Comprehensive in vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms.

Severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2) is the positive-sense RNA virus that causes coronavirus disease 2019 (COVID-19). The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form RNA structures, yet as much as 97% of its 30 kilobases have not been structurally explored. Here, we apply a novel long amplicon strategy to determine the secondary structure of the SARS-CoV-2 RNA genome at single-nucleotide resolution in infected cells. Our in-depth structural analysis reveals networks of well-folded RNA structures throughout Orf1ab and reveals aspects of SARS-CoV-2 genome architecture that distinguish it from other RNA viruses. Evolutionary analysis shows that several features of the SARS-CoV-2 genomic structure are conserved across ß-coronaviruses, and we pinpoint regions of well-folded RNA structure that merit downstream functional analysis. The native, secondary structure of SARS-CoV-2 presented here is a roadmap that will facilitate focused studies on the viral life cycle, facilitate primer design, and guide the identification of RNA drug targets against COVID-19.

The West Nile virus genome harbors essential riboregulatory elements with conserved and host-specific functional roles.

West Nile virus (WNV) is an arthropod-borne, positive-sense RNA virus that poses an increasing global threat due to warming climates and lack of effective therapeutics. Like other enzootic viruses, little is known about how host context affects the structure of the full-length RNA genome. Here, we report a complete secondary structure of the entire WNV genome within infected mammalian and arthropod cell lines. Our analysis affords structural insights into multiple, conserved aspects of flaviviral biology. We show that the WNV genome folds with minimal host dependence, and we prioritize well-folded regions for functional validation using structural homology between hosts as a guide. Using structure-disrupting, antisense locked nucleic acids, we then demonstrate that the WNV genome contains riboregulatory structures with conserved and host-specific functional roles. These results reveal promising RNA drug targets within flaviviral genomes, and they highlight the therapeutic potential of ASO-LNAs as both WNV-specific and pan-flaviviral therapeutic agents.

Abolished frameshifting for predicted structure-stabilizing SARS-CoV-2 mutants: Implications to alternative conformations and their statistical structural analyses.

The SARS-CoV-2 frameshifting element (FSE) has been intensely studied and explored as a therapeutic target for coronavirus diseases including COVID-19. Besides the intriguing virology, this small RNA is known to adopt many length-dependent conformations, as verified by multiple experimental and computational approaches. However, the role these alternative conformations play in the frameshifting mechanism and how to quantify this structural abundance has been an ongoing challenge. Here, we show by DMS and dual-luciferase functional assays that previously predicted FSE mutants (using the RAG graph theory approach) suppress structural transitions and abolish frameshifting. Furthermore, correlated mutation analysis of DMS data by three programs (DREEM, DRACO, and DANCE-MaP) reveals important differences in their estimation of specific RNA conformations, suggesting caution in the interpretation of such complex conformational landscapes. Overall, the abolished frameshifting in three different mutants confirms that all alternative conformations play a role in the pathways of ribosomal transition.

Structure-first identification of RNA elements that regulate dengue virus genome architecture and replication.

The genomes of RNA viruses encode the information required for replication in host cells both in their linear sequence and in complex higher-order structures. A subset of these RNA genome structures show clear sequence conservation, and have been extensively described for well-characterized viruses. However, the extent to which viral RNA genomes contain functional structural elements-unable to be detected by sequence alone-that nonetheless are critical to viral fitness is largely unknown. Here, we devise a structure-first experimental strategy and use it to identify 22 structure-similar motifs across the coding sequences of the RNA genomes for the four dengue virus serotypes. At least 10 of these motifs modulate viral fitness, revealing a significant unnoticed extent of RNA structure-mediated regulation within viral coding sequences. These viral RNA structures promote a compact global genome architecture, interact with proteins, and regulate the viral replication cycle. These motifs are also thus constrained at the levels of both RNA structure and protein sequence and are potential resistance-refractory targets for antivirals and live-attenuated vaccines. Structure-first identification of conserved RNA structure enables efficient discovery of pervasive RNA-mediated regulation in viral genomes and, likely, other cellular RNAs.

SRSF3 shapes the structure of miR-17-92 cluster RNA and promotes selective processing of miR-17 and miR-20a.

MicroRNA (miRNA) biogenesis is tightly regulated to maintain distinct miRNA expression patterns. Almost half of mammalian miRNAs are generated from miRNA clusters, but this process is not well understood. We show here that Serine-arginine rich splicing factor 3 (SRSF3) controls the processing of miR-17-92 cluster miRNAs in pluripotent and cancer cells. SRSF3 binding to multiple CNNC motifs downstream of Drosha cleavage sites within miR-17-92 is required for the efficient processing of the cluster. SRSF3 depletion specifically compromises the processing of two paralog miRNAs, miR-17 and miR-20a. In addition to SRSF3 binding to the CNNC sites, the SRSF3 RS-domain is essential for miR-17-92 processing. SHAPE-MaP probing demonstrates that SRSF3 binding disrupts local and distant base pairing, resulting in global changes in miR-17-92 RNA structure. Our data suggest a model where SRSF3 binding, and potentially its RS-domain interactions, may facilitate an RNA structure that promotes miR-17-92 processing. SRSF3-mediated increase in miR-17/20a levels inhibits the cell cycle inhibitor p21, promoting self-renewal in normal and cancer cells. The SRSF3-miR-17-92-p21 pathway operates in colorectal cancer, linking SRSF3-mediated pri-miRNA processing and cancer pathogenesis.

SHAPE-enabled fragment-based ligand discovery for RNA.

The transcriptome represents an attractive but underused set of targets for small-molecule ligands. Here, we devise a technology that leverages fragment-based screening and SHAPE-MaP RNA structure probing to discover small-molecule fragments that bind an RNA structure of interest. We identified fragments and cooperatively binding fragment pairs that bind to the thiamine pyrophosphate (TPP) riboswitch with millimolar to micromolar affinities. We then used structure-activity relationship information to efficiently design a linked-fragment ligand, with no resemblance to the native ligand, with high ligand efficiency and druglikeness, that binds to the TPP thiM riboswitch with high nanomolar affinity and that modulates RNA conformation during cotranscriptional folding. Principles from this work are broadly applicable, leveraging cooperativity and multisite binding, for developing high-quality ligands for diverse RNA targets.

The m6A landscape of polyadenylated nuclear (PAN) RNA and its related methylome in the context of KSHV replication.

Polyadenylated nuclear (PAN) RNA is a long noncoding transcript involved in Kaposi's sarcoma-associated herpesvirus (KSHV) lytic reactivation and regulation of cellular and viral gene expression. We have previously shown that PAN RNA has dynamic secondary structure and protein binding profiles that can be influenced by epitranscriptomic modifications. N6-methyladenosine (m6A) is one of the most abundant chemical signatures found in viral RNA genomes and virus-encoded RNAs. Here, we combined antibody-independent next-generation mapping with direct RNA sequencing to address the epitranscriptomic status of PAN RNA in KSHV infected cells. We showed that PAN m6A status is dynamic, reaching the highest number of modifications at the late lytic stages of KSHV infection. Using a newly developed method, termed selenium-modified deoxythymidine triphosphate (SedTTP)-reverse transcription (RT) and ligation assisted PCR analysis of m6A (SLAP), we gained insight into the fraction of modification at identified sites. By applying comprehensive proteomic approaches, we identified writers and erasers that regulate the m6A status of PAN, and readers that can convey PAN m6A phenotypic effects. We verified the temporal and spatial subcellular availability of the methylome components for PAN modification by performing confocal microscopy analysis. Additionally, the RNA biochemical probing (SHAPE-MaP) outlined local and global structural alterations invoked by m6A in the context of full-length PAN RNA. This work represents the first comprehensive overview of the dynamic interplay that takes place between the cellular epitranscriptomic machinery and a specific viral RNA in the context of KSHV infected cells.

The Evaluation of SHAPE-MaP RNA Structure Probing Protocols Reveals a Novel Role of Mn2+ in the Detection of 2'-OH Adducts.

Chemical probing, for decades, has been one of the most popular tools for studying the secondary structure of RNA molecules. Recently, protocols for simultaneous analysis of multiple RNAs have been developed, enabling in vivo transcriptome-wide interrogation of the RNA structure dynamics. One of the most popular methods is the selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP). In this study, we describe the evaluation of this protocol by addressing the influence of the reverse transcription enzymes, buffer conditions, and chemical probes on the properties of the cDNA library and the quality of mutational profiling-derived structural signals. Our results reveal a SuperScript IV (SSIV) reverse transcriptase as a more efficient enzyme for mutational profiling of SHAPE adducts and shed new light on the role of Mn2+ cations in the modulation of SSIV readthrough efficiency.

RNA structure probing to characterize RNA-protein interactions on a low abundance pre-mRNA in living cells.

In vivo RNA structure analysis has become a powerful tool in molecular biology, largely due to the coupling of an increasingly diverse set of chemical approaches with high-throughput sequencing. This has resulted in a transition from single target to transcriptome-wide approaches. However, these methods require sequencing depths that preclude studying low abundance targets, which are not sufficiently captured in transcriptome-wide approaches. Here we present a ligation-free method to enrich for low abundance RNA sequences, which improves the diversity of molecules analyzed and results in improved analysis. In addition, this method is compatible with any choice of chemical adduct or read-out approach. We utilized this approach to study an autoregulated event in the pre-mRNA of the splicing factor, muscleblind-like splicing regulator 1 (MBNL1).

Insights into the secondary and tertiary structure of the Bovine Viral Diarrhea Virus Internal Ribosome Entry Site.

The internal ribosome entry site (IRES) RNA of bovine viral diarrhoea virus (BVDV), an economically significant Pestivirus, is required for the cap-independent translation of viral genomic RNA. Thus, it is essential for viral replication and pathogenesis. We applied a combination of high-throughput biochemical RNA structure probing (SHAPE-MaP) and in silico modelling approaches to gain insight into the secondary and tertiary structures of BVDV IRES RNA. Our study demonstrated that BVDV IRES RNA in solution forms a modular architecture composed of three distinct structural domains (I-III). Two regions within domain III are represented in tertiary interactions to form an H-type pseudoknot. Computational modelling of the pseudoknot motif provided a fine-grained picture of the tertiary structure and local arrangement of helices in the BVDV IRES. Furthermore, comparative genomics and consensus structure predictions revealed that the pseudoknot is evolutionarily conserved among many Pestivirus species. These studies provide detailed insight into the structural arrangement of BVDV IRES RNA H-type pseudoknot and encompassing motifs that likely contribute to the optimal functionality of viral cap-independent translation element.

Mapping circular RNA structures in living cells by SHAPE-MaP.

Circular RNAs are produced from back-splicing of exons of precursor mRNAs (pre-mRNAs). The sequences of exons in circular RNAs are identical to their linear cognate mRNAs, but the circular format may confer constraints on their folding and conformation, leading to potentially different functions from their linear RNA cognates. Here, we describe experimental and computational steps that optimize the selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) to probe circular RNA secondary structure at single-nucleotide resolution in living cells.

Using SHAPE-MaP To Model RNA Secondary Structure and Identify 3'UTR Variation in Chikungunya Virus.

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus associated with debilitating arthralgia in humans. RNA secondary structure in the viral genome plays an important role in the lifecycle of alphaviruses; however, the specific role of RNA structure in regulating CHIKV replication is poorly understood. Our previous studies found little conservation in RNA secondary structure between alphaviruses, and this structural divergence creates unique functional structures in specific alphavirus genomes. Therefore, to understand the impact of RNA structure on CHIKV biology, we used SHAPE-MaP to inform the modeling of RNA secondary structure throughout the genome of a CHIKV isolate from the 2013 Caribbean outbreak. We then analyzed regions of the genome with high levels of structural specificity to identify potentially functional RNA secondary structures and identified 23 regions within the CHIKV genome with higher than average structural stability, including four previously identified, functionally important CHIKV RNA structures. We also analyzed the RNA flexibility and secondary structures of multiple 3'UTR variants of CHIKV that are known to affect virus replication in mosquito cells. This analysis found several novel RNA structures within these 3'UTR variants. A duplication in the 3'UTR that enhances viral replication in mosquito cells led to an overall increase in the amount of unstructured RNA in the 3'UTR. This analysis demonstrates that the CHIKV genome contains a number of unique, specific RNA secondary structures and provides a strategy for testing these secondary structures for functional importance in CHIKV replication and pathogenesis.IMPORTANCE Chikungunya virus (CHIKV) is a mosquito-borne RNA virus that causes febrile illness and debilitating arthralgia in humans. CHIKV causes explosive outbreaks but there are no approved therapies to treat or prevent CHIKV infection. The CHIKV genome contains functional RNA secondary structures that are essential for proper virus replication. Since RNA secondary structures have only been defined for a small portion of the CHIKV genome, we used a chemical probing method to define the RNA secondary structures of CHIKV genomic RNA. We identified 23 highly specific structured regions of the genome, and confirmed the functional importance of one structure using mutagenesis. Furthermore, we defined the RNA secondary structure of three CHIKV 3'UTR variants that differ in their ability to replicate in mosquito cells. Our study highlights the complexity of the CHIKV genome and describes new systems for designing compensatory mutations to test the functional relevance of viral RNA secondary structures.

Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure.

The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1 Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.

Predicting RNA SHAPE scores with deep learning.

Secondary structure prediction approaches rely typically on models of equilibrium free energies that are themselves based on in vitro physical chemistry. Recent transcriptome-wide experiments of in vivo RNA structure based on SHAPE-MaP experiments provide important information that may make it possible to extend current in vitro-based RNA folding models in order to improve the accuracy of computational RNA folding simulations with respect to the experimentally measured in vivo RNA secondary structure. Here we present a machine learning approach that utilizes RNA secondary structure prediction results and nucleotide sequence in order to predict in vivo SHAPE scores. We show that this approach has a higher Pearson correlation coefficient with experimental SHAPE scores than thermodynamic folding. This could be an important step towards augmenting experimental results with computational predictions and help with RNA secondary structure predictions that inherently take in-vivo folding properties into account.

Using SHAPE-MaP to probe small molecule-RNA interactions.

RNA is a regulator and catalyst of many cellular processes. Efforts to therapeutically harness RNA began with the discovery of myriad coding and non-coding RNAs and their versatile modes of action. However, due to its dynamic structure and the polar and repetitive nature of its surface, RNA presents a challenging target for drug design. For an RNA to be druggable, it must contain a motif that assumes a nearly fixed and unique conformation that a small molecule can recognize and bind consistently and with high affinity. Hence, reliable methods for determining the secondary and tertiary structures of RNA, and even the features and occupancy of potential drug binding sites are of utmost importance for the effective design of RNA-based therapeutics. Selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) has emerged as such a method, by which RNA secondary structure can be probed at single-nucleotide resolution, under a variety of conditions, and in the presence of RNA-specific small-molecule ligands. In this review, we describe an in-depth protocol for using SHAPE-MaP to characterize RNA-small molecule interactions in cell culture (in cellulo). This method can be applied to transcripts of any size or abundance, and to determine the sites and affinities of small molecule binding, making it an essential and versatile tool for drug discovery.

Visualization of RNA structure models within the Integrative Genomics Viewer.

Analyses of the interrelationships between RNA structure and function are increasingly important components of genomic studies. The SHAPE-MaP strategy enables accurate RNA structure probing and realistic structure modeling of kilobase-length noncoding RNAs and mRNAs. Existing tools for visualizing RNA structure models are not suitable for efficient analysis of long, structurally heterogeneous RNAs. In addition, structure models are often advantageously interpreted in the context of other experimental data and gene annotation information, for which few tools currently exist. We have developed a module within the widely used and well supported open-source Integrative Genomics Viewer (IGV) that allows visualization of SHAPE and other chemical probing data, including raw reactivities, data-driven structural entropies, and data-constrained base-pair secondary structure models, in context with linear genomic data tracks. We illustrate the usefulness of visualizing RNA structure in the IGV by exploring structure models for a large viral RNA genome, comparing bacterial mRNA structure in cells with its structure under cell- and protein-free conditions, and comparing a noncoding RNA structure modeled using SHAPE data with a base-pairing model inferred through sequence covariation analysis.

SHAPE reveals transcript-wide interactions, complex structural domains, and protein interactions across the Xist lncRNA in living cells.

The 18-kb Xist long noncoding RNA (lncRNA) is essential for X-chromosome inactivation during female eutherian mammalian development. Global structural architecture, cell-induced conformational changes, and protein-RNA interactions within Xist are poorly understood. We used selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) to examine these features of Xist at single-nucleotide resolution both in living cells and ex vivo. The Xist RNA forms complex well-defined secondary structure domains and the cellular environment strongly modulates the RNA structure, via motifs spanning one-half of all Xist nucleotides. The Xist RNA structure modulates protein interactions in cells via multiple mechanisms. For example, repeat-containing elements adopt accessible and dynamic structures that function as landing pads for protein cofactors. Structured RNA motifs create interaction domains for specific proteins and also sequester other motifs, such that only a subset of potential binding sites forms stable interactions. This work creates a broad quantitative framework for understanding structure-function interrelationships for Xist and other lncRNAs in cells.

The "Topological Train Ride" of a viral long non-coding RNA.

As the notion of small molecule targeting of regulatory viral and cellular RNAs gathers momentum, understanding their structure, and variations thereof, in the appropriate biological context will play a critical role. This is especially true of the ∼1100-nt polyadenylated nuclear (PAN) long non-coding (lnc) RNA of Kaposi's sarcoma herpesvirus (KSHV), whose interaction with viral and cellular proteins is central to lytic infection. Nuclear accumulation of PAN RNA is mediated via a unique triple helical structure at its 3' terminus (within the expression and nuclear retention element, or ENE) which protects it from deadenylation-dependent decay. Additionally, significant levels of PAN RNA have been reported in both the cytoplasm of KSHV-infected cells and in budding virions, leading us to consider which viral and host proteins might associate with, or dissociate from, this lncRNA during its "journey" through the cell. By combining the power of SHAPE-mutational profiling (SHAPE-MaP) with large scale virus culture facilities of the National Cancer Institute, Frederick MD, Sztuba-Solinska et al. have provide the first detailed description of KSHV PAN nucleoprotein complexes in multiple biological contexts, complementing this by mapping sites of recombinant KSHV proteins on an in vitro-synthesized, polyadenylated counterpart.

A 44-Nucleotide Region in the Chikungunya Virus 3' UTR Dictates Viral Fitness in Disparate Host Cells.

We previously reported that deletion of a 44-nucleotide element in the 3' untranslated region (UTR) of the Chikungunya virus (CHIKV) genome enhances the virulence of CHIKV infection in mice. Here, we find that while this 44-nucleotide deletion enhances CHIKV fitness in murine embryonic fibroblasts in a manner independent of the type I interferon response, the same mutation decreases viral fitness in C6/36 mosquito cells. Further, the fitness advantage conferred by the UTR deletion in mammalian cells is maintained in vivo in a mouse model of CHIKV dissemination. Finally, SHAPE-MaP analysis of the CHIKV 3' UTR revealed this 44-nucleotide element forms a distinctive two-stem-loop structure that is ablated in the mutant 3' UTR without altering additional 3' UTR RNA secondary structures.

Analysis of Telomerase RNA Structure in Physcomitrium patens Indicates Functionally Relevant Transitions Between OPEN and CLOSED Conformations.

Telomerase RNA (TR) conformation determines its function as a template for telomere synthesis and as a scaffold for the assembly of the telomerase nucleoprotein complex. Experimental analyses of TR secondary structure using DMS-Map Seq and SHAPE-Map Seq techniques show its CLOSED conformation as the consensus structure where the template region cannot perform its function. Our data show that the apparent discrepancy between experimental results and predicted TR functional conformation, mostly ignored in published studies, can be explained using data analysis based on single-molecule structure prediction from individual sequencing reads by the recently established DaVinci method. This method results in several clusters of secondary structures reflecting the structural dynamics of TR, possibly related to its multiple functional states. Interestingly, the presumed active (OPEN) conformation of TR corresponds to a minor fraction of TR under in vivo conditions. Therefore, structural polymorphism and dynamic TR transitions between CLOSED and OPEN conformations may be involved in telomerase activity regulation as a switch that functions independently of total TR transcript levels.