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

HIV-1 genotypic resistance testing using single molecule real-time sequencing

BACKGROUND

HIV-1 resistance testing is recommended in clinical management and next-generation sequencing (NGS) methods are now available in many virology laboratories.

OBJECTIVES

To evaluate the diagnostic performance of Long-Read Single Molecule Real-time (SMRT) sequencing (Sequel, PacBio) for HIV-1 polymerase genotyping.

STUDY DESIGN

111 prospective clinical samples (83 plasma and 28 leukocyte-enriched blood fraction) were analyzed for routine HIV-1 resistance genotyping using Sanger sequencing, Vela NGS, and SMRT sequencing. We developed a SMRT sequencing protocol and a bio-informatics pipeline to infer antiretroviral resistance on both haplotype and variant calling approaches.

RESULTS

The polymerase was successfully sequenced by the three platforms in 98 % of plasma RNA samples for viral loads above 4 log copies/mL. The success rate decreased to 83 % using Sanger or Vela sequencing and to 67 % using SMRT sequencing for viral loads of 3 to 4 log copies/mL. Sensitivities of 50 %, 54 % and 61 % were obtained using SMRT, Vela, and Sanger sequencing, respectively, in cellular DNA from patients with prolonged undetectable plasma HIV-1 RNA. Ninety-eight percent of resistance-associated mutations (RAMs) identified with Sanger sequencing were detected using SMRT sequencing. Furthermore, 91 % of RAMs (> 5 % threshold) identified with Vela NGS were detected using SMRT sequencing. RAM quantification using Vela and SMRT sequencing was well correlated (Spearman correlation ρ = 0.82; P < 0.0001).

CONCLUSIONS

SMRT sequencing of the full-length HIV-1 polymerase appeared performant for characterizing HIV-1 genotypic resistance on both RNA and DNA clinical samples. Long-read sequencing is a new tool for mutation haplotyping and resistance analysis.

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.

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.

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.

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.

High frequency of multiclass HCV resistance-associated mutations in patients failing direct-acting antivirals: real-life data

BACKGROUND

Direct-acting antiviral (DAA) therapy has dramatically increased sustained virological response rates in HCV-infected patients. However, resistance-associated substitutions (RAS) interfering with NS3- and NS5A-targeted therapy, still emerge. This real-life study analysed the type and frequency of RAS in rare cases of patients failing DAA regimens in 12 clinical centres in Israel.

METHODS

Blood samples and clinical data from 49 patients who failed various DAAs were collected. RAS identified in the NS3 and NS5A regions by population (Sanger) and next-generation sequencing (NGS) were compared by treatment regimen and HCV subtypes.

RESULTS

The majority (71.4%, 35/49) of patients were infected with the genotype (GT)1b strain, while 12.2% (6/49) carried GT1a and 14.3% GT3a/b (7), GT4a (1) and GT1b/GT3a. RAS were identified in 85.7% (42/49) of failures, of which 90.5% (38/42) were clinically relevant RAS (known to be associated with a specific GT and DAA in patients failing therapy or those with more than twofold change in in vitro replicon assays). The most abundant RAS were 168A/E/Q/G/N/V (32.6%, 16/49) identified in NS3, and 93H/N (61.2%, 30/49), 31I/M/V (34.7%, 17/49) and 30R/H/K (12.2%, 6/49), identified in NS5A. Significantly more clinically relevant RAS were identified in NS5A (82.2%, 37/45) than in NS3 (35.7%, 10/28; P<0.01). While RAS were identified in all GT1a, GT3b and GT4a failures (100%, 10/10), only 71.8% (28/39) of GT1b or GT3a failures had RAS (P=0.09). In four cases, NGS identified additional clinically relevant RAS and in one patient, NGS deciphered coexistence of GT3a and GT1b infections.

CONCLUSIONS

Our findings, together with additional real-life data, will contribute to the optimization of retreatment in DAA failure, when cost-related and suboptimal regimens must be employed.

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.

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.

The dawn of precision medicine in HIV: state of the art of pharmacotherapy

INTRODUCTION

Combination antiretroviral therapy (ART) reduces viral load to under the limit of detection, successfully decreasing HIV-related morbidity and mortality. Due to viral mutations, complex drug combinations and different patient response, there is an increasing demand for individualized treatment options for patients.

AREAS COVERED

This review first summarizes the pharmacokinetic and pharmacodynamic profile of clinical first-line drugs, which serves as guidance for antiretroviral precision medicine. Factors which have influential effects on drug efficacy and thus precision medicine are discussed: patients' pharmacogenetic information, virus mutations, comorbidities, and immune recovery. Furthermore, strategies to improve the application of precision medicine are discussed.

EXPERT OPINION

Precision medicine for ART requires comprehensive information on the drug, virus, and clinical data from the patients. The clinically available genetic tests are a good starting point. To better apply precision medicine, deeper knowledge of drug concentrations, HIV reservoirs, and efficacy associated genes, such as polymorphisms of drug transporters and metabolizing enzymes, are required. With advanced computer-based prediction systems which integrate more comprehensive information on pharmacokinetics, pharmacodynamics, pharmacogenomics, and the clinically relevant information of the patients, precision medicine will lead to better treatment choices and improved disease outcomes.

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.

Comparison between next-generation and Sanger-based sequencing for the detection of transmitted drug-resistance mutations among recently infected HIV-1 patients in Israel, 2000-2014

INTRODUCTION

Transmitted drug-resistance mutations (TDRM) may hamper successful anti-HIV-1 therapy and impact future control of the HIV-1 epidemic. Recently infected, therapy-naïve individuals are best suited for surveillance of such TDRM. In this study, TDRM, detected by next-generation sequencing (NGS) were compared to those identified by Sanger-based population sequencing (SBS) in recently infected HIV-1 patients.

METHODS

Historical samples from 80 recently infected HIV-1 patients, diagnosed between 2000 and 2014, were analysed by MiSeq (NGS) and ABI (SBS). DeepChek-HIV (ABL) was used for interpretation of the results.

RESULTS

Most patients were males (80%); Men who have sex with men (MSM) was the major transmission group (58.8%). Overall, TDRM were detected in 31.3% of patients by NGS and 8.8% by SBS, with SBS TDRM restricted to persons infected with subtype B. All SBS-detected TDRM were identified by NGS. The prevalence of TDRM impacting protease inhibitors (PI), nucleoside reverse transcriptase inhibitors (NRTI) and non-nucleoside reverse transcriptase inhibitors (NNRTI) was 11.3, 26.2 7.5%, respectively, in NGS analyses and 0, 3.8 and 5%, respectively, in SBS analyses. More patients with NGS and SBS TDRM were identified in 2008-2014 (37.2% or 13.9%, respectively) compared to 2000-2007 (24.3% or 2.7%, respectively), and a significantly greater number of these patients had multiple NGS TDRM. The most abundant, albeit, minor-frequency RT TDRM, were the K65R and D67N, while K103N, M184V and T215S were high-frequency mutations. Minor TDRM did not become a major variant in later samples and did not hinder successful treatment.

CONCLUSIONS

NGS can replace SBS for mutation detection and allows for the detection of low-frequency TDRM not identified by SBS. Although rates of TDRM in Israel continued to increase from 2000 to 2014, minor TDRM did not become major species. The need for ongoing surveillance of low-frequency TDRM should be revisited in a larger study.

Performance comparison of next-generation sequencing platforms for determining HIV-1 coreceptor use

The coreceptor used by HIV-1 must be determined before a CCR5 antagonist, part of the arsenal of antiretroviral drugs, is prescribed because viruses that enter cells using the CXCR4 coreceptor are responsible for treatment failure. HIV-1 tropism is also correlated with disease progression and so must be determined for virological studies. Tropism can be determined by next-generation sequencing (NGS), but not all of these new technologies have been fully validated for use in clinical practice. The Illumina NGS technology is used in many laboratories but its ability to predict HIV-1 tropism has not been evaluated while the 454 GS-Junior (Roche) is used for routine diagnosis. The genotypic prediction of HIV-1 tropism is based on sequencing the V3 region and interpreting the results with an appropriate algorithm. We compared the performances of the MiSeq (Illumina) and 454 GS-Junior (Roche) systems with a reference phenotypic assay. We used clinical samples for the NGS tropism predictions and assessed their ability to quantify CXCR4-using variants. The data show that the Illumina platform can be used to detect minor CXCR4-using variants in clinical practice but technical optimization are needed to improve quantification.

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.

Linking short tandem repeat polymorphisms with cytosine modifications in human lymphoblastoid cell lines

Inter-individual variation in cytosine modifications has been linked to complex traits in humans. Cytosine modification variation is partially controlled by single nucleotide polymorphisms (SNPs), known as modified cytosine quantitative trait loci (mQTL). However, little is known about the role of short tandem repeat polymorphisms (STRPs), a class of structural genetic variants, in regulating cytosine modifications. Utilizing the published data on the International HapMap Project lymphoblastoid cell lines (LCLs), we assessed the relationships between 721 STRPs and the modification levels of 283,540 autosomal CpG sites. Our findings suggest that, in contrast to the predominant cis-acting mode for SNP-based mQTL, STRPs are associated with cytosine modification levels in both cis-acting (local) and trans-acting (distant) modes. In local scans within the ±1 Mb windows of target CpGs, 21, 9, and 21 cis-acting STRP-based mQTL were detected in CEU (Caucasian residents from Utah, USA), YRI (Yoruba people from Ibadan, Nigeria), and the combined samples, respectively. In contrast, 139,420, 76,817, and 121,866 trans-acting STRP-based mQTL were identified in CEU, YRI, and the combined samples, respectively. A substantial proportion of CpG sites detected with local STRP-based mQTL were not associated with SNP-based mQTL, suggesting that STRPs represent an independent class of mQTL. Functionally, genetic variants neighboring CpG-associated STRPs are enriched with genome-wide association study (GWAS) loci for a variety of complex traits and diseases, including cancers, based on the National Human Genome Research Institute (NHGRI) GWAS Catalog. Therefore, elucidating these STRP-based mQTL in addition to SNP-based mQTL can provide novel insights into the genetic architectures of complex traits.

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.

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

Limited access to HIV drug resistance testing in low- and middle-income countries impedes clinical decision-making at the individual patient level. An efficient protocol to address this issue must be established to minimize negative therapeutic outcomes for HIV-1-infected individuals in such settings. This is an observational study to ascertain the potential of newer genomic sequencing platforms, such as the Illumina MiSeq instrument, to provide accurate HIV drug resistance genotypes for hundreds of samples simultaneously. Plasma samples were collected from Canadian patients during routine drug resistance testing (n = 759) and from a Ugandan study cohort (n = 349). Amplicons spanning HIV reverse transcriptase codons 90 to 234 were sequenced with both MiSeq sequencing and conventional Sanger sequencing methods. Sequences were evaluated for nucleotide concordance between methods, using coverage and mixture parameters for quality control. Consensus sequences were also analyzed for disparities in the identification of drug resistance mutations. Sanger and MiSeq sequencing was successful for 881 samples (80%) and 892 samples (81%), respectively, with 832 samples having results from both methods. Most failures were for samples with viral loads of <3.0 log10 HIV RNA copies/ml. Overall, 99.3% nucleotide concordance between methods was observed. MiSeq sequencing achieved 97.4% sensitivity and 99.3% specificity in detecting resistance mutations identified by Sanger sequencing. Findings suggest that the Illumina MiSeq platform can yield high-quality data with a high-multiplex "wide" sequencing approach. This strategy can be used for multiple HIV subtypes, demonstrating the potential for widespread individual testing and annual population surveillance in resource-limited settings.

How next-generation sequencing and multiscale data analysis will transform infectious disease management

We should integrate next-generation sequencing data from pathogen specimens with phenotypes from electronic medical records to create quantitative, predictive models of infectious disease. Precision infection control and antimicrobial interventions can address urgent global problems, including healthcare-associated infections and multidrug resistance.

Recent reviews have examined the extent to which routine next-generation sequencing (NGS) on clinical specimens will improve the capabilities of clinical microbiology laboratories in the short term, but do not explore integrating NGS with clinical data from electronic medical records (EMRs), immune profiling data, and other rich datasets to create multiscale predictive models. This review introduces a range of “omics” and patient data sources relevant to managing infections and proposes 3 potentially disruptive applications for these data in the clinical workflow. The combined threats of healthcare-associated infections and multidrug-resistant organisms may be addressed by multiscale analysis of NGS and EMR data that is ideally updated and refined over time within each healthcare organization. Such data and analysis should form the cornerstone of future learning health systems for infectious disease.