Risk of breast cancer with CXCR4-using HIV defined by V3 loop sequencing.

OBJECTIVE: Evaluate the risk of female breast cancer associated with HIV-CXCR4 (X4) tropism as determined by various genotypic measures. METHODS: A breast cancer case-control study, with pairwise comparisons of tropism determination methods, was conducted. From the Women's Interagency HIV Study repository, one stored plasma specimen was selected from 25 HIV-infected cases near the breast cancer diagnosis date and 75 HIV-infected control women matched for age and calendar date. HIV-gp120 V3 sequences were derived by Sanger population sequencing (PS) and 454-pyro deep sequencing (DS). Sequencing-based HIV-X4 tropism was defined using the geno2pheno algorithm, with both high-stringency DS [false-positive rate (3.5) and 2% X4 cutoff], and lower stringency DS (false-positive rate, 5.75 and 15% X4 cutoff). Concordance of tropism results by PS, DS, and previously performed phenotyping was assessed with kappa (κ) statistics. Case-control comparisons used exact P values and conditional logistic regression. RESULTS: In 74 women (19 cases, 55 controls) with complete results, prevalence of HIV-X4 by PS was 5% in cases vs 29% in controls (P = 0.06; odds ratio, 0.14; confidence interval: 0.003 to 1.03). Smaller case-control prevalence differences were found with high-stringency DS (21% vs 36%, P = 0.32), lower stringency DS (16% vs 35%, P = 0.18), and phenotyping (11% vs 31%, P = 0.10). HIV-X4 tropism concordance was best between PS and lower stringency DS (93%, κ = 0.83). Other pairwise concordances were 82%-92% (κ = 0.56-0.81). Concordance was similar among cases and controls. CONCLUSIONS: HIV-X4 defined by population sequencing (PS) had good agreement with lower stringency DS and was significantly associated with lower odds of breast cancer.

"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.

Theoretical and experimental assessment of degenerate primer tagging in ultra-deep applications of next-generation sequencing.

Primer IDs (pIDs) are random oligonucleotide tags used in next-generation sequencing to identify sequences that originate from the same template. These tags are produced by degenerate primers during the reverse transcription of RNA molecules into cDNA. The use of pIDs helps to track the number of RNA molecules carried through amplification and sequencing, and allows resolution of inconsistencies between reads sharing a pID. Three potential issues complicate the above applications. First, multiple cDNAs may share a pID by chance; we found that while preventing any cDNAs from sharing a pID may be unfeasible, it is still practical to limit the number of these collisions. Secondly, a pID must be observed in at least three sequences to allow error correction; as such, pIDs observed only one or two times must be rejected. If the sequencing product contains copies from a high number of RT templates but produces few reads, our findings indicate that rejecting such pIDs will discard a great deal of data. Thirdly, the use of pIDs could influence amplification and sequencing. We examined the effects of several intrinsic and extrinsic factors on sequencing reads at both the individual and ensemble level.

Comparative performances of HIV-1 RNA load assays at low viral load levels: results of an international collaboration.

Low-level viremia during antiretroviral therapy and its accurate measurement are increasingly relevant. Here, we present an international collaboration of 4,221 paired blood plasma viral load (pVL) results from four commercial assays, emphasizing the data with low pVL. The assays compared were the Abbott RealTime assay, the Roche Amplicor assay, and the Roche TaqMan version 1 and version 2 assays. The correlation between the assays was 0.90 to 0.97. However, at a low pVL, the correlation fell to 0.45 to 0.85. The observed interassay concordance was higher when detectability was defined as 200 copies/ml than when it was defined as 50 copies/ml. A pVL of ∼100 to 125 copies/ml by the TaqMan version 1 and version 2 assays corresponded best to a 50-copies/ml threshold with the Amplicor assay. Correlation and concordance between the viral load assays were lower at a low pVL. Clear guidelines are needed on the clinical significance of low-level viremia.

HIV drug resistance detected during low-level viraemia is associated with subsequent virologic failure.

BACKGROUND: The clinical implications of emergent HIV drug resistance on samples with low-level viraemia (LLV <1000 copies/ml) remain unclear. We undertook the present analysis to evaluate the impact of emergent HIV drug resistance at LLV on the risk of subsequent virologic failure. METHODS: One thousand, nine hundred and sixty-five patients had genotype results at LLV. Risk of virologic failure (≥1000 copies/ml) after LLV was evaluated by Kaplan-Meier analysis and Cox proportional hazards regression. Resistance was assessed using the Stanford algorithm or virtual phenotypes. Patients were grouped into four susceptibility categories ('GSS' or 'vPSS') during LLV, corresponding to the number of 'active' drugs prescribed: <1; 1-1.5; 2-2.5; and ≥3. RESULTS: A total of 1702 patients with follow-up on constant therapy were eligible for analysis. Participants excluded due to changing therapy or loss to follow-up before their next observation had mostly similar characteristics to included participants. There was a 'dose-dependent' increase in the hazard ratio for virologic failure with susceptibility categories at LLV. Compared with a GSS of at least 3, hazard ratios for virologic failure were 1.4 for GSS 2-2.5; 2.0 for GSS 1-1.5; and 3.0 for GSS less than 1 (P < 0.001). Numerous sensitivity analyses confirmed these findings. CONCLUSION: Our results demonstrate that emergent HIV drug resistance at LLV is strongly associated with subsequent virologic failure. Furthermore, we uncovered a 'dose-dependent' increase in the hazard ratio for virologic failure with decreasing GSS estimated at the time of LLV. On the basis of these findings, we propose that resistance genotyping be encouraged for HIV-infected individuals on antiretroviral therapy experiencing low-level viraemia.

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.

Results of external quality assessment for proviral DNA testing of HIV tropism in the Maraviroc Switch collaborative study.

The Maraviroc Switch collaborative study (MARCH) is a study in aviremic patients on stable antiretroviral therapy and utilizes population-based sequencing of proviral DNA to determine HIV tropism and susceptibility to maraviroc. An external quality assessment (EQA) program was implemented to ensure competency in assessing the tropism of clinical samples conducted by MARCH laboratories (n = 14). The MARCH EQA has three prestudy phases assessing V3 loop sequencing and tropism determination using the bioinformatic algorithm geno2pheno, which generates a false-positive rate (FPR). DNA sequences with low FPRs are more likely to be from CXCR4-using (X4) viruses. Phase 1 of the EQA involved chromatogram interpretation. Phases 2, 2/3, and 3 involved patient and clonal samples. Clinical samples used in these phases were from treatment-experienced HIV-infected volunteers; 18/20 had viral loads of <50 copies/ml, and 10/15 were CXCR4-tropic on prior phenotyping. All samples were tested in triplicate, and any replicate with a geno2pheno FPR of <10% was designated X4. Performance was deemed adequate if ≤2 R5 and ≤1 X4 specimens were miscalled. For several clinical samples in the EQA, triplicate testing revealed marked DNA variability (FPR range, 0 to 96.7%). Therefore, a consensus-based approach was employed for each sample, i.e., a median FPR across laboratories was used to define sample tropism. Further sequencing analysis showed mixed viral populations in the clinical samples, explaining the differences in tropism predictions. All laboratories passed the EQA after achieving predefined competence thresholds in either of the phase 2 rounds. The use of clinical samples from patients resembling those who were likely to be screened in the MARCH, coupled with triplicate testing, revealed inherent DNA variability that might have been missed if single or duplicate testing and/or clonal samples alone were used. These data highlight the importance of intensive EQA of tropism laboratories before embarking on clinical studies. (This study has been registered at ClinicalTrials.gov under registration no. NCT01384682 [http://www.clinicaltrials.gov/ct2/show/study/NCT01384682?term=NCT01384682&rank=1].).

Use of cellular HIV DNA to predict virologic response to maraviroc: performance of population-based and deep sequencing.

BACKGROUND: A tropism test is required before administration of the antiretroviral drug maraviroc. However, plasma RNA testing is not possible in patients with undetectable plasma viral loads. Here we assess genotypic testing of cellular human immunodeficiency virus (HIV) DNA from peripheral blood mononuclear cells (PBMCs) to predict virologic responses in treatment-experienced patients beginning maraviroc-containing regimens. METHODS: PBMC samples from 181 maraviroc recipients at study entry in MOTIVATE or A4001029 (51% R5 by original Trofile). The V3 loop was amplified in triplicate from cellular HIV DNA, and matching plasma RNA (n = 156). Sequencing was performed using standard population-based methods and next-generation deep sequencing, with tropism assessment as previously defined. RESULTS: Genotypic DNA-based tropism testing from the cellular compartment had 78%-81% sensitivity relative to RNA-based Trofile at the same time point. Cell-based genotypic tropism methods and plasma-based phenotypic and genotypic methods were predictive of virologic response. However, when classifications were discordant, the outcomes favored the plasma predictions over the DNA ones. CONCLUSIONS: Genotypic determination of HIV tropism can be performed using cell-derived viral DNA, and is a predictor of virologic success on maraviroc in therapy-experienced patients. However, the PBMC compartment appears to be a suboptimal predictor compared to plasma.

Development of a novel codon-specific polymerase chain reaction for the detection of CXCR4-utilizing HIV type 1 subtype B.

Insight concerning the switch in HIV-1 coreceptor use will lead to a better understanding of HIV-1 pathogenesis and host-virus dynamics. Predicting CXCR4 utilization by analyzing HIV-1 envelope consensus sequences is highly specific, but minority variants in the viral population are often missed resulting in low sensitivity. Commercial phenotypic assays are costly, and the development of sensitive in-house phenotypic assays to detect CXCR4-using HIV may not be feasible for some laboratories. A sensitive, inexpensive genotyping assay was developed to detect viral sequences associated with CXCR4-utilizing virus (X4). Codon-specific primer pairs were used to detect X4-associated codons at five positions in the HIV-1 envelope V3 loop (11, 13, 24, 25, and 32). Sixty plasma samples from HIV-1-infected individuals were analyzed by consensus sequencing and codon-specific PCR (CS-PCR). Forty-six of these were also phenotyped by Trofile or Enhanced Sensitivity Trofile (ESTA). CS-PCR detected X4 variants 17% more often than 11/24/25 consensus sequencing alone (n=60), 30% more often than Trofile (n=27), and in a limited data set, 16% more often than ESTA (n=19). CS-PCR combined with consensus sequencing had approximately 80% concordance with ESTA.

Factors influencing the sensitivity and specificity of conventional sequencing in human immunodeficiency virus type 1 tropism testing.

Human immunodeficiency virus type 1 (HIV-1) V3 loop sequence can be used to infer viral coreceptor use. The effect of input copy number on population-based sequencing of the V3 loop of HIV-1 was examined through replicate deep and population-based sequencing of samples with known tropism, a heterogeneous clinical sample (624 population-based sequences and 47 deep-sequencing replicates), and a large cohort of clinical samples from phase III clinical trials of maraviroc including the MOTIVATE/A4001029 studies (n = 1,521). Proviral DNA from two independent samples from each of 101 patients from the MOTIVATE/A4001029 studies was also analyzed. Cumulative technical error occurred at a rate of 3 × 10(-4) mismatches/bp, without observed effect on inferred tropism. Increasing PCR replication increased minority species detection with an ~10% minority population detected in 18% of cases using a single replicate at a viral load of 1,072 copies/ml and in 44% of cases using three replicates. The nucleotide prevalence detected by population-based and deep sequencing were highly correlated (Spearman's ρ, 0.73), and the accuracy increased with increasing input copy number (P < 0.001). Triplicate sequencing was able to predict tropism changes in the MOTIVATE/A4001029 studies for both low (P = 0.05) and high (P = 0.02) viral loads. Sequences derived from independently extracted and processed samples of proviral DNA for the same patient were equivalent to replicates from the same extraction (P = 0.45) and had correlated position-specific scoring matrix scores (Spearman's ρ, 0.75; P << 0.001); however, concordance in tropism inference was only 83%. Input copy number and PCR replication are important factors in minority species detection in samples with significant heterogeneity.

Reconstructing the dynamics of HIV evolution within hosts from serial deep sequence data.

At the early stage of infection, human immunodeficiency virus (HIV)-1 predominantly uses the CCR5 coreceptor for host cell entry. The subsequent emergence of HIV variants that use the CXCR4 coreceptor in roughly half of all infections is associated with an accelerated decline of CD4+ T-cells and rate of progression to AIDS. The presence of a 'fitness valley' separating CCR5- and CXCR4-using genotypes is postulated to be a biological determinant of whether the HIV coreceptor switch occurs. Using phylogenetic methods to reconstruct the evolutionary dynamics of HIV within hosts enables us to discriminate between competing models of this process. We have developed a phylogenetic pipeline for the molecular clock analysis, ancestral reconstruction, and visualization of deep sequence data. These data were generated by next-generation sequencing of HIV RNA extracted from longitudinal serum samples (median 7 time points) from 8 untreated subjects with chronic HIV infections (Amsterdam Cohort Studies on HIV-1 infection and AIDS). We used the known dates of sampling to directly estimate rates of evolution and to map ancestral mutations to a reconstructed timeline in units of days. HIV coreceptor usage was predicted from reconstructed ancestral sequences using the geno2pheno algorithm. We determined that the first mutations contributing to CXCR4 use emerged about 16 (per subject range 4 to 30) months before the earliest predicted CXCR4-using ancestor, which preceded the first positive cell-based assay of CXCR4 usage by 10 (range 5 to 25) months. CXCR4 usage arose in multiple lineages within 5 of 8 subjects, and ancestral lineages following alternate mutational pathways before going extinct were common. We observed highly patient-specific distributions and time-scales of mutation accumulation, implying that the role of a fitness valley is contingent on the genotype of the transmitted variant.

Prolonged and substantial discordance in prevalence of raltegravir-resistant HIV-1 in plasma versus PBMC samples revealed by 454 "deep" sequencing.

The evolution of drug resistance mutations in plasma samples is relatively well-characterized. However, the viral population and diversity in other body compartments such as peripheral blood mononuclear cells (PBMC) remains poorly understood. Previous studies have mostly focused on protease and reverse transcriptase drug resistance mutations (DRMs). In this study, we used 454 "deep" sequencing technology to observe and quantify longitudinally the prevalence of resistance mutations associated with the integrase inhibitor, raltegravir, in plasma versus PBMC samples from a San Francisco-based cohort. Four heavily treatment-experienced subjects were monitored in this study over a median of 1.2 years since the initiation of raltegravir-containing regimens. We observed a consistent discordance in the prevalence of DRMs, but not resistance pathway(s), in the plasma versus PBMC viral populations. In the final paired samples that were tested while the subjects were on a raltegravir-containing regimen, DRM prevalence reached 100% in plasma but remained 1% in PBMC on day 177 post-therapy in Subject 3180 (Q148H/G140S), 100% in plasma and 36% in PBMC on day 224 in Subject 3242 (N155H), 78% in plasma and 11-12% in PBMC on day 338 in Subject 3501 (Q148H/G140S), and 100% in plasma and 0% in PBMC on day 197 in Subject 3508 (Y143R). Furthermore, absolute sequence homology comparison between the two compartments revealed that 21% - 99% of PBMC sequences had no match in plasma, whereas 14% - 100% of plasma sequences had no match in PBMC. Overall, our observations suggested that plasma and PBMC hosted drastically different HIV-1 populations even after a prolonged exposure to raltegravir selection pressure.

Technical and regulatory shortcomings of the TaqMan version 1 HIV viral load assay.

BACKGROUND: The lower limit of detection of the original Roche Amplicor HIV plasma viral load (pVL) assay (50 copies/mL) has defined HIV treatment success. The Amplicor assay, however, has been replaced by the Roche TaqMan assay(s). Changes to the limits of detection and calibration have not been validated for clinical utility. Sudden increases in the number of patients with detectable pVL have been reported following the introduction of the TaqMan version 1 assay. METHODS: Between October 2009 and April 2010 all routine pVL samples from British Columbia, Canada, with 40-250 copies/mL by TaqMan were re-tested by Amplicor (N = 1198). Subsequent short-term virological and resistance outcomes were followed in patients with unchanged therapy (N = 279; median 3.2 months follow-up). RESULTS: TaqMan and Amplicor values correlated poorly at low pVL values. Low-level pVL by TaqMan was not associated with impending short-term virological failure; only 17% of patients with 40-250 copies/mL by TaqMan had detectable pVL by Amplicor at follow-up. During the follow-up period only 20% of patients had an increase in pVL by TaqMan (median [IQR]: 80 [36-283] copies/mL). In addition, in ~2.4% of samples pVL was dramatically underestimated by TaqMan due to poor binding of the proprietary TaqMan primers. CONCLUSIONS: The replacement of Amplicor with the TaqMan assay has altered the previously accepted definition of HIV treatment failure without any evidence to support the clinical relevance of the new definition. Given the systematic differences in measurement in the low pVL range the British Columbia HIV treatment guidelines now use a threshold of >250 copies/mL by TaqMan to define treatment failure.

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.

Deep V3 sequencing for HIV type 1 tropism in treatment-naive patients: a reanalysis of the MERIT trial of maraviroc.

BACKGROUND: Deep sequencing is a highly sensitive technique that can detect and quantify the proportion of non-R5 human immunodeficiency virus (HIV) variants, including small minorities, that may emerge and cause virologic failure in patients who receive maraviroc-containing regimens. We retrospectively tested the ability of deep sequencing to predict response to a maraviroc-containing regimen in the Maraviroc versus Efavirenz in Treatment-Naive Patients (MERIT) trial. Results were compared with those obtained using the Enhanced Sensitivity Trofile Assay (ESTA), which is widely used in clinical practice. METHODS: Screening plasma samples from treatment-naive patients who received maraviroc and efavirenz in the MERIT trial were assessed. Samples were extracted, and the V3 region of HIV type 1 glycoprotein 120 was amplified in triplicate and combined in equal quantities before sequencing on a Roche/454 Genome Sequencer-FLX (n = 859). Tropism was inferred from third variable (V3) sequences, with samples classified as non-R5 if ≥2% of the viral population scored ≤3.5 using geno2pheno. RESULTS: Deep sequencing distinguished between responders and nonresponders to maraviroc. Among patients identified as having R5-HIV by deep sequencing, 67% of maraviroc recipients and 69% of efavirenz recipients had a plasma viral load <50 copies/mL at week 48, similar to the ESTA results: 68% and 68%, respectively. CONCLUSIONS: Reanalysis of the MERIT trial using deep V3 loop sequencing indicates that, had patients originally been screened using this method, the maraviroc arm would have likely been found to be noninferior to the efavirenz arm.

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.

Mutations in gp41 are correlated with coreceptor tropism but do not improve prediction methods substantially.

BACKGROUND: The main determinants of HIV-1 coreceptor usage are located in the V3-loop of gp120, although mutations in V2 and gp41 are also known. Incorporation of V2 is known to improve prediction algorithms; however, this has not been confirmed for gp41 mutations. METHODS: Samples with V3 and gp41 genotypes and Trofile assay (Monogram Biosciences, South San Francisco, CA, USA) results were taken from the HOMER cohort (n=444) and from patients screened for the MOTIVATE studies (n=1,916; 859 with maraviroc outcome data). Correlations of mutations with tropism were assessed using Fisher's exact test and prediction models trained using support vector machines. Models were validated by cross-validation, by testing models from one dataset on the other, and by analysing virological outcome. RESULTS: Several mutations within gp41 were highly significant for CXCR4 usage; most strikingly an insertion occurring in 7.7% of HOMER-R5 and 46.3% of HOMER-X4 samples (MOTIVATE 5.7% and 25.2%, respectively). Models trained on gp41 sequence alone achieved relatively high areas under the receiver-operating characteristic curve (AUCs; HOMER 0.713 and MOTIVATE 0.736) that were almost as good as V3 models (0.773 and 0.884, respectively). However, combining the two regions improved predictions only marginally (0.813 and 0.902, respectively). Similar results were found when models were trained on HOMER and validated on MOTIVATE or vice versa. The difference in median log viral load decrease at week 24 between patients with R5 and X4 virus was 1.65 (HOMER 2.45 and MOTIVATE 0.79) for V3 models, 1.59 for gp41-models (2.42 and 0.83, respectively) and 1.58 for the combined predictor (2.44 and 0.86, respectively). CONCLUSIONS: Several mutations within gp41 showed strong correlation with tropism in two independent datasets. However, incorporating gp41 mutations into prediction models is not mandatory because they do not improve substantially on models trained on V3 sequences alone.

"Dynamic range" of inferred phenotypic HIV drug resistance values in clinical practice.

BACKGROUND: 'Virtual' or inferred phenotypes (vPhenotypes) are commonly used to assess resistance to antiretroviral agents in patients failing therapy. In this study, we provide a clinical context for understanding vPhenotype values. METHODS: All HIV-infected persons enrolled in the British Columbia Drug Treatment Program with a baseline plasma viral load (pVL) and follow-up genotypic resistance and pVL results were included up to October 29, 2008 (N = 5,277). Change from baseline pVL was determined as a function of Virco vPhenotype, and the "dynamic range" (defined here by the 10th and 90th percentiles for fold-change in IC50 amongst all patients) was estimated from the distribution of vPhenotye fold-changes across the cohort. RESULTS: The distribution of vPhenotypes from a large cohort of HIV patients who have failed therapy are presented for all available antiretroviral agents. A maximum change in IC50 of at least 13-fold was observed for all drugs. The dideoxy drugs, tenofovir and most PIs exhibited small "dynamic ranges" with values of <4-fold change observed in > 99% of samples. In contrast, zidovudine, lamivudine, emtricitabine and the non-nucleoside reverse transcriptase inihibitors (excluding etravirine) had large dynamic ranges. CONCLUSION: We describe the populational distribution of vPhenotypes such that vPhenotype results can be interpreted relative to other patients in a drug-specific manner.

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 log10 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.