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

Preparation of RNA 3' End Sequencing Libraries of Total and 4-thiouracil Labeled RNA for Simultaneous Measurement of Transcription, RNA Synthesis and Decay in S. cerevisiae

Cellular RNA levels are determined by the rates of RNA transcription from the gene template and subsequent RNA stability. Knowledge about both transcription and RNA decay is, therefore, necessary to interpret RNA levels and gene expression, especially during cellular processes where these parameters change. Numerous experimental strategies have been developed to measure transcription and RNA decay rates. However, to our knowledge, none of those techniques can simultaneously interrogate transcription and RNA decay. The presented protocol allows this and provides a simple approach to simultaneously estimate total RNA levels, transcription and decay rates from the same RNA sample. It is based on brief metabolic labeling of RNA and subsequent concurrent sequencing of polyA⁺ and polyA^- RNA 3' ends. The protocol was developed in S. cerevisiae and should be broadly applicable.

Role of the pre-initiation complex in Mediator recruitment and dynamics

The Mediator complex stimulates the cooperative assembly of a pre-initiation complex (PIC) and recruitment of RNA Polymerase II (Pol II) for gene activation. The core Mediator complex is organized into head, middle, and tail modules, and in budding yeast (Saccharomyces cerevisiae), Mediator recruitment has generally been ascribed to sequence-specific activators engaging the tail module triad of Med2-Med3-Med15 at upstream activating sequences (UASs). We show that yeast lacking Med2-Med3-Med15 are viable and that Mediator and PolII are recruited to promoters genome-wide in these cells, albeit at reduced levels. To test whether Mediator might alternatively be recruited via interactions with the PIC, we examined Mediator association genome-wide after depleting PIC components. We found that depletion of Taf1, Rpb3, and TBP profoundly affected Mediator association at active gene promoters, with TBP being critical for transit of Mediator from UAS to promoter, while Pol II and Taf1 stabilize Mediator association at proximal promoters.

Global analysis of RNA metabolism using bio-orthogonal labeling coupled with next-generation RNA sequencing

Many open questions in RNA biology relate to the kinetics of gene expression and the impact of RNA binding regulatory factors on processing or decay rates of particular transcripts. Steady state measurements of RNA abundance obtained from RNA-seq approaches are not able to separate the effects of transcription from those of RNA decay in the overall abundance of any given transcript, instead only giving information on the (presumed steady-state) abundances of transcripts. Through the combination of metabolic labeling and high-throughput sequencing, several groups have been able to measure both transcription rates and decay rates of the entire transcriptome of an organism in a single experiment. This review focuses on the methodology used to specifically measure RNA decay at a global level. By comparing and contrasting approaches and describing the experimental protocols in a modular manner, we intend to provide both experienced and new researchers to the field the ability to combine aspects of various protocols to fit the unique needs of biological questions not addressed by current methods.

Gaining insight into transcriptome-wide RNA population dynamics through the chemistry of 4-thiouridine

Cellular RNA levels are the result of a juggling act between RNA transcription, processing, and degradation. By tuning one or more of these parameters, cells can rapidly alter the available pool of transcripts in response to stimuli. While RNA sequencing (RNA-seq) is a vital method to quantify RNA levels genome-wide, it is unable to capture the dynamics of different RNA populations at steady-state or distinguish between different mechanisms that induce changes to the steady-state (i.e., altered rate of transcription vs. degradation). The dynamics of different RNA populations can be studied by targeted incorporation of noncanonical nucleosides. 4-Thiouridine (s⁴ U) is a commonly used and versatile RNA metabolic label that allows the study of many properties of RNA metabolism from synthesis to degradation. Numerous experimental strategies have been developed that leverage the power of s⁴ U to label newly transcribed RNA in whole cells, followed by enrichment with activated disulfides or chemistry to induce C mutations at sites of s⁴ U during sequencing. This review presents existing methods to study RNA population dynamics genome-wide using s⁴ U metabolic labeling, as well as a discussion of considerations and challenges when designing s⁴ U metabolic labeling experiments. This article is categorized under: RNA Methods > RNA Analyses in Cells RNA Turnover and Surveillance > Regulation of RNA Stability.

Saccharomyces cerevisiae Metabolic Labeling with 4-thiouracil and the Quantification of Newly Synthesized mRNA As a Proxy for RNA Polymerase II Activity

Global defects in RNA polymerase II transcription might be overlooked by transcriptomic studies analyzing steady-state RNA. Indeed, the global decrease in mRNA synthesis has been shown to be compensated by a simultaneous decrease in mRNA degradation to restore normal steady-state levels. Hence, the genome-wide quantification of mRNA synthesis, independently from mRNA decay, is the best direct reflection of RNA polymerase II transcriptional activity. Here, we discuss a method using non-perturbing metabolic labeling of nascent RNAs in Saccharomyces cerevisiae (S. cerevisiae). Specifically, the cells are cultured for 6 min with a uracil analog, 4-thiouracil, and the labeled newly transcribed RNAs are purified and quantified to determine the synthesis rates of all individual mRNA. Moreover, using labeled Schizosaccharomyces pombe cells as internal standard allows comparing mRNA synthesis in different S. cerevisiae strains. Using this protocol and fitting the data with a dynamic kinetic model, the corresponding mRNA decay rates can be determined.

Global role for coactivator complexes in RNA polymerase II transcription

SAGA and TFIID are related transcription complexes, which were proposed to alternatively deliver TBP at different promoter classes. Recent genome-wide studies in yeast revealed that both complexes are required for the transcription of a vast majority of genes by RNA polymerase II raising new questions about the role of coactivators.

Contributions of regulated transcription and mRNA decay to the dynamics of gene expression

Organisms have acquired sophisticated regulatory networks that control gene expression in response to cellular perturbations. Understanding of the mechanisms underlying the coordinated changes in gene expression in response to external and internal stimuli is a fundamental issue in biology. Recent advances in high-throughput technologies have enabled the measurement of diverse biological information, including gene expression levels, kinetics of gene expression, and interactions among gene expression regulatory molecules. By coupling these technologies with quantitative modeling, we can now uncover the biological roles and mechanisms of gene regulation at the system level. This review consists of two parts. First, we focus on the methods using uridine analogs that measure synthesis and decay rates of RNAs, which demonstrate how cells dynamically change the regulation of gene expression in response to both internal and external cues. Second, we discuss the underlying mechanisms of these changes in kinetics, including the functions of transcription factors and RNA-binding proteins. Overall, this review will help to clarify a system-level view of gene expression programs in cells. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Turnover and Surveillance > Regulation of RNA Stability RNA Methods > RNA Analyses in vitro and In Silico.

Non-invasive measurement of mRNA decay reveals translation initiation as the major determinant of mRNA stability

The cytoplasmic abundance of mRNAs is strictly controlled through a balance of production and degradation. Whereas the control of mRNA synthesis through transcription has been well characterized, less is known about the regulation of mRNA turnover, and a consensus model explaining the wide variations in mRNA decay rates remains elusive. Here, we combine non-invasive transcriptome-wide mRNA production and stability measurements with selective and acute perturbations to demonstrate that mRNA degradation is tightly coupled to the regulation of translation, and that a competition between translation initiation and mRNA decay -but not codon optimality or elongation- is the major determinant of mRNA stability in yeast. Our refined measurements also reveal a remarkably dynamic transcriptome with an average mRNA half-life of only 4.8 min - much shorter than previously thought. Furthermore, global mRNA destabilization by inhibition of translation initiation induces a dose-dependent formation of processing bodies in which mRNAs can decay over time.

A Nuclear Export Block Triggers the Decay of Newly Synthesized Polyadenylated RNA

Genomes are promiscuously transcribed, necessitating mechanisms that facilitate the sorting of RNA for function or destruction. The polyA (pA) tail is one such distinguishing feature, which in the Saccharomyces cerevisiae nucleus is bound by the Nab2p protein, yielding transcript protection. As Nab2p also contacts the main nuclear export factor Mex67p, we asked whether transport kinetics contributes to RNA sorting. Indeed, 3' end sequencing of newly transcribed pA+ RNAs demonstrates that nuclear depletion of Mex67p elicits their instant and global decay. A similar phenotype is evident upon inactivation of other export factors and proportional to the amount of nuclear pA+ RNA. As RNA expression is partially rescued by Nab2p overexpression, we propose that an export block out-titrates Nab2p onto nuclear-retained pA+ RNA, reducing the pool of Nab2p available to protect new transcripts. More generally, we suggest that nuclear RNA decay, negotiated by Nab2p availability, aids in balancing cellular transcript supply with demand.

Distinct patterns of histone acetyltransferase and Mediator deployment at yeast protein-coding genes

Here, Bruzzone et al. explore the relative contributions of the transcriptional coactivators Mediator and two histone acetyltransferase (HAT) complexes, NuA4 and SAGA, to RNA polymerase II association at specific genes and gene classes by rapid nuclear depletion of key complex subunits. They show that the NuA4 HAT Esa1 differentially affects certain groups of genes, whereas the SAGA HAT Gcn5 has a weaker but more uniform effect, and their findings suggest that at least three distinct combinations of coactivator deployment are used to generate moderate or high transcription levels.

The transcriptional coactivators Mediator and two histone acetyltransferase (HAT) complexes, NuA4 and SAGA, play global roles in transcriptional activation. Here we explore the relative contributions of these factors to RNA polymerase II association at specific genes and gene classes by rapid nuclear depletion of key complex subunits. We show that the NuA4 HAT Esa1 differentially affects certain groups of genes, whereas the SAGA HAT Gcn5 has a weaker but more uniform effect. Relative dependence on Esa1 and Tra1, a shared component of NuA4 and SAGA, distinguishes two large groups of coregulated growth-promoting genes. In contrast, we show that the activity of Mediator is particularly important at a separate, small set of highly transcribed TATA-box-containing genes. Our analysis indicates that at least three distinct combinations of coactivator deployment are used to generate moderate or high transcription levels and suggests that each may be associated with distinct forms of regulation.

The conserved RNA recognition motif and C3H1 domain of the Not4 ubiquitin ligase regulate in vivo ligase function

The Ccr4-Not complex controls RNA polymerase II (Pol II) dependent gene expression and proteasome function. The Not4 ubiquitin ligase is a Ccr4-Not subunit that has both a RING domain and a conserved RNA recognition motif and C3H1 domain (referred to as the RRM-C domain) with unknown function. We demonstrate that while individual Not4 RING or RRM-C mutants fail to replicate the proteasomal defects found in Not4 deficient cells, mutation of both exhibits a Not4 loss of function phenotype. Transcriptome analysis revealed that the Not4 RRM-C affects a specific subset of Pol II-regulated genes, including those involved in transcription elongation, cyclin-dependent kinase regulated nutrient responses, and ribosomal biogenesis. The Not4 RING, RRM-C, or RING/RRM-C mutations cause a generalized increase in Pol II binding at a subset of these genes, yet their impact on gene expression does not always correlate with Pol II recruitment which suggests Not4 regulates their expression through additional mechanisms. Intriguingly, we find that while the Not4 RRM-C is dispensable for Ccr4-Not association with RNA Pol II, the Not4 RING domain is required for these interactions. Collectively, these data elucidate previously unknown roles for the conserved Not4 RRM-C and RING domains in regulating Ccr4-Not dependent functions in vivo.

Stochastic system identification without an a priori chosen kinetic model-exploring feasible cell regulation with piecewise linear functions

Kinetic models are at the heart of system identification. A priori chosen rate functions may, however, be unfitting or too restrictive for complex or previously unanticipated regulation. We applied general purpose piecewise linear functions for stochastic system identification in one dimension using published flow cytometry data on E.coli and report on identification results for equilibrium state and dynamic time series. In metabolic labelling experiments during yeast osmotic stress response, we find mRNA production and degradation to be strongly co-regulated. In addition, mRNA degradation appears overall uncorrelated with mRNA level. Comparison of different system identification approaches using semi-empirical synthetic data revealed the superiority of single-cell tracking for parameter identification. Generally, we find that even within restrictive error bounds for deviation from experimental data, the number of viable regulation types may be large. Indeed, distinct regulation can lead to similar expression behaviour over time. Our results demonstrate that molecule production and degradation rates may often differ from classical constant, linear or Michaelis–Menten type kinetics.

Classical cell-regulation models are often imperfectly fitting or even inconsistent with experimental data suggesting inappropriate model assumptions. Martin Hoffmann from Fraunhofer ITEM Regensburg and Jörg Galle from IZBI Leipzig analysed different protein and gene expression data using general purpose piecewise linear functions for system identification. They assessed data corresponding to various experimental techniques for their potential to determine the parameters of their models. Single-cell recordings of expression values over time were most effective for parameter identification. Generally, different and often non-classical cell-regulation models were consistent with the experimental data, even for restrictive error bounds. The authors used virtual treatment experiments to demonstrate that precise knowledge of cell regulation is important for assessing therapy effects. Their findings clearly argue in favour of system identification being performed without an a priori chosen kinetic model.

Cell size control and gene expression homeostasis in single-cells

Growth of a cell and its subsequent division into daughters is a fundamental aspect of all cellular living systems. During these processes, how do individual cells correct size aberrations so that they do not grow abnormally large or small? How do cells ensure that the concentration of essential gene products are maintained at desired levels, in spite of dynamic/stochastic changes in cell size during growth and division? Both these questions have fascinated researchers for over a century. We review how advances in singe-cell technologies and measurements are providing unique insights into these questions across organisms from prokaryotes to human cells. More specifically, diverse strategies based on timing of cell-cycle events, regulating growth, and number of daughters are employed to maintain cell size homeostasis. Interestingly, size homeostasis often results in size optimality - proliferation of individual cells in a population is maximized at an optimal cell size. We further discuss how size-dependent expression or gene-replication timing can buffer concentration of a gene product from cell-to-cell size variations within a population. Finally, we speculate on an intriguing hypothesis that specific size control strategies may have evolved as a consequence of gene-product concentration homeostasis.

Arabidopsis mRNA decay landscape arises from specialized RNA decay substrates, decapping-mediated feedback, and redundancy

The decay of mRNA plays a vital role in modulating mRNA abundance, which, in turn, influences cellular and organismal processes. In plants and metazoans, three distinct pathways carry out the decay of most cytoplasmic mRNAs: The mRNA decapping complex, which requires the scaffold protein VARICOSE (VCS), removes a protective 5' cap, allowing for 5' to 3' decay via EXORIBONUCLEASE4 (XRN4, XRN1 in metazoans and yeast), and both the exosome and SUPPRESSOR OF VCS (SOV)/DIS3L2 degrade RNAs in the 3' to 5' direction. However, the unique biological contributions of these three pathways, and whether they degrade specialized sets of transcripts, are unknown. In Arabidopsis, the participation of SOV in RNA homeostasis is also unclear, because Arabidopsis sov mutants have a normal phenotype. We carried out mRNA decay analyses in wild-type, sov, vcs, and vcs sov seedlings, and used a mathematical modeling approach to determine decay rates and quantify gene-specific contributions of VCS and SOV to decay. This analysis revealed that VCS (decapping) contributes to decay of 68% of the transcriptome, and, while it initiates degradation of mRNAs with a wide range of decay rates, it especially contributes to decay of short-lived RNAs. Only a few RNAs were clear SOV substrates in that they decayed more slowly in sov mutants. However, 4,506 RNAs showed VCS-dependent feedback in sov that modulated decay rates, and, by inference, transcription, to maintain RNA abundances, suggesting that these RNAs might also be SOV substrates. This feedback was shown to be independent of siRNA activity.

Genome-Wide Mapping of Decay Factor-mRNA Interactions in Yeast Identifies Nutrient-Responsive Transcripts as Targets of the Deadenylase Ccr4

The Ccr4 (carbon catabolite repression 4)-Not complex is a major regulator of stress responses that controls gene expression at multiple levels, from transcription to mRNA decay. Ccr4, a “core” subunit of the complex, is the main cytoplasmic deadenylase in Saccharomyces cerevisiae; however, its mRNA targets have not been mapped on a genome-wide scale. Here, we describe a genome-wide approach, RNA immunoprecipitation (RIP) high-throughput sequencing (RIP-seq), to identify the RNAs bound to Ccr4, and two proteins that associate with it, Dhh1 and Puf5. All three proteins were preferentially bound to lowly abundant mRNAs, most often at the 3′ end of the transcript. Furthermore, Ccr4, Dhh1, and Puf5 are recruited to mRNAs that are targeted by other RNA-binding proteins that promote decay and mRNA transport, and inhibit translation. Although Ccr4-Not regulates mRNA transcription and decay, Ccr4 recruitment to mRNAs correlates better with decay rates, suggesting it imparts greater control over transcript abundance through decay. Ccr4-enriched mRNAs are refractory to control by the other deadenylase complex in yeast, Pan2/3, suggesting a division of labor between these deadenylation complexes. Finally, Ccr4 and Dhh1 associate with mRNAs whose abundance increases during nutrient starvation, and those that fluctuate during metabolic and oxygen consumption cycles, which explains the known genetic connections between these factors and nutrient utilization and stress pathways.

Measuring mRNA Decay in Budding Yeast Using Single Molecule FISH

Cellular mRNA levels are determined by the rates of mRNA synthesis and mRNA decay. Typically, mRNA degradation kinetics are measured on a population of cells that are either chemically treated or genetically engineered to inhibit transcription. However, these manipulations can affect the mRNA decay process itself by inhibiting regulatory mechanisms that govern mRNA degradation, especially if they occur on short time-scales. Recently, single molecule fluorescent in situ hybridization (smFISH) approaches have been implemented to quantify mRNA decay rates in single, unperturbed cells. Here, we provide a step-by-step protocol that allows quantification of mRNA decay in single Saccharomyces cerevisiae using smFISH. Our approach relies on fluorescent labeling of single cytoplasmic mRNAs and nascent mRNAs found at active sites of transcription, coupled with mathematical modeling to derive mRNA half-lives. Commercially available, single-stranded smFISH DNA oligonucleotides (smFISH probes) are used to fluorescently label mRNAs followed by the quantification of cellular and nascent mRNAs using freely available spot detection algorithms. Our method enables quantification of mRNA decay of any mRNA in single, unperturbed yeast cells and can be implemented to quantify mRNA turnover in a variety of cell types as well as tissues.

Impact of Methods on the Measurement of mRNA Turnover

The turnover of the RNA molecules is determined by the rates of transcription and RNA degradation. Several methods have been developed to study RNA turnover since the beginnings of molecular biology. Here we summarize the main methods to measure RNA half-life: transcription inhibition, gene control, and metabolic labelling. These methods were used to detect the cellular activity of the mRNAs degradation machinery, including the exo-ribonuclease Xrn1 and the exosome. On the other hand, the study of the differential stability of mature RNAs has been hampered by the fact that different methods have often yielded inconsistent results. Recent advances in the systematic comparison of different method variants in yeast have permitted the identification of the least invasive methodologies that reflect half-lives the most faithfully, which is expected to open the way for a consistent quantitative analysis of the determinants of mRNA stability.

Thiol-linked alkylation of RNA to assess expression dynamics

Gene expression profiling by high-throughput sequencing reveals qualitative and quantitative changes in RNA species at steady state but obscures the intracellular dynamics of RNA transcription, processing and decay. We developed thiol(SH)-linked alkylation for the metabolic sequencing of RNA (SLAM seq), an orthogonal-chemistry-based RNA sequencing technology that detects 4-thiouridine (s4U) incorporation in RNA species at single-nucleotide resolution. In combination with well-established metabolic RNA labeling protocols and coupled to standard, low-input, high-throughput RNA sequencing methods, SLAM seq enabled rapid access to RNA-polymerase-II-dependent gene expression dynamics in the context of total RNA. We validated the method in mouse embryonic stem cells by showing that the RNA-polymerase-II-dependent transcriptional output scaled with Oct4/Sox2/Nanog-defined enhancer activity, and we provide quantitative and mechanistic evidence for transcript-specific RNA turnover mediated by post-transcriptional gene regulatory pathways initiated by microRNAs and N6-methyladenosine. SLAM seq facilitates the dissection of fundamental mechanisms that control gene expression in an accessible, cost-effective and scalable manner.

SAGA mediates transcription from the TATA-like element independently of Taf1p/TFIID but dependent on core promoter structures in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, core promoters of class II genes contain a TATA element, either a TATA box (TATA[A/T]A[A/T][A/G]) or TATA-like element (1 or 2 bp mismatched version of the TATA box). The TATA element directs the assembly of the preinitiation complex (PIC) to ensure accurate transcriptional initiation. It has been proposed the PIC is assembled by two distinct pathways in which TBP is delivered by TFIID or SAGA, leading to the widely accepted model that these complexes mediate transcription mainly from TATA-like element- or TATA box-containing promoters, respectively. Although both complexes are involved in transcription of nearly all class II genes, it remains unclear how efficiently SAGA mediates transcription from TATA-like element-containing promoters independently of TFIID. We found that transcription from the TATA box-containing AGP1 promoter was greatly stimulated in a Spt3p-dependent manner after inactivation of Taf1p/TFIID. Thus, this promoter provides a novel experimental system in which to evaluate SAGA-mediated transcription from TATA-like element(s). We quantitatively measured transcription from various TATA-like elements in the Taf1p-dependent CYC1 promoter and Taf1p-independent AGP1 promoter. The results revealed that SAGA could mediate transcription from at least some TATA-like elements independently of Taf1p/TFIID, and that Taf1p-dependence or -independence is highly robust with respect to variation of the TATA sequence. Furthermore, chimeric promoter mapping revealed that Taf1p-dependence or independence was conferred by the upstream activating sequence (UAS), whereas Spt3p-dependent transcriptional stimulation after inactivation of Taf1p/TFIID was specific to the AGP1 promoter and dependent on core promoter regions other than the TATA box. These results suggest that TFIID and/or SAGA are regulated in two steps: the UAS first specifies TFIID or SAGA as the predominant factor on a given promoter, and then the core promoter structure guides the pertinent factor to conduct transcription in an appropriate manner.

From transcriptional complexity to cellular phenotypes: Lessons from yeast

Pervasive transcription has been reported in many eukaryotic organisms, revealing a highly interleaved transcriptome organization that involves thousands of coding and non-coding RNAs. However, to date, the biological impact of transcriptome complexity is still poorly understood. Here I will review how subtle variations of the transcriptome can lead to divergent cellular phenotypes by fine-tuning both its coding potential and regulation. I will discuss strategies that can be used to link molecular variations with divergent biological outcomes. Finally, I will explore the implication of transcriptional complexity for our understanding of gene expression in the context of cell-to-cell phenotypic variability. Copyright © 2017 John Wiley & Sons, Ltd.

SAGA Is a General Cofactor for RNA Polymerase II Transcription

Prior studies suggested that SAGA and TFIID are alternative factors that promote RNA polymerase II transcription with about 10% of genes in S. cerevisiae dependent on SAGA. We reassessed the role of SAGA by mapping its genome-wide location and role in global transcription in budding yeast. We find that SAGA maps to the UAS elements of most genes, overlapping with Mediator binding and irrespective of previous designations of SAGA or TFIID-dominated genes. Disruption of SAGA through mutation or rapid subunit depletion reduces transcription from nearly all genes, measured by newly-synthesized RNA. We also find that the acetyltransferase Gcn5 synergizes with Spt3 to promote global transcription and that Spt3 functions to stimulate TBP recruitment at all tested genes. Our data demonstrate that SAGA acts as a general cofactor required for essentially all RNA polymerase II transcription and is not consistent with the previous classification of SAGA and TFIID-dominated genes.

Transcription of Nearly All Yeast RNA Polymerase II-Transcribed Genes Is Dependent on Transcription Factor TFIID

Previous studies suggested that expression of most yeast mRNAs is dominated by either transcription factor TFIID or SAGA. We re-examined the role of TFIID by rapid depletion of S. cerevisiae TFIID subunits and measurement of changes in nascent transcription. We find that transcription of nearly all mRNAs is strongly dependent on TFIID function. Degron-dependent depletion of Taf1, Taf2, Taf7, Taf11, and Taf13 showed similar transcription decreases for genes in the Taf1-depleted, Taf1-enriched, TATA-containing, and TATA-less gene classes. The magnitude of TFIID dependence varies with growth conditions, although this variation is similar genome-wide. Many studies have suggested differences in gene-regulatory mechanisms between TATA and TATA-less genes, and these differences have been attributed in part to differential dependence on SAGA or TFIID. Our work indicates that TFIID participates in expression of nearly all yeast mRNAs and that differences in regulation between these two gene categories is due to other properties.

Phosphorylation and nuclear transit modulate the balance between normal function and terminal aggregation of the yeast RNA-binding protein Ssd1

Ssd1 targets mRNAs to daughter cells for translation or to stress granules (SGs) and P-bodies (PBs) for storage or decay. PB components also assist in its nuclear export. If Ssd1 fails to localize to the nucleus, it is targeted to IPOD. IPOD and PB/SG association requires a prion-like domain, whose activity is differentially regulated by Cbk1 phosphorylation.

Yeast Ssd1 is an RNA-binding protein that shuttles between the nucleus and cytoplasm. Ssd1 interacts with its target mRNAs initially during transcription by binding through its N-terminal prion-like domain (PLD) to the C-terminal domain of RNA polymerase II. Ssd1 subsequently targets mRNAs acquired in the nucleus either to daughter cells for translation or to stress granules (SGs) and P-bodies (PBs) for mRNA storage or decay. Here we show that PB components assist in the nuclear export of Ssd1and subsequent targeting of Ssd1 to PB sites in the cytoplasm. In the absence of import into the nucleus, Ssd1 fails to associate with PBs in the cytoplasm but rather is targeted to cytosolic insoluble protein deposits (IPODs). The association of Ssd1 either with IPOD sites or with PB/SG requires the PLD, whose activity is differentially regulated by the Ndr/LATS family kinase, Cbk1: phosphorylation suppresses PB/SG association but enhances IPOD formation. This regulation likely accrues from a phosphorylation-sensitive nuclear localization sequence located in the PLD. The results presented here may inform our understanding of aggregate formation by RBP in certain neurological diseases.

Comprehensive analysis of nucleocytoplasmic dynamics of mRNA in Drosophila cells

Eukaryotic mRNAs undergo a cycle of transcription, nuclear export, and degradation. A major challenge is to obtain a global, quantitative view of these processes. Here we measured the genome-wide nucleocytoplasmic dynamics of mRNA in Drosophila cells by metabolic labeling in combination with cellular fractionation. By mathematical modeling of these data we determined rates of transcription, export and cytoplasmic decay for 5420 genes. We characterized these kinetic rates and investigated links with mRNA features, RNA-binding proteins (RBPs) and chromatin states. We found prominent correlations between mRNA decay rate and transcript size, while nuclear export rates are linked to the size of the 3'UTR. Transcription, export and decay rates are each associated with distinct spectra of RBPs. Specific classes of genes, such as those encoding cytoplasmic ribosomal proteins, exhibit characteristic combinations of rate constants, suggesting modular control. Binding of splicing factors is associated with faster rates of export, and our data suggest coordinated regulation of nuclear export of specific functional classes of genes. Finally, correlations between rate constants suggest global coordination between the three processes. Our approach provides insights into the genome-wide nucleocytoplasmic kinetics of mRNA and should be generally applicable to other cell systems.

All mRNAs start from production in the nucleus, undergo exportation through nuclear pores and finally are degraded in the cytoplasm. A comprehensive characterization of the kinetic rates of all mRNAs is an important prerequisite for a global understanding of the regulation of the transcriptome and the cell. By conducting a time-series experiment and building a mathematical model, we trace the dynamics of mRNAs from the nucleus to the cytoplasm and determine the rates at each kinetic step at transcriptome-wide level. This information allows us to associate mRNA kinetic rates with a wealth of biological features and made some intriguing discoveries. We show mRNA decay is positively linked to transcript length while mRNA export is negatively linked to the length of the 3' UTR. We show binding of splicing factors is associated with faster rates of mRNA export. We provide evidence for global coordination between nuclear export an decay of mRNA. We show genes sharing specific functions tend to have similar nucleoplasmic kinetics, in which ribosomal proteins possessing special kinetic features exclusively stand out. Altogether, our integrated approach to quantitatively determine the rates of kinetic steps on a gene-by-gene basis provides a blueprint to obtain the global understanding of RNA regulation.

The interplay between transcription and mRNA degradation in Saccharomyces cerevisiae

The cellular transcriptome is shaped by both the rates of mRNA synthesis in the nucleus and mRNA degradation in the cytoplasm under a specified condition. The last decade witnessed an exciting development in the field of post-transcriptional regulation of gene expression which underscored a strong functional coupling between the transcription and mRNA degradation. The functional integration is principally mediated by a group of specialized promoters and transcription factors that govern the stability of their cognate transcripts by “marking” them with a specific factor termed “coordinator.” The “mark” carried by the message is later decoded in the cytoplasm which involves the stimulation of one or more mRNA-decay factors, either directly by the “coordinator” itself or in an indirect manner. Activation of the decay factor(s), in turn, leads to the alteration of the stability of the marked message in a selective fashion. Thus, the integration between mRNA synthesis and decay plays a potentially significant role to shape appropriate gene expression profiles during cell cycle progression, cell division, cellular differentiation and proliferation, stress, immune and inflammatory responses, and may enhance the rate of biological evolution.

Multiplexed gene control reveals rapid mRNA turnover

The rates of mRNA synthesis and decay determine the mRNA expression level. The two processes are under coordinated control, which makes the measurements of these rates challenging, as evidenced by the low correlation among the methods of measurement of RNA half-lives. We developed a minimally invasive method, multiplexed gene control, to shut off expression of genes with controllable synthetic promoters. The method was validated by measuring the ratios of the nascent to mature mRNA molecules and by measuring the half-life with endogenous promoters that can be controlled naturally or through inserting short sequences that impart repressibility. The measured mRNA half-lives correlated highly with those obtained with the metabolic pulse-labeling method in yeast. However, mRNA degradation was considerably faster in comparison to previous estimates, with a median half-life of around 2 min. The half-life permits the estimation of promoter-dependent and promoter-independent transcription rates. The dynamical range of the promoter-independent transcription rates was larger than that of the mRNA half-lives. The rapid mRNA turnover and the broad adjustability of promoter-independent transcription rates are expected to have a major impact on stochastic gene expression and gene network behavior.

Analysis of mRNA deadenylation by multi-protein complexes

Poly(A) tails are found at the 3' end of almost every eukaryotic mRNA and are important for the stability of mRNAs and their translation into proteins. Thus, removal of the poly(A) tail, a process called deadenylation, is critical for regulation of gene expression. Most deadenylation enzymes are components of large multi-protein complexes. Here, we describe an in vitro deadenylation assay developed to study the exonucleolytic activities of the multi-protein Ccr4-Not and Pan2-Pan3 complexes. We discuss how this assay can be used with short synthetic RNAs, as well as longer RNA substrates generated using in vitro transcription. Importantly, quantitation of the reactions allows detailed analyses of deadenylation in the presence and absence of accessory factors, leading to new insights into targeted mRNA decay.

Architecture of the RNA polymerase II-Paf1C-TFIIS transcription elongation complex

The conserved polymerase-associated factor 1 complex (Paf1C) plays multiple roles in chromatin transcription and genomic regulation. Paf1C comprises the five subunits Paf1, Leo1, Ctr9, Cdc73 and Rtf1, and binds to the RNA polymerase II (Pol II) transcription elongation complex (EC). Here we report the reconstitution of Paf1C from Saccharomyces cerevisiae, and a structural analysis of Paf1C bound to a Pol II EC containing the elongation factor TFIIS. Cryo-electron microscopy and crosslinking data reveal that Paf1C is highly mobile and extends over the outer Pol II surface from the Rpb2 to the Rpb3 subunit. The Paf1-Leo1 heterodimer and Cdc73 form opposite ends of Paf1C, whereas Ctr9 bridges between them. Consistent with the structural observations, the initiation factor TFIIF impairs Paf1C binding to Pol II, whereas the elongation factor TFIIS enhances it. We further show that Paf1C is globally required for normal mRNA transcription in yeast. These results provide a three-dimensional framework for further analysis of Paf1C function in transcription through chromatin.

mRNA decay: an adaptation tool for the environmental fungal pathogen Cryptococcus neoformans

Fungi are ubiquitous in the environment and humans constantly encounter them in the soil, air, water, and food. The vast majority of these interactions are inconsequential. However, in the context of immunodeficiency precipitated by HIV infection, hematologic malignancy, or transplantation, a small subset of fungi can cause devastating, systemic infection. The most deadly of the opportunistic environmental fungi, Cryptococcus neoformans, is estimated to cause hundreds of thousands of deaths per year, mostly in the context of HIV co-infection. The cellular processes that mediate adaptation to the host environment are of great interest as potential novel therapeutic targets. One such cellular process important for host adaptation is mRNA decay, which mediates the specific degradation of subsets of functionally related mRNAs in response to stressors relevant to pathogenesis, including human core body temperature, carbon limitation, and reactive oxygen stress. Thus, for C. neoformans, host adaptation requires mRNA decay to mediate rapid transcriptome remodeling in the face of stressors encountered in the host. Several nodes of stress-responsive signaling that govern the stress-responsive transcriptome also control the decay rate of mRNAs cleared from the ribosome during stress, suggesting an additional layer of coupling between mRNA synthesis and decay that allows C. neoformans to be a successful pathogen of humans. WIREs RNA 2017, 8:e1424. doi: 10.1002/wrna.1424 For further resources related to this article, please visit the WIREs website.

The link between adjacent codon pairs and mRNA stability

Background

Evidence in diverse organisms suggests that codon optimality is a major determinant of mRNA translation and degradation. Codon optimality is thought to act by modulating the efficiency of ribosome elongation. In Saccharomyces cerevisiae, a recent study has identified 17 adjacent codon pairs that mediate strong inhibition of translation elongation. However, relationships between the inhibitory codon pairs and other aspects of gene expression are unknown.

Results

To gain insights into how the inhibitory codon pairs may affect aspects of gene expression, we utilized existing datasets to conduct genome-scale analyses in S. cerevisiae. Our analysis revealed the following points. First, the inhibitory codon pairs are significantly associated with faster mRNA decay. The association is not solely due to the content of nucleotides, individual codons, or dipeptides encoded by the inhibitory codon pairs. Second, the inhibitory codon pairs cannot fully explain the previously known relationship of codon optimality with mRNA stability, suggesting that optimality of individual codons and properties of adjacent codon pairs both contribute to gene regulation. Finally, although the inhibitory codon pairs are associated with slower mRNA synthesis and protein instability, the associations can be attributed to usage bias in individual codons.

Conclusions

This study suggests an association of inhibitory codon pairs with mRNA stability and thus another layer of complexity in the codon-mediated gene regulation.

Electronic supplementary material

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

Capturing in vivo RNA transcriptional dynamics from the malaria parasite Plasmodium falciparum

To capture the transcriptional dynamics within proliferating cells, methods to differentiate nascent transcription from preexisting mRNAs are desired. One approach is to label newly synthesized mRNA transcripts in vivo through the incorporation of modified pyrimidines. However, the human malaria parasite, Plasmodium falciparum, is incapable of pyrimidine salvage for mRNA biogenesis. To capture cellular mRNA dynamics during Plasmodium development, we engineered parasites that can salvage pyrimidines through the expression of a single bifunctional yeast fusion gene, cytosine deaminase/uracil phosphoribosyltransferase (FCU). We show that expression of FCU allows for the direct incorporation of thiol-modified pyrimidines into nascent mRNAs. Using developmental stage-specific promoters to express FCU-GFP enables the biosynthetic capture and in-depth analysis of mRNA dynamics from subpopulations of cells undergoing differentiation. We demonstrate the utility of this method by examining the transcriptional dynamics of the sexual gametocyte stage transition, a process that is essential to malaria transmission between hosts. Using the pfs16 gametocyte-specific promoter to express FCU-GFP in 3D7 parasites, we found that sexual stage commitment is governed by transcriptional reprogramming and stabilization of a subset of essential gametocyte transcripts. We also measured mRNA dynamics in F12 gametocyte-deficient parasites and demonstrate that the transcriptional program required for sexual commitment and maturation is initiated but likely aborted due to the absence of the PfAP2-G transcriptional regulator and a lack of gametocyte-specific mRNA stabilization. Biosynthetic labeling of Plasmodium mRNAs is incredibly versatile, can be used to measure transcriptional dynamics at any stage of parasite development, and will allow for future applications to comprehensively measure RNA-protein interactions in the malaria parasite.

Kinetic CRAC uncovers a role for Nab3 in determining gene expression profiles during stress

RNA-binding proteins play a key role in shaping gene expression profiles during stress, however, little is known about the dynamic nature of these interactions and how this influences the kinetics of gene expression. To address this, we developed kinetic cross-linking and analysis of cDNAs (χCRAC), an ultraviolet cross-linking method that enabled us to quantitatively measure the dynamics of protein-RNA interactions in vivo on a minute time-scale. Here, using χCRAC we measure the global RNA-binding dynamics of the yeast transcription termination factor Nab3 in response to glucose starvation. These measurements reveal rapid changes in protein-RNA interactions within 1 min following stress imposition. Changes in Nab3 binding are largely independent of alterations in transcription rate during the early stages of stress response, indicating orthogonal transcriptional control mechanisms. We also uncover a function for Nab3 in dampening expression of stress-responsive genes. χCRAC has the potential to greatly enhance our understanding of in vivo dynamics of protein-RNA interactions.Protein RNA interactions are dynamic and regulated in response to environmental changes. Here the authors describe 'kinetic CRAC', an approach that allows time resolved analyses of protein RNA interactions with minute time point resolution and apply it to gain insight into the function of the RNA-binding protein Nab3.

Identification of a transcriptional activation domain in yeast repressor activator protein 1 (Rap1) using an altered DNA-binding specificity variant

Repressor activator protein 1 (Rap1) performs multiple vital cellular functions in the budding yeast Saccharomyces cerevisiae These include regulation of telomere length, transcriptional repression of both telomere-proximal genes and the silent mating type loci, and transcriptional activation of hundreds of mRNA-encoding genes, including the highly transcribed ribosomal protein- and glycolytic enzyme-encoding genes. Studies of the contributions of Rap1 to telomere length regulation and transcriptional repression have yielded significant mechanistic insights. However, the mechanism of Rap1 transcriptional activation remains poorly understood because Rap1 is encoded by a single copy essential gene and is involved in many disparate and essential cellular functions, preventing easy interpretation of attempts to directly dissect Rap1 structure-function relationships. Moreover, conflicting reports on the ability of Rap1-heterologous DNA-binding domain fusion proteins to serve as chimeric transcriptional activators challenge use of this approach to study Rap1. Described here is the development of an altered DNA-binding specificity variant of Rap1 (Rap1AS). We used Rap1AS to map and characterize a 41-amino acid activation domain (AD) within the Rap1 C terminus. We found that this AD is required for transcription of both chimeric reporter genes and authentic chromosomal Rap1 enhancer-containing target genes. Finally, as predicted for a bona fide AD, mutation of this newly identified AD reduced the efficiency of Rap1 binding to a known transcriptional coactivator TFIID-binding target, Taf5. In summary, we show here that Rap1 contains an AD required for Rap1-dependent gene transcription. The Rap1AS variant will likely also be useful for studies of the functions of Rap1 in other biological pathways.

Genome-wide Analysis of RNA Polymerase II Termination at Protein-Coding Genes

At the end of protein-coding genes, RNA polymerase (Pol) II undergoes a concerted transition that involves 3'-processing of the pre-mRNA and transcription termination. Here, we present a genome-wide analysis of the 3'-transition in budding yeast. We find that the 3'-transition globally requires the Pol II elongation factor Spt5 and factors involved in the recognition of the polyadenylation (pA) site and in endonucleolytic RNA cleavage. Pol II release from DNA occurs in a narrow termination window downstream of the pA site and requires the "torpedo" exonuclease Rat1 (XRN2 in human). The Rat1-interacting factor Rai1 contributes to RNA degradation downstream of the pA site. Defects in the 3'-transition can result in increased transcription at downstream genes.

External calibration with Drosophila whole-cell spike-ins delivers absolute mRNA fold changes from human RNA-Seq and qPCR data

Gene expression measurements are typically performed on a fixed-weight aliquot of RNA, which assumes that the total number of transcripts per cell stays nearly constant across all conditions. In cases where this assumption does not hold (e.g., when comparing cell types with different cell sizes) the expression data provide a distorted view of cellular events. Assuming constant numbers of total transcripts, increases in expression of some RNAs must be compensated for by decreases in expression of others. Therefore, we propose calibrating gene expression data to an external reference point, the number of cells in the sample, using whole-cell spike-ins. In a systematic dilution experiment, we mixed varying numbers of human cells with fixed numbers of Drosophila melanogaster cells and scaled the expression levels of the human genes relative to those of the Drosophila genes. This approach restored the original gene expression ratios generated by the dilutions. We then used Drosophila whole-cell spike-ins to uncover non-symmetric gene expression changes, in this case much larger numbers of induced than repressed genes, under perturbations of the human cell line P493–6. Drosophila whole-cell spike-ins are an experimentally and computationally easy and low-priced method to derive mRNA fold changes of absolute abundances from RNA sequencing (RNA-Seq) and quantitative real-time PCR (qPCR) data.

The CELF1 RNA-Binding Protein Regulates Decay of Signal Recognition Particle mRNAs and Limits Secretion in Mouse Myoblasts

We previously identified several mRNAs encoding components of the secretory pathway, including signal recognition particle (SRP) subunit mRNAs, among transcripts associated with the RNA-binding protein CELF1. Through immunoprecipitation of RNAs crosslinked to CELF1 in myoblasts and in vitro binding assays using recombinant CELF1, we now provide evidence that CELF1 directly binds the mRNAs encoding each of the subunits of the SRP. Furthermore, we determined the half-lives of the Srp transcripts in control and CELF1 knockdown myoblasts. Our results indicate CELF1 is a destabilizer of at least five of the six Srp transcripts and that the relative abundance of the SRP proteins is out of balance when CELF1 is depleted. CELF1 knockdown myoblasts exhibit altered secretion of a luciferase reporter protein and are impaired in their ability to migrate and close a wound, consistent with a defect in the secreted extracellular matrix. Importantly, similar defects in wound healing are observed when SRP subunit imbalance is induced by over-expression of SRP68. Our studies support the existence of an RNA regulon containing Srp mRNAs that is controlled by CELF1. One implication is that altered function of CELF1 in myotonic dystrophy may contribute to changes in the extracellular matrix of affected muscle through defects in secretion.

Molecular mechanisms that distinguish TFIID housekeeping from regulatable SAGA promoters

An important distinction is frequently made between constitutively expressed housekeeping genes versus regulated genes. Although generally characterized by different DNA elements, chromatin architecture and cofactors, it is not known to what degree promoter classes strictly follow regulatability rules and which molecular mechanisms dictate such differences. We show that SAGA‐dominated/TATA‐box promoters are more responsive to changes in the amount of activator, even compared to TFIID/TATA‐like promoters that depend on the same activator Hsf1. Regulatability is therefore an inherent property of promoter class. Further analyses show that SAGA/TATA‐box promoters are more dynamic because TATA‐binding protein recruitment through SAGA is susceptible to removal by Mot1. In addition, the nucleosome configuration upon activator depletion shifts on SAGA/TATA‐box promoters and seems less amenable to preinitiation complex formation. The results explain the fundamental difference between housekeeping and regulatable genes, revealing an additional facet of combinatorial control: an activator can elicit a different response dependent on core promoter class.

Enriching s4 U-RNA Using Methane Thiosulfonate (MTS) Chemistry

Metabolic labeling of cellular RNA is a useful approach to study RNA biology. 4-Thiouridine (s⁴ U) is a convenient nucleoside for metabolic labeling because it is cell permeable and is incorporated into newly transcribed RNA, and the sulfur moiety provides a handle for biochemical purification. However, a critical step in the purification of s⁴ U-RNA is the efficiency of the chemistry used to enrich s⁴ U-RNA. Here, we present a protocol for s⁴ U-RNA enrichment that includes efficient and reversible covalent chemistry to biotinylate s⁴ U-RNA using the activated disulfide methane thiosulfonate conjugated to biotin (MTS-biotin), followed by enrichment on streptavidin beads. The efficiency of this chemistry reduces enrichment bias and requires less starting material, thereby expanding the utility of s⁴ U to study RNA biology. © 2016 by John Wiley & Sons, Inc.

RNase J depletion leads to massive changes in mRNA abundance in Helicobacter pylori

Degradation of RNA as an intermediate message between genes and corresponding proteins is important for rapid attenuation of gene expression and maintenance of cellular homeostasis. This process is controlled by ribonucleases that have different target specificities. In the bacterial pathogen Helicobacter pylori, an exo- and endoribonuclease RNase J is essential for growth. To explore the role of RNase J in H. pylori, we identified its putative targets at a global scale using next generation RNA sequencing. We found that strong depletion for RNase J led to a massive increase in the steady-state levels of non-rRNAs. mRNAs and RNAs antisense to open reading frames were most affected with over 80% increased more than 2-fold. Non-coding RNAs expressed in the intergenic regions were much less affected by RNase J depletion. Northern blotting of selected messenger and non-coding RNAs validated these results. Globally, our data suggest that RNase J of H. pylori is a major RNase involved in degradation of most cellular RNAs.

Measures of RNA metabolism rates: Toward a definition at the level of single bonds

We give an overview of experimental and computational methods to estimate RNA metabolism rates genome-wide. We then advocate a local definition of RNA metabolism rate at the level of individual phosphodiester bonds. Rates of formation and disappearance of individual bonds are unambiguously defined, in contrast to rates of complete transcripts. We show that over previous approaches, the recently developed transient transcriptome sequencing (TT-seq) protocol allows for estimation of metabolism rates of individual bonds with least positional bias.

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.

Inference of gene regulation functions from dynamic transcriptome data

To quantify gene regulation, a function is required that relates transcription factor binding to DNA (input) to the rate of mRNA synthesis from a target gene (output). Such a ‘gene regulation function’ (GRF) generally cannot be measured because the experimental titration of inputs and simultaneous readout of outputs is difficult. Here we show that GRFs may instead be inferred from natural changes in cellular gene expression, as exemplified for the cell cycle in the yeast S. cerevisiae. We develop this inference approach based on a time series of mRNA synthesis rates from a synchronized population of cells observed over three cell cycles. We first estimate the functional form of how input transcription factors determine mRNA output and then derive GRFs for target genes in the CLB2 gene cluster that are expressed during G2/M phase. Systematic analysis of additional GRFs suggests a network architecture that rationalizes transcriptional cell cycle oscillations. We find that a transcription factor network alone can produce oscillations in mRNA expression, but that additional input from cyclin oscillations is required to arrive at the native behaviour of the cell cycle oscillator.

DOI:http://dx.doi.org/10.7554/eLife.12188.001

Living cells rely on networks of genes to control their behavior, including how they grow, develop and respond to stress. Genes encode instructions needed to make proteins and other molecules, and much of the control is exerted at the first stage of protein production, known as transcription. During this process, a gene is copied to make molecules known as transcripts that may later be used as templates to make proteins.

Many genes encode proteins that act to regulate transcription. Therefore, an individual gene may receive inputs from other genes, and these inputs affect how much transcript the gene produces, which can be considered as the gene’s output. While these inputs and outputs can often be wired together to form a network, it is less clear exactly how all the different inputs at a gene interact to determine its output. These interactions are known as “gene regulation functions”, and knowing them would be an important step towards understanding gene networks, which would help us to predict how cells will behave in different situations.

Gene regulation functions are difficult to measure directly, so researchers would like to find other ways to assess them indirectly. A recently developed experimental technique called “dynamic transcriptome analysis” seemed promising as it measures both the inputs and outputs of all genes in a cell over time. Hillenbrand et al. used this technique to infer gene regulation functions with one or two inputs in yeast cells.

Comparing these estimates with experimental data from previous studies showed that these inferred gene regulation functions could successfully predict the output of a gene based on its inputs. Hillenbrand et al. then used these estimates to search and model a well-known genetic network that is thought to be part of the molecular clockwork that controls the timing of events that cause a cell to divide.

Currently, the approach used by Hillenbrand et al. treats gene regulation functions like “black boxes”. This means that, while an output can be predicted if the inputs are known, it cannot reveal all of the detailed mechanisms behind it. Gaining insights into the inner workings of these black boxes will require taking more data into account, such as how abundant the proteins that regulate transcription are, where they are located within cells or whether they are active or not. Therefore, the next challenge is to incorporate these kinds of data to gain a fuller picture of how gene networks operate within cells.

DOI:http://dx.doi.org/10.7554/eLife.12188.002

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.

Engineered Covalent Inactivation of TFIIH-Kinase Reveals an Elongation Checkpoint and Results in Widespread mRNA Stabilization

During transcription initiation, the TFIIH-kinase Kin28/Cdk7 marks RNA polymerase II (Pol II) by phosphorylating the C-terminal domain (CTD) of its largest subunit. Here we describe a structure-guided chemical approach to covalently and specifically inactivate Kin28 kinase activity in vivo. This method of irreversible inactivation recapitulates both the lethal phenotype and the key molecular signatures that result from genetically disrupting Kin28 function in vivo. Inactivating Kin28 impacts promoter release to differing degrees and reveals a "checkpoint" during the transition to productive elongation. While promoter-proximal pausing is not observed in budding yeast, inhibition of Kin28 attenuates elongation-licensing signals, resulting in Pol II accumulation at the +2 nucleosome and reduced transition to productive elongation. Furthermore, upon inhibition, global stabilization of mRNA masks different degrees of reduction in nascent transcription. This study resolves long-standing controversies on the role of Kin28 in transcription and provides a rational approach to irreversibly inhibit other kinases in vivo.

The transcription factor ERG recruits CCR4-NOT to control mRNA decay and mitotic progression

Control of mRNA levels, a fundamental aspect in the regulation of gene expression, is achieved through a balance between mRNA synthesis and decay. E26-related gene (Erg) proteins are canonical transcription factors whose previously described functions are confined to the control of mRNA synthesis. Here, we report that ERG also regulates gene expression by affecting mRNA stability and identify the molecular mechanisms underlying this function in human cells. ERG is recruited to mRNAs via interaction with the RNA-binding protein RBPMS, and it promotes mRNA decay by binding CNOT2, a component of the CCR4-NOT deadenylation complex. Transcriptome-wide mRNA stability analysis revealed that ERG controls the degradation of a subset of mRNAs highly connected to Aurora signaling, whose decay during S phase is necessary for mitotic progression. Our data indicate that control of gene expression by mammalian transcription factors may follow a more complex scheme than previously anticipated, integrating mRNA synthesis and degradation.

Divergence of a conserved elongation factor and transcription regulation in budding and fission yeast

Complex regulation of gene expression in mammals has evolved from simpler eukaryotic systems, yet the mechanistic features of this evolution remain elusive. Here, we compared the transcriptional landscapes of the distantly related budding and fission yeast. We adapted the Precision Run-On sequencing (PRO-seq) approach to map the positions of RNA polymerase active sites genome-wide in Schizosaccharomyces pombe and Saccharomyces cerevisiae. Additionally, we mapped preferred sites of transcription initiation in each organism using PRO-cap. Unexpectedly, we identify a pause in early elongation, specific to S. pombe, that requires the conserved elongation factor subunit Spt4 and resembles promoter-proximal pausing in metazoans. PRO-seq profiles in strains lacking Spt4 reveal globally elevated levels of transcribing RNA Polymerase II (Pol II) within genes in both species. Messenger RNA abundance, however, does not reflect the increases in Pol II density, indicating a global reduction in elongation rate. Together, our results provide the first base-pair resolution map of transcription elongation in S. pombe and identify divergent roles for Spt4 in controlling elongation in budding and fission yeast.

The Modifier of Transcription 1 (Mot1) ATPase and Spt16 Histone Chaperone Co-regulate Transcription through Preinitiation Complex Assembly and Nucleosome Organization

Modifier of transcription 1 (Mot1) is a conserved and essential Swi2/Snf2 ATPase that can remove TATA-binding protein (TBP) from DNA using ATP hydrolysis and in so doing exerts global effects on transcription. Spt16 is also essential and functions globally in transcriptional regulation as a component of the facilitates chromatin transcription (FACT) histone chaperone complex. Here we demonstrate that Mot1 and Spt16 regulate a largely overlapping set of genes in Saccharomyces cerevisiae. As expected, Mot1 was found to control TBP levels at co-regulated promoters. In contrast, Spt16 did not affect TBP recruitment. On a global scale, Spt16 was required for Mot1 promoter localization, and Mot1 also affected Spt16 localization to genes. Interestingly, we found that Mot1 has an unanticipated role in establishing or maintaining the occupancy and positioning of nucleosomes at the 5' ends of genes. Spt16 has a broad role in regulating chromatin organization in gene bodies, including those nucleosomes affected by Mot1. These results suggest that the large scale overlap in Mot1 and Spt16 function arises from a combination of both their unique and shared functions in transcription complex assembly and chromatin structure regulation.

Translational Capacity of a Cell Is Determined during Transcription Elongation via the Ccr4-Not Complex

The current understanding of gene expression considers transcription and translation to be independent processes. Challenging this notion, we found that translation efficiency is determined during transcription elongation through the imprinting of mRNAs with Not1, the central scaffold of the Ccr4-Not complex. We determined that another subunit of the complex, Not5, defines Not1 binding to specific mRNAs, particularly those produced from ribosomal protein genes. This imprinting mechanism specifically regulates ribosomal protein gene expression, which in turn determines the translational capacity of cells. We validate our model by SILAC and polysome profiling experiments. As a proof of concept, we demonstrate that enhanced translation compensates for transcriptional elongation stress. Taken together, our data indicate that in addition to defining mRNA stability, components of the Ccr4-Not imprinting complex regulate RNA translatability, thus ensuring global gene expression homeostasis.