HLA and red blood cell antigen genotyping in SARS-CoV-2 convalescent plasma donors
==================================================================================

* William Lemieux
* Josée Perreault
* Gabriel André Leiva-Torres
* Nadia Baillargeon
* Jessica Constanzo Yanez
* Marie-Claire Chevrier
* Lucie Richard
* Antoine Lewin
* Patrick Trépanier

## Abstract

**Introduction** The SARS-CoV-2 pandemic has put significant additional pressure on healthcare systems throughout the world. The identification of at-risk population beyond age, pre-existing medical conditions and socioeconomic status has been the subject of only a small part of the global COVID-19 research so far. To this day, more data is required regarding the association between HLA allele and red blood cell (RBC) antigens’ expression in regard to SARS-CoV-2 infection susceptibility and virus clearance capability, and COVID-19 susceptibility, severity, and duration.

**Methods** The phenotypes for ABO and RhD, and the genotypes for 37 RBC antigens and HLA-A, B, C, DRB1, DQB1 and DPB1 were determined using high throughput platforms (Luminex and Next-generation Sequencing) in 90 Caucasian convalescent plasma donors. The results were compared to expected reference frequencies, local and international databases, and literature.

**Results** The AB group was significantly increased (1.5x, p=0.018) and a non-significant (2.2x, p=0.030) increase was observed for the FY*A allele frequency in the convalescent cohort (N=90) compared to reference frequencies. Some HLA alleles were found significantly overrepresented (HLA-B*44:02, C*05:01, DPB1*04:01, DRB1*04:01 and DRB1*07:01) or underrepresented (A*01:01, B51:01 and DPB1*04:02) in convalescent individuals compared to the local bone marrow registry population.

**Conclusion** Our study of infection-susceptible but non-hospitalized Caucasian COVID-19 patients contributes to the global understanding of host genetic factors associated with SARS-CoV-2 infection susceptibility and severity of the associated disease.

Key words
*   Covid-19
*   SARS-CoV-2
*   Susceptibility
*   Red Blood Cell Antigens
*   Human Leukocyte Antigens

## Introduction

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has caused over 5.3 million deaths worldwide as of December 2021 ([https://covid19.who.int/](https://covid19.who.int/)). Intense work has been done since the beginning of the pandemic to protect the most at-risk populations. Being male, of older age, obese, of a certain ethnic origin, having diabetes, asthma, or many other medical conditions, have all been associated with an increased risk of COVID-19-related complications or death1–5. Genetic factors, such as the expression of angiotensin-converting enzyme (ACE)-related genes, may also play a role in disease severity and could serve as a predictive marker for at-risk populations6. Efforts in getting a better understanding of the molecular gateways for viral entry and interactions with the host’s response have been relentless, focusing on ACE2 expression, HLA, cytokine storms and TLRs7. Such research could lead to an improved understanding of the disease’s susceptibility and severity, of the capacity for virus clearance, of long term symptoms8 and mortality consequent to infection in humans.

Relationships between ABO blood group expression and susceptibility to infection by SARS-CoV-2, as well as the risks of hospitalization and death from COVID-19 have been explored by a few groups and recently reviewed9,10. O group individuals were identified as having a decreased risk of infection compared to other ABO groups, although no differences were observed regarding hospitalization and death rates associated with COVID-1911,12. Several hypotheses to explain this observed link between ABO type and risk of infection have been suggested, such as the presence of anti-A and -B antibodies in O individuals13 and the binding of Receptor-Binding Domain (RBD) of the viral Spike protein to group A antigens14. Given the important immunological deregulation and potential cytokine storm associated with mortality in COVID-19 patients15 and the interesting report from an Italian laboratory of higher rates of direct antiglobulin test (DAT) reactivity in COVID-19 patients16, RBC antigens and related antibodies could be involved in COVID-19 and its clinical manifestations.

Of additional importance is the potential contribution of HLAs, given their central role in regulating immunity against viruses17. Considerable efforts have been made in trying to predict and identify protective or susceptibility-enhancing HLA alleles18,19. Binding assays have shown that HLA-A*02:01 and B*40:01 could preferentially associate with SARS-CoV-2 epitopes20. Using bioinformatic prediction, HLA-A*02:03 and A*31:01 were identified as protective, while A*03:02 was identified as a risk allele21. *In silico* binding affinity studies have shown that HLA-B*46:01 could increase susceptibility to disease, whereas HLA-B*15:03 is associated with protective immunity, and HLA-A*02:02, B*15:03, C*12:03 were the most frequently encountered haplotypes associated with presentation of viral epitopes18. Other retrospective studies aimed at identifying potentially protective or risk-increasing alleles22–28 have been published. Such a global effort requires several geographically diverse laboratories to analyze and share available data in order to generate a comprehensive overview.

In the present study, we hypothesized that associations might exist between certain HLA alleles or blood groups and SARS-CoV-2 susceptibility, as well as the ability to overcome infection without hospitalization, in convalescent plasma donors. We therefore sought to analyze and determine the existence of any trends in the expression of an extended panel of RBC antigens (ABO, RhD and 37 other antigens), and of HLA-A, B, C, DRB1, DQB1, and DPB1 alleles within a Caucasian convalescent plasma donor cohort, in comparison to different reference frequencies (textbook, local and international databases, and literature). The identification of differential patterns of RBC antigen or HLA expression in convalescent individuals, who were infected but were not hospitalized, could contribute to a better understanding of SARS-CoV-2 susceptibility and COVID-19 severity, and may also help interpret the discrepancies observed in the rates of infection in different ethnic groups5,29.

## Material and methods

### Samples

Caucasian convalescent plasmas were randomly chosen from adult participants of the Québec cohort in the CONCOR-1 clinical trial (#[NCT04348656](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT04348656&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom)). All subjects received an official diagnosis of COVID-19 by the Québec Provincial Health Authority after epidemiologic investigation or after confirmation by polymerase chain reaction (PCR) test. All subjects were COVID-19 symptomatic during infection, cleared the infection without hospitalization, and were free of symptoms for at least two weeks before donation. Since COVID-19 diagnostic was given before spring 2021, subjects were presumably infected with the Alpha (B.1.1.7) and Beta (B.1.351) variant, since no case of the Delta (B.1.617.2) variant was yet reported in the province of Quebec30. All subjects were self-identified Caucasians. Our convalescent plasma donor cohort had an average age of 40.4 ± 15.0 years and consisted of 68% males. Donors were not selected with respect to their ABO group. All donors gave consent to participate in this research project, which was approved by the Héma-Québec Research Ethics Committee.

### Phenotyping and genotyping

ABO and RhD phenotype testing were done by serologic detection using the PK7300 from Beckman Coulter, as per the manufacturer’s protocol. The DNA used for genotyping was extracted from the buffy coat of whole blood samples collected in ethylenediaminetetraacetic acid (EDTA) tubes, using QIAamp Blood Mini kit (Qiagen, Hilden, Germany). RBC genotyping was done using the Luminex xMAP® technology with the ID CORE XT platform (Progenika Biopharma-Grifols, Bizkaia, Spain), as per the manufacturer’s protocol, for the following blood group antigens: Rh (C, c, E, e, Cw, hrS, hrB, V, VS), Kell (K, k, Kpa, Kpb, Jsa, Jsb), Kidd (Jka, Jkb), Duffy (Fya, Fyb), MNS (M, N, S, s, U, Mia), Diego (Dia, Dib), Dombrock (Doa, Dob, Hy, Joa), Colton (Coa, Cob), Yt (Yta, Ytb) and Lutheran (Lua, Lub). HLA genotyping was done by next-generation sequencing (NGS) on a MiSeqDx (Illumina, San Diego, CA, USA), using NGSgo®-Ampx v2 kits, and interpreted with NGSengine® v2.21 (both from GenDX, Utrecht, The Netherlands).

### Statistical analyses

For RBC antigen and HLA allele comparisons in Tables 1, 2 and 3, population-wide proportions were assumed to correspond to the Caucasian prevalence estimates (from the Blood Group Antigen FactsBook31 and from the National Marrow Donor Program database32), while the HLA G-group allele frequencies of the 90 participants (180 individual HLA alleles) were estimated along with the Clopper-Pearson 95% confidence interval. Z-tests for two proportions were used to test for statistical significance between populational and observed antigen prevalence. A p value inferior to the Bonferroni correction for multiple comparison per antigen group was considered significant.

View this table:
[Table 1](http://medrxiv.org/content/early/2022/01/21/2021.03.17.21253821/T1)

Table 1 
RBC antigen frequencies for the Caucasian convalescent plasma donors compared to reference frequencies

View this table:
[Table 2](http://medrxiv.org/content/early/2022/01/21/2021.03.17.21253821/T2)

Table 2 
ABO and RhD blood group distributions for the Caucasian convalescent plasma donors compared to reference frequencies

View this table:
[Table 3](http://medrxiv.org/content/early/2022/01/21/2021.03.17.21253821/T3)

Table 3 
Caucasian convalescent donor HLA allele frequency comparison with the NMDP database for HLA*A, B, C, DRB1 and DQB1. The g notation specifies G-groups corresponding to those presented in Gragert, Loren, Abeer Madbouly, John Freeman, and Martin Maiers. Six-locus high resolution HLA haplotype frequencies derived from mixed-resolution DNA typing for the entire US donor registry. *Human Immunology* 74, n° 10.

For HLA genotype comparison in Tables 4 and 5, individuals from Héma-Québec’s bone marrow donor registry were used as reference. In total, 1370 registered individuals typed in high-resolution, from the Montréal and Montérégie regions and of self-reported Caucasian ethnicity, were selected to match with the characteristics of the studied convalescent cohort. Allele frequencies were calculated using the GENE[RATE] tool for HLA-A, B, C, DRB1, DQB1, and DPB133. HLA allelic frequencies were used with a modification of the hierfstat package to calculate the genetic distance (Latter) globally and for each locus between the subjects and the reference population34,35. Standardized residuals of the cohort subjects were calculated to identify alleles with significant differences between the subjects and the registry. Residuals were calculated by considering the subjects’ allele frequencies as the independent variable and the registry frequency as the dependant variable for each locus independently. Frequencies were deemed different at or above a difference of 3 absolute standardized residual. Statistical analyses were performed using the free statistical software R36.

View this table:
[Table 4](http://medrxiv.org/content/early/2022/01/21/2021.03.17.21253821/T4)

Table 4 
Genetic distance comparison between Caucasian convalescent donor HLA allele frequencies (2N=180) and Héma-Québec bone marrow donor registry (2N=2740) for HLA*A, B, C, DPB1, DQB1 and DRB1, per loci, and globally.

View this table:
[Table 5](http://medrxiv.org/content/early/2022/01/21/2021.03.17.21253821/T5)

Table 5 
HLA alleles identified as significantly over- or underrepresented from pairwise comparisons between the convalescent cohort (2N=180) and Héma-Québec bone marrow donor registry (2N=2740), and their potential clinical significance.

## Results

### Red blood cell genotype frequencies

Genotype frequencies for Rh, Kell, MNS, Duffy, Kidd, Diego, Dombrock, Colton, Yt and Lutheran blood groups were determined in each individual, and the resulting predicted phenotypes were compared to the expected Caucasian reference frequencies. The FY\*A/\*A genotype (Fy(a+b-) predicted phenotype) appears to be overrepresented (non-significant, p=0.030) in our convalescent cohort compared to expected frequencies (0.256 vs 0.170, respectively), as presented in Table 1. Incidentally, FY\*A/\*B (Fy(a+b+) predicted phenotype) individuals appear to be trending towards a decreased frequency of 0.400 compared to the expected 0.490 (Table 1), although the trend is not significant (p=0.087). Overall, no antigen group combinations deviated significantly from the expected frequencies.

### ABO and RhD phenotyping

ABO and RhD phenotypes were determined for each convalescent individual and compared to the expected Caucasian reference frequencies for the entire cohort (N=90) and within FY*A/A individuals (N=23), given its non-significant but trending towards overrepresentation (Table 1). Table 2 presenting the ABO phenotyping analysis for the cohort allowed for the identification of a significant (p=0.0178) 2.2x increase for the AB group compared to reference frequencies (0.089 vs 0.040, respectively). While non-significant (p=0.110), the O group trends towards a 0.8x underrepresentation within convalescent individuals compared to expected frequencies (0.356 vs 0.440, respectively). Interestingly, a significant (p=0.028) 3.3x AB group overrepresentation was also found within FY*A individuals (Table 2, FY\*A/\*A). An apparent but non-significant decrease in O group individuals can be observed within the FY\*A/\*A individuals versus the reference frequency (0.304 vs 0.440). No other trend or significant observation was made regarding A and B groups, nor regarding RhD expression.

### HLA allele frequency comparisons

HLA typing was done by NGS in the Caucasian convalescent donors’ cohort and individual allele frequencies were determined. The convalescent cohort allele frequencies (2N=180) were first compared to the NMDP registry frequencies for HLA-A, B, C, DRB1 and DQB1 for the most frequent alleles (Table 3), and no significant differences were identified in the most common alleles in the cohort. The convalescent cohort frequencies for HLA-A, B, C, DRB1, DQB1 and DPB1 were then compared to the Héma-Québec Stem Cell Donor Registry frequencies for Caucasians in the same geographical region (2N=2740). The genetic distance between the convalescent donors’ cohort and the Héma-Québec Registry was calculated for all loci together and on a per-locus basis for HLA-A, B, C, DRB1, DQB1, and DPB1 (Table 4). There was no comparative measurement available for the genetic distance, however the distance is low in all loci compared. Standardized residual analysis was conducted (Table 5 and Figure 1) and led to the identification of alleles that were significantly different between the convalescent cohort (N=90 individuals) and the Héma-Québec Registry (N=1370 individuals). HLA-B*44:02, C*05:01, DPB1*04:01, DRB1*04:01 and DRB1*07:01 were significantly overrepresented within our cohort, while A*01:01, B51:01 and DPB1*04:02 were significantly underrepresented. Finally, the alleles identified as significantly different were used to search previously published data for suspected associations with the SARS-CoV-2 virus infection and COVID-19 disease characteristics; these results are presented in Table 5. Of note, for two of the eight alleles that we have identified (DPB1*04:02 and DRB1*07:01), there was no data in the literature for comparison.

![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/21/2021.03.17.21253821/F1.medium.gif)

[Figure 1](http://medrxiv.org/content/early/2022/01/21/2021.03.17.21253821/F1)

Figure 1 
HLA Allelic frequencies in the convalescent cohort and the Héma-Québec bone marrow donor registry. For each locus, the cohort allele frequencies (horizontal axis) are plotted against the registry frequencies (vertical axis). The dashed line represents identical expression in the cohort and registry. The color of the alleles represents the absolute standardized residuals (|Std. Res.|) to the identity line. The dotted lines are the thresholds of 3 |Std. Res.| away from the identity line.

## Discussion

Our study took a deeper look into the RBC and HLA characteristics of a COVID-19 recovered and non-hospitalized cohort enrolled in the CONCOR-1 convalescent plasma study. We used ABO and RhD automated blood donor testing, and RBC and HLA high throughput genotyping platforms to determine the existence of potential trends regarding the frequencies of ABO, RhD, 37 other RBC antigens and HLA genotypes within the cohort, compared to reference Caucasian populations from textbooks, public databases, and our local bone marrow donor registry. A significant AB blood group overrepresentation was identified, as well as a non-significant trend in FY*A/A individuals. These results suggest a possible involvement of ABO and Duffy red blood cell antigens in SARS-CoV-2 susceptibility and COVID-19 severity, as all these individuals contracted the virus, yet only had mild symptoms, and cleared the infection without hospitalization.

The case of HLA association with disease susceptibility and severity is more complex. Overall, the genetic distance calculated from HLA allele frequencies is low (below 0.01) and suggests the cohort is similar to the reference population chosen. When looking at allele-level frequencies, eight HLA alleles were identified as potential markers between the convalescent cohort and a stem cell registry from the same geographical region by the standardized residuals analysis. We have also considered using 2×2 contingency tables along with Chi Squared tests to evaluate individual allelic differences. While some differences were significant after multiple testing corrections (HLA-A*68:02, B*39:06, B*51:01, B*52:01, C*12:02, C*14:02, DPB1*15:01, and DRB1*13:03), they involved less frequent alleles in the Héma-Québec Registry that were mostly absent from the cohort. Of these alleles, only HLA-B*51:01 was observed in the convalescent cohort, but at a lower frequency than in the Registry. Interestingly, this allele was also significant in the standardized residual analysis. One of the limitations of Chi Squared tests is the need for sufficient observed values (>5), which is not respected in these cases and prevents us from drawing conclusions from this analysis. Further studies and more data will be required to draw conclusion from these hypothetical associations.

One of the limitations of this study concerns the lack of a more diverse stratification of disease severity and the limited sample size for COVID-19 affected individuals. Indeed, our study lacks blood group and HLA data from hospitalized, deceased and asymptomatic COVID-19 patients, especially given that there could be a significant link between ABO and severity37. Additionally, while the historical ABO frequencies of Quebec Caucasian blood donors (internal data) matched that of FactsBook, such unbiased information about the frequencies of other RBC antigens is not currently available, hence the use of FactsBook Caucasian reference frequencies. Our study also does not directly address the major RBC and HLA antigen frequency differences between ethnicities. While our convalescent plasma donation program reflected our donor pool38, the low number of non-Caucasian individuals was insufficient to conduct statistical analysis, which is unfortunate given the importance of understanding the disproportionate impact of COVID-19 on minorities5,29. Nonetheless, the identification of a potential overrepresentation of FY\*A/\*A within Caucasian COVID-19 convalescent individuals and the potential implication of the Duffy blood group could have an impact on future research directions. The trend towards overrepresentation of the Fy(a+b-) predicted phenotype among our COVID-19 convalescent cohort could be explained by the absence of the Fyb antigen, since no significant difference was observed when comparing Fy(a+b+) and Fy(a-b+). In individuals of African descent, the Fy(a-b-) phenotype is caused by a GATA box mutation upstream of the *FY* gene silencing Fyb expression in RBCs39. Given that 67% of African Americans (AA) are Duffy null40, and that Duffy null patients have an increased mortality rate as a consequence of acute lung injury4, some groups have already hypothesized a role for Duffy in COVID-19 AA individuals40. We therefore suggest that Duffy allele identification might be used to select individuals at-risk for COVID-19 complications for further research on associations between COVID-19 and RBC antigens. The involvement of ABO blood groups in COVID-19 has previously been described13,41. The mechanism underlying the association remains elusive, but could be related to circulating natural anti-A and anti-B antibodies, or a low-efficiency furin cleavage in O-group individuals9. A significant overrepresentation of the AB group, and the non-significant but trend towards underrepresentation of the O group within our cohort, appear to be in agreement with other groups’ suggestion that O individuals could be less susceptible to SARS-CoV-2 infection11,12, given that infection was confirmed in every individual in our study. Our sample size does not allow us to draw conclusions as to whether individuals from the AB group are more susceptible to infection and more efficient at clearing infection without hospitalization, or whether this bias is a consequence of the trend towards underrepresentation of O individuals who are less susceptible to infection. It would be interesting to extend our observations within larger cohort analyses that include severely affected patients.

The involvement of HLA alleles in the disease outcome of COVID-19 patients is progressively getting clearer42,43. Various HLA alleles have shown high binding affinity to SARS-CoV-2 peptides,18,44 and trends have been observed in many populations22,23. While our sample size is limited, eight alleles were identified as significantly different between the studied cohort and matched individuals from the Héma-Québec Stem Cell Donor Registry, for which we already had high resolution HLA information. Two of these eight HLA alleles have not been previously identified for their association with COVID-19 in the literature: DPB1*04:02 and DRB1*07:01. The underrepresentation of HLA-A*01:01 is in agreement with its suggested association with high risk and mortality in COVID-19 individuals45,46, and the overrepresentation of HLA-B*44:02 is compatible with an increased susceptibility to infection26. Interestingly, A*01:01, B*44:02 and B*51:01 were predicted as weak binder of SARS-CoV-2 peptides, and none of the other alleles we identified were found to be strong peptide binders47. Given that we found A*01:01 and B51:01 to be underrepresented and B*44:02 overrepresented in our cohort, it is difficult to establish a relationship between these data without considering complete haplotypes or other disease severity groups. The other associations are inconclusive, but should be considered in larger cohorts. Overall, more data is required from as many diverse populations as possible in order to develop a comprehensive view, and to have data readily available for emerging variants and infection waves.

Altogether, we provide additional information regarding the role of RBC antigens and HLA in SARS-CoV-2 susceptibility, and consequential COVID-19 susceptibility, severity, resolution and long-term clinical consequences. More research needs to be done to get a better understanding of potentially at-risk populations, and for the identification of molecular pathways of this virus and its variants.

## Data Availability

None.

## Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

## Acknowledgements

The authors are grateful to the convalescent plasma donors who participated in this study and the Héma-Québec team involved in convalescent donor recruitment, sample collection, as well as the technologists who performed the phenotyping and genotyping workup. The authors thank Jean-Francois Leblanc and Marie-Eve Rhéaume for revising the manuscript. The authors are also thankful for healthcare workers’ dedication worldwide.

P.T. and W.L. designed the research and wrote the manuscript. P.T., L.R. and M-C.C. supervised the research. J.P., L.R., W.L. and J.C.Y. acquired the data. P.T., J.P., L.R., W.L., N.B., A.L. and J.C.Y. analysed the data. P.T., L.R., W.L., G.A.L, N.B., M-C.C and J.C.Y. reviewed and edited the manuscript.

## Footnotes

*   The authors have disclosed no conflicts of interest. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. All donors gave consent to participate in this research project, which was approved by the Héma-Québec Research Ethics Committee.

*   Addition of HLA coverage and overall statistical analysis update.

*   Received March 17, 2021.
*   Revision received January 21, 2022.
*   Accepted January 21, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Williamson, E. J. et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature 584, 430–436 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2521-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

2.  2.Guo, W. et al. Diabetes is a risk factor for the progression and prognosis of COVID-19. Diabetes Metab. Res. Rev. 36, (2020).
    
    

3.  3.Lighter, J. et al. Obesity in Patients Younger Than 60 Years Is a Risk Factor for COVID-19 Hospital Admission. Clin. Infect. Dis. 71, 896–897 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/CID/CIAA415&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

4.  4.Kangelaris, K. N. et al. The Association Between a Darc Gene Polymorphism and Clinical Outcomes in African American Patients With Acute Lung Injury. Chest 141, 1160–1169 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1378/chest.11-1766&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22207676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

5.  5.Raisi-Estabragh, Z. et al. Greater risk of severe COVID-19 in Black, Asian and Minority Ethnic populations is not explained by cardiometabolic, socioeconomic or behavioural factors, or by 25(OH)-vitamin D status: study of 1326 cases from the UK Biobank. J. Public Health 42, 451–460 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/pubmed/fdaa095&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

6.  6.Yamamoto, N. et al. SARS-CoV-2 infections and COVID-19 mortalities strongly correlate with ACE1 I/D genotype. Gene 758, 144944 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.gene.2020.144944&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

7.  7.Debnath, M., Banerjee, M. & Berk, M. Genetic gateways to COVID-19 infection: Implications for risk, severity, and outcomes. FASEB J. 34, 8787–8795 (2020).
    
    

8.  8.Callard, F. & Perego, E. How and why patients made Long Covid. Soc. Sci. Med. 268, 113426 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.socscimed.2020.113426&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

9.  9.Pendu, J. L., Breiman, A., Rocher, J., Dion, M. & Ruvoën-Clouet, N. ABO Blood Types and COVID-19: Spurious, Anecdotal, or Truly Important Relationships? A Reasoned Review of Available Data. Viruses 13, 160 (2021).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

10. 10.Gutiérrez-Valencia, M. et al. ABO blood group and risk of COVID-19 infection and complications: A systematic review and meta-analysis. Transfusion (Paris) n/a,.
    
    

11. 11.Barnkob, M. B. et al. Reduced prevalence of SARS-CoV-2 infection in ABO blood group O. Blood Adv. 4, 4990–4993 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1182/bloodadvances.2020002657&link_type=DOI) 

12. 12.Wu, Y., Feng, Z., Li, P. & Yu, Q. Relationship between ABO blood group distribution and clinical characteristics in patients with COVID-19. Clin. Chim. Acta 509, 220–223 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cca.2020.06.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

13. 13.Goel, R. et al. ABO blood group and COVID-19: a review on behalf of the ISBT COVID-19 working group. Vox Sang. vox.13076 (2021) doi:10.1111/vox.13076.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/vox.13076&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33578447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

14. 14.Wu, S.-C. et al. The SARS-CoV-2 receptor-binding domain preferentially recognizes blood group A. Blood Adv. 5, 1305–1309 (2021).
    
    

15. 15.Hu, B., Huang, S. & Yin, L. The cytokine storm and COVID-19. J. Med. Virol. 93, 250–256 (2021).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

16. 16.Berzuini, A. et al. Red cell–bound antibodies and transfusion requirements in hospitalized patients with COVID-19. Blood 136, 766–768 (2020).
    
    

17. 17.Sanchez-Mazas, A. HLA studies in the context of coronavirus outbreaks. Swiss Med. Wkly. 150, w20248 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4414/smw.2020.20248&link_type=DOI) 

18. 18.Nguyen, A. et al. Human Leukocyte Antigen Susceptibility Map for Severe Acute Respiratory Syndrome Coronavirus 2. J. Virol. 94, e00510–20, /jvi/94/13/JVI.00510-20.atom (2020).
    
    

19. 19.Saulle, I., Vicentini, C., Clerici, M. & Biasin, M. Antigen presentation in SARS-CoV-2 infection: the role of class I HLA and ERAP polymorphisms. Hum. Immunol. 82, 551–560 (2021).
    
    

20. 20.Ahmed, S. F., Quadeer, A. A. & McKay, M. R. Preliminary Identification of Potential Vaccine Targets for the COVID-19 Coronavirus (SARS-CoV-2) Based on SARS-CoV Immunological Studies. Viruses 12, 254 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/vl2030254&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

21. 21.Jp, R.-L. et al. Prediction of SARS-CoV2 spike protein epitopes reveals HLA-associated susceptibility. Res. Prepr. (2020) doi:10.21203/rs.3.rs-25844/v1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21203/rs.3.rs-25844/v1&link_type=DOI) 

22. 22.Wang, W., Zhang, W., Zhang, J., He, J. & Zhu, F. Distribution of HLA allele frequencies in 82 Chinese individuals with coronavirus disease-2019 (COVID-19). HLA 96, 194–196 (2020).
    
    

23. 23.Novelli, A. et al. HLA allele frequencies and susceptibility to COVID-19 in a group of 99 Italian patients. HLA tan.14047 (2020) doi:10.1111/tan.14047.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/tan.14047&link_type=DOI) 

24. 24.Lorente, L. et al. HLA genetic polymorphisms and prognosis of patients with COVID-19. Med. Intensiva S0210569120302667 (2020) doi:10.1016/j.medin.2020.08.004.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.medin.2020.08.004&link_type=DOI) 

25. 25.Warren, R. L. & Birol, I. HLA predictions from the bronchoalveolar lavage fluid and blood samples of eight COVID-19 patients at the pandemic onset. Bioinformatics btaa756 (2020) doi:10.1093/bioinformatics/btaa756.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btaa756&link_type=DOI) 

26. 26.Correale, P. et al. HLA-B*44 and C*01 Prevalence Correlates with Covid19 Spreading across Italy. Int. J. Mol. Sci. 21, 5205 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms21155205&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32717807&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

27. 27.Sakuraba, A., Haider, H. & Sato, T. Population Difference in Allele Frequency of HLA-C*05 and Its Correlation with COVID-19 Mortality. Viruses 12, 1333 (2020).
    
    

28. 28.Langton, D. J. et al. The influence of HLA genotype on the severity of COVID-19 infection. HLA 98, 14–22 (2021).
    
    

29. 29.Tai, D. B. G., Shah, A., Doubeni, C. A., Sia, I. G. & Wieland, M. L. The Disproportionate Impact of COVID-19 on Racial and Ethnic Minorities in the United States. Clin. Infect. Dis. 72, 703–706 (2021).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

30. 30.Données sur les variants du SRAS-CoV-2 au Québec. INSPQ [https://www.inspq.qc.ca/covid-19/donnees/variants](https://www.inspq.qc.ca/covid-19/donnees/variants).
    
    

31. 31.Reid, M. E., Lomas-Francis, C. & Olsson, M. L. The Blood Group Antigen FactsBook. (Academic press, 2012).
    
    

32. 32.Gragert, L., Madbouly, A., Freeman, J. & Maiers, M. Six-locus high resolution HLA haplotype frequencies derived from mixed-resolution DNA typing for the entire US donor registry. Hum. Immunol. 74, 1313–1320 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.humimm.2013.06.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23806270&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

33. 33.Nunes, J. M., Buhler, S., Roessli, D., Sanchez-Mazas, A., & the HLA-net 2013 collaboration. The HLA-net GENE[RATE] pipeline for effective HLA data analysis and its application to 145 population samples from Europe and neighbouring areas. Tissue Antigens 83, 307–323 (2014).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 

34. 34.Weir, B. S. & Goudet, J. A Unified Characterization of Population Structure and Relatedness. Genetics 206, 2085 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6MTA6IjIwNi80LzIwODUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMS8yMS8yMDIxLjAzLjE3LjIxMjUzODIxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

35. 35.Takezaki, N. & Nei, M. Genetic Distances and Reconstruction of Phylogenetic Trees From Microsatellite DNA. Genetics 144, 389 (1996).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6OToiMTQ0LzEvMzg5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDEvMjEvMjAyMS4wMy4xNy4yMTI1MzgyMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

36. 36.R Core Team. R: A Language and Environment for Statistical Computing. in R Foundation for Statistical Computing vol. 1 (2013).
    
    

37. 37.Al-Youha, S. A. et al. The impact of ABO blood groups on clinical outcomes and susceptibility to COVID-19: A retrospective study in an unselected population. Transfusion (Paris) 61, 1631–1641 (2021).
    
    

38. 38.Trépanier, P. et al. Adapting to supply-and-demand emerging trends for antigen-negative red blood cell units. Transfusion (Paris) trf.16285 (2021) doi:10.1111/trf.16285.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/trf.16285&link_type=DOI) 

39. 39.Tournamille, C., Colin, Y., Cartron, J. P. & Le Van Kim, C. Disruption of a GATA motif in the Duffy gene promoter abolishes erythroid gene expression in Duffy–negative individuals. Nat. Genet. 10, 224–228 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng0695-224&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7663520&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F21%2F2021.03.17.21253821.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RA85700022&link_type=ISI) 

40. 40.Hebbel, R. P. & Vercellotti, G. M. SARS-CoV-2 severity in African Americans – A role for Duffy Null? Haematologica 105, 2892 (2020).
    
    

41. 41.Szymanski, J. et al. ABO blood type association with SARS-CoV-2 infection mortality: A single-center population in New York City. Transfusion (Paris) trf.16339 (2021) doi:10.1111/trf.16339.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/trf.16339&link_type=DOI) 

42. 42.Migliorini, F. et al. Association between HLA genotypes and COVID-19 susceptibility, severity and progression: a comprehensive review of the literature. Eur. J. Med. Res. 26, 84 (2021).
    
    

43. 43.Tavasolian, F. et al. HLA, Immune Response, and Susceptibility to COVID-19. Front. Immunol. 11, 3581 (2021).
    
    

44. 44.Tomita, Y., Ikeda, T., Sato, R. & Sakagami, T. Association between HLA gene polymorphisms and mortality of COVID-19: An in silico analysis. Immun. Inflamm. Dis. iid3.358 (2020) doi:10.1002/iid3.358.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/iid3.358&link_type=DOI) 

45. 45.Shkurnikov, M. et al. Association of HLA Class I Genotypes With Severity of Coronavirus Disease-19. Front. Immunol. 12, 423 (2021).
    
    

46. 46.Pisanti, S. et al. Correlation of the two most frequent HLA haplotypes in the Italian population to the differential regional incidence of Covid-19. J. Transl. Med. 18, 352 (2020).
    
    

47. 47.Barquera, R. et al. Binding affinities of 438 HLA proteins to complete proteomes of seven pandemic viruses and distributions of strongest and weakest HLA peptide binders in populations worldwide. HLA 96, 277–298 (2020).
    
    

48. 48.Vishnubhotla, R., Sasikala, M., Ketavarapu, V. & Reddy, D. N. High-resolution HLA genotyping identifies alleles associated with severe COVID-19: A preliminary study from India. Immun. Inflamm. Dis. 9, 1781–1785 (2021).
    
    

49. 49.Langton, D. J. et al. The Influence of Hla Genotype on Susceptibility to, and Severity of, Covid-19 Infection. Medrxiv 2020.12.31.20249081 (2021) doi:10.1101/2020.12.31.20249081.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMi4zMS4yMDI0OTA4MXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDEvMjEvMjAyMS4wMy4xNy4yMTI1MzgyMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)