"Deep" sequencing accuracy and reproducibility using Roche/454 technology for inferring co-receptor usage in HIV-1.

Next generation, "deep", sequencing has increasing applications both clinically and in disparate fields of research. This study investigates the accuracy and reproducibility of "deep" sequencing as applied to co-receptor prediction using the V3 loop of Human Immunodeficiency Virus-1. Despite increasing use in HIV co-receptor prediction, the accuracy and reproducibility of deep sequencing technology, and the factors which can affect it, have received only a limited level of investigation. To accomplish this, repeated deep sequencing results were generated using the Roche GS-FLX (454) from a number of sources including a non-homogeneous clinical sample (N = 47 replicates over 18 deep sequencing runs), and a large clinical cohort from the MOTIVATE and A400129 studies (N = 1521). For repeated measurements of a non-homogeneous clinical sample, increasing input copy number both decreased variance in the measured proportion of non-R5 using virus (p<<0.001 and 0.02 for single replicates and triplicates respectively) and increased measured viral diversity (p<0.001; multiple measures). Detection of sequences with a mean abundance less than 1% abundance showed a 2 fold increase in median coefficient of variation (CV) in repeated measurements of a non-homogeneous clinical sample, and a 2.7 fold increase in CV in the MOTIVATE/A400129 dataset compared to sequences with ≥1% abundance. An unexpected source of error included read position, with low accuracy reads occurring more frequently towards the edge of sequencing regions (p<<0.001). Overall, the primary source of variability was sampling error caused by low input copy number/minority species prevalence, though other sources of error including sequence intrinsic, temporal, and read-position related errors were detected.

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.

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

BACKGROUND: Tropism testing should rule out CXCR4-using HIV before treatment with CCR5 antagonists. Currently, the recombinant phenotypic Trofile assay (Monogram) is most widely utilized; however, genotypic tests may represent alternative methods. METHODS: Independent triplicate amplifications of the HIV gp120 V3 region were made from either plasma HIV RNA or proviral DNA. These underwent standard, population-based sequencing with an ABI3730 (RNA n = 63; DNA n = 40), or "deep" sequencing with a Roche/454 Genome Sequencer-FLX (RNA n = 12; DNA n = 12). Position-specific scoring matrices (PSSMX4/R5) (-6.96 cutoff) and geno2pheno[coreceptor] (5% false-positive rate) inferred tropism from V3 sequence. These methods were then independently validated with a separate, blinded dataset (n = 278) of screening samples from the maraviroc MOTIVATE trials. RESULTS: Standard sequencing of HIV RNA with PSSM yielded 69% sensitivity and 91% specificity, relative to Trofile. The validation dataset gave 75% sensitivity and 83% specificity. Proviral DNA plus PSSM gave 77% sensitivity and 71% specificity. "Deep" sequencing of HIV RNA detected >2% inferred-CXCR4-using virus in 8/8 samples called non-R5 by Trofile, and <2% in 4/4 samples called R5. CONCLUSIONS: Triplicate analyses of V3 standard sequence data detect greater proportions of CXCR4-using samples than previously achieved. Sequencing proviral DNA and "deep" V3 sequencing may also be useful tools for assessing tropism.

HLA-associated immune escape pathways in HIV-1 subtype B Gag, Pol and Nef proteins.

BACKGROUND: Despite the extensive genetic diversity of HIV-1, viral evolution in response to immune selective pressures follows broadly predictable mutational patterns. Sites and pathways of Human Leukocyte-Antigen (HLA)-associated polymorphisms in HIV-1 have been identified through the analysis of population-level data, but the full extent of immune escape pathways remains incompletely characterized. Here, in the largest analysis of HIV-1 subtype B sequences undertaken to date, we identify HLA-associated polymorphisms in the three HIV-1 proteins most commonly considered in cellular-based vaccine strategies. Results are organized into protein-wide escape maps illustrating the sites and pathways of HLA-driven viral evolution. METHODOLOGY/PRINCIPAL FINDINGS: HLA-associated polymorphisms were identified in HIV-1 Gag, Pol and Nef in a multicenter cohort of >1500 chronically subtype-B infected, treatment-naïve individuals from established cohorts in Canada, the USA and Western Australia. At q< or =0.05, 282 codons commonly mutating under HLA-associated immune pressures were identified in these three proteins. The greatest density of associations was observed in Nef (where close to 40% of codons exhibited a significant HLA association), followed by Gag then Pol (where approximately 15-20% of codons exhibited HLA associations), confirming the extensive impact of immune selection on HIV evolution and diversity. Analysis of HIV codon covariation patterns identified over 2000 codon-codon interactions at q< or =0.05, illustrating the dense and complex networks of linked escape and secondary/compensatory mutations. CONCLUSIONS/SIGNIFICANCE: The immune escape maps and associated data are intended to serve as a user-friendly guide to the locations of common escape mutations and covarying codons in HIV-1 subtype B, and as a resource facilitating the systematic identification and classification of immune escape mutations. These resources should facilitate research in HIV epitope discovery and host-pathogen co-evolution, and are relevant to the continued search for an effective CTL-based AIDS vaccine.

Human leukocyte antigen-specific polymorphisms in HIV-1 Gag and their association with viral load in chronic untreated infection.

OBJECTIVE: Selection of specific human leukocyte antigen (HLA)-restricted cytotoxic T-lymphocyte (CTL) escape mutations in key Gag epitopes has been associated with loss of HIV immune control on an individual basis. Here we undertake a population-based identification of HLA-associated polymorphisms in Gag and investigate their relationship with plasma viral load. DESIGN: Cross-sectional analysis of 567 chronically HIV subtype B-infected, treatment-naive individuals. METHODS: HLA class I-associated Gag substitutions were identified using phylogenetically corrected analysis methods featuring a multivariate adjustment for HLA linkage disequilibrium and a q-value correction for multiple tests. Presence of HLA-associated substitutions and markers of HIV disease status were correlated using Spearman's rank test. RESULTS: We have created a gene-wide map of HLA class I-associated substitutions in HIV-1 subtype B Gag. This features 111 HLA-associated substitutions occurring at 51 of 500 Gag codons, more than 50% of which occur within published and/or putative HLA-restricted CTL epitopes. A modest inverse correlation was observed between the total number of HLA-associated Gag polymorphic sites within each individual and plasma viral load in chronic untreated infection (R = -0.17, P < 0.0001), supporting the hypothesis that a broad ability to target Gag in vivo contributes to viral control. A modest positive correlation was observed between the proportion of these sites exhibiting HLA-associated substitutions and plasma viral load (R = 0.09, P = 0.03), consistent with a loss of viremia control with the accumulation of CTL escape mutations. CONCLUSION: Results contribute to our understanding of immune-driven viral adaptation and suggest that the accumulation of CTL escape mutations in Gag results in clinically detectable consequences at the population level. These data have implications for HIV vaccines.

Current V3 genotyping algorithms are inadequate for predicting X4 co-receptor usage in clinical isolates.

OBJECTIVE: Integrating CCR5 antagonists into clinical practice would benefit from accurate assays of co-receptor usage (CCR5 versus CXCR4) with fast turnaround and low cost. DESIGN: Published HIV V3-loop based predictors of co-receptor usage were compared with actual phenotypic tropism results in a large cohort of antiretroviral naive individuals to determine accuracy on clinical samples and identify areas for improvement. METHODS: Aligned HIV envelope V3 loop sequences (n = 977), derived by bulk sequencing were analyzed by six methods: the 11/25 rule; a neural network (NN), two support vector machines, and two subtype-B position specific scoring matrices (PSSM). Co-receptor phenotype results (Trofile Co-receptor Phenotype Assay; Monogram Biosciences) were stratified by CXCR4 relative light unit (RLU) readout and CD4 cell count. RESULTS: Co-receptor phenotype was available for 920 clinical samples with V3 genotypes having fewer than seven amino acid mixtures (n = 769 R5; n = 151 X4-capable). Sensitivity and specificity for predicting X4 capacity were evaluated for the 11/25 rule (30% sensitivity/93% specificity), NN (44%/88%), PSSM(sinsi) (34%/96%), PSSM(x4r5) (24%/97%), SVMgenomiac (22%/90%) and SVMgeno2pheno (50%/89%). Quantitative increases in sensitivity could be obtained by optimizing the cut-off for methods with continuous output (PSSM methods), and/or integrating clinical data (CD4%). Sensitivity was directly proportional to strength of X4 signal in the phenotype assay (P < 0.05). CONCLUSIONS: Current default implementations of co-receptor prediction algorithms are inadequate for predicting HIV X4 co-receptor usage in clinical samples, particularly those X4 phenotypes with low CXCR4 RLU signals. Significant improvements can be made to genotypic predictors, including training on clinical samples, using additional data to improve predictions and optimizing cutoffs and increasing genotype sensitivity.