Minority variants of drug-resistant HIV

HIV-1 Drug Resistance Detected by Next-Generation Sequencing among ART-Naïve Individuals: A Systematic Review and Meta-Analysis

BACKGROUND

There are an increasing number of articles focused on the prevalence and clinical impact of pretreatment HIV drug resistance (PDR) detected by Sanger sequencing (SGS). PDR may contribute to the increased likelihood of virologic failure and the emergence of new resistance mutations. As SGS is gradually replaced by next-generation sequencing (NGS), it is necessary to assess the levels of PDR using NGS in ART-naïve patients systematically. NGS can detect the viral variants (low-abundance drug-resistant HIV-1 variants (LA-DRVs)) of virus quasi-species at levels below 20% that SGS may fail to detect. NGS has the potential to optimize current HIV drug resistance surveillance methods and inform future research directions. As the NGS technique has high sensitivity, it is highly likely that the level of pretreatment resistance would be underestimated using conventional techniques.

METHODS

For the systematic review and meta-analysis, we searched for original studies published in PubMed, Web of Science, Scopus, and Embase before 30 March 2023 that focused exclusively on the application of NGS in the detection of HIV drug resistance. Pooled prevalence estimates were calculated using a random effects model using the 'meta' package in R (version 4.2.3). We described drug resistance detected at five thresholds (>1%, 2%, 5%, 10%, and 20% of virus quasi-species). Chi-squared tests were used to analyze differences between the overall prevalence of PDR reported by SGS and NGS.

RESULTS

A total of 39 eligible studies were selected. The studies included a total of 15,242 ART-naïve individuals living with HIV. The prevalence of PDR was inversely correlated with the mutation detection threshold. The overall prevalence of PDR was 29.74% at the 1% threshold, 22.43% at the 2% threshold, 15.47% at the 5% threshold, 12.95% at the 10% threshold, and 11.08% at the 20% threshold. The prevalence of PDR to INSTIs was 1.22% (95%CI: 0.58-2.57), which is the lowest among the values for all antiretroviral drugs. The prevalence of LA-DRVs was 9.45%. At the 2% and 20% detection threshold, the prevalence of PDR was 22.43% and 11.08%, respectively. Resistance to PIs and INSTIs increased 5.52-fold and 7.08-fold, respectively, in those with a PDR threshold of 2% compared with those with PDR at 20%. However, resistance to NRTIs and NNRTIs increased 2.50-fold and 2.37-fold, respectively. There was a significant difference between the 2% and 5% threshold for detecting HIV drug resistance. There was no statistically significant difference between the results reported by SGS and NGS when using the 20% threshold for reporting resistance mutations.

CONCLUSION

In this study, we found that next-generation sequencing facilitates a more sensitive detection of HIV-1 drug resistance than SGS. The high prevalence of PDR emphasizes the importance of baseline resistance and assessing the threshold for optimal clinical detection using NGS.

Performance characteristics of Allele-Specific PCR (ASPCR) in detecting drug resistance mutations among non-B HIV-1 Variants

Allele-Specific Polymerase Chain Reaction (ASPCR) is an affordable point-mutation assay whose validation could improve the detection of HIV-1 drug resistance mutations (DRMs) in resource-limited settings (RLS). We assessed the performance of ASPCR onforty-four non-B HIV-1 plasma samples from patients who were ARV treated in failure in N'Djamena-Chad. Viral RNA was reverse-transcribed and amplified using LightCycler® FastStart DNA MasterPLUS SYBR Green I. Detection of six major DRMs (K70R, K103N, Y181C, M184V, T215F, T215Y) was evaluated on Roche LightCycler®480 automated system (with dilutions 0.01-100%). ASPCR-results were compared to Sanger-sequencing (gold-standard). Correlations of mutation curves were excellent (R2 >0.97); all DRMs were detected with desirable mutant/wild-type threshold differences (ΔCt≥9) except K70R(ΔCtK70R=6; ΔCtK103N=13; ΔCtM184V=9; ΔCtT215F=12; ΔCtT215Y=12; ΔCtY181C=9) and positive controls were below required thresholds. Also, ASPCR reproducibility on DRMs was assessed by using dilutions of intra-assay and inter-assay coefficient of variations respectively with a threshold of less than 50(i.e.<0.50 variation) which are;: K70R (0.02-0.28 vs. 0.12-0.37), K103N (0.08-0.42 vs. 0.12-0.37), Y181C (0.12-0.39 vs. 0.31-0.37), M184V (0.13-0.39 vs. 0.23-0.42), T215F (0.05-0.43 vs. 0.04-0.45) and T215Y (0.13-0.41 vs. 0.19-0.41). DRM detection-rate by ASPCR vs Sanger was respectively: M184V (63.6% vs. 38.6%); T215F (18.1% vs. 9.1%); T215Y (6.8% vs. 2.3%); K70R (4.5% vs. 2.3%). K103N (22.7% vs. 13.6%); Y181C (13.6% vs. 11.4%). Correlations of mutation curves were excellent (R2 >0.97); all DRMs were detected with desirable mutant/wild-type threshold differences (ΔCt≥9) except K70R(ΔCtK70R=6; ΔCtK103N=13; ΔCtM184V=9; ΔCtT215F=12; ΔCtT215Y=12; ΔCtY181C=9) and positive controls were below required thresholds. Also, ASPCR reproducibility on DRMs was assessed by using dilutions of intra-assay and inter-assay coefficient of variations respectively with a threshold of less than 50(i.e.<0.50 variation) which are;: K70R (0.02-0.28 vs. 0.12-0.37), K103N (0.08-0.42 vs. 0.12-0.37), Y181C (0.12-0.39 vs. 0.31-0.37), M184V (0.13-0.39 vs. 0.23-0.42), T215F (0.05-0.43 vs. 0.04-0.45) and T215Y (0.13-0.41 vs. 0.19-0.41). DRM detection-rate by ASPCR vs Sanger was respectively: M184V (63.6% vs. 38.6%); T215F (18.1% vs. 9.1%); T215Y (6.8% vs. 2.3%); K70R (4.5% vs. 2.3%). K103N (22.7% vs. 13.6%); Y181C (13.6% vs. 11.4%). ASPCR appears more efficient for detecting DRMs on diverse HIV-1 non-B circulating in RLS like Chad.

Application of the ViroKey® SQ FLEX assay for detection of cytomegalovirus antiviral resistance

BACKGROUND

Cytomegalovirus (CMV) is a viral infection which establishes lifelong latency, often reactivating and causing disease in immunosuppressed individuals, including haematopoietic stem cell transplant (HSCT) recipients. Treatment can be problematic due to antiviral resistance which substantially increases the risk of patient mortality. Diagnostic testing capabilities for CMV antiviral resistance in Australia and elsewhere have traditionally relied on gene-specific Sanger sequencing approaches, however, are now being superseded by next generation sequencing protocols.

OBJECTIVE

Provide a snapshot of local mutations and explore the feasibility of the ViroKeyࣨ® SQ FLEX Genotyping Assay (Vela Diagnostics Pty Ltd) by examining sequencing success.

METHOD

Performed sequencing on adult (n = 38) and paediatric (n = 81) plasma samples, over a large range of viral loads (above and below the assay recommended threshold of ≥1,000 International Units (IU)/mL; noting most of our paediatric samples have loads <1,000 IU/mL).

RESULTS

Eleven test runs (including three repeat runs; 14 to 15 samples per run) were conducted, and four runs were deemed valid. The overall individual sample success rate for the four evaluable test runs was 71.2% (42/59 samples); 80.4% (37/46) samples ≥1,000 IU/mL were valid. Ten clinically important antiviral resistance mutations were detected, the most common being A594V in the UL97 gene, found in 6 (5%) samples.

CONCLUSIONS

A range of technical issues were experienced, however with improvement this platform could be a useful addition to routine pathology workflows, providing timely antiviral resistance results for patients undergoing HSCT.

Full-spectrum HIV drug resistance mutation detection by high-resolution complete pol gene sequencing

BACKGROUND

Drug resistance mutation testing is a key element for HIV clinical management, informing effective treatment regimens. However, resistance screening in current clinical practice is limited in reporting linked cross-class resistance mutations and minority variants, both of which may increase the risk of virological failure.

METHODS

To address these limitations, we obtained 358 full-length pol gene sequences from 52 specimens of 20 HIV infected individuals by combining microdroplet amplification, unique molecular identifier (UMI) labeling, and long-read high-throughput sequencing.

RESULTS

We conducted a rigorous assessment of the accuracy of our pipeline for precision drug resistance mutation detection, verifying that a sequencing depth of 35 high-throughput reads achieved complete, error-free pol gene sequencing. We detected 26 distinct drug resistance mutations to Protease Inhibitors (PIs), Nucleoside Reverse Transcriptase Inhibitors (NRTIs), Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs), and Integrase Strand Transfer Inhibitors (INSTIs). We detected linked cross-class drug resistance mutations (PI+NRTI, PI+NNRTI, and NRTI+NNRTI) that confer cross-resistance to multiple drugs in different classes. Fourteen different types of minority mutations were also detected with frequencies ranging from 3.2% to 19%, and the presence of these mutations was verified by Sanger reference sequencing. We detected a putative transmitted drug resistance mutation (TDRM) in one individual that persisted for over seven months from the first sample collected at the acute stage of infection prior to seroconversion.

CONCLUSIONS

Our comprehensive drug resistance mutation profiling can advance clinical practice by reporting mutation linkage and minority variants to better guide antiretroviral therapy options.

Towards Next-Generation Sequencing for HIV-1 Drug Resistance Testing in a Clinical Setting

The HIV genotypic resistance test (GRT) is a standard of care for the clinical management of HIV/AIDS patients. In recent decades, population or Sanger sequencing has been the foundation for drug resistance monitoring in clinical settings. However, the advent of high-throughput or next-generation sequencing has caused a paradigm shift towards the detection and characterization of low-abundance covert mutations that would otherwise be missed by population sequencing. This is clinically significant, as these mutations can potentially compromise the efficacy of antiretroviral therapy, causing poor virologic suppression. Therefore, it is important to develop a more sensitive method so as to reliably detect clinically actionable drug-resistant mutations (DRMs). Here, we evaluated the diagnostic performance of a laboratory-developed, high-throughput, sequencing-based GRT using 103 archived clinical samples that were previously tested for drug resistance using population sequencing. As expected, high-throughput sequencing found all the DRMs that were detectable by population sequencing. Significantly, 78 additional DRMs were identified only by high-throughput sequencing, which is statistically significant based on McNemar's test. Overall, our results complement previous studies, supporting the notion that the two methods are well correlated, and the high-throughput sequencing method appears to be an excellent alternative for drug resistance testing in a clinical setting.

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.

Droplet Microfluidics for High-Resolution Virology

Microfluidics has enabled a new era of cellular and molecular assays due to the small length scales, parallelization, and the modularity of various analysis and actuation functions. Droplet microfluidics, in particular, has been instrumental in providing new tools for biology with its ability to quickly and reproducibly generate drops that act as individual reactors. A notable beneficiary of this technology has been single-cell RNA sequencing, which has revealed new heterogeneities and interactions for the fundamental unit of life. However, viruses far surpass the diversity of cellular life, affect the dynamics of all ecosystems, and are a chronic source of global health crises. Despite their impact on the world, high-throughput and high-resolution viral profiling has been difficult, with conventional methods being limited to population-level averaging, large sample volumes, and few cultivable hosts. Consequently, most viruses have not been identified and studied. Droplet microfluidics holds the potential to address many of these limitations and offers new levels of sensitivity and throughput for virology. This Feature highlights recent efforts that have applied droplet microfluidics to the detection and study of viruses, including for diagnostics, virus-host interactions, and cell-independent virus assays. In combination with traditional virology methods, droplet microfluidics should prove a potent tool toward achieving a better understanding of the most abundant biological species on Earth.

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.

Diversity of HIV-1 subtypes and transmitted drug-resistance mutations among minority HIV-1 variants in a Turkish cohort

BACKGROUND

The World Health Organization (WHO) recommends the surveillance of transmitted drug resistance mutations (TDRMs) to ensure the effectiveness and sustainability of HIV treatment programs.

OBJECTIVE

Our aim was to determine the TDRMs and evaluate the distribution of HIV-1 subtypes using and compared next-generation sequencing (NGS) and Sanger-based sequencing (SBS) in a cohort of 44 antiretroviral treatment-naïve patients.

METHODS

All samples that were referred to the microbiology laboratory for HIV drug resistance analysis between December 2016 and February 2018 were included in the study. After exclusions, 44 treatment-naive adult patients with a viral load of >1000 copies/mL were analyzed. DNA sequencing for reverse transcriptase and protease regions was performed using both DeepChek ABL single round kit and Sanger-based ViroSeq HIV-1 Genotyping System. The mutations and HIV-1 subtypes were analyzed using the Stanford HIVdb version 8.6.1 Genotypic Resistance software, and TDRMs were assessed using the WHO surveillance drug-resistance mutation database. HIV-1 subtypes were confirmed by constructing a maximum-likelihood phylogenetic tree using Los Alamos IQ-Tree software.

RESULTS

NGS identified nucleos(t)ide reverse transcriptase inhibitor (NRTI)-TDRMs in 9.1% of the patients, non-nucleos(t)ide reverse transcriptase inhibitor (NNRTI)-TDRMs in 6.8% of the patients, and protease inhibitor (PI)-TDRMs in 18.2% of the patients at a detection threshold of ≥1%. Using SBS, 2.3% and 6.8% of the patients were found to have NRTI- and NNRTI-TDRMs, respectively, but no major PI mutations were detected. M41L, L74I, K65R, M184V, and M184I related to NRTI, K103N to NNRTI, and N83D, M46I, I84V, V82A, L24I, L90M, I54V to the PI sites were identified using NGS. Most mutations were found in low-abundance (frequency range: 1.0% - 4.7%) HIV-1 variants, except M41L and K103N. The subtypes of the isolates were found as follows; 61.4% subtype B, 18.2% subtype B/CRF02_AG recombinant, 13.6% subtype A, 4.5% CRF43_02G, and 2.3% CRF02_AG. All TDRMs, except K65R, were detected in HIV-1 subtype B isolates.

CONCLUSION

The high diversity of protease site TDRMs in the minority HIV-1 variants and prevalence of CRFs were remarkable in this study. All minority HIV-1 variants were missed by conventional sequencing. TDRM prevalence among minority variants appears to be decreasing over time at our center.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

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.

Detection of human immunodeficiency virus Type 1 phylogenetic clusters with multidrug resistance mutations among 2011 to 2017 blood donors from the highly endemic Northern Brazilian Amazon

BACKGROUND

This study describes transmitted drug resistance (TDR) in blood donors diagnosed with human immunodeficiency virus Type 1 (HIV-1) infection from 2011 to 2017 in three reference public blood centers from the Northern Brazilian Amazon.

STUDY DESIGN AND METHODS

This was a cross-sectional study on HIV-positive blood donors from HEMOAM, Manaus, Amazonas, AM (n = 198); HEMERON, Porto Velho, Rondônia, RO (n = 20); and HEMORAIMA, Boa Vista, Roraima, RR (n = 9). HIV-1 pol sequences (protease, reverse transcriptase) were analyzed for drug resistance mutations (DRMs) using the Calibrated Population Resistance tool (Stanford). TDR/DRM clusters were investigated by phylogenetic analysis after removing positions associated with drug resistance of Subtype B sequences from untreated and treated subjects from Northern Brazil.

RESULTS

Transmitted drug resistance/DRM in blood donors was 11% (25 of 227), all of them from HEMOAM. Most blood donors with TDR/DRM had multiple and similar DRMs. Nonnucleoside reverse transcriptase inhibitor (NNRTI) mutations predominated (10.1%), followed by nucleoside reverse transcriptase inhibitor (NRTI) mutations (5.3%) and protease inhibitor mutations (0.4%). Dual-class NNRTI/NRTI mutations represented 4.8%. Three highly supported Subtype B monophyletic clades mostly composed by individuals from Amazonas with TDR/DRM mutations were identified. The largest transmission cluster contained 10 sequences, eight from HEMOAM and two sequences described previously (one from a treated subject from Amazonas and the other one from Roraima). This cluster was characterized by NRTI (D67N, T69D, T215S/F/L, K219Q) and NNRTI (K101H, K103 N, G190A) mutations. The other two transmission clades comprised only three and two sequences from HEMOAM sharing the E138A NNRTI mutation.

CONCLUSIONS

The identification of transmission clusters of multidrug-resistant viruses in blood donors from Amazonas highlight the need of continued monitoring of TDR/DRM and the importance of pretreatment genotyping in the highly endemic Amazonas state.

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.

Surveillance of transmitted HIV-1 antiretroviral drug resistance in the context of decentralized HIV care in Senegal and the Ebola outbreak in Guinea

Objectives

Disruption in HIV care provision may enhance the development and spread of drug resistance due to inadequate antiretroviral therapy. This study thus determined the prevalence of HIV-1 transmitted drug resistance (TDR) in settings of decentralized therapy and care in Senegal and, the Ebola outbreak in Guinea. Antiretroviral-naïve patients were enrolled following a modified WHO TDR Threshold Survey method, implemented in Senegal (January–March 2015) and Guinea (August–September 2015). Plasma and dried blood spots specimens, respectively from Senegalese (n = 69) and Guinean (n = 50) patients, were collected for direct sequencing of HIV-1 pol genes. The Stanford Calibrated Population Resistance program v6.0 was used for Surveillance Drug Resistance Mutations (SDRMs).

Results

Genotyping was successful from 54/69 (78.2%) and 31/50 (62.0%) isolates. In Senegal, TDR prevalence was 0% (mean duration since HIV diagnosis 4.08 ± 3.53 years). In Guinea, two patients exhibited SDRMs M184V (NRTI), T215F (TAM) and, G190A (NNRTI), respectively. TDR prevalence at this second site, however, could not be ascertained because of low sample size. Phylogenetic inference confirmed CRF02_AG predominance in Senegal (62.96%) and Guinea (77.42%). TDR prevalence in Senegal remains extremely low suggesting improved control measures. Continuous surveillance in both settings is mandatory and, should be done closest to diagnosis/transmission time and with larger sample size.

Focus on recently developed assays for detection of resistance/sensitivity to reverse transcriptase inhibitors

The biology of HIV is rather complex due to high rate of replication, frequent recombination, and introduction of mutations. This gives rise to a number of distinct variants referred as quasispecies. In addition, the latency within reservoir allows the periodic reactivation of virus replication. The rapid replication of HIV allows immune response escape and establishment of resistance to therapy that can be acquired through drug selection and/or transmitted among individuals. This prompted, over the years, the development of a range of assays aimed to determine drug resistance and sensitivity, to be used both in clinical practice and in antiviral research. Reverse transcriptase (RT) inhibitors have an eminent place among the anti-HIV drugs, being constantly present from the beginning until today in the most commonly used antiviral regimens. This mini-review seeks to provide an up-to-date overview of recent efforts in developing even more reliable and simple methods, of both genotypic and phenotypic types, for specifically detecting drug resistance and sensitivity to RT inhibitors.

Sensitive detection of HIV-1 resistance to Zidovudine and impact on treatment outcomes in low- to middle-income countries

Background

Thymidine analogs, namely AZT (Zidovudine or Retrovir™) and d4T (Stavudine or Zerit™) are antiretroviral drugs still employed in over 75% of first line combination antiretroviral therapy (cART) in Kampala, Uganda despite aversion to prescribing these drugs for cART in high income countries due in part to adverse events. For this study, we explored how the continued use of these thymidine analogs in cART could impact emergence of drug resistance and impact on future treatment success in Uganda, a low-income country.

Methods

We examined the drug resistance genotypes by Sanger sequencing of 262 HIV-infected patients failing a first line combined antiretroviral treatment containing either AZT or d4T, which represents approximately 5% of the patients at the Joint Clinical Research Center receiving a AZT or d4T containing treatment. Next generation sequencing (DEEPGEN™HIV) and multiplex oligonucleotide ligation assays (AfriPOLA) were then performed on a subset of patient samples to detect low frequency drug resistant mutations. CD4 cell counts, viral RNA loads, and treatment changes were analyzed in a cohort of treatment success and failures.

Results

Over 80% of patients failing first line AZT/d4T-containing cART had predicted drug resistance to 3TC (Lamivudine) and non-nucleoside RT inhibitors (NNRTIs) in the treatment regimen but only 45% had resistance AZT/d4T associated resistance mutations (TAMs). TAMs were however detected at low frequency within the patients HIV quasispecies (1–20%) in 21 of 34 individuals who were failing first-line AZT-containing cART and lacked TAMs by Sanger. Due to lack of TAMs by Sanger, AZT was typically maintained in second-line therapies and these patients had a low frequency of subsequent virologic success.

Conclusions

Our findings suggest that continued use of AZT and d4T in first-line treatment in low-to-middle income countries may lead to misdiagnosis of HIV-1 drug resistance and possibly enhance a succession of second- and third-line treatment failures.

Electronic supplementary material

The online version of this article (doi: 10.1186/s40249-017-0377-0) contains supplementary material, which is available to authorized users.

Ultradeep sequencing of B and non-B HIV-1 subtypes: Viral diversity and drug resistance mutations before and after one month of antiretroviral therapy in naive patients

BACKGROUND

Ultradeep pyrosequencing technologies permit an assessment of the genetic diversity and the presence and frequency of minority variants in a viral population. The effect of these parameters on the outcome of highly active antiretroviral therapy (HAART) in HIV-infected patients is poorly understood.

OBJECTIVES

The present study used the pyrosequencing Roche 454 prototype assay to determine whether antiretroviral efficacy is correlated with viral diversity and minority drug resistance mutations in HIV-infected treatment-naive patients and to compare assay performance in B and non-B subtypes.

STUDY DESIGN

The study included 30 HIV-1 infected naive patients (20 with subtype non-B and 10 with subtype B). Ultradeep pyrosequencing of protease and reverse transcriptase genes was performed at baseline and 1 month after HAART initiation. Plasma HIV VL was measured at 0 and after 1, 3, and 6 months of HAART.

RESULTS

Pre-HAART minority drug resistance mutations were observed to NRTI in 4 patients, to NNRTI in 6 patients, and to PI in 1 patient; there was no difference in HAART-induced VL decay between patients. Pre-HAART diversity was significantly correlated with the time elapsed since HIV-1 infection diagnosis, but not with the subtype, VL, or CD4 count. Patients with an undetectable VL after 3 months of HAART had a higher pre-HAART diversity. Pre- and post-HAART diversities were not statistically different. There was no difference in assay performance between subtype B and non-B.

CONCLUSIONS

A high pre-HAART viral diversity might have a positive effect on the outcome of HAART. Pre-therapeutic minority drug resistance mutations are uncommon in naive patients.

HIV-1 Drug Resistance by Ultra-Deep Sequencing Following Short Course Zidovudine, Single-Dose Nevirapine, and Single-Dose Tenofovir with Emtricitabine for Prevention of Mother-to-Child Transmission

Supplemental Digital Content is Available in the Text.

Antiretroviral drug resistance following pMTCT strategies remains a significant problem. With rapid advancements in next generation sequencing technologies, there is more focus on HIV drug-resistant variants of low frequency, or the so-called minority variants. In South Africa, AZT monotherapy for pMTCT, similar to World Health Organization option A, has been used since 2008. In 2010, a single dose of co-formulated TDF/FTC was included in the strategy for prevention of resistance conferred by single-dose nevirapine (sd NVP). The study was conducted in KwaZulu-Natal, South Africa, among pMTCT participants who received AZT monotherapy from 14 weeks of gestation, intrapartum AZT and sd NVP, and postpartum sd TDF/FTC. Twenty-six specimens collected at 6 weeks post-delivery were successfully sequenced using 454 ultra-deep sequencing. Non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance was detected in 17 of 26 (65%) patients, 2 (7%) had Thymidine analogue mutations, and 3 (11%) had K65R. Of the 17 patients with NNRTI resistance, 11 (65%) had high-level NNRTI resistance, whereas 6 (35%) had intermediate NNRTI resistance. The levels of NNRTI resistance are much higher than would be expected, given the inclusion of antepartum AZT and postpartum TDF/FTC. This high level of NNRTI resistance could impact future NNRTI-containing treatment for a large proportion of pMTCT-exposed women. The detection of Thymidine analogue mutations highlights the need to understand the clinical impact of these on AZT-containing antiretroviral treatment in women exposed to AZT monotherapy.

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

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

Deep sequencing for HIV-1 clinical management

The emerging HIV-1 resistance epidemic is threatening the impressive global advances in HIV-1 infection treatment and prevention achieved in the last decade. Next-generation sequencing is improving our ability to understand, diagnose and prevent HIV-1 resistance, being increasingly cost-effective and more accessible. However, NGS still faces a number of limitations that need to be addressed to enable its widespread use. Here, we will review the main NGS platforms available for HIV-1 diagnosis, the factors affecting the clinical utility of NGS testing and the evidence supporting -or not- ultrasensitive genotyping over Sanger sequencing for routine HIV-1 diagnosis. Now that global HIV-1 eradication might be within our reach, making NGS accessible also to LMICs has become a priority. Reductions in sequencing costs, particularly in library preparation, and accessibility to low-cost, robust but simplified automated bioinformatic analyses of NGS data will remain essential to end the HIV-1 pandemic.

Pre-treatment minority HIV-1 drug resistance mutations and long term virological outcomes: is prediction possible?

Background

Although the use of highly active antiretroviral therapy in HIV positive individuals has proved to be effective in suppressing the virus to below detection limits of commonly used assays, virological failure associated with drug resistance is still a major challenge in some settings. The prevalence and effect of pre-treatment resistance associated variants on virological outcomes may also be underestimated because of reliance on conventional population sequencing data which excludes minority species. We investigated long term virological outcomes and the prevalence and pattern of pre-treatment minority drug resistance mutations in individuals initiating HAART at a local HIV clinic.

Methods

Patient’s records of viral load results and CD4 cell counts from routine treatment monitoring were used and additional pre-treatment blood samples for Sanger sequencing were obtained. A selection of pre-treatment samples from individuals who experienced virological failure were evaluated for minority resistance associated mutations to 1 % prevalence and compared to individuals who achieved viral suppression.

Results

At least one viral load result after 6 months or more of treatment was available for 65 out of 78 individuals followed for up to 33 months. Twenty (30.8 %) of the 65 individuals had detectable viremia and eight (12.3 %) of them had virological failure (viral load > 1000 RNA copies/ml) after at least 6 months of HAART. Viral suppression, achieved by month 8 to month 13, was followed by low level viremia in 10.8 % of patients and virological failure in one patient after month 20. There was potentially reduced activity to Emtricitabine or Tenofovir in three out of the eight cases in which minority drug resistance associated variants were investigated but detectable viremia occurred in one of these cases while the activity of Efavirenz was generally reduced in all the eight cases.

Conclusions

Early viral suppression was followed by low level viremia for some patients which may be an indication of failure to sustain viral suppression over time. The low level viremia may also be representing early stages of resistance development. The mutation patterns detected in the minority variants showed potential reduced drug sensitivity which highlights their potential to dominate after treatment initiation.

Trial registration

Not applicable.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0628-x) contains supplementary material, which is available to authorized users.

Recent advances in inferring viral diversity from high-throughput sequencing data

Rapidly evolving RNA viruses prevail within a host as a collection of closely related variants, referred to as viral quasispecies. Advances in high-throughput sequencing (HTS) technologies have facilitated the assessment of the genetic diversity of such virus populations at an unprecedented level of detail. However, analysis of HTS data from virus populations is challenging due to short, error-prone reads. In order to account for uncertainties originating from these limitations, several computational and statistical methods have been developed for studying the genetic heterogeneity of virus population. Here, we review methods for the analysis of HTS reads, including approaches to local diversity estimation and global haplotype reconstruction. Challenges posed by aligning reads, as well as the impact of reference biases on diversity estimates are also discussed. In addition, we address some of the experimental approaches designed to improve the biological signal-to-noise ratio. In the future, computational methods for the analysis of heterogeneous virus populations are likely to continue being complemented by technological developments.

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.

Human DDX3 protein is a valuable target to develop broad spectrum antiviral agents

Targeting a host factor essential for the replication of different viruses but not for the cells offers a higher genetic barrier to the development of resistance, may simplify therapy regimens for coinfections, and facilitates management of emerging viral diseases. DEAD-box polypeptide 3 (DDX3) is a human host factor required for the replication of several DNA and RNA viruses, including some of the most challenging human pathogens currently circulating, such as HIV-1, Hepatitis C virus, Dengue virus, and West Nile virus. Herein, we showed for the first time, to our knowledge, that the inhibition of DDX3 by a small molecule could be successfully exploited for the development of a broad spectrum antiviral agent. In addition to the multiple antiviral activities, hit compound 16d retained full activity against drug-resistant HIV-1 strains in the absence of cellular toxicity. Pharmacokinetics and toxicity studies in rats confirmed a good safety profile and bioavailability of 16d. Thus, DDX3 is here validated as a valuable therapeutic target.

Quantifying Next Generation Sequencing Sample Pre-Processing Bias in HIV-1 Complete Genome Sequencing

Genetic analyses play a central role in infectious disease research. Massively parallelized “mechanical cloning” and sequencing technologies were quickly adopted by HIV researchers in order to broaden the understanding of the clinical importance of minor drug-resistant variants. These efforts have, however, remained largely limited to small genomic regions. The growing need to monitor multiple genome regions for drug resistance testing, as well as the obvious benefit for studying evolutionary and epidemic processes makes complete genome sequencing an important goal in viral research. In addition, a major drawback for NGS applications to RNA viruses is the need for large quantities of input DNA. Here, we use a generic overlapping amplicon-based near full-genome amplification protocol to compare low-input enzymatic fragmentation (Nextera™) with conventional mechanical shearing for Roche 454 sequencing. We find that the fragmentation method has only a modest impact on the characterization of the population composition and that for reliable results, the variation introduced at all steps of the procedure—from nucleic acid extraction to sequencing—should be taken into account, a finding that is also relevant for NGS technologies that are now more commonly used. Furthermore, by applying our protocol to deep sequence a number of pre-therapy plasma and PBMC samples, we illustrate the potential benefits of a near complete genome sequencing approach in routine genotyping.

Towards Better Precision Medicine: PacBio Single-Molecule Long Reads Resolve the Interpretation of HIV Drug Resistant Mutation Profiles at Explicit Quasispecies (Haplotype) Level

Development of HIV-1 drug resistance mutations (HDRMs) is one of the major reasons for the clinical failure of antiretroviral therapy. Treatment success rates can be improved by applying personalized anti-HIV regimens based on a patient's HDRM profile. However, the sensitivity and specificity of the HDRM profile is limited by the methods used for detection. Sanger-based sequencing technology has traditionally been used for determining HDRM profiles at the single nucleotide variant (SNV) level, but with a sensitivity of only ≥ 20% in the HIV population of a patient. Next Generation Sequencing (NGS) technologies offer greater detection sensitivity (~ 1%) and larger scope (hundreds of samples per run). However, NGS technologies produce reads that are too short to enable the detection of the physical linkages of individual SNVs across the haplotype of each HIV strain present. In this article, we demonstrate that the single-molecule long reads generated using the Third Generation Sequencer (TGS), PacBio RS II, along with the appropriate bioinformatics analysis method, can resolve the HDRM profile at a more advanced quasispecies level. The case studies on patients' HIV samples showed that the quasispecies view produced using the PacBio method offered greater detection sensitivity and was more comprehensive for understanding HDRM situations, which is complement to both Sanger and NGS technologies. In conclusion, the PacBio method, providing a promising new quasispecies level of HDRM profiling, may effect an important change in the field of HIV drug resistance research.

Rapid identification of pathogens using molecular techniques

Real-time PCR is the traditional face of nucleic acid detection in the diagnostic microbiology laboratory and is now generally regarded as robust enough to be widely adopted. Methods based on nucleic acid detection of this type are bringing increased accuracy to diagnosis in areas where culture is difficult and/or expensive, and these methods are often effective partners to other rapid molecular diagnostic tools such as matrix-assisted laser desorption ionisation-time of flight mass spectrometry (MALDI-TOF MS). This change in practice has particularly affected the recognition of viruses and fastidious or antibiotic-exposed bacteria, but has been also shown to be effective in the recognition of troublesome or specialised phenotypes such as antiviral resistance and transmissible antibiotic resistance in the Enterobacteriaceae. Quantitation and high-intensity sequencing (of multiple whole genomes) has brought new opportunities as well as new challenges to the microbiology community. Diagnostic microbiologists currently training might be expected to deal less with the culture-based techniques of the last half-century than with the high-volume data and complex analyses of the next.

Transmission of HIV drug resistance: lessons from sensitive screening assays

PURPOSE OF REVIEW

The review discusses new technologies for the sensitive detection of HIV drug resistance, with a focus on applications in antiretroviral treatment (ART)-naïve populations.

RECENT FINDINGS

Conventional sequencing is well established for detecting HIV drug resistance in routine care and guides optimal treatment selection in patients starting ART. Access to conventional sequencing is nearly universal in Western countries, but remains limited in Asia, Latin America, and Africa. Technological advances now allow detection of resistance with greatly improved sensitivity compared with conventional sequencing, variably increasing the yield of resistance testing in ART-naïve populations. There is strong cumulative evidence from retrospective studies that sensitive detection of resistant mutants in baseline plasma samples lacking resistance by conventional sequencing more than doubles the risk of virological failure after starting efavirenz-based or nevirapine-based ART.

SUMMARY

Sensitive resistance testing methods are mainly confined to research applications and in this context have provided great insight into the dynamics of drug resistance development, persistence, and transmission. Adoption in care settings is becoming increasingly possible, although important challenges remain. Platforms for diagnostic use must undergo technical improvements to ensure good performance and ease of use, and clinical validation is required to ensure utility.

One-Step Ligation on RNA Amplification for the Detection of Point Mutations

The detection of point mutations is required in the diagnosis of many human diseases. The conformal specificity of DNA ligases was elegantly used to distinguish single-nucleotide mismatches. However, to detect point mutations in RNA retroviruses, conventional ligase-mediated approaches require the reverse transcription of viral genomes before separate ligation and amplification steps. We developed one-step ligation on RNA amplification (LRA) for the direct detection of RNA point mutations. The process combines the ligase-mediated joining of two oligonucleotides and subsequent hot start amplification into a single-tube reaction. We report that modifications to the structure of the oligonucleotide ligation probes improve the rate of ligation and the specificity of mutation detection on RNA. We applied LRA to the detection of a common, clinically relevant HIV-1 reverse transcriptase drug-resistant point mutation, K103N, and compared it with allele-specific PCR and pyrosequencing. LRA achieved a limit of specific quantitation of 1:100 (1%), and a limit of specific detection for mutant K103N RNA transcripts among excess wild-type strands of 1:10,000 (0.01%). LRA also exhibited good detection threshold of 5 × 10(2) copies/μL K103N RNA transcripts. LRA is a novel point mutation detection method, with potential utilization in HIV drug resistance detection and early diagnostics of genetic disorders associated with other infectious diseases and cancer.

Glycans flanking the hypervariable connecting peptide between the A and B strands of the V1/V2 domain of HIV-1 gp120 confer resistance to antibodies that neutralize CRF01_AE viruses

Understanding the molecular determinants of sensitivity and resistance to neutralizing antibodies is critical for the development of vaccines designed to prevent HIV infection. In this study, we used a genetic approach to characterize naturally occurring polymorphisms in the HIV envelope protein that conferred neutralization sensitivity or resistance. Libraries of closely related envelope genes, derived from virus quasi-species, were constructed from individuals infected with CRF01_AE viruses. The libraries were screened with plasma containing broadly neutralizing antibodies, and neutralization sensitive and resistant variants were selected for sequence analysis. In vitro mutagenesis allowed us to identify single amino acid changes in three individuals that conferred resistance to neutralization by these antibodies. All three mutations created N-linked glycosylation sites (two at N136 and one at N149) proximal to the hypervariable connecting peptide between the C-terminus of the A strand and the N-terminus of the B strand in the four-stranded V1/V2 domain β-sheet structure. Although N136 has previously been implicated in the binding of broadly neutralizing monoclonal antibodies, this glycosylation site appears to inhibit the binding of neutralizing antibodies in plasma from HIV-1 infected subjects. Previous studies have reported that the length of the V1/V2 domain in transmitted founder viruses is shorter and possesses fewer glycosylation sites compared to viruses isolated from chronic infections. Our results suggest that vaccine immunogens based on recombinant envelope proteins from clade CRF01_AE viruses might be improved by inclusion of envelope proteins that lack these glycosylation sites. This strategy might improve the efficacy of the vaccines used in the partially successful RV144 HIV vaccine trial, where the two CRF01_AE immunogens (derived from the A244 and TH023 isolates) both possessed glycosylation sites at N136 and N149.

Evaluation of GS Junior and MiSeq next-generation sequencing technologies as an alternative to Trugene population sequencing in the clinical HIV laboratory

Population HIV-1 sequencing is currently the method of choice for the identification and follow-up of HIV-1 antiretroviral drug resistance. It has limited sensitivity and results in a consensus sequence showing the most prevalent nucleotide per position. Moreover concomitant sequencing and interpretation of the results for several samples together is laborious and time consuming. In this study, the practical use of GS Junior and MiSeq bench-top next generation sequencing (NGS) platforms as an alternative to Trugene Sanger-based population sequencing in the clinical HIV laboratory was assessed. DeepChek(®)-HIV TherapyEdge software was used for processing all the protease and reverse transcriptase sequences and for resistance interpretation. Plasma samples from nine HIV-1 carriers, representing the major HIV-1 subtypes in Israel, were compared. The total number of amino acid substitutions identified in the nine samples by GS Junior (232 substitutions) and MiSeq (243 substitutions) was similar and higher than Trugene (181 substitutions), emphasizing the advantage of deep sequencing on population sequencing. More than 80% of the identified substitutions were identical between the GS Junior and MiSeq platforms, most of which (184 of 199) at similar frequency. Low abundance substitutions accounted for 20.9% of the MiSeq and 21.9% of the GS Junior output, the majority of which were not detected by Trugene. More drug resistance mutations were identified by both the NGS platforms, primarily, but not only, at low abundance. In conclusion, in combination with DeepChek, both GS Junior and MiSeq were found to be more sensitive than Trugene and adequate for HIV-1 resistance analysis in the clinical HIV laboratory.

Impact of human immunodeficiency virus type-1 sequence diversity on antiretroviral therapy outcomes

Worldwide circulating HIV-1 genomes show extensive variation represented by different subtypes, polymorphisms and drug-resistant strains. Reports on the impact of sequence variation on antiretroviral therapy (ART) outcomes are mixed. In this review, we summarize relevant published data from both resource-rich and resource-limited countries in the last 10 years on the impact of HIV-1 sequence diversity on treatment outcomes. The prevalence of transmission of drug resistant mutations (DRMs) varies considerably, ranging from 0% to 27% worldwide. Factors such as geographic location, access and availability to ART, duration since inception of treatment programs, quality of care, risk-taking behaviors, mode of transmission, and viral subtype all dictate the prevalence in a particular geographical region. Although HIV-1 subtype may not be a good predictor of treatment outcome, review of emerging evidence supports the fact that HIV-1 genome sequence-resulting from natural polymorphisms or drug-associated mutations-matters when it comes to treatment outcomes. Therefore, continued surveillance of drug resistant variants in both treatment-naïve and treatment-experienced populations is needed to reduce the transmission of DRMs and to optimize the efficacy of the current ART armamentarium.

VirVarSeq: a low-frequency virus variant detection pipeline for Illumina sequencing using adaptive base-calling accuracy filtering

MOTIVATION

In virology, massively parallel sequencing (MPS) opens many opportunities for studying viral quasi-species, e.g. in HIV-1- and HCV-infected patients. This is essential for understanding pathways to resistance, which can substantially improve treatment. Although MPS platforms allow in-depth characterization of sequence variation, their measurements still involve substantial technical noise. For Illumina sequencing, single base substitutions are the main error source and impede powerful assessment of low-frequency mutations. Fortunately, base calls are complemented with quality scores (Qs) that are useful for differentiating errors from the real low-frequency mutations.

RESULTS

A variant calling tool, Q-cpileup, is proposed, which exploits the Qs of nucleotides in a filtering strategy to increase specificity. The tool is imbedded in an open-source pipeline, VirVarSeq, which allows variant calling starting from fastq files. Using both plasmid mixtures and clinical samples, we show that Q-cpileup is able to reduce the number of false-positive findings. The filtering strategy is adaptive and provides an optimized threshold for individual samples in each sequencing run. Additionally, linkage information is kept between single-nucleotide polymorphisms as variants are called at the codon level. This enables virologists to have an immediate biological interpretation of the reported variants with respect to their antiviral drug responses. A comparison with existing SNP caller tools reveals that calling variants at the codon level with Q-cpileup results in an outstanding sensitivity while maintaining a good specificity for variants with frequencies down to 0.5%.

AVAILABILITY

The VirVarSeq is available, together with a user's guide and test data, at sourceforge: http://sourceforge.net/projects/virtools/?source=directory.

Contribution of human immunodeficiency virus type 1 minority variants to reduced drug susceptibility in patients on an integrase strand transfer inhibitor-based therapy

The role of HIV-1 minority variants on transmission, pathogenesis, and virologic failure to antiretroviral regimens has been explored; however, most studies of low-level HIV-1 drug-resistant variants have focused in single target regions. Here we used a novel HIV-1 genotypic assay based on deep sequencing, DEEPGEN (Gibson et al 2014 Antimicrob Agents Chemother 58∶2167) to simultaneously analyze the presence of minority variants carrying mutations associated with reduced susceptibility to protease (PR), reverse transcriptase (RT), and integrase strand transfer integrase inhibitors (INSTIs), as well as HIV-1 coreceptor tropism. gag-p2/NCp7/p1/p6/pol-PR/RT/INT and env/C2V3 PCR products were obtained from twelve heavily treatment-experienced patients experiencing virologic failure while participating in a 48-week dose-ranging study of elvitegravir (GS-US-183-0105). Deep sequencing results were compared with (i) virological response to treatment, (ii) genotyping based on population sequencing, (iii) phenotyping data using PhenoSense and VIRALARTS, and (iv) HIV-1 coreceptor tropism based on the phenotypic test VERITROP. Most patients failed the antiretroviral regimen with numerous pre-existing mutations in the PR and RT, and additionally newly acquired INSTI-resistance mutations as determined by population sequencing (mean 9.4, 5.3, and 1.4 PI- RTI-, and INSTI-resistance mutations, respectively). Interestingly, since DEEPGEN allows the accurate detection of amino acid substitutions at frequencies as low as 1% of the population, a series of additional drug resistance mutations were detected by deep sequencing (mean 2.5, 1.5, and 0.9, respectively). The presence of these low-abundance HIV-1 variants was associated with drug susceptibility, replicative fitness, and coreceptor tropism determined using sensitive phenotypic assays, enhancing the overall burden of resistance to all four antiretroviral drug classes. Further longitudinal studies based on deep sequencing tests will help to clarify (i) the potential impact of minority HIV-1 drug resistant variants in response to antiretroviral therapy and (ii) the importance of the detection of HIV minority variants in the clinical practice.

When can HIV clinical trials detect treatment effects on drug resistance?

Methods of sampling patients for resistance testing, and statistical analyses of HIV drug resistance, have not been standardised in HIV clinical trials. We analysed methods of genotyping and rates of treatment-emergent drug resistance from 27 clinical trials identified from a MEDLINE search. Sample size calculations were conducted using NQUERY software, assuming 5% significance level, 80% power and 1:1 randomisation. The percentage of patients with treatment-emergent IAS-USA mutations after 96 weeks ranged from 1.8% to 9.1% for first-line 2NRTI/NNRTI treatments, 0.6% to 6.3% for first-line 2NRTI/PI/r treatments and 0.0% to 2.0% in switch trials of boosted PIs. The prevalence of drug resistance was higher in trials with no screening for drug resistance at baseline, where the HIV RNA cut-off for genotyping was >50 copies/mL, where patients were tested for drug resistance after discontinuation of treatment, and where follow-up times were 96 weeks or longer. HIV clinical trials could be designed to detect differences in the risk of HIV drug resistance between treatments, as an analysis supporting HIV RNA suppression as the primary endpoint. However, this would require a standardised approach, with intent-to-treat analyses, testing of all samples with HIV RNA>50 copies/mL and genotyping after drug discontinuation.

An international multicenter study on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing

The detection of mutant spectra within the viral quasispecies is critical for therapeutic management of HIV-1 infections. Routine clinical application of ultrasensitive genotyping requires reproducibility and concordance within and between laboratories. The goal of the study was to evaluate a new protocol on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing (454-UDS) in an international multicenter study. Sixteen blinded HIV-1 subtype B samples were provided for 454-UDS as both RNA and cDNA with viral titers of 88,600-573,000 HIV-1 RNA copies/ml. Eight overlapping amplicons spanning protease (PR) codons 10-99 and reverse transcriptase (RT) codons 1-251 were generated using molecular barcoded primers. 454-UDS was performed using the 454 Life Sciences/Roche GS FLX platform. PR and RT sequences were analyzed using 454 Life Sciences Amplicon Variant Analyzer (AVA) software. Quantified variation data were analyzed for intra-laboratory reproducibility and inter-laboratory concordance. Routine population sequencing was performed using the ViroSeq HIV-1 genotyping system. Eleven laboratories and the reference laboratory 454 Life Sciences sequenced the HIV-1 sample set. Data presented are derived from seven laboratories and the reference laboratory since severe study protocol execution errors occurred in four laboratories leading to exclusion. The median sequencing depth across all sites was 1364 reads per position (IQR=809-2065). 100% of the ViroSeq-reported mutations were also detected by 454-UDS. Minority HIV-1 drug resistance mutations, defined as HIV-1 drug resistance mutations identified at frequencies of 1-25%, were only detected by 454-UDS. Analysis of 10 preselected majority and minority mutations were consistently found across sites. The analysis of drug-resistance mutations detected between 1 and 10% demonstrated high intra- and inter-laboratory consistency in frequency estimates for both RNA and prepared cDNA samples, indicating robustness of the method. HIV-1 drug resistance testing using 454 ultra-deep pyrosequencing results in an accurate and highly reproducible, albeit complex, approach to the analysis of HIV-1 mutant spectra, even at frequencies well below those detected by routine population sequencing.

Sensitive deep-sequencing-based HIV-1 genotyping assay to simultaneously determine susceptibility to protease, reverse transcriptase, integrase, and maturation inhibitors, as well as HIV-1 coreceptor tropism

With 29 individual antiretroviral drugs available from six classes that are approved for the treatment of HIV-1 infection, a combination of different phenotypic and genotypic tests is currently needed to monitor HIV-infected individuals. In this study, we developed a novel HIV-1 genotypic assay based on deep sequencing (DeepGen HIV) to simultaneously assess HIV-1 susceptibilities to all drugs targeting the three viral enzymes and to predict HIV-1 coreceptor tropism. Patient-derived gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3 PCR products were sequenced using the Ion Torrent Personal Genome Machine. Reads spanning the 3' end of the Gag, protease (PR), reverse transcriptase (RT), integrase (IN), and V3 regions were extracted, truncated, translated, and assembled for genotype and HIV-1 coreceptor tropism determination. DeepGen HIV consistently detected both minority drug-resistant viruses and non-R5 HIV-1 variants from clinical specimens with viral loads of ≥1,000 copies/ml and from B and non-B subtypes. Additional mutations associated with resistance to PR, RT, and IN inhibitors, previously undetected by standard (Sanger) population sequencing, were reliably identified at frequencies as low as 1%. DeepGen HIV results correlated with phenotypic (original Trofile, 92%; enhanced-sensitivity Trofile assay [ESTA], 80%; TROCAI, 81%; and VeriTrop, 80%) and genotypic (population sequencing/Geno2Pheno with a 10% false-positive rate [FPR], 84%) HIV-1 tropism test results. DeepGen HIV (83%) and Trofile (85%) showed similar concordances with the clinical response following an 8-day course of maraviroc monotherapy (MCT). In summary, this novel all-inclusive HIV-1 genotypic and coreceptor tropism assay, based on deep sequencing of the PR, RT, IN, and V3 regions, permits simultaneous multiplex detection of low-level drug-resistant and/or non-R5 viruses in up to 96 clinical samples. This comprehensive test, the first of its class, will be instrumental in the development of new antiretroviral drugs and, more importantly, will aid in the treatment and management of HIV-infected individuals.

Low-frequency nevirapine (NVP)-resistant HIV-1 variants are not associated with failure of antiretroviral therapy in women without prior exposure to single-dose NVP

BACKGROUND

Low-frequency nevirapine (NVP)-resistant variants have been associated with virologic failure (VF) of initial NVP-based combination antiretroviral therapy (cART) in women with prior exposure to single-dose NVP (sdNVP). We investigated whether a similar association exists in women without prior sdNVP exposure.

METHODS

Pre-cART plasma was analyzed by allele-specific polymerase chain reaction to quantify NVP-resistant mutants in human immunodeficiency virus-infected African women without prior sdNVP who were starting first-line NVP-based cART in the OCTANE/A5208 trial 2. Associations between NVP-resistant mutants and VF or death were determined and compared with published results from women participating in the OCTANE/A5208 trial 1 who had taken sdNVP and initiated NVP-based cART.

RESULTS

Pre-cART NVP-resistant variants were detected in 18% (39/219) of women without prior sdNVP exposure, compared to 45% (51/114) with prior sdNVP exposure (P < .001). Among women without prior sdNVP exposure, 8 of 39 (21%) with NVP-resistant variants experienced VF or death vs 31 of 180 (17%) without such variants (P = .65); this compares with 21 of 51 (41%) vs 9 of 63 (14%) among women with prior exposure (P = .001).

CONCLUSIONS

The risk of VF on NVP-based cART from NVP-resistant variants differs between sdNVP-exposed and -unexposed women. This difference may be driven by drug-resistance mutations emerging after sdNVP exposure that are linked on the same viral genome.

CLINICAL TRIALS REGISTRATION

NCT00089505.

Baseline clinical HIV genotypes are a valid measure of transmitted drug resistance within the treatment-naive population

OBJECTIVE

To examine whether baseline clinical genotypes are equivalent to diagnostic serum genotypes for surveillance of HIV transmitted drug resistance (TDR).

DESIGN

Current HIV TDR surveillance in Canada is conducted through genotyping remnant diagnostic sera from new HIV diagnoses. As part of routine care, baseline genotyping is now conducted on all newly diagnosed HIV infections, with TDR data being generated a second time on the same patients.

METHODS

Surveillance genotyping, on HIV diagnostic serum, was performed on newly diagnosed HIV cases from 2007 to 2010 in Alberta, Canada. All subjects with a baseline clinical genotype result on file, and no evidence of antiretroviral therapy, were studied further. The HIV sequences from diagnosis and from the first clinical genotype were compared according to elapsed time between testing and by evaluating timing of infection based on BED capture enzyme immunoassay (BED-CEIA, abbreviated as BED in this article).

RESULTS

Eighty-seven genotype pairs were available for analysis, most of which were subtype B. The time between genotypes ranged from 0 to 755 days, with a median of 36 days and an interquartile range of 155.25 days. Genetic distance between genotypes varied between 0 and 0.03389 substitutions per site and did not correlate with sampling times. There was a tendency for the genotypes of infections classified as recent by BED to be more similar to their clinical genotypes but this effect was lost when adjusted for elapsed time between tests. There was no difference in the identified drug resistance.

CONCLUSIONS

Baseline clinical genotypes from treatment-naive patients may be used for HIV TDR surveillance.

Molecular characterization of ambiguous mutations in HIV-1 polymerase gene: implications for monitoring HIV infection status and drug resistance

Detection of recent HIV infections is a prerequisite for reliable estimations of transmitted HIV drug resistance (t-HIVDR) and incidence. However, accurately identifying recent HIV infection is challenging due partially to the limitations of current serological tests. Ambiguous nucleotides are newly emerged mutations in quasispecies, and accumulate by time of viral infection. We utilized ambiguous mutations to establish a measurement for detecting recent HIV infection and monitoring early HIVDR development. Ambiguous nucleotides were extracted from HIV-1 pol-gene sequences in the datasets of recent (HIVDR threshold surveys [HIVDR-TS] in 7 countries; n=416) and established infections (1 HIVDR monitoring survey at baseline; n=271). An ambiguous mutation index of 2.04×10^-3 nts/site was detected in HIV-1 recent infections which is equivalent to the HIV-1 substitution rate (2×10^-3 nts/site/year) reported before. However, significantly higher index (14.41×10^-3 nts/site) was revealed with established infections. Using this substitution rate, 75.2% subjects in HIVDR-TS with the exception of the Vietnam dataset and 3.3% those in HIVDR-baseline were classified as recent infection within one year. We also calculated mutation scores at amino acid level at HIVDR sites based on ambiguous or fitted mutations. The overall mutation scores caused by ambiguous mutations increased (0.54×10^-23.48×10^-2/DR-site) whereas those caused by fitted mutations remained stable (7.50-7.89×10^-2/DR-site) in both recent and established infections, indicating that t-HIVDR exists in drug-naïve populations regardless of infection status in which new HIVDR continues to emerge. Our findings suggest that characterization of ambiguous mutations in HIV may serve as an additional tool to differentiate recent from established infections and to monitor HIVDR emergence.

Sensitive assessment of the virologic outcomes of stopping and restarting non-nucleoside reverse transcriptase inhibitor-based antiretroviral therapy

BACKGROUND

Non-nucleoside reverse transcriptase inhibitor (NNRTI)-resistant mutants have been shown to emerge after interruption of suppressive NNRTI-based antiretroviral therapy (ART) using routine testing. The aim of this study was to quantify the risk of resistance by sensitive testing and correlate the detection of resistance with NNRTI concentrations after treatment interruption and virologic responses after treatment resumption.

METHODS

Resistance-associated mutations (RAMs) and NNRTI concentrations were studied in plasma from 132 patients who interrupted suppressive ART within SMART. RAMs were detected by Sanger sequencing, allele-specific PCR, and ultra-deep sequencing. NNRTI concentrations were measured by sensitive high-performance liquid chromatography.

RESULTS

Four weeks after NNRTI interruption, 19/31 (61.3%) and 34/39 (87.2%) patients showed measurable nevirapine (>0.25 ng/ml) or efavirenz (>5 ng/ml) concentrations, respectively. Median eight weeks after interruption, 22/131 (16.8%) patients showed ≥1 NNRTI-RAM, including eight patients with NNRTI-RAMs detected only by sensitive testing. The adjusted odds ratio (OR) of NNRTI-RAM detection was 7.62 (95% confidence interval [CI] 1.52, 38.30; p = 0.01) with nevirapine or efavirenz concentrations above vs. below the median measured in the study population. Staggered interruption, whereby nucleos(t)ide reverse transcriptase inhibitors (NRTIs) were continued for median nine days after NNRTI interruption, did not prevent NNRTI-RAMs, but increased detection of NRTI-RAMs (OR 4.25; 95% CI 1.02, 17.77; p = 0.03). After restarting NNRTI-based ART (n = 90), virologic suppression rates <400 copies/ml were 8/13 (61.5%) with NNRTI-RAMs, 7/11 (63.6%) with NRTI-RAMs only, and 51/59 (86.4%) without RAMs. The ORs of re-suppression were 0.18 (95% CI 0.03, 0.89) and 0.17 (95% CI 0.03, 1.15) for patients with NNRTI-RAMs or NRTI-RAMs only respectively vs. those without RAMs (p = 0.04).

CONCLUSIONS

Detection of resistant mutants in the rebound viremia after interruption of efavirenz- or nevirapine-based ART affects outcomes once these drugs are restarted. Further studies are needed to determine RAM persistence in untreated patients and impact on newer NNRTIs.

Origin of minority drug-resistant HIV-1 variants in primary HIV-1 infection

BACKGROUND

Drug-resistant human immunodeficiency virus type 1 (HIV-1) minority variants (MVs) are present in some antiretroviral therapy (ART)-naive patients. They may result from de novo mutagenesis or transmission. To date, the latter has not been proven.

METHODS

MVs were quantified by allele-specific polymerase chain reaction in 204 acute or recent seroconverters from the Zurich Primary HIV Infection study and 382 ART-naive, chronically infected patients. Phylogenetic analyses identified transmission clusters.

RESULTS

Three lines of evidence were observed in support of transmission of MVs. First, potential transmitters were identified for 12 of 16 acute or recent seroconverters harboring M184V MVs. These variants were also detected in plasma and/or peripheral blood mononuclear cells at the estimated time of transmission in 3 of 4 potential transmitters who experienced virological failure accompanied by the selection of the M184V mutation before transmission. Second, prevalence between MVs harboring the frequent mutation M184V and the particularly uncommon integrase mutation N155H differed highly significantly in acute or recent seroconverters (8.2% vs 0.5%; P < .001). Third, the prevalence of less-fit M184V MVs is significantly higher in acutely or recently than in chronically HIV-1-infected patients (8.2% vs 2.5%; P = .004).

CONCLUSIONS

Drug-resistant HIV-1 MVs can be transmitted. To what extent the origin-transmission vs sporadic appearance-of these variants determines their impact on ART needs to be further explored.

HIV Drug Resistance: Problems and Perspectives

Access to combination antiretroviral treatment (ART) has improved greatly over recent years. At the end of 2011, more than eight million HIV-infected people were receiving ART in low-income and middle-income countries. ART generally works well in keeping the virus suppressed and the patient healthy. However, treatment only works as long as the virus is not resistant against the drugs used. In the last decades, HIV treatments have become better and better at slowing down the evolution of drug resistance, so that some patients are treated for many years without having any resistance problems. However, for some patients, especially in low-income countries, drug resistance is still a serious threat to their health. This essay will review what is known about transmitted and acquired drug resistance, multi-class drug resistance, resistance to newer drugs, resistance due to treatment for the prevention of mother-to-child transmission, the role of minority variants (low-frequency drug-resistance mutations), and resistance due to pre-exposure prophylaxis.

Detection of HIV-1 minority variants containing the K103N drug-resistance mutation using a simple method to amplify RNA targets (SMART)

The simple method for amplifying RNA targets (SMART) was used to detect K103N, a common HIV-1 reverse transcriptase drug-resistance mutation. Novel amplifiable SMART probes served as reporter molecules for RNA sequences that are captured and separated on a microfluidic platform under zero-flow conditions. Assays were performed both off chip and in a microchip reservoir using a modified version of real-time nucleic acid sequence-based amplification, without the noncyclic phase, and 65°C preheat. A total of 6000 copies/mL of the synthetic sequences were detected within 180 minutes of amplification. Although the sensitivity of research platforms is higher, SMART has the potential to offer comparable sensitivity and speed to commercially available viral load and HIV detection kits. Furthermore, SMART uses an inexpensive, practical, and more accurate isothermal exponential amplification technique. The use of molecular beacons resulted in relatively fast real-time detection (<180 minutes); however, they were also shown to hinder the amplification process when compared with end point detection. Finally, SMART probes were used for modeling of K103N concentrations within an unknown sample. Only 1% of the SMART probes was detected within the wild-type population (6 × 10(8) copies/mL). These results establish the groundwork for point-of-care drug resistance and viral load monitoring in clinical samples, which can revolutionize HIV patient care globally.

How bioinformatics influences health informatics: usage of biomolecular sequences, expression profiles and automated microscopic image analyses for clinical needs and public health

ABSTRACT

The currently hyped expectation of personalized medicine is often associated with just achieving the information technology led integration of biomolecular sequencing, expression and histopathological bioimaging data with clinical records at the individual patients' level as if the significant biomedical conclusions would be its more or less mandatory result. It remains a sad fact that many, if not most biomolecular mechanisms that translate the human genomic information into phenotypes are not known and, thus, most of the molecular and cellular data cannot be interpreted in terms of biomedically relevant conclusions. Whereas the historical trend will certainly be into the general direction of personalized diagnostics and cures, the temperate view suggests that biomedical applications that rely either on the comparison of biomolecular sequences and/or on the already known biomolecular mechanisms have much greater chances to enter clinical practice soon. In addition to considering the general trends, we exemplarily review advances in the area of cancer biomarker discovery, in the clinically relevant characterization of patient-specific viral and bacterial pathogens (with emphasis on drug selection for influenza and enterohemorrhagic E. coli) as well as progress in the automated assessment of histopathological images. As molecular and cellular data analysis will become instrumental for achieving desirable clinical outcomes, the role of bioinformatics and computational biology approaches will dramatically grow.

AUTHOR SUMMARY

With DNA sequencing and computers becoming increasingly cheap and accessible to the layman, the idea of integrating biomolecular and clinical patient data seems to become a realistic, short-term option that will lead to patient-specific diagnostics and treatment design for many diseases such as cancer, metabolic disorders, inherited conditions, etc. These hyped expectations will fail since many, if not most biomolecular mechanisms that translate the human genomic information into phenotypes are not known yet and, thus, most of the molecular and cellular data collected will not lead to biomedically relevant conclusions. At the same time, less spectacular biomedical applications based on biomolecular sequence comparison and/or known biomolecular mechanisms have the potential to unfold enormous potential for healthcare and public health. Since the analysis of heterogeneous biomolecular data in context with clinical data will be increasingly critical, the role of bioinformatics and computational biology will grow correspondingly in this process.