HIV-1 tropism dynamics and phylogenetic analysis from longitudinal ultra-deep sequencing data of CCR5- and CXCR4-using variants

Intra-host dynamics and co-receptor usage of HIV-1 quasi-species in vertically infected patients with phenotypic switch

HIV-1 infection through vertical transmission provides a good model to evaluate intra-host viral evolution and allows to gain insight into the dynamics of viral populations. Our aim was to assess the diversity and dynamics of X4- and R5-using HIV-1 variants in vertically infected children who presented a switch in SI/ NSI phenotype in MT-2 cell assays during chronic infection. Through molecular cloning and next generation sequencing of the C2-V5 env fragment, we investigated HIV-1 evolution and co-receptor usage based on V3 loop prediction bioinformatic tools of longitudinal samples obtained from 4 children. In all cases, the phylogenetic relationships were assessed by Maximum-Likelihood trees constructed with MEGA 6.0. In two cases, V3 loop sequences predicted exclusively R5-using and or X4-using strains, while in another two a higher degree of concordance was observed between the phenotypic and genotypic characteristics. In 3 of the 4 cases, C2-V5 env sequences from different time points were intermingled in phylogenetic trees, with no segregation neither by time or tropism. In only one case monophyletic clustering defined groups of sequences with different co-receptor usage. Comparison of amino acid frequency between isolates with SI and NSI phenotype allowed the identification of 9 possible genetic determinants in subtype F C2-V5 region of env associated to SI/ NSI phenotype in these patients, one of which had previously been reported for subtype B. Overall, we found a low degree of correlation between phenotypic and genotypic properties of HIV-1 quasispecies in patients under chronic infection. Whether HIV-1 subtype or other factors influence the evolution of HIV-1 in vivo will require further research.

Dolutegravir reshapes the genetic diversity of HIV-1 reservoirs

Objectives

Better understanding of the dynamics of HIV reservoirs under ART is a critical step to achieve a functional HIV cure. Our objective was to assess the genetic diversity of archived HIV-1 DNA over 48 weeks in blood cells of individuals starting treatment with a dolutegravir-based regimen.

Methods

Eighty blood samples were prospectively and longitudinally collected from 20 individuals (NCT02557997) including: acutely (n = 5) and chronically (n = 5) infected treatment-naive individuals, as well as treatment-experienced individuals who switched to a dolutegravir-based regimen and were either virologically suppressed (n = 5) or had experienced treatment failure (n = 5). The integrase and V3 loop regions of HIV-1 DNA isolated from PBMCs were analysed by pyrosequencing at baseline and weeks 4, 24 and 48. HIV-1 genetic diversity was calculated using Shannon entropy.

Results

All individuals achieved or maintained viral suppression throughout the study. A low and stable genetic diversity of archived HIV quasispecies was observed in individuals starting treatment during acute infection. A dramatic reduction of the genetic diversity was observed at week 4 of treatment in the other individuals. In these patients and despite virological suppression, a recovery of the genetic diversity of the reservoirs was observed up to 48 weeks. Viral variants bearing dolutegravir resistance-associated substitutions at integrase position 50, 124, 230 or 263 were detected in five individuals (n = 5/20, 25%) from all groups except those who were ART-failing at baseline. None of these substitutions led to virological failure.

Conclusions

These data demonstrate that the genetic diversity of the HIV-1 reservoir is reshaped following the initiation of a dolutegravir-based regimen and strongly suggest that HIV-1 can continue to replicate despite successful treatment.

Application of deep sequencing methods for inferring viral population diversity

The first deep sequencing method was announced in 2005. Due to an increasing number of sequencing data and a reduction in the costs of each sequencing dataset, this innovative technique was soon applied to genetic investigations of viral genome diversity in various viruses, particularly RNA viruses. These deep sequencing findings documented viral epidemiology and evolution and provided high-resolution data on the genetic changes in viral populations. Here, we review deep sequencing platforms that have been applied in viral quasispecies studies. Further, we discuss recent deep sequencing studies on viral inter- and intrahost evolution, drug resistance, and humoral immune selection, especially in emerging and re-emerging viruses. Deep sequencing methods are becoming the standard for providing comprehensive results of viral population diversity, and their applications are discussed.

Clonal analysis of HIV-1 genotype and function associated with virologic failure in treatment-experienced persons receiving maraviroc: Results from the MOTIVATE phase 3 randomized, placebo-controlled trials

Detailed clonal phenotypic/genotypic analyses explored viral-escape mechanisms during maraviroc-based therapy in highly treatment-experienced participants from the MOTIVATE trials. To allow real-time assessment of samples while maintaining a blind trial, the first 267 enrolled participants were selected for evaluation. At failure, plasma samples from 20/50 participants (16/20 maraviroc-treated) with CXCR4-using virus and all 38 (13 maraviroc-treated) with CCR5-tropic virus were evaluated. Of those maraviroc-treated participants with CXCR4-using virus at failure, genotypic and phenotypic clonal tropism determinations showed >90% correspondence in 14/16 at Day 1 and 14/16 at failure. Phylogenetic analysis of clonal sequences detected pre-treatment progenitor CXCR4-using virus, or on-treatment virus highly divergent from the Day 1 R5 virus, excluding possible co-receptor switch through maraviroc-mediated evolution. Re-analysis of pre-treatment samples using the enhanced-sensitivity Trofile® assay detected CXCR4-using virus pre-treatment in 16/20 participants failing with CXCR4-using virus. Post-maraviroc reversion of CXCR4-use to CCR5-tropic occurred in 7/8 participants with follow-up, suggesting selective maraviroc inhibition of CCR5-tropic variants in a mixed-tropic viral population, not emergence of de novo mutations in CCR5-tropic virus, as the main virologic escape mechanism. Maraviroc-resistant CCR5-tropic virus was observed in plasma from 5 treated participants with virus displaying reduced maximal percent inhibition (MPI) but no evidence of IC50 change. Env clones with reduced MPI showed 1-5 amino acid changes specific to each V3-loop region of env relative to Day 1. However, transferring on-treatment resistance-associated changes using site-directed mutagenesis did not always establish resistance in Day 1 virus, and key 'signature' mutation patterns associated with reduced susceptibility to maraviroc were not identified. Evolutionary divergence of the CXCR4-using viruses is confirmed, emphasizing natural selection not influenced directly by maraviroc; maraviroc simply unmasks pre-existing lineages by inhibiting the R5 virus. For R5-viral failure, resistance development through drug selection pressure was uncommon and manifested through reduced MPI and with virus strain-specific mutational patterns.

Full-Length Envelope Analyzer (FLEA): A tool for longitudinal analysis of viral amplicons

Next generation sequencing of viral populations has advanced our understanding of viral population dynamics, the development of drug resistance, and escape from host immune responses. Many applications require complete gene sequences, which can be impossible to reconstruct from short reads. HIV env, the protein of interest for HIV vaccine studies, is exceptionally challenging for long-read sequencing and analysis due to its length, high substitution rate, and extensive indel variation. While long-read sequencing is attractive in this setting, the analysis of such data is not well handled by existing methods. To address this, we introduce FLEA (Full-Length Envelope Analyzer), which performs end-to-end analysis and visualization of long-read sequencing data. FLEA consists of both a pipeline (optionally run on a high-performance cluster), and a client-side web application that provides interactive results. The pipeline transforms FASTQ reads into high-quality consensus sequences (HQCSs) and uses them to build a codon-aware multiple sequence alignment. The resulting alignment is then used to infer phylogenies, selection pressure, and evolutionary dynamics. The web application provides publication-quality plots and interactive visualizations, including an annotated viral alignment browser, time series plots of evolutionary dynamics, visualizations of gene-wide selective pressures (such as dN/dS) across time and across protein structure, and a phylogenetic tree browser. We demonstrate how FLEA may be used to process Pacific Biosciences HIV env data and describe recent examples of its use. Simulations show how FLEA dramatically reduces the error rate of this sequencing platform, providing an accurate portrait of complex and variable HIV env populations. A public instance of FLEA is hosted at http://flea.datamonkey.org. The Python source code for the FLEA pipeline can be found at https://github.com/veg/flea-pipeline. The client-side application is available at https://github.com/veg/flea-web-app. A live demo of the P018 results can be found at http://flea.murrell.group/view/P018.

Viral populations constantly evolve and diversify. In this article we introduce a method, FLEA, for reconstructing and visualizing the details of evolutionary changes. FLEA specifically processes data from sequencing platforms that generate reads that are long, but error-prone. To study the evolutionary dynamics of entire genes during viral infection, data is collected via long-read sequencing at discrete time points, allowing us to understand how the virus changes over time. However, the experimental and sequencing process is imperfect, so the resulting data contain not only real evolutionary changes, but also mutations and other genetic artifacts caused by sequencing errors. Our method corrects most of these errors by combining thousands of erroneous sequences into a much smaller number of unique consensus sequences that represent biologically meaningful variation. The resulting high-quality sequences are used for further analysis, such as building an evolutionary tree that tracks and interprets the genetic changes in the viral population over time. FLEA is open source, and is freely available online.

Ultradeep single-molecule real-time sequencing of HIV envelope reveals complete compartmentalization of highly macrophage-tropic R5 proviral variants in brain and CXCR4-using variants in immune and peripheral tissues

Despite combined antiretroviral therapy (cART), HIV+ patients still develop neurological disorders, which may be due to persistent HIV infection and selective evolution in brain tissues. Single-molecule real-time (SMRT) sequencing technology offers an improved opportunity to study the relationship among HIV isolates in the brain and lymphoid tissues because it is capable of generating thousands of long sequence reads in a single run. Here, we used SMRT sequencing to generate ~ 50,000 high-quality full-length HIV envelope sequences (> 2200 bp) from seven autopsy tissues from an HIV+/cART+ subject, including three brain and four non-brain sites. Sanger sequencing was used for comparison with SMRT data and to clone functional pseudoviruses for in vitro tropism assays. Phylogenetic analysis demonstrated that brain-derived HIV was compartmentalized from HIV outside the brain and that the variants from each of the three brain tissues grouped independently. Variants from all peripheral tissues were intermixed on the tree but independent of the brain clades. Due to the large number of sequences, a clustering analysis at three similarity thresholds (99, 99.5, and 99.9%) was also performed. All brain sequences clustered exclusive of any non-brain sequences at all thresholds; however, frontal lobe sequences clustered independently of occipital and parietal lobes. Translated sequences revealed potentially functional differences between brain and non-brain sequences in the location of putative N-linked glycosylation sites (N-sites), V1 length, V3 charge, and the number of V4 N-sites. All brain sequences were predicted to use the CCR5 co-receptor, while most non-brain sequences were predicted to use CXCR4 co-receptor. Tropism results were confirmed by in vitro infection assays. The study is the first to use a SMRT sequencing approach to study HIV compartmentalization in tissues and supports other reports of limited trafficking between brain and non-brain sequences during cART. Due to the long sequence length, we could observe changes along the entire envelope gene, likely caused by differential selective pressure in the brain that may contribute to neurological disease.

High-Sequence Diversity and Rapid Virus Turnover Contribute to Higher Rates of Coreceptor Switching in Treatment-Experienced Subjects with HIV-1 Viremia

Coreceptor switching from CCR5 to CXCR4 is common during chronic HIV-1 infection, but is even more common in individuals who have failed antiretroviral therapy (ART). Prior studies have suggested rapid mutation and/or recombination of HIV-1 envelope (env) genes during coreceptor switching. We compared the functional and genotypic changes in env of viruses from viremic subjects who had failed ART just before and after coreceptor switching and compared those to viruses from matched subjects without coreceptor switching. Analysis of multiple unique functional env clones from each subject revealed extensive diversity at both sample time points and rapid diversification of sequences during the 4-month interval in viruses from both 9 subjects with coreceptor switching and 15 control subjects. Only two subjects had envs with evidence of recombination. Three findings distinguished env clones from subjects with coreceptor switching from controls: (1) lower entry efficiency via CCR5; (2) longer V1/V2 regions; and (3), lower nadir CD4 T cell counts during prior years of infection. Most of these subjects harbored virus with lower replicative capacity associated with protease (PR) and/or reverse transcriptase inhibitor resistance mutations, and the extensive diversification tended to lead either to improved entry efficiency via CCR5 or the gain of entry function via CXCR4. These results suggest that R5X4 or X4 variants emerge from a diverse, low-fitness landscape shaped by chronic infection, multiple ART resistance mutations, the availability of target cells, and reduced entry efficiency via CCR5.

Possible involvement of distinct phylogenetic clusters of HIV-1 variants in the discrepancies between coreceptor tropism predictions based on viral RNA and proviral DNA

Background

The coreceptor tropism testing should be conducted prior to commencing a regimen containing a CCR5 antagonist for treatment of HIV-1 infection. For aviremic patients on long antiretroviral therapy, proviral DNA is often used instead of viral RNA in genotypic tropism testing. However, the tropism predictions from RNA and DNA are sometimes different. We examined the cause of the discrepancies between HIV-1 tropism predictions based on viral RNA and proviral DNA.

Methods

The nucleotide sequence of the env C2V3C3 region was determined using pair samples of plasma RNA and peripheral blood mononuclear cell DNA from 50 HIV-1 subtype B-infected individuals using population-based sequencing. The samples with discrepant tropism assessments between RNA and DNA were further analyzed using deep sequencing, followed by phylogenetic analysis. The tropism was assessed using the algorithm geno2pheno with a false-positive rate cutoff of 10 %.

Results

In population-based sequencing, five of 50 subjects showed discrepant tropism predictions between their RNA and DNA samples: four exhibited R5 tropism in RNA and X4 tropism in DNA, while one exhibited the opposite pattern. In the deep sequencing and phylogenetic analysis, three subjects had single clusters comprising of RNA- and DNA-derived sequences that were a mixture of R5 and X4 sequences. The other two subjects had two and three distinct phylogenetic clusters of sequences, respectively, each of which was dominated by R5 or X4 sequences; sequences of the R5-dominated cluster were mostly found in RNA, while sequences of the X4-dominated cluster were mostly in DNA.

Conclusions

Some of HIV-1 tropism discrepancies between viral RNA and proviral DNA seem to be caused by phylogenetically distinct clusters which resides in plasma and cells in different frequencies. Our findings suggest that the tropism testing using PBMC DNA or deep sequencing may be required when the viral load is not suppressed or rebounds in the course of a CCR5 antagonist-containing regimen.

Deep Sequencing of the HIV-1 env Gene Reveals Discrete X4 Lineages and Linkage Disequilibrium between X4 and R5 Viruses in the V1/V2 and V3 Variable Regions

HIV-1 requires the CD4 receptor and a coreceptor (CCR5 [R5 phenotype] or CXCR4 [X4 phenotype]) to enter cells. Coreceptor tropism can be assessed by either phenotypic or genotypic analysis, the latter using bioinformatics algorithms to predict tropism based on the env V3 sequence. We used the Primer ID sequencing strategy with the MiSeq sequencing platform to reveal the structure of viral populations in the V1/V2 and C2/V3 regions of the HIV-1 env gene in 30 late-stage and 6 early-stage subjects. We also used endpoint dilution PCR followed by cloning of env genes to create pseudotyped virus to explore the link between genotypic predictions and phenotypic assessment of coreceptor usage. We found out that the most stringently sequence-based calls of X4 variants (Geno2Pheno false-positive rate [FPR] of ≤2%) formed distinct lineages within the viral population, and these were detected in 24 of 30 late-stage samples (80%), which was significantly higher than what has been seen previously by using other approaches. Non-X4 lineages were not skewed toward lower FPR scores in X4-containing populations. Phenotypic assays showed that variants with an intermediate FPR (2 to 20%) could be either X4/dual-tropic or R5 variants, although the X4 variants made up only about 25% of the lineages with an FPR of <10%, and these variants carried a distinctive sequence change. Phylogenetic analysis of both the V1/V2 and C2/V3 regions showed evidence of recombination within but very little recombination between the X4 and R5 lineages, suggesting that these populations are genetically isolated.

IMPORTANCE

Primer ID sequencing provides a novel approach to study genetic structures of viral populations. X4 variants may be more prevalent than previously reported when assessed by using next-generation sequencing (NGS) and with a greater depth of sampling than single-genome amplification (SGA). Phylogenetic analysis to identify lineages of sequences with intermediate FPR values may provide additional information for accurately predicting X4 variants by using V3 sequences. Limited recombination occurs between X4 and R5 lineages, suggesting that X4 and R5 variants are genetically isolated and may be replicating in different cell types or that X4/R5 recombinants have reduced fitness.

Evolution of hepatitis C virus in HIV coinfected patients under antiretroviral therapy

Five patients (P) were followed-up for an average of 7.73years after highly active antiretroviral therapy (HAART) initiation. Patients' immune and virological status were determined by periodical CD4+T-cell counts and HIV and HCV viral load. HCV populations were studied using longitudinal high throughput sequence data obtained in parallel by virological and immunological parameters. Two patients (P7, P28) with sub-optimal responses to HAART presented HCV viral loads significantly higher than those recorded for two patients (P1, P18) that achieved good responses to HAART. Interestingly, HCV populations from P7 and P28 displayed a stable phylogenetic structure, whereas HCV populations from P1 and P18showeda significant increase in their phylogenetic structure, followed by a decrease after achieving acceptable CD4+T-cell counts (>500 cell/μl). The fifth patient (P25) presented high HCV viral loads, preserved CD4+T-cell counts from baseline and all along the follow-up, and displayed a constant viral phylogenetic structure. These results strongly suggest that HAART-induced immune recovery induces a decrease in HCV viral load and an increase in the HCV population phylogenetic structure likely reflecting the virus diversification in response to the afresh immune response. The relatively low HCV viral load observed in the HAART responder patients suggests that once HCV is adapted it reaches a maximum number of haplotypes higher than that achieved during the initial stages of the immune response as inferred from the two recovering patients. Future studies using larger number of patients are needed to corroborate these hypotheses.

Evolution of coreceptor utilization to escape CCR5 antagonist therapy

The HIV-1 envelope interacts with coreceptors CCR5 and CXCR4 in a dynamic, multi-step process, its molecular details not clearly delineated. Use of CCR5 antagonists results in tropism shift and therapeutic failure. Here we describe a novel approach using full-length patient-derived gp160 quasispecies libraries cloned into HIV-1 molecular clones, their separation based on phenotypic tropism in vitro, and deep sequencing of the resultant variants for structure-function analyses. Analysis of functionally validated envelope sequences from patients who failed CCR5 antagonist therapy revealed determinants strongly associated with coreceptor specificity, especially at the gp120-gp41 and gp41-gp41 interaction surfaces that invite future research on the roles of subunit interaction and envelope trimer stability in coreceptor usage. This study identifies important structure-function relationships in HIV-1 envelope, and demonstrates proof of concept for a new integrated analysis method that facilitates laboratory discovery of resistant mutants to aid in development of other therapeutic agents.

Drug resistance and tropism as markers of the dynamics of HIV-1 DNA quasispecies in blood cells of heavily pretreated patients who achieved sustained virological suppression

OBJECTIVES

The objective of this study was to address the dynamics of archived resistant quasispecies in cell-associated HIV-1 DNA over time in heavily ART-experienced patients with currently suppressed plasma HIV-1 RNA.

METHODS

Longitudinal ultra-deep sequencing (UDS) analysis of reverse transcriptase, protease and V3 Env regions was performed on blood-cell-associated HIV-1 DNA samples. Drug-resistance-associated mutations (DRAMs) and tropism were interpreted using the ANRS and Geno2Pheno algorithms. We analysed frozen blood cells from patients enrolled in the INNOVE and ANRS 123 ETOILE studies who achieved sustained viral suppression after salvage optimized ART (SOT).

RESULTS

Samples were available at baseline and 6 and ≥12 months after SOT initiation in 10 patients. V3 loop sequences displayed wide intra-individual dynamics over time. Viral variants harbouring DRAMs exhibited three non-exclusive scenarios. First, when SOT exerted the same selective pressure as previous failing regimens, some viral quasispecies still harboured the same DRAMs at the same level as at the time of virological failure. Thus, as DRAMs were mostly associated with the same viral variant, variants with a complete resistance pattern remained archived. Second, some viral variants harbouring DRAMs were no longer detected over time when SOT consisted of new antiretroviral classes or had resistance profiles distinct from those of previous failing regimens. Third, variants with new DRAMs associated with SOT emerged in blood cells during follow-up despite sustained virological control.

CONCLUSIONS

Using longitudinal UDS analysis and focusing on DRAMs and tropism as markers, we demonstrated that, despite sustained virological control, archived HIV-1 DNA quasispecies continued to evolve.

Differential evolution of a CXCR4-using HIV-1 strain in CCR5wt/wt and CCR5∆32/∆32 hosts revealed by longitudinal deep sequencing and phylogenetic reconstruction

Rare individuals homozygous for a naturally-occurring 32 base pair deletion in the CCR5 gene (CCR5∆32/∆32) are resistant to infection by CCR5-using ("R5") HIV-1 strains but remain susceptible to less common CXCR4-using ("X4") strains. The evolutionary dynamics of X4 infections however, remain incompletely understood. We identified two individuals, one CCR5wt/wt and one CCR5∆32/∆32, within the Vancouver Injection Drug Users Study who were infected with a genetically similar X4 HIV-1 strain. While early-stage plasma viral loads were comparable in the two individuals (~4.5-5 log10 HIV-1 RNA copies/ml), CD4 counts in the CCR5wt/wt individual reached a nadir of <20 CD4 cells/mm(3) within 17 months but remained >250 cells/mm(3) in the CCR5∆32/∆32 individual. Ancestral phylogenetic reconstructions using longitudinal envelope-V3 deep sequences suggested that both individuals were infected by a single transmitted/founder (T/F) X4 virus that differed at only one V3 site (codon 24). While substantial within-host HIV-1 V3 diversification was observed in plasma and PBMC in both individuals, the CCR5wt/wt individual's HIV-1 population gradually reverted from 100% X4 to ~60% R5 over ~4 years whereas the CCR5∆32/∆32 individual's remained consistently X4. Our observations illuminate early dynamics of X4 HIV-1 infections and underscore the influence of CCR5 genotype on HIV-1 V3 evolution.