Improved detection of CXCR4-using HIV by V3 genotyping: application of population-based and "deep" sequencing to plasma RNA and proviral DNA

Correlation of the virological response to short-term maraviroc monotherapy with standard and deep-sequencing-based genotypic tropism prediction methods

Genotypic tropism testing methods are emerging as the first step before prescription of the CCR5 antagonist maraviroc (MVC) to HIV-infected patients in Europe. Studies validating genotypic tests have included other active drugs that could have potentially convoluted the effects of MVC. The maraviroc clinical test (MCT) is an in vivo drug sensitivity test based on the virological response to a short-term exposure to MVC monotherapy. Thus, our aim was to compare the results of genotypic tropism testing methods with the short-term virological response to MVC monotherapy. A virological response in the MCT was defined as a ≥ 1-log(10) decrease in HIV RNA or undetectability after 8 days of drug exposure. Seventy-three patients undergoing the MCT were included in this study. We used both standard genotypic methods (n = 73) and deep sequencing (n = 27) on MCT samples at baseline. For the standard methods, the most widely used genotypic algorithms for analyzing the V3 loop sequence, geno2pheno and PSSM, were used. For deep sequencing, the geno2pheno algorithm was used with a false-positive rate cutoff of 3.5. The discordance rates between the standard genotypic methods and the virological response were approximately 20% (including mostly patients without a virological response). Interestingly, these discordance rates were similar to that obtained from deep sequencing (18.5%). The discordance rates between the genotypic methods (tropism assays predictive of the use of the CCR5 coreceptor) and the MCT (in vivo MVC sensitivity assay) indicate that the algorithms used by genotypic methods are still not sufficiently optimized.

Dates of HIV infection can be estimated for seroprevalent patients by coalescent analysis of serial next-generation sequencing data

OBJECTIVE

To reconstruct dates of HIV infection by the coalescent analysis of longitudinal next-generation sequencing (NGS) data.

DESIGN

The coalescent predicts the time that has elapsed since the most recent common ancestor (MRCA) of a population. Because HIV tends to undergo severe bottlenecks upon transmission, the MRCA may be a good predictor of the time of infection. NGS provides an efficient means for performing large-scale clonal sequencing of HIV populations within patients, and the ideal raw material for coalescent analysis.

METHODS

Baseline and follow-up plasma samples were obtained from 19 individuals enrolled into the Montréal Primary HIV Infection cohort. Dates of infection were initially estimated at baseline from nongenetic data (clinical and serological markers and patient questionnaires). HIV RNA was extracted and seven regions of the genome were amplified, subjected to parallel-tagged 454 pyrosequencing, and analyzed using the software package BEAST.

RESULTS

Mean estimates of the time to the MRCA per patient were significantly correlated with nongenetic estimates (Spearman's ρ = 0.65, P = 4.4 × 10(-3)). The median absolute difference between coalescent and nongenetic date estimates was smallest (median 29.4 days) for highly variable regions of the HIV genome such as env V3, and greater (median 114.9 days) for more conserved regions such as pol.

CONCLUSION

This application of NGS represents an important advancement, not only because accurate estimates of dates of infection can be derived retrospectively from archived specimens, but also because each analysis is patient-specific and, therefore, robust to variation in rates of HIV evolution.

A comparison of parallel pyrosequencing and sanger clone-based sequencing and its impact on the characterization of the genetic diversity of HIV-1

BACKGROUND

Pyrosequencing technology has the potential to rapidly sequence HIV-1 viral quasispecies without requiring the traditional approach of cloning. In this study, we investigated the utility of ultra-deep pyrosequencing to characterize genetic diversity of the HIV-1 gag quasispecies and assessed the possible contribution of pyrosequencing technology in studying HIV-1 biology and evolution.

METHODOLOGY/PRINCIPAL FINDINGS

HIV-1 gag gene was amplified from 96 patients using nested PCR. The PCR products were cloned and sequenced using capillary based Sanger fluorescent dideoxy termination sequencing. The same PCR products were also directly sequenced using the 454 pyrosequencing technology. The two sequencing methods were evaluated for their ability to characterize quasispecies variation, and to reveal sites under host immune pressure for their putative functional significance. A total of 14,034 variations were identified by 454 pyrosequencing versus 3,632 variations by Sanger clone-based (SCB) sequencing. 11,050 of these variations were detected only by pyrosequencing. These undetected variations were located in the HIV-1 Gag region which is known to contain putative cytotoxic T lymphocyte (CTL) and neutralizing antibody epitopes, and sites related to virus assembly and packaging. Analysis of the positively selected sites derived by the two sequencing methods identified several differences. All of them were located within the CTL epitope regions.

CONCLUSIONS/SIGNIFICANCE

Ultra-deep pyrosequencing has proven to be a powerful tool for characterization of HIV-1 genetic diversity with enhanced sensitivity, efficiency, and accuracy. It also improved reliability of downstream evolutionary and functional analysis of HIV-1 quasispecies.

Genotypic resistance testing in routine clinical care

PURPOSE OF REVIEW

Genotypic resistance testing has become part of routine clinical management of HIV-infected patients. Focussing on observational studies, this review looks at recent advances in this area.

RECENT FINDINGS

Translation of the nucleotide sequence generated by the resistance test into clinically useful information remains a major challenge. A recent key development is the availability of therapy optimization tools to predict regimens that are most likely to achieve virological suppression. Standard genotypic resistance testing only examines protease and part of reverse transcriptase; as drugs are licensed to further targets, it has become necessary to expand the repertoire for testing. Traditionally, genotypic testing has not been attempted at viral loads less than 1000 copies/ml, but recent studies indicate that major mutations are often detected at much lower levels. Similarly, various methods have been developed for the detection of minority variants including allele-specific PCR, single-genome sequencing, and ultra-deep sequencing.

SUMMARY

The technology and interpretation of genotypic resistance tests is in a phase of rapid development. It remains uncertain which of these developments will become part of routine clinical practice.

Detection of inferred CCR5- and CXCR4-using HIV-1 variants and evolutionary intermediates using ultra-deep pyrosequencing

The emergence of CXCR4-using human immunodeficiency virus type 1 (HIV-1) variants is associated with accelerated disease progression. CXCR4-using variants are believed to evolve from CCR5-using variants, but due to the extremely low frequency at which transitional intermediate variants are often present, the kinetics and mutational pathways involved in this process have been difficult to study and are therefore poorly understood. Here, we used ultra-deep sequencing of the V3 loop of the viral envelope in combination with the V3-based coreceptor prediction tools PSSM_NSI/SI and geno2pheno_[coreceptor] to detect HIV-1 variants during the transition from CCR5- to CXCR4-usage. We analyzed PBMC and serum samples obtained from eight HIV-1-infected individuals at three-month intervals up to one year prior to the first phenotypic detection of CXCR4-using variants in the MT-2 assay. Between 3,482 and 10,521 reads were generated from each sample. In all individuals, V3 sequences of predicted CXCR4-using HIV-1 were detected at least three months prior to phenotypic detection of CXCR4-using variants in the MT-2 assay. Subsequent analysis of the genetic relationships of these V3 sequences using minimum spanning trees revealed that the transition in coreceptor usage followed a stepwise mutational pathway involving sequential intermediate variants, which were generally present at relatively low frequencies compared to the major predicted CCR5- and CXCR4-using variants. In addition, we observed differences between individuals with respect to the number of predicted CXCR4-using variants, the diversity among major predicted CCR5-using variants, and the presence or absence of intermediate variants with discordant phenotype predictions. These results provide the first detailed description of the mutational pathways in V3 during the transition from CCR5- to CXCR4-usage in natural HIV-1 infection.

The first step in the infection of a target cell by human immunodeficiency virus type 1 (HIV-1) is binding of the envelope spike to its receptor CD4 and a coreceptor on the cellular surface. HIV-1 variants present early in the course of infection mainly use the coreceptor CCR5, while virus variants that use CXCR4 can appear later in infection. This change in coreceptor usage is associated with mutations in the third variable (V3) loop of the envelope spike, but has been difficult to study due to the low presence of intermediate variants. Using ultra-deep sequencing, we obtained thousands of sequences of the V3 loop from HIV-1 infected individuals in the year before CXCR4-using variants were first detected, including sequences from almost all intermediate variants. We show that mutations are introduced sequentially in the V3 loop during the evolution from CCR5- to CXCR4-usage. Furthermore, we describe differences and similarities between HIV-1-infected individuals that are related to this change in coreceptor usage, which provides the first detailed overview of this evolutionary process during natural HIV-1 infection.

Genotypic tropism testing by massively parallel sequencing: qualitative and quantitative analysis

BACKGROUND

Inferring viral tropism from genotype is a fast and inexpensive alternative to phenotypic testing. While being highly predictive when performed on clonal samples, sensitivity of predicting CXCR4-using (X4) variants drops substantially in clinical isolates. This is mainly attributed to minor variants not detected by standard bulk-sequencing. Massively parallel sequencing (MPS) detects single clones thereby being much more sensitive. Using this technology we wanted to improve genotypic prediction of coreceptor usage.

METHODS

Plasma samples from 55 antiretroviral-treated patients tested for coreceptor usage with the Monogram Trofile Assay were sequenced with standard population-based approaches. Fourteen of these samples were selected for further analysis with MPS. Tropism was predicted from each sequence with geno2pheno[coreceptor].

RESULTS

Prediction based on bulk-sequencing yielded 59.1% sensitivity and 90.9% specificity compared to the trofile assay. With MPS, 7600 reads were generated on average per isolate. Minorities of sequences with high confidence in CXCR4-usage were found in all samples, irrespective of phenotype. When using the default false-positive-rate of geno2pheno[coreceptor] (10%), and defining a minority cutoff of 5%, the results were concordant in all but one isolate.

CONCLUSIONS

The combination of MPS and coreceptor usage prediction results in a fast and accurate alternative to phenotypic assays. The detection of X4-viruses in all isolates suggests that coreceptor usage as well as fitness of minorities is important for therapy outcome. The high sensitivity of this technology in combination with a quantitative description of the viral population may allow implementing meaningful cutoffs for predicting response to CCR5-antagonists in the presence of X4-minorities.

Deep sequencing to infer HIV-1 co-receptor usage: application to three clinical trials of maraviroc in treatment-experienced patients

BACKGROUND

The Maraviroc versus Optimized Therapy in Viremic Antiretroviral Treatment-Experienced Patients (MOTIVATE) studies compared maraviroc versus placebo in treatment-experienced patients with CCR5-using (R5) human immunodeficiency virus type 1 (HIV-1), screened using the original Trofile assay. A subset with non-R5 HIV infection entered the A4001029 trial. We retrospectively examined the performance of a genotypic tropism assay based on deep sequencing of the HIV env V3 loop in predicting virologic response to maraviroc in these trials.

METHODS

V3 amplicons were prepared from 1827 screening plasma samples and sequenced on a Roche/454 GS-FLX to a depth of >3000 sequences/sample. Samples were considered non-R5 if ≥2% of their viral population scored greater than or equal to -4.75 or ≤3.5 using the PSSM(x4/R5) or geno2pheno algorithms, respectively.

RESULTS

Deep sequencing identified more than twice as many maraviroc recipients as having non-R5 HIV, compared with the original Trofile. With use of genotyping, we determined that 49% of maraviroc recipients with R5 HIV at screening had a week 48 viral load <50 copies/mL versus 26% of recipients with non-R5. Corresponding percentages were 46% and 23% with screening by Trofile. In cases in which screening assays differed, median week 8 log₁₀ copies/mL viral load decrease favored 454. Other parameters predicted by genotyping included likelihood of changing to non-R5 tropism.

CONCLUSIONS

This large study establishes deep V3 sequencing as a promising tool for identifying treatment-experienced individuals who could benefit from CCR5-antagonist-containing regimens.

Detection of quasispecies variants predicted to use CXCR4 by ultra-deep pyrosequencing during early HIV infection

OBJECTIVES

HIV-1 V3 quasispecies was analyzed by ultra-deep pyrosequencing, in early HIV-infected patients, to assess possible correlations between quasispecies diversity, frequency of variants predicted to use CXCR4 and need for early antiretroviral treatment.

METHODS

Twenty patients were retrospectively enrolled: 10 patients (group A) required HAART within 6 months from seroconversion and 10 (group B) remained free of therapy during this period. V3 quasispecies was assessed on plasma viral RNA and in peripheral blood mononuclear cell-associated proviral DNA. Prediction of coreceptor usage was performed by position-specific score matrix analysis.

RESULTS

Variants predicted to use CXCR4 were detected (frequency ≥0.3%) in the plasma of 50% of early infected patients (60% from group A and 40% from group B). Intrapatient frequency of these variants was highly variable (0.3-56.3%). A positive correlation was observed between the proportion of X4 variants and intrapatient quasispecies diversity. Quasispecies diversity and absolute numbers of X4 variants were significantly higher in patients from group A. The analysis of proviral DNA quasispecies, performed in a subgroup of five patients, showed that X4 variants were not detected in patients with RNA frequency below 0.3%, and detected at 3.6% in the patient with 56.3% of X4 plasma variants.

CONCLUSION

Our findings show that X4 variants may be frequently found, at variable intrapatient frequency, in early infected patients, and that quasispecies diversity and absolute numbers of X4 variants are significantly higher in patients undergoing early antiretroviral treatment. Further studies are mandatory to explore the clinical relevance of X4 variants present during early infection with respect to clinical progression and possible therapeutic implications.

HIV-1 dual/mixed tropic isolates show different genetic and phenotypic characteristics and response to maraviroc in vitro

Dual/mixed-tropic HIV-1 strains are predominant in a significative proportion of patients, though few information is available regarding the genetic characteristics, quasispecies composition, and susceptibility against CCR5-antagonists of the primary-isolates. For this reason, we investigated in deep details, both phenotypically and genotypically, the characteristics of 54 HIV-1 primary-isolates obtained from HIV-infected patients. Tropism was assessed by multiple-cycles phenotypic-assay on U87MG-CD4(+)-CCR5(+)-/CXCR4(+)-expressing cells. In vitro selection in PBMCs of X4-tropic viral strains following maraviroc-treatment was also performed. Phenotypic-assay reported pure R5-tropic viruses in 31 (57.4%) isolates, dual/mixed-tropic viruses in 22 (40.7%), and pure X4-tropic virus in only 1 (1.8%). Among dual/mixed-tropic isolates, 12 showed a remarkably higher replication-efficacy in CCR5-expressing cells (R5(+)/X4), and 2 in CXCR4-expressing cells (R5/X4(+)). Genotypic-tropism testing showed a correlation between PSSM-scores, geno2pheno false-positive-rate, and V3-net-charge with both CCR5-usage and syncytium-inducing ability. Moreover, specific gp120- and gp41-mutations were significantly associated with tropism and/or syncytium-inducing ability. Ultra-deep V3-pyrosequencing showed the presence of a swarm of genetically distinct species with a preference for CCR5-coreceptor not only in all pure R5-isolates, but also in 6/7 R5(+)/X4-tropic isolates. In both pure-X4 and R5/X4(+)-isolates, we observed extensive prevalence of X4-using species. In vitro selection-experiments with CCR5-inhibitor maraviroc (up to 2 months) showed no-emergence of X4-tropic variants for all R5- and R5(+)/X4-isolates tested (while X4-virus remained fully-resistant). In conclusion, our study shows that dual/mixed-tropic viruses are constituted by different species, whereby those with characteristics R5(+)/X4 are genotypically and phenotypically similar to the pure-R5 isolates; thus the use of CCR5-antagonists in patients with R5(+)/X4-tropic viruses may be a therapeutic-option that deserves further investigations.

Clinical significance of HIV-1 coreceptor usage

The identification of phenotypically distinct HIV-1 variants with different prevalence during the progression of the disease has been one of the earliest discoveries in HIV-1 biology, but its relevance to AIDS pathogenesis remains only partially understood. The physiological basis for the phenotypic variability of HIV-1 was elucidated with the discovery of distinct coreceptors employed by the virus to infect susceptible cells. The role of the viral phenotype in the variable clinical course and treatment outcome of HIV-1 infection has been extensively investigated over the past two decades. In this review, we summarize the major findings on the clinical significance of the HIV-1 coreceptor usage.

Genotypic inference of HIV-1 tropism using population-based sequencing of V3

BACKGROUND

Prior to receiving a drug from CCR5-antagonist class in HIV therapy, a patient must undergo an HIV tropism test to confirm that his or her viral population uses the CCR5 coreceptor for cellular entry, and not an alternative coreceptor. One approach to tropism testing is to examine the sequence of the V3 region of the HIV envelope, which interacts with the coreceptor.

METHODS

Viral RNA is extracted from blood plasma. The V3 region is amplified in triplicate with nested reverse transcriptase-PCR. The amplifications are then sequenced and analyzed using the software, RE_Call. Sequences are then submitted to a bioinformatic algorithm such as geno2pheno to infer viral tropism from the V3 region. Sequences are inferred to be non-R5 if their geno2pheno false positive rate falls below 5.75%. If any one of the three sequences from a sample is inferred to be non-R5, the patient is unlikely to respond to a CCR5-antagonist.

The evolutionary analysis of emerging low frequency HIV-1 CXCR4 using variants through time--an ultra-deep approach

Large-scale parallel pyrosequencing produces unprecedented quantities of sequence data. However, when generated from viral populations current mapping software is inadequate for dealing with the high levels of variation present, resulting in the potential for biased data loss. In order to apply the 454 Life Sciences' pyrosequencing system to the study of viral populations, we have developed software for the processing of highly variable sequence data. Here we demonstrate our software by analyzing two temporally sampled HIV-1 intra-patient datasets from a clinical study of maraviroc. This drug binds the CCR5 coreceptor, thus preventing HIV-1 infection of the cell. The objective is to determine viral tropism (CCR5 versus CXCR4 usage) and track the evolution of minority CXCR4-using variants that may limit the response to a maraviroc-containing treatment regimen. Five time points (two prior to treatment) were available from each patient. We first quantify the effects of divergence on initial read k-mer mapping and demonstrate the importance of utilizing population-specific template sequences in relation to the analysis of next-generation sequence data. Then, in conjunction with coreceptor prediction algorithms that infer HIV tropism, our software was used to quantify the viral population structure pre- and post-treatment. In both cases, low frequency CXCR4-using variants (2.5-15%) were detected prior to treatment. Following phylogenetic inference, these variants were observed to exist as distinct lineages that were maintained through time. Our analysis, thus confirms the role of pre-existing CXCR4-using virus in the emergence of maraviroc-insensitive HIV. The software will have utility for the study of intra-host viral diversity and evolution of other fast evolving viruses, and is available from http://www.bioinf.manchester.ac.uk/segminator/.