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.

iM-Seeker: a webserver for DNA i-motifs prediction and scoring via automated machine learning.

DNA, beyond its canonical B-form double helix, adopts various alternative conformations, among which the i-motif, emerging in cytosine-rich sequences under acidic conditions, holds significant biological implications in transcription modulation and telomere biology. Despite recognizing the crucial role of i-motifs, predictive software for i-motif forming sequences has been limited. Addressing this gap, we introduce 'iM-Seeker', an innovative computational platform designed for the prediction and evaluation of i-motifs. iM-Seeker exhibits the capability to identify potential i-motifs within DNA segments or entire genomes, calculating stability scores for each predicted i-motif based on parameters such as the cytosine tracts number, loop lengths, and sequence composition. Furthermore, the webserver leverages automated machine learning (AutoML) to effortlessly fine-tune the optimal i-motif scoring model, incorporating user-supplied experimental data and customised features. As an advanced, versatile approach, 'iM-Seeker' promises to advance genomic research, highlighting the potential of i-motifs in cell biology and therapeutic applications. The webserver is freely available at https://im-seeker.org.

A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics.

Plants, as sessile organisms, deploy transcriptional dynamics for adapting to extreme growth conditions such as cold stress. Emerging evidence suggests that chromatin architecture contributes to transcriptional regulation. However, the relationship between chromatin architectural dynamics and transcriptional reprogramming in response to cold stress remains unclear. Here, we apply a chemical-crosslinking assisted proximity capture (CAP-C) method to elucidate the fine-scale chromatin landscape, revealing chromatin interactions within gene bodies closely associated with RNA polymerase II (Pol II) densities across initiation, pausing, and termination sites. We observe dynamic changes in chromatin interactions alongside Pol II activity alterations during cold stress, suggesting local chromatin dynamics may regulate Pol II activity. Notably, cold stress does not affect large-scale chromatin conformations. We further identify a comprehensive promoter-promoter interaction (PPI) network across the genome, potentially facilitating co-regulation of gene expression in response to cold stress. Our study deepens the understanding of chromatin conformation-associated gene regulation in plant response to cold.

Systematic identification of wheat spike developmental regulators by integrated multi-omics, transcriptional network, GWAS, and genetic analyses.

The spike architecture of wheat plays a crucial role in determining grain number, making it a key trait for optimization in wheat breeding programs. In this study, we used a multi-omic approach to analyze the transcriptome and epigenome profiles of the young spike at eight developmental stages, revealing coordinated changes in chromatin accessibility and H3K27me3 abundance during the flowering transition. We constructed a core transcriptional regulatory network (TRN) that drives wheat spike formation and experimentally validated a multi-layer regulatory module involving TaSPL15, TaAGLG1, and TaFUL2. By integrating the TRN with genome-wide association studies, we identified 227 transcription factors, including 42 with known functions and 185 with unknown functions. Further investigation of 61 novel transcription factors using multiple homozygous mutant lines revealed 36 transcription factors that regulate spike architecture or flowering time, such as TaMYC2-A1, TaMYB30-A1, and TaWRKY37-A1. Of particular interest, TaMYB30-A1, downstream of and repressed by WFZP, was found to regulate fertile spikelet number. Notably, the excellent haplotype of TaMYB30-A1, which contains a C allele at the WFZP binding site, was enriched during wheat breeding improvement in China, leading to improved agronomic traits. Finally, we constructed a free and open access Wheat Spike Multi-Omic Database (http://39.98.48.156:8800/#/). Our study identifies novel and high-confidence regulators and offers an effective strategy for dissecting the genetic basis of wheat spike development, with practical value for wheat breeding.

Prediction of DNA i-motifs via machine learning.

i-Motifs (iMs), are secondary structures formed in cytosine-rich DNA sequences and are involved in multiple functions in the genome. Although putative iM forming sequences are widely distributed in the human genome, the folding status and strength of putative iMs vary dramatically. Much previous research on iM has focused on assessing the iM folding properties using biophysical experiments. However, there are no dedicated computational tools for predicting the folding status and strength of iM structures. Here, we introduce a machine learning pipeline, iM-Seeker, to predict both folding status and structural stability of DNA iMs. The programme iM-Seeker incorporates a Balanced Random Forest classifier trained on genome-wide iMab antibody-based CUT&Tag sequencing data to predict the folding status and an Extreme Gradient Boosting regressor to estimate the folding strength according to both literature biophysical data and our in-house biophysical experiments. iM-Seeker predicts DNA iM folding status with a classification accuracy of 81% and estimates the folding strength with coefficient of determination (R2) of 0.642 on the test set. Model interpretation confirms that the nucleotide composition of the C-rich sequence significantly affects iM stability, with a positive correlation with sequences containing cytosine and thymine and a negative correlation with guanine and adenine.

Capture the in vivo intact RNA structurome by CAP-STRUCTURE-seq.

RNA decay serves as a crucial mechanism for maintaining cellular homeostasis and regulating gene expression. Large-scale analyses indicate that altered rates of decay contribute significantly to changes in mRNA levels, with up to half of these changes attributed to decay. The regulation of RNA decay is, at least in part, through structured RNA elements, especially in the non-coding regions of the mRNAs. The development of next-generation sequencing, and in vivo chemical probing techniques has allowed for unprecedented understanding of RNA folding in vivo and genome-wide. To explore the RNA structure elements that are responsible for RNA cleavage, we need to capture the RNA structure before cleavage. In this method, we introduce a new experimental procedure called CAP-STRUCTURE-seq, a modified STRUCTURE-Seq approach combining with the enrichment of in intact mRNAs by the use of terminator exonuclease treatment (5'-Phosphate-Dependent Exonuclease) that digests RNA containing 5-monophosphate ends. This approach is designed to investigate the RNA structure for these intact RNAs, providing a means to study the impact of RNA structure on RNA decay in greater detail. This method can provide insights into the function of RNA structure in RNA decay and help advance our understanding of biological processes.

3' Untranslated Region Structural Elements in CYP24A1 Are Associated With Infantile Hypercalcemia Type 1.

Loss-of-function mutations in the CYP24A1 protein-coding region causing reduced 25 hydroxyvitamin D (25OHD) and 1,25 dihydroxyvitamin D (1,25(OH)2 D) catabolism have been observed in some cases of infantile hypercalcemia type 1 (HCINF1), which can manifest as nephrocalcinosis, hypercalcemia and adult-onset hypercalciuria, and renal stone formation. Some cases present with apparent CYP24A1 phenotypes but do not exhibit pathogenic mutations. Here, we assessed the molecular mechanisms driving apparent HCINF1 where there was a lack of CYP24A1 mutation. We obtained blood samples from 47 patients with either a single abnormality of no obvious cause or a combination of hypercalcemia, hypercalciuria, and nephrolithiasis as part of our metabolic and stone clinics. We used liquid chromatography tandem mass spectrometry (LC-MS/MS) to determine serum vitamin D metabolites and direct sequencing to confirm CYP24A1 genotype. Six patients presented with profiles characteristic of altered CYP24A1 function but lacked protein-coding mutations in CYP24A1. Analysis upstream and downstream of the coding sequence showed single nucleotide variants (SNVs) in the CYP24A1 3' untranslated region (UTR). Bioinformatics approaches revealed that these 3' UTR abnormalities did not result in microRNA silencing but altered the CYP24A1 messenger RNA (mRNA) secondary structure, which negatively impacted translation. Our experiments showed that mRNA misfolding driven by these 3' UTR sequence-dependent structural elements was associated with normal 25OHD but abnormal 1,25(OH)2 D catabolism. Using CRISPR-Cas9 gene editing, we developed an in vitro mutant model for future CYP24A1 studies. Our results form a basis for future studies investigating structure-function relationships and novel CYP24A1 mutations producing a semifunctional protein. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Beyond transcription: compelling open questions in plant RNA biology.

The study of RNAs has become one of the most influential research fields in contemporary biology and biomedicine. In the last few years, new sequencing technologies have produced an explosion of new and exciting discoveries in the field but have also given rise to many open questions. Defining these questions, together with old, long-standing gaps in our knowledge, is the spirit of this article. The breadth of topics within RNA biology research is vast, and every aspect of the biology of these molecules contains countless exciting open questions. Here, we asked 12 groups to discuss their most compelling question among some plant RNA biology topics. The following vignettes cover RNA alternative splicing; RNA dynamics; RNA translation; RNA structures; R-loops; epitranscriptomics; long non-coding RNAs; small RNA production and their functions in crops; small RNAs during gametogenesis and in cross-kingdom RNA interference; and RNA-directed DNA methylation. In each section, we will present the current state-of-the-art in plant RNA biology research before asking the questions that will surely motivate future discoveries in the field. We hope this article will spark a debate about the future perspective on RNA biology and provoke novel reflections in the reader.

Reference genes for quantitative Arabidopsis single molecule RNA fluorescence in situ hybridization.

Subcellular mRNA quantities and spatial distributions are fundamental for driving gene regulatory programmes. Single molecule RNA fluorescence in situ hybridization (smFISH) uses fluorescent probes to label individual mRNA molecules, thereby facilitating both localization and quantitative studies. Validated reference mRNAs function as positive controls and are required for calibration. Here we present selection criteria for the first set of Arabidopsis smFISH reference genes. Following sequence and transcript data assessments, four mRNA probe sets were selected for imaging. Transcript counts per cell, correlations with cell size, and corrected fluorescence intensities were all calculated for comparison. In addition to validating reference probe sets, we present sample preparation steps that can retain green fluorescent protein fluorescence, thereby providing a method for simultaneous RNA and protein detection. In summary, our reference gene analyses, modified protocol, and simplified quantification method together provide a firm foundation for future quantitative single molecule RNA studies in Arabidopsis root apical meristem cells.

G4Atlas: a comprehensive transcriptome-wide G-quadruplex database.

RNA G-quadruplex (rG4) is a vital RNA tertiary structure motif that involves the base pairs on both Hoogsteen and Watson-Crick faces of guanines. rG4 is of great importance in the post-transcriptional regulation of gene expression. Experimental technologies have advanced to identify in vitro and in vivo rG4s across diverse transcriptomes. Building on these recent advances, here we present G4Atlas, the first transcriptome-wide G-quadruplex database, in which we have collated, classified, and visualized transcriptome rG4 experimental data, generated from rG4-seq, chemical profiling and ligand-binding methods. Our comprehensive database includes transcriptome-wide rG4s generated from 82 experimental treatments and 238 samples across ten species. In addition, we have also included RNA secondary structure prediction information across both experimentally identified and unidentified rG4s to enable users to display any potential competitive folding between rG4 and RNA secondary structures. As such, G4Atlas will enable users to explore the general functions of rG4s in diverse biological processes. In addition, G4Atlas lays the foundation for further data-driven deep learning algorithms to examine rG4 structural features.

RNA G-quadruplex structure contributes to cold adaptation in plants.

Nucleotide composition is suggested to infer gene functionality and ecological adaptation of species to distinct environments. However, the underlying biological function of nucleotide composition dictating environmental adaptations is largely unknown. Here, we systematically analyze the nucleotide composition of transcriptomes across 1000 plants (1KP) and their corresponding habitats. Intriguingly, we find that plants growing in cold climates have guanine (G)-enriched transcriptomes, which are prone to forming RNA G-quadruplex structures. Both immunofluorescence detection and in vivo structure profiling reveal that RNA G-quadruplex formation in plants is globally enhanced in response to cold. Cold-responsive RNA G-quadruplexes strongly enhanced mRNA stability, rather than affecting translation. Disruption of individual RNA G-quadruplex promotes mRNA decay in the cold, leading to impaired plant cold response. Therefore, we propose that plants adopted RNA G-quadruplex structure as a molecular signature to facilitate their adaptation to the cold during evolution.

In vivo single-molecule analysis reveals COOLAIR RNA structural diversity.

Cellular RNAs are heterogeneous with respect to their alternative processing and secondary structures, but the functional importance of this complexity is still poorly understood. A set of alternatively processed antisense non-coding transcripts, which are collectively called COOLAIR, are generated at the Arabidopsis floral-repressor locus FLOWERING LOCUS C (FLC)1. Different isoforms of COOLAIR influence FLC transcriptional output in warm and cold conditions2-7. Here, to further investigate the function of COOLAIR, we developed an RNA structure-profiling method to determine the in vivo structure of single RNA molecules rather than the RNA population average. This revealed that individual isoforms of the COOLAIR transcript adopt multiple structures with different conformational dynamics. The major distally polyadenylated COOLAIR isoform in warm conditions adopts three predominant structural conformations, the proportions and conformations of which change after cold exposure. An alternatively spliced, strongly cold-upregulated distal COOLAIR isoform6 shows high structural diversity, in contrast to proximally polyadenylated COOLAIR. A hyper-variable COOLAIR structural element was identified that was complementary to the FLC transcription start site. Mutations altering the structure of this region changed FLC expression and flowering time, consistent with an important regulatory role of the COOLAIR structure in FLC transcription. Our work demonstrates that isoforms of non-coding RNA transcripts adopt multiple distinct and functionally relevant structural conformations, which change in abundance and shape in response to external conditions.

Deep Learning in RNA Structure Studies.

Deep learning, or artificial neural networks, is a type of machine learning algorithm that can decipher underlying relationships from large volumes of data and has been successfully applied to solve structural biology questions, such as RNA structure. RNA can fold into complex RNA structures by forming hydrogen bonds, thereby playing an essential role in biological processes. While experimental effort has enabled resolving RNA structure at the genome-wide scale, deep learning has been more recently introduced for studying RNA structure and its functionality. Here, we discuss successful applications of deep learning to solve RNA problems, including predictions of RNA structures, non-canonical G-quadruplex, RNA-protein interactions and RNA switches. Following these cases, we give a general guide to deep learning for solving RNA structure problems.

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.

A cytosolic pentatricopeptide repeat protein is essential for tapetal plastid development by regulating OsGLK1 transcript levels in rice.

Most plant pentatricopeptide repeat (PPR) proteins localize to and function inside plastids and mitochondria. However, the function of PPRs that only localize to the cytoplasm remains unknown. Here, we demonstrated that the rice (Oryza sativa) PPR protein CYTOPLASM-LOCALIZED PPR1 (OsCPPR1) contributes to pollen development and localizes to the cytoplasm. Knocking down OsCPPR1 led to abnormal plastid development in tapetal cells, prolonged tapetal programmed cell death (PCD) and tapetum degradation, and significantly reduced pollen fertility. Transcriptome analysis revealed that the transcript level of OsGOLDEN-LIKE1 (OsGLK1), which encodes a transcription factor that regulates plastid development and maintenance, was significantly higher in the OsCPPR1 knockdown plants compared to wild-type plants. We further determined that OsCPPR1 downregulates OsGLK1 transcription by directly binding to the single-stranded regions of OsGLK1 mRNAs. Overexpression of OsGLK1 resulted in abnormal tapetum and plastid development, similar to that seen in OsCPPR1 knockdown plants, and suppression of OsGLK1 partially restored pollen fertility in the OsCPPR1 knockdown plants. We therefore conclude that OsCPPR1 suppresses OsGLK1 in the regulation of plastid development and PCD in the tapetum. Our work revealed novel functions for a cytosolic PPR, demonstrating the diverse roles of PPRs in plants and identifying a new regulatory mechanism for regulating pollen development in rice.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry.

BACKGROUND: Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry. RESULTS: Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat. CONCLUSIONS: We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

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.

Targeting the Conserved Stem Loop 2 Motif in the SARS-CoV-2 Genome.

RNA structural elements occur in numerous single-stranded positive-sense RNA viruses. The stem-loop 2 motif (s2m) is one such element with an unusually high degree of sequence conservation, being found in the 3' untranslated region (UTR) in the genomes of many astroviruses, some picornaviruses and noroviruses, and a variety of coronaviruses, including severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2. The evolutionary conservation and its occurrence in all viral subgenomic transcripts imply a key role for s2m in the viral infection cycle. Our findings indicate that the element, while stably folded, can nonetheless be invaded and remodeled spontaneously by antisense oligonucleotides (ASOs) that initiate pairing in exposed loops and trigger efficient sequence-specific RNA cleavage in reporter assays. ASOs also act to inhibit replication in an astrovirus replicon model system in a sequence-specific, dose-dependent manner and inhibit SARS-CoV-2 replication in cell culture. Our results thus permit us to suggest that the s2m element is readily targeted by ASOs, which show promise as antiviral agents. IMPORTANCE The highly conserved stem-loop 2 motif (s2m) is found in the genomes of many RNA viruses, including SARS-CoV-2. Our findings indicate that the s2m element can be targeted by antisense oligonucleotides. The antiviral potential of this element represents a promising start for further research into targeting conserved elements in RNA viruses.

In vivo nuclear RNA structurome reveals RNA-structure regulation of mRNA processing in plants.

BACKGROUND: mRNA processing is critical for gene expression. A challenge in regulating mRNA processing is how to recognize the actual mRNA processing sites, such as splice and polyadenylation sites, when the sequence content is insufficient for this purpose. Previous studies suggested that RNA structure affects mRNA processing. However, the regulatory role of RNA structure in mRNA processing remains unclear. RESULTS: Here, we perform in vivo selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemical profiling on Arabidopsis and generate the in vivo nuclear RNA structure landscape. We find that nuclear mRNAs fold differently from cytosolic mRNAs across translation start and stop sites. Notably, we discover a two-nucleotide single-stranded RNA structure feature upstream of 5' splice sites that is strongly associated with splicing and the selection of alternative 5' splice sites. The regulatory role of this RNA structure feature is further confirmed by experimental validation. Moreover, we find the single-strandedness of branch sites is also associated with 3' splice site recognition. We also identify an RNA structure feature comprising two close-by single-stranded regions that is specifically associated with both polyadenylation and alternative polyadenylation events. CONCLUSIONS: We successfully identify pre-mRNA structure features associated with splicing and polyadenylation at whole-genome scale and validate an RNA structure feature which can regulate splicing. Our study unveils a new RNA structure regulatory mechanism for mRNA processing.