In vivo single-molecule analysis reveals COOLAIR RNA structural diversity

RNA Structure: Function and Application in Plant Biology

RNA orchestrates intricate structures that influence gene expression and protein production in all living organisms, with implications for fundamental biology, medicine, and agriculture. Although extensive research has been conducted on RNA biology, many regulatory mechanisms remain elusive due to the complex and dynamic nature of RNA structures and past technological limitations. Recent advancements in RNA structure technology have revolutionized plant RNA biology research. Here, we review cutting-edge technologies for studying RNA structures in plants and their functional significance in diverse biological processes. Additionally, we highlight the pivotal role of RNA structure in influencing plant growth, development, and responses to environmental stresses. We also discuss the potential evolutionary significance of RNA structure in natural adaptation and crop domestication. Finally, we propose leveraging RNA structure-mediated gene regulation as an innovative strategy to bolster plant resilience against climate change.

Rapid folding of nascent RNA regulates eukaryotic RNA biogenesis

RNA's catalytic, regulatory, or coding potential depends on structure formation. Because base pairing occurs during transcription, early structural states can govern RNA processing events and dictate the formation of functional conformations. These co-transcriptional states remain mostly unknown. Here, we develop co-transcriptional structure tracking (CoSTseq), which detects nascent RNA base pairing within and upon exit from RNA polymerases (Pols) transcriptome wide in living yeast cells. Monitoring each nucleotide's base pairing activity during transcription, CoSTseq reveals predominantly rapid pairing-within 25 bp of transcription after addition to the nascent chain. Moreover, ∼23% of rRNA nucleotides attain their final base pairing state near Pol I, while most other nucleotides must undergo changes in pairing status during later steps of ribosome biogenesis. We show that helicases act immediately to remodel structures across the rDNA locus to facilitate ribosome biogenesis. By contrast, nascent pre-mRNAs attain local structures indistinguishable from mature mRNAs, suggesting that refolding behind elongating ribosomes resembles co-transcriptional folding behind Pol II.

Telomerase RNA evolution: a journey from plant telomeres to broader eukaryotic diversity

Telomeres, essential for maintaining genomic stability, are typically preserved through the action of telomerase, a ribonucleoprotein complex that synthesizes telomeric DNA. One of its two core components, telomerase RNA (TR), serves as the template for this synthesis, and its evolution across different species is both complex and diverse. This review discusses recent advancements in understanding TR evolution, with a focus on plants (Viridiplantae). Utilizing novel bioinformatic tools and accumulating genomic and transcriptomic data, combined with corresponding experimental validation, researchers have begun to unravel the intricate pathways of TR evolution and telomere maintenance mechanisms. Contrary to previous beliefs, a monophyletic origin of TR has been demonstrated first in land plants and subsequently across the broader phylogenetic megagroup Diaphoretickes. Conversely, the discovery of plant-type TRs in insects challenges assumptions about the monophyletic origin of TRs in animals, suggesting evolutionary innovations coinciding with arthropod divergence. The review also highlights key challenges in TR identification and provides examples of how these have been addressed. Overall, this work underscores the importance of expanding beyond model organisms to comprehend the full complexity of telomerase evolution, with potential applications in agriculture and biotechnology.

The conformational space of RNase P RNA in solution

RNA conformational diversity has fundamental biological roles1-5, but direct visualization of its full conformational space in solution has not been possible using traditional biophysical techniques. Using solution atomic force microscopy, a deep neural network and statistical analyses, we show that the ribonuclease P (RNase P) RNA adopts heterogeneous conformations consisting of a conformationally invariant core and highly flexible peripheral structural elements that sample a broad conformational space, with amplitudes as large as 20-60 Å in a multitude of directions, with very low net energy cost. Increasing Mg2+ drives compaction and enhances enzymatic activity, probably by narrowing the conformational space. Moreover, analyses of the correlations and anticorrelations between spatial flexibility and sequence conservation suggest that the functional roles of both the structure and dynamics of key regions are embedded in the primary sequence. These findings reveal the structure-dynamics basis for the embodiment of both enzymatic precision and substrate promiscuity in the RNA component of the RNase P. Mapping the conformational space of the RNase P RNA demonstrates a new general approach to studying RNA structure and dynamics.

Carbon and nitrogen signaling regulate FLOWERING LOCUS C and impact flowering time in Arabidopsis

The timing of flowering in plants is modulated by both carbon (C) and nitrogen (N) signaling pathways. In a previous study, we established a pivotal role of the sucrose-signaling trehalose 6-phosphate pathway in regulating flowering under N-limited short-day conditions. In this work, we show that both wild-type Arabidopsis (Arabidopsis thaliana) plants grown under N-limited conditions and knock-down plants of TREHALOSE PHOSPHATE SYNTHASE 1 induce FLOWERING LOCUS C (FLC) expression, a well-known floral repressor associated with vernalization. When exposed to an extended period of cold, a flc mutant fails to respond to N availability and flowers at the same time under N-limited and full-nutrition conditions. Our data suggest that SUCROSE NON-FERMENTING 1 RELATED KINASE 1-dependent trehalose 6-phosphate-mediated C signaling and a mechanism downstream of N signaling (likely involving NIN-LIKE PROTEIN 7) impact the expression of FLC. Collectively, our data underscore the existence of a multi-factor regulatory system in which the C and N signaling pathways jointly govern the regulation of flowering in plants.

Protocol for mapping RNA G-quadruplex for chromatin-bound RNA using d-rG4-seq

Here, we present a protocol for using d-rG4-seq, a technique for mapping RNA G-quadruplex (rG4) for chromatin-bound RNA. We describe steps for identifying in vivo rG4 structures based on differential sensitivity of rG4 to dimethyl sulfate (DMS) modification, folding in the presence of monovalent cations, K+ versus Li+, and reverse transcriptase (RT) readthrough when folded. We then detail procedures for isolating RNA from the chromatin-bound fractions to enrich for epigenetic regulators and comparing in vitro versus in vivo profiles. For complete details on the use and execution of this protocol, please refer to Lee et al.1.

An interpretable RNA foundation model for exploring functional RNA motifs in plants

The complex 'language' of plant RNA encodes a vast array of biological regulatory elements that orchestrate crucial aspects of plant growth, development and adaptation to environmental stresses. Recent advancements in foundation models (FMs) have demonstrated their unprecedented potential to decipher complex 'language' in biology. In this study, we introduced PlantRNA-FM, a high-performance and interpretable RNA FM specifically designed for plants. PlantRNA-FM was pretrained on an extensive dataset, integrating RNA sequences and RNA structure information from 1,124 distinct plant species. PlantRNA-FM exhibits superior performance in plant-specific downstream tasks. PlantRNA-FM achieves an F1 score of 0.974 for genic region annotation, whereas the current best-performing model achieves 0.639. Our PlantRNA-FM is empowered by our interpretable framework that facilitates the identification of biologically functional RNA sequence and structure motifs, including both RNA secondary and tertiary structure motifs across transcriptomes. Through experimental validations, we revealed translation-associated RNA motifs in plants. Our PlantRNA-FM also highlighted the importance of the position information of these functional RNA motifs in genic regions. Taken together, our PlantRNA-FM facilitates the exploration of functional RNA motifs across the complexity of transcriptomes, empowering plant scientists with capabilities for programming RNA codes in plants.

RNA tertiary structure and conformational dynamics revealed by BASH MaP

The functional effects of an RNA can arise from complex three-dimensional folds known as tertiary structures. However, predicting the tertiary structure of an RNA and whether an RNA adopts distinct tertiary conformations remains challenging. To address this, we developed BASH MaP, a single-molecule dimethyl sulfate (DMS) footprinting method and DAGGER, a computational pipeline, to identify alternative tertiary structures adopted by different molecules of RNA. BASH MaP utilizes potassium borohydride to reveal the chemical accessibility of the N7 position of guanosine, a key mediator of tertiary structures. We used BASH MaP to identify diverse conformational states and dynamics of RNA G-quadruplexes, an important RNA tertiary motif, in vitro and in cells. BASH MaP and DAGGER analysis of the fluorogenic aptamer Spinach reveals that it adopts alternative tertiary conformations which determine its fluorescence states. BASH MaP thus provides an approach for structural analysis of RNA by revealing previously undetectable tertiary structures.

Nexus between RNA conformational dynamics and functional versatility

RNA conformational dynamics is pivotal for functional regulations in biology. RNA can function as versatile as protein but adopts multiple distinct structures. In this review, we provide a focused review of the recent advances in studies of RNA conformational dynamics and address some of the misconceptions about RNA structure and its conformational dynamics. We discuss why the traditional methods for structure determination come up short in describing RNA conformational space. The examples discussed provide illustrations of the structure-based mechanisms of RNAs with diverse roles, including viral, long noncoding, and catalytic RNAs, one of which focuses on the debated area of conformational heterogeneity of an RNA structural element in the HIV-1 genome.

The role of RNA structure in 3' end processing in eukaryotes

Maturation of pre-mRNA into fully functional mRNA involves a series of highly coordinated steps that are essential for eukaryotic gene expression. RNA structure has been found to play regulatory roles in many of these steps, including cleavage, polyadenylation, and termination. Recent advances in structure probing techniques have been instrumental in revealing how nascent transcript conformation contributes to these dynamic, co-transcriptional processes. In this review, we present examples where RNA structure affects accessibility and/or function of key processing enzymes, thereby influencing the efficiency and precision of 3' end processing machinery. We also discuss emerging technologies that could further enhance our understanding of RNA structure mediated regulation of 3' end processing.

RNA ensembles from in vitro to in vivo: Toward predictive models of RNA cellular function

Deepening our understanding of RNA biology and accelerating development of RNA-based therapeutics go hand-in-hand-both requiring a transition from qualitative descriptions of RNA structure to quantitative models capable of predicting RNA behaviors, and from a static to an ensemble view. Ensembles are determined from their free energy landscapes, which define the relative populations of conformational states and the energetic barriers separating them. Experimental determination of RNA ensembles over the past decade has led to powerful predictive models of RNA behavior in vitro. It has also been shown during this time that the cellular environment redistributes RNA ensembles, changing the abundances of functionally relevant conformers relative to in vitro contexts with subsequent functional RNA consequences. However, recent studies have demonstrated that testing models built from in vitro ensembles with highly quantitative measurements of RNA cellular function, aided by emerging computational methodologies, enables predictive modelling of cellular activity and biological discovery.

SVALKA: A Long Noncoding Cis-Natural Antisense RNA That Plays a Role in the Regulation of the Cold Response of Arabidopsis thaliana

RNA plays important roles in the regulation of gene expression in response to environmental stimuli. SVALKA, a long noncoding cis-natural antisense RNA, is a key component of regulating the response to cold temperature in Arabidopsis thaliana. There are three mechanisms through which SVALKA fine tunes the transcriptional response to cold temperatures. SVALKA regulates the expression of the CBF1 (C-Repeat Dehydration Binding Factor 1) transcription factor through a collisional transcription mechanism and a dsRNA and DICER mediated mechanism. SVALKA also interacts with Polycomb Repressor Complex 2 to regulate the histone methylation of CBF3. Both CBF1 and CBF3 are key components of the COLD REGULATED (COR) regulon that direct the plant's response to cold temperature over time, as well as plant drought adaptation, pathogen responses, and growth regulation. The different isoforms of SVALKA and its potential to form dynamic RNA conformations are important features in regulating a complex gene network in concert with several other noncoding RNA. This review will summarize the three mechanisms through which SVALKA participates in gene regulation, describe the ways that dynamic RNA structures support the function of regulatory noncoding RNA, and explore the potential for improving agricultural genetic engineering with a better understanding of the roles of noncoding RNA.

Rapid folding of nascent RNA regulates eukaryotic RNA biogenesis

An RNA's catalytic, regulatory, or coding potential depends on RNA structure formation. Because base pairing occurs during transcription, early structural states can govern RNA processing events and dictate the formation of functional conformations. These co-transcriptional states remain unknown. Here, we develop CoSTseq, which detects nascent RNA base pairing within and upon exit from RNA polymerases (Pols) transcriptome-wide in living yeast cells. By monitoring each nucleotide's base pairing activity during transcription, we identify distinct classes of behaviors. While 47% of rRNA nucleotides remain unpaired, rapid and delayed base pairing - with rates of 48.5 and 13.2 kb -1 of transcribed rDNA, respectively - typically completes when Pol I is only 25 bp downstream. We show that helicases act immediately to remodel structures across the rDNA locus and facilitate ribosome biogenesis. In contrast, nascent pre-mRNAs attain local structures indistinguishable from mature mRNAs, suggesting that refolding behind elongating ribosomes resembles co-transcriptional folding behind Pol II.

Nanopore direct RNA sequencing reveals N6-methyladenosine and polyadenylation landscapes on long non-coding RNAs in Arabidopsis thaliana

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in various biological processes, including stage development in plants. N6-methyladenosine (m6A) modification and polyadenylation are noteworthy regulatory processes that impact transcript functions by modulating their abundance. However, the specific landscapes of m6A modification and polyadenylation on lncRNAs remain largely unexplored. The advent of nanopore direct RNA sequencing (DRS) provides unprecedented opportunities for directly detecting m6A modifications and estimating polyadenine (poly[A]) tail lengths on individual RNA molecules.

RESULTS

Here we utilized nanopore DRS to identify lncRNAs and map the transcriptome-wide m6A modification and polyadenylation landscapes in the model plant Arabidopsis thaliana. Leveraging the Low-abundance Aware Full-length Isoform clusTEr (LAFITE) assembly pipeline, we identified 1149 novel lncRNAs in seventeen nanopore DRS datasets from the wild-type Columbia-0. Through the precise detection of 2381 m6A modification sites on lncRNAs at single-base resolution, we observed that lncRNAs exhibited lower methylation levels compared to protein-coding RNAs, and m6A modification facilitated lncRNA abundance. Additionally, we estimated the poly(A) tail lengths of individual lncRNAs and found that poly(A) tails contributed to lncRNA stability, while their effect was not length-dependent. Furthermore, by comparing lncRNA abundance between 2-week seedlings and 5-week floral buds, we revealed the dynamic expression patterns of lncRNAs during the transition from the vegetative stage to the reproductive stage. These observations provided insights into their potential roles in specific tissues or stages in Arabidopsis, including regulating stage development. Moreover, by integrating information on m6A modification, we unveiled a positive correlation between methylation variances and differential expressions of lncRNAs during stage development.

CONCLUSIONS

These findings highlighted the significance of epigenetic modification and post-transcriptional processing in shaping lncRNA expression and their functions during Arabidopsis stage development, contributing to the growing field of lncRNA research in plants.

CLINICAL TRIAL NUMBER

Not applicable.

Cotranscriptional splicing is required in the cold to produce COOLAIR isoforms that repress Arabidopsis FLC

Plants use seasonal cold to time the transition to reproductive development. Short- and long-term cold exposure is registered via parallel transcriptional shutdown and Polycomb-dependent epigenetic silencing of the Arabidopsis thaliana major flowering repressor locus FLOWERING LOCUS C (FLC). The cold-induced antisense transcripts (COOLAIR) determine the dynamics of FLC transcriptional shutdown, but the thermosensory mechanisms are still unresolved. Here, through a forward genetic screen, we identify a mutation that perturbs cold-induced COOLAIR expression and FLC repression. The mutation is a hypomorphic allele of SUPPRESSORS OF MEC-8 AND UNC-52 1 (SMU1), a conserved subunit of the spliceosomal B complex. SMU1 interacts in vivo with the proximal region of nascent COOLAIR and RNA 3' processing/cotranscriptional regulators and enhances COOLAIR proximal intron splicing to promote specific COOLAIR isoforms. SMU1 also interacts with ELF7, an RNA Polymerase II Associated Factor (Paf1) component and limits COOLAIR transcription. Cold thus changes cotranscriptional splicing/RNA Pol II functionality in an SMU1-dependent mechanism to promote two different isoforms of COOLAIR that lead to reduced FLC transcription. Such cotranscriptional mechanisms are emerging as important regulators underlying plasticity in gene expression.

Regulatory Networks Underlying Plant Responses and Adaptation to Cold Stress

Cold is an important environmental factor limiting plant growth and development. Recent studies have revealed the complex regulatory networks associated with plant responses to cold and identified their interconnections with signaling pathways related to light, the circadian clock, plant hormones, and pathogen defense. In this article, we review recent advances in understanding the molecular basis of cold perception and signal transduction pathways. We also summarize recent developments in the study of cold-responsive growth and flowering. Finally, we propose future directions for the study of long-term cold sensing, RNA secondary structures in response to cold, and the development of cold-tolerant and high-yield crops.

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.

Structural and biophysical dissection of RNA conformational ensembles

RNA's ability to form and interconvert between multiple secondary and tertiary structures is critical to its functional versatility and the traditional view of RNA structures as static entities has shifted towards understanding them as dynamic conformational ensembles. In this review we discuss RNA structural ensembles and their dynamics, highlighting the concept of conformational energy landscapes as a unifying framework for understanding RNA processes such as folding, misfolding, conformational changes, and complex formation. Ongoing advancements in cryo-electron microscopy and chemical probing techniques are significantly enhancing our ability to investigate multiple structures adopted by conformationally dynamic RNAs, while traditional methods such as nuclear magnetic resonance spectroscopy continue to play a crucial role in providing high-resolution, quantitative spatial and temporal information. We discuss how these methods, when used synergistically, can provide a comprehensive understanding of RNA conformational ensembles, offering new insights into their regulatory functions.

Visualizing RNA structure ensembles by single-molecule correlated chemical probing

RNA molecules fold to form complex internal structures. Many of these RNA structures populate ensembles with rheostat-like properties, with each state having a distinct function. Until recently, analysis of RNA structures, especially within cells, was limited to modeling either a single averaged structure or computationally-modeled ensembles. These approaches obscure the intrinsic heterogeneity of many structured RNAs. Single-molecule correlated chemical probing (smCCP) strategies are now making it possible to measure and deconvolute RNA structure ensembles based on efficiently executed chemical probing experiments. Here, we provide an overview of fundamental single-molecule probing principles, review current ensemble deconvolution strategies, and discuss recent applications to diverse biological systems. smCCP is enabling a revolution in understanding how the plasticity of RNA structure is exploited in biological systems to respond to stimuli and alter gene function. The energetics of RNA ensembles are often subtle and a subset can likely be targeted to modulate disease-associated biological processes.

SEISMICgraph: a web-based tool for RNA structure data visualization

In recent years, RNA has been increasingly recognized for its essential roles in biology, functioning not only as a carrier of genetic information but also as a dynamic regulator of gene expression through its interactions with other RNAs, proteins, and itself. Advances in chemical probing techniques have significantly enhanced our ability to identify RNA secondary structures and understand their regulatory roles. These developments, alongside improvements in experimental design and data processing, have greatly increased the resolution and throughput of structural analyses. Here, we introduce SEISMICgraph, a web-based tool designed to support RNA structure research by offering data visualization and analysis capabilities for a variety of chemical probing modalities. SEISMICgraph enables simultaneous comparison of data across different sequences and experimental conditions through a user-friendly interface that requires no programming expertise. We demonstrate its utility by investigating known and putative riboswitches and exploring how RNA modifications influence their structure and binding. SEISMICgraph's ability to rapidly visualize adenine-dependent structural changes and assess the impact of pseudouridylation on these transitions provides novel insights and establishes a roadmap for numerous future applications.

RNA tertiary structure and conformational dynamics revealed by BASH MaP

The functional effects of an RNA can arise from complex three-dimensional folds known as tertiary structures. However, predicting the tertiary structure of an RNA and whether an RNA adopts distinct tertiary conformations remains challenging. To address this, we developed BASH MaP, a single-molecule dimethyl sulfate (DMS) footprinting method and DAGGER, a computational pipeline, to identify alternative tertiary structures adopted by different molecules of RNA. BASH MaP utilizes potassium borohydride to reveal the chemical accessibility of the N7 position of guanosine, a key mediator of tertiary structures. We used BASH MaP to identify diverse conformational states and dynamics of RNA G-quadruplexes, an important RNA tertiary motif, in vitro and in cells. BASH MaP and DAGGER analysis of the fluorogenic aptamer Spinach reveals that it adopts alternative tertiary conformations which determine its fluorescence states. BASH MaP thus provides an approach for structural analysis of RNA by revealing previously undetectable tertiary structures.

Involvement of ribosomal protein L17 and Y-box binding protein 1 in the assembly of hepatitis C virus potentially via their interaction with the 3' untranslated region of the viral genome

The 3' untranslated region (3'UTR) of the hepatitis C virus (HCV) RNA genome, which contains a highly conserved 3' region named the 3'X-tail, plays an essential role in RNA replication and promotes viral IRES-dependent translation. Although our previous work has found a cis-acting element for genome encapsidation within 3'X, there is limited information on the involvement of the 3'UTR in particle formation. In this study, proteomic analyses identified host cell proteins that bind to the 3'UTR containing the 3'X region but not to the sequence lacking the 3'X. Further characterization showed that RNA-binding proteins, ribosomal protein L17 (RPL17), and Y-box binding protein 1 (YBX1) facilitate the efficient production of infectious HCV particles in the virus infection cells. Using small interfering RNA (siRNA)-mediated gene silencing in four assays that distinguish between the various stages of the HCV life cycle, RPL17 and YBX1 were found to be most important for particle assembly in the trans-packaging assay with replication-defective subgenomic RNA. In vitro assays showed that RPL17 and YBX1 bind to the 3'UTR RNA and deletion of the 3'X region attenuates their interaction. Knockdown of RPL17 or YBX1 resulted in reducing the amount of HCV RNA co-precipitating with the viral Core protein by RNA immunoprecipitation and increasing the relative distance in space between Core and double-stranded RNA by confocal imaging, suggesting that RPL17 and YBX1 potentially affect HCV RNA-Core interaction, leading to efficient nucleocapsid assembly. These host factors provide new clues to understanding the molecular mechanisms that regulate HCV particle formation.

IMPORTANCE

Although basic research on the HCV life cycle has progressed significantly over the past two decades, our understanding of the molecular mechanisms that regulate the process of particle formation, in particular encapsidation of the genome or nucleocapsid assembly, has been limited. We present here, for the first time, that two RNA-binding proteins, RPL17 and YBX1, bind to the 3'X in the 3'UTR of the HCV genome, which potentially acts as a packaging signal, and facilitates the viral particle assembly. Our study revealed that RPL17 and YBX1 exert a positive effect on the interaction between HCV RNA and Core protein, suggesting that the presence of both host factors modulate an RNA structure or conformation suitable for packaging the viral genome. These findings help us to elucidate not only the regulatory mechanism of the particle assembly of HCV but also the function of host RNA-binding proteins during viral infection.

LncRNAs: the art of being influential without protein

The plant long noncoding (lnc)RNA field is on the brink of transitioning from large-scale identification of lncRNAs to their functional characterization. Due to the cross-kingdom conservation of interaction types and molecular functions, there is much to be learned from mammalian lncRNA research. Here, we discuss the different molecular processes involving lncRNAs from the regulation of chromatin to splicing. Furthermore, we discuss the lncRNA interactome, which includes proteins, other RNAs, and DNA. We explore and discuss how mammalian lncRNA functionalities could be reflected in similar pathways in plants and hypothesize that several breakthroughs in mammalian research could lead to the discovery of novel plant lncRNA molecular functions. Expanding our knowledge of the biological role of lncRNAs and their multiple applications paves the way for future agricultural applications.

Probing RNA structure and dynamics using nanopore and next generation sequencing

It has become increasingly evident that the structures RNAs adopt are conformationally dynamic; the various structured states that RNAs sample govern their interactions with other nucleic acids, proteins, and ligands to regulate a myriad of biological processes. Although several biophysical approaches have been developed and used to study the dynamic landscape of structured RNAs, technical limitations have limited their application to all classes of RNA due to variable size and flexibility. Recent advances combining chemical probing experiments with next-generation- and direct sequencing have emerged as an alternative approach to exploring the conformational dynamics of RNA. In this review, we provide a methodological overview of the sequencing-based techniques used to study RNA conformational dynamics. We discuss how different techniques have enabled us to better understand the propensity of RNAs from a variety of different classes to sample multiple conformational states. Finally, we present examples of the ways these techniques have reshaped how we think about RNA structure.

Isoform-specific RNA structure determination using Nano-DMS-MaP

RNA structure determination is essential to understand how RNA carries out its diverse biological functions. In cells, RNA isoforms are readily expressed with partial variations within their sequences due, for example, to alternative splicing, heterogeneity in the transcription start site, RNA processing or differential termination/polyadenylation. Nanopore dimethyl sulfate mutational profiling (Nano-DMS-MaP) is a method for in situ isoform-specific RNA structure determination. Unlike similar methods that rely on short sequencing reads, Nano-DMS-MaP employs nanopore sequencing to resolve the structures of long and highly similar RNA molecules to reveal their previously hidden structural differences. This Protocol describes the development and applications of Nano-DMS-MaP and outlines the main considerations for designing and implementing a successful experiment: from bench to data analysis. In cell probing experiments can be carried out by an experienced molecular biologist in 3-4 d. Data analysis requires good knowledge of command line tools and Python scripts and requires a further 3-5 d.

Transcription regulation by long non-coding RNAs: mechanisms and disease relevance

Long non-coding RNAs (lncRNAs) outnumber protein-coding transcripts, but their functions remain largely unknown. In this Review, we discuss the emerging roles of lncRNAs in the control of gene transcription. Some of the best characterized lncRNAs have essential transcription cis-regulatory functions that cannot be easily accomplished by DNA-interacting transcription factors, such as XIST, which controls X-chromosome inactivation, or imprinted lncRNAs that direct allele-specific repression. A growing number of lncRNA transcription units, including CHASERR, PVT1 and HASTER (also known as HNF1A-AS1) act as transcription-stabilizing elements that fine-tune the activity of dosage-sensitive genes that encode transcription factors. Genetic experiments have shown that defects in such transcription stabilizers often cause severe phenotypes. Other lncRNAs, such as lincRNA-p21 (also known as Trp53cor1) and Maenli (Gm29348) contribute to local activation of gene transcription, whereas distinct lncRNAs influence gene transcription in trans. We discuss findings of lncRNAs that elicit a function through either activation of their transcription, transcript elongation and processing or the lncRNA molecule itself. We also discuss emerging evidence of lncRNA involvement in human diseases, and their potential as therapeutic targets.

Decoding plant specialized metabolism: new mechanistic insights

Secondary metabolite (SM) production provides biotic and abiotic stress resistance and enables plants to adapt to the environment. Biosynthesis of these metabolites involves a complex interplay between transcription factors (TFs) and regulatory elements, with emerging evidence suggesting an integral role for chromatin dynamics. Here we review key TFs and epigenetic regulators that govern SM biosynthesis in different contexts. We summarize relevant emerging technologies and results from the model species arabidopsis (Arabidopsis thaliana) and outline aspects of regulation that may also function in food, feed, fiber, oil, or industrial crop plants. Finally, we highlight how effective translation of fundamental knowledge from model to non-model species can benefit understanding of SM production in a variety of ecological, agricultural, and pharmaceutical contexts.

Flowering time: From physiology, through genetics to mechanism

Plant species have evolved different requirements for environmental/endogenous cues to induce flowering. Originally, these varying requirements were thought to reflect the action of different molecular mechanisms. Thinking changed when genetic and molecular analysis in Arabidopsis thaliana revealed that a network of environmental and endogenous signaling input pathways converge to regulate a common set of "floral pathway integrators." Variation in the predominance of the different input pathways within a network can generate the diversity of requirements observed in different species. Many genes identified by flowering time mutants were found to encode general developmental and gene regulators, with their targets having a specific flowering function. Studies of natural variation in flowering were more successful at identifying genes acting as nodes in the network central to adaptation and domestication. Attention has now turned to mechanistic dissection of flowering time gene function and how that has changed during adaptation. This will inform breeding strategies for climate-proof crops and help define which genes act as critical flowering nodes in many other species.

From environmental responses to adaptation: the roles of plant lncRNAs

As sessile organisms, plants are continuously exposed to heterogeneous and changing environments and constantly need to adapt their growth strategies. They have evolved complex mechanisms to recognize various stress factors, activate appropriate signaling pathways, and respond accordingly by reprogramming the expression of multiple genes at the transcriptional, post-transcriptional, and even epigenome levels to tolerate stressful conditions such as drought, high temperature, nutrient deficiency, and pathogenic interactions. Apart from protein-coding genes, long non-coding RNAs (lncRNAs) have emerged as key players in plant adaptation to environmental stresses. They are transcripts larger than 200 nucleotides without protein-coding potential. Still, they appear to regulate a wide range of processes, including epigenetic modifications and chromatin reorganization, as well as transcriptional and post-transcriptional modulation of gene expression, allowing plant adaptation to various environmental stresses. LncRNAs can positively or negatively modulate stress responses, affecting processes such as hormone signaling, temperature tolerance, and nutrient deficiency adaptation. Moreover, they also seem to play a role in stress memory, wherein prior exposure to mild stress enhances plant ability to adapt to subsequent stressful conditions. In this review, we summarize the contribution of lncRNAs in plant adaptation to biotic and abiotic stresses, as well as stress memory. The complex evolutionary conservation of lncRNAs is also discussed and provides insights into future research directions in this field.

Frontiers in plant RNA research in ICAR2023: from lab to innovative agriculture

The recent growth in global warming, soil contamination, and climate instability have widely disturbed ecosystems, and will have a significant negative impact on the growth of plants that produce grains, fruits and woody biomass. To conquer this difficult situation, we need to understand the molecular bias of plant environmental responses and promote development of new technologies for sustainable maintenance of crop production. Accumulated molecular biological data have highlighted the importance of RNA-based mechanisms for plant stress responses. Here, we report the most advanced plant RNA research presented in the 33rd International Conference on Arabidopsis Research (ICAR2023), held as a hybrid event on June 5-9, 2023 in Chiba, Japan, and focused on "Arabidopsis for Sustainable Development Goals". Six workshops/concurrent sessions in ICAR2023 targeted plant RNA biology, and many RNA-related topics could be found in other sessions. In this meeting report, we focus on the workshops/concurrent sessions targeting RNA biology, to share what is happening now at the forefront of plant RNA research.

A long non-coding RNA functions as a competitive endogenous RNA to modulate TaNAC018 by acting as a decoy for tae-miR6206

Increasing evidence indicates a strong correlation between the deposition of cuticular waxes and drought tolerance. However, the precise regulatory mechanism remains elusive. Here, we conducted a comprehensive transcriptome analysis of two wheat (Triticum aestivum) near-isogenic lines, the glaucous line G-JM38 rich in cuticular waxes and the non-glaucous line NG-JM31. We identified 85,143 protein-coding mRNAs, 4,485 lncRNAs, and 1,130 miRNAs. Using the lncRNA-miRNA-mRNA network and endogenous target mimic (eTM) prediction, we discovered that lncRNA35557 acted as an eTM for the miRNA tae-miR6206, effectively preventing tae-miR6206 from cleaving the NAC transcription factor gene TaNAC018. This lncRNA-miRNA interaction led to higher transcript abundance for TaNAC018 and enhanced drought-stress tolerance. Additionally, treatment with mannitol and abscisic acid (ABA) each influenced the levels of tae-miR6206, lncRNA35557, and TaNAC018 transcript. The ectopic expression of TaNAC018 in Arabidopsis also improved tolerance toward mannitol and ABA treatment, whereas knocking down TaNAC018 transcript levels via virus-induced gene silencing in wheat rendered seedlings more sensitive to mannitol stress. Our results indicate that lncRNA35557 functions as a competing endogenous RNA to modulate TaNAC018 expression by acting as a decoy target for tae-miR6206 in glaucous wheat, suggesting that non-coding RNA has important roles in the regulatory mechanisms responsible for wheat stress tolerance.

RNAvigate: efficient exploration of RNA chemical probing datasets

Chemical probing technologies enable high-throughput examination of diverse structural features of RNA, including local nucleotide flexibility, RNA secondary structure, protein and ligand binding, through-space interaction networks, and multistate structural ensembles. Deep understanding of RNA structure-function relationships typically requires evaluating a system under structure- and function-altering conditions, linking these data with additional information, and visualizing multilayered relationships. Current platforms lack the broad accessibility, flexibility and efficiency needed to iterate on integrative analyses of these diverse, complex data. Here, we share the RNA visualization and graphical analysis toolset RNAvigate, a straightforward and flexible Python library that automatically parses 21 standard file formats (primary sequence annotations, per- and internucleotide data, and secondary and tertiary structures) and outputs 18 plot types. RNAvigate enables efficient exploration of nuanced relationships between multiple layers of RNA structure information and across multiple experimental conditions. Compatibility with Jupyter notebooks enables nonburdensome, reproducible, transparent and organized sharing of multistep analyses and data visualization strategies. RNAvigate simplifies and accelerates discovery and characterization of RNA-centric functions in biology.

Analysis of Telomerase RNA Structure in Physcomitrium patens Indicates Functionally Relevant Transitions Between OPEN and CLOSED Conformations

Telomerase RNA (TR) conformation determines its function as a template for telomere synthesis and as a scaffold for the assembly of the telomerase nucleoprotein complex. Experimental analyses of TR secondary structure using DMS-Map Seq and SHAPE-Map Seq techniques show its CLOSED conformation as the consensus structure where the template region cannot perform its function. Our data show that the apparent discrepancy between experimental results and predicted TR functional conformation, mostly ignored in published studies, can be explained using data analysis based on single-molecule structure prediction from individual sequencing reads by the recently established DaVinci method. This method results in several clusters of secondary structures reflecting the structural dynamics of TR, possibly related to its multiple functional states. Interestingly, the presumed active (OPEN) conformation of TR corresponds to a minor fraction of TR under in vivo conditions. Therefore, structural polymorphism and dynamic TR transitions between CLOSED and OPEN conformations may be involved in telomerase activity regulation as a switch that functions independently of total TR transcript levels.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

2'-OH as a universal handle for studying intracellular RNAs

RNA plays pivotal roles in most cellular processes, serving as both the traditional carrier of genetic information and as a key regulator of cellular functions. The advent of chemical technologies has contributed critically to the analysis of cellular RNA structures, functions, and interactions. Many of these methods and molecules involve the utilization of chemically reactive handles in RNAs, either introduced externally or inherent within the polymer itself. Among these handles, the 2'-hydroxyl (2'-OH) group has emerged as an exceptionally well-suited and general chemical moiety for the modification and profiling of RNAs in intracellular studies. In this review, we provide an overview of the recent advancements in intracellular applications of acylation at the 2'-OH group of RNA. We outline progress made in probing RNA structure and interactomes, controlling RNA function, RNA imaging, and analyzing RNA-small molecule interactions, all achieved in living cells through this simple chemical handle on the biopolymer.

G·U base pairing motifs in long non-coding RNAs

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts involved in gene expression regulation and associated with diseases. Despite the unprecedented molecular complexity of these transcripts, recent studies of the secondary and tertiary structure of lncRNAs are starting to reveal the principles of lncRNA structural organization, with important functional implications. It therefore starts to be possible to analyze lncRNA structures systematically. Here, using a set of prototypical and medically-relevant lncRNAs of known secondary structure, we specifically catalogue the distribution and structural environment of one of the first-identified and most frequently occurring non-canonical Watson-Crick interactions, the G·U base pair. We compare the properties of G·U base pairs in our set of lncRNAs to those of the G·U base pairs in other well-characterized transcripts, like rRNAs, tRNAs, ribozymes, and riboswitches. Furthermore, we discuss how G·U base pairs in these targets participate in establishing interactions with proteins or miRNAs, and how they enable lncRNA tertiary folding by forming intramolecular or metal-ion interactions. Finally, by identifying highly-G·U-enriched regions of yet unknown function in our target lncRNAs, we provide a new rationale for future experimental investigation of these motifs, which will help obtain a more comprehensive understanding of lncRNA functions and molecular mechanisms in the future.

Is winter coming? Impact of the changing climate on plant responses to cold temperature

Climate change is causing alterations in annual temperature regimes worldwide. Important aspects of this include the reduction of winter chilling temperatures as well as the occurrence of unpredicted frosts, both significantly affecting plant growth and yields. Recent studies advanced the knowledge of the mechanisms underlying cold responses and tolerance in the model plant Arabidopsis thaliana. However, how these cold-responsive pathways will readjust to ongoing seasonal temperature variation caused by global warming remains an open question. In this review, we highlight the plant developmental programmes that depend on cold temperature. We focus on the molecular mechanisms that plants have evolved to adjust their development and stress responses upon exposure to cold. Covering both genetic and epigenetic aspects, we present the latest insights into how alternative splicing, noncoding RNAs and the formation of biomolecular condensates play key roles in the regulation of cold responses. We conclude by commenting on attractive targets to accelerate the breeding of increased cold tolerance, bringing up biotechnological tools that might assist in overcoming current limitations. Our aim is to guide the reflection on the current agricultural challenges imposed by a changing climate and to provide useful information for improving plant resilience to unpredictable cold regimes.

Causal role of a promoter polymorphism in natural variation of the Arabidopsis floral repressor gene FLC

Noncoding polymorphism frequently associates with phenotypic variation, but causation and mechanism are rarely established. Noncoding single-nucleotide polymorphisms (SNPs) characterize the major haplotypes of the Arabidopsis thaliana floral repressor gene FLOWERING LOCUS C (FLC). This noncoding polymorphism generates a range of FLC expression levels, determining the requirement for and the response to winter cold. The major adaptive determinant of these FLC haplotypes was shown to be the autumnal levels of FLC expression. Here, we investigate how noncoding SNPs influence FLC transcriptional output. We identify an upstream transcription start site (uTSS) cluster at FLC, whose usage is increased by an A variant at the promoter SNP-230. This variant is present in relatively few Arabidopsis accessions, with the majority containing G at this site. We demonstrate a causal role for the A variant at -230 in reduced FLC transcriptional output. The G variant upregulates FLC expression redundantly with the major transcriptional activator FRIGIDA (FRI). We demonstrate an additive interaction of SNP-230 with an intronic SNP+259, which also differentially influences uTSS usage. Combinatorial interactions between noncoding SNPs and transcriptional activators thus generate quantitative variation in FLC transcription that has facilitated the adaptation of Arabidopsis accessions to distinct climates.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

Emerging functional principles of tRNA-derived small RNAs and other regulatory small RNAs

Recent advancements in small RNA sequencing have unveiled a previously hidden world of regulatory small noncoding RNAs (sncRNAs) that extend beyond the well-studied small interfering RNAs, microRNAs, and piwi-interacting RNAs. This exploration, starting with tRNA-derived small RNAs, has led to the discovery of a diverse universe of sncRNAs derived from various longer structured RNAs such as rRNAs, small nucleolar RNAs, small nuclear RNAs, Y RNAs, and vault RNAs, with exciting uncharted functional possibilities. In this perspective, we discuss the emerging functional principles of sncRNAs beyond the well-known RNAi-like mechanisms, focusing on those that operate independent of linear sequence complementarity but rather function in an aptamer-like fashion. Aptamers use 3D structure for specific interactions with ligands and are modulated by RNA modifications and subcellular environments. Given that aptamer-like sncRNA functions are widespread and present in species lacking RNAi, they may represent an ancient functional principle that predates RNAi. We propose a rethinking of the origin of RNAi and its relationship with these aptamer-like functions in sncRNAs and how these complementary mechanisms shape biological processes. Lastly, the aptamer-like function of sncRNAs highlights the need for caution in using small RNA mimics in research and therapeutics, as their specificity is not restricted solely to linear sequence.

Deciphering shared attributes of plant long non-coding RNAs through a comparative computational approach

Over the past decade, long non-coding RNA (lncRNA), which lacks protein-coding potential, has emerged as an essential regulator of the genome. The present study examined 13,599 lncRNAs in Arabidopsis thaliana, 11,565 in Oryza sativa, and 32,397 in Zea mays for their characteristic features and explored the associated genomic and epigenomic features. We found lncRNAs were distributed throughout the chromosomes and the Helitron family of transposable elements (TEs) enriched, while the terminal inverted repeat depleted in lncRNA transcribing regions. Our analyses determined that lncRNA transcribing regions show rare or weak signals for most epigenetic marks except for H3K9me2 and cytosine methylation in all three plant species. LncRNAs showed preferential localization in the nucleus and cytoplasm; however, the distribution ratio in the cytoplasm and nucleus varies among the studied plant species. We identified several conserved endogenous target mimic sites in the lncRNAs among the studied plants. We found 233, 301, and 273 unique miRNAs, potentially targeting the lncRNAs of A. thaliana, O. sativa, and Z. mays, respectively. Our study has revealed that miRNAs, which interact with lncRNAs, target genes that are involved in a diverse array of biological and molecular processes. The miRNA-targeted lncRNAs displayed a strong affinity for several transcription factors, including ERF and BBR-BPC, mutually present in all three plants, advocating their conserved functions. Overall, the present study showed that plant lncRNAs exhibit conserved genomic and epigenomic characteristics and potentially govern the growth and development of plants.

Mutation signature filtering enables high-fidelity RNA structure probing at all four nucleobases with DMS

Chemical probing experiments have transformed RNA structure analysis, enabling high-throughput measurement of base-pairing in living cells. Dimethyl sulfate (DMS) is one of the most widely used structure probing reagents and has played a pivotal role in enabling next-generation single-molecule probing analyses. However, DMS has traditionally only been able to probe adenine and cytosine nucleobases. We previously showed that, using appropriate conditions, DMS can also be used to interrogate base-pairing of uracil and guanines in vitro at reduced accuracy. However, DMS remained unable to informatively probe guanines in cells. Here, we develop an improved DMS mutational profiling (MaP) strategy that leverages the unique mutational signature of N1-methylguanine DMS modifications to enable high-fidelity structure probing at all four nucleotides, including in cells. Using information theory, we show that four-base DMS reactivities convey greater structural information than current two-base DMS and SHAPE probing strategies. Four-base DMS experiments further enable improved direct base-pair detection by single-molecule PAIR analysis, and ultimately support RNA structure modeling at superior accuracy. Four-base DMS probing experiments are straightforward to perform and will broadly facilitate improved RNA structural analysis in living cells.

RNA structure determination: From 2D to 3D

RNA molecules serve a wide range of functions that are closely linked to their structures. The basic structural units of RNA consist of single- and double-stranded regions. In order to carry out advanced functions such as catalysis and ligand binding, certain types of RNAs can adopt higher-order structures. The analysis of RNA structures has progressed alongside advancements in structural biology techniques, but it comes with its own set of challenges and corresponding solutions. In this review, we will discuss recent advances in RNA structure analysis techniques, including structural probing methods, X-ray crystallography, nuclear magnetic resonance, cryo-electron microscopy, and small-angle X-ray scattering. Often, a combination of multiple techniques is employed for the integrated analysis of RNA structures. We also survey important RNA structures that have been recently determined using various techniques.

Long noncoding RNAs contribute to DNA damage resistance in Arabidopsis thaliana

Efficient repair of DNA lesions is essential for the faithful transmission of genetic information between somatic cells and for genome integrity across generations. Plants have multiple, partially redundant, and overlapping DNA repair pathways, probably due to the less constricted germline and the inevitable exposure to light including higher energy wavelengths. Many proteins involved in DNA repair and their mode of actions are well described. In contrast, a role for DNA damage-associated RNA components, evident from many other organisms, is less well understood. Here, we have challenged young Arabidopsis thaliana plants with two different types of genotoxic stress and performed de novo assembly and transcriptome analysis. We identified three long noncoding RNAs (lncRNAs) that are lowly or not expressed under regular conditions but up-regulated or induced by DNA damage. We generated CRISPR/Cas deletion mutants and found that the absence of the lncRNAs impairs the recovery capacity of the plants from genotoxic stress. The genetic loci are highly conserved among world-wide distributed Arabidopsis accessions and within related species in the Brassicaceae group. Together, these results suggest that the lncRNAs have a conserved function in connection with DNA damage and provide a basis for mechanistic analysis of their role.

Long Noncoding RNAs in Plant-Pathogen Interactions

Long noncoding RNAs (lncRNAs) are commonly defined as transcripts that lack protein-coding capacity and are longer than 200 nucleotides. Since the emergence of next-generation sequencing technologies in this century, thousands of lncRNAs have been identified from nearly all living organisms. Notably, various pathogens also express their own lncRNAs in host cells during infection. In plants, many lncRNAs exhibit dynamic expression patterns in response to environmental stimuli, including pathogen attacks. In contrast to well-established methods in identifying such lncRNAs, the current understanding of lncRNAs' functional mechanisms is in its infancy. Some lncRNAs serve as precursors for generating small RNAs or serve as target mimics to sequester functional small RNAs, which have been extensively reviewed in the literature. This review focuses on the emerging evidence supporting that certain lncRNAs function as negative or positive regulators of plant immunity. A common theme is that those regulations rely on specific interactions between lncRNAs and key regulatory proteins. Viroids as single-stranded circular noncoding RNAs provide a handle to investigate how RNA local motifs render interaction specificity between lncRNAs and regulatory proteins. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

RNA covariation at helix-level resolution for the identification of evolutionarily conserved RNA structure

Many biologically important RNAs fold into specific 3D structures conserved through evolution. Knowing when an RNA sequence includes a conserved RNA structure that could lead to new biology is not trivial and depends on clues left behind by conservation in the form of covariation and variation. For that purpose, the R-scape statistical test was created to identify from alignments of RNA sequences, the base pairs that significantly covary above phylogenetic expectation. R-scape treats base pairs as independent units. However, RNA base pairs do not occur in isolation. The Watson-Crick (WC) base pairs stack together forming helices that constitute the scaffold that facilitates the formation of the non-WC base pairs, and ultimately the complete 3D structure. The helix-forming WC base pairs carry most of the covariation signal in an RNA structure. Here, I introduce a new measure of statistically significant covariation at helix-level by aggregation of the covariation significance and covariation power calculated at base-pair-level resolution. Performance benchmarks show that helix-level aggregated covariation increases sensitivity in the detection of evolutionarily conserved RNA structure without sacrificing specificity. This additional helix-level sensitivity reveals an artifact that results from using covariation to build an alignment for a hypothetical structure and then testing the alignment for whether its covariation significantly supports the structure. Helix-level reanalysis of the evolutionary evidence for a selection of long non-coding RNAs (lncRNAs) reinforces the evidence against these lncRNAs having a conserved secondary structure.

Target Landscape of Conserved Plant MicroRNAs and the Complexities of Their Ancient MicroRNA-Binding Sites

In plants, microRNA (miRNA)-target interactions (MTIs) require high complementarity, a feature from which bioinformatic programs have predicted numerous and diverse targets for any given miRNA, promoting the idea of complex miRNA networks. Opposing this is a hypothesis of constrained miRNA specificity, in which functional MTIs are restricted to the few targets whose required expression output is compatible with the expression of the miRNA. To explore these opposing views, the bioinformatic pipeline Targets Ranked Using Experimental Evidence was applied to strongly conserved miRNAs to identity their high-evidence (HE) targets across species. For each miRNA family, HE targets predominantly consisted of homologs from one conserved target gene family (primary family). These primary families corresponded to the known canonical miRNA-target families, validating the approach. Very few additional HE target families were identified (secondary family), and if so, they were likely functionally related to the primary family. Many primary target families contained highly conserved nucleotide sequences flanking their miRNA-binding sites that were enriched in HE homologs across species. A number of these flanking sequences are predicted to form conserved RNA secondary structures that preferentially base pair with the miRNA-binding site, implying that these sites are highly structured. Our findings support a target landscape view that is dominated by the conserved primary target families, with a minority of either secondary target families or non-conserved targets. This is consistent with the constrained hypothesis of functional miRNA specificity, which potentially in part is being facilitated by features beyond complementarity.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.