Concordance among digital gene expression, microarrays, and qPCR when measuring differential expression of microRNAs.

Profiling microRNA (miRNA) expression is of widespread interest given the critical role of miRNAs in many cellular functions. Profiling can be achieved via hybridization-based (microarrays), sequencing-based, or amplification-based (quantitative reverse transcription-PCR, qPCR) technologies. Among these, microarrays face the significant challenge of accurately distinguishing between mature and immature miRNA forms, and different vendors have developed different methods to meet this challenge. Here we measure differential miRNA expression using the Affymetrix, Agilent, and Illumina microarray platforms, as well as qPCR (Applied Biosystems) and ultra high-throughput sequencing (Illumina). We show that the differential expression measurements are more divergent when the three types of microarrays are compared than when the Agilent microarray, qPCR, and sequencing technology measurements are compared, which exhibit a good overall concordance.

Expression of the peroxisome proliferator-activated receptor alpha gene is stimulated by stress and follows a diurnal rhythm.

Peroxisome proliferator-activated receptors (PPARs) are nuclear hormone receptors that can be activated by fatty acids and peroxisome proliferators. The PPAR alpha subtype mediates the pleiotropic effects of these activators in liver and regulates several target genes involved in fatty acid catabolism. In primary hepatocytes cultured in vitro, the PPAR alpha gene is regulated at the transcriptional level by glucocorticoids. We investigated if this hormonal regulation also occurs in the whole animal in physiological situations leading to increased plasma corticosterone levels in rats. We show here that an immobilization stress is a potent and rapid stimulator of PPAR alpha expression in liver but not in hippocampus. The injection of the synthetic glucocorticoid dexamethasone into adult rats produces a similar increase in PPAR alpha expression in liver, whereas the administration of the antiglucocorticoid RU 486 inhibits the stress-dependent stimulation. We conclude that glucocorticoids are major mediators of the stress response. Consistent with this hormonal regulation, hepatic PPAR alpha mRNA and protein levels follow a diurnal rhythm, which parallels that of circulating corticosterone. To test the effects of variations in PPAR alpha expression on PPAR alpha target gene activity, high glucocorticoid-dependent PPAR alpha expression was mimicked in cultured primary hepatocytes. Under these conditions, hormonal stimulation of receptor expression synergizes with receptor activation by WY-14,643 to induce the expression of the PPAR alpha target gene acyl-CoA oxidase. Together, these results show that regulation of the PPAR alpha expression levels efficiently modulates PPAR activator signaling and thus may affect downstream metabolic pathways involved in lipid homeostasis.

Polarity and specific sequence requirements of peroxisome proliferator-activated receptor (PPAR)/retinoid X receptor heterodimer binding to DNA. A functional analysis of the malic enzyme gene PPAR response element.

The malic enzyme (ME) gene is a target for both thyroid hormone receptors and peroxisome proliferator-activated receptors (PPAR). Within the ME promoter, two direct repeat (DR)-1-like elements, MEp and MEd, have been identified as putative PPAR response elements (PPRE). We demonstrate that only MEp and not MEd is able to bind PPAR/retinoid X receptor (RXR) heterodimers and mediate peroxisome proliferator signaling. Taking advantage of the close sequence resemblance of MEp and MEd, we have identified crucial determinants of a PPRE. Using reciprocal mutation analyses of these two elements, we show the preference for adenine as the spacing nucleotide between the two half-sites of the PPRE and demonstrate the importance of the two first bases flanking the core DR1 in 5'. This latter feature of the PPRE lead us to consider the polarity of the PPAR/RXR heterodimer bound to its cognate element. We demonstrate that, in contrast to the polarity of RXR/TR and RXR/RAR bound to DR4 and DR5 elements respectively, PPAR binds to the 5' extended half-site of the response element, while RXR occupies the 3' half-site. Consistent with this polarity is our finding that formation and binding of the PPAR/RXR heterodimer requires an intact hinge T region in RXR while its integrity is not required for binding of the RXR/TR heterodimer to a DR4.

Understanding mutational and selective processes that govern gene duplication in mammals

Abstract : Gene duplication is an essential source of material for the origin of genetic novelties. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with at least ~3600 detectable retrocopies. We find that ~30% of these retrocopies are transcribed, generally in testes. Their transcription often relies on preexisting regulatory elements (or open chromatin) close to their insertion site, which is illustrated by mRNA molecules containing retrocopies fused to their neighboring genes. Retrocopies appear to have been profoundly shaped by selection. Consistently, human retrocopies with an intact open reading (ORF) are more often transcribed than retropseudogenes, which leads to a minimal estimate of 120 functional retrogenes present in our genome. We also performed an analysis of Ka/Ks for human retrocopies. This analysis demonstrates that several intact retrocopies evolved under purifying selection and yields an estimated formation rate of ~1 retrogene per million year in the primate lineage. Using DNA sequencing and evolutionary simulations, we have identified 7 such primate-specific retrogenes that emerged on the lineage leading to humans In therian genomes, we found an excess of retrogenes with X-linked parents. Expression analyses support the idea that this "out of X" movement was driven by natural selection to produce autosomal functional counterparts for X-linked genes, which are silenced during male meiosis. Phylogenetic dating of this "out of X" movement suggests that our sex chromosomes arose about 180 MYA ago and are thus much younger than previously thought. Finally, we have also analyzed young gene duplications (and deletions) that arose by non allelic-homologous recombination and are not fixed in species. Using wild-caught and laboratory animals, we detected thousands of DNA segments that are polymorphic in copy number in mice. These copy number variants were found to profoundly alter the transcriptome of several mouse tissues. Strikingly, their influence on gene expression is not limited to the gene they contain but seems to extend to genes located up to 1.5 million bases away.

Immuno-electron microscopic identification of human estrogen receptor-DNA complexes at the estrogen-responsive element and in the first intron of a Xenopus vitellogenin gene.

Using an extract of nuclei from the estrogen-responsive human breast cancer cell line MCF-7, protein-DNA complexes were assembled in vitro at the 5' end of the Xenopus laevis vitellogenin gene B2 that is normally expressed in liver after estrogen induction. The complexes formed were analyzed by electron microscopy after labeling by the indirect colloidal gold immunological method using a monoclonal antibody specific for the human estrogen receptor. As identified by its interaction with protein A-gold, the antibody was found linked to two protein-DNA complexes, the first localized at the estrogen responsive element of the gene and the second in intron I, thus proving a direct participation of the receptor in these two complexes. The procedure used allows the visualization and rapid localization of specific transcription factors bound in vitro to a promoter or any other gene region.

Estrogen receptor level determines sex-specific in vitro transcription from the Xenopus vitellogenin promoter.

Female-specific expression of the Xenopus laevis vitellogenin gene was reconstituted in vitro by addition of recombinant vaccinia-virus-produced estrogen receptor to nuclear extracts from male livers, in which this gene is silent. Transcription enhancement was at least 30 times and was selectively restricted to vitellogenin templates containing the estrogen-responsive unit. Thus, in male hepatocytes, estrogen receptor is the limiting regulatory factor that in the female liver controls efficient and accurate sex-specific expression of the vitellogenin gene. Furthermore, the Xenopus liver factor B, which is essential in addition to the estrogen receptor for the activation of this gene, was successfully replaced in the Xenopus extract by purified human nuclear factor I, identifying factor B of Xenopus as a functional homolog of this well-characterized human transcription factor.

DNA methylation by CcrM activates the transcription of two genes required for the division of Caulobacter crescentus.

DNA methylation regulates many processes, including gene expression, by superimposing secondary information on DNA sequences. The conserved CcrM enzyme, which methylates adenines in GANTC sequences, is essential to the viability of several Alphaproteobacteria. In this study, we find that Caulobacter crescentus cells lacking the CcrM enzyme accumulate low levels of the two conserved FtsZ and MipZ proteins, leading to a severe defect in cell division. This defect can be compensated by the expression of the ftsZ gene from an inducible promoter or by spontaneous suppressor mutations that promote FtsZ accumulation. We show that CcrM promotes the transcription of the ftsZ and mipZ genes and that the ftsZ and mipZ promoter regions contain a conserved CGACTC motif that is critical to their activities and to their regulation by CcrM. In addition, our results suggest that the ftsZ promoter has the lowest activity when the CGACTC motif is non-methylated, an intermediate activity when it is hemi-methylated and the highest activity when it is fully methylated. The regulation of ftsZ expression by DNA methylation may explain why CcrM is essential in a subset of Alphaproteobacteria.

Members of the PHO1 gene family show limited functional redundancy in phosphate transfer to the shoot, and are regulated by phosphate deficiency via distinct pathways.

The PHO1 family comprises 11 members in Arabidopsis thaliana. In order to decipher the role of these genes in inorganic phosphate (Pi) transport and homeostasis, complementation of the pho1 mutant, deficient in loading Pi to the root xylem, was determined by the expression of the PHO1 homologous genes under the control of the PHO1 promoter. Only PHO1 and the homologue PHO1;H1 could complement pho1. The PHO1;H1 promoter was active in the vascular cylinder of roots and shoots. Expression of PHO1;H1 was very low in Pi-sufficient plants, but was strongly induced under Pi-deficient conditions. T-DNA knock-out mutants of PHO1;H1 neither showed growth defects nor alteration in Pi transport dynamics, or Pi content, compared with wild type. However, the double mutant pho1/pho1;h1 showed a strong reduction in growth and in the capacity to transfer Pi from the root to the shoot compared with pho1. Grafting experiments revealed that phenotypes associated with the pho1 and pho1/pho1;h1 mutants were linked to the lack of gene expression in the root. The increased expression of PHO1;H1 under Pi deficiency was largely controlled by the transcription factor PHR1 and was suppressed by the phosphate analogue phosphite, whereas the increase of PHO1 expression was independent of PHR1 and was not influenced by phosphite. Together, these data reveal that although transfer of Pi to the root xylem vessel is primarily mediated by PHO1, the homologue PHO1;H1 also contributes to Pi loading to the xylem, and that the two corresponding genes are regulated by Pi deficiency by distinct signal transduction pathways.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution.

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Sequence homologies in the region preceding the transcription initiation site of the liver estrogen-responsive vitellogenin and apo-VLDLII genes.

In the liver of oviparous vertebrates vitellogenin gene expression is controlled by estrogen. The nucleotide sequence of the 5' flanking region of the Xenopus laevis vitellogenin genes A1, A2, B1 and B2 has been determined. These sequences have been compared to each other and to the equivalent region of the chicken vitellogenin II and apo-VLDLII genes which are also expressed in the liver in response to estrogen. The homology between the 5' flanking region of the Xenopus genes B1 and B2 is higher than between the corresponding regions of the other closely related genes A1 and A2. Four short blocks of sequence homology which are present at equivalent positions in the vitellogenin genes of both Xenopus laevis and chicken are characterized. A short sequence with two-fold rotational symmetry (GGTCANNNTGACC) was found at similar positions upstream of the five vitellogenin genes and is also present in two copies close to the 5' end of the chicken apo-VLDLII gene. The possible functional significance of this sequence, common to liver estrogen-responsive genes, is discussed.

Differential regulation of Na-K-ATPase isoform gene expression by T3 during rat brain development.

A fetal rat telencephalon organotypic cell culture system was found to reproduce the developmental pattern of Na-K-adenosinetriphosphatase (ATPase) gene expression observed in vivo [Am. J. Physiol. 258 (Cell Physiol. 27): C1062-C1069, 1990]. We have used this culture system to study the effects of triiodothyronine (T3; 0.003-30 nM) on mRNA abundance and basal transcription rates of Na-K-ATPase isoforms. Steady-state mRNA levels were low at culture day 6 (corresponding to the day of birth) but distinct for each isoform alpha 3 much greater than beta 1 = beta 2 greater than alpha 2 greater than alpha 1. At culture day 6, T3 did not modify mRNA abundance of any isoform. At culture day 12 (corresponding to day 7 postnatal), T3 increased the mRNA level of alpha 2 (4- to 7-fold), beta 2 (4- to 5-fold), alpha 1 (3- to 6-fold), and beta 1 (1.5-fold), whereas alpha 3 mRNA levels remained unchanged. Interestingly, the basal transcription rate for each isoform differed strikingly (alpha 2 greater than alpha 1 much greater than beta 1 = beta 2 greater than alpha 3) but remained stable throughout 12 days of culture and was not regulated by T3. Thus we observed an inverse relationship between rate of transcription and rate of mRNA accumulation for each alpha-isoform, suggesting that alpha 1- and alpha 2-mRNA are turning over rapidly whereas alpha 3-mRNA is turning over slowly. Our data indicate that one of the mechanisms by which T3 selectively controls Na-K-ATPase gene expression during brain development in vitro occurs at the posttranscriptional level.

The G0/G1 switch gene 2 is a novel PPAR target gene.

PPARs (peroxisome-proliferator-activated receptors) alpha, beta/delta and gamma are a group of transcription factors that are involved in numerous processes, including lipid metabolism and adipogenesis. By comparing liver mRNAs of wild-type and PPARalpha-null mice using microarrays, a novel putative target gene of PPARalpha, G0S2 (G0/G1 switch gene 2), was identified. Hepatic expression of G0S2 was up-regulated by fasting and by the PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. Surprisingly, the G0S2 mRNA level was highest in brown and white adipose tissue and was greatly up-regulated during mouse 3T3-L1 and human SGBS (Simpson-Golabi-Behmel syndrome) adipogenesis. Transactivation, gel shift and chromatin immunoprecipitation assays indicated that G0S2 is a direct PPARgamma and probable PPARalpha target gene with a functional PPRE (PPAR-responsive element) in its promoter. Up-regulation of G0S2 mRNA seemed to be specific for adipogenesis, and was not observed during osteogenesis or myogenesis. In 3T3-L1 fibroblasts, expression of G0S2 was associated with growth arrest, which is required for 3T3-L1 adipogenesis. Together, these data indicate that G0S2 is a novel target gene of PPARs that may be involved in adipocyte differentiation.

Estrogen-dependent in vitro transcription from the vitellogenin promoter in liver nuclear extracts.

One approach to analyzing the molecular mechanisms of gene expression in vivo is to reconstitute these events in cell-free systems in vitro. Although there is some evidence for tissue-specific transcription in vitro, transcriptionally active extracts that mimic a steroid hormone-dependent enhancement of transcription have not been described. In the study reported here, nuclear extracts of liver from the frog Xenopus laevis were capable of estrogen-dependent induction of a homologous vitellogenin promoter that contained the estrogen-responsive element.

SEPT4, a p53 target gene product, implicates p 53 in the regulation of exocytosis

Résumé : Le cancer, qui est responsable d'un quart des décès en Suisse, exhibe un état cellulaire désordonné, qui lui-même, est la conséquence d'un dérèglement des gènes. Le gène le plus fréquemment altéré, dans les cas de cancers humains, est p53. Ce gène encode un facteur de transcription, impliqué dans la régulation de nombreux gènes impliqués dans le cycle cellulaire, l'apoptose ou la différenciation. Notre laboratoire a récemment identifié seize nouveaux gènes, dont l'expression est régulée par p53, parmi lesquels sept4, su jet de cette thèse. La protéine 5EPT4 appartient à la famille des septines, qui est impliquée dans la cytokinèse. Dans ce travail, nous avons confirmé la régulation de l'expression de sept4 par p53 dans des tissus de souris, et étonnamment, seul un des deux promoteurs du gène sept4 est contrôlé par p53. En outre, l'approche immunohistologique nous a permis de supposer une implication de la protéine SEPT4 dans le mécanisme de l'exocytose. Cette hypothèse a été confirmée par l'interaction de SEPT4 avec la protéine syntaxine 1A, et par son activité inhibitrice sur la sécrétion stimulée. En élargissant l'étude de la protéine SEPT4, nous avons découvert que celle-ci avait comme partenaire fonctionnel, la protéine Pinl, une enzyme qui catalyse l'isomérisation cis-trans du lien peptidique précédant une proline. bans ce contexte, nous avons démontré que l'interaction entre ces deux protéines reposait sur le domaine WW de Pinl, un type de domaine reconnaissant les motifs phosphoséryl-prolyl et phosphothréonyl-prolyl. Ce dernier résultat nous a conduit à examiner la phosphorylation de 5EPT4. Nous avons démontré que la partie N-terminale de SEPT4 était phosphorylée par la kinase Cdk5. La co¬expression de Cdk5 et de SEPT4 stimule la dégradation de SEPT4, indépendamment de la voie du protéasome. Ainsi, l'ensemble de nos observations fournissent l'évidence de l'engagement de la protéine SEPT4 dans la régulation de l'exocytose, et soutiennent le rôle de p53 dans le contrôle de l'exocytose, via SEPT4, ce qui constituerait un nouveau rôle fonctionnel pour ce gardien du génome. Summary: Cancer, which is responsible for a quarter of the deaths in Switzerland, exhibits a disordered cellular state, which itself, is the consequence of an altered state of genes. The most frequently altered gene in human cancer is p53. This gene encodes a transcription factor, implicated in the regulation of numerous genes involved in cell cycle, apoptosis or differentiation. Our laboratory has recently identified sixteen new genes whose expression is regulated by p53, amongst them septin 4, which is the subject of this thesis. The SEPT4 protein belongs to the septin family which is implicated in cytokinesis. In the present work, we have confirmed the regulation of sept4 expression by p53 in mouse tissues, and surprisingly, only one of the two sept4 promoters is regulated by p53. In addition, the immunohistologic approach enabled us to suppose a role of SEPT4 in exocytosis. This assumption was confirmed by the interaction of SEPT4 with syntaxin 1A, and by its inhibiting activity on stimulated secretion. By widening the analysis of SEPT4, we identified Pin1 as an interacting protein. Pin1 is an enzyme which catalyzes the cis-trans isomerization of the peptide bond preceding a proline residue. In this context, we showed that the interaction between these two proteins depends on the WW domain of Pin 1. This domain has been shown to function as a phosphoserine- or phosphothreonine¬binding module. This last result prompted us to examine phosphorylation of SEPT4. We demonstrated that the N-terminal portion of SEPT4 was phosphorylated by the Cdk5 kinase. The co-expression of Cdk5 with 5EPT4 stimulates SEPT4 degradation, independently of the proteasome pathway. Thus, these observations provide evidence for the engagement of SEPT4 in the regulation of exocytosis, and supports the role of p53 in the control of exocytosis, via SEPT4, which constitutes a new functional role for this guardian of the genome.

Membrane mu poly(A) signal and 3' flanking sequences function as a transcription terminator for immunoglobulin-encoding genes.

Developmentally regulated mechanisms involving alternative RNA splicing and/or polyadenylation, as well as transcription termination, are implicated in controlling the levels of secreted mu (mu s), membrane mu (mu m) and delta immunoglobulin (Ig) heavy chain mRNAs during B cell differentiation (mu gene encodes the mu heavy chain). Using expression vectors constructed with genomic DNA segments composed of the mu m polyadenylation signal region, we analyzed poly(A) site utilization and termination of transcription in stably transfected myeloma cells and in murine fibroblast L cells. We found that the gene segment containing the mu m poly(A) signals, along with 536 bp of downstream flanking sequence, acted as a transcription terminator in both myeloma cells and L cell fibroblasts. Neither a 141-bp DNA fragment (which directed efficient polyadenylation at the mu m site), nor the 536-bp flanking nucleotide sequence alone, were sufficient to obtain a similar regulation. This shows that the mu m poly(A) region plays a central role in controlling developmentally regulated transcription termination by blocking downstream delta gene expression. Because this gene segment exhibited the same RNA processing and termination activities in fibroblasts, it appears that these processes are not tissue-specific.