Identification of minority resistance mutations in the HIV-1 integrase coding region using next generation sequencing

Rate of response to initial antiretroviral therapy according to level of pre-existing HIV-1 drug resistance detected by next-generation sequencing in the strategic timing of antiretroviral treatment (START) study

OBJECTIVES

The main objective of this analysis was to evaluate the impact of pre-existing drug resistance by next-generation sequencing (NGS) on the risk of treatment failure (TF) of first-line regimens in participants enrolled in the START study.

METHODS

Stored plasma from participants with entry HIV RNA >1000 copies/mL were analysed using NGS (llumina MiSeq). Pre-existing drug resistance was defined using the mutations considered by the Stanford HIV Drug Resistance Database (HIVDB v8.6) to calculate the genotypic susceptibility score (GSS, estimating the number of active drugs) for the first-line regimen at the detection threshold windows of >20%, >5%, and >2% of the viral population. Survival analysis was conducted to evaluate the association between the GSS and risk of TF (viral load >200 copies/mL plus treatment change).

RESULTS

Baseline NGS data were available for 1380 antiretroviral therapy (ART)-naïve participants enrolled over 2009-2013. First-line ART included a non-nucleoside reverse transcriptase inhibitor (NNRTI) in 976 (71%), a boosted protease inhibitor in 297 (22%), or an integrase strand transfer inhibitor in 107 (8%). The proportions of participants with GSS <3 were 7% for >20%, 10% for >5%, and 17% for the >2% thresholds, respectively. The adjusted hazard ratio of TF associated with a GSS of 0-2.75 versus 3 in the subset of participants with mutations detected at the >2% threshold was 1.66 (95% confidence interval 1.01-2.74; p = 0.05) and 2.32 (95% confidence interval 1.32-4.09; p = 0.003) after restricting the analysis to participants who started an NNRTI-based regimen.

CONCLUSIONS

Up to 17% of participants initiated ART with a GSS <3 on the basis of NGS data. Minority variants were predictive of TF, especially for participants starting NNRTI-based regimens.

Emergence of the G118R Pan-Integrase Resistance Mutation as a Result of Low Compliance to a Dolutegravir-Based cART

HIV-1 resistance towards integrase inhibitors is a potential threat of the success of HIV-1 combination treatment. G118R is a rare drug resistance mutation conferring pan-integrase resistance. Here, we describe the occurrence of G118R in a HIV-1 subtype-B-positive individual with major compliance problems, detected while the patient was on dolutegravir-based cART. We speculate the pre-selection of M184I/V aided the occurrence of G118R in this case, and discuss the robustness of dolutegravir-based therapies.

In vivo drug resistance mutation dynamics from the early to chronic stage of infection in antiretroviral-therapy-naïve HIV-infected men who have sex with men

Human immunodeficiency virus type 1 (HIV) primary drug resistance mutations (DRMs) influence the long-term therapeutic effects of antiretroviral treatment (ART). Drug-resistance genotyping based on polymerase gene sequences obtained by next-generation sequencing (NGS) was performed using samples from 10 ART-naïve HIV-infected men who have sex with men (MSM; P1-P10) from the acute/early to chronic stage of infection. Three of the 10 subjects exhibited the presence of major (abundance, ≥ 20%) viral populations carrying DRM at early/acute stage that later, at the chronic stage, dropped drastically (V106M) or remained highly abundant (E138A). Four individuals exhibited additional DRMs (M46I/L; I47A; I54M, L100V) as HIV minority populations (abundance, 2-20%) that emerged during the chronic stage but ephemerally.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

Prevalence and clinical impact of minority resistant variants in patients failing an integrase inhibitor-based regimen by ultra-deep sequencing

Background

Integrase strand transfer inhibitors (INSTIs) are recommended by international guidelines as first-line therapy in antiretroviral-naive and -experienced HIV-1-infected patients.

Objectives

This study aimed at evaluating the prevalence at failure of INSTI-resistant variants and the impact of baseline minority resistant variants (MiRVs) on the virological response to an INSTI-based regimen.

Methods

Samples at failure of 134 patients failing a raltegravir-containing (n = 65), an elvitegravir-containing (n = 20) or a dolutegravir-containing (n = 49) regimen were sequenced by Sanger sequencing and ultra-deep sequencing (UDS). Baseline samples of patients with virological failure (VF) (n = 34) and of those with virological success (VS) (n = 31) under INSTI treatment were sequenced by UDS. Data were analysed using the SmartGene platform, and resistance was interpreted according to the ANRS algorithm version 27.

Results

At failure, the prevalence of at least one INSTI-resistant variant was 39.6% by Sanger sequencing and 57.5% by UDS, changing the interpretation of resistance in 17/134 (13%) patients. Among 53 patients harbouring at least one resistance mutation detected by both techniques, the most dominant INSTI resistance mutations were N155H (45%), Q148H/K/R (23%), T97A (19%) and Y143C (11%). There was no difference in prevalence of baseline MiRVs between patients with VF and those with VS. MiRVs found at baseline in patients with VF were not detected at failure either in majority or minority mutations.

Conclusions

UDS is more sensitive than Sanger sequencing at detecting INSTI MiRVs at treatment failure. The presence of MiRVs at failure could be important to the decision to switch to other INSTIs. However, there was no association between the presence of baseline MiRVs and the response to INSTI-based therapies in our study.

Viral evolution in the cell-associated HIV-1 DNA during early ART can lead to drug resistance and virological failure in children

Using cell-associated DNA and cell-free RNA of human immunodeficiency virus type-1 (HIV-1), we investigated the role of drug-resistant viral variants that emerged during early antiretroviral therapy (ART) in determining virological outcome. This case-control study compared virologic nonresponder children (two viral loads [VLs] ≥ 200 copies/mL within 2 years of ART) and responder children (two VLs < 200 copies/mL after six months of ART) infected with HIV-1 initiated on nonnucleoside reverse-transcriptase inhibitor (NNRTI)-based ART. The partial reverse-transcriptase gene of HIV-1 in cell-associated DNA was genotyped using next-generation sequencing (NGS; Illumina; threshold 0.5%; at baseline and month six of ART) and in cell-free RNA (concurrently and at virological failure; VL > 1000 copies/mL at ≥ 12 months of ART) using the Sanger method. Among 30 nonresponders and 37 responders, baseline differences were insignificant while adherence, VL, and drug resistance mutations (DRMs) observed at month six differed significantly ( P ≥ 0.05). At month six, NGS estimated a higher number of DRMs compared with Sanger (50% vs 33%; P = 0.001). Among the nonresponders carrying a resistant virus (86.6%) at virological failure, 26% harbored clinically relevant low-frequency DRMs in the cell-associated DNA at month six (0.5%-20%; K103N, G190A, Y181C, and M184I). Plasma VL of > 3 log ₁₀ copies/mL (AOR, 30.4; 95% CI, 3.3-281; P = 0.003) and treatment-relevant DRMs detected in the cell-associated DNA at month six (AOR, 24.2; 95% CI, 2.6-221; P = 0.005) were independently associated with increased risk for early virological failure. Our findings suggest that treatment-relevant DRMs acquired in cell-associated DNA during the first six months of ART can predict virological failure in children initiated on NNRTI-based ART.

Limited Marginal Utility of Deep Sequencing for HIV Drug Resistance Testing in the Age of Integrase Inhibitors

HIV drug resistance genotyping is a critical tool in the clinical management of HIV infections. Although resistance genotyping has traditionally been conducted using Sanger sequencing, next-generation sequencing (NGS) is emerging as a powerful tool due to its ability to detect low-frequency alleles.

HIV drug resistance genotyping is a critical tool in the clinical management of HIV infections. Although resistance genotyping has traditionally been conducted using Sanger sequencing, next-generation sequencing (NGS) is emerging as a powerful tool due to its ability to detect low-frequency alleles. However, the clinical value added from NGS approaches to antiviral resistance testing remains to be demonstrated. We compared the variant detection capacity of NGS versus Sanger sequencing methods for resistance genotyping in 144 drug resistance tests (105 protease-reverse transcriptase tests and 39 integrase tests) submitted to our clinical virology laboratory over a four-month period in 2016 for Sanger-based HIV drug resistance testing. NGS detected all true high-frequency drug resistance mutations (>20% frequency) found by Sanger sequencing, with greater accuracy in one instance of a Sanger-detected false positive. Freely available online NGS variant callers HyDRA and PASeq were superior to Sanger methods for interpretations of allele linkage and automated variant calling. NGS additionally detected low-frequency mutations (1 to 20% frequency) associated with higher levels of drug resistance in 30/105 (29%) protease-reverse transcriptase tests and 4/39 (10%) integrase tests. In clinical follow-up of 69 individuals for a median of 674 days, we did not find a difference in rates of virological failure between individuals with and without low-frequency mutations, although rates of virological failure were higher for individuals with drug-relevant low-frequency mutations. However, all 27 individuals who experienced virological failure reported poor adherence to their drug regimen during the preceding follow-up time, and all 19 who subsequently improved their adherence achieved viral suppression at later time points, consistent with a lack of clinical resistance. In conclusion, in a population with low antiviral resistance emergence, NGS methods detected numerous instances of minor alleles that did not result in subsequent bona fide virological failure due to antiviral resistance.

Clinically relevant thresholds for ultrasensitive HIV drug resistance testing: a multi-country nested case-control study

BACKGROUND

Implementation of ultrasensitive HIV drug resistance tests for routine clinical use is hampered by uncertainty about the clinical relevance of drug-resistant minority variants. We assessed different detection thresholds for pretreatment drug resistance to predict an increased risk of virological failure.

METHODS

We did a case-control study nested within a prospective multicountry cohort. Our study included patients from 12 clinical sites in Kenya, Nigeria, South Africa, Uganda, and Zambia. We defined cases as patients with virological failure (ie, those who had either viral load ≥400 copies per mL at 12 months or had switched to second-line antiretroviral therapy [ART] as a result of virological failure before 12 months) and controls as those with viral suppression (viral load <400 copies per mL at 12 months) on first-line non-nucleoside reverse transcriptase inhibitor-based antiretroviral therapy. We assessed pretreatment drug resistance with Illumina MiSeq next-generation sequencing, using the International Antiviral Society (IAS)-USA mutation list or the Stanford HIV Drug Resistance Database (HIVDB) genotypic sensitivity score. We calculated diagnostic accuracy measures and assessed the odds of virological failure using conditional logistic regression for 1%, 5%, and 10% pretreatment drug resistance detection thresholds, compared with the conventional 20% or more used in Sanger-based sequencing.

FINDINGS

Paired viral load results before ART and at month 12 of follow-up were available from 1896 participants. We identified 178 patients with virological failure and selected 338 matched controls. We excluded 117 patients from pretreatment drug resistance analysis; therefore, 152 cases of virological failure and 247 controls were included in the final analysis. With the IAS-USA mutation list, at a detection threshold of 20% or more in patients with pretreatment drug resistance, the adjusted odds ratio (OR) for virological failure was 9·2 (95% CI 4·2-20·1) compared with those without pretreatment drug resistance. Lowering the threshold resulted in adjusted ORs of virological failure of 6·8 (95% CI 3·3-13·9) at the 10% threshold, 7·6 (3·4-17·1) at the 5% threshold, and 4·5 (2·0-10·2) at the 1% threshold. Lowering the detection threshold from 20% improved the sensitivity (ie, ability to identify cases) from 12% (n=18) to 13% (n=19) at detection threshold 10%, to 15% (n=23) at detection threshold 5%, and to 17% (n=26) at detection threshold 1%, but caused a slight reduction in specificity (ie, ability to identify controls) from 98% (n=241) to 96% (n=238) at the 10% threshold, 96% (n=236) at the 5% threshold, and a larger reduction to 92% (n=227) at the 1% threshold. Diagnostic ORs were 5·4 (95% CI 2·1-13·9) at the 20% threshold, 3·8 (1·7-8·6) at the 10% threshold, 3·8 (1·8-8·1) at the 5% threshold, and 2·3 (1·2-4·2) at the 1% threshold. Use of the Stanford HIVDB genotypic sensitivity scores yielded similar ORs for virological failure, sensitivities, specificities, and diagnostic ORs.

INTERPRETATION

Ultrasensitive resistance testing for pretreatment drug resistance improved identification of people at risk of virological failure; however, this came with a reduction in our ability to identify people with viral suppression, especially at very low thresholds. Further modelling is needed to estimate the optimal trade-off for the 5% and 20% thresholds, balancing improved case finding against unnecessary regimen switching.

FUNDING

The Netherlands Ministry of Foreign Affairs, IrsiCaixa, and European Union.

Performance comparison of deep sequencing platforms for detecting HIV-1 variants in the pol gene

The present study compares the performances of an in-house sequencing protocol developed on MiSeq, the Sanger method, and the 454 GS-FLX for detecting and quantifying drug-resistant mutations (DRMs) in the human immunodeficiency virus polymerase gene (reverse transcriptase [RT] and protease [PR]). MiSeq sequencing identified all the resistance mutations detected by bulk sequencing (n = 84). Both the MiSeq and 454 GS-FLX platforms identified 67 DRMs in the RT and PR regions, but a further 25 DRMs were identified by only one or other of them. Pearson's analysis showed good concordance between the percentage of drug-resistant variants determined by MiSeq and 454 GS-FLX sequencing (ρ = .77, P < .0001). The MiSeq platform is as accurate as the 454 GS-FLX Roche system for determining RT and PR DRMs and could be used for monitoring human immunodeficiency virus type 1 drug resistance.

The Role of HIV-1 Drug-Resistant Minority Variants in Treatment Failure

Human immunodeficiency virus type 1 (HIV-1) drug resistance genotyping is recommended to help in the selection of antiretroviral therapy and to prevent virologic failure. There are several ultrasensitive assays able to detect HIV-1 drug-resistance minority variants (DRMVs) not detectable by standard population sequencing-based HIV genotyping assays. Presence of these DRMVs has been shown to be clinically relevant, but its impact does not appear to be uniform across drug classes. In this review, we summarize key evidence for the clinical impact of DRMVs across drug classes for both antiretroviral treatment-naive and antiretroviral treatment-experienced patients, and highlight areas where more supporting evidence is needed.

Development of a Qualitative Quantitative Polymerase Chain Reaction Test to Identify Patients Failing First-Line Therapy to Non-Nucleotide Reverse Transcriptase Inhibitor

Antiretroviral therapy (ART) can be compromised by selection of drug resistance strains, which can be promoted by lack of adherence during therapy and drug tolerance, and some of these drug-resistant strains can persist for years as minority populations. The K103N drug resistance mutation is selected by the use of non-nucleotide reverse transcriptase inhibitors, including nevirapine or efavirenz (EFV), used in low-income countries. Here we describe the use of a less expensive qualitative point mutation polymerase chain reaction (PMqPCR_K103N) targeting K103N mutation. To validate the use of this methodology, we tested previously sequenced samples from patients treated with highly active ART with viral loads above 2,000 copies/ml and compared the results of our assay with Illumina deep sequencing. Due to its low cost and high specificity, this test is particularly suitable for low-income countries to screen for pretreatment resistance in patients either initiating ART or failing first-line regimens containing EFV.

Clinical and biological insights from viral genome sequencing

Whole-genome sequencing (WGS) of pathogens is becoming increasingly important not only for basic research but also for clinical science and practice. In virology, WGS is important for the development of novel treatments and vaccines, and for increasing the power of molecular epidemiology and evolutionary genomics. In this Opinion article, we suggest that WGS of viruses in a clinical setting will become increasingly important for patient care. We give an overview of different WGS methods that are used in virology and summarize their advantages and disadvantages. Although there are only partially addressed technical, financial and ethical issues in regard to the clinical application of viral WGS, this technique provides important insights into virus transmission, evolution and pathogenesis.

Clinical benefit of dolutegravir in HIV-1 management related to the high genetic barrier to drug resistance

This manuscript reviews the reasons why Integrase inhibitors should now routinely constitute a part of first line antiretroviral therapy for the treatment of HIV disease. The use of these drugs that are generally well tolerated has resulted in far less drug resistance than was the case with most other categories of antiviral compounds. In addition, the integrase inhibitor family of drugs has been less prone to the problem of transmitted drug resistance which is due to a wide variety of substitutions in the HIV genome that can be sexually transmitted from one person to another. However, the use of integrase inhibitors in first line therapy may unfortunately not soon happen in developing country settings where non-nucleoside reverse transcriptase inhibitors continue to be a mainstay of initial therapy, primarily for reasons of cost. As long as this situation continues, problems of drug resistance and transmitted drug resistance will be common in such settings. Current evidence also suggests that the use of dolutegravir as a first line integrase inhibitor may limit development of drug resistance to an extent that exceeds the use of other members of this family of drugs. This may be due to particular patterns of resistance involving dolutegravir, whereby the mutations that are associated with resistance against this compound may actually diminish both HIV replication capacity as well as integrase enzymatic activity in a far-reaching and unique manner. This gives potential hope that the use of dolutegravir in first line therapy could actually form part of the long-sought goal of attainment of a functional cure for HIV disease.