Evidence for a high frequency of simultaneous double-nucleotide substitutions

Science

Volumn 287, Issue 5456, 2000, Pages 1283-1286

  Averof, Michalis ^a   Rokas, Antonis ^b   Wolfe, Kenneth H ^c   Sharp, Paul M ^d  

^a INSTITUTE OF MOLECULAR BIOLOGY AND BIOTECHNOLOGY   (Greece)

^b UNIVERSITY OF EDINBURGH   (United Kingdom)

^c TRINITY COLLEGE DUBLIN   (Ireland)

^d UNIVERSITY OF NOTTINGHAM   (United Kingdom)

Author keywords

[No Author keywords available]

Indexed keywords

ARTICLE; CODON; GENE LOCUS; GENE MUTATION; HUMAN; NONHUMAN; NUCLEOTIDE SEQUENCE; PHYLOGENY; POINT MUTATION; PRIORITY JOURNAL;

ALLELES; AMINO ACID SUBSTITUTION; ANIMALS; CODON; EVOLUTION, MOLECULAR; GLOBINS; HAPLORHINI; HUMANS; MODELS, GENETIC; MUTATION; NUCLEOTIDES; PHYLOGENY; SERINE;

EID: 0034681511     PISSN: 00368075     EISSN: None     Source Type: Journal    
DOI: 10.1126/science.287.5456.1283     Document Type: Article

References (42)

1. G. B. Golding and B. W. Glickman, Proc. Natl. Acad. Sci. U.S.A. 82, 8577 (1985).
2. D. M. Hampsey, J. F. Ernst, J. W. Stewart, F. J. Sherman, J. Mol. Biol. 201, 471 (1988).
3. H. Nakazawa et al., Proc. Natl. Acad. Sci. U.S.A. 91, 360 (1994).
4. K. H. Wolfe and P. M. Sharp, J. Mol. Evol. 37, 441 (1993).
5. D. Mouchiroud, C. Gautier, G. Bernardi, J. Mol. Evol. 40, 107 (1995); J. M. Comeron and M. Kreitman, Genetics 150, 767 (1998); N. G. C. Smith and L. D. Hurst. Genetics 153, 1395 (1999).
6. D. Mouchiroud, C. Gautier, G. Bernardi, J. Mol. Evol. 40, 107 (1995); J. M. Comeron and M. Kreitman, Genetics 150, 767 (1998); N. G. C. Smith and L. D. Hurst. Genetics 153, 1395 (1999).
7. D. Mouchiroud, C. Gautier, G. Bernardi, J. Mol. Evol. 40, 107 (1995); J. M. Comeron and M. Kreitman, Genetics 150, 767 (1998); N. G. C. Smith and L. D. Hurst. Genetics 153, 1395 (1999).
8. R. Grantham, Science 185, 862 (1974).
9. P. M. Sharp and W.-H. Li, J. Mol. Evol. 25, 58 (1987).
10. S. Brenner, Nature 334, 528 (1988).
11. We collected data sets of homologous proteins showing wide phylogenetic conservation (representing diverse eukaryotic and eubacterial lineages) and unambiguous relationships (excluding multigene families, cases of horizontal transfer, concerted evolution). For each data set, conserved protein sequences were obtained from the SWISS-PROT database (release 8/98) using BLAST (22). Protein sequences were aligned using CLUSTALW (23) and were searched for unambiguously aligned sites where serine is absolutely conserved (i.e., present in all available sequences). The corresponding codons were determined from the respective nucleotide sequences, obtained from the GenBank/EMBL (European Molecular Biology Laboratory) database. Changes in serine codon type were determined in the most parsimonious way on the basis of phytogenies (illustrated Fig. 1B), and rates of change were estimated as the number of inferred changes over the time sampled at each site (the sum of all branch lengths). Phylogenetic relationships and times of divergence svere based on published data for the respective species (20, 24) (Fig. 1A). Trees were also constructed from the protein sequences themselves [using the Neighbor-Joining method (25)], and sequences showing an inconsistent phylogenetic placement were eliminated. Because of difficulties in determining their times of divergence, eubacterial and eukaryotic sequences were treated separately and archaebacterial sequences were excluded. In estimating the rates of codon switches, we tried to be conservative, for example, by overestimating times of divergence in cases of uncertainty. The analysis was also carried out on the basis of alternative published phylogenies and the results were always robust (M. Averof et al., data not shown).
12. Y. Diaz-Lazcoz, A. Henaut, P. Vigier, J. L. Risler, J. Mol. Biol. 250, 123 (1995).
13. M. Kimura, in Population Genetics and Molecular Evolution, T. Ohta and K. Aoki, Eds. (Springer-Verlag, Berlin, 1985), pp. 19-39.
14. N. Averof and P. M. Sharp, data not shown.
15. Ubiquitin, one of the most highly conserved proteins known, is encoded by repeats sharing recent common ancestry due to concerted evolution (7); serine switches among repeats have been reported [P. M. Sharp and W.-H. Li, Trends Ecol. Evol. 2, 328 (1987); P. M. Sharp, M. Averof, A. T. Lloyd, G. Matassi, J. F. Peden, Philos. Trans. R. Soc. London Ser. B 349, 241 (1995)]. We have also found serine switches in orthologous genes from the two strains of Helicobacter pylori for which complete genome sequences have been determined [J.-F. Tomb et al., Nature 388, 539 (1997); R. A. Alm et al., Nature 397, 176 (1999)]. Switches were seen at sites where serine is conserved in a range of other species.
16. Ubiquitin, one of the most highly conserved proteins known, is encoded by repeats sharing recent common ancestry due to concerted evolution (7); serine switches among repeats have been reported [P. M. Sharp and W.-H. Li, Trends Ecol. Evol. 2, 328 (1987); P. M. Sharp, M. Averof, A. T. Lloyd, G. Matassi, J. F. Peden, Philos. Trans. R. Soc. London Ser. B 349, 241 (1995)]. We have also found serine switches in orthologous genes from the two strains of Helicobacter pylori for which complete genome sequences have been determined [J.-F. Tomb et al., Nature 388, 539 (1997); R. A. Alm et al., Nature 397, 176 (1999)]. Switches were seen at sites where serine is conserved in a range of other species.
17. Ubiquitin, one of the most highly conserved proteins known, is encoded by repeats sharing recent common ancestry due to concerted evolution (7); serine switches among repeats have been reported [P. M. Sharp and W.-H. Li, Trends Ecol. Evol. 2, 328 (1987); P. M. Sharp, M. Averof, A. T. Lloyd, G. Matassi, J. F. Peden, Philos. Trans. R. Soc. London Ser. B 349, 241 (1995)]. We have also found serine switches in orthologous genes from the two strains of Helicobacter pylori for which complete genome sequences have been determined [J.-F. Tomb et al., Nature 388, 539 (1997); R. A. Alm et al., Nature 397, 176 (1999)]. Switches were seen at sites where serine is conserved in a range of other species.
18. Ubiquitin, one of the most highly conserved proteins known, is encoded by repeats sharing recent common ancestry due to concerted evolution (7); serine switches among repeats have been reported [P. M. Sharp and W.-H. Li, Trends Ecol. Evol. 2, 328 (1987); P. M. Sharp, M. Averof, A. T. Lloyd, G. Matassi, J. F. Peden, Philos. Trans. R. Soc. London Ser. B 349, 241 (1995)]. We have also found serine switches in orthologous genes from the two strains of Helicobacter pylori for which complete genome sequences have been determined [J.-F. Tomb et al., Nature 388, 539 (1997); R. A. Alm et al., Nature 397, 176 (1999)]. Switches were seen at sites where serine is conserved in a range of other species.
19. D. M. Irwin, Nature 336, 429 (1988).
20. P. S. Goldfarb, Nature 336, 429 (1988).
21. E. V. Koonin and A. E. Gorbalenya, Nature 338, 467 (1989).
22. M. M. Miyamoto, B. F. Koop, J. L. Slightom, M. Goodman, Proc. Natl. Acad. Sci. U.S.A. 85, 7627 (1988); W. J. Bailey et al., Mol. Phylogenet. Evol. 1, 97 (1992).
23. M. M. Miyamoto, B. F. Koop, J. L. Slightom, M. Goodman, Proc. Natl. Acad. Sci. U.S.A. 85, 7627 (1988); W. J. Bailey et al., Mol. Phylogenet. Evol. 1, 97 (1992).
24. 2L, where L is the total number of aligned nucleotides (sample size), ObsS and ObsD are the numbers of changes observed as single- and double-nucleotide substitutions, respectively, RealS and RealD are the numbers of changes that have occurred as single- and simultaneous double-nucleotide substitutions, respectively, and ExpD is the number of doublet substitutions (observed at adjacent nucleotides) that have occurred by coincidence of two separate single-nucleotide changes. The difference of ObsD and ExpD was evaluated by a chi-squared test. In mammalian genomes, CG dinucleotides are methylated, which renders them susceptible to mutation (by deamination), yielding TC and CA dinucleotides at relatively high frequencies [A. P. Bird, Nucleic Acids Res. 8, 1499 (1980)]. To prevent sequential CN→CG→TG and NG→CG→CA changes (or parallel CC→TG and CG→CA mutations in independent lineages) from being counted as doublet mutations, CN→TG and NG→CA substitutions were excluded from the analysis. We also examined the rates of substitution of each of the four nucleotides separately; the biases in these rates are not sufficient to affect our results.
25. The rate of double-nucleotide substitution was estimated as 2% of the overall rate of single-nucleotide substitution (RealD/ObsS in Table 2). Average rates of neutral nucleotide substitutions are around 5 per site per Gyr (20, 27), yielding an estimated doublet rate of 0.1 per site per Gyr. The excess of double-nucleotide substitutions (the difference between ExpD and ObsD) is highly significant, with a chisquared value of 14.2 (df = 1, P ≪ 0.005).
26. H. Ochman and A. C. Wilson, J. Mol. Evol. 26, 74 (1987).
27. W.-H. Li, Molecular Evolution (Sinauer, Sunderland, MA, 1997); B. S. Caut, Evolutionary Biology, M. K. Hecht et al., Eds. (Plenum, New York, 1998), pp. 93-120; W. Makalowski and M. S. Boguski, Proc. Natl. Acad. Sci. U.S.A. 95, 9407 (1998).
28. W.-H. Li, Molecular Evolution (Sinauer, Sunderland, MA, 1997); B. S. Caut, Evolutionary Biology, M. K. Hecht et al., Eds. (Plenum, New York, 1998), pp. 93-120; W. Makalowski and M. S. Boguski, Proc. Natl. Acad. Sci. U.S.A. 95, 9407 (1998).
29. W.-H. Li, Molecular Evolution (Sinauer, Sunderland, MA, 1997); B. S. Caut, Evolutionary Biology, M. K. Hecht et al., Eds. (Plenum, New York, 1998), pp. 93-120; W. Makalowski and M. S. Boguski, Proc. Natl. Acad. Sci. U.S.A. 95, 9407 (1998).
30. S. F. Altschul, W. Cish, W. Miller, E. W. Myers, D. J. Lipman, J. Mol. Biol. 215, 403 (1990).
31. J. D. Thompson, D. C. Higgins, T. J. Gibson, Nucleic Acids Res. 22, 4673 (1994).
32. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
33. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
34. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
35. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
36. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
37. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
38. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
39. G. J. Olsen, C. R. Woese, R. Overbeek, J. Bacteriol. 176, 1 (1994); D. F. Feng, G. Cho, R. F. Doolittle, Proc. Natl. Acad. Sci. U.S.A. 94, 13028 (1997); T. M. Embley, R. P. Hirt, Curr. Opin. Genet. Dev. 8, 624 (1998); J. D. Palmer, R. K. Jansen, H. J. Michaels, M. W. Chase, J. R. Manhart, Ann. Mo. Bot. Card. 75, 1180 (1988); K. H. Wolfe, M. Gouy, V. W. Yang, P. M. Sharp, W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 86, 6201 (1989); F. J. Ayala, A. Rzhetsky, F. J. Ayala, Proc. Natl. Acad. Sci. U.S.A. 95, 606 (1998); A. M. A. Aguinaldo et al., Nature 367, 489 (1997); S. Kumar and S. B. Hedges, Nature 392, 917 (1998).
40. N. Saitou and M. Nei, Mol. Biol. Evol 4, 406 (1987).
41. D. Swofford, PAUP (Phylogenetic Analysis Using Parsimony), version 3.0q (Illinois Natural History Survey, Champaign, IL, 1990).
42. This work was initiated in the Genetics Department of Trinity College, University of Dublin, and was continued with the support of C. Delidakis at the IMBB, Crete. It was supported in part by the EPET II programme of the General Secretariat for Research and Technology, Greece, and a UK Biotechnology and Biological Sciences Research Council grant G04905.