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.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression.

The mammalian mitochondrial genome encodes thirteen oxidative phosphorylation system proteins, crucial in aerobic energy transduction. These proteins are translated from 9 monocistronic and 2 bicistronic transcripts, whose native structures remain unexplored, leaving fundamental molecular determinants of mitochondrial gene expression unknown. To address this gap, we developed a mitoDMS-MaPseq approach and used DREEM clustering to resolve the native human mitochondrial mt-mRNA structurome. We gained insights into mt-mRNA biology and translation regulatory mechanisms, including a unique programmed ribosomal frameshifting for the ATP8/ATP6 transcript. Furthermore, absence of the mt-mRNA maintenance factor LRPPRC led to a mitochondrial transcriptome structured differently, with specific mRNA regions exhibiting increased or decreased structuredness. This highlights the role of LRPPRC in maintaining mRNA folding to promote mt-mRNA stabilization and efficient translation. In conclusion, our mt-mRNA folding maps reveal novel mitochondrial gene expression mechanisms, serving as a detailed reference and tool for studying them in different physiological and pathological contexts.

Dissecting telomerase RNA structural heterogeneity in living human cells with DMS-MaPseq.

Telomerase is a specialized reverse transcriptase that uses an intrinsic RNA subunit as the template for telomeric DNA synthesis. Biogenesis of human telomerase requires its RNA subunit (hTR) to fold into a multi-domain architecture that includes the template-containing pseudoknot (t/PK) and the three-way junction (CR4/5). These two hTR domains bind the telomerase reverse transcriptase (hTERT) protein and are thus essential for telomerase catalytic activity. Here, we probe the structure of hTR in living cells using dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) and ensemble deconvolution analysis. Unexpectedly, approximately 15% of the steady state population of hTR has a CR4/5 conformation lacking features required for hTERT binding. Mutagenesis demonstrates that stabilization of the alternative CR4/5 conformation is detrimental to telomerase assembly and activity. We propose that this misfolded portion of the cellular hTR pool is either slowly refolded or degraded. Thus, kinetic traps for RNA folding that have been so well-studied in vitro may also present barriers for assembly of ribonucleoprotein complexes in vivo.

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.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts.

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

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.

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.

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

As an essential post-transcriptional regulator of gene expression, microRNA (miR) levels must be strictly maintained. The biogenesis of many, but not all, miRs 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 pre-miR-31 can provide a mechanism by which its biogenesis is regulated. We determined the solution structure of the precursor element of miR-31 (pre-miR-31) to investigate the role of distinct structural elements in regulating Dicer processing. We found that the presence or absence of mismatches within the helical stem do not strongly influence Dicer processing of the pre-miR. However, both the apical loop size and structure at the Dicing site are key elements for discrimination by Dicer. 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 influence both Dicer binding and processing. Our results enrich our understanding of the active role that RNA structure plays in regulating Dicer processing which has direct implications for control of gene expression.

3D RNA-scaffolded wireframe origami.

Hybrid RNA:DNA origami, in which a long RNA scaffold strand folds into a target nanostructure via thermal annealing with complementary DNA oligos, has only been explored to a limited extent despite its unique potential for biomedical delivery of mRNA, tertiary structure characterization of long RNAs, and fabrication of artificial ribozymes. Here, we investigate design principles of three-dimensional wireframe RNA-scaffolded origami rendered as polyhedra composed of dual-duplex edges. We computationally design, fabricate, and characterize tetrahedra folded from an EGFP-encoding messenger RNA and de Bruijn sequences, an octahedron folded with M13 transcript RNA, and an octahedron and pentagonal bipyramids folded with 23S ribosomal RNA, demonstrating the ability to make diverse polyhedral shapes with distinct structural and functional RNA scaffolds. We characterize secondary and tertiary structures using dimethyl sulfate mutational profiling and cryo-electron microscopy, revealing insight into both global and local, base-level structures of origami. Our top-down sequence design strategy enables the use of long RNAs as functional scaffolds for complex wireframe origami.

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.

An intranasal ASO therapeutic targeting SARS-CoV-2.

The COVID-19 pandemic is exacting an increasing toll worldwide, with new SARS-CoV-2 variants emerging that exhibit higher infectivity rates and that may partially evade vaccine and antibody immunity. Rapid deployment of non-invasive therapeutic avenues capable of preventing infection by all SARS-CoV-2 variants could complement current vaccination efforts and help turn the tide on the COVID-19 pandemic. Here, we describe a novel therapeutic strategy targeting the SARS-CoV-2 RNA using locked nucleic acid antisense oligonucleotides (LNA ASOs). We identify an LNA ASO binding to the 5' leader sequence of SARS-CoV-2 that disrupts a highly conserved stem-loop structure with nanomolar efficacy in preventing viral replication in human cells. Daily intranasal administration of this LNA ASO in the COVID-19 mouse model potently suppresses viral replication (>80-fold) in the lungs of infected mice. We find that the LNA ASO is efficacious in countering all SARS-CoV-2 "variants of concern" tested both in vitro and in vivo. Hence, inhaled LNA ASOs targeting SARS-CoV-2 represents a promising therapeutic approach to reduce or prevent transmission and decrease severity of COVID-19 in infected individuals. LNA ASOs are chemically stable and can be flexibly modified to target different viral RNA sequences and could be stockpiled for future coronavirus pandemics.

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.

Secondary structural ensembles of the SARS-CoV-2 RNA genome in infected cells.

SARS-CoV-2 is a betacoronavirus with a single-stranded, positive-sense, 30-kilobase RNA genome responsible for the ongoing COVID-19 pandemic. Although population average structure models of the genome were recently reported, there is little experimental data on native structural ensembles, and most structures lack functional characterization. Here we report secondary structure heterogeneity of the entire SARS-CoV-2 genome in two lines of infected cells at single nucleotide resolution. Our results reveal alternative RNA conformations across the genome and at the critical frameshifting stimulation element (FSE) that are drastically different from prevailing population average models. Importantly, we find that this structural ensemble promotes frameshifting rates much higher than the canonical minimal FSE and similar to ribosome profiling studies. Our results highlight the value of studying RNA in its full length and cellular context. The genomic structures detailed here lay groundwork for coronavirus RNA biology and will guide the design of SARS-CoV-2 RNA-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.

HnRNP A1/A2 Proteins Assemble onto 7SK snRNA via Context Dependent Interactions.

7SK small nuclear RNA (snRNA) is an abundant and ubiquitously expressed noncoding RNA that functions to modulate the activity of RNA Polymerase II (RNAPII) in part by stabilizing distinct pools of 7SK-protein complexes. Prevailing models suggest that the secondary structure of 7SK is dynamically remodeled within its alternative RNA-protein pools such that its architecture differentially regulates the exchange of cognate binding partners. The nuclear hnRNP A1/A2 proteins influence the biology of 7SK snRNA via processes that require an intact stem loop (SL) 3 domain; however, the molecular details by which hnRNPs assemble onto 7SK snRNA are yet to be described. Here, we have taken an integrated approach to present a detailed description of the 7SK-hnRNP A1 complex. We show that unbound 7SK snRNA adopts at least two major conformations in solution, with significant structural differences localizing to the SL2-3 linker and the base of SL3. Phylogenetic analysis indicates that this same region is the least genetically conserved feature of 7SK snRNA. By performing DMS modifications with the presence of excess protein, we reveal that hnRNP A1 binds with selectivity to SL3 through mechanisms that increase the flexibility of the RNA adjacent to putative binding sites. Calorimetric titrations further validate that hnRNP A1-SL3 assembly is complex with the affinity of discrete binding events modulated by the surrounding RNA structure. To interpret this context-dependent binding phenomenon, we determined a 3D model of SL3 to show that it folds to position minimal hnRNP A1/A2 binding sites (5'-Y/RAG-3') within different local environments. SL3-protein complexes resolved by SEC-MALS-SAXS confirm that up to four hnRNP A1 proteins bind along the entire surface of SL3 via interactions that preserve the overall structural integrity of this domain. In sum, the collective results presented here reveal a specific role for a folded SL3 domain to scaffold hnRNP A1/A2-7SK assembly via mechanisms modulated by the surrounding RNA structure.

DMS-MaPseq for Genome-Wide or Targeted RNA Structure Probing In Vitro and In Vivo.

DMS-MaPseq is a chemical probing method combined with high throughput sequencing used to study RNA structure. Here we present a flexible protocol for adherent and suspension mammalian cells to analyze RNA structure in vitro or in vivo. The protocol provides instruction on either a targeted sequencing of a lncRNA of interest or a transcriptome-wide approach that provides structural data on all expressed RNAs, including lncRNAs. This technique is particularly useful for comparing in vitro and in vivo structure of RNAs, determining how mutations and polymorphisms with phenotypic effects influence RNA structure and analyzing RNA structure across the entire transcriptome.

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.

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.

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.