Comparison of HIV-1 drug-resistance genotyping by ultra-deep sequencing and sanger sequencing using clinical samples

Frequency matters: comparison of drug resistance mutation detection by Sanger and next-generation sequencing in HIV-1

BACKGROUND

Next-generation sequencing (NGS) is gradually replacing Sanger sequencing (SS) as the primary method for HIV genotypic resistance testing. However, there are limited systematic data on comparability of these methods in a clinical setting for the presence of low-abundance drug resistance mutations (DRMs) and their dependency on the variant-calling thresholds.

METHODS

To compare the HIV-DRMs detected by SS and NGS, we included participants enrolled in the Swiss HIV Cohort Study (SHCS) with SS and NGS sequences available with sample collection dates ≤7 days apart. We tested for the presence of HIV-DRMs and compared the agreement between SS and NGS at different variant-calling thresholds.

RESULTS

We included 594 pairs of SS and NGS from 527 SHCS participants. Males accounted for 80.5% of the participants, 76.3% were ART naive at sample collection and 78.1% of the sequences were subtype B. Overall, we observed a good agreement (Cohen's kappa >0.80) for HIV-DRMs for variant-calling thresholds ≥5%. We observed an increase in low-abundance HIV-DRMs detected at lower thresholds [28/417 (6.7%) at 10%-25% to 293/812 (36.1%) at 1%-2% threshold]. However, such low-abundance HIV-DRMs were overrepresented in ART-naive participants and were in most cases not detected in previously sampled sequences suggesting high sequencing error for thresholds <3%.

CONCLUSIONS

We found high concordance between SS and NGS but also a substantial number of low-abundance HIV-DRMs detected only by NGS at lower variant-calling thresholds. Our findings suggest that a substantial fraction of the low-abundance HIV-DRMs detected at thresholds <3% may represent sequencing errors and hence should not be overinterpreted in clinical practice.

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.

Genotypic Resistance Testing of HIV-1 DNA in Peripheral Blood Mononuclear Cells

HIV-1 DNA exists in nonintegrated linear and circular episomal forms and as integrated proviruses. In patients with plasma viremia, most peripheral blood mononuclear cell (PBMC) HIV-1 DNA consists of recently produced nonintegrated virus DNA while in patients with prolonged virological suppression (VS) on antiretroviral therapy (ART), most PBMC HIV-1 DNA consists of proviral DNA produced months to years earlier. Drug-resistance mutations (DRMs) in PBMCs are more likely to coexist with ancestral wild-type virus populations than they are in plasma, explaining why next-generation sequencing is particularly useful for the detection of PBMC-associated DRMs. In patients with ongoing high levels of active virus replication, the DRMs detected in PBMCs and in plasma are usually highly concordant. However, in patients with lower levels of virus replication, it may take several months for plasma virus DRMs to reach detectable levels in PBMCs. This time lag explains why, in patients with VS, PBMC genotypic resistance testing (GRT) is less sensitive than historical plasma virus GRT, if previous episodes of virological failure and emergent DRMs were either not prolonged or not associated with high levels of plasma viremia. Despite the increasing use of PBMC GRT in patients with VS, few studies have examined the predictive value of DRMs on the response to a simplified ART regimen. In this review, we summarize what is known about PBMC HIV-1 DNA dynamics, particularly in patients with suppressed plasma viremia, the methods used for PBMC HIV-1 GRT, and the scenarios in which PBMC GRT has been used clinically.

HIV-1 Drug Resistance Assay Using Ion Torrent Next Generation Sequencing and On-Instrument End-to-End Analysis Software

HIV-1 antiretroviral therapy management requires sequencing the protease, reverse transcriptase, and integrase portions of the HIV-1 pol gene. Most resistance testing is performed with Sanger sequencing, which has limited ability to detect minor variants. Next generation sequencing (NGS) platforms enable variant detection at frequencies as low as 1% allowing for earlier detection of resistance and modification of therapy. Implementation of NGS assays in the clinical laboratory is hindered by complicated assay design, cumbersome wet bench procedures, and the complexity of data analysis and bioinformatics. We developed a complete NGS protocol and companion analysis and reporting pipeline using AmpliSeq multiplex PCR, Ion Torrent S5 XL sequencing, and Stanford's HIVdb resistance algorithm. Implemented as a Torrent Suite software plugin, the pipeline runs automatically after sequencing. An optimum variant frequency threshold of 10% was determined by comparing Sanger sequences of archived samples from ViroSeq testing, resulting in a sensitivity of 98.2% and specificity of 99.0%. The majority (91%) of drug resistance mutations were detected by both Sanger and NGS, with 1.7% only by Sanger and 7.3% only by NGS. Variant calls were highly reproducible and there was no cross-reactivity to VZV, HBV, CMV, EBV, and HCV. The limit of detection was 500 copies/mL. The NGS assay performance was comparable to ViroSeq Sanger sequencing and has several advantages, including a publicly available end-to-end analysis and reporting plugin. The assay provides a straightforward path for implementation of NGS for HIV drug resistance testing in the laboratory setting without additional investment in bioinformatics infrastructure and resources.

HIV Drug Resistance Mutations Detection by Next-Generation Sequencing during Antiretroviral Therapy Interruption in China

Patients with antiretroviral therapy interruption have a high risk of virological failure when re-initiating antiretroviral therapy (ART), especially those with HIV drug resistance. Next-generation sequencing may provide close scrutiny on their minority drug resistance variant. A cross-sectional study was conducted in patients with ART interruption in five regions in China in 2016. Through Sanger and next-generation sequencing in parallel, HIV drug resistance was genotyped on their plasma samples. Rates of HIV drug resistance were compared by the McNemar tests. In total, 174 patients were included in this study, with a median 12 (interquartile range (IQR), 6–24) months of ART interruption. Most (86.2%) of them had received efavirenz (EFV)/nevirapine (NVP)-based first-line therapy for a median 16 (IQR, 7–26) months before ART interruption. Sixty-one (35.1%) patients had CRF07_BC HIV-1 strains, 58 (33.3%) CRF08_BC and 35 (20.1%) CRF01_AE. Thirty-four (19.5%) of the 174 patients were detected to harbor HIV drug-resistant variants on Sanger sequencing. Thirty-six (20.7%), 37 (21.3%), 42 (24.1%), 79 (45.4%) and 139 (79.9) patients were identified to have HIV drug resistance by next-generation sequencing at 20% (v.s. Sanger, p = 0.317), 10% (v.s. Sanger, p = 0.180), 5% (v.s. Sanger, p = 0.011), 2% (v.s. Sanger, p < 0.001) and 1% (v.s. Sanger, p < 0.001) of detection thresholds, respectively. K65R was the most common minority mutation, of 95.1% (58/61) and 93.1% (54/58) in CRF07_BC and CRF08_BC, respectively, when compared with 5.7% (2/35) in CRF01_AE (p < 0.001). In 49 patients that followed-up a median 10 months later, HIV drug resistance mutations at >20% frequency such as K103N, M184VI and P225H still existed, but with decreased frequencies. The prevalence of HIV drug resistance in ART interruption was higher than 15% in the survey. Next-generation sequencing was able to detect more minority drug resistance variants than Sanger. There was a sharp increase in minority drug resistance variants when the detection threshold was below 5%.

PANDAA intentionally violates conventional qPCR design to enable durable, mismatch-agnostic detection of highly polymorphic pathogens

Sensitive and reproducible diagnostics are fundamental to containing the spread of existing and emerging pathogens. Despite the reliance of clinical virology on qPCR, technical challenges persist that compromise their reliability for sustainable epidemic containment as sequence instability in probe-binding regions produces false-negative results. We systematically violated canonical qPCR design principles to develop a Pan-Degenerate Amplification and Adaptation (PANDAA), a point mutation assay that mitigates the impact of sequence variation on probe-based qPCR performance. Using HIV-1 as a model system, we optimized and validated PANDAA to detect HIV drug resistance mutations (DRMs). Ultra-degenerate primers with 3' termini overlapping the probe-binding site adapt the target through site-directed mutagenesis during qPCR to replace DRM-proximal sequence variation. PANDAA-quantified DRMs present at frequency ≥5% (2 h from nucleic acid to result) with a sensitivity and specificity of 96.9% and 97.5%, respectively. PANDAA is an innovative advancement with applicability to any pathogen where target-proximal genetic variability hinders diagnostic development.

HIV-1 drug resistance mutations detection and HIV-1 subtype G report by using next-generation sequencing platform

BACKGROUND

Based on world health organization (WHO) recommend, drug resistance assay should be performed in initial of treatment and after treatment for administering and monitoring of anti-retroviral regime in HIV-1 infected patients.

MATERIAL AND METHOD

NGS analyses were performed on forty-one plasma samples from HIV-1 affected patients using the Sentosa SQ HIV genotyping assay (Vela-Diagnostics, Germany). This system comprises a semi-automated Ion torrent based platform and the sequencing results were analyzed based on ANRS, REGA and Stanford drug resistance algorithms. Phylogenetic analysis was analyzed based on https://comet.lih.lu database as well as MEGA5 Software.

RESULTS

Drug resistances were identified in thirty-three samples (80%) out of forty-one samples. The Phylogenetic analysis results showed that CRF-35AD (94%) and subtypes B (2.4%) and G (2.4%) were dominant subtypes in this study. NRTI and NNRTI associated dominant mutations were M184I/V and K103 N.High-level resistance to lamivudine (3 TC) and Emtricitabine (FTC) were detected in 34.3% of patients while 53.1% were resistant to Efavirenz (EFV) and Nevirapine (NVP). The Protease inhibitor (PI) minor and major mutations were not reported but more than 95% of samples had polymorphisms mutation in K20R, M36I, H69K, L89 M positions. These mutations are subtype dependent and completely are absent in subtype B virus. The secondary mutations were reported in positions of E157Q, S230 N, and T97A of integrase gene and four samples represent low-level resistance to integrase strand transfer inhibitor (INSTI).

CONCLUSIONS

This is the first preliminary evaluation of HIV-1 drug resistance mutation (DRM) by using the Sentosa SQ HIV Genotyping Assay in Iran. The NGS represent a promising tool for the accurate detection of DRMs of CRF-35AD that is dominant subtype in Iranian HIV-1 infected population and for the first time revealed HIV-1 subtype G in Iranian population. In the present study polymorphic mutation in the position of K20R, M36I, H69K, L89 M were properly reported in CRF35AD that is dominant in Iranian HIV patients.

Shall I trust the report? Variable performance of Sanger sequencing revealed by deep sequencing on HIV drug resistance mutation detection

BACKGROUND

The clinical utilisation of deep sequencing in HIV treatment has been hindered due to its unknown correlation with standard Sanger genotyping and the undetermined value of minority drug resistance mutation (DRM) detection.

OBJECTIVES

To compare deep sequencing performance to standard Sanger genotyping with clinical samples, in an effort to delineate the correlation between the results from the two methods and to find the optimal deep sequencing threshold for clinical utilisation.

METHODS

We conducted a retrospective study using stored plasma collected from August 2014 to March 2018 for HIV genotyping with the commercial Sanger genotyping kit. Samples with available Sanger genotyping reports were further deep sequenced. Drug resistance was interpreted according to the Stanford HIV drug resistance database algorithm.

RESULTS

At 15-25% minority detection thresholds, 9-15% cases had underestimated DRMs by Sanger sequencing. The concordance between the Sanger and deep sequencing reports was 68-82% in protease-reverse transcriptase region and 88-97% in integrase region at 5-25% thresholds. The undetected drug resistant minority variants by Sanger sequencing contributed to the lower negative predictive value of Sanger genotyping in cases harbouring DRMs.

CONCLUSIONS

Use of deep sequencing improved detection of antiretroviral resistance mutations especially in cases with virological failure or previous treatment interruption. Deep sequencing with 10-15% detection thresholds may be considered a suitable substitute for Sanger sequencing on antiretroviral DRM detection.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Increased predominance of HIV-1 CRF01_AE and its recombinants in the Philippines

The growth rate of new HIV infections in the Philippines was the fastest of any countries in the Asia-Pacific region between 2010 and 2016. To date, HIV-1 subtyping results in the Philippines have been determined by characterizing only partial viral genome sequences. It is not known whether recombination occurs in the majority of unsequenced genome regions. Near-full-length genome (NFLG) sequences were obtained by amplifying two overlapping half genomes from plasma samples collected between 2015 and 2017 from 23 newly diagnosed infected individuals in the Philippines. Phylogenetic analysis showed that the newly characterized sequences were CRF01_AE (14), subtype B (3), CRF01/B recombinants (5) and a CRF01/CRF07/B recombinant (1). All 14 CRF01_AE formed a tight cluster, suggesting that they were derived from a single introduction. The time to the most recent common ancestor (tMRCA) for CRF01_AE in the Philippines was 1995 (1992-1998), about 10-15 years later than that of CRF01_AE in China and Thailand. All five CRF01/B recombinants showed distinct recombination patterns, suggesting ongoing recombination between the two predominant circulating viruses. The identification of partial CRF07_BC sequences in one CRF01/CRF07/B recombinant, not reported previously in the Philippines, indicated that CRF07_BC may have been recently introduced into that country from China, where CRF07_BC is prevalent. Our results show that the major epidemic strains may have shifted to an increased predominance of CRF01_AE and its recombinants, and that other genotypes such as CRF07_BC may have been introduced into the Philippines.

First evaluation of the Next-Generation Sequencing platform for the detection of HIV-1 drug resistance mutations in Belgium

Introduction

The WHO urges action against the threat posed by HIV drug resistance. It is well known that the sensitivity of Next-Generation Sequencing (NGS) is greater than that of Sanger Sequencing (SS). The objective of this study was to evaluate the performance of the novel NGS HIV-1 drug resistance monitoring system.

Materials & methods

NGS analyses were performed on 67 plasma samples from HIV-1 infected patients using the Sentosa SQ HIV Genotyping Assay from Vela-Dx. This kit was used on a semi-automated Ion Torrent-based platform. Sequences were compared to those obtained by SS. Samples were analysed in the same and in separate runs. Quality controls (QC) were added to control sequencing processes of protease (PRO), reverse transcriptase (RT) and integrase (INT) regions.

Results

Of the 41 analysed samples, 33 (80.5%) had identical drug resistance interpretation reports. Discrepant results were observed for eight samples. Five of them were only detected by NGS and had drug resistance mutations (DRMs) with an allelic frequency below the limit of detection of the SS method (between 6.3 to 20.5%). Two DRMs were only identified using the SS method. The sequences were similar in 98.2% of cases (counting variants as mismatches) and homologous in 99.9% if missed variants. Duplicated samples in a single run were similar in 95.7% (99.9%) of cases. Duplicated samples in two different runs were 98% (100%) homologous. QC results were manually assessed with a score of 340/340 for detection of DRMs in PRO and RT and 100% for INT sequencing.

Conclusions

This is the first preliminary evaluation in Belgium employing the Sentosa SQ HIV Genotyping Assay. The NGS appears to be a promising tool for the detection of DRMs in HIV-1 patients and showed a higher sensitivity compared to SS. Large studies assessing the clinical relevance of low frequency DRMs are needed.

Comparison of an In Vitro Diagnostic Next-Generation Sequencing Assay with Sanger Sequencing for HIV-1 Genotypic Resistance Testing

The ability of next-generation sequencing (NGS) technologies to detect low frequency HIV-1 drug resistance mutations (DRMs) not detected by dideoxynucleotide Sanger sequencing has potential advantages for improved patient outcomes. We compared the performance of an in vitro diagnostic (IVD) NGS assay, the Sentosa SQ HIV genotyping assay for HIV-1 genotypic resistance testing, with Sanger sequencing on 138 protease/reverse transcriptase (RT) and 39 integrase sequences. The NGS assay used a 5% threshold for reporting low-frequency variants. The level of complete plus partial nucleotide sequence concordance between Sanger sequencing and NGS was 99.9%. Among the 138 protease/RT sequences, a mean of 6.4 DRMs was identified by both Sanger and NGS, a mean of 0.5 DRM was detected by NGS alone, and a mean of 0.1 DRM was detected by Sanger sequencing alone. Among the 39 integrase sequences, a mean of 1.6 DRMs was detected by both Sanger sequencing and NGS and a mean of 0.15 DRM was detected by NGS alone. Compared with Sanger sequencing, NGS estimated higher levels of resistance to one or more antiretroviral drugs for 18.2% of protease/RT sequences and 5.1% of integrase sequences. There was little evidence for technical artifacts in the NGS sequences, but the G-to-A hypermutation was detected in three samples. In conclusion, the IVD NGS assay evaluated in this study was highly concordant with Sanger sequencing. At the 5% threshold for reporting minority variants, NGS appeared to attain a modestly increased sensitivity for detecting low-frequency DRMs without compromising sequence accuracy.