Next Article in Journal
Therapeutic Properties of Ayahuasca Components in Ischemia/Reperfusion Injury of the Eye
Next Article in Special Issue
Human iNKT Cells Modulate Macrophage Survival and Phenotype
Previous Article in Journal
Accumulation of Arachidonic Acid, Precursor of Pro-Inflammatory Eicosanoids, in Adipose Tissue of Obese Women: Association with Breast Cancer Aggressiveness Indicators
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomedicines
Volume 10
Issue 5
10.3390/biomedicines10050996
biomedicines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Resilience of Spike-Specific Immunity Induced by COVID-19 Vaccines against SARS-CoV-2 Variants

by
Laura Ballesteros-Sanabria
†,
Hector F. Pelaez-Prestel
†,
Alvaro Ras-Carmona
† and
Pedro A. Reche
*
Laboratory of Immunomedicine, Department of Immunology & O2, School of Medicine, Complutense University of Madrid, Pza Ramon y Cajal, S/N, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
†
These authors have contributed equally to this work.
Biomedicines 2022, 10(5), 996; https://doi.org/10.3390/biomedicines10050996
Submission received: 7 April 2022 / Revised: 20 April 2022 / Accepted: 22 April 2022 / Published: 26 April 2022
(This article belongs to the Special Issue State-of-the-Art Immunology and Immunotherapy in Europe)
Browse Figures
Versions Notes

Abstract

:
The outbreak of SARS-CoV-2 leading to the declaration of the COVID-19 global pandemic has led to the urgent development and deployment of several COVID-19 vaccines. Many of these new vaccines, including those based on mRNA and adenoviruses, are aimed to generate neutralizing antibodies against the spike glycoprotein, which is known to bind to the receptor angiotensin converting enzyme 2 (ACE2) in host cells via the receptor-binding domain (RBD). Antibodies binding to this domain can block the interaction with the receptor and prevent viral entry into the cells. Additionally, these vaccines can also induce spike-specific T cells which could contribute to providing protection against the virus. However, the emergence of new SARS-CoV-2 variants can impair the immunity generated by COVID-19 vaccines if mutations occur in cognate epitopes, precluding immune recognition. Here, we evaluated the chance of five SARS-CoV-2 variants of concern (VOCs), Alpha, Beta, Gamma, Delta and Omicron, to escape spike-specific immunity induced by vaccines. To that end, we examined the impact of the SARS-CoV-2 variant mutations on residues located on experimentally verified spike-specific epitopes, deposited at the Immune Epitope Database, that are targeted by neutralizing antibodies or recognized by T cells. We found about 300 of such B cell epitopes, which were largely overlapping, and could be grouped into 54 B cell epitope clusters sharing ≥ 7 residues. Most of the B cell epitope clusters map in the RBD domain (39 out of 54) and 20%, 50%, 37%, 44% and 57% of the total are mutated in SARS-CoV-2 Alpha, Beta, Gamma, Delta and Omicron variants, respectively. We also found 234 experimentally verified CD8 and CD4 T cell epitopes that were distributed evenly throughout the spike protein. Interestingly, in each SARS-CoV-2 VOC, over 87% and 79% of CD8 and CD4 T cell epitopes, respectively, are not mutated. These observations suggest that SARS-CoV-2 VOCs—particularly the Omicron variant—may be prone to escape spike-specific antibody immunity, but not cellular immunity, elicited by COVID-19 vaccines.

1. Introduction

The detection in Wuhan city of the new coronavirus nCoV-19, now renamed as SARS-CoV-2, responsible for the ongoing coronavirus disease-19 (COVID-19) pandemic, has prompted unprecedented research efforts to characterize the virus and develop vaccines and therapeutics [1]. Coronaviruses (CoVs), such as SARS-CoV, SARS-CoV-2 and MERS-CoV, are enveloped viruses with a positive-sense, single-stranded RNA genome (ssRNA+) belonging to the Betacoronavirus genus [2]. These three CoVs infect humans causing severe respiratory illness, including pneumonia. Health crises caused by SARS-CoV and MERS-CoV in 2002 and 2012, respectively, as well as the current COVID-19 pandemic, are a reminder that new threatening CoVs are likely to emerge in the future. Therefore, it is of interest to know as much as possible about the pathogenesis and immunogenicity of these viruses so that prevention methods, antiviral treatments, such as protease inhibitors [3], antibody therapies [4] and effective vaccines can be developed.
SARS-CoV-2 has a tropism for the respiratory tract, as the spike glycoprotein binds angiotensin converting enzyme 2 (ACE2), which is highly expressed in the lung epithelium. The receptor-binding domain (RBD), a region between residues 319 and 541, is responsible for binding to the N-terminal helix of ACE2, specifically the 70 residues that comprise the receptor-binding motif inside this domain. For interaction with ACE2 to occur, the spike trimer must undergo a conformational change to expose the RBD (open state), which is buried between the trimer protomers in the close state [5]. This interaction triggers cell entry by fusion or endocytosis, which requires proteolytic cleavage of the spike by the host membrane protease TMPRSS2, or furin [6].
The spike glycoprotein plays a major role in viral transmission and thereby most current COVID-19 vaccines, including those based on the mRNA platforms (Pfizer BNT162b2 and Moderna mRNA-1273) and non-replicating viral vectors (AstraZeneca AZD-1222 and Johnson & Johnson Janssen’s Ad26.COV2.S) are directed towards the generation of neutralizing antibodies against the spike glycoprotein of the initial Wuhan sequence deposited in NCBI (NC_045512) in January of 2020. Anti-spike antibodies, particularly those recognizing the RBD, can block the binding with ACE2 and the consequent entry into the cell [7]. Subsequently, mutations in the spike protein may render neutralizing antibodies unable to recognize their cognate epitopes, thus impairing the efficacy of current vaccines. This is the case of the E484K mutation in the RBD, which is associated with the decreased neutralizing capacity of antibodies [8]. Moreover, some mutations in spike can also enhance SARS-CoV-2 pathogenicity and infectivity. For instance, the D614G mutation increases transmissibility while the N501Y mutation enhances the binding to ACE2 [9]. SARS-CoV-2 variants with such mutations pose a risk to public health [10] and are classified as variants of concern (VOCs). Up to February 2022, the currently designated VOCs are the following: B.1.1.7 originated and/or initially expanded in the United Kingdom; B.1.351 emanating from South Africa; P.1 from Brazil; the so-called “double mutant” Indian variant B.1.617.2 and the variant B.1.1.529 (likely originated in Africa) which is now predominant worldwide. These VOCs are labeled by the World Health Organization (WHO) as Alpha, Beta, Gamma, Delta and Omicron, respectively.
In this work, we evaluated to what extent SARS-CoV-2 VOCs could defy the immunity generated by current vaccines. To that end, we studied the variation of experimentally verified spike-specific epitopes deposited in the Immune Epitope Database (IEDB) [11,12] that are recognized by neutralizing antibodies or T cells. We identified 313 unique B cell epitopes—comprising 337 distinct residues of the spike protein, 172 of those in the RBD—that could be grouped into 54 clusters sharing seven or more residues. We found that the number of mutated clusters varies between variants, with Omicron being the VOC with higher variation (31 clusters with mutations). Nevertheless, the number of mutated B cell epitope clusters in SARS-CoV-2 Alpha, Beta, Gamma and Delta is still considerable (11, 27, 20 and 24, respectively). On the other hand, T cell epitopes were much less affected by mutations. In all SARS-CoV-2 VOCs above 87% and 79% of CD8 and CD4 T cell epitopes are not mutated. Overall, these results suggest that SARS-CoV-2 VOCs can elude humoral immunity elicited by current COVID-19 but not cellular immunity.

2. Materials and Methods

2.1. Protein Sequences of SARS-CoV-2 Variants

Protein sequences of SARS-CoV-2 variants were obtained after the genome sequences of the five SARS-CoV-2 VOCs: Alpha, Beta, Gamma, Delta and Omicron (Table 1). The genome sequences of these variants were identified and obtained from Global Initiative on Sharing All Influenza Data (GISAID) EpiCoVTM Database (https://www.gisaid.org/, accessed on 15 February 2022), considering only complete sequences (>29,000 bp) and excluding those entries with low coverage (>5% undefined bases) to discard poor quality sequences. To identify the mutations in the selected variants genomes, the Bioinformatic tool Nextclade v.1.14.0 was used (https://clades.nextstrain.org/, accessed on 17 Frebruary 2022). Nextclade is a Web application that can identify differences between a query and a reference sequence (SARS-CoV-2 isolate Wuhan-Hu-1, NCBI accession NC_045512). Subsequently, protein sequences of membrane spike glycoprotein (NCBI accession YP_009724390.1) for SARS-CoV-2 variants of the study were generated.

2.2. SARS-CoV-2 Spike-Specific B Cell Epitopes Targeted by Neutralizing Antibodies

B cell epitopes on the spike necessary for the subsequent analysis of the impact of the mutations were obtained through relevant searches at the IEDB (https://www.iedb.org/, accessed on 25 February 2022). To identify B cell epitopes, search parameters were restricted to human host and B cell assays that detected positive neutralization activity by antibodies against the spike glycoprotein (P0DTC2) of SARS-CoV-2 (ID:2697049) as epitope source. Both linear and discontinuous epitopes were included in the search. To simplify the analysis of B cell epitope mutations on spike glycoprotein variants, B cell epitopes encompassing between seven and 20 residues were grouped into clusters based on sequence similarity (sharing ≥ 7 residues), using a PYTHON script generated ad-hoc. The script can be obtained from the corresponding author upon written request.

2.3. SARS-CoV-2 Spike-Specific T Cell Epitopes

CD4 and CD8 T cell epitopes mapping on the spike protein were retrieved by querying the IEDB (https://www.iedb.org/, accessed on 25 February 2022). To identify CD8 T cell epitopes on SARS-CoV-2 spike glycoprotein, IEDB search parameters were restricted to infectious disease, human host, positive T cell assay and binding to MHC class I. Only those epitopes with nine residues were considered. To select CD4 T cell epitopes on SARS-CoV-2 spike glycoprotein, search parameters were restricted to infectious disease, human host, positive T cell assay and binding to MHC class II. Only epitopes larger than 11 and shorter or equal to 18 residues were considered.

2.4. Other Procedures

Multiple sequence amino acid alignments of SARS-CoV-2 spike proteins were generated using ClustalW version 2.0 [13]. The complete tertiary (3D) structure of SARS-CoV-2 spike glycoprotein was generated by homology modeling, using SWISS-MODEL (https://swissmodel.expasy.org/, accessed on 09 March 2022) with default settings. PyMOL Molecular Graphics System Version 2.4.1 Schrödinger, LLC was used to generate molecular renderings. Venn diagrams were generated using the nVennR package version 0.2.3 [14]. The percentage of the world population that could respond to CD8 and CD4 T cell epitopes attending to their HLA restriction elements was computed using EPISOPT [15] and the IEDB PPC tool (http://tools.iedb.org/tools/population/iedb_input, accessed on 16 March 2022), respectively.

3. Results

3.1. Identification of Mutations in SARS-CoV-2 Variants

We compared the genome of the reference SARS-CoV-2 with those of five SARS-Co-2 VOCs: Alpha, Beta, Gamma, Delta and Omicron. The selected genome sequences (Table 1) were obtained from the GISAID EpiCoVTM Database (https://www.gisaid.org/, accessed on 15 February 2022) after lineage-based searches. Genome sequence mutations resulting in amino acid changes were identified using Nextclade (details in Material and Methods). We identified mutations throughout the entire genome of SARS-CoV-2 VOCs but hereafter, we will focus on those located in the spike glycoprotein (Table 2).
The number of spike mutations in the studied SARS-CoV-2 variants, including amino acid changes and deletions, ranges from nine in the Delta variant to 24 mutations in the Omicron variant (Figure 1A). In total, the selected variants have 52 different residues that are mutated. Many of these mutations in SARS-CoV-2 VOCs are private but others are shared between different variants (Figure 1B). For example, the D614G mutation is found in all variants. Likewise, all SARS-CoV-2 VOCs but Delta share the N501Y mutation in the RBD. In other cases, the same residue can be affected by different mutations. For instance, Beta and Gamma include the E484K mutation, while in the same residue Delta presents the E484Q mutation and Omicron the E484A mutation.

3.2. Mutations in B Cell Epitopes Targeted by Neutralizing Antibodies

To examine the chance of SARS-CoV-2 variants escaping to spike-specific humoral immunity, we determined the occurrence of mutations on spike residues reported to be recognized by neutralizing antibodies. To that end, we identified and retrieved from the IEDB 313 distinct B cell epitopes, differing in at least one residue, which are recognized by neutralizing antibodies (search criteria in Material and Methods). All these 313 epitopes encompass 337 spike residues, 172 located in the RBD. Given that the spike and the RBD have 1273 and 223 residues, respectively, 26% of the spike and 77% of the RBD appears to be susceptible to some recognition by neutralizing antibodies. However, it is worth noting that there is a large overlap between the collected B cell epitopes and the same spike residues can be found in many distinct B cell epitopes (Figure 2A). This is exemplified by residues K417, E484, F486, N487, Y489, Q493 and N501, which are included in more than 75 different B cell epitopes. In fact, only a handful of residues (18 of 337), which are all located in the RBD (residues 319 and 541), participates in more than 50 distinct B cell epitopes (Figure 2A,B).
A comparison between the 337 SARS-CoV-2 spike residues that are recognized by neutralizing antibodies in all selected SARS-CoV-2 variants is shown in Figure 3. We identified that of 337 residues belonging to B cell epitopes, only 37 residues are mutated in at least one of the SARS-CoV-2 variants (Figure 3). However, it is worth noting that many mutations occur in residues that are frequently recognized by neutralizing antibodies. Thus, of the noted residues K417, E484, F486, N487, Y489, Q493 and N501, which are present in more than 75 different B cell epitopes, K417, E484, Q493 and N501 are affected by mutations in at least one variant.
Those mutations that occurred in SARS-CoV-2 residues frequently targeted by neutralizing antibodies suggest that they could have profound disturbing effects. Thereby, we examined the impact of mutations in SARS-CoV-2 VOCs on humoral immunity considering entire B cell epitopes. To that end, we first clustered together all B cell epitopes sharing at least seven or more residues, thus defining B cell epitope cores. We identified 54 clusters including 132 different B cell epitopes (Supplementary Materials Table S1). The number of B cell epitopes in the clusters varied from 1 to 13 and most of them, 104, were mapped on the RBD. There was still significant overlap between the identified B cell epitope cores, but none shared more than six residues. Subsequently, we analyzed how mutations in SARS-CoV-2 VOCs affected B cell epitope cores. The epitope core that defines each cluster along with the analysis of mutations is provided in Table 3.
We found that 43 B cell epitope cores out of 54 were mutated in at least one SARS-CoV-2 variant. The specific number of mutated B cell epitope cores in Alpha, Beta, Gamma, Delta and Omicron variants was 11, 27, 20, 24 and 31, respectively (Table 3). SARS-CoV-2 VOCs do not generally include more than one or two mutations per B cell epitope. A notable exception is the B cell epitope core defined by cluster L18-V1176, which includes 12 residues and all are mutated in the Gamma variant (Table 3).

3.3. Spike-Specific T Cell Epitopes Mutated in SARS-CoV-2 VOCs

SARS-CoV-2 vaccines can also induce spike-specific T cell responses [16], which can contribute to viral clearance and containment [17]. Thereby, we also analyzed the impact of SARS-CoV-2 mutations on experimentally verified spike-specific T cell epitopes that were identified in the IEDB. The selected T cell epitopes included 152 CD8 and 82 CD4 T cell epitopes, respectively. As CD8 T cell epitopes, we only considered those with a size of nine residues, which is the optimal size for binding to class I HLA molecules [18]. In contrast, the size of CD4 T cell epitopes varied between 11 and 18, as HLA II molecules can fit longer peptides than HLA I molecules and with variable lengths [19]. The selected T cell epitopes are restricted by different HLA molecules and distributed throughout the spike protein (Figure 4).
We verified that all the selected T cell epitopes correspond exactly to the reference SARS-CoV-2 spike protein and examined if they were mutated in SARS-CoV-2 VOCs. The result of the analysis for each individual CD8 and CD4 T cell epitope is provided in Supplementary Materials Tables S2 and S3, respectively, and the percentage of conserved epitopes in each variant is shown in Figure 5. Clearly, T cell epitopes are not affected to a great extent by mutations, and the percentage of mutated CD8 or CD4 T cell epitopes is lower than 21% in all SARS-CoV2 VOCs. The variant with less mutated T cell epitopes is Delta (94.1% of CD8 and 93.9% of CD4 T cell epitopes are mutation free, conserved) while Omicron has the largest percentage of mutated CD8 and CD4 epitopes, yet 87.5% and 79.3% of CD8 and CD4 T cell epitopes, respectively, are conserved.
Moreover, attending to the verified HLA I and HLA II restriction elements of CD8 and CD4 T cell epitopes that are not mutated in any of the SARS-CoV-2 VOCs, the percentage of the world population that will have spike-specific CD8 and CD4 T cells elicited by vaccines will be ≥98%.

4. Discussion

The SARS-CoV-2 outbreak has led to the approval of various emergency vaccines to palliate the COVID-19 pandemic crisis and over half of the world population is currently vaccinated. The most widely used COVID-19 vaccines, and the first to obtain clearance for human use by the FDA and EMA agencies, are two mRNA vaccines, Pfizer BNT-162b2 and MODERNA mRNA-1273, two adenovirus vaccines, AstraZeneca AZD-1222 and Johnson & Johnson Ad26.COV2.S and a subunit vaccine, Novavax NVX-CoV2373. These vaccines were trusted to elicit neutralizing antibodies against the spike protein [19] and have shown efficacy and safety in clinical trials [20]. However, numerous SARS-CoV-2 variants, including the Alpha, Beta, Gamma, Delta and Omicron VOCs, have emerged that could defy the immunity generated by COVID-19 vaccines. In fact, SARS-CoV-2 escape mutants do indeed rise in vitro under the presence of neutralizing antibodies [21] and COVID-19 vaccinated subjects are still susceptible to infection or reinfection [22]. The most relevant example is the SARS-CoV-2 Omicron variant, which is more likely to infect vaccinated persons than other variants [23]. However, re-infection of COVID-19 vaccinated subjects can also be the result of waning immunity [24] or inappropriate immunization rather than the result of viral immune escape through variation. Therefore, it is still a matter of debate to what extent SARS-CoV-2 VOCs challenge the immunity elicited by vaccines.
All SARS-CoV-2 VOCs do indeed carry several mutations in the spike protein (Table 1 and Figure 2), including the D614G mutation, which increases receptor affinity and transmissibility [25]. Actually, this same mutation can be found in most SARS-CoV-2 variants [26,27]. SARS-CoV-2 VOCs include additional spike mutations that have been shown to increase infectivity [28,29,30]. For instance, the SARS-CoV-2 Delta variant includes the L452R and E484Q mutations, which, in combination with the P681R mutation (one of the cut-off sites of TMPRSS2 or furin), enhance the dissociation between S1 and S2, and thus viral infectivity [31]. Likewise, the SARS-CoV-2 Omicron variant carries the P681H mutation, which together with the H655Y and N679K mutations, increases viral infectivity by enhancing the excision of spike [32,33]. SARS-CoV-2 Omicron also has the N501Y and Q498R mutations, which further increases the affinity for ACE2 [34]. Overall, SARS-CoV-2 Omicron is the VOC with more mutations in the spike protein (24) followed by Alpha and Gamma (12), Beta (10) and Delta (9) (Figure 1). In this work, we analyzed epitope legacy experimentation to explain why and to what extent these mutations enable SARS-CoV-2 VOCs to escape the immunity elicited by COVID-19 vaccines.
Immune escape by SARS-CoV-2 variants has been analyzed in other works by pure computational means, assessing the impact of mutations on the spike molecular surface of conformational epitopes that are predicted to be recognized by antibodies [35]. In general, these analyses concluded that SARS-CoV-2 variants do not pose a major concern to antibody-mediated immunity elicited by COVID-19 vaccines. However, these previous works did not take into consideration T cell immunity, nor that immune epitope recognition breadth is limited by patterns of immunodominance. Thus, of all the potential epitopes that are susceptible to recognition, the immune system obstinately chose a handful of them [36,37,38,39,40,41]. Thereby, we assessed the chance of immune evasion by SARS-CoV-2 VOCs by analyzing the impact of mutations in experimentally verified B and T cell epitopes, which remain deposited as legacy experimentation in dedicated databases, such as the IEDB. In this approach, we cannot discard potential errors in the data deposited in the IEDB and we have not differentiated between epitopes determined by different experimental assays. Moreover, it is worth noting that we have not considered that some mutations may have a larger impact than others. In particular, some mutations in B cell epitopes, even outside of the identified cores, could alter the protein tertiary structure and thus have a large impact on conformational epitopes.
By considering experimentally verified spike-specific B cell epitopes, we identified 337 SARS-CoV-2 spike residues that are subject to recognition by neutralizing antibodies and of those only 37 residues were mutated in at least one of the SARS-CoV-2 variants. This result alone would appear to indicate that SARS-CoV-2 VOCs are unlikely to escape the action of neutralizing antibodies. However, some of these residues were much more frequently recognized than others, and a clear pattern of immunodominance emerged around the RBD of the spike protein (Figure 2 and Figure 3). This pattern of immunodominance indicates that few mutations in SARS-CoV-2 VOCs can actually have a large detrimental impact on the recognition by neutralizing antibodies elicited by COVID-19 vaccines. We found that all experimentally verified B cell epitopes could be clustered in just 54 clusters sharing at least seven residues (B cell epitope cores). Interestingly, 39 out of 54 B cell epitope cores lie in the RBD region, and 30 are affected by mutations in SARS-CoV-2 VOCs that could impair recognition by antibodies elicited by COVID-19 vaccines (Table 3 and Figure 6). There is already evidence that some of the mutations found in B cell epitope cores, such as N501Y, E484K and K417N can impair antibody neutralization [42,43,44]. Likewise, the L452R and E484Q mutations, found in the Delta variant, have been associated with loss of neutralizing capacity by cognate antibodies [31]. That such few mutations impair neutralizing antibodies against SARS-CoV-2 VOCs emphasizes the identified pattern of immunodominance. The SARS-CoV-2 variant with more B cell cores affected by mutations is Omicron (31 of 54), followed by Beta and Gamma and Delta (27, 20 and 24 mutated cores, respectively) (Table 3 and Figure 6). In addition, the SARS-CoV-2 Omicron variant has, in general, more mutations per B cell epitope core than any other variant, pointing to a greater chance for immune evasion.
Although the protection provided by COVID-19 vaccines against each variant is a case-by-case issue, the number of mutated B cell epitope cores in SARS-CoV-2 VOCs explains the current data on the effectiveness of COVID-19 vaccines, particularly in terms of neutralizing activity. For instance, in clinical assays, AstraZeneca COVID-19 vaccine has shown an effectiveness of <70% [45] against the SARS-CoV-2 Alpha variant, but only 10–22% against the SARS-CoV-2 Beta variant [46,47], which doubles SARS-CoV-2 Alpha variant in mutated B cell epitope cores. It has also been shown that antibodies elicited by Moderna and Novavax vaccines can neutralize the SARS-CoV-2 Alpha variant [48] but the reported effectiveness against the Beta variant is 49% [49] Likewise, serum from people vaccinated with BNT162b2 (Pfizer vaccine) cannot effectively neutralize Beta and Gamma variants [43]. COVID-19 vaccines’ effectiveness against SARS-CoV-2 Gamma and Delta variants is also similar since, in both variants, mRNA-1273 and BNT162b2 suffer a partial loss of neutralization capacity [50]. Overall, Beta, Gamma and Delta have shown higher resistance to neutralization compared to Alpha [51]. In line with our results, many studies agree that SARS-CoV-2 Omicron is largely resistant to current COVID-19 vaccines [50] and makes it understandable that SARS-CoV-2 Omicron became the dominant variant during mass vaccinations [52]. In line with the epidemiological data, sera from vaccinated individuals appear to contain very low to undetectable levels of neutralizing antibodies against SARS-CoV-2 Omicron [53]. Interestingly, SARS-CoV-2 Omicron does not completely escape COVID-19 vaccines [54], nor any other variants [55]. Moreover, it has been reported that a third booster dose with BNT162b2, AZD-1222, or mRNA-1273 vaccines enhances protection against the SARS-CoV-2 Omicron variant [56,57]. While it cannot be discarded that a third booster dose could lead to neutralizing antibodies, it is likely that the enhanced protection is mediated by T cells.
It is now accepted that mRNA and adenovirus COVID-19 vaccines induce T cell responses and memory [16,58,59], which plays a pivotal role in anti-SARS-CoV-2 immunity [17,60]. Thereby, we also analyzed the impact of mutations in 152 CD8 and 82 CD4 T cell epitopes known to be targets of T cell responses induced by SARS-CoV-2 infection or vaccination. Unlike B cell epitopes, these T cell epitopes are distributed throughout the entire spike protein, without displaying a clear immunodominant region. This is expected as T cell epitopes must be recognized bound to HLA molecules, which are highly polymorphic in the population and allelic variants bind distinct sets of peptides [61]. Compared with B cells, T cell epitopes were affected by mutations to a much lesser extent, and over 79% and 87% of CD4 and CD8 T cells were preserved in each distinct SARS-CoV-2 VOC (Figure 5). Moreover, the diversity of the restriction elements of conserved T cell epitopes guarantees that spike-specific T cell responses are present in the COVID-19 vaccinated population against any SARS-CoV-2 VOC, regardless of their HLA background. These results support that the current evidence indicating that COVID-19 vaccines induce T cell memory against all five SARS-CoV-2 VOCs [62] applies to the entire world population.
The detected preservation of cellular immune responses could also explain why despite the high rate of infection by SARS-CoV-2 Omicron, severity, deaths and hospitalizations have decreased upon the introduction of COVID-19 vaccines [63]. Moreover, given the resilience of T cell immunity to SARS-CoV-2 variants, it would be of great interest to explore new vaccine strategies that favor this type of immunity, including that of eliciting cross-reactive T cells by tetanic vaccines [64,65,66]. Peptide-based vaccines, such as CoVac-1, have also been found to be a valid strategy to induce T cell immunity against SARS-CoV-2 [67]. Moreover, the induction of T cell immunity against conserved epitope peptide fragments shared between SARS-CoV-2 and other human coronaviruses could lead to a universal vaccine against emerging coronavirus. These vaccines could be instrumental in the generation of long-term immunity mediated by memory T cells, overcoming the fact that antibody levels generated by coronavirus infections and COVID-19 vaccines are short lived [58]. Neutralizing antibodies following natural infection with SARS-CoV-2 are diminished after a year while T cell response is still detectable even against variants [68].

5. Conclusions

The pattern of mutations on spike-specific B cell epitopes targeted by neutralizing antibodies supports that SARS-CoV-2 VOCs, particularly SARS-CoV-2 Omicron, can escape humoral adaptive immunity induced by COVID-19 vaccines. In contrast, SARS-CoV-2 VOCs cannot escape T cell immunity mediated by COVID-19 vaccines, as most spike-specific T cell epitopes are conserved. These findings highlight the relevance of T cell immunity induced by COVID-19 vaccines.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines10050996/s1, Supplementary Table S1: B cell epitope clusters in SARS-CoV-2 spike glycoprotein; Supplementary Table S2: Conservation of SARS-CoV-2 spike CD8 T cell epitopes in variants of concern; Supplementary Table S3: Conservation of SARS-CoV-2 spike CD4 T cell epitopes in variants of concern.

Author Contributions

Conceptualization, P.A.R.; Methodology, L.B.-S., A.R.-C. and P.A.R.; Investigation: L.B.-S., H.F.P.-P. and P.A.R.; Writing—Original Draft: L.B.-S., H.F.P.-P. and P.A.R.; Final Writing and Editing: L.B.-S., H.F.P.-P., A.R.-C. and P.A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by a REACT-European Union grant from the Comunidad de Madrid to the ANTICIPA project of Complutense University of Madrid. The funders had no role in study design, data collection and analysis, preparation of the manuscript or decision to publish.Esta actividad está financiada por la Comunidad de Madrid y la Unión Europea, a través del Fondo Europeo de Desarrollo Regional (FEDER), financiado como parte de la respuesta de la Unión a la pandemia de COVID-19.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data is available online as supplementary material at https://www.mdpi.com/article/10.3390/biomedicines10050996/s1.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Marian, A.J. Current state of vaccine development and targeted therapies for COVID-19: Impact of basic science discoveries. Cardiovasc. Pathol. 2021, 50, 107278. [Google Scholar] [CrossRef] [PubMed]
  2. Schoch, C.L.; Ciufo, S.; Domrachev, M.; Hotton, C.L.; Kannan, S.; Khovanskaya, R.; Leipe, D.; Mcveigh, R.; O’Neill, K.; Robbertse, B.; et al. NCBI Taxonomy: A comprehensive update on curation, resources and tools. Database 2020, 2020, baaa062. [Google Scholar] [CrossRef]
  3. Shamsi, A.; Mohammad, T.; Anwar, S.; AlAjmi, M.F.; Hussain, A.; Rehman, M.T.; Islam, A.; Hassan, M.I. Glecaprevir and Maraviroc are high-affinity inhibitors of SARS-CoV-2 main protease: Possible implication in COVID-19 therapy. Biosci. Rep. 2020, 26, 40. [Google Scholar] [CrossRef]
  4. Hwang, Y.C.; Lu, R.M.; Su, S.C.; Chiang, P.Y.; Ko, S.H.; Ke, F.Y.; Liang, K.H.; Hsieh, T.Y.; Wu, H.C. Monoclonal antibodies for COVID-19 therapy and SARS-CoV-2 detection. J. Biomed. Sci. 2022, 29, 1. [Google Scholar] [CrossRef]
  5. Walls, A.C.; Park, Y.J.; Tortorici, M.A.; Wall, A.; McGuire, A.T.; Veesler, D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell 2020, 181, 281–292.e6. [Google Scholar] [CrossRef]
  6. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280.e8. [Google Scholar] [CrossRef] [PubMed]
  7. Ju, B.; Zhang, Q.; Ge, J.; Wang, R.; Sun, J.; Ge, X.; Yu, J.; Shan, S.; Zhou, B.; Song, S.; et al. Human neutralizing antibodies elicited by SARS-CoV-2 infection. Nature 2020, 584, 115–119. [Google Scholar] [CrossRef] [PubMed]
  8. Weisblum, Y.; Schmidt, F.; Zhang, F.; DaSilva, J.; Poston, D.; Lorenzi, J.C.C.; Muecksch, F.; Rutkowska, M.; Hoffmann, H.H.; Michailidis, E.; et al. Escape from neutralizing antibodies 1 by SARS-CoV-2 spike protein variants. Elife 2020, 9, e61312. [Google Scholar] [CrossRef]
  9. Galloway, S.E.; Paul, P.; MacCannell, D.R.; Johansson, M.A.; Brooks, J.T.; MacNeil, A.; Slayton, R.B.; Tong, S.; Silk, B.J.; Armstrong, G.L.; et al. Emergence of SARS-CoV-2 B.1.1.7 Lineage—United States, December 29, 2020–January 12, 2021. MMWR Morb. Mortal Wkly. Rep. 2021, 70, 95–99. [Google Scholar] [CrossRef]
  10. Choi, J.Y.; Smith, D.M. SARS-CoV-2 variants of concern. Yonsei Med. J. 2021, 62, 961–968. [Google Scholar] [CrossRef]
  11. Vita, R.; Overton, J.A.; Greenbaum, J.A.; Ponomarenko, J.; Clark, J.D.; Cantrell, J.R.; Wheeler, D.K.; Gabbard, J.L.; Hix, D.; Sette, A.; et al. The immune epitope database (IEDB) 3.0. Nucleic Acids Res. 2015, 43, D405–D412. [Google Scholar] [CrossRef]
  12. Vita, R.; Mahajan, S.; Overton, J.A.; Dhanda, S.K.; Martini, S.; Cantrell, J.R.; Wheeler, D.K.; Sette, A.; Peters, B. The Immune Epitope Database (IEDB): 2018 update. Nucleic Acids Res. 2019, 47, D339–D343. [Google Scholar] [CrossRef] [Green Version]
  13. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; Mcgettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [Green Version]
  14. Pérez-Silva, J.G.; Araujo-Voces, M.; Quesada, V. NVenn: Generalized, quasi-proportional Venn and Euler diagrams. Bioinform. Oxf. Univ. Press 2018, 34, 2322–2324. [Google Scholar] [CrossRef] [Green Version]
  15. Molero-Abraham, M.; Lafuente, E.M.; Flower, D.R.; Reche, P.A. Selection of conserved epitopes from hepatitis c virus for pan-populational stimulation of T-cell responses. Clin. Dev. Immunol. 2013, 2013, 601943. [Google Scholar] [CrossRef]
  16. Nguyen, T.H.O.; Cohen, C.A.; Rowntree, L.C.; Bull, M.B.; Hachim, A.; Kedzierska, K.; Valkenburg, S.A. T Cells Targeting SARS-CoV-2: By Infection, Vaccination, and Against Future Variants. Front. Med. 2021, 8, 793102. [Google Scholar] [CrossRef]
  17. Niessl, J.; Sekine, T.; Buggert, M. T cell immunity to SARS-CoV-2. Semin. Immunol. 2021, 55, 101505. [Google Scholar] [CrossRef] [PubMed]
  18. Lafuente, E.; Reche, P. Prediction of MHC-Peptide Binding: A Systematic and Comprehensive Overview. Curr. Pharm. Des. 2009, 15, 3209–3220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Kyriakidis, N.C.; López-Cortés, A.; González, E.V.; Grimaldos, A.B.; Prado, E.O. SARS-CoV-2 vaccines strategies: A comprehensive review of phase 3 candidates. npj Vaccines 2021, 6, 28. [Google Scholar] [CrossRef]
  20. Cheng, H.; Peng, Z.; Luo, W.; Si, S.; Mo, M.; Zhou, H.; Xin, X.; Liu, H.; Yu, Y. Efficacy and safety of covid-19 vaccines in phase iii trials: A meta-analysis. Vaccines 2021, 9, 582. [Google Scholar] [CrossRef] [PubMed]
  21. Harvey, W.T.; Carabelli, A.M.; Jackson, B.; Gupta, R.K.; Thomson, E.C.; Harrison, E.M.; Ludden, C.; Reeve, R.; Rambaut, A.; COVID-19 Genomics UK (COG-UK) Consortium; et al. SARS-CoV-2 variants, spike mutations and immune escape. Nat. Rev. Microbiol. 2021, 19, 409–424. [Google Scholar] [CrossRef] [PubMed]
  22. Hacisuleyman, E.; Hale, C.; Saito, Y.; Blachere, N.E.; Bergh, M.; Conlon, E.G.; Schaefer-Babajew, D.J.; DaSilva, J.; Muecksch, F.; Gaebler, C.; et al. Vaccine Breakthrough Infections with SARS-CoV-2 Variants. N. Engl. J. Med. 2021, 384, 2212–2218. [Google Scholar] [CrossRef]
  23. Chaguza, C.; Coppi, A.; Earnest, R.; Ferguson, D.; Kerantzas, N.; Warner, F.; Young, H.P.; Breban, M.I.; Billig, K.; Koch, R.T.; et al. Rapid Emergence of SARS-CoV-2 Omicron Variant Is Associated with an Infection Advantage over Delta in Vaccinated Persons. medRxiv, 2022; online ahead of print. [Google Scholar] [CrossRef]
  24. Levin, E.G.; Lustig, Y.; Cohen, C.; Fluss, R.; Indenbaum, V.; Amit, S.; Doolman, R.; Asraf, K.; Mendelson, E.; Ziv, A.; et al. Waning Immune Humoral Response to BNT162b2 Covid-19 Vaccine over 6 Months. N. Engl. J. Med. 2021, 385, e84. [Google Scholar] [CrossRef]
  25. Volz, E.; Hill, V.; McCrone, J.T.; Price, A.; Jorgensen, D.; O’Toole, Á.; Southgate, J.; Johnson, R.; Jackson, B.; Nascimento, F.; et al. Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity. Cell 2021, 184, 64–75.e11. [Google Scholar] [CrossRef] [PubMed]
  26. Hou, Y.J.; Chiba, S.; Halfmann, P.; Ehre, C.; Kuroda, M.; Dinnon, K.H.; Leist, S.R.; Schäfer, A.; Nakajima, N.; Takahashi, K.; et al. SARS-CoV-2 D614G variant exhibits efficient replication ex vivo and transmission in vivo. Science 2020, 370, 1464–1468. [Google Scholar] [CrossRef]
  27. Korber, B.; Fischer, W.M.; Gnanakaran, S.; Yoon, H.; Theiler, J.; Abfalterer, W.; Hengartner, N.; Giorgi, E.E.; Bhattacharya, T.; Foley, B.; et al. Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus. Cell 2020, 182, 812–827.e19. [Google Scholar] [CrossRef]
  28. Starr, T.N.; Greaney, A.J.; Hilton, S.K.; Ellis, D.; Crawford, K.H.D.; Dingens, A.S.; Navarro, M.; Bowen, J.; Tortorici, M.; Walls, A.; et al. Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding. Cell 2020, 182, 1295–1310.e20. [Google Scholar] [CrossRef] [PubMed]
  29. Teruel, N.; Mailhot, O.; Najmanovich, R.J. Modelling conformational state dynamics and its role on infection for SARS-CoV-2 Spike protein variants. PLoS Comput. Biol. 2021, 17, e1009286. [Google Scholar] [CrossRef]
  30. Nelson, G.; Buzko, O.; Patricia, S.; Niazi, K.; Rabizadeh, S.; Soon-Shiong, P. Molecular dynamic simulation reveals E484K mutation enhances spike RBD-ACE2 affinity and the 1 combination of E484K, K417N and N501Y mutations (501Y.V2 variant) induces conformational 2 change greater than N501Y mutant alone, potentially resulting in an escape mutant. bioRxiv 2021. [Google Scholar] [CrossRef]
  31. Cherian, S.; Potdar, V.; Jadhav, S.; Yadav, P.; Gupta, N.; Das, M.; Rakshit, P.; Singh, S.; Abraham, P.; Panda, S.; et al. SARS-CoV-2 Spike Mutations, L452R, T478K, E484Q and P681R, in the Second Wave of COVID-19 in Maharashtra, India. Microorganisms 2021, 9, 1542. [Google Scholar] [CrossRef] [PubMed]
  32. Hossain, M.G.; Tang, Y.; Akter, S.; Zheng, C. Roles of the polybasic furin cleavage site of spike protein in SARS-CoV-2 replication, pathogenesis, and host immune responses and vaccination. J. Med. Virol. 2021, 94, 1815–1820. [Google Scholar] [CrossRef]
  33. Gong, S.Y.; Chatterjee, D.; Richard, J.; Prévost, J.; Tauzin, A.; Gasser, R.; Bo, Y.; Vézina, D.; Goyette, G.; Gendron-Lepage, G.; et al. Contribution of single mutations to selected SARS-CoV-2 emerging variants spike antigenicity. Virology 2021, 563, 134–145. [Google Scholar] [CrossRef] [PubMed]
  34. Zahradník, J.; Marciano, S.; Shemesh, M.; Zoler, E.; Harari, D.; Chiaravalli, J.; Meyer, B.; Rudich, Y.; Li, C.; Marton, I.; et al. SARS-CoV-2 variant prediction and antiviral drug design are enabled by RBD in vitro evolution. Nat. Microbiol. 2021, 6, 1188–1198. [Google Scholar] [CrossRef]
  35. Barnes, C.O.; Jette, C.A.; Abernathy, M.E.; Dam, K.M.A.; Esswein, S.R.; Gristick, H.B.; Malyutin, A.; Sharaf, N.; Huey-Tubman, K.; Lee, Y.; et al. SARS-CoV-2 neutralizing antibody structures inform therapeutic strategies. Nature 2020, 588, 682–687. [Google Scholar] [CrossRef] [PubMed]
  36. Mathew, N.R.; Angeletti, D. Recombinant Influenza Vaccines: Saviors to Overcome Immunodominance. Front. Immunol. 2020, 10, 2997. [Google Scholar] [CrossRef] [PubMed]
  37. Abbott, R.K.; Crotty, S. Factors in B cell competition and immunodominance. Immunol. Rev. 2020, 296, 120–131. [Google Scholar] [CrossRef]
  38. Angeletti, D.; Yewdell, J.W. Understanding and Manipulating Viral Immunity: Antibody Immunodominance Enters Center Stage. Trends Immunol. 2018, 39, 549–561. [Google Scholar] [CrossRef] [PubMed]
  39. Lehmann, A.A.; Kirchenbaum, G.A.; Zhang, T.; Reche, P.A.; Lehmann, P.V. Deconvoluting the T Cell Response to SARS-CoV-2: Specificity Versus Chance and Cognate Cross-Reactivity. Front. Immunol. 2021, 12, 635942. [Google Scholar] [CrossRef]
  40. Kim, A.R.; Sadegh-Nasseri, S. Determinants of immunodominance for CD4 T cells. Curr. Opin. Immunol. 2015, 34, 9–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Yewdell, J.W. Confronting Complexity: Real-World Immunodominance in Antiviral CD8+ T Cell Responses. Immunity 2006, 25, 533–543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Nonaka, C.K.V.; Franco, M.M.; Gräf, T.; De Lorenzo Barcia, C.A.; De Ávila Mendonça, R.N.; De Sousa, K.A.F.; Neiva, L.; Fosenca, V.; Mendes, A.; de Aguiar, R.; et al. Genomic evidence of SARS-CoV-2 reinfection involving E484K spike mutation, Brazil. Emerg. Infect. Dis. 2021, 27, 1522–1524. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, P.; Nair, M.S.; Liu, L.; Iketani, S.; Luo, Y.; Guo, Y.; Wang, M.; Yu, J.; Zhang, B.; Kwong, P.D.; et al. Antibody resistance of SARS-CoV-2 variants B.1.351 and B.1.1.7. Nature 2021, 593, 130–135. [Google Scholar] [CrossRef] [PubMed]
  44. Hoffmann, M.; Arora, P.; Groß, R.; Seidel, A.; Hörnich, B.F.; Hahn, A.S.; Krüger, N.; Graichen, L.; Hofmann-Winkler, H.; Kempf, A.; et al. SARS-CoV-2 variants B.1.351 and P.1 escape from neutralizing antibodies. Cell 2021, 184, 2384–2393.e12. [Google Scholar] [CrossRef] [PubMed]
  45. Voysey, M.; Clemens, S.A.C.; Madhi, S.A.; Weckx, L.Y.; Folegatti, P.M.; Aley, P.K.; Angus, B.; Baillie, V.L.; Barnabas, S.L.; Bhorat, Q.E.; et al. Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: An interim analysis of four randomised controlled trials in Brazil, South Africa, and the UK. Lancet 2021, 397, 99–111. [Google Scholar] [CrossRef]
  46. Madhi, S.A.; Izu, A.; Pollard, A.J. ChAdOx1 nCoV-19 Vaccine Efficacy against the B.1.351 Variant. Reply. N. Engl. J. Med. 2021, 385, 571–572. [Google Scholar]
  47. Madhi, S.A.; Baillie, V.; Cutland, C.L.; Voysey, M.; Koen, A.L.; Fairlie, L.; Padayachee, S.D.; Dheda, K.; Barnabas, S.L.; Bhorat, Q.E.; et al. Efficacy of the ChAdOx1 nCoV-19 Covid-19 Vaccine against the B.1.351 Variant. N. Engl. J. Med. 2021, 384, 1885–1898. [Google Scholar] [CrossRef] [PubMed]
  48. Shen, X.; Tang, H.; McDanal, C.; Wagh, K.; Fischer, W.; Theiler, J.; Yoon, H.; Li, D.; Haynes, B.; Saunders, K.; et al. SARS-CoV-2 variant B.1.1.7 is susceptible to neutralizing antibodies elicited by ancestral spike vaccines. Cell Host Microbe 2021, 29, 529–539.e3. [Google Scholar] [CrossRef]
  49. Wadman, M. Novavax vaccine delivers 89% efficacy against COVID-19 in U.K.—But is less potent in South Africa. Science 2021, 12, 277. [Google Scholar] [CrossRef]
  50. Rayati Damavandi, A.; Dowran, R.; Al Sharif, S.; Kashanchi, F.; Jafari, R. Molecular variants of SARS-CoV-2: Antigenic properties and current vaccine efficacy. Med. Microbiol. Immunol. 2022. [Google Scholar] [CrossRef]
  51. Mistry, P.; Barmania, F.; Mellet, J.; Peta, K.; Strydom, A.; Viljoen, I.M.; James, W.; Gordon, S.; Pepper, M. SARS-CoV-2 Variants, Vaccines, and Host Immunity. Front. Immunol. 2022, 12, 809244. [Google Scholar] [CrossRef] [PubMed]
  52. World Health Organization. Enhancing response to Omicron SARS-CoV-2 Variant. Available online: https://www.who.int/publications/m/item/enhancing-readiness-for-omicron-(b.1.1.529)-technical-brief-and-priority-actions-for-member-states (accessed on 22 April 2022).
  53. Carreño, J.M.; Alshammary, H.; Tcheou, J.; Singh, G.; Raskin, A.; Kawabata, H.; Sominsky, L.A.; Clark, J.J.; Adelsberg, D.C.; Bielak, D.A.; et al. Activity of convalescent and vaccine serum against SARS-CoV-2 Omicron. Nature 2021, 602, 682–688. [Google Scholar] [CrossRef]
  54. Cele, S.; Jackson, L.; Khoury, D.S.; Khan, K.; Moyo-Gwete, T.; Tegally, H.; San, J.; Cromer, D.; Scheepers, C.; Amoako, D.; et al. Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization. Nature 2021, 602, 654–656. [Google Scholar] [CrossRef]
  55. Fiolet, T.; Kherabi, Y.; Macdonald, C.; Ghosn, J.; Peiffer-smadja, N. Comparing COVID-19 vaccines for their characteristics, efficacy and effectiveness against SARS-CoV-2 and variants of concern: A narrative review. Clin. Microbiol. Infect. 2022, 28, 202–221. [Google Scholar] [CrossRef]
  56. Dejnirattisai, W.; Huo, J.; Zhou, D.; Zahradník, J.; Supasa, P.; Liu, C.; Duyvesteyn, H.; Ginn, H.M.; Mentzer, A.J.; Tuekprakhon, A.; et al. SARS-CoV-2 Omicron-B.1.1.529 leads to widespread escape from neutralizing antibody responses. Cell 2022, 185, 467–484.e15. [Google Scholar] [CrossRef] [PubMed]
  57. Andrews, N.; Stowe, J.; Kirsebom, F.; Toffa, S.; Rickeard, T.; Gallagher, E.; Gower, C.; Kall, M.; Groves, N.; O’Connell, A.M.; et al. Covid-19 Vaccine Effectiveness against the Omicron (B.1.1.529) Variant. N. Engl. J. Med. 2022, 386, 1532–1546. [Google Scholar] [CrossRef] [PubMed]
  58. Ewer, K.J.; Barrett, J.R.; Belij-Rammerstorfer, S.; Sharpe, H.; Makinson, R.; Morter, R.; Flaxman, A.; Wright, D.; Bellamy, D.; Bittaye, M.; et al. T cell and antibody responses induced by a single dose of ChAdOx1 nCoV-19 (AZD1222) vaccine in a phase 1/2 clinical trial. Nat. Med. 2021, 27, 270–278. [Google Scholar] [CrossRef]
  59. Sahin, U.; Muik, A.; Derhovanessian, E.; Vogler, I.; Kranz, L.M.; Vormehr, M.; Baum, A.; Pascal, K.; Quandt, J.; Maurus, D.; et al. COVID-19 vaccine BNT162b1 elicits human antibody and TH1 T cell responses. Nature 2020, 586, 594–599. [Google Scholar] [CrossRef]
  60. Moss, P. The T cell immune response against SARS-CoV-2. Vol. 23, Nature Immunology. Nat. Res. 2022, 23, 186–193. [Google Scholar]
  61. Reche, P.A.; Reinherz, E.L. Sequence variability analysis of human class I and class II MHC molecules: Functional and structural correlates of amino acid polymorphisms. J. Mol. Biol. 2003, 331, 623–641. [Google Scholar] [CrossRef] [Green Version]
  62. Tarke, A.; Coelho, C.H.; Zhang, Z.; Dan, J.M.; Yu, E.D.; Methot, N.; Bloom, N.I.; Goodwin, B.; Phillips, E.; Mallal, S.; et al. SARS-CoV-2 vaccination induces immunological T cell memory able to cross-recognize variants from Alpha to Omicron. Cell 2022, 185, 847–859.e11. [Google Scholar] [CrossRef] [PubMed]
  63. Barandalla, I.; Alvarez, C.; Barreiro, P.; de Mendoza, C.; González-Crespo, R.; Soriano, V. Impact of scaling up SARS-CoV-2 vaccination on COVID-19 hospitalizations in Spain. Int. J. Infect. Dis. 2021, 112, 81–88. [Google Scholar] [CrossRef] [PubMed]
  64. Monereo-Sánchez, J.; Luykx, J.J.; Pinzón-Espinosa, J.; Richard, G.; Motazedi, E.; Westlye, L.T.; Andreassen, O.A.; van der Meer, D. Diphtheria and Tetanus Vaccination History Is Associated With Lower Odds of COVID-19 Hospitalization. Front. Immunol. 2021, 12, 749264. [Google Scholar] [CrossRef]
  65. Mysore, V. Clinical and Translational Article Protective heterologous T cell immunity in COVID-19 induced by the trivalent MMR and Tdap vaccine antigens. Med 2021, 2, 1050–1071.e7. [Google Scholar] [CrossRef]
  66. Reche, P.A. Potential Cross-Reactive Immunity to SARS-CoV-2 From Common Human Pathogens and Vaccines. Front. Immunol. 2020, 11, 586984. [Google Scholar] [CrossRef] [PubMed]
  67. Heitmann, J.S.; Bilich, T.; Tandler, C.; Nelde, A.; Maringer, Y.; Marconato, M.; Reusch, J.; Jäger, S.; Denk, M.; Richter, M.; et al. A COVID-19 peptide vaccine for the induction of SARS-CoV-2 T cell immunity. Nature 2022, 601, 617–622. [Google Scholar] [CrossRef]
  68. Guo, L.; Wang, G.; Wang, Y.; Zhang, Q.; Ren, L.; Gu, X.; Huang, T.; Zhong, J.; Wang, Y.; Wang, X.; et al. SARS-CoV-2-specific antibody and T-cell responses 1 year after infection in people recovered from COVID-19: A longitudinal cohort study. Lancet Microbe, 2022; online ahead of print. [Google Scholar] [CrossRef]
Biomedicines 10 00996 g001 550
Figure 1. Mutations in spike glycoprotein. (A) Total number of mutations (Substitutions in pink and deletions in green) of each spike glycoprotein from SARS-CoV-2 variants; (B) Venn diagram showing the shared mutations (by residue position) between variants.
Figure 1. Mutations in spike glycoprotein. (A) Total number of mutations (Substitutions in pink and deletions in green) of each spike glycoprotein from SARS-CoV-2 variants; (B) Venn diagram showing the shared mutations (by residue position) between variants.
Biomedicines 10 00996 g001
Biomedicines 10 00996 g002 550
Figure 2. Spike residues targeted by neutralizing antibodies. (A) Plot depicting the occurrence/frequency of SARS-CoV-2 spike residues (X-axis) in distinct B cell epitopes targeted by neutralizing antibodies (Y-axis). The receptor-binding domain (RBD) is highlighted in black; (B) Molecular surface of SARS-CoV-2 spike protein with residue recognition frequency by neutralizing antibodies shown in a color gradient from dark blue (not recognized) to green and then to bright red (most frequently recognized). Rendering was generated with PyMOL upon a PDB of the 3D-structure of SARS-CoV-2 spike protein with the occurrence of residues in B cell epitopes assigned to B-factors.
Figure 2. Spike residues targeted by neutralizing antibodies. (A) Plot depicting the occurrence/frequency of SARS-CoV-2 spike residues (X-axis) in distinct B cell epitopes targeted by neutralizing antibodies (Y-axis). The receptor-binding domain (RBD) is highlighted in black; (B) Molecular surface of SARS-CoV-2 spike protein with residue recognition frequency by neutralizing antibodies shown in a color gradient from dark blue (not recognized) to green and then to bright red (most frequently recognized). Rendering was generated with PyMOL upon a PDB of the 3D-structure of SARS-CoV-2 spike protein with the occurrence of residues in B cell epitopes assigned to B-factors.
Biomedicines 10 00996 g002
Biomedicines 10 00996 g003 550
Figure 3. Spike glycoprotein residues recognized by neutralizing antibodies. A multiple sequence alignment of the spike glycoprotein from SARS-CoV-2 and variants of concern (VOCs) was generated and residues reported to be part of B cell epitopes recognized by neutralizing antibodies are shown. Spike residues that are mutated in SARS-CoV-2 VOCs are boxed and mutations highlighted. The receptor-binding domain (RBD) is shaded in cyan. Residue numbering is shown at the top. The bottom bar represents the occurrence of residues in B cell epitopes (1 to 99) in a color gradient from dark blue (1, at least in one B cell epitope) to green and then to bright red (top value is for E484 which is part of 99 B cell epitopes).
Figure 3. Spike glycoprotein residues recognized by neutralizing antibodies. A multiple sequence alignment of the spike glycoprotein from SARS-CoV-2 and variants of concern (VOCs) was generated and residues reported to be part of B cell epitopes recognized by neutralizing antibodies are shown. Spike residues that are mutated in SARS-CoV-2 VOCs are boxed and mutations highlighted. The receptor-binding domain (RBD) is shaded in cyan. Residue numbering is shown at the top. The bottom bar represents the occurrence of residues in B cell epitopes (1 to 99) in a color gradient from dark blue (1, at least in one B cell epitope) to green and then to bright red (top value is for E484 which is part of 99 B cell epitopes).
Biomedicines 10 00996 g003
Biomedicines 10 00996 g004 550
Figure 4. Distribution of spike-specific T cell epitopes. The Figure depicts the location of 152 CD8 (blue bars) and 82 CD4 (orange) T cell epitopes in the primary structure of SARS-CoV-2 spike protein. CD8 and CD4 T cell epitopes were obtained from the Immune Epitope Database resource.
Figure 4. Distribution of spike-specific T cell epitopes. The Figure depicts the location of 152 CD8 (blue bars) and 82 CD4 (orange) T cell epitopes in the primary structure of SARS-CoV-2 spike protein. CD8 and CD4 T cell epitopes were obtained from the Immune Epitope Database resource.
Biomedicines 10 00996 g004
Biomedicines 10 00996 g005 550
Figure 5. CD8 and CD4 T cell epitopes conservation in SARS-CoV-2 variants of concern (VOCs). Percentage of T cell epitopes (CD8 in blue and CD4 in orange) not affected by mutations in SARS-CoV-2 VOCs.
Figure 5. CD8 and CD4 T cell epitopes conservation in SARS-CoV-2 variants of concern (VOCs). Percentage of T cell epitopes (CD8 in blue and CD4 in orange) not affected by mutations in SARS-CoV-2 VOCs.
Biomedicines 10 00996 g005
Biomedicines 10 00996 g006 550
Figure 6. B cell epitope cores affected by mutations in SARS-CoV-2 variants. All 132 B cell epitopes targeted by neutralizing antibodies could be grouped in 54 clusters sharing seven or more residues (B cell epitope cores), and those mutated in each variant are indicated. The length of the bars is proportional to the percentage of mutated residues. Clusters within the receptor-binding domain (RBD) are framed. The residue composition of B cell epitope cores is provided in Table 3.
Figure 6. B cell epitope cores affected by mutations in SARS-CoV-2 variants. All 132 B cell epitopes targeted by neutralizing antibodies could be grouped in 54 clusters sharing seven or more residues (B cell epitope cores), and those mutated in each variant are indicated. The length of the bars is proportional to the percentage of mutated residues. Clusters within the receptor-binding domain (RBD) are framed. The residue composition of B cell epitope cores is provided in Table 3.
Biomedicines 10 00996 g006
Table
Table 1. SARS-CoV-2 sample genome sequences selected as representative lineage variants.
Table 1. SARS-CoV-2 sample genome sequences selected as representative lineage variants.
Virus NameGISAID Accession IDLineageLocationCollection Date
hCoV-19/England/
QEUH-119AACD/2021		EPI_ISL_944882B.1.1.7 (Alpha)United Kingdom/England28 January 2021
hCoV-19/South Africa/N00390/2020EPI_ISL_712081B.1.315 (Beta)South Africa/Eastern Cape8 October 2020
hCoV-19/Brazil/AM-1031/2021EPI_ISL_906075P.1 (Gamma)Brazil/Manaus19 January 2021
hCoV-19/India/TN-SEQ_5119_S219_R1_001/2021EPI_ISL_6033571B.1.617.2 (Delta)India/Tamil Nadu26 July 2021
hCoV-19/Nepal/NPHL-S-263/2021EPI_ISL_7196120B.1.1.529 (Omicron)Nepal/Bagmati26 November 2021
Table
Table 2. Amino acid substitutions and deletions in spike glycoprotein from SARS-CoV-2 variants.
Table 2. Amino acid substitutions and deletions in spike glycoprotein from SARS-CoV-2 variants.
VariantAmino Acid SubstitutionsAmino Acid Deletions
B.1.1.7 (Alpha)S13T, T385N, N501Y, A570D, D614G, P681H, T716I, S982A, D1118HH69-, V70-, Y144-
B.1.315 (Beta)D80A, R246I, K417N, E484K, N501Y, D614G, A701VL241-, L242-, A243-
P.1 (Gamma)L18F, T20N, P26S, D138Y, R190S, K417T, E484K, N501Y, D614G, H655Y, T1027I, V1176F-
B.1.617.2 (Delta)T19R, G142D, L452R, T478K, E484Q, D614G, P681R, P812T, D950N-
B. 1.1.529
(Omicron)		A67V, T95I, G339D, S371L, S373P, S375F, S477N, T478K, E484A, Q493R, G496S, Q498R, N501Y, Y505H, T547K, D614G, H655Y, N679K, P681H, D796Y, N856K, Q954H, N969K, L981F-
-, No amino acid deletions were found in Gamma, Delta and Omicron.
Table
Table 3. Mutations in spike-specific B cell epitope clusters in SARS-CoV-2 variants.
Table 3. Mutations in spike-specific B cell epitope clusters in SARS-CoV-2 variants.
Cluster 1Epitope CoreN° Mutations
αβγδο
S12-R246S12,C15,L18,T19,C136,G142,H146,K147,N149,R24601120
Q14-D253Q14,Y144,Y145,H146,K147,F157,G252,D25310000
V16-D253V16,Y144,K147,R246,P251,G252,D25301000
L18-T716L18,Y144,L242,A243,L244,R246,E484,N501,A701,T71626312
L18-V1176L18,T20,P26,D138,R190,K417,E484,N501,D614,H655,T1027,V1176241224
H69-G261H69,V70,Y144,L242,A243,L244,G26122000
D80-A701D80,D215,K417,E484,N501,D614,A70126423
F140-H245F140,G142,V143,Y145,H146,N148,N149,W152,E154,F157,A243,L244,H24511010
Y144-Y248Y144,Y145,H146,K147,K150,W152,Y24810000
Y144-P251Y144,H146,K147,R246,Y248,L249,P25101000
Y144-D253Y144,Y145,H146,K147,R158,R246,L249,T250,P251,G252,D25311000
Y144-S255Y144,R246,S247,Y248,T250,P251,D253,S254,S25501000
Y145-G261Y145,H146,K147,R246,Y248,P251,G26111000
N334-L441N334,L335,P337,G339,E340,N343,A344,T345,K356,R357,S359,C361,L44100001
G339-L441G339,E340,N343,A344,T345,R346,L44100001
T342-F487T342,R343,Y348,Y446,N447,L449,I465,T467,F48700010
R343-T497R343,N436,N437,L438,K441,V442,G443,G444,P496,T49700000
T345-V445T345,R346,S373,W436,N437,N440,L441,K444,V44500001
T345-V455-N450T345,R346,N440,L441,S443,K444,V445,N45000000
T345-Y449-N450T345,R346,K444,G446,G447,Y449,N45000000
T345-S494T345,R346,A348,S349,N354,K356,Y449,N450,L452,T470,Q493,S49400011
R346-Q498R346,N440,L441,K444,V445,G446,N448,Y449,Q49800001
Y351-R509Y351,A352,W353,N360,L368,A419,V433,Y449,N450,D467,C480,E484,C488,F490,S494,R50901111
W353-R466W353,R355,R357,Y396,K462,F464,R46600000
R355-R466R355,R357,D428,K462,P463,F464,R46600000
R355-I498R355,R457,K462,F464,E465,R466,D467,I46800000
F374-P384F374,S375,T376,F377,K378,C379,P38400001
F374-E516F374,S375,T376,F377,C379,F392,D427,E51600001
R403-R509R403,D405,A419,G502,G504,Y505,R50900001
T415-G502T415,Y421,A475,G476,N487,S494,G50200000
K417-Y489K417,L455,F456,A475,E484,F486,Y48902211
K417-N501K417,Y449,L455,F456,Y489,G496,Q498,T500,N50112203
S443-S494S443,V445,G446,G447,Y449,P499,T50000000
K444-N450-S494K444,V445,G446,G447,Y449,N450,S49400000
K444-I472-S494K444,V445,G446,L452,L455,F456,T470,E471,I472,S49400010
K444-T500K444,V445,G446,N450,Q498,P499,T50000001
V445-G502V445,G446,Q498,P499,T500,N501,G50211102
G446-G447-F490G446,G447,N448,Y449,N450,L452,E484,F49001121
G446-Y449-F490G446,Y449,L452,V483,E484,G485,F49001121
G446-S494G446,Y449,E484,F490,L492,Q493,S49401112
Y449-F490Y449,I472,N481,V483,E484,G485,F486,F49001111
Y449-L492Y449,N481,G482,V483,E484,F490,L49201111
Y449-Q493Y449,E484,G485,F486,Y489,F490,Q49301112
Y449-L452-S494Y449,L452,T470,E484,F490,L492,S49401121
Y449-N481-S494Y449,N481,G482,V483,E484,F490,S49401111
Y449-Y505Y449,Y453,L455,F456,Y489,Q493,Y495,Y50500002
L452-Q490L452,F453,T475,G482,F483,N484,Y486,Q49000012
L455-Y489L455,F456,A475,T478,E484,F486,N487,Y48901122
L455-F490L455,V483,E484,G485,F486,Y489,F49001111
L455-Q493L455,A475,T478,G485,F486,Y489,Q49300012
F456-Q493F456,E484,G485,F486,C488,Y489,Q49301112
Y473-F486Y473,A475,G476,S477,E484,G485,F48601112
N556-D568N556,K557,K558,L560,P561,F562,Q563,D56800000
1 Clusters are named by the first and last residue in the epitope core. In those clusters where these residues coincide, the next residue at the beginning or the end of the core is added to the name.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ballesteros-Sanabria, L.; Pelaez-Prestel, H.F.; Ras-Carmona, A.; Reche, P.A. Resilience of Spike-Specific Immunity Induced by COVID-19 Vaccines against SARS-CoV-2 Variants. Biomedicines 2022, 10, 996. https://doi.org/10.3390/biomedicines10050996

AMA Style

Ballesteros-Sanabria L, Pelaez-Prestel HF, Ras-Carmona A, Reche PA. Resilience of Spike-Specific Immunity Induced by COVID-19 Vaccines against SARS-CoV-2 Variants. Biomedicines. 2022; 10(5):996. https://doi.org/10.3390/biomedicines10050996

Chicago/Turabian Style

Ballesteros-Sanabria, Laura, Hector F. Pelaez-Prestel, Alvaro Ras-Carmona, and Pedro A. Reche. 2022. "Resilience of Spike-Specific Immunity Induced by COVID-19 Vaccines against SARS-CoV-2 Variants" Biomedicines 10, no. 5: 996. https://doi.org/10.3390/biomedicines10050996

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop