Gastroenterology Research and Practice
Volume 2019 (2019), Article ID 7369014, 7 pages
https://doi.org/10.1155/2019/7369014

Research Article

HLA Haplotype Association with Celiac Disease in Albanian Pediatric Patients from Kosovo

Atifete Ramosaj-Morina,¹ Marija Burek Kamenaric,² Mehmedali Azemi,¹ Lidvana Spahiu,¹ Zorana Grubic,² and Renata Zunec²

¹Pediatric Clinic, University Clinical Centre of Kosovo, Rrethi i Spitalit p.n., 10000 Prishtina, Kosovo
²Tissue Typing Centre, Clinical Department of Transfusion Medicine and Transplantation Biology, University Hospital Centre Zagreb, Kispaticeva 12, HR-10000 Zagreb, Croatia

Correspondence should be addressed to Renata Zunec; rzunec@kbc-zagreb.hr

Received 7 March 2019; Accepted 19 May 2019; Published 10 June 2019

Guest Editor: Francesco Valitutti

Copyright © 2019 Atifete Ramosaj-Morina et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

Table 1: The HLA-A, -B, and -DRB1 allele frequencies observed in the Albanian pediatric celiac disease patients from Kosovo () with statistically significant differences in comparison to healthy individuals in the control group ().

Patients
()
Control group
()

AF

AF

HLA-A

01
32
0.2667
26
0.1048
<0.0001

02
26
0.2167
81
0.3266
0.0312

HLA-B

08
35
0.2917
14
0.0565
<0.0001

35
4
0.0333
31
0.1250
0.0093

50
7
0.0583
1
0.0040
0.0012

HLA-DRB1

03
46
0.3833
18
0.0726
<0.0001

07
21
0.1750
17
0.0685
0.0023

11
16
0.1333
55
0.2218
0.0460

13
6
0.0500
36
0.1452
0.0102

16
7
0.0583
34
0.1371
0.0289

Legend (alphabetic order): AF = allele frequency; = number of tested individuals; = number of allelic group occurrence.

Table 2: The HLA-DQA1 and -DQB1 allele frequencies observed in the Albanian pediatric celiac disease patients from Kosovo () with statistically significant differences in comparison to healthy individuals in the control group ().

Patients
()
Control group
()

AF

AF

HLA-DQA1

01:02
16
0.1333
66
0.2661
0.0049

02:01
21
0.1750
17
0.0685
0.0023

05:01
48
0.4000
19
0.0766
<0.0001

05:05
16
0.1333
62
0.2500
0.0115

HLA-DQB1

02:01
48
0.4000
18
0.0726
<0.0001

02:02
22
0.1833
19
0.0766
0.0020

03:01
17
0.1417
66
0.2661
0.0084

05:02
8
0.0667
36
0.1452
0.0341

Legend (alphabetic order): AF = allele frequency; = number of tested individuals; = number of allele occurrence.

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.

Table 3: The HLA-A~B~DRB1~DQA1~DQB1 haplotypes with ≥2 number of haplotype copies in the Albanian pediatric celiac disease patients from Kosovo () and compared with the frequencies of those haplotypes in the control group ().

Haplotype HLA
Patients ()
Control group ()

A~B~DRB1~DQA1~DQB1

HF

HF

01~08~03~05:01~02:01
27.0
0.2250
7.0
0.0282
<0.0001

02~50~07~02:01~02:02
5.0
0.0416
1.0
0.0041
0.0311

02~51~11~05:05~03:01
5.0
0.0416
6.0
0.0242
0.3619

68~44~07~02:01~02:02
4.0
0.0334
0
0
0.0049

24~07~15~01:02~06:02
3.0
0.0250
3.0
0.0121
0.3699

23~18~11~05:05~03:01
3.0
0.0250
0
0
0.0755

68~40~03~05:01~02:01
2.0
0.0167
1.0
0.0041
0.2443

02~13~07~02:01~02:02
2.0
0.0167
1.0
0.0041
0.2443

02~08~03~05:01~02:01
2.0
0.0167
1.0
0.0041
0.2443

03~44~07~02:01~02:02
2.0
0.0167
1.0
0.0041
0.2443

24~08~03~05:01~02:01
2.0
0.0167
0
0
0.1303

32~08~03~05:01~02:01
2.0
0.0167
1.0
0.0041
0.2443

03~38~13~01:03~06:03
2.0
0.0167
3.0
0.0121
0.7237

29~27~11~05:01~02:01
2.0
0.0167
1.0
0.0041
0.2443

23~44~07~02:01~02:02
2.0
0.0167
6.0
0.0242
0.6444

03~07~03~05:01~02:01
2.0
0.0167
0
0
0.1303

Legend (alphabetic order): HF = haplotype frequency; = number of tested individuals; = number of observed haplotypes.

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.

Table 4: The frequency of celiac disease-predisposing HLA-DQ haplotypes and genotypes detected among Albanian pediatric celiac disease patients from Kosovo () and the control group ().

Patients ()
Control group ()

HLA-DQ haplotypes
(%)
(%)

Susceptible haplotypes

DQ2.5cis
48 (40.01)
19 (7.66)
<0.0001

DQ2.5trans
8 (6.66)
3 (1.20)
0.0088

DQ8
5 (4.17)
20 (8.06)
0.1712

Fractional susceptible haplotypes

DQ2.2
14 (11.66)
14 (6.45)
0.0631

DQ2.3
0
1 (0.41)
0.8170

DQ7
17 (14.16)
67 (27.01)
0.0068

Nonsusceptible haplotypes

DQ4
1 (0.83)
6 (2.41)
0.3191

DQ5
15 (12.50)
68 (27.41)
0.0017

DQ6
12 (10.00)
48 (19.35)
0.0252

DQ9
0
2 (0.80)
0.5650

HLA-DQ genotypes

DQ2 and/or DQ8 positive
57 (95.00)
50 (40.32)
<0.0001

    Dual-dosage susceptible genotypes

        DQ2.5/DQ2.2
12 (20.00)
2 (1.61)
0.0005

        DQ2.5/DQ2.5
3 (5.00)
0
0.0738

        DQ2.2/DQ2.2
0
2 (1.61)
0.5616

        DQ8/DQ2.5
3 (5.00)
1 (0.81)
0.1091

        DQ8/DQ2.2
0
0
n/a

        DQ8/DQ8
0
1 (0.81)
0.8144

    Sole-dosage susceptible genotypes

        DQ2.5/DQX
35 (58.33)
19 (15.32)
<0.0001

        DQ2.2/DQX
2 (3.33)
8 (6.45)
0.3902

        DQ8/DQX
2 (3.33)
17 (13.71)
0.0459

DQ2 and/or DQ8 negative
3 (5.00)
74 (59.68)
<0.0001

Legend (alphabetic order): DQ2.5 = HLA-DQA105~DQB102 (HLA-DRB103); DQ8 = HLA-DQA103~DQB103:02 (HLA-DRB104); DQ2.2 = HLA-DQA102~DQB102 (HLA-DRB107); DQ2.3 = HLA-DQA103~DQB102 (HLA-DRB104/09/11); DQ7 = HLA-DQB103:01 (HLA-DRB111/12/X); DQ4, DQ5, DQ6, and DQ9 were assigned if HLA-DQB104, HLA-DQB105, HLA-DQB106, and HLA-DQB103:03 alleles were present; DQX = presence of any other HLA-DQB1 allele than HLA-DQ2 and/or HLA-DQ8; = number of tested individuals; = number of HLA-DQ haplotype/genotype occurrence; n/a = not applicable (due to ).

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 [19–22].

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%) [23–25]. 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).

References

S. Husby, S. Koletzko, I. R. Korponay-Szabó et al., “European Society for Pediatric Gastroenterology, Hepatology, and Nutrition Guidelines for the Diagnosis of Coeliac Disease,” Journal of Pediatric Gastroenterology and Nutrition, vol. 54, no. 1, pp. 136–160, 2012.

K. Mustalahti, C. Catassi, A. Reunanen et al., “The prevalence of celiac disease in Europe: results of a centralized, international mass screening project,” Annals of Medicine, vol. 42, no. 8, pp. 587–595, 2010.

M. Leja, Z. Shums, L. Nikitina-Zake et al., “Prevalence estimation of celiac disease in the general adult population of Latvia using serology and HLA genotyping,” United European Gastroenterology Journal, vol. 3, no. 2, pp. 190–199, 2015.

C. Dubé, A. Rostom, R. Sy et al., “The prevalence of celiac disease in average-risk and at-risk Western European populations: a systematic review,” Gastroenterology, vol. 128, no. 4, pp. S57–S67, 2005.

A. Fasano, I. Berti, T. Gerarduzzi et al., “Prevalence of celiac disease in at-risk and not-at-risk groups in the United States: a large multicenter study,” Archives of Internal Medicine, vol. 163, no. 3, pp. 286–292, 2003.

M. F. Kagnoff, “Celiac disease: pathogenesis of a model immunogenetic disease,” Journal of Clinical Investigation, vol. 117, no. 1, pp. 41–49, 2007.

P. L. Stokes, G. K. T. Holmes, P. Asquith, P. Mackintosh, and W. T. Cooke, “Histocompatibility antigens associated with adult coeliac disease,” The Lancet, vol. 300, no. 7769, pp. 162–164, 1972.

J. Ek, D. Albrechtsen, B. G. Solheim, and E. Thorsby, “Strong association between the HLA-Dw3-related B cell alloantigen -DRw3 and coeliac disease,” Scandinavian Journal of Gastroenterology, vol. 13, no. 2, pp. 229–233, 2010.

S. Trabace, A. Giunta, M. Rosso et al., “HLA-ABC and DR antigens in celiac disease. A study in a pediatric Italian population,” Vox Sanguinis, vol. 46, no. 2, pp. 102–106, 1984.

A. S. Louka, S. Nilsson, M. Olsson et al., “HLA in coeliac disease families: a novel test of risk modification by the “other” haplotype when at least one DQA105-DQB102 haplotype is carried,” Tissue Antigens, vol. 60, no. 2, pp. 147–154, 2002.

N. Bratanic, D. Smigoc Schweiger, A. Mendez, N. Bratina, T. Battelino, and B. Vidan-Jeras, “An influence of HLA-A, B, DR, DQ, and MICA on the occurrence of celiac disease in patients with type 1 diabetes,” Tissue Antigens, vol. 76, no. 3, pp. 208–215, 2010.

G. Trynka, K. A. Hunt, N. A. Bockett et al., “Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease,” Nature Genetics, vol. 43, no. 12, pp. 1193–1201, 2011.

R. P. Anderson, M. J. Henry, R. Taylor et al., “A novel serogenetic approach determines the community prevalence of celiac disease and informs improved diagnostic pathways,” BMC Medicine, vol. 11, no. 1, 2013.

“Revised criteria for diagnosis of coeliac disease. Report of Working Group of European Society of Paediatric Gastroenterology and Nutrition,” Archives of Disease in Childhood, vol. 65, no. 8, pp. 909–911, 1990.

K. Dalva and M. Beksac, “HLA typing with sequence-specific oligonucleotide primed PCR (PCR-SSO) and use of the Luminex™ technology,” Methods in Molecular Biology, vol. 1109, pp. 87–99, 2014.

A. K. Lancaster, R. M. Single, O. D. Solberg, M. P. Nelson, and G. Thomson, “PyPop update—a software pipeline for large-scale multilocus population genomics,” Tissue Antigens, vol. 69, pp. 192–197, 2007.

M. P. Horan, S. Y. Chai, N. Munusamy et al., “High rates of variation in HLA-DQ2/DQ8 testing for coeliac disease: results from an RCPAQAP pilot program,” Journal of Clinical Pathology, vol. 71, no. 10, pp. 900–905, 2018.

P. Margaritte-Jeannin, M. C. Babron, M. Bourgey et al., “HLA-DQ relative risks for coeliac disease in European populations: a study of the European Genetics Cluster on Coeliac Disease,” Tissue Antigens, vol. 63, no. 6, pp. 562–567, 2004.

K. Karell, P. Holopainen, K. Mustalahti, P. Collin, M. Mäki, and J. Partanen, “Not all HLA DR3 DQ2 haplotypes confer equal susceptibility to coeliac disease: transmission analysis in families,” Scandinavian Journal of Gastroenterology, vol. 37, no. 1, pp. 56–61, 2009.

E. Bolognesi, K. Karell, S. Percopo et al., “Additional factor in some HLA DR3/DQ2 haplotypes confers a fourfold increased genetic risk of celiac disease,” Tissue Antigens, vol. 61, no. 4, pp. 308–316, 2003.

L. M. Medrano, B. Dema, A. López-Larios et al., “HLA and celiac disease susceptibility: new genetic factors bring open questions about the HLA influence and gene-dosage effects,” PLoS One, vol. 7, no. 10, article e48403, 2012.

J. Gutierrez-Achury, A. Zhernakova, S. L. Pulit et al., “Fine mapping in the MHC region accounts for 18% additional genetic risk for celiac disease,” Nature Genetics, vol. 47, no. 6, pp. 577-578, 2015.

T. Imanishi, T. Akaza, A. Kimura, K. Tokunaga, and T. Gojobori, “Allele and haplotype frequencies for HLA and complement loci in various ethnic groups,” in Proceedings of the Eleventh International Histocompatibility Workshop and Conference, K. Tsuji, M. Aizawa, and T. Sasazuki, Eds., vol. 1, pp. 1065–1220, Oxford University Press, Yokohama, Japan, 1992.

A. Arnaiz-Villena, M. Karin, N. Bendikuze et al., “HLA alleles and haplotypes in the Turkish population: relatedness to Kurds, Armenians and other Mediterraneans,” Tissue Antigens, vol. 57, no. 4, pp. 308–317, 2001.

A. Canossi, D. Piancatelli, A. Aureli et al., “Correlation between genetic HLA class I and II polymorphisms and anthropological aspects in the Chaouya population from Morocco (Arabic speaking),” Tissue Antigens, vol. 76, no. 3, pp. 177–193, 2010.

F. F. González-Galarza, L. Y. C. Takeshita, E. J. M. Santos et al., “Allele frequency net 2015 update: new features for HLA epitopes, KIR and disease and HLA adverse drug reaction associations,” Nucleic Acids Research, vol. 43, no. D1, pp. D784–D788, 2015.

K. Karell, A. S. Louka, S. J. Moodie et al., “HLA types in celiac disease patients not carrying the DQA105-DQB102 (DQ2) heterodimer: results from the European Genetics Cluster on Celiac Disease,” Human Immunology, vol. 64, no. 4, pp. 469–477, 2003.

M. Krini, G. Chouliaras, M. Kanariou et al., “HLA class II high-resolution genotyping in Greek children with celiac disease and impact on disease susceptibility,” Pediatric Research, vol. 72, no. 6, pp. 625–630, 2012.

R. Žunec, Z. Grubić, Z. Jurčić, M. Peršić, A. Kaštelan, and V. Kerhin-Brkljačić, “HLA-DQ2 heterodimer in the diagnosis of celiac disease,” Biochemia Medica, vol. 14, pp. 119–124, 2004.

L. Koskinen, J. Romanos, K. Kaukinen et al., “Cost-effective HLA typing with tagging SNPs predicts celiac disease risk haplotypes in the Finnish, Hungarian, and Italian populations,” Immunogenetics, vol. 61, no. 4, pp. 247–256, 2009.

A. Basturk, R. Artan, and A. Yilmaz, “The incidence of HLA-DQ2/DQ8 in Turkish children with celiac disease and a comparison of the geographical distribution of HLA-DQ,” Gastroenterology Review, vol. 4, no. 4, pp. 256–261, 2017.

N. Selleski, L. M. Almeida, F. C. d. Almeida, C. B. Pratesi, Y. K. d. M. Nóbrega, and L. Gandolfi, “Prevalence of celiac disease predisposing genotypes, including HLA-DQ2.2 variant, in Brazilian children,” Arquivos de Gastroenterologia, vol. 55, no. 1, pp. 82–85, 2018.