Temporal dynamics of HIV-1 circulating subtypes in distinct exposure categories in southern Brazil
© Almeida et al.; licensee BioMed Central Ltd. 2012
Received: 2 March 2012
Accepted: 7 December 2012
Published: 12 December 2012
The HIV-1 epidemic in Brazil is predominantly driven by subtype B. However, in Brazilian Southern region subtype C prevails and a relatively high AIDS incidence rate is observed. The aim of the present study was to assess the temporal dynamics of HIV-1 subtypes circulating in patients from distinct exposure categories in Southern Brazil. For this purpose 166 HIV-1 samples collected at the years of 1998 (group I) and 2005–2008 (group II) were analyzed.
Analysis of group I revealed statistically significant (p < 0.05) associations between MSM and subtype B as well as between IDU and subtype C; while no statistical significant association between HIV-1 subtypes and exposure category was verified for group II. An overall temporal increase in the prevalence of subtype C and BC recombinants was observed in both HET and MSM populations, accompanied by a proportional decrease in the prevalence of the pure subtype B.
The present study shows an association between HIV subtypes and exposure categories at the middle 1990s in Southern Brazil. Our findings suggest that MSM and IDU populations might have played a major role in the introduction and initial dissemination of subtypes B and C, respectively, in Southern Brazil. This study also suggests a trend towards homogenization of HIV-1 strains across distinct exposure categories as a consequence of an overall increase in the prevalence of subtype C and BC recombinants in both HET and MSM populations.
KeywordsHIV-1 Brazil Subtypes Exposure categories Temporal dynamics
Acquired immune deficiency syndrome
Circulating recombinant form
Injection drug users
Human immunodeficiency virus
Men who have sex with men
Unique recombinant form.
The main hallmark of the HIV-1 is an extraordinary evolution rate, which results in high molecular diversity and dynamism of the AIDS epidemic . HIV-1 is classified in four groups and the group M, currently estimated to infect around 33 million people around the world, is subdivided in 9 subtypes (A–D, F–H, J, and K) and 54 circulating recombinant forms (CRFs) [1–3].
Since 1980, Brazil has registered 608,230 cases of AIDS, representing an overall prevalence of 0.6% in adult population . The HIV-1 subtype B is the predominant variant in most of the Brazilian regions followed by subtype F1, subtype C and a large variety of BF1 and BC recombinant forms [4–10]. The distribution of HIV-1 subtypes, however, is not homogeneous across the country and a distinct HIV-1 molecular epidemiologic scenario is observed in the Southern region. Composed by the states of Rio Grande do Sul, Santa Catarina and Paraná, the Brazilian Southern region shows a remarkably high prevalence of subtype C and BC recombinant forms [11–21].
The HIV-1 molecular epidemiology in Porto Alegre, capital of the Rio Grande do Sul state, is characterized by a high prevalence of subtype C (~30-40%), subtype B (~30-45%) and the circulating recombinant form (CRF) 31_BC (~10-25%), and lower prevalence of unique recombinant forms (URFs) and subtype F1 [12, 16–18, 21]. The southern region of Brazil is not only characterized by a distinct subtype profile, but also by a relative high AIDS incidence rate. This is particularly evident for the city of Porto Alegre, showing the highest AIDS incidence rate among all the Brazilian capitals since 1997 . The reported AIDS incidence in Porto Alegre in 2010 (99.8 cases per 100,000 habitants) was more than five times higher than the mean incidence for the whole country (17.9 cases per 100,000 habitants) . In this city, almost 80% of the new AIDS cases registered in 2009 correspond to heterosexual (HET) individuals, approximately 10% to men who have sex with men (MSM) and another 10% to injection drug users (IDU) .
Some studies indicate a tendency towards an increasing proportion of subtype C infections over time in some cities from the Southern region. The estimated prevalence of subtype C increased from 36% before 1997 to 53% in 2008 in Rio Grande (Rio Grande do Sul state) [11, 23], and from 56% to 78% between 2004 and 2009 in Florianopolis (Santa Catarina state) . Further studies also point to a different spreading of the HIV-1 variants among the exposure categories in southern Brazil. Studies conducted in all the three states of Brazilian Southern region described that subtype C is more frequently seen in the HET population, while subtype B is more common in MSM [14, 20, 23]. One study also found an association between subtype C infections and the use of intravenous drugs among males in the state of Parana . So far, no significant associations between HIV-1 clades and exposure categories were detected in the city of Porto Alegre, the main capital of Southern Brazil, even though an association of subtype B with anal sex practices and a tendency of subtype C with females have been reported [11, 17, 21].
Most of the previous studies that analyzed the temporal dynamics of HIV-1 diversity in Southern Brazil were based on the stratification of cross-sectional samples according to time of HIV diagnosis and without taking into account the exposure categories of the individuals. In this study, we described the temporal trends of HIV-1 clades in patients of different exposure categories living in the city of Porto Alegre, through the analysis of two groups of samples collected at 1998 and 2005–2008.
Clinical and demographic data of the HIV-positive patients from Porto Alegre at two distinct time periods - 1998 (group I) and 2005–2008 (group II)
Group I (n=83)
Group II (n=83)
33 ± 10
35 ± 10
CD4 T cell count
First positive serology for HIV (years)
HIV-1 amplification, sequencing and subtyping
DNA samples were extracted from 200 μl of whole blood using a QIAamp DNA kit (Qiagen Inc., CA, U.S.A.), according to the manufacturer’s protocol. Amplification and sequencing of the PR/RT region was performed as described elsewhere . The sequences generated were ~1,160 nt long and covered the protease (PR) and part of the reverse transcriptase (RT) genes (nucleotides 2253–3413 relative to HXB2).
Nucleotide sequences were aligned using the Clustal X program  and three strategies were used to well characterize the HIV-1 sequences as pure subtype, CRF-like or URF using the pol alignment: i) a Neighbor-Joining (NJ) phylogenetic tree was first built under the Tamura-Nei substitution model in 1000 bootstrapped data sets, as implemented in MEGA program ; ii) all sequences were subsequently subjected to bootscanning with the Simplot software version 3.5.1 , using reference sequences representative of HIV-1 subtypes A1, B, C, and F1 available in Los Alamos database, http://www.hiv.lanl.gov/components/sequence/HIV/search/search.html). Bootstrap values supporting branching with reference sequences were determined in NJ trees constructed using the K2-parameter model , based on 100 re-samplings, with a 300 nt sliding window moving in steps of 10 bases; iii) to better characterize the recombination breakpoints suggested in the previous analyses, the putative recombinants were subjected to informative site analyses as described elsewhere . Consensus sequences used in our analyses were generated from a total of 491 subtype B and 164 subtype C pol Brazilian sequences downloaded from the Los Alamos HIV Sequence Database, as implemented in the DAMBE program . Sequences were submitted to GeneBank under the following accession numbers: [JF487830 - JF487913] and [JQ619540 - JQ619621].
Statistical comparisons among subtypes groups and exposure categories were made using Pearson’s χ2-test with adjusted residues and Fisher’s exact test when appropriate. Statistical analysis was performed using the SPSS 16.0 statistical package and the significance level was set at p<0.05.
Group I was composed by HIV-1 positive samples collected in 1998 from 83 individuals, of which 71.0% were men. Based on medical records, patients were classified into three exposure categories: HET (50.6%), MSM (38.6%) and IDU (10.8%). The mean diagnostic year of patients from group I was 1994 (Table 1). Group II, comprised 83 HIV-1 positive individuals sampled from 2005 to 2008, of which 59.3% were men (Table 1). Regarding the exposure categories, group II showed 73.5% of HET individuals and 26.5% of MSM, the mean diagnostic year being 2005. Most patients of both groups (39,2% and 50%, respectively) presented T CD4 cells counts above 400cell/mm3. All patients included in group II were antiretroviral naïve, while 55% of individuals from group I declared to be under antiretroviral therapy.
According to the phylogenetic, bootscanning, and informative site analyses of the PR/RT region, the 166 HIV-1 samples from Porto Alegre here analyzed were classified as follows: subtype B (n = 65), subtype C (n = 37), CRF31_BC (n = 28), URFs_BC (n = 25), subtype F1 (n = 3), URFs_BF (n = 5), and URF_BCF (n = 3) (Additional file 1: Figure S1). Some of the URFs_BC detected in the present study (n = 13), shared one of the recombination breakpoints with the CRF31_BC and were probably created by the recombination of the CRF31_BC with local subtypes C and B. Since CRF31_BC and the URFs_BC with a CRF31_BC-related structure probably share a common evolutionary history, they were grouped together for further analyses. Other URFs_BC with a mosaic structure not related to the CRF31_BC (n = 12) were allocated in a separate group. The recombination pattern of all URFs and a schematic draw of the CRF31_BC-related structure are provided as Additional file (Additional file 2: Figure S2).
HIV-1 subtype frequencies according to the patient’s exposure category
CRF31_BC and related recombinants
Subtype F1 and URF_BF
CRF31_BC and related recombinants
Subtype F1 and URF_BF/URF_BCF
Analysis of the group II reveals no significant association between HIV-1 clades and the exposure category, even though some differences in the most prevalent clades circulating in HET and MSM individuals were evident. Most HIV-1 infections in the HET population were associated to subtype C (34.4%) and CRF31_BC and related recombinants (36%), followed by subtype B (21.4%) and URFs_BC (6.6%). By contrast, the majority of HIV-1 infections detected in the MSM population were related to subtype B (31.8%), followed by roughly similar frequencies of CRF31_BC and related recombinants (22.7%), URFs_BC (22.7%), and subtype C (18.2%). In both groups, the frequency of subtype F1 and URFs_BF/URFs_BCF was very low (<5%).
The current study compared the HIV-1 molecular epidemiologic scenario at two distinct time periods in the city of Porto Alegre, Southern Brazil. Patients from group I, recruited in 1998 were diagnosed between 1986 to 1998 and were likely infected around the early 1990s. Patients from group II, recruited between 2005 and 2008 were diagnosed between 2002 to 2008, and constitute a picture of the HIV-1 epidemic around the early to middle 2000s. Given the broad scope of these two sample groups, the analyses performed here have a temporal feature, shedding some light to the origins of the singular HIV-1 epidemic in Southern Brazil.
The proposed scenario is consistent with the proposed evolutionary history of subtype B and C epidemics in Brazil. Some studies suggest that the pandemic subtype B clade was introduced in Brazil around 1965–1970, at least 10–15 years earlier than subtype C [30, 31]. Thus, the pandemic subtype B is the most probable HIV-1 clade responsible for the initial AIDS cases described in Brazil, including the city of Porto Alegre, during the early 1980s, that mostly affect the MSM population (Figure 2). The proportion of AIDS cases corresponding to IDU and HET individuals started to growth in Porto Alegre around the middle 1980s (Figure 2), shortly after the estimated introduction of subtype C clade . Of note, the estimated onset date of the CRF31_BC clade (around the late 1980s) , coincides with a peak in the proportion of AIDS cases due to IDU infections in Porto Alegre (Figure 2). Since HIV transmission is very efficient through injecting equipment, super-infections with distinct viral forms are constant among IDU , which may have allowed the emergence of diverse recombinants forms, including the CRF31_BC.
Analysis of the group II reveals that by the middle 2000s most HIV-1 infections in the HET population were associated to subtype C (34%) and CRF31_BC and related recombinants (36%), while subtype B was the most common genetic variant in the MSM population (32%). However, no significant association between HIV-1 clades and the exposure infection risk were detected at this second time point. This coincides with other studies conducted in Porto Alegre between 2002 and 2009, that also have failed to find a significant association between HIV-1 clades and the exposure category [17, 21]. This lack of association may be explained by a progressive intermixing and homogenization of HIV-1 clades between different transmission networks in Porto Alegre over time, characteristic of a mature epidemic.
Comparison between groups I and II revealed a temporal increase in the proportion of subtype C (from 14% to 34%) and CRF31_BC and related recombinants (from 21% to 36%) in HET individuals. This is consistent with a previous coalescent analysis that supports an exponential expansion of subtype C and CRF31_BC clades in Porto Alegre during the 1980s and 1990s . This expansion also coincides with the growing proportion of AIDS cases in IDU and, particularly, in the HET population in Porto Alegre (Figure 2), in agreement with the process of feminization of the AIDS epidemic observed throughout the country . Similar to the AIDS Brazilian epidemic, the present study also observed an increase in heterosexual transmissions (50.6%-73.5%) and in female cases (29.0%- 40.7%) between the two groups of patients (Table 1). It is conceivable that shortly after introduction into IDU population, subtype C may have passed to HET. The early introduction of subtype C into the HET group may have promoted the successful spread of this variant in Porto Alegre, even after the decrease of HIV infections among IDU.
An increase in the proportion of subtype C (from 6% to 18%) and CRF31_BC and related recombinants (from 9% to 23%) between the two time periods was also evident among MSM individuals from Porto Alegre, thus providing evidence that dissemination of subtype C and CRF31_BC was not limited to the HET group. These results support the hypothesis that some associations between HIV-1 strains and sexual exposure categories detected in Brazil may be due to founder effects, rather than to a different efficacy of transmission of subtypes B and C through MSM or HET individuals. Analysis of the MSM population further reveals a significant increase in the frequency of URFs_BC (from 3% to 23%), that could be explained by the gradual dissemination of subtype C into this group coupled to high rates of coinfection and/or superinfection. Bisexual men and homosexual IDU might have played an essential role in this process, promoting a linkage between partially isolated transmission networks.
The overall increase in the prevalence of subtype C and BC recombinants in Porto Alegre was accompanied by a proportional decrease in the prevalence of subtype B in both HET (from 50% to 21%) and MSM (from 72% to 32%) populations. One hypothesis to explain this result is that subtype C displays a higher sexual transmissibility than subtype B. Alternatively, the dissimilar outcome of subtypes B and C in the city of Porto Alegre may be a consequence of differences in the transmission networks that promoted the initial dissemination of each subtype in the city. It is possible that subtype C and the CRF31_BC gained access to large networks of IDU and HET groups before subtype B, and this may have conditioned the subsequent dissemination of those clades among all different exposure categories. The great variation in the relative prevalence of clades B and C across different cities in Southern Brazil favors this second hypothesis. Of note, interpretations of the observed patterns should be considered with caution because of the limited sample size and the absence of IDU in the second casuistic. Future studies including larger number of HIV+ patients within different exposure categories will be necessary to confirm the observed trends.
The data presented here points to a possible introduction of subtype B and C through different transmission networks in the city of Porto Alegre. Subtype B was probably introduced through the MSM individuals while subtype C may have been introduced into the IDU network and rapidly disseminated to the HET group. This study also suggests a trend towards homogenization of HIV-1 strains across distinct exposure categories as there has been an overall increase in the prevalence of subtype C and BC recombinants in both HET and MSM populations. Understanding the mechanism responsible for the expansion and contraction of subtypes C and B, respectively, in Porto Alegre, is of paramount importance to understand the HIV-1 dynamics in this country region.
We wish to thank Vera Bongertz for English revision. We also acknowledge Regina Loureiro for patients’ recruitment and sample collection.
- Tebit DM, Arts EJ: Tracking a century of global expansion and evolution of HIV to drive understanding and to combat disease. Lancet Infect Dis 2011, 11:45–56.PubMedView Article
- World Health Organization, Global report: UNAIDS report on the global AIDS epidemic. 2010. http://www.unaids.org/globalreport/Global_report.htm.
- Ng KT, Ong LY, Takebe Y, Kamarulzaman A, Tee KK: Genome Sequence of a Novel HIV-1 Circulating Recombinant Form 54_01B from Malaysia. J Virol 2012, 86:11405–11406.PubMedView Article
- Brazilian Ministry of health: AIDS epidemiological bulletin (in portuguese). January-june 2011; ano viii, n°01. http://www.aids.gov.br/publicacao/2011/boletim_epidemiologico_2011.
- Guimarães ML, Dos Santos Moreira A, Loureiro R, Galvão-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.PubMedView Article
- Brígido LF, Franco HM, Custódio RM, Oliveira CA, P Ferreira JL, Eira M, Bergel F, Araújo F, Carvalheiro JR, Rodrigues R: Molecular characteristics of HIV type 1 circulating in São Paulo. Brazil. AIDS Res Hum Retroviruses 2005, 21:673–682.View Article
- Pedroso C, Queiroz AT, Alcântara LC, Drexler JF, Diaz RS, Weyll N, Brites C: High prevalence of primary antiretroviral resistance among HIV-1-infected adults and children in Bahia, a northeast state of Brazil. J Acquir Immune Defic Syndr 2007, 45:251–253.PubMedView Article
- Cardoso LP, Queiroz BB, Stefani MM: HIV-1 pol phylogenetic diversity and antiretroviral resistance mutations in treatment naïve patients from Central West Brazil. J Clin Virol 2009, 46:134–139.PubMedView Article
- Machado LF, Ishak MO, Vallinoto AC, Lemos JA, Azevedo VN, Moreira MR, Souza MI, Fernandes LM, Souza LL, Ishak R: Molecular epidemiology of HIV type 1 in northern Brazil: identification of subtypes C and D and the introduction of CRF02_AG in the Amazon region of Brazil. AIDS Res Hum Retroviruses 2009, 25:961–966.PubMedView Article
- Monteiro-Cunha JP, Araujo AF, Santos E, Galvao-Castro B, Alcantara LC: Lack of high-level resistance mutations in HIV type 1 BF recombinant strains circulating in northeast Brazil. AIDS Res Hum Retroviruses 2011, 27:623–631.PubMedView Article
- Soares EA, Martínez AM, Souza TM, Santos AF, Da Hora V, Silveira J, Bastos FI, Tanuri A, Soares MA: HIV-1 subtype C dissemination in southern Brazil. AIDS 2005,19(Suppl 4):S81-S86.PubMedView Article
- Brígido LF, Nunes CC, Oliveira CM, Knoll RK, Ferreira JL, Freitas CA, Alves MA, Dias C, Rodrigues R, Program RC: HIV type 1 subtype C and CB Pol recombinants prevail at the cities with the highest AIDS prevalence rate in Brazil. AIDS Res Hum Retroviruses 2007, 23:1579–1586.PubMedView Article
- Rodrigues R, Scherer LC, Oliveira CM, Franco HM, Sperhacke RD, Ferreira JL, Castro SM, Stella IM, Brigido LF: Low prevalence of primary antiretroviral resistance mutations and predominance of HIV-1 clade C at polymerase gene in newly diagnosed individuals from south Brazil. Virus Res 2006, 116:201–207.PubMedView Article
- Gräf T, Passaes CP, Ferreira LG, Grisard EC, Morgado MG, Bello G, Pinto AR: HIV-1 genetic diversity and drug resistance among treatment naïve patients from Southern Brazil: an association of HIV-1 subtypes with exposure categories. J Clin Virol 2011, 51:186–191.PubMedView Article
- Raboni SM, Almeida SM, Rotta I, Ribeiro CE, Rosario D, Vidal LR, Nogueira MB, Riedel M, Winhescki MG, Ferreira KA, Ellis R: Molecular epidemiology of HIV-1 clades in Southern Brazil. Mem Inst Oswaldo Cruz 2010, 105:1044–1049.PubMedView Article
- Santos AF, Schrago CG, Martinez AM, Mendoza-Sassi R, Silveira J, Sousa TM, Lengruber RB, Soares EA, Sprinz E, Soares MA: Epidemiologic and evolutionary trends of HIV-1 CRF31_BC-related strains in southern Brazil. J Acquir Immune Defic Syndr 2007, 45:328–333.PubMed
- Dias CF, Nunes CC, Freitas IO, Lamego IS, Oliveira IM, Gilli S, Rodrigues R, Brigido LF: High prevalence and association of HIV-1 non-B subtype with specific sexual transmission risk among antiretroviral naïve patients in Porto Alegre, RS, Brazil. Rev Inst Med Trop Sao Paulo 2009, 51:191–196.PubMed
- de Medeiros RM, Junqueira DM, Matte MC, Barcellos NT, Chies JA, Matos Almeida SE: Co-circulation HIV-1 subtypes B, C, and CRF31_BC in a drug-naïve population from Southernmost Brazil: analysis of primary resistance mutations. J Med Virol 2011, 83:1682–1688.PubMedView Article
- Toledo PV, Carvalho DS, Rossi SG, Brindeiro R, de Queiroz-Telles F: Genetic diversity of human immunodeficiency virus-1 isolates in Paraná. Brazil. Braz J Infect Dis 2010, 14:230–236.View Article
- Silva MM, Telles FQ, da Cunha CA, Rhame FS: HIV subtype, epidemiological and mutational correlations in patients from Paraná. Brazil. Braz J Infect Dis 2010, 14:495–501.View Article
- Simon D, Béria JU, Tietzmann DC, Carli R, Stein AT, Lunge VR: Prevalence of HIV-1 subtypes in patients of an urban center in Southern Brazil. Rev Saude Publica 2010, 44:1094–1101.PubMedView Article
- Brazilian Ministry of health: Coordenação DST/AIDS, Brasil. http://www.aids.gov.br/pagina/tabulacao-de-dados.
- Silveira J, Santos AF, Martínez AM, Góes LR, Mendoza-Sassi R, Muniz CP, Tupinambás U, Soares MA, Greco DB: Heterosexual transmission of human immunodeficiency virus type 1 subtype C in southern Brazil. J Clin Virol 2012, 54:36–41.PubMedView Article
- Passaes CP, Bello G, Lorete RS, Matos Almeida SE, Junqueira DM, Veloso VG, Morgado MG, Guimarães 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- Ray S: Simplot [computer program] version 3.5.1. http://sray.med.som.jhmi.edu/scroftware/.
- 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.PubMedView Article
- Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered 2001, 92:371–373.PubMedView Article
- Bello G, Guimarães ML, Morgado MG: Evolutionary history of HIV-1 subtype B and F infections in Brazil. AIDS 2006, 20:763–768.PubMedView Article
- Bello G, Eyer-Silva WA, Couto-Fernandez JC, Guimarães 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.PubMedView Article
- 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.PubMedView Article
- Bergenstrom AM, Abdul-Quader AS: Injection drug use, HIV and the current response in selected low-income and middle-income countries. AIDS 2010,24(Suppl 3):S20-S29.PubMedView Article
- Bello G, Guimarães ML, Passaes CP, Matos 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.PubMedView Article
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.