Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 09 June 2022
Sec. Immunogenetics
This article is part of the Research Topic Progress in understanding the immunogenetic basis of disease susceptibility and outcomes View all 5 articles

HLA Diversity in Saudi Population: High Frequency of Homozygous HLA Alleles and Haplotypes

Aziz Alami Chentoufi,
Aziz Alami Chentoufi1,2*F. Aytül UyarF. Aytül Uyar3Hiba A. ChentoufiHiba A. Chentoufi4Khalid AlzahraniKhalid Alzahrani1Maria PazMaria Paz1Ahmed BahnassyAhmed Bahnassy1Ghaleb ElyamanyGhaleb Elyamany5Assem Elghazaly,Assem Elghazaly1,6
  • 1King Fahad Medical City, Riyadh, Saudi Arabia
  • 2Department of Medical Microbiology, Faculty of Health Sciences, University of Pretoria, Pretoria, South Africa
  • 3Department of Physiology, Istanbul Medical Faculty, Istanbul University, Istanbul, Turkey
  • 4Biomedical Sciences, University of Ottawa, Ottawa, ON, Canada
  • 5Department of Central Military Laboratory and Blood Bank, Prince Sultan Military Medical City, Riyadh, Saudi Arabia
  • 6Department of Hematology/Oncology, King Saud Medical City, Riyadh, Saudi Arabia

Human leukocyte antigens (HLA) diversity has a tremendous impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host–pathogen interactions. Here, we conducted a retrospective study of HLA class I and class II homozygosity at allelic and haplotype levels in unrelated individuals genotyped from 2012 to 2016 in a tertiary hospital in the capital of Saudi Arabia. Among 5,000 individuals, 2,773 individuals meet inclusion criteria and were retrospectively analyzed for HLA-A, -B, -C–DRB1, and -DQB1 homozygosity at allelic and haplotype levels. HLA molecular typing was performed using a commercial reverse sequence-specific oligonucleotide (rSSO) kit. We were able to identify 15 HLA-A, 20 HLA-B, 11 HLA-C, 13 HLA-DRB1, and five HLA-DQB1 homozygous alleles demonstrating a very low genetic diversity in the Saudi population. The highest homozygosity in HLA class I was found in locus C followed by A and B (20.3% > 16.1% > 15.5%; p < 0.001) where the most homozygote alleles were A*02 (9.2%), B*51 and B*50 (5.7% and 3.7%), and C*07, C*06, and C*15 (7.2%, 5.48%, and 3.3%) and in HLA class II, the highest homozygosity was found in locus DQB1 compared to DRB1 (31.71% > 19.2%; p < 0.001), with the most common homozygote alleles being DRB1*07 and DRB1*04 (5.33% and 4.2%) and DQB1*02, DQB1*06, and DQB1*03 (13.55%, 7.92%, and 7.64%). The frequency of finding an individual with one homozygote allele was (24.6%), two homozygote alleles (13.5%), three homozygote alleles (4.7%), four homozygote alleles (3.4%), and five alleles were (4.8%). The most frequent homozygote haplotypes are A*23∼C*06∼B*50∼DRB1*07∼DQB1*02 and A*02∼C*06∼B*50∼DRB1*07∼DQB1*02. This study shows low diversity of both class I and II alleles and haplotypes in the Saudi population, which would have a significant impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host–pathogen interactions.

Introduction

Human leukocyte antigens (HLA) molecules, also known as MHC molecules (major histocompatibility complex), are encoded by the most polymorphic cluster of genes in the human genome located on the short arm of chromosome 6 (Marrack and Kappler, 1986). Classically divided into HLA class I (HLA-A, -B, and -C) and HLA class II (HLA-DR, -DQ, and -DP) with more than 27.000 HLA alleles sequences (Robinson et al., 2020) encoding for molecules critical to immunity, including antigen processing and presentation to immune effector CD4+ and CD8+ T-cells in an HLA-restricted manner, as such, HLA molecules have been extensively studied in the context of organ and hematopoietic cell transplantation and infectious and autoimmune disease associations (Shiina et al., 2004; Dehn et al., 2008; Trowsdale and Knight, 2013; Dendrou et al., 2018).

Up to date, HLA class I and class II alleles and haplotype frequencies were characterized in different populations around the globe including the Saudi population (Maiers et al., 2007; Saw et al., 2015; Creary et al., 2019; Jawdat et al., 2020). Understanding HLA diversity in these populations has a tremendous impact on clinical outcomes of solid organ and hematopoietic cell transplantation where a mismatch in HLA class I and/or class II alleles between donors and recipients represent a major cause for graft failure, rejection, and graft versus host disease (GVHD) (Lee et al., 2007; Dehn et al., 2008; Petersdorf, 2008; Kunze-Schumacher et al., 2014). In addition, HLA diversity presents an important challenge in transfusing medicine where HLA homozygosity represents the major cause of transfusion-associated-graft versus host disease (TA-GVHD) in immunocompetent patients. TA-GvHD is a very rare but fatal complication that occurs following transfusion of lymphocyte-containing blood components such as packed red blood cells and platelets (Treleaven et al., 2011).

The HLA genetic diversity has been shown to play a critical role in developing an effective host response against pathogens, disease progression, and outcome. Indeed, it has been shown that HLA heterozygous individuals have an advantage over homozygous individuals in presenting pathogen-derived epitopes to effector CD8+ and CD4+ T-cells (Doherty and Zinkernagel, 1975; Thursz et al., 1997; Penn et al., 2002). Multiple infectious disease studies have confirmed the negative clinical impact of HLA homozygosity on disease development and outcome (Thursz et al., 1997; Carrington et al., 1999; Burgner et al., 2006). Similarly, it has been shown for non-infectious neoplasms that HLA homozygosity reduces the ability of the immune system to differentiate cancer cells from noncancer cells (Hojjat-Farsangi et al., 2008; Shah et al., 2011).

The Kingdom of Saudi Arabia (KSA) constitutes the majority of the Arab Peninsula in southwest Asia, with a land area of about 2,000,000 km2. According to the 2020 survey, the KSA population is 35 million with an annual growth rate of 2.4. The gender ratio of the Saudi population is 1.04 (General Authority of Statistics, 2016). At least 24.8% of the whole population is less than 15 years and people in the working age who range between 15–64 years old represent 72%. The Saudi Population consists of a multi-ethnic population, of which the largest portion (90% of which) constitutes Arab. KSA is characterized by a large migrant population, which forms 36.3% of its residents, the majority of whom are of Asian and African ethnicities. 26.3% of immigrants are native to India (General Authority of Statistics, 2016). Most of the marriage is between Saudi people, and the overall rate of consanguinity is about 57.7%, one of the highest worldwide (El-Hazmi et al., 1995), a trend that did not show a changeover generation (Warsy et al., 2014).

The aim of this work is to study the prevalence of homozygosity of HLA class I and Class II alleles and haplotypes which would reflect the degree of genetic diversity (or similarity) in the Saudi population. This study would pave the way to better transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host–pathogen interactions in the Saudi population.

Materials and Methods

Cohort Study

This is a retrospective study, using already available HLA typing results of Saudi individuals (HSCT-unrelated donors). The study was performed on 5,000 individuals, who have been HLA-typed for HLA-class I (A, B, and C) and HLA-class II (DR and DQ) in the HLA laboratory at KFMC between January 2012 and June 2016. Only 2,773 unrelated individuals were included in this cohort study representing Saudi Arab ethnicity. Saudi ethnicity was defined as Arabs belonging to known Saudi tribes and citizens of the kingdom of Saudi Arabia. The study population shows that they were in the age range of <25 years: 15%; 25–45 years: 55%; and 45–76 years: 30%. The male: female ratio of patients was 51%:49%. The retrospective study was approved for human investigations at King Fahad Medical Hospital Tertiary Hospital under an institutional review board–approved protocol (IRB# 15-198).

Blood Samples, DNA Extraction, and Quantification

Blood samples were drawn in the ward in 2-ml EDTA tubes and sent to the HLA laboratory for HLA typing at KFMC. For DNA extraction from blood samples, an automated MagNA Pure Compact System (Roche Life Sciences, Indianapolis, United States) was used. The protocol was performed as per the manufacturer’s instructions without modifications. DNA extraction: The extracted DNA was quantified using Nanodrop 2000 C (Thermo Scientific, Wilmington, DE, United States) and normalized to the desired concentration as recommended by the manufacturer’s protocol.

HLA-A, -B, -C, -DRB1, and -DQB1 Genotyping Using rSSO

HLA-A typing was performed using a commercial reverse sequence-specific oligonucleotides kit (SSO; One Lambda, Canoga Park, CA). Briefly, extracted and quantified genomic DNA from the blood has been subjected to a PCR reaction using primers specific for areas of the exon 2 and 3 for HLA class I (HLA-A, B, and C) and the exon 2 for HLA-DRB1 and HLA-DQB1. The PCR product was hybridized with a mixture of beads (n = 100), each one having a specific sequence and color. Then, the beads were acquired on the Luminex platform and analyzed by HLA Fusion Software (One Lambda, Inc., Canoga Park, CA). To make sure that the obtained homozygosity is not false results commonly obtained with rSSO, we confirmed our results by family segregation study and/or sequence-based typing for the concerned homozygous allele.

Statistical Analysis

The HLA class I and class II allele homozygosity frequencies were calculated by dividing the number of homozygote individuals for this allele by the total number of homozygote individuals at that locus. Allele homozygosity population frequency was calculated by dividing the number of homozygote individuals for this allele by the total number of studied individuals. All frequencies were calculated from an Excel spreadsheet counts and using Fisher’s exact test (two-tailed) to calculate the p value (p value < 0.05 is considered significant).

The population analyses and genetic diversity measures were calculated using Arlequin 3.5 software. Allele frequencies and five locus haplotype frequencies were determined by the expectation-maximization (EM) algorithm. For testing the deviation from Hardy Weinberg equilibrium, a Markov chain algorithm was performed (Excoffier and Lischer, 2010). Alleles are said to be in Hardy Weinberg equilibrium when the observed heterozygote frequencies do not significantly differ from the expected frequencies (p > 0.05).

Results

HLA Class I Allele Frequency Distribution of Homozygosity

In our cohort, 15 HLA-A, 20 HLA-B, and 11 HLA-C homozygous alleles were identified (Tables 13). Of the 15 HLA-A homozygous alleles identified, only one with the frequency of homozygosity above 5% was found (Table 1). This allele was HLA-A*02 (9.2%). Similarly, of the 20 HLA-B homozygous alleles identified in this population, only one allele B*51 (5.7%) showed a frequency of homozygosity above 5% (Table 2). Among the 11 HLA-C homozygous alleles identified in this population, two alleles C*07 and C*06 (7.2% and 5.48%) showed a frequency of homozygosity above 5% (Table 3). Comparatively, as shown in (Figure 1A), HLA-C showed the highest frequency of homozygosity followed by HLA-A and HLA- B (20.3% > 16.1% > 15.5%; p < 0.001). In order to analyze the influence of gender on HLA allele frequencies, we have calculated HLA class I and class II allele frequencies in males and females. A comparison between male and female allele frequencies of HLA-A, -B, -C, - DRB1, and DQB1 loci show that there were no significant differences between male and female in term of allele frequencies, demonstrating that gender has a nonsignificant influence on HLA allele frequencies in Saudi population (Supplementary Table S1).

TABLE 1
www.frontiersin.org

TABLE 1. List of HLA-A alleles detected in locus A, the frequency of homozygosity, and allele frequencies.

TABLE 2
www.frontiersin.org

TABLE 2. List of HLA-B alleles detected in locus B and the frequency of homozygosity and allele frequencies.

TABLE 3
www.frontiersin.org

TABLE 3. List of HLA-C alleles detected in locus C and the frequency of homozygosity and allele frequencies.

FIGURE 1
www.frontiersin.org

FIGURE 1. Comparison of the Frequency of homozygosity between HLA class I Loci and HLA class II loci in Saudi population.

HLA Class II Allele Frequency Distribution of Homozygosity

As expected, the HLA class II loci are less diverse than the class I loci in this population. Only a total of 13 HLA-DRB1 and five HLA-DQB1 homozygous alleles (Tables 4, 5) were identified. HLA-DQB Class II homozygous alleles with frequencies above 5% in this population were HLA-DQB1*02, -DQB*06, and -DQB*03 (13.55%, 7.92%, and 7.64%) as shown in Table 4. Only one HLA-DRB1 homozygous allele has shown a frequency above the 5% in the population: HLA-DRB1*07 (5.33%) as shown in Table 5. Comparatively, as shown in Figure 1B, the highest homozygosity in HLA class II is found in locus HLA-DQB1 compared to HLA-DRB1 (31.71% > 19.2%; p < 0.001).

TABLE 4
www.frontiersin.org

TABLE 4. List of HLA-DRB1 alleles detected in locus DRB1 and the frequency of homozygosity and allele frequencies.

TABLE 5
www.frontiersin.org

TABLE 5. List of HLA-DQB1 alleles detected in locus DQB1 and the frequency of homozygosity and allele frequencies.

For each HLA locus, we observed a statistically reduced heterozygosity compared with the expected heterozygosity as shown by Hardy Weinberg equilibrium analysis (Supplementary Table S2). As displayed in Supplementary Table S2, the HWE deviation was observed for the heterozygote component of the HLA-A, -B, -C, -DRB1, and -DQB1 loci data, suggesting that this population may have been positively selected. Excess homozygosity observed in this study is similar to what has been reported previously by Hajeer et al. (2013) and Jawdat et al. (2012).

Class I and Class II Loci and Haplotypes Homozygosity

The frequency of allele homozygosity was overall lower in HLA class I genes than that in class II genes. Among the 2,773 individuals who were fully typed at class I and class II loci, 683 individuals were homozygous at one class I or class II allele (24.6%), 374 individuals (13.5%) were homozygous at two class I and/or HLA class II alleles, 131 individuals (4.7%) were homozygous at three class I and/or class II alleles, 94 individuals (3.4%) were homozygous at four class I and/or HLA class II alleles, and interestingly among the 2,773 individuals who were fully typed at class I and class II loci, 134 (4.8%) were homozygous at the five loci of HLA class I and class II (Table 6). The most frequent homozygote haplotypes are A*23∼C*06∼B*50∼DRB1*07∼DQB1*02 (0.8%) and A*02∼C*06∼B*50∼DRB1*07∼DQB1*02 (0.55%) (Table 7). Saudi population showed higher HLA homozygous haplotype frequency (4.8%) than other populations such as Japanese (1%) and East African (0.96) (Supplementary Table S3) (Kawashima et al., 2012; Peterson et al., 2014).

TABLE 6
www.frontiersin.org

TABLE 6. Frequency to find one, two, three, four, or five homozygote HLA class I and/or alleles in our cohort study.

TABLE 7
www.frontiersin.org

TABLE 7. Frequency of HLA- A ∼ C ∼ B ∼ DRB1∼DQB1 haplotype homozygosity and haplotype frequency in Saudi population.

Discussion

In this study, homozygosity at the HLA class I and class II loci was significantly higher than that in any other population. This low frequency of diversity in the Saudi population would have a significant impact in shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, host–pathogen interaction, and management of recessive mutation expression in the population. HLA class I and II allele frequency distribution in the Saudi population is of prominence in disease association studies, transfusion medicine practices, and vaccine development (Ollier et al., 1985; Sheth et al., 1985; Jawdat et al., 2009; Al-Awwami et al., 2012; Jawdat et al., 2012; Hajeer et al., 2013; Osman et al., 2014). While many genetic studies of the diversity of different populations have been conducted, very few have analyzed HLA class I and class II allele homozygosity (Gladman et al., 1986; Rohowsky-Kochan et al., 1998; Azurdia et al., 1999).

We found here that the Saudi population presents high levels of homozygosity at allele and haplotype levels. The most HLA class I homozygote alleles (>5%) were HLA-A*02, HLA-B*51, and HLA- C*07 and -C*06. HLA-C showed the highest frequency of homozygosity followed by HLA-A and HLA- B (20.3% > 16.1% > 15.5%; p < 0.001). In HLA class II, the most frequent homozygotes alleles were HLA-DRB1*07 and HLA-DQB1*02 and -DQ*06, and -DQ*03. The highest homozygosity in HLA class II is found in locus DQB1 compared to DRB1* (31.71% > 19.2%; p < 0.001). This alleles homozygosity impacted the levels of haplotype homozygosity, with the most homozygous haplotypes in Saudi population being A*23∼C*06∼B*50∼DRB1*07∼DQB1*02 (0.8%) and A*02∼C*06∼B*50∼DRB1*07∼DQB1*02 (0.55%), which are among the most frequent haplotypes in the Saudi population and the Arab population at large (Jawdat et al., 2020; Hajjej et al., 2018; Askar et al., 2013). Hardy Weinberg equilibrium analysis shows an excess of homozygotes in this cohort (Supplementary Table S2). This phenomenon has been previously described (Kanda et al., 2016; Jawdat et al., 2019), which might be explained by the high level of consanguineous marriages among the Saudi population (El-Hazmi et al., 1995).

HLA homozygosity has many advantages in allogeneic bone marrow transplantation such as the limited number of alleles to match and has been shown to give comparable results to haploidentical transplant (Anasetti et al., 1990; Morishima et al., 2020; Fakhoury et al., 2021). In a hematopoietic stem cell transplantation setup, an HLA matching sibling is an ideal donor because of the low risk of GVHD and graft rejection. The incidence of GVHD is proportional to the degree of HLA mismatch (Kanda et al., 2015). The mismatch is in the graft-versus-host direction when a mismatched allele of the donor is homozygous. A mismatched antigen in the GVH direction can be a major target for donor T cells and can cause GVHD, while the mismatch is in the host-versus-graft direction when a mismatched allele of the recipient is homozygous. An antigen in the HVG direction can be a major target for the remaining recipient T cells and can lead to graft rejection (Arima et al., 2018). Interestingly, a study from Japan showed a reduced risk of relapse after HLA matched transplantation in homozygous HLA-C1 patients with acute myeloid leukemia and chronic myeloid leukemia but not with acute lymphoblastic leukemia (Ito et al., 1988).

For these patients with HLA homozygous haplotype, related family donors who are haploidentical with no HLA mismatch in the GVHD direction are willingly available for hematopoietic cell transplantation (HCT). Kanda et al. have recently shown that patients with HLA homozygous haplotypes achieved better transplantation outcomes when the donor is an HLA-matched sibling than with HLA heterozygous haplotypes and when it is needed, transplantation from a haploidentical donor without T cell depletion is a viable option, given the comparable transplant outcomes for heterozygous-to-homozygous HCT and HLA-matched sibling heterozygous HCT (Morishima et al., 2020), suggesting that the Saudi population will have more success rate in HCT haploidentical related donor transplant.

HLA association studies with TA-GVHD have shown that homozygosity was associated with worse clinical outcomes. TA-GVHD in immunocompetent patients occurs mostly when lymphocytes are transfused from an HLA homozygous into an HLA heterozygous patient who shares one haplotype with the donor (Juji et al., 1989; Uchida et al., 1994). It was assumed that recipients are not capable of recognizing donor lymphocytes as foreign because they have similar HLA. This is exactly similar to what happens in mice when T cells are transferred from an inbred homozygous parent to their F1-hybrid offspring (Mitsunaga et al., 1992). Previously, typing HLA before and after the TA-GVHD occurs was so difficult because pancytopenia progresses so rapidly that there are not sufficient lymphocytes for HLA typing. Therefore, with the development of high-throughput HLA typing taking only a few hours using polymerase chain reaction (PCR) amplification technologies, some centers already have started using this procedure (Mitsunaga et al., 1998). The knowledge of the diversity of HLA alleles and haplotypes of the Saudi population can contribute to the understanding of how host genetic factors influence transfusion associated-graft versus host disease (TA-GvHD) risk and transfusion practice in Saudi Arabia. Such knowledge would guide decision-making regarding prophylactic blood irradiation, that is, whether to continue following the current as indicated irradiation policy or it would be safer to shift to the universal irradiation practice.

In conclusion, this work showed a high prevalence of homozygosity of HLA class I and class II alleles and haplotypes in the Saudi population, which reflects the degree of genetic diversity (or similarity) in this population. We believe that these findings would pave the way to better transplantation practices and transfusion-associated graft versus host disease prevention strategies in the Saudi population.

Data Availability Statement

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

Ethics Statement

This retrospective study using participants medical records was reviewed and approved by the KFMC institutional review board and the study was exempt from obtaining informed consent (IRB# 15-198). The study was exempt from obtaining informed consent for each individual as these HSCT donors already consented that their data could be used for research purposes.

Author Contributions

FU (statistical data analysis using Arlequin software and data interpretation); HC ( analyzed data and writing manuscript), KA (experimental design and analyzed data), MP (experimental design and analyzed data), AB (statistical data analysis), GE (data analysis and writing the manuscript), AE (project co-leader, analyzed data, and writing manuscript), AC (project leader, experimental design, analyzed data, and writing manuscript). All authors contributed to the article and approved the submitted version.

Funding

This work was supported by KFMC research grant # 015-019.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors would like to thank the technologists of the Immunology and Pathology Laboratories at KFMC for preparing the patients samples.

Supplementary Material

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

References

Al-Awwami, M., Aljurf, M., Al-Humidan, H., El-Solh, H., and Almeshari, K. (2012). HLA Haplotype Frequencies in saudi arabia for Design of a Saudi Stem-Cell Registry. Hum. Immunol. 73, 106. doi:10.1016/j.humimm.2012.07.211

CrossRef Full Text | Google Scholar

Anasetti, C., Beatty, P. G., Storb, R., Martin, P. J., Mori, M., Sanders, J. E., et al. (1990). Effect of HLA Incompatibility on Graft-Versus-Host Disease, Relapse, and Survival after Marrow Transplantation for Patients with Leukemia or Lymphoma. Hum. Immunol. 29, 79–91. doi:10.1016/0198-8859(90)90071-V

PubMed Abstract | CrossRef Full Text | Google Scholar

Arima, N., Kanda, J., Tanaka, J., Yabe, T., Morishima, Y., Kim, S.-W., et al. (2018). Homozygous HLA-C1 Is Associated with Reduced Risk of Relapse after HLA-Matched Transplantation in Patients with Myeloid Leukemia. Biol. Blood Marrow Transplant. 24, 717–725. doi:10.1016/j.bbmt.2017.11.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Askar, M., Daghstani, J., Thomas, D., Leahy, N., Dunn, P., Claas, F., et al. (2013). 16thIHIW: Global Distribution of Extended HLA Haplotypes. Int. J. Immunogenet. 40, 31–38. doi:10.1111/IJI.12029

PubMed Abstract | CrossRef Full Text | Google Scholar

Azurdia, R. M., Luzzi, G. A., Byren, I., Welsh, K., Wojnarowska, F., Marren, P., et al. (1999). Lichen Sclerosus in Adult Men: A Study of HLA Associations and Susceptibility to Autoimmune Disease. Br. J. Dermatol. 140, 79–83. doi:10.1046/j.1365-2133.1999.02611.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgner, D., Jamieson, S. E., and Blackwell, J. M. (2006). Genetic Susceptibility to Infectious Diseases: Big Is Beautiful, but Will Bigger Be Even Better? Lancet Infect. Dis. 6(10), 653–663. doi:10.1016/S1473-3099(06)70601-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrington, M., Nelson, G. W., Martin, M. P., Kissner, T., Vlahov, D., Goedert, J. J., et al. (1999). HLA and HIV-1: Heterozygote Advantage and B*35 - Cw*04 Disadvantage. Science 283, 1748–1752. doi:10.1126/science.283.5408.1748

PubMed Abstract | CrossRef Full Text | Google Scholar

Creary, L. E., Gangavarapu, S., Mallempati, K. C., Montero-Martín, G., Caillier, S. J., Santaniello, A., et al. (2019). Next-generation Sequencing Reveals New Information about HLA Allele and Haplotype Diversity in a Large European American Population. Hum. Immunol. 80, 807–822. doi:10.1016/j.humimm.2019.07.275

PubMed Abstract | CrossRef Full Text | Google Scholar

Dehn, J., Arora, M., Spellman, S., Setterholm, M., Horowitz, M., Confer, D., et al. (2008). Unrelated Donor Hematopoietic Cell Transplantation: Factors Associated with a Better HLA Match. Biol. Blood Marrow Transplant. 14, 1334–1340. doi:10.1016/j.bbmt.2008.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dendrou, C. A., Petersen, J., Rossjohn, J., and Fugger, L. (2018). HLA Variation and Disease. Nat. Rev. Immunol. 18, 325–339. doi:10.1038/nri.2017.143

PubMed Abstract | CrossRef Full Text | Google Scholar

Doherty, P. C., and Zinkernagel, R. M. (1975). Enhanced Immunological Surveillance in Mice Heterozygous at the H-2 Gene Complex. Nature 256, 50–52. doi:10.1038/256050a0

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Hazmi, M. A., Al-Swailem, A. R., Warsy, A. S., Al-Swailem, A. M., Sulaimani, R., and Al-Meshari, A. A. (1995). Consanguinity Among the Saudi Arabian Population. J. Med. Genet. 32, 623–626. doi:10.1136/JMG.32.8.623

PubMed Abstract | CrossRef Full Text | Google Scholar

Excoffier, L., and Lischer, H. E. L. (2010). Arlequin Suite Ver 3.5: a New Series of Programs to Perform Population Genetics Analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi:10.1111/J.1755-0998.2010.02847.X

PubMed Abstract | CrossRef Full Text | Google Scholar

Fakhoury, H., Alburaq, M., AlSebaiheen, F., Alshubaili, A., and Hajeer, A. H. (2021). Identification of a Novel HLA‐B *18 Variant, HLA‐B *18:01:01:52 , in a Saudi Individual. HLA 97, 359–360. doi:10.1111/tan.14161

PubMed Abstract | CrossRef Full Text | Google Scholar

General Authority of Statistics, (2016). Kingdom of Saudi Arabia. Demography Survey 2021, Demography Survey 2016. Available at: https//Www.Stats.Gov. Sa/Sites/Default/Files/En-Demographic-Research-2016_2.Pdf (Accessed 04 21, 2022).

Google Scholar

Gladman, D. D., Anhorn, K. A. B., Schachter, R. K., and Mervart, H. (1986). HLA Antigens in Psoriatic Arthritis. J. Rheumatol. 13(3), 586–592.

PubMed Abstract | Google Scholar

Hajeer, A. H., Al Balwi, M. A., Aytül Uyar, F., Alhaidan, Y., Alabdulrahman, A., Al Abdulkareem, I., et al. (2013). HLA-A, -B, -C, -DRB1 and -DQB1 Allele and Haplotype Frequencies in Saudis Using Next Generation Sequencing Technique. Tissue Antigens 82, 252–258. doi:10.1111/TAN.12200

PubMed Abstract | CrossRef Full Text | Google Scholar

Hajjej, A., Almawi, W. Y., Arnaiz-Villena, A., Hattab, L., and Hmida, S. (2018). The Genetic Heterogeneity of Arab Populations as Inferred from HLA Genes. PLoS One 13, e0192269. doi:10.1371/journal.pone.0192269

PubMed Abstract | CrossRef Full Text | Google Scholar

Hojjat-Farsangi, M., Jeddi-Tehrani, M., Amirzargar, A. A., Razavi, S. M., Sharifian, R. A., Rabbani, H., et al. (2008). Human Leukocyte Antigen Class II Allele Association to Disease Progression in Iranian Patients with Chronic Lymphocytic Leukemia. Hum. Immunol. 69, 666–674. doi:10.1016/j.humimm.2008.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, K., Yoshida, H., Yanagibashi, K., Shimada, Y., Imamura, M., Tobe, T., et al. (1988). Change of HLA Phenotype in Postoperative Erythroderma. Lancet 331, 413–414. doi:10.1016/S0140-6736(88)91204-4

CrossRef Full Text | Google Scholar

Jawdat, D., Al Saleh, S., Sutton, P., Al Anazi, H., Shubaili, A., Uyar, F. A., et al. (2012). HLA-C Polymorphisms in Two Cohorts of Donors for Bone Marrow Transplantation. Saudi J. Kidney Dis. Transpl. 23, 467–470. doi:10.4103/1319-2442.95705

PubMed Abstract | CrossRef Full Text | Google Scholar

Jawdat, D., Al‐Zahrani, M., Al‐Askar, A., Fakhoury, H., Uyar, F. A., and Hajeer, A. (2019). HLA‐A, B, C, DRB1 and DQB1 Allele and Haplotype Frequencies in Volunteer Bone Marrow Donors from Eastern Region of Saudi Arabia. HLA 94, 49–56. doi:10.1111/TAN.13533

PubMed Abstract | CrossRef Full Text | Google Scholar

Jawdat, D. M., Al Saleh, S., Sutton, P., Al Anazi, H., Shubaili, A., Tamim, H., et al. (2009). Chances of Finding an HLA-Matched Sibling: The Saudi Experience. Biol. Blood Marrow Transplant. 15, 1342–1344. doi:10.1016/j.bbmt.2009.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Jawdat, D., Uyar, F. A., Alaskar, A., Müller, C. R., and Hajeer, A. (2020). HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1 Allele and Haplotype Frequencies of 28,927 Saudi Stem Cell Donors Typed by Next-Generation Sequencing. Front. Immunol. 11. doi:10.3389/fimmu.2020.544768

PubMed Abstract | CrossRef Full Text | Google Scholar

Juji, T., Takahashi, K., Shibata, Y., Ide, H., Sakakibara, T., Ino, T., et al. (1989). Post-transfusion Graft-Versus-Host Disease in Immunocompetent Patients after Cardiac Surgery in Japan. N. Engl. J. Med. 321, 56. doi:10.1056/nejm198907063210116

CrossRef Full Text | Google Scholar

Kanda, J., Ichinohe, T., Fuji, S., Maeda, Y., Ohashi, K., Fukuda, T., et al. (2015). Impact of HLA Mismatch Direction on the Outcome of Unrelated Bone Marrow Transplantation: A Retrospective Analysis from the Japan Society for Hematopoietic Cell Transplantation. Biol. Blood Marrow Transplant. 21, 305–311. doi:10.1016/j.bbmt.2014.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanda, J., Ikegame, K., Fuji, S., Kurokawa, M., Kanamori, H., Fukuda, T., et al. (2016). Haploidentical and Matched Sibling Donor Hematopoietic Cell Transplantation for Patients with HLA-Homozygous Haplotypes. Biol. Blood Marrow Transplant. 22, 2031–2037. doi:10.1016/j.bbmt.2016.07.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawashima, M., Ohashi, J., Nishida, N., and Tokunaga, K. (2012). Evolutionary Analysis of Classical HLA Class I and II Genes Suggests that Recent Positive Selection Acted on DPB1*04∶01 in Japanese Population. PLoS One 7, e46806. doi:10.1371/JOURNAL.PONE.0046806

PubMed Abstract | CrossRef Full Text | Google Scholar

Kunze-Schumacher, H., Blasczyk, R., and Bade-Doeding, C. (2014). Soluble HLA Technology as a Strategy to Evaluate the Impact of HLA Mismatches. J. Immunol. Res. 2014, 1–8. doi:10.1155/2014/246171

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. J., Klein, J., Haagenson, M., Baxter-Lowe, L. A., Confer, D. L., Eapen, M., et al. (2007). High-resolution Donor-Recipient HLA Matching Contributes to the Success of Unrelated Donor Marrow Transplantation. Blood 110, 4576–4583. doi:10.1182/blood-2007-06-097386

PubMed Abstract | CrossRef Full Text | Google Scholar

Maiers, M., Gragert, L., and Klitz, W. (2007). High-resolution HLA Alleles and Haplotypes in the United States Population. Hum. Immunol. 68, 779–788. doi:10.1016/j.humimm.2007.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Marrack, P., and Kappler, J. (1986). The Antigen-specific, Major Histocompatibility Complex-Restricted Receptor on T Cells. Adv. Immunol. 38, 1–30. doi:10.1016/S0065-2776(08)60005-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitsunaga, S., Kuwata, S., Tokunaga, K., Uchikawa, C., Takahashi, K., Akaza, T., et al. (1992). Family Study on HLA-DPB1 Polymorphism: Linkage Analysis with HLA-DR/DQ and Two "new" Alleles. Hum. Immunol. 34, 203–211. doi:10.1016/0198-8859(92)90113-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitsunaga, S., Tokunaga, K., Kashiwase, K., Akaza, T., Tadokoro, K., and Juji, T. (1998). A Nested PCR-RFLP Method for High-Resolution Typing of HLA-A Alleles. Eur. J. Immunogenet. 25, 15–27. doi:10.1046/j.1365-2370.1998.00093.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Morishima, Y., Morishima, S., Murata, M., Arima, N., Uchida, N., Sugio, Y., et al. (2020). Impact of Homozygous Conserved Extended HLA Haplotype on Single Cord Blood Transplantation: Lessons for Induced Pluripotent Stem Cell Banking and Transplantation in Allogeneic Settings. Biol. Blood Marrow Transplant. 26, 132–138. doi:10.1016/j.bbmt.2019.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Ollier, W., Doyle, P., Alonso, A., Awad, J., Williams, E., Gill, D., et al. (1985). HLA Polymorphisms in Saudi Arabs. Tissue Antigens 25, 87–95. doi:10.1111/j.1399-0039.1985.tb00420.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Osman, A. E., Gonzalez-Galarza, F., Mubasher, M., Al-Harthi, H., El Sheikh, N., Berka, N., et al. (2014). HLA-A, -B, -C, -DRB1, and -DQB1 Allele Lineages and Haplotype Frequencies Among Saudis. Iii 6, 1–6. doi:10.4137/III.S16796

CrossRef Full Text | Google Scholar

Penn, D. J., Damjanovich, K., and Potts, W. K. (2002). MHC Heterozygosity Confers a Selective Advantage against Multiple-Strain Infections. Proc. Natl. Acad. Sci. U.S.A. 99, 11260–11264. doi:10.1073/pnas.162006499

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersdorf, E. W. (2008). Optimal HLA Matching in Hematopoietic Cell Transplantation. Curr. Opin. Immunol. 20, 588–593. doi:10.1016/j.coi.2008.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Peterson, T. A., Bielawny, T., Lacap, P., Hardie, R.-A., Daniuk, C., Mendoza, L., et al. (2014). Diversity and Frequencies of HLA Class I and Class II Genes of an East African Population. OJGen 04, 99–124. doi:10.4236/ojgen.2014.42013

CrossRef Full Text | Google Scholar

Robinson, J., Barker, D. J., Georgiou, X., Cooper, M. A., Flicek, P., and Marsh, S. G. E. (2020). IPD-IMGT/HLA Database. Nucleic Acids Res. 48(1), 948–955. doi:10.1093/nar/gkz950

PubMed Abstract | CrossRef Full Text | Google Scholar

Rohowsky-Kochan, C., Skurnick, J., Molinaro, D., and Louria, D. (1998). HLA Antigens Associated with Susceptibility/resistance to HIV-1 Infection. Hum. Immunol. 59(12), 802–815. doi:10.1016/S0198-8859(98)00086-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Saw, W.-Y., Liu, X., Liu, X., Khor, C.-C., Takeuchi, F., Katsuya, T., et al. (2015). Mapping the Genetic Diversity of HLA Haplotypes in the Japanese Populations. Sci. Rep. 5. doi:10.1038/srep17855

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, N., Decker, W. K., Lapushin, R., Xing, D., Robinson, S. N., Yang, H., et al. (2011). HLA Homozygosity and Haplotype Bias Among Patients with Chronic Lymphocytic Leukemia: Implications for Disease Control by Physiological Immune Surveillance. Leukemia 25, 1036–1039. doi:10.1038/leu.2011.30

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheth, K. V., Edwards, J. A., and Godwin, J. T. (1985). Study of the HLA Gene and Antigen Frequency from a Saudi Arabian Hospital. Tissue Antigens 25, 156–162. doi:10.1111/j.1399-0039.1985.tb00430.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shiina, T., Inoko, H., and Kulski, J. K. (2004). An Update of the HLA Genomic Region, Locus Information and Disease Associations: 2004. Tissue Antigens 64, 631–649. doi:10.1111/j.1399-0039.2004.00327.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Thursz, M. R., Thomas, H. C., Greenwood, B. M., and Hill, A. V. S. (1997). Heterozygote Advantage for HLA Class-II Type in Hepatitis B Virus Infection. Nat. Genet. 17, 11–12. doi:10.1038/ng0997-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Treleaven, J., Gennery, A., Marsh, J., Norfolk, D., Page, L., Parker, A., et al. (2011). Guidelines on the Use of Irradiated Blood Components Prepared by the British Committee for Standards in Haematology Blood Transfusion Task Force. Br. J. Haematol. 152, 35–51. doi:10.1111/j.1365-2141.2010.08444.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Trowsdale, J., and Knight, J. C. (2013). Major Histocompatibility Complex Genomics and Human Disease. Annu. Rev. Genom. Hum. Genet. 14, 301–323. doi:10.1146/annurev-genom-091212-153455

PubMed Abstract | CrossRef Full Text | Google Scholar

Uchida, S., Suzuki, K., Akiyama, S., Miyamoto, M., Juji, T., and Fujiwara, M. (1994). Suppressive Effect of Cyclophosphamide on the Progression of Lethal Graft-Versus-Host Disease in Mice-Aa Therapeutic Model of Fatal Post-transfusion GVHD. Ther. Immunol. 1, 313–318.

PubMed Abstract | Google Scholar

Warsy, A., Al-Jaser, M., Albdass, A., Al-Daihan, S., and Alanazi, M. (2014). Is Consanguinity Prevalence Decreasing in Saudis?: A Study in Two Generations. Afr. H. Sci. 14, 314–321. doi:10.4314/AHS.V14I2.5

CrossRef Full Text | Google Scholar

Keywords: HLA, genetic diversity, allele frequency, allele and haplotype homozygosity, TA-GVHD

Citation: Chentoufi AA, Uyar FA, Chentoufi HA, Alzahrani K, Paz M, Bahnassy A, Elyamany G and Elghazaly A (2022) HLA Diversity in Saudi Population: High Frequency of Homozygous HLA Alleles and Haplotypes. Front. Genet. 13:898235. doi: 10.3389/fgene.2022.898235

Received: 17 March 2022; Accepted: 25 April 2022;
Published: 09 June 2022.

Edited by:

Divyansh Agarwal, Massachusetts General Hospital and Harvard Medical School, United States

Reviewed by:

Xiuwen Zheng, AbbVie, United States
Seik-Soon Khor, National Center For Global Health and Medicine, Japan

Copyright © 2022 Chentoufi, Uyar, Chentoufi, Alzahrani, Paz, Bahnassy, Elyamany and Elghazaly. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Aziz Alami Chentoufi, Y2hlbnRvdWZpYXppekBob3RtYWlsLmNvbQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.