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

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.

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.

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.

Splicing-Mediated Autoregulation Modulates Rpl22p Expression in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, splicing is critical for expression of ribosomal protein genes (RPGs), which are among the most highly expressed genes and are tightly regulated according to growth and environmental conditions. However, knowledge of the precise mechanisms by which RPG pre-mRNA splicing is regulated on a gene-by-gene basis is lacking. Here we show that Rpl22p has an extraribosomal role in the inhibition of splicing of the RPL22B pre-mRNA transcript. A stem loop secondary structure within the intron is necessary for pre-mRNA binding by Rpl22p in vivo and splicing inhibition in vivo and in vitro and can rescue splicing inhibition in vitro when added in trans to splicing reactions. Splicing inhibition by Rpl22p may be partly attributed to the reduction of co-transcriptional U1 snRNP recruitment to the pre-mRNA at the RPL22B locus. We further demonstrate that the inhibition of RPL22B pre-mRNA splicing contributes to the down-regulation of mature transcript during specific stress conditions, and provide evidence hinting at a regulatory role for this mechanism in conditions of suppressed ribosome biogenesis. These results demonstrate an autoregulatory mechanism that fine-tunes the expression of the Rpl22 protein and by extension Rpl22p paralog composition according to the cellular demands for ribosome biogenesis.

Ribosomes are responsible for protein production in all living cells, serving as the grounds for the translation of genetic information from RNA to protein. Given the vital role of the ribosome in protein assembly, ribosome components are highly expressed and are subject to tight regulation. Some ribosomal proteins are also known to engage in extra-ribosomal activities. In our study, we demonstrate that the ribosomal protein Rpl22p is able to regulate its own expression by inhibiting the processing of its own RNA transcript, leading to degradation of the RNA. We also show that this self-imposed regulation plays a role in limiting RPL22 transcript levels in specific stress conditions. We suggest that this mechanism may impact the composition of ribosomes by influencing the availability of the Rpl22p paralogs.

The Rtr1p CTD phosphatase autoregulates its mRNA through a degradation pathway involving the REX exonucleases

Rtr1p is a phosphatase that impacts gene expression by modulating the phosphorylation status of the C-terminal domain of the large subunit of RNA polymerase II. Here, we show that Rtr1p is a component of a novel mRNA degradation pathway that promotes its autoregulation through turnover of its own mRNA. We show that the 3'UTR of the RTR1 mRNA contains a cis element that destabilizes this mRNA. RTR1 mRNA turnover is achieved through binding of Rtr1p to the RTR1 mRNP in a manner that is dependent on this cis element. Genetic evidence shows that Rtr1p-mediated decay of the RTR1 mRNA involves the 5'-3' DExD/H-box RNA helicase Dhh1p and the 3'-5' exonucleases Rex2p and Rex3p. Rtr1p and Rex3p are found associated with Dhh1p, suggesting a model for recruiting the REX exonucleases to the RTR1 mRNA for degradation. Rtr1p-mediated decay potentially impacts additional transcripts, including the unspliced BMH2 pre-mRNA. We propose that Rtr1p may imprint its RNA targets cotranscriptionally and determine their downstream degradation mechanism by directing these transcripts to a novel turnover pathway that involves Rtr1p, Dhh1p, and the REX family of exonucleases.

Expression homeostasis during DNA replication

Genome replication introduces a stepwise increase in the DNA template available for transcription. Genes replicated early in S phase experience this increase before late-replicating genes, raising the question of how expression levels are affected by DNA replication. We show that in budding yeast, messenger RNA (mRNA) synthesis rate is buffered against changes in gene dosage during S phase. This expression homeostasis depends on acetylation of H3 on its internal K56 site by Rtt109/Asf1. Deleting these factors, mutating H3K56 or up-regulating its deacetylation, increases gene expression in S phase in proportion to gene replication timing. Therefore, H3K56 acetylation on newly deposited histones reduces transcription efficiency from replicated DNA, complementing its role in guarding genome stability. Our study provides molecular insight into the mechanism maintaining expression homeostasis during DNA replication.

Rpb1 foot mutations demonstrate a major role of Rpb4 in mRNA stability during stress situations in yeast

The RPB1 mutants in the foot region of RNA polymerase II affect the assembly of the complex by altering the correct association of both the Rpb6 and the Rpb4/7 dimer. Assembly defects alter both transcriptional activity as well as the amount of enzyme associated with genes. Here, we show that the global transcriptional analysis of foot mutants reveals the activation of an environmental stress response (ESR), which occurs at a permissive temperature under optimal growth conditions. Our data indicate that the ESR that occurs in foot mutants depends mostly on a global post-transcriptional regulation mechanism which, in turn, depends on Rpb4-mRNA imprinting. Under optimal growth conditions, we propose that Rpb4 serves as a key to globally modulate mRNA stability as well as to coordinate transcription and decay. Overall, our results imply that post-transcriptional regulation plays a major role in controlling the ESR at both the transcription and mRNA decay levels.

Determinants of RNA metabolism in the Schizosaccharomyces pombe genome

To decrypt the regulatory code of the genome, sequence elements must be defined that determine the kinetics of RNA metabolism and thus gene expression. Here, we attempt such decryption in an eukaryotic model organism, the fission yeast S. pombe. We first derive an improved genome annotation that redefines borders of 36% of expressed mRNAs and adds 487 non‐coding RNAs (ncRNAs). We then combine RNA labeling in vivo with mathematical modeling to obtain rates of RNA synthesis and degradation for 5,484 expressed RNAs and splicing rates for 4,958 introns. We identify functional sequence elements in DNA and RNA that control RNA metabolic rates and quantify the contributions of individual nucleotides to RNA synthesis, splicing, and degradation. Our approach reveals distinct kinetics of mRNA and ncRNA metabolism, separates antisense regulation by transcription interference from RNA interference, and provides a general tool for studying the regulatory code of genomes.

Genome-wide quantification of 5'-phosphorylated mRNA degradation intermediates for analysis of ribosome dynamics

Co-translational mRNA degradation is a widespread process in which 5'-3' exonucleolytic degradation follows the last translating ribosome, thus producing an in vivo ribosomal footprint that delimits the 5' position of the mRNA molecule within the ribosome. To study this degradation process and ribosome dynamics, we developed 5PSeq, which is a method that profiles the genome-wide abundance of mRNA degradation intermediates by virtue of their 5'-phosphorylated (5'P) ends. The approach involves targeted ligation of an oligonucleotide to the 5'P end of mRNA degradation intermediates, followed by depletion of rRNA molecules, reverse transcription of 5'P mRNAs and Illumina high-throughput sequencing. 5PSeq can identify translational pauses at rare codons that are often masked when using alternative methods. This approach can be applied to previously extracted RNA samples, and it is straightforward and does not require polyribosome purification or in vitro RNA footprinting. The protocol we describe here can be applied to Saccharomyces cerevisiae and potentially to other eukaryotic organisms. Three days are required to generate 5PSeq libraries.

Ccr4–Not is at the core of the eukaryotic gene expression circuitry

In this mini-review, we summarize our current knowledge about the cross-talk between the different levels of gene expression. We introduce the Ccr4 (carbon catabolite repressed 4)–Not (negative on TATA-less) complex as a candidate to be a master regulator that orchestrates between the different levels of gene expression. An integrated view of the findings about the Ccr4–Not complex suggests that it is involved in gene expression co-ordination. Since the discovery of the Not proteins in a selection for transcription regulators in yeast [Collart and Struhl (1994) Genes Dev. 8, 525–537], the Ccr4–Not complex has been connected to every step of the mRNA lifecycle. Moreover, it has been found to be relevant for appropriate protein folding and quaternary protein structure by being involved in co-translational protein complex assembly.

Tracking Distinct RNA Populations Using Efficient and Reversible Covalent Chemistry

We describe a chemical method to label and purify 4-thiouridine (s(4)U)-containing RNA. We demonstrate that methanethiosulfonate (MTS) reagents form disulfide bonds with s(4)U more efficiently than the commonly used HPDP-biotin, leading to higher yields and less biased enrichment. This increase in efficiency allowed us to use s(4)U labeling to study global microRNA (miRNA) turnover in proliferating cultured human cells without perturbing global miRNA levels or the miRNA processing machinery. This improved chemistry will enhance methods that depend on tracking different populations of RNA, such as 4-thiouridine tagging to study tissue-specific transcription and dynamic transcriptome analysis (DTA) to study RNA turnover.

Tip60 complex binds to active Pol II promoters and a subset of enhancers and co-regulates the c-Myc network in mouse embryonic stem cells

Background

Tip60 (KAT5) is the histone acetyltransferase (HAT) of the mammalian Tip60/NuA4 complex. While Tip60 is important for early mouse development and mouse embryonic stem cell (mESC) pluripotency, the function of Tip60 as reflected in a genome-wide context is not yet well understood.

Results

Gel filtration of nuclear mESCs extracts indicate incorporation of Tip60 into large molecular complexes and exclude the existence of large quantities of “free” Tip60 within the nuclei of ESCs. Thus, monitoring of Tip60 binding to the genome should reflect the behaviour of Tip60-containing complexes. The genome-wide mapping of Tip60 binding in mESCs by chromatin immunoprecipitation (ChIP) coupled with high-throughput sequencing (ChIP-seq) shows that the Tip60 complex is present at promoter regions of predominantly active genes that are bound by RNA polymerase II (Pol II) and contain the H3K4me3 histone mark. The coactivator HAT complexes, Tip60- and Mof (KAT8)-containing (NSL and MSL), show a global overlap at promoters, whereas distinct binding profiles at enhancers suggest different regulatory functions of each essential HAT complex. Interestingly, Tip60 enrichment peaks at about 200 bp downstream of the transcription start sites suggesting a function for the Tip60 complexes in addition to histone acetylation. The comparison of genome-wide binding profiles of Tip60 and c-Myc, a somatic cell reprogramming factor that binds predominantly to active genes in mESCs, demonstrate that Tip60 and c-Myc co-bind at 50–60 % of their binding sites. We also show that the Tip60 complex binds to a subset of bivalent developmental genes and defines a set of mESC-specific enhancer as well as super-enhancer regions.

Conclusions

Our study suggests that the Tip60 complex functions as a global transcriptional co-activator at most active Pol II promoters, co-regulates the ESC-specific c-Myc network, important for ESC self-renewal and cell metabolism and acts at a subset of active distal regulatory elements, or super enhancers, in mESCs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0039-z) contains supplementary material, which is available to authorized users.

Next-generation analysis of gene expression regulation--comparing the roles of synthesis and degradation

Technological advances now enable routine measurement of mRNA and protein abundances, and estimates of their rates of synthesis and degradation that inform on their values and the degree of change in response to stimuli. Importantly, more and more data on time-series experiments are emerging, e.g. of cells responding to stress, enabling first insights into a new dimension of gene expression regulation - its dynamics and how it allows for very different response signals across genes. This review discusses recently published methods and datasets, their impact on what we now know about the relationships between concentrations and synthesis rates of mRNAs and proteins in yeast and mammalian cells, their evolution, and new hypotheses on translation regulatory mechanisms generated by approaches that involve ribosome footprinting.

Measuring Transcription Rate Changes via Time-Course 4-Thiouridine Pulse-Labelling Improves Transcriptional Target Identification

Identifying changes in the transcriptional regulation of target genes from high-throughput studies is important for unravelling molecular mechanisms controlled by a given perturbation. When measuring global transcript levels only, the effect of the perturbation [e.g., transcription factor (TF) overexpression or drug treatment] on its target genes is often obscured by delayed feedback and secondary effects until the changes are fully propagated. As a proof of principle, we show that selective measuring of transcripts that are only synthesised after a perturbation [4-thiouridine (4sU) sequencing (4sU-seq)] is a more sensitive method to identify targets and time-dependent transcriptional responses than global transcript profiling. By metabolically labelling RNA in a time-course setup, we could vastly increase the sensitivity of MYCN target gene detection compared to traditional RNA sequencing. The validity of targets identified by 4sU-seq was demonstrated using chromatin immunoprecipitation sequencing and neuroblastoma microarray tumour data. Here, we describe the methodology, both molecular biology and computational aspects, required to successfully apply this 4sU-seq approach.

Gene looping facilitates TFIIH kinase-mediated termination of transcription

TFIIH is a general transcription factor with kinase and helicase activities. The kinase activity resides in the Kin28 subunit of TFIIH. The role of Kin28 kinase in the early steps of transcription is well established. Here we report a novel role of Kin28 in the termination of transcription. We show that RNAPII reads through a termination signal upon kinase inhibition. Furthermore, the recruitment of termination factors towards the 3′ end of a gene was compromised in the kinase mutant, thus confirming the termination defect. A concomitant decrease in crosslinking of termination factors near the 5′ end of genes was also observed in the kinase-defective mutant. Simultaneous presence of termination factors towards both the ends of a gene is indicative of gene looping; while the loss of termination factor occupancy from the distal ends suggest the abolition of a looped gene conformation. Accordingly, CCC analysis revealed that the looped architecture of genes was severely compromised in the Kin28 kinase mutant. In a looping defective sua7-1 mutant, even the enzymatically active Kin28 kinase could not rescue the termination defect. These results strongly suggest a crucial role of Kin28 kinase-dependent gene looping in the termination of transcription in budding yeast.