Global analysis of RNA secondary structure in two metazoans

Probing RNA structural landscapes across Candida yeast genomes

Understanding the intricate roles of RNA molecules in virulence and host-pathogen interactions can provide valuable insights into combatting infections and improving human health. Although much progress has been achieved in understanding transcriptional regulation during host-pathogen interactions in diverse species, more is needed to know about the structure of pathogen RNAs. This is particularly true for fungal pathogens, including pathogenic yeasts of the Candida genus, which are the leading cause of hospital-acquired fungal infections. Our work addresses the gap between RNA structure and their biology by employing genome-wide structure probing to comprehensively explore the structural landscape of mRNAs and long non-coding RNAs (lncRNAs) in the four major Candida pathogens. Specifically focusing on mRNA, we observe a robust correlation between sequence conservation and structural characteristics in orthologous transcripts, significantly when sequence identity exceeds 50%, highlighting structural feature conservation among closely related species. We investigate the impact of single nucleotide polymorphisms (SNPs) on mRNA secondary structure. SNPs within 5' untranslated regions (UTRs) tend to occur in less structured positions, suggesting structural constraints influencing transcript regulation. Furthermore, we compare the structural properties of coding regions and UTRs, noting that coding regions are generally more structured than UTRs, consistent with similar trends in other species. Additionally, we provide the first experimental characterization of lncRNA structures in Candida species. Most lncRNAs form independent subdomains, similar to human lncRNAs. Notably, we identify hairpin-like structures in lncRNAs, a feature known to be functionally significant. Comparing hairpin prevalence between lncRNAs and protein-coding genes, we find enrichment in lncRNAs across Candida species, humans, and Arabidopsis thaliana, suggesting a conserved role for these structures. In summary, our study offers valuable insights into the interplay between RNA sequence, structure, and function in Candida pathogens, with implications for gene expression regulation and potential therapeutic strategies against Candida infections.

Tailor made: the art of therapeutic mRNA design

mRNA medicine is a new and rapidly developing field in which the delivery of genetic information in the form of mRNA is used to direct therapeutic protein production in humans. This approach, which allows for the quick and efficient identification and optimization of drug candidates for both large populations and individual patients, has the potential to revolutionize the way we prevent and treat disease. A key feature of mRNA medicines is their high degree of designability, although the design choices involved are complex. Maximizing the production of therapeutic proteins from mRNA medicines requires a thorough understanding of how nucleotide sequence, nucleotide modification and RNA structure interplay to affect translational efficiency and mRNA stability. In this Review, we describe the principles that underlie the physical stability and biological activity of mRNA and emphasize their relevance to the myriad considerations that factor into therapeutic mRNA design.

PELOTA and HBS1 suppress co-translational messenger RNA decay in Arabidopsis

Various messenger RNA (mRNA) decay mechanisms play major roles in controlling mRNA quality and quantity in eukaryotic organisms under different conditions. While it is known that the recently discovered co-translational mRNA decay (CTRD), the mechanism that allows mRNAs to be degraded while still being actively translated, is prevalent in yeast, humans, and various angiosperms, the regulation of this decay mechanism is less well studied. Moreover, it is still unclear whether this decay mechanism plays any role in the regulation of specific physiological processes in eukaryotes. Here, by re-analyzing the publicly available polysome profiling or ribosome footprinting and degradome sequencing datasets, we discovered that highly translated mRNAs tend to have lower co-translational decay levels. Based on this finding, we then identified Pelota and Hbs1, the translation-related ribosome rescue factors, as suppressors of co-translational mRNA decay in Arabidopsis. Furthermore, we found that Pelota and Hbs1 null mutants have lower germination rates compared to the wild-type plants, implying that proper regulation of co-translational mRNA decay is essential for normal developmental processes. In total, our study provides further insights into the regulation of CTRD in Arabidopsis and demonstrates that this decay mechanism does play important roles in Arabidopsis physiological processes.

Pervasive downstream RNA hairpins dynamically dictate start-codon selection

Translational reprogramming allows organisms to adapt to changing conditions. Upstream start codons (uAUGs), which are prevalently present in mRNAs, have crucial roles in regulating translation by providing alternative translation start sites1-4. However, what determines this selective initiation of translation between conditions remains unclear. Here, by integrating transcriptome-wide translational and structural analyses during pattern-triggered immunity in Arabidopsis, we found that transcripts with immune-induced translation are enriched with upstream open reading frames (uORFs). Without infection, these uORFs are selectively translated owing to hairpins immediately downstream of uAUGs, presumably by slowing and engaging the scanning preinitiation complex. Modelling using deep learning provides unbiased support for these recognizable double-stranded RNA structures downstream of uAUGs (which we term uAUG-ds) being responsible for the selective translation of uAUGs, and allows the prediction and rational design of translating uAUG-ds. We found that uAUG-ds-mediated regulation can be generalized to human cells. Moreover, uAUG-ds-mediated start-codon selection is dynamically regulated. After immune challenge in plants, induced RNA helicases that are homologous to Ded1p in yeast and DDX3X in humans resolve these structures, allowing ribosomes to bypass uAUGs to translate downstream defence proteins. This study shows that mRNA structures dynamically regulate start-codon selection. The prevalence of this RNA structural feature and the conservation of RNA helicases across kingdoms suggest that mRNA structural remodelling is a general feature of translational reprogramming.

Recent advances in RNA structurome

RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

An essential role for the piRNA pathway in regulating the ribosomal RNA pool in C. elegans

Piwi-interacting RNAs (piRNAs) are RNA effectors with key roles in maintaining genome integrity and promoting fertility in metazoans. In Caenorhabditis elegans loss of piRNAs leads to a transgenerational sterility phenotype. The plethora of piRNAs and their ability to silence transcripts with imperfect complementarity have raised several (non-exclusive) models for the underlying drivers of sterility. Here, we report the extranuclear and transferable nature of the sterility driver, its suppression via mutations disrupting the endogenous RNAi and poly-uridylation machinery, and copy-number amplification at the ribosomal DNA locus. In piRNA-deficient animals, several small interfering RNA (siRNA) populations become increasingly overabundant in the generations preceding loss of germline function, including ribosomal siRNAs (risiRNAs). A concomitant increase in uridylated sense rRNA fragments suggests that poly-uridylation may potentiate RNAi-mediated gene silencing of rRNAs. We conclude that loss of the piRNA machinery allows for unchecked amplification of siRNA populations, originating from abundant highly structured RNAs, to deleterious levels.

The Expression Profiles of mRNAs and lncRNAs in Buffalo Muscle Stem Cells Driving Myogenic Differentiation

Buffalo breeding has become an important branch of the beef cattle industry. Hence, it is of great significance to study buffalo meat production and meat quality. However, the expression profiles of mRNA and long non-coding RNAs (lncRNA) molecules in muscle stem cells (MuSCs) development in buffalo have not been explored fully. We, therefore, performed mRNA and lncRNA expression profiling analysis during the proliferation and differentiation phases of MuSCs in buffalo. The results showed that there were 4,820 differentially expressed genes as well as 12,227 mRNAs and 1,352 lncRNAs. These genes were shown to be enriched in essential biological processes such as cell cycle, p53 signaling pathway, RNA transport and calcium signaling pathway. We also identified a number of functionally important genes, such as MCMC4, SERDINE1, ISLR, LOC102394806, and LOC102403551, and found that interference with MYLPF expression significantly inhibited the differentiation of MuSCs. In conclusion, our research revealed the characteristics of mRNA and lncRNA expression during the differentiation of buffalo MuSCs. This study can be used as an important reference for the study of RNA regulation during muscle development in buffalo.

RNA structure probing uncovers RNA structure-dependent biological functions

RNA molecules fold into complex structures that enable their diverse functions in cells. Recent revolutionary innovations in transcriptome-wide RNA structural probing of living cells have ushered in a new era in understanding RNA functions. Here, we summarize the latest technological advances for probing RNA secondary structures and discuss striking discoveries that have linked RNA regulation and biological processes through interrogation of RNA structures. In particular, we highlight how different long noncoding RNAs form into distinct secondary structures that determine their modes of interactions with protein partners to realize their unique functions. These dynamic structures mediate RNA regulatory functions through altering interactions with proteins and other RNAs. We also outline current methodological hurdles and speculate about future directions for development of the next generation of RNA structure-probing technologies of higher sensitivity and resolution, which could then be applied in increasingly physiologically relevant studies.

The RNA structurome in the asexual blood stages of malaria pathogen plasmodium falciparum

Plasmodium falciparum is a deadly human pathogen responsible for the devastating disease called malaria. In this study, we measured the differential accumulation of RNA secondary structures in coding and non-coding transcripts from the asexual developmental cycle in P. falciparum in human red blood cells. Our comprehensive analysis that combined high-throughput nuclease mapping of RNA structures by duplex RNA-seq, SHAPE-directed RNA structure validation, immunoaffinity purification and characterization of antisense RNAs collectively measured differentially base-paired RNA regions throughout the parasite's asexual RBC cycle. Our mapping data not only aligned to a diverse pool of RNAs with known structures but also enabled us to identify new structural RNA regions in the malaria genome. On average, approximately 71% of the genes with secondary structures are found to be protein coding mRNAs. The mapping pattern of these base-paired RNAs corresponded to all regions of mRNAs, including the 5' UTR, CDS and 3' UTR as well as the start and stop codons. Histone family genes which are known to form secondary structures in their mRNAs and transcripts from genes which are important for transcriptional and post-transcriptional control, such as the unique plant-like transcription factor family, ApiAP2, DNA-/RNA-binding protein, Alba3 and proteins important for RBC invasion and malaria cytoadherence also showed strong accumulation of duplex RNA reads in various asexual stages in P. falciparum. Intriguingly, our study determined stage-specific, dynamic relationships between mRNA structural contents and translation efficiency in P. falciparum asexual blood stages, suggesting an essential role of RNA structural changes in malaria gene expression programs.

Computational search of hybrid human/SARS-CoV-2 dsRNA reveals unique viral sequences that diverge from those of other coronavirus strains

The role of the RNAi/Dicer/Ago system in degrading RNA viruses has been elusive in mammals in the past, which has prompted authors to think that interferon (IFN) synthesis is essential in this clade, relegating the RNAi defense strategy against viral infection as an accessory function. However, recent publications highlight the existence of abundant viral small interference and micro RNAs (VsiRNAs and VmiRNAs) in both cell-line and whole organism based experiments, indicating a contribution of these molecules in host responses and/or viral replication. We explore the theoretical possibility that RNAi triggered by SARS-CoV-2 might degrade some host transcripts in the opposite direction, although this hypothesis seems counterintuitive. The SARS-CoV-2 genome was therefore computationally searched for exact intrapairing within the viral RNA and exact hybrid pairing with the human transcriptome over a minimum of 20 bases in length. Minimal segments of 20-base lengths of SARS-CoV-2 RNA were found based on the theoretical matching with existing complementary strands in the human host transcriptome. Few human genes potentially annealing with SARS-CoV-2 RNA, including mitochondrial deubiquitinase USP30, the subunit of ubiquitin protein ligase complex FBXO21 and two long noncoding RNAs, were retrieved. The hypothesis that viral-originated RNAi might mediate degradation of host transcriptome messages was corroborated by published high throughput sequencing of RNA from infected tissues and cultured cells, clinical observation and phylogenetic comparative analysis, indicating a strong specificity of these SARS-CoV-2 hybrid pairing sequences for human genomes.

RNA architecture influences plant biology

The majority of the genome is transcribed to RNA in living organisms. RNA transcripts can form astonishing arrays of secondary and tertiary structures via Watson-Crick, Hoogsteen, or wobble base pairing. In vivo, RNA folding is not a simple thermodynamic event of minimizing free energy. Instead, the process is constrained by transcription, RNA-binding proteins, steric factors, and the microenvironment. RNA secondary structure (RSS) plays myriad roles in numerous biological processes, such as RNA processing, stability, transportation, and translation in prokaryotes and eukaryotes. Emerging evidence has also implicated RSS in RNA trafficking, liquid-liquid phase separation, and plant responses to environmental variations such as temperature and salinity. At molecular level, RSS is correlated with splicing, polyadenylation, protein synthesis, and miRNA biogenesis and functions. In this review, we summarize newly reported methods for probing RSS in vivo and functions and mechanisms of RSS in plant physiology.

In Vivo and In Vitro Genome-Wide Profiling of RNA Secondary Structures Reveals Key Regulatory Features in Plasmodium falciparum

It is widely accepted that the structure of RNA plays important roles in a number of biological processes, such as polyadenylation, splicing, and catalytic functions. Dynamic changes in RNA structure are able to regulate the gene expression programme and can be used as a highly specific and subtle mechanism for governing cellular processes. However, the nature of most RNA secondary structures in Plasmodium falciparum has not been determined. To investigate the genome-wide RNA secondary structural features at single-nucleotide resolution in P. falciparum, we applied a novel high-throughput method utilizing the chemical modification of RNA structures to characterize these structures. Structural data from parasites are in close agreement with the known 18S ribosomal RNA secondary structures of P. falciparum and can help to predict the in vivo RNA secondary structure of a total of 3,396 transcripts in the ring-stage and trophozoite-stage developmental cycles. By parallel analysis of RNA structures in vivo and in vitro during the Plasmodium parasite ring-stage and trophozoite-stage intraerythrocytic developmental cycles, we identified some key regulatory features. Recent studies have established that the RNA structure is a ubiquitous and fundamental regulator of gene expression. Our study indicate that there is a critical connection between RNA secondary structure and mRNA abundance during the complex biological programme of P. falciparum. This work presents a useful framework and important results, which may facilitate further research investigating the interactions between RNA secondary structure and the complex biological programme in P. falciparum. The RNA secondary structure characterized in this study has potential applications and important implications regarding the identification of RNA structural elements, which are important for parasite infection and elucidating host-parasite interactions and parasites in the environment.

Noncoding RNAs in Glioblastoma: Emerging Biological Concepts and Potential Therapeutic Implications

Simple Summary

Since the completion of the Human Genome Project, noncoding RNAs (ncRNAs) have emerged as an important class of genetic regulators. Several classes of ncRNAs, which include microRNAs (miRNAs), long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and piwi-interacting RNAs (piRNAs), have been shown to play important roles in controlling developmental and disease processes. In this article, we discuss the potential roles of ncRNAs in regulating glioblastoma (GBM) formation and progression as well as potential strategies to exploit the diagnostic and therapeutic potential of ncRNAs in GBM.

Abstract

Noncoding RNAs (ncRNAs) have emerged as a novel class of genomic regulators, ushering in a new era in molecular biology. With the advent of advanced genetic sequencing technology, several different classes of ncRNAs have been uncovered, including microRNAs (miRNAs), long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and piwi-interacting RNAs (piRNAs), which have been linked to many important developmental and disease processes and are being pursued as clinical and therapeutic targets. Molecular phenotyping studies of glioblastoma (GBM), the most common and lethal cancer of the adult brain, revealed that several ncRNAs are frequently dysregulated in its pathogenesis. Additionally, ncRNAs regulate many important aspects of glioma biology including tumour cell proliferation, migration, invasion, apoptosis, angiogenesis, and self-renewal. Here, we present an overview of the biogenesis of the different classes of ncRNAs, discuss their biological roles, as well as their relevance to gliomagenesis. We conclude by discussing potential approaches to therapeutically target the ncRNAs in clinic.

Messenger RNA 5' NAD+ Capping Is a Dynamic Regulatory Epitranscriptome Mark That Is Required for Proper Response to Abscisic Acid in Arabidopsis

Although eukaryotic messenger RNAs (mRNAs) normally possess a 5' end N⁷-methyl guanosine (m⁷G) cap, a non-canonical 5' nicotinamide adenine dinucleotide (NAD⁺) cap can tag certain transcripts for degradation mediated by the NAD⁺ decapping enzyme DXO1. Despite this importance, whether NAD⁺ capping dynamically responds to specific stimuli to regulate eukaryotic transcriptomes remains unknown. Here, we reveal a link between NAD⁺ capping and tissue- and hormone response-specific mRNA stability. In the absence of DXO1 function, transcripts displaying a high proportion of NAD⁺ capping are instead processed into RNA-dependent RNA polymerase 6-dependent small RNAs, resulting in their continued turnover likely to free the NAD⁺ molecules. Additionally, the NAD⁺-capped transcriptome is significantly remodeled in response to the essential plant hormone abscisic acid in a mechanism that is primarily independent of DXO1. Overall, our findings reveal a previously uncharacterized and essential role of NAD⁺ capping in dynamically regulating transcript stability during specific physiological responses.

New RNA Structural Elements Identified in the Coding Region of the Coxsackie B3 Virus Genome

Here we present a set of new structural elements formed within the open reading frame of the virus, which are highly probable, evolutionarily conserved and may interact with host proteins. This work focused on the coding regions of the CVB3 genome (particularly the V4-, V1-, 2C-, and 3D-coding regions), which, with the exception of the cis-acting replication element (CRE), have not yet been subjected to experimental analysis of their structures. The SHAPE technique, chemical modification with DMS and RNA cleavage with Pb²⁺, were performed in order to characterize the RNA structure. The experimental results were used to improve the computer prediction of the structural models, whereas a phylogenetic analysis was performed to check universality of the newly identified structural elements for twenty CVB3 genomes and 11 other enteroviruses. Some of the RNA motifs turned out to be conserved among different enteroviruses. We also observed that the 3'-terminal region of the genome tends to dimerize in a magnesium concentration-dependent manner. RNA affinity chromatography was used to confirm RNA-protein interactions hypothesized by database searches, leading to the discovery of several interactions, which may be important for virus propagation.

Determination of isoform-specific RNA structure with nanopore long reads

Current methods for determining RNA structure with short-read sequencing cannot capture most differences between distinct transcript isoforms. Here we present RNA structure analysis using nanopore sequencing (PORE-cupine), which combines structure probing using chemical modifications with direct long-read RNA sequencing and machine learning to detect secondary structures in cellular RNAs. PORE-cupine also captures global structural features, such as RNA-binding-protein binding sites and reactivity differences at single-nucleotide variants. We show that shared sequences in different transcript isoforms of the same gene can fold into different structures, highlighting the importance of long-read sequencing for obtaining phase information. We also demonstrate that structural differences between transcript isoforms of the same gene lead to differences in translation efficiency. By revealing isoform-specific RNA structure, PORE-cupine will deepen understanding of the role of structures in controlling gene regulation.

RNA structures in alternative splicing and back-splicing

Alternative splicing greatly expands the transcriptomic and proteomic diversities related to physiological and developmental processes in higher eukaryotes. Splicing of long noncoding RNAs, and back- and trans- splicing further expanded the regulatory repertoire of alternative splicing. RNA structures were shown to play an important role in regulating alternative splicing and back-splicing. Application of novel sequencing technologies made it possible to identify genome-wide RNA structures and interaction networks, which might provide new insights into RNA splicing regulation in vitro to in vivo. The emerging transcription-folding-splicing paradigm is changing our understanding of RNA alternative splicing regulation. Here, we review the insights into the roles and mechanisms of RNA structures in alternative splicing and back-splicing, as well as how disruption of these structures affects alternative splicing and then leads to human diseases. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

mRNA adenosine methylase (MTA) deposits m6A on pri-miRNAs to modulate miRNA biogenesis in Arabidopsis thaliana

Recently N⁶-methyladenosine (m⁶A) methylation has emerged as a biological process with significant impact on cellular functions. However, almost all the research regarding m⁶A methylation has been based on mRNAs. In our research, we focus on how m⁶A methylation affects microRNA (miRNA) biogenesis in Arabidopsis. In brief, we show that m⁶A methylation is necessary to maintain proper levels of mature miRNAs as well as their precursors. m⁶A mark affects pri-miRNA secondary structures and affects the recruitment of the Microprocessor to pri-miRNAs. We also demonstrate the interactions of MTA (m⁶A writer) with other proteins involved in miRNA biogenesis, namely RNA Polymerase II and TOUGH. Our study provides evidence of the role played by m⁶A in plant miRNA biogenesis.

In Arabidopsis thaliana, the METTL3 homolog, mRNA adenosine methylase (MTA) introduces N⁶-methyladenosine (m⁶A) into various coding and noncoding RNAs of the plant transcriptome. Here, we show that an MTA-deficient mutant (mta) has decreased levels of microRNAs (miRNAs) but accumulates primary miRNA transcripts (pri-miRNAs). Moreover, pri-miRNAs are methylated by MTA, and RNA structure probing analysis reveals a decrease in secondary structure within stem–loop regions of these transcripts in mta mutant plants. We demonstrate interaction between MTA and both RNA Polymerase II and TOUGH (TGH), a plant protein needed for early steps of miRNA biogenesis. Both MTA and TGH are necessary for efficient colocalization of the Microprocessor components Dicer-like 1 (DCL1) and Hyponastic Leaves 1 (HYL1) with RNA Polymerase II. We propose that secondary structure of miRNA precursors induced by their MTA-dependent m⁶A methylation status, together with direct interactions between MTA and TGH, influence the recruitment of Microprocessor to plant pri-miRNAs. Therefore, the lack of MTA in mta mutant plants disturbs pri-miRNA processing and leads to the decrease in miRNA accumulation. Furthermore, our findings reveal that reduced miR393b levels likely contributes to the impaired auxin response phenotypes of mta mutant plants.

Determinants of translation efficiency in the evolutionarily-divergent protist Trichomonas vaginalis

BACKGROUND

Trichomonas vaginalis, the causative agent of a prevalent urogenital infection in humans, is an evolutionarily divergent protozoan. Protein-coding genes in T. vaginalis are largely controlled by two core promoter elements, producing mRNAs with short 5' UTRs. The specific mechanisms adopted by T. vaginalis to fine-tune the translation efficiency (TE) of mRNAs remain largely unknown.

RESULTS

Using both computational and experimental approaches, this study investigated two key factors influencing TE in T. vaginalis: codon usage and mRNA secondary structure. Statistical dependence between TE and codon adaptation index (CAI) highlighted the impact of codon usage on mRNA translation in T. vaginalis. A genome-wide interrogation revealed that low structural complexity at the 5' end of mRNA followed closely by a highly structured downstream region correlates with TE variation in this organism. To validate these findings, a synthetic library of 15 synonymous iLOV genes was created, representing five mRNA folding profiles and three codon usage profiles. Fluorescence signals produced by the expression of these synonymous iLOV genes in T. vaginalis were consistent with and validated our in silico predictions.

CONCLUSIONS

This study demonstrates the role of codon usage bias and mRNA secondary structure in TE of T. vaginalis mRNAs, contributing to a better understanding of the factors that influence, and possibly regulate, gene expression in this human pathogen.

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.

N6-methyladenosine and RNA secondary structure affect transcript stability and protein abundance during systemic salt stress in Arabidopsis

After transcription, a messenger RNA (mRNA) is further post-transcriptionally regulated by several features including RNA secondary structure and covalent RNA modifications (specifically N6-methyladenosine, m6A). Both RNA secondary structure and m6A have been demonstrated to regulate mRNA stability and translation and have been independently linked to plant responses to soil salinity levels. However, the effect of m6A on regulating RNA secondary structure and the combinatorial interplay between these two RNA features during salt stress response has yet to be studied. Here, we globally identify RNA-protein interactions and RNA secondary structure during systemic salt stress. This analysis reveals that RNA secondary structure changes significantly during salt stress, and that it is independent of global changes in RNA-protein interactions. Conversely, we find that m6A is anti-correlated with RNA secondary structure in a condition-dependent manner, with salt-specific m6A correlated with a decrease in mRNA secondary structure during salt stress. Taken together, we suggest that salt-specific m6A deposition and the associated loss of RNA secondary structure results in increases in mRNA stability for transcripts encoding abiotic stress response proteins and ultimately increases in protein levels from these stabilized transcripts. In total, our comprehensive analyses reveal important post-transcriptional regulatory mechanisms involved in plant long-term salt stress response and adaptation.

The landscape of eukaryotic mRNPs

The proper regulation of mRNA processing, localization, translation, and degradation occurs on mRNPs. However, the global principles of mRNP organization are poorly understood. We utilize the limited, but existing, information available to present a speculative synthesis of mRNP organization with the following key points. First, mRNPs form a compacted structure due to the inherent folding of RNA. Second, the ribosome is the principal mechanism by which mRNA regions are partially decompacted. Third, mRNPs are 50%-80% protein by weight, consistent with proteins modulating mRNP organization, but also suggesting the majority of mRNA sequences are not directly interacting with RNA-binding proteins. Finally, the ratio of mRNA-binding proteins to mRNAs is higher in the nucleus to allow effective RNA processing and limit the potential for nuclear RNA based aggregation. This synthesis of mRNP understanding provides a model for mRNP biogenesis, structure, and regulation with multiple implications.

Transcriptome-wide-scale-predicted dsRNAs potentially involved in RNA homoeostasis are remarkably excluded from genes with no/very low expression in all developmental stages

RNA interference (RNAi) refers to a conserved posttranscriptional mechanism for the degradation of RNA by short dsRNAs. A genome-wide analysis of mRNAs that are complementary to RNAs of variable length that are transcribed from the full transcriptome and susceptible to being loaded onto Argonaute type 2 was performed through computational searches in the Drosophila model. We report the segments of RNAs that are complementary to mRNAs originating from introns, the exons of mRNAs and lncRNAs as a potential source of siRNAs. A full catalogue of the mRNAs that fulfill these criteria is presented, along with the quantification of multiple annealing. The catalogue was assessed for biological validation using three published lists: two for Ago2-associated RNAs and one for dsRNAs isolated from a crude extract. A broad spectrum of mRNAs were found to theoretically form intermolecular segmental dsRNAs, which should qualify them as Dicer/Ago2 substrates if they exist in vivo. These results suggest a genome-wide scale of mRNA homoeostasis via RNAi metabolism and could extend the known roles of canonical miRNAs and hairpin RNAs. The distribution of the genes for which transcripts are engaged in intermolecular segmental pairing is largely lacking in the gene collections defined as showing no expression in each individual developmental stage from early embryos to adulthood. This trend was also observed for the genes showing very low expression from the 8-12-hour embryonic to larval stage 2. This situation was also suggested by the 3 lists generated with minimal 20-, 25- and 30-base pairing lengths.

RNAdt: An online tutorial and data portal for the RNA structurome era

RNA is not only a passive transporter of genetic information, but also a pivotal player in all domains of life. RNA can regulate gene expression because of its involvement in transcription, mRNA modification and processing, and translation. RNA also possesses other intricate functions such as catalysis, ligand sensing, interaction with biomolecules, response to environment stresses, and information storage. The primary structure of RNA is single stranded, but it always folds into complex secondary and tertiary structures owing to base pairing and effects from the cellular environment. The importance of structure has been increasingly recognized in understanding the myriad functions of RNA. After decades of development, there is a wide range of RNA structure probing techniques. The marriage between structure probing and high-throughput sequencing (HTS) especially enables the measurement of RNA structure on a transcriptomic scale, advancing the advent of the RNA structurome era. Dozens of HTS-associated RNA structure probing methods have been published, so it is urgent to provide a user-friendly and easy-to-use resource for users who are perplexed by selecting the most suitable method for their experiments. Motivated by this demand, we collected currently available HTS-associated RNA structure probing methods and then developed RNAdt (freely accessible at http://www.zhounan.org/rnadt). RNAdt can be used as a web-based tutorial to learn fundamental knowledge of HTS-associated RNA structure probing methods. RNAdt can also be used as a data portal to access HTS data sets from previous RNA structurome studies. At the end of this work, we also provided perspectives on future development of RNA structure probing methods. Our study is expected to facilitate RNA structure probing and ultimately elucidate the connection between RNA structure and biological functions.

Co-temporal Force and Fluorescence Measurements Reveal a Ribosomal Gear Shift Mechanism of Translation Regulation by Structured mRNAs

The movement of ribosomes on mRNA is often interrupted by secondary structures that present mechanical barriers and play a central role in translation regulation. We investigate how ribosomes couple their internal conformational changes with the activity of translocation factor EF-G to unwind mRNA secondary structures using high-resolution optical tweezers with single-molecule fluorescence capability. We find that hairpin opening occurs during EF-G-catalyzed translocation and is driven by the forward rotation of the small subunit head. Modulating the magnitude of the hairpin barrier by force shows that ribosomes respond to strong barriers by shifting their operation to an alternative 7-fold-slower kinetic pathway prior to translocation. Shifting into a slow gear results from an allosteric switch in the ribosome that may allow it to exploit thermal fluctuations to overcome mechanical barriers. Finally, we observe that ribosomes occasionally open the hairpin in two successive sub-codon steps, revealing a previously unobserved translocation intermediate.

The regulation properties of RNA secondary structure in alternative splicing

The RNA secondary structure is important for many functional processes in the cell. The secondary and tertiary structures of cellular RNAs are essential for the activity of these molecules in processes such as transcription, splicing, translation, and localization. New high-throughput analytical methods, including next generation sequencing, have allowed for the in-depth characterization of the 'RNA structurome': a new term describing how the RNA structure controls the activity of RNA by itself and how it regulates the expression of genes. In this review, we present many examples of the influence of structural motifs of RNA, long range interactions and global RNA structure on the alternative splicing processes. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

High-throughput determination of RNA structures

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

Effects of Arbovirus Multi-Host Life Cycles on Dinucleotide and Codon Usage Patterns

Arthropod-borne viruses (arboviruses) of vertebrates including dengue, zika, chikungunya, Rift Valley fever, and blue tongue viruses cause extensive morbidity and mortality in humans, agricultural animals, and wildlife across the globe. As obligate intercellular pathogens, arboviruses must be well adapted to the cellular and molecular environment of both their arthropod (invertebrate) and vertebrate hosts, which are vastly different due to hundreds of millions of years of separate evolution. Here we discuss the comparative pressures on arbovirus RNA genomes as a result of a dual host life cycle, focusing on pressures that do not alter amino acids. We summarize what is currently known about arboviral genetic composition, such as dinucleotide and codon usage, and how cyclical infection of vertebrate and invertebrate hosts results in different genetic profiles compared with single-host viruses. To serve as a comparison, we compile what is known about arthropod tRNA, dinucleotide, and codon usages and compare this with vertebrates. Additionally, we discuss the potential roles of genetic robustness in arboviral evolution and how it may vary from other viruses. Overall, both arthropod and vertebrate hosts influence the resulting genetic composition of arboviruses, but a great deal remains to be investigated.

Mapping the dsRNA World

Long double-stranded RNAs (dsRNAs) are abundantly expressed in animals, in which they frequently occur in introns and 3' untranslated regions of mRNAs. Functions of long, cellular dsRNAs are poorly understood, although deficiencies in adenosine deaminases that act on RNA, or ADARs, promote their recognition as viral dsRNA and an aberrant immune response. Diverse dsRNA-binding proteins bind cellular dsRNAs, hinting at additional roles. Understanding these roles is facilitated by mapping the genomic locations that express dsRNA in various tissues and organisms. ADAR editing provides a signature of dsRNA structure in cellular transcripts. In this review, we detail approaches to map ADAR editing sites and dsRNAs genome-wide, with particular focus on high-throughput sequencing methods and considerations for their successful application to the detection of editing sites and dsRNAs.

Inverted repeat structures are associated with essential and highly expressed genes on C. elegans autosome distal arms

Complementary sequences in cellular transcripts base-pair to form double-stranded RNA (dsRNA) structures. Because transposon-derived repeats often give rise to self-complementary sequences, dsRNA structures are prevalent in eukaryotic genomes, typically occurring in gene introns and untranslated regions (UTRs). However, the regulatory impact of double-stranded structures within genes is not fully understood. We used three independent methods to define loci in Caenorhabditis elegans predicted to form dsRNA and correlated these structures with patterns of gene expression, gene essentiality, and genome organization. As previously observed, dsRNA loci are enriched on distal arms of C. elegans autosomes, where genes typically show less conservation and lower overall expression. In contrast, we find that dsRNAs are associated with essential genes on autosome arms, and dsRNA-associated genes exhibit higher-than-expected expression and histone modification patterns associated with transcriptional elongation. Genes with significant repetitive sequence content are also highly expressed, and, thus, observed gene expression trends may relate either to dsRNA structures or to repeat content. Our results raise the possibility that as-yet-undescribed mechanisms promote expression of loci that produce dsRNAs, despite their well-characterized roles in gene silencing.

Revealing lncRNA Structures and Interactions by Sequencing-Based Approaches

Long noncoding RNAs (lncRNAs) have emerged as significant players in almost every level of gene function and regulation. Thus, characterizing the structures and interactions of lncRNAs is essential for understanding their mechanistic roles in cells. Through a combination of (bio)chemical approaches and automated capillary and high-throughput sequencing (HTS), the complexity and diversity of RNA structures and interactions has been revealed in the transcriptomes of multiple species. These methods have uncovered important biological insights into the mechanistic and functional roles of lncRNA in gene expression and RNA metabolism, as well as in development and disease. In this review, we summarize the latest sequencing strategies to reveal RNA structure, RNA-RNA, RNA-DNA, and RNA-protein interactions, and highlight the recent applications of these approaches to map functional lncRNAs. We discuss the advantages and limitations of these strategies, and provide recommendations to further advance methodologies capable of mapping RNA structure and interactions in order to discover new biology of lncRNAs and decipher their molecular mechanisms and implication in diseases.

RNA virus evasion of nonsense-mediated decay

Nonsense-mediated decay (NMD) is a host RNA control pathway that removes aberrant transcripts with long 3' untranslated regions (UTRs) due to premature termination codons (PTCs) that arise through mutation or defective splicing. To maximize coding potential, RNA viruses often contain internally located stop codons that should also be prime targets for NMD. Using an agroinfiltration-based NMD assay in Nicotiana benthamiana, we identified two segments conferring NMD-resistance in the carmovirus Turnip crinkle virus (TCV) genome. The ribosome readthrough structure just downstream of the TCV p28 termination codon stabilized an NMD-sensitive reporter as did a frameshifting element from umbravirus Pea enation mosaic virus. In addition, a 51-nt unstructured region (USR) at the beginning of the TCV 3' UTR increased NMD-resistance 3-fold when inserted into an unrelated NMD-sensitive 3' UTR. Several additional carmovirus 3' UTRs also conferred varying levels of NMD resistance depending on the construct despite no sequence similarity in the analogous region. Instead, these regions displayed a marked lack of RNA structure immediately following the NMD-targeted stop codon. NMD-resistance was only slightly reduced by conversion of 19 pyrimidines in the USR to purines, but resistance was abolished when a 2-nt mutation was introduced downstream of the USR that substantially increased the secondary structure in the USR through formation of a stable hairpin. The same 2-nt mutation also enhanced the NMD susceptibility of a subgenomic RNA expressed independently of the genomic RNA. The conserved lack of RNA structure among most carmoviruses at the 5' end of their 3' UTR could serve to enhance subgenomic RNA stability, which would increase expression of the encoded capsid protein that also functions as the RNA silencing suppressor. These results demonstrate that the TCV genome has features that are inherently NMD-resistant and these strategies could be widespread among RNA viruses and NMD-resistant host mRNAs with long 3' UTRs.

Plant 24-nt reproductive phasiRNAs from intramolecular duplex mRNAs in diverse monocots

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

RNA, Action through Interactions

As transcription of the human genome is quite pervasive, it is possible that many novel functions of the noncoding genome have yet to be identified. Often the noncoding genome's functions are carried out by their RNA transcripts, which may rely on their structures and/or extensive interactions with other molecules. Recent technology developments are transforming the fields of RNA biology from studying one RNA at a time to transcriptome-wide mapping of structures and interactions. Here, we highlight the recent advances in transcriptome-wide RNA interaction analysis. These technologies revealed surprising versatility of RNA to participate in diverse molecular systems. For example, tens of thousands of RNA-RNA interactions have been revealed in cultured cells as well as in mouse brain, including interactions between transposon-produced transcripts and mRNAs. In addition, most transcription start sites in the human genome are associated with noncoding RNA transcribed from other genomic loci. These recent discoveries expanded our understanding of RNAs' roles in chromatin organization, gene regulation, and intracellular signaling.

Analyses of mRNA structure dynamics identify embryonic gene regulatory programs

RNA folding plays a crucial role in RNA function. However, knowledge of the global structure of the transcriptome is limited to cellular systems at steady state, thus hindering the understanding of RNA structure dynamics during biological transitions and how it influences gene function. Here, we characterized mRNA structure dynamics during zebrafish development. We observed that on a global level, translation guides structure rather than structure guiding translation. We detected a decrease in structure in translated regions and identified the ribosome as a major remodeler of RNA structure in vivo. In contrast, we found that 3' untranslated regions (UTRs) form highly folded structures in vivo, which can affect gene expression by modulating microRNA activity. Furthermore, dynamic 3'-UTR structures contain RNA-decay elements, such as the regulatory elements in nanog and ccna1, two genes encoding key maternal factors orchestrating the maternal-to-zygotic transition. These results reveal a central role of RNA structure dynamics in gene regulatory programs.

Adenine Enrichment at the Fourth CDS Residue in Bacterial Genes Is Consistent with Error Proofing for +1 Frameshifts

Beyond selection for optimal protein functioning, coding sequences (CDSs) are under selection at the RNA and DNA levels. Here, we identify a possible signature of “dual-coding,” namely extensive adenine (A) enrichment at bacterial CDS fourth sites. In 99.07% of studied bacterial genomes, fourth site A use is greater than expected given genomic A-starting codon use. Arguing for nucleotide level selection, A-starting serine and arginine second codons are heavily utilized when compared with their non-A starting synonyms. Several models have the ability to explain some of this trend. In part, A-enrichment likely reduces 5′ mRNA stability, promoting translation initiation. However T/U, which may also reduce stability, is avoided. Further, +1 frameshifts on the initiating ATG encode a stop codon (TGA) provided A is the fourth residue, acting either as a frameshift “catch and destroy” or a frameshift stop and adjust mechanism and hence implicated in translation initiation. Consistent with both, genomes lacking TGA stop codons exhibit weaker fourth site A-enrichment. Sequences lacking a Shine–Dalgarno sequence and those without upstream leader genes, that may be more error prone during initiation, have greater utilization of A, again suggesting a role in initiation. The frameshift correction model is consistent with the notion that many genomic features are error-mitigation factors and provides the first evidence for site-specific out of frame stop codon selection. We conjecture that the NTG universal start codon may have evolved as a consequence of TGA being a stop codon and the ability of NTGA to rapidly terminate or adjust a ribosome.

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.

RNA–protein interactions often depend on the recognition of extended RNA elements but the identification of these motifs is challenging. Here, the authors present a global integrated approach to analyze RNA–protein binding landscapes, mapping extended RNA interaction motifs for four RNA-binding proteins.

New Era of Studying RNA Secondary Structure and Its Influence on Gene Regulation in Plants

The dynamic structure of RNA plays a central role in post-transcriptional regulation of gene expression such as RNA maturation, degradation, and translation. With the rise of next-generation sequencing, the study of RNA structure has been transformed from in vitro low-throughput RNA structure probing methods to in vivo high-throughput RNA structure profiling. The development of these methods enables incremental studies on the function of RNA structure to be performed, revealing new insights of novel regulatory mechanisms of RNA structure in plants. Genome-wide scale RNA structure profiling allows us to investigate general RNA structural features over 10s of 1000s of mRNAs and to compare RNA structuromes between plant species. Here, we provide a comprehensive and up-to-date overview of: (i) RNA structure probing methods; (ii) the biological functions of RNA structure; (iii) genome-wide RNA structural features corresponding to their regulatory mechanisms; and (iv) RNA structurome evolution in plants.

Rice In Vivo RNA Structurome Reveals RNA Secondary Structure Conservation and Divergence in Plants

RNA secondary structure plays a critical role in gene regulation. Rice (Oryza sativa) is one of the most important food crops in the world. However, RNA structure in rice has scarcely been studied. Here, we have successfully generated in vivo Structure-seq libraries in rice. We found that the structural flexibility of mRNAs might associate with the dynamics of biological function. Higher N⁶-methyladenosine (m⁶A) modification tends to have less RNA structure in 3' UTR, whereas GC content does not significantly affect in vivo mRNA structure to maintain efficient biological processes such as translation. Comparative analysis of RNA structurome between rice and Arabidopsis revealed that higher GC content does not lead to stronger structure and less RNA structural flexibility. Moreover, we found a weak correlation between sequence and structure conservation of the orthologs between rice and Arabidopsis. The conservation and divergence of both sequence and in vivo RNA structure corresponds to diverse and specific biological processes. Our results indicate that RNA secondary structure might offer a separate layer of selection to the sequence between monocot and dicot. Therefore, our study implies that RNA structure evolves differently in various biological processes to maintain robustness in development and adaptational flexibility during angiosperm evolution.

RNA editing by ADAR1 leads to context-dependent transcriptome-wide changes in RNA secondary structure

Adenosine deaminase acting on RNA 1 (ADAR1) is the master RNA editor, catalyzing the deamination of adenosine to inosine. RNA editing is vital for preventing abnormal activation of cytosolic nucleic acid sensing pathways by self-double-stranded RNAs. Here we determine, by parallel analysis of RNA secondary structure sequencing (PARS-seq), the global RNA secondary structure changes in ADAR1 deficient cells. Surprisingly, ADAR1 silencing resulted in a lower global double-stranded to single-stranded RNA ratio, suggesting that A-to-I editing can stabilize a large subset of imperfect RNA duplexes. The duplexes destabilized by editing are composed of vastly complementary inverted Alus found in untranslated regions of genes performing vital biological processes, including housekeeping functions and type-I interferon responses. They are predominantly cytoplasmic and generally demonstrate higher ribosomal occupancy. Our findings imply that the editing effect on RNA secondary structure is context dependent and underline the intricate regulatory role of ADAR1 on global RNA secondary structure.

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.

A Global View of RNA-Protein Interactions Identifies Post-transcriptional Regulators of Root Hair Cell Fate

The Arabidopsis thaliana root epidermis is comprised of two cell types, hair and nonhair cells, which differentiate from the same precursor. Although the transcriptional programs regulating these events are well studied, post-transcriptional factors functioning in this cell fate decision are mostly unknown. Here, we globally identify RNA-protein interactions and RNA secondary structure in hair and nonhair cell nuclei. This analysis reveals distinct structural and protein binding patterns across both transcriptomes, allowing identification of differential RNA binding protein (RBP) recognition sites. Using these sequences, we identify two RBPs that regulate hair cell development. Specifically, we find that SERRATE functions in a microRNA-dependent manner to inhibit hair cell fate, while also terminating growth of root hairs mostly independent of microRNA biogenesis. In addition, we show that GLYCINE-RICH PROTEIN 8 promotes hair cell fate while alleviating phosphate starvation stress. In total, this global analysis reveals post-transcriptional regulators of plant root epidermal cell fate.

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.

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.

RNA structure, binding, and coordination in Arabidopsis

From the moment of transcription, up through degradation, each RNA transcript is bound by an ever-changing cohort of RNA binding proteins. The binding of these proteins is regulated by both the primary RNA sequence, as well as the intramolecular RNA folding, or secondary structure, of the transcript. Thus, RNA secondary structure regulates many post-transcriptional processes. With the advent of next generation sequencing, several techniques have been developed to generate global landscapes of both RNA-protein interactions and RNA secondary structure. In this review, we describe the current state of the field detailing techniques to globally interrogate RNA secondary structure and/or RNA-protein interaction sites, as well as our current understanding of these features in the transcriptome of the model plant Arabidopsis thaliana. WIREs RNA 2017, 8:e1426. doi: 10.1002/wrna.1426 For further resources related to this article, please visit the WIREs website.

Translation regulation in plants: an interesting past, an exciting present and a promising future

Changes in gene expression are at the core of most biological processes, from cell differentiation to organ development, including the adaptation of the whole organism to the ever-changing environment. Although the central role of transcriptional regulation is solidly established and the general mechanisms involved in this type of regulation are relatively well understood, it is clear that regulation at a translational level also plays an essential role in modulating gene expression. Despite the large number of examples illustrating the critical role played by translational regulation in determining the expression levels of a gene, our understanding of the molecular mechanisms behind such types of regulation has been slow to emerge. With the recent development of high-throughput approaches to map and quantify different critical parameters affecting translation, such as RNA structure, protein-RNA interactions and ribosome occupancy at the genome level, a renewed enthusiasm toward studying translation regulation is warranted. The use of these new powerful technologies in well-established and uncharacterized translation-dependent processes holds the promise to decipher the likely complex and diverse, but also fascinating, mechanisms behind the regulation of translation.

Measuring RNA structure transcriptome-wide with icSHAPE

RNA molecules can be found at the heart of nearly every aspect of gene regulation: from gene expression to protein translation. The ability of RNA molecules to fold into intricate structures guides their function. Chemical methods to measure RNA structure have been part of the RNA biologists toolkit for several decades. These methods, although often cumbersome and difficult to perform on large RNAs, are notable for their accuracy and precision of structural measurements. Recent extension of these methods to transcriptome-wide analyses has opened the door to interrogating the structure of complete RNA molecules inside cells. Within this manuscript we describe the biochemical basis for the methodology behind a novel technology, icSHAPE, which measures RNA flexibility and single-strandedness in RNA. Novel methods such as icSHAPE have greatly expanded our understanding of RNA function and have paved the way to expansive analyses of large groups of RNA structures as they function inside the native environment of the cell.