Genome-wide measurement of RNA folding energies

Leveraging Saccharomyces cerevisiae for ADAR research: From high-yield purification to high-throughput screening and therapeutic applications

Saccharomyces cerevisiae, a model eukaryotic organism with a rich history in research and industry, has become a pivotal tool for studying Adenosine Deaminase Acting on RNA (ADAR) enzymes despite lacking these enzymes endogenously. This chapter reviews the diverse methodologies harnessed using yeast to elucidate ADAR structure and function, emphasizing its role in advancing our understanding of RNA editing. Initially, Saccharomyces cerevisiae was instrumental in the high-yield purification of ADARs, addressing challenges associated with enzyme stability and activity in other systems. The chapter highlights the successful application of yeast in high-throughput screening platforms that identify key structural motifs and substrate preferences of ADARs, showcasing its utility in revealing complex enzyme mechanics. Furthermore, we discuss the development of yeast-based systems to optimize guide RNA sequences for site-directed RNA editing (SDRE), demonstrating how these systems can be employed to refine therapeutic strategies targeting genetic mutations. Additionally, exogenous expression of ADARs from various species in yeast has shed light on enzyme potency and substrate recognition across different temperatures, offering insights into evolutionary adaptations. Overall, Saccharomyces cerevisiae has proven to be an invaluable asset in ADAR research, facilitating significant advances in our understanding of RNA editing mechanisms and therapeutic applications.

Unveiling RNA structure-mediated regulations of RNA stability in wheat

Despite the critical role of mRNA stability in post-transcriptional gene regulation, research on this topic in wheat, a vital agricultural crop, remains unclear. Our study investigated the mRNA decay landscape of durum wheat (Triticum turgidum L. ssp. durum, BBAA), revealing subgenomic asymmetry in mRNA stability and its impact on steady-state mRNA abundance. Our findings indicate that the 3' UTR structure and homoeolog preference for RNA structural motifs can influence mRNA stability, leading to subgenomic RNA decay imbalance. Furthermore, single-nucleotide variations (SNVs) selected for RNA structural motifs during domestication can cause variations in subgenomic mRNA stability and subsequent changes in steady-state expression levels. Our research on the transcriptome stability of polyploid wheat highlights the regulatory role of non-coding region structures in mRNA stability, and how domestication shaped RNA structure, altering subgenomic mRNA stability. These results illustrate the importance of RNA structure-mediated post-transcriptional gene regulation in wheat and pave the way for its potential use in crop improvement.

Exonuclease-assisted enrichment and base resolution analysis of pseudouridine in single-stranded RNA

Pseudouridine (Ψ) is one of the most abundant RNA modifications, playing crucial roles in various biological processes. Identifying Ψ sites is vital for understanding their functions. In this study, we proposed a novel method for identifying Ψ sites with an improved signal-to-noise ratio. This method, called RNA exonuclease-assisted identification of pseudouridine sites (RIPS), combines specific CMC-labeling of Ψ sites with an exonuclease-assisted digestion strategy for the detection of Ψ sites. Utilizing exonuclease XRN1 to digest RNA strands not labeled by CMC, RIPS significantly reduces the background signal from unlabeled strands and enhances the positive signal of Ψ sites labeled by CMC, which terminates exonuclease digestion. As a result, we can enrich Ψ sites and identify them at single-base resolution. Considering the unique functions of single-stranded RNA (ssRNA), we employed RIPS to distinguish Ψ sites in single-stranded and double-stranded regions of RNA. Our results indicated that CMC could specifically label Ψ sites in ssRNA under natural conditions, enabling RIPS to selectively identify Ψ sites in ssRNA, which may facilitate the study on the functions of Ψ sites.

Technological advancements in deciphering RNA-RNA interactions

RNA-RNA interactions (RRIs) can dictate RNA molecules to form intricate higher-order structures and bind their RNA substrates in diverse biological processes. To elucidate the function, binding specificity, and regulatory mechanisms of various RNA molecules, especially the vast repertoire of non-coding RNAs, advanced technologies and methods that globally map RRIs are extremely valuable. In the past decades, many state-of-the-art technologies have been developed for this purpose. This review focuses on those high-throughput technologies for the global mapping of RRIs. We summarize the key concepts and the pros and cons of different technologies. In addition, we highlight the novel biological insights uncovered by these RRI mapping methods and discuss the future challenges for appreciating the crucial roles of RRIs in gene regulation across bacteria, viruses, archaea, and mammals.

Dynamic structural remodeling of LINC01956 enhances temozolomide resistance in MGMT-methylated glioblastoma

The mechanisms underlying stimuli-induced dynamic structural remodeling of RNAs for the maintenance of cellular physiological function and survival remain unclear. Here, we showed that in MGMT promoter-methylated glioblastoma (GBM), the RNA helicase DEAD-box helicase 46 (DDX46) is phosphorylated by temozolomide (TMZ)-activated checkpoint kinase 1 (CHK1), resulting in a dense-to-loose conformational change and an increase in DDX46 helicase activity. DDX46-mediated tertiary structural remodeling of LINC01956 exposes the binding motifs of LINC01956 to the 3' untranslated region of O6-methylguanine DNA methyltransferase (MGMT). This accelerates recruitment of MGMT mRNA to the RNA export machinery and transportation of MGMT mRNA from the nucleus to the cytoplasm, leading to increased MGMT abundance and TMZ resistance. Using patient-derived xenograft (PDX) and tumor organoid models, we found that treatment with the CHK1 inhibitor SRA737abolishes TMZ-induced structural remodeling of LINC01956 and subsequent MGMT up-regulation, resensitizing TMZ-resistant MGMT promoter-methylated GBM to TMZ. In conclusion, these findings highlight a mechanism underlying temozolomide-induced RNA structural remodeling and may represent a potential therapeutic strategy for patients with TMZ-resistant MGMT promoter-methylated GBM.

Identification of RNA structures and their roles in RNA functions

The development of high-throughput RNA structure profiling methods in the past decade has greatly facilitated our ability to map and characterize different aspects of RNA structures transcriptome-wide in cell populations, single cells and single molecules. The resulting high-resolution data have provided insights into the static and dynamic nature of RNA structures, revealing their complexity as they perform their respective functions in the cell. In this Review, we discuss recent technical advances in the determination of RNA structures, and the roles of RNA structures in RNA biogenesis and functions, including in transcription, processing, translation, degradation, localization and RNA structure-dependent condensates. We also discuss the current understanding of how RNA structures could guide drug design for treating genetic diseases and battling pathogenic viruses, and highlight existing challenges and future directions in RNA structure research.

Association Analysis of METTL23 Gene Polymorphisms with Reproductive Traits in Kele Pigs

Methyltransferase-like 23 (METTL23) is a kind of RNA methyltransferase that catalyzes the methylation transfer to the N6-adenosine of RNA, serving as one of the key mediators in this process. However, the METTL23 gene has been poorly researched in pigs. In this study, we investigated the genetic effects of METTL23 single-nucleotide polymorphism(SNPs) on reproductive traits in Kele pigs. The DNA was extracted from 228 healthy multiparous Kele sows, and Sanger sequencing revealed three SNPs, g.4804958 G > T (intron 2), g.4805082 C > T (exon 2), and g.4806821 A > G (exon 3). The polymorphism information content (PIC) for each SNP was 0.264, 0.25, and 0.354, indicating moderate polymorphism (0.25 < PIC < 0.5) and providing genetic information. Linkage disequilibrium analysis showed no strong linkage disequilibrium between the three SNPs. The association analysis revealed that in the SNP g.4804958 G > T individuals with the GG genotype had a significantly higher number of piglets born alive, litter birth weight, number of weaned piglets, and weaning litter weight compared to those with the TT genotype (p < 0.05). Individuals with the GG genotype in the SNP g.4806821 A > G group had significantly higher litter birth weight and average birth weight than those with the AA genotype (p < 0.05). The H4H4 diplotype showed significant effects on the number of piglets born alive, litter birth weight, number of weaned piglets, weaning litter weight, and weaning weight (p < 0.05). Together, the METTL23 gene could be used as a candidate gene for the selection of reproductive traits in Kele pigs.

Revealing Differential RNA Editing Specificity of Human ADAR1 and ADAR2 in Schizosaccharomyces pombe

Adenosine-to-inosine (A-to-I) RNA editing is an important post-transcriptional modification mediated by the adenosine deaminases acting on RNA (ADAR) family of enzymes, expanding the transcriptome by altering selected nucleotides A to I in RNA molecules. Recently, A-to-I editing has been explored for correcting disease-causing mutations in RNA using therapeutic guide oligonucleotides to direct ADAR editing at specific sites. Humans have two active ADARs whose preferences and specificities are not well understood. To investigate their substrate specificity, we introduced hADAR1 and hADAR2, respectively, into Schizosaccharomyces pombe (S. pombe), which lacks endogenous ADARs, and evaluated their editing activities in vivo. Using transcriptome sequencing of S. pombe cultured at optimal growth temperature (30 °C), we identified 483 A-to-I high-confident editing sites for hADAR1 and 404 for hADAR2, compared with the non-editing wild-type control strain. However, these sites were mostly divergent between hADAR1 and hADAR2-expressing strains, sharing 33 common sites that are less than 9% for each strain. Their differential specificity for substrates was attributed to their differential preference for neighboring sequences of editing sites. We found that at the -3-position relative to the editing site, hADAR1 exhibits a tendency toward T, whereas hADAR2 leans toward A. Additionally, when varying the growth temperature for hADAR1- and hADAR2-expressing strains, we observed increased editing sites for them at both 20 and 35 °C, compared with them growing at 30 °C. However, we did not observe a significant shift in hADAR1 and hADAR2's preference for neighboring sequences across three temperatures. The vast changes in RNA editing sites at lower and higher temperatures were also observed for hADAR2 previously in budding yeast, which was likely due to the influence of RNA folding at these different temperatures, among many other factors. We noticed examples of longer lengths of dsRNA around the editing sites that induced editing at 20 or 35 °C but were absent at the other two temperature conditions. We found genes' functions can be greatly affected by editing of their transcripts, for which over 50% of RNA editing sites for both hADAR1 and hADAR2 in S. pombe were in coding sequences (CDS), with more than 60% of them resulting in amino acid changes in protein products. This study revealed the extensive differences in substrate selectivity between the two active human ADARS, i.e., ADAR1 and ADAR2, and provided novel insight when utilizing the two different enzymes for in vivo treatment of human genetic diseases using the RNA editing approach.

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.

Single-molecule quantitation of RNA-binding protein occupancy and stoichiometry defines a role for Yra1 (Aly/REF) in nuclear mRNP organization

RNA-binding proteins (RBPs) interact with mRNA to form supramolecular complexes called messenger ribonucleoprotein (mRNP) particles. These dynamic assemblies direct and regulate individual steps of gene expression; however, their composition and functional importance remain largely unknown. Here, we develop a total internal reflection fluorescence-based single-molecule imaging assay to investigate stoichiometry and co-occupancy of 15 RBPs within mRNPs from Saccharomyces cerevisiae. We show compositional heterogeneity of single mRNPs and plasticity across different growth conditions, with major co-occupants of mRNPs containing the nuclear cap-binding complex identified as Yra1 (1-10 copies), Nab2 (1-6 copies), and Npl3 (1-6 copies). Multicopy Yra1-bound mRNPs are specifically co-occupied by the THO complex and assembled on mRNAs biased by transcript length and RNA secondary structure. Yra1 depletion results in decreased compaction of nuclear mRNPs demonstrating a packaging function. Together, we provide a quantitative framework for gene- and condition-dependent RBP occupancy and stoichiometry in individual nuclear mRNPs.

Tissue-specific temperature dependence of RNA editing levels in zebrafish

BACKGROUND

RNA editing by adenosine deaminase acting on RNA (ADAR) occurs in all metazoans and fulfils several functions. Here, we examined effects of acclimation temperature (27 °C, 18 °C,13 °C) on editing patterns in six tissues of zebrafish (Danio rerio).

RESULTS

Sites and total amounts of editing differed among tissues. Brain showed the highest levels, followed by gill and skin. In these highly edited tissues, decreases in temperatures led to large increases in total amounts of editing and changes in specific edited sites. Gene ontology analysis showed both similarities (e.g., endoplasmic reticulum stress response) and differences in editing among tissues. The majority of edited sites were in transcripts of transposable elements and the 3'UTR regions of protein coding genes. By experimental validation, translation efficiency was directly related to extent of editing of the 3'UTR region of an mRNA.

CONCLUSIONS

RNA editing increases 3'UTR polymorphism and affects efficiency of translation. Such editing may lead to temperature-adaptive changes in the proteome through altering relative amounts of synthesis of different proteins.

Incorporation of noncanonical base Z yields modified mRNA with minimal immunogenicity and improved translational capacity in mammalian cells

Chemically modified mRNAs hold great potential for therapeutic applications in vivo. Currently, the base modification scheme largely preserves the canonical Watson-Crick base pairing, thus missing one mode of mRNA modulation by altering its secondary structure. Here we report the incorporation of base Z (2-aminoadenine) into mRNA to create Z-mRNA with improved translational capacity, decreased cytotoxicity, and drastically reduced immunogenicity compared to the unmodified mRNA in mammalian cells. In particular, the A-to-Z substitution renders modified mRNAs less immunogenic than the state-of-the-art base modification N1-methylpseudouridine (m1ψ) in mouse embryonic fibroblast cells. As a proof of concept, we developed a Z-mRNA-based vaccine against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Antigen-encoding Z-mRNA elicited substantial humoral and cellular immune responses in vivo in mice, albeit with relatively lower efficacy than the state-of-the-art m1ψ-mRNA. Z-mRNA expands the scope of mRNA base modifications toward noncanonical bases and could offer an advantageous platform for mRNA-based therapeutics where minimal immunogenicity is desired.

RNA-RNA competitive interactions: a molecular civil war ruling cell physiology and diseases

The idea that proteins are the main determining factors in the functioning of cells and organisms, and their dysfunctions are the first cause of pathologies, has been predominant in biology and biomedicine until recently. This protein-centered view was too simplistic and failed to explain the physiological and pathological complexity of the cell. About 80% of the human genome is dynamically and pervasively transcribed, mostly as non-protein-coding RNAs (ncRNAs), which competitively interact with each other and with coding RNAs generating a complex RNA network regulating RNA processing, stability, and translation and, accordingly, fine-tuning the gene expression of the cells. Qualitative and quantitative dysregulations of RNA-RNA interaction networks are strongly involved in the onset and progression of many pathologies, including cancers and degenerative diseases. This review will summarize the RNA species involved in the competitive endogenous RNA network, their mechanisms of action, and involvement in pathological phenotypes. Moreover, it will give an overview of the most advanced experimental and computational methods to dissect and rebuild RNA networks.

Natural variation in codon bias and mRNA folding strength interact synergistically to modify protein expression in Saccharomyces cerevisiae

Codon bias and mRNA folding strength (mF) are hypothesized molecular mechanisms by which polymorphisms in genes modify protein expression. Natural patterns of codon bias and mF across genes as well as effects of altering codon bias and mF suggest that the influence of these 2 mechanisms may vary depending on the specific location of polymorphisms within a transcript. Despite the central role codon bias and mF may play in natural trait variation within populations, systematic studies of how polymorphic codon bias and mF relate to protein expression variation are lacking. To address this need, we analyzed genomic, transcriptomic, and proteomic data for 22 Saccharomyces cerevisiae isolates, estimated protein accumulation for each allele of 1,620 genes as the log of protein molecules per RNA molecule (logPPR), and built linear mixed-effects models associating allelic variation in codon bias and mF with allelic variation in logPPR. We found that codon bias and mF interact synergistically in a positive association with logPPR, and this interaction explains almost all the effects of codon bias and mF. We examined how the locations of polymorphisms within transcripts influence their effects and found that codon bias primarily acts through polymorphisms in domain-encoding and 3' coding sequences, while mF acts most significantly through coding sequences with weaker effects from untranslated regions. Our results present the most comprehensive characterization to date of how polymorphisms in transcripts influence protein expression.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

Thermomorphogenesis: Opportunities and challenges in posttranscriptional regulation

Plants exposed to mildly elevated temperatures display morphological and developmental changes collectively termed thermomorphogenesis. This adaptative process has several undesirable consequences to food production, including yield reduction and increased vulnerability to pathogens. Understanding thermomorphogenesis is, thus, critical for understanding how plants will respond to increasingly warmer temperature conditions, such as those caused by climate change. Recently, we have made major advances in that direction, and it has become apparent that plants resource to a broad range of molecules and molecular mechanisms to perceive and respond to increases in environmental temperature. However, most of our efforts have been focused on regulation of transcription and protein abundance and activity, with an important gap encompassing nearly all processes involving RNA (i.e., posttranscriptional regulation). Here, I summarized our current knowledge of thermomorphogenesis involving transcriptional, posttranscriptional, and posttranslational regulation, focused on opportunities and challenges in understanding posttranscriptional regulation-a fertile field for exciting new discoveries.

Identification of exceptionally potent adenosine deaminases RNA editors from high body temperature organisms

The most abundant form of RNA editing in metazoa is the deamination of adenosines into inosines (A-to-I), catalyzed by ADAR enzymes. Inosines are read as guanosines by the translation machinery, and thus A-to-I may lead to protein recoding. The ability of ADARs to recode at the mRNA level makes them attractive therapeutic tools. Several approaches for Site-Directed RNA Editing (SDRE) are currently under development. A major challenge in this field is achieving high on-target editing efficiency, and thus it is of much interest to identify highly potent ADARs. To address this, we used the baker yeast Saccharomyces cerevisiae as an editing-naïve system. We exogenously expressed a range of heterologous ADARs and identified the hummingbird and primarily mallard-duck ADARs, which evolved at 40-42°C, as two exceptionally potent editors. ADARs bind to double-stranded RNA structures (dsRNAs), which in turn are temperature sensitive. Our results indicate that species evolved to live with higher core body temperatures have developed ADAR enzymes that target weaker dsRNA structures and would therefore be more effective than other ADARs. Further studies may use this approach to isolate additional ADARs with an editing profile of choice to meet specific requirements, thus broadening the applicability of SDRE.

Genome-wide sequencing identifies a thermal-tolerance related synonymous mutation in the mussel, Mytilisepta virgata

The roles of synonymous mutations for adapting to stressful thermal environments are of fundamental biological and ecological interests but poorly understood. To study whether synonymous mutations influence thermal adaptation at specific microhabitats, a genome-wide genotype-phenotype association analysis is carried out in the black mussels Mytilisepta virgata. A synonymous mutation of Ubiquitin-specific Peptidase 15 (MvUSP15) is significantly associated with the physiological upper thermal limit. The individuals carrying GG genotype (the G-type) at the mutant locus possess significantly lower heat tolerance compared to the individuals carrying GA and AA genotypes (the A-type). When heated to sublethal temperature, the G-type exhibit higher inter-individual variations in MvUSP15 expression, especially for the mussels on the sun-exposed microhabitats. Taken together, a synonymous mutation in MvUSP15 can affect the gene expression profile and interact with microhabitat heterogeneity to influence thermal resistance. This integrative study sheds light on the ecological importance of adaptive synonymous mutations as an underappreciated genetic buffer against heat stress and emphasizes the importance of integrative studies at a microhabitat scale for evaluating and predicting the impacts of climate change.

Detecting RNA-RNA interactome

The last decade has seen a robust increase in various types of novel RNA molecules and their complexity in gene regulation. RNA molecules play a critical role in cellular events by interacting with other biomolecules, including protein, DNA, and RNA. It has been established that RNA-RNA interactions play a critical role in several biological processes by regulating the biogenesis and function of RNA molecules. Interestingly, RNA-RNA interactions regulate the biogenesis of diverse RNA molecules, including mRNAs, microRNAs, tRNAs, and circRNAs, through splicing or backsplicing. Structured RNAs like rRNA, tRNA, and snRNAs achieve their functional conformation by intramolecular RNA-RNA interactions. In addition, functional consequences of many intermolecular RNA-RNA interactions have been extensively studied in the regulation of gene expression. Hence, it is essential to understand the mechanism and functions of RNA-RNA interactions in eukaryotes. Conventionally, RNA-RNA interactions have been identified through diverse biochemical methods for decades. The advent of high-throughput RNA-sequencing technologies has revolutionized the identification of global RNA-RNA interactome in cells and their importance in RNA structure and function in gene expression regulation. Although these technologies revealed tens of thousands of intramolecular and intermolecular RNA-RNA interactions, we further look forward to future unbiased and quantitative high-throughput technologies for detecting transcriptome-wide RNA-RNA interactions. With the ability to detect RNA-RNA interactome, we expect that future studies will reveal the higher-order structures of RNA molecules and multi-RNA hybrids impacting human health and diseases. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Response and regulatory mechanisms of heat resistance in pathogenic fungi

Abstract

Both the increasing environmental temperature in nature and the defensive body temperature response to pathogenic fungi during mammalian infection cause heat stress during the fungal existence, reproduction, and pathogenic infection. To adapt and respond to the changing environment, fungi initiate a series of actions through a perfect thermal response system, conservative signaling pathways, corresponding transcriptional regulatory system, corresponding physiological and biochemical processes, and phenotypic changes. However, until now, accurate response and regulatory mechanisms have remained a challenge. Additionally, at present, the latest research progress on the heat resistance mechanism of pathogenic fungi has not been summarized. In this review, recent research investigating temperature sensing, transcriptional regulation, and physiological, biochemical, and morphological responses of fungi in response to heat stress is discussed. Moreover, the specificity thermal adaptation mechanism of pathogenic fungi in vivo is highlighted. These data will provide valuable knowledge to further understand the fungal heat adaptation and response mechanism, especially in pathogenic heat-resistant fungi.

Key points

• Mechanisms of fungal perception of heat pressure are reviewed.

• The regulatory mechanism of fungal resistance to heat stress is discussed.

• The thermal adaptation mechanism of pathogenic fungi in the human body is highlighted.

The Bacillaceae-1 RNA motif comprises two distinct classes

Non-coding RNAs are key regulatory players in bacteria. Many computationally predicted non-coding RNAs, however, lack functional associations. An example is the Bacillaceae-1 RNA motif, whose Rfam model consists of two hairpin loops. We find the motif conserved in nine of 13 non-pathogenic strains of the genus Bacillus but only in one pathogenic strain. To elucidate functional characteristics, we studied 118 hits of the Rfam model in 11 Bacillus spp. and found two distinct classes based on the ensemble diversity of their RNA secondary structure and the genomic context concerning the ribosomal RNA (rRNA) cluster. Forty hits are associated with the rRNA cluster, of which all 19 hits upstream flanking of 16S rRNA have a reverse complementary structure of low structural diversity. Fifty-two hits have large ensemble diversity, of which 38 are located between two coding genes. For eight hits in Bacillus subtilis, we investigated public expression data under various conditions and observed either the forward or the reverse complementary motif expressed. Five hits are associated with the rRNA cluster. Four of them are located upstream of the 16S rRNA and are not transcriptionally active, but instead, their reverse complements with low structural diversity are expressed together with the rRNA cluster. The three other hits are located between two coding genes in non-conserved genomic loci. Two of them are independently expressed from their surrounding genes and are structurally diverse. In summary, we found that Bacillaceae-1 RNA motifs upstream flanking of ribosomal RNA clusters tend to have one stable structure with the reverse complementary motif expressed in B. subtilis. In contrast, a subgroup of intergenic motifs has the thermodynamic potential for structural switches.

Advances and Trends in Omics Technology Development

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

Diffuse Ions Coordinate Dynamics in a Ribonucleoprotein Assembly

Proper ionic concentrations are required for the functional dynamics of RNA and ribonucleoprotein (RNP) assemblies. While experimental and computational techniques have provided many insights into the properties of chelated ions, less is known about the energetic contributions of diffuse ions to large-scale conformational rearrangements. To address this, we present a model that is designed to quantify the influence of diffuse monovalent and divalent ions on the dynamics of biomolecular assemblies. This model employs all-atom (non-H) resolution and explicit ions, where effective potentials account for hydration effects. We first show that the model accurately predicts the number of excess Mg²⁺ ions for prototypical RNA systems, at a level comparable to modern coarse-grained models. We then apply the model to a complete ribosome and show how the balance between diffuse Mg²⁺ and K⁺ ions can control the dynamics of tRNA molecules during translation. The model predicts differential effects of diffuse ions on the free-energy barrier associated with tRNA entry and the energy of tRNA binding to the ribosome. Together, this analysis reveals the direct impact of diffuse ions on the dynamics of an RNP assembly.

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.

mSRFR: a machine learning model using microalgal signature features for ncRNA classification

This work presents mSRFR (microalgae SMOTE Random Forest Relief model), a classification tool for noncoding RNAs (ncRNAs) in microalgae, including green algae, diatoms, golden algae, and cyanobacteria. First, the SMOTE technique was applied to address the challenge of imbalanced data due to the different numbers of microalgae ncRNAs from different species in the EBI RNA-central database. Then the top 20 significant features from a total of 106 features, including sequence-based, secondary structure, base-pair, and triplet sequence-structure features, were selected using the Relief feature selection method. Next, ten-fold cross-validation was applied to choose a classifier algorithm with the highest performance among Support Vector Machine, Random Forest, Decision Tree, Naïve Bayes, K-nearest Neighbor, and Neural Network, based on the receiver operating characteristic (ROC) area. The results showed that the Random Forest classifier achieved the highest ROC area of 0.992. Then, the Random Forest algorithm was selected and compared with other tools, including RNAcon, CPC, CPC2, CNCI, and CPPred. Our model achieved a high accuracy of about 97% and a low false-positive rate of about 2% in predicting the test dataset of microalgae. Furthermore, the top features from Relief revealed that the %GA dinucleotide is a signature feature of microalgal ncRNAs when compared to Escherichia coli, Saccharomyces cerevisiae, Arabidopsis thaliana, and Homo sapiens.

E2F4 may be a core transcription factor in the lncRNA-TF regulatory network in cervical cancer

Background

Cervical cancer is the most common gynecological cancer worldwide and is associated with high morbidity and mortality. Despite improvements in therapeutic strategies, the network regulation mechanism remains unclear and the treatment effect is not satisfactory. Therefore, there is a need to continue studying the mechanism of cervical cancer to explore effective gene targets and precise targeted therapy drugs.

Methods

First, three paired tissues (cancer tissues and noncancerous tissues) from patients with cervical squamous cell carcinoma were collected, grouped, and analyzed by microarray. Second, differentially expressed mRNAs (DEMs) and differentially expressed lncRNAs (DELs) (|fold change| ≥ 2 and p < 0.05) between the two groups were screened. For DEMs, functional annotation and pathway analysis were performed using DAVID. Functional prediction of DELs was then performed and their cis‐regulatory and trans‐regulatory networks were explored.

Results

Function prediction of DELs (both up‐regulated and down‐regulated) shows that the highest frequency Cellular Component (CC) item is cytosol, the highest frequency Molecular function (MF) item is mitotic cell cycle and the highest frequency Biological Process (BP) item is protein binding. Through cis‐regulation analysis of DELs, the cis‐regulatory relationship of 96 DELs was predicted. The lncRNA‐trans‐regulation network analysis suggested that E2F4 may be the core transcription factor in the lncRNA‐TF regulatory network in cervical cancer.

Conclusions

The lncRNA‐TF regulatory network plays an important role in the occurrence and progression of cervical cancer, and E2F4 may be a critical transcription factor in the regulatory network.

RNA thermometers and other regulatory elements: Diversity and importance in bacterial pathogenesis

Survival of microorganisms depends to a large extent on environmental conditions and the occupied host. By adopting specific strategies, microorganisms can thrive in the surrounding environment and, at the same time, preserve their viability. Evading the host defenses requires several mechanisms compatible with the host survival which include the production of RNA thermometers to regulate the expression of genes responsible for heat or cold shock as well as of those involved in virulence. Microorganisms have developed a variety of molecules in response to the environmental changes in temperature and even more specifically to the host they invade. Among all, RNA-based regulatory mechanisms are the most common ones, highlighting the importance of such molecules in gene expression control and novel drug development by suitable structure-based alterations. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA in Disease and Development > RNA in Disease RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

An ultra low-input method for global RNA structure probing uncovers Regnase-1-mediated regulation in macrophages

To enable diverse functions and precise regulation, an RNA sequence often folds into complex yet distinct structures in different cellular states. Probing RNA in its native environment is essential to uncovering RNA structures of biological contexts. However, current methods generally require large amounts of input RNA and are challenging for physiologically relevant use. Here, we report smartSHAPE, a new RNA structure probing method that requires very low amounts of RNA input due to the largely reduced artefact of probing signals and increased efficiency of library construction. Using smartSHAPE, we showcased the profiling of the RNA structure landscape of mouse intestinal macrophages upon inflammation, and provided evidence that RNA conformational changes regulate immune responses. These results demonstrate that smartSHAPE can greatly expand the scope of RNA structure-based investigations in practical biological systems, and also provide a research paradigm for the study of post-transcriptional regulation.

Thermal adaptation of mRNA secondary structure: stability versus lability

Macromolecular function commonly involves rapidly reversible alterations in three-dimensional structure (conformation). To allow these essential conformational changes, macromolecules must possess higher order structures that are appropriately balanced between rigidity and flexibility. Because of the low stabilization free energies (marginal stabilities) of macromolecule conformations, temperature changes have strong effects on conformation and, thereby, on function. As is well known for proteins, during evolution, temperature-adaptive changes in sequence foster retention of optimal marginal stability at a species' normal physiological temperatures. Here, we extend this type of analysis to messenger RNAs (mRNAs), a class of macromolecules for which the stability-lability balance has not been elucidated. We employ in silico methods to determine secondary structures and estimate changes in free energy of folding (ΔG_fold) for 25 orthologous mRNAs that encode the enzyme cytosolic malate dehydrogenase in marine mollusks with adaptation temperatures spanning an almost 60 °C range. The change in free energy that occurs during formation of the ensemble of mRNA secondary structures is significantly correlated with adaptation temperature: ΔG_fold values are all negative and their absolute values increase with adaptation temperature. A principal mechanism underlying these adaptations is a significant increase in synonymous guanine + cytosine substitutions with increasing temperature. These findings open up an avenue of exploration in molecular evolution and raise interesting questions about the interaction between temperature-adaptive changes in mRNA sequence and in the proteins they encode.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

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

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

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.

Structural characterization of NORAD reveals a stabilizing role of spacers and two new repeat units

Long non-coding RNAs (lncRNAs) can perform a variety of key cellular functions by interacting with proteins and other RNAs. Recent studies have shown that the functions of lncRNAS are largely mediated by their structures. However, our structural knowledge for most lncRNAS is limited to sequence-based computational predictions. Non-coding RNA activated by DNA damage (NORAD) is an atypical lncRNA due to its abundant expression and high sequence conservation. NORAD regulates genomic stability by interacting with proteins and microRNAs. Previous sequence-based characterization has identified a modular organization of NORAD composed of several NORAD repeat units (NRUs). These units comprise the protein-binding elements and are separated by regular spacers. Here, we experimentally determine for the first time the secondary structure of NORAD using the nextPARS approach. Our results suggest that the spacer regions provide structural stability to NRUs. Furthermore, we uncover two previously unreported NRUs, and determine the core structural motifs conserved across NRUs. Overall, these findings will help to elucidate the function and evolution of NORAD.

Decoding LncRNAs

Non-coding RNAs (ncRNAs) have been considered as unimportant additions to the transcriptome. Yet, in light of numerous studies, it has become clear that ncRNAs play important roles in development, health and disease. Long-ignored, long non-coding RNAs (lncRNAs), ncRNAs made of more than 200 nucleotides have gained attention due to their involvement as drivers or suppressors of a myriad of tumours. The detailed understanding of some of their functions, structures and interactomes has been the result of interdisciplinary efforts, as in many cases, new methods need to be created or adapted to characterise these molecules. Unlike most reviews on lncRNAs, we summarize the achievements on lncRNA studies by taking into consideration the approaches for identification of lncRNA functions, interactomes, and structural arrangements. We also provide information about the recent data on the involvement of lncRNAs in diseases and present applications of these molecules, especially in medicine.

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.

Synonymous variants that disrupt messenger RNA structure are significantly constrained in the human population

BACKGROUND

The role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of "silent" genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking.

RESULTS

We developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population.

CONCLUSIONS

These findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a "Structural Predictivity Index" (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

Transcription Regulation Through Nascent RNA Folding

Folding of RNA into secondary structures through intramolecular base pairing determines an RNA's three-dimensional architecture and associated function. Simple RNA structures like stem loops can provide specialized functions independent of coding capacity, such as protein binding, regulation of RNA processing and stability, stimulation or inhibition of translation. RNA catalysis is dependent on tertiary structures found in the ribosome, tRNAs and group I and II introns. While the extent to which non-coding RNAs contribute to cellular maintenance is generally appreciated, the fact that both non-coding and coding RNA can assume relevant structural states has only recently gained attention. In particular, the co-transcriptional folding of nascent RNA of all classes has the potential to regulate co-transcriptional processing, RNP (ribonucleoprotein particle) formation, and transcription itself. Riboswitches are established examples of co-transcriptionally folded coding RNAs that directly regulate transcription, mainly in prokaryotes. Here we discuss recent studies in both prokaryotes and eukaryotes showing that structure formation may carry a more widespread regulatory logic during RNA synthesis. Local structures forming close to the catalytic center of RNA polymerases have the potential to regulate transcription by reducing backtracking. In addition, stem loops or more complex structures may alter co-transcriptional RNA processing or its efficiency. Several examples of functional structures have been identified to date, and this review provides an overview of physiologically distinct processes where co-transcriptionally folded RNA plays a role. Experimental approaches such as single-molecule FRET and in vivo structural probing to further advance our insight into the significance of co-transcriptional structure formation are discussed.

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.

Specificity of RNA Folding and Its Association with Evolutionarily Adaptive mRNA Secondary Structures

The secondary structure is a fundamental feature of both noncoding and messenger RNAs. However, our understanding of the secondary structure of mRNA, especially that of the coding regions, remains elusive, likely due to translation and the lack of RNA-binding proteins that sustain the consensus structure, such as those that bind to noncoding RNA. Indeed, mRNA has recently been found to adopt diverse alternative structures, the overall functional significance of which remains untested. We hereby approached this problem by estimating the folding specificity, i.e., the probability that a fragment of RNA folds back to the same partner once refolded. We showed that the folding specificity of mRNA is lower than that of noncoding RNA and exhibits moderate evolutionary conservation. Notably, we found that specific rather than alternative folding is likely evolutionarily adaptive since specific folding is frequently associated with functionally important genes or sites within a gene. Additional analysis in combination with ribosome density suggests the ability to modulate ribosome movement as one potential functional advantage provided by specific folding. Our findings revealed a novel facet of the RNA structurome with important functional and evolutionary implications and indicated a potential method for distinguishing the mRNA secondary structures maintained by natural selection from molecular noise.

Analytical ultracentrifuge: an ideal tool for characterization of non-coding RNAs

Analytical ultracentrifugation (AUC) has emerged as a robust and reliable technique for biomolecular characterization with extraordinary sensitivity. AUC is widely used to study purity, conformational changes, biomolecular interactions, and stoichiometry. Furthermore, AUC is used to determine the molecular weight of biomolecules such as proteins, carbohydrates, and DNA and RNA. Due to the multifaceted role(s) of non-coding RNAs from viruses, prokaryotes, and eukaryotes, research aimed at understanding the structure-function relationships of non-coding RNAs is rapidly increasing. However, due to their large size, flexibility, complicated secondary structures, and conformations, structural studies of non-coding RNAs are challenging. In this review, we are summarizing the application of AUC to evaluate the homogeneity, interactions, and conformational changes of non-coding RNAs from adenovirus as well as from Murray Valley, Powassan, and West Nile viruses. We also discuss the application of AUC to characterize eukaryotic long non-coding RNAs, Xist, and HOTAIR. These examples highlight the significant role AUC can play in facilitating the structural determination of non-coding RNAs and their complexes.

The molecular structure of long non-coding RNAs: emerging patterns and functional implications

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts that regulate vital cellular processes and are crucially connected to diseases. Despite their unprecedented molecular complexity, it is emerging that lncRNAs possess distinct structural motifs. Remarkably, the 3D shape and topology of full-length, native lncRNAs have been visualized for the first time in the last year. These studies reveal that lncRNA structures dictate lncRNA functions. Here, we review experimentally determined lncRNA structures and emphasize that lncRNA structural characterization requires synergistic integration of computational, biochemical and biophysical approaches. Based on these emerging paradigms, we discuss how to overcome the challenges posed by the complex molecular architecture of lncRNAs, with the goal of obtaining a detailed understanding of lncRNA functions and molecular mechanisms in the future.

Group II intron as cold sensor for self-preservation and bacterial conjugation

Group II introns are self-splicing ribozymes and mobile genetic elements. Splicing is required for both expression of the interrupted host gene and intron retromobility. For the pRS01 plasmid-encoded Lactococcus lactis group II intron, Ll.LtrB, splicing enables expression of the intron's host relaxase protein. Relaxase, in turn, initiates horizontal transfer of the conjugative pRS01 plasmid and stimulates retrotransposition of the intron. Little is known about how splicing of bacterial group II introns is influenced by environmental conditions. Here, we show that low temperatures can inhibit Ll.LtrB intron splicing. Whereas autocatalysis is abolished in the cold, splicing is partially restored by the intron-encoded protein (IEP). Structure profiling reveals cold-induced disruptions of key tertiary interactions, suggesting that a kinetic trap prevents the intron RNA from assuming its native state. Interestingly, while reduced levels of transcription and splicing lead to a paucity of excised intron in the cold, levels of relaxase mRNA are maintained, partially due to diminished intron-mediated mRNA targeting, allowing intron spread by conjugal transfer. Taken together, this study demonstrates not only the intrinsic cold sensitivity of group II intron splicing and the role of the IEP for cold-stress adaptation, but also maintenance of horizontal plasmid and intron transfer under cold-shock.

Identification and analysis of RNA structural disruptions induced by single nucleotide variants using Riprap and RiboSNitchDB

RNA conformational alteration has significant impacts on cellular processes and phenotypic variations. An emerging genetic factor of RNA conformational alteration is a new class of single nucleotide variant (SNV) named riboSNitch. RiboSNitches have been demonstrated to be involved in many genetic diseases. However, identifying riboSNitches is notably difficult as the signals of RNA structural disruption are often subtle. Here, we introduce a novel computational framework–RIboSNitch Predictor based on Robust Analysis of Pairing probabilities (Riprap). Riprap identifies structurally disrupted regions around any given SNVs based on robust analysis of local structural configurations between wild-type and mutant RNA sequences. Compared to previous approaches, Riprap shows higher accuracy when assessed on hundreds of known riboSNitches captured by various experimental RNA structure probing methods including the parallel analysis of RNA structure (PARS) and the selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE). Further, Riprap detects the experimentally validated riboSNitch that regulates human catechol-O-methyltransferase haplotypes and outputs structurally disrupted regions precisely at base resolution. Riprap provides a new approach to interpreting disease-related genetic variants. In addition, we construct a database (RiboSNitchDB) that includes the annotation and visualization of all presented riboSNitches in this study as well as 24 629 predicted riboSNitches from human expression quantitative trait loci.

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.

RNA Conformation Capture by Proximity Ligation

RNA proximity ligation is a set of molecular biology techniques used to analyze the conformations and spatial proximity of RNA molecules within cells. A typical experiment starts with cross-linking of a biological sample using UV light or psoralen, followed by partial fragmentation of RNA, RNA-RNA ligation, library preparation, and high-throughput sequencing. In the past decade, proximity ligation has been used to study structures of individual RNAs, networks of interactions between small RNAs and their targets, and whole RNA-RNA interactomes, in models ranging from bacteria to animal tissues and whole animals. Here, we provide an overview of the field, highlight the main findings, review the recent experimental and computational developments, and provide troubleshooting advice for new users. In the final section, we draw parallels between DNA and RNA proximity ligation and speculate on possible future research directions.

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.