Bone marrow CD34+ progenitor cells may harbour HIV-DNA even in successfully treated patients

AbstractThe issue about bone marrow hematopoietic progenitor cells harbouring HIV-DNA in infected patients is still under scrutiny. We studied nine HIV-infected individuals undergoing bone marrow aspiration for diagnostic purposes. In all patients, even in those receiving successful antiretroviral therapy for several years, HIV-DNA was detected in purified CD34+ lineage-bone marrow progenitor cells. This finding, although not conclusive due to the low number of patients examined, adds further evidence that current treatment strategies may be insufficient to resolve latent infection in bone marrow CD34+ hematopoietic progenitor cells

Analysis of HIV quasispecies and virological outcome of an HIV D+/R+ kidney–liver transplantation

Introduction: Transplantation among HIV positive patients may be a valuable therapeutic intervention. This study involves an HIV D+/R+ kidney–liver transplantation, where PBMC-associated HIV quasispecies were analyzed in donor and transplant recipients (TR) prior to transplantation and thereafter, together with standard viral monitoring. Methods: The donor was a 54 year of age HIV infected woman: kidney and liver recipients were two HIV infected men, aged 49 and 61. HIV quasispecies in PBMC was analyzed by ultra-deep sequencing of V3 env region. During TR follow-up, plasma HIV-1 RNA, HIV-1 DNA in PBMC, analysis of proviral integration sites and drug-resistance genotyping were performed. Other virological and immunological monitoring included CMV and EBV DNA quantification in blood and CD4 T cell counts. Results: Donor and TR were all ART-HIV suppressed at transplantation. Thereafter, TR maintained a nearly suppressed HIV-1 viremia, but HIV-1 RNA blips and the increase of proviral integration sites in PBMC attested some residual HIV replication. A transient peak in HIV-1 DNA occurred in the liver recipient. No major changes of drug-resistance genotype were detected after transplantation. CMV and EBV transient reactivations were observed only in the kidney recipient, but did not require specific treatment. CD4 counts remained stable. No intermixed quasispecies between donor and TR was observed at transplantation or thereafter. Despite signs of viral evolution in TR, HIV genetic heterogeneity did not increase over the course of the months of follow up. Conclusions: No evidence of HIV superinfection was observed in the donor nor in the recipients. The immunosuppressive treatment administrated to TR did not result in clinical relevant viral reactivations

The genotypic false positive rate determined by V3 population sequencing can predict the burden of HIV-1 CXCR4-using species detected by pyrosequencing

The false-positive rate (FPR) is a percentage-score provided by Geno2Pheno-algorithm indicating the likelihood that a V3-sequence is falsely predicted as CXCR4-using. We evaluated the correlation between FPR obtained by V3 population-sequencing and the burden of CXCR4-using variants detected by V3 ultra-deep sequencing (UDPS) and Enhanced-Sensitivity Trofile assay (ESTA)

HIV Drug Resistance Surveillance Using Pooled Pyrosequencing

BACKGROUND: Surveillance for HIV transmitted drug resistance (TDR) is performed using HIV genotype results from individual specimens. Pyrosequencing, through its massive parallel sequencing ability, can analyze large numbers of specimens simultaneously. Instead of using pyrosequencing conventionally, to sequence a population of viruses within an individual, we interrogated a single combined pool of surveillance specimens to demonstrate that it is possible to determine TDR rates in HIV protease from a population of individuals. METHODOLOGY/PRINCIPAL FINDINGS: The protease region from 96 treatment naïve, HIV+ serum specimens was genotyped using standard Sanger sequencing method. The 462 bp protease amplicons from these specimens were pooled in equimolar concentrations and re-sequenced using the GS FLX Titanium system. The nucleotide (NT) and amino acid (AA) differences from the reference sequence, along with TDR mutations, detected by each method were compared. In the protease sequence, there were 212 nucleotide and 81 AA differences found using conventional sequencing and 345 nucleotide and 168 AA differences using pyrosequencing. All nucleotide and amino acid polymorphisms found at frequencies >/=5% in pyrosequencing were detected using both methods with the rates of variation highly correlated. Using Sanger sequencing, two TDR mutations, M46L and I84V, were each detected as mixtures at a frequency of 1.04% (1/96). These same TDR mutations were detected by pyrosequencing with a prevalence of 0.29% and 0.34% respectively. Phylogenetic analysis established that the detected low frequency mutations arose from the same single specimens that were found to contain TDR mutations by Sanger sequencing. Multiple clinical protease DR mutations present at higher frequencies were concordantly identified using both methods. CONCLUSIONS/SIGNIFICANCE: We show that pyrosequencing pooled surveillance specimens can cost-competitively detect protease TDR mutations when compared with conventional methods. With few modifications, the method described here can be used to determine population rates of TDR in both protease and reverse transcriptase. Furthermore, this pooled pyrosequencing technique may be generalizable to other infectious agents where a survey of DR rates is required

Performance of Ultra-Deep Pyrosequencing in Analysis of HIV-1 pol Gene Variation

INTRODUCTION: Ultra-deep pyrosequencing (UDPS) has been used to detect minority variants within HIV-1 populations. Some aspects of the quality and reproducibility of UDPS have been previously evaluated, but comprehensive studies are still needed. PRINCIPAL FINDING: In this study the UDPS technology (FLX platform) was evaluated by analyzing a 120 base pair fragment of the HIV-1 pol gene from plasma samples from two patients and artificial mixtures of molecular clones. UDPS was performed using an optimized experimental protocol and an in-house data cleaning strategy. Nine samples and mixtures were analyzed and the average number of reads per sample was 19,404 (range 8,858-26,846). The two patient plasma samples were analyzed twice and quantification of viral variants was found to be highly repeatable for variants representing >0.27% of the virus population, whereas some variants representing 0.11-0.27% were detected in only one of the two UDPS runs. Bland-Altman analysis showed that a repeated measurement would have a 95% likelihood to lie approximately within ±0.5 log(10) of the initial estimate. A similar level of agreement was observed for variant frequency estimates in forward vs. reverse sequencing direction, but here the agreement was higher for common variants than for rare variants. UDPS following PCR amplification with alternative primers indicated that some variants may be incorrectly quantified due to primer-related selective amplification. Finally, the in vitro recombination rate during PCR was evaluated using artificial mixtures of clones and was found to be low. The most abundant in vitro recombinant represented 0.25% of all UDPS reads. CONCLUSION: This study demonstrates that this UDPS protocol results in low experimental noise and high repeatability, which is relevant for future research and clinical use of the UDPS technology. The low rate of in vitro recombination suggests that this UDPS system can be used to study genetic variants and mutational linkage

Sentinel mutations in standard population sequencing can predict the presence of HIV-1 reverse transcriptase major mutations detectable only by ultra-deep pyrosequencing

Objectives: This proof-of-concept study aimed to identify whether mutations considered not yet relevant for drug resistance (but located at key drug-resistance positions) can act as 'sentinels' of minority resistant variants in HIV-1 drug-naive patients. Methods: We focused our attention on three reverse transcriptase (RT) mutations (T69S, L210M and K103R) easily detected by standard population sequencing [i.e. the genotypic resistance test (GRT)]. Ultra-deep pyrosequencing (UDPS) of HIV-1 RT was performed using GS-FLX Roche, on plasma RNA from 40 drug-naive patients infected with HIV-1 subtype B without primary resistance detected by GRT. Only RT drug resistance mutations detected at >0.1% in both forward and reverse directions were considered. Associations between GRT sentinel mutations and UDPS drug resistance were assessed using Fisher's exact test. Results: UDPS detected drug resistance mutations in 18/40 drug-naive patients. Patients carrying HIV-1 strains with T69S and L210M by GRT showed a trend to greater infection by minority drug-resistant variants than control patients infected by HIV-1 without these mutations (5/10 and 7/10 versus 3/10; P not significant). No association was found for K103R by GRT. Notably, T69S and L210M(but not K103R or control viruses) were associated with GRT minority drug resistant variants with a prevalence >1% (3/10 and 4/10 versus 0/20 in K103R and controls; P=0.03 and P=0.008, respectively). Moreover, the presence of L210M or T69S viruses by GRT significantly correlated with that of minority thymidine analogue mutations by UDPS (6/20 patients carrying HIV-1 strains with T69S/L210M versus 0/20 patients carrying HIV-1 having K103R or none of these mutations; P=0.03). Conclusions: This proof-of-concept study suggests the existence of genetic markers, detectable by routine testing, potentially acting as sentinel mutations of minority drug resistance. Their identification may help in the selection of patients at high risk of resistance in reservoirs without the necessity of using UDPS

Development of a Low Bias Method for Characterizing Viral Populations Using Next Generation Sequencing Technology

Background: With an estimated 38 million people worldwide currently infected with human immunodeficiency virus (HIV), and an additional 4.1 million people becoming infected each year, it is important to understand how this virus mutates and develops resistance in order to design successful therapies. Methodology/Principal Findings: We report a novel experimental method for amplifying full-length HIV genomes without the use of sequence-specific primers for high throughput DNA sequencing, followed by assembly of full length viral genome sequences from the resulting large dataset. Illumina was chosen for sequencing due to its ability to provide greater coverage of the HIV genome compared to prior methods, allowing for more comprehensive characterization of the heterogeneity present in the HIV samples analyzed. Our novel amplification method in combination with Illumina sequencing was used to analyze two HIV populations: a homogenous HIV population based on the canonical NL4-3 strain and a heterogeneous viral population obtained from a HIV patient's infected T cells. In addition, the resulting sequence was analyzed using a new computational approach to obtain a consensus sequence and several metrics of diversity. Significance: This study demonstrates how a lower bias amplification method in combination with next generation DNA sequencing provides in-depth, complete coverage of the HIV genome, enabling a stronger characterization of the quasispecies present in a clinically relevant HIV population as well as future study of how HIV mutates in response to a selective pressure

Analysis of high-depth sequence data for studying viral diversity: a comparison of next generation sequencing platforms using Segminator II

Background: Next generation sequencing provides detailed insight into the variation present within viral populations, introducing the possibility of treatment strategies that are both reactive and predictive. Current software tools, however, need to be scaled up to accommodate for high-depth viral data sets, which are often temporally or spatially linked. In addition, due to the development of novel sequencing platforms and chemistries, each with implicit strengths and weaknesses, it will be helpful for researchers to be able to routinely compare and combine data sets from different platforms/chemistries. In particular, error associated with a specific sequencing process must be quantified so that true biological variation may be identified. Results: Segminator II was developed to allow for the efficient comparison of data sets derived from different sources. We demonstrate its usage by comparing large data sets from 12 influenza H1N1 samples sequenced on both the 454 Life Sciences and Illumina platforms, permitting quantification of platform error. For mismatches median error rates at 0.10 and 0.12%, respectively, suggested that both platforms performed similarly. For insertions and deletions median error rates within the 454 data (at 0.3 and 0.2%, respectively) were significantly higher than those within the Illumina data (0.004 and 0.006%, respectively). In agreement with previous observations these higher rates were strongly associated with homopolymeric stretches on the 454 platform. Outside of such regions both platforms had similar indel error profiles. Additionally, we apply our software to the identification of low frequency variants. Conclusion: We have demonstrated, using Segminator II, that it is possible to distinguish platform specific error from biological variation using data derived from two different platforms. We have used this approach to quantify the amount of error present within the 454 and Illumina platforms in relation to genomic location as well as location on the read. Given that next generation data is increasingly important in the analysis of drug-resistance and vaccine trials, this software will be useful to the pathogen research community. A zip file containing the source code and jar file is freely available for download from http://www.bioinf.manchester.ac.uk/segminator/

The Evolutionary Analysis of Emerging Low Frequency HIV-1 CXCR4 Using Variants through Time—An Ultra-Deep Approach

Large-scale parallel pyrosequencing produces unprecedented quantities of sequence data. However, when generated from viral populations current mapping software is inadequate for dealing with the high levels of variation present, resulting in the potential for biased data loss. In order to apply the 454 Life Sciences' pyrosequencing system to the study of viral populations, we have developed software for the processing of highly variable sequence data. Here we demonstrate our software by analyzing two temporally sampled HIV-1 intra-patient datasets from a clinical study of maraviroc. This drug binds the CCR5 coreceptor, thus preventing HIV-1 infection of the cell. The objective is to determine viral tropism (CCR5 versus CXCR4 usage) and track the evolution of minority CXCR4-using variants that may limit the response to a maraviroc-containing treatment regimen. Five time points (two prior to treatment) were available from each patient. We first quantify the effects of divergence on initial read k-mer mapping and demonstrate the importance of utilizing population-specific template sequences in relation to the analysis of next-generation sequence data. Then, in conjunction with coreceptor prediction algorithms that infer HIV tropism, our software was used to quantify the viral population structure pre- and post-treatment. In both cases, low frequency CXCR4-using variants (2.5–15%) were detected prior to treatment. Following phylogenetic inference, these variants were observed to exist as distinct lineages that were maintained through time. Our analysis, thus confirms the role of pre-existing CXCR4-using virus in the emergence of maraviroc-insensitive HIV. The software will have utility for the study of intra-host viral diversity and evolution of other fast evolving viruses, and is available from http://www.bioinf.manchester.ac.uk/segminator/

Dynamics of Envelope Evolution in Clade C SHIV-Infected Pig-Tailed Macaques during Disease Progression Analyzed by Ultra-Deep Pyrosequencing

Understanding the evolution of the human immunodeficiency virus type 1 (HIV-1) envelope during disease progression can provide tremendous insights for vaccine development, and simian-human immunodeficiency virus (SHIV) infection of non-human primate provides an ideal platform for such studies. A newly developed clade C SHIV, SHIV-1157ipd3N4, which was able to infect rhesus macaques, closely resembled primary HIV-1 in transmission and pathogenesis, was used to infect several pig-tailed macaques. One of the infected animals subsequently progressed to AIDS, whereas one remained a non-progressor. The viral envelope evolution in the infected animals during disease progression was analyzed by a bioinformatics approach using ultra-deep pyrosequencing. Our results showed substantial envelope variations emerging in the progressor animal after the onset of AIDS. These envelope variations impacted the length of the variable loops and charges of different envelope regions. Additionally, multiple mutations were located at the CD4 and CCR5 binding sites, potentially affecting receptor binding affinity, viral fitness and they might be selected at late stages of disease. More importantly, these envelope mutations are not random since they had repeatedly been observed in a rhesus macaque and a human infant infected by either SHIV or HIV-1, respectively, carrying the parental envelope of the infectious molecular clone SHIV-1157ipd3N4. Moreover, similar mutations were also observed from other studies on different clades of envelopes regardless of the host species. These recurring mutations in different envelopes suggest that there may be a common evolutionary pattern and selection pathway for the HIV-1 envelope during disease progression