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

Long-Read Sequencing with Hierarchical Clustering for Antiretroviral Resistance Profiling of Mixed Human Immunodeficiency Virus Quasispecies

BACKGROUND

HIV infections often develop drug resistance mutations (DRMs), which can increase the risk of virological failure. However, it has been difficult to determine if minor mutations occur in the same genome or in different virions using Sanger sequencing and short-read sequencing methods. Oxford Nanopore Technologies (ONT) sequencing may improve antiretroviral resistance profiling by allowing for long-read clustering.

METHODS

A new ONT sequencing-based method for profiling DRMs in HIV quasispecies was developed and validated. The method used hierarchical clustering of long amplicons that cover regions associated with different types of antiretroviral drugs. A gradient series of an HIV plasmid and 2 plasma samples was prepared to validate the clustering performance. The ONT results were compared to those obtained with Sanger sequencing and Illumina sequencing in 77 HIV-positive plasma samples to evaluate the diagnostic performance.

RESULTS

In the validation study, the abundance of detected quasispecies was concordant with the predicted result with the R2 of > 0.99. During the diagnostic evaluation, 59/77 samples were successfully sequenced for DRMs. Among 18 failed samples, 17 were below the limit of detection of 303.9 copies/μL. Based on the receiver operating characteristic analysis, the ONT workflow achieved an F1 score of 0.96 with a cutoff of 0.4 variant allele frequency. Four cases were found to have quasispecies with DRMs, in which 2 harbored quasispecies with more than one class of DRMs. Treatment modifications were recommended for these cases.

CONCLUSIONS

Long-read sequencing coupled with hierarchical clustering could differentiate the quasispecies resistance profiles in HIV-infected samples, providing a clearer picture for medical care.

Added Value of Next Generation over Sanger Sequencing in Kenyan Youth with Extensive HIV-1 Drug Resistance

HIV-1 drug resistance testing in children and adolescents in low-resource settings is both important and challenging. New (more sensitive) drug resistance testing technologies may improve clinical care, but evaluation of their added value is limited. We assessed the potential added value of using next-generation sequencing (NGS) over Sanger sequencing for detecting nucleoside reverse transcriptase inhibitor (NRTI) and nonnucleoside reverse transcriptase inhibitor (NNRTI) drug resistance mutations (DRMs). Participants included 132 treatment-experienced Kenyan children and adolescents with diverse HIV-1 subtypes and with already high levels of drug resistance detected by Sanger sequencing. We examined overall and DRM-specific resistance and its predicted impact on antiretroviral therapy and evaluated the discrepancy between Sanger sequencing and six NGS thresholds (1%, 2%, 5%, 10%, 15%, and 20%). Depending on the NGS threshold, agreement between the two technologies was 62% to 88% for any DRM, 83% to 92% for NRTI DRMs, and 73% to 94% for NNRTI DRMs, with more DRMs detected at low NGS thresholds. NGS identified 96% to 100% of DRMs detected by Sanger sequencing, while Sanger identified 83% to 99% of DRMs detected by NGS. Higher discrepancy between technologies was associated with higher DRM prevalence. Even in this resistance-saturated cohort, 12% of participants had higher, potentially clinically relevant predicted resistance detected only by NGS. These findings, in a young, vulnerable Kenyan population with diverse HIV-1 subtypes and already high resistance levels, suggest potential benefits of more sensitive NGS over existing technology. Good agreement between technologies at high NGS thresholds supports their interchangeable use; however, the significance of DRMs identified at lower thresholds to patient care should be explored further. IMPORTANCE HIV-1 drug resistance in children and adolescents remains a significant problem in countries facing the highest burden of the HIV epidemic. Surveillance of HIV-1 drug resistance in children and adolescents is an important public health strategy, particularly in resource-limited settings, and yet, it is limited due mostly to cost and infrastructure constraints. Whether newer and more sensitive next-generation sequencing (NGS) adds substantial value beyond traditional Sanger sequencing in detecting HIV-1 drug resistance in real life settings remains an open and debatable question. In this paper, we attempt to address this issue by performing a comprehensive comparison of drug resistance identified by Sanger sequencing and six NGS thresholds. We conducted this study in a well-characterized, vulnerable cohort of children and adolescents living with diverse HIV-1 subtypes in Kenya and, importantly, failing antiretroviral therapy (ART) with already extensive drug resistance. Our findings suggest a potential added value of NGS over Sanger even in this unique cohort.

Drug Resistance Mutations in a Population Before Antiretroviral Therapy Initiation in Northern South Africa

South Africa introduced the "diagnose and treat" universal HIV treatment program in September 2016. This program enables all identified HIV-positive patients to immediately start first-line antiretroviral therapy (ART). However, the presence of drug-resistant (DR) viruses in the drug-naive population complicates the choice of ART. We used next-generation sequencing (NGS) to determine the prevalence and diversity of HIV DR mutations in patients entering HIV treatment programs in northern South Africa. RNA was isolated from plasma of drug-naive HIV-1-infected patients. Using reverse transcriptase polymerase chain reaction, the HIV-1-pol gene comprising the complete protease (PR) and the first 900 bp of reverse transcriptase (RT) was amplified and sequenced on an Illumina MiniSeq platform. Consensus sequences were derived at >20% threshold and at >5% threshold using Geneious PRIME^® software version 2020.1.2. HIV-1 surveillance drug resistance mutations (SDRM) were inferred using Calibrated Population Resistance tool in HIV Drug Resistance Database. Viral subtypes were determined using REGA and RIP genotyping tools. The HIV PR/RT region was successfully sequenced from 241 patients. From these, 23 (9.5%) had at least one SDRM detected at >20% threshold, with a prevalence of 9.5% (n = 18), 3% (n = 7), and 0.4% (n = 1) for non-nucleoside reverse transcriptase inhibitors (NNRTI), nucleoside reverse transcriptase inhibitors (NRTI), and protease inhibitors (PI), respectively. The number of patients with SDRM increased to 31 (12.9%) when minority variants were accounted for at >5% threshold. The most frequent SDRMs based on drug class were; K103N (7.9%-NNRTI), K65R (2.5%-NRTI), and D30N (0.8%-PI). Four cases of dual NRTI/NNRTI mutations were identified. All consensus sequences were subtype C, except three, which were C/A1, C/F1, and C/G recombinants. NGS analysis confirms that individuals entering HIV treatment programs in northern South Africa, habor moderate levels of SDRM, including cases of dual-class drug resistance. Further SDRM studies may be required to better understand resistance in the drug-naive population in the era of "diagnose and treat" in Limpopo Province, South Africa.

Investigation of drug resistance against protease, reverse transcriptase, and integrase inhibitors by next-generation sequencing in HIV-positive patients

Our aim was to investigate the mutations in protease (PR), reverse transcriptase (RT), and integrase (IN) gene regions in human immunodeficiency virus (HIV) using a single amplicon via next-generation sequencing (NGS). The study included plasma samples from 49 HIV-1-positive patients, which were referred for HIV-1 drug resistance testing during 2017. A nested polymerase chain reaction (PCR) was performed after the RNA extraction and one-step reverse transcription stages. The sequencing of the HIV genome in the PR, RT, and IN gene regions was carried out using MiSeq NGS technology. Sanger sequencing (SS) was used to analyze resistance mutations in the PR and RT gene regions using a ViroSeq HIV-1 Genotyping System. PCR products were analyzed with an ABI3500 GeneticAnalyzer (Applied Biosystems). Resistance mutations detected with NGS at frequencies above 20% were identical to the SS results. Resistance to at least one antiretroviral (ARV) drug was 22.4% (11 of 49) with NGS and 10.2% (5 of 49) with SS. At least one low-frequency resistance mutation was detected in 18.3% (9 of 49) of the samples. Low-frequency resistance mutations resulted in virological failure in only one patient. The cost of the analyses was reduced by sample pooling and multiplex analysis using the MiSeq system. This is the first study in Turkey to use NGS technologies for the detection of resistance mutations in all three gene (PR, RT, IN) regions using a single amplicon. Our findings suggest that NGS is more sensitive and cost-effective than the SS method.

Low-Abundance Drug-Resistant HIV-1 Variants in Antiretroviral Drug-Naive Individuals: A Systematic Review of Detection Methods, Prevalence, and Clinical Impact

BACKGROUND

The presence of high-abundance drug-resistant HIV-1 jeopardizes success of antiretroviral therapy (ART). Despite numerous investigations, the clinical impact of low-abundance drug-resistant HIV-1 variants (LA-DRVs) at levels <15%-25% of the virus population in antiretroviral (ARV) drug-naive individuals remains controversial.

METHODS

We systematically reviewed 103 studies assessing prevalence, detection methods, technical and clinical detection cutoffs, and clinical significance of LA-DRVs in antiretroviral drug-naive adults.

RESULTS

In total, 14 919 ARV drug-naive individuals were included. Prevalence of LA-DRVs (ie, proportion of individuals harboring LA-DRVs) was 0%-100%. Technical detection cutoffs showed a 4 log range (0.001%-10%); 42/103 (40.8%) studies investigating the impact of LA-DRVs on ART; 25 studies included only individuals on first-line nonnucleoside reverse transcriptase inhibitor-based ART regimens. Eleven of those 25 studies (44.0%) reported a significantly association between preexisting LA-DRVs and risk of virological failure whereas 14/25 (56.0%) did not.

CONCLUSIONS

Comparability of the 103 studies is hampered by high heterogeneity of the studies' designs and use of different methods to detect LA-DRVs. Thus, evaluating clinical impact of LA-DRVs on first-line ART remains challenging. We, the WHO HIVResNet working group, defined central areas of future investigations to guide further efforts to implement ultrasensitive resistance testing in routine settings.

Prevalence of HIV-1 Integrase Strand Transfer Inhibitor Resistance in Treatment-Naïve Voluntary Counselling and Testing Clients by Population Sequencing and Illumina Next-Generation Sequencing in Taiwan

Purpose

Integrase strand transfer inhibitors (INSTIs) are used as first-line therapy for HIV-1-infected patients. Next-generation sequencing (NGS) can detect low-frequency mutants; however, the clinical value of NGS to detect resistance variants is unknown. This study aimed to evaluate the prevalence of INSTI resistance in southern Taiwan and determine the clinical implications of using NGS to detect integrase region low-level resistant variants.

Patients and Methods

This retrospective cohort study included antiretroviral therapy-naïve HIV-1-infected individuals at Kaohsiung Veterans General Hospital, Taiwan, from 2013 to 2017. Drug-resistance mutations were determined, and an in-house polymerase chain reaction was used for genotyping INSTI resistance. NGS was used to assess INSTI resistance (≧1%), and the results were compared with those from population sequencing. Drug resistance-associated mutations were defined according to the 2019 IAS-USA HIV drug resistance-associated mutations list, and accessory mutations by a Stanford HIVdb score ≥10 to at least one INSTI.

Results

A total of 224 patients were included. Subtype B HIV-1 strains were found in 96% of the individuals and subtype CRF01_AE in 4%. The prevalence rates for nucleoside reverse transcriptase inhibitors, non-nucleoside reverse transcriptase inhibitors, protease inhibitors and INSTI resistance were 4%, 5.8%, 0.4% and 0.9%, respectively. The most common INSTI resistance-associated mutations were G163K (0.4%) and E138A (0.4%). Of the 38 patients diagnosed in 2017 who had both NGS and population sequencing data, none had INSTI resistance-associated mutations by population sequencing; however, NGS detected four more INSTI resistance-associated mutations with low frequencies (G163R 3.25%, S153F 3.21%, S153Y 1.36% and Y143H 2.06%). Two patients with S153F and S153Y low frequencies mutations started INSTI-based highly active antiretroviral therapy, and none had virological failure by week 48.

Conclusion

Our findings showed a low rate of HIV drug resistance to INSTIs (0.9%) in treatment-naïve patients. NGS detected more INSTI resistance-associated mutations at a low frequency. Low-level drug resistance-associated mutations to INSTIs identified by NGS did not have an impact on the treatment response to INSTI-based first-line therapy.

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.

A single variant sequencing method for sensitive and quantitative detection of HIV-1 minority variants

HIV drug resistance is a major threat to achieving long-term viral suppression in HIV-positive individuals. Drug resistant HIV variants, including minority variants, can compromise response to antiretroviral therapy. Many studies have investigated the clinical relevance of drug resistant minority variants, but the level at which minority variants become clinically relevant remains unclear. A combination of Primer-ID and deep sequencing is a promising approach that may quantify minority variants more accurately compared to standard deep sequencing. However, most studies that used the Primer-ID method have analyzed clinical samples directly. Thus, its sensitivity and quantitative accuracy have not been adequately validated using known controls. Here, we constructed defined proportions of artificial RNA and virus quasispecies and measured their relative proportions using the Primer-ID based, quantitative single-variant sequencing (qSVS) assay. Our results showed that minority variants present at 1% of quasispecies were detected reproducibly with minimal variations between technical replicates. In addition, the measured frequencies were comparable to the expected frequencies. These data validate the accuracy and reproducibility of the qSVS assay in quantifying authentic HIV minority variants, and support the use of this approach to examine the impacts of minority HIV variants on virologic response and clinical outcome.

Analysis of unusual and signature APOBEC-mutations in HIV-1 pol next-generation sequences

Introduction

At low mutation-detection thresholds, next generation sequencing (NGS) for HIV-1 genotypic resistance testing is susceptible to artifactual detection of mutations arising from PCR error and APOBEC-mediated G-to-A hypermutation.

Methods

We analyzed published HIV-1 pol Illumina NGS data to characterize the distribution of mutations at eight NGS mutation detection thresholds: 20%, 10%, 5%, 2%, 1%, 0.5%, 0.2%, and 0.1%. At each threshold, we determined proportions of amino acid mutations that were unusual (defined as having a prevalence <0.01% in HIV-1 group M sequences) or signature APOBEC mutations.

Results

Eight studies, containing 855 samples, in the NCBI Sequence Read Archive were analyzed. As detection thresholds were lowered, there was a progressive increase in the proportion of positions with usual and unusual mutations and in the proportion of all mutations that were unusual. The median proportion of positions with an unusual mutation increased gradually from 0% at the 20% threshold to 0.3% at the 1% threshold and then exponentially to 1.3% (0.5% threshold), 6.9% (0.2% threshold), and 23.2% (0.1% threshold). In two of three studies with available plasma HIV-1 RNA levels, the proportion of positions with unusual mutations was negatively associated with virus levels. Although the complete set of signature APOBEC mutations was much smaller than that of unusual mutations, the former outnumbered the latter in one-sixth of samples at the 0.5%, 1%, and 2% thresholds.

Conclusions

The marked increase in the proportion of positions with unusual mutations at thresholds below 1% and in samples with lower virus loads suggests that, at low thresholds, many unusual mutations are artifactual, reflecting PCR error or G-to-A hypermutation. Profiling the numbers of unusual and signature APOBEC pol mutations at different NGS mutation detection thresholds may be useful to avoid selecting a threshold that is too low and poses an unacceptable risk of identifying artifactual mutations.

Development of a new comprehensive HIV-1 genotypic drug resistance assay for all commercially available reverse transcriptase, protease and integrase inhibitors in patients infected with group M HIV-1 strains

Comprehensive PCR assays for the genotypic drug resistance analysis of all HIV-1 antiretroviral agents (reverse transcriptase, protease and integrase inhibitors) are increasingly in demand due to introduction of integrase inhibitors in the first line regimens and the increasing presence of non-B HIV-1 clades around the world. This study focused on the development and evaluation of a new PCR-based assay for the amplification and sequencing of the entire HIV-1 pol region of major circulating group M HIV-1 strains in Europe for genotypic drug resistance analysis. The comprehensive touchdown PCR assay developed in this study utilized HIV-1 RNA extracted from the plasma of blood samples of consenting HIV-1 infected patients in Cyprus, collected from 2017 to 2019. The HIV-1 pol region was amplified by touchdown PCR for both the primary RT-PCR and the secondary PCR steps. Successful PCR amplicons were determined by population DNA sequencing, using the Sanger method and the genotypic drug resistance analysis was performed with the Stanford University HIV Drug Resistance Database Program. The newly developed assay successfully amplified the entire HIV-1 pol region (2844 nucleotides long) of 141 out of 144 samples of group M HIV-1 subtypes and recombinant strains of the Cyprus HIV-1 Transmission Cohort Study (CHICS) isolated from 2017 to 2019 and genotypic analyses were conducted for all currently available HIV-1 reverse transcriptase, protease and integrase inhibitors. The drug resistance, epidemiological and demographic data of these study subjects will be expanded upon in the CHICS (L.G. Kostrikis et al., manuscript in preparation for publication). The newly developed HIV-1 genotypic drug resistance assay would benefit clinical settings, and research focusing on the world-wide spread of HIV-1 drug-resistant strains, especially in geographic regions characterized by polyphyletic HIV-1 infections.

Companion and complementary diagnostics for infectious diseases

INTRODUCTION

Companion diagnostics (CDx) are important in oncology therapeutic decision-making, but specific regulatory-approved CDx for infectious disease treatment are officially lacking. While not approved as CDx, several ID diagnostics are used as CDx. The diagnostics community, manufacturers, and regulatory agencies have made major efforts to ensure that diagnostics for new antimicrobials are available at or near release of new agents.

AREAS COVERED

This review highlights the status of Complementary and companion diagnostic (c/CDx) in the infectious disease literature, with a focus on genotypic antimicrobial resistance testing against pathogens as a class of diagnostic tests.

EXPERT OPINION

CRISPR, sepsis markers, and narrow spectrum antimicrobials, in addition to current and emerging technologies, present opportunities for infectious disease c/CDx. Challenges include slow guideline revision, high costs for regulatory approval, lengthy buy in by agencies, discordant pharmaceutical/diagnostic partnerships, and higher treatment costs. The number of patients and available medications used to treat different infectious diseases is well suited to support competing diagnostic tests. However, newer approaches to treatment (for example, narrow spectrum antibiotics), may be well suited for a small number of patients, i.e. a niche market in support of a CDx. The current emphasis is rapid and point-of-care (POC) diagnostic platforms as well as changes in treatment.

MiDRMpol: A High-Throughput Multiplexed Amplicon Sequencing Workflow to Quantify HIV-1 Drug Resistance Mutations against Protease, Reverse Transcriptase, and Integrase Inhibitors

The detection of drug resistance mutations (DRMs) in minor viral populations is of potential clinical importance. However, sophisticated computational infrastructure and competence for analysis of high-throughput sequencing (HTS) data lack at most diagnostic laboratories. Thus, we have proposed a new pipeline, MiDRMpol, to quantify DRM from the HIV-1 pol region. The gag-vpu region of 87 plasma samples from HIV-infected individuals from three cohorts was amplified and sequenced by Illumina HiSeq2500. The sequence reads were adapter-trimmed, followed by analysis using in-house scripts. Samples from Swedish and Ethiopian cohorts were also sequenced by Sanger sequencing. The pipeline was validated against the online tool PASeq (Polymorphism Analysis by Sequencing). Based on an error rate of <1%, a value of >1% was set as reliable to consider a minor variant. Both pipelines detected the mutations in the dominant viral populations, while discrepancies were observed in minor viral populations. In five HIV-1 subtype C samples, minor mutations were detected at the <5% level by MiDRMpol but not by PASeq. MiDRMpol is a computationally as well as labor efficient bioinformatics pipeline for the detection of DRM from HTS data. It identifies minor viral populations (<20%) of DRMs. Our method can be incorporated into large-scale surveillance of HIV-1 DRM.

Validation of publicly-available software used in analyzing NGS data for HIV-1 drug resistance mutations and transmission networks in a Washington, DC, Cohort

The DC Cohort is an ongoing longitudinal observational study of persons living with HIV. To better understand HIV-1 drug resistance and potential transmission clusters among these participants, we performed targeted, paired-end next-generation sequencing (NGS) of protease, reverse transcriptase and integrase amplicons. We elected to use free, publicly-available software (HyDRA Web, Stanford HIVdb and HIV-TRACE) for data analyses so that laboratory personnel without extensive bioinformatics expertise could use it; making the approach accessible and affordable for labs worldwide. With more laboratories transitioning away from Sanger-based chemistries to NGS platforms, lower frequency drug resistance mutations (DRMs) can be detected, yet their clinical relevance is uncertain. We looked at the impact choice in cutoff percentage had on number of DRMs detected and found an inverse correlation between the two. Longitudinal studies will be needed to determine whether low frequency DRMs are an early indicator of emerging resistance. We successfully validated this pipeline against a commercial pipeline, and another free, publicly-available pipeline. RT DRM results from HyDRA Web were compared to both SmartGene and PASeq Web; using the Mantel test, R2 values were 0.9332 (p<0.0001) and 0.9097 (p<0.0001), respectively. PR and IN DRM results from HyDRA Web were then compared with PASeq Web only; using the Mantel test, R2 values were 0.9993 (p<0.0001) and 0.9765 (p<0.0001), respectively. Drug resistance was highest for the NRTI drug class and lowest for the PI drug class in this cohort. RT DRM interpretation reports from this pipeline were also highly correlative compared to SmartGene pipeline; using the Spearman's Correlation, rs value was 0.97757 (p<0.0001). HIV-TRACE was used to identify potential transmission clusters to better understand potential linkages among an urban cohort of persons living with HIV; more individuals were male, of black race, with an HIV risk factor of either MSM or High-risk Heterosexual. Common DRMs existed among individuals within a cluster. In summary, we validated a comprehensive, easy-to-use and affordable NGS approach for tracking HIV-1 drug resistance and identifying potential transmission clusters within the community.