Impact Factor 7.561 | CiteScore 8.1
More on impact ›

ORIGINAL RESEARCH article

Front. Immunol., 04 December 2020 | https://doi.org/10.3389/fimmu.2020.605688

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

Roberto Littera1,2*†, Marcello Campagna3*†, Silvia Deidda4, Goffredo Angioni5, Selene Cipri1,6, Maurizio Melis3, Davide Firinu3, Simonetta Santus7, Alberto Lai7, Rita Porcella1, Sara Lai1, Stefania Rassu1, Rosetta Scioscia3, Federico Meloni3, Daniele Schirru3, William Cordeddu5, Marta Anna Kowalik8, Maria Serra9, Paola Ragatzu9, Mauro Giovanni Carta3, Stefano Del Giacco3, Angelo Restivo10, Simona Deidda10, Sandro Orrù9, Antonella Palimodde4, Roberto Perra4, Germano Orrù11, Maria Conti12, Cinzia Balestrieri12, Giancarlo Serra12, Simona Onali8, Francesco Marongiu3, Andrea Perra2,8*† and Luchino Chessa2,3,12*†
  • 1Complex Structure of Medical Genetics, R. Binaghi Hospital, Area Socio-Sanitaria Locale (ASSL) Cagliari, Azienda per la Tutela della Salute (ATS) Sardegna, Italy
  • 2Associazione per l’Avanzamento della Ricerca per i Trapianti O.d.V., non profit organisation, Cagliari, Italy
  • 3Department of Medical Sciences and Public Health, University of Cagliari, Cagliari, Italy
  • 4Complex Structure of Pneumology, PO SS Trinità, ASSL Cagliari, ATS Sardegna, Cagliari, Italy
  • 5Complex Structure of Infectious Diseases, PO SS Trinità, ASSL Cagliari ATS Sardegna, Cagliari, Italy
  • 6Department of Medical Biotechnology and Translational Medicine, University of Milan, Milan, Italy
  • 7Unitá di Crisi Locale (UCL) ATS Sardegna, Cagliari, Italy
  • 8Unit of Oncology and Molecular Pathology, Department of Biomedical Sciences, University of Cagliari, Cagliari, Italy
  • 9Medical Genetics, Department of Medical Sciences and Public Health, University of Cagliari, Cagliari, Italy
  • 10Colorectal Surgery Unit, Department of Surgical Science, University of Cagliari, Cagliari, Italy
  • 11Molecular Biology Service Laboratory, Department of Surgical Science, University of Cagliari, Cagliari, Italy
  • 12Liver Unit, Department of Internal Medicine, University Hospital of Cagliari, Cagliari, Italy

Aim: SARS-CoV-2 infection is a world-wide public health problem. Several aspects of its pathogenesis and the related clinical consequences still need elucidation. In Italy, Sardinia has had very low numbers of infections. Taking advantage of the low genetic polymorphism in the Sardinian population, we analyzed clinical, genetic and immunogenetic factors, with particular attention to HLA class I and II molecules, to evaluate their influence on susceptibility to SARS-CoV-2 infection and the clinical outcome.

Method and Materials: We recruited 619 healthy Sardinian controls and 182 SARS-CoV-2 patients. Thirty-nine patients required hospital care and 143 were without symptoms, pauci-symptomatic or with mild disease. For all participants, we collected demographic and clinical data and analyzed the HLA allele and haplotype frequencies.

Results: Male sex and older age were more frequent in hospitalized patients, none of whom had been vaccinated during the previous seasonal flu vaccination campaignes. Compared to the group of asymptomatic or pauci-symptomatic patients, hospitalized patients also had a higher frequency of autoimmune diseases and glucose-6-phosphate-dehydrogenase (G6PDH) deficiency. None of these patients carried the beta-thalassemia trait, a relatively common finding in the Sardinian population. The extended haplotype HLA-A*02:05, B*58:01, C*07:01, DRB1*03:01 [OR 0.1 (95% CI 0–0.6), Pc = 0.015] was absent in all 182 patients, while the HLA-C*04:01 allele and the three-loci haplotype HLA-A*30:02, B*14:02, C*08:02 [OR 3.8 (95% CI 1.8–8.1), Pc = 0.025] were more frequently represented in patients than controls. In a comparison between in-patients and home care patients, the HLA-DRB1*08:01 allele was exclusively present in the hospitalized patients [OR > 2.5 (95% CI 2.7–220.6), Pc = 0.024].

Conclusion: The data emerging from our study suggest that the extended haplotype HLA-A*02:05, B*58:01, C*07:01, DRB1*03:01 has a protective effect against SARS-CoV-2 infection in the Sardinian population. Genetic factors that resulted to have a negative influence on the disease course were presence of the HLA-DRB1*08:01 allele and G6PDH deficiency, but not the beta-thalassemic trait. Absence of influenza vaccination could be a predisposing factor for more severe disease.

Introduction

Outbreak of the coronavirus disease (COVID-19) was declared a pandemic by the World Health Organization (WHO) on 11 March 2020 (1) and since then continues to pose a threat to the sustainability of public health systems worldwide. Although consensus on total mortality and morbidity has yet to be reached, current data show that over 50 million people have contracted the SARS-CoV-2 virus globally, and of these, more than 1,250,000 have died (2).

Although only a few months have passed since the onset of the pandemic, a myriad of attempts have been made to unravel the unknowns surrounding disease pathogenesis and the related clinical consequences. Currently, it is assumed that the mean incubation period is 4–6 days with an estimated range of 2–14 days, although rare cases of up to 27 days have been reported. The large majority of patients have a good prognosis with an asymptomatic or pauci-symptomatic disease course (3). However, about 20% of infected subjects become seriously ill and 4% develop severe and even fatal disease. Clinical features in hospitalized patients include complications such as pneumonia and, in the most severe cases, acute respiratory distress syndrome (ARDS). Devastating neurological and gastrointestinal symptoms have also been reported. Sepsis is another feared complication of COVID-19 that can cause lasting damage to the lungs and other organs. Moreover, heart failure, kidney failure, liver damage, hypercoagulability, septic shock, and multiple-organ failure (MOF) have been shown to precipitate death. The most severe cases of COVID-19 with admission to the Intensive Care Unit (ICU) are generally more frequent in males and the elderly, especially those with comorbidities such as chronic cardiovascular and/or respiratory disease (48).

A plethora of drugs and treatments have been tried around the world in the battle against the pandemic but, so far, none of those tested or included in clinical trials have yielded a significant reduction in the morbidity and mortality rates (9). Contemporarily, the scientific community is striving hard to produce a vaccine but in spite of the many encouraging results, health experts warn that a coronavirus vaccine is unlikely to be available in time to stop the second wave of the epidemic which is expected to peak in European countries, including Italy, this autumn or winter (10).

After China and South Korea, Italy was the first European country to be hit by the explosive pandemic potential of COVID-19, which sparked from the first cases diagnosed in the province of Lodi (11) to a quickly rising toll of 250,000 infected individuals and more than 35,000 deaths. The Northern regions of Italy were the territories with the large majority of cases (12). A decreasing gradient of infection levels from north to south was clearly evident from the beginning of the epidemic, with atypical differences even between the provinces of the most affected regions (13).

Data published by the Italian National Institute of Statistics (ISTAT) and the Italian National Institute of Health (ISS) on 3 August 2020 demonstrated that the number of people who had been infected with the virus was six-fold higher than the total number of individuals who had been identified by positive test results for both nasopharyngeal and oropharyngeal swabs. The seroprevalence survey on SARS-CoV-2 conducted at national level showed that 1 million 482,000 individuals tested positive for anti-SARS-CoV-2 IgG antibodies, corresponding to 2.5% of the resident population. Data analysis confirmed the important differences observed for spread of the infection between the northern and southern regions of Italy, with the Islands of Sardinia and Sicily having the lowest seroprevalence of anti-SARS-CoV-2 IgG antibodies (0.3%) (14). The different routes of transmission as well as the clinical, epidemiological and virologic features of COVID-19 still need to be clarified today, at more than 6 months from the start of the pandemic (15, 16).

While some traditional or home remedies may alleviate symptoms of mild COVID-19, there are still no drugs or treatments that have been shown to prevent or cure the disease. An initially high viral load has been shown to be independently associated with disease severity and could be influenced by host immune responses (17). Based on the experience gained during the previous Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS) epidemics, it is likely that both innate and adaptive host immunity play a role in viral clearance, disease severity and the different clinical manifestations of the disease (18).

Recently published data suggest that appropriate innate and adaptive T-cell-mediated humoral and cellular immune responses could have a critical role in the elimination of the SARS-CoV-2 virus and confer consequent protective immunity, which in most cases coincides with clinical recovery (11, 19, 20). On the other hand, an excessive cell-mediated and dysregulated innate and adaptive immune response can lead to an aggressive inflammatory reaction with the release of large amounts of pro-inflammatory cytokines. This so-called “cytokine storm”—resulting from the excessive production of cytokines by immune cells such as the innate dendritic cells, macrophages, natural killer (NK) cells and the adaptive T and B cells—directly correlates with lung injury, ARDS and MOF including the kidneys and the central nervous system and overall leads to an unfavourable prognosis (5, 11, 2126).

The human leukocyte antigen (HLA) class I and II molecules are actively engaged in immune response mechanisms against invading pathogens. Their pivotal role in viral immune response has already been highlighted in studies performed on the two coronaviruses responsible for SARS and MERS. Also, in SARS-Cov-2 infection, there is evidence of how different HLA alleles of the major histocompatibility complex (MHC) can define individual susceptibility to infection and its spread (19).

Previous studies on different viruses have shown a correlation between the susceptibility and/or severity of disease and the genetic variability of HLA alleles. These alleles are critical components of the viral antigen presentation pathway. In fact, a study carried out on transgenic mice has shown that the DRB1*04:01 allele, which confers susceptibility to certain autoimmune diseases, generates a robust toll-like receptor (TLR) response and eliminates H1N1 infection. Furthermore, following vaccination and exposure to the H1N1 influenza strain, these mice exhibited cross protective immunity also for the H3N2 strain (27).

Interesting studies conducted on the SARS-CoV-1 virus have identified HLA polymorphisms associated with the disease risk in the Asian population. The HLA-B*46:01 allele was observed with a high frequency in a group of patients from Taiwan classified as being “probable SARS cases” from coronavirus infection. This allele was also significantly increased in the group of seriously ill patients (28). In Chinese patients, the HLA-B*07:03 and HLA-DRB1*03:01 alleles closely correlated with SARS-CoV-1 infection. Moreover, the HLA-B*07:03 and HLA-B*60 alleles were observed with a significantly higher frequency in patients in comparison to the expected frequency of the general population (29). In a group of SARS-CoV-1 positive Vietnamese patients, a strong association was observed for HLA-DRB1*12:02 in comparison to the control group (30). Also, the HLA-C*08:01 allele seemed to confer susceptibility to SARS-CoV-1 (31).

Immunogenetic variation in humans is becoming an increasingly important target for clinical diagnosis and therapeutic intervention. Binding between peptide epitopes and HLA proteins significantly contributes to cellular immune response mechanisms in human beings. Several studies describe the pivotal role of peptides in the specificity, magnitude and quality of both humoral and cellular immune responses. This is furthermore supported by the recent emphasis put on the use of peptides in vaccine design and medical diagnostics. In silico studies have greatly facilitated analysis of the binding affinity between all the viral peptides of SARS-CoV-2 and different HLA class I genotypes. It has been shown that the HLA-B*46:01 allele has a low degree of binding affinity, suggesting that subjects with this allele may have a higher risk of developing the more severe forms of COVID-19. On the other hand, the HLA-B*15:03 allele is reported as having the highest binding affinity for viral peptides (32).

Another study describes high binding affinity between SARS-CoV-2 epitopes and the HLA-A*02:06, HLA-B*52:01, and HLA-C*12:02 alleles. In particular, two epitopes displayed strong binding affinity for HLA-A*24:02, HLA-A*02:01, and HLA-A*02:06 (33). Unfortunately, these two studies present mathematical predictions that need to be clinically and experimentally evaluated and confirmed in a suitable study population with appropriate immunogenetic characteristics.

Researchers in Italy set up a study to investigate whether specific class I HLA alleles could at least partially explain the huge differences observed for the spread of SARS-CoV-2 infection between northern and southern Italy. They compared HLA allele prevalence retrieved through the Italian Bone-Marrow Donor Registry with the incidence of SARS-CoV-2 infections in the different geographical regions. It emerged that HLA-A*25, B*08, B*44, B*15:01, B*51, C*01, and C*03 positively correlated with the incidence of SARS-CoV-2 infection, while HLA-B*14, B*18, and B*49 showed an inverse correlation. After applying a multiple regression model to eliminate confounding factors, only the HLA-C*01 and –B*44 alleles, which are present with a higher frequency in the northern regions of Italy, remained positively associated with COVID-19 (34).

Sardinia is an autonomous region of Italy and the second-largest Island in the Mediterranean Sea. The Island has one of the lowest infection and mortality rates in Italy and on 7 November 2020 registered 11,412 positive individuals (infections: 7x10-3 Sardinian inhabitants vs 14x10-3 Italian mainland inhabitants) and 249 deaths (mortality: 2.18% Sardinian inhabitants vs 3.81% Italian mainland inhabitants). Moreover, during the second wave of the epidemic, Sardinia continued to be one of the Italian regions with the lowest transmission rate (Rt) for SARS-CoV-2 (Sardinian Rt index = 1.16 vs Italian Rt index 1.7— http://www.salute.gov.it/portale/nuovocoronavirus/).

The geographic isolation of Sardinia throughout history has limited the genetic diversity of its people. The population with its low levels of genetic polymorphism has made important contributions to genome-wide association studies of complex disease traits. Many disease risk alleles are shared between Sardinian and other populations which makes it an ideal location for the identification of immunogenetic factors potentially involved in resistance or susceptibility to SARS-CoV-2 infection.

Based on these premises, we set up a study to investigate the clinical characteristics and genetic traits of Sardinian patients affected by SARS-CoV-2 infection. Particular emphasis was put on the influence of HLA Class I and II alleles on disease susceptibility and the varying severity of the disease course.

Materials and Methods

A cohort of 182 patients were recruited from 1 June to 1 August 2020. All patients had been diagnosed with SARS-CoV-2 by RT-PCR from nasopharyngeal swab. The patients were divided into two groups according to disease severity. Thirty-nine patients had been admitted to the Covid Unit of the SS.Trinità Hospital in Cagliari with moderate or severe disease (group S) and 143 asymptomatic or pauci-symptomatic patients (group A) were confined to home quarantine. According to WHO classification, patients with severe disease were considered to be those who needed invasive mechanical ventilation or high-flow nasal oxygen in the hospital, while patients classified as having moderate symptoms did not require oxygen. Pauci-symptomatic patients presented mild symptoms such as loss of taste or smell and/or cold or flu-like symptoms (35). Six hundred and nineteen healthy controls, without SARS-CoV-2 infection (negative for RT-PCR from nasopharyngeal swab), were extracted from the Sardinian Voluntary Bone Marrow Donor Registry, which is highly representative of the genetically homogeneous island population of Sardinia, Italy (36, 37). The controls were selected to appropriately represent the male-to-female ratio and genetic frequencies of the population pertaining to the areas of central and southern Sardinia from where the COVID-19 patients were recruited.

The HLA class I and class II allele frequencies (HLA-A, -B, -C, -DRB1) observed in the two groups of patients (Group S and Group A) were compared to those of the 619 unrelated healthy controls.

Ethics Statement

Patients were recruited and enrolled in the study protocol at the Department of Medical Sciences and Public Health of the University of Cagliari, the University Hospital of Cagliari (AOUCA) and the SS.Trinità Hospital of the Sardinian Regional Company for the Protection of Health (ATS Sardegna). Written informed consent was obtained from all patients and controls in accordance with the ethical standards (institutional and national) of the local human research committee. The study protocol, including informed consent procedures, conforms to the ethical guidelines of the Declaration of Helsinki and was approved by the responsible ethics committee (Ethics Committee of the Cagliari University Hospital; date of approval: May, 27, 2020; protocol number GT/2020/10894). Records of written informed consent are kept on file and are included in the clinical record of each patient.

HLA Allele Typing

DNA from nasopharyngeal swab was extracted with the Qiagen QIAamp DNA Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s instructions. Briefly, each swab was washed with PBS (200 µl). The recovered washing buffer was mixed with 20 µl of proteinase k and 200 µl of the AL buffer was added to the mix. After 10 min of incubation at 56°C, the sample was processed according to the kit instructions. Purity and quantification of total DNA were assessed with the NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific).

Patients and controls were typed at high resolution for the alleles at the HLA-A, B, C, and DRB1 loci using a polymerase chain reaction sequence-specific primer (PCR-SSP) method according to the manufacturer’s instructions (Allele-specific PCR-SSP kits: Olerup SSP AB, Stockholm, Sweden).

HLA typing validation was performed by Next Generation Sequencing on 96 randomly selected samples of patients and controls, according to a previously reported method (38). The concordance of the two different methods was 100%.

HLA Haplotype Analysis in Patients and Controls

Four-loci HLA haplotypes are sequences of four HLA alleles localized in four different HLA loci: HLA-A*ai, HLA-B*bj, HLA-C*ch, HLA-DRB1*dk, where ai, bj, ch, dk are integers (or couples of integers in high resolution) and each index i, j, h, k = 1,2 indicates one of the two alleles in the respective HLA-locus. The four-loci HLA haplotypes in each subject are 24 = 16, obtained by combining four alleles localized in four different HLA loci. Analogously, the three-loci HLA haplotypes in each subject are 23 = 8, obtained by combining three alleles localized in three different HLA loci. Finally, the two-loci HLA haplotypes are 22 = 4, formed by two alleles localized in two different HLA loci.

The frequency of each HLA haplotype in a sample of size is computed by dividing the number of times the HLA haplotype appears in that sample by the total number of alleles for each HLA locus.

The HLA analysis in patients and controls was performed using a specific programming code written with R language (R core39). The statistically significant outcomes given by the R code were also checked “by hand” through the direct evaluation of the number of HLA haplotypes in each group of subjects.

Statistical Analysis

The clinical and genetic characteristics of patients with SARS-CoV-2 infection were analyzed using R software version 4.0.2 (39). Mean standard deviation (SD) and interquartile range (IQR) were calculated for all continuous variables; percentages and 95% confidence intervals (95% CI) were computed for categorical data. P values were calculated using two-tailed Fisher’s exact test or Student’s test as appropriate. Only P values below 0.05 were considered to be statistically significant.

To compare the HLA allele and haplotype frequencies of COVID-19 patients with those of the healthy controls we used the two-tailed Fisher’s exact test. For each single HLA allele, Pc was obtained by multiplying the uncorrected P value by the number of alleles observed. The frequencies of the HLA extended (four loci) and partial (three or two loci) haplotypes were obtained analytically using a programming code created with R language.

Cochran’s rule—claiming that a minimum expected frequency of 5 can be regarded as adequate in analyzing tables with more than a single degree of freedom—was used to establish which HLA haplotypes had frequencies sufficiently high to allow for significant comparisons between the two groups of patients and controls.

For each HLA haplotype the corrected P value (Pc) was calculated by multiplying the P value obtained with the two-tailed Fisher’s exact test by the number of tested allelic combinations. Only Pc values lower than 0.05 were considered to be statistically significant (see Supplementary Material).

Results

Main Clinical and Genetic Characteristics of Sardinian COVID-19 Patients

The most relevant clinical and genetic features of patients positive for SARS-CoV-2 are shown in Table 1. Mean age at the time of infection was 53 years (mean ± SD: 53.2 ± 18.1; 95% CI 50.6–55.9; IQR = 27.4), with a prevalence of female patients (61.5%). Moreover, 53% (n = 96) of the patients had an age at onset of ≤ 50 years; 21% (n = 38) of the patients were over 65 years of age. In line with previous clinical studies, the most severe symptoms were present in adults aged 65 years and older [OR 10.1 (95% CI 4.2–25.2), P = 1.5 10-8]. The clinical manifestations were less severe in female patients (n =112). In fact, 28 of the 39 critically ill patients were males [71.8% vs 28.2%, OR 6.1 (95% CI 2.6–14.8), P = 2.2 10-6].

TABLE 1
www.frontiersin.org

Table 1 Comparisons of baseline clinical, genetic and biochemical parameters between COVID-19 patients with asymptomatic/pauci-symptomatic and moderate/severe disease.

Approximately 13% of the examined subjects had been vaccinated against the flu in the last year. This percentage rose to 15% when the 3-year period 2017-2020 was considered. None of the 39 patients with moderate or severe disease had received anti-flu vaccination during this 3-year period, unlike the group of patients with pauci-symptomatic infection or mild symptoms [0% vs 16.1%, OR 0.07 (95% CI 0–0.5), P = 0.002]. During the winter season in Sardinia of 2019, about 14.2% of the general population underwent flu vaccination. An annual flu vaccination is strongly recommended for people at risk of severe outcomes and in 2019 was administered to 46,5% of the population over 65 years of age (40). It is interesting to note that none of the 22 seriously ill patients, who were over the age of 65, had received a flu vaccination in recent years.

About 8% of the subjects examined had ischemic heart disease, 15% arterial hypertension, 10% hypercholesterolemia, 12% autoimmune diseases and 3.3% type 1 diabetes mellitus. Comorbidities were present in both groups of patients. The only significant difference was a higher frequency of autoimmune diseases (rheumatoid arthritis, type I diabetes and autoimmune hepatitis) in group S [28.2% vs 7.7%, OR 4.7 (95% CI 1.9–11.9), P = 0.001]. Type I diabetes mellitus was present with a higher frequency in the group of patients with a severe clinical picture (5.1% vs 2.8%), but considered separately, did not reach statistical significance.

Five percent (5%) of group S subjects took steroidal anti-inflammatory drugs, 6% non-steroidal anti-inflammatory drugs, 10.4% ACE II inhibitors and/or blockers of angiotensin II receptors (5.5%), 13.7% beta blockers or calcium channel blockers. No statistically significant differences were observed when comparing chronic drug intake between the two groups.

Additionally, two very common genetic traits were considered in the Sardinian population: the beta thalassemic trait and the G6PDH enzyme deficiency. Nineteen COVID-19 positive patients (10.4%) were heterozygous for the β0-39 mutation of the beta globin chain. None of the group S patients carried this mutation [0% vs 13.3%, OR 0.2 (95% CI 0.1–0.8), P = 0.015].

Twenty-four patients (13.2%) had G6PDH enzyme deficiency, resulting in a markedly higher frequency in group S patients compared to group A patients [25.6% vs 9.8%, OR 3.2 (95% CI 1.3–7.9), P = 0.015]. The gene encoding G6PDH maps to the long arm of the X chromosome. There are different genetic variants of G6PDH depending on the geographical area. In Sardinia the enzyme deficiency is caused by the Mediterranean variant (G6PDHMed) which is characterized by the absence or marked reduction of enzymatic activity (0%–10% compared to normal enzyme activity) (41). It should be noted that all 10 patients (25.6%) with G6PDH deficiency of the 39 patients of group S were male and therefore lacked or had markedly reduced enzyme activity.

Allele and Haplotype Frequencies Associated With SARS-CoV-2 Infection

The HLA allele frequencies are illustrated in Table 2 and the Supplementary Material. The 364 alleles of the 182 patients infected by SARS-CoV-2 were compared to the 1,238 alleles of the 619 unrelated healthy controls representative of the HLA frequencies in the Sardinian population. No substantial differences were observed for the HLA Class I (HLA-A, HLA-B, HLA-C) and Class II (HLA-DRB1) alleles. Most of the significant differences (P value < 0.05) were lost after adjustment for multiple comparisons (Pc). The only allele that maintained significance was HLA-C*04:01 [OR = 1.8 (95% CI: 1.3–2.4), P = 0.001; Pc = 0.012]. The HLA-B40 antigen, which initially resulted to be significantly associated to SARS-Cov-2 (Table 2), only exhibited the single allele HLA-B*40:02 in the control group, while in patients it was present with two allelic variants: HLA-B*40:02 and HLA-B*40:01. Therefore, even if the HLA-B40 antigen remained weakly significant after correction for multiple-testing, the alleles lost statistical significance (P = 0.016, Pc = 0.24).

TABLE 2
www.frontiersin.org

Table 2 Comparisons of HLA Class I and Class II allele frequencies between patients and controls and patients with an asymptomatic/pauci-symptomatic or moderate/severe disease course.

More important results were obtained for the haplotypes that characterize the Sardinian population (Table 3). Rather surprisingly, none of the patients affected by COVID-19 carried the HLA-A*02:05, B*58:01, C*07:01, DRB1*03:01 extended haplotype [OR 0.1 (95% CI 0–0.6), P = 0.0006, Pc = 0.015]. After breaking down this extended haplotype into partial haplotypes with three or two loci, deriving from the different allelic combinations, it emerged that no patient infected with SARS-CoV-2 had the three-loci extended haplotype HLA-A*02:05, B*58:01, DR*03:01 [OR 0.1 (95% CI 0–0.5), P = 0.0007, Pc = 0.016]. Also, the HLA-A*02:05, B*58:01, C*07:01 three-loci haplotype had a significantly reduced frequency in patients compared to controls [OR 0.4 (95% CI 0.2–0.7), P = 0.0004, Pc = 0.011].

TABLE 3
www.frontiersin.org

Table 3 Comparison of HLA haplotype frequencies between controls and patients.

Conversely, the HLA-A*30:02, B*14:02, C*08:02 three-loci haplotype was overexpressed in patients compared to the controls [OR 3.8 (95% CI 1.8 –8.1), P = 0.0008, Pc = 0.025].

Overall, the analysis of the HLA alleles and haplotype frequencies revealed seven HLA alleles or haplotypes with a protective effect (OR <1) against SARS-CoV-2 infection and five alleles or haplotypes that instead were associated with an increased susceptibility to infection (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1 The HLA alleles and haplotypes which confer susceptibility to or protection against SARS-CoV-2 infection. The figure shows the odds ratios (OR) for the significantly different frequencies observed for HLA alleles and haplotypes (HLA-B*40:02; HLA-C*04:01; HLA-A*30:02, -B*14:02; HLA-A*02:05, -B*58:01; HLA-A*30:02, -C*08:02; HLA-A*02:05, -C*07:01; HLA-A*02:05, -B*58:01, -C*07:01; HLA-A*30:02, -B*14:02, -C*08:02; HLA-A*02:05, -B*58:01, -DRB1*03:01; HLA-A*02:05, -B*18:01, -DRB1*16:01; HLA-A*02:05, -B*58:01, -C*07:01, -DRB1*03:01 and HLA-A*02:05, -B*58:01, -C*07:01, DRB1*16:01) in the two groups of COVID-19 patients and controls. The error bars, computed according to the two-tailed Fisher’s exact test, represent the 95% confidence intervals of the OR. The red line represents the threshold OR = 1 at which the HLA allele and haplotype frequencies in the two groups of patients and controls are equal.

Correlation Between HLA Allele and Haplotype Frequencies and Severity of Clinical Manifestations in SARS-CoV-2 Infection

Comparisons of the allele frequencies between the two groups A and S are shown in Table 4 and the Supplementary Material. The alleles that remained significant after correction for multiple comparisons were HLA-A*23:01 and HLA-DRB1*08:01. These two alleles were exclusively present in patients with a moderate or severe disease course [OR > 2.5 (95% CI 2.7–220.6), P = 0.002, Pc = 0.038 and OR > 2.5 (95% CI 2.7–220.6), P = 0.002; Pc = 0.024, respectively]. However, it must be considered that HLA-A*23:01 is an uncommon allele in the Sardinian population for which it is difficult to assess its true impact on the evolution of the disease. Furthermore, the association with the severe form of the disease could also be distorted by the limited number of subjects who carried this allele and the fact that they all belonged to the same family nucleus.

TABLE 4
www.frontiersin.org

Table 4 Comparison of HLA allele and haplotype frequencies between patients with pauci-symptomatic/mild (Group A) and moderate/severe disease (Group S).

Only the HLA-A*30:02, B*14:02, C*08:02 three-loci haplotype maintained statistical significance after correction of the P values (Supplementary Material). This haplotype strongly correlated with disease severity [OR 10.4 (95% CI 3.2–34.1), P = 0.0007; Pc = 0.008].

Discussion

Investigation of the SARS-CoV-2 infection confirmed that patients over 65 years of age, particularly males, with comorbidities or organ-associated pathologies are at a higher risk of developing a severe, critical or even fatal disease course (42). Consistent with the findings of other authors, type I diabetes mellitus was more common in patients with the severe clinical manifestations of SARS-CoV-2. Concomitant conditions such as hypertension and cardiovascular disease, obesity and/or a pro-inflammatory and pro-coagulative state are all likely to increase the risk of a worse outcome (43). The fact that type I diabetes mellitus was not significantly associated to a worse prognosis in our study can be attributed to the limited number of patients affected by this comorbidity and evidence of a good glucometabolic control. Overall, our findings suggest that autoimmune diseases (rheumatoid arthritis, autoimmune hepatitis, and type 1 diabetes) are associated with the most severe cases.

In autoimmune disorders, impaired regulation of cell-mediated immune response mechanisms triggers an exaggerated release of proinflammatory cytokines and chemokines by T lymphocytes leading to the so-called “cytokine storm” that often complicates the disease course of COVID-19 (5). Indeed, it appears that SARS-CoV-2 positive patients receiving immunosuppressive treatment for an autoimmune pathology have a much lower risk of developing severe clinical manifestations, and so far no complications such as acute respiratory distress syndrome (ARDS) have been reported in these patients (44, 45).

Similar to the data from Wuhan (46), Lombardy (47) and the USA (48), also in our study the female gender was associated with less severe clinical symptoms. In fact, only 28% of the patients who required oxygen therapy or mechanical ventilation were female. Female sex hormones can only partially explain this protection against SARS-CoV-2 infection. An important role can certainly be attributed to the X chromosome and, in particular, to the conservation of interacting gene clusters related to immune response (49). These genes (microRNA, TLR-7, TLR-8, the chemokine receptor CXCR3 and the interleukin-2 receptor subunit gamma IL2RG etc.) are more highly expressed in women than in males (50). Chromosome X inactivation is random in normal females for which the X chromosome inherited from the mother is active in some cells and the X chromosome inherited from the father is active in others. However, some genes on the X chromosome escape X-inactivation and may at least partially explain why women develop higher innate, humoral and cellular immune responses against pathogens and viral infections (49).

Human Leukocyte Antigen (HLA) molecules play a critical role in the battle between host innate and adaptive (cell-mediated and humoral) immune response mechanisms and invading pathogens. Epitope specificity of HLA Class I molecules creates diversity in innate immune response mechanisms and sees NK cells among the main protagonists. This lymphocyte subpopulation regulates its activity through binding of HLA class I molecules with specific activating and inhibitory receptors (killer cell immunoglobulin-like receptors, NKG2D etc.) and thereby mediates innate defense strategies against viral infections (51, 52).

The HLA complex also has a central role in the recognition and presentation of viral antigens to the immune system (CD4+ and CD8+ lymphocytes) generating a wide range of more or less effective responses. The degree of affinity of the different HLA molecules for specific viral peptides of coronaviruses, such as SARS-CoV and MERS-CoV, drastically influences immune response mechanisms and the consequential clinical manifestations of infection (2833).

The HLA is a highly polymorphic system with profound differences in different geographical areas. For this reason, association studies in the literature often report different and sometimes conflicting results. In the study conducted by Ellinghaus and colleagues, analysis of the classical HLA loci (53) was performed to identify which variants could be correlated to the clinical variability of the COVID-19 disease course. All the HLA allelic variants, observed in 750 patients from Lombardy (northern Italy region) and 850 Spanish patients, lost significance after appropriate statistical correction.

Conversely, a study including about 500,000 voluntary donors enrolled in the Italian Bone Marrow Donor Registry (IBMDR) and representative of the entire Italian peninsula, compared the distribution of certain genes crucial to functioning of the immune system to examine whether they could underlie the geographic differences in COVID-19 incidence. According to the results, the prevalence of the HLA-C*01 and HLA-B*44 alleles are predictive of the incidence of COVID-19 in the different Italian regions and therefore confer susceptibility to SARS-CoV-2 infection (34). This association was not observed in the Sardinian population, probably because the HLA-C*01 and HLA-B*44 alleles, while being fairly common in the north eastern areas of Italy, are scarcely represented in the central-south (34). In other populations, the HLA-B*44 allele would seem to exert a protective effect (54), as opposite to what was observed by Correale et al. (34). In our study of the Sardinian population, we found that the HLA-C*04:01 allele conferred susceptibility to SARS-CoV-2 infection.

NK cell function is regulated by a large number of receptors that transmit either activating or inhibitory signals. Killer cell immunoglobulin-like receptors (KIRs) are among the most important receptors used by NK cells to distinguish viruses and other pathogens from healthy self-counterparts. HLA-C is the dominant ligand for KIR on NK cells. All HLA-C allotypes fall into two major groups of KIR epitopes (HLA-C1 and HLA-C2) on the basis of a dimorphism at position 80 of the alpha domain. The allele found associated to SARS-Cov-2 in our study encodes an HLA-C antigen pertaining to the HLA-C2 group which has high binding affinity for the activating KIR2DS1 receptor and its inhibitory homologue KIR2DL1 (55). Upon binding to HLA-C2, the inhibitory signal generated by KIR2DL1 would seem to override target cell-induced activating signals via KIR2DS1. It is plausible that the higher frequency of HLA-C*04:01 molecules observed in our SARS-CoV-2 positive patients creates an imbalance between HLA molecules of group C1 in favor of those of HLA-C2. This could cause a reduction in some of the functional KIR-HLA ligand interactions with a consequent decrease in the ability of NK cells to control infection.

One of the most interesting findings emerging from the analysis of HLA in Sardinian SARS-CoV-2 patients was the total absence of the extended haplotype HLA-A*02:05, B*58:01, C*07:01, DRB1*03:01. This HLA extended haplotype is third in order of frequency in the Sardinian population (36) and would therefore seem to confer protection against SARS-CoV-2 infection. The three-loci haplotype HLA-A*02:05, B*58:01 DRB1*03:01, derived from this extended haplotype, was also completely absent in SARS-CoV-2 patients whereas the other three-loci haplotype HLA-A*02:05, B*58:01 C*07:01 was present with a significantly lower frequency in patients compared to the general Sardinian population. However, patients carrying the latter three-loci haplotype only had mild disease. Therefore, it is likely that both these three-loci haplotypes exert a protective effect against the severe clinical manifestations of SARS-CoV-2 infection.

The extended haplotype HLA-A*02:05, B*58:01, C*07:01, DRB1*03:01 is extremely rare in other populations where it has a frequency of <0.1% (http://www.allelefrequencies.net/). Therefore, it will not be possible to confirm this result in other geographical areas. The Sardinian population still today represents a genetic isolate with a lowly polymorphic HLA system. This low level of polymorphism makes it easier to analyze the role of HLA molecules in the onset and progression of diseases that are strongly influenced by the immune system.

It is interesting to note how the protective effect against SARS-CoV-2 infection is provided by the extended HLA haplotype and not the individual component alleles. This suggests that combinations rather than single HLA class I and II alleles are able to generate interactions within the complex network of the immune system that allow for an efficient fight against infection.

On the contrary, we observed a greater susceptibility to SARS-CoV-2 infection and a severe disease course in subjects who possessed the HLA-A*30:02, B*14:02, C*08:02 three-loci haplotype. Previous studies have shown that the two-loci haplotype HLA-B*14:02, C*08:02 is associated with hypersensitivity reactions to non-nucleoside reverse transcriptase inhibitors (in particular nevirapine) which are used in the treatment of HIV (56), while HLA-A*30:02 is part of an extended haplotype commonly found in Sardinia (HLA-A*30:02, B*18:01, C*05:01, DRB1*03:01) that is associated with autoimmune diseases such as multiple sclerosis and autoimmune type I hepatitis (37, 57). It is likely that these HLA molecules have low binding affinity for the antigenic peptides of SARS-CoV-2, generating immune responses that cannot effectively clear the virus, but instead overreact causing self-damage.

Finally, we need to discuss the HLA-DRB1*08 (HLA-DRB1*08:01) allele since it appears to be the allele associated with the highest risk for severe clinical manifestations in Sardinian SARS-CoV-2 patients. Because this allele confers susceptibility to numerous autoimmune diseases (58), it can be postulated that the more severe clinical manifestations observed in SARS-CoV-2 patients carrying this allele are caused by altered regulation of cell-mediated immune response (5).

It is clear that the associations between HLA and susceptibility to SARS-CoV-2 infection obtained in our study may be difficult to confirm in other populations with a more polymorphic HLA system. However, they provide valuable information on the binding affinity between viral peptides and MHC class I and II antigens and can contribute to the creation of new antiviral molecules and the selection of vaccines with greater immunogenic potential.

Another surprising finding was that none of the 39 patients with severe disease had been vaccinated against the flu which is strongly recommended, particularly for people over 65 years of age. In Sardinia, the average flu vaccination rate for this category of the population is 46.2%. There is growing consensus that influenza vaccination does not increase the risk of an adverse outcome of COVID-19 but instead may provide some protection against the onset of the more severe or critical clinical manifestations (59). The pathogenesis of influenza viruses and coronaviruses depends upon their ability to dock and enter suitable human host cells. In influenza, the two major internal proteins hemagglutinin (HA) and neuraminidase (NA) mediate virus entry whereas in COVID-19 the spike (S) protein is the leading mediator of virus entry and is a primary determinant of cell tropism and pathogenesis. HA and S proteins are both heavily glycosylated which could make them susceptible to similar components of innate and adaptive immunity (60). Further explanation for a protective effect of the flu vaccine could be that the vaccine for H1N1 seems to have a protective effect on other viruses such as H3N2, thanks to a protective cross-response mediated by Toll Like Receptors through MHC II molecules such as DRB1*04:01 (27).

Another very interesting finding is the higher frequency of G6PDH enzyme deficiency in patients with severe symptoms, compared to pauci-symptomatic patients (25.6% vs 9.8%). It should be noted that the G6PDH deficiency found in 25.6% of the 39 patients with severe clinical pictures was characterized by a marked reduction (10%) in normal enzymatic activity. The importance of this finding should not be underestimated, particularly in light of the high frequency of the G6PDHMed enzymatic variant in Sardinia which fluctuates between 10-20% depending on the coastal or inland areas of the island (61). Although it is rather difficult to explain how the absence of enzymatic activity can promote the onset of severe systemic manifestations in subjects affected by COVID-19, previous studies have demonstrated that virus infection induces production of reactive oxygen species (ROS) and reactive nitrogen species (RNS), which can both cause damage to proteins, DNA and cellular components of cells when antioxidant enzyme metabolism is impaired. Since G6PDH deficiency triggers a redox imbalance in the erythrocytes leading to hemolysis and tissue damage as a result of insufficient oxygen transportation, COVID-19 might increase the mortality risk of patients with this deficit (62, 63).

Another finding worthy of mention was that none of the seriously ill patients were carriers of beta-thalassemia. Two recent reports suggest that this genetic trait is likely to exert a protective effect against COVID-19 since Italian regions with a high number of beta-thalassemia carriers such as Sicily, Puglia, and Sardinia have a very low incidence of SARS-Cov-2 infections (64, 65). This hypothesis is supported by the consideration that ORF8 protein and the surface glycoproteins of SARS-Cov-2 were shown capable of binding to porphyrins. Moreover, computational proteomic analysis revealed that some SARS-CoV-2 proteins could attack the heme on the 1-ß chain of hemoglobin, resulting in the dissociation of iron to form porphyrin SARS-Cov-2 viral proteins (66).

In synthesis, this study highlights numerous clinical and genetic factors that are more or less capable of influencing the course of SARS-CoV-2, all of which can lead to the dreaded severe clinical manifestations of the infection (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2 Factors influencing the clinical course of COVID-19 in Sardinian patients. To establish which factors had a major influence on the clinical disease course, we compared asymptomatic/pauci-asymptomatic patients (Group A) to those with moderate or severe symptoms (Group S). Protective and risk factors were plotted against the odds ratio (OR). Age > 65 years, male gender, autoimmune disease, the HLA-A*30:02, B*14:02, C*08:02 haplotype, G6PDH deficiency, and the HLA-DRB1*08:01 allele conferred an increased risk for severe illness, whereas influenza vaccine (FLU vaccine), the HLA-A*02:05, B*58:01, C07:01, DRB1*03:01 haplotype, female gender, and the beta(0)39-thalassemia trait would seem to offer protection.

The new associations with COVID-19 revealed in the present study deserve investigation on larger cohorts of patients. Meanwhile, the information provided should contribute to the optimisation of treatment and management of patients with the disease. Our results also highlight the importance of influenza vaccination in the battle against COVID-19. During the winter months, SARS-CoV-2 infection can initially easily be confused with the flu and lead to a delayed diagnosis which in turn hampers timely contact tracing procedures to contain further spread of the disease and overload of public health systems. Furthermore, this and other studies suggest that the flu vaccine may offer a certain level of protection against the more severe clinical forms of COVID-19.

However, despite the importance of the information provided, our study only represents a small step ahead in the global effort to understand the disease. Therefore, we strongly recommend analysis on larger samples and diverse populations to further unveil the immunological pathways leading to disease severity and death.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of the Azienda Ospedaliero Universitaria di Cagliari. The patients/participants provided their written informed consent to participate in this study.

Author ContributionS

RL, MC, DF, AP, and LC contributed to study conception and design. RL, MC, SD, GA, DF, SS, AL, FM, DS, WC, MK, APa, MCo, CB, GS, AP and LC especially contributed to acquiring data. RL, SC, MM, SL, SR, RS, RP, MS, PR, GO, SO, AP, LC contributed to analyzing data. RL, MC, SC, MM, AP, LC contributed to interpreting data. RL, MC, SC, MM, MGC, SDG, AR, SmD, SOr, RPe, FMa, DF, AP, and LC contributed to enhancing the intellectual content. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by: Fondazione di Sardegna: grant # 2020.2197 Associazione per l’Avanzamento della Ricerca per i Trapianti O.D.V.

Conflict of Interest

The 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.

Acknowledgments

We are grateful to Anna Maria Koopmans for translations, professional writing assistance and preparation of the manuscript. We wish to thank the Sardinian Foundation (Fondazione di Sardegna) and Associazione per l’Avanzamento della Ricerca per i Trapianti (AART-ODV) for financial support to this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2020.605688/full#supplementary-material

Abbreviations

G6PDH, glucose-6-phosphate dehydrogenase; HA, hemagglutinin; NA, neuraminidase; MOF, multiple-organ failure; ARDS, acute respiratory distress syndrome; ICU, Intensive Care Unit; SARS, Severe Acute Respiratory Syndrome; MERS, Middle East Respiratory Syndrome; HLA, human leukocyte antigen; MHC, major histocompatibility complex; TLR, toll-like receptor; RT-PCR, real time quantitative PCR.

References

1. World Health Organization. Virtual press conference on COVID-19 (2020). Available at: https://www.who.int/docs/default-source/coronaviruse/transcripts/who-audio-emergencies-coronavirus-press-conference-full-and-final-11mar2020.vpdf?sfvrsn=cb432bb3_2 (Accessed August 23, 2020).

Google Scholar

2. World Health Organization. (2020). Available at: https://covid19.who.int/ (Accessed November 10, 2020).

Google Scholar

3. Backer JA, Klinkenberg D, Wallinga J. Incubation period of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China, 20–28 January 2020. Euro Surveill (2020) 25(5):2000062. doi: 10.2807/1560-7917.ES.2020.25.5.2000062

CrossRef Full Text | Google Scholar

4. Lai CC, Shih T-P, Ko W-C, Tang H-J, Hsueh P-R. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges. Int J Antimicrob Agents (2020) 55(3):105924. doi: 10.1016/j.ijantimicag.2020.105924

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ye Q, Wang B, Mao J. The pathogenesis and treatment of the `Cytokine Storm’ in COVID-19. J Infect (2020) 80(6):607–13. doi: 10.1016/j.jinf.2020.03.037

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Kochi AN, Tagliari AP, Forleo GB, Fassini GM, Tondo C. Cardiac and arrhythmic complications in patients with COVID-19. J Cardiovasc Electrophysiol (2020) 31(5):1003–8. doi: 10.1111/jce.14479

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Baig AM. Neurological manifestations in COVID-19 caused by SARS-CoV-2. CNS Neurosci Ther (2020) 26(5):499–501. doi: 10.1111/cns.13372

PubMed Abstract | CrossRef Full Text | Google Scholar

8. I-Cheng L, Huo T-I, Huang Y-H. Gastrointestinal and liver manifestations in patients with COVID-19. J Chin Med Assoc (2020) 83(6):521–3. doi: 10.1097/JCMA.0000000000000319

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Italian Medicines Agency. Agenzia Italiana del Farmaco (2020). Available at: https://www.aifa.gov.it/fr/sperimentazioni-cliniche-covid-19 (Accessed August 23, 2020).

Google Scholar

10. Callaway E. The race for coronavirus vaccines: a graphical guide. Eight ways in which scientists hope to provide immunity to SARS-CoV-2. Nature (2020) 580:576–7. doi: 10.1038/d41586-020-01221-y

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Mutti L, Pentimalli F, Baglio G, Maiorano P, Saladino RE, Correale P, et al. Coronavirus. Disease (COVID-19): What Are We Learning in a Country With High Mortality Rate? Front Immunol (2020) 11:1208. doi: 10.3389/fimmu.2020.01208

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Document from the Italian Civil Protection Agency. Protezione Civile (2020). Available at: http://opendatadpc.maps.arcgis.com/apps/opsdashboard/index.html#/b0c68bce2cce478eaac82fe38d4138b1 (Accessed August 23, 2020).

Google Scholar

13. Livingston E, Bucher K. Coronavirus Disease 2019 (COVID-19) in Italy. JAMA (2020) 323:1335. doi: 10.1001/jama.2020.4344

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Italian Institute of Statistics. Istituto Nazionale di Statistica ISTAT (2020). Available at: https://www.istat.it/it/files//2020/08/ReportPrimiRisultatiIndagineSiero.pdf (Accessed August 3, 2020).

Google Scholar

15. Phelan AL, Katz R, Gostin LO. The novel coronavirus originating in Wuhan, China: challenges for global health governance. JAMA (2020) 323(8):709–10. doi: 10.1001/jama.2020.1097

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Li Q, Guan X, Wu P, Wang X, Lei Z, Tong Y, et al. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N Engl J Med (2020) 382:1199–207. doi: 10.1056/NEJMoa2001316

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Shlomai A, Ben-Zvi H, Glusman Bendersky A, Shafran N, Goldberg E, Sklan EH. Nasopharyngeal viral load predicts hypoxemia and disease outcome in admitted COVID-19 patients. Crit Care (2020) 24(1):539. doi: 10.1186/s13054-020-03244-3

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Janice Oh H-L, Ken-En Gan S, Bertoletti A, Tan Y-J. Understanding the T cell immune response in SARS coronavirus infection. Emerg Microbes Infect (2012) 1:e23. doi: 10.1038/emi.2012.26

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Li X, Gen M, Peng Y, Meng L, Lu S. Molecular immune pathogenesis and diagnosis of COVID-19. J Pharm Anal (2020) 10(2):102–8. doi: 10.1016/j.jpha.2020.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Nicoli F, Solis-Soto MT, Paudel D, Marconi P, Gavioli R, Appay V, et al. Age-related decline of de novo T cell responsiveness as a cause of COVID-19 severity. GeroScience (2020) 42:1015–9. doi: 10.1007/s11357-020-00217-w

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Cheng Y, Luo R, Wang K, Zhang M, Wang Z, Dong L, et al. Kidney disease is associated with in-hospital death of patients with COVID-19. Kidney Int (2020) 95:829–38. doi: 10.1016/j.kint.2020.03.005

CrossRef Full Text | Google Scholar

22. Moriguchi T, Hari N, Goto J, Harada D, Sugawara H, Takamino J, et al. A first Case of Meningitis/Encephalitis associated with SARS-Coronavirus-2. Int J Infect Dis (2020) 95:55–8. doi: 10.1016/j.ijid.2020.03.062

CrossRef Full Text | Google Scholar

23. Wu F, Zhao S, Yu B, Chen Y-M, Wang W, Song Z-G, et al. A new coronavirus associated with human respiratory disease in China. Nature (2020) 579:265–9. doi: 10.1038/s41586-020-2008-3

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet (2020) 395:497–506. doi: 10.1016/S0140-6736(20)30183-5

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Xu Z, Shi L, Wang Y, Zhang J, Huang L, Zhang C, et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir Med (2020) 8(4):420–2. doi: 10.1016/S2213-2600(20)30076-X

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Mangalmurti N, Hunter CA. Cytokine Storms: Understanding COVID-19. Immunity (2020) 53(1):19–25. doi: 10.1016/j.immuni.2020.06.017

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Luckey D, Weaver EA, Osborne DG, Billadeu DD, Taneia V. Immunity to Influenza is dependent on MHC II polymorphism: study with 2 HLA transgenic strains. Sci Rep (2019) 9:19061. doi: 10.1038/s41598-019-55503-1

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lin M, Tseng H-K, Trejaut JA, Lee H-L, Loo H-L, Chu C-C, et al. Association of HLA class I with severe acute respiratory syndrome coronavirus infection. BMC Med Genet (2003) 4:9. doi: 10.1186/1471-2350-4-9

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Ng MHL, Lau KM, Li L, Cheng S-H, 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:515–18. doi: 10.1086/421523

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Keicho N, Itoyama S, Kashiwase K, Phi NC, Long HT, Ha D, 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

31. Chen Y-M, Liang S-Y, Shih Y-P, Chen C-Y, Lee Y-M, Chang L, et al. Epidemiological and Genetic Correlates of Severe Acute Respiratory Syndrome Coronavirus Infection in the Hospital with the Highest Nosocomial Infection Rate in Taiwan in 2003. J Clin Microbiol (2006) 44(2):359–65. doi: 10.1128/JCM.44.2.359-365.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Nguyen A, David JK, Made 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

33. Kazuma K, Yujiro T, Kensaku N, Yusuke N. Bioinformatic prediction of potential T cell epitopes for SARS-Cov-2. J Hum Gen (2020) 65:569–75. doi: 10.1038/s10038-020-0771-5

CrossRef Full Text | Google Scholar

34. Correale P, Mutti L, Pentimalli F, Baglio G, Saladino RE, Sileri P, et al. HLA-B*44 and C*01 Prevalence Correlates with Covid19 Spreading across Italy. Int J Mol Sci (2020) 21:5205. doi: 10.3390/ijms21155205

CrossRef Full Text | Google Scholar

35. Clinical Management of COVID-19. Available at: https://apps.who.int/iris/bitstream/handle/10665/332196/WHO-2019-nCoV-clinical-2020.5-eng.pdf?sequence=1&isAllowed=y (Accessed Access May, 27 2020). Interim guidance, 27 May 2020. World Health Organization. WHO reference number: WHO/2019-nCoV/clinical/2020.5.

Google Scholar

36. Contu L, Arras M, Carcassi C, La Nasa G, Mulargia M. HLA structure of the Sardinian population: a haplotype study of 551 families. Tissue Antigens (1992) 40(4):165–74. doi: 10.1111/j.1399-0039.1992.tb02041.x

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Melis M, Littera R, Cocco E, Frau J, Loi S, Congeddu E, et al. Entropy of human leukocyte antigen and killer-cell immunoglobulin-like receptor systems in immune-mediated disorders: A pilot study on multiple sclerosis. PloS One (2019) 4(12):e0226615. doi: 10.1371/journal.pone.0226615

CrossRef Full Text | Google Scholar

38. Closa L, Vidal F, Herrero MJ, Caro JL. Design and Validation of a Multiplex KIR and HLA Class I Geno-typing Method Using Next Generation Sequencing. Front Immunol (2018) 9:2991. doi: 10.3389/fimmu.2018.02991

PubMed Abstract | CrossRef Full Text | Google Scholar

39. R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing (2013). Available at: http://www.R-project.org/.

Google Scholar

40. Italian Ministry of Health. MINISTERO DELLA SALUTE - ISS (2020). Available at: https://www.epicentro.iss.it/influenza/stagione-in-corso (Accessed August).

Google Scholar

41. Beutler E. Glucose-6-phosphate dehydrogenase deficiency: a historical perspective. Blood (2008) 111(1):16–24. doi: 10.1182/blood-2007-04-077412

PubMed Abstract | CrossRef Full Text | Google Scholar

42. 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

43. Apicella M, Campopiano MC, Mantuano M, Mazoni L, Coppelli A, Del Prato S. COVID-19 in people with diabetes: understanding the reasons for worse outcomes. Lancet Diabetes Endocrinol (2020) 8(9):782–92. doi: 10.1016/S2213-8587(20)30238-2

PubMed Abstract | CrossRef Full Text | Google Scholar

44. D’Antiga L. Coronaviruses and Immunosuppressed Patients: The Facts During the Third Epidemic. Liver Transpl (2020) 26(6):832–4. doi: 10.1002/lt.25756

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lleo A, Invernizzi P, Lohes AW, Aghemo A, Carbone M. Management of patients with autoimmune liver disease during COVID-19 pandemic. J Hepatol (2020) 73(2):453–5. doi: 10.1016/j.jhep.2020.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinic features of patients infected with 2019 novel corona- virus in Wuhan, China. Lancet (2020) 395(10223):497–506. doi: 10.1016/S01406736(20)30183-5

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Grasselli G, Zangrillo A, Zanella A, Antonelli M, Cabrini L, Castelli A, et al. Baseline Characteristics and Outcomes of 1591 Patients Infected With SARS-CoV-2 Admitted to ICUs of the Lombardy Region, Italy. JAMA (2020) 323(16):1574–81. doi: 10.1001/jama.2020.5394

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Richardson S, Hirsch JS, Narasimhan M, Crawford JM, McGinn T, Davidson KW, et al. Presenting Characteristics, Comorbidities, and Outcomes Among 5700 Patients Hospitalized With COVID-19 in the New York City Area. JAMA (2020) 323(20):2052–9. doi: 10.1001/jama.2020.6775

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Schurz H, Salie M, Tromp G, Hoal EG, Kinnear CJ, Möller M. The X chromosome and sex-specific effects in infectious disease susceptibility. Hum Genomics (2019) 13(1):2. doi: 10.1186/s40246-018-0185-z

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Carrel L, Willard HF. X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature (2005) 434(7031):400–4. doi: 10.1038/nature03479

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Parham P, Guethlein LA. Genetics of Natural Killer Cells in Human Health,Disease, and Survival. Annu Rev Immunol (2018) 26(36):519–48. doi: 10.1146/annurev-immunol-042617-053149

CrossRef Full Text | Google Scholar

52. Kulkarni S, Martin MP, Carrington M. The Yin and Yang of HLA and KIR in human disease. Semin Immunol (2008) 20(6):343–52. doi: 10.1016/j.smim.2008.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Ellinghaus D, Degenhardt F, Bujanda L, Buti M, Albillos A, Invernizzi P, et al. Genomewide Association Study of Severe Covid-19 with Respiratory Failure. N Engl J Med (2020) 382:727–33. doi: 10.1056/NEJMoa2020283

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Poulton K, Wright P, Hughes P, Savic S, Welberry Smith M, Guiver M, et al. A role for human leucocyte antigens in the susceptibility to SARS-Cov-2 infection observed in transplant patients. Int J Immunogenet (2020) 47(4):324–8. doi: 10.1111/iji.12505

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Brooks AG, Boyington JC, Sun PD. Natural killer cell recognition of HLA class I molecules. Rev Immunogenet (2000) 2(3):433–48.

PubMed Abstract | Google Scholar

56. Littera R, Carcassi C, Masala A, Piano P, Serra P, Ortu F, et al. HLA-dependent hypersensitivity to nevirapine in Sardinian HIV patients. AIDS (2006) 20(12):1621–6. doi: 10.1097/01.aids.0000238408.82947.09

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Littera R, Luchino C, Simona O, Francesco F, Sara L, Luca S. Exploring the Role of Killer Cell Immunoglobulin-Like Receptors and Their HLA Class I Ligands in Autoimmune Hepatitis. PloS One (2016) 11(1):e0146086. doi: 10.1371/journal.pone.0146086

PubMed Abstract | CrossRef Full Text | Google Scholar

58. International MHC and Autoimmunity Genetics Network, Rioux JD, Goyette P, Vyse TJ, Hammarström L, Fernando MM, et al. Mapping of multiple susceptibility variants within the MHC region for 7 immune-mediated diseases. Proc Natl Acad Sci USA (2009) 106(44):18680–5. doi: 10.1073/pnas.0909307106

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Noale M, Trevisan C, Maggi S, Antonelli Incalzi R, Pedone C, Di Bari M, et al. The Association between Influenza and Pneumococcal Vaccinations and SARS-Cov-2 Infection: Data from the EPICOVID19 Web-Based Survey. Vaccines (Basel) (2020) 8(3):E471. doi: 10.3390/vaccines8030471

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Liu Y-M, Mahumud S-A, Chen X, Chen T-H, Liao K-S, Lo J-M. Carbohydrate-Binding Protein from the Edible Lablab Beans Effectively Blocks the Infections of Influenza Viruses and SARS-CoV-2. Cell Rep (2020) 32(6):108016. doi: 10.1016/j.celrep.2020.108016

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Carcassi U, Del Giacco GS, Pintus A, Perpignano G, Locci F, Loy M, et al. HLA antigens and erythrocyte glucose-6-phosphate dehydrogenase (G-6-PD) deficiency in Sardinia. Int Arch Allergy Appl Immunol (1977) 54(3):285–7. doi: 10.1159/000231838

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Wu Y-H, Tseng C-P, Cheng M-L, Ho H-Y, Shih S-R, Chiu DT-Y, et al. Glucose-6-phosphate dehydrogenase deficiency enhances human coronavirus 229E infection. J Infect Dis (2008) 197(6):812–6. doi: 10.1086/528377

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Aydemir D, Ulusu NN. Is glucose-6-phosphate dehydrogenase enzyme deficiency a factor in Coronavirus-19 (COVID-19) infections and deaths? Pathog Glob Health (2020) 114(3):109–10. doi: 10.1080/20477724.2020.1751388

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Lansiaux E, Pébaÿ PP, Picard J-L, Son-Forget J. COVID-19: beta-thalassemia subjects immunised? Med Hypotheses (2020) 142:109827. doi: 10.1016/j.mehy.2020.109827

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Drouin E. Beta-thalassemia may protect against COVID 19. Med Hypotheses (2020) 143:110014. doi: 10.1016/j.mehy.2020.110014

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Wenzhong L, Hualan L. COVID-19: attacks the 1-beta chain of hemoglobin and captures the porphyrin to inhibit human heme metabolism. ChemRxiv (2020). doi: 10.26434/chemrxiv.11938173.v4

CrossRef Full Text | Google Scholar

Keywords: human leukocyte antigen, SARS-CoV-2 infection, immunogenetic background, COVID-19 severity, Sardinian population, alleles, haplotypes, glucose-6-phosphate dehydrogenase

Citation: Littera R, Campagna M, Deidda S, Angioni G, Cipri S, Melis M, Firinu D, Santus S, Lai A, Porcella R, Lai S, Rassu S, Scioscia R, Meloni F, Schirru D, Cordeddu W, Kowalik MA, Serra M, Ragatzu P, Carta MG, Del Giacco S, Restivo A, Deidda S, Orrù S, Palimodde A, Perra R, Orrù G, Conti M, Balestrieri C, Serra G, Onali S, Marongiu F, Perra A and Chessa L (2020) 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. 11:605688. doi: 10.3389/fimmu.2020.605688

Received: 12 September 2020; Accepted: 18 November 2020;
Published: 04 December 2020.

Edited by:

Jukka Partanen, Finnish Red Cross Blood Service, Finland

Reviewed by:

Marja-Liisa Lokki, University of Helsinki, Finland
Riccardo Gavioli, University of Ferrara, Italy

Copyright © 2020 Littera, Campagna, Deidda, Angioni, Cipri, Melis, Firinu, Santus, Lai, Porcella, Lai, Rassu, Scioscia, Meloni, Schirru, Cordeddu, Kowalik, Serra, Ragatzu, Carta, Del Giacco, Restivo, Deidda, Orrù, Palimodde, Perra, Orrù, Conti, Balestrieri, Serra, Onali, Marongiu, Perra and Chessa. 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: Roberto Littera, roby.litter@gmail.com; Marcello Campagna, mam.campagna@gmail.com; Andrea Perra, andrea.perra@unica.it; Luchino Chessa, luchinochessa@unica.it

These authors have contributed equally to this work