High level of HIV-1 drug resistance mutations in patients with unsuppressed viral loads in rural northern South Africa
© The Author(s) 2017
Received: 26 May 2017
Accepted: 19 July 2017
Published: 27 July 2017
Combination antiretroviral therapy (cART) has significantly reduced HIV morbidity and mortality in both developed and developing countries. However, the sustainability of cART may be compromised by the emergence of viral drug resistance mutations (DRM) and the cellular persistence of proviruses carrying these DRM. This is potentially a more serious problem in resource limited settings.
DRM were evaluated in individuals with unsuppressed viral loads after first or multiple lines of cART at two sites in rural Limpopo, South Africa. Seventy-two patients with viral loads of >1000 copies/ml were recruited between March 2014 and December 2015. Complete protease (PR) and partial Reverse Transcriptase (RT) sequences were amplified from both plasma RNA and paired proviral DNA from 35 of these subjects. Amplicons were directly sequenced to determine subtype and DRM using the Stanford HIV Drug Resistance Interpretation algorithm.
Among the 72 samples, 69 could be PCR amplified from RNA and 35 from both RNA and DNA. Sixty-five (94.2%) viruses were subtype C, while one was subtype B (1.4%), one recombinant K/C, one recombinant C/B and one unclassified. Fifty-eight (84%) sequences carried at least one DRM, while 11 (15.9%) displayed no DRM. DRM prevalence according to drug class was: NRTI 60.8% NNRTI 65.2%, and PI 5.8%. The most common DRMs were; M184V (51.7%), K103N (50%), V106M (20.6%), D67N (13.3%), K65R (12%). The frequency of the DRM tracked well with the frequency of use of medications to which the mutations were predicted to confer resistance. Interestingly, a significant number of subjects showed predicted resistance to the newer NNRTIs, etravirine (33%) and rilpivirine (42%), both of which are not yet available in this setting. The proportion of DRM in RNA and DNA were mostly similar with the exception of the thymidine analogue mutations (TAMs) D67N, K70R, K219QE; and K103N which were slightly more prevalent in DNA than RNA. Subjects who had received cART for at least 5 years were more likely to harbour >2 DRM (p < 0.05) compared to those treated for a shorter period. DRM were more prevalent in this rural setting compared to a neighbouring urban setting.
We found a very high prevalence of NRTI and NNRTI DRM in patients from rural Limpopo settings with different durations of treatment. The prevalence was significantly higher than those reported in urban settings in South Africa. The dominance of NNRTI based mutations late in treatment supports the use of PI based regimens for second line treatment in this setting. The slight dominance of TAMs in DNA from infected PBMCs compared to plasma virus requires further studies that should include cART subjects with suppressed virus. Such studies will improve our understanding of the pattern of drug resistance and dynamics of viral persistence in these rural settings.
KeywordsHIV drug resistance Antiretroviral therapy Peripheral blood mononuclear cells Plasma Rural South Africa
HIV drug resistance remains a major threat to the success of scaling up of combination antiretroviral therapy (cART) in developing countries, especially in Africa where about 60% of HIV-infected individuals were on cART in 2012 . The increased availability of cART has led to remarkable treatment results in most programs in sub-Saharan Africa (SSA), leading to a significant reduction of new infections . Despite this success, some patients with incomplete adherence to treatment protocols might fail therapy and develop resistant viruses . One challenge of cART in SSA has been the inability to routinely monitor virological parameters such as viral load (VL) and drug resistance, which are important to guide patient management and choose the right salvage treatments [4, 5]. Recent WHO recommendations stipulate that VL monitoring be utilized as the preferred method to determine treatment failure in SSA . Within the past years, VL monitoring has in fact been successfully implemented as part of care in South Africa .
Sequencing is the most accurate method to track the emergence of DRM and is typically done by Sanger sequencing. However, cost and availability of infrastructure to perform these assays still remain major obstacles . As a consequence, most patients failing therapy remain on first-line therapy for long periods causing an accumulation of mutations, which might lead to the emergence of high-level resistance [8, 9]. This is particularly true for patients who started treatment long before VL was more widely available. Samples from such long term cART patients provide a good source to track previous DRM which can be found only in the cellular compartment.
Genotypic drug resistance studies of subjects with unsuppressed VL have been performed in several SSA countries including South Africa and suggest a high level of DRM as well as some subtype specific mutations [10–16]. For example, previous studies among South Africans reported that infections with HIV-1 subtype C may lead to a rapid emergence of certain non-nucleoside reverse transcriptase (NNRTI) or nucleotide reverse transcriptase inhibitor (NRTI) DRM, which are less common among other subtypes. Specifically, V106M (NNRTI), K65R (NRTI) and N348I are more common among subtype C strains compared to other group M strains [17–19].
The South African treatment protocols are consistent with the 2016 WHO guidelines, which recommend a combination of two NRTIs and one NNRTIs as a first-line regimen and two-NRTI and one protease inhibitor (PI) as second line. According to the April 2014 South African revised guidelines, a preferred first-line regimen could be: tenofovir (TDF)/lamivudine (3TC) + efavirenz (EFV)/nevirapine (NVP), or alternatively abacavir (ABC) + 3TC + EFV; or stavudine (d4T) + 3TC + EFV. For second line treatment, TDF + 3TC + lopinavir/ritonavir (LPV/r); or 3TC + AZT + LPV/r could be used.
Despite the benefits of cART, a significant number of patients still experience treatment failure, associated with the development of DRM due to incomplete adherence to therapy or the presence of transmitted DRM [20, 21]. Given the changing treatment strategies and options, genotypic and phenotypic drug resistance patterns will continue to evolve. Although considerable data on virologic failure in patients with unsuppressed VL on first or second line treatment from South Africa are available, these studies have been mostly done in urban settings [15, 22, 23]. On the contrary, rural settings (making up almost 40% of the South African population) still remain largely understudied. Such rural populations often face significant challenges such as distance to treatment sites and access to virological monitoring. In addition, cultural beliefs often increase stigma and could impact adherence and prevalence of DRM [24, 25]. Two recent studies in Uganda and South Africa reported that the intensity of virological monitoring, frequency and pattern of DRM is influenced by the location (urban or rural) of the treatment site [26, 27]. Such results illustrate that inadequate infrastructure and low per capita income common in rural settings might influence the successful implementation of cART. It is therefore imperative that frequent monitoring of DRM is implemented in order to better understand the relationship between drug resistance development and therapy failures in these rural settings.
The use of integrated proviral DNA rather than RNA isolated from circulating virus to monitor HIV drug resistance in patients failing treatment has been considered as an alternative method for drug resistance genotyping . Sequencing of DNA is also an efficient method to study the archived and persistent viral reservoirs in patients with low or suppressed viral loads . With the current focus on HIV persistence, it is important that studies investigating persisting viral populations during suppression are performed to better understand the viral dynamics and alternative future treatment options.
In this study, we describe the drug resistance profile of HIV-1 infected individuals experiencing unsuppressed VL attending two rural HIV treatment centres; the HIV/AIDS Prevention Group Wellness Clinic (HAPG) in Bela-Bela and the Donald Frazer Hope Clinic (DFHC) in Vhufhuli, both located in the Limpopo Province of northern South Africa. Our findings show a high prevalence of NNRTI and NRTI, but low level PI mutations. DRM were more prevalent in this rural setting than reported for a nearby urban setting in Pretoria. Furthermore, by comparing the DRM (NRTI and NNRTI) in viral RNA as well as proviral DNA among circulating viruses in plasma and peripheral blood compartments in 35 subjects, we found that the DRM in these two compartments were mostly similar with a slight dominance of thymidine analogue mutations (TAMs) in DNA, a result that necessitates further studies with a larger cohort.
Subjects and specimen collection
From July 2013 to December 2015, blood specimens were collected from participating subjects who comprised a random selection of those with unsuppressed viremia at the HIV/AIDS Prevention Group Wellness Clinic (HAPG) in Bela-Bela and the Donald Frazer Hope Clinic (DFHC) in Vhufhuli, in Limpopo Province of rural northern South Africa. All recruited patients met the study criteria of plasma VL > 1000 copies/ml in one of the following ways: (i) two consecutive VL greater than 1000 copies/ml after previous suppression (ii) one VL > 1000 copies/ml after previous suppression followed by a change in treatment (iii) one VL > 1000 copies/ml after 6 months on ART without suppression.
CD4 count, viral load determination and sample processing
VL and CD4+ T cell counts were all performed by the National Health Laboratory Services (NLS) of South Africa. VL were done annually as stipulated by the 2015 South African Treatment Guideline and involved clinic staff usually calling to remind patients of their next VL testing. For drug resistance determination, 5 ml of whole blood was collected from each study participant and transported to the HIV/AIDS and Global Health Research and Training Laboratory at the Department of Microbiology, University of Venda for processing and analyses. On arrival in the laboratory, plasma was separated from cells by centrifuging at 3000 rpm for 5 min and stored at −80 °C in several aliquots. Total cells comprising mostly buffy coats were stored at 4 °C until further processing.
DNA and RNA extractions
Plasma RNA and proviral DNA from total cells were extracted using Qiagen Viral RNA Mini Kit and Qiagen blood DNA Midi Kit (Qiagen, Valencia, CA, USA), respectively, according to the manufacturers’ protocols. Extracted DNA and RNA from these samples were then used for PCR and RT-PCR, respectively.
Polymerase chain reaction and sequencing
A fragment of the pol gene approximately 1.5 kb in size (comprising the complete PR and the first 1200 bp of RT) was amplified using a nested PCR protocol described previously [10, 30]. Briefly, cDNA was prepared from viral RNA and used for a nested PCR . Cycling conditions for cDNA and PCRs are available upon request. Amplicons were purified using the Qiaquick PCR purification kit (Qiagen, Valencia, CA, USA) according to the manufacturers’ protocol.
Direct population based Sanger sequencing was performed on an ABI 3730XL platform using the ABI Prism® BigDye™ Terminator v3.1 ready reaction sequencing kit (Applied Biosystems, Warrington, UK). The resulting forward and reverse nucleotide sequence electro-pherograms were assembled, manually edited and translated into predicted amino acids with SeqMan Pro II software, version 8.0 (DNASTAR, Lasergene, USA). Primer sequences used for amplification and sequencing are available upon request.
Phylogenetic analyses were performed to determine the viral subtypes in relation to reference subtypes from the Los Alamos National Laboratory. Sequence alignment of the gene fragments was performed using MUSCLE in MEGA 6 software . The phylogenetic trees were derived using the Maximum Likelihood Model, which calculates the difference in transition and transversion rates. The reliability of the tree was assessed by bootstrapping of 1000 replicates.
Drug resistance genotyping
For drug resistance genotyping, the newly derived sequences were submitted to the Stanford HIV Drug Resistance Database (http://hivdb.stanford.edu) for examination of drug resistance mutations. Observed predicted DRM were confirmed using the International AIDS Society-USA (IAS) Panel for drug resistance guideline . All sequences associated with high level resistance according to Stanford and IAS algorithms were considered to be resistant.
The Wilcoxon rank sum test was used to compare groups. P-value <0.05 was considered to be significant. Correlation between the rural and urban settings was performed using the Pearson Correlation.
Study participants, clinical data and study sites
Demographic, clinical and virological characteristics of 72 patients on ART experiencing an unsuppressed viral load
Number of patients
Median age (range) in years
Median CD4 count (range) in cells/mm3
Median plasma viral load (range) in copies/ml
ART regimen at time of genotyping
Time on ART (months)
Clinical stage (CDC) (%)
Amplified Genotypes (%)
Others (K/C, C/B and unclassified)
HIV genotyping and drug resistance prediction
Three subjects carried the accessory mutation A62V that has been reported to be widespread among subtype A viruses from countries of the former Soviet Union . A62V often occurs in combination with multi-resistance NRTI mutations K65R and Q151M. In this study, all the A62V mutations described were linked to K65R, but not to Q151M . Out of the five patients who carried multiple thymidine analogue mutations (TAMs) only one had the K65R mutation, suggesting that these TAMs likely resulted from previous AZT or d4T treatment and not exposure to TDF. TAMs were dominated by D67N (13.8%), followed by K219Q/E (12.1%) and K70R (12.1%) with a low occurrence of T215Y/F (6.9%). The T215Y/F mutation seem to be subtype specific as it was recently found as the dominant TAM in 33.6% of non-subtype C infected study subjects in five Central African countries . Interestingly, among the viruses with mutations at this position, tyrosine (Y) was the dominant amino acid (6.9%; T215Y) compared to the complete absence of phenylalanine (F) (0%; T215F; Fig. 1). Five percent of viruses carried an isoleucine (I, encoded by ATT, ATC) mutation at position 215 (215I) which resulted from a single nucleotide changing threonine (ACT, ACC) to isoleucine. The 215I mutation seems to be transitional (intermediate) to T215F (phenylalanine-TTT, TTC).
Comparison of the NRTI and NNRTI DRM prevalence from Limpopo (rural), KwaZulu-Natal (rural) and urban (Pretoria)
p-value (Limpopo vs. Pretoria)b
This study was performed among patients failing treatment in two treatment centres in rural South Africa: the HAPG Wellness Clinic which started cART administration in 2000, 4 years before the roll-out of antiretroviral therapy in the public health sector in South Africa; and the Donald Fraser Clinic which started providing treatment in 2004. Although both sites have been implementing VL and CD4+ T-cell count monitoring, according to national guidelines, and as part of HIV management for several years, the prevalence of DRM was still high. The M184V mutation was the most prevalent DRM, observed in 52% of study subjects. This coincided with the extensive use of 3TC, surpassed only by the use of TDF and EFV in the population. Interestingly, TDF mutations were one of the least common DRM in the study despite its frequent usage, supporting its inclusion as the main backbone during first line treatment.
Although currently used cART is effective for all group M subtypes, the pattern of drug resistance differs between subtypes [10, 36]. Our study confirmed the dominance of K65R and V106M, two mutations which are most common among subtype C strains [36, 37]. In fact, very recently in a cohort of more than 3700 subjects failing cART in several West and Central Africa countries, these two DRM were found to be more prevalent among subtypes C strains compared to CRF02_AG, CRF06cpx, and subtypes A and G [17, 18, 36]. This dominance has been attributed to the genomic sequence of subtype C, which requires only a single nucleotide change at these codon positions to transition to a drug resistant amino acid residue . Similarly, we speculate that the occurrence of the E138A/K/Q mutation could be due to the codon usage of subtype C RT at this position which is GAA versus GAG in other group M strains. In fact, in the Stanford Drug Resistance database, subtype C strains from drug exposed patients have a higher prevalence of the E138A/K/Q mutation compared to other group M strains (8–9.5% vs. 2–5%). Further verification of such subtype-specific differences in the pattern of resistance mutations would clearly be of great significance for treatment management and also for a better understanding of the overall resistance mechanism.
The high occurrence of both K103N and V106M could be a reflection of the wide use of EFV and NVP as part of the first line therapy. In addition, the high levels of V106M mutations could also be attributed to subtype-specificity, since this mutation had been found most frequently among subtype C [17, 36]. It seems quite clear that NNRTIs should not be used for second line treatment in rural South African populations. This conclusion is supported by the fact that a significant number of the patients we studied already show resistance to newer NNRTIs (ETR and RPV), which have not even been introduced in South Africa.
The selection of TAMs such as D67N, M41L and T215Y could be attributed to treatment with the thymidine analogues AZT and d4T, which also cause cross resistance to other NRTIs . Similarly, patients harbouring TAMs also had high level resistance to NNRTIs (EFV and NVP) and 3TC supporting the point that they were likely failing thymidine containing ART at the point when TDF replaced AZT in the national guidelines . TAMs have been reported to be more common among CRF06cpx strains than other HIV-1 subtypes or CRFs [10, 36]. In this study, there was a noticeable absence among subjects of the L210W and very low occurrence of the M41L DRM, both components of the Type 1 TAMs pathway (M41L, L210W and T215Y). Although, this could likely be related to the reduced usage of d4T, recent reports suggest that compared to other dominant HIV-1 genotypes (subtypes A, CRF02_AG and CRF06cpx), subtype C show the lowest prevalence of M41L DRM . The absence of M41L and L210W was initially reported from a subtype C cohort in Botswana where the most common TAM pathway was D67N, K70R and T215Y . In our study, we observed that the four TAMs (D67N, K70R, T215Y and K219Q/E) which make up the Type 2 TAM pathway, occurred more frequently and correlated well with the usage of AZT. Although we could not state categorically the order of emergence of these mutations as this was a cross sectional study, the mutations D67N and K70R were the most frequent single occurring TAMs. Both AZT and d4T have been reported to cause TAMs at an equal rate .
The majority of participants in this study were adults >18 years old. This is not surprising given that treatment of children and adolescents with cART has not been vigorously pursued in SSA until recently. Out of the 13 subjects <18 years old with unsuppressed virus, DRM could be detected in 70% (6 were TAMs and 3 were either NNRTI or other NRTI). DRM has been reported to be highly prevalent among African children with unsuppressed VL . In Uganda and South Africa, DRMs were higher in children from rural compared to urban clinics indicating a better management and/or easier adherence in urban clinics [26, 27]. Due to the low number of children and adolescents in this study, we could not independently confirm this. With the increasing usage of cART in children and adolescents in SSA, more longitudinal studies are required to better understand the pattern of DRM among this group.
The observed high levels of DRM in these rural settings could be linked to several factors. First, the availability of transportation to treatment centers is generally significantly reduced in rural areas compared to urban settings. Secondly, shortage of cART supplies are common occurrences in rural areas. Finally, cultural differences or stigma might affect the open administration of cART especially among women, a phenomenon that is more strongly present in rural areas. Another factor that likely influenced the high rates of DRM in one of the study sites (HAPG) was the inadequate planning and poor follow-up by health authorities. In 2009, at the HAPG privately run facility in Bela-Bela, patients’ healthcare providers were switched as patients were transferred from HAPG to a newly opened nearby Public Health Clinic (PHC). Records dating from the beginning of cART therapy at HAPG in 2004 showed good patient follow-up and less treatment failure. However, after the massive patient transfer in 2009 to the PHC and their return to the HAPG private clinic in 2012, the rate of virologic failure increased significantly among returning patients. Despite the high prevalence of DRM observed in this study, the number of patients with detectable VL, but no DRM was similar to previous reports (16% vs. 19–36% respectively) in some urban centres in South Africa [22, 36, 42, 43].
There are few reports on the prevalence and pattern of DRM in rural South Africa [26, 30, 44]. A recent study compared the frequency of DRM from rural KwaZulu-Natal and urban Pretoria, and found a clear difference between these two sites . In our current study, we confirmed the findings of Rossouw et al. Notably, the number of DRMs at the two sites in rural Limpopo were generally higher than in neighboring urban Pretoria (Table 2). Only two DRMs were substantially more common in Pretoria as compared with the Limpopo population. The mutations V106A/M and Y181C were 13.3 and 4.8% more common in Pretoria than in Limpopo, respectively. Both mutations are associated with exposure to NVP . At the time of sequencing, the Pretoria subjects were much more likely to be exposed to NVP than were the subjects in Limpopo (32.4% vs. 5%). Although subjects in these rural settings would be predicted to respond to the ART regimens that are currently in use, we suggest that the WHO recommendation that includes a PI in the second line treatment and close monitoring be strictly observed in order to reduce the transmission of DRM. This recommendation is supported by the very low level of PI DRMs observed in this study.
We compared DRM from both plasma (RNA) and cellular (proviral DNA) compartments from patients experiencing unsuppressed viral load who had received treatment for varying periods, ranging from 1 to 10 years. Plasma is routinely used for drug resistance testing, because viral RNA reflects the current replicating virus in blood. The prevalence of DRM in RNA and DNA was mostly similar with the exception of the TAMs, which were more prevalent in DNA than RNA. The persistence of some TAMs in DNA could be explained by the formerly common administration of NRTI (AZT and d4T) regimens. Of note is the fact that d4T is gradually being phased out and was not included in any of the current treatment regimens received by the study subjects. Three viral RNA and two proviral DNA derived sequences from patients who had never been treated with ddI, ABC, or TDF harboured a K65R mutation.
It is important to point out that our analyses were done using routine population-based Sanger sequencing that does not allow the detection of minority viral variants which comprise less than 20% of the viral population in a patient sample [46, 47]. Recent methods such as next generation sequencing would most likely detect minor drug resistant variants at levels as low as 1% [47, 48]. Considering the fact that undetected drug resistant minority variants may persist when cART is discontinued or changed, it will be worthwhile to perform comparative RNA/DNA drug resistance studies using more sensitive next generation sequencing . Although advances in cART has improved the lives of HIV infected patients, current regimens have failed to eliminate latent reservoirs, and it is acknowledged that eradication of the virus is not possible with current cART regimens alone . Furthermore, the efficacy of switching drugs for patients experiencing a rebound in VL will depend mainly on the number of DRM persisting in the proviral reservoir, following previous therapeutic failures [50, 51].
Overall, these results clearly show that patients receiving cART, but experiencing therapeutic failure, harbour diverse drug-resistant variants, which in some cases persisted in the cellular blood compartment. Detection of this persisting variants can be improved by using more sensitive methods such as next generation sequencing. The proviral variants observed in this study closely resemble the population found in the plasma compartment (viral RNA) with the exception of some persistent TAMs; confirming the known fact that proviral DNA constitutes a reservoir for drug-resistant variants which might replenish plasma pool during suboptimal therapy.
Two limitations should be noted in our investigation: first, this was a cross-sectional study in which data on adherence to treatment was not available; as a result, it was not possible to associate the observed virologic failure to poor adherence. Second, information on prior exposure to NVP in the early days of mother-to-children transmission (MTCT) prevention efforts was also not available. The use of nevirapine monotherapy could have potentially impacted on the observed level of NNRTI mutations particularly K103N.
In conclusion, this study suggests that there is a high level of DRM in rural South Africa and confirms recent observations from KwaZulu-Natal that DRMs are more prevalent in rural than urban areas. It also confirms that long term ART with no virological monitoring could lead to an accumulation of DRMs especially of the NRTI class. Although plasma and PBMCs presented largely identical DRMs, the later compartment may act as an archive of drug resistant variants, especially TAMs, an observation that needs further studies with a larger sample size including subjects with suppressed virus. Such studies could shed more light on the persistence of DRMs and their impact on HIV treatment and cure.
DMT, PB, DR and MLH designed the study; EME, CM, LM, KMG were involved in sampling and performing experiments; PJ interpreted the data and participated in writing the manuscript. All authors read and approved the manuscript.
We are grateful to all members of the HIV/AIDS & Global Health Research laboratory, Department of Microbiology, University of Venda for their help with sample processing.
The authors declare that they have no competing interests.
Availability of data and materials
Primer and patient sequences are available upon request from the corresponding authors.
Consent for publication
Participants of the study also gave their consent for publication of the study results.
This study was approved by the Health, Safety and Research Ethics Committee of the University of Venda, South Africa (Protocol Number SMNS/13/MBY/01/0625); and permission to access public sector health facilities was obtained from the Limpopo Provincial Department of Health, South Africa. Informed consent was obtained from all study participants prior to demographic data collection and blood draw.
KMG was supported by a stipend from the Division of Global Health, Department of Infectious Diseases, University of Virginia for a Summer Internship Program at the University of Venda, South Africa. DMT was supported by Funds from the Myles H. Thaler Center for AIDS and Human Retrovirus Research. DMT also received partial support through a Carnegie African Diaspora Fellowship Award. EME was supported by the Research and Publications Committee of the University of Venda. POB received grant support from the South African Medical Research Council and the National Research Foundation (GUN 76684). The contents expressed here do not necessarily represent the views of the South African Medical Research Council or the National Research Foundation.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- WHO U, UNAIDS: Global update on HIV treatment 2013: results, impact and opportunities. 2014.Google Scholar
- Tanser F, de Oliveira T, Maheu-Giroux M, Barnighausen T. Concentrated HIV subepidemics in generalized epidemic settings. Curr Opin HIV AIDS. 2014;9(2):115–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Bangsberg DR. Preventing HIV antiretroviral resistance through better monitoring of treatment adherence. J Infect Dis. 2008;197:S272–8.View ArticlePubMedGoogle Scholar
- Gupta RK, Jordan MR, Sultan BJ, Hill A, Davis DH, Gregson J, Sawyer AW, Hamers RL, Ndembi N, Pillay D, 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. Lancet. 2012;380(9849):1250–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Richman DD, Margolis DM, Delaney M, Greene WC, Hazuda D, Pomerantz RJ. The challenge of finding a cure for HIV infection. Science. 2009;323(5919):1304–7.View ArticlePubMedGoogle Scholar
- Peter T, Ellenberger D, Kim AA, Boeras D, Messele T, Roberts T, Stevens W, Jani I, Abimiku A, Ford N, et al. Early antiretroviral therapy initiation: access and equity of viral load testing for HIV treatment monitoring. Lancet Infect Dis. 2017;17(1):e26–9.View ArticlePubMedGoogle Scholar
- Roberts T, Cohn J, Bonner K, Hargreaves S. Scale-up of routine viral load testing in resource-poor settings: current and future implementation challenges. Clin Infect Dis. 2016;62(8):1043–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Goodall RL, Dunn DT, Nkurunziza P, Mugarura L, Pattery T, Munderi P, Kityo C, Gilks C, Kaleebu P, Pillay D, et al. Rapid accumulation of HIV-1 thymidine analogue mutations and phenotypic impact following prolonged viral failure on zidovudine-based first-line ART in sub-Saharan Africa. J Antimicrob Chemother. 2017;72(5):1450–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Boender TS, Kityo CM, Boerma RS, Hamers RL, Ondoa P, Wellington M, Siwale M, Nankya I, Kaudha E, Akanmu AS, et al. Accumulation of HIV-1 drug resistance after continued virological failure on first-line ART in adults and children in sub-Saharan Africa. J Antimicrob Chemother. 2016;71(10):2918–27.View ArticlePubMedGoogle Scholar
- Tebit DM, Sangare L, Makamtse A, Yameogo S, Somlare H, Bado G, Kouldiaty BG, Sathiandee K, Tiba F, Sanou I, et al. HIV drug resistance pattern among HAART-exposed patients with suboptimal virological response in Ouagadougou, Burkina Faso. J Acquir Immune Defic Syndr. 2008;49(1):17–25.View ArticlePubMedGoogle Scholar
- Delaugerre C, Braun J, Charreau I, Delarue S, Nere ML, de Castro N, May T, Marchou B, Simon F, Molina JM, et al. Comparison of resistance mutation patterns in historical plasma HIV RNA genotypes with those in current proviral HIV DNA genotypes among extensively treated patients with suppressed replication. HIV Med. 2012;13(9):517–25.PubMedGoogle Scholar
- Aghokeng AF, Kouanfack C, Eymard-Duvernay S, Butel C, Edoul GE, Laurent C, Koulla-Shiro S, Delaporte E, Mpoudi-Ngole E, Peeters M. Virological outcome and patterns of HIV-1 drug resistance in patients with 36 months’ antiretroviral therapy experience in Cameroon. J Int AIDS Soc. 2013;16:18004.View ArticlePubMedPubMed CentralGoogle Scholar
- Dolling DI, Dunn DT, Geretti AM, Sabin CA, Database UHDR. Study UKCHC: HIV-1 subtype and virological response to antiretroviral therapy: a confirmatory analysis. Clin Infect Dis. 2013;56(1):162–3.View ArticlePubMedGoogle Scholar
- El-Khatib Z, Delong AK, Katzenstein D, Ekstrom AM, Ledwaba J, Mohapi L, Laher F, Petzold M, Morris L, Kantor R. Drug resistance patterns and virus re-suppression among HIV-1 subtype C infected patients receiving non-nucleoside reverse transcriptase inhibitors in South Africa. J AIDS Clin Res. 2011;2(117). doi:10.4172/2155-6113.1000117
- Wallis CL, Aga E, Ribaudo H, Saravanan S, Norton M, Stevens W, Kumarasamy N, Bartlett J, Katzenstein D. Drug susceptibility and resistance mutations after first-line failure in resource limited settings. Clin Infect Dis. 2014;59(5):706–15.View ArticlePubMedPubMed CentralGoogle Scholar
- TenoRes Study G. Global epidemiology of drug resistance after failure of WHO recommended first-line regimens for adult HIV-1 infection: a multicentre retrospective cohort study. Lancet Infect Dis. 2016;16(5):565–75.View ArticleGoogle Scholar
- Brenner B, Turner D, Oliveira M, Moisi D, Detorio M, Carobene M, Marlink RG, Schapiro J, Roger M, Wainberg MA. A V106M mutation in HIV-1 clade C viruses exposed to efavirenz confers cross-resistance to non-nucleoside reverse transcriptase inhibitors. AIDS. 2003;17(1):F1–5.View ArticlePubMedGoogle Scholar
- Li JZ, Paredes R, Ribaudo HJ, Kozal MJ, Svarovskaia ES, Johnson JA, Geretti AM, Metzner KJ, Jakobsen MR, Hullsiek KH, et al. Impact of minority nonnucleoside reverse transcriptase inhibitor resistance mutations on resistance genotype after virologic failure. J Infect Dis. 2013;207(6):893–7.View ArticlePubMedGoogle Scholar
- Brehm JH, Koontz DL, Wallis CL, Shutt KA, Sanne I, Wood R, McIntyre JA, Stevens WS, Sluis-Cremer N, Mellors JW, et al. Frequent emergence of N348I in HIV-1 subtype C reverse transcriptase with failure of initial therapy reduces susceptibility to reverse-transcriptase inhibitors. Clin Infect Dis. 2012;55(5):737–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Nachega JB, Hislop M, Nguyen H, Dowdy DW, Chaisson RE, Regensberg L, Cotton M, Maartens G. Antiretroviral therapy adherence, virologic and immunologic outcomes in adolescents compared with adults in southern Africa. J Acquir Immune Defic Syndr 1999. 2009;51(1):65–71.View ArticleGoogle Scholar
- Tang MW, Shafer RW. HIV-1 antiretroviral resistance: scientific principles and clinical applications. Drugs. 2012;72(9):e1–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Orrell C, Walensky RP, Losina E, Pitt J, Freedberg KA, Wood R. HIV type-1 clade C resistance genotypes in treatment-naive patients and after first virological failure in a large community antiretroviral therapy programme. Antivir Ther. 2009;14(4):523–31.PubMedPubMed CentralGoogle Scholar
- Hoffmann CJ, Charalambous S, Sim J, Ledwaba J, Schwikkard G, Chaisson RE, Fielding KL, Churchyard GJ, Morris L, Grant AD. Viremia, resuppression, and time to resistance in human immunodeficiency virus (HIV) subtype C during first-line antiretroviral therapy in South Africa. Clin Infect Dis. 2009;49(12):1928–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Heckman TG, Somlai AM, Peters J, Walker J, Otto-Salaj L, Galdabini CA, Kelly JA. Barriers to care among persons living with HIV/AIDS in urban and rural areas. AIDS Care. 1998;10(3):365–75.View ArticlePubMedGoogle Scholar
- Wilson DP, Kahn J, Blower SM. Predicting the epidemiological impact of antiretroviral allocation strategies in KwaZulu-Natal: the effect of the urban–rural divide. Proc Natl Acad Sci USA. 2006;103(38):14228–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Rossouw TM, Nieuwoudt M, Manasa J, Malherbe G, Lessells RJ, Pillay S, Danaviah S, Mahasha P, van Dyk G, de Oliveira T. HIV drug resistance levels in adults failing first-line antiretroviral therapy in an urban and a rural setting in South Africa. HIV Med. 2017;18(2):104–14.View ArticlePubMedGoogle Scholar
- Musiime V, Kayiwa J, Kiconco M, Tamale W, Alima H, Mugerwa H, Abwola M, Apilli E, Ahimbisibwe F, Kizito H, et al. Response to antiretroviral therapy of HIV type 1-infected children in urban and rural settings of Uganda. AIDS Res Hum Retrovir. 2012;28(12):1647–57.View ArticlePubMedGoogle Scholar
- Kabamba-Mukadi B, Duquenne A, Henrivaux P, Musuamba F, Ruelle J, Yombi JC, Bodeus M, Vandercam B, Goubau P. HIV-1 proviral resistance mutations: usefulness in clinical practice. HIV Med. 2010;11(8):483–92.PubMedGoogle Scholar
- Gallien S, Charreau I, Nere ML, Mahjoub N, Simon F, de Castro N, Aboulker JP, Molina JM, Delaugerre C. Archived HIV-1 DNA resistance mutations to rilpivirine and etravirine in successfully treated HIV-1-infected individuals pre-exposed to efavirenz or nevirapine. J Antimicrob Chemother. 2015;70(2):562–5.View ArticlePubMedGoogle Scholar
- Bessong PO, Mphahlele J, Choge IA, Obi LC, Morris L, Hammarskjold ML, Rekosh DM. Resistance mutational analysis of HIV type 1 subtype C among rural South African drug-naive patients prior to large-scale availability of antiretrovirals. AIDS Res Hum Retrovir. 2006;22(12):1306–12.View ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Thompson MA, Aberg JA, Hoy JF, Telenti A, Benson C, Cahn P, Eron JJ, Gunthard HF, Hammer SM, Reiss P, et al. Antiretroviral treatment of adult HIV infection: 2012 recommendations of the International Antiviral Society-USA panel. JAMA. 2012;308(4):387–402.View ArticlePubMedGoogle Scholar
- Rimsky L, Vingerhoets J, Van Eygen V, Eron J, Clotet B, Hoogstoel A, Boven K, Picchio G. Genotypic and phenotypic characterization of HIV-1 isolates obtained from patients on rilpivirine therapy experiencing virologic failure in the phase 3 ECHO and THRIVE studies: 48-week analysis. J Acquir Immune Defic Syndr (1999). 2012;59(1):39–46.View ArticleGoogle Scholar
- Kolomeets AN, Varghese V, Lemey P, Bobkova MR, Shafer RW. A uniquely prevalent nonnucleoside reverse transcriptase inhibitor resistance mutation in Russian subtype A HIV-1 viruses. AIDS. 2014;28(17):F1–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Rhee SY, Varghese V, Holmes SP, Van Zyl GU, Steegen K, Boyd MA, Cooper DA, Nsanzimana S, Saravanan S, Charpentier C, et al. Mutational correlates of virological failure in individuals receiving a WHO-recommended tenofovir-containing first-line regimen: an international collaboration. EBioMedicine. 2017;18:225–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Villabona-Arenas CJ, Vidal N, Guichet E, Serrano L, Delaporte E, Gascuel O, Peeters M. In-depth analysis of HIV-1 drug resistance mutations in HIV-infected individuals failing first-line regimens in West and Central Africa. AIDS. 2016;30(17):2577–89.View ArticlePubMedGoogle Scholar
- Brenner BG, Oliveira M, Doualla-Bell F, Moisi DD, Ntemgwa M, Frankel F, Essex M, Wainberg MA. HIV-1 subtype C viruses rapidly develop K65R resistance to tenofovir in cell culture. AIDS. 2006;20(9):F9–13.View ArticlePubMedGoogle Scholar
- Gregson J, Kaleebu P, Marconi VC, van Vuuren C, Ndembi N, Hamers RL, Kanki P, Hoffmann CJ, Lockman S, Pillay D, et al. Occult HIV-1 drug resistance to thymidine analogues following failure of first-line tenofovir combined with a cytosine analogue and nevirapine or efavirenz in sub Saharan Africa: a retrospective multi-centre cohort study. Lancet Infect Dis. 2017;17(3):296–304.View ArticlePubMedGoogle Scholar
- Novitsky V, Wester CW, DeGruttola V, Bussmann H, Gaseitsiwe S, Thomas A, Moyo S, Musonda R, Van WE, Marlink RG, 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 Res Hum Retrovir. 2007;23(7):868–78.View ArticlePubMedGoogle Scholar
- Sarmati L, Nicastri E, Parisi SG, D’Ettorre G, Narciso P, Mancino G, Gallo I, Abbadessa V, Dalle NE, Traina C, et al. Failure of stavudine–lamivudine combination therapy in antiretroviral-naive patients with AZT-like HIV-1 resistance mutations. J Med Virol. 2001;65(4):631–6.View ArticlePubMedGoogle Scholar
- Boerma RS, Sigaloff KC, Akanmu AS, Inzaule S, van Boele Hensbroek M, de Rinke Wit TF, Calis JC. Alarming increase in pretreatment HIV drug resistance in children living in sub-Saharan Africa: a systematic review and meta-analysis. J Antimicrob Chemother. 2017;72(2):365–71.View ArticlePubMedGoogle Scholar
- Wallis CL, Papathanasopolous MA, Fox M, Conradie F, Ive P, Orrell C, Zeinecker J, Sanne I, Wood R, McIntyre J, et al. Low rates of nucleoside reverse transcriptase inhibitor resistance in a well-monitored cohort in South Africa on antiretroviral therapy. Antivir Ther. 2012;17(2):313–20.View ArticlePubMedGoogle Scholar
- Hamers RL, Sigaloff KC, Kityo C, Mugyenyi P, de Wit TF. HIV-1 drug resistance in antiretroviral-naive patients in sub-Saharan Africa. Lancet Infect Dis. 2013;13(3):196–7.View ArticlePubMedGoogle Scholar
- Hamers RL, Derdelinckx I, van VM, Stevens W, de Wit TFR, Schuurman R. The status of HIV-1 resistance to antiretroviral drugs in sub-Saharan Africa. Antivir Ther. 2008;13(5):625–39 [Erratum in: Antivir Ther. 2008;13(8):1129].PubMedGoogle Scholar
- Rhee SY, Gonzales MJ, Kantor R, Betts BJ, Ravela J, Shafer RW. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. 2003;31(1):298–303.View ArticlePubMedPubMed CentralGoogle Scholar
- Gibson RM, Meyer AM, Winner D, Archer J, Feyertag F, Ruiz-Mateos E, Leal M, Robertson DL, Schmotzer CL, Quinones-Mateu ME. 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(4):2167–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Kyeyune F, Gibson RM, Nankya I, Venner C, Metha S, Akao J, Ndashimye E, Kityo CM, Salata RA, Mugyenyi P, et al. Low frequency drug resistance in HIV-infected Ugandans on antiretroviral treatment is associated with regimen failure. Antimicrob Agents Chemother. 2016;60(6):3380–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Palmer S, Boltz V, Martinson N, Maldarelli F, Gray G, McIntyre J, Mellors J, Morris L, Coffin J. Persistence of nevirapine-resistant HIV-1 in women after single-dose nevirapine therapy for prevention of maternal-to-fetal HIV-1 transmission. Proc Natl Acad Sci USA. 2006;103(18):7094–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Kimata JT, Rice AP, Wang J. Challenges and strategies for the eradication of the HIV reservoir. Curr Opin Immunol. 2016;42:65–70.View ArticlePubMedGoogle Scholar
- Devereux HL, Loveday C, Youle M, Sabin CA, Burke A, Johnson MA. Reduction in human immunodeficiency virus type 1 mutations associated with drug resistance after initiating new therapeutic regimens in pretreated patients. J Infect Dis. 2000;181(5):1804–7.View ArticlePubMedGoogle Scholar
- Turriziani O, Boni A, Falasca F, Graziano F, Bucci M, D’Ettorre G, Fantauzzi A, Paoletti F, Massetti AP, Mezzaroma I, et al. Mutational resistance pattern of HIV type 1 in CD14+ monocytes, CD4+ T cells, and plasma from treated patients. AIDS Res Hum Retrovir. 2010;26(6):625–34.View ArticlePubMedGoogle Scholar