Alternative Splicing and Protein Diversity: Plants Versus Animals

Geneapp: A web application for visualizing alternative splicing for biomedicine

Alternative Splicing (AS) is an essential mechanism for eukaryotes. However, the consequences of deleting a single exon can be dramatic for the organism and can lead to cancer in humans. Additionally, alternative 5' and 3' splice sites, which define the boundaries of exons, also play key roles to human disorders. Therefore, Investigating AS events is crucial for understanding the molecular basis of human diseases and developing therapeutic strategies. Workflow for AS event analysis can be sampling followed by data analysis with bioinformatics to identify the different AS events in the control and case samples, data visualization for curation, and selection of relevant targets for experimental validation. The raw output of the analysis software does not favor the inspection of events by bioinformaticians requiring custom scripts for data visualization. In this work, we propose the Geneapp application with three modules: GeneappScript, GeneappServer, and GeneappExplorer. GeneappScript is a wrapper that assists in identifying AS in samples compared in two different approaches, while GeneappServer integrates data from AS analysis already performed by the user. In GeneappExplorer, the user visualizes the previous dataset by exploring AS events in genes with functional annotation. This targeted screens that Geneapp allows to perform helps in the identification of targets for experimental validation to confirm the hypotheses under study. The Geneapp is freely available for non-commercial use at https://geneapp.net to advance research on AS for bioinformatics.

The Many Roads from Alternative Splicing to Cancer: Molecular Mechanisms Involving Driver Genes

Cancer driver genes are either oncogenes or tumour suppressor genes that are classically activated or inactivated, respectively, by driver mutations. Alternative splicing-which produces various mature mRNAs and, eventually, protein variants from a single gene-may also result in driving neoplastic transformation because of the different and often opposed functions of the variants of driver genes. The present review analyses the different alternative splicing events that result in driving neoplastic transformation, with an emphasis on their molecular mechanisms. To do this, we collected a list of 568 gene drivers of cancer and revised the literature to select those involved in the alternative splicing of other genes as well as those in which its pre-mRNA is subject to alternative splicing, with the result, in both cases, of producing an oncogenic isoform. Thirty-one genes fall into the first category, which includes splicing factors and components of the spliceosome and splicing regulators. In the second category, namely that comprising driver genes in which alternative splicing produces the oncogenic isoform, 168 genes were found. Then, we grouped them according to the molecular mechanisms responsible for alternative splicing yielding oncogenic isoforms, namely, mutations in cis splicing-determining elements, other causes involving non-mutated cis elements, changes in splicing factors, and epigenetic and chromatin-related changes. The data given in the present review substantiate the idea that aberrant splicing may regulate the activation of proto-oncogenes or inactivation of tumour suppressor genes and details on the mechanisms involved are given for more than 40 driver genes.

Physiological function and regulation of ascorbate peroxidase isoforms

Ascorbate peroxidase (APX) reduces H2O2 to H2O by utilizing ascorbate as a specific electron donor and constitutes the ascorbate-glutathione cycle in organelles of plants including chloroplasts, cytosol, mitochondria, and peroxisomes. It has been almost 40 years since APX was discovered as an important plant-specific H2O2-scavenging enzyme, during which time many research groups have conducted molecular physiological analyses. It is now clear that APX isoforms function not only just as antioxidant enzymes but also as important factors in intracellular redox regulation through the metabolism of reactive oxygen species. The function of APX isoforms is regulated at multiple steps, from the transcriptional level to post-translational modifications of enzymes, thereby allowing them to respond flexibly to ever-changing environmental factors and physiological phenomena such as cell growth and signal transduction. In this review, we summarize the physiological functions and regulation mechanisms of expression of each APX isoform.

FER-like iron deficiency-induced transcription factor (FIT) accumulates in nuclear condensates

The functional importance of nuclear protein condensation remains often unclear. The bHLH FER-like iron deficiency-induced transcription factor (FIT) controls iron acquisition and growth in plants. Previously described C-terminal serine residues allow FIT to interact and form active transcription factor complexes with subgroup Ib bHLH factors such as bHLH039. FIT has lower nuclear mobility than mutant FITmSS271AA. Here, we show that FIT undergoes a light-inducible subnuclear partitioning into FIT nuclear bodies (NBs). Using quantitative and qualitative microscopy-based approaches, we characterized FIT NBs as condensates that were reversible and likely formed by liquid-liquid phase separation. FIT accumulated preferentially in NBs versus nucleoplasm when engaged in protein complexes with itself and with bHLH039. FITmSS271AA, instead, localized to NBs with different dynamics. FIT colocalized with splicing and light signaling NB markers. The NB-inducing light conditions were linked with active FIT and elevated FIT target gene expression in roots. FIT condensation may affect nuclear mobility and be relevant for integrating environmental and Fe nutrition signals.

Alternative splicing patterns of hnrnp genes in gill tissues of rainbow trout (Oncorhynchus mykiss) during salinity changes

Alternative splicing (AS) plays an important role in various physiological processes in eukaryotes, such as the stress response. However, patterns of AS events remain largely unexplored during salinity acclimation in fishes. In this study, we conducted AS analysis using RNA-seq datasets to explore splicing patterns in the gill tissues of rainbow trout exposed to altered salinity environments, ranging from 0 ‰ (T0) to 30 ‰ (T30). The results revealed 1441, 351, 483, 1051 and 1049 differentially alternatively spliced (DAS) events in 5 pairwise comparisons, including T6 vs. T0, T12 vs. T0, T18 vs. T0, T24 vs. T0, and T30 vs. T0, respectively. These DAS events were derived from 1290, 328, 444, 963 and 948 genes. Enrichment analysis indicated that these DAS genes were related to RNA splicing and processing. Among these, 14 DAS genes were identified as members of the large heterogeneous nuclear RNP (hnRNP) gene family. Alternative 3' splice site (A3SS), exon skipping (SE) and intron retention (RI) events resulted in the fragmentation or even loss of the functional RNA recognition motif (RRM) domains in hnrnpa0, hnrnp1a, hnrnp1b and hnrnpc genes. The incomplete RRM domains would hinder the interactions between hnRNP genes and pre-mRNAs. It would in turn influence the splicing patterns and mRNA stability of downstream target genes in response to salinity changes. The study provides insights into salinity acclimation in gill tissues of rainbow trout and serves as a significant reference on the osmoregulation mechanisms at post-transcription regulation levels in fish.

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).

Impact of alternative splicing on Arabidopsis proteome

Alternative splicing is an important regulatory process in eukaryotes. In plants, the major form of alternative splicing is intron retention. Despite its importance, the global impact of AS on the Arabidopsis proteome has not been investigated. In this study, we address this gap by performing a comprehensive integrated analysis of how changes in AS can affect the Arabidopsis proteome using mutants that disrupt ACINUS and PININ, two evolutionarily conserved alternative splicing factors. We used tandem mass tagging (TMT) with real-time search MS3 (RTS-SPS-MS3) coupled with extensive sample fractionations to achieve very high coverage and accurate protein quantification. We then integrated our proteomic data with transcriptomic data to assess how transcript changes and increased intron retention (IIR) affect the proteome. For differentially expressed transcripts, we have observed a weak to moderate correlation between transcript changes and protein changes. Our studies revealed that some IIRs have no effect on either transcript or protein levels, while some IIRs can significantly affect protein levels. Surprisingly, we found that IIRs have a much smaller effect on increasing protein diversity. Notably, the increased intron retention events detected in the double mutant are also detected in the WT under various biotic or abiotic stresses. We further investigated the characteristics of the retained introns. Our extensive proteomic data help to guide the phenotypic analysis and reveal that collective protein changes contribute to the observed phenotypes of the increased anthocyanin, pale green, reduced growth, and short root observed in the acinus pnn double mutant. Overall, our study provides insight into the intricate regulatory mechanism of intron retention and its impact on protein abundance in plants.

Exploring the anticancer properties of the gum of Ferula gummosa: impact on cytotoxicity, caspase 3/7 activity and apoptosis, and gene expression in SW-480 cells

Ferula gummosa Boiss. is a well-known Iranian endemic plant that grows in the north and northeast regions of Iran. In Iranian traditional medicine, its gum is utilized to treat inflammation, pain, and infections of the gastrointestinal system. However, no studies have been conducted to investigate the anticancer potential of its gum against colorectal cancer cells. This study aimed to identify the chemical components of the gum of F. gummosa and investigate its effects on SW-480 cells. The experiments included MTT, clonogenic, micronucleus formation, acridine orange/ethidium bromide stain, DNA degradation, caspase 3/7 activity assay, and in vitro wound-healing experiment and investigating the expression of BAX, BCL2, MTOR, and PTEN genes. Chemical analysis using GC/MS identified 102 compounds. The gum had a significant cytotoxic effect on SW-480 cells, with an IC50 value of 1.8 µg/ml for 48 hours. The gum induced apoptosis. Microscopic observations revealed a decrease in cell proliferation, as evidenced by nuclear condensation, increased micronucleus formation, and inhibition of colony formation. Additionally, the gum suppressed cell migration, induced the expression of PTEN and BAX, and down-regulated MTOR and BCL2 genes. These findings suggest that Ferula gummosa has strong cytotoxic properties and warrants further investigation.

iFLAS: positive-unlabeled learning facilitates full-length transcriptome-based identification and functional exploration of alternatively spliced isoforms in maize

The advent of full-length transcriptome sequencing technologies has accelerated the discovery of novel splicing isoforms. However, existing alternative splicing (AS) tools are either tailored for short-read RNA-Seq data or designed for human and animal studies. The disparities in AS patterns between plants and animals still pose a challenge to the reliable identification and functional exploration of novel isoforms in plants. Here, we developed integrated full-length alternative splicing analysis (iFLAS), a plant-optimized AS toolkit that introduced a semi-supervised machine learning method known as positive-unlabeled (PU) learning to accurately identify novel isoforms. iFLAS also enables the investigation of AS functions from various perspectives, such as differential AS, poly(A) tail length, and allele-specific AS (ASAS) analyses. By applying iFLAS to three full-length transcriptome sequencing datasets, we systematically identified and functionally characterized maize (Zea mays) AS patterns. We found intron retention not only introduces premature termination codons, resulting in lower expression levels of isoforms, but may also regulate the length of 3'UTR and poly(A) tail, thereby affecting the functional differentiation of isoforms. Moreover, we observed distinct ASAS patterns in two genes within heterosis offspring, highlighting their potential value in breeding. These results underscore the broad applicability of iFLAS in plant full-length transcriptome-based AS research.

Unlocking allelic variation in circadian clock genes to develop environmentally robust and productive crops

MAIN CONCLUSION

Molecular mechanisms of biological rhythms provide opportunities to harness functional allelic diversity in core (and trait- or stress-responsive) oscillator networks to develop more climate-resilient and productive germplasm. The circadian clock senses light and temperature in day-night cycles to drive biological rhythms. The clock integrates endogenous signals and exogenous stimuli to coordinate diverse physiological processes. Advances in high-throughput non-invasive assays, use of forward- and inverse-genetic approaches, and powerful algorithms are allowing quantitation of variation and detection of genes associated with circadian dynamics. Circadian rhythms and phytohormone pathways in response to endogenous and exogenous cues have been well documented the model plant Arabidopsis. Novel allelic variation associated with circadian rhythms facilitates adaptation and range expansion, and may provide additional opportunity to tailor climate-resilient crops. The circadian phase and period can determine adaptation to environments, while the robustness in the circadian amplitude can enhance resilience to environmental changes. Circadian rhythms in plants are tightly controlled by multiple and interlocked transcriptional-translational feedback loops involving morning (CCA1, LHY), mid-day (PRR9, PRR7, PRR5), and evening (TOC1, ELF3, ELF4, LUX) genes that maintain the plant circadian clock ticking. Significant progress has been made to unravel the functions of circadian rhythms and clock genes that regulate traits, via interaction with phytohormones and trait-responsive genes, in diverse crops. Altered circadian rhythms and clock genes may contribute to hybrid vigor as shown in Arabidopsis, maize, and rice. Modifying circadian rhythms via transgenesis or genome-editing may provide additional opportunities to develop crops with better buffering capacity to environmental stresses. Models that involve clock gene‒phytohormone‒trait interactions can provide novel insights to orchestrate circadian rhythms and modulate clock genes to facilitate breeding of all season crops.

Lectin diversity and their positive roles in WSSV replication through regulation of calreticulin expression and inhibiting ALFs expression

In biological evolution, gene duplication (GD) generates new genes to facilitate new functions. C-type lectins (CTLs) in crayfish have been extended by GD to expand their family members. In this study, four CTL genes generated by GD were identified from Procambarus clarkii (PcLec1-4). Among these four genes, PcLec1 can also generate new isoforms with different numbers of tandem repeats through DNA slip mispairing. PcLec1-4 was widely expressed in multiple tissues. The expression levels of PcLec1-4 were upregulated in the intestine of P. clarkii upon white spot syndrome virus (WSSV) challenge at multiple time points. Further analysis indicated that GATA transcription factor regulated PcLec1-4 expression. RNA interference and recombinant PcLec1-4 protein injection experiments suggested that PcLec1-4 promoted the expression of calreticulin (PcCRT) and negatively regulated the expression of antimicrobial peptides, thereby promoting WSSV replication. This study contributes to the understanding of the function of CTLs produced by GD during WSSV invasion in crustaceans.

Bulk and single-cell alternative splicing analyses reveal roles of TRA2B in myogenic differentiation

Alternative splicing (AS) disruption has been linked to disorders of muscle development, as well as muscular atrophy. However, the precise changes in AS patterns that occur during myogenesis are not well understood. Here, we employed isoform long-reads RNA-seq (Iso-seq) and single-cell RNA-seq (scRNA-seq) to investigate the AS landscape during myogenesis. Our Iso-seq data identified 61,146 full-length isoforms representing 11,682 expressed genes, of which over 52% were novel. We identified 38,022 AS events, with most of these events altering coding sequences and exhibiting stage-specific splicing patterns. We identified AS dynamics in different types of muscle cells through scRNA-seq analysis, revealing genes essential for the contractile muscle system and cytoskeleton that undergo differential splicing across cell types. Gene-splicing analysis demonstrated that AS acts as a regulator, independent of changes in overall gene expression. Two isoforms of splicing factor TRA2B play distinct roles in myogenic differentiation by triggering AS of TGFBR2 to regulate canonical TGF-β signalling cascades differently. Our study provides a valuable transcriptome resource for myogenesis and reveals the complexity of AS and its regulation during myogenesis.

Nonsense-mediated mRNA decay pathway in plants under stress: general gene regulatory mechanism and advances

MAIN CONCLUSION

Nonsense-mediated mRNA decay in eukaryotes is vital to cellular homeostasis. Further knowledge of its putative role in plant RNA metabolism under stress is pivotal to developing fitness-optimizing strategies. Nonsense-mediated mRNA decay (NMD), part of the mRNA surveillance pathway, is an evolutionarily conserved form of gene regulation in all living organisms. Degradation of mRNA-bearing premature termination codons and regulation of physiological RNA levels highlight NMD's role in shaping the cellular transcriptome. Initially regarded as purely a tool for cellular RNA quality control, NMD is now considered to mediate various aspects of plant developmental processes and responses to environmental changes. Here we offer a basic understanding of NMD in eukaryotes by explaining the concept of premature termination codon recognition and NMD complex formation. We also provide a detailed overview of the NMD mechanism and its role in gene regulation. The potential role of effectors, including ABCE1, in ribosome recycling during the translation process is also explained. Recent reports of alternatively spliced variants of corresponding genes targeted by NMD in Arabidopsis thaliana are provided in tabular format. Detailed figures are also provided to clarify the NMD concept in plants. In particular, accumulating evidence shows that NMD can serve as a novel alternative strategy for genetic manipulation and can help design RNA-based therapies to combat stress in plants. A key point of emphasis is its function as a gene regulatory mechanism as well as its dynamic regulation by environmental and developmental factors. Overall, a detailed molecular understanding of the NMD mechanism can lead to further diverse applications, such as improving cellular homeostasis in living organisms.

Diurnal regulation of alternative splicing associated with thermotolerance in rice by two glycine-rich RNA-binding proteins

Rice (Oryza sativa L.) production is threatened by global warming associated with extreme high temperatures, and rice heat sensitivity is differed when stress occurs between daytime and nighttime. However, the underlying molecular mechanism are largely unknown. We show here that two glycine-rich RNA binding proteins, OsGRP3 and OsGRP162, are required for thermotolerance in rice, especially at nighttime. The rhythmic expression of OsGRP3/OsGRP162 peaks at midnight, and at these coincident times, is increased by heat stress. This is largely dependent on the evening complex component OsELF3-2. We next found that the double mutant of OsGRP3/OsGRP162 is strikingly more sensitive to heat stress in terms of survival rate and seed setting rate when comparing to the wild-type plants. Interestingly, the defect in thermotolerance is more evident when heat stress occurred in nighttime than that in daytime. Upon heat stress, the double mutant of OsGRP3/OsGRP162 displays globally reduced expression of heat-stress responsive genes, and increases of mRNA alternative splicing dominated by exon-skipping. This study thus reveals the important role of OsGRP3/OsGRP162 in thermotolerance in rice, and unravels the mechanism on how OsGRP3/OsGRP162 regulate thermotolerance in a diurnal manner.

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.

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

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

The many faces of lysine acylation in proteins: Phytohormones as unexplored substrates

Protein post-translational modification (PTM) is a ubiquitous process that occurs in most proteins. Lysine residues containing an ε-amino group are recognized as hotspots for the addition of different chemical groups. Lysine acetylation, extensively studied in histones, serves as an epigenetic hallmark capable of promoting changes in chromatin structure and availability. Acyl groups derived from molecules involved in carbohydrate and lipid metabolisms, such as lactate, succinate and hydroxybutyrate, were identified as lysine modifications of histones and other proteins. Lysine-acyltransferases do not exhibit significant substrate specificity concerning acyl donors. Furthermore, plant hormones harboring acyl groups often form conjugates with free amino acids to regulate their activity and function during plant physiological processes and responses, a process mediated by GH3 enzymes. Besides forming low-molecular weight conjugates, auxins have been shown to covalently modify proteins in bean seeds. Aside from auxins, other phytohormones with acyl groups are unexplored potential substrates for post-translational acylation of proteins. Using MS data searches, we revealed various proteins with lysine residues linked to auxin, abscisic acid, gibberellic acid, jasmonic acid, and salicylic acid. These findings raise compelling questions about the ability of plant hormones harboring carboxyl groups to serve as new candidates for protein acylation and acting in protein PTM and modulation.

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.

Phosphorylation mediated regulation of RNA splicing in plants

For the past two decades, the study of alternative splicing (AS) and its involvement in plant development and stress response has grown in popularity. Only recently however, has the focus shifted to the study of how AS regulation (or lack-thereof) affects downstream mRNA and protein landscapes and how these AS regulatory events impact plant development and stress tolerance. In humans, protein phosphorylation represents one of the predominant mechanisms by which AS is regulated and thus the protein kinases governing these phosphorylation events are of interest for further study. Large-scale phosphoproteomic studies in plants have consistently found that RNA splicing-related proteins are extensively phosphorylated, however, the signaling pathways involved in AS regulation have not been resolved. In this mini-review, we summarize our current knowledge of the three major splicing-related protein kinase families in plants that are suggested to mediate AS phospho-regulation and draw comparisons to their metazoan orthologs. We also summarize and contextualize the phosphorylation events identified as occurring on splicing-related protein families to illustrate the high degree to which splicing-related proteins are modified, placing a new focus on elucidating the impacts of AS at the protein and PTM-level.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

Protein diversification through post-translational modifications, alternative splicing, and gene duplication

Proteins provide the basis for cellular function. Having multiple versions of the same protein within a single organism provides a way of regulating its activity or developing novel functions. Post-translational modifications of proteins, by means of adding/removing chemical groups to amino acids, allow for a well-regulated and controlled way of generating functionally distinct protein species. Alternative splicing is another method with which organisms possibly generate new isoforms. Additionally, gene duplication events throughout evolution generate multiple paralogs of the same genes, resulting in multiple versions of the same protein within an organism. In this review, we discuss recent advancements in the study of these three methods of protein diversification and provide illustrative examples of how they affect protein structure and function.

Differential RNA-Seq Analysis Predicts Genes Related to Terpene Tailoring in Caryopteris × clandonensis

Enzymatic terpene functionalization is an essential part of plant secondary metabolite diversity. Within this, multiple terpene-modifying enzymes are required to enable the chemical diversity of volatile compounds essential in plant communication and defense. This work sheds light on the differentially transcribed genes within Caryopteris × clandonensis that are capable of functionalizing cyclic terpene scaffolds, which are the product of terpene cyclase action. The available genomic reference was subjected to further improvements to provide a comprehensive basis, where the number of contigs was minimized. RNA-Seq data of six cultivars, Dark Knight, Grand Bleu, Good as Gold, Hint of Gold, Pink Perfection, and Sunny Blue, were mapped on the reference, and their distinct transcription profile investigated. Within this data resource, we detected interesting variations and additionally genes with high and low transcript abundancies in leaves of Caryopteris × clandonensis related to terpene functionalization. As previously described, different cultivars vary in their modification of monoterpenes, especially limonene, resulting in different limonene-derived molecules. This study focuses on predicting the cytochrome p450 enzymes underlying this varied transcription pattern between investigated samples. Thus, making them a reasonable explanation for terpenoid differences between these plants. Furthermore, these data provide the basis for functional assays and the verification of putative enzyme activities.

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.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

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.

Cell-type-specific alternative splicing in the Arabidopsis germline

During sexual reproduction in flowering plants, the two haploid sperm cells (SCs) embedded within the cytoplasm of a growing pollen tube are carried to the embryo sac for double fertilization. Pollen development in flowering plants is a dynamic process that encompasses changes at transcriptome and epigenome levels. While the transcriptome of pollen and SCs in Arabidopsis (Arabidopsis thaliana) is well documented, previous analyses have mostly been based on gene-level expression. In-depth transcriptome analysis, particularly the extent of alternative splicing (AS) at the resolution of SC and vegetative nucleus (VN), is still lacking. Therefore, we performed RNA-seq analysis to generate a spliceome map of Arabidopsis SCs and VN isolated from mature pollen grains. Based on our de novo transcriptome assembly, we identified 58,039 transcripts, including 9,681 novel transcripts, of which 2,091 were expressed in SCs and 3,600 in VN. Four hundred and sixty-eight genes were regulated both at gene and splicing levels, with many having functions in mRNA splicing, chromatin modification, and protein localization. Moreover, a comparison with egg cell RNA-seq data uncovered sex-specific regulation of transcription and splicing factors. Our study provides insights into a gamete-specific AS landscape at unprecedented resolution.

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/.

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

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

Calcium signalling in weeds under herbicide stress: An outlook

The continuous use of herbicides for controlling weeds has led to the evolution of resistance to all major herbicidal modes of action globally. Every year, new cases of herbicide resistance are reported. Resistance is still in progress in many species, which must be stopped before it becomes a worldwide concern. Several herbicides are known to cause stressful conditions that resemble plant abiotic stresses. Variation in intracellular calcium (Ca²⁺) concentration is a primary event in a wide range of biological processes in plants, including adaptation to various biotic and abiotic stresses. Ca²⁺ acts as a secondary messenger, connecting various environmental stimuli to different biological processes, especially during stress rejoindering in plants. Even though many studies involving Ca²⁺ signalling in plants have been published, there have been no studies on the roles of Ca²⁺ signalling in herbicide stress response. Hence, this mini-review will highlight the possible sensing and molecular communication via Ca²⁺ signals in weeds under herbicide stress. It will also discuss some critical points regarding integrating the sensing mechanisms of multiple stress conditions and subsequent molecular communication. These signalling responses must be addressed in the future, enabling researchers to discover new herbicidal targets.

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.

Plant Intron-Splicing Efficiency Database (PISE): exploring splicing of ∼1,650,000 introns in Arabidopsis, maize, rice, and soybean from ∼57,000 public RNA-seq libraries

Intron retention is the most common alternative splicing event in plants and plays a crucial role in the responses of plants to environmental signals. Despite a large number of RNA-seq libraries from different treatments and genetic mutants stored in public domains, a resource for querying the intron-splicing ratio of individual intron is still required. Here, we established the first-ever large-scale splicing efficiency database in any organism. Our database includes over 57,000 plant public RNA-seq libraries, comprising 25,283 from Arabidopsis, 17,789 from maize, 10,710 from rice, and 3,974 from soybean, and covers a total of 1.6 million introns in these four species. In addition, we manually curated and annotated all the mutant- and treatment-related libraries as well as their matched controls included in our library collection, and added graphics to display intron-splicing efficiency across various tissues, developmental stages, and stress-related conditions. The result is a large collection of 3,313 treatment conditions and 3,594 genetic mutants for discovering differentially regulated splicing efficiency. Our online database can be accessed at https://plantintron.com/ .

Age-related changes in gonadotropin-releasing hormone (GnRH) splice variants in mouse brain

Gonadotropin-releasing hormone (GnRH) is the primary regulator of the mammalian reproductive axis. We investigated the spatiotemporal expression of GnRH splice variants (V1, V2, and V3) and splicing factors (Srsf7, Srsf9, and Tra-2) in the male mice brain. Further, using in silico tools, we predicted protein structure and the reason for the low translational efficiency of V2 and V3. Messenger RNA levels of GnRH variants and splicing factors were quantified using real-time reverse transcription-polymerase chain reaction at different age groups. Our data show that expression of almost all the variants alters with aging in all the brain regions studied; even in comparison to the hypothalamus, several brain areas were found to have higher expression of these variants. Hypothalamic expression of splicing factors such as Srsf7, Srsf9, and Tra-2 also change with aging. Computational studies have translation repressors site on the V3, which probably reduces its translation efficiency. Also, V2 is an intrinsically disordered protein that might have a regulatory or signaling function. In conclusion, this study provides novel crucial information and multiple starting points for future analysis of GnRH splice variants in the brain.

Long non-coding RNAs reveal new regulatory mechanisms controlling gene expression

A plethora of non-coding RNAs have been found in eukaryotes, notably with the advent of modern sequencing technologies to analyze the transcriptome. Apart from the well-known housekeeping RNA genes (such as the ribosomal RNA or the transfer RNA), many thousands of transcripts detected are not evidently linked to a protein-coding gene. These, so called non-coding RNAs, may code for crucial regulators of gene expression, the small si/miRNAs, for small peptides (translated under specific conditions) or may act as long RNA molecules (antisense, intronic or intergenic long non-coding RNAs or lncRNAs). The lncRNAs interact with members of multiple machineries involved in gene regulation. In this review, we discussed about how plant lncRNAs permitted to discover new regulatory mechanisms acting in epigenetic control, chromatin 3D structure and alternative splicing. These novel regulations diversified the expression patterns and protein variants of target protein-coding genes and are an important element of the response of plants to environmental stresses and their adaptation to changing conditions.

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.

Regulation of alternative splicing by retrograde and light signals converges to control chloroplast proteins

Retrograde signals sent by chloroplasts control transcription in the nucleus. These signals antagonistically converge with light signals to coordinate the expression of genes involved in chloroplast functioning and seedling development. Although significant advances have been made in understanding the molecular interplay between light and retrograde signals at the transcriptional level, little is known about their interconnection at the post-transcriptional level. By using different publicly available datasets, this study addresses the influence of retrograde signaling on alternative splicing and defines the molecular and biological functions of this regulation. These analyses revealed that alternative splicing mimics transcriptional responses triggered by retrograde signals at different levels. First, both molecular processes similarly depend on the chloroplast-localized pentatricopeptide-repeat protein GUN1 to modulate the nuclear transcriptome. Secondly, as described for transcriptional regulation, alternative splicing coupled with the nonsense-mediated decay pathway effectively downregulates expression of chloroplast proteins in response to retrograde signals. Finally, light signals were found to antagonistically control retrograde signaling-regulated splicing isoforms, which consequently generates opposite splicing outcomes that likely contribute to the opposite roles these signals play in controlling chloroplast functioning and seedling development.

Spliceostatin C, a component of a microbial bioherbicide, is a potent phytotoxin that inhibits the spliceosome

Spliceostatin C (SPC) is a component of a bioherbicide isolated from the soil bacterium Burkholderia rinojensis. The chemical structure of SPC closely resembles spliceostatin A (SPA) which was characterized as an anticancer agent and splicing inhibitor. SPC inhibited the growth of Arabidopsis thaliana seedlings with an IC50 value of 2.2 µM. The seedlings exposed to SPC displayed a significant response with decreased root length and number and inhibition of gravitropism. Reverse transcriptase semi-quantitative PCR (RT-sqPCR) analyses of 19 selected genes demonstrated the active impact of SPC on the quality and quantity of transcripts that underwent intron rearrangements as well as up or down expression upon exposure to SPC. Qualitative and quantitative proteomic profiles identified 66 proteins that were significantly affected by SPC treatment. Further proteomics data analysis revealed that spliceostatin C induces hormone-related responses in Arabidopsis seedlings. In silico binding studies showed that SPC binds to a pocket between the SF3B3 and PF5A of the spliceosome.

Plant Single-Cell/Nucleus RNA-seq Workflow

Single-cell transcriptomics technologies allow researchers to investigate how individual cells, in complex multicellular organisms, differentially use their common genomic DNA. In plant biology, these technologies were recently applied to reveal the transcriptomes of various plant cells isolated from different organs and different species and in response to environmental stresses. These first studies support the potential of single-cell transcriptomics technology to decipher the biological function of plant cells, their developmental programs, cell-type-specific gene networks, programs controlling plant cell response to environmental stresses, etc. In this chapter, we provide information regarding the critical steps and important information to consider when developing an experimental design in plant single-cell biology. We also describe the current status of bioinformatics tools used to analyze single-cell RNA-seq datasets and how additional emerging technologies such as spatial transcriptomics and long-read sequencing technologies will provide additional information on the differential use of the genome by plant cells.

Transcriptomic and splicing changes underlying tomato responses to combined water and nutrient stress

Tomato is a horticultural crop of high economic and nutritional value. Suboptimal environmental conditions, such as limited water and nutrient availability, cause severe yield reductions. Thus, selection of genotypes requiring lower inputs is a goal for the tomato breeding sector. We screened 10 tomato varieties exposed to water deficit, low nitrate or a combination of both. Biometric, physiological and molecular analyses revealed different stress responses among genotypes, identifying T270 as severely affected, and T250 as tolerant to the stresses applied. Investigation of transcriptome changes caused by combined stress in roots and leaves of these two genotypes yielded a low number of differentially expressed genes (DEGs) in T250 compared to T270, suggesting that T250 tailors changes in gene expression to efficiently respond to combined stress. By contrast, the susceptible tomato activated approximately one thousand and two thousand genes in leaves and roots respectively, indicating a more generalized stress response in this genotype. In particular, developmental and stress-related genes were differentially expressed, such as hormone responsive factors and transcription factors. Analysis of differential alternative splicing (DAS) events showed that combined stress greatly affects the splicing landscape in both genotypes, highlighting the important role of AS in stress response mechanisms. In particular, several stress and growth-related genes as well as transcription and splicing factors were differentially spliced in both tissues. Taken together, these results reveal important insights into the transcriptional and post-transcriptional mechanisms regulating tomato adaptation to growth under reduced water and nitrogen inputs.

SWAP1-SFPS-RRC1 splicing factor complex modulates pre-mRNA splicing to promote photomorphogenesis in Arabidopsis

Regulation of transcription and pre–messenger RNA (mRNA) splicing is essential for transcript diversity and light-regulated gene expression in plants. Although several transcription factors have been described, only a few splicing factors are known to regulate light signaling pathways. Here, we describe the functional characterization of splicing factor SWAP1, which interacts with two previously characterized splicing factors, SFPS and RRC1, forming a ternary complex. SWAP1 also interacts with photoactivated phytochrome B (phyB). All three single and double mutant combinations display similar hypocotyl lengths under red light. SWAP1 modulates alternative splicing of a large number of genes, and a subset of these genes is coordinately regulated by SFPS, RRC1, and SWAP1. These results highlight the importance of phyB-interacting splicing factors in light-regulated plant development.

Light signals perceived by a group of photoreceptors have profound effects on the physiology, growth, and development of plants. The red/far-red light–absorbing phytochromes (phys) modulate these aspects by intricately regulating gene expression at multiple levels. Here, we report the identification and functional characterization of an RNA-binding splicing factor, SWAP1 (SUPPRESSOR-OF-WHITE-APRICOT/SURP RNA-BINDING DOMAIN-CONTAINING PROTEIN1). Loss-of-function swap1-1 mutant is hyposensitive to red light and exhibits a day length–independent early flowering phenotype. SWAP1 physically interacts with two other splicing factors, (SFPS) SPLICING FACTOR FOR PHYTOCHROME SIGNALING and (RRC1) REDUCED RED LIGHT RESPONSES IN CRY1CRY2 BACKGROUND 1 in a light-independent manner and forms a ternary complex. In addition, SWAP1 physically interacts with photoactivated phyB and colocalizes with nuclear phyB photobodies. Phenotypic analyses show that the swap1sfps, swap1rrc1, and sfpsrrc1 double mutants display hypocotyl lengths similar to that of the respective single mutants under red light, suggesting that they function in the same genetic pathway. The swap1sfps double and swap1sfpsrrc1 triple mutants display pleiotropic phenotypes, including sterility at the adult stage. Deep RNA sequencing (RNA-seq) analyses show that SWAP1 regulates the gene expression and pre–messenger RNA (mRNA) alternative splicing of a large number of genes, including those involved in plant responses to light signaling. A comparative analysis of alternative splicing among single, double, and triple mutants showed that all three splicing factors coordinately regulate the alternative splicing of a subset of genes. Our study uncovered the function of a splicing factor that modulates light-regulated alternative splicing by interacting with photoactivated phyB and other splicing factors.

Comparative Analysis of Environment-Responsive Alternative Splicing in the Inflorescences of Cultivated and Wild Tomato Species

Cultivated tomato (Solanum lycopersicum) is bred for fruit production in optimized environments, in contrast to harsh environments where their ancestral relatives thrive. The process of domestication and breeding has profound impacts on the phenotypic plasticity of plant development and the stress response. Notably, the alternative splicing (AS) of precursor message RNA (pre-mRNA), which is one of the major factors contributing to transcriptome complexity, is responsive to developmental cues and environmental change. To determine a possible association between AS events and phenotypic plasticity, we investigated environment-responsive AS events in the inflorescences of cultivated tomato and its ancestral relatives S. pimpinellifolium. Despite that similar AS frequencies were detected in the cultivated tomato variety Moneymaker and two S. pimpinellifolium accessions under the same growth conditions, 528 genes including splicing factors showed differential splicing in the inflorescences of plants grown in open fields and plastic greenhouses in the Moneymaker variety. In contrast, the two S. pimpinellifolium accessions, LA1589 and LA1781, had 298 and 268 genes showing differential splicing, respectively. Moreover, seven heat responsive genes showed opposite expression patterns in response to changing growth conditions between Moneymaker and its ancestral relatives. Accordingly, there were eight differentially expressed splice variants from genes involved in heat response in Moneymaker. Our results reveal distinctive features of AS events in the inflorescences between cultivated tomato and its ancestral relatives, and show that AS regulation in response to environmental changes is genotype dependent.

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.

Phytophthora effector PSR1 hijacks the host pre-mRNA splicing machinery to modulate small RNA biogenesis and plant immunity

Phytophthora effector PSR1 suppresses small RNA (sRNA)-mediated immunity in plants, but the underlying mechanism remains unknown. Here, we show that Phytophthora suppressor of RNA silencing 1 (PSR1) contributes to the pathogenicity of Phytophthora sojae and specifically binds to three conserved C-terminal domains of the eukaryotic PSR1-Interacting Protein 1 (PINP1). PINP1 encodes PRP16, a core pre-mRNA splicing factor that unwinds RNA duplexes and binds to primary microRNA transcripts and general RNAs. Intriguingly, PSR1 decreased both RNA helicase and RNA-binding activity of PINP1, thereby dampening sRNA biogenesis and RNA metabolism. The PSR1-PINP1 interaction caused global changes in alternative splicing (AS). A total of 5,135 genes simultaneously exhibited mis-splicing in both PSR1-overexpressing and PINP1-silenced plants. AS upregulated many mRNA transcripts that had their introns retained. The high occurrence of intron retention in AS-induced transcripts significantly promoted Phytophthora pathogen infection in Nicotiana benthamiana, and this might be caused by the production of truncated proteins. Taken together, our findings reveal a key role for PINP1 in regulating sRNA biogenesis and plant immunity.

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

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

Genome-Wide Survey and Expression Analyses of Hexokinase Family in Poplar (Populus trichocarpa)

Hexokinase (HXK) family proteins exert critical roles in catalyzing hexose phosphorylation, sugar sensing, and modulation of plant growth and stress adaptation. Nevertheless, a large amount remains unknown about the molecular profile of HXK enzymes in Populus trichocarpa, a woody model tree species. A genome-wide survey of HXK-encoding genes, including phylogenies, genomic structures, exon/intron organization, chromosomal distribution, and conserved features, was conducted, identifying six putative HXK isogenes (PtHXK1-6) in the Populus genome. The evolutionary tree demonstrated that 135 homologous HXKs between 17 plant species were categorized into four major subfamilies (type A, B, C, and D), clustering one plastidic (PtHXK3) and five mitochondrial PtHXKs grouped into type A and B, respectively. The in silico deduction prompted the presence of the conserved sugar-binding core (motif 4), phosphorylation sites (motif 2 and 3), and adenosine-binding domains (motif 7). The transcriptomic sequencing (RNA-seq) and the quantitative real-time PCR (qRT-PCR) assays revealed that three isogenes (PtHXK2, 3, and 6) were abundantly expressed in leaves, stems, and roots, while others appeared to be dominantly expressed in the reproductive tissues. Under the stress exposure, PtHXK2 and 6 displayed a significant induction upon the pathogenic fungi (Fusarium solani) infection and marked promotions by glucose feeding in roots. In contrast, the PtHXK3 and 6 are ABA-responsive genes, following a dose-dependent manner. The comprehensive analyses of the genomic patterns and expression profiling provide theoretical clues and lay a foundation for unraveling the physiological and signaling roles underlying the fine-tuned PtHXKs responding to diverse stressors.

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.

An Arabidopsis pre-RNA processing8a (prp8a) missense allele restores splicing of a subset of mis-spliced mRNAs

Eukaryotic precursor mRNAs often harbor noncoding introns that must be removed prior to translation. Accurate splicing of precursor messenger RNA depends on placement and assembly of small nuclear ribonucleoprotein (snRNP) sub-complexes of the spliceosome. Yeast (Saccharomyces cerevisiae) studies established a role in splice-site selection for PRE-RNA PROCESSING8 (PRP8), a conserved spliceosome scaffolding protein of the U5 snRNP. However, analogous splice-site selection studies in multicellular eukaryotes are lacking. Such studies are crucial for a comprehensive understanding of alternative splicing, which is extensive in plants and animals but limited in yeast. In this work, we describe an Arabidopsis (Arabidopsis thaliana) prp8a mutant that modulates splice-site selection. We isolated prp8a-14 from a screen for suppressors of pex14-6, which carries a splice-site mutation in the PEROXIN14 (PEX14) peroxisome biogenesis gene. To elucidate Arabidopsis PRP8A function in spliceosome fidelity, we combined prp8a-14 with various pex14 splice-site mutations and monitored the double mutants for physiological and molecular consequences of dysfunctional and functional peroxisomes that correspond to impaired and recovered splicing, respectively. prp8a-14 restored splicing and PEX14 function to alleles with mutations in the exonic guanine of the 5'-splice site but did not restore splicing or function to alleles with mutations in the intronic guanine of 5'- or 3'-splice sites. We used RNA-seq to reveal the systemic impact of prp8a-14 and found hundreds of differentially spliced transcripts and thousands of transcripts with significantly altered levels. Among differentially spliced transcripts, prp8a-14 significantly altered 5'- and 3'-splice-site utilization to favor sites resulting in shorter introns. This study provides a genetic platform for probing splicing in plants and hints at a role for plant PRP8 in splice-site selection.

The splicing factor 1-FLOWERING LOCUS M module spatially regulates temperature-dependent flowering by modulating FLOWERING LOCUS T and LEAFY expression

The AtSF1-FLM module spatially controls temperature-dependent flowering by negatively regulating the expression of FT and LFY in the leaf and shoot apex, respectively. Alternative splicing mediated by various splicing factors is important for the regulation of plant growth and development. Our recent reports have shown that a temperature-dependent interaction between Arabidopsis thaliana splicing factor 1 (AtSF1) and FLOWERING LOCUS M (FLM) pre-mRNA introns controls the differential production of FLM-β transcripts at different temperatures, eventually resulting in temperature-responsive flowering. However, the molecular and genetic interactions between the AtSF1-FLM module and floral activator genes remain unknown. Here, we aimed to identify the interactions among AtSF1, FLM, FLOWERING LOCUS T (FT), and LEAFY (LFY) by performing molecular and genetic analyses. FT and TWIN SISTER OF FT (TSF) expression in atsf1-2 mutants significantly increased in the morning and middle of the night at 16 and 23 °C, respectively, under long-day conditions. In addition, ft mutation suppressed the early flowering of atsf1-2 and atsf1-2 flm-3 mutants and masked the temperature response of atsf1-2 flm-3 mutants, suggesting that FT is a downstream target gene of the AtSF1-FLM module. LFY expression significantly increased in the diurnal samples of atsf1-2 mutants and in the shoot apex regions of atsf1-2 ft-10 mutants at different temperatures. The chromatin immunoprecipitation (ChIP) assay revealed that FLM directly binds to the genomic regions of LFY but not of APETALA1 (AP1). Moreover, lfy mutation suppressed the early flowering of flm-3 mutants, suggesting that LFY is another target of the AtSF1-FLM module. Our results reveal that the AtSF1-FLM module spatially modulates temperature-dependent flowering by regulating FT and LFY expressions.

Landscape and regulation of alternative splicing and alternative polyadenylation in a plant pathogenic fungus

Alternative splicing (AS) and alternative polyadenylation (APA) contribute significantly to the regulation of gene expression in higher eukaryotes. Their biological impact in filamentous fungi, however, is largely unknown. Here we combine PacBio Isoform-Sequencing and strand-specific RNA-sequencing of multiple tissues and mutant characterization to reveal the landscape and regulation of AS and APA in Fusarium graminearum. We generated a transcript annotation comprising 51 617 isoforms from 17 189 genes. In total, 4997 and 11 133 genes are alternatively spliced and polyadenylated, respectively. Majority of the AS events alter coding sequences. Unexpectedly, the AS transcripts containing premature-termination codons are not sensitive to nonsense-mediated messenger RNA decay. Unlike in yeasts and animals, distal APA sites have strong signals, but proximal APA isoforms are highly expressed in F. graminearum. The 3'-end processing factors FgRNA15, FgHRP1, and FgFIP1 play roles in promoting proximal APA site usage and intron splicing. A genome-wide increase in intron inclusion and distal APA site usage and downregulation of the spliceosomal and 3'-end processing factors were observed in older and quiescent tissues, indicating intron inclusion and 3'-untranslated region lengthening as novel mechanisms in regulating aging and dormancy in fungi. This study provides new insights into the complexity and regulation of AS and APA in filamentous fungi.

Dynamics and functional roles of splicing factor autoregulation

Non-core spliceosome components are essential, conserved regulators of alternative splicing. They provide concentration-dependent control of diverse pre-mRNAs. Many splicing factors direct unproductive splicing of their own pre-mRNAs through negative autoregulation. However, the impact of such feedback loops on splicing dynamics at the single-cell level remains unclear. Here, we developed a system to quantitatively analyze negative autoregulatory splicing dynamics by splicing factor SRSF1 in response to perturbations in single HEK293 cells. We show that negative autoregulatory splicing provides critical functions for gene regulation, establishing a ceiling of SRSF1 protein concentration, reducing cell-cell heterogeneity in SRSF1 levels, and buffering variation in transcription. Most important, it adapts SRSF1 splicing activity to variations in demand from other pre-mRNA substrates. A minimal mathematical model of autoregulatory splicing explains these experimentally observed features and provides values for effective biochemical parameters. These results reveal the unique functional roles that splicing negative autoregulation plays in homeostatically regulating transcriptional programs.