Wobble splicing reveals the role of the branch point sequence-to-NAGNAG region in 3' tandem splice site selection

Considerations for Generating Humanized Mouse Models to Test Efficacy of Antisense Oligonucleotides

Over the last decades, animal models have become increasingly important in therapeutic drug development and assessment. The use of these models, mainly mice and rats, allow evaluating drugs in the real-organism environment and context. However, several molecular therapeutic approaches are sequence-dependent, and therefore, the humanization of such models is required to assess the efficacy. The generation of genetically modified humanized mouse models is often an expensive and laborious process that may not always recapitulate the human molecular and/or physiological phenotype. In this chapter, we summarize basic aspects to consider before designing and generating humanized models, especially when they are aimed to test antisense-based therapies.

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.

Top-ranked expressed gene transcripts of human protein-coding genes investigated with GTEx dataset

With considerable accumulation of RNA-Seq transcriptome data, we have extended our understanding about protein-coding gene transcript compositions. However, alternatively compounded patterns of human protein-coding gene transcripts would complicate gene expression data processing and interpretation. It is essential to exhaustively interrogate complex mRNA isoforms of protein-coding genes with an unified data resource. In order to investigate representative mRNA transcript isoforms to be utilized as transcriptome analysis references, we utilized GTEx data to establish a top-ranked transcript isoform expression data resource for human protein-coding genes. Distinctive tissue specific expression profiles and modulations could be observed for individual top-ranked transcripts of protein-coding genes. Protein-coding transcripts or genes do occupy much higher expression fraction in transcriptome data. In addition, top-ranked transcripts are the dominantly expressed ones in various normal tissues. Intriguingly, some of the top-ranked transcripts are noncoding splicing isoforms, which imply diverse gene regulation mechanisms. Comprehensive investigation on the tissue expression patterns of top-ranked transcript isoforms is crucial. Thus, we established a web tool to examine top-ranked transcript isoforms in various human normal tissue types, which provides concise transcript information and easy-to-use graphical user interfaces. Investigation of top-ranked transcript isoforms would contribute understanding on the functional significance of distinctive alternatively spliced transcript isoforms.

Nuclear Accumulation of LAP1:TRF2 Complex during DNA Damage Response Uncovers a Novel Role for LAP1

Lamina-associated polypeptide 1 (LAP1) is a nuclear envelope (NE) protein whose function remains poorly characterized. In a recent LAP1 protein interactome study, a putative regulatory role in the DNA damage response (DDR) has emerged and telomeric repeat-binding factor 2 (TRF2), a protein intimately associated with this signaling pathway, was among the list of LAP1 interactors. To gain insights into LAP1's physiological properties, the interaction with TRF2 in human cells exposed to DNA-damaging agents was investigated. The direct LAP1:TRF2 binding was validated in vitro by blot overlay and in vivo by co-immunoprecipitation after hydrogen peroxide and bleomycin treatments. The regulation of this protein interaction by LAP1 phosphorylation was demonstrated by co-immunoprecipitation and mass spectrometry following okadaic acid exposure. The involvement of LAP1 and TRF2 in the DDR was confirmed by their increased nuclear protein levels after bleomycin treatment, evaluated by immunoblotting, as well as by their co-localization with DDR factors at the NE and within the nucleoplasm, assessed by immunocytochemistry. Effectively, we showed that the LAP1:TRF2 complex is established during a cellular response against DNA damage. This work proposes a novel functional role for LAP1 in the DDR, revealing a potential biological mechanism that may be disrupted in LAP1-associated pathologies.

Cloning and functional complementation of ten Schistosoma mansoni phosphodiesterases expressed in the mammalian host stages

Only a single drug against schistosomiasis is currently available and new drug development is urgently required but very few drug targets have been validated and characterised. However, regulatory systems including cyclic nucleotide metabolism are emerging as primary candidates for drug discovery. Here, we report the cloning of ten cyclic nucleotide phosphodiesterase (PDE) genes of S. mansoni, out of a total of 11 identified in its genome. We classify these PDEs by homology to human PDEs. Male worms displayed higher expression levels for all PDEs, in mature and juvenile worms, and schistosomula. Several functional complementation approaches were used to characterise these genes. We constructed a Trypanosoma brucei cell line in which expression of a cAMP-degrading PDE complements the deletion of TbrPDEB1/B2. Inhibitor screens of these cells expressing only either SmPDE4A, TbrPDEB1 or TbrPDEB2, identified highly potent inhibitors of the S. mansoni enzyme that elevated the cellular cAMP concentration. We further expressed most of the cloned SmPDEs in two pde1Δ/pde2Δ strains of Saccharomyces cerevisiae and some also in a specialised strain of Schizosacharomyces pombe. Five PDEs, SmPDE1, SmPDE4A, SmPDE8, SmPDE9A and SmPDE11 successfully complemented the S. cerevisiae strains, and SmPDE7var also complemented to a lesser degree, in liquid culture. SmPDE4A, SmPDE8 and SmPDE11 were further assessed in S. pombe for hydrolysis of cAMP and cGMP; SmPDE11 displayed considerable preferrence for cGMP over cAMP. These results and tools enable the pursuit of a rigorous drug discovery program based on inhibitors of S. mansoni PDEs.

Impact of acceptor splice site NAGTAG motif on exon recognition

Pre-mRNA splicing is an essential step in gene expression, when introns are removed and exons joined by the complex of proteins called spliceosome. Correct splicing requires a precise exon/intron junction definition, which is determined by a consensual donor and acceptor splice site at the 5' and 3' end, respectively. An acceptor splice site (3'ss) consists of highly conserved AG nucleotides in positions E^-2 and E^-1. These nucleotides can appear in tandem, located 3 bp from each other. Then they are referred to as NAGNAG or tandem 3'ss, which can be alternatively spliced. NAG/TAG 3'ss motif abundance is extremely low and cannot be easily explained by just a nucleotide preference in this position. We tested artificial NAG/TAG motif's potential negative effect on exon recognition using a minigene assay. Introducing the NAG/TAG motif into seven different exons revealed no general negative effect on exon recognition. The only observed effect was the partial use of the newly formed distal 3'ss. We can conclude that this motif's extremely low preference in a natural 3'ss is not a consequence of the NAG/TAG motif's negative effect on exon recognition, but more likely the result of other RNA processing aspects, such as an alternative 3'ss choice, decreased 3'ss strength, or incorporating an amber stop codon.

Protein interaction network of alternatively spliced NudCD1 isoforms

NudCD1, also known as CML66 or OVA66, is a protein initially identified as overexpressed in patients with chronic myelogenous leukemia. The mRNA of NudCD1 is expressed in heart and testis of normal tissues, and is overexpressed in several cancers. Previous studies have shown that the expression level of the protein correlates with tumoral phenotype, possibly interacting upstream of the Insulin Growth Factor - 1 Receptor (IGF-1R). The gene encoding the NudCD1 protein consists of 12 exons that can be alternative spliced, leading to the expression of three different isoforms. These isoforms possess a common region of 492 amino acids in their C-terminus region and have an isoform specific N-terminus. To determine the distinct function of each isoforms, we have localised the isoforms within the cells using immunofluorescence microscopy and used a quantitative proteomics approach (SILAC) to identify specific protein interaction partners for each isoforms. Localization studies showed a different subcellular distribution for the different isoforms, with the first isoform being nuclear, while the other two isoforms have distinct cytoplasmic and nuclear location. We found that the different NudCD1 isoforms have unique interacting partners, with the first isoform binding to a putative RNA helicase named DHX15 involved in mRNA splicing.

Coordinated tissue-specific regulation of adjacent alternative 3' splice sites in C. elegans

Adjacent alternative 3′ splice sites, those separated by ≤18 nucleotides, provide a unique problem in the study of alternative splicing regulation; there is overlap of the cis-elements that define the adjacent sites. Identification of the intron's 3′ end depends upon sequence elements that define the branchpoint, polypyrimidine tract, and terminal AG dinucleotide. Starting with RNA-seq data from germline-enriched and somatic cell-enriched Caenorhabditis elegans samples, we identify hundreds of introns with adjacent alternative 3′ splice sites. We identify 203 events that undergo tissue-specific alternative splicing. For these, the regulation is monodirectional, with somatic cells preferring to splice at the distal 3′ splice site (furthest from the 5′ end of the intron) and germline cells showing a distinct shift toward usage of the adjacent proximal 3′ splice site (closer to the 5′ end of the intron). Splicing patterns in somatic cells follow C. elegans consensus rules of 3′ splice site definition; a short stretch of pyrimidines preceding an AG dinucleotide. Splicing in germline cells occurs at proximal 3′ splice sites that lack a preceding polypyrimidine tract, and in three instances the germline-specific site lacks the AG dinucleotide. We provide evidence that use of germline-specific proximal 3′ splice sites is conserved across Caenorhabditis species. We propose that there are differences between germline and somatic cells in the way that the basal splicing machinery functions to determine the intron terminus.

Identification of a novel human LAP1 isoform that is regulated by protein phosphorylation

Lamina associated polypeptide 1 (LAP1) is an integral protein of the inner nuclear membrane that is ubiquitously expressed. LAP1 binds to lamins and chromatin, probably contributing to the maintenance of the nuclear envelope architecture. Moreover, LAP1 also interacts with torsinA and emerin, proteins involved in DYT1 dystonia and X-linked Emery-Dreifuss muscular dystrophy disorder, respectively. Given its relevance to human pathological conditions, it is important to better understand the functional diversity of LAP1 proteins. In rat, the LAP1 gene (TOR1AIP1) undergoes alternative splicing to originate three LAP1 isoforms (LAP1A, B and C). However, it remains unclear if the same occurs with the human TOR1AIP1 gene, since only the LAP1B isoform had thus far been identified in human cells. In silico analysis suggested that, across different species, potential new LAP1 isoforms could be generated by alternative splicing. Using shRNA to induce LAP1 knockdown and HPLC-mass spectrometry analysis the presence of two isoforms in human cells was described and validated: LAP1B and LAP1C; the latter is putatively N-terminal truncated. LAP1B and LAP1C expression profiles appear to be dependent on the specific tissues analyzed and in cultured cells LAP1C was the major isoform detected. Moreover, LAP1B and LAP1C expression increased during neuronal maturation, suggesting that LAP1 is relevant in this process. Both isoforms were found to be post-translationally modified by phosphorylation and methionine oxidation and two LAP1B/LAP1C residues were shown to be dephosphorylated by PP1. This study permitted the identification of the novel human LAP1C isoform and partially unraveled the molecular basis of LAP1 regulation.

The nuclear envelope protein, LAP1B, is a novel protein phosphatase 1 substrate

Protein phosphatase 1 (PP1) binding proteins are quintessential regulators, determining substrate specificity and defining subcellular localization and activity of the latter. Here, we describe a novel PP1 binding protein, the nuclear membrane protein lamina associated polypeptide 1B (LAP1B), which interacts with the DYT1 dystonia protein torsinA. The PP1 binding domain in LAP1B was here identified as the REVRF motif at amino acids 55-59. The LAP1B:PP1 complex can be immunoprecipitated from cells in culture and rat cortex and the complex was further validated by yeast co-transformations and blot overlay assays. PP1, which is enriched in the nucleus, binds to the N-terminal nuclear domain of LAP1B, as shown by immunocolocalization and domain specific binding studies. PP1 dephosphorylates LAP1B, confirming the physiological relevance of this interaction. These findings place PP1 at a key position to participate in the pathogenesis of DYT1 dystonia and related nuclear envelope-based diseases.

Extensive regulation of NAGNAG alternative splicing: new tricks for the spliceosome?

A recent study using massive parallel sequencing demonstrates unequivocally that alternative tandem acceptor splicing is tissue-specifically regulated.

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.

Extensive regulation of NAGNAG alternative splicing: new tricks for the spliceosome?

A recent study using massive parallel sequencing demonstrates unequivocally that alternative tandem acceptor splicing is tissue-specifically regulated.

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.

Alternative splicing at a NAGNAG acceptor site as a novel phenotype modifier

Approximately 30% of alleles causing genetic disorders generate premature termination codons (PTCs), which are usually associated with severe phenotypes. However, bypassing the deleterious stop codon can lead to a mild disease outcome. Splicing at NAGNAG tandem splice sites has been reported to result in insertion or deletion (indel) of three nucleotides. We identified such a mechanism as the origin of the mild to asymptomatic phenotype observed in cystic fibrosis patients homozygous for the E831X mutation (2623G>T) in the CFTR gene. Analyses performed on nasal epithelial cell mRNA detected three distinct isoforms, a considerably more complex situation than expected for a single nucleotide substitution. Structure-function studies and in silico analyses provided the first experimental evidence of an indel of a stop codon by alternative splicing at a NAGNAG acceptor site. In addition to contributing to proteome plasticity, alternative splicing at a NAGNAG tandem site can thus remove a disease-causing UAG stop codon. This molecular study reveals a naturally occurring mechanism where the effect of either modifier genes or epigenetic factors could be suspected. This finding is of importance for genetic counseling as well as for deciding appropriate therapeutic strategies.

Mild disease outcomes associated with premature termination codons can result from at least three different mechanisms, but none of these mechanisms explain the mild phenotype observed in some patients. Subtle differences in alternative transcripts have recently been reported at NAGNAG tandem acceptor motifs, which can be detected in 30% of human genes. We provide the first experimental evidence of premature termination codon removal by alternative splicing at a NAGNAG acceptor splice site. Our study emphasizes the biological significance of such alternative splicing in the context of disease-causing mutations and defines a new phenotype-modifying mechanism that buffers nonsense mutations.

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.

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.

Molecular determinants and evolutionary dynamics of wobble splicing

Alternative splicing at tandem splice sites (wobble splicing) is widespread in many species, but the mechanisms specifying the tandem sites remain poorly understood. Here, we used synaptotagmin I as a model to analyze the phylogeny of wobble splicing spanning more than 300 My of insect evolution. Phylogenetic analysis indicated that the occurrence of species-specific wobble splicing was related to synonymous variation at tandem splice sites. Further mutagenesis experiments demonstrated that wobble splicing could be lost by artificially induced synonymous point mutations due to destruction of splice acceptor sites. In contrast, wobble splicing could not be correctly restored through mimicking an ancestral tandem acceptor by artificial synonymous mutation in in vivo splicing assays, which suggests that artificial tandem splice sites might be incompatible with normal wobble splicing. Moreover, combining comparative genomics with hybrid minigene analysis revealed that alternative splicing has evolved from the 3′ tandem donor to the 5′ tandem acceptor in Culex pipiens, as a result of an evolutionary shift of cis element sequences from 3′ to 5′ splice sites. These data collectively suggest that the selection of tandem splice sites might not simply be an accident of history but rather in large part the result of coevolution between splice site and cis element sequences as a basis for wobble splicing. An evolutionary model of wobble splicing is proposed.

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

Background

Several recent studies indicate that alternative splicing in Arabidopsis and other plants is a common mechanism for post-transcriptional modulation of gene expression. However, few analyses have been done so far to elucidate the functional relevance of alternative splicing in higher plants. Representing a frequent and universal subtle alternative splicing event among eukaryotes, alternative splicing at NAGNAG acceptors contributes to transcriptome diversity and therefore, proteome plasticity. Alternatively spliced NAGNAG acceptors are overrepresented in genes coding for proteins with RNA-recognition motifs (RRMs). As SR proteins, a family of RRM-containing important splicing factors, are known to be extensively alternatively spliced in Arabidopsis, we analyzed alternative splicing at NAGNAG acceptors in SR and SR-related genes.

Results

In a comprehensive analysis of the Arabidopsis thaliana genome, we identified 6,772 introns that exhibit a NAGNAG acceptor motif. Alternative splicing at these acceptors was assessed using available EST data, complemented by a sequence-based prediction method. Of the 36 identified introns within 30 SR and SR-related protein-coding genes that have a NAGNAG acceptor, we selected 15 candidates for an experimental analysis of alternative splicing under several conditions. We provide experimental evidence for 8 of these candidates being alternatively spliced. Quantifying the ratio of NAGNAG-derived splice variants under several conditions, we found organ-specific splicing ratios in adult plants and changes in seedlings of different ages. Splicing ratio changes were observed in response to heat shock and most strikingly, cold shock. Interestingly, the patterns of differential splicing ratios are similar for all analyzed genes.

Conclusion

NAGNAG acceptors frequently occur in the Arabidopsis genome and are particularly prevalent in SR and SR-related protein-coding genes. A lack of extensive EST coverage can be compensated by using the proposed sequence-based method to predict alternative splicing at these acceptors. Our findings indicate that the differential effects on NAGNAG alternative splicing in SR and SR-related genes are organ- and condition-specific rather than gene-specific.

Widespread and subtle: alternative splicing at short-distance tandem sites

Alternative splicing at donor or acceptor sites located just a few nucleotides apart is widespread in many species. It results in subtle changes in the transcripts and often in the encoded proteins. Several of these tandem splice events contribute to the repertoire of functionally different proteins, whereas many are neutral or deleterious. Remarkably, some of the functional events are differentially spliced in tissues or developmental stages, whereas others exhibit constant splicing ratios, indicating that function is not always associated with differential splicing. Stochastic splice site selection seems to play a major role in these processes. Here, we review recent progress in understanding functional and evolutionary aspects as well as the mechanism of splicing at short-distance tandem sites.

Assessing the fraction of short-distance tandem splice sites under purifying selection

Many alternative splice events result in subtle mRNA changes, and most of them occur at short-distance tandem donor and acceptor sites. The splicing mechanism of such tandem sites likely involves the stochastic selection of either splice site. While tandem splice events are frequent, it is unknown how many are functionally important. Here, we use phylogenetic conservation to address this question, focusing on tandems with a distance of 3-9 nucleotides. We show that previous contradicting results on whether alternative or constitutive tandem motifs are more conserved between species can be explained by a statistical paradox (Simpson's paradox). Applying methods that take biases into account, we found higher conservation of alternative tandems in mouse, dog, and even chicken, zebrafish, and Fugu genomes. We estimated a lower bound for the number of alternative sites that are under purifying (negative) selection. While the absolute number of conserved tandem motifs decreases with the evolutionary distance, the fraction under selection increases. Interestingly, a number of frameshifting tandems are under selection, suggesting a role in regulating mRNA and protein levels via nonsense-mediated decay (NMD). An analysis of the intronic flanks shows that purifying selection also acts on the intronic sequence. We propose that stochastic splice site selection can be an advantageous mechanism that allows constant splice variant ratios in situations where a deviation in this ratio is deleterious.