RNA levers and switches controlling viral gene expression.

RNA viruses are diverse and abundant pathogens that are responsible for numerous human diseases. RNA viruses possess relatively compact genomes and have therefore evolved multiple mechanisms to maximize their coding capacities, often by encoding overlapping reading frames. These reading frames are then decoded by mechanisms such as alternative splicing and ribosomal frameshifting to produce multiple distinct proteins. These solutions are enabled by the ability of the RNA genome to fold into 3D structures that can mimic cellular RNAs, hijack host proteins, and expose or occlude regulatory protein-binding motifs to ultimately control key process in the viral life cycle. We highlight recent findings focusing on less conventional mechanisms of gene expression and new discoveries on the role of RNA structures.

A ribosome-bound quality control complex triggers degradation of nascent peptides and signals translation stress.

The conserved transcriptional regulator heat shock factor 1 (Hsf1) is a key sensor of proteotoxic and other stress in the eukaryotic cytosol. We surveyed Hsf1 activity in a genome-wide loss-of-function library in Saccaromyces cerevisiae as well as ~78,000 double mutants and found Hsf1 activity to be modulated by highly diverse stresses. These included disruption of a ribosome-bound complex we named the Ribosome Quality Control Complex (RQC) comprising the Ltn1 E3 ubiquitin ligase, two highly conserved but poorly characterized proteins (Tae2 and Rqc1), and Cdc48 and its cofactors. Electron microscopy and biochemical analyses revealed that the RQC forms a stable complex with 60S ribosomal subunits containing stalled polypeptides and triggers their degradation. A negative feedback loop regulates the RQC, and Hsf1 senses an RQC-mediated translation-stress signal distinctly from other stresses. Our work reveals the range of stresses Hsf1 monitors and elucidates a conserved cotranslational protein quality control mechanism.

Author Correction: Determination of RNA structural diversity and its role in HIV-1 RNA splicing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Determination of RNA structural diversity and its role in HIV-1 RNA splicing.

Human immunodeficiency virus 1 (HIV-1) is a retrovirus with a ten-kilobase single-stranded RNA genome. HIV-1 must express all of its gene products from a single primary transcript, which undergoes alternative splicing to produce diverse protein products that include structural proteins and regulatory factors1,2. Despite the critical role of alternative splicing, the mechanisms that drive the choice of splice site are poorly understood. Synonymous RNA mutations that lead to severe defects in splicing and viral replication indicate the presence of unknown cis-regulatory elements3. Here we use dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) to investigate the structure of HIV-1 RNA in cells, and develop an algorithm that we name 'detection of RNA folding ensembles using expectation-maximization' (DREEM), which reveals the alternative conformations that are assumed by the same RNA sequence. Contrary to previous models that have analysed population averages4, our results reveal heterogeneous regions of RNA structure across the entire HIV-1 genome. In addition to confirming that in vitro characterized5 alternative structures for the HIV-1 Rev responsive element also exist in cells, we discover alternative conformations at critical splice sites that influence the ratio of transcript isoforms. Our simultaneous measurement of splicing and intracellular RNA structure provides evidence for the long-standing hypothesis6-8 that heterogeneity in RNA conformation regulates splice-site use and viral gene expression.

A Stress Response that Monitors and Regulates mRNA Structure Is Central to Cold Shock Adaptation.

Temperature influences the structural and functional properties of cellular components, necessitating stress responses to restore homeostasis following temperature shift. Whereas the circuitry controlling the heat shock response is well understood, that controlling the E. coli cold shock adaptation program is not. We found that during the growth arrest phase (acclimation) that follows shift to low temperature, protein synthesis increases, and open reading frame (ORF)-wide mRNA secondary structure decreases. To identify the regulatory system controlling this process, we screened for players required for increased translation. We identified a two-member mRNA surveillance system that enables recovery of translation during acclimation: RNase R assures appropriate mRNA degradation and the Csps dynamically adjust mRNA secondary structure to globally modulate protein expression level. An autoregulatory switch in which Csps tune their own expression to cellular demand enables dynamic control of global translation. The universality of Csps in bacteria suggests broad utilization of this control mechanism.

KREH2 helicase represses ND7 mRNA editing in procyclic-stage Trypanosoma brucei by opposite modulation of canonical and 'moonlighting' gRNA utilization creating a proposed mRNA structure.

Unknown factors regulate mitochondrial U-insertion/deletion (U-indel) RNA editing in procyclic-form (PCF) and bloodstream-form (BSF) T. brucei. This editing, directed by anti-sense gRNAs, creates canonical protein-encoding mRNAs and may developmentally control respiration. Canonical editing by gRNAs that specify protein-encoding mRNA sequences occurs amid massive non-canonical editing of unclear sources and biological significance. We found PCF-specific repression at a major early checkpoint in mRNA ND7, involving helicase KREH2-dependent opposite modulation of canonical and non-canonical 'terminator' gRNA utilization. Terminator-programmed editing derails canonical editing and installs proposed repressive structure in 30% of the ND7 transcriptome. BSF-to-PCF differentiation in vitro recreated this negative control. Remarkably, KREH2-RNAi knockdown relieved repression and increased editing progression by reverting canonical/terminator gRNA utilization. ND7 transcripts lacking early terminator-directed editing in PCF exhibited similar negative editing control along the mRNA sequence, suggesting global modulation of gRNA utilization fidelity. The terminator is a 'moonlighting' gRNA also associated with mRNA COX3 canonical editing, so the gRNA transcriptome seems multifunctional. Thus, KREH2 is the first identified repressor in developmental editing control. This and our prior work support a model whereby KREH2 activates or represses editing in a stage and substrate-specific manner. KREH2's novel dual role tunes mitochondrial gene expression in either direction during development.

Structure of pre-miR-31 reveals an active role in Dicer-TRBP complex processing.

As an essential posttranscriptional regulator of gene expression, microRNA (miRNA) levels must be strictly maintained. The biogenesis of many miRNAs is mediated by trans-acting protein partners through a variety of mechanisms, including remodeling of the RNA structure. miR-31 functions as an oncogene in numerous cancers, and interestingly, its biogenesis is not known to be regulated by protein-binding partners. Therefore, the intrinsic structural properties of the precursor element of miR-31 (pre-miR-31) can provide a mechanism by which its biogenesis is regulated. We determined the solution structure of pre-miR-31 to investigate the role of distinct structural elements in regulating processing by the Dicer-TRBP complex. We found that the presence or absence of mismatches within the helical stem does not strongly influence Dicer-TRBP processing of the pre-miRNAs. However, both the apical loop size and structure at the Dicing site are key elements for discrimination by the Dicer-TRBP complex. Interestingly, our NMR-derived structure reveals the presence of a triplet of base pairs that link the Dicer cleavage site and the apical loop. Mutational analysis in this region suggests that the stability of the junction region strongly influences processing by the Dicer-TRBP complex. Our results enrich our understanding of the active role that RNA structure plays in regulating miRNA biogenesis, which has direct implications for the control of gene expression.

Trypanosome RNA helicase KREH2 differentially controls non-canonical editing and putative repressive structure via a novel proposed 'bifunctional' gRNA in mRNA A6.

U-insertion/deletion (U-indel) RNA editing in trypanosome mitochondria is directed by guide RNAs (gRNAs). This editing may developmentally control respiration in bloodstream forms (BSF) and insect procyclic forms (PCF). Holo-editosomes include the accessory RNA Editing Substrate Binding Complex (RESC) and RNA Editing Helicase 2 Complex (REH2C), but the specific proteins controlling differential editing remain unknown. Also, RNA editing appears highly error prone because most U-indels do not match the canonical pattern. However, despite extensive non-canonical editing of unknown functions, accurate canonical editing is required for normal cell growth. In PCF, REH2C controls editing fidelity in RESC-bound mRNAs. Here, we report that KREH2, a REH2C-associated helicase, developmentally controls programmed non-canonical editing, including an abundant 3' element in ATPase subunit 6 (A6) mRNA. The 3' element sequence is directed by a proposed novel regulatory gRNA. In PCF, KREH2 RNAi-knockdown up-regulates the 3' element, which establishes a stable structure hindering element removal by canonical initiator-gRNA-directed editing. In BSF, KREH2-knockdown does not up-regulate the 3' element but reduces its high abundance. Thus, KREH2 differentially controls extensive non-canonical editing and associated RNA structure via a novel regulatory gRNA, potentially hijacking factors as a 'molecular sponge'. Furthermore, this gRNA is bifunctional, serving in canonical CR4 mRNA editing whilst installing a structural element in A6 mRNA.

Recurrent emergence of Klebsiella pneumoniae carbapenem resistance mediated by an inhibitory ompK36 mRNA secondary structure.

Outer membrane porins in Gram-negative bacteria facilitate antibiotic influx. In Klebsiella pneumoniae, modifications in the porin OmpK36 are implicated in increasing resistance to carbapenems. An analysis of large K. pneumoniae genome collections, encompassing major healthcare-associated clones, revealed the recurrent emergence of a synonymous cytosine-to-thymine transition at position 25 (25c > t) in ompK36. We show that the 25c > t transition increases carbapenem resistance through depletion of OmpK36 from the outer membrane. The mutation attenuates K. pneumoniae in a murine pneumonia model, which accounts for its limited clonal expansion observed by phylogenetic analysis. However, in the context of carbapenem treatment, the 25c > t transition tips the balance toward treatment failure, thus accounting for its recurrent emergence. Mechanistically, the 25c > t transition mediates an intramolecular messenger RNA (mRNA) interaction between a uracil encoded by 25t and the first adenine within the Shine-Dalgarno sequence. This specific interaction leads to the formation of an RNA stem structure, which obscures the ribosomal binding site thus disrupting translation. While mutations reducing OmpK36 expression via transcriptional silencing are known, we uniquely demonstrate the repeated selection of a synonymous ompK36 mutation mediating translational suppression in response to antibiotic pressure.

Web-based platform for analysis of RNA folding from high throughput chemical probing data.

RNA structures play critical roles in regulating gene expression across all domains of life and viruses. Chemical probing methods coupled with massively parallel sequencing have revolutionized the RNA structure field by enabling the assessment of many structures in their native, physiological context. Previously, we developed Dimethyl-Sulfate-based Mutational Profiling and Sequencing (DMS-MaPseq), which uses DMS to label the Watson-Crick face of open and accessible adenine and cytosine bases in the RNA. We used this approach to determine the genome-wide structures of HIV-1 and SARS-CoV-2 in infected cells, which permitted uncovering new biology and identifying therapeutic targets. Due to the simplicity and ease of the experimental procedure, DMS-MaPseq has been adopted by labs worldwide. However, bioinformatic analysis remains a substantial hurdle for labs that often lack the necessary infrastructure and computational expertise. Here we present a modern web-based interface that automates the analysis of chemical probing profiles from raw sequencing files (http://rnadreem.org). The availability of a simple web-based platform for DMS-MaPseq analysis will dramatically expand studies of RNA structure and aid in the design of structure-based therapeutics.

Sigma factor dependent translational activation in Bacillus subtilis.

Sigma factors are an important class of bacterial transcription factors that lend specificity to RNA polymerases by binding to distinct promoter elements for genes in their regulons. Here we show that activation of the general stress sigma factor, σB, in Bacillus subtilis paradoxically leads to dramatic induction of translation for a subset of its regulon genes. These genes are translationally repressed when transcribed by the housekeeping sigma factor, σA, owing to extended RNA secondary structures as determined in vivo using DMS-MaPseq. Transcription from σB-dependent promoters ablates the secondary structures and activates translation, leading to dual induction. Translation efficiencies between σB- and σA-dependent RNA isoforms can vary by up to 100-fold, which in multiple cases exceeds the magnitude of transcriptional induction. These results highlight the role of long-range RNA folding in modulating translation and demonstrate that a transcription factor can regulate protein synthesis beyond its effects on transcript levels.

Viral RNA structure analysis using DMS-MaPseq.

RNA structure is critically important to RNA viruses in every part of the replication cycle. RNA structure is also utilized by DNA viruses in order to regulate gene expression and interact with host factors. Advances in next-generation sequencing have greatly enhanced the utility of chemical probing in order to analyze RNA structure. This review will cover some recent viral RNA structural studies using chemical probing and next-generation sequencing as well as the advantages of dimethyl sulfate (DMS)-mutational profiling and sequencing (MaPseq). DMS-MaPseq is a robust assay that can easily modify RNA in vitro, in cell and in virion. A detailed protocol for whole-genome DMS-MaPseq from cells transfected with HIV-1 and the structure of TAR as determined by DMS-MaPseq is presented. DMS-MaPseq has the ability to answer a variety of integral questions about viral RNA, including how they change in different environments and when interacting with different host factors.

DMS-MaPseq for genome-wide or targeted RNA structure probing in vivo.

Coupling of structure-specific in vivo chemical modification to next-generation sequencing is transforming RNA secondary structure studies in living cells. The dominant strategy for detecting in vivo chemical modifications uses reverse transcriptase truncation products, which introduce biases and necessitate population-average assessments of RNA structure. Here we present dimethyl sulfate (DMS) mutational profiling with sequencing (DMS-MaPseq), which encodes DMS modifications as mismatches using a thermostable group II intron reverse transcriptase. DMS-MaPseq yields a high signal-to-noise ratio, can report multiple structural features per molecule, and allows both genome-wide studies and focused in vivo investigations of even low-abundance RNAs. We apply DMS-MaPseq for the first analysis of RNA structure within an animal tissue and to identify a functional structure involved in noncanonical translation initiation. Additionally, we use DMS-MaPseq to compare the in vivo structure of pre-mRNAs with their mature isoforms. These applications illustrate DMS-MaPseq's capacity to dramatically expand in vivo analysis of RNA structure.

Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo.

RNA has a dual role as an informational molecule and a direct effector of biological tasks. The latter function is enabled by RNA's ability to adopt complex secondary and tertiary folds and thus has motivated extensive computational and experimental efforts for determining RNA structures. Existing approaches for evaluating RNA structure have been largely limited to in vitro systems, yet the thermodynamic forces which drive RNA folding in vitro may not be sufficient to predict stable RNA structures in vivo. Indeed, the presence of RNA-binding proteins and ATP-dependent helicases can influence which structures are present inside cells. Here we present an approach for globally monitoring RNA structure in native conditions in vivo with single-nucleotide precision. This method is based on in vivo modification with dimethyl sulphate (DMS), which reacts with unpaired adenine and cytosine residues, followed by deep sequencing to monitor modifications. Our data from yeast and mammalian cells are in excellent agreement with known messenger RNA structures and with the high-resolution crystal structure of the Saccharomyces cerevisiae ribosome. Comparison between in vivo and in vitro data reveals that in rapidly dividing cells there are vastly fewer structured mRNA regions in vivo than in vitro. Even thermostable RNA structures are often denatured in cells, highlighting the importance of cellular processes in regulating RNA structure. Indeed, analysis of mRNA structure under ATP-depleted conditions in yeast shows that energy-dependent processes strongly contribute to the predominantly unfolded state of mRNAs inside cells. Our studies broadly enable the functional analysis of physiological RNA structures and reveal that, in contrast to the Anfinsen view of protein folding whereby the structure formed is the most thermodynamically favourable, thermodynamics have an incomplete role in determining mRNA structure in vivo.

Causal signals between codon bias, mRNA structure, and the efficiency of translation and elongation.

Ribosome profiling data report on the distribution of translating ribosomes, at steady-state, with codon-level resolution. We present a robust method to extract codon translation rates and protein synthesis rates from these data, and identify causal features associated with elongation and translation efficiency in physiological conditions in yeast. We show that neither elongation rate nor translational efficiency is improved by experimental manipulation of the abundance or body sequence of the rare AGG tRNA. Deletion of three of the four copies of the heavily used ACA tRNA shows a modest efficiency decrease that could be explained by other rate-reducing signals at gene start. This suggests that correlation between codon bias and efficiency arises as selection for codons to utilize translation machinery efficiently in highly translated genes. We also show a correlation between efficiency and RNA structure calculated both computationally and from recent structure probing data, as well as the Kozak initiation motif, which may comprise a mechanism to regulate initiation.

Untangling the pseudoknots of SARS-CoV-2: Insights into structural heterogeneity and plasticity.

Since the onset of the COVID-19 pandemic, one productive area of research has focused on the intricate two- and three-dimensional structures taken on by SARS-CoV-2's RNA genome. These structures control essential viral processes, making them tempting targets for therapeutic intervention. This review focuses on two such structured regions, the frameshift stimulation element (FSE), which controls the translation of viral protein, and the 3' untranslated region (3' UTR), which is thought to regulate genome replication. For the FSE, we discuss its canonical pseudoknot's threaded and unthreaded topologies, as well as the diversity of competing two-dimensional structures formed by local and long-distance base pairing. For the 3' UTR, we review the evidence both for and against the formation of its replication-enabling pseudoknot.

The Folding of Germ Granule mRNAs Controls Intermolecular Base Pairing in Germ Granules and Maintains Normal Fly Development.

Drosophila germ granules enrich mRNAs critical for fly development. Within germ granules, mRNAs form multi-transcript clusters marked by increased mRNA concentration, creating an elevated potential for intermolecular base pairing. However, the type and abundance of intermolecular base pairing in mRNA clusters is poorly characterized. Using single-molecule super-resolution microscopy, chemical probing for base accessibility, phase separation assays, and simulations, we demonstrated that mRNAs remain well-folded upon localization to germ granules. While most base pairing is intramolecular, mRNAs still display the ability for intermolecular base pairing, facilitating clustering without high sequence complementarity or significant melting of secondary structure. This base pairing among mRNAs is driven by scattered and discontinuous stretches of bases appearing on the surface of folded RNAs, providing multivalency to clustering but exhibits low probability for sustained interactions. Notably, engineered germ granule mRNAs with exposed GC-rich complementary sequences (CSs) presented within stable stem loops induce sustained base pairing in vitro and enhanced intermolecular interactions in vivo. However, the presence of these stem loops alone disrupts fly development, and the addition of GC-rich CSs exacerbates this phenotype. Although germ granule mRNAs contain numerous GC-rich CSs capable of stable intermolecular base pairing, they are primarily embedded by RNA folding. This study emphasizes the role of RNA folding in controlling the type and abundance of intermolecular base pairing, thereby preserving the functional integrity of mRNAs within the germ granules.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression.

The human mitochondrial genome encodes crucial oxidative phosphorylation system proteins, pivotal for aerobic energy transduction. They are translated from nine monocistronic and two bicistronic transcripts whose native structures remain unexplored, posing a gap in understanding mitochondrial gene expression. In this work, we devised the mitochondrial dimethyl sulfate mutational profiling with sequencing (mitoDMS-MaPseq) method and applied detection of RNA folding ensembles using expectation-maximization (DREEM) clustering to unravel the native mitochondrial messenger RNA (mt-mRNA) structurome in wild-type (WT) and leucine-rich pentatricopeptide repeat-containing protein (LRPPRC)-deficient cells. Our findings elucidate LRPPRC's role as a holdase contributing to maintaining mt-mRNA folding and efficient translation. mt-mRNA structural insights in WT mitochondria, coupled with metabolic labeling, unveil potential mRNA-programmed translational pausing and a distinct programmed ribosomal frameshifting mechanism. Our data define a critical layer of mitochondrial gene expression regulation. These mt-mRNA folding maps provide a reference for studying mt-mRNA structures in diverse physiological and pathological contexts.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA.

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA.

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.