Genome-wide probing RNA structure with the modified DMS-MaPseq in Arabidopsis

Parallel degradome-seq and DMS-MaPseq substantially revise the miRNA biogenesis atlas in Arabidopsis

MicroRNAs (miRNAs) are produced from highly structured primary transcripts (pri-miRNAs) and regulate numerous biological processes in eukaryotes. Due to the extreme heterogeneity of these structures, the initial processing sites of plant pri-miRNAs and the structural rules that determine their processing have been predicted for many miRNAs but remain elusive for others. Here we used semi-active DCL1 mutants and advanced degradome-sequencing strategies to accurately identify the initial processing sites for 147 of 326 previously annotated Arabidopsis miRNAs and to illustrate their associated pri-miRNA cleavage patterns. Elucidating the in vivo RNA secondary structures of 73 pri-miRNAs revealed that about 95% of them differ from in silico predictions, and that the revised structures offer clearer interpretation of the processing sites and patterns. Finally, DCL1 partners Serrate and HYL1 could synergistically and independently impact processing patterns and in vivo RNA secondary structures of pri-miRNAs. Together, our work sheds light on the precise processing mechanisms of plant pri-miRNAs.

Structured 3' UTRs destabilize mRNAs in plants

BACKGROUND

RNA secondary structure (RSS) can influence the regulation of transcription, RNA processing, and protein synthesis, among other processes. 3' untranslated regions (3' UTRs) of mRNA also hold the key for many aspects of gene regulation. However, there are often contradictory results regarding the roles of RSS in 3' UTRs in gene expression in different organisms and/or contexts.

RESULTS

Here, we incidentally observe that the primary substrate of miR159a (pri-miR159a), when embedded in a 3' UTR, could promote mRNA accumulation. The enhanced expression is attributed to the earlier polyadenylation of the transcript within the hybrid pri-miR159a-3' UTR and, resultantly, a poorly structured 3' UTR. RNA decay assays indicate that poorly structured 3' UTRs could promote mRNA stability, whereas highly structured 3' UTRs destabilize mRNA in vivo. Genome-wide DMS-MaPseq also reveals the prevailing inverse relationship between 3' UTRs' RSS and transcript accumulation in the transcriptomes of Arabidopsis, rice, and even human. Mechanistically, transcripts with highly structured 3' UTRs are preferentially degraded by 3'-5' exoribonuclease SOV and 5'-3' exoribonuclease XRN4, leading to decreased expression in Arabidopsis. Finally, we engineer different structured 3' UTRs to an endogenous FT gene and alter the FT-regulated flowering time in Arabidopsis.

CONCLUSIONS

We conclude that highly structured 3' UTRs typically cause reduced accumulation of the harbored transcripts in Arabidopsis. This pattern extends to rice and even mammals. Furthermore, our study provides a new strategy of engineering the 3' UTRs' RSS to modify plant traits in agricultural production and mRNA stability in biotechnology.

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.

Rock, scissors, paper: How RNA structure informs function

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

High-temperature adaptation of an OsNRT2.3 allele is thermoregulated by small RNAs

Climate change negatively affects crop yield, which hinders efforts to reach agricultural sustainability and food security. Here, we show that a previously unidentified allele of the nitrate transporter gene OsNRT2.3 is required to maintain high yield and high nitrogen use efficiency under high temperatures. We demonstrate that this tolerance to high temperatures in rice accessions harboring the HTNE-2 (high temperature resistant and nitrogen efficient-2) alleles from enhanced translation of the OsNRT2.3b mRNA isoform and the decreased abundance of a unique small RNA (sNRT2.3-1) derived from the 5' untranslated region of OsNRT2.3. sNRT2.3-1 binds to the OsNRT2.3a mRNA in a temperature-dependent manner. Our findings reveal that allelic variation in the 5' untranslated region of OsNRT2.3 leads to an increase in OsNRT2.3b protein levels and higher yield during high-temperature stress. Our results also provide a breeding strategy to produce rice varieties with higher grain yield and lower N fertilizer input suitable for a sustainable agriculture that is resilient against climate change.

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 Key Molecular Regulator, RNA G-Quadruplex and Its Function in Plants

RNA structure plays key roles in plant growth, development, and adaptation. One of the complex RNA structures is the RNA G-quadruplex (RG4) where guanine-rich sequences are folded into two or more layers of G-quartets. Previous computational predictions of RG4 revealed that it is widespread across the whole transcriptomes in many plant species, raising the hypothesis that RG4 is likely to be an important regulatory motif in plants. Recently, with the advances in both high-throughput sequencing and cell imaging technologies, RG4 can be detected in living cells as well as at the genome-wide scale. Here, we provide a comprehensive review of recent developments in new methods for detecting RG4 in plants. We also summarize the new functions of RG4 in regulating plant growth and development. We then discuss the possible role of RG4 in adapting to environmental conditions along with evolutionary perspectives.

Probing in vivo RNA Structure With Optimized DMS-MaPseq in Rice

RNA transcripts form various secondary and tertiary structures that have a wide range of regulatory functions. Several methods have been developed to profile in vivo RNA secondary structure in prokaryotes and eukaryotes. These methods, such as dimethyl sulfate (DMS) mutational profiling with high-throughput sequencing (DMS-MaPseq), couple small chemical-mediated RNA modifications with next-generation sequencing. DMS-MaPseq, a powerful method for genome-wide and target-specific RNA secondary structure profiling, has been applied in yeast, mammals, Drosophila, and Arabidopsis thaliana, but not in crops. Here, we used DMS-MaPseq to conduct a target-specific and genome-wide profile of in vivo RNA secondary structure in rice (Oryza sativa). The DMS treatment conditions were optimized for rice leaf and root tissues. To increase the sequencing depth and coverage of low-abundance transcripts in genome-wide DMS-MaPseq, we used streptavidin-biotin depletion to reduce the abundance of highly expressed chloroplast transcripts during library construction. The resulting target-specific and genome-wide rice DMS-MaPseq data were of high quality and reproducibility. Furthermore, we used DMS-MaPseq to profile the in vivo RNA secondary structure of an OsmiR399 target region located at 5'UTR of OsPHO2, which participates in rice phosphate homeostasis. An unfolded RNA structure downstream of miRNA target site was observed in predicted in vivo RNA secondary structure, reminiscence of the TAM (Target Adjacent nucleotide Motif) involved in mRNA structure-mediated regulation in miRNA cleavage. Our study optimized DMS-MaPseq for probing in vivo RNA secondary structure in rice, facilitating the study of RNA structure-mediated regulations in crops.

A natural single-nucleotide polymorphism variant in sulfite reductase influences sulfur assimilation in maize

Plants absorb sulfur from the environment and assimilate it into suitable forms for the biosynthesis of a broad range of molecules. Although the biochemical pathway of sulfur assimilation is known, how genetic differences contribute to natural variation in sulfur assimilation remains poorly understood. Here, using a genome-wide association study, we uncovered a single-nucleotide polymorphism (SNP) variant in the sulfite reductase (SiR) gene that was significantly associated with SiR protein abundance in a maize natural association population. We also demonstrated that the synonymous C to G base change at SNP69 may repress translational activity by altering messenger RNA secondary structure, which leads to reduction in ZmSiR protein abundance and sulfur assimilation activity. Population genetic analyses showed that the SNP69C allele was likely a variant occurring after the initial maize domestication and accumulated with the spread of maize cultivation from tropical to temperate regions. This study provides the first evidence that genetic polymorphisms in the exon of ZmSiR could influence the protein abundance through a posttranscriptional mechanism and in part contribute to natural variation in sulfur assimilation. These findings provide a prospective target to improve maize varieties with proper sulfur nutrient levels assisted by molecular breeding and engineering.

Arabidopsis retains vertebrate-type telomerase accessory proteins via a plant-specific assembly

The recent discovery of the bona-fide telomerase RNA (TR) from plants reveals conserved and unique secondary structure elements and the opportunity for new insight into the telomerase RNP. Here we examine how two highly conserved proteins previously implicated in Arabidopsis telomere maintenance, AtPOT1a and AtNAP57 (dyskerin), engage plant telomerase. We report that AtPOT1a associates with Arabidopsis telomerase via interaction with TERT. While loss of AtPOT1a does not impact AtTR stability, the templating domain is more accessible in pot1a mutants, supporting the conclusion that AtPOT1a stimulates telomerase activity but does not facilitate telomerase RNP assembly. We also show, that despite the absence of a canonical H/ACA binding motif within AtTR, dyskerin binds AtTR with high affinity and specificity in vitro via a plant specific three-way junction (TWJ). A core element of the TWJ is the P1a stem, which unites the 5′ and 3′ ends of AtTR. P1a is required for dyskerin-mediated stimulation of telomerase repeat addition processivity in vitro, and for AtTR accumulation and telomerase activity in vivo. The deployment of vertebrate-like accessory proteins and unique RNA structural elements by Arabidopsis telomerase provides a new platform for exploring telomerase biogenesis and evolution.

RNA architecture influences plant biology

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

Light-Dependent Translation Change of Arabidopsis psbA Correlates with RNA Structure Alterations at the Translation Initiation Region

mRNA secondary structure influences translation. Proteins that modulate the mRNA secondary structure around the translation initiation region may regulate translation in plastids. To test this hypothesis, we exposed Arabidopsis thaliana to high light, which induces translation of psbA mRNA encoding the D1 subunit of photosystem II. We assayed translation by ribosome profiling and applied two complementary methods to analyze in vivo RNA secondary structure: DMS-MaPseq and SHAPE-seq. We detected increased accessibility of the translation initiation region of psbA after high light treatment, likely contributing to the observed increase in translation by facilitating translation initiation. Furthermore, we identified the footprint of a putative regulatory protein in the 5' UTR of psbA at a position where occlusion of the nucleotide sequence would cause the structure of the translation initiation region to open up, thereby facilitating ribosome access. Moreover, we show that other plastid genes with weak Shine-Dalgarno sequences (SD) are likely to exhibit psbA-like regulation, while those with strong SDs do not. This supports the idea that changes in mRNA secondary structure might represent a general mechanism for translational regulation of psbA and other plastid genes.

RNA structure probing to characterize RNA-protein interactions on a low abundance pre-mRNA in living cells

In vivo RNA structure analysis has become a powerful tool in molecular biology, largely due to the coupling of an increasingly diverse set of chemical approaches with high-throughput sequencing. This has resulted in a transition from single target to transcriptome-wide approaches. However, these methods require sequencing depths that preclude studying low abundance targets, which are not sufficiently captured in transcriptome-wide approaches. Here we present a ligation-free method to enrich for low abundance RNA sequences, which improves the diversity of molecules analyzed and results in improved analysis. In addition, this method is compatible with any choice of chemical adduct or read-out approach. We utilized this approach to study an autoregulated event in the pre-mRNA of the splicing factor, muscleblind-like splicing regulator 1 (MBNL1).

RNA G-quadruplex structures exist and function in vivo in plants

BACKGROUND

Guanine-rich sequences are able to form complex RNA structures termed RNA G-quadruplexes in vitro. Because of their high stability, RNA G-quadruplexes are proposed to exist in vivo and are suggested to be associated with important biological relevance. However, there is a lack of direct evidence for RNA G-quadruplex formation in living eukaryotic cells. Therefore, it is unclear whether any purported functions are associated with the specific sequence content or the formation of an RNA G-quadruplex structure.

RESULTS

Using rG4-seq, we profile the landscape of those guanine-rich regions with the in vitro folding potential in the Arabidopsis transcriptome. We find a global enrichment of RNA G-quadruplexes with two G-quartets whereby the folding potential is strongly influenced by RNA secondary structures. Using in vitro and in vivo RNA chemical structure profiling, we determine that hundreds of RNA G-quadruplex structures are strongly folded in both Arabidopsis and rice, providing direct evidence of RNA G-quadruplex formation in living eukaryotic cells. Subsequent genetic and biochemical analyses show that RNA G-quadruplex folding is able to regulate translation and modulate plant growth.

CONCLUSIONS

Our study reveals the existence of RNA G-quadruplex in vivo and indicates that RNA G-quadruplex structures act as important regulators of plant development and growth.

Secondary Structure of Chloroplast mRNAs In Vivo and In Vitro

mRNA secondary structure can influence gene expression, e.g., by influencing translation initiation. The probing of in vivo mRNA secondary structures is therefore necessary to understand what determines the efficiency and regulation of gene expression. Here, in vivo mRNA secondary structure was analyzed using dimethyl sulfate (DMS)-MaPseq and compared to in vitro-folded RNA. We used an approach to analyze specific, full-length transcripts. To test this approach, we chose low, medium, and high abundant mRNAs. We included both monocistronic and multicistronic transcripts. Because of the slightly alkaline pH of the chloroplast stroma, we could probe all four nucleotides with DMS. The structural information gained was evaluated using the known structure of the plastid 16S rRNA. This demonstrated that the results obtained for adenosines and cytidines were more reliable than for guanosines and uridines. The majority of mRNAs analyzed were less structured in vivo than in vitro. The in vivo secondary structure of the translation initiation region of most tested genes appears to be optimized for high translation efficiency.