key: cord-0852794-a3mgjo07
authors: Ahmed, Syed Faraz; Sohail, Muhammad Saqib; Quadeer, Ahmed Abdul; McKay, Matthew R.
title: Identification of Potential SARS-CoV-2 CD8(+) T Cell Escape Mutants
date: 2022-03-31
journal: Vaccines (Basel)
DOI: 10.3390/vaccines10040542
sha: 39eb796881d710748c253ff3256331002b67ef05
doc_id: 852794
cord_uid: a3mgjo07

Memory SARS-CoV-2-specific CD8(+) T cell responses induced upon infection or COVID-19 vaccination have been important for protecting against severe COVID-19 disease while being largely robust against variants of concern (VOCs) observed so far. However, T cell immunity may be weakened by genetic mutations in future SARS-CoV-2 variants that lead to widespread T cell escape. The capacity for SARS-CoV-2 mutations to escape memory T cell responses requires comprehensive experimental investigation, though this is prohibited by the large number of SARS-CoV-2 mutations that have been observed. To guide targeted experimental studies, here we provide a screened list of potential SARS-CoV-2 T cell escape mutants. These mutants are identified as candidates for T cell escape as they lie within CD8(+) T cell epitopes that are commonly targeted in individuals and are predicted to abrogate HLA–peptide binding.

SARS-CoV-2-specific neutralizing antibodies and T cells are generated in response to COVID-19 infection and vaccines. The level of immunity that they provide is challenged by the evolution of SARS-CoV-2 and the rise of new genetic variants such as Omicron (see [1] for details on SARS-CoV-2 variants). This presents a key issue for neutralizing antibodies for which genetic mutations enabling escape have been widely reported (e.g., [2] [3] [4] ). For SARS-CoV-2-specific T cell responses, multiple studies have suggested that they are largely robust to mutations observed in the variants of concern (VOCs) so far [5] [6] [7] [8] . However, in a small fraction of recovered and vaccinated individuals, loss of CD8 + T cell responses has been reported against variants such as Omicron [8, 9] . Moreover, for some specific CD8 + T cell epitopes, it has been reported that mutations have weakened or abrogated T cell recognition in some individuals [10, 11] . These recent immunological studies [8] [9] [10] [11] indicate at least some capacity for SARS-CoV-2 to escape T cell responses.

CD8 + T cell responses present an important weapon against COVID-19, having been associated with limiting disease severity in infected [12] as well as vaccinated individuals [13] .

Mutations that facilitate T cell escape could potentially compromise the protection provided by COVID-19 vaccines and infection-induced immunity against severe disease [14, 15] . However, the degree to which circulating SARS-CoV-2 viruses may evade CD8 + T cell responses is not yet well understood and requires detailed experimental investigation. Here, we report results to aid such studies by identifying and screening mutations observed among sequenced SARS-CoV-2 viruses that may facilitate CD8 + T cell escape.

We obtained SARS-CoV-2 CD8 + T cell epitope data from the dashboard reported in [5] (https://www.mckayspcb.com/SARS2TcellEpitopes/; accessed on 1 February 2022). This dashboard provides up-to-date information on SARS-CoV-2 epitopes targeted by T cells in COVID-19 recovered individuals. The downloaded data comprised a total of 753 distinct HLA-specific CD8 + T cell epitopes.

SARS-CoV-2 sequence data were obtained from the GISAID database (https://www. gisaid.org/) accessed on 6 September 2021. Only the complete (full genome) sequences derived from human hosts with high coverage were downloaded. All (2,436,276) downloaded sequences were aligned to the SARS-CoV-2 reference genome (GenBank: NC_045512.2) using MAFFT [16] . The aligned genomic sequences were then translated using an in-house code to obtain the protein sequences. The positions of the open reading frames provided with the reference sequence were used to identify the respective protein regions of the full genome.

For each SARS-CoV-2 epitope, we screened the processed sequence data at the epitope location and defined an epitope mutant as any epitope sequence which differed from the reported sequence of that epitope. Epitope mutants comprising any ambiguous amino acid were discarded. The count of each identified epitope mutant, i.e., the number of sequences comprising the epitope mutant, was also recorded.

To understand if an identified epitope mutant impacts the binding of the associated epitope with its cognate HLA allele, we used the well-known peptide-HLA binding prediction method, NetMHCpan4.1 [17] . All predictions were made using the default parameters of the prediction method.

To identify potential CD8 + T cell escape mutants, we analyzed mutations lying within the 753 experimentally determined SARS-CoV-2 epitopes based on~2.5 million SARS-CoV-2 sequences (Materials and Methods). This analysis revealed over 35,000 epitope mutants across proteins. Comprehensive experimental testing of these mutants for their potential to evade CD8 + T cell responses is prohibitive, necessitating further screening. A more refined set was obtained by considering those mutants that are associated with epitopes that are commonly targeted globally and which thereby have the potential to result in widespread T cell escape. A list of 20 such immunoprevalent SARS-CoV-2 CD8 + T cell epitopes was reported in [5] , which are collectively estimated to be targeted by at least 85% of the global population based on their HLA alleles association [18] . Individually, the global coverage of the HLA alleles associated with the immunoprevalent epitopes ranged from 13% to 39%. An important feature of these epitopes is that T cell responses against them were recorded in multiple studies across different geographical locations and in more than half of the COVID-19 recovered individuals tested [5] . These 20 immunoprevalent T cell epitopes have been remarkably conserved so far and do not comprise any mutation specific to VOCs observed to date [7] . Restricted to this immunoprevalent epitope set, the number of HLA-specific epitope mutants was reduced to 914.

To determine the likely escape mutations within immunoprevalent epitopes, a further screening step can be performed using tools to predict the impact of mutations on peptide-HLA binding (a prerequisite for T cell recognition [19] ). For this purpose, we employed NetMHCpan4.1 [17] , a computational method whose predictions have lined up well with SARS-CoV-2 experimental studies describing CD8 + T cell epitopes [20] , as well as those reporting HLA binding of specific T cell epitopes [10, 11] .

Examination of each mutant of the immunoprevalent epitopes revealed that the majority (744/914) were predicted to have no effect on HLA binding, while very few increased binding (4/914). However, almost one-fifth (166/914) of the epitope mutants were predicted to decrease binding, and hence, these presented candidates for enabling escape from T cell responses against commonly targeted SARS-CoV-2 CD8 + T cell epitopes. We recommend that these 166 epitope mutants (listed in Supplementary Material Table S1 ) be considered for experimental investigation. This may involve conducting assays such as targeted HLA binding or HLA ligand elution assays, in combination with standard T cell assays such as multimer qualitative binding, enzyme-linked immunospot (ELISPOT), or intracellular cytokine staining (ICS) interferon-γ (IFN-γ) release, and activation-induced markers (AIM) assays.

Among the identified 166 epitope mutants, most have been observed so far at only low frequencies within the global SARS-CoV-2 sequence data. Only 83 epitope mutants have been observed five times or more (listed in Table 1 ). Among these, about half (42/83) are associated with HLA-A*02:01, the most prevalent HLA allele globally. Of the 83 epitope mutants, 16 are derived from the spike (S) protein, which could be particularly important given that most COVID-19 vaccines target the S antigen. These also include a mutant of a notable immunoprevalent epitope, YLQPRTFLL, which has been shown to be targeted during COVID-19 infection in more than 10 different cohorts so far [5] . The remaining 67/83 epitope mutants are derived from SARS-CoV-2 proteins besides S, which are of interest for T cell responses induced by inactivated COVID-19 vaccines or prior infection. 1 Reported as immunoprevalent in [5] . The numbers in the subscript before and after the epitope sequence are the start and end indices of the epitope within the corresponding protein. 2 An epitope mutant represents an experimentally determined CD8 + T cell epitope with at least one mutation (shown in bold) observed in its sequence. The listed epitope mutants are predicted (using the default parameters of NetMHCpan4.1) to decrease HLA binding and occur at least five times in SARS-CoV-2 global sequence data (~2.5 million). The epitope mutants belonging to a specific protein are grouped together. Within each group, the epitope mutants are presented in descending order according to the number of times they are observed in the global sequence data.

The approach that we have adopted, which screens potential T cell epitope mutants based on their ability to evade peptide-HLA binding, is a conservative strategy. This is because a loss of peptide-HLA binding is irreparable since such a peptide would not be expressed on the cell surface for recognition by T cells. It is also possible that epitope mutants that do not decrease HLA binding may still result in T cell escape, mediated by other mechanisms such as abrogation of T cell receptor (TCR) binding with the peptide-HLA complex [21] . However, T cell escape via this means may be recovered by TCR recognition by other T cell clonotypes present in an individual [22] . Thus, while T cell escape through abrogation of TCR binding should not be ignored, the epitope mutants that we have identified to decrease peptide-HLA binding appear to be clear candidates for conferring T cell escape.

SARS-CoV-2 may evolve mutations to broadly escape T cell responses. While memory T cell responses induced by prior infection or vaccination have so far mainly remained robust [5] [6] [7] [8] , some evidence of T cell escape has been reported [8] [9] [10] [11] . Continued monitoring of the capacity of SARS-CoV-2 to broadly evade T cell responses is required, as widespread T cell escape could compromise population-level immunity and have significant clinical and public health implications. To aid such monitoring efforts, our study considered SARS-CoV-2-specific CD8 + T cell responses and sought to identify potential SARS-CoV-2 CD8 + T cell escape mutants.

A comprehensive experimental evaluation of all potential SARS-CoV-2 CD8 + T cell escape mutants would involve testing all relevant mutations observed in available genetic sequence data, which is prohibitive as it would involve evaluating tens of thousands of mutations within known SARS-CoV-2 epitopes. To help address this challenge, we performed a computational analysis to facilitate targeted experimental screening of T cell escape mutants. The main outcome of our study, given in Table 1 , is a selective list of 83 observed SARS-CoV-2 mutations that lie within commonly targeted CD8 + T cell epitopes and have been observed in the SARS-CoV-2 global sequence data. An important feature of these epitope mutants is that they are expected to decrease peptide-HLA binding, which is a prerequisite for T cell recognition. Escape from these epitopes could be problematic, as it could affect vaccine and infection-induced T cell immunity in a large fraction of the population.

Overall, our computational study leverages available SARS-CoV-2 sequence and immunological data to direct attention towards commonly targeted T cell epitope mutants that emerge as clear candidates to confer escape from the vaccine or infection-induced responses. These epitope mutants are recommended for immediate experimental investigation.

Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/vaccines10040542/s1, Table S1 (included as a separate file): List of SARS-CoV-2 immunoprevalent CD8 + T cell epitope mutants predicted to decrease HLA binding. 

The authors confirm that the data supporting the findings of this study are available within the article or its Supplementary Materials.

SARS-CoV-2 infection: New molecular, phylogenetic, and pathogenetic Insights. Efficacy of current vaccines and the potential risk of variants

New SARS-CoV-2 variants-clinical, public health, and vaccine implications

Considerable escape of SARS-CoV-2 Omicron to antibody neutralization

Sensitivity of infectious SARS-CoV-2 B.1.1.7 and B.1.351 variants to neutralizing antibodies

Landscape of epitopes targeted by T cells in 852 individuals recovered from COVID-19: Meta-analysis, immunoprevalence, and web platform

Impact of SARS-CoV-2 variants on the total CD4+ and CD8+ T cell reactivity in infected or vaccinated individuals. Cell Rep

SARS-CoV-2 T cell responses elicited by COVID-19 vaccines or infection are expected to remain robust against Omicron

T cell responses to SARS-CoV-2 spike cross-recognize Omicron

T-cell reactivity to the SARS-CoV-2 Omicron variant is preserved in most but not all individuals

The COVID-19 Genomics UK (COG-UK) consortium; et al. Emergence of immune escape at dominant SARS-CoV-2 killer T-cell epitope

SARS-CoV-2 mutations in MHC-I-restricted epitopes evade CD8+ T cell responses

Antigen-specific adaptive immunity to SARS-CoV-2 in acute COVID-19 and associations with age and disease severity

Early T cell and binding antibody responses are associated with COVID-19 RNA vaccine efficacy onset

The T cell immune response against SARS-CoV-2

Population impact of SARS-CoV-2 variants with enhanced transmissibility and/or partial immune escape

MAFFT multiple sequence alignment software Version 7: Improvements in performance and usability

NetMHCpan-4.1 and NetMHCIIpan-4.0: Improved predictions of MHC antigen presentation by concurrent motif deconvolution and integration of MS MHC eluted ligand data

Predicting population coverage of T-cell epitope-based diagnostics and vaccines

T cell receptor recognition of self and foreign antigens in the induction of autoimmunity

In silico T cell epitope identification for SARS-CoV-2: Progress and perspectives

T cell epitope predictions

A molecular basis for the control of preimmune escape variants by HIV-Specific CD8+ T cells

We acknowledge the contributions of all the researchers, scientists and technical staff involved in the open sharing of SARS-CoV-2 sequence data through the GISAID platform.

The authors have filed for patent protection for various aspects of SARS-CoV-2 T cell epitope and vaccine design work.