Evolution of human immunodeficiency virus type 1 coreceptor usage during antiretroviral Therapy: a Bayesian approach

Long-term antiretroviral therapy mitigates mortality and morbidity independent of HIV tropism: 18 years follow-up in a women's cohort

OBJECTIVE

CXCR4 (X4)-tropic HIV-1 was found previously to herald CD4 + cell depletion and disease progression in individuals who were antiretroviral-naive or took combination antiretroviral therapy (cART) for less than 5 years. We updated this finding by investigating whether the deleterious effect of X4-tropic strains is mitigated by long-term cART.

DESIGN

We examined morbidity and mortality in relation to HIV-1 tropism and cART in 529 participants followed up to 18 years in the Women's Interagency HIV Study; 91% were women of color.

METHODS

Plasma-derived HIV-1 tropism was determined genotypically.

RESULTS

We categorized participants according to the number of visits reported on cART after initiation. Group 1: three or less visits, 74% of these participants reporting no cART; group 2: at least four visits and less than 70% of visits on cART; group 3: at least 70% of visits on cART. AIDS mortality rates for participants in each group with X4 virus compared with those with R5 virus exclusively were, respectively: 62 vs. 40% ( P  = 0.0088); 23% vs. 22% [nonsignificant (NS)]; 7% vs. 14% (NS). Kaplan-Meier curves showed accelerated progression to AIDS death or AIDS-defining illness in participants with three or less cART visits and X4 viruses ( P  = 0.0028) but no difference in progression rates stratified by tropism in other groups. Logistic regression found that HIV-1 suppression for at least 10 semiannual visits (≥5 years total) mitigated X4 tropism's deleterious effect on mortality, controlling for maximal viral load, and CD4 + nadir.

CONCLUSION

Long-term cART markedly mitigated the deleterious effect of X4 viruses on AIDS morbidity and mortality. Mitigation was correlated with duration of viral suppression, supporting HIV-1 suppression as a crucial goal.

Contribution of the HIV-1 Envelope Glycoprotein to AIDS Pathogenesis and Clinical Progression

In the absence of antiviral therapy, HIV-1 infection progresses to a wide spectrum of clinical manifestations that are the result of an entangled contribution of host, immune and viral factors. The contribution of these factors is not completely established. Several investigations have described the involvement of the immune system in the viral control. In addition, distinct HLA-B alleles, HLA-B27, -B57-58, were associated with infection control. The combination of these elements and antiviral host restriction factors results in different clinical outcomes. The role of the viral proteins in HIV-1 infection has been, however, less investigated. We will review contributions dedicated to the pathogenesis of HIV-1 infection focusing on studies identifying the function of the viral envelope glycoprotein (Env) in the clinical progression because of its essential role in the initial events of the virus life-cycle. Some analysis showed that inefficient viral Envs were dominant in non-progressor individuals. These poorly-functional viral proteins resulted in lower cellular activation, viral replication and minor viral loads. This limited viral antigenic production allows a better immune response and a lower immune exhaustion. Thus, the properties of HIV-1 Env are significant in the clinical outcome of the HIV-1 infection and AIDS pathogenesis.

The Characteristics of the HIV-1 Env Glycoprotein Are Linked With Viral Pathogenesis

The understanding of HIV-1 pathogenesis and clinical progression is incomplete due to the variable contribution of host, immune, and viral factors. The involvement of viral factors has been investigated in extreme clinical phenotypes from rapid progressors to long-term non-progressors (LTNPs). Among HIV-1 proteins, the envelope glycoprotein complex (Env) has been concentrated on in many studies for its important role in the immune response and in the first steps of viral replication. In this study, we analyzed the contribution of 41 Envs from 24 patients with different clinical progression rates and viral loads (VLs), LTNP-Elite Controllers (LTNP-ECs); Viremic LTNPs (vLTNPs), and non-controller individuals contemporary to LTNPs or recent, named Old and Modern progressors. We studied the Env expression, the fusion and cell-to-cell transfer capacities, as well as viral infectivity. The sequence and phylogenetic analysis of Envs were also performed. In every functional characteristic, the Envs from subjects with viral control (LTNP-ECs and vLTNPs) showed significant lower performance compared to those from the progressor individuals (Old and Modern). Regarding sequence analysis, the variable loops of the gp120 subunit of the Env (i.e., V2, V4, and mainly V5) of the progressor individuals showed longer and more glycosylated sequences than controller subjects. Therefore, HIV-1 Envs from virus of patients presenting viremic control and the non-progressor clinical phenotype showed poor viral functions and shorter sequences, whereas functional Envs were associated with virus of patients lacking virological control and with progressor clinical phenotypes. These correlations support the role of Env genotypic and phenotypic characteristics in the in vivo HIV-1 infection and pathogenesis.

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.

Cladograms with Path to Event (ClaPTE): a novel algorithm to detect associations between genotypes or phenotypes using phylogenies

BACKGROUND

Associations between genotype and phenotype provide insight into the evolution of pathogenesis, drug resistance, and the spread of pathogens between hosts. However, common ancestry can lead to apparent associations between biologically unrelated features. The novel method Cladograms with Path to Event (ClaPTE) detects associations between character-pairs (either a pair of mutations or a mutation paired with a phenotype) while adjusting for common ancestry, using phylogenetic trees.

METHODS

ClaPTE tests for character-pairs changing close together on the phylogenetic tree, consistent with an associated character-pair. ClaPTE is compared to three existing methods (independent contrasts, mixed model, and likelihood ratio) to detect character-pair associations adjusted for common ancestry. Comparisons utilize simulations on gene trees for: HIV Env, HIV promoter, and bacterial DnaJ and GuaB; and case studies for Oseltamavir resistance in Influenza, and for DnaJ and GuaB. Simulated data include both true-positive/associated character-pairs, and true-negative/not-associated character-pairs, used to assess type I (frequency of p-values in true-negatives) and type II (sensitivity to true-positives) error control.

RESULTS AND CONCLUSIONS

ClaPTE has competitive sensitivity and better type I error control than existing methods. In the Influenza/Oseltamavir case study, ClaPTE reports no new permissive mutations but detects associations between adjacent (in primary sequence) amino acid positions which other methods miss. In the DnaJ and GuaB case study, ClaPTE reports more frequent associations between positions both from the same protein family than between positions from different families, in contrast to other methods. In both case studies, the results from ClaPTE are biologically plausible.

Natural Scrub Typhus Antibody Suppresses HIV CXCR4(X4) Viruses

Viral load generally rises in HIV-infected individuals with a concomitant infection, but falls markedly in some individuals with scrub typhus (ST), a common Asian rickettsial infection. ST infection appears to shift the viral population from CXCR4-using (X4) to CCR5-utilizing (R5) strains, and there is evidence of cross-reactivity between ST-specific antibodies and HIV-1. We examined the mechanism of ST suppression of HIV by measuring the effects of ST infection on X4 and R5 viruses in vivo and in vitro, and assessing the relative contributions of antibodies and chemokines to the inhibitory effect. In vivo, a single scrub typhus plasma infusion markedly reduced the subpopulation of HIV-1 viruses using the X4 co-receptor in all 8 recipients, and eliminated X4 viruses 6 patients. In vitro, the 14 ST sera tested all inhibited the replication of an X4 but not an R5 virus. This inhibitory effect was maintained if ST sera were depleted of chemokines but was lost upon removal of antibodies. Sera from STinfected mice recognized a target that co-localized with X4 HIV gp120 in immunofluorescent experiments. These in vivo and in vitro data suggest that acute ST infection generates cross-reactive antibodies that produce potent suppression of CXCR4- but not CCR5-using HIV-1 viruses. ST suppression of HIV replication could reveal novel mechanisms that could be exploited for vaccination strategies, as well as aid in the development of fusion inhibitors and other new therapeutic regimens. This also appears to be the first instance where one pathogen is neutralized by antibody produced in response to infection by a completely unrelated organism.

Emergence, dominance, and possible decline of CXCR4 chemokine receptor usage during the course of HIV infection

Binding to a chemokine receptor, either CCR5 or CXCR4, by the gp120 glycoprotein is an essential step in the pathway by which HIV enters host cells. Recently, CCR5 antagonists have been developed that obstruct binding of CCR5 by gp120, thus inhibiting host cell entry. Resistance to such CCR5 antagonists may emerge, however, through the selection of viral strains capable of utilizing CXCR4 receptors. This study explores the evolutionary context of emergence, and in many cases decline, of dominant CXCR4-usage (X4) during disease progression within a number of individuals. Of seven individuals exhibiting a switch to dominant CXCR4 usage, such dominance is transient in five of them with CCR5-usage (R5) re-emerging to dominate the viral population later in disease progression. Three individuals conform to documented X4 transience in that the re-emergence of R5 dominance is an outgrowth from the predominant R5 strain. However, in two individuals we observe a novel pathway for R5 re-emergence in that R5 strains emerge to dominate late in disease progression through continued evolution of the X4 population. This suggests that the molecular mechanism of such switches between R5 and X4-usage is strain specific and that no single mechanism is shared between individuals. These findings have implications for the understanding of the mechanisms of potential emergence of resistance to CCR5 antagonists through use of the CXCR4 receptor and support the importance to have an appropriately optimized background therapy for use with entry inhibitors and, as for all HAART, to monitor drug resistance in a comprehensive manner.

Impact of CCR5delta32 host genetic background and disease progression on HIV-1 intrahost evolutionary processes: efficient hypothesis testing through hierarchical phylogenetic models

The interplay between C-C chemokine receptor type 5 (CCR5) host genetic background, disease progression, and intrahost HIV-1 evolutionary dynamics remains unclear because differences in viral evolution between hosts limit the ability to draw conclusions across hosts stratified into clinically relevant populations. Similar inference problems are proliferating across many measurably evolving pathogens for which intrahost sequence samples are readily available. To this end, we propose novel hierarchical phylogenetic models (HPMs) that incorporate fixed effects to test for differences in dynamics across host populations in a formal statistical framework employing stochastic search variable selection and model averaging. To clarify the role of CCR5 host genetic background and disease progression on viral evolutionary patterns, we obtain gp120 envelope sequences from clonal HIV-1 variants isolated at multiple time points in the course of infection from populations of HIV-1-infected individuals who only harbored CCR5-using HIV-1 variants at all time points. Presence or absence of a CCR5 wt/Δ32 genotype and progressive or long-term nonprogressive course of infection stratify the clinical populations in a two-way design. As compared with the standard approach of analyzing sequences from each patient independently, the HPM provides more efficient estimation of evolutionary parameters such as nucleotide substitution rates and d(N)/d(S) rate ratios, as shown by significant shrinkage of the estimator variance. The fixed effects also correct for nonindependence of data between populations and results in even further shrinkage of individual patient estimates. Model selection suggests an association between nucleotide substitution rate and disease progression, but a role for CCR5 genotype remains elusive. Given the absence of clear d(N)/d(S) differences between patient groups, delayed onset of AIDS symptoms appears to be solely associated with lower viral replication rates rather than with differences in selection on amino acid fixation.

Evolution and recombination of genes encoding HIV-1 drug resistance and tropism during antiretroviral therapy

Characterization of residual plasma virus during antiretroviral therapy (ART) is a high priority to improve understanding of HIV-1 pathogenesis and therapy. To understand the evolution of HIV-1 pol and env genes in viremic patients under selective pressure of ART, we performed longitudinal analyses of plasma-derived pol and env sequences from single HIV-1 genomes. We tested the hypotheses that drug resistance in pol was unrelated to changes in coreceptor usage (tropism), and that recombination played a role in evolution of viral strains. Recombinants were identified by using Bayesian and other computational methods. High-level genotypic resistance was seen in approximately 70% of X4 and R5 strains during ART. There was no significant association between resistance and tropism. Each patient displayed at least one recombinant encompassing env and representing a change in predicted tropism. These data suggest that, in addition to mutation, recombination can play a significant role in shaping HIV-1 evolution.

Modulation of the virus-receptor interaction by mutations in the V5 loop of feline immunodeficiency virus (FIV) following in vivo escape from neutralising antibody

Background

In the acute phase of infection with feline immunodeficiency virus (FIV), the virus targets activated CD4+ T cells by utilising CD134 (OX40) as a primary attachment receptor and CXCR4 as a co-receptor. The nature of the virus-receptor interaction varies between isolates; strains such as GL8 and CPGammer recognise a "complex" determinant on CD134 formed by cysteine-rich domains (CRDs) 1 and 2 of the molecule while strains such as PPR and B2542 require a more "simple" determinant comprising CRD1 only for infection. These differences in receptor recognition manifest as variations in sensitivity to receptor antagonists. In this study, we ask whether the nature of the virus-receptor interaction evolves in vivo.

Results

Following infection with a homogeneous viral population derived from a pathogenic molecular clone, a quasispecies emerged comprising variants with distinct sensitivities to neutralising antibody and displaying evidence of conversion from a "complex" to a "simple" interaction with CD134. Escape from neutralising antibody was mediated primarily by length and sequence polymorphisms in the V5 region of Env, and these alterations in V5 modulated the virus-receptor interaction as indicated by altered sensitivities to antagonism by both anti-CD134 antibody and soluble CD134.

Conclusions

The FIV-receptor interaction evolves under the selective pressure of the host humoral immune response, and the V5 loop contributes to the virus-receptor interaction. Our data are consistent with a model whereby viruses with distinct biological properties are present in early versus late infection and with a shift from a "complex" to a "simple" interaction with CD134 with time post-infection.

CD4 T-cell hyperactivation and susceptibility to cell death determine poor CD4 T-cell recovery during suppressive HAART

BACKGROUND

The failure to increase CD4 T-cell counts in some HAART-treated HIV-infected patients with satisfactory virological responses has been related to low CD4 T-cell production, high turnover and death. However, the relative contribution of these factors is still unclear, strongly limiting the definition of appropriate therapeutic strategies for these patients.

METHODS

A cross-sectional study was designed to evaluate the contribution of thymic activity, microbial translocation, cellular activation and death to CD4 T-cell recovery. We included 230 HIV-infected individuals on suppressive HAART (>2 years); 95 of them were considered 'discordant' (CD4 T-cell count <350 cells/mul) and 135 were considered 'concordant'. Comparative and logistic regression analyses were performed.

RESULTS

Discordant patients showed higher levels of activated [human leukocyte antigen (HLA)-DRCD95 and CD38CD45RA] cells in both the CD8 and CD4 T-cell compartments. Notably, the most significant differences were observed in CD4 T cells. Discordant patients showed lower naive CD4 T-cell production (CD45RACD31 cells), higher spontaneous ex-vivo CD4 T-cell death and higher plasma levels of soluble CD14. Multivariate analysis showed that activation and death of CD4 T cells, along with nadir CD4 T-cell counts, were the only predictive factors for poor immune recovery. Moreover, the low correlations found between CD4 T-cell activation or death with thymic output and bacterial translocation suggest that additional factors modulate cellular activation and death and, in turn, CD4 T-cell recovery.

CONCLUSION

CD4 T-cell repopulation during HAART is determined by CD4 T-cell activation and death. Therefore, strategies aimed to reduce these parameters should be envisaged to treat discordant patients.

Two-way Bayesian hierarchical phylogenetic models: An application to the co-evolution of gp120 and gp41 during and after enfuvirtide treatment

Enfuvirtide (ENF) is a fusion inhibitor that prevents the entry of HIV virions into target cells. Studying the characteristics of viral evolution during treatment and after a treatment interruption can lend insight into the mechanisms of viral evolution and fitness. Although interruption of anti-HIV therapy often results in rapid emergence of an archived "wild-type" virus population, previous work from our group indicates that when only ENF is interrupted, viral gp41 continues to evolve forward and resistance mutations are lost due to back-mutation and remodeling of the envelope protein. To examine the co-evolution of gp120 and gp41 during ENF interruption we extend the Bayesian Hierarchical Phylogenetic model (HPM). Current HPMs enforce conditional independence across all outcomes while biologically all gene regions within a patient should return the same tree unless recombination confers an evolutionary selective advantage. A two-way-interaction HPM is proposed that provides middle ground between these two extremes and allows us to test for differences in evolutionary pressures across gene regions in multiple patients simultaneously. When the model is applied to a well-characterized cohort of HIV-infected patients interrupting ENF we find that across patients, the virus continued to evolve forward in both gene regions. Overall, the hypothesis of independence over dependence between the gene regions is supported. Models that allow for the examination of co-evolution over time will be increasingly important as more therapeutic classes are developed, each of which may impact other through novel and complex mechanisms.

Reappearance of an 11-year-old sequence in an HIV-1 infected patient during treatment interruption

HIV-1 from a patient with multi-drug resistant virus was identified as wild type during treatment interruption. The aim of the study was to describe how the viral population is affected by treatment interruptions and use phylogeny to reconstruct the evolutionary pattern. 15 samples covering 13 y and 2 treatment interruptions were analysed in both pol and env. The wild type virus found in the sample from the second treatment interruption in 2002 had not been present as a dominant population since 1994. Phylogeny showed that the 2002 sample was more closely related to wild type sequences than to other sequences sampled in 2002. This indicated that the wild type virus was caused by recruitment from the viral archives rather than reversion of previously circulating resistant strains. A few weeks after re-initiated treatment, virus showed full resistance, indicating that resistant virus was present as a subpopulation and reselected due to higher fitness in the presence of drugs. Phylogeny of env showed that CCR5 and CXCR4 viruses coexist in the patient. In conclusion, the study showed that at all times during infection, virus is archived in the cells and can be recruited when the surrounding environment changes and the archived virus is more fit.

Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics

Kingman's coalescent process opens the door for estimation of population genetics model parameters from molecular sequences. One paramount parameter of interest is the effective population size. Temporal variation of this quantity characterizes the demographic history of a population. Because researchers are rarely able to choose a priori a deterministic model describing effective population size dynamics for data at hand, nonparametric curve-fitting methods based on multiple change-point (MCP) models have been developed. We propose an alternative to change-point modeling that exploits Gaussian Markov random fields to achieve temporal smoothing of the effective population size in a Bayesian framework. The main advantage of our approach is that, in contrast to MCP models, the explicit temporal smoothing does not require strong prior decisions. To approximate the posterior distribution of the population dynamics, we use efficient, fast mixing Markov chain Monte Carlo algorithms designed for highly structured Gaussian models. In a simulation study, we demonstrate that the proposed temporal smoothing method, named Bayesian skyride, successfully recovers "true" population size trajectories in all simulation scenarios and competes well with the MCP approaches without evoking strong prior assumptions. We apply our Bayesian skyride method to 2 real data sets. We analyze sequences of hepatitis C virus contemporaneously sampled in Egypt, reproducing all key known aspects of the viral population dynamics. Next, we estimate the demographic histories of human influenza A hemagglutinin sequences, serially sampled throughout 3 flu seasons.

Continued evolution in gp41 after interruption of enfuvirtide in subjects with advanced HIV type 1 disease

Careful examination of viral dynamics during antiretroviral treatment can provide important insights into HIV pathogenesis. We examined viral evolution during and after enfuvirtide (T-20) treatment with an objective of defining the characteristics of viral evolution during advanced HIV immunodeficiency. Specifically, we examined the viral quasispecies from nine patients who experienced incomplete viral suppression on an enfuvirtide-based regimen, and who subsequently interrupted enfuvirtide while remaining on a stable optimized background regimen (enfuvirtide "partial treatment interruption"). On average, eight clones were sequenced from three time points: pre-enfuvirtide, post-enfuvirtide failure, and post-enfuvirtide interruption. We utilized a Bayesian hierarchical phylogenetic model to assess the evolutionary relatedness of the virus that emerged over time. In most subjects, interruption of enfuvirtide was associated with continued viral evolution as enfuvirtide mutations waned; in no subject did we find clear evidence supporting the reemergence of archived wild-type variants (Bayes factor=30.2, p=0.01). Evidence supporting ongoing viral evolution was particularly strong in subjects whose virus remained diverse or became more diverse during enfuvirtide therapy. In contrast to observations when all drugs are interrupted, loss of resistance during enfuvirtide interruption is most likely due to ongoing viral evolution (and back-mutation), rather than emergence of an archived virus.

Genetic and stochastic influences on the interaction of human immunodeficiency virus type 1 and cytotoxic T lymphocytes in identical twins

Human immunodeficiency virus type 1 (HIV-1) evolves in vivo under selective pressure from CD8+ T-lymphocyte (CTL) responses, which are in turn determined by host and viral genetic factors, such as restricting major histocompatibility complex molecules and the available viral epitope sequences. However, CTL are derived stochastically through the random gene rearrangements to produce T-cell receptors (TCR), and the relative impact of genetic versus stochastic processes on CTL targeting of HIV and immune-driven viral evolution is unclear. Here we evaluate identical twins infected with HIV-1 as neonates from a common blood transfusion, with subsequently similar environmental exposures, thereby allowing controlled comparisons of CTL targeting and viral evolution. Seventeen years after infection, their CTL targeting of HIV-1 was remarkably similar. In contrast, their overall TCR profiles were highly dissimilar, and a dominant epitope was recognized by distinctly different TCR in each twin. Furthermore, their viral epitopes had diverged, and there was ongoing viral phylogenetic divergence between the twins between 12 and 17 years after infection. These results indicate that while CTL targeting is predominately genetically determined, stochastic influences render the interaction of HIV-1 and host immunity, and therefore viral escape and CTL efficacy, unpredictable.