HIV-1 Coreceptor Usage Assessment by Ultra-Deep Pyrosequencing and Response to Maraviroc

Background

Maraviroc is an HIV entry inhibitor that alters the conformation of CCR5 and is poorly efficient in patients infected by viruses that use CXCR4 as an entry coreceptor. The goal of this study was to assess the capacity of ultra-deep pyrosequencing (UDPS) and different data analysis approaches to characterize HIV tropism at baseline and predict the therapeutic outcome on maraviroc treatment.

Methods

113 patients with detectable HIV-1 RNA on HAART were treated with maraviroc. The virological response was assessed at months 1, 3 and 6. The sequence of the HIV V3 loop was determined at baseline and prediction of maraviroc response by different software and interpretation algorithms was analyzed.

Results

UDPS followed by analysis with the Pyrotrop software or geno2pheno algorithm provided better prediction of the response to maraviroc than Sanger sequencing. We also found that the H34Y/S substitution in the V3 loop was the strongest individual predictor of maraviroc response, stronger than substitutions at positions 11 or 25 classically used in interpretation algorithms.

Conclusions

UDPS is a powerful tool that can be used with confidence to predict maraviroc response in HIV-1-infected patients. Improvement of the predictive value of interpretation algorithms is possible and our results suggest that adding the H34S/Y substitution would substantially improve the performance of the 11/25/charge rule.     

Tracking Changes in HIV-1 Envelope Quasispecies Using DNA Heteroduplex Analysis

A DNA heteroduplex tracking assay (HTA) using single-stranded probes is described. This assay provides a rapid means of resolving genetic variants coamplified by PCR and of measuring the level of particular variants in complex populations. To confidently detect minor quasispecies changes, the importance of maximizing template input into nested PCR (nPCR) and of duplicating nPCR and HTA to ensure correct population sampling is highlighted. The sensitivity of detection of rare variants within a genetically mixed population using single-stranded DNA probes is shown to be 1:500. The effects of nucleotide substitution at different locations on heteroduplex electrophoretic mobility are used to illustrate the limits of HTA for mutation detection. This simple assay may be used to track the evolution of HIV as well as to address issues of contamination and transmission.     

Ultra-Deep Pyrosequencing (UDPS) Data Treatment to Study Amplicon HCV Minor Variants

We have investigated the reliability and reproducibility of HCV viral quasispecies quantification by ultra-deep pyrosequencing (UDPS) methods. Our study has been divided in two parts. First of all, by UDPS sequencing of clone mixes samples we have established the global noise level of UDPS and fine tuned a data treatment workflow previously optimized for HBV sequence analysis. Secondly, we have studied the reproducibility of the methodology by comparing 5 amplicons from two patient samples on three massive sequencing platforms (FLX+, FLX and Junior) after applying the error filters developed from the clonal/control study. After noise filtering the UDPS results, the three replicates showed the same 12 polymorphic sites above 0.7%, with a mean CV of 4.86%. Two polymorphic sites below 0.6% were identified by two replicates and one replicate respectively. A total of 25, 23 and 26 haplotypes were detected by GS-Junior, GS-FLX and GS-FLX+. The observed CVs for the normalized Shannon entropy (Sn), the mutation frequency (Mf), and the nucleotidic diversity (Pi) were 1.46%, 3.96% and 3.78%. The mean absolute difference in the two patients (5 amplicons each), in the GS-FLX and GS-FLX+, were 1.46%, 3.96% and 3.78% for Sn, Mf and Pi. No false polymorphic site was observed above 0.5%.Our results indicate that UDPS is an optimal alternative to molecular cloning for quantitative study of HCV viral quasispecies populations, both in complexity and composition. We propose an UDPS data treatment workflow for amplicons from the RNA viral quasispecies which, at a sequencing depth of at least 10,000 reads per strand, enables to obtain sequences and frequencies of consensus haplotypes above 0.5% abundance with no erroneous mutations, with high confidence, resistant mutants as minor variants at the level of 1%, with high confidence that variants are not missed, and highly confident measures of quasispecies complexity.     

Transmission of Single HIV-1 Genomes and Dynamics of Early Immune Escape Revealed by Ultra-Deep Sequencing

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3–6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses – using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2–7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.     

HIV-1 Quasispecies Delineation by Tag Linkage Deep Sequencing

Trade-offs between throughput, read length, and error rates in high-throughput sequencing limit certain applications such as monitoring viral quasispecies. Here, we describe a molecular-based tag linkage method that allows assemblage of short sequence reads into long DNA fragments. It enables haplotype phasing with high accuracy and sensitivity to interrogate individual viral sequences in a quasispecies. This approach is demonstrated to deduce ∼2000 unique 1.3 kb viral sequences from HIV-1 quasispecies in vivo and after passaging ex vivo with a detection limit of ∼0.005% to ∼0.001%. Reproducibility of the method is validated quantitatively and qualitatively by a technical replicate. This approach can improve monitoring of the genetic architecture and evolution dynamics in any quasispecies population.     

Improved Detection of Rare HIV-1 Variants using 454 Pyrosequencing

454 pyrosequencing, a massively parallel sequencing (MPS) technology, is often used to study HIV genetic variation. However, the substantial mismatch error rate of the PCR required to prepare HIV-containing samples for pyrosequencing has limited the detection of rare variants within viral populations to those present above ~1%. To improve detection of rare variants, we varied PCR enzymes and conditions to identify those that combined high sensitivity with a low error rate. Substitution errors were found to vary up to 3-fold between the different enzymes tested. The sensitivity of each enzyme, which impacts the number of templates amplified for pyrosequencing, was shown to vary, although not consistently across genes and different samples. We also describe an amplicon-based method to improve the consistency of read coverage over stretches of the HIV-1 genome. Twenty-two primers were designed to amplify 11 overlapping amplicons in the HIV-1 clade B gag-pol and env gp120 coding regions to encompass 4.7 kb of the viral genome per sample at sensitivities as low as 0.01-0.2%.     

Dynamic Correlation between Intrahost HIV-1 Quasispecies Evolution and Disease Progression

Quantifying the dynamics of intrahost HIV-1 sequence evolution is one means of uncovering information about the interaction between HIV-1 and the host immune system. In the chronic phase of infection, common dynamics of sequence divergence and diversity have been reported. We developed an HIV-1 sequence evolution model that simulated the effects of mutation and fitness of sequence variants. The amount of evolution was described by the distance from the founder strain, and fitness was described by the number of offspring a parent sequence produces. Analysis of the model suggested that the previously observed saturation of divergence and decrease of diversity in later stages of infection can be explained by a decrease in the proportion of offspring that are mutants as the distance from the founder strain increases rather than due to an increase of viral fitness. The prediction of the model was examined by performing phylogenetic analysis to estimate the change in the rate of evolution during infection. In agreement with our modeling, in 13 out of 15 patients (followed for 3–12 years) we found that the rate of intrahost HIV-1 evolution was not constant but rather slowed down at a rate correlated with the rate of CD4+ T-cell decline. The correlation between the dynamics of the evolutionary rate and the rate of CD4+ T-cell decline, coupled with our HIV-1 sequence evolution model, explains previously conflicting observations of the relationships between the rate of HIV-1 quasispecies evolution and disease progression.     

HIV-1 Tropism Dynamics and Phylogenetic Analysis from Longitudinal Ultra-Deep Sequencing Data of CCR5- and CXCR4-Using Variants

Objective

Coreceptor switch from CCR5 to CXCR4 is associated with HIV disease progression. The molecular and evolutionary mechanisms underlying the CCR5 to CXCR4 switch are the focus of intense recent research. We studied the HIV-1 tropism dynamics in relation to coreceptor usage, the nature of quasispecies from ultra deep sequencing (UDPS) data and their phylogenetic relationships.

Methods

Here, we characterized C2-V3-C3 sequences of HIV obtained from 19 patients followed up for 54 to 114 months using UDPS, with further genotyping and phylogenetic analysis for coreceptor usage. HIV quasispecies diversity and variability as well as HIV plasma viral load were measured longitudinally and their relationship with the HIV coreceptor usage was analyzed. The longitudinal UDPS data were submitted to phylogenetic analysis and sampling times and coreceptor usage were mapped onto the trees obtained.

Results

Although a temporal viral genetic structuring was evident, the persistence of several viral lineages evolving independently along the infection was statistically supported, indicating a complex scenario for the evolution of viral quasispecies. HIV X4-using variants were present in most of our patients, exhibiting a dissimilar inter- and intra-patient predominance as the component of quasispecies even on antiretroviral therapy. The viral populations from some of the patients studied displayed evidences of the evolution of X4 variants through fitness valleys, whereas for other patients the data favored a gradual mode of emergence.

Conclusions

CXCR4 usage can emerge independently, in multiple lineages, along the course of HIV infection. The mode of emergence, i.e. gradual or through fitness valleys seems to depend on both virus and patient factors. Furthermore, our analyses suggest that, besides becoming dominant after population-level switches, minor proportions of X4 viruses might exist along the infection, perhaps even at early stages of it. The fate of these minor variants might depend on both viral and host factors.     

PCR-Induced Transitions Are the Major Source of Error in Cleaned Ultra-Deep Pyrosequencing Data

Background

Ultra-deep pyrosequencing (UDPS) is used to identify rare sequence variants. The sequence depth is influenced by several factors including the error frequency of PCR and UDPS. This study investigated the characteristics and source of errors in raw and cleaned UDPS data.

Results

UDPS of a 167-nucleotide fragment of the HIV-1 SG3Δenv plasmid was performed on the Roche/454 platform. The plasmid was diluted to one copy, PCR amplified and subjected to bidirectional UDPS on three occasions. The dataset consisted of 47,693 UDPS reads. Raw UDPS data had an average error frequency of 0.30% per nucleotide site. Most errors were insertions and deletions in homopolymeric regions. We used a cleaning strategy that removed almost all indel errors, but had little effect on substitution errors, which reduced the error frequency to 0.056% per nucleotide. In cleaned data the error frequency was similar in homopolymeric and non-homopolymeric regions, but varied considerably across sites. These site-specific error frequencies were moderately, but still significantly, correlated between runs (r = 0.15–0.65) and between forward and reverse sequencing directions within runs (r = 0.33–0.65). Furthermore, transition errors were 48-times more common than transversion errors (0.052% vs. 0.001%; p<0.0001). Collectively the results indicate that a considerable proportion of the sequencing errors that remained after data cleaning were generated during the PCR that preceded UDPS.

Conclusions

A majority of the sequencing errors that remained after data cleaning were introduced by PCR prior to sequencing, which means that they will be independent of platform used for next-generation sequencing. The transition vs. transversion error bias in cleaned UDPS data will influence the detection limits of rare mutations and sequence variants.     

Reconciling Longitudinal Naive T-Cell and TREC Dynamics during HIV-1 Infection

Naive T cells in untreated HIV-1 infected individuals have a reduced T-cell receptor excision circle (TREC) content. Previous mathematical models have suggested that this is due to increased naive T-cell division. It remains unclear, however, how reduced naive TREC contents can be reconciled with a gradual loss of naive T cells in HIV-1 infection. We performed longitudinal analyses in humans before and after HIV-1 seroconversion, and used a mathematical model to investigate which processes could explain the observed changes in naive T-cell numbers and TRECs during untreated HIV-1 disease progression. Both CD4⁺ and CD8⁺ naive T-cell TREC contents declined biphasically, with a rapid loss during the first year and a much slower loss during the chronic phase of infection. While naive CD8⁺ T-cell numbers hardly changed during follow-up, naive CD4⁺ T-cell counts continually declined. We show that a fine balance between increased T-cell division and loss in the peripheral naive T-cell pool can explain the observed short- and long-term changes in TRECs and naive T-cell numbers, especially if T-cell turnover during the acute phase is more increased than during the chronic phase of infection. Loss of thymic output, on the other hand, does not help to explain the biphasic loss of TRECs in HIV infection. The observed longitudinal changes in TRECs and naive T-cell numbers in HIV-infected individuals are most likely explained by a tight balance between increased T-cell division and death, suggesting that these changes are intrinsically linked in HIV infection.     

Dynamics of HIV-1 Assembly and Release

Assembly and release of human immunodeficiency virus (HIV) occur at the plasma membrane of infected cells and are driven by the Gag polyprotein. Previous studies analyzed viral morphogenesis using biochemical methods and static images, while dynamic and kinetic information has been lacking until very recently. Using a combination of wide-field and total internal reflection fluorescence microscopy, we have investigated the assembly and release of fluorescently labeled HIV-1 at the plasma membrane of living cells with high time resolution. Gag assembled into discrete clusters corresponding to single virions. Formation of multiple particles from the same site was rarely observed. Using a photoconvertible fluorescent protein fused to Gag, we determined that assembly was nucleated preferentially by Gag molecules that had recently attached to the plasma membrane or arrived directly from the cytosol. Both membrane-bound and cytosol derived Gag polyproteins contributed to the growing bud. After their initial appearance, assembly sites accumulated at the plasma membrane of individual cells over 1–2 hours. Assembly kinetics were rapid: the number of Gag molecules at a budding site increased, following a saturating exponential with a rate constant of ∼5×10⁻³ s⁻¹, corresponding to 8–9 min for 90% completion of assembly for a single virion. Release of extracellular particles was observed at ∼1,500±700 s after the onset of assembly. The ability of the virus to recruit components of the cellular ESCRT machinery or to undergo proteolytic maturation, or the absence of Vpu did not significantly alter the assembly kinetics.     

Antiviral properties of combination peptides of HIV-1 Rev NLS and NES

Nuclear translocation signal has been identified as a mediator of protein shuttling between nuclear and cytoplasm. Here we report that the combination of peptides from nuclear localization signal (NLS) and nuclear export signal (NES) of HIV-1 Rev have an antiviral activity against the Herpes virus of turkey and Marek's disease virus serotype 1.     

Exploring the Dynamic Core Microbiome of Plaque Microbiota during Head-and-Neck Radiotherapy Using Pyrosequencing

Radiotherapy is the primary treatment modality used for patients with head-and-neck cancers, but inevitably causes microorganism-related oral complications. This study aims to explore the dynamic core microbiome of oral microbiota in supragingival plaque during the course of head-and-neck radiotherapy. Eight subjects aged 26 to 70 were recruited. Dental plaque samples were collected (over seven sampling time points for each patient) before and during radiotherapy. The V1–V3 hypervariable regions of bacterial 16S rRNA genes were amplified, and the high-throughput pyrosequencing was performed. A total of 140 genera belonging to 13 phyla were found. Four phyla (Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria) and 11 genera (Streptococcus, Actinomyces, Veillonella, Capnocytophaga, Derxia, Neisseria, Rothia, Prevotella, Granulicatella, Luteococcus, and Gemella) were found in all subjects, supporting the concept of a core microbiome. Temporal variation of these major cores in relative abundance were observed, as well as a negative correlation between the number of OTUs and radiation dose. Moreover, an optimized conceptual framework was proposed for defining a dynamic core microbiome in extreme conditions such as radiotherapy. This study presents a theoretical foundation for exploring a core microbiome of communities from time series data, and may help predict community responses to perturbation as caused by exposure to ionizing radiation.