Luke C. Swenson

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.

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≤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 ∼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≤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.

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.

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.

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.

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 PSSMNSI/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.

Performance of HIV-1 Drug Resistance Testing at Low-Level Viremia and Its Ability to Predict Future Virologic Outcomes and Viral Evolution in Treatment-Naive Individuals

BACKGROUND Low-level viremia (LLV; human immunodeficiency virus [HIV-1] RNA 50-999 copies/mL) occurs frequently in patients receiving antiretroviral therapy (ART), but there are few or no data available demonstrating that HIV-1 drug resistance testing at a plasma viral load (pVL) <1000 copies/mL provides potentially clinically useful information. Here, we assess the ability to perform resistance testing by genotyping at LLV and whether it is predictive of future virologic outcomes in patients beginning ART. METHODS Resistance testing by genotyping at LLV was attempted on 4915 plasma samples from 2492 patients. A subset of previously ART-naive patients was analyzed who achieved undetectable pVL and subsequently rebounded with LLV (n = 212). A genotypic sensitivity score (GSS) was calculated based on therapy and resistance testing results by genotyping, and stratified according to number of active drugs. RESULTS Eighty-eight percent of LLV resistance assays produced useable sequences, with higher success at higher pVL. Overall, 16 of 212 (8%) patients had pretherapy resistance. Thirty-eight of 196 (19%) patients without pretherapy resistance evolved resistance to 1 or more drug classes, primarily the nucleoside reverse transcriptase (14%) and/or nonnucleoside reverse transcriptase (9%) inhibitors. Patients with resistance at LLV (GSS <3) had a 2.1-fold higher risk of virologic failure (95% confidence interval, 1.2- to 3.7-fold) than those without resistance (P = .007). Progressively lower GSS scores at LLV were associated with a higher increase in pVL over time (P < .001). Acquisition of additional resistance mutations to a new class of antiretroviral drugs during LLV was not found in a subset of patients. CONCLUSIONS Routine HIV-1 genotyping of LLV samples can be performed with a reasonably high success rate, and the results appear predictive of future virologic outcomes.

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.

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.

Comparative Performances of HIV-1 RNA Load Assays at Low Viral Load Levels: Results of an International Collaboration

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