Research | Open | Published:
The HIV-1 epidemic in Bolivia is dominated by subtype B and CRF12_BF "family" strains
Virology Journalvolume 9, Article number: 19 (2012)
Molecular epidemiological studies of HIV-1 in South America have revealed the occurrence of subtypes B, F1 and BF1 recombinants. Even so, little information concerning the HIV-1 molecular epidemiology in Bolivia is available. In this study we performed phylogenetic analyses from samples collected in Bolivia at two different points in time over a 10 year span. We analyzed these samples to estimate the trends in the HIV subtype and recombinant forms over time.
Materials and methods
Fifty one HIV-1 positive samples were collected in Bolivia over two distinct periods (1996 and 2005). These samples were genetically characterized based on partial pol protease/reverse transcriptase (pr/rt) and env regions. Alignment and neighbor-joining (NJ) phylogenetic analyses were established from partial env (n = 37) and all pol sequences using Mega 4. The remaining 14 env sequences from 1996 were previously characterized based on HMA-env (Heteroduplex mobility assay). The Simplot v.3.5.1 program was used to verify intragenic recombination, and SplitsTree 4.0 was employed to confirm the phylogenetic relationship of the BF1 recombinant samples.
Phylogenetic analysis of both env and pol regions confirmed the predominance of "pure" subtype B (72.5%) samples circulating in Bolivia and revealed a high prevalence of BF1 genotypes (27.5%). Eleven out of 14 BF1 recombinants displayed a mosaic structure identical or similar to that described for the CRF12_BF variant, one sample was classified as CRF17_BF, and two others were F1pol/Benv. No "pure" HIV-1 subtype F1 or B" variant of subtype B was detected in the present study. Of note, samples characterized as CRF12_BF-related were depicted only in 2005.
HIV-1 genetic diversity in Bolivia is mostly driven by subtype B followed by BF1 recombinant strains from the CRF12_BF "family". No significant temporal changes were detected between the mid-1990s and the mid-2000s for subtype B (76.2% vs 70.0%) or BF1 recombinant (23.8% vs 30.0%) samples from Bolivia.
Latin America contains the third highest estimated number of people living with HIV in the world . The HIV epidemic in countries from the Southern Cone is characterized by the co-circulation of HIV-1 subtypes B, F1, and BF1 recombinant forms [2–4]. Despite such overall similarity, the mosaic patterns of the HIV-1 BF1 Circulanting Recombinant Forms (CRFs_BF) and Unique Recombinant Forms (URFs_BF) display some important differences among countries.
The HIV-1 BF1 epidemics in Argentina [5–9], Chile , Paraguay  and Uruguay  are dominated by the CRF12_BF and URFs_BF sharing CRF12-related structures, in other words sharing some of CRF12_BF breakpoints. In addition, new local CRFs_BF with a CRF12-related structure have been identified within each country, including CRF17_BF in Argentina and Paraguay [9, 11], CRF38_BF in Uruguay , and CRF44_BF in Chile . By contrast, the HIV-1 BF1 Brazilian epidemic is driven mainly by a myriad of different URFs_BF and at least five distinct CRFs_BF [14–19], unrelated to CRF12_BF, which usually show no common mosaic recombinant profile among one another. Another important difference is that some "pure" full-length subtype F1 genomes have been described in Brazil , while the frequency of such full-length subtype F1 viruses or even non-recombinant subtype F1 pol sequences is almost absent in other South American countries [8, 12, 21, 22].
Bolivia is a landlocked country in central South America, bordered by Brazil to the north and east, Paraguay and Argentina to the south, and Chile and Peru to the west. The total number of HIV/AIDS cases reported in Bolivia was 5,184 from 1984 through March 2010 , with an estimated country-wide underreporting rate of over 70% . UNAIDS estimates that approximately 12,000 people in Bolivia were living with HIV in 2009, and of those cases, 500 to 1500 were new infections . These estimates demonstrate that the Bolivian epidemic is still potentially growing. Very scarce data on HIV epidemiology in Bolivia are available, with only two major previous studies conducted in the 1990s, which revealed the occurrence of subtypes B, F1 and BF1 recombinant forms [22, 25], similar to what has been described for other countries from the Southern Cone [2–13]. In these previous studies, however, most samples were genetically characterized on the basis of HMA env analysis, and very few (n = 14) HIV-1 nucleotide sequences from Bolivian seropositive patients are available at the Los Alamos HIV database, from those only three were full length genomes that reveal the circulation of CRF12_BF  and CRF17_BF.
The objective of the present study was to assess the trends in the HIV molecular epidemiology in Bolivia based on the phylogenetic analysis of HIV-1 pol and env sequences at two different time points (1996 and 2005).
In order to assess the trends of the HIV-1 molecular epidemiology in Bolivia, 62 samples were collected at two points in time, 1996 (30 samples) and 2005 (32 samples), and submitted to amplification and sequencing of both the pol (pr/rt) and env genomic regions. Eleven samples (nine from 1996 and two from 2005) failed to be amplified or sequenced from the proviral DNA in one or both genomic regions and were excluded from the study. The remaining 51 samples showed a predominance of the B subtype (72.5%) and BF1 recombinants (27.5%) (Figure 1). Of note, no pure subtype F1 genomes (F1pol/F1env) were detected. Fourteen highly supported (> 95%) monophyletic clusters of at least two sequences could be detected in the pol (Figure 1a) and 11 in env region (Figure 2b). Some clusters were composed of sequences retrieved from individuals with a known direct epidemiological link, such as sexual partners (clusters 2a, 5 and 6a), and a mother/child pair (cluster 7). No apparent direct epidemiological relationships were detected for the remaining 10 clusters (1, 2, 3, 4, 6, 8, 9, 10, 11, 12). Among the 12 monophyletic clusters, nine displayed a Bpol/Benv profile (1-9), two were BF1pol/F1env (10, 11), and one was F1pol/Benv (12) (Figure 1).
The 14 BF1 recombinant samples identified in this study displayed the following recombination profiles: BF1pol/F1env (n = 11), BF1pol/Benv (n = 1) and F1pol/Benv (n = 2) (Figure 1). A more detailed evaluation of the 12 BF1 pol recombinant samples was performed using bootscan and SplitsTree analyses. From these, five samples displayed a recombination profile in pr/rt identical to the CRF12_BF, being classified as CRF12_BF-like (96BOLP29, 96BOSC14, 05BO06, 05BO15, 05BO18), and one sample exhibited a CRF17_BF-like mosaic pattern (96BOLP23) (Figure 2).
No epidemiological data were available to confirm the proposed close genetic relationship between the CRF17_BF (BO02.BO119) strain available in the Los Alamos Sequence Database and the BF recombinant genome described in the present study. The remaining six BF1 pol samples segregate in two highly supported monophyletic clusters and exhibited two different CRF12_BF-related (F1/B/F1/B) mosaic profiles (Figure 2). In the first cluster (cluster 10), which is composed of samples 05BO09 and 05BO20, the first (F1/B) and second (B/F1) recombination breakpoints at pol coincide with those observed in the CRF12_BF, while the third one (F1/B) is displaced downstream. In the second cluster (cluster 11), which is composed of samples 05BO03, 05BO16, 05BO22 and 05BO24, only the third (F1/B) recombination breakpoint at pol coincides with those observed in the CRF12_BF, while the first (F1/B) and the second (B/F1) recombination breakpoints at pol were displaced upstream in comparison to those detected in CRF12_BF (Figure 3). None of the CRFs_BF pol mosaic patterns described in Brazil (CRF28, CRF29, CRF39, and CRF40), Uruguay (CRF38), or Chile (CRF44) was observed for the HIV-1 samples obtained from Bolivia.
Table 1summarizes the subtype distribution and recombination profiles of the 51 HIV-1 samples from Bolivia according to the time period. No significant temporal change in the HIV-1 subtype distribution was observed in Bolivia between 1996 and 2005. During this period the prevalence of subtype B decreased from 76.2% to 70%, while the prevalence of BF1 recombinants increased from 23.8% to 30%.
The present study represents the most comprehensive analysis of the HIV-1 molecular diversity in Bolivia performed to date. From a total of 51 samples investigated in both pol (pr/rt) and env (C2-V3) regions, 72.5% of cases were classified as "pure" subtype B, and 27.5% were classified as BF1 recombinants. Eleven out of 14 BF1 recombinants displayed a mosaic structure identical or similar to that described for the CRF12_BF strain, and the remaining two BF1 recombinants displayed a F1(pol)/B(env) mosaic structure. Subtype F1 (pol)/F1 (env) strains were not detected in the present study.
The high frequency of BF1 recombinants with CRF12_BF-like or CRF12_BF-related structures and the absence of non-recombinant subtype F1 strains found in Bolivia reveal a scenario closely related to that found in Argentina, Chile, Paraguay and Uruguay [6–12], and different from that found in other neighboring countries, such as Brazil and Peru [17–19, 22]. This pattern suggests a more intense HIV-1 influx into Bolivia from the southern border (Argentina and/or Paraguay) than from the eastern (Peru) or western/northern (Brazil) frontiers. It is interesting to note that none of the subtype B samples studied presented the typical Brazilian motif GWGR in the top of the v3 loop (B" variant of subtype B), which is highly (40%) represented in Brazil . The reason for such potential unequal viral genetic influxes into Bolivia from neighboring countries remains unknown and deserves further investigation.
In addition to the widely disseminated CRF12_BF variant, other CRFs_BF variants with a CRF12_BF-related structure were also described as circulating in Argentina and Paraguay (CRF17_BF), Uruguay (CRF38_BF) and Chile (CRF44_BF) [9–13]. One sequence with a CRF17_BF-like recombinant structure was also detected in our study, confirming the circulation of this CRF in Bolivia since 1996. Of note, we also detected one highly supported cluster (cluster 11) of four BF1 sequences displaying the same CRF12_BF-related (pol)/F1 (env) mosaic profile, obtained from epidemiologically unlinked HIV-1 seropositive patients from Bolivia. Such a result may represent a new local Bolivian CRF_BF variant. Full-length genome analysis of these four BF1 sequences will be necessary to confirm this hypothesis.
It is interesting to note that all BF1 recombinants with a CRF12_BF-related structure detected in Bolivia displayed larger subtype B segments at pol than CRF12_BF. A similar finding was also observed in Argentina [5–9] and Chile . The generation of CRF12_BF-related recombinants with larger subtype F1 fragments in those countries is unexpected because of the scarcity (or complete absence) of pure subtype F1 variants outside of Brazil. It has been proposed that the CRF12_BF strain and other CRFs_BF and URFs_BF variants with related structures circulating in South America, as well as other groups of related recombinants should be considered components of the same "family" of recombinants [3, 27].
Since its emergence around the early 1980s, the CRF12_BF variant has displayed an initial exponential expansion in Argentina and Uruguay . Such a fast rate of spread and the high prevalence (> 50%) of CRF12_BF and CRF12_BF-related mosaic forms in those countries could be interpreted as evidence of the high viral fitness and/or transmissibility of the CRF12 "family" of recombinants. However, much lower prevalence (< 25%) of CRF12_BF-like and CRF12_BF-related variants has been found in Bolivia [this study], Chile  and Paraguay . Our results did not depict a significant increase in the proportion of those mosaic HIV-1 forms in Bolivia between 1996 and 2005, nevertheless all samples from 2005 were collected in only one city: La Paz. This finding suggests that the rate of expansion and the resulting prevalence of CRF12_BF "family" variants in South America are country-specific and might be dependent on non-virological factors. It has been also proposed that the prevalence of CRF12_BF and related BF1 strains is higher in countries with more extensive intravenous drug epidemics, possibly due to a founder effect and the rapid spread of HIV in these groups .
In summary, our results demonstrate that HIV-1 genetic diversity in Bolivia is mostly driven by subtype B followed by BF1 recombinant strains from the CRF12_BF "family" and no significant temporal changes in the subtype distributions were detected between the mid-1990s and the mid-2000s.
Materials and methods
A set of 62 samples were collected in Bolivia in 1996 (30 samples) and 2005 (32 samples). Demographic, laboratory and clinical data of Bolivian subjects collected in 1996 were previously described . Samples from 1996 were collected in three distinct Bolivia cities [La Paz (n = 9), Cochabamba (n = 5) and Santa Cruz (n = 16)]. The second cohort was composed exclusively of patients living in La Paz, and most of them were under antiretroviral therapy. All participants signed the informed consent form, and the study was approved by the Bolivian Ministry of Health.
HIV-1 DNA amplification and sequencing
Genomic DNA was extracted from whole blood by the phenol/chloroform method from samples collected in 1996, as described . For samples collected in 2005, the QIAamp DNA kit was used (Qiagen Inc., CA, U.S.A.) according to the manufacturer's protocol.
Amplification of HIV-1 env and pol regions was performed using an in-house nested PCR protocol method as described elsewhere . Samples collected in 1996 were previously subtyped based on the HMA env-gp120 C2-V3 region and 11 of them were also sequenced and phylogenetically analyzed as described . Partial protease/reverse transcriptase (pr/rt) sequences were generated for the samples from both periods, and new env sequences were obtained only from samples collected in 2005. No HIV-1 seropositive individuals from the first time period were undergoing antiretroviral therapy at the time of blood collection; however, more than a half of the individuals in the second group were being treated.
HIV-1 subtype classification
HIV-1 subtype classification was performed by phylogenetic analysis with the neighbor-joining method using the Tamura-Nei model as implemented in Mega 4.0 . Those samples that present divergent clustering in the phylogenetic analysis were investigated for recombination by bootscan analysis as implemented in the Simplot version 3.5.1 (http://www.med.jhu.edu/deptmed/sray/download/), using representative HIV-1 sequences from group M. Bootstrap values supporting branching with reference sequences were determined in neighbor-joining trees constructed using the Kimura 2-parameter model, based on 100 re-samplings, with a 200 nucleotide sliding window moving in steps of 20 bases. The SplitsTree program version 4  was employed to confirm the phylogenetic relationship of the recombinant samples using the NeighborNet based on the pairwise distance estimated by F84 parameter model.
GenBank accession numbers: JN710748-JN710798 (env-gp120 C2-V3 region), JN710799-JN710830 (pr/rt).
UNAIDS Global Report. [http://www.unaids.org/globalreport/Global_report.htm,accessed in 2011]
Morgado MG, Guimaraes ML, Galvao-Castro B: HIV-1 polymorphism: a challenge for vaccine development - a review. Mem Inst Oswaldo Cruz. 2002, 97: 143-150.
Zhang M, Foley B, Schultz AK, Macke JP, Bulla I, Stanke M, Morgenstern B, Korber B: Leitner: The role of recombination in the emergence of a complex and dynamic HIV epidemic. Retrovirology. 2010, 7: 25-10.1186/1742-4690-7-25.
Hemelaar J, Gouws E, Ghys PD, Osmanov S: Global trends in molecular epidemiology of HIV-1 during 2000-2007. AIDS. 2011, 25: 679-689. 10.1097/QAD.0b013e328342ff93.
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-F47. 10.1097/00002030-200110190-00002.
Thomson MM, Delgado E, Herrero I, Villahermosa ML, Vázquez-de Parga E, Cuevas MT, Carmona R, Medrano L, Pérez-Alvarez L, Cuevas L, Nájera 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.
Quarleri JF, Rubio A, Carobene M, Turk G, Vignoles M, Harrigan RP, Montaner JS, Salomón H, Gómez-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.
Dilernia DA, Gomez AM, Lourtau L, Marone R, Losso MH, Salomón H, Gómez-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 (10): 1201-1207. 10.1089/aid.2007.0068.
Aulicino PC, Bello G, Guimaraes ML, Ruchansky D, Rocco C, Mangano A, Morgado MG, Sen L: Longitudinal analysis of HIV-1 BF1 recombinant strains in vertically infected children from Argentina reveals a decrease in CRF12_BF pol gene mosaic patterns and high diversity of BF unique recombinant forms. Infect Genet Evol. 2011, 11: 349-357. 10.1016/j.meegid.2010.11.008.
Ríos M, Delgado E, Pérez-Alvarez L, Fernández J, Gálvez P, de Parga EV, Yung V, Thomson MM, Nájera R: Antiretroviral drug resistance and phylogenetic diversity of HIV-1 in Chile. J Med Virol. 2007, 79: 647-656. 10.1002/jmv.20881.
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. Revista Da Sociedade Brasileira De Medicina Tropical. 2008, 41: 225-231. 10.1590/S0037-86822008000300001.
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.
Delgado E, Ríos M, Fernández J, Pérez-Alvarez L, Nájera R, Thomson MM: Identification of a new HIV type 1 BF intersubtype circulating recombinant form (CRF44_BF) in Chile. AIDS Res Hum Retroviruses. 2010, 26: 821-826. 10.1089/aid.2010.0006.
De Sá Filho DJ, Sucupira MC, Caseiro 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.
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.
Sanabani SS, Pastena ER, Neto WK, Martinez VP, Sabino EC: Characterization and frequency of a newly identified HIV-1 BF1 intersubtype circulating recombinant form in Sao Paulo, Brazil. Virol J. 2010, 7: 74-10.1186/1743-422X-7-74.
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.
de Souza AC, de Oliveira CM, Rodrigues CL, Silva SA, Levi JE: Short communication: molecular characterization of HIV type 1 BF pol recombinants from Sao Paulo, Brazil. AIDS Res Hum Retroviruses. 2008, 24: 1521-1525. 10.1089/aid.2008.0089.
Guimarães ML, Couto-Fernandez JC, Eyer-Silva Wde A, Teixeira SL, Chequer-Fernandez SL, Morgado MG: Analysis of HIV-1 BF pr/rt recombinant strains from Rio de Janeiro/Brazil reveals multiple unrelated mosaic structures. Infect Genet Evol. 2010, 10: 1094-1100. 10.1016/j.meegid.2010.07.001.
Sanabani S, Kleine Neto W, 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: 368-377. 10.1016/j.meegid.2006.01.003.
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.
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.
USAID HIV/AIDS Health profile. Available from http://www.usaid.gov/our_work/global_health/aids/Countries/lac/bolivia_profile.pdf
Protto JP, Schaaf D, Suarez MF, Darras C: Epidemiological environment and response to the HIV epidemic in Bolivia. Rev Panam Salud Publica. 2008, 23: 288-294. 10.1590/S1020-49892008000400012.
Velarde-Dunois KG, Guimarães ML, La Fuente C, Andrade R, Arévalo R, Pantoja S, Mariscal R, Sandoval R, Iriarte F, Chamón V, Melgar ML, Carvajal R, Morgado MG: Molecular characterization of human immunodeficiency virus type 1-infected individuals from bolivia reveals the presence of two distinct genetic subtypes B and F. AIDS Res Hum Retroviruses. 2000, 16: 1921-1926. 10.1089/08892220050195883.
Morgado MG, Guimarães ML, Neves Júnior I, dos Santos VG, Linhares-de-Carvalho MI, Castello-Branco LR, Bastos FI, Castilho EA, Galvão-Castro B, Bongertz V: Molecular epidemiology of HIV in Brazil: polymorphism of the antigenically distinct HIV-1 B subtype strains. Mem Inst Oswaldo Cruz. 1998, 93 (3): 383-386. 10.1590/S0074-02761998000300020.
Thomson MM, Nájera R: Molecular epidemiology of HIV-1 variants in the global AIDS pandemic: an update. AIDS Rev. 2005, 7 (4): 210-224. Review
Bello G, Aulicino PC, Ruchansky D, Guimarães ML, Lopez-Galindez C, Casado C, Chiparelli H, Rocco C, Mangano A, Sen L, Morgado MG: Phylodynamics of HIV-1 circulating recombinant forms 12_BF and 38_BF in Argentina and Uruguay. Retrovirology. 2010, 7: 22-24. 10.1186/1742-4690-7-22.
Bello G, Passaes CP, Guimarães 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.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23: 254-267.
We wish to thank Dr Gonzalo Bello for critically reviewing the manuscript. We also acknowledge Dr Vera Bongertz for English revision. This study was supported by the Brazilian CNPq/PROSUL grant 490585/2007-8.
The authors declare that they have no competing interests.
M.L.G., K.G.V.D and M.G.M conceived and designed the study. D.S and K.G.V.D were responsible for patients' recruitment and sample collection in Bolivia. K.G.V.D, M.L.G. and M.G.M. performed the characterization of the samples from 1996 in the env region. M.L.G. conducted the characterization of the pol region from both studied periods. M.L.G performed the phylogenetic analyses. M.L.G. and M.G.M wrote the first draft and all authors contributed to the final version of the paper.