THE SEARCH FOR AN ASSOCIATION OF HLA ALLELES AND COVID-19 RELATED MORTALITY IN THE RUSSIAN POPULATION ===================================================================================================== * Valery Cheranev * Irina Bulusheva * Valery Vechorko * Dmitriy Korostin * Denis Rebrikov ## Abstract HLA genes play a pivotal role in an immune response via the presentation of pathogen peptides in a complex on the surface of cells of a host organism. Here, we studied the association of class I and class II genes with the severity of COVID-19 infection and HLA allele variants. We performed high-resolution sequencing of class I and class II HLA genes using the sample population of 147 patients who died of COVID-19 and statistically compared our results with the frequencies of the HLA genotypes in a control population of 270 samples. The obtained data demonstrated that 51:05 and 15:18 alleles from locus B* are statistically significantly associated with COVID-19 severity, while C*14:02 allele correlates with the probability of death from COVID-19 for patients without comorbidities. ## Introduction In 2020, the human population faced an epidemic that affected almost all countries to a certain degree. It was caused by a new virus from the Coronaviridae family SARS-CoV-2 [1]. The first outbreak of the disease was registered in Wuhan (China) and the disease rapidly spread over the other countries. Currently, the number of registered cases exceeds 53 million, 1.3 million being lethal [2]. In January 2020, the World Health Organization (WHO) declared the outbreak to be a pandemic and termed the disease caused by the virus COVID-19. One of the essential factors facilitating the rapid spread of COVID-19 is a difference in the severity of a disease and its manifestations widely ranging in individual patients [3, 4, 5]. The accumulated data on SARS-CoV-2 suggest some correlation between the disease outcome and sex, age, and concurrent diseases [6]. Disease severity can be also associated with individual genetic characteristics of a patient [7, 8]. One of the genetic predictors for severity and an outcome of the disease could be class I and II HLA genes encoding the proteins of the major histocompatibility complex (MHC). The main function of these molecules is presenting an antigen on the surface of the plasma membrane for their recognition by immune cells. Therefore, the alleles encoding the aminoacid sequence in MHC directly affect the susceptibility to certain diseases. Based on the data on the alleles of class I and II HLA genes and the protein composition of SARS-CoV-2, we analyzed the level of the affinity of MHC binding to all possible viral epitopes [9]. The lowest predicted level of the interaction with viral antigens belonged to the protein encoded by B*46:01 allele while the highest level belonged to the B*15:03 allele. Thus, we proceeded with a retrospective analysis aimed to reveal any association between the identified alleles and disease severity [10-15]. Previously, the correlation between HLA alleles and the disease in the test and control groups was performed using statistical analysis by the 2-level procedure of identification of significant alleles [16]. On the first stage, the distribution of all alleles in patients was compared to the distribution of the reference population in order to choose several candidate alleles. After that, candidate alleles were compared between two groups of patients using the Bonferroni method. The correlation between HLA and the disease was studied at the background of concurrent diseases [17, 18, 19]. In our work, we compared the frequencies of different alleles of class I (A, B, C) and class II (DRB1, DQB1) HLA genes between the groups of patients and healthy donors. The main goal of the work is to identify and validate the alleles significantly related to various COVID-19 outcomes. ## Background ### Materials and methods #### Ethics statement This study conformed to the principles of the Declaration of Helsinki. The appropriate institutional review board approval for this study was obtained from the Ethics Committee at Pirogov Medical University. All patients provided written informed consent for the collection of samples, subsequent analysis, and publication thereof. #### Biomaterial collection The biomaterial consisting of whole venous blood was collected in EDTA-coated tubes used as a preservative. #### gDNA isolation gDNA was isolated from 100 ul of whole venous blood with Proba-Mch-Maks reagent kit using automated dosing station дТстрИм (DT stream) (DNA Technology, Russia). Quality check of the isolated DNA was performed using gel electrophoresis in agarose gel, the concentration was measured using Qubit 2 fluorimeter with Qubit dsDNA BR Assay kit (ThermoFisher Scientific, USA). #### HLA library preparation The preparation of amplicon libraries for HLA high-resolution genotyping was performed using HLA Expert kit (DNA Technology, Russia), following the manufacturer’s protocol. #### Sequencing Sequencing was performed with Illumina MiSeq using MiSeq Reagent Kit v3 (600-cycle), according to the manufacturer’s instructions. #### HLA high-resolution genotyping The analysis of fastq files was performed with HLA-Expert software v.2.0 (DNA Technology, Russia) according to the manufacturer’s protocol. #### Statistical Analysis Statistical analysis was performed using Arlequin v.3.5.2.2 [16] and Pearson’s goodness-of-fit test. We created several scripts allowing for estimating the diversity of each gene and differences in frequencies of the individual alleles in groups. We also created a script that cleared an input table containing patients’ data from errors and transformed the names of HLA alleles, according to a unified syntax ([https://github.com/genomecenter/HLA_article](https://github.com/genomecenter/HLA_article)). We created a script that generated an input file containing patients’ data for Arlequin. ## Results As a study population, we used a collection of venous blood samples and clinical data from 147 patients (group III) that had died in Moscow clinics from coronavirus disease COVID-19 during the pandemic in 2020. The positive test for COVID-19 was confirmed using qPCR kits in Moscow clinical diagnostic laboratories that collected the biomaterial. Clinical features of the group III is presented in Table 1. As the control sample population, we used 270 venous blood samples collected from the members of the National Registry of Bone Marrow Donors by the Pirogov Medical University at the beginning of 2020. All members of the control sample population were surveyed and stratified into two groups: those who hadn’t suffered from COVID-19 during our study (245 samples, group 1) and those who had recovered from COVID-19. View this table: [Table 1.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T1) Table 1. Clinical data in the cohort of deceased patients with COVID-19. The data from HLA high-resolution genotyping for each patient included the information on both alleles of A*, B*, C*, DRB1*, DQB1* genes in the HLA histocompatibility complex (Group I-III HLA genotypes are presented in Supplemental Table 1). We performed this analysis for every gene and for each allele of a gene that had significant differences. The distribution of allele frequencies over 5 loci from three groups is presented in Figure 1. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/24/2020.12.22.20248695/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/F1) Figure 1. **The histogram of the frequencies of alleles of class I (A) and class II (B) HLA loci. The horizontal axis, the value of the 1st level allele. The vertical axis, its frequency in a sample population. Group I is denoted by yellow bars, group II is denoted by orange bars, group III is denoted by green bars. Starts denote the alleles that showed statistically significant associations**. For three groups, we used a method of distances upon pairwise comparison [21] using the fixation index FST which is a special case of the Wright F-statistics based on the correlation of randomly chosen alleles within a population with the whole population which is equivalent to the ratio of genetic diversity due to the differences in allele frequency in populations [22]. Results are presented in Table 2. Zero and negative values of FST mean that there was no genetic stratification between populations Group I and III are similar, group II is equidistant from both of them. View this table: [Table 2.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T2) Table 2. The results obtained by the distance method upon pairwise comparisons of groups I-III. After that, we performed the population Assignment Test (Figure 2). For each sample population, we created a probability graph of the genotype within the sample from the population compared to the other populations. The results obtained are similar to those obtained from a pairwise comparison. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/24/2020.12.22.20248695/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/F2) Figure 2. **Results from the population Assignment Test. For each group (I – A, II – B, III – C), we created the genotype probability plot within the sample from the population compared to all other populations**. Alternatively, we used Pearson’s goodness-of-fit test that enabled a separate analysis of each gene. We analysed the 2-level alleles for each gene. The null hypothesis stated that HLA did not affect the divergence in the groups. We chose a p-value of 0.05 as a threshold value. We compiled a contingency matrix for each gene using one of four methods (Table 3). Although considering both alleles from each patient as independent ones (V2) is incorrect from the point of view of statistical analysis, we decided to perform this analysis in order to compare the results. View this table: [Table 3.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T3) Table 3. Three ways to collect patient’s data to compile a contingency matrix for each gene. Analysis of the combined groups: I+II vs III (i.e. healthy or recovered donors vs patients who died from COVID-19), II+III vs I (i.e. all who suffered from COVID-19 vs healthy donors) didn’t reveal any significant differences. Separate analysis of three groups revealed statistically significant effects of HLA-C* allele combinations (Table 4). As a multiple comparison correction, we used the Bonferroni method [23] with the significance level of 0.05 Locus C* significantly affects the group stratification. To determine significant alleles for each gene based on the way of compiling a contingency matrix, we created a special script. The results obtained for V3 are presented in Table 5. For alleles we considered only once, we chose the V4 method, as the results from the V1 method should be considered in the iterative dynamics due to the random choice of an allele. The results obtained from V4 are presented in Table 6. View this table: [Table 4.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T4) Table 4. Comparison of groups I, II, and III of samples based on Pearson’s goodness-of-fit test. View this table: [Table 5.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T5) Table 5. **The list of statistically significant combinations of HLA gene alleles that define the differences between groups I, II, and III following the V3 method**. View this table: [Table 6.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T6) Table 6. **The list of statistically significant combinations of HLA gene alleles that define the differences between groups I, II, and III following the V4 method**. To determine the significant alleles, we analysed the other group combinations, according to the V4 method. In particular, the group of deceased patients was divided into subgroups, according to their age (age not exceeding 65 or exceeding 65 at the time of death) and comorbidities. The group designations and sizes are presented in Table 7. The results obtained from the other combinations are presented in View this table: [Table 7.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T7) Table 7. The designations and number of the groups subjected to analysis. Table 8. We should note that we had to exclude some patients from the sample population (Table 1) due to the lack of full clinical data. The distribution of allele frequencies of the locus C* in the subgroups IIIA-IIID is presented in Figure 3. View this table: [Table 8.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T8) Table 8. Comparison of the group combinations analyzed by V4 method. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/24/2020.12.22.20248695/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/F3) Figure 3. **The distribution of allele frequencies of locus C* in the patients who died from COVID-19. The subgroups are indicated in the diagram**. ## Discussion Thus, locus C* statistically significantly affects the groups of patients who didn’t suffer from COVID-19, recovered and died (Table 4); estimating the significance of allele combinations reveal statistically significant associations with COVID-19 severity (Table 5). Locus B 51:05 and 15:18 alleles show statistically significant association with the COVID-19 severity (Table 6), while the p-value of C*14:02 and DQB1*06:01 alleles are close to the threshold value. Upon stratification of the patients who died from COVID-19 according to their age, C*14:02 allele was significant (Table 8). Importantly, the manipulations with group stratification based on various parameters allow determining a criterion that would make the stratification significant. This is related to a limited size of the experimental sample population. Therefore, we are planning to confirm the results obtained in this work using a larger sample that is being collected at the moment. We hope to obtain conclusive evidence regarding the C*14:02 and DQB1*06:01 alleles. We performed a comparative analysis of the previously published results on the search for the association between HLA and COVID-19 severity (Table 9). Statistically significant and close to statistically significant B*51, C*14:02, and DQB1*06:01 alleles we singled out in our work turned out to be identified in other works as well. Some divergences in the results may be accounted for distinct frequencies of HLA alleles within the particular ethnic groups. View this table: [Table 9.](http://medrxiv.org/content/early/2020/12/24/2020.12.22.20248695/T9) Table 9. Comparison of the results on alleles for COVID-19 obtained by other investigators ## Conclusions Locus B*51:05 and 15:18 alleles show statistically significant association with COVID-19 severity upon comparing the sample populations of patients in Russia who died from coronavirus infection with the group of bone marrow donors. They are more frequent in patients who recovered from COVID-19. The C*14:02 allele revealed a statistically significant correlation with the increased probability of death from COVID-19 of patients without comorbidities. ## Supporting information Supplemental Table 1 [[supplements/248695_file05.xlsx]](pending:yes) ## Data Availability All received data are available upon request ## Funding This research was funded by grant 𝒩º075-15-2019-1789 from the Ministry of Science and Higher Education of the Russian Federation allocated to the Center for Precision Genome Editing and Genetic Technologies for Biomedicine. ## Conflicts of Interest The authors declare no conflict of interest. ## Acknowledgments * Received December 22, 2020. * Revision received December 22, 2020. * Accepted December 23, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. 1.Meyerowitz-Katz G, Merone L. A systematic review and meta-analysis of published research data on COVID-19 infection fatality rates. Int J Infect Dis. 2020 Sep 29; 101:138–148. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F24%2F2020.12.22.20248695.atom) 2. 2.WHO. Coronavirus disease (COVID-19) Weekly Epidemiological Update and Weekly Operational Update: Weekly Epidemiological Update 13 November 2020 [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports) (2020). 3. 3.Mao, Ren, et al. “Manifestations and prognosis of gastrointestinal and liver involvement in patients with COVID-19: a systematic review and meta- analysis. ” The lancet Gastroenterology & hepatology (2020). 4. 4.Spinato, Giacomo, et al. “Alterations in smell or taste in mildly symptomatic outpatients with SARS-CoV-2 infection.” Jama (2020). 5. 5.Vetter, P. et al. Clinical features of COVID-19. Br. Med. J. 369, m1470 (2020). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjkvYXByMTdfMTYvbTE0NzAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMi8yNC8yMDIwLjEyLjIyLjIwMjQ4Njk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 6. 6.Zhou, F. et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet 395, 1054–1062 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30566-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F24%2F2020.12.22.20248695.atom) 7. 7.Ellinghaus, D. et al. Genomewide association study of severe COVID-19 with respiratory failure. N. Engl. J.Med. [https://doi.org/10.1056/NEJMoa2020283](https://doi.org/10.1056/NEJMoa2020283) (2020). 8. 8.Cascella M, Rajnik M, Cuomo A, et al. Features, Evaluation, and Treatment of Coronavirus. [Updated 2020 Oct 4]. In: StatPearls [Internet]. Treasure Island (FL): StatPearls Publishing; 2020 Jan-. Available from: [https://www.ncbi.nlm.nih.gov/books/NBK554776/](https://www.ncbi.nlm.nih.gov/books/NBK554776/) 9. 9.Nguyen, A., David, J. K., Maden, S. K., Wood, M. A., Weeder, B. R., Nellore, A., & Thompson, R. F. (2020). Human leukocyte antigen susceptibility map for SARS-CoV-2. Journal of virology. 10. 10.Iturrieta-Zuazo, I., Rita, C. G., García-Soidán, A., de Malet Pintos-Fonseca, A., Alonso-Alarcón, N., Pariente-Rodríguez, R., Tejeda-Velarde, A., Serrano-Villar, S., Castañer-Alabau, J. L., & Nieto-Gañán, I. (2020). Possible role of HLA class-I genotype in SARS-CoV-2 infection and progression: A pilot study in a cohort of Covid-19 Spanish patients. Clinical immunology (Orlando, Fla.), 219, 108572. 11. 11.Warren, R. L., & Birol, I. (2020). Retrospective in silico HLA predictions from COVID-19 patients reveal alleles associated with disease prognosis. medRxiv : the preprint server for health sciences, 2020.10.27.20220863. 12. 12.Pisanti, S., Deelen, J., Gallina, A. M., Caputo, M., Citro, M., Abate, M., Sacchi, N., Vecchione, C., & Martinelli, R. (2020). Correlation of the two most frequent HLA haplotypes in the Italian population to the differential regional incidence of Covid-19. Journal of translational medicine, 18(1), 352. 13. 13.Warren, R. L., & Birol, I. (2020). HLA predictions from the bronchoalveolar lavage fluid samples of five patients at the early stage of the Wuhan seafood market COVID-19 outbreak. ArXiv, 2004.07108v3. 14. 14.Wang, W., Zhang, W., Zhang, J., He, J., & Zhu, F. (2020). Distribution of HLA allele frequencies in 82 Chinese individuals with coronavirus disease-2019 (COVID-19). HLA, 96(2), 194–196. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/tan.13941&link_type=DOI) 15. 15.Barquera, R., Collen, E., Di, D., Buhler, S., Teixeira, J., Llamas, B., Nunes, J. M., & Sanchez-Mazas, A. (2020). 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(3), 277–298. 16. 16.Goverdhan, S. V., Howell, M. W., Mullins, R. F., Osmond, C., Hodgkins, P. R., Self, J., … & Lotery, A. J. (2005). Association of HLA class I and class II polymorphisms with age-related macular degeneration. Investigative ophthalmology & visual science, 46(5), 1726–1734. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo5OiI0Ni81LzE3MjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMi8yNC8yMDIwLjEyLjIyLjIwMjQ4Njk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 17. 17.Teixeira, S. L. M., De Sá, N. B. R., Campos, D. P., Coelho, A. B., Guimaraes, M. L., Leite, T. C. N. F., … & Morgado, M. G. (2014). Association of the HLA-B* 52 allele with non-progression to AIDS in Brazilian HIV-1-infected individuals. Genes & Immunity, 15(4), 256–262. 18. 18.Karnes, J. H., Shaffer, C. M., Bastarache, L., Gaudieri, S., Glazer, A. M., Steiner, H. E., … & Roden, D. M. (2017). Comparison of HLA allelic imputation programs. PLoS One, 12(2), e0172444. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0172444&link_type=DOI) 19. 19.Ding, S. J., Zhang, Y., Zhang, X. M., Jiang, X. L., Pang, B., Song, Y. H., … & Yu, X. J. (2016). Correlation Between HLA-A, B and DRB1 Alleles and Severe Fever with Thrombocytopenia Syndrome. PLoS neglected tropical diseases, 10(10), e0005076. 20. 20.Excoffier, L., & Lischer, H. E. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular ecology resources, 10(3), 564–567. 21. 21.Beerli, P. (2005). Pairwise distance methods. Computational Evolucionary Biology, No BSC5936-Fall, 1–7. 22. 22.Holsinger, K. E., & Weir, B. S. (2009). Genetics in geographically structured populations: defining, estimating and interpreting F ST. Nature Reviews Genetics, 10(9), 639–650. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg2611&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19687804&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F24%2F2020.12.22.20248695.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269063200013&link_type=ISI) 23. 23.Bland, J. M., & Altman, D. G. (1995). Multiple significance tests: the Bonferroni method. Bmj, 310(6973), 170. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjEyOiIzMTAvNjk3My8xNzAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMi8yNC8yMDIwLjEyLjIyLjIwMjQ4Njk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 24. 24.Wang, F., Huang, S., Gao, H., Zhou, Y., Lai, C., Li, Z., … & Tang, Q. (2020). Initial Whole Genome Sequencing and Analysis of the Host Genetic Contribution to COVID-19 Severity and Susceptibility. medRxiv. 25. 25.Nguyen, A., David, J. K., Maden, S. K., Wood, M. A., Weeder, B. R., Nellore, A., & Thompson, R. F. (2020). Human leukocyte antigen susceptibility map for SARS-CoV-2. Journal of virology. 26. 26.Lin, M., Tseng, H. K., Trejaut, J. A., Lee, H. L., Loo, J. H., Chu, C. C., … & Lin, R. Y. (2003). Association of HLA class I with severe acute respiratory syndrome coronavirus infection. BMC Medical Genetics, 4(1), 9. 27. 27.Ng, M. H., Lau, K. M., Li, L., Cheng, S. H., Chan, W. Y., Hui, P. K., … & Sung, J. J. (2004). Association of human-leukocyte-antigen class I (B* 0703) and class II (DRB1* 0301) genotypes with susceptibility and resistance to the development of severe acute respiratory syndrome. Journal of Infectious Diseases, 190(3), 515–518. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/421523&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15243926&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F24%2F2020.12.22.20248695.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222549500014&link_type=ISI) 28. 28.Liu, D., Qiu, Y., Zha, Y., Li, W., Li, D., & Wu, T. (2018). Association of HLA class?and class?genes with severe acute respiratory syndrome in the northern Chinese population. 29. 29.Sanchez-Mazas, A. (2020). HLA studies in the context of coronavirus outbreaks. Swiss Medical Weekly, 150(1516).