Stress-induced nuclear RNA degradation pathways regulate yeast bromodomain factor 2 to promote cell survival

Identification of acidic stress-responsive genes and acid tolerance engineering in Synechococcus elongatus PCC 7942

Cyanobacteria are excellent autotrophic photosynthetic chassis employed in synthetic biology, and previous studies have suggested that they have alkaline tolerance but low acid tolerance, significantly limiting their productivity as photosynthetic chassis and necessitating investigations into the acid stress resistance mechanism. In this study, differentially expressed genes were obtained by RNA sequencing-based comparative transcriptomic analysis under long-term acidic stress conditions and acidic shock treatment, in the model cyanobacterium Synechococcus elongatus PCC 7942. A pathway enrichment analysis revealed the upregulated and downregulated pathways during long-term acidic and shock stress treatment. The subsequent single gene knockout and phenotype analysis showed that under acidic stress conditions, the strains with chlL, chlN, pex, synpcc7942_2038, synpcc7942_1890, or synpcc7942_2547 knocked out grew worse than the wild type, suggesting their involvement in acid tolerance. This finding was further confirmed by introducing the corresponding genes back into the knockout mutant individually. Moreover, individual overexpression of the chlL and chlN genes in the wild type successfully improved the tolerance of S. elongatus PCC 7942 to acidic stress. This work successfully identified six genes involved in acidic stress responses, and overexpressing chIL or chIN individually successfully improved acid tolerance in S. elongatus PCC 7942, providing valuable information to better understand the acid resistance mechanism in S. elongatus PCC 7942 and novel insights into the robustness and tolerance engineering of cyanobacterial chassis. KEY POINTS: • DEGs were identified by RNA-seq based transcriptomics analysis in response to acidic stress in S. elongatus PCC 7942. • Six genes were identified to be involved in acid tolerance in S. elongatus PCC 7942. • Overexpression of chIL or chIN individually successfully improved the acid tolerance of S. elongatus PCC 7942.

Mechanisms and regulation of spliceosome-mediated pre-mRNA splicing in Saccharomyces cerevisiae

Pre-mRNA splicing, the removal of introns and ligation of flanking exons, is a crucial step in eukaryotic gene expression. The spliceosome, a macromolecular complex made up of five small nuclear RNAs (snRNAs) and dozens of proteins, assembles on introns via a complex pathway before catalyzing the two transesterification reactions necessary for splicing. All of these steps have the potential to be highly regulated to ensure correct mRNA isoform production for proper cellular function. While Saccharomyces cerevisiae (yeast) has a limited set of intron-containing genes, many of these genes are highly expressed, resulting in a large number of transcripts in a cell being spliced. As a result, splicing regulation is of critical importance for yeast. Just as in humans, yeast splicing can be influenced by protein components of the splicing machinery, structures and properties of the pre-mRNA itself, or by the action of trans-acting factors. It is likely that further analysis of the mechanisms and pathways of splicing regulation in yeast can reveal general principles applicable to other eukaryotes. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Splicing factor Prp18p promotes genome-wide fidelity of consensus 3'-splice sites

The fidelity of splice site selection is critical for proper gene expression. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is challenging considering the low complexity of the 3'SS consensus sequence YAG. Here, we show that absence of the Prp18p splicing factor results in genome-wide activation of alternative 3'SS in S. cerevisiae, including highly unusual non-YAG sequences. Usage of these non-canonical 3'SS in the absence of Prp18p is enhanced by upstream poly(U) tracts and by their potential to interact with the first intronic nucleoside, allowing them to dock in the spliceosome active site instead of the normal 3'SS. The role of Prp18p in 3'SS fidelity is facilitated by interactions with Slu7p and Prp8p, but cannot be fulfilled by Slu7p, identifying a unique role for Prp18p in 3'SS fidelity. This fidelity function is synergized by the downstream proofreading activity of the Prp22p helicase, but is independent from another late splicing helicase, Prp43p. Our results show that spliceosomes exhibit remarkably relaxed 3'SS sequence usage in the absence of Prp18p and identify a network of spliceosomal interactions centered on Prp18p which are required to promote the fidelity of the recognition of consensus 3'SS sequences.

Spliceosomal mutations decouple 3' splice site fidelity from cellular fitness

The fidelity of splice site selection is thought to be critical for proper gene expression and cellular fitness. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is a daunting task considering the low complexity of the 3'SS consensus sequence YAG. Here we show that inactivating the near-essential splicing factor Prp18p results in a global activation of alternative 3'SS, many of which harbor sequences that highly diverge from the YAG consensus, including some highly unusual non-AG 3'SS. We show that the role of Prp18p in 3'SS fidelity is promoted by physical interactions with the essential splicing factors Slu7p and Prp8p and synergized by the proofreading activity of the Prp22p helicase. Strikingly, structure-guided point mutations that disrupt Prp18p-Slu7p and Prp18p-Prp8p interactions mimic the loss of 3'SS fidelity without any impact on cellular growth, suggesting that accumulation of incorrectly spliced transcripts does not have a major deleterious effect on cellular viability. These results show that spliceosomes exhibit remarkably relaxed fidelity in the absence of Prp18p, and that new 3'SS sampling can be achieved genome-wide without a major negative impact on cellular fitness, a feature that could be used during evolution to explore new productive alternative splice sites.

Genetic suppressors of Δgrx3 Δgrx4, lacking redundant multidomain monothiol yeast glutaredoxins, rescue growth and iron homeostasis

Saccharomyces cerevisiae Grx3 and Grx4 are multidomain monothiol glutaredoxins that are redundant with each other. They can be efficiently complemented by heterologous expression of their mammalian ortholog, PICOT, which has been linked to tumor development and embryogenesis. PICOT is now believed to act as a chaperone distributing Fe-S clusters, although the first link to iron metabolism was observed with its yeast counterparts. Like PICOT, yeast Grx3 and Grx4 reside in the cytosol and nucleus where they form unusual Fe-S clusters coordinated by two glutaredoxins with CGFS motifs and two molecules of glutathione. Depletion or deletion of Grx3/Grx4 leads to functional impairment of virtually all cellular iron-dependent processes and loss of cell viability, thus making these genes the most upstream components of the iron utilization system. Nevertheless, the Δgrx3/4 double mutant in the BY4741 genetic background is viable and exhibits slow but stable growth under hypoxic conditions. Upon exposure to air, growth of the double deletion strain ceases, and suppressor mutants appear. Adopting a high copy-number library screen approach, we discovered novel genetic interactions: overexpression of ESL1, ESL2, SOK1, SFP1 or BDF2 partially rescues growth and iron utilization defects of Δgrx3/4. This genetic escape from the requirement for Grx3/Grx4 has not been previously described. Our study shows that even a far-upstream component of the iron regulatory machinery (Grx3/4) can be bypassed, and cellular networks involving RIM101 pH sensing, cAMP signaling, mTOR nutritional signaling, or bromodomain acetylation, may confer the bypassing activities.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

Stress-induced inhibition of mRNA export triggers RNase III-mediated decay of the BDF2 mRNA

The expression of bromodomain-containing proteins that regulate chromatin structure and accessibility must be tightly controlled to ensure the appropriate regulation of gene expression. In the yeast S. cerevisiae, Bromodomain Factor 2 (BDF2) expression is extensively regulated post-transcriptionally during stress by RNase III-mediated decay (RMD), which is triggered by cleavage of the BDF2 mRNA in the nucleus by the RNase III homolog Rnt1p. Previous studies have shown that RMD-mediated down-regulation of BDF2 is hyperactivated in osmotic stress conditions, yet the mechanisms driving the enhanced nuclear cleavage of BDF2 RNA under these conditions remain unknown. Here, we show that RMD hyperactivation can be detected in multiple stress conditions that inhibit mRNA export, and that Rnt1p remains primarily localized in the nucleus during salt stress. We show that globally inhibiting mRNA nuclear export by anchoring away mRNA biogenesis or export factors out of the nucleus can recapitulate RMD hyperactivation in the absence of stress. RMD hyperactivation requires Rnt1p nuclear localization but does not depend on the BDF2 gene endogenous promoter, and its efficiency is affected by the structure of the stem-loop cleaved by Rnt1p. Because multiple stress conditions have been shown to mediate global inhibition of mRNA export, our results suggest that the hyperactivation of RMD is primarily the result of the increased nuclear retention of the BDF2 mRNA during stress.

Rrp6 Moonlights in an RNA Exosome-Independent Manner to Promote Cell Survival and Gene Expression during Stress

The nuclear RNA exosome is essential for RNA processing and degradation. Here, we show that the exosome nuclear-specific subunit Rrp6p promotes cell survival during heat stress through the cell wall integrity (CWI) pathway, independently of its catalytic activity or association with the core exosome. Rrp6p exhibits negative genetic interactions with the Slt2/Mpk1p or Paf1p elongation factors required for expression of CWI genes during stress. Overexpression of Rrp6p or of its catalytically inactive or exosome-independent mutants can partially rescue the growth defect of the mpk1Δ mutant and stimulates expression of the Mpk1 p target gene FKS2. The rrp6Δ and mpk1Δ mutants show similarities in deficient expression of CWI genes during heat shock, and overexpression of the CWI gene HSP150 can rescue the stress-induced lethality of the mpk1Δrp6Δ mutant. These results demonstrate that Rrp6p moonlights independently from the exosome to ensure proper expression of CWI genes and to promote cell survival during stress.

Wang et al. show that Rrp6 functions with the Slt2/Mpk1 and Paf1 elongation factors for the proper expression of CWI genes during heat stress. The role of Rrp6p in promoting heat-stress-induced gene expression does not require Rrp6 catalytic activity or interaction with the nuclear RNA exosome.

Robust mapping of polyadenylated and non-polyadenylated RNA 3' ends at nucleotide resolution by 3'-end sequencing

3'-end poly(A)+ sequencing is an efficient and economical method for global measurement of mRNA levels and alternative poly(A) site usage. A common method involves oligo(dT)19V reverse-transcription (RT)-based library preparation and high-throughput sequencing with a custom primer ending in (dT)19. While the majority of library products have the first sequenced nucleotide reflect the bona fide poly(A) site (pA), a substantial fraction of sequencing reads arise from various mis-priming events. These can result in incorrect pA site calls anywhere from several nucleotides downstream to several kilobases upstream from the bona fide pA site. While these mis-priming events can be mitigated by increasing annealing stringency (e.g. increasing temperature from 37 °C to 42 °C), they still persist at an appreciable level (∼10%) and computational methods must be used to prevent artifactual calls. Here we present a bioinformatics workflow for precise mapping of poly(A)+ 3' ends and handling of artifacts due to oligo(dT) mis-priming and sample polymorphisms. We test pA site calling with three different read mapping programs (STAR, BWA, and BBMap), and show that the way in which each handles terminal mismatches and soft clipping has a substantial impact on identifying correct pA sites, with BWA requiring the least post-processing to correct artifacts. We demonstrate the use of this pipeline for mapping pA sites in the model eukaryote S. cerevisiae, and further apply this technology to non-polyadenylated transcripts by employing in vitro polyadenylation prior to library prep (IVP-seq). As proof of principle, we show that a fraction of tRNAs harbor CCU 3' tails instead of the canonical CCA tail, and globally identify 3' ends of splicing intermediates arising from inefficiently spliced transcripts.

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

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

Rice In Vivo RNA Structurome Reveals RNA Secondary Structure Conservation and Divergence in Plants

RNA secondary structure plays a critical role in gene regulation. Rice (Oryza sativa) is one of the most important food crops in the world. However, RNA structure in rice has scarcely been studied. Here, we have successfully generated in vivo Structure-seq libraries in rice. We found that the structural flexibility of mRNAs might associate with the dynamics of biological function. Higher N⁶-methyladenosine (m⁶A) modification tends to have less RNA structure in 3' UTR, whereas GC content does not significantly affect in vivo mRNA structure to maintain efficient biological processes such as translation. Comparative analysis of RNA structurome between rice and Arabidopsis revealed that higher GC content does not lead to stronger structure and less RNA structural flexibility. Moreover, we found a weak correlation between sequence and structure conservation of the orthologs between rice and Arabidopsis. The conservation and divergence of both sequence and in vivo RNA structure corresponds to diverse and specific biological processes. Our results indicate that RNA secondary structure might offer a separate layer of selection to the sequence between monocot and dicot. Therefore, our study implies that RNA structure evolves differently in various biological processes to maintain robustness in development and adaptational flexibility during angiosperm evolution.

Survival in Quiescence Requires the Euchromatic Deployment of Clr4/SUV39H by Argonaute-Associated Small RNAs

Quiescence (G0) is a ubiquitous stress response through which cells enter reversible dormancy, acquiring distinct properties including reduced metabolism, resistance to stress, and long life. G0 entry involves dramatic changes to chromatin and transcription of cells, but the mechanisms coordinating these processes remain poorly understood. Using the fission yeast, here, we track G0-associated chromatin and transcriptional changes temporally and show that as cells enter G0, their survival and global gene expression programs become increasingly dependent on Clr4/SUV39H, the sole histone H3 lysine 9 (H3K9) methyltransferase, and RNAi proteins. Notably, G0 entry results in RNAi-dependent H3K9 methylation of several euchromatic pockets, prior to which Argonaute1-associated small RNAs from these regions emerge. Overall, our data reveal another function for constitutive heterochromatin proteins (the establishment of the global G0 transcriptional program) and suggest that stress-induced alterations in Argonaute-associated sRNAs can target the deployment of transcriptional regulatory proteins to specific sequences.

Bdf1 Bromodomains Are Essential for Meiosis and the Expression of Meiotic-Specific Genes

Bromodomain and Extra-terminal motif (BET) proteins play a central role in transcription regulation and chromatin signalling pathways. They are present in unicellular eukaryotes and in this study, the role of the BET protein Bdf1 has been explored in Saccharomyces cerevisiae. Mutation of Bdf1 bromodomains revealed defects on both the formation of spores and the meiotic progression, blocking cells at the exit from prophase, before the first meiotic division. This phenotype is associated with a massive deregulation of the transcription of meiotic genes and Bdf1 bromodomains are required for appropriate expression of the key meiotic transcription factor NDT80 and almost all the Ndt80-inducible genes, including APC complex components. Bdf1 notably accumulates on the promoter of Ndt80 and its recruitment is dependent on Bdf1 bromodomains. In addition, the ectopic expression of NDT80 during meiosis partially bypasses this dependency. Finally, purification of Bdf1 partners identified two independent complexes with Bdf2 or the SWR complex, neither of which was required to complete sporulation. Taken together, our results unveil a new role for Bdf1 –working independently from its predominant protein partners Bdf2 and the SWR1 complex–as a regulator of meiosis-specific genes.

Chromatin modifying proteins play a central role in transcription regulation and chromatin signalling. In this study we investigated the functional role of the bromodomains of the chromatin protein Bdf1 during yeast gametogenesis. Our results show that the bromodomains of Bdf1 are essential for meiotic progression and the formation of mature spores. Bdf1 bromodomains are required for the expression of key meiotic genes and the master regulator NDT80. Forced expression of NDT80 can partially rescue the formation of spores when Bdf1 bromodomains are mutated. The results presented here indicate that Bdf1 forms two exclusive complexes, with Bdf2 or with the SWR complex. However, none of these complexes are required for sporulation progression. To conclude, our findings suggest that Bdf1 is a new regulator of the meiotic transcription program and of the expression of the master regulator NDT80.

Splicing-Mediated Autoregulation Modulates Rpl22p Expression in Saccharomyces cerevisiae

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

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

The regulation and functions of the nuclear RNA exosome complex

The RNA exosome complex is the most versatile RNA-degradation machine in eukaryotes. The exosome has a central role in several aspects of RNA biogenesis, including RNA maturation and surveillance. Moreover, it is emerging as an important player in regulating the expression levels of specific mRNAs in response to environmental cues and during cell differentiation and development. Although the mechanisms by which RNA is targeted to (or escapes from) the exosome are still not fully understood, general principles have begun to emerge, which we discuss in this Review. In addition, we introduce and discuss novel, previously unappreciated functions of the nuclear exosome, including in transcription regulation and in the maintenance of genome stability.

Regulation of mRNA Levels by Decay-Promoting Introns that Recruit the Exosome Specificity Factor Mmi1

In eukaryotic cells, inefficient splicing is surprisingly common and leads to the degradation of transcripts with retained introns. How pre-mRNAs are committed to nuclear decay is unknown. Here, we uncover a mechanism by which specific intron-containing transcripts are targeted for nuclear degradation in fission yeast. Sequence elements within these "decay-promoting" introns co-transcriptionally recruit the exosome specificity factor Mmi1, which induces degradation of the unspliced precursor and leads to a reduction in the levels of the spliced mRNA. This mechanism negatively regulates levels of the RNA helicase DDX5/Dbp2 to promote cell survival in response to stress. In contrast, fast removal of decay-promoting introns by co-transcriptional splicing precludes Mmi1 recruitment and relieves negative expression regulation. We propose that decay-promoting introns facilitate the regulation of gene expression. Based on the identification of multiple additional Mmi1 targets, including mRNAs, long non-coding RNAs, and sn/snoRNAs, we suggest a general role in RNA regulation for Mmi1 through transcript degradation.

Widespread alternative and aberrant splicing revealed by lariat sequencing

Alternative splicing is an important and ancient feature of eukaryotic gene structure, the existence of which has likely facilitated eukaryotic proteome expansions. Here, we have used intron lariat sequencing to generate a comprehensive profile of splicing events in Schizosaccharomyces pombe, amongst the simplest organisms that possess mammalian-like splice site degeneracy. We reveal an unprecedented level of alternative splicing, including alternative splice site selection for over half of all annotated introns, hundreds of novel exon-skipping events, and thousands of novel introns. Moreover, the frequency of these events is far higher than previous estimates, with alternative splice sites on average activated at ∼3% the rate of canonical sites. Although a subset of alternative sites are conserved in related species, implying functional potential, the majority are not detectably conserved. Interestingly, the rate of aberrant splicing is inversely related to expression level, with lowly expressed genes more prone to erroneous splicing. Although we validate many events with RNAseq, the proportion of alternative splicing discovered with lariat sequencing is far greater, a difference we attribute to preferential decay of aberrantly spliced transcripts. Together, these data suggest the spliceosome possesses far lower fidelity than previously appreciated, highlighting the potential contributions of alternative splicing in generating novel gene structures.