Skip to main content


Front. Immunol., 03 February 2022
Sec. Viral Immunology

HLA-C*04:01 Affects HLA Class I Heterozygosity and Predicted Affinity to SARS-CoV-2 Peptides, and in Combination With Age and Sex of Armenian Patients Contributes to COVID-19 Severity

Anahit Hovhannisyan1,2*, Vergine Madelian3, Sevak Avagyan3, Mihran Nazaretyan3, Armine Hyussyan3, Alina Sirunyan3, Rubina Arakelyan4, Zorayr Manukyan5, Levon Yepiskoposyan2, Karine R. Mayilyan1† and Frieda Jordan3*†
  • 1Institute of Molecular Biology, National Academy of Sciences, Yerevan, Armenia
  • 2Russian-Armenian University, Yerevan, Armenia
  • 3Armenian Bone Marrow Donor Registry Charitable Trust, Yerevan, Armenia
  • 4ClinSoft Armenia, Yerevan, Armenia
  • 5ClinStat Group, Lexington, MA, United States

The novel SARS-CoV-2 coronavirus infection has become a global health concern, causing the COVID-19 pandemic. The disease symptoms and outcomes depend on the host immunity, in which the human leukocyte antigen (HLA) molecules play a distinct role. The HLA alleles have an inter-population variability, and understanding their link to the COVID-19 in an ethnically distinct population may contribute to personalized medicine. The present study aimed at detecting associations between common HLA alleles and COVID-19 susceptibility and severity in Armenians. In 299 COVID-19 patients (75 asymptomatic, 102 mild/moderate, 122 severe), the association between disease severity and classic HLA-I and II loci was examined. We found that the advanced age, male sex of patients, and sex and age interaction significantly contributed to the severity of the disease. We observed that an age-dependent effect of HLA-B*51:01 carriage [odds ratio (OR)=0.48 (0.28-0.80), Pbonf <0.036] is protective against severe COVID-19. Contrary, the HLA-C*04:01 allele, in a dose-dependent manner, was associated with a significant increase in the disease severity [OR (95% CI) =1.73 (1.20-2.49), Pbonf <0.021] and an advancing age (P<0.013). The link between HLA-C*04:01 and age was secondary to a stronger association between HLA-C*04:01 and disease severity. However, HLA-C*04:01 exerted a sex-dependent differential distribution between clinical subgroups [females: P<0.0012; males: P=0.48]. The comparison of HLA-C*04:01 frequency between subgroups and 2,781 Armenian controls revealed a significant incidence of HLA-C*04:01 deficiency in asymptomatic COVID-19. HLA-C*04:01 homozygous genotype in patients blueprinted a decrease in heterozygosity of HLA-B and HLA class-I loci. In HLA-C*04:01 carriers, these changes translated to the SARS-CoV-2 peptide presentation predicted inefficacy by HLA-C and HLA class-I molecules, simultaneously enhancing the appropriate HLA-B potency. In patients with clinical manifestation, due to the high prevalence of HLA-C*04:01, these effects provided a decrease of the HLA class-I heterozygosity and an ability to recognize SARS-CoV-2 peptides. Based on our observations, we developed a prediction model involving demographic variables and HLA-C*04:01 allele for the identification of potential cases with the risk of hospitalization (the area under the curve (AUC) = 86.2%) or severe COVID-19 (AUC =71%).


The novel coronavirus disease COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has emerged at unprecedented scale and ferocity, putting enormous pressure on health systems globally. As of 31 August 2021, according to the World Health Organization, the cumulative number of coronavirus-infected patients identified is 216,867,420 (1). In terms of the clinical course of the disease, its spectrum can range from asymptomatic to severe manifestation, with acute respiratory failure leading to death in some cases (2, 3). The early reports on the COVID-19 clinical characteristics indicate that up to 20% of patients required admission to critical care, such as intensive care unit (ICU), and prolonged support by mechanical ventilation in case of respiratory failure. In-hospital mortality rate was as high as 28% (3, 4). Preceding health conditions, age and sex of individuals, environmental and socioeconomic factors are known to affect both the susceptibility and the severity of COVID-19 (5, 6). Meanwhile, knowledge on the influence of the host genetic profile on the susceptibility and outcome of the disease is being continuously generated.

HLA molecules play a pivotal role in the development and regulation of immune response and are shown to be associated with the predisposition to multiple viral diseases, such as HIV (7, 8), HBV, HCV (9), CMV (10), and SARS-CoV-1 (1114). Several studies have reported a link between HLA alleles and COVID-19 infection in different populations (1520), some of which also indicated the association between an HLA-C allele and the severe clinical course of the disease (17, 19, 21, 22). Nevertheless, our understanding of the HLA genetic contribution to the COVID-19 clinical manifestation and severity remains insufficient and requires further exploration. Furthermore, because of HLA genetic diversity across human populations, the association between the disease severity and specific HLA alleles may vary from one ethnic/ancestral group to another. Therefore, HLA landscape of each population should be also considered in COVID-19 association studies (23).

Here, we report the results of the HLA classic loci association analyses with severity of COVID-19 infection in the Armenians, a population with its unique HLA variations developed over long-term genetic isolation (24). Our findings implicate the HLA-C*04:01 allele in the severe course of the disease in a dose-dependent manner, and the HLA-B*51:01 carriage in the protection against the disease but with a lesser extent.

Materials and Methods


The study included 299 SARS-CoV-2 infected patients. During April-May 2020, blood and nasopharyngeal swab samples from patients were collected in six hospitals of the Health Ministry of the Republic of Armenia: National Center for Infection Diseases, Surb Grigor Lusavorich, Surb Astvatsamayr and Artashat Medical Centres, Maralik and Akhuryan Mother and Child Health Centers. The diagnosis was performed on nasopharyngeal specimens by detection of SARS-CoV-2 RNA, using the following kits: SARS-CoV-2/SARS-CoV Multiplex REAL-TIME PCR Detection Kit (LLC “NPF DNA-Technology”, Cat. No. R3-P436-23/9EU), LiliF™ COVID-19 Real-time RT-PCR Kit (iNtRON Biotechnology, Cat. No. IPH21505.50) and Real-time-fluorescent-RT-PCR (BGI, Cat. No. HW5105). Patients were classified into three groups according to their disease symptoms. The severe disease group (n=122) comprised patients with severe respiratory problems, who required critical care at ICU and the use of oxygen supplementation and/or mechanical ventilation. In the second group, namely mild (in figures) symptomatic group (n=102), were included patients hospitalized during the same period with various mild and moderate COVID-19 symptoms, such as fever, loss of taste, body aches, fatigue, etc., but not requiring ICU admission. The asymptomatic group (n=75) consisted of hospital staff at the COVID-19 treating centers who tested positive for the virus but had no symptoms.

To compare the allelic frequency profiles of patient groups with those specific to the general population, we used HLA data of 2,781 Armenian individuals who were registered donors at the Armenian Bone Marrow Donor Registry (ABMDR) before the pandemic. This group was used as a control cohort for age-genotype and allele frequency analyses. Informed consent for genetic studies was obtained from all subjects. The study protocol was approved by the Ethics Committee of Yerevan State Medical University after Mkhitar Heratsi (Meeting of 04.06.2020, Protocol N 7-3/20), and carried out following the Declaration of Helsinki.

HLA Sequencing

Blood samples (3 ml, EDTA) were sent out to HistoGenetics LLC laboratory (Ossining, NY) for HLA sequencing of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 loci (25). Next-generation sequencing was performed using Illumina MiSeq and PacBio platforms.

The SARS-CoV-2 Peptides Binding Affinity by HLA Class I Molecules

A number of SARS-CoV-2 peptides binding ability for each HLA-A, HLA-B, HLA-C, and total HLA class I molecules was calculated from the data provided by Nguyen et al. (26). Based on the availability of HLA alleles represented in the Armenian population and, particularly, in this dataset, we were able to predict the number of SARS-CoV-2 peptides binding to the HLA-A for 281, HLA-A for 276, and HLA-C and HLA class I molecules for 234 and 203 patients, respectively.


A Hardy-Weinberg equilibrium (HWE) of HLA genotypes was calculated with PyPop (Python for Population genomics) v.0.7.0 software (27). The HLA allele frequency data for world populations were obtained from the Allele Frequency Net Database (28).

Statistical analyses were performed using SPSS program v. 17.0 and STATISTICA v.6.1 StatSoft Inc. (Tulsa, USA). For general statistics, Fisher’s exact, χ2 or logistic regression tests were used for comparison of frequencies. Kolmogorov-Smirnov test with Lilliefors correction was applied for inspection of normality of the continuous data distribution. Spearman’s rank or Pearson’s r tests were used for correlation analyses based on the data distribution. Two-tailed Student’s t and Mann-Whitney tests were used for comparison of two groups, while ANOVA and Kruskal–Wallis tests were used for multiple groups based on the distribution of a continuous variable. Linear regression analysis was applied for continuous variable(s) to test the difference between several groups, and/or for conditional analysis of covariates. In a linear regression of the HLA class I affinity to the SARS-CoV-2 peptides, four outliers were omitted, as they had a value more than 1.5 times the interquartile range from either 25th percentile or 75th percentile of the group.

To test the association of age and sex with the COVID-19 severity characteristics as well as with the HLA allele genotypes, unconditional linear and ordinal logistic regression analyses were performed, respectively. For association analyses of each common HLA allele (freq > 0.05%) with disease severity, unconditional ordinal logistic regressions were applied in additive mode:


where dose is the number of the HLA allele of interest in the genotype of the given HLA locus (or the number of diverse HLA class I alleles in all three loci, when testing the effect of HLA class I heterozygosity). Thus, allelic-doses 0, 1, and 2 correspond to genotype with homozygote absence, heterozygote and homozygote presence of a given HLA allele (i.e., HLA genotypes: 0/0, 0/1, and 1/1, respectively), while for HLA class I heterozygosity number of alleles ranged between 3 (homozygote for all three HLA class I loci) and 6 (heterozygote for HLA-A, -B and –C loci). P values of nominally significant associations were adjusted for multiple comparisons using Bonferroni correction (Pbonf). Since age and sex, as well as sex and age interaction were significantly different between subgroups of cases, subsequent conditional ordinal logistic regression analyses were performed with the following model:


where OR is calculated from the probability of having given or much severe course of disease, dose is the number of the HLA allele of interest in the genotype of the given HLA locus, θasymptomatic is β of the first order intercept and θmild is β of the second order intercept. The models fit were assessed by Akaike`s information Criterion (AIC) (29) and Bayesian information Criterion (BIC). In case of HLA-C*04:01 these variables were 556.9 and 582.4, respectively, while for HLA-B*51:01 carriage – 561.0 and 586.9, respectively. For the contextualization of potential genetic mechanisms of HLA allele association with the COVID-19 severity, the allelic frequencies in asymptomatic, mild/moderate, and severe patient groups were compared with those of the controls from the general population applying the χ2 test.

HLA haplotype prediction and their association with the disease manifestation characteristics were performed using Bringing ImmunoGenomics Data Analysis Workflow Gaps (BIGDAWG) R package (30).

We developed a risk prediction model for detecting severe COVID-19 cases admitted to ICU, as well as for the identification of hospitalized vs. non-hospitalized COVID-19 cases. For this, we merged a risk prediction models of unconditional ordinal logistic regressions of severity characteristics against demographic and genetic predictors. The final association model was the following:


where OR is calculated from the probability of having given or much severe COVID-19 using the following coefficients: βmale vs. female =2.300; βage(male) = 0.049 per year; βage(female) = 0.091 per year; βHLA=0.545 per copy of HLA-C*04:01 allele in genotype.

The receiver operating characteristics (ROC) curve was plotted using SPSS program v.17.0. The area under the curve (AUC) was used to evaluate the ability of the model to discriminate between severe and non-severe cases as well as between asymptomatic and cases with clinical symptoms of COVID-19 disease.


Demographic Characteristics of Patients

The demographic characteristics of the controls and the patients from each group are summarized in Table 1. The huge disparity across cohorts in demographics was based on the diverse aim in sample collection (i.e. relatively young controls were registered donors at the Armenian Bone Marrow Donor Registry), the hidden course of the disease in asymptomatic cases (all asymptomatic cases in the study were hospital workers of non-retiremental age with the majority of female staff members, who routinely underwent SARS-CoV-2 RNA test) and the well-documented age association of COVID-19 severity. Particularly, in unconditional ordinal logistic regression analyses, the severity score of COVID-19 infection was found to be significantly correlated with the advanced age (P<9.5E-18; see Figure S1A) but not to the sex of the patients (P>0.87). Despite the latter non-significant association, sex of individuals showed disparity across clinical subgroups (χ2 = 8.5; P<0.014) (Figure S1B). Apparently, there was a significant age and sex interaction as well. Particularly, when age and sex interaction term was added to the model, all the three variables, i.e., age (P<1E-20), sex, and age and sex interaction were significantly associated with severity score, which unveiled ~10-fold higher relative OR for males versus females: [OR ± 95% confidence interval (CI) = 9.98 (1.66-59.97), P<0.012] (Table S1). OR of age effect per year in females versus males was equal to 1.042 (95% CI: 1.011-1.075, P<0.0075).


Table 1 Demographic characteristics of patient subgroups and controls.

Associations of HLA Alleles With Susceptibility and Severity of COVID-19

In the total group of patients, out of the eight investigated loci, only the HLA-A genotype was subtly departed from the HWE (P<0.048). Using data of 299 patients and 2,781 population control subjects, we first checked for the association of HLA common alleles with susceptibility to COVID-19, and did not find any significant association (Table S2).

Three HLA alleles, HLA-C*04:01 (P<0.0035), B*35:01 (P<0.047), and HLA-B*51:01 (P<0.0075) were nominally associated with the disease severity in an unconditional additive mode (Table 2). After applying Bonferroni’s correction, only the association of the HLA-C*04:01 allele remained significant (Pbonf<0.021), while HLA-B*51:01 showed a trend for an association. None of the two and three allelic common haplotypes showed a significant association with the disease severity characteristics (data not shown), indicating an independent nature of HLA-C*04:01 and HLA-B*51:01 associations.


Table 2 Unconditional ordinal logistic regressions for the association of the genotypes of each common allele (freq > 0.05%) of the HLA class I and class II loci with the disease severity characteristics in the additive model of each common allele.

Patients homozygous for the HLA-C*04:01 allele were identified in the mild/moderate and severe subgroups (Figure 1A), which resulted in a dose-dependent increase of the risk for severe clinical manifestation (Figure 1C). In contrast, only one person from each subgroup was homozygote for HLA-B*51:01 (Figure 1B), causing an inflation of statistical error. Therefore, we applied a dominant model for reinspection of this association. Our results indicated that the presence of HLA-B*51:01 allele conferred protection against the severe clinical manifestation of COVID-19 infection (P<0.0051; Figure 1D), which survived the correction for multiple testing (Pbonf <0.036).


Figure 1 The association of HLA class I alleles with COVID-19 severity. (A, B) Allelic distribution of HLA-C*04:01 (A) and HLA-B*51:01 (B) in COVID-19 patient subgroups with asymptomatic (n=75), mild/moderate (n=102), severe (n=122) manifestation of the disease. (C, D) COVID-19 disease severity risk associated odds ratio (OR) of HLA-C*04:01 (C) and HLA-B*51:01 (D) variations in the ordinal logistic regression analyses. Error bars represent 95% confidence intervals (95% CI) of OR. Data of HLA-C genotypes are the result of the ordinal logistic regression analyses using codominant model for visualization purposes. (E, F) Age distribution in the COVID-19 patients and controls from the general population associated with the allelic-load of HLA-C*04:01 (E) and HLA-B*51:01 carriage (F). In (E, F) P-values were derived from ANOVA tests to analyze differential genotype effect separately in cases and controls. In (C–F) “0” denotes any other allele different from the allele of interest. In (E, F) control groups refers to 2,781 individuals who were registered donors at the Armenian Bone Marrow Donor Registry (ABMDR) before the pandemic.

HLA-C*04:01 exerted a sex-dependent differential distribution between subgroups of COVID-19 cases. Particularly, in females but not in males [HLA-C*04:01 per copy: OR (95% CI) = 1.21(0.72-2.02), P =0.48], carrying this allele, appeared to predispose to a relatively severe outcome in dose-dependent manner [HLA-C*04:01 per copy in females: OR (95% CI) = 2.15 (1.44-4.39), P = 0.0012] (Figure S2). Besides, we observed an association of the allelic-load of HLA-C*04:01 and the presence of HLA-B*51:01 with the age of patients (P<0.013 and P<0.0045, respectively; Figures 1E, F). Particularly, the relationship between the age and HLA-C*04:01 mirrored the association between advancing age and the severity score of COVID-19 (Table 1), while that of HLA-B*51:01 was in the opposite direction. We additionally performed the conditional analyses controlling for the age and sex of the patients. Unlike HLA-B*51:01 (P>0.15), the relationship between HLA-C*04:01 and the disease severity did not meaningfully change, reaching the significance threshold of suggestive association (P<0.063). Notably, mirroring a sex-related differential distribution between clinical subgroups, the HLA-C*04:01 dose-dependent association remained significant in females (P<0.014) but not in male patients (P=0.96; see Table S3). We found that the association of HLA-C*04:01 genotype with the age in the COVID-19 cases was secondary to the association between the genotype and the COVID-19 severity characteristics (HLA-C*04:01 association with age conditioned for the disease severity, P=0.13). Moreover, a similar analysis in the control cohort with 2,781 individuals from the general population showed no correlation between allelic-load of HLA-C*04:01 or HLA*B51:01 and age (Figures 1E, F). Thus, HLA-C*04:01 association with advancing age was specific to the patients and rather evolved from the allele precipitation in subgroups of mild/moderate and severe COVID-19 cases with a relatively higher average age (Table 1 and Figure S1A).

To elucidate the association mechanism of the HLA-C and COVID-19 severity, we compared the frequency distribution of the HLA-C*04:01 allele in the three clinical subgroups with a representative cohort of the Armenian population. In comparison with controls, a non-significant difference in HLA-C*04:01 allele frequency was observed for the COVID-19 patients with mild/moderate (P=0.93) and severe (P=0.15) symptoms composing the group of hospitalized inpatients (Figure S3A). In contrast, asymptomatic patients with positive SARS-CoV-2 had significantly lower HLA-C*04:01 frequency (P<0.012). Besides, we found that the frequency of HLA-B*51:01 in discrete subsets of patients was not significantly different from that of the general population (0.063 ≤ P < 0.73; Figure S3B). Thus, our results are consistent across allelic and genotype analyses and indicate that the lack of HLA-C*04:01 is protective against a manifestation of the disease symptoms in Armenian patients rather than its load confers a severe course of the disease.

Testing HLA Class I-Related Theories on COVID-19 Severity

Searching for the biological consequences of detected HLA-C*04:01 association with clinical manifestation, we tested a possible contribution of this allele to the HLA class I heterozygosity and number of SARS-CoV-2 peptides presenting ability, which reportedly are protective against progression of COVID-19 (31). For this goal, we first tested if those HLA class I variables significantly correlate with COVID-19 severity characteristics in this data set. Even though there was no statistical differences for separate HLA-A, HLA-B and HLA-C loci between subgroups (P=0.25, P=0.48, P=0.23, respectively, for HLA-A, -B, and -C), we observed a significant association between low HLA class I heterozygosity and the disease severity [unconditional likelihood ratio test global P< 0.035; Figure 2A], which overcame the covarying effect of age and sex [global P<0.037, in comparison with heterozygote genotype for 0 or 1 HLA class I locus, heterozygotes for two HLA class I loci: OR (95% CI) = 0.23 (0.05-1.025), P=0.054; heterozygotes for all three HLA class I loci: OR (95% CI) = 0.19 (0.046-0.77), P = 0.020]. While testing the hypothesis on the association between COVID-19 severity and predicted low number of HLA-recognizable SARS-CoV-2 peptides (26, 31), HLA-A, HLA-B and HLA-C loci showed no difference between subgroups (Mann-Whitney test P>0.24; see Figures S4A–C). Thus, we could not confirm the findings of Iturrieta-Zuazo et al. (31) on separate HLA class I molecules. As to the general HLA class I affinity, we found a better theoretical capacity to bind SARS-CoV-2 peptides in asymptomatic patients compared to mild/moderate and severe cases (Figure 2B).


Figure 2 Testing the theories on the association between HLA class I loci and COVID-19 severity. (A) HLA class I heterozygosity was a predictor of the disease severity. Error bars represent 95% confidence intervals (95% CI) of OR. (B) Differential HLA class I affinity to SARS-CoV-2 peptides was observed in asymptomatic patients (n=45) compared with mild/moderate (n=67) to severe (n=87) cases (P-value was driven by two-tailed t-test) for whom affinity data from the recent in silico study were available (26). (C) Effect of HLA-C*04:01 allelic load on HLA-A (n=281), HLA-B (n=276), HLA-C (n=234) binding capacity to SARS-CoV-2 peptides (P-value was driven by Kruskal–Wallis test). (D) HLA-C*04:01 allele carriage was associated with an attenuated average HLA class I binding affinity to SARS-CoV-2 peptides (P-value was driven by a two-tailed t-test of the data from 78 carriers vs 121 non-carriers). In (B–D), box-plots indicate median and interquartile ranges of data distribution without four outliers.

Strikingly, both of those HLA class I variables, i.e., heterozygosity and SARS-CoV-2 peptide presenting ability, showed stronger associations to the HLA-C*04:01 allele than to the clinical manifestation of COVID-19. First, the HLA-C*04:01 homozygosis was simultaneously associated with low HLA class I and HLA-B heterozygosity (P<4.2E-12 and P<0.0039, respectively; Figure S5). Second, the allelic-load of HLA-C*04:01 increased and decreased, respectively, the median number of SARS-CoV-2 peptides affinity to HLA-B and HLA-C molecules (Figure 2C). Lastly, those counter-effects of HLA-C*04:01-load culminated in a relatively low HLA class I-binding ability in the HLA-C*04:01 carriers (Figure 2D).

As in other populations, the LD between HLA-B and HLA-C loci was the highest across HLA classical genes in our sample (R2 = 0.98). To test if any HLA-C allele homozygosity is a sufficient genetic determinant to drive HLA I and HLA-B low heterozygosity observed in severe COVID-19 due to the high LD between HLA-B and HLA-C, we have performed similar analyses with HLA-C*07:01 (next frequent HLA-C allele in our cohort, see Table 2). We found that although HLA-C*07:01 homozygosity significantly affects HLA-B heterozygosity (P<0.043; cf. HLA-C*04:01: P<0.0056), but in overall statistics, this HLA-C*07:01 effect is not significant (P=0.13; cf. HLA-C*04:01: P<0.0039, see Figure S5B). Besides, due to relatively low frequency than HLA-C*04:01, and non-significant association with COVID-19 (P>0.36), it does not preclude a significant association between a scanty HLA class I heterozygosity and COVID-19 (after controlling for HLA-C*07:01, P<0.035). These comparative data indicate that the association between HLA-C*04:01 and COVID-19 is so substantial that it also drives a decrease in HLA I heterozygosity in severe cases.

In conjunction with the results of an in silico analysis that the HLA-C*04:01 allele is one of the least efficient for SARS-CoV-2 peptide presentation (26), our aforementioned findings suggested that the scanty HLA class I heterozygosity and SARS-CoV-2 peptide presenting ability in COVID-19 cases with clinical symptoms may be driven by the stronger association of the HLA-C*04:01 with the disease severity. Therefore, we performed a couple of reciprocal conditional analyses. In an ordinal logistic regression, we found that the HLA-C*04:01 association with COVID-19 severity retained its significance (P<0.033), despite controlling for HLA class I heterozygosity. On the other hand, the significant diversity in the SARS-CoV-2 peptide presenting ability of the HLA class I molecules across severity groups was constrained (P>0.17). Simultaneously, HLA-B (P<0.030), HLA-C (P<3.3E-07) and HLA class I (P<0.0038) binding ability showed a similar pattern of the association to the HLA-C*04:01 allelic–load, as in unconditional analyses. These results implicate HLA-C*04:01 allele in COVID-19 severity with a biological consequence affecting HLA class I qualitative and quantitative properties, and serving as a genetic determinant for the SARS-CoV-2 peptide presentation quantitative trait (QT) beyond its effect on the HLA-C locus.

Predictors of Hospitalization and ICU Admission of Armenian COVID-19

Based on our findings, we further developed a risk prediction model for detecting severe COVID-19 cases admitted to the ICU, as well as for the identification of hospitalized vs. non-hospitalized COVID-19 cases by merging risk prediction models of unconditional ordinal logistic regressions of severity characteristics against demographic (i.e., age, sex, and sex and age interaction) and HLA-C*04:01 predictors (see Material and Methods). The risk scores better separated hospitalized patients with mild to severe symptoms from asymptomatic COVID-19 subjects than severe versus non-severe cases (Figure 3A), as the individual age was the dominant variable (Figure S6). Nevertheless, the contribution of HLA-C*04:01 allele explained 15.4% of the variation in log-transformed OR score (Figure 3B). Thus, the AUC of a risk prediction model for a severe course of the disease that includes HLA-C*04:01, age, sex, and sex and age interaction was 0.710 (95% CI = 0.653-0.768) with a significant difference between severe and non-severe cases (P=6,26E-10; Figure 3C and Figure S7). The resulting OR for predicting severe cases or cases admitted to ICU was relatively moderate [OR=4.22 (2.57-6.92); sensitivity = 0.705, specificity = 0.638]. The same model, however, was more efficient in indicating hospitalization of cases with AUC equal to 0.862 (95% CI =0.815-0.909), with a sensitivity of 0.821 and specificity of 0.720. The threshold of the plot was calculated from the optimal sensitivity and specificity in the AUC, which revealed a significant difference between hospitalized and non-hospitalized cases (P=6,54E-21; Figure 3D). The model was able to discriminate higher OR between hospitalized and non-hospitalized cases [OR (95% CI) = 11.83 (6.43-21.75), Figure 3E].


Figure 3 The model characteristics for prediction of the risk scores between severe (ICU-admitted) and non-severe, as well as hospitalized and non-hospitalized COVID-19 cases. (A) Distribution of risk scores in different subsets of COVID-19 cases based on the clinical picture. (B) Distribution of risk scores in different subsets of COVID-19 cases based on the genomic number of HLA-C*04:01 allele. P-value was driven from Pearson correlation test. (C, D) Receiver operating characteristic (ROC) curve involving the contribution of HLA-C*04:01, sex, and sex-dependent age effects for discrimination of severe vs. non-severe (C) and hospitalized vs. non-hospitalized (D) COVID-19 cases. True positive rate (sensitivity) was plotted against false positive rate (1-specificity). The areas under the curves are presented for each outcome. (E) Histogram plot of hospitalized and non-hospitalized cases versus distribution of ln(OR). The dotted red line represents the threshold of the plot obtained from the AUC (Area Under Curve) with optimal sensitivity (0.821) and specificity (0.720). A significant difference was observed between risk and non-risk cases (OR (95% CI) = 11.83 [6.43-21.75], P=6,54E-21).


Based on the conjoint approach of HLA classical loci genotype and allelic distribution analyses between Armenian population controls and COVID-19 cases, as well as between patients with different clinical phenotypes, we discovered an age-related protective effect of the HLA-B*51:01 carriage against COVID-19 severity. In contrast, we found that HLA-C*04:01 allelic-load contributes to the risk for the severe COVID-19 independently from age and classical HLA variables. However, the HLA-C*04:01 association with COVID-19 manifestation is rather in concert with female gender and older age. Along with the expected significant decrease in cumulative heterozygosity of HLA class I loci, we report a previously undescribed decrease in heterozygosity of HLA-B locus blueprinted by the HLA-C*04:01 homozygous genotype. These changes in HLA heterogeneity may translate to the SARS-CoV-2 peptide presentation relative inefficacy by HLA-C and HLA class I molecules in HLA-C*04:01 carriers, simultaneously enhancing the appropriate potency of HLA-B. In patients with mild/moderate to severe COVID-19, due to high prevalence of HLA-C*04:01, these effects provide a decrease of the HLA class I loci heterozygosity and a down-modulation of the predicted HLA class I ability to recognize SARS-CoV-2 peptides. Given these findings, we have developed a risk prediction model for the early identification of Armenian COVID-19 cases with a potential to be hospitalized or, subsequently, to be admitted to intensive care units. The genomic number of HLA-C*04:01 allele and demographic variables (i.e. age and sex) compose the model, in which >15% potential share in the detection of cases with adverse clinical phenotypes belongs to HLA-C*04:01.

Our finding on the association between HLA-C*04:01 genomic quantity and COVID-19 severity is consistent with some previous reports of cross-sectional studies, involving Indian (32), German (21) and USA (multiethnic) (21, 22) cohorts. Consistently to our results, Warren and Birol (22) also report that the high prevalence of HLA-C*04:01 allele may predispose to severe COVID-19 in patients with a white ethnic background, female gender and advanced age (≥65 years old). This phenomenon of interaction with sex is not unusual for the HLA and COVID-19 association, as it can appear even on haplotype level as well (33). Nevertheless, most studies land partial or no support to our findings (16, 17, 19, 34, 35). Thus, the first genome-wide association study (GWAS) reports a trend for HLA-C*04:01 association to the disease severity in the Italian population with a marginal significance, and a significantly lower number of HLA class I bound SARS‐CoV‐2 peptides in the most severe group of the Italian and Spanish patients (35). Another study indicates HLA-C*04:01 association with susceptibility but not COVID-19 severity in the Sardinian/Italian population (34). Despite this allelic variability conferring COVID‐19 or its clinical heterogeneity across ethnicities, the inability of the identified alleles to effectively bind the SARS‐CoV‐2 viral peptides and to trigger an adequate immune response among others may be one of the main causative mechanisms for severe disease (20, 26, 31, 32, 36).

HLA class I molecules present antigens to cytotoxic T lymphocytes (CTLs) and natural killer (NK) cells. HLA-C, being the least potent in antigen presentation to CTLs (37), works extremely well in controlling and modulating NK cell immune response (38, 39). The HLA-C mediated protection, harnessing NK cells, is relevant to a number of viral infections (4042), including SARS‐CoV‐1 infection (43). In spite of genetic association data, a similar role of HLA-C in severe COVID-19 is yet unknown. An impediment for a direct approximation of genetic findings are transcriptional and epigenetic changes seen in severe cases, such as two differentially methylated CpG sites regulating HLA-C locus (44), and HLA-C alternative transcription and RNA splicing events (45). Another obstacle is a high incidence of lymphopenia (46), in particular, the reduced number of NK cells in the peripheral blood of COVID-19 patients, which is related to the severe stage of infection (4750). Moreover, it is already established that the pathogenesis of severe SARS‐CoV‐2 infection involves a hyper-function of CD4+ and CD8+ T cells and NK cell hyporesponsiveness (49).

Given the key role of HLA-C molecules in antigen presentation to KIR receptors on NK cells, HLA-C*04:01 molecules with (i) a poor affinity to the SARS-CoV-2 peptides (26), (ii) a low cell surface expression ability [due to bearing HLA-C 3UTR’ target site that miR-148a binds (51)], (iii) and differential avidity to its cognate receptors (i.e. high for the inhibiting KIR (KIR2DL1) than the activating KIR [functional KIR2DS4 (52)], may cause NK cell hyporesponsiveness since the early phase of infection (49). As NK cells are at the forefront of antiviral response, particularly at that of the respiratory system (53), NK hyporesponsiveness may delay SARS‐CoV‐2 viral recognition and, subsequently, generate an inadequate immune response, supporting the virus propagation and lung damage, such as seen in severe pneumonia with a genetically predicted HLA-C low expression (54). A series of genetic studies with a Sardinian/Italian sample potentially support this hypothesis, by reporting a significantly high prevalence of both, HLA-C*04:01 allele and NK-cell inhibitor KIR2DL1 receptor, in COVID-19 cases (34, 55).

As the avidity of KIR2DS4 receptor is limited to several HLA class I alleles (52), the other mechanism may involve HLA-C*04:01 antigen presentation to this receptor. At the late phase of the COVID-19 progression, in the frame of cytokine storm and lymphopenia (46, 56), the lung residing infected monocytes might recruit a new portion of NK cells from blood (48). The specific induction of HLA-C expression on differentiated macrophages in the lung of COVID-19 HLA-C*04:01 carriers (57) may trigger CTL and NK cells activation via KIR2DS4 (52), and result in those cells exhaustion. Of note, the two parts of our hypothesis are not mutually exclusive and rather discriminate differential spatial and temporal responses of the HLA-C restricted CTL and NK cells. Remarkably, considerable literature is supportive of the latter interpretation, as several reports from different populations with high prevalence of the functional KIR2DS4*001 gene (58, 59) delineate a risk association of HLA-A*11:01, HLA-C*05:01, HLA-C*04:01, and HLA-C*16:01 to the severe or fatal COVID-19 (16, 17, 19, 21, 22, 34, 6064), and all of those alleles are recognized by the functional KIR2DS4. From the other side of the same coin, HLA-B enhanced affinity to the SARS-CoV-2 viral peptides, coinciding with the carriage of HLA-C*04:01, may over-activate CD8+ T cells in severe patients, since a significantly greater number of CTLs responses are HLA-B restricted (65).

Thus, the results of our study strongly suggest a putative role of HLA-C genetic variation in the development of a specific immune response to COVID-19. The exact role of HLA-C*04:01 to the COVID-19 severity merits further study. However, screening for HLA-C genotype in patients from the populations, in which HLA-C*04:01 is a risk allele for the severe COVID-19, in addition to demographic variables, may support strategic clinical management of the disease.

Data Availability Statement

The data presented in the study are deposited in the Figshare repository, with the following link

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of Yerevan State Medical University after Mkhitar Heratsi (Meeting of 04.06.2020, Protocol N 7-3/20). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author Contributions

FJ initiated the study. VM designed and developed the project. AHo, KM, and LY performed the data analyses and interpretation of the results. KM provided data analysis methodology and the results conceptualization. ZM and RA performed ABMDR database data analyses. AHy, SA, MN, and AS collected the data. AHo and KM performed statistics, designed the figures and wrote the manuscript. All authors contributed to the final version of the manuscript.


The research was supported in part by funds from Dr. George and Carolann Najarian, Armenian Medical Society of California. AnH and LY were supported by the 21AG-1F025 Grant from State Committee Science of the Ministry of Education and Science of the Republic of Armenia. AnH, LY, and KRM acknowledge the grant E17 provided by the State Committee Science of the Ministry of Education and Science of the Republic of Armenia.

Conflict of Interest

RA and ZM declare that they work at the ClinSoft Armenia, Yerevan, Armenia, and ClinStat Group, Lexington, MA.

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

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.


We are grateful to the Ministry of Health of the Republic of Armenia, Yerevan State Medical University, and HistoGenetics LLC.

Supplementary Material

The Supplementary Material for this article can be found online at:


1. World Health Organization. (2021). Available at: (Accessed August 31, 2021).

Google Scholar

2. Chen N, Zhou M, Dong X, Qu J, Gong F, Han Y, et al. Epidemiological and Clinical Characteristics of 99 Cases of 2019 Novel Coronavirus Pneumonia in Wuhan, China: A Descriptive Study. Lancet (2020) 395(10223):507–13. doi: 10.1016/S0140-6736(20)30211-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Weiss P, Murdoch DR. Clinical Course and Mortality Risk of Severe COVID-19. Lancet (2020) 395:1014–5. doi: 10.1016/S0140-6736(20)30633-4

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical Course and Risk Factors for Mortality of Adult Inpatients With COVID-19 in Wuhan, China: A Retrospective Cohort Study. Lancet (2020) 395(10229):1054–62. doi: 10.1016/S0140-6736(20)30566-3

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Liu X, Zhou H, Zhou Y, Wu X, Zhao Y, Lu Y, et al. Risk Factors Associated With Disease Severity and Length of Hospital Stay in COVID-19 Patients. J Infect (2020) 81(1):e95–7. doi: 10.1016/j.jinf.2020.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Zheng Z, Peng F, Xu B, Zhao J, Liu H, Peng J, et al. Risk Factors of Critical & Mortal COVID-19 Cases: A Systematic Literature Review and Meta-Analysis. J Infect (2020) 81(2):e16–25. doi: 10.1016/j.jinf.2020.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Carrington M, Nelson GW, Martin MP, Kissner T, Vlahov D, Goedert JJ, et al. HLA and HIV-1: Heterozygote Advantage and B* 35-Cw* 04 Disadvantage. Science (1999) 283(5408):1748–52. doi: 10.1126/science.283.5408.1748

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Matthews PC, Koyanagi M, Kløverpris HN, Harndahl M, Stryhn A, Akahoshi T, et al. Differential Clade-Specific HLA-B* 3501 Association With HIV-1 Disease Outcome Is Linked to Immunogenicity of a Single Gag Epitope. J Virol (2012) 86(23):12643–54. doi: 10.1128/JVI.01381-12

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Singh R, Kaul R, Kaul A, Khan K. A Comparative Review of HLA Associations With Hepatitis B and C Viral Infections Across Global Populations. WJG (2007) 13(12):1770–87. doi: 10.3748/wjg.v13.i12.1770

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Hassan J, O’Neill D, Honari B, De Gascun C, Connell J, Keogan M, et al. Cytomegalovirus Infection in Ireland: Seroprevalence, HLA Class I Alleles, and Implications. Medicine (2016) 95(6):e2735. doi: 10.1097/MD.0000000000002735

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lin M, Tseng HK, Trejaut JA, Lee HL, Loo JH, Chu CC, et al. Association of HLA Class I With Severe Acute Respiratory Syndrome Coronavirus Infection. BMC Med Genet (2003) 4(1):1–7. doi: 10.1186/1471-2350-4-9

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Ng MH, Lau KM, Li L, Cheng SH, Chan WY, Hui PK, et al. Association of Human-Leukocyte-Antigen Class I (B* 0703) and Class II (DRB1* 0301) Genotypes With Susceptibility and Resistance to the Development of Severe Acute Respiratory Syndrome. J Infect Dis (2004) 190(3):515–8. doi: 10.1086/421523

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Keicho N, Itoyama S, Kashiwase K, Phi NC, Long HT, Van Ban V, et al. Association of Human Leukocyte Antigen Class II Alleles With Severe Acute Respiratory Syndrome in the Vietnamese Population. Hum Immunol (2009) 70(7):527–31. doi: 10.1016/j.humimm.2009.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Sanchez-Mazas A. HLA Studies in the Context of Coronavirus Outbreaks. Swiss Med Wkly (2020) 150:w20248. doi: 10.4414/smw.2020.20248

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Novelli A, Andreani M, Biancolella M, Liberatoscioli L, Passarelli C, Colona VL, et al. HLA Allele Frequencies and Susceptibility to COVID-19 in a Group of 99 Italian Patients. Hla (2020) 96(5):610–4. doi: 10.1111/tan.14047

PubMed Abstract | CrossRef Full Text | Google Scholar

16. 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) 96(2):194–6. doi: 10.1111/tan.13941

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang F, Huang S, Gao R, Zhou Y, Lai C, Li Z, et al. Initial Whole-Genome Sequencing and Analysis of the Host Genetic Contribution to COVID-19 Severity and Susceptibility. Cell Discov (2020) 6(1):1–16. doi: 10.1038/s41421-020-00231-4

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Yung YL, Cheng CK, Chan HY, Xia JT, Lau KM, Wong RS, et al. Association of HLA-B22 Serotype With SARS-CoV-2 Susceptibility in Hong Kong Chinese Patients. Hla (2021) 97(2):127–32. doi: 10.1111/tan.14135

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Khor SS, Omae Y, Nishida N, Sugiyama M, Kinoshita N, Suzuki T, et al. HLA-A*11:01:01:01, HLA-C*12:02:02:01-HLA-B*52:01:02:02, Age and Sex Are Associated With Severity of Japanese COVID-19 With Respiratory Failure. Front Immunol (2021) 12:658570. doi: 10.3389/fimmu.2021.658570

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Shkurnikov M, Nersisyan S, Jankevic T, Galatenko A, Gordeev I, Vechorko V, et al. Association of HLA Class I Genotypes With Severity of Coronavirus Disease-19. Front Immunol (2021) 12:641900. doi: 10.3389/fimmu.2021.641900

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Weiner 3J, Suwalski P, Holtgrewe M, Rakitko A, Thibeault C, Müller M, et al. Increased Risk of Severe Clinical Course of COVID-19 in Carriers of HLA-C* 04:01. EClinicalMedicine (2021) 40:101099. doi: 10.1016/j.eclinm.2021.101099

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Warren RL, Birol I. HLA Alleles Measured From COVID-19 Patient Transcriptomes Reveal Associations With Disease Prognosis in a New York Cohort. PeerJ (2021) 9:e12368. doi: 10.7717/peerj.12368

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Fricke-Galindo I, Falfán-Valencia R. Genetics Insight for COVID-19 Susceptibility and Severity: A Review. Front Immunol (2021) 12:622176. doi: 10.3389/fimmu.2021.622176

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Matevosyan L, Chattopadhyay S, Madelian V, Avagyan S, Nazaretyan M, Hyussian A, et al. HLA-A, HLA-B, and HLA-DRB1 Allele Distribution in a Large Armenian Population Sample. Tissue Antigens (2011) 78(1):21–30. doi: 10.1111/j.1399-0039.2011.01668.x

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Cereb N, Kim HR, Ryu J, Yang SY. Advances in DNA Sequencing Technologies for High Resolution HLA Typing. Hum Immunol (2015) 76(12):923–7. doi: 10.1016/j.humimm.2015.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Nguyen A, David JK, Maden SK, Wood MA, Weeder BR, Nellore A, et al. Human Leukocyte Antigen Susceptibility Map for Severe Acute Respiratory Syndrome Coronavirus 2. J Virol (2020) 94(13):e00510-20. doi: 10.1128/JVI.00510-20

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lancaster AK, Single RM, Solberg OD, Nelson MP, Thomson G. PyPop Update–a Software Pipeline for Large-Scale Multilocus Population Genomics. Tissue Antigens (2007) 69(1):192–7. doi: 10.1111/j.1399-0039.2006.00769.x

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Gonzalez-Galarza FF, Christmas S, Middleton D, Jones AR. Allele Frequency Net: A Database and Online Repository for Immune Gene Frequencies in Worldwide Populations. Nucleic Acids Res (2010) 39(suppl_1):D913–919. doi: 10.1093/nar/gkq1128

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Akaike H. A New Look at the Statistical Model Identification. IEEE Trans (1974) 19(6):716–23. doi: 10.1109/TAC.1974.1100705

CrossRef Full Text | Google Scholar

30. Pappas DJ, Marin W, Hollenbach JA, Mack SJ. Bridging Immunogenomic Data Analysis Workflow Gaps (BIGDAWG): An Integrated Case-Control Analysis Pipeline. Hum Immunol (2016) 77(3):283–7. doi: 10.1016/j.humimm.2015.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Iturrieta-Zuazo I, Rita CG, García-Soidán A, de Malet Pintos-Fonseca A, Alonso-Alarcón N, Pariente-Rodríguez R, et al. Possible Role of HLA Class-I Genotype in SARS-CoV-2 Infection and Progression: A Pilot Study in a Cohort of Covid-19 Spanish Patients. Clin Immunol (2020) 219:108572. doi: 10.1016/j.clim.2020.108572

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Vishnubhotla R, Sasikala M, Ketavarapu V, Reddy DN. High-Resolution HLA Genotyping Identifies Alleles Associated With Severe COVID-19: A Preliminary Study From India. Immun Inflamm Dis (2021) 9(4):1781–85. doi: 10.1002/iid3.481

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Langton D, Bourke S, Lie BA, Reiff G, Natu S, Darlay R, et al. The Influence of HLA Genotype on the Severity of COVID-19 Infection. HLA (2021) 98:14–22. doi: 10.1111/tan.14284

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Littera R, Campagna M, Deidda S, Angioni G, Cipri S, Melis M, et al. Human Leukocyte Antigen Complex and Other Immunogenetic and Clinical Factors Influence Susceptibility or Protection to SARS-CoV-2 Infection and Severity of the Disease Course The Sardinian Experience. Front Immunol (2020) 11:605688. doi: 10.3389/fimmu.2020.605688

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Ellinghaus D, Degenhardt F, Bujanda L, Buti M. Genomewide Association Study of Severe Covid-19 With Respiratory Failure. NEJM (2021) 383(16):1522–34. doi: 10.1056/NEJMoa2020283

CrossRef Full Text | Google Scholar

36. Tavasolian F, Rashidi M, Hatam GR, Jeddi M, Hosseini AZ, Mosawi SH, et al. HLA, Immune Response, and Susceptibility to COVID-19. Front Immunol (2021) 11:601886. doi: 10.3389/fimmu.2020.601886

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Falk CS, Schendel DJ. HLA-C Revisited. Immunol Res (1997) 16(2):203–14. doi: 10.1007/BF02786363

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Colonna M, Borsellino G, Falco M, Ferrara GB, Strominger JL. HLA-C Is the Inhibitory Ligand That Determines Dominant Resistance to Lysis by NK1-And NK2-Specific Natural Killer Cells. PNAS (1993) 90(24):12000–4. doi: 10.1073/pnas.90.24.12000

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Blais ME, Dong T, Rowland-Jones S. HLA-C as a Mediator of Natural Killer and T-Cell Activation: Spectator or Key Player? Immunol (2011) 133(1):1–7. doi: 10.1111/j.1365-2567.2011.03422.x

CrossRef Full Text | Google Scholar

40. Khakoo SI, Thio CL, Martin MP, Brooks CR, Gao X, Astemborski J, et al. HLA and NK Cell Inhibitory Receptor Genes in Resolving Hepatitis C Virus Infection. Science (2004) 305(5685):872–4. doi: 10.1126/science.1097670

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Rauch A, Laird R, McKinnon E, Telenti A, Furrer H, Weber R, et al. Influence of Inhibitory Killer Immunoglobulin-Like Receptors and Their HLA-C Ligands on Resolving Hepatitis C Virus Infection. Tissue Antigens (2007) 69:237–40. doi: 10.1111/j.1399-0039.2006.773_4.x

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Alter G, Martin MP, Teigen N, Carr WH, Suscovich TJ, Schneidewind A, et al. Differential Natural Killer Cell–Mediated Inhibition of HIV-1 Replication Based on Distinct KIR/HLA Subtypes. J Exp Med (2007) 204(12):3027–36. doi: 10.1084/jem.20070695

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Wang SF, Chen KH, Chen M, Li WY, Chen YJ, Tsao CH, et al. Human-Leukocyte Antigen Class I Cw 1502 and Class II DR 0301 Genotypes Are Associated With Resistance to Severe Acute Respiratory Syndrome (SARS) Infection. Viral Immunol (2011) 24(5):421–6. doi: 10.1089/vim.2011.0024

PubMed Abstract | CrossRef Full Text | Google Scholar

44. de Moura MC, Davalos V, Planas-Serra L, Alvarez-Errico D, Arribas C, Ruiz M, et al. Epigenome-Wide Association Study of COVID-19 Severity With Respiratory Failure. EBioMedicine (2021) 66:103339. doi: 10.1016/j.ebiom.2021.103339

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Monaghan SF, Fredericks AM, Jentzsch MS, Cioffi WG, Cohen M, Fairbrother WG, et al. Deep RNA Sequencing of Intensive Care Unit Patients With COVID-19. medRxiv (2021). doi: 10.1101/2021.01.11.21249276

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Tavakolpour S, Rakhshandehroo T, Wei EX, Rashidian M. Lymphopenia During the COVID-19 Infection: What It Shows and What Can be Learned. Immunol Lett (2020) 225:31. doi: 10.1016/j.imlet.2020.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Masselli E, Vaccarezza M, Carubbi C, Pozzi G, Presta V, Mirandola P, et al. NK Cells: A Double Edge Sword Against SARS-CoV-2. Adv Biol Regul (2020) 77:100737. doi: 10.1016/j.jbior.2020.100737

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Vabret N, Britton GJ, Gruber C, Hegde S, Kim J, Kuksin M, et al. Immunology of COVID-19: Current State of the Science. Immunity (2020) 52(6):910–41. doi: 10.1016/j.immuni.2020.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Wang F, Hou H, Luo Y, Tang G, Wu S, Huang M, et al. The Laboratory Tests and Host Immunity of COVID-19 Patients With Different Severity of Illness. JCI Insight (2020) 5(10):e137799. doi: 10.1172/jci.insight.137799

CrossRef Full Text | Google Scholar

50. Wang F, Nie J, Wang H, Zhao Q, Xiong Y, Deng L, et al. Characteristics of Peripheral Lymphocyte Subset Alteration in COVID-19 Pneumonia. J Infect Dis (2020) 221(11):1762–9. doi: 10.1093/infdis/jiaa150

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Kulkarni S, Savan R, Qi Y, Gao X, Yuki Y, Bass SE, et al. Differential microRNA Regulation of HLA-C Expression and Its Association With HIV Control. Nature (2011) 472(7344):495–8. doi: 10.1038/nature09914

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Graef T, Moesta AK, Norman PJ, Abi-Rached L, Vago L, Older Aguilar AM, 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(11):2557–72. doi: 10.1084/jem.20091010

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Tang Y, Li X, Wang M, Zou Q, Zhao S, Sun B, et al. Increased Numbers of NK Cells, NKT-Like Cells, and NK Inhibitory Receptors in Peripheral Blood of Patients With Chronic Obstructive Pulmonary Disease. Clin Dev Immunol (2013) 2013:721782. doi: 10.1155/2013/721782

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Khadzhieva MB, Kuzovlev AN, Salnikova LE. Pneumonia: Host Susceptibility and Shared Genetics With Pulmonary Function and Other Traits. Clin Exp Immunol (2019) 198(3):367–80. doi: 10.1111/cei.13367

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Littera R, Chessa L, Deidda S, Angioni G, Campagna M, Lai S, et al. Natural Killer-Cell Immunoglobulin-Like Receptors Trigger Differences in Immune Response to SARS-CoV-2 Infection. PloS One (2021) 16(8):e0255608. doi: 10.1371/journal.pone.0255608

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Soleimanian S, Yaghobi R. Harnessing Memory NK Cell to Protect Against COVID-19. Front Pharmacol (2020) 11:1309. doi: 10.3389/fphar.2020.01309

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Schaefer MR, Williams M, Kulpa DA, Blakely PK, Yaffee AQ, Collins KL. A Novel Trafficking Signal Within the HLA-C Cytoplasmic Tail Allows Regulated Expression Upon Differentiation of Macrophages. J Immunol (2008) 180(12):7804–17. doi: 10.4049/jimmunol.180.12.7804

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Middleton D, Meenagh A, Gourraud PA. KIR Haplotype Content at the Allele Level in 77 Northern Irish Families. Immunogenetics (2007) 59(2):145–58. doi: 10.1007/s00251-006-0181-7

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Sim MJ, Rajagopalan S, Altmann DM, Boyton RJ, Sun PD, Long EO. Human NK Cell Receptor KIR2DS4 Detects a Conserved Bacterial Epitope Presented by HLA-C. PNAS (2019) 116(26):12964–73. doi: 10.1073/pnas.1903781116

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Gutiérrez-Bautista JF, Rodriguez-Nicolas A, Rosales-Castillo A, López-Ruz MA, Martín-Casares AM, Fernández-Rubiales A, et al. Study of HLA-A, -B, -C, -DRB1 and -DQB1 Polymorphisms in COVID-19 Patients. J Microbiol Immunol Infect (2021). doi: 10.1016/j.jmii.2021.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

61. 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. doi: 10.3390/v12111333

CrossRef Full Text | Google Scholar

62. Lorente L, Martín MM, Franco A, Barrios Y, Cáceres JJ, Solé-Violán J, et al. HLA Genetic Polymorphisms and Prognosis of Patients With COVID-19. Med Intensiva (2021) 45(2):96–103. doi: 10.1016/j.medin.2020.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Toyoshima Y, Nemoto K, Matsumoto S, Nakamura Y, Kiyotani K. SARS-CoV-2 Genomic Variations Associated With Mortality Rate of COVID-19. J Hum Genet (2020) 65(12):1075–82. doi: 10.1038/s10038-020-0808-9

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Tomita Y, Ikeda T, Sato R, Sakagami T. Association Between HLA Gene Polymorphisms and Mortality of COVID-19: An In Silico Analysis. Immun Inflamm Dis (2020) 8(4):684–94. doi: 10.1002/iid3.358

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Kiepiela P, Leslie AJ, Honeyborne I, Ramduth D, Thobakgale C, Chetty S, et al. Dominant Influence of HLA-B in Mediating the Potential Co-Evolution of HIV and HLA. Nature (2004) 432(7018):769–75. doi: 10.1038/nature03113

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: COVID-19 severity, HLA classical genes, Armenian population, HLA-I heterozygosity, affinity to SARS-CoV-2, HLA-C*04:01, HLA-B*51:01, severity risk modelling

Citation: Hovhannisyan A, Madelian V, Avagyan S, Nazaretyan M, Hyussyan A, Sirunyan A, Arakelyan R, Manukyan Z, Yepiskoposyan L, Mayilyan KR and Jordan F (2022) HLA-C*04:01 Affects HLA Class I Heterozygosity and Predicted Affinity to SARS-CoV-2 Peptides, and in Combination With Age and Sex of Armenian Patients Contributes to COVID-19 Severity. Front. Immunol. 13:769900. doi: 10.3389/fimmu.2022.769900

Received: 02 September 2021; Accepted: 13 January 2022;
Published: 03 February 2022.

Edited by:

Vijayakumar Velu, Emory University, United States

Reviewed by:

Anusmita Sahoo, Emory University, United States
Nanda Kishore Routhu, Emory University, United States

Copyright © 2022 Hovhannisyan, Madelian, Avagyan, Nazaretyan, Hyussyan, Sirunyan, Arakelyan, Manukyan, Yepiskoposyan, Mayilyan and Jordan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Anahit Hovhannisyan,; Frieda Jordan,

†These authors shared senior authorship