Timing the Ancestor of the HIV-1 Pandemic Strains

American Association for the Advancement of Science (AAAS) - Tập 288 Số 5472 - Trang 1789-1796 - 2000
Bette Korber1,2, Mark Muldoon3,2, James Theiler1, Feng Gao4, R. Gupta1, Alan S. Lapedes1,2, Beatrice H. Hahn4, Steven M. Wolinsky5, Tanmoy Bhattacharya1
1Los Alamos National Laboratory, Los Alamos, NM 87545, USA
2Santa Fe Institute, Santa Fe, NM 87501, USA
3Department of Mathematics, University of Manchester Institute of Technology, Manchester M60 1QD, UK.
4Department of Medicine and Microbiology, University of Alabama at Birmingham, Birmingham, AL 35294, USA.
5Division of Infectious Diseases, Department of Medicine, Northwestern University, Chicago, IL 60611, USA.

Tóm tắt

HIV-1 sequences were analyzed to estimate the timing of the ancestral sequence of the main group of HIV-1, the strains responsible for the AIDS pandemic. Using parallel supercomputers and assuming a constant rate of evolution, we applied maximum-likelihood phylogenetic methods to unprecedented amounts of data for this calculation. We validated our approach by correctly estimating the timing of two historically documented points. Using a comprehensive full-length envelope sequence alignment, we estimated the date of the last common ancestor of the main group of HIV-1 to be 1931 (1915–41). Analysis of a gag gene alignment, subregions of envelope including additional sequences, and a method that relaxed the assumption of a strict molecular clock also supported these results.

Từ khóa


Tài liệu tham khảo

10.1126/science.287.5453.607

Peeters M., et al., AIDS 10, 625 (1989);

Peeters M., et al., AIDS 6, 447 (1992);

Gao F., et al., Nature 397, 436 (1999);

10.1128/JVI.74.1.529-534.2000

. S. Souquiere et al. paper presented at the 7th Conference on Retroviruses and Opportunistic Infections San Francisco 2000 (www.retroconference.org/).

Hirsch V., et al., Nature 339, 389 (1989);

Gao F., et al., Nature 358, 495 (1992);

; F. Gao et al. J. Virol. 68 7433 (1994);

Chen Z., et al., J. Virol. 70, 3617 (1996);

Chen A., et al., J. Virol. 71, 3953 (1997).

De Leys R., et al., J. Virol. 64, 1207 (1990);

Charneau P., et al., Virology 205, 247 (1994);

Gurtler L., Lancet 348, 176 (1996).

Simon F., et al., Nature Med. 4, 1032 (1998).

D. L. Robertson et al. in Human Retroviruses and AIDS 1999 C. Kuiken et al. Eds. (Los Alamos National Laboratory Los Alamos NM in press) (available at hiv-web.lanl.gov);

10.1126/science.288.5463.55d

Meyer A., et al., Med. Trop. 51, 53 (1991);

Heymann D., Szczeniowski M., Esteves K., Br. Med. Bull. 54, 693 (1998);

Cohen J., Science 277, 312 (1997).

Voevodin A., et al., Virology 238, 212 (1997);

Slattery J. P., Franchini G., Gessain A., Genome Res. 9, 525 (1999).

Heneine W., et al., Nature Med. 4, 403 (1998);

10.1128/JVI.73.11.9619-9624.1999

Report from the Joint United Nations Programme on HIV/AIDS Global HIV/AIDS epidemic update 1999 available at www.unaids.org/publications/.

Zhu T., et al., Nature 391, 594 (1998).

Jonassen T., et al., Virology 231, 43 (1997).

10.1126/science.6189183

10.1126/science.6601823

Leitner T., Escanilla D., Franzen C., Uhlen M., Albert J., Proc. Natl. Acad. Sci. U.S.A. 93, 10864 (1996);

Leitner T., Kumar S., Albert J., J. Virol. 71, 4761 (1997) ;

Leitner T., Albert J., Proc. Natl. Acad. Sci. U.S.A. 96, 10752 (1999).

Sharp P. M., Nature 336, 315 (1988);

Smith T. F., Srinivasan A., Schochetman G., Marcus M., Myers G., Nature 333, 573 (1988);

Kasper P., et al., AIDS Res. Hum. Retroviruses 11, 1197 (1995).

10.1126/science.280.5371.1868

10.1007/BF01734359

; Cladistics 5 164 (1989);

Olsen G. J., Matsuda H., Hagstrom R., Overbeek R., Comput. Appl. Biosci. 10, 41 (1994).

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 /).

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 1996) pp. 407–514.

D. M. Hillis B. K. Mable C. Moritz in Molecular Systematics 2nd ed. D. M. Hillis C. Moritz B. K. Mable Eds. (Sinauer Sunderland MA (1996) pp. 515–543.

10.1093/oxfordjournals.molbev.a025892

The alignment was based on those provided in Human Retroviruses and AIDS B. Korber et al. Eds. (Los Alamos National Laboratory Los Alamos NM 1998).

A complete description of the alignments details of the phylogenetic analysis the sequence alignments and the links new code written for this study are provided at www.santafe.edu/btk/science-paper/bette.html.

Shankarappa R., et al., J. Virol. 73, 10489 (1999).

10.1126/science.272.5261.537

; S. Ganeshan et al.

Wolinsky S., J. Virol. 71, 663 (1997);

Markham R. B., et al., Proc. Natl. Acad. Sci. U.S.A. 95, 12568 (1998);

Bagnarelli P., et al., J. Virol. 73, 3764 (1999).

Furtado M. R., et al., J. Virol. 69, 2092 (1995);

10.1126/science.271.5255.1582

10.1126/science.278.5341.1295

10.1038/8394

Gunthard H., et al., J. Virol. 73, 9404 (1999);

10.1056/NEJM199905273402101

Furtado M. R., N. Engl. J. Med. 340, 1614 (1999).

10.1126/science.253.5018.390

Grassly N., Harvery P., Holmes E., Genetics 151, 427 (1999).

Robertson D. L., Hahn B. H., Sharp P. M., Mol. Evol. 40, 249 (1995);

Heyndrickx L., et al., J. Virol. 74, 363 (2000).

Kuiken C. L., et al., AIDS 10, 31 (1996).

Salemi M., et al., Proc. Natl. Acad. Sci. U.S.A. 96, 13253 (1999);

Smith P., Simmons D., J. Virol. 73, 5787 (1999);

Bollyky P., Holmes E., Mol. Evol. 49, 130 (1999).

Yang Z., J. Mol. Evol. 39, 105 (1994);

; J. Mol. Evol. 39 306 (1994); J. Mol. Evol. 42 587 (1996).

10.1126/science.8171318

Yang Z., Goldman N., Friday A., Mol. Biol. Evol. 11, 316 (1994);

Huelsenbeck J., Mol. Biol. Evol. 12, 843 (1995);

Kuhner M., Felsenstein J., Mol. Biol. Evol. 11, 459 (1994).

Yang Z., Genetics 139, 993 (1995).

10.1126/science.276.5310.227

Because it is not possible to test all possible tree configurations tree-building programs use heuristics to estimate the best tree and the final tree is dependent on the input order of sequences. To optimize the final trees we randomized the input order of the sequences five to seven times until the best maximum-likelihood scores were very similar (1). Given the number of taxa we included and consequently the combinatorially vast potential for different branching orders we do not expect our trees to be optimal solutions. Limited testing of the final timing estimates based on different input orders of sequences did not significantly affect our calculations of the timing of divergence from a common ancestor. We also compared the likelihood of the data under different evolutionary models (1) and over 100 maximum-likelihood trees were run in the course of this study.

Yang Z., Mol. Biol. Evol. 10, 1396 (1993).

D. L. Swofford PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods) (Sinauer Sunderland MA 1999).

We also tested other aspects of our evolutionary model. We found that the assignment of base frequencies by means of the phylogenetic trees gave consistently better results than empirical base frequencies. The REV model performed better than an F84 model (2) which only includes rate parameters for transitions and transversions instead of for each pair of bases. Also for the envelope gene analyses the improvement in the log-likelihood scores comparing the REV model with a uniform rate of evolution at all sites to the REV model with rate variation at different sites estimated by the maximum-likelihood method was many times larger than the number of positions (1) justifying the increase in parameters (3).

Korber B., Sharp P., Ho D., Nature 400, 326 (1999);

Goudsmit J., Lukoshov V., Nature 400, 325 (1999).

B. Korber et al. data not shown.

Courgnaud V., et al., Virology 247, 41 (1998).

Wangroongsarb Y., et al., Southeast Asian J. Trop. Med. Public Health 16, 517 (1985);

Smith D., Lancet 335, 781 (1990);

Mason C., et al., J. Acquir. Immune Defic. Syndr. Hum. Retrovirol. 19, 165 (1996);

Bunnell R., et al., AIDS 13, 509 (1999).

C. Kuiken et al. Am. J. Epidemiol. in press.

10.1128/jvi.70.6.3331-3338.1996

Subbarao S., et al., AIDS Res. Human Retroviruses 14, 319 (1998).

Gao F., et al., J. Virol. 70, 7013 (1996);

Carr J. K., et al., J. Virol. 70, 5935 (1996).

Gottlieb M. S., et al., N. Engl. J. Med. 305, 1425 (1981);

Gottlieb M. S., et al., Morb. Mortal. Wkly. Rep. 30, 250 (1981).

Selik R., Haverkos H., Curran J., Am. J. Med. 76, 493 (1984);

Pape J., N. Engl. J. Med. 309, 945 (1983);

Gazzolo L., N. Engl. J. Med. 311, 1252 (1984).

E. Hooper. The River (Little Brown Boston 1999). See pp. 77–82 and 440–443 for discussion of early cases in the United States and Haiti and pp. 550 791 and 1009 for a discussion of the number of primate kidneys required to make OPV.

Chevret S., et al., J. Epidemiol. Commun. Health 46, 582 (1992).

Li W.-H., Tanimura M., Sharp P., Mol. Biol. Evol. 5, 313 (1988);

; T. Gojobori et al. Proc. Natl. Acad. Sci. U.S.A. 340 1605 4108 (1990); J. Kelly Genet. Res. 64 1 1994.

Liitsola K., et al., AIDS 12, 1907 (1998).

A record of the ages of chimpanzees from Camp Lindi used for research noted a range from <1 to 10 years with more than 80% less than 4 years old (S. Plotkin personal communication; data taken from the laboratory notes of F. Deinhardt).

M. Grmek. History of AIDS Emergence and Origin of a Modern Pandemic (Princeton Univ. Press Princeton NJ 1990) chaps. 10 and 15.

10.1089/088922200309548

We thank D. Pollock T. Leitner and B. Bruno for suggestions concerning phylogenetics maximum likelihood and estimating the error on time of sampling; G. Shaw for suggesting the 1959 control; S. Wain-Hobson and G. Myers for clarifying discussions on the interpretation and limitations of these results; B. Foley and C. Kuiken for numerous helpful discussions; and K. Rock and J. Shepard for technical support. G. Olsen and J. Thorne generously supplied source code and helped us interpret their work. The research of the Los Alamos authors was supported under internal funds from the Delphi Project S.W. and B.K. were supported by NIH (RO1-HD37356) B.K. and M.M. were supported through the Pediatric AIDS Foundation and an anonymous foundation supplied further support for S.W. B.H.H. was supported by grants NO1 AI 85338 RO1 AI 44596 and RO1 AI 40951 from NIH.