Abstract
Immunodominance describes a phenomenon whereby the immune system consistently targets only a fraction of the available Ag pool derived from a given pathogen. In the case of CD8+ T cells, these constrained epitope-targeting patterns are linked to HLA class I expression and determine disease progression. Despite the biological importance of these predetermined response hierarchies, little is known about the factors that control immunodominance in vivo. In this study, we conducted an extensive analysis of CD8+ T cell responses restricted by a single HLA class I molecule to evaluate the mechanisms that contribute to epitope-targeting frequency and antiviral efficacy in HIV-1 infection. A clear immunodominance hierarchy was observed across 20 epitopes restricted by HLA-B*42:01, which is highly prevalent in populations of African origin. Moreover, in line with previous studies, Gag-specific responses and targeting breadth were associated with lower viral load set-points. However, peptide–HLA-B*42:01 binding affinity and stability were not significantly linked with targeting frequencies. Instead, immunodominance correlated with epitope-specific usage of public TCRs, defined as amino acid residue–identical TRB sequences that occur in multiple individuals. Collectively, these results provide important insights into a potential link between shared TCR recruitment, immunodominance, and antiviral efficacy in a major human infection.
Introduction
CD8+ T cells play a central role in the immune response to viruses (1). However, clear-cut differences exist between distinct specificities in terms of antiviral efficacy (2, 3). Moreover, epitope-targeting patterns are often predetermined within a hierarchy of immunodominance across restriction elements (4–12). These biologically imposed limitations can therefore dictate the outcome of certain viral infections.
In the case of HIV-1, disease progression is strongly affected by the expression of particular HLA alleles (13, 14). Although several mechanisms have been proposed to explain this association (13), one key factor is the relative ability of epitope-specific CD8+ T cell populations to kill HIV-infected cells. For example, Gag-specific responses typically contain HIV-1 more effectively in vivo (15–19) and in vitro (20, 21) compared with those that target other viral proteins. In addition, protective HLA class I molecules, such as HLA-B*27 and HLA-B*57 (22–24), restrict immunodominant Gag-specific responses that select viral escape variants with impaired replicative capacity (22, 25–27).
A number of factors can influence immunodominance (2, 4, 6, 11). Ag presentation itself is the culmination of several upstream events, such as the kinetics of protein expression, the abundance of protein delivered into the cytoplasm, intracytoplasmic proteosomal cleavage, translocation by TAP, peptide loading onto MHC, ERAP1/2 trimming, and transport to the cell surface (28). Peptide–MHC binding affinity and stability determine the subsequent availability of Ag over time, whereas TCR avidity and the frequency of naive precursors govern the size of the potentially responsive T cell pool (7, 9, 10, 12). The phenomenon of immunodomination, whereby certain responses are subordinated in the presence of particular high-frequency responses (29), also plays a role. Nonetheless, despite the complexity of these multifaceted processes, the end result is a largely predictable pattern of immunodominance for any given virus.
Emerging studies highlight a key role for the TCR repertoire as an independent determinant of antiviral efficacy in multiple systems (30–35). Although the process of V(D)J rearrangement can theoretically generate 1015–20 distinct TCRs (36), extreme biases exist during recombination, thymic selection, naive T cell recruitment, and subsequent clonal expansion (6, 10, 36–38). These biases can ultimately generate identical or “public” epitope-specific TCRs in multiple individuals (38). However, the role of TCR bias in relation to immunodominance remains ill-defined, with murine studies yielding apparently contradictory data (39, 40).
To assess the potential influence of TCR publicity on CD8+ T cell immunodominance and antiviral efficacy in a human viral infection, we conducted an extensive analysis of HLA-B*42:01–restricted responses directed against an array of different epitopes derived from HIV-1. The prevalence of HLA-B*42:01 in most populations of African origin, combined with the substantial repertoire of associated viral peptides, enabled a large-scale study controlled for the restriction element across multiple individuals with chronic HIV-1 infection. Consequently, we were able to detect statistically meaningful correlations between these parameters.
Materials and Methods
Study subjects
The study cohort comprised 2093 female adults with chronic, antiretroviral therapy–naive C-clade HIV-1 infection, recruited from five cohorts: Durban, South Africa (14, 17, 25, 41); Gaborone, Botswana (42); Bloemfontein, South Africa (43); Kimberley, South Africa (44); and Thames Valley, U.K. (45). A total of 246 HLA-B*42:01+ individuals with documented proviral DNA sequences, CD4+ T cell counts, plasma viral loads, and four-digit HLA-genotyping data was identified within the entire cohort, from which 181 were screened with overlapping peptides (OLPs) to map HIV-specific CD8+ T cell responses in IFN-γ ELISPOT assays. Informed consent was obtained from all participants. The following Institutional Review Boards approved the study: University of KwaZulu-Natal; University of the Free State, South Africa; Health Research Development Committee, Botswana Ministry of Health, Botswana; Office of Human Research Administration, Harvard School of Public Health; and University of Oxford.
IFN-γ ELISPOT
Virus-specific CD8+ T cell responses across the whole HIV-1 proteome were determined for 1009 C-clade–infected subjects via direct ex vivo IFN-γ ELISPOT analysis. Ags comprised 410 OLPs based on the C-clade consensus (2001) arranged in a matrix system with 11–12 peptides/pool. Responses to matrix pools were deconvoluted by subsequent testing with the individual 18-mer peptides contained in each pool (15). Associations between HLA-B*42:01 expression and OLP targeting were calculated from a total of 181 HLA-B*42:01+ individuals. Percent targeting frequency (immunodominance rank) calculations for HLA-B*42:01–restricted epitopes required the exclusion of 27 individuals coexpressing HLA-B*07:02/39:10/42:02/81:01 to eliminate alternative presentation by other B7 superfamily members (46). Viral load set-point calculations versus targeting of HLA-B*42:01–restricted OLPs were based exclusively on data from the Durban cohort (n = 126) to minimize the influence of external factors.
HLA class I typing
Four-digit HLA-A/B/C genotyping was performed using real-time reverse sequence-specific oligonucleotide kits (Dynal), as described previously (47).
Epitope mapping and HLA restriction
Epitope mapping and recognition assays were performed as described previously (16, 45). Single allele–matched B lymphoblastoid cell lines (BLCLs) were used to determine HLA restriction (16, 48).
Tetramers
Tetrameric peptide–HLA complexes were generated and used as described previously (16, 49). Samples were acquired using an LSR II flow cytometer (BD Biosciences) and analyzed with FlowJo version 8.8.6 (TreeStar). Events were gated serially on singlets, lymphocytes, live cells, and CD3+ T cells prior to analysis in bivariate tetramer versus CD8 plots.
Proviral DNA sequencing
Sequences from Gag (n = 1857), Pol (n = 1052), and Nef (n = 1327) were generated by extraction of genomic DNA from PBMCs and amplification via nested PCRs, as described previously (17, 44, 47). Purified products were sequenced using the Big Dye Terminator Ready Reaction (Life Technologies) (41, 50). Vif, Rev, and Env sequences were available from 255 subjects (51).
Peptide–HLA binding and stability assays
Peptide–HLA stability was measured as described previously (52). Peptide–HLA binding affinity was measured via AlphaScreen technology (PerkinElmer) (53).
TCR clonotyping
Clonotypic analysis of Ag-specific CD8+ T cell populations restricted by HLA-B*42:01 was performed as described previously (54). Briefly, viable HIV-specific tetramer+ CD8+ T cell populations (n = 48) were sorted by flow cytometry at >98% purity directly into RNAlater. A median of 2120 cells (25th percentile = 1068 cells, 75th percentile = 4340 cells) was sorted per population, with a median response magnitude of 0.60% (25th percentile = 0.3%, 75th percentile = 1.4%). The number of sorted cells correlated with response magnitude (r = +0.8, p < 0.0001). Unbiased amplification of all expressed TRB gene products was conducted using a template-switch anchored RT-PCR with a 3′ C region primer. Amplicons were subcloned, sampled, sequenced (n = 3592), and analyzed as described previously (55). IMGT nomenclature is used in this report (56).
Statistical analysis
Associations between HLA-B*42:01 expression, HIV-1 polymorphisms, and 18-mer peptide (OLP) responses were determined as described previously (44, 57). The Dunn multiple-comparisons test was used to compare median viral loads between OLP responders and nonresponders, the number of OLP responses between different HIV-1 proteins (breadth), and the number of different TCR clonotypes between responses. The two-tailed Mann–Whitney U test was used to compare viral loads and CD4+ T cell counts between individuals carrying wild-type sequences or escape polymorphisms within the Gag-RM9 epitope. All p values were calculated using GraphPad Prism version 6.0c (GraphPad). Differences in Gag-RM9 sequence polymorphisms between HLA-B*42:01+ and HLA-B*42:01− subjects were calculated using the Fisher exact test. Correlations between percent targeting frequency (IFN-γ OLP ELISPOT responders) and either peptide–HLA binding affinities (IC50 KD value) or peptide–HLA binding half-lives (hours) were calculated using the Spearman rank test. Clonotypic data were normalized as described previously (58–60). Briefly, all samples were rarefied down to the lowest estimated coverage (88%) prior to calculations of TCR sharing. This process was repeated 10,000 times, and mean publicity scores were used for statistical analysis.
Results
HLA-B*42:01–restricted responses conform to a strict immunodominance hierarchy
Initially, we used IFN-γ ELISPOT assays to screen 1009 antiretroviral therapy–naive C-clade–infected female individuals, 181 (18%) of whom carried the HLA-B*42:01 allele, for responses to a panel of 410 OLPs (15, 16, 61) spanning the entire HIV-1 proteome (Table I). This panel comprised seven HLA-B*42:01–restricted epitopes already listed in the Los Alamos Immunology Database (http://www.hiv.lanl.gov) (62). In addition, we identified 13 novel epitopes via statistical associations between recognition of a particular 18-mer OLP and the expression of HLA-B*42:01 (q < 0.05). This set of 20 HLA-B*42:01–restricted HIV-1 epitopes was used for further analysis.
Next, we ranked the response hierarchy among these 20 epitopes based on targeting frequency (immunodominance rank). A total of 154 HLA-B*42:01+ individuals was included in this analysis, stratified for lack of HLA-B*07:02/39:10/42:02/81:01 coexpression. Targeting frequency correlated with response magnitude (r = +0.96, p < 0.0001; Supplemental Fig. 1) and conformed to a clear pattern of immunodominance (Fig. 1A).
Identification of immunodominant and subdominant HLA-B*42:01–restricted epitopes. (A) Targeting frequencies for “protein–18-mer peptide-optimal epitope name” OLPs in HLA-B*42:01+ individuals, stratified for lack of HLA-B*07:02/39:10/42:02/81:01 coexpression (n = 154). Responses were determined via IFN-γ ELISPOT assays. (B) Correlation between HLA-B*42:01 status and targeting of OLP-275 (KVVPRRKAKIIKDYGKQM) in IFN-γ ELISPOT assays. Significance was determined using the Fisher exact test. (C) Identification of the optimal HLA-B*42:01–restricted epitope Int-IM9 (IIKDYGKQM) via peptide truncations in IFN-γ ELISPOT assays. (D) Confirmation of HLA-B*42:01 as the restriction element for Int-IM9 via peptide pulsing of BLCLs partially HLA matched to donor R014 (A*30:01/33:01, B*42:01/57:03, Cw17:01/18:01). Autologous or mismatched BLCLs were used as positive and negative controls, respectively. (E) Unequivocal confirmation of HLA-B*42:01 as the restriction element for Int-IM9 via cognate tetramer staining of responding PBMCs. Similar data are shown for three other novel epitopes: Gag-RM9 (RPGGKKHYM), Vif-HI10 (HPKVSSEVHI), and Env-IF9 (IPRRIRQGF). An HLA-mismatched tetramer was used as the negative control.
To validate the novel HLA-B*42:01–restricted CD8+ T cell epitopes, we first mapped Int-IM9 (IIKDYGKQM) as the optimal peptide within OLP Int-275 (KVVPRRKAKIIKDYGKQM), which was targeted at a frequency of 53% (Fig. 1B–D). Of note, Int-IM9 is the first example of a B7 superfamily–restricted epitope that does not contain proline at position 2, which acts as the primary anchor residue. We then used HLA-B*42:01 tetramers (16, 63) to rapidly and unequivocally identify an additional three optimal epitopes from OLP Gag-3 (Gag-RM9, RPGGKKHYM), OLP Vif-407 (Vif-HI10, HPKVSSEVHI), and OLP Env-401 (Env-IF9, IPRRIRQGF) (Fig. 1E, Table I). Finally, we demonstrated that 10 of the novel epitopes listed in Table I bound strongly to HLA-B*42:01 (KD, 2–82 nM; binding half-life, 1.1–22.4 h) (Supplemental Table I).
The identification of 20 HLA-B*42:01–restricted epitopes with predictable targeting patterns provided a unique opportunity to probe the biological impact and mechanistic basis of immunodominance patterns in HIV-1 infection.
Gag-specific responses are associated with lower viral load set-points
To establish the efficacy of each HLA-B*42:01–restricted response, we compared viral load set-points in HLA-B*42:01+ responders and nonresponders (Fig. 2). Five of the six Gag-specific responses showed a trend toward lower median viral loads, which was statistically significant for the Gag-RM9 epitope after multiple-comparisons analysis (p = 0.002). In contrast, three of the six responses directed against accessory/regulatory proteins (two Nef-specific and one Vif-specific) showed a trend toward higher median viral loads, although no significant associations were observed after correction for multiple comparisons. Discordant associations with viral load set-point were also observed for response breadth (Fig. 3). In particular, only the number of Gag-specific responses was linked to immune control in the context of HLA-B*42:01 restriction, confirming previous analyses across the entire cohort (15). Thus, epitope-specific differences discriminate HLA-B*42:01–restricted responses with respect to immune control of HIV-1 replication.
Subdominant HLA-B*42:01–restricted Gag-specific responses are associated with lower viral loads. Viral load set-points were compared across OLP-specific CD8+ T cell responses targeting Gag (A), Pol (B), Nef (C), or Vpr/Vif/Env (D) in HLA-B*42:01+ responders and nonresponders, stratified for lack of HLA-B*07:02/39:10/42:02/81:01 coexpression (n = 126). Horizontal lines indicate median values. All p values < 0.1, by the Dunn multiple-comparisons test, are shown.
Discordant viral load associations with HLA-B*42:01–restricted protein-specific responses. Viral load set-points were compared across groups making 0, 1, 2, and >2 protein-specific responses directed against Gag (A), Pol (B), Nef (C), or Vpr/Vif/Env (D) in the cohort of HLA-B*42:01+ individuals, stratified for lack of HLA-B*07:02/39:10/42:02/81:01 coexpression (n = 126). Horizontal lines indicate median values. All p values < 0.05, by the Dunn multiple-comparisons test, are shown.
Gag-specific selection pressure is associated with loss of antiviral efficacy
The ability of certain CD8+ T cell responses to exert antiviral selection pressure has previously been associated with immune efficacy (15–17, 44). To examine this phenomenon across HLA-B*42:01–restricted specificities, we initially focused on the Gag-RM9 response, which exhibited the strongest association with viremic control (Fig. 2A). Escape mutations were selected in 92% of HLA-B*42:01+ subjects compared with only 54% of HLA-B*42:01− individuals (p = 2 × 10−35) (Fig. 4A). These polymorphisms were associated with higher viral load set-points (median HIV RNA copies/ml plasma = 3,860 versus 18,550, p = 0.003) (Fig. 4B) and lower CD4+ T cell counts (median CD4+ T cells/μl = 467 versus 328, p = 0.03) (Fig. 4C) in HLA-B*42:01+ individuals. Given the absence of such associations in HLA-B*42:01− subjects, it is likely that these effects operate via evasion of the Gag-RM9 response. Moreover, a beneficial effect on viral load set-point in Gag-RM9 responders carrying the wild-type epitope was observed across 60 HLA-B*42:01+ individuals for whom IFN-γ ELISPOT and viral sequence data were available from the same time point (median HIV RNA copies/ml plasma = 387 versus 2470, p = 0.04) (Fig. 4D). These findings suggest that the Gag-RM9 response contributes to viremic control most effectively in the absence of viral escape, but also to some extent in the presence of mutations that compromise CD8+ T cell recognition (Fig. 4E).
HIV-1 sequence polymorphisms in Gag-RM9 are CD8+ T cell escape mutations associated with loss of immune control. (A) HLA-B*42:01+ (red) and HLA-B*42:01− (green) associations with mutations in Gag-RM9 (departure from wild-type). Significance was determined using the Fisher exact test. Viral load set-points (B) and CD4+ T cell counts (C) in HLA-B*42:01+ (red) and HLA-B*42:01− (green) individuals carrying either wild-type Gag-RM9 (RPGGKKHYM) or any mutation in this epitope. The x-axis shows the number of sequences analyzed in each case. (D) Viral load set-points in 60 HLA-B*42:01+ individuals carrying either wild-type Gag-RM9 or any mutation in this epitope stratified for responder or nonresponder status. (E) Impact of the commonly selected H28R, M30R, and M30K variants on CD8+ T cell recognition in IFN-γ ELISPOT assays. Data from subject R019 (A*02:01/30:01, B*35:01/42:01, Cw*16:01/17:01) are shown. (F) HLA-B*42:01 binding affinities and half-lives (stability) for the peptides shown in (E). In (B)–(D), horizontal lines indicate median values presented by 25th percentile boxes and ranges (SD); significance was determined using the Mann–Whitney U test.
In further analyses, we examined the relationship between viral sequences across the whole HIV-1 genome and HLA-B*42:01 expression in a total of 1867 individuals. Four polymorphisms in the Gag-RM9 epitope were selected by HLA-B*42:01 (H28R, H28S, M30R, and M30K) (Table II). These mutations resulted in reduced CD8+ T cell recognition (Fig. 4E), at least in part through decreased peptide–HLA-B*42:01 binding (Fig. 4F). In total, we identified 20 polymorphic sites located in 9 of the 20 HLA-B*42:01–restricted epitopes (q < 0.05), including the novel Gag-GF9 epitope (Table II). Of the 10 top-ranking immunodominant HLA-B*42:01–restricted responses, all six targeting epitopes in Gag, Pol, or Vpr showed evidence of selection pressure on the virus. In contrast, only one of the four responses targeting epitopes in Nef or Vif showed evidence of selection pressure. These data provide further evidence of variable antiviral efficacy across the 20 HLA-B*42:01–restricted responses, most notably those targeting immunodominant epitopes.
TCR bias in the HLA-B*42:01–restricted repertoire
The presence of public TCR clonotypes, defined on the basis of TRB amino acid sequence identity across multiple individuals responding to the same epitope (36, 38, 64), has previously been linked with CD8+ T cell efficacy in SIV (33, 34), HIV (30, 65), and CMV (35) infection. To examine this phenomenon in the setting of HLA-B*42:01–restricted responses to HIV-1, we sorted viable tetramer+ CD8+ T cell populations (n = 48) specific for six of the seven most frequently targeted epitopes (Gag-TL9, Int-IM9, Vpr-FL9, Nef-TL10, Vif-HI10, and Int-LI9) (Supplemental Fig. 2) and sequenced a total of 3592 constituent TCR clones across eight donors/specificity (Tables III, IV, V, VI, VII, and VIII). Clonotypic data for the Gag-RM9 response were only available from one subject (Supplemental Table II) and could not be included in further analyses.
For the top two ranking immunodominant responses, Gag-TL9 (Table III) and Int-IM9 (Table IV), public TCRs were identified in six of eight and four of eight subjects, respectively. The most prevalent examples in terms of recurrence were TRBV12-3/4/CASSFSKNTEAF/TRBJ1-1 (Gag-TL9) and TRBV9/CASSVDKGGTDTQY/TRBJ2-3 (Int-IM9), each of which were shared by four of eight individuals. Public clonotypes were also identified for Vpr-FL9, ranked third in the immunodominance hierarchy (Table V). Each specificity was characterized by diverse, but distinct, patterns of TRBV usage (Fig. 5A). Moreover, within each specificity, almost all public clonotypes used the same TRBV gene (Fig. 5B). For example, the four Int-IM9–specific public clonotypes all expressed TRBV9 (Table IV), suggesting “hard-wired” germline-encoded Ag recognition.
Clonality and TCR bias in HLA-B*42:01–restricted epitope-specific responses. (A) TRBV usage is depicted for six HLA-B*42:01–restricted epitope-specific responses. The number of individuals analyzed is indicated in each case. (B) Percentage of responders to Gag-TL9, Int-IM9, or Vpr-FL9 stratified for the presence of public TCRs. Significance was determined using the Fisher exact test. (C) Number of different TCR clonotypes detected in each of the six epitope-specific responses. All p values < 0.1, by the Dunn multiple-comparisons test, are shown. (D) Correlation between clonality and response magnitude. Significance was determined using the Spearman rank test.
Overall, we identified 11 different public clonotypes, all of which were confined exclusively to the top three ranked immunodominant specificities (Gag-TL9, Int-IM9, and Vpr-FL9) (Tables III–V). In contrast, no public clonotypes were present within the Nef-TL10, Vif-HI10, and Int-LI9 specificities, ranked fifth to seventh in the immunodominance hierarchy (Tables VI–VIII). These observations suggest that clonotypic publicity is a feature of immunodominant HIV-specific CD8+ T cell responses restricted by HLA-B*42:01. The Vpr-FL9–specific CD8+ T cell populations were the most polyclonal (Fig. 5C), but no correlation was found between clonality and response magnitude (p = 0.2) (Fig. 5D). Almost identical results were obtained after normalization procedures that account for differences in sampling depth between individuals (data not shown).
Public clonotypes are involved in immunodominant CD8+ T cell responses
A stable peptide–MHC interaction is a prerequisite for immunogenicity, but other factors also contribute (12). Based on our observations above, we hypothesized that TCR bias may play a role in immunodominance. Accordingly, we compared measures of public TCR occurrence with epitope targeting across all HLA-B*42:01+ individuals (n = 154) (Fig. 6A). Strong correlations were observed for immunodominance rank (r = −0.94, p < 0.0001) (Fig. 6B) and targeting frequency (r = +0.92, p = 0.03) (Fig. 6C), both of which held whether we compared the frequency of individuals sharing any public TCR or the frequency of individuals sharing the most common public TCR. The latter correlation was even more marked after normalization for sampling depth (r = +0.94, p = 0.005; data not shown). Although we deliberately analyzed an identical number of CD8+ T cell populations from an identical number of individuals for each specificity (n = 8), this result is important because it indicates that the detected association is not an artifact of differential sequence coverage. An additional correlation was observed for response magnitude (Fig. 6D). Moreover, exclusively private responses were significantly less dominant than responses incorporating public TCRs (p = 0.04) (Fig. 6E).
Epitope-targeting frequency correlates with shared TCR recruitment. (A) Targeting frequencies for six optimal epitopes in HLA-B*42:01+ individuals, stratified for lack of HLA-B*07:02/39:10/42:02/81:01 coexpression (n = 154). (B) Immunodominance rank, shown as most targeted ranked 1, versus the percentage of individuals sharing any public TCR. Significance was determined using the Spearman rank test. (C) Targeting frequencies for the same six epitopes versus the percentage of individuals sharing the most frequent public TCR. Significance was determined using the Spearman rank test. (D) Response magnitude, including responders and nonresponders (mean), versus the percentage of individuals sharing the most frequent public TCR. Significance was determined using the Spearman rank test. (E) Targeting frequencies for all six epitopes stratified for the presence or absence of public TCRs. Significance was determined using the Mann–Whitney U test.
Type 1 bias, defined by TRBV sharing alone, was not directly linked to immunodominance (p = 0.2; data not shown). Similarly, there was no association between type IV bias, defined by near-identical (disparity < 2 aa) CDR3 sequences (36), and targeting frequency (p = 0.7; data not shown). It is also notable that response magnitude did not correlate with the frequency of public TCRs (p = 0.74; data not shown). Thus, public TCR recruitment is linked to increased targeting frequency independently of nonidentical bias and response magnitude.
Peptide–HLA-B*42:01 binding contributes minimally to immunodominance
To assess the relative impact of TCR bias with respect to other factors that shape immunodominance (6, 16, 66), we examined peptide–HLA-B*42:01 binding affinity and stability (52, 53) for 16 of the epitopes described in this study (Supplemental Table I). The two immunodominant epitopes, Gag-TL9 and Int-IM9, bound HLA-B*42:01 with affinities of KD = 82 and 19 nM, respectively, and corresponding half-lives of 3.6 and 13.5 h (Fig. 7A, 7B). Across all epitopes, however, no significant correlations were detected between these binding parameters and immunodominance (Fig. 7C, 7D). Nonetheless, peptide–HLA-B*42:01 binding half-life correlated positively with targeting frequency when Gag-TL9 was excluded from the analysis (r = +0.53, p = 0.04; data not shown). Neither binding affinity nor stability correlated with response magnitude (p = 0.42 and p = 0.32, respectively; data not shown). Thus, peptide binding alone cannot explain the observed patterns of immunodominance in this system.
Peptide–HLA-B*42:01 binding affinities and half-lives for 16 optimal epitopes. HLA-B*42:01 binding affinities (A) and half-lives (B) for the Gag-TL9 (TPQDLNTML) and Int-IM9 (IIKDYGKQM) peptides. Peptide binding affinities (C) and half-lives (D) for the optimal epitopes versus targeting frequencies for the 16 corresponding OLPs in HLA-B*42:01+ individuals, stratified for lack of HLA-B*07:02/39:10/42:02/81:01 coexpression (n = 154). Significance was determined using the Spearman rank test.
Discussion
In this study, we conducted an extensive analysis of CD8+ T cell responses restricted by a single HLA class I molecule to evaluate the mechanisms that contribute to immunodominance and antiviral efficacy in HIV-1 infection. Epitope-targeting frequencies conformed to a clear hierarchy across 20 different HLA-B*42:01–restricted specificities, within which the presence and multiplicity of responses directed against Gag were consistently associated with lower viral load set-points. Immunodominance patterns correlated with epitope-specific public TCR usage, but no clear role was identified for peptide–HLA-B*42:01 binding affinity or stability in the observed response profile. Collectively, these results suggest that the available TCR repertoire can influence immunodominance patterns and CD8+ T cell efficacy in chronic HIV-1 infection.
Initially, we used OLPs spanning the entire viral proteome to screen a large cohort of individuals with chronic C-clade infection for HIV-specific CD8+ T cell responses. This comprehensive and unbiased strategy allowed the identification of novel and previously defined HLA-B*42:01–restricted epitopes that adhered to a strict pattern of immunodominance, thereby enabling more detailed downstream analyses. It is notable that the accurate quantification of certain responses can be compromised by this approach, especially when the optimal epitope resides in the central part of the corresponding OLP (18, 67). Nonetheless, we found a strong correlation between optimal peptide-specific and OLP-defined responses (r = +0.85, p < 0.0001; data not shown). Although this finding justifies the use of OLP screening, it remains conceivable that minor discrepancies could negatively influence the observed correlation between response magnitude and targeting frequency.
Subsequent analyses revealed a significant association between HLA-B*42:01–restricted Gag-specific CD8+ T cell responses and lower viral load set-points, consistent with the notion that protein-targeting patterns influence immune efficacy. Mechanistically, this may be a function of the rapid processing kinetics and relative abundance of Gag-derived epitopes on the cell surface, enabling the elimination of infected targets prior to the production of viral progeny (57, 68, 69). In addition, Gag targeting may be beneficial due to the fitness costs incurred by viral escape mutations within highly conserved, functionally constrained regions of the viral proteome (25–27, 70, 71). However, not all escape mutations reduce viral replicative capacity, as exemplified by the Gag-RM9 variants associated with higher viral loads in our study. This observation is consistent with a recent analysis of the M30R polymorphism, which was shown to increase viral replicative capacity and precipitate disease progression in association with HLA-B*42:01 (72).
The detected association between TCR publicity and immunodominance is particularly striking given the vast recombinatorial diversity that shapes the available repertoire in each individual (36, 38, 64, 65). Across a total of 48 ex vivo datasets, however, TCR bias was apparent at multiple levels. In particular, each epitope-specific CD8+ T cell population displayed distinct patterns of TRBV usage. Moreover, the extent of TCR sharing varied as a function of epitope specificity. In line with previous studies, clonal diversity per se did not correlate with immunodominance (73). However, the most frequently targeted epitopes mobilized cognate repertoires characterized by the presence of high-frequency public TCRs. It is established that public TCRs arise in the naive repertoire as a function of convergent recombination, whereby production frequency is dictated on a probabilistic basis by the number of rearrangements, nucleotide additions, and amino acid codons that can generate a particular sequence (64, 74). The subsequent recruitment of these TCRs into the epitope-specific memory repertoire is governed by Ag avidity and clonal proliferation (10, 55). Accordingly, our observations suggest that commonly targeted epitopes are structurally compatible with frequently generated TCRs, at least in the context of HLA-B*4201. In this scenario, public TCRs populate a recurrence “hot spot” in the naive repertoire, enabling immunodominant epitopes to initiate and maintain CD8+ T cell responses in multiple individuals. Consistent with this proposition, naive precursor frequencies are known to influence immunodominance patterns (6, 7, 9, 75, 76). However, the extent to which frequently produced TCRs contribute to the overall precursor pool for any given epitope remains to be determined.
It is noteworthy that the public TCRs detected for each of the three immunodominant specificities were largely constrained by a distinct TRBV segment. This suggests a determinative role for the germline-encoded CDR1 and CDR2 loops. One possibility is that specific residues encoded by these TRBV genes interact with HLA-B*42:01, potentially influencing thymic selection to skew naive CD8+ T cell frequencies toward particular specificities (77). Alternatively, germline-encoded recognition of the bound peptide may contribute to immunodominance patterns (78). Immunodominant epitope-specific TCRs constructed almost entirely from germline DNA have been described previously and may represent an evolutionarily conserved mechanism to combat ancient pathogens that have coevolved with the human race (37, 64, 79, 80). Although this is unlikely to apply directly in the case of HIV-1, it is intriguing to speculate that structural homology with such epitopes may inadvertently underlie the immunodominance patterns described in this report.
In contrast to the association between immunodominance and TCR bias controlled within the framework of a single HLA class I molecule, peptide binding to HLA-B*42:01 contributed little to the observed epitope-targeting frequencies in this study. However, it is important to note that other factors, such as Ag-processing efficiency (81, 82), kinetics (57, 69, 83), and protein abundance (84, 85), play a key role alongside a requirement for sufficient peptide–MHC binding (9, 86). Further studies will therefore be necessary to assess the contribution of TCR recruitment in relation to these well-defined determinants of immunodominance.
In summary, we present clear evidence linking epitope-targeting frequencies to TCR bias. Although the extent to which this phenomenon applies across other systems remains to be defined, our data suggest that Ag-specific repertoire studies will be important for a full understanding of both natural and vaccine-induced immune responses against intracellular pathogens.
Disclosures
The authors have no financial conflicts of interest.
Footnotes
This work was supported by the Wellcome Trust (to D.A.P. and P.J.G.) and National Institutes of Health Grant R01 AI46995. H.N.K. is funded by the Danish Agency for Science, Technology and Innovation via Grant 12-132295, The Lundbeck Foundation via Grant R151-2013-14624, and The MAERSK Foundation for Medical Improvement. D.A.P. is a Wellcome Trust Senior Investigator.
The online version of this article contains supplemental material.
Abbreviations used in this article:
- BLCL
- B lymphoblastoid cell line
- OLP
- overlapping peptide.
- Received April 3, 2014.
- Accepted February 25, 2015.
- Copyright © 2015 by The American Association of Immunologists, Inc.