RNA Structure Probing Research Articles

RASP v2.0: an updated atlas for RNA structure probing data.

RNA molecules function in numerous biological processes by folding into intricate structures. Here we present RASP v2.0, an updated database for RNA structure probing data featuring a substantially expanded collection of datasets along with enhanced online structural analysis functionalities. Compared to the previous version, RASP v2.0 includes the following improvements: (i) the number of RNA structure datasets has increased from 156 to 438, comprising 216 transcriptome-wide RNA structure datasets, 141 target-specific RNA structure datasets, and 81 RNA-RNA interaction datasets, thereby broadening species coverage from 18 to 24, (ii) a deep learning-based model has been implemented to impute missing structural signals for 59 transcriptome-wide RNA structure datasets with low structure score coverage, significantly enhancing data quality, particularly for low-abundance RNAs, (iii) three new online analysis modules have been deployed to assist RNA structure studies, including missing structure score imputation, RNA secondary and tertiary structure prediction, and RNA binding protein (RBP) binding prediction. By providing a resource of much more comprehensive RNA structure data, RASP v2.0 is poised to facilitate the exploration of RNA structure-function relationships across diverse biological processes. RASP v2.0 is freely accessible at http://rasp2.zhanglab.net/.

Template switching enables chemical probing of native RNA structures.

RNAs are often studied in non-native sequence contexts to facilitate structural studies. However, seemingly innocuous changes to an RNA sequence may perturb the native structure and generate inaccurate or ambiguous structural models. To facilitate the investigation of native RNA secondary structure by selective 2' hydroxyl acylation analyzed by primer extension (SHAPE), we engineered an approach that couples minimal enzymatic steps to RNA chemical probing and mutational profiling (MaP) reverse transcription (RT) methods - a process we call template switching and mutational profiling (Switch-MaP). In Switch-MaP, RT templates and additional library sequences are added post-probing through ligation and template switching, capturing reactivities for every nucleotide. For a candidate SAM-I riboswitch, we compared RNA structure models generated by the Switch-MaP approach to those of traditional primer-based MaP, including RNAs with or without appended structure cassettes. Primer-based MaP masked reactivity data in the 5' and 3' ends of the RNA, producing ambiguous ensembles inconsistent with the conserved SAM-I riboswitch secondary structure. Structure cassettes enabled unambiguous modeling of an aptamer-only construct but introduced non-native interactions in the full length riboswitch. In contrast, Switch-MaP provided reactivity data for all nucleotides in each RNA and enabled unambiguous modeling of secondary structure, consistent with the conserved SAM-I fold. Switch-MaP is a straightforward alternative approach to primer-based and cassette-based chemical probing methods that precludes primer masking and the formation of alternative secondary structures due to non-native sequence elements.

Sequential structure probing of cotranscriptional RNA folding intermediates.

Cotranscriptional RNA folding pathways typically involve the sequential formation of folding intermediates. Existing methods for cotranscriptional RNA structure probing map the structure of nascent RNA in the context of a terminally arrested transcription elongation complex. Consequently, the rearrangement of RNA structures as nucleotides are added to the transcript can be inferred but is not assessed directly. To address this limitation, we have developed linked- m ultipoint T ranscription E longation C omplex RNA structure prob ing (TECprobe-LM), which assesses the cotranscriptional rearrangement of RNA structures by sequentially positioning E. coli RNAP at two or more points within a DNA template so that nascent RNA can be chemically probed. We validated TECprobe-LM by measuring known folding events that occur within the E. coli signal recognition particle RNA, Clostridium beijerinckii pfl ZTP riboswitch, and Bacillus cereus crcB fluoride riboswitch folding pathways. Our findings establish TECprobe-LM as a strategy for detecting cotranscriptional RNA folding events directly using chemical probing.

Probing RNA structure and dynamics using cutting-edge biophysical techniques

Probing RNA structure and dynamics using cutting-edge biophysical techniques

Phase separation and viral factories: unveiling the physical processes supporting RNA packaging in dsRNA viruses.

Understanding of the physicochemical properties and functions of biomolecular condensates has rapidly advanced over the past decade. More recently, many RNA viruses have been shown to form cytoplasmic replication factories, or viroplasms, via phase separation of their components, akin to numerous cellular membraneless organelles. Notably, diverse viruses from the Reoviridae family containing 10-12 segmented double-stranded RNA genomes induce the formation of viroplasms in infected cells. Little is known about the inner workings of these membraneless cytoplasmic inclusions and how they may support stoichiometric RNA assembly in viruses with segmented RNA genomes, raising questions about the roles of phase separation in coordinating viral genome packaging. Here, we discuss how the molecular composition of viroplasms determines their properties, highlighting the interplay between RNA structure, RNA remodelling, and condensate self-organisation. Advancements in RNA structural probing and theoretical modelling of condensates can reveal the mechanisms through which these ribonucleoprotein complexes support the selective enrichment and stoichiometric assembly of distinct viral RNAs.

SHAPEwarp-web: sequence-agnostic search for structurally homologous RNA regions across databases of chemical probing data.

RNA molecules perform a variety of functions in cells, many of which rely on their secondary and tertiary structures. Chemical probing methods coupled with high-throughput sequencing have significantly accelerated the mapping of RNA structures, and increasingly large datasets of transcriptome-wide RNA chemical probing data are becoming available. Analogously to what has been done for decades in the protein world, this RNA structural information can be leveraged to aid the discovery of structural similarity to a known RNA (or RNA family), which, in turn, can inform about the function of transcripts. We have previously developed SHAPEwarp, a sequence-agnostic method for the search of structurally homologous RNA segments in a database of reactivity profiles derived from chemical probing experiments. In its original implementation, however, SHAPEwarp required substantial computational resources, even for moderately sized databases, as well as significant Linux command line know-how. To address these limitations, we introduce here SHAPEwarp-web, a user-friendly web interface to rapidly query large databases of RNA chemical probing data for structurally similar RNAs. Aside from featuring a completely rewritten core, which speeds up by orders of magnitude the search inside large databases, the web server hosts several high-quality chemical probing databases across multiple species. SHAPEwarp-web is available fromhttps://shapewarp.incarnatolab.com.

An RNA pseudoknot mediates toxin translation and antitoxin inhibition

Type I toxin-antitoxin systems (T1TAs) are bipartite bacterial loci encoding a growth-inhibitory toxin and an antitoxin small RNA (sRNA). In many of these systems, the transcribed toxin mRNA is translationally inactive, but becomes translation-competent upon ribonucleolytic processing. The antitoxin sRNA targets the processed mRNA to inhibit its translation. This two-level control mechanism prevents cotranscriptional translation of the toxin and allows its synthesis only when the antitoxin is absent. Contrary to this, we found that the timP mRNA of the timPR T1TA locus does not undergo enzymatic processing. Instead, the full-length timP transcript is both translationally active and can be targeted by the antitoxin TimR. Thus, tight control in this system relies on a noncanonical mechanism. Based on the results from in vitro binding assays, RNA structure probing, and cell-free translation experiments, we suggest that timP mRNA adopts mutually exclusive structural conformations. The active form uniquely possesses an RNA pseudoknot structure which is essential for translation initiation. TimR preferentially binds to the active conformation, which leads to pseudoknot destabilization and inhibited translation. Based on this, we propose a model in which "structural processing" of timP mRNA enables tight inhibition by TimR in nonpermissive conditions, and TimP synthesis only upon TimR depletion.

A new reagent for in vivo structure probing of RNA G and U residues that improves RNA structure prediction alone and combined with DMS.

A key to understanding the roles of RNA in regulating gene expression is knowing their structures in vivo. One way to obtain this information is through probing the structures of RNA with chemicals. To probe RNA structure directly in cells, membrane-permeable reagents that modify the Watson-Crick (WC) face of unpaired nucleotides can be used. Although dimethyl sulfate (DMS) has led to substantial insight into RNA structure, it has limited nucleotide specificity in vivo, with WC face reactivity only at adenine (A) and cytosine (C) at neutral pH. The reagent 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) was recently shown to modify the WC face of guanine (G) and uracil (U). Although useful at lower concentrations in experiments that measure chemical modifications by reverse transcription (RT) stops, at higher concentrations necessary for detection by mutational profiling (MaP), EDC treatment leads to degradation of RNA. Here, we demonstrate EDC-stimulated degradation of RNA in Gram-negative and Gram-positive bacteria. In an attempt to overcome these limitations, we developed a new carbodiimide reagent, 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide methiodide (ETC), which we show specifically modifies unpaired Gs and Us in vivo without substantial degradation of RNA. We establish ETC as a probe for MaP and optimize the RT conditions and computational analysis in Escherichia coli Importantly, we demonstrate the utility of ETC as a probe for improving RNA structure prediction both alone and with DMS.

Structural atlas of human primary microRNAs generated by SHAPE-MaP

MicroRNA (miRNA) maturation is critically dependent on structural features of primary transcripts (pri-miRNAs). However, the scarcity of determined pri-miRNA structures has limited our understanding of miRNA maturation. Here, we employed selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP), a high-throughput RNA structure probing method, to unravel the secondary structures of 476 high-confidence human pri-miRNAs. Our SHAPE-based structures diverge substantially from those inferred solely from computation, particularly in the apical loop and basal segments, underlining the need for experimental data in RNA structure prediction. By comparing the structures with high-throughput processing data, we determined the optimal structural features of pri-miRNAs. The sequence determinants are influenced substantially by their structural contexts. Moreover, we identified an element termed the bulged GWG motif (bGWG) with a 3' bulge in the lower stem, which promotes processing. Our structure-function mapping better annotates the determinants of pri-miRNA processing and offers practical implications for designing small hairpin RNAs and predicting the impacts of miRNA mutations.

RNAvigate: efficient exploration of RNA chemical probing datasets.

Chemical probing technologies enable high-throughput examination of diverse structural features of RNA, including local nucleotide flexibility, RNA secondary structure, protein and ligand binding, through-space interaction networks, and multistate structural ensembles. Deep understanding of RNA structure-function relationships typically requires evaluating a system under structure- and function-altering conditions, linking these data with additional information, and visualizing multilayered relationships. Current platforms lack the broad accessibility, flexibility and efficiency needed to iterate on integrative analyses of these diverse, complex data. Here, we share the RNA visualization and graphical analysis toolset RNAvigate, a straightforward and flexible Python library that automatically parses 21 standard file formats (primary sequence annotations, per- and internucleotide data, and secondary and tertiary structures) and outputs 18 plot types. RNAvigate enables efficient exploration of nuanced relationships between multiple layers of RNA structure information and across multiple experimental conditions. Compatibility with Jupyter notebooks enables nonburdensome, reproducible, transparent and organized sharing of multistep analyses and data visualization strategies. RNAvigate simplifies and accelerates discovery and characterization of RNA-centric functions in biology.

Single-cell probing of RNA structure.

Single-cell probing of RNA structure.

RNA structure: implications in viral infections and neurodegenerative diseases

AbstractRNA is an intermediary between DNA and protein, a catalyzer of biochemical reactions, and a regulator of genes and transcripts. RNA structures are essential for complicated functions. Recent years have witnessed rapid advancements in RNA secondary structure probing techniques. These technological strides provided comprehensive insights into RNA structures, which significantly contributed to our understanding of diverse cellular regulatory processes, including gene regulation, epigenetic regulation, and post-transactional regulation. Meanwhile, they have facilitated the creation of therapeutic tools for tackling human diseases. Despite their therapeutic applications, RNA structure probing methods also offer a promising avenue for exploring the mechanisms of human diseases, potentially providing the key to overcoming existing research constraints and obtaining the in-depth information necessary for a deeper understanding of disease mechanisms.

RNA Secondary Structure Modeling Following the IPANEMAP Workflow.

The structure of RNA molecules and their complexes are crucial for understanding biology at the molecular level. Resolving these structures holds the key to understanding their manifold structure-mediated functions ranging from regulating gene expression to catalyzing biochemical processes. Predicting RNA secondary structure is a prerequisite and a key step to accurately model their three dimensional structure. Although dedicated modelling software are making fast and significant progresses, predicting an accurate secondary structure from the sequence remains a challenge. Their performance can be significantly improved by the incorporation of experimental RNA structure probing data. Many different chemical and enzymatic probes have been developed; however, only one set of quantitative data can be incorporated as constraints for computer-assisted modelling. IPANEMAP is a recent workflow based on RNAfold that can take into account several quantitative or qualitative data sets to model RNA secondary structure. This chapter details the methods for popular chemical probing (DMS, CMCT, SHAPE-CE, and SHAPE-Map) and the subsequent analysis and structure prediction using IPANEMAP.

Observation of coordinated RNA folding events by systematic cotranscriptional RNA structure probing

RNA begins to fold as it is transcribed by an RNA polymerase. Consequently, RNA folding is constrained by the direction and rate of transcription. Understanding how RNA folds into secondary and tertiary structures therefore requires methods for determining the structure of cotranscriptional folding intermediates. Cotranscriptional RNA chemical probing methods accomplish this by systematically probing the structure of nascent RNA that is displayed from an RNA polymerase. Here, we describe a concise, high-resolution cotranscriptional RNA chemical probing procedure called variable length Transcription Elongation Complex RNA structure probing (TECprobe-VL). We demonstrate the accuracy and resolution of TECprobe-VL by replicating and extending previous analyses of ZTP and fluoride riboswitch folding and mapping the folding pathway of a ppGpp-sensing riboswitch. In each system, we show that TECprobe-VL identifies coordinated cotranscriptional folding events that mediate transcription antitermination. Our findings establish TECprobe-VL as an accessible method for mapping cotranscriptional RNA folding pathways.

An intronic RNA element modulates Factor VIII exon-16 splicing.

Pathogenic variants in the human Factor VIII (F8) gene cause Hemophilia A (HA). Here, we investigated the impact of 97 HA-causing single-nucleotide variants on the splicing of 11 exons from F8. For the majority of F8 exons, splicing was insensitive to the presence of HA-causing variants. However, splicing of several exons, including exon-16, was impacted by variants predicted to alter exonic splicing regulatory sequences. Using exon-16 as a model, we investigated the structure-function relationship of HA-causing variants on splicing. Intriguingly, RNA chemical probing analyses revealed a three-way junction structure at the 3'-end of intron-15 (TWJ-3-15) capable of sequestering the polypyrimidine tract. We discovered antisense oligonucleotides (ASOs) targeting TWJ-3-15 partially rescue splicing-deficient exon-16 variants by increasing accessibility of the polypyrimidine tract. The apical stem loop region of TWJ-3-15 also contains two hnRNPA1-dependent intronic splicing silencers (ISSs). ASOs blocking these ISSs also partially rescued splicing. When used in combination, ASOs targeting both the ISSs and the region sequestering the polypyrimidine tract, fully rescue pre-mRNA splicing of multiple HA-linked variants of exon-16. Together, our data reveal a putative RNA structure that sensitizes F8 exon-16 to aberrant splicing.

Variable eIF4E-binding sites and their synergistic effect on cap-independent translation in a novel IRES of wheat yellow mosaic virus RNA2 isolates

Some viral proteins are translated cap-independently via the internal ribosome entry site (IRES), which maintains conservative characteristic among different isolates of the same virus species. However, IRES activity showed a 7-fold variance in RNA2 of wheat yellow mosaic virus (WYMV) HC and LYJN isolates in this study. Based on RNA structure probing and mutagenesis assay, the loosened middle stem of H1 and the hepta-nucleotide top loop of H2 in the LYJN isolate synergistically ensured higher IRES activity than that in the HC isolate. In addition, the conserved top loop of H1 ensured basic IRES activity in HC and LYJN isolates. WYMV RNA2 5′-UTR specifically interacted with the wheat eIF4E, accomplished by the top loop of H1 in the HC isolate or the top loop of H1 and H2 in the LYJN isolate. The high IRES activity of the WYMV RNA2 LYJN isolate was regulated by two eIF4E-binding sites, which showed a synergistic effect mediated by the proximity of the H1 and H2 top loops owing to the flexibility of the middle stem in H1. This report presents a novel evolution pattern of IRES, which altered the number of eIF4E-binding sites to regulate IRES activity.

In vivo-like nearest neighbor parameters improve prediction of fractional RNA base-pairing in cells.

We conducted a thermodynamic analysis of RNA stability in Eco80 artificial cytoplasm, which mimics in vivo conditions, and compared it to transcriptome-wide probing of mRNA. Eco80 contains 80% of Escherichia coli metabolites, with biological concentrations of metal ions, including 2 mM free Mg2+ and 29 mM metabolite-chelated Mg2+. Fluorescence-detected binding isotherms (FDBI) were used to conduct a thermodynamic analysis of 24 RNA helices and found that these helices, which have an average stability of -12.3 kcal/mol, are less stable by ΔΔGo37 ∼1 kcal/mol. The FDBI data was used to determine a set of Watson-Crick free energy nearestneighbor parameters (NNPs), which revealed that Eco80 reduces the stability of three NNPs. This information was used to adjust the NN model using the RNAstructure package. The in vivo-like adjustments have minimal effects on the prediction of RNA secondary structures determined in vitro and in silico, but markedly improve prediction of fractional RNA base pairing in E. coli, as benchmarked with our in vivo DMS and EDC RNA chemical probing data. In summary, our thermodynamic and chemical probing analyses of RNA helices indicate that RNA secondary structures are less stable in cells than in artificially stable in vitro buffer conditions.

Genome-wide probing of eukaryotic nascent RNA structure elucidates cotranscriptional folding and its antimutagenic effect

The transcriptional intermediates of RNAs fold into secondary structures with multiple regulatory roles, yet the details of such cotranscriptional RNA folding are largely unresolved in eukaryotes. Here, we present eSPET-seq (Structural Probing of Elongating Transcripts in eukaryotes), a method to assess the cotranscriptional RNA folding in Saccharomyces cerevisiae. Our study reveals pervasive structural transitions during cotranscriptional folding and overall structural similarities between nascent and mature RNAs. Furthermore, a combined analysis with genome-wide R-loop and mutation rate approximations provides quantitative evidence for the antimutator effect of nascent RNA folding through competitive inhibition of the R-loops, known to facilitate transcription-associated mutagenesis. Taken together, we present an experimental evaluation of cotranscriptional folding in eukaryotes and demonstrate the antimutator effect of nascent RNA folding. These results suggest genome-wide coupling between the processing and transmission of genetic information through RNA folding.

In vivo structure probing of RNA in Archaea: novel insights into the ribosome structure of Methanosarcina acetivorans.

Structure probing combined with next-generation sequencing (NGS) has provided novel insights into RNA structure-function relationships. To date, such studies have focused largely on bacteria and eukaryotes, with little attention given to the third domain of life, archaea. Furthermore, functional RNAs have not been extensively studied in archaea, leaving open questions about RNA structure and function within this domain of life. With archaeal species being diverse and having many similarities to both bacteria and eukaryotes, the archaea domain has the potential to be an evolutionary bridge. In this study, we introduce a method for probing RNA structure in vivo in the archaea domain of life. We investigated the structure of ribosomal RNA (rRNA) from Methanosarcina acetivorans, a well-studied anaerobic archaeal species, grown with either methanol or acetate. After probing the RNA in vivo with dimethyl sulfate (DMS), Structure-seq2 libraries were generated, sequenced, and analyzed. We mapped the reactivity of DMS onto the secondary structure of the ribosome, which we determined independently with comparative analysis, and confirmed the accuracy of DMS probing in M. acetivorans Accessibility of the rRNA to DMS in the two carbon sources was found to be quite similar, although some differences were found. Overall, this study establishes the Structure-seq2 pipeline in the archaea domain of life and informs about ribosomal structure within M. acetivorans.

The oxidative stress response, in particular the katY gene, is temperature-regulated in Yersinia pseudotuberculosis.

Pathogenic bacteria, such as Yersinia pseudotuberculosis encounter reactive oxygen species (ROS) as one of the first lines of defense in the mammalian host. In return, the bacteria react by mounting an oxidative stress response. Previous global RNA structure probing studies provided evidence for temperature-modulated RNA structures in the 5'-untranslated region (5'-UTR) of various oxidative stress response transcripts, suggesting that opening of these RNA thermometer (RNAT) structures at host-body temperature relieves translational repression. Here, we systematically analyzed the transcriptional and translational regulation of ROS defense genes by RNA-sequencing, qRT-PCR, translational reporter gene fusions, enzymatic RNA structure probing and toeprinting assays. Transcription of four ROS defense genes was upregulated at 37°C. The trxA gene is transcribed into two mRNA isoforms, of which the most abundant short one contains a functional RNAT. Biochemical assays validated temperature-responsive RNAT-like structures in the 5'-UTRs of sodB, sodC and katA. However, they barely conferred translational repression in Y. pseudotuberculosis at 25°C suggesting partially open structures available to the ribosome in the living cell. Around the translation initiation region of katY we discovered a novel, highly efficient RNAT that was primarily responsible for massive induction of KatY at 37°C. By phenotypic characterization of catalase mutants and through fluorometric real-time measurements of the redox-sensitive roGFP2-Orp1 reporter in these strains, we revealed KatA as the primary H2O2 scavenger. Consistent with the upregulation of katY, we observed an improved protection of Y. pseudotuberculosis at 37°C. Our findings suggest a multilayered regulation of the oxidative stress response in Yersinia and an important role of RNAT-controlled katY expression at host body temperature.