Gene expression in IFN-gamma-activated murine macrophages

Macrophages are critical for natural immunity and play a central role in specific acquired immunity. The IFN-gamma activation of macrophages derived from A/J or BALB/c mice yielded two different patterns of antiviral state in murine hepatitis virus 3 infection, which were related to a down-regulation of the main virus receptor. Using cDNA hybridization to evaluate mRNA accumulation in the cells, we were able to identify several genes that are differently up- or down-regulated by IFN-gamma in A/J (267 and 266 genes, respectively, up- and down-regulated) or BALB/c (297 and 58 genes, respectively, up- and down-regulated) mouse macrophages. Macrophages from mice with different genetic backgrounds behave differently at the molecular level and comparison of the patterns of non-activated and IFN-gamma-activated A/J or BALB/c mouse macrophages revealed, for instance, an up-regulation and a down-regulation of genes coding for biological functions such as enzymatic reactions, nucleic acid synthesis and transport, protein synthesis, transport and metabolism, cytoskeleton arrangement and extracellular matrix, phagocytosis, resistance and susceptibility to infection and tumors, inflammation, and cell differentiation or activation. The present data are reported in order to facilitate future correlation of proteomic/transcriptomic findings as well as of results obtained from a classical approach for the understanding of biological phenomena. The possible implication of the role of some of the gene products relevant to macrophage biology can now be further scrutinized. In this respect, a down-regulation of the main murine hepatitis virus 3 receptor gene was detected only in IFN-gamma-activated macrophages of resistant mice.

A conceptual and practical overview of cDNA microarray technology: implications for basic and clinical sciences

cDNA microarray is an innovative technology that facilitates the analysis of the expression of thousands of genes simultaneously. The utilization of this methodology, which is rapidly evolving, requires a combination of expertise from the biological, mathematical and statistical sciences. In this review, we attempt to provide an overview of the principles of cDNA microarray technology, the practical concerns of the analytical processing of the data obtained, the correlation of this methodology with other data analysis methods such as immunohistochemistry in tissue microarrays, and the cDNA microarray application in distinct areas of the basic and clinical sciences.

Pim-1 kinase inhibits the activation of reporter gene expression in Elk-1 and c-Fos reporting systems but not the endogenous gene expression: an artifact of the reporter gene assay by transient co-transfection

We have studied the molecular mechanism and signal transduction of pim-1, an oncogene encoding a serine-threonine kinase. This is a true oncogene which prolongs survival and inhibits apoptosis of hematopoietic cells. In order to determine whether the effects of Pim-1 occur by regulation of the mitogen-activated protein kinase pathway, we used a transcriptional reporter assay by transient co-transfection as a screening method. In this study, we found that Pim-1 inhibited the Elk-1 and NFkappaB transcriptional activities induced by activation of the mitogen-activated protein kinase cascade in reporter gene assays. However, Western blots showed that the induction of Elk-1-regulated expression of endogenous c-Fos was not affected by Pim-1. The phosphorylation and activation of neither Erk1/2 nor Elk-1 was influenced by Pim-1. Also, in the gel shift assay, the pattern of endogenous NFkappaB binding to its probe was not changed in any manner by Pim-1. These data indicate that Pim-1 does not regulate the activation of Erk1/2, Elk-1 or NFkappaB. These contrasting results suggest a pitfall of the transient co-transfection reporter assay in analyzing the regulation of transcription factors outside of the chromosome context. It ensures that results from reporter gene expression assay should be verified by study of endogenous gene expression.

Gene expression profiling of clinical stages II and III breast cancer

Clinical stage (CS) is an established indicator of breast cancer outcome. In the present study, a cDNA microarray platform containing 692 genes was used to identify molecular differences between CSII and CSIII disease. Tumor samples were collected from patients with CSII or CSIII breast cancer, and normal breast tissue was collected from women without invasive cancer. Seventy-eight genes were deregulated in CSIII tumors and 22 in CSII tumors when compared to normal tissue, and 20 of them were differentially expressed in both CSII and CSIII tumors. In addition, 58 genes were specifically altered in CSIII and expression of 6 of them was tested by real time RT-PCR in another cohort of patients with CSII or CSIII breast cancer and in women without cancer. Among these genes, MAX, KRT15 and S100A14, but not APOBEC3G or KRT19, were differentially expressed on both CSIII and CSII tumors as compared to normal tissue. Increased HMOX1 levels were detected only in CSIII tumors and may represent a molecular marker of this stage. A clear difference in gene expression pattern occurs at the normal-to-cancer transition; however, most of the differentially expressed genes are deregulated in tumors of both CS (II and III) compared to normal breast tissue.

Hypothyroidism decreases proinsulin gene expression and the attachment of its mRNA and eEF1A protein to the actin cytoskeleton of INS-1E cells

The actions of thyroid hormone (TH) on pancreatic beta cells have not been thoroughly explored, with current knowledge being limited to the modulation of insulin secretion in response to glucose, and beta cell viability by regulation of pro-mitotic and pro-apoptotic factors. Therefore, the effects of TH on proinsulin gene expression are not known. This led us to measure: a) proinsulin mRNA expression, b) proinsulin transcripts and eEF1A protein binding to the actin cytoskeleton, c) actin cytoskeleton arrangement, and d) proinsulin mRNA poly(A) tail length modulation in INS-1E cells cultured in different media containing: i) normal fetal bovine serum - FBS (control); ii) normal FBS plus 1 µM or 10 nM T3, for 12 h, and iii) FBS depleted of TH for 24 h (Tx). A decrease in proinsulin mRNA content and attachment to the cytoskeleton were observed in hypothyroid (Tx) beta cells. The amount of eEF1A protein anchored to the cytoskeleton was also reduced in hypothyroidism, and it is worth mentioning that eEF1A is essential to attach transcripts to the cytoskeleton, which might modulate their stability and rate of translation. Proinsulin poly(A) tail length and cytoskeleton arrangement remained unchanged in hypothyroidism. T3 treatment of control cells for 12 h did not induce any changes in the parameters studied. The data indicate that TH is important for proinsulin mRNA expression and translation, since its total amount and attachment to the cytoskeleton are decreased in hypothyroid beta cells, providing evidence that effects of TH on carbohydrate metabolism also include the control of proinsulin gene expression.

Arginine induces GH gene expression by activating NOS/NO signaling in rat isolated hemi-pituitaries

The amino acid arginine (Arg) is a recognized secretagogue of growth hormone (GH), and has been shown to induce GH gene expression. Arg is the natural precursor of nitric oxide (NO), which is known to mediate many of the effects of Arg, such as GH secretion. Arg was also shown to increase calcium influx in pituitary cells, which might contribute to its effects on GH secretion. Although the mechanisms involved in the effects of Arg on GH secretion are well established, little is known about them regarding the control of GH gene expression. We investigated whether the NO pathway and/or calcium are involved in the effects of Arg on GH gene expression in rat isolated pituitaries. To this end, pituitaries from approximately 170 male Wistar rats (~250 g) were removed, divided into two halves, pooled (three hemi-pituitaries) and incubated or not with Arg, as well as with different pharmacological agents. Arg (71 mM), the NO donor sodium nitroprusside (SNP, 1 and 0.1 mM) and a cyclic guanosine monophosphate (cGMP) analogue (8-Br-cGMP, 1 mM) increased GH mRNA expression 60 min later. The NO acceptor hemoglobin (0.3 µM) blunted the effect of SNP, and the combined treatment with Arg and L-NAME (a NO synthase (NOS) inhibitor, 55 mM) abolished the stimulatory effect of Arg on GH gene expression. The calcium channel inhibitor nifedipine (3 µM) also abolished Arg-induced GH gene expression. The present study shows that Arg directly induces GH gene expression in hemi-pituitaries isolated from rats, excluding interference from somatostatinergic neurons, which are supposed to be inhibited by Arg. Moreover, the data demonstrate that the NOS/NO signaling pathway and calcium mediate the Arg effects on GH gene expression.

Probe-level linear model fitting and mixture modeling results in high accuracy detection of differential gene expression

Affiliation: Institut de recherche en immunologie et en cancérologie, Université de Montréal

Regulation of gene expression in the dinoflagellate Lingulodinium polyedrum

Les dinoflagellés sont des eucaryotes unicellulaires que l’on retrouve autant en eau douce qu’en milieu marin. Ils sont particulièrement connus pour causer des fleurs d’algues toxiques nommées ‘marée-rouge’, ainsi que pour leur symbiose avec les coraux et pour leur importante contribution à la fixation du carbone dans les océans. Au point de vue moléculaire, ils sont aussi connus pour leur caractéristiques nucléaires uniques, car on retrouve généralement une quantité immense d’ADN dans leurs chromosomes et ceux-ci sont empaquetés et condensés sous une forme cristalline liquide au lieu de nucléosomes. Les gènes encodés par le noyau sont souvent présents en multiples copies et arrangés en tandem et aucun élément de régulation transcriptionnelle, y compris la boite TATA, n’a encore été observé. L’organisation unique de la chromatine des dinoflagellés suggère que différentes stratégies sont nécessaires pour contrôler l’expression des gènes de ces organismes. Dans cette étude, j’ai abordé ce problème en utilisant le dinoflagellé photosynthétique Lingulodinium polyedrum comme modèle. L. polyedrum est d’un intérêt particulier, car il a plusieurs rythmes circadiens (journalier). À ce jour, toutes les études sur l’expression des gènes lors des changements circadiens ont démontrées une régulation à un niveau traductionnel. Pour mes recherches, j’ai utilisé les approches transcriptomique, protéomique et phosphoprotéomique ainsi que des études biochimiques pour donner un aperçu de la mécanique de la régulation des gènes des dinoflagellés, ceci en mettant l’accent sur l’importance de la phosphorylation du système circadien de L. polyedrum. L’absence des protéines histones et des nucléosomes est une particularité des dinoflagellés. En utilisant la technologie RNA-Seq, j’ai trouvé des séquences complètes encodant des histones et des enzymes modifiant les histones. L polyedrum exprime donc des séquences conservées codantes pour les histones, mais le niveau d’expression protéique est plus faible que les limites de détection par immunodétection de type Western. Les données de séquençage RNA-Seq ont également été utilisées pour générer un transcriptome, qui est une liste des gènes exprimés par L. polyedrum. Une recherche par homologie de séquences a d’abord été effectuée pour classifier les transcrits en diverses catégories (Gene Ontology; GO). Cette analyse a révélé une faible abondance des facteurs de transcription et une surprenante prédominance, parmi ceux-ci, des séquences à domaine Cold Shock. Chez L. polyedrum, plusieurs gènes sont répétés en tandem. Un alignement des séquences obtenues par RNA-Seq avec les copies génomiques de gènes organisés en tandem a été réalisé pour examiner la présence de transcrits polycistroniques, une hypothèse formulée pour expliquer le manque d’élément promoteur dans la région intergénique de la séquence de ces gènes. Cette analyse a également démontré une très haute conservation des séquences codantes des gènes organisés en tandem. Le transcriptome a également été utilisé pour aider à l’identification de protéines après leur séquençage par spectrométrie de masse, et une fraction enrichie en phosphoprotéines a été déterminée comme particulièrement bien adapté aux approches d’analyse à haut débit. La comparaison des phosphoprotéomes provenant de deux périodes différentes de la journée a révélée qu’une grande partie des protéines pour lesquelles l’état de phosphorylation varie avec le temps est reliées aux catégories de liaison à l’ARN et de la traduction. Le transcriptome a aussi été utilisé pour définir le spectre des kinases présentes chez L. polyedrum, qui a ensuite été utilisé pour classifier les différents peptides phosphorylés qui sont potentiellement les cibles de ces kinases. Plusieurs peptides identifiés comme étant phosphorylés par la Casein Kinase 2 (CK2), une kinase connue pour être impliquée dans l’horloge circadienne des eucaryotes, proviennent de diverses protéines de liaison à l’ARN. Pour évaluer la possibilité que quelques-unes des multiples protéines à domaine Cold Shock identifiées dans le transcriptome puissent moduler l’expression des gènes de L. polyedrum, tel qu’observé chez plusieurs autres systèmes procaryotiques et eucaryotiques, la réponse des cellules à des températures froides a été examinée. Les températures froides ont permis d’induire rapidement un enkystement, condition dans laquelle ces cellules deviennent métaboliquement inactives afin de résister aux conditions environnementales défavorables. Les changements dans le profil des phosphoprotéines seraient le facteur majeur causant la formation de kystes. Les phosphosites prédits pour être phosphorylés par la CK2 sont la classe la plus fortement réduite dans les kystes, une découverte intéressante, car le rythme de la bioluminescence confirme que l’horloge a été arrêtée dans le kyste.

Role of histone methylation in the regulation of COX-2, iNOS, and mPGES-1 gene expression in human chondrocytes: Implication for Osteoarthritis

L'arthrose (OA) est une maladie articulaire dégénérative, classée comme la forme la plus fréquente au monde. Elle est caractérisée par la dégénérescence du cartilage articulaire, l’inflammation de la membrane synoviale, et le remodelage de l’os sous-chondral. Ces changements structurels et fonctionnels sont dues à de nombreux facteurs. Les cytokines, les prostaglandines (PG), et les espèces réactives de l'oxygène sont les principaux médiateurs impliqués dans la pathophysiologie de l'OA. L'interleukine-1β (IL-1β) est une cytokine pro-inflammatoire majeure qui joue un rôle crucial dans l'OA. L'IL-1β induit l'expression de la cyclooxygénase-2 (COX-2), la microsomale prostaglandine E synthase-1 (mPGES-1), la synthase inductible de l'oxyde nitrique (iNOS), ainsi que leurs produits la prostaglandine E2 (PGE2) et l'oxyde nitrique (NO). Ce sont des médiateurs essentiels de la réponse inflammatoire au cours de l'OA qui contribuent aux mécanismes des douleurs, de gonflement, et de destruction des tissus articulaires. Les modifications épigénétiques jouent un rôle très important dans la régulation de l’expression de ces gènes pro-inflammatoires. Parmi ces modifications, la méthylation/ déméthylation des histones joue un rôle critique dans la régulation des gènes. La méthylation/ déméthylation des histones est médiée par deux types d'enzymes: les histones méthyltransférases (HMT) et les histones déméthylases (HDM) qui favorisent l’activation et/ou la répression de la transcription. Il est donc nécessaire de comprendre les mécanismes moléculaires qui contrôlent l’expression des gènes de la COX-2, la mPGES-1, et l’iNOS. L'objectif de cette étude est de déterminer si la méthylation/déméthylation des histones contribute à la régulation de l’expression des gènes COX-2, mPGES-1, et iNOS dans des chondrocytes OA humains induits par l'IL-1β. Nous avons montré que la méthylation de la lysine K4 de l'histone H3 (H3K4) par SET-1A contribue à l’activation des gènes COX-2 et iNOS dans les chondrocytes humains OA induite par l'IL-1β. Nous avons également montré que la lysine K9 de l’histone H3 (H3K9) est déméthylée par LSD1, et que cette déméthylation contribue à l’expression de la mPGES-1 induite par IL-1β dans les chondrocytes humains OA. Nous avons aussi trouvé que les niveaux d'expression des enzymes SET-1A et LSD1 sont élevés au niveau du cartilage OA. Nos résultats montrent, pour la première fois, l'implication de la méthylation/ déméthylation des histones dans la régulation de l’expression des gènes COX-2, mPGES-1, et iNOS. Ces données suggèrent que ces mécanismes pourraient être une cible potentielle pour une intervention pharmacologique dans le traitement de la physiopathologie de l'OA.

GABAA and GABAB Receptor Gene Expression and Functional Regulation During Pancreatic Regeneration and Insulin Secretion in Rats

The present study demonstrate the functional alterations of the GABAA and GABAB receptors and the gene expression during the regeneration of pancreas following partial pancreatectomy. The role of these receptors in insulin secretion and pancreatic DNA synthesis using the specific agonists and antagonists also are studied in vitro. The alterations of GABAA and GABAR receptor function and gene expression in the brain stem, crebellum and hypothalamus play an important role in the sympathetic regulation of insulin secretion during pancreatic regeneration. Previous studies have given much information linking functional interaction between GABA and the peripheral nervous system. The involvement of specific receptor subtypes functional regulation during pancreatic regeneration has not given emphasis and research in this area seems to be scarce. We have observed a decreased GABA content, down regulation of GABAA receptors and an up regulation of GABAB receptors in the cerebral cortex, brain stem and hypothalamus. Real Time-PCR analysis confirmed the receptor data in the brain regions. These alterations in the GABAA and GABAB receptors of the brain are suggested to govern the regenerative response and growth regulation of the pancreas through sympathetic innervation. In addition, receptor binding studies and Real Time-PCR analysis revealed that during pancreatic regeneration GABAA receptors were down regulated and GABAB receptors were up regulated in pancreatic islets. This suggests an inhibitory role for GABAA receptors in islet cell proliferation i.e., the down regulation of this receptor facilitates proliferation. Insulin secretion study during 1 hour showed GABA has inhibited the insulin secretion in a dose dependent manner in normal and hyperglycaemic conditions. Bicuculline did not antagonize this effect. GABAA agonist, muscimol inhibited glucose stimulated insulin secretion from pancreatic islets except in the lowest concentration of 1O-9M in presence of 4mM glucose.Musclmol enhanced insulin secretion at 10-7 and 10-4M muscimol in presence of 20mM glucose- 4mM glucose represents normal and 20mM represent hyperglycaemic conditions. GABAB agonist, baclofen also inhibited glucose induced insulin secretion and enhanced at the concentration of 1O-5M at 4mM glucose and at 10-9M baclofen in presence of 20mM glucose. This shows a differential control of the GABAA and GABAB receptors over insulin release from the pancreatic islets. During 24 hours in vitro insulin secretion study it showed that low concentration of GABA has inhibited glucose stimulated insulin secretion from pancreatic islets. Muscimol, the GABAA agonist, inhibited the insulin secretion but, gave an enhanced secretion of insulin in presence of 4mM glucose at 10-7 , 10-5 and 1O-4M muscimol. But in presence of 20mM glucose muscimol significantly inhibited the insulin secretion. GABAB agonist, baclofen also inhibited glucose induced insulin secretion in presence of both 4mM and 20mM glucose. This shows the inhibitory role of GABA and its specific receptor subtypes over insulin synthesis from pancreatic bete-islets. In vitro DNA synthesis studies showed that activation of GABAA receptor by adding muscimol, a specific agonist, inhibited islet DNA synthesis. Also, the addition of baclofen, a specific agonist of GABAB receptor resulted in the stimulation of DNA synthesis.Thus the brain and pancreatic GABAA and GABAB receptor gene expression differentially regulates pancreatic insulin secretion and islet cell proliferation during pancreatic regeneration. This will have immense clinical significance in therapeutic applications in the management of Diabetes mellitus.

“Unmixing” Tissue Gene Expression Signatures from Tumor Biopsies

Emergent molecular measurement methods, such as DNA microarray, qRTPCR, and many others, offer tremendous promise for the personalized treatment of cancer. These technologies measure the amount of specific proteins, RNA, DNA or other molecular targets from tumor specimens with the goal of “fingerprinting” individual cancers. Tumor specimens are heterogeneous; an individual specimen typically contains unknown amounts of multiple tissues types. Thus, the measured molecular concentrations result from an unknown mixture of tissue types, and must be normalized to account for the composition of the mixture. For example, a breast tumor biopsy may contain normal, dysplastic and cancerous epithelial cells, as well as stromal components (fatty and connective tissue) and blood and lymphatic vessels. Our diagnostic interest focuses solely on the dysplastic and cancerous epithelial cells. The remaining tissue components serve to “contaminate” the signal of interest. The proportion of each of the tissue components changes as a function of patient characteristics (e.g., age), and varies spatially across the tumor region. Because each of the tissue components produces a different molecular signature, and the amount of each tissue type is specimen dependent, we must estimate the tissue composition of the specimen, and adjust the molecular signal for this composition. Using the idea of a chemical mass balance, we consider the total measured concentrations to be a weighted sum of the individual tissue signatures, where weights are determined by the relative amounts of the different tissue types. We develop a compositional source apportionment model to estimate the relative amounts of tissue components in a tumor specimen. We then use these estimates to infer the tissuespecific concentrations of key molecular targets for sub-typing individual tumors. We anticipate these specific measurements will greatly improve our ability to discriminate between different classes of tumors, and allow more precise matching of each patient to the appropriate treatment

Juvenile-mature wood transition in pine: correlation between wood properties and candidate gene expression profiles

There is a strong desire to exploit transcriptomics data from model species for the genetic improvement of non-model crops. Here, we use gene expression profiles from the commercial model Pinus taeda to identify candidate genes implicated in juvenile-mature wood transition in the non-model relative, P. sylvestris. Re-analysis of 'public domain' SAGE data from xylem tissues of P. taeda revealed 283 mature-abundant and 396 juvenile-abundant tags (P < 0.01), of which 70 and 137, respectively matched to genes with known function. Based on sequence similarity, we then isolated 16 putative homologues of genes that in P. taeda exhibited widest divergence in expression between juvenile and mature samples. Candidate expression levels in P. sylvestris were almost invariably differential between juvenile and mature woody tissue samples among two cohorts of five trees collected from the same seed source and selected for genetic uniformity by genetic distance analysis. However, the direction of differential expression was not always consistent with that described in the original P. taeda SAGE data. Correlation was observed between gene expression and juvenile-mature wood anatomical characteristics by OPLS analysis. Four candidates (alpha-tubulin, porin MIP1, lipid transfer protein and aquaporin like protein) apparently had greatest influence on the wood traits measured. Speculative function of these genes in relation to juvenile-mature wood transition is briefly explored. Thus, we demonstrate the feasibility of exploiting SAGE data from a model species to identify consistently differentially expressed candidates in a related non-model species.

Comparison of the global gene expression profiles produced by methylparaben, n-butylparaben and 17 beta-oestradiol in MCF7 human breast cancer cells

Since the alkyl esters of p-hydroxybenzoic acid (parabens) can be measured intact in the human breast and possess oestrogenic properties, it has been suggested that they could contribute to an aberrant burden of oestrogen signalling in the human breast and so play a role in the rising incidence of breast cancer. However, although parabens have been shown to regulate a few single genes (reporter genes, pS2, progesterone receptor) in a manner similar to that of 17 beta-oestradiol, the question remains as to the full extent of the similarity in the overall gene profile induced in response to parabens compared with 17 beta-oestradiol. The GE-Amersham CodeLink 20 K human expression microarray system was used to profile the expression of 19881 genes in MCF7 human breast cancer cells following a 7-day exposure to 5 x 10(-4) m methylparaben, 10(-5) m n-butylparaben and 10(-8) m 17 beta-oestradiol. At these concentrations, the parabens gave growth responses in MCF7 cells of similar magnitude to 17 beta-oestradiol. The study identified genes which are upregulated or downregulated to a similar extent by methylparaben, n-butylparaben and 17 beta-oestradiol. However, the majority of genes were not regulated in the same way by all three treatments. Some genes responded differently to parabens from 17 beta-oestradiol, and furthermore, differences in expression of some genes could be detected even between the two individual parabens. Therefore, although parabens possess oestrogenic properties, their mimicry in terms of global gene expression patterns is not perfect and differences in gene expression profiles could result in consequences to the cells that are not identical to those following exposure to 17 beta-oestradiol. Copyright (c) 2006 John Wiley & Sons, Ltd.

Development and application of in vivo expression technology (IVET) for analysing microbial gene expression in complex environments

Establishing the mechanisms by which microbes interact with their environment, including eukaryotic hosts, is a major challenge that is essential for the economic utilisation of microbes and their products. Techniques for determining global gene expression profiles of microbes, such as microarray analyses, are often hampered by methodological restraints, particularly the recovery of bacterial transcripts (RNA) from complex mixtures and rapid degradation of RNA. A pioneering technology that avoids this problem is In Vivo Expression Technology (IVET). IVET is a 'promoter-trapping' methodology that can be used to capture nearly all bacterial promoters (genes) upregulated during a microbe-environment interaction. IVET is especially useful because there is virtually no limit to the type of environment used (examples to date include soil, oomycete, a host plant or animal) to select for active microbial promoters. Furthermore, IVET provides a powerful method to identify genes that are often overlooked during genomic annotation, and has proven to be a flexible technology that can provide even more information than identification of gene expression profiles. A derivative of IVET, termed resolvase-IVET (RIVET), can be used to provide spatio-temporal information about environment-specific gene expression. More recently, niche-specific genes captured during an IVET screen have been exploited to identify the regulatory mechanisms controlling their expression. Overall, IVET and its various spin-offs have proven to be a valuable and robust set of tools for analysing microbial gene expression in complex environments and providing new targets for biotechnological development.

Dietary alpha-tocopherol affects differential gene expression in rat testes

Gene-chip technology was employed to study the effect of dietary vitamin E (VE) on gene expression in rat testes. Male albino rats were fed with either a diet deficient in VE or a standard diet containing VE. Differential gene expression was monitored at five individual time-points over a period of 14 months with all animals individually pro. led. Low VE intake resulted in the consistent upregulation of 7-dehydrocholesterol reductase and GATA binding protein 4, both involved in testosterone synthesis. Cyclin D3, important in cell cycle progression and Wilms tumor 1, related to cancer development, were also up-regulated in the vitamin E deficient animals. This study demonstrates that low dietary VE intake has long-term effects on gene expression in the testes. Our data provides insights into the possible molecular mechanisms underlying the beneficial effects of vitamin E on the male reproductive organ.