Drug Resistance Mutations Against Protease, Reverse Transcriptase and Integrase Inhibitors in People Living With HIV-1 Receiving Boosted Protease Inhibitors in South Africa

The South African national combination antiretroviral therapy (cART) roll-out program started in 2006, with over 4.4 million people accessing treatment since it was first introduced. HIV-1 drug resistance can hamper the success of cART. This study determined the patterns of HIV-1 drug-resistance associated mutations (RAMs) in People Living with HIV-1 (PLHIV-1). Receiving first (for children below 3 years of age) and second-line (for adults) cART regimens in South Africa. During 2017 and 2018, 110 patients plasma samples were selected, 96 samples including those of 17 children and infants were successfully analyzed. All patients were receiving a boosted protease inhibitor (bPI) as part of their cART regimen. The viral sequences were analyzed for RAMs through genotypic resistance testing. We performed genotypic resistance testing (GRT) for Protease inhibitors (PIs), Reverse transcriptase inhibitors (RTIs) and Integrase strand transfer inhibitors (InSTIs). Viral sequences were subtyped using REGAv3 and COMET. Based on the PR/RT sequences, HIV-1 subtypes were classified as 95 (99%) HIV-1 subtype C (HIV-1C) while one sample as 02_AG. Integrase sequencing was successful for 89 sequences, and all the sequences were classified as HIV-1C (99%, 88/89) except one sequence classified CRF02_AG, as observed in PR/RT. Of the 96 PR/RT sequences analyzed, M184V/I (52/96; 54%) had the most frequent RAM nucleoside reverse transcriptase inhibitor (NRTI). The most frequent non-nucleoside reverse transcriptase inhibitor (NNRTI) RAM was K103N/S (40/96, 42%). Protease inhibitor (PI) RAMs M46I and V82A were present in 12 (13%) of the sequences analyzed. Among the InSTI major RAM two (2.2%) sequences have Y143R and T97A mutations while one sample had T66I. The accessory RAM E157Q was identified in two (2.2%). The data indicates that the majority of the patients failed on bPIs didn’t have any mutation; therefore adherence could be major issue in these groups of individuals. We propose continued viral load monitoring for better management of infected PLHIV.

Introduction

Exceptional improvements in combination antiretroviral therapy (cART) regimens have changed HIV/AIDS from a deadly pandemic to a chronic and manageable disease (Trickey et al., 2017). cART has made significant contributions to reducing the rates of morbidity and mortality in people living with HIV (PLHIV) and has led to better management of infection at an individual level, not only in high-income countries but also in low- and middle-income countries (Hightower et al., 2011; UNAIDS, 2017). South Africa’s national HIV cART program was introduced in 2006, with a public health approach (UNAIDS, 2018). Besides problems related to adherence, the development and spread of drug resistance have constantly challenged the long-term management of PLHIV in public health settings, where patients are often monitored using clinical or immunological parameters (Rousseau et al., 2006).

In accordance with the World Health Organization’s (WHO) guidelines, the recommended first-line cART in South Africa consists of a non-nucleoside reverse transcriptase inhibitor (NNRTI) backboned regimen of efavirenz (EFV), combined with two nucleoside reverse transcriptase inhibitors (NRTIs), namely lamivudine (3TC) and either tenofovir disoproxil fumarate (TDF), for adults, or abacavir (ABC), for children. The recommended second-line cART consists of the NRTIs zidovudine (AZT) and 3TC and a ritonavir-boosted (/r) protease inhibitor (PI), usually lopinavir (LPV/r) which was revised to atazanavir (ATV/r) in 2017 (Meintjes et al., 2017). The WHO guidelines also recommend the PI lopinavir co-formulated with ritonavir (lopinavir/ritonavir [LPV/r]) in a four-to-one ratio in first-line cART for children younger than 3 years, based on its superiority when compared with nevirapine (NVP), regardless of previous NVP exposure to prevent mother-to-child HIV transmission (Meintjes et al., 2017).

In vitro studies on PI-naïve PLHIV-1 infected with HIV-1 subtype C (HIV-1C) viruses, have indicated wide variations in their respective susceptibility to the PIs LPV/r and ATV/r (Sutherland et al., 2016). Observational studies from sub-Saharan Africa have shown a 14–32% prevalence of virological failure to second-line boosted PI- (bPI) based cART (Ajose et al., 2012; Sigaloff et al., 2012). In South Africa, reports of drug resistance patterns in patients receiving bPIs are scarce (Collier et al., 2017). With this study, we aimed to identify the pattern of acquired drug resistance mutations (DRMs) among PLHIV in South Africa receiving bPI second-line cART. Furthermore, we characterized the presence of primary integrase strand-transfer inhibitor (InSTI) DRMs in this specific population.

Materials and Methods

Ethics Statement

The study was approved by the Health Research Ethics Committee of Stellenbosch University, South Africa (N15/08/071). The study was conducted according to the ethical guidelines and principles of the Declaration of Helsinki 2013, the South African Guidelines for Good Clinical Practice and the Medical Research Council Ethical Guidelines for Research. A waiver of written informed consent was awarded to conduct sequence analyses on these samples by the Health Research Ethics Committee of Stellenbosch University, South Africa.

Viral Load

HIV-1 Viral load was performed using the Abbott m2000sp and the Abbott m2000rt analyzers (Abbott laboratories, Abbott Park, IL, United States). RNA was isolated from patient samples according to the manufacturer’s instructions using the Abbott RealTime HIV-1 amplification reagent Kit.

Study Design

HIV-1-positive patient samples were obtained randomly, without any knowledge of drug-resistance patterns, from the diagnostic section at the Division of Medical Virology, Stellenbosch University, and the South African National Health Laboratory Services (NHLS). Samples were collected between March 2017 and February 2018. We excluded patient samples with no previous cART regimen history and patients receiving first-line cART treatment regimen. Demographic and clinical information such as age, cART regimen, and viral load measurement (Table 1). Patients had their samples submitted for HIV-1 genotypic resistance testing to the NHLS. The NHLS provides routine genotypic antiretroviral drug resistance testing for clinics from the Western Cape, Gauteng and Eastern Cape provinces.

TABLE 1

Table 1. Characteristics and patterns of mutations in 96 patients at the time of treatment failure.

We included samples from children (aged below 16 years) suspected of failing on bPI – which is used as first-line therapy in children – and adults suspected of experiencing virological failure on a bPI second-line cART regimen, for which treatment information, as provided by the physicians, was available. The treatment history was collected retrospectively. The selection consisted of plasma samples (n = 96) obtained from patients receiving bPIs cART, according to the South African national cART guidelines (Meintjes et al., 2017). These patients are eligible for InSTI treatment consideration when PI mutations are present.

Genotypic Resistance Testing

We performed genotypic resistance testing using viral RNA extracted from plasma. The HIV-1 protease and reverse transcriptase gene fragments were PCR-amplified using a slightly modified protocol as previously described by us (Jacobs et al., 2008). Briefly, HIV-1 protease and reverse transcriptase first-round cDNA synthesis through reverse transcription was done using amplification primers HIV-PR outer 50prot1 (5′-TAA TTT TTT AGG GAA GAT CTG GCC TTC C-3′) and HIVRT outer Mj4 (5′-CTG TTA GTG CTT TGG TTC CTC T–3′), position 2085-2109 and 3399-3420 of the HXB2 reference numbering, with an expected fragment size of approximately 1314 base pairs (bps) (Plantier et al., 2006). For second-round PCR amplification, primers 50prot2 (5′-TCA GAG CAG ACC AGA GCC AAC AGC CCC A–3′) and NE13 (5′-CCT ACT AAC TTC TGT ATG TCA TTG ACA GTC CAG CT–3′), position 2136–2163 and 3334–3300 of the HXB2 reference numbering, with an expected fragment size of approximately 1300 bps, were used (Plantier et al., 2006). The integrase gene fragment amplification steps were performed as previously described by us (Brado et al., 2018). Sequencing reactions were performed as previously described by us (Brado et al., 2018). As part of quality control, each of the viral sequences was inferred on a phylogenetic tree in order to eliminate possible contamination. DRMs were interpreted using the Stanford University HIV Drug Resistance Database version 8.7¹. Subtyping was carried out using REGAv3 and COMET (Pineda-Peña et al., 2013). Phylogenetic analysis was carried out using MEGAv6.

Results

We included patients receiving bPI as part of their cART regimens. We confirmed the successful amplification and Sanger sequencing of the protease, reverse transcriptase and integrase gene fragments of the HIV-1 polymerase gene for 96 samples. In the Integrase region seven sequences did not pass the quality control. Those sequences were excluded from the final analyses. Among the patients, 4% (4/96) of the samples did not indicate their patient file as being either male or female. Hence, they were classified as unknown. The ages ranged from 2 to 66 years.

Seventeen (n = 17; 18%) of the patients were 16 years or younger. Of these patients, three (n = 3; 3%) were female, 12 (13%) were male, while for two (2%) the gender had not been disclosed by the physician. The NRTI cART regimen combinations administered were ABC plus 3TC (n = 6; 6% of patients), AZT plus 3TC (n = 6; 6%), stavudine [d4T] plus 3TC (n = 1; 1%), TDF plus emtricitabine [FTC] (n = 2; 2%), and TDF plus AZT (n = 1; 1%). Fourteen (15%) had received LPV/r, two (2%) had received ATV/r and one (1%) had received darunavir (DRV/r) as their bPIs.

We had a total of 76 adults; 55 (57%) were women, and 22 (23%) were men, while with two (2%) the gender was not disclosed by the physician. The most commonly used NRTI combination was AZT plus 3TC (n = 56; 58%), compared with those receiving ABC plus 3TC (n = 11; 11%) and TDF plus 3TC (n = 8; 8%). Other regimens given were d4T plus 3TC (n = 2; 2%), TDF plus FTC (n = 2; 2%), TDF plus etravirine [ETR] (n = 1; 1%) and AZT plus TDF (n = 1; 1%). Fifty-eight (60%) had received LPV/r and six (6%) had received ATV/r as their bPIs and three (3%) are currently receiving DRV/r.

HIV-1 Subtyping

Subtyping was carried out using REGAv3 and COMET. While REGAv3 provide subtyping for all the RT/PR sequences as well as IN sequences, COMET failed to subtype the IN sequences as the majority of the IN sequences were typed as CPZ. Therefore we used HIV-1 BLAST to identify the nearest subtype for the IN sequences. Based on the PR/RT sequences (n = 96), 99% were identified as HIV-1C while one as 02_AG (PT405ZA). Based on the IN sequences (n = 89), patient samples PT405ZA identified as 02_AG while 88 (99%) samples as HIV-1C. The subtyping data is presented in Supplementary Table S1. We also performed phylogenetic analysis to identify any specific clusters. The neighbor-joining phylogenetic tree did not identify any specific cluster of a transmission network (Figure 1).

FIGURE 1

Figure 1. Neighbor-joining phylogenetic tree. The phylogenetic analysis is carried out using the PR/RT sequences. Bootstrap support > 70% was shown.

NRTI Resistance-Associated Mutations

Table 1 shows the number of resistance-associated mutations (RAMs) observed among the 96 sequences analyzed. We observed M184V/I as the most prevalent NRTI mutation. M184V/I was detected in 52 (54%) patients suspected of failing cART. Of these, 36 (69%) patients were receiving a combination of AZT plus 3TC, compared with 13 (14%) patients receiving an ABC plus 3TC combination. M184V/I was also found in two (4%) patients receiving TDF plus 3TC, compared with those receiving FTC plus TDF (n = 1; 2%). L74V was detected in five (5%) patients – three (3%) patients receiving ABC plus 3TC, and two (2%) patients receiving AZT plus 3TC. The K65R/N mutation occurred in five (5%) patients; K65R occurred in two (2%) patients receiving AZT plus 3TC, and in one (1%) patient receiving TDF plus FTC. The Y115F mutation occurred in five (5%) patients. However, Y115F occurred more often in patients receiving AZT plus 3TC (n = 4; 4%); it occured in one (1%) patient receiving ABC plus 3TC. The A62V mutation occurred in one (1%) patient receiving AZT plus 3TC. V75I occurred in one (1%) patient receiving AZT plus 3TC.

Thymidine analog mutations (TAMs) were grouped into TAMs1 and TAMs2. The most frequent TAMs observed were K70R/E in 10 (10%) patients receiving AZT plus 3TC, in three (3%) patients receiving ABC plus 3TC, and in one (1%) patient receiving FTC plus TDF. D67N occurred in five (5%) patients receiving AZT plus 3TC and in three (3%) patients receiving ABC plus 3TC. D67N occurred in one (1%) patient receiving both TDF plus 3TC and an FTC plus TDF-based regimen. The K219E/Q mutation occurred in six (6%) patients receiving AZT plus 3TC, and in two (2%) patients receiving FTC plus TDF. The M41L mutation occurred in three (3%) patients receiving AZT plus 3TC, and in one (1%) patient receiving 3TC plus ATV/r. T215Y occurred in two (2%) patients receiving AZT plus 3TC.

NNRTI Resistance-Associated Mutations

Table 1 shows the number of NNRTI RAMs observed. We observed that the K103N/S mutation occurred in 41 (43%) of those patients failing cART. Of the patients with this mutation, 24 (59%) patients were receiving AZT plus 3TC, and five (5%) patients were receiving ABC plus 3TC. P225H occurred in 10 (10%) patients receiving AZT plus 3TC, in three (3%) patients receiving ABC plus 3TC, and in two (2%) patients receiving TDF plus 3TC. V106M occurred in 10 (10%) patients receiving AZT plus 3TC, and in two (2%) patients receiving ABC plus 3TC. The patients might have received the EFV or NVP based regiments as their first-line treatment though the past treatment history was not available.

Y188L occurred in five (5%) patients receiving AZT plus 3TC, and in two (2%) patients receiving TDF plus ATV/r, and in one (1%) patient receiving ABC plus 3TC. G190G/A occurred in six (6%) patients receiving AZT plus 3TC, and in three (3%) patients receiving ABC plus 3TC. K101EP occurred in three (3%) patients receiving AZT plus 3TC, and in one (1%) patient receiving ABC plus 3TC.

E138QGA occurred in four (4%) patients receiving AZT plus 3TC, in three (3%) patients receiving ABC plus 3TC, and in one (1%) patient receiving TDF plus 3TC. V108I occurred in one (1%) patient receiving AZT plus 3TC. H221Y and M230L occurred in one (1%) patient each; both these patients were receiving AZT plus 3TC.

PI Resistance-Associated Mutations

Of the 96 patients, 75 (81%) were receiving the LPV/r-containing regimen, followed by the ATV/r- (n = 8; 8%) and DRV/r- (n = 4; 4%) containing regimens. We identified 18 (18%) patients with major PI RAMs. Of those, a substantial majority of 16 (89%) patients were receiving LPV/r as their bPI regimen, while two (11%) patients were receiving ATV/r (Table 1). The most common major PI RAMs observed were M46I and V82A (n = 12; 12%); I54V (n = 10; 10%); I84V and L76V (n = 7; 7%); I47A/V (n = 3; 3%); I50L/V (n = 2; 2%); and V32I (n = 2; 2%) (Table 1).

InSTIs Resistance-Associated Mutations

We successfully sequence 89 INI samples. In our cohort, we identified major InSTI RAMs in patients who had received the DRV/r-containing regimen. Two (2%) patient had a viral sequences with Y143R major InSTI mutation in combination with the accessory T97A mutation, which confers high-level resistance to raltegravir (RAL), intermediate resistance to elvitegravir (EVG), and potential low-level resistance to bictegravir and dolutegravir (DTG). A patient receiving AZT, 3TC, EFV, d4T and LPV/r had a virus sequence with the T66I mutation (Table 1). The mutation identified and classified as an ‘accessory’ integrase, E157Q, occurred in two (2%) patients. One patient was receiving EFV, TDF and 3TC and the other AZT, 3TC and ATV/r. The level of drug resistance against all the cARTs is presented in Figure 2. The observed RAMs in patient receiving bPI based treatment regimen shows high-level resistance was demonstrated in 17 (17%) and 14 (14%) of PLHIV against LPV/r and (ATV/r), respectively, while seven (7%) showed intermediate cross-resistance to DRV/r. Despite off-treatment with NNRTIs, more than half of the patients were shown to have high-level resistance to NVP (57%, n = 56) and EFV (56%, n = 55). High-level resistance were also observed with patients receiving NRTI-based treatment to FTC (60%, n = 59) and 3TC (60%, n = 59).

FIGURE 2

Figure 2. The observed resistance in patients receiving bPIs as part of their cART. High-level resistance was shown in 17 (17%) and 14 (14%) of PLHIV against LPV/r and (ATV/r), respectively, while seven (7%) showed intermediate cross-resistance to DRV/r. Despite off-treatment with NNRTIs, more than half of the patients were shown to have high-level resistance to NVP (57%, n = 56) and EFV (56%, n = 55).

Discussion

In this study, we analyzed 96 HIV-1 RT/PR and 89 IN sequences from PI-experienced and InSTI-naïve patients for major RAMs. The patients were being treated with a bPI second-line cART regimen and were suspected of virological failure. Second-line cART consisted of two NRTIs, backboned by a PI, if previously treated with an NNRTI-based regimen, and vice versa.

As expected, major DRMs against NRTIs and NNRTIs were present at a rate of 65% (n = 65) and 62% (n = 62), respectively. Despite being on a bPI, only 18% (n = 18) of our study sequences harbored major PI RAMs. This is in line with a previous study conducted in Sweden, where it was predicted that HIV-1 subtype C would be more prone to failure in bPIs (Amanda et al., 2016). We identified 27 (27%) sequences not showing any DRM against the drug classes mentioned above and therefore could indicate a problem of poor adherence, rather than the selection of resistant variants.

M184V, the most common NRTI RAM, occurred more frequently in patients receiving AZT plus 3TC, in comparison with patients receiving the ABC plus 3TC regimen. Our findings correspond with previous studies conducted in South Africa with PLHIV, showing M184V/I as the most prevalent NRTI mutation (Marconi et al., 2008; Van Zyl et al., 2011, 2013; Wallis et al., 2011; Lombaard et al., 2016; Neogi et al., 2016; Steegen et al., 2016a, b; Rossouw et al., 2017; Penrose et al., 2018). The K65R and Y115F RAMs occurred more frequently in patients receiving AZT plus 3TC, rather than in patients receiving ABC plus 3TC. TAMS 1 and 2 pathway mutations occurred more frequently in patients receiving the AZT plus 3TC cART regimen but were low in patients receiving ABC plus 3TC, TDF plus 3TC, and FTC plus TDF. The most frequent TAM was K70R/E, which occurred mostly in patients receiving AZT plus 3TC, as opposed to ABC plus 3TC and FTC plus TDF. In our study, M184V, L74V, K65R, and Y115F were the most common major NRTI RAMs in patients receiving LPV/r as their bPIs.

The K103N mutation occurred at a higher frequency in patients receiving AZT plus 3TC or ABC plus 3TC than in those receiving TDF plus 3TC, and 3TC plus d4T. The high rate of K103N RAM is also well documented and has been observed in several previous studies^–(Bronze et al., 2012; Steegen et al., 2017). The group receiving AZT plus 3TC or ABC plus 3TC showed the highest rates of NNRTI mutations such as P225H, V106M, E138A/G/K/Q, G190A/S, and Y188L occurred most frequently in patients receiving AZT plus 3TC or ABC plus 3TC. The presence of NNRTI RAMs when they were off the NNRTI indicating either the carryover NNRTI RAMs from the failed first-line therapy or could be the transmitted DRMs in children from the mother due to vertical transmission prophylaxis.

The majority of our patients were receiving LPV/r as their bPIs. The M46I and V82A RAMs were the most common mutations observed in patients receiving LPV/r compared with ATV/r. We identified a low frequency of M46I and V82A in patients receiving ATV/r as their bPIs. The group receiving AZT plus 3TC or ABC plus 3TC showed the highest rate of PIs such as I54V, I84V, L76V, I47A/V, I50L/V, and V32I. Our findings are in agreement with Neogi et al.’s (2016) study where major PI RAMs were observed in 5% of patients; among them, V82A 65% (28/43), I54V 63% (27/43), L76V 23% (10/43), and L90M 16% (7/43) were the most frequent (Neogi et al., 2016). Our findings also showed M46I and V82A RAMs as the most prevalent major PI RAMs. Furthermore, these findings are in agreement with Rossouw et al. (2015) that observed both V82A and M46I has the most common mutation in infected children receiving PI-based cART. Chimbetete et al. (2018) observed similar results, with all three drug classes showing their DRMs at similar rates. The most common PI RAM reported by Chimbetete et al. (2018) was M46I 28 (33%), followed by I50V 18 (21%) and V82A 18 (21%). We observed more high-level resistance to patients receiving LPV/r compared with ATV/r and DRV/r. Van Zyl et al. (2013) findings are consistent with ours, as they also identified more high-level resistance in patients receiving LPV/r compared with those receiving ATV/r and DRV/r.

We analyzed the integrase gene for the presence of treatment-compromising polymorphisms and DRMs against InSTIs. We observed the presence of Y143R in combination with T97A in one of our patients receiving ABC, 3TC, LPV/r. N155H, Q148H/R/K, and Y143R/C/H are the three major recognized pathways of genotypic resistance against InSTIs (Doyle et al., 2015). We confirmed Y143R in our study and this mutation in combination with T97A also impaired EVG susceptibility and showed possible low-level resistance. Furthermore, our data suggest that EVG activity is compromised in the presence of any RAL RAM, in this case Y143R. We also identified the presence of E157Q in 2 (2%) patients. The presence of this mutation is concerning, as these mutations are associated with potential low-level resistance to both RAL and EVG. In a previous study conducted by Brado et al. (2018), we also found E157Q on HIV-1-infected South African sequences retrieved from the HIV database. Viruses having E157Q were found to be associated with treatment failure of a DTG-containing regimen (Brado et al., 2018). A study has shown that eight patients who had the E157Q mutation and were initiated with DTG-based therapy did not experience viremia suppression below detection level (Neogi et al., 2018).

Furthermore, we identified the presence of T66I mutation in 1 (1%) patient. T66I confers low-level resistance to RAL and high-level resistance to EVG. The low prevalence of DRMs to InSTIs in our cohort should not be underestimated. RAMs against InSTI raises the question about the positioning of DTG in the treatment guidelines for South Africa. Previous studies have shown that a considerable minority of patients develop cross-resistance to DTG after exposure to RAL and EVG; resistance to DTG has not yet been reported in patients from South Africa (Fourati et al., 2015; Mesplède and Wainberg, 2015). As DTG was proposed as the first-line drug, it is essential to conduct studies in real-life clinical settings to identify the efficacy of DTG as limited viral load monitoring and without drug resistance genotyping may compromise the next-generation InSTIs to be used.

Our study had some limitations that merit comments. First, the sample size was small compared to the total number of patients who are receiving cART in South Africa. However, to the best of our knowledge, there has not been any statistical study that reports on the number of patients receiving second-line cART from South Africa. Second, the majority of our patients were receiving LPV/r as their bPIs compared to other bPI regimens. We cannot tell for certain whether the patients having RAL resistance according to the sequence have had access to an RAL-based treatment regimen. Finally, we did not have any adherence data for these patients and the DRM data were only based on population sequencing, therefore we could not detect minor mutations below 20% of the population.

Conclusion

We identified patterns of RAMs against reverse transcriptase inhibitors and PIs from patients suspected of failing on the South African second-line national cART program. Very low or no primary InSTI RAMs were detected in second-line failure patients. The majority of them had M184V mutations that could have been carried over from the first-line cART. Given the negative effect of M184V mutations on viral fitness, it is more plausible to the recycling of 3TC/FTC in second-line cART maintains the presence of M184V. The non-identification of any RAMs in one-third of the patients and the presence of PI RAM in only one-fifth of the patients indicate that the failure may not be due to RAM, but might be due to adherence. Given the limited cART drug availability and high public health burden, we propose for genotypic resistance testing should be performed before switching to InSTIs-based regimen in our setting. This will not only detect treatment failure earlier but will also identify poor treatment adherence. Data generated from this study can assist in the development of cART guidelines for patients who experience treatment failure in resource-limited settings where genotyping is not available. Studies that address operational issues, such as the optimal use of treatment monitoring tools, should be a research priority.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Ethics Statement

The study was approved by the Health Research Ethics Committee of Stellenbosch University, South Africa (N15/08/071). The study was conducted according to the ethical guidelines and principles of the Declaration of Helsinki 2013, the South African Guidelines for Good Clinical Practice and the Medical Research Council Ethical Guidelines for Research. A waiver of written informed consent was awarded to conduct sequence analyses on these samples by the Health Research Ethics Committee of Stellenbosch University, South Africa.

Author Contributions

GJ and UN conceptualized and designed the study. AO performed the laboratory experiments, detailed the sequence analyses, and wrote the first draft of the manuscript. SM helped with demographic data and sample collection. UN performed the sequencing experiments and the sequence analyses. DB, RC, and KS helped with manuscript proofreading and editing. All authors read and approved the final manuscript.

Funding

This study was funded by the National Research Foundation of South Africa, the Poliomyelitis Research Foundation of South Africa, and the NHLS Research Trust. We thank Stellenbosch University and the Harry Crossley Foundation for additional financial support. UN is supported by a grant from the Swedish Research Council (2017-01330).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We acknowledge Shambhu Prasad Ganeshappa Aralaguppe for performing the sequencing experiments. We also acknowledge Mathilda Claassen for sample collection.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00438/full#supplementary-material

TABLE S1 | The subtyping data is presented in Supplementary Table S1. Subtyping was carried out using REGAv3 and COMET. Based on the PR/RT sequences (n = 96), 99% were identified as HIV-1C while one as 02_AG (PT405ZA). Based on the IN sequences (n = 89), patient samples PT405ZA identified as 02_AG while 88 (99%) samples as HIV-1C.

Footnotes

1. ^ https://hivdb.stanford.edu/

References

Ajose, O., Mookerjee, S., Mills, E. J., Boulle, A., and Ford, N. (2012). Treatment outcomes of patients on second-line antiretroviral therapy in resource-limited settings: a systematic review and meta-analysis. AIDS 26, 929–938. doi: 10.1097/QAD.0b013e328351f5b2

PubMed Abstract | CrossRef Full Text | Google Scholar

Amanda, H., Svedhem, V., Singh, K., Sonnerborg, A., and Neogi, U. (2016). Virological failure in patients with HIV1 subtype C receiving antiretroviral therapy: an analysis of a prospective national cohort in Sweden. Lancet HIV 3, 1–18. doi: 10.1186/s40945-017-0033-39

PubMed Abstract | CrossRef Full Text | Google Scholar

Brado, D., Obasa, A. E., Ikomey, G. M., Cloete, R., Singh, K., Engelbrecht, S., et al. (2018). Analyses of HIV-1 integrase sequences prior to South African national HIV-treatment program and availability of integrase inhibitors in Cape Town, South Africa. Sci. Rep. 8:4709. doi: 10.1038/s41598-018-22914-22915

CrossRef Full Text | Google Scholar

Bronze, M., Steegen, K., Wallis, C. L., De Wolf, H., Papathanasopoulos, M. A., Van Houtte, M., et al. (2012). ART-A Consortium. HIV-1 phenotypic reverse transcriptase inhibitor drug resistance test interpretation is not dependent on the subtype of the virus backbone. PLoS One 7:e34708. doi: 10.1371/journal.pone.0034708

PubMed Abstract | CrossRef Full Text | Google Scholar

Chimbetete, C., Katzenstein, D., Shamu, T., Spoerri, A., Estill, J., Egger, M., et al. (2018). HIV-1 drug resistance and third-line therapy outcomes in patients failing second-line therapy in Zimbabwe. Open Forum Infect. Dis. 5, 1–8. doi: 10.1093/ofid/ofy005

PubMed Abstract | CrossRef Full Text | Google Scholar

Collier, D., Iwuji, C., Derache, A., de Oliveira, T., Okesola, N., Calmy, A., et al. (2017). Virological outcomes of second-line protease inhibitor-based treatment for human immunodeficiency virus type 1 in a high-prevalence Rural South African setting: a competing-risks prospective cohort analysis. Clin. Infect. Dis. 64, 1006–1016. doi: 10.1093/cid/cix015

PubMed Abstract | CrossRef Full Text | Google Scholar

Doyle, T., Dunn, D. T., Ceccherini-Silberstein, F., De Mendoza, C., Garcia, F., Smit, E., et al. (2015). Integrase inhibitor (INI) genotypic resistance in treatment-naive and raltegravir-experienced patients infected with diverse HIV-1 clades. J. Antimicrob. Chemother. 70, 3080–3086. doi: 10.1093/jac/dkv243

PubMed Abstract | CrossRef Full Text | Google Scholar

Fourati, S., Charpentier, C., Amiel, C., Morand-Joubert, L., Reigadas, S., Trabaud, M. A., et al. (2015). Cross-resistance to elvitegravir and dolutegravir in 502 patients failing on raltegravir: a French national study of raltegravir-experienced HIV-1-infected patients. J. Antimicrob. Chemother. 70, 1507–1512. doi: 10.1093/jac/dku535

PubMed Abstract | CrossRef Full Text | Google Scholar

Hightower, K. E., Wang, R., Deanda, F., Johns, B. A., Weaver, K., Shen, Y., et al. (2011). Dolutegravir (S/GSK1349572) exhibits significantly slower dissociation than raltegravir and elvitegravir from wild-type and integrase inhibitor-resistant HIV-1 integrase-DNA complexes. Antimicrob. Agents Chemother. 55, 4552–4559. doi: 10.1128/AAC.00157-111

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobs, G. B., Laten, A., van Rensburg, E. J., Bodem, J., Weissbrich, B., Rethwilm, A., et al. (2008). Phylogenetic diversity and low level antiretroviral resistance mutations in HIV type 1 treatment-naive patients from Cape Town, South Africa. AIDS Res. Hum. Retrovirus. 24, 1009–1012. doi: 10.1089/aid.2008.0028

PubMed Abstract | CrossRef Full Text | Google Scholar

Lombaard, J., Bunupuradah, T., Flynn, P. M., Ramapuram, J., Ssali, F., Crauwels, H., et al. (2016). Rilpivirine as a treatment for HIV-infected antiretroviral-naïve adolescents: week 48 safety, efficacy, virology and pharmacokinetics. Pediatr. Infect. Dis. J. 35, 1215–1221. doi: 10.1097/inf.0000000000001275

PubMed Abstract | CrossRef Full Text | Google Scholar

Marconi, V. C., Sunpath, H., Lu, Z., Gordon, M., Koranteng-Apeagyei, K., Hampton, J., et al. (2008). South africa resistance cohort study team. Prevalence of HIV-1 drug resistance after failure of a first highly active antiretroviral therapy regimen in KwaZulu Natal, South Africa. Clin. Infect. Dis. 46, 1589–1597. doi: 10.1086/587109

PubMed Abstract | CrossRef Full Text | Google Scholar

Meintjes, G., Moorhouse, M. A., Carmona, S., Davies, N., Dlamini, S., van Vuuren, C., et al. (2017). Adult antiretroviral therapy guidelines 2017. Southern Afr. J. HIV Med. 18:24. doi: 10.4102/sajhivmed.v18i1.776

PubMed Abstract | CrossRef Full Text | Google Scholar

Mesplède, T., and Wainberg, M. A. (2015). Resistance against integrase strand transfer inhibitors and relevance to HIV persistence. Viruses 7, 3703–3718. doi: 10.3390/v7072790

PubMed Abstract | CrossRef Full Text | Google Scholar

Neogi, U., Engelbrecht, S., Claassen, M., Jacobs, G. B., van Zyl, G., Preiser, W., et al. (2016). Mutational heterogeneity in p6 gag late assembly (L) domains in HIV-1 subtype C viruses from South Africa. AIDS Res. Hum. Retrovirus. 32, 80–84. doi: 10.1089/aid.2015.0266

PubMed Abstract | CrossRef Full Text | Google Scholar

Neogi, U., Singh, K., Aralaguppe, S. G., Rogers, L. C., Njenda, D. T., Sarafianos, S. G., et al. (2018). Ex-vivo antiretroviral potency of newer integrase strand transfer inhibitors cabotegravir and bictegravir in HIV type 1 non-B subtypes. AIDS 32, 469–476. doi: 10.1097/QAD.0000000000001726

PubMed Abstract | CrossRef Full Text | Google Scholar

Penrose, K. J., Brumme, C. J., Scoulos-Hanson, M., Hamanishi, K., Gordon, K., Viana, R. V., et al. (2018). Frequent cross-resistance to rilpivirine among subtype C HIV-1 from first-line antiretroviral therapy failures in South Africa. Antivir. Chem. Chemother. 26:204020661876298. doi: 10.1177/2040206618762985

PubMed Abstract | CrossRef Full Text | Google Scholar

Pineda-Peña, A., Faria, N. R., Imbrechts, S., Libin, P., Abecasis, A. B., Deforche, K., et al. (2013). V. A. automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes: performance evaluation of the new REGA version 3 and seven other tools. Infect. Genet. Evol. 19, 337–348. doi: 10.1016/j.meegid.2013.04.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Plantier, J. C., Lemée, V., Nabias, R., Sire, J. M., and Simon, F. (2006). Generic screening test for HIV infection. AIDS 20, 1345–1347. doi: 10.1097/01.aids.0000232251.66533.ca

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossouw, T. M., Feucht, U. D., Melikian, G., van Dyk, G., Thomas, W., du Plessis, N. M., et al. (2015). Factors associated with the development of drug resistance mutations in HIV-1 infected children failing protease inhibitor-based antiretroviral therapy in South Africa. PLoS One 10:e0133452. doi: 10.1371/journal.pone.0133452

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossouw, T. M., Nieuwoudt, M., Manasa, J., Malherbe, G., Lessells, R. J., Pillay, S., et al. (2017). HIV drug resistance levels in adults failing first-line antiretroviral therapy in an urban and a rural setting in South Africa. HIV Med. 18, 104–114. doi: 10.1111/hiv.12400

PubMed Abstract | CrossRef Full Text | Google Scholar

Rousseau, C. M., Birditt, B. A., McKay, A. R., Stoddard, J. N., Lee, T. C., McLaughlin, S., et al. (2006). Large-scale amplification, cloning and sequencing of near full-length HIV-1 subtype C genomes. J. Virol. Methods 136, 118–125. doi: 10.1016/j.jviromet.2006.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Sigaloff, K. C., Hamers, R. L., Wallis, C. L., Kityo, C., Siwale, M., Ive, P., et al. (2012). Second-line antiretroviral treatment successfully resuppresses drug-resistant HIV-1 after first-line Failure: prospective cohort in sub-saharan Africa. J. Infect. Dis. 205, 1739–1744. doi: 10.1093/infdis/jis261

PubMed Abstract | CrossRef Full Text | Google Scholar

Steegen, K., Bronze, M., Papathanasopoulos, M. A., van Zyl, G., Goedhals, D., Van Vuuren, C., et al. (2016a). Prevalence of antiretroviral drug resistance in patients who are not responding to protease inhibitor – based treatment: results from the first national survey in South Africa. J. Infect. Dis. 214, 1826–1830. doi: 10.1093/infdis/jiw491

PubMed Abstract | CrossRef Full Text | Google Scholar

Steegen, K., Bronze, M., Papathanasopoulos, M. A., van Zyl, G., Goedhals, D., Variava, E., et al. (2017). HIV-1 antiretroviral drug resistance patterns in patients failing NNRTI-based treatment: results from a national survey in South Africa. J. Antimicrob. Chemother. 72, 210–219. doi: 10.1093/jac/dkw358

PubMed Abstract | CrossRef Full Text | Google Scholar

Steegen, K., Carmona, S., Bronze, M., Papathanasopoulos, M. A., van Zyl, G., Goedhals, D., et al. (2016b). Moderate levels of pre-treatment HIV-1 antiretroviral drug resistance detected in the first south african national survey. PLoS One 11:e0166305. doi: 10.1371/journal.pone.0166305

PubMed Abstract | CrossRef Full Text | Google Scholar

Sutherland, K., Collier, D. A., Claiborne, D. T., Prince, J. L., Deymier, M. J., Goldstein, R. A., et al. (2016). Wide variation in susceptibility of transmitted/founder HIV-1 subtype C Isolates to protease inhibitors and association with in vitro replication efficiency. Sci. Rep. 6:38153. doi: 10.1038/srep38153

PubMed Abstract | CrossRef Full Text | Google Scholar

Trickey, A., May, M. T., Vehreschild, J. J., Obel, N., Gill, M. J., Crane, H. M., et al. (2017). Survival of HIV-positive patients starting antiretroviral therapy between 1996 and 2013: a collaborative analysis of cohort studies. Lancet HIV 4:e00349-6. doi: 10.1016/S2352-3018(17)30066-30068

CrossRef Full Text | Google Scholar

UNAIDS (2017). Data 2017. Programme on HIV/AIDS. Geneva: UNAIDS.

Google Scholar

UNAIDS (2018). Country Report South Africa. Geneva: UNAIDS.

Google Scholar

Van Zyl, G. U., Liu, T. F., Claassen, M., Engelbrecht, S., de Oliveira, T., Preiser, W., et al. (2013). 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 8:e67188. doi: 10.1371/journal.pone.0067188

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Zyl, G. U., van der Merwe, L., Claassen, M., Zeier, M., and Preiser, W. (2011). Antiretroviral resistance patterns and factors associated with resistance in adult patients failing NNRTI-based regimens in the Western Cape, South Africa. J. Med. Virol. 83, 1764–1769. doi: 10.1002/jmv.22189

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallis, C. L., Mellors, J. W., Venter, W. D. F., Sanne, I., and Stevens, W. (2011). Protease Inhibitor resistance is uncommon in HIV-1 subtype C infected patients on failing second-line lopinavir/r-containing antiretroviral therapy in South Africa. AIDS Res. Treat. Hindawi 2011:769627. doi: 10.1155/2011/769627

PubMed Abstract | CrossRef Full Text | Google Scholar