Alternative splicing of pre-messenger RNAs in plants in the genomic era

Diversification and evolution of serine/arginine-rich splicing factors in plants

In plants, serine/arginine-rich splicing factors (SRSFs) play an essential role in pre-mRNA splicing, influencing growth and stress response. The present study investigated the evolution and diversification of SR proteins in plants by phylogenetic analysis. First, we identified 1650 SR proteins from 132 plant species. Large-scale phylogenetic analysis divided these SR proteins into nine subfamilies: SR, RS, SC, SCL, RSZ, RS2Z, SR45, SR45a, and a newly discovered RSZ-like subfamily. The analysis also suggested that SR proteins of animals and plants originated from an ancient SR gene of a common ancestor of eukaryotes and showed that RSZ-like subfamily diverged in fern. We also found the presence of the newly discovered RSZ-like subfamily in ferns, gymnosperms, and angiosperms as a sister group of SR45 and SR45a, indicating a common ancestor. In addition, we detected deletion of the N-terminal RS domain in a few members of SR45 and SR45a of plants. These findings provide novel insights into the phylogeny of the SR gene family in plants, emphasizing their crucial function as a splicing factor.

The splicing auxiliary factor OsU2AF35a enhances thermotolerance via protein separation and promoting proper splicing of OsHSA32 pre-mRNA in rice

Heat stress significantly impacts global rice production, highlighting the critical need to understand the genetic basis of heat resistance in rice. U2AF (U2 snRNP auxiliary factor) is an essential splicing complex with critical roles in recognizing the 3'-splice site of precursor messenger RNAs (pre-mRNAs). The U2AF small subunit (U2AF35) can bind to the 3'-AG intron border and promote U2 snRNP binding to the branch-point sequences of introns through interaction with the U2AF large subunit (U2AF65). However, the functions of U2AF35 in plants are poorly understood. In this study, we discovered that the OsU2AF35a gene was vigorously induced by heat stress and could positively regulate rice thermotolerance during both the seedling and reproductive growth stages. OsU2AF35a interacts with OsU2AF65a within the nucleus, and both of them can form condensates through liquid-liquid phase separation (LLPS) following heat stress. The intrinsically disordered regions (IDR) are accountable for their LLPS. OsU2AF35a condensation is indispensable for thermotolerance. RNA-seq analysis disclosed that, subsequent to heat treatment, the expression levels of several genes associated with water deficiency and oxidative stress in osu2af35a-1 were markedly lower than those in ZH11. In accordance with this, OsU2AF35a is capable of positively regulating the oxidative stress resistance of rice. The pre-mRNAs of a considerable number of genes in the osu2af35a-1 mutant exhibited defective splicing, among which was the OsHSA32 gene. Knocking out OsHSA32 significantly reduced the thermotolerance of rice, while overexpressing OsHSA32 could partially rescue the heat sensitivity of osu2af35a-1. Together, our findings uncovered the essential role of OsU2AF35a in rice heat stress response through protein separation and regulating alternative pre-mRNA splicing.

Nanopore-based full-length transcriptome sequencing of the skin in Pseudopleuronectes yokohamae identifies novel antimicrobial peptide genes

The marbled flounder (Pseudopleuronectes yokohamae) is highly esteemed for its exceptional nutritional value and delicious taste. However, this species has extremely limited transcriptome data, which can offer priceless information for disease protection. In the study, the skin transcriptomic sequencing of P. yokohamae revealed 7.72 GB of clean data using the Nanopore sequencing platform. The results revealed 30,498 transcripts of functional annotations in the P. yokohamae transcriptome. All transcripts were searched in eight functional databases. A total of 10,337 ORFs were obtained, of which 6081 complete ORFs accounted for 58.83% of all predicted CDS. Moreover, 10,195 SSRs were detected. Meanwhile, the non-pecific immunity pathways were investigated for better understanding of the immunological reaction in P. yokohamae, and seven innate immune pathways were identified. The innate-immune related genes were highly expressed in the NOD-like receptor signaling pathway, followed by the C-type lectin receptor signaling pathway, Toll-like receptor signaling pathway, RIG-I-like receptor signaling pathway and Cytosolic DNA-sensing pathway. In this study, four families of antimicrobial peptides (AMPs) in P. yokohamae were analysed for the first time, including piscidins, hepcidins, liver-expressed antimicrobial peptide and defensins. Seven AMPs, including Pypleurocidin-like WF3, Pypleurocidin-like WFX, Pyhepcidin 1, Pyhepcidin-like 1, PyLEAP-2, Pybeta-defensin and Pybeta-defensin-like 1, were further identified. The seven AMPs showed a highly identity in their cDNA and genomic structures and an inducible expression pattern preferable to skin in response to pathogens. The transcriptomic data and investigation of AMPs from P. yokohamae promote a deeper awareness of fish mucosal immunity and provide information in the prevention of fish diseases.

Genome-Scale Identification of Wild Soybean Serine/Arginine-Rich Protein Family Genes and Their Responses to Abiotic Stresses

Serine/arginine-rich (SR) proteins mostly function as splicing factors for pre-mRNA splicing in spliceosomes and play critical roles in plant development and adaptation to environments. However, detailed study about SR proteins in legume plants is still lacking. In this report, we performed a genome-wide investigation of SR protein genes in wild soybean (Glycine soja) and identified a total of 31 GsSR genes from the wild soybean genome. The analyses of chromosome location and synteny show that the GsSRs are unevenly distributed on 15 chromosomes and are mainly under the purifying selection. The GsSR proteins can be phylogenetically classified into six sub-families and are conserved in evolution. Prediction of protein phosphorylation sites indicates that GsSR proteins are highly phosphorylated proteins. The protein-protein interaction network implies that there exist numerous interactions between GsSR proteins. We experimentally confirmed their physical interactions with the representative SR proteins of spliceosome-associated components such as U1-70K or U2AF35 by yeast two-hybrid assays. In addition, we identified various stress-/hormone-responsive cis-acting elements in the promoter regions of these GsSR genes and verified their expression patterns by RT-qPCR analyses. The results show most GsSR genes are highly expressed in root and stem tissues and are responsive to salt and alkali stresses. Splicing analysis showed that the splicing patterns of GsSRs were in a tissue- and stress-dependent manner. Overall, these results will help us to further investigate the biological functions of leguminous plant SR proteins and shed new light on uncovering the regulatory mechanisms of plant SR proteins in growth, development, and stress responses.

A global survey of bicarbonate stress-induced pre-mRNA alternative splicing in soybean via integrative analysis of Iso-seq and RNA-seq

Alternative splicing (AS) plays important roles in modulating environmental stress responses in plants. However, little is known about the functions of bicarbonate-induced AS in cultivated soybean (Glycine max L. Merr.). In this study, we combined PacBio isoform sequencing (Iso-seq) and Illumina RNA sequencing (RNA-seq) to elucidate the bicarbonate-induced AS events in soybean root and leaf tissues. Compared to RNA-seq, Iso-seq identified more novel genes and transcripts, as well as more AS events, indicating that Iso-seq is more efficient in AS detection. Combining these two technologies, we found that intron retention (IR) is the most frequent AS event type. We identified a total of 913 and 1974 bicarbonate stress-responsive differentially alternative spliced genes (DAGs) in soybean leaves and roots respectively, from our RNA-seq results. Additionally, we determined a transcription factor (GmNTL9) and a splicing factor (GmRSZ22), and validated their roles in bicarbonate stress response by AS. Overall, our study opens an avenue for evaluating plant AS regulatory networks, and the obtained global landscape of alternative splicing provides valuable insights into the AS-mediated bicarbonate-responsive mechanisms in plant species.

Rare occurrence of cryptic 5' splice sites by downstream 3' splice site/exon boundary mutations in a heavy-ion-induced egy1-4 allele of Arabidopsis thaliana

Pre-mRNA splicing is a fundamental process in eukaryotic gene expression, and the mechanism of intron definition, involving the recognition of the canonical GU (5'-splice site) and AG (3'-splice site) dinucleotides by splicing factors, has been postulated for most cases of splicing initiation in plants. Splice site mutations have played crucial roles in unraveling the mechanism of pre-mRNA splicing in planta. Typically, splice site mutations abolish splicing events or activate one or more cryptic splice sites surrounding the mutated region. In this report, we investigated the splicing pattern of the EGY1 gene in an Ar-ion-induced egy1-4 allele of Arabidopsis thaliana. egy1-4 has an AG-to-AC mutation in the 3'-end of intron 3, along with 4-bp substitutions and a 5-bp deletion in adjacent exon 4. RT-PCR, cDNA cloning, and amplicon sequencing analyses of EGY1 revealed that while most wild-type EGY1 mRNAs had a single splicing pattern, egy1-4 mRNAs had multiple splicing defects. Almost half of EGY1 transcripts showed 'intron retention' at intron 3, while the other half exhibited activation of 3' cryptic splice sites either upstream or downstream of the original 3'-splice site. Unexpectedly, around 8% of EGY1 transcripts in egy1-4 exhibited activation of cryptic 5'-splice sites positioned upstream of the authentic 5'-splice site of intron 3. Whole genome resequencing of egy1-4 indicated that it has no other known impactful mutations. These results may provide a rare, but real case of activation of cryptic 5'-splice sites by downstream 3'-splice site/exon mutations in planta.

Alternative splicing of a potato disease resistance gene maintains homeostasis between growth and immunity

Plants possess a robust and sophisticated innate immune system against pathogens and must balance growth with rapid pathogen detection and defense. The intracellular receptors with nucleotide-binding leucine-rich repeat (NLR) motifs recognize pathogen-derived effector proteins and thereby trigger the immune response. The expression of genes encoding NLR receptors is precisely controlled in multifaceted ways. The alternative splicing (AS) of introns in response to infection is recurrently observed but poorly understood. Here we report that the potato (Solanum tuberosum) NLR gene RB undergoes AS of its intron, resulting in 2 transcriptional isoforms, which coordinately regulate plant immunity and growth homeostasis. During normal growth, RB predominantly exists as an intron-retained isoform RB_IR, encoding a truncated protein containing only the N-terminus of the NLR. Upon late blight infection, the pathogen induces intron splicing of RB, increasing the abundance of RB_CDS, which encodes a full-length and active R protein. By deploying the RB splicing isoforms fused with a luciferase reporter system, we identified IPI-O1 (also known as Avrblb1), the RB cognate effector, as a facilitator of RB AS. IPI-O1 directly interacts with potato splicing factor StCWC15, resulting in altered localization of StCWC15 from the nucleoplasm to the nucleolus and nuclear speckles. Mutations in IPI-O1 that eliminate StCWC15 binding also disrupt StCWC15 re-localization and RB intron splicing. Thus, our study reveals that StCWC15 serves as a surveillance facilitator that senses the pathogen-secreted effector and regulates the trade-off between RB-mediated plant immunity and growth, expanding our understanding of molecular plant-microbe interactions.

Relevance and regulation of alternative splicing in plant secondary metabolism: current understanding and future directions

The secondary metabolism of plants is an essential life process enabling organisms to navigate various stages of plant development and cope with ever-changing environmental stresses. Secondary metabolites, abundantly found in nature, possess significant medicinal value. Among the regulatory mechanisms governing these metabolic processes, alternative splicing stands out as a widely observed post-transcriptional mechanism present in multicellular organisms. It facilitates the generation of multiple mRNA transcripts from a single gene by selecting different splicing sites. Selective splicing events in plants are widely induced by various signals, including external environmental stress and hormone signals. These events ultimately regulate the secondary metabolic processes and the accumulation of essential secondary metabolites in plants by influencing the synthesis of primary metabolites, hormone metabolism, biomass accumulation, and capillary density. Simultaneously, alternative splicing plays a crucial role in enhancing protein diversity and the abundance of the transcriptome. This paper provides a summary of the factors inducing alternative splicing events in plants and systematically describes the progress in regulating alternative splicing with respect to different secondary metabolites, including terpenoid, phenolic compounds, and nitrogen-containing compounds. Such elucidation offers critical foundational insights for understanding the role of alternative splicing in regulating plant metabolism and presents novel avenues and perspectives for bioengineering.

AtSNU13 modulates pre-mRNA splicing of RBOHD and ALD1 to regulate plant immunity

Pre-mRNA splicing is a significant step for post-transcriptional modifications and functions in a wide range of physiological processes in plants. Human NHP2L binds to U4 snRNA during spliceosome assembly; it is involved in RNA splicing and mediates the development of human tumors. However, no ortholog has yet been identified in plants. Therefore, we report At4g12600 encoding the ortholog NHP2L protein, and AtSNU13 associates with the component of the spliceosome complex; the atsnu13 mutant showed compromised resistance in disease resistance, indicating that AtSNU13 is a positive regulator of plant immunity. Compared to wild-type plants, the atsnu13 mutation resulted in altered splicing patterns for defense-related genes and decreased expression of defense-related genes, such as RBOHD and ALD1. Further investigation shows that AtSNU13 promotes the interaction between U4/U6.U5 tri-snRNP-specific 27 K and the motif in target mRNAs to regulate the RNA splicing. Our study highlights the role of AtSNU13 in regulating plant immunity by affecting the pre-mRNA splicing of defense-related genes.

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.

An alternative 3' splice site of PeuHKT1;3 improves the response to salt stress through enhancing affinity to K+ in Populus

Alternative splicing (AS) serves as a crucial post-transcriptional regulator in plants that contributes to the resistance to salt stress. However, the underlying mechanism is largely unknown. In this research, we identified an important AS transcript in Populus euphratica, PeuHKT1:3a, generated by alternative 3' splice site splicing mode that resulted in the removal of 252 bases at the 5' end of the first exon in PeuHKT1:3. Protein sequence comparison showed that the site of AS occurred in PeuHKT1:3 is located at a crucial Ser residue within the first pore-loop domain, which leads to inefficient K+ transport in HKT I-type transporters. Expressing PeuHKT1;3a in an axt3 mutant yeast strain can effectively compensate for the lack of intracellular K+, whereas the expression of PeuHKT1;3 cannot yield the effect. Furthermore, in transgenic Arabidopsis and poplar plants, it was observed that lines expressing PeuHKT1;3a exhibited greater salt tolerance compared to those expressing the PeuHKT1;3 strain. Analysis of ion content and flux demonstrated that the transgenic PeuHKT1;3a line exhibited significantly higher K+ content compared to the PeuHKT1;3 line, while there was no significant difference in Na+ content. In conclusion, our findings revealed that AS can give rise to novel variants of HKT I-type proteins in P. euphratica with modified K+ selectivity to keep a higher K+/Na+ ratio to enhanced salt tolerance.

Positioning of pyrimidine motifs around cassette exons defines their PTB-dependent splicing in Arabidopsis

Alternative splicing (AS) is a complex process that generates transcript variants from a single pre-mRNA and is involved in numerous biological functions. Many RNA-binding proteins are known to regulate AS; however, little is known about the underlying mechanisms, especially outside the mammalian clade. Here, we show that polypyrimidine tract binding proteins (PTBs) from Arabidopsis thaliana regulate AS of cassette exons via pyrimidine (Py)-rich motifs close to the alternative splice sites. Mutational studies on three PTB-dependent cassette exon events revealed that only some of the Py motifs in this region are critical for AS. Moreover, in vitro binding of PTBs did not reflect a motif's impact on AS in vivo. Our mutational studies and bioinformatic investigation of all known PTB-regulated cassette exons from A. thaliana and human suggested that the binding position of PTBs relative to a cassette exon defines whether its inclusion or skipping is induced. Accordingly, exon skipping is associated with a higher frequency of Py stretches within the cassette exon, and in human also upstream of it, whereas exon inclusion is characterized by increased Py motif occurrence downstream of said exon. Enrichment of Py motifs downstream of PTB-activated 5' splice sites is also seen for PTB-dependent intron removal and alternative 5' splice site events from A. thaliana, suggesting this is a common step of exon definition. In conclusion, the position-dependent AS regulatory mechanism by PTB homologs has been conserved during the separate evolution of plants and mammals, while other critical features, in particular intron length, have considerably changed.

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

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

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

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.

Quantitative Phosphoproteomic Analysis Provides Insights into the Sodium Bicarbonate Responsiveness of Glycine max

Sodium bicarbonate stress caused by NaHCO3 is one of the most severe abiotic stresses affecting agricultural production worldwide. However, little attention has been given to the molecular mechanisms underlying plant responses to sodium bicarbonate stress. To understand phosphorylation events in signaling pathways triggered by sodium bicarbonate stress, TMT-labeling-based quantitative phosphoproteomic analyses were performed on soybean leaf and root tissues under 50 mM NaHCO3 treatment. In the present study, a total of 7856 phosphopeptides were identified from cultivated soybeans (Glycine max L. Merr.), representing 3468 phosphoprotein groups, in which 2427 phosphoprotein groups were newly identified. These phosphoprotein groups contained 6326 unique high-probability phosphosites (UHPs), of which 77.2% were newly identified, increasing the current soybean phosphosite database size by 43.4%. Among the phosphopeptides found in this study, we determined 67 phosphopeptides (representing 63 phosphoprotein groups) from leaf tissue and 554 phosphopeptides (representing 487 phosphoprotein groups) from root tissue that showed significant changes in phosphorylation levels under sodium bicarbonate stress (fold change >1.2 or <0.83, respectively; p < 0.05). Localization prediction showed that most phosphoproteins localized in the nucleus for both leaf and root tissues. GO and KEGG enrichment analyses showed quite different enriched functional terms between leaf and root tissues, and more pathways were enriched in the root tissue than in the leaf tissue. Moreover, a total of 53 different protein kinases and 7 protein phosphatases were identified from the differentially expressed phosphoproteins (DEPs). A protein kinase/phosphatase interactor analysis showed that the interacting proteins were mainly involved in/with transporters/membrane trafficking, transcriptional level regulation, protein level regulation, signaling/stress response, and miscellaneous functions. The results presented in this study reveal insights into the function of post-translational modification in plant responses to sodium bicarbonate stress.

The Abundant and Unique Transcripts and Alternative Splicing of the Artificially Autododecaploid London Plane (Platanus × acerifolia)

Transcription and alternative splicing (AS) are now appreciated in plants, but few studies have examined the effects of changing ploidy on transcription and AS. In this study, we showed that artificially autododecaploid plants of London plane (Platanus × acerifolia (Aiton) Willd) had few flowers relative to their hexaploid progenitors. Transcriptome analysis based on full-length Oxford Nanopore Technologies (ONTs) and next-generation sequencing (NGS) revealed that the increased ploidy level in P. × acerifolia led to more transcribed isoforms, accompanied by an increase in the number of isoforms per gene. The functional enrichment of genes indicated that novel genes transcribed specifically in the dodecaploids may have been highly correlated with the ability to maintain genome stability. The dodecaploids showed a higher number of genes with upregulated differentially expressed genes (DEGs) compared with the hexaploid counterpart. The genome duplication of P. × acerifolia resulted mainly in the DEGs involved in basic biological pathways. It was noted that there was a greater abundance of alternative splicing (AS) events and AS genes in the dodecaploids compared with the hexaploids in P. × acerifolia. In addition, a significant difference between the structure and expression of AS events between the hexaploids and dodecaploids of Platanus was found. Of note, some DEGs and differentially spliced genes (DSGs) related to floral transition and flower development were consistent with the few flower traits in the dodecaploids of P. × acerifolia. Collectively, our findings explored the difference in transcription and AS regulation between the hexaploids and dodecaploids of P. × acerifolia and gained new insight into the molecular mechanisms underlying the few-flower phenotype of P. × acerifolia. These results contribute to uncovering the regulatory role of transcription and AS in polyploids and breeding few-flower germplasms.

Sequence variations affect the 5' splice site selection of plant introns

Introns are noncoding sequences spliced out of pre-mRNAs by the spliceosome to produce mature mRNAs. The 5' ends of introns mostly begin with GU and have a conserved sequence motif of AG/GUAAGU that could base-pair with the core sequence of U1 snRNA of the spliceosome. Intriguingly, ∼ 1% of introns in various eukaryotic species begin with GC. This occurrence could cause misannotation of genes; however, the underlying splicing mechanism is unclear. We analyzed the sequences around the intron 5' splice site (ss) in Arabidopsis (Arabidopsis thaliana) and found sequences at the GC intron ss are much more stringent than those of GT introns. Mutational analysis at various positions of the intron 5' ss revealed that although mutations impair base pairing, different mutations at the same site can have different effects, suggesting that steric hindrance also affects splicing. Moreover, mutations of 5' ss often activate a hidden ss nearby. Our data suggest that the 5' ss is selected via a competition between the major ss and the nearby minor ss. This work not only provides insights into the splicing mechanism of intron 5' ss but also improves the accuracy of gene annotation and the study of the evolution of intron 5' ss.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

Catsnap: a user-friendly algorithm for determining the conservation of protein variants reveals extensive parallelisms in the evolution of alternative splicing

Understanding the evolutionary conservation of complex eukaryotic transcriptomes significantly illuminates the physiological relevance of alternative splicing (AS). Examining the evolutionary depth of a given AS event with ordinary homology searches is generally challenging and time-consuming. Here, we present Catsnap, an algorithmic pipeline for assessing the conservation of putative protein isoforms generated by AS. It employs a machine learning approach following a database search with the provided pair of protein sequences. We used the Catsnap algorithm for analyzing the conservation of emerging experimentally characterized alternative proteins from plants and animals. Indeed, most of them are conserved among other species. Catsnap can detect the conserved functional protein isoforms regardless of the AS type by which they are generated. Notably, we found that while the primary amino acid sequence is maintained, the type of AS determining the inclusion or exclusion of protein regions varies throughout plant phylogenetic lineages in these proteins. We also document that this phenomenon is less seen among animals. In sum, our algorithm highlights the presence of unexpectedly frequent hotspots where protein isoforms recurrently arise to carry physiologically relevant functions. The user web interface is available at https://catsnap.cesnet.cz/.

The CmMYB3 transcription factors isolated from the Chrysanthemum morifolium regulate flavonol biosynthesis in Arabidopsis thaliana

Chrysanthemum morifolium MYB3 factors are transcriptional activators for the regulation of flavonol biosynthesis. Flavonol was not only the critical secondary metabolite participating in the growth and development of plants but also the main active ingredient in medicinal chrysanthemum. However, few pieces of research revealed the transcriptional regulation of flavonol biosynthesis in Chrysanthemum morifolium. Here, we isolated two CmMYB3 transcription factors (CmMYB3a and CmMYB3b) from the capitulum of Chrysanthemum morifolium cv 'Hangju'. According to the sequence characteristics, the CmMYB3a and CmMYB3b belonged to the R2R3-MYB subgroup 7, whose members were often reported to regulate flavonol biosynthesis positively. CmMYB3a and CmMYB3b factors were identified to localize in the nucleus by subcellular localization assay. Besides, both of them have obvious transcriptional self-activation activity in their C-terminal. After the overexpression of CmMYB3 genes in Nicotiana benthamiana and Arabidopsis thaliana, the flavonol contents in plants were increased, and the expression of AtCHS, AtCHI, AtF3H, and AtFLS genes in A. thaliana was also improved. Interestingly, the CmMYB3a factor had stronger functions in improving flavonol contents and related gene expression levels than CmMYB3b. The interaction analysis between transcription factors and promoters suggested that CmMYB3 could bind and activate the promoters of CmCHI and CmFLS genes in C. morifolium, and CmMYB3a also functioned more powerfully. Overall, these results indicated that CmMYB3a and CmMYB3b work as transcriptional activators in controlling flavonol biosynthesis.

Evidence for the role of transcription factors in the co-transcriptional regulation of intron retention

BACKGROUND

Alternative splicing is a widespread regulatory phenomenon that enables a single gene to produce multiple transcripts. Among the different types of alternative splicing, intron retention is one of the least explored despite its high prevalence in both plants and animals. The recent discovery that the majority of splicing is co-transcriptional has led to the finding that chromatin state affects alternative splicing. Therefore, it is plausible that transcription factors can regulate splicing outcomes.

RESULTS

We provide evidence for the hypothesis that transcription factors are involved in the regulation of intron retention by studying regions of open chromatin in retained and excised introns. Using deep learning models designed to distinguish between regions of open chromatin in retained introns and non-retained introns, we identified motifs enriched in IR events with significant hits to known human transcription factors. Our model predicts that the majority of transcription factors that affect intron retention come from the zinc finger family. We demonstrate the validity of these predictions using ChIP-seq data for multiple zinc finger transcription factors and find strong over-representation for their peaks in intron retention events.

CONCLUSIONS

This work opens up opportunities for further studies that elucidate the mechanisms by which transcription factors affect intron retention and other forms of splicing.

AVAILABILITY

Source code available at https://github.com/fahadahaf/chromir.

Alternative splicing in ABA signaling during seed germination

Seed germination is an essential step in a plant's life cycle. It is controlled by complex physiological, biochemical, and molecular mechanisms and external factors. Alternative splicing (AS) is a co-transcriptional mechanism that regulates gene expression and produces multiple mRNA variants from a single gene to modulate transcriptome diversity. However, little is known about the effect of AS on the function of generated protein isoforms. The latest reports indicate that alternative splicing (AS), the relevant mechanism controlling gene expression, plays a significant role in abscisic acid (ABA) signaling. In this review, we present the current state of the art about the identified AS regulators and the ABA-related changes in AS during seed germination. We show how they are connected with the ABA signaling and the seed germination process. We also discuss changes in the structure of the generated AS isoforms and their impact on the functionality of the generated proteins. Also, we point out that the advances in sequencing technology allow for a better explanation of the role of AS in gene regulation by more accurate detection of AS events and identification of full-length splicing isoforms.

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

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

Comprehensive profiling of epigenetic modifications in fast-growing Moso bamboo shoots

The fast growth of Moso bamboo (Phyllostachys edulis) shoots is caused by the rapid elongation of each internode. However, the key underlying cellular processes and epigenetic mechanisms remain largely unexplored. We used microscopy and multi-omics approaches to investigate two regions (bottom and middle) of the 18th internode from shoots of two different heights (2 and 4 m). We observed that internode cells become longer, and that lignin biosynthesis and glycosyltransferase family 43 (GT43) genes are substantially upregulated with shoot height. Nanopore direct RNA sequencing (DRS) revealed a higher N6-methyladenine (m6A) modification rate in 2-m shoots than in 4-m shoots. In addition, different specific m6A modification sites were enriched at different growth stages. Global DNA methylation profiling indicated that DNA methylation levels are higher in 4-m shoots than in 2-m shoots. We also detected shorter poly(A) tail lengths (PALs) in 4-m shoots compared with 2-m shoots. Genes showing differential PAL were mainly enriched in the functional terms of protein translation and vesicle fusion. An association analysis between PALs and DNA methylation strongly suggested that gene body CG methylation levels are positively associated with PAL. This study provides valuable information to better understand post-transcriptional regulations responsible for fast-growing shoots in Moso bamboo.

Alternative Splicing of TaHsfA2-7 Is Involved in the Improvement of Thermotolerance in Wheat

High temperature has severely affected plant growth and development, resulting in reduced production of crops worldwide, especially wheat. Alternative splicing (AS), a crucial post-transcriptional regulatory mechanism, is involved in the growth and development of eukaryotes and the adaptation to environmental changes. Previous transcriptome data suggested that heat shock transcription factor (Hsf) TaHsfA2-7 may form different transcripts by AS. However, it remains unclear whether this post-transcriptional regulatory mechanism of TaHsfA2-7 is related to thermotolerance in wheat (Triticum aestivum). Here, we identified a novel splice variant, TaHsfA2-7-AS, which was induced by high temperature and played a positive role in thermotolerance regulation in wheat. Moreover, TaHsfA2-7-AS is predicted to encode a small truncated TaHsfA2-7 isoform, retaining only part of the DNA-binding domain (DBD). TaHsfA2-7-AS is constitutively expressed in various tissues of wheat. Notably, the expression level of TaHsfA2-7-AS is significantly up-regulated by heat shock (HS) during flowering and grain-filling stages in wheat. Further studies showed that TaHsfA2-7-AS was localized in the nucleus but lacked transcriptional activation activity. Ectopic expression of TaHsfA2-7-AS in yeast exhibited improved thermotolerance. Compared to non-transgenic plants, overexpression of TaHsfA2-7-AS in Arabidopsis results in enhanced tolerance to heat stress. Simultaneously, we also found that TaHsfA1 is directly involved in the transcriptional regulation of TaHsfA2-7 and TaHsfA2-7-AS. In summary, our findings demonstrate the function of TaHsfA2-7-AS splicing variant in response to heat stress and establish a link between regulatory mechanisms of AS and the improvement of thermotolerance in wheat.

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.

The RNA helicase UAP56 and the E3 ubiquitin ligase COP1 coordinately regulate alternative splicing to repress photomorphogenesis in Arabidopsis

Light is a key environmental signal that regulates plant growth and development. While posttranscriptional regulatory mechanisms of gene expression include alternative splicing (AS) of pre-messenger RNA (mRNA) in both plants and animals, how light signaling affects AS in plants is largely unknown. Here, we identify DExD/H RNA helicase U2AF65-associated protein (UAP56) as a negative regulator of photomorphogenesis in Arabidopsis thaliana. UAP56 is encoded by the homologs UAP56a and UAP56b. Knockdown of UAP56 led to enhanced photomorphogenic responses and diverse developmental defects during vegetative and reproductive growth. UAP56 physically interacts with the central light signaling repressor constitutive photomorphogenic 1 (COP1) and U2AF65. Global transcriptome analysis revealed that UAP56 and COP1 co-regulate the transcription of a subset of genes. Furthermore, deep RNA-sequencing analysis showed that UAP56 and COP1 control pre-mRNA AS in both overlapping and distinct manners. Ribonucleic acid immunoprecipitation assays showed that UAP56 and COP1 bind to common small nuclear RNAs and mRNAs of downstream targets. Our study reveals that both UAP56 and COP1 function as splicing factors that coordinately regulate AS during light-regulated plant growth and development.

Transcriptome and Proteome Co-Profiling Offers an Understanding of Pre-Harvest Sprouting (PHS) Molecular Mechanisms in Wheat (Triticum aestivum)

While wheat (Triticum aestivum L.) is a widely grown and enjoyed crop, the diverse and complex global situation and climate are exacerbating the instability of its supply. In particular, pre-harvest sprouting (PHS) is one of the major abiotic stresses that frequently occurs due to irregular climate conditions, causing serious damage to wheat and its quality. In this study, transcriptomic analysis with RNA-seq and proteomic analysis with LC-MS/MS were performed in PHS-treated spikes from two wheat cultivars presenting PHS sensitivity and tolerance, respectively. A total of 13,154 differentially expressed genes (DEGs) and 706 differentially expressed proteins (DEPs) were identified in four comparison groups between the susceptible/tolerant cultivars. Gene function and correlation analysis were performed to determine the co-profiled genes and proteins affected by PHS treatment. In the functional annotation of each comparative group, similar functions were confirmed in each cultivar under PHS treatment; however, in Keumgang PHS+7 (K7) vs. Woori PHS+7 (W7), functional annotations presented clear differences in the "spliceosome" and "proteasome" pathways. In addition, our results indicate that alternative splicing and ubiquitin-proteasome support the regulation of germination and seed dormancy. This study provides an advanced understanding of the functions involved in transcription and translation related to PHS mechanisms, thus enabling specific proposals for the further analysis of germination and seed dormancy mechanisms and pathways in wheat.

A novel full-length transcriptome resource from multiple immune-related tissues in turbot (Scophthalmus maximus) using Pacbio SMART sequencing

Turbot (Scophthalmus maximus) is an important cold-water economic fish. However, the production and development of turbot industry has been constantly hindered by the frequent occurrence of some diseases. Lacking full-length transcriptome for turbot limits immune gene discoveries and gene structures analysis. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of turbot with PacBio Sequel platform. In this study, a total of 31.7 Gb high quality data were generated with the average subreads length of 2618 bp. According to the presence of 5' and 3' primers as well as poly (A) tails, FL (Full-length) and NFL (Non-full-length) isoforms were obtained. Meanwhile, we identified 32,003 non-redundant transcripts, 76.02% of which was novel isoforms of known genes. In addition, 12,176 alternative splicing (AS) events, 6614 polyadenylation (APA) events, 1905 transcription factors, and 2703 lncRNAs were identified. This work is a comprehensive report on the full-length transcriptome of immune-related tissues of turbot, and it also provides valuable molecular resources for future research on the adaptation mechanisms and functional genomics of turbot.

Transcriptomic and alternative splicing analyses provide insights into the roles of exogenous salicylic acid ameliorating waxy maize seedling growth under heat stress

BACKGROUND

Salicylic acid (SA) is a phytohormone which works to regulate the abiotic stress response of plants. However, the molecular mechanism by which SA mediates heat tolerance in waxy maize (Zea mays L. sinsensis Kulesh) remains unknown.

RESULTS

Two varieties of waxy maize seedlings, heat-tolerant 'Yunuo7' (Y7) and heat-sensitive 'Suyunuo5' (S5), were pretreated with SA prior to heat stress (HTS). After treatment, physiological and transcriptomic changes were analyzed. Compared with HTS, the exogenous application of SA enhanced the shoot dry weight, the activities of antioxidant enzymes (e.g., SOD, POD, CAT and APX), and the concentration of endogenous phytohormones (e.g., SA, ABA, IAA, GA3), while decreased the MDA content. Transcriptome analysis showed that the number of differentially expressed genes (DEGs) identified in the control (CK) vs HTS and HTS vs HTS + SA comparisons were more in S5 than in Y7. HTS induced the downregulation of genes involved in photosynthesis and the upregulation of genes encoding heat shock transcription factors (HSFs) and heat shock proteins (HSPs). Compared with HTS, SA pretreatment reversed the expression of 5 photosynthesis-related genes, 26 phytohormone-related genes, and all genes encoding HSFs and HSPs in S5. Furthermore, the number of alternative splicing (AS) events increased under HTS treatment for both varieties, while decreased under SA pretreatment of S5. Differentially spliced genes (DSGs) showed little overlap with DEGs, and DEGs and DSGs differed significantly in functional enrichment.

CONCLUSIONS

Physiological and transcriptional together indicated that HTS and SA pretreatment had a greater effect on S5 than Y7. Additionally, it appears that transcriptional regulation and AS work synergistically to enhance thermotolerance in heat-sensitive waxy maize. Our study revealed the regulatory effects and underlying molecular mechanisms of SA on waxy maize seedling under HTS.

Alternative splicing of Arabidopsis G6PD5 recruits NADPH-producing OPPP reactions to the endoplasmic reticulum

Glucose-6-phosphate dehydrogenase is the rate-limiting enzyme of the oxidative pentose-phosphate pathway (OPPP). The OPPP mainly provides NADPH and sugar-phosphate building blocks for anabolic pathways and is present in all eukaryotes. In plant cells, the irreversible part of the OPPP is found in several compartments. Among the isoforms catalyzing the first OPPP step in Arabidopsis, G6PD1 to G6PD4 target plastids (with G6PD1 being also directed to peroxisomes), whereas G6PD5 and G6PD6 operate in the cytosol. We noticed that alternative splice forms G6PD5.4 and G6PD5.5 encode N-terminally extended proteoforms. Compared to G6PD5.1, RT-PCR signals differed and fluorescent reporter fusions expressed in Arabidopsis protoplasts accumulated in distinct intracellular sites. Co-expression with organelle-specific markers revealed that the G6PD5.4 and G6PD5.5 proteoforms label different subdomains of the endoplasmic reticulum (ER), and analysis of C-terminal roGFP fusions showed that their catalytic domains face the cytosol. In g6pd5-1 g6pd6-2 mutant protoplasts lacking cytosolic G6PDH activity, the ER-bound proteoforms were both active and thus able to form homomers. Among the Arabidopsis 6-phosphogluconolactonases (catalyzing the second OPPP step), we noticed that isoform PGL2 carries a C-terminal CaaX motif that may be prenylated for membrane attachment. Reporter-PGL2 fusions co-localized with G6PD5.4 in ER subdomains, which was abolished by Cys-to-Ser exchange in the ²⁵⁶CSIL motif. Among the Arabidopsis 6-phosphogluconate dehydrogenases (catalyzing the third OPPP step), S-acylated peptides were detected for all three isoforms in a recent palmitoylome, with dual cytosolic/peroxisomal PGD2 displaying three sites. Co-expression of GFP-PGD2 diminished crowding of OFP-G6PD5.4 at the ER, independent of PGL2's presence. Upon pull-down of GFP-G6PD5.4, not only unlabeled PGD2 and PGL2 were enriched, but also enzymes that depend on NADPH provision at the ER, indicative of physical interaction with the OPPP enzymes. When membrane-bound G6PD5.5 and 5.4 variants were co-expressed with KCR1 (ketoacyl-CoA reductase, involved in fatty acid elongation), ATR1 (NADPH:cytochrome-P450 oxidoreductase), or pulled C4H/CYP73A5 (cinnamate 4-hydroxylase) as indirectly (via ATR) NADPH-dependent cytochrome P450 enzyme, co-localization in ER subdomains was observed. Thus, alternative splicing of G6PD5 can direct the NADPH-producing OPPP reactions to the cytosolic face of the ER, where they may operate as membrane-bound metabolon to support several important biosynthetic pathways of plant cells.

RNA biology takes root in plant systems

Advances in RNA biology such as RNAi, CRISPR, and the first mRNA vaccine represent the enormous potential of RNA research to address current problems. Additionally, plants are a diverse and undeniably essential resource for life threatened by climate change, loss of arable land, and pollution. Different aspects of RNA such as its processing, modification and structure are intertwined with plant development, physiology and stress response. This report details the findings of researchers around the world during the 23rd Penn State Symposium in Plant Biology with a focus in RNA biology.

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.

Comprehensive profiling of alternative splicing landscape during secondary dormancy in oilseed rape (Brassica napus L.)

Alternative splicing is a general mechanism that regulates gene expression at the post-transcriptional level, which increases the transcriptomic diversity. Oilseed rape (Brassica napus L.), one of the main oil crops worldwide, is prone to secondary dormancy. However, how alternative splicing landscape of oilseed rape seed changes in response to secondary dormancy is unknown. Here, we analyzed twelve RNA-seq libraries from varieties "Huaiyou-SSD-V1" and "Huaiyou-WSD-H2" which exhibited high (> 95%) and low (< 5%) secondary dormancy potential, respectively, and demonstrated that alternative splicing changes led to a significant increase with the diversity of the transcripts in response to secondary dormancy induction via polyethylene glycol 6000 (PEG6000) treatment. Among the four basic alternative splicing types, intron retention dominates, and exon skipping shows the rarest frequency. A total of 8% of expressed genes had two or more transcripts after PEG treatment. Further analysis revealed that global isoform expression percentage variations in alternative splicing in differently expressed genes (DEGs) is more than three times as much as those in non-DEGs, suggesting alternative splicing change is associated with the transcriptional activity change in response to secondary dormancy induction. Eventually, 342 differently spliced genes (DSGs) associated with secondary dormancy were identified, five of which were validated by RT-PCR. The number of the overlapped genes between DSGs and DEGs associated with secondary dormancy was much less than that of either DSGs or DEGs, suggesting that DSGs and DEGs may independently regulates secondary dormancy. Functional annotation analysis of DSGs revealed that spliceosome components are overrepresented among the DSGs, including small nuclear ribonucleoprotein particles (snRNPs), serine/arginine-rich (SR) proteins, and other splicing factors. Thus, it is proposed that the spliceosome components could be exploited to reduce secondary dormancy potential in oilseed rape.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01314-8.

Moderate Soil Drying-Induced Alternative Splicing Provides a Potential Novel Approach for the Regulation of Grain Filling in Rice Inferior Spikelets

Poor grain filling of inferior spikelets, especially in some large-panicle rice varieties, is becoming a major limitation in breaking the ceiling of rice production. In our previous studies, we proved that post-anthesis moderate soil drying (MD) was an effective way to promote starch synthesis and inferior grain filling. As one of the most important regulatory processes in response to environmental cues and at different developmental stages, the function of alternative splicing (AS) has not yet been revealed in regulating grain filling under MD conditions. In this study, AS events at the most active grain-filling stage were identified in inferior spikelets under well-watered control (CK) and MD treatments. Of 16,089 AS events, 1840 AS events involving 1392 genes occurred differentially between the CK and MD treatments, many of which function on spliceosome, ncRNA metabolic process, starch, and sucrose metabolism, and other functions. Some of the splicing factors and starch synthesis-related genes, such as SR protein, hnRNP protein, OsAGPL2, OsAPS2, OsSSIVa, OsSSIVb, OsGBSSII, and OsISA1 showed differential AS changes under MD treatment. The expression of miR439f and miR444b was reduced due to an AS event which occurred in the intron where miRNAs were located in the MD-treated inferior spikelets. On the contrary, OsAGPL2, an AGPase encoding gene, was alternatively spliced, resulting in different transcripts with or without the miR393b binding site, suggesting a potential mechanism for miRNA-mediated gene regulation on grain filling of inferior spikelets in response to MD treatment. This study provides some new insights into the function of AS on the MD-promoted grain filling of inferior spikelets, and potential application in agriculture to increase rice yields by genetic approaches.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.

Genome-Wide Identification and Characterization of CPR5 Genes in Gossypium Reveals Their Potential Role in Trichome Development

Trichomes protect plants against insects, microbes, herbivores, and abiotic damages and assist seed dispersal. The function of CPR5 genes have been found to be involved in the trichome development but the research on the underlying genetic and molecular mechanisms are extremely limited. Herein, genome wide identification and characterization of CPR5 genes was performed. In total, 26 CPR5 family members were identified in Gossypium species. Phylogenetic analysis, structural characteristics, and synteny analysis of CPR5s showed the conserved evolution relationships of CPR5. The promoter analysis of CPR5 genes revealed hormone, stress, and development-related cis-elements. Gene ontology (GO) enrichment analysis showed that the CPR5 genes were largely related to biological regulation, developmental process, multicellular organismal process. Protein-protein interaction analysis predicted several trichome development related proteins (SIM, LGO, and GRL) directly interacting with CPR5 genes. Further, nine putative Gossypium-miRNAs were also identified, targeting Gossypium CPR5 genes. RNA-Seq data of G. arboreum (with trichomes) and G. herbaceum (with no trichomes) was used to perform the co-expression network analysis. GheCPR5.1 was identified as a hub gene in a co-expression network analysis. RT-qPCR of GheCPR5.1 gene in different tissues suggests that this gene has higher expressions in the petiole and might be a key candidate involved in the trichome development. Virus induced gene silencing of GheCPR5.1 (Ghe02G17590) confirms its role in trichome development and elongation. Current results provide proofs of the possible role of CPR5 genes and provide preliminary information for further studies of GheCPR5.1 functions in trichome development.

Relevance and Regulation of Alternative Splicing in Plant Heat Stress Response: Current Understanding and Future Directions

Alternative splicing (AS) is a major mechanism for gene expression in eukaryotes, increasing proteome diversity but also regulating transcriptome abundance. High temperatures have a strong impact on the splicing profile of many genes and therefore AS is considered as an integral part of heat stress response. While many studies have established a detailed description of the diversity of the RNAome under heat stress in different plant species and stress regimes, little is known on the underlying mechanisms that control this temperature-sensitive process. AS is mainly regulated by the activity of splicing regulators. Changes in the abundance of these proteins through transcription and AS, post-translational modifications and interactions with exonic and intronic cis-elements and core elements of the spliceosomes modulate the outcome of pre-mRNA splicing. As a major part of pre-mRNAs are spliced co-transcriptionally, the chromatin environment along with the RNA polymerase II elongation play a major role in the regulation of pre-mRNA splicing under heat stress conditions. Despite its importance, our understanding on the regulation of heat stress sensitive AS in plants is scarce. In this review, we summarize the current status of knowledge on the regulation of AS in plants under heat stress conditions. We discuss possible implications of different pathways based on results from non-plant systems to provide a perspective for researchers who aim to elucidate the molecular basis of AS under high temperatures.

A rust fungus effector directly binds plant pre-mRNA splice site to reprogram alternative splicing and suppress host immunity

Alternative splicing (AS) is a crucial post-transcriptional regulatory mechanism in plant resistance. However, whether and how plant pathogens target splicing in their host remains mostly unknown. For example, although infection by Puccinia striiformis f. sp. tritici (Pst), a pathogenic fungus that severely affects the yield of wheat worldwide, has been shown to significantly influence the levels of alternatively spliced transcripts in the host, the mechanisms that govern this process, and its functional consequence have not been examined. Here, we identified Pst_A23 as a new Pst arginine-rich effector that localizes to host nuclear speckles, nuclear regions enriched in splicing factors. We demonstrated that transient expression of Pst_A23 suppresses plant basal defence dependent on the Pst_A23 nuclear speckle localization and that this protein plays an important role in virulence, stable silencing of which improves wheat stripe rust resistance. Remarkably, RNA-Seq data revealed that AS patterns of 588 wheat genes are altered in Pst_A23-overexpressing lines compared to control plants. To further examine the direct relationship between Pst_A23 and AS, we confirmed direct binding between two RNA motifs predicted from these altered splicing sites and Pst_A23 in vitro. The two RNA motifs we chose occur in the cis-element of TaXa21-H and TaWRKY53, and we validated that Pst_A23 overexpression results in decreased functional transcripts of TaXa21-H and TaWRKY53 while silencing of TaXa21-H and TaWRKY53 impairs wheat resistance to Pst. Overall, this represents formal evidence that plant pathogens produce 'splicing' effectors, which regulate host pre-mRNA splicing by direct engagement of the splicing sites, thereby interfering with host immunity.

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.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro)

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.

XAP5 CIRCADIAN TIMEKEEPER specifically modulates 3' splice site recognition and is important for circadian clock regulation partly by alternative splicing of LHY and TIC

Pre-mRNA splicing is an essential step during gene expression, which takes place in the spliceosome, a large dynamic ribonucleoprotein complex assembled in a stepwise manner. During the last decade, several spliceosomal mutants were functionally identified to cause a lengthened circadian period by introducing intron retention defects into circadian clock genes in Arabidopsis. However, the spliceosomal components that play opposite roles in the circadian period via alternative 3' splice site (Alt 3'ss) are largely unknown. Here, we demonstrated that XCT (XAP5 CIRCADIAN TIMEKEEPER) is a key spliceosomal component associated with multiple splicing factors. Moreover, genome-wide analysis revealed that inactivation of XCT particularly results in defects in Alt 3'ss recognition by RNA sequencing. Further analysis indicated that a strong alteration in the 3' splice sites of LHY and TIC partly accounts for the shortened circadian period of the xct mutant. Therefore, our results demonstrated that mutations in XCT shortened the circadian period partly by alternative splicing of LHY and TIC particularly in 3' splice site recognition, which provides new insight into the link between alternative splicing and the circadian clock.

Alternative Splicing Dynamics During the Lifecycle of Salvia miltiorrhiza Root Revealed the Fine Tuning in Root Development and Ingredients Biosynthesis

Alternative splicing (AS) is an essential post-transcriptional process that enhances the coding and regulatory potential of the genome, thereby strongly influencing multiple plant physiology processes, such as metabolic biosynthesis. To explore how AS affects the root development and synthesis of tanshinones and phenolic acid pathways in Salvia miltiorrhiza roots, we investigated the dynamic landscape of AS events in S. miltiorrhiza roots during an annual life history. Temporal profiling represented a distinct temporal variation of AS during the entire development stages, showing the most abundant AS events at the early seedling stage (ES stage) and troughs in 45 days after germination (DAG) and 120 DAG. Gene ontology (GO) analysis indicated that physiological and molecular events, such as lateral root formation, gravity response, RNA splicing regulation, and mitogen-activated protein kinase (MAPK) cascade, were greatly affected by AS at the ES stage. AS events were identified in the tanshinones and phenolic acids pathways as well, especially for the genes for the branch points of the pathways as SmRAS and SmKSL1. Fifteen Ser/Arg-rich (SR) proteins and eight phosphokinases (PKs) were identified with high transcription levels at the ES stage, showing their regulatory roles for the high frequency of AS in this stage. Simultaneously, a co-expression network that includes 521 highly expressed AS genes, SRs, and PKs, provides deeper insight into the mechanism for the variable programming of AS.

Transcriptome profiling based on Illumina- and SMRT-based RNA-seq reveals circadian regulation of key pathways in flower bud development in walnut

Flower bud development is a defining feature of walnut, which contributes to the kernel yield, yield stability, fruit quality and commodity value. However, little is known about the mechanism of the flower bud development in walnut. Here, the stages of walnut female flower bud development were divided into five period (P01-05) by using histological observation. They were further studied through PacBio Iso-Seq and RNA-seq analysis. Accordingly, we obtained 52,875 full-length transcripts, where 4,579 were new transcripts, 3,065 were novel genes, 1,437 were consensus lncRNAs and 20,813 were alternatively spliced isoforms. These transcripts greatly improved the current genome annotation and enhanced our understanding of the walnut transcriptome. Next, RNA sequencing of female flower buds at five periods revealed that circadian rhythm-plant was commonly enriched along with the flower bud developmental gradient. A total of 14 differentially expressed genes (DEGs) were identified, and six of them were confirmed by real-time quantitative analysis. Additionally, six and two differentially expressed clock genes were detected to be regulated by AS events and lncRNAs, respectively. All these detected plant circadian genes form a complex interconnected network to regulate the flower bud development. Thus, investigation of key genes associated with the circadian clock could clarify the process of flower bud development in walnut.

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.

Regulation of alternative splicing in response to temperature variation in plants

Plants have evolved numerous molecular strategies to cope with perturbations in environmental temperature, and to adjust growth and physiology to limit the negative effects of extreme temperature. One of the strategies involves alternative splicing of primary transcripts to encode alternative protein products or transcript variants destined for degradation by nonsense-mediated decay. Here, we review how changes in environmental temperature-cold, heat, and moderate alterations in temperature-affect alternative splicing in plants, including crops. We present examples of the mode of action of various temperature-induced splice variants and discuss how these alternative splicing events enable favourable plant responses to altered temperatures. Finally, we point out unanswered questions that should be addressed to fully utilize the endogenous mechanisms in plants to adjust their growth to environmental temperature. We also indicate how this knowledge might be used to enhance crop productivity in the future.