key: cord-294429-isivkz8b
authors: Grifoni, Alba; Weiskopf, Daniela; Ramirez, Sydney I.; Mateus, Jose; Dan, Jennifer M.; Moderbacher, Carolyn Rydyznski; Rawlings, Stephen A.; Sutherland, Aaron; Premkumar, Lakshmanane; Jadi, Ramesh S.; Marrama, Daniel; de Silva, Aravinda M.; Frazier, April; Carlin, Aaron; Greenbaum, Jason A.; Peters, Bjoern; Krammer, Florian; Smith, Davey M.; Crotty, Shane; Sette, Alessandro
title: Targets of T cell responses to SARS-CoV-2 coronavirus in humans with COVID-19 disease and unexposed individuals
date: 2020-05-20
journal: Cell
DOI: 10.1016/j.cell.2020.05.015
sha: 
doc_id: 294429
cord_uid: isivkz8b

Summary Understanding adaptive immunity to SARS-CoV-2 is important for vaccine development, interpreting coronavirus disease 2019 (COVID-19) pathogenesis, and calibration of pandemic control measures. Using HLA class I and II predicted peptide ‘megapools’, circulating SARS-CoV-2−specific CD8+ and CD4+ T cells were identified in ∼70% and 100% of COVID-19 convalescent patients, respectively. CD4+ T cell responses to spike, the main target of most vaccine efforts, were robust and correlated with the magnitude of the anti-SARS-CoV-2 IgG and IgA titers. The M, spike and N proteins each accounted for 11-27% of the total CD4+ response, with additional responses commonly targeting nsp3, nsp4, ORF3a and ORF8, among others. For CD8+ T cells, spike and M were recognized, with at least eight SARS-CoV-2 ORFs targeted. Importantly, we detected SARS-CoV-2−reactive CD4+ T cells in ∼40-60% of unexposed individuals, suggesting cross-reactive T cell recognition between circulating ‘common cold’ coronaviruses and SARS-CoV-2.

COVID-19 is a world-wide emergency. The first cases occurred in December 2019 and now more than 240,000 deaths and 3,000,000 cases of SARS-CoV-2 infection have been reported globally as of May 1 st , (Dong et al., 2020; Wu and McGoogan, 2020) . Vaccines against SARS-CoV-2 are just beginning development (Amanat and Krammer, 2020; Thanh Le et al., 2020) . An understanding of human T cell responses to SARS-CoV2 is lacking, due to the rapid emergence of the pandemic. There is an urgent need for foundational information about T cell responses to this virus.

The first steps for such an understanding are the ability to quantify the virus-specific CD4 + and CD8 + T cells. Such knowledge is of immediate relevance, as it will provide insights into immunity and pathogenesis of SARS-CoV-2 infection, and the same knowledge will assist vaccine design and evaluation of candidate vaccines. Estimations of immunity are also central to epidemiological model calibration of future social distancing pandemic control measures (Kissler et al., 2020) . Such projections are dramatically different depending on whether SARS-CoV-2 infection creates substantial immunity, and whether any crossreactive immunity exists between SARS-CoV-2 and circulating seasonal 'common cold' human coronaviruses. Definition and assessment of human antigen-specific SARS-CoV2 T cell responses are best made with direct ex vivo T cell assays using broad-based epitope pools and assays capable of detecting T cells of any cytokine polarization. Herein, we have completed such an assessment with blood samples from COVID-19 patients.

There is also great uncertainty about whether adaptive immune responses to SARS-CoV-2 are protective or pathogenic, or whether both scenarios can occur depending on timing, composition, or magnitude of the adaptive immune response. Hypotheses range the full gamut (Peeples, 2020) , based on available clinical data from severe acute respiratory disease syndrome (SARS) or MERS (Alshukairi et al., 2018; Wong et al., 2004; Zhao et al., 2017) or animal model data with SARS in mice (Zhao et al., 2016; Zhao et al., 2010; Zhao et al., 2009) , SARS in NHPs (Liu et al., 2019; Takano et al., 2008) or FIPV in cats (Vennema et al., 1990) . Protective immunity, immunopathogenesis, and vaccine development for COVID-19 are each briefly discussed below, related to introducing the importance of defining T cell responses to SARS-CoV-2.

Based on data from SARS patients in [2003] [2004] (caused by SARS-CoV, the most closely related human betacoronavirus to SARS-CoV-2), and based on the fact that most acute viral infections result in development of protective immunity (Sallusto et al., 2010) , a likely possibility has been that substantial CD4 + T cell, CD8 + T cell, and neutralizing antibody responses develop to SARS-CoV-2 and all contribute to clearance of the acute infection; and, as a corollary, some of the T and B cells are retained long-term (i.e., multiple years) as immunological memory and protective immunity against SARS-CoV-2 infection (Guo et al., 2020b; Li et al., 2008) . However, a contrarian viewpoint is also legitimate. While most acute infections result in the development of protective immunity, available data for human coronaviruses suggest the possibility that substantive adaptive immune responses can fail to occur (Choe et al., 2017; Okba et al., 2019; Zhao et al., 2017) and robust protective immunity can fail to develop (Callow et al., 1990) . A failure to develop protective immunity could occur due to a T cell and/or antibody response of insufficient magnitude or durability, with the neutralizing antibody response being dependent on the CD4 + T cell response (Crotty, 2019; Zhao et al., 2016) . Thus, there is urgent need to understand the magnitude and composition of the human CD4 + and CD8 + T cell responses to SARS-CoV-2. If natural infection with SARS-CoV-2 elicits potent CD4 + and CD8 + T cell responses commonly associated with protective antiviral immunity, COVID-19 is a strong candidate for rapid vaccine development.

Immunopathogenesis in COVID-19 is a serious concern (Cao, 2020; Peeples, 2020) . It is most likely that an early CD4 + and CD8 + T cell response against SARS-CoV-2 is protective, but an early response is difficult to generate because of efficient innate immune evasion mechanisms of SARS-CoV-2 in humans (Blanco-Melo et al., 2020) . Immune evasion by SARS-CoV-2 is likely exacerbated by reduced myeloid cell antigen presenting cell (APC) function or availability in the elderly (Zhao et al., 2011) . In such cases, it is conceivable that late T cell responses may instead amplify pathogenic inflammatory outcomes in the presence of sustained high viral loads in the lungs, by multiple hypothetical possible mechanisms (Guo et al., 2020a; Li et al., 2008; Liu et al., 2019) . Critical (ICU) and fatal outcomes are associated with elevated levels of inflammatory cytokines and chemokines, including IL-6 (Giamarellos- Bourboulis et al., 2020; Wong et al., 2004; Zhou et al., 2020) Vaccine development against acute viral infections classically focuses on vaccine-elicited recapitulation of the type of protective immune response elicited by natural infection. Such foundational knowledge is currently missing for COVID-19, including how the balance and the phenotypes of responding cells vary as a function of disease course and severity. Such knowledge can guide selection of vaccine strategies most likely to elicit protective immunity against SARS-CoV-2. Furthermore, knowledge of the T cell responses to COVID-19 can guide selection of appropriate immunological endpoints for COVID-19 candidate vaccine clinical trials, which are already starting.

Limited information is also available about which SARS-CoV-2 proteins are recognized by human T cell immune responses. In some infections, T cell responses are strongly biased towards certain viral proteins, and the targets can vary substantially between CD4 + and CD8 + T cells (Moutaftsi et al., 2010; Tian et al., 2019) . Knowledge of SARS-CoV-2 proteins and epitopes recognized by human T cell responses is of immediate relevance, as it will allow for monitoring of COVID-19 immune responses in laboratories worldwide. Epitope knowledge will also assist candidate vaccine design and facilitate evaluation of vaccine candidate immunogenicity. Almost all of the current COVID-19 vaccine candidates are focused on the spike protein.

A final key issue to consider in the study of SARS-CoV-2 immunity is whether some degree of cross-reactive coronavirus immunity exists in a fraction of the human population, and whether this might influence susceptibility to COVID-19 disease. This issue is also relevant for vaccine development, as cross-reactive immunity could influence responsiveness to candidate vaccines (Andrews et al., 2015) .

In sum, the ability to measure and understand the human CD4 + and CD8 + T cell responses to SARS-CoV-2 infection is a major knowledge gap currently impeding COVID-19 vaccine development, interpretation of COVID-19 disease pathogenesis, and calibration of future social distancing pandemic control measures.

We recently predicted SARS-CoV-2 T cell epitopes utilizing the Immune Epitope Database and Analysis Resource (IEDB) Vita et al., 2019) . Utilizing bioinformatic approaches, we identified specific peptides in SARS-CoV-2 with increased probability of being T cell targets (Grifoni et al., 2020) . We previously developed the megapool (MP) approach to allow simultaneous testing of large numbers of epitopes. By this technique, numerous epitopes are solubilized, pooled and relyophilized to avoid cell toxicity problems . These MPs have been used in human T cell studies of a number of indications, including allergies (Hinz et al., 2016) , tuberculosis , tetanus, pertussis (Bancroft et al., 2016; da Silva Antunes et al., 2017) and dengue virus, for both CD4 + and CD8 + T cell epitopes (Grifoni et al., 2017; Weiskopf et al., 2015) . Here, we generated MPs based on predicted SARS-CoV-2 epitopes. Specifically, one MP corresponds to 221 predicted HLA class II CD4 + T cell epitopes (Grifoni et al., 2020) covering all proteins in the viral genome, apart from the spike (S) antigen (CD4_R MP). The prediction strategy utilized is geared to capture ~50% of the total response (Dhanda et al., 2018; Paul et al., 2015) and was designed and validated to predict dominant epitopes independently of ethnicity and HLA polymorphism. This approach takes advantage of the extensive cross-reactivity and repertoire overlap between different HLA class II loci and allelic variants to predict promiscuous epitopes, capable of binding many of the most common HLA class II prototypic specificities (Greenbaum et al., 2011; O'Sullivan et al., 1991; Sidney et al., 2010a, b; Southwood et al., 1998) .

For the spike protein, to ensure that all T cell reactivity against this important antigen can be detected, we generated a separate MP covering the entire antigen with 253 15-mer peptides overlapping by 10-residues (MP_S, Table S1). As stated above, the MP used to probe the non-spike regions is expected to capture ~50% of the total response. The use of overlapping peptides spanning entire ORFs instead allows for a more complete characterization, but also requires more cells. This factor should be kept in mind in terms of comparison of the magnitude of the CD4 + T cell responses to those pools.

In the case of CD8 epitopes, since the overlap between different HLA class I allelic variants and loci is more limited to specific groups of alleles, or supertypes (Sidney et al., 2008) , we targeted a set of the 12 most prominent HLA class I A and B alleles, which together allow broad coverage (>85%) of the general population. Two class I MPs were synthesized based on epitope predictions for those 12 most common HLA A and B alleles (Grifoni et al., 2020) , which collectively encompass 628 predicted HLA class I CD8 + T cell epitopes from the entire SARS-CoV-2 proteome (CD8 MP-A and MP-B).

To test for the generation of SARS-CoV-2 CD4 + and CD8 + T cell responses following infection, we initially recruited 20 adult patients who had recovered from COVID-19 disease ( Table 1) . We also utilized PBMC and plasma samples from local healthy control donors collected in 2015-2018 (see STAR Methods). Blood samples were collected at 20-35 days post-symptoms onset from nonhospitalized COVID-19 patients who were no longer symptomatic. SARS-CoV-2 infection was determined by swab test viral PCR during the acute phase of the infection. Verification of SARS-CoV-2 exposure was attempted both by lateral flow serology and SARS-CoV-2 spike protein receptor binding domain (RBD) ELISA (Stadlbauer et al., 2020) , using plasma from the convalescence stage blood draw. Most patients were confirmed positive by lateral flow Ig tests ( Table 1) . All patients were confirmed COVID-19 cases by SARS-CoV-2 RBD ELISA (Fig. 1, S1 ). All cases were IgG positive; anti-RBD IgM and IgA was also detected in the large majority of cases ( Fig. 1, S1) .

We defined a 21-color flow cytometry panel of mononuclear leukocyte lineage and phenotypic markers (Table S2) to broadly assess the immunological cellular profile of recovered COVID-19 patients (Fig. 1, S2) . The frequency of CD3 + cells was slightly increased in recovered COVID-19 patients relative to non-exposed controls, while no significant differences overall were observed in the frequencies of CD4 + or CD8 + T cells between the two groups. Frequencies of CD19 + cells were somewhat decreased, while no differences were observed in the frequencies of CD3 -CD19cells or CD14 + CD16monocytes ( Fig. 1, S2) . No evidence of general lymphopenia was observed in the convalescing patients, consistent with the literature. Next, we utilized the SARS-CoV-2 MPs to probe CD4 + and CD8 + T cell responses.

We utilized T cell receptor− (TCR) dependent Activation Induced Marker (AIM) assays to identify and quantify SARS-CoV-2−specific CD4 + T cells in recovered COVID-19 patients. Initial definition and assessment of human antigen-specific SARS-CoV-2 T cell responses are best made with direct ex vivo T cell assays using broad-based epitope pools, such as MPs, and assays capable of detecting T cells of unknown cytokine polarization and functional attributes. AIM assays are cytokine-independent assays to identify antigen-specific CD4 + T cells Reiss et al., 2017) . AIM assays have been successfully used to identify virus-specific, vaccine-specific, or tuberculosisspecific CD4 + T cells in a range of studies (Dan et al., 2019; Dan et al., 2016; Herati et al., 2017; Morou et al., 2019) . We stimulated PBMCs from 10 COVID-19 cases and 11 healthy controls (SARS-CoV-2 unexposed, collected in 2015-2018) with a spike MP (MP_S) and the class II MP covering the remainder of the SARS-CoV-2 orfeome ("non-spike", MP CD4_R). A CMV MP was used as a positive control, while DMSO was used as the negative control (Fig. 2, S3) . SARS-CoV-2 spike-specific CD4 + T cell responses (OX40 + CD137 + ) were detected in 100% of COVID-19 cases (p < 0.0001 vs. unexposed donors spike MP, Fig. 2A -B. p = 0.002 vs. DMSO control, Fig. 2C) . CD4 + T cell responses to the remainder of the SARS-CoV-2 orfeome were also detected in 100% of COVID-19 cases (p < 0.0079 vs. unexposed donors non-spike MP, Fig. 2A -B. p = 0.002, non-spike vs. DMSO control, Fig. 2C ). The magnitude of the SARS-CoV-2−specific CD4 + T cell responses measured was similar to that of the CMV MP (Fig. S3C) . The concordance between SARS-CoV-2−specific CD4 + T cell measurements in independent experiments was high (p < 0.0002, Fig S3D) . To assess functionality and polarization of the SARS-CoV-2−specific CD4 + T cell response, we measured cytokines secreted in response to MP stimulation. The SARS-CoV-2−specific CD4 + T cells were functional, as the cells produced IL-2 in response to non-spike and spike MPs (Fig. 2D) . Polarization of the cells appeared to be a classical T H 1 type, as substantial IFNγ was produced (Fig. 2E) , while little to no IL-4, Il-5, IL-13, or IL-17α was expressed ( Fig. S3G-J) .

Thus, recovered COVID-19 patients consistently generated a substantial CD4 + T cell response against SARS-CoV-2. Similar conclusions were reached using stimulation index as the metric ( Fig.  S3E-F) . In terms of total CD4 + T cell response per donor ( Fig. 2A) , on average ~50% of the detected response was directed against the spike protein, and ~50% was directed against the MP representing the remainder of the SARS-CoV-2 orfeome ( Fig. 2A) . This is of significance, since the SARS-CoV-2 spike protein is a key component of the vast majority of candidate COVID-19 vaccines under development. Of note, given the nature of the MP_R peptide predictions, the actual CD4 + T cell response to be ascribed to non-spike ORFs was likely to be higher, addressed in further experiments below.

To measure SARS-CoV-2−specific CD8 + T cells in the recovered COVID-19 patients, we utilized two complementary methodologies, AIM assays and intracellular cytokine staining (ICS). The two SARS-CoV-2 class I MPs were used, CD8-A and CD8-B, with CMV MP and DMSO serving as positive and negative controls, respectively (Fig. 3, S4) . CD8 + T cell responses were detected by AIM (CD69 + CD137 + ) in 70% of COVID-19 cases (p < 0.0011 vs. unexposed donors "CD8 total", Fig. 3A -B. p = 0.002, CD8-A or CD8-B vs. DMSO control, Fig. S4B ). MP CD8-A contains spike epitopes, among epitopes to other proteins. The magnitude of the SARS-CoV-2 reactive CD8 + T cell responses measured by AIM was somewhat lower than the CMV MP (Fig. S4C) . Similar conclusions were reached using stimulation index ( Fig. S3D-E) .

Independently, ICS assays detected IFNγ + SARS-CoV-2−specific CD8 + T cells in the majority of COVID-19 cases ( Fig. 3C-D) . The majority of IFNγ + cells co-expressed granzyme B (Fig. 3D-E) . A substantial fraction of the IFNγ + cells expressed TNF, but not IL-10 (Fig. 3D) . Thus, the majority of recovered COVID-19 patients generated a CD8 + T cell response against SARS-CoV-2.

Relationship between SARS-CoV-2− − − −specific CD4 + T cell responses and IgG and IgA titers Most protective antibody responses are dependent on CD4 + T cell help. Therefore, we assessed whether stronger SARS-CoV-2−specific CD4 + T cell responses were associated with higher antibody titers in COVID-19 cases. Given that spike is the primary target of SARS neutralizing antibodies, we examined spike-specific CD4 + T cells. Spike-specific CD4 + T cell responses correlated well with the magnitude of the anti-spike RBD IgG titers (R = 0.81. p < 0.0001. Fig. 4A . Similar results were obtained using stimulation index, Fig. S5A ). The non-spike SARS-CoV-2−specific CD4 + T cell response did not correlate as well with anti-spike RBD IgG titers (Fig. 4B, S5B ), consistent with a common requirement for intramolecular CD4 + T cell help . Anti-spike IgA titers also correlated with spikespecific CD4 + T cells (p < 0.0002, Fig. S5 ). Thus, COVID-19 patients make anti-spike RBD antibody responses commensurate with the magnitude of their spike-specific CD4 + T cell response. We then assessed the relationship between the CD4 + and CD8 + T cell responses to SARS-CoV-2. SARS-CoV-2−specific CD4 + and CD8 + T cell responses were well correlated (R=0.62. p = 0.0025, Fig. 4C and S5) . Thus, antibody, CD4 + and CD8 + T cell responses to SARS-CoV-2 were generally well correlated.

While spike− and non-spike−specific CD4 + T cell responses were detectable in all COVID-19 cases, cells were also detected in unexposed individuals ( Fig. 3A-B) . These responses were statistically significant for non-spike-specific CD4 + T cell reactivity (non-spike, p = 0.039. Spike, p = 0.067. Fig. 5A -B). Non-spike−specific CD4 + T cell responses were above the limit of detection in 50% of donors based on SI (Fig. S3E) . All of the donors were recruited between 2015-2018, excluding any possibility of exposure to SARS-CoV-2. Four human coronaviruses are known causes of seasonal 'common cold' upper respiratory tract infections: HCoV-OC43, HCoV-HKU1, HCoV-NL63, and HCoV-229E. We tested the SARS-CoV-2 unexposed donors for seroreactivity to HCoV-OC43 and HCoV-NL63 as a representative betacoronavirus and alphacoronavirus, respectively. All donors were IgG seropositive to HCoV-OC43 and HCoV-NL63 RBD, to varying degrees ( Fig. 5C) , consistent with the endemic nature of these viruses (Gorse et al., 2010; Huang et al., 2020; Severance et al., 2008) . We therefore examined whether these represented true pan-coronavirus T cells capable of recognizing SARS-CoV-2 epitopes.

A most pressing, yet unresolved, set of issues in understanding SARS-CoV-2 immune responses is what antigens are targeted by CD4 + and CD8 + T cells, whether the corresponding antigens are the same or different, and how do they reflect the antigens currently considered for COVID-19 vaccine development. We synthesized sets of overlapping peptides spanning the entire sequence of SARS-CoV-2, and pooled them separately so that each pool would represent one antigen (with the exception of nsp3, for which two pools were made. Table S1 ).

In the case of CD4 + T cell responses, no obvious pattern of antigen specificity was observed based on SARS-CoV-2 genome organization; however, coronaviruses increase protein synthesis of certain ORFs in infected cells via subgenomic RNAs. Accounting for the relative abundance of subgenomic RNAs (Fig. 6A ) (Irigoyen et al., 2016; Snijder et al., 2003; Xie et al., 2020) , the ORFs were re-ordered based on predicted protein abundance (Fig. 6B) . A clear hierarchy of SARS-CoV-2−specific CD4 + T cell targets was then apparent, with the majority of the CD4 + T cell response in COVID-19 cases directed against highly expressed SARS-CoV-2 ORFs spike, M, and N. On average, these antigens accounted for 27%, 21% and 11% of the total CD4 + T cell response, respectively. Most COVID-19 cases also had CD4 + T cells specific for SARS-CoV-2 nsp3, nsp4 and ORF8 (Fig. 6B) , on average each accounting for ~5% of the total CD4 + T cell response (Fig. 6C) . E, ORF6, hypothetical ORF10, and nsp1 are all small antigens (or potentially not expressed, in the case of ORF10) and were most likely predominantly unrecognized as a result. These results are somewhat unexpected, because data for other coronaviruses, from 27 different studies curated in the IEDB, reported that spike accounted for nearly two-thirds of reported CD4 + T cell reactivity (Table S3) . N accounted for most of the remaining epitopes in the published literature, although human N-specific CD4 + T cell responses were not observed in one of the most comprehensive studies of human SARS-CoV-1 T cell responses (Li et al., 2008) . Coronavirus M has not previously been described as a prominent target of CD4 + T cell responses (Table S3 ). In sum, these results, fully scanning the SARS2 orfeome, demonstrate a pattern of robust and diverse SARS-CoV-2−specific CD4 + T cell reactivity in convalescing COVID-19 cases that correlated largely with predicted viral protein abundance in infected cells.

When examining the non-exposed donors, the pattern of CD4 + T cell targets changed. While S was still a relatively prominent target (23% of total, on average), there was no, or marginal, reactivity against SARS-CoV-2 N and M. Among donors with detectable CD4 + T cells, a shift in reactivity was observed towards SARS-CoV-2 nsp14 (25%), nsp4 (15%) and nsp6 (14%) (Fig 6B-C) . SARS-CoV-2−reactive CD4 + T cells were detected in at least six different unexposed donors, demonstrating that the crossreactivity is relatively widely distributed (Fig. S6A) .

Having scanned the full SARS-CoV-2 orfeome for CD4 + T cell reactivity in multiple donors, it was possible to assess whether the epitope prediction MP approach successfully enriched for SARS-CoV-2 epitopes targeted by human CD4 + T cells. When the total reactivity observed with the CD4_R MP was plotted versus the sum total of all antigen pools (excluding spike, given that spike predictions were not included in the CD4_R MP), a significant correlation was observed (p < 0.0002, Fig S6C) . The single MP-R captured ~50% (44% +/-range 28-80%) of the non-spike response per COVID-19 donor, demonstrating the success of the prediction approach, which, as mentioned above, was devised to attempt to capture approximately 50% of the total response (Dhanda et al., 2018; Paul et al., 2015) .

In the case of CD8 + T cell responses, the data in the literature from other coronaviruses (57 different studies curated in the IEDB, Table S3 ) reported spike accounting for 50% and N accounting for 36% of the defined epitopes. In a large study of human SARS-CoV-1 responses, spike was reported as essentially the only target of CD8 + T cell responses (Li et al., 2008) , while in a study of MERS CD8 + T cells, responses were noted for spike, N and a pool of M/E peptides (Zhao et al., 2017) . Few epitopes have been reported from other coronavirus antigens (Table S3) . Here, we scanned the full SARS-CoV-2 orfeome for CD8 + T cell recognition. Our data indicate a somewhat different pattern of immunodominance for SARS-CoV-2 CD8 + T cell reactivity ( Fig. 6D-E) , with spike protein accounting for ~26% of the reactivity, and N accounting for ~12%. Significant reactivity in COVID-19 recovered subjects was derived from other antigens, such as M (22%), nsp6 (15%), ORF8 (10%) and ORF3a (7%) (Fig. 6D-E) . In unexposed donors, SARS-CoV-2−reactive CD8 + T cells were detected in at least four different donors (Fig. S7) , with less clear targeting of specific SARS-CoV-2 proteins than was observed for CD4 + T cells, suggesting that coronavirus CD8 + T cell crossreactivity exists but is less widespread than CD4 + T cell crossreactivity.

There is a critical need for foundational knowledge about T cell responses to SARS-CoV-2. Here we report functional validation of predicted epitopes when arranged in epitope "megapools", utilizing PBMCs derived from convalescing COVID-19 cases. The experiments also used protein-specific peptide pools to determine which SARS-CoV-2 proteins are the predominant targets of human SARS-CoV-2−specific CD4 + and CD8 + T cells generated during COVID-19 disease. Importantly, we utilized the exact same series of experimental techniques with blood samples from healthy control donors (PBMCs collected in the 2015-2018 time frame), and substantial crossreactive coronavirus T cell memory was observed.

Our results demonstrate that the epitope MPs are reagents well suited to analyze and detect SARS-CoV-2-specific T cell responses with limited sample material. We also developed and tested peptide pools corresponding to each of the 25 proteins encoded in the SARS-CoV-2 genome. Data from both the epitope MPs and protein peptide pool experiments can be interpreted in the context of previously reported T cell response immunodominance patterns observed for other coronaviruses, particularly the SARS and MERS viruses, which have been studied in humans, HLA-transgenic mice, wild-type mice and other species. In the case of CD4 + T cell responses, data for other coronaviruses found that spike accounted for nearly two-thirds of reported CD4 + T cell reactivity, with N and M accounting for limited reactivity, and no reactivity in one large study of human SARS-CoV-1 responses (Li et al., 2008) . Our SARS-CoV-2 data reveal that the pattern of immunodominance in COVID-19 is different. In particular, M, spike and N proteins were clearly co-dominant, each recognized by 100% of COVID-19 cases studied here. Significant CD4 + T cell responses were also directed against nsp3, nsp4, ORF3s, ORF7a, nsp12 and ORF8. These data suggest that a candidate COVID-19 vaccine consisting only of SARS-CoV-2 spike would be capable of eliciting SARS-CoV-2−specific CD4 + T cell responses of similar representation to that of natural COVID-19 disease, but the data also indicate that there are many potential CD4 + T cell targets in SARS-CoV-2 and inclusion of additional SARS-CoV-2 structural antigens such as M and N would better mimic the natural SARS-CoV-2−specific CD4 + T cell response observed in mild to moderate COVID-19 disease.

Regarding SARS-CoV-2 CD8 + T cell responses, the pattern of immunodominance found here differed from the literature for other coronaviruses. However, stringent comparisons are not possible, as some earlier studies were not similarly comprehensive and did not utilize the same experimental strategy. The spike protein was a target of human SARS-CoV-2 CD8 + T cell responses, but it is not dominant. SARS-CoV-2 M was just as strongly recognized, and significant reactivity was noted for other antigens, mostly nsp6, ORF3a, and N, which comprised nearly 50% of the total CD8 + T cell response, on average. Thus, these data indicate that candidate COVID-19 vaccines endeavoring to elicit CD8 + T cell responses against the spike protein will be eliciting a relatively narrow CD8 + T cell response compared to the natural CD8 + T cell response observed in mild to moderate COVID-19 disease. An optimal vaccine CD8 + T cell response to SARS-CoV-2 might benefit from additional class I epitopes, such as the ones derived from the M, nsp6, ORF3a and/or N.

There have been concerns regarding vaccine enhancement of disease by certain candidate COVID-19 vaccine approaches, via antibody-dependent enhancement (ADE) or development of a T H 2 responses (Peeples, 2020) . Herein, we saw predominant T H 1 responses in convalescing COVID-19 cases, with little to no T H 2 cytokines. Clearly more studies are required, but the data here appear to predominantly represent a classical T H 1 response to SARS-CoV-2.

While it was important to identify antigen-specific T cell responses in COVID-19 cases, it is also of great interest to understand whether crossreactive immunity exists between coronaviruses to any degree. A key step in developing that understanding is to examine antigen-specific CD4 + and CD8 + T cells in COVID-19 cases and in unexposed healthy controls, utilizing the exact same antigens and series of experimental techniques. CD4 + T cell responses were detected in 40-60% of unexposed individuals. This may be reflective of some degree of crossreactive, preexisting immunity to SARS-CoV-2 in some, but not all, individuals. Whether this immunity is relevant in influencing clinical outcomes is unknown-and cannot be known without T cell measurements before and after SARS-CoV-2 infection of individuals-but it is tempting to speculate that the crossreactive CD4 + T cells may be of value in protective immunity, based on SARS mouse models (Zhao et al., 2016) . Clear identification of the crossreactive peptides, and their sequence homology relation to other coronaviruses, requires deconvolution of the positive peptide pools, which is not feasible with the cell numbers presently available, and time frame of the present study.

Regarding the value of crossreactive T cells, influenza (flu) immunology in relationship to pandemics may be instructive. In the context of the 2009 H1N1 influenza pandemic, preexisting T cell immunity existed in the adult population, which focused on the more conserved internal influenza viral proteins (Greenbaum et al., 2009) . The presence of crossreactive T cells was found to correlate with less severe disease (Sridhar et al., 2013; Wilkinson et al., 2012) . The frequent availability of crossreactive memory T cell responses might have been one factor contributing to the lesser severity of the H1N1 flu pandemic (Hancock et al., 2009) . Cross-reactive immunity to influenza strains has been modeled to be a critical influencer of susceptibility to newly emerging, potentially pandemic, influenza strains (Gostic et al., 2016) . Given the severity of the ongoing COVID-19 pandemic, it has been modeled that any degree of crossprotective coronavirus immunity in the population could have a very substantial impact on the overall course of the pandemic, and the dynamics of the epidemiology for years to come (Kissler et al., 2020) .

Caveats of this study include the sample size and the focus on non-hospitalized COVID-19 cases. Sample size was limited by expediency. The focus on non-hospitalized cases of COVID-19 is a strength, in that these donors had uncomplicated disease of moderate duration, and thus it was encouraging that substantial CD4 + T cell and antibody responses were detected in all cases, and CD8 + T cell responses in the majority of cases. Complementing these data with MP T cell data from acute patients and patients with complicated disease course will also be of clear value, as will studies on the longevity of SARS-CoV-2 immunological memory. Additionally, lack of detailed information on common cold history or matched blood samples pre-exposure to SARS-CoV-2 prevents conclusions regarding the abundance of crossreactive coronavirus T cells before exposure to SARS-CoV-2 and any potential protective efficacy of such cells. Finally, full epitope mapping in the future will add important detailed resolution of the human coronavirus-specific T cell responses.

In sum, we measured SARS-CoV-2−specific CD4 + and CD8 + T cells responses in COVID-19 cases. Using multiple experimental approaches, SARS-CoV-2−specific CD4 + T cell and antibody responses were observed in all COVID-19 cases, and CD8 + T cell responses were observed in most. Importantly, pre-existing SARS-CoV-2−crossreactive T cell responses were observed in healthy donors, indicating some potential for pre-existing immunity in the human population. ORF mapping of T cell specificities revealed valuable targets for incorporation in candidate vaccine development, and revealed distinct specificity patterns between COVID-19 cases and unexposed healthy controls. Figure S4 .

Total MP responses per donor were used in each case ("Non-spike" + "spike" (CD4_R + MP_S) for CD4 + T cells, CD8_A + CD8_B for CD8 + T cells). Statistical comparisons were performed using Spearman correlation. See also Figure S5 . (B) SARS-CoV-2 antigen-specific CD4 + T cells (AIM + , OX40 + CD137 + ) quantified by stimulation index, using a peptide pool for each viral protein (with two exceptions, see Table S1 ). COVID-19 cases (top, in blue. n=10) and unexposed donors (bottom, in white. n=10).

(D) SARS-CoV-2 antigen-specific CD4 + T cells (AIM + , OX40 + CD137 + ) quantified by stimulation index, using a peptide pool for each viral protein (with two exceptions, see Table S1 ). COVID-19 cases (top, in red. n=10) and unexposed donors (bottom, in grey. n=10). Table S6 . Representative gating of CD3 + T cells, CD19 + B cells, CD3 -CD19cells, CD4 + T cells, CD8 + T cells and CD14 + monocytes from donor PBMCs is shown. Briefly, mononuclear cells were gated out of all events followed by subsequent singlet gating. Live cells are gated as Zombie UV -. Cells were then gated as CD19-PE-Cy5 + , CD3-buv395 + or CD19 -CD3cells. T cells were further subdivided into either CD8-buv805 + or CD4-PerCPefluor710 + populations. CD3 -CD19cells were defined as CD56-PE-Dazzle bright NK cells, CD56 dim CD-16buv737 + NK cells or CD56monocytes. Monocytes were further classified on differential expression of CD14-bv510 and CD16. (C) CMV-specific CD8 + T cells as percentage of AIM + (CD69 + CD137 + ) CD8 + T cells after stimulation of PBMCs with CMV peptide pool. Data were background subtracted against DMSO negative control and are shown with geometric mean and geometric standard deviation. Samples were from unexposed donors ("Unexposed", n=11) and recovered COVID-19 patients ("COVID-19", n=10). (D-E) Stimulation index quantitation of AIM + (CD69 + CD137 + ) CD8 + T cells; the same samples as in Figure 2 and Fig S4C were analyzed. Statistical comparisons across cohorts were performed with the Mann-Whitney test, while paired sample comparisons were performed with the Wilcoxon test. **p<0.01; ***p<0.001. ns not significant. (E) Correlation between SARS-CoV-2−specific CD4 + T cells and SARS-CoV-2−specific CD8 + T cells, using stimulation index. Total MP responses per donor were used in each case ("Non-spike" + "spike" (CD4_R + MP_S) for CD4 + T cells, CD8-A + CD8-B for CD8 + T cells). Statistical comparisons were performed using Spearman correlation. Figure 6 . (A) The same data as Figure 6B , but with each unexposed donor color coded. (B) The same experiment as Figure 6B , but with SARS-CoV-2−specific CD4 + T cells measured as percentage of AIM + (OX40 + CD137 + ) CD4 + T cells, after background subtraction. COVID-19 cases (top, in blue. n=10) and unexposed donors (bottom, in white. n=10). (C) Correlation of SARS-CoV-2−specific CD4 + T cells detected using the epitope prediction approach (CD4_R MP) compared against the sum total of all antigen pools of overlapping peptides (excluding spike), run with samples from the same donors in two different experiment series. Dotted line indicates 1:1 concordance. Statistical comparison was performed using Spearman correlation. Figure 6 . (A) The same data as Figure 6D , but with each unexposed donor color coded. (B) The same experiment as Figure 6D , but with SARS-CoV-2−specific CD8 + T cells measured as percentage of AIM + (CD69 + CD137 + ) CD8 + T cells, after background subtraction. COVID-19 cases (top, in red. n=10) and unexposed donors (bottom, in grey. n=10).

Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Dr. Alessandro Sette (alex@lji.org).

Aliquots of synthesized sets of peptides utilized in this study will be made available upon request. There are restrictions to the availability of the peptide reagents due to cost and limited quantity.

The published article includes all data generated or analyzed during this study, and summarized in the accompanying tables, figures and Supplemental materials.

Healthy Unexposed Donors. Samples from healthy adult donors were obtained by the La Jolla Institute for Immunology (LJI) Clinical Core or provided by a commercial vendor (Carter Blood Care) for prior, unrelated studies between early 2015 and early 2018. These samples were considered to be from unexposed controls, given that SARS-CoV-2 emerged as a novel pathogen in late 2019, more than one year after the collection of any of these samples. These donors were considered healthy in that they had no known history of any significant systemic diseases, including, but not limited to, autoimmune disease, diabetes, kidney or liver disease, congestive heart failure, malignancy, coagulopathy, hepatitis B or C, or HIV. An overview of the characteristics of these unexposed donors is provided in Table 1 .

The LJI Institutional Review Board approved the collection of these samples (LJI; VD-112). At the time of enrollment in the initial studies, all individual donors provided informed consent that their samples could be used for future studies, including this study.

Convalescent COVID-19 Donors. The Institutional Review Boards of the University of California, San Diego (UCSD; 200236X) and La Jolla Institute (LJI; VD-214) approved blood draw protocols for convalescent donors. All human subjects were assessed for capacity using a standardized and approved assessment. Subjects deemed to have capacity voluntarily gave informed consent prior to being enrolled in the study. Individuals did not receive compensation for their participation.

Study inclusion criteria included subjects over the age of 18 years, regardless of disease severity, race, ethnicity, gender, pregnancy or nursing status, who were willing and able to provide informed consent, or with a legal guardian or representative willing and able to provide informed consent when the participant could not personally do so. Study exclusion criteria included lack of willingness or ability to provide informed consent, or lack of an appropriate legal guardian or representative to provide informed consent.

Blood from convalescent donors was obtained at a UC San Diego Health clinic. Blood was collected in acid citrate dextrose (ACD) tubes and stored at room temperature prior to processing for PBMC isolation and plasma collection. A separate serum separator tube (SST) was collected from each donor. Samples were de-identified prior to analysis. Other efforts to maintain the confidentiality of participants included referring to specimens and other records via an assigned, coded identification number.

Prior to enrollment in the study, donors were asked to provide proof of positive testing for SARS-CoV-2, and screened for clinical history and/or epidemiological risk factors consistent with the World Health Organization (WHO) or Centers for Disease Control and Prevention (CDC) case definitions of COVID-19 or Persons Under Investigation (PUI) (https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/surveillanceand-case-definitions, https://www.cdc.gov/coronavirus/2019-nCoV/hcp/clinical-criteria.html). Per CDC and WHO guidance, clinical features consistent with COVID-19 included subjective or measured fever, signs or symptoms of lower respiratory tract illness (e.g. cough or dyspnea). Epidemiologic risk factors included close contact with a laboratory-confirmed case of SARS-CoV-2 within 14 days of symptom onset or a history of travel to an area with a high rate of COVID-19 cases within 14 days of symptom onset.

Disease severity was defined as mild, moderate, severe or critical based on a modified version of the WHO interim guidance, "Clinical management of severe acute respiratory infection when COVID-19 is suspected" (WHO Reference Number: WHO/2019-nCoV/clinical/2020.4). Mild disease was defined as an uncomplicated upper respiratory tract infection (URI) with potential non-specific symptoms (e.g. fatigue, fever, cough with or without sputum production, anorexia, malaise, myalgia, sore throat, dyspnea, nasal congestion, headache; rarely diarrhea, nausea and vomiting) that did not require hospitalization. Moderate disease was defined as the presence of lower respiratory tract disease or pneumonia without the need for supplemental oxygen, without signs of severe pneumonia, or a URI requiring hospitalization (including observation admission status). Severe disease was defined as severe lower respiratory tract infection or pneumonia with fever plus any one of the following: tachypnea (respiratory rate > 30 breaths per minute), respiratory distress, or oxygen saturation less than 93% on room air. Critical disease was defined as the need for ICU admission or the presence of acute respiratory distress syndrome (ARDS), sepsis, or septic shock, as defined in the WHO guidance document.

Convalescent donors were screened for symptoms prior to scheduling blood draws, and had to be symptom-free and approximately 3 weeks out from symptom onset at the time of the initial blood draw. Following enrollment, whole blood from convalescent donors was run on a colloidal-gold immunochromatographic 'lateral flow' assay to evaluate for prior exposure to SARS-CoV-2. This assay detects IgM or IgG antibodies directed against recombinant SARS-CoV-2 antigen labeled with a colloidal gold tracer (20/20 BioResponse CoronaCheck). Ninety percent of convalescent donors tested positive for IgM or IgG to SARS-CoV-2 by this assay (Table 1) .

Convalescent donors were California residents, who were either referred to the study by a health care provider or self-referred. The majority (75%) of donors had a known sick contact with COVID-19 or suspected exposure to SARS-CoV-2 ( Table 1) . The most common symptoms reported were cough, fatigue, fever, anosmia, and dyspnea. Seventy percent of donors experienced mild illness. Donors were asked to self-report any known medical illnesses. Of note, 65% of these individuals had no known underlying medical illnesses.

Epitope MegaPool (MP) design and preparation. SARS-CoV-2 virus-specific CD4 and CD8 peptides were synthesized as crude material (A&A, San Diego, CA), resuspended in DMSO, pooled and sequentially lyophilized as previously reported . SARS-CoV-2 epitopes were predicted using the protein sequences derived from the SARS-CoV-2 reference (GenBank: MN908947) and IEDB analysis-resource as previously described Grifoni et al., 2020) . Specifically, CD4 SARS-CoV-2 epitope prediction was carried out using a previously described approach in Tepitool resource in IEDB Paul et al., 2016) , to select peptides with median consensus percentile ≤ 20, similar to what was previously described, but removing the resulting spike glycoprotein epitopes from this prediction (CD4-R (remainder) "Non-spike" MP, n=221). This approach takes advantage of the extensive cross-reactivity and repertoire overlap between different HLA class II loci and allelic variants to predict promiscuous epitopes, capable of binding across the most common HLA class II prototypic specificities (Greenbaum et al., 2011; O'Sullivan et al., 1991; Sidney et al., 2010a, b; Southwood et al., 1998) . The algorithm utilizes predictions for seven common HLA-DR alleles (DRB1*03:01, DRB1*07:01, DRB1*15:01, DRB3*01:01, DRB3*02:02, DRB4*01:01 and DRB5*01:01) empirically determined to allow coverage of diverse populations and for different pathogens and antigen systems (Dhanda et al., 2018; Paul et al., 2015) .

To investigate in-depth spike-specific CD4 T cells, 15-mer peptides (overlapping by 10 amino acids) spanning the entire antigen have been synthesized and pooled separately (CD-4 S (spike) MP, n=253).

In the case of CD8 epitopes, since the overlap between different HLA class I allelic variants and loci is more limited to specific groups of alleles, or supertypes (Sidney et al., 2008) , we targeted a set of the 12 most prominent HLA class I A and B alleles (A*01:01, A*02:01, A*03:01, A*11:01, A*23:01, A*24:02, B*07:02, B*08:01, B*35:01, B*40:01, B*44:02, B*44:03), which have been shown to allow broad coverage of the general population. CD8 SARS-CoV-2 epitope prediction was performed as previously reported, using NetMHC pan EL 4.0 algorithm (Jurtz et al., 2017) for the top 12 more frequent HLA alleles and selecting the top 1 percentile predicted epitope per HLA allele clustered with nested/overlap reduction (Grifoni et al., 2020) . The 628 predicted CD8 epitopes were split in two CD8 MPs containing 314 peptides each (CD8-A and CD8-B). The CMV MP is a pool of previously reported class I and class II epitopes .

Protein peptide pools. In the case of the protein pools, peptides of 15 amino acid length overlapping by 10 spanning each entire protein sequence were tested in a single MP (6-253 peptides per pool). Table  S1 lists the number of peptides pooled for each of the viral proteins. Upon request we are prepared to make these MP available to the scientific community for use in a diverse set of investigations.

For all samples whole blood was collected in ACD tubes (COVID-19 donors) or heparin coated blood bag (healthy unexposed donors). Whole blood was then centrifuged for 15 min at 1850 rpm to separate the cellular fraction and plasma. The plasma was then carefully removed from the cell pellet and stored at -20C.

Peripheral blood mononuclear cells (PBMC) were isolated by density-gradient sedimentation using Ficoll-Paque (Lymphoprep, Nycomed Pharma, Oslo, Norway) as previously described (Weiskopf et al., 2013) . Isolated PBMC were cryopreserved in cell recovery media containing 10% DMSO (Gibco), supplemented with 10% heat inactivated fetal bovine serum, depending on the processing laboratory, (FBS; Hyclone Laboratories, Logan UT) and stored in liquid nitrogen until used in the assays.

SARS-CoV-2 Receptor Binding Domain (RBD) protein was obtained courtesy of Florian Krammer and Peter Kim (Stadlbauer et al., 2020) . Corning 96-well half-area plates (ThermoFisher 3690) were coated with 1µg/mL SARS-CoV-2 RBD overnight at 4°C. ELISA protocol generally followed that of the Krammer lab, which previously demonstrated specificity (Stadlbauer et al., 2020) . Plates were blocked the next day with 3% milk (Skim Milk Powder ThermoFisher LP0031 by weight/volume) in Phosphate Buffered Saline (PBS) containing 0.05% Tween-20 (ThermoScientific J260605-AP) for 2 hours at room temperature. Plasma was then added to the plates and incubated for 1.5 hours at room temperature. Prior to plasma addition to the plates, plasma was heat inactivated at 56°C for 30-60 minutes. Plasma was diluted in 1% milk in 0.05% PBS-Tween 20 starting at a 1:3 dilution and diluting each sample at by 1:3. Plates were then washed 5 times with 0.05% PBS-Tween 20. Secondary antibodies were diluted in 1% milk in 0.05% Tween-20 and incubated for 1 hour. For IgG, anti-human IgG peroxidase antibody produced in goat (Sigma A6029) was used at a 1:5000 dilution. For IgM, anti-human IgM peroxidase antibody produced in goat (Sigma A6907) was used at a 1:10,000 dilution. For IgA, anti-human IgA horseradish peroxidase antibody (Hybridoma Reagent Laboratory HP6123-HRP) was used at a 1:1,000 dilution. Plates were washed 5 times with 0.05% PBS-Tween 20. Plates were developed with TMB Substrate Kit (ThermoScientific 34021) for 15 minutes at room temperature. The reaction was stopped with 2M sulfuric acid. Plates were read on a Spectramax Plate Reader at 450 nm using SoftMax Pro, and ODs were background subtracted. A positive control standard was created by pooling plasma from six convalescing COVID-19 patients. Positive control standard was run on each plate and was used to calculate titers (relative units) for all samples using non-linear regression interpolations, done to quantify the amount of anti-RBD IgG, anti-RBD IgM, and anti-RBD IgA present in each specimen. Titers were plotted for each specimen and compared to COVID-19 negative specimens. As a second analytical approach, Area under the curve was also calculated for each specimen to compare COVID-19 to negative specimens, using a baseline of 0.05 for peak calculations.

An in-house ELISA at UNC was performed by coating with recombinant S RBD antigens (SARS-CoV-2, SARS-CoV, OC43-CoV and NL63-CoV) in TBS for 1 h at 37 C. After blocking, we added 1:20 diluted serum and incubated at 37°C for 1 h. Antigen-specific antibodies (Ig) were measured at 405 nm by using alkaline phosphatase conjugated goat anti-human IgG, IgA and IgM Abs and 4-Nitrophenyl phosphate.

Direct ex vivo PBMC immune cell phenotyping. For the surface stain, 1x10 6 PBMCs were resuspended in 100 µl PBS with 2% FBS (FACS buffer) and stained with antibody cocktail for 1 hour at 4°C in the dark. Following surface staining, cells were washed twice with FACS buffer. Cells were then fixed/permeabilized for 40min at 4C in the dark using the eBioscience FoxP3 transcription factor buffer kit (ThermoFisher Scientific, Waltham, MA). Following fixation/permeabilization, cells were washed twice with 1x permeabilization buffer, resuspended in 100 µl permeabilization buffer and stained with intracellular/intranuclear antibodies for 1 hour at 4ºC in the dark. Samples were washed twice with 1x permeabilization buffer following staining. After the final wash, cells were resuspended in 200µl FACS buffer. All samples were acquired on a BD FACSymphony cell sorter (BD Biosciences, San Diego, CA). A list of antibodies used in this panel can be found in Table S2 .

T cell stimulations. For all flow cytometry assays of stimulated T cells, cryopreserved cells were thawed by diluting them in 10 ml complete RPMI 1640 with 5% human AB serum (Gemini Bioproducts) in the presence of benzonase [20ul/10mL]. All samples were acquired on a ZE5 Cell analyzer (Bio-rad laboratories), and analyzed with FlowJo software (Tree Star, San Carlos, CA).

Activation induced cell marker assay. Cells were cultured for 24 hours in the presence of SARS-CoV-2 specific MPs [1 µg/ml] or 10 µg/mL PHA in 96-wells U bottom plates at 1x10 6 PBMC per well. A stimulation with an equimolar amount of DMSO was performed as negative control, phytohemagglutinin (PHA, Roche, 1µg/ml) and stimulation with a combined CD4 and CD8 cytomegalovirus MP (CMV, 1µg/ml) were included as positive controls. Supernatants were harvested at 24 hours post-stimulation for multiplex detection of cytokines. Antibodies used in the AIM assay are listed in Table S4 . AIM assays shown in Figure 2 and 3 and AIM assays shown in Figure 6 had five COVID-19 donors in common and nine Unexposed donors. Full raw data is listed in Table S6 .

Intracellular cytokine staining assay. For the intracellular cytokine staining, PBMC were cultured in the presence of SARS-CoV-2 specific MPs [1 µg/ml] for 9 hours. Golgi-Plug containing brefeldin A (BD Biosciences, San Diego, CA) and monensin (Biolegend, San Diego, CA) were added 3 hours into the culture. Cells were then washed and surface stained for 30 minutes on ice, fixed with 1% of paraformaldehyde (Sigma-Aldrich, St. Louis, MO) and kept at 4 o C overnight. Antibodies used in the ICS assay are listed in Table S5 . The gates applied for the identification of IFNγ, GzB, TNFα, or IL-10 production on the total population of CD8 + T cells were defined according to the cells cultured with DMSO for each individual.

Supernatants were collected from 24-hour stimulation cultures of the AIM assays and stored in 96 well plates at -20 °C. Cytokines in cell culture supernatants of the same samples used for AIM were quantified using a human Th cytokine panel (13-plex) kit (LEGENDplex, Biolegend) according to the manufacturer's instruction. Supernatants were mixed with beads coated with capture antibodies specific for IL-5, IL-13, IL-2, IL-6, IL-9, IL-10, IFNγ, TNFα, IL-17a, IL-17F, IL-4, IL-21 and IL-22 and incubated on a 96 well filter plate for 2 hours. Beads were washed and incubated with biotin-labelled detection antibodies for 1 hour, followed by a final incubation with streptavidin-PE. Beads were analyzed by flow cytometry using a FACS Canto cytometer. Analysis was performed using the LEGENDplex analysis software v8.0, which distinguishes between the 13 different analytes on basis of bead size and internal dye.

To identify coronavirus epitopes and associated references, the IEDB was searched (on April 16, 2020) utilizing the following queries. A first query was run to identify references associated with class I restricted CD8 epitopes, which utilized the criteria settings "Antigen": Organism = Coronavirus (taxonomy ID 11118); "Assay": Positive assays only; "Assay": T cell assay; "MHC restriction" = MHC Class II; no parameters were defined for "Host" or "Disease." This query identified 57 references, which are listed and displayed under the "References" tab on the results page.

A second query was run to identify references associated with class II restricted CD4 epitopes which utilized the criteria settings "Antigen": Organism = Coronavirus (taxonomy ID 11118); "Assay": Positive assays only; "Assay": T cell assay; "MHC restriction" = MHC Class II; no parameters were defined for "Host" or "Disease." This query identified 27 references, which are listed and displayed under the "References" tab on the results page.

A third query was run to specifically capture epitopes and map them back to the antigen of origin using the setting; "Antigen": Organism = Coronavirus (taxonomy ID 11118); "Assay": Positive assays only; "Assay": T cell assay; no parameters were defined for "MHC restriction," "Host" or "Disease." Results were exported as csv files, and then examined in Excel to tabulate the number of CD4 and CD8 epitopes recognized in humans, mice, transgenic mice and other hosts associated with each respective antigen.

Data and statistical analyses were done in FlowJo 10 and GraphPad Prism 8.4, unless otherwise stated. The statistical details of the experiments are provided in the respective figure legends. Data plotted in linear scale were expressed as Mean + Standard Deviation (SD). Data plotted in logarithmic scales were expressed as Geometric Mean + Geometric Standard Deviation (SD). Correlation analyses were performed using Spearman, while Mann-Whitney or Wilcoxon tests were applied for unpaired or paired comparisons, respectively. Details pertaining to significance are also noted in the respective legends. T cell data have been calculated as background subtracted data or stimulation index. Background subtracted data were derived by subtracting the percentage of AIM + cells after SARS-CoV-2 stimulation from the DMSO stimulation. Stimulation Index was calculated instead by dividing the percentage of AIM+ cells after SARS-CoV-2 stimulation with the percentage of AIM + cells derived from DMSO stimulation. If the AIM + cells percentage after DMSO stimulation was equal to 0, the minimum value across each cohort was used. When two stimuli were combined together, the percentage of AIM + cells after SARS-CoV-2 stimulation was combined and either subtracted twice or divided by twice the value of the percentage of AIM+ cells derived from DMSO stimulation. Additional data analysis techniques are described in the Methods sections above. 

High Prevalence of MERS-CoV Infection in Camel Workers in Saudi Arabia

SARS-CoV-2 Vaccines: Status Report

Immune history profoundly affects broadly protective B cell responses to influenza

Th1 versus Th2 T cell polarization by whole-cell and acellular childhood pertussis vaccines persists upon re-immunization in adolescence and adulthood

SARS-CoV-2 launches a unique transcriptional signature from in vitro, ex vivo, and in vivo systems

The time course of the immune response to experimental coronavirus infection of man

COVID-19: immunopathology and its implications for therapy

Automatic Generation of Validated Specific Epitope Sets

MERS-CoV Antibody Responses 1 Year after Symptom Onset

T Follicular Helper Cell Biology: A Decade of Discovery and Diseases

Definition of Human Epitopes Recognized in Tetanus Toxoid and Development of an Assay Strategy to Detect Ex Vivo Tetanus CD4+ T Cell Responses

Recurrent group A Streptococcus tonsillitis is an immunosusceptibility disease involving antibody deficiency and aberrant TFH cells

A Cytokine-Independent Approach To Identify Antigen-Specific Human Germinal Center T Follicular Helper Cells and Rare Antigen-Specific CD4+ T Cells in Blood

Predicting HLA CD4 Immunogenicity in Human Populations. Front Immunol 9

IEDB-AR: immune epitope database-analysis resource in 2019

An interactive web-based dashboard to track COVID-19 in real time

Complex Immune Dysregulation in COVID-19 Patients with Severe Respiratory Failure

Prevalence of Antibodies to Four Human Coronaviruses Is Lower in Nasal Secretions than in Serum

Potent protection against H5N1 and H7N9 influenza via childhood hemagglutinin imprinting

Functional classification of class II human leukocyte antigen (HLA) molecules reveals seven different supertypes and a surprising degree of repertoire sharing across supertypes

Pre-existing immunity against swine-origin H1N1 influenza viruses in the general human population

Global Assessment of Dengue Virus-Specific CD4+ T Cell Responses in Dengue-Endemic Areas

A Sequence Homology and Bioinformatic Approach Can Predict Candidate Targets for Immune Responses to SARS-CoV-2

Cardiovascular Implications of Fatal Outcomes of Patients With Coronavirus Disease 2019 (COVID-19)

Long-Term Persistence of IgG Antibodies in SARS-CoV Infected Healthcare Workers. medRxiv

Cross-reactive antibody responses to the 2009 pandemic H1N1 influenza virus

Cytokine-Independent Detection of Antigen-Specific Germinal Center T Follicular Helper Cells in Immunized Nonhuman Primates Using a Live Cell Activation-Induced Marker Technique

Successive annual influenza vaccination induces a recurrent oligoclonotypic memory response in circulating T follicular helper cells

Lack of allergy to timothy grass pollen is not a passive phenomenon but associated with the allergen-specific modulation of immune reactivity

A systematic review of antibody mediated immunity to coronaviruses: antibody kinetics, correlates of protection, and association of antibody responses with severity of disease. medRxiv

High-Resolution Analysis of Coronavirus Gene Expression by RNA Sequencing and Ribosome Profiling

NetMHCpan-4.0: Improved Peptide-MHC Class I Interaction Predictions Integrating Eluted Ligand and Peptide Binding Affinity Data

Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period

T cell responses to whole SARS coronavirus in humans

A Quantitative Analysis of Complexity of Human Pathogen-Specific CD4 T Cell Responses in Healthy M. tuberculosis Infected South Africans

Antispike IgG causes severe acute lung injury by skewing macrophage responses during acute SARS-CoV infection

Altered differentiation is central to HIV-specific CD4(+) T cell dysfunction in progressive disease

Uncovering the interplay between CD8, CD4 and antibody responses to complex pathogens

On the interaction of promiscuous antigenic peptides with different DR alleles. Identification of common structural motifs

Sensitive and Specific Detection of Low-Level Antibody Responses in Mild Middle East Respiratory Syndrome Coronavirus Infections

Development and validation of a broad scheme for prediction of HLA class II restricted T cell epitopes

TepiTool: A Pipeline for Computational Prediction of T Cell Epitope Candidates

News Feature: Avoiding pitfalls in the pursuit of a COVID-19 vaccine

Comparative analysis of activation induced marker (AIM) assays for sensitive identification of antigen-specific CD4 T cells

From vaccines to memory and back

Selective CD4+ T cell help for antibody responses to a large viral pathogen: deterministic linkage of specificities

Development of a nucleocapsid-based human coronavirus immunoassay and estimates of individuals exposed to coronavirus in a U.S. metropolitan population

HLA class I supertypes: a revised and updated classification

Divergent motifs but overlapping binding repertoires of six HLA-DQ molecules frequently expressed in the worldwide human population

Five HLA-DP molecules frequently expressed in the worldwide human population share a common HLA supertypic binding specificity

Unique and conserved features of genome and proteome of SARScoronavirus, an early split-off from the coronavirus group 2 lineage

Several common HLA-DR types share largely overlapping peptide binding repertoires

Cellular immune correlates of protection against symptomatic pandemic influenza

SARS-CoV-2 Seroconversion in Humans: A Detailed Protocol for a Serological Assay, Antigen Production, and Test Setup

Antibody-dependent enhancement occurs upon re-infection with the identical serotype virus in feline infectious peritonitis virus infection

The COVID-19 vaccine development landscape

Human T Cell Response to Dengue Virus Infection

Early death after feline infectious peritonitis virus challenge due to recombinant vaccinia virus immunization

The Immune Epitope Database (IEDB): 2018 update

Comprehensive analysis of dengue virus-specific responses supports an HLA-linked protective role for CD8+ T cells

Human CD8+ T-Cell Responses Against the 4 Dengue Virus Serotypes Are Associated With Distinct Patterns of Protein Targets

Preexisting influenza-specific CD4+ T cells correlate with disease protection against influenza challenge in humans

Plasma inflammatory cytokines and chemokines in severe acute respiratory syndrome

Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72 314 Cases From the Chinese Center for Disease Control and Prevention

An Infectious cDNA Clone of SARS-CoV-2

Recovery from the Middle East respiratory syndrome is associated with antibody and T-cell responses

Age-related increases in PGD(2) expression impair respiratory DC migration, resulting in diminished T cell responses upon respiratory virus infection in mice

Airway Memory CD4(+) T Cells Mediate Protective Immunity against Emerging Respiratory Coronaviruses

T cell responses are required for protection from clinical disease and for virus clearance in severe acute respiratory syndrome coronavirus-infected mice

Evasion by stealth: inefficient immune activation underlies poor T cell response and severe disease in SARS-CoV-infected mice

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

Epitope pools detect CD4 + and CD8 + T cells in 100 and 70% of convalescent COVID patients • T cell responses are focused not only on spike but also on M, N and other ORFs • T cell reactivity to SARS-CoV-2 epitopes is also detected in non-exposed individuals An analysis of immune cell responses to SARS-CoV-2

We would like to thank Cheryl Kim, Director of the LJI flow cytometry core facility for outstanding expertise. We thank Prof. Peter Kim, Abigail Powell, PhD, and colleagues (Stanford) for RBD protein synthesized from Prof. Florian Krammer (Mt. Sinai) constructs. J.M. was supported by Ph.D. student fellowships from the Departamento Administrativo de Ciencia, Tecnología e Innovación (COLCIENCIAS) and Pontificia Universidad Javeriana (Convocatoria 727 Doctorados Nacionales). This work was funded by the NIH NIAID under awards AI42742 (Cooperative Centers for Human Immunology) (S.C., 

The authors declare no competing interests.