Phylogenetically widespread alternative splicing at unusual GYNGYN donors
© Hiller et al.; licensee BioMed Central Ltd. 2006
Received: 3 April 2006
Accepted: 25 July 2006
Published: 25 July 2006
Splice donor sites have a highly conserved GT or GC dinucleotide and an extended intronic consensus sequence GTRAGT that reflects the sequence complementarity to the U1 snRNA. Here, we focus on unusual donor sites with the motif GYNGYN (Y stands for C or T; N stands for A, C, G, or T).
While only one GY functions as a splice donor for the majority of these splice sites in human, we provide computational and experimental evidence that 110 (1.3%) allow alternative splicing at both GY donors. The resulting splice forms differ in only three nucleotides, which results mostly in the insertion/deletion of one amino acid. However, we also report the insertion of a stop codon in four cases. Investigating what distinguishes alternatively from not alternatively spliced GYNGYN donors, we found differences in the binding to U1 snRNA, a strong correlation between U1 snRNA binding strength and the preferred donor, over-represented sequence motifs in the adjacent introns, and a higher conservation of the exonic and intronic flanks between human and mouse. Extending our genome-wide analysis to seven other eukaryotic species, we found alternatively spliced GYNGYN donors in all species from mouse to Caenorhabditis elegans and even in Arabidopsis thaliana. Experimental verification of a conserved GTAGTT donor of the STAT3 gene in human and mouse reveals a remarkably similar ratio of alternatively spliced transcripts in both species.
In contrast to alternative splicing in general, GYNGYN donors in addition to NAGNAG acceptors enable subtle protein variations.
Given the rather limited number of human genes , alternative splicing is believed to be a major mechanism to bridge the gap between the gene and protein number [2, 3]. Most human multi-exon genes express more than one splice variant . Protein isoforms, produced by alternative splicing, can differ in various aspects, including ligand binding affinity, signaling activity, protein domain composition, subcellular localization, and protein half-life . In coordination with nonsense-mediated mRNA decay, alternatively spliced transcripts can be degraded rapidly, providing a regulation and fine-tuning mechanism of the adjustment of the protein level .
The skipping of an exon is the most frequent alternative splice event, followed by alternative splice donor and acceptor sites . Such splice events often result in large effects for the proteins, for example, by deleting functional units like protein domains [8, 9] or transmembrane helices [10, 11]. On the other hand, alternative splicing also allows the production of many very similar protein isoforms. The most frequent of these subtle events is the alternative splicing at NAGNAG or tandem acceptors . In the NAGNAG motif (N stands for A, C, G or T/U; throughout the paper we write T instead of U also when referring to an RNA sequence), we have termed the upstream acceptor the E acceptor (since the downstream NAG becomes exonic in case of splicing at this site) and the downstream one the I acceptor (since the whole tandem becomes intronic). This splice acceptor motif frequently allows the selection of one of the two AGs in the splice process, resulting in the insertion/deletion (indel) of the I acceptor NAG in mRNAs, preferably if both Ns are either A, C, or T [13–15]. Despite the rather simple genomic structure, these NAG indels lead to a surprisingly high diversity at the protein level. Depending on the sequence of the up- and downstream exon and the phase of the intron, eight different single amino acid indels, the exchange of a dipeptide for an unrelated amino acid, and the indel of a stop codon are possible . These subtle protein changes can result in functional differences for the respective protein isoforms [15–18].
The recognition of donor and acceptor splice sites is different. While the acceptor AG and its preceding polypyrimidine tract is recognized by the U2AF heterodimer , the donor site has an extended consensus sequence AG|GTRAGT (| is the splice site, R stands for A or G), that is bound by base pairing to the 5' end of the U1 snRNA . However, two donor sites that are only three nucleotides (nt) apart would result in overlapping U1 snRNA binding sites and the GTNGTN motif differs from the donor consensus sequence at the two conserved positions +4 and +5. According to the consensus, an alternative usage of the GT dinucleotide 4 nucleotides downstream is much more likely but results in a frameshift and thus a dramatic change of the protein if the donor is located in the coding sequence (CDS).
Here we investigate whether alternative splicing at a GT or GC donor dinucleotide 3nt up- or downstream is possible. This type of alternative splicing requires a GYNGYN donor motif (Y stands for C or T) and is of interest because it would result in similar subtle protein changes like at NAGNAG tandem acceptors and thus increase the proteome plasticity. We found expressed sequence tag (EST) and/or mRNA evidence for alternative splicing at 110 human GYNGYN tandem donors and confirm the existence of both splice forms by RT-PCR experiments in seven cases. We report the occurrence of alternative splicing at GYNGYN tandem donors in six other animals and a plant. Analyzing the GYNGYN motifs that do and do not allow alternative splicing, we found significant differences in the stability of the U1 snRNA binding, conserved exonic and intronic flanks between human and mouse, and over-represented sequence motifs in the intronic flanks.
Alternative splicing at tandem donor sites
Human tandem donor sites divided into the four different GYNGYN patterns
Splice donor pattern
Number and % of tandem donors*
Number and % of confirmed donors
By searching dbEST and the human mRNAs from GenBank, we identified experimental evidence for alternative splicing at 110 (1.3% of 8,550) tandem donors (in the following we term these tandem donors 'confirmed') (Table 1; Additional data file 1). We term the remaining 8,440 donors 'unconfirmed' with the notion that they are enriched in GYNGYN donors that are not functional. The percentage of confirmed tandem donors is considerably higher for GTNGTN (2%) and GTNGCN (1.6%) patterns. No confirmed GCNGCN donor was found, presumably because this motif is very rare and because the weaker GC donor requires a more stringent sequence context. Since ESTs are random high-throughput samples from the transcriptome, spurious or mis-spliced entries may pollute dbEST, especially if the EST number for a particular locus is high [21, 22]. However, the likelihood of splicing errors decreases if the respective splice event is represented by more than one EST and/or if the EST ratio between alternative splice forms is not extreme. From the 110 confirmed tandems, 50 (45%) have at least two ESTs and 19 (17%) have at least five ESTs for e as well as i transcripts. Likewise, in 85 cases (77%) the minor splice form is confirmed by more than 1% of the ESTs that are spliced at the tandem donor, and in 49 cases (45%) this fraction is at least 5%. Thus, although we cannot exclude that some confirmations of GYNGYN tandem donors represent rare errors of the splice machinery, the majority seems to comprise real alternative splice events.
A or G is strongly preferred at intron position +3 for standard donor sites GTN, while T and C have lower frequencies . We classified the confirmed GTNGTN donors according to their pattern into three groups: GTRGTR (R = A or G); GTTGTR, GTRGTT or GTTGTT; and GTCGTN or GTNGTC. The GTRGTR pattern is clearly preferred as 86% (70 of 81) of the confirmed GTNGTN donors belong to this group. A smaller fraction has one or two T at the N-positions (8 of 81, 10%) and the third group is very rare, with only three cases. These findings indicate that the common splicing machinery is operating at these sites. For GTNGCN and GCNGTN donors, we found very similar results: 21 of 29 (72%) have R at both N-positions and two (7%) one T. In addition, we found the exceptional pattern GTAGCC six times (21%).
In some cases it has been reported that single nucleotide polymorphisms (SNPs) in the vicinity of donor sites lead to a shift in the splice site [24–26]. To check if there is a general trend that confirmed GYNGYNs might be influenced by SNPs in their genomic flanks, thus giving rise to allele-specific splice forms , we selected all SNPs from dbSNP that are mapped to the 100 nucleotide context up- and downstream of these tandem donors. We found that 64 (58%) of the confirmed GYNGYNs do not have an annotated SNP in this 206 nucleotide region. As a control we randomly selected 500 unconfirmed GYNGYNs and found that 56% (279 of 500) do not have a SNP in this 206 nucleotide region. Thus, we conclude that most of the confirmed tandems are not associated with allele-specific splice forms.
Experimental verification of alternative GYNGYN splicing
Experimental verification of human GYNGYN donors
e < i
e > i
e > i
e < i
e < i
e > i
e < i
Differences in U1 snRNA binding for confirmed and unconfirmed GYNGYN donors
The U1 snRNA determines the donor site by base pairing with the mRNA . To define the strength of a donor site, we calculated: the average free energy of U1 snRNA binding; the average number of base pairs between donor sites and U1 snRNA ; and the maximum entropy scores . In general, the e donor of confirmed GTNGTNs has a higher strength compared to the i donor (Additional data file 2). In agreement with that, the e donor is annotated in 73% (59 of 81) of the confirmed GTNGTN donors in RefSeq. Furthermore, the e donor is represented by an average of 233 ESTs, which is about tenfold higher than the average of 24 ESTs for the i donor. These findings can be explained with a stronger consensus sequence downstream of a standard GT donor compared to the three upstream positions (Figure 2a). For GTNGCN and GCNGTN tandems, we have to distinguish between the GT and GC donor site since GT is stronger than GC (Additional data file 2). Consistently, of the 29 confirmed GTNGCN and GCNGTN tandems, the GT donor is annotated in 23 cases (79%) in RefSeq and the splicing at the GT donor is represented by an average of 116 ESTs compared to the average of four ESTs for GC donors.
Characteristics of U1 snRNA binding to human confirmed and unconfirmed GYNGYN donors
Free energy (kcal/mol)*
Number of base pairs*
Maximum entropy score†
Unconfirmed, e annotated
Unconfirmed, i annotated
We assumed that the strength of both donors might be a criterion to distinguish functional from non-functional tandem donors. To test this experimentally, we selected nine unconfirmed GTNGTNs with a low free energy for both donor sites for experimental verification. As for confirmed GTNGTNs, RT-PCR products were directly sequenced and the sequencing traces were inspected for overlapping sequences. For none of the nine candidates, we found evidence for alternative splicing at the tandem donor, suggesting that the majority of unconfirmed GTNGTNs is presumably not alternatively spliced. However, our direct sequencing approach does not exclude that the alternative transcript is expressed at a low frequency.
We conclude that: stable U1 binding is necessary but not sufficient for alternative tandem donor splicing; the currently confirmed GYNGYN represent a large fraction of all functional tandem donors; and, in contrast to NAGNAG acceptors , alternatively spliced GYNGYNs are not easily predictable.
Confirmed tandem donors have over-represented motifs in their intron flanks
Since the free energy of U1 binding seems not to be the only discriminative criterion, we searched for other differences between confirmed and unconfirmed GTNGTNs. The regulation of alternative splicing often involves auxiliary exonic and intronic splice enhancer and silencer elements (abbreviated ESE, ESS, ISE, and ISS, respectively) that are bound by trans-acting RNA-binding proteins like serine/arginine rich (SR) proteins and hnRNPs [33–36]. Previous computational studies followed by experimental verification identified 238 hexamers as ESEs , 2,060 octamers as ESEs and 1,019 octamers as ESSs , and 133 hexamers as ISE motifs in the vicinity of donor sites . We used these motifs to compare their average frequency between both groups. The 100 nucleotide exonic flanks of confirmed GTNGTNs are indistinguishable from unconfirmed ones when comparing the frequency of the 238 ESE hexamers (average of 10 ESEs per exon flank for both groups) and have a slight but not significantly higher frequency of ESE octamers (average 7.2 versus 6.5). The ESS frequency is slightly, but not significantly, lower for exon flanks of confirmed tandem donors (average 1.3 versus 2). However, we found a significantly higher frequency of ISE motifs in the 100 nucleotide intron flanks for confirmed GTNGTNs (average 10 versus 8, t-test: P = 0.0174). We repeated this analysis using a shorter exonic/intronic context (50 nt) and found consistent results (data not shown).
To find out if specific ISE hexamers are statistically over-represented, we used a resampling strategy. We randomly sampled 10,000 sets, each comprising 81 intron flanks from unconfirmed GTNGTNs. We estimated the P-value as the fraction of random sets with a higher frequency of a given ISE hexamer compared to the observed frequency in confirmed tandem donors. CGGGGT is the only one among the 133ISE motifs that is significantly over-represented in the vicinity of confirmed GTNGTN donors as all 10,000 random sets have a lower frequency (P < 1/10,000 × 133 = 0.0133 to correct for multiple testing). To find out if other sequence motifs are over-represented in the intron flanks of confirmed tandem donors, we repeated this procedure with tetramers. A word length of 4 nucleotides was chosen to account for the rather small data set. We only compared the 119 tetramers that occur at least with the expected frequency in the intron flanks of confirmed GTNGTNs. We found a significant overrepresentation for GGGT and CGGG (both have a higher frequency in only two random sets, P < 3/10,000 × 119 = 0.0357), while the tetramer GGGG has a corrected P-value slightly above the 0.05 threshold (higher frequency in five random sets, P < 0.0714). Since both GGGT and CGGG are substrings of the over-represented ISE CGGGGT, no new sequence motifs were found. The common feature of the over-represented sequence motifs is the G triplet. Interestingly, this motif occurs in 82 of the 133 ISEs  and is a known splice enhancer . Since both splice sites of confirmed GTNGTNs are weaker compared to the annotated splice site of unconfirmed ones (Table 3), the G triplets might simply be associated with weak GTNGTNs. To exclude this possibility, we compared the average GGG frequency with unconfirmed GTNGTNs having a low U1 binding potential for both e and i donor (average free energy -3 kcal/mol for the e donor, -2.2 for the i donor) and still found an over-representation in the intron flanks of confirmed GTNGTNs (average 4.4 versus 2.6 G triplets per intron flank). Since this triplet was found to be more frequent in shorter introns , we divided our confirmed and unconfirmed datasets into short and long introns using 200 nucleotide as a cut-off. Consistently, the GGG is more frequent in the flanks of short as well as long introns with confirmed GTNGTNs (average 8.3 versus 4.4 G triplets per short intron, average 3.4 vs 2.7 per long intron). We also observed a noticeable higher ISE hexamer frequency in the intronic flanks of confirmed GTNGCN and GCNGTN tandems (average 12.7 versus 10.8, not significant), but only a slightly higher frequency of ESE hexamers and octamers in their exon flanks (data not shown). Specifically, G triplets are also more frequent in the intronic flanks of confirmed GTNGCN and GCNGTN donors compared to unconfirmed ones (average 5 versus 4.1). Thus, the occurrence of G triplets is another discriminating criterion between confirmed and unconfirmed tandem donors.
Protein impact of alternative splicing at GYNGYN donors
Of the 81 confirmed GTNGTNs, 72 (89%) are located downstream of a coding exon; thus, alternative splicing at these sites results in 3 nucleotide indels into the coding sequence. The effect for the protein depends on the phase of the intron as well as the sequence of the GTNGTN and the upstream/downstream exon. In intron phase 0 (intron location between two codons) the GTN of the i donor is inserted/deleted and codes for a valine. In intron phase 1 and 2 (location between the first and second codon position, respectively), three different events are possible: indel of a single amino acid; exchange of a dipeptide and a different amino acid; and indel of a stop codon. Of the 72 GTNGTNs, 37 (51%) are located in phase 0, thus a valine indel is the most frequent event at the protein level. Of the 28 (39%) GTNGTNs in phase 1, 18 result in single amino acid events (14 times glycine, 2 times arginine, 2 times serine), 8 exchange a dipeptide and an unrelated amino acid and in two cases splicing at the i donor creates a stop codon. The 7 (10%) confirmed GTNGTNs in phase 2 are interesting since they either result in indels of rare amino acids (three times tryptophan, one cysteine, one tyrosine) or insert/delete a stop codon in two cases.
Thus, alternative splicing at four tandem donors has a more drastic impact on the proteins by the indel of a stop codon (phase 1: FAM65A, NM_024519, exon 21; BRSK2, NM_003957, exon 19; phase 2: ABC1, NM_022070, exon 8; KLHL5, NM_001007075, exon 3). In two cases (KLHL5, ABC1) the splice form with the premature stop codon is a clear candidate for nonsense-mediated mRNA decay , which potentially results in a down-regulation of the protein level. For the other two cases, the tandem donor affects the last intron of the transcript, thus the stop codon-containing splice variant should be translated into a protein with a shortened carboxyl terminus (FAM65A 33 residue deletion at the carboxyl terminus, BRSK2 7 residue deletion).
The protein impact of the 28 GTNGCN and GCNGTN tandems that are located within the CDS comprise 22 single amino acid indels (intron phase 0: alanine, valine; phase 1: arginine, glycine, serine; phase 2: arginine, glutamine, leucine) and 6 dipeptide exchanges (all in phase 1). Compared to the protein events for GTNGTNs, the indel of alanine, glutamine and leucine is only observed for tandems with a GC donor.
Next, we compared the frequency of single amino acid events in phase 1 and 2 for confirmed and as a control for unconfirmed GTNGTNs. While only 42% (495 of 1,180) unconfirmed tandem donors in phase 1 result in a single residue indel, this percentage is significantly higher for confirmed tandems (64%, 18 of 28, Fisher's exact test: P = 0.02). The small number of phase 2 tandems does not allow a significant result, although the same trend is visible (100%, 5 of 5 confirmed tandems; 76%, 431 of 566 unconfirmed tandems; leaving stop codon events out). These findings argue for a preference to insert/delete only single amino acids, presumably because this is a less dramatic event compared to dipeptide exchanges. However, we cannot exclude the possibility that this result is an indirect consequence of a sequence bias of the GTNGTN motif and its context for confirmed tandems that primarily aims at a more stable U1 snRNA binding. For GTNGCN and GCNGTN tandems, only phase 2 donors result in a high percentage of single amino acid indels (100%, 5 of 5 confirmed; 81%, 781 of 962 unconfirmed) while phase 1 events show no bias (40%, 4 of 10 confirmed; 48%, 524 of 1103 unconfirmed).
Tandem donors in seven other species
GTNGTN, GTNGCN and GCNGTN donors in eight investigated species
Number of donors*
GTNGCN and GCNGTN
Number of confirmed GTNGTN donors divided into three groups according to their motif
GTRGTT + GTTGTR + GTTGTT
GTCGTN + GTNGTC
Conservation of exonic and intronic flanks in mouse
Having observed several alternative GTNGTN splice events in human and mouse, we found conservation of the GTNGTN motif for 53 (65.4%) of the 81 human confirmed GTNGTNs. To assess whether this percentage is high or not, we counted GTNGTN conservation for the 3,909 unconfirmed tandems (162 of the 4,071 have no orthologous locus in mouse) and found a very similar percentage of 65.5% (2,561 of 3,909). The fraction of tandem donors that have a completely identical GTNGTN pattern in mouse is also equal: 40 of 81 (49.4%) confirmed, 1,939 of 3,909 (49.6%) unconfirmed. Thus, there is no evidence for a general selection pressure to maintain confirmed tandem donors since the divergence of the human-mouse ancestor.
However, a considerable fraction (10 of 53 (19%)) of the conserved and confirmed human GTNGTNs is also confirmed in mouse. For example, the GTAGTT donor of intron 21 of STAT3 is conserved in mouse and both e and i transcripts are supported by mouse ESTs. As in humans, we performed RT-PCR in several mouse tissues to further support the EST-derived confirmation of alternative splice events (Figure 3b). We found experimental evidence for alternative splicing at the Stat3 tandem donor in all investigated tissues and observed a strikingly similar trace pattern in human and mouse (Figure 3, e+i). Accordingly, the ratio of e and i transcripts estimated by the EST/mRNA counts are virtually identical (e transcripts 57 of 74 (77%)) human ESTs vs 55 of 69 (79.7%) mouse ESTs). To accurately quantify the ratio of e and i transcripts in one selected tissue, we subcloned the RT-PCR product, sequenced individual clones and found a remarkable agreement in the transcript ratio: 82.8% of the human clones indicate splicing at the e donor, which is almost equal to 85.3% in the mouse (Figure 3, e:i). Interestingly, this tandem donor is conserved in several other mammals and the e:i ratio is very similar (9:2 ESTs for rat, 12:3 ESTs for cow, 9:1 ESTs for dog). This indicates that, in addition to the tandem donor, regulatory elements may be conserved.
We report the occurrence of alternative splice donor usage for GTNGTN, GTNGCN, and GCNGTN motifs in eight investigated eukaryotic species. Apart from our experimental verification of seven human and one mouse GYNGYN donors, several lines of evidence indicate that the majority of observed events is attributable to real alternative splicing. Firstly, numerous GTNGTNs are confirmed by multiple ESTs/mRNAs and for several of these events both e and i transcripts are deposited in the RefSeq database. Secondly, the existence of orthologous tandem donors that are confirmed in two or more species makes EST artifacts or database errors unlikely. Thirdly, these GTNGTN donors have a higher conservation of the exonic and intronic flanking regions, a situation that is typical for conserved alternative splice events [45, 46, 48, 49]. Fourthly, all of the six investigated human individuals express e and i transcripts for STAT3, thus excluding the possibility of allele-specific instead of alternative splicing . Finally, by manual examination of all human confirmed GYNGYNs, we excluded the existence of paralogs or processed pseudogenes that could mimic alternative splicing at a tandem donor.
NAGNAG acceptors for seven species
Number of acceptors*
Confirmed NAGNAG acceptor
Confirmed HAGHAG acceptor
Although only a fraction of the tandem donors is confirmed, we found features that distinguish confirmed from unconfirmed ones. Since the non-annotated donor of unconfirmed tandems does not allow a sufficiently stable binding to the U1 snRNA, the other donor is used exclusively in the splice process. For confirmed tandem donors, both sites allow a stable binding to U1 snRNA. However, in most of the confirmed cases one donor has a better strength and this results in its preferred usage as measured by the EST ratio between both transcripts. The second discriminative feature is the overabundance of G triplets in the intronic flanks of confirmed GTNGTNs, especially for introns shorter than 200 nt. This triplet is the core of many known ISE motifs [39, 40] and was demonstrated to function in splice site definition . Interestingly, in the human alpha-globin gene, GGG elements were shown to exert their effect by binding to the nucleotides 8-10 (5'-CCT-3') of the U1 snRNA . We have searched for over-represented tetramers and found a significantly higher frequency of CGGG and GGGT. Strikingly, the nucleotides 7-11 of U1 snRNA are 5'-ACCTG-3'. The CGGG as well as the GGGT motifs are complementary to this part of U1; thus, it is tempting to speculate that these motifs bind to U1 snRNA with four instead of three base pairs. Since CGGG and GGGT are more frequent in the intronic flanks of confirmed tandem donors, they may be involved in alternative splicing at these donor sites. If U1 snRNA is a critical factor, we do not expect much variation in splicing between tissues since U1 is ubiquitously expressed in high amounts. Consistent with this notion, the seven experimentally investigated tandem donors exhibit similar e:i transcript ratios in all tissues.
Most confirmed GYNGYNs have a low free energy of U1 snRNA binding to both the e and i donor, suggesting that the U1 snRNA can stably bind to both sites. However, there are a few exceptions where one donor is much stronger than the other one in a confirmed tandem (Figure 4a). The mechanism of splicing at these sites remains unclear but there are several hypotheses that might guide future investigations. For example, it has been reported that U6 snRNA rather than U1 snRNA determines a donor site in the human FGFR1 gene . Moreover, there is evidence that splicing can occur without U1 snRNA binding to the donor site [54, 55]. Furthermore, other protein factors can influence the splice site choice and/or (de)stabilize U1 snRNA binding [56, 57]. We believe that a further experimental investigation of confirmed tandem splice donors may help to elucidate further details of the splicing process.
Previously, we found that the impact of SNPs in NAGNAG acceptors on alternative splicing can be accurately predicted . Therefore, it would be interesting to check if similar statements are possible for SNPs in GYNGYN donors. In principle, a SNP in close proximity to an unconfirmed GYNGYN donor might increase the base pairing capability to the U1 snRNA for the alternative donor, thus enabling alternative splicing. SNPs that affect a confirmed tandem donor might weaken U1 binding for one donor and result in the exclusive usage of the other. During the SNP mapping, we found two SNPs in the GTNGTN motif of human confirmed tandem donors. For exon 4 of FAM3B (NM_206964), the verified SNP rs417708 results in two alleles, GTGGTA and GCGGTA. For the C allele, the free energy of the i and e donor is -3.4 and -6.2 kcal/mol, respectively, while this value is more balanced for the T allele (-5.3 for i and -5.9 kcal/mol for the e donor). This agrees well with the prediction of the splice site analysis server [58, 59]. Thus, it is likely that only the T allele produces two splice forms. The second SNP (rs11672749) is especially interesting, since it affects the tandem donor of exon 5 of the maternally imprinted PEG3 gene (NM_006210, GTGGTG and GTGGGG alleles). Both homozygous G genotypes and heterozygous genotypes with a maternally inactivated T allele will result in an exclusive splicing at the i donor.
Higher eukaryotes typically express multiple transcripts and proteins from a single gene. A prominent mechanism is alternative splicing as about 74% of the human multi-exon genes express more than one splice variant . Protein isoforms can also be expressed from paralogous genes. Large gene families are observed to have a reduced frequency of alternative splicing, consistent with the idea that the variability of those gene products comes from the divergence of the gene copies . While most research focused on large changes introduced by alternative splicing, it is becoming clear that there is a surprisingly high number of very similar protein isoforms. There are several mechanisms to introduce subtle protein changes. The most widespread type is alternative splicing at NAGNAG acceptors [12, 15]. Furthermore, very similar mutually exclusive exons can lead to similar but functionally different proteins . Here, we found that alternative splicing at GYNGYN donor sites occurs in all eight investigated species. Despite not as frequent as confirmed NAGNAG acceptors, the diverse protein changes further contribute to the plasticity of these proteomes. Confirmed tandem donors and acceptors are able to insert 12 of the 20 different amino acids by single amino acid events and the dipeptide exchanges are even more diverse. Further flexibility comes from the simultaneous use of a GYNGYN donor and a NAGNAG acceptor for one intron (Figure 1b). Such an example is intron 9 of BRUNOL4 (NM_020180), for which we found 14 e-E, 3 i-E and 6 e-I transcripts in dbEST that result in protein forms with a GPA, AA, or GP peptide, respectively.
Despite many GYNGYN donors, we found only a minority that allows alternative splicing. Nevertheless, among the human confirmed and evolutionary conserved tandem donors we found a significant fraction to be confirmed in other species. Moreover, the splicing pattern of the STAT3 GTNGTN donor is strikingly equal in human and mouse. In light of the discussion about functional versus non-functional alternative splicing [21, 62], this is a strong indication that these alternative splicing events are not splicing noise. Consistently, such subtle changes by alternative splicing may result in functional differences for the two proteins. An arginine insertion between two zinc fingers results in a human glucocorticoid receptor isoform (NM_001018075, exon 3, GTAGTG) with an activity reduced to 48% [63, 64]. Interestingly, this tandem donor is also conserved and confirmed in mouse. A similar subtle 6 nucleotide shift at a GTAAATGT donor of ALDH18A1 results in an isoform that is insensitive to ornithine inhibition . Furthermore, there are at least four reported cases of functional differences by alternative NAGNAG splicing [15–18]. Thus, subtle alternative splice events are interesting candidates for further research, especially since several of them occur in known disease genes .
Materials and methods
Identification of GYNGYN donors and NAGNAG acceptors
We downloaded from the UCSC Genome Browser  the Human genome assembly (hg17, May 2004) as well as RefSeq annotation (refGene.txt.gz, November 2005). We discarded transcripts with an erroneous open reading frame or ambiguous characters in their sequence. If more than one entry with the same RefSeq ID exists, we selected the transcript with the highest number of exons. From the transcripts, we extracted a list of unique genomic positions of donor sites. We checked the presence of a GT/GC dinucleotide at the annotated donor position and 3 nucleotides up- or downstream. We repeated this procedure for the other six eukaryotes having a RefSeq annotation in the UCSC Genome Browser. The respective genome assemblies are mm6 (March 2005) for M. musculus, rn3 (June 2003) for R. norvegicus, galGal2 (February 2004) for G. gallus, danRer2 (June 2004) for D. rerio, dm2 (April 2004) for D. melanogaster, and ce2 (March 2004) for C. elegans. For A. thaliana, we used the genome assembly (NCBI, build 5.0) to build a list of donor sites according to the CDS feature annotations in GenBank format and screened this list for donor sites with a GYNGYN sequence pattern.
For checking if a GYNGYN donor is confirmed (at least one EST/mRNA for the e as well as i transcript), we compiled two search strings: 30 nucleotides from the upstream exon and 30 nucleotides from the downstream exon for the i transcript; and 30 nucleotides from the upstream exon, the GYN of the i donor and 30 nucleotides from the downstream exon for the e transcript. Then, we used BLAST against all ESTs and mRNAs for the respective species, allowing at most one mismatch or one gap but demanding exact identity for the region 27-33 for i-transcripts and 27-36 for e-transcripts. The ESTs were downloaded from the UCSC Genome Browser (est.fa.gz, November 2005). mRNAs were downloaded from GenBank at the same date as the ESTs. For acceptor sites with a NAGNAG pattern, we repeated the analysis using analogous procedures and the same data described above.
Conservation analysis in mouse
Human-mouse genomic alignments (hg17-mm6) were downloaded from the UCSC Genome Browser (vsMm6/axtNet, March 2005). We used the genomic position of human and mouse donor sites to select the respective alignment chain. From the alignments, we determined whether a human GTNGTN donor is conserved (there is also a GTNGTN motif in mouse) or completely identical. For the per-position identity computation, we considered the alignment part up to 100 positions upstream and downstream. For each position, we counted how often there is identity between human and mouse (Nid), and how often there is a mismatch or gap (Nmm). The per-position identity value is Nid/(Nid+Nmm). Alignment positions with an 'N' aligned to a nucleotide were ignored.
To find tandem donors that are orthologous and confirmed in human and mouse, we used BLAST with the human-confirmed search strings against the search strings of the mouse-confirmed GTNGTNs. Furthermore, we used BLAST with the human-confirmed search strings against the mouse ESTs and mRNAs. Using the UCSC and Ensembl genome browser, we manually checked each hit with an E-value of less than 1e-3 for being alternatively spliced in both species and for having a true orthologous relationship.
Strength of a donor site
We extracted a 9 nucleotide genomic context (3 nucleotides upstream to 3 nucleotides downstream of the GYN) for the e and i donor of confirmed and unconfirmed GYNGYNs. The free energy and number of base pairs with the U1 snRNA were computed according to  with the Splice-site Analyzer tool . The score according to the maximum entropy model  was computed using MaxEntScan .
We extracted the genomic sequence 100 nucleotides upstream (exonic) and 100 nucleotides downstream (intronic) of GTNGTN donor motifs. To identify over-represented ISE hexamer motifs, we used a resampling procedure to estimate the P-value for a higher frequency in the intronic flanks of confirmed GTNGTNs. To this end we randomly sampled 10,000 sets of 81 intronic flanks of unconfirmed GTNGTNs and computed the frequency for each of the 133 ISE motifs in the 10,000 random sets. The P-value for one ISE is the fraction of random sets with a higher frequency compared to the observed frequency for confirmed GTNGTNs. To correct for multiple testing, each P-value is multiplied by 133. For the general search for over-represented motifs, we decided to use tetramers (word length 4 nt) instead of hexamers since the dataset of the confirmed tandems is rather small. Since we were searching for over-represented motifs in the intronic flanks of confirmed GTNGTNs, we expected that such motifs occur at least with the expected frequency under a null model and with a significant higher frequency compared to the flanks of unconfirmed GTNGTNs. There are 97 overlapping tetramers in a 100 nucleotide sequence, thus we analyzed a total of 81 × 97 = 7,857 tetramer occurrences. For complete random sequences, each tetramer should occur 7857/256 = 30.7 times. Since intron sequences are not random, we found a total of 119 tetramers that occur 30 times or more in the flanks of confirmed GTNGTNs. For these 119 tetramers, we repeated the procedure described above but multiplied the P-value by 119.
Experimental verification of alternative splicing at tandem donors
Eight genes with multiple EST evidence for alternative splicing at a tandem donor were analyzed by RT-PCR in different tissues by using cDNA from multiple tissue cDNA panels (BD Clontech Germany, Heidelberg, Germany) as PCR templates. Primers were designed for the exons flanking the tandem donor with distances to these donors that allow reliable amplification and sequencing. PCR was performed in a total volume of 25 μl using ReadyToGo PCR beads (GE Healthcare Europe, Munich, Germany) with 5 pmoles of each primer and 1 μl of cDNA. Cycling conditions were 94°C for 30 s followed by 35 cycles with 94°C for 20 s, 57°C for 30 s and 72°C for 30 s, followed by a final extension at 72°C for 10 minutes. Amplified fragments were precipitated with ethanol and ammonium acetate, washed with ethanol and sequenced using DyeTerminator chemistry (Applied Biosystems, Foster city, USA) and the respective PCR primers on a 3730 xl DNA Analyzer (Applied Biosystems). Genes and their primer sequences were: TOM1 (AGTTTGACATGTTTGCGCTG, GCAGCCTTAACACCAGAGGA); STAT3 (GCCATCTTGAGCACTAAGCC, GGTTCAGCACCTTCACCATT); ANAPC4 (AGATGCTGCAGGAATCGAAG, CTGGCTTTTGCAAACACTGA); RBM10 (AGGCTGGATCAGCAGACACT, TCCCTCTTAGAACCCTTGGC); ANGPT1 (ACAAGGAAGAGTTGGACACC, GGGATTTCCAAAACCCATTT); SEMA5B (AGCACGTCCTGTGGCATC, GTCCTCGTCTCGGTCCTTCT); CXorf44 (GAGGGCAGGACTATGGGAG, AAATACTTCTCCTTCATAGCGGA); and LTBP1 (GGACCTGTATTTGTCAAGCCA, TAATGCAGTGTCCTGCTCCA). In addition, Stat3 of M. musculus was analyzed in the respective Clontech mouse tissue panel by amplifying and sequencing the homologous region with the oligonucleotides GCCATCCTAAGCACAAAGCC and GGCTCAGCACCTTCACCGTT. All sequence traces have been deposited in the NCBI Trace Archive (TIs human: 1166719658-1166720385; mouse; 1166879453-1166879628). To estimate the relative amounts of e and i transcripts of human and mouse STAT3, we cloned the respective amplicons into pCR2.1-TOPO (Invitrogen, Karlsruhe, Germany) and propagated the clones in Escherichia coli TOP10 cells. Plasmids were isolated from several isolated clones and their inserts were sequenced using universal M13 primers.
The same strategy was applied to a set of genes with unconfirmed tandem splice donors. Genes and their primers were: A2ML1 (NM_144670, exon 17, ACTTTCCTCAGCCCCTCATT, AGTGCAGAAACTCATCGCCT); TMEM63C (NM_020431, exon 1, GTGCTGAGGACGCAAATCA, CATCTCCAAGGAAGTCTCCG); RNPC3 (NM_017619, exon 13, ACCGGGTGAACCAAACTGTA, AGCTGTTACGCACAGTTCCA); GOLGA3 (NM_005895, exon 5, CACCCCCTATATGGTCAACG, CACGACTGCTTCAGGGTGT); ART5 (NM_053017, exon 2, GCCCCTATACAGGCCTTCTC, ATTGCAACACCGTTCAATCA); KIAA1853 (NM_194286, exon 8, CCCTCAAGCTGTGAGAGCAG, TGGTGAAGGAGTTCCCTGAA); K5B (NM_173352, exon 5, ACAACAACCGCTACCTGGAC, AATCTCCACATCCAGGGAAA); SAMD4B (NM_018028, exon 15, AACAGCATGCCCAGTCAGA, CTCAGCAGAGATCCCTCGAC); and LOC221711 (NM_194299, exon 9, TTCAGATGTTGGATTCCTTCC, TTTTTCATCCGCTGGTTTTC).
To facilitate further experimental and computational studies of tandem splicesites, we recently developed a database, TassDB , which provides large collections of GYNGYN donors and NAGNAG acceptors of eight species.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is an Excel spreadsheet containing data on all human confirmed GYNGYN splice donor sites identified in this study. Additional data file 2 is an Excel spreadsheet providing data about the strength of e and i human GYNGYN donors. Additional data file 3 is an Excel spreadsheet listing all confirmed GYNGYN donors for seven species. Additional data file 4 is an Excel spreadsheet presenting information about selected sequence traces that exemplify the experimental verification of GYNGYN donors. Additional data file 5 is a Word file describing the control experiments for detecting minor splice forms by direct sequencing of RT-PCR products.
We thank Gene Yeo for providing the ISE hexamer list and Anke Busch for critical reading of the manuscript. The skillful technical assistance of Beate Szafranski and Ivonne Görlich is gratefully acknowledged. This work was supported by grants from the German Ministry of Education and Research to SS (01GS0426) and to MP (01GR0504, 0313652D) as well as from the Deutsche Forschungsgemeinschaft (SFB604-02) to MP.
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.PubMedView ArticleGoogle Scholar
- Maniatis T, Tasic B: Alternative pre-mRNA splicing and proteome expansion in metazoans. Nature. 2002, 418: 236-243. 10.1038/418236a.PubMedView ArticleGoogle Scholar
- Graveley BR: Alternative splicing: increasing diversity in the proteomic world. Trends Genet. 2001, 17: 100-107. 10.1016/S0168-9525(00)02176-4.PubMedView ArticleGoogle Scholar
- Johnson JM, Castle J, Garrett-Engele P, Kan Z, Loerch PM, Armour CD, Santos R, Schadt EE, Stoughton R, Shoemaker DD: Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science. 2003, 302: 2141-2144. 10.1126/science.1090100.PubMedView ArticleGoogle Scholar
- Stamm S, Ben-Ari S, Rafalska I, Tang Y, Zhang Z, Toiber D, Thanaraj TA, Soreq H: Function of alternative splicing. Gene. 2005, 344: 1-20. 10.1016/j.gene.2004.10.022.PubMedView ArticleGoogle Scholar
- Lewis BP, Green RE, Brenner SE: Evidence for the widespread coupling of alternative splicing and nonsense-mediated mRNA decay in humans. Proc Natl Acad Sci USA. 2003, 100: 189-192. 10.1073/pnas.0136770100.PubMedPubMed CentralView ArticleGoogle Scholar
- Stamm S, Zhu J, Nakai K, Stoilov P, Stoss O, Zhang MQ: An alternative-exon database and its statistical analysis. DNA Cell Biol. 2000, 19: 739-756. 10.1089/104454900750058107.PubMedView ArticleGoogle Scholar
- Kriventseva EV, Koch I, Apweiler R, Vingron M, Bork P, Gelfand MS, Sunyaev S: Increase of functional diversity by alternative splicing. Trends Genet. 2003, 19: 124-128. 10.1016/S0168-9525(03)00023-4.PubMedView ArticleGoogle Scholar
- Resch A, Xing Y, Modrek B, Gorlick M, Riley R, Lee C: Assessing the impact of alternative splicing on domain interactions in the human proteome. J Proteome Res. 2004, 3: 76-83. 10.1021/pr034064v.PubMedView ArticleGoogle Scholar
- Xing Y, Xu Q, Lee C: Widespread production of novel soluble protein isoforms by alternative splicing removal of transmembrane anchoring domains. FEBS Lett. 2003, 555: 572-578. 10.1016/S0014-5793(03)01354-1.PubMedView ArticleGoogle Scholar
- Hiller M, Huse K, Platzer M, Backofen R: Creation and disruption of protein features by alternative splicing - a novel mechanism to modulate function. Genome Biol. 2005, 6: R58-10.1186/gb-2005-6-7-r58.PubMedPubMed CentralView ArticleGoogle Scholar
- Hiller M, Huse K, Szafranski K, Jahn N, Hampe J, Schreiber S, Backofen R, Platzer M: Widespread occurrence of alternative splicing at NAGNAG acceptors contributes to proteome plasticity. Nat Genet. 2004, 36: 1255-1257. 10.1038/ng1469.PubMedView ArticleGoogle Scholar
- Akerman M, Mandel-Gutfreund Y: Alternative splicing regulation at tandem 3' splice sites. Nucleic Acids Res. 2006, 34: 23-31. 10.1093/nar/gkj408.PubMedPubMed CentralView ArticleGoogle Scholar
- Hiller M, Huse K, Szafranski K, Jahn N, Hampe J, Schreiber S, Backofen R, Platzer M: Single-nucleotide polymorphisms in NAGNAG acceptors are highly predictive for variations of alternative splicing. Am J Hum Genet. 2006, 78: 291-302. 10.1086/500151.PubMedPubMed CentralView ArticleGoogle Scholar
- Tadokoro K, Yamazaki-Inoue M, Tachibana M, Fujishiro M, Nagao K, Toyoda M, Ozaki M, Ono M, Miki N, Miyashita T, Yamada M: Frequent occurrence of protein isoforms with or without a single amino acid residue by subtle alternative splicing: the case of Gln in DRPLA affects subcellular localization of the products. J Hum Genet. 2005, 50: 382-394. 10.1007/s10038-005-0261-9.PubMedView ArticleGoogle Scholar
- Lorkovic ZJ, Lehner R, Forstner C, Barta A: Evolutionary conservation of minor U12-type spliceosome between plants and humans. RNA. 2005, 11: 1095-1107. 10.1261/rna.2440305.PubMedPubMed CentralView ArticleGoogle Scholar
- Vogan KJ, Underhill DA, Gros P: An alternative splicing event in the Pax-3 paired domain identifies the linker region as a key determinant of paired domain DNA-binding activity. Mol Cell Biol. 1996, 16: 6677-6686.PubMedPubMed CentralView ArticleGoogle Scholar
- Condorelli G, Bueno R, Smith RJ: Two alternatively spliced forms of the human insulin-like growth factor I receptor have distinct biological activities and internalization kinetics. J Biol Chem. 1994, 269: 8510-8516.PubMedGoogle Scholar
- Wu S, Romfo CM, Nilsen TW, Green MR: Functional recognition of the 3' splice site AG by the splicing factor U2AF35. Nature. 1999, 402: 832-835. 10.1038/45996.PubMedView ArticleGoogle Scholar
- Zhuang Y, Weiner AM: A compensatory base change in U1 snRNA suppresses a 5' splice site mutation. Cell. 1986, 46: 827-835. 10.1016/0092-8674(86)90064-4.PubMedView ArticleGoogle Scholar
- Sorek R, Shamir R, Ast G: How prevalent is functional alternative splicing in the human genome?. Trends Genet. 2004, 20: 68-71. 10.1016/j.tig.2003.12.004.PubMedView ArticleGoogle Scholar
- Kan Z, States D, Gish W: Selecting for functional alternative splices in ESTs. Genome Res. 2002, 12: 1837-1845. 10.1101/gr.764102.PubMedPubMed CentralView ArticleGoogle Scholar
- Abril JF, Castelo R, Guigo R: Comparison of splice sites in mammals and chicken. Genome Res. 2005, 15: 111-119. 10.1101/gr.3108805.PubMedPubMed CentralView ArticleGoogle Scholar
- Barbaux S, Niaudet P, Gubler MC, Grunfeld JP, Jaubert F, Kuttenn F, Fekete CN, Souleyreau-Therville N, Thibaud E, Fellous M, McElreavy K: Donor splice-site mutations in WT1 are responsible for Frasier syndrome. Nat Genet. 1997, 17: 467-470. 10.1038/ng1297-467.PubMedView ArticleGoogle Scholar
- Valentonyte R, Hampe J, Huse K, Rosenstiel P, Albrecht M, Stenzel A, Nagy M, Gaede KI, Franke A, Haesler R, et al: Sarcoidosis is associated with a truncating splice site mutation in BTNL2. Nat Genet. 2005, 37: 357-364. 10.1038/ng1519.PubMedView ArticleGoogle Scholar
- Costa DB, Lozovatsky L, Gallagher PG, Forget BG: A novel splicing mutation of the alpha-spectrin gene in the original hereditary pyropoikilocytosis kindred. Blood. 2005, 106: 4367-4369. 10.1182/blood-2005-05-1813.PubMedPubMed CentralView ArticleGoogle Scholar
- Nembaware V, Wolfe KH, Bettoni F, Kelso J, Seoighe C: Allele-specific transcript isoforms in human. FEBS Lett. 2004, 577: 233-238. 10.1016/j.febslet.2004.10.018.PubMedView ArticleGoogle Scholar
- Carmel I, Tal S, Vig I, Ast G: Comparative analysis detects dependencies among the 5' splice-site positions. RNA. 2004, 10: 828-840. 10.1261/rna.5196404.PubMedPubMed CentralView ArticleGoogle Scholar
- Yeo G, Burge CB: Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals. J Comput Biol. 2004, 11: 377-394. 10.1089/1066527041410418.PubMedView ArticleGoogle Scholar
- Sorek R, Lev-Maor G, Reznik M, Dagan T, Belinky F, Graur D, Ast G: Minimal conditions for exonization of intronic sequences: 5' splice site formation in alu exons. Mol Cell. 2004, 14: 221-231. 10.1016/S1097-2765(04)00181-9.PubMedView ArticleGoogle Scholar
- Roca X, Sachidanandam R, Krainer AR: Determinants of the inherent strength of human 5' splice sites. RNA. 2005, 11: 683-698. 10.1261/rna.2040605.PubMedPubMed CentralView ArticleGoogle Scholar
- Bi J, Xia H, Li F, Zhang X, Li Y: The effect of U1 snRNA binding free energy on the selection of 5' splice sites. Biochem Biophys Res Commun. 2005, 333: 64-69. 10.1016/j.bbrc.2005.05.078.PubMedView ArticleGoogle Scholar
- Blencowe BJ: Exonic splicing enhancers: mechanism of action, diversity and role in human genetic diseases. Trends Biochem Sci. 2000, 25: 106-110. 10.1016/S0968-0004(00)01549-8.PubMedView ArticleGoogle Scholar
- Ladd AN, Cooper TA: Finding signals that regulate alternative splicing in the post-genomic era. Genome Biol. 2002, 3: reviews0008.1-0008.16. 10.1186/gb-2002-3-11-reviews0008.View ArticleGoogle Scholar
- Graveley BR: Sorting out the complexity of SR protein functions. RNA. 2000, 6: 1197-1211. 10.1017/S1355838200000960.PubMedPubMed CentralView ArticleGoogle Scholar
- Pozzoli U, Sironi M: Silencers regulate both constitutive and alternative splicing events in mammals. Cell Mol Life Sci. 2005, 62: 1579-1604. 10.1007/s00018-005-5030-6.PubMedView ArticleGoogle Scholar
- Fairbrother WG, Yeh RF, Sharp PA, Burge CB: Predictive identification of exonic splicing enhancers in human genes. Science. 2002, 297: 1007-1013. 10.1126/science.1073774.PubMedView ArticleGoogle Scholar
- Zhang XH, Chasin LA: Computational definition of sequence motifs governing constitutive exon splicing. Genes Dev. 2004, 18: 1241-1250. 10.1101/gad.1195304.PubMedPubMed CentralView ArticleGoogle Scholar
- Yeo G, Hoon S, Venkatesh B, Burge CB: Variation in sequence and organization of splicing regulatory elements in vertebrate genes. Proc Natl Acad Sci USA. 2004, 101: 15700-15705. 10.1073/pnas.0404901101.PubMedPubMed CentralView ArticleGoogle Scholar
- McCullough AJ, Berget SM: An intronic splicing enhancer binds U1 snRNPs to enhance splicing and select 5' splice sites. Mol Cell Biol. 2000, 20: 9225-9235. 10.1128/MCB.20.24.9225-9235.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- McCullough AJ, Berget SM: G triplets located throughout a class of small vertebrate introns enforce intron borders and regulate splice site selection. Mol Cell Biol. 1997, 17: 4562-4571.PubMedPubMed CentralView ArticleGoogle Scholar
- Maquat LE: Nonsense-mediated mRNA decay: splicing, translation and mRNP dynamics. Nat Rev Mol Cell Biol. 2004, 5: 89-99. 10.1038/nrm1310.PubMedView ArticleGoogle Scholar
- Reddy R, Henning D, Busch H: Pseudouridine residues in the 5'-terminus of uridine-rich nuclear RNA I (U1 RNA). Biochem Biophys Res Commun. 1981, 98: 1076-1083. 10.1016/0006-291X(81)91221-3.PubMedView ArticleGoogle Scholar
- Thomas J, Lea K, Zucker-Aprison E, Blumenthal T: The spliceosomal snRNAs of Caenorhabditis elegans. Nucleic Acids Res. 1990, 18: 2633-2642.PubMedPubMed CentralView ArticleGoogle Scholar
- Sorek R, Ast G: Intronic sequences flanking alternatively spliced exons are conserved between human and mouse. Genome Res. 2003, 13: 1631-1637. 10.1101/gr.1208803.PubMedPubMed CentralView ArticleGoogle Scholar
- Sugnet CW, Kent WJ, Ares M, Haussler D: Transcriptome and genome conservation of alternative splicing events in humans and mice. Pac Symp Biocomput. 2004, 66-77.Google Scholar
- Sorek R, Shemesh R, Cohen Y, Basechess O, Ast G, Shamir R: A non-EST-based method for exon-skipping prediction. Genome Res. 2004, 14: 1617-1623. 10.1101/gr.2572604.PubMedPubMed CentralView ArticleGoogle Scholar
- Baek D, Green P: Sequence conservation, relative isoform frequencies, and nonsense-mediated decay in evolutionarily conserved alternative splicing. Proc Natl Acad Sci USA. 2005, 102: 12813-12818. 10.1073/pnas.0506139102.PubMedPubMed CentralView ArticleGoogle Scholar
- Yeo GW, Van Nostrand E, Holste D, Poggio T, Burge CB: Identification and analysis of alternative splicing events conserved in human and mouse. Proc Natl Acad Sci USA. 2005, 102: 2850-2855. 10.1073/pnas.0409742102.PubMedPubMed CentralView ArticleGoogle Scholar
- Kazan K: Alternative splicing and proteome diversity in plants: the tip of the iceberg has just emerged. Trends Plant Sci. 2003, 8: 468-471. 10.1016/j.tplants.2003.09.001.PubMedView ArticleGoogle Scholar
- Iida K, Seki M, Sakurai T, Satou M, Akiyama K, Toyoda T, Konagaya A, Shinozaki K: Genome-wide analysis of alternative pre-mRNA splicing in Arabidopsis thaliana based on full-length cDNA sequences. Nucleic Acids Res. 2004, 32: 5096-5103. 10.1093/nar/gkh845.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang H, Blumenthal T: Functional analysis of an intron 3' splice site in Caenorhabditis elegans. RNA. 1996, 2: 380-388.PubMedPubMed CentralGoogle Scholar
- Brackenridge S, Wilkie AO, Screaton GR: Efficient use of a 'dead-end' GA 5' splice site in the human fibroblast growth factor receptor genes. EMBO J. 2003, 22: 1620-1631. 10.1093/emboj/cdg163.PubMedPubMed CentralView ArticleGoogle Scholar
- Crispino JD, Blencowe BJ, Sharp PA: Complementation by SR proteins of pre-mRNA splicing reactions depleted of U1 snRNP. Science. 1994, 265: 1866-1869.PubMedView ArticleGoogle Scholar
- Crispino JD, Sharp PA: A U6 snRNA:pre-mRNA interaction can be rate-limiting for U1-independent splicing. Genes Dev. 1995, 9: 2314-2323.PubMedView ArticleGoogle Scholar
- Chen JY, Stands L, Staley JP, Jackups RR, Latus LJ, Chang TH: Specific alterations of U1-C protein or U1 small nuclear RNA can eliminate the requirement of Prp28p, an essential DEAD box splicing factor. Mol Cell. 2001, 7: 227-232. 10.1016/S1097-2765(01)00170-8.PubMedView ArticleGoogle Scholar
- Hastings ML, Krainer AR: Pre-mRNA splicing in the new millennium. Curr Opin Cell Biol. 2001, 13: 302-309. 10.1016/S0955-0674(00)00212-X.PubMedView ArticleGoogle Scholar
- Automated Splice Site Analyses. [https://splice.cmh.edu/]
- Nalla VK, Rogan PK: Automated splicing mutation analysis by information theory. Hum Mutat. 2005, 25: 334-342. 10.1002/humu.20151.PubMedView ArticleGoogle Scholar
- Kopelman NM, Lancet D, Yanai I: Alternative splicing and gene duplication are inversely correlated evolutionary mechanisms. Nat Genet. 2005, 37: 588-589. 10.1038/ng1575.PubMedView ArticleGoogle Scholar
- Echard A, Opdam FJ, de Leeuw HJ, Jollivet F, Savelkoul P, Hendriks W, Voorberg J, Goud B, Fransen JA: Alternative splicing of the human Rab6A gene generates two close but functionally different isoforms. Mol Biol Cell. 2000, 11: 3819-3833.PubMedPubMed CentralView ArticleGoogle Scholar
- Chern TM, van Nimwegen E, Kai C, Kawai J, Carninci P, Hayashizaki Y, Zavolan M: A simple physical model predicts small exon length variations. PLoS Genet. 2006, 2: e45-10.1371/journal.pgen.0020045.PubMedPubMed CentralView ArticleGoogle Scholar
- Ray DW, Davis JR, White A, Clark AJ: Glucocorticoid receptor structure and function in glucocorticoid-resistant small cell lung carcinoma cells. Cancer Res. 1996, 56: 3276-3280.PubMedGoogle Scholar
- Rivers C, Levy A, Hancock J, Lightman S, Norman M: Insertion of an amino acid in the DNA-binding domain of the glucocorticoid receptor as a result of alternative splicing. J Clin Endocrinol Metab. 1999, 84: 4283-4286. 10.1210/jc.84.11.4283.PubMedView ArticleGoogle Scholar
- Hu CA, Lin WW, Obie C, Valle D: Molecular enzymology of mammalian Delta1-pyrroline-5-carboxylate synthase. Alternative splice donor utilization generates isoforms with different sensitivity to ornithine inhibition. J Biol Chem. 1999, 274: 6754-6762. 10.1074/jbc.274.10.6754.PubMedView ArticleGoogle Scholar
- Hinrichs AS, Karolchik D, Baertsch R, Barber GP, Bejerano G, Clawson H, Diekhans M, Furey TS, Harte RA, Hsu F, et al: The UCSC Genome Browser Database: update 2006. Nucleic Acids Res. 2006, D590-D598. 10.1093/nar/gkj144. 34 Database
- Splice-site Analyzer Tool. [http://ast.bioinfo.tau.ac.il/SpliceSiteFrame.htm]
- MaxEntScan Splice Site Scores. [http://genes.mit.edu/burgelab/maxent/Xmaxentscan_scoreseq.html]
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14: 1188-1190. 10.1101/gr.849004.PubMedPubMed CentralView ArticleGoogle Scholar
- TAndem Splice Site DataBase. [http://helios.informatik.uni-freiburg.de/TassDB/]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.