Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

Fast rate of evolution in alternatively spliced coding regions of mammalian genes

Ekaterina O Ermakova12*, Ramil N Nurtdinov1 and Mikhail S Gelfand12

Author Affiliations

1 Department of Bioengineering and Bioinformatics, Moscow State University, Vorob'evy gory, 1-73, 119992, Moscow, Russia

2 Research and Training Center "Bioinformatics", Institute for Information Transmission Problems, Russian Academy of Sciences, Bolshoi Karetny per. 19, 127994, Moscow, Russia

For all author emails, please log on.

BMC Genomics 2006, 7:84  doi:10.1186/1471-2164-7-84


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/7/84


Received:21 December 2005
Accepted:18 April 2006
Published:18 April 2006

© 2006 Ermakova et al; licensee 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.

Abstract

Background

At least half of mammalian genes are alternatively spliced. Alternative isoforms are often genome-specific and it has been suggested that alternative splicing is one of the major mechanisms for generating protein diversity in the course of evolution. Another way of looking at alternative splicing is to consider sequence evolution of constitutive and alternative regions of protein-coding genes. Indeed, it turns out that constitutive and alternative regions evolve in different ways.

Results

A set of 3029 orthologous pairs of human and mouse alternatively spliced genes was considered. The rate of nonsynonymous substitutions (dN), the rate of synonymous substitutions (dS), and their ratio (ω = dN/dS) appear to be significantly higher in alternatively spliced coding regions compared to constitutive regions. When N-terminal, internal and C-terminal alternatives are analysed separately, C-terminal alternatives appear to make the main contribution to the observed difference. The effects become even more pronounced in a subset of fast evolving genes.

Conclusion

These results provide evidence of weaker purifying selection and/or stronger positive selection in alternative regions and thus one more confirmation of accelerated evolution in alternative regions. This study corroborates the theory that alternative splicing serves as a testing ground for molecular evolution.

Background

Alternative splicing is a major mechanism for generating functional and evolutionary diversity of proteins in mammals [1,2], for a review see [3]. Indeed, alternative splicing allows for generation of novel proteins without sacrificing old ones [2]. If a new isoform proves to be beneficial, its fraction increases by subtle regulatory changes. On the other hand, unlike gene duplication, alternative splicing does not lead to dramatic changes in protein concentrations. Moreover, it has been demonstrated that duplicated genes are rarely alternatively spliced compared to singletons [4,5].

There are good reasons to believe that some key mutational events driving evolution might reside in introns, untranslated regions (UTRs) and/or nontranscribed regulatory regions [6-8]. A large fraction of alternative splicing events occur in untranslated regions [9]. Nevertheless, most studies of molecular evolution have focused on the analysis of protein coding regions, as these data are simpler to obtain and are more amenable to functional interpretation.

From this point of view, alternative regions of genes occupy an intermediate position. Alternatively spliced regions are often evolutionary young: indeed, about a half of genes in human-mouse orthologous pairs have species-specific isoforms [2,10]. In many respects, constitutive and alternative regions are organized in different ways. Alternative human splice sites are on the average weaker than constitutive ones [11]. Non-canonical GC-AG introns tend to be alternative [12]. Among human exons conserved in mouse, about 77% of alternative cassette exons are flanked on both sides by long conserved intronic sequences, compared to only 17% of the constitutive exons [13]. Overall, statistical and evolutionary features of constitutive and alternative exons are sufficiently different to provide for computational recognition of these exons [14-16].

In several early studies it has been observed that patterns of nucleotide substitutions are different in alternative and constitutive coding regions. Iida and Akashi [17] analysed 26 pairs of alternatively spliced human genes and their non-human mammalian orthologs and demonstrated that synonymous divergence was lower and the nonsynonymous divergence was higher in alternative regions compared to constitutive regions. Evidence for diversifying selection was observed in alternative regions of CD45 [18], whereas the reduced rate of synonymous substitutions in an alternative region of BRCA1 [19] was assigned to purifying selection due to exonic splicing enhancer sites [20]. Recently, lower synonymous divergence in alternative exons compared to constitutive ones was demonstrated in a large-scale study of human, chimpanzee, mouse and rat genes [21].

Here we analyze evolutionary patterns in a set of 3029 pairs of orthologous human and mouse genes. We consider all types of alternative splicing and analyze separately 5'-, internal, and 3'-regions of genes, as well as faster and slower evolving genes.

Results

We considered 3029 alternatively spliced human genes and their mouse orthlogs (Figure 1). The sample was divided into three bins of equal size with respect to nucleotide identity in coding regions. Nucleotide alignments of coding regions were sliced into constitutive (C) and alternative (A) fragments. Alternative fragments were further sorted into N-terminal (AN), internal (AI), and C-terminal (AC). To reveal the general pattern of evolution in these regions, we estimated the amino acid identity (Id), the nonsynonymous substitution rate (dN), the synonymous substitution rate (dS), and ω = dN/dS (Table 1). Global meta-alignments (concatenated alignment fragments across all considered genes and across three rate bins, see Methods) of five types (C, A, AN, AI, AC) were used. We performed 2000 bootstrap replications to evaluate the robustness of our estimates (Figures 2, 3, 4, 5).

thumbnailFigure 1. The data flow through the analysis pipeline.

thumbnailFigure 2. Amino acid identity in alignments of constitutive and alternative coding regions of 3029 alternatively spliced genes classified by the overall evolution rate. C – constitutive regions. A – alternative regions. AN – N-terminal alternative regions. AI – internal alternative regions. AC – C-terminal alternative regions. Five columns on the left show results for three equal bins of slowly evolving genes, genes evolving with medium speed, and fast evolving genes. Five columns on the right correspond to the total sample of 3029 alternatively spliced genes. To estimate the average amino acid identity, we used global meta-alignments (see Methods). The boxplots represent the results for 2000 bootstrap replications.

thumbnailFigure 3. Nonsynonymous substitution rate (dN) in alignments of constitutive and alternative coding regions of 3029 alternatively spliced genes genes classified by the overall evolution rate. Notation and layout as in Figure 2.

thumbnailFigure 4. Synonymous substitution rate (dS) in alignments of constitutive and alternative coding regions of 3029 alternatively spliced genes classified by the overall evolution rate. Notation and layout as in Figure 2.

thumbnailFigure 5. Selection measure ω = dN/dS in alignments of constitutive and alternative coding regions of 3029 alternatively spliced genes classified by the overall evolution rate. Notation and layout as in Figure 2.

Table 1. Substitution rates in constitutive and alternative regions of human-mouse concatenated alignments

It turned out that Id(A)<Id(C), dN(A)>dN(C), and ω(A)>ω(C) for alternatively spliced genes irrespective of the rate of evolution. These results show that negative selection is weaker and/or positive selection is stronger in alternative regions and thus confirm that alternatively spliced coding regions are hotspots of molecular evolution. Unexpectedly, dN and ω rise dramatically at the C-terminal alternative regions (Figures 3 and 5).

The pattern of synonymous substitutions is more complicated, as it depends on the rate of evolution (Figure 4). The general pattern is that dS in alternative regions increases in the 5' to 3' direction. In genes evolving at the medium rate, dS(AN)<dS(AI)<dS(AC), whereas in fast evolving genes dS(AI)>dS(AC).

For control, we considered N-terminal and C-terminal constitutively spliced regions and performed similar analysis. All computed evolutionary parameters were the same as in the constitutive regions in general (data not shown). Thus the observed difference cannot be explained simply by faster evolution at gene termini.

Next, we analyzed individual gene pairs. For each of 2358 genes with the total lengths of constitutive and of alternative regions both exceeding 80 bp, ω for constitutive (ωC) and alternative (ωA) regions was calculated separately. Figure 6a represents the distribution of the difference (ωCA). The distribution is skewed, showing that ω tends to be greater in alternative regions. We used the chi-squared test to compare the distributions of |ωCA| in the case ωCA and in the case ωCA. The null hypothesis that the distributions were the same was rejected at the significance level 10-7. When N-terminal, internal, and C-terminal alternatives were considered separately, the effect was the same (Figure 6b–d). The null hypothesis that the distribution of |ωCAN| was symmetrical was rejected at the significance level 10-9, the one that the distribution of |ωCAI| was symmetrical, at the significance level 10-3, and the one that the distribution of |ωCAC| was symmetrical, at the significance level 10-2. Therefore, the detailed analysis of individual genes confirmed the observations made on concatenated alignments.

thumbnailFigure 6. Distributions of ωCA, ωCAN, ωCAI, ωCAC. a. The distribution of ωCA for 2358 genes with the total length of human-mouse alignments of constitutive regions and of alternative regions both exceeding 80 bp shows that ω tends to be larger in alternative regions. In particular, there are 23 genes with ωCA<-0.8 and no genes with ωCA>0.8. b, c, d. The distributions of ωCAN, ωCAI, ωCAC for genes with long N-terminal (1674 genes, top, ωCAN), internal (976 genes, middle, ωCAI), and C-terminal (110 genes, bottom, ωCAC) alternative regions, respectively, show that ωA* tends to be larger than ωc in all types of alternative regions. the grey shadows are symmetrical to the left parts of the histograms.

Discussion

Evolutionary patterns in different functional regions are known to be significantly different. Conserved genes are duplicated relatively more often [22], although shortly after duplication the evolutionary rate might increase, as the purifying selection is weaker [23,24], and the selection pattern in the two copies may be different [25]. Duret and Mouchiroud [26] observed lower nonsynonymous divergence in genes expressed in multiple tissues when compared to genes with more limited expression patterns, whereas the synonymous substitution rate was roughly the same. Similarly, Pál, Papp and Hurst [27] demonstrated that highly expressed genes tend to be more conserved then genes expressed at a lower level. Our results are consistent with these observations if one assumes that constitutive regions are expressed in more tissues and at a higher rate that alternative ones. Indeed, the former assumption holds for genes with isoforms having clear tissue specific expression pattern, whereas the latter holds for genes with all isoforms expressed ubiquitously.

Young gene regions tend to evolve faster. Several studies [23-25,28] demonstrated post-duplicational relaxation of purifying selection in paralogs. Our results provided evidence of stronger positive selection and/or weaker purifying selection in alternative gene regions.

One possible explanation for our observations could be that the data are contaminated by non-functional isoforms (hence relaxation of purifying selection takes place). We do not believe that to be the case for the following reasons. Firstly, these regions were conserved between human and mouse at a sufficiently high similarity level of 70% nucleotide identity. Secondly, the observed pattern of increased dN level in alternative regions was the most pronounced in 3' (C-terminal) regions, that are the most reliable as regards gene recognition and have higher EST coverage due to polyA-primed ESTs.

As we considered only alternatives derived from RefSeq proteins, we could miss some alternatives and thus label a fraction of the alternative regions as constitutive. However, that could only contaminate the constitutive sample with alternative regions and thus blur the observed differences, but not create any spurious effect.

Recently Xing and Lee [21] observed similar rate of non-synonymous substitutions in alternative and constitutive regions whereas the rate of synonymous substitutions was lower in alternative regions, especially in tissue-specific exons [29]. One possible explanation for that was based on the assumption that conserved alternative exons contain more candidate splicing enhancer sites than constitutive ones [15]. As such sites could be expected to be conserved, like in BRCA1 [19,20], that could lead to higher conservation of synonymous codon positions in alternative regions compared to constitutive ones. However, this explanation seems to be incorrect, since, although indeed dS is lower in splicing enhancers, the fractions of constitutive and alternative regions covered by splicing enhancers are the same [30], and if the RNA selection pressure is the same in alternative and constitutive regions, it cannot distort the measurement of ω [31].

However, this effect has not been observed in our study, and the substitution rates differ from those in [21]. As our results are consistent and statistically significant for all classes of genes (fast, medium and slow evolving) and all gene regions (N-terminal, internal, C-terminal), and do not seem to be caused by contamination, there should be other reasons for this discrepancy. One of them could be the fact that we considered all types of alternatives, as opposed to only cassette exons in other studies. We also considered short alternative regions, its skipping or inclusion might be regulated "outside". Another one could be the use of different methods to calculate the rates of evolution. We used our own implementation of the first method of Ina [32] here, as we needed a tool for very long alignments (~ 3·106 bp), whereas Xing and Lee [21,29] used a maximum likelihood method implemented in the PAML package [33]. On the other hand, we considered only RefSeq isoforms and did not distinguish between the minor and major isoform alternatives.

An explanation for our finding could be that the total length of regulatory sites experiencing purifying selection is still small compared to the total length of alternative regions. The pattern of substitutions in insects is less consistent [34]: in N-terminal alternatives, the synonymous rate is higher than in constant regions, whereas in internal alternatives, there are more amino acid substitutions, similar to our observations here.

Conclusion

Overall, this study corroborates the idea that alternative splicing serves as a testing ground for molecular evolution. Several lines of evidence confirm this hypothesis: (i) alternatively spliced isoforms are often evolutionary young both in mammals [2,10] and in insects [35]; (ii) the rate of nonsynonymous substitutions is higher in alternative regions compared to constitutive ones (this study), (iii) constitutive exons in genes with genome-specific alternative splicing evolve faster than constitutive regions in genes with conserved structure [36] (cf. a similar observation for duplicated genes [23-25]), (iv) many young (rodent-specific, missing in human and pig as an outgroup) exons are alternatively spliced and tend to have ω>1 in the mouse-rat comparison [30], and (v) the frequency of nonsynonymous SNPs in human genes is higher in alternative regions than in constitutive regions [37].

Methods

Definitions

In an alternatively spliced gene, constitutive regions are defined as the ones that are always exonic and coding, and alternative regions as the ones that are either coding or spliced out. An exon can be either completely constitutive, or completely alternative, or non-coding, or consist of constitutive, alternative and non-coding regions.

A local meta-alignment is a concatenate of all alignment fragments of a fixed type (for example, coding alternative regions) for one particular gene.

A global meta-alignment is the concatenate of local meta-alignments of a fixed type for all genes of a fixed group (for example, for all fast-evolving genes).

Data

Human and mouse mRNA sequences were taken from the NCBI RefSeq database [38] and orthologs were identified and aligned as described previously [39].

Overall, 12356 pairs of orthologous human and mouse genes were considered. The data flow through the analysis pipeline is shown in Figure 1. Out of 12356 human genes, 5754 genes had more than one protein isoform in the EDAS database of alternatively spliced genes [40]. These proteins were mapped to the human-mouse mRNA alignments using Pro-Frame [41] and the results were parsed using a set of Perl scripts that identified constitutive and alternative coding fragments of the human genes and their reading frames. Alternatives confirmed by protein sequences and read in a single frame were identified in 3079 genes. We further restricted the dataset to 3029 pairs with more than 70% nucleotide identity between the human and mouse genes, as we doubted that the other ones were reliable.

2358 genes were selected for individual substitution rate analysis. These were the ones with both the total length of the human-mouse alignment length of the alternative regions and that of the constitutive regions exceeding 80 base pairs.

Data classification

We grouped genes with comparable average substitution rates and formed three bins of equal size: slow, medium-speed, and fast-evolving genes.

We also considered alternative coding regions corresponding to protein N-terminal, middle, and C-terminal parts separately.

Estimation of substitution rates

The transitional to transversional substitution rate ratio R, as well as thenumbers of synonymous (dS) and nonsynonymous (dN) substitutions per site were estimated by the Ina method I [32]. Unlike maximum likelihood methods, it afforded considerable results for very long alignments (~ 3·106 bp) and it proved to be fast enough to allow bootstrap resampling. We used our own implementation of this method (a set of Perl scripts).

Bootstrapping

To evaluate the robustness of the estimates for evolutionary parameters of the global meta-alignments, we used bootstrapping to form 2000 alighments of the same length for reach global meta-alignment and estimated amino-acid identity, dN, dS, and ω = dN/dS.

Authors' contributions

RN provided the EDAS database. EE and RN analysed the data. MG designed the project. EE and MG wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

We are grateful to I. King Jordan for alignments of human and mouse mRNAs and to Georgii Bazykin, Alexey Kondrashov, Dmitry Malko, Andrei Mironov, Dmitri Petrov, and Vasily Ramensky for useful discussions.

This study was partially supported by grants from the Ludwig Institute of Cancer Research (CDRF RBO-1268), the Howard Hughes Medical Institute (55000309), the Russian Fund of Basic Research (04-04-49440), the Russian Academy of Sciences (program "Molecular and Cellular Biology"), and the Russian Science Support Fund.

References

  1. 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. PubMed Abstract | Publisher Full Text OpenURL

  2. 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. PubMed Abstract | Publisher Full Text OpenURL

  3. Gelfand MS: Computational analysis of alternative splicing. In Handbook of Computational Molecular Biology, Chapman & Hall/CRC Computer & Information Science Series. Volume 9. Edited by Alluru S. New York: Chapman & Hall/CRC; 2005. OpenURL

  4. Kopelman NM, Lancet D, Yanai I: Alternative splicing and gene duplication are inversely correlated evolutionary mechanisms.

    Nature Genet 2005, 37:588-589. PubMed Abstract | Publisher Full Text OpenURL

  5. Su Z, Wang J, Yu J, Huang X, Gu X: Evolution of alternative splicing after gene duplication.

    Genome Res 2006, 16:182-189. PubMed Abstract | Publisher Full Text OpenURL

  6. Mazumder B, Seshadri V, Fox PL: Translational control by the 3'-UTR: the ends specify the means.

    Trends Biochem Sci 2003, 28:91-98. PubMed Abstract | Publisher Full Text OpenURL

  7. Pagani F, Baralle FE: Genomic variants in exons and introns: identifying the splicing spoilers.

    Nat Rev Genet 2004, 5:389-396. PubMed Abstract | Publisher Full Text OpenURL

  8. Lynch M, Scofield DG, Hong X: The evolution of transcription-initiation sites.

    Mol Biol Evol 2005, 22:1137-1146. PubMed Abstract | Publisher Full Text OpenURL

  9. Mironov AA, Fickett JW, Gelfand MS: Frequent alternative splicing of human genes.

    Genome Res 1999, 9:1288-1293. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Nurtdinov RN, Artamonova II, Mironov AA, Gelfand MS: Low conservation of alternative splicing patterns in the human and mouse genomes.

    Hum Mol Genet 2003, 12:1313-1320. PubMed Abstract | Publisher Full Text OpenURL

  11. Clark F, Thanaraj TA: Categorization and characterization of transcript-confirmed constitutively and alternatively spliced introns and exons from human.

    Hum Mol Genet 2001, 11:451-464. Publisher Full Text OpenURL

  12. Burset M, Seledtsov IA, Solovyev VV: Analysis of canonical and non-canonical splice sites in mammalian genomes.

    Nucleic Acids Res 2000, 28:4364-4375. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Sorek R, Ast G: Intronic sequences flanking alternatively spliced exons are conserved between human and mouse.

    Genome Res 2003, 13:1631-1637. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. 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 U S A 2005, 102:2850-2855. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Xing Y, Lee C: Assessing the application of Ka/Ks ratio test to alternatively spliced exons.

    Bioinformatics 2005, 21:3701-3797. PubMed Abstract | Publisher Full Text OpenURL

  17. Iida K, Akashi H: A test of translational selection at 'silent' sites in the human genome: base composition comparisons in alternatively spliced genes.

    Gene 2000, 261:93-105. PubMed Abstract | Publisher Full Text OpenURL

  18. Filip LC, Mundy NI: Rapid evolution by positive Darwinian selection in the extracellular domain of the abundant lymphocyte protein CD45 in primates.

    Mol Biol Evol 2004, 21:1504-1511. PubMed Abstract | Publisher Full Text OpenURL

  19. Hurst LD, Pál C: Evidence for purifying selection acting on silent sites in BRCA1.

    Trends Genet 2001, 17:62-65. PubMed Abstract | Publisher Full Text OpenURL

  20. Orban TI, Olah E: Purifying selection on silent sites – a constraint from splicing regulation?

    Trends Genet 2001, 17:252-253. PubMed Abstract | Publisher Full Text OpenURL

  21. Xing Y, Lee C: Evidence of functional selection pressure for alternative splicing events that accelerate evolution of protein subsequences.

    Proc Natl Acad Sci U S A 2005, 102:13526-13531. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Davis JC, Petrov DA: Preferential duplication of conserved proteins in eukaryotic genomes.

    PLoS Biol 2004, 2:318-326. Publisher Full Text OpenURL

  23. Kondrashov FA, Rogozin IB, Wolf YI, Koonin EV: Selection in the evolution of gene duplications.

    Genome Biol 2002, 3:RESEARCH0008. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Conant GC, Wagner A: Asymmetric sequence divergence of duplicate genes.

    Genome Res 2003, 13:2052-2058. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Zhang P, Gu Z, Li WH: Different evolutionary patterns between young duplicate genes in the human genome.

    Genome Biol 2003, 4:R56. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  26. Duret L, Mouchiroud D: Determinants of substitution rates in mammalian genes: expression pattern affects selection intensity but not mutation rate.

    Mol Biol Evol 2000, 17:68-74. PubMed Abstract | Publisher Full Text OpenURL

  27. Pal C, Papp B, Hurst LD: Highly expressed genes in yeast evolve slowly.

    Genetics 2001, 158:927-931. PubMed Abstract | Publisher Full Text OpenURL

  28. Jordan IK, Wolf YI, Koonin EV: Duplicated genes evolve slower than singletons despite the initial rate increase.

    BMC Evol Biol 2004, 4:22. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Xing Y, Lee CJ: Protein Modularity of Alternatively Spliced Exons Is Associated with Tissue-Specific Regulation of Alternative Splicing.

    PLoS Genet 2005, 1:e34. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Parmley JL, Chamary JV, Hurst LD: Evidence for purifying selection against synonymous mutations in mammalian exonic splicing enhancers.

    Mol Biol Evol 2006, 23:301-309. PubMed Abstract | Publisher Full Text OpenURL

  31. Xing Y, Lee C: Can RNA selection pressure distort the measurement of Ka/Ks?

    Gene, in press.

    2006, Feb 16

    OpenURL

  32. Ina Y: New methods for estimating the numbers of synonymous and nonsynonymous substitutions.

    J Mol Evol 1995, 40:190-226. PubMed Abstract | Publisher Full Text OpenURL

  33. Yang Z: PAML: a program package for phylogenetics analysis by maximum likelyhood.

    Comput Appl Biosci 1997, 13:555-556. PubMed Abstract OpenURL

  34. Ermakova EO, Malko DB, Gelfand MS: Evolutionary differences between alternative and constant protein-coding regions of alternatively spliced Drosophila genes.

    Biofizika, in press. OpenURL

  35. Malko DB, Makeev VJ, Mironov AA, Gelfand MS: Evolution of the exon-intron structure and alternative splicing in fruit flies and malarial mosquito genomes.

    Genome Res 2006, 16(4):505-9.

    PubMed Abstract | Publisher Full Text OpenURL

  36. 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. PubMed Abstract | Publisher Full Text OpenURL

  37. Ramensky V, Neverov AD, Nurtdinov RN, Mironov AA, Gelfand MS: Human SNPs and alternative splicing.

    Proceedings of the Second International Moscow Conference on Computational Molecular Biology: 18–21 July 2005; Moscow 2005, 317-318. OpenURL

  38. NCBI Reference Sequence (RefSeq) [http://www.ncbi.nlm.nih.gov/RefSeq/] webcite

  39. Jordan IK, Kondrashov FA, Rogozin IB, Tatusov RL, Wolf YI, Koonin EV: Constant relative rate of protein evolution and detection of functional diversification among bacterial, archaeal and eukaryotic proteins.

    Genome Biol 2001, 2:research0053.1-0053.9. BioMed Central Full Text OpenURL

  40. EDAS: EST-Derived Alternative Splicing Database [http://www.belozersky.msu.ru/edas/] webcite

  41. Mironov AA, Novichkov PS, Gelfand MS: Pro-Frame: similarity-based gene recognition in eukaryotic DNA sequences with errors.

    Bioinformatics 2001, 17:13-15. PubMed Abstract | Publisher Full Text OpenURL