Accumulation of GC donor splice signals in mammals
© Churbanov et al; licensee BioMed Central Ltd. 2008
Received: 07 July 2008
Accepted: 09 July 2008
Published: 09 July 2008
The GT dinucleotide in the first two intron positions is the most conserved element of the U2 donor splice signals. However, in a small fraction of donor sites, GT is replaced by GC. A substantial enrichment of GC in donor sites of alternatively spliced genes has been observed previously in human, nematode and Arabidopsis, suggesting that GC signals are important for regulation of alternative splicing. We used parsimony analysis to reconstruct evolution of donor splice sites and inferred 298 GT > GC conversion events compared to 40 GC > GT conversion events in primate and rodent genomes. Thus, there was substantive accumulation of GC donor splice sites during the evolution of mammals. Accumulation of GC sites might have been driven by selection for alternative splicing.
This article was reviewed by Jerzy Jurka and Anton Nekrutenko. For the full reviews, please go to the Reviewers' Reports section.
In vertebrates, most of the protein-coding genes are interrupted by multiple introns that are removed at the donor and acceptor splice sites so that the adjacent exons are spliced. This process is mediated by an elaborate molecular machine, the spliceosome that consists of 5 snRNPs (small nuclear ribonucleoprotein particles) along with numerous less stably associated proteins, and is conserved throughout the eukaryotic world [1–3]. The U2 spliceosome (the major eukaryotic spliceosome) interacts with specific parts of the intron and the flanking exons to ensure accurate and efficient splicing . The nucleotides at the intron termini and the adjacent nucleotides in the exons are involved in these interactions and comprise the splicing signal. The (A/C)AG|GT(A/G)AGT consensus sequence (the exon|intron boundary is shown by the vertical streak and the first two nucleotides of the intron are underlined) at the donor splice signal is complementary to the 5' end of U1 snRNA, and this interaction is believed to be the major requirement for splicing [5–7].
The GT dinucleotide in the first two intron positions is the most conserved element of the U2 donor splice signal. However, in a small fraction of donor sites (<1%), GT is replaced by GC; in these cases, the rest of the nucleotides in the donor signal adhere more closely to the consensus sequence, apparently, compensating for the T to C substitution that is unfavorable for splicing [8–10]. This rare class of donor splice signals has been implicated in alternative splicing [9, 11, 12]. For example, the conserved C at the +2 position of the 10th intron of the let-2 gene which encodes one of the collagen isoforms is essential for developmentally regulated alternative splicing in the nematode C. elegans. Replacement of the GC donor signal with a moderate or strong GT signal abolishes splicing regulation and leads to excessive usage of exon 10 of let-2 in embryos . Generally, a substantial enrichment of GC donor signals in alternatively spliced genes has been observed in human, C. elegans and Arabidopsis [9, 11, 12].
Pairwise comparisons of GC splicing signals in the nematodes Caenorhabditis elegans and C. briggsae suggested that GC donor signals are not evolutionarily conserved in nematodes: among the 26 C. elegans GC-AG introns, only 5 had a GC-AG counterpart in C. briggsae . Frequent switching between GT-AG and GC-AG introns has been reported for 5 vertebrate genomes . We were interested in exploring the genome-wide evolutionary dynamics of the donor splice sites and, in particular, sought to determine whether there might be a trend toward depletion or accumulation of GC.
The GC>GT and GT>GC conversion events were reconstructed using maximum parsimony (Figure 1). Unexpectedly, we observed a pronounced excess of GT>GC conversion over GC>GT conversion that is indicative of accumulation of GC donor splice sites in both primate and rodent genomes (Figure 1). The trend is stronger in the rodent lineage than in the primate lineage (Figure 1), an observation that is consistent with the overall fast genome evolution in rodents . The apparent accumulation of GC donor signals was further supported by the analysis of the terminal branches of the tree although the excess of GC>GT conversions in macaque compared to human (Figure 1) could be caused by sequencing errors and/or population polymorphism. The observed excess of GT>GC conversions was robust with respect to the composition of the outgroup species set (Additional file 2).
The observed excess of GT>GC conversions hardly can be explained by a nucleotide substitution bias. It has been repeatedly shown that mammalian genomes have a tendency to become more AT-rich [17–19]. Even if one assumes that, due to unknown reasons, this trend is reversed in the donor sites so that T to C substitutions are twice as frequent as C to T substitutions, such a bias would not account for the observed excess of GT>GC conversion events (P < 10-10 according to the χ2 test for pooled conversion events).
Considering that mutational bias did not seem to be a plausible cause of the observed accumulation of GC donor sites, it seems most likely that this trend has to do with the involvement of GC sites in alternative splicing that is widespread and essential in mammals [9, 11, 12]. As GT>GC conversion can substantially alter the pattern of alternative splicing , these changes might become beneficial and eventually would be fixed in the population. Thus, positive selection could be a plausible explanation for the observed accumulation of GC in donor sites. However, an even more plausible scenario would involve evolution of a strong splice site context that would allow neutral fixation of GC sites. The neutrally fixed GC sites, then, could be recruited for alternative splicing and thus would become subject to purifying selection forbidding the reverse GC>GT conversion.
Reviewer's report 1: Jerzy Jurka, Genetic Information Research Institute
This paper reports the relatively straightforward observation that GC donor splice sites tend to accumulate during evolution of mammals. However, the suggestion of selection for alternative splicing would be more convincing if the authors could demonstrate the GC accumulation separately in AT-rich and GC-rich genomic regions in mammals, where the dynamics of GT replacement by GC may be different.
We appreciate the suggestion that the dynamics of GC accumulation could depend on the base composition in the respective regions of the mammalian genomes.
GT and GC donor splice sites and conversion events in AT-rich and GC-rich regionsa
Splice sites/conversion events
%A+T > 50%
%A+T < 50%
Total number of donor splice sites(%)
Reviewer's report 2: Anton Nekrutenko, Pennsylvania State University
In this discovery note authors point out accumulation of non- canonical GC donor splice signals in mammals, against the previously observed nucleotide substitution bias. They provide a convincing explanation suggesting that GT->GC conversion may be beneficial for mammals as it creates additional possibilities for alternative splicing events. In my opinion these observations provide a platform for launching more detailed investigation of alternative splicing through comparative genomics and raise numerous interesting question (e.g., are any of the GC sites overlap with known SNPs?). Thus publication of this note will appeal to a broad evolutionary genomics community.
On a technical side the authors used a rather complex procedure for retrieving splice sites from TBA alignments. Instead, this can be easily and quickly done using Galaxy system:
as explained here:
We appreciate the reviewer pointing out the utility of the Galaxy platform and hope to exploit Galaxy in future genome analyses.
The research of IBR and EVK is supported by the Department of Health and Human Services intramural program (NIH, National Library of Medicine). The research of AC and SWH is supported by an NIH National Library of Medicine K-22 award (K22LM008794), and via private funding from the Research Institute for Children, New Orleans.
- Jurica MS, Moore MJ: Pre-mRNA splicing: awash in a sea of proteins. Mol Cell. 2003, 12 (1): 5-14. 10.1016/S1097-2765(03)00270-3.PubMedView ArticleGoogle Scholar
- Nilsen TW: The spliceosome: the most complex macromolecular machine in the cell?. Bioessays. 2003, 25 (12): 1147-1149. 10.1002/bies.10394.PubMedView ArticleGoogle Scholar
- Collins L, Penny D: Complex spliceosomal organization ancestral to extant eukaryotes. Mol Biol Evol. 2005, 22 (4): 1053-1066. 10.1093/molbev/msi091.PubMedView ArticleGoogle Scholar
- Patel AA, Steitz JA: Splicing double: insights from the second spliceosome. Nat Rev Mol Cell Biol. 2003, 4 (12): 960-970. 10.1038/nrm1259.PubMedView ArticleGoogle Scholar
- Rosbash M, Seraphin B: Who's on first? The U1 snRNP-5' splice site interaction and splicing. Trends Biochem Sci. 1991, 16 (5): 187-190. 10.1016/0968-0004(91)90073-5.PubMedView ArticleGoogle Scholar
- Du H, Rosbash M: The U1 snRNP protein U1C recognizes the 5' splice site in the absence of base pairing. Nature. 2002, 419 (6902): 86-90. 10.1038/nature00947.PubMedView ArticleGoogle Scholar
- Carmel I, Tal S, Vig I, Ast G: Comparative analysis detects dependencies among the 5' splice-site positions. Rna. 2004, 10 (5): 828-840. 10.1261/rna.5196404.PubMedPubMed CentralView ArticleGoogle Scholar
- Burset M, Seledtsov IA, Solovyev VV: Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Res. 2000, 28 (21): 4364-4375. 10.1093/nar/28.21.4364.PubMedPubMed CentralView ArticleGoogle Scholar
- Thanaraj TA, Clark F: Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions. Nucleic Acids Res. 2001, 29 (12): 2581-2593. 10.1093/nar/29.12.2581.PubMedPubMed CentralView ArticleGoogle Scholar
- Abril JF, Castelo R, Guigo R: Comparison of splice sites in mammals and chicken. Genome Res. 2005, 15 (1): 111-119. 10.1101/gr.3108805.PubMedPubMed CentralView ArticleGoogle Scholar
- Farrer T, Roller AB, Kent WJ, Zahler AM: Analysis of the role of Caenorhabditis elegans GC-AG introns in regulated splicing. Nucleic Acids Res. 2002, 30 (15): 3360-3367. 10.1093/nar/gkf465.PubMedPubMed CentralView ArticleGoogle Scholar
- Campbell MA, Haas BJ, Hamilton JP, Mount SM, Buell CR: Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis. BMC Genomics. 2006, 7: 327-10.1186/1471-2164-7-327.PubMedPubMed CentralView ArticleGoogle Scholar
- Blanchette M, Kent WJ, Riemer C, Elnitski L, Smit AF, Roskin KM, Baertsch R, Rosenbloom K, Clawson H, Green ED, Haussler D, Miller W: Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res. 2004, 14 (4): 708-715. 10.1101/gr.1933104.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuhn RM, Karolchik D, Zweig AS, Trumbower H, Thomas DJ, Thakkapallayil A, Sugnet CW, Stanke M, Smith KE, Siepel A, Rosenbloom KR, Rhead B, Raney BJ, Pohl A, Pedersen JS, Hsu F, Hinrichs AS, Harte RA, Diekhans M, Clawson H, Bejerano G, Barber GP, Baertsch R, Haussler D, Kent WJ: The UCSC genome browser database: update 2007. Nucleic Acids Res. 2007, 35 (Database issue): D668-73. 10.1093/nar/gkl928.PubMedPubMed CentralView ArticleGoogle Scholar
- Rogozin IB, Wolf YI, Carmel L, Koonin EV: Ecdysozoan clade rejected by genome-wide analysis of rare amino acid replacements. Mol Biol Evol. 2007, 24 (4): 1080-1090. 10.1093/molbev/msm029.PubMedView ArticleGoogle Scholar
- Thomas JW, Touchman JW, Blakesley RW, Bouffard GG, Beckstrom-Sternberg SM, Margulies EH, Blanchette M, Siepel AC, Thomas PJ, McDowell JC, Maskeri B, Hansen NF, Schwartz MS, Weber RJ, Kent WJ, Karolchik D, Bruen TC, Bevan R, Cutler DJ, Schwartz S, Elnitski L, Idol JR, Prasad AB, Lee-Lin SQ, Maduro VV, Summers TJ, Portnoy ME, Dietrich NL, Akhter N, Ayele K, Benjamin B, Cariaga K, Brinkley CP, Brooks SY, Granite S, Guan X, Gupta J, Haghighi P, Ho SL, Huang MC, Karlins E, Laric PL, Legaspi R, Lim MJ, Maduro QL, Masiello CA, Mastrian SD, McCloskey JC, Pearson R, Stantripop S, Tiongson EE, Tran JT, Tsurgeon C, Vogt JL, Walker MA, Wetherby KD, Wiggins LS, Young AC, Zhang LH, Osoegawa K, Zhu B, Zhao B, Shu CL, De Jong PJ, Lawrence CE, Smit AF, Chakravarti A, Haussler D, Green P, Miller W, Green ED: Comparative analyses of multi-species sequences from targeted genomic regions. Nature. 2003, 424 (6950): 788-793. 10.1038/nature01858.PubMedView ArticleGoogle Scholar
- Pozdniakov MA, Rogozin IB, Babenko VN, Kolchanov NA: Neighboring base effect on emergence of spontaneous mutations in human pseudogenes. Dokl Akad Nauk. 1997, 356 (4): 566-568.PubMedGoogle Scholar
- Petrov DA, Hartl DL: Patterns of nucleotide substitution in Drosophila and mammalian genomes. Proc Natl Acad Sci U S A. 1999, 96 (4): 1475-1479. 10.1073/pnas.96.4.1475.PubMedPubMed CentralView ArticleGoogle Scholar
- Jiang C, Zhao Z: Directionality of point mutation and 5-methylcytosine deamination rates in the chimpanzee genome. BMC Genomics. 2006, 7: 316-10.1186/1471-2164-7-316.PubMedPubMed CentralView ArticleGoogle 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.