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.