Concordance between two phenotypic assays and ultradeep pyrosequencing for determining HIV-1 tropism

A Review of Next Generation Sequencing Methods and its Applications in Laboratory Diagnosis

Next-generation sequencing (NGS) is a new technology used to detect the sequence of DNA and RNA and to detect mutations or variations of significance. NGS generates large quantities of sequence data within a short time duration. The various types of sequencing includes Sanger Sequencing, Pyrosequencing, Sequencing by Synthesis (Illumina), Ligation (SoLID), Single molecule Fluorescent Sequencing (Helicos), Single molecule Real time Sequencing (Pacbio), Semiconductor sequencing (Ion torrent technology), Nanopore sequencing and fourth generation sequencing. These methods of sequencing have been modified and improved over the years such that it has become cost effective and accessible to diagnostic laboratories. Management of Outbreaks, rapid identification of bacteria, molecular case finding, taxonomy, detection of the zoonotic agents and guiding prevention strategies in HIV outbreaks are just a few of the many applications of Next Generation sequencing in clinical microbiology.

A highly sensitive method for the detection of recombinant PERV-A/C env RNA using next generation sequencing technologies

Several xenogenic cell-based therapeutic products are currently under development around the world for the treatment of human diseases. Porcine islet cell products for treating human diabetes are a typical example. Since porcine cells possess endogenous retrovirus (PERV), which can replicate in human cells in vitro, the potential transmission of PERV has raised concerns in the development of these products. Four subgroups of infectious PERV have been identified, namely PERV-A, -B, -C, and recombinant PERV-A/C. Among them, PERV-A/C shows a high titre and there was a paper reported that an incidence of PERV-A/C viremia was increased in diseased pigs; thus, it would be important to monitor the emergence of PERV-A/C after transplantation of porcine products. In this study, we developed a highly sensitive method for the detection of PERV-A/C using next generation sequencing (NGS) technologies. A model PERV-C spiked with various doses of PERV-A/C were amplified by RT-PCR and the amplicons were analysed by NGS. We found that the NGS analysis allowed the detection of PERV-A/C at the abundance ratios of 1% and 0.1% with true positive rates of 100% and 57%, respectively, indicating that it would be useful for the rapid detection of PERV-A/C emergence after transplantation of porcine products.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

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

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

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.

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.

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.

Frequent detection of CXCR4-using viruses among Brazilian blood donors with HIV-1 long-standing infection and unknown clinical stage: Analysis of massive parallel sequencing data

The determination of viral tropism is critically important and highly recommended to guide therapy with the CCR5 antagonist, which does not inhibit the effect of X4-tropic viruses. Here, we report the prevalence of HIV-1×4 HIV strains in 84 proviral DNA massively parallel sequencing “MPS” data from well-defined non-recently infected first-time Brazilian blood donors. The MPS data covering the entire V3 region of the env gene was extracted from our recently generated HIV-1 genomes sequenced by a paired-end protocol (Illumina). Of the 84 MPS data samples, 63 (75%) were derived from donors with long-standing infection and 21 (25%) were lacking stage information. HIV‐1 tropism was inferred using Geno2pheno (g2p) _[454] algorithm (FPR=1%, 2.5%, and 3.75%). Among the 84 data samples for which tropism was defined by g2p_2.5%, 13 (15.5%) participants had detectable CXCR4-using viruses in their MPS reads. Mixed infections with R5 and X4 were observed in 11.9% of the study subjects and minority X4 viruses were detected in 7 (8.3%) of participants. Nine of the 63 (14.3%) subjects with LS infection were predicted by g2p _2.5% to harbor proviral CXCR4-using viruses. Our findings of a high proportion of blood donors (15.5%) harboring CXCR4-using viruses in PBMCs may indicate that this phenomenon is common. These findings may have implications for clinical and therapeutic aspects and may benefit individuals who plan to receive CCR5 antagonists.

Position-specific automated processing of V3 env ultra-deep pyrosequencing data for predicting HIV-1 tropism

HIV-1 coreceptor usage must be accurately determined before starting CCR5 antagonist-based treatment as the presence of undetected minor CXCR4-using variants can cause subsequent virological failure. Ultra-deep pyrosequencing of HIV-1 V3 env allows to detect low levels of CXCR4-using variants that current genotypic approaches miss. However, the computation of the mass of sequence data and the need to identify true minor variants while excluding artifactual sequences generated during amplification and ultra-deep pyrosequencing is rate-limiting. Arbitrary fixed cut-offs below which minor variants are discarded are currently used but the errors generated during ultra-deep pyrosequencing are sequence-dependant rather than random. We have developed an automated processing of HIV-1 V3 env ultra-deep pyrosequencing data that uses biological filters to discard artifactual or non-functional V3 sequences followed by statistical filters to determine position-specific sensitivity thresholds, rather than arbitrary fixed cut-offs. It allows to retain authentic sequences with point mutations at V3 positions of interest and discard artifactual ones with accurate sensitivity thresholds.

HIV-1 Coreceptor Usage Assessment by Ultra-Deep Pyrosequencing and Response to Maraviroc

Background

Maraviroc is an HIV entry inhibitor that alters the conformation of CCR5 and is poorly efficient in patients infected by viruses that use CXCR4 as an entry coreceptor. The goal of this study was to assess the capacity of ultra-deep pyrosequencing (UDPS) and different data analysis approaches to characterize HIV tropism at baseline and predict the therapeutic outcome on maraviroc treatment.

Methods

113 patients with detectable HIV-1 RNA on HAART were treated with maraviroc. The virological response was assessed at months 1, 3 and 6. The sequence of the HIV V3 loop was determined at baseline and prediction of maraviroc response by different software and interpretation algorithms was analyzed.

Results

UDPS followed by analysis with the Pyrotrop software or geno2pheno algorithm provided better prediction of the response to maraviroc than Sanger sequencing. We also found that the H34Y/S substitution in the V3 loop was the strongest individual predictor of maraviroc response, stronger than substitutions at positions 11 or 25 classically used in interpretation algorithms.

Conclusions

UDPS is a powerful tool that can be used with confidence to predict maraviroc response in HIV-1-infected patients. Improvement of the predictive value of interpretation algorithms is possible and our results suggest that adding the H34S/Y substitution would substantially improve the performance of the 11/25/charge rule.

A diagnostic HIV-1 tropism system based on sequence relatedness

Key clinical studies for HIV coreceptor antagonists have used the phenotyping-based Trofile test. Meanwhile various simpler-to-do genotypic tests have become available that are compatible with standard laboratory equipment and Web-based interpretation tools. However, these systems typically analyze only the most prominent virus sequence in a specimen. We present a new diagnostic HIV tropism test not needing DNA sequencing. The system, XTrack, uses physical properties of DNA duplexes after hybridization of single-stranded HIV-1 env V3 loop probes to the clinical specimen. Resulting "heteroduplexes" possess unique properties driven by sequence relatedness to the reference and resulting in a discrete electrophoretic mobility. A detailed optimization process identified diagnostic probe candidates relating best to a large number of HIV-1 sequences with known tropism. From over 500 V3 sequences representing all main HIV-1 subtypes (Los Alamos database), we obtained a small set of probes to determine the tropism in clinical samples. We found a high concordance with the commercial TrofileES test (84.9%) and the Web-based tool Geno2Pheno (83.0%). Moreover, the new system reveals mixed virus populations, and it was successful on specimens with low virus loads or on provirus from leukocytes. A replicative phenotyping system was used for validation. Our data show that the XTrack test is favorably suitable for routine diagnostics. It detects and dissects mixed virus populations and viral minorities; samples with viral loads (VL) of <200 copies/ml are successfully analyzed. We further expect that the principles of the platform can be adapted also to other sequence-divergent pathogens, such as hepatitis B and C viruses.

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.

Evolution of HIV-1 quasispecies and coreceptor use in cell reservoirs of patients on suppressive antiretroviral therapy

OBJECTIVES

To track changes in the V3 env region of HIV-1 quasispecies and determine virus coreceptor use in cell reservoirs of patients on long-term suppressive antiretroviral therapy (ART).

PATIENTS AND METHODS

Ten patients whose plasma viraemia had been suppressed for a median of 5.5 years were followed for 5 years. The V3 env regions of viruses in peripheral blood mononuclear cells were analysed by ultra-deep sequencing (UDS). HIV-1 tropism was predicted using the geno2pheno 5.75 algorithm and a phenotypic assay.

RESULTS

The UDS and phenotypic assay data were concordant for predicting HIV-1 tropism. CXCR4-using viruses detected by UDS accounted for 14.7%-76.5% of the virus populations in samples from five patients at enrolment. Five patients harboured pure R5 virus populations and no X4 viruses emerged during the 5 years. The selection pressures estimated by the dN/dS ratio were acting on the V3 region to produce diversification of the quasispecies in CXCR4-infected patients and purification of the quasispecies in R5-infected patients on effective ART.

CONCLUSIONS

UDS showed that the virus quasispecies in cell reservoirs of patients on long-term suppressive ART continued to evolve. CXCR4-using variants became more diversified. Analysis of the selection pressures on the virus quasispecies could provide a clearer picture of virus persistence in patients on effective ART.

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.

Algorithm-based prediction of HIV-1 subtype D coreceptor use

We compared the coreceptor tropism-predicting performance of a specific genotypic algorithm for HIV-1 subtype D and that of the geno2pheno algorithm with different cutoffs. The D-specific algorithm and geno2pheno with a false-positivity rate cutoff of 2.5% had the same concordance with the phenotypic determination. The geno2pheno algorithm with a false-positivity rate cutoff of 2.5%, more sensitive but slightly less specific, seems to be an appropriate alternative.

Current assays for HIV-1 diagnostics and antiretroviral therapy monitoring: challenges and possibilities

In 2011, there were over 34 million people living with HIV infections, placing a heavy burden on public health sectors. HIV infection is a lifelong threat that cannot be prevented by vaccination or cured by antiretroviral drugs. The infected patients rely on daily antiretroviral therapy to suppress HIV viral replication. Hence, it is important to diagnose HIV infections as early as possible and to monitor the efficacy of antiretroviral therapy every 3–6 months. Different immunoassays detecting HIV antigens and antibodies have been modified to offer better sensitivity and more rapid diagnosis. Several clinical and virological parameters, including CD4+ cell counts, viral load and drug resistance mutations, are also used for treatment monitoring. Many molecular assay optimizations are now being utilized to improve patient care. This review will focus on the most updated HIV diagnostic assays, as well as discussing the upcoming possibilities of other advanced technologies.

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.

The genotypic false positive rate determined by V3 population sequencing can predict the burden of HIV-1 CXCR4-using species detected by pyrosequencing

Objective

The false-positive rate (FPR) is a percentage-score provided by Geno2Pheno-algorithm indicating the likelihood that a V3-sequence is falsely predicted as CXCR4-using. We evaluated the correlation between FPR obtained by V3 population-sequencing and the burden of CXCR4-using variants detected by V3 ultra-deep sequencing (UDPS) and Enhanced-Sensitivity Trofile assay (ESTA).

Methods

54 HIV-1 B-subtype infected-patients (all maraviroc-naïve), with viremia >10,000copies/ml, were analyzed. HIV-tropism was assessed by V3 population-sequencing, UDPS (considering variants with >0.5% prevalence), and ESTA.

Results

By UDPS, CCR5-using variants were detected in 53/54 patients, irrespective of FPR values, and their intra-patient prevalence progressively increased by increasing the FPR obtained by V3 population-sequencing (rho = 0.75, p = 5.0e-8). Conversely, the intra-patient prevalence of CXCR4-using variants in the 54 patients analyzed progressively decreased by increasing the FPR (rho = −0.61; p = 9.3e-6). Indeed, no CXCR4-using variants were detected in 13/13 patients with FPR>60. They were present in 7/18 (38.8%) patients with FPR 20–60 (intra-patient prevalence range: 2.1%–18.4%), in 5/7 (71.4%) with FPR 10–20, in 4/6 (66.7%) with FPR 5–10, and in 10/10(100%) with FPR<5 (intra-patient prevalence range: 12.1%–98.1%).

Conclusions

FPR by V3 population-sequencing can predict the burden of CXCR4-using variants. This information can be used to optimize the management of tropism determination in clinical practice. Due to its low cost and short turnaround time, V3 population-sequencing may represent the most feasible test for HIV-1 tropism determination. More sensitive methodologies (as UDPS) might be useful when V3 population-sequencing provides a FPR >20 (particularly in the range 20–60), allowing a more careful identification of patients harboring CXCR4-using variants.

Next-generation sequencing to assess HIV tropism

PURPOSE OF REVIEW

Deep sequencing of the V3 region of the HIV envelope gene can detect minority non-R5 variants in patients with high sensitivity and specificity. As next-generation sequencing approaches have matured, the clinical utility of deep sequencing for HIV tropism has entered the clinic. Accurate and sensitive tropism testing is essential for successful treatment with the CCR5 antagonist class of antiretrovirals.

RECENT FINDINGS

This review will focus on five aspects of next-generation sequencing for assessing HIV tropism: some background on the necessity of deep sequencing versus other tropism methods; the methodological process of 454 sequencing and analysis; other next-generation sequencing technologies; the diagnostic performance of deep sequencing relative to other tropism assays; and the use of deep sequencing in clinical practice.

SUMMARY

This method has emerged quickly as both a research and clinical tool because of its high concordance with commonly used phenotypic tropism assays and its ability to predict virological response to CCR5 antagonist-containing regimens.

Phenotyping methods for determining HIV tropism and applications in clinical settings

PURPOSE OF REVIEW

HIV-1 enters CD4-expressing cells via one or both of the chemokine receptors CCR5 and CXCR4. Specific CCR5 antagonists are now in clinical use, but only for CCR5-tropic viruses. Hence, several methods have been developed for assessing HIV-1 tropism in patients who are candidates for CCR5 antagonists. This article reviews current data on phenotypic assays of tropism.

RECENT FINDINGS

Phenotypic assays are still used as reference, although genotypic methods have improved. The main advantages of phenotypic assays are their great sensitivity for detecting minor CXCR4-using variants and their capacity to assess non-B subtypes of HIV-1. Clinical trials of maraviroc have, thus, relied on the phenotypic determination of HIV-1 tropism. However, new genotypic approaches that are more sensitive for minor CXCR4-using variants, notably ultra-deep pyrosequencing, are now challenging phenotypic assays. Nevertheless, phenotypic assays are essential for improving genotypic algorithms for determining HIV-1 tropism as well as for assessing the resistance of R5-tropic viruses to CCR5 antagonists.

SUMMARY

HIV-1 tropism should be determined before using CCR5 antagonists. Phenotypic recombinant assays are still the benchmark tests for characterizing HIV-1 tropism as their great sensitivity enables them to detect minor CXCR4-using variants of both B and non-B HIV-1 subtypes.

HIV-1 dynamics and coreceptor usage in Maraviroc-treated patients with ongoing replication

There is evidence that HIV-1 evolution under maraviroc (MVC) pressure can lead to the selection of either X4-tropic variants and/or R5-tropic, MVC-resistant isolates. However, the viral dynamics of HIV-1 variants in patients with virological failure (VF) on MVC-containing regimens remain poorly studied. Here, we investigated the V3 loop evolution of HIV-1 on MVC in relation to coreceptor usage and the nature of HIV-1 quasispecies before MVC therapy using bulk population sequences and ultradeep sequencing. The majority of patients had no detectable minority X4 variant at baseline. The evolution of tropism was followed up until VF and showed three possibilities for viral evolution in these patients: emergence of preexisting X4 variants, de novo selection of R5 variants presenting V3 loop mutations, or replication of R5 variants without selection of known mutations.

Correlation between genotypic (V3 population sequencing) and phenotypic (Trofile ES) methods of characterizing co-receptor usage of HIV-1 from 200 treatment-naïve HIV patients screened for Study A4001078

Assessment of HIV-1 co-receptor usage is essential to identify patients who are likely to respond to maraviroc (MVC)-containing regimens. Co-receptor usage of plasma virus from all treatment-naïve patients screened for a MVC clinical trial was assessed using phenotypic and genotypic methodologies to evaluate concordance between testing methods and to assess the quantity of CXCR4-using (non-R5) virus in samples giving discordant results. Co-receptor usage was prospectively measured using the enhanced sensitivity Trofile assay (Trofile ES) to screen patients for enrollment in Study A4001078. Population and deep sequencing methodologies were utilized retrospectively to analyze all screening samples, with co-receptor usage determined using the geno2pheno algorithm. Concordance between methods was explored using descriptive statistics. The quantity of non-R5 virus in all samples was measured using deep sequencing. Trofile ES and matched genotype results were obtained for 199screening samples. Concordance of Trofile ES with population genotyping (5.75% false-positive rate [FPR]) and deep sequencing (3.5% FPR; 2% non-R5 threshold) was 91.7% and 89.6%, respectively. Population genotype data were available for all samples with non-reportable Trofile ES results; the distribution of co-receptor usage in this set was consistent with that in the overall population: 75% (12/16) R5 and 25% (4/16) non-R5. The majority of samples contained non-R5 plasma HIV-1 RNA estimated at either <1 log(10) (62.0%) or ⩾4 log(10) (30.5%) copies/mL; the absolute amount of detectable non-R5 virus remained stable between screening and baseline visits. Samples originally classified as non-R5 by Trofile ES but R5 by population sequencing had a relatively low absolute amount of non-R5 virus (mean 2.1%, median 0.1%). The determination of co-receptor usage using either Trofile ES or genotyping methodologies showed similar frequencies of R5 and non-R5 virus in this treatment-naïve study population. For both concordant and discordant samples, population sequencing appropriately identified R5 samples with low levels of non-R5-using virus.

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

HIV-1 coreceptor tropism assays are required to rule out the presence of CXCR4-tropic (non-R5) viruses prior treatment with CCR5 antagonists. Phenotypic (e.g., Trofile™, Monogram Biosciences) and genotypic (e.g., population sequencing linked to bioinformatic algorithms) assays are the most widely used. Although several next-generation sequencing (NGS) platforms are available, to date all published deep sequencing HIV-1 tropism studies have used the 454™ Life Sciences/Roche platform. In this study, HIV-1 co-receptor usage was predicted for twelve patients scheduled to start a maraviroc-based antiretroviral regimen. The V3 region of the HIV-1 env gene was sequenced using four NGS platforms: 454™, PacBio® RS (Pacific Biosciences), Illumina®, and Ion Torrent™ (Life Technologies). Cross-platform variation was evaluated, including number of reads, read length and error rates. HIV-1 tropism was inferred using Geno2Pheno, Web PSSM, and the 11/24/25 rule and compared with Trofile™ and virologic response to antiretroviral therapy. Error rates related to insertions/deletions (indels) and nucleotide substitutions introduced by the four NGS platforms were low compared to the actual HIV-1 sequence variation. Each platform detected all major virus variants within the HIV-1 population with similar frequencies. Identification of non-R5 viruses was comparable among the four platforms, with minor differences attributable to the algorithms used to infer HIV-1 tropism. All NGS platforms showed similar concordance with virologic response to the maraviroc-based regimen (75% to 80% range depending on the algorithm used), compared to Trofile (80%) and population sequencing (70%). In conclusion, all four NGS platforms were able to detect minority non-R5 variants at comparable levels suggesting that any NGS-based method can be used to predict HIV-1 coreceptor usage.

Challenges and opportunities in estimating viral genetic diversity from next-generation sequencing data

Many viruses, including the clinically relevant RNA viruses HIV (human immunodeficiency virus) and HCV (hepatitis C virus), exist in large populations and display high genetic heterogeneity within and between infected hosts. Assessing intra-patient viral genetic diversity is essential for understanding the evolutionary dynamics of viruses, for designing effective vaccines, and for the success of antiviral therapy. Next-generation sequencing (NGS) technologies allow the rapid and cost-effective acquisition of thousands to millions of short DNA sequences from a single sample. However, this approach entails several challenges in experimental design and computational data analysis. Here, we review the entire process of inferring viral diversity from sample collection to computing measures of genetic diversity. We discuss sample preparation, including reverse transcription and amplification, and the effect of experimental conditions on diversity estimates due to in vitro base substitutions, insertions, deletions, and recombination. The use of different NGS platforms and their sequencing error profiles are compared in the context of various applications of diversity estimation, ranging from the detection of single nucleotide variants (SNVs) to the reconstruction of whole-genome haplotypes. We describe the statistical and computational challenges arising from these technical artifacts, and we review existing approaches, including available software, for their solution. Finally, we discuss open problems, and highlight successful biomedical applications and potential future clinical use of NGS to estimate viral diversity.

[Methods for determination of HIV tropism and their clinical use]

Determination of HIV-1 tropism is mandatory before using CCR5 antagonists in clinical practice. One drug of this class, maraviroc, has been approved for the treatment of HIV infection. The phenotypic assay, TrofileTM, was clinically validated in the clinical development program of maraviroc and has been widely used to select candidates for maraviroc therapy. Phenotypic tests, however, have the disadvantage of being complex, are costly and time-consuming, and their accessibility is limited, which hampers their routine use in clinical diagnosis. Genotypic assays, based on sequencing the third hypervariable (V3 loop) of the viral gene env, interpreted according to various genotypic bioinformatic tools, such as geno2pheno and PSSM, are faster and cheaper than phenotypic assays, and are also more accessible. In retrospective analyses of the maraviroc pivotal trials, genotypic methods using either conventional ("bulk") or deep-sequencing technology predicted virologic response to maraviroc similarly to phenotypic assays and are now included within several European recommendations to guide the clinical use of CCR5 antagonists.

Discordance in HIV-1 co-receptor use prediction by different genotypic algorithms and phenotype assay: intermediate profile in relation to concordant predictions

Concordant and discordant genotypic predictions of HIV-1 co-receptor tropism were analyzed. V3 region was sequenced from plasma samples of patients screened for R5 tropism by the Trofile® assay, before CCR5 antagonist prescription. Ten tools including geno2pheno, PSSM, an "11/25" and "net charge" rule, and other published algorithms were used. Patients were grouped according to concordance or discordance between tools and Trofile® result. Trofile® tropism reports from 50 patient samples were R5 in 38 and Dual/Mixed (DM) in 12. Prediction with the genotypic tools were concordant for 23 R5 samples, and discordant for the 15 other ones. From Trofile® DM strains were concordant in 6 and discordant in 6. V3 sequences were not clearly distinct between R5 and DM strains, except a greater diversity in the later. Discordances were found with any tool or combination of them, so that no one can be proposed as better than the others. Predictive values of each algorithm were similar and rather good (efficacy ranged from 74% to 84%), but the rate of non-confirmed prediction is greater when compelling the results of all tools with each individual sample. The mean of quantitative values obtained with one tool when another tool give the opposite prediction were different from those obtained when all tools agree with that prediction. The two discordant groups were often not distinguishable from each other. These results suggest that viruses giving discordant prediction with bioinformatic tools could be functionally distinct and/or in a different evolutionary state compared to those with concordant prediction.

Geno2pheno[454]: a Web server for the prediction of HIV-1 coreceptor usage from next-generation sequencing data

Inferring HIV-1 coreceptor usage from a genotype is becoming more and more important for the appropriate treatment of long-term patients. While results are already encouraging where standard bulk-nucleic acid sequencing methods are used, they are limited with respect to the detection of minor variants. In contrast, next-generation sequencing methods (ultradeep sequencing, pyrosequencing) are capable of sequencing virus quasispecies at very low quantities. However, as well as being very expensive, these methods generate vast amounts of data such that sequence analysis has to be automated by computer assistance. Here, we describe the geno2pheno[454] system which handles all processing and prediction steps involved in the prediction of coreceptor usage from massively parallel sequencing data. The system is split into a JAVA preprocessor which is run locally on the client side and a Web server which generates the prediction results. Predictions are based on the same prediction method as used in the geno2pheno[coreceptor] tool.

High-throughput sequencing of PCR products tagged with universal primers using 454 life sciences systems

This unit describes a method to convert PCR products (amplicons) flanked by universal M13 primers into a library for use on all 454 Sequencing Systems (454 Life Sciences, a Roche Company). This is especially useful for simultaneous sequencing and analysis of large numbers of amplicons or for the detection of minor variations within the amplified products. The method described here involves preparing a library of DNA with specific primers containing adaptor sequences recognized by the GS Junior System sequencing process. The data from the sequencing run are processed and analyzed by 454 Life Sciences software. This approach allows for multiplexing of a large number of amplicons to streamline processing and analysis. Any pre-existing universally tagged amplicon, primer set, or plasmid with M13 sequences flanking the cloning site can be used in this protocol.

Drug safety evaluation of maraviroc for the treatment of HIV infection

INTRODUCTION

Maraviroc is the only C-chemokine receptor 5 (CCR5) antagonist approved for the treatment of infection with HIV. This article reviews the safety and efficacy of maraviroc in the treatment of HIV infection.

AREAS COVERED

The PubMed database was searched using the keywords 'maraviroc' and 'HIV'. In addition, conference proceedings from CROI, IAS and EACS meetings were searched for maraviroc clinical trials. The PubMed search revealed one Phase IIb - III clinical trial in treatment-naive HIV(+) patients (MERIT) and three Phase IIb - III randomized clinical trials (RCTs) in treatment-experienced patients (MOTIVATE 1 and 2, A4001029). All RCTs showed an excellent safety profile for maraviroc in the treatment of HIV-1 infection. However, long-term (> 3 years) safety data generated on maraviroc therapy are still scarce. Based on the findings from RCTs so far, no relevant toxicities and co-morbidities such as coronary heart disease or hepatotoxicity have been described. The overall CD4(+) cell count increase resulting from a maraviroc-containing regimen appears to be higher than those seen with other antiretroviral regimens. However, the significance remains controversial. To date, maraviroc has shown a potent and durable virological efficacy profile for the treatment of HIV-1 infection. The only use of maraviroc depends on pretreatment testing for CCR5 tropism.

EXPERT OPINION

Maraviroc is a generally safe and well-tolerated medication for the treatment of HIV-1 infection with a unique mechanism of action. Long-term (i.e., > 5 years) risks are not known and have to be carefully monitored.

RNA virus population diversity: implications for inter-species transmission

RNA viruses are notorious for rapidly generating genetically diverse populations during a single replication cycle, and the implications of this mutant population, often referred to as quasispecies, can be vast. Previous studies have linked RNA virus genetic variability to changes in viral pathogenesis, the ability to adapt to a host during infection, and to the acquisition of mechanisms required to switch hosts entirely. However, these initial studies are just the beginning. With the development of next generation technologies, groups will be able to dig deeper into the sequence space that is generated during an RNA virus infection and more clearly understand the development, role, and consequences of viral genetic diversity.

Frequency of CXCR4-using viruses in primary HIV-1 infections using ultra-deep pyrosequencing

We used ultra-deep pyrosequencing and the Toulouse Tropism Test phenotypic assay to determine the prevalence of CXCR4-using viruses in 21 patients with primary HIV-1 infections. We found X4-containing virus populations in 9% of patients by ultra-deep pyrosequencing using position-specific scoring matrices (PSSM(X4/R5)) or geno2pheno(5.75) and in 14% using the combined 11/25 and net charge rule. The phenotypic assay identified 9% of CXCR4-using viruses. This confirms that R5 viruses are predominant in primary HIV-1 infections.

Genotypic prediction of HIV-1 subtype D tropism

Background

HIV-1 subtype D infections have been associated with rapid disease progression and phenotypic assays have shown that CXCR4-using viruses are very prevalent. Recent studies indicate that the genotypic algorithms used routinely to assess HIV-1 tropism may lack accuracy for non-B subtypes. Little is known about the genotypic determinants of HIV-1 subtype D tropism.

Results

We determined the HIV-1 coreceptor usage for 32 patients infected with subtype D by both a recombinant virus phenotypic entry assay and V3-loop sequencing to determine the correlation between them. The sensitivity of the Geno2pheno10 genotypic algorithm was 75% and that of the combined 11/25 and net charge rule was 100% for predicting subtype D CXCR4 usage, but their specificities were poor (54% and 68%). We have identified subtype D determinants in the V3 region associated with CXCR4 use and built a new simple genotypic rule for optimizing the genotypic prediction of HIV-1 subtype D tropism. We validated this algorithm using an independent GenBank data set of 67 subtype D V3 sequences of viruses of known phenotype. The subtype D genotypic algorithm was 68% sensitive and 95% specific for predicting X4 viruses in this data set, approaching the performance of genotypic prediction of HIV-1 subtype B tropism.

Conclusion

The genotypic determinants of HIV-1 subtype D coreceptor usage are slightly different from those for subtype B viruses. Genotypic predictions based on a subtype D-specific algorithm appear to be preferable for characterizing coreceptor usage in epidemiological and pathophysiological studies.