Specific and global regulation of mRNA stability during osmotic stress in Saccharomyces cerevisiae

Controlling gene expression in response to stress

Acute stress puts cells at risk, and rapid adaptation is crucial for maximizing cell survival. Cellular adaptation mechanisms include modification of certain aspects of cell physiology, such as the induction of efficient changes in the gene expression programmes by intracellular signalling networks. Recent studies using genome-wide approaches as well as single-cell transcription measurements, in combination with classical genetics, have shown that rapid and specific activation of gene expression can be accomplished by several different strategies. This article discusses how organisms can achieve generic and specific responses to different stresses by regulating gene expression at multiple stages of mRNA biogenesis from chromatin structure to transcription, mRNA stability and translation.

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

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

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

Diverse environmental stresses elicit distinct responses at the level of pre-mRNA processing in yeast

Gene expression in eukaryotic cells is profoundly influenced by the post-transcriptional processing of mRNAs, including the splicing of introns in the nucleus and both nuclear and cytoplasmic degradation pathways. These processes have the potential to affect both the steady-state levels and the kinetics of changes to levels of intron-containing transcripts. Here we report the use of a splicing isoform-specific microarray platform to investigate the effects of diverse stress conditions on pre-mRNA processing. Interestingly, we find that diverse stresses cause distinct patterns of changes at this level. The responses we observed are most dramatic for the RPGs and can be categorized into three major classes. The first is characterized by accumulation of RPG pre-mRNA and is seen in multiple types of amino acid starvation regimes; the magnitude of splicing inhibition correlates with the severity of the stress. The second class is characterized by a rapid decrease in both pre- and mature RPG mRNA and is seen in many stresses that inactivate the TORC1 kinase complex. These decreases depend on nuclear turnover of the intron-containing pre-RNAs. The third class is characterized by a decrease in RPG pre-mRNA, with only a modest reduction in the mature species; this response is observed in hyperosmotic and cation-toxic stresses. We show that casein kinase 2 (CK2) makes important contributions to the changes in pre-mRNA processing, particularly for the first two classes of stress responses. In total, our data suggest that complex post-transcriptional programs cooperate to fine-tune expression of intron-containing transcripts in budding yeast.

Coupled evolution of transcription and mRNA degradation

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

Regulatory mechanisms and networks couple the different phases of gene expression

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

Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast

To obtain rates of mRNA synthesis and decay in yeast, we established dynamic transcriptome analysis (DTA). DTA combines non-perturbing metabolic RNA labeling with dynamic kinetic modeling. DTA reveals that most mRNA synthesis rates are around several transcripts per cell and cell cycle, and most mRNA half-lives range around a median of 11 min. DTA can monitor the cellular response to osmotic stress with higher sensitivity and temporal resolution than standard transcriptomics. In contrast to monotonically increasing total mRNA levels, DTA reveals three phases of the stress response. During the initial shock phase, mRNA synthesis and decay rates decrease globally, resulting in mRNA storage. During the subsequent induction phase, both rates increase for a subset of genes, resulting in production and rapid removal of stress-responsive mRNAs. During the recovery phase, decay rates are largely restored, whereas synthesis rates remain altered, apparently enabling growth at high salt concentration. Stress-induced changes in mRNA synthesis rates are predicted from gene occupancy with RNA polymerase II. DTA-derived mRNA synthesis rates identified 16 stress-specific pairs/triples of cooperative transcription factors, of which seven were known. Thus, DTA realistically monitors the dynamics in mRNA metabolism that underlie gene regulatory systems.

Heat shock response in yeast involves changes in both transcription rates and mRNA stabilities

We have analyzed the heat stress response in the yeast Saccharomyces cerevisiae by determining mRNA levels and transcription rates for the whole transcriptome after a shift from 25 °C to 37 °C. Using an established mathematical algorithm, theoretical mRNA decay rates have also been calculated from the experimental data. We have verified the mathematical predictions for selected genes by determining their mRNA decay rates at different times during heat stress response using the regulatable tetO promoter. This study indicates that the yeast response to heat shock is not only due to changes in transcription rates, but also to changes in the mRNA stabilities. mRNA stability is affected in 62% of the yeast genes and it is particularly important in shaping the mRNA profile of the genes belonging to the environmental stress response. In most cases, changes in transcription rates and mRNA stabilities are homodirectional for both parameters, although some interesting cases of antagonist behavior are found. The statistical analysis of gene targets and sequence motifs within the clusters of genes with similar behaviors shows that both transcriptional and post-transcriptional regulons apparently contribute to the general heat stress response by means of transcriptional factors and RNA binding proteins.

Genomic-wide methods to evaluate transcription rates in yeast

Gene transcription is a dynamic process in which the desired amount of an mRNA is obtained by the equilibrium between its transcription (TR) and degradation (DR) rates. The control mechanism at the RNA polymerase level primarily causes changes in TR. Despite their importance, TRs have been rarely measured. In the yeast Saccharomyces cerevisiae, we have implemented two techniques to evaluate TRs: run-on and chromatin immunoprecipitation of RNA polymerase II. These techniques allow the discrimination of the relative importance of TR and DR in gene regulation for the first time in a eukaryote.

Global estimation of mRNA stability in yeast

Turnover of mRNA is an important level of gene regulation. Individual mRNAs have different intrinsic stabilities. Moreover, mRNA stability changes dynamically with conditions such as hormonal stimulation or cellular stress. While accurate methods exist to measure the half-life of an individual transcript, global methods to estimate mRNA turnover have limitations in terms of resolution in time and precision. We describe and compare two complementary approaches to estimating global transcript stability: (1) direct measurement of decay rates; (2) indirect estimation of turnover from determination of mRNA synthesis rates and steady-state levels. Since the two approaches have distinct strengths yet confer different cellular perturbations, it is valuable to consider results obtained with both methods. The practical aspects of the chapter are written from a yeast perspective; the general considerations hold true for all eukaryotes, however.

Cellular stress induces cytoplasmic RNA granules in fission yeast

Severe stress causes plant and animal cells to form large cytoplasmic granules containing RNA and proteins. Here, we demonstrate the existence of stress-induced cytoplasmic RNA granules in Schizosaccharomyces pombe. Homologs to several known protein components of mammalian processing bodies and stress granules are found in fission yeast RNA granules. In contrast to mammalian cells, poly(A)-binding protein (Pabp) colocalizes in stress-induced granules with decapping protein. After glucose deprivation, protein kinase A (PKA) is required for accumulation of Pabp-positive granules and translational down-regulation. This is the first demonstration of a role for PKA in RNA granule formation. In mammals, the translation initiation protein eIF2α is a key regulator of formation of granules containing poly(A)-binding protein. In S. pombe, nonphosphorylatable eIF2α does not block but delays granule formation and subsequent clearance after exposure to hyperosmosis. At least two separate pathways in S. pombe appear to regulate stress-induced granules: pka1 mutants are fully proficient to form granules after hyperosmotic shock; conversely, eIF2α does not affect granule formation in glucose starvation. Further, we demonstrate a Pka1-dependent link between calcium perturbation and RNA granules, which has not been described earlier in any organism.

A complete set of nascent transcription rates for yeast genes

The amount of mRNA in a cell is the result of two opposite reactions: transcription and mRNA degradation. These reactions are governed by kinetics laws, and the most regulated step for many genes is the transcription rate. The transcription rate, which is assumed to be exercised mainly at the RNA polymerase recruitment level, can be calculated using the RNA polymerase densities determined either by run-on or immunoprecipitation using specific antibodies. The yeast Saccharomyces cerevisiae is the ideal model organism to generate a complete set of nascent transcription rates that will prove useful for many gene regulation studies. By combining genomic data from both the GRO (Genomic Run-on) and the RNA pol ChIP-on-chip methods we generated a new, more accurate nascent transcription rate dataset. By comparing this dataset with the indirect ones obtained from the mRNA stabilities and mRNA amount datasets, we are able to obtain biological information about posttranscriptional regulation processes and a genomic snapshot of the location of the active transcriptional machinery. We have obtained nascent transcription rates for 4,670 yeast genes. The median RNA polymerase II density in the genes is 0.078 molecules/kb, which corresponds to an average of 0.096 molecules/gene. Most genes have transcription rates of between 2 and 30 mRNAs/hour and less than 1% of yeast genes have >1 RNA polymerase molecule/gene. Histone and ribosomal protein genes are the highest transcribed groups of genes and other than these exceptions the transcription of genes is an infrequent phenomenon in a yeast cell.

Turnover of AU-rich-containing mRNAs during stress: a matter of survival

Cells undergo various adaptive measures in response to stress. Among these are specific changes in the posttranscriptional regulation of various genes. In particular, the turnover of mRNA is modified to either increase or decrease the abundance of certain target messages. Some of the best-studied mRNAs that are affected by stress are those that contain adenine/uridine-rich elements (AREs) in their 3'-untranslated regions. ARE-containing mRNAs are involved in many important cellular processes and are normally labile, but in response to stress they are differentially regulated through the concerted efforts of ARE-binding proteins (AUBPs) such as HuR, AUF1, tristetraprolin, BRF1, and KSRP, along with microRNA-mediated effects. Additionally, the fate of ARE-containing mRNAs is modified by inducing their localization to stress granules or mRNA processing bodies. Coordination of these various mechanisms controls the turnover of ARE-containing mRNAs, and thereby enables proper responses to cellular stress. In this review, we discuss how AUBPs regulate their target mRNAs in response to stress, along with the involvement of cytoplasmic granules in this process.

Toward a genomic view of the gene expression program regulated by osmostress in yeast

Osmostress triggers profound adaptive changes in the physiology of the cell with a great impact on gene expression. Saccharomyces cerevisiae has served as an instructive model system to unravel the complexity of the stress response at the transcriptional level. The main signal transduction pathways like the HOG (high osmolarity glycerol) MAP kinase cascade or the protein kinase A pathway regulate multiple specific transcription factors to accomplish large changes in the expression pattern of the genome. Transcription profiling and proteomic studies give us an idea about the impact of osmostress on gene expression and the overall protein composition. Recent genome wide location studies for several transcription factors and signaling kinases involved in the transcriptional stress response shed light on the genomic organization of the osmostress response at the level of the dynamic association of regulators with chromatin. Finally, global surveys of mRNA stability complete our picture of the mechanisms underlying the massive reprogramming of global gene expression, which leads to efficient adaptation to osmotic stress.

Defects in the secretory pathway and high Ca2+ induce multiple P-bodies

mRNA is sequestered and turned over in cytoplasmic processing bodies (PBs), which are induced by various cellular stresses. Unexpectedly, in Saccharomyces cerevisiae, mutants of the small GTPase Arf1 and various secretory pathway mutants induced a significant increase in PB number, compared with PB induction by starvation or oxidative stress. Exposure of wild-type cells to osmotic stress or high extracellular Ca(2+) mimicked this increase in PB number. Conversely, intracellular Ca(2+)-depletion strongly reduced PB formation in the secretory mutants. In contrast to PB induction through starvation or osmotic stress, PB formation in secretory mutants and by Ca(2+) required the PB components Pat1 and Scd6, and calmodulin, indicating that different stressors act through distinct pathways. Consistent with this hypothesis, when stresses were combined, PB number did not correlate with the strength of the translational block, but rather with the type of stress encountered. Interestingly, independent of the stressor, PBs appear as spheres of approximately 40-100 nm connected to the endoplasmic reticulum (ER), consistent with the idea that translation and silencing/degradation occur in a spatially coordinated manner at the ER. We propose that PB assembly in response to stress occurs at the ER and depends on intracellular signals that regulate PB number.

The HOG pathway dictates the short-term translational response after hyperosmotic shock

Cellular responses to environmental changes occur on different levels. We investigated the translational response of yeast cells after mild hyperosmotic shock by isolating mRNA associated with multiple ribosomes (polysomes) followed by array analysis. Globally, recruitment of preexisting mRNAs to ribosomes (translational response) is faster than the transcriptional response. Specific functional groups of mRNAs are recruited to ribosomes without any corresponding increase in total mRNA. Among mRNAs under strong translational up-regulation upon shock, transcripts encoding membrane-bound proteins including hexose transporters were enriched. Similarly, numerous mRNAs encoding cytoplasmic ribosomal proteins run counter to the overall trend of down-regulation and are instead translationally mobilized late in the response. Surprisingly, certain transcriptionally induced mRNAs were excluded from ribosomal association after shock. Importantly, we verify, using constructs with intact 5' and 3' untranslated regions, that the observed changes in polysomal mRNA are reflected in protein levels, including cases with only translational up-regulation. Interestingly, the translational regulation of the most highly osmostress-regulated mRNAs was more strongly dependent on the stress-activated protein kinases Hog1 and Rck2 than the transcriptional regulation. Our results show the importance of translational control for fine tuning of the adaptive responses.

Defects in the secretory pathway and high Ca2+ induce multiple P-bodies

mRNA is sequestered and turned over in cytoplasmic processing bodies (PBs), which are induced by various cellular stresses. Unexpectedly, in Saccharomyces cerevisiae, mutants of the small GTPase Arf1 and various secretory pathway mutants induced a significant increase in PB number, compared with PB induction by starvation or oxidative stress. Exposure of wild-type cells to osmotic stress or high extracellular Ca(2+) mimicked this increase in PB number. Conversely, intracellular Ca(2+)-depletion strongly reduced PB formation in the secretory mutants. In contrast to PB induction through starvation or osmotic stress, PB formation in secretory mutants and by Ca(2+) required the PB components Pat1 and Scd6, and calmodulin, indicating that different stressors act through distinct pathways. Consistent with this hypothesis, when stresses were combined, PB number did not correlate with the strength of the translational block, but rather with the type of stress encountered. Interestingly, independent of the stressor, PBs appear as spheres of approximately 40-100 nm connected to the endoplasmic reticulum (ER), consistent with the idea that translation and silencing/degradation occur in a spatially coordinated manner at the ER. We propose that PB assembly in response to stress occurs at the ER and depends on intracellular signals that regulate PB number.

Major role for mRNA stability in shaping the kinetics of gene induction

Background

mRNA levels in cells are determined by the relative rates of RNA production and degradation. Yet, to date, most analyses of gene expression profiles were focused on mechanisms which regulate transcription, while the role of mRNA stability in modulating transcriptional networks was to a large extent overlooked. In particular, kinetic waves in transcriptional responses are usually interpreted as resulting from sequential activation of transcription factors.

Results

In this study, we examined on a global scale the role of mRNA stability in shaping the kinetics of gene response. Analyzing numerous expression datasets we revealed a striking global anti-correlation between rapidity of induction and mRNA stability, fitting the prediction of a kinetic mathematical model. In contrast, the relationship between kinetics and stability was less significant when gene suppression was analyzed. Frequently, mRNAs that are stable under standard conditions were very rapidly down-regulated following stimulation. Such effect cannot be explained even by a complete shut-off of transcription, and therefore indicates intense modulation of RNA stability.

Conclusion

Taken together, our results demonstrate the key role of mRNA stability in determining induction kinetics in mammalian transcriptional networks.

The bidirectional cytomegalovirus immediate/early promoter is regulated by Hog1 and the stress transcription factors Sko1 and Hot1 in yeast

The work presented here intends to address the question of whether the immediate/early promoter of cytomegalovirus (CMV), which is widely used for expressing transgenes in eukaryotic cells, yields a constitutive expression of the transgenes under stress conditions in Saccharomyces cerevisiae cells. This information would also be relevant because in the tetracycline-regulated expression (tetO) system, which is one of the first choices for studying gene function from yeast to human cells, the CMV promoter controls the expression of the tetO transactivator. We found that the CMV promoter in yeast cells is bidirectionally induced by osmotic stress and in glycerol media. The mitogen-activated protein (MAP) kinase Hog1 controls CMV activation by osmotic stress through the ATF/CRE-related transcription factor Sko1 and the yeast osmostress factor Hot1. Our results indicate that the CMV and tetO expression systems respond to external signals and this should be considered before using these systems in yeast. Moreover, our results also suggest that CMV could be regulated by the intracellular glucose concentration in human cells.

Comparing transcription rate and mRNA abundance as parameters for biochemical pathway and network analysis

The cells adapt to extra- and intra-cellular signals by dynamic orchestration of activities of pathways in the biochemical networks. Dynamic control of the gene expression process represents a major mechanism for pathway activity regulation. Gene expression has thus been routinely measured, most frequently at steady-state mRNA abundance level using micro-array technology. The results are widely used in statistical inference of the structures of underlying biochemical networks, with the assumption that functionally related genes exhibit similar dynamic profiles. Steady-state mRNA abundance, however, is a composite of two factors: transcription rate and mRNA degradation rate. The question being asked here is therefore whether steady-state mRNA abundance or any of two factors is a more informative measurement target for studying network dynamics. The yeast S. cerevisiae was used as model organism and transcription rate was chosen out of the two factors in this study, because genome-wide determination of transcription rates has been reported for several physiological processes in this species. Our strategy is to test which one is a better measurement of functional relatedness between genes. The analysis was performed on those S. cerevisiae genes that have bacterial orthologs as identified by reciprocal BLAST analysis, so that functional relatedness of a gene pair can be measured by the frequency at which their bacterial orthologs co-occur in the same operon in the collection of bacterial genomes. It is found that transcription rate data is generally a better parameter for functional relatedness than steady state mRNA abundance, suggesting transcription rate data is more informative to use in deciphering the logics used by the cells in dynamic regulation of biochemical network behaviors. The significance of this finding for network and systems biology, as well as biomedical research in general, is discussed.

Kaposi's sarcoma-associated herpesvirus ORF57 protein binds and protects a nuclear noncoding RNA from cellular RNA decay pathways

The control of RNA stability is a key determinant in cellular gene expression. The stability of any transcript is modulated through the activity of cis- or trans-acting regulatory factors as well as cellular quality control systems that ensure the integrity of a transcript. As a result, invading viral pathogens must be able to subvert cellular RNA decay pathways capable of destroying viral transcripts. Here we report that the Kaposi's sarcoma-associated herpesvirus (KSHV) ORF57 protein binds to a unique KSHV polyadenylated nuclear RNA, called PAN RNA, and protects it from degradation by cellular factors. ORF57 increases PAN RNA levels and its effects are greatest on unstable alleles of PAN RNA. Kinetic analysis of transcription pulse assays shows that ORF57 protects PAN RNA from a rapid cellular RNA decay process, but ORF57 has little effect on transcription or PAN RNA localization based on chromatin immunoprecipitation and in situ hybridization experiments, respectively. Using a UV cross-linking technique, we further demonstrate that ORF57 binds PAN RNA directly in living cells and we show that binding correlates with function. In addition, we define an ORF57-responsive element (ORE) that is necessary for ORF57 binding to PAN RNA and sufficient to confer ORF57-response to a heterologous intronless beta-globin mRNA, but not its spliced counterparts. We conclude that ORF57 binds to viral transcripts in the nucleus and protects them from a cellular RNA decay pathway. We propose that KSHV ORF57 protein functions to enhance the nuclear stability of intronless viral transcripts by protecting them from a cellular RNA quality control pathway.

Multilayered control of gene expression by stress-activated protein kinases

Stress-activated protein kinases (SAPKs) are key elements for intracellular signalling networks that serve to respond and adapt to extracellular changes. Exposure of yeast to high osmolarity results in the activation of p38-related SAPK, Hog1, which is essential for reprogramming the gene expression capacity of the cell by regulation of several steps of the transcription process. At initiation, active Hog1 not only directly phosphorylates several transcription factors to alter their activities, but also associates at stress-responsive promoters through such transcription factors. Once at the promoters, Hog1 serves as a platform to recruit general transcription factors, chromatin-modifying activities and RNA Pol II. In addition, the SAPK pathway has a role in elongation. At the stress-responsive ORFs, Hog1 recruits the RSC chromatin-remodelling complex to modify nucleosome organization. Several SAPKs from yeast to mammals have maintained some of the regulatory abilities of Hog1. Thus, elucidating the control of gene expression by the Hog1 SAPK should help to understand how eukaryotic cells implement a massive and rapid change on their transcriptional capacity in response to adverse conditions.