HIV-1 genetic diversity and antiretroviral drug resistance among individuals from Roraima state, northern Brazil

Abstract

The HIV-1 epidemic in Brazil has spread towards the Northern country region, but little is known about HIV-1 subtypes and prevalence of HIV strains with resistance mutations to antiretrovirals in some of the Northern states. HIV-1 protease (PR) and reverse transcriptase (RT) sequences were obtained from 73 treatment-naive and -experienced subjects followed between 2013 and 2014 at a public health reference unit from Roraima, the northernmost Brazilian state. The most prevalent HIV-1 clade observed in the study population was the subtype B (91%), followed by subtype C (9%). Among 12 HIV-1 strains from treatment-naïve patients, only one had a transmitted drug resistance mutation for NNRTI. Among 59 treatment-experienced patients, 12 (20%) harbored HIV-1 strains with acquired drug resistance mutations (ADRM) that reduce the susceptibility to two classes of antiretroviral drugs (NRTI and NNRTI or NRTI and PI), and five (8%) harbored HIV-1 strains with ADRM that reduced susceptibility to only one class of antiretroviral drugs (NNRTI or PI). No patients harboring HIV strains with reduced susceptibility to all three classes of antiretroviral drugs were detected. A substantial fraction of treatment-experienced patients with (63%) and without (70%) ADRM had undetectable plasma viral loads (<40 copies/ml) at the time of sampling. Among treatment-experienced with plasma viral loads above 2,000 copies/ml, 44% displayed no ADRM. This data showed that the HIV-1 epidemic in Roraima displayed a much lower level of genetic diversity and a lower prevalence of ADRM than that described in other Brazilian states.

About the Authors:

André de Lima Guerra Corado

Contributed equally to this work with: André de Lima Guerra Corado, Gonzalo Bello

* E-mail: [email protected] (ALGC); [email protected] (FN)

Affiliation: Instituto Leônidas e Maria Deane, FIOCRUZ, Manaus, Amazonas, Brazil

Gonzalo Bello

Contributed equally to this work with: André de Lima Guerra Corado, Gonzalo Bello

Affiliation: Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, FIOCRUZ, Rio de Janeiro, Brazil

Renato Augusto Carvalho Leão

Affiliation: Laboratório Central de Saúde Pública de Roraima (LACEN/RR), Boa Vista, Roraima, Brazil

Fabiana Granja

Affiliation: Laboratório de Biologia Molecular, Universidade Federal de Roraima (UFRR), Boa Vista, Roraima, Brazil

Felipe Gomes Naveca

* E-mail: [email protected] (ALGC); [email protected] (FN)

Affiliation: Instituto Leônidas e Maria Deane, FIOCRUZ, Manaus, Amazonas, Brazil

Introduction

In the early of 1980s the first case of HIV-1 infection was officially diagnosed in Brazil and by the end of 2014, about 780,000 people were living with HIV/AIDS in this country [1]. During the first 15 years of the epidemic, most AIDS cases were concentrated in the capital cities of the Southern and Southeastern regions, but since 1995 there was a clear expansion of the AIDS cases to the capital and inner cities of the Northeastern, Central-Western and Northern Brazilian regions [1]. HIV-1 subtype B accounts for most (70%) of HIV infections in Brazil, although other genetic variants including subtype F1, subtype C, and intersubtype recombinants (B/F and B/C) can also be found at relatively high prevalences around the country [2–8].

Nearly half (52%) of Brazilian people living with HIV/AIDS are under antiretroviral therapy (ART) [1]. This has substantially reduced the AIDS-related morbidity and mortality, but also has led to a substantial increase in the emergence and transmission of drug resistant viruses in the Brazilian population [9]. Nation-wide studies of HIV-infected treatment-naive individuals estimate that the overall rate of transmitted drug resistance mutations (TDRM) to nucleoside reverse transcriptase inhibitors (NRTI), non-nucleoside reverse transcriptase inhibitors (NNRTI) and/or protease inhibitors (PI) in Brazil increased from 6.6% to 12.2% over the past 10–15 years [2–5,10]. In the same period, the prevalence of acquired drug resistance mutations (ADRM) to NRTI, NNRTI and/or PI remained persistently high (>80%) among ART-experienced patients from Brazil [11–19].

Brazil is a continental sized country and information about the prevalence of HIV-1 subtypes and antiviral drug resistance mutations in some locations away from the most populated urban centers, particularly in the Northern Brazilian region, are scarce. Roraima is the northernmost and least populated state of Brazil, with about 500,000 inhabitants in 2015, mostly concentrated (65%) in the state capital Boa Vista. Roraima is bordered by Venezuela to the north and northwest, Guyana to the east, the state of Amazonas to the south and southwest and the state of Para to the southeast. By June 2015, 1,821 AIDS cases had been cumulatively reported in Roraima since the first identification of AIDS in the state in 1988, most of them among heterosexuals (56%) and male (66%) individuals [1].

Roraima represents an interesting location for molecular surveillance studies of the HIV epidemic because of its geographic localization and AIDS incident rate trends. Roraima maintains a very intense movement of people with Venezuela and Guyana owing to activities such as tourism, sexual exploitation and mining [20], nowadays increased by the economic crisis in Venezuela, supporting for the possibility of an intense HIV flux with those neighboring countries. Moreover, the incidence rate of AIDS cases in Roraima (29.2/100,000 inhabitants) and its capital Boa Vista (37.8/100,000 inhabitants) was well above the national mean (19.7 cases per 100.000 inhabitants) in 2014 [1]. Roraima currently figures among the top five states with the country's largest incidence rates of AIDS cases, supporting the dynamic nature of the HIV epidemic in this Brazilian state.

The only HIV molecular epidemiologic study performed in Roraima revealed that subtype B was the predominant viral strain (86%), followed by BF recombinant variants (14%) [8]. That study, however, analyzed a very low number of patients (n = 36) and with no description about the prevalence of drug resistance mutations among the study population. The objective of this study was to perform a cross-sectional study of the HIV-1 subtypes and rate of TDRM and ADRM detected among HIV-infected patients from different regions of the Roraima state.

Materials and methods

Study patients

We collected 5ml of peripheral blood samples from 73 HIV-1-infected persons that attended the Public Health Central Laboratory in Boa Vista/Roraima (LACEN-RR), a reference unit from the Brazilian Ministry of Health that receives samples for monitoring of CD4+ T cell count and plasma viral load, between January 2013 and February 2014. Inclusion criteria were adult patients (≥18 years of age) of either sex with a recent or chronic diagnosis of HIV-1 infection, and of any treatment status. Our sampling represents more than 10% of the HIV-infected population in Roraima in this period and covered eight out of 15 municipalities in the Roraima state. All patients that participated in this study were informed of the procedures and signed the informed consent. This study was approved by the Ethics Research Committee of the Federal University of Roraima (COEP/UFRR) (CAAE number: 15629013.8.0000.5302).

HIV-1 DNA amplification and sequencing

Blood samples collected in Boa Vista were transported to the Instituto Leônidas e Maria Deane, the FIOCRUZ unit in Manaus, for HIV amplification and sequencing. Genomic DNA was extracted from 200 μl of total blood using the automated system QIAcube, with QIAamp DNABlood Kit (QIAGEN, Germany), following manufacturer’s recommendations. The entire protease (PR) and part of the reverse transcriptase (RT) (nucleotides 2253–3272 relative to the HXB2 clone) gene was amplified and sequenced as previously described [21]. Sequences were assembled using the Geneious software package version 7.1.4 (Biomatters Limited, Auckland, New Zealand).

HIV-1 subtyping

HIV-1 subtypes were initially determined with the REGA HIV-1 Subtyping Tool 3.0 software [22] and later confirmed by phylogenetic and bootscanning analyses with HIV-1 reference sequences of subtypes A–D, F–H, J, and K retrieved from the Los Alamos HIV database (www.hiv.lanl.gov). Maximum-likelihood (ML) phylogenetic trees were reconstructed with the PhyML 3.0 program [23,24] to verify the clustering of Brazilian HIV-1 sequences with subtypes, using the approximate likelihood-ratio test (aLRT) [25] to estimate the reliability of the obtained tree topology. All sequences were also subjected to bootscanning analyses using the SimPlot 3.5.1 software [26] to identify possible recombination breakpoints. Bootstrap values supporting branching with HIV-1 reference sequences were determined in Neighbor Joining (NJ) trees, based on 100 resamplings, with a 250 nt sliding window moving in steps of 10 bases.

HIV-1 drug resistance analyses

The TDRM and ADRM were identified using the Calibrated Population Resistance (CPR) tool [27] (http://cpr.stanford.edu/cpr.cgi) and the Genotypic Resistance Interpretation Algorithm of the HIVdb program [28] (http://sierra2.stanford.edu/sierra/servlet/JSierra), respectively, both available through the Stanford University HIV Drug Resistance Database. The HIVdb program was also used to infer the resistance profile of the HIV-1 sequences.

Nucleotide sequence accession numbers

HIV-1 pol (PR/RT) sequences from Roraima were deposited in GenBank under accession numbers KX443015-KX443087.

Results

The epidemiological characteristics of the study population are described in Table 1 (individual data for each patient is described in the S1 Table). Among the 73 individuals included in this study, 12 were drug naïve and 61 were under ART. The period of diagnosis of the HIV infection ranged between 1995 and 2013. Most patients were males (57.5%) and the overall median age was 37 years (21–57 years range), with similar values for both males (37 years, 24–57 years range) and females (37 years, 21–56 years range). More male patients declared to be heterosexual (57%) than homosexual (33%) or bisexual (10%), whereas only one female declared to be homosexual. The major risk factor for HIV infection identified was unprotected sex (81%), followed by iatrogenic transmission/intravenous drug use (12%), while for 7% of the patients the route of transmission was unknown. The median of CD4 T cell counts was 448 cells/μl (74–1,822 cells/μl range). The CD4 T cell counts was <350 cells/μl for 36%, between 350 and 500 cells/μl for 22%, and >500 cells/μl for 33% of the participants. The median of plasma viral load was 54,614 copies/ml (<40–1,383,017 copies/ml range). Plasma viral loads were below limit of detection (40 copies/ml) for 58%, between 51 and 10,000 copies/ml for 19%, and over 10,000 copies/ml for 16% of the participants. Most patients (79%) were from the capital, Boa Vista. The remaining ones were from other municipalities of the Roraima state (16%) or neighboring cities outside Roraima (5%), including one patient from Manaus (Amazonas state), two patients from St. Helena de Uairén (Venezuela) and one patient from Lethem (Guyana).

[Figure omitted. See PDF.]

Table 1. Epidemiological and clinical characteristics of HIV-1-infected patients attended at the Public Health Central Laboratory in Boa Vista, Roraima.

https://doi.org/10.1371/journal.pone.0173894.t001

All patients recruited had HIV-1 pol (PR/RT) genes amplified and sequenced, and nearly all were characterized as subtype B (91%), except one classified as subtype C (9%) (Fig 1). Among the 12 treatment-naïve patients analyzed, only one harbored a TDRM to NNRTI (G190A), whereas no TDRM to NRTI or PI were detected in the study population. Two out of 61 HIV-1 PR/RT sequences obtained from patients under ART were excluded due to the presence of APOBEC3GF-mediated G to A mutations. Among the remaining 59 treatment-experienced patients, 17 (29%) harbored HIV-1 strains with ADRM that conferred some level of resistance to NNRTI, NRTI and/or PI (Table 2). Among these, 14 (24%) patients had ADRM to NNRTI (being K103N and G190A the most common), 12 (20%) patients had ADRM to NRTI (being M184V the most frequently found) and three (5%) patients had ADRM to PI. Five (8%) patients showed ADRM to only one class of antiretroviral drugs (NNRTI or PI), 10 (17%) patients showed ADRM to both NRTI and NNRTI, and two (3%) patients showed ADRM to both NRTI and PI. No patients showed ADRM to all three classes of antiretroviral drugs. The median plasma viral load of patients without ADRM (2,917 copies/ml) was much lower than of those with ADRM (130,360 copies/ml range), although the proportion of patients without ADRM that displayed undetectable plasma viral loads (70%) was similar to the corresponding proportion of patients with ADRM (63%) and the prevalence of ADRM among patients with detectable plasma viral loads (32%) was only slightly higher than among patients with undetectable viral loads (25%). HIV-1 sequences containing DRM were widely dispersed in the phylogenetic tree among sequences with no DRM (Fig 1). Only two highly supported (aLRT > 0.90) pairs of related sequences with ADRM were identified (Fig 1), but sequences within each pair displayed a different set of ADRM (Table 2).

[Figure omitted. See PDF.]

Table 2. List of mutations known to confer reduced susceptibility to antiretroviral agents among ART-experienced HIV-1-infected patients attended at the Public Health Central Laboratory in Boa Vista, Roraima.

https://doi.org/10.1371/journal.pone.0173894.t002

[Figure omitted. See PDF.]

Fig 1. Maximum likelihood phylogenetic tree of HIV-1 pol (PR/RT) sequences obtained from ART naive and experienced patients from the Roraima state.

HIV-1 pol sequences from Roraima (red terminal branches) were combined with reference sequences (black terminal branches) of all HIV-1 group M subtypes (A–D, F–H, J, and K). Shaded boxes indicate the position of HIV-1 subtypes detected in Roraima. Asterisks indicate the position of sequences with ADRM (dark) and TDRM (blue). Only aLRT values at the branch of subtype B, subtype C, and the two highly supported clades of sequences with ADRM are shown. Trees are rooted at midpoint and the branch lengths are drawn to scale with the bar at the bottom indicating nucleotide substitutions per site.

https://doi.org/10.1371/journal.pone.0173894.g001

Discussion

The current analysis showed that the HIV-1 epidemic in the northernmost Brazilian state of Roraima is almost exclusively driven by the subtype B. Among the 73 HIV-1 pol (PR/RT) gene sequences analyzed, 72 were classified as subtype B and one as subtype C. A previous study analyzed 36 HIV-1 pol (PR/RT) gene sequences from Roraima and classified 31 sequences as subtype B and five sequences as BF recombinants [8]. When both studies are combined, we obtained that the predominant HIV-1 clade circulating in Roraima was subtype B (94%), followed by BF recombinants (5%) and subtype C (1%).

These results indicate that the HIV-1 molecular epidemiologic profile observed in Roraima is much more homogenous than that described in any other Brazilian state, even from the Northern region [8,18,29–33], but comparable to that described in the neighboring country of Venezuela [34,35]. Roraima is one of the most isolated Brazilian states, being connected by land only to the state of Amazonas. On the other hand, Boa Vista (the capital of the Roraima state) is located close to the Venezuelan and Guyanese cities of St. Helena de Uairén (230 km) and Lethem (114 km), respectively. Notably, among the 73 patients included in our study and that attended the Public Health Central Laboratory in Boa Vista, two patients were from St. Helena de Uairén and one patient from Lethem, reflecting the intense human mobility across the Roraima borders with Venezuela and Guyana. Consistent with this, we recently described a singular high prevalence of non-pandemic subtype B variants of Caribbean origin in Roraima [36]. These data suggest that the HIV-1 epidemic in Roraima is more closely related to epidemics in neighboring South American countries (Venezuela and Guyana) than to epidemics in other Brazilian regions.

Among the 12 HIV-1 pol sequences from drug naïve patients here analyzed, only one displayed a TDRM to NNRTI, whereas no TDRM to NRTI or PI were detected. A much higher level of ADRM (28%), by contrast, was observed among subjects under ART included in our study. Most HIV-1 isolates with ADRM harbored some level of resistance to both NRTI and NNRTI (60%), whereas others displayed resistance to only NNRTI (23%), to both PI and NRTI (12%), or to only PI (6%). No patients with ADRM to all three classes of antiretroviral drugs were detected in this study. The phylogenetic analysis reveals no highly supported clades of HIV-1 sequences with the same set of ADRM, thus supporting that DRM observed among patients under ART most probably evolved independently after initiation of therapy and were not directly transmitted among subjects.

The prevalence of ADRM here detected in Roraima (28%) was relatively low compared to the higher levels of resistance to NRTI, NNRTI and/or PI (>80%) reported in other Brazilian states from all country regions [11–19]. It could be argued that most previous studies of ADRM conducted in Brazil analyzed populations of HIV-1-infected individuals failing ART and that the relative low rate of ADRM observed in Roraima may be explained by a low frequency of virological failure among treated patients. However, the estimated prevalence of ADRM in Roraima remains comparatively lower to that described in other Brazilian states even if we only considered those ART-experienced patients with detectable viral load (32%), or those patients with plasma viral loads >1,000 copies/ml (56%) that is the viral load cut-off for HIV resistance genotyping according to the guidelines of the Brazilian Ministry of Health (http://www.aids.gov.br/pagina/2012/51642).

It is also interesting to note that although the median plasma viral load of patients without ADRM (2,917 copies/ml) was much lower than of those with ADRM (130,360 copies/ml range), not all patients harboring ADRM displayed evidence of treatment failure. Indeed, a significant proportion of patients with ADRM (63%) displayed undetectable plasma viral loads at the time of sampling, and the prevalence of ADRM among patients with detectable plasma viral loads (32%) was only slightly higher than among patients with undetectable viral loads (25%). On the other hand, among the ART-experienced patients with plasma viral loads >1,000 copies/ml, 44% harbored HIV strains fully susceptible to antiretrovirals. Thus, a significant proportion of the apparent therapeutic failures among HIV-1-infected patients from Roraima could probably be attributed to problems of sub-optimal therapeutic adherence rather than to selection of drug-resistance strains.

Although this study suggests a low prevalence of both TDRM and ADRM among HIV-infected individuals from the Roraima state, these results should be interpreted with caution due to the relatively low size of the study populations. New surveillance studies of both untreated and ART-treated patients will be necessary to obtain more accurate estimates of the prevalence of TDRM and ADRM in this Brazilian state. The application of more specific questionnaires will be necessary to quantify the relative contribution of therapeutic success and sub-optimal therapeutic adherence to the low prevalence of ADRM observed in Roraima. Finally, information about previous and current ART regimens used by the patients will be also of paramount importance to fully explain the high rate of ADRM observed among patients with undetectable plasma viremia.

In conclusion, our study indicates that the HIV-1 epidemic in Roraima is mostly dominated by the B subtype and that the epidemiological relevance of other subtypes and recombinants forms is much lower than that described in other Brazilian states. Our findings also revealed a relatively high frequency of undetectable viral load and a relatively low prevalence of ADRM among treatment-experienced individuals from Roraima. As the HIV-1 epidemic and the use of ART spreads towards the innermost parts of Brazil, continued molecular surveillance studies will be essential to monitor temporal changes in the prevalence of drug resistant strains among treatment naïve and experienced patients in those areas.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Viral load, CD4 information and genbank access number of the sequences of each patient of HIV-1-infected patients attended at the Public Health Central Laboratory in Boa Vista, Roraima.

https://doi.org/10.1371/journal.pone.0173894.s001

(DOCX)

Acknowledgments

The authors thank the Program for Technological Development in Tools for Health-PDTIS-FIOCRUZ and the Instituto Lêonidas e Maria Deane, for use of its facilities (ILMD sequencing platform), the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes), and the Governo do Estado de Roraima e Laboratório Central de Saúde Pública do Estado de Roraima (LACEN-RR). The authors also thanks to Dr Vera Bongertz for critical review of the manuscript.

Author Contributions

1. Conceptualization: ALGC FGN FG.

2. Data curation: ALGC RACL FGN.

3. Formal analysis: ALGC GB FGN.

4. Funding acquisition: FGN.

5. Investigation: ALGC FGN GB RACL FG.

6. Methodology: ALGC GB FGN.

7. Project administration: ALGC FGN.

8. Resources: ALGC FGN GB RACL FG.

9. Supervision: ALGC GB FGN.

10. Visualization: GB.

11. Writing – original draft: GB.

12. Writing – review & editing: ALGC FGN GB.

Citation: Corado AdLG, Bello G, Leão RAC, Granja F, Naveca FG (2017) HIV-1 genetic diversity and antiretroviral drug resistance among individuals from Roraima state, northern Brazil. PLoS ONE 12(3): e0173894. https://doi.org/10.1371/journal.pone.0173894

References

1. Brazilian (2015) Ministry of Health. AIDS Epidemiological Bulletin [in Portuguese] January-June 2015; Ano IV, n° 01 http://wwwaidsgovbr/sites/default/files/anexos/publicacao/2015/58534/boletim_aids_11_2015_web_pdf_19105pdf.

2. Brindeiro RM, Diaz RS, Sabino EC, Morgado MG, Pires IL, Brigido L, et al. (2003) Brazilian Network for HIV Drug Resistance Surveillance (HIV-BResNet): a survey of chronically infected individuals. Aids 17: 1063–1069. pmid:12700457

3. Inocencio LA, Pereira AA, Sucupira MC, Fernandez JC, Jorge CP, Souza DF, et al. (2009) Brazilian Network for HIV Drug Resistance Surveillance: a survey of individuals recently diagnosed with HIV. J Int AIDS Soc 12: 20. pmid:19765271

4. Alencar CS, Sabino EC, Carvalho SM, Leao SC, Carneiro-Proietti AB, Capuani L, et al. (2013) HIV genotypes and primary drug resistance among HIV-seropositive blood donors in Brazil: role of infected blood donors as sentinel populations for molecular surveillance of HIV. J Acquir Immune Defic Syndr 63: 387–392. pmid:23507660

5. de Moraes Soares CM, Vergara TR, Brites C, Brito JD, Grinberg G, Caseiro MM, et al. (2014) Prevalence of transmitted HIV-1 antiretroviral resistance among patients initiating antiretroviral therapy in Brazil: a surveillance study using dried blood spots. J Int AIDS Soc 17: 19042. pmid:25249214

6. Stefani MM, Pereira GA, Lins JA, Alcantara KC, Silveira AA, Viegas AA, et al. (2007) Molecular screening shows extensive HIV-1 genetic diversity in Central West Brazil. J Clin Virol 39: 205–209. pmid:17537671

7. Pessoa R, Loureiro P, Esther Lopes M, Carneiro-Proietti AB, Sabino EC, Busch MP, et al. (2016) Ultra-Deep Sequencing of HIV-1 near Full-Length and Partial Proviral Genomes Reveals High Genetic Diversity among Brazilian Blood Donors. PLoS One 11: e0152499. pmid:27031505

8. da Costa CM, Costa de Oliveira CM, Chehuan de Melo YF, Delatorre E, Bello G, Couto-Fernandez JC (2016) High HIV-1 Genetic Diversity in Patients from Northern Brazil. AIDS Res Hum Retroviruses 32: 918–922. pmid:27091699

9. Hacker MA, Kaida A, Hogg RS, Bastos FI (2007) The first ten years: achievements and challenges of the Brazilian program of universal access to HIV/AIDS comprehensive management and care, 1996–2006. Cad Saude Publica 23 Suppl 3: S345–359.

10. Sprinz E, Netto EM, Patelli M, Lima JS, Furtado JJ, da Eira M, et al. (2009) Primary antiretroviral drug resistance among HIV type 1-infected individuals in Brazil. AIDS Res Hum Retroviruses 25: 861–867. pmid:19689190

11. Tanuri A, Caridea E, Dantas MC, Morgado MG, Mello DL, Borges S, et al. (2002) Prevalence of mutations related to HIV-1 antiretroviral resistance in Brazilian patients failing HAART. J Clin Virol 25: 39–46. pmid:12126720

12. Bahia F, Pedroso C, Netto EM, Figueiredo R, Pinto Neto L, Brites C (2004) Evaluation of the genotypic pattern of HIV-1 resistance in AIDS patients failing antiretroviral therapy. Braz J Infect Dis 8: 281–289. pmid:15565258

13. Couto-Fernandez JC, Silva-de-Jesus C, Veloso VG, Rachid M, Gracie RS, Chequer-Fernandez SL, et al. (2005) Human immunodeficiency virus type 1 (HIV-1) genotyping in Rio de Janeiro, Brazil: assessing subtype and drug-resistance associated mutations in HIV-1 infected individuals failing highly active antiretroviral therapy. Mem Inst Oswaldo Cruz 100: 73–78.

14. Cavalcanti AM, Lacerda HR, Brito AM, Pereira S, Medeiros D, Oliveira S (2007) Antiretroviral resistance in individuals presenting therapeutic failure and subtypes of the human immunodeficiency virus type 1 in the Northeast Region of Brazil. Mem Inst Oswaldo Cruz 102: 785–792. pmid:17992369

15. Cardoso LP, Stefani MM (2009) High level of multidrug resistance mutations in HIV type 1 pol gene and resistance-associated mutations to enfuvirtide (T-20) among antiretroviral-experienced patients from central Brazil. AIDS Res Hum Retroviruses 25: 943–950. pmid:19792869

16. Toledo PV, Carvalho DS, Romagnoli L, Marcinko G, Cunha CA, Souza MN, et al. (2010) HIV-1 genotypic resistance profile of patients failing antiretroviral therapy in Parana, Brazil. Braz J Infect Dis 14: 360–371. pmid:20963321

17. Figueiredo Saad F, Kaminami Morimoto H, Wiechmann SL, Bonametti AM, Ueda LT, Matsuo T, et al. (2010) Frequency and diversity of human immunodeficiency virus type 1 mutations associated with antiretroviral resistance among patients from Southern Brazil failing highly active antiretroviral therapy (HAART). Int J Mol Med 26: 585–593. pmid:20818500

18. Lopes CA, Soares MA, Falci DR, Sprinz E (2015) The Evolving Genotypic Profile of HIV-1 Mutations Related to Antiretroviral Treatment in the North Region of Brazil. Biomed Res Int 2015: 738528. pmid:26543866

19. Diaz RS, Inocencio LA, Sucupira MC, Pereira AA, Hunter J, Ferreira JE, et al. (2015) The Virological and Immunological Characteristics of the HIV-1-Infected Population in Brazil: From Initial Diagnosis to Impact of Antiretroviral Use. PLoS One 10: e0139677. pmid:26509441

20. Diniz AMA, dos Santos RO (2005) O vertiginoso crescimento populacional de Roraima e seus impactos socioambientais (in Portuguese). Caderno de Geografia 15: 23–44.

21. Bello G, Passaes CP, Guimaraes ML, Lorete RS, Matos Almeida SE, Medeiros RM, et al. (2008) Origin and evolutionary history of HIV-1 subtype C in Brazil. AIDS 22: 1993–2000. pmid:18753928

22. de Oliveira T, Deforche K, Cassol S, Salminen M, Paraskevis D, Seebregts C, et al. (2005) An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics 21: 3797–3800. pmid:16076886

23. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59: 307–321. pmid:20525638

24. Guindon S, Lethiec F, Duroux P, Gascuel O (2005) PHYML Online—a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res 33: W557–559. pmid:15980534

25. Anisimova M, Gascuel O (2006) Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol 55: 539–552. pmid:16785212

26. Ray S Simplot v2.5.0. http://wwwmedjhuedu/deptmed/sray/download/.

27. Gifford RJ, Liu TF, Rhee SY, Kiuchi M, Hue S, Pillay D, et al. (2009) The calibrated population resistance tool: standardized genotypic estimation of transmitted HIV-1 drug resistance. Bioinformatics 25: 1197–1198. pmid:19304876

28. Liu TF, Shafer RW (2006) Web resources for HIV type 1 genotypic-resistance test interpretation. Clin Infect Dis 42: 1608–1618. pmid:16652319

29. Dos Anjos Silva L, Divino F, da Silva Rego MO, Lima Lopes IG, Nobrega Costa CM, da Silva Pereira FC, et al. (2016) HIV-1 Genetic Diversity and Transmitted Drug Resistance in Antiretroviral Treatment-Naive Individuals from Amapa State, Northern Brazil. AIDS Res Hum Retroviruses 32: 373–376. pmid:26529282

30. Carvalho BC, Cardoso LP, Damasceno S, Stefani MM (2011) Moderate prevalence of transmitted drug resistance and interiorization of HIV type 1 subtype C in the inland North State of Tocantins, Brazil. AIDS Res Hum Retroviruses 27: 1081–1087. pmid:21417758

31. Cunha LK, Kashima S, Amarante MF, Haddad R, Rodrigues ES, Silva KL, et al. (2012) Distribution of human immunodeficiency virus type 1 subtypes in the State of Amazonas, Brazil, and subtype C identification. Braz J Med Biol Res 45: 104–112. pmid:22249428

32. Macêdo O, Ferreira LM, Lopes CAF, de Sousa RCM, Araujo JRM, Vasconcelos PFC (2012) Distribution of HIV-1 subtypes in patients with HAARTtherapeutic failure in the States of Pará and Amazonas, Brazil: 2002 to 2006. Rev Pan-Amaz Saude 3: 11–16.

33. Machado LF, Ishak MO, Vallinoto AC, Lemos JA, Azevedo VN, Moreira MR, et al. (2009) 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 25: 961–966. pmid:19795985

34. Delgado E, Leon-Ponte M, Villahermosa ML, Cuevas MT, Deibis L, Echeverria G, et al. (2001) Analysis of HIV type 1 protease and reverse transcriptase sequences from Venezuela for drug resistance-associated mutations and subtype classification: a UNAIDS study. AIDS Res Hum Retroviruses 17: 753–758. pmid:11429115

35. Rangel HR, Garzaro DJ, Torres JR, Castro J, Suarez JA, Naranjo L, et al. (2009) Prevalence of antiretroviral drug resistance among treatment-naive and treated HIV-infected patients in Venezuela. Mem Inst Oswaldo Cruz 104: 522–525. pmid:19547882

36. Divino F, de Lima Guerra Corado A, Gomes Naveca F, Stefani MM, Bello G (2016) High Prevalence and Onward Transmission of Non-Pandemic HIV-1 Subtype B Clades in Northern and Northeastern Brazilian Regions. PLoS One 11: e0162112. pmid:27603317