T-cell immunity plays an important role in augmenting vaccine-mediated protection against SARS-CoV-2, particularly in the face of emerging variants such as Omicron1,2, and is likely to be critically important in individuals who fail to generate efficient neutralizing antibody responses due to underlying disease 3,4. Following the emergence of Omicron, studies have demonstrated that most T-cell epitopes are conserved 5,6, and the majority of vaccine recipients maintain T-cell immunity.7-9 Nevertheless, evidence in convalescent individuals has demonstrated the capacity for variant strains to evade T cell responses10. In this study we sought to investigate the impact of changes in Omicron on a uniquely immunodominant spike-encoded CD8+ T-cell epitope.
Our group recruited 59 COVID-19 convalescent participants during the initial 2020 wave of the pandemic in Australia. All participants were typed for human leukocyte antigens (HLA, Supplementary Table 1) and SARS-CoV-2-specific T cells were expanded from peripheral blood mononuclear cells (PBMC) using SARS-CoV-2-encoded antigens 11. In this analysis, we noted striking immunodominant CD8+ T cell responses to both spike (S-1 pool) and ORF3A in HLA-A*29:02+ participants (Figure 1A, Supplementary Figure 1). We detected a median S-1 pool-specific CD8+ T-cell response of 44.6% in HLA-A*29:02+ participants compared to 0.85% in HLA-A*29:02- participants (Figure 1B). The ORF3A-specific CD8+ T-cell response was detected at a median of 31.2% in HLA-A*29:02+ participants compared to 1.47% in HLA-A*29:02- participants. Using an overlapping peptide matrix, we identified that the spike-specific T-cell response targeted spike matrix pools 5, 6 and 23 (Figure 1C), corresponding to overlapping peptides GFNCYFPLQSYGFQP and YFPLQSYGFQPTNGV (Supplementary Figure 2A). Peptide minimization and HLA restriction (Supplementary Figure 2B) revealed the immunodominant response to be encoded by a HLA-A*29:02-restricted epitope, YFPLQSYGF (Figure 1D). Similar analysis revealed an immunodominant response in matrix pools 7, 8 and 15 from ORF3a (Figure 1E), corresponding to overlapping peptides VVLHSYFTSDYYQLY and SYFTSDYYQLYSTQL (Supplementary Figure 2C), which contain the previously published HLA-A*29:02-restricted epitope YFTSDYYQLY (Figure 1F) 12.
Analysis of the amino acid sequences of the HLA-A*29:02 epitopes in SARS-CoV-2 variants revealed that while the ORF3A-encoded epitope was conserved, the Omicron variant contained amino acid substitutions at positions 5 and 8 in the spike-encoded epitope (Figure 1G). Although not anchor residues for A*29:02, the substitutions, particularly the non-conserved glutamine (Q) to arginine (R) change at position 5, have the potential to significantly impact T-cell receptor (TCR) engagement. To address this, we generated T cells specific for the YFPLQSYGF epitope from three COVID-19-convalescent participants and performed intracellular cytokine analysis using 10-fold serial peptide titration with both the cognate epitope and Omicron variant. Despite strong functional avidity directed towards the cognate peptide, we saw poor reactivity against the Omicron variant (Figure 1H, Supplementary Figure 3).
We next sought to determine if individuals vaccinated against SARS-CoV-2 would show a similar pattern in their response to YFPLQSYGF. From a cohort of 104 participants, recruited following vaccination with either Vaxzevria or Comirnaty, we identified three participants with HLA-A*29:02 (Supplementary Table 2). PBMC from these individuals were isolated 28 days after their second vaccine dose, then stimulated with a spike peptide pool, the YFPLQSYGF epitope or the Omicron variant. We assessed expression of IFN-g and TNF by CD8+ T cells using intracellular cytokine analysis to identify low-frequency T-cell responses. Controls were unstimulated PBMC (no peptide) and PBMC incubated with a cell stimulation cocktail (CSC). Although low in abundance, IFN-g+TNF+ CD8+ T cell responses against the spike overlapping peptide pools (0.006%, 0.009%, and 0.002% above no peptide) and the YFPLQSYGF epitope (0.046%, 0.012%, and 0.006% above no peptide) were detected in PBMC from all three individuals (Figure 2A). Similar observations were evident following analysis of IFN-g alone (Supplementary Figure 4A). Conversely, CD8+ T cells capable of recognizing the Omicron variant were not detected. To validate these observations, we generated T cells from these three participants by stimulating with the spike peptide pools and culturing in the presence of IL-2 for 14 days. In addition, we set up cultures from a fourth individual in our vaccine cohort who is HLA-A*29:11+. Spike-specific CD8+ T cell responses were observed in all four T cell expansions (Figure 2B, Supplementary Figure 4B). All four cultures were dominated by T cells specific for YFPLQSYGF, but displayed no recognition of the Omicron variant. In the three HLA-A*29:02+ cultures, the frequency of the YFPLQSYGF-specific CD8+IFN-g+TNF+ response was comparable to the response against whole spike protein, demonstrating similar immunodominance to that seen in convalescent individuals. To determine if we could generate variant-specific T cells, we stimulated PBMC with either YFPLRSYSF or the Omicron variant and assessed functional avidity after 14 days in culture. While the Omicron epitope failed to induce the expansion of T cells in these donors, YFPLRSYSF-specific T cells were detected in cultures from all four vaccine recipients (Figure 2C&D).
Given the differences in the T cell response between YFPLRSYSF and its variant, we wanted to understand the impact of the Omicron mutation on peptide presentation by HLA-A*29:02. YFPLRSYSF peptide was refolded with HLA-A*29:02 and then crystallised to understand the basis of peptide presentation (Supplementary Table 3). The HLA-A*29:02 molecule adopted the canonical fold of HLA molecules with YFPLRSYSF (10.3390/ijms22010068) (Supplementary Figure 5A). The primary anchors of the YFPLRSYSF peptide at P2 and P9 are both phenylalanine residues that are favored for HLA-A*29:02 binding (Figure 2E, Supplementary Figure 6) and in line with the thermal stability observed for the YFP-HLA-A*29:02 complex (Supplementary Table 4). The P5-Gln acts as a secondary anchor with its side chain buried in the center of the peptide binding cleft (Figure 2E), forming hydrogen bond with the Arg114 as well as hydrophobic interaction with Tyr99, Gln70, Leu156 and Met97 (Supplementary Figure 5B). The P7-Tyr is half buried with its aromatic side chain placed in between the a2-helix and the backbone of the peptide. The peptide has four residue side chains that are solvent exposed (Figure 2E), namely the P1-Tyr, P4-Leu, P6-Ser and P8-Gly that could potentially be contacted by TCRs.
The two residues of the peptide that are mutated in Omicron are at positions 5 and 8. While the substitution at P8 from a glycine to a serine is unlikely to impact peptide presentation and conformation (Figure 2F), the mutation at position 5 will have an impact on the conformation of the peptide. The P5-Gln is buried in the cleft of HLA-A*29:02, and modeling of a P5-Arg shows a steric clash of the large and charged side chain (Figure 2F) with residues in the binding cleft due to the presence of large amino acid residues Met97, Tyr99 and Arg114 that occupy the b-sheet floor of HLA-A*29:02. The P5 mutation in Omicron is likely to change peptide and could also weaken its binding to HLA-A*29:02.
To investigate the impact of amino acid changes on another spike-encoded epitope, we assessed T cell responses to the HLA-B*07:02-restricted epitope SPRRARSVA13 in five HLA-B*07:02-positive vaccine recipients (Supplementary Table 2). This epitope contains mutations detectable in multiple variants that alter the HLA-B*07:02 anchor residue at P2 (Supplementary Table 5). While these changes ablated the ability of T cells from vaccinated individuals to recognize the variant peptides, other spike-specific T cell responses were detected, including those against the conserved HLA-A*02:01–restricted epitope YLQPRTFLL12 (Figure 2G. Supplementary Figure 4C), demonstrating the complex nature of T-cell immunodominance in HLA-distinct individuals.
Despite the potentially large number of antigenic targets in the spike protein12-14, our observations suggest that the immunodominance profile in HLA-A29-positive individuals is associated with a skewed response towards a single dominant epitope. While it remains to be determined if this could impact susceptibility in vaccinated individuals, it was clear that changes in the amino acid sequence ablated T cell activation. Although HLA-A29 was present in only 6.8% and 3.8% of our convalescent and vaccinated cohort respectively, HLA-A29 frequencies have been reported as high as 24% in some populations in Africa (allelefrequencies.net, Supplementary Table 6), which could render a high proportion of these populations susceptible to evasion of CD8+ T-cell immunity following exposure to Omicron. These analyses highlight the potential impact of genetic variation in newly emerging SARS-CoV-2 variants. Demonstrating that ongoing genetic variation may lead to escape from cellular immune control, especially in ethnic groups in which HLA alleles impacted by epitope mutations are dominant.