Cells respond to endogenous and exogenous cues, through the dynamic alteration of their gene expression program. In recent years, a tremendous development of genome-wide methodologies allows the precise and comprehensive description of transcriptome changes in different conditions. In most transcriptomic studies, microarray hybridization or high-throughput sequencing are used to quantify RNA levels from a total steady-state RNA fraction. Transcriptional changes under a specific perturbation can display a wide range of possible outcomes, with either specific gene expression changes or a large spectrum of genes being either up- or downregulated. Gene expression is the result of a fine-tuned equilibrium — or steady-state — between RNA synthesis by RNA polymerases and other processes affecting RNA levels. RNA polymerase II transcription, including its three distinct phases (initiation, elongation, and termination), is highly and intricately associated with mRNA processing, cytoplasmic export, translation, and degradation.
Several recent studies demonstrated that mRNAs synthesis and decay are coupled mechanisms and showed that transcriptional effects upon mutation or under stimuli can be overlooked when quantifying total steady-state RNA. First, the detection of transcriptional changes through the analyses of steady-state levels of mRNA always depends on mRNAs half-life. Once the perturbation is introduced, the steady-state levels of mRNAs with long half-lives will be much less affected than those of mRNAs with short half-lives. Therefore, the detectability of the changes in RNA synthesis is strongly biased in favor of short-lived transcripts, while the analysis of longer-lived mRNA species might fail to reveal dynamic changes in transcription rate. Second, several reports have shown that, both in yeast and mammals, global changes in transcription might be overlooked when analyzing the steady-state levels of mRNA. This is likely due to the mechanisms that link mRNA synthesis and degradation resulting in mRNA buffering. This prompted the development of new protocols to quantify mRNA synthesis uncoupled from degradation, through the analysis of newly transcribed mRNA. In recent years, several alternatives have been presented, including global run-on sequencing (GRO-seq)1, and native elongation transcript sequencing (NET-seq)2,3. Here, we are presenting a protocol initially developed in mammalian cells4,5,6 and then adapted to yeast7,8,9,10,11, which is based on RNA labeling with a thiolated nucleoside or base analog, 4-thiouridine (4sU) or 4-thiouracil (4tU), respectively.
This method specifically purifies newly transcribed RNA from the cells in which RNA are pulse-labeled with 4sU with virtually no interference in the cell homeostasis. Hence, once the cells are exposed to 4sU, the molecule is rapidly uptaken, phosphorylated to 4sU-triphosphate, and incorporated in RNAs being transcribed. Once pulse-labeled, it is possible to extract total cellular RNA (corresponding to steady-state levels of RNA), and, subsequently, the 4sU-labeled RNA fraction is thiol-specifically modified, leading to the formation of a disulfide bond between biotin and the newly transcribed RNA4,5. However, 4sU can only be uptaken by the cells expressing a nucleoside transporter, like the human equilibrative nucleoside transporter (hENT1), preventing its immediate use in budding or fission yeast. While one could express hENT1 in either S. pombe or S. cerevisiae, an easier approach can be achieved using the modified base 4tU, since yeast cells can take up 4tU, without the need of expression of a nucleoside transporter10,11,12,13. In fact, the metabolism of 4tU requires the activity of the enzyme uracil phosphoribosyltransferase (UPRT). In several organisms, including yeast but not mammals, UPRT is essential for a pyrimidine salvage pathway, recycling uracil to uridine monophosphate.
An important bias in transcriptomic studies can be introduced by the normalization between different samples analyzed in parallel. Indeed, many deviating factors can affect the comparative analysis of the transcriptome of mutant and wild-type strains: the efficiency of cell lysis, differences in the extraction and recovery of RNA, and variances in scanner calibration for microarray analyses, among others. As discussed above, such variations can be particularly misleading when global effects on RNA polymerase II transcription are expected. An elegant mean to accurately compare mRNA synthesis rates between different samples was designed by using the distantly related fission yeast Schizosaccharomyces pombe as an internal standard. For that, a fixed number of labeled S. pombe cells is added to the S. cerevisiae samples, either wild-type or mutant cells, prior to cell lysis and RNA extraction10. Subsequently, both steady-state and newly synthesized RNAs from S. pombe and S. cerevisiae are quantified either by RT-qPCR or via the use of microarray chips or high-throughput sequencing10. Combining these data with kinetic modeling, absolute rates of mRNA synthesis and decay in budding yeast can be measured.
In the framework of this manuscript, we will show how the analysis of newly transcribed RNA allowed to reveal a global role for the coactivator complexes SAGA and TFIID in RNA polymerase II transcription in budding yeast14,15,16. Importantly, past studies quantified steady-state mRNA levels in S. cerevisiae and suggested that SAGA plays a predominant function on a limited set of yeast genes which are strongly affected by mutations in SAGA but relatively resistant to TFIID mutations17,18,19. Surprisingly, the SAGA enzymatic activities were shown to act on the whole transcribed genome, suggesting a broader role for this co-activator in RNA polymerase II transcription. Decreased RNA polymerase II recruitment at most expressed genes was observed upon the inactivation of SAGA or TFIID, suggesting that these coactivators work together at most genes. Hence, the quantification of newly transcribed mRNA revealed that SAGA and TFIID are required for the transcription of nearly all genes by RNA polymerase II14,15,16. The implementation of compensatory mechanisms emerges as a way for the cells to cope with a global decrease in mRNA synthesis which is buffered by a simultaneous global decrease in mRNA degradation. SAGA adds to the list of factors having a global effect on RNA polymerase II transcription, such as RNA Pol II subunits10, the Mediator coactivator complex20, the general transcription factor TFIIH21,22, and indirectly, elements of the mRNA degradation machinery9,10,23. Such compensatory events were universally observed in SAGA mutants, accounting for the modest and limited changes in steady-state mRNA levels despite a global and severe decrease in mRNA synthesis14. Similar analyses were also performed in a BRE1 deletion strain, resulting in a complete loss of histone H2B ubiquitination. Interestingly, a much milder but consistent global effect on RNA polymerase II transcription could be detected in the absence of Bre1, indicating that metabolic labeling of newly transcribed RNA in yeast can detect and quantify a wide range of changes in mRNA synthesis rates.