HIV-1 infected monozygotic twins: a tale of two outcomes
© Tazi et al; licensee BioMed Central Ltd. 2011
Received: 16 December 2009
Accepted: 8 March 2011
Published: 8 March 2011
Replicate experiments are often difficult to find in evolutionary biology, as this field is inherently an historical science. However, viruses, bacteria and phages provide opportunities to study evolution in both natural and experimental contexts, due to their accelerated rates of evolution and short generation times. Here we investigate HIV-1 evolution by using a natural model represented by monozygotic twins infected synchronically at birth with an HIV-1 population from a shared blood transfusion source. We explore the evolutionary processes and population dynamics that shape viral diversity of HIV in these monozygotic twins.
Despite the identical host genetic backdrop of monozygotic twins and the identical source and timing of the HIV-1 inoculation, the resulting HIV populations differed in genetic diversity, growth rate, recombination rate, and selection pressure between the two infected twins.
Our study shows that the outcome of evolution is strikingly different between these two "replicates" of viral evolution. Given the identical starting points at infection, our results support the impact of random epigenetic selection in early infection dynamics. Our data also emphasize the need for a better understanding of the impact of host-virus interactions in viral evolution.
Over the past few decades, evolutionary biology has had an increasing impact on biomedical research [1–3]. Evolutionary theory can address pertinent questions related to the control of infectious diseases and particularly to pathogen virulence . RNA viruses serve as exciting models for testing this theory because of their potential for rapid evolution. For instance, a striking feature of the Human Immunodeficiency Virus 1 (HIV-1) is the rapid population dynamics resulting in a high degree of genetic diversity within and between infected individuals [5, 6]. The virus can then capitalize on this genetic diversity to evade a host immune response [7, 8]. Nevertheless, host-virus coevolution is also important with respect to disease progression [9, 10]. Viral evolution can be strongly shaped through antiviral pressure applied by the host's immune system [11, 12].
A few complications in the study of host-virus coevolution are the different genetic backgrounds of the host and the genetic diversity of the infecting viral population. The model of monozygotic twins infected synchronically through exposure to the same viral population, therefore, provides an opportunity to examine the population dynamics of HIV while holding a number of confounding variables constant. Here we exploit HIV-1 as a model system  to examine the population genetic processes and epigenetic influence on viral evolution. We collected nucleotide sequence data from the genes Protease (pro), Reverse Transcriptase (rt), and Envelope (env) in each monozygotic twin. Using these data, we estimated relevant population genetic parameters, including recombination rates, genetic diversity, growth rates, and selection pressure to compare between the two host individuals. While these twins have identical genetic backgrounds and were infected with the same source of HIV (a common transfusion source collected from one donor and administered simultaneously to both twins at birth), they have remarkable differences in their clinical courses. Twin A is almost normal in terms of immune system (total average of CD4 T cell counts: 860 cells/μl for twin A; and 319 cells/μl for twin B) and growth, whereas twin B is almost 5 years delayed in terms of size and had many complications including a rare neoplasm.
Recent studies focusing on monozygotic twins are conflicting on whether HIV evolution can be predictable [14, 15] or unpredictable [16, 17], related to the immunological repertoire recruitment in these hosts. Those studies examined immune responses, viral evolution, and disease outcome in monozygotic twins infected simultaneously with the same virus. They showed that immune selection driven by dominant sequences in each host could contribute to specific pathways of HIV-1 evolution. The same natural model is also investigated here; however we focus on the evolutionary processes acting on HIV-1 evolution in identical twins by characterizing the population genetic parameters, including positive selection in the viral population. If the viral evolution is predictable, we expect to see similar population dynamics in both twins. However, if there is significant epigenetic impact on viral evolution, then we would see differences in the population dynamics and associated population genetic parameters.
Results and Discussion
Estimates of genetic diversity (θ), recombination (r and C), and growth (g) for pro, rt and env in each twin.
0.44 [0.23 - 0.92]
0.07 [0.02 - 0.17]
0.03 [0.005 - 0.16]
324 [222 - 426]
0.26 [0.14 - 0.6]
0.23 [0.1 - 0.38]
0.06 [0.01 - 0.23]
177 [86 - 277]
0.2 [0.13 - 0.33]
0.03 [0.01 - 0.07]
0.01 [0.001 - 0.02]
42 [19 - 66]
0.2 [0.09 - 0.71]
0.1 [0.003 - 0.47]
0.02 [0.0003 - 0.33]
175 [65 - 480]
0.12 [0.07 - 0.24]
0.14 [0.04 - 0.25]
0.02 [0.003 - 0.06]
81 [27 - 165]
0.1 [0.06 - 0.16]
0.02 [0.01 - 0.07]
0.002 [0.001 - 0.01]
13 [(-18) - 41]
Log-likelihood values and parameter estimates (ω, p, and n) for pro, rt and env in each twin.
Selection analyses for pro, rt and env in each twin.
Our goal was to investigate HIV-1 evolution in identical twins infected synchronically at birth with the same blood transfusion. We found compelling distinctions between the viral populations from each twin with respect to their population dynamics, phylogenetic structure, growth rates, recombination rates, genetic diversity, and selection pressures. These results were unexpected due to a combination of having identical starting points with respect to both the infecting viral population and the host genetic background. That coupled with the seemingly limited pathways of evolution for both immune evasion and evolution of drug resistance [14, 23] would lead one to predict similar patterns of genetic diversity and dynamics in viral populations resulting in similar clinical outcomes. Instead, we found higher growth rates, higher genetic diversity, and higher recombination in rt in the healthier twin A compared to twin B. We also found sites under diversifying selection in env in twin A whereas twin B had only one site under selection in rt (in PAML analysis). Thus, the higher genetic diversity and higher number of selected sites in env appear to be associated with slower disease progression, results concordant with that found in a broad study of disease progression in infants . Similarly, the twins differed in their population dynamics and these differed by gene region. The rt and env regions showed the viral population in the healthier twin A with higher levels of genetic diversity throughout the history of infection even when there were significant shifts in overall levels of diversity. On the other hand, pro showed the viral population in twin B with a gradual increase in diversity post drug therapy with a more rapid increase in twin A that was delayed by the same time period as the delay in the RT inhibitor (2 years). This result suggests that the shape of the response to drug therapy in terms of the HIV population diversity might be diagnostic of future disease progression, but further study with larger sample sizes are needed to better test this response as predictive of disease progression. Nevertheless, these twins clearly show very different responses to infection.
This difference in viral population dynamics is concordant with the differences observed in the clinical courses in each twin. The immune system in twin A shows CD4 T cells at an almost normal rate. The immune system in twin B is depressed, hence no strong selective pressure is acting upon its virus population to evolve fast . All these results combined provide strong evidence that, at least in this case, the replicate evolutionary experiment did not result in an identical outcome demonstrating the importance of selective response to random epigenetic factors impacting disease progression .
Indeed, some studies in monozygotic twins revealed increasing epigenetic differences with age [25, 26]. Additionally, there is a clear potential for founding effects upon infection , even in the context of a blood transfusion as the viral population in a shared blood donation is certainly reduced in genetic diversity and number compared to infectious virus from an infected individual. The combination of a reduced effective population size coupled with strong selective pressure is a key ingredient for founder effects [28, 29], resulting in populations with very different characteristics as evident here in both the population dynamics and immunology. Clearly, the early impact of founder effects and epigenetic factors on viral population dynamics has diversifying impact over time as the viral populations undergo independent evolution - even in the face of similar genetic selection pressures, identical genetic starting points, and identical host genetic backgrounds (Figure 1). This epigenetic drift during development can be either stochastic (especially when impacted by genetic drift) or determined by environmental factors . Host-virus interactions in early stages HIV infection are presumed to have a large impact on the disease course and viral evolution [31–34], yet they are exceptionally difficult to study because researchers are typically not able to design experiments to investigate viral dynamics at infection. Our study capitalizes on the infection of monozygotic twins through a common contaminated blood transfusion to demonstrate that even more complicated epigenetic factors need to be taken into account in developing hypotheses associated with genetic diversity, population dynamics, selection pressure and their association with disease progression.
We used monozygotic twins infected at birth from the same blood transfusion contaminated with HIV-1 to study the association of population genetic and phylogenetic diversity with disease progression and clinical outcomes. We documented that these twins had very different clinical outcomes with twin A being relatively healthy compared to twin B. Associated with this slower disease progression in twin A, we found phylogenetic differences, higher growth rates, and higher genetic diversity in the HIV population and higher recombination rates in rt. We also found differences in population dynamics across all three gene regions. These differences suggest that epigenetic factors are important in disease progression and can impact viral genetic population dynamics.
Subjects A and B are monozygotic twins. Due to intrapartum blood loss, they simultaneously received transfusion shortly after birth from a common blood transfusion source in 1983. The parents were notified that the twins were infected with HIV when they were 2 years old and the twins were declared HIV positive in the same year (1985). At that time, they were both asymptomatic. No superinfection was reported during the first five years of age of the twins. Also the twins did not receive additional transfusions since their birth. The twins came to NIH when they were 5 years old (1988). Since then, they have been closely monitored at NIH. The twins started receiving nucleoside analog reverse transcriptase inhibitors at ages 5 (Twin B) and 7 (Twin A). In 1988, twin B started AZT therapy because he was sicker than his brother. Twin A started ddI therapy two years later (1990). Twin A is normal in terms of immune system and growth. In contrast, twin B is 5 years delayed in terms of size and he presented a rare neoplasm, more likely due to HIV infection. Blood samples and clinical information were obtained 15 years after their birth at the NIH Clinical Center, Bethesda, MD (1998). The twins subsequently showed extremely different clinical courses. Samples were collected under protocols approved by the institutional review boards at the National Cancer Institute. Written informed consent was obtained from participants. This research complied with all relevant federal guidelines and institutional policies. Permission was obtained from the parents of the patients for utilization of samples and data for research, analysis and publication.
Peripheral blood mononuclear cells (PBMCs) were isolated from blood samples by Ficoll-Hypaque gradient density at 650 g for 30 min, and the cell pellets were stored in liquid nitrogen until use. DNA was extracted from PBMCs using the PureGene genomic DNA isolation kit (Gentra Systems, Minneapolis, MN) and prepared at optimal concentration for nested polymerase chain reaction (PCR) method. Nested PCR amplification of the V1-V5 region of envelope coding gene, the protease gene, and a part of reverse transcriptase gene was performed with the Expand High Fidelity PCR System (Roche Applied Science, Indianapolis, IN) in a 50 μ l reaction containing 1× Expand High Fidelity buffer 3, 200 μ M dNTPs, 2 mM MgCl2, 400 nM primers, and 1.75 U of Expand High Fidelity PCR System enzyme mix. Primer sets used to amplify a 1.1-Kbp fragment, encompassing the env V1-V5 region were: +6559 (sense) 5'-GGGATCAAAGCCTAAAGCCA-3' and -7648 (antisense) 5'- ACTTCTCCAATTGTCCCTCA-3' in a first round reaction; +6586 (sense) 5'- AATTAACCCCACTCTGTGTTA-3' and -7627 (antisense) 5'- TATCTCCTCCTCCAGGTCTGA-3' in a second round reaction. Amplification of the HIV-1 pro was done with the following primer sets: +2165 (sense) 5'- CAGAAGAGAGCTTCAGGTTTGGG-3' and -2588 (antisense) 5'- ACTTTTGGGCCATCCATTCCTGGY-3' in a first round reaction; +2208 (sense) 5'-TCAGAAGCAGGAGCCGATAGAC-3' and -2550 (antisense) 5'-TGGTACAGTCTCAATAGGACTAATGGG-3' in a second round reaction. The part of rt was amplified with: +1882 (sense) 5'-GAAGCAATGAGCCAAGTAACAAAT-3' and -3544 (antisense) 5'-GATATGTCCATTGGCCTTGCCCCT-3' in a first round reaction; +1966 (sense) 5'-TTCAATTGTGGCAAAGAAGGGCAC-3' and -3501 (antisense) 5'-TAAGTCTTTTGATGGGTCATAATA-3' in a second round reaction. Each round of PCR consisted of 25 cycles, with the initial denaturation at 94°C for 2 min, followed by 25 cycles of denaturation at 94°C for 15 s, annealing at 50°C for 30 s, and extension at 72°C for 1 min, with the final extension at 72°C for 7 min. The PCR products were purified with the QIA quick PCR purification kit (QIAGEN, Valencia, CA), and then cloned into pCR2.1-TOPO vector using the TOPO TA Cloning Kit (Invitrogen, Carlsbad, CA). The clones for each gene were sequenced using the ABI BigDye Terminator v3.1 Ready Reaction Cycle Sequencing Kit (Applied Biosystems, Foster City, CA). Sequences were then analyzed on an ABI PRISM 3130 × l Genetic Analyzer (Applied Biosystems, Foster City, CA). GenBank accession numbers for the sequences determined in this study are GQ118464 to GQ118640.
Phylogenetic analysis and network reconstruction
Sequences from each twin were aligned using MAFFT v5.3  and manually edited using MacClade 4.05 . Phylogenetic trees were estimated using the maximum likelihood approach  with nodal support assessed via bootstrapping (1,000 pseudoreplicates) , as implemented in PHYML . We also estimated phylogenies using Bayesian methods  coupled with Markov Chain Monte Carlo (BMCMC) inference as implemented in MrBayes v3.1.2 . Model selection for these analyses followed the procedure outlined by Posada and Buckley  and implemented in ModelTest v3.6 (using PAUP*)  under the Akaike Information Criterion (AIC) . Two independent BMCMC analyses were run, each consisting of four MCMC chains (one cold and three heated). Each Markov chain started from a random tree and ran for 2.0 × 107 cycles, sampling every 1000th generation. In order to confirm that the chains converged and mixed well, we monitored the likelihood scores and compared means and variances of all likelihood parameters from independent runs using Tracer v1.4.1 . The relatedness of the twins' infections was assessed by adding 20 HIV-1 reference sequences to the phylogenetic analysis: 11 lab reference strains (HIVU455, HIVELI, HIVSC, HIVHXB2, HIVJRFL, HIVJRCSF, HIVD31, HIVBRVA, HIVMN, HIVNY5 and HIVSF2) and the nine closest sequence matches (Sequence identity: 85%-95%) identified by BLAST analyses (DQ410259, AF025761, AY247221, DQ313246, AY247218, AY835759, EU184128, HQ110637 and DQ410044). Evolutionary relationships among the sequences in each twin were also assessed using the method of statistical parsimony , as implemented in the software package TCS v1.21 . This approach allows for the visualization of evolutionary relationships as a network instead of a bifurcating tree which is often more appropriate for viral populations that recombine.
Population genetic parameters
Genetic diversity (θ = 4N e μ; where N e is the effective population size and μ is the mutation rate per site), recombination rates (r = c/μ; where c is the recombination chance per site), and growth rates (g) were estimated for each set of genes using the maximum likelihood coalescent approach implemented in LAMARC v2.0.2 . Three independent runs were performed for each gene in order to assess the reproducibility of the different parameter estimates. Furthermore, an estimate of the diversity generated per recombination (C) was obtained by multiplying r and θ .
HIV-1 past population dynamics in each twin was inferred using the Bayesian skyline plot  and relaxed clock (lognormal) models  in BEAST v1.4.2 . The rate of substitution was calibrated based on the age of infection (15 years). The hyperparameter m (number of grouped intervals) was set up to 1/4 of the sequences in each dataset. All the Bayesian MCMC output generated by BEAST was then analyzed in Tracer v1.4.1  to test for convergence and mixing and plot population demographics.
The extent of natural selection was inferred by estimating the ratio of nonsynonymous to synonymous substitutions (ω = dN/dS) per site and per gene using the codon-based nested models M1 (neutral)/M2 (selection), and M7 (beta)/M8 (beta and ω), as implemented in PAML v4 . Model likelihoods were compared using a Likelihood Ratio Test (LRT) to determine the best-fit model. The Bayes empirical Bayes approach was applied to identify the potential sites under diversifying selection as indicated by a posterior probability (pP > 0.95) . In order to confirm the PAML analysis, additional ω estimates were also calculated in DataMonkey . Positively and negatively selected sites were identified using the REL approach . Recombination was taken into account by screening for recombination breakpoints using GARD and allowing each recombinant fragment to have its own phylogenetic tree .
The authors would like to thank study participants for their help in completing this study. This work was supported by NIH grant R01-GM66276 (KAC, DP) and Brigham Young University, and by the Spanish Ministry of Science and Education grant number BIO2007-61411 to DP. This project has been funded in part with federal funds from the National Cancer Institute, National Institutes of Health, under Contract No. HHSN261200800001E and by the Intramural Research Program of the National Institute of Allergy and Infectious Diseases. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.
- Ewald PW, Sussman JB, Distler MT, Libel C, Chammas WP, Dirita VJ, Salles CA, Vicente AC, Heitmann I, Cabello F: Evolutionary control of infectious disease: prospects for vectorborne and waterborne pathogens. Mem Inst Oswaldo Cruz. 1998, 93: 567-576. 10.1590/S0074-02761998000500002.View ArticlePubMedGoogle Scholar
- Gammelgaard A: Evolutionary biology and the concept of disease. Med Health Care Philos. 2000, 3: 109-116. 10.1023/A:1009999502884.View ArticlePubMedGoogle Scholar
- Ellison PT: Evolutionary biology for doctors. Science. 2009, 325: 1207-10.1126/science.1179152.View ArticleGoogle Scholar
- Frank SA, Schmid-Hempel P: Mechanisms of pathogenesis and the evolution of parasite virulence. J Evol Biol. 2008, 21: 396-404. 10.1111/j.1420-9101.2007.01480.x.View ArticlePubMedGoogle Scholar
- Coffin JM: Genetic diversity and evolution of retroviruses. Curr Top Microbiol Immunol. 1992, 176: 143-164.PubMedGoogle Scholar
- Wei X, Ghosh SK, Taylor ME, Johnson VA, Emini EA, Deutsch P, Lifson JD, Bonhoeffer S, Nowak MA, Hahn BH, et al: Viral dynamics in human immunodeficiency virus type 1 infection. Nature. 1995, 373: 117-122. 10.1038/373117a0.View ArticlePubMedGoogle Scholar
- Goulder PJ, Watkins DI: HIV and SIV CTL escape: implications for vaccine design. Nat Rev Immunol. 2004, 4: 630-640. 10.1038/nri1417.View ArticlePubMedGoogle Scholar
- Fauci AS, Johnston MI, Dieffenbach CW, Burton DR, Hammer SM, Hoxie JA, Martin M, Overbaugh J, Watkins DI, Mahmoud A, et al: HIV vaccine research: the way forward. Science. 2008, 321: 530-532. 10.1126/science.1161000.View ArticlePubMedGoogle Scholar
- Pantaleo G, Fauci AS: Immunopathogenesis of HIV infection. Annu Rev Microbiol. 1996, 50: 825-854. 10.1146/annurev.micro.50.1.825.View ArticlePubMedGoogle Scholar
- Letvin NL, Walker BD: Immunopathogenesis and immunotherapy in AIDS virus infections. Nat Med. 2003, 9: 861-866. 10.1038/nm0703-861.View ArticlePubMedGoogle Scholar
- Moore CB, John M, James IR, Christiansen FT, Witt CS, Mallal SA: Evidence of HIV-1 adaptation to HLA-restricted immune responses at a population level. Science. 2002, 296: 1439-1443. 10.1126/science.1069660.View ArticlePubMedGoogle Scholar
- Yang OO, Sarkis PT, Ali A, Harlow JD, Brander C, Kalams SA, Walker BD: Determinant of HIV-1 mutational escape from cytotoxic T lymphocytes. J Exp Med. 2003, 197: 1365-1375. 10.1084/jem.20022138.View ArticlePubMedPubMed CentralGoogle Scholar
- Crandall KA: The Evolution of HIV. 1999, Baltimore, MD: The Johns Hopkins University PressGoogle Scholar
- Draenert R, Allen TM, Liu Y, Wrin T, Chappey C, Verrill CL, Sirera G, Eldridge RL, Lahaie MP, Ruiz L, et al: Constraints on HIV-1 evolution and immunodominance revealed in monozygotic adult twins infected with the same virus. J Exp Med. 2006, 203: 529-539. 10.1084/jem.20052116.View ArticlePubMedPubMed CentralGoogle Scholar
- O'Connor DH, Burton DR: Immune responses and HIV: a little order from the chaos. J Exp Med. 2006, 203: 501-503.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang OO, Church J, Kitchen CM, Kilpatrick R, Ali A, Geng Y, Killian MS, Sabado RL, Ng H, Suen J, et al: Genetic and stochastic influences on the interaction of human immunodeficiency virus type 1 and cytotoxic T lymphocytes in identical twins. J Virol. 2005, 79: 15368-15375. 10.1128/JVI.79.24.15368-15375.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu XG, Lichterfeld M, Williams KL, Martinez-Picado J, Walker BD: Random T-cell receptor recruitment in human immunodeficiency virus type 1 (HIV-1)-specific CD8+ T cells from genetically identical twins infected with the same HIV-1 strain. J Virol. 2007, 81: 12666-12669. 10.1128/JVI.01450-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Posada D, Crandall KA, Holmes EC: Recombination in evolutionary genomics. Annu Rev Genet. 2002, 36: 75-97. 10.1146/annurev.genet.36.040202.111115.View ArticlePubMedGoogle Scholar
- Schinazi RF, Larder BA, Mellors JW: Mutations in retroviral genes associated with drug resistance. Int Antiviral News. 1997, 5: 129-142.Google Scholar
- Nozawa M, Suzuki Y, Nei M: Reliabilities of identifying positive selection by the branch-site and the site-prediction methods. Proc Natl Acad Sci USA. 2009, 106: 6700-6705. 10.1073/pnas.0901855106.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Z, Nielsen R, Goldman N: In defense of statistical methods for detecting positive selection. Proc Natl Acad Sci USA. 2009, 106: E95-10.1073/pnas.0904550106.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Crandall KA, Kelsey CR, Imamichi H, Salzman NP: Parallel evolution of drug resistance in HIV: failure of nonsynonymous/synonymous substitution rate ratio to detect selection. Mol Biol Evol. 1999, 16: 372-382.View ArticlePubMedGoogle Scholar
- Carvajal-Rodriguez A, Posada D, Perez-Losada M, Keller E, Abrams EJ, Viscidi RP, Crandall KA: Disease progression and evolution of the HIV-1 env gene in 24 infected infants. Infect Genet Evol. 2008, 8: 110-120. 10.1016/j.meegid.2007.10.009.View ArticlePubMedGoogle Scholar
- Fraga MF, Ballestar E, Paz MF, Ropero S, Setien F, Ballestar ML, Heine-Suner D, Cigudosa JC, Urioste M, Benitez J, et al: Epigenetic differences arise during the lifetime of monozygotic twins. Proc Natl Acad Sci USA. 2005, 102: 10604-10609. 10.1073/pnas.0500398102.View ArticlePubMedPubMed CentralGoogle Scholar
- Wong AH, Gottesman II, Petronis A: Phenotypic differences in genetically identical organisms: the epigenetic perspective. Hum Mol Genet. 2005, 14: 11-18. 10.1093/hmg/ddi116.View ArticleGoogle Scholar
- Leigh Brown AJ: Analysis of HIV-1 env gene sequences reveals evidence for a low effective number in the viral population. Proc Natl Acad Sci USA. 1997, 94: 1862-1865. 10.1073/pnas.94.5.1862.View ArticleGoogle Scholar
- Templeton AR: Modes of speciation and inferences based on genetic distances. Evolution. 1980, 34: 719-729. 10.2307/2408026.View ArticleGoogle Scholar
- Carson HL: Increased genetic variance after a population bottleneck. Trends Ecol Evol. 1990, 5: 228-230. 10.1016/0169-5347(90)90137-3.View ArticlePubMedGoogle Scholar
- Poulsen P, Esteller M, Vaag A, Fraga MF: The epigenetic basis of twin discordance in age-related diseases. Pediatr Res. 2007, 61: 38-42. 10.1203/pdr.0b013e31803c7b98.View ArticleGoogle Scholar
- Keele BF, Giorgi EE, Salazar-Gonzalez JF, Decker JM, Pham KT, Salazar MG, Sun C, Grayson T, Wang S, Li H, et al: Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proc Natl Acad Sci USA. 2008, 105: 7552-7557. 10.1073/pnas.0802203105.View ArticlePubMedPubMed CentralGoogle Scholar
- Hay C, Rosenberg E: Immunologic response to HIV. AIDS Clin Care. 1998, 10: 1-3.PubMedGoogle Scholar
- Streeck H, Brumme ZL, Anastario M, Cohen KW, Jolin JS, Meier A, Brumme CJ, Rosenberg ES, Alter G, Allen TM, et al: Antigen load and viral sequence diversification determine the functional profile of HIV-1-specific CD8+ T cells. PLoS Med. 2008, 5: e100-10.1371/journal.pmed.0050100.View ArticlePubMedPubMed CentralGoogle Scholar
- Brumme ZL, Brumme CJ, Carlson J, Streeck H, John M, Eichbaum Q, Block BL, Baker B, Kadie C, Markowitz M, et al: Marked epitope- and allele-specific differences in rates of mutation in human immunodeficiency type 1 (HIV-1) Gag, Pol, and Nef cytotoxic T-lymphocyte epitopes in acute/early HIV-1 infection. J Virol. 2008, 82: 9216-9227. 10.1128/JVI.01041-08.View ArticlePubMedPubMed CentralGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.View ArticlePubMedPubMed CentralGoogle Scholar
- Maddison DR, Maddison WP: MacClade 4: Analysis of Phylogeny and Character Evolution, version 4.0. 2000, Sinauer AssociatesGoogle Scholar
- Felsenstein J: Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981, 17: 368-376. 10.1007/BF01734359.View ArticlePubMedGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Larget B, Miller RE, Ronquist F: Potential applications and pitfalls of Bayesian inference of phylogeny. Syst Biol. 2002, 51: 673-688. 10.1080/10635150290102366.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Posada D, Buckley TR: Model selection and model averaging in phylogenetics: Advantages of akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst Biol. 2004, 53: 793-808. 10.1080/10635150490522304.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Hirotugu A: A new look at the statistical model identification. IEEE Transact Automatic Control. 1974, 19: 716-723. 10.1109/TAC.1974.1100705.View ArticleGoogle Scholar
- Rambaut A, Drummond AJ: Tracer: MCMC trace analysis tool. 1.4.1 edition. 2009, Edinburgh: Institute of Evolutionary Biology, [http://tree.bio.ed.ac.uk/software/tracer/]Google Scholar
- Templeton AR, Crandall KA, Sing CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132: 619-633.PubMedPubMed CentralGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMedGoogle Scholar
- Kuhner MK: LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006, 22: 768-770. 10.1093/bioinformatics/btk051.View ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005, 22: 1185-1192. 10.1093/molbev/msi103.View ArticlePubMedGoogle Scholar
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.View ArticlePubMedPubMed CentralGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.View ArticlePubMedPubMed CentralGoogle Scholar
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148: 929-936.PubMedPubMed CentralGoogle Scholar
- Kosakovsky Pond SL, Frost SDW: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, 21: 2531-2533. 10.1093/bioinformatics/bti320.View ArticleGoogle Scholar
- Kosakovsky Pond SL, Frost SDW: Not So Different After All: A Comparison of Methods for Detecting Amino Acid Sites Under Selection. Mol Biol Evol. 2005, 22: 1208-1222. 10.1093/molbev/msi105.View ArticlePubMedGoogle Scholar
- Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SDW: Automated Phylogenetic Detection of Recombination Using a Genetic Algorithm. Mol Biol Evol. 2006, 23: 1891-1901. 10.1093/molbev/msl051.View ArticlePubMedGoogle Scholar