Background

Routine HIV-1 antiretroviral drug resistance testing for patients failing NNRTI-based regimens is not recommended in resource-limited settings. Therefore, surveys are required to monitor resistance profiles in patients failing ART.

Methods

A cross-sectional survey was conducted amongst patients failing NNRTI-based regimens in the public sector throughout South Africa. Virological failure was defined as two consecutive HIV-1 viral load results >1000 RNA copies/mL. Pol sequences were obtained using RT–PCR and Sanger sequencing and submitted to Stanford HIVdb v7.0.1.

Results

A total of 788 sequences were available for analysis. Most patients failed a tenofovir-based NRTI backbone (74.4%) in combination with efavirenz (82.1%) after median treatment duration of 36 months. K103N (48.9%) and V106M (34.9%) were the most common NNRTI mutations. Only one-third of patients retained full susceptibility to second-generation NNRTIs such as etravirine (36.5%) and rilpivirine (27.3%). After M184V/I (82.7%), K65R was the most common NRTI mutation (45.8%). The prevalence of K65R increased to 57.5% in patients failing a tenofovir regimen without prior stavudine exposure. Cross-resistance to NRTIs was often observed, but did not seem to affect the predicted activity of zidovudine as 82.9% of patients remained fully susceptible to this drug.

Conclusions

The introduction of tenofovir-based first-line regimens has dramatically increased the prevalence of K65R mutations in the HIV-1-infected South African population. However, most patients failing tenofovir-based regimens remained fully susceptible to zidovudine. Based on these data, there is currently no need to change either the recommended first- or second-line ART regimens in South Africa.

Introduction

It is well established that the introduction of ART has greatly reduced morbidity and mortality in HIV-1-infected patients. In sub-Saharan Africa, the number of AIDS-related deaths fell by 42% between 2004 and 2014.1 By July 2015, 3.1 million South Africans were on ART.2 Over the years the South African national treatment guidelines have changed: prior to 2010, a first-line regimen consisting of stavudine or zidovudine in combination with lamivudine and nevirapine or efavirenz was recommended with a preferred switch to zidovudine, didanosine and ritonavir-boosted lopinavir for second-line treatment,3 whereas the 2010 guidelines advised starting patients on a tenofovir-based regimen, with a preferred switch to zidovudine, lamivudine and ritonavir-boosted lopinavir or tenofovir, lamivudine and ritonavir-boosted lopinavir for those patients failing a stavudine- or zidovudine-based first-line regimen.4 The latest guidelines emphasize the use of fixed-dose combinations in first-line treatment.5,6 These changes are expected to have an effect over time on the observed HIV-1 drug resistance profiles.

With increasing ART coverage and longer treatment, the proportion of patients experiencing ART failure, with associated drug resistance will increase and might jeopardize the long-term success of ART programmes. WHO recommends an HIV drug resistance surveillance strategy for resource-constrained settings, including a combination of monitoring early warning indicators of poor programme performance, surveillance of transmitted drug resistance and evaluation of acquired HIV-1 drug resistance.7

In South Africa drug resistance surveillance activities are conducted in combination with limited HIV-1 drug resistance testing for diagnostic use5,6 and research activities, which makes for a diverse dataset. A total of 18 studies, conducted between 2000 and 2012, have described the HIV drug resistance patterns in South Africans failing first-line treatment regimens. Six studies exclusively described HIV-1 drug resistance profiles of patients failing stavudine- or zidovudine-based regimens.8–13 In six additional studies, most patients failed stavudine- or zidovudine-based regimens, but some patients failed other regimens such as tenofovir- or didanosine-based regimens.14–19 Four studies addressed the development of mutations after exposure to tenofovir-based regimens.20–23 The final two studies did not specify the details of the ART regimens used.24,25 The majority of studies reported the M184V/I mutation as the most prevalent NRTI mutation. The presence of at least one thymidine analogue mutation (TAM) varied substantially across the various studies, ranging from none to 55%; this variation is likely due to the different treatment regimens, duration of treatment and frequency of viral load testing. The introduction of tenofovir in first-line regimens has significantly increased the frequency of the cross-resistance mutation K65R.20–23 Skhosana et al.21 reported that the use of tenofovir increased the risk of developing K65R by almost 5-fold, compared with the use of stavudine. For NNRTI mutations, there is general consensus in these reports that K103N was present in more than half of the patients. The prevalence of Y181C and V106M/A was largely dependent on the use of nevirapine or efavirenz in each population.11,12,21,23

The national treatment guidelines only recommend HIV drug resistance testing for patients failing second-line treatment.5,6 In the absence of routine HIV-1 antiretroviral drug resistance testing in first-line failures, we embarked on the first national cross-sectional survey of patients failing NNRTI-based treatment to provide information on the current resistance patterns within the South African population.

Methods

Survey design

A cross-sectional survey was designed using probability proportional to size sampling. This method ensures that samples were proportionally collected from each of the nine South African provinces, based on the number of patients receiving ART in each region. Only tier 2 and 3 facilities with at least 500 patients on ART were considered, to ensure adequate sampling options. At the time of the study, these facilities treated ∼85% of patients on ART. The number of patients on ART was obtained from the South African District Health Information System. Any patient ≥18 years old failing an NNRTI-based regimen who gave signed informed consent was eligible to be included in the study. Virological failure was defined as two consecutive HIV-1 VL (HIVVL) measurements >1000 RNA copies/mL after being exposed to an NNRTI-based regimen for at least 6 months.

Data and sample collection

Healthcare professionals at selected sentinel sites informed eligible patients about the survey, obtained written informed consent and completed a study questionnaire. Samples were collected between February 2013 and October 2014 from 91 healthcare facilities across 37 districts and nine provinces (Table S1 and Figure S1, available as Supplementary data at JAC Online). The samples were then sent to one of three participating laboratories: Tygerberg Hospital, Stellenbosch; Universitas, Bloemfontein; or Charlotte Maxeke Johannesburg Academic Hospital, Johannesburg. Upon receipt, samples were centrifuged and plasma was stored at −80°C until processing. Each laboratory used their own validated population-based in-house genotyping method to obtain pol sequences.26–28 The laboratories at Tygerberg Hospital and Charlotte Maxeke Johannesburg Academic Hospital have successfully participated in at least one international external quality assurance panel per year for at least 5 consecutive years. An inter-laboratory comparison panel was run between the Universitas laboratory and the Charlotte Maxeke Johannesburg Academic Hospital laboratory to ensure all results were comparable, regardless of the testing laboratory.

Sequence analysis

All obtained nucleotide sequences were collated at the Charlotte Maxeke Johannesburg Academic Hospital for analysis. The Stanford HIVdb v7.0.1 tool (http://sierra2.stanford.edu/sierra/servlet/JSierra) was used to identify drug resistance mutations and generate predicted resistance profiles. Predicted antiretroviral drug resistance profiles were categorized as being susceptible (including potential low-level resistance), intermediate resistance (including low-level resistance) or high-level drug resistance.

Sequences were aligned using Muscle29 and PhyML v 3.130 and FastTree31,32 were used to construct a phylogenetic tree to ensure that no duplicate patient samples were included in the analyses.

Subtyping

A multiple alignment of the 788 nucleotide sequences with the 2010 subtype reference sequences from HIV-1 subtypes A–K, and several Circulating Recombinant Forms (http://www.hiv.lanl.gov/), were generated using Clustal X (version 2.1) (http://www.clustal.org) to determine subtypes The generated Clustal X alignment was used to construct a Neighbour-Joining phylogenetic tree in MEGA 7 (http://www.megasoftware.net) with the Kimura two-parameter model. The stability of the nodes was assessed by bootstrap analysis (1000 replicates), and bootstrap values >70% were considered significant.

In addition, pol subtyping was performed using Rega HIV subtyping tool v2.0 (http://bioafrica.net/rega-genotype/html/subtypinghiv.html). Any sequence that was assigned to be a non-C subtype was also analysed using the NCBI genotyping tool http://www.ncbi.nlm.nih.gov/projects/genotyping/ and the jumping profile HMM tool http://jphmm.gobics.de/submission_hiv.html to confirm the subtype.

Statistical analysis

Statistical analyses were performed using GraphPad Prism 6 and SAS v9.3. Differences between groups for non-parametric data were calculated using the two-sided Wilcoxon–Mann–Whitney test. Two-sided Fisher's exact tests were performed to assess proportional differences between groups for gender and mutation prevalence. P < 0.05 was considered statistically significant. The 95% CIs were calculated using the modified Wald method. Log-binomial regression analysis was used to test the association between different mutations and antiretroviral drug pressure. The relative risks (RRs) and corresponding 95% CIs are presented. Mutations identified as being significantly associated with drug pressure in the univariate analysis (P < 0.2) and a priori variables of importance were included in multivariate models. Models were adjusted (aRR) for CD4 count, HIVVL, total time on ART and time on current ART where possible, which was dependent on sample sizes.

Sequence data

The pol nucleotide sequences were submitted to GenBank using Sequin v9.50 (www.ncbi.nlm.nih.gov/Sequin) and are available under accession numbers KU127587–KU128374.

Ethics statement

This study was conducted with ethics clearance by the Research on Human Subjects (Medical) Committee at the University of the Witwatersrand (Clearance Number M120254). Written informed consent was obtained from all participants prior to sample collection.

Results

A total of 1033 samples were collected from eight provinces (no samples could be obtained from the Northern Cape province). Only 788 samples were included for analysis. The main reason for exclusion was the inability of obtaining a PCR product (n = 151, 14.6%). Other reasons for exclusion from the study were: volunteers <18 years of age (n = 35, 3.4%); <6 months of ART exposure (n = 16, 1.5%); treatment with a non-standard regimen at time of failure (n = 11, 1.1%); inadequate completion of informed consent (n = 10, 1.0%); poor sample integrity (n = 8, 0.8%); <6 months of NNRTI exposure (n = 6, 0.6%); or failure to obtain a sequence (n = 3, 0.3%). In addition, five sequences (0.5%) were shown to be duplicate samples as per phylogenetic analysis. The quality of samples for which no PCR product could be obtained may have been compromised by prolonged transport time or, more likely, had a low HIVVL. The median HIVVL was significantly lower for the samples that failed PCR (4.0 log copies/mL) versus those that could be reported (4.7 log copies/mL, P < 0.0001). The median time between the latest HIVVL and the date of sample collection for resistance testing was significantly higher for the samples for which no result could be obtained (58 days) versus those that could be reported (46 days, P = 0.0059). Due to budgetary reasons, virological failure at time of sample collection could not be confirmed in the laboratory prior to HIV drug resistance testing.

Most patients were female (n = 511, 64.8%) and the median age was 36 years. Patients had been exposed to ART for a median duration of 36 months (IQR: 19–59). The latest HIVVL was taken 47 days (median, IQR: 28–84) prior to the study sample, with a median HIVVL of 4.7 log RNA copies/mL (IQR: 4.2–5.2). A summary of other baseline clinical characteristics is presented in Table 1.

Table 1.

Demographic and baseline clinical characteristics for patients failing NNRTI-based regimens throughout South Africa

VariableNumber of observationsPercentage of observationsMedianIQRMinimum–maximum
Gender (female)51164.8
Age788100.03631–4218–89
Time between diagnosis and study sample (months)67785.95228–827–277
Baseline CD4 (cells/mm3)64481.79334–1720–1111
Latest CD4 (cells/mm3)72892.417370–3051–1020
Time between latest CD4 and study sample (days)72592.010635–2100–3005
Latest HIVVL (log copies/mL)76296.74.74.2–5.23.0–8.0
Time between latest HIVVL and study sample (days)75796.14728–840–540
Previous HIVVL (log copies/mL)65883.54.53.8–5.01.3–6.7
Time between previous HIVVL and latest HIVVL (days)63480.5163106–30014–2818
Time on current regimen76497.02010–371–168
Total time on NNRTI regimen (months)75495.73519–587–168
Total time on ART (months)76196.63619–597–168
NRTI backbone
 TDF + 3TC41252.3
 TDF + FTC17422.1
 d4T + 3TC11614.7
 ZDV + 3TC557.0
 ABC + 3TC222.8
 other91.1
NNRTI backbone
 EFV64782.1
 NVP14117.9
VariableNumber of observationsPercentage of observationsMedianIQRMinimum–maximum
Gender (female)51164.8
Age788100.03631–4218–89
Time between diagnosis and study sample (months)67785.95228–827–277
Baseline CD4 (cells/mm3)64481.79334–1720–1111
Latest CD4 (cells/mm3)72892.417370–3051–1020
Time between latest CD4 and study sample (days)72592.010635–2100–3005
Latest HIVVL (log copies/mL)76296.74.74.2–5.23.0–8.0
Time between latest HIVVL and study sample (days)75796.14728–840–540
Previous HIVVL (log copies/mL)65883.54.53.8–5.01.3–6.7
Time between previous HIVVL and latest HIVVL (days)63480.5163106–30014–2818
Time on current regimen76497.02010–371–168
Total time on NNRTI regimen (months)75495.73519–587–168
Total time on ART (months)76196.63619–597–168
NRTI backbone
 TDF + 3TC41252.3
 TDF + FTC17422.1
 d4T + 3TC11614.7
 ZDV + 3TC557.0
 ABC + 3TC222.8
 other91.1
NNRTI backbone
 EFV64782.1
 NVP14117.9

HIVVL, viral load; TDF, tenofovir; 3TC, lamivudine; FTC, emtricitabine; d4T, stavudine; ZDV, zidovudine; ABC, abacavir; EFV, efavirenz; NVP, nevirapine.

Table 1.

Demographic and baseline clinical characteristics for patients failing NNRTI-based regimens throughout South Africa

VariableNumber of observationsPercentage of observationsMedianIQRMinimum–maximum
Gender (female)51164.8
Age788100.03631–4218–89
Time between diagnosis and study sample (months)67785.95228–827–277
Baseline CD4 (cells/mm3)64481.79334–1720–1111
Latest CD4 (cells/mm3)72892.417370–3051–1020
Time between latest CD4 and study sample (days)72592.010635–2100–3005
Latest HIVVL (log copies/mL)76296.74.74.2–5.23.0–8.0
Time between latest HIVVL and study sample (days)75796.14728–840–540
Previous HIVVL (log copies/mL)65883.54.53.8–5.01.3–6.7
Time between previous HIVVL and latest HIVVL (days)63480.5163106–30014–2818
Time on current regimen76497.02010–371–168
Total time on NNRTI regimen (months)75495.73519–587–168
Total time on ART (months)76196.63619–597–168
NRTI backbone
 TDF + 3TC41252.3
 TDF + FTC17422.1
 d4T + 3TC11614.7
 ZDV + 3TC557.0
 ABC + 3TC222.8
 other91.1
NNRTI backbone
 EFV64782.1
 NVP14117.9
VariableNumber of observationsPercentage of observationsMedianIQRMinimum–maximum
Gender (female)51164.8
Age788100.03631–4218–89
Time between diagnosis and study sample (months)67785.95228–827–277
Baseline CD4 (cells/mm3)64481.79334–1720–1111
Latest CD4 (cells/mm3)72892.417370–3051–1020
Time between latest CD4 and study sample (days)72592.010635–2100–3005
Latest HIVVL (log copies/mL)76296.74.74.2–5.23.0–8.0
Time between latest HIVVL and study sample (days)75796.14728–840–540
Previous HIVVL (log copies/mL)65883.54.53.8–5.01.3–6.7
Time between previous HIVVL and latest HIVVL (days)63480.5163106–30014–2818
Time on current regimen76497.02010–371–168
Total time on NNRTI regimen (months)75495.73519–587–168
Total time on ART (months)76196.63619–597–168
NRTI backbone
 TDF + 3TC41252.3
 TDF + FTC17422.1
 d4T + 3TC11614.7
 ZDV + 3TC557.0
 ABC + 3TC222.8
 other91.1
NNRTI backbone
 EFV64782.1
 NVP14117.9

HIVVL, viral load; TDF, tenofovir; 3TC, lamivudine; FTC, emtricitabine; d4T, stavudine; ZDV, zidovudine; ABC, abacavir; EFV, efavirenz; NVP, nevirapine.

The implementation of the various ART guidelines is reflected in the use of both the NRTI backbone and the choice of NNRTI (Table 1). Over half of the patients were treated with a tenofovir and lamivudine backbone and an additional 22.1% of patients received a tenofovir and emtricitabine backbone at the time of resistance testing. The proportion of emtricitabine use might be under-represented as some healthcare workers consider emtricitabine and lamivudine to be the same. Overall, 304 patients (38.6%) had been exposed to stavudine at some stage during their ART; 118 of these patients (15.0%) were still receiving a stavudine regimen at the time of sampling. When analysing the NRTI backbone data, a distinction was made between those currently failing a tenofovir-based regimen without prior exposure to stavudine (‘tenofovir-only group’, n = 424); those failing a tenofovir-based regimen with recorded prior exposure to stavudine (‘tenofovir, prior stavudine group’, n = 162); those failing a stavudine-based regimen (‘stavudine group’, n = 118); patients failing a zidovudine-based regimen (‘zidovudine group’, n = 55) and finally patients failing any other regimens (‘other group’, n = 29). Patients failing a tenofovir-only regimen had a significantly higher median HIVVL and lower CD4 count compared with stavudine-treated patients, zidovudine-treated patients and those failing tenofovir with known prior exposure to stavudine (Table 2).

Table 2.

Comparison of CD4 count, HIVVL measurements and time on ART between different treatment groups for patients failing an NNRTI-based ART regimen throughout South Africa

Treatment group
CD4 count (cells/mm3)
HIVVL (log copies/mL)
Total time on ART (months)
Reference groupmedianPmedianPmedianP
TDF only (n = 424)1414.825
TDF prior d4T1880.00224.70.012662<0.0001
d4T210<0.00014.5<0.000152<0.0001
ZDV249<0.00014.50.010453<0.0001
other171NS4.7NS29NS
TDF prior d4T (n = 162)1884.762 
d4T210NS4.5NS520.0077
ZDV249NS4.5NS53NS
other171NS4.7NS29<0.0001
d4T (n = 118)210 4.5 52 
ZDV249NS4.5NS53NS
other171NS4.7NS290.0004
ZDV (n = 55)249 4.5 53 
other (n = 29)1710.02364.7NS290.0093
EFV (n = 647)163 4.8 34 
NVP (n = 141)2390.00104.50.0012450.0004
Treatment group
CD4 count (cells/mm3)
HIVVL (log copies/mL)
Total time on ART (months)
Reference groupmedianPmedianPmedianP
TDF only (n = 424)1414.825
TDF prior d4T1880.00224.70.012662<0.0001
d4T210<0.00014.5<0.000152<0.0001
ZDV249<0.00014.50.010453<0.0001
other171NS4.7NS29NS
TDF prior d4T (n = 162)1884.762 
d4T210NS4.5NS520.0077
ZDV249NS4.5NS53NS
other171NS4.7NS29<0.0001
d4T (n = 118)210 4.5 52 
ZDV249NS4.5NS53NS
other171NS4.7NS290.0004
ZDV (n = 55)249 4.5 53 
other (n = 29)1710.02364.7NS290.0093
EFV (n = 647)163 4.8 34 
NVP (n = 141)2390.00104.50.0012450.0004

TDF, tenofovir; d4T, stavudine; ZDV, zidovudine; EFV, efavirenz; NVP, nevirapine; NS, not significant.

Table 2.

Comparison of CD4 count, HIVVL measurements and time on ART between different treatment groups for patients failing an NNRTI-based ART regimen throughout South Africa

Treatment group
CD4 count (cells/mm3)
HIVVL (log copies/mL)
Total time on ART (months)
Reference groupmedianPmedianPmedianP
TDF only (n = 424)1414.825
TDF prior d4T1880.00224.70.012662<0.0001
d4T210<0.00014.5<0.000152<0.0001
ZDV249<0.00014.50.010453<0.0001
other171NS4.7NS29NS
TDF prior d4T (n = 162)1884.762 
d4T210NS4.5NS520.0077
ZDV249NS4.5NS53NS
other171NS4.7NS29<0.0001
d4T (n = 118)210 4.5 52 
ZDV249NS4.5NS53NS
other171NS4.7NS290.0004
ZDV (n = 55)249 4.5 53 
other (n = 29)1710.02364.7NS290.0093
EFV (n = 647)163 4.8 34 
NVP (n = 141)2390.00104.50.0012450.0004
Treatment group
CD4 count (cells/mm3)
HIVVL (log copies/mL)
Total time on ART (months)
Reference groupmedianPmedianPmedianP
TDF only (n = 424)1414.825
TDF prior d4T1880.00224.70.012662<0.0001
d4T210<0.00014.5<0.000152<0.0001
ZDV249<0.00014.50.010453<0.0001
other171NS4.7NS29NS
TDF prior d4T (n = 162)1884.762 
d4T210NS4.5NS520.0077
ZDV249NS4.5NS53NS
other171NS4.7NS29<0.0001
d4T (n = 118)210 4.5 52 
ZDV249NS4.5NS53NS
other171NS4.7NS290.0004
ZDV (n = 55)249 4.5 53 
other (n = 29)1710.02364.7NS290.0093
EFV (n = 647)163 4.8 34 
NVP (n = 141)2390.00104.50.0012450.0004

TDF, tenofovir; d4T, stavudine; ZDV, zidovudine; EFV, efavirenz; NVP, nevirapine; NS, not significant.

As per national guidelines, efavirenz is the preferred NNRTI with 82.1% of the patients receiving efavirenz at the time of sample collection (Table 1). Patients failing an efavirenz-based regimen presented with a higher HIVVL and lower CD4 count compared with those failing a nevirapine-based regimen (Table 2). For further analysis the treatment groups were split into efavirenz-based (n = 647) versus nevirapine-based (n = 141).

Significant differences were observed for the total time on ART between the various groups. Patients failing a tenofovir-based regimen without prior exposure to stavudine had been treated for an average of 25 months versus 62, 52 and 53 months for tenofovir with prior exposure to stavudine, stavudine and zidovudine-treated patients, respectively (P < 0.0001), which is reflective of the earlier introduction of stavudine and zidovudine in the national treatment guidelines. Likewise, patients treated with efavirenz at the time of sample collection had been treated for a shorter time compared with those that were still on a nevirapine-based regimen (34 versus 45 months; Table 2).

Sixteen patients (2.0%) had exposure to a PI during their treatment history; however, they were receiving an NNRTI-based regimen at the time of drug resistance testing.

Most patients (n = 773, 98.1%) were infected with HIV-1 subtype C, which is the predominant subtype circulating in South Africa. The remaining 15 patients were infected with a C–D recombinant (n = 3), subtype A1 (n = 2), subtype B (n = 2), A1–C recombinant (n = 2), CRF01_AE (n = 1), CRF02_AG (n = 1), subtype G (n = 1) and A1–D recombinant (n = 1) virus. For two patients, the pol subtype could not be determined.

Completely susceptible resistance profiles across the three drug classes (PI, NRTI and NNRTI) were observed in only 31 of 788 patients (3.9%, 95% CI = 2.8%–5.6%). This included nine patients who presented with intermediate resistance caused by the presence of one or more polymorphisms (T74S, Q58E and E138A), which are common in HIV-1 subtype C, without the need for drug pressure.

Most patients presented with NNRTI resistance (752 of 788, 95.4%, 95% CI = 93.7%–96.7%). Overall, K103N was the most commonly observed NNRTI mutation (n = 385, 48.9%), followed by V106M (n = 275, 34.9%), Y181C (n = 206, 26.1%) and G190A (n = 169, 21.4%). The prevalence of these and other NNRTI mutations is presented in Figure 1. Detailed analysis of the resistance mutations in relation to the ART regimens revealed that the K103N mutation was more commonly observed in patients failing an efavirenz-based regimen n = 342, 52.9% versus n = 44, 30.5%, aRR = 0.51, 95% CI = 0.38–0.68, P < 0.0001). Similarly, V106M was selected twice as often under efavirenz pressure (n = 247, 38.2% versus n = 28, 19.9%, aRR = 0.57, 95% CI = 0.39–0.82, P = 0.0025). In contrast the Y181C mutation was found in 54.6% (n = 77) and 20.0% (n = 129) of patients exposed to nevirapine and efavirenz, respectively (aRR = 2.4, 95% CI = 1.99–2.89, P < 0.0001). P225H is known to be selected under efavirenz pressure, which is reflected in the difference in prevalence between the two groups: 98 patients (15.1%) in the efavirenz group compared with only three patients (2.1%) in the nevirapine group (aRR = 0.15, 95% CI = 0.05–0.47, P = 0.0012).
Prevalence of NNRTI mutations observed in 788 South African patients failing efavirenz (EFV)- or nevirapine (NVP)-based NNRTI regimens.
Figure 1.

Prevalence of NNRTI mutations observed in 788 South African patients failing efavirenz (EFV)- or nevirapine (NVP)-based NNRTI regimens.

More than 90% of the patients failing NNRTI-based regimens have lost any susceptibility to the first-generation NNRTI drugs, nevirapine and efavirenz. Only one-third of these patients remain fully susceptible to the second-generation NNRTIs such as etravirine (36.5%) and rilpivirine (27.3%; Figure S2). Patients exposed to efavirenz at the time of resistance testing were twice as likely to remain fully susceptible to etravirine (n = 258, 39.9%) versus those failing a nevirapine-based regimen (n = 29, 20.6%, aRR = 0.48, 95% CI = 0.35–0.66, P < 0.0001). Similar differences were observed for rilpivirine where 194 patients (30.0%) failing efavirenz-based regimens remained fully susceptible to rilpivirine versus only 20 (14.2%) nevirapine-exposed patients (aRR = 0.42, 95% CI = 0.27–0.66, P < 0.0001).

Most patients presented with at least one NRTI mutation (719 of 788, 91.2%, 95% CI = 89.1%–93.0%). The complete NRTI mutation profile is presented in Figure 2. The M184V/I mutation was detected in >80% of patients (n = 652, 82.7%). The next most prevalent mutation was K65R (n = 361, 45.8%). The prevalence of K65R in patients failing a tenofovir-based regimen, with or without prior exposure to stavudine, was 51.2% and 57.5%, respectively. The K65R mutation was significantly more common in the tenofovir-only group (n = 244, 57.5%) compared with the stavudine group (n = 18, 15.3%, aRR = 0.31, 95% CI = 0.19–0.50, P < 0.0001) and the zidovudine group (n = 6, 10.9%, aRR = 0.23, 95% CI = 0.11–0.49, P < 0.0002). The Y115F mutation, which is commonly selected under abacavir pressure, was most commonly observed in the other group (27.6%), which is not surprising as 23 of 29 patients in this group failed an abacavir-based regimen. However, this mutation was also present in 24.5% of the patients in the tenofovir-only group. The L74IV mutation, commonly selected under abacavir pressure, was often detected in the other group (34.5%). At least one in five patients presented with ≥1 TAM, regardless of the treatment group, with the zidovudine group showing the highest prevalence of ≥1 TAM (43.6%) and ≥3 TAMs (20.0%). Patients in the zidovudine group were more likely to present with ≥3 TAMs (n = 11, 20%) compared with those patients in the tenofovir-only group (n = 12, 2.8%, aRR = 6.25, 95% CI = 2.74–14.3, P < 0.0001). The Q151M mutation was rarely observed across the total study population (n = 12, 1.5%).
Prevalence of NRTI mutations observed in 788 South African patients failing NNRTI-based regimens. TDF, tenofovir; d4T, stavudine; ZDV, zidovudine.
Figure 2.

Prevalence of NRTI mutations observed in 788 South African patients failing NNRTI-based regimens. TDF, tenofovir; d4T, stavudine; ZDV, zidovudine.

The above-mentioned mutations lead to predicted resistance profiles (Figure 3). Tenofovir remained fully active in 34.8% of all participants, but only one in four patients failing a tenofovir-based regimen remained fully susceptible to tenofovir. Cross-resistance to most NRTIs was frequently observed due to the high prevalence of K65R and relatively frequent detection of TAMs. However, in most cases, zidovudine is not affected by this cross-resistance as is demonstrated by 82.9% of patients who remain fully susceptible to zidovudine. This number increased to 92.2% in patients failing a tenofovir-based regimen, without prior exposure to stavudine. However, only 1 in 10 patients remained fully susceptible to abacavir, which is a concern for those patients who cannot be switched to zidovudine due to anaemia.
NRTI drug resistance profiles in 788 South African patients failing NNRTI-based regimens. Group 1: all patients failing an NNRTI-based regimen (n = 788). Group 2: patients failing a tenofovir-based regimen (n = 424). Group 3: patients failing a tenofovir-based regimen with prior exposure to stavudine (n = 162). Group 4: patients failing a stavudine-based regimen (n = 118). Group 5: patients failing a zidovudine-based regimen (n = 55). Group 6: patients failing other regimens (n = 29). Resistance calls: susceptible (S; white), intermediate resistance (IR; grey) and high-level resistance (HLR; black). 3TC, lamivudine; FTC, emtricitabine; ABC, abacavir; ZDV, zidovudine; d4T, stavudine; ddI, didanosine; TDF, tenofovir.
Figure 3.

NRTI drug resistance profiles in 788 South African patients failing NNRTI-based regimens. Group 1: all patients failing an NNRTI-based regimen (n = 788). Group 2: patients failing a tenofovir-based regimen (n = 424). Group 3: patients failing a tenofovir-based regimen with prior exposure to stavudine (n = 162). Group 4: patients failing a stavudine-based regimen (n = 118). Group 5: patients failing a zidovudine-based regimen (n = 55). Group 6: patients failing other regimens (n = 29). Resistance calls: susceptible (S; white), intermediate resistance (IR; grey) and high-level resistance (HLR; black). 3TC, lamivudine; FTC, emtricitabine; ABC, abacavir; ZDV, zidovudine; d4T, stavudine; ddI, didanosine; TDF, tenofovir.

One in five patients (n = 150, 19.0%) presented with the T74S PI polymorphism, causing intermediate resistance to ritonavir-boosted nelfinavir. The polymorphisms Q58E and L23I were detected in seven and one patient, causing intermediate resistance to ritonavir-boosted tipranavir and ritonavir-boosted nelfinavir, respectively. Six patients (0.8%) presented with at least one major PI mutation; however, only one of these patients was amongst the 16 who had recorded exposure to a ritonavir-boosted lopinavir-based regimen. Five of the six patients presented with a single major PI mutation: M46L (n = 2), I54V (n = 1), V82A (n = 1) and L90M (n = 1). The last patient presented with a combination of three major PI mutations (I54V, V82A and L90M).

Discussion

The first national survey amongst South African patients failing an NNRTI-based regimen in the public sector, which was conducted between February 2013 and October 2014, confirmed good adherence to the current treatment guidelines, as more than half of the patients were failing a regimen based on a tenofovir and lamivudine backbone, which has been the recommendation since 2010.4 In addition, almost one in four patients were failing the fixed-dose combination (tenofovir, emtricitabine and efavirenz), which has been implemented in the public sector since 2013.5 Adherence to the national treatment guidelines was also reflected in the common use of efavirenz (82.1%).

The frequency of completely drug susceptible virus (3.9%) was shown to be substantially lower in this survey compared with previously published studies from South Africa where the prevalence of susceptible virus ranged from 4.7% to 66.2%.8–20,24 This might indicate that the level of adherence to ART among patients has improved over time as fully drug susceptible virus is only expected in patients that have an extremely low adherence to their drug regimens.33 Alternatively, healthcare workers have refined the selection criteria to identify true treatment failures.

K103N remains the most common NNRTI mutation, previous studies have shown that Y181C is more common in patients exposed to nevirapine, whereas patients exposed to efavirenz more often present with V106M/A;11,12,21,23 these findings have been confirmed in this survey.

The M184V/I mutation remains the most common NRTI mutation, which is known to develop rapidly after lamivudine or emtricitabine exposure; however, lamivudine is recycled in second-line regimens to maintain this mutation due to the negative effect on viral fitness.34–36 The second most common NRTI mutation was K65R, this mutation is known to be rapidly selected under tenofovir pressure in HIV-1 subtype C.37–40 Although this mutation causes cross-resistance to most NRTIs, it is also known to have a crippling effect on viral fitness, comparable to M184V.36,41 Overall, at least one in five patients presented with ≥1 TAM, although the frequency remained highest in patients failing zidovudine-based regimens. The Q151M mutation was rarely observed and was only seen in patients failing stavudine- and zidovudine-based regimens, or patients that had previously been exposed to stavudine. The Y115F mutation, which is commonly selected under abacavir pressure,42 was most commonly observed in the other group (27.6%), which is not surprising as 23 of 29 patients in this group failed an abacavir-based regimen. However, this mutation was also present in 24.5% of the patients in the tenofovir-only group. In vitro experiments have shown that this mutation is often selected after the selection of K65R.43 Similarly, in the population surveyed, 84.8% of patients who presented with Y115F in the tenofovir-only group also presented with the K65R mutation.

Due to the relatively high prevalence of K65R and TAMs, most patients presented with NRTI cross-resistance patterns; however, 82.9% of these patients remained fully susceptible to zidovudine. The frequency of zidovudine susceptibility increased to 92.2% in patients who were failing a tenofovir-based regimen without prior exposure to stavudine. This finding indicates the currently available second-line regimens are still suitable for most patients failing an NNRTI-based regimen. There is, however, a concern in patients for whom zidovudine might be contra-indicated due to anaemia, as only 10% of patients remained fully susceptible to abacavir. However, a recent randomized controlled trial has shown that a PI-based second-line regimen in combination with an NRTI backbone was superior to PI monotherapy and non-inferior to ritonavir-boosted lopinavir combined with raltegravir.44 Furthermore, the finding of substantial residual activity of the NRTI backbone, despite the presence of extensive NRTI cross-resistance, was remarkable.44 In addition, this study found that patients with multiple resistance mutations at first-line failure were more likely to become suppressed when exposed to second-line regimens, as compared with those without resistance mutations. This study supported prior data from the SECOND-LINE study45 and observational cohort data,46 which showed that ritonavir-boosted lopinavir-based second-line regimen success is independent of NRTI resistance. This observation can likely be explained by the fact that those patients who do develop resistance have better adherence compared with those without mutations. Therefore, most patients who are empirically switched to a PI-based second-line regimen are likely to suppress, provided good treatment adherence is observed.47

The use of second-generation NNRTIs in third and salvage regimens should be addressed with caution as only one-third of the patients failing NNRTI-based regimens remained fully susceptible to rilpivirine and/or etravirine. Since routine HIV-1 drug resistance testing is not recommended after first-line failure, these NNRTI mutations and resistance profiles might be underestimated when resistance testing is performed at second-line failure, when some of these mutations might have become minority populations, and therefore undetectable by standard genotyping assays, due to the lack of exposure to NNRTIs during second-line treatment.

Despite the careful design of this national HIV drug resistance survey, a cross-sectional survey is known to have some limitations. Only healthcare facilities with at least 500 patients on ART were considered as a collection site; however, these facilities are responsible for 85% of the South African population on ART. Data on the total number of patients failing ART at the time of study in each of the participating facilities were unavailable; therefore, it is unknown which proportion of patients failing ART were surveyed, and baseline characteristics between patients included in the study could not be compared with the overall population failing ART. The study was dependent on the willingness of healthcare workers to obtain written consent from patients, complete study questionnaires and collect specimens. These activities place a burden on clinics that are often already understaffed. This resulted in a relatively prolonged period of sample collection to reach the desired sample size, with some areas, being under-represented by healthcare facilities unable to participate (Table S1). Additionally, the data quality was dependent on the completeness and accuracy of the questionnaires. Treatment history may have been incomplete due to patients transferring from one healthcare facility to another and no direct questions on prevention of mother-to-child transmission exposure were posed. Some bias could have been introduced if healthcare workers were more or less likely to include patients that were suspected of being non-adherent. This potential selection bias could not be verified due to the large number of collection sites. Patients were not systematically sampled at the clinic level; however, we do not believe that it is possible for clinicians to select failing patients based upon demographic or clinical characteristics that would bias the proportion of HIV drug resistance mutations.

Finally, resistance mutations were detected by population-based Sanger sequencing, the sensitivity of this method may have underestimated the prevalence of minority drug resistance mutations, particularly in patients who were non-compliant, or have been exposed to multiple antiretroviral drugs over time. It is known that drug resistance mutations can often not be picked up in the bulk viral population using population-based sequencing, without the presence of antiretroviral drug pressure. It is therefore possible that archived mutations were missed in this analysis.

Overall, given the total sample size of the survey we are confident that these results represent a close reflection of the situation of antiretroviral drug resistance in HIV-1-infected patients in care in the South African ART programme. Based on the findings of this national survey there is currently no evidence to change the recommended first- or second-line regimens.

Funding

This work was funded by The Global Fund to Fight AIDS, Tuberculosis and Malaria (SAF-H-RTC, www.theglobalfund.org). The funder had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

Transparency declarations

None to declare.

Disclaimer

The views described herein are the views of the researchers involved, and do not represent the views or opinions of The Global Fund to Fight AIDS, Tuberculosis and Malaria; nor is there any approval or authorization of this material, express or implied, by The Global Fund to Fight AIDS, Tuberculosis and Malaria.

Supplementary data

Table S1 and Figures S1 and S2 are available as Supplementary data at JAC Online (http://jac.oxfordjournals.org/).

Acknowledgements

Preliminary data from this manuscript were presented as a poster at the International AIDS Society Conference, Vancouver, Canada, 2015 (TUPEB288).

We would like to thank all study participants, healthcare workers and regional ART programme coordinators who made this study possible. In addition, we would like to show our gratitude to E. Letsoalo, A. Lukhwareni, C. Mutuku, L. Skhosana, H. Zwane, A. Bester and M. Claassen, the staff at the HIV drug resistance testing laboratories, who ensured that all samples were processed in a timely manner. In addition, we thank L. Noble for assisting with data capture.

References

2

UNAIDS
.
UNAIDS announces that the goal of 15 million people on life-saving HIV treatment by 2015 has been met nine months ahead of schedule
. .

3

South African Antiretroviral Treatment Guidelines
. .

5

Guidelines for the Use of Antiretroviral Agents in HIV-1-Infected Adults and Adolescents
. .

6

National Consolidated Guidelines for the Prevention of Mother-to-Child Transmission of HIV (PMTCT) and the Management of HIV in Children, Adolescents and Adults
. .

7

Bennett
DE
,
Myatt
M
,
Bertagnolio
S
et al. .
Recommendations for surveillance of transmitted HIV drug resistance in countries scaling up antiretroviral treatment
.
Antivir Ther
2008
;
13
Suppl 2
:
25
–
36
.

8

Barth
RE
,
Wensing
AM
,
Tempelman
HA
et al. .
Rapid accumulation of nonnucleoside reverse transcriptase inhibitor-associated resistance: evidence of transmitted resistance in rural South Africa
.
AIDS
2008
;
22
:
2210
–
2
.

9

El-Khatib
Z
,
DeLong
AK
,
Katzenstein
D
et al. .
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
:
pii=1000117
.

10

Hoffmann
CJ
,
Charalambous
S
,
Sim
J
et al. .
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
:
1928
–
35
.

11

Orrell
C
,
Walensky
RP
,
Losina
E
et al. .
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
:
523
–
31
.

12

Wallis
CL
,
Mellors
JW
,
Venter
WD
et al. .
Varied patterns of HIV-1 drug resistance on failing first-line antiretroviral therapy in South Africa
.
J Acquir Immune Defic Syndr
2010
;
53
:
480
–
4
.

13

Wallis
CL
,
Papathanasopolous
MA
,
Fox
M
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
:
313
–
20
.

14

Dlamini
JN
,
Hu
Z
,
Ledwaba
J
et al. .
Genotypic resistance at viral rebound among patients who received lopinavir/ritonavir-based or efavirenz-based first antiretroviral therapy in South Africa
.
J Acquir Immune Defic Syndr
2011
;
58
:
304
–
8
.

15

Marconi
VC
,
Sunpath
H
,
Lu
Z
et al. .
Prevalence of HIV-1 drug resistance after failure of a first highly active antiretroviral therapy regimen in KwaZulu Natal, South Africa
.
Clin Infect Dis
2008
;
46
:
1589
–
97
.

16

Sigaloff
KC
,
Ramatsebe
T
,
Viana
R
et al. .
Accumulation of HIV drug resistance mutations in patients failing first-line antiretroviral treatment in South Africa
.
AIDS Res Hum Retroviruses
2012
;
28
:
171
–
5
.

17

Singh
A
,
Sunpath
H
,
Green
TN
et al. .
Drug resistance and viral tropism in HIV-1 subtype C-infected patients in KwaZulu-Natal, South Africa: implications for future treatment options
.
J Acquir Immune Defic Syndr
2011
;
58
:
233
–
40
.

18

van Zyl
GU
,
van der Merwe
L
,
Claassen
M
et al. .
Antiretroviral resistance patterns and factors associated with resistance in adult patients failing NNRTI-based regimens in the Western Cape, South Africa
.
J Med Virol
2011
;
83
:
1764
–
9
.

19

Manasa
J
,
Lessells
RJ
,
Skingsley
A
et al. .
High-levels of acquired drug resistance in adult patients failing first-line antiretroviral therapy in a rural HIV treatment programme in KwaZulu-Natal, South Africa
.
PLoS One
2013
;
8
:
e72152
.

20

Hoffmann
CJ
,
Ledwaba
J
,
Li
JF
et al. .
Resistance to tenofovir-based regimens during treatment failure of subtype C HIV-1 in South Africa
.
Antivir Ther
2013
;
18
:
915
–
20
.

21

Skhosana
L
,
Steegen
K
,
Bronze
M
et al. .
High prevalence of the K65R mutation in HIV-1 subtype C infected patients failing tenofovir-based first-line regimens in South Africa
.
PLoS One
2015
;
10
:
e0118145
.

22

Sunpath
H
,
Wu
B
,
Gordon
M
et al. .
High rate of K65R for antiretroviral therapy-naive patients with subtype C HIV infection failing a tenofovir-containing first-line regimen
.
AIDS
2012
;
26
:
1679
–
84
.

23

Van Zyl
GU
,
Liu
TF
,
Claassen
M
et al. .
Trends in genotypic HIV-1 antiretroviral resistance between 2006 and 2012 in South African patients receiving first- and second-line antiretroviral treatment regimens
.
PLoS One
2013
;
8
:
e67188
.

24

Pillay
V
,
Ledwaba
J
,
Hunt
G
et al. .
Antiretroviral drug resistance surveillance among drug-naive HIV-1-infected individuals in Gauteng Province, South Africa in 2002 and 2004
.
Antivir Ther
2008
;
13
Suppl 2
:
101
–
7
.

25

El-Khatib
Z
,
Ekstrom
AM
,
Ledwaba
J
et al. .
Viremia and drug resistance among HIV-1 patients on antiretrovial treatment: a cross-sectional study in Soweto, South Africa
.
AIDS
2010
;
24
:
1679
–
87
.

26

Bronze
M
,
Moodley
K
,
Steegen
K
et al. .
Development and validation of a one-step genotyping protocol encompassing HIV-1 protease and reverse transcriptase for patient antiretroviral drug resistance monitoring
.
First International Conference of the African Society for Laboratory Medicine, Cape Town, South Africa
.

27

Claassen
M
,
Van Zyl
GU
.
The improvement, validation and quality assurance of in-house HIV-1 antiretroviral resistance tests for practical use in resistance monitoring of patients receiving antiretroviral combination therapy at Tygerberg Hospital Family Clinic
. .

28

Lukhwareni
A
,
Bronze
M
,
Ledwaba
J
et al. .
Validation of a robust in-house genotyping assay for the detection of HIV-1 antiretroviral drug resistance
.
Pathpoint Conference, Cape Town, South Africa
.

29

Edgar
RC
.
MUSCLE: multiple sequence alignment with high accuracy and high throughput
.
Nucleic Acids Res
2004
;
32
:
1792
–
7
.

30

Guindon
S
,
Gascuel
O
.
A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood
.
Syst Biol
2003
;
52
:
696
–
704
.

31

Price
MN
,
Dehal
PS
,
Arkin
AP
.
FastTree: computing large minimum evolution trees with profiles instead of a distance matrix
.
Mol Biol Evol
2009
;
26
:
1641
–
50
.

32

Price
MN
,
Dehal
PS
,
Arkin
AP
.
FastTree 2—approximately maximum-likelihood trees for large alignments
.
PLoS One
2010
;
5
:
e9490
.

33

Perno
CF
,
Ceccherini-Silberstein
F
,
De Luca
A
et al. .
Virologic correlates of adherence to antiretroviral medications and therapeutic failure
.
J Acquir Immune Defic Syndr
2002
;
31
Suppl 3
:
S118
–
22
.

34

Turner
D
,
Brenner
BG
,
Routy
JP
et al. .
Rationale for maintenance of the M184V resistance mutation in human immunodeficiency virus type 1 reverse transcriptase in treatment experienced patients
.
New Microbiol
2004
;
27
:
31
–
9
.

35

Wainberg
MA
.
The impact of the M184V substitution on drug resistance and viral fitness
.
Expert Rev Anti Infect Ther
2004
;
2
:
147
–
51
.

36

White
KL
,
Margot
NA
,
Wrin
T
et al. .
Molecular mechanisms of resistance to human immunodeficiency virus type 1 with reverse transcriptase mutations K65R and K65R+M184V and their effects on enzyme function and viral replication capacity
.
Antimicrob Agents Chemother
2002
;
46
:
3437
–
46
.

37

Brenner
BG
,
Oliveira
M
,
Doualla-Bell
F
et al. .
HIV-1 subtype C viruses rapidly develop K65R resistance to tenofovir in cell culture
.
AIDS
2006
;
20
:
F9
–
13
.

38

Coutsinos
D
,
Invernizzi
CF
,
Moisi
D
et al. .
A template-dependent dislocation mechanism potentiates K65R reverse transcriptase mutation development in subtype C variants of HIV-1
.
PLoS One
2011
;
6
:
e20208
.

39

Tang
MW
,
Rhee
SY
,
Bertagnolio
S
et al. .
Nucleoside reverse transcriptase inhibitor resistance mutations associated with first-line stavudine-containing antiretroviral therapy: programmatic implications for countries phasing out stavudine
.
J Infect Dis
2013
;
207
Suppl 2
:
S70
–
7
.

40

Theys
K
,
Vercauteren
J
,
Snoeck
J
et al. .
HIV-1 subtype is an independent predictor of reverse transcriptase mutation K65R in HIV-1 patients treated with combination antiretroviral therapy including tenofovir
.
Antimicrob Agents Chemother
2013
;
57
:
1053
–
6
.

41

Weber
J
,
Chakraborty
B
,
Weberova
J
et al. .
Diminished replicative fitness of primary human immunodeficiency virus type 1 isolates harboring the K65R mutation
.
J Clin Microbiol
2005
;
43
:
1395
–
400
.

42

Miller
V
,
Ait-Khaled
M
,
Stone
C
et al. .
HIV-1 reverse transcriptase (RT) genotype and susceptibility to RT inhibitors during abacavir monotherapy and combination therapy
.
AIDS
2000
;
14
:
163
–
71
.

43

Stone
C
,
Ait-Khaled
M
,
Craig
C
et al. .
Human immunodeficiency virus type 1 reverse transcriptase mutation selection during in vitro exposure to tenofovir alone or combined with abacavir or lamivudine
.
Antimicrob Agents Chemother
2004
;
48
:
1413
–
5
.

44

Paton
NI
,
Kityo
C
,
Hoppe
A
et al. .
Assessment of second-line antiretroviral regimens for HIV therapy in Africa
.
N Engl J Med
2014
;
371
:
234
–
47
.

45

Boyd
MA
,
Kumarasamy
N
,
Moore
CL
et al. .
Ritonavir-boosted lopinavir plus nucleoside or nucleotide reverse transcriptase inhibitors versus ritonavir-boosted lopinavir plus raltegravir for treatment of HIV-1 infection in adults with virological failure of a standard first-line ART regimen (SECOND-LINE): a randomised, open-label, non-inferiority study
.
Lancet
2013
;
381
:
2091
–
9
.

46

Sigaloff
KC
,
Hamers
RL
,
Wallis
CL
et al. .
Second-line antiretroviral treatment successfully resuppresses drug-resistant HIV-1 after first-line failure: prospective cohort in Sub-Saharan Africa
.
J Infect Dis
2012
;
205
:
1739
–
44
.

47

Paton
NI
,
Kityo
C
,
Thompson
J
et al. .
Impact of NRTI cross-resistance on second-line PI+NRTI therapy outcomes in Africa
. In:
Abstracts of the Twenty-second Conference on Retroviruses and Opportunistic Infections, Seattle, WA, USA, 2015
.
Abstract 119. Foundation for Retrovirology and Human Health, Alexandria, VA, USA
.

Supplementary data