Genome-Wide Approaches for RNA Structure Probing

RNA structure: implications in viral infections and neurodegenerative diseases

RNA 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.

Advances and Trends in Omics Technology Development

The human history has witnessed the rapid development of technologies such as high-throughput sequencing and mass spectrometry that led to the concept of “omics” and methodological advancement in systematically interrogating a cellular system. Yet, the ever-growing types of molecules and regulatory mechanisms being discovered have been persistently transforming our understandings on the cellular machinery. This renders cell omics seemingly, like the universe, expand with no limit and our goal toward the complete harness of the cellular system merely impossible. Therefore, it is imperative to review what has been done and is being done to predict what can be done toward the translation of omics information to disease control with minimal cell perturbation. With a focus on the “four big omics,” i.e., genomics, transcriptomics, proteomics, metabolomics, we delineate hierarchies of these omics together with their epiomics and interactomics, and review technologies developed for interrogation. We predict, among others, redoxomics as an emerging omics layer that views cell decision toward the physiological or pathological state as a fine-tuned redox balance.

Lead-seq: transcriptome-wide structure probing in vivo using lead(II) ions

The dynamic conformation of RNA molecules within living cells is key to their function. Recent advances in probing the RNA structurome in vivo, including the use of SHAPE (Selective 2'-Hydroxyl Acylation analyzed by Primer Extension) or kethoxal reagents or DMS (dimethyl sulfate), provided unprecedented insights into the architecture of RNA molecules in the living cell. Here, we report the establishment of lead probing in a global RNA structuromics approach. In order to elucidate the transcriptome-wide RNA landscape in the enteric pathogen Yersinia pseudotuberculosis, we combined lead(II) acetate-mediated cleavage of single-stranded RNA regions with high-throughput sequencing. This new approach, termed 'Lead-seq', provides structural information independent of base identity. We show that the method recapitulates secondary structures of tRNAs, RNase P RNA, tmRNA, 16S rRNA and the rpsT 5'-untranslated region, and that it reveals global structural features of mRNAs. The application of Lead-seq to Y. pseudotuberculosis cells grown at two different temperatures unveiled the first temperature-responsive in vivo RNA structurome of a bacterial pathogen. The translation of candidate genes derived from this approach was confirmed to be temperature regulated. Overall, this study establishes Lead-seq as complementary approach to interrogate intracellular RNA structures on a global scale.

Repeat-associated RNA structure and aberrant splicing

Over 30 hereditary disorders attributed to the expansion of microsatellite repeats have been identified. Despite variant nucleotide content, number of consecutive repeats, and different locations in the genome, many of these diseases have pathogenic RNA gain-of-function mechanisms. The repeat-containing RNAs can form structures in vitro predicted to contribute to the disease through assembly of intracellular RNA aggregates termed foci. The expanded repeat RNAs within these foci sequester RNA binding proteins (RBPs) with important roles in the regulation of RNA metabolism, most notably alternative splicing (AS). These deleterious interactions lead to downstream alterations in transcriptome-wide AS directly linked with disease symptoms. This review summarizes existing knowledge about the association between the repeat RNA structures and RBPs as well as the resulting aberrant AS patterns, specifically in the context of myotonic dystrophy. The connection between toxic, structured RNAs and dysregulation of AS in other repeat expansion diseases is also discussed. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

UPA-seq: prediction of functional lncRNAs using differential sensitivity to UV crosslinking

While a large number of long noncoding RNAs (lncRNAs) are transcribed from the genome of higher eukaryotes, systematic prediction of their functionality has been challenging due to the lack of conserved sequence motifs or structures. Assuming that some lncRNAs function as large ribonucleoprotein complexes and thus are easily crosslinked to proteins upon UV irradiation, we performed RNA-seq analyses of RNAs recovered from the aqueous phase after UV irradiation and phenol-chloroform extraction (UPA-seq). As expected, the numbers of UPA-seq reads mapped to known functional lncRNAs were remarkably reduced upon UV irradiation. Comparison with ENCODE eCLIP data revealed that lncRNAs that exhibited greater decreases upon UV irradiation preferentially associated with proteins containing prion-like domains (PrLDs). Fluorescent in situ hybridization (FISH) analyses revealed the nuclear localization of novel functional lncRNA candidates, including one that accumulated at the site of transcription. We propose that UPA-seq provides a useful tool for the selection of lncRNA candidates to be analyzed in depth in subsequent functional studies.

High-throughput in vivo mapping of RNA accessible interfaces to identify functional sRNA binding sites

Herein we introduce a high-throughput method, INTERFACE, to reveal the capacity of contiguous RNA nucleotides to establish in vivo intermolecular RNA interactions for the purpose of functional characterization of intracellular RNA. INTERFACE enables simultaneous accessibility interrogation of an unlimited number of regions by coupling regional hybridization detection to transcription elongation outputs measurable by RNA-seq. We profile over 900 RNA interfaces in 71 validated, but largely mechanistically under-characterized, Escherichia coli sRNAs in the presence and absence of a global regulator, Hfq, and find that two-thirds of tested sRNAs feature Hfq-dependent regions. Further, we identify in vivo hybridization patterns that hallmark functional regions to uncover mRNA targets. In this way, we biochemically validate 25 mRNA targets, many of which are not captured by typically tested, top-ranked computational predictions. We additionally discover direct mRNA binding activity within the GlmY terminator, highlighting the information value of high-throughput RNA accessibility data.

Mapping RNA accessibility is valuable for identifying functional/regulatory RNA regions. Here the authors introduce INTERFACE, an intracellular method that quantifies antisense hybridization efficacy of any number of RNA regions simultaneously via a transcriptional elongation output, measurable via RNA-seq

Ribothrypsis, a novel process of canonical mRNA decay, mediates ribosome-phased mRNA endonucleolysis

mRNAs transmit the genetic information that dictates protein production and are a nexus for numerous pathways that regulate gene expression. The prevailing view of canonical mRNA decay is that it is mediated by deadenylation and decapping followed by exonucleolysis from the 3' and 5' ends. By developing Akron-seq, a novel approach that captures the native 3' and 5' ends of capped and polyadenylated RNAs, respectively, we show that canonical human mRNAs are subject to repeated cotranslational and ribosome-phased endonucleolytic cuts at the exit site of the mRNA ribosome channel, in a process that we term ribothrypsis. We uncovered RNA G quadruplexes among likely ribothrypsis triggers and show that ribothrypsis is a conserved process. Strikingly, we found that mRNA fragments are abundant in living cells and thus have important implications for the interpretation of experiments, such as RNA-seq, that rely on the assumption that mRNAs exist largely as full-length molecules in vivo.

Production of HIV-1 vif mRNA Is Modulated by Natural Nucleotide Variations and SLSA1 RNA Structure in SA1D2prox Genomic Region

Genomic RNA of HIV-1 contains localized structures critical for viral replication. Its structural analysis has demonstrated a stem-loop structure, SLSA1, in a nearby region of HIV-1 genomic splicing acceptor 1 (SA1). We have previously shown that the expression level of vif mRNA is considerably altered by some natural single-nucleotide variations (nSNVs) clustering in SLSA1 structure. In this study, besides eleven nSNVs previously identified by us, we totally found nine new nSNVs in the SLSA1-containing sequence from SA1, splicing donor 2, and through to the start codon of Vif that significantly affect the vif mRNA level, and designated the sequence SA1D2prox (142 nucleotides for HIV-1 NL4-3). We then examined by extensive variant and mutagenesis analyses how SA1D2prox sequence and SLSA1 secondary structure are related to vif mRNA level. While the secondary structure and stability of SLSA1 was largely changed by nSNVs and artificial mutations introduced to restore the original NL4-3 form from altered ones by nSNVs, no clear association of the two SLSA1 properties with vif mRNA level was observed. In contrast, when naturally occurring SA1D2prox sequences that contain multiple nSNVs were examined, we attained significant inverse correlation between the vif level and SLSA1 stability. These results may suggest that SA1D2prox sequence adapts over time, and also that the altered SA1D2prox sequence, SLSA1 stability, and vif level are mutually related. In total, we show here that the entire SA1D2prox sequence and SLSA1 stability critically contribute to the modulation of vif mRNA level.

RNA Regulations and Functions Decoded by Transcriptome-wide RNA Structure Probing

RNA folds into intricate structures that are crucial for its functions and regulations. To date, a multitude of approaches for probing structures of the whole transcriptome, i.e., RNA structuromes, have been developed. Applications of these approaches to different cell lines and tissues have generated a rich resource for the study of RNA structure–function relationships at a systems biology level. In this review, we first introduce the designs of these methods and their applications to study different RNA structuromes. We emphasize their technological differences especially their unique advantages and caveats. We then summarize the structural insights in RNA functions and regulations obtained from the studies of RNA structuromes. And finally, we propose potential directions for future improvements and studies.

Interpreting Reverse Transcriptase Termination and Mutation Events for Greater Insight into the Chemical Probing of RNA

Chemical probing has the power to provide insight into RNA conformation in vivo and in vitro, but interpreting the results depends on methods to detect the chemically modified nucleotides. Traditionally, the presence of modified bases was inferred from their ability to halt reverse transcriptase during primer extension and the locations of termination sites observed by electrophoresis or sequencing. More recently, modification-induced mutations have been used as a readout for chemical probing data. Given the variable propensity for mismatch incorporation and read-through with different reverse transcriptases, we examined how termination and mutation events compare to each other in the same chemical probing experiments. We found that mutations and terminations induced by dimethyl sulfate probing are both specific for methylated bases, but these two measures have surprisingly little correlation and represent largely nonoverlapping indicators of chemical modification data. We also show that specific biases for modified bases depend partly on local sequence context and that different reverse transcriptases show different biases toward reading a modification as a stop or a mutation. These results support approaches that incorporate analysis of both termination and mutation events into RNA probing experiments.

A combined sequence and structure based method for discovering enriched motifs in RNA from in vivo binding data

RNA binding proteins (RBPs) play an important role in regulating many processes in the cell. RBPs often recognize their RNA targets in a specific manner. In addition to the RNA primary sequence, the structure of the RNA has been shown to play a central role in RNA recognition by RBPs. In recent years, many experimental approaches, both in vitro and in vivo, were developed and employed to identify and characterize RBP targets and extract their binding specificities. In vivo binding techniques, such as CrossLinking and ImmunoPrecipitation (CLIP)-based methods, enable the characterization of protein binding sites on RNA targets. However, these methods do not provide information regarding the structural preferences of the protein. While methods to obtain the structure of RNA are available, inferring both the sequence and the structure preferences of RBPs remains a challenge. Here we present SMARTIV, a novel computational tool for discovering combined sequence and structure binding motifs from in vivo RNA binding data relying on the sequences of the target sites, the ranking of their binding scores and their predicted secondary structure. The combined motifs are provided in a unified representation that is informative and easy for visual perception. We tested the method on CLIP-seq data from different platforms for a variety of RBPs. Overall, we show that our results are highly consistent with known binding motifs of RBPs, offering additional information on their structural preferences.