BioMed Research International
Volume 2021 (2021), Article ID 5567332, 9 pages
https://doi.org/10.1155/2021/5567332

Research Article

High Detection Rate of HIV Drug Resistance Mutations among Patients Who Fail Combined Antiretroviral Therapy in Manaus, Brazil

Yury Oliveira Chaves,¹ Flávio Ribeiro Pereira,² Rebeca de Souza Pinheiro,¹ Diego Rafael Lima Batista,¹ Antônio Alcirley da Silva Balieiro,¹ Marcus Vinicius Guimarães de Lacerda,1,2 Paulo Afonso Nogueira,1,2 and Monick Lindenmeyer Guimarães³

¹Laboratório de Diagnóstico e Controle de Doenças Infecciosas na Amazônia (DCDIA), Instituto Leônidas e Maria Deane (ILMD)–Fiocruz Amazônia, Manaus CEP 69.057-070, Brazil
²Fundação de Medicina Tropical Dr. Heitor Vieira Dourado, Instituto de Pesquisa Clínica Carlos Borborema, Manaus, Amazonas CEP 69040-000, Brazil
³Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, FIOCRUZ, Rio de Janeiro, RJ CEP 21040-900, Brazil

Correspondence should be addressed to Yury Oliveira Chaves; yurychaves@gmail.com

Received 9 February 2021; Revised 10 May 2021; Accepted 27 May 2021; Published 9 June 2021

Academic Editor: Yanjin Zhang

Copyright © 2021 Yury Oliveira Chaves et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Virologic failure may occur because of poor treatment adherence and/or viral drug resistance mutations (DRM). In Brazil, the northern region exhibits the worst epidemiological scenarios for the human immunodeficiency virus (HIV). Thus, this study is aimed at investigating the genetic diversity of HIV-1 and DRM in Manaus. The cross-sectional study included people living with HIV on combined antiretroviral therapy and who had experienced virological failure during 2018-2019. Sequencing of the protease/reverse transcriptase (PR/RT) and C2V3 of the viral envelope gp120 (Env) regions was analyzed to determine subtypes/variants of HIV-1, DRMs, and tropism. Ninety-two individuals were analyzed in the study. Approximately 72% of them were male and 74% self-declared as heterosexual. Phylogenetic inference (PR/RT-Env) showed that most sequences were subtype, followed by BF1 or mosaic genomes and few F1 and sequences. Among the variants of subtype at PR/RT, 84.3% were pandemic (), and 15.7% were Caribbean (). The DRMs most frequent were M184I/V (82.9%) for nucleoside reverse transcriptase inhibitors (NRTI), K103N/S (63.4%) for nonnucleoside reverse transcriptase inhibitor (NNRTI), and V82A/L/M (7.3%) for protease inhibitors (PI). DRM analysis depicted high levels of resistance for lamivudine and efavirenz in over 82.9% of individuals; although, low (7.7%) cross-resistance to etravirine was observed. A low level of resistance to protease inhibitors was found and included patients that take atazanavir/ritonavir (16.6%) and lopinavir (11.1%), which confirms that these antiretrovirals can be used—for most individuals. The thymidine analog mutations-2 (TAM-2) resistance pathway was higher in than in . Similar results from other Brazilian studies regarding HIV drug resistance were observed; however, we underscore a need for additional studies regarding subtype variants. Molecular epidemiology studies are an important tool for monitoring the prevalence of HIV drug resistance and can influence the public health policies.

1. Introduction

Since 1996, Brazil has offered free combined antiretroviral therapy(cART) to all people living with HIV/AIDS (PLWHA) regardless of the TCD4+ cell count—and has adopted innovative actions over the years, mainly in primary care [1, 2]. In cART, the drugs used act on different viral targets, and this therapy has brought significant advances in the treatment of PLWHA, such as a fast reduction of the viral load and their more extended maintenance, followed by a significant increase in the TCD4 lymphocyte counts [3, 4]. This has led to better quality of life, with individuals achieving immunological stability and a consequent reduction in the incidence of opportunistic infections and HIV/AIDS-related deaths. Another benefit that stems from the suppression of HIV replication is the reduced possibility of HIV transmission [5, 6].

Treatment failure in HIV therapy can be related to the presence of drug resistance mutations, drug intolerance, and low adherence to therapeutic regimens causing suboptimal ARV dosage, as well as anthropogenic aspects [7, 8] such as physical, psychological, social, or cultural aspects [9]. In addition, the use of ARVs with a low genetic barrier for resistance in the first therapeutic regimen may favor the rapid selection of DRMs [10]. Virological failure is actually defined by the presence of after at least six months of cART; hence, HIV-1 genotyping is recommended in order to verify the presence of acquired DRMs [11–13]. Depending on the type of DRM, crossresistance to another ARV of the same class could be verified, as reported for thymidine analog mutations (TAMs). TAMs confers crossresistance to nucleoside analogues, potentially affecting first cART regimens, as may prevent the use of another NRTI in subsequent regimens, and could result in the replacement of the ARV class. [14–17].

Despite advances in HIV treatment in Brazil, the North and northeastern regions registered an increase in the AIDS detection rate between 2009 and 2019. In the Amazonas state, detection arose from 31.7 to 34.8 per 100,000 habitants in this period, despite of a 4.5% decrease in the AIDS mortality rate (6.7 to 6.4 per 100,000 habitants). In 2019, Manaus occupied third place in the mortality coefficient rankings (10.5 per 100,000 habitants) and fifth place in AIDS detection rate rankings (54.7 per 100,000 habitants) [18].

The genetic diversity of HIV-1 found in the northern region is little distinct from other regions in Brazil [19], characterized mainly by the prevalence of subtype (63-92%), subtype F1 (0-14%), subtype (0-6%), BF1 (4-18%), and (0-8%) [20–24]. Contrasting to other Brazilian regions, high prevalence of variant (14.4%) was verified in this region [25]. The studies about subtype variants have been focusing in epidemiological aspects; however, the possible impact in the pathogenesis or acquisition of DRM was not investigated.

This problematic epidemiological scenario, along with little data concerning HIV-1 diversity and DRMs in Manaus, led us to investigate these topics in individuals with virological failure ongoing cART schemes.

2. Material and Methods

2.1. Study Subject Characteristics and Ethical Aspects

This cross-sectional study included 100 individuals with positive HIV-1 serology, who were older than 18 years of age, undergoing antiretroviral therapy for at least six months and presenting virologic failure ( of plasma), receiving a medical request for HIV-1 genotyping following the guidelines of the Brazilian Ministry of Health at the time of the study. All individuals were included in the study and had their samples collected between 2018 and 2019.. All subjects were treated at the Tropical Medicine Foundation-Heitor Vieira Dourado (FMT-HVD), Manaus, Amazonas state, Brazil. Demographic and clinical data were collected using a structured questionnaire. A sample of 4 ml of whole blood was collected by venous pulse, and the samples were processed for plasma separation and stored at -80°C until RNA extraction. The study was approved by the Ethical Review Board, IOC/FIOCRUZ under CAAE 87171018.4.0000.5248 protocol, the FMT-HVD under approval CAAE 87171018.4.3001.0005, and all subjects gave written informed consent in accordance with the Declaration of Helsinki.

2.2. RNA and RT-PCR Extraction

Extraction of genetic materials was performed using the QIAamp RNA mini kit (Qiagen, Germany) following the manufacturer’s instructions. This was followed by reverse transcription of RNA to obtain cDNA using MMLV enzyme (Invitrogen, Carlsbad, CA) protocol with specific external reverse PCR primers for each gene target.

2.3. Amplification, Purification, and Sequencing of PR/RT and C2-V3

The pol region was amplified by nested PCR with initial primer G17S (AAAAAGGGCTGTTGTTGGTGGAATGTGGA)/MMRTR6 (TTTTACATTAGTGTGGG) and MMRTR5 (TAAATTTGATATGTCCATTG)/MMRT10 (CAGGCTAATTTTTTAGGGAA) generating a final fragment of 1478 bp (positions 2077-3574 relative to HXB2). For the env region, we used the external primers ED5 (ATGGGATCAAAGCCTAAAAGCCATGTG)/ED12 (AGTGCTTCCTGCTGCTCCCAAGAACCCAAG) and ED31(CCTCAGCCATTACACAGGCCTGTCCAAAG)/ED33 (TTACAGTAGAAAAATTCCCCTC) generating a final fragment of 563 bp (positions 6816-7380 relative to HXB2), as described by Delwart et al. and Delatorre et al. [26, 27]. The cycling protocol for both study regions was 94°C for 2 minutes for initial denaturation, followed by 35 cycles of 94°C for 30 seconds, 55°C for 30 seconds, and 72°C for 2 minutes, with a final extension of 10 minutes at 72°C.

The amplified PCR products were purified using the Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare, United Kingdom) following the manufacturer’s instructions. The products were subsequently subjected to sequencing using ABI BigDye Terminator v.3.1 Cycle Sequencing Ready (Applied Biosystems, United States). The samples were sequenced in the ABI Prism 3130 automatic sequencer (Applied Biosystems).

2.4. Sequence Analysis

The sequences were edited using SeqMan software of the DNASTAR 4.0 program and aligned by Clustal W in the MEGA 6.0 program together with corresponding reference sequences of the different HIV-1 subtypes obtained from the Los Alamos database (http://hiv.lanl.gov). The final pol (partial PR/RT) and env (partial gp120) alignments contained a 1260 nucleotide (nt) (positions 2254-3514 relative to the HXB2) and a 397 nt fragment (positions 6858-7225 relative to the HXB2), respectively.

The phylogenetic analysis was performed based on the neighbor-joining method and the Tamura-Nei substitution model. The bootstrap test with 1,000 replicates was used to estimate the confidence level of the branching on the phylogenetic tree. Possible recombinant sequences were detected via bootscan analysis, which was performed using the Simplot 3.5.1 program. These analyses were performed on a sliding window of 200 nt moving in 10 nt steps. For the characterization of the BBR variant in the C2V3 region, an amino acid sequence was obtained by the translation of the env alignment using MEGA, and all sequences showing the GWG motif at the top of loop V3 were considered. The Brazilian HIV-1 B pol sequences were aligned using the Clustal W in the MEGA program, with reference B (Pandemic or Caribbean) sequences, and the positions associated with DRMs were removed. The analysis was performed using the maximum likelihood (ML) method in the PhyML program, and the FigTree v1.4.0 program was used for confection and visualization of the phylogenetic trees.

Analysis of PR/RT resistance mutations was performed for each sequence on the Stanford HIV Drug Resistance Database website (http://hivdb.stanford.edu/). All major resistance mutations from the updated list on this site were considered (https://hivdb.stanford.edu/pages/download/resistanceMutations_handout.pdf and https://www.iasusa.org/resources/hiv-drug-resistance-mutations/). The susceptibility of each sequence to the ARV was also determined using the HIVdb Program from the Stanford HIV Drug Resistance Database. However, concerning the resistance level to the ARV only the high-level resistance was applied.

The genotypic prediction was used to infer the use of CCR5 and/or CXCR4 coreceptors. The prediction was performed based on the amino acid sequence of the V3 region of the gp120 of the envelope of each sequence using the program Geno2pheno (available at https://coreceptor.geno2pheno.org/), with a false positive cut-off rate of 10% according to the European Tropism Test guidelines.

2.5. Statistical Analyses

The data were analyzed using the GraphPad Prism software version 5 and the software, with a significance level set at . Statistical comparisons between and within groups were made using Fisher’s exact test, the Mann–Whitney test, and the Kruskal-Wallis test when appropriate.

3. Results

3.1. Epidemiological and Clinical Data

Of the 100 individuals enrolled in the study, we were able to amplifiy 95% at the PR/RT region and to sequence 92% of them. From those 92 individuals, approximately, 72% were male, with an average age of 40 (31-46) years, and approximately, 95% were residents of Manaus, the capital of the Amazonas state. Regarding sexual orientation, 74% declared themselves heterosexual, 19% as homosexual, and 6% as bisexual. The median viral load (log) was 4.4 (IQR 4-5), and the median CD4+ count was 168 cell/mm³ (IQR 50.5-342.5) at the time of sample collection. The patients from this cohort presented a mean of six years (IQR 4-7) receiving cATR, 57.6% of them were receiving their first treatment, and 42.4% had already undergone treatment changes (TARV). From the 92 (PR/RT) analyzed individuals, only 78 (84.8%) sequences were determined in the envelope region. Thus, according to the phylogenetic inference of the combined pol (RT/PR) and env regions together, it was possible to classify these 78 HIV-1 (pol/env) sequences as 68.5% subtype /, F1/F1 (1.1%), and / (1.1%), and 14.1%were mosaic genomes, for which 13% BF1 and 1.1% . However, for the remaining 14 sequences (15.2%) obtained only in the pol region, 14.2% were subtype , and 1% was subtype . Regarding the HIV-1 subtype variants, 84.3% were classified as pandemic and 15.7% as Caribbean in the pol region, and 5.1% were Brazilian variants () in env. The bioinformatic analyses were able to predict that the viral tropism found in the sequences of the env region demonstrated that about the majority (66.7%) of the 78 studied individuals sequenced at env region had the R5 virus at the time of the sample collection.

Among the 92 studied individuals, DRMs were observed in 82 sequences (89%), from those 91.5% presented DRMs associated with NRTI () and 84.1% with NNRTI () (Figure 1). All patients were taking NRTIs at the time the genotyping test was done, with 3TC being used by 92.6% of the individuals () and TDF by 84.1% (). Concerning NNRTIs, all individuals had already used EFV; however, only 63.4% (n =52/82) were currently using this ARV at the collection time (Table 1). The most frequently detected DRMs were M184I/V (68/82; 82.9%), K70E/R (16/82; 19.5%), and T215F/Y (17/82; 20.7%) to NRTI and K103N/S (51/82; 62.1%), P225H (15/82; 18.2%), and V106A/I/M (11/82; 13.4%) to NNRTI (Figure 1). Thymidine analog mutations (TAM) were detected in 35.4% of the individuals [17.1% (14/82) TAM-1 (including M41L, L210W, and T215Y) and 20.7% (17/82) TAM-2 (including D67N, K70R, T215F, and K219Q/E)].

(a)

(b)

(c)

(d)

Figure 1: Frequency of HIV-1 drug resistance mutations according to inhibitor class and TAMs detected in 82 individuals. (a) DRM to NRTI: nucleoside reverse transcriptase inhibitors (91.5%, ). (b) DRM to NNRTI: nonnucleoside reverse transcriptase inhibitor (84.1%, ). (c) DRM to PI: protease inhibitors (17%, n =14/82). (d) Thymidine analog mutations TAM [TAM-1 (), TAM-2 ()].

Table 1: Prevalence of ARV used at any time during treatment and current in use, prevalence of the most prevalent DRMs to each ARV class in relation to the ARV current in use, and percentage of individuals with high ARV resistance related to the ARV used.

ARV
ARV used, (%)
ARV current in use, (%)
Drug resistance mutation, (%)
ARV high resistance level, (%)

M184I/V
K70E/R
T215F/Y

Total NRTI
82/82 (100)
82/82 (100)
68/82 (82.9)
16/82 (19.5)
17/82 (20.7)
70/82 (85.3)

3TC
82/82 (100)
76/82 (92.6)

68/82 (82.9)

TDF
70/82 (85.3)
69/82 (84.1)

6/70 (8.5)

AZT
41/82 (50.0)
5/82 (6.0)

9/41 (21.9)

K103N/S
V106A/I/M
P225H

Total NNRTI
71/82 (86.5)
48/82 (58.5)
51/71 (71.8)
11/71 (15.4)
15/71 (21.1)
64/71 (90.1)

EFV
71/82 (86.5)
52/82 (63.4)

64/71 (90.1)

V82A/L/M
I54L/M/V
M46L/I

Total PI
46/82 (50.0)
26/82 (31.7)
5/46 (10.8)
5/46 (10.8)
4/46 (8.6)
6/46(13.0)

LPV
18/82 (21.9)
1/82 (1.2)

2/18 (11.1)

ATZ/r
24/82 (29.2)
19/82 (23.1)

4/24 (16.6)

DRV
6/82 (7.3)
4/82 (4.8)

0 (0)

ARVs: 3TC: lamivudine; EFV: efavirenz; TDF: tenofovir; AZT: zidovudine; AZT/r: atazanavir/ritonavir; DRV: darunavir; LPV: lopinavir; PI: protease inhibitors. ARV resistance inputted according to high resistance level in the Stanford HIV Drug Resistance database (https://hivdb.stanford.edu/).

Only half of the individuals had used PI (), only 17% of individuals presented DRMs associated with protease inhibitors (), and the most frequent mutations were V82A/L/M (6/82; 7.3%) and I54L/M/V (5/82; 6%), L90M (4/82; 4.8%), and M46L/I (4/82; 4.8%) (Figure 1). ATZ/r was the most ARV used at the time of the genotyping test (23.1%, ) (Table 1). Concerning the resistance mutations, the most prevalent inter-ARV class combinations were M184/K103 and M184/K103/P225 (12.1% each) and M184/K70/K103 (8.5%). Multidrug resistance mutations were verified in 72% to NNRTI/NRTI, 5% NRTI/PI, and 3% NNRTI/NRTI/PI.

High resistance level to NRTI class was verified in 85.3% (), being most of them 82.9% () among those under 3TC use, and the M184I/V mutation was presented in all of them. However, only 6 from those 70 individuals that were taking TDF (8.5%) had high resistance level against TDF, despite its high use (Table 1), in the same way as high resistance level was detected just in 6 from those 46 individuals that were taking PI (13.0%). TDF remained acceptable as for being used for most of the individuals studied (91.5%), as well as the DRV to all individuals (Table 1). In regard to NNRTI, high ARV resistance level was detected in approximately 90.0% of them, representing 64 of those 71 individuals who were taking EFV, and mutations K103N/S (71.8) and P225H (21.1%) were the most prevalent among them. Among those individuals under use of NRTIs, crossover DRM confers high resistance level to abacavir in 47.5% of them, as well as for the 71 individuals under use of NNRTIs, and high resistance of 66.4% was verified to nevirapine, followed to 39.4% to rilpavirine, 19.7% to doravirine, and only 7.7% to etravirine (data not shown).

While analyzing the resistance profile of individuals with BCAR and BPAN variants, the heat map showed that the most frequent mutations in subtype variants were M184I/V (NRTI) and K103N/S (NNRTI), and these also presented similar prevalence in both groups (Figure 2(a)). Despite the low frequencies of TAMs in the general group below 20% (Figure 1), when we compared them among and , it was possible to observe that the TAM-2 mutational pattern was more present () in the variants (Figure 2(b)) and TAM-1 in the variants ().

(a)

(b)

Figure 2: Comparison of DRMs among subtype variants: () and (). (a) Heat map of the different classes of antiretroviral drugs: protease inhibitors (PI), nucleoside reverse transcriptase inhibitors (NRTI) and nonnucleoside reverse transcriptase inhibitors (NNRTI), and the frequency of drug-resistant mutations. (b) Profile of mutation pathways associated with thymidine—TAMs (TAM-1 () and () and TAM-2 () and () between subtype variants). Statistical analysis using the Mann–Whitney test: TAM-1 () TAM-2 ().

4. Discussion

Even though the northern region has presented one of the worst HIV epidemiological scenarios in Brazil in the last few years, there has been scarce surveillance of epidemiological and molecular resistance. In the present study regarding HIV-1 individuals failing cART individuals, the epidemiological data, genetic HIV-1 diversity,, and frequency of DRMs and their ARV used were analyzed. We observed both high prevalence of individuals with virologic failure without DRMs and also elevated frequency of DRMs, which can lead to an increase of TDRM (transmitted DRMs) and in HIV transmission in the region, as has been recently observed [28].

The general male-to-female ratio (M/F) detected was 2.5, in agreement with a 2.3 rate observed in 2018 in Brazil on a national level [18], and in 2011 in the state of Amazonas during the 2001-2012 HIV/AIDS epidemic [29]. Our results show a higher prevalence of heterosexual individuals. This category represents 71% of the study males, in contrast with the 31.4% observed nationally from 2007 to June 2019 in males over 13 years [18] .

Our results show a high prevalence of HIV-1 subtype and a very low presence of subtypes F1 and , as previously reported in other studies in the state of Amazonas [21, 22, 30, 31]. An increased number of HIV BF1 recombinant forms (13%) were verified when compared to other studies based on the pol region (4.3%-10%), mainly because of the inclusion of the env region in the analysis. Caribbean variants of subtype pol viruses were detected in the high rate of 15.7% when compared to other Brazilian regions, which goes in agreement with the 22% previously report at the same Amazon region [32]. A low prevalence (5.1%) of the variant of HIV-1 subtype (env) was observed [33], in contrast with 20-60% found in studies from other regions in the country [34–37].

The majority of the patients in our study had T CD4 counts below 200 cells/mm³ even when under a longer period of cART. This could be the result of a long virologic failure period until HIV genotyping and ARV replacement and also due to a late arrival of the individuals to care and consequently late HIV diagnosis. This may be demonstrated by the low CD4 count when analyzing the immune status of patients presenting for the first time, where half of patients showed CD4+ T cell counts below 350/mm³ and one-third below 200/mm³ at the initial HIV diagnosis [38]. These results also indicate a significant advance in the progression of the infection to AIDS in these individuals.

Concerning the viral tropism based in the sequences of the env region, we found that about 33% of the viruses presented X4 tropism; when considering only HIV-1 subtype env (), the X4 virus was verified in 38% of the individuals, which is very similar to previous reports of 40% of subtype virus switching during disease progression [39, 40]. On the other hand, 66.7% of the analyzed samples in our study still had R5 viral tropism and could be eligible for rescue therapies involving the use of maraviroc, as described by Bon et al. [41].

The studied individuals show high prevalence of DRMs associated with to NRTI (91%) and to NNRTI (84%). Other Brazilian studies have documented an increase in NNRTI DRMs over time, as well as the decrease in the rate of DRMs for NRTI [38, 42]. The increase of NNRTI resistance is mainly correlated to the presence of the K103N mutation, which was detected in 71.8% of individuals that use NNRTI, as also previously described [24]. The prevalence of K103N/S in our study (63.4%) was higher than the mean of 40% detected for the five largest Brazilian cities (São Paulo, Rio de Janeiro, Espirito Santo-Southeast region and Bahia and Ceará-Northeast region) and also in Belém, another city in the northern region [24, 38]. An elevated frequency of M184V (75%) was also detected herein, as well as in the study conducted in Belém [24]; this has ranged from 54.4% in São Paulo to 81.2% in Bahia [38]. Despite the low impact of the M184V alone to NRTI resistance, when taken together to other NRTI DRM, an increased level of resistance could be verified as shown in this study and in others [43–45].

A higher prevalence of TAM 1 was verified herein, as well as a lower prevalence of TAM-2, similar to those reported by Lopes et al. in the northern region [24]. It is important to point out that some resistance mutations observed in our and in other cross-sectional studies conducted with individuals failing cART could result from the transmitted virus since, up until now, genotyping prior to the HIV treatment is required only in a few cases, such as for pregnant women in Brazil, for example.

We observed high resistance to 3TC and EFV, and we have also demonstrated that the NRTIs (TDF and AZT) still presented high viability of use for most individuals. Due to the presence of mutations, such as L100I, K101P, Y181C, M230L, observed in of individuals failing in the use of EFV, the viability of using ETR in a rescue scheme is reduced for few individuals (data not shown) [46].

The high resistance levels found for these ARVs () reflect their high use and the genetic barrier to resistance of these ARVs [47]. Although some studies suggest that the viability of using first-line ARV drugs may be related to the combination of drugs with a low and/or high genetic barrier, consequently, the choice of drugs with a low barrier may favor the rapid selection of DRMs causing therapeutic failure [10]. The studied group has little experience with PIs, taken together to the fact that boosted PIs have a high genetic barrier, which could explain the low prevalence of major PI mutations observed. Our results demonstrated a high prevalence of NRTI/NNRTI dual class ARV NNRTI over the years in several studies, and resistance to NNRTI and NRTI is the most common forms of TDR [48–50] Due to the increased prevalence of resistance to NRTI and NNRTI-based regimens, ARVs focusing new targets on the HIV replication cycle were developed and implemented at the therapy. On a clinical trial comparing once-daily dolutegravir versus twice-daily raltegravir in antiretroviral-experienced individuals, the superior virological effect was verified in the dolutegravir group [51]. Dolutegravir has a high genetic barrier to resistance, which can be effective for HIV, both for first-line and rescue schemes. The study by Cahn et al. subsidized the World Health Organization to recommend the use of dolutegravir in the first therapeutic regimen since 2016, and it was implemented in Brazil in 2017. After the genotyping performed at the present study, various of the individuals started to use this ARV.

These factors can contribute to the DRM potentiating therapeutic failure and increase HIV transmission and TDRM prevalence in the region. Thus, the clinical and constant viral load monitoring of PLWHA is essential in order to reduce the risk of DRM accumulation and thus rapidly detect therapeutic failure, in order to preserve susceptibility to subsequent therapeutic schemes [15, 52]. However, protease inhibitors remain a good rescue therapy with less than 13% of the associated resistance verified, despite their use in 50% of individuals.

No significant epidemiological differences were seen among individuals infected by subtype variants, as previously verified among specific viral subtypes [53, 54]. Our results showed that the TAM-2 pathway was more prevalent in the variants, while TAM-1 was greater in , as shown previously in the northern region in blood donors [28]. It is important to point out that the TAM pathway verified could result from the specific use of thymidine analogs used during the treatment of these individuals. However, specific TAMs pathways were verified in some viral subtypes [55, 56], and these pathways are potentially related to the difficulty in achieving viral suppression [57]. However, further studies are needed to corroborate the specific characteristics concerning HIV-1 subtype variants ( and ) detected herein.

5. Conclusion

This is the most comprehensive analysis of the molecular epidemiology of HIV-1 and acquired DRMs in Manaus. Herein, we verified a high prevalence of the HIV-1 subtype , a significant difference in the TAM pathway in relation to variants ( and ), as well as an increase in recombinant BF1. We also verified that, despite the high resistance to 3TC and EFV, the low crossresistance of ARV of the same classes may favor rapid replication control after the change of viable ARV. Therefore, HIV-1 surveillance at local sites is essential for understanding the different obstacles for treatment adherence and rapid detection of therapeutic failure. This information will help the National Health Service redirect its strategies to achieve greater HIV/AIDS control.

Data Availability

All nucleotide sequences are available from the GenBank database (https://www.ncbi.nlm.nih.gov/genbank/) under accession numbers POL (MW545333-MW545424) and ENV (MW582429-MW582506).

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

MLG conceived and designed the experiments. YOC, RSP, DRLB, and MLG carried out the experiments. YOC, RSP, MLG, and FPR analyzed the data. MLG and PAN contributed reagents/materials/analysis tools. FPR and MVGL conducted patient recruitment and follow-up. YOC and MLG wrote the first draft.

All authors revised and approved the final manuscript.

Acknowledgments

The study is supported In collaboration with the Tropical Medicine Foundation Heitor Vieira Dourado and the Carlos Borborema Research Institute in Manaus, Amazonas, Brazil and the financial support from the Laboratory of AIDS and Molecular Immunology, Oswaldo Cruz Institute, FIOCRUZ, Rio de Janeiro, Rio de Janeiro, Brazil, the authors are grateful to the Genomic Platform DNA Sequencing–RPT01A and RPT01H (Rede de Plataformas Tecnológicas Fiocruz). “This study was partially supported by the Coordination for the Improvement of Higher Education Personnel Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-CAPES)-Finance Code 001.” MLG has a CNPq fellowship number 305919/2018-0.

References

A. E. H. V. Brasil, Ministério da Saúde, Secretaria de Vigilância em Saúde, and Departamento de DST, Cuidado integral às pessoas que vivem com HIV pela Atenção Básica Manual para a equipe multiprofissional, 2015.

A. E. H. V. Brasil, Ministério da Saúde, Secretaria de Vigilância em Saúde, and Departamento de DST, Agenda Estratégica para ampliação do acesso e cuidado integral das populações-chave em HIV, hepatites virais e outras infecções sexualmente transmissíveis, 2018.

C. J. Hoffmann, M. Schomaker, M. P. Fox et al., “CD4 count slope and mortality in HIV-infected patients on antiretroviral therapy,” JAIDS Journal of Acquired Immune Deficiency Syndromes, vol. 63, no. 1, pp. 34–41, 2013.

M. S. Saag, C. A. Benson, R. T. Gandhi et al., “Antiretroviral drugs for treatment and prevention of HIV infection in adults,” JAMA, vol. 320, no. 4, pp. 379–396, 2018.

J. Fokam, S. M. Sosso, B. Yagai et al., “Viral suppression in adults, adolescents and children receiving antiretroviral therapy in Cameroon: adolescents at high risk of virological failure in the era of ‘test and treat,” AIDS Research and Therapy, vol. 16, no. 1, pp. 10–17, 2019.

M. S. Cohen, Y. Q. Chen, M. McCauley et al., “Antiretroviral therapy for the prevention of HIV-1 transmission,” New England Journal of Medicine, vol. 375, no. 9, pp. 830–839, 2016.

J. B. Nachega, V. C. Marconi, G. U. van Zyl et al., “HIV treatment adherence, drug resistance, virologic failure: evolving concepts,” Infectious Disorders-Drug Targets, vol. 11, no. 2, pp. 167–174, 2011.

Y. Yu, D. Luo, X. Chen, Z. Huang, M. Wang, and S. Xiao, “Medication adherence to antiretroviral therapy among newly treated people living with HIV,” BMC Public Health, vol. 18, no. 1, 2018.

P. Ü. Calvetti, G. R. M. Giovelli, G. J. C. Gauer, and J. F. D. de Moraes, “Psychosocial factors associated with adherence to treatment and quality of life in people living with HIV/AIDS in Brazil,” Jornal Brasileiro de Psiquiatria, vol. 63, no. 1, pp. 8–15, 2014.

C. Katusiime, P. Ocama, and A. Kambugu, “Basis of selection of first and second line highly active antiretroviral therapy for hiv/aids on genetic barrier to resistance: a literature review,” African Health Sciences, vol. 14, no. 3, pp. 679–681, 2014.

A. D. L. G. Corado, G. Bello, R. A. C. Leão, F. Granja, and F. G. Naveca, “HIV-1 genetic diversity and antiretroviral drug resistance among individuals from Roraima state, northern Brazil,” PLoS One, vol. 12, no. 3, p. e0173894, 2017.

S. Sanguansittianant, N. Nooroon, P. Phaengchomduan, and P. Ammaranond, “Trends in prevalence of HIV-1 drug resistance in Thailand 2009-2010,” Journal of Clinical Laboratory Analysis., vol. 27, no. 5, pp. 346–353, 2013.

R. K. Gupta, A. Hill, A. W. Sawyer et al., “Virological monitoring and resistance to first-line highly active antiretroviral therapy in adults infected with HIV-1 treated under WHO guidelines: a systematic review and meta-analysis,” The Lancet Infectious Diseases., vol. 9, no. 7, pp. 409–417, 2009.

A. E. H. V. Brasil, Ministério da Saúde, Secretaria de Vigilância em Saúde, and Departamento de DST, Recomendações para terapia anti-retroviral em adultos infectados pelo HIV, 2008.

K. C. E. Sigaloff, T. Ramatsebe, R. Viana, T. F. R. de Wit, C. L. Wallis, and W. S. Stevens, “Accumulation of HIV drug resistance mutations in patients failing first-line antiretroviral treatment in South Africa,” AIDS Research and Human Retroviruses., vol. 28, no. 2, pp. 171–175, 2012.

D. Thirunavukarasu, V. Udhaya, H. S. Iqbal, and T. Umaarasu, “Patterns of HIV-1 drug-resistance mutations among patients failing first-line antiretroviral treatment in South India,” Journal of the International Association of Providers of AIDS Care, vol. 15, no. 3, pp. 261–268, 2013.

B.M da S.S de em S.D. de D de C.C e I S T (DCCI), Brasil, Ministério da Saúde, Secretaria de Vigilância em Saúde, and Departamento de DST, Manual Técnico para Avaliação de Exames de Genotipagem do HIV, 2019.

A. E. H. V. Brasil, Ministério da Saúde, Secretaria de Vigilância em Saúde, and Departamento de DST, Boletim Epidemiológico HIV/Aids |2019, 2019.

T. Gräf and A. R. Pinto, “The increasing prevalence of HIV-1 subtype C in southern Brazil and its dispersion through the continent,” Virology, vol. 435, no. 1, pp. 170–178, 2013.

B. C. Carvalho, L. P. V. Cardoso, S. Damasceno, and M. M. D. A. Stefani, “Moderate prevalence of transmitted drug resistance and interiorization of HIV type 1 subtype C in the inland north state of Tocantins, Brazil,” AIDS Research and Human Retroviruses., vol. 27, no. 10, pp. 1081–1087, 2011.

L. K. H. Cunha, S. Kashima, M. F. C. Amarante et al., “Distribution of human immunodeficiency virus type 1 subtypes in the state of Amazonas, Brazil, and subtype C identification,” Brazilian Journal of Medical and Biological Research., vol. 45, no. 2, pp. 104–112, 2012.

C. M. Da Costa, C. M. Costa De Oliveira, Y. F. Chehuan De Melo, E. Delatorre, G. Bello, and J. C. Couto-Fernandez, “High HIV-1 genetic diversity in patients from northern Brazil,” AIDS Research and Human Retroviruses., vol. 32, no. 9, pp. 918–922, 2016.

L. Dos Anjos Silva, F. Divino, M. O. Da Silva Rêgo et al., “HIV-1 genetic diversity and transmitted drug resistance in antiretroviral treatment-naive individuals from Amapá state, northern Brazil,” AIDS Research and Human Retroviruses., vol. 32, no. 4, pp. 373–376, 2016.

C. A. F. Lopes, M. A. Soares, D. R. Falci, and E. Sprinz, “The evolving genotypic profile of HIV-1 mutations related to antiretroviral treatment in the north region of Brazil,” BioMed Research International, vol. 2015, 11 pages, 2015.

F. Divino, A. L. G. De Corado, F. G. Naveca, M. M. A. Stefani, and G. Bello, “High prevalence and onward transmission of non-pandemic HIV-1 subtype B clades in northern and northeastern brazilian regions,” PLoS One, vol. 11, no. 9, p. e0162112, 2016.

E. Delwart, E. Shpaer, J. Louwagie et al., “Genetic relationships determined by a DNA heteroduplex mobility assay: analysis of HIV-1 env genes,” Science, vol. 262, no. 5137, pp. 1257–1261, 1993.

E. Delatorre, C. Silva-De-Jesus, J. C. Couto-Fernandez, J. H. Pilotto, and M. G. Morgado, “High HIV-1 diversity and prevalence of transmitted drug resistance among antiretroviral-naive HIV-infected pregnant women from Rio de Janeiro, Brazil,” AIDS Research and Human Retroviruses., vol. 33, no. 1, pp. 68–73, 2017.

M. A. Esashika Crispim, M. N. da Guarda Reis, N. Fraiji, G. Bello, and M. M. A. Stefani, “Detection of human immunodeficiency virus type 1 phylogenetic clusters with multidrug resistance mutations among 2011 to 2017 blood donors from the highly endemic northern Brazilian Amazon,” Transfusion, vol. 59, no. 8, pp. 2593–2601, 2019.

R. D. S. M. De Oliveira, A. S. Benzaken, V. Saraceni, and M. Sabidó, “Hiv/aids epidemic in the state of amazonas: Characteristics and trends from 2001 to 2012,” Revista da Sociedade Brasileira de Medicina Tropical, vol. 48, suppl 1, pp. 70–78, 2015.

M. A. E. Crispim, M. N. D. G. Reis, C. Abrahim et al., “Homogenous HIV-1 subtype B from the Brazilian Amazon with infrequent diverse BF1 recombinants, subtypes F1 and C among blood donors,” PLOS ONE, vol. 14, no. 9, p. e0221151, 2019.

S. D. De Andrade, M. Sabidó, W. M. Monteiro, A. S. Benzaken, and A. Tanuri, “Drug resistance in antiretroviral-naive children newly diagnosed with HIV-1 in Manaus, Amazonas,” Journal of Antimicrobial Chemotherapy, vol. 72, no. 6, pp. 1774–1783, 2017.

A. C. P. Vicente, K. Otsuki, N. B. Silva et al., “The HIV epidemic in the Amazon Basin is driven by prototypic and recombinant HIV-1 subtypes B and F,” Journal of Acquired Immune Deficiency Syndromes, vol. 23, no. 4, pp. 327–331, 2000.

M. G. Morgado, M. L. Guimarães, C. B. G. Gripp et al., “Molecular epidemiology of HIV-1 in Brazil: high prevalence of HIV-1 subtype B and identification of an HIV-1 subtype D infection in the city of Rio de Janeiro, Brazil,” Journal of Acquired Immune Deficiency Syndromes and Human Retrovirology, vol. 18, no. 5, pp. 488–494, 1998.

V. F. Pimentel, M. G. Morgado, G. Bello et al., “Temporal trends and molecular epidemiology of HIV type 1 infection in Rio de Janeiro, Brazil,” AIDS Research and Human Retroviruses., vol. 29, no. 12, pp. 1553–1561, 2013.

K. E. Potts, M. L. Kalish, T. Lott et al., “Genetic heterogeneity of the V3 region of the HIV-1 envelope glycoprotein in Brazil,” AIDS, vol. 7, no. 9, pp. 1191–1198, 1993.

M. G. MORGADO, E. C. SABINO, E. G. SHPAER et al., “V3 region polymorphisms in HIV-1 from Brazil: prevalence of subtype B strains divergent from North American/European prototype and detection of subtype F,” AIDS Research and Human Retroviruses, vol. 10, no. 5, pp. 569–576, 1994.

D. M. Junqueira, R. M. de Medeiros, T. C. N. Ferreira Leite et al., “Detection of the B″-GWGR variant in the southernmost region of Brazil: unveiling the complexity of the human immunodeficiency virus-1 subtype B epidemic,” Memorias do Instituto Oswaldo Cruz, vol. 108, no. 6, pp. 735–740, 2013.

R. S. Diaz, L. A. Inocêncio, M. C. A. Sucupira et al., “The virological and immunological characteristics of the HIV-1-infected population in Brazil: from initial diagnosis to impact of antiretroviral use,” PLoS One, vol. 10, no. 10, p. e0139677, 2015.

M. Mild, A. Kvist, J. Esbjörnsson, I. Karlsson, E. M. Fenyö, and P. Medstrand Patrik, “Differences in molecular evolution between switch (R5 to R5X4/X4-tropic) and non-switch (R5-tropic only) HIV-1 populations during infection,” Infection, Genetics and Evolution, vol. 10, no. 3, pp. 356–364, 2010.

C. Tscherning, A. Alaeus, R. Fredriksson et al., “Differences in chemokine coreceptor usage between genetic subtypes of HIV-1,” Virology, vol. 241, no. 2, pp. 181–188, 1998.

I. Bon, A. Clò, M. Borderi et al., “Prevalence of R5 strains in multi-treated HIV subjects and impact of new regimens including maraviroc in a selected group of patients with CCR5-tropic HIV-1 infection,” International Journal of Infectious Diseases, vol. 17, no. 10, pp. e875–e882, 2013.

H. Duani, A. W. Aleixo, and U. Tupinambás, “Trends and predictors of HIV-1 acquired drug resistance in Minas Gerais, Brazil: 2002–2012,” Brazilian Journal of Infectious Diseases, vol. 21, no. 2, pp. 148–154, 2017.

P. R. Harrigan, C. Stone, P. Griffin et al., “Resistance profile of the human immunodeficiency virus type 1 reverse transcriptase inhibitor abacavir (1592U89) after monotherapy and combination therapy,” The Journal of Infectious Diseases, vol. 181, no. 3, pp. 912–920, 2000.

L. Bacheler, S. Jeffrey, G. Hanna et al., “Genotypic correlates of phenotypic resistance to efavirenz in virus isolates from patients failing nonnucleoside reverse transcriptase inhibitor therapy,” Journal of Virology, vol. 75, no. 11, pp. 4999–5008, 2001.

G. L. Melikian, S. Y. Rhee, V. Varghese et al., “Non-nucleoside reverse transcriptase inhibitor (NNRTI) cross-resistance: implications for preclinical evaluation of novel NNRTIs and clinical genotypic resistance testing,” Journal of Antimicrobial Chemotherapy, vol. 69, no. 1, pp. 12–20, 2014.

J. V. Madruga, P. Cahn, B. Grinsztejn et al., “Efficacy and safety of TMC125 (etravirine) in treatment-experienced HIV-1-infected patients in DUET-1: 24-week results from a randomised, double-blind, placebo-controlled trial,” Lancet, vol. 370, no. 9581, pp. 29–38, 2007.

M. W. Tang and R. W. Shafer, “HIV-1 Antiretroviral Resistance,” Drugs, vol. 72, no. 9, pp. e1–e25, 2012.

S. Y. Rhee, J. L. Blanco, M. R. Jordan et al., “Geographic and temporal trends in the molecular epidemiology and genetic mechanisms of transmitted HIV-1 drug resistance: an individual-patient- and sequence-level meta-analysis,” PLoS Medicine, vol. 12, no. 4, 2015.

R. K. Gupta, M. R. Jordan, B. J. Sultan et al., “Global trends in antiretroviral resistance in treatment-naive individuals with HIV after rollout of antiretroviral treatment in resource-limited settings: a global collaborative study and meta-regression analysis,” The Lancet., vol. 380, no. 9849, pp. 1250–1258, 2012.

A. M. Wensing, V. Calvez, F. Ceccherini-Silberstein et al., “2019 update of the drug resistance mutations in HIV-1,” Topics in antiviral medicine., vol. 27, no. 3, pp. 111–121, 2019.

P. Cahn, A. L. Pozniak, H. Mingrone et al., “Dolutegravir versus raltegravir in antiretroviral-experienced, integrase- inhibitor-naive adults with HIV: week 48 results from the randomised, double- blind, non-inferiority SAILING study,” The Lancet, vol. 382, no. 9893, pp. 700–708, 2013.

A. De Luca, Z. J. Sidumo, G. Zanelli et al., “Accumulation of HIV-1 drug resistance in patients on a standard thymidine analogue-based first line antiretroviral therapy after virological failure: implications for the activity of next-line regimens from a longitudinal study in Mozambique,” BMC Infectious Diseases, vol. 17, no. 1, 2017.

C. M. Venner, I. Nankya, F. Kyeyune et al., “Infecting HIV-1 subtype predicts disease progression in women of sub-Saharan Africa,” eBioMedicine, vol. 13, pp. 305–314, 2016.

N. Kiwanuka, O. Laeyendecker, M. Robb et al., “Effect of human immunodeficiency virus type 1 (HIV-1) subtype on disease progression in persons from Rakai, Uganda, with incident HIV-1 infection,” The Journal of Infectious Diseases., vol. 197, no. 5, pp. 707–713, 2008.

V. Novitsky, C. W. Wester, V. DeGruttola et al., “The reverse transcriptase 67N 70R 215Y genotype is the predominant TAM pathway associated with virologic failure among HIV type 1C-infected adults treated with ZDV/ddI-containing HAART in southern Africa,” AIDS Research and Human Retroviruses., vol. 23, no. 7, pp. 868–878, 2007.

A. G. Marcelin, C. Delaugerre, M. Wirden et al., “Thymidine analogue reverse transcriptase inhibitors resistance mutations profiles and association to other nucleoside reverse transcriptase inhibitors resistance mutations observed in the context of virological failure,” Journal of Medical Virology, vol. 72, no. 1, pp. 162–165, 2004.

G. J. Hanna, V. A. Johnson, D. R. Kuritzkes et al., “Patterns of resistance mutations selected by treatment of human immunodeficiency virus type 1 infection with zidovudine, didanosine, and nevirapine,” The Journal of Infectious Diseases, vol. 181, no. 3, pp. 904–911, 2000.