Massive parallel-sequencing-based hydroxyl radical probing of RNA accessibility

Navigating the landscape of epitranscriptomics and host immunity

RNA modifications, also termed epitranscriptomic marks, encompass chemical alterations to individual nucleotides, including processes such as methylation and editing. These marks contribute to a wide range of biological processes, many of which are related to host immune system defense. The functions of immune-related RNA modifications can be categorized into three main groups: regulation of immunogenic RNAs, control of genes involved in innate immune response, and facilitation of adaptive immunity. Here, we provide an overview of recent research findings that elucidate the contributions of RNA modifications to each of these processes. We also discuss relevant methods for genome-wide identification of RNA modifications and their immunogenic substrates. Finally, we highlight recent advances in cancer immunotherapies that aim to reduce cancer cell viability by targeting the enzymes responsible for RNA modifications. Our presentation of these dynamic research avenues sets the stage for future investigations in this field.

Probing the dynamic RNA structurome and its functions

RNA is a key regulator of almost every cellular process, and the structures adopted by RNA molecules are thought to be central to their functions. The recent fast-paced evolution of high-throughput sequencing-based RNA structure mapping methods has enabled the rapid in vivo structural interrogation of entire cellular transcriptomes. Collectively, these studies are shedding new light on the long underestimated complexity of the structural organization of the transcriptome - the RNA structurome. Moreover, recent analyses are challenging the view that the RNA structurome is a static entity by revealing how RNA molecules establish intricate networks of alternative intramolecular and intermolecular interactions and that these ensembles of RNA structures are dynamically regulated to finely tune RNA functions in living cells. This new understanding of how RNA can shape cell phenotypes has important implications for the development of RNA-targeted therapeutic strategies.

Multi-model predictive analysis of RNA solvent accessibility based on modified residual attention mechanism

Predicting RNA solvent accessibility using only primary sequence data can be regarded as sequence-based prediction work. Currently, the established studies for sequence-based RNA solvent accessibility prediction are limited due to the available number of datasets and black box prediction. To improve these issues, we first expanded the available RNA structures and then developed a sequence-based model using modified attention layers with different receptive fields to conform to the stem-loop structure of RNA chains. We measured the improvement with an extended dataset and further explored the model's interpretability by analysing the model structures, attention values and hyperparameters. Finally, we found that the developed model regarded the pieces of a sequence as templates during the training process. This work will be helpful for researchers who would like to build RNA attribute prediction models using deep learning in the future.

Predicting RNA solvent accessibility from multi-scale context feature via multi-shot neural network

Knowledge of RNA solvent accessibility has recently become attractive due to the increasing awareness of its importance for key biological process. Accurately predicting the solvent accessibility of RNA is crucial for understanding its 3D structure and biological function. In this study, we develop a novel computational method, termed M²pred, for accurately predicting the solvent accessibility of RNA from sequence-based multi-scale context feature. In M²pred, three single-view features, i.e., base-pairing probabilities, position-specific frequency matrix, and a binary one-hot encoding, are first generated as three feature sources, and immediately concatenated to engender a super feature. Secondly, for the super feature, the matrix-format features of each nucleotide are extracted using an initialized sliding window technique, and regularly stacked into a cube-format feature. Then, using multi-scale context feature extraction strategy, a pyramid feature constructed of contextual feature of four scales related to target nucleotides is extracted from the cube-format feature. Finally, a customized multi-shot neural network framework, which is equipped with four different scales of receptive fields mainly integrating several residual attention blocks, is designed to dig discrimination information from the contextual pyramid feature. Experimental results demonstrate that the proposed M²pred achieve a high prediction performance and outperforms existing state-of-the-art prediction methods of RNA solvent accessibility.

Probing RNA structures and functions by solvent accessibility: an overview from experimental and computational perspectives

Characterizing RNA structures and functions have mostly been focused on 2D, secondary and 3D, tertiary structures. Recent advances in experimental and computational techniques for probing or predicting RNA solvent accessibility make this 1D representation of tertiary structures an increasingly attractive feature to explore. Here, we provide a survey of these recent developments, which indicate the emergence of solvent accessibility as a simple 1D property, adding to secondary and tertiary structures for investigating complex structure–function relations of RNAs.

High-throughput mapping of RNA solvent accessibility at the single-nucleotide resolution by RtcB ligation between a fixed 5'-OH-end linker and unique 3'-P-end fragments from hydroxyl radical cleavage

Given the challenges for the experimental determination of RNA tertiary structures, probing solvent accessibility has become increasingly important to gain functional insights. Among various chemical probes developed, backbone-cleaving hydroxyl radical is the only one that can provide unbiased detection of all accessible nucleotides. However, the readouts have been based on reverse transcription (RT) stop at the cleaving sites, which are prone to false positives due to PCR amplification bias, early drop-off of reverse transcriptase, and the use of random primers in RT reaction. Here, we introduced a fixed-primer method called RL-Seq by performing RtcB Ligation (RL) between a fixed 5'-OH-end linker and unique 3'-P-end fragments from hydroxyl radical cleavage prior to high-throughput sequencing. The application of this method to E. coli ribosomes confirmed its ability to accurately probe solvent accessibility with high sensitivity (low required sequencing depth) and accuracy (strong correlation to structure-derived values) at the single-nucleotide resolution. Moreover, a near-perfect correlation was found between the experiments with and without using unique molecular identifiers, indicating negligible PCR biases in RL-Seq. Further improvement of RL-Seq and its potential transcriptome-wide applications are discussed.

"Bind, cleave and leave": multiple turnover catalysis of RNA cleavage by bulge-loop inducing supramolecular conjugates

Antisense sequence-specific knockdown of pathogenic RNA offers opportunities to find new solutions for therapeutic treatments. However, to gain a desired therapeutic effect, the multiple turnover catalysis is critical to inactivate many copies of emerging RNA sequences, which is difficult to achieve without sacrificing the sequence-specificity of cleavage. Here, engineering two or three catalytic peptides into the bulge-loop inducing molecular framework of antisense oligonucleotides achieved catalytic turnover of targeted RNA. Different supramolecular configurations revealed that cleavage of the RNA backbone upon sequence-specific hybridization with the catalyst accelerated with increase in the number of catalytic guanidinium groups, with almost complete demolition of target RNA in 24 h. Multiple sequence-specific cuts at different locations within and around the bulge-loop facilitated release of the catalyst for subsequent attacks of at least 10 further RNA substrate copies, such that delivery of only a few catalytic molecules could be sufficient to maintain knockdown of typical RNA copy numbers. We have developed fluorescent assay and kinetic simulation tools to characterise how the limited availability of different targets and catalysts had restrained catalytic reaction progress considerably, and to inform how to accelerate the catalytic destruction of shorter linear and larger RNAs even further.

Chemical and Enzymatic Probing of Viral RNAs: From Infancy to Maturity and Beyond

RNA molecules are key players in a variety of biological events, and this is particularly true for viral RNAs. To better understand the replication of those pathogens and try to block them, special attention has been paid to the structure of their RNAs. Methods to probe RNA structures have been developed since the 1960s; even if they have evolved over the years, they are still in use today and provide useful information on the folding of RNA molecules, including viral RNAs. The aim of this review is to offer a historical perspective on the structural probing methods used to decipher RNA structures before the development of the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) methodology and to show how they have influenced the current probing techniques. Actually, these technological breakthroughs, which involved advanced detection methods, were made possible thanks to the development of next-generation sequencing (NGS) but also to the previous works accumulated in the field of structural RNA biology. Finally, we will also discuss how high-throughput SHAPE (hSHAPE) paved the way for the development of sophisticated RNA structural techniques.

Single-sequence and profile-based prediction of RNA solvent accessibility using dilated convolutional neural network

MOTIVATION

RNA solvent accessibility, similar to protein solvent accessibility, reflects the structural regions that are accessible to solvents or other functional biomolecules, and plays an important role for structural and functional characterization. Unlike protein solvent accessibility, only a few tools are available for predicting RNA solvent accessibility despite the fact that millions of RNA transcripts have unknown structures and functions. Also, these tools have limited accuracy. Here, we have developed RNAsnap2 that uses a dilated convolutional neural network with a new feature, based on predicted base-pairing probabilities from LinearPartition.

RESULTS

Using the same training set from the recent predictor RNAsol, RNAsnap2 provides an 11% improvement in median Pearson Correlation Coefficient (PCC) and 9% improvement in mean absolute errors for the same test set of 45 RNA chains. A larger improvement (22% in median PCC) is observed for 31 newly deposited RNA chains that are non-redundant and independent from the training and the test sets. A single-sequence version of RNAsnap2 (i.e. without using sequence profiles generated from homology search by Infernal) has achieved comparable performance to the profile-based RNAsol. In addition, RNAsnap2 has achieved comparable performance for protein-bound and protein-free RNAs. Both RNAsnap2 and RNAsnap2 (SingleSeq) are expected to be useful for searching structural signatures and locating functional regions of non-coding RNAs.

AVAILABILITY AND IMPLEMENTATION

Standalone-versions of RNAsnap2 and RNAsnap2 (SingleSeq) are available at https://github.com/jaswindersingh2/RNAsnap2. Direct prediction can also be made at https://sparks-lab.org/server/rnasnap2. The datasets used in this research can also be downloaded from the GITHUB and the webserver mentioned above.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Assaying RNA structure with LASER-Seq

Chemical probing methods are crucial to our understanding of the structure and function of RNA molecules. The majority of chemical methods used to probe RNA structure report on Watson–Crick pairing, but tertiary structure parameters such as solvent accessibility can provide an additional layer of structural information, particularly in RNA-protein complexes. Herein we report the development of Light Activated Structural Examination of RNA by high-throughput sequencing, or LASER-Seq, for measuring RNA structure in cells with deep sequencing. LASER relies on a light-generated nicotinoyl nitrenium ion to form covalent adducts with the C8 position of adenosine and guanosine. Reactivity is governed by the accessibility of C8 to the light-generated probe. We compare structure probing by RT-stop and mutational profiling (MaP), demonstrating that LASER can be integrated with both platforms for RNA structure analyses. We find that LASER reactivity correlates with solvent accessibility across the entire ribosome, and that LASER can be used to rapidly survey for ligand binding sites in an unbiased fashion. LASER has a particular advantage in this last application, as it readily modifies paired nucleotides, enabling the identification of binding sites and conformational changes in highly structured RNA.

Excess primer degradation by Exo I improves the preparation of 3' cDNA ligation-based sequencing libraries

RNA sequencing library construction using single-stranded ligation of a DNA adapter to 3' ends of cDNAs often produces primer-adapter byproducts, which compete with cDNA-adapter ligation products during library amplification and, therefore, reduces the number of informative sequencing reads. We find that Escherichia coli Exo I digestion efficiently and selectively removes surplus reverse transcription primer and thereby reduces the primer-adapter product contamination in 3' cDNA ligation-based sequencing libraries, including small RNA libraries, which are typically similar in size to the primer-adapter products. We further demonstrate that Exo I treatment does not lead to trimming of the cDNA 3' end when duplexed with the RNA template. Exo I digestion is easy to perform and implement in other protocols and could facilitate a more widespread use of 3' cDNA ligation for sequencing-based applications.

Computational modeling of RNA 3D structure based on experimental data

RNA molecules are master regulators of cells. They are involved in a variety of molecular processes: they transmit genetic information, sense cellular signals and communicate responses, and even catalyze chemical reactions. As in the case of proteins, RNA function is dictated by its structure and by its ability to adopt different conformations, which in turn is encoded in the sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore the majority of known RNAs remain structurally uncharacterized. To address this problem, predictive computational methods were developed based on the accumulated knowledge of RNA structures determined so far, the physical basis of the RNA folding, and taking into account evolutionary considerations, such as conservation of functionally important motifs. However, all theoretical methods suffer from various limitations, and they are generally unable to accurately predict structures for RNA sequences longer than 100-nt residues unless aided by additional experimental data. In this article, we review experimental methods that can generate data usable by computational methods, as well as computational approaches for RNA structure prediction that can utilize data from experimental analyses. We outline methods and data types that can be potentially useful for RNA 3D structure modeling but are not commonly used by the existing software, suggesting directions for future development.

The evolution of RNA structural probing methods: From gels to next-generation sequencing

RNA molecules are important players in all domains of life and the study of the relationship between their multiple flexible states and the associated biological roles has increased in recent years. For several decades, chemical and enzymatic structural probing experiments have been used to determine RNA structure. During this time, there has been a steady improvement in probing reagents and experimental methods, and today the structural biologist community has a large range of tools at its disposal to probe the secondary structure of RNAs in vitro and in cells. Early experiments used radioactive labeling and polyacrylamide gel electrophoresis as read-out methods. This was superseded by capillary electrophoresis, and more recently by next-generation sequencing. Today, powerful structural probing methods can characterize RNA structure on a genome-wide scale. In this review, we will provide an overview of RNA structural probing methodologies from a historical and technical perspective. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Methods > RNA Analyses in vitro and In Silico RNA Methods > RNA Analyses in Cells.

Profiling the Nucleobase and Structure Selectivity of Anticancer Drugs and other DNA Alkylating Agents by RNA Sequencing

Drugs that covalently modify DNA are components of most chemotherapy regimens, often serving as first-line treatments. Classically, the reactivity and selectivity of DNA alkylating agents has been determined in vitro with short oligonucleotides. A statistically sound analysis of sequence preferences of alkylating agents is untenable with serial analysis methods because of the combinatorial explosion of sequence possibilities. Next-generation sequencing (NGS) is ideally suited for the broad characterization of sequence or structure selectivities because it analyzes many sequences at once. Herein, NGS is used to report on the chemoselectivity of alkylating agents on RNA and this technology is applied to the previously uncharacterized alkylating agent trimethylsilyl diazomethane.

Time-Resolved Hydroxyl Radical Footprinting of RNA with X-Rays

RNA footprinting by hydroxyl radical cleavage provides 'snapshots' of RNA tertiary structure or protein interactions that bury the RNA backbone. Generation of hydroxyl radicals with a high-flux synchrotron X-ray beam provides analysis on a short timescale (5-100 msec), which enables the structures of folding intermediates or other transient conformational states to be determined in biochemical solutions or cells. This article provides protocols for using synchrotron beamlines for hydroxyl radical footprinting. © 2018 by John Wiley & Sons, Inc.

Directed hydroxyl radical probing reveals Upf1 binding to the 80S ribosomal E site rRNA at the L1 stalk

Upf1 is an SF1-family RNA helicase that is essential for the nonsense-mediated decay (NMD) process in eukaryotes. While Upf1 has been shown to interact with 80S ribosomes, the molecular details of this interaction were unknown. Using purified recombinant proteins and high-throughput sequencing combined with Fe-BABE directed hydroxyl radical probing (HTS-BABE) we have characterized the interaction between Upf1 and the yeast 80S ribosome. We identify the 1C domain of Upf1, an alpha-helical insertion in the RecA helicase core, to be essential for ribosome binding, and determine that the L1 stalk of 25S rRNA is the binding site for Upf1 on the ribosome. Using the cleavage sites identified by hydroxyl radical probing and high-resolution structures of both yeast Upf1 and the human 80S ribosome, we provide a model of a Upf1:80S structure. Our model requires that the L1 stalk adopt an open configuration as adopted by an un-rotated, or classical-state, ribosome. Our results shed light on the interaction between Upf1 and the ribosome, and suggest that Upf1 may specifically engage a classical-state ribosome during translation.

Structure-seq2: sensitive and accurate genome-wide profiling of RNA structure in vivo

RNA serves many functions in biology such as splicing, temperature sensing, and innate immunity. These functions are often determined by the structure of RNA. There is thus a pressing need to understand RNA structure and how it changes during diverse biological processes both in vivo and genome-wide. Here, we present Structure-seq2, which provides nucleotide-resolution RNA structural information in vivo and genome-wide. This optimized version of our original Structure-seq method increases sensitivity by at least 4-fold and improves data quality by minimizing formation of a deleterious by-product, reducing ligation bias, and improving read coverage. We also present a variation of Structure-seq2 in which a biotinylated nucleotide is incorporated during reverse transcription, which greatly facilitates the protocol by eliminating two PAGE purification steps. We benchmark Structure-seq2 on both mRNA and rRNA structure in rice (Oryza sativa). We demonstrate that Structure-seq2 can lead to new biological insights. Our Structure-seq2 datasets uncover hidden breaks in chloroplast rRNA and identify a previously unreported N¹-methyladenosine (m¹A) in a nuclear-encoded Oryza sativa rRNA. Overall, Structure-seq2 is a rapid, sensitive, and unbiased method to probe RNA in vivo and genome-wide that facilitates new insights into RNA biology.

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.

Managing the sequence-specificity of antisense oligonucleotides in drug discovery

All drugs perturb the expression of many genes in the cells that are exposed to them. These gene expression changes can be divided into effects resulting from engaging the intended target and effects resulting from engaging unintended targets. For antisense oligonucleotides, developments in bioinformatics algorithms, and the quality of sequence databases, allow oligonucleotide sequences to be analyzed computationally, in terms of the predictability of their interactions with intended and unintended RNA targets. Applying these tools enables selection of sequence-specific oligonucleotides where no- or only few unintended RNA targets are expected. To evaluate oligonucleotide sequence-specificity experimentally, we recommend a transcriptomics protocol where two or more oligonucleotides targeting the same RNA molecule, but with entirely different sequences, are evaluated together. This helps to clarify which changes in cellular RNA levels result from downstream processes of engaging the intended target, and which are likely to be related to engaging unintended targets. As required for all classes of drugs, the toxic potential of oligonucleotides must be evaluated in cell- and animal models before clinical testing. Since potential adverse effects related to unintended targeting are sequence-dependent and therefore species-specific, in vitro toxicology assays in human cells are especially relevant in oligonucleotide drug discovery.

Dawn of the in vivo RNA structurome and interactome

RNA is one of the most fascinating biomolecules in living systems given its structural versatility to fold into elaborate architectures for important biological functions such as gene regulation, catalysis, and information storage. Knowledge of RNA structures and interactions can provide deep insights into their functional roles in vivo For decades, RNA structural studies have been conducted on a transcript-by-transcript basis. The advent of next-generation sequencing (NGS) has enabled the development of transcriptome-wide structural probing methods to profile the global landscape of RNA structures and interactions, also known as the RNA structurome and interactome, which transformed our understanding of the RNA structure-function relationship on a transcriptomic scale. In this review, molecular tools and NGS methods used for RNA structure probing are presented, novel insights uncovered by RNA structurome and interactome studies are highlighted, and perspectives on current challenges and potential future directions are discussed. A more complete understanding of the RNA structures and interactions in vivo will help illuminate the novel roles of RNA in gene regulation, development, and diseases.

Structural transitions during large ribosomal subunit maturation analyzed by tethered nuclease structure probing in S. cerevisiae

Yeast large ribosomal subunit (LSU) precursors are subject to substantial changes in protein composition during their maturation due to coordinated transient interactions with a large number of ribosome biogenesis factors and due to the assembly of ribosomal proteins. These compositional changes go along with stepwise processing of LSU rRNA precursors and with specific rRNA folding events, as revealed by recent cryo-electron microscopy analyses of late nuclear and cytoplasmic LSU precursors. Here we aimed to analyze changes in the spatial rRNA surrounding of selected ribosomal proteins during yeast LSU maturation. For this we combined a recently developed tethered tertiary structure probing approach with both targeted and high throughput readout strategies. Several structural features of late LSU precursors were faithfully detected by this procedure. In addition, the obtained data let us suggest that early rRNA precursor processing events are accompanied by a global transition from a flexible to a spatially restricted rRNA conformation. For intermediate LSU precursors a number of structural hallmarks could be addressed which include the fold of the internal transcribed spacer between 5.8S rRNA and 25S rRNA, the orientation of the central protuberance and the spatial organization of the interface between LSU rRNA domains I and III.

RNA structure through multidimensional chemical mapping

The discoveries of myriad non-coding RNA molecules, each transiting through multiple flexible states in cells or virions, present major challenges for structure determination. Advances in high-throughput chemical mapping give new routes for characterizing entire transcriptomes in vivo, but the resulting one-dimensional data generally remain too information-poor to allow accurate de novo structure determination. Multidimensional chemical mapping (MCM) methods seek to address this challenge. Mutate-and-map (M2), RNA interaction groups by mutational profiling (RING-MaP and MaP-2D analysis) and multiplexed •OH cleavage analysis (MOHCA) measure how the chemical reactivities of every nucleotide in an RNA molecule change in response to modifications at every other nucleotide. A growing body of in vitro blind tests and compensatory mutation/rescue experiments indicate that MCM methods give consistently accurate secondary structures and global tertiary structures for ribozymes, ribosomal domains and ligand-bound riboswitch aptamers up to 200 nucleotides in length. Importantly, MCM analyses provide detailed information on structurally heterogeneous RNA states, such as ligand-free riboswitches that are functionally important but difficult to resolve with other approaches. The sequencing requirements of currently available MCM protocols scale at least quadratically with RNA length, precluding general application to transcriptomes or viral genomes at present. We propose a modify-cross-link-map (MXM) expansion to overcome this and other current limitations to resolving the in vivo 'RNA structurome'.

Comparative and integrative analysis of RNA structural profiling data: current practices and emerging questions

BACKGROUND

Structure profiling experiments provide single-nucleotide information on RNA structure. Recent advances in chemistry combined with application of high-throughput sequencing have enabled structure profiling at transcriptome scale and in living cells, creating unprecedented opportunities for RNA biology. Propelled by these experimental advances, massive data with ever-increasing diversity and complexity have been generated, which give rise to new challenges in interpreting and analyzing these data.

RESULTS

We review current practices in analysis of structure profiling data with emphasis on comparative and integrative analysis as well as highlight emerging questions. Comparative analysis has revealed structural patterns across transcriptomes and has become an integral component of recent profiling studies. Additionally, profiling data can be integrated into traditional structure prediction algorithms to improve prediction accuracy.

CONCLUSIONS

To keep pace with experimental developments, methods to facilitate, enhance and refine such analyses are needed. Parallel advances in analysis methodology will complement profiling technologies and help them reach their full potential.

The RNA Epistructurome: Uncovering RNA Function by Studying Structure and Post-Transcriptional Modifications

A large fraction of higher metazoan genomes transcribe RNA molecules whose functions extend far beyond carrying instructions for protein synthesis. Although RNA is apparently a simple molecule, the ways in which it performs many of its functions have remained highly elusive for decades. As learned from studying ribosomal and transfer RNAs, two of the key features influencing the function of RNA are its structure and post-transcriptional modifications. A deep understanding of RNA function therefore requires rapid and straightforward approaches to study the complex and intricate landscape of RNA structures and modifications. In this review we summarize and discuss the most recent methods and findings in the field of RNA biology, with an eye toward new frontiers and open questions.

Metrics for rapid quality control in RNA structure probing experiments

MOTIVATION

The diverse functionalities of RNA can be attributed to its capacity to form complex and varied structures. The recent proliferation of new structure probing techniques coupled with high-throughput sequencing has helped RNA studies expand in both scope and depth. Despite differences in techniques, most experiments face similar challenges in reproducibility due to the stochastic nature of chemical probing and sequencing. As these protocols expand to transcriptome-wide studies, quality control becomes a more daunting task. General and efficient methodologies are needed to quantify variability and quality in the wide range of current and emerging structure probing experiments.

RESULTS

We develop metrics to rapidly and quantitatively evaluate data quality from structure probing experiments, demonstrating their efficacy on both small synthetic libraries and transcriptome-wide datasets. We use a signal-to-noise ratio concept to evaluate replicate agreement, which has the capacity to identify high-quality data. We also consider and compare two methods to assess variability inherent in probing experiments, which we then utilize to evaluate the coverage adjustments needed to meet desired quality. The developed metrics and tools will be useful in summarizing large-scale datasets and will help standardize quality control in the field.

AVAILABILITY AND IMPLEMENTATION

The data and methods used in this article are freely available at: http://bme.ucdavis.edu/aviranlab/SPEQC_software CONTACT: saviran@ucdavis.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Data-directed RNA secondary structure prediction using probabilistic modeling

Structure dictates the function of many RNAs, but secondary RNA structure analysis is either labor intensive and costly or relies on computational predictions that are often inaccurate. These limitations are alleviated by integration of structure probing data into prediction algorithms. However, existing algorithms are optimized for a specific type of probing data. Recently, new chemistries combined with advances in sequencing have facilitated structure probing at unprecedented scale and sensitivity. These novel technologies and anticipated wealth of data highlight a need for algorithms that readily accommodate more complex and diverse input sources. We implemented and investigated a recently outlined probabilistic framework for RNA secondary structure prediction and extended it to accommodate further refinement of structural information. This framework utilizes direct likelihood-based calculations of pseudo-energy terms per considered structural context and can readily accommodate diverse data types and complex data dependencies. We use real data in conjunction with simulations to evaluate performances of several implementations and to show that proper integration of structural contexts can lead to improvements. Our tests also reveal discrepancies between real data and simulations, which we show can be alleviated by refined modeling. We then propose statistical preprocessing approaches to standardize data interpretation and integration into such a generic framework. We further systematically quantify the information content of data subsets, demonstrating that high reactivities are major drivers of SHAPE-directed predictions and that better understanding of less informative reactivities is key to further improvements. Finally, we provide evidence for the adaptive capability of our framework using mock probe simulations.

Probing the structure of ribosome assembly intermediates in vivo using DMS and hydroxyl radical footprinting

The assembly of the Escherichia coli ribosome has been widely studied and characterized in vitro. Despite this, ribosome biogenesis in living cells is only partly understood because assembly is coupled with transcription, modification and processing of the pre-ribosomal RNA. We present a method for footprinting and isolating pre-rRNA as it is synthesized in E. coli cells. Pre-rRNA synthesis is synchronized by starvation, followed by nutrient upshift. RNA synthesized during outgrowth is metabolically labeled to facilitate isolation of recent transcripts. Combining this technique with two in vivo RNA probing methods, hydroxyl radical and DMS footprinting, allows the structure of nascent RNA to be probed over time. Together, these can be used to determine changes in the structures of ribosome assembly intermediates as they fold in vivo.

Predicting RNA secondary structures from sequence and probing data

RNA secondary structures have proven essential for understanding the regulatory functions performed by RNA such as microRNAs, bacterial small RNAs, or riboswitches. This success is in part due to the availability of efficient computational methods for predicting RNA secondary structures. Recent advances focus on dealing with the inherent uncertainty of prediction by considering the ensemble of possible structures rather than the single most stable one. Moreover, the advent of high-throughput structural probing has spurred the development of computational methods that incorporate such experimental data as auxiliary information.

High-Throughput Nuclease Probing of RNA Structures Using FragSeq

High-throughput sequencing of cDNA (RNA-Seq) can be used to generate nuclease accessibility data for many distinct transcripts in the same mixture simultaneously. Such assays accelerate RNA structure analysis and provide researchers with new technologies to tackle biological questions on a transcriptome-wide scale. FragSeq is an experimental assay for transcriptome-wide RNA structure probing using RNA-Seq, coupled with data analysis tools that allow quantitative determination of nuclease accessibility at single-base resolution. We provide a practical guide to designing and carrying out FragSeq experiments and data analysis.

Probing Xist RNA Structure in Cells Using Targeted Structure-Seq

The long non-coding RNA (lncRNA) Xist is a master regulator of X-chromosome inactivation in mammalian cells. Models for how Xist and other lncRNAs function depend on thermodynamically stable secondary and higher-order structures that RNAs can form in the context of a cell. Probing accessible RNA bases can provide data to build models of RNA conformation that provide insight into RNA function, molecular evolution, and modularity. To study the structure of Xist in cells, we built upon recent advances in RNA secondary structure mapping and modeling to develop Targeted Structure-Seq, which combines chemical probing of RNA structure in cells with target-specific massively parallel sequencing. By enriching for signals from the RNA of interest, Targeted Structure-Seq achieves high coverage of the target RNA with relatively few sequencing reads, thus providing a targeted and scalable approach to analyze RNA conformation in cells. We use this approach to probe the full-length Xist lncRNA to develop new models for functional elements within Xist, including the repeat A element in the 5’-end of Xist. This analysis also identified new structural elements in Xist that are evolutionarily conserved, including a new element proximal to the C repeats that is important for Xist function.

To do their jobs, many RNAs need to fold into structures (through base-paring). We were interested in the conformation of a specific mammalian RNA, Xist, when it is inside a cell. Xist is a very large non-coding RNA (lncRNA), that is >17,000 nt long. Xist is particularly important because it is one of the first lncRNAs to be discovered, and turns genes off across an entire chromosome. To figure out how Xist RNA is folded in mouse cells, we developed a new approach, Targeted Structure-Seq, to examine the conformation of large RNAs like Xist. Using computer modeling, we identified parts of Xist that are base paired into RNA duplexes. We also determined which parts of the Xist RNA are likely to be structured. This work provides a new tool for studying the secondary structure of any large RNA, and helps us understand what the important pieces of Xist look like while Xist does its work in the cell.

Understanding DNA under oxidative stress and sensitization: the role of molecular modeling

DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.

The potential of the riboSNitch in personalized medicine

RNA conformation plays a significant role in stability, ligand binding, transcription, and translation. Single nucleotide variants (SNVs) have the potential to disrupt specific structural elements because RNA folds in a sequence-specific manner. A riboSNitch is an element of RNA structure with a specific function that is disrupted by an SNV or a single nucleotide polymorphism (SNP; or polymorphism; SNVs occur with low frequency in the population, <1%). The riboSNitch is analogous to a riboswitch, where binding of a small molecule rather than mutation alters the structure of the RNA to control gene regulation. RiboSNitches are particularly relevant to interpreting the results of genome-wide association studies (GWAS). Often GWAS identify SNPs associated with a phenotype mapping to noncoding regions of the genome. Because a majority of the human genome is transcribed, significant subsets of GWAS SNPs are putative riboSNitches. The extent to which the transcriptome is tolerant of SNP-induced structure change is still poorly understood. Recent advances in ultra high-throughput structure probing begin to reveal the structural complexities of mutation-induced structure change. This review summarizes our current understanding of SNV and SNP-induced structure change in the human transcriptome and discusses the importance of riboSNitch discovery in interpreting GWAS results and massive sequencing projects.

Current perspectives on RNA secondary structure probing

The range of roles played by structured RNAs in biological systems is vast. At the same time as we are learning more about the importance of RNA structure, recent advances in reagents, methods and technology mean that RNA secondary structural probing has become faster and more accurate. As a result, the capabilities of laboratories that already perform this type of structural analysis have increased greatly, and it has also become more widely accessible. The present review summarizes established and recently developed techniques. The information we can derive from secondary structural analysis is assessed, together with the areas in which we are likely to see exciting developments in the near future.

SHAPE Selection (SHAPES) enrich for RNA structure signal in SHAPE sequencing-based probing data

Selective 2' Hydroxyl Acylation analyzed by Primer Extension (SHAPE) is an accurate method for probing of RNA secondary structure. In existing SHAPE methods, the SHAPE probing signal is normalized to a no-reagent control to correct for the background caused by premature termination of the reverse transcriptase. Here, we introduce a SHAPE Selection (SHAPES) reagent, N-propanone isatoic anhydride (NPIA), which retains the ability of SHAPE reagents to accurately probe RNA structure, but also allows covalent coupling between the SHAPES reagent and a biotin molecule. We demonstrate that SHAPES-based selection of cDNA-RNA hybrids on streptavidin beads effectively removes the large majority of background signal present in SHAPE probing data and that sequencing-based SHAPES data contain the same amount of RNA structure data as regular sequencing-based SHAPE data obtained through normalization to a no-reagent control. Moreover, the selection efficiently enriches for probed RNAs, suggesting that the SHAPES strategy will be useful for applications with high-background and low-probing signal such as in vivo RNA structure probing.

Reproducible Analysis of Sequencing-Based RNA Structure Probing Data with User-Friendly Tools

RNA structure-probing data can improve the prediction of RNA secondary and tertiary structure and allow structural changes to be identified and investigated. In recent years, massive parallel sequencing has dramatically improved the throughput of RNA structure probing experiments, but at the same time also made analysis of the data challenging for scientists without formal training in computational biology. Here, we discuss different strategies for data analysis of massive parallel sequencing-based structure-probing data. To facilitate reproducible and standardized analysis of this type of data, we have made a collection of tools, which allow raw sequencing reads to be converted to normalized probing values using different published strategies. In addition, we also provide tools for visualization of the probing data in the UCSC Genome Browser and for converting RNA coordinates to genomic coordinates and vice versa. The collection is implemented as functions in the R statistical environment and as tools in the Galaxy platform, making them easily accessible for the scientific community. We demonstrate the usefulness of the collection by applying it to the analysis of sequencing-based hydroxyl radical probing data and comparing different normalization strategies.

Probing small non-coding RNAs structures

The diverse roles of RNAs depend on their ability to fold so as to form biologically functional structures. Thus, understanding the function of a given RNA molecule often requires experimental analysis of its secondary structure by in vitro RNA probing, which is more accurate than using prediction programs only. This chapter presents in vitro RNA probing protocols that we routinely use, from RNA transcript production and purification to RNA structure determination using enzymatic (RNases T1, T2, and V1) and chemical (DMS, CMCT, kethoxal, and Pb(2+)) probing performed on both unlabeled and end-labeled RNAs.

Rational experiment design for sequencing-based RNA structure mapping

Structure mapping is a classic experimental approach for determining nucleic acid structure that has gained renewed interest in recent years following advances in chemistry, genomics, and informatics. The approach encompasses numerous techniques that use different means to introduce nucleotide-level modifications in a structure-dependent manner. Modifications are assayed via cDNA fragment analysis, using electrophoresis or next-generation sequencing (NGS). The recent advent of NGS has dramatically increased the throughput, multiplexing capacity, and scope of RNA structure mapping assays, thereby opening new possibilities for genome-scale, de novo, and in vivo studies. From an informatics standpoint, NGS is more informative than prior technologies by virtue of delivering direct molecular measurements in the form of digital sequence counts. Motivated by these new capabilities, we introduce a novel model-based in silico approach for quantitative design of large-scale multiplexed NGS structure mapping assays, which takes advantage of the direct and digital nature of NGS readouts. We use it to characterize the relationship between controllable experimental parameters and the precision of mapping measurements. Our results highlight the complexity of these dependencies and shed light on relevant tradeoffs and pitfalls, which can be difficult to discern by intuition alone. We demonstrate our approach by quantitatively assessing the robustness of SHAPE-Seq measurements, obtained by multiplexing SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry in conjunction with NGS. We then utilize it to elucidate design considerations in advanced genome-wide approaches for probing the transcriptome, which recently obtained in vivo information using dimethyl sulfate (DMS) chemistry.

Chemical probing of RNA with the hydroxyl radical at single-atom resolution

While hydroxyl radical cleavage is widely used to map RNA tertiary structure, lack of mechanistic understanding of strand break formation limits the degree of structural insight that can be obtained from this experiment. Here, we determine how individual ribose hydrogens of sarcin/ricin loop RNA participate in strand cleavage. We find that substituting deuterium for hydrogen at a ribose 5'-carbon produces a kinetic isotope effect on cleavage; the major cleavage product is an RNA strand terminated by a 5'-aldehyde. We conclude that hydroxyl radical abstracts a 5'-hydrogen atom, leading to RNA strand cleavage. We used this approach to obtain structural information for a GUA base triple, a common tertiary structural feature of RNA. Cleavage at U exhibits a large 5' deuterium kinetic isotope effect, a potential signature of a base triple. Others had noted a ribose-phosphate hydrogen bond involving the G 2'-OH and the U phosphate of the GUA triple, and suggested that this hydrogen bond contributes to backbone rigidity. Substituting deoxyguanosine for G, to eliminate this hydrogen bond, results in a substantial decrease in cleavage at G and U of the triple. We conclude that this hydrogen bond is a linchpin of backbone structure around the triple.

Insights into RNA structure and function from genome-wide studies

A comprehensive understanding of RNA structure will provide fundamental insights into the cellular function of both coding and non-coding RNAs. Although many RNA structures have been analysed by traditional biophysical and biochemical methods, the low-throughput nature of these approaches has prevented investigation of the vast majority of cellular transcripts. Triggered by advances in sequencing technology, genome-wide approaches for probing the transcriptome are beginning to reveal how RNA structure affects each step of protein expression and RNA stability. In this Review, we discuss the emerging relationships between RNA structure and the regulation of gene expression.