Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Minireviews
    • JVI Classic Spotlights
    • Archive
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About JVI
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Journal of Virology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Minireviews
    • JVI Classic Spotlights
    • Archive
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About JVI
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Genetic Diversity and Evolution

Human Leukocyte Antigen Susceptibility Map for Severe Acute Respiratory Syndrome Coronavirus 2

Austin Nguyen, Julianne K. David, Sean K. Maden, Mary A. Wood, Benjamin R. Weeder, Abhinav Nellore, Reid F. Thompson
Tom Gallagher, Editor
Austin Nguyen
aComputational Biology Program, Oregon Health & Science University, Portland, Oregon, USA
bDepartment of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Austin Nguyen
Julianne K. David
aComputational Biology Program, Oregon Health & Science University, Portland, Oregon, USA
bDepartment of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Julianne K. David
Sean K. Maden
aComputational Biology Program, Oregon Health & Science University, Portland, Oregon, USA
bDepartment of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sean K. Maden
Mary A. Wood
aComputational Biology Program, Oregon Health & Science University, Portland, Oregon, USA
cPortland VA Research Foundation, Portland, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Mary A. Wood
Benjamin R. Weeder
aComputational Biology Program, Oregon Health & Science University, Portland, Oregon, USA
bDepartment of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Benjamin R. Weeder
Abhinav Nellore
aComputational Biology Program, Oregon Health & Science University, Portland, Oregon, USA
bDepartment of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
dDepartment of Surgery, Oregon Health & Science University, Portland, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Abhinav Nellore
Reid F. Thompson
aComputational Biology Program, Oregon Health & Science University, Portland, Oregon, USA
bDepartment of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
eDepartment of Radiation Medicine, Oregon Health & Science University, Portland, Oregon, USA
fDepartment of Medical Informatics and Clinical Epidemiology, Oregon Health & Science University, Portland, Oregon, USA
gDivision of Hospital and Specialty Medicine, VA Portland Healthcare System, Portland, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Reid F. Thompson
Tom Gallagher
Loyola University Chicago
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/JVI.00510-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Genetic variability across the three major histocompatibility complex (MHC) class I genes (human leukocyte antigen A [HLA-A], -B, and -C genes) may affect susceptibility to and severity of the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for coronavirus disease 2019 (COVID-19). We performed a comprehensive in silico analysis of viral peptide-MHC class I binding affinity across 145 HLA-A, -B, and -C genotypes for all SARS-CoV-2 peptides. We further explored the potential for cross-protective immunity conferred by prior exposure to four common human coronaviruses. The SARS-CoV-2 proteome was successfully sampled and was represented by a diversity of HLA alleles. However, we found that HLA-B*46:01 had the fewest predicted binding peptides for SARS-CoV-2, suggesting that individuals with this allele may be particularly vulnerable to COVID-19, as they were previously shown to be for SARS (M. Lin, H.-T. Tseng, J. A. Trejaut, H.-L. Lee, et al., BMC Med Genet 4:9, 2003, https://bmcmedgenet.biomedcentral.com/articles/10.1186/1471-2350-4-9). Conversely, we found that HLA-B*15:03 showed the greatest capacity to present highly conserved SARS-CoV-2 peptides that are shared among common human coronaviruses, suggesting that it could enable cross-protective T-cell-based immunity. Finally, we reported global distributions of HLA types with potential epidemiological ramifications in the setting of the current pandemic.

IMPORTANCE Individual genetic variation may help to explain different immune responses to a virus across a population. In particular, understanding how variation in HLA may affect the course of COVID-19 could help identify individuals at higher risk from the disease. HLA typing can be fast and inexpensive. Pairing HLA typing with COVID-19 testing where feasible could improve assessment of severity of viral disease in the population. Following the development of a vaccine against SARS-CoV-2, the virus that causes COVID-19, individuals with high-risk HLA types could be prioritized for vaccination.

INTRODUCTION

Recently, a new strain of betacoronavirus (severe acute respiratory syndrome coronavirus 2, or SARS-CoV-2) emerged as a global pathogen, prompting the World Health Organization in January 2020 to declare an international public health emergency (1). In the large coronavirus family, comprising enveloped positive-strand RNA viruses, SARS-CoV-2 is the seventh encountered strain that causes respiratory disease in humans (2) ranging from mild—the common cold—to severe—disease caused by the zoonotic Middle East respiratory syndrome coronavirus (MERS-CoV) and severe acute respiratory syndrome coronavirus (SARS-CoV). As of April 2020, there are over one million presumed or confirmed cases of coronavirus disease 19 (COVID-19) worldwide, with total deaths exceeding 50,000 (3). While age and many comorbidities, including cardiovascular and pulmonary disease, appear to increase the severity and mortality of COVID-19 (4–9), approximately 80% of infected individuals have mild symptoms (10). As with SARS-CoV (11, 12) and MERS-CoV (13, 14), children seem to have low susceptibility to the disease (15–17); despite infection rates similar to those seen with adults (18), only 5.9% of pediatric cases are severe or critical, possibly due to lower binding ability of the ACE2 receptor in children or generally higher levels of antiviral antibodies (19). Other similarities (20–22), including genomic (23, 24) and immune system response (25–33) similarities, between SARS-CoV-2 and other coronaviruses (34), especially SARS-CoV and MERS-CoV, are topics of ongoing active research, results of which may inform an understanding of the severity of infection (35) and improve the ongoing work of immune landscape profiling (36–40) and vaccine discovery (28, 37, 41–48).

Human leukocyte antigen (HLA) alleles, which are critical components of the viral antigen presentation pathway, have been shown in previous studies to confer differential viral susceptibility and severity of disease. For instance, disease caused by the closely related SARS-CoV (23, 24) shows increased severity among individuals with the HLA-B*46:01 genotype (49). Associations between HLA genotype and disease severity extend broadly to several other unrelated viruses. For example, in human immunodeficiency virus 1 (HIV-1), certain HLA types (e.g., HLA-A*02:05) may reduce risk of seroconversion (50), and in dengue virus, certain HLA alleles (e.g., HLA-A*02:07 and HLA-B*51) are associated with increased secondary disease severity among ethnic Thais (51).

While the details of the clinical picture of the COVID-19 pandemic continue to emerge, there remain substantial unanswered questions regarding the role of individual genetic variability in the immune response against SARS-CoV-2 (51). We hypothesize that individual HLA genotypes may differentially induce the T-cell mediated antiviral response and could potentially alter the course of disease and its transmission. In this study, we performed a comprehensive in silico analysis of viral peptide-major histocompatibility complex (MHC) class I binding affinity across 145 different HLA types for the entire SARS-CoV-2 proteome.

RESULTS

To explore the potential for a given HLA allele to produce an antiviral response, we assessed the HLA binding affinity of all possible 8-mers to 12-mers from the SARS-CoV-2 proteome (n = 48,395 unique peptides). We then removed from further consideration 16,138 peptides that were not predicted to enter the MHC class I antigen processing pathway via proteasomal cleavage. For the remaining 32,257 peptides, we repeated binding affinity predictions for a total of 145 different HLA types, and we show here the SARS-CoV-2-specific distribution of per-allele proteome presentation (predicted binding affinity threshold of <500 nM) (Fig. 1; see also Table S1 in the supplemental material). Importantly, we note that the putative capacity for SARS-CoV-2 antigen presentation is unrelated to the HLA allelic frequency in the population (Fig. 1). We identify HLA-B*46:01 as the HLA allele with the fewest predicted binding peptides for SARS-CoV-2. We performed the same analyses for the closely related SARS-CoV proteome (see Fig. S1 in the supplemental material) and similarly note that HLA-B*46:01 was predicted to present the fewest SARS-CoV peptides, in keeping with previous clinical data associating this allele with severe disease (49).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Distribution of HLA allelic presentation of 8- to 12-mers from the SARS-CoV-2 proteome. At right, the number of peptides (see Table S1) that putatively bind to each of 145 HLA alleles is shown as a series of horizontal bars, with dark and light shading indicating the number of tightly (<50 nM) and loosely (<500 nM) binding peptides, respectively, and with green, orange, and purple representing HLA-A, -B, and -C alleles, respectively. Alleles are sorted in descending order based on the number of peptides that they bind (<500 nM). The corresponding estimated allelic frequency in the global population is also shown (left), with the length of each horizontal bar indicating absolute frequency in the population.

To assess the potential for cross-protective immunity conferred by prior exposure to common human coronaviruses (i.e., HKU1, OC43, NL63, and 229E), we next sought to characterize the conservation of the SARS-CoV-2 proteome across diverse coronavirus subgenera to identify highly conserved linear epitopes. After aligning reference proteome sequence data for 5 essential viral components (ORF1ab, S, E, M, and N proteins) across 34 distinct alpha- and betacoronaviruses, including all known human coronaviruses, we identified 48 highly conserved amino acid sequence spans (see Data File S1 in the supplemental material). Acknowledging the challenges to inferring cross-protective immunity among closely related peptides, we confined our attention exclusively to identical peptide matches. Among the conserved sequences, 44 SARS-CoV-2 sequences would each be anticipated to generate at least one 8- to 12-mer linear peptide epitope also present within at least one other common human coronavirus (Fig. 2; see also Table S2). In total, 564 such 8- to 12-mer peptides were found to share 100% identity with corresponding OC43, HKU1, NL63, and 229E sequences (467, 460, 179, and 157 peptides, respectively) (Table S3).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Amino acid sequence conservation of four linear peptide example sequences from three human coronavirus proteins. Protein sequence alignments are shown for nucleocapsid (N), membrane (M), and ORF1ab polyprotein (helicase) across all five known human betacoronaviruses (SARS-CoV-2, SARS-CoV, HKU1, OC43, and MERS-CoV) and two known human alphacoronaviruses (229E and NL63). Each row in the three depicted sequence alignments corresponds to the protein sequence from the indicated coronavirus, with the starting coordinate of the viral protein sequence shown at left and position coordinates of the overall alignment displayed above. Blue shading indicates the extent of sequence identity, with the darkest blue shading indicating a 100% match for that amino acid across all sequences. The four red-highlighted sequences correspond to highly conserved peptides ≥8 amino acids in length (PRWYFYYLGTGP, WSFNPETN, QPPGTGKSH, and VYTACSHAAVDALCEKA, see Table S2).

For the subset of these potentially cross-protective peptides that are anticipated to be generated via the MHC class I antigen processing pathway, we performed binding affinity predictions across 145 different HLA-A, -B, and -C alleles (see Data File S3). As described above, we demonstrated the SARS-CoV-2-specific distribution of per-allele presentation for these conserved peptides. We found that alleles HLA-A*02:02, HLA-B*15:03, and HLA-C*12:03 were the top presenters of conserved peptides. Conversely, we note that 56 different HLA alleles demonstrated no appreciable binding affinity (<500 nM) to any of the conserved SARS-CoV-2 peptides, suggesting a concomitant lack of potential for cross-protective immunity from other human coronaviruses. We note, in particular, that HLA-B*46:01 was among these alleles. We note also that the putative capacity for conserved peptide presentation is unrelated to the HLA allelic frequency in the population (Fig. 3). Moreover, we see no appreciable global correlation between conservation of the SARS-CoV-2 proteome and its predicted MHC binding affinity, suggesting a lack of selective pressure for or against the capacity to present coronavirus epitopes (P = 0.27 [Fisher’s exact test]; see Fig. S2).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Distribution of HLA allelic presentations of highly conserved human coronavirus peptides with potential to elicit cross-protective immunity to COVID-19. At right, the number of conserved peptides (see Table S3) that putatively bind to a subset of 89 HLA alleles is shown as a series of horizontal bars, with dark and light shading indicating the number of tightly (<50 nM) and loosely (<500 nM) binding peptides, respectively, and with green, orange, and purple representing HLA-A, -B, and -C alleles, respectively. Alleles are sorted in descending order based on the number of peptides they are anticipated to present (binding affinity, <500 nM). The corresponding allelic frequency in the global population is also shown (left), with the length of each horizontal bar indicating absolute frequency in the population.

We were further interested in whether certain regions of the SARS-CoV-2 proteome showed differential presentation by the MHC class I pathway. Accordingly, we surveyed the distribution of antigen presentation capacity across the entire proteome, highlighting its most conserved peptide sequences (Fig. 4). Throughout the entire proteome, HLA-A and HLA-C alleles exhibited the relatively largest and smallest capacities to present SARS-CoV-2 antigens, respectively. However, each of the three major class I genes exhibited very similar patterns of peptide presentation across the proteome (Fig. S3). We additionally note that peptide presentation appears to be independent of estimated time of peptide production during viral life cycle, with indistinguishable levels of peptide presentation of both early and late SARS-CoV-2 peptides (Fig. S4).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Distribution of allelic presentation of conserved 8- to 12-mers across the entire SARS-CoV-2 proteome for all HLA alleles and individually for HLA-A, HLA-B, and HLA-C (first, second, third, and fourth plots from top, respectively) with dark and light shading indicating the number of tightly (<50 nM) and loosely (<500 nM) binding peptides, respectively. Positions are derived from a concatenation of coding sequences (CDSs) as indicated in the bottom panel. Tightly binding peptides are confined to ORF1ab. The sequence begins with only the last 12 amino acids of ORF1a because all but the last four amino acids of ORF1a are contained in ORF1ab, and we considered binding peptides up to 12 amino acids (AA) in length.

Given the global nature of the current COVID-19 pandemic, we sought to describe population-level distributions of the HLA alleles most (and least) capable of generating a repertoire of SARS-CoV-2 epitopes in support of a T-cell-based immune response. While we present global maps of individual HLA allele frequencies for the full set of 145 different alleles studied here (Data File S2), we specifically highlight the global distributions of the three best-presenting (A*02:02, B*15:03, and C*12:03) and three of the worst-presenting (A*25:01, B*46:01, and C*01:02) HLA-A, -B, and -C alleles (Fig. 5). Note that all allelic frequencies are aggregated by country but that they implicitly reflect the distribution of HLA data available on the Allele Frequency Net Database (52).

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Global HLA allele frequency distribution heat maps for six HLA-A, -B, and -C alleles. The leftmost panels show the global allele frequency distributions by country for three representative alleles (HLA-A*02:02, HLA-B*15:03, and HLA-C*12:03) with the predicted capacities to present the greatest repertoire of epitopes from the SARS-CoV-2 proteome (21.1%, 19.1%, and 7.9% of presentable epitopes, respectively). The rightmost panels show the global allele frequency distributions by country for three representative alleles (HLA-A*25:01, HLA-B*46:01, and HLA-C*01:02) with the lowest predicted levels epitope presentation from the SARS-CoV-2 proteome (0.2%, 0%, and 0% of presentable epitopes, respectively). Heat map coloring corresponds to the individual HLA allele frequency within each country, ranging from lowest (white/yellow) to highest (red) frequency as indicated in the legend below each map.

Finally, we acknowledge that nearly all individuals have two HLA-A/B/C haplotypes constituting as few as three but as many as six distinct alleles, potentially buffering against the lack of presentation from a single poorly presenting allele. We sought to determine whether allele-specific variability in SARS-CoV-2 presentation extends to full HLA haplotypes and to whole individual HLA genotypes. For six representative alleles with the highest (HLA-A*02:02, HLA-B*15:03, and HLA-C*12:03) and lowest (HLA-A*25:01, HLA-B*46:01, and HLA-C*01:02) predicted capacity for SARS-CoV-2 epitope presentation, these differences remain significant at the haplotype level, albeit with wide variability in presentation among different haplotypes (Fig. 6). Haplotype-level data for all 145 alleles are included in Fig. S5 and Data File S2. We then identified 3,382 individuals with full HLA genotype data and noted wide variability in their capacity to present peptides from the SARS-CoV-2 proteome, albeit with a small minority of individuals at either extreme (Fig. S6).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Distributions of SARS-CoV-2 peptide presentation across HLA haplotypes. The leftmost panels show the distributions of SARS-CoV-2 peptide presentation capacity for haplotypes containing one of three representative HLA alleles (HLA-A*02:02, HLA-B*15:03, and HLA-C*12:03) with the greatest predicted repertoire of epitopes from the SARS-CoV-2 proteome. The rightmost panels show the distributions of SARS-CoV-2 peptide presentation capacity for haplotypes containing one of three representative alleles (HLA-A*25:01, HLA-B*46:01, and HLA-C*01:02) with the lowest predicted levels of epitope presentation from the SARS-CoV-2 proteome. Black and gray bars represent full and partial haplotypes, respectively. Blue and red dashed lines represent the percentages of presented SARS-CoV-2 peptides for the indicated allele itself (blue) and its global population frequency weighted average presentation across its observed haplotypes (red).

DISCUSSION

To the best of our knowledge, this was the first study to evaluate per-allele viral proteome presentation across a wide range of HLA alleles using MHC-peptide binding affinity predictors. This report also introduces the relationship between coronavirus sequence conservation and MHC class I antigen presentation. We show that individual HLA, haplotype, and full-genotype variability likely influence the capacity to respond to SARS-CoV-2 infection, and we note certain alleles in particular (e.g., HLA-B*46:01) that could be associated with more-severe infection, as previously shown with SARS-CoV (49). Indeed, we further compare SARS-CoV and SARS-CoV-2 peptide presentation data and note a high degree of similarity between the two across HLA types. Finally, this is the first report to present global distributions of HLA types and haplotypes with potential epidemiological ramifications in the setting of the current pandemic. We found that, in general, there is no correlation between the HLA allelic frequency in the population and allelic capacity to bind SARS-CoV or SARS-CoV-2 peptides, irrespective of the estimated timing of peptide production during the viral replication cycle. While we are not aware of any studies explicitly reporting the relationship between human coronavirus epitope abundance and immune response, there are vaccinia virus data that suggest that early peptide antigens are more likely to generate CD8+ T-cell responses whereas antibody and CD4+ T-cell responses are more likely to target later mRNA expression with higher peptide abundance in the virion (53).

We note, however, several limitations to our work. First and foremost, while we note that a few of our binding affinity predictions were borne out in experimentally validated SARS-CoV peptides (see Table S4 in the supplemental material), we acknowledge that ours was a study performed entirely in silico. As we are unable to obtain individual-level HLA typing and clinical outcome data for any real-world COVID-19 populations at this time, the data presented are theoretical in nature and are subject to many of the same limitations implicit in the MHC binding affinity prediction tool(s) upon which it is based. As such, we are unable to assess the relative importance of HLA type compared to known disease-modifying risk factors such as age and clinical comorbidities (4–9). We further note that peptide-MHC binding affinity is limited in its utility as a predictor of subsequent T-cell responses (54–56), and we did not study T-cell responses here. As such, we are ill-equipped to explore phenomena such as original antigenic sin (57–59), where prior exposure to a closely related infection(s) might trigger T-cell anergy (60–62) or immunopathogenesis (63) in the setting of a novel infection. We explored only a limited set of 145 well-studied HLA alleles but note that this analysis could be performed across a wider diversity of genotypes (48). Additionally, we did not assess genotypic heterogeneity or in vivo evolution of SARS-CoV-2, which could modify the repertoire of viral epitopes presented or could otherwise modulate virulence in an HLA-independent manner (64, 65) (https://nextstrain.org/ncov). We also did not address the potential for individual-level genetic variation in other proteins (e.g., angiotensin converting enzyme 2 [ACE2] or transmembrane serine protease 2 [TMPRSS2], essential host proteins for SARS-CoV-2 priming and cell entry [66]) to modulate the host-pathogen interface.

Unless and until the findings we present here are clinically validated, they should not be employed for any clinical purposes. However, we do at this juncture recommend integrating HLA testing into clinical trials and pairing HLA typing with COVID-19 testing where feasible to more rapidly develop and deploy a predictor(s) of viral severity in the population and, potentially, to tailor future vaccination strategies to genotypically at-risk populations. This approach may have additional implications for the management of a broad array of other viruses.

MATERIALS AND METHODS

Sequence retrieval and alignments.Full polyprotein 1ab (ORF1ab), spike (S) protein, membrane (M) protein, envelope (E) protein, and nucleocapsid (N) protein sequences were obtained for each of 34 distinct but representative alpha and betacoronaviruses from broad genus and subgenus distributions, including all known human coronaviruses (i.e., SARS-CoV, SARS-CoV-2, MERS-CoV, HKU1, OC43, NL63, and 229E). FASTA-formatted protein sequence data (the full accession number list is available in Table S5 in the supplemental material) were retrieved from the National Center of Biotechnology Information (NCBI) (67). For each of the protein classes (i.e., ORF1ab, S, M, E, and N), all 34 coronavirus sequences were aligned using the Clustal Omega v1.2.4 multisequence aligner tool employing the following parameters: sequence type [Protein], output alignment format [clustal_num], dealign [false], mBed-like clustering guide-tree [true], mBed-like clustering iteration [true], number of combined iterations 0, maximum guide tree iterations [-1], and maximum HMM iterations [-1] (68). For the purposes of estimating time of viral peptide production, we classified ORF1a and ORF1b peptides as “early” whereas all other peptides produced by subgenomic mRNAs were classified as “late” (69, 70).

Conserved peptide assessment.Aligned sequences were imported into Jalview v. 2.1.1 (71) with automated generation of the following alignment annotations: (i) sequence consensus, calculated as the percentage of the modal residue per column; (ii) sequence conservation (0 to 11), measured as a numerical index reflecting conservation of amino acid physicochemical properties in the alignment; (iii) alignment quality (0 to 1), measured as a normalized sum of BLOSUM62 ratios for all residues at each position; and (iv) occupancy, calculated as the number of aligned residues (not including gaps) for each position. In all cases, sequence conservation was assessed for each of the following three groups: only human-infecting coronavirus sequences (n = 7), all betacoronavirus sequences (n = 16), and all alpha- and betacoronavirus sequences combined (n = 34). Aligned SARS-CoV-2 sequences and all annotations were manually exported for subsequent analysis. Conserved human coronavirus peptides were defined as those with a length of ≥8 consecutive amino acids, each showing agreement with SARS-CoV-2 sequences and ≥4 other human coronavirus sequences with the consensus sequence (Table S2). For each of these conserved peptides, we also assessed the component number of 8- to 12-mers sharing identical amino acid sequence between SARS-CoV-2 and each of the four other common human coronaviruses (i.e., OC43, HKU1, NL63, and 229E) (Table S3). For all peptides, human, beta, and combined conservation scores were obtained using a custom R v.3.6.2 script representing mean sequence conservation (minus gap penalties where relevant) (see https://github.com/pdxgx/covid19).

Peptide-MHC class I binding affinity predictions.FASTA-formatted input protein sequences from the entire SARS-CoV-2 and SARS-CoV proteomes were obtained from the NCBI RefSeq database (67) under accession numbers NC_045512.2 and NC_004718.3. We kmerized each of these sequences into 8- to 12-mers to assess MHC class I-peptide binding affinity across the entire proteome. MHC class I binding affinity predictions were performed using 145 different HLA alleles for which global allele frequency data were available as described previously (72) (see Table S5) with netMHCpan v4.0 (73) using the ‘-BA’ option to include binding affinity predictions and the ‘-l’ option to specify peptides 8 to 12 amino acids in length (Table S1). Binding affinity was not predicted for peptides containing the character ‘|’ in their sequences. Additional MHC class I binding affinity predictions were performed on all 66 MHCflurry-supported alleles (–list-supported-alleles; Table S6) using both MHCnuggets 2.3.2 (74) and MHCflurry 1.4.3 (75) (see Tables S7, S8, and S9 and Fig. S7 to S10 in the supplemental material). We further cross-referenced these lists of peptides with existing experimentally validated SARS-CoV epitopes present in the Immune Epitope Database (Table S4) (76). We then performed consensus binding affinity predictions for the 66 supported alleles shared by all three tools by taking the union set of alleles and filtering for peptide-allele pairs matching the union set of alleles. For the SARS-CoV-specific and SARS-CoV-2-specific distributions of per-allele proteome presentation, we exclude all peptide-allele pairs with >500 nM predicted binding. In all cases, we used the netchop v3.0 (77) “C-term” model with a cleavage threshold of 0.1 to further remove any peptides that were not predicted to undergo canonical MHC class I antigen processing via proteasomal cleavage (of the peptide’s C terminus).

Global HLA allele and haplotype frequencies.HLA-A, -B, and -C allele and haplotype frequency data were obtained from the Allele Frequency Net Database (52) for 805 distinct populations pertaining to 101 different countries and 2,628 distinct major/minor (4-digit) alleles, corresponding to 20,478 distinct haplotypes (https://github.com/pdxgx/covid19). We also identified full HLA genotype data for 3,382 individuals whose HLA types were confined to the 145 HLA alleles studied here. Population allele and haplotype frequency data were aggregated by country as a mean of all constituent population allele or haplotype frequencies weighted by sample size of the population but not accounting for the representative ethnic demographic size of the population. Global allele frequency maps were generated using the rworldmap v1.3-6 package (78), with total global allele and haplotype frequency estimates calculated as the mean of per-country allele and haplotype frequencies, weighted by each country’s population in 2005.

Data availability.Source code is available at https://github.com/pdxgx/covid19 under the Massachusetts Institute of Technology (MIT) license. Data File S4 can be found at https://github.com/pdxgx/covid19/blob/master/supporting_data/Appendix_4.zip.

ACKNOWLEDGMENTS

We thank Christopher Loo and Jeffrey Barnet for their critical readings of the manuscript. We thank Jonah Sacha and Paul Spellman for their helpful discussions.

The contents do not necessarily represent the views of the U.S. Department of Veterans Affairs or the United States Government.

R.F.T. was supported by the U.S. Department of Veterans Affairs under award number 1IK2CX002049-01 and by the Sunlin & Priscilla Chou Foundation.

FOOTNOTES

    • Received 23 March 2020.
    • Accepted 14 April 2020.
    • Accepted manuscript posted online 17 April 2020.
  • Supplemental material is available online only.

  • Copyright © 2020 Nguyen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    World Health Organization. 2020. Statement on the second meeting of the International Health Regulations (2005) Emergency Committee regarding the outbreak of novel coronavirus (2019-nCoV). WHO (World Health Organization), Geneva, Switzerland.
  2. 2.↵
    1. Zhu N,
    2. Zhang D,
    3. Wang W,
    4. Li X,
    5. Yang B,
    6. Song J,
    7. Zhao X,
    8. Huang B,
    9. Shi W,
    10. Lu R,
    11. Niu P,
    12. Zhan F,
    13. Ma X,
    14. Wang D,
    15. Xu W,
    16. Wu G,
    17. Gao GF,
    18. Tan W
    ; China Novel Coronavirus Investigating and Research Team. 2020. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med 382:727–733. doi:10.1056/NEJMoa2001017.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Roser M,
    2. Ritchie HR,
    3. Ortiz-Ospina E,
    4. Hasell J
    . 2020. Coronavirus disease (COVID-19)—statistics and research. https://ourworldindata.org/coronavirus.
  4. 4.↵
    1. Caramelo F,
    2. Ferreira N,
    3. Oliveiros B
    . 2020. Estimation of risk factors for COVID-19 mortality-preliminary results. medRxiv doi:10.1101/2020.02.24.20027268.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Jain V,
    2. Yuan J-M
    . 2020. Systematic review and meta-analysis of predictive symptoms and comorbidities for severe COVID-19 infection. medRxiv doi:10.1101/2020.03.15.20035360.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Yang X,
    2. Yu Y,
    3. Xu J,
    4. Shu H,
    5. Xia J,
    6. Liu H,
    7. Wu Y,
    8. Zhang L,
    9. Yu Z,
    10. Fang M,
    11. Yu T,
    12. Wang Y,
    13. Pan S,
    14. Zou X,
    15. Yuan S,
    16. Shang Y
    . 24 February 2020, posting date. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. Lancet Respir Med doi:10.1016/S2213-2600(20)30079-5.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Wang D,
    2. Hu B,
    3. Hu C,
    4. Zhu F,
    5. Liu X,
    6. Zhang J,
    7. Wang B,
    8. Xiang H,
    9. Cheng Z,
    10. Xiong Y,
    11. Zhao Y,
    12. Li Y,
    13. Wang X,
    14. Peng Z
    . 2020. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan, China. JAMA 323:1061–1069. doi:10.1001/jama.2020.1585.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Zhou F,
    2. Yu T,
    3. Du R,
    4. Fan G,
    5. Liu Y,
    6. Liu Z,
    7. Xiang J,
    8. Wang Y,
    9. Song B,
    10. Gu X,
    11. Guan L,
    12. Wei Y,
    13. Li H,
    14. Wu X,
    15. Xu J,
    16. Tu S,
    17. Zhang Y,
    18. Chen H,
    19. Cao B
    . 2020. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet 395:1054–1062. doi:10.1016/S0140-6736(20)30566-3.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Guan W-J,
    2. Liang W-H,
    3. Zhao Y,
    4. Liang H-R,
    5. Chen Z-S,
    6. Li Y-M,
    7. Liu X-Q,
    8. Chen R-C,
    9. Tang C-L,
    10. Wang T
    , et al. 2020. Comorbidity and its impact on 1,590 patients with COVID-19 in China: a nationwide analysis. medRxiv doi:10.1101/2020.02.25.20027664.
    OpenUrlCrossRef
  10. 10.↵
    Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. 2020. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China. Zhonghua Liu Xing Bing Xue Za Zhi 41:145–151. (In Chinese.) doi:10.3760/cma.j.issn.0254-6450.2020.02.003.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Lau JTF,
    2. Lau M,
    3. Kim JH,
    4. Wong E,
    5. Tsui H-Y,
    6. Tsang T,
    7. Wong TW
    . 2004. Probable secondary infections in households of SARS patients in Hong Kong. Emerg Infect Dis 10:235–243. doi:10.3201/eid1002.030626.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    1. Denison MR
    . 2004. Severe acute respiratory syndrome coronavirus pathogenesis, disease and vaccines: an update. Pediatr Infect Dis J 23:S207–S214. doi:10.1097/01.inf.0000144666.95284.05.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Thabet F,
    2. Chehab M,
    3. Bafaqih H,
    4. AlMohaimeed S
    . 2015. Middle East respiratory syndrome coronavirus in children. Saudi Med J 36:484–343. doi:10.15537/smj.2015.4.10243.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Al-Tawfiq JA,
    2. Kattan RF,
    3. Memish ZA
    . 2016. Middle East respiratory syndrome coronavirus disease is rare in children: an update from Saudi Arabia. World J Clin Pediatr 5:391–396. doi:10.5409/wjcp.v5.i4.391.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Cao Q,
    2. Chen Y-C,
    3. Chen C-L,
    4. Chiu C-H
    . 2020. SARS-CoV-2 infection in children: transmission dynamics and clinical characteristics. J Formos Med Assoc 119:670–673. doi:10.1016/j.jfma.2020.02.009.
    OpenUrlCrossRef
  16. 16.↵
    1. Tang A,
    2. Xu W,
    3. Shen M,
    4. Chen P,
    5. Li G,
    6. Liu Y,
    7. Liu L
    . 2020. A retrospective study of the clinical characteristics of COVID-19 infection in 26 children. medRxiv doi:10.1101/2020.03.08.20029710.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Lu X,
    2. Zhang L,
    3. Du H,
    4. Zhang J,
    5. Li YY,
    6. Qu J,
    7. Zhang W,
    8. Wang Y,
    9. Bao S,
    10. Li Y,
    11. Wu C,
    12. Liu H,
    13. Liu D,
    14. Shao J,
    15. Peng X,
    16. Yang Y,
    17. Liu Z,
    18. Xiang Y,
    19. Zhang F,
    20. Silva RM,
    21. Pinkerton KE,
    22. Shen K,
    23. Xiao H,
    24. Xu S,
    25. Wong GWK
    , Chinese Pediatric Novel Coronavirus Study Team. 18 March 2020, posting date. SARS-CoV-2 infection in children. N Engl J Med doi:10.1056/NEJMc2005073.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Bi Q,
    2. Wu Y,
    3. Mei S,
    4. Ye C,
    5. Zou X,
    6. Zhang Z,
    7. Liu X,
    8. Wei L,
    9. Truelove SA,
    10. Zhang T,
    11. Gao W,
    12. Cheng C,
    13. Tang X,
    14. Wu X,
    15. Wu Y,
    16. Sun B,
    17. Huang S,
    18. Sun Y,
    19. Zhang J,
    20. Ma T,
    21. Lessler J,
    22. Feng T
    . 2020. Epidemiology and transmission of COVID-19 in Shenzhen China: analysis of 391 cases and 1,286 of their close contacts. medRxiv doi:10.1101/2020.03.03.20028423.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Dong Y,
    2. Mo X,
    3. Hu Y,
    4. Qi X,
    5. Jiang F,
    6. Jiang Z,
    7. Tong S
    . 16 March 2020, posting date Epidemiological characteristics of 2143 pediatric patients with 2019 coronavirus disease in China. Pediatrics doi:10.1542/peds.2020-0702.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Ashour HM,
    2. Elkhatib WF,
    3. Rahman MM,
    4. Elshabrawy HA
    . 2020. Insights into the recent 2019 novel coronavirus (SARS-CoV-2) in light of past human coronavirus outbreaks. Pathogens 9:E186. doi:10.3390/pathogens9030186.
    OpenUrlCrossRef
  21. 21.↵
    1. Yang Y,
    2. Peng F,
    3. Wang R,
    4. Guan K,
    5. Jiang T,
    6. Xu G,
    7. Sun J,
    8. Chang C
    . 3 March 2020, posting date. The deadly coronaviruses: the 2003 SARS pandemic and the 2020 novel coronavirus epidemic in China. J Autoimmun doi:10.1016/j.jaut.2020.102434.
    OpenUrlCrossRef
  22. 22.↵
    1. Ge Y,
    2. Tian T,
    3. Huang S,
    4. Wan F,
    5. Li J,
    6. Li S,
    7. Yang H,
    8. Hong L,
    9. Wu N,
    10. Yuan E,
    11. Cheng L,
    12. Lei Y,
    13. Shu H,
    14. Feng X,
    15. Jiang Z,
    16. Chi Y,
    17. Guo X,
    18. Cui L,
    19. Xiao L,
    20. Li Z,
    21. Yang C,
    22. Miao Z,
    23. Tang H,
    24. Chen L,
    25. Zeng H,
    26. Zhao D,
    27. Zhu F,
    28. Shen X,
    29. Zeng J
    . 2020. A data-driven drug repositioning framework discovered a potential therapeutic agent targeting COVID-19. bioRxiv doi:10.1101/2020.03.11.986836.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Ceraolo C,
    2. Giorgi FM
    . 2020. Genomic variance of the 2019-nCoV coronavirus. J Med Virol 92:522–528. doi:10.1002/jmv.25700.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Lu R,
    2. Zhao X,
    3. Li J,
    4. Niu P,
    5. Yang B,
    6. Wu H,
    7. Wang W,
    8. Song H,
    9. Huang B,
    10. Zhu N,
    11. Bi Y,
    12. Ma X,
    13. Zhan F,
    14. Wang L,
    15. Hu T,
    16. Zhou H,
    17. Hu Z,
    18. Zhou W,
    19. Zhao L,
    20. Chen J,
    21. Meng Y,
    22. Wang J,
    23. Lin Y,
    24. Yuan J,
    25. Xie Z,
    26. Ma J,
    27. Liu WJ,
    28. Wang D,
    29. Xu W,
    30. Holmes EC,
    31. Gao GF,
    32. Wu G,
    33. Chen W,
    34. Shi W,
    35. Tan W
    . 2020. Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet 395:565–574. doi:10.1016/S0140-6736(20)30251-8.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Prompetchara E,
    2. Ketloy C,
    3. Palaga T
    . 2020. Immune responses in COVID-19 and potential vaccines: lessons learned from SARS and MERS epidemic. Asian Pac J Allergy Immunol 38:1–9. doi:10.12932/AP-200220-0772.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Zheng Z,
    2. Monteil VM,
    3. Maurer-Stroh S,
    4. Yew CW,
    5. Leong C,
    6. Mohd-Ismail NK,
    7. Arularasu SC,
    8. Chow VTK,
    9. Pin RLT,
    10. Mirazimi A,
    11. Hong W,
    12. Tan Y-J
    . 2020. Monoclonal antibodies for the S2 subunit of spike of SARS-CoV cross-react with the newly-emerged SARS-CoV-2. bioRxiv doi:10.1101/2020.03.06.980037.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Yang Y,
    2. Shen C,
    3. Li J,
    4. Yuan J,
    5. Yang M,
    6. Wang F,
    7. Li G,
    8. Li Y,
    9. Xing L,
    10. Peng L,
    11. Wei J,
    12. Cao M,
    13. Zheng H,
    14. Wu W,
    15. Zou R,
    16. Li D,
    17. Xu Z,
    18. Wang H,
    19. Zhang M,
    20. Zhang Z,
    21. Liu L,
    22. Liu Y
    . 2020. Exuberant elevation of IP-10, MCP-3 and IL-1ra during SARS-CoV-2 infection is associated with disease severity and fatal outcome. medRxiv doi:10.1101/2020.03.02.20029975.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Ahmed SF,
    2. Quadeer AA,
    3. McKay MR
    . 2020. Preliminary identification of potential vaccine targets for the COVID-19 coronavirus (SARS-CoV-2) based on SARS-CoV immunological studies. Viruses 12:E254. doi:10.3390/v12030254.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Li G,
    2. Fan Y,
    3. Lai Y,
    4. Han T,
    5. Li Z,
    6. Zhou P,
    7. Pan P,
    8. Wang W,
    9. Hu D,
    10. Liu X,
    11. Zhang Q,
    12. Wu J
    . 2020. Coronavirus infections and immune responses. J Med Virol 92:424–432. doi:10.1002/jmv.25685.
    OpenUrlCrossRef
  30. 30.↵
    1. Lokugamage KG,
    2. Schindewolf C,
    3. Menachery VD
    . 2020. SARS-CoV-2 sensitive to type I interferon pretreatment. bioRxiv doi:10.1101/2020.03.07.982264.
    OpenUrlCrossRef
  31. 31.↵
    1. Wang C,
    2. Li W,
    3. Drabek D,
    4. Okba NMA,
    5. van Haperen R,
    6. Albert DM,
    7. van Kuppeveld FJM,
    8. Haagmans BL,
    9. Grosveld F,
    10. Bosch B-J
    . 2020. A human monoclonal antibody blocking SARS-CoV-2 infection. bioRxiv doi:10.1101/2020.03.11.987958.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Lv H,
    2. Wu NC,
    3. Tsang O-Y,
    4. Yuan M,
    5. Ranawaka AP,
    6. Leung WS,
    7. So RTY,
    8. Chan JMC,
    9. Yip GK,
    10. Chik TSH,
    11. Wang Y,
    12. Choi CYC,
    13. Lin Y,
    14. Ng WW,
    15. Zhao J,
    16. Poon LLM,
    17. Malik Peiris JS,
    18. Wilson IA,
    19. Mok CKP
    . 2020. Cross-reactive antibody response between SARS-CoV-2 and SARS-CoV infections. bioRxiv doi:10.1101/2020.03.15.993097.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Yuan M,
    2. Wu NC,
    3. Zhu X,
    4. Lee C-C,
    5. So RTY,
    6. Lv H,
    7. Mok CKP,
    8. Wilson IA
    . 3 April 2020, posting date. A highly conserved cryptic epitope in the receptor-binding domains of SARS-CoV-2 and SARS-CoV. Science doi:10.1126/science.abb7269.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Tan Y,
    2. Schneider T,
    3. Leong M,
    4. Aravind L,
    5. Zhang D
    . 2020. Novel immunoglobulin domain proteins provide insights into evolution and pathogenesis mechanisms of SARS-related coronaviruses. bioRxiv doi:10.1101/2020.03.04.977736.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Tetro JA
    . 2020. Is COVID-19 receiving ADE from other coronaviruses? Microbes Infect 22:72–73. doi:10.1016/j.micinf.2020.02.006.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Thevarajan I,
    2. Nguyen THO,
    3. Koutsakos M,
    4. Druce J,
    5. Caly L,
    6. van de Sand CE,
    7. Jia X,
    8. Nicholson S,
    9. Catton M,
    10. Cowie B,
    11. Tong S,
    12. Lewin S,
    13. Kedzierska K
    . 2020. Breadth of concomitant immune responses underpinning viral clearance and patient recovery in a non-severe case of COVID-19. medRxiv doi:10.1101/2020.02.20.20025841.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Zhu J,
    2. Kim J,
    3. Xiao X,
    4. Wang Y,
    5. Luo D,
    6. Chen R,
    7. Xu L,
    8. Zhang H,
    9. Xiao G,
    10. Zhan X,
    11. Wang T,
    12. Xie Y
    . 2020. Profiling the immune vulnerability landscape of the 2019 novel coronavirus. bioRxiv doi:10.1101/2020.02.08.939553.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Dong C,
    2. Ni L,
    3. Ye F,
    4. Chen M-L,
    5. Feng Y,
    6. Deng Y-Q,
    7. Zhao H,
    8. Wei P,
    9. Ge J,
    10. Li X,
    11. Sun L,
    12. Wang P,
    13. Liang P,
    14. Guo H,
    15. Wang X,
    16. Qin C-F,
    17. Chen F
    . 2020. Characterization of anti-viral immunity in recovered individuals infected by SARS-CoV-2. medRxiv doi:10.1101/2020.03.17.20036640.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Okba NMA,
    2. Muller MA,
    3. Li W,
    4. Wang C,
    5. GeurtsvanKessel CH,
    6. Corman VM,
    7. Lamers MM,
    8. Sikkema RS,
    9. de Bruin E,
    10. Chandler FD,
    11. Yazdanpanah Y,
    12. Le Hingrat Q,
    13. Descamps D,
    14. Houhou-Fidouh N,
    15. Reusken CBE,
    16. Bosch B-J,
    17. Drosten C,
    18. Koopmans MPG,
    19. Haagmans BL
    . 2020. SARS-CoV-2 specific antibody responses in COVID-19 patients. medRxiv doi:10.1101/2020.03.18.20038059.
    OpenUrlCrossRef
  40. 40.↵
    1. Chen G,
    2. Wu D,
    3. Guo W,
    4. Cao Y,
    5. Huang D,
    6. Wang H,
    7. Wang T,
    8. Zhang X,
    9. Chen H,
    10. Yu H,
    11. Zhang X,
    12. Zhang M,
    13. Wu S,
    14. Song J,
    15. Chen T,
    16. Han M,
    17. Li S,
    18. Luo X,
    19. Zhao J,
    20. Ning Q
    . 27 March 2020, posting date. Clinical and immunologic features in severe and moderate coronavirus disease 2019. J Clin Invest doi:10.1172/JCI137244.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Tahir Ul Qamar M,
    2. Rehman A,
    3. Ashfaq UA,
    4. Awan MQ,
    5. Fatima I,
    6. Shahid F,
    7. Chen L-L
    . 2020. Designing of a next generation multiepitope based vaccine (MEV) against SARS-COV-2: immunoinformatics and in silico approaches. bioRxiv doi:10.1101/2020.02.28.970343.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Fast E,
    2. Chen B
    . 2020. Potential T-cell and B-cell epitopes of 2019-nCoV. bioRxiv doi:10.1101/2020.02.19.955484.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Abdelmageed MI,
    2. Abdelmoneim AH,
    3. Mustafa MI,
    4. Elfadol NM,
    5. Murshed NS,
    6. Shantier SW,
    7. Makhawi AM
    . 2020. Design of multi epitope-based peptide vaccine against E protein of human COVID-19: an immunoinformatics approach. bioRxiv doi:10.1101/2020.02.04.934232.
    OpenUrlCrossRef
  44. 44.↵
    1. Baruah V,
    2. Bose S
    . 2020. Immunoinformatics-aided identification of T cell and B cell epitopes in the surface glycoprotein of 2019-nCoV. J Med Virol 92:495–500. doi:10.1002/jmv.25698.
    OpenUrlCrossRef
  45. 45.↵
    1. Robson B
    . 2020. Computers and viral diseases. Preliminary bioinformatics studies on the design of a synthetic vaccine and a preventative peptidomimetic antagonist against the SARS-CoV-2 (2019-nCoV, COVID-19) coronavirus. Comput Biol Med 119:103670. doi:10.1016/j.compbiomed.2020.103670.
    OpenUrlCrossRef
  46. 46.↵
    1. Dhama K,
    2. Sharun K,
    3. Tiwari R,
    4. Dadar M,
    5. Malik YS,
    6. Singh KP,
    7. Chaicumpa W
    , 18 March 2020, posting date. COVID-19, an emerging coronavirus infection: advances and prospects in designing and developing vaccines, immunotherapeutics, and therapeutics. Hum Vaccin Immunother doi:10.1080/21645515.2020.1735227.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Grifoni A,
    2. Sidney J,
    3. Zhang Y,
    4. Scheuermann RH,
    5. Peters B,
    6. Sette A
    . 16 March 2020, posting date. A sequence homology and bioinformatic approach can predict candidate targets for immune responses to SARS-CoV-2. Cell Host Microbe doi:10.1016/j.chom.2020.03.002.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Campbell KM,
    2. Steiner G,
    3. Wells DK,
    4. Ribas A,
    5. Kalbasi A
    . 2020. Prediction of SARS-CoV-2 epitopes across 9360 HLA class I alleles. bioRxiv doi:10.1101/2020.03.30.016931.
    OpenUrlCrossRef
  49. 49.↵
    1. Lin M,
    2. Tseng H-K,
    3. Trejaut JA,
    4. Lee H-L,
    5. Loo J-H,
    6. Chu C-C,
    7. Chen P-J,
    8. Su Y-W,
    9. Lim KH,
    10. Tsai Z-U,
    11. Lin R-Y,
    12. Lin R-S,
    13. Huang C-H
    . 2003. Association of HLA class I with severe acute respiratory syndrome coronavirus infection. BMC Med Genet 4:9. doi:10.1186/1471-2350-4-9.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. MacDonald KS,
    2. Fowke KR,
    3. Kimani J,
    4. Dunand VA,
    5. Nagelkerke NJ,
    6. Ball TB,
    7. Oyugi J,
    8. Njagi E,
    9. Gaur LK,
    10. Brunham RC,
    11. Wade J,
    12. Luscher MA,
    13. Krausa P,
    14. Rowland-Jones S,
    15. Ngugi E,
    16. Bwayo JJ,
    17. Plummer FA
    . 2000. Influence of HLA supertypes on susceptibility and resistance to human immunodeficiency virus type 1 infection. J Infect Dis 181:1581–1589. doi:10.1086/315472.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Stephens HAF,
    2. Klaythong R,
    3. Sirikong M,
    4. Vaughn DW,
    5. Green S,
    6. Kalayanarooj S,
    7. Endy TP,
    8. Libraty DH,
    9. Nisalak A,
    10. Innis BL,
    11. Rothman AL,
    12. Ennis FA,
    13. Chandanayingyong D
    . 2002. HLA-A and -B allele associations with secondary dengue virus infections correlate with disease severity and the infecting viral serotype in ethnic Thais. Tissue Antigens 60:309–318. doi:10.1034/j.1399-0039.2002.600405.x.
    OpenUrlCrossRefPubMedWeb of Science
  52. 52.↵
    1. González-Galarza FF,
    2. Takeshita LYC,
    3. Santos EJM,
    4. Kempson F,
    5. Maia MHT,
    6. da Silva ALS,
    7. Luiz Teles A,
    8. Ghattaoraya GS,
    9. Alfirevic A,
    10. Jones AR,
    11. Middleton D
    . 2015. Allele frequency net 2015 update: new features for HLA epitopes, KIR and disease and HLA adverse drug reaction associations. Nucleic Acids Res 43:D784–D788. doi:10.1093/nar/gku1166.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Sette A,
    2. Grey H,
    3. Oseroff C,
    4. Peters B,
    5. Moutaftsi M,
    6. Crotty S,
    7. Assarsson E,
    8. Greenbaum J,
    9. Kim Y,
    10. Kolla R,
    11. Tscharke D,
    12. Koelle D,
    13. Johnson RP,
    14. Blum J,
    15. Head S,
    16. Sidney J
    . 2009. Definition of epitopes and antigens recognized by vaccinia specific immune responses: their conservation in variola virus sequences, and use as a model system to study complex pathogens. Vaccine 27(Suppl 6):G21–G26. doi:10.1016/j.vaccine.2009.10.011.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Sibener LV,
    2. Fernandes RA,
    3. Kolawole EM,
    4. Carbone CB,
    5. Liu F,
    6. McAffee D,
    7. Birnbaum ME,
    8. Yang X,
    9. Su LF,
    10. Yu W,
    11. Dong S,
    12. Gee MH,
    13. Jude KM,
    14. Davis MM,
    15. Groves JT,
    16. Goddard WA,
    17. Heath JR,
    18. Evavold BD,
    19. Vale RD,
    20. Garcia KC
    . 2018. Isolation of a structural mechanism for uncoupling T cell receptor signaling from peptide-MHC binding. Cell 174:672–687.e27. doi:10.1016/j.cell.2018.06.017.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Gálvez J,
    2. Gálvez JJ,
    3. García-Peñarrubia P
    . 2019. Is TCR/pMHC affinity a good estimate of the T-cell response? An answer based on predictions from 12 phenotypic models. Front Immunol 10:349. doi:10.3389/fimmu.2019.00349.
    OpenUrlCrossRef
  56. 56.↵
    1. Zehn D,
    2. Lee SY,
    3. Bevan MJ
    . 2009. Complete but curtailed T-cell response to very low-affinity antigen. Nature 458:211–214. doi:10.1038/nature07657.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    1. Vatti A,
    2. Monsalve DM,
    3. Pacheco Y,
    4. Chang C,
    5. Anaya J-M,
    6. Gershwin ME
    . 2017. Original antigenic sin: a comprehensive review. J Autoimmun 83:12–21. doi:10.1016/j.jaut.2017.04.008.
    OpenUrlCrossRef
  58. 58.↵
    1. Park MS,
    2. Kim JI,
    3. Park S,
    4. Lee I,
    5. Park M-S
    . 2016. Original antigenic sin response to RNA viruses and antiviral immunity. Immune Netw 16:261–270. doi:10.4110/in.2016.16.5.261.
    OpenUrlCrossRef
  59. 59.↵
    1. Lambert P-H,
    2. Liu M,
    3. Siegrist C-A
    . 2005. Can successful vaccines teach us how to induce efficient protective immune responses? Nat Med 11:S54–S62. doi:10.1038/nm1216.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Singh RAK,
    2. Rodgers JR,
    3. Barry MA
    . 2002. The role of T cell antagonism and original antigenic sin in genetic immunization. J Immunol 169:6779–6786. doi:10.4049/jimmunol.169.12.6779.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Weiskopf D,
    2. Angelo MA,
    3. de Azeredo EL,
    4. Sidney J,
    5. Greenbaum JA,
    6. Fernando AN,
    7. Broadwater A,
    8. Kolla RV,
    9. De Silva AD,
    10. de Silva AM,
    11. Mattia KA,
    12. Doranz BJ,
    13. Grey HM,
    14. Shresta S,
    15. Peters B,
    16. Sette A
    . 2013. Comprehensive analysis of dengue virus-specific responses supports an HLA-linked protective role for CD8 T cells. Proc Natl Acad Sci U S A 110:E2046–E2053. doi:10.1073/pnas.1305227110.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Klenerman P,
    2. Zinkernagel RM
    . 1998. Original antigenic sin impairs cytotoxic T lymphocyte responses to viruses bearing variant epitopes. Nature 394:482–485. doi:10.1038/28860.
    OpenUrlCrossRefPubMed
  63. 63.↵
    1. Rothman AL
    . 2011. Immunity to dengue virus: a tale of original antigenic sin and tropical cytokine storms. Nat Rev Immunol 11:532–543. doi:10.1038/nri3014.
    OpenUrlCrossRefPubMed
  64. 64.↵
    1. Shen Z,
    2. Xiao Y,
    3. Kang L,
    4. Ma W,
    5. Shi L,
    6. Zhang L,
    7. Zhou Z,
    8. Yang J,
    9. Zhong J,
    10. Yang D,
    11. Guo L,
    12. Zhang G,
    13. Li H,
    14. Xu Y,
    15. Chen M,
    16. Gao Z,
    17. Wang J,
    18. Ren L,
    19. Li M
    . 4 March 2020, posting date. Genomic diversity of SARS-CoV-2 in coronavirus disease 2019 patients. Clin Infect Dis doi:10.1093/cid/ciaa203.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Tang X,
    2. Wu C,
    3. Li X,
    4. Song Y,
    5. Yao X,
    6. Wu X,
    7. Duan Y,
    8. Zhang H,
    9. Wang Y,
    10. Qian Z,
    11. Cui J,
    12. Lu J
    . 3 March 2020, posting date. On the origin and continuing evolution of SARS-CoV-2. Nat Sci Rev doi:10.1093/nsr/nwaa036.
    OpenUrlCrossRef
  66. 66.↵
    1. Hoffmann M,
    2. Kleine-Weber H,
    3. Schroeder S,
    4. Krüger N,
    5. Herrler T,
    6. Erichsen S,
    7. Schiergens TS,
    8. Herrler G,
    9. Wu N-H,
    10. Nitsche A,
    11. Müller MA,
    12. Drosten C,
    13. Pöhlmann S
    . 4 March 2020, posting date. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell doi:10.1016/j.cell.2020.02.052.
    OpenUrlCrossRefPubMed
  67. 67.↵
    1. O'Leary NA,
    2. Wright MW,
    3. Brister JR,
    4. Ciufo S,
    5. Haddad D,
    6. McVeigh R,
    7. Rajput B,
    8. Robbertse B,
    9. Smith-White B,
    10. Ako-Adjei D,
    11. Astashyn A,
    12. Badretdin A,
    13. Bao Y,
    14. Blinkova O,
    15. Brover V,
    16. Chetvernin V,
    17. Choi J,
    18. Cox E,
    19. Ermolaeva O,
    20. Farrell CM,
    21. Goldfarb T,
    22. Gupta T,
    23. Haft D,
    24. Hatcher E,
    25. Hlavina W,
    26. Joardar VS,
    27. Kodali VK,
    28. Li W,
    29. Maglott D,
    30. Masterson P,
    31. McGarvey KM,
    32. Murphy MR,
    33. O'Neill K,
    34. Pujar S,
    35. Rangwala SH,
    36. Rausch D,
    37. Riddick LD,
    38. Schoch C,
    39. Shkeda A,
    40. Storz SS,
    41. Sun H,
    42. Thibaud-Nissen F,
    43. Tolstoy I,
    44. Tully RE,
    45. Vatsan AR,
    46. Wallin C,
    47. Webb D,
    48. Wu W,
    49. Landrum MJ,
    50. Kimchi A, et al
    . 2016. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res 44:D733–D745. doi:10.1093/nar/gkv1189.
    OpenUrlCrossRefPubMed
  68. 68.↵
    1. Sievers F,
    2. Wilm A,
    3. Dineen D,
    4. Gibson TJ,
    5. Karplus K,
    6. Li W,
    7. Lopez R,
    8. McWilliam H,
    9. Remmert M,
    10. Söding J,
    11. Thompson JD,
    12. Higgins DG
    . 2011. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7:539. doi:10.1038/msb.2011.75.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Baranov PV,
    2. Henderson CM,
    3. Anderson CB,
    4. Gesteland RF,
    5. Atkins JF,
    6. Howard MT
    . 2005. Programmed ribosomal frameshifting in decoding the SARS-CoV genome. Virology 332:498–510. doi:10.1016/j.virol.2004.11.038.
    OpenUrlCrossRefPubMedWeb of Science
  70. 70.↵
    1. Taiaroa G,
    2. Rawlinson D,
    3. Featherstone L,
    4. Pitt M,
    5. Caly L,
    6. Druce J,
    7. Purcell D,
    8. Harty L,
    9. Tran T,
    10. Roberts J,
    11. Scott N,
    12. Catton M,
    13. Williamson D,
    14. Coin L,
    15. Duchene S
    . 2020. Direct RNA sequencing and early evolution of SARS-CoV-2. bioRxiv doi:10.1101/2020.03.05.976167.
    OpenUrlAbstract/FREE Full Text
  71. 71.↵
    1. Waterhouse AM,
    2. Procter JB,
    3. Martin DMA,
    4. Clamp M,
    5. Barton GJ
    . 2009. Jalview Version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics 25:1189–1191. doi:10.1093/bioinformatics/btp033.
    OpenUrlCrossRefPubMedWeb of Science
  72. 72.↵
    1. Wood MA,
    2. Paralkar M,
    3. Paralkar MP,
    4. Nguyen A,
    5. Struck AJ,
    6. Ellrott K,
    7. Margolin A,
    8. Nellore A,
    9. Thompson RF
    . 2018. Population-level distribution and putative immunogenicity of cancer neoepitopes. BMC Cancer 18:414. doi:10.1186/s12885-018-4325-6.
    OpenUrlCrossRef
  73. 73.↵
    1. Jurtz V,
    2. Paul S,
    3. Andreatta M,
    4. Marcatili P,
    5. Peters B,
    6. Nielsen M
    . 2017. NetMHCpan-4.0: improved peptide–MHC class I interaction predictions integrating eluted ligand and peptide binding affinity data. J Immunol 199:3360–3368. doi:10.4049/jimmunol.1700893.
    OpenUrlAbstract/FREE Full Text
  74. 74.↵
    1. Shao XM,
    2. Bhattacharya R,
    3. Huang J,
    4. Sivakumar IKA,
    5. Tokheim C,
    6. Zheng L,
    7. Hirsch D,
    8. Kaminow B,
    9. Omdahl A,
    10. Bonsack M,
    11. Riemer AB,
    12. Velculescu VE,
    13. Anagnostou V,
    14. Pagel KA,
    15. Karchin R
    . 2020. High-throughput prediction of MHC class I and II neoantigens with MHCnuggets. Cancer Immunol Res 8:396–408. doi:10.1158/2326-6066.CIR-19-0464.
    OpenUrlAbstract/FREE Full Text
  75. 75.↵
    1. O'Donnell TJ,
    2. Rubinsteyn A,
    3. Bonsack M,
    4. Riemer AB,
    5. Laserson U,
    6. Hammerbacher J
    . 2018. MHCflurry: open-source class I MHC binding affinity prediction. Cell Syst 7:129–132.e4. doi:10.1016/j.cels.2018.05.014.
    OpenUrlCrossRef
  76. 76.↵
    1. Vita R,
    2. Mahajan S,
    3. Overton JA,
    4. Dhanda SK,
    5. Martini S,
    6. Cantrell JR,
    7. Wheeler DK,
    8. Sette A,
    9. Peters B
    . 2019. The Immune Epitope Database (IEDB): 2018 update. Nucleic Acids Res 47:D339–D343. doi:10.1093/nar/gky1006.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Nielsen M,
    2. Lundegaard C,
    3. Lund O,
    4. Keşmir C
    . 2005. The role of the proteasome in generating cytotoxic T-cell epitopes: insights obtained from improved predictions of proteasomal cleavage. Immunogenetics 57:33–41. doi:10.1007/s00251-005-0781-7.
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.↵
    1. South A
    . 2011. rworldmap: a new R package for mapping global data. R J 3:35–43. doi:10.32614/RJ-2011-006.
    OpenUrlCrossRef
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Human Leukocyte Antigen Susceptibility Map for Severe Acute Respiratory Syndrome Coronavirus 2
Austin Nguyen, Julianne K. David, Sean K. Maden, Mary A. Wood, Benjamin R. Weeder, Abhinav Nellore, Reid F. Thompson
Journal of Virology Jun 2020, 94 (13) e00510-20; DOI: 10.1128/JVI.00510-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Journal of Virology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Human Leukocyte Antigen Susceptibility Map for Severe Acute Respiratory Syndrome Coronavirus 2
(Your Name) has forwarded a page to you from Journal of Virology
(Your Name) thought you would be interested in this article in Journal of Virology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Human Leukocyte Antigen Susceptibility Map for Severe Acute Respiratory Syndrome Coronavirus 2
Austin Nguyen, Julianne K. David, Sean K. Maden, Mary A. Wood, Benjamin R. Weeder, Abhinav Nellore, Reid F. Thompson
Journal of Virology Jun 2020, 94 (13) e00510-20; DOI: 10.1128/JVI.00510-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

COVID-19
HLA
MHC class I
SARS-CoV-2
coronavirus

Related Articles

Cited By...

About

  • About JVI
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #Jvirology

@ASMicrobiology

       

 

JVI in collaboration with

American Society for Virology

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Print ISSN: 0022-538X; Online ISSN: 1098-5514