Using machine learning and big data to explore the drug resistance landscape in HIV.

Drug resistance mutations (DRMs) appear in HIV under treatment pressure. DRMs are commonly transmitted to naive patients. The standard approach to reveal new DRMs is to test for significant frequency differences of mutations between treated and naive patients. However, we then consider each mutation individually and cannot hope to study interactions between several mutations. Here, we aim to leverage the ever-growing quantity of high-quality sequence data and machine learning methods to study such interactions (i.e. epistasis), as well as try to find new DRMs. We trained classifiers to discriminate between Reverse Transcriptase Inhibitor (RTI)-experienced and RTI-naive samples on a large HIV-1 reverse transcriptase (RT) sequence dataset from the UK (n ≈ 55, 000), using all observed mutations as binary representation features. To assess the robustness of our findings, our classifiers were evaluated on independent data sets, both from the UK and Africa. Important representation features for each classifier were then extracted as potential DRMs. To find novel DRMs, we repeated this process by removing either features or samples associated to known DRMs. When keeping all known resistance signal, we detected sufficiently prevalent known DRMs, thus validating the approach. When removing features corresponding to known DRMs, our classifiers retained some prediction accuracy, and six new mutations significantly associated with resistance were identified. These six mutations have a low genetic barrier, are correlated to known DRMs, and are spatially close to either the RT active site or the regulatory binding pocket. When removing both known DRM features and sequences containing at least one known DRM, our classifiers lose all prediction accuracy. These results likely indicate that all mutations directly conferring resistance have been found, and that our newly discovered DRMs are accessory or compensatory mutations. Moreover, apart from the accessory nature of the relationships we found, we did not find any significant signal of further, more subtle epistasis combining several mutations which individually do not seem to confer any resistance

HIV-1 drug resistance mutations emerging on darunavir therapy in PI-naive and -experienced patients in the UK

BACKGROUND: Darunavir is considered to have a high genetic barrier to resistance. Most darunavir-associated drug resistance mutations (DRMs) have been identified through correlation of baseline genotype with virological response in clinical trials. However, there is little information on DRMs that are directly selected by darunavir in clinical settings. OBJECTIVES: We examined darunavir DRMs emerging in clinical practice in the UK. PATIENTS AND METHODS: Baseline and post-exposure protease genotypes were compared for individuals in the UK Collaborative HIV Cohort Study who had received darunavir; analyses were stratified for PI history. A selection analysis was used to compare the evolution of subtype B proteases in darunavir recipients and matched PI-naive controls. RESULTS: Of 6918 people who had received darunavir, 386 had resistance tests pre- and post-exposure. Overall, 2.8% (11/386) of these participants developed emergent darunavir DRMs. The prevalence of baseline DRMs was 1.0% (2/198) among PI-naive participants and 13.8% (26/188) among PI-experienced participants. Emergent DRMs developed in 2.0% of the PI-naive group (4 mutations) and 3.7% of the PI-experienced group (12 mutations). Codon 77 was positively selected in the PI-naive darunavir cases, but not in the control group. CONCLUSIONS: Our findings suggest that although emergent darunavir resistance is rare, it may be more common among PI-experienced patients than those who are PI-naive. Further investigation is required to explore whether codon 77 is a novel site involved in darunavir susceptibility

Associations between baseline characteristics, CD4 cell count response and virological failure on first-line efavirenz + tenofovir + emtricitabine for HIV

Objectives: The aim of this study was to investigate associations between baseline characteristics and CD4 cell count response on first-line antiretroviral therapy and risk of virological failure (VF) with or without drug resistance. Methods: We conducted an analysis of UK Collaborative HIV Cohort data linked to the UK HIV Drug Resistance Database. Inclusion criteria were viral sequence showing no resistance prior to initiation of first-line efavirenz + tenofovir disoproxil fumarate + emtricitabine and virological suppression within 6 months. Outcomes of VF (≥200 copies/mL) with or without drug resistance were assessed using a competing risks approach fitted jointly with a model for CD4 cell count recovery. Hazard ratios for each VF outcome were estimated for baseline CD4 cell count and viral load and characteristics of CD4 cell count response using latent variables on a standard normal scale. Results: A total of 3640 people were included with 338 VF events; corresponding viral sequences were available in 134 with ≥1 resistance mutation in 36. VF with resistance was associated with lower baseline CD4 (0.30, 0.09-0.62), lower CD4 recovery (0.04, 0.00-0.17) and higher CD4 variability (4.40, 1.22-12.68). A different pattern of associations was observed for VF without resistance, but the strength of these results was less consistent across sensitivity analyses. Cumulative incidence of VF with resistance was estimated to be <2% at 3 years for baseline CD4 ≥350 cells/μL. Conclusion: Lower baseline CD4 cell count and suboptimal CD4 recovery are associated with VF with drug resistance. People with low CD4 cell count before ART or with suboptimal CD4 recovery on treatment should be a priority for regimens with high genetic barrier to resistance

Risk factors and outcomes for the Q151M and T69 insertion HIV-1 resistance mutations in historic UK data

BACKGROUND: The prevalence of HIV-1 resistance to antiretroviral therapies (ART) has declined in high-income countries over recent years, but drug resistance remains a substantial concern in many low and middle-income countries. The Q151M and T69 insertion (T69i) resistance mutations in the viral reverse transcriptase gene can reduce susceptibility to all nucleoside/tide analogue reverse transcriptase inhibitors, motivating the present study to investigate the risk factors and outcomes associated with these mutations. METHODS: We considered all data in the UK HIV Drug Resistance Database for blood samples obtained in the period 1997-2014. Where available, treatment history and patient outcomes were obtained through linkage to the UK Collaborative HIV Cohort study. A matched case-control approach was used to assess risk factors associated with the appearance of each of the mutations in ART-experienced patients, and survival analysis was used to investigate factors associated with viral suppression. A further analysis using matched controls was performed to investigate the impact of each mutation on survival. RESULTS: A total of 180 patients with Q151M mutation and 85 with T69i mutation were identified, almost entirely from before 2006. Occurrence of both the Q151M and T69i mutations was strongly associated with cumulative period of virological failure while on ART, and for Q151M there was a particular positive association with use of stavudine and negative association with use of boosted-protease inhibitors. Subsequent viral suppression was negatively associated with viral load at sequencing for both mutations, and for Q151M we found a negative association with didanosine use but a positive association with boosted-protease inhibitor use. The results obtained in these analyses were also consistent with potentially large associations with other drugs. Analyses were inconclusive regarding associations between the mutations and mortality, but mortality was high for patients with low CD4 at detection. CONCLUSIONS: The Q151M and T69i resistance mutations are now very rare in the UK. Our results suggest that good outcomes are possible for people with these mutations. However, in this historic sample, viral load and CD4 at detection were important factors in determining prognosis

An association between K65R and HIV-1 subtype C viruses in patients treated with multiple NRTIs

Objectives: HIV-1 subtype C might have a greater propensity to develop K65R mutations in patients with virological failure compared with other subtypes. However, the strong association between viral subtype and confounding factors such as exposure groups and ethnicity affects the calculation of this propensity. We exploited the diversity of viral subtypes within the UK to undertake a direct comparative analysis. Patients and methods: We analysed only sequences with major IAS-defined mutations from patients with virological failure. Prevalence of K65R was related to subtype and exposure to the NRTIs that primarily select for this mutation (tenofovir, abacavir, didanosine and stavudine). A multivariate logistic regression model quantified the effect of subtype on the prevalence of K65R, adjusting for previous and current exposure to all four specified drugs. Results: Subtype B patients ( n  =   3410) were mostly MSM (78%) and those with subtype C ( n  =   810) were mostly heterosexual (82%). K65R was detected in 7.8% of subtype B patients compared with 14.2% of subtype C patients. The subtype difference in K65R prevalence was observed irrespective of NRTI exposure and K65R was frequently selected by abacavir, didanosine and stavudine in patients with no previous exposure to tenofovir. Multivariate logistic regression confirmed that K65R was significantly more common in subtype C viruses (adjusted OR = 2.02, 95% CI = 1.55-2.62, P  <   0.001). Conclusions: Patients with subtype C HIV-1 have approximately double the frequency of K65R in our database compared with other subtypes. The exact clinical implications of this finding need to be further elucidated

Virological outcomes of boosted protease inhibitor-based first-line ART in subjects harbouring thymidine analogue-associated mutations as the sole form of transmitted drug resistance

OBJECTIVES: In subjects with transmitted thymidine analogue mutations (TAMs), boosted PIs (PI/b) are often chosen to overcome possible resistance to the NRTI backbone. However, data to guide treatment selection are limited. Our aim was to obtain firmer guidance for clinical practice using real-world cohort data. METHODS: We analysed 1710 subjects who started a PI/b in combination with tenofovir or abacavir plus emtricitabine or lamivudine, and compared their virological outcomes with those of 4889 patients who started an NNRTI (predominantly efavirenz), according to the presence of ≥1 TAM as the sole form of transmitted drug resistance. RESULTS: Participants with ≥1 TAM comprised predominantly MSM (213 of 269, 79.2%), subjects of white ethnicity (206 of 269, 76.6%) and HIV-1 subtype B infections (234 of 269, 87.0%). Most (203 of 269, 75.5%) had singleton TAMs, commonly a revertant of T215Y or T215F (112 of 269, 41.6%). Over a median of 2.5 years of follow-up, 834 of 6599 (12.6%) subjects experienced viraemia (HIV-1 RNA >50 copies/mL). The adjusted HR for viraemia was 2.17 with PI/b versus NNRTI-based therapy (95% CI 1.88–2.51; P < 0.001). Other independent predictors of viraemia included injecting drug use, black ethnicity, higher viral load and lower CD4 cell count at baseline, and receiving abacavir instead of tenofovir. Resistance showed no overall impact (adjusted HR 0.77 with ≥1 TAM versus no resistance; 95% CI 0.54–1.10; P = 0.15). CONCLUSION: In this cohort, patients harbouring ≥1 TAM as the sole form of transmitted drug resistance gained no apparent virological advantage from starting first-line ART with a PI/b

Molecular Epidemiology of HIV-1 Subtype B Reveals Heterogeneous Transmission Risk: Implications for Intervention and Control

Background: The impact of HIV pre-exposure prophylaxis (PrEP) depends on infections averted by protecting vulnerable individuals as well as infections averted by preventing transmission by those who would have been infected if not receiving PrEP. Analysis of HIV phylogenies reveals risk factors for transmission, which we examine as potential criteria for allocating PrEP. Methods: We analyzed 6912 HIV-1 partial pol sequences from men who have sex with men (MSM) in the United Kingdom combined with global reference sequences and patient-level metadata. Population genetic models were developed that adjust for stage of infection, global migration of HIV lineages, and changing incidence of infection through time. Models were extended to simulate the effects of providing susceptible MSM with PrEP. Results: We found that young age <25 years confers higher risk of HIV transmission (relative risk = 2.52 [95% confidence interval, 2.32-2.73]) and that young MSM are more likely to transmit to one another than expected by chance. Simulated interventions indicate that 4-fold more infections can be averted over 5 years by focusing PrEP on young MSM. Conclusions: Concentrating PrEP doses on young individuals can avert more infections than random allocation

Non-disclosed men who have sex with men in UK HIV transmission networks: phylogenetic analysis of surveillance data.

BACKGROUND: Patients who do not disclose their sexuality, including men who do not disclose same-sex behaviour, are difficult to characterise through traditional epidemiological approaches such as interviews. Using a recently developed method to detect large networks of viral sequences from time-resolved trees, we localised non-disclosed men who have sex with men (MSM) in UK transmission networks, gaining crucial insight into the behaviour of this group. METHODS: For this phylogenetic analysis, we obtained HIV pol sequences from the UK HIV Drug Resistance Database (UKRDB), a central repository for resistance tests done as part of routine clinical care throughout the UK. Sequence data are linked to demographic and clinical data held by the UK Collaborative HIV Cohort study and the national HIV/AIDS reporting system database. Initially, we reconstructed maximum likelihood phylogenies from these sequences, then sequences were selected for time-resolved analysis in BEAST if they were clustered with at least one other sequence at a genetic distance of 4·5% or less with support of at least 90%. We used time-resolved phylogenies to create networks by linking together nodes if sequences shared a common ancestor within the previous 5 years. We identified potential non-disclosed MSM (pnMSM), defined as self-reported heterosexual men who clustered only with men. We measured the network position of pnMSM, including betweenness (a measure of connectedness and importance) and assortativity (the propensity for nodes sharing attributes to link). FINDINGS: 14 405 individuals were in the network, including 8452 MSM, 1743 heterosexual women and 1341 heterosexual men. 249 pnMSM were identified (18·6% of all clustered heterosexual men) in the network. pnMSM were more likely to be black African (p<0·0001), less likely to be infected with subtype B (p=0·006), and were slightly older (p=0·002) than the MSM they clustered with. Mean betweenness centrality was lower for pnMSM than for MSM (1·31, 95% CI 0·48-2·15 in pnMSM vs 2·24, 0·98-3·51 in MSM; p=0·002), indicating that pnMSM were in peripheral positions in MSM clusters. Assortativity by risk group was higher than expected (0·037 vs -0·037, p=0·01) signifying that pnMSM were linked to each other. We found that self-reported heterosexual men were more likely to link MSM and heterosexual women than heterosexual women were to link MSM and heterosexual men (Fisher's exact test p=0·0004; OR 2·24) but the number of such transmission chains was small (only 54 in total vs 32 in women). INTERPRETATION: pnMSM are a subgroup distinct from both MSM and from heterosexual men. They are more likely to choose sexual partners who are also pnMSM and might exhibit lower-risk sexual behaviour than MSM (eg, choosing low-risk partners or consistently using condoms). Heterosexual men are the group most likely to be diagnosed with late-stage disease (ie, low CD4 counts) and non-disclosed MSM might put female partners at higher risk than heterosexual men because non-disclosed MSM have male partners. Hence, pnMSM require specific consideration to ensure they are included in public health interventions. FUNDING: National Institutes of Health

Continuation of emtricitabine/lamivudine within combination antiretroviral therapy following detection of the M184V/I HIV-1 resistance mutation

OBJECTIVES: The aim of the study was to investigate whether lamivudine (3TC) or emtricitabine (FTC) use following detection of M184V/I is associated with better virological outcomes. METHODS: We identified people with viruses harbouring the M184V/I mutation in UK multicentre data sets who had treatment change/initiation within 1 year. We analysed outcomes of viral suppression (< 200 HIV-1 RNA copies/mL) and appearance of new major drug resistance mutations (DRMs) using Cox and Poisson models, with stratification by new drug regimen (excluding 3TC/FTC) and Bayesian implementation, and estimated the effect of 3TC/FTC adjusted for individual and viral characteristics. RESULTS: We included 2597 people with the M184V/I resistance mutation, of whom 665 (25.6%) were on 3TC and 458 (17.6%) on FTC. We found a negative adjusted association between 3TC/FTC use and viral suppression [hazard ratio (HR) 0.84; 95% credibility interval (CrI) 0.71-0.98]. On subgroup analysis of individual drugs, there was no evidence of an association with viral suppression for 3TC (n = 184; HR 0.94; 95% CrI 0.73-1.15) or FTC (n = 454; HR 0.99; 95% CrI 0.80-1.19) amongst those on tenofovir-containing regimens, but we estimated a reduced rate of viral suppression for people on 3TC amongst those without tenofovir use (n = 481; HR 0.71; 95% CrI 0.54-0.90). We found no association between 3TC/FTC and detection of any new DRM (overall HR 0.92; 95% CrI 0.64-1.18), but found inconclusive evidence of a lower incidence rate of new DRMs (overall incidence rate ratio 0.69; 95% CrI 0.34-1.11). CONCLUSIONS: We did not find evidence that 3TC or FTC use is associated with an increase in viral suppression, but it may reduce the appearance of additional DRMs in people with M184V/I. 3TC was associated with reduced viral suppression amongst people on regimens without tenofovir

A phylotype-based analysis highlights the role of drug-naive HIV-positive individuals in the transmission of antiretroviral resistance in the UK

International audienceOBJECTIVE:Antiretroviral-naive HIV-positive individuals contribute to the transmission of drug-resistant viruses, compromising first-line therapy. Using phylogenetic inference, we quantified the proportion of transmitted drug-resistance originating from a treatment-naive source.METHODS:Using a novel phylotype-based approach, 24 550 HIV-1 subtype B partial pol gene sequences from the UK HIV Drug Resistance database were analysed. Ongoing transmission of drug resistance amongst HIV-positive individuals was identified as phylotypes of at least three sequences with at least one shared drug resistance mutation, a maximum intra-clade genetic distance of 4.0% and a basal branch support at least 90%. The time of persistence of the transmission chains was estimated using a fast least-squares molecular clock inference approach.RESULTS:Around 70% of transmitted drug-resistance had a treatment-naive source. The most commonly transmitted mutations were L90M in the protease gene and K103N, T215D and T215S in reverse transcriptase. Reversion to wild type occurred at a low frequency and drug-independent reservoirs of resistance have persisted for up to 13 years.CONCLUSION:These results illustrate the impact of viral fitness on the establishment of resistance reservoirs and support the notion that earlier diagnoses and treatment of HIV infections are warranted for counteracting the spread of antiretroviral resistance. Phylotype-based phylogenetic inference is an attractive approach for the routine surveillance of transmitted drug resistance in HIV as well as in other pathogens for which genotypic resistance data are available