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.

Increasingly successful highly active antiretroviral therapy delays the emergence of new HLA class I-associated escape mutations in HIV-1.

BACKGROUND: HLA class I-restricted cytotoxic T lymphocytes and highly active antiretroviral therapy (HAART) exert strong selective pressures on human immunodeficiency virus type 1 (HIV-1), leading to escape mutations compromising virologic control. Immune responses continue to shape HIV-1 evolution after HAART initiation, but the extent and rate at which this occurs remain incompletely quantified. Here, we characterize the incidence and clinical correlates of HLA-associated evolution in HIV-1 Pol after HAART initiation in a large, population-based observational cohort. METHODS: British Columbia HAART Observational, Medical Evaluation and Research cohort participants with available HLA class I types and longitudinal posttherapy protease/reverse transcriptase sequences were studied (n = 619; median, 5 samples per patient and 5.2 years of follow-up). HLA-associated polymorphisms were defined according to published reference lists. Rates and correlates of immune-mediated HIV-1 evolution were investigated using multivariate Cox proportional hazard models incorporating baseline and time-dependent plasma viral load and CD4 response data. RESULTS: New HLA-associated escape events were observed in 269 (43%) patients during HAART and occurred at 49 of 63 (78%) investigated immune-associated sites in Pol. In time-dependent analyses adjusting for baseline factors, poorer virologic, but not immunologic, response to HAART was associated with increased risk of immune escape of 1.9-fold per log(10) viral load increment (P < .0001). Reversion of escape mutations following HAART initiation was extremely rare. CONCLUSIONS: HLA-associated HIV-1 evolution continues during HAART to an extent that is inversely related to the virologic success of therapy. Minimizing the degree of immune escape could represent a secondary benefit of effective HAART.

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.

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.

Phylogenetic analysis of population-based and deep sequencing data to identify coevolving sites in the nef gene of HIV-1.

Rapidly evolving viruses such as HIV-1 display extensive sequence variation in response to host-specific selection, while simultaneously maintaining functions that are critical to replication and infectivity. This apparent conflict between diversifying and purifying selection may be resolved by an abundance of epistatic interactions such that the same functional requirements can be met by highly divergent sequences. We investigate this hypothesis by conducting an extensive characterization of sequence variation in the HIV-1 nef gene that encodes a highly variable multifunctional protein. Population-based sequences were obtained from 686 patients enrolled in the HOMER cohort in British Columbia, Canada, from which the distribution of nonsynonymous substitutions in the phylogeny was reconstructed by maximum likelihood. We used a phylogenetic comparative method on these data to identify putative epistatic interactions between residues. Two interactions (Y120/Q125 and N157/S169) were chosen to further investigate within-host evolution using HIV-1 RNA extractions from plasma samples from eight patients. Clonal sequencing confirmed strong linkage between polymorphisms at these sites in every case. We used massively parallel pyrosequencing (MPP) to reconstruct within-host evolution in these patients. Experimental error associated with MPP was quantified by performing replicates at two different stages of the protocol, which were pooled prior to analysis to reduce this source of variation. Phylogenetic reconstruction from these data revealed correlated substitutions at Y120/Q125 or N157/S169 repeated across multiple lineages in every host, indicating convergent within-host evolution shaped by epistatic interactions.

Validation of a Genotype-Independent Hepatitis C Virus Near-Whole Genome Sequencing Assay.

Despite the effectiveness of direct-acting antiviral agents in treating hepatitis C virus (HCV), cases of treatment failure have been associated with the emergence of resistance-associated substitutions. To better guide clinical decision-making, we developed and validated a near-whole-genome HCV genotype-independent next-generation sequencing strategy. HCV genotype 1-6 samples from direct-acting antiviral agent treatment-naïve and -treated HCV-infected individuals were included. Viral RNA was extracted using a NucliSens easyMAG and amplified using nested reverse transcription-polymerase chain reaction. Libraries were prepared using Nextera XT and sequenced on the Illumina MiSeq sequencing platform. Data were processed by an in-house pipeline (MiCall). Nucleotide consensus sequences were aligned to reference strain sequences for resistance-associated substitution identification and compared to NS3, NS5a, and NS5b sequence data obtained from a validated in-house assay optimized for HCV genotype 1. Sequencing success rates (defined as achieving >100-fold read coverage) approaching 90% were observed for most genotypes in samples with a viral load >5 log10 IU/mL. This genotype-independent sequencing method resulted in >99.8% nucleotide concordance with the genotype 1-optimized method, and 100% agreement in genotype assignment with paired line probe assay-based genotypes. The assay demonstrated high intra-run repeatability and inter-run reproducibility at detecting substitutions above 2% prevalence. This study highlights the performance of a freely available laboratory and bioinformatic approach for reliable HCV genotyping and resistance-associated substitution detection regardless of genotype.

Transmission of drug-resistant HIV-1 from an infected individual to a caregiver.

OBJECTIVE: To describe a suspected case of intrafamilial transmission of drug-resistant HIV-1 from an infected individual to his caregiver. METHODS: HIV-1 protease, reverse transcriptase, nef and gp41 DNA sequences were determined from plasma samples after RNA extraction and nested RT-PCR. Phylogenetic trees were generated using neighbour-joining methods. RESULTS: HIV-1 infection was diagnosed in an individual (index) following the death of her HIV-infected adopted son (source). The index case cared for her son in the 9 months prior to his death. No obvious instances of contact with bodily secretions or other traditional risk factors for HIV transmission were identified, although the index case had mild eczema. At time of death, the source case's HIV plasma viral load was > 1,000,000 HIV copies/ml. Genotyping of the protease and reverse transcriptase of both the source and index samples revealed similar drug-resistance mutations despite the index case being antiretroviral-naive. Both source and index V3 genotypes were consistent with syncytium-inducing virus. The mean pairwise amino acid identity between the source and index samples was 97.0%, 99.2%, 92.6% and 94.6% in protease, reverse transcriptase, nef and gp41, respectively. The source and index sequences clustered together on phylogenetic trees when compared with 50 randomly selected sequences from British Columbia. CONCLUSION: Although exceptionally rare, our case demonstrates transmission of drug-resistant HIV-1 from an infected individual to his caregiver in the absence of traditional risk factors. It should be noted that the source had a very high viral load during the terminal phase of his illness.

Emergence of drug resistance is associated with an increased risk of death among patients first starting HAART.

BACKGROUND: The impact of the emergence of drug-resistance mutations on mortality is not well characterized in antiretroviral-naïve patients first starting highly active antiretroviral therapy (HAART). Patients may be able to sustain immunologic function with resistant virus, and there is limited evidence that reduced sensitivity to antiretrovirals leads to rapid disease progression or death. We undertook the present analysis to characterize the determinants of mortality in a prospective cohort study with a median of nearly 5 y of follow-up. The objective of this study was to determine the impact of the emergence of drug-resistance mutations on survival among persons initiating HAART. METHODS AND FINDINGS: Participants were antiretroviral therapy naïve at entry and initiated triple combination antiretroviral therapy between August 1, 1996, and September 30, 1999. Marginal structural modeling was used to address potential confounding between time-dependent variables in the Cox proportional hazard regression models. In this analysis resistance to any class of drug was considered as a binary time-dependent exposure to the risk of death, controlling for the effect of other time-dependent confounders. We also considered each separate class of mutation as a binary time-dependent exposure, while controlling for the presence/absence of other mutations. A total of 207 deaths were identified among 1,138 participants over the follow-up period, with an all cause mortality rate of 18.2%. Among the 679 patients with HIV-drug-resistance genotyping done before initiating HAART, HIV-drug resistance to any class was observed in 53 (7.8%) of the patients. During follow-up, HIV-drug resistance to any class was observed in 302 (26.5%) participants. Emergence of any resistance was associated with mortality (hazard ratio: 1.75 [95% confidence interval: 1.27, 2.43]). When we considered each class of resistance separately, persons who exhibited resistance to non-nucleoside reverse transcriptase inhibitors had the highest risk: mortality rates were 3.02 times higher (95% confidence interval: 1.99, 4.57) for these patients than for those who did not exhibit this type of resistance. CONCLUSIONS: We demonstrated that emergence of resistance to non-nucleoside reverse transcriptase inhibitors was associated with a greater risk of subsequent death than was emergence of protease inhibitor resistance. Future research is needed to identify the particular subpopulations of men and women at greatest risk and to elucidate the impact of resistance over a longer follow-up period.

Rates of antiretroviral resistance among HIV-infected patients with and without a history of injection drug use.

BACKGROUND: There exist concerns regarding the potential for elevated rates of antiretroviral resistance among HIV-infected injection drug users (IDUs) prescribed highly active antiretroviral therapy (HAART), however, no population-based study has examined if IDUs have elevated rates of antiretroviral resistance in comparison to non-IDUs. OBJECTIVE: To evaluate the time to the development of antiretroviral resistance among antiretroviral-naive patients with and without a history of injection drug use. METHODS: In British Columbia there is a province-wide HIV/AIDS treatment program that provides antiretrovirals free of charge. We examined all antiretroviral-naive patients initiating HAART between 1 August 1996 and 30 September 2000 and who were followed to 31 March 2002. The main outcome measure was the time to class-specific antiretroviral resistance. Cumulative antiretroviral resistance rates among IDUs and non-IDUs were evaluated using Kaplan-Meier methods and relative hazards were estimated using Cox regression. RESULTS: Overall, 1191 antiretroviral-naive patients initiated HAART during the study period. Resistance mutations were observed in 298 (25%) subjects during the first 30 months of HAART. In comparison with non-IDUs, the risk of protease inhibitor resistance [relative hazard (RH), 0.9; 95% confidence interval (CI), 0.5-1.6] and non-nucleoside reverse transcriptase inhibitor resistance (RH, 1.5; 95% CI, 1.0-2.2) were similar among IDUs, and there were no differences in the rates of resistance to the sub-classes of nucleoside reverse transcriptase inhibitors. CONCLUSIONS: Resistance to all major classes of antiretrovirals were similar among IDUs and non-IDUs after 30 months of follow-up. These findings should help to allay fears that prescribing HAART to IDUs may result in elevated rates of resistance.

Influence of polymorphisms within the CX3CR1 and MDR-1 genes on initial antiretroviral therapy response.

OBJECTIVE: Single nucleotide polymorphisms (SNP) in the genes encoding the human CX3CR1 chemokine receptor and the P-glycoprotein multidrug transporter have been associated with accelerated disease progression in untreated individuals and implicated in therapeutic response, respectively. This retrospective study assessed the influence of SNP in the CX3CR1 and MDR-1 genes on initial virological and immunological response in 461 HIV-infected, antiretroviral-naive individuals initiating antiretroviral therapy in British Columbia, Canada. METHODS: CX3CR1 and MDR-1 SNP were determined by PCR amplification of human DNA from plasma, followed by DNA sequencing. Time to virological success [time to HIV plasma viral load (pVL) < or = 500 copies/ml], virological failure (subsequent time to the second of two consecutive pVL > or = 500) and immunological failure (time to the second consecutive CD4 cell count below baseline) were analyzed by Kaplan-Meier methods. RESULTS: Frequencies of CX3CR1 amino acid haplotypes were 249V 280T (0.75), 249I 280M (0.15), and 249I 280T (0.1). Frequencies of MDR-1 nucleotide polymorphisms were 3435C (0.47) and 3435T (0.53). There was no effect detected for SNP in CX3CR1 or MDR-1 on time to virological success, nor of CX3CR1 and MDR-1 SNP on time to virological and immunological failure, respectively ( P > 0.1). There was a trend to earlier virological failure in the MDR-1 3435C/C genotype group ( P = 0.07), and a statistically significant trend to earlier immunological failure in individuals with the CX3CR1 249I polymorphism ( P = 0.02). These remained significant after correcting for baseline age, sex, pVL, CD4 cell count, type of therapy, and adherence ( P < or = 0.05). CONCLUSION: Polymorphisms in MDR-1 and CX3CR1 may be associated with accelerated virological and immunological therapy failure, respectively.

Clinical and immunological impact of HIV envelope V3 sequence variation after starting initial triple antiretroviral therapy.

BACKGROUND: The HIV-1 envelope third variable loop (V3 loop) is an important determinant of viral phenotype and co-receptor usage. We wished to determine the impact of specific V3 genotypes associated with viral phenotype and co-receptor usage on response to initial triple antiretroviral therapy. METHODS: Pre-therapy plasma samples from the HOMER cohort of 1191 antiretroviral-naive, HIV-infected adults who initiated triple therapy in British Columbia, Canada between August 1996 and September 1999 were genotyped for V3 loop sequence. V3 sequences were dichotomized by the presence or absence of positively charged residues at codons 11 and/or 25 (an '11/25' genotype). Neural network (NN) and Position Specific Scoring Matrix (PSSM) approaches were used as alternative V3 sequence interpretation methods. The association of V3 genotypes with clinical endpoints was assessed over a median of 43 months of follow up. RESULTS: One-hundred and eighteen (10.9%) of the 1085 isolates successfully genotyped for V3 displayed the 11/25 genotype. In multivariate analyses, this genotype was associated with a more rapid CD4 decline [risk ratio, (RR), 1.38; P = 0.012] and earlier mortality (RR, 1.70; P = 0.027), despite comparable viral load suppression below 500 HIV RNA copies/ml. We observed no influence of the 11/25 genotype on time to viral rebound or the development of drug resistance. PSSM-based sequence categories were similarly predictive of outcomes. NN sequence categories were not associated with any endpoints. CONCLUSION: The 11/25 genotype of the HIV V3 loop is an independent predictor of poor immunological response and more rapid mortality even after starting triple antiretroviral therapy. These results may prove to be useful for the clinical management of HIV-infected individuals.

Prevalence and clinical implications of insertions in the HIV-1 p6Gag N-terminal region in drug-naive individuals initiating antiretroviral therapy.

We assessed the prevalence and clinical impact of insertions within the HIV-1 p6Gag proline-rich (PTAP) region on initial antiretroviral therapy response in 461 HIV-infected, drug-naive individuals initiating therapy in British Columbia, Canada between June 1996 and August 1998. HIV p6Gag insertions were detected by nested RT-PCR of extracted patient plasma followed by direct DNA sequencing. Insertions were observed in 70 of 423 successfully genotyped samples (16.5%). HIV p6Gag insertions were significantly associated with a lower baseline CD4 cell count (P<<0.05) and the presence of basic amino acids at key positions in the HIV envelope V3 loop linked to a syncytium-inducing phenotype (P<<0.05). After adjusting for baseline factors, no effect of HIV p6Gag insertions was observed on time to virological success (pVL < or = 500 copies/ml), virological failure (subsequent confirmed pVL > or = 500 copies/ml) or immunological failure (confirmed CD4 count below baseline), as evaluated by Kaplan-Meier methods and Cox proportional hazard regression (P>0.1). The data suggest that HIV p6Gag insertions are not exclusively related to drug resistance and may not influence response to antiretroviral therapy, but may be linked to sequence variations in the HIV envelope.

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.

Predictors of HIV drug-resistance mutations in a large antiretroviral-naive cohort initiating triple antiretroviral therapy.

OBJECTIVE: The objective of this study was to systematically characterize the incidence and determinants of antiretroviral resistance in the HOMER (Highly Active Antiretroviral Therapy [HAART] Observational Medical Evaluation and Research) cohort of 1191 human immunodeficiency virus-infected, antiretroviral-naive adults initiating HAART in British Columbia, Canada. METHODS: All plasma samples with plasma virus loads (pVLs) >1000 copies/mL collected during the first 30 months of follow-up were genotyped for drug resistance. The primary outcome measure was time to the first detection of major drug-resistance mutation(s). Cox proportional hazard regression was used to identify factors significantly associated with the detection of drug-resistance mutations. RESULTS: Drug-resistance mutations were detected in 298 subjects (25%). Factors significantly associated with detection of drug-resistance mutations included high baseline pVL (multivariate hazard ratio [HR], 1.59; P<.001) and adherence (estimated using prescription-refill data and/or untimed plasma drug-concentration measurements). When compared with subjects with low (0%-<20%) prescription-refill percentages, subjects at an elevated risk of harboring drug-resistance mutations were those with relatively high but imperfect prescription-refill percentages (80%-<90%; multivariate HR, 4.15; P<.001) and those with essentially perfect (>/=95%) refill percentages but with 2 plasma drug concentrations below the steady-state trough concentration minus 1 standard deviation (multivariate HR, 4.57; P<.001). Initial use of nonnucleoside reverse-transcriptase inhibitor-based HAART was significantly associated with multiclass drug resistance (multivariate HR, 1.84; P=.001). CONCLUSION: High baseline pVLs and substantial but imperfect levels of adherence were major predictors of antiretroviral resistance.

Antiretroviral resistance among HIV-infected persons who have died in British Columbia, in the era of modern antiretroviral therapy.

BACKGROUND: The prevalence of antiretroviral resistance among persons enrolled in the centralized HIV/AIDS Drug Treatment Program in British Columbia, Canada, who had died between July 1997 and December 2001, was investigated, to determine the degree to which antiretroviral resistance contributed to mortality. METHODS: During this period, 637 deaths had occurred. The last plasma sample obtained during therapy was genotyped retrospectively for treated individuals who had died of a nonaccidental cause. Samples with plasma human immunodeficiency virus (HIV) loads <500 copies/mL were not genotyped. Drug resistance among 1220 living HIV-infected persons who had experienced virologic therapy failure during the study period also was examined. RESULTS: Of 554 individuals who had died of nonaccidental causes, 58 (10.4%) were antiretroviral naive, and 99 (17.9%) had very brief exposure to antiretroviral therapy (median, 2 months). The majority of isolates from the remaining 397 individuals harbored either no major resistance mutations or represented samples with plasma HIV suppression of <500 copies/mL. Resistance to >/=1, >/=2, or 3 drug classes was observed in 76%, 42%, and 11% of individuals, respectively, in the group of 1220 living individuals experiencing virologic therapy failure, compared with only 44%, 23%, and 5% of individuals, respectively, who had died (P<.001). CONCLUSION: Only a relatively low prevalence of multidrug resistance was observed in this cohort, indicating that the exhaustion of treatment options because of drug resistance was not a significant contributor to mortality.