Regulatory impact of RNA secondary structure across the Arabidopsis transcriptome

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.

miRNA-like secondary structures in maize (Zea mays) genes and transposable elements correlate with small RNAs, methylation, and expression

RNA molecules carry information in their primary sequence and also their secondary structure. Secondary structure can confer important functional information, but it is also a signal for an RNAi-like host epigenetic response mediated by small RNAs (smRNAs). In this study, we used two bioinformatic methods to predict local secondary structures across features of the maize genome, focusing on small regions that had similar folding properties to pre-miRNA loci. We found miRNA-like secondary structures to be common in genes and most, but not all, superfamilies of RNA and DNA transposable elements (TEs). The miRNA-like regions map to a higher diversity of smRNAs than regions without miRNA-like structure, explaining up to 27% of variation in smRNA mapping for some TE superfamilies. This mapping bias is more pronounced among putatively autonomous TEs relative to nonautonomous TEs. Genome-wide, miRNA-like regions are also associated with elevated methylation levels, particularly in the CHH context. Among genes, those with miRNA-like secondary structure are 1.5-fold more highly expressed, on average, than other genes. However, these genes are also more variably expressed across the 26 nested association mapping founder lines, and this variability positively correlates with the number of mapping smRNAs. We conclude that local miRNA-like structures are a nearly ubiquitous feature of expressed regions of the maize genome, that they correlate with higher smRNA mapping and methylation, and that they may represent a trade-off between functional requirements and the potentially negative consequences of smRNA production.

Translation initiation at AUG and non-AUG triplets in plants

In plants and other eukaryotes, precise selection of translation initiation site (TIS) on mRNAs shapes the proteome in response to cellular events or environmental cues. The canonical translation of mRNAs initiates at a 5' proximal AUG codon in a favorable context. However, the coding and non-coding regions of plant genomes contain numerous unannotated alternative AUG and non-AUG TISs. Determining how and why these unexpected and prevalent TISs are activated in plants has emerged as an exciting research area. In this review, we focus on the selection of plant TISs and highlight studies that revealed previously unannotated TISs used in vivo via comparative genomics and genome-wide profiling of ribosome positioning and protein N-terminal ends. The biological signatures of non-AUG TIS-initiated open reading frames (ORFs) in plants are also discussed. We describe what is understood about cis-regulatory RNA elements and trans-acting eukaryotic initiation factors (eIFs) in the site selection for translation initiation by featuring the findings in plants along with supporting findings in non-plant species. The prevalent, unannotated TISs provide a hidden reservoir of ORFs that likely help reshape plant proteomes in response to developmental or environmental cues. These findings underscore the importance of understanding the mechanistic basis of TIS selection to functionally annotate plant genomes, especially for crops with large genomes.

SeedTransNet: a directional translational network revealing regulatory patterns during seed maturation and germination

Seed maturation is the developmental process that prepares the embryo for the desiccated waiting period before germination. It is associated with a series of physiological changes leading to the establishment of seed dormancy, seed longevity, and desiccation tolerance. We studied translational changes during seed maturation and observed a gradual reduction in global translation during seed maturation. Transcriptome and translatome profiling revealed specific reduction in the translation of thousands of genes. By including previously published data on germination and seedling establishment, a regulatory network based on polysome occupancy data was constructed: SeedTransNet. Network analysis predicted translational regulatory pathways involving hundreds of genes with distinct functions. The network identified specific transcript sequence features suggesting separate translational regulatory circuits. The network revealed several seed maturation-associated genes as central nodes, and this was confirmed by specific seed phenotypes of the respective mutants. One of the regulators identified, an AWPM19 family protein, PM19-Like1 (PM19L1), was shown to regulate seed dormancy and longevity. This putative RNA-binding protein also affects the translational regulation of its target mRNA, as identified by SeedTransNet. Our data show the usefulness of SeedTransNet in identifying regulatory pathways during seed phase transitions.

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.

Small and intermediate size structural RNAs in the unicellular parasite Cryptosporidium parvum as revealed by sRNA-seq and comparative genomics

Small and intermediate-size noncoding RNAs (sRNAs and is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Cryptosporidium parvum, an obligate zoonotic protist parasite responsible for the diarrhoeal disease cryptosporidiosis. Using Illumina sequencing of a small RNA library, a systematic identification of novel small and is-ncRNAs was performed in C. parvum excysted sporozoites. A total of 79 novel is-ncRNA candidates, including antisense, intergenic and intronic is-ncRNAs, were identified, including 7 new small nucleolar RNAs (snoRNAs). Expression of select novel is-ncRNAs was confirmed by RT-PCR. Phylogenetic conservation was analysed using covariance models (CMs) in related Cryptosporidium and apicomplexan parasite genome sequences. A potential new type of small ncRNA derived from tRNA fragments was observed. Overall, a deep profiling analysis of novel is-ncRNAs in C. parvum and related species revealed structural features and conservation of these novel is-ncRNAs. Covariance models can be used to detect is-ncRNA genes in other closely related parasites. These findings provide important new sequences for additional functional characterization of novel is-ncRNAs in the protist pathogen C. parvum.

Two evolutionarily duplicated domains individually and post-transcriptionally control SWEET expression for phloem transport

Cell type-specific gene expression is critical for the specialized functions within multicellular organisms. In Arabidopsis, SWEET11 and SWEET12 sugar transporters are specifically expressed in phloem parenchyma (PP) cells and are responsible for sucrose efflux from the PP, the first step of a two-step apoplasmic phloem-loading strategy that initiates the long-distance transport of sugar from leaves to nonphotosynthetic sink tissues. However, we know nothing about what determines the PP cell-specific expression of these SWEETs. Sequence deletions, histochemical β-glucuronidase (GUS) analysis, cross-sectioning, live-cell imaging, and evolutionary analysis were used to elucidate domains responsible for PP specificity, while a Förster resonance energy transfer (FRET) sensor-based transport assay was used to determine whether substrate specificity coevolved with PP specificity. We identified two domains in the Arabidopsis SWEET11 coding sequence that, along with its promoter (including 5' UTR), regulate PP-specific expression at the post-transcriptional level, probably involving RNA-binding proteins. This mechanism is conserved among vascular plants but independent of transport substrate specificity. We conclude that two evolutionarily duplicated coding sequence domains are essential and individually sufficient for PP-specific expression of SWEET11. We also provide a crucial experimental tool to study PP physiology and development.

Novel insights into the pervasive role of RNA structure in post-transcriptional regulation of gene expression in plants

RNA folding is an intrinsic property of RNA that serves a key role in every step of post-transcriptional regulation of gene expression, from RNA maturation to translation in plants. Recent developments of genome-wide RNA structure profiling methods have transformed research in this area enabling focus to shift from individual molecules to the study of tens of thousands of RNAs. Here, we provide a comprehensive review of recent advances in the field. We discuss these new insights of RNA structure functionality within the context of post-transcriptional regulation including mRNA maturation, translation, and RNA degradation in plants. Notably, we also provide an overview of how plants exhibit different RNA structures in response to environmental changes.

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 initiation of RNA interference (RNAi) in plants

When an mRNA enters into the RNA degradation pathway called RNA interference (RNAi), it is cleaved into small interfering RNAs (siRNAs) that then target complementary mRNAs for destruction. The consequence of entry into RNAi is mRNA degradation, post-transcriptional silencing and in some cases transcriptional silencing. RNAi functions as a defense against transposable element and virus activity, and in plants, RNAi additionally plays a role in development by regulating some genes. However, it is unknown how specific transcripts are selected for RNAi, and how most genic mRNAs steer clear. This Current Opinion article explores the key question of how RNAs are selected for entry into RNAi, and proposes models that enable the cell to distinguish between transcripts to translate versus destroy.

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.

Eating microRNAs: pharmacological opportunities for cross-kingdom regulation and implications in host gene and gut microbiota modulation

Cross-kingdom communication via non-coding RNAs is a recent discovery. Exogenous microRNAs (exog-miRNAs) mainly enter the host via the diet. Generally considered unstable in the gastrointestinal tract, some exogenous RNAs may resist these conditions, especially if transported in extracellular vesicles. They could then reach the intestines and more probably exert a regulatory effect. We give an overview of recent discoveries concerning dietary miRNAs, possible ways of enhancing their resistance to food processing and gut conditions, their transport in extracellular vesicles (animal- and plant-origin) and possible biological effects on recipient cells after ingestion. We critically focus on what we believe are the most relevant data for future pharmacological development of dietary miRNAs as therapeutic agents. Finally, we discuss the miRNA-mediated cross-kingdom regulation between diet, host and the gut microbiota. We conclude that, despite many obstacles and challenges, extracellular miRNAs are serious candidates to be targeted pharmacologically for development of new therapeutic agents.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Genome-scale molecular principles of mRNA half-life regulation in yeast

Precise control of protein and messenger RNA (mRNA) degradation is essential for cellular metabolism and homeostasis. Controlled and specific degradation of both molecular species necessitates their engagements with the respective degradation machineries; this engagement involves a disordered/unstructured segment of the substrate traversing the degradation tunnel of the machinery and accessing the catalytic sites. However, while molecular factors influencing protein degradation have been extensively explored on a genome scale, and in multiple organisms, such a comprehensive understanding remains missing for mRNAs. Here, we analyzed multiple genome-scale experimental yeast mRNA half-life data in light of experimentally derived mRNA secondary structures and protein binding data, along with high-resolution X-ray crystallographic structures of the RNase machines. Results unraveled a consistent genome-scale trend that mRNAs comprising longer terminal and/or internal unstructured segments have significantly shorter half-lives; the lengths of the 5'-terminal, 3'-terminal, and internal unstructured segments that affect mRNA half-life are compatible with molecular structures of the 5' exo-, 3' exo-, and endoribonuclease machineries. Sequestration into ribonucleoprotein complexes elongates mRNA half-life, presumably by burying ribonuclease engagement sites under oligomeric interfaces. After gene duplication, differences in terminal unstructured lengths, proportions of internal unstructured segments, and oligomerization modes result in significantly altered half-lives of paralogous mRNAs. Side-by-side comparison of molecular principles underlying controlled protein and mRNA degradation in yeast unravels their remarkable mechanistic similarities and suggests how the intrinsic structural features of the two molecular species, at two different levels of the central dogma, regulate their half-lives on genome scale.

Translational gene regulation in plants: A green new deal

The molecular machinery for protein synthesis is profoundly similar between plants and other eukaryotes. Mechanisms of translational gene regulation are embedded into the broader network of RNA-level processes including RNA quality control and RNA turnover. However, over eons of their separate history, plants acquired new components, dropped others, and generally evolved an alternate way of making the parts list of protein synthesis work. Research over the past 5 years has unveiled how plants utilize translational control to defend themselves against viruses, regulate translation in response to metabolites, and reversibly adjust translation to a wide variety of environmental parameters. Moreover, during seed and pollen development plants make use of RNA granules and other translational controls to underpin developmental transitions between quiescent and metabolically active stages. The economics of resource allocation over the daily light-dark cycle also include controls over cellular protein synthesis. Important new insights into translational control on cytosolic ribosomes continue to emerge from studies of translational control mechanisms in viruses. Finally, sketches of coherent signaling pathways that connect external stimuli with a translational response are emerging, anchored in part around TOR and GCN2 kinase signaling networks. These again reveal some mechanisms that are familiar and others that are different from other eukaryotes, motivating deeper studies on translational control in plants. This article is categorized under: Translation > Translation Regulation RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Structure-seq2 probing of RNA structure upon amino acid starvation reveals both known and novel RNA switches in Bacillus subtilis

RNA structure influences numerous processes in all organisms. In bacteria, these processes include transcription termination and attenuation, small RNA and protein binding, translation initiation, and mRNA stability, and can be regulated via metabolite availability and other stresses. Here we use Structure-seq2 to probe the in vivo RNA structurome of Bacillus subtilis grown in the presence and absence of amino acids. Our results reveal that amino acid starvation results in lower overall dimethyl sulfate (DMS) reactivity of the transcriptome, indicating enhanced protection owing to protein binding or RNA structure. Starvation-induced changes in DMS reactivity correlated inversely with transcript abundance changes. This correlation was particularly pronounced in genes associated with the stringent response and CodY regulons, which are involved in adaptation to nutritional stress, suggesting that RNA structure contributes to transcript abundance change in regulons involved in amino acid metabolism. Structure-seq2 accurately reported on four known amino acid-responsive riboswitches: T-box, SAM, glycine, and lysine riboswitches. Additionally, we discovered a transcription attenuation mechanism that reduces yfmG expression when amino acids are added to the growth medium. We also found that translation of a leader peptide (YfmH) encoded just upstream of yfmG regulates yfmG expression. Our results are consistent with a model in which a slow rate of yfmH translation caused by limitation of the amino acids encoded in YfmH prevents transcription termination in the yfmG leader region by favoring formation of an overlapping antiterminator structure. This novel RNA switch offers a way to simultaneously monitor the levels of multiple amino acids.

On the Way to Understanding the Interplay between the RNA Structure and Functions in Cells: A Genome-Wide Perspective

RNAs adopt specific structures in order to perform their biological activities. The structure of RNA is an important layer of gene expression regulation, and can impact a plethora of cellular processes, starting with transcription, RNA processing, and translation, and ending with RNA turnover. The development of high-throughput technologies has enabled a deeper insight into the sophisticated interplay between the structure of the cellular transcriptome and the living cells environment. In this review, we present the current view on the RNA structure in vivo resulting from the most recent transcriptome-wide studies in different organisms, including mammalians, yeast, plants, and bacteria. We focus on the relationship between the mRNA structure and translation, mRNA stability and degradation, protein binding, and RNA posttranscriptional modifications.

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.

Intact RNA structurome reveals mRNA structure-mediated regulation of miRNA cleavage in vivo

MicroRNA (miRNA)-mediated cleavage is involved in numerous essential cellular pathways. miRNAs recognize target RNAs via sequence complementarity. In addition to complementarity, in vitro and in silico studies have suggested that RNA structure may influence the accessibility of mRNAs to miRNA-induced silencing complexes (miRISCs), thereby affecting RNA silencing. However, the regulatory mechanism of mRNA structure in miRNA cleavage remains elusive. We investigated the role of in vivo RNA secondary structure in miRNA cleavage by developing the new CAP-STRUCTURE-seq method to capture the intact mRNA structurome in Arabidopsis thaliana. This approach revealed that miRNA target sites were not structurally accessible for miRISC binding prior to cleavage in vivo. Instead, we found that the unfolding of the target site structure plays a key role in miRISC activity in vivo. We found that the single-strandedness of the two nucleotides immediately downstream of the target site, named Target Adjacent nucleotide Motif, can promote miRNA cleavage but not miRNA binding, thus decoupling target site binding from cleavage. Our findings demonstrate that mRNA structure in vivo can modulate miRNA cleavage, providing evidence of mRNA structure-dependent regulation of biological processes.

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.

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.

HiPR: High-throughput probabilistic RNA structure inference

Recent high-throughput structure-sensitive genome-wide sequencing-based assays have enabled large-scale studies of RNA structure, and robust transcriptome-wide computational prediction of individual RNA structures across RNA classes from these assays has potential to further improve the prediction accuracy. Here, we describe HiPR, a novel method for RNA structure prediction at single-nucleotide resolution that combines high-throughput structure probing data (DMS-seq, DMS-MaPseq) with a novel probabilistic folding algorithm. On validation data spanning a variety of RNA classes, HiPR often increases accuracy for predicting RNA structures, giving researchers new tools to study RNA structure.

5'-UTR and ORF elements, as well as the 3'-UTR regulate the translation of Cyclin

Mammalian oocyte maturation is wholly dependent on the translation of accumulated maternal transcripts. Therefore, measuring the translation of specific genes, especially Ccnb1 and Ccnb2, which are key regulators of the oocyte cell cycle in mice, is essential to monitor oocyte cell cycle progression. For this purpose, almost all previous research has used a reporter construct containing the 3'-untranslated region (UTR) of Ccnb. It is based on the concept that the 3'-UTR is the main modulator of translation. Here, we investigated the expression pattern of Renilla luciferase (RL) reporters combining the 5'-UTR and/or open reading frame (ORF) as well as the 3'-UTR (RL-3', 5'-RL-3', RL-ORF-3', and 5'-RL-ORF-3') of Ccnb1 and Ccnb2 in somatic cells and mouse oocytes. The addition of the 5'-UTR and/or ORF of Ccnb altered the expression of the RL-3' reporter in HEK293T cells and mouse oocytes. The ORF tended to suppress RL expression, whereas the 5'-UTR enhanced the expression in most cases. The increased rate in expression was the highest when only the 3'-UTR of Ccnb1 (RL-3') was used, whereas the 5'-RL-ORF-3' reporter showed a relatively lower increase during oocyte maturation. For Ccnb2, the RL-ORF-3' reporter showed the largest increase, and other reporters exhibited a similar increase in expression during oocyte maturation. Results show that the expression of these genes is modulated not only by the 3'-UTR but also by the 5'-UTR and ORF. Therefore, special caution should be taken when using only the 3'-UTR to monitor the expression of specific genes.

Tissue-specific changes in the RNA structurome mediate salinity response in Arabidopsis

Little is known concerning the effects of abiotic factors on in vivo RNA structures. We applied Structure-seq to assess the in vivo mRNA structuromes of Arabidopsis thaliana under salinity stress, which negatively impacts agriculture. Structure-seq utilizes dimethyl sulfate reactivity to identify As and Cs that lack base-pairing or protection. Salt stress refolded transcripts differentially in root versus shoot, evincing tissue specificity of the structurome. Both tissues exhibited an inverse correlation between salt stress-induced changes in transcript reactivity and changes in abundance, with stress-related mRNAs showing particular structural dynamism. This inverse correlation is more pronounced in mRNAs wherein the mean reactivity of the 5'UTR, CDS, and 3'UTR concertedly change under salinity stress, suggesting increased susceptibility to abundance control mechanisms in transcripts exhibiting this phenomenon, which we name "concordancy." Concordant salinity-induced increases in reactivity were notably observed in photosynthesis genes, thereby implicating mRNA structural loss in the well-known depression of photosynthesis by salt stress. Overall, changes in secondary structure appear to impact mRNA abundance, molding the functional specificity of the transcriptome under stress.

Seed-Stored mRNAs that Are Specifically Associated to Monosomes Are Translationally Regulated during Germination

The life cycle of many organisms includes a quiescent stage, such as bacterial or fungal spores, insect larvae, or plant seeds. Common to these stages is their low water content and high survivability during harsh conditions. Upon rehydration, organisms need to reactivate metabolism and protein synthesis. Plant seeds contain many mRNAs that are transcribed during seed development. Translation of these mRNAs occurs during early seed germination, even before the requirement of transcription. Therefore, stored mRNAs are postulated to be important for germination. How these mRNAs are stored and protected during long-term storage is unknown. The aim of this study was to investigate how mRNAs are stored in dry seeds and whether they are indeed translated during seed germination. We investigated seed polysome profiles and the mRNAs and protein complexes that are associated with these ribosomes in seeds of the model organism Arabidopsis (Arabidopsis thaliana). We showed that most stored mRNAs are associated with monosomes in dry seeds; therefore, we focus on monosomes in this study. Seed ribosome complexes are associated with mRNA-binding proteins, stress granule, and P-body proteins, which suggests regulated packing of seed mRNAs. Interestingly, ∼17% of the mRNAs that are specifically associated with monosomes are translationally up-regulated during seed germination. These mRNAs are transcribed during seed maturation, suggesting a role for this developmental stage in determining the translational fate of mRNAs during early germination.

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.

Coding or Noncoding, the Converging Concepts of RNAs

Technological advances over the past decade have unraveled the remarkable complexity of RNA. The identification of small peptides encoded by long non-coding RNAs (lncRNAs) as well as regulatory functions mediated by non-coding regions of mRNAs have further complicated our understanding of the multifaceted functions of RNA. In this review, we summarize current evidence pointing to dual roles of RNA molecules defined by their coding and non-coding potentials. We also discuss how the emerging roles of RNA transform our understanding of gene expression and evolution.

Conserved Secondary Structures in Viral mRNAs

RNA secondary structure in untranslated and protein coding regions has been shown to play an important role in regulatory processes and the viral replication cycle. While structures in non-coding regions have been investigated extensively, a thorough overview of the structural repertoire of protein coding mRNAs, especially for viruses, is lacking. Secondary structure prediction of large molecules, such as long mRNAs remains a challenging task, as the contingent of structures a sequence can theoretically fold into grows exponentially with sequence length. We applied a structure prediction pipeline to Viral Orthologous Groups that first identifies the local boundaries of potentially structured regions and subsequently predicts their functional importance. Using this procedure, the orthologous groups were split into structurally homogenous subgroups, which we call subVOGs. This is the first compilation of potentially functional conserved RNA structures in viral coding regions, covering the complete RefSeq viral database. We were able to recover structural elements from previous studies and discovered a variety of novel structured regions. The subVOGs are available through our web resource RNASIV (RNA structure in viruses; http://rnasiv.bio.wzw.tum.de).

Directions for research and training in plant omics: Big Questions and Big Data

A key remit of the NSF-funded "Arabidopsis Research and Training for the 21st Century" (ART-21) Research Coordination Network has been to convene a series of workshops with community members to explore issues concerning research and training in plant biology, including the role that research using Arabidopsis thaliana can play in addressing those issues. A first workshop focused on training needs for bioinformatic and computational approaches in plant biology was held in 2016, and recommendations from that workshop have been published (Friesner et al., Plant Physiology, 175, 2017, 1499). In this white paper, we provide a summary of the discussions and insights arising from the second ART-21 workshop. The second workshop focused on experimental aspects of omics data acquisition and analysis and involved a broad spectrum of participants from academics and industry, ranging from graduate students through post-doctorates, early career and established investigators. Our hope is that this article will inspire beginning and established scientists, corporations, and funding agencies to pursue directions in research and training identified by this workshop, capitalizing on the reference species Arabidopsis thaliana and other valuable plant systems.

A Model-Based Approach for Identifying Functional Intergenic Transcribed Regions and Noncoding RNAs

With advances in transcript profiling, the presence of transcriptional activities in intergenic regions has been well established. However, whether intergenic expression reflects transcriptional noise or activity of novel genes remains unclear. We identified intergenic transcribed regions (ITRs) in 15 diverse flowering plant species and found that the amount of intergenic expression correlates with genome size, a pattern that could be expected if intergenic expression is largely nonfunctional. To further assess the functionality of ITRs, we first built machine learning models using Arabidopsis thaliana as a model that accurately distinguish functional sequences (benchmark protein-coding and RNA genes) and likely nonfunctional ones (pseudogenes and unexpressed intergenic regions) by integrating 93 biochemical, evolutionary, and sequence-structure features. Next, by applying the models genome-wide, we found that 4,427 ITRs (38%) and 796 annotated ncRNAs (44%) had features significantly similar to benchmark protein-coding or RNA genes and thus were likely parts of functional genes. Approximately 60% of ITRs and ncRNAs were more similar to nonfunctional sequences and were likely transcriptional noise. The predictive framework established here provides not only a comprehensive look at how functional, genic sequences are distinct from likely nonfunctional ones, but also a new way to differentiate novel genes from genomic regions with noisy transcriptional activities.

Complex Relationships between Chromatin Accessibility, Sequence Divergence, and Gene Expression in Arabidopsis thaliana

Variation in regulatory DNA is thought to drive phenotypic variation, evolution, and disease. Prior studies of regulatory DNA and transcription factors across animal species highlighted a fundamental conundrum: Transcription factor binding domains and cognate binding sites are conserved, while regulatory DNA sequences are not. It remains unclear how conserved transcription factors and dynamic regulatory sites produce conserved expression patterns across species. Here, we explore regulatory DNA variation and its functional consequences within Arabidopsis thaliana, using chromatin accessibility to delineate regulatory DNA genome-wide. Unlike in previous cross-species comparisons, the positional homology of regulatory DNA is maintained among A. thaliana ecotypes and less nucleotide divergence has occurred. Of the ∼50,000 regulatory sites in A. thaliana, we found that 15% varied in accessibility among ecotypes. Some of these accessibility differences were associated with extensive, previously unannotated sequence variation, encompassing many deletions and ancient hypervariable alleles. Unexpectedly, for the majority of such regulatory sites, nearby gene expression was unaffected. Nevertheless, regulatory sites with high levels of sequence variation and differential chromatin accessibility were the most likely to be associated with differential gene expression. Finally, and most surprising, we found that the vast majority of differentially accessible sites show no underlying sequence variation. We argue that these surprising results highlight the necessity to consider higher-order regulatory context in evaluating regulatory variation and predicting its phenotypic consequences.

Beyond transcription factors: roles of mRNA decay in regulating gene expression in plants

Gene expression is typically quantified as RNA abundance, which is influenced by both synthesis (transcription) and decay. Cytoplasmic decay typically initiates by deadenylation, after which decay can occur through any of three cytoplasmic decay pathways. Recent advances reveal several mechanisms by which RNA decay is regulated to control RNA abundance. mRNA can be post-transcriptionally modified, either indirectly through secondary structure or through direct modifications to the transcript itself, sometimes resulting in subsequent changes in mRNA decay rates. mRNA abundances can also be modified by tapping into pathways normally used for RNA quality control. Regulated mRNA decay can also come about through post-translational modification of decapping complex subunits. Likewise, mRNAs can undergo changes in subcellular localization (for example, the deposition of specific mRNAs into processing bodies, or P-bodies, where stabilization and destabilization occur in a transcript- and context-dependent manner). Additionally, specialized functions of mRNA decay pathways were implicated in a genome-wide mRNA decay analysis in Arabidopsis. Advances made using plants are emphasized in this review, but relevant studies from other model systems that highlight RNA decay mechanisms that may also be conserved in plants are discussed.

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.

MicroRNAs from plants to animals, do they define a new messenger for communication?

MicroRNAs (miRNAs), a class of single-stranded non-coding RNA of about 22 nucleotides, are potent regulators of gene expression existing in both plants and animals. Recent studies showed that plant miRNAs could enter mammalian bloodstream via gastrointestinal tract, through which access a variety of tissues and cells of recipients to exert therapeutic effects. This intriguing phenomenon indicates that miRNAs of diet/plant origin may act as a new class of bioactive ingredients communicating with mammalian systems. In this review, in order to pinpoint the reason underlying discrepancies of miRNAs transmission from diet/plant to animals, the pathways that generate miRNAs and machineries involved in the functions of miRNAs in both kingdoms were outlined and compared. Then, the current controversies concerning cross-kingdom regulations and the potential mechanisms responsible for absorption and transfer of diet/plant-derived miRNAs were interpreted. Furthermore, the hormone-like action of miRNAs and the intricate interplay between miRNAs and hormones were implicated. Finally, how these findings may impact nutrition and medicine were briefly discussed.

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 spatial analysis of Arabidopsis natural variants implicates 5'UTR splicing of LATE ELONGATED HYPOCOTYL in responses to temperature

How plants perceive and respond to temperature remains an important question in the plant sciences. Temperature perception and signal transduction may occur through temperature-sensitive intramolecular folding of primary mRNA transcripts. Recent studies suggested a role for retention of the first intron in the 5'UTR of the clock component LATE ELONGATED HYPOCOTYL (LHY) in response to changes in temperature. Here, we identified a set of haplotypes in the LHY 5'UTR, examined their global spatial distribution, and obtained evidence that haplotype can affect temperature-dependent splicing of LHY transcripts. Correlations of haplotype spatial distributions with global bioclimatic variables and altitude point to associations with annual mean temperature and temperature fluctuation. Relatively rare relict type accessions correlate with lower mean temperature and greater temperature fluctuation and the spatial distribution of other haplotypes may be informative of evolutionary processes driving colonization of ecosystems. We propose that haplotypes may possess distinct 5'UTR pre-mRNA folding thermodynamics and/or specific biological stabilities based around the binding of trans-acting RNA splicing factors, a consequence of which is scalable splicing sensitivity of a central clock component that is likely tuned to specific temperature environments.

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.

Discovering Putative Peptides Encoded from Noncoding RNAs in Ribosome Profiling Data of Arabidopsis thaliana

Most noncoding RNAs are considered by their expression at low levels and as having a limited phylogenetic distribution in the cytoplasm, indicating that they may be only involved in specific biological processes. However, recent studies showed the protein-coding potential of ncRNAs, indicating that they might be a source of some special proteins. Although there are increasing noncoding RNAs identified to be able to code proteins, it is challenging to distinguish coding RNAs from previously annotated ncRNAs, and to detect the proteins from their translation. In this article, we designed a pipeline to identify these noncoding RNAs in Arabidopsis thaliana from three NCBI GEO data sets with coding potential and predict their translation products. 31 311 noncoding RNAs were predicted to be translated into peptides, and they showed lower conservation rate than common proteins. In addition, we built an interaction network between these peptides and annotated Arabidopsis proteins using BIPS, which included 69 peptides from noncoding RNAs. Peptides in the interaction network showed different characteristics from other noncoding RNA-derived peptides, and they participated in several crucial biological processes, such as photorespiration and stress-responses. All the information of putative ncPEPs and their interaction with proteins predicted above are finally integrated in a database, PncPEPDB ( http://bis.zju.edu.cn/PncPEPDB ). These results showed that peptides derived from noncoding RNAs may play important roles in noncoding RNA regulation, which provided another hypothesis that noncoding RNA may regulate the metabolism via their translation products.

RNA secondary structure profiling in zebrafish reveals unique regulatory features

Background

RNA is known to play diverse roles in gene regulation. The clues for this regulatory function of RNA are embedded in its ability to fold into intricate secondary and tertiary structure.

Results

We report the transcriptome-wide RNA secondary structure in zebrafish at single nucleotide resolution using Parallel Analysis of RNA Structure (PARS). This study provides the secondary structure map of zebrafish coding and non-coding RNAs. The single nucleotide pairing probabilities of 54,083 distinct transcripts in the zebrafish genome were documented. We identified RNA secondary structural features embedded in functional units of zebrafish mRNAs. Translation start and stop sites were demarcated by weak structural signals. The coding regions were characterized by the three-nucleotide periodicity of secondary structure and display a codon base specific structural constrain. The splice sites of transcripts were also delineated by distinct signature signals. Relatively higher structural signals were observed at 3’ Untranslated Regions (UTRs) compared to Coding DNA Sequence (CDS) and 5’ UTRs. The 3′ ends of transcripts were also marked by unique structure signals. Secondary structural signals in long non-coding RNAs were also explored to better understand their molecular function.

Conclusions

Our study presents the first PARS-enabled transcriptome-wide secondary structure map of zebrafish, which documents pairing probability of RNA at single nucleotide precision. Our findings open avenues for exploring structural features in zebrafish RNAs and their influence on gene expression.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4497-0) contains supplementary material, which is available to authorized users.

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.

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.

Target RNA Secondary Structure Is a Major Determinant of miR159 Efficacy

In plants, microRNA (miRNA)-target complementarity has long been considered the predominant factor determining the silencing outcome of the miRNA-target interaction, although the efficacy of such interactions have rarely been appraised in plants. Here, we perform in planta silencing efficacy assays on seven Arabidopsis MYB genes, all of which contain conserved miR159-binding sites of analogous complementarity. These genes were found to be differentially silenced by miR159; MYB81, MYB97, MYB101, MYB104, and DUO1 were all poorly silenced, whereas MYB33 and MYB65 were strongly silenced. Curiously, this is consistent with previous genetic analysis defining MYB33 and MYB65 as the major functional targets of miR159. Neither the free energy of miR159-target complementarity, nor miRNA binding site accessibility, as determined by flanking region AU content, could fully explain the discrepancy of miR159 silencing efficacy. Instead, we found that MYB33 and MYB65 were both predicted to contain a distinctive RNA secondary structure abutting the miR159 binding site. The structure is composed of two stem-loops (SLs) that are predicted to form in MYB33/65 homologs of species as evolutionary distant as gymnosperms. Functional analysis found that the RNA structure in MYB33 correlated with strong silencing efficacy; introducing mutations to disrupt either SL attenuated miR159 efficacy, while introducing complementary mutations to restore the SLs, but not the sequence, restored strong miR159-mediated silencing. Therefore, it appears that this RNA secondary structure demarcates MYB33/65 as sensitive targets of miR159, which underpins the narrow functional specificity of Arabidopsis miR159.

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.