|Home | About | Journals | Submit | Contact Us | Français|
The free antiretroviral therapy (ART) program in India has scaled up to register second largest number of people living with HIV/AIDS across the globe. To assess the effectiveness of current first-line regimen we estimated virological suppression on completion of 1 year of ART. The study describes the correlates of virological failure (VF) and multinucleoside reverse transcriptase inhibitor (NRTI) drug resistance mutations (DRMs).
In this cross-sectional study conducted between June and August 2014, consecutive adults from 4 State sponsored ART clinics of western India were recruited for plasma viral load screening at 12±2 months of ART initiation. Individuals with plasma viral load >1000copies/mL were selected for HIV drug resistance (HIVDR) genotyping. Logistic regression analyses were performed to assess factors associated with VF and multi-NRTI resistance mutations. Criteria adopted for multi-NRTI resistance mutation were either presence of K65R or 3 or more thymidine analog mutations (TAMs) or presence of M184V along with 2 TAMs.
Of the 844 study participants, virological suppression at 1 year was achieved in 87.7% of individuals. Factors significantly associated with VF (P<0.005) were 12 months CD4 count of ≤100cells/μL (adjusted OR −7.11), low reported adherence (adjusted OR −4.44), and those living without any partner (adjusted OR −1.98). In patients with VF, the prevalence of non-nucleoside reverse transcriptase inhibitor (NNRTI) DRM (78.75%) were higher as compared to NRTI (58.75%). Multi-NRTI DRMs were present in 32.5% of sequences and were significantly associated with CD4 count of ≤100cells/μL at baseline (adjusted OR −13.00) and TDF-based failing regimen (adjusted OR −20.43). Additionally, low reported adherence was negatively associated with multi-NRTI resistance (adjusted OR −0.11, P=0.015). K65R mutation was significantly associated with tenofovir (TDF)-based failing regimen (P<0.001).
The study supports early linkage of HIV-infected individuals to the program for ART initiation, adherence improvement, and introduction of viral load monitoring. With recent introduction of TDF-based regimen, the emergence of K65R needs to be monitored closely among HIV-1 subtype C-infected Indian population.
India has an estimated population of 2.11 million (1.71–2.64 million) HIV-infected individuals with adult (15–49 years) prevalence of 0.26% in 2015. Initiated on April 1, 2004, the program of provision of free antiretroviral therapy (ART) to the eligible individuals, has scaled up considerably to register second largest number of people living with HIV/AIDS (PLHA) across the globe. Compared to 2000 and 2007, the Indian National AIDS Control Program (NACP) has, respectively, achieved 66% and 32% decline in new cases of HIV infection in 2015. With maturing of the program it is essential to assess the effectiveness of ART by estimating virological suppression among individuals initiated on first-line regimen. Concern over emergence and spread of HIV drug resistance (HIVDR), also necessitates resistance surveillance studies to monitor treatment outcome.
The national program in India implements public health approach for initiation of first-line ART comprising of 2 nucleoside/nucleotide analog reverse transcriptase inhibitors (NRTI) and a non-nucleoside reverse transcriptase inhibitor (NNRTI). The NRTI options available in the program are zidovudine (AZT), abacavir (ABC), lamivudine (3TC), stavudine (d4T), didanosine (DDI), and emtricitabine (FTC), whereas tenofovir (TDF) is the only nucleotide reverse transcriptase inhibitor. The NNRTI options currently available are efavirenz (EFV) and nevirapine (NVP).
In developed nations, viral load (VL) monitoring is a part of standard treatment protocol, whereas our program still relies on immunological and clinical monitoring, though targeted viral load is offered for confirmation of treatment failure. In the absence of periodic plasma viral load estimation, detection of treatment failure is delayed which may result in sequential accumulation of resistance mutations, especially thymidine analog mutations (TAMs). A large study conducted in Sub-Saharan Africa and Southeast Asia indicated virological failure (VF) rate of 11.1% at 12 months of ART. A recent study from southern India among 599 participants recorded cumulative VF incidence of 13.2% in the first year, however drug resistance outcome data were unavailable. A prior study from 2 ART clinics in southern and western India documented viral load suppression (<1000copies/mL) in 75% and 64.6% patients, respectively, at the end of 12 months.
In the absence of routine virological monitoring, HIVDR outcome studies from India at 12 months of ART are limited by sample size, different inclusion criteria and variable duration of ART treatment before genotyping.[9–13] In one of the prior large study involving genotyping of 138 patients with failure from South India, M184V and Y181C emerged as most common NRTI and NNRTI mutations, respectively. Sinha et al also reported overall drug resistance mutation (DRM) prevalence of 93.8% among 128 individuals from North India with failure of first-line ART. There is a need to have systematically collected data as recommended by WHO to assess effectiveness of first-line ART regimen at 12 months of its initiation.
Conventionally, mutations like K69 Insertion, Q151M complex and multiple TAMs reduce susceptibility to all currently available NRTI. Additionally K65R mutation selected by TDF, ABC, d4T, and DDI impart resistance to all NRTI except AZT. In subtype C-infected Indian population, K69 Insertion and Q151M are seldom reported and common mutations responsible for multi-NRTI resistance includes, K65R and multiple TAMs with or without M184V. M184V is selected by 3TC/FTC and reduces susceptibility to these drugs by 100-fold. In combination with TAMs, M184V also reduces susceptibility to ABC and DDI.[18,20] As NRTI forms an important backbone of both first and second-line regimen, there is a need to estimate the proportion of individuals developing multi-NRTI resistance mutations at 12 months of ART.
Therefore the primary objective of this study was to ascertain the rate of population level virological suppression and factors associated with failure at 12 months of ART. The secondary objective was to characterize HIVDR mutation pattern with special attention to multi-NRTI resistance mutations.
This cross-sectional study was carried out between June and August 2014 at 4 State sponsored ART clinics of Pune city in western India. Consecutive ART naive adults (age >18 years) were recruited on completion of 12±2 months of ART for viral load screening. All ART clinics followed the National guidelines wherein, ART was initiated at CD4 count of ≤350cells/μL or else irrespective of CD4 count for individuals in WHO clinical stage III and IV. Individuals with history of exposure to antiretrovirals (ARV) outside program were excluded, except in cases of females, where ARV were offered for prevention of mother to child transmission. The primary study outcome was VF at the end of 1 year of ART initiation and development of resistance mutations was secondary outcome.
The details of participant's age, gender, marital status, ARV treatment history, and immunological profile were recorded from the data maintained at respective clinics. Married and those in “Living-in” relationship were considered to be living with partner. Unmarried, divorced, and separated couples were considered as living without partner. A counselor elicited history of prior ART exposure and adherence to the medication at the time of sample collection. The drug adherence was assessed at the time of monthly drug refill by calculating the proportion of doses missed since last visit. An average of recent 3 months adherence was taken into account for analysis. Ten milliliters of peripheral blood was collected in EDTA collection vials for CD4 cell count, viral load testing, and drug resistance genotyping. Plasma viral load estimation was performed on m2000RT Abbott Real Time HIV-1 assay (Abbott Molecular, Inc., Des Plaines, IL). Individuals with VF, defined as single plasma HIV-1 viral load of more than 1000copies/mL were subjected to HIVDR genotyping.
Genotyping to characterize DRMs was performed by an In-house population-based sequencing protocol, validated previously. Bidirectional sequencing of complete protease (PR, 1–99 amino acids) and partial reverse transcriptase (RT, 1–256 amino acids) were performed on ABI 3730XL DNA Sequencer (Applied Biosystems, Inc., Foster City, CA). The quality of PR/RT sequences were screened using Sequence Quality Analysis Tool (SQUAT) and any pair of sequences with genetic distance of less than 0.99 or more than 15 were excluded. Resistance patterns were determined using Stanford University HIVdb Genotypic Interpretation Algorithm version 7.0 dated February 27, 2014. Predicted susceptibility for each NRTI and NNRTI was calculated by adding up penalty score associated with each DRM in a given sequence. A total score of 14 or less was considered susceptible and a score of 15 or more was considered resistant. Phylogenetic trees were constructed by the maximum likelihood method based on general time reversible model, using MEGA 6.0. Multi-NRTI resistance mutations were defined as either presence of K65R or 3 or more TAMs or presence of M184V along with 2 TAMs.
The demographic, clinical, and biological characteristics of study participants were summarized by medians and interquartile range (IQR) for continuous variables and by proportions for categorical variables. Differences between groups were compared using the Mann–Whitney U test for continuous data and Pearson Chi-square test for categorical data. Univariable logistic regression analysis was performed to assess factors associated with VF and multi-NRTI resistance mutations. Variables with significant association (P<0.05) with outcome in univariable regression analysis were included in multivariable logistic regression model and were adjusted for age and gender. Adjusted odds ratio (aOR) with 95% confidence intervals (95% CIs) were calculated. All significant tests were 2-sided and “P” value of <0.05 was considered statistically significant. Data were entered in a spreadsheet (MS Excel 2010) and analyzed using SPSS version 17.0.
Minimum sample size calculated for our study, assuming proportion of individuals with VF of 15%, precision of 2.5% and alpha error as 5%, was 784. Anticipating nonparticipation by 10%, we targeted recruitment of 863 individuals for our study.
The study was approved by Ethics committee of National AIDS Research Institute (ICMR) and the individuals from each ART clinics were recruited voluntarily after obtaining written informed consent.
A total 844 individuals consented to participate at 12±2 months of ART initiation as shown in study flow chart (Fig. (Fig.1).1). Compared to baseline median CD4 count of 213cells/μL (IQR: 116–309), significant increase in median CD4 count at 12 months of ART (379cells/μL, IQR: 256–526) was observed (Table (Table1).1). At the time of ART initiation 681 (80.7%) and 160 (19.0%) individuals were on AZT- and TDF-based regimen, respectively. Toxicity was the main reason for which 43.7% of individual underwent drug substitution in ART during first year.
Study flow chart showing participation of 844 people living with HIV/AIDS (PLHA) at 12±2 months of ART.
Demographic and clinical characteristics of study participants.
On completion of 12 months of ART, successful suppression of plasma viral load to less than 1000copies/mL was achieved in 87.7% (95% CI: 85.5–89.9%) of individuals. Among 740 virologically suppressed individuals, 58 (7.83%) patients satisfied the criteria for immunological failure. Of these 50 (6.75%) had CD4 cell count of below baseline value suggestive of immuno-virological discordance and 17 (2.28%) had 12 months CD4 count of ≤100cells/μL. Nine patients in common satisfied both the criteria of 12 months CD4 count of <100cells/μL as well as below baseline value. Out of 104 individuals with VF, only 32 satisfied the immunological failure criteria. First 85 plasma samples from VF patients were processed for HIVDR genotyping due to budgetary constraints and 80 were sequenced successfully.
The factors significantly associated (P<0.005) with VF (Table (Table2)2) were 12 months CD4 count of <100cells/μL (aOR −7.11; 95% CI: 3.10–16.31), <95% reported adherence (aOR −4.44; 95% CI: 2.74–7.18) and those living without any partner (OR −1.98; 95% CI: 1.24–3.14).
Logistic regression analyses of factors associated with virological failure in (n=844) individuals at completion of 12±2 months of first-line antiretroviral therapy.
The HIV-1 subtype C was the most predominant subtype, seen in 97.5% of sequences. The phylogenetic tree of 80 partial pol gene sequences along with reference sequences retrieved from Los Alamos database is shown in Supplementary Figure. A single sequence (KR816094) clustered with subtype A1 and another sequence (KR816077) was separated away from Indian subtype C cluster suggestive of recombination. Near full-length genomic sequencing of this isolate (KR816077) confirmed unique recombination pattern of CRF01_AE and subtype C.
Of these 80 successfully genotyped sequences, 17 (21.25%) sequences did not show any DRM. The prevalence of DRMs were 58.75%, 78.75%, and 1.25% for NRTIs, NNRTIs, and protease inhibitors (PIs), respectively (Fig. (Fig.2).2). The most common NRTI and NNRTI mutation were M184V/I (51.25%) and K103N (36.25%), respectively. Only 1 individual had major PI resistance mutation L90M and 5 had minor PI resistance mutations L89M, V77I, L63P, H69K/R/Q, M36I, K20I/M/R/T, G16E, and L10V/I. Out of 80, 17.5% sequences showed at least 1 TAMs and there was no significant difference in the mutations arising from TAM-1 (16.25%) and TAM-2 (21.25%) pathway.
Prevalence of nucleoside reverse transcriptase inhibitors (NRTIs) and nonnucleoside reverse transcriptase inhibitors (NNRTIs) drug resistance mutations among individuals with virological failures at 12 months of ART initiation.
Based on Stanford resistance scoring system (Fig. (Fig.3),3), among NRTI, the predicted susceptibility of AZT was highest (86.25%) followed by TDF (68.75%). Comparison of mutation pattern among individuals failing on TDF- and AZT-based regimen (Supplementary Table) showed statistically significant association of K65R mutation with TDF-based regimen (P<0.001). Among NNRTI, the predicted susceptibility for NVP as well as EFV was only 23.75%. Cross-resistance to second-generation NNRTI, namely etravirine (ETR) and rilpivirine (RPV) was seen in 51.25% and 60% of sequences, respectively.
Predicted efficacy of nucleoside reverse transcriptase inhibitors (NRTI) and nonnucleoside reverse transcriptase inhibitors (NNRTI) following virological failure, based on Stanford genotype susceptibility scoring.
Factors significantly associated with presence of multi-NRTI resistance mutations (Table (Table3)3) were CD4 count of ≤100cells/μL at baseline (aOR −13; 95% CI: 2.26–74.77) and TDF-based failing regimen (aOR −20.43, 95% CI: 4.21–99.14). In addition low reported adherence (<95%) was negatively associated with multi-NRTI resistance (aOR −0.046; 95% CI: 0.02–0.66).
Logistic regression analyses of determinants of multi-NRTI drug resistance mutations (DRM) in individuals with virological failure (n=80) at 12 months of ART initiation.
Since 2013, WHO has recommended plasma viral load testing as preferred method to monitor individuals on ART.[25,27] To our knowledge, this is the largest study from India describing virological response and drug resistance outcome after 12 months of initiation of first-line ART. Among 844 participants, who were part of free ART program, we report successful virological suppression in 87.6%, which is well above WHO recommended target of 70%. Prior studies with limited sample size reported VF rate of 19%, 25%, and 35.4%, respectively, at ART clinics of Pune, Chennai, and Mumbai.[9,13] Similar VF rate of 12% was reported in Kigali, Rwanda. McMahon et al carried out the first systematic review to quantify population-level viral suppression 12 months after ART initiation in low and middle income countries (LMIC). In nine cohorts with viral suppression threshold of <1000copies/mL, the review showed suppression in 83.5% (95% CI: 77.8–88.4; n=3192) of the combined on-treatment populations. Present study reports better virological suppression, that is, in 87.7% of individuals as compared to other LMIC.
In this study, 58 (7.83%) virologically suppressed individuals failed to show adequate immunological recovery at 12 months of ART. Of these, 17 (2.28%) of individuals had 12 months CD4 count of ≤100cells/μL and 50 (6.75%) had fall in CD4 cell count to baseline or below. Zoufaly et al have indicated increased risk of developing non-AIDS morbidity and mortality among those individuals who failed to achieve CD4 count of above baseline value. In absence of virological monitoring, individuals with immunovirological discordance would be missed and misclassified as treatment failure. Also out of 104 individuals with VF, only 32 meet the WHO immunological failure criteria, indicating the need for plasma viral load testing for early identification of treatment failure.
In the current setting, low adherence to ART (<90%), low CD4 count (<100cells/μL) and living without a partner emerged as important predictors of VF among PLHA. Several prior studies have indicated suboptimal adherence as important factor for VF in resource-limited settings.[8,33,34] Our study showed that VF was more likely among PLHA who are living alone rather than those living with their partners. This highlights the need for strengthening social support measures in country where HIV/AIDS is still associated with as stigma.
The overall prevalence of DRM in individuals with VF was 78.25%, with NNRTI DRMs leading the list. A prior multicountry WHO HIV drug resistance survey (n=269) reported overall prevalence of 72.1% in patients failing therapy at 12 months with resistance to NNRTI and NRTI as, 69.5% and 62.5%, respectively. The most common NRTI mutations reported were M184V and K65R whereas K103N/S and V106A/M were common NNRTI resistance mutations. The pattern was similar to as observed by other investigators in resource-limited settings.[9,33,35] Thus, the pattern of DRMs is predictable at 12 months of ART in individuals with access to virological monitoring.
In this study, 21.25% sequences from individuals with VL >1000copies/mL, did not contain any DRM, suggestive of unreported treatment interruption. VF in the absence of DRM has been reported from India and elsewhere.[13,28] Failure without DRM also support the need of virological monitoring complemented with drug resistance genotyping or else 1 in 5 patients with VF may be subjected to costly second-line regimen unnecessarily.
A prior large study involving 6 Sub-Saharan African countries revealed prevalence of resistance to NNRTI and NRTI as 60.6% and 57.8%, respectively. Higher prevalence of NNRTI resistance mutations (79%) in our study can be explained by its low genetic barrier. Earlier, Neogi et al performed temporal analysis of primary drug resistance in India and reported rising trend of NNRTI resistance. Importantly, 51.25% and 60% of sequences respectively also indicated low to high level of resistance to second-generation NNRTI, ETR and rilpivirine (RPV). These NNRTI analogs are not yet introduced in the program and emergence of cross-resistance will reduce their future utility. High cross-resistance to ETR and RPV has been reported previously.[38,39]
TAMs are important in resource-limited settings as they impart resistance to all available NRTIs and the degree of cross-resistance depends on specific mutations and total number of TAMs. In presence of TAMs, M184V also causes additional low-level resistance to DDI and ABC. In this study only 5% of sequences from VF patients had ≥3 TAMs, due to which the susceptibility to AZT was preserved. WHO guidelines also recommends use of AZT in the second line, if TDF is utilized in failing first-line regimen. Accumulation of TAMs in case of immunological monitoring jeopardize NRTI backbone leaving PI as only active ingredient in second line. Introduction of newer class of ARV in second-line regimen can be averted if failure is detected early by virological monitoring. Although recent trials have indicated promising results of simplified PI monotherapy, the success of this strategy ultimately depends on viral load monitoring.
In resource-limited settings NRTI analogs are utilized in first as well as second-line regimen. Due to low prevalence of multi drug mutation such as T69Ins or Q151M, we adopted criteria for multi-NRTI DRMs as either presence of K65R or presence of 2 TAMs along with M184V or presence of 3 or more TAMs. Multivariable analysis indicated that low CD4 count at baseline and TDF-based failing regimen were significantly associated with multi-NRTI DRMs. Recently a multicentre retrospective cohort study has reported strong association of low pre-ART CD4 count (<100cells/μL) with development of TDF resistance. The median CD4 count of study participants at ART initiation was 213cells/μL (IQR: 116–309) as against the program recommendation of 350cells/μL for ART initiation. Though early diagnosis and linkage to ART program is challenging in resource-limited settings, initiating ART at higher CD4 count will likely to improve virological suppression and resistance outcome.[27,44,45] A similar study from Asia found multi-NRTI resistant associated mutations (RAMs) in 37% of the patients, wherein, multi-NRTI RAMs were defined as presence of either Q151M; 69Ins; 2 TAMs; or M184V+1 TAM. In our study, with modified definition, multi-NRTI DRMs were present in 32.5% of sequences.
The emergence of K65R among individuals failing TDF-based regimen is primarily responsible for multi-NRTI resistance. In this study, 52.94% of individuals who were exposed to TDF-developed K65R mutation. Though TDF-based regimen is associated with higher rate of viral suppression as compared to AZT, one must be cautious of the emergence of K65R mutation among individuals failing TDF-based regimen. These results have implications on program as recently we adopted fixed-dose combination of TDF, 3TC, and EFV as primary initiating regimen based on WHO recommendations.[26,27] In a recent study from South Africa, authors have concluded that patients failing on a TDF-containing regimen were almost 5 times more likely to present with a K65R mutation compared to d4T-exposed patients. The selection of K65R is known to be facilitated in subtype C, which is a predominant circulating subtype in India. The lower genetic barrier in subtype C for K65R may be attributed to enzymatic pausing arising at the end of poly-adenine stretches. In addition, recent reports have suggested antagonism between K65R and TAMs, indicating that both pathways are unlikely to occur simultaneously.
Surprisingly among individuals with VF, higher reported adherence (> 95%) was associated with multi-NRTI DRMs. As we captured adherence data of recent 3 months before failure, it is possible that these individuals had periods of suboptimal adherence initially. In a failing individual high level of adherence provides environment for selective drug pressure. Thus in absence of virological monitoring, continuation of failing regimen despite higher adherence may further lead to development of DRMs. Nevertheless intensive efforts should be made to ensure optimal adherence among individuals initiating ART.
Our study does have few limitations. Firstly, pretreatment drug resistance was not assessed due to financial constraints, which may influence the final outcome. However, recent studies from India indicates transmitted drug resistance to be less than 5%.[54–56] In this study, VF was diagnosed by single plasma viral load estimation at 12±2 months. Thought 2 tests are preferable, most studies resort to single tests as cost-saving measure. Use of population-based genotyping method in this study may underestimate the prevalence of resistance. Finally, it should be recorded that the ART clinics involved in this study were located in large city, catering for urban population and therefore, findings may not be nationally representative.
To conclude, implementation of viral load monitoring into the program is a long-standing priority. Though, this study achieved WHO recommended target for viral load suppression, the program should focus on timely ART initiation and optimal adherence. Psychosocial support of partner is essential for improvement in virological outcome among PLHA. The study highlights the importance of AZT as preferred NRTI option in second-line ART due to selection of K65R by TDF. With recent introduction of fixed-dose combination of TDF, 3TC, and EFV in national program, the emergence of K65R need to be monitored closely among HIV-1 subtype C-infected Indian population.
The authors would like to thank our counselors Girish Rahane, Narayan Panchal, and Ayesha Momin, staff of Model Colony NACO ART Centre and all participants for valuable assistance.
Abbreviations: 3TC = lamivudine, ABC = abacavir, ART = antiretroviral therapy, ARV = antiretrovirals, AZT = zidovudine, d4T = stavudine, DDI = didanosine, DRMs = drug resistance mutations, EFV = efavirenz, ETR = etravirine, FTC = emtricitabine, NACP = National AIDS Control Program, India, NNRTI = non-nucleoside reverse transcriptase inhibitor, NRTI = nucleoside reverse transcriptase inhibitor, NVP = nevirapine, PLHA = people living with HIV/AIDS, RPV = rilpivirine, TAMs = thymidine analog mutations, TDF = tenofovir, VF = virological failure.
Funding: This work was supported by World Health Organization (India) and an intramural grant from Indian Council of Medical Research, New Delhi.
Authors’ contribution: RRG, BBR, SKK, MVG, and RSP designed and planned the research. MVG, NG, DBK, SS, and RRG supported the study implementation and data collection. SKK, DNC, SSK, and SNK performed the lab experiments. MVG, SKK, SG, RJ, SSS, and RRG designed the database and performed the analysis. SKK and MVG drafted the first manuscript. All authors commented on design, data interpretation, critically reviewed the article draft, read and approved the final manuscript. SKK and MVG contributed equally to this work.
The GenBank accession numbers obtained for the newly submitted pol gene sequences are KR816018–KR816097.
The authors have no conflicts of interest to disclose.
Supplemental Digital Content is available for this article.