HIV drug resistance testing by high-multiplex "wide" sequencing on the MiSeq instrument

Optimization of genetic distance threshold for inferring the CRF01_AE molecular network based on next-generation sequencing

Introduction

HIV molecular network based on genetic distance (GD) has been extensively utilized. However, the GD threshold for the non-B subtype differs from that of subtype B. This study aimed to optimize the GD threshold for inferring the CRF01_AE molecular network.

Methods

Next-generation sequencing data of partial CRF01_AE pol sequences were obtained for 59 samples from 12 transmission pairs enrolled from a high-risk cohort during 2009 and 2014. The paired GD was calculated using the Tamura-Nei 93 model to infer a GD threshold range for HIV molecular networks.

Results

2,019 CRF01_AE pol sequences and information on recent HIV infection (RHI) from newly diagnosed individuals in Shenyang from 2016 to 2019 were collected to construct molecular networks to assess the ability of the inferred GD thresholds to predict recent transmission events. When HIV transmission occurs within a span of 1-4 years, the mean paired GD between the sequences of the donor and recipient within the same transmission pair were as follow: 0.008, 0.011, 0.013, and 0.023 substitutions/site. Using these four GD thresholds, it was found that 98.9%, 96.0%, 88.2%, and 40.4% of all randomly paired GD values from 12 transmission pairs were correctly identified as originating from the same transmission pairs. In the real world, as the GD threshold increased from 0.001 to 0.02 substitutions/site, the proportion of RHI within the molecular network gradually increased from 16.6% to 92.3%. Meanwhile, the proportion of links with RHI gradually decreased from 87.0% to 48.2%. The two curves intersected at a GD of 0.008 substitutions/site.

Discussion

A suitable range of GD thresholds, 0.008-0.013 substitutions/site, was identified to infer the CRF01_AE molecular transmission network and identify HIV transmission events that occurred within the past three years. This finding provides valuable data for selecting an appropriate GD thresholds in constructing molecular networks for non-B subtypes.

Comparison of HIV-1 DNA load measurements in blood and in relation to successful proviral sequencing

OBJECTIVE

HIV DNA sequencing is now routinely used for HIV-infected individuals on antiretroviral therapy (ART) with or without partial genotypic history. Successful amplification of HIV pol gene has yet to be correlated with HIV DNA levels. Here, we assessed the relationship between HIV DNA load and sequencing results.

METHODS

We analyzed three different qPCR measurements of total (LTR and LTR-gag) and integrated (Alu-LTR) HIV DNA in blood samples collected from viremic as well as virally suppressed HIV-infected individuals on ART. HIV DNA levels were compared to HIV DNA Sanger sequencing and clinical and therapeutic parameters.

RESULTS

Among the 135 individuals analyzed for HIV DNA measurements and sequencing, all three HIV DNA measurements were associated with HIV DNA Sanger sequencing results. A threshold of around 2 and 1.5 log copies/million leukocytes of total HIV DNA was identified for LTR and LTR-gag qPCRs, respectively. Integrated HIV DNA positivity was also associated with successful sequencing. We further compared HIV DNA measurement techniques in an extended cohort of 312 individuals and showed that all measurements correlated between the different techniques, regardless of the HIV-1 subtypes analyzed. However, higher detection rates were observed with LTR (96%) compared to LTR-gag (86%) and Alu-LTR (59%) qPCRs. Duration of virological control on ART and CD4 nadir were the main determinants of HIV reservoir size.

CONCLUSIONS

HIV DNA measurement is associated with Sanger sequencing success, regardless of the technique used. In a clinical setting, Application of HIV DNA quantification before sequencing should be further evaluated.

Optimizing Antiretroviral Therapy in Heavily ART-Experienced Patients with Multi-Class Resistant HIV-1 Using Proviral DNA Genotypic Resistance Testing

Resistance to multiple antiretroviral drugs among people living with HIV (PLWH) can result in a high pill burden, causing toxicity and drug interactions. Thus, the goal is to simplify treatment regimens while maintaining effectiveness. However, former resistance analysis data may not be current or complete. The use of proviral DNA genotyping may assist in selecting appropriate treatment options. A retrospective study was carried out on individuals belonging to the Cologne HIV cohort with a resistance history to two or more antiretroviral (ARV) classes and on non-standard antiretroviral therapy (ART). Patients required former viral RNA and a recent proviral DNA resistance test to be available prior to the switch to ART. Potential discrepancies between resistance test results obtained through RNA and proviral DNA methods and the consequent virological and clinical outcomes following ART adjustments were analyzed. Out of 1250 patients, 35 were eligible for inclusion in this study. The median length of known HIV infection was 27 years, and the median duration of ART was 22 years. Of the 35 participants, 16 had received all five ARV classes. Based on proviral DNA genotyping results, ART was simplified in 17 patients. At the last follow-up examination after changing therapy, 15 patients had HIV RNA <50 copies/mL (median 202 days, range 21-636). The mean number of pills per day decreased from eight to three, and the median intake frequency decreased from two to one time/day (ranges 1-2). Our study supports the use of proviral DNA genotyping as a safe strategy for switching to simplified ART regimens. However, the lack of extensive research on the advantages of proviral DNA genotyping makes it challenging to fully assess its benefits in terms of treatment selection.

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.

HIV-1 resistance genotyping by ultra-deep sequencing and 6-month virological response to first-line treatment

OBJECTIVES

To evaluate the routine use of the Sentosa ultra-deep sequencing (UDS) system for HIV-1 polymerase resistance genotyping in treatment-naïve individuals and to analyse the virological response (VR) to first-line antiretroviral treatment.

METHODS

HIV drug resistance was determined on 237 consecutive samples from treatment-naïve individuals using the Sentosa UDS platform with two mutation detection thresholds (3% and 20%). VR was defined as a plasma HIV-1 virus load <50 copies/mL after 6 months of treatment.

RESULTS

Resistance to at least one antiretroviral drug with a mutation threshold of 3% was identified in 29% and 16% of samples according to ANRS and Stanford algorithms, respectively. The ANRS algorithm also revealed reduced susceptibility to at least one protease inhibitor (PI) in 14.3% of samples, to one reverse transcriptase inhibitor in 12.7%, and to one integrase inhibitor (INSTI) in 5.1%. For a mutation threshold of 20%, resistance was identified in 24% and 13% of samples according to ANRS and Stanford algorithms, respectively. The 6 months VR was 87% and was similar in the 58% of patients given INSTI-based treatment, in the 16% given PI-based treatment and in the 9% given NNRTI-based treatment. Multivariate analysis indicated that the VR was correlated with the baseline HIV virus load and resistance to at least one PI at both 3% and 20% mutation detection thresholds (ANRS algorithm).

CONCLUSIONS

The Vela UDS platform is appropriate for determining antiretroviral resistance in patients on a first-line antiretroviral treatment. Further studies are needed on the use of UDS for therapeutic management.

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.

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.

Analytical Assessment of the Vela Diagnostics NGS Assay for HIV Genotyping and Resistance Testing: The Apulian Experience

Drug-resistance monitoring is one of the hardest challenges in HIV management. Next-generation sequencing (NGS) technologies speed up the detection of drug resistance, allowing the adjustment of antiretroviral therapy and enhancing the quality of life of people living with HIV. Recently, the NGS Sentosa^® SQ HIV Genotyping Assay (Vela Diagnostics) received approval for in vitro diagnostics use. This work is the first Italian evaluation of the performance of the Vela Diagnostics NGS platform, assessed with 420 HIV-1 clinical samples. A comparison with Sanger sequencing performance is also reported, highlighting the advantages and disadvantages of the Sentosa^® NGS assay. The precision of the technology was studied with reference specimens, while intra- and inter-assay reproducibility were evaluated for selected clinical samples. Vela Diagnostics’ NGS assay reached an 87% success rate through 30 runs of analysis in a real-world clinical context. The concordance with Sanger sequencing outcomes was equal to 97.2%. Several detected mismatches were due to NGS’s superior sensitivity to low-frequency variants. A high accuracy was observed in testing reference samples. Repeatability and reproducibility assays highlighted the good performance of the NGS platform. Beyond a few technical issues that call for further optimization, the key improvement will be a better balance between costs and processing speed. Once these issues have been solved, the Sentosa^® SQ HIV Genotyping Assay will be the way forward for HIV resistance testing.

Genetic Diversity and Acquired Drug Resistance Mutations Detected by Deep Sequencing in Virologic Failures among Antiretroviral Treatment Experienced Human Immunodeficiency Virus-1 Patients in a Pastoralist Region of Ethiopia

Purpose

This study was conducted to investigate the drug resistance mutations and genetic diversity of HIV-1 in ART experienced patients in South Omo, Ethiopia.

Patients and Methods

A cross-sectional study conducted on 253 adult patients attending ART clinics for ≥6 months in South Omo. Samples with VL ≥1000 copies/mL were considered as virological failures (VF) and their reverse transcriptase gene codons 90–234 were sequenced using Illumina MiSeq. MinVar was used for the identification of the subtypes and drug resistance mutations. Phylogenetic tree was constructed by neighbor-joining method using the maximum likelihood model.

Results

The median duration of ART was 51 months and 18.6% (47/253) of the patients exhibited VF. Of 47 viraemic patients, the genome of 41 were sequenced and subtype C was dominant (87.8%) followed by recombinant subtype BC (4.9%), M-09-CPX (4.9) and BF1 (2.4%). Of 41 genotyped subjects, 85.4% (35/41) had at least one ADR mutation. Eighty-one percent (33/41) of viraemic patients harbored NRTI resistance mutations, and 48.8% (20/41) were positive for NNRTI resistance mutations, with 43.9% dual resistance mutations. Among NRTI resistance mutations, M184V (73.2%), K219Q (63.4%) and T215 (56.1%) complex were the most mutated positions, while the most common NNRTI resistance mutations were K103N (24.4%), K101E, P225H and V108I 7.5% each. Active tuberculosis (aOR=13, 95% CI= 3.46–29.69), immunological failure (aOR=3.61, 95% CI=1.26–10.39), opportunistic infections (aOR=8.39, 95% CI= 1.75–40.19), and poor adherence were significantly associated with virological failure, while rural residence (aOR 2.37; 95% CI: 1.62–9.10, P= 0.05), immunological failures (aOR 2.37; 95% CI: 1.62–9.10, P= 0.05) and high viral load (aOR 16; 95% CI: 5.35 51.59, P <0.001) were predictors of ADR mutation among the ART experienced and viraemic study subjects.

Conclusion

The study revealed considerable prevalence of VF and ADR mutation with the associated risk indicators. Regular virological monitoring and drug resistance genotyping methods should be implemented for better ART treatment outcomes of the nation.

Kinetics of archived M184V mutation in treatment-experienced virally suppressed HIV-infected patients

BACKGROUND

We aimed to assess the kinetics of drug-resistant viral variants (DRVs) harboring the M184V mutation in the proviral DNA of long-term virally suppressed patients, and factors associated with DRV persistence.

METHODS

HIV-DNA from blood cells stored in 2019 and 2016 was sequenced using both Sanger and ultradeep sequencing (SS and UDS, with a detection threshold of 1%) in ART-treated patients with HIV-RNA <50 copies/mL for at least 5 years, with past M184V mutation documented in HIV-RNA.

RESULTS

Among the 79 tested patients, by combining SS and UDS, the M184V was found to be absent in 26/79 (33%) patients (M184V- patients), and persisted in 53/79 (67%) (M184V+ patients). The M184V+ patients had a longer history of ART, a lower CD4 nadir and higher pretherapeutic HIV-RNA. Among the 37 patients with viral sequences assessed by UDS, the proportion of M184V+ DRVs significantly decreased between 2016 and 2019 (40% versus 14%, p=0.005). The persistence of M184V was associated with the duration and level of HIV-RNA replication under 3TC/FTC (p=0.0009 and p=0.009, respectively).

CONCLUSION

While it decreased over time in HIV-DNA, the M184V mutation was more frequently persistent in the HIV-DNA of more experienced patients with longer past replication under 3TC/FTC.

Laboratory Optimization Tweaks for Sanger Sequencing in a Resource-Limited Setting

Despite various challenges that hinder the implementation of high-tech molecular methods in resource-limited settings, we have been able to implement and achieve International Organization for Standardization 15189:2012 accreditation for genotypic HIV drug resistance testing in our facility. At the Center for Human Virology and Genomics, Nigerian Institute of Medical Research, Nigeria has recorded a high sequencing success rate and good quality sequence data. This was achieved by optimizing laboratory processes from 2008 to the current date. We have optimized sample preparation, RT-PCR, several post-PCR processes, and the cycle sequencing to improve the sensitivity of amplification even with limited plasma samples and low viral copy numbers. The optimized workflow maximizes output, minimizes reagent wastage, and achieves substantial cost savings without compromising the quality of the sequence data. Our performance at our last external quality assurance program is a testimonial to the efficiency of the workflow. For the 5-sample panel, each with 67-68 mutation points evaluated, we scored 100% for all 5 specimens. Our optimized laboratory workflow is thus documented to support laboratories and to help researchers achieve excellent results the first time and eliminate contamination while minimizing the wastage of costly sequencing reagents.

From Sequence Data to Patient Result: A Solution for HIV Drug Resistance Genotyping With Exatype, End to End Software for Pol-HIV-1 Sanger Based Sequence Analysis and Patient HIV Drug Resistance Result Generation

Introduction:

With the rapid scale-up of antiretroviral therapy (ART) to treat HIV infection, there are ongoing concerns regarding probable emergence and transmission of HIV drug resistance (HIVDR) mutations. This scale-up has to lead to an increased need for routine HIVDR testing to inform the clinical decision on a regimen switch. Although the majority of wet laboratory processes are standardized, slow, labor-intensive data transfer and subjective manual sequence interpretation steps are still required to finalize and release patient results. We thus set out to validate the applicability of a software package to generate HIVDR patient results from raw sequence data independently.

Methods:

We assessed the performance characteristics of Hyrax Bioscience’s Exatype (a sequence data to patient result, fully automated sequence analysis software, which consolidates RECall, MEGA X and the Stanford HIV database) against the standard method (RECall and Stanford database). Exatype is a web-based HIV Drug resistance bioinformatic pipeline available at sanger.exatype.com. To validate the exatype, we used a test set of 135 remnant HIV viral load samples at the National HIV Reference Laboratory (NHRL).

Result:

We analyzed, and successfully generated results of 126 sequences out of 135 specimens by both Standard and Exatype software. Result production using Exatype required minimal hands-on time in comparison to the Standard (6 computation-hours using the standard method versus 1.5 Exatype computation-hours). Concordance between the 2 systems was 99.8% for 311,227 bases compared. 99.7% of the 0.2% discordant bases, were attributed to nucleotide mixtures as a result of the sequence editing in Recall. Both methods identified similar (99.1%) critical antiretroviral resistance-associated mutations resulting in a 99.2% concordance of resistance susceptibility interpretations. The Base-calling comparison between the 2 methods had Cohen’s kappa (0.97 to 0.99), implying an almost perfect agreement with minimal base calling variation. On a predefined dataset, RECall editing displayed the highest probability to score mixtures accurately 1 vs. 0.71 and the lowest chance to inaccurately assign mixtures to pure nucleotides (0.002–0.0008). This advantage is attributable to the manual sequence editing in RECall.

Conclusion:

The reduction in hands-on time needed is a benefit when using the Exatype HIV DR sequence analysis platform and result generation tool. There is a minimal difference in base calling between Exatype and standard methods. Although the discrepancy has minimal impact on drug resistance interpretation, allowance of sequence editing in Exatype as RECall can significantly improve its performance.

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

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

Genotypic Methods for HIV Drug Resistance Monitoring: The Opportunities and Challenges Faced by China

AIDS is a globalized infectious disease. In 2014, UNAIDS launched a global project of "90-90-90" to end the HIV epidemic by 2030. The second and third 90 require 90% of HIV-1 infected individuals receiving antiretroviral therapy (ART) and durable virological suppression. However, wide use of ART will greatly increase the emergence and spreading of HIV drug resistance and current HIV drug resistance test (DRT) assays in China are seriously lagging behind, hindering to achieve virological suppression. Therefore, recommending an appropriate HIV DRT method is critical for HIV routine surveillance and prevention in China. In this review, we summarized the current existing HIV drug resistance genotypic testing methods around the world and discussed the advantages and disadvantages of these methods.

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.

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 Program for Next-Generation Sequencing-Based HIV Drug Resistance Testing: Logistical Considerations

Next-generation sequencing (NGS) is likely to become the new standard method for HIV drug resistance (HIVDR) genotyping. Despite the significant advances in the development of wet-lab protocols and bioinformatic data processing pipelines, one often-missing critical component of an NGS HIVDR assay for clinical use is external quality assessment (EQA). EQA is essential for ensuring assay consistency and laboratory competency in performing routine biomedical assays, and the rollout of NGS HIVDR tests in clinical practice will require an EQA. In September 2019, the 2nd International Symposium on NGS HIVDR was held in Winnipeg, Canada. It convened a multidisciplinary panel of experts, including research scientists, clinicians, bioinformaticians, laboratory biologists, biostatisticians, and EQA experts. A themed discussion was conducted on EQA strategies towards such assays during the symposium. This article describes the logistical challenges identified and summarizes the opinions and recommendations derived from these discussions, which may inform the development of an inaugural EQA program for NGS HIVDR in the near future.

A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing

Washington, DC continues to experience a generalized HIV-1 epidemic. We characterized the local phylodynamics of HIV-1 in DC using next-generation sequencing (NGS) data. Viral samples from 68 participants from 2016 through 2017 were sequenced and paired with epidemiological data. Phylogenetic and network inferences, drug resistant mutations (DRMs), subtypes and HIV-1 diversity estimations were completed. Haplotypes were reconstructed to infer transmission clusters. Phylodynamic inferences based on the HIV-1 polymerase (pol) and envelope genes (env) were compared. Higher HIV-1 diversity (n.s.) was seen in men who have sex with men, heterosexual, and male participants in DC. 54.0% of the participants contained at least one DRM. The 40-49 year-olds showed the highest prevalence of DRMs (22.9%). Phylogenetic analysis of pol and env sequences grouped 31.9-33.8% of the participants into clusters. HIV-TRACE grouped 2.9-12.8% of participants when using consensus sequences and 9.0-64.2% when using haplotypes. NGS allowed us to characterize the local phylodynamics of HIV-1 in DC more broadly and accurately, given a better representation of its diversity and dynamics. Reconstructed haplotypes provided novel and deeper phylodynamic insights, which led to networks linking a higher number of participants. Our understanding of the HIV-1 epidemic was expanded with the powerful coupling of HIV-1 NGS data with epidemiological data.

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.

Performance evaluation of the Vela Dx Sentosa next-generation sequencing system for HIV-1 DNA genotypic resistance

BACKGROUND

Patients on antiretroviral therapy could benefit from HIV-1 DNA resistance genotyping for exploring virological failure with low viral load or to guide treatment simplification. Few new generation sequencing data are available.

OBJECTIVE

To check that the automated deep sequencing Sentosa platform (Vela DX) detected minority resistant variants well enough for HIV DNA genotyping.

STUDY DESIGN

We evaluated the Sentosa SQ HIV genotyping assay with automated extraction on 40 DNA longitudinal samples from treatment-experienced patients by comparison with Sanger sequencing. HIV drug resistance was interpreted using the ANRS algorithm (v29) at the threshold of 20 % and 3 %.

RESULTS

The Sentosa SQ HIV genotyping assay was 100 % successful to amplify and sequence PR and RT and 86 % to amplify and sequence IN when the HIV DNA load was >2.5 log copies/million cells. The Sentosa and Sanger sequencing were concordant for predicting PR-RT resistance at the threshold of 20 % in 14/18 samples successfully sequenced. A higher level of resistance was predicted by Sentosa in three samples and by Sanger in one sample. The prevalence of resistance was 7 % to PI, 59 % to NRTI, 31 % to NNRTI and 20 % to integrase inhibitors using the Sentosa SQ genotyping assay at the threshold of 3 %. Seven additional mutations <20 % were detected using the Sentosa assay.

CONCLUSION

Automated DNA extraction and sequencing using the Sentosa SQ HIV genotyping assay accurately predicted HIV DNA drug resistance by comparison with Sanger. Prospective studies are needed to evaluate the clinical interest of HIV DNA genotyping.

Phylogenetic surveillance of travel-related Zika virus infections through whole-genome sequencing methods

In 2018, the World Health Organization identified the Zika virus (ZIKV) as a pathogen that should be prioritized for public health research due to its epidemic potential. In this study, whole-genome sequencing (WGS) of travel-acquired ZIKV infections was used to examine the limitations of phylogenetic analysis. WGS and phylogenetic analysis were performed to investigate geographic clustering of samples from five Canadians with travel-acquired ZIKV infections and to assess the limitations of phylogenetic analysis of ZIKV sequences using a phylogenetic cluster approach. Genomic variability of ZIKV samples was assessed and for context, compared with hepatitis C virus (HCV) samples. Phylogenetic analysis confirmed the suspected region of ZIKV infection for one of five samples and one sample failed to cluster with sequences from its suspected country of infection. Travel-acquired ZIKV samples depicted low genomic variability relative to HCV samples. A floating patristic distance threshold classified all pre-2000 ZIKV sequences into separate clusters, while only Cambodian, Peruvian, Malaysian, and South Korean sequences were similarly classifiable. While phylogenetic analysis of ZIKV data can identify the broad geographical region of ZIKV infection, ZIKV’s low genomic variability is likely to limit precise interpretations of phylogenetic analysis of the origins of travel-related cases.

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.

Impact of Next-generation Sequencing Defined Human Immunodeficiency Virus Pretreatment Drug Resistance on Virological Outcomes in the ANRS 12249 Treatment-as-Prevention Trial

BACKGROUND

Previous studies in human immunodeficiency virus (HIV)-positive individuals on thymidine analogue backbone antiretroviral therapy (ART) with either nevirapine or efavirenz have suggested poorer virological outcomes in the presence of pretreatment drug resistance (PDR). We assessed the impact of PDR on virological suppression (VS; <50 copies/mL) in individuals prescribed primarily tenofovir/emtricitabine/efavirenz in rural KwaZulu-Natal within a treatment-as-prevention trial.

METHODS

Among 1557 HIV-positive individuals who reported no prior ART at study entry and provided plasma samples, 1328 individuals with entry viral load (VL) >1000 copies/mL had next-generation sequencing (NGS) of the HIV pol gene with MiSeq technology. Results were obtained for 1148 individuals, and the presence of PDR was assessed at 5% and 20% detection thresholds. Virological outcome was assessed using Cox regression in 837 of 920 ART initiators with at least 1 follow-up VL after ART initiation.

RESULTS

PDR prevalence was 9.5% (109/1148) and 12.8% (147/1148) at 20% and 5% thresholds, respectively. After a median of 1.36 years (interquartile range, 0.91-2.13), mostly on fixed-dose combination tenofovir/emtricitabine/efavirenz, presence of both nonnucleoside reverse transcriptase inhibitor (NNRTI)/nucleoside reverse transcriptase inhibitor PDR vs no PDR was associated with longer time to VS (adjusted hazard ratio [aHR], 0.32; 95% confidence interval [CI], 0.12-0.86), while there was no difference between those with only NNRTI PDR vs no PDR (aHR, 1.05; 95% CI, 0.82-1.34) at the 5% threshold. Similar differences were observed for mutations detected at the 20% threshold, although without statistical significance.

CONCLUSIONS

NGS uncovered a high prevalence of PDR among participants enrolled in trial clinics in rural KwaZulu-Natal. Dual-class PDR to a mainly tenofovir/emtricitabine/efavirenz regimen was associated with poorer VS. However, there was no impact of NNRTI PDR alone.

CLINICAL TRIALS TEGISTRATION

NCT01509508; South African National Clinical Trials Register: DOH-27-0512-3974.

HIV-1 genotypic resistance testing using the Vela automated next-generation sequencing platform

Objectives

To evaluate the diagnostic performance of the Vela next-generation sequencing (NGS) system in conjunction with the Sentosa SQ HIV Genotyping Assay for genotyping HIV-1.

Methods

Plasma RNA was extracted and templates prepared with the Sentosa SX instrument before sequencing the HIV-1 polymerase on the Sentosa SQ301 Sequencer (PGM IonTorrent). The Vela NGS System was compared with direct sequencing and the 454 GS-FLX (Roche) and MiSeq (Illumina) systems for genotypic resistance testing on clinical samples.

Results

The Vela NGS system detected majority resistance mutations in subtype B and CRF02-AG samples at 500 copies/mL and minority variants with a sensitivity of 5% at 100 000 copies/mL. The Vela NGS system and direct sequencing identified resistance mutations with 97% concordance in 46 clinical samples. Vela identified 1/20 of the 1%-5% mutations identified by 454, 5/12 of the 5%-20% mutations and 60/61 of the >20% mutations. Vela identified 3/14 of the 1%-5% mutations identified by MiSeq, 0/2 of the 5%-20% mutations and 47/47 of the >20% mutations. The resistance mutation quantifications by Vela and 454 were concordant (bias: 2.31%), as were those by Vela and MiSeq (bias: 1.06%).

Conclusions

The Vela NGS system provides automated nucleic acid extraction, PCR reagent distribution, library preparation and bioinformatics analysis. The analytical performance was very good when compared with direct sequencing, but was less sensitive than two other NGS platforms for detecting minority variants.

A MiSeq-HyDRA platform for enhanced HIV drug resistance genotyping and surveillance

Conventional HIV drug resistance (HIVDR) genotyping utilizes Sanger sequencing (SS) methods, which are limited by low data throughput and the inability of detecting low abundant drug resistant variants (LADRVs). Here we present a next generation sequencing (NGS)-based HIVDR typing platform that leverages the advantages of Illumina MiSeq and HyDRA Web. The platform consists of a fully validated sample processing protocol and HyDRA web, an open web portal that allows automated customizable NGS-based HIVDR data processing. This platform was characterized and validated using a panel of HIV-spiked plasma representing all major HIV-1 subtypes, pedigreed plasmids, HIVDR proficiency specimens and clinical specimens. All examined major HIV-1 subtypes were consistently amplified at viral loads of ≥1,000 copies/ml. The gross error rate of this platform was determined at 0.21%, and minor variations were reliably detected down to 0.50% in plasmid mixtures. All HIVDR mutations identifiable by SS were detected by the MiSeq-HyDRA protocol, while LADRVs at frequencies of 1~15% were detected by MiSeq-HyDRA only. As compared to SS approaches, the MiSeq-HyDRA platform has several notable advantages including reduced cost and labour, and increased sensitivity for LADRVs, making it suitable for routine HIVDR monitoring for both patient care and surveillance purposes.

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

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

Predicted antiviral activity of tenofovir versus abacavir in combination with a cytosine analogue and the integrase inhibitor dolutegravir in HIV-1-infected South African patients initiating or failing first-line ART

Objectives

The WHO recently recommended the use of a new first-line ART containing dolutegravir. We investigated the efficacy of NRTI backbones (tenofovir or abacavir with a cytosine analogue) in low- and middle-income countries where there is significant prior exposure to antiretrovirals and drug resistance to NRTIs.

Methods

Within the treatment-as-prevention study in South Africa, we selected participants with available next-generation sequencing (NGS) data for the HIV-1 pol gene at trial entry; they were either ART initiators (n = 1193) or already established on ART (n = 94). NGS of the HIV-1 pol gene was carried out using MiSeq technology; reverse transcriptase drug resistance mutations (DRMs) were detected at 5% (DRM5%) and 20% (DRM20%) for all 1287 participants. Genotypic susceptibility was assessed using the Stanford HIVDB resistance interpretation algorithm.

Results

NRTI DRM20% and DRM5% were detected among 5/1193 (0.4%) and 9/1193 (0.8%) of ART initiators, respectively. There was tenofovir exposure in 73/94 (77.7%) of those established on ART, with full susceptibility to abacavir in 57/94 (60.6%) and 56/94 (59.6%) for DRM20% and DRM5%, respectively, while 67/94 (71.3%) and 64/94 (68.1%) were fully susceptible to tenofovir, respectively. The differences between tenofovir and abacavir were not statistically significant at the 20% or 5% variant level (P = 0.16 and 0.29, respectively). NGS detection of variants at the 5% level increased detection of K65R in both naive and treated groups. One of 607 integrase sequences carried a DRM20% (Q148R).

Conclusions

Dolutegravir with a cytosine analogue plus tenofovir or abacavir appears to have similar efficacy in South Africans naive to ART. NGS should be considered in HIV drug resistance surveillance.

Diagnosis of Human Immunodeficiency Virus Infection

HIV diagnostics have played a central role in the remarkable progress in identifying, staging, initiating, and monitoring infected individuals on life-saving antiretroviral therapy. They are also useful in surveillance and outbreak responses, allowing for assessment of disease burden and identification of vulnerable populations and transmission "hot spots," thus enabling planning, appropriate interventions, and allocation of appropriate funding. HIV diagnostics are critical in achieving epidemic control and require a hybrid of conventional laboratory-based diagnostic tests and new technologies, including point-of-care (POC) testing, to expand coverage, increase access, and positively impact patient management. In this review, we provide (i) a historical perspective on the evolution of HIV diagnostics (serologic and molecular) and their interplay with WHO normative guidelines, (ii) a description of the role of conventional and POC testing within the tiered laboratory diagnostic network, (iii) information on the evaluations and selection of appropriate diagnostics, (iv) a description of the quality management systems needed to ensure reliability of testing, and (v) strategies to increase access while reducing the time to return results to patients. Maintaining the central role of HIV diagnostics in programs requires periodic monitoring and optimization with quality assurance in order to inform adjustments or alignment to achieve epidemic control.

Characterization of minority HIV-1 drug resistant variants in the United Kingdom following the verification of a deep sequencing-based HIV-1 genotyping and tropism assay

BACKGROUND

The widespread global access to antiretroviral drugs has led to considerable reductions in morbidity and mortality but, unfortunately, the risk of virologic failure increases with the emergence, and potential transmission, of drug resistant viruses. Detecting and quantifying HIV-1 drug resistance has therefore become the standard of care when designing new antiretroviral regimens. The sensitivity of Sanger sequencing-based HIV-1 genotypic assays is limited by its inability to identify minority members of the quasispecies, i.e., it only detects variants present above ~ 20% of the viral population, thus, failing to detect minority variants below this threshold. It is clear that deep sequencing-based HIV-1 genotyping assays are an important step change towards accurately monitoring HIV-infected individuals.

METHODS

We implemented and verified a clinically validated HIV-1 genotyping assay based on deep sequencing (DEEPGEN™) in two clinical laboratories in the United Kingdom: St. George's University Hospitals Healthcare NHS Foundation Trust (London) and at NHS Lothian (Edinburgh), to characterize minority HIV-1 variants in 109 plasma samples from ART-naïve or -experienced individuals.

RESULTS

Although subtype B HIV-1 strains were highly prevalent (44%, 48/109), most individuals were infected with non-B subtype viruses (i.e., A1, A2, C, D, F1, G, CRF02_AG, and CRF01_AE). DEEPGEN™ was able to accurately detect drug resistance-associated mutations not identified using standard Sanger sequencing-based tests, which correlated significantly with patient's antiretroviral treatment histories. A higher proportion of minority PI-, NRTI-, and NNRTI-resistance mutations was detected in NHS Lothian patients compared to individuals from St. George's, mainly M46I/L and I50 V (associated with PIs), D67 N, K65R, L74I, M184 V/I, and K219Q (NRTIs), and L100I (NNRTIs). Interestingly, we observed an inverse correlation between intra-patient HIV-1 diversity and CD4⁺ T cell counts in the NHS Lothian patients.

CONCLUSIONS

This is the first study evaluating the transition, training, and implementation of DEEPGEN™ between three clinical laboratories in two different countries. More importantly, we were able to characterize the HIV-1 drug resistance profile (including minority variants), coreceptor tropism, subtyping, and intra-patient viral diversity in patients from the United Kingdom, providing a rigorous foundation for basing clinical decisions on highly sensitive and cost-effective deep sequencing-based HIV-1 genotyping assays in the country.

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.

Selective resistance profiles emerging in patient-derived clinical isolates with cabotegravir, bictegravir, dolutegravir, and elvitegravir

Background

Integrase strand transfer inhibitors (INSTIs) are recommended for first-line HIV therapy based on their relatively high genetic barrier to resistance. Although raltegravir (RAL) and elvitegravir (EVG) resistance profiles are well-characterized, resistance patterns for dolutegravir (DTG), bictegravir (BIC), and cabotegravir (CAB) remain largely unknown. Here, in vitro drug selections compared the development of resistance to DTG, BIC, CAB, EVG and RAL using clinical isolates from treatment-naïve primary HIV infection (PHI) cohort participants (n = 12), and pNL4.3 recombinant strains encoding patient-derived Integrase with (n = 5) and without (n = 5) the E157Q substitution.

Results

Patient-derived viral isolates were serially passaged in PHA-stimulated cord blood mononuclear cells in the presence of escalating concentrations of INSTIs over the course of 36–46 weeks. Drug resistance arose more rapidly in primary clinical isolates with EVG (12/12), followed by CAB (8/12), DTG (8/12) and BIC (6/12). For pNL4.3 recombinant strains encoding patient-derived integrase, the comparative genetic barrier to resistance was RAL > EVG > CAB > DTG and BIC. The E157Q substitution in integrase delayed the advent of resistance to INSTIs. With EVG, T66I/A, E92G/V/Q, T97A or R263K (n = 16, 3, 2 and 1, respectively) arose by weeks 8–16, followed by 1–4 accessory mutations, conferring high-level resistance (> 100-fold) by week 36. With DTG and BIC, solitary R263K (n = 27), S153F/Y (n = 7) H51Y (n = 2), Q146 R (n = 3) or S147G (n = 1) mutations conferred low-level (< 3-fold) resistance at weeks 36–46. Similarly, most CAB selections (n = 18) resulted in R263K, S153Y, S147G, H51Y, or Q146L solitary mutations. However, three CAB selections resulted in Q148R/K followed by secondary mutations conferring high-level cross-resistance to all INSTIs. EVG-resistant viruses (T66I/R263K, T66I/E157Q/R263K, and S153A/R263K) retained residual susceptibility when switched to DTG, BIC or CAB, losing T66I by week 27. Two EVG-resistant variants developed resistance to DTG, BIC and CAB through the additional acquisition of E138A/Q148R and S230N, respectively. One EVG-resistant variant (T66I) acquired L74M/G140S/S147G, L74M/E138K/S147G and H51Y with DTG CAB and BIC, respectively.

Conclusions

Second generation INSTIs show a higher genetic barrier to resistance than EVG and RAL. The potency of CAB was lower than BIC and DTG. The development of Q148R/K with CAB can result in high-level cross-resistance to all INSTIs.

Evaluating the accuracy and sensitivity of detecting minority HIV-1 populations by Illumina next-generation sequencing

The accuracy and sensitivity of deep sequencing were assessed using viral standards (pNL4-3 and pLAI.2) of both DNA and RNA. The sequencing accuracy did not depend on the type of nucleic acid, but critically depended on the number of reads and threshold of sensitivity to minor viral populations. With coverage of more than 236 reads, the accuracy of viral RNA sequencing was equal to or exceeded 99.9%, with a sensitivity threshold to minor nucleotides of 20%. When the sensitivity threshold was below 1%, reduced accuracy dynamics were clearly visible even when the coverage was massive (more than 9.000 reads). It was found that the floating sensitivity threshold allowed the sequencing accuracy to be maintained at an acceptable level in cases of low coverage (less than 1.500-2.000) of reads. These results indicate the quality that can be expected with a specific number of reads and sensitivity threshold. Deep sequencing is a very powerful tool that can significantly improve the value of study results, but despite its superior performance, it should be used with caution regarding its sensitivity to minor populations below 1%.

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.

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

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

Next-Generation Human Immunodeficiency Virus Sequencing for Patient Management and Drug Resistance Surveillance

High-quality, simplified, and low-cost human immunodeficiency virus (HIV) drug resistance tests that are able to provide timely actionable HIV resistance data at individual, population, and programmatic levels are needed to confront the emerging drug-resistant HIV epidemic. Next-generation sequencing technologies embedded in automated cloud-computing analysis environments are ideally suited for such endeavor. Whereas NGS can reduce costs over Sanger sequencing, automated analysis pipelines make NGS accessible to molecular laboratories regardless of the available bioinformatic skills. They can also produce highly structured, high-quality data that could be examined by healthcare officials and program managers on a real-time basis to allow timely public health action. Here we discuss the opportunities and challenges of such an approach.

Prevalence and clinical impacts of HIV-1 intersubtype recombinants in Uganda revealed by near-full-genome population and deep sequencing approaches

OBJECTIVES

HIV-1 subtypes A1 and D cocirculate in a rural community in Mbarara, Uganda. This study examines HIV-1 intersubtype recombination in this community under a full-genome sequencing context. We aim to estimate prevalence, examine time trends, and test for clinical correlates and outcomes associated with intersubtype recombinants.

METHODS

Near-full-genome HIV-1 Sanger sequence data were collected from plasma samples of 504 treatment-naïve individuals, who then received protease inhibitor or nonnucleoside reverse transcriptase inhibitor-containing regimens and were monitored for up to 7.5 years. Subtypes were inferred by Los Alamos Recombinant Identification Program (RIP) 3.0 and compared with Sanger/REGA and MiSeq/RIP. 'Nonrecombinants' and 'recombinants' infections were compared in terms of pretherapy viral load, CD4 cell count, posttherapy time to virologic suppression, virologic rebound, first CD4 rise above baseline and sustained CD4 recovery.

RESULTS

Prevalence of intersubtype recombinants varied depending on the genomic region examined: gag (15%), prrt (11%), int (8%), vif (10%), vpr (2%), vpu (9%), GP120 (8%), GP41 (18%), and nef (4%). Of the 200 patients with near-full-genome data, prevalence of intersubtype recombination was 46%; the most frequently observed recombinant was A1-D (25%). Sanger/REGA and MiSeq/RIP yielded generally consistent results. Phylogenetic tree revealed most recombinants did not share common ancestors. No temporal trend was observed (all P > 0.1). Subsequent subtype switches were detected in 27 of 143 (19%) study participants with follow-up sequences. Nonrecombinant versus recombinants infections were not significantly different in any pre nor posttherapy clinical correlates examined (all P > 0.2).

CONCLUSION

Intersubtype recombination was highly prevalent (46%) in Uganda if the entire HIV genome was considered, but was neither associated with clinical correlates nor therapy outcomes.

Promises and pitfalls of Illumina sequencing for HIV resistance genotyping

Genetic sequencing ("genotyping") plays a critical role in the modern clinical management of HIV infection. This virus evolves rapidly within patients because of its error-prone reverse transcriptase and short generation time. Consequently, HIV variants with mutations that confer resistance to one or more antiretroviral drugs can emerge during sub-optimal treatment. There are now multiple HIV drug resistance interpretation algorithms that take the region of the HIV genome encoding the major drug targets as inputs; expert use of these algorithms can significantly improve to clinical outcomes in HIV treatment. Next-generation sequencing has the potential to revolutionize HIV resistance genotyping by lowering the threshold that rare but clinically significant HIV variants can be detected reproducibly, and by conferring improved cost-effectiveness in high-throughput scenarios. In this review, we discuss the relative merits and challenges of deploying the Illumina MiSeq instrument for clinical HIV genotyping.

HIV-1 drug resistance and resistance testing

The global scale-up of antiretroviral (ARV) therapy (ART) has led to dramatic reductions in HIV-1 mortality and incidence. However, HIV drug resistance (HIVDR) poses a potential threat to the long-term success of ART and is emerging as a threat to the elimination of AIDS as a public health problem by 2030. In this review we describe the genetic mechanisms, epidemiology, and management of HIVDR at both individual and population levels across diverse economic and geographic settings. To describe the genetic mechanisms of HIVDR, we review the genetic barriers to resistance for the most commonly used ARVs and describe the extent of cross-resistance between them. To describe the epidemiology of HIVDR, we summarize the prevalence and patterns of transmitted drug resistance (TDR) and acquired drug resistance (ADR) in both high-income and low- and middle-income countries (LMICs). We also review to two categories of HIVDR with important public health relevance: (i) pre-treatment drug resistance (PDR), a World Health Organization-recommended HIVDR surveillance metric and (ii) and pre-exposure prophylaxis (PrEP)-related drug resistance, a type of ADR that can impact clinical outcomes if present at the time of treatment initiation. To summarize the implications of HIVDR for patient management, we review the role of genotypic resistance testing and treatment practices in both high-income and LMIC settings. In high-income countries where drug resistance testing is part of routine care, such an understanding can help clinicians prevent virological failure and accumulation of further HIVDR on an individual level by selecting the most efficacious regimens for their patients. Although there is reduced access to diagnostic testing and to many ARVs in LMIC, understanding the scientific basis and clinical implications of HIVDR is useful in all regions in order to shape appropriate surveillance, inform treatment algorithms, and manage difficult cases.

Low-Frequency Drug Resistance in HIV-Infected Ugandans on Antiretroviral Treatment Is Associated with Regimen Failure

Most patients failing antiretroviral treatment in Uganda continue to fail their treatment regimen even if a dominant drug-resistant HIV-1 genotype is not detected. In a recent retrospective study, we observed that approximately 30% of HIV-infected individuals in the Joint Clinical Research Centre (Kampala, Uganda) experienced virologic failure with a susceptible HIV-1 genotype based on standard Sanger sequencing. Selection of minority drug-resistant HIV-1 variants (not detectable by Sanger sequencing) under antiretroviral therapy pressure can lead to a shift in the viral quasispecies distribution, becoming dominant members of the virus population and eventually causing treatment failure. Here, we used a novel HIV-1 genotyping assay based on deep sequencing (DeepGen) to quantify low-level drug-resistant HIV-1 variants in 33 patients failing a first-line antiretroviral treatment regimen in the absence of drug-resistant mutations, as screened by standard population-based Sanger sequencing. Using this sensitive assay, we observed that 64% (21/33) of these individuals had low-frequency (or minority) drug-resistant variants in the intrapatient HIV-1 population, which correlated with treatment failure. Moreover, the presence of these minority HIV-1 variants was associated with higher intrapatient HIV-1 diversity, suggesting a dynamic selection or fading of drug-resistant HIV-1 variants from the viral quasispecies in the presence or absence of drug pressure, respectively. This study identified low-frequency HIV drug resistance mutations by deep sequencing in Ugandan patients failing antiretroviral treatment but lacking dominant drug resistance mutations as determined by Sanger sequencing methods. We showed that these low-abundance drug-resistant viruses could have significant consequences for clinical outcomes, especially if treatment is not modified based on a susceptible HIV-1 genotype by Sanger sequencing. Therefore, we propose to make clinical decisions using more sensitive methods to detect minority HIV-1 variants.

Genetic Variability of HIV-1 for Drug Resistance Assay Development

A hybridization-based point-of-care (POC) assay for HIV-1 drug resistance would be useful in low- and middle-income countries (LMICs) where resistance testing is not routinely available. The major obstacle in developing such an assay is the extreme genetic variability of HIV-1. We analyzed 27,203 reverse transcriptase (RT) sequences from the Stanford HIV Drug Resistance Database originating from six LMIC regions. We characterized the variability in a 27-nucleotide window surrounding six clinically important drug resistance mutations (DRMs) at positions 65, 103, 106, 181, 184, and 190. The number of distinct codons at each DRM position ranged from four at position 184 to 11 at position 190. Depending on the mutation, between 11 and 15 of the 24 flanking nucleotide positions were variable. Nonetheless, most flanking sequences differed from a core set of 10 flanking sequences by just one or two nucleotides. Flanking sequence variability was also lower in each LMIC region compared with overall variability in all regions. We also describe an online program that we developed to perform similar analyses for mutations at any position in RT, protease, or integrase.

HIV Whole-Genome Sequencing Now: Answering Still-Open Questions

Diversity, evolution, and epidemiology of HIV are directly relevant to HIV transmission and pathogenesis; hence, they play a key role in antiretroviral treatment and vaccine design. Global HIV whole-genome sequencing would provide a treasure chest of data to answer many questions still open in these fields. An article by Berg et al. in this issue of theJournal of Clinical Microbiologydescribes a universal strategy to amplify and sequence heterogeneous HIV whole genomes (M. G. Berg, J. Yamaguchi, E. Alessandri-Gradt, R. W. Tell, J.-C. Plantier, and C. A. Brennan, J Clin Microbiol 54:868-882, 2016,http://dx.doi.org/10.1128/JCM.02479-15).

HIV-1 Drug Resistance Mutations: Potential Applications for Point-of-Care Genotypic Resistance Testing

The increasing prevalence of acquired and transmitted HIV-1 drug resistance is an obstacle to successful antiretroviral therapy (ART) in the low- and middle-income countries (LMICs) hardest hit by the HIV-1 pandemic. Genotypic drug resistance testing could facilitate the choice of initial ART in areas with rising transmitted drug resistance (TDR) and enable care-providers to determine which individuals with virological failure (VF) on a first- or second-line ART regimen require a change in treatment. An inexpensive near point-of-care (POC) genotypic resistance test would be useful in settings where the resources, capacity, and infrastructure to perform standard genotypic drug resistance testing are limited. Such a test would be particularly useful in conjunction with the POC HIV-1 viral load tests that are currently being introduced in LMICs. A POC genotypic resistance test is likely to involve the use of allele-specific point mutation assays for detecting drug-resistance mutations (DRMs). This study proposes that two major nucleoside reverse transcriptase inhibitor (NRTI)-associated DRMs (M184V and K65R) and four major NNRTI-associated DRMs (K103N, Y181C, G190A, and V106M) would be the most useful for POC genotypic resistance testing in LMIC settings. One or more of these six DRMs was present in 61.2% of analyzed virus sequences from ART-naïve individuals with intermediate or high-level TDR and 98.8% of analyzed virus sequences from individuals on a first-line NRTI/NNRTI-containing regimen with intermediate or high-level acquired drug resistance. The detection of one or more of these DRMs in an ART-naïve individual or in a individual with VF on a first-line NRTI/NNRTI-containing regimen may be considered an indication for a protease inhibitor (PI)-containing regimen or closer virological monitoring based on cost-effectiveness or country policy.