The small glycine-rich RNA binding protein AtGRP7 promotes floral transition in Arabidopsis thaliana

Structural insights to the RRM-domain of the glycine-rich RNA-binding protein from Sorghum bicolor and its role in cold stress tolerance in E. coli

Sorghum bicolor Glycine-rich RNA-binding protein (SbGRBP), exhibit the ability to bind both single-stranded and double-stranded DNA. The expression of SbGRBP is regulated by heat stress, with the protein localizing to the nucleus and cytosol. The present study delves into the structure and ssDNA binding ability of its truncated version (SbGRBP1-119) which lacks glycine rich domain (GR). This protein has the ability to bind ssDNA Using Nuclear Magnetic Resonance (NMR) spectroscopy, we have revealed the secondary structure of SbGRBP1-119, highlighting the typical configuration of GRBPs with four β-sheets and two α-helices. Notably, we found two additional α-helices at the N-terminal region that seem to interact with ssDNA, a novel observation for GRBPs. Key residues crucial for ssDNA binding were identified, suggesting a specific interaction with the oligonucleotide sequence 5'-TTCTGG-3'. Preliminary assays hinted that SbGRBP1-119 might bolster E. coli resilience to cold stress, indicating a potential chaperone-like role under stress conditions. This study sheds light on the structural basis of SbGRBP1-119's interaction with nucleic acids, deepening our understanding about the role of GRBPs' in RNA metabolism and regulation.

Assessing the Role of AtGRP7 Arginine 141, a Target of Dimethylation by PRMT5, in Flowering Time Control and Stress Response

PROTEIN ARGININE METHYLTRANSFERASES (PRMTs) catalyze arginine (R) methylation that is critical for transcriptional and post-transcriptional gene regulation. In Arabidopsis, PRMT5 that catalyzes symmetric R dimethylation is best characterized. PRMT5 mutants are late-flowering and show altered responses to environmental stress. Among PRMT5 targets are Arabidopsis thaliana GLYCINE RICH RNA BINDING PROTEIN 7 (AtGRP7) and AtGRP8 that promote the transition to flowering. AtGRP7 R141 has been shown to be modified by PRMT5. Here, we tested whether this symmetric dimethylation of R141 is important for AtGRP7's physiological role in flowering time control. We constructed AtGRP7 mutant variants with non-methylable R141 (R141A, R141K). Genomic clones containing these variants complemented the late-flowering phenotype of the grp7-1 mutant to the same extent as wild-type AtGRP7. Furthermore, overexpression of AtGRP7 R141A or R141K promoted flowering similar to overexpression of the wild-type protein. Thus, flowering time does not depend on R141 and its modification. However, germination experiments showed that R141 contributes to the activity of AtGRP7 in response to abiotic stress reactions mediated by abscisic acid during early development. Immunoprecipitation of AtGRP7-GFP in the prmt5 background revealed that antibodies against dimethylated arginine still recognized AtGRP7, suggesting that additional methyltransferases may be responsible for modification of AtGRP7.

Revealing the Arabidopsis AtGRP7 mRNA binding proteome by specific enhanced RNA interactome capture

BACKGROUND

The interaction of proteins with RNA in the cell is crucial to orchestrate all steps of RNA processing. RNA interactome capture (RIC) techniques have been implemented to catalogue RNA- binding proteins in the cell. In RIC, RNA-protein complexes are stabilized by UV crosslinking in vivo. Polyadenylated RNAs and associated proteins are pulled down from cell lysates using oligo(dT) beads and the RNA-binding proteome is identified by quantitative mass spectrometry. However, insights into the RNA-binding proteome of a single RNA that would yield mechanistic information on how RNA expression patterns are orchestrated, are scarce.

RESULTS

Here, we explored RIC in Arabidopsis to identify proteins interacting with a single mRNA, using the circadian clock-regulated Arabidopsis thaliana GLYCINE-RICH RNA-BINDING PROTEIN 7 (AtGRP7) transcript, one of the most abundant transcripts in Arabidopsis, as a showcase. Seedlings were treated with UV light to covalently crosslink RNA and proteins. The AtGRP7 transcript was captured from cell lysates with antisense oligonucleotides directed against the 5'untranslated region (UTR). The efficiency of RNA capture was greatly improved by using locked nucleic acid (LNA)/DNA oligonucleotides, as done in the enhanced RIC protocol. Furthermore, performing a tandem capture with two rounds of pulldown with the 5'UTR oligonucleotide increased the yield. In total, we identified 356 proteins enriched relative to a pulldown from atgrp7 mutant plants. These were benchmarked against proteins pulled down from nuclear lysates by AtGRP7 in vitro transcripts immobilized on beads. Among the proteins validated by in vitro interaction we found the family of Acetylation Lowers Binding Affinity (ALBA) proteins. Interaction of ALBA4 with the AtGRP7 RNA was independently validated via individual-nucleotide resolution crosslinking and immunoprecipitation (iCLIP). The expression of the AtGRP7 transcript in an alba loss-of-function mutant was slightly changed compared to wild-type, demonstrating the functional relevance of the interaction.

CONCLUSION

We adapted specific RNA interactome capture with LNA/DNA oligonucleotides for use in plants using AtGRP7 as a showcase. We anticipate that with further optimization and up scaling the protocol should be applicable for less abundant transcripts.

Mapping protein-RNA binding in plants with individual-nucleotide-resolution UV cross-linking and immunoprecipitation (plant iCLIP2)

Despite crucial roles of RNA-binding proteins (RBPs) in plant physiology and development, methods for determining their transcriptome-wide binding landscape are less developed than those used in other model organisms. Cross-linking and immunoprecipitation (CLIP) methods (based on UV-mediated generation of covalent bonds between RNAs and cognate RBPs in vivo, purification of the cross-linked complexes and identification of the co-purified RNAs by high-throughput sequencing) have been applied mainly in mammalian cells growing in monolayers or in translucent tissue. We have developed plant iCLIP2, an efficient protocol for performing individual-nucleotide-resolution CLIP (iCLIP) in plants, tailored to overcome the experimental hurdles posed by plant tissue. We optimized the UV dosage to efficiently cross-link RNA and proteins in plants and expressed epitope-tagged RBPs under the control of their native promoters in loss-of-function mutants. We select epitopes for which nanobodies are available, allowing stringent conditions for immunopurification of the RNA-protein complexes to be established. To overcome the inherently high RNase content of plant cells, RNase inhibitors are added and the limited RNA fragmentation step is modified. We combine the optimized isolation of RBP-bound RNAs with iCLIP2, a streamlined protocol that greatly enhances the efficiency of library preparation for high-throughput sequencing. Plant researchers with experience in molecular biology and handling of RNA can complete this iCLIP2 protocol in ~5 d. Finally, we describe a bioinformatics workflow to determine targets of Arabidopsis RBPs from iCLIP data, covering all steps from downloading sequencing reads to identifying cross-linking events ( https://github.com/malewins/Plant-iCLIPseq ), and present the R/Bioconductor package BindingSiteFinder to extract reproducible binding sites ( https://bioconductor.org/packages/release/bioc/html/BindingSiteFinder.html ).

Arabidopsis thaliana GLYCINE RICH RNA-BINDING PROTEIN 7 interaction with its iCLIP target LHCB1.1 correlates with changes in RNA stability and circadian oscillation

The importance of RNA-binding proteins (RBPs) for plant responses to environmental stimuli and development is well documented. Insights into the portfolio of RNAs they recognize, however, clearly lack behind the understanding gathered in non-plant model organisms. Here, we characterize binding of the circadian clock-regulated Arabidopsis thaliana GLYCINE-RICH RNA-BINDING PROTEIN 7 (AtGRP7) to its target transcripts. We identified novel RNA targets from individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) data using an improved bioinformatics pipeline that will be broadly applicable to plant RBP iCLIP data. 2705 transcripts with binding sites were identified in plants expressing AtGRP7-GFP that were not recovered in plants expressing an RNA-binding dead variant or GFP alone. A conserved RNA motif enriched in uridine residues was identified at the AtGRP7 binding sites. NMR titrations confirmed the preference of AtGRP7 for RNAs with a central U-rich motif. Among the bound RNAs, circadian clock-regulated transcripts were overrepresented. Peak abundance of the LHCB1.1 transcript encoding a chlorophyll-binding protein was reduced in plants overexpressing AtGRP7 whereas it was elevated in atgrp7 mutants, indicating that LHCB1.1 was regulated by AtGRP7 in a dose-dependent manner. In plants overexpressing AtGRP7, the LHCB1.1 half-life was shorter compared to wild-type plants whereas in atgrp7 mutant plants, the half-life was significantly longer. Thus, AtGRP7 modulates circadian oscillations of its in vivo binding target LHCB1.1 by affecting RNA stability.

Phase separation of GRP7 facilitated by FERONIA-mediated phosphorylation inhibits mRNA translation to modulate plant temperature resilience

Changes in ambient temperature profoundly affect plant growth and performance. Therefore, the molecular basis of plant acclimation to temperature fluctuation is of great interest. In this study, we discovered that GLYCINE-RICH RNA-BINDING PROTEIN 7 (GRP7) contributes to cold and heat tolerance in Arabidopsis thaliana. We found that exposure to a warm temperature rapidly induces GRP7 condensates in planta, which can be reversed by transfer to a lower temperature. Cell biology and biochemical assays revealed that GRP7 undergoes liquid-liquid phase separation (LLPS) in vivo and in vitro. LLPS of GRP7 in the cytoplasm contributes to the formation of stress granules that recruit RNA, along with the translation machinery component eukaryotic initiation factor 4E1 (eIF4E1) and the mRNA chaperones COLD SHOCK PROTEIN 1 (CSP1) and CSP3, to inhibit translation. Moreover, natural variations in GRP7 affecting the residue phosphorylated by the receptor kinase FERONIA alter its capacity to undergo LLPS and correlate with the adaptation of some Arabidopsis accessions to a wider temperature range. Taken together, our findings illustrate the role of translational control mediated by GRP7 LLPS to confer plants with temperature resilience.

Epigenetic variation in early and late flowering plants of the rubber-producing Russian dandelion Taraxacum koksaghyz provides insights into the regulation of flowering time

The Russian dandelion (Taraxacum koksaghyz) grows in temperate zones and produces large amounts of poly(cis-1,4-isoprene) in its roots, making it an attractive alternative source of natural rubber. Most T. koksaghyz plants require vernalization to trigger flower development, whereas early flowering varieties that have lost their vernalization dependence are more suitable for breeding and domestication. To provide insight into the regulation of flowering time in T. koksaghyz, we induced epigenetic variation by in vitro cultivation and applied epigenomic and transcriptomic analysis to the resulting early flowering plants and late flowering controls, allowing us to identify differences in methylation patterns and gene expression that correlated with flowering. This led to the identification of candidate genes homologous to vernalization and photoperiodism response genes in other plants, as well as epigenetic modifications that may contribute to the control of flower development. Some of the candidate genes were homologous to known floral regulators, including those that directly or indirectly regulate the major flowering control gene FT. Our atlas of genes can be used as a starting point to investigate mechanisms that control flowering time in T. koksaghyz in greater detail and to develop new breeding varieties that are more suited to domestication.

Post-transcriptional regulation of grain weight and shape by the RBP-A-J-K complex in rice

RNA-binding proteins (RBPs) are components of the post-transcriptional regulatory system, but their regulatory effects on complex traits remain unknown. Using an integrated strategy involving map-based cloning, functional characterizations, and transcriptomic and population genomic analyses, we revealed that RBP-K (LOC_Os08g23120), RBP-A (LOC_Os11g41890), and RBP-J (LOC_Os10g33230) encode proteins that form an RBP-A-J-K complex that negatively regulates rice yield-related traits. Examinations of the RBP-A-J-K complex indicated RBP-K functions as a relatively non-specific RBP chaperone that enables RBP-A and RBP-J to function normally. Additionally, RBP-J most likely affects GA pathways, resulting in considerable increases in grain and panicle lengths, but decreases in grain width and thickness. In contrast, RBP-A negatively regulates the expression of genes most likely involved in auxin-regulated pathways controlling cell wall elongation and carbohydrate transport, with substantial effects on the rice grain filling process as well as grain length and weight. Evolutionarily, RBP-K is relatively ancient and highly conserved, whereas RBP-J and RBP-A are more diverse. Thus, the RBP-A-J-K complex may represent a typical functional model for many RBPs and protein complexes that function at transcriptional and post-transcriptional levels in plants and animals for increased functional consistency, efficiency, and versatility, as well as increased evolutionary potential. Our results clearly demonstrate the importance of RBP-mediated post-transcriptional regulation for the diversity of complex traits. Furthermore, rice grain yield and quality may be enhanced by introducing various complete or partial loss-of-function mutations to specific RBP genes using clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 technology and by exploiting desirable natural tri-genic allelic combinations at the loci encoding the components of the RBP-A-J-K complex through marker-assisted selection.

Are Arabidopsis thaliana plants able to recover from exposure to gamma radiation? A molecular perspective

Most plant research focuses on the responses immediately after exposure to ionizing irradiation (IR). However, it is as important to investigate how plants recover after exposure since this has a profound effect on future plant growth and development and hence on the long-term consequences of exposure to stress. This study aimed to investigate the IR-induced responses after exposure and during recovery by exposing 1-week old A. thaliana seedlings to gamma dose rates ranging from 27 to 103.7 mGy/h for 2 weeks and allowing them to recover for 4 days. A high-throughput RNAsequencing analysis was carried out. An enrichment of GO terms related to the metabolism of hormones was observed both after irradiation and during recovery at all dose rates. While plants exposed to the lowest dose rate activate defence responses after irradiation, they recover from the IR by resuming normal growth during the recovery period. Plants exposed to the intermediate dose rate invest in signalling and defence after irradiation. During recovery, in the plants exposed to the highest dose rate, fundamental metabolic processes such as photosynthesis and RNA modification were still affected. This might lead to detrimental effects in the long-term or in the next generations of those irradiated plants.

SPINDLY mediates O-fucosylation of hundreds of proteins and sugar-dependent growth in Arabidopsis

The recent discovery of SPINDLY (SPY)-catalyzed protein O-fucosylation revealed a novel mechanism for regulating nucleocytoplasmic protein functions in plants. Genetic evidence indicates the important roles of SPY in diverse developmental and physiological processes. However, the upstream signal controlling SPY activity and the downstream substrate proteins O-fucosylated by SPY remain largely unknown. Here, we demonstrated that SPY mediates sugar-dependent growth in Arabidopsis (Arabidopsis thaliana). We further identified hundreds of O-fucosylated proteins using lectin affinity chromatography followed by mass spectrometry. All the O-fucosylation events quantified in our proteomic analyses were undetectable or dramatically decreased in the spy mutants, and thus likely catalyzed by SPY. The O-fucosylome includes mostly nuclear and cytosolic proteins. Many O-fucosylated proteins function in essential cellular processes, phytohormone signaling, and developmental programs, consistent with the genetic functions of SPY. The O-fucosylome also includes many proteins modified by O-linked N-acetylglucosamine (O-GlcNAc) and by phosphorylation downstream of the target of rapamycin (TOR) kinase, revealing the convergence of these nutrient signaling pathways on key regulatory functions such as post-transcriptional/translational regulation and phytohormone responses. Our study identified numerous targets of SPY/O-fucosylation and potential nodes of crosstalk among sugar/nutrient signaling pathways, enabling future dissection of the signaling network that mediates sugar regulation of plant growth and development.

LncRNA FLAIL affects alternative splicing and represses flowering in Arabidopsis

How the noncoding genome affects cellular functions is a key biological question. A particular challenge is to distinguish the effects of noncoding DNA elements from long noncoding RNAs (lncRNAs) that coincide at the same loci. Here, we identified the flowering-associated intergenic lncRNA (FLAIL) in Arabidopsis through early flowering flail mutants. Expression of FLAIL RNA from a different chromosomal location in combination with strand-specific RNA knockdown characterized FLAIL as a trans-acting RNA molecule. FLAIL directly binds to differentially expressed target genes that control flowering via RNA-DNA interactions through conserved sequence motifs. FLAIL interacts with protein and RNA components of the spliceosome to affect target mRNA expression through co-transcriptional alternative splicing (AS) and linked chromatin regulation. In the absence of FLAIL, splicing defects at the direct FLAIL target flowering gene LACCASE 8 (LAC8) correlated with reduced mRNA expression. Double mutant analyses support a model where FLAIL-mediated splicing of LAC8 promotes its mRNA expression and represses flowering. Our study suggests lncRNAs as accessory components of the spliceosome that regulate AS and gene expression to impact organismal development.

Genome-wide screening and characterization of long noncoding RNAs involved in flowering/bolting of Lactuca sativa

BACKGROUND

Lettuce (Lactuca sativa L.) is considered the most important vegetable in the leafy vegetable group. However, bolting affects quality, gives it a bitter taste, and as a result makes it inedible. Bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components. Although bolting/flowering responsive genes have been identified in most sensitive and non-sensitive species, non-coding RNA molecules like long non-coding RNAs (lncRNAs) have not been investigated in lettuce. Hence, in this study, potential long non-coding RNAs that regulate flowering /bolting were investigated in two lettuce strains S24 (resistant strain) and S39 (susceptible strain) in different flowering times to better understand the regulation of lettuce bolting mechanism. For this purpose, we used two RNA-seq datasets to discover the lncRNA transcriptome profile during the transition from vegetative to reproductive phase.

RESULTS

For identifying unannotated transcripts in these datasets, a 7-step pipeline was employed to filter out these transcripts and terminate with 293 novel lncRNAs predicted by PLncPRO and CREMA. These transcripts were then utilized to predict cis and trans flowering-associated targets and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Computational predictions of target gene function showed the involvement of putative flowering-related genes and enrichment of the floral regulators FLC, CO, FT, and SOC1 in both datasets. Finally, 17 and 18 lncRNAs were proposed as competing endogenous target mimics (eTMs) for novel and known lncRNA miRNAs, respectively.

CONCLUSION

Overall, this study provides new insights into lncRNAs that control the flowering time of plants known for bolting, such as lettuce, and opens new windows for further study.

Selection of favorable alleles of genes controlling flowering and senescence improves malt barley quality

Malt barley (Hordeum vulgare L.) is an important cash crop with stringent grain quality standards. Timing of the switch from vegetative to reproductive growth and timing of whole-plant senescence and nutrient remobilization are critical for cereal grain yield and quality. Understanding the genetic variation in genes associated with these developmental traits can streamline genotypic selection of superior malt barley germplasm. Here, we determined the effects of allelic variation in three genes encoding a glycine-rich RNA-binding protein (HvGR-RBP1) and two NAC transcription factors (HvNAM1 and HvNAM2) on malt barley agronomics and quality using previously developed markers for HvGR-RBP1 and HvNAM1 and a novel marker for HvNAM2. Based on a single-nucleotide polymorphism (SNP) in the first intron, the utilized marker differentiates NAM2 alleles of low-grain protein variety 'Karl' and of higher protein variety 'Lewis'. We demonstrate that the selection of favorable alleles for each gene impacts heading date, senescence timing, grain size, grain protein concentration, and malt quality. Specifically, combining 'Karl' alleles for the two NAC genes with the 'Lewis' HvGR-RBP1 allele extends grain fill duration, increases the percentage of plump kernels, decreases grain protein, and provides malt quality stability. Molecular markers for these genes are therefore highly useful tools in malt barley breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01331-7.

The MYB33, MYB65, and MYB101 transcription factors affect Arabidopsis and potato responses to drought by regulating the ABA signaling pathway

Drought is one of the main climate threats limiting crop production. Potato is one of the four most important food crop species worldwide and is sensitive to water shortage. The CBP80 gene was shown to affect Arabidopsis and potato responses to drought by regulating the level of microRNA159 and, consequently, the levels of the MYB33 and MYB101 transcription factors (TFs). Here, we show that three MYB TFs, MYB33, MYB65, and MYB101, are involved in plant responses to water shortage. Their downregulation in Arabidopsis causes stomatal hyposensitivity to abscisic acid (ABA), leading to reduced tolerance to drought. Transgenic Arabidopsis and potato plants overexpressing these genes, with a mutated recognition site in miR159, show hypersensitivity to ABA and relatively high tolerance to drought conditions. Thus, the MYB33, MYB65, and MYB101 genes may be potential targets for innovative breeding to obtain crops with relatively high tolerance to drought.

Alternative splicing as a key player in the fine-tuning of the immunity response in Arabidopsis

Plants, like animals, are constantly exposed to abiotic and biotic stresses, which often inhibit plant growth and development, and cause tissue damage, disease, and even plant death. Efficient and timely response to stress requires appropriate co- and posttranscriptional reprogramming of gene expression. Alternative pre-mRNA splicing provides an important layer of this regulation by controlling the level of factors involved in stress response and generating additional protein isoforms with specific features. Recent high-throughput studies have revealed that several defence genes undergo alternative splicing that is often affected by pathogen infection. Despite extensive work, the exact mechanisms underlying these relationships are still unclear, but the contribution of alternative protein isoforms to the defence response and the role of regulatory factors, including components of the splicing machinery, have been established. Modulation of gene expression in response to stress includes alternative splicing, chromatin remodelling, histone modifications, and nucleosome occupancy. How these processes affect plant immunity is mostly unknown, but these facets open new regulatory possibilities. Here we provide an overview of the current state of knowledge and recent findings regarding the growing importance of alternative splicing in plant response to biotic stress.

Machine-Learning-Based Genome-Wide Association Studies for Uncovering QTL Underlying Soybean Yield and Its Components

A genome-wide association study (GWAS) is currently one of the most recommended approaches for discovering marker-trait associations (MTAs) for complex traits in plant species. Insufficient statistical power is a limiting factor, especially in narrow genetic basis species, that conventional GWAS methods are suffering from. Using sophisticated mathematical methods such as machine learning (ML) algorithms may address this issue and advance the implication of this valuable genetic method in applied plant-breeding programs. In this study, we evaluated the potential use of two ML algorithms, support-vector machine (SVR) and random forest (RF), in a GWAS and compared them with two conventional methods of mixed linear models (MLM) and fixed and random model circulating probability unification (FarmCPU), for identifying MTAs for soybean-yield components. In this study, important soybean-yield component traits, including the number of reproductive nodes (RNP), non-reproductive nodes (NRNP), total nodes (NP), and total pods (PP) per plant along with yield and maturity, were assessed using a panel of 227 soybean genotypes evaluated at two locations over two years (four environments). Using the SVR-mediated GWAS method, we were able to discover MTAs colocalized with previously reported quantitative trait loci (QTL) with potential causal effects on the target traits, supported by the functional annotation of candidate gene analyses. This study demonstrated the potential benefit of using sophisticated mathematical approaches, such as SVR, in a GWAS to complement conventional GWAS methods for identifying MTAs that can improve the efficiency of genomic-based soybean-breeding programs.

The Potential Role of RNA Structure in Crop Molecular Breeding

The continually growing human population creates a concomitantly increasing demand for nutritious crops with high yields. Advances in high throughput sequencing technologies have revealed the genetic architecture of major crops. This includes extensive information enabling comprehensive genetic markers for breeding selection, new gene discoveries, and novel gene regulatory strategies for crop editing. RNA structure is an important type of genetic feature, essential for post-transcriptional regulation of gene expression. Here, we summarize recent advances in genome-wide RNA structure studies in crops and review the associated RNA structure-mediated regulation of gene expression. We also discuss the functional importance of those single nucleotide variations that induce large RNA structure disparities. Lastly, we discuss the potential role of RNA structure in crop molecular breeding.

Arabidopsis apoplastic fluid contains sRNA- and circular RNA-protein complexes that are located outside extracellular vesicles

Previously, we have shown that apoplastic wash fluid (AWF) purified from Arabidopsis leaves contains small RNAs (sRNAs). To investigate whether these sRNAs are encapsulated inside extracellular vesicles (EVs), we treated EVs isolated from Arabidopsis leaves with the protease trypsin and RNase A, which should degrade RNAs located outside EVs but not those located inside. These analyses revealed that apoplastic RNAs are mostly located outside and are associated with proteins. Further analyses of these extracellular RNAs (exRNAs) revealed that they include both sRNAs and long noncoding RNAs (lncRNAs), including circular RNAs (circRNAs). We also found that exRNAs are highly enriched in the posttranscriptional modification N6-methyladenine (m6A). Consistent with this, we identified a putative m6A-binding protein in AWF, GLYCINE-RICH RNA-BINDING PROTEIN 7 (GRP7), as well as the sRNA-binding protein ARGONAUTE2 (AGO2). These two proteins coimmunoprecipitated with lncRNAs, including circRNAs. Mutation of GRP7 or AGO2 caused changes in both the sRNA and lncRNA content of AWF, suggesting that these proteins contribute to the secretion and/or stabilization of exRNAs. We propose that exRNAs located outside of EVs mediate host-induced gene silencing, rather than RNA located inside EVs.

New paradigms in cell adaptation: decades of discoveries on the CrRLK1L receptor kinase signalling network

Receptor-like kinases (RLKs), which constitute the largest receptor family in plants, are essential for perceiving and relaying information about various environmental stimuli. Tremendous progress has been made in the past few decades towards elucidating the mechanisms of action of several RLKs, with emerging paradigms pointing to their roles in cell adaptations. Among these paradigms, Catharanthus roseus receptor-like kinase 1-like (CrRLK1L) proteins and their rapid alkalinization factor (RALF) peptide ligands have attracted much interest. In particular, FERONIA (FER) is a CrRLK1L protein that participates in a wide array of physiological processes associated with RALF signalling, including cell growth and monitoring cell wall integrity, RNA and energy metabolism, and phytohormone and stress responses. Here, we analyse FER in the context of CrRLK1L members and their ligands in multiple species. The FER working model raises many questions about the role of CrRLK1L signalling networks during cell adaptation. For example, how do CrRLK1Ls recognize various RALF peptides from different organisms to initiate specific phosphorylation signal cascades? How do RALF-FER complexes achieve their specific, sometimes opposite, functions in different cell types? Here, we summarize recent major findings and highlight future perspectives in the field of CrRLK1L signalling networks.

Nonsense-mediated mRNA decay modulates Arabidopsis flowering time via the SET DOMAIN GROUP 40-FLOWERING LOCUS C module

The nonsense-mediated mRNA decay (NMD) surveillance system clears aberrant mRNAs from the cell, thus preventing the accumulation of truncated proteins. Although loss of the core NMD proteins UP-FRAMESHIFT1 (UPF1) and UPF3 leads to late flowering in Arabidopsis, the underlying mechanism remains elusive. Here, we showed that mutations in UPF1 and UPF3 cause temperature- and photoperiod-independent late flowering. Expression analyses revealed high FLOWERING LOCUS C (FLC) mRNA levels in upf mutants; in agreement with this, the flc mutation strongly suppressed the late flowering of upf mutants. Vernalization accelerated flowering of upf mutants in a temperature-independent manner. FLC transcript levels rose in wild-type plants upon NMD inhibition. In upf mutants, we observed increased enrichment of H3K4me3 and reduced enrichment of H3K27me3 in FLC chromatin. Transcriptome analyses showed that SET DOMAIN GROUP 40 (SDG40) mRNA levels increased in upf mutants, and the SDG40 transcript underwent NMD-coupled alternative splicing, suggesting that SDG40 affects flowering time in upf mutants. Furthermore, NMD directly regulated SDG40 transcript stability. The sdg40 mutants showed decreased H3K4me3 and increased H3K27me3 levels in FLC chromatin, flowered early, and rescued the late flowering of upf mutants. Taken together, these results suggest that NMD epigenetically regulates FLC through SDG40 to modulate flowering time in Arabidopsis.

The CBP/p300 histone acetyltransferases function as plant-specific MEDIATOR subunits in Arabidopsis

In eukaryotes, MEDIATOR is a conserved multi-subunit complex that links transcription factors and RNA polymerase II and that thereby facilitates transcriptional initiation. Although the composition of MEDIATOR has been well studied in yeast and mammals, relatively little is known about the composition of MEDIATOR in plants. By affinity purification followed by mass spectrometry, we identified 28 conserved MEDIATOR subunits in Arabidopsis thaliana, including putative MEDIATOR subunits that were not previously validated. Our results indicated that MED34, MED35, MED36, and MED37 are not Arabidopsis MEDIATOR subunits, as previously proposed. Our results also revealed that two homologous CBP/p300 histone acetyltransferases, HAC1 and HAC5 (HAC1/5) are in fact plant-specific MEDIATOR subunits. The MEDIATOR subunits MED8 and MED25 (MED8/25) are partially responsible for the association of MEDIATOR with HAC1/5, MED8/25 and HAC1/5 co-regulate gene expression and thereby affect flowering time and floral development. Our in vitro observations indicated that MED8 and HAC1 form liquid-like droplets by phase separation, and our in vivo observations indicated that these droplets co-localize in the nuclear bodies at a subset of nuclei. The formation of liquid-like droplets is required for MED8 to interact with RNA polymerase II. In summary, we have identified all of the components of Arabidopsis MEDIATOR and revealed the mechanism underlying the link of histone acetylation and transcriptional regulation.

The Landscape of RNA-Protein Interactions in Plants: Approaches and Current Status

RNAs transmit information from DNA to encode proteins that perform all cellular processes and regulate gene expression in multiple ways. From the time of synthesis to degradation, RNA molecules are associated with proteins called RNA-binding proteins (RBPs). The RBPs play diverse roles in many aspects of gene expression including pre-mRNA processing and post-transcriptional and translational regulation. In the last decade, the application of modern techniques to identify RNA-protein interactions with individual proteins, RNAs, and the whole transcriptome has led to the discovery of a hidden landscape of these interactions in plants. Global approaches such as RNA interactome capture (RIC) to identify proteins that bind protein-coding transcripts have led to the identification of close to 2000 putative RBPs in plants. Interestingly, many of these were found to be metabolic enzymes with no known canonical RNA-binding domains. Here, we review the methods used to analyze RNA-protein interactions in plants thus far and highlight the understanding of plant RNA-protein interactions these techniques have provided us. We also review some recent protein-centric, RNA-centric, and global approaches developed with non-plant systems and discuss their potential application to plants. We also provide an overview of results from classical studies of RNA-protein interaction in plants and discuss the significance of the increasingly evident ubiquity of RNA-protein interactions for the study of gene regulation and RNA biology in plants.

Combined effects of a glycine-rich RNA-binding protein and a NAC transcription factor extend grain fill duration and improve malt barley agronomic performance

Two key barley genes independently control anthesis and senescence timing, enabling the manipulation of grain fill duration, grain size/plumpness, and grain protein concentration. Plant developmental processes such as flowering and senescence have direct effects on cereal yield and quality. Previous work highlighted the importance of two tightly linked genes encoding a glycine-rich RNA-binding protein (HvGR-RBP1) and a NAC transcription factor (HvNAM1), controlling barley anthesis timing, senescence, and percent grain protein. Varieties that differ in HvGR-RBP1 expression, 'Karl'(low) and 'Lewis'(high), also differ in sequence 1 KB upstream of translation start site, including an ~ 400 bp G rich insertion in the 5'-flanking region of the 'Karl' allele, which could disrupt gene expression. To improve malt quality, the (low-grain protein, delayed-senescence) 'Karl' HvNAM1 allele was introgressed into Montana germplasm. After several seasons of selection, the resulting germplasm was screened for the allelic combinations of HvGR-RBP1 and HvNAM1, finding lines combining 'Karl' alleles for both genes (-/-), lines combining 'Lewis' (functional, expressed) HvGR-RBP1 with 'Karl' HvNAM1 alleles ( ±), and lines combining 'Lewis' alleles for both genes (+ / +). Field experiments indicate that the functional ('Lewis,'  +) HvGR-RBP1 allele is associated with earlier anthesis and with slightly shorter plants, while the 'Karl' (-) HvNAM1 allele delays maturation. Genotypes carrying the ± allele combination therefore had a significantly (3 days) extended grain fill duration, leading to a higher percentage of plump kernels, slightly enhanced test weight, and lower grain protein concentration when compared to the other allele combinations. Overall, our data suggest an important function for HvGR-RBP1 in the control of barley reproductive development and set the stage for a more detailed functional analysis of this gene.

Transcriptomic analysis of poco1, a mitochondrial pentatricopeptide repeat protein mutant in Arabidopsis thaliana

BACKGROUND

Flowering is a crucial stage during plant development. Plants may respond to unfavorable conditions by accelerating reproductive processes like flowering. In a recent study, we showed that PRECOCIOUS1 (POCO1) is a mitochondrial pentatricopeptide repeat (PPR) protein involved in flowering time and abscisic acid (ABA) signaling in Arabidopsis thaliana. Here, we use RNA-seq data to investigate global gene expression alteration in the poco1 mutant.

RESULTS

RNA-seq analysis was performed during different developmental stages for wild-type and poco1 plants. The most profound differences in gene expression were found when wild-type and poco1 plants of the same developmental stage were compared. Coverage analysis confirmed the T-DNA insertion in POCO1, which was concomitant with truncated transcripts. Many biological processes were found to be enriched. Several flowering-related genes such as FLOWERING LOCUS T (FT), which may be involved in the early-flowering phenotype of poco1, were differentially regulated. Numerous ABA-associated genes, including the core components of ABA signaling such as ABA receptors, protein phosphatases, protein kinases, and ABA-responsive element (ABRE) binding proteins (AREBs)/ABRE-binding factors (ABFs) as well as important genes for stomatal function, were mostly down-regulated in poco1. Drought and oxidative stress-related genes, including ABA-induced stress genes, were differentially regulated. RNA-seq analysis also uncovered differentially regulated genes encoding various classes of transcription factors and genes involved in cellular signaling. Furthermore, the expression of stress-associated nuclear genes encoding mitochondrial proteins (NGEMPs) was found to be altered in poco1. Redox-related genes were affected, suggesting that the redox state in poco1 might be altered.

CONCLUSION

The identification of various enriched biological processes indicates that complex regulatory mechanisms underlie poco1 development. Differentially regulated genes associated with flowering may contribute to the early-flowering phenotype of poco1. Our data suggest the involvement of POCO1 in the early ABA signaling process. The down-regulation of many ABA-related genes suggests an association of poco1 mutation with the ABA signaling deficiency. This condition further affects the expression of many stress-related, especially drought-associated genes in poco1, consistent with the drought sensitivity of poco1. poco1 mutation also affects the expression of genes associated with the cellular regulation, redox, and mitochondrial perturbation.

RALF1-FERONIA complex affects splicing dynamics to modulate stress responses and growth in plants

The environmentally responsive signaling pathways that link global transcriptomic changes through alternative splicing (AS) to plant fitness remain unclear. Here, we found that the interaction of the extracellular rapid alkalinization FACTOR 1 (RALF1) peptide with its receptor FERONIA (FER) triggered a rapid and massive RNA AS response by interacting with and phosphorylating glycine-rich RNA binding protein7 (GRP7) to elevate GRP7 nuclear accumulation in Arabidopsis thaliana. FER-dependent GRP7 phosphorylation enhanced its mRNA binding ability and its association with the spliceosome component U1-70K to enable splice site selection, modulating dynamic AS. Genetic reversal of a RALF1-FER-dependent splicing target partly rescued mutants deficient in GRP7. AS of GRP7 itself induced nonsense-mediated decay feedback to the RALF1-FER-GRP7 module, fine-tuning stress responses, and cell growth. The RALF1-FER-GRP7 module provides a paradigm for regulatory mechanisms of RNA splicing in response to external stimuli.

Genome-Wide Transcriptomic Analysis Reveals a Regulatory Network of Oxidative Stress-Induced Flowering Signals Produced in Litchi Leaves

Litchi is an important subtropical fruit tree that requires an appropriately low temperature to trigger floral initiation. Our previous studies have shown that reactive oxygen species (ROS) are involved in litchi flowering. To identify oxidative stress-induced flowering related genes in leaves, ‘Nuomici’ potted trees were grown at medium low-temperature conditions (18/13 °C for day/night, medium-temperature). The trees were treated with the ROS generator methyl viologen dichloride hydrate (MV) as the MV-generated ROS treatment (MM, medium-temperature plus MV) and water as the control treatment (M, medium-temperature plus water). Sixteen RNA-sequencing libraries were constructed, and each library generated more than 5,000,000 clean reads. A total of 517 differentially expressed genes (DEGs) were obtained. Among those DEGs, plant hormone biosynthesis and signal transduction genes, ROS-specific transcription factors, such as AP2/ERF and WRKY genes, stress response genes, and flowering-related genes FLOWERING LOCUS T1 (FT1) and FLOWERING LOCUS T2 (FT2) were significantly enriched. Then, as a confirmatory experiment, the potted trees were uniformly sprayed with MV, N,N’-dimethylthiourea (DMTU, ROS scavenger) plus MV, and water at medium-temperature. The results showed that the MV-generated ROS promoted flowering and changed related gene expression, but these effects were repressed by DMTU treatment. The results of our studies indicate that ROS could promote flowering and partly bypass chilling for litchi flowering.

KHZ1 and KHZ2, novel members of the autonomous pathway, repress the splicing efficiency of FLC pre-mRNA in Arabidopsis

As one of the most important events during the life cycle of flowering plants, the floral transition is of crucial importance for plant propagation and requires the precise coordination of multiple endogenous and external signals. There have been at least four flowering pathways (i.e. photoperiod, vernalization, gibberellin, and autonomous) identified in Arabidopsis. We previously reported that two Arabidopsis RNA-binding proteins, KHZ1 and KHZ2, redundantly promote flowering. However, the underlying mechanism was unclear. Here, we found that the double mutant khz1 khz2 flowered late under both long-day and short-day conditions, but responded to vernalization and gibberellin treatments. The late-flowering phenotype was almost completely rescued by mutating FLOWERING LOCUS C (FLC) and fully rescued by overexpressing FLOWERING LOCUS T (FT). Additional experiments demonstrated that the KHZs could form homodimers or interact to form heterodimers, localized to nuclear dots, and repressed the splicing efficiency of FLC pre-mRNA. Together, these data indicate that the KHZs could promote flowering via the autonomous pathway by repressing the splicing efficiency of FLC pre-mRNA.

CLIP and RNA interactome studies to unravel genome-wide RNA-protein interactions in vivo in Arabidopsis thaliana

Post-transcriptional regulation makes an important contribution to adjusting the transcriptome to environmental changes in plants. RNA-binding proteins are key players that interact specifically with mRNAs to co-ordinate their fate. While the regulatory interactions between proteins and RNA are well understood in animals, until recently little information was available on the global binding landscape of RNA-binding proteins in higher plants. This is not least due to technical challenges in plants. In turn, while numerous RNA-binding proteins have been identified through mutant analysis and homology-based searches in plants, only recently a full compendium of proteins with RNA-binding activity has been experimentally determined for the reference plant Arabidopsis thaliana. State-of-the-art techniques to determine RNA-protein interactions genome-wide in animals are based on the covalent fixation of RNA and protein in vivo by UV light. This has only recently been successfully applied to plants. Here, we present practical considerations on the application of UV irradiation based methods to comprehensively determine in vivo RNA-protein interactions in Arabidopsis thaliana, focussing on individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) and mRNA interactome capture.

Regulation of Flowering Time by the RNA-Binding Proteins AtGRP7 and AtGRP8

The timing of floral initiation is a tightly controlled process in plants. The circadian clock regulated glycine-rich RNA-binding protein (RBP) AtGRP7, a known regulator of splicing, was previously shown to regulate flowering time mainly by affecting the MADS-box repressor FLOWERING LOCUS C (FLC). Loss of AtGRP7 leads to elevated FLC expression and late flowering in the atgrp7-1 mutant. Here, we analyze genetic interactions of AtGRP7 with key regulators of the autonomous and the thermosensory pathway of floral induction. RNA interference- mediated reduction of the level of the paralogous AtGRP8 in atgrp7-1 further delays floral transition compared of with atgrp7-1. AtGRP7 acts in parallel to FCA, FPA and FLK in the branch of the autonomous pathway (AP) comprised of RBPs. It acts in the same branch as FLOWERING LOCUS D, and AtGRP7 loss-of-function mutants show elevated levels of dimethylated lysine 4 of histone H3, a mark for active transcription. In addition to its role in the AP, AtGRP7 acts in the thermosensory pathway of flowering time control by regulating alternative splicing of the floral repressor FLOWERING LOCUS M (FLM). Overexpression of AtGRP7 selectively favors the formation of the repressive isoform FLM-β. Our results suggest that the RBPs AtGRP7 and AtGRP8 influence MADS-Box transcription factors in at least two different pathways of flowering time control. This highlights the importance of RBPs to fine-tune the integration of varying cues into flowering time control and further strengthens the view that the different pathways, although genetically separable, constitute a tightly interwoven network to ensure plant reproductive success under changing environmental conditions.

Genome-wide analysis of alternative splicing divergences between Brassica hexaploid and its parents

Compared with its parents, Brassica hexaploid underwent significant AS changes, which may provide diversified gene expression regulation patterns and could enhance its adaptability during evolution Polyploidization is considered a significant evolution force that promotes species formation. Alternative splicing (AS) plays a crucial role in multiple biological processes during plant growth and development. To explore the effects of allopolyploidization on the AS patterns of genes, a genome-wide AS analysis was performed by RNA-seq in Brassica hexaploid and its parents. In total, we found 7913 (27540 AS events), 14447 (70179 AS events), and 13205 (60804 AS events) AS genes in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. A total of 920 new AS genes were discovered in Brassica hexaploid. There were 56 differently spliced genes between Brassica hexaploid and its parents. In addition, most of the alternative 5' splice sites were located 4 bp upstream of the dominant 5' splice sites, and most of the alternative 3' splice sites were located 3 bp downstream of the dominant 3' splice sites in Brassica hexapliod, which was similar to B. carinata. Furthermore, we cloned and sequenced all amplicons from the RT-PCR products of GRP7/8, namely, Bol045859, Bol016025 and Bol02880. The three genes were found to produce AS transcripts in a new way. The AS patterns of genes were diverse between Brassica hexaploid and its parents, including the loss and gain of AS events. Allopolyploidization changed alternative splicing sites of pre-mRNAs in Brassica hexaploid, which brought about alterations in the sequences of transcripts. Our study provided novel insights into the AS patterns of genes in allopolyploid plants, which may provide a reference for the study of polyploidy adaptability.

Gene regulatory network and abundant genetic variation play critical roles in heading stage of polyploidy wheat

BACKGROUND

The extensive adaptability of polyploidy wheat is attributed to its complex genome, and accurately controlling heading stage is a prime target in wheat breeding process. Wheat heading stage is an essential growth and development processes since it starts at a crucial point in the transition from vegetative phase to reproductive phase.

MAIN BODY

Heading stage is mainly decided by vernalization, photoperiod, hormone (like gibberellic acid, GA), and earliness per se (Eps). As a polyploidy species, common wheat possesses the abundant genetic variation, such as allelic variation, copy number variation etc., which have a strong effect on regulation of wheat growth and development. Therefore, understanding genetic manipulation of heading stage is pivotal for controlling the heading stage in wheat. In this review, we summarized the recent advances in the genetic regulatory mechanisms and abundant variation in genetic diversity controlling heading stage in wheat, as well as the interaction mechanism of different signals and the contribution of different genetic variation. We first summarized the genes involved in vernalization, photoperoid and other signals cross-talk with each other to control wheat heading stage, then the abundant genetic variation related to signal components associated with wheat heading stage was also elaborated in detail.

CONCLUSION

Our knowledge of the regulatory network of wheat heading can be used to adjust the duration of the growth phase for the purpose of acclimatizing to different geographical environments.

RNA Splicing of FLC Modulates the Transition to Flowering

Flowering is a critical stage of plant development and is closely correlated with seed production and crop yield. Flowering transition is regulated by complex genetic networks in response to endogenous and environmental signals. FLOWERING LOCUS C (FLC) is a central repressor in the flowering transition of Arabidopsis thaliana. The regulation of FLC expression is well studied at transcriptional and post-transcriptional levels. A subset of antisense transcripts from FLC locus, collectively termed cold-induced long antisense intragenic RNAs (COOLAIR), repress FLC expression under cold exposure. Recent studies have provided important insights into the alternative splicing of COOLAIR and FLC sense transcripts in response to developmental and environmental cues. Herein, at the 20th anniversary of FLC functional identification, we summarise new research advances in the alternative splicing of FLC sense and antisense transcripts that regulates flowering.

Identification of putative flowering genes and transcription factors from flower de novo transcriptome dataset of tuberose (Polianthes tuberosa L.)

Polianthes tuberosa is commercially popular because of their economic importance in floriculture for cut and loose flowers and in perfume industry because of the unique fragrance. Despite its commercial importance, no ready-to-use transcript sequence information is available in the public database. We have sequenced the RNA obtained from tuberose flowers using the Illumina HiSeq. 2000 platform and have carried out a de novo analysis of the transcriptome data. The de novo assembly generated 11,100 transcripts. These transcripts represent a total of 7876 unigenes that were considered for downstream analysis. These 7876 unigenes, which was further annotated using blast2go and KEGG pathways, were also assigned. Tuberose transcripts were also assigned to metabolic pathways using the Kyoto Encyclopedia of Genes and Genomes database to determine their biochemical functions. 4591 of the tuberose transcripts matched to genes in KEGG pathways and 66 transcripts were mapped to the Flavonoid biosynthesis pathway. 21 flowering genes have been identified in this tuberose transcriptome. Transcription factor analysis helped in the identification of a large number of transcripts similar to key genes in the flowering regulation network of Arabidopsis thaliana. Among the transcription factors identified “NAC” which is associated with plant stress response represented the most abundant category followed by APETALA2 (AP2)/ethylene-responsive element binding proteins (EREBPs) which plays various role in floral organ identity and respond to different biotic and abiotic stress.

A Glycine-Rich RNA-Binding Protein, CsGR-RBP3, Is Involved in Defense Responses Against Cold Stress in Harvested Cucumber (Cucumis sativus L.) Fruit

Plant glycine-rich RNA-binding proteins (GR-RBPs) have been shown to play important roles in response to abiotic stresses in actively proliferating organs such as young plants, root tips, and flowers, but their roles in chilling responses of harvested fruit remains largely unknown. Here, we investigated the role of CsGR-RBP3 in the chilling response of cucumber fruit. Pre-storage cold acclimation at 10°C (PsCA) for 3 days significantly enhanced chilling tolerance of cucumber fruit compared with the control fruit that were stored at 5°C. In the control fruit, only one of the six cucumber CsGR-RBP genes, CsGR-RBP2, was enhanced whereas the other five, i.e., CsGR-RBP3, CsGR-RBP4, CsGR-RBP5, CsGR-RBP-blt801, and CsGR-RBP-RZ1A were not. However, in the fruit exposed to PsCA before storage at 5°C, CsGR-RBP2 transcript levels were not obviously different from those in the controls, whereas the other five were highly upregulated, with CsGR-RBP3 the most significantly induced. Treatment with endogenous ABA and NO biosynthesis inhibitors, tungstate and L-nitro-arginine methyl ester, respectively, prior to PsCA treatment, clearly downregulated CsGR-RBP3 expression and significantly aggravated chilling injury. These results suggest a strong connection between CsGR-RBP3 expression and chilling tolerance in cucumber fruit. Transient expression in tobacco suggests CsGR-RBP3 was located in the mitochondria, implying a role for CsGR-RBP3 in maintaining mitochondria-related functions under low temperature. Arabidopsis lines overexpressing CsGR-RBP3 displayed faster growth at 23°C, lower electrolyte leakage and higher Fv/Fm ratio at 0°C, and higher survival rate at -20°C, than wild-type plants. Under cold stress conditions, Arabidopsis plants overexpressing CsGR-RBP3 displayed lower reactive oxygen species levels, and higher catalase and superoxide dismutase expression and activities, compared with the wild-type plants. In addition, overexpression of CsGR-RBP3 significantly upregulated nine Arabidopsis genes involved in defense responses to various stresses, including chilling. These results strongly suggest CsGR-RBP3 plays a positive role in defense against chilling stress.

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.

Computational exploration of cis-regulatory modules in rhythmic expression data using the "Exploration of Distinctive CREs and CRMs" (EDCC) and "CRM Network Generator" (CNG) programs

Understanding the effect of cis-regulatory elements (CRE) and clusters of CREs, which are called cis-regulatory modules (CRM), in eukaryotic gene expression is a challenge of computational biology. We developed two programs that allow simple, fast and reliable analysis of candidate CREs and CRMs that may affect specific gene expression and that determine positional features between individual CREs within a CRM. The first program, "Exploration of Distinctive CREs and CRMs" (EDCC), correlates candidate CREs and CRMs with specific gene expression patterns. For pairs of CREs, EDCC also determines positional preferences of the single CREs in relation to each other and to the transcriptional start site. The second program, "CRM Network Generator" (CNG), prioritizes these positional preferences using a neural network and thus allows unbiased rating of the positional preferences that were determined by EDCC. We tested these programs with data from a microarray study of circadian gene expression in Arabidopsis thaliana. Analyzing more than 1.5 million pairwise CRE combinations, we found 22 candidate combinations, of which several contained known clock promoter elements together with elements that had not been identified as relevant to circadian gene expression before. CNG analysis further identified positional preferences of these CRE pairs, hinting at positional information that may be relevant for circadian gene expression. Future wet lab experiments will have to determine which of these combinations confer daytime specific circadian gene expression.

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.

TAF15b, involved in the autonomous pathway for flowering, represses transcription of FLOWERING LOCUS C

TATA-binding protein-associated factors (TAFs) are general transcription factors within the transcription factor IID (TFIID) complex, which recognizes the core promoter of genes. In addition to their biochemical function, it is known that several TAFs are involved in the regulation of developmental processes. In this study, we found that TAF15b affects flowering time, especially through the autonomous pathway (AP) in Arabidopsis. The mutant taf15b shows late flowering compared with the wild type plant during both long and short days, and vernalization accelerates the flowering time of taf15b. In addition, taf15b shows strong upregulation of FLOWERING LOCUS C (FLC), a flowering repressor in Arabidopsis, and the flc taf15b double mutant completely offsets the late flowering of taf15b, indicating that TAF15b is a typical AP gene. The taf15b mutant also shows increased transcript levels of COOLAIR, an antisense transcript of FLC. Consistently, chromatin immunoprecipitation (ChIP) analyses showed that the TAF15b protein is enriched around both sense and antisense transcription start sites of the FLC locus. In addition, co-immunoprecipitation showed that TAF15b interacts with RNA polymerase II (Pol II), while ChIP showed increased enrichment of the phosphorylated forms, both serine 2 (Ser2) and Ser5, of the C-terminal domain of Pol II at the FLC locus, which is indicative of transcriptional elongation. Finally, taf15b showed higher enrichment of the active histone marker, H3K4me3, on FLC chromatin. Taken together, our results suggest that TAF15b affects flowering time through transcriptional repression of FLC in Arabidopsis.

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.

A photo-responsive F-box protein FOF2 regulates floral initiation by promoting FLC expression in Arabidopsis

Floral initiation is regulated by various genetic pathways in response to light, temperature, hormones and developmental status; however, the molecular mechanisms underlying the interactions between different genetic pathways are not fully understood. Here, we show that the photoresponsive gene FOF2 (F-box of flowering 2) negatively regulates flowering. FOF2 encodes a putative F-box protein that interacts specifically with ASK14, and its overexpression results in later flowering under both long-day and short-day photoperiods. Conversely, transgenic plants expressing the F-box domain deletion mutant of FOF2 (FOF2ΔF), or double loss of function mutant of FOF2 and FOL1 (FOF2-LIKE 1) present early flowering phenotypes. The late flowering phenotype of the FOF2 overexpression lines is suppressed by the flc-3 loss-of-function mutation. Furthermore, FOF2 mRNA expression is regulated by autonomous pathway gene FCA, and the repressive effect of FOF2 in flowering can be overcome by vernalization. Interestingly, FOF2 expression is regulated by light. The protein level of FOF2 accumulates in response to light, whereas it is degraded under dark conditions via the 26S proteasome pathway. Our findings suggest a possible mechanistic link between light conditions and the autonomous floral promotion pathway in Arabidopsis.

Structural disorder in plant proteins: where plasticity meets sessility

Plants are sessile organisms. This intriguing nature provokes the question of how they survive despite the continual perturbations caused by their constantly changing environment. The large amount of knowledge accumulated to date demonstrates the fascinating dynamic and plastic mechanisms, which underpin the diverse strategies selected in plants in response to the fluctuating environment. This phenotypic plasticity requires an efficient integration of external cues to their growth and developmental programs that can only be achieved through the dynamic and interactive coordination of various signaling networks. Given the versatility of intrinsic structural disorder within proteins, this feature appears as one of the leading characters of such complex functional circuits, critical for plant adaptation and survival in their wild habitats. In this review, we present information of those intrinsically disordered proteins (IDPs) from plants for which their high level of predicted structural disorder has been correlated with a particular function, or where there is experimental evidence linking this structural feature with its protein function. Using examples of plant IDPs involved in the control of cell cycle, metabolism, hormonal signaling and regulation of gene expression, development and responses to stress, we demonstrate the critical importance of IDPs throughout the life of the plant.

Comparative proteomic analysis of eggplant (Solanum melongena L.) heterostylous pistil development

Heterostyly is a common floral polymorphism, but the proteomic basis of this trait is still largely unexplored. In this study, self- and cross-pollination of L-morph and S-morph flowers and comparison of embryo sac development in eggplant (Solanum melongena L.) suggested that lower fruit set from S-morph flowers results from stigma-pollen incompatibility. To explore the molecular mechanism underlying heterostyly development, we conducted isobaric tags for relative and absolute quantification (iTRAQ) proteomic analysis of eggplant pistils for L- and S-morph flowers. A total of 5,259 distinct proteins were identified during heterostyly development. Compared S-morph flowers with L-morph, we discovered 57 and 184 differentially expressed proteins (DEPs) during flower development and maturity, respectively. Quantitative real time polymerase chain reactions were used for nine genes to verify DEPs from the iTRAQ approach. During flower development, DEPs were mainly involved in morphogenesis, biosynthetic processes, and metabolic pathways. At flower maturity, DEPs primarily participated in biosynthetic processes, metabolic pathways, and the formation of ribosomes and proteasomes. Additionally, some proteins associated with senescence and programmed cell death were found to be upregulated in S-morph pistils, which may lead to the lower fruit set in S-morph flowers. Although the exact roles of these related proteins are not yet known, this was the first attempt to use an iTRAQ approach to analyze proteomes of heterostylous eggplant flowers, and these results will provide insights into biochemical events taking place during the development of heterostyly.

A Global View of RNA-Protein Interactions Identifies Post-transcriptional Regulators of Root Hair Cell Fate

The Arabidopsis thaliana root epidermis is comprised of two cell types, hair and nonhair cells, which differentiate from the same precursor. Although the transcriptional programs regulating these events are well studied, post-transcriptional factors functioning in this cell fate decision are mostly unknown. Here, we globally identify RNA-protein interactions and RNA secondary structure in hair and nonhair cell nuclei. This analysis reveals distinct structural and protein binding patterns across both transcriptomes, allowing identification of differential RNA binding protein (RBP) recognition sites. Using these sequences, we identify two RBPs that regulate hair cell development. Specifically, we find that SERRATE functions in a microRNA-dependent manner to inhibit hair cell fate, while also terminating growth of root hairs mostly independent of microRNA biogenesis. In addition, we show that GLYCINE-RICH PROTEIN 8 promotes hair cell fate while alleviating phosphate starvation stress. In total, this global analysis reveals post-transcriptional regulators of plant root epidermal cell fate.

Identification of reference genes for quantitative expression analysis using large-scale RNA-seq data of Arabidopsis thaliana and model crop plants

In quantitative gene expression analysis, normalization using a reference gene as an internal control is frequently performed for appropriate interpretation of the results. Efforts have been devoted to exploring superior novel reference genes using microarray transcriptomic data and to evaluating commonly used reference genes by targeting analysis. However, because the number of specifically detectable genes is totally dependent on probe design in the microarray analysis, exploration using microarray data may miss some of the best choices for the reference genes. Recently emerging RNA sequencing (RNA-seq) provides an ideal resource for comprehensive exploration of reference genes since this method is capable of detecting all expressed genes, in principle including even unknown genes. We report the results of a comprehensive exploration of reference genes using public RNA-seq data from plants such as Arabidopsis thaliana (Arabidopsis), Glycine max (soybean), Solanum lycopersicum (tomato) and Oryza sativa (rice). To select reference genes suitable for the broadest experimental conditions possible, candidates were surveyed by the following four steps: (1) evaluation of the basal expression level of each gene in each experiment; (2) evaluation of the expression stability of each gene in each experiment; (3) evaluation of the expression stability of each gene across the experiments; and (4) selection of top-ranked genes, after ranking according to the number of experiments in which the gene was expressed stably. Employing this procedure, 13, 10, 12 and 21 top candidates for reference genes were proposed in Arabidopsis, soybean, tomato and rice, respectively. Microarray expression data confirmed that the expression of the proposed reference genes under broad experimental conditions was more stable than that of commonly used reference genes. These novel reference genes will be useful for analyzing gene expression profiles across experiments carried out under various experimental conditions.

Proteomic Response of Hordeum vulgare cv. Tadmor and Hordeum marinum to Salinity Stress: Similarities and Differences between a Glycophyte and a Halophyte

Response to a high salinity treatment of 300 mM NaCl was studied in a cultivated barley Hordeum vulgare Syrian cultivar Tadmor and in a halophytic wild barley H. marinum. Differential salinity tolerance of H. marinum and H. vulgare is underlied by qualitative and quantitative differences in proteins involved in a variety of biological processes. The major aim was to identify proteins underlying differential salinity tolerance between the two barley species. Analyses of plant water content, osmotic potential and accumulation of proline and dehydrin proteins under high salinity revealed a relatively higher water saturation deficit in H. marinum than in H. vulgare while H. vulgare had lower osmotic potential corresponding with high levels of proline and dehydrins. Analysis of proteins soluble upon boiling isolated from control and salt-treated crown tissues revealed similarities as well as differences between H. marinum and H. vulgare. The similar salinity responses of both barley species lie in enhanced levels of stress-protective proteins such as defense-related proteins from late-embryogenesis abundant family, several chaperones from heat shock protein family, and others such as GrpE. However, there have also been found significant differences between H. marinum and H. vulgare salinity response indicating an active stress acclimation in H. marinum while stress damage in H. vulgare. An active acclimation to high salinity in H. marinum is underlined by enhanced levels of several stress-responsive transcription factors from basic leucine zipper and nascent polypeptide-associated complex families. In salt-treated H. marinum, enhanced levels of proteins involved in energy metabolism such as glycolysis, ATP metabolism, and photosynthesis-related proteins indicate an active acclimation to enhanced energy requirements during an establishment of novel plant homeostasis. In contrast, changes at proteome level in salt-treated H. vulgare indicate plant tissue damage as revealed by enhanced levels of proteins involved in proteasome-dependent protein degradation and proteins related to apoptosis. The results of proteomic analysis clearly indicate differential responses to high salinity and provide more profound insight into biological mechanisms underlying salinity response between two barley species with contrasting salinity tolerance.

Comparative proteomic analysis of the shoot apical meristem in maize between a ZmCCT-associated near-isogenic line and its recurrent parent

The ZmCCT, one of the most important genes affecting photoperiod response, delays flowering under long-day conditions in maize (Zea mays). In this study we used the isobaric tags for relative and absolute quantification (iTRAQ) technique-based proteomics approach to identify differentially expressed proteins between a near-isogenic line (NIL) and its recurrent parent, contrasting in alleles of ZmCCT. A total of 5,259 distinct proteins were identified. Among them, 386 proteins were differentially expressed between NIL-cml line (ZmCCT-positive) and H4 line (ZmCCT-negative). Functional categorization showed that the differentially proteins were mainly involved in energy production, photosynthesis, signal transduction, and cell organization and biogenesis. Our results showed that during shoot apical meristem (SAM) development cell division proteins, carbohydrate metabolism–related proteins, and flower inhibition-related proteins were more abundant in the ZmCCT-positive line than the ZmCCT-negative line. These results, taken together with morphological observations, showed that the effect of ZmCCT on flowering might be caused by its effect on one or all of these biological processes. Although the exact roles of these putative related proteins remain to be examined, our results obtained using the proteomics approach lead to a better understanding of the photoperiodicity mechanism in maize plants.

Integrating roots into a whole plant network of flowering time genes in Arabidopsis thaliana

Molecular data concerning the involvement of roots in the genetic pathways regulating floral transition are lacking. In this study, we performed global analyses of the root transcriptome in Arabidopsis in order to identify flowering time genes that are expressed in the roots and genes that are differentially expressed in the roots during the induction of flowering. Data mining of public microarray experiments uncovered that about 200 genes whose mutations are reported to alter flowering time are expressed in the roots (i.e. were detected in more than 50% of the microarrays). However, only a few flowering integrator genes passed the analysis cutoff. Comparison of root transcriptome in short days and during synchronized induction of flowering by a single 22-h long day revealed that 595 genes were differentially expressed. Enrichment analyses of differentially expressed genes in root tissues, gene ontology categories, and cis-regulatory elements converged towards sugar signaling. We concluded that roots are integrated in systemic signaling, whereby carbon supply coordinates growth at the whole plant level during the induction of flowering. This coordination could involve the root circadian clock and cytokinin biosynthesis as a feed forward loop towards the shoot.

Genomic dissection of plant development and its impact on thousand grain weight in barley through nested association mapping

Flowering time is a key agronomic trait that plays an important role in crop yield. There is growing interest in dissecting the developmental subphases of flowering to better understand and fine-tune plant development and maximize yield. To do this, we used the wild barley nested association mapping (NAM) population HEB-25, comprising 1420 BC1S3 lines, to map quantitative trait loci (QTLs) controlling five developmental traits, plant height, and thousand grain weight. Genome-wide association studies (GWAS) enabled us to locate a total of 89 QTLs that genetically regulate the seven investigated traits. Several exotic QTL alleles proved to be highly effective and potentially useful in barley breeding. For instance, thousand grain weight was increased by 4.5 g and flowering time was reduced by 9.3 days by substituting Barke elite QTL alleles for exotic QTL alleles at the denso/sdw1 and the Ppd-H1 loci, respectively. We showed that the exotic allele at the semi-dwarf locus denso/sdw1 can be used to increase grain weight since it uncouples the negative correlation between shoot elongation and the ripening phase. Our study demonstrates that nested association mapping of HEB-25 can help unravel the genetic regulation of plant development and yield formation in barley. Moreover, since we detected numerous useful exotic QTL alleles in HEB-25, we conclude that the introgression of these wild barley alleles into the elite barley gene pool may enable developmental phases to be specifically fine-tuned in order to maximize thousand grain weight and, potentially, yield in the long term.