Dynamic profiling of mRNA turnover reveals gene-specific and system-wide regulation of mRNA decay

The NEDD4-binding protein N4BP1 degrades mRNA substrates through the coding sequence independent of nonsense-mediated decay

3'UTRs are recognized for their role in regulating mRNA turnover while the turnover of a specific group of mRNAs mediated by coding sequences (CDSs) remains poorly understood. N4BP1 is a critical inflammatory regulator in vivo with a molecular mechanism that is not yet clearly defined. Our study reveals that N4BP1 efficiently degrades its mRNA targets via CDS rather than the 3'-UTR. This CDS-dependent mRNA turnover mechanism appears to be a general feature of N4BP1, as evidenced by testing multiple mRNA substrates, such as Fos-C, Fos-B, Jun-B, and C-X-C motif chemokine ligand 1. Detailed mapping of the motif identified a crucial 33-nt (289-322) sequence near the 5'-end of Fos-C-CDS, where the presence of polyC is necessary for N4BP1-mediated degradation. Functional studies involving domain deletion and point mutations showed that both the K homology and N4BP1, YacP-like nuclease domains are essential for N4BP1 to restrict mRNA substrates. The function of N4BP1 in mRNA turnover is not dependent on nonsense-mediated decay as it efficiently restricts mRNA substrates even in cells deficient in UPF1, UPF3A, and UPF3B. Additionally, the function of N4BP1 is not reliant on LUC7L3 despite its known association with this protein. Our findings suggest that N4BP1 acts as an endoribonuclease to degrade mRNA substrates primarily through CDSs containing a C-rich motif.

Proteome Dynamics in iPSC-Derived Human Dopaminergic Neurons

Dopaminergic neurons participate in fundamental physiological processes and are the cell type primarily affected in Parkinson's disease. Their analysis is challenging due to the intricate nature of their function, involvement in diverse neurological processes, and heterogeneity and localization in deep brain regions. Consequently, most of the research on the protein dynamics of dopaminergic neurons has been performed in animal cells ex vivo. Here we use iPSC-derived human mid-brain-specific dopaminergic neurons to study general features of their proteome biology and provide datasets for protein turnover and dynamics, including a human axonal translatome. We cover the proteome to a depth of 9409 proteins and use dynamic SILAC to measure the half-life of more than 4300 proteins. We report uniform turnover rates of conserved cytosolic protein complexes such as the proteasome and map the variable rates of turnover of the respiratory chain complexes in these cells. We use differential dynamic SILAC labeling in combination with microfluidic devices to analyze local protein synthesis and transport between axons and soma. We report 105 potentially novel axonal markers and detect translocation of 269 proteins between axons and the soma in the time frame of our analysis (120 h). Importantly, we provide evidence for local synthesis of 154 proteins in the axon and their retrograde transport to the soma, among them several proteins involved in RNA editing such as ADAR1 and the RNA helicase DHX30, involved in the assembly of mitochondrial ribosomes. Our study provides a workflow and resource for the future applications of quantitative proteomics in iPSC-derived human neurons.

Mitochondrial protein heterogeneity stems from the stochastic nature of co-translational protein targeting in cell senescence

A decline in mitochondrial function is a hallmark of aging and neurodegenerative diseases. It has been proposed that changes in mitochondrial morphology, including fragmentation of the tubular mitochondrial network, can lead to mitochondrial dysfunction, yet the mechanism of this loss of function is unclear. Most proteins contained within mitochondria are nuclear-encoded and must be properly targeted to the mitochondria. Here, we report that sustained mRNA localization and co-translational protein delivery leads to a heterogeneous protein distribution across fragmented mitochondria. We find that age-induced mitochondrial fragmentation drives a substantial increase in protein expression noise across fragments. Using a translational kinetic and molecular diffusion model, we find that protein expression noise is explained by the nature of stochastic compartmentalization and that co-translational protein delivery is the main contributor to increased heterogeneity. We observed that cells primarily reduce the variability in protein distribution by utilizing mitochondrial fission-fusion processes rather than relying on the mitophagy pathway. Furthermore, we are able to reduce the heterogeneity of the protein distribution by inhibiting co-translational protein targeting. This research lays the framework for a better understanding of the detrimental impact of mitochondrial fragmentation on the physiology of cells in aging and disease.

Strategies to reduce the risks of mRNA drug and vaccine toxicity

mRNA formulated with lipid nanoparticles is a transformative technology that has enabled the rapid development and administration of billions of coronavirus disease 2019 (COVID-19) vaccine doses worldwide. However, avoiding unacceptable toxicity with mRNA drugs and vaccines presents challenges. Lipid nanoparticle structural components, production methods, route of administration and proteins produced from complexed mRNAs all present toxicity concerns. Here, we discuss these concerns, specifically how cell tropism and tissue distribution of mRNA and lipid nanoparticles can lead to toxicity, and their possible reactogenicity. We focus on adverse events from mRNA applications for protein replacement and gene editing therapies as well as vaccines, tracing common biochemical and cellular pathways. The potential and limitations of existing models and tools used to screen for on-target efficacy and de-risk off-target toxicity, including in vivo and next-generation in vitro models, are also discussed.

Mutations in yeast Pcf11, a conserved protein essential for mRNA 3' end processing and transcription termination, elicit the Environmental Stress Response

The Pcf11 protein is an essential subunit of the large complex that cleaves and polyadenylates eukaryotic mRNA precursor. It has also been functionally linked to gene-looping, termination of RNA Polymerase II (Pol II) transcripts, and mRNA export. We have examined a poorly characterized but conserved domain (amino acids 142-225) of the Saccharomyces cerevisiae  Pcf11 and found that while it is not needed for mRNA 3' end processing or termination downstream of the poly(A) sites of protein-coding genes, its presence improves the interaction with Pol II and the use of transcription terminators near gene promoters. Analysis of genome-wide Pol II occupancy in cells with Pcf11 missing this region, as well as Pcf11 mutated in the Pol II CTD Interacting Domain, indicates that systematic changes in mRNA expression are mediated primarily at the level of transcription. Global expression analysis also shows that a general stress response, involving both activation and suppression of specific gene sets known to be regulated in response to a wide variety of stresses, is induced in the two pcf11 mutants, even though cells are grown in optimal conditions. The mutants also cause an unbalanced expression of cell wall-related genes that does not activate the Cell Wall Integrity pathway but is associated with strong caffeine sensitivity. Based on these findings, we propose that Pcf11 can modulate the expression level of specific functional groups of genes in ways that do not involve its well-characterized role in mRNA 3' end processing.

Dynamics of mRNA fate during light stress and recovery: from transcription to stability and translation

Transcript stability is an important determinant of its abundance and, consequently, translational output. Transcript destabilisation can be rapid and is well suited for modulating the cellular response. However, it is unclear the extent to which RNA stability is altered under changing environmental conditions in plants. We previously hypothesised that recovery-induced transcript destabilisation facilitated a phenomenon of rapid recovery gene downregulation (RRGD) in Arabidopsis thaliana (Arabidopsis) following light stress, based on mathematical calculations to account for ongoing transcription. Here, we test this hypothesis and investigate processes regulating transcript abundance and fate by quantifying changes in transcription, stability and translation before, during and after light stress. We adapt syringe infiltration to apply a transcriptional inhibitor to soil-grown plants in combination with stress treatments. Compared with measurements in juvenile plants and cell culture, we find reduced stability across a range of transcripts encoding proteins involved in RNA binding and processing. We also observe light-induced destabilisation of transcripts, followed by their stabilisation during recovery. We propose that this destabilisation facilitates RRGD, possibly in combination with transcriptional shut-off that was confirmed for HSP101, ROF1 and GOLS1. We also show that translation remains highly dynamic over the course of light stress and recovery, with a bias towards transcript-specific increases in ribosome association, independent of changes in total transcript abundance, after 30 min of light stress. Taken together, we provide evidence for the combinatorial regulation of transcription and stability that occurs to coordinate translation during light stress and recovery in Arabidopsis.

The circular logic of mRNA homeostasis

Eukaryotic cells rely upon dynamic, multifaceted regulation at each step of RNA biogenesis to maintain mRNA pools and ensure normal protein synthesis. Studies in budding yeast indicate a buffering phenomenon that preserves global mRNA levels through the reciprocal balancing of RNA synthesis rates and mRNA decay. In short, changes in transcription impact the efficiency of mRNA degradation and defects in either nuclear or cytoplasmic mRNA degradation are somehow sensed and relayed to control a compensatory change in mRNA transcription rates. Here, we review current views on molecular mechanisms that might explain this apparent bidirectional sensing process that ensures homeostasis of the stable mRNA pool.

Simultaneous estimation of gene regulatory network structure and RNA kinetics from single cell gene expression

Cells respond to environmental and developmental stimuli by remodeling their transcriptomes through regulation of both mRNA transcription and mRNA decay. A central goal of biology is identifying the global set of regulatory relationships between factors that control mRNA production and degradation and their target transcripts and construct a predictive model of gene expression. Regulatory relationships are typically identified using transcriptome measurements and causal inference algorithms. RNA kinetic parameters are determined experimentally by employing run-on or metabolic labeling (e.g. 4-thiouracil) methods that allow transcription and decay rates to be separately measured. Here, we develop a deep learning model, trained with single-cell RNA-seq data, that both infers causal regulatory relationships and estimates RNA kinetic parameters. The resulting in silico model predicts future gene expression states and can be perturbed to simulate the effect of transcription factor changes. We acquired model training data by sequencing the transcriptomes of 175,000 individual Saccharomyces cerevisiae cells that were subject to an external perturbation and continuously sampled over a one hour period. The rate of change for each transcript was calculated on a per-cell basis to estimate RNA velocity. We then trained a deep learning model with transcriptome and RNA velocity data to calculate time-dependent estimates of mRNA production and decay rates. By separating RNA velocity into transcription and decay rates, we show that rapamycin treatment causes existing ribosomal protein transcripts to be rapidly destabilized, while production of new transcripts gradually slows over the course of an hour. The neural network framework we present is designed to explicitly model causal regulatory relationships between transcription factors and their genes, and shows superior performance to existing models on the basis of recovery of known regulatory relationships. We validated the predictive power of the model by perturbing transcription factors in silico and comparing transcriptome-wide effects with experimental data. Our study represents the first step in constructing a complete, predictive, biophysical model of gene expression regulation.

Evaluating single-cell variability in proteasomal decay

Gene expression is a stochastic process that leads to variability in mRNA and protein abundances even within an isogenic population of cells grown in the same environment. This variation, often called gene-expression noise, has typically been attributed to transcriptional and translational processes while ignoring the contributions of protein decay variability across cells. Here we estimate the single-cell protein decay rates of two degron GFPs in Saccharomyces cerevisiae using time-lapse microscopy. We find substantial cell-to-cell variability in the decay rates of the degron GFPs. We evaluate cellular features that explain the variability in the proteasomal decay and find that the amount of 20s catalytic beta subunit of the proteasome marginally explains the observed variability in the degron GFP half-lives. We propose alternate hypotheses that might explain the observed variability in the decay of the two degron GFPs. Overall, our study highlights the importance of studying the kinetics of the decay process at single-cell resolution and that decay rates vary at the single-cell level, and that the decay process is stochastic. A complex model of decay dynamics must be included when modeling stochastic gene expression to estimate gene expression noise.

Evolutionary trade-off and mutational bias could favor transcriptional over translational divergence within paralog pairs

How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.

Clinical Delivery of Circular RNA: Lessons Learned from RNA Drug Development

Circular RNAs (circRNA) represent a distinct class of covalently closed-loop RNA molecules, which play diverse roles in regulating biological processes and disease states. The enhanced stability of synthetic circRNAs compared to their linear counterparts has recently garnered considerable research interest, paving the way for new therapeutic applications. While clinical circRNA technology is still in its early stages, significant advancements in mRNA technology offer valuable insights into its potential future applications. Two primary obstacles that must be addressed are the development of efficient production methods and the optimization of delivery systems. To expedite progress in this area, this review aims to provide an overview of the current state of knowledge on circRNA structure and function, outline recent techniques for synthesizing circRNAs, highlight key delivery strategies and applications, and discuss the current challenges and future prospects in the field of circRNA-based therapeutics.

mRNA 5' terminal sequences drive 200-fold differences in expression through effects on synthesis, translation and decay

mRNA regulatory sequences control gene expression at multiple levels including translation initiation and mRNA decay. The 5' terminal sequences of mRNAs have unique regulatory potential because of their proximity to key post-transcriptional regulators. Here we have systematically probed the function of 5' terminal sequences in gene expression in human cells. Using a library of reporter mRNAs initiating with all possible 7-mer sequences at their 5' ends, we find an unexpected impact on transcription that underlies 200-fold differences in mRNA expression. Library sequences that promote high levels of transcription mirrored those found in native mRNAs and define two basic classes with similarities to classic Initiator (Inr) and TCT core promoter motifs. By comparing transcription, translation and decay rates, we identify sequences that are optimized for both efficient transcription and growth-regulated translation and stability, including variants of terminal oligopyrimidine (TOP) motifs. We further show that 5' sequences of endogenous mRNAs are enriched for multi-functional TCT/TOP hybrid sequences. Together, our results reveal how 5' sequences define two general classes of mRNAs with distinct growth-responsive profiles of expression across synthesis, translation and decay.

Mitotic checkpoint gene expression is tuned by codon usage bias

The mitotic checkpoint (also called spindle assembly checkpoint, SAC) is a signaling pathway that safeguards proper chromosome segregation. Correct functioning of the SAC depends on adequate protein concentrations and appropriate stoichiometries between SAC proteins. Yet very little is known about the regulation of SAC gene expression. Here, we show in the fission yeast Schizosaccharomyces pombe that a combination of short mRNA half-lives and long protein half-lives supports stable SAC protein levels. For the SAC genes mad2+ and mad3+ , their short mRNA half-lives are caused, in part, by a high frequency of nonoptimal codons. In contrast, mad1+ mRNA has a short half-life despite a higher frequency of optimal codons, and despite the lack of known RNA-destabilizing motifs. Hence, different SAC genes employ different strategies of expression. We further show that Mad1 homodimers form co-translationally, which may necessitate a certain codon usage pattern. Taken together, we propose that the codon usage of SAC genes is fine-tuned to ensure proper SAC function. Our work shines light on gene expression features that promote spindle assembly checkpoint function and suggests that synonymous mutations may weaken the checkpoint.

Nonsense-mediated RNA decay: an emerging modulator of malignancy

Nonsense-mediated RNA decay (NMD) is a highly conserved RNA turnover pathway that selectively degrades RNAs harbouring truncating mutations that prematurely terminate translation, including nonsense, frameshift and some splice-site mutations. Recent studies show that NMD shapes the mutational landscape of tumours by selecting for mutations that tend to downregulate the expression of tumour suppressor genes but not oncogenes. This suggests that NMD can benefit tumours, a notion further supported by the finding that mRNAs encoding immunogenic neoantigen peptides are typically targeted for decay by NMD. Together, this raises the possibility that NMD-inhibitory therapy could be of therapeutic benefit against many tumour types, including those with a high load of neoantigen-generating mutations. Complicating this scenario is the evidence that NMD can also be detrimental for many tumour types, and consequently tumours often have perturbed NMD. NMD may suppress tumour generation and progression by degrading subsets of specific normal mRNAs, including those encoding stress-response proteins, signalling factors and other proteins beneficial for tumours, as well as pro-tumour non-coding RNAs. Together, these findings suggest that NMD-modulatory therapy has the potential to provide widespread therapeutic benefit against diverse tumour types. However, whether NMD should be stimulated or repressed requires careful analysis of the tumour to be treated.

Regulation of both transcription and RNA turnover contribute to germline specification

The nuanced mechanisms driving primordial germ cells (PGC) specification remain incompletely understood since genome-wide transcriptional regulation in developing PGCs has previously only been defined indirectly. Here, using SLAMseq analysis, we determined genome-wide transcription rates during the differentiation of embryonic stem cells (ESCs) to form epiblast-like (EpiLC) cells and ultimately PGC-like cells (PGCLCs). This revealed thousands of genes undergoing bursts of transcriptional induction and rapid shut-off not detectable by RNAseq analysis. Our SLAMseq datasets also allowed us to infer RNA turnover rates, which revealed thousands of mRNAs stabilized and destabilized during PGCLC specification. mRNAs tend to be unstable in ESCs and then are progressively stabilized as they differentiate. For some classes of genes, mRNA turnover regulation collaborates with transcriptional regulation, but these processes oppose each other in a surprisingly high frequency of genes. To test whether regulated mRNA turnover has a physiological role in PGC development, we examined three genes that we found were regulated by RNA turnover: Sox2, Klf2 and Ccne1. Circumvention of their regulated RNA turnover severely impaired the ESC-to-EpiLC and EpiLC-to-PGCLC transitions. Our study demonstrates the functional importance of regulated RNA stability in germline development and provides a roadmap of transcriptional and post-transcriptional regulation during germline specification.

Global SLAM-seq for accurate mRNA decay determination and identification of NMD targets

Gene expression analysis requires accurate measurements of global RNA degradation rates, earlier problematic with methods disruptive to cell physiology. Recently, metabolic RNA labeling emerged as an efficient and minimally invasive technique applied in mammalian cells. Here, we have adapted SH-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) for a global mRNA stability study in yeast using 4-thiouracil pulse-chase labeling. We assign high-confidence half-life estimates for 67.5% of expressed ORFs, and measure a median half-life of 9.4 min. For mRNAs where half-life estimates exist in the literature, their ranking order was in good agreement with previous data, indicating that SLAM-seq efficiently classifies stable and unstable transcripts. We then leveraged our yeast protocol to identify targets of the nonsense-mediated decay (NMD) pathway by measuring the change in RNA half-lives, instead of steady-state RNA level changes. With SLAM-seq, we assign 580 transcripts as putative NMD targets, based on their measured half-lives in wild-type and upf3Δ mutants. We find 225 novel targets, and observe a strong agreement with previous reports of NMD targets, 61.2% of our candidates being identified in previous studies. This indicates that SLAM-seq is a simpler and more economic method for global quantification of mRNA half-lives. Our adaptation for yeast yielded global quantitative measures of the NMD effect on transcript half-lives, high correlation with RNA half-lives measured previously with more technically challenging protocols, and identification of novel NMD regulated transcripts that escaped prior detection.

StaRTrEK:in silico estimation of RNA half-lives from genome-wide time-course experiments without transcriptional inhibition

Background

Gene expression is the result of the balance between transcription and degradation. Recent experimental findings have shown fine and specific regulation of RNA degradation and the presence of various molecular machinery purposely devoted to this task, such as RNA binding proteins, non-coding RNAs, etc. A biological process can be studied by measuring time-courses of RNA abundance in response of internal and/or external stimuli, using recent technologies, such as the microarrays or the Next Generation Sequencing devices. Unfortunately, the picture provided by looking only at the transcriptome abundance may not gain insight into its dynamic regulation. By contrast, independent simultaneous measurement of RNA expression and half-lives could provide such valuable additional insight. A computational approach to the estimation of RNAs half-lives from RNA expression time profiles data, can be a low-cost alternative to its experimental measurement which may be also affected by various artifacts.

Results

Here we present a computational methodology, called StaRTrEK (STAbility Rates ThRough Expression Kinetics), able to estimate half-life values basing only on genome-wide gene expression time series without transcriptional inhibition. The StaRTrEK algorithm makes use of a simple first order kinetic model and of a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$l_1$$\end{document}l1-norm regularized least square optimization approach to find its parameter values. Estimates provided by StaRTrEK are validated using simulated data and three independent experimental datasets of two short (6 samples) and one long (48 samples) time-courses.

Conclusions

We believe that our algorithm can be used as a fast valuable computational complement to time-course experimental gene expression studies by adding a relevant kinetic property, i.e. the RNA half-life, with a strong biological interpretation, thus providing a dynamic picture of what is going in a cell during the biological process under study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04730-x.

Modulation of RNA stability regulates gene expression in two opposite ways: through buffering of RNA levels upon global perturbations and by supporting adapted differential expression

The steady state levels of RNAs, often referred to as expression levels, result from a well-balanced combination of RNA transcription and decay. Alterations in RNA levels will therefore result from tight regulation of transcription rates, decay rates or both. Here, we explore the role of RNA stability in achieving balanced gene expression and present genome-wide RNA stabilities in Drosophila melanogaster male and female cells as well as male cells depleted of proteins essential for dosage compensation. We identify two distinct RNA-stability mediated responses involved in regulation of gene expression. The first of these responds to acute and global changes in transcription and thus counteracts potentially harmful gene mis-expression by shifting the RNA stability in the direction opposite to the transcriptional change. The second response enhances inter-individual differential gene expression by adjusting the RNA stability in the same direction as a transcriptional change. Both mechanisms are global, act on housekeeping as well as non-housekeeping genes and were observed in both flies and mammals. Additionally, we show that, in contrast to mammals, modulation of RNA stability does not detectably contribute to dosage compensation of the sex-chromosomes in D. melanogaster.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3′-untranslated region (3′-UTR) bioparts are limited. Thus, transcript 3′-ends require further investigation to understand the underlying regulatory role and applications of the 3′-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3′-UTR regulatory functions and to provide a diverse collection of tunable 3′-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3′-end positions revealed multiple 3′-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3′-UTR bioparts is advantageous over promoter- or 5′-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3′-UTR engineering in synthetic biology applications.

Genetic variants associated mRNA stability in lung

Background

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs).

Results

Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08405-y.

Combinations of slow-translating codon clusters can increase mRNA half-life in Saccharomyces cerevisiae

The presence of a single cluster of nonoptimal codons was found to decrease a transcript's half-life through the interaction of the ribosome-associated quality control machinery with stalled ribosomes in Saccharomyces cerevisiae The impact of multiple nonoptimal codon clusters on a transcript's half-life, however, is unknown. Using a kinetic model, we predict that inserting a second nonoptimal cluster near the 5' end can lead to synergistic effects that increase a messenger RNA's (mRNA's) half-life in S. cerevisiae Specifically, the 5' end cluster suppresses the formation of ribosome queues, reducing the interaction of ribosome-associated quality control factors with stalled ribosomes. We experimentally validate this prediction by introducing two nonoptimal clusters into three different genes and find that their mRNA half-life increases up to fourfold. The model also predicts that in the presence of two clusters, the cluster closest to the 5' end is the primary determinant of mRNA half-life. These results suggest the "translational ramp," in which nonoptimal codons are located near the start codon and increase translational efficiency, may have the additional biological benefit of allowing downstream slow-codon clusters to be present without decreasing mRNA half-life. These results indicate that codon usage bias plays a more nuanced role in controlling cellular protein levels than previously thought.

Replication stress inhibits synthesis of histone mRNAs in yeast by removing Spt10p and Spt21p from the histone promoters

Proliferating cells coordinate histone and DNA synthesis to maintain correct stoichiometry for chromatin assembly. Histone mRNA levels must be repressed when DNA replication is inhibited to prevent toxicity and genome instability due to free non-chromatinized histone proteins. In mammalian cells, replication stress triggers degradation of histone mRNAs, but it is unclear if this mechanism is conserved from other species. The aim of this study was to identify the histone mRNA decay pathway in the yeast Saccharomyces cerevisiae and determine the mechanism by which DNA replication stress represses histone mRNAs. Using reverse transcription-quantitative PCR and chromatin immunoprecipitation–quantitative PCR, we show here that histone mRNAs can be degraded by both 5′ → 3′ and 3′ → 5′ pathways; however, replication stress does not trigger decay of histone mRNA in yeast. Rather, replication stress inhibits transcription of histone genes by removing the histone gene–specific transcription factors Spt10p and Spt21p from histone promoters, leading to disassembly of the preinitiation complexes and eviction of RNA Pol II from histone genes by a mechanism facilitated by checkpoint kinase Rad53p and histone chaperone Asf1p. In contrast, replication stress does not remove SCB-binding factor transcription complex, another activator of histone genes, from the histone promoters, suggesting that Spt10p and Spt21p have unique roles in the transcriptional downregulation of histone genes during replication stress. Together, our data show that, unlike in mammalian cells, replication stress in yeast does not trigger decay of histone mRNAs but inhibits histone transcription.

The RNA-binding protein Puf5 contributes to buffering of mRNA upon chromatin-mediated changes in nascent transcription

Gene expression involves regulation of chromatin structure and transcription, as well as processing of the transcribed mRNA. While there are feedback mechanisms, it is not clear whether these include crosstalk between chromatin architecture and mRNA decay. To address this, we performed a genome-wide genetic screen using a Saccharomyces cerevisiae strain harbouring the H3K56A mutation, which is known to perturb chromatin structure and nascent transcription. We identified Puf5 (also known as Mpt5) as essential in an H3K56A background. Depletion of Puf5 in this background leads to downregulation of Puf5 targets. We suggest that Puf5 plays a role in post-transcriptional buffering of mRNAs, and support this by transcriptional shutoff experiments in which Puf5 mRNA targets are degraded slower in H3K56A cells compared to wild-type cells. Finally, we show that post-transcriptional buffering of Puf5 targets is widespread and does not occur only in an H3K56A mutant, but also in an H3K4R background, which leads to a global increase in nascent transcription. Our data suggest that Puf5 determines the fate of its mRNA targets in a context-dependent manner acting as an mRNA surveillance hub balancing deregulated nascent transcription to maintain physiological mRNA levels.

Genetic and pharmacological evidence for kinetic competition between alternative poly(A) sites in yeast

Most eukaryotic mRNAs accommodate alternative sites of poly(A) addition in the 3’ untranslated region in order to regulate mRNA function. Here, we present a systematic analysis of 3’ end formation factors, which revealed 3’UTR lengthening in response to a loss of the core machinery, whereas a loss of the Sen1 helicase resulted in shorter 3’UTRs. We show that the anti-cancer drug cordycepin, 3’ deoxyadenosine, caused nucleotide accumulation and the usage of distal poly(A) sites. Mycophenolic acid, a drug which reduces GTP levels and impairs RNA polymerase II (RNAP II) transcription elongation, promoted the usage of proximal sites and reversed the effects of cordycepin on alternative polyadenylation. Moreover, cordycepin-mediated usage of distal sites was associated with a permissive chromatin template and was suppressed in the presence of an rpb1 mutation, which slows RNAP II elongation rate. We propose that alternative polyadenylation is governed by temporal coordination of RNAP II transcription and 3’ end processing and controlled by the availability of 3’ end factors, nucleotide levels and chromatin landscape.

A prion-like domain of Tpk2 catalytic subunit of protein kinase A modulates P-body formation in response to stress in budding yeast

Low complexity regions are involved in the assembly and disassembly of P-bodies (PBs). Saccharomyces cerevisiae contains three genes encoding the protein kinase A (PKA) catalytic subunit: TPK1, TPK2 and TPK3. Tpk2 and Tpk3 isoforms localize to PBs upon glucose starvation showing different mechanisms and kinetics of accumulation. In contrast to the other two isoforms, Tpk2 harbors a glutamine-rich prion-like domain (PrLD) at the N-terminus. Here we show that the appearance of Tpk2 foci in response to glucose starvation, heat stress or stationary phase was dependent on its PrLD. Moreover, the PrLD of Tpk2 was necessary for efficient PB and stress granule aggregation during stress conditions and in quiescent cells. Deletion of PrLD does not affect the in vitro and in vivo kinase activity of Tpk2 or its interaction with the regulatory subunit Bcy1. We present evidence that the PrLD of Tpk2 serves as a scaffold domain for PB assembly in a manner that is independent of Pat1 phosphorylation by PKA. In addition, a mutant strain where Tpk2 lacks PrLD showed a decrease of turnover of mRNA during glucose starvation. This work therefore provides new insight into the mechanism of stress-induced cytoplasmic mRNP assembly, and the role of isoform specific domains in the regulation of PKA catalytic subunit specificity and dynamic localization to cytoplasmic RNPs granules.

Genome-scale molecular principles of mRNA half-life regulation in yeast

Precise control of protein and messenger RNA (mRNA) degradation is essential for cellular metabolism and homeostasis. Controlled and specific degradation of both molecular species necessitates their engagements with the respective degradation machineries; this engagement involves a disordered/unstructured segment of the substrate traversing the degradation tunnel of the machinery and accessing the catalytic sites. However, while molecular factors influencing protein degradation have been extensively explored on a genome scale, and in multiple organisms, such a comprehensive understanding remains missing for mRNAs. Here, we analyzed multiple genome-scale experimental yeast mRNA half-life data in light of experimentally derived mRNA secondary structures and protein binding data, along with high-resolution X-ray crystallographic structures of the RNase machines. Results unraveled a consistent genome-scale trend that mRNAs comprising longer terminal and/or internal unstructured segments have significantly shorter half-lives; the lengths of the 5'-terminal, 3'-terminal, and internal unstructured segments that affect mRNA half-life are compatible with molecular structures of the 5' exo-, 3' exo-, and endoribonuclease machineries. Sequestration into ribonucleoprotein complexes elongates mRNA half-life, presumably by burying ribonuclease engagement sites under oligomeric interfaces. After gene duplication, differences in terminal unstructured lengths, proportions of internal unstructured segments, and oligomerization modes result in significantly altered half-lives of paralogous mRNAs. Side-by-side comparison of molecular principles underlying controlled protein and mRNA degradation in yeast unravels their remarkable mechanistic similarities and suggests how the intrinsic structural features of the two molecular species, at two different levels of the central dogma, regulate their half-lives on genome scale.

Metabolic RNA labeling for probing RNA dynamics in bacteria

Metabolic labeling of RNAs with noncanonical nucleosides that are chemically active, followed by chemoselective conjugation with imaging probes or enrichment tags, has emerged as a powerful method for studying RNA transcription and degradation in eukaryotes. However, metabolic RNA labeling is not applicable for prokaryotes, in which the complexity and distinctness of gene regulation largely remain to be explored. Here, we report 2'-deoxy-2'-azidoguanosine (AzG) as a noncanonical nucleoside compatible with metabolic labeling of bacterial RNAs. With AzG, we develop AIR-seq (azidonucleoside-incorporated RNA sequencing), which enables genome-wide analysis of transcription upon heat stress in Escherichia coli. Furthermore, AIR-seq coupled with pulse-chase labeling allows for global analysis of bacterial RNA degradation. Finally, we demonstrate that RNAs of mouse gut microbiotas can be metabolically labeled with AzG in living animals. The AzG-enabled metabolic RNA labeling should find broad applications in studying RNA biology in various bacterial species.

The Arabidopsis thaliana mRNA decay factor PAT1 functions in osmotic stress responses and decaps ABA-responsive genes

mRNA decapping plays essential roles in regulating gene expression during cellular reprogramming in response to developmental and environmental cues. The evolutionarily conserved PAT1 proteins activate decapping by binding mRNA, recruiting other decapping components, and promoting processing body (PB) assembly. Arabidopsis encodes 3 PAT proteins: PAT1, PATH1, and PATH2. Here, we report that only pat1 mutants exhibit hypersensitivity to ABA and that transcripts of ABA-responsive genes, but not those of ABA biosynthesis genes, persist longer in these mutants. The pat1 mutants also exhibit increased resistance to drought stress and resistance to Pythium irregulare. This is supported by assays showing that PAT1 functions specifically in decapping of the canonical ABA-responsive gene COR15A. In summary, PAT1 protein mediates decay of ABA-responsive genes and, thus, regulates stress responses.

Global Analysis of the Human RNA Degradome Reveals Widespread Decapped and Endonucleolytic Cleaved Transcripts

RNA decay is an important regulatory mechanism for gene expression at the posttranscriptional level. Although the main pathways and major enzymes that facilitate this process are well defined, global analysis of RNA turnover remains under-investigated. Recent advances in the application of next-generation sequencing technology enable its use in order to examine various RNA decay patterns at the genome-wide scale. In this study, we investigated human RNA decay patterns using parallel analysis of RNA end-sequencing (PARE-seq) data from XRN1-knockdown HeLa cell lines, followed by a comparison of steady state and degraded mRNA levels from RNA-seq and PARE-seq data, respectively. The results revealed 1103 and 1347 transcripts classified as stable and unstable candidates, respectively. Of the unstable candidates, we found that a subset of the replication-dependent histone transcripts was polyadenylated and rapidly degraded. Additionally, we identified 380 endonucleolytically cleaved candidates by analyzing the most abundant PARE sequence on a transcript. Of these, 41.4% of genes were classified as unstable genes, which implied that their endonucleolytic cleavage might affect their mRNA stability. Furthermore, we identified 1877 decapped candidates, including HSP90B1 and SWI5, having the most abundant PARE sequences at the 5′-end positions of the transcripts. These results provide a useful resource for further analysis of RNA decay patterns in human cells.

Salt-Induced Stability of SR1/CAMTA3 mRNA Is Mediated by Reactive Oxygen Species and Requires the 3' End of Its Open Reading Frame

Soil salinity, a prevalent abiotic stress, causes enormous losses in global crop yields annually. Previous studies have shown that salt stress-induced reprogramming of gene expression contributes to the survival of plants under this stress. However, mechanisms regulating gene expression in response to salt stress at the posttranscriptional level are not well understood. In this study, we show that salt stress increases the level of Signal Responsive 1 (SR1) mRNA, a member of signal-responsive Ca2+/calmodulin-regulated transcription factors, by enhancing its stability. We present multiple lines of evidence indicating that reactive oxygen species generated by NADPH oxidase activity mediate salt-induced SR1 transcript stability. Using mutants impaired in either nonsense-mediated decay, XRN4 or mRNA decapping pathways, we show that neither the nonsense-mediated mRNA decay pathway, XRN4 nor the decapping of SR1 mRNA is required for its decay. We analyzed the salt-induced accumulation of eight truncated versions of the SR1 coding region (∼3 kb) in the sr1 mutant background. This analysis identified a 500-nt region at the 3' end of the SR1 coding region to be required for the salt-induced stability of SR1 mRNA. Potential mechanisms by which this region confers SR1 transcript stability in response to salt are discussed.

In Vivo RNAi-Mediated eIF3m Knockdown Affects Ribosome Biogenesis and Transcription but Has Limited Impact on mRNA-Specific Translation

Translation is an essential biological process, and dysregulation is associated with a range of diseases including ribosomopathies, diabetes, and cancer. Here, we examine translation dysregulation in vivo using RNAi to knock down the m-subunit of the translation initiation factor eIF3 in the mouse liver. Transcriptome sequencing, ribosome profiling, whole proteome, and phosphoproteome analyses show that eIF3m deficiency leads to the transcriptional response and changes in cellular translation that yield few detectable differences in the translation of particular mRNAs. The transcriptional response fell into two main categories: ribosome biogenesis (increased transcription of ribosomal proteins) and cell metabolism (alterations in lipid, amino acid, nucleic acid, and drug metabolism). Analysis of ribosome biogenesis reveals inhibition of rRNA processing, highlighting decoupling of rRNA synthesis and ribosomal protein gene transcription in response to eIF3m knockdown. Interestingly, a similar reduction in eIF3m protein levels is associated with induction of the mTOR pathway in vitro but not in vivo. Overall, this work highlights the utility of a RNAi-based in vivo approach for studying the regulation of mammalian translation in vivo.

Using TTchem-seq for profiling nascent transcription and measuring transcript elongation

The dynamics of transcription can be studied genome wide by high-throughput sequencing of nascent and newly synthesized RNA. 4-thiouridine (4SU) labeling in vivo enables the specific capture of such new transcripts, with 4SU residues being tagged by biotin linkers and captured using streptavidin beads before library production and high-throughput sequencing. To achieve high-resolution profiles of transcribed regions, an RNA fragmentation step before biotin tagging was introduced, in an approach known as transient transcriptome sequencing (TT-seq). We recently introduced a chemical approach for RNA fragmentation that we refer to as TTchem-seq. We describe how TTchem-seq can be used in combination with transient inhibition of early elongation using the reversible CDK9 inhibitor, 5,6-dichlorobenzimidazole 1-β-D-ribofuranoside (DRB), to measure RNA polymerase II (RNAPII) elongation rates in vivo, a technique we call DRB/TTchem-seq. Here, we provide detailed protocols for carrying out TTchem-seq and DRB/TTchem-seq, including computational analysis. Experiments and data analysis can be performed over a period of 10-13 d and require molecular biology and bioinformatics skills.

Two roles for the yeast transcription coactivator SAGA and a set of genes redundantly regulated by TFIID and SAGA

Deletions within genes coding for subunits of the transcription coactivator SAGA caused strong genome-wide defects in transcription and SAGA-mediated chromatin modifications. In contrast, rapid SAGA depletion produced only modest transcription defects at 13% of protein-coding genes - genes that are generally more sensitive to rapid TFIID depletion. However, transcription of these 'coactivator-redundant' genes is strongly affected by rapid depletion of both factors, showing the overlapping functions of TFIID and SAGA at this gene set. We suggest that this overlapping function is linked to TBP-DNA recruitment. The remaining 87% of expressed genes that we term 'TFIID-dependent' are highly sensitive to rapid TFIID depletion and insensitive to rapid SAGA depletion. Genome-wide mapping of SAGA and TFIID found binding of both factors at many genes independent of gene class. Promoter analysis suggests that the distinction between the gene classes is due to multiple components rather than any single regulatory factor or promoter sequence motif.

Alterations of signaling pathways in response to chemical perturbations used to measure mRNA decay rates in yeast

Assessing variations in mRNA stability typically involves inhibiting transcription either globally or in a gene-specific manner. Alternatively, mRNA pulse-labeling strategies offer a means to calculate mRNA stability without inhibiting transcription. However, key stress-responsive cell signaling pathways, which affect mRNA stability, may themselves be perturbed by the approaches used to measure mRNA stability, leading to artifactual results. Here, we have focused on common strategies to measure mRNA half-lives in yeast and determined that commonly used transcription inhibitors thiolutin and 1,10 phenanthroline inhibit TORC1 signaling, PKC signaling, and partially activate HOG signaling. Additionally, 4-thiouracil (4tU), a uracil analog used in mRNA pulse-labeling approaches, modestly induces P-bodies, mRNA-protein granules implicated in storage and decay of nontranslating mRNA. Thiolutin also induces P-bodies, whereas phenanthroline has no effect. Doxycycline, which controls "Tet On/Tet Off" regulatable promoters, shows no impact on the above signaling pathways or P-bodies. In summary, our data argues that broad-acting transcriptional inhibitors are problematic for determining mRNA half-life, particularly if studying the impacts of the TORC1, HOG, or PKC pathway on mRNA stability. Regulatable promoter systems are a preferred approach for individual mRNA half-life studies, with 4tU labeling representing a good approach to global mRNA half-life analysis, despite modestly inducing P-bodies.

Mip6 binds directly to the Mex67 UBA domain to maintain low levels of Msn2/4 stress-dependent mRNAs

RNA-binding proteins (RBPs) participate in all steps of gene expression, underscoring their potential as regulators of RNA homeostasis. We structurally and functionally characterize Mip6, a four-RNA recognition motif (RRM)-containing RBP, as a functional and physical interactor of the export factor Mex67. Mip6-RRM4 directly interacts with the ubiquitin-associated (UBA) domain of Mex67 through a loop containing tryptophan 442. Mip6 shuttles between the nucleus and the cytoplasm in a Mex67-dependent manner and concentrates in cytoplasmic foci under stress. Photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation experiments show preferential binding of Mip6 to mRNAs regulated by the stress-response Msn2/4 transcription factors. Consistent with this binding, MIP6 deletion affects their export and expression levels. Additionally, Mip6 interacts physically and/or functionally with proteins with a role in mRNA metabolism and transcription such as Rrp6, Xrn1, Sgf73, and Rpb1. These results reveal a novel role for Mip6 in the homeostasis of Msn2/4-dependent transcripts through its direct interaction with the Mex67 UBA domain.

Attenuated Codon Optimality Contributes to Neural-Specific mRNA Decay in Drosophila

Tissue-specific mRNA stability is important for cell fate and physiology, but the mechanisms involved are not fully understood. We found that zygotic mRNA stability in Drosophila correlates with codon content: optimal codons are enriched in stable transcripts associated with metabolic functions like translation, while non-optimal codons are enriched in unstable transcripts, including those associated with neural development. Bioinformatic analyses and reporter assays revealed that similar codons stabilize or destabilize mRNAs in the nervous system and other tissues, but the link between codon content and stability is attenuated in the nervous system. We confirmed that optimal codons are decoded by abundant tRNAs while non-optimal codons are decoded by less abundant tRNAs in embryos and in the nervous system. We conclude that codon optimality is a general determinant of zygotic mRNA stability, and attenuation of codon optimality allows trans-acting factors to exert greater influence over mRNA decay in the nervous system.

Burow et al. report that codon optimality is a general determinant of zygotic mRNA stability in Drosophila embryos, but the link between codons and stability is weak in the nervous system. Bioinformatics, reporter transcript assays, and tRNA quantitation show that the attenuation of codon optimality establishes neuralspecific mRNA decay.

Big data in yeast systems biology

Systems biology uses computational and mathematical modeling to study complex interactions in a biological system. The yeast Saccharomyces cerevisiae, which has served as both an important model organism and cell factory, has pioneered both the early development of such models and modeling concepts, and the more recent integration of multi-omics big data in these models to elucidate fundamental principles of biology. Here, we review the advancement of big data technologies to gain biological insight in three aspects of yeast systems biology: gene expression dynamics, cellular metabolism and the regulation network between gene expression and metabolism. The role of big data and complementary modeling approaches, including the expansion of genome-scale metabolic models and machine learning methodologies, are discussed as key drivers in the rapid advancement of yeast systems biology.

Codon stabilization coefficient as a metric to gain insights into mRNA stability and codon bias and their relationships with translation

The codon stabilization coefficient (CSC) is derived from the correlation between each codon frequency in transcripts and mRNA half-life experimental data. In this work, we used this metric as a reference to compare previously published Saccharomyces cerevisiae mRNA half-life datasets and investigate how codon composition related to protein levels. We generated CSCs derived from nine studies. Four datasets produced similar CSCs, which also correlated with other independent parameters that reflected codon optimality, such as the tRNA abundance and ribosome residence time. By calculating the average CSC for each gene, we found that most mRNAs tended to have more non-optimal codons. Conversely, a high proportion of optimal codons was found for genes coding highly abundant proteins, including proteins that were only transiently overexpressed in response to stress conditions. We also used CSCs to identify and locate mRNA regions enriched in non-optimal codons. We found that these stretches were usually located close to the initiation codon and were sufficient to slow ribosome movement. However, in contrast to observations from reporter systems, we found no position-dependent effect on the mRNA half-life. These analyses underscore the value of CSCs in studies of mRNA stability and codon bias and their relationships with protein expression.

Metabolic response to Parkinson's disease recapitulated by the haploinsufficient diploid yeast cells hemizygous for the adrenodoxin reductase gene

Adrenodoxin reductase, a widely conserved mitochondrial P450 protein, catalyses essential steps in steroid hormone biosynthesis and is highly expressed in the adrenal cortex. The yeast adrenodoxin reductase homolog, Arh1p, is involved in cytoplasmic and mitochondrial iron homeostasis and is required for activity of enzymes containing an Fe-S cluster. In this paper, we investigated the response of yeast to the loss of a single copy of ARH1, an oxidoreductase of the mitochondrial inner membrane, which is among the few mitochondrial proteins that is essential for viability in yeast. The phenotypic, transcriptional, proteomic, and metabolic landscape indicated that Saccharomyces cerevisiae successfully adapted to this loss, displaying an apparently dosage-insensitive cellular response. However, a considered investigation of transcriptional regulation in ARH1-impaired yeast highlighted that a significant hierarchical reorganisation occurred, involving the iron assimilation and tyrosine biosynthetic processes. The interconnected roles of the iron and tyrosine pathways, coupled with oxidative processes, are of interest beyond yeast since they are involved in dopaminergic neurodegeneration associated with Parkinson's disease. The identification of similar responses in yeast, albeit preliminary, suggests that this simple eukaryote could have potential as a model system for investigating the regulatory mechanisms leading to the initiation and progression of early disease responses in humans.

Condition-Specific Modeling of Biophysical Parameters Advances Inference of Regulatory Networks

Large-scale inference of eukaryotic transcription-regulatory networks remains challenging. One underlying reason is that existing algorithms typically ignore crucial regulatory mechanisms, such as RNA degradation and post-transcriptional processing. Here, we describe InfereCLaDR, which incorporates such elements and advances prediction in Saccharomyces cerevisiae. First, InfereCLaDR employs a high-quality Gold Standard dataset that we use separately as prior information and for model validation. Second, InfereCLaDR explicitly models transcription factor activity and RNA half-lives. Third, it introduces expression subspaces to derive condition-responsive regulatory networks for every gene. InfereCLaDR’s final network is validated by known data and trends and results in multiple insights. For example, it predicts long half-lives for transcripts of the nucleic acid metabolism genes and members of the cytosolic chaperonin complex as targets of the proteasome regulator Rpn4p. InfereCLaDR demonstrates that more biophysically realistic modeling of regulatory networks advances prediction accuracy both in eukaryotes and prokaryotes.

Role of DCP1-DCP2 complex regulated by viral and host microRNAs in DNA virus infection

The DCP1-DCP2 complex can regulate the antiviral immunity of animals by the decapping of retrovirus RNAs and the suppression of RNAi during RNA virus infection. However, the influence of DCP1-DCP2 complex on DNA virus infection and the regulation of DCP1-DCP2 complex by microRNAs (miRNAs) remain unclear. In this study, the role of miRNA-regulated DCP1-DCP2 complex in DNA virus infection was characterized. Our results showed that the DCP1-DCP2 complex played a positive role in the infection of white spot syndrome virus (WSSV), a DNA virus of shrimp. In the DCP1-DCP2 complex, the N-terminal regulatory domain of DCP2 was interacted with the EVH1 domain of DCP1. Furthermore, shrimp miRNA miR-87 inhibited WSSV infection by targeting the host DCP2 gene and viral miRNA WSSV-miR-N46 took a negative effect on WSSV replication by targeting the host DCP1 gene. Therefore, our study provided novel insights into the underlying mechanism of DCP1-DCP2 complex and its regulation by miRNAs in virus-host interactions. IMPORTANCE: During RNA virus infection, the DCP1-DCP2 complex can play important roles in the animal antiviral immunity by decapping retrovirus RNAs and suppressing RNAi. In the present study, the findings indicated that the silencing of DCP1 and DCP2 inhibited the infection of WSSV, a DNA virus of shrimp, suggesting that the DCP1-DCP2 complex facilitated DNA virus infection. Due to the suppressive role of the DCP1-DCP2 complex in shrimp RNAi against WSSV infection, the DCP1-DCP2 complex could promote WSSV infection in shrimp. The results showed that WSSV-miR-N46 and shrimp miR-87 could respectively suppress the expressions of DCP1 and DCP2 to affect virus infection. Therefore, our study contributed novel aspects of the DCP1-DCP2 complex and its regulation by miRNAs in virus-host interactions.

Methods for the analysis of transcriptome dynamics

The transcriptome is the complete set of transcripts in a cell or tissue and includes ribosomal RNA (rRNA), messenger RNA (mRNA), transfer RNA (tRNA), and regulatory noncoding RNA. At steady-state, the transcriptome results from a compensatory variation of the transcription and decay rate to maintain the RNA concentration constant. RNA transcription constitutes the first stage in gene expression, and thus is a major and primary mode of gene expression control. Nevertheless, regulation of RNA decay is also a key factor in gene expression control, involving either selective RNA stabilization or enhanced degradation. Transcriptome analysis allows the identification of gene expression alterations, providing new insights regarding the pathways and mechanisms involved in physiological and pathological processes. Upon perturbation of cell homeostasis, rapid changes in gene expression are required to adapt to new conditions. Thus, to better understand the regulatory mechanisms associated with gene expression alterations, it is vital to acknowledge the relative contribution of RNA synthesis and decay to the transcriptome. To the toxicology field, the study of gene expression regulation mechanisms can help identify the early and mechanistic relevant cellular events associated with a particular response. This review aims to provide a critical comparison of the available methods used to analyze the contribution of RNA transcription and decay to gene expression dynamics. Notwithstanding, an integration of the data obtained is necessary to understand the entire repercussions of gene transcription changes at a system-level. Thus, a brief overview of the methods available for the integration and analysis of the data obtained from transcriptome analysis will also be provided.

Ksp1-dependent phosphorylation of eIF4G modulates post-transcriptional regulation of specific mRNAs under glucose deprivation conditions

Post-transcriptional regulation is an important mechanism for modulating gene expression and is performed by numerous mRNA-binding proteins. To understand the mechanisms underlying post-transcriptional regulation, we investigated the phosphorylation status of 32 mRNA-binding proteins under glucose deprivation conditions in Saccharomyces cerevisiae. We identified 17 glucose-sensitive phosphoproteins and signal pathways implicated in their phosphorylation. Notably, phosphorylation of the eukaryotic translation initiation factor 4G (eIF4G) was regulated by both the Snf1/AMPK pathway and the target of rapamycin complex 1 (TORC1) pathway. The serine/threonine protein kinase Ksp1 has previously been suggested to be a downstream effector of TORC1, but its detailed function has rarely been discussed. We identified that Snf1/AMPK and TORC1 signalings converge on Ksp1, which phosphorylates eIF4G under glucose deprivation conditions. Ksp1-dependent phosphorylation of eIF4G regulates the degradation of specific mRNAs (e.g. glycolytic mRNAs and ribosomal protein mRNAs) under glucose deprivation conditions likely through the recruitment of Dhh1. Taken together, our results suggest that Ksp1 functions as a novel modulator of post-transcriptional regulation in yeast.

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

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

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

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

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

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

Brr6 plays a role in gene recruitment and transcriptional regulation at the nuclear envelope

Correlation between transcriptional regulation and positioning of genes at the nuclear envelope is well established in eukaryotes, but the mechanisms involved are not well understood. We show that brr6-1, a mutant of the essential yeast envelope transmembrane protein Brr6p, impairs normal positioning and expression of the PAB1 and FUR4-GAL1,10,7 loci. Similarly, expression of a dominant negative nucleoplasmic Brr6 fragment in wild-type cells reproduced many of the brr6-1 effects. Histone chromatin immunoprecipitation (ChIP) experiments showed decreased acetylation at the key histone H4K16 residue in the FUR4-GAL1,10,7 region in brr6-1. Importantly, blocking deacetylation significantly suppressed selected brr6-1 phenotypes. ChIPseq with FLAG-tagged Brr6 fragments showed enrichment at FUR4 and several other genes that showed striking changes in brr6-1 RNAseq data. These associations depended on a Brr6 putative zinc finger domain. Importantly, artificially tethering the GAL1 locus to the envelope suppressed the brr6-1 effects on GAL1 and FUR4 expression and increased H4K16 acetylation between GAL1 and FUR4 in the mutant. Together these results argue that Brr6 interacts with chromatin, helping to maintain normal chromatin architecture and transcriptional regulation of certain loci at the nuclear envelope.

Global role for coactivator complexes in RNA polymerase II transcription

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