Intersubtype differences in the effect of a rare p24 gag mutation on HIV-1 replicative fitness

Changing selection on amino acid substitutions in Gag protein between major HIV-1 subtypes

Amino acid preferences at a protein site depend on the role of this site in protein function and structure as well as on external constraints. All these factors can change in the course of evolution, making amino acid propensities of a site time-dependent. When viral subtypes divergently evolve in different host subpopulations, such changes may depend on genetic, medical, and sociocultural differences between these subpopulations. Here, using our previously developed phylogenetic approach, we describe sixty-nine amino acid sites of the Gag protein of human immunodeficiency virus type 1 (HIV-1) where amino acids have different impact on viral fitness in six major subtypes of the type M. These changes in preferences trigger adaptive evolution; indeed, 32 (46 per cent) of these sites experienced strong positive selection at least in one of the subtypes. At some of the sites, changes in amino acid preferences may be associated with differences in immune escape between subtypes. The prevalence of an amino acid in a protein site within a subtype is only a poor predictor for whether this amino acid is preferred in this subtype according to the phylogenetic analysis. Therefore, attempts to identify the factors of viral evolution from comparative genomics data should integrate across multiple sources of information.

Genotypic and Mechanistic Characterization of Subtype-Specific HIV Adaptation to Host Cellular Immunity

The extent to which viral genetic context influences HIV adaptation to human leukocyte antigen (HLA) class I-restricted immune pressures remains incompletely understood. The Ugandan HIV epidemic, where major pandemic group M subtypes A1 and D cocirculate in a single host population, provides an opportunity to investigate this question. We characterized plasma HIV RNA gag, pol, and nef sequences, along with host HLA genotypes, in 464 antiretroviral-naive individuals chronically infected with HIV subtype A1 or D. Using phylogenetically informed statistical approaches, we identified HLA-associated polymorphisms and formally compared their strengths of selection between viral subtypes. A substantial number (32%) of HLA-associated polymorphisms identified in subtype A1 and/or D had previously been reported in subtype B, C, and/or circulating recombinant form 01_AE (CRF01_AE), confirming the shared nature of many HLA-driven escape pathways regardless of viral genetic context. Nevertheless, 34% of the identified HLA-associated polymorphisms were significantly differentially selected between subtypes A1 and D. Experimental investigation of select examples of subtype-specific escape revealed distinct underlying mechanisms with important implications for vaccine design: whereas some were attributable to subtype-specific sequence variation that influenced epitope-HLA binding, others were attributable to differential mutational barriers to immune escape. Overall, our results confirm that HIV genetic context is a key modulator of viral adaptation to host cellular immunity and highlight the power of combined bioinformatic and mechanistic studies, paired with knowledge of epitope immunogenicity, to identify appropriate viral regions for inclusion in subtype-specific and universal HIV vaccine strategies.IMPORTANCE The identification of HIV polymorphisms reproducibly selected under pressure by specific HLA alleles and the elucidation of their impact on viral function can help identify immunogenic viral regions where immune escape incurs a fitness cost. However, our knowledge of HLA-driven escape pathways and their functional costs is largely limited to HIV subtype B and, to a lesser extent, subtype C. Our study represents the first characterization of HLA-driven adaptation pathways in HIV subtypes A1 and D, which dominate in East Africa, and the first statistically rigorous characterization of differential HLA-driven escape across viral subtypes. The results support a considerable impact of viral genetic context on HIV adaptation to host HLA, where HIV subtype-specific sequence variation influences both epitope-HLA binding and the fitness costs of escape. Integrated bioinformatic and mechanistic characterization of these and other instances of differential escape could aid rational cytotoxic T-lymphocyte-based vaccine immunogen selection for both subtype-specific and universal HIV vaccines.

PTAP motif duplication in the p6 Gag protein confers a replication advantage on HIV-1 subtype C

HIV-1 subtype C (HIV-1C) may duplicate longer amino acid stretches in the p6 Gag protein, leading to the creation of an additional Pro-Thr/Ser-Ala-Pro (PTAP) motif necessary for viral packaging. However, the biological significance of a duplication of the PTAP motif for HIV-1 replication and pathogenesis has not been experimentally validated. In a longitudinal study of two different clinical cohorts of select HIV-1 seropositive, drug-naive individuals from India, we found that 8 of 50 of these individuals harbored a mixed infection of viral strains discordant for the PTAP duplication. Conventional and next-generation sequencing of six primary viral quasispecies at multiple time points disclosed that in a mixed infection, the viral strains containing the PTAP duplication dominated the infection. The dominance of the double-PTAP viral strains over a genetically similar single-PTAP viral clone was confirmed in viral proliferation and pairwise competition assays. Of note, in the proximity ligation assay, double-PTAP Gag proteins exhibited a significantly enhanced interaction with the host protein tumor susceptibility gene 101 (Tsg101). Moreover, Tsg101 overexpression resulted in a biphasic effect on HIV-1C proliferation, an enhanced effect at low concentration and an inhibitory effect only at higher concentrations, unlike a uniformly inhibitory effect on subtype B strains. In summary, our results indicate that the duplication of the PTAP motif in the p6 Gag protein enhances the replication fitness of HIV-1C by engaging the Tsg101 host protein with a higher affinity. Our results have implications for HIV-1 pathogenesis, especially of HIV-1C.

Subtype-Specific Differences in Gag-Protease-Driven Replication Capacity Are Consistent with Intersubtype Differences in HIV-1 Disease Progression

There are marked differences in the spread and prevalence of HIV-1 subtypes worldwide, and differences in clinical progression have been reported. However, the biological reasons underlying these differences are unknown. Gag-protease is essential for HIV-1 replication, and Gag-protease-driven replication capacity has previously been correlated with disease progression. We show that Gag-protease replication capacity correlates significantly with that of whole isolates (r = 0.51; P = 0.04), indicating that Gag-protease is a significant contributor to viral replication capacity. Furthermore, we investigated subtype-specific differences in Gag-protease-driven replication capacity using large well-characterized cohorts in Africa and the Americas. Patient-derived Gag-protease sequences were inserted into an HIV-1 NL4-3 backbone, and the replication capacities of the resulting recombinant viruses were measured in an HIV-1-inducible reporter T cell line by flow cytometry. Recombinant viruses expressing subtype C Gag-proteases exhibited substantially lower replication capacities than those expressing subtype B Gag-proteases (P < 0.0001); this observation remained consistent when representative Gag-protease sequences were engineered into an HIV-1 subtype C backbone. We identified Gag residues 483 and 484, located within the Alix-binding motif involved in virus budding, as major contributors to subtype-specific replicative differences. In East African cohorts, we observed a hierarchy of Gag-protease-driven replication capacities, i.e., subtypes A/C < D < intersubtype recombinants (P < 0.0029), which is consistent with reported intersubtype differences in disease progression. We thus hypothesize that the lower Gag-protease-driven replication capacity of subtypes A and C slows disease progression in individuals infected with these subtypes, which in turn leads to greater opportunity for transmission and thus increased prevalence of these subtypes.IMPORTANCE HIV-1 subtypes are unevenly distributed globally, and there are reported differences in their rates of disease progression and epidemic spread. The biological determinants underlying these differences have not been fully elucidated. Here, we show that HIV-1 Gag-protease-driven replication capacity correlates with the replication capacity of whole virus isolates. We further show that subtype B displays a significantly higher Gag-protease-mediated replication capacity than does subtype C, and we identify a major genetic determinant of these differences. Moreover, in two independent East African cohorts we demonstrate a reproducible hierarchy of Gag-protease-driven replicative capacity, whereby recombinants exhibit the greatest replication, followed by subtype D, followed by subtypes A and C. Our data identify Gag-protease as a major determinant of subtype differences in disease progression among HIV-1 subtypes; furthermore, we propose that the poorer viral replicative capacity of subtypes A and C may paradoxically contribute to their more efficient spread in sub-Saharan Africa.

The PTAP sequence duplication in HIV-1 subtype C Gag p6 in drug-naive subjects of India and South Africa

Background

HIV-1 subtype C demonstrates several biological properties distinct from other viral subtypes. One such variation is the duplication of PTAP motif in p6 Gag. PTAP motif is a key player in viral budding. Here, we studied the prevalence of PTAP motif duplication in subtype C viral strains in a longitudinal study.

Methods

In a prospective follow-up study, 65 HIV-1 seropositive drug-naive subjects were monitored in two different clinical cohorts of India for 2 years with repeated sampling at 6-month intervals. The viral RNA was extracted from plasma, the gag segment was amplified and sequenced. From a subset of viral isolates the sequences of pol, env and LTR were sequenced. Using HIV-1 gag amino acid sequences available from public databases and additional sequences derived from the Indian and South-African cohorts, we examined the nature of PTAP motif duplication in subtype C.

Results

In 16% (8 of 50) of the primary viral strains of India, we identified a sequence duplication of the PTAP motif in Gag p6. The length of the sequence duplication varied from 6 to 14 amino acids in the viral isolates but remained fixed within a subject over a period of 24–36 month follow-up. In the duplicated motif, the core PTAP motif was invariable, but the flanking residues were highly variable. In an acute phase clinical cohort of South Africa, in a subset of 75 subjects, we found the presence of the PTAP duplication at a frequency of 29.3%. An analysis of the gag sequences from the extant databases showed that unlike other subtypes of HIV-1, subtype C has a natural propensity to generate the PTAP motif duplication at a significantly higher frequency and of greater length. Additionally, the global prevalence of PTAP duplication in subtype C appears to be increasing progressively over the past 30 years.

Conclusion

We showed that in subtype C, the duplication of the PTAP motif in p6 Gag involves sequence stretches of greater length, and at a much higher frequency as compared to other HIV-1 subtypes. Given that subtype C naturally lacks the Alix binding motif, the acquisition of an additional PTAP motif may confer replication advantage on this HIV-1 subtype. Further investigation is warranted to examine the significance of PTAP motif duplication on the replicative fitness of HIV-1.

Reversion and T cell escape mutations compensate the fitness loss of a CD8+ T cell escape mutant in their cognate transmitted/founder virus

Immune escape mutations that revert back to the consensus sequence frequently occur in newly HIV-1-infected individuals and have been thought to render the viruses more fit. However, their impact on viral fitness and their interaction with other immune escape mutations have not been evaluated in the background of their cognate transmitted/founder (T/F) viral genomes. To precisely determine the role of reversion mutations, we introduced reversion mutations alone or together with CD8⁺ T cell escape mutations in their unmodified cognate T/F viral genome and determined their impact on viral fitness in primary CD4⁺ T cells. Two reversion mutations, V247I and I64T, were identified in Gag and Tat, respectively, but neither had measurable effect on the fitness of their cognate T/F virus. The V247I and G248A mutations that were detected before and concurrently with the potent T cell escape mutation T242N, respectively, were selected by early T cell responses. The V247I or the G248A mutation alone partially restored the fitness loss caused by the T242N mutation. Together they could fully restore the fitness of the T242N mutant to the T/F level. These results demonstrate that the fitness loss caused by a T cell escape mutation could be compensated by preexisting or concurrent reversion and other T cell escape mutations. Our findings indicate that the overall viral fitness is modulated by the complex interplay among T cell escape, compensatory and reversion mutations to maintain the balance between immune escape and viral replication capacity.

Impact of HLA selection pressure on HIV fitness at a population level in Mexico and Barbados

Previous studies have demonstrated that effective cytotoxic T lymphocyte (CTL) responses drive the selection of escape mutations that reduce viral replication capacity (VRC). Escape mutations, including those with reduced VRC, can be transmitted and accumulate in a population. Here we compared two antiretroviral therapy (ART)-naive HIV clade B-infected cohorts, in Mexico and Barbados, in which the most protective HLA alleles (HLA-B*27/57/58:01/81:01) are differentially expressed, at 8% and 34%, respectively. Viral loads were significantly higher in Mexico than in Barbados (median, 40,774 versus 14,200; P < 0.0001), and absolute CD4⁺ T-cell counts were somewhat lower (median, 380/mm³ versus 403/mm³; P = 0.007). We tested the hypothesis that the disparate frequencies of these protective HLA alleles would be associated with a higher VRC at the population level in Mexico. Analysis of VRC in subjects in each cohort, matched for CD4⁺ T-cell count, revealed that the VRC was indeed higher in the Mexican cohort (mean, 1.13 versus 1.03; P = 0.0025). Although CD4 counts were matched, viral loads remained significantly higher in the Mexican subjects (P = 0.04). This VRC difference was reflected by a significantly higher frequency in the Barbados cohort of HLA-B*27/57/58:01/81:01-associated Gag escape mutations previously shown to incur a fitness cost on the virus (P = 0.004), a difference between the two cohorts that remained statistically significant even in subjects not expressing these protective alleles (P = 0.01). These data suggest that viral set points and disease progression rates at the population level may be significantly influenced by the prevalence of protective HLA alleles such as HLA-B*27/57/58:01/81:01 and that CD4 count-based guidelines to initiate antiretroviral therapy may need to be modified accordingly, to optimize the effectiveness of treatment-for-prevention strategies and reduce HIV transmission rates to the absolute minimum.

IMPORTANCE Immune control of HIV at an individual level is strongly influenced by the HLA class I genotype. HLA class I molecules mediating effective immune control, such as HLA-B*27 and HLA-B*57, are associated with the selection of escape mutants that reduce viral replicative capacity. The escape mutants selected in infected patients can be transmitted and affect the viral load and CD4 count in the recipient. These findings prompt the hypothesis that the frequency of protective alleles in a population may affect viral set points and rates of disease progression in that population. These studies in Mexico and Barbados, where the prevalence rates of protective HLA alleles are 8% and 34%, respectively, support this hypothesis. These data suggest that antiretroviral therapy (ART) treatment-for-prevention strategies will be less successful in populations such as those in Mexico, where viral loads are higher for a given CD4 count. Consideration may therefore usefully be given to ART initiation at higher absolute CD4 counts in such populations to optimize the impact of ART for prevention.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

HIV-1 conserved-element vaccines: relationship between sequence conservation and replicative capacity

To overcome the problem of HIV-1 variability, candidate vaccine antigens have been designed to be composed of conserved elements of the HIV-1 proteome. Such candidate vaccines could be improved with a better understanding of both HIV-1 evolutionary constraints and the fitness cost of specific mutations. We evaluated the in vitro fitness cost of 23 mutations engineered in the HIV-1 subtype B Gag-p24 Center-of-Tree (COT) protein through fitness competition assays. While some mutations at conserved sites exacted a high fitness cost, as expected under the assumption that the most conserved residue confers the highest fitness, there was no overall strong relationship between sequence conservation and replicative capacity. By comparing sites that have evolved since the beginning of the epidemic to those that have remain unchanged, we found that sites that have evolved over time were more likely to correspond to HLA-associated sites and that their mutation had limited fitness costs. Our data showed no transcendent link between high conservation and high fitness cost, indicating that merely focusing on conserved segments of HIV-1 would not be sufficient for a successful vaccine strategy. Nonetheless, a subset of sites exacted a high fitness cost upon mutation--these sites have been under selective pressure to change since the beginning of the epidemic but have proved virtually nonmutable and could constitute preferred targets for vaccine design.