RNA Structural Elements Research Articles

Functional RNA structures as targets in emerging arboviruses

Emerging arthropod-borne viruses, such as Zika virus (ZIKV) and Chikungunya virus (CHIKV), represent a significant and increasing threat to global health. Despite the expanding prevalence of their vectors, the Aedes sp.mosquitoes, and the potential for major epidemics, there are currently no specific antiviral compounds or vaccines available for either viral pathogen. The positive-strand genomes of ZIKV and CHIKV, members of the flavivirus and alphavirus genera respectively, contain functional, structured cis-acting RNA elements. We are investigating a range of approaches for targeting these RNA elements and analysing the effect on virus replication at different stages of their life cycle. Gaining high resolution structural data is essential prior to targeting RNA elements, and consequently, we mapped RNA structural elements within the ZIKV 5’ genome region using a combination of biochemical SHAPE probing, thermodynamic predictions and phylogenetic analysis. We are currently validating our structural data by analysis of mutant phenotypes in a reverse genetic system. Using antisense locked nucleic acid oligonucleotides (antisense-LNA), we demonstrated that functional RNA elements in CHIKV can be specifically targeted – inhibiting replication in both sub-genomic replicon and infectious virus systems. Surface plasmon resonance confirmed that the antisense-LNA binds to a specific stem-loop target with aKd of 310 nM and has an IC50 of 35 nM in the sub-genomic replicon system. In future work, we aim to investigate selection of RNA-aptamers against CHIKV and target ZIKV genomic stem-loops using antisense-LNAs.

Exoribonuclease-Resistant RNAs Exist within both Coding and Noncoding Subgenomic RNAs.

Many viruses produce protein-coding and noncoding subgenomic RNAs (sgRNAs) that are critical for infection. A recently discovered pathway for viral sgRNA production uses exoribonuclease-resistant RNAs (xrRNAs), discrete folded RNA elements that block the processive exoribonucleolytic degradation of RNA. xrRNAs are widespread in animal-infecting flaviviruses but had been found only in three members of the plant virus genus Dianthovirus Also, xrRNAs had been found only in the 3' untranslated regions (3'UTRs) of viral RNAs, where they produce noncoding sgRNAs. The degree to which xrRNA elements exist in other viruses, the conservation of their ring-like fold, and the ability of xrRNAs to operate in diverse contexts were unknown. Using computational tools and biochemical assays, we discovered xrRNA elements pervading two large families of plant-infecting RNA viruses, demonstrating their importance and widespread utility. Comparison of the sequences and functional requirements suggests that all adopt the characteristic ring-like fold. Unexpectedly, many of these newly discovered xrRNAs are located in intergenic regions rather than 3´UTRs, and some are associated with the 5' ends of subgenomic RNAs that encode viral proteins. This suggests that xrRNAs are involved in the production of both coding and noncoding subgenomic RNAs and can operate as part of broader mechanisms to regulate RNA levels and protein expression. These discoveries expand the potential roles for xrRNAs and suggest that xrRNAs may represent a more general strategy for RNA maturation and maintenance than previously known.IMPORTANCE During infection, viruses often produce subgenomic RNAs (sgRNAs) that either serve as the template for protein synthesis or act as "riboregulators" that interact with and influence the viral and cellular machinery. Recently, a mechanism for producing sgRNAs was found that depends on the presence of specifically structured RNA elements (xrRNAs). However, the degree to which this mechanism is used, where the elements are found, their structural diversity, and what types of sgRNAs are produced by this pathway were unclear. This article describes the discovery of these structured RNA elements in two large families of plant viruses and shows that they are used to produce both protein-coding sgRNAs and "riboregulatory" RNAs. These discoveries provide evidence that xrRNA-based RNA maturation pathways may be more widespread than previously anticipated and that they are involved in producing a variety of RNAs of diverse functions.

ScanFold: an approach for genome-wide discovery of local RNA structural elements-applications to Zika virus and HIV.

In addition to encoding RNA primary structures, genomes also encode RNA secondary and tertiary structures that play roles in gene regulation and, in the case of RNA viruses, genome replication. Methods for the identification of functional RNA structures in genomes typically rely on scanning analysis windows, where multiple partially-overlapping windows are used to predict RNA structures and folding metrics to deduce regions likely to form functional structure. Separate structural models are produced for each window, where the step size can greatly affect the returned model. This makes deducing unique local structures challenging, as the same nucleotides in each window can be alternatively base paired. We are presenting here a new approach where all base pairs from analysis windows are considered and weighted by favorable folding. This results in unique base pairing throughout the genome and the generation of local regions/structures that can be ranked by their propensity to form unusually thermodynamically stable folds. We applied this approach to the Zika virus (ZIKV) and HIV-1 genomes. ZIKV is linked to a variety of neurological ailments including microcephaly and Guillain–Barré syndrome and its (+)-sense RNA genome encodes two, previously described, functionally essential structured RNA regions. HIV, the cause of AIDS, contains multiple functional RNA motifs in its genome, which have been extensively studied. Our approach is able to successfully identify and model the structures of known functional motifs in both viruses, while also finding additional regions likely to form functional structures. All data have been archived at the RNAStructuromeDB (www.structurome.bb.iastate.edu), a repository of RNA folding data for humans and their pathogens.

Ribosome-induced RNA conformational changes in a viral 3\u2032-UTR sense and regulate translation levels

Structured RNA elements, programmed RNA conformational changes, and interactions between different RNA domains underlie many modes of regulating gene expression, mandating studies to understand the foundational principles that govern these phenomena. Exploring the structured 3′ untranslated region (UTR) of a viral RNA, we discovered that different contexts of the 3′-UTR confer different abilities to enhance translation of an associated open reading frame. In one context, ribosome-induced conformational changes in a ‘sensor’ RNA domain affect a separate RNA ‘functional’ domain, altering translation efficiency. The structure of the entire 3′-UTR reveals that structurally distinct domains use a spine of continuously stacked bases and a strut-like linker to create a conduit for communication within the higher-order architecture. Thus, this 3′-UTR RNA illustrates how RNA can use programmed conformational changes to sense the translation status of an upstream open reading frame, then create a tuned functional response by communicating that information to other RNA elements.

Integrative Analysis of Zika Virus Genome RNA Structure Reveals Critical Determinants of Viral Infectivity

Zika virus (ZIKV) strains can be classified into the ancestral African and contemporary Asian lineages, with the latter responsible for recent epidemics associated with neurological conditions. To understand how Asian strains lead to exacerbated disease, a crucial step is identifying genomic variations that affect infectivity and pathogenicity. Here we use two high-throughput sequencing approaches to assess RNA secondary structures and intramolecular RNA-RNA interactions invivo for the RNA genomes of Asian and African ZIKV lineages. Our analysis identified functional RNA structural elements and a functional long-range intramolecular interaction specific for the Asian epidemic strains. Mutants that disrupt this extended RNA interaction between the 5' UTR and the E protein coding region reduce virus infectivity, which is partially rescued with compensatory mutants, restoring this RNA-RNA interaction. These findings illuminate the structural basis of ZIKV regulation and provide a resource for the discovery of RNA structural elements important for ZIKV infection.

RNA Structure Elements Conserved between Mouse and 59 Other Vertebrates.

In this work, we present a computational screen conducted for functional RNA structures, resulting in over 100,000 conserved RNA structure elements found in alignments of mouse (mm10) against 59 other vertebrates. We explicitly included masked repeat regions to explore the potential of transposable elements and low-complexity regions to give rise to regulatory RNA elements. In our analysis pipeline, we implemented a four-step procedure: (i) we screened genome-wide alignments for potential structure elements using RNAz-2, (ii) realigned and refined candidate loci with LocARNA-P, (iii) scored candidates again with RNAz-2 in structure alignment mode, and (iv) searched for additional homologous loci in mouse genome that were not covered by genome alignments. The 3’-untranslated regions (3’-UTRs) of protein-coding genes and small noncoding RNAs are enriched for structures, while coding sequences are depleted. Repeat-associated loci make up about 95% of the homologous loci identified and are, as expected, predominantly found in intronic and intergenic regions. Nevertheless, we report the structure elements enriched in specific genome elements, such as 3’-UTRs and long noncoding RNAs (lncRNAs). We provide full access to our results via a custom UCSC genome browser trackhub freely available on our website (http://rna.tbi.univie.ac.at/trackhubs/#RNAz).

NMR resonance assignments for the SAM/SAH-binding riboswitch RNA bound to S-adenosylhomocysteine.

Riboswitches are structured RNA elements in the 5'-untranslated regions of bacterial mRNAs that are able to control the transcription or translation of these mRNAs in response to the specific binding of small molecules such as certain metabolites. Riboswitches that bind with high specificity to either S-adenosylmethionine (SAM) or S-adenosylhomocysteine (SAH) are widespread in bacteria. Based on differences in secondary structure and sequence these riboswitches can be grouped into a number of distinct classes. X-ray structures for riboswitch RNAs in complex with SAM or SAH established a structural basis for understanding ligand recognition and discrimination in many of these riboswitch classes. One class of riboswitches-the so-called SAM/SAH riboswitch class-binds SAM and SAH with similar affinity. However, this class of riboswitches is structurally not yet characterized and the structural basis for its unusual bispecificity is not established. In order to understand the ligand recognition mode that enables this riboswitch to bind both SAM and SAH with similar affinities, we are currently determining its structure in complex with SAH using NMR spectroscopy. Here, we present the NMR resonance assignment of the SAM/SAH binding riboswitch (env9b) in complex with SAH as a prerequisite for a solution NMR-based high-resolution structure determination.

Global pairwise RNA interaction landscapes reveal core features of protein recognition

RNA–protein interactions permeate biology. Transcription, translation, and splicing all hinge on the recognition of structured RNA elements by RNA-binding proteins. Models of RNA–protein interactions are generally limited to short linear motifs and structures because of the vast sequence sampling required to access longer elements. Here, we develop an integrated approach that calculates global pairwise interaction scores from in vitro selection and high-throughput sequencing. We examine four RNA-binding proteins of phage, viral, and human origin. Our approach reveals regulatory motifs, discriminates between regulated and non-regulated RNAs within their native genomic context, and correctly predicts the consequence of mutational events on binding activity. We design binding elements that improve binding activity in cells and infer mutational pathways that reveal permissive versus disruptive evolutionary trajectories between regulated motifs. These coupling landscapes are broadly applicable for the discovery and characterization of protein–RNA recognition at single nucleotide resolution.

Polymorphism studies on microRNA targetome of thalassemia.

Thalassemia is one of the most prevalent hemoglobin disorders. It is caused by the decreased or absent synthesis of one globin chain that leads to moderate to severe hemolytic anemia in clinical complications. Some genetic factors cause these phenotypic variations by the alteration of gene expression. MicroRNAs (miRNAs) are post-transcriptional regulators in gene expression. Therefore, variations in 3'-untranslated region (3'-UTR) of target genes may affect gene expression. It is of interest to evaluate the impact of noncoding SNPs in thalassemia related genes on miRNA: mRNA interactions in the severity of thalassemia. Polymorphisms that alter miRNA: mRNA interactions were predicted using PolymiRTS and Mirsnpscore tools. Then, the effect of predicted target SNPs on thermodynamic stability, local RNA structure and regulatory elements was investigated using RNAhybrid, RNAsnp and RegulomeDB, respectively. The molecular functions and the Biological process of candidate genes were extracted and interaction network was created. Forty-six SNPs were predicted to affect 188 miRNA interactions. These results suggest that 3'-UTR SNP may affect gene expression and cause phenotypic variation in thalassemia patients.

Extensive reshaping of bacterial operons by programmed mRNA decay.

Bacterial operons synchronize the expression of multiple genes by placing them under the control of a shared promoter. It was previously shown that polycistronic transcripts can undergo differential RNA decay, leaving some genes within the polycistron more stable than others, but the extent of regulation by differential mRNA decay or its evolutionary conservation remains unknown. Here, we find that a substantial fraction of E. coli genes display non-uniform mRNA stoichiometries despite being coded from the same operon. We further show that these altered operon stoichiometries are shaped post-transcriptionally by differential mRNA decay, which is regulated by RNA structures that protect specific regions in the transcript from degradation. These protective RNA structures are generally coded within the protein-coding regions of the regulated genes and are frequently evolutionarily conserved. Furthermore, we provide evidence that differences in ribosome densities across polycistronic transcript segments, together with the conserved structural RNA elements, play a major role in the differential decay process. Our results highlight a major role for differential mRNA decay in shaping bacterial transcriptomes.

RNA Characterization by Solid-State NMR Spectroscopy.

The structures of RNAs, which play critical roles in various biological processes, provide important clues and insights into the biological functions of these molecules. However, RNA structure determination remains a challenging topic. In recent years, magic-angle-spinning solid-state NMR (MAS SSNMR) has emerged as an alternative technique for structural and dynamic characterization of RNA. MAS SSNMR has been successfully applied to provide atomic-level structural information about several RNA molecules and RNA-protein complexes. In this Minireview, we give an overview of recent progress in the field of MAS SSNMR based RNA structural characterization, and introduce sample preparation strategies and SSNMR spectroscopic techniques that have been incorporated to identify RNA structural elements. We also highlight a few impressive examples of RNAs that have been investigated extensively by SSNMR. Finally, we briefly discuss future technical trends in the use of MAS SSNMR to facilitate RNA structure determination.

Structural divergence creates new functional features in alphavirus genomes.

Alphaviruses are mosquito-borne pathogens that cause human diseases ranging from debilitating arthritis to lethal encephalitis. Studies with Sindbis virus (SINV), which causes fever, rash, and arthralgia in humans, and Venezuelan equine encephalitis virus (VEEV), which causes encephalitis, have identified RNA structural elements that play key roles in replication and pathogenesis. However, a complete genomic structural profile has not been established for these viruses. We used the structural probing technique SHAPE-MaP to identify structured elements within the SINV and VEEV genomes. Our SHAPE-directed structural models recapitulate known RNA structures, while also identifying novel structural elements, including a new functional element in the nsP1 region of SINV whose disruption causes a defect in infectivity. Although RNA structural elements are important for multiple aspects of alphavirus biology, we found the majority of RNA structures were not conserved between SINV and VEEV. Our data suggest that alphavirus RNA genomes are highly divergent structurally despite similar genomic architecture and sequence conservation; still, RNA structural elements are critical to the viral life cycle. These findings reframe traditional assumptions about RNA structure and evolution: rather than structures being conserved, alphaviruses frequently evolve new structures that may shape interactions with host immune systems or co-evolve with viral proteins.

Know Your Enemy: Successful Bioinformatic Approaches to Predict Functional RNA Structures in Viral RNAs.

Structured RNA elements may control virus replication, transcription and translation, and their distinct features are being exploited by novel antiviral strategies. Viral RNA elements continue to be discovered using combinations of experimental and computational analyses. However, the wealth of sequence data, notably from deep viral RNA sequencing, viromes, and metagenomes, necessitates computational approaches being used as an essential discovery tool. In this review, we describe practical approaches being used to discover functional RNA elements in viral genomes. In addition to success stories in new and emerging viruses, these approaches have revealed some surprising new features of well-studied viruses e.g., human immunodeficiency virus, hepatitis C virus, influenza, and dengue viruses. Some notable discoveries were facilitated by new comparative analyses of diverse viral genome alignments. Importantly, comparative approaches for finding RNA elements embedded in coding and non-coding regions differ. With the exponential growth of computer power we have progressed from stem-loop prediction on single sequences to cutting edge 3D prediction, and from command line to user friendly web interfaces. Despite these advances, many powerful, user friendly prediction tools and resources are underutilized by the virology community.

In Vivo Visualization and Function Probing of Transport mRNPs Using Injected FIT Probes.

Fluorogenic hybridization methods, such as the use of FIT probes, enable the in vivo detection of specific mRNAs transcribed from their endogenous, genetically nonmodified loci. Here, we describe the design, synthesis and injection of nuclease resistant FIT probes into developing Drosophila oocytes to detect endogenous localizing mRNAs as wells as to probe function of structural RNA elements.

Probing of RNA structures in a positive sense RNA virus reveals selection pressures for structural elements.

In single stranded (+)-sense RNA viruses, RNA structural elements (SEs) play essential roles in the infection process from replication to encapsidation. Using selective 2′-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) and covariation analysis, we explore the structural features of the third genome segment of cucumber mosaic virus (CMV), RNA3 (2216 nt), both in vitro and in plant cell lysates. Comparing SHAPE-Seq and covariation analysis results revealed multiple SEs in the coat protein open reading frame and 3′ untranslated region. Four of these SEs were mutated and serially passaged in Nicotiana tabacum plants to identify biologically selected changes to the original mutated sequences. After passaging, loop mutants showed partial reversion to their wild-type sequence and SEs that were structurally disrupted by mutations were restored to wild-type-like structures via synonymous mutations in planta. These results support the existence and selection of virus open reading frame SEs in the host organism and provide a framework for further studies on the role of RNA structure in viral infection. Additionally, this work demonstrates the applicability of high-throughput chemical probing in plant cell lysates and presents a new method for calculating SHAPE reactivities from overlapping reverse transcriptase priming sites.

Structural and functional conservation of cis-acting RNA elements in coronavirus 5'-terminal genome regions

Structure predictions suggest a partial conservation of RNA structure elements in coronavirus terminal genome regions. Here, we determined the structures of stem-loops (SL) 1 and 2 of two alphacoronaviruses, human coronavirus (HCoV) 229E and NL63, by RNA structure probing and studied the functional relevance of these putative cis-acting elements. HCoV-229E SL1 and SL2 mutants generated by reverse genetics were used to study the effects on viral replication of single-nucleotide substitutions predicted to destabilize the SL1 and SL2 structures. The data provide conclusive evidence for the critical role of SL1 and SL2 in HCoV-229E replication and, in some cases, revealed parallels with previously characterized betacoronavirus SL1 and SL2 elements. Also, we were able to rescue viable HCoV-229E mutants carrying replacements of SL2 with equivalent betacoronavirus structural elements. The data obtained in this study reveal a remarkable degree of structural and functional conservation of 5′-terminal RNA structural elements across coronavirus genus boundaries.

Functional 5' UTR mRNA structures in eukaryotic translation regulation and how to find them.

RNA molecules can fold into intricate shapes that can provide an additional layer of control of gene expression beyond that of their sequence. In this Review, we discuss the current mechanistic understanding of structures in 5' untranslated regions (UTRs) of eukaryotic mRNAs and the emerging methodologies used to explore them. These structures may regulate cap-dependent translation initiation through helicase-mediated remodelling of RNA structures and higher-order RNA interactions, as well as cap-independent translation initiation through internal ribosome entry sites (IRESs), mRNA modifications and other specialized translation pathways. We discuss known 5' UTR RNA structures and how new structure probing technologies coupled with prospective validation, particularly compensatory mutagenesis, are likely to identify classes of structured RNA elements that shape post-transcriptional control of gene expression and the development of multicellular organisms.

Trick or TREAT: A Scary-Good New Approach for Single-Molecule mRNA Decay Analysis

In this issue of Molecular Cell, Horvathova etal. (2017) have developed a powerful approach to single-molecule assessment of RNA decay in living cells by exploiting the ability of flavivirus RNA structural elements to trap XRN1 decay intermediates in dual-labeled reporter constructs.

The DEAD-Box Protein CYT-19 Uses Arginine Residues in Its C-Tail To Tether RNA Substrates.

DEAD-box proteins are nonprocessive RNA helicases that play diverse roles in cellular processes. The Neurospora crassa DEAD-box protein CYT-19 promotes mitochondrial group I intron splicing and functions as a general RNA chaperone. CYT-19 includes a disordered, arginine-rich "C-tail" that binds RNA, positioning the helicase core to capture and unwind nearby RNA helices. Here we probed the C-tail further by varying the number and positions of arginines within it. We found that removing sets of as few as four of the 11 arginines reduced RNA unwinding activity (kcat/KM) to a degree equivalent to that seen upon removal of the C-tail, suggesting that a minimum or "threshold" number of arginines is required. In addition, a mutant with 16 arginines displayed RNA unwinding activity greater than that of wild-type CYT-19. The C-tail modifications impacted unwinding only of RNA helices within constructs that included an adjacent helix or structured RNA element that would allow C-tail binding, indicating that the helicase core remained active in the mutants. In addition, changes in RNA unwinding efficiency of the mutants were mirrored by changes in functional RNA affinity, as determined from the RNA concentration dependence of ATPase activity, suggesting that the C-tail functions primarily to increase RNA affinity. Interestingly, the salt concentration dependence of RNA unwinding activity is unaffected by C-tail composition, suggesting that the C-tail uses primarily hydrogen bonding, not electrostatic interactions, to bind double-stranded RNA. Our results provide insights into how an unstructured C-tail contributes to DEAD-box protein activity and suggest parallels with other families of RNA- and DNA-binding proteins.

Structures of riboswitch RNA reaction states by mix-and-inject XFEL serial crystallography

Riboswitches are structural RNA elements that are generally located in the 5' untranslated region of messenger RNA. During regulation of gene expression, ligand binding to the aptamer domain of a riboswitch triggers a signal to the downstream expression platform. A complete understanding of the structural basis of this mechanism requires the ability to study structural changes over time. Here we use femtosecond X-ray free electron laser (XFEL) pulses to obtain structural measurements from crystals so small that diffusion of a ligand can be timed to initiate a reaction before diffraction. We demonstrate this approach by determining four structures of the adenine riboswitch aptamer domain during the course of a reaction, involving two unbound apo structures, one ligand-bound intermediate, and the final ligand-bound conformation. These structures support a reaction mechanism model with at least four states and illustrate the structural basis of signal transmission. The three-way junction and the P1 switch helix of the two apo conformers are notably different from those in the ligand-bound conformation. Our time-resolved crystallographic measurements with a 10-second delay captured the structure of an intermediate with changes in the binding pocket that accommodate the ligand. With at least a 10-minute delay, the RNA molecules were fully converted to the ligand-bound state, in which the substantial conformational changes resulted in conversion of the space group. Such notable changes in crystallo highlight the important opportunities that micro- and nanocrystals may offer in these and similar time-resolved diffraction studies. Together, these results demonstrate the potential of 'mix-and-inject' time-resolved serial crystallography to study biochemically important interactions between biomacromolecules and ligands, including those that involve large conformational changes.