Limitations of a molecular clock applied to considerations of the origin of HIV-1

Science

Volumn 280, Issue 5371, 1998, Pages 1868-1871

  Korber, Bette ^a,b   Theiler, James ^c   Wolinsky, Steven ^d  

^a LOS ALAMOS NATIONAL LABORATORY   (United States)

^b SANTA FE INSTITUTE   (United States)

^c LOS ALAMOS NATIONAL LABORATORY   (United States)

^d NORTHWESTERN UNIVERSITY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

HUMAN IMMUNODEFICIENCY VIRUS VACCINE; VIRUS ENVELOPE PROTEIN;

ACQUIRED IMMUNE DEFICIENCY SYNDROME; EVOLUTION; HUMAN IMMUNODEFICIENCY VIRUS 1; HUMAN IMMUNODEFICIENCY VIRUS INFECTION; NONHUMAN; PHYLOGENY; PRIORITY JOURNAL; SHORT SURVEY; VIRAL GENETICS; VIRUS ENVELOPE; VIRUS RECOMBINANT;

HUMAN IMMUNODEFICIENCY VIRUS; HUMAN IMMUNODEFICIENCY VIRUS 1;

EID: 0032546809     PISSN: 00368075     EISSN: None     Source Type: Journal    
DOI: 10.1126/science.280.5371.1868     Document Type: Short Survey

References (78)

1. J. M. Coffin, Science 267, 483 (1995); B. D. Preston, B. J. Poiesz, L. A. Loeb, ibid. 242, 1168 (1988); K. Bebenek, J. Abbotts, S. H. Wilson, T. A. Kunkel, J. Biol. Chem. 268, 10324 (1993).
2. J. M. Coffin, Science 267, 483 (1995); B. D. Preston, B. J. Poiesz, L. A. Loeb, ibid. 242, 1168 (1988); K. Bebenek, J. Abbotts, S. H. Wilson, T. A. Kunkel, J. Biol. Chem. 268, 10324 (1993).
3. J. M. Coffin, Science 267, 483 (1995); B. D. Preston, B. J. Poiesz, L. A. Loeb, ibid. 242, 1168 (1988); K. Bebenek, J. Abbotts, S. H. Wilson, T. A. Kunkel, J. Biol. Chem. 268, 10324 (1993).
4. E. C. Sabino et al., J. Virol. 68, 6340 (1994); D. L. Robertson, B. H. Hahn, P. M. Sharp, Mol. Evol. 40, 249 (1995); F. E. McCutchan, M. O. Salminen, J. K. Carr, D. S. Burke, AIDS 10, S13 (1996); M. Cornelissen et al., J. Virol. 70, 8209 (1996).
5. E. C. Sabino et al., J. Virol. 68, 6340 (1994); D. L. Robertson, B. H. Hahn, P. M. Sharp, Mol. Evol. 40, 249 (1995); F. E. McCutchan, M. O. Salminen, J. K. Carr, D. S. Burke, AIDS 10, S13 (1996); M. Cornelissen et al., J. Virol. 70, 8209 (1996).
6. E. C. Sabino et al., J. Virol. 68, 6340 (1994); D. L. Robertson, B. H. Hahn, P. M. Sharp, Mol. Evol. 40, 249 (1995); F. E. McCutchan, M. O. Salminen, J. K. Carr, D. S. Burke, AIDS 10, S13 (1996); M. Cornelissen et al., J. Virol. 70, 8209 (1996).
7. E. C. Sabino et al., J. Virol. 68, 6340 (1994); D. L. Robertson, B. H. Hahn, P. M. Sharp, Mol. Evol. 40, 249 (1995); F. E. McCutchan, M. O. Salminen, J. K. Carr, D. S. Burke, AIDS 10, S13 (1996); M. Cornelissen et al., J. Virol. 70, 8209 (1996).
8. D. L. Robertson, F. Gao, B. H. Hahn, P. M. Sharp, in Human Retroviruses and AIDS, B. Korber et al., Eds. (Los Alamos National Laboratory, Los Alamos, NM, 1997), part III, pp. 25-30, http://hiv-web.lanl.gov/ HTML/reviews/hahn.html.
9. S. M. Wolinsky et al., Science 255, 1134 (1992); C. Wills, G. Myers, A. Farmer, AIDS Res. Hum. Retroviruses 12, 1383 (1996).
10. S. M. Wolinsky et al., Science 255, 1134 (1992); C. Wills, G. Myers, A. Farmer, AIDS Res. Hum. Retroviruses 12, 1383 (1996).
11. P. M. Sharp, D. L. Robertson, F. Gao, B. H. Hahn, AIDS 8, S27 (1994).
12. F. Clavel et al., Nature 324, 691 (1986); M. Grez et al., J. Virol. 68, 2161 (1994); B. Korber, E. Allen, A. Farmer, G. Myers, AIDS 9 (suppl. A), S5 (1995); P. J. Kanki, M. Peeters, A. Gueye-Ndiaye, ibid. 11 (suppl. B), S33 (1997).
13. F. Clavel et al., Nature 324, 691 (1986); M. Grez et al., J. Virol. 68, 2161 (1994); B. Korber, E. Allen, A. Farmer, G. Myers, AIDS 9 (suppl. A), S5 (1995); P. J. Kanki, M. Peeters, A. Gueye-Ndiaye, ibid. 11 (suppl. B), S33 (1997).
14. F. Clavel et al., Nature 324, 691 (1986); M. Grez et al., J. Virol. 68, 2161 (1994); B. Korber, E. Allen, A. Farmer, G. Myers, AIDS 9 (suppl. A), S5 (1995); P. J. Kanki, M. Peeters, A. Gueye-Ndiaye, ibid. 11 (suppl. B), S33 (1997).
15. F. Clavel et al., Nature 324, 691 (1986); M. Grez et al., J. Virol. 68, 2161 (1994); B. Korber, E. Allen, A. Farmer, G. Myers, AIDS 9 (suppl. A), S5 (1995); P. J. Kanki, M. Peeters, A. Gueye-Ndiaye, ibid. 11 (suppl. B), S33 (1997).
16. P. Chameau et al., Virology 205, 247 (1994); I. Loussert-Ajaka et al., J. Virol. 69, 5640 (1995).
17. P. Chameau et al., Virology 205, 247 (1994); I. Loussert-Ajaka et al., J. Virol. 69, 5640 (1995).
18. T. Huet et al., Nature 345, 356 (1990).
19. D. L. Swofford, G. J. Olsen, P. J. Waddell, D. M. Hillis, in Molecular Systematics, D. M. Hillis, C. Moritz, B. K. Mable, Eds. (Sinauer, Sunderland, MA, ed. 2, 1996), pp. 407-514.
20. J. Louwagie et al., AIDS 7, 769 (1993); Workshop Report from the European Commission (DG XII, INCO-DC) and the Joint United Nations Programme on HIV/AIDS, AIDS 11, 17 (1997).
21. J. Louwagie et al., AIDS 7, 769 (1993); Workshop Report from the European Commission (DG XII, INCO-DC) and the Joint United Nations Programme on HIV/AIDS, AIDS 11, 17 (1997).
22. T. Leitner, B. Korber, D. Robertson, F. Gao, B. Hahn, in (3), pp. 19-26, http://hiv-web.lanl.gov/HTML/ reviews/Thomas.html.
23. F. E. McCutchan et al., in The Proceedings of the Fifth Annual HIV Dynamics and Evolution Meeting, Santa Fe, NM, 18 to 20 April 1998, B. Korber, Ed. (The Santa Fe Institute, Santa Fe, NM, 1998), vol. 5, p. 23, www.santafe.edu/hiv5/mccutchan.html; M. Peeters et al., ibid., p. 73, www.santafe.edu/hiv5/ peeters.html.
24. F. E. McCutchan et al., in The Proceedings of the Fifth Annual HIV Dynamics and Evolution Meeting, Santa Fe, NM, 18 to 20 April 1998, B. Korber, Ed. (The Santa Fe Institute, Santa Fe, NM, 1998), vol. 5, p. 23, www.santafe.edu/hiv5/mccutchan.html; M. Peeters et al., ibid., p. 73, www.santafe.edu/hiv5/ peeters.html.
25. W.-H. Li, M. Tanimura, P. M. Sharp, Mol. Biol. Evol. 5, 313 (1988); P. Sharp and W.-H. Li, Nature 333, 315 (1988); T. Gojobori, E. N. Moriyama, M. Kimura. Proc. Natl. Acad. Sci. U.S.A. 87, 10015 (1990).
26. W.-H. Li, M. Tanimura, P. M. Sharp, Mol. Biol. Evol. 5, 313 (1988); P. Sharp and W.-H. Li, Nature 333, 315 (1988); T. Gojobori, E. N. Moriyama, M. Kimura. Proc. Natl. Acad. Sci. U.S.A. 87, 10015 (1990).
27. W.-H. Li, M. Tanimura, P. M. Sharp, Mol. Biol. Evol. 5, 313 (1988); P. Sharp and W.-H. Li, Nature 333, 315 (1988); T. Gojobori, E. N. Moriyama, M. Kimura. Proc. Natl. Acad. Sci. U.S.A. 87, 10015 (1990).
28. T. Leitner et al., Proc. Natl. Acad. Sci. U.S.A. 93, 10864 (1996).
29. "Absolute, true, and mathematical time, of itself, and from its own nature, flows equably without relation to anything external." Isaac Newton, Principia (1687), F. Cajori, transl. (Univ. of California Press, Berkeley, CA, 1947), p. 6; "Wrongagain, Isaac," paraphrased by P. Ginsparg (personal communication), from Albert Einstein, Special Theory of Relativity (1905).
30. "Absolute, true, and mathematical time, of itself, and from its own nature, flows equably without relation to anything external." Isaac Newton, Principia (1687), F. Cajori, transl. (Univ. of California Press, Berkeley, CA, 1947), p. 6; "Wrongagain, Isaac," paraphrased by P. Ginsparg (personal communication), from Albert Einstein, Special Theory of Relativity (1905).
31. "One hand only. Finest of fine black darts. It advances by fits and starts. No tick. Leaps from dot to dot with so lightning a leap that but for its new position it had not stirred. Whole nights may pass as but a fraction of a second or any intermediate lapse of time soever before it flings itself from one degree to the next." Samuel Beckett, author's translation, III Seen, III Said (Grove, New York, 1981), p. 76.
32. D. M. Hillis, B. K. Mable, C. Moritz, in (9), pp. 515-543.
33. The use of the distance to a common M group ancestral node to estimate the year of sampling of a sequence gives a root mean square error of 8.7 years. Given that the virus that causes AIDS was identified only 15 years ago, this is not much predictive power. This calculation was specifically based on a linear regression analysis of genetic distance : between contemporary sequences and the estimated ancestral origin of the HIV M group, shown in Fig. 1, plotted as a function of the year of sampling. The linear correlation coefficient (Pearson's r) is 0.282 for the full set of data in Fig. 1, with a Student's t probability of 0.002, that is, a highly significant but weak correlation. Using a regression analysis to estimate the rate of divergence and a bootstrap approach to . estimate the error (45), we calculated annual divergence rates of 0.00169 per year (95% confidence interval of 0.0011 to 0.0022). These rates are lower than numbers given in the past (approaching 1% per year) [C. Kuiken and B. Korber, AIDS 8, S73 (1994)], perhaps because these calculations were based on sequences of the full-length envelope protein gp160, including the conserved regions, and typically shorter variable domains of envelope have been studied.
34. M. O. Salminen et al., J. Virol. 71, 2647 (1997); R. S. Diaz et al., ibid. 69, 3273 (1995).
35. M. O. Salminen et al., J. Virol. 71, 2647 (1997); R. S. Diaz et al., ibid. 69, 3273 (1995).
36. T. Zhu et al., Nature 391, 594 (1998).
37. The branch location of the ZR.59 sample, between the B, D, and F clades, could have been the result of ZR.59 being a mosaic genome; however, the branch length was extremely short between the ancestral node of these clades and the ZR.59 sequence, indicating that it was simply close to the ancestral sequence. There was no evidence that the ZR.59 sequence was a mosaic when it was analyzed carefully for that possibility with a sliding window approach. However, evidence for a preclade recombination event emerged from the study. The B and D subtypes are generally the most closely associated subtypes throughout the genome. But in the fragments of the genome where sample ZR.59 was sequenced, B and F subtype sequences were most closely associated. [Additional evidence for early recombination events is that the subtypes cluster with different patterns of association in phylogenetic trees on the basis of different regions of the genome (11)].
38. S. Wolinsky et al., Science 272, 537 (1996); S. Ganeshan, R. Dickover, B. Korber, Y. Bryson, S. Wolinsky, J. Virol. 71, 663 (1997); S.-L. Liu et al., ibid., p. 4284; E. Delwart et al., ibid., p. 7498; V. V. Lukashov et al., ibid. 69, 6911 (1995).
39. S. Wolinsky et al., Science 272, 537 (1996); S. Ganeshan, R. Dickover, B. Korber, Y. Bryson, S. Wolinsky, J. Virol. 71, 663 (1997); S.-L. Liu et al., ibid., p. 4284; E. Delwart et al., ibid., p. 7498; V. V. Lukashov et al., ibid. 69, 6911 (1995).
40. S. Wolinsky et al., Science 272, 537 (1996); S. Ganeshan, R. Dickover, B. Korber, Y. Bryson, S. Wolinsky, J. Virol. 71, 663 (1997); S.-L. Liu et al., ibid., p. 4284; E. Delwart et al., ibid., p. 7498; V. V. Lukashov et al., ibid. 69, 6911 (1995).
41. S. Wolinsky et al., Science 272, 537 (1996); S. Ganeshan, R. Dickover, B. Korber, Y. Bryson, S. Wolinsky, J. Virol. 71, 663 (1997); S.-L. Liu et al., ibid., p. 4284; E. Delwart et al., ibid., p. 7498; V. V. Lukashov et al., ibid. 69, 6911 (1995).
42. S. Wolinsky et al., Science 272, 537 (1996); S. Ganeshan, R. Dickover, B. Korber, Y. Bryson, S. Wolinsky, J. Virol. 71, 663 (1997); S.-L. Liu et al., ibid., p. 4284; E. Delwart et al., ibid., p. 7498; V. V. Lukashov et al., ibid. 69, 6911 (1995).
43. M. A. Nowak et al., Nature 375, 606 (1995).
44. A. S. Perelson, A. U. Neumann, M. Markowitz, J. M. Leonard, D. D. Ho, Science 271, 1582 (1996).
45. M. Eigen, Trends Microbiol. 4, 216 (1996).
46. P. Borrow et al., Nature Med. 3, 205 (1997); D. A. Price et al., Proc. Natl. Acad. Sci. U.S.A. 94, 1890 (1997).
47. P. Borrow et al., Nature Med. 3, 205 (1997); D. A. Price et al., Proc. Natl. Acad. Sci. U.S.A. 94, 1890 (1997).
48. H. F. Gunthard et al., J. Virol. 72, 2422 (1998); J. Hammond, B. A. Larder, R. F. Schinazi, J. W. Mellors, in (3), pp. 206-249.
49. H. F. Gunthard et al., J. Virol. 72, 2422 (1998); J. Hammond, B. A. Larder, R. F. Schinazi, J. W. Mellors, in (3), pp. 206-249.
50. R. Doms and J. P. Moore, in (3), p. 1, http://hiv-web.lanl.gov/HTML/reviews/doms.html.
51. P. Dighe, B. Korber, B. Foley, in (3), p. 74, http://hiv-web.lanl.gov/HTML/97compendium.html.
52. B. T. M. Korber, K. MacInnes, R. Smith, G. Myers, J. Virol. 68, 6730 (1994); C. L. Kuiken, B. T. Foley, E. Guzman, B. T. M. Korber, in Molecular Evolution of HIV, K. Crandall, Ed. (Johns Hopkins Univ. Press, Baltimore, MD, in press).
53. B. T. M. Korber, K. MacInnes, R. Smith, G. Myers, J. Virol. 68, 6730 (1994); C. L. Kuiken, B. T. Foley, E. Guzman, B. T. M. Korber, in Molecular Evolution of HIV, K. Crandall, Ed. (Johns Hopkins Univ. Press, Baltimore, MD, in press).
54. H. Cao et al., J. Virol. 71, 8615 (1997); G. Ferrari, et al., Proc. Natl. Acad. Sci. U.S.A. 94, 1396 (1997).
55. H. Cao et al., J. Virol. 71, 8615 (1997); G. Ferrari, et al., Proc. Natl. Acad. Sci. U.S.A. 94, 1396 (1997).
56. Using a method that determines the best fit line for given data points including error estimates in both dimensions, one can determine how different assumptions of error influence the slope of a line [W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. Flannery, Numerical Recipes in C: The Art of Scientific Computing (Cambridge Univ. Press, Cambridge, ed. 2, 1992), section 15.3, p. 666]. Use of the Poisson error as the fixed error for the distance and variation of the time error yields a range of quite different slopes (and as a consequence, the estimated year of origin of the M group) for the data shown in Fig. 1. This ranges from the extreme estimates of 1917 (no error in the time axis) to 1972 (at the limit of large error in the time axis, the equivalent to a fit with no error in the distance axis). If one assumes that the year of sampling could be confounded by a small pool of latently infected cells, carrying virus that originated some years before sampling, the predicted year of origin is intermediate. For example, if one assumes a 6-year error in the time axis, the estimated year of an ancestral sequence of the M group is 1942 (44).
57. J. P. Moore and D. D. Ho, AIDS 9 (suppl. A), S117 (1995); M. P. D'Souza, D. Livnat, J. A. Bradac, S. H. Bridges, Infect. Dis. 175, 1056 (1997); S. Beddows et al., J. Gen. Virol. 79, 77 (1998).
58. J. P. Moore and D. D. Ho, AIDS 9 (suppl. A), S117 (1995); M. P. D'Souza, D. Livnat, J. A. Bradac, S. H. Bridges, Infect. Dis. 175, 1056 (1997); S. Beddows et al., J. Gen. Virol. 79, 77 (1998).
59. J. P. Moore and D. D. Ho, AIDS 9 (suppl. A), S117 (1995); M. P. D'Souza, D. Livnat, J. A. Bradac, S. H. Bridges, Infect. Dis. 175, 1056 (1997); S. Beddows et al., J. Gen. Virol. 79, 77 (1998).
60. P. Berman et al., J. Virol. 66, 4464 (1992); P. Berman et al., J. Infect Dis. 176, 384 (1997); "FDA authorizes first full testing for HIV vaccine," New York Times, 4 June 1998, p. A1.
61. P. Berman et al., J. Virol. 66, 4464 (1992); P. Berman et al., J. Infect Dis. 176, 384 (1997); "FDA authorizes first full testing for HIV vaccine," New York Times, 4 June 1998, p. A1.
62. P. Berman et al., J. Virol. 66, 4464 (1992); P. Berman et al., J. Infect Dis. 176, 384 (1997); "FDA authorizes first full testing for HIV vaccine," New York Times, 4 June 1998, p. A1.
63. B. Korber et al., AIDS Res. Hum. Retroviruses 8, 1461 (1992); S. D. Rencher, T. D. Lockey, K. S. Slobod, J. L. Hurwitz, ibid. 13, 527 (1997).
64. B. Korber et al., AIDS Res. Hum. Retroviruses 8, 1461 (1992); S. D. Rencher, T. D. Lockey, K. S. Slobod, J. L. Hurwitz, ibid. 13, 527 (1997).
65. R. Connor et al., ibid. 72, 1552 (1998).
66. Los Alamos, NM (1997). Theoretical Biology and Biophysics, Los Alamos National Laboratory, http:// hiv-web.lanl.gov/HTML/ENV-index.html.
67. J. K. Carr et al., J. Virol. 70, 5935 (1996); F. Gao et al., ibid., p. 7013.
68. J. K. Carr et al., J. Virol. 70, 5935 (1996); F. Gao et al., ibid., p. 7013.
69. F. Gao et al., ibid. 70, 1651 (1996); F. McCutchan et al., AIDS Res. Hum. Retroviruses 14, 329 (1998).
70. F. Gao et al., ibid. 70, 1651 (1996); F. McCutchan et al., AIDS Res. Hum. Retroviruses 14, 329 (1998).
71. Two-letter country codes indicating the origins of the viruses included in the envelope tree in Fig. 1 are as follows: Malawi (MW), Burundi (BU), Ethiopia (ET), Brazil (BR), Central African Republic (CF), Thailand (TH), Uganda (UG), Kenya (KE), Rwanda (RW), Zaire (ZR), United States (US), Australia (AU), France (FR), and Haiti (HT).
72. G. J. Olsen et al., Comput. Appl. Biosci. 10, 41 (1994).
73. FastDNAml and DNArates were written by Gary Olsen and colleagues at the Ribosomal Database Project (RDP) at the University of Illinois at Urbana-Champaign, available by anonymous ftp from ftp:// rdp.life.uiuc.edu/.
74. To create this tree, we generated a preliminary maximum likelihood tree which was then used as input for the DNArates program (47) that produced the maximum likelihood estimate of the rate of nucleotide substitution at every position in the alignment, allowing nine categories of rates. The preliminary tree and rates and categories were used to create a new tree, allowing global rearrangements and adjusting of branch lengths. Maximum likelihood was used to generate the tree to get the most accurate branch lengths possible (9, 11). It took 3 weeks of computer time on a SPARC workstation to generate this tree, so bootstrap analysis was done by using the neighbor joining package of Phylip with the maximum likelihood model to calculate the genetic distances (43).
75. J. Felsenstein, Evolution 39, 783 (1985); Cladistics 5, 164 (1989).
76. J. Felsenstein, Evolution 39, 783 (1985); Cladistics 5, 164 (1989).
77. B. Korber and J. Theiler, unpublished observations.
78. We thank C. Kuiken, T. Leitner, K. Rock, and B. Foley for helpful discussions, R. Thakallapally and A. Ploss for their help with the preparation of this manuscript, and P. Ginsparg for suggestions concerning the text. B.K. was supported by interagency agreements between the U.S. Department of Energy (DOE) and NIH, 3-Y01-AI-70001-13 and Y1-A1-4058-03; J.T. by the DOE; and S.M.W. by grants RO1 -HD-31756 and PO1 -AI-42454 from NIH, and a gift from an anonymous foundation.