Modeling HIV quasispecies evolutionary dynamics

Background

During the HIV infection several quasispecies of the virus arise, which are able to use different coreceptors, in particular the CCR5 and CXCR4 coreceptors (R5 and X4 phenotypes, respectively). The switch in coreceptor usage has been correlated with a faster progression of the disease to the AIDS phase. As several pharmaceutical companies are starting large phase III trials for R5 and X4 drugs, models are needed to predict the co-evolutionary and competitive dynamics of virus strains.

Results

We present a model of HIV early infection which describes the dynamics of R5 quasispecies and a model of HIV late infection which describes the R5 to X4 switch. We report the following findings: after superinfection (multiple infections at different times) or coinfection (simultaneous infection by different strains), quasispecies dynamics has time scales of several months and becomes even slower at low number of CD4+ T cells. Phylogenetic inference of chemokine receptors suggests that viral mutational pathway may generate a large variety of R5 variants able to interact with chemokine receptors different from CXCR4. The decrease of CD4+ T cells, during AIDS late stage, can be described taking into account the X4-related Tumor Necrosis Factor dynamics.

Conclusion

The results of this study bridge the gap between the within-patient and the inter-patients (i.e. world-wide) evolutionary processes during HIV infection and may represent a framework relevant for modeling vaccination and therapy.

     

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.     

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₁₀ 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.     

HIV tropism and decreased risk of breast cancer

Background

During the first two decades of the U.S. AIDS epidemic, and unlike some malignancies, breast cancer risk was significantly lower for women with human immunodeficiency virus (HIV) infection compared to the general population. This deficit in HIV-associated breast cancer could not be attributed to differences in survival, immune deficiency, childbearing or other breast cancer risk factors. HIV infects mononuclear immune cells by binding to the CD4 molecule and to CCR5 or CXCR4 chemokine coreceptors. Neoplastic breast cells commonly express CXCR4 but not CCR5. In vitro, binding HIV envelope protein to CXCR4 has been shown to induce apoptosis of neoplastic breast cells. Based on these observations, we hypothesized that breast cancer risk would be lower among women with CXCR4-tropic HIV infection.

Methods and Findings

We conducted a breast cancer nested case-control study among women who participated in the WIHS and HERS HIV cohort studies with longitudinally collected risk factor data and plasma. Cases were HIV-infected women (mean age 46 years) who had stored plasma collected within 24 months of breast cancer diagnosis and an HIV viral load ≥500 copies/mL. Three HIV-infected control women, without breast cancer, were matched to each case based on age and plasma collection date. CXCR4-tropism was determined by a phenotypic tropism assay. Odds ratios (OR) and 95% confidence intervals (CI) for breast cancer were estimated by exact conditional logistic regression. Two (9%) of 23 breast cancer cases had CXCR4-tropic HIV, compared to 19 (28%) of 69 matched controls. Breast cancer risk was significantly and independently reduced with CXCR4 tropism (adjusted odds ratio, 0.10, 95% CI 0.002–0.84) and with menopause (adjusted odds ratio, 0.08, 95% CI 0.001–0.83). Adjustment for CD4⁺ cell count, HIV viral load, and use of antiretroviral therapy did not attenuate the association between infection with CXCR4-tropic HIV and breast cancer.

Conclusions

Low breast cancer risk with HIV is specifically linked to CXCR4-using variants of HIV. These variants are thought to exclusively bind to and signal through a receptor that is commonly expressed on hyperplastic and neoplastic breast duct cells. Additional studies are needed to confirm these observations and to understand how CXCR4 might reduce breast cancer risk.     

Sexually-Transmitted/Founder HIV-1 Cannot Be Directly Predicted from Plasma or PBMC-Derived Viral Quasispecies in the Transmitting Partner

Objective

Characterization of HIV-1 sequences in newly infected individuals is important for elucidating the mechanisms of viral sexual transmission. We report the identification of transmitted/founder viruses in eight pairs of HIV-1 sexually-infected patients enrolled at the time of primary infection (“recipients”) and their transmitting partners (“donors”).

Methods

Using a single genome-amplification approach, we compared quasispecies in donors and recipients on the basis of 316 and 376 C2V5 env sequences amplified from plasma viral RNA and PBMC-associated DNA, respectively.

Results

Both DNA and RNA sequences indicated very homogeneous viral populations in all recipients, suggesting transmission of a single variant, even in cases of recent sexually transmitted infections (STIs) in donors (n = 2) or recipients (n = 3). In all pairs, the transmitted/founder virus was derived from an infrequent variant population within the blood of the donor. The donor variant sequences most closely related to the recipient sequences were found in plasma samples in 3/8 cases and/or in PBMC samples in 6/8 cases. Although donors were exclusively (n = 4) or predominantly (n = 4) infected by CCR5-tropic (R5) strains, two recipients were infected with highly homogeneous CXCR4/dual-mixed-tropic (X4/DM) viral populations, identified in both DNA and RNA. The proportion of X4/DM quasispecies in donors was higher in cases of X4/DM than R5 HIV transmission (16.7–22.0% versus 0–2.6%), suggesting that X4/DM transmission may be associated with a threshold population of X4/DM circulating quasispecies in donors.

Conclusions

These suggest that a severe genetic bottleneck occurs during subtype B HIV-1 heterosexual and homosexual transmission. Sexually-transmitted/founder virus cannot be directly predicted by analysis of the donor’s quasispecies in plasma and/or PBMC. Additional studies are required to fully understand the traits that confer the capacity to transmit and establish infection, and determine the role of concomitant STIs in mitigating the genetic bottleneck in mucosal HIV transmission.     

Differential Role of Autophagy in CD4 T Cells and Macrophages during X4 and R5 HIV-1 Infection

Background

HIV-1 can infect and replicate in both CD4 T cells and macrophages. In these cell types, HIV-1 entry is mediated by the binding of envelope glycoproteins (gp120 and gp41, Env) to the receptor CD4 and a coreceptor, principally CCR5 or CXCR4, depending on the viral strain (R5 or X4, respectively). Uninfected CD4 T cells undergo X4 Env-mediated autophagy, leading to their apoptosis, a mechanism now recognized as central to immunodeficiency.

Methodology/Principal Findings

We demonstrate here that autophagy and cell death are also induced in the uninfected CD4 T cells by HIV-1 R5 Env, while autophagy is inhibited in productively X4 or R5-infected CD4 T cells. In contrast, uninfected macrophages, a preserved cell population during HIV-1 infection, do not undergo X4 or R5 Env-mediated autophagy. Autophagosomes, however, are present in macrophages exposed to infectious HIV-1 particles, independently of coreceptor use. Interestingly, we observed two populations of autophagic cells: one highly autophagic and the other weakly autophagic. Surprisingly, viruses could be detected in the weakly autophagic cells but not in the highly autophagic cells. In addition, we show that the triggering of autophagy in macrophages is necessary for viral replication but addition of Bafilomycin A1, which blocks the final stages of autophagy, strongly increases productive infection.

Conclusions/Significance

Taken together, our data suggest that autophagy plays a complex, but essential, role in HIV pathology by regulating both viral replication and the fate of the target cells.     

Hybrid Approach for Predicting Coreceptor Used by HIV-1 from Its V3 Loop Amino Acid Sequence

Background

HIV-1 infects the host cell by interacting with the primary receptor CD4 and a coreceptor CCR5 or CXCR4. Maraviroc, a CCR5 antagonist binds to CCR5 receptor. Thus, it is important to identify the coreceptor used by the HIV strains dominating in the patient. In past, a number of experimental assays and in-silico techniques have been developed for predicting the coreceptor tropism. The prediction accuracy of these methods is excellent when predicting CCR5(R5) tropic sequences but is relatively poor for CXCR4(X4) tropic sequences. Therefore, any new method for accurate determination of coreceptor usage would be of paramount importance to the successful management of HIV-infected individuals.

Results

The dataset used in this study comprised 1799 R5-tropic and 598 X4-tropic third variable (V3) sequences of HIV-1. We compared the amino acid composition of both types of V3 sequences and observed that certain types of residues, e.g., Asparagine and Isoleucine, were preferred in R5-tropic sequences whereas residues like Lysine, Arginine, and Tryptophan were preferred in X4-tropic sequences. Initially, Support Vector Machine-based models were developed using amino acid composition, dipeptide composition, and split amino acid composition, which achieved accuracy up to 90%. We used BLAST to discriminate R5- and X4-tropic sequences and correctly predicted 93.16% of R5- and 75.75% of X4-tropic sequences. In order to improve the prediction accuracy, a Hybrid model was developed that achieved 91.66% sensitivity, 81.77% specificity, 89.19% accuracy and 0.72 Matthews Correlation Coefficient. The performance of our models was also evaluated on an independent dataset (256 R5- and 81 X4-tropic sequences) and achieved maximum accuracy of 84.87% with Matthews Correlation Coefficient 0.63.

Conclusion

This study describes a highly efficient method for predicting HIV-1 coreceptor usage from V3 sequences. In order to provide a service to the scientific community, a webserver HIVcoPred was developed (http://www.imtech.res.in/raghava/hivcopred/) for predicting the coreceptor usage.     

In-depth characterization of viral isolates from plasma and cells compared with plasma circulating quasispecies in early HIV-1 infection

Background

The use of in vitro models to unravel the phenotypic characteristics of circulating viral variants is key to understanding HIV-1 pathogenesis but limited by the availability of primary viral isolates from biological samples. However, overall in vivo genetic variability of HIV-1 within a subject may not be reflected in the viable viral population obtained after isolation. Although several studies have tried to determine whether viral populations expanded in vitro are representative of in vivo findings, the answer remains unclear due to the reduced number of clonal sequences analyzed or samples compared. In order to overcome previous experimental limitations, here we applied Deep Pyrosequencing (DPS) technology in combination with phenotypic experiments to analyze and compare with unprecedented detail the composition of viral isolates and in vivo quasispecies.

Methodology/Principal Findings

We amplified by DPS HIV-1 genomic regions covering gag, protease, integrase and env-V3 to characterize paired isolates from plasma and peripheral blood mononuclear cells and compare them with total plasma viral RNA in four recently HIV-1 infected subjects. Our study demonstrated the presence of unique haplotypes scattered between sample types with conservation of major variants. In addition, no differences in intra- and inter-population encoded protein variability were found between the different types of isolates or when these were compared to plasma viral RNA within subjects. Additionally, in vitro experiments demonstrated phenotypic similarities in terms of replicative capacity and co-receptor usage between viral isolates and plasma viral RNA.

Conclusion

This study is the first in-depth comparison and characterization of viral isolates from different sources and plasma circulating quasispecies using DPS in recently HIV-1 infected subjects. Our data supports the use of primary isolates regardless of their plasma or cellular origin to define genetic variability and biological traits of circulating HIV-1 quasispecies.     

Virus genome dynamics under different propagation pressures: reconstruction of whole genome haplotypes of west nile viruses from NGS data

Background

Extensive focus is placed on the comparative analyses of consensus genotypes in the study of West Nile virus (WNV) emergence. Few studies account for genetic change in the underlying WNV quasispecies population variants. These variants are not discernable in the consensus genome at the time of emergence, and the maintenance of mutation-selection equilibria of population variants is greatly underestimated. The emergence of lineage 1 WNV strains has been studied extensively, but recent epidemics caused by lineage 2 WNV strains in Hungary, Austria, Greece and Italy emphasizes the increasing importance of this lineage to public health. In this study we explored the quasispecies dynamics of minority variants that contribute to cell-tropism and host determination, i.e. the ability to infect different cell types or cells from different species from Next Generation Sequencing (NGS) data of a historic lineage 2 WNV strain.

Results

Minority variants contributing to host cell membrane association persist in the viral population without contributing to the genetic change in the consensus genome. Minority variants are shown to maintain a stable mutation-selection equilibrium under positive selection, particularly in the capsid gene region.

Conclusions

This study is the first to infer positive selection and the persistence of WNV haplotype variants that contribute to viral fitness without accompanying genetic change in the consensus genotype, documented solely from NGS sequence data. The approach used in this study streamlines the experimental design seeking viral minority variants accurately from NGS data whilst minimizing the influence of associated sequence error.     

X4 Human immunodeficiency virus type 1 gp120 promotes human hepatic stellate cell activation and collagen I expression through interactions with CXCR4

Background & Aims

Patients coinfected with HIV-1 and HCV develop more rapid liver fibrosis than patients monoinfected with HCV. HIV RNA levels correlate with fibrosis progression implicating HIV directly in the fibrotic process. While activated hepatic stellate cells (HSCs) express the 2 major HIV chemokine coreceptors, CXCR4 and CCR5, little is known about the pro-fibrogenic effects of the HIV-1 envelope protein, gp120, on HSCs. We therefore examined the in vitro impact of X4 gp120 on HSC activation, collagen I expression, and underlying signaling pathways and examined the in vivo expression of gp120 in HIV/HCV coinfected livers.

Methods

Primary human HSCs and LX-2 cells, a human HSC line, were challenged with X4 gp120 and expression of fibrogenic markers assessed by qRT-PCR and Western blot +/− either CXCR4-targeted shRNA or anti-CXCR4 neutralizing antibody. Downstream intracellular signaling pathways were evaluated with Western blot and pre-treatment with specific pathway inhibitors. Gp120 immunostaining was performed on HIV/HCV coinfected liver biopsies.

Results

X4 gp 120 significantly increased expression of alpha-smooth muscle actin (a-SMA) and collagen I in HSCs which was blocked by pre-incubation with either CXCR4-targeted shRNA or anti-CXCR4 neutralizing antibody. Furthermore, X4 gp120 promoted Extracellular signal-regulated kinase (ERK) 1/2 phosphorylation and pretreatment with an ERK inhibitor attenuated HSC activation and collagen I expression. Sinusoidal staining for gp120 was evident in HIV/HCV coinfected livers.

Conclusions

X4 HIV-1 gp120 is pro-fibrogenic through its interactions with CXCR4 on activated HSCs. The availability of small molecule inhibitors to CXCR4 make this a potential anti-fibrotic target in HIV/HCV coinfected patients.     

Genetic analysis of HIV-1 subtypes in Nairobi, Kenya

Background

Genetic analysis of a viral infection helps in following its spread in a given population, in tracking the routes of infection and, where applicable, in vaccine design. Additionally, sequence analysis of the viral genome provides information about patterns of genetic divergence that may have occurred during viral evolution.

Objective

In this study we have analyzed the subtypes of Human Immunodeficiency Virus -1 (HIV-1) circulating in a diverse sample population of Nairobi, Kenya.

Methodology

69 blood samples were collected from a diverse subject population attending the Aga Khan University Hospital in Nairobi, Kenya. Total DNA was extracted from peripheral blood mononuclear cells (PBMCs), and used in a Polymerase Chain Reaction (PCR) to amplify the HIV gag gene. The PCR amplimers were partially sequenced, and alignment and phylogenetic analysis of these sequences was performed using the Los Alamos HIV Database.

Results

Blood samples from 69 HIV-1 infected subjects from varying ethnic backgrounds were analyzed. Sequence alignment and phylogenetic analysis showed 39 isolates to be subtype A, 13 subtype D, 7 subtype C, 3 subtype AD and CRF01_AE, 2 subtype G and 1 subtype AC and 1 AG. Deeper phylogenetic analysis revealed HIV subtype A sequences to be highly divergent as compared to subtypes D and C.

Conclusion

Our analysis indicates that HIV-1 subtypes in the Nairobi province of Kenya are dominated by a genetically diverse clade A. Additionally, the prevalence of highly divergent, complex subtypes, intersubtypes, and the recombinant forms indicates viral mixing in Kenyan population, possibly as a result of dual infections.     

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.     

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.     

Monocytes contribute to differential immune pressure on R5 versus X4 HIV through the adipocytokine visfatin/NAMPT

Background

The immune system exerts a diversifying selection pressure on HIV through cellular, humoral and innate mechanisms. This pressure drives viral evolution throughout infection. A better understanding of the natural immune pressure on the virus during infection is warranted, given the clinical interest in eliciting and sustaining an immune response to HIV which can help to control the infection. We undertook to evaluate the potential of the novel HIV-induced, monocyte-derived factor visfatin to modulate viral infection, as part of the innate immune pressure on viral populations.

Results

We show that visfatin is capable of selectively inhibiting infection by R5 HIV strains in macrophages and resting PBMC in vitro, while at the same time remaining indifferent to or even favouring infection by X4 strains. Furthermore, visfatin exerts a direct effect on the relative fitness of R5 versus X4 infections in a viral competition setup. Direct interaction of visfatin with the CCR5 receptor is proposed as a putative mechanism for this differential effect. Possible in vivo relevance of visfatin induction is illustrated by its association with the dominance of CXCR4-using HIV in the plasma.

Conclusions

As an innate factor produced by monocytes, visfatin is capable of inhibiting infections by R5 but not X4 strains, reflecting a potential selective pressure against R5 viruses.     

Analysis of Ultra-Deep Pyrosequencing and Cloning Based Sequencing of the Basic Core Promoter/Precore/Core Region of Hepatitis B Virus Using Newly Developed Bioinformatics Tools

Aims

The aims of this study were to develop bioinformatics tools to explore ultra-deep pyrosequencing (UDPS) data, to test these tools, and to use them to determine the optimum error threshold, and to compare results from UDPS and cloning based sequencing (CBS).

Methods

Four serum samples, infected with either genotype D or E, from HBeAg-positive and HBeAg-negative patients were randomly selected. UDPS and CBS were used to sequence the basic core promoter/precore region of HBV. Two online bioinformatics tools, the “Deep Threshold Tool” and the “Rosetta Tool” (http://hvdr.bioinf.wits.ac.za/tools/), were built to test and analyze the generated data.

Results

A total of 10952 reads were generated by UDPS on the 454 GS Junior platform. In the four samples, substitutions, detected at 0.5% threshold or above, were identified at 39 unique positions, 25 of which were non-synonymous mutations. Sample #2 (HBeAg-negative, genotype D) had substitutions in 26 positions, followed by sample #1 (HBeAg-negative, genotype E) in 12 positions, sample #3 (HBeAg-positive, genotype D) in 7 positions and sample #4 (HBeAg-positive, genotype E) in only four positions. The ratio of nucleotide substitutions between isolates from HBeAg-negative and HBeAg-positive patients was 3.5∶1. Compared to genotype E isolates, genotype D isolates showed greater variation in the X, basic core promoter/precore and core regions. Only 18 of the 39 positions identified by UDPS were detected by CBS, which detected 14 of the 25 non-synonymous mutations detected by UDPS.

Conclusion

UDPS data should be approached with caution. Appropriate curation of read data is required prior to analysis, in order to clean the data and eliminate artefacts. CBS detected fewer than 50% of the substitutions detected by UDPS. Furthermore it is important that the appropriate consensus (reference) sequence is used in order to identify variants correctly.     

Genetic heterogeneity of hepatitis C virus in association with antiviral therapy determined by ultra-deep sequencing

Background and Aims

The hepatitis C virus (HCV) invariably shows wide heterogeneity in infected patients, referred to as a quasispecies population. Massive amounts of genetic information due to the abundance of HCV variants could be an obstacle to evaluate the viral genetic heterogeneity in detail.

Methods

Using a newly developed massive-parallel ultra-deep sequencing technique, we investigated the viral genetic heterogeneity in 27 chronic hepatitis C patients receiving peg-interferon (IFN) α2b plus ribavirin therapy.

Results

Ultra-deep sequencing determined a total of more than 10 million nucleotides of the HCV genome, corresponding to a mean of more than 1000 clones in each specimen, and unveiled extremely high genetic heterogeneity in the genotype 1b HCV population. There was no significant difference in the level of viral complexity between immediate virologic responders and non-responders at baseline (p = 0.39). Immediate virologic responders (n = 8) showed a significant reduction in the genetic complexity spanning all the viral genetic regions at the early phase of IFN administration (p = 0.037). In contrast, non-virologic responders (n = 8) showed no significant changes in the level of viral quasispecies (p = 0.12), indicating that very few viral clones are sensitive to IFN treatment. We also demonstrated that clones resistant to direct-acting antivirals for HCV, such as viral protease and polymerase inhibitors, preexist with various abundances in all 27 treatment-naïve patients, suggesting the risk of the development of drug resistance against these agents.

Conclusion

Use of the ultra-deep sequencing technology revealed massive genetic heterogeneity of HCV, which has important implications regarding the treatment response and outcome of antiviral therapy.     

Phylogenetic Relatedness of Circulating HIV-1C Variants in Mochudi,Botswana

Background

Determining patterns of HIV transmission is increasingly important for the most efficient use of modern prevention interventions. HIV phylogeny can provide a better understanding of the mechanisms underlying HIV transmission networks in communities.

Methods

To reconstruct the structure and dynamics of a local HIV/AIDS epidemic, the phylogenetic relatedness of HIV-1 subtype C env sequences obtained from 785 HIV-infected community residents in the northeastern sector of Mochudi, Botswana, during 2010–2013 was estimated. The genotyping coverage was estimated at 44%. Clusters were defined based on relatedness of HIV-1C env sequences and bootstrap support of splits.

Results

The overall proportion of clustered HIV-1C env sequences was 19.1% (95% CI 17.5% to 20.8%). The proportion of clustered sequences from Mochudi was significantly higher than the proportion of non-Mochudi sequences that clustered, 27.0% vs. 14.7% (p = 5.8E-12; Fisher exact test). The majority of clustered Mochudi sequences (90.1%; 95% CI 85.1% to 93.6%) were found in the Mochudi-unique clusters. None of the sequences from Mochudi clustered with any of the 1,244 non-Botswana HIV-1C sequences. At least 83 distinct HIV-1C variants, or chains of HIV transmission, in Mochudi were enumerated, and their sequence signatures were reconstructed. Seven of 20 genotyped seroconverters were found in 7 distinct clusters.

Conclusions

The study provides essential characteristics of the HIV transmission network in a community in Botswana, suggests the importance of high sampling coverage, and highlights the need for broad HIV genotyping to determine the spread of community-unique and community-mixed viral variants circulating in local epidemics. The proposed methodology of cluster analysis enumerates circulating HIV variants and can work well for surveillance of HIV transmission networks. HIV genotyping at the community level can help to optimize and balance HIV prevention strategies in trials and combined intervention packages.     

Initial Fitness Recovery of HIV-1 Is Associated with Quasispecies Heterogeneity and Can Occur without Modifications in the Consensus Sequence

Background

Fitness recovery of HIV-1 “in vitro” was studied using viral clones that had their fitness decreased as a result of plaque-to-plaque passages.

Principal Findings

After ten large population passages, the viral populations showed an average increase of fitness, although with wide variations among clones. While 5 clones showed significant fitness increases, 3 clones showed increases that were only marginally significant (p<0.1), and 4 clones did not show any change. Fitness recovery was not accompanied by an increase in p24 production, but was associated with an increase in viral titer. Few mutations (an average of 2 mutations per genome) were detected in the consensus nucleotide sequence of the entire genome in all viral populations. Five of the populations did not fix any mutation, and three of them displayed marginally significant fitness increases, illustrating that fitness recovery can occur without detectable alterations of the consensus genomic sequence. The investigation of other possible viral factors associated with the initial steps of fitness recovery, showed that viral quasispecies heterogeneity increased between the initial clones and the passaged populations. A direct statistical correlation between viral heterogeneity and viral fitness was obtained.

Conclusions

Thus, the initial fitness recovery of debilitated HIV-1 clones was mediated by an increase in quasispecies heterogeneity. This observation, together with the invariance of the consensus sequence despite fitness increases demonstrates the relevance of quasispecies heterogeneity in the evolution of HIV-1 in cell culture.     

Multiplex RT-PCR amplification of HIV genes to create a completely autologous DC-based immunotherapy for the treatment of HIV infection

Background

Effective therapy for HIV-infected individuals remains an unmet medical need. Promising clinical trials with dendritic cell (DC)-based immunotherapy consisting of autologous DC loaded with autologous virus have been reported, however, these approaches depend on large numbers of HIV virions to generate sufficient doses for even limited treatment regimens.

Methodology/Principal Findings

The present study describes a novel approach for RT-PCR amplification of HIV antigens. Previously, RT-PCR amplification of autologous viral sequences has been confounded by the high mutation rate of the virus which results in unreliable primer-template binding. To resolve this problem we developed a multiplex RT-PCR strategy that allows reliable strain-independent amplification of highly polymorphic target antigens from any patient and requires neither viral sequence data nor custom-designed PCR primers for each individual. We demonstrate the application of our RT-PCR process to amplify translationally-competent RNA encoding regions of Gag, Vpr, Rev and Nef. The products amplified using this method represent a complex mixture of autologous antigens encoded by viral quasispecies. We further demonstrate that DCs electroporated with in vitro-transcribed HIV RNAs are capable of stimulating poly-antigen-specific CD8+ T cell responses in vitro.

Conclusion/Significance

This study describes a strategy to overcome patient to patient viral diversity enabling strain-independent RT-PCR amplification of RNAs encoding sequence divergent quasispecies of Gag, Vpr, Rev and Nef from small volumes of infectious plasma. The approach allows creation of a completely autologous therapy that does not require advance knowledge of the HIV genomic sequences, does not have yield limitations and has no intact virus in the final product. The simultaneous use of autologous viral antigens and DCs may provoke broad patient-specific immune responses that could potentially induce effective control of viral loads in the absence of conventional antiretroviral drug therapy.