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

Untranslated yet indispensable-UTRs act as key regulators in the environmental control of gene expression

To survive and thrive in a dynamic environment, plants must continuously monitor their surroundings and adjust their development and physiology accordingly. Changes in gene expression underlie these developmental and physiological adjustments, and are traditionally attributed to widespread transcriptional reprogramming. Growing evidence, however, suggests that post-transcriptional mechanisms also play a vital role in tailoring gene expression to a plant's environment. Untranslated regions (UTRs) act as regulatory hubs for post-transcriptional control, harbouring cis-elements that affect an mRNA's processing, localization, translation, and stability, and thereby tune the abundance of the encoded protein. Here, we review recent advances made in understanding the critical function UTRs exert in the post-transcriptional control of gene expression in the context of a plant's abiotic environment. We summarize the molecular mechanisms at play, present examples of UTR-controlled signalling cascades, and discuss the potential that resides within UTRs to render plants more resilient to a changing climate.

Transcriptomic analysis reveals candidate genes associated with salinity stress tolerance during the early vegetative stage in fababean genotype, Hassawi-2

Abiotic stresses are a significant constraint to plant production globally. Identifying stress-related genes can aid in the development of stress-tolerant elite genotypes and facilitate trait and crop manipulation. The primary aim of this study was to conduct whole transcriptome analyses of the salt-tolerant faba bean genotype, Hassawi-2, under different durations of salt stress (6 h, 12 h, 24 h, 48 h, and 72 h) at the early vegetative stage, to better understand the molecular basis of salt tolerance. After de novo assembly, a total of 140,308 unigenes were obtained. The up-regulated differentially expressed genes (DEGs) were 2380, 2863, 3057, 3484, and 4820 at 6 h, 12 h, 24 h, 48 h, and 72 h of salt stress, respectively. Meanwhile, 1974, 3436, 2371, 3502, and 5958 genes were downregulated at 6 h, 12 h, 24 h, 48 h, and 72 h of salt stress, respectively. These DEGs encoded various regulatory and functional proteins, including kinases, plant hormone proteins, transcriptional factors (TFs) basic helix-loop-helix (bHLH), Myeloblastosis (MYB), and (WRKY), heat shock proteins (HSPs), late embryogenesis abundant (LEA) proteins, dehydrin, antioxidant enzymes, and aquaporin proteins. This suggests that the faba bean genome possesses an abundance of salinity resistance genes, which trigger different adaptive mechanisms under salt stress. Some selected DEGs validated the RNA sequencing results, thus confirming similar gene expression levels. This study represents the first transcriptome analysis of faba bean leaves subjected to salinity stress offering valuable insights into the mechanisms governing salt tolerance in faba bean during the vegetative stage. This comprehensive investigation enhances our understanding of precise gene regulatory mechanisms and holds promise for the development of novel salt-tolerant faba bean salt-tolerant cultivars.

In vivo structure probing of RNA in Archaea: novel insights into the ribosome structure of Methanosarcina acetivorans

Structure probing combined with next-generation sequencing (NGS) has provided novel insights into RNA structure-function relationships. To date, such studies have focused largely on bacteria and eukaryotes, with little attention given to the third domain of life, archaea. Furthermore, functional RNAs have not been extensively studied in archaea, leaving open questions about RNA structure and function within this domain of life. With archaeal species being diverse and having many similarities to both bacteria and eukaryotes, the archaea domain has the potential to be an evolutionary bridge. In this study, we introduce a method for probing RNA structure in vivo in the archaea domain of life. We investigated the structure of ribosomal RNA (rRNA) from Methanosarcina acetivorans, a well-studied anaerobic archaeal species, grown with either methanol or acetate. After probing the RNA in vivo with dimethyl sulfate (DMS), Structure-seq2 libraries were generated, sequenced, and analyzed. We mapped the reactivity of DMS onto the secondary structure of the ribosome, which we determined independently with comparative analysis, and confirmed the accuracy of DMS probing in M. acetivorans Accessibility of the rRNA to DMS in the two carbon sources was found to be quite similar, although some differences were found. Overall, this study establishes the Structure-seq2 pipeline in the archaea domain of life and informs about ribosomal structure within M. acetivorans.

Transcriptome-wide probing reveals RNA thermometers that regulate translation of glycerol permease genes in Bacillus subtilis

RNA structure regulates bacterial gene expression by several distinct mechanisms via environmental and cellular stimuli, one of which is temperature. While some genome-wide studies have focused on heat shock treatments and the subsequent transcriptomic changes, soil bacteria are less likely to experience such rapid and extreme temperature changes. Though RNA thermometers (RNATs) have been found in 5' untranslated leader regions (5' UTRs) of heat shock and virulence-associated genes, this RNA-controlled mechanism could regulate other genes as well. Using Structure-seq2 and the chemical probe dimethyl sulfate (DMS) at four growth temperatures ranging from 23°C to 42°C, we captured a dynamic response of the Bacillus subtilis transcriptome to temperature. Our transcriptome-wide results show RNA structural changes across all four temperatures and reveal nonmonotonic reactivity trends with increasing temperature. Then, focusing on subregions likely to contain regulatory RNAs, we examined 5' UTRs to identify large, local reactivity changes. This approach led to the discovery of RNATs that control the expression of glpF (glycerol permease) and glpT (glycerol-3-phosphate permease); expression of both genes increased with increased temperature. Results with mutant RNATs indicate that both genes are controlled at the translational level. Increased import of glycerols at high temperatures could provide thermoprotection to proteins.

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.

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.

Biological solution conditions and flanking sequence modulate LLPS of RNA G-quadruplex structures

Guanine-rich regions of DNA or RNA can form structures with two or more consecutive G-quartets called G-quadruplexes (GQ). Recent studies reveal the potential for these structures to aggregate in vitro. Here, we report effects of in vivo concentrations of additives-amino acids, nucleotides, and crowding agents-on the structure and solution behavior of RNAs containing GQ-forming sequences. We found that cytosine nucleotides destabilize a model GQ structure at biological salt concentrations, while free amino acids and other nucleotides do not do so to a substantial degree. We also report that the tendency of folded GQs to form droplets or to aggregate depends on the nature of flanking sequence and the presence of additives. Notably, in the presence of biological amounts of polyamines, flanking regions on the 5'-end of the RNA drive more droplet-like phase separation, while flanking regions on the 3'-end, as well as both the 5'- and 3'-ends, induce more condensed, granular structures. Finally, we provide an example of a biological sequence in the presence of polyamines and show that crowders such as PEG and dextran can selectively cause its phase separation. These findings have implications for the participation of GQS in LLPS in vivo.

RNA biology takes root in plant systems

Advances in RNA biology such as RNAi, CRISPR, and the first mRNA vaccine represent the enormous potential of RNA research to address current problems. Additionally, plants are a diverse and undeniably essential resource for life threatened by climate change, loss of arable land, and pollution. Different aspects of RNA such as its processing, modification and structure are intertwined with plant development, physiology and stress response. This report details the findings of researchers around the world during the 23rd Penn State Symposium in Plant Biology with a focus in RNA biology.

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.

Genome-wide analysis of the in vivo tRNA structurome reveals RNA structural and modification dynamics under heat stress

RNA structure plays roles in myriad cellular events including transcription, translation, and RNA processing. Genome-wide analyses of RNA secondary structure in vivo by chemical probing have revealed critical structural features of mRNAs and long ncRNAs. Here, we examine the in vivo secondary structure of a small RNA class, tRNAs. Study of tRNA structure is challenging because tRNAs are heavily modified and strongly structured. We introduce "tRNA structure-seq," a new workflow that accurately determines in vivo secondary structures of tRNA. The workflow combines dimethyl sulfate (DMS) probing, ultra-processive RT, and mutational profiling (MaP), which provides mutations opposite DMS and natural modifications thereby allowing multiple modifications to be identified in a single read. We applied tRNA structure-seq to E. coli under control and stress conditions. A leading folding algorithm predicts E. coli tRNA structures with only ∼80% average accuracy from sequence alone. Strikingly, tRNA structure-seq, by providing experimental restraints, improves structure prediction under in vivo conditions to ∼95% accuracy, with more than 14 tRNAs predicted completely correctly. tRNA structure-seq also quantifies the relative levels of tRNAs and their natural modifications at single nucleotide resolution, as validated by LC-MS/MS. Our application of tRNA structure-seq yields insights into tRNA structure in living cells, revealing that it is not immutable but has dynamics, with partial unfolding of secondary and tertiary tRNA structure under heat stress that is correlated with a loss of tRNA abundance. This method is applicable to other small RNAs, including those with natural modifications and highly structured regions.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

Background

Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”

Results

We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.

Conclusion

We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02656-4.

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.

Global N6-Methyladenosine Profiling Revealed the Tissue-Specific Epitranscriptomic Regulation of Rice Responses to Salt Stress

N6-methyladenosine (m6A) methylation represents a new layer of the epitranscriptomic regulation of plant development and growth. However, the effects of m6A on rice responses to environmental stimuli remain unclear. In this study, we performed a methylated-RNA immunoprecipitation sequencing analysis and compared the changes in m6A methylation and gene expression in rice under salt stress conditions. Salt stress significantly increased the m6A methylation in the shoots (p value < 0.05). Additionally, 2537 and 2304 differential m6A sites within 2134 and 1997 genes were identified in the shoots and roots, respectively, under salt stress and control conditions. These differential m6A sites were largely regulated in a tissue-specific manner. A unique set of genes encoding transcription factors, antioxidants, and auxin-responsive proteins had increased or decreased m6A methylation levels only in the shoots or roots under salt stress, implying m6A may mediate salt tolerance by regulating transcription, ROS homeostasis, and auxin signaling in a tissue-specific manner. Integrating analyses of m6A modifications and gene expression changes revealed that m6A changes regulate the expression of genes controlling plant growth, stress responses, and ion transport under saline conditions. These findings may help clarify the regulatory effects of m6A modifications on rice salt tolerance.

RNA multimerization as an organizing force for liquid-liquid phase separation

RNA interactions are exceptionally strong and highly redundant. As such, nearly any two RNAs have the potential to interact with one another over relatively short stretches, especially at high RNA concentrations. This is especially true for pairs of RNAs that do not form strong self-structure. Such phenomena can drive liquid-liquid phase separation, either solely from RNA-RNA interactions in the presence of divalent or organic cations, or in concert with proteins. RNA interactions can drive multimerization of RNA strands via both base-pairing and tertiary interactions. In this article, we explore the tendency of RNA to form stable monomers, dimers, and higher order structures as a function of RNA length and sequence through a focus on the intrinsic thermodynamic, kinetic, and structural properties of RNA. The principles we discuss are independent of any specific type of biomolecular condensate, and thus widely applicable. We also speculate how external conditions experienced by living organisms can influence the formation of nonmembranous compartments, again focusing on the physical and structural properties of RNA. Plants, in particular, are subject to diverse abiotic stresses including extreme temperatures, drought, and salinity. These stresses and the cellular responses to them, including changes in the concentrations of small molecules such as polyamines, salts, and compatible solutes, have the potential to regulate condensate formation by melting or strengthening base-pairing. Reversible condensate formation, perhaps including regulation by circadian rhythms, could impact biological processes in plants, and other organisms.

Studies on Viroid Shed Light on the Role of RNA Three-Dimensional Structural Motifs in RNA Trafficking in Plants

RNAs play essential roles in various biological processes. Mounting evidence has demonstrated that RNA subcellular localization and intercellular/systemic trafficking govern their functions in coordinating plant growth at the organismal level. While numerous types of RNAs (i.e., mRNAs, small RNAs, rRNAs, tRNAs, and long noncoding RNAs) have been found to traffic in a non-cell-autonomous fashion within plants, the underlying regulatory mechanism remains unclear. Viroids are single-stranded circular noncoding RNAs, which entirely rely on their RNA motifs to exploit cellular machinery for organelle entry and exit, cell-to-cell movement through plasmodesmata, and systemic trafficking. Viroids represent an excellent model to dissect the role of RNA three-dimensional (3D) structural motifs in regulating RNA movement. Nearly two decades of studies have found multiple RNA 3D motifs responsible for viroid nuclear import as well as trafficking across diverse cellular boundaries in plants. These RNA 3D motifs function as "keys" to unlock cellular and subcellular barriers and guide RNA movement within a cell or between cells. Here, we summarize the key findings along this line of research with implications for future studies on RNA trafficking in plants.

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.

Long Tracts of Guanines Drive Aggregation of RNA G-Quadruplexes in the Presence of Spermine

G-Quadruplexes (GQs) are compact, stable structures in DNA and RNA comprised of two or more tiers of quartets whose G-rich motif of tracts of two or more G's occurs commonly within genomes and transcriptomes. While thermodynamically stable in vitro, these structures remain difficult to study in vivo. One approach to understanding GQ in vivo behavior is to test whether conditions and molecules found in cells facilitate their folding. Polyamines are biogenic polycations that interact with RNA. Among common polyamines, spermine contains the highest charge and is found in eukaryotes, making it a good candidate for association with high-charge density nucleic acid structures like GQs. Using a variety of techniques, including ultraviolet-detected thermal denaturation, circular dichroism, size exclusion chromatography, and confocal microscopy, on an array of quadruplex sequence variants, we find that eukaryotic biological concentrations of spermine induce microaggregation of three-tiered G-rich sequences, but not of purely two-tiered structures, although higher spermine concentrations induce aggregation of even these. The formation of microaggregates can also be induced by addition of as little as a single G to a two-tiered structure; moreover, they form at biological temperatures, are sensitive to salt, and can form in the presence of at least some flanking sequence. Notably, GQ aggregation is not observed under prokaryotic-like conditions of no spermine and higher NaCl concentrations. The sequence, polyamine, and salt specificity of microaggregation reported herein have implications for the formation and stability of G-rich nucleic acid aggregates in vivo and for functional roles for understudied GQ sequences with only two quadruplex tiers.

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

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

RNA 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.

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

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

N6-methyladenosine and RNA secondary structure affect transcript stability and protein abundance during systemic salt stress in Arabidopsis

After transcription, a messenger RNA (mRNA) is further post-transcriptionally regulated by several features including RNA secondary structure and covalent RNA modifications (specifically N6-methyladenosine, m6A). Both RNA secondary structure and m6A have been demonstrated to regulate mRNA stability and translation and have been independently linked to plant responses to soil salinity levels. However, the effect of m6A on regulating RNA secondary structure and the combinatorial interplay between these two RNA features during salt stress response has yet to be studied. Here, we globally identify RNA-protein interactions and RNA secondary structure during systemic salt stress. This analysis reveals that RNA secondary structure changes significantly during salt stress, and that it is independent of global changes in RNA-protein interactions. Conversely, we find that m6A is anti-correlated with RNA secondary structure in a condition-dependent manner, with salt-specific m6A correlated with a decrease in mRNA secondary structure during salt stress. Taken together, we suggest that salt-specific m6A deposition and the associated loss of RNA secondary structure results in increases in mRNA stability for transcripts encoding abiotic stress response proteins and ultimately increases in protein levels from these stabilized transcripts. In total, our comprehensive analyses reveal important post-transcriptional regulatory mechanisms involved in plant long-term salt stress response and adaptation.