Deep sequencing for HIV-1 clinical management

Next-Generation Sequencing Methods for Near-Real-Time Molecular Epidemiology of HIV and HCV

The World Health Organisation has set targets of reducing the transmission of new hepatitis C (HCV) infections by 90%, and ending human immunodeficiency virus-1 (HIV) as a public health threat, by 2030. To achieve this, efficient and timely viral surveillance, and effective public health interventions, are required. Traditional epidemiological methods are largely dependent on the recognition of incident cases with symptomatic illness; acute HIV and HCV infections are commonly asymptomatic, which may lead to delays in the recognition of such new infections. Instead, for these viruses, molecular epidemiology may improve the detection of, and response to, clusters of viral transmission. Molecular epidemiology using historical datasets has highlighted key populations that may have benefitted from a timely intervention. Similar analyses performed on contemporary samples are needed to underpin the 2030 targets, but this requires the generation of a cohesive dataset of viral genome sequences in near-real-time. To generate such data, methodologies harnessing next-generation sequencing (NGS) should be utilised. Here we discuss the opportunity presented by NGS for public health surveillance of HIV and HCV, and discuss three methods that can generate sequences for such analysis. These include full-length genome amplification, utilised for analysis of HCV in the research space; tiling PCR, which was the method of choice for many diagnostic laboratories in the SARS-CoV-2 pandemic; and bait-capture hybridisation, which has been utilised in local HIV outbreaks. These techniques could be applied for near-real-time HIV and HCV surveillance, informing public health strategies that will be key to achieving 2030 targets.

Comparison of Different HIV-1 Resistance Interpretation Tools for Next-Generation Sequencing in Italy

BACKGROUND

Next-generation sequencing (NGS) is gradually replacing Sanger sequencing for HIV genotypic drug resistance testing (GRT). This work evaluated the concordance among different NGS-GRT interpretation tools in a real-life setting.

METHODS

Routine NGS-GRT data were generated from viral RNA at 11 Italian laboratories with the AD4SEQ HIV-1 Solution v2 commercial kit. NGS results were interpreted by the SmartVir system provided by the kit and by two online tools (HyDRA Web and Stanford HIVdb). NGS-GRT was considered valid when the coverage was >100 reads (100×) at each PR/RT/IN resistance-associated position listed in the HIVdb 9.5.1 algorithm.

RESULTS

Among 629 NGS-GRT, 75.2%, 74.2%, and 70.9% were valid according to SmartVir, HyDRA Web, and HIVdb. Considering at least two interpretation tools, 463 (73.6%) NGS-GRT had a valid coverage for resistance analyses. The proportion of valid samples was affected by viremia <10,000-1000 copies/mL and non-B subtypes. Mutations at an NGS frequency >10% showed fair concordance among different interpretation tools.

CONCLUSION

This Italian survey on NGS resistance testing suggests that viremia levels and HIV subtype affect NGS-GRT coverage. Within the current routine method for NGS-GRT, only mutations with frequency >10% seem reliably detected across different interpretation tools.

The relevance of ultradeep sequencing for low HIV-1 viral loads and proviruses in the clinical setting

Improving the therapeutic management of HIV-positive persons is a major public health issue and includes better detection of drug resistance mutations (DRMs). The aim of this study was (i) to explore DRMs in HIV-1-positive persons presenting a blood viral load (VL) < 1000 genomes copies (gc)/mL, including the analyze of cerebrospinal fluid (CSF) and HIV-DNA from peripheral blood mononuclear cells using ultradeep sequencing (UDS) and (ii), to evaluate how these DRMs could influence the clinical practices. For each patient (n = 12), including five low-VL patients (i.e., <1000 gc/mL), HIV-1 UDS targeting the protease, reverse transcriptase and integrase genes was performed on plasma, proviral DNA, and CSF when available. Sequencing discordances or failures were mostly found in samples from low-VL patients. A 5% UDS cut-off allowed to increase the sensitivity to detect DRMs in different compartments, excepted in CSF. Patients with the highest viral quasispecies heterogeneity were naïve of treatment or presented a medical history suggesting low selection pressure or virological failures. When analyzing compartmentalization and following-up patients: low-frequency variants (LFVs) were responsible for 47% (n = 8) and 76% (n = 13) of changes in drug resistance interpretation, respectively. In such cases, we conclude that UDS is a robust technique, which still could be improved by increase the RNA and/or DNA extraction in low-VL samples to detect LFVs. Further studies are needed to define the impact of LFVs on antiretroviral treatments. At last, when considering a DRM, the use of mutational load would probably be more suitable than frequencies.

Full-Spectrum Surveillance of Pre-Treatment HIV Drug Resistance in Southeastern China

HIV drug resistance compromises the ability of anti-retroviral therapy (ART) to suppress viral replication, resulting in treatment failure. This study investigates the prevalence of pre-treatment drug resistance (PDR) in newly diagnosed individuals in a prosperous city (Wenzhou) in Southeastern China. A cross-sectional investigation was carried out among 473 newly diagnosed ART-naive HIV-1-infected individuals between January and December 2022. The protease-reverse transcriptase (PR-RT) region and integrase (IN) region of HIV-1 were amplified by two separately nested PCRs, followed by sequencing. Drug resistance mutations (DRMs) and drug resistance to nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), protease inhibitors (PIs) and integrase strand transfer inhibitors (INSTIs) were analyzed. The PDR prevalence was 6.5% [95% CI: 4.4-9.1] for any anti-retroviral drug, 0.9% [95% CI: 0.3-2.3] for NRTIs, 4.1% [95% CI: 2.5-6.5] for NNRTIs, 1.8% [95% CI: 0.8-3.6] for PIs and 0.5% [95% CI: 0.1-1.8] for INSTIs. According to the subtyping results of the PR-RT region, 11 different subtypes and 31 unique recombinant forms (URFs) were found. CRF07_BC was the dominant subtype (53.7%, 233/434), followed by CRF01_AE (25.3%, 110/434). V179D (1.6%) and K103N (1.4%) were the most predominant types of NNRTI DRMs. Q58E (1.2%) and M184V (0.7%) were the most frequent PI DRMs and NRTI DRMs, respectively. The INSTI-related DRMs Y143S (causes high-level resistance to RAL) and G163K (causes low-level resistance to EVG and RAL) were found in one patient each. Given the relatively high PDR prevalence of NNRTI (4.1%), non-NNRTI-based ART may be preferred in the future. It is recommended to include genotypic resistance testing before starting ART in regions where feasible.

Validation of an HIV whole genome sequencing method for HIV drug resistance testing in an Australian clinical microbiology laboratory

Detection of HIV drug resistance (HIVDR) is vital to successful anti-retroviral therapy (ART). HIVDR testing to determine drug-resistance mutations is routinely performed in Australia to guide ART choice in newly diagnosed people living with HIV or in cases of treatment failure. In 2022, our clinical microbiology laboratory sought to validate a next-generation sequencing (NGS)-based HIVDR assay to replace the previous Sanger-sequencing (SS)-based ViroSeq. NGS solutions for HIVDR offer higher throughput, lower costs and higher sensitivity for variant detection. We sought to validate the previously described low-cost probe-based NGS method (veSEQ-HIV) for whole-genome recovery and HIVDR-testing in a diagnostic setting. veSEQ-HIV displayed 100% and 98% accuracy in major and minor mutation detection, respectively, and 100% accuracy of subtyping (provided > 1000 mapped reads were obtained). Pairwise comparison exhibited low inter-and intrarun variability across the whole-genome (Jaccard index [J] = 0.993; J = 0.972) and the Pol gene (J = 0.999; J = 0.999), respectively. veSEQ-HIV met all our pre-set criteria based on WHO recommendations and successfully replaced ViroSeq in our laboratory. Scaling-down veSEQ-HIV to a limited batch size and sequencing on Illumina iSeq. 100, allowed easy implementation of the assay into the workflow of a small sequencing laboratory with minimal staff and equipment and the ability to meet clinically relevant test turn-around times. As HIVDR-testing moves from SS- to NGS-based methods and new ART drugs come to market (particularly those with targets outside the Pol region), whole-genome recovery using veSEQ-HIV provides a robust, cost-effective and "future-proof" NGS method for HIVDR-testing.

Impact of pretreatment low-abundance HIV-1 drug resistance on virological failure after 1 year of antiretroviral therapy in China

OBJECTIVES

To assess the impact of pretreatment low-abundance HIV drug-resistant variants (LA-DRVs) on virological outcomes among ART-naive HIV-1-infected Chinese people who initiated ART.

METHODS

A nested case-control study was conducted among HIV-1-infected individuals who had pretreatment drug resistance (PDR) genotypic results. Cases were defined as individuals with virological failure (HIV-1 RNA viral load ≥1000 copies/mL) after 1 year of ART, and controls were individuals from the same cohort whose viral load was less than 1000 copies/mL. Next-generation sequencing was used to identify low-abundance PDR mutations at detection thresholds of 10%, 2% and 1%. The mutant load was calculated by multiplying the abundance of HIV-1 drug-resistant variants by the pretreatment viral load. The impact of pretreatment low-abundance mutations on virological failure was estimated in logistic regression models.

RESULTS

Participants (43 cases and 100 controls) were included in this study for the analysis. The proportion of participants with PDR was higher in cases than in controls at different detection thresholds (44.2% versus 22.0%, P = 0.007 at 10% threshold; 58.1% versus 31.0%, P = 0.002 at 2% threshold; 90.7% versus 69.0%, P = 0.006 at 1% threshold). Compared with participants without PDR, participants with ≥10% detectable PDR mutations were associated with an increased risk of virological failure (adjusted OR 8.0, 95% CI 2.4-26.3, P = 0.001). Besides this, individuals with pretreatment LA-DRVs (2%-9% abundance range) had 5-fold higher odds of virological failure (adjusted OR 5.0, 95% CI 1.3-19.6, P = 0.021). Furthermore, LA-DRVs at 2%-9% abundance resistant to NRTIs and mutants with abundance of ≥10% resistant to NNRTIs had a 4-fold and 8-fold risk of experiencing virological failure, respectively. It was also found that a mutant load of more than 1000 copies/mL was predictive of virological failure (adjusted OR 7.2, 95% CI 2.5-21.1, P = 0.0003).

CONCLUSIONS

Low-abundance PDR mutations ranging from 2% to 9% of abundance can increase the risk of virological failure. Further studies are warranted to define a clinically relevant threshold of LA-DRVs and the role of NRTI LA-DRVs.

Sensitive HIV-1 DNA Pol Next-Generation Sequencing for the Characterisation of Archived Antiretroviral Drug Resistance

Modern HIV-1 treatment effectively suppresses viral amplification in people living with HIV. However, the persistence of HIV-1 DNA as proviruses integrated into the human genome remains the main barrier to achieving a cure. Next-generation sequencing (NGS) offers increased sensitivity for characterising archived drug resistance mutations (DRMs) in HIV-1 DNA for improved treatment options. In this study, we present an ultra-sensitive targeted PCR assay coupled with NGS and a robust pipeline to characterise HIV-1 DNA DRMs from buffy coat samples. Our evaluation supports the use of this assay for Pan-HIV-1 analyses with reliable detection of DRMs across the HIV-1 Pol region. We propose this assay as a new valuable tool for monitoring archived HIV-1 drug resistance in virologically suppressed individuals, especially in clinical trials investigating novel therapeutic approaches.

Molecular transmission network of pretreatment drug resistance among human immunodeficiency virus-positive individuals and the impact of virological failure on those who received antiretroviral therapy in China

Objectives

We investigated the prevalence of pretreatment drug resistance (PDR), the molecular transmission network among HIV-positive individuals, and the impact of virological failure on those who received antiretroviral therapy (ART) in China.

Methods

Based on the World Health Organization (WHO) surveillance guidelines for PDR, a baseline survey and follow-up were conducted in 2018 and 2021, respectively. Demographic information and plasma samples were obtained from all participants. HIV pol gene region sequences were used to analyze the PDR and molecular transmission networks using the Stanford HIV database algorithm and HIV-TRACE, respectively. This study assessed the odds ratios (OR) of PDR to virological failure (viral load ≥ 50 copies/mL) after 3 years of ART using multivariable logistic regression.

Results

Of the 4,084 individuals, 370 (9.1%) had PDR. The prevalence of PDR to non-nucleoside reverse transcriptase inhibitors (5.2%) was notably higher than that to nucleoside reverse transcriptase inhibitors (0.7%, p < 0.001), protease inhibitors (3.0%, p < 0.001), and multidrug resistance (0.3%, p < 0.001). A total of 1,339 (32.8%) individuals from 361 clusters were enrolled in the molecular transmission network. Of the 361 clusters, 22 included two or more individuals with PDR. The prevalence of virological failure among HIV-positive individuals after 3 years of ART without PDR, those with PDR to Chinese listed drugs, and those with PDR to other drugs was 7.9, 14.3, and 12.6%, respectively. Compared with that in HIV-positive individuals without PDR, virological failure after 3 years of ART was significantly higher (OR: 2.02, 95% confidence interval (CI): 1.25–3.27) and not significantly different (OR: 1.72, 95% CI: 0.87–3.43) in individuals with PDR to Chinese listed drugs and those with PDR to other drugs, respectively. Missed doses in the past month were significantly associated with virological failure (OR, 2.82; 95% CI: 4.08–5.89).

Conclusion

The overall prevalence of PDR was close to a high level and had an impact on virological failure after 3 years of ART. Moreover, HIV drug-resistant strains were transmitted in the molecular transmission network. These results illustrate the importance of monitoring PDR and ensuring virological suppression through drug adherence.

Low-frequency HIV-1 drug resistance mutations in antiretroviral naïve individuals in Botswana

BACKGROUND

Individuals living with human immunodeficiency virus (HIV) who experience virological failure (VF) after combination antiretroviral therapy (cART) initiation may have had low-frequency drug resistance mutations (DRMs) at cART initiation. There are no data on low-frequency DRMs among cART-naïve HIV-positive individuals in Botswana.

METHODS

We evaluated the prevalence of low-frequency DRMs among cART-naïve individuals previously sequenced using Sanger sequencing. The generated pol amplicons were sequenced by next-generation sequencing.

RESULTS

We observed low-frequency DRMs (detected at <20% in 33/103 (32%) of the successfully sequenced individuals, of whom four also had mutations detected at >20%. K65R was the most common low-frequency DRM detected in 8 individuals. Eighty-two of the 103 individuals had follow-up viral load data while on cART. Twenty-seven of the 82 individuals harbored low-frequency DRMs. Only 12 of 82 individuals experienced VF. The following low-frequency DRMs were observed in four individuals experiencing VF: K65R, K103N, V108I, and Y188C. No statistically significant difference was observed in the prevalence of low-frequency DRMs between individuals experiencing VF (4/12) and those not experiencing VF (23/70) (P = .97). However, individuals with non-nucleoside reverse transcriptase inhibitors-associated low-frequency DRMs were 2.68 times more likely to experience VF (odds ratio, 2.68; 95% confidential interval, 0.4-13.9) compared with those without (P = .22).

CONCLUSION

Next-generation sequencing was able to detect low-frequency DRMs in this cohort in Botswana, but these DRMs did not contribute significantly to VF.

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.

Advances in multiplex molecular detection technologies for harmful algae

As the eutrophication of natural water bodies becomes more and more serious, the frequency of outbreaks of harmful algal blooms (HABs) mainly formed by harmful algae also increases. HABs have become a global ecological problem that poses a serious threat to human health and food safety. Therefore, it is extremely important to establish methods that can rapidly detect harmful algal species for early warning of HABs. The traditional morphology-based identification method is inefficient and inaccurate. In recent years, the rapid development of molecular biology techniques has provided new ideas for the detection of harmful algae and has become a research hotspot. The current molecular detection methods for harmful algal species mainly include fluorescence in situ hybridization, sandwich hybridization, and quantitative PCR (qPCR), but all of these methods can only detect single harmful algal species at a time. The establishment of methods for the simultaneous detection of multiple harmful algal species has become a new trend in the development of molecular detection technology because various harmful algal species may coexist in the natural water environment. The established molecular techniques for multiple detections of harmful algae mainly include gene chip, multiplex PCR, multiplex qPCR, massively parallel sequencing, antibody chip, and multiple isothermal amplification. This review mainly focuses on the principles, advantages and disadvantages, application progress, and application prospects of these multiple detection technologies, aiming at providing effective references not only for the fisheries but also for economic activities, environment, and human health.

Detection of archived lamivudine-associated resistance mutations in virologically suppressed, lamivudine-experienced HIV-infected adults by different genotyping techniques (GEN-PRO study)

BACKGROUND

Previously selected lamivudine resistance-associated mutations (RAMs) may remain archived within the proviral HIV-DNA.

OBJECTIVES

To evaluate the ability of proviral DNA genotyping to detect lamivudine RAMs in HIV-1 virologically suppressed participants; the correlation between Sanger and next generation sequencing (NGS); and predictive factors for detection of lamivudine RAMs in proviral DNA.

METHODS

Cross-sectional study of participants on stable antiretroviral therapy and suppressed for ≥1 year. Analysis of proviral DNA was performed by Sanger sequencing in whole blood and by NGS in PBMCs.

RESULTS

We analysed samples from 102 subjects (52 with and 50 without lamivudine RAMs in historical plasma RNA-genotypes). Among participants with previous lamivudine resistance, Sanger sequencing detected RAMs in 26.9%. Detection rates significantly increased using NGS: 47.9%, 64.6%, 75% and 87.5% with the 20%, 10%, 5% and 1% thresholds, respectively. As for participants without historical lamivudine resistance, Sanger detected the RAMs in 1/49 (2%), and NGS (5% threshold) in 8/45 (17.8%). Multivariate models fitted to the whole population revealed that having a history of lamivudine resistance was a risk factor for detection of lamivudine RAMs by NGS. Among participants with historical lamivudine resistance, multivariate analysis showed that a longer time since HIV diagnosis was associated with persistence of archived mutations by NGS at thresholds of >10% [OR 1.10 (95% CI: 1.00-1.24)] and >5% [OR 1.16 (95% CI: 1.02-1.32)].

CONCLUSIONS

Proviral DNA Sanger sequencing does not detect the majority of historical lamivudine RAMs. NGS increases the sensitivity of detection at lower thresholds, although the relevance of these minority populations with lamivudine RAMs needs further evaluation.

Next generation sequencing based in-house HIV genotyping method: validation report

Background

HIV genotyping has had a significant impact on the care and treatment of HIV/AIDS. At a clinical level, the test guides physicians on the choice of treatment regimens. At the surveillance level, it informs policy on consolidated treatment guidelines and microbial resistance control strategies. Until recently, the conventional test has utilized the Sanger sequencing (SS) method. Unlike Next Generation Sequencing (NGS), SS is limited by low data throughput and the inability of detecting low abundant drug-resistant variants. NGS can improve sensitivity and quantitatively identify low-abundance variants; in addition, it has the potential to improve efficiency as well as lowering costs when samples are batched. Despite the NGS benefits, its utilization in clinical drug resistance profiling is faced with mixed reactions. These are largely based on a lack of a consensus regarding the quality control strategy. Nonetheless, transitional views suggest validating the method against the gold-standard SS. Therefore, we present a validation report of an NGS-based in-house HIV genotyping method against the SS method in Uganda.

Results

Since there were no established proficiency test panels for NGS-based HIV genotyping, 15 clinical plasma samples for routine care were utilized. The use of clinical samples allowed for accuracy and precision studies. The workflow involved four main steps; viral RNA extraction, targeted amplicon generation, amplicon sequencing and data analysis. Accuracy of 98% with an average percentage error of 3% was reported for the NGS based assay against the SS platform demonstrating similar performance. The coefficient of variation (CV) findings for both the inter-run and inter-personnel precision showed no variability (CV ≤ 0%) at the relative abundance of ≥ 20%. For both inter-run and inter-personnel, a variation that affected the precision was observed at 1% frequency. Overall, for all the frequencies, CV registered a small range of (0–2%).

Conclusion

The NGS-based in-house HIV genotyping method fulfilled the minimum requirements that support its utilization for drug resistance profiling in a clinical setting of a low-income country. For more inclusive quality control studies, well-characterized wet panels need to be established.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12981-021-00390-8.

Primary resistance to integrase strand transfer inhibitors in Spain using ultrasensitive HIV-1 genotyping

BACKGROUND

Transmission of resistance mutations to integrase strand transfer inhibitors (INSTIs) in HIV-infected patients may compromise the efficacy of first-line antiretroviral regimens currently recommended worldwide. Continued surveillance of transmitted drug resistance (TDR) is thus warranted.

OBJECTIVES

We evaluated the rates and effects on virological outcomes of TDR in a 96 week prospective multicentre cohort study of ART-naive HIV-1-infected subjects initiating INSTI-based ART in Spain between April 2015 and December 2016.

METHODS

Pre-ART plasma samples were genotyped for integrase, protease and reverse transcriptase resistance using Sanger population sequencing or MiSeq™ using a ≥ 20% mutant sensitivity cut-off. Those present at 1%-19% of the virus population were considered to be low-frequency variants.

RESULTS

From a total of 214 available samples, 173 (80.8%), 210 (98.1%) and 214 (100.0%) were successfully amplified for integrase, reverse transcriptase and protease genes, respectively. Using a Sanger-like cut-off, the overall prevalence of any TDR, INSTI-, NRTI-, NNRTI- and protease inhibitor (PI)-associated mutations was 13.1%, 1.7%, 3.8%, 7.1% and 0.9%, respectively. Only three (1.7%) subjects had INSTI TDR (R263K, E138K and G163R), while minority variants with integrase TDR were detected in 9.6% of subjects. There were no virological failures during 96 weeks of follow-up in subjects harbouring TDR as majority variants.

CONCLUSIONS

Transmitted INSTI resistance remains rare in Spain and, to date, is not associated with virological failure to first-line INSTI-based regimens.

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%.

High-throughput sequencing reveals a high prevalence of pretreatment HIV-1 drug resistance in Sweden

OBJECTIVES

HIV-1 pretreatment drug resistance (PDR) is a global concern. Our aim was to evaluate high-throughput sequencing (HTS) for HIV-1 resistance testing and describe PDR in Sweden, where 75% of diagnosed individuals are foreign-born.

DESIGN

Cross-sectional study.

METHODS

Individuals entering HIV-1 care in Sweden 2017 to March 2019 (n = 400) were included if a viremic sample was available (n = 220). HTS was performed using an in-house assay. Drug resistance mutations (DRMs) (based on Stanford HIV DB vs. 8.7) at levels 1-5%, 5-19% and at least 20% of the viral population were described. Results from HTS and routine Sanger sequencing were compared.

RESULTS

HTS was successful in 88% of patients, 92% when viral load was at least 1000 copies/ml. DRMs at any level in protease and/or reverse transcriptase were detected in 95 individuals (49%), whereas DRMs at least 20% in 35 (18%) individuals. DRMs at least 20% correlated well to findings in routine Sanger sequencing. Protease/reverse transcriptase (PR/RT) DRMs at least 20% were predicted by treatment exposure; adjusted OR 9.28 (95% CI 2.24-38.43; P = 0.002) and origin in Asia; adjusted OR 20.65 (95% CI 1.66-256.24; P = 0.02). Nonnucleoside reverse transcriptase inhibitor (NNRTI) DRMs at least 20% were common (16%) and over-represented in individuals originating from sub-Saharan Africa or Asia. Low-level integrase strand transfer inhibitor (INSTI) DRMs less than 20% were detected in 15 individuals (8%) with no association with INSTI exposure.

CONCLUSION

Our HTS can efficiently detect PDR and findings of DRMs at least 20% compare well to routine Sanger sequencing. The high prevalence of PDR was because of NNRTI DRMs and associated with migration from areas with emerging PDR.

Enhancing care for people living with HIV: current and future monitoring approaches

INTRODUCTION

Antiretroviral therapy (ART) is the most significant advance in the medical management of HIV-1 infection. Given the fact that HIV cannot be eradicated from the body, ART has to be indefinitely maintained. New approaches need to be defined for monitoring HIV-infected individuals (PLWHIV), including clinical, virologic, immunological parameters and also ways to collect individual points of view and quality of life.

AREAS COVERED

We discuss which tests may be used to improve the management of PLWHIV and respond to a comprehensive health demand.

EXPERT OPINION

Viral load and CD4 counts are well-validated outcome measures and we still need them, but they do not completely depict the health status of PLWHIV. We need to better understand and to apply to clinical practice what happens in sanctuaries, what is the role of HIV DNA, what is the meaning of low-level viremia. Most of these questions do not yet have a definitive response. Further, we need to understand how to modify these variables in order to improve outcomes.Similar points may be raised for immunological measures and for tests exploring the tolerability of drugs. The goal must be the evolution from a viro/immunologic-based to a comprehensive quality-of-health-based evaluation of PLWHIV.

Near point-of-care, point-mutation test to detect drug resistance in HIV-1: a validation study in a Mexican cohort

OBJECTIVE

Pretreatment HIV-drug resistance (PDR, HIVDR) to non-nucleoside reverse transcriptase inhibitors (NNRTIs) is increasing globally. NNRTIs continue to be used as first-line antiretroviral therapy (ART) in some communities due to the cost of dolutegravir-based ART or dolutegravir-associated adverse events. A simplified version of the oligonucleotide ligation assay (OLA) - 'OLA-Simple' - is a low-cost, near point-of-care assay that provides ready-to-use lyophilized reagents and reports HIVDR mutations as colored lines on lateral flow strips. Our objective was to design and validate OLA-Simple for a Mexican cohort.

DESIGN

OLA-Simple probes to detect K65R, K103N/S, Y181C, M184V, and G190A were optimized for HIV Mexican sequences. Sixty clinical plasma specimens were analyzed by OLA-Simple by technicians blinded to Illumina-MiSeq sequences, and HIVDR results were compared.

METHODS

Plasma RNA was tested using OLA-Simple kits. OLA-Simple lateral flow strips were read by in-house software, and were classified as mutant or wild-type at each codon. The comparison of results by OLA-Simple and Miseq was used to generate receiver-operating characteristic curves.

RESULTS

OLA-Simple PCR amplified 59 of 60 specimens and successfully genotyped 287 of 295 codons, with eight of 295 (2.7%) indeterminate results. Compared to MiSeq, OLA-Simple gave five of 295 (1.7%) false-positive and four of 295 (1.4%) false-negative results. Excluding indeterminate results, OLA-Simple classified mutant with an accuracy of 97.4 and 98.8% when using thresholds at 10 and 25% mutant within an individual's HIV quasispecies, respectively.

CONCLUSIONS

Compared to MiSeq, OLA-Simple detected HIVDR with high sensitivity and accuracy. OLA-Simple could expand access to affordable and rapid HIVDR testing to guide appropriate ART choices in populations using NNRTI-based ART.

Multi-Laboratory Comparison of Next-Generation to Sanger-Based Sequencing for HIV-1 Drug Resistance Genotyping

Next-generation sequencing (NGS) is increasingly used for HIV-1 drug resistance genotyping. NGS methods have the potential for a more sensitive detection of low-abundance variants (LAV) compared to standard Sanger sequencing (SS) methods. A standardized threshold for reporting LAV that generates data comparable to those derived from SS is needed to allow for the comparability of data from laboratories using NGS and SS. Ten HIV-1 specimens were tested in ten laboratories using Illumina MiSeq-based methods. The consensus sequences for each specimen using LAV thresholds of 5%, 10%, 15%, and 20% were compared to each other and to the consensus of the SS sequences (protease 4-99; reverse transcriptase 38-247). The concordance among laboratories' sequences at different thresholds was evaluated by pairwise sequence comparisons. NGS sequences generated using the 20% threshold were the most similar to the SS consensus (average 99.6% identity, range 96.1-100%), compared to 15% (99.4%, 88.5-100%), 10% (99.2%, 87.4-100%), or 5% (98.5%, 86.4-100%). The average sequence identity between laboratories using thresholds of 20%, 15%, 10%, and 5% was 99.1%, 98.7%, 98.3%, and 97.3%, respectively. Using the 20% threshold, we observed an excellent agreement between NGS and SS, but significant differences at lower thresholds. Understanding how variation in NGS methods influences sequence quality is essential for NGS-based HIV-1 drug resistance genotyping.

Dry Panels Supporting External Quality Assessment Programs for Next Generation Sequencing-Based HIV Drug Resistance Testing

External quality assessment (EQA) is a keystone element in the validation and implementation of next generation sequencing (NGS)-based HIV drug resistance testing (DRT). Software validation and evaluation is a critical element in NGS EQA programs. While the development, sharing, and adoption of wet lab protocols is coupled with the increasing access to NGS technology worldwide, rendering it easy to produce NGS data for HIV-DRT, bioinformatic data analysis remains a bottleneck for most of the diagnostic laboratories. Several computational tools have been made available, via free or commercial sources, to automate the conversion of raw NGS data into an actionable clinical report. Although different software platforms yield equivalent results when identical raw NGS datasets are analyzed for variations at higher abundance, discrepancies arise when variations at lower frequencies are considered. This implies that validation and performance assessment of the bioinformatics tools applied in NGS HIV-DRT is critical, and the origins of the observed discrepancies should be determined. Well-characterized reference NGS datasets with ground truth on the genotype composition at all examined loci and the exact frequencies of HIV variations they may harbor, so-called dry panels, would be essential in such cases. The strategic design and construction of such panels are challenging but imperative tasks in support of EQA programs for NGS-based HIV-DRT and the validation of relevant bioinformatics tools. Here, we present criteria that can guide the design of such dry panels, which were discussed in the Second International Winnipeg Symposium themed for EQA strategies for NGS HIVDR assays.

Measurement error and variant-calling in deep Illumina sequencing of HIV

MOTIVATION

Next-generation deep sequencing of viral genomes, particularly on the Illumina platform, is increasingly applied in HIV research. Yet, there is no standard protocol or method used by the research community to account for measurement errors that arise during sample preparation and sequencing. Correctly calling high and low-frequency variants while controlling for erroneous variants is an important precursor to downstream interpretation, such as studying the emergence of HIV drug-resistance mutations, which in turn has clinical applications and can improve patient care.

RESULTS

We developed a new variant-calling pipeline, hivmmer, for Illumina sequences from HIV viral genomes. First, we validated hivmmer by comparing it to other variant-calling pipelines on real HIV plasmid datasets. We found that hivmmer achieves a lower rate of erroneous variants, and that all methods agree on the frequency of correctly called variants. Next, we compared the methods on an HIV plasmid dataset that was sequenced using Primer ID, an amplicon-tagging protocol, which is designed to reduce errors and amplification bias during library preparation. We show that the Primer ID consensus exhibits fewer erroneous variants compared to the variant-calling pipelines, and that hivmmer more closely approaches this low error rate compared to the other pipelines. The frequency estimates from the Primer ID consensus do not differ significantly from those of the variant-calling pipelines.

AVAILABILITY AND IMPLEMENTATION

hivmmer is freely available for non-commercial use from https://github.com/kantorlab/hivmmer.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Next-Generation Sequencing for HIV Drug Resistance Testing: Laboratory, Clinical, and Implementation Considerations

Higher accessibility and decreasing costs of next generation sequencing (NGS), availability of commercial kits, and development of dedicated analysis pipelines, have allowed an increasing number of laboratories to adopt this technology for HIV drug resistance (HIVDR) genotyping. Conventional HIVDR genotyping is traditionally carried out using population-based Sanger sequencing, which has a limited capacity for reliable detection of variants present at intra-host frequencies below a threshold of approximately 20%. NGS has the potential to improve sensitivity and quantitatively identify low-abundance variants, improving efficiency and lowering costs. However, some challenges exist for the standardization and quality assurance of NGS-based HIVDR genotyping. In this paper, we highlight considerations of these challenges as related to laboratory, clinical, and implementation of NGS for HIV drug resistance testing. Several sources of variation and bias occur in each step of the general NGS workflow, i.e., starting material, sample type, PCR amplification, library preparation method, instrument and sequencing chemistry-inherent errors, and data analysis options and limitations. Additionally, adoption of NGS-based HIVDR genotyping, especially for clinical care, poses pressing challenges, especially for resource-poor settings, including infrastructure and equipment requirements and cost, logistic and supply chains, instrument service availability, personnel training, validated laboratory protocols, and standardized analysis outputs. The establishment of external quality assessment programs may help to address some of these challenges and is needed to proceed with NGS-based HIVDR genotyping adoption.

Impact of the mutational load on the virological response to a first-line rilpivirine-based regimen

OBJECTIVES

To determine how the load of rilpivirine-resistant variants (mutational load) influences the virological response (VR) of HIV-1-infected patients to a rilpivirine-based first-line regimen.

PATIENTS AND METHODS

Four hundred and eighty-nine patients infected with HIV-1 whose reverse transcriptase gene had been successfully resistance genotyped using next-generation sequencing were given a first-line regimen containing rilpivirine. Variables associated with the VR at 12 months were identified using a logistic model. The results were used to build a multivariate model for each mutational load threshold and the R2 variations were analysed to identify the mutational load threshold that best predicted the VR.

RESULTS

The mutational load at baseline was the only variable linked to the VR at 12 months (P  < 0.01). The VR at 12 months decreased from 96.9% to 83.4% when the mutational load was >1700 copies/mL and to 50% when the mutational load was > 9000 copies/mL. The threshold of 9000 copies/mL was associated with the VR at 12 months with an OR of 36.7 (95% CI 4.7-285.1). The threshold of 1700 copies/mL was associated with the VR at 12 months with an OR of 7.2 (95% CI 1.4-36.8).

CONCLUSIONS

There is quantifiable evidence that determining a mutational load threshold can be used to identify those patients on a first-line regimen containing rilpivirine who are at risk of virological failure. The clinical management of HIV-infected patients can be improved by evaluating the frequency of mutant variants at a threshold of < 20% together with the plasma HIV-1 viral load at the time of resistance genotyping.

Are We Ready for NGS HIV Drug Resistance Testing? The Second "Winnipeg Consensus" Symposium

HIV drug resistance is a major global challenge to successful and sustainable antiretroviral therapy. Next-generation sequencing (NGS)-based HIV drug resistance (HIVDR) assays enable more sensitive and quantitative detection of drug-resistance-associated mutations (DRMs) and outperform Sanger sequencing approaches in detecting lower abundance resistance mutations. While NGS is likely to become the new standard for routine HIVDR testing, many technical and knowledge gaps remain to be resolved before its generalized adoption in regular clinical care, public health, and research. Recognizing this, we conceived and launched an international symposium series on NGS HIVDR, to bring together leading experts in the field to address these issues through in-depth discussions and brainstorming. Following the first symposium in 2018 (Winnipeg, MB Canada, 21-22 February, 2018), a second "Winnipeg Consensus" symposium was held in September 2019 in Winnipeg, Canada, and was focused on external quality assurance strategies for NGS HIVDR assays. In this paper, we summarize this second symposium's goals and highlights.

External Quality Assessment for Next-Generation Sequencing-Based HIV Drug Resistance Testing: Unique Requirements and Challenges

Over the past decade, there has been an increase in the adoption of next generation sequencing (NGS) technologies for HIV drug resistance (HIVDR) testing. NGS far outweighs conventional Sanger sequencing as it has much higher throughput, lower cost when samples are batched and, most importantly, significantly higher sensitivities for variants present at low frequencies, which may have significant clinical implications. Despite the advantages of NGS, Sanger sequencing remains the gold standard for HIVDR testing, largely due to the lack of standardization of NGS-based HIVDR testing. One important aspect of standardization includes external quality assessment (EQA) strategies and programs. Current EQA for Sanger-based HIVDR testing includes proficiency testing where samples are sent to labs and the performance of the lab conducting such assays is evaluated. The current methods for Sanger-based EQA may not apply to NGS-based tests because of the fundamental differences in their technologies and outputs. Sanger-based genotyping reports drug resistance mutations (DRMs) data as dichotomous, whereas NGS-based HIVDR genotyping also reports DRMs as numerical data (percent abundance). Here we present an overview of the need to develop EQA for NGS-based HIVDR testing and some unique challenges that may be encountered.

Evaluation of Vela Diagnostics HIV-1 genotyping assay on an automated next generation sequencing platform

BACKGROUND

Antiretroviral drug resistance testing is an integral part of the management of patients infected with HIV. The traditional Sanger sequencing method is capable of detecting drug resistant mutations (DRMs) that make up at least 10-15% of the viral quasispecies population. Newer next generation sequencing technologies have a greater sensitivity for the detection of minority variant DRMs down to around 1% of the population.

OBJECTIVES

Here NGS sequencing on the Vela Diagnostics automated next generation sequencing platform was evaluated and compared to the currently used Sanger sequencing method.

STUDY DESIGN

Sequences from both methods were obtained from a total of 79 patients, with a range of subtypes (CRF01_AE, A1/G, A1/CRF01_AE, A1/CRF02_AG, A1, A, B, C, CRF01_AG, CRF 06_CPX, D, G, B/G, CRF 57_BC/C, G/CRF 02_AG and CRF 14_BG/G) and viral loads (2.43-7 log10 copies/ml).

RESULTS

A high concordance was seen between the two methods for subtyping (96%) and majority variant detection (97.9%). NGS sequencing detected more variants and DRMs than Sanger sequencing. Of the 76 patient samples 86% (n = 66) had identical drug resistance reports. From the ten discrepant reports, nine had extra DRMs detected by NGS sequencing and all discrepancies were seen for NRTI and NNRTI antiviral resistance.

CONCLUSIONS

This study demonstrated a good performance of the NGS method for HIV-1 genotyping compared to the Sanger sequencing method for detection of majority variants, however the reproducibility for the detection of minority variants was sub-optimal. Adoption of an NGS sequencing approach has the potential to improve the clinical management of HIV-infected patients.

Performance comparison of next generation sequencing analysis pipelines for HIV-1 drug resistance testing

Next generation sequencing (NGS) is a trending new standard for genotypic HIV-1 drug resistance (HIVDR) testing. Many NGS HIVDR data analysis pipelines have been independently developed, each with variable outputs and data management protocols. Standardization of such analytical methods and comparison of available pipelines are lacking, yet may impact subsequent HIVDR interpretation and other downstream applications. Here we compared the performance of five NGS HIVDR pipelines using proficiency panel samples from NIAID Virology Quality Assurance (VQA) program. Ten VQA panel specimens were genotyped by each of six international laboratories using their own in-house NGS assays. Raw NGS data were then processed using each of the five different pipelines including HyDRA, MiCall, PASeq, Hivmmer and DEEPGEN. All pipelines detected amino acid variants (AAVs) at full range of frequencies (1~100%) and demonstrated good linearity as compared to the reference frequency values. While the sensitivity in detecting low abundance AAVs, with frequencies between 1~20%, is less a concern for all pipelines, their specificity dramatically decreased at AAV frequencies <2%, suggesting that 2% threshold may be a more reliable reporting threshold for ensured specificity in AAV calling and reporting. More variations were observed among the pipelines when low abundance AAVs are concerned, likely due to differences in their NGS read quality control strategies. Findings from this study highlight the need for standardized strategies for NGS HIVDR data analysis, especially for the detection of minority HIVDR variants.

OLA-Simple: A software-guided HIV-1 drug resistance test for low-resource laboratories

Background

HIV drug resistance (HIVDR) testing can assist clinicians in selecting treatments. However, high complexity and cost of genotyping assays limit routine testing in settings where HIVDR prevalence has reached high levels.

Methods

The oligonucleotide ligation assay (OLA)-Simple kit was developed for detection of HIVDR against first-line non-nucleoside/nucleoside reverse transcriptase inhibitors and validated on 672 codons (168 specimens) from subtypes A, B, C, D, and AE. The kit uses dry reagents to facilitate assay setup, lateral flow devices for visual HIVDR detections, and in-house software with an interface for guiding users and analyzing results.

Findings

HIVDR analysis of specimens by OLA-Simple compared to Sanger sequencing revealed 99.6 ± 0.3% specificity and 98.2 ± 0.9% sensitivity, and compared to high-sensitivity assays, 99.6 ± 0.6% specificity and 86.2 ± 2.5% sensitivity, with 2.6 ± 0.9% indeterminate results. OLA-Simple was performed more rapidly compared to Sanger sequencing (<4 h vs. 35–72 h). Forty-one untrained volunteers blindly tested two specimens each with 96.8 ± 0.8% accuracy.

Interpretation

OLA-Simple compares favorably with HIVDR genotyping by Sanger and sensitive comparators. Instructional software enabled inexperienced, first-time users to perform the assay with high accuracy. The reduced complexity, cost, and training requirements of OLA-Simple could improve access to HIVDR testing in low-resource settings and potentially allow same-day selection of appropriate antiretroviral therapy.

Fund

USA National Institutes of Health R01; the Clinical and Retrovirology Research Core and the Molecular Profiling and Computational Biology Core of the UW CFAR; Seattle Children's Research Institute; UW Holloman Innovation Challenge Award; Pilcher Faculty Fellowship.

Prospective Evaluation of the Vela Diagnostics Next-Generation Sequencing Platform for HIV-1 Genotypic Resistance Testing

Genotypic antiretroviral drug resistance testing is a critical component of the global efforts to control the HIV-1 epidemic. This study investigates the semiautomated, next-generation sequencing (NGS)-based Vela Diagnostics Sentosa SQ HIV-1 Genotyping Assay in a prospective cohort of HIV-1-infected patients. Two-hundred sixty-nine samples were successfully sequenced by both NGS and Sanger sequencing. Among the 261 protease/reverse transcriptase (PR/RT) sequences, a mean of 0.37 drug resistance mutations were identified by both Sanger and NGS, 0.08 by NGS alone, and 0.03 by Sanger alone. Among the 50 integrase sequences, a mean of 0.3 drug resistance mutations were detected by both Sanger and NGS, and 0.08 by NGS alone. NGS estimated higher levels of drug resistance to one or more antiretroviral drugs for 6.5% of PR/RT sequences and 4.0% of integrase sequences, whereas Sanger estimated higher levels of drug resistance for 3.8% of PR/RT sequences. Although the samples successfully sequenced by the Sentosa SQ HIV Genotyping Assay demonstrated similar predicted resistance compared with Sanger, 44% of Sentosa runs failed quality control requiring 17 additional runs. This semi-automated NGS-based assay may aid in HIV-1 genotypic drug resistance testing, though numerous quality control issues were observed when this platform was used in a clinical laboratory setting. With additional refinement, the Sentosa SQ HIV-1 Genotyping Assay may contribute to the global efforts to control HIV-1.

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.

Dynamic interactions between RNA viruses and human hosts unravelled by a decade of next generation sequencing

BACKGROUND

Next generation sequencing (NGS) methods have significantly contributed to a paradigm shift in genomic research for nearly a decade now. These methods have been useful in studying the dynamic interactions between RNA viruses and human hosts.

SCOPE OF THE REVIEW

In this review, we summarise and discuss key applications of NGS in studying the host - pathogen interactions in RNA viral infections of humans with examples.

MAJOR CONCLUSIONS

Use of NGS to study globally relevant RNA viral infections have revolutionized our understanding of the within host and between host evolution of these viruses. These methods have also been useful in clinical decision-making and in guiding biomedical research on vaccine design.

GENERAL SIGNIFICANCE

NGS has been instrumental in viral genomic studies in resolving within-host viral genomic variants and the distribution of nucleotide polymorphisms along the full-length of viral genomes in a high throughput, cost effective manner. In the future, novel advances such as long read, single molecule sequencing of viral genomes and simultaneous sequencing of host and pathogens may become the standard of practice in research and clinical settings. This will also bring on new challenges in big data analysis.

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.

Surveillance and Diagnosis of West Nile Virus in the Face of Flavivirus Cross-Reactivity

West Nile Virus (WNV) is an arthropod-borne flavivirus whose zoonotic cycle includes both mosquitoes and birds as amplifiers and humans and horses as dead-end hosts. In recent years WNV has been spreading globally and is currently endemic in Africa, The Middle East, India, Australia, central and southern Europe, and the Americas. Integrated surveillance schemes and environmental data aim to detect viral circulation and reduce the risk of infection for the human population emphasizing the critical role for One Health principles in public health. Approximately 20% of WNV infected patients develop West Nile Fever while in less than 1%, infection results in West Nile Neurological Disease. Currently, the diagnosis of WNV infection is primarily based on serology, since molecular identification of WNV RNA is unreliable due to the short viremia. The recent emergence of Zika virus epidemic in America and Asia has added another layer of complexity to WNV diagnosis due to significant cross-reactivity between several members of the Flaviviridae family such as Zika, dengue, Usutu, and West Nile viruses. Diagnosis is especially challenging in persons living in regions with flavivirus co-circulation as well as in travelers from WNV endemic countries traveling to Zika or dengue infected areas or vise-versa. Here, we review the recent studies implementing WNV surveillance of mosquitoes and birds within the One Health initiative. Furthermore, we discuss the utility of novel molecular methods, alongside traditional molecular and serological methods, in WNV diagnosis and epidemiological research.

Bioinformatic data processing pipelines in support of next-generation sequencing-based HIV drug resistance testing: the Winnipeg Consensus

INTRODUCTION

Next-generation sequencing (NGS) has several advantages over conventional Sanger sequencing for HIV drug resistance (HIVDR) genotyping, including detection and quantitation of low-abundance variants bearing drug resistance mutations (DRMs). However, the high HIV genomic diversity, unprecedented large volume of data, complexity of analysis and potential for error pose significant challenges for data processing. Several NGS analysis pipelines have been developed and used in HIVDR research; however, the absence of uniformity in data processing strategies results in lack of consistency and comparability of outputs from different pipelines. To fill this gap, an international symposium on bioinformatic strategies for NGS-based HIVDR testing was held in February 2018 in Winnipeg, Canada, convening laboratory scientists, bioinformaticians and clinicians involved in four recently developed, publicly available NGS HIVDR pipelines. The goal of this symposium was to establish a consensus on effective bioinformatic strategies for NGS data management and its use for HIVDR reporting.

DISCUSSION

Essential functionalities of an NGS HIVDR pipeline were divided into five analytic blocks: (1) NGS read quality control (QC)/quality assurance (QA); (2) NGS read alignment and reference mapping; (3) HIV variant calling and variant QC; (4) NGS HIVDR reporting; and (5) extended data applications and additional considerations for data management. The consensuses reached among the participants on all major aspects of these blocks are summarized here. They encompass not only recommended data management and analysis strategies, but also detailed bioinformatic approaches that help ensure accuracy of the derived HIVDR analysis outputs for both research and potential clinical use.

CONCLUSIONS

While NGS is being adopted more broadly in HIVDR testing laboratories, data processing is often a bottleneck hindering its generalized application. The proposed standardization of NGS read QC/QA, read alignment and reference mapping, variant calling and QC, HIVDR reporting and relevant data management strategies in this "Winnipeg Consensus" may serve as a starting guideline for NGS HIVDR data processing that informs the refinement of existing pipelines and those yet to be developed. Moreover, the bioinformatic strategies presented here may apply more broadly to NGS data analysis of microbes harbouring significant genomic diversity.

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.

Sanger and Next Generation Sequencing Approaches to Evaluate HIV-1 Virus in Blood Compartments

The implementation of antiretroviral treatment combined with the monitoring of drug resistance mutations improves the quality of life of HIV-1 positive patients. The drug resistance mutation patterns and viral genotypes are currently analyzed by DNA sequencing of the virus in the plasma of patients. However, the virus compartmentalizes, and different T cell subsets may harbor distinct viral subsets. In this study, we compared the patterns of HIV distribution in cell-free (blood plasma) and cell-associated viruses (peripheral blood mononuclear cells, PBMCs) derived from ART-treated patients by using Sanger sequencing- and Next-Generation sequencing-based HIV assay. CD4⁺CD45RA⁻RO⁺ memory T-cells were isolated from PBMCs using a BD FACSAria instrument. HIV pol (protease and reverse transcriptase) was RT-PCR or PCR amplified from the plasma and the T-cell subset, respectively. Sequences were obtained using Sanger sequencing and Next-Generation Sequencing (NGS). Sanger sequences were aligned and edited using RECall software (beta v3.03). The Stanford HIV database was used to evaluate drug resistance mutations. Illumina MiSeq platform and HyDRA Web were used to generate and analyze NGS data, respectively. Our results show a high correlation between Sanger sequencing and NGS results. However, some major and minor drug resistance mutations were only observed by NGS, albeit at different frequencies. Analysis of low-frequency drugs resistance mutations and virus distribution in the blood compartments may provide information to allow a more sustainable response to therapy and better disease management.

HIV drug resistance patterns in pregnant women using next generation sequence in Mozambique

Background

Few data on HIV resistance in pregnancy are available from Mozambique, one of the countries with the highest HIV toll worldwide. Understanding the patterns of HIV drug resistance in pregnant women might help in tailoring optimal regimens for prevention of mother to child transmission of HIV (pMTCT) and antenatal care.

Objectives

To describe the frequency and characteristics of HIV drug resistance mutations (HIVDRM) in pregnant women with virological failure at delivery, despite pMTCT or antiretroviral therapy (ART).

Methods

Samples from HIV-infected pregnant women from a rural area in southern Mozambique were analysed. Only women with HIV-1 RNA >400c/mL at delivery were included in the analysis. HIVDRM were determined using MiSeq® (detection threshold 1%) at the first antenatal care (ANC) visit and at the time of delivery.

Results

Ninety and 60 samples were available at the first ANC visit and delivery, respectively. At first ANC, 97% of the women had HIV-1 RNA>400c/mL, 39% had CD4+ counts <350 c/mm³ and 30% were previously not on ART. Thirteen women (14%) had at least one HIVDRM of whom 70% were not on previous ART. Eight women (13%) had at least one HIVDRM at delivery. Out of 37 women with data available from the two time points, 8 (21%) developed at least one new HIVDRM during pMTCT or ART. Twenty seven per cent (53/191), 32% (44/138) and 100% (5/5) of the mutations that were present at enrolment, delivery and that emerged during pregnancy, respectively, were minority mutations (frequency <20%).

Conclusions

Even with ultrasensitive HIV-1 genotyping, less than 20% of women with detectable viremia at delivery had HIVDRM before initiating pMTCT or ART. This suggests that factors other than pre-existing resistance, such as lack of adherence or interruptions of the ANC chain, are also relevant to explain lack of virological suppression at the time of delivery in women receiving antiretrovirals drugs during pregnancy.

Whole-exome sequencing improves genetic testing accuracy in pulmonary artery hypertension

Sanger sequencing, the traditional “gold standard” for mutation detection, has been wildly used in genetic testing of pulmonary artery hypertension (PAH). However, with the advent of whole-exome sequencing (WES), few studies have compared the accuracy of WES and Sanger sequencing in routine genetic testing of PAH. PAH individuals were enrolled from Fu Wai Hospital and Shanghai Pulmonary Hospital. WES was used to analyze DNA samples from 120 PAH patients whose bone morphogenetic protein receptor type 2 (BMPR2) mutation statuses had been previously studied using Sanger sequencing. The Sanger sequencing and WES agreement was 98.3% (118/120) with near-perfect agreement (κ coefficient = 0.848). There was no significant difference between the two methods on the McNemar–Bowker test (P > 0.05). Twenty-one BMPR2 mutation carriers and 99 non-carriers were detected by Sanger sequencing. Among the 21 BMPR2 carriers detected by Sanger sequencing, one variant (c.1040 T > A) was not found by WES. Among the 99 BMPR2 non-carriers, WES detected an extra mutation carrier (c.76 + 1 G > C) missed by Sanger sequencing. Both false-positive and false-negative results were highly conserved and were re-analyzed by Sanger sequencing. WES improved the accuracy of Sanger sequencing and detected 1% (1/99) false-positive and 4.8% (1/21) false-negative results of Sanger sequencing. No false-positive and false-negative results of WES were identified in our analysis. WES is non-inferior to Sanger sequencing and may play a more important role in genetic testing of PAH patients in the future.

A Bioinformatic Pipeline for Monitoring of the Mutational Stability of Viral Drug Targets with Deep-Sequencing Technology

The efficient development of antiviral drugs, including efficient antiviral small interfering RNAs (siRNAs), requires continuous monitoring of the strict correspondence between a drug and the related highly variable viral DNA/RNA target(s). Deep sequencing is able to provide an assessment of both the general target conservation and the frequency of particular mutations in the different target sites. The aim of this study was to develop a reliable bioinformatic pipeline for the analysis of millions of short, deep sequencing reads corresponding to selected highly variable viral sequences that are drug target(s). The suggested bioinformatic pipeline combines the available programs and the ad hoc scripts based on an original algorithm of the search for the conserved targets in the deep sequencing data. We also present the statistical criteria for the threshold of reliable mutation detection and for the assessment of variations between corresponding data sets. These criteria are robust against the possible sequencing errors in the reads. As an example, the bioinformatic pipeline is applied to the study of the conservation of RNA interference (RNAi) targets in human immunodeficiency virus 1 (HIV-1) subtype A. The developed pipeline is freely available to download at the website http://virmut.eimb.ru/. Brief comments and comparisons between VirMut and other pipelines are also presented.

Cryptic Microheteroresistance Explains Mycobacterium tuberculosis Phenotypic Resistance

RATIONALE

Minority drug-resistant Mycobacterium tuberculosis subpopulations can be associated with phenotypic resistance but are poorly detected by Sanger sequencing or commercial molecular diagnostic assays.

OBJECTIVES

To determine the role of targeted next-generation sequencing in resolving these minor variant subpopulations.

METHODS

We used single molecule overlapping reads (SMOR), a targeted next-generation sequencing approach that dramatically reduces sequencing error, to analyze primary cultured isolates phenotypically resistant to rifampin, fluoroquinolones, or aminoglycosides, but for which Sanger sequencing found no resistance-associated variants (RAVs) within respective resistance-determining regions (study group). Isolates also underwent single-colony selection on antibiotic-containing agar, blinded to sequencing results. As a positive control, isolates with multiple colocalizing chromatogram peaks were also analyzed (control group).

MEASUREMENTS AND MAIN RESULTS

Among 61 primary culture isolates (25 study group and 36 control group), SMOR described 66 (49%) and 45 (33%) of 135 total heteroresistant RAVs at frequencies less than 5% and less than 1% of the total mycobacterial population, respectively. In the study group, SMOR detected minor resistant variant subpopulations in 80% (n = 20/25) of isolates with no Sanger-identified RAVs (median subpopulation size, 1.0%; interquartile range, 0.2-3.9%). Single-colony selection on drug-containing media corroborated SMOR results for 90% (n = 18/20) of RAV-containing specimens, and the absence of RAVs in 60% (n = 3/5) of isolates. Overall, Sanger sequencing was concordant with SMOR for 77% (n = 53/69) of macroheteroresistant (5-95% total population), but only 5% of microheteroresistant (<5%) subpopulations (n = 3/66) across both groups.

CONCLUSIONS

Cryptic minor variant mycobacterial subpopulations exist below the resolving capability of current drug susceptibility testing methodologies, and may explain an important proportion of false-negative resistance determinations.