Skip to main content

Sequence-based prediction of SARS-CoV-2 vaccine targets using a mass spectrometry-based bioinformatics predictor identifies immunogenic T cell epitopes

Abstract

Background

The ongoing COVID-19 pandemic has created an urgency to identify novel vaccine targets for protective immunity against SARS-CoV-2. Early reports identify protective roles for both humoral and cell-mediated immunity for SARS-CoV-2.

Methods

We leveraged our bioinformatics binding prediction tools for human leukocyte antigen (HLA)-I and HLA-II alleles that were developed using mass spectrometry-based profiling of individual HLA-I and HLA-II alleles to predict peptide binding to diverse allele sets. We applied these binding predictors to viral genomes from the Coronaviridae family and specifically focused on T cell epitopes from SARS-CoV-2 proteins. We assayed a subset of these epitopes in a T cell induction assay for their ability to elicit CD8+ T cell responses.

Results

We first validated HLA-I and HLA-II predictions on Coronaviridae family epitopes deposited in the Virus Pathogen Database and Analysis Resource (ViPR) database. We then utilized our HLA-I and HLA-II predictors to identify 11,897 HLA-I and 8046 HLA-II candidate peptides which were highly ranked for binding across 13 open reading frames (ORFs) of SARS-CoV-2. These peptides are predicted to provide over 99% allele coverage for the US, European, and Asian populations. From our SARS-CoV-2-predicted peptide-HLA-I allele pairs, 374 pairs identically matched what was previously reported in the ViPR database, originating from other coronaviruses with identical sequences. Of these pairs, 333 (89%) had a positive HLA binding assay result, reinforcing the validity of our predictions. We then demonstrated that a subset of these highly predicted epitopes were immunogenic based on their recognition by specific CD8+ T cells in healthy human donor peripheral blood mononuclear cells (PBMCs). Finally, we characterized the expression of SARS-CoV-2 proteins in virally infected cells to prioritize those which could be potential targets for T cell immunity.

Conclusions

Using our bioinformatics platform, we identify multiple putative epitopes that are potential targets for CD4+ and CD8+ T cells, whose HLA binding properties cover nearly the entire population. We also confirm that our binding predictors can predict epitopes eliciting CD8+ T cell responses from multiple SARS-CoV-2 proteins. Protein expression and population HLA allele coverage, combined with the ability to identify T cell epitopes, should be considered in SARS-CoV-2 vaccine design strategies and immune monitoring.

Background

Coronaviruses are positive-sense single-stranded RNA viruses that have occasionally emerged from zoonotic sources to infect human populations [1]. Most coronavirus infections cause mild respiratory symptoms. However, some recent coronavirus infections have resulted in serious morbidity and mortality, including the severe acute respiratory syndrome coronavirus (SARS-CoV) [2,3,4], Middle East respiratory syndrome coronavirus (MERS-CoV) [5, 6], and SARS-CoV-2, which are responsible for the current worldwide pandemic, COVID-19. These three viruses belong to the genus Betacoronaviridae [1]. SARS-CoV was identified in South China in 2002, and its global spread led to 8096 cases and 774 deaths [7]. The first case of MERS-CoV emerged in 2012 in Saudi Arabia, and since then, a total of 2494 cases and 858 associated deaths have been reported [6]. In contrast to the more limited scope of these other coronavirus infections, SARS-CoV-2, which emerged in Wuhan, China, at the end of December 2019, has resulted in 9,129,146 cases, including 473,797 deaths globally as of June 24, 2020 [8]. The rapid spread of SARS-CoV-2 has resulted in the World Health Organization declaring a global pandemic. Thus, there is an urgent need for effective vaccines and antiviral treatments against SARS-CoV-2 to reduce the spread of this highly infectious agent.

The genome of SARS-CoV-2 spans 30 kb in length and encodes for 13 open reading frames (ORFs), including four structural proteins. These structural proteins are the spike protein (S), the membrane protein (M), the envelope protein (E), and the nucleocapsid protein (N). In addition, there are over 20 non-structural proteins that account for all the proteins involved in the transcription and replication of the virus [9]. All encoded proteins of the virus are potential candidates for developing vaccines to induce robust T cell immunity.

SARS-CoV and SARS-CoV-2 share 76% amino acid identity across the genome [10, 11]. This high degree of sequence similarity allows us to leverage the previous research on protective immune responses to SARS-CoV to aid in vaccine development for SARS-CoV-2 [12,13,14,15]. Both humoral and cellular immune responses have been shown to be important in host responses to SARS-CoV [16]. Antibody responses generated against the S and the N proteins have shown to protect from SARS-CoV infection in mice and have been detected in SARS-CoV and SARS-CoV-2-infected patients [17,18,19,20]. However, the antibody responses detected against the S protein were undetectable in patients 6 years post-recovery [21]. In addition, higher titers of antibodies have been found in more severe clinical cases of viral infection suggesting that a robust antibody response alone may be insufficient for controlling SARS-CoV [22] and SARS-CoV-2 [23,24,25] infection.

Together with B cell immunity, T cell responses seem important in the immune response’s control of SARS-CoV and are also likely important for the control of SARS-CoV-2. In mice, studies have shown that adoptive transfer of SARS-CoV-specific memory CD8+ T cells provided protection against a lethal SARS-CoV infection in aged mice and that adoptive transfer of effector CD4+ and CD8+ T cells to immunodeficient or young mice expedited virus clearance and improved survival [26]. Both CD4+ and CD8+ T cell responses have also been detected in SARS-CoV [16, 27] and SARS-CoV-2-infected patients [28,29,30]. Additionally, SARS-CoV specific memory CD8+ T cells have been found to persist for up to 11 years post-infection in patients who recovered from SARS [31]. These viral specific CD8+ T cells can be cytotoxic and can kill virally infected cells to reduce disease severity [16]. In addition to having effector functions, CD4+ T cells can promote the production of virus-specific antibodies by activating T-dependent B cells. Given the wealth of data from SARS-CoV, the homology between the SARS-CoV-2 and SARS-CoV, as well as emerging data from SARS-CoV-2 [28, 32], T cell immunity likely plays a critical role in providing protection against SARS-CoV-2.

Here, we utilized mass spectrometry (MS)-based HLA-I and HLA-II epitope binding prediction tools [33, 34] to identify SARS-CoV-2 epitopes recognized by CD4+ and CD8+ T cells. These binding predictors were trained on high-quality mono-allelic HLA immunopeptidome data generated via MS. The use of MS for the identification of MHC peptide ligandome yields an extensive and relatively unbiased population of naturally processed and presented MHC binding peptides in vivo. Unlike traditional binding assays which rely on chemical synthesis and a priori knowledge of peptides and ligands to be assayed, MS uses natural peptide-MHC complexes which are subject to the endogenous processing and presentation pathways within the cell. Additionally, the use of engineered mono-allelic cell lines avoids dependence on in silico deconvolution techniques and allows for allele coverage to be expanded in a targeted manner. In Abelin et al. [34], we demonstrated that improved HLA-II binding prediction led to improved immunogenicity prediction by validating this approach on a data set of immune responses to a diverse collection of pathogens and allergens [35, 36].

With this approach, we generated binding predictors for 74 HLA-I and 83 HLA-II alleles (Additional file 1: Table S1 and Additional file 2: Table S2). Alleles selected for data collection were prioritized to maximize population coverage. Here, we specifically validated the binding predictors utilizing Coronaviridae family peptides that had been assayed for T cell reactivity or MHC binding from the Virus Pathogen Resource (ViPR) database [34]. The ViPR database integrates viral pathogen data from internally curated data, researcher submissions, and data from various external sources. Specifically, experimentally determined epitopes were derived from the Immune Epitope Database (IEDB) [37]. Compared with a recent study with a similar aim [38], our approach provides epitope predictions to a wider set of alleles, which were characterized using high-quality mono-allelic MS data. These include high-frequency alleles from diverse populations as well as lower frequency alleles, leading to an expansive set of bioinformatically predicted SARS-CoV-2 epitopes.

We used our HLA-I and HLA-II binding predictors to predict the binding potential of peptide sequences from across the entire SARS-CoV-2 genome for a broad set of HLA-I and HLA-II alleles, covering the vast majority of USA, European, and Asian populations (Additional file 3: Table S3). We additionally confirm that a subset of these epitopes can raise specific CD8+ T cell responses in T cell induction assays using donor PBMCs. Furthermore, we interrogate publicly available proteomic data and demonstrate that the relative expression of SARS-CoV-2 proteins in virally infected cells vary significantly, pointing to another parameter that should be considered in vaccine design to induce cellular immunity. Epitopes predicted to have a high likelihood of binding to multiple HLA-I and HLA-II alleles and exhibit high expression in infected human cells are promising vaccine candidates to elicit T cell responses against SARS-CoV-2.

Methods

Analysis of Coronaviridae family T cell epitopes from ViPR

Experimentally determined epitopes for the Coronaviridae family for human hosts were retrieved from the ViPR database (https://www.viprbrc.org/; accessed March 5, 2020) [39]. To build a validation dataset, both positives and negatives for T cell assays and MHC binding assays were obtained. Only assays associated with alleles identified with at least four-digit resolution and supported by our predictors (Additional file 1: Table S1) were included in this analysis. Positive calls were prioritized: peptide-allele pairs were classified as positive if a given peptide-allele pair was assayed multiple times by a specific assay type and was determined to be positive in any single one of the assays. Specifically, the priority was given by the following order: Positive-High > Positive-Intermediate > Positive-Low > Positive > Negative (e.g., a peptide allele pairing that was assayed three times with the results Positive-High, Positive, and Negative were assigned a Positive-High result). Of note, alternative approaches such as prioritizing negative assay results, or random choice in cases of multiple results, yielded very similar results.

Binding prediction for ViPR Coronaviridae family T cell epitopes

Peptide-HLA-I allele pairs in the ViPR validation dataset were scored using our HLA-I binding predictor, a neural network trained on mono-allelic MS data [33]. Similarly, peptide-HLA-II allele pairs in the ViPR validation dataset were scored using our HLA-II binding predictor, a recently published convolutional neural network-based model also trained on mono-allelic MS data [34]. We scored all 12–20mers contained within a given assay peptide with the HLA-II binding predictor and took the maximum score as the representative binding score for the assay peptide. In vitro MHC binding assays, which represent the vast majority of the ViPR dataset, do not require endogenous processing and presentation for a positive binding result. Since our binding predictor, which is trained on naturally processed and presented ligands observed via MS, is also implicitly learning these endogenous processing rules, we score all potential ligands within an assayed peptide (rather than just the full-length assay peptide itself) to account for this distinction.

Retrieval of SARS-CoV-2 sequence

The GenBank reference sequence for SARS-CoV-2 (accession: NC_045512.2, https://www.ncbi.nlm.nih.gov/nuccore/NC_045512) was used for this study. All twelve annotated open reading frames (ORF1a, ORF1b, S, ORF3a, E, M, ORF6, ORF7a, ORF7b, ORF8, N, and ORF10) were considered as sources of potential epitopes. In addition, due to its high expression level in recently published proteomic datasets [40,41,42], ORF9b, as annotated by UniProt (P0DTD2, https://www.uniprot.org/proteomes/UP000464024), was also used for epitope predictions.

Identification of HLA-I epitopes and prioritization by population coverage

To identify candidate HLA-I epitopes, we exhaustively scored all possible 8–12mer peptide sequences from SARS-CoV-2 with our HLA-I binding predictor [33] for 74 alleles, including 21 HLA-A alleles, 35 HLA-B alleles, and 18 HLA-C alleles. Peptide-allele pairs were assigned a percent rank by comparing their binding scores to those of 1,000,000 reference peptides (selected from a partition of the human proteome that had not been used for model training) for the same respective allele. Peptide-allele pairs that scored in the top 1% of the scores of these reference peptides were considered strong potential binders.

Since a vaccine should ideally benefit a large fraction of the population, these top-ranking peptides were then prioritized based on expected population coverage (allele frequencies obtained from [43]), given all the alleles each peptide was expected to bind to (i.e., all the alleles for which the peptide scored in the top 1%). The estimate of population coverage for each peptide was calculated as:

$$ \mathrm{Coverage}=1-{\Pi}_{\mathrm{loci}}\ {\left(1-{\Sigma}_{\mathrm{locus}\ \mathrm{alleles}}\ {f}_{\mathrm{allele},\mathrm{avg}}\right)}^2 $$

where fallele,avg is the (unweighted) average allele frequency across the US, European, and Asian-Pacific Islander (API) populations (which is intended to represent an approximation of a global population average, focusing on the populations most affected by the pandemic, and the cumulative product is taken across the three HLA-I loci: HLA-A, HLA-B, and HLA-C). The cumulative product itself represents the chance that an individual in the population does not express any one of the contained alleles; hence, the complement describes the probability that at least one is present.

The USA population allele frequency is calculated as the following weighted average of a few subpopulations: 0.623*EUR + 0.133*AFA + 0.068*API + 0.176*HIS, where EUR = European, AFA = African American, API = Asian-Pacific Islander, and HIS = Hispanic populations. These subpopulation frequencies are based on data from the US Census Bureau [44], accommodating for slight variations in different tables and information of mixed races. For alleles where AFA, HIS, or API population frequencies were not available, the US population allele frequency values were set to match EUR. Missing API allele frequency values were conservatively imputed with 0 for our analyses.

We then constructed two types of ranked lists of HLA-I epitopes by coverage. The first ranks all SARS-CoV-2 epitopes by their absolute coverage, such that peptides predicted to bind similar collections of alleles would be ranked similarly (Additional file 4: Table S4). This approach provides the full list of predicted class I epitopes sorted by the expected coverage for each peptide, with the generous assumption that every binding prediction is correct.

The second type of list, referred to as a “disjoint” list, is constructed in an iterative fashion where the peptide with the greatest coverage is selected first, and then, the coverage for the remaining epitopes is updated to nullify contributions from any alleles that have already been selected (Additional file 5: Table S5). Disjoint lists were generated for M, N, and S proteins (the most highly expressed structural proteins) individually, instead of across the entire SARS-CoV-2 genome, to provide protein-level prioritizations. This approach produces a parsimonious list of peptides that is designed to maximize cumulative population coverage with the fewest number of selections.

Identification of HLA-II epitopes and prioritization by population coverage

To identify HLA-II epitopes, we used our HLA-II binding predictor [34] to score all 12–20mer sequences in the SARS-CoV-2 proteome to predict both binding potential and the likely binding core within each 12–20mer. Scoring was performed across all supported HLA-II alleles, consisting of 46 HLA-DR alleles, 17 HLA-DP alleles, and 20 HLA-DQ alleles (Additional file 2: Table S2).

Peptide-allele pairs were assigned a percent rank by comparing their binding scores to those of 100,000 reference peptides (as before, sampled from a partition of the human proteome that was held out from training). Pairs scoring in the top 1% were deemed likely to bind. Additionally, we define the “epitope” of 12–20mers to be the predicted binding core within the sequence. As such, overlapping 12–20mers with the same predicted binding core for a given allele would constitute a single epitope.

Table 1 Summary of the HLA-I and HLA-II epitopes predicted across the 13 SARS-CoV-2 ORFs and their validation

Additionally, we prioritized predicted HLA-II binding 25mers in SARS-CoV-2 by population coverage, given the desire to design vaccines that are effective broadly across the global population. To do this, we associated each 25mer with all subsequences that were likely binders and calculated the population coverage of the corresponding HLA-II alleles. Given a collection of alleles, we calculated the coverage as described in the previous section, the only difference being the cumulative product is taken across the following four HLA-II loci: HLA-DRB1, HLA-DRB3/4/5, HLA-DP, and HLA-DQ. HLA-II allele frequencies were obtained from [43] and Allele Frequency Net Database [45].

As with HLA-I, two types of sorted lists of predicted binding sequences were generated. The first type ranks every predicted SARS-CoV-2 25mers by absolute coverage provided by the HLA-II alleles to which a constituent subsequence is expected to bind (Additional file 6: Table S6). The second type of ranking was again performed for predicted binders in M, N, and S proteins individually, using disjoint coverage, to maximize cumulative population coverage with a parsimonious list of peptides (Additional file 7: Table S7).

Comparison of predicted epitopes to the human proteome

Eight to 12mer sequences (corresponding to predicted HLA-I epitopes), 9mer sequences (corresponding to predicted HLA-II binding cores), and 25mer sequences (corresponding to predicted HLA-II sequences that bound multiple alleles) from SARS-CoV-2 were compared against all subsequences of the same length from the human proteome, using UCSC Genome Browser genes with hg19 annotation of the human genome and its protein coding transcripts (63,691 entries) [46]. Exact matches were identified, flagged in Additional file 4: Table S4, and omitted from the disjoint coverage ranking analysis to avoid prioritizing peptides that may inadvertently induce an autoimmune response. No exact matches were found for the predicted HLA-II binding cores or 25mer sequences.

T cell inductions and assessment of peptide-MHC-positive T cell responses

Human PBMCs from HLA-A02:01-positive human donors were isolated using Ficoll separation from apheresis material (AllCells, USA). Twenty-three SARS-CoV-2 peptides predicted to be strong binders to HLA-A02:01 were pooled by similar binding potential, with up to 6 peptides per pool. The selected peptides represent high ranking peptides predicted to bind HLA-A02:01 from across the S, N, M, E, and ORF1ab proteins, avoiding sequences also prioritized by Grifoni et al. [38]. Specifically, three of the 23 peptide sequences were chosen from ViPR and also scored highly in our predictions; the remaining 20 were the top-ranking peptides by our prediction for the abovementioned proteins. When comparing our 20 SARS-CoV-2-predicted epitopes with the ViPR dataset, 8 epitopes were previously assayed and confirmed as HLA-A02:01 binders. Of these 8, two were reported positive in a T cell assay in ViPR and two were negative. PBMCs were incubated with peptide pools, matured, and cultured in the presence of IL-7 and IL-15 (CellGenix GmbH, Germany) to promote T cell growth. Cells were then harvested, and the frequency of CD8+ T cells specific to peptide-MHC (pMHC) was assayed using combinatorial coding of pMHC multimers [43].

pMHC multimers were prepared as described previously [47, 48]. Briefly, biotinylated HLA-A02:01 monomers loaded with UV cleavable peptides were exchanged under UV light with SARS-CoV-2-predicted peptides. The streptavidin-labeled fluorophores PE, APC, BV421 (Biolegend, Inc., USA), BV650, and BUV395 (BD Biosciences, USA) were added to UV-exchanged monomers to create fluorescently labeled multimer reagents.

Harvested cells were then stained with LIVE/DEAD Fixable Near-IR Dead Cell Stain Kit for 633 or 635 nm excitation (Life Technologies Corporation, USA); anti-CD4 FITC, anti-CD14 FITC, anti-CD16 FITC, and anti-CD19 FITC (BD Biosciences, USA); and anti-CD8 AF700 (Biolegend Inc., USA). Only live CD8+ T cells staining for both fluorochromes of the relevant pMHC multimers were considered positive [47]. Samples were analyzed on FACS LSR Fortessa X20 cytometers (BD Biosciences), and data was analyzed using FlowJo (TreeStar).

Analysis of publicly available SARS-CoV-2 proteomic datasets

SARS-CoV-2 proteomic datasets were downloaded from the PRIDE repository (Bojkova et al. [40]: PXD017710, https://www.ebi.ac.uk/pride/archive/projects/PXD017710; Bezstarosti et al. [41]: PXD018760, https://www.ebi.ac.uk/pride/archive/projects/PXD018760; Davidson et al. [42]: PXD018241, https://www.ebi.ac.uk/pride/archive/projects/PXD018241). In these studies, either Caco-2 human colorectal adenocarcinoma cells [40] or Vero E6 African green monkey kidney epithelial cells [41, 42] were subject to infection with SARS-CoV-2. Tandem mass spectra (MS/MS) acquired with data-dependent acquisition (DDA) were interpreted using Spectrum Mill MS Proteomics software package v7.0 pre-release (Agilent Technologies). Cysteine carbamidomethylation was selected as a fixed modification. Methionine oxidation, asparagine deamidation, protein N-termini acetylation, peptide N-terminal glutamine to pyroglutamic acid, and peptide N-terminal cysteine pyro-carbamidomethylation were selected as variable modifications. For the dataset from Bojkova et al. [40] which employed isobaric mass tags, TMT11 was added as a fixed modification to peptide N-termini and lysines, and 13C6-15N2-TMT11-lysine and 13C6-15N4-arginine were added as variable modifications. All datasets were searched against the SARS-CoV-2 proteome (UniProtKB, 28 April 2020, 14 entries) concatenated to databases containing either the Homo sapiens proteome ([40], UCSC Genome Browser hg19 annotation, 63,691 entries) or the Chlorocebus sabaeus proteome ([41, 42] UniProtKB, 9229 entries). Precursor and fragment mass tolerances were set as described in each manuscript, or as 20 ppm when not specified. Database search results were exported as a list of peptide-spectrum matches (PSMs) with a target-decoy-based false discovery rate (FDR) estimation of 1%. Individual fractions from each study were combined into a single list. To perform spectral counting, PSMs assigned to a single SARS-CoV-2 protein were counted, with ORF1a and ORF1ab treated as a single protein group. Peptides matched to both a host and SARS-CoV-2 protein were discarded. Spectral counts were normalized to the length of each protein, and the maximum value within each dataset was set to 100%.

Results

Bioinformatics predictor validation for viral epitopes using ViPR

We first sought to validate the ability of our predictors to identify epitopes from genomes of the Coronaviridae family. Since SARS-CoV-2 only emerged recently, specific data on SARS-CoV-2 peptide MHC binding and immunogenic epitopes are currently limited. However, other viruses from the Coronaviridae family have been studied thoroughly, specifically MERS-CoV and SARS-CoV. The latter has significant sequence homology to SARS-CoV-2 [12]. We therefore sought to leverage previously tested epitopes from across the Coronaviridae family to validate our predictions of viral peptides, with special interest in peptide sequences that exactly matched protein sequences of the novel SARS-CoV-2 virus. To that end, we used the ViPR database [39], which lists the results of T cell immunogenicity and MHC-peptide binding assays for both HLA-I and HLA-II alleles for viral pathogen epitopes. We used all assays of Coronaviridae family viruses with human hosts from ViPR as our validation dataset. Assays that did not have an associated four-digit HLA allele or were associated with an allele our models did not support were omitted (see Additional file 1: Table S1 and Additional file 2: Table S2 for a list of supported alleles).

For HLA-I, within the validation dataset, there were a total of 4445 unique peptide-HLA allele pairs that were assayed for MHC binding, using variations of (1) cellular MHC or purified MHC, (2) a direct or competitive assay, and (3) measurement by fluorescence or radioactivity. Two additional peptide-MHC allele pairs were confirmed via X-ray crystallography. Depending on the study from which the data was collected, peptide-MHC allele pairs were either defined in ViPR simply as “Negative” and “Positive” for binding, or with a more granular scale of positivity: low, intermediate, and high. We assigned peptide-MHC allele pairs with multiple measurements with the highest MHC binding detected across the replicates (see the “Methods” section).

We then applied our HLA-I binding predictor to the peptide-MHC allele pairs in the validation dataset and compared the computed HLA-I percent ranks of these pairs with the reported MHC binding assay results (Additional file 8: Table S8). A low percent rank value corresponds to high likelihood of binding (e.g., a peptide with a percent rank of 1% scores among the top 1% in a reference set of random peptides). The percent ranks of peptide-MHC allele pairs that had a binary “Positive” result in the MHC binding assay were significantly lower than pairs with a “Negative” result. Further, in the more granular positive results, stronger assay results (low < intermediate < high) were associated with increasingly lower percent ranks (Fig. 1a). In addition, the two peptide-MHC alleles that were confirmed by X-ray crystallography were predicted as very likely binders, with low percent rank scores of 0.07% and 0.30%. Although our HLA-I binding predictor was initially built with the purpose of supporting neoantigen prediction in cancer, this analysis shows that it can be successfully applied to coronavirus proteomes. We evaluated our predictor by performing a Precision-Recall analysis, demonstrating the tradeoff between accurate calling of positive binders and the fraction of true binders that are detected (Fig. 1b).

Fig. 1
figure1

Binding predictions for both peptide-HLA-I and HLA-II pairs from ViPR correlate with their reported assay results. a The log10(percent rank) of scored peptide-HLA-I allele pairs, versus their ViPR reported MHC binding assay result (either binary Negative/Positive or the scaled Negative/Positive-Low/Positive-Intermediate/Positive-High, based on the reported value). In total, there were 4445 peptide-HLA-I allele pairs in the ViPR dataset we obtained (see the “Methods” section). Black lines indicate median values. b A modified precision-recall analysis for the HLA-I binding prediction of ViPR data, in which we demonstrate the fraction of the true positives out of all called positives (precision, indicated by color) and the fraction of detected true positives out of all true positives (recall, indicated x-axis value) as a fraction of the log10(percent-rank) threshold (y-axis value). Red dashed line indicates an example of log10(percent-rank) threshold of − 1, corresponding to the line in a. Area under the precision-recall curve is indicated. c The log10(percent rank) of scored peptide-HLA-II allele pairs, versus their ViPR reported MHC binding assay result (Negative+Positive-Low/Positive-Intermediate/Positive-High, based on the reported value). In total, there were 259 peptide-HLA-II allele pairs. Black lines indicate median values. d A modified precision-recall analysis for the HLA-II binding prediction of ViPR data, in which we demonstrate the fraction of the true positives out of all called positives (precision, indicated by color) and the fraction of detected true positives out of all true positives (recall, indicated x-axis value) as a fraction of the log10(percent-rank) threshold (y-axis value). Red dashed line indicates an example of log10(percent-rank) threshold of − 1, corresponding to the line in c. Area under the precision-recall curve is indicated

Assays of T cell reactivity (e.g., interferon-gamma ELISpots, tetramers), which are stricter measures for T cell immunogenicity to epitopes, were performed in significantly lower numbers compared with MHC binding assays. For HLA-I, the overlap between peptide-MHC allele pairs for which we had a prediction (supported alleles) and pairs with a reported T cell assay consisted of only 32 pairs, of which 23 had a positive result. We did not detect differences in the percent ranks across the positive and negative groups; however, sample sizes are extremely small. In addition, for HLA-I epitopes, the validation dataset only contained T cell assay results for peptide-MHC allele pairs that had a positive result in a binding assay, suggesting a highly biased pool of epitopes selected for testing, as also reflected in the high rate of positive T cell assay results. Indeed, the high rate of positive MHC binding assays compared to what would be expected for completely randomly selected peptides also implies that peptides expected to bind based on prediction or prior data were prioritized for testing (or negative results were under-reported). This underlying bias in peptides assayed is important to keep in mind in evaluating the binding predictor performance on this validation dataset. An even more dramatic difference in scores for positives versus negatives could be expected had random peptides been selected for testing.

In addition to the identification of targets for CD8+ T cells, we have recently demonstrated the ability to predict HLA-II binders [34], allowing us to target CD4+ T cell responses which could be harnessed for SARS-CoV-2 vaccines. These CD4+ responses can potentially bolster both T cell immunity and enhance humoral immunity [49].

In a similar fashion to the HLA-I analysis, we scored all Coronaviridae family peptide-MHC allele pairs with supported HLA-II alleles in ViPR using our HLA-II binding predictor [34] (Additional file 9: Table S9). There were 259 unique peptide-MHC allele pairs assayed by MHC binding assays in the ViPR validation dataset for HLA-II. As before, we compared their percent rank with their reported “best” (in the case of multiple measurements) MHC binding assay result. This comparison could not be performed with the “Negative” pairs as an independent group since there was only one negative result in the validation dataset for HLA-II. The low negative counts may be due to under-reporting of negative assay results or biased selection of the peptides to be assayed. Therefore, we merged the “Negative” and “Positive-Low” groups into one group and compared their percent ranks with either the “Positive-Intermediate” or the “Positive-High” groups (Fig. 1c). This analysis revealed a trend similar to that observed with HLA-I predictions, indicating that stronger MHC binding assay results are associated with a lower predicted percent rank for HLA-II binders, as we expect for a robust predictor. We also evaluated our HLA-II binding predictor by performing a precision-recall analysis (Fig. 1d). The area under the precision-recall curve (AUC) indicated only a small advantage to our predictor over a random guess, which is explained by the heavy bias towards peptides with positive HLA-II binding assay results. Similar to the HLA-I T cell assays, there were too few recorded HLA-II T cell assays in our validation dataset to determine percent rank differences between peptide-HLA II allele pairs testing positive and negative. Together, these findings further corroborate the validity of our epitope predictors, as peptide-MHC allele pairs with positive results in binding assays consistently have lower percent ranks (better scores) by both our HLA-I and HLA-II MHC binding predictors.

Epitope prediction for SARS-CoV-2

We harnessed our HLA binding prediction ability to identify the peptides most relevant to the generation of SARS-CoV-2 T cell responses. We first performed the analysis for HLA-I peptide binding and computed the likelihood of each peptide of lengths 8–12 amino acids from the 13 SARS-CoV-2 ORFs to bind to any HLA-I allele in our database. We then calculated the percent rank of each peptide-MHC allele pair by comparing their binding scores to those of a set of reference peptides; putative binders were identified as sequences predicted to bind to a given allele with a percent rank of 1% or lower (Fig. 2 a–c).

Fig. 2
figure2

A schematic demonstrating our approach to identify SARS-CoV-2 T cell epitopes and their validation. (a) A diagram of the SARS-CoV-2 virus, listing example proteins. (b) Applying our HLA-I and HLA-II binding predictors to the 13 annotated ORFs of SARS-CoV-2 (including ORF9b). (c) Both HLA-I and HLA-II epitopes are ranked by their likelihood to bind a particular HLA allele. (d) Epitopes shared between SARS-CoV-2 and other coronaviruses which were previously assayed are used for validation. (e) A description of a T cell induction assay to assess immunogenicity of select epitopes

By this metric, we detected a total of 11,897 unique SARS-CoV-2 peptides that were predicted to bind at least one HLA-I allele (Additional file 4: Table S4). Sixteen of these peptides overlapped with a subsequence of the human proteome and were marked for considerations of potential autoimmunity (see the “Methods” section, Additional file 4: Table S4).

Unlike HLA-I, which has a closed binding groove that constrains bound peptide lengths to approximately 8 to 12 amino acids, peptides binding HLA-II have a wider length distribution (up to 30 amino acids or even longer) since the HLA-II binding groove is open at both ends. Peptides bind with a 9-amino acid subsequence (termed the binding core) occupying the HLA-II binding groove, with any flanking sequence overhanging the edges of the molecule. We consider a group of peptides that differ in the flanking regions but share a common binding core as a single epitope. Using the HLA-II predictor, we identified 3372 unique binding cores that are predicted to bind at least one HLA-II allele with a percent rank score of 1% or lower (Table 1). The majority of predicted peptide-MHC allele pairs are from ORF1a and ORF1ab, primarily driven by the length of these ORFs. In addition, ORF1a and ORF1ab have very similar sequences, with over 18,000 identical binding peptide-HLA-I allele pairs predicted for both ORFs. We therefore opted to exclude redundant predictions and only reported unique pairs (see * in Table 1). Similarly, all HLA-II predicted epitopes from ORF1a were covered by those reported for ORF1ab.

To test the validity of the SARS-CoV-2-predicted peptide-HLA pairs, we looked for peptide sequences in the Coronaviridae portion of the ViPR database which exactly matched SARS-CoV-2 peptide sequences (Fig. 2d). A total of 374 HLA-I peptide-MHC allele pairs from SARS-CoV-2 had both a percent rank lower than 1% by our predictor and were found in the HLA-I MHC binding validation dataset. Strikingly, of these HLA-I peptide-MHC allele pairs, 333 (89%) had a positive assay result. As a comparison, we also tested for overlap between epitopes predicted to have low likelihood of MHC binding (percent rank 50% or higher) and the validation dataset. Thirty-seven peptide-MHC allele pairs overlapped between these sets, of which 36 (97.2%) had a negative assay result, as predicted. Further, we sought to determine whether our highly predicted SARS-CoV-2 peptide-HLA-I allele pairs (percent rank lower than 1%) would be validated by reported T cell assay results. Despite the significantly smaller number of peptide-MHC allele pairs that were tested for T cell reactivity in the validation dataset, 10 assayed pairs were also highly predicted by our HLA-I binding predictor. Nine out of these 10 (90%) predicted pairs had a positive result to the T cell assay. No low-scoring pairs (percent rank of 50% or above) were reported in the validation dataset. These findings demonstrate the validity of our prediction for peptide-HLA-I allele pairs for SARS-CoV-2 epitopes. Notably, while our algorithms are not trained on T cell reactivity data and are aimed at peptide-MHC binding, for the few examples in ViPR for which T cell reactivity assay results were reported, we were able to show our highly scoring peptide-MHC allele pairs are indeed immunogenic in the vast majority of cases.

For HLA-II peptide-MHC allele pairs, only a single HLA-II peptide-MHC allele pair had both a percent rank lower than 1% and was reported in the validation dataset; this single pair (from the envelope protein) had a “Positive-High” assay result.

Immunogenicity of HLA-A02:01-predicted SARS-CoV-2 epitopes

Our binding prediction algorithms predict the likelihood of an epitope to be presented by a specific HLA allele, but do not directly predict the ability of a T cell receptor to recognize the epitope presented by the MHC molecule. Due to the process of central tolerance, which deletes T cells that could cross-react with peptides from self-antigens, not every epitope that is a strong MHC binder will elicit a T cell response [50]. Therefore, there is a need to further validate high affinity MHC binding peptides in T cell assays as previously described (Fig. 2e) [51,52,53,54]. To address the immunogenicity of a subset of highly predicted MHC binding peptides, we synthesized 23 highly predicted HLA-A02:01 binding epitopes from each of the following SARS-CoV-2 proteins: S, M, N, E, and ORF1ab (Fig. 3a). Of these, 20 were selected solely due to being highly predicted SARS-CoV-2 epitopes while the additional three are also highly predicted but were chosen from ViPR. Pools of these peptides were cultured with PBMCs from three human donors, and the predicted epitopes were considered immunogenic if they elicited a T cell response as detected by binding to pMHC multimers for HLA-A02:01 in at least one of three donors.

Fig. 3
figure3

Experimental Validation of HLA-A02:01 predicted epitopes from SARS-CoV-2 in human T cell induction assays. a 23 peptides that were predicted to be high binders to HLA-A02:01 were synthesized and assayed in T cell inductions using PBMCs from three human donors. Three epitopes marked with asterisk were chosen based on ViPR, while the remaining 20 were chosen solely based on the predictor score. For our assay, epitopes were considered to be immunogenic if at least one donor raised a T cell response to the peptide as determined by pMHC multimer technology. ViPR confirmation refers to identical sequences from SARS-CoV confirmed via either MHC binding or T cell assays. b Flow cytometry plots of pMHC multimer staining from representative immunogenic SARS-CoV-2 epitopes. Multimer positive populations are circled in red, with the frequency of multimer positive CD8+ T cells shown in the upper right-hand corner of each plot

Overall, we detected CD8+ T cell responses in at least one donor for 11 of the 23 highly predicted epitopes in our assay (Fig. 3a, b). Fifteen of the 20 epitopes (75%) selected solely based on the prediction score were either reported positive in ViPR for MHC binding (n = 8) or had T cell reactivity in our assay (n = 8) in association with HLA-A02:01. Importantly, 10 of these 20 predicted epitopes (50%) tested positive for T cell reactivity, either in our assay (n = 8) or in previous reports of SARS-CoV (n = 2), confirming that our binding predictor can identify epitopes that are immunogenic. We were thus able to identify eight novel epitopes not previously reported in ViPR that were recognized by specific CD8+ T cells in donor PBMCs. The responses were generally robust, with nine of the 11 epitopes positive in our assay being recognized by specific CD8+ T cell responses in at least two donors (Fig. 3a), and encouragingly, every ORF from SARS-CoV-2 that was assayed had at least one peptide that led to a T cell response (Fig. 3a). Taken together, these data show that many novel SARS-CoV-2 epitopes that were predicted to be strong binders from our HLA-I binding predictor were found to be immunogenic.

Population coverage of peptides predicted to bind multiple HLA-I and HLA-II alleles

We sought to prioritize peptides from the M, N, and S proteins that are predicted to provide broad coverage for the US, European, and Asian-Pacific Islander populations based on the prevalence of MHC alleles in these populations [43]. We found that a subset of the peptides was predicted to bind a broad set of either HLA-I or HLA-II alleles. For each protein, we determined that a small number of peptide sequences provide saturating coverage for the US, European, and Asian-Pacific Islander populations, with > 99% population coverage achieved with selected 8–12mer epitopes for HLA-I, and > 95% population coverage achieved with selected 25mer sequences for HLA-II, respectively (Fig. 4a, b). Even if the generous assumption that all peptide-MHC allele pairs for which a given peptide scores in the top 1% are indeed immunogenic is not fully upheld, this finding could facilitate the design of a parsimonious, broadly effective vaccine to induce broad T cell immunity.

Fig. 4
figure4

Few predicted multi-allele binding epitopes from individual SARS-CoV-2 proteins can achieve broad population coverage. a Cumulative HLA-I coverage for USA, EUR, and API populations versus the number of included prioritized HLA-I epitopes for M, N, and S proteins, respectively. See Additional file 5: Table S5 for the peptide sequences corresponding to each panel. b Cumulative HLA-II coverage for each population versus the number of included prioritized HLA-II 25mers for M, N, and S proteins, respectively. See Additional file 7: Table S7 for the peptide sequences corresponding to each panel

Leveraging proteomic data to infer relative viral protein abundance

In addition to peptide-MHC binding, another important consideration in the design of a potential SARS-CoV-2 vaccine is the degree of viral protein expression in infected host cells. In order to determine the relative abundance of SARS-CoV-2 proteins, we analyzed three publicly available proteomic datasets that acquired unbiased LC-MS/MS on tryptic digestions of SARS-CoV-2-infected host cells [40,41,42]. Relative abundance of the viral proteins was estimated by spectral counting, a semi-quantitative approach whereby peptide-spectrum matches are counted, and totals are compared across proteins (Additional file 10: Table S10) [55, 56]. This analysis demonstrated the significantly wide range of expression levels of the SARS-CoV-2 proteins. Specifically, it confirmed that the N protein is the most abundant viral protein across all three datasets following SARS-CoV-2 infection (Fig. 5). This finding is corroborated by reports of N-derived peptides being detected in gargle solution samples from COVID-19 patients [57]. Furthermore, the N protein has been used as a biomarker for diagnosing patients infected with the SARS-CoV virus [58]. On the other hand, based solely on genomic information, ORF10 might be considered a potential target for vaccine development. However, there is very little proteomic and transcriptomic evidence that ORF10 is actually expressed in SARS-CoV-2-infected cells [42, 59]. These findings emphasize the value of considering SARS-CoV-2 protein expression levels in addition to HLA binding predictions and the immunogenicity of these epitopes in vaccine design strategies.

Fig. 5
figure5

Analysis of publicly available proteomic datasets demonstrates relative SARS-CoV-2 protein expression levels. Three datasets examining the proteomic response to SARS-CoV-2 infection were re-analyzed, and protein abundance was estimated by spectral counts normalized to protein length. Any annotated ORF not shown in the figure was not detected in these proteomic studies. Across all three studies, the nucleocapsid protein (N) is the most abundant SARS-CoV-2 protein in infected cells

Discussion

In this work, we demonstrated the utility and validity of our HLA-I and HLA-II binding prediction algorithms to the Coronaviridae virus family, and specifically to SARS-CoV-2. We use our validated predictors trained on mono-allelic MS data for both HLA-I and HLA-II binders, which potentially could be leveraged to induce both long-term CD4+ and CD8+ T cell immunity against the virus. Specifically, our HLA-II predictor, which has been trained on a large set of mono-allelic MS data and has been shown to identify immunogenic epitopes, is used here to identify high-quality SARS-CoV-2 CD4+ epitopes that may contribute to both cellular and humoral immunity [34] (Additional file 6: Table S6). Our database of supported HLA-I and HLA-II alleles provides us with the ability to not only identify many peptide-MHC allele pairs, but to generate a narrow list of peptides with many potential HLA pairings that could be presented by the entire US, European, and Asian-Pacific Islander populations. By applying these algorithms to previously assayed peptide-MHC allele pairs in ViPR, we were able to demonstrate excellent concordance between our binding predictions and the results of the binding assays for both HLA-I and HLA-II epitopes. We leveraged the homology within the Coronaviridae family to demonstrate that an exceedingly high portion (~ 90%) of our high-ranking SARS-CoV-2 peptide-MHC allele pairs for which validation was available was indeed confirmed to bind the predicted MHC allele.

We also experimentally confirmed that our binding predictors can identify epitopes that are immunogenic and can lead to CD8+ T cell responses to multiple SARS-CoV-2 proteins in donor PBMCs. It is plausible that the significant fraction of epitopes we experimentally confirmed (50% of highly predicted, tested epitopes) is only an underestimate for overall immunogenicity, since PBMCs from only three donors were used in this initial experiment. While an immunogenicity rate of 50% for epitopes predicted based on HLA binding is encouraging, it is challenging to compare with previous studies due to differences in T cell induction protocols, prediction algorithms and their prediction thresholds, and epitope etiology [53, 60, 61]. Pre-prints released during the revision of this manuscript [60, 61] studied CD8+ T cell responses to a small set of epitopes from the S protein in COVID-19 patient or healthy donor PBMCs. The study by Shomuradova et al. [60] showed that in their cohort, following stimulation, two HLA-A02:01 epitopes (YLQPRTFLL and RLQSLQTYV) distinguished CD8+ T cell responses in COVID-19 patients from healthy donor samples collected before or during the pandemic. Their data also show that 10 additional peptides sporadically induce responses in 2–3 individuals from either the healthy donor or the COVID-19 patient cohorts (total of 31 individuals). We found that 11 of the 12 peptides are highly ranked by our predictor (percent rank < 1%), with the two peptides specifically immunogenic in patients scoring at 0.1 and 0.17% rank, respectively. Interestingly, in our T cell assay, the HLA-A02:01 epitope RLQSLQTYV elicited a CD8+ T cell response in all three healthy donors.

While focusing on a small set of two COVID-19 patients and one healthy donor, Chour et al. [61] demonstrate that five HLA-A02:01-restricted epitopes from the S protein elicited a CD8+ T cell response in all three individuals. Although these epitopes were not tested in our immunogenicity assay, all five (100%) were predicted highly by our HLA-I binding predictor (percent rank < 1%). These studies provide independent confirmation to our algorithm’s ability to not only predict HLA-I binding, but also identify immunogenic peptides.

Though we did not perform T cell assays to evaluate the immunogenicity of the HLA-II predicted epitopes, such analysis would be valuable, especially given the importance of CD4+ T cells in both the cellular and humoral anti-viral response. We thus propose that a combination of B and T cell epitopes could provide long-lasting immunity from SARS-CoV-2 or mitigate the severity of disease when protection is partial.

We therefore concluded that using our HLA binding predictors to predict T cell epitopes from the ORFs of SARS-CoV-2 provides a novel and large set of high-quality T cell vaccine targets for the virus. In comparison to the recent publication by Grifoni et al. [38], we provide a large number of predicted epitopes, which we attribute to two differences in our approaches: (a) we opted to use a less stringent prediction cutoff based on the concordance of our 1% rank cutoff with the previously reported epitopes in ViPR and (b) we provide predictions for alleles from a wider range of population frequencies. This approach better covers non-white populations and provides predictions for rare alleles based on models trained on mono-allelic MS data as opposed to the extrapolation required by pan-allele predictors for alleles lacking sufficient data. This approach allows us to better prioritize vaccine candidates, as well as provide researchers with predictions to investigate materials from individuals with less frequent alleles. The differences in training data lead to algorithms that prioritize substantially different epitope sets, even for well-studied alleles. For instance, of the 91 HLA-A02:01 epitopes nominated by Grifoni et al., only 48 overlap with the top 91 from our predictor. In addition, we provide not only bioinformatics validation based on previously reported T cell epitopes and MHC binding peptides from other viruses from the Coronaviridae family in ViPR, but also experimentally validated, novel SARS-CoV-2 T cell epitopes.

The selection of target sequences can be further guided by levels of protein expression, predicted population coverage, and degree of sequence conservation. First, designing therapeutics against predicted epitopes is only effective if the proteins containing those epitopes are expressed at high enough levels for efficient antigen processing and presentation to take place. Therefore, it is crucial that protein expression be considered when selecting therapeutic targets. Second, prioritization of epitopes that are predicted to bind multiple alleles could provide coverage to significant fractions of the population, while including few epitopes in the vaccine. Lastly, during the viral spread and expansion through the population, genomic modifications are acquired, generating sequence diversity among the SARS-CoV-2 variants. This diversity may allow evasion of immune pressure, and therefore, it is important to prioritize epitopes that are conserved across the SARS-CoV-2 variants [62]. Novel tools which enable restricting peptides to conserved regions have recently become available: https://covidep.ust.hk [63] and http://cov-glue.cvr.gla.ac.uk/#/home [64].

Limiting epitope selection to highly expressed proteins, epitopes predicted to bind multiple high-frequency HLA alleles, and conserved viral sequences restricts the number of potential epitopes. However, the breadth of the list we provide increases the likelihood of identifying many high-quality, highly expressed epitopes. The epitopes characterized here, combined with insights on SARS-CoV-2 protein expression along with further efforts to confirm immunogenicity, can provide pre-clinical validation of epitopes that may be vaccine candidates to induce strong cellular immunity.

Conclusions

In summary, our work provides the most extensive set of both CD4+ and CD8+ T cell epitopes that are spanning the entire SARS-CoV-2 genome and binding a wide set of HLA-I and HLA-II alleles. Combining this epitope list with consideration of protein expression levels, population coverage and viral sequence conservation will lead to generation of a short list of vaccine epitope candidates that are likely immunogenic in the majority of the population. Our predicted list of CD4+ and CD8+ T cell epitopes will complement B cell epitopes and serve as a resource for the scientific community to generate potent SARS-CoV-2 vaccine epitopes and generate long-lasting T cell immunity.

Availability of data and materials

All data generated or analyzed during this study are included in this published article, its supplementary information files, or the following external sources: SARS-CoV-2 reference sequences used in this study were obtained from GenBank (accession: NC_045512.2, https://www.ncbi.nlm.nih.gov/nuccore/NC_045512) and ORF9b, as annotated by UniProt (P0DTD2, https://www.uniprot.org/proteomes/UP000464024). The method has been described in depth in Abelin et al., Immunity 2017 [33] and Abelin et al., Immunity 2019 [34]. Analogous online tools to the ones deployed here are available at (a) http://hlathena.tools/ for HLA I epitopes which is maintained by the Broad institute, and described in Sarkizova et al., Nature Biotechnology 2019 alongside with the associated data [65], and (b) https://neonmhc2.org/ for HLA II epitopes [34]. Mono-allelic MS data used for the training of our HLA-II binding predictor is also available as part of Abelin et al., Immunity 2019 [34]. SARS-CoV-2 proteomic datasets were downloaded from the PRIDE repository (Bojkova et al. [40]: PXD017710, https://www.ebi.ac.uk/pride/archive/projects/PXD017710; Bezstarosti et al. [41]: PXD018760, https://www.ebi.ac.uk/pride/archive/projects/PXD018760; Davidson et al. [42]: PXD018241, https://www.ebi.ac.uk/pride/archive/projects/PXD018241). Custom R and Python scripts used in generation of supplementary tables and figures are included (Additional files 11, 12 and 13).

Abbreviations

HLA:

Human leukocyte antigen

MERS-CoV:

Middle East respiratory syndrome coronavirus

MHC:

Major histocompatibility complex

MS:

Mass spectrometry

PBMC:

Peripheral blood mononuclear cells

PSM:

Peptide-spectrum matches

SARS-CoV:

Severe acute respiratory syndrome coronavirus

SARS-CoV-2:

Severe acute respiratory syndrome coronavirus 2

USA:

United States of America

ViPR:

Virus Pathogen Resource

WHO:

World Health Organization

References

  1. 1.

    Cui J, Li F, Shi Z-L. Origin and evolution of pathogenic coronaviruses. Nat Rev Microbiol. 2019;17(3):181–92.

    CAS  PubMed  Google Scholar 

  2. 2.

    Gu J, Gong E, Zhang B, Zheng J, Gao Z, Zhong Y, et al. Multiple organ infection and the pathogenesis of SARS. J Exp Med. 2005;202(3):415–24.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Rota PA, Oberste MS, Monroe SS, Nix WA, Campagnoli R, Icenogle JP, et al. Characterization of a novel coronavirus associated with severe acute respiratory syndrome. Science. 2003;300(5624):1394–9.

    CAS  PubMed  Google Scholar 

  4. 4.

    Ksiazek TG, Erdman D, Goldsmith CS, Zaki SR, Peret T, Emery S, et al. A novel coronavirus associated with severe acute respiratory syndrome. N Engl J Med. 2003;348(20):1953–66.

    CAS  PubMed  Google Scholar 

  5. 5.

    WHO MERS-CoV Global summary and assessment of risk. World Health Organization; 2018. Available from: https://www.who.int/csr/disease/coronavirus_infections/risk-assessment-august-2018.pdf.

  6. 6.

    WHO | Middle East respiratory syndrome coronavirus (MERS-CoV). WHO. World Health Organization; [cited 2020 Mar 25]. Available from: http://www.who.int/emergencies/mers-cov/en/.

  7. 7.

    WHO | Update 49 - SARS case fatality ratio, incubation period. WHO. World Health Organization; [cited 2020 Mar 25]. Available from: https://www.who.int/csr/sars/archive/2003_05_07a/en/.

  8. 8.

    Coronavirus disease (COVID-19) Pandemic. [cited 2020 Mar 25]. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019.

  9. 9.

    Chan JF-W, Kok K-H, Zhu Z, Chu H, To KK-W, Yuan S, et al. Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan. Emerg Microbes Infect. 2020;9(1):221–36.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Zhou P, Yang X-L, Wang X-G, Hu B, Zhang L, Zhang W, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020;579(7798):270–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Lu R, Zhao X, Li J, Niu P, Yang B, Wu H, et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet. 2020;395(10224):565. 74.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Ahmed SF, Quadeer AA, McKay MR. Preliminary identification of potential vaccine targets for the COVID-19 coronavirus (SARS-CoV-2) based on SARS-CoV immunological studies. Viruses. 2020 25;12(3).

  13. 13.

    Zheng M, Song L. Novel antibody epitopes dominate the antigenicity of spike glycoprotein in SARS-CoV-2 compared to SARS-CoV. Cellular & Molecular Immunology. 2020;17(5):536–8.

    CAS  Google Scholar 

  14. 14.

    Wrapp D, Wang N, Corbett KS, Goldsmith JA, Hsieh C-L, Abiona O, et al. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science. 2020;367(6483):1260–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Walls AC, Park Y-J, Tortorici MA, Wall A, McGuire AT, Veesler D. Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell. 2020;181(2):281–92. e6.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Li G, Fan Y, Lai Y, Han T, Li Z, Zhou P, et al. Coronavirus infections and immune responses. J Med Virol. 2020;92(4):424–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Yang Z-Y, Kong W-P, Huang Y, Roberts A, Murphy BR, Subbarao K, et al. A DNA vaccine induces SARS coronavirus neutralization and protective immunity in mice. Nature. 2004;428(6982):561–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Graham RL, Becker MM, Eckerle LD, Bolles M, Denison MR, Baric RS. A live, impaired-fidelity coronavirus vaccine protects in an aged, immunocompromised mouse model of lethal disease. Nat Med. 2012;18(12):1820–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Wang J, Wen J, Li J, Yin J, Zhu Q, Wang H, et al. Assessment of immunoreactive synthetic peptides from the structural proteins of severe acute respiratory syndrome coronavirus. Clin Chem. 2003;49(12):1989–96.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Liu X, Shi Y, Li P, Li L, Yi Y, Ma Q, et al. Profile of antibodies to the nucleocapsid protein of the severe acute respiratory syndrome (SARS)-associated coronavirus in probable SARS patients. Clin Diagn Lab Immunol. 2004;11(1):227–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Tang F, Quan Y, Xin Z-T, Wrammert J, Ma M-J, Lv H, et al. Lack of peripheral memory B cell responses in recovered patients with severe acute respiratory syndrome: a six-year follow-up study. J Immunol. 2011;186(12):7264–8.

    CAS  PubMed  Google Scholar 

  22. 22.

    Zhang L, Zhang F, Yu W, He T, Yu J, Yi CE, et al. Antibody responses against SARS coronavirus are correlated with disease outcome of infected individuals. J Med Virol. 2006;78(1):1–8.

    CAS  PubMed  Google Scholar 

  23. 23.

    Zhao J, Yuan Q, Wang H, et al. Antibody responses to SARS-CoV-2 in patients of novel coronavirus disease 2019 [published online ahead of print, 2020 Mar 28]. Clin Infect Dis. 2020. https://doi.org/10.1093/cid/ciaa344.

  24. 24.

    Okba NMA, Müller MA, Li W, Wang C, GeurtsvanKessel CH, Corman VM, et al. Early release - severe acute respiratory syndrome coronavirus 2−specific antibody responses in coronavirus disease 2019 patients - Volume 26, Number 7—July 2020 - Emerging Infectious Diseases journal - CDC. [cited 2020 May 11]; Available from: https://wwwnc.cdc.gov/eid/article/26/7/20-0841_article.

  25. 25.

    Cao X. COVID-19: immunopathology and its implications for therapy. Nat Rev Immunol. 2020;20(5):269–70.

    CAS  PubMed  Google Scholar 

  26. 26.

    Channappanavar R, Zhao J, Perlman S. T-cell-mediated immune response to respiratory coronaviruses. Immunol Res 2014;59(0):118–128.

  27. 27.

    Li CK, Wu H, Yan H, Ma S, Wang L, Zhang M, et al. T cell responses to whole SARS coronavirus in humans. J Immunol. 2008;181(8):5490–500.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Thevarajan I, Nguyen THO, Koutsakos M, Druce J, Caly L, van de Sandt CE, et al. Breadth of concomitant immune responses prior to patient recovery: a case report of non-severe COVID-19. Nat Med. 2020;16:1–3.

    Google Scholar 

  29. 29.

    Ni L, Ye F, Cheng M-L, Feng Y, Deng Y-Q, Zhao H, et al. Detection of SARS-CoV-2-specific humoral and cellular immunity in COVID-19 convalescent individuals. Immunity. 2020;S1074761320301813.

  30. 30.

    Weiskopf D, Schmitz KS, Raadsen MP, Grifoni A, Okba NMA, Endeman H, et al. Phenotype and kinetics of SARS-CoV-2–specific T cells in COVID-19 patients with acute respiratory distress syndrome. Science Immunology. 2020 [cited 2020 Jul 17];5(48). Available from: https://immunology.sciencemag.org/content/5/48/eabd2071.

  31. 31.

    Ng O-W, Chia A, Tan AT, Jadi RS, Leong HN, Bertoletti A, et al. Memory T cell responses targeting the SARS coronavirus persist up to 11 years post-infection. Vaccine. 2016;34(17):2008–14.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Grifoni A, Weiskopf D, Ramirez SI, Mateus J, Dan JM, Moderbacher CR, et al. Targets of T cell responses to SARS-CoV-2 coronavirus in humans with COVID-19 disease and unexposed individuals. Cell. 2020;S0092867420306103.

  33. 33.

    Abelin JG, Keskin DB, Sarkizova S, Hartigan CR, Zhang W, Sidney J, et al. Mass spectrometry profiling of HLA-associated peptidomes in mono-allelic cells enables more accurate epitope prediction. Immunity. 2017;46(2):315–26.

  34. 34.

    Abelin JG, Harjanto D, Malloy M, Suri P, Colson T, Goulding SP, et al. Defining HLA-II ligand processing and binding rules with mass spectrometry enhances cancer epitope prediction. Immunity. 2019;51(4):766–779.e17.

  35. 35.

    Archila LLD, Kwok WW. Tetramer-guided epitope mapping: a rapid approach to identify HLA-restricted T-cell epitopes from composite allergens. Methods Mol Biol. 2017;1592:199–209.

    CAS  PubMed  Google Scholar 

  36. 36.

    Yang J, James EA, Huston L, Danke NA, Liu AW, Kwok WW. Multiplex mapping of CD4 T cell epitopes using class II tetramers. Clin Immunol. 2006;120(1):21–32.

    CAS  PubMed  Google Scholar 

  37. 37.

    Vita R, Mahajan S, Overton JA, Dhanda SK, Martini S, Cantrell JR, et al. The Immune Epitope Database (IEDB): 2018 update. Nucleic Acids Res. 2019 08;47(D1):D339–43.

  38. 38.

    Grifoni A, Sidney J, Zhang Y, Scheuermann RH, Peters B, Sette A. A sequence homology and Bioinformatic approach can predict candidate targets for immune responses to SARS-CoV-2. Cell Host Microbe. 2020;12.

  39. 39.

    Pickett BE, Sadat EL, Zhang Y, Noronha JM, Squires RB, Hunt V, et al. ViPR: an open bioinformatics database and analysis resource for virology research. Nucleic Acids Res. 2012;40(Database issue):D593–8.

    CAS  PubMed  Google Scholar 

  40. 40.

    Bojkova D, Klann K, Koch B, Widera M, Krause D, Ciesek S, et al. Proteomics of SARS-CoV-2-infected host cells reveals therapy targets. Nature. 2020;583(7816):469–72.

    CAS  PubMed  Google Scholar 

  41. 41.

    Bezstarosti K, Lamers MM, Haagmans BL, Demmers JAA. Targeted proteomics for the detection of SARS-CoV-2 proteins. bioRxiv. 2020;2020.04.23.057810.

  42. 42.

    Davidson, A.D., Williamson, M.K., Lewis, S. et al. Characterisation of the transcriptome and proteome of SARS-CoV-2 reveals a cell passage induced in-frame deletion of the furin-like cleavage site from the spike glycoprotein. Genome Med 12, 68 (2020). https://doi.org/10.1186/s13073-020-00763-0.

  43. 43.

    Maiers M, Gragert L, Klitz W. High-resolution HLA alleles and haplotypes in the United States population. Hum Immunol. 2007;68(9):779–88.

    CAS  PubMed  Google Scholar 

  44. 44.

    U.S. Census Bureau QuickFacts: United States. [cited 2020 Jun 24]. Available from: https://www.census.gov/quickfacts/fact/table/US/PST045219.

  45. 45.

    González-Galarza FF, Takeshita LYC, Santos EJM, Kempson F, Maia MHT, da Silva ALS, et al. Allele frequency net 2015 update: new features for HLA epitopes, KIR and disease and HLA adverse drug reaction associations. Nucleic Acids Res. 2015;43(Database issue):D784–8.

    PubMed  Google Scholar 

  46. 46.

    Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12(6):996–1006.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Hadrup SR, Bakker AH, Shu CJ, Andersen RS, van Veluw J, Hombrink P, et al. Parallel detection of antigen-specific T-cell responses by multidimensional encoding of MHC multimers. Nat Methods. 2009;6(7):520–6.

    CAS  PubMed  Google Scholar 

  48. 48.

    Toebes M, Coccoris M, Bins A, Rodenko B, Gomez R, Nieuwkoop NJ, et al. Design and use of conditional MHC class I ligands. Nat Med. 2006;12(2):246–51.

    CAS  PubMed  Google Scholar 

  49. 49.

    Zhao J, Zhao J, Mangalam AK, Channappanavar R, Fett C, Meyerholz DK, et al. Airway memory CD4+ T cells mediate protective immunity against emerging respiratory coronaviruses. Immunity. 2016;44(6):1379–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Xing Y, Hogquist KA. T-cell tolerance: central and peripheral. Cold Spring Harb Perspect Biol. 2012;4(6):a006957.

    PubMed  PubMed Central  Google Scholar 

  51. 51.

    Parkhurst MR, Robbins PF, Tran E, Prickett TD, Gartner JJ, Jia L, et al. Unique neoantigens arise from somatic mutations in patients with gastrointestinal cancers. Cancer Discov. 2019;9(8):1022–35.

    PubMed  PubMed Central  Google Scholar 

  52. 52.

    Yadav M, Jhunjhunwala S, Phung QT, Lupardus P, Tanguay J, Bumbaca S, et al. Predicting immunogenic tumour mutations by combining mass spectrometry and exome sequencing. Nature. 2014;515(7528):572–6.

    CAS  PubMed  Google Scholar 

  53. 53.

    Stronen E, Toebes M, Kelderman S, van Buuren MM, Yang W, van Rooij N, et al. Targeting of cancer neoantigens with donor-derived T cell receptor repertoires. Science. 2016;352(6291):1337–41.

    CAS  PubMed  Google Scholar 

  54. 54.

    Ali M, Foldvari Z, Giannakopoulou E, Böschen M-L, Strønen E, Yang W, et al. Induction of neoantigen-reactive T cells from healthy donors. Nat Protoc. 2019;14(6):1926–43.

    CAS  PubMed  Google Scholar 

  55. 55.

    Liu H, Sadygov RG, Yates JR. A model for random sampling and estimation of relative protein abundance in shotgun proteomics. Anal Chem. 2004;76(14):4193–201.

    CAS  PubMed  Google Scholar 

  56. 56.

    Zybailov B, Mosley AL, Sardiu ME, Coleman MK, Florens L, Washburn MP. Statistical analysis of membrane proteome expression changes in Saccharomyces cerevisiae. J Proteome Res. 2006;5(9):2339–47.

    CAS  PubMed  Google Scholar 

  57. 57.

    Ihling C, Tänzler D, Hagemann S, Kehlen A, Hüttelmaier S, Arlt C, et al. Mass spectrometric identification of SARS-CoV-2 proteins from gargle solution samples of COVID-19 patients. J Proteome Res. 2020 [cited 2020 Jul 17]; Available from: https://doi.org/10.1021/acs.jproteome.0c00280.

  58. 58.

    Shi Y, Yi Y, Li P, Kuang T, Li L, Dong M, et al. Diagnosis of severe acute respiratory syndrome (SARS) by detection of SARS coronavirus Nucleocapsid antibodies in an antigen-capturing enzyme-linked Immunosorbent assay. J Clin Microbiol. 2003;41(12):5781–2.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Kim D, Lee J-Y, Yang J-S, Kim JW, Kim VN, Chang H. The architecture of SARS-CoV-2 transcriptome. Cell. 2020 [cited 2020 May 9]; Available from: http://www.sciencedirect.com/science/article/pii/S0092867420304062.

  60. 60.

    Shomuradova AS, Vagida MS, Sheetikov SA, Zornikova KV, Kiryukhin D, Titov A, et al. SARS-CoV-2 epitopes are recognized by a public and diverse repertoire of human T-cell receptors. 2020 May [cited 2020 Jun 24]. Available from: http://medrxiv.org/lookup/doi/10.1101/2020.05.20.20107813.

  61. 61.

    Chour W, Xu AM, Ng AHC, Choi J, Xie J, Yuan D, et al. Shared antigen-specific CD8+ T cell responses against the SARS-COV-2 spike protein in HLA A*02:01 COVID-19 participants. 2020 May [cited 2020 Jun 24]. Available from: http://medrxiv.org/lookup/doi/10.1101/2020.05.04.20085779.

  62. 62.

    Júnior IJM, Polveiro RC, Souza GM, Bortolin DI, Sassaki FT, Lima ATM. The global population of SARS-CoV-2 is composed of six major subtypes. bioRxiv. 2020;2020.04.14.040782.

  63. 63.

    Ahmed SF, Quadeer AA, McKay MR. COVIDep: a web-based platform for real-time reporting of vaccine target recommendations for SARS-CoV-2. Nat Protoc. 2020;17:1–2.

    Google Scholar 

  64. 64.

    Singer J, Gifford R, Cotten M, Robertson D. CoV-GLUE: a web application for tracking SARS-CoV-2 genomic variation. Life Sciences; 2020 Jun [cited 2020 Jun 24]. Available from: https://www.preprints.org/manuscript/202006.0225/v1.

  65. 65.

    Sarkizova S, Klaeger S, Le PM, Li LW, Oliveira G, Keshishian H, et al. A large peptidome dataset improves HLA class I epitope prediction across most of the human population. Nat Biotechnol. 2020;38(2):199–209.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors wish to thank the Research and Development department at BioNTech US (formerly Neon Therapeutics, Inc.) for the team efforts and the wide base of knowledge making this study possible. Specifically, the authors wish to thank Yvonne Ware and Jonathan McGee from the peptide chemistry team for their dedication and hard work on manufacturing the T cell assay peptides in record time. The authors wish to thank Prerna Suri for her assistance in cell culture work for the T cell assay. The authors would also like to thank John Welle for his contribution to graphically visualizing the SARS-CoV-2 virus.

Funding

BioNTech US (formerly Neon Therapeutics, Inc.)

Author information

Affiliations

Authors

Contributions

A.P.: conceptualization, formal analysis, investigation, visualization, writing—original draft; D.H.: conceptualization, formal analysis, investigation, visualization, writing—original draft; M.M.: conceptualization, formal analysis, investigation, visualization, writing—original draft; C.M.A: conceptualization, methodology, investigation, writing—original draft; D.A.R.: methodology, investigation, writing—original draft; D.L.: methodology; MMvB: methodology; T.A.A.: supervision, writing—review and editing; M.S.R: supervision, writing—review and editing; L.S.: supervision, conceptualization, writing—original draft; R.B.G: supervision, funding acquisition, writing—review and editing. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Asaf Poran or Dewi Harjanto.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

AP, DH, MM, CMA, DAR, DL, MMvB, TAA, MSR, LS, and RBG are all current or past employees and shareholders of BioNTech US (formerly Neon Therapeutics, Inc.).

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

HLA-I Alleles Covered by Binding Predictors. List of the 74 alleles covered by the HLA-I binding predictor. (CSV 896 bytes)

Additional file 2: Table S2.

HLA-II Alleles Covered by Binding Predictor. List of the 83 alleles covered by the HLA-II binding predictor. (CSV 1 kb)

Additional file 3: Table S3.

HLA-I and HLA-II Allele Population Frequencies. HLA-I and HLA-II allele frequencies for USA, European (EUR), Asian Pacific Islander (API), African American (AFA), and Hispanic (HIS) populations.

Additional file 4: Table S4.

Predicted HLA-I Binders Ranked by Population Coverage. Table of all predicted HLA-I binders and their associated allele coverage, including additional indicators for overlap with the human proteome or overlap with the ViPR dataset used. (CSV 1448 kb)

Additional file 5: Table S5.

Broadly Binding HLA-I Peptides. The top HLA-I predicted binders from each of the three SARS-CoV-2 proteins: spike, nucleocapsid and membrane with the broadest cumulative allele coverage. (CSV 19 kb)

Additional file 6: Table S6.

SARS-CoV-2 25mers Ranked by HLA-II Population Coverage. Table of all SARS-CoV-2-derived 25mers containing at least 3 predicted HLA-II binders as subsequences. (CSV 1268 kb)

Additional file 7: Table S7.

Broadly Binding HLA-II 25mers. The top HLA-II predicted binders from each of the three SARS-CoV-2 proteins: spike, nucleocapsid. (CSV 16 kb)

Additional file 8: Table S8.

binding prediction of ViPR HLA-I epitopes. The peptide-HLA alleles pairs from the ViPR database which belong to the Coronaviridae family and have a human host that had been scored using our HLA-I binding predictor. (CSV 522 kb)

Additional file 9: Table S9.

binding prediction of ViPR HLA-II epitopes. The peptide-HLA alleles pairs from the ViPR database which belong to the Coronaviridae family and have a human host that had been scored using our HLA-II binding predictor. (CSV 39 kb)

Additional file 10: Table S10.

spectral counts from published SARS-CoV-2 proteomic datasets. MS/MS spectra assigned to peptides from SARS-CoV-2 proteins were tallied across datasets, divided by protein length, and normalized within each dataset to generate Fig. 5.

Additional file 11.

Custom Python script for HLA-I. This Python scripts can be used to generate HLA-I supplementary tables. (PY 5 kb)

Additional file 12.

Custom Python script for HLA-II. This Python scripts can be used to generate HLA-II related supplementary tables. (PY 8 kb)

Additional file 13.

Custom R script for figure plotting. This R script can be used to plot the precision-recall analyses from Fig. 1, reproduce Fig. 5, and to produce Table 1. (R 9 kb)

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Poran, A., Harjanto, D., Malloy, M. et al. Sequence-based prediction of SARS-CoV-2 vaccine targets using a mass spectrometry-based bioinformatics predictor identifies immunogenic T cell epitopes. Genome Med 12, 70 (2020). https://doi.org/10.1186/s13073-020-00767-w

Download citation

Keywords

  • COVID-19
  • SARS-CoV-2 T cell epitopes
  • Computational biology
  • HLA-I binding prediction
  • HLA-II binding prediction
  • T cell assay
  • Vaccine design