Alternative splicing in ABA signaling during seed germination

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

2007; Gehring and Roignant, 2021). In two transesterification reactions catalyzed by splicosome, introns are removed, and exons are linked. In several steps, the spliceosome recognizes conserved sequences at exon-intron junctions and undergoes precise assembly and remodeling (Meyer et al., 2015) ( Figure 1A). The spliceosome functions through numerous RNA and protein interactions and its activity is coordinated by multiple cis-acting elements and trans-acting factors (Punzo et al., 2020;Lin and Zhu, 2021) ( Figure 1B). The constitutive and alternative splicing is distinguished. AS events can be grouped into five general splicing events such as intron retention (IR), exon skipping (ES), mutually exclusive exons, alternative 5′ donor sites (A5′SS) or alternative 3′ acceptor sites (A3′SS) (Kathare and Huq, 2021) ( Figure 2). All introns are excised from pre-mRNA in constitutive splicing, and mature mRNA is assembled from exons. This produces one type of mRNA and generates only one type of protein. While in alternative splicing, introns and exons are removed or included from pre-mRNA in various ways, producing many mRNA isoforms. During exon skipping, the exon is removed from the pre-mRNA together with introns. The intron retention appears when both 5′ and 3′ ends of the specific intron are wrongly recognized, producing mRNA containing that additional intronic sequence. In mutually exclusive exons, only one exon from the cluster of exons remains in the mRNA, while the others are spliced out. When an alternative 5′ splicing site is generated, called a donor site as well, it causes cutting out of the fragment of exon from its 3′ end and nearby intron. The generation of an alternative splice site 3′, which is called the acceptor site as well, causes the cutting out of the intron with the 5′ fragment of the downstream exon (Syed et al., 2012;Chaudhary et al., 2019b). Estimating the number of AS events is demanding because AS is a dynamic process that changes quickly and depends on the type of the cell or tissue, the developmental stage, and environmental conditions (Chaudhary et al., 2019b). However, recent sequencing technologies allowed us much better to Scheme representing the splicing process and its regulation. (A) Steps of spliceosome assembly on precursor messenger RNA (pre-mRNA). Precise pre-mRNA splicing requires four constitutive motifs present in the intronic sequence: the 5′ splice site (GU), the 3′ splice site (AG), the branch point sequence (BPS), and the polypyrimidine tract (Py). The intron excision by the spliceosome is divided into several main stages. In the early step of splicing the E complex is generated. Firstly, splicing factor U1 snRNP U1 small nuclear ribonucleoprotein (U1 snRNP) attaches to the 5′ splice site. Next, SPLICING FACTOR 1 (SF1) recognizes the branch point sequence and U2 snRNP auxiliary factor (U2AF) in the polypyrimidine tract. The A complex is created when U2 snRNP displaces SF1. Then, tri-snRNP U4/U5-U6 combines with complex A to form an inactive B complex. The spliceosome is subjected to ATP-dependent structural changes and two U1 and U4 snRNPs detach leading to spliceosome catalytic activation and forming of an active B act /B* complex. Next, intron excision is catalyzed in two subsequent transesterification reactions. The first reaction occurs in the C complex and the second in the C*complex. In the C*complex exons are joined together and as a result, is formed a post-spliceosomal complex with ligated exons (P complex). The splicing process is completed when the intron is twisted into a loop and excised; the exons are linked together to form mature mRNA, whereas snRNPs are recycled. (B) Role of cis-regulatory sequences and trans-acting splicing regulators in the spliceosome activity. The cis-regulatory sequences include intronic splice enhancers (ISEs) and silencers (ISSs), exonic splice enhancers (ESEs), and silencers (ESSs). To those sequences binds trans-acting RNA binding proteins (RBPs) which promote or inhibit splicing assembly. Serine/arginine (SR) proteins are activators of the splicing that bind ISEs/ESEs and promote exon joining through recruitment of U1 small nuclear ribonucleoprotein (snRNP) to the 5′ splice site and the U2 snRNP to the branch point. On the contrary, the heterogeneous nuclear ribonucleoproteins (hnRNPs) act as splicing inhibitors, attach to ISSs/ESSs, and repress splice site recognition (Gehring and Roignant, 2021). evaluate genes that are found to undergo AS. It is reported that around 95% of multiexon genes of humans and even up to 83% of plant intron-containing genes are subjected to AS (Pan et al., 2008). Plant genes are shorter than animals' and contain fewer exons and shorter introns. These differences in the structure of genes may be the reason for the different number of AS events in plants and animals (Barbazuk et al., 2008;Chang et al., 2017). The most common event of AS in plants is intron retention. In the case of Arabidopsis, this mechanism accounts for approximately 45%-56% of all AS events. Surprisingly, in plants, only 8% of AS are exon skipping events, in contrast to animals, where this AS event is the most common and represents around 58% of all AS events Reddy et al., 2013;Chaudhary et al., 2019b). Proteins generated due to AS may have altered properties or acquire new biological functions, while losing functional domains results in non-functional proteins (Chaudhary et al., 2019a). The change in the proportion of the functional to nonfunctional isoforms fine-tunes the developmental processes and allows the plants quickly react to the changing environmental conditions Ohtani and Wachter, 2019). Alternative splicing also regulates the cell's transcripts level by associating nonsense-mediated mRNA decay (NMD). NMD is the degradation of transcripts containing a premature stop codon (PTC), which causes premature translation termination. The alternatively spliced mRNA isoforms may have a PTC in their nucleotide sequence and then truncated proteins are made from such transcripts. Usually, these proteins are non-functional. Therefore, those aberrant mRNAs are directed and destroyed in the NMD pathway (Drechsel et al., 2013;Nasif et al., 2017). AS has been shown to regulate gene expression through association with the microRNA (miRNA) biogenesis process (Li and Yu, 2021). miRNAs are small and non-coding RNA molecules whose length is around 21-25 nucleotides. They repress the expression of the genes by degrading or inhibiting mRNA translation (Bartel, 2004). AS may occur during primary miRNA transcript (pri-miRNA) processing to mature miRNA. This generates changes in the structure of the pri-miRNA. As a result, the number of functional miRNAs in the cell changes, affecting the expression level of their target genes (Bielewicz et al., 2013;Zhang et al., 2022a). In plants, AS controls the key developmental processes such as growth, circadian clock or flowering time, and adaptation to external conditions, including abiotic stress or pathogen attack (Muhammad et al., 2022). Increasing evidence indicates that AS also regulates seed germination by modulating abscisic acid (ABA) signal transduction (Zhang et al., 2016;Narsai et al., 2017). ABA is the main phytohormone, which through its antagonism with gibberellins (GA), inhibits the process of seed germination (Finch-Savage and Footitt, 2017;Sano and Marion-Poll, 2021). At the seed germination stage, the relevant genes from the ABA signaling pathway, transcription factors, and spliceosome components are alternatively spliced (Table 1).
The alternative splicing modulates ABA signaling components during seed germination HYPERSENSITIVE TO ABA1 (HAB1) phosphatase is one of the better described components of the ABA signaling pathway (Saez et al., 2004;Fujii et al., 2009). HAB1 binds and dephosphorylates SNF1-RELATED PROTEIN KINASE 2 (SnRK2) kinases, especially SnRK2.6 (OPEN STOMATA 1 (OST1), and inhibits the transduction of ABA signal (Umezawa et al., 2009;Vlad et al., 2009). The HAB1 belongs to the PP2C phosphatase family, but only in the case of HAB1 the role of its AS isoforms at the germination stage is known . As a result of AS of HAB1, two protein isoforms are formed that differ in length. The HAB1.1 is four-exonic, and its translation produces a full-length protein. The second isoform, HAB1.2 is created due to intron retention between the third and fourth exon, and this intron sequence contains PTC. This variant produces a truncated protein with a partial catalytic domain because of lacking 105 amino acids at the C-terminal end. This defection changes the phosphatase activity of these two AS variants of HAB1. The HAB1.1 and HAB1.2 can bind to SnRK2, but only the long isoform HAB1.1 can inactivate it by dephosphorylation, which prevents further ABA signal transduction. Notably, the expression level of both HAB1 variants is high in seeds, but what is important is that the regulation of seed Scheme representing constitutive splicing and major alternative splicing events: constitutive splicing, mutually exclusive exons, exon skipping, alternative 5′ splice site, alternative 3′ splice site, and intron retention. Exons are shown as boxes and introns as black lines. Red lines illustrate distinct alternative splicing events (Gurnari et al., 2021). germination in the presence of ABA depends on the HAB1.2/ HAB1.1 isoform ratio. A decrease in the level of HAB1.1 transcripts and an increase in HAB1.2 turn on the ABA signal transduction. On the contrary, an increase in HAB1.1 and a reduction in HAB1.2 turn off the ABA pathway. Studies during seed germination in the presence of ABA using Arabidopsis mutants revealed the role played by HAB1 isoforms in this process. The hab1-1 insertional mutant is hypersensitive to ABA. However, the overexpression of HAB1.1 in the hab1-1 background leads to an ABA-insensitive phenotype, and transgenic seeds can germinate, in contrast to hab1-1. On the other hand, overexpression of HAB1.2 in the hab1-1 background does not complement the germination phenotype to the level of WT. All this indicates that those two HAB1 variants generated in the AS process act antagonistically during seed germination. The role of HAB1.1 is to promote germination, whereas the role of HAB1.2 is to inhibit it . Studies with Arabidopsis mutants in the genes that encode the splicing factors RNA binding motif protein (RMB25), small nuclear ribonucleoprotein E (SmEb), SNW/Skiinteracting protein (SKIP), and ACINUS revealed their involvement in the regulation of AS HAB1 (Cheng et al., 2017;Bi et al., 2021;Hong et al., 2021;Zhang et al., 2022b) (Figure 3). Mutants in these genes are hypersensitive to ABA during seed germination and show an altered expression pattern of HAB1.1 and HAB1.2 isoforms. The RNA binding motif protein (RMB25) contains PWI and RRM motifs important in RNA metabolism and splicing (Zhan et al., 2015). RMB25 interacts with other components of the splicing machinery and controls proper exon junctions (Carlson et al., 2017). ABA activates RBM25 gene expression, and its translation product, the RBM25 protein, attaches to the last intron in the pre-mRNA of HAB1 and generates two HAB1 isoforms, HAB1.1 and HAB1.2. The loss of function of RBM25 increases the HAB1.2/ HAB1.1 isoform ratio in the rbm25 mutant compared to WT Cheng et al., 2017). The SmEb protein is one of the seven Sm core proteins (SmB,SmD3,SmE,SmF,SmG,and SmD1) that are part of the pre-mRNA assembly machinery (Reddy, 2007;Schwer et al., 2016). SmEb contributes to regulating ABA signaling during seed germination through alternative splicing of HAB1 and activating ABA-positive regulators such as the transcription factors ABI3, ABI4, and ABI5 (Hong et al., 2021). Studies using smeb-1 ABI3 splicing during seed germination is regulated by DRT111 and SUA that act in the same pathway. ABI3 homologous in wheat, rice, pea, tomato, flax and alfalfa undergo alternative splicing resulting in isoforms with changed functions. Sugliani et al., 2010;Punzo et al., 2020McKibbin et al., 2002Wilkinson et al., 2005;Fan et al., 2007;Gagete et al., 2009;Gao et al., 2013;Wang et al., 2018;Lalanne et (Zhang et al., 2022b). In response to ABA in the skip-1 mutant, the number of transcripts of functional isoforms of positive ABA regulators PYL7, PYL8, PYL10, ABSCISIC ACID RESPONSIVE ELEMENT-BINDING FACTOR 1 (ABF1), ABF2, and ABI5 are reduced. On the other hand, the level of incorrect intron retention splicing form in negative regulators such as ABI1 and HAB1 is increased. Significantly, SKIP has previously been shown to interact with the SERINE/ARGININE-rich 45 (SR45) spliceosomal protein and control the circadian cycle through the AS of circadian clock genes (Wang et al., 2012). SKIP also controls the flowering time of plants through the AS and is responsible for plant tolerance to salt stress (Feng et al., 2015;Cui et al., 2017). Recently, the homolog of the mammalian apoptotic inducer of nuclear chromatin condensation (ACINUS) was identified in Arabidopsis (Bi et al., 2021). ACINUS, together with the serine-rich domain RNA-binding protein 1 (RNPS1) and the Sin3A-related protein 18 (SAP18), forms a complex of proteins associated with apoptosis and splicing (ASAP) (Schwerk et al., 2003). ACINUS contains an RNA recognition motif (RRM) required for RNA binding. A similar motif present in ACINUS in Arabidopsis has been identified in its distant paralog, another splicing factor called PININ (PNN). Through this motif, PININ also binds to RNPS1 and SAP18 and forms an alternative PSAP complex (Murachelli et al., 2012). Moreover, these two splice factors are also components of the exon junction complex (EJC), which is deposited onto pre-mRNAs to regulate the export of mRNA, translation, and NMD process (Tange et al., 2005;Schlautmann and Gehring, 2020).
Research in Arabidopsis showed that ACINUS and PININ are redundant negative regulators of ABA responses. Knock-out of both ACINUS and PININ results in an abnormal plant phenotype. Double mutant acinus-2 pinin-1 is dwarf and exhibits higher ABA sensitivity and delayed germination and flowering time dynamics (Bi et al., 2021). Transcriptomic analysis of the acinus-2 pinin-1 double mutant exhibited ACINUS and PININ function in the AS mainly by promoting IR events. This type of AS was identified in HAB1 pre-mRNA. In the mutant seedlings, non-functional HAB1.2 was over-accumulated, while the short HAB1.1 isoform was downregulated (Bi et al., 2021). Moreover, not spliced out intron has also been identified in another ABA negative regulator, the ABA HYPERSENSITIVE 1/CAP BINDING PROTEIN 80 (ABH1/CBP80) gene. ABH1/CBP80 and CBP20 are creating the CBC complex (Hugouvieux et al., 2001). Mutations in these genes confer ABA hypersensitive phenotype during seed germination (Hugouvieux et al., 2001;Hugouvieux et al., 2002;Papp et al., 2004;Jäger et al., 2011;Daszkowska-Golec et al., 2013;Daszkowska-Golec et al., 2017). The CBC complex co-transcriptionally interacts with the 5′-end cap structure of pre-mRNA and plays the role of an important splicing regulator (Laubinger et al., 2008;Raczynska et al., 2010). In plants, it has been shown that CBC together with the SERRATE (SE) protein promotes the interaction between the U1 small nuclear ribonucleoprotein (U1 snRNA) factor with the donor splice site of the first intron in the pre-mRNA (Laubinger et al., 2008). Studies with the use of single cbp20 and cbp80(abh1) mutants as well as cbp20/80 double mutant in Arabidopsis showed that a mutation in the CBP80 gene generates more AS events than Scheme showing ABA-dependent splice factors that regulate alternative splicing of genes involved in the regulation of seed germination. Red arrows indicate ABA activation. Black arrows represent alternative splicing promotion. Curved arrows indicate autoregulation of alternative splicing.
Sybilska and Daszkowska-Golec 10.3389/fpls.2023.1144990 Frontiers in Plant Science frontiersin.org mutation in CBP20, which suggests that larger subunit CBP80 may play a decisive role in AS (Raczynska et al., 2010). In the acinus-2 pinin-1 double mutant, the intron-containing isoform of ABH1.2, which produces a non-functional protein, was significantly accumulated, while the splicing form of ABH1.1 was decreased (Bi et al., 2021). Together, all these studies indicate that the observed ABA-hypersensitive phenotype of the mutants in the RMB25, SmEb, SKIP, and ACINUS may relate to the accumulation of defective isoforms of ABA-negative regulators such as HAB1. Thus, these splice factors turn on or off the ABA signal transduction to fine-tune the seed germination by regulating the functional AS isoforms level in the cell.
Transcription factors that regulate seed germination undergo ABArelated alternative splicing AS regulates gene expression by producing various TF isoforms with altered or switched-off DNA binding abilities. This can influence the transcription activation of specific genes and change the transcription rate. One of the better known TFs that undergo AS in plants is ABSCISIC ACID INSENSITIVE 3 (ABI3)/ VIVIPAROUS 1 (VP1). ABI3 plays a significant role in the maturation and germination of the seeds and positively modulates ABA signaling by activating ABA-responsive genes (Lopez-Molina et al., 2002;Finkelstein et al., 2008). Expression of ABI3 is seed-specific and induced by ABA (Parcy et al., 1994;Liu et al., 2013). In Arabidopsis, seeds of abi3 mutant are insensitive to ABA during germination and highly germinate, indicating that ABI3 negatively regulates seed germination in the presence of ABA (Finkelstein et al., 2008). ABI3 induces the expression of MIR159, which encodes the small non-coding RNA miR159. miR159 can regulate ABA signal transduction by targeting two transcription factors MYB DOMAIN PROTEIN 33 (MYB33) and MYB DOMAIN PROTEIN 101 (MYB101), which act as positive regulators of ABA signaling (Reyes and Chua, 2007). Moreover, ABI3 regulates the expression of ABI4 and ABI5, which results in transcription activation of LATE EMBRYOGENESIS ABUNDANT 1 (EM1) and LATE EMBRYOGENESIS ABUNDANT 6 (EM6) genes and leads the inhibition of seed germination (Söderman et al., 2000;Lopez-Molina et al., 2002). Recently, it has been shown that ABI3 can also promote seed dormancy by inducing ABA biosynthesis . ABI3 binds to the promoter of the REVERSAL of RDO5 1 (ODR1), which is a regulator of ABA biosynthesis, and inhibits its expression. This disrupts the interaction of the bHLH57 transcription factor with ODR1 and as a result, allows bHLH57 to activate the transcription of genes from the ABA biosynthetic pathway, such as 9-CIS-EPOXYCAROTENOID DIOXYGENASE 6 (NCED6) and 9-CIS-EPOXYCAROTENOID DIOXYGENASE 9 (NCED9) . However, as well alternative splicing of ABI3 may affect the seed germination process by altering the function of this protein. AS of ABI3 was detected in monocots and dicot plants (McKibbin et al., 2002;Fan et al., 2007;Gagete et al., 2009). ABI3 has four conserved and functional domains in its structure (A1 and B1, B2, B3) (Bies-Etheve et al., 1999;Kurup et al., 2000). The acidic A1 domain is located at the N-terminus and activates gene transcription (McCarty et al., 1991). The B1 domain is necessary for ABI3 interaction with other basic leucine zipper (bZIP) family transcription factors, like ABI5, bZIP10, bZIP25, and TRAB1 (Hobo et al., 1999;Ezcurra et al., 2000;Nakamura et al., 2001;Lara et al., 2003Llorca et al., 2014. In the B2 domain, there is a nuclear localization signal. This domain is also required to transactivate genes containing ABA-responsive elements (ABRE) or G-box motifs (Ezcurra et al., 2000). The B3 domain has DNA binding activity and binds an RY motif (CATGCA), specific for the promoters of genes expressed in the seeds (Suzuki et al., 1997). For the first time, AS of ABI3 was examined in pea (Pisum sativum) (Gagete et al., 2009). In this research were found seven transcripts of PsABI3 the expression of which was detected during seed development, with the highest level mainly in the middle part of the embryogenesis process. The full-length PsABI3 protein containing all functional domains is generated from the PsABI3-1 transcript. Functional PsABI3 protein is also produced from the PsABI3-3 variant, which lacks a fragment located between the A1 and A2 acidic domains. This region contains many amino acids, i.e., proline, serine, and threonine, but its loss does not affect the activity of PsABI3. The PsABI3-1 and PsABI3-3 isoforms can interact with ABI5 due to the presence of the B1 domain. In the case of the PsABI3-2 transcript, the A2 domain, the region between the two acidic domains, and a fragment of the B1 domain are missing, which generates an inactive protein incapable of interacting with ABI5. However, the presence of the B3 domain in PsABI3-2 may cause the ability of this protein to compete with active PsABI3-1 and PsABI3-3 in binding to promoters of target genes and thus prevent transcription of these genes. The PsABI3-5 protein is active despite having only the A1 domain but with an added short amino acid sequence. Notably, the presence of the A1 domain is insufficient for the biological activity of PsABI3. For this reason, the additional amino acid sequence is likely to cause the activity of the PsABI3-5 protein. The remaining three PsABI3-4, PsABI3-6, and PsABI3-7 transcripts produce truncated proteins because of containing a premature termination codon (PTC). Those proteins cannot bind to the promoter sequence of target genes due to the loss of the B3 domain (Gagete et al., 2009). In the model plant Arabidopsis, two ABI3 transcripts have been identified: ABI3-a and ABI3-b. The pre-mRNA of ABI3-a contains a cryptic intron which in translation produces a full-length protein with four functional domains, A1, B1, B2, and B3, which are important for DNA binding. The lack of a cryptic intron in the pre-mRNA of ABI3-b causes the translated protein not to have the B2 and B3 domains and therefore is defective. Analysis of the expression pattern of these AS isoforms revealed that both are expressed in seeds during their development. An increased level of ABI3-a transcript is observed during the entire seed development process. In contrast, an increase in the expression of the ABI3-b isoform is visible in the final stage of seed maturation (Sugliani et al., 2010). In Arabidopsis, it is known that alternative splicing of ABI3 is regulated by the conserved splicing factor SUPPRESSOR OF ABI3-5 (SUA) (Sugliani et al., 2010). SUA is a homolog of the human RNA BINDING MOTIF PROTEIN 5 (RBM5), which binds to the splicing factor U2AF 65 and is responsible for 3′ splice site recognition (Bonnal et al., 2008). The same in Arabidopsis, SUA also interacts with the U2AF 65 subunit and, inhibiting excising of the cryptic intron, promotes the accumulation of the functional ABI3-a isoform. However, the sua-1 mutant has a relatively high level of ABI3-a transcripts indicating that other splicing factors can replace the function of the SUA in AS ABI3 during seed maturation (Sugliani et al., 2010). In the case of tomato (Solanum lycopersicum), the functional ABI3 protein is formed from the S1ABI3-F pre-mRNA (Gao et al., 2013). However, the absence of 90 bp in the pre-mRNA near the region encoding the functional B1 domain results in the formation of a truncated SlABI3-T protein. Interestingly, although S1ABI3-T has all four conserved ABI3 domains, its functioning is disturbed, suggesting the region's important role near the B1 domain for the interaction with ABI5. Only the functional SlABI3-F isoform overexpression in tomato seeds causes hypersensitivity to ABA. This may be caused by the shortening of the SlABI3-T isoform close to the functional B1 domain, which results in a lower binding ability with ABI5 and consequently affects the efficiency of transcription of the ABI5 target genes and disrupts the ABA signaling. Both S1ABI3 transcripts are upregulated mainly during seed imbibition, with the expression of the S1ABI3-F isoform significantly predominant at the germination stage. However, the use of various phytohormones, as well as abiotic stresses, changed the accumulation ratio of the shortened and full-length isoforms. Also, the expression of seed-specific ABA-dependent genes SOMNUS (SlSOM), LATE EMBRYOGENESIS ABUNDANT 1 (SlEM1), and LATE EMBRYOGENESIS ABUNDANT 6 (SlEM6) is regulated differently by both isoforms. The S1ABI3-F increases and S1ABI3-T reduces their expression. Interestingly, studies using transgenic plants with ectopic overexpression of S1ABI3 isoforms in leaves have shown that, in this case, both isoforms function similarly and induce the expression of these ABA-dependent seed specific genes in the leaves. It indicates that the function of S1ABI3-F and S1ABI3-T depends on the tissue in which they are expressed (Gao et al., 2013). The AS of ABI3 has also been shown in flax (Linum usitatissimum) . In this plant, ABI3 produces three transcripts LuABI3-1, LuABI3-2, LuABI3-3. Protein containing all functional domains is produced from those transcripts in the translation process. However, only LuABI3-1 and LuABI3-2 generate functional proteins. The intron retention in the LuABI3-3 pre-mRNA causes the generated protein not to exhibit biological activity; like the research under SlABI3 isoforms in tomatoes and in flax transgenic plants overexpressing LuABI3-1 or LuABI3-2 affects the expression of seed-specific genes and thus regulates seed germination . The latest research on the three stages of seed development of alfalfa (Medicago truncatula) showed the presence of three ABI3 splice isoforms (Lalanne et al., 2021). The first is the full-length MtABI3 protein named SPLICING FORM 1 (SF1), the transcript of which contains nine exons and all functional domains. The SF2 isoform transcript is truncated, contains eight exons, and lacks the A1 activation domain. The pre-mRNA of the third isoform SF3 contains six exons, the A1 domain, and truncated B1 and B3 domains but also lacks the B2 domain. The loss of functional domains in two isoforms caused both to have different cell functions. Moreover, each of these isoforms had a different, specific expression pattern during various phases of seed development and a different and specific pattern of expression depending on the type of tested tissue, i.e., embryo, endosperm, and seed coat. Analysis of the differential expressed genes (DEG) in hairy roots after the ectopic expression of individual splicing forms SF1, SF2, and SF3 in the background of the abi3 mutant revealed that only 41 genes are expressed commonly for all three MtABI3 isoforms. However, 791 specific DEGs for SF1, 357 for SF2, and 177 for SF3 were identified as involved in different biological processes. Specific DEGs in SF1 were associated with photosynthesis, important storage proteins, the photosynthetic PSII complex, and the development of important storage proteins during the seed development process. In the case of SF2, the identified DEGs were associated with the plasma membrane intrinsic protein (PIP) and the degradation of fatty acids and the cell wall. On the other hand, the processes involving differential expressed genes specific to the third SF3 isoform concerned the modification of the cell wall and the activity of pectin methyltransferases. Moreover, alternative splicing of ABSCISIC ACID INSENSITIVE3/ VIVIPAROUS1 (ABI3/VP1) homologs has also been reported in monocots such as wheat (Triticum aestivum) and rice (Oryza sativa) (McKibbin et al., 2002;Wilkinson et al., 2005;Fan et al., 2007). In both species, multiple ABI3 isoforms have been identified that were generated from multiple different AS events, which usually resulted in an open reading frame (ORF) shift and generated a stop codon. Non-functional shortened proteins from aberrantly spliced transcripts are created, which is usually associated with lower seed quality. As a result, the seeds of these crops show reduced primary dormancy as increased pre-harvest sprouting (PHS).
In Arabidopsis, it is known that ABI3 directly interacts with transcription factor ABI5 via the B1 domain and promotes transcription of its target gene by binding to ABRE elements in the DNA (Nakamura et al., 2001). The ABI5 homolog of Arabidopsis in rice OsABI5 is alternatively spliced at the seed germination stage generating two variants of transcripts, the OsABI5-1 and OsABI5-2 (Zou et al., 2007). As a result of the exon skipping the event in the pre-mRNA OsABI5, the translated proteins differ in length in 10 amino acids localized after the bZIP domain. The OsABI5-2 is an elongated protein version; AS changes the binding activity of OsABI5 isoforms with ABI3. The ability to bind to ABI3 of both isoforms is preserved. However, the OsABI5-2 isoform binding is much stronger. This influences the ABI5-ABI3 complex formation efficiency, suggesting different roles of those two ABI5 isoforms in regulating the ABA signaling under changing environmental conditions (Zou et al., 2007). As previously mentioned, two alternatively spliced isoforms of ABI5 have also been identified in Arabidopsis in recent studies of the skip-1 mutant in response to ABA during seed germination (Zhang et al., 2022b).
The PHYTOCHROME INTERACTING FACTORS (PIFs) family consists of the eight proteins belonging to the helix-loophelix (bHLH) TFs that modulate the response of plants to light and act as negative light regulators in the germination process (Leivar and Quail, 2011;Jeong and Choi, 2013;Leivar and Monte, 2014;Yang et al., 2020). PIF1 and PIF6 directly interact with Phytochrome B (PhyB) (Khanna et al., 2004;Narsai et al., 2017). It was revealed that in the presence of light, phytochromes, which are far-red and red-light receptors, are activated and bind to the transcription factor PHYTOCHROME INTERACTING FACTOR 1 (PIF1). This leads to the degradation of PIF1 in ubiquitin and the 26S proteasome pathway and, as a result, enables the seed to germinate (Oh et al., 2004;Shen et al., 2005). In contrast, in the absence of light, the PIF1 protein promotes GA degradation and inhibits the expression of genes from the GA biosynthetic pathway. At the same time, it increases the levels of ABA in seeds by promoting the expression of ABA biosynthesis genes and reducing the expression of genes involved in ABA catabolism (Finkelstein et al., 2008;Oh et al., 2009;Gabriele et al., 2010). In addition to affecting the metabolism of ABA and GA, PIF1 also directly binds to the ABI3 and ABI5 genes from the ABA signal transduction pathway and GA INSENSITIVE DWARF1A (GID1A) and REPRESSOR OF GA1-3 1 (RGA1) genes from the GA signaling pathway (Oh et al., 2007;Oh et al., 2009;Kim et al., 2016a). PIF1 connects the light pathway with the ABA signaling pathway and binds to ABI3 to activate the expression of the negative seed germination regulator SOMNUS (SOM) to inhibit germination of the seeds (Dong et al., 2008;Kim et al., 2008;Park et al., 2011). Regarding the PIF6 transcription factor, it has been shown that it undergoes AS during seed germination. The expression of PIF6 transcripts is particularly observed in dry seeds, and at the initial stage of germination, there is a strong reduction in its expression level. This expression pattern of PIF6 is contrasted with the PIF1, whose transcripts accumulate only at the imbibition stage (Penfield et al., 2010). As a result of the third exon skipping the event, two PIF6 transcripts are generated, encoding two different proteins, PIF6-a and PIF6-b. The PIF6-a variant is a fulllength protein, while the PIF6-b isoform is truncated and lost the helix-loop-helix bHLH domain essential for binding this transcription factor with DNA. Only seeds overexpressing PIF6-b germinate more intensively compared to the germination level of WT, while PIF6-a overexpression does not affect the germination ratio. This shows that AS caused those two PIF6 isoforms to affect the seed germination rate differently. It indicates the antagonistic role of PIF6-a and PIF6-b in the germination process. It is worth mentioning that those two isoforms of PIF6 do not always act differently. Overexpression of PIF6-a or PIF6-b inhibited the elongation of the seedling hypocotyl in R-light, which suggests a similar functioning of these isoforms in seedling development (Penfield et al., 2010). It was revealed that AS pre-mRNA PIF6 is driven by the polypyrimidine tract-binding proteins (PTB), which are the component of the splicosome machinery. PTBs belong to heterogeneous nuclear ribonucleoproteins I (hnRNP I) and act as splicing suppressors (Sawicka et al., 2008;Dai et al., 2022). Arabidopsis genome has three PTB genes (PTB1, PTB2, and PTB3) and all of them undergo alternate splicing to generate two transcripts SPLICING VARIANT I (SPI) and SPLICING VARIANT II (SPII) (Stauffer et al., 2010). The pre-mRNA of SPI is formed into a protein with the correct length, while the SPII has PTC in its nucleotide sequence, which directs this transcript to the NMD pathway. This example shows the role of AS coupling with the NMD mechanism in regulating the level of functional proteins. PTB variants can mutually regulate their splicing but also show the autoregulation ability to balance AS events. Studies in Arabidopsis have shown that AS of PIF6 is regulated only by PTB1 and PTB2 (Rühl et al., 2012). An opposite pattern of accumulation of PIF6 splice isoforms was observed in mutants overexpressing or with knockdown of PTB1 or PTB2. In the case of mutants with PTB1 or PTB2 overexpression and simultaneous overexpression of both PTBs, it led to the accumulation of the correct PIF6-a isoform. In contrast, the knockdown of single PTB1 or PTB2 or both PTBs promoted the truncated PIF6-b isoform. This contrary accumulation pattern of alternatively spliced PIF6 isoforms has been associated with the altered germination ability of PTB mutants in the presence of ABA. Mutants overexpressing PTB1 increased the levels of PIF6-a and delayed seed germination, suggesting that this protein acts as a negative germination regulator in the presence of ABA. In addition, it was confirmed by the observation that the accumulation of non-functional PIF6-b protein in a double mutant with switched off both PTB1 and PTB2 genes increased the seed germination potential in response to ABA. The results of these studies indicate the role of PTBs in the AS of PIF6. Also, other studies in germinating Arabidopsis seeds identified AS of PIF6, in which four differently expressed PIF6 transcripts were detected (Narsai et al., 2017). Moreover, PhyB, which interacts with PIF6, also undergoes AS (Khanna et al., 2004;Narsai et al., 2017). Two PhyB transcripts have been identified that display isoform switch, which suggests that these isoforms may be important at different stages of seed germination. The expression of the first isoform is high in dry seeds and gradually decreases during stratification, while a gradual, high increase in expression is again observed upon exposure to light. In contrast, the expression of the second isoform is the highest during stratification (Narsai et al., 2017). PhyB interacts with the splicing factor DNA-DAMAGE REPAIR/TOLERATION PROTEIN111/SPLICING FACTOR FOR PHYTOCHROME SIGNALING (DRT111/SFPS), which plays a role in photomorphogenesis (Xin et al., 2017). It shows that DRT111 is a hub connecting light and ABA signaling to regulate seed germination (Punzo et al., 2020). The drt111 mutant was observed to have altered expression of genes from the ABA signaling pathway, including ABI3 and the light signaling pathway e.g. PIF1 and PIF6. Analysis of drt111-2/sua-2 double mutants and drt111-2/sua-2/abi3-5 triple mutants showed an interaction between DRT111 and SUA. This indicates that DRT111 controls ABI3 splicing via SUA. Variation in the ABI3 AS isoforms was observed in the single drt111-2 mutant, where only the non-functional ABI3b isoform was increased in dry seeds. Interestingly, at the imbibition stage, both ABI3-a and ABI3-b pre-mRNA amounts increased, but ABI3-b levels rise four times more than WT.
A key regulator of seed dormancy, DOG1 undergoes alternative splicing and interferes with ABA signaling components Interestingly, the main dormancy regulator DELAY OF GERMINATION 1 (DOG1) has been shown to undergo AS, and recent studies have revealed that DOG1 also interacts with the ABA signaling components (Nakabayashi et al., 2015;Dekkers et al., 2016;Neé et al., 2017). As a result of the alternative 3′ and 5′ splice sites selection AS events, five different DOG1 transcripts are generated: DOG1-a, DOG1-b, DOG1-g, DOG1-d, and DOG1-ϵ. These transcripts are translated into three types of proteins because the same protein is produced from the DOG1-b, DOG1-g, and DOG1-ϵ transcripts (Bentsink et al., 2006;Nakabayashi et al., 2015). Although these proteins vary in length, they are all biologically functional. The role of these splice variants of DOG1 is unknown. However, current studies suggest they play a synergistic role in seed dormancy (Nakabayashi et al., 2015). Analysis of the phenotype of mutants with overexpression of the individual single DOG1 isoform in the dog1 mutant background showed that the seeds do not accumulate the DOG1 protein, resulting in no dormant seeds. However, deep dormancy is observed when the mutant overexpresses at least two DOG1 isoforms. This phenotype has been associated with the accumulation of the DOG1 protein, which increases with the number of overexpressed isoforms. This shows that at least two DOG1 isoforms are required to induce seed dormancy. The expression of DOG1 increases with seed development and reaches the highest levels till halfway through the seed maturation process and then decreases until the lowest level in the dry seeds. The rise in the expression of DOG1 during seed maturation is influenced by low temperature and ABA treatment (Chiang et al., 2011;Kendall et al., 2011). The drop in temperature causes a higher amount of DOG1 protein in the cells, which promotes seed dormancy (Zaretskaya et al., 2022). Analysis of the dog1-1 mutant revealed its non-dormant phenotype, which is similar to ABA-deficient and signaling mutants. The dog1-2 cyp707a2-1 double mutant exhibited reduced dormancy compared to the ABA catabolism single mutant cyp707a2. The high ABA levels in the dog1-2 cyp707a2-1 did not lead to seed dormancy (Nakabayashi et al., 2012). However, high levels of DOG1 in the non-dormant ABA-deficient mutant aba1 cannot inhibit the germination of the seeds (Bentsink et al., 2006). Interestingly, recent reports indicate that DOG1 enhances ABA signaling through the interaction and inhibition of PP2C ABA HIPERSENSITIVE GERMINATION1 (AHG1) and ABA HYPERSENSITIVE GERMINATION 3 (AHG3) phosphatases characteristic for seed development and germination (Yoshida et al., 2006;Nishimura et al., 2007;Neé et al., 2017;Nishimura et al., 2018). Studies using mutants in the AHG1 and AHG3, and DOG1 genes have shown that the single ahg1 and ahg3 mutants are significantly less dormant than the ahg1/ahg3 double mutant. On the other hand, double mutants dog1/ahg1 or dog1/ahg3 can still germinate, and only the loss of the function of all three genes causes seed dormancy. This indicates that AGH1 and AHG3 phosphatases may partially take over their functions and that they are downstream from DOG1 (Neé et al., 2017;Nishimura et al., 2018). AHG3 is regulated in the ABA signaling pathway dependent on PYRABACTIN RESISTANCE/PYR-LIKE PROTEIN/ABA RECEPTOR REGULATORY COMPONENTS (PYR/PYL RCAR) receptors and a separate pathway regulated by DOG1. In contrast, AHG1 phosphatase is regulated exclusively by DOG1).
Regulation of ABA-related changes in alternative splicing by the spliceosomal machinery during seed germination Interestingly, it has been shown that spliceosomal components are not only involved in the splicing mechanism but also undergo AS. Serine/arginine-rich (SR) and SR-like proteins belong to the conserved family of pre-mRNA splicing factors. The SR proteins influence the choice of the splicing site and regulate spliceosome assembly by binding to the cis-regulatory sequences in the pre-mRNA (Zahler et al., 1992;Reddy, 2007). In their structure, the SR proteins have one or two RRM motifs at the N-terminus responsible for RNA recognition and an arginine/serine-rich RS domain at the C-terminus allowing for the interaction with other proteins (Wu and Maniatis, 1993;Manley and Tacke, 1996;Graveley et al., 2001;Reddy, 2004;Day et al., 2012). The SR45 protein is an SR-like protein with an additional RS domain located at the N-terminus (Tanabe et al., 2007). In Arabidopsis, SR45 is spliced at various stages of plant development or in response to external stresses (Ali et al., 2007;Zhang et al., 2014;Filichkin et al., 2015). AS for SR45, using alternative 3′ splice sites leads to the production of two isoforms. AS variant SR45.1 is a full-length transcript, while SR45.2 is a truncated version. These isoforms encode very similar proteins differing only in eight amino acids, which are only present in the SR45.1 isoform. Despite this, both isoforms are functional, but surprisingly, each has different functions. The first isoform, SR45.1, plays a primary function in flower petal development, while SR45.2 is in root growth (Ali et al., 2007;Zhang and Mount, 2009;Zhang et al., 2014). Moreover, SR45 has been shown to act in the ABA pathway mediated by glucose (Glc) (Carvalho et al., 2010). Glc induces the expression of ABA biosynthesis and signaling genes, which increases endogenous ABA levels and influences ABA signal transduction (Dekkers et al., 2008;Huang et al., 2016). Such a link between Glc-ABA was observed in the sr45 mutant, which is hypersensitive to glucose and ABA at the seedling stage (Carvalho et al., 2010). The analysis of the ABI3 and ABI5 transcripts levels in the sr45 mutant after treatment with ABA and glucose showed that these transcripts increased significantly compared to the wild type only after glucose treatment. This indicates that the SR45 protein has a role in the negative regulation of ABI3 and ABI5 transcription after glucose treatment. In contrast, reversing the hypersensitivity to Glc of mutant sr45 in complemented lines with SR45.1 and SR45.2 reveals that these isoforms have the same function (Carvalho et al., 2010). Other research studies have shown that the regulatory role of SR45 in the pathway linking sugar signaling with ABA is to control the level of SNF1-RELATED PROTEIN KINASE 1 (SnRK1) in the cell by modulating the AS of INOSITOL POLYPHOSPHATE 5-PHOSPHATASE (5PTase13), which is involved in the proteasomal degradation of SnRK1 (Ananieva et al., 2008;Rodrigues et al., 2013). Analysis of the Arabidopsis seedling transcriptome identified over 4,000 RNAs that directly or indirectly bind to SR45 (Xing et al., 2015). The Gene Ontology (GO) of SR45-related RNA enabled it to indicate the main biological processes in which they are engaged and, through this, confirmed the role of SR45 in processes such as flowering, various aspects of plant development, and RNA splicing. Importantly, the highly enriched term GO responded to abiotic stresses and ABA. This showed that SR45 is involved in the ABA signaling pathway. Moreover, the SR45 knockout affects genes' transcription and the AS from the ABA signaling pathway (Xing et al., 2015). In the case of the HAB1 gene in the sr45 mutant, the transcription of this gene was reduced, but the accumulation of the AS spliced variant HAB1.2 with IR was promoted in the control conditions and in the presence of ABA (Xing et al., 2015). In Arabidopsis, AS isoforms of SR45 also function in response to salt stress (Albaqami et al., 2019). Knock-out of the SR45 resulted in increased salt stress sensitivity at different plant life cycle moments, including seed germination. In the presence of salt stress, sr45 mutant genes related to the SALT OVERLY SENSITIVE (SOS) pathway and ABA signaling were differentially expressed and/or exhibited changed splicing patterns. Interestingly, only the SR45.1 long isoform reversed the salt-sensitive phenotype of the sr45 mutant as well as altered the expression of these genes and splicing patterns. The obtained results indicate that SR45 can positively regulate the response to salt stress in Arabidopsis; in this response, only the longer isoform SR45.1 is involved (Albaqami et al., 2019). SR45a is another splicing regulator that also undergoes AS and, together with the smaller subunit CAP-BINDING PROTEIN 20 (CBP20) of the CAP-BINDING COMPLEX (CBC), cooperates in the salt stress response . SR45a, similar to SR45, has two RS domains, and because of AS, due to the loss of the C-terminal domain, two isoforms with altered functions are created (Tanabe et al., 2009;. SR45a-1a is a full-length functional protein variant, whereas the isoform SR45a-1b is truncated and incapable of interacting with other splicosome elements. However, SR45a-1b plays an important regulatory role. Both SR45a isoforms have the RS domain located at the N-terminal region, so both can interact with the CBP20 protein. In response to salt stress, the complex of SR45a-1a with CBP20 is created, and the binding of those two proteins is enhanced by SR45a-1b. This formed complex regulates the gene expression of the salt-responsive genes and modulates their AS. At the germination stage, the loss of function of SR45a confers tolerance to salinity. At the same time, overexpression of the SR45a-1a or SR45a-1b isoform increased the salinity sensitivity, suggesting a negative regulation of salt stress by the SR45a spliced variants . The latest research also indicates the functionality of SR45 isoforms in seed germination in Arabidopsis. The study of the Arabidopsis transcriptome dynamics during seed germination identified 620 genes that showed isoform variation. Moreover, 612 genes were differentially expressed, including SR45 (Narsai et al., 2017).

Transcriptome analyses reveal the involvement of alternative splicing in seeds
Although the amount of RNA-seq research on AS in plants is constantly growing, more data about seeds and the germination process still needs to be collected, especially in monocots. An RNAseq study of germinating barley (Hordeum vulgare) embryos from four varieties identified 2,200 and 3,900 AS transcripts 24 and 48 h after imbibition (Zhang et al., 2016). The alternative 3′ splice site was the most common AS event, which amounted to 34% and 45% of all AS events. This observation was interesting because the most frequent event in plants is intron retention (Reddy et al., 2012). This research showed that AS process changes over time in germinating seeds. The number of detected transcripts after 24 h and 48 h was similar but not identical. Half of the transcripts present in the embryos after 24 h were no longer detectable at 48 h, but other new transcripts were produced at this time. Evaluating the biological functions of these identified AS genes showed their possible engagement in protein synthesis, energy and metabolism of carbon, and post-transcriptional regulation. Moreover, among AS genes, genes acting in the signaling pathways of various phytohormones have been identified. At the 24 and 48 h of germination, splicing variants of genes from the auxin, cytokinin, and the ABA signaling pathway, such as AUXIN-RESPONSIVE PROTEIN IAA (AUS/IAA) [alternative 3′ splice site (A3S)], A-TYPE RESPONSE REGULATOR (A-ARR) (IR), or SnRK2 (A3S, IR), was detected. However, the A3S alternatively spliced variant of the A-ARR gene was present only for 24 h and the IR isoform of the AUX/ IAA gene after 48 h. Exceptionally high levels of coexpressed genes had auxin efflux carrier, auxin-responsive AUX/IAA proteins, and SnRK2 from the ABA signaling pathway. This indicates that this seed germination stage involves phytohormones such as auxin and ABA. Therefore, the newly-created splice variants may differently modulate various signaling pathways and biological processes and thus affect seed germination. AS was also reported in rice seeds in response to stress conditions . In low oxygen conditions, 1,741 differentially expressed AS were identified. Interestingly, only 5% of AS genes were as well differentially expressed. That result reveals the independent transcription role of alternative splicing in gene regulation. A similar relation was found in the study of the Arabidopsis transcriptome, where only 6% of the identified genes that underwent AS were also differentially expressed (Srinivasan et al., 2016). In 14 and 20 days after pollination, 4,723 and 4,494 AS genes were identified, respectively. Those genes were mainly involved in processes related to mRNA catabolism. In total, 8,927 AS events were detected, of which 88% were newly identified in Arabidopsis. Among them was found a new alternative 3′ splice site (A3P) event of the FUSCA3 (FUS3) gene which codes for a B3 domain transcription factor that acts as a regulator of seed development (Gazzarrini et al., 2004;Tsai and Gazzarrini, 2012). From this AS variant a truncated protein was created due to reading frame shift and formation of a premature stop codon within the B3 domain. This may influence FUS3 functionality and change its DNA binding activity. Although the function of this splice isoform is not determined, it is worth noting that the A3P FUS3 was expressed explicitly at 20 DAP, in contrast to the full-length isoform which was downregulated. The high frequency of the AS process in seeds was also confirmed by RNA-seq studies in dicots (Aghamirzaie et al., 2013;Ruan et al., 2018;Feng et al., 2019;Liu et al., 2022). In the case of peanut (Arachis hypogaea) seeds, 30 days after flowering, 20,213 genes underwent AS, which is 49.69% of all expressed genes. In contrast, on the 50th day after flowering, 19,534 AS genes were identified, corresponding to 48.30% of all expressed genes. 92,483 and 85,562 AS events were detected at these two-time points, respectively. It is worth noting that in the seeds was the most significant number of AS isoforms compared to the root and leaf tissue. It was observed that a considerable percentage of genes that underwent AS was involved in the metabolism of fatty acids (approx. 60%), which is an important process for the development of oil seeds. This suggests that alternative splicing may be a crucial mechanism regulating this process. Also, the new splice variants of Aradu.5N10F, Araip.84LR8, and Araip.92Q2X genes encoding fatty acid desaturase (FAD) were detected. These isoforms were generated as a result of a novel exon presence. However, in the case of the Araip.92Q2X, an additional isoform was also produced which, except for an extra exon, had as well a 3′ alternative exon end (3′-AE). All those newly detected AS events did not disrupt the translation process and resulted in the formation of correct proteins (Ruan et al., 2018). In the soybean (Glycine max) case, 217,371 different AS transcript isoforms from 47,331 AS genes were identified at the stage of embryogenesis (Aghamirzaie et al., 2013). Moreover, RNA-seq studies of germinating seeds showed high AS process dynamics. In cotton (Gossypium australe), the greatest differences in AS were observed in comparison of the 5th and 30th hour of germination (Feng et al., 2019). Similarly, RNAseq studies in Arabidopsis showed variable accumulation of AS isoforms at different phases of seed development and in various tissues, with the most significant changes when comparing seeds after imbibition with dry seeds (Narsai et al., 2017). The AS pattern was recently investigated in seeds of high and low secondary dormancy varieties of oilseed rape (Brassica napus) after polyethylene glycol 6000 (PEG6000) treatment to induce secondary dormancy. After treatment with PEG, 5,136 genes that underwent AS were identified, accounting for 8% of all expressed genes. From these genes, 11,408 splicing variants have been detected, of which over 50% were generated due to intron retention. Among the differentially alternatively spliced genes, a large group consisted of genes encoding spliceosome components, e.g., small nuclear ribonucleoprotein molecules (snRNPs), serine/ arginine-rich proteins (SR) whose function in Arabidopsis has been confirmed in the process of embryo development (JANUS), seed germination (RZ-1C, RBM25) or response to stress conditions (SCL30A, U1-70K, SR45A). Moreover, within the differentially alternatively spliced genes, 342 genes played a role in the secondary dormancy, with only seven of them also significantly expressed. Interestingly, the expression of genes that underwent AS and were also differentially expressed was more than three times higher than the rest of the differentially expressed genes, indicating a correlation between changes in the AS process with changes at the transcription level. However, comparing the number of differentially alternatively spliced genes with differentially expressed genes showed that the number of unique genes is greater than the number of common genes and that suggests an independent action of AS process and transcription in seed dormancy of oilseed rape (Liu et al., 2022).
The RNA-seq data combined with proteomics analysis complements the information about the functions of generated splicing isoforms and may enrich the knowledge of the mechanisms regulating the seed germination process. This approach was used in the study of maize (Zea mays) seed germination in the presence of a high concentration of NaCl. It was noted that during germination, the reaction to the salt stress of the embryo and the endosperm is different. In the salt stress conditions, the embryo could germinate, whereas the germination of seeds with endosperm was inhibited. However, in the presence of ABA, germination was strongly reduced in both cases. These interesting observations suggest that, unlike the endosperm, the embryo is not responsible for the perception of the NaCl, but the signal transduction of the ABA is active. The transcriptomic and proteomics analyses showed significant differences in posttranscriptional events and translational regulation between the embryo and the endosperm. Comparable amounts of differentially expressed posttranscriptional events (DPTE) were detected in the embryo and the endosperm tissue, respectively -14,570 events of 4,529 genes in the embryo and 15,023 events of 4,756 genes in the endosperm. These events included processes that generate different mRNA isoforms, such as the alternative start of transcription, alternative splicing, and alternative polyadenylation. In contrast, few differentially expressed genes were detected in both examined tissues, just 57 in the embryo and 148 in the endosperm. However, proteomics analysis showed the presence of 243 proteins with increased and 419 proteins with reduced expression in the embryo after NaCl treatment. In the case of endosperm, these amounts amounted to 389 and 234, respectively. Interestingly, among the proteins that were produced as a result of DPTE, 11 splicing factors were identified, and their expression at the protein level was differently regulated in the embryo and endosperm. Together, these studies indicate that posttranscriptional events, including AS, may be crucial for seed germination in response to salt stress and that the different expression patterns of splicing factors can be connected with the opposite reaction to NaCl of the embryo and endosperm in maize .

Third-generation sequencing (TGS) in plant transcriptome studies
So far, genome-wide analysis of AS has been primarily based on high-throughput sequencing technology using the RNA-seq Illumina short-read approach. However, detecting various isoforms, repeated sequences, and transposable elements using Illumina is difficult. This is because it is impossible to distinguish different transcripts with identical exons, so reconstructing fulllength alternative transcripts is demanding (Bernard et al., 2014). TGS like Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) overcome these limitations due to the use of SMRT sequencing technology which allows the sequencing of single molecules in real time . The long-read RNA sequencing technologies can sequence full-length transcripts and obtain long cDNA or RNA reads. This allows for identifying new isoforms of transcripts, protein-coding or non-protein-coding genes, and alternative polyadenylation sites. Long reads also remarkably improve the plants' transcriptome annotation (Abdel-Ghany et al., 2016;Zhang et al., 2019;Cui et al., 2020). This together makes PacBio and ONT more often used to analyze plant transcriptomes. TGS has been recently applied in transcriptomic studies of many plant species such as maize (Zea mays) , rice (Oryza sativa) , sorghum (Sorghum bicolor) (Abdel-Ghany et al., 2016), rapeseed (Brassica napus) (Yao et al., 2020), strawberry (Fragaria vesca) , bamboo (Phyllostachys edulis) , Arabidopsis , and tea plant (Camellia sinensis) (Qiao et al., 2019). A maize multi-tissue analysis using PacBio revealed 111,151 transcripts, of which 57% represent a novel, also specific for the tissue isoforms. By discovering new genes, isoforms, and lncRNAs, this method more than doubled the number of alternative transcripts and updated the annotation of the maize genome . To increase assembly quality and reduce the error rate, hybrid sequencing combining NGS and TGS approaches is applied. In the transcriptome study by PacBio and Illumina of adzuki bean (Vigna angularis), two varieties during seed germination under drought stress, 2,457 gene loci, and 46,177 isoforms were newly discovered, and the quality of transcript was significantly improved (Zhu et al., 2020).

Concluding remarks
Alternative splicing (AS) is a crucial co-transcriptional mechanism that regulates gene expression in response to internal cues and environmental signals. Advances in high-throughput sequencing (HTS) technology showed plants' increasing frequency of AS events. However, there still needs to be more information regarding the factors that regulate AS and the effects of AS on mRNA levels and protein functioning. Considering the link between AS and seed germination, further studies are needed to uncover the complex regulatory networks of seed germination and reveal parts of splice variants. This opens another exciting area for research. The obtained knowledge can be used to regulate gene expression to develop crops with new and desirable functional traits.

Author contributions
ES wrote the first draft of the manuscript and prepared the figures and table. AD-G contributed to the conception and design of the manuscript and revised it critically for important intellectual content. All authors contributed to the article and approved the submitted version.

Funding
This work was supported by the National Science Centre, Poland project SONATA BIS10 "(QUEST) Quest for climatesmart barley -the multilayered genomic study of CBC function in ABA signaling" (2020/38/E/NZ9/00346).

Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.