Person:

Goulder, Philip J.

Background: Despite high potential for HIV-1 genetic variation, the emergence of some mutations is constrained by fitness costs, and may be associated with compensatory amino acid (AA) co-variation. To characterize the interplay between Cytotoxic T Lymphocyte (CTL)-mediated pressure and HIV-1 evolutionary pathways, we investigated AA co-variation in Gag sequences obtained from 449 South African individuals chronically infected with HIV-1 subtype C. Methodology/Principal Findings: Individuals with CTL responses biased toward Gag presented lower viral loads than individuals with under-represented Gag-specific CTL responses. Using methods that account for founder effects and HLA linkage disequilibrium, we identified 35 AA sites under Human Leukocyte Antigen (HLA)-restricted CTL selection pressure and 534 AA-to-AA interactions. Analysis of two-dimensional distances between co-varying residues revealed local stabilization mechanisms since 40% of associations involved neighboring residues. Key features of our co-variation analysis included sites with a high number of co-varying partners, such as HLA-associated sites, which had on average 55% more connections than other co-varying sites. Conclusions/Significance: Clusters of co-varying AA around HLA-associated sites (especially at typically conserved sites) suggested that cooperative interactions act to preserve the local structural stability and protein function when CTL escape mutations occur. These results expose HLA-imprinted HIV-1 polymorphisms and their interlinked mutational paths in Gag that are likely due to opposite selective pressures from host CTL-mediated responses and viral fitness constraints.

Retroviruses pack multiple genes into relatively small genomes by encoding several genes in the same genomic region with overlapping reading frames. Both sense and antisense HIV-1 transcripts contain open reading frames for known functional proteins as well as numerous alternative reading frames (ARFs). At least some ARFs have the potential to encode proteins of unknown function, and their antigenic properties can be considered as cryptic epitopes (CEs). To examine the extent of active immune response to virally encoded CEs, we analyzed human leukocyte antigen class I–associated polymorphisms in HIV-1 gag, pol, and nef genes from a large cohort of South Africans with chronic infection. In all, 391 CEs and 168 conventional epitopes were predicted, with the majority (307; 79%) of CEs derived from antisense transcripts. In further evaluation of CD8 T cell responses to a subset of the predicted CEs in patients with primary or chronic infection, both sense- and antisense-encoded CEs were immunogenic at both stages of infection. In addition, CEs often mutated during the first year of infection, which was consistent with immune selection for escape variants. These findings indicate that the HIV-1 genome might encode and deploy a large potential repertoire of unconventional epitopes to enhance vaccine-induced antiviral immunity.

High-resolution HLA typing plays a central role in many areas of immunology, such as in identifying immunogenetic risk factors for disease, in studying how the genomes of pathogens evolve in response to immune selection pressures, and also in vaccine design, where identification of HLA-restricted epitopes may be used to guide the selection of vaccine immunogens. Perhaps one of the most immediate applications is in direct medical decisions concerning the matching of stem cell transplant donors to unrelated recipients. However, high-resolution HLA typing is frequently unavailable due to its high cost or the inability to re-type historical data. In this paper, we introduce and evaluate a method for statistical, in silico refinement of ambiguous and/or low-resolution HLA data. Our method, which requires an independent, high-resolution training data set drawn from the same population as the data to be refined, uses linkage disequilibrium in HLA haplotypes as well as four-digit allele frequency data to probabilistically refine HLA typings. Central to our approach is the use of haplotype inference. We introduce new methodology to this area, improving upon the Expectation-Maximization (EM)-based approaches currently used within the HLA community. Our improvements are achieved by using a parsimonious parameterization for haplotype distributions and by smoothing the maximum likelihood (ML) solution. These improvements make it possible to scale the refinement to a larger number of alleles and loci in a more computationally efficient and stable manner. We also show how to augment our method in order to incorporate ethnicity information (as HLA allele distributions vary widely according to race/ethnicity as well as geographic area), and demonstrate the potential utility of this experimentally. A tool based on our approach is freely available for research purposes at http://microsoft.com/science.

HIV-1 specific HLA-B-restricted CD8+ T cell responses differ from HLA-C-restricted responses in antiviral effectiveness. To investigate possible reasons for these differences, we characterized the frequency and polyfunctionality of immmunodominant HLA-B57/B5801- and HLA-Cw07-restricted CD8+ T cells occurring concurrently in nine study subjects assessing IFN-γ, TNF-α, IL-2, MIP-1β, and CD107a by flow cytometry and analyzed sequence variation in targeted epitopes. HLA-B57/5801 and HLA-Cw07 restricted CD8+ T cells did not differ significantly in polyfunctionality (p = 0.84). Possession of three or more functions correlated positively with CD4+ T cell counts (r = 0.85; p = 0.006) and monofunctional CD8+ T cells inversely correlated with CD4 cell counts (r = −0.79; p = 0.05). There were no differences in polyfunctionality of CD8+ T cells specific to wildtype versus mutated epitopes. These results suggest that loss of polyfunctionality and increase in monofunctional HIV-1-specific CD8+ T cells are associated with disease progression independent of restricting HLA allele. Furthermore, sequence variation does not appear to significantly impact CD8+ T cell polyfunctionality in chronic HIV-1 infection.

Background: A dominance of Gag-specific CD8+ T cell responses is significantly associated with a lower viral load in individuals with chronic, untreated clade C human immunodeficiency virus type 1 (HIV-1) infection. This association has not been investigated in terms of Gag-specific CD4+ T cell responses, nor have clade C HIV-1–specific CD4+ T cell epitopes, likely a vital component of an effective global HIV-1 vaccine, been identified. Methodology/Principal Findings: Intracellular cytokine staining was conducted on 373 subjects with chronic, untreated clade C infection to assess interferon-gamma (IFN-γ) responses by CD4+ T cells to pooled Gag peptides and to determine their association with viral load and CD4 count. Gag-specific IFN-γ–producing CD4+ T cell responses were detected in 261/373 (70%) subjects, with the Gag responders having a significantly lower viral load and higher CD4 count than those with no detectable Gag response (p<0.0001 for both parameters). To identify individual peptides targeted by HIV-1–specific CD4+ T cells, separate ELISPOT screening was conducted on CD8-depleted PBMCs from 32 chronically infected untreated subjects, using pools of overlapping peptides that spanned the entire HIV-1 clade C consensus sequence, and reconfirmed by flow cytometry to be CD4+ mediated. The ELISPOT screening identified 33 CD4+ peptides targeted by 18/32 patients (56%), with 27 of the 33 peptides located in the Gag region. Although the breadth of the CD4+ responses correlated inversely with viral load (p = 0.015), the magnitude of the response was not significantly associated with viral load. Conclusions/Significance: These data indicate that in chronic untreated clade C HIV-1 infection, IFN-γ–secreting Gag-specific CD4+ T cell responses are immunodominant, directed at multiple distinct epitopes, and associated with viral control.

HIV avoids elimination by cytotoxic T-lymphocytes (CTLs) through the evolution of escape mutations. Although there is mounting evidence that these escape pathways are broadly consistent among individuals with similar human leukocyte antigen (HLA) class I alleles, previous population-based studies have been limited by the inability to simultaneously account for HIV codon covariation, linkage disequilibrium among HLA alleles, and the confounding effects of HIV phylogeny when attempting to identify HLA-associated viral evolution. We have developed a statistical model of evolution, called a phylogenetic dependency network, that accounts for these three sources of confounding and identifies the primary sources of selection pressure acting on each HIV codon. Using synthetic data, we demonstrate the utility of this approach for identifying sites of HLA-mediated selection pressure and codon evolution as well as the deleterious effects of failing to account for all three sources of confounding. We then apply our approach to a large, clinically-derived dataset of Gag p17 and p24 sequences from a multicenter cohort of 1144 HIV-infected individuals from British Columbia, Canada (predominantly HIV-1 clade B) and Durban, South Africa (predominantly HIV-1 clade C). The resulting phylogenetic dependency network is dense, containing 149 associations between HLA alleles and HIV codons and 1386 associations among HIV codons. These associations include the complete reconstruction of several recently defined escape and compensatory mutation pathways and agree with emerging data on patterns of epitope targeting. The phylogenetic dependency network adds to the growing body of literature suggesting that sites of escape, order of escape, and compensatory mutations are largely consistent even across different clades, although we also identify several differences between clades. As recent case studies have demonstrated, understanding both the complexity and the consistency of immune escape has important implications for CTL-based vaccine design. Phylogenetic dependency networks represent a major step toward systematically expanding our understanding of CTL escape to diverse populations and whole viral genes.

Background: HLA class-I alleles differ in their ability to control HIV replication through cell-mediated immune responses. No consistent associations have been found between the breadth of Cytotoxic T Lymphocytes (CTL) responses and the control of HIV-1, and it is unknown whether the size or distribution of the viral proteome-wide epitope repertoire, i.e., the intrinsic ability to present fewer, more or specific viral epitopes, could affect clinical markers of disease progression. Methodology/Principal Findings: We used an epitope prediction model to identify all epitope motifs in a set of 302 HIV-1 full-length proteomes according to each individual's HLA (Human Leukocyte Antigen) genotype. The epitope repertoire, i.e., the number of predicted epitopes per HIV-1 proteome, varied considerably between HLA alleles and thus among individual proteomes. In a subgroup of 270 chronically infected individuals, we found that lower viral loads and higher CD4 counts were associated with a larger predicted epitope repertoire. Additionally, in Gag and Rev only, more epitopes were restricted by alleles associated with low viral loads than by alleles associated with higher viral loads. Conclusions/Significance: This comprehensive analysis puts forth the epitope repertoire as a mechanistic component of the multi-faceted HIV-specific CTL response. The favorable impact on markers of disease status of the propensity to present more HLA binding peptides and specific proteins gives impetus to vaccine design strategies that seek to elicit responses to a broad array of HIV-1 epitopes, and suggest a particular focus on Gag.