Global analysis of the RNA-protein interaction and RNA secondary structure landscapes of the Arabidopsis nucleus

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.

LncRNA DANA1 promotes drought tolerance and histone deacetylation of drought responsive genes in Arabidopsis

Although many long noncoding RNAs have been discovered in plants, little is known about their biological function and mode of action. Here we show that the drought-induced long intergenic noncoding RNA DANA1 interacts with the L1p/L10e family member protein DANA1-INTERACTING PROTEIN 1 (DIP1) in the cell nucleus of Arabidopsis, and both DANA1 and DIP1 promote plant drought resistance. DANA1 and DIP1 increase histone deacetylase HDA9 binding to the CYP707A1 and CYP707A2 loci. DIP1 further interacts with PWWP3, a member of the PEAT complex that associates with HDA9 and has histone deacetylase activity. Mutation of DANA1 enhances CYP707A1 and CYP707A2 acetylation and expression resulting in impaired drought tolerance, in agreement with dip1 and pwwp3 mutant phenotypes. Our results demonstrate that DANA1 is a positive regulator of drought response and that DANA1 works jointly with the novel chromatin-related factor DIP1 on epigenetic reprogramming of the plant transcriptome during the response to drought.

Long-range RNA structures in the human transcriptome beyond evolutionarily conserved regions

RNA structure has been increasingly recognized as a critical player in the biogenesis and turnover of many transcripts classes. In eukaryotes, the prediction of RNA structure by thermodynamic modeling meets fundamental limitations due to the large sizes and complex, discontinuous organization of eukaryotic genes. Signatures of functional RNA structures can be found by detecting compensatory substitutions in homologous sequences, but a comparative approach is applicable only within conserved sequence blocks. Here, we developed a computational pipeline called PHRIC, which is not limited to conserved regions and relies on RNA contacts derived from RNA in situ conformation sequencing (RIC-seq) experiments. It extracts pairs of short RNA fragments surrounded by nested clusters of RNA contacts and predicts long, nearly perfect complementary base pairings formed between these fragments. In application to a panel of RIC-seq experiments in seven human cell lines, PHRIC predicted ~12,000 stable long-range RNA structures with equilibrium free energy below -15 kcal/mol, the vast majority of which fall outside of regions annotated as conserved among vertebrates. These structures, nevertheless, show some level of sequence conservation and remarkable compensatory substitution patterns in other clades. Furthermore, we found that introns have a higher propensity to form stable long-range RNA structures between each other, and moreover that RNA structures tend to concentrate within the same intron rather than connect adjacent introns. These results for the first time extend the application of proximity ligation assays to RNA structure prediction beyond conserved regions.

Formation, function, and pathology of RNP granules

Ribonucleoprotein (RNP) granules are diverse membrane-less organelles that form through multivalent RNA-RNA, RNA-protein, and protein-protein interactions between RNPs. RNP granules are implicated in many aspects of RNA physiology, but in most cases their functions are poorly understood. RNP granules can be described through four key principles. First, RNP granules often arise because of the large size, high localized concentrations, and multivalent interactions of RNPs. Second, cells regulate RNP granule formation by multiple mechanisms including posttranslational modifications, protein chaperones, and RNA chaperones. Third, RNP granules impact cell physiology in multiple manners. Finally, dysregulation of RNP granules contributes to human diseases. Outstanding issues in the field remain, including determining the scale and molecular mechanisms of RNP granule function and how granule dysfunction contributes to human disease.

A systematic benchmark of machine learning methods for protein-RNA interaction prediction

RNA-binding proteins (RBPs) are central actors of RNA post-transcriptional regulation. Experiments to profile-binding sites of RBPs in vivo are limited to transcripts expressed in the experimental cell type, creating the need for computational methods to infer missing binding information. While numerous machine-learning based methods have been developed for this task, their use of heterogeneous training and evaluation datasets across different sets of RBPs and CLIP-seq protocols makes a direct comparison of their performance difficult. Here, we compile a set of 37 machine learning (primarily deep learning) methods for in vivo RBP-RNA interaction prediction and systematically benchmark a subset of 11 representative methods across hundreds of CLIP-seq datasets and RBPs. Using homogenized sample pre-processing and two negative-class sample generation strategies, we evaluate methods in terms of predictive performance and assess the impact of neural network architectures and input modalities on model performance. We believe that this study will not only enable researchers to choose the optimal prediction method for their tasks at hand, but also aid method developers in developing novel, high-performing methods by introducing a standardized framework for their evaluation.

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.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

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.

Genomic sequences and RNA binding proteins predict RNA splicing efficiency in various single-cell contexts

MOTIVATION

The RNA splicing efficiency is of high interest for both understanding the regulatory machinery of gene expression and estimating the RNA velocity in single cells. However, its genomic regulation and stochasticity across contexts remain poorly understood.

RESULTS

Here, by leveraging the recent RNA velocity tool, we estimated the relative splicing efficiency across a variety of single-cell RNA-Seq data sets. We further extracted large sets of genomic features and 120 RNA binding protein features and found they are highly predictive to relative RNA splicing efficiency across multiple tissues and organs on human and mouse. This predictive power brings promise to reveal the complexity of RNA processing and to enhance the analysis of single-cell transcription activities.

AVAILABILITY AND IMPLEMENTATION

In order to ensure reproducibility, all preprocessed data sets and scripts used for the prediction and figure generation are publicly available at https://doi.org/10.5281/zenodo.6513669.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pan-cancer analysis of necroptosis-related gene signature for the identification of prognosis and immune significance

BACKGROUND

Necroptosis is a novel programmed cell death mode independent on caspase. A number of studies have revealed that the induction of necroptosis could act as an alternative therapeutic strategy for drug-resistant tumors as well as affect tumor immune microenvironment.

METHODS

Gene expression profiles and clinical data were downloaded from XENA-UCSC (including The Cancer Genome Atlas and Genotype-Tissue Expression), Gene Expression Omnibus, International Cancer Genome Consortium and Chinese Glioma Genome Atlas. We used non-negative matrix factorization method to conduct tumor classification. The least absolute shrinkage and selection operator regression was applied to establish risk models, whose prognostic effectiveness was examined in both training and testing sets with Kaplan-Meier analysis, time-dependent receiver operating characteristic curves as well as uni- and multi-variate survival analysis. Principal Component Analysis, t-distributed Stochastic Neighbor Embedding and Uniform Manifold Approximation and Projection were conducted to check the risk group distribution. Gene Set Enrichment Analyses, immune infiltration analysis based on CIBERSORT, EPIC, MCPcounter, ssGSEA and ESTIMATE, gene mutation and drug sensitivity between the risk groups were also taken into consideration.

RESULTS

There were eight types of cancer with at least ten differentially expressed necroptosis-related genes which could influence patients' prognosis, namely, adrenocortical carcinoma (ACC), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), acute myeloid leukemia (LAML), brain lower grade glioma (LGG), pancreatic adenocarcinoma (PAAD), liver hepatocellular carcinoma (LIHC), skin cutaneous melanoma (SKCM) and thymoma (THYM). Patients could be divided into different clusters with distinct overall survival in all cancers above except for LIHC. The risk models could efficiently predict prognosis of ACC, LAML, LGG, LIHC, SKCM and THYM patients. LGG patients from high-risk group had a higher infiltration level of M2 macrophages and cancer-associated fibroblasts. There were more CD8+ T cells, Th1 cells and M1 macrophages in low-risk SKCM patients' tumor microenvironment. Gene mutation status and drug sensitivity are also different between low- and high-risk groups in the six cancers.

CONCLUSIONS

Necroptosis-related genes can predict clinical outcomes of ACC, LAML, LGG, LIHC, SKCM and THYM patients and help to distinguish immune infiltration status for LGG and SKCM.

High-Throughput Evolutionary Comparative Analysis of Long Intergenic Noncoding RNAs in Multiple Organisms

Comparative genomic and transcriptomic analyses can help prioritize and facilitate the functional analysis of long noncoding RNAs (lncRNAs). Evolinc-II is a bioinformatic pipeline that automates comparative analyses, searching for sequence and structural conservation for thousands of lncRNAs at once. In addition, Evolinc-II takes a phylogenetic approach to infer key evolutionary events that may have occurred during the emergence of each query lncRNA. Here, we describe how to use command line or GUI (CyVerse's Discovery Environment) versions of Evolinc-II to identify lncRNA homologs and prioritize them for functional analysis.

A Conserved Long Intergenic Non-coding RNA Containing snoRNA Sequences, lncCOBRA1, Affects Arabidopsis Germination and Development

Long non-coding RNAs (lncRNAs) are an increasingly studied group of non-protein coding transcripts with a wide variety of molecular functions gaining attention for their roles in numerous biological processes. Nearly 6,000 lncRNAs have been identified in Arabidopsis thaliana but many have yet to be studied. Here, we examine a class of previously uncharacterized lncRNAs termed CONSERVED IN BRASSICA RAPA (lncCOBRA) transcripts that were previously identified for their high level of sequence conservation in the related crop species Brassica rapa, their nuclear-localization and protein-bound nature. In particular, we focus on lncCOBRA1 and demonstrate that its abundance is highly tissue and developmental specific, with particularly high levels early in germination. lncCOBRA1 contains two snoRNAs domains within it, making it the first sno-lincRNA example in a non-mammalian system. However, we find that it is processed differently than its mammalian counterparts. We further show that plants lacking lncCOBRA1 display patterns of delayed germination and are overall smaller than wild-type plants. Lastly, we identify the proteins that interact with lncCOBRA1 and propose a novel mechanism of lincRNA action in which it may act as a scaffold with the RACK1A protein to regulate germination and development, possibly through a role in ribosome biogenesis.

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.

Two evolutionarily duplicated domains individually and post-transcriptionally control SWEET expression for phloem transport

Cell type-specific gene expression is critical for the specialized functions within multicellular organisms. In Arabidopsis, SWEET11 and SWEET12 sugar transporters are specifically expressed in phloem parenchyma (PP) cells and are responsible for sucrose efflux from the PP, the first step of a two-step apoplasmic phloem-loading strategy that initiates the long-distance transport of sugar from leaves to nonphotosynthetic sink tissues. However, we know nothing about what determines the PP cell-specific expression of these SWEETs. Sequence deletions, histochemical β-glucuronidase (GUS) analysis, cross-sectioning, live-cell imaging, and evolutionary analysis were used to elucidate domains responsible for PP specificity, while a Förster resonance energy transfer (FRET) sensor-based transport assay was used to determine whether substrate specificity coevolved with PP specificity. We identified two domains in the Arabidopsis SWEET11 coding sequence that, along with its promoter (including 5' UTR), regulate PP-specific expression at the post-transcriptional level, probably involving RNA-binding proteins. This mechanism is conserved among vascular plants but independent of transport substrate specificity. We conclude that two evolutionarily duplicated coding sequence domains are essential and individually sufficient for PP-specific expression of SWEET11. We also provide a crucial experimental tool to study PP physiology and development.

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.

Length variation in short tandem repeats affects gene expression in natural populations of Arabidopsis thaliana

The genetic basis for the fine-tuned regulation of gene expression is complex and ultimately influences the phenotype and thus the local adaptation of natural populations. Short tandem repeats (STRs) consisting of repetitive DNA motifs have been shown to regulate gene expression. STRs are variable in length within a population and serve as a heritable, but semi-reversible, reservoir of standing genetic variation. For sessile organisms, such as plants, STRs could be of major importance in fine-tuning gene expression as a response to a shifting local environment. Here, we used a transcriptome dataset from natural accessions of Arabidopsis thaliana to investigate population-wide gene expression patterns in light of genome-wide STR variation. We empirically modeled gene expression as a response to the STR length within and around the gene and demonstrated that an association between gene expression and STR length variation is unequivocally present in the sampled population. To support our model, we explored the promoter activity in a transcriptional regulator involved in root hair formation and provided experimentally determined causality between coding sequence length variation and promoter activity. Our results support a general link between gene expression variation and STR length variation in A. thaliana.

Dynamic changes in RNA-protein interactions and RNA secondary structure in mammalian erythropoiesis

Two features of eukaryotic RNA molecules that regulate their post-transcriptional fates are RNA secondary structure and RNA-binding protein (RBP) interaction sites. However, a comprehensive global overview of the dynamic nature of these sequence features during erythropoiesis has never been obtained. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approach to reveal the global landscape of RNA secondary structure and RBP-RNA interaction sites and the dynamics of these features during this important developmental process. We identify dynamic patterns of RNA secondary structure and RBP binding throughout the process and determine a set of corresponding protein-bound sequence motifs along with their dynamic structural and RBP-binding contexts. Finally, using these dynamically bound sequences, we identify a number of RBPs that have known and putative key functions in post-transcriptional regulation during mammalian erythropoiesis. In total, this global analysis reveals new post-transcriptional regulators of mammalian blood cell development.

RNA structure probing uncovers RNA structure-dependent biological functions

RNA molecules fold into complex structures that enable their diverse functions in cells. Recent revolutionary innovations in transcriptome-wide RNA structural probing of living cells have ushered in a new era in understanding RNA functions. Here, we summarize the latest technological advances for probing RNA secondary structures and discuss striking discoveries that have linked RNA regulation and biological processes through interrogation of RNA structures. In particular, we highlight how different long noncoding RNAs form into distinct secondary structures that determine their modes of interactions with protein partners to realize their unique functions. These dynamic structures mediate RNA regulatory functions through altering interactions with proteins and other RNAs. We also outline current methodological hurdles and speculate about future directions for development of the next generation of RNA structure-probing technologies of higher sensitivity and resolution, which could then be applied in increasingly physiologically relevant studies.

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.

Transcription Regulation Through Nascent RNA Folding

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

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

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

BACKGROUND

mRNA processing is critical for gene expression. A challenge in regulating mRNA processing is how to recognize the actual mRNA processing sites, such as splice and polyadenylation sites, when the sequence content is insufficient for this purpose. Previous studies suggested that RNA structure affects mRNA processing. However, the regulatory role of RNA structure in mRNA processing remains unclear.

RESULTS

Here, we perform in vivo selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemical profiling on Arabidopsis and generate the in vivo nuclear RNA structure landscape. We find that nuclear mRNAs fold differently from cytosolic mRNAs across translation start and stop sites. Notably, we discover a two-nucleotide single-stranded RNA structure feature upstream of 5' splice sites that is strongly associated with splicing and the selection of alternative 5' splice sites. The regulatory role of this RNA structure feature is further confirmed by experimental validation. Moreover, we find the single-strandedness of branch sites is also associated with 3' splice site recognition. We also identify an RNA structure feature comprising two close-by single-stranded regions that is specifically associated with both polyadenylation and alternative polyadenylation events.

CONCLUSIONS

We successfully identify pre-mRNA structure features associated with splicing and polyadenylation at whole-genome scale and validate an RNA structure feature which can regulate splicing. Our study unveils a new RNA structure regulatory mechanism for mRNA processing.

Messenger RNA 5' NAD+ Capping Is a Dynamic Regulatory Epitranscriptome Mark That Is Required for Proper Response to Abscisic Acid in Arabidopsis

Although eukaryotic messenger RNAs (mRNAs) normally possess a 5' end N⁷-methyl guanosine (m⁷G) cap, a non-canonical 5' nicotinamide adenine dinucleotide (NAD⁺) cap can tag certain transcripts for degradation mediated by the NAD⁺ decapping enzyme DXO1. Despite this importance, whether NAD⁺ capping dynamically responds to specific stimuli to regulate eukaryotic transcriptomes remains unknown. Here, we reveal a link between NAD⁺ capping and tissue- and hormone response-specific mRNA stability. In the absence of DXO1 function, transcripts displaying a high proportion of NAD⁺ capping are instead processed into RNA-dependent RNA polymerase 6-dependent small RNAs, resulting in their continued turnover likely to free the NAD⁺ molecules. Additionally, the NAD⁺-capped transcriptome is significantly remodeled in response to the essential plant hormone abscisic acid in a mechanism that is primarily independent of DXO1. Overall, our findings reveal a previously uncharacterized and essential role of NAD⁺ capping in dynamically regulating transcript stability during specific physiological responses.

Processing of coding and non-coding RNAs in plant development and environmental responses

Precursor RNAs undergo extensive processing to become mature RNAs. RNA transcripts are subjected to 5′ capping, 3′-end processing, splicing, and modification; they also form dynamic secondary structures during co-transcriptional and post-transcriptional processing. Like coding RNAs, non-coding RNAs (ncRNAs) undergo extensive processing. For example, secondary small interfering RNA (siRNA) transcripts undergo RNA processing, followed by further cleavage to become mature siRNAs. Transcriptome studies have revealed roles for co-transcriptional and post-transcriptional RNA processing in the regulation of gene expression and the coordination of plant development and plant–environment interactions. In this review, we present the latest progress on RNA processing in gene expression and discuss phased siRNAs (phasiRNAs), a kind of germ cell-specific secondary small RNA (sRNA), focusing on their functions in plant development and environmental responses.

The increasing diversity and complexity of the RNA-binding protein repertoire in plants

Post-transcriptional regulation has far-reaching implications on the fate of RNAs. It is gaining increasing momentum as a critical component in adjusting global cellular transcript levels during development and in response to environmental stresses. In this process, RNA-binding proteins (RBPs) are indispensable chaperones that naturally bind RNA via one or multiple globular RNA-binding domains (RBDs) changing the function or fate of the bound RNAs. Despite the technical challenges faced in plants in large-scale studies, several hundreds of these RBPs have been discovered and elucidated globally over the past few years. Recent discoveries have more than doubled the number of proteins implicated in RNA interaction, including identification of RBPs lacking classical RBDs. This review will discuss these new emerging classes of RBPs, focusing on the current state of the RBP repertoire in Arabidopsis thaliana, including the diverse functional roles derived from quantitative studies implicating RBPs in abiotic stress responses. Notably, this review highlights that 836 RBPs are enriched as Arabidopsis RBPs while 1865 can be classified as candidate RBPs. The review will also outline outstanding areas within this field that require addressing to advance our understanding and potential biotechnological applications of RBPs.

mRNA adenosine methylase (MTA) deposits m6A on pri-miRNAs to modulate miRNA biogenesis in Arabidopsis thaliana

In Arabidopsis thaliana, the METTL3 homolog, mRNA adenosine methylase (MTA) introduces N⁶-methyladenosine (m⁶A) into various coding and noncoding RNAs of the plant transcriptome. Here, we show that an MTA-deficient mutant (mta) has decreased levels of microRNAs (miRNAs) but accumulates primary miRNA transcripts (pri-miRNAs). Moreover, pri-miRNAs are methylated by MTA, and RNA structure probing analysis reveals a decrease in secondary structure within stem-loop regions of these transcripts in mta mutant plants. We demonstrate interaction between MTA and both RNA Polymerase II and TOUGH (TGH), a plant protein needed for early steps of miRNA biogenesis. Both MTA and TGH are necessary for efficient colocalization of the Microprocessor components Dicer-like 1 (DCL1) and Hyponastic Leaves 1 (HYL1) with RNA Polymerase II. We propose that secondary structure of miRNA precursors induced by their MTA-dependent m⁶A methylation status, together with direct interactions between MTA and TGH, influence the recruitment of Microprocessor to plant pri-miRNAs. Therefore, the lack of MTA in mta mutant plants disturbs pri-miRNA processing and leads to the decrease in miRNA accumulation. Furthermore, our findings reveal that reduced miR393b levels likely contributes to the impaired auxin response phenotypes of mta mutant plants.

Mechanisms and Regulation of RNA Condensation in RNP Granule Formation

Ribonucleoprotein (RNP) granules are RNA-protein assemblies that are involved in multiple aspects of RNA metabolism and are linked to memory, development, and disease. Some RNP granules form, in part, through the formation of intermolecular RNA-RNA interactions. In vitro, such trans RNA condensation occurs readily, suggesting that cells require mechanisms to modulate RNA-based condensation. We assess the mechanisms of RNA condensation and how cells modulate this phenomenon. We propose that cells control RNA condensation through ATP-dependent processes, static RNA buffering, and dynamic post-translational mechanisms. Moreover, perturbations in these mechanisms can be involved in disease. This reveals multiple cellular mechanisms of kinetic and thermodynamic control that maintain the proper distribution of RNA molecules between dispersed and condensed forms.

mRNA adenosine methylase (MTA) deposits m6A on pri-miRNAs to modulate miRNA biogenesis in Arabidopsis thaliana

Recently N⁶-methyladenosine (m⁶A) methylation has emerged as a biological process with significant impact on cellular functions. However, almost all the research regarding m⁶A methylation has been based on mRNAs. In our research, we focus on how m⁶A methylation affects microRNA (miRNA) biogenesis in Arabidopsis. In brief, we show that m⁶A methylation is necessary to maintain proper levels of mature miRNAs as well as their precursors. m⁶A mark affects pri-miRNA secondary structures and affects the recruitment of the Microprocessor to pri-miRNAs. We also demonstrate the interactions of MTA (m⁶A writer) with other proteins involved in miRNA biogenesis, namely RNA Polymerase II and TOUGH. Our study provides evidence of the role played by m⁶A in plant miRNA biogenesis.

In Arabidopsis thaliana, the METTL3 homolog, mRNA adenosine methylase (MTA) introduces N⁶-methyladenosine (m⁶A) into various coding and noncoding RNAs of the plant transcriptome. Here, we show that an MTA-deficient mutant (mta) has decreased levels of microRNAs (miRNAs) but accumulates primary miRNA transcripts (pri-miRNAs). Moreover, pri-miRNAs are methylated by MTA, and RNA structure probing analysis reveals a decrease in secondary structure within stem–loop regions of these transcripts in mta mutant plants. We demonstrate interaction between MTA and both RNA Polymerase II and TOUGH (TGH), a plant protein needed for early steps of miRNA biogenesis. Both MTA and TGH are necessary for efficient colocalization of the Microprocessor components Dicer-like 1 (DCL1) and Hyponastic Leaves 1 (HYL1) with RNA Polymerase II. We propose that secondary structure of miRNA precursors induced by their MTA-dependent m⁶A methylation status, together with direct interactions between MTA and TGH, influence the recruitment of Microprocessor to plant pri-miRNAs. Therefore, the lack of MTA in mta mutant plants disturbs pri-miRNA processing and leads to the decrease in miRNA accumulation. Furthermore, our findings reveal that reduced miR393b levels likely contributes to the impaired auxin response phenotypes of mta mutant plants.

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.

The landscape of eukaryotic mRNPs

The proper regulation of mRNA processing, localization, translation, and degradation occurs on mRNPs. However, the global principles of mRNP organization are poorly understood. We utilize the limited, but existing, information available to present a speculative synthesis of mRNP organization with the following key points. First, mRNPs form a compacted structure due to the inherent folding of RNA. Second, the ribosome is the principal mechanism by which mRNA regions are partially decompacted. Third, mRNPs are 50%-80% protein by weight, consistent with proteins modulating mRNP organization, but also suggesting the majority of mRNA sequences are not directly interacting with RNA-binding proteins. Finally, the ratio of mRNA-binding proteins to mRNAs is higher in the nucleus to allow effective RNA processing and limit the potential for nuclear RNA based aggregation. This synthesis of mRNP understanding provides a model for mRNP biogenesis, structure, and regulation with multiple implications.

Read Mapping and Transcript Assembly: A Scalable and High-Throughput Workflow for the Processing and Analysis of Ribonucleic Acid Sequencing Data

Next-generation RNA-sequencing is an incredibly powerful means of generating a snapshot of the transcriptomic state within a cell, tissue, or whole organism. As the questions addressed by RNA-sequencing (RNA-seq) become both more complex and greater in number, there is a need to simplify RNA-seq processing workflows, make them more efficient and interoperable, and capable of handling both large and small datasets. This is especially important for researchers who need to process hundreds to tens of thousands of RNA-seq datasets. To address these needs, we have developed a scalable, user-friendly, and easily deployable analysis suite called RMTA (Read Mapping, Transcript Assembly). RMTA can easily process thousands of RNA-seq datasets with features that include automated read quality analysis, filters for lowly expressed transcripts, and read counting for differential expression analysis. RMTA is containerized using Docker for easy deployment within any compute environment [cloud, local, or high-performance computing (HPC)] and is available as two apps in CyVerse's Discovery Environment, one for normal use and one specifically designed for introducing undergraduates and high school to RNA-seq analysis. For extremely large datasets (tens of thousands of FASTq files) we developed a high-throughput, scalable, and parallelized version of RMTA optimized for launching on the Open Science Grid (OSG) from within the Discovery Environment. OSG-RMTA allows users to utilize the Discovery Environment for data management, parallelization, and submitting jobs to OSG, and finally, employ the OSG for distributed, high throughput computing. Alternatively, OSG-RMTA can be run directly on the OSG through the command line. RMTA is designed to be useful for data scientists, of any skill level, interested in rapidly and reproducibly analyzing their large RNA-seq data sets.

Using RNA Affinity Purification Followed by Mass Spectrometry to Identify RNA-Binding Proteins (RBPs)

RNA-binding proteins (RBPs) perform key functions in posttranscriptional regulation, adding complexity to the RNA life cycle. RNA interactome capture techniques have been applied to various organisms of interest and detected hundreds of RBPs, some with uncharacterized functions. However, even in many well-studied organisms, the primary sequence motif for most RBPs remains unknown. Here, we describe a 3-day protocol where users couple an RNA sequence of interest that is known to be bound by an RBP(s) with agarose beads, incubate the now tagged RNA sequence with protein lysate, and then pull down the proteins bound to the RNA. Subsequent mass spectrometry allows users to profile the RNA sequence-interacting proteome and pick out any enriched proteins as RBPs of interest. This protocol allows researchers to match sequences to their RBPs and even often identify novel RBPs or new functions for known RBPs.

Emerging roles of RNA-binding proteins in plant development

RNA-binding proteins (RBPs) influence the fate of target RNAs via direct interactions. During transcription, RBPs and interacting partners are recruited to and modify transcripts, after which they may also participate in critical steps to generate functional RNA. RBP-RNA interactions govern post-transcriptional processing of RNA, consequently regulating gene expression in a spatio-temporal manner. In plants, an increasing number of proteins have been classified as RBPs, many of which have been shown to function as key players in diverse developmental processes. However, a comprehensive understanding of how RBPs function, which RNAs are targeted, and where RBP-RNA interactions occur within plant cells is lacking. Here, we discuss recent findings in the field and newly defined roles for RBPs in plant growth and development. We also describe the mechanistic effects of RBPs on target RNA metabolism and translation.

Genome-Wide Identification of Long Noncoding RNAs and Their Responses to Salt Stress in Two Closely Related Poplars

Long noncoding RNAs (lncRNAs) are involved in various biological regulatory processes, but their roles in plants resistance to salt stress remain largely unknown. To systematically explore the characteristics of lncRNAs and their roles in plant salt responses, we conducted strand-specific RNA-sequencing of four tissue types with salt treatments in two closely related poplars (Populus euphratica and Populus alba var. pyramidalis), and a total of 10,646 and 10,531 lncRNAs were identified, respectively. These lncRNAs showed significantly lower values in terms of length, expression, and expression correction than with mRNA. We further found that about 40% and 60% of these identified lncRNAs responded to salt stress with tissue-specific expression patterns across the two poplars. Furthermore, lncRNAs showed weak evolutionary conservation in sequences and exhibited diverse regulatory styles; in particular, tissue- and species-specific responses to salt stress varied greatly in two poplars, for example, 322 lncRNAs were found highly expressed in P. euphratica but not in P. alba var. pyramidalis and 3,425 lncRNAs were identified to be species-specific in P. euphratica in response to salt stress. Moreover, tissue-specific expression of lncRNAs in two poplars were identified with predicted target genes included Aux/IAA, NAC, MYB, involved in regulating plant growth and the plant stress response. Taken together, the systematic analysis of lncRNAs between sister species enhances our understanding of the characteristics of lncRNAs and their roles in plant growth and salt response.

Transcription Is Just the Beginning of Gene Expression Regulation: The Functional Significance of RNA-Binding Proteins to Post-transcriptional Processes in Plants

Plants have developed sophisticated mechanisms to compensate and respond to ever-changing environmental conditions. Research focus in this area has recently shifted towards understanding the post-transcriptional mechanisms that contribute to RNA transcript maturation, abundance and function as key regulatory steps in allowing plants to properly react and adapt to these never-ending shifts in their environments. At the center of these regulatory mechanisms are RNA-binding proteins (RBPs), the functional mediators of all post-transcriptional processes. In plants, RBPs are becoming increasingly appreciated as the critical modulators of core cellular processes during development and in response to environmental stimuli. With the majority of research on RBPs and their functions historically in prokaryotic and mammalian systems, it has more recently been unveiled that plants have expanded families of conserved and novel RBPs compared with their eukaryotic counterparts. To better understand the scope of RBPs in plants, we present past and current literature detailing specific roles of RBPs during stress response, development and other fundamental transition periods. In this review, we highlight examples of complex regulation coordinated by RBPs with a focus on the diverse mechanisms of plant RBPs and the unique processes they regulate. Additionally, we discuss the importance for additional research into understanding global interactions of RBPs on a systems and network-scale, with genome mining and annotation providing valuable insight for potential uses in improving crop plants in order to maintain high-level production in this era of global climate change.

Beyond CLIP: advances and opportunities to measure RBP-RNA and RNA-RNA interactions

RNA is an essential player in almost all biological processes, and has an ever-growing number of roles in regulating cellular growth and organization. RNA functions extend far beyond just coding for proteins and RNA has been shown to function in signaling events, chromatin organization and transcriptional regulation. Dissecting how the complex network of RNA-binding proteins (RBPs) and regulatory RNAs interact with their substrates within the cell is a real, but exciting, challenge for the RNA community. Investigating these biological questions has fueled the development of new quantitative technologies to measure how RNA and RBPs interact both locally and on a global scale. In this review, we provide an assessment of available approaches to enable researchers to select the protocol most applicable for their experimental question.

RNA structure maps across mammalian cellular compartments

RNA structure is intimately connected to each step of gene expression. Recent advances have enabled transcriptome-wide maps of RNA secondary structure, termed RNA structuromes. However, previous whole-cell analyses lacked the resolution to unravel the landscape and also the regulatory mechanisms of RNA structural changes across subcellular compartments. Here we reveal the RNA structuromes in three compartments — chromatin, nucleoplasm and cytoplasm in human and mouse cells. The cytotopic structuromes substantially expand RNA structural information, and enable detailed investigation of the central role of RNA structure in linking transcription, translation, and RNA decay. We develop a resource to visualize the interplay of RNA-protein interactions, RNA modifications, and RNA structure, and predict both direct and indirect reader proteins of RNA modifications. We also validate a novel role of the RNA binding protein LIN28A as an N6-methyladenosine modification “anti-reader”. Our results highlight the dynamic nature of RNA structures and its functional significance in gene regulation.

Computational analysis of alternative splicing in plant genomes

Computational analyses play crucial roles in characterizing splicing isoforms in plant genomes. In this review, we provide a survey of computational tools used in recently published, genome-scale splicing analyses in plants. We summarize the commonly used software and pipelines for read mapping, isoform reconstruction, isoform quantification, and differential expression analysis. We also discuss methods for analyzing long reads and the strategies to combine long and short reads in identifying splicing isoforms. We review several tools for characterizing local splicing events, splicing graphs, coding potential, and visualizing splicing isoforms. We further discuss the procedures for identifying conserved splicing isoforms across plant species. Finally, we discuss the outlook of integrating other genomic data with splicing analyses to identify regulatory mechanisms of AS on genome-wide scale.

Using Protein Interaction Profile Sequencing (PIP-seq) to Identify RNA Secondary Structure and RNA-Protein Interaction Sites of Long Noncoding RNAs in Plants

From the moment of transcription, RNA molecules are continuously bound by RNA-binding proteins (RBPs). While the majority of research has focused on how these RBPs regulate posttranscriptional gene regulation of messenger RNAs (mRNAs), the majority of cellular RNAs do not code for proteins, such as ribosomal RNAs, transfer RNAs, and microRNAs. Since these RNAs do not code for protein, their function is mainly determined by their interactions with RBPs as well as their intramolecular base pairing, or RNA secondary structure. One class of noncoding RNAs termed long noncoding RNAs (lncRNAs) have recently become the subject of intense research interest. To study the function of lncRNAs in eukaryotic cells, it is important to examine both their interactions with RBPs as well as their RNA secondary structure. Protein interaction profile sequencing (PIP-seq) is a genome-wide method that uses structure-specific ribonucleases (RNases) to identify regions of double-stranded and single-stranded RNA as well as regions that are protected from these RNases, which represent sites of RBP binding. This method is a very powerful way to examine RNA-protein interactions and RNA secondary structure of all lncRNAs expressed in cells and tissues and can reveal cell-type-specific or tissue-specific patterns of these RNA features. Here, we give a step-by-step account of performing this technique for comprehensively analyzing RNA-protein interactions and RNA secondary structure in plant transcriptomes.

An Overview of Methodologies in Studying lncRNAs in the High-Throughput Era: When Acronyms ATTACK!

The discovery of pervasive transcription in eukaryotic genomes provided one of many surprising (and perhaps most surprising) findings of the genomic era and led to the uncovering of a large number of previously unstudied transcriptional events. This pervasive transcription leads to the production of large numbers of noncoding RNAs (ncRNAs) and thus opened the window to study these diverse, abundant transcripts of unclear relevance and unknown function. Since that discovery, recent advances in high-throughput sequencing technologies have identified a large collection of ncRNAs, from microRNAs to long noncoding RNAs (lncRNAs). Subsequent discoveries have shown that many lncRNAs play important roles in various eukaryotic processes; these discoveries have profoundly altered our understanding of the regulation of eukaryotic gene expression. Although the identification of ncRNAs has become a standard experimental approach, the functional characterization of these diverse ncRNAs remains a major challenge. In this chapter, we highlight recent progress in the methods to identify lncRNAs and the techniques to study the molecular function of these lncRNAs and the application of these techniques to the study of plant lncRNAs.

Genome-wide RNA structurome reprogramming by acute heat shock globally regulates mRNA abundance

The heat shock response is crucial for organism survival in natural environments. RNA structure is known to influence numerous processes related to gene expression, but there have been few studies on the global RNA structurome as it prevails in vivo. Moreover, how heat shock rapidly affects RNA structure genome-wide in living systems remains unknown. We report here in vivo heat-regulated RNA structuromes. We applied Structure-seq chemical [dimethyl sulfate (DMS)] structure probing to rice (Oryza sativa L.) seedlings with and without 10 min of 42 °C heat shock and obtained structural data on >14,000 mRNAs. We show that RNA secondary structure broadly regulates gene expression in response to heat shock in this essential crop species. Our results indicate significant heat-induced elevation of DMS reactivity in the global transcriptome, revealing RNA unfolding over this biological temperature range. Our parallel Ribo-seq analysis provides no evidence for a correlation between RNA unfolding and heat-induced changes in translation, in contrast to the paradigm established in prokaryotes, wherein melting of RNA thermometers promotes translation. Instead, we find that heat-induced DMS reactivity increases correlate with significant decreases in transcript abundance, as quantified from an RNA-seq time course, indicating that mRNA unfolding promotes transcript degradation. The mechanistic basis for this outcome appears to be mRNA unfolding at both 5' and 3'-UTRs that facilitates access to the RNA degradation machinery. Our results thus reveal unexpected paradigms governing RNA structural changes and the eukaryotic RNA life cycle.

The nucleotides they are a-changin': function of RNA binding proteins in post-transcriptional messenger RNA editing and modification in Arabidopsis

During and after transcription, the fate of an RNA molecule is almost entirely directed by the cohorts of interacting RNA-binding proteins (RBPs). RBPs regulate all stages of the life cycle of a messenger RNA (mRNA) molecule, including splicing, polyadenylation, transport out of the nucleus, RNA stability, and translation. In addition to these functions, RBPs can function to modify or edit the sequences encoded by the RNA. While the sequence for each transcript is determined in the genome, by the time an RNA reaches its final fate, the sequence may have been edited, where one nucleotide is converted to another, or modified, where a chemical group, or sometimes others moieties, are covalently linked to a nucleotide base. These changes to the RNA sequence have major consequences on the function of the RNA. Additionally, variation in the levels of the RBPs that perform the editing or modification can drastically affect the fitness of an organism. Here, we review RBPs that are known to edit or modify RNA ribonucleotides, focusing on the RNA editing ability of the pentatricopeptide repeat (PPR) proteins and the RBPs that modify adenosine to N⁶- methyladenosine.

New Era of Studying RNA Secondary Structure and Its Influence on Gene Regulation in Plants

The dynamic structure of RNA plays a central role in post-transcriptional regulation of gene expression such as RNA maturation, degradation, and translation. With the rise of next-generation sequencing, the study of RNA structure has been transformed from in vitro low-throughput RNA structure probing methods to in vivo high-throughput RNA structure profiling. The development of these methods enables incremental studies on the function of RNA structure to be performed, revealing new insights of novel regulatory mechanisms of RNA structure in plants. Genome-wide scale RNA structure profiling allows us to investigate general RNA structural features over 10s of 1000s of mRNAs and to compare RNA structuromes between plant species. Here, we provide a comprehensive and up-to-date overview of: (i) RNA structure probing methods; (ii) the biological functions of RNA structure; (iii) genome-wide RNA structural features corresponding to their regulatory mechanisms; and (iv) RNA structurome evolution in plants.

Nuclear Transcriptomes at High Resolution Using Retooled INTACT

Isolated nuclei provide access to early steps in gene regulation involving chromatin as well as transcript production and processing. Here, we describe transfer of the isolation of nuclei from tagged specific cell types (INTACT) to the monocot rice (Oryza sativa L.). The purification of biotinylated nuclei was redesigned by replacing the outer nuclear-envelope-targeting domain of the nuclear tagging fusion (NTF) protein with an outer nuclear-envelope-anchored domain. This modified NTF was combined with codon-optimized Escherichia coli BirA in a single T-DNA construct. We also developed inexpensive methods for INTACT, T-DNA insertion mapping, and profiling of the complete nuclear transcriptome, including a ribosomal RNA degradation procedure that minimizes pre-ribosomal RNA (pre-rRNA) transcripts. A high-resolution comparison of nuclear and steady-state poly(A)+ transcript populations of seedling root tips confirmed the capture of pre-messenger RNA (pre-mRNA) and exposed distinctions in diversity and abundance of the nuclear and total transcriptomes. This retooled INTACT can enable high-resolution monitoring of the nuclear transcriptome and chromatin in specific cell types of rice and other species.

Systems Approaches to Map In Vivo RNA–Protein Interactions in Arabidopsis thaliana

Proteins that specifically interact with mRNAs orchestrate mRNA processing steps all the way from transcription to decay. Thus, these RNA-binding proteins represent an important control mechanism to double check which proportion of nascent pre-mRNAs is ultimately available for translation into distinct proteins. Here, we discuss recent progress to obtain a systems-level understanding of in vivo RNA–protein interactions in the reference plant Arabidopsis thaliana using protein-centric and RNA-centric methods as well as combined protein binding site and structure probing.

A Guide to the Chloroplast Transcriptome Analysis Using RNA-Seq

Since its first use in plants in 2007, high-throughput RNA sequencing (RNA-Seq) has generated a vast amount of data for both model and nonmodel species. Organellar transcriptomes, however, are virtually always overlooked at the data analysis step. We therefore developed ChloroSeq, a bioinformatic pipeline aimed at facilitating the systematic analysis of chloroplast RNA metabolism, and we provide here a step-by-step user's manual. Following the alignment of quality-controlled data to the genome of interest, ChloroSeq measures genome expression level along with splicing and RNA editing efficiencies. When used in combination with the Tuxedo suite (TopHat and Cufflinks), ChloroSeq allows the simultaneous analysis of organellar and nuclear transcriptomes, opening the way to a better understanding of nucleus-organelle cross talk. We also describe the use of R commands to produce publication-quality figures based on ChloroSeq outputs. The effectiveness of the pipeline is illustrated through analysis of an RNA-Seq dataset covering the transition from growth to maturation to senescence of Arabidopsis thaliana leaves.