HLA-A, -B, -C, -DRB1 and -DQB1 allele and haplotype frequencies in Saudis using next generation sequencing technique

High-resolution KIR and HLA genotyping in three Chinese ethnic minorities reveals distinct origins

Polymorphism of killer-cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands impacts the effector activity of cytotoxic NK cell and T cell subsets. Therefore, understanding the extent and implications of KIR and HLA class I genetic polymorphism across various populations is important for immunological and medical research. In this study, we conducted a high-resolution investigation of KIR and HLA class I diversity in three distinct Chinese ethnic minority populations. We studied the She, Yugur, and Tajik, and compared them with the Zhejiang Han population (Zhe), which represents the majority Southern Han ethnicity. Our findings revealed that the Tajik population exhibited the most diverse KIR copy number, allele, and haplotype diversity among the four populations. This diversity aligns with their proposed ancestral origin, closely resembling that of Iranian populations, with a relatively higher presence of KIR-B genes, alleles, and haplotypes compared with the other Chinese populations. The Yugur population displayed KIR distributions similar to those of the Tibetans and Southeast Asians, whereas the She population resembled the Zhe and other East Asians, as confirmed by genetic distance analysis of KIR. Additionally, we identified 12.9% of individuals across the three minority populations as having KIR haplotypes characterized by specific gene block insertions or deletions. Genetic analysis based on HLA alleles yielded consistent results, even though there were extensive variations in HLA alleles. The observed variations in KIR interactions, such as higher numbers of 2DL1-C2 interactions in Tajik and Yugur populations and of 2DL3-C1 interactions in the She population, are likely shaped by demographic and evolutionary mechanisms specific to their local environments. Overall, our findings offer valuable insights into the distribution of KIR and HLA diversity among three distinct Chinese ethnic minority populations, which can inform future clinical and population studies.

Race, ethnicity, ancestry, and aspects that impact HLA data and matching for transplant

Race, ethnicity, and ancestry are terms that are often misinterpreted and/or used interchangeably. There is lack of consensus in the scientific literature on the definition of these terms and insufficient guidelines on the proper classification, collection, and application of this data in the scientific community. However, defining groups for human populations is crucial for multiple healthcare applications and clinical research. Some examples impacted by population classification include HLA matching for stem-cell or solid organ transplant, identifying disease associations and/or adverse drug reactions, defining social determinants of health, understanding diverse representation in research studies, and identifying potential biases. This article describes aspects of race, ethnicity and ancestry information that impact the stem-cell or solid organ transplantation field with particular focus on HLA data collected from donors and recipients by donor registries or transplant centers.

HLA pharmacogenetic markers of drug hypersensitivity from the perspective of the populations of the Greater Middle East

Specific HLA associations with drug hypersensitivity may vary between geographic regions and ethnic groups. There are little to no data related to HLA-drug hypersensitivity on populations who reside in the Greater Middle East (GME), a vast region spanning from Morocco in the west to Pakistan in the east. In this review, the authors intended to summarize the significant HLA alleles associated with hypersensitive drug reactions induced by different drugs, as have been found in different populations, and to summarize the prevalence of these alleles in the specific and diverse populations of the GME. For example, HLA-B*57:01 allele prevalence, associated with abacavir-induced hypersensitivity, ranges from 1% to 3%, and HLA-DPB1*03:01 prevalence, associated with aspirin-induced asthma, ranges from 10% to 14% in the GME population. Studying pharmacogenomic associations in the ethnic groups of the GME may allow the discovery of new associations, confirm ones found with a low evidence rate and enable cost-effectiveness analysis of allele screening before drug use.

Distribution of HLA-B Alleles and Haplotypes in Qatari: Recommendation for Establishing Pharmacogenomic Markers Screening for Drug Hypersensitivity

Human leukocyte antigen (HLA) proteins are present at the cellular surface of antigen-presenting cells and play a crucial role in the adaptive immune response. Class I genes, specifically certain HLA-B alleles, are associated with adverse drug reactions (ADRs) and are used as pharmacogenetic markers. Although ADRs are a common causes of hospitalization and mortality, the data on the prevalence of HLA-B pharmacogenetics markers in Arab countries are scarce. In this study, we investigated the frequencies of major HLA-B pharmacogenomics markers in the Qatari population. Next-generation sequencing data from 1,098 Qatari individuals were employed for HLA-B typing using HLA-HD version 1.4.0 and IPD-IMGT/HLA database. In addition, HLA-B pharmacogenetics markers were obtained from the HLA Adverse Drug Reaction Database. In total, 469 major HLA-B pharmacogenetic markers were identified, with HLA-B*51:01 being the most frequent pharmacogenetic marker (26.67%) in the Qatari population. Moreover, HLA-B*51:01 is associated with phenytoin- and clindamycin-induced ADRs. The second most frequent pharmacogenetic marker was the HLA-B*58:01 allele (6.56%), which is associated with allopurinol-induced ADRs. The third most frequent pharmacogenetic marker was the HLA-B*44:03 allele, which is associated with phenytoin-induced ADRs. The establishment of a pharmacogenetics screening program in Qatar for cost effective interventions aimed at preventing drug-induced hypersensitivity can be aided by the highly prevalent HLA-B pharmacogenetic markers detected here.

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

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.

Identification of 8-Digit HLA-A, -B, -C, and -DRB1 Allele and Haplotype Frequencies in Koreans Using the One Lambda AllType Next-Generation Sequencing Kit

Background

Recent studies have successfully implemented next-generation sequencing (NGS) in HLA typing. We performed HLA NGS in a Korean population to estimate HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies up to an 8-digit resolution, which might be useful for an extended application of HLA results.

Methods

A total of 128 samples collected from healthy unrelated Korean adults, previously subjected to Sanger sequencing for 6-digit HLA analysis, were used. NGS was performed for HLA-A, -B, -C, and -DRB1 using the AllType NGS kit (One Lambda, West Hills, CA, USA), Ion Torrent S5 platform (Thermo Fisher Scientific, Waltham, MA, USA), and Type Steam Visual NGS analysis software (One Lambda).

Results

Eight HLA alleles showed frequencies of ≥10% in the Korean population, namely, A*24:02:01:01 (19.5%), A*33:03:01 (15.6%), A*02:01:01:01 (14.5%), A*11:01:01:01 (13.3%), B*15:01:01:01 (10.2%), C*01:02:01 (19.9%), C*03:04:01:02 (11.3%), and DRB1*09:01:02 (10.2%). Nine previous 6-digit HLA alleles were further identified as two or more 8-digit HLA alleles. Of these, eight alleles (A*24:02:01, B*35:01:01, B*40:01:02, B*55:02:01, B*58:01:01, C*03:02:02, C*07:02:01, and DRB1*07:01:01) were identified as two 8-digit HLA alleles, and one allele (B*51:01:01) was identified as three 8-digit HLA alleles. The most frequent four-loci haplotype was HLA-A*33:03:01-B*44:03:01:01-C*14:03-DRB1*13:02:01.

Conclusions

We identified 8-digit HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies in a healthy Korean population using NGS. These new data can be used as a representative Korean data for further disease-related HLA type analysis.

Association between HLA genotype and ferritin levels in COVID-19 infection: a study of a Saudi cohort

BACKGROUND

An elevation in serum ferritin levels through an unknown mechanism was observed in COVID-19 infected patients. This study examined the association between patients' HLA genotype and serum ferritin level modulation and also assessed the effect of serum ferritin levels on infection severity/mortality.

MATERIALS AND METHODS

One hundred and thirty-six COVID-19 Saudi patients were divided into two groups according to their ferritin levels: group 1 (<500 ng/mL, N = 67) and group 2 (>500 ng/mL, N = 69). HLA genotyping (class I and II) was carried out through the rPCR-SSO method.

RESULTS

High serum ferritin levels were associated with a significant increase in infection severity, as 75% of ICU patients showed high levels of ferritin compared to 43% of patients with moderate symptoms, p = .002. This elevation indicated a gender skew in that 56% of the infected male patients displayed high ferritin levels compared to 36.6% of the female patients, p = .03. In terms of mortality, 74% of patients with fatal outcomes had a high level of serum ferritin compared to 47% of recovered patients, p = .039. There was a significant difference in the HLA frequency between the two groups, with a predominance of HLA-A*01 in the low-ferritin group (19.4 vs. 6.5%, p = .002, p_c = .016) and predominance of C*03 in the high ferritin group (10.9 vs. 3%, p = .047, p_c = .27).

CONCLUSION

High serum ferritin levels are associated with an increase in COVID-19 severity, which may be affected by HLA polymorphism.

Segregation Analysis of Genotyped and Family-Phased, Long Range MHC Classical Class I and Class II Haplotypes in 5 Families With Type 1 Diabetes Proband in the United Arab Emirates

The classical Human Leucocyte Antigen (HLA) class II haplotypes of the Major Histocompatibility Complex (MHC) that are associated with type 1 diabetes (T1D) were identified in five families from the United Arab Emirates (UAE). Segregation analyses were performed on these 5 families with the disease, 3 with one child and 2 with 2 children diagnosed with T1D. Three HLA-DR4 haplotypes were identified: HLA- DRB1^∗04:01:01-DQB1^∗03:02:01:01; HLA- DRB1^∗04:02:01- DQB1^∗03:02:01; and HLA -DRB1^∗04:05:01-DQB1^∗02:02:01:02. All have previously been identified to be associated with T1D in studies of the Arabian population. In the 10 parents from the 5 families, 9 had at least one HLA-DR4 and HLA-DR3 haplotype which potentially increases the risk of T1D. Of these 9 parents, 3 were heterozygous for HLA-DR4/HLA-DR3 and one was homozygous for HLA-DR3. Two haplotypes that were identified here extend to the HLA class I region were previously designated AH8.2 (HLA -A^∗26-B^∗08-DRB1^∗03) and AH50.2 (HLA -C^∗06-B^∗50-DRB1^∗03:01-DQ^∗02) and associated with diabetes in neighboring North Indian populations. This study provides examples of MHC haplotype analysis in pedigrees to improve our understanding of the genetics of T1D in the understudied population of the UAE.

A Novel Allele-Specific PCR Protocol for the Detection of the HLA-C*03:02 Allele, a Pharmacogenetic Marker, in Vietnamese Kinh People

Background

Allopurinol, a common anti-hyperuricemia drug, is well known as an inducer of severe cutaneous adverse drug reactions (SCARs). One of the most well-defined risk factors of allopurinol-induced SCARs is the presence of polymorphic alleles of human leukocyte antigen (HLA) genes, such as HLA-B*58:01 and HLA-C*03:02 alleles. There is no commercial test or published in-house protocol for the specific detection of the HLA-C*03:02 allele. In this article, we established for the first time a simple allele-specific (AS) PCR method to identify HLA-C*03:02 allele carriers, and at the same time, determine their zygosities.

Methods

A two-step AS-PCR protocol, using four primer sets, was designed to specifically amplify and differentiate the HLA-C*03:02 allele from 17 other HLA-C alleles found in Vietnamese people. The protocol was validated with PCR-sequencing-based typing (SBT) of 100 samples of unknown genotypes.

Results

The PCR protocol can detect the HLA-C*03:02 allele and determine the zygosity. The results of this protocol were highly consistent with those of the SBT (ĸ = 0.98, p < 0.001). Regarding the specific detection of the HLA-C*03:02 allele, the PCR protocol had a sensitivity of 100% (95% CI: 91.61-100%) and specificity of 98.3% (95% CI: 90.9-99.7%). The protocol was used to determine the distribution of the HLA-C*03:02 allele in 810 unrelated Vietnamese Kinh people, 14.2% of which were HLA-C*03:02 carriers, the allele frequency was 7.5%.

Conclusion

A novel AS-PCR protocol with a sensitivity of 100% for the detection of the HLA-C*03:02 allele was established. The protocol can be used for personalized treatment with allopurinol in order to minimize the risk of SCARs in Vietnamese people as well as in other Asian populations with similar genetic characteristics.

High-resolution HLA allele and haplotype frequencies of the Saudi Arabian population based on 45,457 individuals and corresponding stem cell donor matching probabilities

We estimated HLA allele and haplotype frequencies of the Saudi Arabian population from a sample of 45,457 registered stem cell donors. The most frequent HLA alleles were A*02:01g (18.5%), C*06:02g (16.1%), B*51:01g (14.1%), DRB1*07:01g (16.2%), DQB1*02:01g (30.5%), and DPB1*04:01g (33.6%). The most frequent 5-locus haplotypes were A*02:05g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.73%), A*02:01g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.66%), and A*26:01g~C*07:02g~B*08:01g~DRB1*03:01g~DQB1*02:01g (1.38%). Furthermore, we used the calculated haplotype frequencies to estimate stem cell donor matching probabilities for Saudi Arabian donor and patient populations under various matching requirements. These results are relevant for strategic donor registry planning in the Kingdom of Saudi Arabia.

Association of HLA-DR-DQ alleles, haplotypes, and diplotypes with type 1 diabetes in Saudis

AIMS

Type 1 diabetes (T1D) is an autoimmune disease that affects many children worldwide. Genetic factors and environmental triggers play crucial interacting roles in the aetiology. This study aimed to assess the contribution of HLA-DRB1-DQA1-DQB1 alleles, haplotypes, and genotypes to the risk of T1D among Saudis.

METHODS

A total of 222 children with T1D and 342 controls were genotyped for HLA-DRB1, -DQA1, and -DQB1 using reverse sequence-specific oligonucleotide (rSSO) Lab Type high definition (HD) kits. Alleles, haplotypes, and diplotypes were compared between cases and controls using the SAS statistical package.

RESULTS

DRB1*03:01-DQA1*05:01-DQB1*02:01 (32.4%; OR = 3.68; P_c < .0001), DRB1*04:05-DQA1*03:02-DQB1*03:02 (6.6%; OR = 6.76; P_c < .0001), DRB1*04:02-DQA1*03:01-DQB1*03:02 (6.0%; OR = 3.10; P_c = .0194), DRB1*04:01-DQA1*03:01-DQB1*03:02 (3.7%; OR = 4.22; P_c = .0335), and DRB1*04:05-DQA1*03:02-DQB1*02:02 (2.7%; OR = 6.31; P_c = .0326) haplotypes were significantly increased in cases compared to controls, whereas DRB1*07:01-DQA1*02:01-DQB1*02:02 (OR = 0.41; P_c = .0001), DRB1*13:01-DQA1*01:03-DQB1*06:03 (OR = 0.05; P_c < .0001), DRB1*15:01-DQA1*01:02-DQB1*06:02 (OR = 0.03; P_c < .0001), and DRB1*11:01-DQA1*05:05-DQB1*03:01 (OR = 0.07; P_c = .0291) were significantly decreased. Homozygous DRB1*03:01-DQA1*05:01-DQB1*02:01 genotypes and combinations of DRB1*03:01-DQA1*05:01-DQB1*02:01 with DRB1*04:05-DQA1*03:02-DQB1*03:02, DRB1*04:02-DQA1*03:01-DQB1*03:02, and DRB1*04:01-DQA1*03:01-DQB1*03:02 were significantly increased in cases than controls. Combinations of DRB1*03:01-DQA1*05:01-DQB1*02:01 with DRB1*07:01-DQA1*02:01-DQB1*02:02 and DRB1*13:02-DQA1*01:02-DQB1*06:04 showed low OR values but did not remain significantly decreased after Bonferroni correction.

CONCLUSIONS

HLA-DRB1-DQA1-DQB1 alleles, haplotypes, and diplotypes in Saudis with T1D are not markedly different from those observed in Western and Middle-Eastern populations but are quite different than those of East Asians.

Diversity of HLA Class I and Class II blocks and conserved extended haplotypes in Lacandon Mayans

Here we studied HLA blocks and haplotypes in a group of 218 Lacandon Maya Native American using a high-resolution next generation sequencing (NGS) method. We assessed the genetic diversity of HLA class I and class II in this population, and determined the most probable ancestry of Lacandon Maya HLA class I and class II haplotypes. Importantly, this Native American group showed a high degree of both HLA homozygosity and linkage disequilibrium across the HLA region and also lower class II HLA allelic diversity than most previously reported populations (including other Native American groups). Distinctive alleles present in the Lacandon population include HLA-A*24:14 and HLA-B*40:08. Furthermore, in Lacandons we observed a high frequency of haplotypes containing the allele HLA-DRB1*04:11, a relatively frequent allele in comparison with other neighboring indigenous groups. The specific demographic history of the Lacandon population including inbreeding, as well as pathogen selection, may have elevated the frequencies of a small number of HLA class II alleles and DNA blocks. To assess the possible role of different selective pressures in determining Native American HLA diversity, we evaluated the relationship between genetic diversity at HLA-A, HLA-B and HLA-DRB1 and pathogen richness for a global dataset and for Native American populations alone. In keeping with previous studies of such relationships we included distance from Africa as a covariate. After correction for multiple comparisons we did not find any significant relationship between pathogen diversity and HLA genetic diversity (as measured by polymorphism information content) in either our global dataset or the Native American subset of the dataset. We found the expected negative relationship between genetic diversity and distance from Africa in the global dataset, but no relationship between HLA genetic diversity and distance from Africa when Native American populations were considered alone.

HLA Haplotype Frequencies and Genetic Profiles of the Kuwaiti Population

OBJECTIVE

The aim of this study was to assess the HLA haplotype frequencies and genetic profiles of the Kuwaiti population.

MATERIALS AND METHODS

Whole venous blood was obtained from 595 healthy, unrelated Kuwaiti volunteers. The study population was genotyped for HLA class I (HLA-A, HLA-B, and HLA-C) and class II (HLA-DRB1 and HLA-DQB1) loci using sequence-specific oligonucleotide (SSO) probe-based hybridization and high-resolution HLA genotyping. Haplotype frequencies were estimated using an implementation of the expectation maximization algorithm that resolves both phase and allelic ambiguity. The Kuwaiti population was compared with other populations from the US National Marrow Donor Program (NMDP), by running a principal component analysis (PCA) on the relevant haplotype frequencies.

RESULTS

The most common HLA class I alleles in Kuwait were HLA-A*02:01g, HLA-C*06:02g, and HLA-B*50:01g with frequencies of 16, 14, and 12%, respectively. The most common HLA class II alleles in Kuwait were HLA-DQB1*02:01g and HLA-DRB1*07:01 with frequencies of 29.7 and 16.5%, respectively. The most common Kuwaiti haplotype observed was HLA-A*02:01g∼HLA-C*06:02g∼HLA-B*50:01g∼HLA-DRB1*07:01∼HLA-DQB1*02:01g at a frequency of 2.3%. The PCA demonstrated close genetic proximity of the Kuwaiti population with Middle Eastern, Southeast Asian, and North African populations in the NMDP.

CONCLUSION

Identifying the haplotype diversity in the Kuwaiti population will contribute to the selection of an HLA-match for HSCT, disease associations, pharmacogenomics, and knowledge of pop-ulation HLA diversity.

Diagnostic Utility of Human Leukocyte Antigen B*15:02 Screening in Severe Carbamazepine Hypersensitivity Syndrome

Background

Despite many studies suggesting an association between human leukocyte antigen (HLA)-B*15:02 and carbamazepine (CBZ)-induced severe cutaneous adverse drug reactions essentially toxic epidermal necrolysis (TEN) and Stevens-Johnson syndrome (SJS), the evidence of association in different populations and the degree of association remain uncertain.

Materials and Methods

The primary analysis was based on population control studies. Data were pooled by means of a random-effects model, and sensitivity, specificity, positive and negative likelihood ratios (LR+ and LR-), diagnostic odds ratios (DOR), and areas under the summary receiver operating characteristic curve (AUC) were calculated.

Results

In 23 population control studies, HLA-B*15:02 was measured in 373 patients with CBZ-induced TEN/SJS and 3452 patients without CBZ-induced TEN/SJS. The pooled sensitivity, specificity, LR+, LR-, DOR, and AUC were 0.67 (95% confidence interval [CI] = 0.63-0.72), 0.98 (95% CI = 0.98-0.99), 19.73 (95% CI = 10.54-36.92), 0.34 (95% CI = 0.23-0.49), 71.38 (95% CI = 34.89-146.05), and 0.96 (95% CI = 0.92-0.98), respectively. Subgroup analyses for Han Chinese, Thai, and Malaysian populations yielded similar findings. Specifically, racial/ethnic subgroup analyses revealed similar findings with respect to DOR for Han Chinese (99.28; 95% CI = 22.20-443.88), Thai (61.01; 95% CI = 23.05-161.44), and Malaysian (30; 95% CI = 7.08-126.68) populations, which are similar to the pooled DOR for the relationship between the HLA-B*15:02 allele and CBZ-induced TEN/SJS across all populations (71.38; 95% CI = 34.89-146.05).

Conclusions

The present study reveals that CBZ is the leading cause of TEN/SJS in many countries. Screening of HLA-B*15:02 may help patients to prevent the occurrence of CBZ-induced TEN/SJS, especially in populations with a higher (≥5%) risk allele frequency.

HLA-A, B, C, DRB1 and DQB1 allele and haplotype frequencies in volunteer bone marrow donors from Eastern Region of Saudi Arabia

We analyzed HLA allele and haplotype frequencies from donors from the Eastern region of Saudi Arabia. A cross-sectional study was performed on 2405 bone marrow donors from the Eastern region. HLA typing was carried out by sequencing. The most common HLA allele groups were HLA-A*02:01:01G (11.08%), A*01:01:01G (10.40%), HLA-B*52:01:01G (8.79%), B*18:01:01G (8.07%), HLA-C*04:01:01G (17.88%), C*12:03:01G (10.23%), HLA-DRB1*10:01 (14.89%), DRB1*03:01:01G (14.10%), HLA-DQB1*02:01:01G (24.53%) and DQB1*05:01:01G (20.17%). The most frequent HLA-A^~ C^~ B^~ DRB1^~ DQB1 haplotypes were HLA-A*01:03^~ C*15:05:01G^~ B*73:01^~ DRB1*10:01:01^~ DQB1*05:01:01G (3.11%) and HLA-A*01:01:01G^~ C*12:02:01G^~ B*52:01:01G^~ DRB1*15:02:01^~ DQB1*06:01:01G (2.25%). When comparing the allele and haplotype frequencies of the Eastern regions' population to those from the Central region we found significant differences in several allele frequencies including A*01:01:01G (P ≤ 0.0001), B*52:01:01G (P ≤ 0.0001), B*18:01:01G (P = 0.0001), C*12:03:01G (P < 0.0001), DRB1*10:01:01 (P < 0.0001) and DQB1*05:01:01G (P < 0.0001). Our data confirms genetic heterogeneity among the Saudi population.

HLA class I allele lineages and haplotype frequencies in Arabs of the United Arab Emirates

The high degree of polymorphism of the HLA system provides suitable genetic markers to study the diversity and migration of different world populations and is beneficial for forensic identification, anthropology, transplantation and disease associations. Although the United Arab Emirates (UAE) population of about nine million people is heterogeneous, information is limited for the HLA class I allele and haplotype frequencies of the Bedouin ethnic group. We performed low-resolution PCR-SSP genotyping of three HLA class I loci at HLA-A, -B and -C for 95 unrelated healthy Bedouins from the cities of Al Ain and Abu Dhabi in the UAE. A total of 54 HLA allele lineages were detected; the most frequent low-resolution allele lineages at each HLA locus were A*02 (0.268), B*51 (0.163) and C*07 (0.216). The inferred estimates for the two most frequent HLA-A and HLA-B haplotypes were HLA-A*02 ~ HLA-B*50 (0.070) and HLA-A*02 ~ HLA-B*51 (0.051), and the most frequent 3-locus haplotype was HLA-A*02 ~ HLA-B*50 ~ HLA-C*06 (0.068). The HLA allele lineage frequencies of the UAE Arabs were compared to those previously reported for 70 other world populations, and a strong genetic similarity was detected between the UAE Arabs and the Saudi Arabians from the west with evidence of a limited gene flow between the UAE Arabs and Pakistani across the Gulf from the east, and the UAE Arabs and Omani from the south of the Gulf Peninsula.

Next-Generation Sequencing Based HLA Typing: Deciphering Immunogenetic Aspects of Sarcoidosis

Unraveling of the HLA-related immunogenetic basis of several immune disorders is complex due to the extensive HLA polymorphism and strong linkage-disequilibrium between HLA loci. A lack of in phase sequence information, a relative deficiency of high resolution genotyping including non-coding regions and ambiguous haplotype assignment make it difficult to compare findings across association studies and to attribute a causal role to specific HLA alleles/haplotypes in disease susceptibility and modification of disease phenotypes. Earlier, historical antibody and DNA-based methods of HLA typing, primarily of low resolution at antigen/alellic group levels, yielded "indicative" findings which were partially improved by high-resolution DNA-based typing. Only recently, next-generation sequencing (NGS) approaches based on deep-sequencing of the complete HLA genes combined with bioinformatics tools began to provide the access to complete information at an allelic level. Analyzing HLA with NGS approaches, therefore, promises to provide further insight in the etiopathogenesis of several immune disorders in which HLA associations have been implicated. These range from coeliac disease and rheumatological conditions to even more complex disorders, such as type-1 diabetes, systemic lupus erythematosus and sarcoidosis. A systemic disease of unknown etiology, sarcoidosis has previously been associated with numerous HLA variants and also other gene polymorphisms, often in linkage with the HLA region. To date, the biological significance of these associations has only partially been defined. Therefore, more precise assignments of HLA alleles/haplotypes using NGS approaches could help to elucidate the exact role of HLA variation in the multifaceted etiopathogenesis of sarcoidosis, including epigenetic mechanisms. NGS-based HLA analyses may be also relevant for defining variable clinical phenotypes and for predicting the disease course or the response to current/plausible novel therapies.

The genetic heterogeneity of Arab populations as inferred from HLA genes

This is the first genetic anthropology study on Arabs in MENA (Middle East and North Africa) region. The present meta-analysis included 100 populations from 36 Arab and non-Arab communities, comprising 16,006 individuals, and evaluates the genetic profile of Arabs using HLA class I (A, B) and class II (DRB1, DQB1) genes. A total of 56 Arab populations comprising 10,283 individuals were selected from several databases, and were compared with 44 Mediterranean, Asian, and sub-Saharan populations. The most frequent alleles in Arabs are A*01, A*02, B*35, B*51, DRB1*03:01, DRB1*07:01, DQB1*02:01, and DQB1*03:01, while DRB1*03:01-DQB1*02:01 and DRB1*07:01-DQB1*02:02 are the most frequent class II haplotypes. Dendrograms, correspondence analyses, genetic distances, and haplotype analysis indicate that Arabs could be stratified into four groups. The first consists of North Africans (Algerians, Tunisians, Moroccans, and Libyans), and the first Arabian Peninsula cluster (Saudis, Kuwaitis, and Yemenis), who appear to be related to Western Mediterraneans, including Iberians; this might be explained for a massive migration into these areas when Sahara underwent a relatively rapid desiccation, starting about 10,000 years BC. The second includes Levantine Arabs (Palestinians, Jordanians, Lebanese, and Syrians), along with Iraqi and Egyptians, who are related to Eastern Mediterraneans. The third comprises Sudanese and Comorians, who tend to cluster with Sub-Saharans. The fourth comprises the second Arabian Peninsula cluster, made up of Omanis, Emiratis, and Bahrainis. It is noteworthy that the two large minorities (Berbers and Kurds) are indigenous (autochthonous), and are not genetically different from "host" and neighboring populations. In conclusion, this study confirmed high genetic heterogeneity among present-day Arabs, and especially those of the Arabian Peninsula.

HLA class II polymorphism in Saudi patients with multiple sclerosis

Several studies have investigated the association of different HLA antigens with multiple sclerosis (MS). However, only few studies have considered the association of high-resolution HLA type and MS with none yet from Saudi Arabia. The aim of this study was to investigate the association of HLA class II alleles with MS in the Saudi population. We used next-generation sequencing to investigate HLA association with MS. This study was conducted at King Abdulaziz Medical City in Riyadh, Saudi Arabia. We found that several HLA-DRB1 and DQB1 alleles were associated with MS. These alleles included HLA-DRB1*15:01 (odds ratio [OR]: 3.01; 95%, confidence interval [CI]: 1.68-5.54; P = .0001), HLA-DQB1*02:01 (OR: 1.76; 95% CI: 1.20-2.58; P = .0022), HLA-DQB1*06:02 (OR: 3.52; 95% CI: 1.87-6.86; P < .0001), and HLA-DQB1*06:03 (OR: 2.42; 95% CI: 1.16-5.25; P = 0.01). Interestingly, HLA-DRB1*15:01 was associated with increased risk of previous relapses. In addition, HLA-DRB1*15:01 and HLA-DQB1*06:02 were found to be associated with lower vitamin D levels. This study provides insights on the association of different HLA alleles with clinical characteristics and outcome of MS among Saudis. These insights can have future implications for the clinical management of MS based on the patient genetic profile.

Association of human leukocyte antigen DRB1 with anti-cyclic citrullinated peptide autoantibodies in Saudi patients with rheumatoid arthritis

BACKGROUND

The genetic association between human leukocyte antigen (HLA)-DRB1 alleles and the risk of development of autoantibodies has been investigated, but there are few studies from the Gulf region.

OBJECTIVES

To investigate the association between the HLA-DRB1 shared epitope and the risk for development of autoantibodies in rheumatoid arthritis (RA) patients in a Saudi population.

DESIGN

Analytical cross-sectional study.

SETTING

Tertiary care hospital in Riyadh, Saudi Arabia.

PATIENTS AND METHODS

We enrolled consecutive Saudi RA patients attending the rheumatology clinic between January and April 2015. Previously published data on HLA typing on unmatched healthy controls were used for comparison. HLA typing was performed using sequence-specific oligonucleotide probes (SSOP). Rheumatoid factor (RF), anti-cyclic citrullinated peptide (anti-CCP) antibodies, and antinuclear antibodies (ANA) were also measured. Logistic regression analysis was used to study the autoantibodies as possible explanatory variables for the presence of the HLA-DRB1 shared epitope.

MAIN OUTCOME MEASURE(S)

The association between the presence of the shared epitope and the risk of developing anti-CCP antibodies, ANA, and RF.

RESULTS

In 76 patients with RA, carrying the shared epitope was associated with a significantly higher risk of having RA [OR=2.65, 95% CI (1.42-4.94), P=.0009]. However, only HLA-DRB1*04:05 was significantly as.sociated with RA [OR=3.73, 95% CI (1.61-8.96), Pc=.016]. In the logistic regression analysis, only anti-CCP was significantly associated with the shared epitope [OR=14.51, 95% CI (1.53-137.49), P=.02].

CONCLUSIONS

Our analysis indicates that the presence of the HLA-DRB1 shared epitope is strongly associated with the development of anti-CCP antibodies in Saudi patients with RA.

LIMITATIONS

A larger sample size is needed to confirm our finding.

Association of human leukocyte antigen class II alleles with severe Middle East respiratory syndrome-coronavirus infection

BACKGROUND:

Middle East Respiratory Syndrome (MERS) is a disease of the lower respiratory tract and is characterized by high mortality. It is caused by a beta coronavirus (CoV) referred to as MERS-CoV. Majority of MERS-CoV cases have been reported from Saudi Arabia.

AIM:

We investigated the human leukocyte antigen (HLA) Class II alleles in patients with severe MERS who were admitted in our Intensive Care Unit.

METHODS:

A total of 23 Saudi patients with severe MERS-CoV infection were typed for HLA class II, results were compared with those of 161 healthy controls.

RESULTS:

Two HLA class II alleles were associated with the disease; HLA-DRB1*11:01 and DQB1*02:02, but not with the disease outcome.

CONCLUSIONS:

Our results suggest that the HLA-DRB1*11:01 and DQB1*02:02 may be associated with susceptibility to MERS.

The impact of next-generation sequencing technologies on HLA research

In the past decade, the development of next-generation sequencing (NGS) has paved the way for whole-genome analysis in individuals. Research on the human leukocyte antigen (HLA), an extensively studied molecule involved in immunity, has benefitted from NGS technologies. The HLA region, a 3.6-Mb segment of the human genome at 6p21, has been associated with more than 100 different diseases, primarily autoimmune diseases. Recently, the HLA region has received much attention because severe adverse effects of various drugs are associated with particular HLA alleles. Owing to the complex nature of the HLA genes, classical direct sequencing methods cannot comprehensively elucidate the genomic makeup of HLA genes. Thus far, several high-throughput HLA-typing methods using NGS have been developed. In HLA research, NGS facilitates complete HLA sequencing and is expected to improve our understanding of the mechanisms through which HLA genes are modulated, including transcription, regulation of gene expression and epigenetics. Most importantly, NGS may also permit the analysis of HLA-omics. In this review, we summarize the impact of NGS on HLA research, with a focus on the potential for clinical applications.

Three new HLA-C alleles (HLA-C*14:02:13, HLA-C*15:72 and HLA-C*15:74) in Saudi bone marrow donors

Three new HLA-C alleles were identified by sequence-based typing method (SBT) in donors for the Saudi Bone Marrow Donor Registry (SBMDR). HLA-C*14:02:13 differs from HLA-C*14:02:01 by a silent G to A substitution at nucleotide position 400 in exon 2, where lysine at position 66 remains unchanged. HLA-C*15:72 differs from HLA-C*15:22 by a nonsynonymous C to A substitution at nucleotide position 796 in exon 3, resulting in an amino acid change from phenylalanine to leucine at position 116. HLA-C*15:74 differs from HLA-C*15:08 by a nonsynonymous C to T substitution at nucleotide position 914 in exon 3, resulting in an amino acid change from arginine to tryptophan at position 156.

Screening panel-reactive antibody negative, single-antigen positive: a case report

Sensitized patients remain a challenge for successful transplant. Virtual crossmatch is used to determine the presence or absence of donor-specific antibodies. A 60-year-old woman with a negative screening for panel-reactive antibodies (PRA) received an A*11, A*68 type kidney with a negative anti-human globulin/complement-dependent cytotoxicity (AHG-CDC) crossmatch. Her transplant course was complicated by delayed graft function, and she required hemodialysis. On day 8 after receiving the transplant, she had a kidney biopsy that showed features of antibody-mediated rejection/severe acute tubular necrosis, which was treated by plasmapheresis for 5 sessions and intravenous immunoglobulin (2 g/kg). Her serum level of creatinine decreased from 6.7 to 3.6 mg/dL (600-320 μmol/L). Panel-reactive antibody by Luminex was repeated and again was negative. Single-antigen detection was tried next. Surprisingly, A*11:02 came up positive with a mean fluorescence intensity of 9500. High-resolution donor HLA type was A*68:01 and A*11:01. A*11:02 is not part of the screening Luminex PRA whereas the 11:01 allele is. Serologically, HLA-A11 has 2 defined splits, A11.1 and A11.2, which encode A*11:01 and A*11:02, respectively. In this case, the A*11:02 antibody does not seem to be responsible for the increasing creatinine level. However, if the donor had been A*11:02, a humeral rejection would have occurred and been missed by a virtual crossmatch. Thus virtual crossmatch may not work at all times. Screening for PRA by single antigens is suggested even in PRA-negative cases, if only virtual crossmatch is to be used.

The relevance of HLA sequencing in population genetics studies

Next generation sequencing (NGS) is currently being adapted by different biotechnological platforms to the standard typing method for HLA polymorphism, the huge diversity of which makes this initiative particularly challenging. Boosting the molecular characterization of the HLA genes through efficient, rapid, and low-cost technologies is expected to amplify the success of tissue transplantation by enabling us to find donor-recipient matching for rare phenotypes. But the application of NGS technologies to the molecular mapping of the MHC region also anticipates essential changes in population genetic studies. Huge amounts of HLA sequence data will be available in the next years for different populations, with the potential to change our understanding of HLA variation in humans. In this review, we first explain how HLA sequencing allows a better assessment of the HLA diversity in human populations, taking also into account the methodological difficulties it introduces at the statistical level; secondly, we show how analyzing HLA sequence variation may improve our comprehension of population genetic relationships by facilitating the identification of demographic events that marked human evolution; finally, we discuss the interest of both HLA and genome-wide sequencing and genotyping in detecting functionally significant SNPs in the MHC region, the latter having also contributed to the makeup of the HLA molecular diversity observed today.

Two novel alleles HLA-A*02:433 and HLA-A*02:434 identified in Saudi bone marrow donors using sequence-based typing

In this report, we present two novel HLA-A alleles: HLA-A*02:433 and HLA-A*02:434. These alleles were identified by sequence-based typing method (SBT), in two donors for the Saudi Bone Marrow Donor Registry (SBMDR). Allele A*02:433 is identical to A*02:05:01G except for a G to A substitution at nucleotide position 449 in exon 2. This substitution results in glycine to serine substitution at position 83. Whereas, allele A*02:434 is identical to A*02:01:01G except for a C to A substitution at nucleotide position 245 in exon 2, which results in phenylalanine to threonine substitution at position 15. The generation of both alleles appears to be the result of nucleotide point mutation involving 02:01:01 and 02:05:01.