Global coordination of transcriptional control and mRNA decay during cellular differentiation

Membrane and organelle rearrangement during ascospore formation in budding yeast

SUMMARYIn ascomycete fungi, sexual spores, termed ascospores, are formed after meiosis. Ascospore formation is an unusual cell division in which daughter cells are created within the cytoplasm of the mother cell by de novo generation of membranes that encapsulate each of the haploid chromosome sets created by meiosis. This review describes the molecular events underlying the creation, expansion, and closure of these membranes in the budding yeast, Saccharomyces cerevisiae. Recent advances in our understanding of the regulation of gene expression and the dynamic behavior of different membrane-bound organelles during this process are detailed. While less is known about ascospore formation in other systems, comparison to the distantly related fission yeast suggests that the molecular events will be broadly similar throughout the ascomycetes.

Exploring the role of transcriptional and post-transcriptional processes in mRNA co-expression

Co-expression of two or more genes at the single-cell level is usually associated with functional co-regulation. While mRNA co-expression-measured as the correlation in mRNA levels-can be influenced by both transcriptional and post-transcriptional events, transcriptional regulation is typically considered dominant. We review and connect the literature describing transcriptional and post-transcriptional regulation of co-expression. To enhance our understanding, we integrate four datasets spanning single-cell gene expression data, single-cell promoter activity data and individual transcript half-lives. Confirming expectations, we find that positive co-expression necessitates promoter coordination and similar mRNA half-lives. Surprisingly, negative co-expression is favored by differences in mRNA half-lives, contrary to initial predictions from stochastic simulations. Notably, this association manifests specifically within clusters of genes. We further observe a striking compensation between promoter coordination and mRNA half-lives, which additional stochastic simulations suggest might give rise to the observed co-expression patterns. These findings raise intriguing questions about the functional advantages conferred by this compensation between distal kinetic steps.

A Trans-Omics Comparison Reveals Common Gene Expression Strategies in Four Model Organisms and Exposes Similarities and Differences between Them

The ultimate goal of gene expression regulation is on the protein level. However, because the amounts of mRNAs and proteins are controlled by their synthesis and degradation rates, the cellular amount of a given protein can be attained by following different strategies. By studying omics data for six expression variables (mRNA and protein amounts, plus their synthesis and decay rates), we previously demonstrated the existence of common expression strategies (CESs) for functionally related genes in the yeast Saccharomyces cerevisiae. Here we extend that study to two other eukaryotes: the yeast Schizosaccharomyces pombe and cultured human HeLa cells. We also use genomic data from the model prokaryote Escherichia coli as an external reference. We show that six-variable profiles (6VPs) can be constructed for every gene and that these 6VPs are similar for genes with similar functions in all the studied organisms. The differences in 6VPs between organisms can be used to establish their phylogenetic relationships. The analysis of the correlations among the six variables supports the hypothesis that most gene expression control occurs in actively growing organisms at the transcription rate level, and that translation plays a minor role. We propose that living organisms use CESs for the genes acting on the same physiological pathways, especially for those belonging to stable macromolecular complexes, but CESs have been modeled by evolution to adapt to the specific life circumstances of each organism.

The Fission Yeast RNA-Binding Protein Meu5 Is Involved in Outer Forespore Membrane Breakdown during Spore Formation

In Schizosaccharomyces pombe, the spore wall confers strong resistance against external stress. During meiosis II, the double-layered intracellular forespore membrane (FSM) forms de novo and encapsulates the nucleus. Eventually, the inner FSM layer becomes the plasma membrane of the spore, while the outer layer breaks down. However, the molecular mechanism and biological significance of this membrane breakdown remain unknown. Here, by genetic investigation of an S. pombe mutant (E22) with normal prespore formation but abnormal spores, we showed that Meu5, an RNA-binding protein known to bind to and stabilize more than 80 transcripts, is involved in this process. We confirmed that the E22 mutant does not produce Meu5 protein, while overexpression of meu5⁺ in E22 restores the sporulation defect. Furthermore, electron microscopy revealed that the outer membrane of the FSM persisted in meu5∆ spores. Investigation of the target genes of meu5⁺ showed that a mutant of cyc1⁺ encoding cytochrome c also showed a severe defect in outer FSM breakdown. Lastly, we determined that outer FSM breakdown occurs coincident with or after formation of the outermost Isp3 layer of the spore wall. Collectively, our data provide novel insights into the molecular mechanism of spore formation.

mmi1 and rep2 mRNAs are novel RNA targets of the Mei2 RNA-binding protein during early meiosis in Schizosaccharomyces pombe

The RNA-binding protein Mei2 is crucial for meiosis in Schizosaccharomyces pombe. In mei2 mutants, pre-meiotic S-phase is blocked, along with meiosis. Mei2 binds a long non-coding RNA (lncRNA) called meiRNA, which is a 'sponge RNA' for the meiotic inhibitor protein Mmi1. The interaction between Mei2, meiRNA and Mmi1 protein is essential for meiosis. But mei2 mutants have stronger and different phenotypes than meiRNA mutants, since mei2Δ arrests before pre-meiotic S, while the meiRNA mutant arrests after pre-meiotic S but before meiosis. This suggests Mei2 may bind additional RNAs. To identify novel RNA targets of Mei2, which might explain how Mei2 regulates pre-meiotic S, we used RNA immunoprecipitation and cross-linking immunoprecipitation. In addition to meiRNA, we found the mRNAs for mmi1 (which encodes Mmi1) and for the S-phase transcription factor rep2 There were also three other RNAs of uncertain relevance. We suggest that at meiotic initiation, Mei2 may sequester rep2 mRNA to help allow pre-meiotic S, and then may bind both meiRNA and mmi1 mRNA to inactivate Mmi1 at two levels, the protein level (as previously known), and also the mRNA level, allowing meiosis. We call Mei2-meiRNA a 'double sponge' (i.e. binding both an mRNA and its encoded protein).

Linear Regression Links Transcriptomic Data and Cellular Raman Spectra

Raman microscopy is an imaging technique that has been applied to assess molecular compositions of living cells to characterize cell types and states. However, owing to the diverse molecular species in cells and challenges of assigning peaks to specific molecules, it has not been clear how to interpret cellular Raman spectra. Here, we provide firm evidence that cellular Raman spectra and transcriptomic profiles of Schizosaccharomyces pombe and Escherichia coli can be computationally connected and thus interpreted. We find that the dimensions of high-dimensional Raman spectra and transcriptomes measured by RNA sequencing can be reduced and connected linearly through a shared low-dimensional subspace. Accordingly, we were able to predict global gene expression profiles by applying the calculated transformation matrix to Raman spectra, and vice versa. Highly expressed non-coding RNAs contributed to the Raman-transcriptome linear correspondence more significantly than mRNAs in S. pombe. This demonstration of correspondence between cellular Raman spectra and transcriptomes is a promising step toward establishing spectroscopic live-cell omics studies.

Gene Expression of Pneumocystis murina after Treatment with Anidulafungin Results in Strong Signals for Sexual Reproduction, Cell Wall Integrity, and Cell Cycle Arrest, Indicating a Requirement for Ascus Formation for Proliferation

The echinocandins are a class of antifungal agents that target β-1,3-d-glucan (BG) biosynthesis. In the ascigerous Pneumocystis species, treatment with these drugs depletes the ascus life cycle stage, which contains BG, but large numbers of forms which do not express BG remain in the infected lungs. In the present study, the gene expression profiles of Pneumocystis murina were compared between infected, untreated mice and mice treated with anidulafungin for 2 weeks to understand the metabolism of the persisting forms. Almost 80 genes were significantly up- or downregulated. Like other fungi exposed to echinocandins, genes associated with sexual replication, cell wall integrity, cell cycle arrest, and stress comprised the strongest upregulated signals in P. murina from the treated mice. The upregulation of the P. murina β-1,3-d-glucan endohydrolase and endo-1,3-glucanase was notable and may explain the disappearance of the existing asci in the lungs of treated mice since both enzymes can degrade BG. The biochemical measurement of BG in the lungs of treated mice and fluorescence microscopy with an anti-BG antibody supported the loss of BG. Downregulated signals included genes involved in cell replication, genome stability, and ribosomal biogenesis and function and the Pneumocystis-specific genes encoding the major surface glycoproteins (Msg). These studies suggest that P. murina attempted to undergo sexual replication in response to a stressed environment and was halted in any type of proliferative cycle, likely due to a lack of BG. Asci appear to be a required part of the life cycle stage of Pneumocystis, and BG may be needed to facilitate progression through the life cycle via sexual replication.

Silencing of the lncRNA Zeb2-NAT facilitates reprogramming of aged fibroblasts and safeguards stem cell pluripotency

Aging imposes a barrier to somatic cell reprogramming through poorly understood mechanisms. Here, we report that fibroblasts from old mice express higher levels of Zeb2, a transcription factor that activates epithelial-to-mesenchymal transition. Synthesis of Zeb2 protein is controlled by a natural antisense transcript named Zeb2-NAT. We show that transfection of adult fibroblasts with specific LNA Gapmers induces a robust downregulation of Zeb2-NAT transcripts and Zeb2 protein and enhances the reprogramming of old fibroblasts into pluripotent cells. We further demonstrate that Zeb2-NAT expression is precociously activated by differentiation stimuli in embryonic stem (ES) cells. By knocking down Zeb2-NAT, we were able to maintain ES cells challenged with commitment signals in the ground state of pluripotency. In conclusion, our study identifies a long noncoding RNA that is overlapping and antisense to the Zeb2 locus as a target for rejuvenation strategies.

The efficiency of somatic cell reprogramming is lowered by ageing. Here the authors show that the transcription factor Zeb2 and its long non-coding RNA Zeb2-NAT are expressed at high levels in older fibroblasts and their inhibition increases reprogramming efficiency.

TORC1 and TORC2 converge to regulate the SAGA co-activator in response to nutrient availability

Gene expression regulation is essential for cells to adapt to changes in their environment. Co-activator complexes have well-established roles in transcriptional regulation, but less is known about how they sense and respond to signaling cues. We have previously shown that, in fission yeast, one such co-activator, the SAGA complex, controls gene expression and the switch from proliferation to differentiation in response to nutrient availability. Here, using a combination of genetic, biochemical, and proteomic approaches, we show that SAGA responds to nutrients through the differential phosphorylation of its Taf12 component, downstream of both the TORC1 and TORC2 pathways. Taf12 phosphorylation increases early upon starvation and is controlled by the opposing activities of the PP2A phosphatase, which is activated by TORC1, and the TORC2-activated Gad8^AKT kinase. Mutational analyses suggest that Taf12 phosphorylation prevents cells from committing to differentiation until starvation reaches a critical level. Overall, our work reveals that SAGA is a direct target of nutrient-sensing pathways and has uncovered a mechanism by which TORC1 and TORC2 converge to control gene expression and cell fate decisions.

Rapid Recovery Gene Downregulation during Excess-Light Stress and Recovery in Arabidopsis

Stress recovery may prove to be a promising approach to increase plant performance and, theoretically, mRNA instability may facilitate faster recovery. Transcriptome (RNA-seq, qPCR, sRNA-seq, and PARE) and methylome profiling during repeated excess-light stress and recovery was performed at intervals as short as 3 min. We demonstrate that 87% of the stress-upregulated mRNAs analyzed exhibit very rapid recovery. For instance, HSP101 abundance declined 2-fold every 5.1 min. We term this phenomenon rapid recovery gene downregulation (RRGD), whereby mRNA abundance rapidly decreases promoting transcriptome resetting. Decay constants (k) were modeled using two strategies, linear and nonlinear least squares regressions, with the latter accounting for both transcription and degradation. This revealed extremely short half-lives ranging from 2.7 to 60.0 min for 222 genes. Ribosome footprinting using degradome data demonstrated RRGD loci undergo cotranslational decay and identified changes in the ribosome stalling index during stress and recovery. However, small RNAs and 5'-3' RNA decay were not essential for recovery of the transcripts examined, nor were any of the six excess light-associated methylome changes. We observed recovery-specific gene expression networks upon return to favorable conditions and six transcriptional memory types. In summary, rapid transcriptome resetting is reported in the context of active recovery and cellular memory.

4-Thiouridine Labeling to Analyze mRNA Turnover in Schizosaccharomyces pombe

Traditionally, the half-lives of mRNAs were measured after inhibition of transcription to allow decay of the preexisting population. The protocol presented here is a more recently developed strategy in which mRNA turnover is analyzed by measuring the decline in levels of newly synthesized RNA labeled with 4-thiouridine (4sU) during a brief pulse. After RNA extraction, the 4sU is biotinylated and the labeled species are purified using streptavidin beads. DNA microarrays can then be used to compare this population with total RNA, allowing half-lives to be calculated.

Analysis of RNA Metabolism in Fission Yeast

Here we focus on the biogenesis and function of messenger RNA (mRNA) in fission yeast cells. Following a general introduction that also briefly touches on other classes of RNA, we provide an overview of methods used to analyze mRNAs throughout their life cycles.

Predicted RNA Binding Proteins Pes4 and Mip6 Regulate mRNA Levels, Translation, and Localization during Sporulation in Budding Yeast

In response to starvation, diploid cells of Saccharomyces cerevisiae undergo meiosis and form haploid spores, a process collectively referred to as sporulation. The differentiation into spores requires extensive changes in gene expression. The transcriptional activator Ndt80 is a central regulator of this process, which controls many genes essential for sporulation. Ndt80 induces ∼300 genes coordinately during meiotic prophase, but different mRNAs within the NDT80 regulon are translated at different times during sporulation. The protein kinase Ime2 and RNA binding protein Rim4 are general regulators of meiotic translational delay, but how differential timing of individual transcripts is achieved was not known. This report describes the characterization of two related NDT80-induced genes, PES4 and MIP6, encoding predicted RNA binding proteins. These genes are necessary to regulate the steady-state expression, translational timing, and localization of a set of mRNAs that are transcribed by NDT80 but not translated until the end of meiosis II. Mutations in the predicted RNA binding domains within PES4 alter the stability of target mRNAs. PES4 and MIP6 affect only a small portion of the NDT80 regulon, indicating that they act as modulators of the general Ime2/Rim4 pathway for specific transcripts.

Ribosome Profiling for the Analysis of Translation During Yeast Meiosis

Ribosome profiling provides a genome-wide view of translation with unprecedented resolution. Application of this approach to fission and budding yeast revealed widespread regulation of translational efficiency, translation of short open reading frames on unannotated transcripts, and frequent translation of open reading frames in 5' leader sequences. We present here a detailed protocol for the application of ribosome profiling to meiotic fission yeast cells, although the approach should be easily adapted to budding yeast.

Analysis of the association between codon optimality and mRNA stability in Schizosaccharomyces pombe

Background

Recent experiments have shown that codon optimality is a major determinant of mRNA stability in Saccharomyces cerevisiae and that this phenomenon may be conserved in Escherichia coli and some metazoans, although work in Neurospora crassa is not consistent with this model.

Results

We examined the association between codon optimality and mRNA stability in the fission yeast Schizosaccharomyces pombe. Our analysis revealed the following points. First, we observe a genome-wide association between codon optimality and mRNA stability also in S. pombe, suggesting evolutionary conservation of the phenomenon. Second, in both S. pombe and S. cerevisiae, mRNA synthesis rates are also correlated at the genome-wide analysis with codon optimality, suggesting that the long-appreciated association between codon optimality and mRNA abundance is due to regulation of both mRNA synthesis and degradation. However, when we examined correlation of codon optimality and either mRNA half-lives or synthesis rates controlling for mRNA abundance, codon optimality was still positively correlated with mRNA half-lives in S. cerevisiae, but the association was no longer significant for mRNA half-lives in S. pombe or for synthesis rates in either organism. This illustrates how only the pairwise analysis of multiple correlating variables may limit these types of analyses. Finally, in S. pombe, codon optimality is associated with known DNA/RNA sequence motifs that are associated with mRNA production/stability, suggesting these two features have been under similar selective pressures for optimal gene expression.

Conclusions

Consistent with the emerging body of studies, this study suggests that the association between codon optimality and mRNA stability may be a broadly conserved phenomenon. It also suggests that the association can be explained at least in part by independent adaptations of codon optimality and other transcript features for elevated expression during evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3237-6) contains supplementary material, which is available to authorized users.

Cell types differ in global coordination of splicing and proportion of highly expressed genes

Balance in the transcriptome is regulated by coordinated synthesis and degradation of RNA molecules. Here we investigated whether mammalian cell types intrinsically differ in global coordination of gene splicing and expression levels. We analyzed RNA-seq transcriptome profiles of 8 different purified mouse cell types. We found that different cell types vary in proportion of highly expressed genes and the number of alternatively spliced transcripts expressed per gene, and that the cell types that express more variants of alternatively spliced transcripts per gene are those that have higher proportion of highly expressed genes. Cell types segregated into two clusters based on high or low proportion of highly expressed genes. Biological functions involved in negative regulation of gene expression were enriched in the group of cell types with low proportion of highly expressed genes, and biological functions involved in regulation of transcription and RNA splicing were enriched in the group of cell types with high proportion of highly expressed genes. Our findings show that cell types differ in proportion of highly expressed genes and the number of alternatively spliced transcripts expressed per gene, which represent distinct properties of the transcriptome and may reflect intrinsic differences in global coordination of synthesis, splicing, and degradation of RNA molecules.

Cuf2 Is a Transcriptional Co-Regulator that Interacts with Mei4 for Timely Expression of Middle-Phase Meiotic Genes

The Schizosaccharomyces pombe cuf2+ gene encodes a nuclear regulator that is required for timely activation and repression of several middle-phase genes during meiotic differentiation. In this study, we sought to gain insight into the mechanism by which Cuf2 regulates meiotic gene expression. Using a chromatin immunoprecipitation approach, we demonstrate that Cuf2 is specifically associated with promoters of both activated and repressed target genes, in a time-dependent manner. In case of the fzr1+ gene whose transcription is positively affected by Cuf2, promoter occupancy by Cuf2 results in a concomitant increased association of RNA polymerase II along its coding region. In marked contrast, association of RNA polymerase II with chromatin decreases when Cuf2 negatively regulates target gene expression such as wtf13+. Although Cuf2 operates through a transcriptional mechanism, it is unable to perform its function in the absence of the Mei4 transcription factor, which is a member of the conserved forkhead protein family. Using coimmunoprecipitation experiments, results showed that Cuf2 is a binding partner of Mei4. Bimolecular fluorescence complementation experiments brought further evidence that an association between Cuf2 and Mei4 occurs in the nucleus. Analysis of fzr1+ promoter regions revealed that two FLEX-like elements, which are bound by the transcription factor Mei4, are required for chromatin occupancy by Cuf2. Together, results reported here revealed that Cuf2 and Mei4 co-regulate the timely expression of middle-phase genes during meiosis.

Reconsidering plant memory: Intersections between stress recovery, RNA turnover, and epigenetics

Plants grow in dynamic environments where they can be exposed to a multitude of stressful factors, all of which affect their development, yield, and, ultimately, reproductive success. Plants are adept at rapidly acclimating to stressful conditions and are able to further fortify their defenses by retaining memories of stress to enable stronger or more rapid responses should an environmental perturbation recur. Indeed, one mechanism that is often evoked regarding environmental memories is epigenetics. Yet, there are relatively few examples of such memories; neither is there a clear understanding of their duration, considering the plethora of stresses in nature. We propose that this field would benefit from investigations into the processes and mechanisms enabling recovery from stress. An understanding of stress recovery could provide fresh insights into when, how, and why environmental memories are created and regulated. Stress memories may be maladaptive, hindering recovery and affecting development and potential yield. In some circumstances, it may be advantageous for plants to learn to forget. Accordingly, the recovery process entails a balancing act between resetting and memory formation. During recovery, RNA metabolism, posttranscriptional gene silencing, and RNA-directed DNA methylation have the potential to play key roles in resetting the epigenome and transcriptome and in altering memory. Exploration of this emerging area of research is becoming ever more tractable with advances in genomics, phenomics, and high-throughput sequencing methodology that will enable unprecedented profiling of high-resolution stress recovery time series experiments and sampling of large natural populations.

AnGeLi: A Tool for the Analysis of Gene Lists from Fission Yeast

Genome-wide assays and screens typically result in large lists of genes or proteins. Enrichments of functional or other biological properties within such lists can provide valuable insights and testable hypotheses. To systematically detect these enrichments can be challenging and time-consuming, because relevant data to compare against query gene lists are spread over many different sources. We have developed AnGeLi (Analysis of Gene Lists), an intuitive, integrated web-tool for comprehensive and customized interrogation of gene lists from the fission yeast, Schizosaccharomyces pombe. AnGeLi searches for significant enrichments among multiple qualitative and quantitative information sources, including gene and phenotype ontologies, genetic and protein interactions, numerous features of genes, transcripts, translation, and proteins such as copy numbers, chromosomal positions, genetic diversity, RNA polymerase II and ribosome occupancy, localization, conservation, half-lives, domains, and molecular weight among others, as well as diverse sets of genes that are co-regulated or lead to the same phenotypes when mutated. AnGeLi uses robust statistics which can be tailored to specific needs. It also provides the option to upload user-defined gene sets to compare against the query list. Through an integrated data submission form, AnGeLi encourages the community to contribute additional curated gene lists to further increase the usefulness of this resource and to get the most from the ever increasing large-scale experiments. AnGeLi offers a rigorous yet flexible statistical analysis platform for rich insights into functional enrichments and biological context for query gene lists, thus providing a powerful exploratory tool through which S. pombe researchers can uncover fresh perspectives and unexpected connections from genomic data. AnGeLi is freely available at: www.bahlerlab.info/AnGeLi.

Role of Ccr4-Not complex in heterochromatin formation at meiotic genes and subtelomeres in fission yeast

BACKGROUND

Heterochromatin is essential for chromosome segregation, gene silencing and genome integrity. The fission yeast Schizosaccharomyces pombe contains heterochromatin at centromeres, subtelomeres, and mating type genes, as well as at small islands of meiotic genes dispersed across the genome. This heterochromatin is generated by partially redundant mechanisms, including the production of small interfering RNAs (siRNAs) that are incorporated into the RITS protein complex (RNAi-Induced Transcriptional Silencing). The assembly of heterochromatin islands requires the function of the RNA-binding protein Mmi1, which recruits RITS to its mRNA targets and to heterochromatin islands. In addition, Mmi1 directs its targets to an exosome-dependent RNA elimination pathway.

RESULTS

Ccr4-Not is a conserved multiprotein complex that regulates gene expression at multiple levels, including RNA degradation and translation. We show here that Ccr4-Not is recruited by Mmi1 to its RNA targets. Surprisingly, Ccr4 and Caf1 (the mRNA deadenylase catalytic subunits of the Ccr4-Not complex) are not necessary for the degradation or translation of Mmi1 RNA targets, but are essential for heterochromatin integrity at Mmi1-dependent islands and, independently of Mmi1, at subtelomeric regions. Both roles require the deadenylase activity of Ccr4 and the Mot2/Not4 protein, a ubiquitin ligase that is also part of the complex. Genetic evidence shows that Ccr4-mediated silencing is essential for normal cell growth, indicating that this novel regulation is physiologically relevant. Moreover, Ccr4 interacts with components of the RITS complex in a Mmi1-independent manner.

CONCLUSIONS

Taken together, our results demonstrate that the Ccr4-Not complex is required for heterochromatin integrity in both Mmi1-dependent and Mmi1-independent pathways.

Kinetics effects and modeling of mRNA turnover

Broader comprehension of gene expression regulatory mechanisms can be gained from a global analysis of how transcription and degradation are coordinated to orchestrate complex cell responses. The role of messenger RNA (mRNA) turnover modulation in gene expression levels has become increasingly recognized. From such perspective, in this review we briefly illustrate how a simple but effective mathematical model of mRNA turnover and some experimental findings, may together shed light on the molecular mechanisms underpinning the major role of mRNA decay rates in shaping the kinetics of gene activation and repression.

Systematic analysis of the role of RNA-binding proteins in the regulation of RNA stability

mRNA half-lives are transcript-specific and vary over a range of more than 100-fold in eukaryotic cells. mRNA stabilities can be regulated by sequence-specific RNA-binding proteins (RBPs), which bind to regulatory sequence elements and modulate the interaction of the mRNA with the cellular RNA degradation machinery. However, it is unclear if this kind of regulation is sufficient to explain the large range of mRNA stabilities. To address this question, we examined the transcriptome of 74 Schizosaccharomyces pombe strains carrying deletions in non-essential genes encoding predicted RBPs (86% of all such genes). We identified 25 strains that displayed changes in the levels of between 4 and 104 mRNAs. The putative targets of these RBPs formed biologically coherent groups, defining regulons involved in cell separation, ribosome biogenesis, meiotic progression, stress responses and mitochondrial function. Moreover, mRNAs in these groups were enriched in specific sequence motifs in their coding sequences and untranslated regions, suggesting that they are coregulated at the posttranscriptional level. We performed genome-wide RNA stability measurements for several RBP mutants, and confirmed that the altered mRNA levels were caused by changes in their stabilities. Although RBPs regulate the decay rates of multiple regulons, only 16% of all S. pombe mRNAs were affected in any of the 74 deletion strains. This suggests that other players or mechanisms are required to generate the observed range of RNA half-lives of a eukaryotic transcriptome.

Coupling mRNA synthesis and decay

What has been will be again, what has been done will be done again; there is nothing new under the sun. -Ecclesiastes 1:9 (New International Version) Posttranscriptional regulation of gene expression has an important role in defining the phenotypic characteristics of an organism. Well-defined steps in mRNA metabolism that occur in the nucleus-capping, splicing, and polyadenylation-are mechanistically linked to the process of transcription. Recent evidence suggests another link between RNA polymerase II (Pol II) and a posttranscriptional process that occurs in the cytoplasm-mRNA decay. This conclusion appears to represent a conundrum. How could mRNA synthesis in the nucleus and mRNA decay in the cytoplasm be mechanistically linked? After a brief overview of mRNA processing, we will review the recent evidence for transcription-coupled mRNA decay and the possible involvement of Snf1, the Saccharomyces cerevisiae ortholog of AMP-activated protein kinase, in this process.

Functional significance of nuclear export and mRNA binding of meiotic regulator Spo5 in fission yeast

Background

Meiotic cells undergo two rounds of nuclear division and generate gametes. Previous studies have indicated that a number of transcription factors modulate the transcriptome in successive waves during meiosis and spore formation in fission yeast. However, the mechanisms underlying the post-transcriptional regulation in meiosis are not fully understood. The fission yeast spo5⁺ gene encodes a meiosis-specific RNA-binding protein, which is required for the progression of meiosis II and spore formation. However, the target RNA molecules of Spo5 are yet to be identified. Characterization of meiosis-specific RNA-binding proteins will provide insight into how post-transcriptional regulation influence gene expression during sexual differentiation.

Results

To assess the functional significance of RNA-recognition motifs (RRMs) of Spo5, we constructed a series of new spo5 truncated mutants and previously reported spo5 missense mutants. In addition, we isolated novel spo5 missense mutants. The phenotypic characteristics of these mutants indicated that the RRMs are essential for both the localization and function of the protein. Interestingly, Spo5 is exported from the nucleus to the cytoplasm via the Rae1-dependent mRNA export pathway, but is unlikely to be involved in global mRNA export. Furthermore, cytoplasmic localization of Spo5 is important for its function, which suggests the involvement of Spo5 in post-transcriptional regulation. We identified pcr1⁺ mRNA as one of the critical targets of Spo5. The pcr1⁺ gene encodes an activating transcription factor/cAMP response element binding (ATF/CREB) transcription factor family. Among the four family members, namely Pcr1, Atf1, Atf21, and Atf31, only the mRNA encoding Pcr1 binds to Spo5.

Conclusions

Spo5 is exported from the nucleus with mRNAs via the Rae1-dependent pathway. RRMs are necessary for this process and also for the function of Spo5 after the nuclear export. Spo5 appears to influence the activity of pcr1⁺ mRNA, and the mechanism of how Spo5 stimulates the mRNA to promote the progression of meiosis II and spore formation remains an intriguing question for future research.

Contributions of transcription and mRNA decay to gene expression dynamics of fission yeast in response to oxidative stress

The cooperation of transcriptional and post-transcriptional levels of control to shape gene regulation is only partially understood. Here we show that a combination of two simple and non-invasive genomic techniques, coupled with kinetic mathematical modeling, afford insight into the intricate dynamics of RNA regulation in response to oxidative stress in the fission yeast Schizosaccharomyces pombe. This study reveals a dominant role of transcriptional regulation in response to stress, but also points to the first minutes after stress induction as a critical time when the coordinated control of mRNA turnover can support the control of transcription for rapid gene regulation. In addition, we uncover specialized gene expression strategies associated with distinct functional gene groups, such as simultaneous transcriptional repression and mRNA destabilization for genes encoding ribosomal proteins, delayed mRNA destabilization with varying contribution of transcription for ribosome biogenesis genes, dominant roles of mRNA stabilization for genes functioning in protein degradation, and adjustment of both transcription and mRNA turnover during the adaptation to stress. We also show that genes regulated independently of the bZIP transcription factor Atf1p are predominantly controlled by mRNA turnover, and identify putative cis-regulatory sequences that are associated with different gene expression strategies during the stress response. This study highlights the intricate and multi-faceted interplay between transcription and RNA turnover during the dynamic regulatory response to stress.

The translational landscape of fission-yeast meiosis and sporulation

Sexual development in the fission yeast Schizosaccharomyces pombe culminates in meiosis and sporulation. We used ribosome profiling to investigate the translational landscape of this process. We show that the translation efficiency of hundreds of genes is regulated in complex patterns, often correlating with changes in RNA levels. Ribosome-protected fragments show a three-nucleotide periodicity that identifies translated sequences and their reading frame. Using this property, we identified 46 novel translated genes and found that 24% of non-coding RNAs are actively translated. We also detected 19 nested antisense genes, in which both DNA strands encode translated mRNAs. Finally, we identified 1,735 translated upstream ORFs in leader sequences. In contrast with Saccharomyces cerevisiae, sexual development in S. pombe is not accompanied by large increases in upstream ORF use, suggesting that this is an organism-specific adaptation and not a general feature of developmental processes.

SnapShot-Seq: a method for extracting genome-wide, in vivo mRNA dynamics from a single total RNA sample

mRNA synthesis, processing, and destruction involve a complex series of molecular steps that are incompletely understood. Because the RNA intermediates in each of these steps have finite lifetimes, extensive mechanistic and dynamical information is encoded in total cellular RNA. Here we report the development of SnapShot-Seq, a set of computational methods that allow the determination of in vivo rates of pre-mRNA synthesis, splicing, intron degradation, and mRNA decay from a single RNA-Seq snapshot of total cellular RNA. SnapShot-Seq can detect in vivo changes in the rates of specific steps of splicing, and it provides genome-wide estimates of pre-mRNA synthesis rates comparable to those obtained via labeling of newly synthesized RNA. We used SnapShot-Seq to investigate the origins of the intrinsic bimodality of metazoan gene expression levels, and our results suggest that this bimodality is partly due to spillover of transcriptional activation from highly expressed genes to their poorly expressed neighbors. SnapShot-Seq dramatically expands the information obtainable from a standard RNA-Seq experiment.

Red5 and three nuclear pore components are essential for efficient suppression of specific mRNAs during vegetative growth of fission yeast

Zinc-finger domains are found in many nucleic acid-binding proteins in both prokaryotes and eukaryotes. Proteins carrying zinc-finger domains have important roles in various nuclear transactions, including transcription, mRNA processing and mRNA export; however, for many individual zinc-finger proteins in eukaryotes, the exact function of the protein is not fully understood. Here, we report that Red5 is involved in efficient suppression of specific mRNAs during vegetative growth of Schizosaccharomyces pombe. Red5, which contains five C3H1-type zinc-finger domains, localizes to the nucleus where it forms discrete dots. A red5 point mutation, red5-2, results in the upregulation of specific meiotic mRNAs in vegetative mutant red5-2 cells; northern blot data indicated that these meiotic mRNAs in red5-2 cells have elongated poly(A) tails. RNA-fluorescence in situ hybridization results demonstrate that poly(A)⁺ RNA species accumulate in the nucleolar regions of red5-deficient cells. Moreover, Red5 genetically interacts with several mRNA export factors. Unexpectedly, three components of the nuclear pore complex also suppress a specific set of meiotic mRNAs. These results indicate that Red5 function is important to meiotic mRNA degradation; they also suggest possible connections among selective mRNA decay, mRNA export and the nuclear pore complex in vegetative fission yeast.

Functional characterization of Upf1 targets in Schizosaccharomyces pombe

Nonsense-mediated mRNA decay (NMD) is a highly conserved mechanism of mRNA degradation. NMD eliminates mRNAs containing premature termination codons (PTCs), preventing the production of truncated proteins with possible deleterious effects. However, there is mounting evidence that NMD factors, like Upf1, Upf2 and Upf3, participate in general regulation of gene expression, affecting the expression of genes lacking PTCs. We have used the fission yeast Schizosaccharomyces pombe to identify mRNAs directly regulated by NMD. Using a combination of genetic and biochemical approaches, we have defined a population of fission yeast mRNAs specifically regulated by Upf1. We show that other components of the Upf complex, Upf2 and Upf3, are required for binding of Upf1 to its RNA targets and for the proper response of fission yeast to oxidative stress. Finally, we investigated the physiological importance of this phenomenon, and demonstrate that the Upf1-dependent downregulation of some of its direct targets is necessary for normal resistance to oxidative stress.

Eukaryotic mRNA decay: methodologies, pathways, and links to other stages of gene expression

mRNA concentration depends on the balance between transcription and degradation rates. On both sides of the equilibrium, synthesis and degradation show, however, interesting differences that have conditioned the evolution of gene regulatory mechanisms. Here, we discuss recent genome-wide methods for determining mRNA half-lives in eukaryotes. We also review pre- and posttranscriptional regulons that coordinate the fate of functionally related mRNAs by using protein- or RNA-based trans factors. Some of these factors can regulate both transcription and decay rates, thereby maintaining proper mRNA homeostasis during eukaryotic cell life.

The fate of the messenger is pre-determined: a new model for regulation of gene expression

Recent years have seen a rise in publications demonstrating coupling between transcription and mRNA decay. This coupling most often accompanies cellular processes that involve transitions in gene expression patterns, for example during mitotic division and cellular differentiation and in response to cellular stress. Transcription can affect the mRNA fate by multiple mechanisms. The most novel finding is the process of co-transcriptional imprinting of mRNAs with proteins, which in turn regulate cytoplasmic mRNA stability. Transcription therefore is not only a catalyst of mRNA synthesis but also provides a platform that enables imprinting, which coordinates between transcription and mRNA decay. Here we present an overview of the literature, which provides the evidence of coupling between transcription and decay, review the mechanisms and regulators by which the two processes are coupled, discuss why such coupling is beneficial and present a new model for regulation of gene expression. This article is part of a Special Issue entitled: RNA Decay mechanisms.

Widespread promoter-mediated coordination of transcription and mRNA degradation

Background

Previous work showed that mRNA degradation is coordinated with transcription in yeast, and in several genes the control of mRNA degradation was linked to promoter elements through two different mechanisms. Here we show at the genomic scale that the coordination of transcription and mRNA degradation is promoter-dependent in yeast and is also observed in humans.

Results

We first demonstrate that swapping upstream cis-regulatory sequences between two yeast species affects both transcription and mRNA degradation and suggest that while some cis-regulatory elements control either transcription or degradation, multiple other elements enhance both processes. Second, we show that adjacent yeast genes that share a promoter (through divergent orientation) have increased similarity in their patterns of mRNA degradation, providing independent evidence for the promoter-mediated coupling of transcription to mRNA degradation. Finally, analysis of the differences in mRNA degradation rates between mammalian cell types or mammalian species suggests a similar coordination between transcription and mRNA degradation in humans.

Conclusions

Our results extend previous studies and suggest a pervasive promoter-mediated coordination between transcription and mRNA degradation in yeast. The diverse genes and regulatory elements associated with this coordination suggest that it is generated by a global mechanism of gene regulation and modulated by gene-specific mechanisms. The observation of a similar coupling in mammals raises the possibility that coupling of transcription and mRNA degradation may reflect an evolutionarily conserved phenomenon in gene regulation.

Dynamic changes in translational efficiency are deduced from codon usage of the transcriptome

Translation of a gene is assumed to be efficient if the supply of the tRNAs that translate it is high. Yet high-abundance tRNAs are often also at high demand since they correspond to preferred codons in genomes. Thus to fully model translational efficiency one must gauge the supply-to-demand ratio of the tRNAs that are required by the transcriptome at a given time. The tRNAs’ supply is often approximated by their gene copy number in the genome. Yet neither the demand for each tRNA nor the extent to which its concentration changes across environmental conditions has been extensively examined. Here we compute changes in the codon usage of the transcriptome across different conditions in several organisms by inspecting conventional mRNA expression data. We find recurring dynamics of codon usage in the transcriptome in multiple stressful conditions. In particular, codons that are translated by rare tRNAs become over-represented in the transcriptome in response to stresses. These results raise the possibility that the tRNA pool might dynamically change upon stress to support efficient translation of stress-transcribed genes. Alternatively, stress genes may be typically translated with low efficiency, presumably due to lack of sufficient evolutionary optimization pressure on their codon usage.

A self-organized model for cell-differentiation based on variations of molecular decay rates

Systemic properties of living cells are the result of molecular dynamics governed by so-called genetic regulatory networks (GRN). These networks capture all possible features of cells and are responsible for the immense levels of adaptation characteristic to living systems. At any point in time only small subsets of these networks are active. Any active subset of the GRN leads to the expression of particular sets of molecules (expression modes). The subsets of active networks change over time, leading to the observed complex dynamics of expression patterns. Understanding of these dynamics becomes increasingly important in systems biology and medicine. While the importance of transcription rates and catalytic interactions has been widely recognized in modeling genetic regulatory systems, the understanding of the role of degradation of biochemical agents (mRNA, protein) in regulatory dynamics remains limited. Recent experimental data suggests that there exists a functional relation between mRNA and protein decay rates and expression modes. In this paper we propose a model for the dynamics of successions of sequences of active subnetworks of the GRN. The model is able to reproduce key characteristics of molecular dynamics, including homeostasis, multi-stability, periodic dynamics, alternating activity, differentiability, and self-organized critical dynamics. Moreover the model allows to naturally understand the mechanism behind the relation between decay rates and expression modes. The model explains recent experimental observations that decay-rates (or turnovers) vary between differentiated tissue-classes at a general systemic level and highlights the role of intracellular decay rate control mechanisms in cell differentiation.

Functional and topological characterization of transcriptional cooperativity in yeast

Background

Many cellular programs are regulated through the integration of specific transcriptional signals originated from external stimuli, being cooperation between transcription factors a key feature in this process. In this work, we studied how transcriptional cooperativity in yeast is aimed at integrating different regulatory inputs rather than controlling particular cellular functions from a organizational, evolutionary and functional point of view.

Findings

Our results showed that cooperative transcription factor pairs co-evolve and are essential for the life of the cell. When organized into a layered regulatory network, we observed that cooperative transcription factors were preferentially placed in the middle layers, which highlights a role in regulatory signal integration. We also observed significant co-activity and co-evolution between members of the same cooperative pairs, but a lack of common co-expression profile.

Conclusions

Our results suggest that transcriptional cooperativity has a specific role within the regulatory control scheme of the cell, focused in the amplification and integration of cellular signals rather than control of particular cellular functions. This information can be used for better characterization of regulatory interactions between transcription factors, aimed at determining the spatial and temporal control of gene expression.

Cuf2 is a novel meiosis-specific regulatory factor of meiosis maturation

Background

Meiosis is the specialized form of the cell cycle by which diploid cells produce the haploid gametes required for sexual reproduction. Initiation and progression through meiosis requires that the expression of the meiotic genes is precisely controlled so as to provide the correct gene products at the correct times. During meiosis, four temporal gene clusters are either induced or repressed by a cascade of transcription factors.

Principal Findings

In this report a novel copper-fist-type regulator, Cuf2, is shown to be expressed exclusively during meiosis. The expression profile of the cuf2⁺ mRNA revealed that it was induced during middle-phase meiosis. Both cuf2⁺ mRNA and protein levels are unregulated by copper addition or starvation. The transcription of cuf2⁺ required the presence of a functional mei4⁺ gene encoding a key transcription factor that activates the expression of numerous middle meiotic genes. Microscopic analyses of cells expressing a functional Cuf2-GFP protein revealed that Cuf2 co-localized with both homologous chromosomes and sister chromatids during the meiotic divisions. Cells lacking Cuf2 showed an elevated and sustained expression of several of the middle meiotic genes that persisted even during late meiosis. Moreover, cells carrying disrupted cuf2Δ/cuf2Δ alleles displayed an abnormal morphology of the forespore membranes and a dramatic reduction of spore viability.

Significance

Collectively, the results revealed that Cuf2 functions in the timely repression of the middle-phase genes during meiotic differentiation.

Computational inference of mRNA stability from histone modification and transcriptome profiles

Histone modifications play important roles in regulating eukaryotic gene expression and have been used to model expression levels. Here, we present a regression model to systematically infer mRNA stability by comparing transcriptome profiles with ChIP-seq of H3K4me3, H3K27me3 and H3K36me3. The results from multiple human and mouse cell lines show that the inferred unstable mRNAs have significantly longer 3′Untranslated Regions (UTRs) and more microRNA binding sites within 3′UTR than the inferred stable mRNAs. Regression residuals derived from RNA-seq, but not from GRO-seq, are highly correlated with the half-lives measured by pulse-labeling experiments, supporting the rationale of our inference. Whereas, the functions enriched in the inferred stable and unstable mRNAs are consistent with those from pulse-labeling experiments, we found the unstable mRNAs have higher cell-type specificity under functional constraint. We conclude that the systematical use of histone modifications can differentiate non-expressed mRNAs from unstable mRNAs, and distinguish stable mRNAs from highly expressed ones. In summary, we represent the first computational model of mRNA stability inference that compares transcriptome and epigenome profiles, and provides an alternative strategy for directing experimental measurements.

Comparative dynamic transcriptome analysis (cDTA) reveals mutual feedback between mRNA synthesis and degradation

To monitor eukaryotic mRNA metabolism, we developed comparative dynamic transcriptome analysis (cDTA). cDTA provides absolute rates of mRNA synthesis and decay in Saccharomyces cerevisiae (Sc) cells with the use of Schizosaccharomyces pombe (Sp) as an internal standard. cDTA uses nonperturbing metabolic labeling that supersedes conventional methods for mRNA turnover analysis. cDTA reveals that Sc and Sp transcripts that encode orthologous proteins have similar synthesis rates, whereas decay rates are fivefold lower in Sp, resulting in similar mRNA concentrations despite the larger Sp cell volume. cDTA of Sc mutants reveals that a eukaryote can buffer mRNA levels. Impairing transcription with a point mutation in RNA polymerase (Pol) II causes decreased mRNA synthesis rates as expected, but also decreased decay rates. Impairing mRNA degradation by deleting deadenylase subunits of the Ccr4–Not complex causes decreased decay rates as expected, but also decreased synthesis rates. Extended kinetic modeling reveals mutual feedback between mRNA synthesis and degradation that may be achieved by a factor that inhibits synthesis and enhances degradation.

Learning transcriptional regulation on a genome scale: a theoretical analysis based on gene expression data

The recent advent of high-throughput microarray data has enabled the global analysis of the transcriptome, driving the development and application of computational approaches to study transcriptional regulation on the genome scale, by reconstructing in silico the regulatory interactions of the gene network. Although there are many in-depth reviews of such 'reverse-engineering' methodologies, most have focused on the practical aspect of data mining, and few on the biological problem and the biological relevance of the methodology. Therefore, in this review, from a biological perspective, we used a set of yeast microarray data as a working example, to evaluate the fundamental assumptions implicit in associating transcription factor (TF)-target gene expression levels and estimating TFs' activity, and further explore cooperative models. Finally we confirm that the detailed transcription mechanism is overly-complex for expression data alone to reveal, nevertheless, future network reconstruction studies could benefit from the incorporation of context-specific information, the modeling of multiple layers of regulation (e.g. micro-RNA), or the development of approaches for context-dependent analysis, to uncover the mechanisms of gene regulation.

Measurement of genome-wide RNA synthesis and decay rates with Dynamic Transcriptome Analysis (DTA)

Standard transcriptomics measures total cellular RNA levels. Our understanding of gene regulation would be greatly improved if we could measure RNA synthesis and decay rates on a genome-wide level. To that end, the Dynamic Transcriptome Analysis (DTA) method has been developed. DTA combines metabolic RNA labeling with standard transcriptomics to measure RNA synthesis and decay rates in a precise and non-perturbing manner. Here, we present the open source R/Bioconductor software package DTA. It implements all required bioinformatics steps that allow the accurate absolute quantification and comparison of RNA turnover.

Repression of meiotic genes by antisense transcription and by Fkh2 transcription factor in Schizosaccharomyces pombe

In S. pombe, about 5% of genes are meiosis-specific and accumulate little or no mRNA during vegetative growth. Here we use Affymetrix tiling arrays to characterize transcripts in vegetative and meiotic cells. In vegetative cells, many meiotic genes, especially those induced in mid-meiosis, have abundant antisense transcripts. Disruption of the antisense transcription of three of these mid-meiotic genes allowed vegetative sense transcription. These results suggest that antisense transcription represses sense transcription of meiotic genes in vegetative cells. Although the mechanism(s) of antisense mediated transcription repression need to be further explored, our data indicates that RNAi machinery is not required for repression. Previously, we and others used non-strand specific methods to study splicing regulation of meiotic genes and concluded that 28 mid-meiotic genes are spliced only in meiosis. We now demonstrate that the "unspliced" signal in vegetative cells comes from the antisense RNA, not from unspliced sense RNA, and we argue against the idea that splicing regulates these mid-meiotic genes. Most of these mid-meiotic genes are induced in mid-meiosis by the forkhead transcription factor Mei4. Interestingly, deletion of a different forkhead transcription factor, Fkh2, allows low levels of sense expression of some mid-meiotic genes in vegetative cells. We propose that vegetative expression of mid-meiotic genes is repressed at least two independent ways: antisense transcription and Fkh2 repression.

Programmed fluctuations in sense/antisense transcript ratios drive sexual differentiation in S. pombe

Strand-specific RNA sequencing of S. pombe reveals a highly structured programme of ncRNA expression at over 600 loci. Functional investigations show that this extensive ncRNA landscape controls the complex programme of sexual differentiation in S. pombe.

The model eukaryote S. pombe features substantial numbers of ncRNAs many of which are antisense regulatory transcripts (ARTs), ncRNAs expressed on the opposing strand to coding sequences.

Individual ARTs are generated during the mitotic cycle, or at discrete stages of sexual differentiation to downregulate the levels of proteins that drive and coordinate sexual differentiation.

Antisense transcription occurring from events such as bidirectional transcription is not simply artefactual ‘chatter', it performs a critical role in regulating gene expression.

Regulation of the RNA profile is a principal control driving sexual differentiation in the fission yeast Schizosaccharomyces pombe. Before transcription, RNAi-mediated formation of heterochromatin is used to suppress expression, while post-transcription, regulation is achieved via the active stabilisation or destruction of transcripts, and through at least two distinct types of splicing control (Mata et al, 2002; Shimoseki and Shimoda, 2001; Averbeck et al, 2005; Mata and Bähler, 2006; Xue-Franzen et al, 2006; Moldon et al, 2008; Djupedal et al, 2009; Amorim et al, 2010; Grewal, 2010; Cremona et al, 2011).

Around 94% of the S. pombe genome is transcribed (Wilhelm et al, 2008). While many of these transcripts encode proteins (Wood et al, 2002; Bitton et al, 2011), the majority have no known function. We used a strand-specific protocol to sequence total RNA extracts taken from vegetatively growing cells, and at different points during a time course of sexual differentiation. The resulting data redefined existing gene coordinates and identified additional transcribed loci. The frequency of reads at each of these was used to monitor transcript abundance.

Transcript levels at 6599 loci changed in at least one sample (G-statistic; False Discovery Rate <5%). 4231 (72.3%), of which 4011 map to protein-coding genes, while 809 loci were antisense to a known gene. Comparisons between haploid and diploid strains identified changes in transcript levels at over 1000 loci.

At 354 loci, greater antisense abundance was observed relative to sense, in at least one sample (putative antisense regulatory transcripts—ARTs). Since antisense mechanisms are known to modulate sense transcript expression through a variety of inhibitory mechanisms (Faghihi and Wahlestedt, 2009), we postulated that the waves of antisense expression activated at different stages during meiosis might be regulating protein expression.

To ask whether transcription factors that drive sense-transcript levels influenced ART production, we performed RNA-seq of a pat1.114 diploid meiosis in the absence of the transcription factors Atf21 and Atf31 (responsible for late meiotic transcription; Mata et al, 2002). Transcript levels at 185 ncRNA loci showed significant changes in the knockout backgrounds. Although meiotic progression is largely unaffected by removal of Atf21 and Atf31, viability of the resulting spores was significantly diminished, indicating that Atf21- and Atf31-mediated events are critical to efficient sexual differentiation.

If changes to relative antisense/sense transcript levels during a particular phase of sexual differentiation were to regulate protein expression, then the continued presence of the antisense at points in the differentiation programme where it would normally be absent should abolish protein function during this phase. We tested this hypothesis at four loci representing the three means of antisense production: convergent gene expression, improper termination and nascent transcription from an independent locus. Induction of the natural antisense transcripts that opposed spo4⁺, spo6⁺ and dis1⁺ (Figures 3 and 7) in trans from a heterologous locus phenocopied a loss of function of the target protein. ART overexpression decreased Dis1 protein levels. Antisense transcription opposing spk1⁺ originated from improper termination of the sense ups1⁺ transcript on the opposite strand (Figure 3B, left locus). Expression of either the natural full-length ups1⁺ transcript or a truncated version, restricted to the portion of ups1⁺ overlapping spk1⁺ (Figure 3, orange transcripts) in trans from a heterologous locus phenocopied the spk1.Δ differentiation deficiency. Convergent transcription from a neighbouring gene on the opposing strand is, therefore, an effective mechanism to generate RNAi-mediated (below) silencing in fission yeast. Further analysis of the data revealed, for many loci, substantial changes in UTR length over the course of meiosis, suggesting that UTR dynamics may have an active role in regulating gene expression by controlling the transcriptional overlap between convergent adjacent gene pairs.

The RNAi machinery (Grewal, 2010) was required for antisense suppression at each of the dis1, spk1, spo4 and spo6 loci, as antisense to each locus had no impact in ago1.Δ, dcr1.Δ and rdp1.Δ backgrounds. We conclude that RNAi control has a key role in maintaining the fidelity of sexual differentiation in fission yeast. The histone H3 methyl transferase Clr4 was required for antisense control from a heterologous locus.

Thus, a significant portion of the impact of ncRNA upon sexual differentiation arises from antisense gene silencing. Importantly, in contrast to the extensively characterised ability of the RNAi machinery to operate in cis at a target locus in S. pombe (Grewal, 2010), each case of gene silencing generated here could be achieved in trans by expression of the antisense transcript from a single heterologous locus elsewhere in the genome.

Integration of an antibiotic marker gene immediately downstream of the dis1⁺ locus instigated antisense control in an orientation-dependent manner. PCR-based gene tagging approaches are widely used to fuse the coding sequences of epitope or protein tags to a gene of interest. Not only do these tagging approaches disrupt normal 3′UTR controls, but the insertion of a heterologous marker gene immediately downstream of an ORF can clearly have a significant impact upon transcriptional control of the resulting fusion protein. Thus, PCR tagging approaches can no longer be viewed as benign manipulations of a locus that only result in the production of a tagged protein product.

Repression of Dis1 function by gene deletion or antisense control revealed a key role this conserved microtubule regulator in driving the horsetail nuclear migrations that promote recombination during meiotic prophase.

Non-coding transcripts have often been viewed as simple ‘chatter', maintained solely because evolutionary pressures have not been strong enough to force their elimination from the system. Our data show that phenomena such as improper termination and bidirectional transcription are not simply interesting artifacts arising from the complexities of transcription or genome history, but have a critical role in regulating gene expression in the current genome. Given the widespread use of RNAi, it is reasonable to anticipate that future analyses will establish ARTs to have equal importance in other organisms, including vertebrates.

These data highlight the need to modify our concept of a gene from that of a spatially distinct locus. This view is becoming increasingly untenable. Not only are the 5′ and 3′ ends of many genes indistinct, but that this lack of a hard and fast boundary is actively used by cells to control the transcription of adjacent and overlapping loci, and thus to regulate critical events in the life of a cell.

Strand-specific RNA sequencing of S. pombe revealed a highly structured programme of ncRNA expression at over 600 loci. Waves of antisense transcription accompanied sexual differentiation. A substantial proportion of ncRNA arose from mechanisms previously considered to be largely artefactual, including improper 3′ termination and bidirectional transcription. Constitutive induction of the entire spk1⁺, spo4⁺, dis1⁺ and spo6⁺ antisense transcripts from an integrated, ectopic, locus disrupted their respective meiotic functions. This ability of antisense transcripts to disrupt gene function when expressed in trans suggests that cis production at native loci during sexual differentiation may also control gene function. Consistently, insertion of a marker gene adjacent to the dis1⁺ antisense start site mimicked ectopic antisense expression in reducing the levels of this microtubule regulator and abolishing the microtubule-dependent ‘horsetail' stage of meiosis. Antisense production had no impact at any of these loci when the RNA interference (RNAi) machinery was removed. Thus, far from being simply ‘genome chatter', this extensive ncRNA landscape constitutes a fundamental component in the controls that drive the complex programme of sexual differentiation in S. pombe.

Widespread cotranslational formation of protein complexes

Most cellular processes are conducted by multi-protein complexes. However, little is known about how these complexes are assembled. In particular, it is not known if they are formed while one or more members of the complexes are being translated (cotranslational assembly). We took a genomic approach to address this question, by systematically identifying mRNAs associated with specific proteins. In a sample of 31 proteins from Schizosaccharomyces pombe that did not contain RNA–binding domains, we found that ∼38% copurify with mRNAs that encode interacting proteins. For example, the cyclin-dependent kinase Cdc2p associates with the rum1 and cdc18 mRNAs, which encode, respectively, an inhibitor of Cdc2p kinase activity and an essential regulator of DNA replication. Both proteins interact with Cdc2p and are key cell cycle regulators. We obtained analogous results with proteins with different structures and cellular functions (kinesins, protein kinases, transcription factors, proteasome components, etc.). We showed that copurification of a bait protein and of specific mRNAs was dependent on the presence of the proteins encoded by the interacting mRNAs and on polysomal integrity. These results indicate that these observed associations reflect the cotranslational interaction between the bait and the nascent proteins encoded by the interacting mRNAs. Therefore, we show that the cotranslational formation of protein–protein interactions is a widespread phenomenon.

Most proteins do not function in isolation. Instead, they associate with other proteins to form complexes. Little is known about the assembly of protein complexes within cells. One possibility is that proteins are completely synthesised before they bind to each other. An alternative is that proteins attach to each other as they are being translated in the ribosome (called cotranslational assembly). To investigate if cells use cotranslational assembly to form complexes, we identified mRNAs associated with specific proteins. The expectation is that if protein A binds to protein B as protein B is being translated, A will associate indirectly to the mRNA encoding B. Indeed, we found that for ∼40% of proteins (out of a sample of over 30) this was the case. Proteins associated with a small number of mRNAs, most of which encoded known or predicted interacting proteins. We found examples of this phenomenon in proteins with different functions and structures, indicating that cotranslational assembly is widespread. Cotranslational assembly might be required for certain proteins to associate, or it might be important in cases where the early formation of a protein complex is beneficial, such as when a protein is toxic or unstable unless bound to a partner.

Transcriptome kinetics is governed by a genome-wide coupling of mRNA production and degradation: a role for RNA Pol II

Transcriptome dynamics is governed by two opposing processes, mRNA production and degradation. Recent studies found that changes in these processes are frequently coordinated and that the relationship between them shapes transcriptome kinetics. Specifically, when transcription changes are counter-acted with changes in mRNA stability, transient fast-relaxing transcriptome kinetics is observed. A possible molecular mechanism underlying such coordinated regulation might lay in two RNA polymerase (Pol II) subunits, Rpb4 and Rpb7, which are recruited to mRNAs during transcription and later affect their degradation in the cytoplasm. Here we used a yeast strain carrying a mutant Pol II which poorly recruits these subunits. We show that this mutant strain is impaired in its ability to modulate mRNA stability in response to stress. The normal negative coordinated regulation is lost in the mutant, resulting in abnormal transcriptome profiles both with respect to magnitude and kinetics of responses. These results reveal an important role for Pol II, in regulation of both mRNA synthesis and degradation, and also in coordinating between them. We propose a simple model for production-degradation coupling that accounts for our observations. The model shows how a simple manipulation of the rates of co-transcriptional mRNA imprinting by Pol II may govern genome-wide transcriptome kinetics in response to environmental changes.

Organisms alter genes expression programs in response to changes in their environment. Such programs can specify fast induction, slow relaxation, oscillations, etc. Conceivably these kinetic outputs may depend on proper orchestration of the various phases of gene expression, including transcription, translation, and mRNA decay. In particular, in the transcriptomes of a broad range of species, fast mRNA “spikes” appear to result from surprisingly “pressing the gas and the brakes” together, i.e. by activating both transcription and degradation of same transcripts. A recently discovered molecular mechanism, in which subunits of RNA polymerase II (Pol II) associate to mRNAs during transcription and control their decay, could explain how such transcription-decay counter-action works. Yet, how such potential coupling responds to physiological conditions and how it shapes transcriptome kinetics remain unknown. Here we used a minimalist mutation in yeast RNA Pol II that is defective in the above mechanism in order to show that Pol II governs the ability of the cell to modulate mRNA decay in stress and, most importantly, that Pol II is essential for appropriate coupling between mRNA production and degradation. We further show that this transcription-decay coupling is responsible for shaping the transcriptome kinetic profiles under changing environmental conditions.

Coupled evolution of transcription and mRNA degradation

mRNA levels are determined by the balance between transcription and mRNA degradation, and while transcription has been extensively studied, very little is known regarding the regulation of mRNA degradation and its coordination with transcription. Here we examine the evolution of mRNA degradation rates between two closely related yeast species. Surprisingly, we find that around half of the evolutionary changes in mRNA degradation were coupled to transcriptional changes that exert opposite effects on mRNA levels. Analysis of mRNA degradation rates in an interspecific hybrid further suggests that opposite evolutionary changes in transcription and in mRNA degradation are mechanistically coupled and were generated by the same individual mutations. Coupled changes are associated with divergence of two complexes that were previously implicated both in transcription and in mRNA degradation (Rpb4/7 and Ccr4-Not), as well as with sequence divergence of transcription factor binding motifs. These results suggest that an opposite coupling between the regulation of transcription and that of mRNA degradation has shaped the evolution of gene regulation in yeast.

Regulatory mechanisms and networks couple the different phases of gene expression

Gene expression comprises multiple stages, from transcription to protein degradation. Although much is known about the regulation of each stage separately, an understanding of the regulatory coupling between the different stages is only beginning to emerge. For example, there is a clear crosstalk between translation and transcription, and the localization and stability of an mRNA in the cytoplasm could already be determined during transcription in the nucleus. We review a diversity of mechanisms discovered in recent years that couple the different stages of gene expression. We then speculate on the functional and evolutionary significance of this coupling and suggest certain systems-level functionalities that might be optimized via the various coupling modes. In particular, we hypothesize that coupling is often an economic strategy that allows biological systems to respond robustly and precisely to genetic and environmental perturbations.

Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells

Cellular RNA levels are determined by the interplay of RNA production, processing and degradation. However, because most studies of RNA regulation do not distinguish the separate contributions of these processes, little is known about how they are temporally integrated. Here we combine metabolic labeling of RNA at high temporal resolution with advanced RNA quantification and computational modeling to estimate RNA transcription and degradation rates during the response of mouse dendritic cells to lipopolysaccharide. We find that changes in transcription rates determine the majority of temporal changes in RNA levels, but that changes in degradation rates are important for shaping sharp 'peaked' responses. We used sequencing of the newly transcribed RNA population to estimate temporally constant RNA processing and degradation rates genome wide. Degradation rates vary significantly between genes and contribute to the observed differences in the dynamic response. Certain transcripts, including those encoding cytokines and transcription factors, mature faster. Our study provides a quantitative approach to study the integrative process of RNA regulation.