Abstract

Genetic predisposition to celiac disease (CD) is strongly associated with the presence of HLA alleles in the individual genotype encoding HLA-DQ2 and/or HLA-DQ8 heterodimers. The main aim of this study was to analyze the HLA-A, -B, -DRB1, and -DQ allele and five-locus haplotype frequencies in 60 Albanian pediatric CD patients and 124 non-CD children from Kosovo. The most prevalent haplotype in patients was the ancestral AH 8.1 haplotype present in 22.5% of the cases compared to 2.8% of the controls (). Additionally, two other haplotypes were also overrepresented in patients (HLA-A02~B50~DRB107~DQA102:01~DQB102:02 and HLA-A68~B44~DRB107~DQA102:01~DQB102:02). Analysis showed that 95.0% of CD patients and 43.3% of controls were carriers of HLA-DQ2 and/or HLA-DQ8 heterodimers. The most frequent CD-predisposing HLA-DQ haplotypes in patients were HLA-DQ2.5 (46.7%) and HLA-DQ2.2 (11.6%), while the most prevalent genotypes were HLA-DQ2.5/DQX (58.3%) and HLA-DQ2.5/DQ2.2 (20.0%). The frequency of the HLA-DQ8 heterodimer among CD patients (4.2%) compared to the control group (8.1%) was without statistical significance. The given data demonstrate differences in the distribution of HLA haplotypes among Albanian CD patients from Kosovo in comparison to other European and non-European populations, as well as provide additional population data to supplement the thus far undisputed importance of the role of HLA-DQ2 and HLA-DQ8 heterodimers in the development of CD.

1. Introduction

Celiac disease (CD) is an autoimmune systemic disorder characterized by chronic inflammation of the small intestinal mucosa triggered by gluten and related prolamine uptake that occurs in genetically susceptible individuals [1]. Europe is historically considered a geographical area with a high incidence of CD, with a prevalence of 1%, which may be higher in Northern European countries [2, 3]. Some studies from Scandinavian countries and the United Kingdom population tended to show a higher prevalence of CD of approximately 1.0-1.5% [4]. Although it was thought that some countries, including the United States, were exempt from this disease, it has recently been shown that it has a similar prevalence as in Europe, 0.5-1.0% of the general population [5].

Susceptibility to CD and its activation and perpetuation involve a combination of environmental and genetic factors through some immunological mechanisms [6]. The involvement of Human Leukocyte Antigens (HLA) in CD pathogenesis was first described in the 1970s [7]. Over the subsequent years, the specific alleles that underlie the described associations became clear: CD is strongly associated with HLA-DQ, which encodes HLA-DQ2 and HLA-DQ8 heterodimers.

At least 90% of the patients with CD are positive for the HLA-DQ2 heterodimer in cis formation (HLA-DQ2cis) encoded by HLA-DQA105:01 (α-chain) and DQB102:01 (β-chain) alleles on a HLA-DRB103 haplotype, although the HLA-DQ2 heterodimer may also be encoded in trans position (HLA-DQ2trans) with the HLA-DQA105:05 allele, usually on HLA-DRB111, DRB112, DRB113:03, and DRB113:05 haplotypes, and the HLA-DQB102:02 allele, usually on a HLA-DRB107 haplotype [8, 9], but in some populations on HLA-DRB104 haplotypes as well. The majority of the remaining CD patients carry the HLA-DQ8 heterodimer formed by one α-chain and one β-chain encoded with HLA-DQA103:01 and HLA-DQB103:02 alleles on the HLA-DRB104 haplotypes [10] and sporadically on HLA-DRB108 haplotypes.

Most of the studies exploring the linkage of HLA genes with CD are focused on HLA class II genes encoding HLA-DQ heterodimers, while very few studies investigate extended HLA haplotypes [11]. In contrast to a relatively conserved HLA-DR3 haplotype, the linkage between HLA loci forming the HLA-DR7 and HLA-DR11 haplotypes are not so strong, and since they can form the HLA-DQ2.5trans heterodimer, we wanted to explore if there is a difference in HLA haplotype distribution among CD patients in comparison to healthy subjects.

The importance of genetic factors in the pathogenesis of CD is well documented by many studies. It is clear that CD rarely develops in the absence of HLA-DQ2 and/or HLA-DQ8 heterodimers, and that the predisposing HLA-DQ2 and HLA-DQ8 subtypes are necessary, but not sufficient for causing the disease [12]. Anderson et al. suggested that a combination of HLA typing and confirmatory serology could reduce the number of unnecessary endoscopies as well as the number of false negatives and/or positive diagnoses [13]. In 2012, the role of HLA in the diagnosis of CD was firmly recognized, which resulted in important changes in diagnostic criteria and the inclusion of HLA typing in CD diagnostic guidelines [1].

As far as we know, no previous study has focused on the frequency of CD-predisposing HLA genotypes in affected and nonaffected individuals in Kosovar Albanian children. Thus, the main aim of this research was to analyze the HLA-A~B~DRB1~DQA1~DQB1 haplotype distribution as well as the frequency of CD-predisposing HLA-DQ genotypes in Albanian pediatric CD and non-CD subjects from Kosovo.

2. Materials and Methods

2.1. Patients and Control Group

Sixty pediatric patients with CD (40 females and 20 males) at the University Clinical Centre of Kosovo (UCCK) and 124 unrelated age- and gender-matched children (64 females and 60 males) without a history of autoimmune diseases were included in this study. The patients’ age at time of diagnosis ranged from 17 months to 18 years, with a mean age of 5.5 years (). The control group age at the time of sample collection ranged from 1 to 18 years with the mean age of 8.7 years (). The CD diagnosis was achieved according to the criteria of the European Society of Pediatric Gastroenterology and Nutrition [14] and the revised guidelines of the European Society of Pediatric Gastroenterology, Hepatology, and Nutrition [1]. The study was approved by the Ethics Committee of the UCCK (approved February 12th, 2013) and by the Ethics Committee of the University Hospital Centre (UHC) Zagreb (approved March 28th, 2013) and is in accordance with the Helsinki Declaration. Written informed consent was obtained from parents or caregivers of all children prior to their enrollment in the study.

2.2. HLA Typing

Three milliliters of EDTA blood samples were collected from each child included in the study. Blood samples were stored frozen at -30°C until all samples were collected. DNA isolation and HLA typing were performed in the Tissue Typing Centre, Department of Transfusion Medicine and Transplantation Biology, UHC Zagreb. Genomic DNA was isolated from whole blood with the MagNA Pure Compact Instrument using the corresponding commercially available MagNA Pure Compact Nucleic Acid Isolation Kit I (Roche Diagnostics GmbH, Mannheim, Germany). The final DNA concentration was adjusted to 50 ng/μl. All CD patients and controls were typed for HLA-A, -B, -DRB1, -DQA1, and -DQB1 applying the standard polymerase chain reaction–sequence-specific oligonucleotide probing (PCR-SSOP) method using the commercially available Immucor LIFECODES HLA-SSO typing kit (Immucor Transplant Diagnostics Inc., Stamford, USA), based on the hybridization of a labeled single-stranded PCR product to SSO probes [15]. After hybridization, the samples were processed and quantified on the Luminex 200 flow analyzer (Luminex Corporation, Austin, USA) and the resulting data were analyzed with Match It! DNA Software, version 1.2.4. The HLA-DQA1 and -DQB1 alleles determined at medium resolution level were assigned to the most common corresponding four-digit HLA allele which was additionally confirmed when necessary by the standard polymerase chain reaction–sequence-specific primer (PCR-SSP) high-resolution protocol (CareDx, Olerup SSP AB, Sweden).

2.3. Statistical Analysis

The observed HLA-A, -B, -DRB1, -DQA1, and -DQB1 frequencies in both research groups and five-locus haplotype estimates were calculated using the PyPop (Python for Population genetics, version 0.7.0; http://www.pypop.org) program [16].

Regarding the HLA-DQ profile of the subjects, they were categorized according to the presence of susceptible haplotypes (DQ2.5cis, DQ2.5trans, and DQ8), fractional susceptible haplotypes (DQ2.2, DQ2.3, and DQ7), or nonsusceptible haplotypes (DQ4, DQ5, DQ6, and DQ9). Additionally, based on the copy number of the CD-predisposing HLA-DQ2 and/or HLA-DQ8 alleles in each individual, dual-dosage susceptible genotypes (DQ2/DQ2, DQ2/DQ8, and DQ8/DQ8), sole-dosage susceptible genotypes (DQ2/DQX or DQ8/DQX), and nonsusceptible (DQ2 and/or DQ8 negative) genotypes were defined. Fisher’s exact test was applied for comparisons between the HLA allelic groups, haplotypes, and genotypes of CD patients and those of the control group, and calculated differences were considered statistically significant if the value was less than 0.05.

3. Results

3.1. HLA Allele Frequencies

The observed HLA allele frequencies of HLA-A, -B, -DRB1, -DQA1, and -DQB1 loci with a statistically significant difference between the Albanian pediatric CD patients and the control group from Kosovo are shown in Table 1 and Table 2.

A total of 15 HLA-A alleles were found in the CD patient group, among which the most frequent were HLA-A01 (26.7%), -A02 (21.7%), and -A03 (10.0%). A statistically significant difference was observed for HLA-A01 (), present with a higher frequency among CD patients compared to the control group, while the HLA-A02 allele had a lower frequency among CD patients ().

At the most polymorphic HLA-B loci with 21 different alleles detected in the CD patient group, seven alleles showed a : HLA-B07 (8.3%), -B08 (29.2%), -B18 (6.7%), -B38 (5.0%), -B44 (8.3%), -B50 (5.8%), and -B51 (10.8%). A significantly higher frequency of HLA-B08 () and HLA-B50 () as well as a significantly lower frequency of HLA-B35 () was detected among CD patients in comparison to the control group.

Out of the 12 different alleles observed at the HLA-DRB1 locus among CD patients, the three most frequent were HLA-DRB103 (38.3%), -DRB107 (17.5%), and -DRB111 (13.3%). Additionally, HLA-DRB103 and -DRB107 are alleles with a much higher frequency among CD patients than in the control group ( and , respectively), as opposed to HLA-DRB111, -DRB113, and -DRB116 with a significantly lower frequency among CD patients (, , and , respectively).

At the HLA-DQ locus, 10 different HLA-DQA1 and 12 different HLA-DQB1 alleles were observed in CD subjects. The most frequent HLA-DQA1 alleles were HLA-DQA101:02 (13.3%), -DQA102:01 (17.5%), -DQA105:01 (40.0%), and -DQA105:05 (13.3%), with a statistically significant difference and higher frequency for HLA-DQA102:01 () and -DQA105:01 () in comparison to the control group. On the other hand, the frequency of the HLA-DQA101:02 () and -DQA105:05 () alleles is lower among CD patients. HLA-DQB102:01 (40.0%), -DQB102:02 (18.3%), and -DQB103:01 (14.2%) were the most frequent alleles of the HLA-DQB1 locus in the CD patient group, and a significantly higher frequency in the CD patient group in comparison to controls was observed just for HLA-DQB102:01 () and -DQB102:02 (). On the other hand, the frequencies of HLA-DQB103:01 and -DQB105:02 alleles were significantly lower among CD patients ( and , respectively) in comparison to the control group.

3.2. HLA-A~B~DRB1~DQA1~DQB1 Haplotype Frequencies

A total of 71 different five-locus haplotypes (HLA-A~B~DRB1~DQA1~DQB1) were detected in the CD patient group with a total of 55 haplotypes that occurred only once. The most prevalent haplotype observed was HLA-A01~B08~DRB103~DQA105:01~DQB102:01 (22.5%) followed by HLA-A02~B50~DRB107~DQA102:01~DQB102:02 and HLA-A02~B51~DRB111~DQA105:05~DQB103:01 with a frequency of 4.1%, each. When comparing the CD patients’ haplotype frequency results to the frequencies of 156 different five-locus haplotypes observed in the control group (Table 3), a statistically significant difference was calculated for the HLA-A01~B08~DRB103~DQA105:01~DQB102:01 () and HLA-A02~B50~DRB107~DQA102:01~DQB102:02 () haplotypes, both more frequent among CD patients. Also, the HLA-A68~B44~DRB107~DQA102:01~DQB102:02 haplotype was observed with a high frequency (3.3%) among the CD patient group, which resulted in a significant difference () due to a 0% frequency of this haplotype in the control group.

3.3. Frequencies of CD-Predisposing HLA-DQ Haplotypes and Genotypes

The presence of the HLA-DQ2 and/or HLA-DQ8 heterodimers was detected in 57/60 (95.0%) CD patients and in 50/124 (40.3%) control subjects. The detailed CD-predisposing HLA-DQ haplotype and genotype distribution among CD patients and a comparison with the frequency of those haplotypes/genotypes in the control group is presented in Table 4. The most frequent CD-predisposing HLA-DQ haplotypes found in patients were HLA-DQ2.5 (cis and trans conformation) and HLA-DQ2.2 with a frequency of 46.7% and 11.6%, respectively. Consequently, the most prevalent genotype among CD patients was HLA-DQ2.5/DQX with a frequency of 58.3%, followed by the HLA-DQ2.5/DQ2.2 genotype with a frequency of 20.0%. When compared with the control group, both genotypes show a statistically significantly higher frequency among CD patients ( and , respectively). Three patients were homozygous for HLA-DQ2.5 (having a genotype positive for two copies of the HLA-DQ2.5 heterodimer), with no presence of this genotype in controls. The frequency of the HLA-DQ8 haplotype among CD patients is low (4.2%), and there is no significant difference compared to the control group (8.1%). Looking at the HLA-DQ8-related genotypes, 3 (5.0%) CD patients and 1 (0.8%) individual from the control group carried HLA-DQ8/DQ2.5, while 2 (3.3%) CD patients, but as many as 17 (13.7%) individuals from the control group, carried the HLA-DQ8/DQX genotype (). Only one individual in the control group had the HLA-DQ8 genotype alone (DQ8/DQ8). Three CD patients who are not carriers of any of the susceptible variants (HLA-DQ2.5, HLA-DQ2.2, or HLA-DQ8) had the HLA-DQ7/DQ7, HLA-DQ7/DQ5, and HLA-DQ5/DQ6 genotypes.

4. Discussion

A strong association between the CD and HLA-DQ allelic groups encoding for HLA-DQ2 and HLA-DQ8 heterodimers is well-described, and these specific HLA heterodimers are usually seen in more than 90.0% of the patients with CD [17]. However, due to the huge polymorphism of the HLA system and the different HLA risk levels of CD among populations [18], it is important to perform the analysis in each specific population separately. This is the first analysis of the HLA-A, -B, -DRB1, -DQA1, and -DQB1 allele frequencies and the HLA-DQ2 and HLA-DQ8 heterodimer frequencies in 60 Albanian children from Kosovo diagnosed with CD. The limitation of the study is the small sample size, although the post hoc sample size calculation showed that the number of 60 CD patients was enough for the obtained results with a confidence interval of 85%.

The analysis of the HLA-A, -B, -DRB1, and -DQ frequencies among CD patients and compared to the control group revealed a pronounced increase of the HLA-A01, -B08, -B50, -DRB103, -DRB107, -DQA102:01, -DQA105:01, -DQB102:01, and DQB102:02 alleles. On the other hand, HLA-A02, -B35, -DRB111, -DRB113, -DRB116, -DQA101, DQB105, and -DQB103:01 were significantly less present among CD patients than in controls. The higher frequency of these HLA alleles among CD patients is also reflected in the haplotype frequency. The top-ranked haplotype observed with a high frequency of 22.5% was the ancestral AH 8.1 haplotype (HLA-A01~B08~DRB103~DQA105:01~DQB102:01), which was by contrast present in 2.8% of controls. These results are in concordance with results showing that DQ2 is an absolute requirement for the development of CD, but the presence of the well-known ancestral haplotype AH 8.1 and additional genetic factors in the HLA class I region induce an increased risk of CD [1922].

What is also interesting are the two HLA alleles with a surprisingly high frequency in the CD patient group: HLA-B50 and HLA-DRB107, with frequencies that were more than four times higher and three times higher in comparison to controls, respectively. Consequently, the second most frequent five-locus haplotype observed in our CD patients was HLA-A02~B50~DRB107~DQA102:01~DQB102:02. This haplotype has not been previously reported as a haplotype associated with CD, and it is mostly associated with a probable Euro-Asiatic origin, having been reported in Mongolians (HF: 3.2%), in the Chaouya population from Morocco (HF: 2.9%), Turks and Kurds (HF: 1.3%), Spaniards (HF: 1.2%), and Italians (HF: 0.5%) [2325]. The second uncommon and unexpectedly frequent haplotype, ranked fourth among our CD patients, was HLA-A68~B44~DRB107~DQA102:01~DQB102:02. This haplotype was not observed in the control group, and this five-locus haplotype has not been reported in any European population so far, but only in the Sri Lanka Colombo population with a frequency of 7.0% [26].

These results raise new questions: what is the distribution of this five-locus HLA haplotype in other neighboring populations as well as in other populations of European origin? Furthermore, is the presence of this haplotype specific just for Albanian CD patients from Kosovo, or is this haplotype present among CD patients from other countries?.

The results of this research revealed that 95.0% of the Albanian pediatric CD patients from Kosovo were HLA-DQ2- and/or HLA-DQ8-positive compared to the 40.3% of positive individuals in the control group. The given results are similar to those reported from other studies, and at the same time they confirm the variability of HLA-DQ2 and HLA-DQ8 heterodimer frequencies among CD patients from different populations. One large European collaborative study comprising 1008 CD patients from Finland, France, Italy, Norway, Sweden, and the UK showed that HLA-DQ2 and HLA-DQ8 heterodimer frequencies were higher in the Northern European population (Finland—96.0%, Norway + Sweden—96.6%, UK—95.7%) than in the Southern European population (France—93.4% and Italy—89.5%) [27]. Since Kosovo is situated in Southeastern Europe, our result of 95.0% DQ2- and DQ8-positive patients is in concordance with this observation. Furthermore, our results are very similar to results from Greece where a single centre study found that 95.8% of pediatric CD patients were HLA-DQ2- and/or HLA-DQ8-positive in comparison to 32.5% of healthy individuals [28]. On the other hand, a study from Croatia reports a much higher frequency (98.0%) of patients with CD who were carriers of HLA-DQ2 and HLA-DQ8 heterodimers [29].

Epidemiological HLA studies have shown that the HLA-DQ gene dose has a strong quantitative effect on the magnitude of gluten-specific T cell responses, and these individuals have the highest risk of developing CD [10]. In the present study, 18 (30.0%) patients had a dual dosage of HLA-DQ-susceptible genotypes; 12 were HLA-DQ2.5/2.2 heterozygous, 3 HLA-DQ2.5 homozygous, and 3 were HLA-DQ2/DQ8 heterozygous. The majority of HLA-DQ2-positive patients were found to be homozygous or heterozygous for the HLA-DRB103~DQA105:01~DQB102:01 haplotype (HLA-DR3~DQ2), so the likelihood of HLA-DQ2 positivity (91.6%) is in line with reports from the European populations (90.0-95.0%) [18, 23]. Interestingly, the incidence of the HLA-DQ8 heterodimer (4.2%) alone, in double dose or in combination with other heterodimers, was lower in our CD cases than in controls, but with no statistically significant difference. Those result are in line with the studies from different European populations reporting the incidence of the HLA-DQ8 heterodimer among CD patients from 2.8% to 7.9% [27, 29, 30], although there are also studies reporting much higher HLA-DQ8 heterodimer frequencies, even up to 25.0% [28, 31].

In our cohort, 5.0% (3/60) of CD patients were HLA-DQ2- and/or HLA-DQ8-negative, which is comparable to the data originating from different countries reporting the percentage of HLA-DQ2/DQ8-negative CD patients between 0 and 10.0% [27, 29, 30, 32]. Two out of three HLA-DQ2/DQ8-negative CD patients were positive for the DQA1 part of the DQ2 heterodimer, carrying the HLA-DRB111~DQA105:05~DQB103:01 haplotype. The third patient was negative for all CD-predisposing DQ alleles (HLA haplotypes: HLA-DRB113~DQA101:03~DQB106:03 and HLA-DRB115~DQA101:02~DQB105:02), with positive serology and partial villous atrophy who responded to a gluten-free diet.

In conclusion, this study supports the idea that possibly, in addition to the well-known association of CD with HLA class II alleles, extended HLA haplotypes might be considered as a potential genetic risk factor. The given results also provide population data of Albanian CD patients from Kosovo and support the importance of HLA-DQ2 and HLA-DQ8 heterodimers in the development of CD.

Data Availability

The HLA data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

The abstract of this paper was presented at the 33rd European Immunogenetics and Histocompatibility Conference as a poster presentation with interim findings. The poster’s abstract was published in “Special Issue: Abstracts for the 33rd European Immunogenetics and Histocompatibility Conference, Lisbon, Portugal, May 8–11, 2019” in HLA Immune Response Genetics, 93(5): 249-413. https://onlinelibrary.wiley.com/doi/full/10.1111/tan.13518.

Conflicts of Interest

The authors declare that there is no conflict of interest.

Authors’ Contributions

Atifete Ramosaj-Morina and Marija Burek Kamenaric contributed equally to this manuscript.

Acknowledgments

A special acknowledgment goes to the Tissue Typing Centre, Clinical Department of Transfusion Medicine and Transplantation Biology, University Hospital Centre Zagreb for their help with HLA typing for patients from Kosovo. We also want to thank the nurses of the Department of Gastroenterology, Pediatric Clinic, University Clinical Centre of Kosovo and Main Family Medical Centre, Prishtina, Kosovo, for their help during the blood sample collection for patients and controls. This work was supported in part by the Ministry of Education, Science, and Technology, Republic of Kosovo (grant number 2347-04/12/2014).