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.

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.

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.

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.

Evidence of differential HLA class I-mediated viral evolution in functional and accessory/regulatory genes of HIV-1.

Despite the formidable mutational capacity and sequence diversity of HIV-1, evidence suggests that viral evolution in response to specific selective pressures follows generally predictable mutational pathways. Population-based analyses of clinically derived HIV sequences may be used to identify immune escape mutations in viral genes; however, prior attempts to identify such mutations have been complicated by the inability to discriminate active immune selection from virus founder effects. Furthermore, the association between mutations arising under in vivo immune selection and disease progression for highly variable pathogens such as HIV-1 remains incompletely understood. We applied a viral lineage-corrected analytical method to investigate HLA class I-associated sequence imprinting in HIV protease, reverse transcriptase (RT), Vpr, and Nef in a large cohort of chronically infected, antiretrovirally naïve individuals. A total of 478 unique HLA-associated polymorphisms were observed and organized into a series of "escape maps," which identify known and putative cytotoxic T lymphocyte (CTL) epitopes under selection pressure in vivo. Our data indicate that pathways to immune escape are predictable based on host HLA class I profile, and that epitope anchor residues are not the preferred sites of CTL escape. Results reveal differential contributions of immune imprinting to viral gene diversity, with Nef exhibiting far greater evidence for HLA class I-mediated selection compared to other genes. Moreover, these data reveal a significant, dose-dependent inverse correlation between HLA-associated polymorphisms and HIV disease stage as estimated by CD4(+) T cell count. Identification of specific sites and patterns of HLA-associated polymorphisms across HIV protease, RT, Vpr, and Nef illuminates regions of the genes encoding these products under active immune selection pressure in vivo. The high density of HLA-associated polymorphisms in Nef compared to other genes investigated indicates differential HLA class I-driven evolution in different viral genes. The relationship between HLA class I-associated polymorphisms and lower CD4(+) cell count suggests that immune escape correlates with disease status, supporting an essential role of maintenance of effective CTL responses in immune control of HIV-1. The design of preventative and therapeutic CTL-based vaccine approaches could incorporate information on predictable escape pathways.

HIV-1 drug resistance surveillance using dried whole blood spots.

BACKGROUND: Field-friendly methods for HIV drug resistance (HIVDR) surveillance in resource-limited regions are urgently needed. Despite evidence that dried blood spots (DBS) are suitable for HIV serology, viral load and CD4+ T-cell enumeration, no study has evaluated DBS for HIVDR genotyping. We assessed the feasibility of genotyping HIV-1 from field-collected DBS stored under challenging environmental conditions. METHODS: We prospectively collected specimens from newly diagnosed, treatment-naive HIV-positive subjects in Mexico. Whole blood was spotted onto filter cards, air dried at ambient temperature and stored with desiccant at 37 degrees C and 85% humidity for 3 months. Genotypes obtained from DBS-extracted nucleic acids using an in-house nested reverse transcription-PCR method were compared to genotypes derived from matched plasma. RESULTS: Genotypes from 103 phylogenetically matched plasma and DBS were compared. In total, 90.1% of all DBS specimens could be amplified in either the region of HIV protease or the region of reverse transcriptase. Failure to amplify from DBS did not correlate with low plasma viral loads. Between paired specimens, the median nucleotide similarity was 99.95%. In the nine specimens with drug resistance mutations, all differences between pairs were partial discordances. Mutations identified in plasma were found in the majority of replicate DBS amplifications. CONCLUSION: The results suggest that genotypes obtained from DBS are equivalent to those from plasma. DBS are a promising public health tool for HIVDR surveillance of treatment-naive subjects, especially in regions where specimens might be exposed to severe environmental conditions and where logistical difficulties could prevent timely specimen processing. More studies are needed to validate DBS for patient monitoring.

A 21-base pair insertion/duplication at codon 69 of the HIV type 1 reverse transcriptase in a patient undergoing multiple nucleoside therapy.

We describe the selection of a previously unreported 21-base pair insertion following codon 69 of the HIV-1 reverse transcriptase (RT) from a patient undergoing multiple nucleoside analogue therapy. This insertion was a direct duplication of the preceding 21 bases of HIV-RT, and was selected in a background of NRTI-resistance mutations including substitutions at RT codons 41, 67, 184, 210, and 215. Longitudinal genotypic and phenotypic resistance tests performed before and after selection of the insertion suggested that the insertion conferred an additional decrease in susceptibility to some nucleoside analogues, most notably didanosine, stavudine, abacavir, and tenofovir. However, phenotypic analysis of an insertion-containing site-directed mutant constructed in an HIV-1 HXB2 background revealed no direct association between the 21-base pair insertion and decreased susceptibility to NRTIs, suggesting that the insert requires the context of the patients' virus in order to confer resistance. These observations may offer new insight into the relative contribution of HIV-RT codon 69 insertion mutations to antiretroviral resistance.

Prevalence and clinical impacts of HIV-1 intersubtype recombinants in Uganda revealed by near-full-genome population and deep sequencing approaches.

OBJECTIVES: HIV-1 subtypes A1 and D cocirculate in a rural community in Mbarara, Uganda. This study examines HIV-1 intersubtype recombination in this community under a full-genome sequencing context. We aim to estimate prevalence, examine time trends, and test for clinical correlates and outcomes associated with intersubtype recombinants. METHODS: Near-full-genome HIV-1 Sanger sequence data were collected from plasma samples of 504 treatment-naïve individuals, who then received protease inhibitor or nonnucleoside reverse transcriptase inhibitor-containing regimens and were monitored for up to 7.5 years. Subtypes were inferred by Los Alamos Recombinant Identification Program (RIP) 3.0 and compared with Sanger/REGA and MiSeq/RIP. 'Nonrecombinants' and 'recombinants' infections were compared in terms of pretherapy viral load, CD4 cell count, posttherapy time to virologic suppression, virologic rebound, first CD4 rise above baseline and sustained CD4 recovery. RESULTS: Prevalence of intersubtype recombinants varied depending on the genomic region examined: gag (15%), prrt (11%), int (8%), vif (10%), vpr (2%), vpu (9%), GP120 (8%), GP41 (18%), and nef (4%). Of the 200 patients with near-full-genome data, prevalence of intersubtype recombination was 46%; the most frequently observed recombinant was A1-D (25%). Sanger/REGA and MiSeq/RIP yielded generally consistent results. Phylogenetic tree revealed most recombinants did not share common ancestors. No temporal trend was observed (all P > 0.1). Subsequent subtype switches were detected in 27 of 143 (19%) study participants with follow-up sequences. Nonrecombinant versus recombinants infections were not significantly different in any pre nor posttherapy clinical correlates examined (all P > 0.2). CONCLUSION: Intersubtype recombination was highly prevalent (46%) in Uganda if the entire HIV genome was considered, but was neither associated with clinical correlates nor therapy outcomes.

Utilization of a liver allograft from a hepatitis B surface antigen positive donor.

With today's donor organ shortage, enhanced efforts must be made to utilize organs that previously would have been declined. We report a 26-year-old man with chronic hepatitis B virus (HBV) and hepatitis C virus (HCV) coinfection who received a liver transplant from an HBsAg-positive donor. HBV viremia (6,281,185 copies/ml) was seen early posttransplant despite lamivudine prophylaxis, but became negative with addition of adefovir. Virologic analysis revealed predominantly donor HBV strain immediately posttransplant. At 5 months there was an elevation of liver enzymes accompanied by histologic evidence of hepatitis. At this time, HCV-RNA was positive but HBV DNA was undetectable. Treatment with pegylated interferon and ribavirin resulted in sustained clearance of HCV RNA. Two years posttransplant, the patient has normal liver biochemistry and HCV and HBV viral load are undetectable with persistence of HBsAg. Our experience suggests that with effective antiviral therapy, the use of HBsAg seropositive donors is feasible in selected circumstances.

HIV VprR77Q mutation does not influence clinical response of individuals initiating highly active antiretroviral therapy.

VprR77Q has been associated with long-term nonprogressive (LTNP) HIV infection. We wished to investigate the prevalence, clinical correlates, and effect on treatment response of VprR77Q in a cohort of antiretroviral- naïve individuals initiating highly active antiretroviral therapy (HAART). Baseline plasma samples from 728 subjects were genotyped using RT-PCR and direct DNA sequencing. Cox proportional hazards regression was used to model the effects of VprR77Q on virologic and immunologic responses, and survival following initiation of HAART, over a median 4.5 years follow-up. We found that 308 subjects (42.3%) harbored VprR77Q alone or in combination with another amino acid, while 420 (57.7%) harbored an amino acid other than Q. A cross-sectional analysis found no correlation between R77Q and baseline plasma viral load (pVL), CD4 count, diagnosis of AIDS, or sociodemographic characteristics including age, gender, and history of injection drug use (p > 0.1). In multivariate analyses, no significant associations between VprR77Q and initial pVL and CD4 responses to HAART or survival following initiation of treatment were observed (p > 0.1). The high prevalence and the lack of association with pretherapy clinical parameters in this cohort argue against an association of R77Q with LTNP status. These results do not support an association between R77Q and HAART response.

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.

Rare mutations at codon 103 of HIV-1 reverse transcriptase can confer resistance to non-nucleoside reverse transcriptase inhibitors.

BACKGROUND: The K103N mutation in HIV-1 reverse transcriptase (RT) confers high-level resistance to current non-nucleoside reverse transcriptase inhibitors (NNRTI). The prevalence and resistance profile of HIV-1 with other substitutions at RT codon 103 is less well documented. METHODS: K103 substitutions among over 70,000 clinical samples submitted for routine antiretroviral resistance testing at two independent centres were examined. Phenotypic resistance profiles of isolates harboring rare K103 variants in the absence of known NNRTI-associated resistance mutations were retrieved from Virco's correlative genotype/phenotype database. Genotyped samples with known treatment histories were retrieved from the British Columbia Centre for Excellence in HIV/AIDS database. Site-directed mutants containing K103 variants were constructed and phenotyped. RESULTS: K103N, R and S were observed in 29, 1.8, and 0.9% of Virco isolates and in 16, 1.5 and 0.4% of British Columbia isolates. K103T/Q/H substitutions were observed only rarely (<0.2%). The prevalence of unusual codon 103 substitutions remained stable over 5 years, except K103S, which increased over fourfold in both datasets. K103R/Q-containing clinical isolates remained phenotypically susceptible to NNRTI, whereas K103S/T/H-containing isolates showed over 10-fold decreased NNRTI susceptibility. Among patients with a known treatment history, K103S/T/H were observed primarily in individuals failing NNRTI-containing regimens. Site-directed mutants confirmed decreased susceptibility to NNRTI in K103S/T/H-containing recombinants. CONCLUSION: Variants at HIV RT codon 103 other than K103N are observed relatively rarely in clinical isolates, but K103 S, T and H confer decreased susceptibility to NNRTI. These data are relevant for interpretive genotype algorithms and in the design of assays specific to RT codon 103 mutations.

No association between GB virus-C viremia and virological or immunological failure after starting initial antiretroviral therapy.

INTRODUCTION: Co-infection with GBV-C ('Hepatitis G' virus) appears to be associated with slower disease progression in HIV-infected, untreated individuals. We wished to determine whether detection of GBV-C RNA was associated with differential response to HIV therapy in a population-based cohort of 461 individuals initiating antiretroviral therapy between June 1996 and August 1998, in British Columbia, Canada. METHODS: The presence of GBV-C RNA in plasma was identified by nested RT-PCR, using detection of HIV RNA as a positive control. Time to virological success [achieving HIV plasma viral load (pVL) < or = 500 copies/ml], virological failure (subsequent confirmed pVL > 500 copies/ml) and immunological failure (confirmed CD4 cell count below baseline) were assessed by Kaplan-Meier methods and Cox proportional hazard regression. RESULTS: Of the 441 individuals for whom results were available, 90 (20.4%) had detectable plasma GBV-C RNA. GBV-C RNA was significantly associated with a lower HIV pVL at baseline (P = 0.004). In univariate and multivariate Cox models, GBV-C RNA positive and negative individuals did not differ with respect to time to virological success [risk ratio (RR), 0.98; 95% confidence interval (CI), 0.75-1.27], time to virological failure (RR, 1.10; 95% CI, 0.74-1.65), or time to immunological failure (RR, 1.09; 95% CI, 0.73-1.63). There was no correlation between detection of GBV-C RNA and mutations in the human chemokine receptors CCR5 and CX CR1, or HIV viral tropism as predicted by the HIV envelope sequence (P > 0.1). CONCLUSION: GBV-C viremia is relatively common in individuals seeking treatment for HIV infection; however, it does not appear to have any effect on initial antiretroviral therapy response.

Prevalence and response to antiretroviral therapy of non-B subtypes of HIV in antiretroviral-naive individuals in British Columbia.

In North America, the B subtype of the major group (M) of HIV-1 predominates. Phylogenetic analysis of HIV reverse transcriptase and protease sequences isolated from 479 therapy-naive patients, first seeking treatment in British Columbia between June 1997 and August 1998, revealed a prevalence of 4.4% non-B virus. A range of different subtypes was identified, including one subtype A, 11 C, two D, five CRF01_AE, and one sample that could not be reliably subtyped. Baseline CD4 courts were significantly lower in individuals harbouring the non-B subtypes (P = 0.02), but baseline viral loads were similar (P = 0.80). In this study, individuals infected with non-B variants did not have a significantly different virological response to therapy after up to 18 months.

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.

Limited evolution of inferred HIV-1 tropism while viremia is undetectable during standard HAART therapy.

BACKGROUND: HIV patients on suppressive antiretroviral therapy have undetectable viremia making it impossible to screen plasma HIV tropism if regimen change is required during suppression. We investigated the prevalence and predictors of tropism switch from CCR5-using ("R5") to non-CCR5-using ("non-R5") before and after viral suppression in the initially therapy-naïve HOMER cohort from British Columbia, Canada. METHODS: We compared pre-therapy and post-suppression viral genotypic tropism in patients who initiated on PI/NNRTI-based antiretroviral regimens between 1996-1999 (n = 462). Virologic suppression was defined as having two consecutive viral loads of <500 copies/mL, which was the sensitivity limit of most viral load assays at the time. Viral tropism was inferred by V3-loop-population-sequencing and geno2pheno[coreceptor] with cutoff at 5.75% false positive rate (FPR). RESULTS: When virologic suppression was defined as two-consecutive viral loads <500 copies/mL, 34 (9%) of the 397 patients with pre-therapy R5-virus switched to non-R5 at viral load rebound after a median of 19 months (IQR 8-41 months) of undetectable viremia. Duration of viral load suppression was not a predictor of switch, but lower CD4 count during suppression (median 400 versus 250 cells/mL) and an increased prevalence of pre-therapy non-R5 HIV by "deep" sequencing (median 0.2% versus 3.2%) were independently associated with switch (p = 0.03 and p<0.0001, respectively). CONCLUSION: R5-to-non-R5 tropism switches in plasma virus after undetectable viremia were relatively rare events especially among patients with higher CD4 counts during virologic suppression. Our study supports the use of pre-suppression tropism results if maraviroc is being considered during virologic suppression in this subgroup of patients.

HLA class I sequence-based typing using DNA recovered from frozen plasma.

We describe a rapid, reliable and cost-effective method for intermediate-to-high-resolution sequence-based HLA class I typing using frozen plasma as a source of genomic DNA. The plasma samples investigated had a median age of 8.5 years. Total nucleic acids were isolated from matched frozen PBMC (~2.5 million) and plasma (500 µl) samples from a panel of 25 individuals using commercial silica-based kits. Extractions yielded median [IQR] nucleic acid concentrations of 85.7 [47.0-130.0]ng/µl and 2.2 [1.7-2.6]ng/µl from PBMC and plasma, respectively. Following extraction, ~1000 base pair regions spanning exons 2 and 3 of HLA-A, -B and -C were amplified independently via nested PCR using universal, locus-specific primers and sequenced directly. Chromatogram analysis was performed using commercial DNA sequence analysis software and allele interpretation was performed using a free web-based tool. HLA-A, -B and -C amplification rates were 100% and chromatograms were of uniformly high quality with clearly distinguishable mixed bases regardless of DNA source. Concordance between PBMC and plasma-derived HLA types was 100% at the allele and protein levels. At the nucleotide level, a single partially discordant base (resulting from a failure to call both peaks in a mixed base) was observed out of >46,975 bases sequenced (>99.9% concordance). This protocol has previously been used to perform HLA class I typing from a variety of genomic DNA sources including PBMC, whole blood, granulocyte pellets and serum, from specimens up to 30 years old. This method provides comparable specificity to conventional sequence-based approaches and could be applied in situations where cell samples are unavailable or DNA quantities are limiting.

Population-based sequencing of the V3-loop can predict the virological response to maraviroc in treatment-naive patients of the MERIT trial.

BACKGROUND: MERIT was a randomized trial comparing maraviroc (MVC) + Combivir versus efavirenz (EFV) + Combivir in drug-naive patients screened as having R5 HIV-1 by the original Trofile assay (OTA). We retrospectively evaluated treatment response after rescreening for viral tropism using population-based V3-loop sequencing. METHODS: HIV env V3-loop was amplified in triplicate using reverse transcriptase-polymerase chain reaction from stored screening plasma and sequenced. Automated base calling was performed using custom software (RECall) and tropism inferred by geno2pheno (5.75% false-positive rate). Tropism results by genotype were compared with those of OTA and Enhanced Sensitivity Trofile assay (ESTA), where all results were available (n = 876). RESULTS: Approximately 8% of patients screened as having R5 virus by OTA were classified as having non-R5 virus by V3-loop genotyping. These patients were less likely to have early or sustained week-48 treatment response to MVC, but not EFV. When restricted to patients with R5 virus by genotype, reanalysis of the primary study endpoint (plasma viral load <50 copies/mL at week 48) showed noninferiority of MVC twice daily to EFV (67% vs. 68%). Rescreening by genotype and ESTA had 84% concordance; patients receiving MVC twice daily rescreened as having R5 virus had greater than 1 log10 copies per milliliter decrease in viral load over those rescreened as having non-R5 virus. Where genotype and ESTA screening results were discordant outcomes were similar. CONCLUSIONS: The exclusion of ∼8% of patients with CXCR4-using virus by population-based sequencing would likely have resulted in noninferior responses in the MVC twice-daily and EFV arms. Rescreening by ESTA and population-based sequencing predicted similar virological response.