Next Article in Journal
Diagnostic Performances of Different Genome Amplification Assays for the Detection of Swine Vesicular Disease Virus in Relation to Genomic Lineages That Circulated in Italy
Next Article in Special Issue
Demystifying Excess Immune Response in COVID-19 to Reposition an Orphan Drug for Down-Regulation of NF-κB: A Systematic Review
Previous Article in Journal
Neurological Findings in Children without Congenital Microcephaly Exposed to Zika Virus in Utero: A Case Series Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 11
10.3390/v12111333
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Population Difference in Allele Frequency of HLA-C*05 and Its Correlation with COVID-19 Mortality

by
Atsushi Sakuraba
1,*,
Haider Haider
1 and
Toshiro Sato
2,3
1
Section of Gastroenterology, Hepatology, and Nutrition, University of Chicago Medicine, Chicago, IL 60637, USA
2
Department of Organoid Medicine, Keio University School of Medicine, Tokyo 160-8582, Japan
3
Corona Virus Task Force, Keio University School of Medicine, Tokyo 160-8582, Japan
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(11), 1333; https://doi.org/10.3390/v12111333
Submission received: 30 October 2020 / Revised: 13 November 2020 / Accepted: 19 November 2020 / Published: 20 November 2020
(This article belongs to the Special Issue Virus-Induced Cytokine Storms)
Browse Figures
Versions Notes

Abstract

:
Background: coronavirus disease 2019 (COVID-19) causes severe illness including cytokine storms, but mortality among countries differs largely. In the present study, we investigated the association between human leukocyte antigen (HLA) class I, which plays a major role in susceptibility to viral infections, and the mortality of COVID-19. Methods: data of allele frequencies of HLA-A, -B and -C and COVID-19 mortality were obtained for 74 countries from the Allele Frequency Net Database and worldometer.info. Association between allele frequency of each HLA and mortality was assessed by linear regression followed by multivariable regression. Subsequently, association of HLA-C*05 to its receptor KIR2DS4fl, expressed on natural killer (NK) cells, and differential mortality to historic pandemics were analyzed. Results: HLA-A*01, -B*07, -B*08, -B*44 and -C*05 were significantly associated with the risk of deaths (adjusted p = 0.040, 0.00081, 0.047, 0.0022, 0.00032, respectively), but only HLA-C*05 remained statistically significant (p = 0.000027) after multivariable regression. A 1% increase in the allele frequency of HLA-C*05 was associated with an increase of 44 deaths/million. Countries with different mortality could be categorized by the distribution of HLA-C*05 and its receptor KIR2DS4fl, which in combination cause NK cell-induced hyperactive immune response. Countries with similar ethnic and/or geographic background responded in a similar pattern to each pandemic. Conclusions: we demonstrated that allele frequency of HLA-C*05 and the distribution pattern with its receptor KIR2DS4fl strongly correlated with COVID-19 mortality. Host genetic variance of innate immunity may contribute to the difference in mortality among various countries and further investigation using patient samples is warranted.

1. Introduction

The symptoms of coronavirus disease 2019 (COVID-19), which is the acute respiratory syndrome caused by the virus, and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), vary from asymptomatic to critically ill [1]. It has recently come to attention that the number of infections and deaths among countries worldwide differ largely. In China, where COVID-19 was initially identified, the infection and deaths are 58 and 3 per million population, respectively. Most Asian countries show rates in the range similar to China. In many European countries, these numbers are much greater; e.g., the United Kingdom (UK), Italy and Spain, have rates in the range of 3000–6000 and 300–800, respectively. The difference in the infection rate can be partly explained by demographic background, healthcare and governmental policies, and approach to diagnostic testing and containment of infection. However, host factors contributing to the difference in mortality due to COVID-19 remain largely unknown.
Elderly people are known to be more prone to the risks of COVID-19. The case fatality risk was reported to be 1.38% for all cases, but was substantially greater in patients ≥60 as compared to <60 (6.4% vs. 0.32%) [2]. Other risk factors include male sex, smoking, diabetes and obesity. In the United States of America (USA) and the UK, minority ethnicity groups had a higher incidence of deaths [3]. While structural inequalities contribute to difference in deaths between races in these countries, host genetic differences in susceptibility to, especially severe, COVID-19 among populations worldwide likely exist.
A majority of patients admitted for severe COVID-19 experience complications such as thromboembolism, multiorgan failure and acute respiratory distress syndrome (ARDS). Elevated levels of inflammatory cytokines such as interleukin (IL)-2R, macrophage inflammatory protein and IL-6, characteristic of cytokine storm syndrome have been found in these patients [4,5,6]. This has rationalized the use of immunomodulatory agents, tocilizumab and dexamethasone, in severe or critically ill COVID-19 patients [4,7].
Major histocompatibility complex (MHC) class I molecules are found on the cell surface of all nucleated cells. The main function of MHC class I molecules are to display intracellular peptide fragments to cytotoxic T cells and natural killer (NK) cells, which in turn induce an immune response [8]. Virus-infected host cells utilize this mechanism to trigger immune surveillance and eliminate infected cells, but the immune system can cause autoimmunity or overproduction of cytokines when regulation is impaired [9]. In humans, there are three MHC class I genes, namely human leukocyte antigen (HLA)-A, HLA-B and HLA-C. HLA loci are known as a leading genetic candidate for infectious disease susceptibility. Studies have shown that patients with certain HLA alleles are at risk for viral infections including SARS [10]. The frequency and composition of HLA among different human populations vary largely which is explained by diversifying selection, balancing selection and pathogen-driven selection [11,12]. In the present study, we investigated the association between the allele frequency of class I MHC genes, HLA-A, -B and -C, obtained from the Allele Frequency Net Database and the worldwide mortality due to COVID-19.

2. Materials and Methods

The data of allele frequencies of HLA-A, -B and -C genes among countries were obtained from the Allele Frequency Net Database (http://allelefrequencies.net) [13] and for mortality (deaths per one million population) were obtained from Worldometer (https://www.worldometers.info). Worldometer collects data from individual countries’ communication channels and provides global COVID-19 live statistics. It has been widely used in literature [14]. Both databases were searched on 28 May 2020. For the purpose of this study, we limited the search at the allele group level (first field). There were 79 countries with data of allele frequency of at least one HLA-A, B or C. When multiple datapoints of an allele frequency were available for a country, the weighted mean was calculated according to the sample size. Data for death/1 million population due to COVID-19 were available for 213 countries. The final dataset included 74 countries with data available for both HLA and deaths.
The analysis of the association between the allele frequency of class I MHC genes and the mortality due to COVID-19 demonstrated that HLA-C*05:01 had the strongest association. Literature has demonstrated that peptides bound to HLA-C*05:01 are recognized by an activating killer cell immunoglobulin-like receptor (KIR), KIR2DS4 [8,15]. Thus, we studied the association between the carrier frequency of the functional form of KIR2DS4, KIR2DS4fl and HLA-C*05. Data of frequency of carriers of KIR2DS4fl were obtained from the Allele Frequency Net Database [13] and also by a search of the published literature [16,17,18,19,20].
Based on the results of HLA-C*05/KIR2DS4 association, we compared the mortality due to COVID-19 among countries to those of historic global pandemics, the 1918–1920 Spanish flu and 2009 swine flu [21,22], to examine whether there is a pattern suggestive of pathogen-derived epistatic balancing selection [23].

Statistical Analysis

The association between the allele frequency of each HLA gene and mortality was assessed by linear regression. Multiplicity adjusted p-values following Bonferroni multiple comparison testing was done. An adjusted p-value of <0.05 was considered statistically significant. Multivariate regression with genes with a p-value of <0.05 was performed. We then applied backward elimination with a threshold of 0.05 for selection of the most predictive variable(s). Differences were considered to be statistically significant when the p-value was <0.05. Data were analyzed by EZR (Saitama Medical Center, Jichi Medical University, Saitama, Japan) [24], which is a graphical user interface for R (The R Foundation for Statistical Computing, version 2.13.0, Vienna, Austria), and GraphPad Prism version 6 (San Diego, CA, USA).

3. Results

Data of 21, 34 and 13 HLA-A, -B and -C allele frequencies, respectively, as well as cases of death due to COVID-19 were available for 74 countries (Figure S1). Countries were distributed evenly among continents. There were 10, 25, 25, 6 and 8 countries from Africa, Asia/Oceania, Europe, North/Central America and South America, respectively.

3.1. Association between Worldwide Allele Frequency of HLAs and Mortality to COVID-19

The association of each allele frequency and number of deaths due to COVID-19 per one million population was analyzed by linear regression for all HLA alleles (Table 1). After adjusting for multiplicity, HLA-A*01, -B*07, -B*08, -B*44 and -C*05 were significantly associated with the risk of deaths (adjusted p = 0.040, 0.00081, 0.047, 0.0022, 0.00032, respectively). Multivariate regression with backward elimination demonstrated that HLA-C*05 was the only variable that was significantly associated with the risk of deaths (p = 0.000027). Figure 1 shows the linear regression model of allele frequency of HLA-C*05 and the risk of deaths per million population (adjusted R2 0.37, p = 0.0000047). It is estimated that for each 1% increase in allele frequency of HLA-C*05, mortality increased by 44 deaths per 1 million population.

3.2. Association between HLA-C*05 and Its Receptor, KIR2DS4

A recent report found that peptides bound to HLA-C*05:01 are recognized by one of the activating KIR, KIR2DS4 [8]. There are two major alleles in KIR2DS4, the full-length receptor, KIR2DS4-fl, and KIR2DS4-del that harbors a 22-bp deletion and loses membrane-anchoring [25]. KIR2DS4-fl serves as a potent activating KIR, while KIR2DS4-del is devoid of NK cell activation potential [8]. Intriguingly, KIR2DS4-del is highly associated with HLA-C*05, and both alleles are common in European populations [8], suggesting that the effective NK cell activation by these ligand-receptor pairs drove the epistasis evolution [8]. We obtained data of the frequency of carriers of KIR2DS4fl among 16 countries and plotted them against the allele frequency of HLA-C*05 and deaths per 1 million population (Figure 2). Countries could be separated into HLA-C*05 high/KIR2DS4fl medium countries that exhibited high mortality (France, Italy and Spain), HLA-C*05 medium/KIR2DS4fl medium countries that demonstrated medium mortality (Argentina, Brazil, Croatia, Cuba, Iran, Poland and USA) and HLA-C*05 low/KIR2DS4fl high countries that demonstrated low mortality (India, Japan, Mexico, Saudi Arabia, Senegal and Tunisia).

3.3. Distinct Pattern of Response to Historic and COVID-19 Pandemics

Based upon the notion that historic pandemics have posed evolutional selection pressure on HLA alleles [26], we next compared the mortality due to COVID-19 to those of two major historic pandemics, the 1918–1920 Spanish flu and 2009 swine flu [21]. With the exception of a few outliers, countries with similar ethnic and/or geographic background showed similar degree of mortality to each pandemic (Figure 3). European countries had medium mortality to Spanish flu and swine flu, but were categorized into two groups based on COVID-19 high (Spain, Italy, France, etc.) and medium (Germany, Austria, Switzerland, etc.) mortality. Most countries in the American continent (USA, Canada, Mexico, Brazil, etc.) were clustered as Spanish flu medium/COVID-19 medium/swine flu high mortality. African countries (Kenya, Ghana, Nigeria, etc.) were categorized as Spanish flu high/COVID-19 low/swine flu low. A majority of Asian countries (Japan, Philippines, China, etc.) were clustered as low mortality to all three pandemics. Noteworthy, categorization of countries was near identical to that of HLA-C*05/KIR2DS4fl. Countries with similar ethnic and/or geographic background responded in a similar pattern to each global pandemic indicative of evolutionary epistatic and balancing selection precipitated by these pandemics.

4. Discussion

In the present study, we investigated whether a host genetic factor is associated with the risk of death due to COVID-19. Using the allele frequency data of HLA-A, -B and -C in 74 countries, we have identified five HLAs, namely HLA-A*01, -B*07, -B*08, -B*44 and -C*05, as potential candidates that increase mortality. Amongst them, HLA-C*05 was the most influential in increasing the risk of death due to COVID-19. Our results suggest that the distribution of certain HLAs at the population level account for the disparate susceptibility and mortality to COVID-19 across the world.
The HLA system harbors genes that are crucial in modulating the immune system against pathogens. HLA alleles or haplotypes are related with susceptibility to infections caused by various pathogens [10]. The frequency of HLA alleles varies among different populations or races and has been reported to influence the susceptibility to viral respiratory infections. SARS coronavirus infection, which caused an outbreak in 2003, was extremely contagious and had a high rate of mortality [27]. Taiwanese studies demonstrated associations of the HLA-B*46:01 allele with severity and the HLA-C*15:02 and -DRB1*03:01 alleles with resistance against SARS [10], however, this was not reproduced in another cohort [28]. Keicho et al. showed that the HLA-DRB1*12:02 allele was found at a higher incidence in SARS patients than in controls [29]. Furthermore, the HLA phenotype has been shown to control immune reactivity against the influenza virus [30]. Wang et al. showed that HLA-C*07:29 and B*15:2 were found more frequently in Chinese COVID-19 patients than in controls, however, their study included only 99 patients [31].
Our study has identified HLA-C*05 as the allele most strongly correlating with the risk of deaths due to COVID-19 at a global level. HLA-C molecules load nonamer peptides and function as ligands for the KIR [32]. KIR are expressed on natural killer (NK) cells and T-cells, and classified as activating or inhibitory type by the presence of inhibitory signal domain in their cytoplasmic tail [33,34]. Sim et al. reported that KIR2DS4, expressed on NK cells, detects a conserved bacterial epitope presented by HLA-C*05:01 and exerts a potent activation signal [8]. KIR2DS4 are known to recognize HLA-C-bound viral peptides, so KIR2DS4+ NK cells might also recognize viral peptides presented by HLA-C*05 [8]. The inverse correlation between the frequency of HLA-C*05 and the KIR2DS4fl alleles raised a hypothesis that this HLA-KIR pair induces immune overactivation and subsequently causes negative selection. We found that populations with the highest mortalities (France, Italy and Spain) frequently carried HLA-C*05 and KIR2DS4fl, while countries with lower mortality had HLA-C*05:01 low-medium/KIR2DS4fl high. Thus, patients with a HLA-C*05 and KIR2DS4fl pair, which is most frequently seen in high mortality countries, may be prone to an excess cytokine response and suffer from hypercytokinemia, characteristic of elderly patients with severe COVID-19. Of note, COVID-19 is associated with Kawasaki disease like multisystem inflammatory syndrome in children in Europe and USA, suggesting that hypercytokinemia is a hallmark of severe COVID-19 infection regardless of age [35]. Cytokine storm and HLA–KIR interaction have also been observed with severe influenza virus infections [36] and with Spanish flu which was most fatal in young adults [37]. Furthermore, by analyzing mortality to historic pandemics, we found that countries have distinct patterns of response to each pandemic. Meanwhile, HLA–KIR interaction and their selection through human evolution has been associated with the risk of preeclampsia and reproductive outcome [38,39]. These results suggested that population-level HLA-KIR diversity undergoes epistasis and pathogen-mediated selection and strongly influences the mortality to pandemics.
Limitations of our study include that we used database, but not patient-level samples for assessment of HLA. The Allele Frequency Net database has been widely used to compare immune gene frequencies among different countries, but they may not accurately represent the population of a country. The spread and casualty of SARS-CoV-2 within a country can be assessed by several methods, but we chose deaths as our outcome as it is thought that most countries would universally test a hospitalized patient with interstitial pneumonitis to rule out COVID-19. However, testing criteria for COVID-19 vary and there may be underreporting especially in hard-hit countries during the peak of infection. Furthermore, there may be other factors that may influence the spread or outcome of COVID-19 among countries, such as mask-compliance rates, lockdown policies, healthcare provisions and ICU beds per capita. We are aware that our analysis is preliminary as the COVID-19 pandemic is still progressing and the infection may not have fully penetrated in some countries. Furthermore, there is a risk for second and third waves as we have seen in past viral pandemics and according to observance of countries/areas that lifted lockdowns. While our results suggest the presence of a correlation between HLA-C*05 and death due to COVID-19, prospective genetic studies will be required to study causal relationships [40]. During the preparation of this paper, the first genome-wide association study (GWAS) study of severe COVID-19 with respiratory failure was reported [41]. This study failed to show significant allele association of classical HLA loci, however, they only used genome data from Spanish and Italian populations and used blood donors with unknown COVID-19 status as controls, which may have precluded the detection of weaker correlations among classical HLA loci. On the other hand, a preprint study of Chinese patients with COVID-19 reported that TMEM189-UBE2V1, involved in the IL-1 signaling pathway, and the HLA-A*11:01, -B*51:01 and -C*14:02 alleles were significantly associated with severe disease or worse outcome [42]. This suggests that there may be other risk and/or protective HLA(s) or HLA haplotypes among different races.
Ever since the beginning of the COVID-19 pandemic, researchers and the community have been perplexed by the large variation in the fatality seen among countries. We studied the association between class I MHC, HLA-A, -B and -C, and the risk of deaths due to COVID-19 at a worldwide level and have identified HLA-C*05 as a potential candidate that may influence mortality. The results of our study demonstrated a correlation, so additional research of the functional properties of HLA-C*05 and a subsequent study to assess causal relationship are warranted. If our results are upheld, then the situation will represent a substantial opportunity for modifying the strategy to protect against or treat COVID-19 amongst countries with varying degree of risks.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/11/1333/s1, Figure S1: Allele frequencies of class I HLAs among countries in each continent.

Author Contributions

A.S., conceptualization; methodology; writing and editing. H.H., analysis of data; approval of final draft. T.S., conceptualization (equal); editing and approval of final draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

All authors have no conflicts of interest directly relevant to the content of this article.

References

  1. Wang, F.; Hou, H.; Luo, Y.; Tang, G.; Wu, S.; Huang, M.; Liu, W.; Zhu, Y.; Lin, Q.; Mao, L.; et al. The laboratory tests and host immunity of COVID-19 patients with different severity of illness. JCI Insight 2020, 5. [Google Scholar] [CrossRef]
  2. Verity, R.; Okell, L.C.; Dorigatti, I.; Winskill, P.; Whittaker, C.; Imai, N.; Cuomo-Dannenburg, G.; Thompson, H.; Walker, P.G.T.; Fu, H.; et al. Estimates of the severity of coronavirus disease 2019: A model-based analysis. Lancet Infect. Dis. 2020, 20, 669–677. [Google Scholar] [CrossRef]
  3. Anyane-Yeboa, A.; Sato, T.; Sakuraba, A. Racial disparities in COVID-19 deaths reveal harsh truths about structural inequality in America. J. Intern. Med. 2020. [Google Scholar] [CrossRef]
  4. Zhao, M. Cytokine storm and immunomodulatory therapy in COVID-19: Role of chloroquine and anti-IL-6 monoclonal antibodies. Int. J. Antimicrob. Agents 2020, 55, 105982. [Google Scholar] [CrossRef]
  5. Vabret, N.; Britton, G.J.; Gruber, C.; Hegde, S.; Kim, J.; Kuksin, M.; Levantovsky, R.; Malle, L.; Moreira, A.; Park, M.D.; et al. Immunology of COVID-19: Current State of the Science. Immunity 2020, 52, 910–941. [Google Scholar] [CrossRef]
  6. Mangalmurti, N.; Hunter, C.A. Cytokine Storms: Understanding COVID-19. Immunity 2020, 53, 19–25. [Google Scholar] [CrossRef]
  7. Toniati, P.; Piva, S.; Cattalini, M.; Garrafa, E.; Regola, F.; Castelli, F.; Franceschini, F.; Airò, P.; Bazzani, C.; Beindorf, E.-A.; et al. Tocilizumab for the treatment of severe COVID-19 pneumonia with hyperinflammatory syndrome and acute respiratory failure: A single center study of 100 patients in Brescia, Italy. Autoimmun. Rev. 2020, 19, 102568. [Google Scholar] [CrossRef]
  8. Sim, M.J.W.; Rajagopalan, S.; Altmann, D.M.; Boyton, R.J.; Sun, P.D.; Long, E.O. Human NK cell receptor KIR2DS4 detects a conserved bacterial epitope presented by HLA-C. Proc. Natl. Acad. Sci. USA 2019, 116, 12964–12973. [Google Scholar] [CrossRef] [Green Version]
  9. Coperchini, F.; Chiovato, L.; Croce, L.; Magri, F.; Rotondi, M. The cytokine storm in COVID-19: An overview of the involvement of the chemokine/chemokine-receptor system. Cytokine Growth Factor Rev. 2020, 53, 25–32. [Google Scholar] [CrossRef]
  10. Sun, Y.; Xi, Y. Association between HLA Gene Polymorphism and the Genetic Susceptibility of SARS Infection. In HLA and Associated Important Diseases; InTech: London, UK, 2014. [Google Scholar]
  11. Sanchez-Mazas, A. A review of HLA allele and SNP associations with highly prevalent infectious diseases in human populations. Swiss Med. Wkly. 2020, 150, w20214. [Google Scholar]
  12. Prugnolle, F.; Manica, A.; Charpentier, M.; Guégan, J.F.; Guernier, V.; Balloux, F. Pathogen-Driven Selection and Worldwide HLA Class I Diversity. Curr. Biol. 2005, 15, 1022–1027. [Google Scholar] [CrossRef] [Green Version]
  13. Gonzalez-Galarza, F.F.; McCabe, A.; Santos, E.; Jones, J.; Takeshita, L.; Ortega-Rivera, N.D.; Cid-Pavon, G.M.D.; Ramsbottom, K.; Ghattaoraya, G.; Alfirevic, A.; et al. Allele frequency net database (AFND) 2020 update: Gold-standard data classification, open access genotype data and new query tools. Nucleic Acids Res. 2020, 48, D783–D788. [Google Scholar] [CrossRef]
  14. Worldometer COVID-19 Coronavirus Pandemic. Available online: https://www.worldometers.info/coronavirus/ (accessed on 2 October 2020).
  15. Goodson-Gregg, F.J.; Krepel, S.A.; Anderson, S.K. Tuning of human NK cells by endogenous HLA-C expression. Immunogenetics 2020, 72, 205–215. [Google Scholar] [CrossRef] [Green Version]
  16. Bontadini, A.; Testi, M.; Cuccia, M.; Martinetti, M.; Carcassi, C.; Chiesa, A.; Cosentini, E.; Dametto, E.; Frison, S.; Iannone, A.; et al. Distribution of killer cell immunoglobulin-like receptors genes in the Italian Caucasian population. J. Transl. Med. 2006, 4, 44. [Google Scholar] [CrossRef] [Green Version]
  17. Tajik, N.; Shahsavar, F.; Mousavi, T.; Radjabzadeh, M.F. Distribution ofKIRgenes in the Iranian population. Tissue Antigens 2009, 74, 22–31. [Google Scholar] [CrossRef]
  18. Yawata, M.; Yawata, N.; McQueen, K.; Cheng, N.; Guethlein, L.; Rajalingam, R.; Shilling, H.; Parham, P. Predominance of group A KIR haplotypes in Japanese associated with diverse NK cell repertoires of KIR expression. Immunogenetics 2002, 54, 543–550. [Google Scholar] [CrossRef]
  19. Jiang, K.; Zhu, F.M.; Lv, Q.F.; Yan, L.X. Distribution of killer cell immunoglobulin-like receptor genes in the Chinese Han population. Tissue Antigens 2005, 65, 556–563. [Google Scholar] [CrossRef]
  20. Whang, D.H.; Park, H.; Yoon, J.A.; Park, M.H. Haplotype analysis of killer cell immunoglobulin-like receptor genes in 77 Korean families. Hum. Immunol. 2005, 66, 146–154. [Google Scholar] [CrossRef]
  21. Johnson, N.P.; Mueller, J. Updating the accounts: Global mortality of the 1918–1920 “Spanish” influenza pandemic. Bull. Hist. Med. 2002, 76, 105–115. [Google Scholar] [CrossRef]
  22. 2009 Swine Flu Pandemic by Country. Available online: https://en.wikipedia.org/wiki/2009_swine_flu_pandemic_by_country (accessed on 18 May 2020).
  23. Sabouri Ghannad, M.; Hajilooi, M.; Solgi, G. HLA-KIR Interactions and Immunity to Viral Infections. Res. Mol. Med. 2014, 2, 1–20. [Google Scholar] [CrossRef]
  24. Kanda, Y. Investigation of the freely available easy-to-use software ‘EZR’ for medical statistics. Bone Marrow Transpl. 2013, 48, 452–458. [Google Scholar] [CrossRef] [Green Version]
  25. Graef, T.; Moesta, A.K.; Norman, P.J.; Abi-Rached, L.; Vago, L.; Older Aguilar, A.M.; Gleimer, M.; Hammond, J.A.; Guethlein, L.A.; Bushnell, D.A.; et al. KIR2DS4 is a product of gene conversion with KIR3DL2 that introduced specificity for HLA-A*11 while diminishing avidity for HLA-C. J. Exp. Med. 2009, 206, 2557–2572. [Google Scholar] [CrossRef] [Green Version]
  26. Barreiro, L.B.; Quintana-Murci, L. From evolutionary genetics to human immunology: How selection shapes host defence genes. Nat. Rev. Genet. 2010, 11, 17–30. [Google Scholar] [CrossRef]
  27. Peiris, J.S.M.; Yuen, K.Y.; Osterhaus, A.D.M.E.; Stöhr, K. The Severe Acute Respiratory Syndrome. N. Engl. J. Med. 2003, 349, 2431–2441. [Google Scholar] [CrossRef] [Green Version]
  28. Yuan, F.F.; Velickovic, Z.; Ashton, L.J.; Dyer, W.B.; Geczy, A.F.; Dunckley, H.; Lynch, G.W.; Sullivan, J.S. Influence of HLA gene polymorphisms on susceptibility and outcome post infection with the SARS-CoV virus. Virol. Sin. 2014, 29, 128–130. [Google Scholar] [CrossRef] [Green Version]
  29. Keicho, N.; Itoyama, S.; Kashiwase, K.; Phi, N.C.; Long, H.T.; Ha, L.D.; Ban, V.V.; Hoa, B.K.; Hang, N.T.L.; Hijikata, M.; et al. Association of human leukocyte antigen class II alleles with severe acute respiratory syndrome in the Vietnamese population. Hum. Immunol. 2009, 70, 527–531. [Google Scholar] [CrossRef]
  30. Tussey, L.G.; Rowland-Jones, S.; Zheng, T.S.; Androlewicz, M.J.; Cresswell, P.; Frelinger, J.A.; McMichael, A.J. Different MHC class I alleles compete for presentation of overlapping viral epitopes. Immunity 1995, 3, 65–77. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, W.; Zhang, W.; Zhang, J.; He, J.; Zhu, F. Distribution of HLA allele frequencies in 82 Chinese individuals with coronavirus disease-2019 (COVID-19). HLA 2020. [Google Scholar] [CrossRef]
  32. Parham, P. Immunogenetics of killer-cell immunoglobulin-like receptors. Tissue Antigens 2003, 62, 194–200. [Google Scholar] [CrossRef]
  33. Yawata, M.; Yawata, N.; Draghi, M.; Little, A.M.; Partheniou, F.; Parham, P. Roles for HLA and KIR polymorphisms in natural killer cell repertoire selection and modulation of effector function. J. Exp. Med. 2006, 203, 633–645. [Google Scholar] [CrossRef] [Green Version]
  34. Rajaram, S.; Canaday, L.M.; Ochayon, D.E.; Rangel, K.M.; Ali, A.; Gyurova, I.E.; Krishnamurthy, D.; Fletcher, J.S.; Reighard, S.D.; Cox, A.; et al. The Promise and Peril of Natural Killer Cell Therapies in Pulmonary Infection. Immunity 2020, 52, 887–889. [Google Scholar] [CrossRef]
  35. Whittaker, E.; Bamford, A.; Kenny, J.; Kaforou, M.; Jones, C.E.; Shah, P.; Ramnarayan, P.; Fraisse, A.; Miller, O.; Davies, P.; et al. Clinical Characteristics of 58 Children with a Pediatric Inflammatory Multisystem Syndrome Temporally Associated With SARS-CoV-2. JAMA 2020, 324, 259–269. [Google Scholar] [CrossRef]
  36. La, D.; Czarnecki, C.; El-Gabalawy, H.; Kumar, A.; Meyers, A.F.A.; Bastien, N.; Simonsen, J.N.; Plummer, F.A.; Luo, M. Enrichment of Variations in KIR3DL1/S1 and KIR2DL2/L3 among H1N1/09 ICU Patients: An Exploratory Study. PLoS ONE 2011, 6, e29200. [Google Scholar] [CrossRef]
  37. Kobasa, D.; Jones, S.M.; Shinya, K.; Kash, J.C.; Copps, J.; Ebihara, H.; Hatta, Y.; Hyun Kim, J.; Halfmann, P.; Hatta, M.; et al. Aberrant innate immune response in lethal infection of macaques with the 1918 influenza virus. Nature 2007, 445, 319–323. [Google Scholar] [CrossRef]
  38. Zwerling, A.; Behr, M.A.; Verma, A.; Brewer, T.F.; Menzies, D.; Pai, M. The BCG World Atlas: A Database of Global BCG Vaccination Policies and Practices. PLoS Med. 2011, 8, e1001012. [Google Scholar] [CrossRef] [Green Version]
  39. Hiby, S.E.; Walker, J.J.; O’Shaughnessy, K.M.; Redman, C.W.G.; Carrington, M.; Trowsdale, J.; Moffett, A. Combinations of Maternal KIR and Fetal HLA-C Genes Influence the Risk of Preeclampsia and Reproductive Success. J. Exp. Med. 2004, 200, 957–965. [Google Scholar] [CrossRef]
  40. Ayoub, F.; Sato, T.; Sakuraba, A. Football and COVID-19 risk: Correlation is not causation. Clin. Microbiol. Infect. 2020. [Google Scholar] [CrossRef]
  41. Ellinghaus, D.; Degenhardt, F.; Bujanda, L.; Buti, M.; Albillos, A.; Invernizzi, P.; Fernández, J.; Prati, D.; Baselli, G.; Asselta, R.; et al. Genomewide Association Study of Severe Covid-19 with Respiratory Failure. N. Engl. J. Med. 2020. [Google Scholar] [CrossRef]
  42. Wang, F.; Huang, S.; Gao, H.; Zhou, Y.; Lai, C.; Li, Z.; Xian, W.; Qian, X.; Li, Z.; Huang, Y.; et al. Initial Whole Genome Sequencing and Analysis of the Host Genetic Contribution to COVID-19 Severity and Susceptibility. medRxiv 2020. [Google Scholar] [CrossRef]
Viruses 12 01333 g001 550
Figure 1. Worldwide allele frequency of HLA-C*05 and risk of deaths due to COVID-19. Linear regression showed a strong correlation between the allele frequency of HLA-C*05 and number of deaths per 1 million population (adjusted R2 0.37, p = 0.0000047).
Figure 1. Worldwide allele frequency of HLA-C*05 and risk of deaths due to COVID-19. Linear regression showed a strong correlation between the allele frequency of HLA-C*05 and number of deaths per 1 million population (adjusted R2 0.37, p = 0.0000047).
Viruses 12 01333 g001
Viruses 12 01333 g002 550
Figure 2. Association between allele frequency of HLA-C*05, proportion of individuals carrying KIR2DS4fl and risk of deaths due to COVID-19. Countries with high, medium and low mortality could be categorized into three distinct groups with a specific pattern of HLA-C*05/KIR2DS4fl expression. KIR: killer cell immunoglobulin-like receptor.
Figure 2. Association between allele frequency of HLA-C*05, proportion of individuals carrying KIR2DS4fl and risk of deaths due to COVID-19. Countries with high, medium and low mortality could be categorized into three distinct groups with a specific pattern of HLA-C*05/KIR2DS4fl expression. KIR: killer cell immunoglobulin-like receptor.
Viruses 12 01333 g002
Viruses 12 01333 g003 550
Figure 3. Pattern of mortality to Spanish flu, swine flu and COVID-19 pandemics. Countries with high, medium and low mortality from each pandemic were categorized into five distinct groups. With the exception of a few outliers, African, American, Asian and European countries were clustered, demonstrating a similar pattern of mortality to each pandemic.
Figure 3. Pattern of mortality to Spanish flu, swine flu and COVID-19 pandemics. Countries with high, medium and low mortality from each pandemic were categorized into five distinct groups. With the exception of a few outliers, African, American, Asian and European countries were clustered, demonstrating a similar pattern of mortality to each pandemic.
Viruses 12 01333 g003
Table
Table 1. Results of linear regression and adjusted p values. HLA: human leukocyte antigen.
Table 1. Results of linear regression and adjusted p values. HLA: human leukocyte antigen.
HLAAdjusted R2p ValueAdjusted p ValueHLAAdjusted R2p ValueAdjusted p ValueHLAAdjusted R2p ValueAdjusted p Value
A*010.160.000590.040 *B*070.240.0000120.00081 *C*01−0.00160.341
A*020.0170.151B*080.160.000700.047 *C*02−0.0130.511
A*030.100.00460.31B*130.0420.0571C*03−0.0240.861
A*110.0200.131B*140.150.00150.099C*04−0.0120.491
A*230.0320.0901B*150.0260.101C*050.370.00000470.00032 *
A*240.0320.0781B*18−0.0110.591C*06−0.0190.631
A*250.0220.151B*27−0.00610.441C*070.180.00230.15
A*26−0.0130.691B*35−0.00940.551C*08−0.000440.331
A*290.0660.0211B*37−0.0140.671C*12−0.0220.841
A*300.00810.221B*38−0.00570.421C*150.0640.0511
A*310.0120.581B*39−0.00790.481C*16−0.0290.901
A*320.0140.181B*400.0300.0961C*170.0210.191
A*330.110.00380.26B*41−0.00270.361C*180.110.091
A*340.0300.121B*420.0570.0541
A*360.0110.231B*440.220.0000320.0022 *
A*43−0.0220.531B*450.0270.121
A*660.0410.0781B*460.0420.121
A*68−0.00940.531B*470.0120.221
A*69−0.00730.421B*480.00490.271
A*740.0510.0481B*49−0.0180.861
A*800.0200.201B*50−0.00620.431
B*510.0120.181
B*520.0720.0201
B*530.0300.101
B*540.0100.261
B*55−0.0160.821
B*560.0200.161
B*570.0330.0761
B*580.0910.00890.60
B*590.0600.191
B*670.0530.141
B*73−0.0290.831
B*810.0210.241
B*820.0300.281
* indicates p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sakuraba, A.; Haider, H.; Sato, T. Population Difference in Allele Frequency of HLA-C*05 and Its Correlation with COVID-19 Mortality. Viruses 2020, 12, 1333. https://doi.org/10.3390/v12111333

AMA Style

Sakuraba A, Haider H, Sato T. Population Difference in Allele Frequency of HLA-C*05 and Its Correlation with COVID-19 Mortality. Viruses. 2020; 12(11):1333. https://doi.org/10.3390/v12111333

Chicago/Turabian Style

Sakuraba, Atsushi, Haider Haider, and Toshiro Sato. 2020. "Population Difference in Allele Frequency of HLA-C*05 and Its Correlation with COVID-19 Mortality" Viruses 12, no. 11: 1333. https://doi.org/10.3390/v12111333

APA Style

Sakuraba, A., Haider, H., & Sato, T. (2020). Population Difference in Allele Frequency of HLA-C*05 and Its Correlation with COVID-19 Mortality. Viruses, 12(11), 1333. https://doi.org/10.3390/v12111333

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

Article Metrics

Back to TopTop