Hostname: page-component-78c5997874-mlc7c Total loading time: 0 Render date: 2024-11-05T04:24:17.443Z Has data issue: false hasContentIssue false

How alternative splicing changes the properties of plant proteins

Published online by Cambridge University Press:  01 July 2022

Ivan Kashkan
Affiliation:
Laboratory of Hormonal Regulations in Plants, Institute of Experimental Botany, Czech Academy of Sciences, Prague, Czech Republic Functional Genomics and Proteomics of Plants, Central European Institute of Technology and National Centre for Biomolecular Research, Masaryk University, Brno 62500, Czech Republic
Ksenia Timofeyenko
Affiliation:
Laboratory of Hormonal Regulations in Plants, Institute of Experimental Botany, Czech Academy of Sciences, Prague, Czech Republic Functional Genomics and Proteomics of Plants, Central European Institute of Technology and National Centre for Biomolecular Research, Masaryk University, Brno 62500, Czech Republic
Kamil Růžička*
Affiliation:
Laboratory of Hormonal Regulations in Plants, Institute of Experimental Botany, Czech Academy of Sciences, Prague, Czech Republic
*
Author for correspondence: K. Růžička, E-mail: [email protected]

Abstract

Most plant primary transcripts undergo alternative splicing (AS), and its impact on protein diversity is a subject of intensive investigation. Several studies have uncovered various mechanisms of how particular protein splice isoforms operate. However, the common principles behind the AS effects on protein function in plants have rarely been surveyed. Here, on the selected examples, we highlight diverse tissue expression patterns, subcellular localization, enzymatic activities, abilities to bind other molecules and other relevant features. We describe how the protein isoforms mutually interact to underline their intriguing roles in altering the functionality of protein complexes. Moreover, we also discuss the known cases when these interactions have been placed inside the autoregulatory loops. This review is particularly intended for plant cell and developmental biologists who would like to gain inspiration on how the splice variants encoded by their genes of interest may coordinately work.

Type
Review
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press in association with The John Innes Centre

1. Introduction

In eukaryotes, almost all protein-coding primary transcripts are interrupted by introns, regions that are not translated and are removed by the process called splicing. The position and length of most introns, their 5 and 3 splice sites, are constitutively defined by the adjoining sequence context. However, some of them can be dynamically adjusted by the process of alternative splicing (AS) (Kelemen et al., Reference Kelemen, Convertini, Zhang, Wen, Shen, Falaleeva and Stamm2013; Lee & Rio, Reference Lee and Rio2015; Reddy et al., Reference Reddy, Marquez, Kalyna and Barta2013; Wang & Burge, Reference Wang and Burge2008). Most AS events occur in common patterns, called AS types or modes (Figure 1a). The most frequent plant AS type is intron retention (Figure 1a), which comprises ~60% of the AS events described (Chaudhary et al., Reference Chaudhary, Jabre, Reddy, Staiger and Syed2019a; Marquez et al., Reference Marquez, Brown, Simpson, Barta and Kalyna2012; Reddy et al., Reference Reddy, Marquez, Kalyna and Barta2013) and includes a particular subtype called exitrons, introns with a protein-coding ability present inside exons, holding ~4% (Marquez et al., Reference Marquez, Höpfler, Ayatollahi, Barta and Kalyna2015). Selection of the alternative 5 or 3 splice sites (Figure 1a) represents ~8 and ~16%, respectively, of Arabidopsis thaliana events. The simultaneous change of both 5 and 3 splice sites, such as exon inclusion or skipping (Figure 1a), corresponds collectively to ~8% of events. The remaining AS types, such as mutual exclusion of exons, rarely occur in plants (Chaudhary et al., Reference Chaudhary, Jabre, Reddy, Staiger and Syed2019a; Filichkin et al., Reference Filichkin, Priest, Megraw and Mockler2015; Marquez et al., Reference Marquez, Brown, Simpson, Barta and Kalyna2012; Martín et al., Reference Martín, Márquez, Mantica, Duque and Irimia2021).

Nearly every plant gene produces, besides the main, usually longest and most expressed, canonical transcript, one or more alternative mRNAs (Marquez et al., Reference Marquez, Brown, Simpson, Barta and Kalyna2012; Zhu et al., Reference Zhu, Chen, Ye, Shi, Ma, Yang, Cao, Zhang, Yoshida, Fernie, Fan, Wen, Zhou, Liu, Fan, Gao, Zhang, Hao, Xiao, Liu and Zhang2017). Although a significant portion of the alternative transcripts seems to be functionally neutral (Tress et al., Reference Tress, Abascal and Valencia2017a; Reference Tress, Abascal and Valencia2017b), it was shown that AS affects the function of countless individual genes (Chaudhary et al., Reference Chaudhary, Khokhar, Jabre, Reddy, Byrne, Wilson and Syed2019b; Filichkin et al., Reference Filichkin, Priest, Megraw and Mockler2015; Kelemen et al., Reference Kelemen, Convertini, Zhang, Wen, Shen, Falaleeva and Stamm2013; Martín et al., Reference Martín, Márquez, Mantica, Duque and Irimia2021; Reddy et al., Reference Reddy, Marquez, Kalyna and Barta2013; Staiger & Brown, Reference Staiger and Brown2013; Szakonyi & Duque, Reference Szakonyi and Duque2018). It has been previously insightfully reviewed how AS is carried out at the (pre-)mRNA level in plants (Reddy et al., Reference Reddy, Marquez, Kalyna and Barta2013), how AS impacts the processing and function of long non-coding RNAs (Fonouni-Farde et al., Reference Fonouni-Farde, Ariel and Crespi2021) and how AS patterns differ between plants and animals (Chaudhary et al., Reference Chaudhary, Khokhar, Jabre, Reddy, Byrne, Wilson and Syed2019b; Martín et al., Reference Martín, Márquez, Mantica, Duque and Irimia2021). The spectrum of physiological and developmental processes related to particular AS events was also discussed from various points of view (Carvalho et al., Reference Carvalho, Feijão and Duque2013; Shang et al., Reference Shang, Cao and Ma2017; Staiger & Brown, Reference Staiger and Brown2013; Szakonyi & Duque, Reference Szakonyi and Duque2018). However, in contrast to several comprehensive reviews in the animal field (Kelemen et al., Reference Kelemen, Convertini, Zhang, Wen, Shen, Falaleeva and Stamm2013; Stamm et al., Reference Stamm, Ben-Ari, Rafalska, Tang, Zhang, Toiber, Thanaraj and Soreq2005), the general functional consequences of AS exerted on the protein level have not been properly summarised in plants. Emphasising the cellular and developmental aspects, we have assembled the most representative and well-characterised AS events from Arabidopsis thaliana and other model systems to illustrate the general mechanistic principles in which the plant splice isoforms appear to coordinately work (Figure 1b).

Fig. 1. A summary of the prominent in-depth characterised AS events in plants. (a) A scheme of common AS types. (b) The list of selected in-depth characterised AS isoforms in plants. Prefixes prior to the protein symbols denote St, Solanum tuberosum, Cr, Catharanthus roseus, Sl, Solanum lycopersicum, no prefix indicates Arabidopsis thaliana. The AS diagrams show exons as boxes, introns as lines; coding sequences are coloured in blue, untranslated regions in grey, and coding regions with alternative reading frames as hatched. NLS, nuclear localization signal; TF, transcription factor. Abbreviations in the Evidence column designate the main experimental approaches supporting the existence of respective splice protein isoforms or proposed molecular models: WB, the native variants detected with western blot; Rep, AS reporter; PS, mRNA associates with polysomes; NX, phenotypes complemented by cDNA expressed under natural promoter; OX, phenotypes conferred by cDNA overexpression; sRi, transcript-specific RNAi; NA, biochemical or microscopic assays performed with cDNAs stably expressed in the native plant systems; HA, in vivo assays performed in heterologous systems, IV, in vitro assays (such as experiments with purified recombinant proteins).

2. Diverse tissue-specific roles of splice isoforms

Transcriptomic studies have demonstrated that a substantial number of splice isoforms show differential expression in various cell types (Klepikova et al., Reference Klepikova, Kasianov, Gerasimov, Logacheva and Penin2016; Li et al., Reference Li, Yamada, Han, Ohler and Benfey2016; Martín et al., Reference Martín, Márquez, Mantica, Duque and Irimia2021). Only limited functional evidence underlies their specific expression patterns and/or distinct ability to rescue mutant loss-of-function phenotypes. The gene determining floral organ size BIG PETAL (BPE) encodes a petal-specific BPEp transcript with the last intron retained, in addition to the canonical BPEub with a uniform expression in all organs (Figures 1b and 2a) (Szécsi et al., 2006). The C-terminal sequence exclusively encoded by BPEp is required for the interaction with the AUXIN RESPONSE FACTOR 8 (ARF8), a component of auxin signalling cascade involved in floral organ development. The arf8 loss-of-function mutants phenocopy the petal defects of the bpe knockout mutants. It was therefore assumed that BPEp, in contrast to the bona fide BPEub protein, is required for petal development (Varaud et al., Reference Varaud, Brioudes, Szécsi, Leroux, Brown, Perrot-Rechenmann and Bendahmane2011; Zhang et al., Reference Zhang, Min, Holappa, Walcher‐Chevillet, Duan, Donaldson, Kong and Kramer2020). ARF8 seems to undergo tissue-specific AS as well. The relative levels of the alternative ARF8.4 variant (showing the in-frame retention of the eighth intron along with the alternative 5 site in the last intron and premature stop codon) were recognised as elevated in flowers (Figures 1b and 2a). The overexpression of the ARF8.4 cDNA, but not other ARF8 isoforms, reverts stamen elongation defects associated with the arf8 knockout mutation (Ghelli et al., Reference Ghelli, Brunetti, Napoli, De Paolis, Cecchetti, Tsuge, Serino, Matsui, Mele, Rinaldi, Palumbo, Barozzi, Costantino and Cardarelli2018). The essential regulator of the early embryogenesis AUXIN RESPONSE FACTOR 5/MONOPTEROS (ARF5/MP) produces an alternative transcript denoted MP11ir with the last (eleventh) intron retained, leading to the protein truncation (Figure 1b). MP11ir cDNA complements ovule-specific defects conferred by the mp/arf5 loss-of-function mutation, in contrast to the canonical ARF5 variant that also rescues the remaining, post-embryonic mp/arf5 defects. Both isoforms show rather similar expression pattern, but it was speculated that the truncated MP11ir protein might be instrumental for an auxin-independent activation of the downstream transcriptional pathways in the early ovule development (Figure 2a; Cucinotta et al., Reference Cucinotta, Cavalleri, Guazzotti, Astori, Manrique, Bombarely, Oliveto, Biffo, Weijers, Kater and Colombo2021).

Fig. 2. Schematic outline of main functional changes of proteins caused by AS in plants. (a) Organ-specific occurrence and functionality of representative splice isoforms. Underlined gene symbols denote the known ability to complement diverse phenotypes. While BPEub, ZIFL1.1 and YUC4.1 generally show uniform expression across organs, the BPEp and YUC4.2 transcripts are enriched in flowers and ZIFL1.3 in leaves. In contrast to ARF8.1, ARF8.4 (presumably elevated in flowers) reverts only flower-specific defects of arf8 knockouts. MP rescues all mp loss-of-function phenotypes, while MP11ir complements only those related to ovule development. Both SR45 variants show practically overlapping expression, but SR45.1 reverts floral organ phenotypes, while SR45.2 rescues root elongation defects conferred by the sr45 loss-of-function mutation. (b) Effect of AS on the subcellular localization of splice isoforms: RAD52-1A (nucleus) and RAD52-1B (mitochondria), RAD52-2A (nucleoplasm and chloroplasts) and RAD52-2B (chloroplasts), TTL1− (peroxisomes) and TTL2− (cytosol), YUC4.1 (together with TAA in the cytosol) and YUC4.2 (together with TAR2 on the outer side of the endoplasmic reticulum), ZIFL1.1 (tonoplast) and ZIFL1.3 (plasma membrane). BES1-S is observed in the nucleus and cytoplasm, but the exclusively nuclearly localized BES1-L is able to displace BES1-S to the nucleus. On the contrary, BRC1aL is seen in the nucleus and BRC1aS in the cytoplasm, when BRC1aS directly prevents nuclear targeting of BRC1aL. PIN7a shows lower mobility within the plasma membrane (and presumably higher transporting capacity in planta) than PIN7b, however, both isoforms associate and can directly affect the mobility of the other protein. NLS, nuclear localization signal; MTS, mitochondrial targeting sequence; TM, transmembrane domain(s); PTS, peroxisomal targeting signal. (c) Examples of common effects of AS on protein function in plants. In case of negative interaction, alternative isoform competes with the canonical protein in the dimer and inhibits its binding to DNA (IDD14, CCA1, SGR5, CO, FLM) or abolishes the formation of the complex required for the catalytic conversion of the metabolite (SGD). Among documented positive interactions, the truncated alternative SR45a-1b isoform enhances the association of the canonical SR45a-1a variant with CBP20, a subunit of the mRNA cap-binding complex. In chloroplasts, the RCA variants coordinately function under different light conditions. (d) A scheme of the positive regulatory loop involving AS of HAB1 in ABA signalling, both positive and negative regulatory loops implicated in AS of JAZ10 in the jasmonate (JA) transduction cascade and the immune activation mediated by AS of CPK28. On the schemes on (b) to (d), the amino acid regions shared between both splice variants are coloured in orange and the regions modified by AS in yellow.

The overexpression of the canonical transcript encoding the membrane transporter ZINC-INDUCED FACILITATOR-LIKE 1 (ZIFL1.1, expressed ubiquitously) rescues several auxin-related defects associated with the zifl1 knockout mutation. Nevertheless, the predominantly leaf-specific ZIFL1.3 variant (originating from the 3 alternative splice site in the 14th intron) (Figures 1b and 2a) reverts the subset of phenotypes linked with the abscisic acid (ABA) and drought response (Remy et al., Reference Remy, Cabrito, Baster, Batista, Teixeira, Friml, Sá-Correia and Duque2013). Similarly, the gene encoding the regulator of AS called SERIN-ARGININE-RICH PROTEIN 45 (SR45) is regulated by the choice of the alternative 3 site in the sixth intron, which removes a short amino acid motif with a critical phosphorylation site (Figure 1b). Both transcripts are expressed across most tissues at comparable levels. However, only the full-length SR45.1 variant can rescue the narrow petal phenotypes observed in the sr45 loss-of-function mutants. SR45.2, on the contrary, complements exclusively their root elongation defects (Figure 2a) (Zhang & Mount, Reference Zhang and Mount2009; Zhang et al., Reference Zhang, Mo, Garrett and Cooper2014).

AS of ARF5, ARF8 and BPE appears to modify the respective targets in a cell-specific manner to tune the tissue or organ identity. In the case of ZIFL1 and SR45, AS intriguingly changes the fundamental functional outcomes of the resulting proteins, and it would be exciting to explore further the mechanistic principles underlying these findings.

3. Differential subcellular localization

One of the most noticeable features of AS is the capability to change the subcellular localization of the protein (Figure 2b). The RADIATION SENSITIVE 52 (RAD52-1) gene, required for the homology-dependent DNA double-strand break repair, encodes a RAD52-1A isoform with the last intron retained (Figure 1b). RAD52-1A was shown to be localized to the nucleoplasm, while the full-length RAD52-1B isoform is targeted to mitochondria (Figure 2b). The concurrent event in the near RAD52-2 paralog gives rise to the RAD52-2A variant, present in both nucleoplasm and chloroplast, and RAD52-2B, detected exclusively in chloroplasts (Figure 2b; Samach et al., Reference Samach, Melamed-Bessudo, Avivi-Ragolski, Pietrokovski and Levy2011). Thus, AS controls the delivery of the RAD52 proteins (and DNA repair) between semiautonomous organelles and nuclei. TRANSTHYRETIN-LIKE (TTL), a protein required for the synthesis of allantoin, is modified by the choice of the alternative 3 site of the last intron (Figure 1b). The resulting isoforms, TTL1− and TTL2−, are localized in peroxisomes and cytoplasm, respectively (Figure 2b) (Lamberto et al., Reference Lamberto, Percudani, Gatti, Folli and Petrucco2010). In addition, the full-length ZIFL1.1 protein (see above) is localized to the tonoplast, while the truncated isoform ZIFL1.3 is targeted to the plasma membrane (Figure 2b; Remy et al., Reference Remy, Cabrito, Baster, Batista, Teixeira, Friml, Sá-Correia and Duque2013).

YUCCA 4 (YUC4) encodes a rate-limiting factor required for auxin biosynthesis. The full-length YUC4.1 protein is localized in the cytosol. The alternative (and flower specific) YUC4.2 variant originates from the transcript showing retention of the last intron (Figure 1b). This region encodes a transmembrane domain that holds the YUC4.2 protein at the cytosolic side of the endoplasmic reticulum (Figure 2b). Both isoforms are catalytically active (Kriechbaumer et al., Reference Kriechbaumer, Wang, Hawes and Abell2012). Later, these observations were placed into a broader context. TRYPTOPHAN AMINOTRANSFERASE OF ARABIDOPSIS 1 (TAA1) and its paralog TRYPTOPHAN AMINOTRANSFERASE RELATED 2 (TAR2) show metabolic activity directly upstream of YUC and reside in the cytosol and on the endoplasmic reticulum, respectively. Moreover, other YUC paralogs (unprocessed by AS) as well localize either in the cytosol or on the endoplasmic reticulum (Figure 2b). Hence, TAA1/TAR2 and the YUC4 isoforms closely associate in both compartments to convert tryptophan to auxin in the formed metabolons (Hrtyan et al., Reference Hrtyan, Šliková, Hejátko and Růžička2015; Reference Kriechbaumer, Botchway and HawesKriechbaumer et al., 2017).

The spatial non-overlapping detachment of the splice isoforms inside the cell apparently indicates that they function independently. Although it seems that there can be cases where AS leads, for example, to deactivation of the protein by its deposition in a different compartment (see also Jiang et al., Reference Jiang, Zhang and Wang2015 and Nicolas et al., Reference Nicolas, Rodríguez-Buey, Franco-Zorrilla and Cubas2015 below), the RAD52-1/2 and YUC4 AS outcomes resemble products of separate genes, analogous to gene duplication.

4. Mutually dependent subcellular distribution

In contrast to the examples where protein isoforms act likely independently, a large part of studies reveals that splice isoforms mutually affect each other’s function. For instance, they can coordinately influence their subcellular localization by direct molecular interaction (Figure 2b). The longer BES1-L variant of the transcriptional factor BES1 (BRI1 EMS SUPPRESSOR 1), involved in brassinosteroid signalling, contains two bipartite nuclear localization signals (NLS), which promote retention of the protein in the nucleus. The shorter BES1-S isoform (designated as canonical due to the sequence conservation) with the alternative transcription initiation codon in the second exon lacks the N-terminal part, including the first NLS, and is observed in both nucleus and cytoplasm (Figure 1b). When BES1-S was co-expressed together with BES1-L, it was detected only in the nucleus, probably due to the dimerisation with BES1-L (Figure 2b). The BES1-mediated relocation to the nucleus was also shown for another component of brassinosteroid signalling, BRASSINAZOLE-RESISTANT 1 (BZR1), which otherwise displays dual cytoplasmic and nuclear localization as well. Both BES1-L and BES1-S isoforms are probably functional. However, the overexpression of BES1-L, but not BES1-S, leads to the phenotypes associated with brassinosteroid (Jiang et al., Reference Jiang, Zhang and Wang2015). The gene BRANCHED1a (BRC1a) codes for a TCP (bHLH) transcription factor, which is processed into two isoforms in the Solanum genus. The nuclear-localized canonical BRC1aL variant with the retained in-frame first intron (more potent at inducing ectopic defects when overexpressed) carries a transcription activation domain on the C-terminus (Figure 1b). Splicing of the first intron in the shorter BRC1aS transcript leads to the replacement of the activation domain by a frame-shifted amino acid sequence and prevents the nuclear targeting of BRC1aS from the cytoplasm. The co-expression of BRC1aL and BRC1aS results in their dimerisation and a partial shift of BRC1aL to the cytosol, along with the decreased ability to induce the reporter-monitored BRC1a-dependent transcription (Figure 2b; Nicolas et al., Reference Nicolas, Rodríguez-Buey, Franco-Zorrilla and Cubas2015).

PIN7 (PIN-FORMED 7), an auxin efflux carrier polarly localized on the plasma membrane, is encoded by two major transcripts. The shorter PIN7b is generated by the choice of an alternative 5 splice site in the first intron (Figure 1b). The resulting protein lacks a 4-amino acid stretch inside the large internal hydrophilic loop (Hrtyan et al., Reference Hrtyan, Šliková, Hejátko and Růžička2015). The longer PIN7a variant, expressed under native promoter, rescues the tropic bending responses and other defects associated with the PIN7 locus, even leading to exaggerated phenotypes. In contrast, PIN7b is almost inactive when expressed alone. Both isoforms show the comparable capability of transporting auxin in a heterologous system and similar subcellular localization in the native tissues. However, tracking with the fluorescence recovery after photobleaching revealed that PIN7a shows lower lateral mobility within the plasma membrane than PIN7b. Moreover, PIN7a and PIN7b form homo- and heterodimers and show the rates of lateral mobility dropping closer to intermediate values when co-expressed (Figure 2b). Consistently, PIN7b reverts the exaggerated tropic response conferred by PIN7a, phenocopying that of the wild-type PIN7 allele (Kashkan et al., Reference Kashkan, Timofeyenko, Kollárová and Růžička2020; Reference Kashkan, Hrtyan, Retzer, Humpolíčková, Jayasree, Filepová, Vondráková, Simon, Rombaut, Jacobs, Frilander, Hejátko, Friml, Petrášek and Růžička2021).

On the outlined examples, the localization overlap marks the likely area where the splice isoforms interact and influence each other’s presence in the given spot. The external cues can tune the resulting activity of the AS products population in the cell. For example, the BRC1a transcript ratios can change following various environmental stimuli (light conditions, decapitation, hormonal treatment) (Nicolas et al., Reference Nicolas, Rodríguez-Buey, Franco-Zorrilla and Cubas2015). Similarly, the levels of PIN7b or both BES1 isoforms can be changed by the application of the respective hormone, likely compensating the response to the growth regulator (Jiang et al., Reference Jiang, Zhang and Wang2015; Kashkan et al., Reference Kashkan, Timofeyenko, Kollárová and Růžička2020; Reference Kashkan, Hrtyan, Retzer, Humpolíčková, Jayasree, Filepová, Vondráková, Simon, Rombaut, Jacobs, Frilander, Hejátko, Friml, Petrášek and Růžička2021).

5. Competitive inhibitory effects

The truncated alternative isoforms commonly show the ability to interfere with the canonical proteins. This has been particularly explored on transcription factors, which tend to form homo- or heterodimers (Seo et al., Reference Seo, Hong, Kim and Park2011a). CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) is a transcriptional factor involved in circadian regulation and cold acclimation. In contrast to the full-length CCA1α, the alternative CCA1β isoform, arising from the alternative initiation codon in the fourth exon, lacks the MYB-type DNA-binding motif (Figure 1b). This can abolish the homodimerisation of CCA1α (and also outcompetes the CCA1α paralog LATE ELONGATED HYPOCOTYL (LHY), dimerising with CCA1α as well), preventing it from binding to the promoters of selected downstream target genes (Figure 2c). Accordingly, the simultaneous presence of the 35S:CCA1β transgene can suppress the 35S:CCA1α overexpression phenotypes. Moreover, the overexpression of a single CCA1α or CCA1β shows opposite effects on the transcription of internal circadian rhythm markers and on the survival rates during cold acclimation (Seo et al., Reference Seo, Park, Lim, Kim, Lee, Baldwin and Park2012).

IDD14 (INDERMINATE DOMAIN 14) is a bHLH transcription factor involved in various morphogenetic processes. Analogously to CCA1, IDD14 encodes an alternative IDD14β isoform arising from the retention of the first intron and alternative initiation codon in the second exon (Figure 1b). Due to the missing DNA binding domain at the N-terminus, IDD14β is inactive. However, it heterodimerises with the canonical IDD14α isoform and inhibits its ability to bind the promoter of the downstream target genes, including QQS (QUA-QUINE STARCH), a factor responsible for the starch degradation (Figure 2c). During cold stress, the proportion of IDD14β increases and the QQS expression is reduced, leading to the elevated starch content (a general indicator of cold acclimation), and these effects can be reverted by the IDD14α overexpression (Seo et al., Reference Seo, Kim, Ryu, Jeong and Park2011b). Furthermore, an AS-mediated mechanism of heat-mediated shoot tropic response was proposed by Kim et al. (Reference Kim, Ryu, Baek and Park2016). A close paralog of IDD14, SHOOT GRAVITROPISM 5 (SGR5 or IDD15) shows a virtually identical isoform interaction scheme, including analogous AS type (Figures 1b and 2c). The sgr5 knockouts show defects in shoot gravitropism, and this phenotype can indeed be rescued by the overexpression of the canonical SGR5α isoform at ambient temperature. The expression of the inhibitory SGR5β isoform increases with growing temperature. In accord with the proposed model, wild type shows reduced shoot gravitropism at increased temperature, while the shoots of the sgr5 knockouts overexpressing the sole SGR5α isoform display practically normal gravitropic bending response even under elevated heat conditions (Kim et al., Reference Kim, Ryu, Baek and Park2016).

CONSTANS or B-BOX DOMAIN PROTEIN 1 (CO or BBX1) is a transcription factor that regulates photoperiodic flowering by controlling the integrator gene FLOWERING LOCUS T (FT). Due to the retention of the only intron and premature stop codon presence, the alternative COβ variant lacks the C-terminal CCT domain responsible for binding DNA (and several other proteins interacting with FT) (Figure 1b). COβ heterodimerises with the full-length COα isoform and prevents it from binding DNA (Figure 2c). Moreover, the presence of COβ in the dimer appears to promote the COα degradation by HOS1 (HIGH EXPRESSION OF OSMOTICALLY RESPONSIVE GENES 1) and COP1 (CONSTITUTIVE PHOTOMORPHOGENIC 1), CO-destabilising ubiquitin E3 ligases, but inhibits its binding to the CO-stabilising E3 ubiquitin ligase FKF1 (F-BOX 1). Thus, the overall COα levels seem to be negatively regulated during the night (HOS1) or in the morning (COP1). In the late afternoon, COα can be temporarily preserved (FKF1), being even protected itself from binding to COβ. The diurnally elevated levels of COα can thereby promote flowering during the long day conditions (Gil et al., Reference Gil, Park, Lee, Park, Han, Kwon, Seo, Jung and Park2017). A similar functional model has also recently been hypothesised for the CO ortholog (Huang et al., Reference Huang, Lin and Wu2022; Jiao & Meyerowitz, Reference Jiao and Meyerowitz2010; Job et al., Reference Job, Yadukrishnan, Bursch, Datta and Johansson2018).

Ultimately, a similar isoform interplay was shown for FLOWERING LOCUS M (FLM, see the scheme of mutual exon exclusion on Figure 1b), a MADS-box transcription factor involved in the regulation of flowering at increased temperature (Lee et al., Reference Lee, Ryu, Chung, Pose, Kim, Schmid and Ahn2013; Posé et al., Reference Posé, Verhage, Ott, Yant, Mathieu, Angenent, Immink and Schmid2013), and parallelised by the FLM paralog MADS AFFECTING FLOWERING 2 (MAF2) (Airoldi et al., Reference Airoldi, McKay and Davies2015). Among splice variants, FLM-δ does not bind DNA but competes with the functional FLM-β isoform for the interaction with the SHORT VEGETATIVE PHASE (SVP) protein, a co-repressor of flowering (Figure 2c). While the levels of FLM-β decrease with the growing temperature, the amounts of FLM-δ rise, releasing the block on the downstream transcripts required for early flowering and the downstream developmental response (Lee et al., Reference Lee, Ryu, Chung, Pose, Kim, Schmid and Ahn2013; Posé et al., Reference Posé, Verhage, Ott, Yant, Mathieu, Angenent, Immink and Schmid2013). The whole mechanism is perhaps more complicated. Further research revealed that a sole decrease of the FLM-β levels is sufficient to induce early flowering (Capovilla et al., Reference Capovilla, Symeonidi, Wu and Schmid2017; John et al., Reference John, Olas and Mueller-Roeber2021; Lutz et al., Reference Lutz, Posé, Pfeifer, Gundlach, Hagmann, Wang, Weigel, Mayer, Schmid and Schwechheimer2015; Reference Lutz, Nussbaumer, Spannagl, Diener, Mayer and Schwechheimer2017; Sureshkumar et al., Reference Sureshkumar, Dent, Seleznev, Tasset and Balasubramanian2016), and the FLM-β amounts at the elevated temperature appear to be lowered by the preferential production of other transcripts that are subsequently degraded by non-sense mediated decay (NMD) (Sureshkumar et al., Reference Sureshkumar, Dent, Seleznev, Tasset and Balasubramanian2016).

Besides interfering with the DNA-binding activity, the dominant-negative alternative isoforms were demonstrated to abolish the catalytic activity of the canonical variants of metabolic enzymes. STRICTOSIDINE β-d-GLUCOSIDASE (SGD) is involved in the synthesis of the cytotoxic monoterpene indole alkaloids in Catharanthus roseus. The alternative variant shSGD lacks a large part of the C-terminal sequence, including NLS, resulting from the retention of the last intron and a premature stop codon (Figure 1b). In contrast to the canonical SGD variant, shSGD is catalytically inactive and unable to self-interact. However, it can heterodimerise with SGD and even disrupts the high-molecular complexes formed by SGD in vitro (Figure 2c). shSGD thereby directly inhibits the enzymatic activity of SGD and affects the synthesis of the relevant monoterpene indole alkaloids in planta. In contrast to the nucleus-resided SGD variant, shSGD shows a dual nuclear and cytosolic localization. In the bimolecular fluorescence complementation interaction assays, it binds also THAS1, another nuclear enzyme involved in further steps of the alkaloid synthesis which normally complexes with SGD. Moreover, shSGD can recruit THAS1 to the cytosol, even when co-expressed with the canonical SGD variant (Carqueijeiro et al., Reference Carqueijeiro, Koudounas, Dugé de Bernonville, Sepúlveda, Mosquera, Bomzan, Oudin, Lanoue, Besseau, Lemos Cruz, Kulagina, Stander, Eymieux, Burlaud-Gaillard, Blanchard, Clastre, Atehortùa, St-Pierre, Giglioli-Guivarc'h, Papon, Nagegowda, O'Connor and Courdavault2021).

A high number of studies illustrate how minor truncated isoforms can interfere with the activity of the full-length proteins. This mode of action seems to be common in most eukaryotes (Jangi & Sharp, Reference Jangi and Sharp2014; Seo et al., Reference Seo, Hong, Kim and Park2011a). Removal of protein domains by AS typically reduces the number of interaction partners at least in half in animal systems (Rodriguez et al., Reference Rodriguez, Pozo, di Domenico, Vazquez and Tress2020; Yang et al., Reference Yang, Coulombe-Huntington, Kang, Sheynkman, Hao, Richardson, Sun, Yang, Shen, Murray, Spirohn, Begg, Duran-Frigola, MacWilliams, Pevzner, Zhong, Trigg, Tam, Ghamsari, Sahni and Vidal2016). Mathematical models of regulatory network motifs indicate that gene expression systems containing dominant-negative factors (here alternative isoforms) show faster response times following the signal stimulus. It can thereby represent a potent adaptation to the changing external or developmental cues (Alon, Reference Alon2007; Jangi & Sharp, Reference Jangi and Sharp2014).

6. Various manners of cooperative action of splice isoforms

Occasionally, the interaction of the canonical and alternative variant(s) can lead to a complex functional interaction in the expressed isoform assemblage. AS of DOG1 (DELAY OF GERMINATION 1), a regulator of seed dormancy, leads to five mRNAs, eventually producing three proteins. If overexpressed, they complement the dog1 loss-of-function phenotypes. However, when these cDNAs were expressed under the natural promoter alone, the resulting proteins were degraded rapidly, failing to restore the dog1 dormancy defects of the mutant entirely (Nakabayashi et al., Reference Nakabayashi, Bartsch, Ding and Soppe2015), or at least moderately (Cyrek et al., Reference Cyrek, Fedak, Ciesielski, Guo, Sliwa, Brzezniak, Krzyczmonik, Pietras, Kaczanowski and Liu2016). The expression of two or more isoforms stabilises by unknown mechanism the subsequent DOG1 accumulation in the nucleus and can rescue the dog1 knockout phenotypes (Nakabayashi et al., Reference Nakabayashi, Bartsch, Ding and Soppe2015).

The gene encoding the MITOGEN-ACTIVATED PROTEIN KINASE 13 (MPK13) gives rise to a truncated alternative transcript MPK13_I4 with the fourth intron retained (Figure 1b). In contrast to the canonical MPK13_Full isoform, MPK13_I4 lacks a part of the kinase domain. MPK13_I4 alone does not show the typical (auto)phosphorylation activity, nor the ability to interact with the MKK6 (MITOGEN-ACTIVATED PROTEIN KINASE 6) acting upstream of MPK13. However, adding the recombinant MPK13_I4 protein into the in vitro reaction mixture enhances the activation of MPK13_Full by MKK6 (Lin et al., Reference Lin, Matsuoka, Sasayama and Nanmori2010). Similarly, the alternative truncated isoform SR45a-1b (different from the SR45 protein above), resulting from the cryptic fifth exon, cannot interact with another core spliceosome component U1-70K due to the lack of the essential C-terminal RNA-binding RS domain. However, it remains partially functional in the salt-stress response linked with the SR45a protein and enhances the formation of the complex of the full-length SR45a-1a isoform and the CBP20 cap-binding protein, along with the regulation of AS of numerous salt-stress related genes (Li et al., Reference Li, Tang, Bassham and Howell2021).

The first step of the chloroplast fixation of CO2 in the Calvin cycle is co-regulated by AS of nuclearly encoded RUBISCO ACTIVASE (RCA). The longer RCAα (or RCAL ) and the shorter RCAβ (RCAS ) transcripts differ in the choice of the 5 splice site in the last intron, leading to the frameshift and protein truncation (Figure 1b; Werneke et al., Reference Werneke, Chatfield and Ogren1989). In multiple species, including Arabidopsis thaliana, both proteins activate Rubisco in vitro. However, the truncated RCAβ lacks cysteine residues required for the perception of fluctuating ADP levels (or of changed redox conditions) and bypasses the feedback loop reacting on the shortage of ATP occurring at the decreased light intensities (Figure 2c; Shen et al., Reference Shen, Orozco and Ogren1991; Zhang & Portis, Reference Zhang and Portis1999). The joint action of both isoforms is a part of adaptation to light conditions: the lines harboring RCAα alone are unable to reach the wild-type rates of Rubisco activation under saturating light conditions. In contrast, the lines carrying exclusively RCAβ show the steadily elevated Rubisco activity, regardless of high- or low-light conditions used. Only the lines containing both RCAα and RCAβ in the rca knockout mutant background display the Rubisco activation dynamics similar to wild type (Zhang et al., Reference Zhang, Kallis, Ewy and Portis2002). It was also shown that the expression levels of both isoforms can be regulated by altering the external temperature or during heat acclimation, and that both isoforms seem to be responsible for different photosynthetic activities under (heat) stress conditions (reviewed in Carvalho et al., Reference Carvalho, Feijão and Duque2013).

Hu et al. (Reference Hu, Mesihovic, Jiménez‐Gómez, Röth, Gebhardt, Bublak, Bovy, Scharf, Schleiff and Fragkostefanakis2020) performed a set of protoplast assays to explore the function of two out of seven variants of the HEAT SHOCK TRANSCRIPTIONAL FACTOR A2 (HsfA2) in tomato. HsfA2-I is the longest isoform. It possesses both nuclear export and NLSs and shuttles between the nucleus and cytoplasm. HsfA2-II carries a cryptic intron towards its 3 terminus, which removes the C-terminal nuclear export signal (Figure 1b). The HsfA2-II protein exhibits a predominantly nuclear localization and decreased protein stability. In contrast to HsfA2-I, HsfA2-II shows a limited ability to interact with the Hsp17.4-CII (HEAT SHOCK PROTEIN 17.4-CII), required for its deposition in the heat shock granules. However, both isoforms can induce transcription of the heat-shock responsive genes. The comparison of the allele polymorphisms further supported the scheme that HsfA2-I can be stored in the heat stress granules over a longer time period and re-used in case of repeated heat exposure, while HsfA2-II can be rather involved in the immediate heat-stress response (Hu et al., Reference Hu, Mesihovic, Jiménez‐Gómez, Röth, Gebhardt, Bublak, Bovy, Scharf, Schleiff and Fragkostefanakis2020).

Several models of how splice variants may coordinately interact have been proposed. The mechanism propounded for the DOG1 protein variants can draw up a situation when multiple isoforms are synergistically required for the correct activity of the resulting protein population. Such systems act as a sign-sensitive filter, which creates a response delay and buffers irregular (stochastic) weak signals, responding only to pronounced stimuli (Alon, Reference Alon2007). Systems containing positive autoregulation, such as MPK13 and SR45a, show slower response time or result in an increased signal variability within the examined cell population, depending on the strength of the input signal (Alon, Reference Alon2007; Jangi & Sharp, Reference Jangi and Sharp2014). The elementary functions of the RCA and HsfA2 are equivalent, but they are adapted to different external cues. This can improve the system robustness in the changing conditions (Alon, Reference Alon2007; Jangi & Sharp, Reference Jangi and Sharp2014). Moreover, the evolutional analysis revealed that the two RCA proteins are encoded by separate genes in some species, conceptually similar to some variants with diverse subcellular localization discussed above (Huang et al., Reference Huang, Lin and Wu2022; Nagarajan & Gill, Reference Nagarajan and Gill2018).

7. Complex autoregulatory circuits tuning splice isoform activity

Several studies uncovered that the splice isoforms participate in positive and negative regulatory loops. These findings integrate the previously outlined basic schemes and illustrate the envisaged complexity of the AS-mediated pathways. The main HsfA2-I isoform of Arabidopsis thaliana contains only two exons, in contrast to the situation in tomato (see above). A mild heat stress under 37°C activates a short cryptic exon splitting the only canonical intron 1 into intron 1a and 1b (Figure 1b). The intervening short sequence introduces a premature stop codon, and the resulting HsfA2-II transcript is eventually not translated, being probably subjected to NMD (Sugio et al., Reference Sugio, Dreos, Aparicio and Maule2009). Under a severe temperature pulse, up to 45°C, the 1a intron incorporates into mRNA as well and gives rise to a leucine-rich motif in the nascent amino acid sequence within the translated HsfA2-III isoform (Figure 1b). The resulting short protein lacks the dimerisation and the C-terminal transactivation domain. However, it contains a partially truncated DNA-binding motif and can bind the heat-shock elements in its own promoter, further promoting HsfA2 expression under extreme heat conditions. HsfA2 thus represents an example of a positive autoregulatory loop (Liu et al., Reference Liu, Sun, Liu, Liu, Du, Wang and Qi2013).

Another positive autoregulatory loop was described for HAB1.2, a truncated isoform resulting from the retention of the last intron of the gene coding for the HAB1 (HYPERSENSITIVE TO ABA 1) phosphatase, a negative regulator of the ABA signalling pathway (Figure 1b). HAB1.2 is ABA inducible and binds the downstream protein kinase OST1 (OPEN STOMATA 1), a positive regulator of ABA response, without the ability to dephosphorylate it. The overexpression of the canonical HAB1.1 transcript in the hab1-1 knockouts leads to the increased resistance to ABA, while HAB1.2 confers the hypersensitivity even exceeding the hab1-1 phenotypes. RBM25 (RNA-BINDING PROTEIN 25), a core regulator of AS, directly binds the last intron of the HAB1 transcript (Figure 2d). Accordingly, the rbm25 loss-of-function mutants show hypersensitivity to ABA and enhanced intron retention rates in several genes, particularly in HAB1. Hence, it was proposed that ABA increases the HAB1.2/HAB1.1 expression ratio with the contribution of RBM25 to keep the ABA signal transduction active (Wang et al., Reference Wang, Ji, Yuan, Wang, Su, Yao, Zhao and Li2015; Zhan et al., Reference Zhan, Qian, Cao, Wu, Yang, Guan, Gu, Wang, Okusolubo, Dunn, Zhu and Zhu2015).

Numerous RNA-binding factors show the ability to bind their own transcript to induce AS, leading to the production of the variants that apparently remain untranslated, thus turning off their own expression (Schöning et al., Reference Schöning, Streitner, Meyer, Gao and Staiger2008; Hartmann et al., Reference Hartmann, Wießner and Wachter2018; Quesada et al., Reference Quesada, Macknight, Dean and Simpson2003). A complex negative auto-regulatory loop, occurring arguably at the protein level, was described for JAZ10 (JASMONATE ZIM DOMAIN PROTEIN 10), a major transcriptional repressor of the nuclear-located jasmonate signalling pathway (Figure 2d). The JAZ10.3 and JAZ10.4 variants are produced by the choice of the alternative 5 site in the last and, respectively, second last intron of the JAZ10 primary transcript (Figure 1b). This results in the partial (JAZ10.3) or complete (JAZ10.4) removal of the conserved Jas motif that under normal conditions binds the MYC2 bHLH transcriptional factors to repress jasmonate-dependent signalling. The Jas motif is recognised by COI1 (CORONATINE INSENSITIVE 1), a F-box protein serving as jasmonate co-receptor, which in the presence of the hormone targets JAZ10 for ubiquitination, leading to the derepression of the MYC2 factors and triggering the downstream response. JAZ10.4 is practically unable to interact with COI1, while the ability of JAZ10.3 to bind COI1 is impaired only partially. Thus, both JAZ10.3 and JAZ10.4 repressors show increased stability following the jasmonate treatment (Chung & Howe, Reference Chung and Howe2009).

Interestingly, the crystallographic studies revealed that JAZ10.4 could bind MYC2 transcription factors even stronger than major JAZ10.1 due to the presence of a cryptic MYC2-interacting domain (CMID) located on its N-terminus (Zhang et al., Reference Zhang, Ke, Zhang, Chen, Sugimoto, Howe, Xu, Zhou, He and Melcher2017). Moreover, while the levels of JAZ10.1 can gradually drop due to the COI1-mediated degradation, the JAZ10.4 expression is induced by the jasmonate treatment. Hence, MYC2 factors, initially derepressed by degradation of JAZ10.1, are bound by JAZ10.4 through the CMID domain and return to the repressed state, attenuating the excessive jasmonate response by a negative feedback loop (Moreno et al., Reference Moreno, Shyu, Campos, Patel, Chung, Yao, He and Howe2013; Zhang et al., Reference Zhang, Ke, Zhang, Chen, Sugimoto, Howe, Xu, Zhou, He and Melcher2017; Figure 2d).

MEDIATOR TRANSCRIPTIONAL COACTIVATOR 25 (MED25), a part of the multimeric Mediator complex, directly binds MYC2 to promote the jasmonate response by recruiting RNA polymerase II to the promoters of the jasmonate-responsive genes (Chung & Howe, Reference Chung and Howe2009; Howe et al., Reference Howe, Major and Koo2018; Yan et al., Reference Yan, Stolz, Chételat, Reymond, Pagni, Dubugnon and Farmer2007; Zhang et al., Reference Zhang, Ke, Zhang, Chen, Sugimoto, Howe, Xu, Zhou, He and Melcher2017). Upon the MYC2 repression by the JAZ proteins, MED25 associates with the jasmonate-inducible PRP39a and PRP40a (PRE-MRNA-PROCESSING FACTOR39a and 40a) splicing factors. They together interact with the JAZ10 primary transcript and shift AS towards the production of the canonical JAZ10.1 mRNA, preventing the excessive desensitisation of jasmonate signalling (Wu et al., Reference Wu, Deng, Zhai, Zhao, Chen and Li2020). Altogether, it seems that AS of JAZ10 can be tuned by both positive and negative feedback loops (Figure 2d).

A thorough experimental effort unraveled the mechanisms accompanying the AS of CALCIUM-DEPENDENT PROTEIN KINASE 28 (CPK28), a negative regulator of plant innate immunity. CPK28 phosphorylates a key positive regulator of plant immunity BIK1 (BOTRYTIS-INDUCED KINASE 1), causing its degradation and attenuation of the downstream immune response. In the absence of the signal associated with the pathogen infection, IMMUNOREGULATORY RNA-BINDING PROTEIN (IRR) is phosphorylated and binds the CPK28 pre-mRNA, activating the preferential splicing of the long, fully functional CPK28 isoform to keep the immunogenic pathways inactive (Figure 1b). Following the immune activation by plant elicitor peptides (Peps), dephosphorylated IRR dissociates from the CPK28 primary transcript, which leads to the preferential expression of the CPK28-RI mRNA with the last three introns retained. CPK28-RI lacks two Ca2+-binding EF-hand domains and shows a severely impaired kinase activity, failing to phosphorylate BIK1. That leads, in turn, to the stabilisation of BIK1 and derepression of the Peps-triggered immune response (Figure 2d) (Dressano et al., Reference Dressano, Weckwerth, Poretsky, Takahashi, Villarreal, Shen, Schroeder, Briggs and Huffaker2020).

8. Limits of our knowledge, future directions

Despite the relatively limited number of elaborated studies, several molecular models have been proposed to manifest the diverse roles of splice isoforms in plants. In essence, they can operate either independently or in a joint manner. Independently acting proteins tend to show different tissue-specific expression or subcellular localization. Here, AS can fundamentally change protein roles (represented by ZIFL1 and SR45) or in effect substitute gene duplication (RAD51, RAD52, YUC4, also RCA). Interaction of splice isoforms, in its turn, represents a level of functional regulation, repressing or modifying the activity of the final protein product(s). The splice isoforms can sometimes mutually influence their subcellular localization (BES1, BRC1a, PIN7). Various examples of positive or coordinated modes of action have been also shown (DOG1, MPK13, SR45a, RCA). Nonetheless, a high number of reports demonstrated functional mechanisms involving dominant-negative (competitive) interaction, particularly on DNA-binding transcription factors (CCA1, IDD14, SGR5, CO, FLM) and also on a metabolic enzyme (SGD). That may perhaps reflect the high rate of intron retention observed in plants (Marquez et al., Reference Marquez, Brown, Simpson, Barta and Kalyna2012). Functionally, it was associated with a rapid reaction to various stimuli. Accordingly, it was revealed that the negative (auto)regulation is the most common network motif in the organismal signalling pathways (Alon, Reference Alon2007; Jangi & Sharp, Reference Jangi and Sharp2014; Lee et al., Reference Lee, Rinaldi, Robert, Odom, Bar-Joseph, Gerber, Hannett, Harbison, Thompson, Simon, Zeitlinger, Jennings, Murray, Gordon, Ren, Wyrick, Tagne, Volkert, Fraenkel, Gifford and Young2002). The feedback loops with negative autoregulation (reported for HsfA2, HAB1, CPK28, and particularly for JAZ10) are not much explored, but they likely accompany many or most of the proposed interaction modes. Works of Shikata et al. (Reference Shikata, Hanada, Ushijima, Nakashima, Suzuki and Matsushita2014) or Huang et al. (Reference Huang, Lin and Wu2022) showed a possible large-scale biology-based direction, how to identify such loops and to integrate them among other signalling pathways.

Surprisingly, the extent to which AS produces the physiologically relevant protein-coding transcripts remains highly debated (Blencowe, Reference Blencowe2017; Tress et al., Reference Tress, Abascal and Valencia2017a; Reference Tress, Abascal and Valencia2017b). Depending on the experimental approach (e.g., transcript association with polysomes, proteomics or evolutional conservation), the predicted share of functionally relevant AS events ranges from far negligible amounts (Abascal et al., Reference Abascal, Ezkurdia, Rodriguez-Rivas, Rodriguez, del Pozo, Vázquez, Valencia and Tress2015; Tress et al., Reference Tress, Abascal and Valencia2017a) to almost half of all expressed transcripts in human and Arabidopsis thaliana (Reixachs-Solé et al., Reference Reixachs-Solé, Ruiz-Orera, Albà and Eyras2020; Weatheritt et al., Reference Weatheritt, Sterne-Weiler and Blencowe2016; Yu et al., Reference Yu, Tian, Yu and Jiao2016). Additionally, many AS events can be specifically activated following specific external stimuli or in small cell groups within particular tissues (Kelemen et al., Reference Kelemen, Convertini, Zhang, Wen, Shen, Falaleeva and Stamm2013; Martín et al., Reference Martín, Márquez, Mantica, Duque and Irimia2021; Reddy et al., Reference Reddy, Marquez, Kalyna and Barta2013; Rodriguez et al., Reference Rodriguez, Pozo, di Domenico, Vazquez and Tress2020). It can be thereby often challenging to confirm their exact functional context in the controlled laboratory condition.

Moreover, a few additional methodological issues have been pointed out, particularly for plant model systems. Current gene schemes, including their protein-coding regions, are based mainly on algorithmic predictions. Hence, many annotated transcripts may not immediately code for proteins, exerting their role at the RNA level. These can also be intermediary products from various stages of mRNA maturation, subject of NMD, or even experimental artefacts. Furthermore, the actual open reading frames can largely differ from the predicted ones as well (Brown et al., Reference Brown, Simpson, Marquez, Gadd, Barta and Kalyna2015). In this context, it is, for example, discussed whether the CCA1 transcripts indeed code for authentic proteins (Brown et al., Reference Brown, Simpson, Marquez, Gadd, Barta and Kalyna2015; Seo et al., Reference Seo, Park, Lim, Kim, Lee, Baldwin and Park2012; Zhang et al., Reference Zhang, Liu, Yuan, Li, Wang, Xu and Xie2021). We have summarised the experimental evidence underlying the natural presence of the outlined protein isoforms, reinforcing the proposed molecular models (Figure 1b). Ideally, the immunoblotting (and the complementation test) has been suggested as solid proof. Additionally, perhaps AS reporter, the association of the transcript with the polyribosome, individually with other indirect data, can serve as a piece of good evidence supporting the authenticity of the protein variant (Brown et al., Reference Brown, Simpson, Marquez, Gadd, Barta and Kalyna2015; Chaudhary et al., Reference Chaudhary, Jabre, Reddy, Staiger and Syed2019a; Kanno et al., Reference Kanno, Venhuizen, Wen, Lin, Chiou, Kalyna, Matzke and Matzke2018; Kashkan et al., Reference Kashkan, Timofeyenko, Kollárová and Růžička2020). Hence, the combined high- and low-scale experimental effort may continuously clear out the current mysteries of the physiological relevance and the most common modi operandi of AS.

Acknowledgements

We thank Elena Zemlyanskaya for her comments on the manuscript. We apologise to the authors whose significant contribution was not discussed due to space limitations.

Financial support

This work was supported by the Ministry of Education, Youth and Sports of the Czech Republic (CZ.02.1.01/0.0/0.0/16_019/0000738) to K.R.

Conflict of interest

The authors declare no conflicts of interest.

Authorship contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Data availability statement

No new data or code are presented in this paper.

Footnotes

I.K. and K.T. contributed equally to this work.

References

Abascal, F., Ezkurdia, I., Rodriguez-Rivas, J., Rodriguez, J. M., del Pozo, A., Vázquez, J., Valencia, A., & Tress, M. L. (2015). Alternatively spliced homologous exons have ancient origins and are highly expressed at the protein level. PLOS Computational Biology, 11, e1004325.CrossRefGoogle ScholarPubMed
Airoldi, C. A., McKay, M., & Davies, B. (2015). MAF2 is regulated by temperature-dependent splicing and represses flowering at low temperatures in parallel with FLM. PLoS One, 10, e0126516.CrossRefGoogle ScholarPubMed
Alon, U. (2007). Network motifs: Theory and experimental approaches. Nature Reviews Genetics, 8, 450461.CrossRefGoogle ScholarPubMed
Blencowe, B. J. (2017). The relationship between alternative splicing and proteomic complexity. Trends in Biochemical Sciences, 42, 407408.CrossRefGoogle ScholarPubMed
Brown, J. W. S., Simpson, C. G., Marquez, Y., Gadd, G. M., Barta, A., & Kalyna, M. (2015). Lost in translation: Pitfalls in deciphering plant alternative splicing transcripts. The Plant Cell, 27, 20832087.CrossRefGoogle ScholarPubMed
Capovilla, G., Symeonidi, E., Wu, R., & Schmid, M. (2017). Contribution of major FLM isoforms to temperature-dependent flowering in Arabidopsis thaliana. Journal of Experimental Botany, 68, 51175127.CrossRefGoogle ScholarPubMed
Carqueijeiro, I., Koudounas, K., Dugé de Bernonville, T., Sepúlveda, LJ., Mosquera, A., Bomzan, DP., Oudin, A., Lanoue, A., Besseau, S., Lemos Cruz, P., Kulagina, N., Stander, EA., Eymieux, S., Burlaud-Gaillard, J., Blanchard, E., Clastre, M., Atehortùa, L., St-Pierre, B., Giglioli-Guivarc'h, N., Papon, N., Nagegowda, DA., O'Connor, SE., Courdavault, V. (2021). Alternative splicing creates a pseudo-strictosidine β-d-glucosidase modulating alkaloid synthesis in Catharanthus roseus. Plant Physiology, 185, 836856.CrossRefGoogle Scholar
Carvalho, R. F., Feijão, C. V., & Duque, P. (2013). On the physiological significance of alternative splicing events in higher plants. Protoplasma, 250, 639650.CrossRefGoogle ScholarPubMed
Chaudhary, S., Jabre, I., Reddy, A. S. N., Staiger, D., & Syed, N. H. (2019a). Perspective on alternative splicing and proteome complexity in plants. Trends in Plant Science, 24, 496506.CrossRefGoogle Scholar
Chaudhary, S., Khokhar, W., Jabre, I., Reddy, A. S. N., Byrne, L. J., Wilson, C. M., & Syed, N. H. (2019b). Alternative splicing and protein diversity: Plants versus animals. Frontiers in Plant Science, 10, 708.CrossRefGoogle Scholar
Chung, H. S., & Howe, G. A. (2009). A critical role for the TIFY motif in repression of Jasmonate signaling by a stabilized splice variant of the JASMONATE ZIM-domain protein JAZ10 in Arabidopsis . The Plant Cell, 21, 131145.CrossRefGoogle ScholarPubMed
Cucinotta, M., Cavalleri, A., Guazzotti, A., Astori, C., Manrique, S., Bombarely, A., Oliveto, S., Biffo, S., Weijers, D., Kater, M. M., Colombo, L. (2021). Alternative splicing generates a MONOPTEROS isoform required for ovule development. Current Biology, 31, 892899.e3.CrossRefGoogle ScholarPubMed
Cyrek, M., Fedak, H., Ciesielski, A., Guo, Y., Sliwa, A., Brzezniak, L., Krzyczmonik, K., Pietras, Z., Kaczanowski, S., Liu, F., et al. (2016). Seed dormancy in Arabidopsis is controlled by alternative polyadenylation of DOG1 . Plant Physiology, 170, 947955.CrossRefGoogle ScholarPubMed
Dressano, K., Weckwerth, P. R., Poretsky, E., Takahashi, Y., Villarreal, C., Shen, Z., Schroeder, J. I., Briggs, S. P., & Huffaker, A. (2020). Dynamic regulation of pep-induced immunity through post-translational control of defence transcript splicing. Nature Plants, 6, 10081019.CrossRefGoogle ScholarPubMed
Filichkin, S., Priest, H. D., Megraw, M., & Mockler, T. C. (2015). Alternative splicing in plants: Directing traffic at the crossroads of adaptation and environmental stress. Current Opinion in Plant Biology, 24, 125135.CrossRefGoogle ScholarPubMed
Fonouni-Farde, C., Ariel, F., & Crespi, M. (2021). Plant long noncoding RNAs: New players in the field of post-transcriptional regulations. Non-Coding RNA, 7, 12.CrossRefGoogle ScholarPubMed
Ghelli, R., Brunetti, P., Napoli, N., De Paolis, A., Cecchetti, V., Tsuge, T., Serino, G., Matsui, M., Mele, G., Rinaldi, G., Palumbo, GA., Barozzi, F., Costantino, P., Cardarelli, M. (2018). A newly identified flower-specific splice variant of AUXIN RESPONSE FACTOR8 regulates stamen elongation and endothecium lignification in Arabidopsis. The Plant Cell, 30, 620637.CrossRefGoogle ScholarPubMed
Gil, K.-E., Park, M.-J., Lee, H.-J., Park, Y.-J., Han, S.-H., Kwon, Y.-J., Seo, P. J., Jung, J.-H., & Park, C.-M. (2017). Alternative splicing provides a proactive mechanism for the diurnal CONSTANS dynamics in Arabidopsis photoperiodic flowering. The Plant Journal, 89, 128140.CrossRefGoogle ScholarPubMed
Hartmann, L., Wießner, T., & Wachter, A. (2018). Subcellular compartmentation of alternatively spliced transcripts defines SERINE/ARGININE-RICH PROTEIN30 expression. Plant Physiology, 176, 28862903.CrossRefGoogle ScholarPubMed
Howe, G. A., Major, I. T., & Koo, A. J. (2018). Modularity in Jasmonate signaling for multistress resilience. Annual Review of Plant Biology, 69, 387415.CrossRefGoogle ScholarPubMed
Hrtyan, M., Šliková, E., Hejátko, J., & Růžička, K. (2015). RNA processing in auxin and cytokinin pathways. Journal of Experimental Botany, 66, 48974912.CrossRefGoogle ScholarPubMed
Hu, Y., Mesihovic, A., Jiménez‐Gómez, J. M., Röth, S., Gebhardt, P., Bublak, D., Bovy, A., Scharf, K., Schleiff, E., & Fragkostefanakis, S. (2020). Natural variation in HsfA2 pre‐mRNA splicing is associated with changes in thermotolerance during tomato domestication. New Phytologist, 225, 12971310.CrossRefGoogle ScholarPubMed
Huang, C.-K., Lin, W.-D., & Wu, S.-H. (2022). An improved repertoire of splicing variants and their potential roles in Arabidopsis photomorphogenic development. Genome Biology, 23, 50.CrossRefGoogle ScholarPubMed
Jangi, M., & Sharp, P. A. (2014). Building robust transcriptomes with master splicing factors. Cell, 159, 487498.CrossRefGoogle ScholarPubMed
Jiang, J., Zhang, C., & Wang, X. (2015). A recently evolved isoform of the transcription factor BES1 promotes Brassinosteroid signaling and development in Arabidopsis thaliana . The Plant Cell, 27, 361374.CrossRefGoogle ScholarPubMed
Jiao, Y., & Meyerowitz, E. M. (2010). Cell‐type specific analysis of translating RNAs in developing flowers reveals new levels of control. Molecular Systems Biology, 6, 419.CrossRefGoogle ScholarPubMed
Job, N., Yadukrishnan, P., Bursch, K., Datta, S., & Johansson, H. (2018). Two B-box proteins regulate Photomorphogenesis by oppositely modulating HY5 through their diverse C-terminal domains. Plant Physiology, 176, 29632976.CrossRefGoogle ScholarPubMed
John, S., Olas, J. J., & Mueller-Roeber, B. (2021). Regulation of alternative splicing in response to temperature variation in plants. Journal of Experimental Botany, 72, 61506163.CrossRefGoogle ScholarPubMed
Kanno, T., Venhuizen, P., Wen, T.-N., Lin, W.-D., Chiou, P., Kalyna, M., Matzke, A. J. M., & Matzke, M. (2018). PRP4KA, a putative spliceosomal protein kinase, is important for alternative splicing and development in Arabidopsis thaliana . Genetics, 210, 12671285.CrossRefGoogle ScholarPubMed
Kashkan, I., Hrtyan, M., Retzer, K., Humpolíčková, J., Jayasree, A., Filepová, R., Vondráková, Z., Simon, S., Rombaut, D., Jacobs, T. B., Frilander, M. J., Hejátko, J., Friml, J., Petrášek, J., & Růžička, K. (2021). Mutually opposing activity of PIN7 splicing isoforms is required for auxin‐mediated tropic responses in Arabidopsis thaliana. New Phytologist, 233, 329343.CrossRefGoogle ScholarPubMed
Kashkan, I., Timofeyenko, K., Kollárová, E., & Růžička, K. (2020). In vivo reporters for visualizing alternative splicing of hormonal genes. Plants, 9, 868.CrossRefGoogle ScholarPubMed
Kelemen, O., Convertini, P., Zhang, Z., Wen, Y., Shen, M., Falaleeva, M., & Stamm, S. (2013). Function of alternative splicing. Gene, 514, 130.CrossRefGoogle ScholarPubMed
Kim, J., Ryu, J. Y., Baek, K., & Park, C. (2016). High temperature attenuates the gravitropism of inflorescence stems by inducing SHOOT GRAVITROPISM 5 alternative splicing in Arabidopsis rabidopsis . New Phytologist, 209, 265279.CrossRefGoogle Scholar
Klepikova, A. V., Kasianov, A. S., Gerasimov, E. S., Logacheva, M. D., & Penin, A. A. (2016). A high resolution map of the Arabidopsis thaliana developmental transcriptome based on RNA-seq profiling. The Plant Journal, 88, 10581070.CrossRefGoogle ScholarPubMed
Kriechbaumer, V., Botchway, S. W., & Hawes, C. (2017). Localization and interactions between Arabidopsis auxin biosynthetic enzymes in the TAA/YUC-dependent pathway. Journal of Experimental Botany, 68, 41954207.Google Scholar
Kriechbaumer, V., Wang, P., Hawes, C., & Abell, B. M. (2012). Alternative splicing of the auxin biosynthesis gene YUCCA4 determines its subcellular compartmentation: YUCCA4 and auxin biosynthesis. The Plant Journal, 70, 292302.CrossRefGoogle ScholarPubMed
Lamberto, I., Percudani, R., Gatti, R., Folli, C., & Petrucco, S. (2010). Conserved alternative splicing of Arabidopsis transthyretin-like determines protein localization and S-Allantoin synthesis in peroxisomes. The Plant Cell, 22, 15641574.CrossRefGoogle ScholarPubMed
Lee, T. I., Rinaldi, N. J., Robert, F., Odom, D. T., Bar-Joseph, Z., Gerber, G. K., Hannett, N. M., Harbison, C. T., Thompson, C. M., Simon, I., Zeitlinger, J., Jennings, E. G., Murray, H. L., Gordon, D. B., Ren, B., Wyrick, J. J., Tagne, J. B., Volkert, T. L., Fraenkel, E., Gifford, D. K., Young, R. A. 2002. Transcriptional regulatory networks in saccharomyces cerevisiae. Science, 298, 799804.CrossRefGoogle ScholarPubMed
Lee, Y., & Rio, D. C. (2015). Mechanisms and regulation of alternative pre-mRNA splicing. Annual Review of Biochemistry, 84, 291323.CrossRefGoogle ScholarPubMed
Lee, J. H., Ryu, H.-S., Chung, K. S., Pose, D., Kim, S., Schmid, M., & Ahn, J. H. (2013). Regulation of temperature-responsive flowering by MADS-box transcription factor repressors. Science, 342, 628632.CrossRefGoogle ScholarPubMed
Li, Z., Tang, J., Bassham, D. C., & Howell, S. H. (2021). Daily temperature cycles promote alternative splicing of RNAs encoding SR45a, a splicing regulator in maize. Plant Physiology, 186, 13181335.CrossRefGoogle ScholarPubMed
Li, S., Yamada, M., Han, X., Ohler, U., & Benfey, P. N. (2016). High-resolution expression map of the Arabidopsis root reveals alternative splicing and lincRNA regulation. Developmental Cell, 39, 508522.CrossRefGoogle ScholarPubMed
Lin, W.-Y., Matsuoka, D., Sasayama, D., & Nanmori, T. (2010). A splice variant of Arabidopsis mitogen-activated protein kinase and its regulatory function in the MKK6–MPK13 pathway. Plant Science, 178, 245250.CrossRefGoogle Scholar
Liu, J., Sun, N., Liu, M., Liu, J., Du, B., Wang, X., & Qi, X. (2013). An autoregulatory loop controlling Arabidopsis HsfA2 expression: Role of heat shock-induced alternative splicing. Plant Physiology, 162, 512521.CrossRefGoogle ScholarPubMed
Lutz, U., Nussbaumer, T., Spannagl, M., Diener, J., Mayer, K. F., & Schwechheimer, C. (2017). Natural haplotypes of FLM non-coding sequences fine-tune flowering time in ambient spring temperatures in Arabidopsis. eLife, 6, e22114.CrossRefGoogle ScholarPubMed
Lutz, U., Posé, D., Pfeifer, M., Gundlach, H., Hagmann, J., Wang, C., Weigel, D., Mayer, K. F. X., Schmid, M., & Schwechheimer, C. (2015). Modulation of ambient temperature-dependent flowering in Arabidopsis thaliana by natural variation of FLOWERING LOCUS M. PLoS Genetics, 11, e1005588.CrossRefGoogle ScholarPubMed
Marquez, Y., Brown, J. W. S., Simpson, C., Barta, A., & Kalyna, M. (2012). Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis . Genome Research, 22, 11841195.CrossRefGoogle ScholarPubMed
Marquez, Y., Höpfler, M., Ayatollahi, Z., Barta, A., & Kalyna, M. (2015). Unmasking alternative splicing inside protein-coding exons defines exitrons and their role in proteome plasticity. Genome Research, 25, 9951007.CrossRefGoogle ScholarPubMed
Martín, G., Márquez, Y., Mantica, F., Duque, P., & Irimia, M. (2021). Alternative splicing landscapes in Arabidopsis thaliana across tissues and stress conditions highlight major functional differences with animals. Genome Biology, 22, 35.CrossRefGoogle ScholarPubMed
Moreno, J. E., Shyu, C., Campos, M. L., Patel, L. C., Chung, H. S., Yao, J., He, S. Y., & Howe, G. A. (2013). Negative feedback control of jasmonate signaling by an alternative splice variant of JAZ10. Plant Physiology, 162, 10061017.CrossRefGoogle ScholarPubMed
Nagarajan, R., & Gill, K. S. (2018). Evolution of rubisco activase gene in plants. Plant Molecular Biology, 96, 6987.CrossRefGoogle ScholarPubMed
Nakabayashi, K., Bartsch, M., Ding, J., & Soppe, W. J. J. (2015). Seed dormancy in Arabidopsis requires self-Binding ability of DOG1 protein and the presence of multiple isoforms generated by alternative splicing. PLoS Genetics, 11, e1005737.CrossRefGoogle ScholarPubMed
Nicolas, M., Rodríguez-Buey, M. L., Franco-Zorrilla, J. M., & Cubas, P. (2015). A recently evolved alternative splice site in the BRANCHED1a gene controls potato plant architecture. Current Biology, 25, 17991809.CrossRefGoogle ScholarPubMed
Posé, D., Verhage, L., Ott, F., Yant, L., Mathieu, J., Angenent, G. C., Immink, R. G. H., & Schmid, M. (2013). Temperature-dependent regulation of flowering by antagonistic FLM variants. Nature, 503, 414417.CrossRefGoogle ScholarPubMed
Quesada, V., Macknight, R., Dean, C., & Simpson, G. G. (2003) Autoregulation of FCA pre-mRNA processing controls Arabidopsis owering time. The EMBO Journal, 22, 31423152.CrossRefGoogle Scholar
Reddy, A. S. N., Marquez, Y., Kalyna, M., & Barta, A. (2013). Complexity of the alternative splicing landscape in plants. The Plant Cell, 25, 36573683.CrossRefGoogle ScholarPubMed
Reixachs-Solé, M., Ruiz-Orera, J., Albà, M. M., & Eyras, E. (2020). Ribosome profiling at isoform level reveals evolutionary conserved impacts of differential splicing on the proteome. Nature Communications, 11, 1768.CrossRefGoogle ScholarPubMed
Remy, E., Cabrito, T. R., Baster, P., Batista, R. A., Teixeira, M. C., Friml, J., Sá-Correia, I., & Duque, P. (2013). A major facilitator superfamily transporter plays a dual role in polar auxin transport and drought stress tolerance in Arabidopsis . The Plant Cell, 25, 901926.CrossRefGoogle Scholar
Rodriguez, J. M., Pozo, F., di Domenico, T., Vazquez, J., & Tress, M. L. (2020). An analysis of tissue-specific alternative splicing at the protein level. PLoS Computational Biology, 16, e1008287.CrossRefGoogle ScholarPubMed
Samach, A., Melamed-Bessudo, C., Avivi-Ragolski, N., Pietrokovski, S., & Levy, A. A. (2011). Identification of plant RAD52 homologs and characterization of the Arabidopsis thaliana RAD52 -like genes. The Plant Cell, 23, 42664279.CrossRefGoogle ScholarPubMed
Schöning, J. C., Streitner, C., Meyer, I. M., Gao, Y., & Staiger, D. (2008). Reciprocal regulation of glycine-rich RNA-binding proteins via an interlocked feedback loop coupling alternative splicing to nonsense-mediated decay in Arabidopsis. Nucleic Acids Research, 36, 69776987.CrossRefGoogle ScholarPubMed
Seo, P. J., Hong, S.-Y., Kim, S.-G., & Park, C.-M. (2011a). Competitive inhibition of transcription factors by small interfering peptides. Trends in Plant Science, 16, 541549.CrossRefGoogle Scholar
Seo, P. J., Kim, M. J., Ryu, J.-Y., Jeong, E.-Y., & Park, C.-M. (2011b). Two splice variants of the IDD14 transcription factor competitively form nonfunctional heterodimers which may regulate starch metabolism. Nature Communications, 2, 303.CrossRefGoogle Scholar
Seo, P. J., Park, M.-J., Lim, M.-H., Kim, S.-G., Lee, M., Baldwin, I. T., & Park, C.-M. (2012). A self-regulatory circuit of CIRCADIAN CLOCK-ASSOCIATED1 underlies the circadian clock regulation of temperature responses in Arabidopsis . The Plant Cell, 24, 24272442.CrossRefGoogle ScholarPubMed
Shang, X., Cao, Y., & Ma, L. (2017). Alternative splicing in plant genes: A means of regulating the environmental fitness of plants. International Journal of Molecular Sciences, 18, 432.CrossRefGoogle Scholar
Shen, J. B., Orozco, E. M., & Ogren, W. L. (1991). Expression of the two isoforms of spinach ribulose 1,5-bisphosphate carboxylase activase and essentiality of the conserved lysine in the consensus nucleotide-binding domain. Journal of Biological Chemistry, 266, 89638968.CrossRefGoogle ScholarPubMed
Shikata, H., Hanada, K., Ushijima, T., Nakashima, M., Suzuki, Y., & Matsushita, T. (2014). Phytochrome controls alternative splicing to mediate light responses in Arabidopsis . Proceedings of the National Academy of Sciences, 111, 1878118786.CrossRefGoogle ScholarPubMed
Staiger, D., & Brown, J. W. S. (2013). Alternative splicing at the intersection of biological timing, development, and stress responses. The Plant Cell, 25, 36403656.CrossRefGoogle ScholarPubMed
Stamm, S., Ben-Ari, S., Rafalska, I., Tang, Y., Zhang, Z., Toiber, D., Thanaraj, T. A., & Soreq, H. (2005). Function of alternative splicing. Gene, 344, 120.CrossRefGoogle ScholarPubMed
Sugio, A., Dreos, R., Aparicio, F., & Maule, A. J. (2009). The cytosolic protein response as a subcomponent of the wider heat shock response in Arabidopsis . The Plant Cell, 21, 642654.CrossRefGoogle ScholarPubMed
Sureshkumar, S., Dent, C., Seleznev, A., Tasset, C., & Balasubramanian, S. (2016). Nonsense-mediated mRNA decay modulates FLM-dependent thermosensory flowering response in Arabidopsis. Nature Plants, 2, 16055.CrossRefGoogle ScholarPubMed
Szakonyi, D., & Duque, P. (2018). Alternative splicing as a regulator of early plant development. Frontiers in Plant Science, 9, 1174.CrossRefGoogle ScholarPubMed
Szécsi, J., Joly, C., Bordji, K., Varaud, E., Cock, J. M., Dumas, C., & Bendahmane, M. (2006). BIGPETALp, a bHLH transcription factor is involved in the control ofArabidopsis petal size. EMBO Journal, 25, 39123920.CrossRefGoogle Scholar
Tress, M. L., Abascal, F., & Valencia, A. (2017a). Alternative splicing may not be the key to proteome complexity. Trends in Biochemical Sciences, 42, 98110.CrossRefGoogle Scholar
Tress, M. L., Abascal, F., & Valencia, A. (2017b). Most alternative isoforms are not functionally important. Trends in Biochemical Sciences, 42, 408410.CrossRefGoogle Scholar
Varaud, E., Brioudes, F., Szécsi, J., Leroux, J., Brown, S., Perrot-Rechenmann, C., & Bendahmane, M. (2011). AUXIN RESPONSE FACTOR8 regulates Arabidopsis petal growth by interacting with the bHLH transcription FACTOR BIGPETALp. The Plant Cell, 23, 973983.CrossRefGoogle ScholarPubMed
Wang, Z., & Burge, C. B. (2008). Splicing regulation: From a parts list of regulatory elements to an integrated splicing code. RNA, 14, 802813.CrossRefGoogle Scholar
Wang, Z., Ji, H., Yuan, B., Wang, S., Su, C., Yao, B., Zhao, H., & Li, X. (2015). ABA signalling is fine-tuned by antagonistic HAB1 variants. Nature Communications, 6, 8138.CrossRefGoogle ScholarPubMed
Weatheritt, R. J., Sterne-Weiler, T., & Blencowe, B. J. (2016). The ribosome-engaged landscape of alternative splicing. Nature Structural & Molecular Biology, 23, 11171123.CrossRefGoogle ScholarPubMed
Werneke, J. M., Chatfield, J. M., & Ogren, W. L. (1989) Alternative mRNA splicing generates the two ribulosebisphosphate carboxylase/oxygenase activase polypeptides in spinach and Arabídopsis. 1, 815825.Google Scholar
Wu, F., Deng, L., Zhai, Q., Zhao, J., Chen, Q., & Li, C. (2020). Mediator subunit MED25 couples alternative splicing of JAZ genes with fine-tuning of Jasmonate signaling. The Plant Cell, 32, 429448.CrossRefGoogle ScholarPubMed
Yan, Y., Stolz, S., Chételat, A., Reymond, P., Pagni, M., Dubugnon, L., & Farmer, E. E. (2007). A downstream mediator in the growth repression limb of the jasmonate pathway. The Plant Cell, 19, 24702483.CrossRefGoogle ScholarPubMed
Yang, X., Coulombe-Huntington, J., Kang, S., Sheynkman, G. M., Hao, T., Richardson, A., Sun, S., Yang, F., Shen, Y. A., Murray, R. R., Spirohn, K., Begg, B. E., Duran-Frigola, M., MacWilliams, A., Pevzner, S. J., Zhong, Q., Trigg, S. A., Tam, S., Ghamsari, L., Sahni, N., Vidal, M. (2016). Widespread expansion of protein interaction capabilities by alternative splicing. Cell, 164, 805817.CrossRefGoogle ScholarPubMed
Yu, H., Tian, C., Yu, Y., & Jiao, Y. (2016). Transcriptome survey of the contribution of alternative splicing to proteome diversity in Arabidopsis thaliana. Molecular Plant, 9, 749752.CrossRefGoogle ScholarPubMed
Zhan, X., Qian, B., Cao, F., Wu, W., Yang, L., Guan, Q., Gu, X., Wang, P., Okusolubo, T. A., Dunn, S. L., Zhu, J.-K., & Zhu, J. (2015). An Arabidopsis PWI and RRM motif-containing protein is critical for pre-mRNA splicing and ABA responses. Nature Communications, 6, 8139.CrossRefGoogle ScholarPubMed
Zhang, N., Kallis, R. P., Ewy, R. G., & Portis, A. R. (2002). Light modulation of rubisco in Arabidopsis requires a capacity for redox regulation of the larger rubisco activase isoform. Proceedings of the National Academy of Sciences, 99, 33303334.CrossRefGoogle ScholarPubMed
Zhang, F., Ke, J., Zhang, L., Chen, R., Sugimoto, K., Howe, G. A., Xu, H. E., Zhou, M., He, S. Y., & Melcher, K. (2017). Structural insights into alternative splicing-mediated desensitization of jasmonate signaling. Proceedings of the National Academy of Sciences, 114, 17201725.CrossRefGoogle ScholarPubMed
Zhang, S., Liu, H., Yuan, L., Li, X., Wang, L., Xu, X., & Xie, Q. (2021). Recognition of CCA1 alternative protein isoforms during temperature acclimation. Plant Cell Reports, 40, 421432.CrossRefGoogle ScholarPubMed
Zhang, R., Min, Y., Holappa, L. D., Walcher‐Chevillet, C. L., Duan, X., Donaldson, E., Kong, H., & Kramer, E. M. (2020). A role for the auxin response factors ARF6 and ARF8 homologs in petal spur elongation and nectary maturation in aquilegia . New Phytologist, 227, 13921405.CrossRefGoogle ScholarPubMed
Zhang, X.-N., Mo, C., Garrett, W. M., & Cooper, B. (2014). Phosphothreonine 218 is required for the function of SR45.1 in regulating flower petal development in Arabidopsis . Plant Signaling & Behavior, 9, e29134.CrossRefGoogle ScholarPubMed
Zhang, X.-N., & Mount, S. M. (2009). Two alternatively spliced isoforms of the Arabidopsis SR45 protein have distinct roles during normal plant development. Plant Physiology, 150, 14501458.CrossRefGoogle ScholarPubMed
Zhang, N., & Portis, A. R. (1999). Mechanism of light regulation of rubisco: A specific role for the larger rubisco activase isoform involving reductive activation by thioredoxin-f. Proceedings of the National Academy of Sciences, 96, 94389443.CrossRefGoogle ScholarPubMed
Zhu, F.-Y., Chen, M.-X., Ye, N.-H., Shi, L., Ma, K.-L., Yang, J.-F., Cao, Y.-Y., Zhang, Y., Yoshida, T., Fernie, A. R., Fan, G. Y., Wen, B., Zhou, R., Liu, T. Y., Fan, T., Gao, B., Zhang, D., Hao, G. F., Xiao, S., Liu, Y. G., Zhang, J. (2017). Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings. The Plant Journal, 91, 518533.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. A summary of the prominent in-depth characterised AS events in plants. (a) A scheme of common AS types. (b) The list of selected in-depth characterised AS isoforms in plants. Prefixes prior to the protein symbols denote St, Solanum tuberosum, Cr, Catharanthus roseus, Sl, Solanum lycopersicum, no prefix indicates Arabidopsis thaliana. The AS diagrams show exons as boxes, introns as lines; coding sequences are coloured in blue, untranslated regions in grey, and coding regions with alternative reading frames as hatched. NLS, nuclear localization signal; TF, transcription factor. Abbreviations in the Evidence column designate the main experimental approaches supporting the existence of respective splice protein isoforms or proposed molecular models: WB, the native variants detected with western blot; Rep, AS reporter; PS, mRNA associates with polysomes; NX, phenotypes complemented by cDNA expressed under natural promoter; OX, phenotypes conferred by cDNA overexpression; sRi, transcript-specific RNAi; NA, biochemical or microscopic assays performed with cDNAs stably expressed in the native plant systems; HA, in vivo assays performed in heterologous systems, IV, in vitro assays (such as experiments with purified recombinant proteins).

Figure 1

Fig. 2. Schematic outline of main functional changes of proteins caused by AS in plants. (a) Organ-specific occurrence and functionality of representative splice isoforms. Underlined gene symbols denote the known ability to complement diverse phenotypes. While BPEub, ZIFL1.1 and YUC4.1 generally show uniform expression across organs, the BPEp and YUC4.2 transcripts are enriched in flowers and ZIFL1.3 in leaves. In contrast to ARF8.1, ARF8.4 (presumably elevated in flowers) reverts only flower-specific defects of arf8 knockouts. MP rescues all mp loss-of-function phenotypes, while MP11ir complements only those related to ovule development. Both SR45 variants show practically overlapping expression, but SR45.1 reverts floral organ phenotypes, while SR45.2 rescues root elongation defects conferred by the sr45 loss-of-function mutation. (b) Effect of AS on the subcellular localization of splice isoforms: RAD52-1A (nucleus) and RAD52-1B (mitochondria), RAD52-2A (nucleoplasm and chloroplasts) and RAD52-2B (chloroplasts), TTL1− (peroxisomes) and TTL2− (cytosol), YUC4.1 (together with TAA in the cytosol) and YUC4.2 (together with TAR2 on the outer side of the endoplasmic reticulum), ZIFL1.1 (tonoplast) and ZIFL1.3 (plasma membrane). BES1-S is observed in the nucleus and cytoplasm, but the exclusively nuclearly localized BES1-L is able to displace BES1-S to the nucleus. On the contrary, BRC1aL is seen in the nucleus and BRC1aS in the cytoplasm, when BRC1aS directly prevents nuclear targeting of BRC1aL. PIN7a shows lower mobility within the plasma membrane (and presumably higher transporting capacity in planta) than PIN7b, however, both isoforms associate and can directly affect the mobility of the other protein. NLS, nuclear localization signal; MTS, mitochondrial targeting sequence; TM, transmembrane domain(s); PTS, peroxisomal targeting signal. (c) Examples of common effects of AS on protein function in plants. In case of negative interaction, alternative isoform competes with the canonical protein in the dimer and inhibits its binding to DNA (IDD14, CCA1, SGR5, CO, FLM) or abolishes the formation of the complex required for the catalytic conversion of the metabolite (SGD). Among documented positive interactions, the truncated alternative SR45a-1b isoform enhances the association of the canonical SR45a-1a variant with CBP20, a subunit of the mRNA cap-binding complex. In chloroplasts, the RCA variants coordinately function under different light conditions. (d) A scheme of the positive regulatory loop involving AS of HAB1 in ABA signalling, both positive and negative regulatory loops implicated in AS of JAZ10 in the jasmonate (JA) transduction cascade and the immune activation mediated by AS of CPK28. On the schemes on (b) to (d), the amino acid regions shared between both splice variants are coloured in orange and the regions modified by AS in yellow.

Author comment: How alternative splicing changes the properties of plant proteins — R0/PR1

Comments

Editor Prague, 4 November 2021

Quantitative Plant Biology

Dear Editor,

Following the invitation to the Quantitative Plant Biology communicated with Dr. Olivier Hamant, please find the submission of the review article manuscript entitled How alternative splicing changes the properties of plant proteins.

The research on alternative splicing (AS) represents a dynamically expanding research field and the role of AS in plants has been repeatedly summarized from various points of view. Interestingly, none of the previous contributions attempted to conceptualize how plant AS isoforms mechanistically work, similar to the seminal review articles from the animal field,.

Respecting the uniqueness of model plant systems, we attempted to fill this gap. We gathered the prominent, well-characterized AS events that arose recently in the literature and assembled them in a similar manner as in the referred animal reviews. We thus believe that our study would serve as a comprehensive guide for any researcher interested in the general role of AS in plants.

Please see below the following experts suggested as potential referees. Though, we would be grateful for, if you do not send this manuscript to Gordon Simpson (University of Dundee, UK).

Dr. Maria Kalyna

University of Natural Resources and Life Sciences, Austria

[email protected]

Expertise: molecular mechanisms of alternative splicing, transcriptomics

Assoc. Prof. Misato Ohtani

University of Tokyo, Japan

[email protected]

Expertise: RNA processing in plant development

Prof. John Brown

University of Dundee, UK

[email protected]

Expertise: alternative splicing in plants, transcriptomics

Dr. Paula Duque

University of Lisbon, Portugal

[email protected]

Expertise: alternative splicing in plants

Prof. Artur Jarmolowski

Adam Mickiewicz University, Poland

[email protected]

Expertise: alternative splicing and RNA processing in plants

Dr. Craig Simpson

The James Hutton Institute, UK

[email protected]

Expertise: mechanisms of alternative splicing in plants

Once more, we thank for the considering our contribution and the opportunity to publish in Quantitative Plant Biology.

We apologize for the delay with preparing the manuscript.

With very best wishes,

Kamil Ruzicka, Ivan Kashkan and Ksenia Timofeyenko

Institute of Experimental Botany

Czech Academy of Sciences

Rozvojová 263

165 02 Praha 6 - Lysolaje

Czechia

Review: How alternative splicing changes the properties of plant proteins — R0/PR2

Conflict of interest statement

Reviewer declares none.

Comments

Comments to Author: This review attempts to give a comprehensive summary related to alternative splicing (AS) in plants. This is a timely and meaningful work which can be of interest to the readers of QPB. Overall the review is well written and the structure is clear. I have, however, few main comments aiming to make the review more accessible to a wide audience. A major missing part is some basic descriptions of AS in plants for non-initiated readers who lack knowledge about AS. While the author described a lot studies in details, most sections lack a clear summary and take home message. In the present form, this makes it difficult for readers to capture the key points of each section.

In the following, I provide several broad areas to improve this review article towards being more reader friendly.

1. The concept of alternative splicing

This review directly expounds the impact of AS on the transcriptome at the beginning. But it would be nice to provide an accurate definition of AS. Although biologists are generally familiar with alternative splicing, definitions and general mechanisms would still be necessary for the initial part of the review (even if citing authoritative reviews on AS as done in the introduction). Accurate descriptions can help beginners to understand alternative splicing very well.

2. Types (Modes) and mechanisms of alternative splicing

Following on the description of AS; it would be nice to see a description of the assembly mechanism of alternative splicing in the process of gene expression. The authors could first introduce the concept of AS, and then describe how alternative splicing is spliced. My suggestion is to describe the assembly process of AS according to the different assembly modes of AS (ie. exon skipping, mutually exclusive exons, alternative donor site, and so on), and draw a schematic diagram for each AS mode, which can help readers understand the process of AS more clearly.

3. The author's summary is missing in each section

To make the review go beyond the description of previous work on AS, it would be nice to have a summary of these studies combined with your own understanding as take home message. I absolutely believe that the authors have read a lot of articles to prepare this review. However, a summary of each section is lacking. This could be done also by including some future research prospects for each section.

4. Tissue-specific roles of splice isoforms

The section 2 (Diverse tissue-specific roles of splice isoforms) introduces the tissue specificity of splice isoforms. The author only describes the specificity of many splice isoforms in floral organs, I think that splice isoforms in other tissues or organs ( for example in root system(Li et al. 2016)) could also be introduced in this section.

5. Alternative splicing under stressful conditions

Adaptation of plant for stressful condition (ie. cold, drought and radiation) is one of the focal points (also include defense against biotic stress). There is big interest about roles of AS under different conditions for botanists/experimentalists. Some more general overview of research on this aspect (Li et al. 2013, Yang et al. 2014, Shang et al. 2017, Punzo et al. 2020) could be included to generalize the end of section 7. This section lists already some studies related to roles of AS under light and heat stimulation, but could include or generalize to other stresses.

6. I could suggest that the author integrate together the section 3 and section 4 (subcellular localization), as well as section 5 and section 6 (interaction between splice isoforms).

7. Overall the The conclusions are somehow redundant with each section and could be better and not refined and include a view of the current limitations in AS studies and and prospect for future research and future questions in the field.. A list of the 5 most outstanding questions in the field of plant AS could for example be included as outlook in the conclusion.

Review: How alternative splicing changes the properties of plant proteins — R0/PR3

Conflict of interest statement

Reviewer declares none.

Comments

Comments to Author: Kashkan et al. reviewed the effects of alternative splicing (AS) on protein functionality, based on the experimental evidences. I of course agree with the importance of AS to regulate plant cell activity/function through changes in protein function/localization/activity, and this review article should be important to make the readers aware the diversified roles of AS-based protein regulation. However, unfortunately, the current version of manuscript is simply “a list of examples”. There is no clear and original perspective of advantages and/or physiological significance of AS-mediated regulation of protein functions in plants, I felt.

1. For making this review article more attractive for general readers, I suggest to reconstruct the examples of AS: in the current manuscript, the description was sometimes mixed and jumped, between molecular effects on proteins and physiological aspects, such as stress response (response to exogenous stimuli). It is so confusing; the cases shown in section 7 can be also mentioned in section 2-6, for example.

Could you re-organize the examples of AS genes/proteins perfectly based on the molecular effects? I think the best way is to follow your graphic abstract as the systematic grouping of cases for effects of AS on proteins. Or, more simply, you can insert a systematic (logical) view of molecular effect by AS (based on graphic abstract) as a summary, in the beginning of manuscript (just after Introduction). It helps the readers to understand your concept more clearly.

2. Please put the discussion on physiological significance of AS-mediated regulation of protein functions, before (or within) the section of Conclusion. No original perspective was given in your manuscript. In the discussion, you can also mention the hidden effects of AS for changing protein functionality, such as translational efficiency, mRNA stability, and mRNA localization, and which parts are possible “plant-specific aspects”.

3. In the abstract, you said “but the common underlying functional principles have surprisingly never been generalized.” What did you mean by “common underlying functional principles”? I think 1) the mode of action for AS is well studied, and its basis is widely common among eukaryotes (i.e. specific trans- and cis-factors can act as regulators), and 2) the resultant changes in protein structures/localization by AS are well summarized already; for example, the differential protein structure can change protein functionality apparently. The complicated circuits of AS regulation including splicing factors are also well known in the sex determination of insects. It should be better to make clear what is “the common underlying functional principles” in this review.

4. In Figure 1: the information in the “Mechanism” column was not well unified. For example, “tissue specificity” is not actually mechanism and not at the same level with “subcellular localization.” Rather, “tissue specificity” should be put in the same line with “temperature sensitive ” and “light intensity sensitive”, as their responsiveness. As well, you have to spell out NLS in the caption.

Recommendation: How alternative splicing changes the properties of plant proteins — R0/PR4

Comments

Comments to Author: Dear Dr Růžička,

Many thanks for submitting to Quantitative Plant Biology.

Your review article addresses an important and interesting topic for readers of QPB, however the reviewers have concerns on the aims, structure and organization of this review article. I recommend therefore the submisison of new revised version.

After reading the review article myself I concur with the reviewers that the article at the moment reads somehow as a list of examples but clear links between them or clear take home messages are still lacking. The article could be improved by reshaping some sections and by making the structure clearer (see detailed comments by the reviewers), for example by having a first section on the bases and mechanisms of AS in plants. From there you could highlight what the main questions in the field currently are, this would yield your graphical abstract. I find that the graphical abstract is a good idea but it is not optimally used, while it could be nicely used to guide non-expert readers. This section would lay down a clear plan for the next sections. Each section would consist of several examples and a clear take home message.

I would suggest to add a “Contribution of AS to quantitative biology” section before the outlook/conclusion to show the relevance of AS links for quantitative plant biology (the aim of the journal). When reading the article I got curious about the role and importance of AS in linking genomics and phenotypes. For example what is the relevance of AS for quantitative traits (important for plant breeding), for QTL mapping/GWAs, and for allele specific expression (for important phenotypes such as abiotic stress tolerance, yield, heterosis effects). This section would allow to generalize why AS is important for crop breeding or plant evolution in general. At present most examples seem to stem from A. thaliana, but it is not always clearly stated. Maybe this section would highlight the relevance of AS in crops or in evolution of plant species in general.

The final conclusion/outlook section could highlight few major unanswered questions in the field that need to be addressed with the increasing number of large genomic, transcriptomic, metabolomic, proteomic and phenotypic data. This is one of the main aim of QPB and it would be great to get a sense of what the authors consider important for the future of big data analysis to tackle the relevance/importance/mechanisms of AS in plants.

Many thanks for submitting your work to QPB, and we look forward to receive a revised version of your article.

Sincerely,

Prof. Aurélien Tellier

Decision: How alternative splicing changes the properties of plant proteins — R0/PR5

Comments

No accompanying comment.

Author comment: How alternative splicing changes the properties of plant proteins — R1/PR6

Comments

11 March 2022

Dear Aurélien and Olivier,

Please find enclosed our revised manuscript entitled How alternative splicing changes the properties of plant proteins (QPB-21-0060).

We have discussed the points raised and revised our manuscript accordingly. We believe that the current form of the manuscript meets better the journal scope and appears to be more attractive for potential readers. The detailed response to the editor and referees’ comments is attached.

We greatly appreciate the interest of Quantitative Plant Biology in our work. We also thank the editorial team for the generousity in conceiving the revision of the manuscript and the kindly granted increase of the word limit.

With very best wishes,

Kamil Ruzicka, Ivan Kashkan and Ksenia Timofeyenko

Institute of Experimental Botany

Czech Academy of Sciences

Rozvojová 263

165 02 Praha 6 - Lysolaje

Czechia

Review: How alternative splicing changes the properties of plant proteins — R1/PR7

Comments

Comments to Author: The authors have addressed my comments, I thus recommend the acceptance of this manuscript for the publication after a few small changes.

Although the authors provide full names for some important AS isoforms or genes, there are still many genes that have only gene symbols, not full names. It would be easy to read and check if the authors provide a table including the full and abbreviated names of all AS isoforms mentioned in this manuscript.

It would be more precise for two points as follow:

Line 60. different change to differential

Line 383. tissue change to tissue-specific

Review: How alternative splicing changes the properties of plant proteins — R1/PR8

Comments

Comments to Author: The revised manuscript seemed to be improved, and now it better matches what QPB is looking for in a review article. Before acceptance, I would like to ask 2 minor points,

1. The authors are using “Arabidopsis” to refer to Arabidopsis thaliana. AS events can vary widely among plant species, so in my sense, the authors should clearly limit the specie to Arabidopsis thaliana. At least when it first appears, please write the species name as “Arabidopsis thaliana.”

2. Please mention about “exitron” in the introduction briefly; it is very important form of intron, and it has been shown that the high rate of exitron is a very unique feature of plants as you know. It unnatural that there is no mention of exitron when discussing AS in plants.

Recommendation: How alternative splicing changes the properties of plant proteins — R1/PR9

Comments

Comments to Author:

Dear authors,

I apologize for the time to get the reviews in this round. I concur with the reviewers that the paper has been greatly improved, and gained a lot in clarity. Both reviewers recommend publication pending some last minor changes.

Reviewer 1 suggests to make sure that full names of genes/isoforms are written somewhere (in text or table) and not only the gene symbols.

Reviewer 2 suggests to be precise when to use Arabidopsis or Arabidopsis thaliana, and to mention exitron briefly in the introduction.

Providing you address these minor changes, the manuscript will be accepted for publication.

I look forward to receive your revised version soon.

Best regards,

Aurelien Tellier

Decision: How alternative splicing changes the properties of plant proteins — R1/PR10

Comments

No accompanying comment.

Author comment: How alternative splicing changes the properties of plant proteins — R2/PR11

Comments

Dear Aurélien and Olivier,

Please find enclosed our second revision of the manuscript How alternative splicing changes the properties of plant proteins (QPB-21-0060). We have corrected the remaining minor points as requested.

We thank you for the opportunity to contribute to QPB and for the helpful editorial and reviewer comments. We hope our work will be of interest to many prospective readers!

With very best wishes,

Kamil Ruzicka, Ivan Kashkan and Ksenia Timofeyenko

Laboratory of Hormonal Regulations in Plants

Prague, 29 April 2022

Recommendation: How alternative splicing changes the properties of plant proteins — R2/PR12

Comments

No accompanying comment.

Decision: How alternative splicing changes the properties of plant proteins — R2/PR13

Comments

No accompanying comment.