Alternative splicing at NAGNAG acceptors in Arabidopsis thaliana SR and SR-related protein-coding genes

Alternative splicing responses to salt stress in Glycyrrhiza uralensis revealed by global profiling of transcriptome RNA-seq datasets

Excessive reactive oxygen species stress due to salinity poses a significant threat to the growth of Glycyrrhiza uralensis Fisch. To adapt to salt stress, G. uralensis engages in alternative splicing (AS) to generate a variety of proteins that help it withstand the effects of salt stress. While several studies have investigated the impact of alternative splicing on plants stress responses, the mechanisms by which AS interacts with transcriptional regulation to modulate the salt stress response in G. uralensis remain poorly understood. In this study, we utilized high-throughput RNA sequencing data to perform a comprehensive analysis of AS events at various time points in G. uralensis under salt stress, with exon skipping (SE) being the predominant AS type. KEGG enrichment analysis was performed on the different splicing genes (DSG), and pathways associated with AS were significantly enriched, including RNA transport, mRNA surveillance, and spliceosome. This indicated splicing regulation of genes, resulting in AS events under salt stress conditions. Moreover, plant response to salt stress pathways were also enriched, such as mitogen-activated protein kinase signaling pathway - plant, flavonoid biosynthesis, and oxidative phosphorylation. We focused on four differentially significant genes in the MAPK pathway by AS and qRT-PCR analysis. The alternative splicing type of MPK4 and SnRK2 was skipped exon (SE). ETR2 and RbohD were retained intron (RI) and alternative 5'splice site (A5SS), respectively. The expression levels of isoform1 of these four genes displayed different but significant increases in different tissue sites and salt stress treatment times. These findings suggest that MPK4, SnRK2, ETR2, and RbohD in G. uralensis activate the expression of isoform1, leading to the production of more isoform1 protein and thereby enhancing resistance to salt stress. These findings suggest that salt-responsive AS directly and indirectly governs G. uralensis salt response. Further investigations into AS function and mechanism during abiotic stresses may offer novel references for bolstering plant stress tolerance.

CXIP4 depletion causes early lethality and pre-mRNA missplicing in Arabidopsis

Zinc knuckle (ZCCHC) motif-containing proteins are present in unicellular and multicellular eukaryotes and most ZCCHC proteins with known functions participate in the metabolism of various classes of RNA, such as mRNAs, ribosomal RNAs, and microRNAs. The Arabidopsis (Arabidopsis thaliana) genome encodes 69 ZCCHC-containing proteins, but the functions of most remain unclear. One of these proteins is CAX-INTERACTING PROTEIN 4 (CXIP4), which has been classified as a PTHR31437 family member, along with human SREK1-interacting protein 1 (SREK1IP1), which is thought to function in pre-mRNA splicing and RNA methylation. Metazoan SREK1IP1-like and plant CXIP4-like proteins only share a ZCCHC motif, and their functions remain almost entirely unknown. We studied two loss-of-function alleles of Arabidopsis CXIP4, the first mutations in PTHR31437 family genes described to date: cxip4-1 is likely null and shows early lethality, and cxip4-2 is hypomorphic and viable, with pleiotropic morphological defects. The cxip4-2 mutant exhibited deregulation of defense genes and upregulation of transcription factor encoding genes, some of which might explain its developmental defects. This mutant also exhibited increased intron retention events, and the specific functions of misspliced genes, such as those involved in "gene silencing by DNA methylation" and "mRNA polyadenylation factor" suggest that CXIP4 has additional functions. The CXIP4 protein localizes to the nucleus in a pattern resembling nuclear speckles, which are rich in splicing factors. Therefore, CXIP4 is required for plant survival and proper development, and mRNA maturation.

Comprehensive study of serine/arginine-rich (SR) gene family in rice: characterization, evolution and expression analysis

As important regulators of alternative splicing (AS) events, serine/arginine (SR)-rich proteins play indispensable roles in the growth and development of organisms. Until now, the study of SR genes has been lacking in plants. In the current study, we performed genome-wide analysis on the SR gene family in rice. A total of 24 OsSR genes were phylogenetically classified into seven groups, corresponding to seven subfamilies. The OsSR genes' structures, distribution of conserved domains, and protein tertiary structure of OsSR were conserved within each subfamily. The synteny analysis revealed that segmental duplication events were critical for the expansion of OsSR gene family. Moreover, interspecific synteny revealed the distribution of orthologous SR gene pairs between rice and Arabidopsis, sorghum, wheat, and maize. Among all OsSR genes, 14 genes exhibited NAGNAG acceptors, and only four OsSR genes had AS events on the NAGNAG acceptors. Furthermore, the distinct tissue-specific expression patterns of OsSR genes showed that these genes may function in different developmental stages in rice. The AS patterns on the same OsSR gene were variable among the root, stem, leaf, and grains at different filling stages, and some isoforms could only be detected in one or a few of tested tissues. Meanwhile, our results showed that the expression of some OsSR genes changed dramatically under ABA, GA, salt, drought, cold or heat treatment, which were related to the wide distribution of corresponding cis-elements in their promoter regions, suggesting their specific roles in stress and hormone response. This research facilitates our understanding of SR gene family in rice and provides clues for further exploration of the function of OsSR genes.

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

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

Research progress on the roles of lncRNAs in plant development and stress responses

Long non-coding RNAs (lncRNAs) are RNAs of more than 200 nucleotides in length that are not (or very rarely) translated into proteins. In eukaryotes, lncRNAs regulate gene expression at the transcriptional, post-transcriptional, and epigenetic levels. lncRNAs are categorized according to their genomic position and molecular mechanism. This review summarized the characteristics and mechanisms of plant lncRNAs involved in vegetative growth, reproduction, and stress responses. Our discussion and model provide a theoretical basis for further studies of lncRNAs in plant breeding.

m6A modification of U6 snRNA modulates usage of two major classes of pre-mRNA 5' splice site

Alternative splicing of messenger RNAs is associated with the evolution of developmentally complex eukaryotes. Splicing is mediated by the spliceosome, and docking of the pre-mRNA 5' splice site into the spliceosome active site depends upon pairing with the conserved ACAGA sequence of U6 snRNA. In some species, including humans, the central adenosine of the ACAGA box is modified by N6 methylation, but the role of this m6A modification is poorly understood. Here, we show that m6A modified U6 snRNA determines the accuracy and efficiency of splicing. We reveal that the conserved methyltransferase, FIONA1, is required for Arabidopsis U6 snRNA m6A modification. Arabidopsis fio1 mutants show disrupted patterns of splicing that can be explained by the sequence composition of 5' splice sites and cooperative roles for U5 and U6 snRNA in splice site selection. U6 snRNA m6A influences 3' splice site usage. We generalise these findings to reveal two major classes of 5' splice site in diverse eukaryotes, which display anti-correlated interaction potential with U5 snRNA loop 1 and the U6 snRNA ACAGA box. We conclude that U6 snRNA m6A modification contributes to the selection of degenerate 5' splice sites crucial to alternative splicing.

Missplicing suppressor alleles of Arabidopsis PRE-MRNA PROCESSING FACTOR 8 increase splicing fidelity by reducing the use of novel splice sites

Efficient splicing requires a balance between high-fidelity splice-site (SS) selection and speed. In Saccharomyces cerevisiae, Pre-mRNA processing factor 8 (Prp8) helps to balance precise SS selection and rapid, efficient intron excision and exon joining. argonaute1-52 (ago1-52) and incurvata13 (icu13) are hypomorphic alleles of the Arabidopsis thaliana genes ARGONAUTE1 (AGO1) and AUXIN RESISTANT6 (AXR6) that harbor point mutations creating a novel 3'SS and 5'SS, respectively. The spliceosome recognizes these novel SSs, as well as the intact genuine SSs, producing a mixture of wild-type and aberrant mature mRNAs. Here, we characterized five novel mutant alleles of PRP8 (one of the two Arabidopsis co-orthologs of yeast Prp8), naming these alleles morphology of ago1-52 suppressed5 (mas5). In the mas5-1 background, the spliceosome preferentially recognizes the intact genuine 3'SS of ago1-52 and 5'SS of icu13. Since point mutations that damage genuine SSs make the spliceosome prone to recognizing cryptic SSs, we also tested alleles of four genes carrying damaged genuine SSs, finding that mas5-1 did not suppress their missplicing. The mas5-1 and mas5-3 mutations represent a novel class of missplicing suppressors that increase splicing fidelity by hampering the use of novel SSs, but do not alter general pre-mRNA splicing.

Nucleotides in both donor and acceptor splice sites are responsible for choice in NAGNAG tandem splice sites

Among alternative splicing events in the human transcriptome, tandem NAGNAG acceptor splice sites represent an appreciable proportion. Both proximal and distal NAG can be used to produce two splicing isoforms differing by three nucleotides. In some cases, the upstream exon can be alternatively spliced as well, which further increases the number of possible transcripts. In this study, we showed that NAG choice in tandem splice site depends considerably not only on the concerned acceptor, but also on the upstream donor splice site sequence. Using an extensive set of experiments with systematically modified two-exonic minigene systems of AFAP1L2 or CSTD gene, we recognized the third and fifth intronic upstream donor splice site position and the tandem acceptor splice site region spanning from -10 to +2, including NAGNAG itself, as the main drivers. In addition, competition between different branch points and their composition were also shown to play a significant role in NAG choice. All these nucleotide effects appeared almost additive, which explained the high variability in proximal versus distal NAG usage.

BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq

BACKGROUND

The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants.

RESULTS

A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts - BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427-433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20-28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage.

CONCLUSION

A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.

BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq

BACKGROUND

The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants.

RESULTS

A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts - BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427-433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20-28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage.

CONCLUSION

A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.

Coordinated Regulation of Pre-mRNA Splicing by the SFPS-RRC1 Complex to Promote Photomorphogenesis

Light signals perceived by the phytochrome (phy) family of photoreceptors control gene expression at both transcriptional and posttranscriptional levels to promote photomorphogenesis. Recently, we identified a factor called SPLICING FACTOR FOR PHYTOCHROME SIGNALING (SFPS) that directly interacts with the photoreceptor phyB and regulates pre-mRNA splicing in Arabidopsis (Arabidopsis thaliana). To identify SFPS-interacting proteins, we performed an immunoprecipitation followed by a mass spectrometry and identified the Ser/Arg-like protein REDUCED RED-LIGHT RESPONSES IN CRY1CRY2 BACKGROUND1 (RRC1). Genetic analyses revealed that the sfps-2 rrc1-3 phenotypes are similar to those of the single mutants, suggesting that RRC1 and SFPS might function together. RNA sequence analyses of rrc1-3 identified a large number of genes whose pre-mRNA splicing is altered under dark and light conditions. Comparison of the sequence data revealed a subset of common genes coregulated by SFPS and RRC1 under dark and light conditions. Similar to SFPS, RRC1 also interacts with phyB, colocalizes in nuclear photobodies, and regulates light-dependent pre-mRNA splicing of a subset of genes. Taken together, these data suggest that although SFPS and RRC1 can regulate distinct subsets of genes, they also form a complex and coordinately control pre-mRNA splicing of a subset of genes involved in light signaling and circadian clock pathways to promote photomorphogenesis.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

An eukaryotic elongation factor 2 from Medicago falcata (MfEF2) confers cold tolerance

BACKGROUND

An eukaryotic translation elongation factor-2 (eEF-2) plays an important role in protein synthesis, however, investigation on its role in abiotic stress responses is limited. A cold responsive eEF2 named as MfEF2 was isolated from yellow-flowered alfalfa [Medicago sativa subsp. falcata (L.) Arcang, thereafter M. falcata], a forage legume with great cold tolerance, and transgenic tobacco (Nicotiana tabacum L.) plants overexpressing MfEF2 were analyzed in cold tolerance and proteomic profiling was conducted under low temperature in this study.

RESULTS

MfEF2 transcript was induced and peaked at 24 h and remained at the high level during cold treatment up to 96 h. Overexpression of MfEF2 in trasngenic tobacco plants resulted in enhanced cold tolerance. Compared to the wild type, transgenic plants showed higher survival rate after freezing treatment, higher levels of net photosynthetic rate (A), maximum photochemical efciency of photosystem (PS) II (Fv/Fm) and nonphotochemical quenching (NPQ) and lower levels of ion leakage and reactive oxygen species (ROS) production after chilling treatment. iTRAQ-based quantitative proteomic analysis identified 336 differentially expressed proteins (DEPs) from leaves of one transgenic line versus the wild type after chilling treatment for 48 h. GO and KEGG enrichment were conducted for analysis of the major biological process, cellular component, molecular function, and pathways of the DEPs involving in. It is interesting that many down-regulated DEPs were grouped into "photosynthesis" and "photosynthesis-antenna", such as subunits of PSI and PSII as well as light harvesting chlorophyll protein complex (LHC), while many up-regulated DEPs were grouped into "spliceosome".

CONCLUSIONS

The results suggest that MfEF2 confers cold tolerance through regulating hundreds of proteins synthesis under low temperature conditions. The elevated cold tolerance in MfEF2 transgenic plants was associated with downregulation of the subunits of PSI and PSII as well as LHC, which leads to reduced capacity for capturing sunlight and ROS production for protection of plants, and upregulation of proteins involving in splicesome, which promotes alternative splicing of pre-mRNA under low temperature.

Identification and analysis of alternative splicing events in Phaseolus vulgaris and Glycine max

BACKGROUND

The vast diversification of proteins in eukaryotic cells has been related with multiple transcript isoforms from a single gene that result in alternative splicing (AS) of primary transcripts. Analysis of RNA sequencing data from expressed sequence tags and next generation RNA sequencing has been crucial for AS identification and genome-wide AS studies. For the identification of AS events from the related legume species Phaseolus vulgaris and Glycine max, 157 and 88 publicly available RNA-seq libraries, respectively, were analyzed.

RESULTS

We identified 85,570 AS events from P. vulgaris in 72% of expressed genes and 134,316 AS events in 70% of expressed genes from G. max. These were categorized in seven AS event types with intron retention being the most abundant followed by alternative acceptor and alternative donor, representing ~75% of all AS events in both plants. Conservation of AS events in homologous genes between the two species was analyzed where an overrepresentation of AS affecting 5'UTR regions was observed for certain types of AS events. The conservation of AS events was experimentally validated for 8 selected genes, through RT-PCR analysis. The different types of AS events also varied by relative position in the genes. The results were consistent in both species.

CONCLUSIONS

The identification and analysis of AS events are first steps to understand their biological relevance. The results presented here from two related legume species reveal high conservation, over ~15-20 MY of divergence, and may point to the biological relevance of AS.

Functions of long intergenic non-coding (linc) RNAs in plants

Whole transcriptome analyses in many organisms have revealed that most transcribed RNAs do not encode proteins. These non-coding RNAs likely contribute to the regulation of gene expression during the development of multicellular organisms. In eukaryotes, the roles of small RNAs, one class of non-coding RNAs, in transcriptional and post-transcriptional regulation have been well characterized. However, the functions of a second class of non-coding RNAs, long intergenic noncoding (linc) RNAs, are relatively unknown, especially in plants. Recent advances in RNA-seq and tiling microarray technologies have revealed the presence of many lincRNAs across plant species. This review focuses on the functions of lincRNAs that have been recently reported in plants. One of the most well characterized functions of lincRNAs is to epigenetically regulate gene expression by recruiting proteins for chromosome modification to specific loci. Second, lincRNAs are known to inhibit the physical interaction between microRNAs (miRNAs) and their target mRNAs thus controling protein levels of the target mRNAs. Lastly, lincRNAs control alternative splicing by binding and sequestering the proteins required for alternative splicing.

CoSpliceNet: a framework for co-splicing network inference from transcriptomics data

Background

Alternative splicing has been proposed to increase transcript diversity and protein plasticity in eukaryotic organisms, but the extent to which this is the case is currently unclear, especially with regard to the diversification of molecular function. Eukaryotic splicing involves complex interactions of splicing factors and their targets. Inference of co-splicing networks capturing these types of interactions is important for understanding this crucial, highly regulated post-transcriptional process at the systems level.

Results

First, several transcript and protein attributes, including coding potential of transcripts and differences in functional domains of proteins, were compared between splice variants and protein isoforms to assess transcript and protein diversity in a biological system. Alternative splicing was shown to increase transcript and function-related protein diversity in developing Arabidopsis embryos. Second, CoSpliceNet, which integrates co-expression and motif discovery at splicing regulatory regions to infer co-splicing networks, was developed. CoSpliceNet was applied to temporal RNA sequencing data to identify candidate regulators of splicing events and predict RNA-binding motifs, some of which are supported by prior experimental evidence. Analysis of inferred splicing factor targets revealed an unexpected role for the unfolded protein response in embryo development.

Conclusions

The methods presented here can be used in any biological system to assess transcript diversity and protein plasticity and to predict candidate regulators, their targets, and RNA-binding motifs for splicing factors. CoSpliceNet is freely available at http://delasa.github.io/co-spliceNet/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3172-6) contains supplementary material, which is available to authorized users.

Merging Absolute and Relative Quantitative PCR Data to Quantify STAT3 Splice Variant Transcripts

Human signal transducer and activator of transcription 3 (STAT3) is one of many genes containing a tandem splicing site. Alternative donor splice sites 3 nucleotides apart result in either the inclusion (S) or exclusion (ΔS) of a single residue, Serine-701. Further downstream, splicing at a pair of alternative acceptor splice sites result in transcripts encoding either the 55 terminal residues of the transactivation domain (α) or a truncated transactivation domain with 7 unique residues (β). As outlined in this manuscript, measuring the proportions of STAT3's four spliced transcripts (Sα, Sβ, ΔSα and ΔSβ) was possible using absolute qPCR (quantitative polymerase chain reaction). The protocol therefore distinguishes and measures highly similar splice variants. Absolute qPCR makes use of calibrator plasmids and thus specificity of detection is not compromised for the sake of efficiency. The protocol necessitates primer validation and optimization of cycling parameters. A combination of absolute qPCR and efficiency-dependent relative qPCR of total STAT3 transcripts allowed a description of the fluctuations of STAT3 splice variants' levels in eosinophils treated with cytokines. The protocol also provided evidence of a co-splicing interdependence between the two STAT3 splicing events. The strategy based on a combination of the two qPCR techniques should be readily adaptable to investigation of co-splicing at other tandem splicing sites.

It's a bit over, is that ok? The subtle surplus from tandem alternative splicing

Tandem alternative splice sites (TASS) form a defined class of alternative splicing and give rise to mRNA insertion/deletion variants with only small size differences. Previous work has confirmed evolutionary conservation of TASS elements while many cases show only low tissue specificity of isoform ratios. We pinpoint stochasticity and noise as important methodological issues for the dissection of TASS isoform patterns. Resolving such uncertainties, a recent report showed regulation in a cell culture system, with shifts of alternative splicing isoform ratios dependent on cell density. This novel type of regulation affects not only multiple TASS isoforms, but also other alternative splicing classes, in a concerted manner. Here, we discuss how specific regulatory network architectures may be realized through the novel regulation type and highlight the role of differential isoform functions as a key step in order to better understand the functional role of TASS.

nagnag: Identification and quantification of NAGNAG alternative splicing using RNA-Seq data

Regulation of proteome diversity by alternative splicing has been widely demonstrated in plants and animals. NAGNAG splicing, which was recently defined as a tissue specific event, results in the production of two distinct isoforms that are distinguished by three nucleotides (NAG) as a consequence of the intron proximal or distal to the splice site. Since the NAGNAG mechanism is not well characterized, tools for the identification and quantification of NAGNAG splicing events remain under-developed. Here we report nagnag, an R-based NAGNAG splicing detection tool, which accurately identifies and quantifies NAGNAG splicing events using RNA-Seq. Overall, nagnag produces user-friendly visualization reports and highlights differences between the DNA/RNA/protein across the identified isoforms of the reported gene. The package is available on https://sourceforge.net/projects/nagnag/files/; or http://genome.sdau.edu.cn/research/software/nagnag.html.

Analysis of pollen-specific alternative splicing in Arabidopsis thaliana via semi-quantitative PCR

Alternative splicing enables a single gene to produce multiple mRNA isoforms by varying splice site selection. In animals, alternative splicing of mRNA isoforms between cell types is widespread and supports cellular differentiation. In plants, at least 20% of multi-exon genes are alternatively spliced, but the extent and significance of tissue-specific splicing is less well understood, partly because it is difficult to isolate cells of a single type. Pollen is a useful model system to study tissue-specific splicing in higher plants because pollen grains contain only two cell types and can be collected in large amounts without damaging cells. Previously, we identified pollen-specific splicing patterns by comparing RNA-Seq data from Arabidopsis pollen and leaves. Here, we used semi-quantitative PCR to validate pollen-specific splicing patterns among genes where RNA-Seq data analysis indicated splicing was most different between pollen and leaves. PCR testing confirmed eight of nine alternative splicing patterns, and results from the ninth were inconclusive. In four genes, alternative transcriptional start sites coincided with alternative splicing. This study highlights the value of the low-cost PCR assay as a method of validating RNA-Seq results.

Long noncoding RNA modulates alternative splicing regulators in Arabidopsis

Alternative splicing (AS) of pre-mRNA represents a major mechanism underlying increased transcriptome and proteome complexity. Here, we show that the nuclear speckle RNA-binding protein (NSR) and the AS competitor long noncoding RNA (or ASCO-lncRNA) constitute an AS regulatory module. AtNSR-GFP translational fusions are expressed in primary and lateral root (LR) meristems. Double Atnsr mutants and ASCO overexpressors exhibit an altered ability to form LRs after auxin treatment. Interestingly, auxin induces a major change in AS patterns of many genes, a response largely dependent on NSRs. RNA immunoprecipitation assays demonstrate that AtNSRs interact not only with their alternatively spliced mRNA targets but also with the ASCO-RNA in vivo. The ASCO-RNA displaces an AS target from an NSR-containing complex in vitro. Expression of ASCO-RNA in Arabidopsis affects the splicing patterns of several NSR-regulated mRNA targets. Hence, lncRNA can hijack nuclear AS regulators to modulate AS patterns during development.

Physiological state co-regulates thousands of mammalian mRNA splicing events at tandem splice sites and alternative exons

Thousands of tandem alternative splice sites (TASS) give rise to mRNA insertion/deletion variants with small size differences. Recent work has concentrated on the question of biological relevance in general, and the physiological regulation of TASS in particular. We have quantitatively studied 11 representative TASS cases in comparison to one mutually exclusive exon case and two cassette exons (CEs) using a panel of human and mouse tissues, as well as cultured cell lines. Tissues show small but significant differences in TASS isoform ratios, with a variance 4- to 20-fold lower than seen for CEs. Remarkably, in cultured cells, all studied alternative splicing (AS) cases showed a cell-density-dependent shift of isoform ratios with similar time series profiles. A respective genome-wide co-regulation of TASS splicing was shown by next-generation mRNA sequencing data. Moreover, data from human and mouse organs indicate that this co-regulation of TASS occurs in vivo, with brain showing the strongest difference to other organs. Together, the results indicate a physiological AS regulation mechanism that functions almost independently from the splice site context and sequence.

A subtle alternative splicing event gives rise to a widely expressed human RNase k isoform

Subtle alternative splicing leads to the formation of RNA variants lacking or including a small number of nucleotides. To date, the impact of subtle alternative splicing phenomena on protein biosynthesis has been studied in frame-preserving incidents. On the contrary, mRNA isoforms derived from frame-shifting events were poorly studied and generally characterized as non-coding. This work provides evidence for a frame-shifting subtle alternative splicing event which results in the production of a novel protein isoform. We applied a combined molecular approach for the cloning and expression analysis of a human RNase κ transcript (RNase κ-02) which lacks four consecutive bases compared to the previously isolated RNase κ isoform. RNase κ-02 mRNA is expressed in all human cell lines tested end encodes the synthesis of a 134-amino-acid protein by utilizing an alternative initiation codon. The expression of RNase κ-02 in the cytoplasm of human cells was verified by Western blot and immunofluorescence analysis using a specific polyclonal antibody developed on the basis of the amino-acid sequence difference between the two protein isoforms. The results presented here show that subtle changes during mRNA splicing can lead to the expression of significantly altered protein isoforms.

Genome-wide study of NAGNAG alternative splicing in Arabidopsis

NAGNAG alternative splicing is one type of alternative splicing in mammals and plants. There are two opposite arguments regarding the mechanism of this NAGNAG event, i.e. whether splice variation is controllable by the cell or is just biological noise. In this paper, we systematically investigated NAGNAG acceptors in Arabidopsis thaliana using both cDNA/EST and RNA-Seq data. We identified 9,473 NAGNAG motifs, including 529 cDNA/EST-confirmed NAGNAG acceptors. A nomenclature tree for this type of alternative splicing was defined based on the cDNA/EST validation, location in the exon, sequence and expression level. Low expression of some NAGNAG motifs was observed in various tissues or pathogen-infected samples, indicating the existence of background splicing. Tissue-specific or treatment-specific differences in the dynamic profiles suggest that some NAGNAG acceptors are highly regulated.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

Alternative splicing of transcription factors in plant responses to low temperature stress: mechanisms and functions

Transcription factors play a central role in the gene regulatory networks that mediate various aspects of plant developmental processes and responses to environmental changes. Therefore, their activities are elaborately regulated at multiple steps. In particular, accumulating evidence illustrates that post-transcriptional control of mRNA metabolism is a key molecular scheme that modulates the transcription factor activities in plant responses to temperature fluctuations. Transcription factors have a modular structure consisting of distinct protein domains essential for DNA binding, dimerization, and transcriptional regulation. Alternative splicing produces multiple proteins having different structural domain compositions from a single transcription factor gene. Recent studies have shown that alternative splicing of some transcription factor genes generates small interfering peptides (siPEPs) that negatively regulate the target transcription factors via peptide interference (PEPi), constituting self-regulatory circuits in plant cold stress response. A number of splicing factors, which are involved in RNA binding, splice site selection, and spliceosome assembly, are also affected by temperature fluctuations, supporting the close association of alternative splicing of transcription factors with plant responses to low temperatures. In this review, we summarize recent progress on the temperature-responsive alternative splicing of transcription factors in plants with emphasis on the siPEP-mediated PEPi mechanism.

Deciphering the plant splicing code: experimental and computational approaches for predicting alternative splicing and splicing regulatory elements

Extensive alternative splicing (AS) of precursor mRNAs (pre-mRNAs) in multicellular eukaryotes increases the protein-coding capacity of a genome and allows novel ways to regulate gene expression. In flowering plants, up to 48% of intron-containing genes exhibit AS. However, the full extent of AS in plants is not yet known, as only a few high-throughput RNA-Seq studies have been performed. As the cost of obtaining RNA-Seq reads continues to fall, it is anticipated that huge amounts of plant sequence data will accumulate and help in obtaining a more complete picture of AS in plants. Although it is not an onerous task to obtain hundreds of millions of reads using high-throughput sequencing technologies, computational tools to accurately predict and visualize AS are still being developed and refined. This review will discuss the tools to predict and visualize transcriptome-wide AS in plants using short-reads and highlight their limitations. Comparative studies of AS events between plants and animals have revealed that there are major differences in the most prevalent types of AS events, suggesting that plants and animals differ in the way they recognize exons and introns. Extensive studies have been performed in animals to identify cis-elements involved in regulating AS, especially in exon skipping. However, few such studies have been carried out in plants. Here, we review the current state of research on splicing regulatory elements (SREs) and briefly discuss emerging experimental and computational tools to identify cis-elements involved in regulation of AS in plants. The availability of curated alternative splice forms in plants makes it possible to use computational tools to predict SREs involved in AS regulation, which can then be verified experimentally. Such studies will permit identification of plant-specific features involved in AS regulation and contribute to deciphering the splicing code in plants.

Mining Arabidopsis thaliana RNA-seq data with Integrated Genome Browser reveals stress-induced alternative splicing of the putative splicing regulator SR45a

PREMISE OF THE STUDY

High-throughput sequencing of cDNA libraries prepared from diverse samples (RNA-seq) can reveal genome-wide changes in alternative splicing. Using RNA-seq data to assess splicing at the level of individual genes requires the ability to visualize read alignments alongside genomic annotations. To meet this need, we added RNA-seq visualization capability to Integrated Genome Browser (IGB), a free desktop genome visualization tool. To illustrate this capability, we present an in-depth analysis of abiotic stresses and their effects on alternative splicing of SR45a (AT1G07350), a putative splicing regulator from Arabidopsis thaliana.

METHODS

cDNA libraries prepared from Arabidopsis plants that were subjected to heat and dehydration stresses were sequenced on an Illumina GAIIx sequencer, yielding more than 511 million high-quality 75-base, single-end sequence reads. Reads were aligned onto the reference genome and visualized in IGB.

KEY RESULTS

Using IGB, we confirmed exon-skipping alternative splicing in SR45a. Exon-skipped variant AT1G07350.1 encodes full-length SR45a protein with intact RS and RNA recognition motifs, while nonskipped variant AT1G07350.2 lacks the C-terminal RS region due to a frameshift in the alternative exon. Heat and drought stresses increased both transcript abundance and the proportion of exon-skipped transcripts encoding the full-length protein. We identified new splice sites and observed frequent intron retention flanking the alternative exon.

CONCLUSIONS

This study underlines the importance of visual inspection of RNA-seq alignments when investigating alternatively spliced genes. We showed that heat and dehydration stresses increase overall abundance of SR45a mRNA while also increasing production of transcripts encoding the full-length SR45a protein relative to other splice variants.

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ∼15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.

Plant serine/arginine-rich proteins: roles in precursor messenger RNA splicing, plant development, and stress responses

Global analyses of splicing of precursor messenger RNAs (pre-mRNAs) have revealed that alternative splicing (AS) is highly pervasive in plants. Despite the widespread occurrence of AS in plants, the mechanisms that control splicing and the roles of splice variants generated from a gene are poorly understood. Studies on plant serine/arginine-rich (SR) proteins, a family of highly conserved proteins, suggest their role in both constitutive splicing and AS of pre-mRNAs. SR proteins have a characteristic domain structure consisting of one or two RNA recognition motifs at the N-terminus and a C-terminal RS domain rich in arginine/serine dipeptides. Plants have many more SR proteins compared to animals including several plant-specific subfamilies. Pre-mRNAs of plant SR proteins are extensively alternatively spliced to increase the transcript complexity by about six-fold. Some of this AS is controlled in a tissue- and development-specific manner. Furthermore, AS of SR pre-mRNAs is altered by various stresses, raising the possibility of rapid reprogramming of the whole transcriptome by external signals through regulation of the splicing of these master regulators of splicing. Most SR splice variants contain a premature termination codon and are degraded by up-frameshift 3 (UPF3)-mediated nonsense-mediated decay (NMD), suggesting a link between NMD and regulation of expression of the functional transcripts of SR proteins. Limited functional studies with plant SRs suggest key roles in growth and development and plant responses to the environment. Here, we discuss the current status of research on plant SRs and some promising approaches to address many unanswered questions about plant SRs.

Strict 3' splice site sequence requirements for U2 snRNP recruitment after U2AF binding underlie a genetic defect leading to autoimmune disease

We report that the 3' splice site associated with the alternatively spliced exon 6 of the Fas receptor CD95 displays strict sequence requirements and that a mutation that disrupts this particular sequence arrangement leads to constitutive exon 6 skipping in a patient suffering from autoimmune lymphoproliferative syndrome (ALPS). Specifically, we find an absolute requirement for RCAG/G at the 3' splice site (where R represents purine, and / indicates the intron/exon boundary) and the balance between exon inclusion and skipping is exquisitely sensitive to single nucleotide variations in the uridine content of the upstream polypyrimidine (Py)-tract. Biochemical experiments revealed that the ALPS patient mutation reduces U2 snRNP recruitment to the 3' splice site region and that this effect cannot be explained by decreased interaction with the U2 snRNP Auxiliary Factor U2AF, whose 65- and 35-kDa subunits recognize the Py-tract and 3' splice site AG, respectively. The effect of the mutation, which generates a tandem of two consecutive AG dinucleotides at the 3' splice site, can be suppressed by increasing the distance between the AGs, mutating the natural 3' splice site AG or increasing the uridine content of the Py-tract at a position distal from the 3' splice site. The suppressive effects of these additional mutations correlate with increased recruitment of U2 snRNP but not with U2AF binding, again suggesting that the strict architecture of Fas intron 5 3' splice site region is tuned to regulate alternative exon inclusion through modulation of U2 snRNP assembly after U2AF binding.

Prevalence of alternative splicing choices in Arabidopsis thaliana

BACKGROUND

Around 14% of protein-coding genes of Arabidopsis thaliana genes from the TAIR9 genome release are annotated as producing multiple transcript variants through alternative splicing. However, for most alternatively spliced genes in Arabidopsis, the relative expression level of individual splicing variants is unknown.

RESULTS

We investigated prevalence of alternative splicing (AS) events in Arabidopsis thaliana using ESTs. We found that for most AS events with ample EST coverage, the majority of overlapping ESTs strongly supported one major splicing choice, with less than 10% of ESTs supporting the minor form. Analysis of ESTs also revealed a small but noteworthy subset of genes for which alternative choices appeared with about equal prevalence, suggesting that for these genes the variant splicing forms co-occur in the same cell types. Of the AS events in which both forms were about equally prevalent, more than 80% affected untranslated regions or involved small changes to the encoded protein sequence.

CONCLUSIONS

Currently available evidence from ESTs indicates that alternative splicing in Arabidopsis occurs and affects many genes, but for most genes with documented alternative splicing, one AS choice predominates. To aid investigation of the role AS may play in modulating function of Arabidopsis genes, we provide an on-line resource (ArabiTag) that supports searching AS events by gene, by EST library keyword search, and by relative prevalence of minor and major forms.

Sequence features involved in the mechanism of 3' splice junction wobbling

BACKGROUND

Alternative splicing is an important mechanism mediating the diversified functions of genes in multicellular organisms, and such event occurs in around 40-60% of human genes. Recently, a new splice-junction wobbling mechanism was proposed that subtle modifications exist in mRNA maturation by alternatively choosing at 5'- GTNGT and 3'- NAGNAG, which created single amino acid insertion and deletion isoforms.

RESULTS

By browsing the Alternative Splicing Database information, we observed that most 3' alternative splice site choices occur within six nucleotides of the dominant splice site and the incidence significantly decreases further away from the dominant acceptor site. Although a lower frequency of alternative splicing occurs within the intronic region (alternative splicing at the proximal AG) than in the exonic region (alternative splicing at the distal AG), alternative AG sites located within the intronic region show stronger potential as the acceptor. These observations revealed that the choice of 3' splice sites during 3' splicing junction wobbling could depend on the distance between the duplicated AG and the branch point site (BPS). Further mutagenesis experiments demonstrated that the distance of AG-to-AG and BPS-to-AG can greatly influence 3' splice site selection. Knocking down a known alternative splicing regulator, hSlu7, failed to affect wobble splicing choices.

CONCLUSION

Our results implied that nucleotide distance between proximal and distal AG sites has an important regulatory function. In this study, we showed that occurrence of 3' wobble splicing occurs in a distance-dependent manner and that most of this wobble splicing is probably caused by steric hindrance from a factor bound at the neighboring tandem motif sequence.

Identification and characterization of NAGNAG alternative splicing in the moss Physcomitrella patens

BACKGROUND

Alternative splicing (AS) involving tandem acceptors that are separated by three nucleotides (NAGNAG) is an evolutionarily widespread class of AS, which is well studied in Homo sapiens (human) and Mus musculus (mouse). It has also been shown to be common in the model seed plants Arabidopsis thaliana and Oryza sativa (rice). In one of the first studies involving sequence-based prediction of AS in plants, we performed a genome-wide identification and characterization of NAGNAG AS in the model plant Physcomitrella patens, a moss.

RESULTS

Using Sanger data, we found 295 alternatively used NAGNAG acceptors in P. patens. Using 31 features and training and test datasets of constitutive and alternative NAGNAGs, we trained a classifier to predict the splicing outcome at NAGNAG tandem splice sites (alternative splicing, constitutive at the first acceptor, or constitutive at the second acceptor). Our classifier achieved a balanced specificity and sensitivity of >or= 89%. Subsequently, a classifier trained exclusively on data well supported by transcript evidence was used to make genome-wide predictions of NAGNAG splicing outcomes. By generation of more transcript evidence from a next-generation sequencing platform (Roche 454), we found additional evidence for NAGNAG AS, with altogether 664 alternative NAGNAGs being detected in P. patens using all currently available transcript evidence. The 454 data also enabled us to validate the predictions of the classifier, with 64% (80/125) of the well-supported cases of AS being predicted correctly.

CONCLUSION

NAGNAG AS is just as common in the moss P. patens as it is in the seed plants A. thaliana and O. sativa (but not conserved on the level of orthologous introns), and can be predicted with high accuracy. The most informative features are the nucleotides in the NAGNAG and in its immediate vicinity, along with the splice sites scores, as found earlier for NAGNAG AS in animals. Our results suggest that the mechanism behind NAGNAG AS in plants is similar to that in animals and is largely dependent on the splice site and its immediate neighborhood.

Comparison of methods for quantification of subtle splice variants

Alternative splicing is capable of generating multiple mRNA variants from a single gene and is hence a key mediator of molecular diversity generated at the transcript level. Consequently, delivering quantitative information on the fractions of splice variants is essential for the understanding of their biological roles. Here we compare techniques for subtle splice variant quantification that are able to resolve length differences as small as one nucleotide: PAGE with ethidium-bromide densitometry, pyrosequencing, and CE-LIF. We give comprehensive descriptions of assay designs and calibration procedures and present an evaluation of these methods in terms of accuracy, reproducibility and applicability. We also examined template concentrations and reverse transcription-coupled PCR conditions as potential cause of biased results as they were observed for extreme low template concentrations and/or PCR amplicons with size differences of 195 nt. As proof of concept, we determine the splice ratios of variants differing by 3 and 12 nt in five human tissues. We demonstrate that CE-LIF is the most precise and also the most labor- and time-efficient method.

Comparative analysis indicates that alternative splicing in plants has a limited role in functional expansion of the proteome

Background

Alternative splicing (AS) is a widespread phenomenon in higher eukaryotes but the extent to which it leads to functional protein isoforms and to proteome expansion at large is still a matter of debate. In contrast to animal species, for which AS has been studied extensively at the protein and functional level, protein-centered studies of AS in plant species are scarce. Here we investigate the functional impact of AS in dicot and monocot plant species using a comparative approach.

Results

Detailed comparison of AS events in alternative spliced orthologs from the dicot Arabidopsis thaliana and the monocot Oryza sativa (rice) revealed that the vast majority of AS events in both species do not result from functional conservation. Transcript isoforms that are putative targets for the nonsense-mediated decay (NMD) pathway are as likely to contain conserved AS events as isoforms that are translated into proteins. Similar results were obtained when the same comparison was performed between the two more closely related monocot species rice and Zea mays (maize).

Genome-wide computational analysis of functional protein domains encoded in alternatively and constitutively spliced genes revealed that only the RNA recognition motif (RRM) is overrepresented in alternatively spliced genes in all species analyzed. In contrast, three domain types were overrepresented in constitutively spliced genes. AS events were found to be less frequent within than outside predicted protein domains and no domain type was found to be enriched with AS introns. Analysis of AS events that result in the removal of complete protein domains revealed that only a small number of domain types is spliced-out in all species analyzed. Finally, in a substantial fraction of cases where a domain is completely removed, this domain appeared to be a unit of a tandem repeat.

Conclusion

The results from the ortholog comparisons suggest that the ability of a gene to produce more than one functional protein through AS does not persist during evolution. Cross-species comparison of the results of the protein-domain oriented analyses indicates little correspondence between the analyzed species. Based on the premise that functional genetic features are most likely to be conserved during evolution, we conclude that AS has only a limited role in functional expansion of the proteome in plants.

Accurate prediction of NAGNAG alternative splicing

Alternative splicing (AS) involving NAGNAG tandem acceptors is an evolutionarily widespread class of AS. Recent predictions of alternative acceptor usage reported better results for acceptors separated by larger distances, than for NAGNAGs. To improve the latter, we aimed at the use of Bayesian networks (BN), and extensive experimental validation of the predictions. Using carefully constructed training and test datasets, a balanced sensitivity and specificity of ≥92% was achieved. A BN trained on the combined dataset was then used to make predictions, and 81% (38/47) of the experimentally tested predictions were verified. Using a BN learned on human data on six other genomes, we show that while the performance for the vertebrate genomes matches that achieved on human data, there is a slight drop for Drosophila and worm. Lastly, using the prediction accuracy according to experimental validation, we estimate the number of yet undiscovered alternative NAGNAGs. State of the art classifiers can produce highly accurate prediction of AS at NAGNAGs, indicating that we have identified the major features of the ‘NAGNAG-splicing code’ within the splice site and its immediate neighborhood. Our results suggest that the mechanism behind NAGNAG AS is simple, stochastic, and conserved among vertebrates and beyond.