Next Article in Journal
Dopamine: A New Player in the Pathogenesis of Diabetic Retinopathy?
Previous Article in Journal
StUBC13, a Ubiquitin-Conjugating Enzyme, Positively Regulates Salt and Osmotic Stresses in Potato
Previous Article in Special Issue
Bulk and Single-Cell RNA Sequencing Elucidate the Etiology of Severe COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sex-Specific HLA Alleles Contribute to the Modulation of COVID-19 Severity

by
Serena Spartano
1,2,†,
Maria Vittoria Faggiano
1,2,†,
Giovanna Guidi
3,4,
Pino D’Ambrosio
1,2,
Alessandro Vaisfeld
1,2,
Agnese Novelli
1,2,
Salvatore Falqui
1,2,
Antonella Cingolani
3,4,
Lorenza Lambertenghi
5,
Alessandro Visentin
5,
Annamaria Azzini
5,
Elda Righi
5,
Enrico Maria Trecarichi
6,
Maria Mazzitelli
6,
Silvano Coletti
7,
Jan Mous
7,
Thomas W. Rademacher
8,
Carlo Torti
3,4,6,
Evelina Tacconelli
5,9,
Massimo Fantoni
3,4,
Roberto Cauda
3,4 and
Francesco Danilo Tiziano
1,2,*
add Show full author list remove Hide full author list
1
Section of Genomic Medicine, Department of Life Sciences and Public Health, Catholic University of Sacred Heart, 00168 Rome, Italy
2
Complex Unit of Medical Genetics, Department of Laboratory and Infectivologic Sciences, Policlinico “A. Gemelli” Foundation, 00168 Rome, Italy
3
Section of Infectious Diseases, Department of Safety and Bioethics, Catholic University of Sacred Heart, 00168 Rome, Italy
4
Unit of Infectious Diseases, Department of Medical and Surgical Sciences, Policlinico “A. Gemelli” Foundation, 00168 Rome, Italy
5
Infectious Diseases Division, Department of Diagnostics and Public Health, University of Verona, 37134 Verona, Italy
6
Infectious and Tropical Diseases Unit, Department of Medical and Surgical Sciences, University “Magna Graecia”—“R. Dulbecco”, 88100 Catanzaro, Italy
7
Chelonia SA, 6900 Lugano, Switzerland
8
Division of Infection and Immunity, University College London and Middlesex University, London WC1E 6BT, UK
9
Clinical Research Unit, German Center for Infectious Diseases, Tübingen University, 72074 Tübingen, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(23), 13198; https://doi.org/10.3390/ijms252313198
Submission received: 26 September 2024 / Revised: 5 November 2024 / Accepted: 20 November 2024 / Published: 8 December 2024
(This article belongs to the Special Issue Genomic Studies Related to SARS-CoV-2 and Other Common Pathogens)

Abstract

:
Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection, responsible for Coronavirus Disease 2019 (COVID-19), exhibits a spectrum of clinical manifestations, ranging from asymptomatic to severe pulmonary dysfunction or death. The variability in COVID-19 severity has largely been attributed to the host’s genetic characteristics, suggesting a polygenic genetic architecture, without significant strong evidence of sex-related genetic differences. In this Italian retrospective case–control study, we investigated the association between COVID-19 severity (severe vs. asymptomatic/oligosymptomatic healed individuals) and HLA gene variants, analyzed by next-generation sequencing (NGS). We identified significant HLA alleles (according to the conventional nomenclature), SNPs and haplotypes in the HLA-B, -C, -F, -DQA1, -DRB1, and -DRB5 genes associated with COVID-19 severity. Interestingly, these variants showed biological sex-related effects. Also, we identified specific haplotypes associated with COVID-19 severity that are shared by different conventional HLA alleles, indicated here as “super-haplotypes”. These haplotypes had a biological sex-specific impact on disease severity and markedly increased the risk of severe COVID-19 compared to the conventional HLA alleles (odds ratio of up to 15). Our data suggest that the revision of the current HLA nomenclature may help to identify variants with a stronger effect on disease susceptibility and that association studies could benefit from the stratification of patients by biological sex. If replicated in other disease models, these findings could help to define the functional diversity in immune response between sexes, also based on the HLA system. Finally, due to the global pandemic’s mortality rate, we hypothesize here that SARS-CoV-2 may have acted as a natural selection trigger, leading to a drift in HLA allelic frequencies in the general population.

1. Introduction

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic has had some prominent characteristics, including extremely high transmissibility, relatively high viral mutability, and variability in individual response to infection [1]. The severity of Coronavirus Disease 2019 (COVID-19) has varied from completely asymptomatic forms to very severe pulmonary dysfunction with respiratory insufficiency or death [2]. The mortality rate of COVID-19 has been very variable in different studies and populations, likely based on viral genotype, standards of care for the clinical assistance of patients, vaccination status and individual characteristics of the host [3]. Among the negative prognostic factors were male biological sex, older age and co-morbidities, as confirmed in a recent meta-analysis [2]. However, co-morbidities correlate more with the risk of death, rather than with disease severity [2].
Regarding the correlation between the genetic background of the host and COVID-19 severity, different case–control genome-wide association studies have been performed. In a recent meta-analysis, 49 SNPs were significantly associated with severe forms of COVID-19 [4].
Some studies have concentrated more on the specific role of the Human Leukocyte Antigen (HLA) system as a severity modulator of COVID-19, identifying some potentially relevant HLA alleles [5]. Some alleles, such as HLA-A02, may reduce the risk of COVID-19, while others, like HLA-C04, could increase severity and mortality [6]. In detail, HLA-C* 04:01 has been identified as a genetic predictor of a severe clinical course of COVID-19, leading to a 2-fold increase in the risk of intubation [7]. In an Italian cohort, carriers of the HLA-B*44 allele had a lower risk of developing severe COVID-19 symptoms [8]; in a Dutch study, the HLA-B*15:03 allele was associated with a reduced risk of infection [3]. Finally, it has been shown that HLA-DRB1*04:01 is more frequent among asymptomatic compared to symptomatic subjects, suggesting a protective role [5]. In a recent review, Hoseinnezhad et al. evaluated the role of HLA genetic variants in COVID-19 susceptibility, severity, and mortality in different populations. These authors suggested that the distribution of HLA alleles in the different populations may account, at least partially, for the diverse outcomes of SARS-CoV-2 infection observed in different geographic areas [6]. To our knowledge, in none of the studies mentioned above were cases and controls stratified according to biological sex.
In the present study, we performed an Italian multicentric case–control study of the association of HLA gene variants with the severity of COVID-19. The study design was prospective–retrospective; subjects healed from COVID-19 were prospectively recruited and retrospectively assigned to case (severe) or control (oligo-asymptomatic) arms. HLA genes (listed in Supplementary Table S1) were genotyped by next-generation sequencing (NGS) on gDNA extracted from blood samples collected on Dried Blood Spot (DBS) cards.

2. Results

2.1. Population Characteristics

We recruited 679 subjects (493 cases), with a case–control ratio of 2.6. The demographic characteristics and inclusion criteria are reported in Table 1. Males represented 67.55% and 42.47% of cases and controls, respectively (p < 0.0001). The median age of the cases was 59.24 years, and that of the controls was 52.55 (p < 0.0001).

2.2. HLA-DRB1 Is the Most Polymorphic Gene

To assess the intrinsic genetic variability of HLA genes in our population, we related number of identified variants and size of the sequenced region per each gene. As shown in Figure 1, HLA-DRB1 displayed the highest number of variants/sequenced base pairs (bp). HLA-A, -B, and -C displayed a similar polymorphic rate; HLA-E, -F, and -G appeared to be the least variable (ranging between 0.018 and 0.043 variants/bp); most class II genes displayed a similar range of variation.

2.3. HLA-A, -B, -C, -DPA1, -DPB1, -DQB1, -DRB1 Alleles Are Associated with COVID-19 Severity in a Biological Sex-Specific Manner

For conventional HLA genotyping and nomenclature (HLA-alleles, ImMunoGeneTics project, IMGT), Ref. [9] data were analyzable for 670 individuals. Results are reported in Table 2 and detailed in Supplemental Materials. The alleles indicated below in HLA-A, -B, -C -DPA1, -DPB1, -DQB1 and -DRB1 were associated with COVID-19 severity.
The stratification of the cohorts by biological sex led to the identification of HLA alleles with a sex-specific effect; the same allele might be indifferent, protective or of risk, in one or the other sex (Table 2), or even non-significant in the whole population (see HLA-C*12 or *15 as examples in Table 2).

2.4. Single HLA-SNPs Are Associated with COVID-19 Severity

We identified 5578 Single-Nucleotide Polymorphisms at the HLA locus (HLA-SNPs), 118 significantly associated with COVID-19 severity (FDR ≤ 0.05). Among the significant variants, we prioritized only those in the coding regions and the intronic predicted to affect exon splicing, irrespective of the allelic frequency. We identified 28 variants that were confirmed by Sanger Sequencing (see Table 3, Supplemental Figure S1 and Supplemental Materials). Among the class I genes, HLA-B, -C, and -F genes were significantly associated with the severe disease phenotype, and among the class II, HLA-DQA1, -DRB1, and -DRB5. The significant HLA-SNPs, their relative odds ratio (OR), 95% confidence intervals (CI) and p-values are reported in Table 3. The OR values for the single variants ranged from 1.02 to 3.62 for risk SNPs, and from 0.42 to 0.77 for the protective.

2.5. Age and Biological Sex Impact on the Significance of HLA-SNPs on COVID-19

Due to the significant difference between cases and controls in age and biological sex, we included these two variables as co-variates in a multinomial logistic regression model; this analysis led to a reduction in the number of significant HLA-SNPs. The remaining are indicated in bold in Table 3. According to this model, HLA-DQA1 variants lost significance.
We analyzed the single or combined effect of each HLA-SNP by comparing allelic and genotypic frequencies, respectively, by chi-squared test. This approach was used to evaluate whether the COVID-19 severity risk was influenced by the homozygous or heterozygous state for each variant. To assess the possible isoallelic influence of the single HLA-SNPs (i.e., if SNPs present in the two alleles might have an additive effect on COVID-19 severity risk), we used a log-additive model. Both analyses were conducted in the whole population of cases or controls, or stratified by biological sex. We finally compared the allelic frequency of each HLA-SNP in males and females, independent if cases or controls (column labeled as M vs. F in Table 4).
From the analysis of the two cohorts, HLA-B, -C, -DRB1 and -DRB5 were significantly associated with COVID-19 severity, whereas HLA-F and -DQB1 lost significance. In most cases (except for few variants), the COVID-19 risk was not influenced by the heterozygous or homozygous state, and we did not find consistent evidence of isoallelic influence for the single SNPs, as suggested by the OR values.
The stratification of the cohorts by biological sex led to two relevant observations:
(1) Risk genes were different in males and females (Table 4).
As an example, HLA-C was associated with the phenotype in males only, HLA-B in females.
(2) Allelic and genotypic frequencies of HLA-SNPs were significantly different in the two sexes when stratified by cases and controls but not in the whole study population, as expected for autosomal loci, ruling out possible bias in patients’ enrollment (Table 4).

2.6. “Super-Haplotypes” Are Shared Among Different HLA Alleles

To evaluate the potential combinatorial cis effect of the significant HLA-SNPs, i.e., whether variants located in the same allele might reciprocally interact, we inferred haplotypes and gametic phases per each HLA gene [10]. Haplotypes were arbitrarily numbered from 1 (haploptype ID #), and are listed in Supplementary Table S2A–C. We identified some haplotypes that were associated with COVID-19 severity (listed in Table 5) in a biological sex-specific manner; the risk haplotypes were different in males or females, or might display opposite behavior in biological sexes.
As examples, HLA-B haplotype#3 and #5 were protective in males, of risk in females, non-significant in the whole population; HLA-C haplotype#9 and #17 were significant in females only (risk haplotypes). Also, in these analyses, the haplotypic frequencies in the whole population were not significantly different between males and females, ruling out possible sample biases.
To identify the conventional HLA alleles bearing each significant haplotype, we aligned the nucleotidic positions characterizing each haplotype with the corresponding positions in the HLA cDNA sequences contained in the HLA-IMGT dictionary. Surprisingly, we found that:
(1)
the same haplotype was present in different HLA-alleles;
(2)
only a small number of haplotypes were found in the dictionary (see Supporting Information).
Following these observations, we propose here to define “super-haplotypes” (HLA-SH), those shared by different HLA alleles.

3. Discussion

Most variability of COVID-19 has been attributed to the genetic characteristics of the host, although the association studies performed so far have not revealed major genetic loci or variants with strong modulating effect, suggesting a polygenic genetic architecture; additionally, there is no clear evidence of biological sex-related differences, except for very rare variants, located in the X-chromosome [11].
To evaluate the specific role of HLA genes in the modulation of COVID-19 severity, we have conceived the present study as retrospective with prospective recruitment of patients healed from COVID-19, with the aim to clearly identify dichotomic groups (a-/oligosymptomatic vs. severe subjects) and maximize the chance of detecting relevant predisposing variants in the HLA locus. On the other hand, this strategy prevented us from studying the extreme side of the spectrum of the disease, the patients who died of COVID-19.
Regarding the choice of the experimental technology, we opted for the NGS of the exons and the flanking intronic regions of HLA genes (excluding known pseudogenes); this approach provided the unquestionable advantages of reducing costs/sample. On the other hand, the huge increase in the number of variants identified, including the intronic, has enormously increased the complexity of result interpretation; it is conceivable that most SNPs are irrelevant to the phenotype or are simply in linkage disequilibrium with others that are functionally relevant. In the present study, to facilitate data interpretation, we focused on variants located in the coding regions only, although some intronic variants were significantly associated with COVID-19 severity. This approach might have led to the loss of variants with a splicing-modulating effect; it is well known that alternative splicing is a common event in the HLA genes, although the biological meaning and the regulation of the inclusion/exclusion of specific exons are largely unknown.
Our data have revealed an unexpected variable degree of polymorphism of each HLA gene; the evolutionary reason for a lesser extent of variability of HLA-E, -F, and -G, up to 8-fold lower than the other class I genes, is unknown (see Figure 1). These characteristics are even more evident in the class II genes, with some genes being extremely variable, more so than HLA-A, -B and -C. We wonder whether this could be related to intrinsic characteristics of the HLA locus or could rather be attributed to a negative selection on the variation of some genes, related to the specific function in the immune system.
Regarding the association with disease severity, we first identified significantly associated HLA alleles and found those reported in Table 2. Among the class I genes, HLA-A, -B, and -C were associated with COVID-19 severity, among the class II genes, HLA-DPA1, -DPB1, -DQB1, -DRB1. The significant alleles conferred a moderate increase/decrease in the OR of severe forms, up to 2-fold. None of these alleles was among the previously identified ones, reported in the literature; this discrepancy could be related to different inclusion criteria, technical differences in genotyping or different allelic frequencies in the studied populations. Differently from previously published studies, we evaluated the biological sex-specific effects of HLA alleles, which led to a consistent improvement in patient risk assessment (Table 2).
Subsequently, we analyzed our dataset using a classical genetic case–control association model. This approach led to the identification of significant HLA-SNPs and to consistently different results than the previous approach; the association of some genes lost significance (HLA-A and all class II genes except -DRB5), whereas other genes became significantly associated with COVID-19, namely HLA-F, -DQA1, and -DRB5. Biological sex stratification led to even more striking results; risk genes differed in males and females according to a quite clear dichotomous model (Table 4).
To address the role of the cis effect of the significant variants in each gene, we first inferred the gametic phase of these variants, and identified the HLA-SH, as indicated in the results. Importantly, these haplotypes have an opposite effect in the two sexes. Only a few had a biological sex-independent impact. The identification of HLA-SH markedly improved the risk assessment compared to the conventional genotyping, since the differential risk of severe forms of COVID-19 was increased up to 5-fold between males and females, as in the case of HLA-SH#1 of the HLA-DQA1 gene that is of high risk in females.
Besides the translational impact of our findings, in terms of individual risk assessment for infected subjects, which may still be useful during the current epidemic phase of COVID-19, there are several lessons we have learned from the present study. First, the identification of HLA-SH strongly suggests the need to revise the standard nomenclature of HLA alleles and the technical approaches to HLA genotyping, through validation in exploratory scenarios.
Our study has also revealed the limitations of the gDNA-based approach in terms of interpretation of the HLA variability in the general population, of the combinatorial effect of cis-variants in the same HLA gene and of possible isoallelic influence of the homologous gene in trans. The most used NGS platforms are based on short read sequencing that do not allow the definition of the gametic phase of the single variants. The spreading of long-read sequencing of gDNA may help address this point and define the combinatorial effect of the different HLA gene variants, in cis on the same chromosome. On the other hand, this approach will be unlikely to help define the role of intronic variants, namely in the alternative splicing modulation, which may play a key role in the immune system modulation. As an example, it has been shown that the skipping of exon 4 of HLA-G leads to the production of a soluble secreted isoform that apparently does not have an antigen presentation function but instead has an immunomodulatory effect, via the regulation of Natural Killer (NK) cells [12]. This phenomenon is less known in other HLA genes but, in principle, a variable splicing pattern can produce isoforms with different functions than those classically attributed to the MHC molecules [13]. We propose implementing, at least in an exploratory phase, HLA genotyping on cDNA from whole blood. To our knowledge, there are few such studies reported in the literature [14,15,16]. Even though this approach does not allow the identification of intronic variants, it provides experimental data on the gametic phase of the coding SNPs and on the qualification and quantification of the alternative mRNA isoform repertoire in vivo.
From the immunology point of view, we provide here evidence of biological sex-specific differences in disease susceptibility related to HLA genes in humans. It is well known that the immune response is different in males and females, and that there are clear differences in the epidemiology of some infectious and autoimmune diseases [17]. Some of the differences in immune response have been very recently attributed to the effect of sex hormones, namely androgens [18]. However, we could not find in the literature studies including the differential analysis of HLA genotype according to biological sex. Similarly, we could not find raw data of published association studies on autoimmune or other HLA-related diseases, including sex-specific allelic frequencies, to perform meta-analyses. It could be of extreme interest to assess whether our findings are replicated in other disease models, beyond COVID-19, since this could shed some additional hints on the functional diversity in the immune system function between males and females.
Our study identified a specific HLA-SNP and haplotype of the non-classical class I gene HLA-F, associated with disease severity. More information is needed on the role of this HLA gene in human disease, even though relatively recent, preliminary data suggest the involvement in viral infection [19]. To our knowledge, HLA-F has been mainly related to the embryo tolerance [20] and is involved in immune system regulation through the modulation of NK response [19]. In our population, HLA-F is among the least variable genes; a single variant, having an allelic frequency of 0.04 in the global population and 0.06 in female cases, confers a 5-fold increased risk of severe COVID-19. The functional role of this variant and, more in general, of HLA-F in disease will require future additional dedicated studies.
As a final remark, from the point of view of the classical genetics of populations, our study has provided novel potentially straightforward insights. As expected, when analyzing our population as a whole (without the stratification by cases or controls), allelic and genotypic frequencies of HLA alleles, HLA-SNPs and HLA-SH are not different between males and females, as expected for autosomal loci in the absence of selective pressure. When stratifying by cases and controls, we observed a statistically significant difference in frequencies between males and females; we can hypothesize an “acute” biological sex-specific natural selection event based on HLA genotype. Even if we could not analyze samples from individuals who died of COVID-19, we can hypothesize the existence of HLA “deadly” alleles that were lost in the general population. If this hypothesis is correct, the huge number of COVID-19 victims worldwide (in the order of several million, underestimated due to the lack of punctual data from many countries) may have led to a consistent drift in the HLA allelic frequencies in the general population. A similar phenomenon may have occurred during the Spanish flu pandemics of 1918–1920, with an estimated number of deaths of 50M people worldwide out of an estimated population of 1.8 billion persons [21].
Among the differences between the two pandemics is that most victims of the Spanish flu were subjects of reproductive age (18–35 years) [21], thus likely changing allelic frequencies by a consistent extent of phenotypically relevant genes over the next generations. In the case of COVID-19, most victims were over reproductive age.
We are well aware that the monogenic prediction model of COVID-19 severity is an over-simplification. However, we wonder whether COVID-19 may have led to the decanalization (i.e., to the perturbation of the gene–environment homeostasis due to large effect-size environmental changes or genetic variants) [22] of the immune response to some environmental triggers, mediated by the HLA locus. A similar mechanism has been evoked to explain the increase in type II diabetes in Greenlandic populations, in response to an increase in carbohydrate intake with diet and the presence of a single TBC1D4 SNP (rs61736969) [23]. To address this point, due to the key role of HLA genes in modulating the susceptibility to most diseases, it will be advisable to monitor possible changes in the epidemiology of immune-mediated conditions over the next few years.
Our study presents some limitations. During the enrollment phase, the epidemiology of the pandemics changed consistently, due to the drift of SARS-CoV-2 viral variants. During the first part of the enrollment, patients may have been mainly infected by the wild-type and alpha strains, while during 2022, the main circulating variant was the BA (Omicron variant) [24]. Due to the healed status of the patients, we could not genotype viral variants in our cohorts; viral variant assessment was not routinely performed at the time of the infection. Additionally, starting from November 2021, the epidemiology of COVID-19 changed when vaccines became available, leading to a marked slowdown in recruitment rate.
In conclusion, our study provides some hints that will likely open new perspectives for further studies in understanding HLA-mediated immune modulation. Starting from this evidence, new scenarios could be opened in understanding HLA-mediated immune modulation in response to environmental HLA effect among biological sexes, beyond the SARS-CoV2 pandemics.
The identification of individual genetic fingerprints, predictive of the severity of disease in the case of exposure to pathogens, might contribute, during interpandemic times, to provide operational procedures for the early therapeutic management protocols and vaccination campaigns.

4. Materials and Methods

4.1. Study Design

We conducted a retrospective case–control association study aimed at investigating the role of HLA variants in COVID-19 severity. Healed subjects were retrospectively assigned to case (severe) or control (oligo-/asymptomatics) arms and prospectively enrolled between 7 December 2020 and 12 May 2022.

4.2. Recruiting Centers

The study was conducted in 3 Italian University Hospitals, located in Rome, Verona and Catanzaro, and was approved by the ethics committees of the three institutions. Cases were recruited among healed patients (subjects who tested positive for the presence of SARS-CoV-2 genome and subsequently tested negative in at least two repeated samplings), in follow-up at the post-COVID outpatients’ units; controls were recruited among the oligo-/asymptomatic relatives.
The enrollment criteria were as follows:
Cases: Hospitalized patients who displayed at least one of the following clinical findings:
  • Oxygen saturation (SaO2) < 94% on room air; <90% in presence of known chronic hypoxic conditions or receiving chronic supplemental oxygen;
  • Respiratory rate > 24 breaths/min.
Controls: Subjects who did not require hospitalization and who displayed few or no signs/symptoms associated with COVID-19, without a relevant impairment in clinical conditions or quality of life.
Sample size and study power were calculated by Genetic Power Calculator [25,26]. We have conservatively considered a minimum allele frequency (MAF) of 5% for each genetic variant in the population and a prevalence of the severe phenotype of about 10%. Assuming a population study of 500 subjects/arm and ratio 1:1 cases:controls, the study would have achieved a power of 80% (at an α value of 0.05) for odds ratios > 2. The expected prevalence of Acute Respiratory Distress Syndrome (ARDS, 14–15% of infected subjects) and the choice of dichotomous groups (ARDS vs. a-/oligosymptomatic) were expected to provide an increase in the study power.
Study variables were age, biological sex, and severity of COVID-19 symptoms.

4.3. Molecular Methods

Blood samples were collected on Dried Blood Spots (DBS) from finger-sticking. We opted for this solution in order to simplify blood sample collection, storage and delivery from the recruiting centers.
gDNA extraction from DBS was performed as previously reported [27], using three 3.5 mm punches per sample and the DNA Blood Spot LH kit (Perkin Elmer, Shelton, CT, USA).
NGS barcoded libraries were prepared using the AmpliSeq approach. Briefly, we designed an AmpliSeq panel (https://www.ampliseq.com, accessed on 1 May 2020) for the HLA genes listed in Supplemental Table S1 along with their relative accession numbers (IAD_199394_182), including 819 amplicons. The coverage of the coding region of each gene was 100%, plus we included at least 50 bp of the intronic region flanking each exon.
For library preparation and purification, we used the Ion AmpliSeqTM Library Kit plus (ThermoFisher Scientific, Waltham, MA, USA); for barcoding, we used the IonCode Plates (ThermoFisher Scientific, Waltham, MA, USA), 96 barcoded libraries/run. For purification, we used the AMPure XP Beads (Beckman Coulter, Brea, CA, USA), according to the AmpliSeq protocol instructions. Library preparation was performed in a semi-automatized pipeline, using a home-developed script for the Sciclone NGS Instrument (Perkin Elmer, Shelton, CT, USA).
Following purification, for library concentration balance and pool preparation, we used the Ion Library Equalizer kit (Thermofisher, Waltham, MA, USA), according to the manufacturer’s protocol. The ION Chef Instrument (ThermoFisher Scientific, Waltham, MA, USA) was used for emulsion PCR preparation and chip loading, using the ION Ion PI Hi-Q Chef Kit (ThermoFisher Scientific, Waltham, MA, USA). Sequencing was performed by the ION Proton Instrument, according to the manufacturer’s protocols. Ninety-six samples were run per chip. Basecalling, alignment to GRCh37/hg19 genome, and generation of the variant caller file (VCF) were performed by the resident Torrent Suite version 5.6. Variants (HLA-SNPs) were annotated by the ION-Reporter suite (version 5.8, https://ionreporter.thermofisher.com, accessed on 1 January 2022). We find a total number of 5578 different SNPs (see Supporting Information File) in our cohort, using the default Ion-Reporter parameters.

4.4. NGS Metrics

Quality metrics for NGS data were the following: per base Q30 and average read length between 125 and 275 bp (based on panel characteristics). Only samples with average depth ≥ 30× and uniformity ≥ 80% were analyzed. Samples that did not meet these quality standards were re-run; multiple BAM files of the same sample were merged using Samtools (version 1.13) [28].
The average number of reads per sample was 401,781 ± 295,638, the average depth 4921 ± 2942×, the mean uniformity, 82.51% ± 26.47%.

4.5. Sanger Sequencing

Significant variants, located in the coding region or predicted to affect exon splicing, were confirmed by Sanger sequencing.
Briefly, we designed the specific primer pairs, flanking the specific SNP/SNPs of interest, reported in Supporting Information, Table S4. For all PCR reactions, we used the 2× PCR Taq® DNA Polymerase (Promega Co., Madison, WI, USA), 10 µM of each primer and 30 ng of gDNA, in a final volume of 12 µL. Thermal cycling conditions were as follows: 95 °C, 5′; 35× (95 °C, 45″; variable annealing temperature indicated in Supporting Information Table S4, 45″; 72 °C, 30″); 72 °C, 5′. PCR products were analyzed by agarose gel (1.5% in TAE 1×) standard electrophoresis, and stained with ethidium bromide.
PCR products were purified by exonuclease and alkaline phosphatase (ExoSAP-IT, ThermoFisher Scientific, Waltham, MA, USA) before sequencing, by incubation at 37 °C, 15′ and enzyme heat inactivation (80 °C, 15′). Sequencing reactions were performed by the BigDye Terminator v3.1 Cycle Sequencing Kit (ThermoFisher Scientific, Waltham, MA, USA), in a final volume of 10 µL, according to the manufacturer’s protocol. Sequencing reaction products were purified using the BigDye XTerminator Purification Kit (ThermoFisher Scientific, Waltham, MA, USA), according to the manufacturer’s protocol. Capillary electrophoresis was performed on the Applied Biosystems 3130 Genetic Analyzer, and analyzed by the sequencing analysis software (Version 7.1, ThermoFisher Scientific, Waltham, MA, USA). The samples were electrophoresed, and the resulting data were analyzed to obtain sequence reads (Supplementary Figure S1).

4.6. Bioinformatics and Statistical Analysis

For canonical HLA-genotyping, we used HLA-VBSeq and the IMGT/HLA database21 (release3.15.0) [9].
An in-house Python3.0 script was used to process in batch HLA-VBSeq input files, to write relative output files and to perform statistical analysis of canonical genotyping results.
The association analysis was based on generalized linear models and performed using R4.1.2 (https://cran.r-project.org, accessed on 1 April 2022), a framework and coding language widely used for data analyses. We specifically made use of SNPassoc (https://cran.r-project.org/web/packages/SNPassoc/index.html, accessed on 1 April 2022), an R package (i.e., a collection of functions and other R objects) comprising multiple tools for genetic association studies. For each SNP, single alleles were set as independent variables; we computed the OR associated to the mere presence of the allele, irrespective of the dosage (heterozygous or homozygous) on the basis of the allelic frequency, by chi squared testing on contingency tables. We then assessed for the possible incremental effect of zero, one or two copies of the SNP (i.e., wild-type homozygous, heterozygous, and variant homozygous) in a log-additive model.
Subsequently, we compiled a logistic regression model in IBM SPSS 21.0.0.0 (https://www.ibm.com/products/spss-statistics, accessed on 1 April 2022) setting the genotypes at each HLA gene, age and biological sex as co-variates, to evaluate any combined effect given by age and/or sex. We used a false-discovery-rate (FDR) control for multiple testing adjustment (FDR < 0.05) after calculating nominal p-values for each locus, as customary for genetic association studies. In brief, FDR is the expected proportion of “discoveries” (rejected null hypotheses) that are false (i.e., incorrect rejections of the null).
Haplotypes were predicted using statistically significant SNPs as inputs into Arlequin v3.0, designed for inferring gametic phase starting from NGS output data [24]. Association analysis between haplotypes and severity was performed by chi squared testing on haplotypic frequencies.
Sex-specific statistical analyses were performed by means of chi-squared testing on allelic frequencies contingency tables for each HLA allele, SNP and haplotype. In addition, chi-squared testing was also performed for each allele on contingency tables built by dividing the whole population by sex rather than phenotype, in order to assess possible differences in allele distribution between sexes. All table operations subsequent to the first R analysis were performed using Python3.0 (https://www.python.org/download/releases/3.0/, accessed on 1 April 2022) custom scripts, making use of the “CSV” (Comma-Separated Values, Python built-in), and “SciPy” (https://scipy.org/, accessed on 1 April 2022) libraries. These analyses include computation of OR and CI for HLA alleles, SNPs and haplotypes, in the whole population and in the sex-stratified ones.
A Python3.0 script was also created to compare HLA-SH and .fasta sequences of conventional HLA-alleles, available on the IMGT website (https://www.ebi.ac.uk/ipd/imgt/hla/, accessed on 1 September 2022), to identify correspondences between HLA-SH and HLA alleles.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms252313198/s1.

Author Contributions

F.D.T. and A.V. (Alessandro Vaisfeld) conceptualized and designed the study. G.G. collected the data, while M.V.F., A.N. and S.F. performed laboratory experiments. S.S. and P.D. performed data-analysis, developed bio-informatic tools and drafted the manuscript. A.C., L.L., A.V. (Alessandro Visentin), A.A., E.R., E.M.T., M.M., C.T., E.T., M.F. and R.C. were in charge of patient recruitment and clinical evaluation of patients, and revised the manuscript; S.C., J.M. and T.W.R. provided counseling and critically revised data and manuscript. F.D.T. supervised the project, is responsible for data, and wrote the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The present study was funded by Emergex Vaccines Holding Ltd. as an independent research grant.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committees of the three institutions involved (Catholic University of Sacred Heart/Policlinico “A. Gemelli” Foundation, University of Verona and University “Magna Graecia”, ID3490, approved on 9 October 2020). Clinicaltrials.gov registration number: NCT05711836.

Informed Consent Statement

For all subjects, written informed consent was collected.

Data Availability Statement

Data available in the Supporting Data Values file.

Acknowledgments

We are very grateful to patients and families for their participation in the study.

Conflicts of Interest

Silvano Coletti and Jan Mous were employed by the company Chelonia SA. 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.

References

  1. Zhu, Z.; Lian, X.; Su, X.; Wu, W.; Marraro, G.A.; Zeng, Y. From SARS and MERS to COVID-19: A brief summary and comparison of severe acute respiratory infections caused by three highly pathogenic human coronaviruses. Respir. Res. 2020, 21, 224. [Google Scholar] [CrossRef] [PubMed]
  2. Williamson, E.J.; Walker, A.J.; Bhaskaran, K.; Bacon, S.; Bates, C.; Morton, C.E.; Curtis, H.J.; Mehrkar, A.; Evans, D.; Inglesby, P.; et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature 2020, 584, 430–436. [Google Scholar] [CrossRef] [PubMed]
  3. Fricke-Galindo, I.; Falfán-Valencia, R. Genetics Insight for COVID-19 Susceptibility and Severity: A Review. Front. Immunol. 2021, 12, 622176. [Google Scholar] [CrossRef] [PubMed]
  4. Pairo-Castineira, E.; Rawlik, K.; Bretherick, A.D.; Qi, T.; Wu, Y.; Nassiri, I.; McConkey, G.A.; Zechner, M.; Klaric, L.; Griffiths, F.; et al. GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19. Nature 2023, 619, E61. [Google Scholar] [CrossRef] [PubMed]
  5. Langton, D.J.; Bourke, S.C.; Lie, B.A.; Reiff, G.; Natu, S.; Darlay, R.; Burn, J.; Echevarria, C. The influence of HLA genotype on the severity of COVID-19 infection. HLA 2021, 98, 14–22. [Google Scholar] [CrossRef]
  6. Hoseinnezhad, T.; Soltani, N.; Ziarati, S.; Behboudi, E.; Mousavi, M.J. The Role of HLA Genetic Variants in COVID-19 Susceptibility, Severity, and Mortality: A Global Review. J. Clin. Lab. Anal. 2024, 38, e25005. [Google Scholar] [CrossRef]
  7. January, W., 3rd; Phillip, S.; Manuel, H.; Charlotte, T.; Müller, M.; Dimitri, P.; Claudia, Q.; Krüger, U.; Elisa, T.H.; Lena, L.; et al. HLA-C* 04:01 is a Genetic Risk Allele for Severe Course of COVID-19. MedRxiv 2020. [CrossRef]
  8. Correale, P.; Mutti, L.; Pentimalli, F.; Baglio, G.; Saladino, R.E.; Sileri, P.; Giordano, A. HLA-B*44 and C*01 Prevalence Correlates with Covid19 Spreading across Italy. Int. J. Mol. Sci. 2020, 21, 5205. [Google Scholar] [CrossRef]
  9. Barker, D.J.; Maccari, G.; Georgiou, X.; A Cooper, M.; Flicek, P.; Robinson, J.; E Marsh, S.G. The IPD-IMGT/HLA Database. Nucleic Acids Res. 2023, 51, D1053–D1060. [Google Scholar] [CrossRef]
  10. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinform. 2007, 1, 117693430500100. [Google Scholar] [CrossRef]
  11. Markov, P.V.; Ghafari, M.; Beer, M.; Lythgoe, K.; Simmonds, P.; Stilianakis, N.I.; Katzourakis, A. The evolution of SARS-CoV-2. Nat. Rev. Microbiol. 2023, 21, 361–379. [Google Scholar] [CrossRef] [PubMed]
  12. Butler-Laporte, G.; Povysil, G.; Kosmicki, J.A.; Cirulli, E.T.; Drivas, T.; Furini, S.; Saad, C.; Schmidt, A.; Olszewski, P.; Korotko, U.; et al. Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative. PLoS Genet. 2022, 18, e1010367. [Google Scholar] [CrossRef] [PubMed]
  13. Park, G.M.; Lee, S.; Park, B.; Kim, E.; Shin, J.; Cho, K.; Ahn, K. Soluble HLA-G generated by proteolytic shedding inhibits NK-mediated cell lysis. Biochem. Biophys. Res. Commun. 2004, 313, 606–611. [Google Scholar] [CrossRef] [PubMed]
  14. Shimizu, M.; Kuroda, Y.; Uchida, M.; Takada, S.; Kamada, H.; Takahashi, D.; Nakajima, F.; Miyata, S.; Igarashi, S.; Satake, M. A new HLA-C allele with an alternative splice site in exon 3: HLA-C*03:23N. HLA 2020, 95, 555–560. [Google Scholar] [CrossRef] [PubMed]
  15. Cole, C.; Byrne, A.; Adams, M.; Volden, R.; Vollmers, C. Complete characterization of the human immune cell transcriptome using accurate full-length cDNA sequencing. Genome Res. 2020, 30, 589–601. [Google Scholar] [CrossRef] [PubMed]
  16. Yamamoto, F.; Suzuki, S.; Mizutani, A.; Shigenari, A.; Ito, S.; Kametani, Y.; Kato, S.; Fernandez-Viña, M.; Murata, M.; Morishima, S.; et al. Capturing Differential Allele-Level Expression and Genotypes of All Classical HLA Loci and Haplotypes by a New Capture RNA-Seq Method. Front. Immunol. 2020, 11, 941. [Google Scholar] [CrossRef]
  17. Johansson, T.; Koskela, S.; Yohannes, D.A.; Partanen, J.; Saavalainen, P. Targeted RNA-Based Oxford Nanopore Sequencing for Typing 12 Classical HLA Genes. Front. Genet. 2021, 12, 635601. [Google Scholar] [CrossRef]
  18. Klein, S.L.; Flanagan, K.L. Sex differences in immune responses. Nat. Rev. Immunol. 2016, 16, 626–638. [Google Scholar] [CrossRef]
  19. Lakshmikanth, T.; Consiglio, C.; Sardh, F.; Forlin, R.; Wang, J.; Tan, Z.; Barcenilla, H.; Rodriguez, L.; Sugrue, J.; Noori, P.; et al. Immune system adaptation during gender-affirming testosterone treatment. Nature 2024, 633, 155–164. [Google Scholar] [CrossRef]
  20. Lin, A.; Yan, W.-H. The Emerging Roles of Human Leukocyte Antigen-F in Immune Modulation and Viral Infection. Front. Immunol. 2019, 10, 964. [Google Scholar] [CrossRef]
  21. Nagamatsu, T.; Fujii, T.; Matsumoto, J.; Yamashita, T.; Kōzuma, S.; Taketani, Y. Human Leukocyte Antigen F Protein is Expressed in the Extra-Villous Trophoblasts but Not on the Cell Surface of them. Am. J. Reprod. Immunol. 2006, 56, 172–177. [Google Scholar] [CrossRef] [PubMed]
  22. Taubenberger, J.K.; Morens, D.M. 1918 Influenza: The mother of all pandemics. Emerg. Infect. Dis. 2006, 12, 15–22. [Google Scholar] [CrossRef] [PubMed]
  23. Timpson, N.J.; Greenwood, C.M.T.; Soranzo, N.; Lawson, D.J.; Richards, J.B. Genetic architecture: The shape of the genetic contribution to human traits and disease. Nat. Rev. Genet. 2018, 19, 110–124. [Google Scholar] [CrossRef] [PubMed]
  24. Moltke, I.; Grarup, N.; Jørgensen, M.E.; Bjerregaard, P.; Treebak, J.T.; Fumagalli, M.; Korneliussen, T.S.; Andersen, M.A.; Nielsen, T.S.; Krarup, N.T.; et al. A common Greenlandic TBC1D4 variant confers muscle insulin resistance and type 2 diabetes. Nature 2014, 512, 190–193. [Google Scholar] [CrossRef] [PubMed]
  25. Purcell, S.; Cherny, S.S.; Sham, P.C. Genetic Power Calculator: Design of linkage and association genetic mapping studies of complex traits. Bioinformatics 2003, 19, 149–150. [Google Scholar] [CrossRef]
  26. Hong, E.P.; Park, J.W. Sample size and statistical power calculation in genetic association studies. Genom. Inform. 2012, 10, 117–122. [Google Scholar] [CrossRef]
  27. Abiusi, E.; Vaisfeld, A.; Fiori, S.; Novelli, A.; Spartano, S.; Faggiano, M.V.; Giovanniello, T.; Angeloni, A.; Vento, G.; Santoloci, R.; et al. Experience of a 2-year spinal muscular atrophy NBS pilot study in Italy: Towards specific guidelines and standard operating procedures for the molecular diagnosis. J. Med. Genet. 2023, 60, 697–705. [Google Scholar] [CrossRef]
  28. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
Figure 1. (A): Size of the coding region of each gene (bp) vs. # of variants. (B): Ratio between # of variants and amount of sequenced gDNA.
Figure 1. (A): Size of the coding region of each gene (bp) vs. # of variants. (B): Ratio between # of variants and amount of sequenced gDNA.
Ijms 25 13198 g001
Table 1. Inclusion criteria for cases/controls and sample distribution.
Table 1. Inclusion criteria for cases/controls and sample distribution.
Inclusion CriteriaM
# (%)
F
# (%)
Total
# (%)
Median Age
Years (±IQR)
CasesOxygen saturation:
<94% on room air
<90% with known chronic hypoxic conditions
or
Respiratory rate > 24/min
333 (67.55)160 (32.45)493 (73)59.24 (11.45)
ControlsFew or no signs/symptoms associated with COVID-19
or
No relevant impairment in quality of life
86 (42.47)100 (57.53)186 (27)52.55 (16.73)
M: males; F: females; IQR: interquartile range.
Table 2. HLA alleles significantly associated with COVID-19 severity. The two cohorts were also stratified according to the biological sex.
Table 2. HLA alleles significantly associated with COVID-19 severity. The two cohorts were also stratified according to the biological sex.
GeneAllele #Whole PopulationMalesFemalesM vs. F
# of AllelesORCIp# of AllelesORCIp# of AllelesORCIpp
Cases
(n = 970)
Ctrls
(n = 284)
LowUpCases
(n = 654)
Ctrls
(n = 172)
LowUpCases
(n = 316)
Ctrls
(n = 198)
LowUp
A*24151421.4412.070.05104151.983.51.120.0247271.111.840.660.81.00
B*1435260.50.290.830.0119160.290.580.155*10-4161012.260.6610.48
B*5517140.450.220.930.031060.431.20.150.18780.541.510.190.350.25
C*0822170.470.250.90.0211120.230.530.15*10-41151.394.070.660.730.73
C*12120510.870.611.230.4392211.181.950.710.6128300.540.940.660.040.20
C*1515110.50.231.110.08660.260.80.080.03951.133.430.6610.10
DPA1*016362211.2811.640.054231081.091.540.770.712131131.562.240.660.020.75
DPA1*023181410.790.621.020.07219600.941.340.660.899810.660.950.660.030.64
DPB1*1723130.670.331.330.251831.595.480.460.635100.30.90.660.050.68
DPB1*416440.380.091.520.16130.090.830.010.04311.8918.280.660.970.5
DQB1*052291080.750.570.980.04148510.691.010.480.0781570.851.270.570.490.26
DRB1*1021111.910.2216.420.551180.350.890.140.041032.127.820.660.380.79
HLA alleles in our populations were represented as absolute number (#) and percentage (%). No statistically significant differences were seen in HLA allele distribution between sexes in the whole population, ruling out possible sampling bias (see the p-value in M vs. F column). Different groups were compared by χ2 test. The protective/risk effect of each HLA allele was evaluated by the odds ratio (OR) and lower (Low) and upper (Up) confidence interval (CI). Differences were considered statistically significant for p < 0.05. The significant HLA alleles, in at least one category, are indicated in bold in coloured cells.
Table 3. HLA-SNPs significantly associated with COVID-19 severity in chi squared testing.
Table 3. HLA-SNPs significantly associated with COVID-19 severity in chi squared testing.
GenegDNA (chr6:)dbSNP IDcDNAa.a. ChangeExonEffectORCIp
LowUp
Class I
HLA-B31324086G>Crs709054c.477C>Gp.Ala159=3R1.020.472.230.04
31324145A>Crs9266150c.418T>Gp.Tyr140Asp4R1.861.113.110.01
31324523T>Crs1131214c.285A>Gp.Ala95=2P0.440.240.820.01
31324526CTG>GGAn.a.c.280CAG>TCCp.Gln94Ser2P0.480.270.850.01
31324535GTA>CATn.ac.271TAC>ATGp.Tyr91Met2P0.50.290.880.02
31324539A>Trs1131202c.269T>Ap.Ile90Thr2P0.440.240.820.01
31324542T>Crs1131201c.266A>Gp.Gln89Arg2P0.420.230.770.01
HLA-C31237124T>Crs1130838c.1087T>Cp.Thr363Ser7R3.621.051.770.02
31237162C>Grs35708511c.1049C>Gp.Cys350Ser7R1.361.031.760.03
31237773T>Crs1130947c.985T>Cp.Thr329Ala5R1.851.442.37<10−4
31237774T>Grs41540512c.984T>Gp.Val328=5R1.851.442.37<10−4
31237779C>Trs146911342c.979C>Tp.Val327Leu5R1.851.521.10.02
31237858T>Crs34794906c.900T>Cp.Pro300=5R1.411.061.90.02
31237987C>Trs41556321c.895C>Tp.Glu299Ter4R1.481.052.090.02
31238009TT>CGrs1131014c.872TT>CGp.Gln291Pro4R1.311.031.660.02
31238027C>Ars41540117c.855C>Ap.Met285Ile4R1.471.062.040.01
31238029T>Crs2308622c.853T>Cp.Met285Leu4R1.281.011.620.04
31238138C>Trs1050320c.744C>Tp.Gln248His4R1.341.011.770.03
31239449C>Grs28626310c.270C>Gp.Lys90Asn2P0.620.470.82<10−4
31239501G>Trs1050409c.218G>Tp.Ala73Glu2R1.341.001.810.05
HLA-F29694777G>Ars17875384c.1153G>Ap.Arg385Gln7R3.621.2810.250.01
Class II
HLA-DQA132609806C>Trs707952c.388C>Tp.Thr130Ile3P0.720.520.990.04
32609952T>Crs707949c.534T>Cp.Phe179Leu3P0.690.490.960.02
32609969T>Grs707963c.551T>Gp.Asp184Glu3P0.680.490.950.02
32609974T>Grs707962c.556T>Gp.Ile186Ser3P0.660.470.920.01
HLA-DRB132551939GG>ATrs16822512c.317CC>ATp. Thr106Asn2P0.770.610.970.02
32552092A>Trs1059569c.174T>Ap.Phe55Tyr2R1.521.521.080.03
HLA-DRB532487256G>Ars112401921c.543C>Tp.Asp181=3P0.480.250.940.03
OR, lower (Low) and upper (Up) CI and p-values are indicated. Bolded SNPs in coloured cells were significant also in a multinomial logistic regression model, including age and biological sex as co-variates. gDNA: precise genomic position according to the hg19 version of the reference genome assembly (https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_000001405.13/ accessed on 1 July 2022); dbSNP ID: SNP reference number as of dbSNP database (https://www.ncbi.nlm.nih.gov/snp/ accessed on 1st July 2022); R: risk allele; P: protective allele
Table 4. Results of HLA-SNP multivariate analysis. The two cohorts were also stratified according to the biological sex.
Table 4. Results of HLA-SNP multivariate analysis. The two cohorts were also stratified according to the biological sex.
GenegDNA Coordinates
(dbSNP rsID)
Log Additive ModelGenotypic FrequenciesAllelic FrequenciesM vs. F
GlobalMalesFemalesGlobalMalesFemales
ORLowUppORLowUppORLowUppORLowUpPORLowUppORLowUppORLowerUpperpp
HLA-B31324086G>C
(rs709054)
10.52.230.051.090.781.520.671.110.691.770.141.150.71.880.611.120.871.460.391.630.761.530.711.17.411.630.680.23
31324145A>C
(rs9266150)
1.91.13.110.011.813.10.042.214.80.051.50.593.820.4821.63.410.011.840.933.670.101.680.684.130.290.07
31324523T>C
(rs1131214)
0.40.20.820.010.50.20.90.020.620.261.470.220.20.10.730.010.40.20.80.010.560.251.260.160.20.10.750.010.73
31324526CTG>GGA
(n.a.)
0.50.30.850.020.790.511.230.290.670.31.520.360.30.10.70.010.50.30.850.020.610.291.320.190.30.10.80.010.88
31324535GTA>CAT
(n.a)
0.50.30.880.030.50.30.90.030.710.321.580.380.20.10.780.010.50.30.880.020.630.31.370.280.30.10.80.010.8
31324539A>T
(rs1131202)
0.40.20.820.010.50.20.90.020.690.261.470.320.20.70.730.010.40.20.810.010.560.251.260.160.20.10.750.010.73
31324542T>C
(rs1131201)
0.40.20.770.010.40.20.80.010.540.241.250.150.20.70.730.010.40.20.7710-30.50.231.10.100.20.10.750.010.65
HLA-C31237124T>C
(rs1130838)
1.41.11.770.022.91.65.310−32.91.36.30.012.230.826.050.121.41.11.790.021.61.22.36*10-31.10.751.630.690.43
31237162C>G
(rs35708511)
1.311.760.032.91.55.13*10-32.81.16.80.031.980.715.50.191.311.730.031.61.12.37*10-31.010.691.4910.17
31237773T>C
(rs1130947)
1.91.42.37<10−41.280.861.920.241.5413.020.060.960.541.1711.311.30.031.61.12.20.010.970.681.380.850.11
31237774T>G
(rs41540512)
1.91.42.37<10−41.280.861.920.241.5413.020.060.960.541.1711.311.30.031.61.12.20.010.970.681.380.850.11
31237779C>T
(rs146911342)
1.51.12.10.011.51.12.20.022.21.33.84*10-31.320.772.280.341.51.12.10.0121.33.34*10-31.140.721.820.630.56
31237858T>C
(rs34794906)
1.41.11.90.021.71.22.35*10-31.450.92.340.121.660.992.80.061.41.11.910.011.512.30.051.310.852.020.230.28
31237987C>T
(rs41556321)
1.51.12.090.021.81.22.63*10-31.91.13.30.011.420.862.460.271.51.12.050.021.71.12.80.031.210.741.970.460.47
31238009TT>CG
(rs1131014)
1.311.660.021.71.22.42*10-31.560.982.470.071.480.92.440.131.41.11.80.011.51.12.20.021.160.791.70.490.06
31238027C>A
(rs41540117)
1.51.12.040.021.81.22.63*10-31.81.13.10.031.540.872.670.161.51.12.120.011.81.12.90.021.270.782.060.390.49
31238029T>C
(rs2308622)
1.311.620.041.61.12.40.021.7212.970.071.070.61.930.881.311.660.021.61.12.20.011.050.741.50.780.36
31238138C>T
(rs1050320)
1.311.770.041.80.12.52*10-32.21.43.62*10-31.060.631.780.891.220.961.560.111.411.970.050.980.691.4210.18
31239449C>G
(rs28626310)
0.60.50.82<10−40.50.40.7<10−40.50.30.73*10-30.80.481.320.430.10.40.810-30.50.30.7<10−40.870.571.320.520.39
HLA-F29694777G>A
(rs17875384)
3.61.310.30.011.880.814.330.142.810.6412.280.185.11.1230.023.61.310.20.012.90.6712.530.194.91.121.90.020.62
HLA-DQA132609806C>T
(rs707952)
0.70.50.990.050.660.41.080.100.730.361.490.440.640.31.380.320.60.49.430.020.670.381.20.200.590.321.10.100.47
32609952T>C
(rs707949)
0.70.50.960.030.740.531.040.090.70.441.110.140.670.41.130.140.760.571.020.070.760.511.140.200.710.431.120.150.58
32609969T>G
(rs707963)
0.70.50.950.020.740.521.030.080.740.461.180.220.620.371.040.080.760.561.010.060.80.531.190.280.670.431.060.100.71
32609974T>G
(rs707962)
0.70.50.920.010.710.50.990.050.720.451.150.180.60.40.980.050.70.60.980.040.780.521.170.230.640.411.010.060.66
HLA-DRB132551939GG>AT
(rs16822512)
0.80.60.970.030.830.591.150.270.60.410.050.890.541.470.700.80.60.970.030.720.511.020.660.750.511.110.160.59
32552092A>T
(rs1059569)
1.51.13.070.051.81.22.90.012.21.14.40.031.170.62.30.731.612.460.041.910.993.680.051.30.692.450.430.62
HAL-DRB532487256G>A
(rs112401921)
0.50.30.940.030.40.20.90.050.20.60.16*10-30.840.223.2110.40.20.870.020.30.10.64*10-30.950.33.0310.40
The two trans (i.e., in each copy of the single HLA genes) alleles were included as covariates in a log-additive model. Significant HLA-SNPs were analyzed by genotypic (chi squared testing from 2 × 3 contingency table) and allelic frequencies (2 × 2 table chi squared testing), in the whole population or stratified by biological sex. As for conventional alleles (Table 2), no statistically significant differences were seen in SNP distribution between sexes in the whole population, ruling out possible sampling bias (see the p-value in M vs. F column). The table reported odds ratio (OR), lower (Low) and upper (Up) confidence intervals (CI) and p-values for each analysis. Bolded SNPs in coloured cells were significant.
Table 5. HLA super haplotypes associated with COVID-19 severity in the whole population: males and females.
Table 5. HLA super haplotypes associated with COVID-19 severity in the whole population: males and females.
HLA GeneHaplotype ID #Whole PopulationMalesFemalesM vs. F
# of AllelesORCIp# of AllelesORCIp# of AllelesORCIpp
Cases n = 964Ctrls n = 396LowUpCases n = 654Ctrls n = 186LowUpCases n = 310Ctrls n = 210LowUp
-B36332700.890.71.150.414111440.490.340.723*10-42221261.681.162.440.010.75
5184751.010.751.361121480.650.440.960.0463271.731.062.820.040.2
92319.651.371.740.011113.160.4124.670.42120---0.010.23
1553112.041.053.940.043552.050.795.30.191862.10.825.370.170.9
-C4114730.590.430.822*10-380400.510.330.772*10-334330.660.391.110.150.47
866231.190.731.950.5642170.680.381.230.262462.851.157.110.030.36
1773231.330.822.160.3046150.860.471.580.752782.411.075.410.040.71
19240.20.041.120.11140.070.010.630.0110---10.28
44140.10.010.910.04110.280.024.550.9203---0.130.32
-F13443.581.2610.160.012022.90.6712.530.221424.921.1121.870.040.74
-DQA117842991.411.071.870.025031580.590.380.920.022811414.742.947.654*10-110.27
2111560.790.561.120.2174350.550.350.850.0137211.220.692.150.590.32
-DRB116652571.21.540.940.144301360.711.010.490.062351212.33.361.5810-40.68
222180.490.260.920.0415110.370.170.830.02770.670.231.940.640.25
31931120.630.480.8310-3117620.440.30.6310-576501.040.691.570.940.21
-DRB5113130.40.180.880.03680.210.070.610-3750.950.33.0310.53
No statistically significant differences were seen in HLA haplotype distribution between sexes in the whole population, ruling out possible sampling bias (see the p-value in M vs. F column). The table reports odds ratio (OR), lower (Low) and upper (Up) confidence intervals (CI) and p-values for each analysis. Bolded SNPs in coloured cells were significant.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Spartano, S.; Faggiano, M.V.; Guidi, G.; D’Ambrosio, P.; Vaisfeld, A.; Novelli, A.; Falqui, S.; Cingolani, A.; Lambertenghi, L.; Visentin, A.; et al. Sex-Specific HLA Alleles Contribute to the Modulation of COVID-19 Severity. Int. J. Mol. Sci. 2024, 25, 13198. https://doi.org/10.3390/ijms252313198

AMA Style

Spartano S, Faggiano MV, Guidi G, D’Ambrosio P, Vaisfeld A, Novelli A, Falqui S, Cingolani A, Lambertenghi L, Visentin A, et al. Sex-Specific HLA Alleles Contribute to the Modulation of COVID-19 Severity. International Journal of Molecular Sciences. 2024; 25(23):13198. https://doi.org/10.3390/ijms252313198

Chicago/Turabian Style

Spartano, Serena, Maria Vittoria Faggiano, Giovanna Guidi, Pino D’Ambrosio, Alessandro Vaisfeld, Agnese Novelli, Salvatore Falqui, Antonella Cingolani, Lorenza Lambertenghi, Alessandro Visentin, and et al. 2024. "Sex-Specific HLA Alleles Contribute to the Modulation of COVID-19 Severity" International Journal of Molecular Sciences 25, no. 23: 13198. https://doi.org/10.3390/ijms252313198

APA Style

Spartano, S., Faggiano, M. V., Guidi, G., D’Ambrosio, P., Vaisfeld, A., Novelli, A., Falqui, S., Cingolani, A., Lambertenghi, L., Visentin, A., Azzini, A., Righi, E., Trecarichi, E. M., Mazzitelli, M., Coletti, S., Mous, J., Rademacher, T. W., Torti, C., Tacconelli, E., ... Tiziano, F. D. (2024). Sex-Specific HLA Alleles Contribute to the Modulation of COVID-19 Severity. International Journal of Molecular Sciences, 25(23), 13198. https://doi.org/10.3390/ijms252313198

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

Article Metrics

Back to TopTop