Genome analysis: RNA recognition motif (RRM) and K homology (KH) domain RNA-binding proteins from the flowering plant Arabidopsis thaliana

Divergence in regulatory mechanisms of GR-RBP genes in different plants under abiotic stress

The IVa subfamily of glycine-rich proteins (GRPs) comprises a group of glycine-rich RNA binding proteins referred to as GR-RBPa here. Previous studies have demonstrated functions of GR-RBPa proteins in regulating stress response in plants. However, the mechanisms responsible for the differential regulatory functions of GR-RBPa proteins in different plant species have not been fully elucidated. In this study, we identified and comprehensively studied a total of 34 GR-RBPa proteins from five plant species. Our analysis revealed that GR-RBPa proteins were further classified into two branches, with proteins in branch I being relatively more conserved than those in branch II. When subjected to identical stresses, these genes exhibited intensive and differential expression regulation in different plant species, corresponding to the enrichment of cis-acting regulatory elements involving in environmental and internal signaling in these genes. Unexpectedly, all GR-RBPa genes in branch I underwent intensive alternative splicing (AS) regulation, while almost all genes in branch II were only constitutively spliced, despite having more introns. This study highlights the complex and divergent regulations of a group of conserved RNA binding proteins in different plants when exposed to identical stress conditions. These species-specific regulations may have implications for stress responses and adaptations in different plant species.

A K homology (KH) domain protein identified by a forward genetic screen affects bundle sheath anatomy in Arabidopsis thaliana

Because of their photosynthetic capacity, leaves function as solar panels providing the basis for the growth of the entire plant. Although the molecular mechanisms of leaf development have been well studied in model dicot and monocot species, a lot of information is still needed about the interplay of the genes that regulate cell division and differentiation and thereby affect the photosynthetic performance of the leaf. We were specifically interested in understanding the differentiation of mesophyll and bundle sheath cells in Arabidopsis thaliana and aimed to identify genes that are involved in determining bundle sheath anatomy. To this end, we established a forward genetic screen by using ethyl methanesulfonate (EMS) for mutagenizing a reporter line expressing a chloroplast-targeted green fluorescent protein (sGFP) under the control of a bundle sheath-specific promoter. Based on the GFP fluorescence phenotype, numerous mutants were produced, and by pursuing a mapping-by-sequencing approach, the genomic segments containing mutated candidate genes were identified. One of the lines with an enhanced GFP fluorescence phenotype (named ELEVATED BUNDLE SHEATH CELLS SIGNAL 1 [ebss1]) was selected for further study, and the responsible gene was verified by CRISPR/Cas9-based mutagenesis of candidate genes located in the mapped genomic segment. The verified gene, At2g25970, encodes a K homology (KH) domain-containing protein.

Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD

Advances in affordable transcriptome sequencing combined with better exon and gene prediction has motivated many to compare transcription across the tree of life. We develop a mathematical framework to calculate complexity and compare transcript models. Structural features, i.e. intron retention (IR), donor/acceptor site variation, alternative exon cassettes, alternative 5'/3' UTRs, are compared and the distance between transcript models is calculated with nucleotide level precision. All metrics are implemented in a PyPi package, TranD and output can be used to summarize splicing patterns for a transcriptome (1GTF) and between transcriptomes (2GTF). TranD output enables quantitative comparisons between: annotations augmented by empirical RNA-seq data and the original transcript models; transcript model prediction tools for longread RNA-seq (e.g. FLAIR versus Isoseq3); alternate annotations for a species (e.g. RefSeq vs Ensembl); and between closely related species. In C. elegans, Z. mays, D. melanogaster, D. simulans and H. sapiens, alternative exons were observed more frequently in combination with an alternative donor/acceptor than alone. Transcript models in RefSeq and Ensembl are linked and both have unique transcript models with empirical support. D. melanogaster and D. simulans, share many transcript models and long-read RNAseq data suggests that both species are under-annotated. We recommend combined references.

Populus euphratica GRP2 Interacts with Target mRNAs to Negatively Regulate Salt Tolerance by Interfering with Photosynthesis, Na+, and ROS Homeostasis

The transcription of glycine-rich RNA-binding protein 2 (PeGRP2) transiently increased in the roots and shoots of Populus euphratica (a salt-resistant poplar) upon initial salt exposure and tended to decrease after long-term NaCl stress (100 mM, 12 days). PeGRP2 overexpression in the hybrid Populus tremula × P. alba '717-1B4' (P. × canescens) increased its salt sensitivity, which was reflected in the plant's growth and photosynthesis. PeGRP2 contains a conserved RNA recognition motif domain at the N-terminus, and RNA affinity purification (RAP) sequencing was developed to enrich the target mRNAs that physically interacted with PeGRP2 in P. × canescens. RAP sequencing combined with RT-qPCR revealed that NaCl decreased the transcripts of PeGRP2-interacting mRNAs encoding photosynthetic proteins, antioxidative enzymes, ATPases, and Na+/H+ antiporters in this transgenic poplar. Specifically, PeGRP2 negatively affected the stability of the target mRNAs encoding the photosynthetic proteins PETC and RBCMT; antioxidant enzymes SOD[Mn], CDSP32, and CYB1-2; ATPases AHA11, ACA8, and ACA9; and the Na+/H+ antiporter NHA1. This resulted in (i) a greater reduction in Fv/Fm, YII, ETR, and Pn; (ii) less pronounced activation of antioxidative enzymes; and (iii) a reduced ability to maintain Na+ homeostasis in the transgenic poplars during long-term salt stress, leading to their lowered ability to tolerate salinity stress.

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.

Genome-wide identification of RNA recognition motif (RRM1) in Brassica rapa and functional analysis of RNA-binding protein (BrRBP) under low-temperature stress

BACKGROUND

The RNA recognition motif (RRM) is primarily engaged in the processing of mRNA and rRNA following gene transcription as well as the regulation of RNA transport; it is critical in preserving RNA stability.

RESULTS

In this study, we identified 102 members of the RRM1 gene family in Brassica rapa, which were dispersed across 10 chromosomes with the ninth chromosome being the most extensively distributed. The RRM1 gene family members of Brassica rapa and Arabidopsis thaliana were grouped into 14 subclades (I-XIV) using phylogenetic analysis. Moreover, the results of transcriptome analysis and RT-qPCR indicated that the expression of Brapa05T000840 was upregulated in the cultivars 'Longyou 7' and 'Longyou 99' following exposure to cold stress at a temperature of 4 °C for 24 h. The levels of expression in the leaves and growth cones of the 'Longyou 7' variety were found to be significantly higher than those observed in the 'Longyou 99' variety under conditions of low temperature and NaCl stress. It illustrates the involvement of the RRM1 gene in the physiological response to both low temperature and salt stress. In addition, it was observed that the survival rate of transgenic BrRBP (Brapa05T000840) Arabidopsis thaliana plants was notably higher compared to that of wild-type plants when subjected to varying durations of low temperature treatment. Furthermore, the expression of the BrRBP gene in transgenic plants exhibited an upward trend as the duration of low temperature treatment increased, reaching its peak at 24 h. The in-vivo enzymatic activity of reactive oxygen species-scavenging enzymes were found to be significantly elevated in comparison to wild-type plants, suggesting that the BrRBP gene may enhance the cold tolerance of Arabidopsis thaliana.

CONCLUSIONS

This study offers a significant foundation for comprehending the regulation mechanism of the RRM1 gene family in winter Brassica rapa subjected to cold stress, as well as for finding key genes associated with cold resistance.

In silico identification, characterization, and expression analysis of RNA recognition motif (RRM) containing RNA-binding proteins in Aedes aegypti

RNA-binding proteins (RBPs) are the proteins that bind RNAs and regulate their functioning. RBPs in mosquitoes are gaining attention due to their ability to bind flaviviruses and regulate their replication and transmission. Despite their relevance, RBPs in mosquitoes are not explored much. In this study, we screened the whole genome of Aedes aegypti, the primary vector of several pathogenic viruses, and identified the proteins containing RNA recognition motif (RRM), the most abundant protein domain in eukaryotes. Using several in silico strategies, a total of 135 RRM-containing RBPs were identified in Ae. aegypti. The proteins were characterized based on their available annotations and the sequence similarity with Drosophila melanogaster. Ae. aegypti RRM-containing RBPs included serine/arginine-rich (SR) proteins, polyadenylate-binding proteins (PABP), heteronuclear ribonucleoproteins (hnRNP), small nuclear ribonucleoproteins (snRNP), splicing factors, eukaryotic initiation factors, transformers, and nucleolysins. Phylogenetic analysis revealed that the proteins and the domain organization are conserved among Ae. aegypti, Bombyx mori, and Drosophila melanogaster. However, the gene length and the intron-exon organization varied across the insect species. Expression analysis of the genes encoding RBPs using publicly available RNA sequencing data for different developmental time points of the mosquito life cycle starting from the ovary and eggs up to the adults revealed stage-specific expression with several genes preferentially expressed in early embryonic stages and blood-fed female ovaries. This is the first database for the Ae. aegypti RBPs that can serve as the reference base for future investigations. Stage-specific genes can be further explored to determine their role in mosquito growth and development with a focus on developing novel mosquito control strategies.

Cauliflower mosaic virus protein P6 is a multivalent node for RNA granule proteins and interferes with stress granule responses during plant infection

Biomolecular condensation is a multipurpose cellular process that viruses use ubiquitously during their multiplication. Cauliflower mosaic virus replication complexes are condensates that differ from those of most viruses, as they are nonmembranous assemblies that consist of RNA and protein, mainly the viral protein P6. Although these viral factories (VFs) were described half a century ago, with many observations that followed since, functional details of the condensation process and the properties and relevance of VFs have remained enigmatic. Here, we studied these issues in Arabidopsis thaliana and Nicotiana benthamiana. We observed a large dynamic mobility range of host proteins within VFs, while the viral matrix protein P6 is immobile, as it represents the central node of these condensates. We identified the stress granule (SG) nucleating factors G3BP7 and UBP1 family members as components of VFs. Similarly, as SG components localize to VFs during infection, ectopic P6 localizes to SGs and reduces their assembly after stress. Intriguingly, it appears that soluble rather than condensed P6 suppresses SG formation and mediates other essential P6 functions, suggesting that the increased condensation over the infection time-course may accompany a progressive shift in selected P6 functions. Together, this study highlights VFs as dynamic condensates and P6 as a complex modulator of SG responses.

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.

Arabidopsis SRPKII family proteins regulate flowering via phosphorylation of SR proteins and effects on gene expression and alternative splicing

Alternative splicing of pre-mRNAs is crucial for plant growth and development. Serine/arginine-rich (SR) proteins are a conserved family of RNA-binding proteins that are critical for both constitutive and alternative splicing. However, how phosphorylation of SR proteins regulates gene transcription and alternative splicing during plant development is poorly understood. We found that the Arabidopsis thaliana L. SR protein-specific kinase II family proteins (SRPKIIs) play an important role in plant development, including flowering. SRPKIIs regulate the phosphorylation status of a subset of specific SR proteins, including SR45 and SC35, which subsequently mediates their subcellular localization. A phospho-dead SR45 mutant inhibits the assembly of the apoptosis-and splicing-associated protein complex and thereby upregulates the expression of FLOWERING LOCUS C (FLC) via epigenetic modification. The splicing efficiency of FLC introns was significantly increased in the shoot apex of the srpkii mutant. Transcriptomic analysis revealed that SRPKIIs regulate the alternative splicing of c. 400 genes, which largely overlap with those regulated by SR45 and SC35-SCL family proteins. In summary, we found that Arabidopsis SRPKIIs specifically affect the phosphorylation status of a subset SR proteins and regulate the expression and alternative splicing of FLC to control flowering time.

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.

The flowering time regulator FLK controls pathogen defense in Arabidopsis thaliana

Plant disease resistance is a complex process that is maintained in an intricate balance with development. Increasing evidence indicates the importance of posttranscriptional regulation of plant defense by RNA binding proteins. In a genetic screen for suppressors of Arabidopsis (Arabidopsis thaliana) accelerated cell death 6-1 (acd6-1), a small constitutive defense mutant whose defense level is grossly in a reverse proportion to plant size, we identified an allele of the canonical flowering regulatory gene FLOWERING LOCUS K HOMOLOGY DOMAIN (FLK) encoding a putative protein with triple K homology (KH) repeats. The KH repeat is an ancient RNA binding motif found in proteins from diverse organisms. The relevance of KH-domain proteins in pathogen resistance is largely unexplored. In addition to late flowering, the flk mutants exhibited decreased resistance to the bacterial pathogen Pseudomonas syringae and increased resistance to the necrotrophic fungal pathogen Botrytis cinerea. We further found that the flk mutations compromised basal defense and defense signaling mediated by salicylic acid (SA). Mutant analysis revealed complex genetic interactions between FLK and several major SA pathway genes. RNA-seq data showed that FLK regulates expression abundance of some major defense- and development-related genes as well as alternative splicing of a number of genes. Among the genes affected by FLK is ACD6, whose transcripts had increased intron retentions influenced by the flk mutations. Thus, this study provides mechanistic support for flk suppression of acd6-1 and establishes that FLK is a multifunctional gene involved in regulating pathogen defense and development of plants.

Two environmental signal-driven RNA metabolic processes: Alternative splicing and translation

Plants live in fixed locations and have evolved adaptation mechanisms that integrate multiple responses to various environmental signals. Among the different components of these response pathways, receptors/sensors represent nodes that recognise environmental signals. Additionally, RNA metabolism plays an essential role in the regulation of gene expression and protein synthesis. With the development of RNA biotechnology, recent advances have been made in determining the roles of RNA metabolism in response to different environmental signals-especially the roles of alternative splicing and translation. In this review, we discuss recent progress in research on how the environmental adaptation mechanisms in plants are affected at the posttranscriptional level. These findings improve our understanding of the mechanism through which plants adapt to environmental changes by regulating the posttranscriptional level and are conducive for breeding stress-tolerant plants to cope with dynamic and rapidly changing environments.

The Phytophthora sojae effector PsFYVE1 modulates immunity-related gene expression by targeting host RZ-1A protein

Oomycete pathogens secrete numerous effectors to manipulate plant immunity and promote infection. However, relatively few effector types have been well characterized. In this study, members of an FYVE domain-containing protein family that are highly expanded in oomycetes were systematically identified, and one secreted protein, PsFYVE1, was selected for further study. PsFYVE1 enhanced Phytophthora capsici infection in Nicotiana benthamiana and was necessary for Phytophthora sojae virulence. The FYVE domain of PsFYVE1 had PI3P-binding activity that depended on four conserved amino acid residues. Furthermore, PsFYVE1 targeted RNA-binding proteins RZ-1A/1B/1C in N. benthamiana and soybean (Glycine max), and silencing of NbRZ-1A/1B/1C genes attenuated plant immunity. NbRZ-1A was associated with the spliceosome complex that included three important components, glycine-rich RNA-binding protein 7 (NbGRP7), glycine-rich RNA-binding protein 8 (NbGRP8), and a specific component of the U1 small nuclear ribonucleoprotein complex (NbU1-70K). Notably, PsFYVE1 disrupted NbRZ-1A-NbGRP7 interaction. RNA-seq and subsequent experimental analysis demonstrated that PsFYVE1 and NbRZ-1A not only modulated pre-mRNA alternative splicing (AS) of the necrotic spotted lesions 1 (NbNSL1) gene, but also co-regulated transcription of hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (NbHCT), ethylene insensitive 2 (NbEIN2), and sucrose synthase 4 (NbSUS4) genes, which participate in plant immunity. Collectively, these findings indicate that the FYVE domain-containing protein family includes potential uncharacterized effector types and also highlight that plant pathogen effectors can regulate plant immunity-related genes at both AS and transcription levels to promote disease.

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.

The rice blast fungus SR protein 1 regulates alternative splicing with unique mechanisms

Serine/arginine-rich (SR) proteins are well known as splicing factors in humans, model animals and plants. However, they are largely unknown in regulating pre-mRNA splicing of filamentous fungi. Here we report that the SR protein MoSrp1 enhances and suppresses alternative splicing in a model fungal plant pathogen Magnaporthe oryzae. Deletion of MoSRP1 caused multiple defects, including reduced virulence and thousands of aberrant alternative splicing events in mycelia, most of which were suppressed or enhanced intron splicing. A GUAG consensus bound by MoSrp1 was identified in more than 94% of the intron or/and proximate exons having the aberrant splicing. The dual functions of regulating alternative splicing of MoSrp1 were exemplified in enhancing and suppressing the consensus-mediated efficient splicing of the introns in MoATF1 and MoMTP1, respectively, which both were important for mycelial growth, conidiation, and virulence. Interestingly, MoSrp1 had a conserved sumoylation site that was essential to nuclear localization and enhancing GUAG binding. Further, we showed that MoSrp1 interacted with a splicing factor and two components of the exon-joining complex via its N-terminal RNA recognition domain, which was required to regulate mycelial growth, development and virulence. In contrast, the C-terminus was important only for virulence and stress responses but not for mycelial growth and development. In addition, only orthologues from Pezizomycotina species could completely rescue defects of the deletion mutants. This study reveals that the fungal conserved SR protein Srp1 regulates alternative splicing in a unique manner.

Dwarf and High Tillering1 represses rice tillering through mediating the splicing of D14 pre-mRNA

Strigolactones (SLs) constitute a class of plant hormones that regulate many aspects of plant development, including repressing tillering in rice (Oryza sativa). However, how SL pathways are regulated is still poorly understood. Here, we describe a rice mutant dwarf and high tillering1 (dht1), which exhibits pleiotropic phenotypes (such as dwarfism and increased tiller numbers) similar to those of mutants defective in SL signaling. We show that DHT1 encodes a monocotyledon-specific hnRNP-like protein that acts as a previously unrecognized intron splicing factor for many precursor mRNAs (pre-mRNAs), including for the SL receptor gene D14. We find that the dht1 (DHT1I232F) mutant protein is impaired in its stability and RNA binding activity, causing defective splicing of D14 pre-mRNA and reduced D14 expression, and consequently leading to the SL signaling-defective phenotypes. Overall, our findings deepen our understanding of the functional diversification of hnRNP-like proteins and establish a connection between posttranscriptional splicing and SL signaling in the regulation of plant development.

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.

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.

Phase separation of HRLP regulates flowering time in Arabidopsis

RNA binding proteins mediate posttranscriptional RNA metabolism and play regulatory roles in many developmental processes in eukaryotes. Despite their known effects on the floral transition from vegetative to reproductive growth in plants, the underlying mechanisms remain largely obscure. Here, we show that a hitherto unknown RNA binding protein, hnRNP R-LIKE PROTEIN (HRLP), inhibits cotranscriptional splicing of a key floral repressor gene FLOWERING LOCUS C (FLC). This, in turn, facilitates R-loop formation near FLC intron I to repress its transcription, thereby promoting the floral transition in Arabidopsis thaliana. HRLP, together with the splicing factor ARGININE/SERINE-RICH 45, forms phase-separated nuclear condensates with liquid-like properties, which is essential for HRLP function in regulating FLC splicing, R-loop formation, and RNA Polymerase II recruitment. Our findings reveal that inhibition of cotranscriptional splicing of FLC by nuclear HRLP condensates constitutes the molecular basis for down-regulation of FLC transcript levels to ensure the reproductive success of Arabidopsis.

Genome-Wide Characterization and Expression Analysis of KH Family Genes Response to ABA and SA in Arabidopsis thaliana

K-homologous (KH) family is a type of nucleic acid-binding protein containing the KH domain and has been found to affect splicing and transcriptional regulation. However, KH family genes haven’t been investigated in plant species systematically. In this study, we identified 30 genes that belonged to the KH family based on HMM of the KH domain in Arabidopsis thaliana. Phylogenetic tree analysis showed that the KH family is grouped into three subgroups. Synteny analysis showed that AtKH9 and AtKH29 have the conserved synteny relationship between A. thaliana and the other five species. The AtKH9 and AtKH29 were located in the cytoplasm and nucleus. The seed germination rates of the mutants atkh9 and atkh29 were higher than wild-type after abscisic acid (ABA) and salicylic acid (SA) treatments. In addition, the expression of ABA-related genes, such as ABRE-binding factor 2 (ABF2), ABRE-binding factor 4 (ABF4), and delta 1-pyrroline-5-carboxylate synthase (P5CS), and an SA-related gene pathogenesis-related proteins b (PR1b) were downregulated after ABA and SA treatments, respectively. These results suggested that atkh9 and atkh29 mutants inhibit the effect of ABA and SA on seed germination. In conclusion, our results provide valuable information for further exploration of the function of KH family genes and propose directions and ideas for the identification and characterization of KH family genes in other plants.

Overexpression of cassava RSZ21b enhances drought tolerance in Arabidopsis

Drought is one of the major environmental constraints affecting crop productivity. Plants have to adjust their developmental and physiological processes to cope with drought. We previously identified 18 cassava serine/arginine-rich (SR) proteins that had a pivotal role in alternative splicing in response to environmental stress. However, functional characterization of SR proteins is rarely explored. Here, we characterized the RSZ subfamily gene MeRSZ21b in cassava. The RSZ21b belongs to the RSZ subfamily, which was widely distributed in major crops and was highly conserved. Quantitative RT-PCR assay showed that the expression of MeRSZ21b was significantly induced by drought. Moreover, overexpression of MeRSZ21b in Arabidopsis was hypersensitive to abscisic acid (ABA) in the phases of seed germination and post-germination seedling growth. Meantime, MeRSZ21b overexpression lines were resistant to sorbitol treatment, and quickly closed the stomata when compared with Col-0 under drought condition. Importantly, overexpression of MeRSZ21b resulted in improved drought tolerance through modulating ABA-dependent signaling. Therefore, our findings refine our knowledge of the SR protein-coding genes and provide novel insights for enhancing plant resistance to environmental stress.

The RNA recognition motif-containing protein UBA2c prevents early flowering by promoting transcription of the flowering repressor FLM in Arabidopsis

FLOWERING LOCUS M (FLM) is a well-known MADS-box transcription factor that is required for preventing early flowering under low temperatures in Arabidopsis thaliana. Alternative splicing of FLM is involved in the regulation of temperature-responsive flowering. However, how the basic transcript level of FLM is regulated is largely unknown. Here, we conducted forward genetic screening and identified a previously uncharacterized flowering repressor gene, UBA2c. Genetic analyses indicated that UBA2c represses flowering at least by promoting FLM transcription. We further demonstrated that UBA2c directly binds to FLM chromatin and facilitates FLM transcription by inhibiting histone H3K27 trimethylation, a histone marker related to transcriptional repression. UBA2c encodes a protein containing two putative RNA recognition motifs (RRMs) and one prion-like domain (PrLD). We found that UBA2c forms speckles in the nucleus and that both the RRMs and PrLD are required not only for forming the nuclear speckles but also for the biological function of UBA2c. These results identify a previously unknown flowering repressor and provide insights into the regulation of flowering time.

The mRNA-binding proteome of a critical phase transition during Arabidopsis seed germination

Arabidopsis thaliana seed germination is marked by extensive translational control at two critical phase transitions. The first transition refers to the start of hydration, the hydration translational shift. The second shift, the germination translational shift (GTS) is the phase between testa rupture and radicle protrusion at which the seed makes the all or nothing decision to germinate. The mechanism behind the translational regulation at these phase transitions is unknown. RNA binding proteins (RBPs) are versatile players in the post-transcriptional control of messenger RNAs (mRNAs) and as such candidates for regulating translation during seed germination. Here, we report the mRNA binding protein repertoire of seeds during the GTS. Thirty seed specific RBPs and 22 dynamic RBPs were identified during the GTS, like the putative RBP Vacuolar ATPase subunit A and RBP HSP101. Several stress granule markers were identified in this study, which suggests that seeds are prepared to quickly adapt the translation of specific mRNAs in response to changes in environmental conditions during the GTS. Taken together this study provides a detailed insight into the world of RBPs during seed germination and their possible regulatory role during this developmentally regulated process.

RNA-Binding Proteins: The Key Modulator in Stress Granule Formation and Abiotic Stress Response

To cope with abiotic environmental stress, plants rapidly change their gene expression transcriptionally and post-transcriptionally, the latter by translational suppression of selected proteins and the assembly of cytoplasmic stress granules (SGs) that sequester mRNA transcripts. RNA-binding proteins (RBPs) are the major players in these post-transcriptional processes, which control RNA processing in the nucleus, their export from the nucleus, and overall RNA metabolism in the cytoplasm. Because of their diverse modular domain structures, various RBP types dynamically co-assemble with their targeted RNAs and interacting proteins to form SGs, a process that finely regulates stress-responsive gene expression. This review summarizes recent findings on the involvement of RBPs in adapting plants to various abiotic stresses via modulation of specific gene expression events and SG formation. The relationship of these processes with the stress hormone abscisic acid (ABA) is discussed.

The Role of Pumilio RNA Binding Protein in Plants

Eukaryotic organisms have a posttranscriptional/translational regulation system for the control of translational efficiency. RNA binding proteins (RBPs) have been known to control target genes. One type of protein, Pumilio (Pum)/Puf family RNA binding proteins, show a specific binding of 3′ untranslational region (3′ UTR) of target mRNA and function as a post-transcriptional/translational regulator in eukaryotic cells. Plant Pum protein is involved in development and biotic/abiotic stresses. Interestingly, Arabidopsis Pum can control target genes in a sequence-specific manner and rRNA processing in a sequence-nonspecific manner. As shown in in silico Pum gene expression analysis, Arabidopsis and rice Pum genes are responsive to biotic/abiotic stresses. Plant Pum can commonly contribute to host gene regulation at the post-transcriptional/translational step, as can mammalian Pum. However, the function of plant Pum proteins is not yet fully known. In this review, we briefly summarize the function of plant Pum in defense, development, and environmental responses via recent research and bioinformatics data.

Multi-resBind: a residual network-based multi-label classifier for in vivo RNA binding prediction and preference visualization

BACKGROUND

Protein-RNA interactions play key roles in many processes regulating gene expression. To understand the underlying binding preference, ultraviolet cross-linking and immunoprecipitation (CLIP)-based methods have been used to identify the binding sites for hundreds of RNA-binding proteins (RBPs) in vivo. Using these large-scale experimental data to infer RNA binding preference and predict missing binding sites has become a great challenge. Some existing deep-learning models have demonstrated high prediction accuracy for individual RBPs. However, it remains difficult to avoid significant bias due to the experimental protocol. The DeepRiPe method was recently developed to solve this problem via introducing multi-task or multi-label learning into this field. However, this method has not reached an ideal level of prediction power due to the weak neural network architecture.

RESULTS

Compared to the DeepRiPe approach, our Multi-resBind method demonstrated substantial improvements using the same large-scale PAR-CLIP dataset with respect to an increase in the area under the receiver operating characteristic curve and average precision. We conducted extensive experiments to evaluate the impact of various types of input data on the final prediction accuracy. The same approach was used to evaluate the effect of loss functions. Finally, a modified integrated gradient was employed to generate attribution maps. The patterns disentangled from relative contributions according to context offer biological insights into the underlying mechanism of protein-RNA interactions.

CONCLUSIONS

Here, we propose Multi-resBind as a new multi-label deep-learning approach to infer protein-RNA binding preferences and predict novel interactions. The results clearly demonstrate that Multi-resBind is a promising tool to predict unknown binding sites in vivo and gain biology insights into why the neural network makes a given prediction.

Silencing of StRIK in potato suggests a role in periderm related to RNA processing and stress

BACKGROUND

The periderm is a protective barrier crucial for land plant survival, but little is known about genetic factors involved in its development and regulation. Using a transcriptomic approach in the cork oak (Q. suber) periderm, we previously identified an RS2-INTERACTING KH PROTEIN (RIK) homologue of unknown function containing a K homology (KH)-domain RNA-binding protein, as a regulatory candidate gene in the periderm.

RESULTS

To gain insight into the function of RIK in the periderm, potato (S. tuberosum) tuber periderm was used as a model: the full-length coding sequence of RIK, hereafter referred to as StRIK, was isolated, the transcript profile analyzed and gene silencing in potato performed to analyze the silencing effects on periderm anatomy and transcriptome. The StRIK transcript accumulated in all vegetative tissues studied, including periderm and other suberized tissues such as root and also in wounded tissues. Downregulation of StRIK in potato by RNA interference (StRIK-RNAi) did not show any obvious effects on tuber periderm anatomy but, unlike Wild type, transgenic plants flowered. Global transcript profiling of the StRIK-RNAi periderm did show altered expression of genes associated with RNA metabolism, stress and signaling, mirroring the biological processes found enriched within the in silico co-expression network of the Arabidopsis orthologue.

CONCLUSIONS

The ubiquitous expression of StRIK transcript, the flower associated phenotype and the differential expression of StRIK-RNAi periderm point out to a general regulatory role of StRIK in diverse plant developmental processes. The transcriptome analysis suggests that StRIK might play roles in RNA maturation and stress response in the periderm.

The Brachypodium distachyon cold-acclimated plasma membrane proteome is primed for stress resistance

In order to survive subzero temperatures, some plants undergo cold acclimation (CA) where low, nonfreezing temperatures, and/or shortened day lengths allow cold-hardening and survival during subsequent freeze events. Central to this response is the plasma membrane (PM), where low temperature is perceived and cellular homeostasis must be preserved by maintaining membrane integrity. Here, we present the first PM proteome of cold-acclimated Brachypodium distachyon, a model species for the study of monocot crops. A time-course experiment investigated CA-induced changes in the proteome following two-phase partitioning PM enrichment and label-free quantification by nano-liquid chromatography-mass spectrophotometry. Two days of CA were sufficient for membrane protection as well as an initial increase in sugar levels and coincided with a significant change in the abundance of 154 proteins. Prolonged CA resulted in further increases in soluble sugars and abundance changes in more than 680 proteins, suggesting both a necessary early response to low-temperature treatment, as well as a sustained CA response elicited over several days. A meta-analysis revealed that the identified PM proteins have known roles in low-temperature tolerance, metabolism, transport, and pathogen defense as well as drought, osmotic stress, and salt resistance suggesting crosstalk between stress responses, such that CA may prime plants for other abiotic and biotic stresses. The PM proteins identified here present keys to an understanding of cold tolerance in monocot crops and the hope of addressing economic losses associated with modern climate-mediated increases in frost events.

PRPI-SC: an ensemble deep learning model for predicting plant lncRNA-protein interactions

BACKGROUND

Plant long non-coding RNAs (lncRNAs) play vital roles in many biological processes mainly through interactions with RNA-binding protein (RBP). To understand the function of lncRNAs, a fundamental method is to identify which types of proteins interact with the lncRNAs. However, the models or rules of interactions are a major challenge when calculating and estimating the types of RBP.

RESULTS

In this study, we propose an ensemble deep learning model to predict plant lncRNA-protein interactions using stacked denoising autoencoder and convolutional neural network based on sequence and structural information, named PRPI-SC. PRPI-SC predicts interactions between lncRNAs and proteins based on the k-mer features of RNAs and proteins. Experiments proved good results on Arabidopsis thaliana and Zea mays datasets (ATH948 and ZEA22133). The accuracy rates of ATH948 and ZEA22133 datasets were 88.9% and 82.6%, respectively. PRPI-SC also performed well on some public RNA protein interaction datasets.

CONCLUSIONS

PRPI-SC accurately predicts the interaction between plant lncRNA and protein, which plays a guiding role in studying the function and expression of plant lncRNA. At the same time, PRPI-SC has a strong generalization ability and good prediction effect for non-plant data.

Plant Stress Granules: Trends and Beyond

Stress granules (SGs) are dynamic membrane-less condensates transiently assembled through liquid-liquid phase separation (LLPS) in response to stress. SGs display a biphasic architecture constituted of core and shell phases. The core is a conserved SG fraction fundamental for its assembly and consists primarily of proteins with intrinsically disordered regions and RNA-binding domains, along with translational-related proteins. The shell fraction contains specific SG components that differ among species, cell type, and developmental stage and might include metabolic enzymes, receptors, transcription factors, untranslated mRNAs, and small molecules. SGs assembly positively correlates with stalled translation associated with stress responses playing a pivotal role during the adaptive cellular response, post-stress recovery, signaling, and metabolic rewire. After stress, SG disassembly releases mRNA and proteins to the cytoplasm to reactivate translation and reassume cell growth and development. However, under severe stress conditions or aberrant cellular behavior, SG dynamics are severely disturbed, affecting cellular homeostasis and leading to cell death in the most critical cases. The majority of research on SGs has focused on yeast and mammals as model organism. Nevertheless, the study of plant SGs has attracted attention in the last few years. Genetics studies and adapted techniques from other non-plant models, such as affinity capture coupled with multi-omics analyses, have enriched our understanding of SG composition in plants. Despite these efforts, the investigation of plant SGs is still an emerging field in plant biology research. In this review, we compile and discuss the accumulated progress of plant SGs regarding their composition, organization, dynamics, regulation, and their relation to other cytoplasmic foci. Lastly, we will explore the possible connections among the most exciting findings of SGs from mammalian, yeast, and plants, which might help provide a complete view of the biology of plant SGs in the future.

Mapping of Functional Subdomains in the atALKBH9B m6A-Demethylase Required for Its Binding to the Viral RNA and to the Coat Protein of Alfalfa Mosaic Virus

N ⁶-methyladenosine (m⁶A) modification is a dynamically regulated RNA modification that impacts many cellular processes and pathways. This epitranscriptomic methylation relies on the participation of RNA methyltransferases (referred to as "writers") and demethylases (referred to as "erasers"), respectively. We previously demonstrated that the Arabidopsis thaliana protein atALKBH9B showed m⁶A-demethylase activity and interacted with the coat protein (CP) of alfalfa mosaic virus (AMV), causing a profound impact on the viral infection cycle. To dissect the functional activity of atALKBH9B in AMV infection, we performed a protein-mapping analysis to identify the putative domains required for regulating this process. In this context, the mutational analysis of the protein revealed that the residues between 427 and 467 positions are critical for in vitro binding to the AMV RNA. The atALKBH9B amino acid sequence showed intrinsically disordered regions (IDRs) located at the N-terminal part delimiting the internal AlkB-like domain and at the C-terminal part. We identified an RNA binding domain containing an RGxxxRGG motif that overlaps with the C-terminal IDR. Moreover, bimolecular fluorescent experiments allowed us to determine that residues located between 387 and 427 are critical for the interaction with the AMV CP, which should be critical for modulating the viral infection process. Finally, we observed that atALKBH9B deletions of either N-terminal 20 residues or the C-terminal's last 40 amino acids impede their accumulation in siRNA bodies. The involvement of the regions responsible for RNA and viral CP binding and those required for its localization in stress granules in the viral cycle is discussed.

Plant RNA Binding Proteins as Critical Modulators in Drought, High Salinity, Heat, and Cold Stress Responses: An Updated Overview

Plant abiotic stress responses are tightly regulated by different players at multiple levels. At transcriptional or post-transcriptional levels, several RNA binding proteins (RBPs) regulate stress response genes through RNA metabolism. They are increasingly recognized as critical modulators of a myriad of biological processes, including stress responses. Plant RBPs are heterogeneous with one or more conservative RNA motifs that constitute canonical/novel RNA binding domains (RBDs), which can bind to target RNAs to determine their regulation as per the plant requirements at given environmental conditions. Given its biological significance and possible consideration as a potential tool in genetic manipulation programs to improve key agronomic traits amidst frequent episodes of climate anomalies, studies concerning the identification and functional characterization of RBP candidate genes are steadily mounting. This paper presents a comprehensive overview of canonical and novel RBPs and their functions in major abiotic stresses including drought, heat, salt, and cold stress conditions. To some extent, we also briefly describe the basic motif structure of RBPs that would be useful in forthcoming studies. Additionally, we also collected RBP genes that were modulated by stress, but that lacked functional characterization, providing an impetus to conduct further research.

Identification of rye B chromosome-associated peptides by mass spectrometry

B chromosomes (Bs) are supernumerary dispensable components of the standard genome (A chromosomes, As) that have been found in many eukaryotes. So far, it is unkown whether the B-derived transcripts translate to proteins or if the host proteome is changed due to the presence of Bs. Comparative mass spectrometry was performed using the protein samples isolated from shoots of rye plants with and without Bs. We aimed to identify B-associated peptides and analyzed the effects of Bs on the total proteome. Our comparative proteome analysis demonstrates that the presence of rye Bs affects the total proteome including different biological function processes. We found 319 of 16 776 quantified features in at least three out of five +B plants but not in 0B plants; 31 of 319 features were identified as B-associated peptide features. According to our data mining, one B-specific protein fragment showed similarity to a glycine-rich RNA binding protein which differed from its A-paralogue by two amino acid insertions. Our result represents a milestone in B chromosome research, because this is the first report to demonstrate the existence of Bs changing the proteome of the host.

An interactome analysis reveals that Arabidopsis thaliana GRDP2 interacts with proteins involved in post-transcriptional processes

The Arabidopsis thaliana glycine-rich domain protein 2 (AtGRDP2) gene encodes a protein of unknown function that is involved in plant growth and salt stress tolerance. The AtGRDP2 protein (787 aa, At4g37900) is constituted by three domains: a DUF1399 located at the N-terminus, a potential RNA Recognition Motif (RRM) in the central region, and a short glycine-rich domain at the C-terminus. Herein, we analyzed the subcellular localization of AtGRDP2 protein as a GFP translational fusion and found it was localized in the cytosol and the nucleus of tobacco leaf cells. Truncated versions of AtGRDP2 showed that the DUF1399 or the RRM domains were sufficient for nuclear localization. In addition, we performed a yeast two-hybrid split-ubiquitin assay (Y2H) to identify potential interactors for AtGRDP2 protein. The Y2H assay identified proteins associated with RNA binding functions such as PABN3 (At5g65260), EF-1α (At1g07920), and CL15 (At3g25920). Heterodimeric associations in planta between AtGRDP2 and its interactors were carried out by Bimolecular Fluorescence Complementation (BiFC) assays. The data revealed heterodimeric interactions between AtGRDP2 and PABN3 in the nucleus and AtGRDP2 with EF-1α in the cytosol, while AtGRDP2-CL15 associations occurred only in the chloroplasts. Finally, functional characterization of the protein-protein interaction regions revealed that both DUF1399 and RRM domains were key for heterodimerization with its interactors. The AtGRDP2 interaction with these proteins in different compartments suggests that this glycine-rich domain protein is involved in post-transcriptional processes.

OsCRP1, a Ribonucleoprotein Gene, Regulates Chloroplast mRNA Stability That Confers Drought and Cold Tolerance

Chloroplast ribonucleoproteins (cpRNPs) are nuclear-encoded and highly abundant proteins that are proposed to function in chloroplast RNA metabolism. However, the molecular mechanisms underlying the regulation of chloroplast RNAs involved in stress tolerance are poorly understood. Here, we demonstrate that CHLOROPLAST RNA-BINDING PROTEIN 1 (OsCRP1), a rice (Oryza sativa) cpRNP gene, is essential for stabilization of RNAs from the NAD(P)H dehydrogenase (NDH) complex, which in turn enhances drought and cold stress tolerance. An RNA-immunoprecipitation assay revealed that OsCRP1 is associated with a set of chloroplast RNAs. Transcript profiling indicated that the mRNA levels of genes from the NDH complex significantly increased in the OsCRP1 overexpressing compared to non-transgenic plants, whereas the pattern in OsCRP1 RNAi plants were opposite. Importantly, the OsCRP1 overexpressing plants showed a higher cyclic electron transport (CET) activity, which is essential for elevated levels of ATP for photosynthesis. Additionally, overexpression of OsCRP1 resulted in significantly enhanced drought and cold stress tolerance with higher ATP levels compared to wild type. Thus, our findings suggest that overexpression of OsCRP1 stabilizes a set of mRNAs from genes of the NDH complex involved in increasing CET activity and production of ATP, which consequently confers enhanced drought and cold tolerance.

The KH domain facilitates the substrate specificity and unwinding processivity of DDX43 helicase

The K-homology (KH) domain is a nucleic acid-binding domain present in many proteins. Recently, we found that the DEAD-box helicase DDX43 contains a KH domain in its N-terminus; however, its function remains unknown. Here, we purified recombinant DDX43 KH domain protein and found that it prefers binding ssDNA and ssRNA. Electrophoretic mobility shift assay and NMR revealed that the KH domain favors pyrimidines over purines. Mutational analysis showed that the GXXG loop in the KH domain is involved in pyrimidine binding. Moreover, we found that an alanine residue adjacent to the GXXG loop is critical for binding. Systematic evolution of ligands by exponential enrichment, chromatin immunoprecipitation-seq, and cross-linking immunoprecipitation-seq showed that the KH domain binds C-/T-rich DNA and U-rich RNA. Bioinformatics analysis suggested that the KH domain prefers to bind promoters. Using 15N-heteronuclear single quantum coherence NMR, the optimal binding sequence was identified as TTGT. Finally, we found that the full-length DDX43 helicase prefers DNA or RNA substrates with TTGT or UUGU single-stranded tails and that the KH domain is critically important for sequence specificity and unwinding processivity. Collectively, our results demonstrated that the KH domain facilitates the substrate specificity and processivity of the DDX43 helicase.

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 Arabidopsis KH-domain protein FLOWERING LOCUS Y delays flowering by upregulating FLOWERING LOCUS C family members

We identified FLY as a previously uncharacterized RNA-binding-family protein that controls flowering time by positively regulating the expression of FLC clade members. The ability of flowering plants to adjust the timing of the floral transition based on endogenous and environmental signals contributes to their adaptive success. In Arabidopsis thaliana, the MADS-domain protein FLOWERING LOCUS C (FLC) and the FLC clade members FLOWERING LOCUS M/MADS AFFECTING FLOWERING1 (FLM/MAF1), MAF2, MAF3, MAF4, and MAF5 form nuclear complexes that repress flowering under noninductive conditions. However, how FLM/MAF genes are regulated requires further study. Using a genetic strategy, we showed that the previously uncharacterized K-homology (KH) domain protein FLOWERING LOCUS Y (FLY) modulates flowering time. The fly-1 knockout mutant and FLY artificial microRNA knockdown line flowered earlier than the wild type under long- and short-day conditions. The knockout fly-1 allele, a SALK T-DNA insertion mutant, contains an ~ 110-kb genomic deletion induced by T-DNA integration. FLC clade members were downregulated in the fly-1 mutants and FLY artificial microRNA knockdown line, whereas the level of the FLC antisense transcript COOLAIR was similar to that of the wild type. Our results identify FLY as a regulator that affects flowering time through upregulation of FLC clade members.

Identification of the wheat C3H gene family and expression analysis of candidates associated with seed dormancy and germination

C3H zinc finger transcription factors play important roles in managing various biotic/abiotic stresses in Aarabidopsis, rice, and maize. The functions of these factors in wheat, however, remain largely unclear. We identified 88 TaC3H genes that were divided into four subfamilies in this analysis. Gene structure and conserved domain analyses indicate that most members of the same subfamily have similar structures. A total of 76 paralogous and 48 orthologous pairs were identified and Ka/Ks values were used to analyze replication relationships amongst wheat, rice, and Arabidopsis. Gene ontology (GO) annotation analysis showed that most TaC3H genes possessed molecular functions, while transcriptome results showed that the 88 TaC3H genes responded to water imbibition. Microarray data for 53 TaC3H genes were obtained and heat maps were generated; these results indicate that these genes are expressed in 13 wheat tissues. Subcellular localization prediction analysis indicates that most TaC3H genes are located in the nucleus. Promoter analysis indicates that most TaC3H genes contained cis-elements including ABRE, GARE-motif, and MBS, indicating that these can respond to various biotic/abiotic stresses. Transcriptome data and quantitative real-time PCR analysis of wheat cultivars with contrasting seed dormancy phenotypes show that five genes TaC3H4/-18/-37/-51/-72 were very likely involved in seed dormancy and germination. Exogenous ABA treatment further indicated that these five genes were responsive to ABA, suggesting that there may be a crosstalk between these genes and ABA signaling pathway in controlling seed dormancy and germination. These results provide a theoretical basis for subsequent studies on TaC3H gene function and also contribute to studies on the C3H gene in other species.

Natural variation in glycine-rich region of Brassica oleracea cold shock domain protein 5 (BoCSDP5) is associated with low temperature tolerance

BACKGROUND

Low temperature (LT) or cold stress is a major environmental stress that seriously affects plant growth and development, limiting crop productivity. Cold shock domain proteins (CSDPs), which are present in most living organism, are involved in RNA metabolisms influencing abiotic stress tolerance.

OBJECTIVE

The aims of this study are to identify target gene for LT-tolerance, like CSDPs, characterize genetics, and develop molecular marker distinguishing LT-tolerance in cabbage (Brassica oleracea var. capitata).

METHODS

Semi-quantitative RT-PCR or qRT-PCR was used in gene expression study. LT-tolerance was determined by electrolyte leakage and PCR with allelic specific primers.

RESULTS

Allelic variation was found in BoCSDP5 coding sequence (CDs) between LT-tolerant (BN106 and BN553) and -susceptible inbred lines (BN107 and BN554). LT-tolerant inbred lines contained variant type of BoCSDP5 (named as BoCSDP5v) which encodes extra CCHC zinc finger domain at C-terminus. Association of LT-tolerance with BoCSDP5v was confirmed by electrolyte leakage and segregation using genetic population derived from BN553 and BN554 cross. Allelic variation in BoCSDP5 gene does not influence the rate of gene expression, but produces different proteins with different number of CCHC zinc finger domains. LT-tolerance marker designed on the basis of polymorphism between BoCSDP5 and BoCSDP5v was confirmed with samples used in previous B. oleracea CIRCADIAN CLOCK ASSOCIATED 1 (BoCCA1) marker validation.

CONCLUSIONS

LT-tolerant allele (BoCSDP5v) is dominant and independent of CBF pathway, and sufficient to generate molecular markers to identify LT-tolerant cabbage when it is used in combination with another marker, like BoCCA1-derived one. Production and analysis of overexpressing plants of BoCSDP1, BoCSDP3, BoCSDP5 and BoCSDP5v will be required for elucidating the function of CCHC zinc finger domains in LT-tolerance.

Identification and expression analysis of chloroplast ribonucleoproteins (cpRNPs) in Arabidopsis and rice

Chloroplast ribonucleoproteins (cpRNPs) are implicated in splicing, editing, and stability control of chloroplast RNAs as well as in regulating development and stress tolerance. To facilitate a comprehensive understanding of their functions, we carried out a genome-wide identification, curation, and phylogenetic analysis of cpRNP genes in Oryza sativa (rice) and Arabidopsis thaliana (Arabidopsis). Ten cpRNP genes were identified in each of Arabidopsis and rice genomes based on the presence of two RRM (RNA-recognition motif) domains and an N-terminal chloroplast targeting signal peptide in the predicted proteins. These proteins are localized to chloroplasts. Gene expression analysis revealed that cpRNP genes have differential tissue expression patterns and some cpRNP genes are induced by abiotic stresses such as cold, heat, and drought. Taken together, our study provides a comprehensive annotation of the cpRNP gene family and their expression patterns in Arabidopsis and rice which will facilitate further studies on their roles in plant growth and stress responses.

SPLICING FACTOR1 Is Important in Chloroplast Development under Cold Stress

RNA SPLICING FACTOR1 (SF1) is responsible for recognizing the branch point site (BPS) sequence in introns and is critical for pre-mRNA splicing. In Arabidopsis (Arabidopsis thaliana), splicing factor1 (AtSF1) has been shown to retain the conserved function, but it is unexpected that null atsf1 mutants are viable. Here, we identified an allele of atsf1, named suppressor of thf1-4 (sot4), from suppressor screening for leaf variegation of thylakoid formation1 The sot4 mutant resulting from the G-to-R mutation at the highly conserved 198th amino acid residue within the functionally unknown domain exhibits leaf virescence associated with less accumulation of mature plastid ribosomal RNA, particularly under cold stress. Interestingly, the same point mutation in yeast Saccharomyces cerevisiae MUD synthetic-lethal 5p (SF1/Msl5p) also causes hypersensitivity to coldness and a low splicing activity for the introns with suboptimal BPS sequences. Transcriptomic profiling and reverse-transcription quantitative PCR analyses showed that expression of many genes were up- or downregulated in atsf1 via insufficient intron splicing. Our search for a BPS consensus from the retained introns in atsf1 transcriptomes, combined with RNA electrophoresis mobility shift assays, revealed that AtSF1 directly binds to the BPS consensus containing 5'-CU(U/A)AU-3'. Taken together, our data provide insight into a role for AtSF1 in regulating intron splicing efficiency, which helps plants acclimate to coldness.

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.

Emerging Roles of RNA-Binding Proteins in Seed Development and Performance

Seed development, dormancy, and germination are key physiological events that are not only important for seed generation, survival, and dispersal, but also contribute to agricultural production. RNA-binding proteins (RBPs) directly interact with target mRNAs and fine-tune mRNA metabolism by governing post-transcriptional regulation, including RNA processing, intron splicing, nuclear export, trafficking, stability/decay, and translational control. Recent studies have functionally characterized increasing numbers of diverse RBPs and shown that they participate in seed development and performance, providing significant insight into the role of RBP-mRNA interactions in seed processes. In this review, we discuss recent research progress on newly defined RBPs that have crucial roles in RNA metabolism and affect seed development, dormancy, and germination.

Transcriptome profiling of raspberry (Rubus idaeus Var. Amira) in response to infection by tomato ringspot virus (ToRSV)

Raspberry (Rubus sp.) is a berries fruit with an ongoing agricultural and commercial interest due to its high contents of flavonoids and nutrients beneficial for human health. The growing demand for raspberries is facing great challenges associated mainly with the dispersal of diseases, which produces a decrease in productivity and fruit quality. A broad range of genomic resources is available for other Rosaceae species; however, genomic resources for species of the Rubus genus are still limited. Here, we characterize the transcriptome of the Rubus idaeus (Var. Amira) in order to 1) provide clues in the transcriptional changes of R. idaeus against tomato ringspot virus (ToRSV); and 2) generate genomic resources for this economically important species. We generate more than 200 million sequencing reads from two mRNA samples of raspberry, infected and not infected by ToRSV, using Illumina technology. After de novo assembly, we obtained 68,853 predicted protein-coding sequences of which 71.3% and 61.3% were annotated using Gene Ontology and Pfam databases, respectively. Moreover, we find 2,340 genes with differential expression between raspberries infected and not infected by ToRSV. Analysis of these genes shows functional enrichments of the oxidation-reduction process, cell wall biogenesis, terpene synthase activity, and lyase activity. These genes could be involved in the raspberry immune response through the interaction of different metabolic pathways; however, this statement needs further investigations. Up-regulation of genes encoding terpene synthases, multicopper oxidases, laccases, and beta-glucosidases might suggest that these enzymes appear to be the predominant transcriptome immune response of R. idaeus against ToRSV. Furthermore, we identify thousands of molecular markers (i.e., SSRs and SNPs), increasing considerably the genomic resources currently available for raspberries. This study is the first report on investigating the transcriptional changes of R. idaeus against ToRSV.

Molecular Characterization, Expression Pattern and Function Analysis of Glycine-Rich Protein Genes Under Stresses in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Plant Glycine-rich proteins (GRP), a superfamily with a glycine-rich domain, play an important role in various stresses such as high or low temperature stress and drought stress. GRP genes have been studied in many plants, but seldom in Chinese cabbage (Brassica rapa L. ssp. pekinensis). In this study, a total of 64 GRP genes were identified in Chinese cabbage by homology comparative analysis. The physical and chemical characteristics predicted by ProtParam tool revealed that 62.5% of BrGRPs were alkaline, 53.1% were stable, and 79.7% were hydrophilic. Conserved domain analysis by MEME and TBtools showed that 64 BrGRPs contained 20 of the same conserved motifs, based on which BrGRPs were classified into five main classes and four subclasses in class IV to clarify their evolutionary relationship. Our results demonstrated that The BrGRP genes were located on ten chromosomes and in three different subgenomes of Chinese cabbage, and 43 pairs of orthologous GRP genes were found between Chinese cabbage and Arabidopsis. According to the transcriptome data, 64 BrGRP genes showed abnormal expression under high temperature stress, 52 under low temperature stress, 39 under drought stress, and 23 responses to soft rot. A large number of stress-related cis-acting elements, such as DRE, MYC, MYB, and ABRE were found in their promoter regions by PlantCare, which corresponded with differential expressions. Two BrGRP genes-w546 (Bra030284) and w1409 (Bra014000), both belonging to the subfamily Subclass IVa RBP-GRP (RNA binding protein-glycine rich protein), were up-regulated under 150 mmol⋅L-1 NaCl stress in Chinese cabbage. However, the overexpressed w546 gene could significantly inhibit seed germination, while w1409 significantly accelerated seed germination under 100 mmol⋅L-1 NaCl or 300 mmol⋅L-1 mannitol stresses. In short, most BrGRP genes showed abnormal expression under adversity stress, and some were involved in multiple stress responses, suggesting a potential capacity to resist multiple biotic and abiotic stresses, which is worthy of further study. Our study provides a systematic investigation of the molecular characteristics and expression patterns of BrGRP genes and promotes for further work on improving stress resistance of Chinese cabbage.

Transcript isoform sequencing reveals widespread promoter-proximal transcriptional termination in Arabidopsis

RNA polymerase II (RNAPII) transcription converts the DNA sequence of a single gene into multiple transcript isoforms that may carry alternative functions. Gene isoforms result from variable transcription start sites (TSSs) at the beginning and polyadenylation sites (PASs) at the end of transcripts. How alternative TSSs relate to variable PASs is poorly understood. Here, we identify both ends of RNA molecules in Arabidopsis thaliana by transcription isoform sequencing (TIF-seq) and report four transcript isoforms per expressed gene. While intragenic initiation represents a large source of regulated isoform diversity, we observe that ~14% of expressed genes generate relatively unstable short promoter-proximal RNAs (sppRNAs) from nascent transcript cleavage and polyadenylation shortly after initiation. The location of sppRNAs correlates with the position of promoter-proximal RNAPII stalling, indicating that large pools of promoter-stalled RNAPII may engage in transcriptional termination. We propose that promoter-proximal RNAPII stalling-linked to premature transcriptional termination may represent a checkpoint that governs plant gene expression.

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.