Genomic era analyses of RNA secondary structure and RNA-binding proteins reveal their significance to post-transcriptional regulation in plants

Plant RNA-binding proteins: Phase separation dynamics and functional mechanisms underlying plant development and stress responses

RNA-binding proteins (RBPs) accompany RNA from synthesis to decay, mediating every aspect of RNA metabolism and impacting diverse cellular and developmental processes in eukaryotes. Many RBPs undergo phase separation along with their bound RNA to form and function in dynamic membraneless biomolecular condensates for spatiotemporal coordination or regulation of RNA metabolism. Increasing evidence suggests that phase-separating RBPs with RNA-binding domains and intrinsically disordered regions play important roles in plant development and stress adaptation. Here, we summarize the current knowledge about how dynamic partitioning of RBPs into condensates controls plant development and enables sensing of experimental changes to confer growth plasticity under stress conditions, with a focus on the dynamics and functional mechanisms of RBP-rich nuclear condensates and cytoplasmic granules in mediating RNA metabolism. We also discuss roles of multiple factors, such as environmental signals, protein modifications, and N6-methyladenosine RNA methylation, in modulating the phase separation behaviors of RBPs, and highlight the prospects and challenges for future research on phase-separating RBPs in crops.

A transient in planta editing assay identifies specific binding of the splicing regulator PTB as a prerequisite for cassette exon inclusion

The dynamic interaction of RNA-binding proteins (RBPs) with their target RNAs contributes to the diversity of ribonucleoprotein (RNP) complexes that are involved in a myriad of biological processes. Identifying the RNP components at high resolution and defining their interactions are key to understanding their regulation and function. Expressing fusions between an RBP of interest and an RNA editing enzyme can result in nucleobase changes in target RNAs, representing a recent addition to experimental approaches for profiling RBP/RNA interactions. Here, we have used the MS2 protein/RNA interaction to test four RNA editing proteins for their suitability to detect target RNAs of RBPs in planta. We have established a transient test system for fast and simple quantification of editing events and identified the hyperactive version of the catalytic domain of an adenosine deaminase (hADARcd) as the most suitable editing enzyme. Examining fusions between homologs of polypyrimidine tract binding proteins (PTBs) from Arabidopsis thaliana and hADARcd allowed determining target RNAs with high sensitivity and specificity. Moreover, almost complete editing of a splicing intermediate provided insight into the order of splicing reactions and PTB dependency of this particular splicing event. Addition of sequences for nuclear localisation of the fusion protein increased the editing efficiency, highlighting this approach's potential to identify RBP targets in a compartment-specific manner. Our studies have established the editing-based analysis of interactions between RBPs and their RNA targets in a fast and straightforward assay, offering a new system to study the intricate composition and functions of plant RNPs in vivo.

Plant phase extraction: A method for enhanced discovery of the RNA-binding proteome and its dynamics in plants

RNA-binding proteins (RBPs) play critical roles in posttranscriptional gene regulation. Current methods of systematically profiling RBPs in plants have been predominantly limited to proteins interacting with polyadenylated (poly(A)) RNAs. We developed a method called plant phase extraction (PPE), which yielded a highly comprehensive RNA-binding proteome (RBPome), uncovering 2,517 RBPs from Arabidopsis (Arabidopsis thaliana) leaf and root samples with a highly diverse array of RNA-binding domains. We identified traditional RBPs that participate in various aspects of RNA metabolism and a plethora of nonclassical proteins moonlighting as RBPs. We uncovered constitutive and tissue-specific RBPs essential for normal development and, more importantly, revealed RBPs crucial for salinity stress responses from a RBP-RNA dynamics perspective. Remarkably, 40% of the RBPs are non-poly(A) RBPs that were not previously annotated as RBPs, signifying the advantage of PPE in unbiasedly retrieving RBPs. We propose that intrinsically disordered regions contribute to their nonclassical binding and provide evidence that enzymatic domains from metabolic enzymes have additional roles in RNA binding. Taken together, our findings demonstrate that PPE is an impactful approach for identifying RBPs from complex plant tissues and pave the way for investigating RBP functions under different physiological and stress conditions at the posttranscriptional level.

Investigating the NRAS 5' UTR as a target for small molecules

Neuroblastoma RAS (NRAS) is an oncogene that is deregulated and highly mutated in cancers including melanomas and acute myeloid leukemias. The 5' untranslated region (UTR) (5' UTR) of the NRAS mRNA contains a G-quadruplex (G4) that regulates translation. Here we report a novel class of small molecule that binds to the G4 structure located in the 5' UTR of the NRAS mRNA. We used a small molecule microarray screen to identify molecules that selectively bind to the NRAS-G4 with submicromolar affinity. One compound inhibits the translation of NRAS in vitro but showed only moderate effects on the NRAS levels in cellulo. Rapid Amplification of cDNA Ends and RT-PCR analysis revealed that the predominant NRAS transcript does not possess the G4 structure. Thus, although NRAS transcripts lack a G4 in many cell lines the concept of targeting folded regions within 5' UTRs to control translation remains a highly attractive strategy.

Toward a systems view on RNA-binding proteins and associated RNAs in plants: Guilt by association

RNA-binding proteins (RBPs) have a broad impact on most biochemical, physiological, and developmental processes in a plant's life. RBPs engage in an on-off relationship with their RNA partners, accompanying virtually every stage in RNA processing and function. While the function of a plethora of RBPs in plant development and stress responses has been described, we are lacking a systems-level understanding of components in RNA-based regulation. Novel techniques have substantially enlarged the compendium of proteins with experimental evidence for binding to RNAs in the cell, the RNA-binding proteome. Furthermore, ribonomics methods have been adapted for use in plants to profile the in vivo binding repertoire of RBPs genome-wide. Here, we discuss how recent technological achievements have provided novel insights into the mode of action of plant RBPs at a genome-wide scale. Furthermore, we touch upon two emerging topics, the connection of RBPs to phase separation in the cell and to extracellular RNAs. Finally, we define open questions to be addressed to move toward an integrated understanding of RBP function.

BPL3 binds the long non-coding RNA nalncFL7 to suppress FORKED-LIKE7 and modulate HAI1-mediated MPK3/6 dephosphorylation in plant immunity

RNA-binding proteins (RBPs) participate in a diverse set of biological processes in plants, but their functions and underlying mechanisms in plant-pathogen interactions are largely unknown. We previously showed that Arabidopsis thaliana BPA1-LIKE PROTEIN3 (BPL3) belongs to a conserved plant RBP family and negatively regulates reactive oxygen species (ROS) accumulation and cell death under biotic stress. In this study, we demonstrate that BPL3 suppresses FORKED-LIKE7 (FL7) transcript accumulation and raises levels of the cis-natural antisense long non-coding RNA (lncRNA) of FL7 (nalncFL7). FL7 positively regulated plant immunity to Phytophthora capsici while nalncFL7 negatively regulated resistance. We also showed that BPL3 directly binds to and stabilizes nalncFL7. Moreover, nalncFL7 suppressed accumulation of FL7 transcripts. Furthermore, FL7 interacted with HIGHLY ABA-INDUCED PP2C1 (HAI1), a type 2C protein phosphatase, and inhibited HAI1 phosphatase activity. By suppressing HAI1 activity, FL7 increased the phosphorylation levels of MITOGEN-ACTIVATED PROTEIN KINASE 3 (MPK3) and MPK6, thus enhancing immunity responses. BPL3 and FL7 are conserved in all plant species tested, but the BPL3-nalncFL7-FL7 cascade was specific to the Brassicaceae. Thus, we identified a conserved BPL3-nalncFL7-FL7 cascade that coordinates plant immunity.

Identification of Pri-miRNA Stem-Loop Interacting Proteins in Plants Using a Modified Version of the Csy4 CRISPR Endonuclease

Regulation at the RNA level by RNA-binding proteins (RBPs) and microRNAs (miRNAs) is key to coordinating eukaryotic gene expression. In plants, the importance of miRNAs is highlighted by severe developmental defects in mutants impaired in miRNA biogenesis. MiRNAs are processed from long primary-microRNAs (pri-miRNAs) with internal stem-loop structures by endonucleolytic cleavage. The highly structured stem-loops constitute the basis for the extensive regulation of miRNA biogenesis through interaction with RBPs. However, trans-acting regulators of the biogenesis of specific miRNAs are largely unknown in plants. Therefore, we exploit an RNA-centric approach based on modified versions of the conditional CRISPR nuclease Csy4* to pull down interactors of the Arabidopsis pri-miR398b stem-loop (pri-miR398b-SL) in vitro. We designed three epitope-tagged versions of the inactive Csy4* for the immobilization of the protein together with the pri-miR398b-SL bait on high affinity matrices. After incubation with nucleoplasmic extracts from Arabidopsis and extensive washing, pri-miR398b-SL, along with its specifically bound proteins, were released by re-activating the cleavage activity of the Csy4* upon the addition of imidazole. Co-purified proteins were identified via quantitative mass spectrometry and data sets were compared. In total, we identified more than 400 different proteins, of which 180 are co-purified in at least two out of three independent Csy4*-based RNA pulldowns. Among those, the glycine-rich RNA-binding protein AtRZ-1a was identified in all pulldowns. To analyze the role of AtRZ-1a in miRNA biogenesis, we determined the miR398 expression level in the atrz-1a mutant. Indeed, the absence of AtRZ-1a caused a decrease in the steady-state level of mature miR398 with a concomitant reduction in pri-miR398b levels. Overall, we show that our modified Csy4*-based RNA pulldown strategy is suitable to identify new trans-acting regulators of miRNA biogenesis and provides new insights into the post-transcriptional regulation of miRNA processing by plant RBPs.

Glycine-rich RNA-binding cofactor RZ1AL is associated with tomato ripening and development

Tomato ripening is a complex and dynamic process coordinated by many regulatory elements, including plant hormones, transcription factors, and numerous ripening-related RNAs and proteins. Although recent studies have shown that some RNA-binding proteins are involved in the regulation of the ripening process, understanding of how RNA-binding proteins affect fruit ripening is still limited. Here, we report the analysis of a glycine-rich RNA-binding protein, RZ1A-Like (RZ1AL), which plays an important role in tomato ripening, especially fruit coloring. To analyze the functions of RZ1AL in fruit development and ripening, we generated knockout cr-rz1al mutant lines via the CRISPR/Cas9 gene-editing system. Knockout of RZ1AL reduced fruit lycopene content and weight in the cr-rz1al mutant plants. RZ1AL encodes a nucleus-localized protein that is associated with Cajal-related bodies. RNA-seq data demonstrated that the expression levels of genes that encode several key enzymes associated with carotenoid biosynthesis and metabolism were notably downregulated in cr-rz1al fruits. Proteomic analysis revealed that the levels of various ribosomal subunit proteins were reduced. This could affect the translation of ripening-related proteins such as ZDS. Collectively, our findings demonstrate that RZ1AL may participate in the regulation of carotenoid biosynthesis and metabolism and affect tomato development and fruit ripening.

Conserved long-range base pairings are associated with pre-mRNA processing of human genes

The ability of nucleic acids to form double-stranded structures is essential for all living systems on Earth. Current knowledge on functional RNA structures is focused on locally-occurring base pairs. However, crosslinking and proximity ligation experiments demonstrated that long-range RNA structures are highly abundant. Here, we present the most complete to-date catalog of conserved complementary regions (PCCRs) in human protein-coding genes. PCCRs tend to occur within introns, suppress intervening exons, and obstruct cryptic and inactive splice sites. Double-stranded structure of PCCRs is supported by decreased icSHAPE nucleotide accessibility, high abundance of RNA editing sites, and frequent occurrence of forked eCLIP peaks. Introns with PCCRs show a distinct splicing pattern in response to RNAPII slowdown suggesting that splicing is widely affected by co-transcriptional RNA folding. The enrichment of 3'-ends within PCCRs raises the intriguing hypothesis that coupling between RNA folding and splicing could mediate co-transcriptional suppression of premature pre-mRNA cleavage and polyadenylation.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.

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.

Comparative analysis of RNA-binding proteomes under Arabidopsis thaliana-Pst DC3000-PAMP interaction by orthogonal organic phase separation

RNA-binding proteins (RBPs) are pivotal participants in post-transcriptional gene regulation. They interact with RNA directly to perform several post-transcriptional RNA regulatory functions or direct metabolic processes. Despite the essential importance, the understanding of plant RBPs is elementary, which derives mainly from other kingdoms via bioinformatic extrapolation or mRNA-binding proteins captured through UV crosslinked method. Recently, orthogonal organic phase separation (OOPS) method for RBP identification has been used in mammals and Escherichia coli. And plentiful RBPs were enriched without molecular tagging or capture of polyadenylated RNA in an unbiased way. In our study, OOPS was conducted on Arabidopsis and 468 RBPs were discovered including 244 putative RBPs. There were 17 peroxidases in 232 RBPs with enzymatic activities. In addition, Arabidopsis thaliana-Pst DC3000-chitinpentaose interaction system was chosen to explore whether OOPS can be used to dig specific RBPs under special physiological conditions. Eighty-four differential RBPs in this system were found and some of them involved in reactive oxygen species (ROS) metabolic pathway. These results showed OOPS can be applied to plants successfully and would be a useful method to identify RBPomes and specific RBPs.

Discovering the RNA-Binding Proteome of Plant Leaves with an Improved RNA Interactome Capture Method

RNA-binding proteins (RBPs) play a crucial role in regulating RNA function and fate. However, the full complement of RBPs has only recently begun to be uncovered through proteome-wide approaches such as RNA interactome capture (RIC). RIC has been applied to various cell lines and organisms, including plants, greatly expanding the repertoire of RBPs. However, several technical challenges have limited the efficacy of RIC when applied to plant tissues. Here, we report an improved version of RIC that overcomes the difficulties imposed by leaf tissue. Using this improved RIC method in Arabidopsis leaves, we identified 717 RBPs, generating a deep RNA-binding proteome for leaf tissues. While 75% of these RBPs can be linked to RNA biology, the remaining 25% were previously not known to interact with RNA. Interestingly, we observed that a large number of proteins related to photosynthesis associate with RNA in vivo, including proteins from the four major photosynthetic supercomplexes. As has previously been reported for mammals, a large proportion of leaf RBPs lack known RNA-binding domains, suggesting unconventional modes of RNA binding. We anticipate that this improved RIC method will provide critical insights into RNA metabolism in plants, including how cellular RBPs respond to environmental, physiological and pathological cues.

Sugar Signaling and Post-transcriptional Regulation in Plants: An Overlooked or an Emerging Topic?

Plants are autotrophic organisms that self-produce sugars through photosynthesis. These sugars serve as an energy source, carbon skeletons, and signaling entities throughout plants' life. Post-transcriptional regulation of gene expression plays an important role in various sugar-related processes. In cells, it is regulated by many factors, such as RNA-binding proteins (RBPs), microRNAs, the spliceosome, etc. To date, most of the investigations into sugar-related gene expression have been focused on the transcriptional level in plants, while only a few studies have been conducted on post-transcriptional mechanisms. The present review provides an overview of the relationships between sugar and post-transcriptional regulation in plants. It addresses the relationships between sugar signaling and RBPs, microRNAs, and mRNA stability. These new items insights will help to reach a comprehensive understanding of the diversity of sugar signaling regulatory networks, and open onto new investigations into the relevance of these regulations for plant growth and development.

Genomic survey of RNA recognition motif (RRM) containing RNA binding proteins from barley (Hordeum vulgare ssp. vulgare)

One of the major mechanisms of post-transcriptional gene regulation is achieved by proteins bearing well-defined sequence motifs involved in 'RNA binding'. In eukaryotes, RNA binding proteins (RBPs) are key players of RNA metabolism that includes synthesis, processing, editing, modifying, transport, storage and stability of RNA. In plants, the family of RBPs is vastly expanded compared to other eukaryotes including humans. In this study we identified 363 RBPs in the barley genome. Gene ontology enrichment analysis of barley RBPs indicated these proteins were in all the major cellular compartments and associated with key biological processes including translation, splicing, seed development and stress signaling. Members with the classical RNA binding motifs such as the RNA recognition motif (RRM), KH domain, Helicase, CRM, dsRNA and Pumilio were identified in the repertoire of barley RBPs. Similar to Arabidopsis, the RRM containing RBPs were the most abundant in barley genome. In-depth analysis of the RRM containing proteins - polyA binding proteins, Ser/Arg rich proteins and Glycine-rich RBPs were undertaken. Reanalysis of the proteome dataset of various stages during barley malting identified 38 RBPs suggesting an important role for these proteins during the malting process. This survey provides a systematic analysis of barley RBPs and serves as the basis for the further functional characterization of this important family of proteins.

Wheat differential gene expression induced by different races of Puccinia triticina

Puccinia triticina, the causal agent of wheat leaf rust, causes significant losses in wheat yield and quality each year worldwide. During leaf rust infection, the host plant recognizes numerous molecules, some of which trigger host defenses. Although P. triticina reproduces clonally, there is still variation within the population due to a high mutation frequency, host specificity, and environmental adaptation. This study explores how wheat responds on a gene expression level to different P. triticina races. Six P. triticina races were inoculated onto a susceptible wheat variety and samples were taken at six days post inoculation, just prior to pustule eruption. RNA sequence data identified 63 wheat genes differentially expressed between the six races. A time course, conducted over the first seven days post inoculation, was used to examine the expression pattern of 63 genes during infection. Forty-seven wheat genes were verified to have differential expression. Three common expression patterns were identified. In addition, two genes were associated with race specific gene expression. Differential expression of an ER molecular chaperone gene was associated with races from two different P. triticina lineages. Also, differential expression in an alanine glyoxylate aminotransferase gene was associated with races with virulence shifts for leaf rust resistance genes.

Plant Ribonomics: Proteins in Search of RNA Partners

Research into the regulation of gene expression underwent a shift from focusing on DNA-binding proteins as key transcriptional regulators to RNA-binding proteins (RBPs) that come into play once transcription has been initiated. RBPs orchestrate all RNA-processing steps in the cell. To obtain a global view of in vivo targets, the RNA complement associated with particular RBPs is determined via immunoprecipitation of the RBP and subsequent identification of bound RNAs via RNA-seq. Here, we describe technical advances in identifying RBP in vivo targets and their binding motifs. We provide an up-to-date view of targets of nucleocytoplasmic RBPs collected in arabidopsis. We also discuss current experimental limitations and provide an outlook on how the approaches may advance our understanding of post-transcriptional networks.

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.

Adaptation of iCLIP to plants determines the binding landscape of the clock-regulated RNA-binding protein AtGRP7

Background

Functions for RNA-binding proteins in orchestrating plant development and environmental responses are well established. However, the lack of a genome-wide view of their in vivo binding targets and binding landscapes represents a gap in understanding the mode of action of plant RNA-binding proteins. Here, we adapt individual nucleotide resolution crosslinking and immunoprecipitation (iCLIP) genome-wide to determine the binding repertoire of the circadian clock-regulated Arabidopsis thaliana glycine-rich RNA-binding protein AtGRP7.

Results

iCLIP identifies 858 transcripts with significantly enriched crosslink sites in plants expressing AtGRP7-GFP that are absent in plants expressing an RNA-binding-dead AtGRP7 variant or GFP alone. To independently validate the targets, we performed RNA immunoprecipitation (RIP)-sequencing of AtGRP7-GFP plants subjected to formaldehyde fixation. Of the iCLIP targets, 452 were also identified by RIP-seq and represent a set of high-confidence binders. AtGRP7 can bind to all transcript regions, with a preference for 3′ untranslated regions. In the vicinity of crosslink sites, U/C-rich motifs are overrepresented. Cross-referencing the targets against transcriptome changes in AtGRP7 loss-of-function mutants or AtGRP7-overexpressing plants reveals a predominantly negative effect of AtGRP7 on its targets. In particular, elevated AtGRP7 levels lead to damping of circadian oscillations of transcripts, including DORMANCY/AUXIN ASSOCIATED FAMILY PROTEIN2 and CCR-LIKE. Furthermore, several targets show changes in alternative splicing or polyadenylation in response to altered AtGRP7 levels.

Conclusions

We have established iCLIP for plants to identify target transcripts of the RNA-binding protein AtGRP7. This paves the way to investigate the dynamics of posttranscriptional networks in response to exogenous and endogenous cues.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1332-x) contains supplementary material, which is available to authorized users.

RNA-Binding Proteins Revisited - The Emerging Arabidopsis mRNA Interactome

RNA-protein interaction is an important checkpoint to tune gene expression at the RNA level. Global identification of proteins binding in vivo to mRNA has been possible through interactome capture - where proteins are fixed to target RNAs by UV crosslinking and purified through affinity capture of polyadenylated RNA. In Arabidopsis over 500 RNA-binding proteins (RBPs) enriched in UV-crosslinked samples have been identified. As in mammals and yeast, the mRNA interactomes came with a few surprises. For example, a plethora of the proteins caught on RNA had not previously been linked to RNA-mediated processes, for example proteins of intermediary metabolism. Thus, the studies provide unprecedented insights into the composition of the mRNA interactome, highlighting the complexity of RNA-mediated processes.

Post-transcriptional regulation of photosynthetic genes is a key driver of C4 leaf ontogeny

C₄ photosynthesis allows highly efficient carbon fixation that originates from tightly regulated anatomical and biochemical modifications of leaf architecture. Recent studies showed that leaf transcriptome modifications during leaf ontogeny of closely related C₃ (Tarenaya hassleriana) and C₄ (Gynandropsis gynandra) species within the Cleomaceae family existed but they did not identify any dedicated transcriptional networks or factors specifically driving C₄ leaf ontogeny. RNAseq analysis provides a steady-state quantification of whole-cell mRNAs but does not allow any discrimination between transcriptional and post-transcriptional processes that may occur simultaneously during leaf ontogeny. Here we use exon-intron split analysis (EISA) to determine the extent to which transcriptional and post-transcriptional processes are involved in the regulation of gene expression between young and expanded leaves in both species. C₄-specific changes in post-transcriptional regulation were observed for genes involved in the Calvin-Benson cycle and some photosystem components but not for C₄ core-cycle genes. Overall, this study provides an unbiased genome-wide insight into the post-transcriptional mechanisms that regulate gene expression through the control of mRNA levels and could be central to the onset of C₄ photosynthesis. This mechanism is cytosolic which implies cell-specific modifications of mRNA stability. Understanding this mechanism may be crucial when aiming to transform C₃ crops into C₄ crops.

In Planta Determination of the mRNA-Binding Proteome of Arabidopsis Etiolated Seedlings

RNA binding proteins (RBPs) control the fate and expression of a transcriptome. Despite this fundamental importance, our understanding of plant RBPs is rudimentary, being mainly derived via bioinformatic extrapolation from other kingdoms. Here, we adapted the mRNA-protein interactome capture method to investigate the RNA binding proteome in planta. From Arabidopsis thaliana etiolated seedlings, we captured more than 700 proteins, including 300 with high confidence that we have defined as the At-RBP set. Approximately 75% of these At-RBPs are bioinformatically linked with RNA biology, containing a diversity of canonical RNA binding domains (RBDs). As no prior experimental RNA binding evidence exists for the majority of these proteins, their capture now authenticates them as RBPs. Moreover, we identified protein families harboring emerging and potentially novel RBDs, including WHIRLY, LIM, ALBA, DUF1296, and YTH domain-containing proteins, the latter being homologous to animal RNA methylation readers. Other At-RBP set proteins include major signaling proteins, cytoskeleton-associated proteins, membrane transporters, and enzymes, suggesting the scope and function of RNA-protein interactions within a plant cell is much broader than previously appreciated. Therefore, our foundation data set has provided an unbiased insight into the RNA binding proteome of plants, on which future investigations into plant RBPs can be based.

Chilling-responsive mechanisms in halophyte Puccinellia tenuiflora seedlings revealed from proteomics analysis

Alkali grass (Puccinellia tenuiflora), a monocotyledonous perennial halophyte species, is a good pasture with great nutritional value for livestocks. It can thrive under low temperature in the saline-alkali soil of Songnen plain in northeastern China. In the present study, the chilling-responsive mechanism in P. tenuiflora leaves was investigated using physiological and proteomic approaches. After treatment of 10°C for 10 and 20days, photosynthesis, biomass, contents of osmolytes and antioxidants, and activities of reactive oxygen species scavenging enzymes were analyzed in leaves of 20-day-old seedlings. Besides, 89 chilling-responsive proteins were revealed from proteomic analysis. All the results highlighted that the growth of seedlings was inhibited due to chilling-decreased enzymes in photosynthesis, carbohydrate metabolism, and energy supplying. The accumulation of osmolytes (i.e., proline, soluble sugar, and glycine betaine) and enhancement of ascorbate-glutathione cycle and glutathione peroxidase/glutathione S-transferase pathway in leaves could minimize oxidative damage of membrane and other molecules under the chilling conditions. In addition, protein synthesis and turnover in cytoplasm and chloroplast were altered to cope with the chilling stress. This study provides valuable information for understanding the chilling-responsive and cross-tolerant mechanisms in monocotyledonous halophyte plant species.

Regulatory effects of cotranscriptional RNA structure formation and transitions

RNAs, which play significant roles in many fundamental biological processes of life, fold into sophisticated and precise structures. RNA folding is a dynamic and intricate process, which conformation transition of coding and noncoding RNAs form the primary elements of genetic regulation. The cellular environment contains various intrinsic and extrinsic factors that potentially affect RNA folding in vivo, and experimental and theoretical evidence increasingly indicates that the highly flexible features of the RNA structure are affected by these factors, which include the flanking sequence context, physiochemical conditions, cis RNA-RNA interactions, and RNA interactions with other molecules. Furthermore, distinct RNA structures have been identified that govern almost all steps of biological processes in cells, including transcriptional activation and termination, transcriptional mutagenesis, 5'-capping, splicing, 3'-polyadenylation, mRNA export and localization, and translation. Here, we briefly summarize the dynamic and complex features of RNA folding along with a wide variety of intrinsic and extrinsic factors that affect RNA folding. We then provide several examples to elaborate RNA structure-mediated regulation at the transcriptional and posttranscriptional levels. Finally, we illustrate the regulatory roles of RNA structure and discuss advances pertaining to RNA structure in plants. WIREs RNA 2016, 7:562-574. doi: 10.1002/wrna.1350 For further resources related to this article, please visit the WIREs website.

αCP binding to a cytosine-rich subset of polypyrimidine tracts drives a novel pathway of cassette exon splicing in the mammalian transcriptome

Alternative splicing (AS) is a robust generator of mammalian transcriptome complexity. Splice site specification is controlled by interactions of cis-acting determinants on a transcript with specific RNA binding proteins. These interactions are frequently localized to the intronic U-rich polypyrimidine tracts (PPT) located 5′ to the majority of splice acceptor junctions. αCPs (also referred to as polyC-binding proteins (PCBPs) and hnRNPEs) comprise a subset of KH-domain proteins with high affinity and specificity for C-rich polypyrimidine motifs. Here, we demonstrate that αCPs promote the splicing of a defined subset of cassette exons via binding to a C-rich subset of polypyrimidine tracts located 5′ to the αCP-enhanced exonic segments. This enhancement of splice acceptor activity is linked to interactions of αCPs with the U2 snRNP complex and may be mediated by cooperative interactions with the canonical polypyrimidine tract binding protein, U2AF65. Analysis of αCP-targeted exons predicts a substantial impact on fundamental cell functions. These findings lead us to conclude that the αCPs play a direct and global role in modulating the splicing activity and inclusion of an array of cassette exons, thus driving a novel pathway of splice site regulation within the mammalian transcriptome.

Transcriptome-Wide Identification of RNA Targets of Arabidopsis SERINE/ARGININE-RICH45 Uncovers the Unexpected Roles of This RNA Binding Protein in RNA Processing

Plant SR45 and its metazoan ortholog RNPS1 are serine/arginine-rich (SR)-like RNA binding proteins that function in splicing/postsplicing events and regulate diverse processes in eukaryotes. Interactions of SR45 with both RNAs and proteins are crucial for regulating RNA processing. However, in vivo RNA targets of SR45 are currently unclear. Using RNA immunoprecipitation followed by high-throughput sequencing, we identified over 4000 Arabidopsis thaliana RNAs that directly or indirectly associate with SR45, designated as SR45-associated RNAs (SARs). Comprehensive analyses of these SARs revealed several roles for SR45. First, SR45 associates with and regulates the expression of 30% of abscisic acid (ABA) signaling genes at the postsplicing level. Second, although most SARs are derived from intron-containing genes, surprisingly, 340 SARs are derived from intronless genes. Expression analysis of the SARs suggests that SR45 differentially regulates intronless and intron-containing SARs. Finally, we identified four overrepresented RNA motifs in SARs that likely mediate SR45's recognition of its targets. Therefore, SR45 plays an unexpected role in mRNA processing of intronless genes, and numerous ABA signaling genes are targeted for regulation at the posttranscriptional level. The diverse molecular functions of SR45 uncovered in this study are likely applicable to other species in view of its conservation across eukaryotes.

Rational experiment design for sequencing-based RNA structure mapping

Structure mapping is a classic experimental approach for determining nucleic acid structure that has gained renewed interest in recent years following advances in chemistry, genomics, and informatics. The approach encompasses numerous techniques that use different means to introduce nucleotide-level modifications in a structure-dependent manner. Modifications are assayed via cDNA fragment analysis, using electrophoresis or next-generation sequencing (NGS). The recent advent of NGS has dramatically increased the throughput, multiplexing capacity, and scope of RNA structure mapping assays, thereby opening new possibilities for genome-scale, de novo, and in vivo studies. From an informatics standpoint, NGS is more informative than prior technologies by virtue of delivering direct molecular measurements in the form of digital sequence counts. Motivated by these new capabilities, we introduce a novel model-based in silico approach for quantitative design of large-scale multiplexed NGS structure mapping assays, which takes advantage of the direct and digital nature of NGS readouts. We use it to characterize the relationship between controllable experimental parameters and the precision of mapping measurements. Our results highlight the complexity of these dependencies and shed light on relevant tradeoffs and pitfalls, which can be difficult to discern by intuition alone. We demonstrate our approach by quantitatively assessing the robustness of SHAPE-Seq measurements, obtained by multiplexing SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry in conjunction with NGS. We then utilize it to elucidate design considerations in advanced genome-wide approaches for probing the transcriptome, which recently obtained in vivo information using dimethyl sulfate (DMS) chemistry.

The evolution of Momordica cyclic peptides

Cyclic proteins have evolved for millions of years across all kingdoms of life to confer structural stability over their acyclic counterparts while maintaining intrinsic functional properties. Here, we show that cyclic miniproteins (or peptides) from Momordica (Cucurbitaceae) seeds evolved in species that diverged from an African ancestor around 19 Ma. The ability to achieve head-to-tail cyclization of Momordica cyclic peptides appears to have been acquired through a series of mutations in their acyclic precursor coding sequences following recent and independent gene expansion event(s). Evolutionary analysis of Momordica cyclic peptides reveals sites that are under selection, highlighting residues that are presumably constrained for maintaining their function as potent trypsin inhibitors. Molecular dynamics of Momordica cyclic peptides in complex with trypsin reveals site-specific residues involved in target binding. In a broader context, this study provides a basis for selecting Momordica species to further investigate the biosynthesis of the cyclic peptides and for constructing libraries that may be screened against evolutionarily related serine proteases implicated in human diseases.

Transcriptome-wide ribonuclease-mediated protein footprinting to identify RNA-protein interaction sites

RNA-binding proteins (RBPs) are intimately involved in all aspects of RNA processing and regulation and are linked to neurodegenerative diseases and cancer. Therefore, investigating the relationship between RBPs and their RNA targets is critical for a broader understanding of post-transcriptional regulation in normal and disease processes. The majority of approaches to study RNA-protein interactions interrogate only individual RBPs. However, there are hundreds of these proteins encoded in the human genome, and each cell type expresses a different repertoire, greatly limiting the ability of current methods to capture the global landscape of RNA-protein interactions. To address this gap, we and others have recently developed methods to globally identify regions of RNAs that are bound by proteins in an unbiased manner. Here, we describe a detailed protocol for performing our ribonuclease-mediated protein footprint sequencing approach, termed protein interaction profile sequencing (PIP-seq). In this protocol, RNA-protein interactions are stabilized by cross-linking, and unbound regions are digested with ribonucleases (RNases), leaving only the protein-bound regions intact. To control for RNase insensitive regions, proteins are first denatured and degraded, then protein-depleted RNAs are subjected to RNase treatment. After high-throughput sequencing of the remaining fragments, peak calling is performed to identify protein-protected sites (PPSs). We describe the application of this protocol to a human embryonic kidney cell line (HEK293T) and perform basic quality control, reproducibility, and benchmarking analyses. Finally, we delineate the landscape of protein-interactions in HEK293T cells, underscoring the value of this approach. Future applications of this method to study the dynamics of RNA-protein interactions in developmental and disease processes will help to further uncover the role of RBPs in post-transcriptional regulation.

IRE mRNA riboregulators use metabolic iron (Fe(2+)) to control mRNA activity and iron chemistry in animals

A family of noncoding RNAs bind Fe(2+) to increase protein synthesis. The structures occur in messenger RNAs encoding animal proteins for iron metabolism. Each mRNA regulatory sequence, ∼30 ribonucleotides long, is called an IRE (Iron Responsive Element), and folds into a bent, A-RNA helix with a terminal loop. Riboregulatory RNAs, like t-RNAs, r-RNAs micro-RNAs, etc. contrast with DNA, since single-stranded RNA can fold into a variety of complex, three-dimensional structures. IRE-RNAs bind two types of proteins: (1) IRPs which are protein repressors, sequence-related to mitochondrial aconitases. (2) eIF-4F, which bind ribosomes and enhances general protein biosynthesis. The competition between IRP and eIF-4F binding to IRE-RNA is controlled by Fe(2+)-induced changes in the IRE-RNA conformation. Mn(2+), which also binds to IRE-RNA in solution, is a convenient experimental proxy for air-sensitive Fe(2+) studies of in vitro protein biosynthesis and protein binding. However, only Fe(2+) has physiological effects on protein biosynthesis directed by IRE-mRNAs. The structures of the IRE-RNA riboregulators is known indirectly from effects of base substitutions on function, from solution NMR of the free RNA, and of X-ray crystallography of the IRE-RNA-IRP repressor complex. However, the inability to date, to crystallize the free IRE-RNA, and the dissociation of the IRE-RNA-IRP complex when metal binds, have hampered direct identification and characterization of the RNA-metal binding sites. The high conservation of the primary sequence in IRE-mRNA control elements has facilitated their identification and analysis of metal-assisted riboregulator function. Expansion of RNA search analyses beyond primary will likely reveal other, metal-dependent families of mRNA riboregulators.

Insights into RNA structure and function from genome-wide studies

A comprehensive understanding of RNA structure will provide fundamental insights into the cellular function of both coding and non-coding RNAs. Although many RNA structures have been analysed by traditional biophysical and biochemical methods, the low-throughput nature of these approaches has prevented investigation of the vast majority of cellular transcripts. Triggered by advances in sequencing technology, genome-wide approaches for probing the transcriptome are beginning to reveal how RNA structure affects each step of protein expression and RNA stability. In this Review, we discuss the emerging relationships between RNA structure and the regulation of gene expression.

The impact of RNA structure on coding sequence evolution in both bacteria and eukaryotes

Background

Many studies have found functional RNA secondary structures are selectively conserved among species. But, the effect of RNA structure selection on coding sequence evolution remains unknown. To address this problem, we systematically investigated the relationship between nucleotide conservation level and its structural sensitivity in four model organisms, Escherichia coli, yeast, fly, and mouse.

Results

We define structurally sensitive sites as those with putative local structure-disruptive mutations. Using both the Mantel-Haenszel procedure and association test, we found structurally sensitive nucleotide sites evolved more slowly than non-sensitive sites in all four organisms. Furthermore, we observed that this association is more obvious in highly expressed genes and region near the start codon.

Conclusion

We conclude that structurally sensitive sites in mRNA sequences normally have less nucleotide divergence in all species we analyzed. This study extends our understanding of the impact of RNA structure on coding sequence evolution, and is helpful to the development of a codon model with RNA structure information.

Computational analysis of conserved RNA secondary structure in transcriptomes and genomes

Transcriptomics experiments and computational predictions both enable systematic discovery of new functional RNAs. However, many putative noncoding transcripts arise instead from artifacts and biological noise, and current computational prediction methods have high false positive rates. I discuss prospects for improving computational methods for analyzing and identifying functional RNAs, with a focus on detecting signatures of conserved RNA secondary structure. An interesting new front is the application of chemical and enzymatic experiments that probe RNA structure on a transcriptome-wide scale. I review several proposed approaches for incorporating structure probing data into the computational prediction of RNA secondary structure. Using probabilistic inference formalisms, I show how all these approaches can be unified in a well-principled framework, which in turn allows RNA probing data to be easily integrated into a wide range of analyses that depend on RNA secondary structure inference. Such analyses include homology search and genome-wide detection of new structural RNAs.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

Arabidopsis mRNA secondary structure correlates with protein function and domains

RNAs fold into intricate structures that are determined by specific base pairing interactions encoded within their primary sequences. Recently, a number of transcriptome-wide studies have suggested that RNA secondary structure is a potent cis-acting regulator of numerous post-transcriptional processes in viruses and eukaryotes. However, the need for experimentally-based structure determination methods has not been well addressed. Here, we show that the regulatory significance of Arabidopsis RNA secondary structure is revealed specifically through high-throughput, sequencing-based, structure mapping data, not by computational prediction. Additionally, we find that transcripts with similar levels of secondary structure in their UTRs (5' or 3') or CDS tend to encode proteins with coherent functions. Finally, we reveal that portions of mRNAs encoding predicted protein domains are significantly more structured than those specifying inter-domain regions. In total, our findings show the utility of high-throughput, sequencing-based, structure-mapping approaches and suggest that mRNA folding regulates protein maturation and function.