- Research article
- Open Access
A comprehensive survey of human polymorphisms at conserved splice dinucleotides and its evolutionary relationship with alternative splicing
© Shimada et al; licensee BioMed Central Ltd. 2010
- Received: 30 October 2009
- Accepted: 30 April 2010
- Published: 30 April 2010
Alternative splicing (AS) is a key molecular process that endows biological functions with diversity and complexity. Generally, functional redundancy leads to the generation of new functions through relaxation of selective pressure in evolution, as exemplified by duplicated genes. It is also known that alternatively spliced exons (ASEs) are subject to relaxed selective pressure. Within consensus sequences at the splice junctions, the most conserved sites are dinucleotides at both ends of introns (splice dinucleotides). However, a small number of single nucleotide polymorphisms (SNPs) occur at splice dinucleotides. An intriguing question relating to the evolution of AS diversity is whether mutations at splice dinucleotides are maintained as polymorphisms and produce diversity in splice patterns within the human population. We therefore surveyed validated SNPs in the database dbSNP located at splice dinucleotides of all human genes that are defined by the H-Invitational Database.
We found 212 validated SNPs at splice dinucleotides (sdSNPs); these were confirmed to be consistent with the GT-AG rule at either allele. Moreover, 53 of them were observed to neighbor ASEs (AE dinucleotides). No significant differences were observed between sdSNPs at AE dinucleotides and those at constitutive exons (CE dinucleotides) in SNP properties including average heterozygosity, SNP density, ratio of predicted alleles consistent with the GT-AG rule, and scores of splice sites formed with the predicted allele. We also found that the proportion of non-conserved exons was higher for exons with sdSNPs than for other exons.
sdSNPs are found at CE dinucleotides in addition to those at AE dinucleotides, suggesting two possibilities. First, sdSNPs at CE dinucleotides may be robust against sdSNPs because of unknown mechanisms. Second, similar to sdSNPs at AE dinucleotides, those at CE dinucleotides cause differences in AS patterns because of the arbitrariness in the classification of exons into alternative and constitutive type that varies according to the dataset. Taking into account the absence of differences in sdSNP properties between those at AE and CE dinucleotides, the increased proportion of non-conserved exons found in exons flanked by sdSNPs suggests the hypothesis that sdSNPs are maintained at the splice dinucleotides of newly generated exons at which negative selection pressure is relaxed.
- Splice Site
- Splice Pattern
- Flank Exon
- Negative Selection Pressure
Pre-mRNA splicing of eukaryotes requires three basic signals (splicing motifs) for the recognition of introns. The splicing motifs are the 5' intron end (donor) and the 3' intron end (acceptor), and the branch site. The splicing motifs at the 5' and 3' splice sites, known as "ag|GTragt" ("|" is the splice junction; "r" is a or g) and "(y)12-17nAG|g" ("y" is c or t; "n" is a, t, g or c; and subscript indicates the repeat number) [1, 2]. A human expressed sequence tag-based study showed that 99.24% and 0.69% of introns are flanked by GT-AG and GC-AG dinucleotides (splice dinucleotides), respectively . Other types of splice dinucleotides are also found in the human genome; these are AT-AC(0.05%) and others (0.02%) . Irrespectively of these variations at the splice dinucleotides, there are two well-studied splicing mechanisms . One mechanism utilizes a major spliceosome, an assembly of five small nuclear ribonucleoprotein particles (U1, U2, U4, U5, and U6 snRNP). The other mechanism uses the minor spliceosome, consisting of U11, U12, U4atac, U5, and U6atac, instead. Thus, most of exons are flanked by the virtually "invariant" GT and AG dinucleotides (splice dinucleotides) . Other additional splicing motifs, such as enhancers and silencers located in exons and introns have vast variety in motif signals and locations, but contribute to splicing fidelity.
Some single base-pair substitutions occurring at the "invariant" splice dinucleotides cause alteration of splice patterns and are associated with serious diseases, for example, NF1 , GSTM4 , cyclin D1 , NUDT1(MTH1) , and LDLR , (for review, see [10, 11]). The Human Genome Mutation Database (HGMD) at the Institute of Medical Genetics in Cardiff [12, 13] has annotated a total of 9267 entries for mutations in the vicinity of splice sites, which include 2362, 756, 1199, and 1355 entries for mutations at splice dinucleotides at sites '+1(G)', '+2(T)', '-2(A)', and '-1(G)', respectively . The databases DBASS5  and DBASS3  contain 431 and 283 details of aberrant splice sites, respectively , which are generated as a result of disease-causing mutations in humans.
While other single base-pair substitutions at splicing dinucleotides are known to be maintained as single nucleotide polymorphisms (SNPs) in human populations [17, 18], the question is: "what determines whether a single base-pair substitution at a splice dinucleotide will be maintained as a SNP or eliminated from the population?"
To address this question, we evaluated SNPs at splice dinucleotides (sdSNPs) in the context of selective pressure in the course of evolution. Generally, functional constraints on exons differ between alternatively spliced exons (ASEs) and constitutively spliced exons (CSEs). ASEs are subject to relaxed negative selective pressure, which is suggested by their significantly higher Ka/Ks values compared with other exons . This relaxation is the most fundamental conceptual constituent of exon creation via alternative splicing (AS) and was first proposed at the time of the discovery of the exon/intron structure in the 1970s . Recent studies revealed that AS is an important mechanism for creating new exons [21–24] and that accelerated accumulation of SNPs at additional splicing motifs after gene duplication enhances exon generation .
Current advances in genome informatics and comparative genomics demonstrate that ASEs can be sub-divided into two contrasting categories. When ASEs are classified as conserved or non-conserved in exon structure, low synonymous rates are characteristic of conserved ASEs but not of those with non-conserved exonic structure . Moreover, when they are classified as boundary-shifting (complex) ASEs or non-boundary shifting (simple) ASEs (those of the former type change the exon/intron boundaries of the flanking exons whereas those of the latter type do not), complex ASEs are under stronger selection pressure at the amino acid level but less pressure at the RNA level than CSEs, while reverse trends were observed in simple ASEs [27, 28]. These opposite evolutionary effects between different AS patterns have been discussed as a key role of AS in the 'switch-like' regulation of gene expression .
If it is supposed that sdSNPs are related to variation within populations in the regulation of gene expression through alterations in splicing patterns, the evolutionary profiles of sdSNP may differ between sdSNPs flanking ASEs (sdSNPs at AE dinucleotides) and those flanking CSEs (sdSNPs at CE dinucleotides).
Here, we extracted the sdSNPs from all human genes, and evaluated them by comparing flanking exon properties between ASEs or CSEs. Each group was subsequently divided into three subgroups according to exon conservation status between human and mouse, namely transcript-conserved, genome-conserved, and non-conserved (See Methods for details of the criteria used). We found that sdSNPs exist in the human genome with high allele frequencies, and with no significant difference in flanking exon properties between ASEs and CSEs. Moreover, we also found that these sdSNPs are prone to be maintained at the splice dinucleotides of newly generated exons. These results suggest that sdSNPs are associated with relaxed selection pressure for newly generated exons.
Validated SNPs at splice sites
SNP distribution and density at splice dinucleotides
Comparison of sdSNP properties between ASEs and CSEs
Properties of sdSNP
#HITs accordance with GT-AG rule3)
Splice score with matched allele4)
Splice score Δ two alleles5)
Conservation of exons flanked by SNPs
Allele prediction and matching to the GT-AG rule
Number of transcript-SNP pairs (number of SNPs)
Comparison of splice site scores
The observation that more predicted alleles followed the GT-AG rule at AE dinucleotides (35/47, 74%) than at CE dinucleotides (100/175, 57%) may suggest that recognition of CSE is robust against SNPs at their dinucleotides, while recognition of ASE is sensitive to SNPs (Table 1). We checked whether exon strength of ASEs was generally weaker than that of CSE, and whether this explains the observation.
We did not find any evidence to suggest that multiple sdSNPs interact with each other to maintain exon strength. First, exons flanked by multiple sdSNPs and introns with a sdSNP at both ends were not observed. Second, a splice dinucleotide located in two consecutive sdSNPs (rs3667 and rs1130151) was observed; however, these did not affect our analysis because the predicted haplotype did not follow the GT-AG rule, i.e., the AA haplotype was predicted at the 3' splice site.
The SNP dataset includes validated SNPs only (see Methods). The heterozygosities of these selected SNPs indicate that most of them are common SNPs (Additional file 1, Tables S1, Additional file 1, Table S2); the average heterozygosities of sdSNPs in AE and CE dinucleotides were 0.24 (S.E. 0.03) and 0.19 (S.E. 0.03), respectively (Table 1); the difference was not significant (Aspin-Welch's t = 1.15, df = 60, P > 0.05). The SNP validation and heterozygosity data therefore suggest that SNPs observed in sdSNPs are not mutations related to serious disease but polymorphisms that are maintained in the human population.
All comparisons of various properties of the sdSNPs, including average heterozygosity, SNP density, the ratio of predicted alleles satisfying the GT-AG rule, and splice scores using the predicted allele, show no significant differences between ASEs and CSEs (Table 1). This suggests that sdSNPs are maintained regardless of whether the flanking exon is of the ASE or CSE type.
We mapped 212 sdSNPs in the human genome by searching the total human transcriptome. It is worth noting that 159 (75%) of these sdSNPs mapped to CE dinucleotides, and 53 sdSNPs (25%) mapped to AE dinucleotides (Process 3 in Figure 2). The 53 sdSNPs mapped to AE dinucleotides suggests that these sdSNPs are potential candidates for AS polymorphisms, which are SNPs that alter the splicing pattern . In contrast to sdSNPs at AE dinucleotides, it is an intriguing discovery that sdSNPs exist at CE dinucleotides: this means that exons adjacent to these SNPs would be recognized constitutively regardless of allele. One explanation for the existence of sdSNPs at CE dinucleotides is that they represent an inconsistency between the definition of CSE based on publicly available human transcript data and on diversity in exon usage, which is affected by cell or tissue type, developmental stage, external stimuli [38, 39], and genetic variation [40–44]. The number of ASEs is increasing as detection methods develop [29, 45, 46]. Because we defined CSEs based on publicly available human transcript data, our CSE category may include some ASEs if defined according to deep sequencing technology. Moreover, the number of specimens used for the current full-length cDNA data is smaller than that for SNP surveys in the dbSNP database. This may also lead to a difference between our definition of ASEs and that after taking polymorphism into consideration. Thus, CSEs are defined arbitrarily and vary with the amount of transcript data, which means that the ASEs in our dataset are likely to be highly alternatively spliced but that the CSEs are not necessarily constitutive. Comparison between genome and transcript sequences using common numerous specimens may resolve this uncertainty.
Another explanation is that CSE is robust against SNPs at their splice dinucleotides because of an unknown mechanism that can work even if the GT-AG rule is partially violated. U1-independent splicing may be an example that supports this explanation [47, 48]. Although the detailed mechanism of U1-independent splicing is unclear and further research is needed, hF1γ intron 9, whose 5' splice site sequences are different from the consensus sequence at positions -3 and +5, spliced out without U1 snRNP. In this way, the 5' splice sites of these CSEs may be relaxed from selective pressure, which leads to maintenance of the sdSNPs at CE dinucleotides.
Our observation that the proportion of non-conserved exons was higher than that of conserved exons when the exons were flanked by sdSNPs suggests a relationship between exons with sdSNPs and newly generated exons (Figure 3). A question arises whether maintained sdSNPs lead to the generation of new exons, or whether newly generated exons bring about the maintenance of new flanking sdSNPs. Previous studies reported a lower degree of conservation of ASEs than CSEs and constructed models in which exons were generated as ASEs [21, 49, 50]. A previous study using the same datasets as the present study also showed that a limited proportion (6%) of AS variants was conserved between human and mouse . Taking account of these, it is unlikely that as many CSEs are generated as ASEs. Consequently, our finding of no significant difference between ASEs and CSEs regarding sdSNPs suggests only a low possibility that sdSNPs lead to the generation of new exons. Instead, newly generated exons can be considered to serve additional functions and to be subject to relaxed selective pressure, in accord with the relaxation that occurs with duplicated genes (i.e., evolution by gene duplication) . Our results suggest that newly generated exons allow the maintenance of flanking sdSNPs because of relaxed selective pressure on these exons.
We found 212 validated sdSNPs in the dbSNP database that were consistent with the GT-AG rule at either allele. Proportion of non-conserved exon was higher when the exons were flanked by sdSNPs. We found no significant difference in the properties of sdSNPs between the two types of flanking exon (ASEs and CSEs). This includes average heterozygosity, SNP density, the proportion of predicted alleles satisfying the GT-AG rule, and scores for splice sites formed with predicted alleles.
The sdSNPs flanking ASEs can be explained by SNPs that alter the splicing pattern by their alleles, while the existence of sdSNPs flanking CSEs requires an explanation for constitutive splicing even when an expressed allele violates the GT-AG rule. Although this may be explained by the arbitrariness of the definition of the CSE, which varies among dataset, it suggests that a novel splicing mechanism makes CSEs robust against sdSNPs, such as U1-independent splicing.
Taking account of previous studies suggesting an exon generation model through ASEs, there is little likelihood that sdSNPs provide as much opportunities to generate CSEs as ASEs. Our finding of similar properties of sdSNPs between ASEs and CSEs does not suggest the possibility that sdSNPs are exon generators. Instead, the increased proportion of non-conserved exons found in exons flanked by sdSNPs suggests the hypothesis that sdSNPs are maintained at the splice dinucleotides of newly generated exons whose negative selection pressures are relaxed.
Determination of exon/intron structure using H-InvDB
We used a comprehensive annotation resource for all human transcripts released by the H-Invitational Database 4.6 (H-InvDB_4.6 ), which includes gene structure, transcript mapping location, and splice isoforms . These data from H-InvDB_4.6 include all of the publicly available 120 558 human mRNAs mapped onto the human genome sequence (NCBI build 36.1) , forming 34 699 gene clusters, each of which we regarded as a putative "locus" (Figure 1A).
The H-Invitational consortium discriminates ASEs from CSEs using the H-Inv transcripts as previously described [32, 33]. The ASEs obtained include the following types: cassette, internal acceptor, internal donor, mutually excusive, and retained intron. H-Inv transcripts with the same gene structure were grouped, and a representative AS variant (RASV) was selected from each group [32, 54].
From all H-Inv transcripts, 80 997 were identified as AS variants and 41 289 were selected as RASVs ; these contained 107 606 ASEs and 298 102 CSEs (a total of 405 708 exons). We determined splice dinucleotides of all exons included in the filtered H-Inv transcripts using the genome position annotated in H-InvDB_4.6 (Figure 1A). We calculated the number of splice sites in the human genome by mapping all the RASV exons into the human genome.
Extraction of SNPs at splice dinucleotides from dbSNP build 127
We searched for SNPs located at splice dinucleotides using the downloaded data of dbSNP build 127 . The downloaded SNP data were mapped onto transcript positions (Figure 1A). We classified these retrieved SNPs located at splice dinucleotides into the following two groups according to the determination of ASEs and CSEs by H-DBAS: splice sites adjacent to ASEs; and those of CSEs, [32, 33]. We selected validated SNPs from dbSNP to avoid inserting sequencing errors into our dataset. Such errors are prone to occur at the G base of 3' splice dinucleotides because of suppression of G after incorporation of A . We removed the SNP-transcript pairs in which neither allele followed the GT-AG rule (Process 3 in Figure 2).
Calculation of SNP density
To calculate SNP densities in splice dinucleotides of ASEs and CSEs, we obtained the total length of each splice dinucleotide region in the genome by mapping splice dinucleotides of all RASVs given by H-DBAS [32, 33] onto the human reference genome sequence (NCBI build 36.1) . For the comparison with average SNP density in the intron region, we defined the total intron length as the region obtained by subtracting the total length of exons from all representative transcripts in H-InvDB_4.6. We excluded the region when an intron overlapped with an exon of another representative transcript. Thus, we counted validated SNPs within the total intron region and calculated the total length of the intron.
Allele prediction of SNPs at splice dinucleotides and inspection of agreement with the GT-AG rule
Since splice dinucleotides occur at both ends of introns, transcript sequences do not contain splice dinucleotides. To predict the SNP allele at splice dinucleotides for each transcript, we searched for another SNP that was most closely located at the exon neighboring the splice dinucleotides in every transcript. We extracted alleles of these two SNPs from published BAC clone sequences used for the human genome project . We presumed the combination of alleles of these two SNPs was a haplotype, assuming LD between them. We also presumed that the other combination of the other alleles at the two SNPs formed another haplotype. According to the allele of the exonic SNP on a transcript, we determined which haplotype corresponded to each of the transcripts (Figure 1B). In the case of Figure 1B, transcript 1 was estimated as the same haplotype as the BAC sequence because the sequence at the SNP1 position on the BAC sequence is A, and transcript 1 also has A at the same position. Consequently, because the SNP2 position of the BAC sequence has G, the SNP2 allele of the transcript 1 was estimated as G. Transcript 2 was estimated as a different haplotype from the BAC sequence because of the difference in sequence at the SNP1 position, which led to the prediction that the SNP2 allele is different from the haplotype of the BAC sequence. In this case there are four possible combinations of the two alleles (table in Figure 1B). Among these, however, two combinations (A at SNP1 and G at SNP2, T at SNP1 and C at SNP2) are expected to be haplotypes, assuming LD between the two SNPs.
For all transcripts, we examined whether the allele of each sdSNP predicted from the haplotype of the transcript agreed with the expectation from the GT-AG rule.
Maximum entropy scores
To estimate the strengths of 5' and 3' splice junctions in exons adjacent to SNPs at their splice dinucleotides, we calculated their splice site scores using the splice site models of MaxEntScan introduced by Yeo and Burge . MaxEntScan is based on an approach to modeling short sequences that takes into account adjacent and nonadjacent dependencies using large datasets of human splice sites, based on the principle of maximum entropy. These splice site models assign a log-odds ratio as a MaXEnt Score to 9-mer (5' splice site) and 23-mer (3' splice site) sequences. The higher the assigned score, the higher the probability that the sequence is a true splice site. We downloaded the Perl scripts of MaxEntScan via the Internet  and applied this method to both splice junction sequences taken from BAC clones used for the human reference genome and replaced SNP sites with another allele.
Inspection of exon mapping
Given two alleles of a SNP on a splice dinucleotide, in some cases there are two possible exon/intron boundaries that agree with the GT-AG rule (Figure 1C). One of these is the original boundary defined by H-InvDB ("Boundary-1 Allele-G" in Figure 1C), while the other is shifted one base pair ("Boundary-2 Allele-A" in Figure 1C). In this case, we selected the boundary that included the highest splice site score between two alleles' scores. We eliminated SNP-transcript pairs scoring negative values for both alleles when we summarized splice scores and compared them between ASEs and CSEs (Process 4 in Figure 2).
Classification of sdSNPs by conservation status between human and mouse genomes
For the sdSNPs that were classified according to the type of flanking exons as ASEs or CSEs (Process 2 in Figure 2), filtering out those for which both alleles did not satisfy the GT-AG rule (Process 3 in Figure 2), we subclassified according to the conservation status of flanking exons between human and mouse genomes (Processes 6 and 7 in Figure 2). We used exon classification as performed in a previous study  and published through H-DBAS [32, 33]; these references contain a detailed description of the method. Briefly, the criterion of the classification is as follows , human and mouse genome sequences that mapped with exons of their transcripts were aligned. If a human exon was conserved with respect to that of the mouse counterpart exon in the genome alignment with coverage and identity greater or equal in the threshold of 70% and 60%, respectively, it was defined as transcript-conserved. If a human exonic region aligned to the mouse genomic region exceeded these thresholds where no counterpart mouse exon was mapped, it was defined as genome-conserved. If a human exonic region in the genome did not align with the mouse genome using the thresholds stated above, the exon was defined as non-conserved. Because this conservation status was defined by genome position, the conservation status of a given sdSNP was determined uniquely even if multiple transcripts included the exon flanking the sdSNP (Processes 6 and 7 in Figure 2).
We thank Yumi Yamaguchi-Kabata and Ryuzo Matsumoto for their planning and advice on computer programming, and Craig Gough and Akila Mayeda for comments on an earlier version of this manuscript. We also thank all the other members of the Integrated Database and Systems Biology Team of BIRC, AIST and former members of the Integrated Database Team of JBIRC for technical assistance and for providing gene structure and annotation data. This research was financially supported by the Ministry of Economy, Trade and Industry of Japan (METI) and the Japan Biological Informatics Consortium (JBIC).
- Chasin LA: Searching for splicing motifs. Alternative Splicing in the Postgenomic Era. Edited by: Blencowe BJ, Graveley BR. 2007, New York: Landes Biosceince and Springer Science+Business Media, 85-106. full_text. [Back N, Cohen IR, Lajha A, Lambris JD, Paoletti R (Series Editors): Advances in Experimental Medicine and Biology, vol 623.]View ArticleGoogle Scholar
- Srebrow A, Kornblihtt AR: The connection between splicing and cancer. J Cell Sci. 2006, 119: 2635-2641. 10.1242/jcs.03053.View ArticlePubMedGoogle Scholar
- Burset M, Seledtsov IA, Solovyev VV: Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Res. 2000, 28: 4364-4375. 10.1093/nar/28.21.4364.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu Q, Krainer AR: AT-AC pre-mRNA splicing mechanisms and conservation of minor introns in voltage-gated ion channel genes. Mol Cell Biol. 1999, 19: 3225-3236.PubMed CentralView ArticlePubMedGoogle Scholar
- Ars E, Kruyer H, Morell M, Pros E, Serra E, Ravella A, Estivill X, Lazaro C: Recurrent mutations in the NF1 gene are common among neurofibromatosis type 1 patients. J Med Genet. 2003, 40: e82-10.1136/jmg.40.6.e82.PubMed CentralView ArticlePubMedGoogle Scholar
- Denson J, Xi Z, Wu Y, Yang W, Neale G, Zhang J: Screening for inter-individual splicing differences in human GSTM4 and the discovery of a single nucleotide substitution related to the tandem skipping of two exons. Gene. 2006, 379: 148-155. 10.1016/j.gene.2006.05.012.View ArticlePubMedGoogle Scholar
- Knudsen KE, Diehl JA, Haiman CA, Knudsen ES: Cyclin D1: polymorphism, aberrant splicing and cancer risk. Oncogene. 2006, 25: 1620-1628. 10.1038/sj.onc.1209371.View ArticlePubMedGoogle Scholar
- Oda H, Taketomi A, Maruyama R, Itoh R, Nishioka K, Yakushiji H, Suzuki T, Sekiguchi M, Nakabeppu Y: Multi-forms of human MTH1 polypeptides produced by alternative translation initiation and single nucleotide polymorphism. Nucleic Acids Res. 1999, 27: 4335-4343. 10.1093/nar/27.22.4335.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu H, Tucker HM, Grear KE, Simpson JF, Manning AK, Cupples LA, Estus S: A common polymorphism decreases low-density lipoprotein receptor exon 12 splicing efficiency and associates with increased cholesterol. Hum Mol Genet. 2007, 16: 1765-1772. 10.1093/hmg/ddm124.PubMed CentralView ArticlePubMedGoogle Scholar
- Cartegni L, Chew SL, Krainer AR: Listening to silence and understanding nonsense: exonic mutations that affect splicing. Nat Rev Genet. 2002, 3: 285-298. 10.1038/nrg775.View ArticlePubMedGoogle Scholar
- Garcia-Blanco MA, Baraniak AP, Lasda EL: Alternative splicing in disease and therapy. Nat Biotechnol. 2004, 22: 535-546. 10.1038/nbt964.View ArticlePubMedGoogle Scholar
- The Human Genome Mutation Database (HGMD) at the Institute of Medical Genetics in Cardiff. [http://www.hgmd.cf.ac.uk/ac/index.php]
- Stenson P, Mort M, Ball E, Howells K, Phillips A, Thomas N, Cooper D: The Human Gene Mutation Database: 2008 update. Genome Med. 2009, 1: 13-10.1186/gm13.PubMed CentralView ArticlePubMedGoogle Scholar
- The Human Gene Mutation Database (HGMD®): Statistics for Splicing Mutations. [http://www.hgmd.cf.ac.uk/ac/hoho2.php]
- Buratti E, Chivers M, Kralovicova J, Romano M, Baralle M, Krainer AR, Vorechovsky I: Aberrant 5' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization. Nucleic Acids Res. 2007, 35: 4250-4263. 10.1093/nar/gkm402.PubMed CentralView ArticlePubMedGoogle Scholar
- Vorechovsky I: Aberrant 3' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization. Nucleic Acids Res. 2006, 34: 4630-4641. 10.1093/nar/gkl535.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Platzer M, Hiller M, Szafranski K, Jahn N, Hampe J, Schreiber S, Backofen R, Huse K: Sequencing errors or SNPs at splice-acceptor guanines in dbSNP?. Nat Biotech. 2006, 24: 1068-1070. 10.1038/nbt0906-1068b.View ArticleGoogle Scholar
- Xing Y, Lee C: Evidence of functional selection pressure for alternative splicing events that accelerate evolution of protein subsequences. Proc Natl Acad Sci USA. 2005, 102: 13526-13531. 10.1073/pnas.0501213102.PubMed CentralView ArticlePubMedGoogle Scholar
- Gilbert W: Why genes in pieces?. Nature. 1978, 271: 501-10.1038/271501a0.View ArticlePubMedGoogle Scholar
- Modrek B, Lee CJ: Alternative splicing in the human, mouse and rat genomes is associated with an increased frequency of exon creation and/or loss. Nat Genet. 2003, 34: 177-180. 10.1038/ng1159.View ArticlePubMedGoogle Scholar
- Alekseyenko AV, Kim N, Lee CJ: Global analysis of exon creation versus loss and the role of alternative splicing in 17 vertebrate genomes. RNA. 2007, 13: 661-670. 10.1261/rna.325107.PubMed CentralView ArticlePubMedGoogle Scholar
- Sorek R: The birth of new exons: Mechanisms and evolutionary consequences. RNA. 2007, 13: 1603-1608. 10.1261/rna.682507.PubMed CentralView ArticlePubMedGoogle Scholar
- Corvelo A, Eyras E: Exon creation and establishment in human genes. Genome Biol. 2008, 9: R141-10.1186/gb-2008-9-9-r141.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Z, Zhou L, Wang P, Liu Y, Chen X, Hu L, Kong X: Divergence of exonic splicing elements after gene duplication and the impact on gene structures. Genome Biol. 2009, 10: R120-10.1186/gb-2009-10-11-r120.PubMed CentralView ArticlePubMedGoogle Scholar
- Plass M, Eyras E: Differentiated evolutionary rates in alternative exons and the implications for splicing regulation. BMC Evol Biol. 2006, 6: 50-10.1186/1471-2148-6-50.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen FC, Chaw SM, Tzeng YH, Wang SS, Chuang TJ: Opposite evolutionary effects between different alternative splicing patterns. Mol Biol Evol. 2007, 24: 1443-1446. 10.1093/molbev/msm072.View ArticlePubMedGoogle Scholar
- Chen FC, Chuang TJ: Different alternative splicing patterns are subject to opposite selection pressure for protein reading frame preservation. BMC Evol Biol. 2007, 7: 179-10.1186/1471-2148-7-179.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509.PubMed CentralView ArticlePubMedGoogle Scholar
- H-Invitational Database (H-InvDB) an integrated database of human genes and transcripts. [http://hinv.jp/]
- Yamasaki C, Murakami K, Fujii Y, Sato Y, Harada E, Takeda J-i, Taniya T, Sakate R, Kikugawa S, Shimada M, Tanino M, Koyanagi KO, Barrero RA, Gough C, Chun HW, Habara T, Hanaoka H, Hayakawa Y, Hilton PB, Kaneko Y, Kanno M, Kawahara Y, Kawamura T, Matsuya A, Nagata N, Nishikata K, Noda AO, Nurimoto S, Saichi N, Sakai H, Sanbonmatsu R, Shiba R, Suzuki M, Takabayashi K, Takahashi A, Tamura T, Tanaka M, Tanaka S, Todokoro F, Yamaguchi K, Yamamoto N, Okido T, Mashima J, Hashizume A, Jin L, Lee KB, Lin YC, Nozaki A, Sakai K, Tada M, Miyazaki S, Makino T, Ohyanagi H, Osato N, Tanaka N, Suzuki Y, Ikeo K, Saitou N, Sugawara H, O'Donovan C, Kulikova T, Whitfield E, Halligan B, Shimoyama M, Twigger S, Yura K, Kimura K, Yasuda T, Nishikawa T, Akiyama Y, Motono C, Mukai Y, Nagasaki H, Suwa M, Horton P, Kikuno R, Ohara O, Lancet D, Eveno E, Graudens E, Imbeaud S, Debily MA, Hayashizaki Y, Amid C, Han M, Osanger A, Endo T, Thomas MA, Hirakawa M, Makalowski W, Nakao M, Kim NS, Yoo HS, De Souza SJ, Bonaldo Mde F, Niimura Y, Kuryshev V, Schupp I, Wiemann S, Bellgard M, Shionyu M, Jia L, Thierry-Mieg D, Thierry-Mieg J, Wagner L, Zhang Q, Go M, Minoshima S, Ohtsubo M, Hanada K, Tonellato P, Isogai T, Zhang J, Lenhard B, Kim S, Chen Z, Hinz U, Estreicher A, Nakai K, Makalowska I, Hide W, Tiffin N, Wilming L, Chakraborty R, Soares MB, Chiusano ML, Auffray C, Yamaguchi-Kabata Y, Itoh T, Hishiki T, Fukuchi S, Nishikawa K, Sugano S, Nomura N, Tateno Y, Imanishi T, Gojobori T: The H-Invitational Database (H-InvDB), a comprehensive annotation resource for human genes and transcripts. Nucleic Acids Res. 2008, 36: D793-799.PubMedGoogle Scholar
- Takeda J-i, Suzuki Y, Nakao M, Barrero RA, Koyanagi KO, Jin L, Motono C, Hata H, Isogai T, Nagai K, Otsuki T, Kuryshev V, Shionyu M, Yura K, Go M, Thierry-Mieg J, Thierry-Mieg D, Wiemann S, Nomura N, Sugano S, Gojobori T, Imanishi T: Large-scale identification and characterization of alternative splicing variants of human gene transcripts using 56 419 completely sequenced and manually annotated full-length cDNAs. Nucleic Acids Res. 2006, 34: 3917-3928. 10.1093/nar/gkl507.PubMed CentralView ArticlePubMedGoogle Scholar
- Takeda J-i, Suzuki Y, Nakao M, Kuroda T, Sugano S, Gojobori T, Imanishi T: H-DBAS: Alternative splicing database of completely sequenced and manually annotated full-length cDNAs based on H-Invitational. Nucleic Acids Res. 2007, 35: D104-109. 10.1093/nar/gkl854.PubMed CentralView ArticlePubMedGoogle Scholar
- Human-transcriptome DataBase for Alternative Splicing (H-DBAS). [http://h-invitational.jp/h-dbas/]
- VarySysDB a human genetic polymorphism database based on all H-InvDB transcripts. [http://h-invitational.jp/varygene/home.htm]
- Shimada MK, Matsumoto R, Hayakawa Y, Sanbonmatsu R, Gough C, Yamaguchi-Kabata Y, Yamasaki C, Imanishi T, Gojobori T: VarySysDB: a human genetic polymorphism database based on all H-InvDB transcripts. Nucleic Acids Res. 2009, 37: D810-815. 10.1093/nar/gkn798.PubMed CentralView ArticlePubMedGoogle Scholar
- Nembaware V, Wolfe KH, Bettoni F, Kelso J, Seoighe C: Allele-specific transcript isoforms in human. FEBS Letters. 2004, 577: 233-238. 10.1016/j.febslet.2004.10.018.View ArticlePubMedGoogle Scholar
- Hiller M: Bioinformatics analyses of alternative splicing: Non-EST based prediction, influence of secondary structures and tandem splice sites. 2008, Saarbrücken: VDM Verlag Dr. MullerGoogle Scholar
- Stamm S: Signals and their transduction pathways regulating alternative splicing: a new dimension of the human genome. Hum Mol Genet. 2002, 11: 2409-2416. 10.1093/hmg/11.20.2409.View ArticlePubMedGoogle Scholar
- Hull J, Campino S, Rowlands K, Chan MS, Copley RR, Taylor MS, Rockett K, Elvidge G, Keating B, Knight J, Kwiatkowski D: Identification of common genetic variation that modulates alternative splicing. PLoS Genet. 2007, 3: e99-10.1371/journal.pgen.0030099.PubMed CentralView ArticlePubMedGoogle Scholar
- Kwan T, Benovoy D, Dias C, Gurd S, Provencher C, Beaulieu P, Hudson TJ, Sladek R, Majewski J: Genome-wide analysis of transcript isoform variation in humans. Nat Genet. 2008, 40: 225-231. 10.1038/ng.2007.57.View ArticlePubMedGoogle Scholar
- Coulombe-Huntington J, Lam KCL, Dias C, Majewski J: Fine-scale variation and genetic determinants of alternative splicing across individuals. PLoS Genetics. 2009, 5: e1000766-10.1371/journal.pgen.1000766.PubMed CentralView ArticlePubMedGoogle Scholar
- ElSharawy A, Hundrieser B, Brosch M, Wittig M, Huse K, Platzer M, Becker A, Simon M, Rosenstiel P, Schreiber S, Krawczak M, Hampe J: Systematic evaluation of the effect of common SNPs on pre-mRNA splicing. Human Mutation. 2009, 30: 625-632. 10.1002/humu.20906.View ArticlePubMedGoogle Scholar
- Nembaware V, Lupindo B, Schouest K, Spillane C, Scheffler K, Seoighe C: Genome-wide survey of allele-specific splicing in humans. BMC Genomics. 2008, 9: 265-10.1186/1471-2164-9-265.PubMed CentralView ArticlePubMedGoogle Scholar
- Blencowe BJ, Ahmad S, Lee LJ: Current-generation high-throughput sequencing: deepening insights into mammalian transcriptomes. Genes Dev. 2009, 23: 1379-1386. 10.1101/gad.1788009.View ArticlePubMedGoogle Scholar
- Kim N, Lee C: Bioinformatics detection of alternative splicing. Bioinformatics: Data, Sequence Analysis and Evolution. 2008, Totowa, NJ; New York: Humana Press, 179-197. [Keith JM (Series Editor): Methods in Molecular Biology, vol 452.]View ArticleGoogle Scholar
- Fukumura K, Taniguchi I, Sakamoto H, Ohno M, Inoue K: U1-independent pre-mRNA splicing contributes to the regulation of alternative splicing. Nucleic Acids Res. 2009, 37: 1907-1914. 10.1093/nar/gkp050.PubMed CentralView ArticlePubMedGoogle Scholar
- Fukumura K, Inoue K: Role and mechanism of U1-independent pre-mRNA splicing in the regulation of alternative splicing. RNA Biol. 2009, 6: 395-398. 10.4161/rna.6.4.9318.View ArticlePubMedGoogle Scholar
- Cusack BP, Wolfe KH: Changes in alternative splicing of human and mouse genes are accompanied by faster evolution of constitutive exons. Mol Biol Evol. 2005, 22: 2198-2208. 10.1093/molbev/msi218.View ArticlePubMedGoogle Scholar
- Wang W, Zheng H, Yang S, Yu H, Li J, Jiang H, Su J, Yang L, Zhang J, McDermott J, Samudrala R, Wang J, Yang H, Yu J, Kristiansen K, Wong GK: Origin and evolution of new exons in rodents. Genome Res. 2005, 15: 1258-1264. 10.1101/gr.3929705.PubMed CentralView ArticlePubMedGoogle Scholar
- Takeda J-i, Suzuki Y, Sakate R, Sato Y, Seki M, Irie T, Takeuchi N, Ueda T, Nakao M, Sugano S, Gojobori T, Imanishi T: Low conservation and species-specific evolution of alternative splicing in humans and mice: comparative genomics analysis using well-annotated full-length cDNAs. Nucleic Acids Res. 2008, 36: 6386-6395. 10.1093/nar/gkn677.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohno S, Wolf U, Atkin NB: Evolution from fish to mammals by gene duplication. Hereditas. 1968, 59: 169-187. 10.1111/j.1601-5223.1968.tb02169.x.View ArticlePubMedGoogle Scholar
- Entrez Genome (NCBI). [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/sites/entrez?db=genome]
- Kimura K, Wakamatsu A, Suzuki Y, Ota T, Nishikawa T, Yamashita R, Yamamoto Ji, Sekine M, Tsuritani K, Wakaguri H, Ishii S, Sugiyama T, Saito K, Isono Y, Irie R, Kushida N, Yoneyama T, Otsuka R, Kanda K, Yokoi T, Kondo H, Wagatsuma M, Murakawa K, Ishida S, Ishibashi T, Takahashi-Fujii A, Tanase T, Nagai K, Kikuchi H, Nakai K, Isogai T, Sugano S: Diversification of transcriptional modulation: Large-scale identification and characterization of putative alternative promoters of human genes. Genome Res. 2006, 16: 55-65. 10.1101/gr.4039406.PubMed CentralView ArticlePubMedGoogle Scholar
- dbSNP NCBI database of Single Nucleotide Polymorphism. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/projects/SNP/]
- 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.View ArticlePubMedGoogle Scholar
- A simple perl scripts to run MaxEntScan. [http://genes.mit.edu/burgelab/maxent/download/]
- Scheme of genomic comparison between RASVs and mouse cDNAs (H-Inv full-length cDNAs only). [http://h-invitational.jp/h-dbas/document/scheme_genomic_comparison.html]
- Eng L, Coutinho G, Nahas S, Yeo G, Tanouye R, Babaei M, Dök T, Burge C, Gatti RA: Nonclassical splicing mutations in the coding and noncoding regions of the ATM gene: Maximum entropy estimates of splice junction strengths. Human Mutation. 2004, 23: 67-76. 10.1002/humu.10295.View ArticlePubMedGoogle Scholar
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.