Use of four next-generation sequencing platforms to determine HIV-1 coreceptor tropism

Evolution of Multiple Domains of the HIV-1 Envelope Glycoprotein during Coreceptor Switch with CCR5 Antagonist Therapy

HIV-1 uses CD4 as a receptor and chemokine receptors CCR5 and/or CXCR4 as coreceptors. CCR5 antagonists are a class of antiretrovirals used to inhibit viral entry. Phenotypic prediction algorithms such as Geno2Pheno are used to assess CCR5 antagonist eligibility, for which the V3 region is screened. However, there exist scenarios where the algorithm cannot give an accurate prediction of tropism. The current study examined coreceptor shift of HIV-1 from CCR5-tropic strains to CXCR4-tropic or dual-tropic strains among five subjects in a clinical trial of the CCR5 antagonist vicriviroc. Envelope gene amplicon libraries were constructed and subjected to next-generation sequencing, as well as single-clone sequencing and functional analyses. Approximately half of the amplified full-length single envelope-encoding clones had no significant activity for infection of cells expressing high levels of CD4 and CCR5 or CXCR4. Functional analysis of 9 to 21 individual infectious clones at baseline and at the time of VF were used to construct phylogenetic trees and sequence alignments. These studies confirmed that specific residues and the overall charge of the V3 loop were the major determinants of coreceptor use, in addition to specific residues in other domains of the envelope protein in V1/V2, V4, C3, and C4 domains that may be important for coreceptor shift. These results provide greater insight into the viral genetic determinants of coreceptor shift. IMPORTANCE This study is novel in combining single-genome sequence analysis and next-generation sequencing to characterize HIV-1 quasispecies. The work highlights the importance of mutants present at frequencies of 1% or less in development of drug resistance. This study highlights a critical role of specific amino acid substitutions outside V3 that contribute to coreceptor shift as well as important roles of the V1/V2, V4, C3, and C4 domain residues.

Reduced and highly diverse peripheral HIV-1 reservoir in virally suppressed patients infected with non-B HIV-1 strains in Uganda

BACKGROUND

Our understanding of the peripheral human immunodeficiency virus type 1 (HIV-1) reservoir is strongly biased towards subtype B HIV-1 strains, with only limited information available from patients infected with non-B HIV-1 subtypes, which are the predominant viruses seen in low- and middle-income countries (LMIC) in Africa and Asia.

RESULTS

In this study, blood samples were obtained from well-suppressed ART-experienced HIV-1 patients monitored in Uganda (n = 62) or the U.S. (n = 50), with plasma HIV-1 loads < 50 copies/ml and CD4+ T-cell counts > 300 cells/ml. The peripheral HIV-1 reservoir, i.e., cell-associated HIV-1 RNA and proviral DNA, was characterized using our novel deep sequencing-based EDITS assay. Ugandan patients were slightly younger (median age 43 vs 49 years) and had slightly lower CD4+ counts (508 vs 772 cells/ml) than U.S. individuals. All Ugandan patients were infected with non-B HIV-1 subtypes (31% A1, 64% D, or 5% C), while all U.S. individuals were infected with subtype B viruses. Unexpectedly, we observed a significantly larger peripheral inducible HIV-1 reservoir in U.S. patients compared to Ugandan individuals (48 vs. 11 cell equivalents/million cells, p < 0.0001). This divergence in reservoir size was verified measuring proviral DNA (206 vs. 88 cell equivalents/million cells, p < 0.0001). However, the peripheral HIV-1 reservoir was more diverse in Ugandan than in U.S. individuals (8.6 vs. 4.7 p-distance, p < 0.0001).

CONCLUSIONS

The smaller, but more diverse, peripheral HIV-1 reservoir in Ugandan patients might be associated with viral (e.g., non-B subtype with higher cytopathicity) and/or host (e.g., higher incidence of co-infections or co-morbidities leading to less clonal expansion) factors. This highlights the need to understand reservoir dynamics in diverse populations as part of ongoing efforts to find a functional cure for HIV-1 infection in LMICs.

Application of Next Generation Sequencing in Laboratory Medicine

The rapid development of next-generation sequencing (NGS) technology, including advances in sequencing chemistry, sequencing technologies, bioinformatics, and data interpretation, has facilitated its wide clinical application in precision medicine. This review describes current sequencing technologies, including short- and long-read sequencing technologies, and highlights the clinical application of NGS in inherited diseases, oncology, and infectious diseases. We review NGS approaches and clinical diagnosis for constitutional disorders; summarize the application of U.S. Food and Drug Administration-approved NGS panels, cancer biomarkers, minimal residual disease, and liquid biopsy in clinical oncology; and consider epidemiological surveillance, identification of pathogens, and the importance of host microbiome in infectious diseases. Finally, we discuss the challenges and future perspectives of clinical NGS tests.

Molecular epidemiology in the HIV and SARS-CoV-2 pandemics

PURPOSE OF REVIEW

The aim of this review was to compare and contrast the application of molecular epidemiology approaches for the improved management and understanding of the HIV versus SARS-CoV-2 epidemics.

RECENT FINDINGS

Molecular biology approaches, including PCR and whole genome sequencing (WGS), have become powerful tools for epidemiological investigation. PCR approaches form the basis for many high-sensitivity diagnostic tests and can supplement traditional contact tracing and surveillance strategies to define risk networks and transmission patterns. WGS approaches can further define the causative agents of disease, trace the origins of the pathogen, and clarify routes of transmission. When coupled with clinical datasets, such as electronic medical record data, these approaches can investigate co-correlates of disease and pathogenesis. In the ongoing HIV epidemic, these approaches have been effectively deployed to identify treatment gaps, transmission clusters and risk factors, though significant barriers to rapid or real-time implementation remain critical to overcome. Likewise, these approaches have been successful in addressing some questions of SARS-CoV-2 transmission and pathogenesis, but the nature and rapid spread of the virus have posed additional challenges.

SUMMARY

Overall, molecular epidemiology approaches offer unique advantages and challenges that complement traditional epidemiological tools for the improved understanding and management of epidemics.

Compartmentalization and Clonal Amplification of HIV-1 in the Male Genital Tract Characterized Using Next-Generation Sequencing

Compartmentalization of HIV-1 between the systemic circulation and the male genital tract may have a substantial impact on which viruses are available for sexual transmission to new hosts. We studied compartmentalization and clonal amplification of HIV-1 populations between the blood and the genital tract from 10 antiretroviral-naive men using Illumina MiSeq with a PrimerID approach. We found evidence of some degree of compartmentalization in every study participant, unlike previous studies, which collectively showed that only ∼50% of analyzed individuals exhibited compartmentalization of HIV-1 lineages between the male genital tract (MGT) and blood. Using down-sampling simulations, we determined that this disparity can be explained by differences in sampling depth in that had we sequenced to a lower depth, we would also have found compartmentalization in only ∼50% of the study participants. For most study participants, phylogenetic trees were rooted in blood, suggesting that the male genital tract reservoir is seeded by incoming variants from the blood. Clonal amplification was observed in all study participants and was a characteristic of both blood and semen viral populations. We also show evidence for independent viral replication in the genital tract in the individual with the most severely compartmentalized HIV-1 populations. The degree of clonal amplification was not obviously associated with the extent of compartmentalization. We were also unable to detect any association between history of sexually transmitted infections and level of HIV-1 compartmentalization. Overall, our findings contribute to a better understanding of the dynamics that affect the composition of virus populations that are available for transmission.IMPORTANCE Within an individual living with HIV-1, factors that restrict the movement of HIV-1 between different compartments-such as between the blood and the male genital tract-could strongly influence which viruses reach sites in the body from which they can be transmitted. Using deep sequencing, we found strong evidence of restricted HIV-1 movements between the blood and genital tract in all 10 men that we studied. We additionally found that neither the degree to which particular genetic variants of HIV-1 proliferate (in blood or genital tract) nor an individual's history of sexually transmitted infections detectably influenced the degree to which virus movements were restricted between the blood and genital tract. Last, we show evidence that viral replication gave rise to a large clonal amplification in semen in a donor with highly compartmentalized HIV-1 populations, raising the possibility that differential selection of HIV-1 variants in the genital tract may occur.

Single-molecular real-time deep sequencing reveals the dynamics of multi-drug resistant haplotypes and structural variations in the hepatitis C virus genome

While direct-acting antivirals (DAAs) for hepatitis C virus (HCV) have dramatically progressed, patients still suffer from treatment failures. For the radical eradication of HCV, a deeper understanding of multiple resistance-associated substitutions (RASs) at the single-clone level is essential. To understand HCV quasispecies and their dynamics during DAA treatment, we applied single-molecule real-time (SMRT) deep sequencing on sera from 12 patients with genotype-1b HCV infections with DAA treatment failures, both pre- and post-treatment. We identified >3.2 kbp sequences between NS3 and NS5A genes of 187,539 clones in total, classifying into haplotype codes based on the linkage of seven RAS loci. The number of haplotype codes during the treatment, per sample, significantly decreased from 14.67 ± 9.12 to 6.58 ± 7.1, while the number of nonsynonymous codons on the seven RAS loci, per clone, significantly increased from 1.50 ± 0.92 to 3.64 ± 0.75. In five cases, the minority multi-drug resistant haplotypes at pre-treatment were identical to the major haplotypes at relapse. Moreover, various structural variations (SVs) were detected and their dynamics analysed. These results suggest that SMRT deep sequencing is useful for detecting minority haplotypes and SVs, and to evaluate the dynamics of viral genomes at the single-clone level.

Companion and complementary diagnostics for infectious diseases

INTRODUCTION

Companion diagnostics (CDx) are important in oncology therapeutic decision-making, but specific regulatory-approved CDx for infectious disease treatment are officially lacking. While not approved as CDx, several ID diagnostics are used as CDx. The diagnostics community, manufacturers, and regulatory agencies have made major efforts to ensure that diagnostics for new antimicrobials are available at or near release of new agents.

AREAS COVERED

This review highlights the status of Complementary and companion diagnostic (c/CDx) in the infectious disease literature, with a focus on genotypic antimicrobial resistance testing against pathogens as a class of diagnostic tests.

EXPERT OPINION

CRISPR, sepsis markers, and narrow spectrum antimicrobials, in addition to current and emerging technologies, present opportunities for infectious disease c/CDx. Challenges include slow guideline revision, high costs for regulatory approval, lengthy buy in by agencies, discordant pharmaceutical/diagnostic partnerships, and higher treatment costs. The number of patients and available medications used to treat different infectious diseases is well suited to support competing diagnostic tests. However, newer approaches to treatment (for example, narrow spectrum antibiotics), may be well suited for a small number of patients, i.e. a niche market in support of a CDx. The current emphasis is rapid and point-of-care (POC) diagnostic platforms as well as changes in treatment.

Long-term evolution of transmitted CXCR4-using HIV-1 under effective antiretroviral therapy

OBJECTIVE

To study the long-term evolution of the transmitted CXCR4-using viruses. CCR5-using viruses (R5 viruses) predominate during primary HIV-1 infections (PHI) while CXCR4-using viruses are isolated in less than 10% of PHI.

DESIGN

Six patients infected with an R5X4 virus, detected by a sensitive phenotypic assay during PHI, were matched with six patients infected with a pure R5 virus for sex, Fiebig stage, time of antiretroviral initiation and duration of follow-up.

METHODS

We used MiSeq ultra-deep sequencing to determine the composition of the virus quasispecies during PHI and at the end of follow-up (median time of follow-up: 12.5 years).

RESULTS

X4 viruses were detected by genetic analysis in three of six samples from the R5X4 group, accounting for 1.3-100% of the virus quasispecies, during PHI, and in four of six samples (accounting for 6.7-100%) at the end of follow-up. No X4 virus was detected in the R5 group during PHI and in only one patient (accounting for 1.2%) at the end of follow-up. The complexity of the virus quasispecies at the stage of PHI was higher in the R5X4 group than in the R5 group. Complexity increased from PHI to the end of follow-up in the R5 group but remained stable in the R5X4 group.

CONCLUSION

CXCR4-using viruses persisted in the peripheral blood mononuclear cells of several patients on suppressive antiretroviral therapy for a median duration of 12.5 years after PHI. The genetic complexity of HIV-1 evolved differently post-PHI in patients infected with R5X4 viruses from those infected with R5 viruses.

Exploring the hepatitis C virus genome using single molecule real-time sequencing

Single molecular real-time (SMRT) sequencing, also called third-generation sequencing, is a novel sequencing technique capable of generating extremely long contiguous sequence reads. While conventional short-read sequencing cannot evaluate the linkage of nucleotide substitutions distant from one another, SMRT sequencing can directly demonstrate linkage of nucleotide changes over a span of more than 20 kbp, and thus can be applied to directly examine the haplotypes of viruses or bacteria whose genome structures are changing in real time. In addition, an error correction method (circular consensus sequencing) has been established and repeated sequencing of a single-molecule DNA template can result in extremely high accuracy. The advantages of long read sequencing enable accurate determination of the haplotypes of individual viral clones. SMRT sequencing has been applied in various studies of viral genomes including determination of the full-length contiguous genome sequence of hepatitis C virus (HCV), targeted deep sequencing of the HCV NS5A gene, and assessment of heterogeneity among viral populations. Recently, the emergence of multi-drug resistant HCV viruses has become a significant clinical issue and has been also demonstrated using SMRT sequencing. In this review, we introduce the novel third-generation PacBio RSII/Sequel systems, compare them with conventional next-generation sequencers, and summarize previous studies in which SMRT sequencing technology has been applied for HCV genome analysis. We also refer to another long-read sequencing platform, nanopore sequencing technology, and discuss the advantages, limitations and future perspectives in using these third-generation sequencers for HCV genome analysis.

Next-generation sequencing of HIV-1 single genome amplicons

The analysis of HIV-1 sequences has helped understand the viral molecular epidemiology, monitor the development of antiretroviral drug resistance, and design candidate vaccines. The introduction of single genome amplification (SGA) has been a major advancement in the field, allowing for the characterization of multiple sequences per patient while preserving linkage among polymorphisms in the same viral genome copy. Sequencing of SGA amplicons is performed by capillary Sanger sequencing, which presents low throughput, requires a high amount of template, and is highly sensitive to template/primer mismatching. In order to meet the increasing demand for HIV-1 SGA amplicon sequencing, we have developed a platform based on benchtop next-generation sequencing (NGS) (IonTorrent) accompanied by a bioinformatics pipeline capable of running on computer resources commonly available at research laboratories. During assay validation, the NGS-based sequencing of 10 HIV-1 env SGA amplicons was fully concordant with Sanger sequencing. The field test was conducted on plasma samples from 10 US Navy and Marine service members with recent HIV-1 infection (sampling interval: 2005–2010; plasma viral load: 5,884–194,984 copies/ml). The NGS analysis of 101 SGA amplicons (median: 10 amplicons/individual) showed within-individual viral sequence profiles expected in individuals at this disease stage, including individuals with highly homogeneous quasispecies, individuals with two highly homogeneous viral lineages, and individuals with heterogeneous viral populations. In a scalability assessment using the Ion Chef automated system, 41/43 tested env SGA amplicons (95%) multiplexed on a single Ion 318 chip showed consistent gene-wide coverage >50×. With lower sample requirements and higher throughput, this approach is suitable to support the increasing demand for high-quality and cost-effective HIV-1 sequences in fields such as molecular epidemiology, and development of preventive and therapeutic strategies.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Universal Target Capture of HIV Sequences From NGS Libraries

Background: Global surveillance of viral sequence diversity is needed to keep pace with the constant evolution of HIV. Recent next generation sequencing (NGS) methods have realized the goal of sequencing circulating virus directly from patient specimens. Yet, a simple, universal approach that maximizes sensitivity and sequencing capacity remains elusive. Here we present a novel HIV enrichment strategy to yield near complete genomes from low viral load specimens. Methodology: A non-redundant biotin-labeled probe set (HIV-xGen; n = 652) was synthesized to tile all HIV-1 (groups M, N, O, and P) and HIV-2 (A and B) strains. Illumina Nextera barcoded libraries of either gene-specific or randomly primed cDNA derived from infected plasma were hybridized to probes in a single pool and unbound sequences were washed away. Captured viral cDNA was amplified by Illumina adaptor primers, sequenced on a MiSeq, and NGS reads were demultiplexed for alignment with CLC Bio software. Results: HIV-xGen probes selectively captured and amplified reads spanning the entirety of the HIV phylogenetic tree. HIV sequences clearly present in unenriched libraries of specimens but previously not observed due to high host background levels, insufficient sequencing depth or the extent of multiplexing, were now enriched by >1,000-fold. Thus, xGen selection not only substantially increased the depth of existing sequence, but also extended overall genome coverage by an average of 40%. We characterized 50 new, diverse HIV strains from clinical specimens and demonstrated a viral load cutoff of approximately log 3.5 copies/ml for full length coverage. Genome coverage was <20% for 5/10 samples with viral loads <log 3.5 copies/ml and >90% for 35/40 samples with higher viral loads. Conclusions: Characterization of >20 complete genomes at a time is now possible from a single probe hybridization and MiSeq run. With the versatility to capture all HIV strains and the sensitivity to detect low titer specimens, HIV-xGen will serve as an important tool for monitoring HIV sequence diversity.

geno2pheno[ngs-freq]: a genotypic interpretation system for identifying viral drug resistance using next-generation sequencing data

Identifying resistance to antiretroviral drugs is crucial for ensuring the successful treatment of patients infected with viruses such as human immunodeficiency virus (HIV) or hepatitis C virus (HCV). In contrast to Sanger sequencing, next-generation sequencing (NGS) can detect resistance mutations in minority populations. Thus, genotypic resistance testing based on NGS data can offer novel, treatment-relevant insights. Since existing web services for analyzing resistance in NGS samples are subject to long processing times and follow strictly rules-based approaches, we developed geno2pheno[ngs-freq], a web service for rapidly identifying drug resistance in HIV-1 and HCV samples. By relying on frequency files that provide the read counts of nucleotides or codons along a viral genome, the time-intensive step of processing raw NGS data is eliminated. Once a frequency file has been uploaded, consensus sequences are generated for a set of user-defined prevalence cutoffs, such that the constructed sequences contain only those nucleotides whose codon prevalence exceeds a given cutoff. After locally aligning the sequences to a set of references, resistance is predicted using the well-established approaches of geno2pheno[resistance] and geno2pheno[hcv]. geno2pheno[ngs-freq] can assist clinical decision making by enabling users to explore resistance in viral populations with different abundances and is freely available at http://ngs.geno2pheno.org.

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.

Blood CXCR3+ CD4 T Cells Are Enriched in Inducible Replication Competent HIV in Aviremic Antiretroviral Therapy-Treated Individuals

We recently demonstrated that lymph nodes (LNs) PD-1⁺/T follicular helper (Tfh) cells from antiretroviral therapy (ART)-treated HIV-infected individuals were enriched in cells containing replication competent virus. However, the distribution of cells containing inducible replication competent virus has been only partially elucidated in blood memory CD4 T-cell populations including the Tfh cell counterpart circulating in blood (cTfh). In this context, we have investigated the distribution of (1) total HIV-infected cells and (2) cells containing replication competent and infectious virus within various blood and LN memory CD4 T-cell populations of conventional antiretroviral therapy (cART)-treated HIV-infected individuals. In the present study, we show that blood CXCR3-expressing memory CD4 T cells are enriched in cells containing inducible replication competent virus and contributed the most to the total pool of cells containing replication competent and infectious virus in blood. Interestingly, subsequent proviral sequence analysis did not indicate virus compartmentalization between blood and LN CD4 T-cell populations, suggesting dynamic interchanges between the two compartments. We then investigated whether the composition of blood HIV reservoir may reflect the polarization of LN CD4 T cells at the time of reservoir seeding and showed that LN PD-1⁺ CD4 T cells of viremic untreated HIV-infected individuals expressed significantly higher levels of CXCR3 as compared to CCR4 and/or CCR6, suggesting that blood CXCR3-expressing CD4 T cells may originate from LN PD-1⁺ CD4 T cells. Taken together, these results indicate that blood CXCR3-expressing CD4 T cells represent the major blood compartment containing inducible replication competent virus in treated aviremic HIV-infected individuals.

A heterogeneous human immunodeficiency virus-like particle (VLP) formulation produced by a novel vector system

First identified as the etiological agent behind Acquired Immunodeficiency Syndrome (AIDS) in the early 1980s, HIV-1 has continued to spread into a global pandemic and major public health concern. Despite the success of antiretroviral therapy at reducing HIV-1 viremia and preventing the dramatic CD4+ T-cell collapse, infected individuals remain HIV positive for life. Unfortunately, it is increasingly clear that natural immunity is not, and may never be, protective against this pathogen. Therefore, efficacious vaccine interventions, which can either prevent infection or eradicate the latent viral reservoir and effect cure, are a major medical priority. Here we describe the development of a safe vaccine platform, currently being utilized in on-going prophylactic and therapeutic preclinical studies and consisting of highly heterogeneous virus-like particle formulations that represent the virus diversity within infected individuals. These VLPs contain no 5'LTR, no functional integrase, and have a severely mutated stem loop 1-thereby preventing any potential reverse transcription, integration, and RNA packaging. Furthermore, we demonstrate that these VLPs are morphologically identical to wild-type virus with polyvalent Env in a functional form. Finally, we show that the VLPs are antigenic and capable of generating strong immune recall responses.

Higher sequence diversity in the vaginal tract than in blood at early HIV-1 infection

In the majority of cases, human immunodeficiency virus type 1 (HIV-1) infection is transmitted through sexual intercourse. A single founder virus in the blood of the newly infected donor emerges from a genetic bottleneck, while in rarer instances multiple viruses are responsible for systemic infection. We sought to characterize the sequence diversity at early infection, between two distinct anatomical sites; the female reproductive tract vs. systemic compartment. We recruited 72 women from Uganda and Zimbabwe within seven months of HIV-1 infection. Using next generation deep sequencing, we analyzed the total genetic diversity within the C2-V3-C3 envelope region of HIV-1 isolated from the female genital tract at early infection and compared this to the diversity of HIV-1 in plasma. We then compared intra-patient viral diversity in matched cervical and blood samples with three or seven months post infection. Genetic analysis of the C2-V3-C3 region of HIV-1 env revealed that early HIV-1 isolates within blood displayed a more homogeneous genotype (mean 1.67 clones, range 1–5 clones) than clones in the female genital tract (mean 5.7 clones, range 3–10 clones) (p<0.0001). The higher env diversity observed within the genital tract compared to plasma was independent of HIV-1 subtype (A, C and D). Our analysis of early mucosal infections in women revealed high HIV-1 diversity in the vaginal tract but few transmitted clones in the blood. These novel in vivo finding suggest a possible mucosal sieve effect, leading to the establishment of a homogenous systemic infection.

During chronic HIV-1 infection, high viral diversity can be found in the blood and semen of donors. However, a single HIV-1 clone establishes productive infection in the recipient following heterosexual transmission. To investigate the genetic bottleneck occurring at the earliest stages of heterosexual HIV-1 transmission, we characterized the HIV-1 envelope sequence diversity at very early and early stages of infection in the female reproductive tract and matched plasma samples from a cohort of Ugandan and Zimbabwean women. A more diverse viral population was observed in the endocervical swab samples compared to plasma. Endocervical samples harbored a larger number of viral clones, while in the majority of plasma samples only a single clone was present early in infection. Interestingly, these observations were independent of HIV-1 subtype, hormonal contraceptive use or the number of sex acts and partners. Furthermore, in the cases of higher HIV-1 diversity in the blood during early infection, faster CD4 T cell decline were observed during chronic disease suggesting faster disease progression. Our findings provide novel in vivo evidence for the existence of an intra-patient genetic bottleneck restricting the HIV-1 from the vaginal tract to the blood during early heterosexual HIV-1 transmission.

HIV-1 envelope sequence-based diversity measures for identifying recent infections

Identifying recent HIV-1 infections is crucial for monitoring HIV-1 incidence and optimizing public health prevention efforts. To identify recent HIV-1 infections, we evaluated and compared the performance of 4 sequence-based diversity measures including percent diversity, percent complexity, Shannon entropy and number of haplotypes targeting 13 genetic segments within the env gene of HIV-1. A total of 597 diagnostic samples obtained in 2013 and 2015 from recently and chronically HIV-1 infected individuals were selected. From the selected samples, 249 (134 from recent versus 115 from chronic infections) env coding regions, including V1-C5 of gp120 and the gp41 ectodomain of HIV-1, were successfully amplified and sequenced by next generation sequencing (NGS) using the Illumina MiSeq platform. The ability of the four sequence-based diversity measures to correctly identify recent HIV infections was evaluated using the frequency distribution curves, median and interquartile range and area under the curve (AUC) of the receiver operating characteristic (ROC). Comparing the median and interquartile range and evaluating the frequency distribution curves associated with the 4 sequence-based diversity measures, we observed that the percent diversity, number of haplotypes and Shannon entropy demonstrated significant potential to discriminate recent from chronic infections (p<0.0001). Using the AUC of ROC analysis, only the Shannon entropy measure within three HIV-1 env segments could accurately identify recent infections at a satisfactory level. The env segments were gp120 C2_1 (AUC = 0.806), gp120 C2_3 (AUC = 0.805) and gp120 V3 (AUC = 0.812). Our results clearly indicate that the Shannon entropy measure represents a useful tool for predicting HIV-1 infection recency.

Application of viromics: a new approach to the understanding of viral infections in humans

This review is focused at exploring the strengths of modern technology driven data compiled in the areas of virus gene sequencing, virus protein structures and their implication to viral diagnosis and therapy. The information for virome analysis (viromics) is generated by the study of viral genomes (entire nucleotide sequence) and viral genes (coding for protein). Presently, the study of viral infectious diseases in terms of etiopathogenesis and development of newer therapeutics is undergoing rapid changes. Currently, viromics relies on deep sequencing, next generation sequencing (NGS) data and public domain databases like GenBank and unique virus specific databases. Two commonly used NGS platforms: Illumina and Ion Torrent, recommend maximum fragment lengths of about 300 and 400 nucleotides for analysis respectively. Direct detection of viruses in clinical samples is now evolving using these methods. Presently, there are a considerable number of good treatment options for HBV/HIV/HCV. These viruses however show development of drug resistance. The drug susceptibility regions of the genomes are sequenced and the prediction of drug resistance is now possible from 3 public domains available on the web. This has been made possible through advances in the technology with the advent of high throughput sequencing and meta-analysis through sophisticated and easy to use software and the use of high speed computers for bioinformatics. More recently NGS technology has been improved with single-molecule real-time sequencing. Here complete long reads can be obtained with less error overcoming a limitation of the NGS which is inherently prone to software anomalies that arise in the hands of personnel without adequate training. The development in understanding the viruses in terms of their genome, pathobiology, transcriptomics and molecular epidemiology constitutes viromics. It could be stated that these developments will bring about radical changes and advancement especially in the field of antiviral therapy and diagnostic virology.

Intrasubtype B HIV-1 Superinfection Correlates with Delayed Neutralizing Antibody Response

Understanding whether the neutralizing antibody (NAb) response impacts HIV-1 superinfection and how superinfection subsequently modulates the NAb response can help clarify correlates of protection from HIV exposures and better delineate pathways of NAb development. We examined associations between the development of NAb and the occurrence of superinfection in a well-characterized, antiretroviral therapy (ART)-naive, primary infection cohort of men who have sex with men. Deep sequencing was applied to blood plasma samples from the cohort to detect cases of superinfection. We compared the NAb activity against autologous and heterologous viruses between 10 participants with intrasubtype B superinfection and 19 monoinfected controls, matched to duration of infection and risk behavior. Three to 6 months after primary infection, individuals who would later become superinfected had significantly weaker NAb activity against tier 1 subtype B viruses (P = 0.003 for SF-162 and P = 0.017 for NL4-3) and marginally against autologous virus (P = 0.054). Lower presuperinfection NAb responses correlated with weaker gp120 binding and lower plasma total IgG titers. Soon after superinfection, the NAb response remained lower, but between 2 and 3 years after primary infection, NAb levels strengthened and reached those of controls. Superinfecting viruses were typically not susceptible to neutralization by presuperinfection plasma. These observations suggest that recently infected individuals with a delayed NAb response against primary infecting and tier 1 subtype B viruses are more susceptible to superinfection.IMPORTANCE Our findings suggest that within the first year after HIV infection, a relatively weak neutralizing antibody response against primary and subtype-specific neutralization-sensitive viruses increases susceptibility to superinfection in the face of repeated exposures. As natural infection progresses, the immune response strengthens significantly in some superinfected individuals. These findings will inform HIV vaccine design by providing testable correlates of protection from initial HIV infection.

Evolution of selective-sequencing approaches for virus discovery and virome analysis

•

Description of virus enrichment techniques for metagenomics based virome analysis.

•

Usefulness of recently developed virome capture sequencing techniques.

•

Perspective on negative and positive selection approaches for virome analysis.

Recent advances in sequencing technologies have transformed the field of virus discovery and virome analysis. Once mostly confined to the traditional Sanger sequencing based individual virus discovery, is now entirely replaced by high throughput sequencing (HTS) based virus metagenomics that can be used to characterize the nature and composition of entire viromes. To better harness the potential of HTS for the study of viromes, sample preparation methodologies use different approaches to exclude amplification of non-viral components that can overshadow low-titer viruses. These virus-sequence enrichment approaches mostly focus on the sample preparation methods, like enzymatic digestion of non-viral nucleic acids and size exclusion of non-viral constituents by column filtration, ultrafiltration or density gradient centrifugation. However, recently a new approach of virus-sequence enrichment called virome-capture sequencing, focused on the amplification or HTS library preparation stage, was developed to increase the ability of virome characterization. This new approach has the potential to further transform the field of virus discovery and virome analysis, but its technical complexity and sequence-dependence warrants further improvements. In this review we discuss the different methods, their applications and evolution, for selective sequencing based virome analysis and also propose refinements needed to harness the full potential of HTS for virome analysis.

Emergence of CXCR4-tropic HIV-1 variants followed by rapid disease progression in hemophiliac slow progressors

OBJECTIVE

The association between emergence of CXCR4-tropic HIV-1 variants (X4 variants) and disease progression of HIV-1 infection has been reported. However, it is not known whether the emergence of X4 variants is the cause or result of HIV-1 disease progression. We tried to answer this question.

DESIGN

HIV-1 env sequences around the V3 region were analyzed in serially stocked samples in order to determine whether X4 variants emerged before or after the fall in CD4+ T-cell count.

METHODS

The study subjects were five HIV-1-infected hemophiliac slow progressors. Deep sequencing around the HIV-1 env V3 region was conducted in duplicate. Tropism was predicted by geno2pheno [coreceptor] 2.5 with cutoff value of false positive ratio at <5%. When X4 variant was identified in the latest stocked sample before the introduction of antiretroviral therapy, we checked viral genotype in previously stocked samples to determine the time of emergence of X4 variants.

RESULTS

Emergence of X4 variants was noted in two of the five patients when their CD4+ T-cell counts were still high. The rate of decrease of CD4+ T-cell count or of rise of HIV-1 load accelerated significantly after the emergence of X4 variants in these two cases. Phylogenetic analysis showed that these X4 variants emerged from CCR5-tropic HIV-1 viruses with several amino acid changes in the V3 region.

CONCLUSIONS

The emergence of X4 variants preceded HIV-1 disease progression in two hemophiliac slow progressors.

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.

Statistical method to compare massive parallel sequencing pipelines

Background

Today, sequencing is frequently carried out by Massive Parallel Sequencing (MPS) that cuts drastically sequencing time and expenses. Nevertheless, Sanger sequencing remains the main validation method to confirm the presence of variants. The analysis of MPS data involves the development of several bioinformatic tools, academic or commercial. We present here a statistical method to compare MPS pipelines and test it in a comparison between an academic (BWA-GATK) and a commercial pipeline (TMAP-NextGENe®), with and without reference to a gold standard (here, Sanger sequencing), on a panel of 41 genes in 43 epileptic patients. This method used the number of variants to fit log-linear models for pairwise agreements between pipelines. To assess the heterogeneity of the margins and the odds ratios of agreement, four log-linear models were used: a full model, a homogeneous-margin model, a model with single odds ratio for all patients, and a model with single intercept. Then a log-linear mixed model was fitted considering the biological variability as a random effect.

Results

Among the 390,339 base-pairs sequenced, TMAP-NextGENe® and BWA-GATK found, on average, 2253.49 and 1857.14 variants (single nucleotide variants and indels), respectively. Against the gold standard, the pipelines had similar sensitivities (63.47% vs. 63.42%) and close but significantly different specificities (99.57% vs. 99.65%; p < 0.001). Same-trend results were obtained when only single nucleotide variants were considered (99.98% specificity and 76.81% sensitivity for both pipelines).

Conclusions

The method allows thus pipeline comparison and selection. It is generalizable to all types of MPS data and all pipelines.

Electronic supplementary material

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

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.

Infection of rhesus macaques with a pool of simian immunodeficiency virus with the envelope genes from acute HIV-1 infections

Background

New simian–human immunodeficiency chimeric viruses with an HIV-1 env (SHIVenv) are critical for studies on HIV pathogenesis, vaccine development, and microbicide testing. Macaques are typically exposed to single CCR5-using SHIVenv which in most instances does not reflect the conditions during acute/early HIV infection (AHI) in humans. Instead of individual and serial testing new SHIV constructs, a pool of SHIVenv_B derived from 16 acute HIV-1 infections were constructed using a novel yeast-based SHIV cloning approach and then used to infect macaques.

Results

Even though none of the 16 SHIVenvs contained the recently reported mutations in env genes that could significantly enhance their binding affinity to RhCD4, one SHIVenv (i.e. SHIVenv_B3-PRB926) established infection in macaques exposed to this pool. AHI SHIVenv_B viruses as well as their HIVenv_B counterparts were analyzed for viral protein content, function, and fitness to identify possible difference between SHIVenv_B3-PRB926 and the other 15 SHIVenvs in the pool. All of the constructs produced SHIV or HIV chimeric with wild type levels of capsid (p27 and p24) content, reverse transcriptase (RT) activity, and expressed envelope glycoproteins that could bind to cell receptors CD4/CCR5 and mediate virus entry. HIV-1env_B chimeric viruses were propagated in susceptible cell lines but the 16 SHIVenv_B variants showed only limited replication in macaque peripheral blood mononuclear cells (PBMCs) and 174×CEM.CCR5 cell line. AHI chimeric viruses including HIVenv_B3 showed only minor variations in cell entry efficiency and kinetics as well as replicative fitness in human PBMCs. Reduced number of N-link glycosylation sites and slightly greater CCR5 affinity/avidity was the only distinguishing feature of env_B3 versus other AHI env’s in the pool, a feature also observed in the HIV establishing new infections in humans.

Conclusion

Despite the inability to propagate in primary cells and cell lines, a pool of 16 SHIVenv viruses could establish infection but only one virus, SHIVenv_B3 was isolated in the macaque and then shown to repeatedly infected macaques. This SHIVenv_B3 virus did not show any distinct phenotypic property from the other 15 SHIVenv viruses but did have the fewest N-linked glycosylation sites.

Electronic supplementary material

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

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.

HIV-1 drug resistance and resistance testing

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

Deep Sequencing of the HIV-1 env Gene Reveals Discrete X4 Lineages and Linkage Disequilibrium between X4 and R5 Viruses in the V1/V2 and V3 Variable Regions

HIV-1 requires the CD4 receptor and a coreceptor (CCR5 [R5 phenotype] or CXCR4 [X4 phenotype]) to enter cells. Coreceptor tropism can be assessed by either phenotypic or genotypic analysis, the latter using bioinformatics algorithms to predict tropism based on the env V3 sequence. We used the Primer ID sequencing strategy with the MiSeq sequencing platform to reveal the structure of viral populations in the V1/V2 and C2/V3 regions of the HIV-1 env gene in 30 late-stage and 6 early-stage subjects. We also used endpoint dilution PCR followed by cloning of env genes to create pseudotyped virus to explore the link between genotypic predictions and phenotypic assessment of coreceptor usage. We found out that the most stringently sequence-based calls of X4 variants (Geno2Pheno false-positive rate [FPR] of ≤2%) formed distinct lineages within the viral population, and these were detected in 24 of 30 late-stage samples (80%), which was significantly higher than what has been seen previously by using other approaches. Non-X4 lineages were not skewed toward lower FPR scores in X4-containing populations. Phenotypic assays showed that variants with an intermediate FPR (2 to 20%) could be either X4/dual-tropic or R5 variants, although the X4 variants made up only about 25% of the lineages with an FPR of <10%, and these variants carried a distinctive sequence change. Phylogenetic analysis of both the V1/V2 and C2/V3 regions showed evidence of recombination within but very little recombination between the X4 and R5 lineages, suggesting that these populations are genetically isolated.

IMPORTANCE

Primer ID sequencing provides a novel approach to study genetic structures of viral populations. X4 variants may be more prevalent than previously reported when assessed by using next-generation sequencing (NGS) and with a greater depth of sampling than single-genome amplification (SGA). Phylogenetic analysis to identify lineages of sequences with intermediate FPR values may provide additional information for accurately predicting X4 variants by using V3 sequences. Limited recombination occurs between X4 and R5 lineages, suggesting that X4 and R5 variants are genetically isolated and may be replicating in different cell types or that X4/R5 recombinants have reduced fitness.

Multiplexed next-generation sequencing and de novo assembly to obtain near full-length HIV-1 genome from plasma virus

Analysing the HIV-1 near full-length genome (HIV-NFLG) facilitates new understanding into the diversity of virus population dynamics at individual or population level. In this study we developed a simple but high-throughput next generation sequencing (NGS) protocol for HIV-NFLG using clinical specimens and validated the method against an external quality control (EQC) panel. Clinical specimens (n=105) were obtained from three cohorts from two highly conserved HIV-1C epidemics (India and Ethiopia) and one diverse epidemic (Sweden). Additionally an EQC panel (n=10) was used to validate the protocol. HIV-NFLG was performed amplifying the HIV-genome (Gag-to-nef) in two fragments. NGS was performed using the Illumina HiSeq2500 after multiplexing 24 samples, followed by de novo assembly in Iterative Virus Assembler or VICUNA. Subtyping was carried out using several bioinformatics tools. Amplification of HIV-NFLG has 90% (95/105) success-rate in clinical specimens. NGS was successful in all clinical specimens (n=45) and EQA samples (n=10) attempted. The mean error for mutations for the EQC panel viruses were <1%. Subtyping identified two as A1C recombinant. Our results demonstrate the feasibility of a simple NGS-based HIV-NFLG that can potentially be used in the molecular surveillance for effective identification of subtypes and transmission clusters for operational public health intervention.

Rapid Sequencing of Complete env Genes from Primary HIV-1 Samples

The ability to study rapidly evolving viral populations has been constrained by the read length of next-generation sequencing approaches and the sampling depth of single-genome amplification methods. Here, we develop and characterize a method using Pacific Biosciences' Single Molecule, Real-Time (SMRT®) sequencing technology to sequence multiple, intact full-length human immunodeficiency virus-1 env genes amplified from viral RNA populations circulating in blood, and provide computational tools for analyzing and visualizing these data.

Evolution of coreceptor utilization to escape CCR5 antagonist therapy

The HIV-1 envelope interacts with coreceptors CCR5 and CXCR4 in a dynamic, multi-step process, its molecular details not clearly delineated. Use of CCR5 antagonists results in tropism shift and therapeutic failure. Here we describe a novel approach using full-length patient-derived gp160 quasispecies libraries cloned into HIV-1 molecular clones, their separation based on phenotypic tropism in vitro, and deep sequencing of the resultant variants for structure-function analyses. Analysis of functionally validated envelope sequences from patients who failed CCR5 antagonist therapy revealed determinants strongly associated with coreceptor specificity, especially at the gp120-gp41 and gp41-gp41 interaction surfaces that invite future research on the roles of subunit interaction and envelope trimer stability in coreceptor usage. This study identifies important structure-function relationships in HIV-1 envelope, and demonstrates proof of concept for a new integrated analysis method that facilitates laboratory discovery of resistant mutants to aid in development of other therapeutic agents.

Highly Accurate Structure-Based Prediction of HIV-1 Coreceptor Usage Suggests Intermolecular Interactions Driving Tropism

HIV-1 entry into host cells is mediated by interactions between the V3-loop of viral glycoprotein gp120 and chemokine receptor CCR5 or CXCR4, collectively known as HIV-1 coreceptors. Accurate genotypic prediction of coreceptor usage is of significant clinical interest and determination of the factors driving tropism has been the focus of extensive study. We have developed a method based on nonlinear support vector machines to elucidate the interacting residue pairs driving coreceptor usage and provide highly accurate coreceptor usage predictions. Our models utilize centroid-centroid interaction energies from computationally derived structures of the V3-loop:coreceptor complexes as primary features, while additional features based on established rules regarding V3-loop sequences are also investigated. We tested our method on 2455 V3-loop sequences of various lengths and subtypes, and produce a median area under the receiver operator curve of 0.977 based on 500 runs of 10-fold cross validation. Our study is the first to elucidate a small set of specific interacting residue pairs between the V3-loop and coreceptors capable of predicting coreceptor usage with high accuracy across major HIV-1 subtypes. The developed method has been implemented as a web tool named CRUSH, CoReceptor USage prediction for HIV-1, which is available at http://ares.tamu.edu/CRUSH/.

A Follow-Up of the Multicenter Collaborative Study on HIV-1 Drug Resistance and Tropism Testing Using 454 Ultra Deep Pyrosequencing

BACKGROUND

Ultra deep sequencing is of increasing use not only in research but also in diagnostics. For implementation of ultra deep sequencing assays in clinical laboratories for routine diagnostics, intra- and inter-laboratory testing are of the utmost importance.

METHODS

A multicenter study was conducted to validate an updated assay design for 454 Life Sciences' GS FLX Titanium system targeting protease/reverse transcriptase (RTP) and env (V3) regions to identify HIV-1 drug-resistance mutations and determine co-receptor use with high sensitivity. The study included 30 HIV-1 subtype B and 6 subtype non-B samples with viral titers (VT) of 3,940-447,400 copies/mL, two dilution series (52,129-1,340 and 25,130-734 copies/mL), and triplicate samples. Amplicons spanning PR codons 10-99, RT codons 1-251 and the entire V3 region were generated using barcoded primers. Analysis was performed using the GS Amplicon Variant Analyzer and geno2pheno for tropism. For comparison, population sequencing was performed using the ViroSeq HIV-1 genotyping system.

RESULTS

The median sequencing depth across the 11 sites was 1,829 reads per position for RTP (IQR 592-3,488) and 2,410 for V3 (IQR 786-3,695). 10 preselected drug resistant variants were measured across sites and showed high inter-laboratory correlation across all sites with data (P<0.001). The triplicate samples of a plasmid mixture confirmed the high inter-laboratory consistency (mean% ± stdev: 4.6 ±0.5, 4.8 ±0.4, 4.9 ±0.3) and revealed good intra-laboratory consistency (mean% range ± stdev range: 4.2-5.2 ± 0.04-0.65). In the two dilutions series, no variants >20% were missed, variants 2-10% were detected at most sites (even at low VT), and variants 1-2% were detected by some sites. All mutations detected by population sequencing were also detected by UDS.

CONCLUSIONS

This assay design results in an accurate and reproducible approach to analyze HIV-1 mutant spectra, even at variant frequencies well below those routinely detectable by population sequencing.

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.

Multiregion deep sequencing of hepatitis C virus: An improved approach for genetic relatedness studies

Hepatitis C virus (HCV) is a major public health problem that affects more than 180 million people worldwide. Identification of HCV transmission networks is of critical importance for disease control. HCV related cases are often difficult to identify due to the characteristic long incubation period and lack of symptoms during the acute phase of the disease, making it challenging to link related cases to a common source of infection. Additionally, HCV transmission chains are difficult to trace back since viral variants from epidemiologically linked cases are genetically related but rarely identical. Genetic relatedness studies primarily rely on information obtained from the rapidly evolving HCV hypervariable region 1 (HVR1). However, in some instances, the rapid divergence of this region can lead to loss of genetic links between related isolates, which represents an important challenge for outbreak investigations and genetic relatedness studies. Sequencing of multiple and longer sub-genomic regions has been proposed as an alternative to overcome the limitations imposed by the rapid molecular evolution of the HCV HVR1. Additionally, conventional molecular approaches required to characterize the HCV intra-host genetic variation are laborious, time-consuming, and expensive while providing limited information about the composition of the viral population. Next generation sequencing (NGS) approaches enormously facilitate the characterization of the HCV intra-host population by detecting rare variants at much lower frequencies. Thus, NGS approaches using multiple sub-genomic regions should improve the characterization of the HCV intra-host population. Here, we explore the usefulness of multiregion sequencing using a NGS platform for genetic relatedness studies among HCV cases.

A Pan-HIV Strategy for Complete Genome Sequencing

Molecular surveillance is essential to monitor HIV diversity and track emerging strains. We have developed a universal library preparation method (HIV-SMART [i.e.,switchingmechanismat 5' end ofRNAtranscript]) for next-generation sequencing that harnesses the specificity of HIV-directed priming to enable full genome characterization of all HIV-1 groups (M, N, O, and P) and HIV-2. Broad application of the HIV-SMART approach was demonstrated using a panel of diverse cell-cultured virus isolates. HIV-1 non-subtype B-infected clinical specimens from Cameroon were then used to optimize the protocol to sequence directly from plasma. When multiplexing 8 or more libraries per MiSeq run, full genome coverage at a median ∼2,000× depth was routinely obtained for either sample type. The method reproducibly generated the same consensus sequence, consistently identified viral sequence heterogeneity present in specimens, and at viral loads of ≤4.5 log copies/ml yielded sufficient coverage to permit strain classification. HIV-SMART provides an unparalleled opportunity to identify diverse HIV strains in patient specimens and to determine phylogenetic classification based on the entire viral genome. Easily adapted to sequence any RNA virus, this technology illustrates the utility of next-generation sequencing (NGS) for viral characterization and surveillance.

Utility of Metagenomic Next-Generation Sequencing for Characterization of HIV and Human Pegivirus Diversity

Given the dynamic changes in HIV-1 complexity and diversity, next-generation sequencing (NGS) has the potential to revolutionize strategies for effective HIV global surveillance. In this study, we explore the utility of metagenomic NGS to characterize divergent strains of HIV-1 and to simultaneously screen for other co-infecting viruses. Thirty-five HIV-1-infected Cameroonian blood donor specimens with viral loads of >4.4 log₁₀ copies/ml were selected to include a diverse representation of group M strains. Random-primed NGS libraries, prepared from plasma specimens, resulted in greater than 90% genome coverage for 88% of specimens. Correct subtype designations based on NGS were concordant with sub-region PCR data in 31 of 35 (89%) cases. Complete genomes were assembled for 25 strains, including circulating recombinant forms with relatively limited data available (7 CRF11_cpx, 2 CRF13_cpx, 1 CRF18_cpx, and 1 CRF37_cpx), as well as 9 unique recombinant forms. HPgV (formerly designated GBV-C) co-infection was detected in 9 of 35 (25%) specimens, of which eight specimens yielded complete genomes. The recovered HPgV genomes formed a diverse cluster with genotype 1 sequences previously reported from Ghana, Uganda, and Japan. The extensive genome coverage obtained by NGS improved accuracy and confidence in phylogenetic classification of the HIV-1 strains present in the study population relative to conventional sub-region PCR. In addition, these data demonstrate the potential for metagenomic analysis to be used for routine characterization of HIV-1 and identification of other viral co-infections.

HIV-1 Tropism Determines Different Mutation Profiles in Proviral DNA

In order to establish new infections HIV-1 particles need to attach to receptors expressed on the cellular surface. HIV-1 particles interact with a cell membrane receptor known as CD4 and subsequently with another cell membrane molecule known as a co-receptor. Two major different co-receptors have been identified: C-C chemokine Receptor type 5 (CCR5) and C-X-C chemokine Receptor type 4 (CXCR4) Previous reports have demonstrated cellular modifications upon HIV-1 binding to its co-receptors including gene expression modulations. Here we investigated the effect of viral binding to either CCR5 or CXCR4 co-receptors on viral diversity after a single round of reverse transcription. CCR5 and CXCR4 pseudotyped viruses were used to infect non-stimulated and stimulated PBMCs and purified CD4 positive cells. We adopted the SOLiD methodology to sequence virtually the entire proviral DNA from all experimental infections. Infections with CCR5 and CXCR4 pseudotyped virus resulted in different patterns of genetic diversification. CCR5 virus infections produced extensive proviral diversity while in CXCR4 infections a more localized substitution process was observed. In addition, we present pioneering results of a recently developed method for the analysis of SOLiD generated sequencing data applicable to the study of viral quasi-species. Our findings demonstrate the feasibility of viral quasi-species evaluation by NGS methodologies. We presented for the first time strong evidence for a host cell driving mechanism acting on the HIV-1 genetic variability under the control of co-receptor stimulation. Additional investigations are needed to further clarify this question, which is relevant to viral diversification process and consequent disease progression.

Clinical value of ultradeep HIV-1 genotyping and tropism testing in late presenters with advanced disease

OBJECTIVE

This article aims to investigate if the detection of preexisting drug-resistant minority variant (DRMV) and/or X4 HIV-1 variants could improve the efficacy of first-line combined antiretroviral therapy (ART) in late presenters.

DESIGN

Post-hoc, combined analysis of two open-label, prospective, randomized clinical trials comparing first-line ART with efavirenz (EFV) vs. ritonavir-boosted protease inhibitor (PI/r)-based regimens in ART-naive, HIV-1-infected patients, with CD4 T-cell counts less than 100 cells/μl and wild-type HIV-1 by bulk sequencing.

METHODS

Pre-ART samples were reanalyzed for the presence of DRMVs and X4 HIV-1 using 454 sequencing. Kaplan-Meier curves and Cox regression were used to evaluate the association between X4 HIV and DRMVs and risk of virological failure.

RESULTS

From 141 evaluable patients, 57 received EFV, and 84 received PI/r, including first-line ART. Median pre-ART CD4 T-cell counts and HIV-1 RNA levels were 39 cells/μl and 257 424 copies/ml, respectively; 35.5% of patients had X4 HIV variants. Detection of DRMVs leading to an ART-specific cumulative HIVdb score of at least 10 increased the risk of virological failure in patients initiating EFV [log-rank P = 0.048, hazard ratio = 4.3 (95% confidence interval: 0.8, 25.0), P = 0.074], but not in those starting PI/r. Presence of X4 HIV did not affect virological outcomes, but was associated with impaired CD4 T-cell count recovery over 2 years (214 vs. 315 cells/μl with X4 vs. R5 HIV-1 tropism, respectively, P = 0.017).

CONCLUSION

Accounting for preexisting DRMVs may improve the outcomes of first-line nonnucleoside reverse transcriptase inhibitor-based ART in late presenters with advanced immune suppression. Presence of X4 HIV-1 at diagnosis predicts impaired immune restoration under ART.

Towards allele-level human leucocyte antigens genotyping - assessing two next-generation sequencing platforms: Ion Torrent Personal Genome Machine and Illumina MiSeq

Human leucocyte antigens (HLA) typing has been a challenge due to extreme polymorphism of the HLA genes and limitations of the current technologies and protocols used for their characterization. Recently, next-generation sequencing techniques have been shown to be a well-suited technology for the complete characterization of the HLA genes. However, a comprehensive assessment of the different platforms for HLA typing, describing the limitations and advantages of each of them, has not been presented. We have compared the Ion Torrent Personal Genome Machine (PGM) and Illumina MiSeq, currently the two most frequently used platforms for diagnostic applications, for a number of metrics including total output, quality score per position across the reads and error rates after alignment which can all affect the accuracy of HLA genotyping. For this purpose, we have used one homozygous and three heterozygous well-characterized samples, at HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1. The total output of bases produced by the MiSeq was higher, and they have higher quality scores and a lower overall error rate than the PGM. The MiSeq also has a higher fidelity when sequencing through homopolymer regions up to 9 bp in length. The need to set phase between distant polymorphic sites was more readily achieved with MiSeq using paired-end sequencing of fragments that are longer than those obtained with PGM. Additionally, we have assessed the workflows of the different platforms for complexity of sample preparation, sequencer operation and turnaround time. The effects of data quality and quantity can impact the genotyping results; having an adequate amount of good quality data to analyse will be imperative for confident HLA genotyping. The overall turnaround time can be very comparable between the two platforms; however, the complexity of sample preparation is higher with PGM, while the actual sequencing time is longer with MiSeq.

Primer ID Validates Template Sampling Depth and Greatly Reduces the Error Rate of Next-Generation Sequencing of HIV-1 Genomic RNA Populations

Validating the sampling depth and reducing sequencing errors are critical for studies of viral populations using next-generation sequencing (NGS). We previously described the use of Primer ID to tag each viral RNA template with a block of degenerate nucleotides in the cDNA primer. We now show that low-abundance Primer IDs (offspring Primer IDs) are generated due to PCR/sequencing errors. These artifactual Primer IDs can be removed using a cutoff model for the number of reads required to make a template consensus sequence. We have modeled the fraction of sequences lost due to Primer ID resampling. For a typical sequencing run, less than 10% of the raw reads are lost to offspring Primer ID filtering and resampling. The remaining raw reads are used to correct for PCR resampling and sequencing errors. We also demonstrate that Primer ID reveals bias intrinsic to PCR, especially at low template input or utilization. cDNA synthesis and PCR convert ca. 20% of RNA templates into recoverable sequences, and 30-fold sequence coverage recovers most of these template sequences. We have directly measured the residual error rate to be around 1 in 10,000 nucleotides. We use this error rate and the Poisson distribution to define the cutoff to identify preexisting drug resistance mutations at low abundance in an HIV-infected subject. Collectively, these studies show that >90% of the raw sequence reads can be used to validate template sampling depth and to dramatically reduce the error rate in assessing a genetically diverse viral population using NGS.

IMPORTANCE

Although next-generation sequencing (NGS) has revolutionized sequencing strategies, it suffers from serious limitations in defining sequence heterogeneity in a genetically diverse population, such as HIV-1 due to PCR resampling and PCR/sequencing errors. The Primer ID approach reveals the true sampling depth and greatly reduces errors. Knowing the sampling depth allows the construction of a model of how to maximize the recovery of sequences from input templates and to reduce resampling of the Primer ID so that appropriate multiplexing can be included in the experimental design. With the defined sampling depth and measured error rate, we are able to assign cutoffs for the accurate detection of minority variants in viral populations. This approach allows the power of NGS to be realized without having to guess about sampling depth or to ignore the problem of PCR resampling, while also being able to correct most of the errors in the data set.

Two-Year Follow-Up of Macaques Developing Intermittent Control of the Human Immunodeficiency Virus Homolog Simian Immunodeficiency Virus SIVmac251 in the Chronic Phase of Infection

Off-therapy control of viremia by HIV-infected individuals has been associated with two likely players: a restricted viral reservoir and an efficient cell-mediated immune response. We previously showed that a combination of highly suppressive antiretroviral therapy and two experimental drugs, i.e., auranofin and buthionine sulfoximine, was able to reduce the viral reservoir, elicit efficient cell-mediated antiviral responses, and induce intermittent posttherapy viral load control in chronically SIVmac251-infected macaques. We here show that the macaques that had received this drug combination and then stopped antiretroviral therapy were also able to maintain low numbers of activated CD4+ T cells at viral rebound. Moreover, these macaques consistently displayed low-level simian immunodeficiency virus (SIV) diversity, which was in line with the strong and broadly reactive cell-mediated immune responses against conserved Gag antigens. Extended follow-up showed that the two macaques that had received the complete drug combination remained healthy and did not develop AIDS in 2 years of follow-up after therapy suspension. This disease-free survival is longer than twice the average time of progression to AIDS in SIVmac251-infected rhesus macaques. These results suggest that limited numbers of activated T cells at viral rebound and subsequent development of broadly reactive cell-mediated responses may be interrelated in reducing the viral reservoir.

IMPORTANCE

The HIV reservoir in CD4+ T cells represents one main obstacle to HIV eradication. Recent studies, however, show that a drastic reduction of this reservoir is insufficient for inducing a functional cure of AIDS. In the present work, we thoroughly studied and subjected to long-term follow-up two macaques showing intermittent control of the virus following suspension of antiretroviral therapy plus an experimental antireservoir treatment, i.e., the gold salt auranofin and the investigational chemotherapeutic agent buthionione sulfoximine (BSO). We found that these drugs were able to decrease the number of activated CD4+ T cells, which are preferential targets for HIV infection. Then, efficient immune responses against the virus were developed in the macaques, which remained healthy during 2 years of follow-up. This result may furnish another building block for future attempts to cure HIV/AIDS.

From benchtop to desktop: important considerations when designing amplicon sequencing workflows

Amplicon sequencing has been the method of choice in many high-throughput DNA sequencing (HTS) applications. To date there has been a heavy focus on the means by which to analyse the burgeoning amount of data afforded by HTS. In contrast, there has been a distinct lack of attention paid to considerations surrounding the importance of sample preparation and the fidelity of library generation. No amount of high-end bioinformatics can compensate for poorly prepared samples and it is therefore imperative that careful attention is given to sample preparation and library generation within workflows, especially those involving multiple PCR steps. This paper redresses this imbalance by focusing on aspects pertaining to the benchtop within typical amplicon workflows: sample screening, the target region, and library generation. Empirical data is provided to illustrate the scope of the problem. Lastly, the impact of various data analysis parameters is also investigated in the context of how the data was initially generated. It is hoped this paper may serve to highlight the importance of pre-analysis workflows in achieving meaningful, future-proof data that can be analysed appropriately. As amplicon sequencing gains traction in a variety of diagnostic applications from forensics to environmental DNA (eDNA) it is paramount workflows and analytics are both fit for purpose.

Targeted single molecule sequencing methodology for ovarian hyperstimulation syndrome

Background

One of the most significant issues surrounding next generation sequencing is the cost and the difficulty assembling short read lengths. Targeted capture enrichment of longer fragments using single molecule sequencing (SMS) is expected to improve both sequence assembly and base-call accuracy but, at present, there are very few examples of successful application of these technologic advances in translational research and clinical testing. We developed a targeted single molecule sequencing (T-SMS) panel for genes implicated in ovarian response to controlled ovarian hyperstimulation (COH) for infertility.

Results

Target enrichment was carried out using droplet-base multiplex polymerase chain reaction (PCR) technology (RainDance®) designed to yield amplicons averaging 1 kb fragment size from candidate 44 loci (99.8% unique base-pair coverage). The total targeted sequence was 3.18 Mb per sample. SMS was carried out using single molecule, real-time DNA sequencing (SMRT® Pacific Biosciences®), average raw read length = 1178 nucleotides, 5% of the amplicons >6000 nucleotides). After filtering with circular consensus (CCS) reads, the mean read length was 3200 nucleotides (97% CCS accuracy). Primary data analyses, alignment and filtering utilized the Pacific Biosciences® SMRT portal. Secondary analysis was conducted using the Genome Analysis Toolkit for SNP discovery l and wANNOVAR for functional analysis of variants. Filtered functional variants 18 of 19 (94.7%) were further confirmed using conventional Sanger sequencing. CCS reads were able to accurately detect zygosity. Coverage within GC rich regions (i.e.VEGFR; 72% GC rich) was achieved by capturing long genomic DNA (gDNA) fragments and reading into regions that flank the capture regions. As proof of concept, a non-synonymous LHCGR variant captured in two severe OHSS cases, and verified by conventional sequencing.

Conclusions

Combining emulsion PCR-generated 1 kb amplicons and SMRT DNA sequencing permitted greater depth of coverage for T-SMS and facilitated easier sequence assembly. To the best of our knowledge, this is the first report combining emulsion PCR and T-SMS for long reads using human DNA samples, and NGS panel designed for biomarker discovery in OHSS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1451-2) contains supplementary material, which is available to authorized users.

Deep sequencing analysis of viral infection and evolution allows rapid and detailed characterization of viral mutant spectrum

Motivation: The study of RNA virus populations is a challenging task. Each population of RNA virus is composed of a collection of different, yet related genomes often referred to as mutant spectra or quasispecies. Virologists using deep sequencing technologies face major obstacles when studying virus population dynamics, both experimentally and in natural settings due to the relatively high error rates of these technologies and the lack of high performance pipelines. In order to overcome these hurdles we developed a computational pipeline, termed ViVan (Viral Variance Analysis). ViVan is a complete pipeline facilitating the identification, characterization and comparison of sequence variance in deep sequenced virus populations.

Results: Applying ViVan on deep sequenced data obtained from samples that were previously characterized by more classical approaches, we uncovered novel and potentially crucial aspects of virus populations. With our experimental work, we illustrate how ViVan can be used for studies ranging from the more practical, detection of resistant mutations and effects of antiviral treatments, to the more theoretical temporal characterization of the population in evolutionary studies.

Availability and implementation: Freely available on the web at http://www.vivanbioinfo.org

Contact: nshomron@post.tau.ac.il

Supplementary information:Supplementary data are available at Bioinformatics online.

The evolution of HIV-1 interactions with coreceptors and mannose C-type lectin receptors

The phenotype of human immunodeficiency virus type 1 (HIV-1) commonly evolves between and within infected individuals, at virus transmission, and during disease progression. This evolution includes altered interactions between the virus and its coreceptors, i.e., chemokine receptors, as well as mannose C-type lectin receptors (CLRs). Transmitted/founder viruses are predominantly restricted to CCR5, whereas the subsequent intrapatient evolution of HIV-1 coreceptor use during progressive disease can be subdivided into two distinct pathways. Accordingly, the CCR5-restricted virus population is either gradually replaced by virus variants able to use CXCR4 or evolves toward an altered, more flexible use of CCR5. Despite a strong dependency on these coreceptors for host cell entry, HIV-1 also interacts with other cell surface molecules during target cell attachment, including the CLRs. The virus interaction with the CLRs may result either in the efficient transfer of virus to CD4(+) T cells or in the degradation of the virus in endosomal compartments. The determinants of the diverse outcomes depend on which CLR is engaged and also on the glycan makeup of the envelope glycoproteins, which may evolve with the strength of the immune pressure during the disease course. With the current clinical introduction of CCR5 antagonists and the development of additional entry inhibitors, knowledge on the evolution and baseline characteristics of HIV-1 interactions with coreceptor and CLR interactions may play important roles for individualized and optimized treatment strategies. This review summarizes our current understanding of the evolution of HIV-1 interactions with these receptors.

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.