Sensitive detection of HIV-1 resistance to Zidovudine and impact on treatment outcomes in low- to middle-income countries
Infectious Diseases of Poverty volume 6, Article number: 163 (2017)
Thymidine analogs, namely AZT (Zidovudine or Retrovir™) and d4T (Stavudine or Zerit™) are antiretroviral drugs still employed in over 75% of first line combination antiretroviral therapy (cART) in Kampala, Uganda despite aversion to prescribing these drugs for cART in high income countries due in part to adverse events. For this study, we explored how the continued use of these thymidine analogs in cART could impact emergence of drug resistance and impact on future treatment success in Uganda, a low-income country.
We examined the drug resistance genotypes by Sanger sequencing of 262 HIV-infected patients failing a first line combined antiretroviral treatment containing either AZT or d4T, which represents approximately 5% of the patients at the Joint Clinical Research Center receiving a AZT or d4T containing treatment. Next generation sequencing (DEEPGEN™HIV) and multiplex oligonucleotide ligation assays (AfriPOLA) were then performed on a subset of patient samples to detect low frequency drug resistant mutations. CD4 cell counts, viral RNA loads, and treatment changes were analyzed in a cohort of treatment success and failures.
Over 80% of patients failing first line AZT/d4T-containing cART had predicted drug resistance to 3TC (Lamivudine) and non-nucleoside RT inhibitors (NNRTIs) in the treatment regimen but only 45% had resistance AZT/d4T associated resistance mutations (TAMs). TAMs were however detected at low frequency within the patients HIV quasispecies (1–20%) in 21 of 34 individuals who were failing first-line AZT-containing cART and lacked TAMs by Sanger. Due to lack of TAMs by Sanger, AZT was typically maintained in second-line therapies and these patients had a low frequency of subsequent virologic success.
Our findings suggest that continued use of AZT and d4T in first-line treatment in low-to-middle income countries may lead to misdiagnosis of HIV-1 drug resistance and possibly enhance a succession of second- and third-line treatment failures.
Please see Additional file 1 for translations of the abstract into the five official working languages of the United Nations.
Uganda has an estimated 1.6 million people living with human immunodeficiency virus type 1 (HIV-1) with 140 000 new diagnoses and 63 000 AIDS-related deaths per year, which has significant socioeconomic consequences on the nation’s population . The advent of combination antiretroviral therapy (cART) in Uganda resulted in a marked reduction in morbidity and mortality among HIV-infected individuals [2, 3], associated with stabilization and potential recovery of immune function , and even helping with prevention of HIV-1 transmission . Although 28 antiretroviral drugs from six different classes are currently approved for treatment , patients in Uganda and other low-to-middle income countires (LIMCs) rarely have access to once a day, all-in-one drug combination pills of the most effective combined antiretroviral therapies (cARTs), e.g., Triumeq, Stribild, Complera or Atripla. Rather, antiretroviral drugs are commonly administered as multiple pill cART, or in co-formulations of generic cART such as TRIOMUNE (Nevirapine[NVP] + Stavudine[d4T] + Lamivudine[3TC]) or DUOVIR-E or -N (Efavirenz[EFV] or NVP + Zidovudine[AZT] + 3TC). Although treatment with d4T in LIMCs is less common in the past 5–7 years , AZT is still heavily utilized in first-line treatment regimens, even though the use of all thymidine analogues have largely been abandoned in high-income-countries (HICs) due to side effects such as anemia, neutropenia, and neuropathy (the latter commonly related to d4T) [7,8,9].
The HIV drug resistance surveillance program conducted by the World Health Organization (WHO) in 2012 concluded that 72% of virologic failures during antiretroviral therapy (ART) in LIMCs were directly associated with HIV-1 drug resistance . In high income countries, HIV-1 with primary drug resistant mutations are observed at greater than 80% of the patients when initially failing cART which increases to greater 90% if the patient remains on the failing therapy for another 6 months  in most settings around the world, HIV-1 drug resistance has been monitored in Uganda using standard Sanger sequencing-based HIV-1 genotyping assays [12,13,14,15]. During virologic failure of an AZT or d4T containing treatment regimen, there is commonly a timed order in the appearance of specific thymidine analog resistance mutations (TAMs), i.e. K70R often appears early in the patient’s HIV reverse transcriptase gene followed by L210 W, T215F/Y, M41 L and D67N [16, 17]. All of these TAMs have been readily detected in HIV-infected Ugandans failing AZT or d4T containing cART [12, 13]. This step-by-step appearance of TAMs appears to be directly proportional to the level of drug resistance and/or the viral replicative fitness cost associated with the accumulation of these mutations [18,19,20]. Due to the high fitness cost of these mutations, HIV-1 within a patient will revert back to wild type HIV-1 in the absence of treatment [21, 22], but these drug resistance mutations may still be observed at a low frequency in the patient’s viral quasispecies. Recent findings suggest low-frequency HIV-1 drug resistance mutations are readily observed in Ugandan patients failing treatment , i.e., below the limit of Sanger sequence detection of 20 to 30% of the viral population [24,25,26]. The clinical significance of these low-abundance HIV-1 drug resistant variants remains controversial [23, 27,28,29,30,31]; however, treatment with specific antiretroviral drugs (e.g., AZT) is often associated with the appearance of the respective drug resistance mutations at low frequency (e.g., TAMs) especially if subsequent treatment failure ensues [27, 32,33,34,35,36].
In this study, we identified 50 HIV-infected individuals receiving first-line cART with or without AZT or d4T and experiencing virologic failure while harboring viruses with dominant FTC/3TC resistance mutations (e.g., M184I/V) or NNRTI resistance mutations, but lacking any TAMs based on Sanger sequencing. As expected for Uganda, the majority of these patients were retained on AZT or d4T and received a treatment boost by adding a protease inhibitor while removing or retaining the NNRTI. We used both a deep sequencing-based HIV-1 genotyping assay (DEEPGEN™HIV)  and an oligonucleotide ligation assay (AfriPOLA ) to quantify both dominant and low frequency HIV-1 drug resistant variants in these HIV-infected Ugandan patients. To better understand the impact of minority TAM mutations on patient outcome, CD4+ T-cell counts and plasma HIV RNA loads were examined prior to and following the initial Sanger sequencing drug resistance test.
We originally identified 262 HIV-infected individuals failing first-line AZT- or d4T-containing cART at the Joint Clinical Research Center (JCRC; Kampala, Uganda). Plasma-derived RNA specimens from a subset of these patients (n = 50) were selected to be analyzed using the AfriPOLA  and DEEPGEN™HIV  assays. Clinical and virological data was obtained from the patient care database at the JCRC under IRB approval (EM10-07) for HIV-1 drug resistance testing. At the JCRC, virologic failure during treatment is defined as plasma HIV-1 RNA load above 2000 copies/ml and/or a loss of CD4+ T-cells to counts below 250 cells/ml. Sanger sequencing-based HIV-1 genotyping assays had been performed on all patient samples at the time of treatment failure (Table 1).
HIV-1 genotyping based on Sanger sequencing
RT-PCR products, corresponding to the HIV-1 reverse transcriptase coding regions (pol-PR-RT) of HIV-1, were sequenced in the WHO and NIH VQA-certified Case Western Reserve University/Center for AIDS Research Uganda Laboratory at the JCRC as part of routine HIV-1 drug resistance testing . HIV-1 sequences were interpreted, and drug resistance profiles generated, based on the HIVdb Program Genotypic Resistance Interpretation Algorithm from the Stanford University HIV Drug Resistance Database (http://hivdb.stanford.edu). Phylogenetic analysis was used to predict HIV-1 subtype as previously described .
Oligonucleotide ligation assay for HIV-1 genotyping of the pol-RT coding region (AfriPOLA)
The POLA assay has been previously defined, developed, optimized, and patented (US Patent# 9487839) to measure HIV-1 drug resistance (even at frequencies as low as 1%) in patients receiving cART . For AfriPOLA, the oligonucleotides were of the same design as described , but were altered to allow for coverage of subtype A, B, C, and D HIV-1 strains and detect the following NRTI resistant mutations: K65R, D67N, K70R, L74 V, Y115F, M184 V, L201 W, T215Y and K219Q. Ligase discrimination reactions (LDRs) using the oligonucleotide sets were performed on RT-PCR products corresponding to the HIV-1 RT as described . Following the LDR, only oligonucleotides that annealed and ligated opposite to a wild type or drug resistant mutation site on the PCR products were distinguished by subsequent hybridization to Luminex MTAG beads, via an anti-TAG bead specific sequence on 5′ end of the ligated produce and a biotin tag on 3’end. Populated MTAG beads were then measured using a Magpix instrument (Luminex) following addition of Streptavidin-R-phycoerythrin conjugate. Samples were run in triplicate and assay cutoff values were determined by calculating the 95% confidence interval (CI) for the signal obtained from the solvent template alone.
HIV-1 genotyping based on deep sequencing of the pol-RT coding region
An RT-PCR product corresponding to the pol-RT coding region of HIV-1 was sequenced using the HIV-1 genotyping assay DEEPGEN™HIV as described . The clinical DEEPGEN™HIV assay was performed in the CLIA/CAP-accredited University Hospitals Translational Laboratory under a good laboratory practice framework. Sequence reads were mapped and aligned against sample-specific reference sequences constructed for the pol-RT genomic region using the DEEPGEN Software Tool Suite . In this study, minority variants were defined as amino acid substitutions detected at ≥1% (based on the intrinsic error rate of the sequencing system ) and <20% of the virus population (based on the lower limit of detection by Sanger population sequencing [24,25,26, 39, 40]).
Results are expressed as median values, standard deviations, and confidence intervals. Two population proportional tests (one tailed) were used to compare drug regimen and percentages of drug resistance. ANOVA was used to compare longitudinal measurements of plasma viral loads and CD4+ T-cell counts. Pearson correlation coefficients were calculated for all comparison studies. All differences with a P value of <0.05 were considered statistically significant. All statistical analyses were performed using GraphPad Prism v.6.0b (GraphPad Software, La Jolla, CA) unless otherwise specified. pol-PR/RT nucleotide sequences obtained by deep sequencing in this study have been submitted to the Los Alamos National Laboratory HIV-DB Next Generation Sequence Archive (http://www.hiv.lanl.gov/content/sequence/HIV/NextGenArchive /Gibson2017).
Thymidine analog resistance in patients receiving cART with single d4T or AZT pills or as a single pill formulation
At the JCRC, over 5000 patients are on or have received AZT or d4T as part of first-line cART. In first-line cART, AZT is most frequently prescribed, more recently as Combivir™ or the Cipla Inc. generic, Durovirv (i.e., both in combination with 3TC). d4T is commonly prescribed QD as part of the Cipla generic single pill Triomune™ combination with 3TC and NVP. For this study, we first explored emergence of TAMs in 262 patients failing cART where AZT or d4T was administered as generic or brand name pill(s) alone or co-formulated with other ARVs.
With d4T- and AZT-containing cART, drug resistance was observed in over 80% of virologic and/or immunologic treatment failures, defined as having 1+ primary DRMs conferring resistance to any of the three drugs in the first-line treatment regimen (Fig. 1). However, based on Sanger sequencing less than 50% of patients had primary TAMs in the RT coding region (Fig. 1). About 3% to 8% of HIV-1 drug resistance genotypes (n = 13) in all three cART groups described in Fig. 1 had a multi-NRTI drug resistance pattern (combinations of A62V, 69 insertion, V75I, F77 L, F116Y, or Q151M). Similar percentages of AZT resistance (~45%) were observed regardless of infecting subtype in these patients (A, C, and D). NNRTI or 3TC resistance was observed in over 75% of patients receiving AZT- or d4T-containing cART. TAMs in absence of any NRTI or NNRTI resistance mutations were observed in only two of these 262 (0.8%) patient samples. There were no statistically significant differences (based on ANOVA) in the frequencies of resistance mutations to NNRTI (range 77–81%), 3TC (range 72–79%), or thymidine analogs (range 41–50%) among patients failing different AZT-containing treatment regimens: (i) AZT + 3TC + NNRTI administered as separate pills (BID/QD/QD, respectively) [n = 140], (ii) Triomune (d4T + 3TC + NVP) (one pill combination QD) [n = 53], or (iii) AZT + 3TC administered as Combivir or Duovir BID + NNRTI (NVP or EFV QD) [n = 69].
Detection of low frequency TAMs in patients failing AZT-containing cART but with AZT-susceptible HIV-1 genotypes based on Sanger sequencing
From the cohort described above, we selected 35 HIV-infected Ugandan patients treated with AZT- or d4T-containing cART and experiencing virologic failure without dominant AZT or d4T resistance detected by Sanger sequencing-based HIV-1 genotyping (AZT/d4T group). As controls, we evaluated 15 patients failing cART for which AZT or d4T were not prescribed (control group). The distribution of HIV-1 subtypes among these patients was consistent with the region: subtype A (20/50, 40%), D (21/50, 42%), B (1/50, 2%), C (2/50, 4%) or CRF01_AE (3/50, 6%) (Table 1). Despite apparently lacking TAMs, 49 of 50 (98%) patients had resistance to 3TC and FTC based on the presence of M184 V or M184I mutations. Five individuals had additional drug resistance mutations, such as A62V, T69 N/A, K70E, and/or L74 V/I, which have been associated with multi-NRTI resistance .
AfriPOLA is an assay used to measure low frequency TAMs and other select drug resistance mutations (>1% in the patient’s quasispecies) in the HIV-1 pol gene of patients failing an AZT/d4T (with or without dominant TAMs detected by Sanger sequencing) (Fig. 2). Resistance mutations at both low (1% to 20%) or high frequencies (>20%) within patients were then used to predict maximum resistances to specific antiretroviral drugs (Fig. 3) with the HIVdb algorithm. Using AfriPOLA, we probed for 9 specific HIV-1 drug resistance mutations: K65R, D67N, K70R, L74 V, Y115F, T215Y, K219Q, L210 W, and M184 V in all 50 HIV-infected individuals. Detection of M41 L by AfriPOLA resulted in low to unreportable signal, most likely as a result of oligonucleotide binding constraints during ligase discrimination. Enhanced detection of TAMs by AfriPOLA resulted in 19 of 34 (56%) patients in the AZT/d4T group identified as harboring low-frequency AZT resistance, despite being originally reported as AZT susceptible based on Sanger sequencing (Fig. 2a). Using AfriPOLA, an average of two TAMs at low frequencies were detected per patient in this group. Similarly, TAMs were detected in the control group by AfriPOLA (7/15) assay. When comparing specific TAMs, we detected D67N in 11 patients by AfriPOLA not detected by Sanger sequencing, as well as K70R in 18 patients, L210 W in 5, T215Y in 8, and K219Q in 8 patients by AfriPOLA but not detected by Sanger (Fig. 2b).
In 94% (47/50) of patients, AfriPOLA identified at least one additional drug resistance-associated mutation that was not reported by Sanger sequencing. Percentage of patients infected with viruses carrying the K65R, L74 V, and Y115F mutations at any frequency (>1%) increased from 0 to 18%, 6 to 44%, and 6 to 40%, respectively (Fig. 2b). AfriPOLA detected mutations at low frequency and predicted resistance to ABC (K65R, L74 V, and Y115F) in all four patients receiving ABC and originally determined as ABC sensitive by Sanger sequencing and interpretive algorithms (Figs. 2a and 3). Type 2 TAMs (M41 L, L210 W, and T215Y) or K65R were detected only at low frequency by AfriPOLA in two patients (DR0279–08 and DR-0319-08, respectively) treated with a TDF-containing regimen. High level TDF resistance may have emerged in one TDF treated patient (DR0242-08) based on the low frequency detection of type 1 TAMs (M41 L, L210 W, and T215Y) by AfriPOLA (Fig. 2a and b). Another four patients treated with TDF-containing regimens did not harbor TDF resistance (dominant or low frequency) despite high level resistance to cytidine analogy (3TC/FTC) and to NNRIs (EFV/NVP).
To confirm the AfriPOLA results, we selected 15 samples of the 50 clinical specimens analyzed by Sanger sequencing and AfriPOLA to be tested by the deep sequencing-based DEEPGEN™HIV assay. Using DEEPGEN™HIV we confirmed the presence of dominant mutations (>20% frequency) identified by Sanger sequencing and AfriPOLA, including the highly prevalent M184 V mutation (Fig. 2c). To directly compare the AfriPOLA results with DEEPGEN a box plot was created comparing all amino acid sequencing frequencies to AfriPOLA positive (>1%) or negative results (<1%) (Fig. 2d). Due to some variability in oligonucleotide annealing efficiency to the divergent DNA templates of each patient, the absolute mutations frequency can vary between patient samples but the limit of mutation detection by AfriPOLA was set to two standard deviations above the assay cutoff . We obtained concordant results for 59 of mutations using DEEPGEN and AfriPOLA. In cases of undetectable mutation frequencies by AfriPOLA but detectable mutations frequencies by DEEPGEN (18 disconcordant results), the frequency detected by DEEPGEN was typically <10% which may have corresponded to a read by AfriPOLA just below the limit of detection. In summary, AfriPOLA and DEEPGEN™HIV predicted resistance to AZT and d4T in 27 of 35 (77%) patients failing an AZT-containing treatment regimen, while Sanger sequencing predicted resistance to AZT and d4T in only 1 of the 35 (2.9%) patient samples.
We compared the HIV-1 genotyping data obtained using Sanger sequencing, AfriPOLA, and DEEPGEN™HIV (Fig. 4). Concordance among the different HIV-1 genotyping assays varied depending on the RT site (amino acid position); however, results from AfriPOLA and DEEPGEN™HIV showed stronger correlation with each other than compared with the less sensitive Sanger sequencing method (Fig. 4a). DEEPGEN™HIV and AfriPOLA showed a strong significant correlation detecting the low frequency TAMs as well as the other drug resistance mutations (r = 0.88, P < 0.0001, Pearson coefficient correlation) while no correlation was observed between any of the high sensitive assays and Sanger sequencing (Fig. 4b). In most cases, non-concordant results between DEEPGEN™HIV and AfriPOLA were most common when the HIV-1 drug resistance mutations were detected at the lowest frequency detection (<5%) by DEEPGEN™HIV. Reduced sensitivity for the lowest frequency mutations by AfriPOLA relates to limited signal on the Luminex beads in comparison to background. For both AfriPOLA and DEEPGEN™HIV, we are monitoring over 1000 HIV templates when the viral RNA in plasma is >10 000 copies/ml. Low frequency mutations (<5%) may not be as accurate if viral RNA is below 10 000 copies/ml of plasma derived from the viral RNA in plasma (i.e. 15 of 50 samples; Table 1).
Patient outcome post Sanger drug resistance testing
Patients failing cART and with resistance to all three drugs including AZT (n = 157) would typically start a new second line treatment regimen, which likely explains for the significant reduction in plasma HIV-1 RNA loads and increase in CD4+ T cells 3 months post drug resistance testing (Fig. 5a). Again, the majority of these patients would initiate a PI-based treatment regimen following detection of resistance to a first line AZT/d4T-containing treatment regimen. In contrast, treatment failure in the absence of AZT resistance but with resistance to the NNRTI and 3TC in the regimen (n = 99) led to no statistically significant difference in the high viral load or low CD4 cell counts before and after Sanger-based drug resistance testing (Fig. 5b). In a majority of these cases, patients would remain on AZT or d4T while the physician at the JCRC might intensify the treatment regimen by adding new drugs. Based on data presented herein, we predict that these patients would have low level TAMs that might influence treatment outcomes. Again, our findings indicate that 94% (47/50) of patients failing AZT- (or d4T-) without TAMs detected by Sanger did harbor low frequency TAMs detected by AfriPOLA and/or DEEPGEN. Next, we specifically examined the patient subset receiving but failing an AZT/d4T-containing cART (n = 27 with available follow up data) where we tested for and found low frequency TAMs. Viral loads and CD4 T cell counts remained unchanged over the 90 days prior to and 90 days post drug resistance testing (Fig. 5c). However, there was trend for decreases in viral load and increase CD4 T cell counts following 90 days post testing which relates a change in treatment regimen over this first year (post 90 days). In the subset cohort presented in Fig. 5c where we now know that low frequency TAMs were present, the physicians received the Sanger drug resistance report indicating the absence of TAMs. In 21 of 27 patients AZT (or d4T) was maintained in the salvage treatment regimen. Thus, complete viral suppression may never be achieved and subsequent treatment failure may be associated with new drug resistance mutations which further limits antiretroviral drugs for third line treatment.
Zidovudine (AZT or Retrovir™) was the first FDA approved drug for treatment of HIV infections . Within a year of AZT approval and widespread use, AZT-resistant HIV was readily identified in almost all patients receiving AZT for longer than 6 months . Mutations in HIV-1 reverse transcriptase conferring AZT resistance have proven more complicated than resistance to any other nucleoside RT inhibitor in terms of (i) sheer number of primary mutations, (ii) complex pattern of these mutations, (iii) temporal emergence of these mutations during treatment failure, (iv) high genetic barrier and high relative fitness costs of AZT resistance, (v) initial incompatibility of AZT and 3TC resistance mutations, and (vi) unique mechanism(s) of AZT resistance by RT [16, 44,45,46,47,48,49,50,51]. All of these aspects of AZT resistance actually improve effectiveness of AZT in combination therapy and, until recently, AZT was still the favored NRTI in first line cART regimens in LMICs. Estimates suggest that the accumulative AZT treatment years (>25 million years) is greater than any other drug in human history. Nonetheless, drawbacks of AZT are also clear and include required BID dosing and long-term adverse events such as anemia, neutropenia, hepatotoxicity, cardiomyopathy, peripheral neuropathy, and myopathy . Nonetheless, AZT is still used as a generic or as part of Combivir (AZT + 3TC) across Africa and LIMCs. In contrast, TDF + FTC are the favored NRTIs for the vast majority of first line treatment regimens, especially in HICs.
At the JCRC in Kampala, Uganda, we treat approximately 10 000 patients and over 70% have received a first line AZT- or d4T-containing cART. Virologic failures are slightly more common in this LIMC population (~10%/year) of first line treatment failure) than observed in cART treated patients in HICs [53,54,55,56,57]. In 2013, we published a 10-year retrospective study on cART treatment failures and drug resistance of patients attending the JCRC clinics  and reported that HIV-1 drug resistance was only observed in 82% of treatment failures, significantly less than 95% of drug resistance genotypes observed with treatment failure in HICs [13, 56, 57]. Over 60% of these patients failing treatment and lacking drug resistance by Sanger did indeed harbor drug resistance mutations at a frequency of 1 to 20% in the viral quasispecies using a deep sequencing-based HIV-1 genotyping assay (DEEPGEN™HIV) . If these patients did not switch treatment regimen, they continued to fail ongoing treatment despite lack of HIV-1 drug resistance based on Sanger sequencing . Whereas the impact of these minority drug resistance mutations on treatment outcomes in HICs has been very controversial [27,28,29,30,31, 58], clinicians working in LMICs have now recognized that a negative drug resistance profile (by Sanger sequencing) during treatment failures may not rule out the presence of HIV-1 drug resistance [13, 55, 57]. Based on these studies/concerns, the WHO-certified HIV-1 drug resistance testing laboratory at the JCRC has now adopted DEEPGEN™HIV  as a routine test for HIV-1 drug resistance. However, cost of the required equipment, software, and expertise for deep-sequencing based drug resistance tests can be prohibitive for many LIMCs. For that reason, here we used AfriPOLA, an affordable multiplex oligonucleotide assay that can reproducibly detect drug resistant mutations >1% within a patient, and which is effective for subtypes A, C, and D that are responsible for >98% of HIV-1 infections in Africa [59, 60].
Drug resistance to at least one drug in the regimen was observed in 82% of JCRC patients failing first line cART involving AZT or d4T. However, majority of this drug resistance was related to the M184 V mutation conferring 3TC/FTC resistance or various NNRTI mutations whereas only 45% of these same patients had TAMs conferring resistance to AZT or d4T. Second line treatment was provided but AZT or d4T was commonly retained in the regimen if resistance was not detected by Sanger. In this study, majority of these patients lacking TAMs (but harboring 3TC and NNRTI resistant virus) had low frequency TAMs as determined by AfriPOLA and/or DEEPGEN™HIV analyses. In addition, our assays also detected minority mutations associated with resistance to TDF, ABC, ddI, NVP, and EFV (e.g., K65R, L74 V, Y115F, M184 V) that were not identified by Sanger sequencing. These findings suggest that low frequency TAMs (>1%) in the HIV-1 quasispecies could impact treatment outcomes in continued first line or salvage treatments involving AZT or d4T. Moreover, AfriPOLA-based drug resistance testing in Uganda or in other LMICs may be an affordable and easily implemented alternative to Sanger sequencing, which would also provide valuable information on low frequency drug resistance mutations.
As we discussed in a previous study , antiretroviral treatment failure is common in LIMCs which is largely due poor adherence as is the case in HICs. However, behavior-based decisions to stop medications are more common in HIC whereas poor adherence in LIMCs is often related to limited access to clinical centers, high travel costs, intermittent antiretroviral drug supply, and interruptions of funding programs, which all in turn limit access to antiretroviral treatment . Often, interruptions in treatment due to these issues are also associated with missed clinic visits and prolonged absence of treatment, which is rarely recorded in the patient charts. Unlike the HIC study showing an increased in the prevalence of drug resistance mutations in patients remaining on a failing treatment (ref), high percentages of HIV drug resistance is observed upon treatment failure at the JCRC in Uganda when testing is performed on patients with frequent clinic visits (every 2–3 months). However, HIV drug resistance is often absent in those patients missing 1 to 3 visits (3–9 months), a gap which appears to relate a concurrent break in cART. We have shown in this and previous studies  that cART treatment failure without drug resistance at the JCRC typically correlates with the presence of low frequency drug resistance mutations in the patient’s HIV quasispecies. During this extended time between clinic visits (>6 months), the dominance of drug resistant HIV-1 has faded and replaced by the more fit wild type HIV in the patient. Due to the lack of drug resistance based on Sanger genotyping, inability to test for low frequency mutations, and limited treatment options, first line treatment is typically resumed leading to treatment failure . We are currently completing a study involving collection of GIS (geographical information systems) information, adherence rates, and clinic visit frequency to determine factors related to treatment failure. We have previously described that patients typically miss clinic visit during poor adherence, leading to treatment failure and drug resistance.
In summary, dominant AZT and d4T resistance was only observed in 45% of patients failing an AZT- or d4T-containing treatment regimen in Uganda. Slower emergence of TAMs as compared to 3TC and NNRTI drug resistance mutations during virologic failure has been described in terms of historic treatments involving AZT in HICs [51, 61]; however, to our knowledge no studies have documented this low percentage of dominant TAMs upon failure to AZT-containing cART in LMICs. Appearance of TAMs is more hampered by higher fitness costs and higher genetic barriers than appearance of NNRTI resistance mutations which have minimal replicative fitness costs [16, 43, 45,46,47,48,49,50]. Larder et al. reported that appearance of M184 V in AZT + 3TC combination treatments delayed emergence of TAMs due to synergistic loss in fitness and some incompatibility in the resistance mechanisms . Patients failing AZT-containing cART but lacking TAMs by Sanger did have low frequency TAMs when using a deep sequencing or a quantitative oligonucleotide ligation assay. Why low frequency TAMs are found in one half while dominant TAMs are found in the other half of patients failing AZT-/d4T-containing cART remains unanswered. We suspect that low frequency resistance to AZT but dominant resistance to 3TC/FTC and the NNRTI may arise if patients stop only AZT treatment in the first line cART. Although not significant due to the limited sample size, in patients where AZT was administered with 3TC as single pill BID (i.e. Combivir or generic Durovir), there was higher percentages of patients failing with dominant TAMs than in patients failing treatment when AZT was prescribed as a separate pill in a cART regimen. Further studies based on ultrasensitive HIV-1 genotyping assays, such as DEEPGEN™HIV and AfriPOLA, are needed to elucidate the lingering effect of these minority HIV-1 resistance variants in patients from LMICs. Unfortunately, the fact remains that AZT (or zidovudine mostly as a generic) is still commonly used throughout the world due to low cost and despite the lack of AZT use in first line regimens in high income countries. Given the lower AZT adherence, toxicity, and higher frequency of failures, most of the various international aid organizations have dropped the use of AZT and d4T for the TDF-containing regimens. However, disparities in treatment still exist when considering HICs commonly utilize dolutegravir in place of efavirenz in first line regimens.
combination antiretroviral therapy
High income countries
Lower middle-income countries
Non-nucleoside RT inhibitors
Thymidine Analog Mutations
World Health Organization
UNAIDS. The HIV and AIDS Uganda Country Progress Report 2014. Geneva, Switzerland. 2014; http://www.unaids.org/sites/default/files/country/documents/UGA_narrative_report_2015.pdf.
Paredes R, Clotet B. Clinical management of HIV-1 resistance. Antivir Res. 2010;85:245–65.
Zolopa AR. The evolution of HIV treatment guidelines: current state-of-the-art of ART. Antivir Res. 2010;85:241–4.
Cohen J. Breakthrough of the year. HIV treatment as prevention. Science. 2011;334:1628.
Menendez-Arias L. Molecular basis of human immunodeficiency virus type 1 drug resistance: overview and recent developments. Antivir Res. 2013;98:93–120.
WHO. Consolidated guidelines on the use of antiretroviral drugs for treating and preventing HIV infection. 2013. http://appswhoint/iris/bitstream/10665/85321/1/9789241505727_engpdf?ua=1.
Pujades-Rodriguez M, Dantony E, Pinoges L, Ecochard R, Etard JF, Carrillo-Casas E, et al. Toxicity associated with stavudine dose reduction from 40 to 30 mg in first-line antiretroviral therapy. PLoS One. 2011;6:e28112.
Menezes CN, Maskew M, Sanne I, Crowther NJ, Raal FJ. A longitudinal study of stavudine-associated toxicities in a large cohort of south African HIV infected subjects. BMC Infect Dis. 2011;11:244.
Castelnuovo B, Kiragga A, Kamya MR, Manabe Y. Stavudine toxicity in women is the main reason for treatment change in a 3-year prospective cohort of adult patients started on first-line antiretroviral treatment in Uganda. J Acquir Immune Defic Syndr. 2011;56:59–63.
WHO. HIV drug resistance report. HIV/AIDS Prorgramme. 2012.
Cozzi-Lepri A, Phillips AN, Ruiz L, Clotet B, Loveday C, Kjaer J, et al. Evolution of drug resistance in HIV-infected patients remaining on a virologically failing combination antiretroviral therapy regimen. AIDS. 2007;21:721–32.
Richard N, Juntilla M, Abraha A, Demers K, Paxinos E, Galovich J, et al. High prevalence of antiretroviral resistance in treated Ugandans infected with non-subtype B human immunodeficiency virus type 1. AIDS Res Hum Retrovir. 2004;20:355–64.
Kyeyune F, Nankya I, Metha S, Akao J, Ndashimye E, Tebit DM, et al. Treatment failure and drug resistance is more frequent in HIV-1 subtype D versus subtype A-infected Ugandans over a 10-year study period. AIDS. 2013;27:1899–909.
Reynolds SJ, Laeyendecker O, Nakigozi G, Gallant JE, Huang W, Hudelson SE, et al. Antiretroviral drug susceptibility among HIV-infected adults failing antiretroviral therapy in Rakai, Uganda. AIDS Res Hum Retroviruses. 2012;28:1739–44.
Ssemwanga D, Kapaata A, Lyagoba F, Magambo B, Nanyonjo M, Mayanja BN, et al. Low drug resistance levels among drug-naive individuals with recent HIV type 1 infection in a rural clinical cohort in southwestern Uganda. AIDS Res Hum Retrovir. 2012;28:1784–7.
Arts EJ, Quinones-Mateu ME, Albright JL, Marois JP, Hough C, Gu Z, et al. 3′-Azido-3′-deoxythymidine (AZT) mediates cross-resistance to nucleoside analogs in the case of AZT-resistant human immunodeficiency virus type 1 variants. J Virol. 1998;72:4858–65.
Larder BA. 3′-Azido-3′-deoxythymidine resistance suppressed by a mutation conferring human immunodeficiency virus type 1 resistance to nonnucleoside reverse transcriptase inhibitors. Antimicrob Agents Chemother. 1992;36:2664–9.
Quinones-Mateu ME, Arts EJ. HIV-1 fitness: implications for drug resistance, disease progression, and global epidemic evolution. In: Kuiken C, Foley B, Hahn B, Marx P, McCutchan F, Mellors J, Wolinsky S, Korber B, editors. HIV sequence compendium 2001. Los Alamos: Theoretical Biology and Biophysics Group, Los Alamos National Laboratory; 2001. p. 134–70.
Quinones-Mateu ME, Arts EJ. Virus fitness: concept, quantification, and application to HIV population dynamics. Curr Top Microbiol Immunol. 2006;299:83–140.
Kellam P, Boucher CA, Tijnagel JM, Larder BA. Zidovudine treatment results in the selection of human immunodeficiency virus type 1 variants whose genotypes confer increasing levels of drug resistance. J Gen Virol. 1994;75(Pt 2):341–51.
Paquet AC, Baxter J, Weidler J, Lie Y, Lawrence J, Kim R, et al. Differences in reversion of resistance mutations to wild-type under structured treatment interruption and related increase in replication capacity. PLoS One. 2011;6:e14638.
Hu Z, Giguel F, Hatano H, Reid P, Lu J, Kuritzkes DR. Fitness comparison of thymidine analog resistance pathways in human immunodeficiency virus type 1. J Virol. 2006;80:7020–7.
Kyeyune F, Gibson RM, Nankya I, Venner C, Metha S, Akao J, et al. Low-frequency drug resistance in HIV-infected Ugandans on antiretroviral treatment is associated with regimen failure. Antimicrob Agents Chemother. 2016;60:3380–97.
Larder BA, Kohli A, Kellam P, Kemp SD, Kronick M, Henfrey RD. Quantitative detection of HIV-1 drug resistance mutations by automated DNA sequencing. Nature. 1993;365:671–3.
Church JD, Jones D, Flys T, Hoover D, Marlowe N, Chen S, et al. Sensitivity of the ViroSeq HIV-1 genotyping system for detection of the K103N resistance mutation in HIV-1 subtypes a, C, and D. J Mol Diagn. 2006;8:430–2. quiz 527
Leitner T, Halapi E, Scarlatti G, Rossi P, Albert J, Fenyo EM, et al. Analysis of heterogeneous viral populations by direct DNA sequencing. BioTechniques. 1993;15:120–7.
Li JZ, Kuritzkes DR. Clinical implications of HIV-1 minority variants. Clin Infect Dis. 2013;56:1667–74.
Stekler JD, Ellis GM, Carlsson J, Eilers B, Holte S, Maenza J, et al. Prevalence and impact of minority variant drug resistance mutations in primary HIV-1 infection. PLoS One. 2011;6:e28952.
Li JZ, Paredes R, Ribaudo HJ, Kozal MJ, Svarovskaia ES, Johnson JA, et al. Impact of minority nonnucleoside reverse transcriptase inhibitor resistance mutations on resistance genotype after virologic failure. J Infect Dis. 2013;207:893–7.
Li JZ, Paredes R, Ribaudo HJ, Svarovskaia ES, Kozal MJ, Hullsiek KH, et al. Relationship between minority nonnucleoside reverse transcriptase inhibitor resistance mutations, adherence, and the risk of virologic failure. AIDS. 2012;26:185–92.
Codoner FM, Pou C, Thielen A, Garcia F, Delgado R, Dalmau D, et al. Added value of deep sequencing relative to population sequencing in heavily pre-treated HIV-1-infected subjects. PLoS One. 2011;6:e19461.
Gianella S, Richman DD. Minority variants of drug-resistant HIV. J Infect Dis. 2010;202:657–66.
Li JZ, Paredes R, Ribaudo HJ, Svarovskaia ES, Metzner KJ, Kozal MJ, et al. Low-frequency HIV-1 drug resistance mutations and risk of NNRTI-based antiretroviral treatment failure: a systematic review and pooled analysis. JAMA. 2011;305:1327–35.
Swenson LC, Mo T, Dong WW, Zhong X, Woods CK, Thielen A, et al. Deep V3 sequencing for HIV type 1 tropism in treatment-naive patients: a reanalysis of the MERIT trial of maraviroc. Clin Infect Dis. 2011;53:732–42.
Swenson LC, Mo T, Dong WW, Zhong X, Woods CK, Jensen MA, et al. Deep sequencing to infer HIV-1 co-receptor usage: application to three clinical trials of maraviroc in treatment-experienced patients. J Infect Dis. 2011;203:237–45.
Goodman DD, Zhou Y, Margot NA, McColl DJ, Zhong L, Borroto-Esoda K, et al. Low level of the K103N HIV-1 above a threshold is associated with virological failure in treatment-naive individuals undergoing efavirenz-containing therapy. AIDS. 2011;25:325–33.
Gibson RM, Meyer AM, Winner D, Archer J, Feyertag F, Ruiz-Mateos E, et al. Sensitive deep sequencing-based HIV-1 genotyping assay to simultaneously determine susceptibility to protease, reverse transcriptase, Integrase, and maturation inhibitors, as well as HIV-1 Coreceptor tropism. Antimicrob Agents Chemother. 2014;58:2167–85.
Lalonde MS, Arts EJ. DNA suspension arrays: silencing discrete artifacts for high-sensitivity applications. PLoS One. 2010;5:e15476.
Halvas EK, Aldrovandi GM, Balfe P, Beck IA, Boltz VF, Coffin JM, et al. Blinded, multicenter comparison of methods to detect a drug-resistant mutant of human immunodeficiency virus type 1 at low frequency. J Clin Microbiol. 2006;44:2612–4.
Grant RM, Kuritzkes DR, Johnson VA, Mellors JW, Sullivan JL, Swanstrom R, et al. Accuracy of the TRUGENE HIV-1 genotyping kit. J Clin Microbiol. 2003;41:1586–93.
Wensing AM, Calvez V, Gunthard HF, Johnson VA, Paredes R, Pillay D, et al. 2015 update of the drug resistance mutations in HIV-1. Top Antivir Med. 2015;23:132–41.
Fischl MA, Richman DD, Griego MH, Gottlieb MS, Volberding PA, Laskin OL, et al. The efficacy of azidothymidine (AZT) in the treatment of patients with AIDS and AIDS-related complex. N Engl J Med. 1987;317:185–91.
Larder BA, Kemp SD. Multiple mutations in HIV-1 reverse transcriptase confer high-level resistance to zidovudine (AZT). Science. 1989;246:1155–8.
Larder BA, Darby G, Richman DD. HIV with reduced sensitivity to zidovudine (AZT) isolated during prolonged therapy. Science. 1989;243:1731–4.
St Clair MH, Martin JL, Tudor-Williams G, Bach MC, Vavro CL, King DM, et al. Resistance to ddI and sensitivity to AZT induced by a mutation in HIV-1 reverse transcriptase. Science. 1991;253:1557–9.
Kellam P, Boucher CA, Larder BA. Fifth mutation in human immunodeficiency virus type 1 reverse transcriptase contributes to the development of high-level resistance to zidovudine. Proc Natl Acad Sci U S A. 1992;89:1934–8.
Larder BA, Kemp SD, Harrigan PR. Potential mechanism for sustained antiretroviral efficacy of AZT-3TC combination therapy. Science. 1995;269:696–9.
Meyer PR, Matsuura SE, Mian AM, So AG, Scott WA. A mechanism of AZT resistance: an increase in nucleotide-dependent primer unblocking by mutant HIV-1 reverse transcriptase. Mol Cell. 1999;4:35–43.
Arion D, Kaushik N, McCormick S, Borkow G, Parniak MA. Phenotypic mechanism of HIV-1 resistance to 3′-azido-3′-deoxythymidine (AZT): increased polymerization processivity and enhanced sensitivity to pyrophosphate of the mutant viral reverse transcriptase. Biochemistry. 1998;37:15908–17.
Arion D, Sluis-Cremer N, Parniak MA. Mechanism by which phosphonoformic acid resistance mutations restore 3′-azido-3′-deoxythymidine (AZT) sensitivity to AZT-resistant HIV-1 reverse transcriptase. J Biol Chem. 2000;275:9251–5.
Rey D, Hughes M, Pi JT, Winters M, Merigan TC, Katzenstein DA. HIV-1 reverse transcriptase codon 215 mutation in plasma RNA: immunologic and virologic responses to zidovudine. The AIDS Clinical Trials Group study 175 virology team. J Acquir Immune Defic Syndr Hum Retrovirol. 1998;17:203–8.
Margolis AM, Heverling H, Pham PA, Stolbach A. A review of the toxicity of HIV medications. J Med Toxicol. 2014;10:26–39.
Braitstein P, Brinkhof MW, Dabis F, Schechter M, Boulle A, Miotti P, et al. Mortality of HIV-1-infected patients in the first year of antiretroviral therapy: comparison between low-income and high-income countries. Lancet. 2006;367:817–24.
Cescon AM, Cooper C, Chan K, Palmer AK, Klein MB, Machouf N, et al. Factors associated with virological suppression among HIV-positive individuals on highly active antiretroviral therapy in a multi-site Canadian cohort. HIV Med. 2011;12:352–60.
Pillay P, Ford N, Shubber Z, Ferrand RA. Outcomes for efavirenz versus nevirapine-containing regimens for treatment of HIV-1 infection: a systematic review and meta-analysis. PLoS One. 2013;8:e68995.
Hamers RL, Sigaloff KC, Wensing AM, Wallis CL, Kityo C, Siwale M, et al. Patterns of HIV-1 drug resistance after first-line antiretroviral therapy (ART) failure in 6 sub-Saharan African countries: implications for second-line ART strategies. Clin Infect Dis. 2012;54:1660–9.
Amoroso A, Etienne-Mesubi M, Edozien A, Ojoo S, Sheneberger R, Obiefune M, et al. Treatment outcomes of recommended first-line antiretroviral regimens in resource-limited clinics. J Acquir Immune Defic Syndr. 2012;60:314–20.
Kagan RM, Johnson EP, Siaw M, Biswas P, Chapman DS, Su Z, et al. A genotypic test for HIV-1 tropism combining Sanger sequencing with ultradeep sequencing predicts virologic response in treatment-experienced patients. PLoS One. 2012;7:e46334.
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.
Rhee SY, Blanco JL, Jordan MR, Taylor J, Lemey P, Varghese V, 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 Med. 2015;12:e1001810.
Katzenstein DA, Hammer SM, Hughes MD, Gundacker H, Jackson JB, Fiscus S, et al. The relation of virologic and immunologic markers to clinical outcomes after nucleoside therapy in HIV-infected adults with 200 to 500 CD4 cells per cubic millimeter. AIDS Clinical Trials Group study 175 virology study team. N Engl J Med. 1996;335:1091–8.
We thank the JCRC patients for participating in this study as well as physicians, nurses, and counselors who continue to provide the best care treatment at the JCRC. This research was presented in part at the 25th Canadian Conference on HIV/AIDS Research, Winnipeg, Canada (2016).
F.K. was funded by a scholarship from the NIH Fogarty International Center grants D43-TW000011 and D43-TW009780. E.J.A. was funded by NIAID/NIH AI49170 and holds the Canada Research Chair in HIV-1 Pathogenesis and Viral Control. M.E.Q-M was partially supported by the CWRU/UH Center for AIDS Research (P30 AI036219) and by funding from University Hospitals Cleveland Medical Center (UHCMC) for the University Hospitals Translational Laboratory (UHTL).
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Clinical and virological data were obtained from patient care database at the JCRC under IRB approval (EM10-07) for HIV-1 drug resistance testing.
Consent for publication
E.J.A developed and patented the technology used in AfriPOLA. R.M.G and M.E.Q-M developed DEEPGEN™HIV. There are no further patents, products in development or marketed products to declare.
About this article
Cite this article
Gibson, R.M., Nickel, G., Crawford, M. et al. Sensitive detection of HIV-1 resistance to Zidovudine and impact on treatment outcomes in low- to middle-income countries. Infect Dis Poverty 6, 163 (2017). https://doi.org/10.1186/s40249-017-0377-0
- Antiretroviral treatment
- Drug resistance