Plant serine/arginine-rich proteins: roles in precursor messenger RNA splicing, plant development, and stress responses

Alternative splicing responses to salt stress in Glycyrrhiza uralensis revealed by global profiling of transcriptome RNA-seq datasets

Excessive reactive oxygen species stress due to salinity poses a significant threat to the growth of Glycyrrhiza uralensis Fisch. To adapt to salt stress, G. uralensis engages in alternative splicing (AS) to generate a variety of proteins that help it withstand the effects of salt stress. While several studies have investigated the impact of alternative splicing on plants stress responses, the mechanisms by which AS interacts with transcriptional regulation to modulate the salt stress response in G. uralensis remain poorly understood. In this study, we utilized high-throughput RNA sequencing data to perform a comprehensive analysis of AS events at various time points in G. uralensis under salt stress, with exon skipping (SE) being the predominant AS type. KEGG enrichment analysis was performed on the different splicing genes (DSG), and pathways associated with AS were significantly enriched, including RNA transport, mRNA surveillance, and spliceosome. This indicated splicing regulation of genes, resulting in AS events under salt stress conditions. Moreover, plant response to salt stress pathways were also enriched, such as mitogen-activated protein kinase signaling pathway - plant, flavonoid biosynthesis, and oxidative phosphorylation. We focused on four differentially significant genes in the MAPK pathway by AS and qRT-PCR analysis. The alternative splicing type of MPK4 and SnRK2 was skipped exon (SE). ETR2 and RbohD were retained intron (RI) and alternative 5'splice site (A5SS), respectively. The expression levels of isoform1 of these four genes displayed different but significant increases in different tissue sites and salt stress treatment times. These findings suggest that MPK4, SnRK2, ETR2, and RbohD in G. uralensis activate the expression of isoform1, leading to the production of more isoform1 protein and thereby enhancing resistance to salt stress. These findings suggest that salt-responsive AS directly and indirectly governs G. uralensis salt response. Further investigations into AS function and mechanism during abiotic stresses may offer novel references for bolstering plant stress tolerance.

SR45a plays a key role in enhancing cotton resistance to Verticillium dahliae by alternative splicing of immunity genes

Alternative splicing (AS) of pre-mRNAs increases the diversity of transcriptome and proteome and plays fundamental roles in plant development and stress responses. However, the prevalent changes in AS events and the regulating mechanisms of plants in response to pathogens remain largely unknown. Here, we show that AS changes are an important mechanism conferring cotton immunity to Verticillium dahliae (Vd). GauSR45a, encoding a serine/arginine-rich RNA binding protein, was upregulated expression and underwent AS in response to Vd infection in Gossypium australe, a wild diploid cotton species highly resistant to Vd. Silencing GauSR45a substantially reduced the splicing ratio of Vd-induced immune-associated genes, including GauBAK1 (BRI1-associated kinase 1) and GauCERK1 (chitin elicitor receptor kinase 1). GauSR45a binds to the GAAGA motif that is commonly found in the pre-mRNA of genes essential for PTI, ETI, and defense. The binding between GauSR45a and the GAAGA motif in the pre-mRNA of BAK1 was enhanced by two splicing factors of GauU2AF35B and GauU1-70 K, thereby facilitating exon splicing; silencing either AtU2AF35B or AtU1-70 K decreased the resistance to Vd in transgenic GauSR45a Arabidopsis. Overexpressing the short splicing variant of BAK1GauBAK1.1 resulted in enhanced Verticillium wilt resistance rather than the long one GauBAK1.2. Vd-induced far more AS events were in G. barbadense (resistant tetraploid cotton) than those in G. hirsutum (susceptible tetraploid cotton) during Vd infection, indicating resistance divergence in immune responses at a genome-wide scale. We provided evidence showing a fundamental mechanism by which GauSR45a enhances cotton resistance to Vd through global regulation of AS of immunity genes.

γ-Aminobutyric acid (GABA) priming alleviates acid-aluminum toxicity to roots of creeping bentgrass via enhancements in antioxidant defense and organic metabolites remodeling

MAIN CONCLUSION

γ-Aminobutyric acid alleviates acid-aluminum toxicity to roots associated with enhanced antioxidant metabolism as well as accumulation and transportation of citric and malic acids. Aluminum (Al) toxicity has become the main limiting factor for crop growth and development in acidic soils and is further being aggravated worldwide due to continuous industrial pollution. The current study was designed to examine effects of GABA priming on alleviating acid-Al toxicity in terms of root growth, antioxidant defense, citrate and malate metabolisms, and extensive metabolites remodeling in roots under acidic conditions. Thirty-seven-day-old creeping bentgrass (Agrostis stolonifera) plants were used as test materials. Roots priming with or without 0.5 mM GABA for 3 days were cultivated in standard nutrient solution for 15 days as control or subjected to nutrient solution containing 5 mM AlCl3·6H2O for 15 days as acid-Al stress treatment. Roots were sampled for determinations of root characteristics, physiological and biochemical parameters, and metabolomics. GABA priming significantly alleviated acid-Al-induced root growth inhibition and oxidative damage, despite it promoted the accumulation of Al in roots. Analysis of metabolomics showed that GABA priming significantly increased accumulations of organic acids, amino acids, carbohydrates, and other metabolites in roots under acid-Al stress. In addition, GABA priming also significantly up-regulated key genes related to accumulation and transportation of malic and citric acids in roots under acid-Al stress. GABA-regulated metabolites participated in tricarboxylic acid cycle, GABA shunt, antioxidant defense system, and lipid metabolism, which played positive roles in reactive oxygen species scavenging, energy conversion, osmotic adjustment, and Al ion chelation in roots.

Regulation of plant resistance to salt stress by the SnRK1-dependent splicing factor SRRM1L

Most splicing factors are extensively phosphorylated but their physiological functions in plant salt resistance are still elusive. We found that phosphorylation by SnRK1 kinase is essential for SRRM1L nuclear speckle formation and its splicing factor activity in plant cells. In Arabidopsis, loss-of-function of SRRM1L leads to the occurrence of alternative pre-mRNA splicing events and compromises plant resistance to salt stress. In Arabidopsis srrm1l mutant line, we identified an intron-retention Nuclear factor Y subunit A 10 (NFYA10) mRNA variant by RNA-Seq and found phosphorylation-dependent RNA-binding of SRRM1L is indispensable for its alternative splicing activity. In the wild-type Arabidopsis, salt stress can activate SnRK1 to phosphorylate SRRM1L, triggering enrichment of functional NFYA10.1 variant to enhance plant salt resistance. By contrast, the Arabidopsis srrm1l mutant accumulates nonfunctional NFYA10.3 variant, sensitizing plants to salt stress. In summary, this work deciphered the molecular mechanisms and physiological functions of SnRK1-SRRM1L-NFYA10 module, shedding light on a regulatory pathway to fine-tune plant adaptation to abiotic stress at the post-transcriptional and post-translational levels.

SCR106 splicing factor modulates abiotic stress responses by maintaining RNA splicing in rice

Plants employ sophisticated molecular machinery to fine-tune their responses to growth, developmental, and stress cues. Gene expression influences plant cellular responses through regulatory processes such as transcription and splicing. Pre-mRNA is alternatively spliced to increase the genome coding potential and further regulate expression. Serine/arginine-rich (SR) proteins, a family of pre-mRNA splicing factors, recognize splicing cis-elements and regulate both constitutive and alternative splicing. Several studies have reported SR protein genes in the rice genome, subdivided into six subfamilies based on their domain structures. Here, we identified a new splicing factor in rice with an RNA recognition motif (RRM) and SR-dipeptides, which is related to the SR proteins, subfamily SC. OsSCR106 regulates pre-mRNA splicing under abiotic stress conditions. It localizes to the nuclear speckles, a major site for pre-mRNA splicing in the cell. The loss-of-function scr106 mutant is hypersensitive to salt, abscisic acid, and low-temperature stress, and harbors a developmental abnormality indicated by the shorter length of the shoot and root. The hypersensitivity to stress phenotype was rescued by complementation using OsSCR106 fused behind its endogenous promoter. Global gene expression and genome-wide splicing analysis in wild-type and scr106 seedlings revealed that OsSCR106 regulates its targets, presumably through regulating the alternative 3'-splice site. Under salt stress conditions, we identified multiple splice isoforms regulated by OsSCR106. Collectively, our results suggest that OsSCR106 is an important splicing factor that plays a crucial role in accurate pre-mRNA splicing and regulates abiotic stress responses in plants.

Nanopore Direct RNA Sequencing Reveals the Short-Term Salt Stress Response in Maize Roots

Transcriptome analysis, relying on the cutting-edge sequencing of cDNA libraries, has become increasingly prevalent within functional genome studies. However, the dependence on cDNA in most RNA sequencing technologies restricts their ability to detect RNA base modifications. To address this limitation, the latest Oxford Nanopore Direct RNA Sequencing (ONT DRS) technology was employed to investigate the transcriptome of maize seedling roots under salt stress. This approach aimed to unveil both the RNA transcriptional profiles and alterations in base modifications. The analysis of the differential expression revealed a total of 1398 genes and 2223 transcripts that exhibited significant variation within the maize root system following brief exposure to salt stress. Enrichment analyses, such as the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway assessments, highlighted the predominant involvement of these differentially expressed genes (DEGs) in regulating ion homeostasis, nitrogen metabolism, amino acid metabolism, and the phytohormone signaling pathways. The protein-protein interaction (PPI) analysis showed the participation of various proteins related to glycolytic metabolism, nitrogen metabolism, amino acid metabolism, abscisic acid signaling, and the jasmonate signaling pathways. It was through this intricate molecular network that these proteins collaborated to safeguard root cells against salt-induced damage. Moreover, under salt stress conditions, the occurrence of variable shear events (AS) in RNA modifications diminished, the average length of poly(A) tails underwent a slight decrease, and the number of genes at the majority of the variable polyadenylation (APA) sites decreased. Additionally, the levels of N5-methylcytosine (m5C) and N6-methyladenosine (m6A) showed a reduction. These results provide insights into the mechanisms of early salt tolerance in maize.

MOS4-associated complex contributes to proper splicing and suppression of ER stress under long-term heat stress in Arabidopsis

Plants are often exposed not only to short-term (S-) but also to long-term (L-)heat stress over several consecutive days. A few Arabidopsis mutants defective in L-heat tolerance have been identified, but the molecular mechanisms are less understood for this tolerance than for S-heat stress tolerance. To elucidate the mechanisms of the former, we used a forward genetic screen for sensitive to long-term heat (sloh) mutants and isolated sloh3 and sloh63. The mutants were hypersensitive to L- but not to S-heat stress, and sloh63 was also hypersensitive to salt stress. We identified the causal genes, SLOH3 and SLOH63, both of which encoded splicing-related components of the MOS4-associated complex (MAC). This complex is widely conserved in eukaryotes and has been suggested to interact with spliceosomes. Both genes were induced by L-heat stress in a time-dependent manner, and some abnormal splicing events were observed in both mutants under L-heat stress. In addition, endoplasmic reticulum (ER) stress and subsequent unfolded protein response occurred in both mutants under L-heat stress and were especially prominent in sloh63, suggesting that enhanced ER stress is due to the salt hypersensitivity of sloh63. Splicing inhibitor pladienolide B led to concentration-dependent disturbance of splicing, decreased L-heat tolerance, and enhanced ER stress. These findings suggest that maintenance of precise mRNA splicing under L-heat stress by the MAC is important for L-heat tolerance and suppressing ER stress in Arabidopsis.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

Comprehensive study of serine/arginine-rich (SR) gene family in rice: characterization, evolution and expression analysis

As important regulators of alternative splicing (AS) events, serine/arginine (SR)-rich proteins play indispensable roles in the growth and development of organisms. Until now, the study of SR genes has been lacking in plants. In the current study, we performed genome-wide analysis on the SR gene family in rice. A total of 24 OsSR genes were phylogenetically classified into seven groups, corresponding to seven subfamilies. The OsSR genes' structures, distribution of conserved domains, and protein tertiary structure of OsSR were conserved within each subfamily. The synteny analysis revealed that segmental duplication events were critical for the expansion of OsSR gene family. Moreover, interspecific synteny revealed the distribution of orthologous SR gene pairs between rice and Arabidopsis, sorghum, wheat, and maize. Among all OsSR genes, 14 genes exhibited NAGNAG acceptors, and only four OsSR genes had AS events on the NAGNAG acceptors. Furthermore, the distinct tissue-specific expression patterns of OsSR genes showed that these genes may function in different developmental stages in rice. The AS patterns on the same OsSR gene were variable among the root, stem, leaf, and grains at different filling stages, and some isoforms could only be detected in one or a few of tested tissues. Meanwhile, our results showed that the expression of some OsSR genes changed dramatically under ABA, GA, salt, drought, cold or heat treatment, which were related to the wide distribution of corresponding cis-elements in their promoter regions, suggesting their specific roles in stress and hormone response. This research facilitates our understanding of SR gene family in rice and provides clues for further exploration of the function of OsSR genes.

Evolutionary Study of Protein Short Tandem Repeats in Protein Families

Tandem repeats in proteins are patterns of residues repeated directly adjacent to each other. The evolution of these repeats can be assessed by using groups of homologous sequences, which can help pointing to events of unit duplication or deletion. High pressure in a protein family for variation of a given type of repeat might point to their function. Here, we propose the analysis of protein families to calculate protein short tandem repeats (pSTRs) in each protein sequence and assess their variability within the family in terms of number of units. To facilitate this analysis, we developed the pSTR tool, a method to analyze the evolution of protein short tandem repeats in a given protein family by pairwise comparisons between evolutionarily related protein sequences. We evaluated pSTR unit number variation in protein families of 12 complete metazoan proteomes. We hypothesize that families with more dynamic ensembles of repeats could reflect particular roles of these repeats in processes that require more adaptability.

Complexity and regulation of age-dependent alternative splicing in Brachypodium distachyon

Alternative splicing (AS) is a gene regulatory mechanism that generates multiple transcripts of the same gene precursor by the spliceosome complex, promoting messenger RNA complexity and proteome diversity. Although AS is extensively studied in response to environmental stresses, whether it mediates age-dependent development and how it is adjusted by growth transitions are largely unknown. Here, we comprehensively explored the AS landscape at different development stages in the grass model plant Brachypodium (Brachypodium distachyon). We identified abundant coding genes and non-coding transcripts subject to dynamic AS regulation during juvenile, adult, and reproductive transitions. Moreover, we revealed that SC35-LIKE SPLICING FACTOR 33 (SCL33), a serine/arginine-rich splicing factor in spliceosomes, plays a redundant and antagonistic role with its putative paralog, SCL33L, in regulating intron assembly across distinct developmental stages. In addition, we determined global AS variations in microRNA156 (miR156)-overproducing plants, in which growth transitions are delayed, and found that SPLs were regulated by miR156 in intron retention alteration in addition to mRNA clearance and translation inhibition manners. Finally, we demonstrated a complex regulatory process of age-dependent AS events in B. distachyon that were coincidently or separately regulated by miR156 and SCL33/SCL33L. These results illustrate a substantial -machinery of AS that mediates phase transitions in plants.

Hydropriming and Nutripriming of Bread Wheat Seeds Improved the Flour's Nutritional Value of the First Unprimed Offspring

Seed hydropriming or nutripriming has been used for wheat biofortification. Previously, the untreated S1 offspring of bread wheat S0 seeds hydro- and nutriprimed with FeSO₄.7H₂O and/or ZnSO₄.7H₂O showed improved yield relative to the offspring of untreated S0 seeds. We hypothesize that such improvement would have its origin in the higher quality of S1 seeds resulting from plants whose seeds were primed. In this work, we characterised biochemically the whole-wheat flour of unprimed S1 offspring whose S0 seeds were hydro- and nutriprimed with Fe and/or Zn and compared it to the offspring of untreated S0 seeds (control). We identified and quantified 16 free amino acids and five soluble sugars per offspring using high-performance liquid chromatography and the Association of Official Analytical Chemists (AOAC) methods. The most abundant amino acids were glutamic acid and glutamine, proline, and glycine, presenting their highest contents in the offspring of seeds nutriprimed with 8 ppm Zn (0.351 mmol∙g^-1), 8 ppm Fe + 8 ppm Zn (0.199 mmol∙g^-1), and (0.135 mmol∙g^-1), respectively. The highest contents of glucose (1.91 mg∙g^-1 sample), ash (24.90 g∙kg^-1 dry matter, DM), and crude protein (209.70 g∙kg^-1 DM) were presented by the offspring resulting from 4 ppm Fe + 4 ppm Zn, 8 ppm Zn, and 8 ppm Fe + 8 ppm Zn, respectively. The highest total starch content (630.10 g∙kg^-1 DM) was detected in the offspring of seeds soaked in 8 ppm Fe. The nutritional value of the flour of the S1 offspring resulting from nutripriming was significantly higher than the control. Overall, the novelty of our research is that seed priming can improve the quality of the wheat grain and flour, at least till the first offspring, without the need to repeat the presowing treatment. Beyond the study of subsequent generations, the unravelling of transgenerational mechanisms underlying the biochemical improvement of the offspring is approached.

Insights into sweet potato SR proteins: from evolution to species-specific expression and alternative splicing

SR proteins from sweet potato have conserved functional domains and similar gene structures as that of Arabidopsis and rice in general. However, expression patterns and alternative splicing regulations of SR genes from different species have changed under stresses. Novel alternative splicing regulations were found in sweet potato SR genes. Serine/arginine-rich (SR) proteins play important roles in plant development and stress response by regulating the pre-mRNA splicing process. However, SR proteins have not been identified so far from an important crop sweet potato. Through bioinformatics analysis, our study identified 24 SR proteins from sweet potato, with comprehensively analyzing of protein characteristics, gene structure, chromosome localization, and cis-acting elements in promotors. Salt, heat, and mimic drought stresses triggered extensive but different expressional regulations on sweet potato SR genes. Interestingly, heat stress caused the most active disturbances in both gene transcription and pre-mRNA alternative splicing (AS). Tissue and species-specific transcriptional and pre-mRNA AS regulations in response to stresses were found in sweet potato, in comparison with Arabidopsis and rice. Moreover, novel patterns of pre-mRNA alternative splicing were found in SR proteins from sweet potato. Our study provided an insight into similarities and differences of SR proteins in different plant species from gene sequences to gene structures and stress responses, indicating SR proteins may regulate their downstream genes differently between different species and tissues by varied transcriptional and pre-mRNA AS regulations.

Regulatory Network of Serine/Arginine-Rich (SR) Proteins: The Molecular Mechanism and Physiological Function in Plants

Serine/arginine-rich (SR) proteins are a type of splicing factor. They play significant roles in constitutive and alternative pre-mRNA splicing, and are involved in post-splicing activities, such as mRNA nuclear export, nonsense-mediated mRNA decay, mRNA translation, and miRNA biogenesis. In plants, SR proteins function under a complex regulatory network by protein–protein and RNA–protein interactions between SR proteins, other splicing factors, other proteins, or even RNAs. The regulatory networks of SR proteins are complex—they are regulated by the SR proteins themselves, they are phosphorylated and dephosphorylated through interactions with kinase, and they participate in signal transduction pathways, whereby signaling cascades can link the splicing machinery to the exterior environment. In a complex network, SR proteins are involved in plant growth and development, signal transduction, responses to abiotic and biotic stresses, and metabolism. Here, I review the current status of research on plant SR proteins, construct a model of SR proteins function, and ask many questions about SR proteins in plants.

An oligonucleotide/oligosaccharide-binding-fold protein enhances the alternative splicing event producing thylakoid membrane-bound ascorbate peroxidase in Nicotiana tabacum

The stromal and thylakoid membrane-bound ascorbate peroxidase isoforms are produced by the alternative splicing event of the 3'-terminal region of the APXII gene in spinach (Spinacia oleracea) and tobacco (Nicotiana tabacum), but not in Arabidopsis (Arabidopsis thaliana). However, all alternative splicing variants were detected in APXII gene-transformed Arabidopsis, indicating the occurrence of its regulatory mechanisms in Arabidopsis. The efficiency of this alternative splicing event in producing thylakoid membrane-bound ascorbate peroxidase mRNA is regulated by a splicing regulatory cis element, but trans splicing regulatory factor(s) for alternative splicing remain unclear. To identify this factor, we conducted a forward genetic screen using Arabidopsis in combination with a luciferase reporter system to evaluate the alternative splicing efficiency of thylakoid membrane-bound ascorbate peroxidase mRNA production. We isolated 9 mutant lines that showed low efficiency of the AS in producing thylakoid membrane-bound ascorbate peroxidase mRNA compared with that in the control plants. From one mutant [APXII alternative splicing inhibition (apsi1)], the causal gene responsible for the phenotype, AT5G38890 (oligonucleotide/oligosaccharide-binding-fold protein, APSI1), was identified. The levels of thylakoid membrane-bound ascorbate peroxidase mRNA from the transformed APXII gene decreased and increased in APSI1 knockout and APSI1-overexpressing plants, respectively. APSI1 was localized to the nucleus and specifically bound to the splicing regulatory cis element sequence. Tobacco plants that disrupted the closest homologs of APSI1 showed low levels of endogenous thylakoid membrane-bound ascorbate peroxidase mRNA. These results indicate that APSI1 is an enhancing component of the alternative splicing event of APXII.

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.

Analysis of Whole-Transcriptome RNA-Seq Data Reveals the Involvement of Alternative Splicing in the Drought Response of Glycyrrhiza uralensis

Alternative splicing (AS) is a post-transcriptional regulatory mechanism that increases protein diversity. There is growing evidence that AS plays an important role in regulating plant stress responses. However, the mechanism by which AS coordinates with transcriptional regulation to regulate the drought response in Glycyrrhiza uralensis remains unclear. In this study, we performed a genome-wide analysis of AS events in G. uralensis at different time points under drought stress using a high-throughput RNA sequencing approach. We detected 2,479 and 2,764 AS events in the aerial parts (AP) and underground parts (UP), respectively, of drought-stressed G. uralensis. Of these, last exon AS and exon skipping were the main types of AS. Overall, 2,653 genes undergoing significant AS regulation were identified from the AP and UP of G. uralensis exposed to drought for 2, 6, 12, and 24 h. Gene Ontology analyses indicated that AS plays an important role in the regulation of nitrogen and protein metabolism in the drought response of G. uralensis. Notably, the spliceosomal pathway and basal transcription factor pathway were significantly enriched with differentially spliced genes under drought stress. Genes related to splicing regulators in the AP and UP of G. uralensis responded to drought stress and underwent AS under drought conditions. In summary, our data suggest that drought-responsive AS directly and indirectly regulates the drought response of G. uralensis. Further in-depth studies on the functions and mechanisms of AS during abiotic stresses will provide new strategies for improving plant stress resistance.

OsSCL30 overexpression reduces the tolerance of rice seedlings to low temperature, drought and salt

Rice is one of the main food crops for the world population. Various abiotic stresses, such as low temperature, drought, and high salinity, affect rice during the entire growth period, determining its yield and quality, and even leading to plant death. In this study, by constructing overexpression vectors D-163 + 1300:OsSCL30 and D-163 + 1300-AcGFP:OsSCL30-GFP, the mechanism of action of OsSCL30 in various abiotic stresses was explored. Bioinformatics analysis showed that OsSCL30 was located on the chromosome 12 of rice Nipponbare, belonging to the plant-specific SCL subfamily of the SR protein family. The 1500 bp section upstream of the open reading frame start site contains stress-related cis-acting elements such as ABRE, MYC, and MYB. Under normal conditions, the expression of OsSCL30 was higher in leaves and leaf sheaths. The results of reverse transcription polymerase chain reaction showed that the expression of OsSCL30 decreased after low temperature, drought and salt treatment. In root cells OsSCL30 was localized in the nuclei. The results of the rice seedling tolerance and recovery tests showed that overexpression of OsSCL30 diminished the resistance to low temperature, drought and salt stresses in transgenic rice and resulted in larger accumulation of reactive oxygen species. This study is of great significance for exploring the response mechanisms of SR proteins under abiotic stresses.

Quantitative phosphoproteomics reveals the role of wild soybean GsSnRK1 as a metabolic regulator under drought and alkali stresses

Drought and alkali stresses cause detrimental effects on plant growth and development. SnRK1 protein kinases act as key energy and stress sensors by phosphorylation-mediated signaling in the regulation of plant defense reactions against adverse environments. To understand SnRK1-dependent phosphorylation events in signaling pathways triggered by abiotic factors, we employed quantitative phosphoproteomics to compare the global changes in phosphopeptides and phosphoproteins in 2kinm mutant Arabidopsis (SnRK1.1 T-DNA knockout and SnRK1.2 knockdown by β-estradiol-induced RNAi) complemented with wild soybean GsSnRK1(wt) or dominant negative mutant GsSnRK1(K49M) in response to drought and alkali stresses. Among 4014 phosphopeptides (representing 2380 phosphoproteins) identified in this study, we finalized 74 phosphopeptides (representing 61 phosphoproteins), and 75 phosphopeptides (representing 57 phosphoproteins) showing significant changes in phosphorylation levels under drought and alkali treatments respectively. Function enrichment and protein-protein interaction analyses indicated that the differentially-expressed phosphoproteins (DPs) under drought and alkali stresses were mainly involved in signaling transduction, stress response, carbohydrate and energy metabolism, transport and membrane trafficking, RNA splicing and processing, DNA binding and gene expression, and protein synthesis/folding/degradation. These results provide assistance to identify bona fide and novel SnRK1 phosphorylation substrates and shed new light on the biological functions of SnRK1 kinase in responses to abiotic stresses. SIGNIFICANCE: These results provide assistance to identify novel SnRK1 phosphorylation substrates and regulatory proteins, and shed new light on investigating the potential roles of reversible phosphorylation in plant responses to abiotic stresses.

Genome-Wide Profiling of Alternative Splicing and Gene Fusion during Rice Black-Streaked Dwarf Virus Stress in Maize (Zea mays L.)

Rice black-streaked dwarf virus (RBSDV) causes maize rough dwarf disease (MRDD), which is a viral disease that significantly affects maize yields worldwide. Plants tolerate stress through transcriptional reprogramming at the alternative splicing (AS), transcriptional, and fusion gene (FG) levels. However, it is unclear whether and how AS and FG interfere with transcriptional reprogramming in MRDD. In this study, we performed global profiling of AS and FG on maize response to RBSDV and compared it with transcriptional changes. There are approximately 1.43 to 2.25 AS events per gene in maize infected with RBSDV. GRMZM2G438622 was only detected in four AS modes (A3SS, A5SS, RI, and SE), whereas GRMZM2G059392 showed downregulated expression and four AS events. A total of 106 and 176 FGs were detected at two time points, respectively, including six differentially expressed genes and five differentially spliced genes. The gene GRMZM2G076798 was the only FG that occurred at two time points and was involved in two FG events. Among these, 104 GOs were enriched, indicating that nodulin-, disease resistance-, and chloroplastic-related genes respond to RBSDV stress in maize. These results provide new insights into the mechanisms underlying post-transcriptional and transcriptional regulation of maize response to RBSDV stress.

Genome-Wide Characterization of Serine/Arginine-Rich Gene Family and Its Genetic Effects on Agronomic Traits of Brassica napus

Serine/arginine-rich (SR) proteins are indispensable factors for RNA splicing, and they play important roles in development and abiotic stress responses. However, little information on SR genes in Brassica napus is available. In this study, 59 SR genes were identified and classified into seven subfamilies: SR, SCL, RS2Z, RSZ, RS, SR45, and SC. In each subfamily, the genes showed relatively conserved structures and motifs, but displayed distinct expression patterns in different tissues and under abiotic stress, which might be caused by the varied cis-acting regulatory elements among them. Transcriptome datasets from Pacbio/Illumina platforms showed that alternative splicing of SR genes was widespread in B. napus and the majority of paralogous gene pairs displayed different splicing patterns. Protein-protein interaction analysis indicated that SR proteins were involved in the regulation of the whole lifecycle of mRNA, from synthesis to decay. Moreover, the association mapping analysis suggested that 12 SR genes were candidate genes for regulating specific agronomic traits, which indicated that SR genes could affect the development and hence influence the important agronomic traits of B. napus. In summary, this study provided elaborate information on SR genes in B. napus, which will aid further functional studies and genetic improvement of agronomic traits in B. napus.

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 pattern of DNA methylation alteration, and its association with the changes of gene expression and alternative splicing during phosphate starvation in tomato

DNA methylation is changed and associates with gene expression alterations in plant response to phosphate starvation (Pi-), a common stress that impacts plant growth and productivity. However, in the horticultural model species Solanum lycopersicum (tomato), the dynamics of DNA methylation and its relationship with changes in gene transcription and alternative splicing (AS) under Pi- are unknown. Here, we performed integrative methylome and transcriptome analyses of tomato seedlings under Pi-deficient and -sufficient conditions. We found Pi- caused a slight increase in the overall methylation level, with millions of differentially methylated cytosines (DmCs) and a few hundred differentially methylated regions (DMRs). We also identified thousands of differentially expressed (DE) and differential AS (DAS) genes induced by Pi-, and found that DmCs were more abundant in non-expressed genes than in DE or DAS genes. Moreover, DNA methylation alterations weakly correlated with transcription changes but not with DAS events, and hyper-CHH-DMRs overlapping with transposable elements (TEs) were enriched in a subset of Pi starvation response (PSR) genes. We propose that changes in DNA methylation may be associated with the differential expression of some PSR genes, but that most of these changes probably control the expression of nearby TEs, rather than directly affecting the transcription or AS of PSR genes. Besides, the pattern of methylation changes upon Pi- may largely be shaped by TE distributions. Together, our study provides comprehensive insights into the association of DNA methylation with gene transcription and AS under Pi- in tomato and may contribute to unveiling novel roles of epigenetic mechanisms in plant stress response.

Overexpression of SCL30A from cassava (Manihot esculenta) negatively regulates salt tolerance in Arabidopsis

Soil salinity is a significant threat to sustainable agricultural production. Plants must adjust their developmental and physiological processes to deal with environmental salt conditions. We previously identified 18 serine-arginine-rich (SR) proteins from cassava (Manihot esculenta Crantz) that play pivotal roles in alternative splicing when encountering the external stress condition. However, functional characterisation of SR proteins is less reported in cassava, which is an important staple crop in the world. In the current study, we found that the expression of cassava spliceosomal component 35-like 30A (MeSCL30A) was significantly induced in response to drought and salt stress. The MeSCL30A overexpressing lines were also obtained in Arabidopsis thaliana L., which flowered earlier when compared with Col-0. Moreover, the MeSCL30A overexpressing lines were hypersensitive to salt and drought stress with lower germination and greening rate in comparison to Col-0. Importantly, soil-grown overexpression lines exhibited salt sensitivity through modulating the reactive oxygen species homeostasis and negatively regulating the gene expression that involved in ionic stress pathway. Therefore, these findings refined the SR protein-coding genes and provided novel insights for enhancing the resistance to environmental stress in plant.

The SR Splicing Factors: Providing Perspectives on Their Evolution, Expression, Alternative Splicing, and Function in Populus trichocarpa

Serine/arginine-rich (SR) proteins are important splicing factors in plant development and abiotic/hormone-related stresses. However, evidence that SR proteins contribute to the process in woody plants has been lacking. Using phylogenetics, gene synteny, transgenic experiments, and RNA-seq analysis, we identified 24 PtSR genes and explored their evolution, expression, and function in Popolus trichocarpa. The PtSR genes were divided into six subfamilies, generated by at least two events of genome triplication and duplication. Notably, they were constitutively expressed in roots, stems, and leaves, demonstrating their fundamental role in P. trichocarpa. Additionally, most PtSR genes (~83%) responded to at least one stress (cold, drought, salt, SA, MeJA, or ABA), and, especially, cold stress induced a dramatic perturbation in the expression and/or alternative splicing (AS) of 18 PtSR genes (~75%). Evidentially, the overexpression of PtSCL30 in Arabidopsis decreased freezing tolerance, which probably resulted from AS changes of the genes (e.g., ICE2 and COR15A) critical for cold tolerance. Moreover, the transgenic plants were salt-hypersensitive at the germination stage. These indicate that PtSCL30 may act as a negative regulator under cold and salt stress. Altogether, this study sheds light on the evolution, expression, and AS of PtSR genes, and the functional mechanisms of PtSCL30 in woody plants.

Genome-wide identification of cassava MeRboh genes and functional analysis in Arabidopsis

Plant respiratory burst oxidase homolog (Rboh) gene family encodes NADPH oxidases, and plays important roles in the production of reactive oxygen species (ROS), plant signaling, growth and stress responses. Cassava is an important starchy crops in tropical region. Environmental stresses, such as drought, pathogen, have caused great yield loss. The mechanisms of stress response are little known in MeRBOH family of cassava. Investigation of Rboh genes response to disease may provide a clue for clarification the disease resistance mechanisms. In this study, eight MeRboh genes were identified from the cassava genome. Comparisons of gene structure, protein motifs, and a phylogenetic tree showed conservation of Rboh gene families in cassava, Arabidopsis and rice. Transcript levels of most MeRboh genes increased following treatment with a pathogen, Xanthomonas axonopodis pv. manihotis, or with phytohormones salicylic acid or jasmonic acid. Analysis of cis-acting elements also indicated that MeRboh genes could response to light, hormone, abiotic and biotic stress. Prediction of miRNA target and post-translation modification sites of MeRboh suggested possible regulations of miRNA and protein phosphorylation; and transient expression of MeRboh in cassava protoplasts confirmed their localization on plasma membrane. Expression of MeRbohB, MeRbohF partially complemented PAMP responses in Arabidopsis rboh mutants, including the expression of PTI marker FRK1, ROS production, peroxide accumulation and callose deposition. It suggesting that MeRbohB and MeRbohF may participate in the PTI pathway and contributed to ROS production triggered by pathogens. Moreover, overexpression of MeRbohB and MeRbohF enhanced the resistance of Arabidopsis against Pseudomonas syringae pv. tomato DC3000. Together, these results suggest the evolutionary conservation of MeRboh gene family and their important role in the immune response and in regulating the plant disease resistance, providing a foundation for revealing molecular mechanisms of cassava disease resistance.

Interpreting machine learning models to investigate circadian regulation and facilitate exploration of clock function

The circadian clock is an internal molecular 24-h timer that is critical to life on Earth. We describe a series of artificial intelligence (AI)– and machine learning (ML)–based approaches that enable more cost-effective analysis and insight into circadian regulation and function. Throughout the manuscript, we illuminate what is inside the ML “black box” via explanation or interpretation of predictive ML models. Using this interpretation of our models, we derive biological insights into why a prediction was made, alongside accurate predictions. Most innovatively, we use only DNA sequence features for accurate circadian gene expression prediction. Using explainable AI, we define possible, responsible regulatory elements as we make these predictions; this critically requires no prior knowledge of regulatory elements.

The circadian clock is an important adaptation to life on Earth. Here, we use machine learning to predict complex, temporal, and circadian gene expression patterns in Arabidopsis. Most significantly, we classify circadian genes using DNA sequence features generated de novo from public, genomic resources, facilitating downstream application of our methods with no experimental work or prior knowledge needed. We use local model explanation that is transcript specific to rank DNA sequence features, providing a detailed profile of the potential circadian regulatory mechanisms for each transcript. Furthermore, we can discriminate the temporal phase of transcript expression using the local, explanation-derived, and ranked DNA sequence features, revealing hidden subclasses within the circadian class. Model interpretation/explanation provides the backbone of our methodological advances, giving insight into biological processes and experimental design. Next, we use model interpretation to optimize sampling strategies when we predict circadian transcripts using reduced numbers of transcriptomic timepoints. Finally, we predict the circadian time from a single, transcriptomic timepoint, deriving marker transcripts that are most impactful for accurate prediction; this could facilitate the identification of altered clock function from existing datasets.

Normal, novel or none: versatile regulation from alternative splicing

Pre-mRNA splicing is a vital step in the posttranscriptional regulation of gene expression. Splicing is catalyzed by the spliceosome, a multidalton RNA-protein complex, through two successive transesterifications to yield mature mRNAs. In Arabidopsis, more than 61% of all transcripts from intron-containing genes are alternatively spliced, thereby resulting in transcriptome and subsequent proteome diversities for cellular processes. Moreover, it is estimated that more alternative splicing (AS) events induced by adverse stimuli occur to confer stress tolerance. Recently, increasing AS variants encoding normal or novel proteins, or degraded by nonsense-mediated decay (NMD) and their corresponding splicing factors or regulators acting at the posttranscriptional level have been functionally characterized. This review comprehensively summarizes and highlights the advances in our understanding of the biological functions and underlying mechanisms of AS events and their regulators in Arabidopsis and provides prospects for further research on AS in crops.

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.

Daily temperature cycles promote alternative splicing of RNAs encoding SR45a, a splicing regulator in maize

Elevated temperatures enhance alternative RNA splicing in maize (Zea mays) with the potential to expand the repertoire of plant responses to heat stress. Alternative RNA splicing generates multiple RNA isoforms for many maize genes, and here we observed changes in the pattern of RNA isoforms with temperature changes. Increases in maximum daily temperature elevated the frequency of the major modes of alternative splices (AS), in particular retained introns and skipped exons. The genes most frequently targeted by increased AS with temperature encode factors involved in RNA processing and plant development. Genes encoding regulators of alternative RNA splicing were themselves among the principal AS targets in maize. Under controlled environmental conditions, daily changes in temperature comparable to field conditions altered the abundance of different RNA isoforms, including the RNAs encoding the splicing regulator SR45a, a member of the SR45 gene family. We established an "in protoplast" RNA splicing assay to show that during the afternoon on simulated hot summer days, SR45a RNA isoforms were produced with the potential to encode proteins efficient in splicing model substrates. With the RNA splicing assay, we also defined the exonic splicing enhancers that the splicing-efficient SR45a forms utilize to aid in the splicing of model substrates. Hence, with rising temperatures on hot summer days, SR45a RNA isoforms in maize are produced with the capability to encode proteins with greater RNA splicing potential.

Analysis of Whole Transcriptome RNA-seq Data Reveals Many Alternative Splicing Events in Soybean Roots under Drought Stress Conditions

Alternative splicing (AS) is a common post-transcriptional regulatory mechanism that modulates gene expression to increase proteome diversity. Increasing evidence indicates that AS plays an important role in regulating plant stress responses. However, the mechanism by which AS coordinates with transcriptional regulation to regulate drought responses in soybean remains poorly understood. In this study, we performed a genome-wide analysis of AS events in soybean (Glycine max) roots grown under various drought conditions using the high-throughput RNA-sequencing method, identifying 385, 989, 1429, and 465 AS events that were significantly differentially spliced under very mild drought stress, mild drought stress, severe drought stress, and recovery after severe drought conditions, respectively. Among them, alternative 3' splice sites and skipped exons were the major types of AS. Overall, 2120 genes that experienced significant AS regulation were identified from these drought-treated root samples. Gene Ontology term analysis indicated that the AS regulation of binding activity has vital roles in the drought response of soybean root. Notably, the genes encoding splicing regulatory factors in the spliceosome pathway and mRNA surveillance pathway were enriched according to the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Splicing regulatory factor-related genes in soybean root also responded to drought stress and were alternatively spliced under drought conditions. Taken together, our data suggest that drought-responsive AS acts as a direct or indirect mode to regulate drought response of soybean roots. With further in-depth research of the function and mechanism of AS in the process of abiotic stress, these results will provide a new strategy for enhancing stress tolerance of plants.

Alternative splicing and its regulatory role in woody plants

Alternative splicing (AS) is an important post-transcriptional process to enhance proteome diversity in eukaryotic organisms. In plants, numerous reports have primarily focused on AS analysis in model plant species or herbaceous plants, leading to a notable lack of research on AS in woody plants. More importantly, emerging evidence indicates that many important traits, including wood formation and stress resistance, in woody plants are controlled by AS. In this review article, we summarize the current progress of all kinds of AS studies in different tree species at various stages of development and in response to various stresses, revealing the significant role played by AS in woody plants, as well as the similar properties and differential regulation within their herbaceous counterparts. Furthermore, we propose several potential strategies to facilitate the functional characterization of splicing factors in woody plants and evaluate a general pipeline for the systematic characterization of splicing isoforms in a complex AS regulatory network. The utilization of genetic studies and high-throughput omics integration approaches to analyze AS genes and splicing factors is likely to further advance our understanding of AS modulation in woody plants.

Phytophthora Effectors Modulate Genome-wide Alternative Splicing of Host mRNAs to Reprogram Plant Immunity

Alternative splicing (AS) of pre-mRNAs increases transcriptome and proteome diversity, regulates gene expression through multiple mechanisms, and plays important roles in plant development and stress responses. However, the prevalence of genome-wide plant AS changes during infection and the mechanisms by which pathogens modulate AS remain poorly understood. Here, we examined the global AS changes in tomato leaves infected with Phytophthora infestans, the infamous Irish famine pathogen. We show that more than 2000 genes exhibiting significant changes in AS are not differentially expressed, indicating that AS is a distinct layer of transcriptome reprogramming during plant-pathogen interactions. Furthermore, our results show that P. infestans subverts host immunity by repressing the AS of positive regulators of plant immunity and promoting the AS of susceptibility factors. To study the underlying mechanism, we established a luminescence-based AS reporter system in Nicotiana benthamiana to screen pathogen effectors modulating plant AS. We identified nine splicing regulatory effectors (SREs) from 87 P. infestans effectors. Further studies revealed that SRE3 physically binds U1-70K to manipulate the plant AS machinery and subsequently modulates AS-mediated plant immunity. Our study not only unveils genome-wide plant AS reprogramming during infection but also establishes a novel AS screening tool to identify SREs from a wide range of plant pathogens, providing opportunities to understand the splicing regulatory mechanisms through which pathogens subvert plant immunity.

The interplay between microRNA and alternative splicing of linear and circular RNAs in eleven plant species

MOTIVATION

MicroRNA (miRNA) and alternative splicing (AS)-mediated post-transcriptional regulation has been extensively studied in most eukaryotes. However, the interplay between AS and miRNAs has not been explored in plants. To our knowledge, the overall profile of miRNA target sites in circular RNAs (circRNA) generated by alternative back splicing has never been reported previously. To address the challenge, we identified miRNA target sites located in alternatively spliced regions of the linear and circular splice isoforms using the up-to-date single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) and Illumina sequencing data in eleven plant species.

RESULTS

In total, we identified 399 401 and 114 574 AS events from linear and circular RNAs, respectively. Among them, there were 64 781 and 41 146 miRNA target sites located in linear and circular AS region, respectively. In addition, we found 38 913 circRNAs to be overlapping with 45 648 AS events of its own parent isoforms, suggesting circRNA regulation of AS of linear RNAs by forming R-loop with the genomic locus. Here, we present a comprehensive database of miRNA targets in alternatively spliced linear and circRNAs (ASmiR) and a web server for deposition and identification of miRNA target sites located in the alternatively spliced region of linear and circular RNAs. This database is accompanied by an easy-to-use web query interface for meaningful downstream analysis. Plant research community can submit user-defined datasets to the web service to search AS regions harboring small RNA target sites. In conclusion, this study provides an unprecedented resource to understand regulatory relationships between miRNAs and AS in both gymnosperms and angiosperms.

AVAILABILITY AND IMPLEMENTATION

The readily accessible database and web-based tools are available at http://forestry.fafu.edu.cn/bioinfor/db/ASmiR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Plant-derived smoke enhances plant growth through ornithine-synthesis pathway and ubiquitin-proteasome pathway in soybean

To investigate the mechanism of promotive effect of plant-derived smoke on the soybean growth, a gel-free/label-free proteomics was performed. Smoke solutions were irrigated on soybean or supplied simultaneously with flooding stress. Morphological and physiological analyses were performed for the confirmation of proteomic result. Metabolomic change was investigated to correlate proteomic change with metabolism regulation. Under normal condition, the length of root including hypocotyl increased in soybean treated with 2000 ppm plant-derived smoke within 4 days, as well as nitric oxide content. Proteins related to protein synthesis especially arginine metabolism were altered; metabolites related to amino acid, carboxylic acids, and sugars were mostly altered. Integrated analysis of omics data indicated that plant-derived smoke regulated nitrogen‑carbon transformation through ornithine synthesis pathway and promoted soybean normal growth. Under flooding, the number of lateral roots increased with root tip degradation in soybean treated with smoke solutions. Proteins related to ubiquitin-proteasome pathway were altered and led to sacrifice-for-survival-mechanism-driven degradation of root tip in soybean, which enabled accumulation of metabolites and guaranteed lateral root development during soybean recovery after flooding. These findings suggest that plant-derived smoke improves early stage of growth in soybean with regulation of ornithine-synthesis pathway and ubiquitin-proteasome pathway. BIOLOGICAL SIGNIFICANCE: Plant-derived smoke plays a key role in crop growth, however, the understanding of soybean in response to smoke treatment remains premature. Therefore, gel-free/label-free proteomic analysis was used for comprehensive study on the dual effect of smoke to soybean under normal and flooding conditions. Under normal condition, plant-derived smoke regulated nitrogen‑carbon transformation through ornithine synthesis pathway and resulted in the increase of the length of root including hypocotyl in soybean within 4 days. Under flooding condition, plant-derived smoke induced inhibition of ubiquitin-proteasome pathway and led to sacrifice-for-survival-mechanism-driven degradation of root tip in soybean, which enabled accumulation of metabolites and promoted lateral root development during soybean recovery after flooding.

Insights into the SAM Synthetase Gene Family and Its Roles in Tomato Seedlings under Abiotic Stresses and Hormone Treatments

S-Adenosyl-L-methionine (SAM) is a key enzyme involved in many important biological processes, such as ethylene and polyamine biosynthesis, transmethylation, and transsulfuration. Here, the SAM synthetase (SAMS) gene family was studied in ten different plants (Arabidopsis, tomato, eggplant, sunflower, Medicago truncatula, soybean, rice, barley, Triticum urartu and sorghum) with respect to its physical structure, physicochemical characteristics, and post-transcriptional and post-translational modifications. Additionally, the expression patterns of SAMS genes in tomato were analyzed based on a real-time quantitative PCR assay and an analysis of a public expression dataset. SAMS genes of monocots were more conserved according to the results of a phylogenetic analysis and the prediction of phosphorylation and glycosylation patterns. SAMS genes showed differential expression in response to abiotic stresses and exogenous hormone treatments. Solyc01g101060 was especially expressed in fruit and root tissues, while Solyc09g008280 was expressed in leaves. Additionally, our results revealed that exogenous BR and ABA treatments strongly reduced the expression of tomato SAMS genes. Our research provides new insights and clues about the role of SAMS genes. In particular, these results can inform future functional analyses aimed at revealing the molecular mechanisms underlying the functions of SAMS genes in plants.

Brassica Rapa SR45a Regulates Drought Tolerance via the Alternative Splicing of Target Genes

The emerging evidence has shown that plant serine/arginine-rich (SR) proteins play a crucial role in abiotic stress responses by regulating the alternative splicing (AS) of key genes. Recently, we have shown that drought stress enhances the expression of SR45a (also known as SR-like 3) in Brassica rapa. Herein, we unraveled the hitherto unknown functions of BrSR45a in drought stress response by comparing the phenotypes, chlorophyll a fluorescence and splicing patterns of the drought-responsive genes of Arabidopsis BrSR45a overexpressors (OEs), homozygous mutants (SALK_052345), and controls (Col-0). Overexpression and loss of function did not result in aberrant phenotypes; however, the overexpression of BrSR45a was positively correlated with drought tolerance and the stress recovery rate in an expression-dependent manner. Moreover, OEs showed a higher drought tolerance index during seed germination (38.16%) than the control lines. Additionally, the overexpression of BrSR45a induced the expression of the drought stress-inducible genes RD29A, NCED3, and DREB2A under normal conditions. To further illustrate the molecular linkages between BrSR45a and drought tolerance, we investigated the AS patterns of key drought-tolerance and BrSR45a interacting genes in OEs, mutants, and controls under both normal and drought conditions. The splicing patterns of DCP5, RD29A, GOLS1, AKR, U2AF, and SDR were different between overexpressors and mutants under normal conditions. Furthermore, drought stress altered the splicing patterns of NCED2, SQE, UPF1, U4/U6-U5 tri-snRNP-associated protein, and UPF1 between OEs and mutants, indicating that both overexpression and loss of function differently influenced the splicing patterns of target genes. This study revealed that BrSR45a regulates the drought stress response via the alternative splicing of target genes in a concentration-dependent manner.

The Features and Regulation of Co-transcriptional Splicing in Arabidopsis

Precursor mRNA (pre-mRNA) splicing is essential for gene expression in most eukaryotic organisms. Previous studies from mammals, Drosophila, and yeast show that the majority of splicing events occurs co-transcriptionally. In plants, however, the features of co-transcriptional splicing (CTS) and its regulation still remain largely unknown. Here, we used chromatin-bound RNA sequencing to study CTS in Arabidopsis thaliana. We found that CTS is widespread in Arabidopsis seedlings, with a large proportion of alternative splicing events determined co-transcriptionally. CTS efficiency correlated with gene expression level, the chromatin landscape and, most surprisingly, the number of introns and exons of individual genes, but is independent of gene length. In combination with enhanced crosslinking and immunoprecipitation sequencing analysis, we further showed that the hnRNP-like proteins RZ-1B and RZ-1C promote efficient CTS globally through direct binding, frequently to exonic sequences. Notably, this general effect of RZ-1B/1C on splicing promotion is mainly observed at the chromatin level, not at the mRNA level. RZ-1C promotes CTS of multiple-exon genes in association with its binding to regions both proximal and distal to the regulated introns. We propose that RZ-1C promotes efficient CTS of genes with multiple exons through cooperative interactions with many exons, introns, and splicing factors. Our work thus reveals important features of CTS in plants and provides methodologies for the investigation of CTS and RNA-binding proteins in plants.

Alternative splicing in tea plants was extensively triggered by drought, heat and their combined stresses

Drought and heat stresses can influence the expressions of genes, and thereby affect the growth and development of plants. Alternative splicing (AS) of genes plays crucial roles through increasing transcriptome diversity in plant stress responses. Tea plants, widely cultivated in the tropics and subtropics, are often simultaneously exposed to drought and heat stresses. In the present study, we performed a global transcriptome of tea leaves treated with drought, heat or their combination. In total, 19,019, 20,025 and 20,253 genes underwent AS in response to drought (DT), heat (HT) and their combined stress (HD), respectively, of which 12,178, 11,912 and 14,413 genes differentially spliced in response to DT, HT and HD, respectively. Also, 2,447 specific differentially spliced genes (DSGs) were found only in response to HD. All DSGs accounted for  48% of the annotated genes in tea tree genome. Comparison of DSGs and differentially expressive genes (DEGs) showed that the proportions of HT and HD-induced DSGs were 13.4% and 9.2%, while the proportion of DT increased to 28.1%. Moreover, the DEG-DSG overlapped genes tended to be enriched in a wide large of pathways in response to DT. The results indicated that the AS of genes in tea leaves was extensively triggered by drought, heat and their combined stresses. In addition, the AS enhanced the transcriptome adaption in response to drought and heat stresses, and the AS also provoked specific molecular functions in response to drought and heat synergy stress. The study might have practical significance for molecular genetic breeding of tea plants with stress resistance.

Novel Amino Acid Assembly in the Silk Tubes of Arid-Adapted Segestriid Spiders

We investigated in different sites inside or outside the Namib Desert the amino acids composition of the protein material forming the tube silk of Ariadna spiders. These spiders belong to the primitive Segestriidae family and spend their life inside vertical silk burrows dug within the sandy and gravelly soil of arid areas. The silks, previously purified by solubilization in hexafluoroisopropanol, were subjected to partial or total acid hydrolysis. Partial hydrolyzed samples, analyzed by mass spectrometry (matrix assisted laser desorption/ionization and electrospray), led to relevant information on the amino acid sequences in the proteins. The free amino acids formed by complete hydrolysis were derivatized with the Marfey's reagent and characterized by electrospray mass spectrometry. The reconstruction of the amino acids highlights a homogeneous plan in the chemical structure of all the analyzed silks. Eight amino acids constituting the primary structure of the proteins were identified. Alanine and glycine are the most abundant ones, with a prevalence of alanine, constituting together at least 61% of the chemical composition of the protein material, differently from what occurs in known spidroins. High percentages of proline, serine and threonine and low percentages of leucine complete the peculiarity of these proteins. The purified silks were also characterized by Fourier-transform Infrared Spectroscopy and their thermal properties were investigated by differential scanning calorimetry. The comparison of the silk tubes among the various Namibian populations, carried out through a multivariate statistical analysis, shows significant differences in their amino acid assembly possibly due to habitat features.

Alternative splicing is a Sorghum bicolor defense response to fungal infection

This study provides new insights that alternative splicing participates with transcriptional control in defense responses to Colletotrichum sublineola in sorghum In eukaryotic organisms, alternative splicing (AS) is an important post-transcriptional mechanism to generate multiple transcript isoforms from a single gene. Protein variants translated from splicing isoforms may have altered molecular characteristics in signal transduction and metabolic activities. However, which transcript isoforms will be translated into proteins and the biological functions of the resulting proteoforms are yet to be identified. Sorghum is one of the five major cereal crops, but its production is severely affected by fungal diseases. For example, sorghum anthracnose caused by Colletotrichum sublineola greatly reduces grain yield and biomass production. In this study, next-generation sequencing technology was used to analyze C. sublineola-inoculated sorghum seedlings compared with mock-inoculated control. It was identified that AS regulation may be as important as traditional transcriptional control during defense responses to fungal infection. Moreover, several genes involved in flavonoid and phenylpropanoid biosynthetic pathways were found to undergo multiple AS modifications. Further analysis demonstrated that non-conventional targets of both 5'- and 3'-splice sites were alternatively used in response to C. sublineola infection. Splicing factors were also affected at both transcriptional and post-transcriptional levels. As the first transcriptome report on C. sublineola infected sorghum, our work also suggested that AS plays crucial functions in defense responses to fungal invasion.

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

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

In-silico prediction of novel genes responsive to drought and salinity stress tolerance in bread wheat (Triticum aestivum)

Common wheat (Triticum aestivum) is the most widely grown cereal crop and is cultivated extensively in dry regions. Water shortage, resulting from either drought or salinity, leads to slow growth and loss of wheat yield. In order to predict new genes responsive to the drought and salt stresses in wheat, 6,717 expressed sequence tags (ESTs), expressed in drought and salinity stress conditions were collected from the National Center for Biotechnology Information (NCBI). The downloaded ESTs were clustered and assembled into 354 contigs; 14 transcription factor families in 29 contigs were identified. In addition, 119 contigs were organized in five enzyme classes. Biological functions were obtained for only 324 of the 354 contigs using gene ontology. In addition, using Kyoto Encyclopedia of Genes and Genomes database, 191 metabolic pathways were identified. The remaining contigs were used for further analysis and the search for new genes responsive to drought and salt stresses. These contigs were mapped on the International Wheat Genome Sequencing Consortium RefSeq v1.0 assembly, the most complete version of the reference sequence of the bread wheat variety Chinese Spring. They were found to have from one to three locations on the subgenomes A, B, and D. Full-length gene sequences were designed for these contigs, which were further validated using promoter analysis. These predicted genes may have applications in molecular breeding programs and wheat drought and salinity research.

Highly ABA-Induced 1 (HAI1)-Interacting protein HIN1 and drought acclimation-enhanced splicing efficiency at intron retention sites

The Highly ABA-Induced 1 (HAI1) protein phosphatase is a central component of drought-related signaling. A screen for HAI1-interacting proteins identified HAI1-Interactor 1 (HIN1), a nuclear protein of unknown function which could be dephosphorylated by HAI1 in vitro. HIN1 colocalization and interaction with serine-arginine rich (SR) splicing factors and appearance of nuclear speckle-localized HIN1 during low water potential (ψ_w) stress suggested a pre-mRNA splicing-related function. RNA sequencing of Arabidopsis Col-0 wild type identified more than 500 introns where moderate severity low ψ_w altered intron retention (IR) frequency. Surprisingly, nearly 90% of these had increased splicing efficiency (decreased IR) during stress. For one-third of these introns, ectopic HIN1 expression (35S:HIN1) in unstressed plants mimicked the increased splicing efficiency seen in stress-treated wild type. HIN1 bound to a GAA-repeat, Exonic Splicing Enhancer-like RNA motif enriched in flanking sequence around HIN1-regulated introns. Genes with stress and HIN1-affected splicing efficiency were enriched for abiotic stress and signaling-related functions. The 35S:HIN1 plants had enhanced growth maintenance during low ψ_w, while hin1 mutants had reduced growth, further indicating the role of HIN1 in drought response. HIN1 is annotated as an MYB/SANT domain protein but has limited homology to other MYB/SANT proteins and is not related to known yeast or metazoan RNA-binding proteins or splicing regulators. Together these data identify HIN1 as a plant-specific RNA-binding protein, show a specific effect of drought acclimation to promote splicing efficiency of IR-prone introns, and also discover HAI1-HIN1 interaction and dephosphorylation that connects stress signaling to splicing regulation.

The Arabidopsis splicing regulator SR45 confers salt tolerance in a splice isoform-dependent manner

Functions of most splice isoforms that are generated by alternative splicing are unknown. We show that two splice variants that encode proteins differing in only eight amino acids have distinct functions in a stress response. Serine/arginine-rich (SR) and SR-like proteins, a conserved family of RNA binding proteins across eukaryotes, play important roles in pre-mRNA splicing and other post-transcriptional processes. Pre-mRNAs of SR and SR-like proteins undergo extensive alternative splicing in response to diverse stresses and produce multiple splice isoforms. However, the functions of most splice isoforms remain elusive. Alternative splicing of pre-mRNA of Arabidopsis SR45, which encodes an SR-like splicing regulator, generates two isoforms (long-SR45.1 and short-SR45.2). The proteins encoded by these two isoforms differ in eight amino acids. Here, we investigated the role of SR45 and its splice variants in salt stress tolerance. The loss of SR45 resulted in enhanced sensitivity to salt stress and changes in expression and splicing of genes involved in regulating salt stress response. Interestingly, only the long isoform (SR45.1) rescued the salt-sensitive phenotype as well as the altered gene expression and splicing patterns in the mutant. These results suggest that SR45 positively regulates salt tolerance. Furthermore, only the long isoform is required for SR45-mediated salt tolerance.

Genome-Wide Analysis of Serine/Arginine-Rich Protein Family in Wheat and Brachypodium distachyon

By regulating the pre-mRNA splicing of other genes and themselves, plant serine/arginine-rich (SR) proteins play important roles in development and in response to abiotic stresses. Presently, the functions of most plant SR protein genes remain unclear. Wheat (Triticumaestivum) and Brachypodiumdistachyon are closely related species. In this study, 40 TaSR and 18 BdSR proteins were identified respectively, and they were classified into seven subfamilies: SR, RS, SCL, RSZ, RS2Z, SC35, and SR45. Similar to Arabidopsis and rice SR protein genes, most TaSR and BdSR protein genes are expressed extensively. Surprisingly, real-time polymerase chain reaction (RT-PCR) analyses showed that no alternative splicing event was found in TaSR protein genes, and only six BdSR protein genes are alternatively spliced genes. The detected alternatively spliced BdSR protein genes and transcripts are much fewer than in Arabidopsis, rice, maize, and sorghum. In the promoter regions, 92 development-related, stress-related, and hormone-related cis-elements were detected, indicating their functions in development and in response to environmental stresses. Meanwhile, 19 TaSR and 16 BdSR proteins were predicted to interact with other SR proteins or non-SR proteins, implying that they are involved in other functions in addition to modulating pre-mRNA splicing as essential components of the spliceosome. These results lay a foundation for further analyses of these genes.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

Chemical Modulation of Alternative Splicing for Molecular-Target Identification by Potential Genetic Control in Agrochemical Research

Alternative splicing (AS), the process of removing introns from pre-mRNA and the rearrangement of exons to produce several types of mature transcripts, is a remarkable step preceding protein synthesis. In particular, it has now been conclusively shown that up to ∼95% of genes are alternatively spliced to generate a complex and diverse proteome in eukaryotic organisms. Consequently, AS is one of the determinants of the functional repertoire of cells. Many studies have revealed that AS in plants can be regulated by cell type, developmental stage, environmental stress, and the circadian clock. Moreover, increasing amounts of evidence reveal that chemical compounds can affect various steps during splicing to induce major effects on plant physiology. Hence, the chemical modulation of AS can serve as a good strategy for molecular-target identification in attempts to potentially control plant genetics. However, the kind of mechanisms involved in the chemical modulation of AS that can be used in agrochemical research remain largely unknown. This review introduces recent studies describing the specific roles AS plays in plant adaptation to environmental stressors and in the regulation of development. We also discuss recent advances in small molecules that induce alterations of AS and the possibility of using this strategy in agrochemical-target identification, giving a new direction for potential genetic control in agrochemical research.

Comparative analysis of transcriptome in two wheat genotypes with contrasting levels of drought tolerance

Drought tolerance is a complex trait that is governed by multiple genes. The study presents differential transcriptome analysis between drought-tolerant (Triticum aestivum Cv. C306) and drought-sensitive (Triticum aestivum Cv. WL711) genotypes, using Affymetrix GeneChip® Wheat Genome Array. Both genotypes exhibited diverse global transcriptional responses under control and drought conditions. Pathway analysis suggested significant induction or repression of genes involved in secondary metabolism, nucleic acid synthesis, protein synthesis, and transport in C306, as compared to WL711. Significant up- and downregulation of transcripts for enzymes, hormone metabolism, and stress response pathways were observed in C306 under drought. The elevated expression of plasma membrane intrinsic protein 1 and downregulation of late embryogenesis abundant in the leaf tissues could play an important role in delayed wilting in C306. The other regulatory genes such as MT, FT, AP2, SKP1, ABA2, ARF6, WRKY6, AOS, and LOX2 are involved in defense response in C306 genotype. Additionally, transcripts with unknown functions were identified as differentially expressed, which could participate in drought responses.

Identification of splice variant of OsGBF1 in Oryza sativa ssp. indica genotypes under salinity stress

G-box-binding factors are plant transcription factors, involved in a wide range of biological processes including abiotic stress responses. In this study, we analyzed the expression of OsGBF1 during salt stress in two contrasting Oryza sativa spp. indica genotypes, Rasi and Tellahamsa. Two-day-old seedlings were exposed to NaCl stress under two different conditions. One set was exposed to 100 mM NaCl before transferring to 250 mM (induction stress), while another set was transferred directly to 250 mM (shock stress). During early induction stress, OsGBF1 was up-regulated in Rasi when compared to Tellahamsa. We cloned full-length OsGBF1 from these two genotypes, and analyzed the sequences. Our analysis indicated the presence of transcript variants, which are designated as OsGBF1a and 1b. OsGBF1b variant retained introns, which lead to the generation of premature termination codon. OsGBF1b transcript levels were not significantly different at 12-h of induction stress in Tellahamsa and Rasi. At 24-h of shock stress, OsGBF1b was up-regulated in both genotypes and the transcript was more in Rasi. Since, OsGBF1a and 1b are predicted to be splice variants, we examined expression pattern of OsSKIP, a splicing factor and component of the spliceosome. In induction stress, OsSKIP was up-regulated at 12- and 24-h in Rasi when compared to Tellahamsa. On the contrary, at 24-h shock stress, OsSKIP was down-regulated in Rasi when compared to Tellahamsa. It is possible that OsSKIP expression was increased in Rasi during induction stress for accurate splicing that could aid in tolerance. This is the first report on OsGBF1 splice variant and the variant could have specific functions linked to stress tolerance in rice.