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

gCUP: rapid GPU-based HIV-1 co-receptor usage prediction for next-generation sequencing

SUMMARY

Next-generation sequencing (NGS) has a large potential in HIV diagnostics, and genotypic prediction models have been developed and successfully tested in the recent years. However, albeit being highly accurate, these computational models lack computational efficiency to reach their full potential. In this study, we demonstrate the use of graphics processing units (GPUs) in combination with a computational prediction model for HIV tropism. Our new model named gCUP, parallelized and optimized for GPU, is highly accurate and can classify >175 000 sequences per second on an NVIDIA GeForce GTX 460. The computational efficiency of our new model is the next step to enable NGS technologies to reach clinical significance in HIV diagnostics. Moreover, our approach is not limited to HIV tropism prediction, but can also be easily adapted to other settings, e.g. drug resistance prediction.

AVAILABILITY AND IMPLEMENTATION

The source code can be downloaded at http://www.heiderlab.de

CONTACT

d.heider@wz-straubing.de.

A simple structure-based model for the prediction of HIV-1 co-receptor tropism

Background

Human Immunodeficiency Virus 1 enters host cells through interaction of its V3 loop (which is part of the gp120 protein) with the host cell receptor CD4 and one of two co-receptors, namely CCR5 or CXCR4. Entry inhibitors binding the CCR5 co-receptor can prevent viral entry. As these drugs are only available for CCR5-using viruses, accurate prediction of this so-called co-receptor tropism is important in order to ensure an effective personalized therapy. With the development of next-generation sequencing technologies, it is now possible to sequence representative subpopulations of the viral quasispecies.

Results

Here we present T-CUP 2.0, a model for predicting co-receptor tropism. Based on our recently published T-CUP model, we developed a more accurate and even faster solution. Similarly to its predecessor, T-CUP 2.0 models co-receptor tropism using information of the electrostatic potential and hydrophobicity of V3-loops. However, extracting this information from a simplified structural vacuum-model leads to more accurate and faster predictions. The area-under-the-ROC-curve (AUC) achieved with T-CUP 2.0 on the training set is 0.968±0.005 in a leave-one-patient-out cross-validation. When applied to an independent dataset, T-CUP 2.0 has an improved prediction accuracy of around 3% when compared to the original T-CUP.

Conclusions

We found that it is possible to model co-receptor tropism in HIV-1 based on a simplified structure-based model of the V3 loop. In this way, genotypic prediction of co-receptor tropism is very accurate, fast and can be applied to large datasets derived from next-generation sequencing technologies. The reduced complexity of the electrostatic modeling makes T-CUP 2.0 independent from third-party software, making it easy to install and use.

Deep sequencing: becoming a critical tool in clinical virology

Population (Sanger) sequencing has been the standard method in basic and clinical DNA sequencing for almost 40 years; however, next-generation (deep) sequencing methodologies are now revolutionizing the field of genomics, and clinical virology is no exception. Deep sequencing is highly efficient, producing an enormous amount of information at low cost in a relatively short period of time. High-throughput sequencing techniques have enabled significant contributions to multiples areas in virology, including virus discovery and metagenomics (viromes), molecular epidemiology, pathogenesis, and studies of how viruses to escape the host immune system and antiviral pressures. In addition, new and more affordable deep sequencing-based assays are now being implemented in clinical laboratories. Here, we review the use of the current deep sequencing platforms in virology, focusing on three of the most studied viruses: human immunodeficiency virus (HIV), hepatitis C virus (HCV), and influenza virus.

Improved prediction of salvage antiretroviral therapy outcomes using ultrasensitive HIV-1 drug resistance testing

BACKGROUND

The clinical relevance of ultrasensitive human immunodeficiency virus type 1 (HIV-1) genotypic resistance testing in antiretroviral treatment (ART)-experienced individuals remains unknown.

METHODS

This was a retrospective, multicentre, cohort study in ART-experienced, HIV-1-infected adults who initiated salvage ART including, at least 1 ritonavir-boosted protease inhibitor, raltegravir or etravirine. Presalvage ART Sanger and 454 sequencing of plasma HIV-1 were used to generate separate genotypic sensitivity scores (GSS) using the HIVdb, ANRS, and REGA algorithms. Virological failure (VF) was defined as 2 consecutive HIV-1 RNA levels ≥200 copies/mL at least 12 weeks after salvage ART initiation, whereas subjects remained on the same ART. The ability of Sanger and 454-GSS to predict VF was assessed by receiver operating characteristic (ROC) curves and survival analyses.

RESULTS

The study included 132 evaluable subjects; 28 (21%) developed VF. Using HIVdb, 454 predicted VF better than Sanger sequencing in the ROC curve analysis (area under the curve: 0.69 vs 0.60, Delong test P = .029). Time to VF was shorter for subjects with 454-GSS < 3 vs 454-GSS ≥ 3 (Log-rank P = .003) but not significantly different between Sanger-GSS < 3 and ≥3. Factors independently associated with increased risk of VF in multivariate Cox regression were a 454-GSS < 3 (HR = 4.6, 95 CI, [1.5, 14.0], P = .007), and the number of previous antiretrovirals received (HR = 1.2 per additional drug, 95 CI, [1.1, 1.3], P = .001). Equivalent findings were obtained with the ANRS and REGA algorithms.

CONCLUSIONS

Ultrasensitive HIV-1 genotyping improves GSS-based predictions of virological outcomes of salvage ART relative to Sanger sequencing. This may improve the clinical management of ART-experienced subjects living with HIV-1.

CLINICAL TRIALS REGISTRATION

NCT01346878.

Targeted deep sequencing of HIV-1 using the IonTorrentPGM platform

The characterization of mixed HIV-1 populations is a key question in clinical and basic research settings. This can be achieved through targeted deep sequencing (TDS), where next-generation sequencing is used to examine in depth a sub-genomic region of interest. This study explores the suitability of IonTorrent PGM(LifeTechnologies) for the TDS-based analysis of HIV-1 evolution. Using laboratory reagents and primary specimens sampled at pre-peak viremia the error rates from misincorporation and in vitro recombination were <0.5%. The sequencing error rate was 2- to 3-fold higher in/around homopolymeric tracts, and could be discerned from true polymorphism using bidirectional sequencing. The limit of detection of complex variants was further lowered by using haplotyping. The application of this system was illustrated on primary samples from an individual infected with HIV-1 followed from pre-peak viremia through six months post-acquisition. TDS provided an augmented view of the extent of genetic diversity, the covariation among polymorphisms, the evolutionary pathways, and the boundaries of the mutational space explored by the viral swarm. Based on its performance, the system can be applied for the characterization of minor viral variants in support of studies of viral evolution, which can inform the rational design of the next generation of vaccines and therapeutics.

Genetic variability of porcine circovirus 2 in vaccinating and non-vaccinating commercial farms

Vaccines against porcine circovirus 2 (PCV2) are now widely used to control the diseases caused by the virus. Although the vaccines protect pigs against the disease, they do not lead to sterilizing immunity and therefore infections with PCV2 continue in farms. It is expected that, due to its high evolutionary rate, PCV2 can adapt quickly to environmental pressures such as vaccination. The goal of this study was to elucidate the molecular variation of PCV2 in relation to vaccination. PCV2 variability was investigated from samples of infected pigs from five farms where vaccination had never been applied and two farms where pigs had been vaccinated for at least 2 years. For the genetic analysis, full PCV2 genomes were amplified and subsequently pooled by vaccination status from serum of eight vaccinated, infected pigs and 16 non-vaccinated, infected pigs. Variability of viral populations was quantified using next-generation sequencing and subsequent bioinformatics analysis. The number of segregating sites was similar in the non-vaccinated (n=109) and vaccinated pools (n=96), but the distribution of these sites in the genome differed. Most notably, in the capsid gene, the number of segregating sites was observed only in the non-vaccinated population. Based on the structural analysis, it is expected that some low-frequency amino acids result in biologically low-fit viruses. On the contrary, D294 in replicase represents a novel amino acid which was dominant and unique in the vaccinated pool. This work showed that variable PCV2 populations were circulating in commercial farms, and that this variability was different in samples obtained from vaccinating and non-vaccinating farms.

Virological response to short-course maraviroc monotherapy does not predict viral tropism in HIV-1-infected treatment-naive patients

OBJECTIVES

We aimed to evaluate whether virological response to a short course of maraviroc monotherapy could predict HIV-1 tropism.

METHODS

A clinical trial was performed in HIV-1 treatment-naive patients infected with R5- or non-R5-tropic virus determined using the Trofile(®) assay, with >1000 HIV-1 RNA copies/mL. Maraviroc was administered for 10 days. Viral load was measured at baseline and days 4, 7, 10 and 28. The main outcome measurement was the decline in HIV-1 RNA at day 10. The trial was registered in the ClinicalTrials.gov database (NCT01060618; TROPISMVC).

RESULTS

Forty patients [30 R5 and 10 dual/mixed (D/M)] were recruited. There was a significant decrease in HIV-1 RNA after 10 days of maraviroc treatment in patients with R5-tropic virus (median 1.52 log10 RNA copies/mL; 95% CI 1.23-1.63; P < 0.0001), but also in patients with D/M-tropic virus (median 1.62 log(10) RNA copies/mL; 95% CI 0.33-1.88; P = 0.00024). The difference in the HIV-1 RNA decrease (-0.16 log(10) RNA copies/mL; 95% CI -0.53 to 0.22) was not significant (P = 0.410). A decrease >0.5 log(10) RNA copies/mL was found in 96.3% of patients with R5-tropic virus and in 70% of patients with D/M-tropic virus (P = 0.052). The differences were not significant when a decline of 1 log(10) RNA copies/mL was considered (92.6% versus 70%; P = 0.11).

CONCLUSIONS

Treatment-naive patients infected with R5- or D/M-tropic virus have similar virological responses to a short course of maraviroc monotherapy. This clinical test thus cannot be used as a surrogate marker of viral tropism in this population.

Comparison of illumina and 454 deep sequencing in participants failing raltegravir-based antiretroviral therapy

BACKGROUND

The impact of raltegravir-resistant HIV-1 minority variants (MVs) on raltegravir treatment failure is unknown. Illumina sequencing offers greater throughput than 454, but sequence analysis tools for viral sequencing are needed. We evaluated Illumina and 454 for the detection of HIV-1 raltegravir-resistant MVs.

METHODS

A5262 was a single-arm study of raltegravir and darunavir/ritonavir in treatment-naïve patients. Pre-treatment plasma was obtained from 5 participants with raltegravir resistance at the time of virologic failure. A control library was created by pooling integrase clones at predefined proportions. Multiplexed sequencing was performed with Illumina and 454 platforms at comparable costs. Illumina sequence analysis was performed with the novel snp-assess tool and 454 sequencing was analyzed with V-Phaser.

RESULTS

Illumina sequencing resulted in significantly higher sequence coverage and a 0.095% limit of detection. Illumina accurately detected all MVs in the control library at ≥0.5% and 7/10 MVs expected at 0.1%. 454 sequencing failed to detect any MVs at 0.1% with 5 false positive calls. For MVs detected in the patient samples by both 454 and Illumina, the correlation in the detected variant frequencies was high (R2 = 0.92, P<0.001). Illumina sequencing detected 2.4-fold greater nucleotide MVs and 2.9-fold greater amino acid MVs compared to 454. The only raltegravir-resistant MV detected was an E138K mutation in one participant by Illumina sequencing, but not by 454.

CONCLUSIONS

In participants of A5262 with raltegravir resistance at virologic failure, baseline raltegravir-resistant MVs were rarely detected. At comparable costs to 454 sequencing, Illumina demonstrated greater depth of coverage, increased sensitivity for detecting HIV MVs, and fewer false positive variant calls.

HIV-1 tropism testing and clinical management of CCR5 antagonists: Quebec review and recommendations

HIV-1 tropism assays play a crucial role in determining the response to CCR5 receptor antagonists. Initially, phenotypic tests were used, but limited access to these tests prompted the development of alternative strategies. Recently, genotyping tropism has been validated using a Canadian technology in clinical trials investigating the use of maraviroc in both experienced and treatment-naive patients. The present guidelines review the evidence supporting the use of genotypic assays and provide recommendations regarding tropism testing in daily clinical management.

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

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

Clinicopathological relevance of BRAF mutations in human cancer

BRAF represents one of the most frequently mutated protein kinase genes in human tumours. The mutation is commonly tested in pathology practice. BRAF mutation is seen in melanoma, papillary thyroid carcinoma (including papillary thyroid carcinoma arising from ovarian teratoma), ovarian serous tumours, colorectal carcinoma, gliomas, hepatobiliary carcinomas and hairy cell leukaemia. In these cancers, various genetic aberrations of the BRAF proto-oncogene, such as different point mutations and chromosomal rearrangements, have been reported. The most common mutation, BRAF V600E, can be detected by DNA sequencing and immunohistochemistry on formalin fixed, paraffin embedded tumour tissue. Detection of BRAF V600E mutation has the potential for clinical use as a diagnostic and prognostic marker. In addition, a great deal of research effort has been spent in strategies inhibiting its activity. Indeed, recent clinical trials involving BRAF selective inhibitors exhibited promising response rates in metastatic melanoma patients. Clinical trials are underway for other cancers. However, cutaneous side effects of treatment have been reported and therapeutic response to cancer is short-lived due to the emergence of several resistance mechanisms. In this review, we give an update on the clinical pathological relevance of BRAF mutation in cancer. It is hoped that the review will enhance the direction of future research and assist in more effective use of the knowledge of BRAF mutation in clinical practice.

Incidence and prevalence of intrasubtype HIV-1 dual infection in at-risk men in the United States

BACKGROUND

Human immunodeficiency virus type 1 (HIV-1) dual infection (DI) has been associated with decreased CD4 T-cell counts and increased viral loads; however, the frequency of intrasubtype DI is poorly understood. We used ultradeep sequencing (UDS) to estimate the frequency of DI in a primary infection cohort of predominantly men who have sex with men (MSM).

METHODS

 HIV-1 genomes from longitudinal blood samples of recently infected, therapy-naive participants were interrogated with UDS. DI was confirmed when maximum sequence divergence was excessive and supported by phylogenetic analysis. Coinfection was defined as DI at baseline; superinfection was monoinfection at baseline and DI at a later time point.

RESULTS

 Of 118 participants, 7 were coinfected and 10 acquired superinfection. Superinfection incidence rate was 4.96 per 100 person-years (95% confidence interval [CI], 2.67-9.22); 6 occurred in the first year and 4 in the second. Overall cumulative prevalence of intrasubtype B DI was 14.4% (95% CI, 8.6%-22.1%). Primary HIV-1 incidence was 4.37 per 100 person-years (95% CI, 3.56-5.36).

CONCLUSIONS

 Intrasubtype DI was frequent and comparable to primary infection rates among MSM in San Diego; however, superinfection rates declined over time. DI is likely an important component of the HIV epidemic dynamics, and development of stronger immune responses to the initial infection may protect from superinfection.

Genotypic analysis of the V3 region of HIV from virologic nonresponders to maraviroc-containing regimens reveals distinct patterns of failure

Changes in HIV tropism from R5 to non-R5 or development of drug resistance is often associated with virologic failure in patients treated with maraviroc, a CCR5 antagonist. We sought to examine changes in HIV envelope sequences and inferred tropism in patients who did not respond to maraviroc-based regimens. We selected 181 patients who experienced early virologic failure on maraviroc-containing therapy in the MOTIVATE trials. All patients had R5 HIV by the original Trofile assay before entry. We used population-based sequencing methods and the geno2pheno algorithm to examine changes in tropism and V3 sequences at the time of failure. Using deep sequencing, we assessed whether V3 sequences observed at failure emerged from preexisting subpopulations. From population genotyping data at failure, 90 patients had R5 results, and 91 had non-R5 results. Of the latter group, the geno2pheno false-positive rate (FPR) value fell from a median of 20 at screening to 1.1 at failure. By deep sequencing, the median percentage of non-R5 variants in these patients rose from 1.4% to 99.5% after a median of 4 weeks on maraviroc. In 70% of cases, deep sequencing could detect a pretreatment CXCR4-using subpopulation, which emerged at failure. Overall, there were two distinct patterns of failure of maraviroc. Patients failing with R5 generally had few V3 substitutions and low non-R5 prevalence by deep sequencing. Patients with non-R5 HIV who were failing developed very-high-prevalence non-R5 HIV (median, 99%) and had very low geno2pheno values.

Entry inhibitors and their use in the treatment of HIV-1 infection

Entry of HIV into target cells is a complex, multi-stage process involving sequential attachment and CD4 binding, coreceptor binding, and membrane fusion. HIV entry inhibitors are a complex group of drugs with multiple mechanisms of action depending on the stage of the viral entry process they target. Two entry inhibitors are currently approved for the treatment of HIV-infected patients. Maraviroc, a CCR5 antagonist, blocks interactions between the viral envelope proteins and the CCR5 coreceptor. Enfuvirtide, a fusion inhibitor, disrupts conformational changes in gp41 that drive membrane fusion. A wide array of additional agents are in various stages of development. This review covers the entry inhibitors and their use in the treatment of HIV-infected patients.

Sensitive cell-based assay for determination of human immunodeficiency virus type 1 coreceptor tropism

CCR5 antagonists are a powerful new class of antiretroviral drugs that require a companion assay to evaluate the presence of CXCR4-tropic (non-R5) viruses prior to use in human immunodeficiency virus (HIV)-infected individuals. In this study, we have developed, characterized, verified, and prevalidated a novel phenotypic test to determine HIV-1 coreceptor tropism (VERITROP) based on a sensitive cell-to-cell fusion assay. A proprietary vector was constructed containing a near-full-length HIV-1 genome with the yeast uracil biosynthesis (URA3) gene replacing the HIV-1 env coding sequence. Patient-derived HIV-1 PCR products were introduced by homologous recombination using an innovative yeast-based cloning strategy. The env-expressing vectors were then used in a cell-to-cell fusion assay to determine the presence of R5 and/or non-R5 HIV-1 variants within the viral population. Results were compared with (i) the original version of Trofile (Monogram Biosciences, San Francisco, CA), (ii) population sequencing, and (iii) 454 pyrosequencing, with the genotypic data analyzed using several bioinformatics tools, i.e., the 11/24/25 rule, Geno2Pheno (2% to 5.75%, 3.5%, or 10% false-positive rate [FPR]), and webPSSM. VERITROP consistently detected minority non-R5 variants from clinical specimens, with an analytical sensitivity of 0.3%, with viral loads of ≥1,000 copies/ml, and from B and non-B subtypes. In a pilot study, a 73.7% (56/76) concordance was observed with the original Trofile assay, with 19 of the 20 discordant results corresponding to non-R5 variants detected using VERITROP and not by the original Trofile assay. The degree of concordance of VERITROP and Trofile with population and deep sequencing results depended on the algorithm used to determine HIV-1 coreceptor tropism. Overall, VERITROP showed better concordance with deep sequencing/Geno2Pheno at a 0.3% detection threshold (67%), whereas Trofile matched better with population sequencing (79%). However, 454 sequencing using Geno2Pheno at a 10% FPR and 0.3% threshold and VERITROP more accurately predicted the success of a maraviroc-based regimen. In conclusion, VERITROP may promote the development of new HIV coreceptor antagonists and aid in the treatment and management of HIV-infected individuals prior to and/or during treatment with this class of drugs.

Clinical implications of HIV-1 minority variants

Technologic advances in human immunodeficiency virus type 1 (HIV-1) sequencing have revolutionized the study of antiretroviral drug resistance and are increasingly moving from the laboratory to clinical practice. These techniques are able to detect HIV-1 drug resistance mutations present at low frequencies not detectable by current HIV-1 genotyping assays. For a number of commonly used antiretroviral medications, such as nonnucleoside reverse transcriptase inhibitors, the detection of these drug-resistant minority variants significantly increases the risk of treatment failure. The level of evidence, however, is insufficient to determine the impact of HIV-1 minority variants for several other classes of antiretroviral medications. Clinicians should be aware of the novel technologies that are moving into routine clinical use and the clinical implications of HIV-1 minority variants. Additional studies are needed to determine the optimal platform for clinical application of these new technologies and to provide guidance to clinicians on the type and frequency of clinically important HIV-1 minority variants.