Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings

Integration of structural study and machine learning to elucidate the RNA-SFs interaction atlas in eukaryotic cells

Alternative splicing (AS) occupies a central position in plant growth and development, stress response, and animal growth and disease processes. Mutations in SF (splicing factor) trigger aberrant AS activities that disrupt these fine biological processes. Although cryo electron microscopy (cryoEM) technology has successfully revealed the fine structure of multiple spliceosomes, the dynamic and complex network of RNA-SFs remains to be fully resolved. This review summarizes the binding patterns of RNA and SFs through machine learning's powerful computational capabilities, the deep structural analysis using cryoEM, and experimental validation of RNA protein binding. Connect RNA protein interaction experiments, high-resolution imaging capabilities of cryoEM, and powerful analytical capabilities of machine learning to jointly construct a detailed RNA-SFs interaction map, forming a powerful toolkit. These knowledge help us better understand the complexity and working mechanisms of biological systems. This article not only has profound significance in revealing the molecular mechanisms of diseases and developing multi-target efficient drugs but also provides in-depth insights into molecular breeding and plant resistance enhancement.

Alternative Splicing of BnABF4L Mediates Response to Abiotic Stresses in Rapeseed (Brassica napus L.)

ABRE BINDING FACTOR 4 (ABF4) is a pivotal regulatory gene in the abscisic acid (ABA) signaling pathway, and changes in its expression levels can modulate the plant's stress resistance. To further explore the specific regulatory mechanisms of alternative splicing (AS) in the ABA signaling pathway and to identify new breakthroughs for breeding high stress-resistant varieties of Brassica napus, we identified 17 homologous genes of ABF4 in the genome. Utilizing bioinformatics techniques, we analyzed their motifs, conserved domains, and cis-acting elements of their promoters. Through transcriptome data from the stress-tolerant dwarf strain ndf2 and its parental line 3529, we uncovered a significantly differentially expressed ABF4 gene, which we named BnABF4L. Subsequently, we analyzed the AS events of BnABF4L under normal growth conditions and different abiotic stresses, as well as the impact of different transcript variants' 5' untranslated region (5'UTR) on gene translation. BnABF4L undergoes alternative 3' splice site (A3SS) selection to produce three transcripts (V1-V3) with divergent 5'UTRs. While V1 translation is suppressed by upstream ORFs (uORFs), V2/V3 exhibit enhanced translational efficiency. Under stress, ndf2 shifts splicing toward V3, circumventing uORF-mediated repression to upregulate stress-adapted isoforms. We validated the inhibitory effect of upstream open reading frames (uORFs) on protein-coding open reading frame (pORFs) and, based on the collective experimental results, proposed the flexible regulatory mechanism of AS events of BnABF4L in response to stress. Our findings provide new insights for future studies on stress resistance in rapeseed as well as for research on the regulation of alternative splicing mechanisms in the ABA signaling pathway.

Rice extracellular vesicles send defense proteins into fungus Rhizoctonia solani to reduce disease

The exchange of molecular information across kingdoms is crucial for the survival of both plants and their pathogens. Recent research has identified that plants transfer their small RNAs and microRNAs into fungal pathogens to suppress infection. However, whether and how plants send defense proteins into pathogens remains unknown. Here, we report that rice (Oryza sativa) plants package defense proteins into extracellular vesicles (EVs) and deliver them to the fungal pathogen Rhizoctonia solani. These EVs, enriched with host defense proteins, are internalized by the fungal cells. Reducing the transfer of host defense proteins via EVs results in increased disease susceptibility. Furthermore, the overexpression of host defense proteins in either rice plants or the fungal cells reduced the infection. Therefore, plants use EVs to send defense proteins into fungal pathogens, thereby combating infection. This mechanism represents a form of protein exchange between plants and pathogens, which contributes to reducing crop diseases.

Splicing and expression dynamics of SR genes in hot pepper (Capsicum annuum): regulatory diversity and conservation under stress

In this study, we identified and characterized 23 genes encoding serine/arginine-rich (SR) protein in hot pepper (Capsicum annuum), named CaSR here. These CaSR proteins are grouped into seven subfamilies. Phylogenetic analysis revealed a high degree of similarity between CaSRs and their homologous proteins in other plants. Promoter regions of SR proteins are enriched with elements relating to light response, stress, hormone signaling, and plant growth. Notably, transcription levels of several proteins, including CaSR33, CaSR34, and CaSR34a, were upregulated by salt, drought, and cold stresses, suggesting potential roles of these proteins in stress tolerance. We also observed an increase of CaSR transcript population resulting from alternative splicing (AS) regulation, mainly intron retention. AS patterns of CaSR genes varied among tissues. Higher AS intensity was found in the RS subfamily, while some genes in the RSZ subfamily showed no AS regulation under the conditions used here. Interestingly, a cross-species comparative study with Arabidopsis (Arabidopsis thaliana) and tomato (Solanum lycopersicum) showed that many AS events impact the region which codes the RNA recognition motif (RRM) domain of the protein, indicating a conserved regulatory mechanism of SR proteins in plants. Our findings reveal the functional diversity and evolutionary conservation of SR proteins in hot pepper and highlight AS as a mechanism enhancing plant adaptability, providing insights for future stress-resistant crop development.

Alternative splicing of VRF1 acts as a molecular switch to regulate stress-induced early flowering

Plants frequently evade extreme environmental stress by initiating early flowering, yet the underlying mechanisms remain incompletely understood. Here, through extensive mutant screening, we identify a vegetative growth to reproductive growth transition factor (vrf1) mutant, which exhibits a deficiency in drought escape. Alternative splicing of VRF1 generates four isoforms, of which two encode functional proteins, VRF1-AS1 and VRF1-AS3. The proportions of VRF1-AS1 and VRF1-AS3 are modulated by environmental factors, serving as a molecular switch mediating the transition from tolerance to early flowering, irrespective of their absolute abundance. VRF1-AS1 and VRF1-AS3 competitively bind to MKK1, resulting in MKK1 phosphorylation at different sites, which opens two distinct regulatory pathways. The role of VRF1 is conserved across various conditions, and 66 Arabidopsis ecotypes suggest its central function in stress avoidance through premature flowering. In summary, our findings show that plants respond precisely and effectively to dynamic environmental changes by modulating their isoform ratios.

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

The Zea mays PeptideAtlas: A New Maize Community Resource

This study presents the Maize PeptideAtlas resource (www.peptideatlas.org/builds/maize) to help solve questions about the maize proteome. Publicly available raw tandem mass spectrometry (MS/MS) data for maize collected from ProteomeXchange were reanalyzed through a uniform processing and metadata annotation pipeline. These data are from a wide range of genetic backgrounds and many sample types and experimental conditions. The protein search space included different maize genome annotations for the B73 inbred line from MaizeGDB, UniProtKB, NCBI RefSeq, and for the W22 inbred line. 445 million MS/MS spectra were searched, of which 120 million were matched to 0.37 million distinct peptides. Peptides were matched to 66.2% of proteins in the most recent B73 nuclear genome annotation. Furthermore, most conserved plastid- and mitochondrial-encoded proteins (NCBI RefSeq annotations) were identified. Peptides and proteins identified in the other B73 genome annotations will improve maize genome annotation. We also illustrate the high-confidence detection of unique W22 proteins. N-terminal acetylation, phosphorylation, ubiquitination, and three lysine acylations (K-acetyl, K-malonyl, and K-hydroxyisobutyryl) were identified and can be inspected through a PTM viewer in PeptideAtlas. All matched MS/MS-derived peptide data are linked to spectral, technical, and biological metadata. This new PeptideAtlas is integrated in MaizeGDB with a peptide track in JBrowse.

Integrative analyses of long and short-read RNA sequencing reveal the spliced isoform regulatory network of seedling growth dynamics in upland cotton

The polyploid genome of cotton has significantly increased the transcript complexity. Recent advances in full-length transcript sequencing are now widely used to characterize the complete landscape of transcriptional events. Such studies in cotton can help us to explore the genetic mechanisms of the cotton seedling growth. Through long-read single-molecule RNA sequencing, this study compared the transcriptomes of three yield contrasting genotypes of upland cotton. Our analysis identified different numbers of spliced isoforms from 31,166, 28,716, and 28,713 genes in SJ48, Z98, and DT8 cotton genotypes, respectively, most of which were novel compared to previous cotton reference transcriptomes, and showed significant differences in the number of exon structures and coding sequence length due to intron retention. Quantification of isoform expression revealed significant differences in expression in the root and leaf of each genotype. An array of key isoform target genes showed protein kinase or phosphorylation functions, and their protein interaction network contained most of the circadian oscillator proteins. Spliced isoforms from the GIGANTEA (GI) protien were differentially regulated in each genotype and might be expected to regulate translational activities, including the sequence and function of target proteins. In addition, these spliced isoforms generate diurnal expression profiles in cotton leaves, which may alter the transcriptional regulatory network of seedling growth. Silencing of the novel spliced GI isoform Gh_A02G0645_N17 significantly affected biomass traits, contributed to variable growth, and increased transcription of the early flowering pathway gene ELF in cotton. Our high-throughput hybrid sequencing results will be useful to dissect functional differences among spliced isoforms in the polyploid cotton genome.

A cotton endoreduplication gene, GaTOP6B, regulates trichome branching development

Trichomes are specialized epidermal structures that protect plants from biotic and abiotic stresses by synthesizing, storing, and secreting defensive compounds. This study investigates the role of the Gossypium arboreum DNA topoisomerase VI subunit B gene (GaTOP6B) in trichome development and branching. Sequence alignment revealed a high similarity between GaTOP6B and AtTOP6B, suggesting a conserved function in trichome regulation. Although AtTOP6B acts as a positive regulator of trichome development, functional analyses showed contrasting effects: Virus-induced gene silencing (VIGS) of GaTOP6B in cotton increased trichome density, while its overexpression in Arabidopsis decreased trichome density but enhanced branching. This demonstrates that GaTOP6B negatively regulates trichome number, indicating species-specific roles in trichome initiation and branching between cotton and Arabidopsis. Overexpression of the GaTOP6B promotes jasmonic acid synthesis, which in turn inhibits the G1/S or G2/M transitions, stalling the cell cycle. On the other hand, it suppresses brassinolide synthesis and signaling while promoting cytokinin degradation, further inhibiting mitosis. These hormonal interactions facilitate the transition of cells from the mitotic cycle to the endoreduplication cycle. As the level of endoreduplication increases, trichomes develop an increased number of branches. These findings highlight GaTOP6B's critical role as a regulator of trichome development, providing new genetic targets for improving cotton varieties in terms of enhanced adaptability and resilience.

Evolution and stress response potential of the plant splicing factor U1C

Alternative splicing is a crucial process in multicellular eukaryote, facilitated by the assembly of spliceosomal complexes comprising numerous small ribonucleoproteins. At an early stage, U1C is thought to be required for 5' splice site recognition and base pairing. However, a systematic analysis of the U1C gene family in response to developmental cues and stress conditions has not yet been conducted in plants. This study identified 114 U1C genes in 72 plant species using basic bioinformatics analyses. Phylogenetic analysis was used to compare gene and protein structures, promoter motifs, and tissue- and stress-specific expression levels, revealing their functional commonalities or diversity in response to developmental cues, such as embryonic expression, or stress treatments, including drought and heat. Fluorescence quantitative expression analysis showed that U1C gene expression changed under salt, low temperature, drought, and Cd stress in rice seedlings. However, gene expression in shoots and roots was not consistent under different stress conditions, suggesting a complex regulatory mechanism. This research provides foundational insights into the U1C gene family's role in plant development and stress responses, highlighting potential targets for future studies.

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.

Alternative splicing as a driver of natural variation in abscisic acid response

Abscisic acid (ABA) is a crucial player in plant responses to the environment. It accumulates under stress, activating downstream signaling to implement molecular responses that restore homeostasis. Natural variance in ABA sensitivity remains barely understood, and the ABA pathway has been mainly studied at the transcriptional level, despite evidence that posttranscriptional regulation, namely, via alternative splicing, contributes to plant stress tolerance. Here, we identified the Arabidopsis accession Kn-0 as less sensitive to ABA than the reference Col-0, as shown by reduced effects of the hormone on seedling establishment, root branching, and stomatal closure, as well as by decreased induction of ABA marker genes. An in-depth comparative transcriptome analysis of the ABA response in the two variants revealed lower expression changes and fewer genes affected for the least ABA-sensitive ecotype. Notably, Kn-0 exhibited reduced levels of the ABA-signaling SnRK2 protein kinases and lower basal expression of ABA-reactivation genes, consistent with our finding that Kn-0 contains less endogenous ABA than Col-0. ABA also markedly affected alternative splicing, primarily intron retention, with Kn-0 being less responsive regarding both the number and magnitude of alternative splicing events, particularly exon skipping. We find that alternative splicing introduces a more ecotype-specific layer of ABA regulation and identify ABA-responsive splicing changes in key ABA pathway regulators that provide a functional and mechanistic link to the differential sensitivity of the two ecotypes. Our results offer new insight into the natural variation of ABA responses and corroborate a key role for alternative splicing in implementing ABA-mediated stress responses.

Identification of key genes for triacylglycerol biosynthesis and storage in herbaceous peony (Paeonia lactifolra Pall.) seeds based on full-length transcriptome

BACKGROUND

The herbaceous peony (Paeonia lactiflora Pall.) is extensively cultivated in China due to its root being used as a traditional Chinese medicine known as 'Radix Paeoniae Alba'. In recent years, it has been discovered that its seeds incorporate abundant unsaturated fatty acids, thereby presenting a potential new oilseed plant. Surprisingly, little is known about the full-length transcriptome sequencing of Paeonia lactiflora, limiting research into its gene function and molecular mechanisms.

RESULTS

A total of 484,931 Reads of Inserts (ROI) sequences and 1,455,771 full-Length non-chimeric reads (FLNC) sequences were obtained for CDS prediction, TF analysis, SSR analysis and lncRNA identification. In addition, gene function annotation and gene structure analysis were performed. A total of 4905 transcripts were related to lipid metabolism biosynthesis pathway, belonging to 28 enzymes. We use these data to identify 10 oleosin (OLE) and 5 diacylglycerol acyltransferase (DGAT) gene members after de-redundancy. The analysis of physicochemical properties and secondary structure showed them similarity in gene family respectively. The phylogenetic analysis showed that the distribution of OLE and DGAT family members was roughly the same as that of Arabidopsis. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed expression changes in different seed development stages, and showed a trend of increasing and then decreasing.

CONCLUSION

In summary, these results provide new insights into the molecular mechanism of triacylglycerol (TAG) biosynthesis and storage during the seedling stage in Paeonia lactiflora. It provides theoretical references for selecting and breeding oil varieties and understanding the functions of oil storage as well as lipid synthesis related genes in Paeonia lactiflora.

E3 SUMO ligase SIZ1 splicing variants localize and function according to external conditions

SIZ1 (SAP and MIZ1) is a member of the Siz/PIAS-type RING family of E3 SUMO (small ubiquitin-related modifier) ligases that play key roles in growth, development, and stress responses in plant and animal systems. Nevertheless, splicing variants of SIZ1 have not yet been characterized. Here, we identified four splicing variants of Arabidopsis (Arabidopsis thaliana) SIZ1, which encode three different protein isoforms. The SIZ1 gene encodes an 873-amino acid (aa) protein. Among the four SIZ1 splicing variants (SSVs), SSV1 and SSV4 encode identical 885 aa proteins; SSV2 encodes an 832 aa protein; and SSV3 encodes an 884 aa protein. SSV2 mainly localized to the plasma membrane, whereas SIZ1, SSV1/SSV4, and SSV3 localized to the nucleus. Interestingly, SIZ1 and all SSVs exhibited similar E3 SUMO ligase activities and preferred SUMO1 and SUMO2 for their E3 ligase activity. Transcript levels of SSV2 were substantially increased by heat treatment, while those of SSV1, SSV3, and SSV4 transcripts were unaffected by various abiotic stresses. SSV2 directly interacted with and sumoylated cyclic nucleotide-gated ion channel 6 (CNGC6), a positive thermotolerance regulator, enhancing the stability of CNGC6. Notably, transgenic siz1-2 mutants expressing SSV2 exhibited greater heat stress tolerance than wild-type plants, whereas those expressing SIZ1 were sensitive to heat stress. Furthermore, transgenic cngc6 plants overaccumulating a mutated mCNGC6 protein (K347R, a mutation at the sumoylation site) were sensitive to heat stress, similar to the cngc6 mutants, while transgenic cngc6 plants overaccumulating CNGC6 exhibited restored heat tolerance. Together, we propose that alternative splicing is an important mechanism that regulates the function of SSVs during development or under adverse conditions, including heat stress.

Transcription and splicing variations of SR genes accompany with genome-wide accumulation of long-introns in pine

Most of mRNAs in Eukaryote were matured after the removal of introns in their pre-mRNA transcripts. Serine/arginine-rich (SR) proteins are a group of splicing regulators regulating the splicing processes globally. Expressions of SR proteins themselves were extensively regulated, at both transcription and splicing levels, under different environmental conditions, specially heat stress conditions. The pine genome is characterized by super-long and easily methylated introns in a large number of genes that derived from the extensive accumulation of transposons (TEs). Here, we identified and analyzed the phylogenetic characteristics of 24 SR proteins and their encoding genes from the pine genome. Then we explored transcription and pre-mRNA splicing expression patterns of SR genes in P. massoniana seedlings under normal and heat stress temperature conditions. Our results showed that the transcription patterns of SR genes in pine exhibited significant changes compared to other plant species, and these changes were not strictly correlated with the intron length and DNA methylation intensity of the SR genes. Interestingly, none of the long introns of SR genes underwent alternative splicing (AS) in our experiment. Furthermore, the intensity of AS regulation may be related to the potential DNA methylation intensity of SR genes. Taken together, this study explores for the first time the characteristics of significant variations in the transcription and splicing patterns of SR proteins in a plant species with an over-accumulation of super-long introns.

Plant genome information facilitates plant functional genomics

MAIN CONCLUSION

In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Proteogenomics-based functional genome research: approaches, applications, and perspectives in plants

Proteogenomics (PG) integrates the proteome with the genome and transcriptome to refine gene models and annotation. Coupled with single-cell (SC) assays, PG effectively distinguishes heterogeneity among cell groups. Affiliating spatial information to PG reveals the high-resolution circuitry within SC atlases. Additionally, PG can investigate dynamic changes in protein-coding genes in plants across growth and development as well as stress and external stimulation, significantly contributing to the functional genome. Here we summarize existing PG research in plants and introduce the technical features of various methods. Combining PG with other omics, such as metabolomics and peptidomics, can offer even deeper insights into gene functions. We argue that the application of PG will represent an important font of foundational knowledge for plants.

Physiology and transcriptomics highlight the underlying mechanism of sunflower responses to drought stress and rehydration

Drought can adversely influence the crop growth and production. Accordingly, sunflowers have strong adaptability to drought; hence, we conducted analyses for sunflower seedlings with drought stress and rehydration drought acclimation through physiological measurements and transcriptomics. It showed that drought can cause the accumulation of ROS and enhance the activity of antioxidant enzymes and the content of osmolytes. After rehydration, the contents of ROS and MDA were significantly reduced concomitant with increased antioxidant activity and osmotic adjustment. Totally, 2,589 DEGs were identified among treatments. Functional enrichment analysis showed that DEGs were mainly involved in plant hormone signal transduction, MAPK signaling, and biosynthesis of secondary metabolites. Comparison between differentially spliced genes and DEGs indicated that bHLH025, NAC53, and SINAT3 may be pivotal genes involved in sunflower drought resistance. Our results not only highlight the underlying mechanism of drought stress and rehydration in sunflower but also provide a theoretical basis for crop genetic breeding.

Advances in alternative splicing identification: deep learning and pantranscriptome

In plants, alternative splicing is a crucial mechanism for regulating gene expression at the post-transcriptional level, which leads to diverse proteins by generating multiple mature mRNA isoforms and diversify the gene regulation. Due to the complexity and variability of this process, accurate identification of splicing events is a vital step in studying alternative splicing. This article presents the application of alternative splicing algorithms with or without reference genomes in plants, as well as the integration of advanced deep learning techniques for improved detection accuracy. In addition, we also discuss alternative splicing studies in the pan-genomic background and the usefulness of integrated strategies for fully profiling alternative splicing.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Eukaryotic splicing machinery in the plant-virus battleground

Plant virual infections are mainly caused by plant-virus parasitism which affects ecological communities. Some viruses are highly pathogen specific that can infect only specific plants, while some can cause widespread harm, such as tobacco mosaic virus (TMV) and cucumber mosaic virus (CMV). After a virus infects the host, undergoes a series of harmful effects, including the destruction of host cell membrane receptors, changes in cell membrane components, cell fusion, and the production of neoantigens on the cell surface. Therefore, competition between the host and the virus arises. The virus starts gaining control of critical cellular functions of the host cells and ultimately affects the fate of the targeted host plants. Among these critical cellular processes, alternative splicing (AS) is an essential posttranscriptional regulation process in RNA maturation, which amplify host protein diversity and manipulates transcript abundance in response to plant pathogens. AS is widespread in nearly all human genes and critical in regulating animal-virus interactions. In particular, an animal virus can hijack the host splicing machinery to re-organize its compartments for propagation. Changes in AS are known to cause human disease, and various AS events have been reported to regulate tissue specificity, development, tumour proliferation, and multi-functionality. However, the mechanisms underlying plant-virus interactions are poorly understood. Here, we summarize the current understanding of how viruses interact with their plant hosts compared with humans, analyze currently used and putative candidate agrochemicals to treat plant-viral infections, and finally discussed the potential research hotspots in the future. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Multi-Omics Pipeline and Omics-Integration Approach to Decipher Plant's Abiotic Stress Tolerance Responses

The present day's ongoing global warming and climate change adversely affect plants through imposing environmental (abiotic) stresses and disease pressure. The major abiotic factors such as drought, heat, cold, salinity, etc., hamper a plant's innate growth and development, resulting in reduced yield and quality, with the possibility of undesired traits. In the 21st century, the advent of high-throughput sequencing tools, state-of-the-art biotechnological techniques and bioinformatic analyzing pipelines led to the easy characterization of plant traits for abiotic stress response and tolerance mechanisms by applying the 'omics' toolbox. Panomics pipeline including genomics, transcriptomics, proteomics, metabolomics, epigenomics, proteogenomics, interactomics, ionomics, phenomics, etc., have become very handy nowadays. This is important to produce climate-smart future crops with a proper understanding of the molecular mechanisms of abiotic stress responses by the plant's genes, transcripts, proteins, epigenome, cellular metabolic circuits and resultant phenotype. Instead of mono-omics, two or more (hence 'multi-omics') integrated-omics approaches can decipher the plant's abiotic stress tolerance response very well. Multi-omics-characterized plants can be used as potent genetic resources to incorporate into the future breeding program. For the practical utility of crop improvement, multi-omics approaches for particular abiotic stress tolerance can be combined with genome-assisted breeding (GAB) by being pyramided with improved crop yield, food quality and associated agronomic traits and can open a new era of omics-assisted breeding. Thus, multi-omics pipelines together are able to decipher molecular processes, biomarkers, targets for genetic engineering, regulatory networks and precision agriculture solutions for a crop's variable abiotic stress tolerance to ensure food security under changing environmental circumstances.

Alternative Splicing Analysis Revealed the Role of Alpha-Linolenic Acid and Carotenoids in Fruit Development of Osmanthus fragrans

Alternative splicing refers to the process of producing different splicing isoforms from the same pre-mRNA through different alternative splicing events, which almost participates in all stages of plant growth and development. In order to understand its role in the fruit development of Osmanthus fragrans, transcriptome sequencing and alternative splicing analysis was carried out on three stages of O. fragrans fruit (O. fragrans "Zi Yingui"). The results showed that the proportion of skipping exon events was the highest in all three periods, followed by a retained intron, and the proportion of mutually exclusive exon events was the lowest and most of the alternative splicing events occurred in the first two periods. The results of enrichment analysis of differentially expressed genes and differentially expressed isoforms showed that alpha-Linolenic acid metabolism, flavonoid biosynthesis, carotenoid biosynthesis, photosynthesis, and photosynthetic-antenna protein pathways were significantly enriched, which may play an important role in the fruit development of O. fragrans. The results of this study lay the foundation for further study of the development and maturation of O. fragrans fruit and further ideas for controlling fruit color and improving fruit quality and appearance.

Plant serine/arginine-rich proteins: versatile players in RNA processing

MAIN CONCLUSION

Serine/arginine-rich (SR) proteins participate in RNA processing by interacting with precursor mRNAs or other splicing factors to maintain plant growth and stress responses. Alternative splicing is an important mechanism involved in mRNA processing and regulation of gene expression at the posttranscriptional level, which is the main reason for the diversity of genes and proteins. The process of alternative splicing requires the participation of many specific splicing factors. The SR protein family is a splicing factor in eukaryotes. The vast majority of SR proteins' existence is an essential survival factor. Through its RS domain and other unique domains, SR proteins can interact with specific sequences of precursor mRNA or other splicing factors and cooperate to complete the correct selection of splicing sites or promote the formation of spliceosomes. They play essential roles in the composition and alternative splicing of precursor mRNAs, providing pivotal functions to maintain growth and stress responses in animals and plants. Although SR proteins have been identified in plants for three decades, their evolutionary trajectory, molecular function, and regulatory network remain largely unknown compared to their animal counterparts. This article reviews the current understanding of this gene family in eukaryotes and proposes potential key research priorities for future functional studies.

PlaASDB: a comprehensive database of plant alternative splicing events in response to stress

BACKGROUND

Alternative splicing (AS) is a co-transcriptional regulatory mechanism of plants in response to environmental stress. However, the role of AS in biotic and abiotic stress responses remains largely unknown. To speed up our understanding of plant AS patterns under different stress responses, development of informative and comprehensive plant AS databases is highly demanded.

DESCRIPTION

In this study, we first collected 3,255 RNA-seq data under biotic and abiotic stresses from two important model plants (Arabidopsis and rice). Then, we conducted AS event detection and gene expression analysis, and established a user-friendly plant AS database termed PlaASDB. By using representative samples from this highly integrated database resource, we compared AS patterns between Arabidopsis and rice under abiotic and biotic stresses, and further investigated the corresponding difference between AS and gene expression. Specifically, we found that differentially spliced genes (DSGs) and differentially expressed genes (DEG) share very limited overlapping under all kinds of stresses, suggesting that gene expression regulation and AS seemed to play independent roles in response to stresses. Compared with gene expression, Arabidopsis and rice were more inclined to have conserved AS patterns under stress conditions.

CONCLUSION

PlaASDB is a comprehensive plant-specific AS database that mainly integrates the AS and gene expression data of Arabidopsis and rice in stress response. Through large-scale comparative analyses, the global landscape of AS events in Arabidopsis and rice was observed. We believe that PlaASDB could help researchers understand the regulatory mechanisms of AS in plants under stresses more conveniently. PlaASDB is freely accessible at http://zzdlab.com/PlaASDB/ASDB/index.html .

Alternative Splicing: From Abiotic Stress Tolerance to Evolutionary Genomics

The post-transcriptional regulation of gene expression, in particular alternative splicing (AS) events, substantially contributes to the complexity of eukaryotic transcriptomes and proteomes [...]

Global Analysis of Dark- and Heat-Regulated Alternative Splicing in Arabidopsis

Alternative splicing (AS) is one of the major post-transcriptional regulation mechanisms that contributes to plant responses to various environmental perturbations. Darkness and heat are two common abiotic factors affecting plant growth, yet the involvement and regulation of AS in the plant responses to these signals remain insufficiently examined. In this study, we subjected Arabidopsis seedlings to 6 h of darkness or heat stress and analyzed their transcriptome through short-read RNA sequencing. We revealed that both treatments altered the transcription and AS of a subset of genes yet with different mechanisms. Dark-regulated AS events were found enriched in photosynthesis and light signaling pathways, while heat-regulated AS events were enriched in responses to abiotic stresses but not in heat-responsive genes, which responded primarily through transcriptional regulation. The AS of splicing-related genes (SRGs) was susceptible to both treatments; while dark treatment mostly regulated the AS of these genes, heat had a strong effect on both their transcription and AS. PCR analysis showed that the AS of the Serine/Arginine-rich family gene SR30 was reversely regulated by dark and heat, and heat induced the upregulation of multiple minor SR30 isoforms with intron retention. Our results suggest that AS participates in plant responses to these two abiotic signals and reveal the regulation of splicing regulators during these processes.

SWATH-MS-based proteogenomic analysis reveals the involvement of alternative splicing in poplar upon lead stress

Alternative splicing (AS) regulates gene expression and increases proteomic diversity for the fine tuning of stress responses in plants, but the exact mechanism through which AS functions in plant stress responses is not thoroughly understood. Here, we investigated how AS functions in poplar (Populus trichocarpa), a popular plant for bioremediation, in response to lead (Pb) stress. Using a proteogenomic analysis, we determine that Pb stress induced alterations in AS patterns that are characterized by an increased use of nonconventional splice sites and a higher abundance of Pb-responsive splicing factors (SFs) associated with Pb-responsive transcription factors. A strong Pb(II)-inducible chaperone protein, PtHSP70, that undergoes AS was further characterized. Overexpression of its two spliced isoforms, PtHSP70-AS1 and PtHSP70-AS2, in poplar and Arabidopsis significantly enhances the tolerance to Pb. Further characterization shows that both isoforms can directly bind to Pb(II), and PtHSP70-AS2 exhibits 10-fold higher binding capacities and a greater increase in expression under Pb stress, thereby reducing cellular toxicity through Pb(II) extrusion and conferring Pb tolerance. AS of PtHSP70 is found to be regulated by PtU1-70K, a Pb(II)-inducible core SF involved in 5'-splice site recognition. Because the same splicing pattern is also found in HSP70 orthologs in other plant species, AS of HSP70 may be a common regulatory mechanism to cope with Pb(II) toxicity. Overall, we have revealed a novel post-transcriptional machinery that mediates heavy metal tolerance in diverse plant species. Our findings offer new molecular targets and bioengineering strategies for phytoremediation and provide new insight for future directions in AS research.

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 Integration of Data from Different Long-Read Sequencing Platforms Enhances Proteoform Characterization in Arabidopsis

The increasing availability of massive omics data requires improving the quality of reference databases and their annotations. The combination of full-length isoform sequencing (Iso-Seq) with short-read transcriptomics and proteomics has been successfully used for increasing proteoform characterization, which is a main ongoing goal in biology. However, the potential of including Oxford Nanopore Technologies Direct RNA Sequencing (ONT-DRS) data has not been explored. In this paper, we analyzed the impact of combining Iso-Seq- and ONT-DRS-derived data on the identification of proteoforms in Arabidopsis MS proteomics data. To this end, we selected a proteomics dataset corresponding to senescent leaves and we performed protein searches using three different protein databases: AtRTD2 and AtRTD3, built from the homonymous transcriptomes, regarded as the most complete and up-to-date available for the species; and a custom hybrid database combining AtRTD3 with publicly available ONT-DRS transcriptomics data generated from Arabidopsis leaves. Our results show that the inclusion and combination of long-read sequencing data from Iso-Seq and ONT-DRS into a proteogenomic workflow enhances proteoform characterization and discovery in bottom-up proteomics studies. This represents a great opportunity to further investigate biological systems at an unprecedented scale, although it brings challenges to current protein searching algorithms.

Nanopore long-read RNAseq reveals transcriptional variations in citrus species

The number of studies on plant transcriptomes using ONT RNAseq technology is rapidly increasing in recent. It is a powerful method to decipher transcriptomic complexity, particularly alternative splicing (AS) event detection. Citrus plants are the most important widely grown fruit crops. Exploring different AS events in citrus contributes to transcriptome improvement and functional genome study. Here, we performed ONT RNAseq in 9 species (Atalantia buxifolia, Citrus clementina, C. grandis, C. ichangensis, C. reticulata, C. sinensis, Clausena lansium, Fortunella hindsii, and Poncirus trifoliata), accompanied with Illumina sequencing. Non-redundant full-length isoforms were identified between 41,957 and 76,974 per species. Systematic analysis including different types of isoforms, number of isoforms per gene locus, isoform distribution, ORFs and lncRNA prediction and functional annotation were performed mainly focused on novel isoforms, unraveling the capability of novel isoforms detection and characterization. For AS events prediction, A3, RI, and AF were overwhelming types across 9 species. We analyzed isoform similarity and evolutionary relationships in all species. We identified that multiple isoforms derived from orthologous single copy genes among different species were annotated as enzymes, nuclear-related proteins or receptors. Isoforms with extending sequences on 5', 3', or both compared with reference genome were filtered out to provide information for transcriptome improvement. Our results provide novel insight into comprehending complex transcriptomes in citrus and valuable information for further investigation on the function of genes with diverse isoforms.

Comparison of the Phenotypic Performance, Molecular Diversity, and Proteomics in Transgenic Rice

The extent of molecular diversity and differentially expressed proteins (DEPs) in transgenic lines provide valuable information to understand the phenotypic performance of transgenic crops compared with their parents. Here, we compared the differences in the phenotypic variation of twelve agronomic and end-use quality traits, the extent of microsatellite diversity, and DEPs of a recurrent parent line with three transgenic rice restorer lines carrying either CRY1C gene on chromosome 11 or CRY2A gene on chromosome 12 or both genes. The three transgenic lines had significantly smaller stem borer infestation than the recurrent parent without showing significant differences among most agronomic traits, yield components, and end-use quality traits. Using 512 microsatellite markers, the three transgenic lines inherited 2.9-4.3% of the Minghui 63 donor genome and 96.3-97.1% of the CH891 recurrent parent genome. As compared with the recurrent parent, the number of upregulated and down-regulated proteins in the three transgenic lines varied from 169 to 239 and from 131 to 199, respectively. Most DEPs were associated with the secondary metabolites biosynthesis transport and catabolism, carbohydrate transport and metabolism, post-translational modification, and signal transduction mechanisms. Although several differentially expressed proteins were observed between transgenic rice and its recurrent parent, the differences may not have been associated with grain yield and most other phenotypic traits in transgenic rice.

Comparative Analysis of Alternative Splicing in Two Contrasting Apple Cultivars Defense against Alternaria alternata Apple Pathotype Infection

Alternaria blotch disease, caused by the Alternaria alternata apple pathotype (A. alternata AP), is one of the most serious fungal diseases in apples. Alternative splicing (AS), one of the pivotal post-transcriptional regulatory mechanisms, plays essential roles in various disease resistance responses. Here, we performed RNA-Seq for two apple cultivars (resistant cultivar 'Jonathan' (J) and susceptible cultivar 'Starking Delicious' (SD)) infected by A. alternata AP to further investigate their AS divergence. In total, 1454, 1780, 1367 and 1698 specifically regulated differential alternative splicing (DAS) events were detected in J36, J72, SD36 and SD72 groups, respectively. Retained intron (RI) was the dominant AS pattern. Conformably, 642, 764, 585 and 742 uniquely regulated differentially spliced genes (DSGs) were found during A. alternata AP infection. Comparative analysis of AS genes in differential splicing and expression levels suggested that only a small proportion of DSGs overlapped with differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis demonstrated that the DSGs were significantly enriched at multiple levels of gene expression regulation. Briefly, the specific AS was triggered in apple defense against A. alternata AP. Therefore, this study facilitates our understanding on the roles of AS regulation in response to A. alternata AP infection in apples.

Integrated Genomic and Transcriptomic Elucidation of Flowering in Garlic

Commercial cultivars of garlic are sterile, and therefore efficient breeding of this crop is impossible. Recent restoration of garlic fertility has opened new options for seed production and hybridization. Transcriptome catalogs were employed as a basis for garlic genetic studies, and in 2020 the huge genome of garlic was fully sequenced. We provide conjoint genomic and transcriptome analysis of the regulatory network in flowering garlic genotypes. The genome analysis revealed phosphatidylethanolamine-binding proteins (PEBP) and LEAFY (LFY) genes that were not found at the transcriptome level. Functions of TFL-like genes were reduced and replaced by FT-like homologs, whereas homologs of MFT-like genes were not found. The discovery of three sequences of LFY-like genes in the garlic genome and confirmation of their alternative splicing suggest their role in garlic florogenesis. It is not yet clear whether AsLFY1 acts alone as the "pioneer transcription factor" or AsLFY2 also provides these functions. The presence of several orthologs of flowering genes that differ in their expression and co-expression network advocates ongoing evolution in the garlic genome and diversification of gene functions. We propose that the process of fertility deprivation in garlic cultivars is based on the loss of transcriptional functions of the specific genes.

Genome-scale analysis of Arabidopsis splicing-related protein kinase families reveals roles in abiotic stress adaptation

Nearly 60 - 80 % of intron-containing plant genes undergo alternative splicing in response to either stress or plant developmental cues. RNA splicing is performed by a large ribonucleoprotein complex called the spliceosome in conjunction with associated subunits such as serine arginine (SR) proteins, all of which undergo extensive phosphorylation. In plants, there are three main protein kinase families suggested to phosphorylate core spliceosome subunits and related splicing factors based on orthology to human splicing-related kinases: the SERINE/ARGININE PROTEIN KINASES (SRPK), ARABIDOPSIS FUS3 COMPLEMENT (AFC), and Pre-mRNA PROCESSING FACTOR 4 (PRP4K) protein kinases. To better define the conservation and role(s) of these kinases in plants, we performed a genome-scale analysis of the three families across photosynthetic eukaryotes, followed by extensive transcriptomic and bioinformatic analysis of all Arabidopsis thaliana SRPK, AFC, and PRP4K protein kinases to elucidate their biological functions. Unexpectedly, this revealed the existence of SRPK and AFC phylogenetic groups with distinct promoter elements and patterns of transcriptional response to abiotic stress, while PRP4Ks possess no phylogenetic sub-divisions, suggestive of functional redundancy. We also reveal splicing-related kinase families are both diel and photoperiod regulated, implicating different orthologs as discrete time-of-day RNA splicing regulators. This foundational work establishes a number of new hypotheses regarding how reversible spliceosome phosphorylation contributes to both diel plant cell regulation and abiotic stress adaptation in plants.

A crosstalk of circadian clock and alternative splicing under abiotic stresses in the plants

The circadian clock is an internal time-keeping mechanism that synchronizes the physiological adaptation of an organism to its surroundings based on day and night transition in a period of 24 h, suggesting the circadian clock provides fitness by adjusting environmental constrains. The circadian clock is driven by positive and negative elements that regulate transcriptionally and post-transcriptionally. Alternative splicing (AS) is a crucial transcriptional regulator capable of generating large numbers of mRNA transcripts from limited numbers of genes, leading to proteome diversity, which is involved in circadian to deal with abiotic stresses. Over the past decade, AS and circadian control have been suggested to coordinately regulate plant performance under fluctuating environmental conditions. However, only a few reports have reported the regulatory mechanism of this complex crosstalk. Based on the emerging evidence, this review elaborates on the existing links between circadian and AS in response to abiotic stresses, suggesting an uncovered regulatory network among circadian, AS, and abiotic stresses. Therefore, the rhythmically expressed splicing factors and core clock oscillators fill the role of temporal regulators participating in improving plant growth, development, and increasing plant tolerance against abiotic stresses.

Dynamic changes in the transcriptome landscape of Arabidopsis thaliana in response to cold stress

Plants must reprogram gene expression to adapt constantly changing environmental temperatures. With the increased occurrence of extremely low temperatures, the negative effects on plants, especially on growth and development, from cold stress are becoming more and more serious. In this research, strand-specific RNA sequencing (ssRNA-seq) was used to explore the dynamic changes in the transcriptome landscape of Arabidopsis thaliana exposed to cold temperatures (4°C) at different times. In total, 7,623 differentially expressed genes (DEGs) exhibited dynamic temporal changes during the cold treatments. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the DEGs were enriched in cold response, secondary metabolic processes, photosynthesis, glucosinolate biosynthesis, and plant hormone signal transduction pathways. Meanwhile, long non-coding RNAs (lncRNAs) were identified after the assembly of the transcripts, from which 247 differentially expressed lncRNAs (DElncRNAs) and their potential target genes were predicted. 3,621 differentially alternatively spliced (DAS) genes related to RNA splicing and spliceosome were identified, indicating enhanced transcriptome complexity due to the alternative splicing (AS) in the cold. In addition, 739 cold-regulated transcription factors (TFs) belonging to 52 gene families were identified as well. This research analyzed the dynamic changes of the transcriptome landscape in response to cold stress, which reveals more complete transcriptional patterns during short- and long-term cold treatment and provides new insights into functional studies of that how plants are affected by cold 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.

Identification of Protein Isoforms Using Reference Databases Built from Long and Short Read RNA-Sequencing

Alternative splicing can lead to distinct protein isoforms. These can have different functions in specific cells and tissues or in different developmental stages. In this study, we explored whether transcripts assembled from long read, nanopore-based, direct RNA-sequencing (RNA-seq) could improve the identification of protein isoforms in human K562 cells. By comparing with Illumina-based short read RNA-seq, we showed that a large proportion of Ensembl transcripts (5949/14,326) and genes expressing alternatively spliced transcripts (486/2981) identified with long direct reads were missed by short paired-end reads. By co-analyzing proteomic and transcriptomic data, we also showed that some peptides (826/35,976), proteins (262/3215), and protein isoforms arising from distinct transcript variants (574/1212) identified with isoform-specific peptides via custom long-read-based databases were missed in Illumina-derived databases. Finally, we generated unequivocal peptide evidence for a set of protein isoforms and showed that long read, direct RNA-seq allows the discovery of novel protein isoforms not already in reference databases or custom databases built from short read RNA-seq data. Our analysis highlights the benefits of long read RNA-seq data in the generation of reference databases to increase tandem mass spectrometry (MS/MS) identification of protein isoforms.

Ubiquitination and Ubiquitin-Like Modifications as Mediators of Alternative Pre-mRNA Splicing in Arabidopsis thaliana

Alternative splicing (AS) is a common post-transcriptional regulatory process in eukaryotes. AS has an irreplaceable role during plant development and in response to environmental stress as it evokes differential expression of downstream genes or splicing factors (e.g., serine/arginine-rich proteins). Numerous studies have reported that loss of AS capacity leads to defects in plant growth and development, and induction of stress-sensitive phenotypes. A role for post-translational modification (PTM) of AS components has emerged in recent years. These modifications are capable of regulating the activity, stability, localization, interaction, and folding of spliceosomal proteins in human cells and yeast, indicating that PTMs represent another layer of AS regulation. In this review, we summarize the recent reports concerning ubiquitin and ubiquitin-like modification of spliceosome components and analyze the relationship between spliceosome and the ubiquitin/26S proteasome pathway in plants. Based on the totality of the evidence presented, we further speculate on the roles of protein ubiquitination mediated AS in plant development and environmental response.

Combination analysis of single-molecule long-read and Illumina sequencing provides insights into the anthocyanin accumulation mechanism in an ornamental grass, Pennisetum setaceum cv. Rubrum

Combination analysis of single-molecule long-read and Illumina sequencing provide full-length transcriptome information and shed new light on the anthocyanin accumulation mechanism of Pennisetum setaceum cv. 'Rubrum'. Pennisetum setaceum cv. 'Rubrum' is an ornamental grass with purple leaves widely used in landscaping. However, the current next-generation sequencing (NGS) transcriptome information of this species is not satisfactory due to the difficulties in obtaining full-length transcripts. Furthermore, the molecular mechanisms of anthocyanin accumulation in P. setaceum have not been thoroughly studied. In this study, we used PacBio full-length transcriptome sequencing (SMRT) combined with NGS technology to build and improve the transcriptomic datasets and reveal the molecular mechanism of anthocyanin accumulation in P. setaceum cv. 'Rubrum'. Therefore, 280,413 full-length non-chimeric reads sequences were obtained using the SMRT technology. We obtained 97,450 high-quality non-redundant transcripts and identified 5352 alternative splicing events. In addition, 93,066 open reading frames (ORFs), including 57,457 full ORFs and 2910 long non-coding RNA (lncRNAs) were screened out. Furthermore, 10,795 differentially expressed genes were identified using NGS. We also explored key genes, synthesis pathways, and detected lncRNA involved in anthocyanin accumulation, providing new insights into anthocyanin accumulation in P. setaceum cv. 'Rubrum'. To our best knowledge, we provided the full-length transcriptome information of P. setaceum cv. 'Rubrum' for the first time. The results of this study will provide baseline information for gene function studies and pave the way for future P. setaceum cv. 'Rubrum' breeding projects.

Emerging Function of Ecotype-Specific Splicing in the Recruitment of Commensal Microbiome

In recent years, host-microbiome interactions in both animals and plants has emerged as a novel research area for studying the relationship between host organisms and their commensal microbial communities. The fitness advantages of this mutualistic interaction can be found in both plant hosts and their associated microbiome, however, the driving forces mediating this beneficial interaction are poorly understood. Alternative splicing (AS), a pivotal post-transcriptional mechanism, has been demonstrated to play a crucial role in plant development and stress responses among diverse plant ecotypes. This natural variation of plants also has an impact on their commensal microbiome. In this article, we review the current progress of plant natural variation on their microbiome community, and discuss knowledge gaps between AS regulation of plants in response to their intimately related microbiota. Through the impact of this article, an avenue could be established to study the biological mechanism of naturally varied splicing isoforms on plant-associated microbiome assembly.

SKIP Regulates ABA Signaling through Alternative Splicing in Arabidopsis

Abscisic acid (ABA) plays key roles in plant development and responses to abiotic stresses. A wide number of transcriptional and posttranslational regulatory mechanisms of ABA signaling are known; however, less is known about the regulatory roles of alternative splicing. In this work, we found that SKIP, a splicing factor, positively regulates ABA signaling. SKIP binds to the pre-mRNA of ABA signaling-related genes, such as PYL7, PYL8, ABI1, HAB1 and ABI5, to regulate their splicing. The precursor mRNA alternative splicing of several PYL receptors, PP2C phosphatases and ABF transcriptional factors is disrupted by the skip-1 mutation. The abnormal alternative splicing in skip-1 represses the expression of ABA-positive regulators, including PYLs and ABFs, and activates the expression of ABA-negative regulators, such as PP2Cs, which confers ABA hyposensitive phenotype of skip-1. We also found that ABA-mediated genome-wide alternative splicing and differential gene expression are changed by the skip-1 mutation. The number of the differential splicing events is increased by skip-1; however, the number of differential expressed genes in response to ABA is reduced by skip-1. Our results reveal a principle on how a splicing factor regulates ABA signaling and ABA-mediated genome-wide alternative splicing.

The Importance of a Genome-Wide Association Analysis in the Study of Alternative Splicing Mutations in Plants with a Special Focus on Maize

Alternative splicing is an important mechanism for regulating gene expressions at the post-transcriptional level. In eukaryotes, the genes are transcribed in the nucleus to produce pre-mRNAs and alternative splicing can splice a pre-mRNA to eventually form multiple different mature mRNAs, greatly increasing the number of genes and protein diversity. Alternative splicing is involved in the regulation of various plant life activities, especially the response of plants to abiotic stresses and is also an important process of plant growth and development. This review aims to clarify the usefulness of a genome-wide association analysis in the study of alternatively spliced variants by summarizing the application of alternative splicing, genome-wide association analyses and genome-wide association analyses in alternative splicing, as well as summarizing the related research progress.

ZmPP2C26 Alternative Splicing Variants Negatively Regulate Drought Tolerance in Maize

Serine/threonine protein phosphatase 2C (PP2C) dephosphorylates proteins and plays crucial roles in plant growth, development, and stress response. In this study, we characterized a clade B member of maize PP2C family, i.e., ZmPP2C26, that negatively regulated drought tolerance by dephosphorylating ZmMAPK3 and ZmMAPK7 in maize. The ZmPP2C26 gene generated ZmPP2C26L and ZmPP2C26S isoforms through untypical alternative splicing. ZmPP2C26S lost 71 amino acids including an MAPK interaction motif and showed higher phosphatase activity than ZmPP2C26L. ZmPP2C26L directly interacted with, dephosphorylated ZmMAPK3 and ZmMAPK7, and localized in chloroplast and nucleus, but ZmPP2C26S only dephosphorylated ZmMAPK3 and localized in cytosol and nucleus. The expression of ZmPP2C26L and ZmPP2C26 was significantly inhibited by drought stress. Meanwhile, the maize zmpp2c26 mutant exhibited enhancement of drought tolerance with higher root length, root weight, chlorophyll content, and photosynthetic rate compared with wild type. However, overexpression of ZmPP2C26L and ZmPP2C26S significantly decreased drought tolerance in Arabidopsis and rice with lower root length, chlorophyll content, and photosynthetic rate. Phosphoproteomic analysis revealed that the ZmPP2C26 protein also altered phosphorylation level of proteins involved in photosynthesis. This study provides insights into understanding the mechanism of PP2C in response to abiotic stress.

ABA Mediates Plant Development and Abiotic Stress via Alternative Splicing

Alternative splicing (AS) exists in eukaryotes to increase the complexity and adaptability of systems under biophysiological conditions by increasing transcriptional and protein diversity. As a classic hormone, abscisic acid (ABA) can effectively control plant growth, improve stress resistance, and promote dormancy. At the transcriptional level, ABA helps plants respond to the outside world by regulating transcription factors through signal transduction pathways to regulate gene expression. However, at the post-transcriptional level, the mechanism by which ABA can regulate plant biological processes by mediating alternative splicing is not well understood. Therefore, this paper briefly introduces the mechanism of ABA-induced alternative splicing and the role of ABA mediating AS in plant response to the environment and its own growth.