The Transcriptomes of Two Heritable Cell Types Illuminate the Circuit Governing Their Differentiation

The differentiation of cells into distinct cell types, each of which is heritable for many generations, underlies many biological phenomena. White and opaque cells of the fungal pathogen Candida albicans are two such heritable cell types, each thought to be adapted to unique niches within their human host. To systematically investigate their differences, we performed strand-specific, massively-parallel sequencing of RNA from C. albicans white and opaque cells. With these data we first annotated the C. albicans transcriptome, finding hundreds of novel differentially-expressed transcripts. Using the new annotation, we compared differences in transcript abundance between the two cell types with the genomic regions bound by a master regulator of the white-opaque switch (Wor1). We found that the revised transcriptional landscape considerably alters our understanding of the circuit governing differentiation. In particular, we can now resolve the poor concordance between binding of a master regulator and the differential expression of adjacent genes, a discrepancy observed in several other studies of cell differentiation. More than one third of the Wor1-bound differentially-expressed transcripts were previously unannotated, which explains the formerly puzzling presence of Wor1 at these positions along the genome. Many of these newly identified Wor1-regulated genes are non-coding and transcribed antisense to coding transcripts. We also find that 5′ and 3′ UTRs of mRNAs in the circuit are unusually long and that 5′ UTRs often differ in length between cell-types, suggesting UTRs encode important regulatory information and that use of alternative promoters is widespread. Further analysis revealed that the revised Wor1 circuit bears several striking similarities to the Oct4 circuit that specifies the pluripotency of mammalian embryonic stem cells. Additional characteristics shared with the Oct4 circuit suggest a set of general hallmarks characteristic of heritable differentiation states in eukaryotes

Diverse Forms of RPS9 Splicing Are Part of an Evolving Autoregulatory Circuit

Ribosomal proteins are essential to life. While the functions of ribosomal protein-encoding genes (RPGs) are highly conserved, the evolution of their regulatory mechanisms is remarkably dynamic. In Saccharomyces cerevisiae, RPGs are unusual in that they are commonly present as two highly similar gene copies and in that they are over-represented among intron-containing genes. To investigate the role of introns in the regulation of RPG expression, we constructed 16 S. cerevisiae strains with precise deletions of RPG introns. We found that several yeast introns function to repress rather than to increase steady-state mRNA levels. Among these, the RPS9A and RPS9B introns were required for cross-regulation of the two paralogous gene copies, which is consistent with the duplication of an autoregulatory circuit. To test for similar intron function in animals, we performed an experimental test and comparative analyses for autoregulation among distantly related animal RPS9 orthologs. Overexpression of an exogenous RpS9 copy in Drosophila melanogaster S2 cells induced alternative splicing and degradation of the endogenous copy by nonsense-mediated decay (NMD). Also, analysis of expressed sequence tag data from distantly related animals, including Homo sapiens and Ciona intestinalis, revealed diverse alternatively-spliced RPS9 isoforms predicted to elicit NMD. We propose that multiple forms of splicing regulation among RPS9 orthologs from various eukaryotes operate analogously to translational repression of the alpha operon by S4, the distant prokaryotic ortholog. Thus, RPS9 orthologs appear to have independently evolved variations on a fundamental autoregulatory circuit

Using RNA-seq to determine the transcriptional landscape and the hypoxic response of the pathogenic yeast Candida parapsilosis

<p>Abstract</p> <p>Background</p> <p><it>Candida parapsilosis </it>is one of the most common causes of <it>Candida </it>infection worldwide. However, the genome sequence annotation was made without experimental validation and little is known about the transcriptional landscape. The transcriptional response of <it>C. parapsilosis </it>to hypoxic (low oxygen) conditions, such as those encountered in the host, is also relatively unexplored.</p> <p>Results</p> <p>We used next generation sequencing (RNA-seq) to determine the transcriptional profile of <it>C. parapsilosis </it>growing in several conditions including different media, temperatures and oxygen concentrations. We identified 395 novel protein-coding sequences that had not previously been annotated. We removed > 300 unsupported gene models, and corrected approximately 900. We mapped the 5' and 3' UTR for thousands of genes. We also identified 422 introns, including two introns in the 3' UTR of one gene. This is the first report of 3' UTR introns in the Saccharomycotina. Comparing the introns in coding sequences with other species shows that small numbers have been gained and lost throughout evolution. Our analysis also identified a number of novel transcriptional active regions (nTARs). We used both RNA-seq and microarray analysis to determine the transcriptional profile of cells grown in normoxic and hypoxic conditions in rich media, and we showed that there was a high correlation between the approaches. We also generated a knockout of the <it>UPC2 </it>transcriptional regulator, and we found that similar to <it>C. albicans</it>, Upc2 is required for conferring resistance to azole drugs, and for regulation of expression of the ergosterol pathway in hypoxia.</p> <p>Conclusion</p> <p>We provide the first detailed annotation of the <it>C. parapsilosis </it>genome, based on gene predictions and transcriptional analysis. We identified a number of novel ORFs and other transcribed regions, and detected transcripts from approximately 90% of the annotated protein coding genes. We found that the transcription factor Upc2 role has a conserved role as a major regulator of the hypoxic response in <it>C. parapsilosis </it>and <it>C. albicans</it>.</p

Large-Scale Evidence for Conservation of NMD Candidature Across Mammals

BACKGROUND: Alternatively-spliced (AS) forms can vary protein function, intracellular localization and post-translational modifications. AS coupled with mRNA nonsense-mediated decay (NMD) can also control the transcript abundance. Here, we have investigated the genome-scale conservation of alternatively-spliced NMD candidates (AS-NMD candidates), in mammals. METHODOLOGY/PRINCIPAL FINDINGS: We mapped>12 million cDNA/EST library transcripts, comprising pooled data from both older and next-generation sequencing techniques, against genomic sequences to annotate AS-NMD candidates generated by in-frame premature termination codons (PTCs), in the human, mouse, rat and cow genomes. In these genomes, we found populations of genes that harbour AS-NMD candidates, varying in number from approximately 149 to 2,051 genes. We discovered that a highly-significant proportion (27%-35%) of AS-NMD candidate genes in mouse, rat and cow, also have human orthologs targeted for NMD. Intron retention was the most abundant type of AS-NMD, ranging from 43% to 67% of genes harbouring an AS-NMD candidate. Groupings of AS-NMD candidate genes either with or without intron retentions also have highly significant AS-NMD conservation, indicating that the trend is not due primarily to conservation of intron retentions. As a subset, the AS-NMD intron retentions are distinguished from non-retained introns by higher GC content, and codon usage similar to the usage in protein-coding sequences. This indicates that most of these alternatively spliced sequences have coded for proteins in the recent evolutionary past. In general, the AS-NMD candidate genes showed a similar pattern of Gene Ontology functional category enrichments in all four species. Genes linked to nucleic-acid interaction and apoptosis, and involved in pathways linked with cancer, were the most common. Finally, we mapped the AS-NMD candidates to mass spectrometry-derived proteomics data, and gathered evidence of truncated polypeptides for at least 10% of all human AS-NMD candidate transcripts. CONCLUSIONS/SIGNIFICANCE: In summary, our analysis provides strong statistical evidence for conservation of functional AS-NMD candidature across Mammalia for a large subset of genes. However, because codon usage of AS-NMD intron retentions is similar to the usage in exons, it is difficult to de-couple conservation of AS-NMD-based regulation from conservation for protein-coding ability, for intron retentions

Heritability in the Efficiency of Nonsense-Mediated mRNA Decay in Humans

BACKGROUND: In eukaryotes mRNA transcripts of protein-coding genes in which an intron has been retained in the coding region normally result in premature stop codons and are therefore degraded through the nonsense-mediated mRNA decay (NMD) pathway. There is evidence in the form of selective pressure for in-frame stop codons in introns and a depletion of length three introns that this is an important and conserved quality-control mechanism. Yet recent reports have revealed that the efficiency of NMD varies across tissues and between individuals, with important clinical consequences. PRINCIPAL FINDINGS: Using previously published Affymetrix exon microarray data from cell lines genotyped as part of the International HapMap project, we investigated whether there are heritable, inter-individual differences in the abundance of intron-containing transcripts, potentially reflecting differences in the efficiency of NMD. We identified intronic probesets using EST data and report evidence of heritability in the extent of intron expression in 56 HapMap trios. We also used a genome-wide association approach to identify genetic markers associated with intron expression. Among the top candidates was a SNP in the DCP1A gene, which forms part of the decapping complex, involved in NMD. CONCLUSIONS: While we caution that some of the apparent inter-individual difference in intron expression may be attributable to different handling or treatments of cell lines, we hypothesize that there is significant polymorphism in the process of NMD, resulting in heritable differences in the abundance of intronic mRNA. Part of this phenotype is likely to be due to a polymorphism in a decapping enzyme on human chromosome 3

Nonsense-mediated mRNA decay controls the changes in yeast ribosomal protein pre-mRNAs levels upon osmotic stress

The expression of ribosomal protein (RP) genes requires a substantial part of cellular transcription, processing and translation resources. Thus, the RP expression must be tightly regulated in response to conditions that compromise cell survival. In Saccharomyces cerevisiae cells, regulation of the RP gene expression at the transcriptional, mature mRNA stability and translational levels during the response to osmotic stress has been reported. Reprogramming global protein synthesis upon osmotic shock includes the movement of ribosomes from RP transcripts to stress-induced mRNAs. Using tiling arrays, we show that osmotic stress yields a drop in the levels of RP pre-mRNAs in S. cerevisiae cells. An analysis of the tiling array data, together with transcription rates data, shows a poor correlation, indicating that the drop in the RP pre-mRNA levels is not merely a result of the lowered RP transcription rates. A kinetic study using quantitative RT-PCR confirmed the decrease in the levels of several RP-unspliced transcripts during the first 15 minutes of osmotic stress, which seems independent of MAP kinase Hog1. Moreover, we found that the mutations in the components of the nonsense-mediated mRNA decay (NMD), Upf1, Upf2, Upf3 or in exonuclease Xrn1, eliminate the osmotic stress-induced drop in RP pre-mRNAs. Altogether, our results indicate that the degradation of yeast RP unspliced transcripts by NMD increases during osmotic stress, and suggest that this might be another mechanism to control RP synthesis during the stress response

An RNA Transport System in Candida albicans Regulates Hyphal Morphology and Invasive Growth

Localization of specific mRNAs is an important mechanism through which cells achieve polarity and direct asymmetric growth. Based on a framework established in Saccharomyces cerevisiae, we describe a She3-dependent RNA transport system in Candida albicans, a fungal pathogen of humans that grows as both budding (yeast) and filamentous (hyphal and pseudohyphal) forms. We identify a set of 40 mRNAs that are selectively transported to the buds of yeast-form cells and to the tips of hyphae, and we show that many of the genes encoded by these mRNAs contribute to hyphal development, as does the transport system itself. Although the basic system of mRNA transport is conserved between S. cerevisiae and C. albicans, we find that the cargo mRNAs have diverged considerably, implying that specific mRNAs can easily move in and out of transport control over evolutionary timescales. The differences in mRNA cargos likely reflect the distinct selective pressures acting on the two species

Cryptic Transcription Mediates Repression of Subtelomeric Metal Homeostasis Genes

Nonsense-mediated mRNA decay (NMD) prevents the accumulation of transcripts bearing premature termination codons. Here we show that Saccharomyces cerevisiae NMD mutants accumulate 5′–extended RNAs (CD-CUTs) of many subtelomeric genes. Using the subtelomeric ZRT1 and FIT3 genes activated in response to zinc and iron deficiency, respectively, we show that transcription of these CD-CUTs mediates repression at the bona fide promoters, by preventing premature binding of RNA polymerase II in conditions of metal repletion. Expression of the main ZRT1 CD-CUT is controlled by the histone deacetylase Rpd3p, showing that histone deacetylases can regulate expression of genes through modulation of the level of CD-CUTs. Analysis of binding of the transcriptional activator Zap1p and insertion of transcriptional terminators upstream from the Zap1p binding sites show that CD-CUT transcription or accumulation also interferes with binding of the transcriptional activator Zap1p. Consistent with this model, overexpressing Zap1p or using a constitutively active version of the Aft1p transcriptional activator rescues the induction defect of ZRT1 and FIT3 in NMD mutants. These results show that cryptic upstream sense transcription resulting in unstable transcripts degraded by NMD controls repression of a large number of genes located in subtelomeric regions, and in particular of many metal homeostasis genes

Transcriptome Analysis of the Model Protozoan, Tetrahymena thermophila, Using Deep RNA Sequencing

Background: The ciliated protozoan Tetrahymena thermophila is a well-studied single-celled eukaryote model organism for cellular and molecular biology. However, the lack of extensive T. thermophila cDNA libraries or a large expressed sequence tag (EST) database limited the quality of the original genome annotation. Methodology/Principal Findings: This RNA-seq study describes the first deep sequencing analysis of the T. thermophila transcriptome during the three major stages of the life cycle: growth, starvation and conjugation. Uniquely mapped reads covered more than 96 % of the 24,725 predicted gene models in the somatic genome. More than 1,000 new transcribed regions were identified. The great dynamic range of RNA-seq allowed detection of a nearly six order-of-magnitude range of measurable gene expression orchestrated by this cell. RNA-seq also allowed the first prediction of transcript untranslated regions (UTRs) and an updated (larger) size estimate of the T. thermophila transcriptome: 57 Mb, or about 55 % of the somatic genome. Our study identified nearly 1,500 alternative splicing (AS) events distributed over 5.2 % of T. thermophila genes. This percentage represents a two order-of-magnitude increase over previous EST-based estimates in Tetrahymena. Evidence of stage-specific regulation of alternative splicing was also obtained. Finally, our study allowed us to completely confirm about 26.8 % of the genes originally predicted by the gene finder, to correct coding sequence boundaries an

A Critical Analysis of Atoh7 (Math5) mRNA Splicing in the Developing Mouse Retina

The Math5 (Atoh7) gene is transiently expressed during retinogenesis by progenitors exiting mitosis, and is essential for ganglion cell (RGC) development. Math5 contains a single exon, and its 1.7 kb mRNA encodes a 149-aa polypeptide. Mouse Math5 mutants have essentially no RGCs or optic nerves. Given the importance of this gene in retinal development, we thoroughly investigated the possibility of Math5 mRNA splicing by Northern blot, 3′RACE, RNase protection assays, and RT-PCR, using RNAs extracted from embryonic eyes and adult cerebellum, or transcribed in vitro from cDNA clones. Because Math5 mRNA contains an elevated G+C content, we used graded concentrations of betaine, an isostabilizing agent that disrupts secondary structure. Although ∼10% of cerebellar Math5 RNAs are spliced, truncating the polypeptide, our results show few, if any, spliced Math5 transcripts exist in the developing retina (<1%). Rare deleted cDNAs do arise via RT-mediated RNA template switching in vitro, and are selectively amplified during PCR. These data differ starkly from a recent study (Kanadia and Cepko 2010), which concluded that the vast majority of Math5 and other bHLH transcripts are spliced to generate noncoding RNAs. Our findings clarify the architecture of the Math5 gene and its mechanism of action. These results have implications for all members of the bHLH gene family, for any gene that is alternatively spliced, and for the interpretation of all RT-PCR experiments