Phylodynamics of HIV-1 Circulating Recombinant Forms 12_BF and 38_BF in Argentina and Uruguay
© Bello et al; licensee BioMed Central Ltd. 2010
Received: 21 December 2009
Accepted: 22 March 2010
Published: 22 March 2010
Although HIV-1 CRF12_BF and CRF38_BF are two epidemiologically important recombinant lineages circulating in Argentina and Uruguay, little is known about their population dynamics.
A total of 120 "CRF12_BF-like" and 20 "CRF38_BF-like" pol recombinant sequences collected in Argentina and Uruguay from 1997 to 2009 were subjected to phylogenetic and Bayesian coalescent-based analyses to estimate evolutionary and demographic parameters.
Phylogenetic analyses revealed that CRF12_BF viruses from Argentina and Uruguay constitute a single epidemic with multiple genetic exchanges among countries; whereas circulation of the CRF38_BF seems to be confined to Uruguay. The mean estimated substitution rate of CRF12_BF at pol gene (2.5 × 10-3 substitutions/site/year) was similar to that previously described for subtype B. According to our estimates, CRF12_BF and CRF38_BF originated at 1983 (1978-1988) and 1986 (1981-1990), respectively. After their emergence, the CRF12_BF and CRF38_BF epidemics seem to have experienced a period of rapid expansion with initial growth rates of around 1.2 year-1 and 0.9 year-1, respectively. Later, the rate of spread of these CRFs_BF seems to have slowed down since the mid-1990s.
Our results suggest that CRF12_BF and CRF38_BF viruses were generated during the 1980s, shortly after the estimated introduction of subtype F1 in South America (~1975-1980). After an initial phase of fast exponential expansion, the rate of spread of both CRFs_BF epidemics seems to have slowed down, thereby following a demographic pattern that resembles those previously reported for the HIV-1 epidemics in Brazil, USA, and Western Europe.
The AIDS epidemic in South America is caused by multiple HIV-1 group M subtypes including subtypes B, F1, and C, in addition to BF1 and BC recombinant forms. The BF1 recombinants represent the most widespread genetic form after subtype B and reach a high prevalence (10%-50%) in countries from the Southern Cone (Argentina, Brazil, Chile, Paraguay, and Uruguay) [1–14].
Genetic characterization of BF1 recombinants in South America revealed some important differences across countries. Although four distinct BF1 circulating recombinant forms (CRFs) have been described in Brazil to date (CRF28_BF, CRF29_BF, CRF39_BF, and CRF40_BF) [15, 16], the Brazilian BF1 epidemic is largely dominated by a variety of unique recombinants forms (URFs) that do not share a common recombinant ancestor [7, 10, 17–19]. In contrast, the Argentine BF1 epidemic comprises the widespread CRF12_BF and several URFs with a CRF12-related structure [6, 8, 20, 21]. The molecular epidemiology of HIV-1 in Uruguay is not so well characterized, but two previous studies suggested that BF1 recombinants circulating in this country are similar to those described in Argentina [3, 20]. Very recently, a novel CRF38_BF1 was described among Uruguayan HIV-1 isolates, indicating that other BF1 recombinants besides CRF12_BF have gained epidemic importance in this country .
Previous studies performed by our group suggest that the HIV-1 subtype F1 and BF1 epidemics in South America were initiated after the introduction of a single F1 strain into Brazil between the middle and late 1970s [22–24]. After its introduction, this founder subtype F1 strain probably recombined with the local subtype B virus generating the large diversity of CRFs_BF1 and URFs_BF1 currently observed in the continent . Based on monophyletic clustering and coincident recombination breakpoints, it was suggested that most BF1 recombinants circulating in Argentina and Uruguay derived from a common recombinant ancestor [21, 24, 25].
To date, however, very little is known about the evolutionary history and epidemic potential of the diverse BF1 recombinants that have expanded in the South American population. Only one previous study was conducted on a small number (n = 40) of CRF12_BF-like vpu sequences from a vertically infected population in Argentina . This study estimated the age of the most recent common ancestor (MRCA) of those CRF12_BF-like viruses between 1981 and 1996, and further suggests an extremely rapid spread of the CRF12_BF-like recombinant viruses, compatible with the demographic pattern of explosive population growth observed in this pediatric population at the start of the epidemic.
The objective of the present study was to reconstruct the evolutionary and demographic history of the CRF12_BF circulating in Argentina and Uruguay through the analysis of a large data set (n = 120) of CRF12_BF-like pol sequences recovered from adults and children living in both countries. In addition, we also analyzed a small data set (n = 20) of CRF38_BF-like pol sequences to estimate the age and demographic history of the CRF38_BF epidemic spreading in Uruguay. This data represent an excellent opportunity to explore potential CRF-specific and regional-specific differences in the patterns of HIV-1 epidemic growth in South America.
HIV-1 CRF12_BF and CRF38_BF data sets.
Characterization of "CRF-like" recombinant profiles
Two strategies were used to characterize the HIV-1 pol sequences used in the present study as CRF12_BF-like or CRF38_BF-like recombinants:
1) First, the recombination breakpoints of each sequence were identified by Bootscanning using Simplot version 3.5.1 . Bootstrap values supporting branching with reference sequences were determined in Neighbor-Joining (NJ) trees constructed using the K2-P  nucleotide substitution model, based on 100 re-samplings, with a 200 bp sliding window moving in steps of 20 bases. Individual query sequences were compared to representative reference sequences of HIV-1 subtypes A1, B, C, and F1. Sequences were considered to have a "CRF-like" profile if recombination sites exactly match those identified in CRF12_BF and CRF38_BF reference sequences.
2) Second, Bayesian and Maximum Likelihood (ML) phylogenetic trees for the final pol alignment including all CRF-like sequences were built to confirm the overall topology and strong support of each CRF clade. Phylogenetic trees were constructed under the GTR  nucleotide substitution model, with a gamma-distribution model of among site rate heterogeneity and a proportion of invariable sites (GTR+I+Γ) selected using the Modeltest program . A Bayesian phylogeny was estimated using MrBayes . Two runs of four chains each were run for 50 × 106 generations, with a burn-in of 5 × 106 generations. Convergence of parameters was assessed by calculating the Effective Sample Size (ESS) using TRACER v1.4 , after excluding an initial 10% for each run. All parameter estimates for each run showed ESS values >100. ML trees were reconstructed with PhyML  using an online web server . Heuristic tree searches were performed using the SPR branch-swapping algorithm, and the approximate likelihood-ratio test (aLRT) based on a Shimodaira-Hasegawa-like procedure was used as a statistical test to calculate branch support. Trees were visualized with the FigTree v1.1.2 program (available at http://tree.bio.ed.ac.uk/software/figtree/).
Estimation of evolutionary rates, dates, and demographic history
The evolutionary rate (μ, units are nucleotide substitutions per site per year), the age of the most recent common ancestor (Tmrca, years), and the mode and rate (r, years-1) of population growth for the CRF12_BF and CRF38_BF strains were estimated using BEAST v1.4.7 [36, 37]. Evolutionary and demographic parameters of CRF12_BF were estimated under a chronological time-scale employing the dates of sample collection. The low number of CRF38_BF sequences analyzed (i.e., 20 sequences), however, was not sufficient to obtain an accurate estimate of the evolutionary rate of this lineage. Therefore, the rates of evolution at pol (PR/RT) gene previously estimated for other HIV-1 group M subtypes (1.5 × 10-3 - 2.5 × 10-3 substitutions/site/year) [38–41] were incorporated as a prior probability distribution in the analysis of this CRF. Estimations of evolutionary and demographic parameters involved two steps. First, the Bayesian skyline plot method  was used to estimate μ, the Tmrca, and the change in effective population size through time. Second, two different demographic models for each data set were compared: exponential and logistic growth; and estimates of the population growth rate were then obtained under the model that provided the best fit to the demographic signal in each data set. Model comparisons in a Bayesian framework were performed by calculating the Bayes Factor (BF)  with TRACER v1.4. Analyses were performed using the GTR+I+Γ nucleotide substitution model under either strict or uncorrelated Lognormal relaxed  molecular clock models. Two separate MCMC chains were run for 10-50 × 106 generations, with a burn-in of 1-5 × 106. BEAST output was analysed using TRACER v1.4, with uncertainty in parameter estimates reflected in the 95% Highest Probability Density (HPD) intervals. Convergence of parameters was assessed through the ESS, with all parameter estimates for each run showing ESS values >100. A graphical representation of the effective number of infections through time was generated by using programs TRACER v1.4 and Prism 4 (GraphPad Software). Posterior trees samples from BEAST runs were summarized using TreeAnnotator v1.4.7 (available from http://beast.bio.ed.ac.uk) to generate time-scaled maximum clade credibility trees.
A total of 115 pol sequences (Argentina = 91, Uruguay = 24) with a "CRF12_BF-like" recombination profile, and 17 pol sequences (Uruguay) with a "CRF38_BF-like" recombinant pattern were identified by Bootscanning analyses. These sequences were aligned with reference sequences of CRF12_BF, CRF38_BF and Brazilian CRFs_BF1, and analyzed using Bayesian and ML approaches. Both phylogenetic approaches showed that the CRF12_BF-like and CRF38_BF-like pol sequences segregated with their respective CRF reference sequences in two well supported monophyletic groups characterized by unique recombination profiles, confirming the common ancestry of each CRF (Fig. 1). Of note, Simplot analysis suggests that the CRF38_BF presents a more complex BF1 mosaic pattern at the PR/RT genomic region than that previously described , characterized by the presence of small subtype F1 fragments between positions 2640 and 3020 relative to HXB2 (Fig. 1). More detailed analysis of the pol genomic region should be performed in order to determine the precise mosaic structure of the CRF38_BF at that region.
Within the CRF12_BF clade, two strongly supported Uruguayan subclusters comprising four (cluster UY-1) and six (cluster UY-2) viruses were identified in both Bayesian (posterior probability [PP] > 0.80) (Fig. 1) and ML (aLTR > 0.70) phylogenetic trees (data not shown). Most (61%) CRF12_BF Uruguayan sequences, however, were randomly interspersed among Argentine sequences, which provides evidence against the existence of a specific Argentine or Uruguayan CRF12_BF lineage. This contrasts with the circulation of CRF38_BF which seems to be restricted to Uruguay, as no strains with a CRF38_BF-like structure were identified after analysis of more than 300 BF1 recombinant pol sequences from Argentina (74 unpublished sequences and 249 sequences retrieved from the Los Alamos HIV Sequence Database).
Bayesian estimates of evolutionary parameters of the HIV-1 CRF12_BF and CRF38_BF epidemics.
2.4 × 10-3
(1.9 × 10-3-2.9 × 10-3)
2.4 × 10-3
(1.8 × 10-3-3.1 × 10-3)
2.4 × 10-3
(1.9 × 10-3-2.8 × 10-3)
2.5 × 10-3
(1.9 × 10-3-3.0 × 10-3)
1.8 × 10-3
(1.5 × 10-3-2.2 × 10-3)
1.9 × 10-3
(1.5 × 10-3-2.3 × 10-3)
1.8 × 10-3
(1.5 × 10-3-2.1 × 10-3)
1.8 × 10-3
(1.4 × 10-3-2.3 × 10-3)
Bayes Factors (BF) between exponential (Exp) and logistic (Log) growth demographic models for the HIV-1 CRF12_BF and CRF38_BF pol data sets.
Evidence against H0b
Exp (H0) vs Log (H1) Strict clock
Exp (H0) vs Log (H1) Relaxed clock
Log Strict (H0) vs Relaxed (H1) clock
Exp (H0) vs Log (H1) Strict clock
Exp (H0) vs Log (H1) Relaxed clock
Log Strict (H0) vs Relaxed (H1) clock
Bayesian estimates of population dynamic parameters of the HIV-1 CRF12_BF and CRF38_BF epidemics.
We have performed an extensive study of the evolutionary history and population dynamics of the two most prevalent HIV-1 BF1 recombinant strains (CRF12_BF and CRF38_BF) spreading in South America. The CRF12_BF was the first CRF to be described in South America among samples isolated from Argentina and Uruguay , and is responsible for a significant part of the HIV-epidemics in those countries. Our phylogenetic analyses of recombinant pol sequences with a CRF12-like structure showed for the first time that CRF12_BF viruses spreading in Argentina and Uruguay constitute a single epidemic with evidences of multiple genetic exchanges among countries. In contrast, circulation of the recently described CRF38_BF  seems to be restricted to Uruguay as no strains with a CRF38_BF-like pol mosaic pattern were observed after screening of a large number (n > 300) of BF1 recombinant sequences from Argentina.
According to our estimates, the mean Tmrca of CRF12_BF and CRF38_BF was 1982-1983 and 1985-1986, respectively; only a few years later than the mean estimated onset date of subtype F1 spread in Brazil (1976-1978) [23, 24]. These mean estimates are fully consistent with epidemiological data that revealed that BF1 recombinants related to the CRF12_BF have been in circulation in Argentine children since the mid 1980s . In agreement with this study, we identified five Argentine children born between 1987 and 1989 infected by CRF12-like BF viruses, and one Uruguayan child born in 1987 infected by a CRF38-like BF strain.
Our study showed that the substitution rate of CRF12_BF viruses at pol gene (2.5 × 10-3 subst./site/year) is similar to rates previously described for HIV-1 subtypes B, C, and CRF31_BC [38, 40, 41], pointing to no major differences in evolution rate among HIV-1 strains circulating in South America. By contrast, some differences in epidemic growth rate may have existed across distinct HIV-1 variants circulating in South America. The estimated initial growth rate of the CRF12_BF epidemic in Argentina and Uruguay was higher than that described for Brazilian subtype B, C, and F1 epidemics [23, 47], and similar to that recently described for the CRF31_BC viruses circulating in southern Brazil  (Table 4); supporting the notion of a fast initial spread of CRF12_BF-like viruses . The initial expansion rate of CRF38_BF also seems to be higher than those described for subtypes B and F1, although the 95% HPD interval of such an estimate was quite large (Table 4).
The rapid initial growth rate of CRF12_BF, CRF31_BC, and CRF38_BF epidemics in South America may indicate that these CRFs displayed a higher fitness and/or transmissibility than parental HIV-1 subtypes. Alternatively, variations in the initial growth rate of HIV-1 variants may reflect differences in the susceptible populations that characterized the initial spread of each strain. Of note, the BF recombinants have been associated with injection drug use populations in Argentina and Uruguay [3, 5], and introduction of a BF1 recombinant virus in such highly connected networks may explain the emergence and rapid initial dissemination of the CRF_BF1 viruses.
After the initial period of fast exponential growth, the expansion rate of the CRF12_BF and CRF38_BF epidemics slowed down since the mid 1990s. The same demographic pattern was described for HIV-1 epidemics in Brazil, USA, and some European countries [23, 38, 47, 49–51]. Such a recent decline in the growth rate of these HIV-1 epidemics may be the consequence of adequate prevention campaigns implemented after the official report of the first AIDS cases in the early 1980s, and/or the result of a saturation of high-risk transmission networks that are loosely connected with low-risk subgroups that exhibit modest levels of HIV-1 infection .
The overall demographic pattern of the CRF12_BF epidemic in Argentina and Uruguay observed in this study is similar with that previously described for the CRF12_BF epidemic in Argentine children ; and our mean estimates of Tmrca and population growth rate fell within the 95% HPD interval of the previous estimates (Table 4). The confidence intervals of these new estimates, however, were much narrower than those previously obtained (Table 4), indicating that the use of a larger data set have substantially improved the accuracy of parameter estimation. In agreement with this idea, the small CRF38 data set was associated with much larger confidence intervals of demographic parameter estimates than those obtained for the larger CRF12_BF data set. Thus, more precise estimations of demographic parameters of CRF38_BF will require the analysis of a larger number of sequences.
Our results suggest that CRF12_BF and CRF38_BF viruses were rapidly generated after the introduction of subtype F1 into South America and have been circulating in the continent over the last 25 years. Both CRFs seem to have spread exponentially at a fast rate during the 1980s and the early 1990s; but the rate of spread of the CRF12_BF and CRF38_BF viruses slowed down since the mid 1990s. Despite similar emergence dates and demographic histories, CRF12_BF have been widely disseminated in Argentina and Uruguay, whereas the CRF38_BF circulation was found limited to Uruguay. Determination of the factors that have shaped the pattern and rate of spread of distinct HIV-1 variants in South America is of paramount importance to understand the epidemic potential of these variants.
We wish to thank Dr Edward C. Holmes and Dr Vera Bongertz for critical reviewing of the manuscript. This study was supported by the Brazilian CNPq/PROSUL grant 490585/2007-8. Work in the "Laboratorio de Biología Celular y Retrovirus" was supported by the "Agencia Nacional de Promoción Científica y Tecnológica" Award PICT 2004-25830, and CONICET Award PIP 6057. Work in the "Centro Nacional de Microbiología" was supported by grants SAF 2005/03833 and SAF 2007/61036 of the "Plan Nacional del SIDA", and by the "Red Tematica Cooperativa de Investigación en SIDA (Red de grupos 173)" of the "Fondo de Investigaciones Sanitarias de la Seguridad Social (FISss)". GB was funded by a fellowship from the Brazilian CNPq/FIOCRUZ "Pesquisador Visitante" Program.
- Thomson MM, Villahermosa ML, Vazquez-de-Parga E, Cuevas MT, Delgado E, Manjon N, Medrano L, Perez-Alvarez L, Contreras G, Carrillo MG, Salomón H, Nájera R: Widespread circulation of a B/F intersubtype recombinant form among HIV-1-infected individuals in Buenos Aires, Argentina. AIDS. 2000, 14: 897-899. 10.1097/00002030-200005050-00020.View ArticlePubMed
- Guimaraes ML, dos Santos Moreira A, Loureiro R, Galvao-Castro B, Morgado MG: High frequency of recombinant genomes in HIV type 1 samples from Brazilian southeastern and southern regions. AIDS Res Hum Retroviruses. 2002, 18: 1261-1269. 10.1089/088922202320886307.View ArticlePubMed
- Hierholzer J, Montano S, Hoelscher M, Negrete M, Hierholzer M, Avila MM, Carrillo MG, Russi JC, Vinoles J, Alava A, Acosta ME, Gianella A, Andrade R, Sanchez JL, Carrion G, Sanchez JL, Russell K, Robb M, Birx D, McCutchan F, Carr JK: Molecular Epidemiology of HIV Type 1 in Ecuador, Peru, Bolivia, Uruguay, and Argentina. AIDS Res Hum Retroviruses. 2002, 18: 1339-1350. 10.1089/088922202320935410.View ArticlePubMed
- Teixeira SL, Bastos FI, Telles PR, Hacker MA, Brigido LF, de FOCA, Bongertz V, Morgado MG: HIV-1 infection among injection and ex-injection drug users from Rio de Janeiro, Brazil: prevalence, estimated incidence and genetic diversity. J Clin Virol. 2004, 31: 221-226. 10.1016/j.jcv.2004.03.016.View ArticlePubMed
- Espinosa A, Vignoles M, Carrillo MG, Sheppard H, Donovan R, Peralta LM, Rossi D, Radulich G, Salomon H, Weissenbacher M: Intersubtype BF recombinants of HIV-1 in a population of injecting drug users in Argentina. J Acquir Immune Defic Syndr. 2004, 36: 630-636. 10.1097/00126334-200405010-00012.View ArticlePubMed
- Quarleri JF, Rubio A, Carobene M, Turk G, Vignoles M, Harrigan RP, Montaner JS, Salomon H, Gomez-Carrillo M: HIV type 1 BF recombinant strains exhibit different pol gene mosaic patterns: descriptive analysis from 284 patients under treatment failure. AIDS Res Hum Retroviruses. 2004, 20: 1100-1107. 10.1089/aid.2004.20.1100.View ArticlePubMed
- Brigido LF, Franco HM, Custodio RM, Oliveira CA, JL PF, Eira M, Bergel F, Araujo F, Carvalheiro JR, Rodrigues R: Molecular characteristics of HIV type 1 circulating in Sao Paulo, Brazil. AIDS Res Hum Retroviruses. 2005, 21: 673-682. 10.1089/aid.2005.21.673.View ArticlePubMed
- Gomez-Carrillo M, Pampuro S, Duran A, Losso M, Harris DR, Read JS, Duarte G, De Souza R, Soto-Ramirez L, Salomon H: Analysis of HIV type 1 diversity in pregnant women from four Latin American and Caribbean countries. AIDS Res Hum Retroviruses. 2006, 22: 1186-1191. 10.1089/aid.2006.22.1186.View ArticlePubMed
- Dilernia DA, Gomez AM, Lourtau L, Marone R, Losso MH, Salomon H, Gomez-Carrillo M: HIV type 1 genetic diversity surveillance among newly diagnosed individuals from 2003 to 2005 in Buenos Aires, Argentina. AIDS Res Hum Retroviruses. 2007, 23: 1201-1207. 10.1089/aid.2007.0068.View ArticlePubMed
- Brennan CA, Brites C, Bodelle P, Golden A, Hackett J, Holzmayer V, Swanson P, Vallari A, Yamaguchi J, Devare S, Pedroso C, Ramos A, Badaro R: HIV-1 strains identified in Brazilian blood donors: significant prevalence of B/F1 recombinants. AIDS Res Hum Retroviruses. 2007, 23: 1434-1441. 10.1089/aid.2007.0121.View ArticlePubMed
- Rios M, Delgado E, Perez-Alvarez L, Fernandez J, Galvez P, de Parga EV, Yung V, Thomson MM, Najera R: Antiretroviral drug resistance and phylogenetic diversity of HIV-1 in Chile. J Med Virol. 2007, 79: 647-656. 10.1002/jmv.20881.View ArticlePubMed
- de Sa-Filho DJ, Soares Mda S, Candido V, Gagliani LH, Cavaliere E, Diaz RS, Caseiro MM: HIV type 1 pol gene diversity and antiretroviral drug resistance mutations in Santos, Brazil. AIDS Res Hum Retroviruses. 2008, 24: 347-353. 10.1089/aid.2007.0203.View ArticlePubMed
- Aguayo N, Laguna-Torres VA, Villafane M, Barboza A, Sosa L, Chauca G, Carrion G, Coenca B, Perez J, Galeano A, Bautista CT, Sanchez JL, Carr JK, Kochel T: Epidemiological and molecular characteristics of HIV-1 infection among female commercial sex workers, men who have sex with men and people living with AIDS in Paraguay. Rev Soc Bras Med Trop. 2008, 41: 225-231. 10.1590/S0037-86822008000300001.View ArticlePubMed
- Ruchansky D, Casado C, Russi JC, Arbiza JR, Lopez-Galindez C: Identification of a new HIV Type 1 circulating recombinant form (CRF38_BF1) in Uruguay. AIDS Res Hum Retroviruses. 2009, 25: 351-356. 10.1089/aid.2008.0248.View ArticlePubMed
- De Sa Filho DJ, Sucupira MC, Casiero MM, Sabino EC, Diaz RS, Janini LM: Identification of two HIV type 1 circulating recombinant forms in Brazil. AIDS Res Hum Retroviruses. 2006, 22: 1-13. 10.1089/aid.2006.22.1.View ArticlePubMed
- Guimaraes ML, Eyer-Silva WA, Couto-Fernandez JC, Morgado MG: Identification of two new CRF_BF in Rio de Janeiro State, Brazil. AIDS. 2008, 22: 433-435. 10.1097/QAD.0b013e3282f47ad0.View ArticlePubMed
- Thomson MM, Sierra M, Tanuri A, May S, Casado G, Manjon N, Najera R: Analysis of near full-length genome sequences of HIV type 1 BF intersubtype recombinant viruses from Brazil reveals their independent origins and their lack of relationship to CRF12_BF. AIDS Res Hum Retroviruses. 2004, 20: 1126-1133. 10.1089/aid.2004.20.1126.View ArticlePubMed
- Sa Filho DJ, Sanabani S, Diaz RS, Munerato P, Brunstein A, Fusuma E, Sabino EC, Janini LM: Analysis of full-length human immunodeficiency virus type 1 genome reveals a variable spectrum of subtypes B and f recombinants in Sao Paulo, Brazil. AIDS Res Hum Retroviruses. 2005, 21: 145-151. 10.1089/aid.2005.21.145.View ArticlePubMed
- Sanabani S, Neto WK, Kalmar EM, Diaz RS, Janini LM, Sabino EC: Analysis of the near full length genomes of HIV-1 subtypes B, F and BF recombinant from a cohort of 14 patients in Sao Paulo, Brazil. Infect Genet Evol. 2006, 6 (5): 368-77. 10.1016/j.meegid.2006.01.003.View ArticlePubMed
- Carr JK, Avila M, Gomez Carrillo M, Salomon H, Hierholzer J, Watanaveeradej V, Pando MA, Negrete M, Russell KL, Sanchez J, Birx DL, Andrade R, Vinoles J, McCutchan FE: Diverse BF recombinants have spread widely since the introduction of HIV-1 into South America. AIDS. 2001, 15: F41-47. 10.1097/00002030-200110190-00002.View ArticlePubMed
- Thomson MM, Delgado E, Herrero I, Villahermosa ML, Vazquez-de Parga E, Cuevas MT, Carmona R, Medrano L, Perez-Alvarez L, Cuevas L, Najera R: Diversity of mosaic structures and common ancestry of human immunodeficiency virus type 1 BF intersubtype recombinant viruses from Argentina revealed by analysis of near full-length genome sequences. J Gen Virol. 2002, 83: 107-119.View ArticlePubMed
- Bello G, Guimaraes ML, Morgado MG: Evolutionary history of HIV-1 subtype B and F infections in Brazil. AIDS. 2006, 20: 763-768. 10.1097/01.aids.0000216377.84313.52.View ArticlePubMed
- Bello G, Eyer-Silva WA, Couto-Fernandez JC, Guimaraes ML, Chequer-Fernandez SL, Teixeira SL, Morgado MG: Demographic history of HIV-1 subtypes B and F in Brazil. Infect Genet Evol. 2007, 7: 263-270. 10.1016/j.meegid.2006.11.002.View ArticlePubMed
- Aulicino PC, Bello G, Rocco C, Romero H, Mangano A, Morgado MG, Sen L: Description of the First Full-Length HIV Type 1 Subtype F1 Strain in Argentina: Implications for the Origin and Dispersion of This Subtype in South America. AIDS Res Hum Retroviruses. 2007, 23: 1176-1182. 10.1089/aid.2007.0038.View ArticlePubMed
- Sierra M, Thomson MM, Rios M, Casado G, Castro RO, Delgado E, Echevarria G, Munoz M, Colomina J, Carmona R, Vega Y, Parga EV, Medrano L, Pérez-Alvarez L, Contreras G, Nájera R: The analysis of near full-length genome sequences of human immunodeficiency virus type 1 BF intersubtype recombinant viruses from Chile, Venezuela and Spain reveals their relationship to diverse lineages of recombinant viruses related to CRF12_BF. Infect Genet Evol. 2005, 5: 209-217. 10.1016/j.meegid.2004.07.010.View ArticlePubMed
- Aulicino PC, Holmes EC, Rocco C, Mangano A, Sen L: Extremely rapid spread of human immunodeficiency virus type 1 BF recombinants in Argentina. J Virol. 2007, 81: 427-429. 10.1128/JVI.01403-06.PubMed CentralView ArticlePubMed
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMed
- Ray S: Simplot v2.5.0. [http://sray.med.som.jhmi.edu/SCRoftware/simplot/]
- Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.View ArticlePubMed
- Tavaré S: Some probabilistic and statistical problems in the analysis of DNA sequences. p. Some mathematical questions in biology--DNA sequence analysis Providence (RI): American Mathematical Society. Edited by: Miura RM. 1986, 57-86.
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMed
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMed
- Rambaut A, Drummond A: Tracer v1.4. 2007, [http://beast.bio.ed.ac.uk/Tracer]
- 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 ArticlePubMed
- Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online--a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33: W557-559. 10.1093/nar/gki352.PubMed CentralView ArticlePubMed
- Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W: Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics. 2002, 161: 1307-1320.PubMed CentralPubMed
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMed
- Hue S, Pillay D, Clewley JP, Pybus OG: Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups. Proc Natl Acad Sci USA. 2005, 102: 4425-4429. 10.1073/pnas.0407534102.PubMed CentralView ArticlePubMed
- Salemi M, de Oliveira T, Ciccozzi M, Rezza G, Goodenow MM: High-Resolution Molecular Epidemiology and Evolutionary History of HIV-1 Subtypes in Albania. PLoS ONE. 2008, 3: e1390-10.1371/journal.pone.0001390.PubMed CentralView ArticlePubMed
- Bello G, Passaes CP, Guimaraes ML, Lorete RS, Matos Almeida SE, Medeiros RM, Alencastro PR, Morgado MG: Origin and evolutionary history of HIV-1 subtype C in Brazil. AIDS. 2008, 22: 1993-2000. 10.1097/QAD.0b013e328315e0aa.View ArticlePubMed
- Passaes CP, Bello G, Lorete RS, Matos Almeida SE, Junqueira DM, Veloso VG, Morgado MG, Guimaraes ML: Genetic characterization of HIV-1 BC recombinants and evolutionary history of the CRF31_BC in Southern Brazil. Infect Genet Evol. 2009, 9: 474-482. 10.1016/j.meegid.2009.01.008.View ArticlePubMed
- 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 ArticlePubMed
- Suchard MA, Weiss RE, Sinsheimer JS: Bayesian selection of continuous-time Markov chain evolutionary models. Mol Biol Evol. 2001, 18: 1001-1013.View ArticlePubMed
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.PubMed CentralView ArticlePubMed
- Worobey M, Gemmel M, Teuwen DE, Haselkorn T, Kunstman K, Bunce M, Muyembe JJ, Kabongo JM, Kalengayi RM, Van Marck E, Gilbert MT, Wolinsky SM: Direct evidence of extensive diversity of HIV-1 in Kinshasa by 1960. Nature. 2008, 455: 661-664. 10.1038/nature07390.PubMed CentralView ArticlePubMed
- Guimaraes ML, Vicente AC, Otsuki K, da Silva RF, Francisco M, da Silva FG, Serrano D, Morgado MG, Bello G: Close phylogenetic relationship between Angolan and Romanian HIV-1 subtype F1 isolates. Retrovirology. 2009, 6: 39-10.1186/1742-4690-6-39.PubMed CentralView ArticlePubMed
- Bello G, Guimaraes ML, Passaes CP, Almeida SE, Veloso VG, Morgado MG: Short Communication: Evidences of Recent Decline in the Expansion Rate of the HIV Type 1 Subtype C and CRF31_BC Epidemics in Southern Brazil. AIDS Res Hum Retroviruses. 2009, 25: 1065-1069. 10.1089/aid.2009.0106.View ArticlePubMed
- Gomez Carrillo M, Avila M, Hierholzer J, Pando M, Martinez PL, McCutchan FE, Carr JK: Mother-to-child HIV type 1 transmission in Argentina: BF recombinants have predominated in infected children since the mid-1980s. AIDS Res Hum Retroviruses. 2002, 18: 477-483. 10.1089/088922202317406619.View ArticlePubMed
- Robbins KE, Lemey P, Pybus OG, Jaffe HW, Youngpairoj AS, Brown TM, Salemi M, Vandamme AM, Kalish ML: U.S. Human immunodeficiency virus type 1 epidemic: date of origin, population history, and characterization of early strains. J Virol. 2003, 77: 6359-6366. 10.1128/JVI.77.11.6359-6366.2003.PubMed CentralView ArticlePubMed
- Walker PR, Pybus OG, Rambaut A, Holmes EC: Comparative population dynamics of HIV-1 subtypes B and C: subtype-specific differences in patterns of epidemic growth. Infect Genet Evol. 2005, 5: 199-208. 10.1016/j.meegid.2004.06.011.View ArticlePubMed
- Paraskevis D, Magiorkinis E, Magiorkinis G, Sypsa V, Paparizos V, Lazanas M, Gargalianos P, Antoniadou A, Panos G, Chrysos G, Sambatakou H, Karafoulidou A, Skoutelis A, Kordossis T, Koratzanis G, Theodoridou M, Daikos GL, Nikolopoulos G, Pybus OG, Hatzakis A, Multicentre Study on HIV Heterogeneity: Increasing prevalence of HIV-1 subtype A in Greece: estimating epidemic history and origin. J Infect Dis. 2007, 196: 1167-1176. 10.1086/521677.View ArticlePubMed
- Blower S: Behaviour change and stabilization of seroprevalence levels in communities of injecting drug users: correlation or causation?. J Acquir Immune Defic Syndr. 1991, 4: 920-923.PubMed
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.