Skip to main content

Contribution of HLA class I (A, B, C) and HLA class II (DRB1, DQA1, DQB1) alleles and haplotypes in exploring ethnic origin of central Tunisians

Abstract

Background

Estimation of HLA (Human leukocyte Antigen) alleles’ frequencies in populations is essential to explore their ethnic origin. Anthropologic studies of central Tunisian population were rarely reported. Then, in this work, we aimed to explore the origin of central Tunisian population using HLA alleles and haplotypes frequencies.

Methods

HLA class I (A, B, C) and HLA class II (DRB1, DQA1, DQB1) loci genotyping of 272 healthy unrelated organ donors was performed by Polymerase Chain Reaction-Sequence Specific Oligonucleotide (PCR-SSO). We compared central Tunisians with other populations (Arabs, Berbers, Mediterraneans, Europeans, Africans, etc.) using alleles and haplotypes frequencies, genetic distances, Neighbour-Joining dendrogram and correspondence analysis.

Results

Among the 19 HLA A alleles, the 26 HLA B alleles, the 13 HLA C alleles, the 15 HLA DRB1 alleles, the 6 HLA DQA1 alleles and the 5 HLA DQB1 alleles identified in the studied population, HLA A*02 (22.8%), HLA B*50 (13.1%), HLA C*06 (21.8%), HLA DRB1*07 (17.8%), HLA DQA1*01 (32.1%) and HLA DQB1*03 (31.6%) were the most frequent alleles. The extended haplotypes HLA A*02-B*50-C*06-DRB1*07-DQA1*02-DQB1*02 (1.97%) was the most frequent HLA six-loci haplotype.

Conclusion

Central Tunisians were very close to other Tunisian populations, to Iberians and North Africans. They were rather distant from sub-Saharan populations and eastern Mediterraneans especially Arabs although the strong cultural and religious impact of Arabs in this population.

Peer Review reports

Introduction

The HLA system is a genetic system located on the short arm of chromosome 6, containing up to 35,000 alleles, registered in the IMGT/HLA database [1, 2]. HLA system codes for surface proteins that initiate the immune response. HLA antigens recognize the foreign peptide, present it to T lymphocytes and thus activate humoral and/or cellular response. However, HLA system is highly polymorphic, enabling the presentation of an enormous repertoire of peptides. This great genetic variability makes this system a genetic predisposition factor for certain diseases, and an obstacle to graft survival in some organ transplantations. Besides, estimation of HLA alleles frequencies in general population is essential in anthropological studies [2,3,4]. In fact, the high polymorphism of HLA system makes it a valuable tool for classifying populations, exploring their ethnic origins and their immigration patterns. HLA genetic studies have been used in many anthropological studies including African anthropological studies [4,5,6,7,8,9,10,11,12,13,14,15,16].

The main origin of present-day Tunisians is the native Berbers. But the other civilizations that settled in Tunisia also had an impact on the heritage of Tunisians. Indeed, Tunisia, including the central Tunisian region, has seen a succession of several civilizations. It was invaded by the Phoenicians which emigrated from the Middle East to the Mediterranean coast between 1000 and 500 BC. In 815 BC they created Carthage. After that, Romans ruined Carthage in 146 BC and colonized Tunisia. Tunisia was then attacked by the Vandals in the fifth century. Muslims came to Tunisia in the seventh century to found Kairouan; the first Islamic city in North Africa. The Muslim invasion of Tunisia continued between the seventh and fifteenth centuries, with the arrival of Arab tribes from Egypt and the Middle East [17]. Tunisia was then invaded by the Ottoman (Turks) in 1574 AD, followed by French colonization until independence in 1956 [18]. Like all regions of Tunisia, the central region is characterized by great ethnic diversity. Our study focuses on the central Tunisian region especially Sousse, Kairouan, Kasserine and Sidi Bouzid governorates (Fig. 1). Central region of Tunisia comprises nearly 19% of all population according to official census in 2014. Inhabitants of central Tunisia are mainly Arab-speaking populations and Berbers.

Fig. 1
figure 1

Map showing the central Tunisian governorates included in this study

The aim of this study was to explore the most likely origin of the population of central Tunisia which was the starting point (especially Kairouan) for the spread of Arabo-muslim civilization in North Africa. This investigation was made using HLA alleles frequencies and comparison with other Tunisian, Mediterranean, sub-Saharan and Arab populations (Table 1).

Table 1 Populations included in the present study

Materials and methods

Study subjects

Our study included retrospectively 272 healthy unrelated donors from solid organ donors’ database of hematology’s laboratory in Sahloul Hospital of Sousse, Tunisia. All individuals included in the study were donors of solid organs and were from the center of Tunisia namely Sousse governorate, Kairouan governorate, Kasserine governorate and Sidi Bouzid governorate. No ethnic or linguistic criterion was used to include subjects in this study. To participate in the study, informed consent was obtained from all participants. This study was approved by the local ethics committee of Sahloul university Hospital. All data and identities of patients were processed with strict confidentiality.

Laboratory analysis

HLA-A, HLA-B, HLA-C, HLA-DR, and HLA-DQ genotyping

EDTA anticoagulated venous blood samples were used for HLA genotyping. Extraction of genomic DNA from these samples was performed using the “QIAamp DNA mini kit” (Qiagen, Hilden, Germany). Gene typing of the HLA-A, HLA-B, HLA-C, HLA-DR, and HLA-DQ loci was performed by the Polymerase chain reaction-sequence specific oligonucleotide (PCR-SSO) technique with the LABType™ SSO typing kits (One lambda Inc, USA). Interpretation was performed on the HLA FUSION™ software.

Statistical analysis

In the present study, HLA alleles of the six HLA loci; HLA-A, -B, -C, -DRB1, -DQA1 and -DQB1 frequencies were determined using the gene counting method. Frequencies of haplotypes were estimated through the maximum-likelihood approach from the related genotypic data using the expectation–maximization algorithm [37, 38] via The Arlequin v3.5.2.2 program [39]. To evaluate Hardy–Weinberg equilibrium (HWE) at each locus, a modified Markov-chain random walk approach with 100,172 steps was applied.

Linkage disequilibrium (LD) between specific alleles at distinct loci, along with the significance level (p < 0.05) for 2 X 2 comparisons and the relative linkage disequilibrium (RLD or D’), were determined as previously outlined [40]. For phylogenetic tree analysis (dendrograms), the Neighbour-Joining (NJ) method [41] was utilized along with standard genetic distances [42]. The DISPAN software containing GNKDST and TREEVIEW programs [43, 44] facilitated these analyses and the tree construction. Furthermore, a three-dimensional correspondence analysis and its two-dimensional representation were performed with R language using the packages Factoshiny and FactoMineR [45]. This analysis provided an overview of the population relationships in the context of HLA allele frequencies.

Results

Alleles frequencies in the studied population

In this study of central Tunisian population, the distribution of HLA genotypes class I and class II were in Hardy Weinberg equilibrium (Table 2). The frequencies of HLA A, B, C, DQA1, DRB1 and DQB1 of Tunisians from the center are illustrated in Table 3. Of the 19 HLA A alleles, A*02 (22.8%), A*01 (13.4%) and A*30 (12.7%) were the most frequent alleles in central Tunisians. They were also the most frequent alleles in Moroccans [35, 46] and in Spanish people [10]. These alleles especially A*01 and A*02 were also in high frequencies in other Tunisian, north African and Iberian populations [6, 7, 11, 12, 16, 21, 47]. As for the 26 HLA B alleles identified, HLA-B*50 (13.1%), HLA-B*44 (9.2%) and HLA-B*51 (9%) were the most frequent alleles. HLA-B*50 and HLA-B*51 are common alleles in Mediterraneans [12, 16, 46]. HLA-B44 and HLA-B50 were also among the most prevalent alleles in the Metalsa Berbers Moroccan population [35]. HLA-B*44 allele was also the most frequent in Tunisian Berbers [9], in the Tunisians-B [21], Spaniards Basques [47] and in high frequencies in Algerians [12], French Basques[48] and in Swiss [49]. Among the 13 HLA C alleles identified, HLA-C*06 (21.8%), HLA-C*07 (18.1%) and HLA-C*04 (13.1%) were the most frequent alleles. In the same way, the three HLA C alleles were the most frequent in North Tunisians [8]. The same result was found in Algerians [12], Moroccans [35] and Spaniards [47].

Table 2 Hardy–Weinberg equilibrium and heterozygosity
Table 3 HLA-A, -B, -DRB1 and -DQB1 allele frequencies

As regards the HLA class II alleles, 15 HLA DRB1 alleles were identified. Among them, HLA-DRB1*07 (17.8%), HLA-DRB1*04 (16.1%) and HLA-DRB1*03 (13.5%) were the most frequent alleles. HLA-DRB1*07 was also the most frequent allele in Berbers [9] and other Tunisians [6, 16]. HLA-DRB1*03 was also observed in high frequencies in many Mediterranean populations [7, 12, 46, 47]. HLA-DRB1*04 was observed in high frequency in Moroccans [46] and in ghanouchians [6]. Of the 6 HLA DQA1 alleles identified, HLA-DQA1*01 (32.1%) and HLA-DQA1*05 (30%) were the most frequent alleles. These two alleles were also the most frequent HLA DQA1 alleles in Spaniards, Basques [47], in Algerians [12] and in Moroccans [13, 46]. HLA-DQA1*05 and HLA-DQA1*01 were also, in addition of HLA-DQA1*02, the most frequent alleles in North Tunisians[8]. In addition, of the 5 HLA DQB1 alleles identified, HLA-DQB1*03 (31.6%) and HLA-DQB1*02 (29%) were the most frequent alleles. These two alleles were also frequent in many other Mediterranean populations [6, 7, 9, 11, 16, 46, 47].

Allelic comparison between Tunisians and other populations

Population of central Tunisia were compared with other Mediterranean, Arab-speaking and worldwide populations, using generic HLA-DRB1 and HLA-DQB1 allele frequencies data.

The allelic comparison was done at the levels of Neighbor-Joining (NJ); (Fig. 2 and 3), correspondence analysis (Fig. 4), and standard genetic distances (SGD); (Table 4).

Fig. 2
figure 2

Neighbor-Joining dendrogram showing relatedness between central Tunisians and other populations using generic genotyping of HLA-DQB1 and HLA-DRB1. Standard genetic distances (SGD) between populations were determined. Populations’ data are from references detailed in Table 1. Bootstrap values from 1.000 replicates are shown

Fig. 3
figure 3

Neighbor-Joining dendrogram showing relatedness between central Tunisians and other populations using generic genotyping of HLA-A and HLA-B. Standard genetic distances (SGD) between populations were determined. Populations’ data are from references detailed in Table 1. Bootstrap values from 1.000 replicates are shown

Fig. 4
figure 4

Correspondence analysis (bi-dimensional representation) based on the standard genetic distances, showing a global view of the relationship between central Tunisians and other Mediterranean populations according to HLA DRB1, and HLA DQB1 allele frequencies data

Table 4 Standard genetic distances (SGD) between central Tunisians and other populations

NJ dendrogram

We compared generic HLA-DRB1 and HLA-DQB1 alleles’ distribution between central Tunisian population and other populations. Comparison was done with NJ dendrogram based on SGD (Fig. 2). Results from NJ dendrogram show four clusters. The first cluster grouped North Africans (Algerians, Moroccans, Tunisians Berbers and Libyans), Europeans (Spanish, Portuguese, French, Belgians, Germans and Austrians) and Saudi Arabians. The second cluster comprised Eastern Mediterraneans (Palestinians, Lebanese, Turks and Greeks), Moroccan Jews and Italians. The third cluster grouped together some Iberian populations (Cantabrians and basques) and some Tunisian populations (Tunisian population A and Tunisia Ghannouch). The fourth cluster contained Sub-Saharan African populations and Tunisia Jerba Berbers. Macedonians formed an out-group.

We compared also generic HLA-A and HLA-B alleles’ distribution between central Tunisian population and other populations, with NJ dendogram based on SGD (Fig. 3). NJ dendogram shows four clusters. The first cluster grouped Europeans (Cantabrians, Portuguese, French, Belgians, Germans, Austrians and basques), Greeks and Macedonians. The second cluster comprised North Africans (Moroccans and Tunisians), sudanese and Saudi Arab. Eastern Mediterraneans (Palestinians, and Jordanians) form the third cluster. The fourth cluster grouped together some Sub-Saharan African populations and Tunisians from Ghannouch.

SGD comparison

Calculation of SGD between central Tunisian population and other populations included in this study was based on generic HLA DRB1 and HLA DQB1 allele frequencies (Table 4). SGD results show that central Tunisians were closer to Western Mediterranean populations than to Eastern Mediterranean populations. In fact, Arab Gabesian Tunisians had the closest genetic distance (5 × 10−3) to central Tunisian population followed by Gabesian Tunisians (1.1 × 10−2), Portuguese (1.1 × 10−2), Belgians (1.3 × 10−2), Tunisians-B (1.5 × 10−2), Moroccans from Settat Chaouya (1.5 × 10−2), Northern Tunisians (1.6 × 10−2), Algerians from Oran (1.7 × 10−2), Moroccans from Atlantic Coast of Chaouya (1.8 × 10−2) and Tunisian Berbers from Matmata (1.9 × 10−2). Thus, Tunisians from the center seem to be distant from Eastern Mediterranean populations including Arab populations (Lebanese (4.7 × 10−2) and Palestinians (2.6 × 10−2)).

Correspondence analysis

Correspondence analysis using generic HLA DRB1 and HLA DQB1 shows four clusters (Fig. 4). The first combines Eastern Mediterraneans (Turks, Greeks, Macedonians, Palestinians and Lebanese), Moroccan Jews and Italians. The second grouped together Europeans (Cantabrians, Portuguese, French, Belgians, Austrians and Germans). The third comprised Western Mediterraneans (Algerians, Moroccans, Tunisians, Libyans and Spanish) and Saudi Arabians. The fourth cluster contained Sub-Saharan African populations, Tunisia Jerba Berber and Tunisians from Matmata Berber. Tunisians from Ghannouch and basques form an out-group.

HLA A, B, DRB1, DQA1 and DQB1 linkage disequilibrium (LD)

The study of central Tunisian HLA haplotypes (Table 5) allows their comparison with other haplotypes found in other populations. Table 5 represents HLA (A, B) and HLA (DRB1, DQB1) two-loci haplotypes with significant LD (p < 0.05 in all cases) in central Tunisians. The most frequent HLA two-loci haplotypes in central Tunisians were DRB1*07–DQB1*02 (14.94%), followed by DRB1*03–DQB1*02 (11.78%) and DRB1*11–DQB1*03 (11.78%). DRB1*07–DQB1*02 haplotype was also the most frequent two-loci haplotype in southern Tunisians (18.02%)[7], in Arab Gabes Tunisians (15%), Tunisians from Sousse (17.7%) [16] and in Tunisian Berbers(16.03%)[9]. This two loci haplotype was also present in many other Mediterranean populations at high frequencies namely Moroccans (12.6%) [46] Ghannouchians(16.46%) and Spaniards (17.3%)[47].

Table 5 HLA class I (A, B) and class II (DRB1, DQB1) two-loci haplotypes with significant linkage disequilibrium (P < 0.05 in all cases) in Central Tunisians

DRB1*03–DQB1*02, known as an Iberian paleo-North African haplotype, was among the most frequent two-loci haplotype in central Tunisian population (11.78%). It was also common in many other west Mediterranean populations. Indeed it was frequent in Tunisian Berbers (11.26%) [9], Tunisians-B (14.08%) [21], Algerians(11.3%)[12], Moroccans (15.3%)[46], Spaniards(13.6%)and Basques (17.5%) [47]. Concerning HLA (A, B) two-loci haplotypes, the most frequent in the population of the Tunisian center were A2-B51 (4.77%) and the A2-B50 haplotype (4.41%). These two-loci haplotypes were also frequent in Berbers, other Tunisian and Mediterranean populations [6, 7, 9, 16, 47].

Concerning HLA (DRB1, DQA1, DQB1) three loci haplotypes, the most frequent were HLA DRB1*07-DQA1*02-DQB1*02 (15.15%) and HLA DRB1*03-DQA1*05-DQB1*02 (12.1%). This was concordant with Tunisian and Mediterranean studies. HLA DRB1*07-DQA1*02-DQB1*02 and HLA DRB1*03-DQA1*05-DQB1*02 were among the most frequent three loci haplotypes in Northern Tunisians [8], Moroccans [46], Algerians [12], Spaniards and Basques [47] (Table 6).

Table 6 HLA (DRB1, DQA1, DQB1) three-loci haplotypes with significant linkage disequilibrium (P < 0.05 in all cases) in Central Tunisians

HLA class I and class II extended haplotype analysis

Extended A-B-DRB1-DQB1 haplotypes and their frequencies in central Tunisian population are represented in Table 7. The most frequent A-B-DRB1-DQB1 extended haplotype in the population of the center of Tunisia was A*02-B*50-DRB1*07-DQB1*02 (3.94%) followed by A*02-B*51-DRB1*13-DQB1-06 (1.97%). HLA A*02-B*50-DRB1*07-DQB1*02 haplotype was also the most frequent haplotype in Tunisians (from the North and the south of Tunisia) (2.2%) [16, 50], in Tunisian Berbers (8.1%) [9], Southern Tunisians(3.2%)[7], in Gabesians(2.6%)[19], in Ghanouchians (1.8%) [6], in Tunisians-B (1.2%) [21], in Spaniards (1.2%) [47], in Turks (1.3%) [51] and Moroccan Jews (2%) [25]. The second most frequent four loci-haplotype A*02-B*51-DRB1*13-DQB1*06 (1.97%), another Ibero-berber haplotype [46], was also frequent in Moroccans [46],and Spaniards [10].

Table 7 The most frequent four-loci haplotype in Central Tunisians

The most frequent extended six-Loci haplotypes (A, B, C, DRB1, DQA1, DQB1) and their frequencies are represented in Table 8. The extended haplotypes HLA A*02 -B*50- C*06-DRB1*07-DQA1*02-DQB1*02 (1.97%) and HLA A*02-B*51-C*16-DRB1*07-DQA1*02-DQB1*02 (1.64%) were the most frequent HLA six-loci haplotypes. The first haplotype was also the most frequent in Moroccans from Chouaya[13] and among the most frequent HLA six-loci haplotypes in Northern Tunisians[8]. However, the second extended haplotype HLA A*02-B*51-C*16-DRB1*07-DQA1*02-DQB1*02 was absent in Iberian, Mediterranean and Tunisian populations [6,7,8, 11,12,13, 32, 46].

Table 8 The most frequent six-loci haplotype in Central Tunisians

Discussion

To the best of our knowledge, this is the first anthropologic published study of the population of the central Tunisian region in terms of HLA alleles and haplotypes, and the relationship of this population to other populations. For this study, we did not only use allele and haplotype frequencies, but also we used methods of anthropological and evolutionary analysis (genetic distance, correspondence analysis, dendrograms, etc.).

Central Tunisians, other Tunisians and North Africans

The study of genetic distances and allelic frequencies has shown that the central Tunisian population is closely related to the southern Tunisians, especially the population of Gabes [5, 19], and to northern Tunisians [8, 21]. This shows that central Tunisians are not distinct from other Tunisians. Indeed, all regions share almost the same historical events. However, the population of Jerba Island was genetically rather distant from the central Tunisian region. Moreover, this genetic distance from the Jerba region was also found in other North Africans. This could be explained by the low level of exogamy within the Jerba Berber population (less than 6%) [5]. Besides, according to Khodjet el khil et al., the Jerba population is rather close to some European populations [52]. Then, central Tunisians are closely related to most of Tunisians. Nevertheless, we have found in our population an extended haplotype HLA A*02-B*51-C*16-DRB1*07-DQA1*02-DQB1*02 which was absent in other Tunisian and Mediterranean populations [6,7,8, 11,12,13, 32, 46]. Thus, this haplotype could be a characteristic of the population of central Tunisia and reflects the important mixing of populations in this region.

Furthermore, central Tunisians were much closer to western Mediterraneans, particularly North Africans, than to the eastern Mediterraneans. This was confirmed by NJ dendrogram, corresponding analysis, allele frequencies and haplotypes frequencies. These results were also found in the other Tunisian populations [5, 7, 8, 19, 21]. The relatedness found in our study between the population of central Tunisia and North African populations especially Algeria and Morocco can be explained by almost similar historical events. Indeed, the original inhabitants of North Africa (including Tunisians) were the Berbers. These countries were then colonized by the Phoenicians around 1000 BC, then the Romans after Punic wars (264–266 BC). Later, there were Arabs and Muslims conquest in seventh century AD, and important Bedouin immigration in the eleventh century. Then, there was immigration of Andalusians and Negroid slaves. This similar succession of civilizations in each of the countries of North Africa explains the links of relatedness between these different countries [50, 53, 54].

Central Tunisians and Iberians

NJ dendrogram, genetic distances, corresponding analysis, alleles and haplotypes comparison have shown that the closest populations to central Tunisians were Western Mediterranean populations especially North Africans and Iberians. This result was also found in most of Tunisian regions [5, 7, 8, 16, 19, 21, 50]. The relatedness between central Tunisians, North Africans and Iberians is explained by many historical events. In fact, North African Berbers (including Tunisians and Central Tunisians) were forced to immigrate to the northern coast of the Mediterranean (Spain, Portugal, etc.) due to hyper-arid climatic conditions. This occurred likely in 10,000–4000 BC [27]. On the other hand, North Africans and Iberians were colonized by nearly the same civilizations. Both were invaded by Phoenicians, Romans, Germans, Arabs and Muslims. Besides, Muslim invasion of the Iberians was launched mainly from North Africa and led by North African Muslims and Berbers. The latter settled in Spain specifically for 8 centuries [7, 55]. These events and observations confirm the gene flow and admixture between Iberians and North Africans (including Central Tunisians) and support their relatedness.

Central Tunisians, Berbers and Arabs

This genetic study of HLA system has revealed that central Tunisians, Algerians, Moroccans, North Tunisians and Southern Tunisians are related to Berbers. This result is not surprising. It is consistent with geography and ancestry. Then, Tunisians and North Africans of the present days are not genetically different from Berbers [9, 16]. However, the genetic impact of Eastern Mediterranean populations, particularly Arabs, on the population of central Tunisia was not significant. Indeed, this study showed that the central Tunisian population is closer to western Mediterranean populations than to eastern ones. Our population was distant from Arabs especially Lebanese and Palestinians. This could be explained by the low influx of Middle East Arabs compared with the settled Berbers [7]. This low genetic contribution of Arabs in Tunisians can be also explained by the low level of mixture between Berber tribes and Arab tribes [9]. In fact, during invasion of Tunisia by the two Arab Tribes (Beni Hilal and Beni Souleim), Berbers were obliged to take refuge in mountains for fear of persecution by the invaders [7, 9, 50]. In addition, the barriers of language, religion and traditions between Arabs and Berbers were additional factors in this low admixture. Nevertheless, although genetic flow from Arabs was low in Tunisia and North Africa, social and cultural effect was very significant evidenced by the adoption of Arab speaking language and Islam religion in all of North African countries. However, Saudi Population was also genetically quite close to central region of Tunisia. This can be probably explained by the by Arab immigration during the Islamic conquests.

Central Tunisians, Blacks and sub-Saharans

Tunisian Blacks are more present in southern Tunisia than in northern Tunisia. Their origins lie mainly in the area extending from Lake of Chad to West Africa [7]. The second origin of Tunisian Blacks was Arab invasion of North African region. However, our genetic analysis of the HLA system has shown that central Tunisians are distant from sub-Saharan populations. This is consistent with the results of other studies which have found that Tunisians, even southern Tunisians, were genetically distant from sub-Saharan Africans [6, 7, 9, 16, 50]. This is probably due to the high endogamy within these black populations, where inter-ethnic marriages were rare because of cultural barriers. This has reduced admixture between sub-Saharan populations and other Tunisians, and limited their genetic contribution to the Tunisian genetic pool [7].

This study does have its limitations, namely the generic typing of HLA alleles, given the lack of resources. However, the strengths of this study are the typing of the 6 HLA-A, -B, -C, -DRB1, -DQA1, DQB1 loci. Indeed, most anthropological studies based on the HLA system have not studied the HLA-C and HLA-DQA1 loci. Moreover, to our knowledge, this is the first anthropological study of the population of the Tunisian center.

In conclusion, according to this study, central Tunisians were very close to other Tunisian populations, including the Berbers. They were also close to North Africans and Iberians. However, central Tunisians were genetically rather distant from eastern Mediterranean populations and sub-Saharans, especially Arabs. This result is somewhat surprising, although it is consistent with other Tunisian studies. Indeed, despite the various Arab invasions of Tunisia and the strong cultural and religious impact, there was no significant Arab contribution to the Tunisian gene pool.

Availability of data and materials

For Data evaluation please contact Dr. Ben Bnina Amène at; amenebenbnina@gmail.com.

References

  1. Barker DJ, Maccari G, Georgiou X, Cooper MA, Flicek P, Robinson J, et al. The IPD-IMGT/HLA Database. Nucleic Acids Res. 2023;51(D1):D1053-60.

    Article  CAS  PubMed  Google Scholar 

  2. Shiina T, Hosomichi K, Inoko H, Kulski JK. The HLA genomic loci map: expression, interaction, diversity and disease. J Hum Genet. 2009;54(1):15–39.

    Article  CAS  PubMed  Google Scholar 

  3. Riley E, Olerup O. HLA polymorphisms and evolution. Immunol Today. 1992;13(9):333–5.

    Article  CAS  PubMed  Google Scholar 

  4. Vina MAF, Hollenbach JA, Lyke KE, Sztein MB, Maiers M, Klitz W, et al. Tracking human migrations by the analysis of the distribution of HLA alleles, lineages and haplotypes in closed and open populations. Philos Trans R Soc B Biol Sci. 2012;367(1590):820–9.

    Article  CAS  Google Scholar 

  5. Abdennaji Guenounou B, Loueslati BY, Buhler S, Hmida S, Ennafaa H, Khodjet-Elkhil H, et al. HLA class II genetic diversity in southern Tunisia and the Mediterranean area. Int J Immunogenet. 2006;33(2):93–103.

    Article  CAS  PubMed  Google Scholar 

  6. Hajjej A, Hmida S, Kaabi H, Dridi A, Jridi A, El Gaaled A, et al. HLA genes in Southern Tunisians (Ghannouch area) and their Relationship with other Mediterraneans. Eur J Med Genet. 2006;49(1):43–56.

    Article  CAS  PubMed  Google Scholar 

  7. Hajjej A, Almawi WY, Hattab L, El-Gaaied A, Hmida S. The investigation of the origin of Southern Tunisians using HLA genes. J Hum Genet. 2017;62(3):419–29.

    Article  CAS  PubMed  Google Scholar 

  8. Ayed K, Ayed-Jendoubi S, Sfar I, Labonne MP, Gebuhrer L. HLA class-I and HLA class-II phenotypic, gene and haplotypic frequencies in Tunisians by using molecular typing data. Tissue Antigens. 2004;64(4):520–32.

    Article  CAS  PubMed  Google Scholar 

  9. Hajjej A, Sellami MH, Kaabi H, Hajjej G, El-Gaaied A, Boukef K, et al. HLA class I and class II polymorphisms in Tunisian Berbers. Ann Hum Biol. 2011;38(2):156–64.

    Article  CAS  PubMed  Google Scholar 

  10. Muro M, Marı́n L, Torı́o A, Moya-Quiles MR, Minguela A, Rosique-Roman J, et al. HLA polymorphism in the murcia population (Spain): in the cradle of the archaeologic Iberians. Hum Immunol. 2001;62(9):910–21.

    Article  CAS  PubMed  Google Scholar 

  11. Sanchez-Velasco P, Gomez-Casado E, Martinez-Laso J, Moscoso J, Zamora J, Lowy E, et al. HLA alleles in isolated populations from North Spain: Origin of the Basques and the ancient Iberians: Sanchez-Velasco et al : HLA in North Spain populations. Tissue Antigens. 2003;61(5):384–92.

    Article  CAS  PubMed  Google Scholar 

  12. Arnaiz-Villena A, Benmamar D, Alvarez M, Diaz-Campos N, Varela P, Gomez-Casado E, Martinez-Laso J. HLA allele and haplotype frequencies in Algerians Relatedness to Spaniards and Basques. Hum Immunol. 1995;43(4):259–68. https://doi.org/10.1016/0198-8859(95)00024-x.

    Article  CAS  PubMed  Google Scholar 

  13. Canossi A, Piancatelli D, Aureli A, Oumhani K, Ozzella G, DelBeato T, et al. Correlation between genetic HLA class I and II polymorphisms and anthropological aspects in the Chaouya population from Morocco (Arabic speaking). Tissue Antigens. 2010;76(3):177–93.

    Article  CAS  PubMed  Google Scholar 

  14. Renquin J, Sanchez-Mazas A, Halle L, Rivalland S, Jaeger G, Mbayo K, et al. HLA class II polymorphism in Aka Pygmies and Bantu Congolese and a reassessment of HLA-DRB1 African diversity. Tissue Antigens. 2001;58(4):211–22.

    Article  CAS  PubMed  Google Scholar 

  15. Rendine S, Ferrero NM, Sacchi N, Costa C, Pollichieni S, Amoroso A. Estimation of human leukocyte antigen class I and class II high-resolution allele and haplotype frequencies in the Italian population and comparison with other European populations. Hum Immunol. 2012;73(4):399–404.

    Article  CAS  PubMed  Google Scholar 

  16. Hajjej A, Almawi WY, Hattab L, Hmida S. Anthropological analysis of Tunisian populations as inferred from HLA class I and class II genetic diversity: A meta-analysis. Immunol Lett. 2017;1(185):12–26.

    Article  Google Scholar 

  17. Stearns PN, Langer WL, editors. The Encyclopedia of World History: Ancient, Medieval, and Modern, Chronologically Arranged. 6th edition. Boston: Houghton Mifflin Harcourt; 2001. p. 1243.

  18. Sellier J, Sellier A. Atlas des peuples d’Orient: Moyen-Orient, Caucase, Asie centrale. La Découverte; 2004. p. 206.

    Google Scholar 

  19. Hajjej A, Hajjej G, Almawi WY, Kaabi H, El-Gaaied A, Hmida S. HLA class I and class II polymorphism in a population from south-eastern Tunisia (Gabes Area): HLA genes in southern Tunisians. Int J Immunogenet. 2011;38(3):191–9.

    Article  CAS  PubMed  Google Scholar 

  20. Mahjoub S, Mehri S, Ghazouani E, Ouarda F, Boussada R, Zaroui A, et al. HLA class II polymorphisms in Tunisian patients with dilated cardiomyopathy. Tissue Antigens. 2010;75(6):679–83.

    Article  CAS  PubMed  Google Scholar 

  21. Hajjej A, Kaabi H, Sellami MH, Dridi A, Jeridi A, El borgi W, et al. The contribution of HLA class I and II alleles and haplotypes to the investigation of the evolutionary history of Tunisians. Tissue Antigens. 2006;68(2):153–62.

    Article  CAS  PubMed  Google Scholar 

  22. Mack S, Tsai Y, Sanchez-Mazas A, Erlich H. 13th International Histocompatibility Workshop Anthropology/Human Genetic Diversity Joint Report - Chapter 3: Anthropology/human genetic diversity population reports. In: Hansen J, editor. Immunobiology of the Human MHC: Proceedings of the 13th International Histocompatibility Workshop and Conference. Seattle: IHWG Press; 2006. p. 580–652.

  23. Arnaiz-Villena A, Elaiwa N, Silvera C, Rostom A, Moscoso J, Gómez-Casado E, et al. The origin of Palestinians and their genetic relatedness with other Mediterranean populations. Hum Immunol. 2001;62(9):889–900.

    CAS  PubMed  Google Scholar 

  24. Elamin NE, Elbashir MI, Elkhidir IM, Elghazali G, Blasczyk R, Horn PA. Identification of four new HLA-Cw alleles in the Sudanese population. Tissue Antigens. 2007;69(3):270–2.

    Article  CAS  PubMed  Google Scholar 

  25. Gonzalez-Galarza FF, McCabe A, Santos EJ, Jones J, Takeshita LY, Ortega-Rivera ND, Del Cid-Pavon GM, Ramsbottom K, Ghattaoraya GS, Alfirevic A, Middleton D, Jones AR. Allele frequency net database (AFND) 2020 update: gold-standard data classification, open access genotype data and new query tools. Nucleic Acid Research. 2020;48:D783-8.

    CAS  Google Scholar 

  26. Jawdat D, Al-Zahrani M, Al-Askar A, Fakhoury H, Uyar FA, Hajeer A. HLA-A, B, C, DRB1 and DQB1 allele and haplotype frequencies in volunteer bone marrow donors from Eastern Region of Saudi Arabia. HLA. 2019;94(1):49–56.

    Article  CAS  PubMed  Google Scholar 

  27. Arnaiz-Villena A, Iliakis P, González-Hevilla M, Longás J, Gómez-Casado E, Sfyridaki K, et al. The origin of Cretan populations as determined by characterization of HLA alleles. Tissue Antigens. 1999;53(3):213–26.

    Article  CAS  PubMed  Google Scholar 

  28. Pascual M, Nieto A, López-Nevot MA, Ramal L, Matarán L, Caballero A, et al. Rheumatoid arthritis in southern Spain: toward elucidation of a unifying role of the HLA class II region in disease predisposition. Arthritis Rheum. 2001;44(2):307–14.

    Article  CAS  PubMed  Google Scholar 

  29. Modiano D, Luoni G, Petrarca V, Sodiomon Sirima B, De Luca M, Simporé J, et al. HLA class I in three West African ethnic groups: genetic distances from sub-Saharan and Caucasoid populations. Tissue Antigens. 2001;57(2):128–37.

    Article  CAS  PubMed  Google Scholar 

  30. Arnaiz-Villena A, Dimitroski K, Pacho A, Moscoso J, Gómez-Casado E, Silvera-Redondo C, et al. HLA genes in Macedonians and the sub-Saharan origin of the Greeks. Tissue Antigens. 2001;57(2):118–27.

    Article  CAS  PubMed  Google Scholar 

  31. Raguénès O, Mercier B, Clèdes J, Whebe B, Férec C. HLA class II typing and idiopathic IgA nephropathy (IgAN): DQB1*0301, a possible marker of unfavorable outcome. Tissue Antigens. 1995;45(4):246–9.

    Article  PubMed  Google Scholar 

  32. Galgani A, Mancino G, Martínez-Labarga C, Cicconi R, Mattei M, Amicosante M, et al. HLA-A, -B and -DRB1 allele frequencies in Cyrenaica population (Libya) and genetic relationships with other populations. Hum Immunol. 2013;74(1):52–9.

    Article  CAS  PubMed  Google Scholar 

  33. Samaha H, Rahal EA, Abou-Jaoude M, Younes M, Dacchache J, Hakime N. HLA class II allele frequencies in the Lebanese population. Mol Immunol. 2003;39(17–18):1079–81.

    Article  CAS  PubMed  Google Scholar 

  34. Elbjeirami WM, Abdel-Rahman F, Hussein AA. Probability of Finding an HLA-Matched Donor in Immediate and Extended Families: The Jordanian Experience The Jordanian Experience. Biol Blood Marrow Transplant 1 févr. 2013;19(2):221–6.

    Article  Google Scholar 

  35. Piancatelli D, Canossi A, Aureli A, Oumhani K, Del Beato T, Di Rocco M, et al. Human leukocyte antigen-A, -B, and -Cw polymorphism in a Berber population from North Morocco using sequence-based typing. Tissue Antigens. 2004;63(2):158–72.

    Article  CAS  PubMed  Google Scholar 

  36. Kijak GH, Walsh AM, Koehler RN, Moqueet N, Eller LA, Eller M, et al. HLA class I allele and haplotype diversity in Ugandans supports the presence of a major east African genetic cluster. Tissue Antigens. 2009;73(3):262–9.

    Article  CAS  PubMed  Google Scholar 

  37. Dempster AP, Laird NM, Rubin DM. Maximum likelihood for incomplete data via algorithm. J R Stat Soc. 1977;39:1–38.

    MathSciNet  Google Scholar 

  38. Excoffier L, Slatkin M. Maximum likelihood estimation of molecular haplotype frequencies in a diploid population. Mol Biol Evol. 1995;12:921–7. https://doi.org/10.1093/oxfordjournals.molbev.a040269.

    Article  CAS  PubMed  Google Scholar 

  39. Schneider S, Kueffer JM, Roessli D, Excoffier L. Arlequin: A Software Environment for the Analysis of Population Genetics Data. Geneva: Genetics and Biometry Lab; 1996.

    Google Scholar 

  40. Imanishi T, Akaza T, Kimura A, Tokunaga K, Gojobori T. Estimation of allele and haplotype frequencies for HLA and complement loci. In: Tsuji K, Aizawa M, Sasazuki T, eds. HLA 1991, Vol. I. Oxford: Oxford University Press. 1992: 76–9. In.

  41. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.

    CAS  PubMed  Google Scholar 

  42. Nei M. Genetic Distance between Populations. Am Nat. 1972;106(949):283–92.

    Article  Google Scholar 

  43. Nei M. Analysis of Gene Diversity in Subdivided Populations. Proc Natl Acad Sci. 1973;70(12):3321–3.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  44. Nei M, Tajima F, Tateno Y. Accuracy of estimated phylogenetic trees from molecular data. J Mol Evol 1 mars. 1983;19(2):153–70.

    Article  CAS  Google Scholar 

  45. Audigier V, Husson F, Josse J. MIMCA: multiple imputation for categorical variables with multiple correspondence analysis. Stat Comput. 2017;27(2):501–18.

    Article  MathSciNet  Google Scholar 

  46. Gómez-Casado E, Del Moral P, Martínez-Laso J, García-Gómez A, Allende L, Silvera-Redondo C, et al. HLA genes in Arabic-speaking Moroccans: close relatedness to Berbers and Iberians: HLA genes in Moroccans. Tissue Antigens. 2000;55(3):239–49.

    Article  PubMed  Google Scholar 

  47. Martinez-Laso J, Juan D, Martinez-Quiles N, Gomez-Casado E, Cuadrado E, Arnaiz-Villena A. The contribution of the HLA-A, -B, -C and -DR, -DQ DNA typing to the study of the origins of Spaniards and Basques. Tissue Antigens. 1995;45(4):237–45.

    Article  CAS  PubMed  Google Scholar 

  48. Grimaldi MC, Crouau-Roy B, Amoros JP, Cambon-Thomsen A, Carcassi C, Orru S, et al. West Mediterranean islands (Corsica, Balearic islands, Sardinia) and the Basque population: contribution of HLA class I molecular markers to their evolutionary history: HLA class I and the west Mediterranean islands. Tissue Antigens. 2001;58(5):281–92.

    Article  CAS  PubMed  Google Scholar 

  49. Clayton J, Lonjou C, editors. Allele and Haplotype frequencies for HLA loci in various ethnic groups. Genetic diversity of HLA. Functional and medical implications. Paris: Charron D; 1997. p. 665–82.

  50. Hajjej A, Almawi WY, Hattab L, El-Gaaied A, Hmida S. HLA Class I and Class II Alleles and Haplotypes Confirm the Berber Origin of the Present Day Tunisian Population Goldberg AC, éditeur. PLOS ONE. 2015;10(8):e0136909.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Arnaiz-Villena A, Karin M, Bendikuze N, Gomez-Casado E, Moscoso J, Silvera C, et al. HLA alleles and haplotypes in the Turkish population: relatedness to Kurds, Armenians and other Mediterraneans: HLA in Turks. Tissue Antigens. 2001;57(4):308–17.

    Article  CAS  PubMed  Google Scholar 

  52. Khil HKE, Marrakchi RT, Loueslati BY, Langaney A, Fellous M, Elgaaied AB. Y chromosome microsatellite variation in three populations of Jerba Island (Tunisia). Ann Hum Genet. 2001;65(3):263–70.

    Article  Google Scholar 

  53. Sanchez-Mazas A. The Berbers of North Africa: Genetic relationships according to HLA and other polymorphisms. Chapter 4. In: In: Arnaiz-Villena A, ed. Prehistoric Iberia: Genetics, Anthropology, and Linguistics. New York: Kluwer-Plenum; 2000. p. 65–77.

    Google Scholar 

  54. Renfrew C. Before Babel: Speculation on the origins of linguistic diversity. Cam Archaeol. 1991;1:3–23.

    Article  Google Scholar 

  55. Fischer, T. in The International Geography (ed. Mill, H. R.) 368–377 (Appleton and Company, New York, NY and London, UK, 1920). In.

Download references

Acknowledgements

We are deeply grateful to all those who contributed to the success of this research project.

Funding

No funding.

Author information

Authors and Affiliations

Authors

Contributions

ABB, wrote the manuscript performed analysis on all samples, and interpreted data. YA  provided statical analysis, helped in results interpretetion  and acted as corresponding author.  YEB, SC  contributed to the writing of the manuscript and discussed the results. NBL, MM, ZM, AM provided HLA Typing and helped in results analysis.  WS  and  DZ edited manuscript, completed data collection. AB responsible for the research design and provided editorial advice.  All Authors have contributed to the drafting and the critical revision of the article. The final version is approved by all authors.

Corresponding author

Correspondence to Amri Yessine.

Ethics declarations

Ethics approval and consent to participate

This study was carried out under the principles of theDeclaration of Helsinki developed by the World Medical Association and approved by the Human Ethics committee of of Sahloul university, Sousse- Tunisia.

Informed consent to participate was obtained from all participants. All data and identities of patients were processed with strict confidentiality.

Consent for publication

Not Applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ben Bnina, A., Yessine, A., El Bahri, Y. et al. Contribution of HLA class I (A, B, C) and HLA class II (DRB1, DQA1, DQB1) alleles and haplotypes in exploring ethnic origin of central Tunisians. BMC Med Genomics 17, 65 (2024). https://doi.org/10.1186/s12920-024-01821-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-024-01821-x

Keywords