Heterogeneity in Palaeolithic Population Continuity and Neolithic Expansion in North Africa

Human Y chromosome haplogroup L1-M22 traces Neolithic expansion in West Asia and supports the Elamite and Dravidian connection

West and South Asian populations profoundly influenced Eurasian genetic and cultural diversity. We investigate the genetic history of the Y chromosome haplogroup L1-M22, which, while prevalent in these regions, lacks in-depth study. Robust Bayesian analyses of 165 high-coverage Y chromosomes favor a West Asian origin for L1-M22 ∼20.6 thousand years ago (kya). Moreover, this haplogroup parallels the genome-wide genetic ancestry of hunter-gatherers from the Iranian Plateau and the Caucasus. We characterized two L1-M22 harboring population groups during the Early Holocene. One expanded with the West Asian Neolithic transition. The other moved to South Asia ∼8-6 kya but showed no expansion. This group likely participated in the spread of Dravidian languages. These South Asian L1-M22 lineages expanded ∼4-3 kya, coinciding with the Steppe ancestry introduction. Our findings advance the current understanding of Eurasian historical dynamics, emphasizing L1-M22's West Asian origin, associated population movements, and possible linguistic impacts.

Understanding the genomic heterogeneity of North African Imazighen: from broad to microgeographical perspectives

The strategic location of North Africa has led to cultural and demographic shifts, shaping its genetic structure. Historical migrations brought different genetic components that are evident in present-day North African genomes, along with autochthonous components. The Imazighen (plural of Amazigh) are believed to be the descendants of autochthonous North Africans and speak various Amazigh languages, which belong to the Afro-Asiatic language family. However, the arrival of different human groups, especially during the Arab conquest, caused cultural and linguistic changes in local populations, increasing their heterogeneity. We aim to characterize the genetic structure of the region, using the largest Amazigh dataset to date and other reference samples. Our findings indicate microgeographical genetic heterogeneity among Amazigh populations, modeled by various admixture waves and different effective population sizes. A first admixture wave is detected group-wide around the twelfth century, whereas a second wave appears in some Amazigh groups around the nineteenth century. These events involved populations with higher genetic ancestry from south of the Sahara compared to the current North Africans. A plausible explanation would be the historical trans-Saharan slave trade, which lasted from the Roman times to the nineteenth century. Furthermore, our investigation shows that assortative mating in North Africa has been rare.

The Impact of Recent Demography on Functional Genetic Variation in North African Human Groups

The strategic location of North Africa has made the region the core of a wide range of human demographic events, including migrations, bottlenecks, and admixture processes. This has led to a complex and heterogeneous genetic and cultural landscape, which remains poorly studied compared to other world regions. Whole-exome sequencing is particularly relevant to determine the effects of these demographic events on current-day North Africans' genomes, since it allows to focus on those parts of the genome that are more likely to have direct biomedical consequences. Whole-exome sequencing can also be used to assess the effect of recent demography in functional genetic variation and the efficacy of natural selection, a long-lasting debate. In the present work, we use newly generated whole-exome sequencing and genome-wide array genotypes to investigate the effect of demography in functional variation in 7 North African populations, considering both cultural and demographic differences and with a special focus on Amazigh (plur. Imazighen) groups. We detect genetic differences among populations related to their degree of isolation and the presence of bottlenecks in their recent history. We find differences in the functional part of the genome that suggest a relaxation of purifying selection in the more isolated groups, allowing for an increase of putatively damaging variation. Our results also show a shift in mutational load coinciding with major demographic events in the region and reveal differences within and between cultural and geographic groups.

North and East African mitochondrial genetic variation needs further characterization towards precision medicine

INTRODUCTION

Mitochondria are maternally inherited cell organelles with their own genome, and perform various functions in eukaryotic cells such as energy production and cellular homeostasis. Due to their inheritance and manifold biological roles in health and disease, mitochondrial genetics serves a dual purpose of tracing the history as well as disease susceptibility of human populations across the globe. This work requires a comprehensive catalogue of commonly observed genetic variations in the mitochondrial DNAs for all regions throughout the world. So far, however, certain regions, such as North and East Africa have been understudied.

OBJECTIVES

To address this shortcoming, we have created the most comprehensive quality-controlled North and East African mitochondrial data set to date and use it for characterizing mitochondrial genetic variation in this region.

METHODS

We compiled 11 published cohorts with novel data for mitochondrial genomes from 159 Sudanese individuals. We combined these 641 mitochondrial sequences with sequences from the 1000 Genomes (n = 2504) and the Human Genome Diversity Project (n = 828) and used the tool haplocheck for extensive quality control and detection of in-sample contamination, as well as Nanopore long read sequencing for haplogroup validation of 18 samples.

RESULTS

Using a subset of high-coverage mitochondrial sequences, we predict 15 potentially novel haplogroups in North and East African subjects and observe likely phylogenetic deviations from the established PhyloTree reference for haplogroups L0a1 and L2a1.

CONCLUSION

Our findings demonstrate common hitherto unexplored variants in mitochondrial genomes of North and East Africa that lead to novel phylogenetic relationships between haplogroups present in these regions. These observations call for further in-depth population genetic studies in that region to enable the prospective use of mitochondrial genetic variation for precision medicine.

The genomic history of the indigenous people of the Canary Islands

The indigenous population of the Canary Islands, which colonized the archipelago around the 3rd century CE, provides both a window into the past of North Africa and a unique model to explore the effects of insularity. We generate genome-wide data from 40 individuals from the seven islands, dated between the 3rd-16rd centuries CE. Along with components already present in Moroccan Neolithic populations, the Canarian natives show signatures related to Bronze Age expansions in Eurasia and trans-Saharan migrations. The lack of gene flow between islands and constant or decreasing effective population sizes suggest that populations were isolated. While some island populations maintained relatively high genetic diversity, with the only detected bottleneck coinciding with the colonization time, other islands with fewer natural resources show the effects of insularity and isolation. Finally, consistent genetic differentiation between eastern and western islands points to a more complex colonization process than previously thought.

Whole mitogenomes reveal that NW Africa has acted both as a source and a destination for multiple human movements

Despite being enclosed between the Mediterranean Sea and the Sahara Desert, North Africa has been the scenario of multiple human migrations that have shaped the genetic structure of its present-day populations. Despite its richness, North Africa remains underrepresented in genomic studies. To overcome this, we have sequenced and analyzed 264 mitogenomes from the Algerian Chaoui-speaking Imazighen (a.k.a. Berbers) living in the Aurès region. The maternal genetic composition of the Aurès is similar to Arab populations in the region, dominated by West Eurasian lineages with a moderate presence of M1/U6 North African and L sub-Saharan lineages. When focusing on the time and geographic origin of the North African specific clades within the non-autochthonous haplogroups, different geographical neighboring regions contributed to the North African maternal gene pool during time periods that could be attributed to previously suggested admixture events in the region, since Paleolithic times to recent historical movements such as the Arabization. We have also observed the role of North Africa as a source of geneflow mainly in Southern European regions since Neolithic times. Finally, the present work constitutes an effort to increase the representation of North African populations in genetic databases, which is key to understand their history.

Identifying signatures of positive selection in human populations from North Africa

Because of its location, North Africa (NA) has witnessed continuous demographic movements with an impact on the genomes of present-day human populations. Genomic data describe a complex scenario with varying proportions of at least four main ancestry components: Maghrebi, Middle Eastern-, European-, and West-and-East-African-like. However, the footprint of positive selection in NA has not been studied. Here, we compile genome-wide genotyping data from 190 North Africans and individuals from surrounding populations, investigate for signatures of positive selection using allele frequencies and linkage disequilibrium-based methods and infer ancestry proportions to discern adaptive admixture from post-admixture selection events. Our results show private candidate genes for selection in NA involved in insulin processing (KIF5A), immune function (KIF5A, IL1RN, TLR3), and haemoglobin phenotypes (BCL11A). We also detect signatures of positive selection related to skin pigmentation (SLC24A5, KITLG), and immunity function (IL1R1, CD44, JAK1) shared with European populations and candidate genes associated with haemoglobin phenotypes (HPSE2, HBE1, HBG2), other immune-related (DOCK2) traits, and insulin processing (GLIS3) traits shared with West and East African populations. Finally, the SLC8A1 gene, which codifies for a sodium-calcium exchanger, was the only candidate identified under post-admixture selection in Western NA.

Genetic analysis based on 15 autosomal short tandem repeats (STRs) in the Chaouia population, western center Morocco, and genetic relationships with worldwide populations

The complex demographic history of human populations in North Africa has resulted in a high degree of genetic heterogeneity across the region. However, little is known about the pattern of these genetic variations in its current populations. The present study provides new data on the genetic background of Chaouia, an Arabic-speaking North African population in the western center of Morocco. A random sample of 150 unrelated healthy individuals from Chaouia was assessed using the AmpFLSTR Identifiler kit. The most polymorphic markers were D21S11 and D18S51, with 23 and 22 alleles, respectively. After Bonferroni's correction, two loci (TH01 and D18S51) deviated from Hardy-Weinberg equilibrium. The phylogeny analysis separated North African populations into northeastern and northwestern groups. The Chaouia population was clustered with northwestern Africans. It was the closest to the Berbers of Azrou. The Chaouia shared close genetic affinities with populations from North Africa, the Middle East, and Europe, particularly Iberians, and to a lesser extent with sub-Saharan populations. The pattern of genetic admixture varied across North African populations without a clear correlation between their geographic (northeastern or northwestern) or linguistic identities (Arab or Berber), however, genetic heterogeneity among Berbers was observed. These findings suggest that the diversity observed in North African populations extends geographical and linguistic boundaries. It is further linked to each population's unique and complex demographic history. Human North African population genetics seems to present an intriguing landscape for future studies in the region and its surrounding populations to trace the origins of the genetic heterogeneity observed in these populations.

Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

The genomic analysis of current-day North African populations reveals the existence of trans-Saharan migrations with different origins and dates

The Sahara Desert has acted as a barrier to human gene-flow between the northern and central parts of Africa since its aridification. Nonetheless, some contacts between both sides of the desert have occurred throughout history, mainly driven by commercial activity. Part of this was the infamous trans-Saharan slave trade, which forcedly brought peoples from south of the Sahara to North Africa from Roman times until the nineteenth century. Although historical records exist, the genetic aspects of these trans-Saharan migrations have not been deeply studied. In the present study, we assess the genetic influence of trans-Saharan migrations in current-day North Africa and characterize its amount, geographical origin, and dates. We confirm the heterogeneous and generally low-frequency presence of genomic segments of sub-Saharan origin in present-day North Africans acquired in recent historical times, and we show evidence of at least two admixture events: one dated around the thirteenth-fourteenth centuries CE between North Africans and a Western-sub-Saharan-like source similar to current-day Senegambian populations, and another one dated around the seventeenth century CE involving Tunisians and an Eastern-sub-Saharan-like source related to current-day south-Sudan and Kenyan populations. Time and location coincide with the peak of trans-Saharan slave-trade activity between Western African empires and North African powers, and are also concordant with the possibility of continuous recent south-to-north gene-flow. These findings confirm the trans-Saharan human genetic contacts, providing new and precise evidence about its possible dates and geographical origins, which are pivotal to understanding the genomic composition of an underrepresented region such as North Africa.

Developing CIRdb as a catalog of natural genetic variation in the Canary Islanders

The current inhabitants of the Canary Islands have a unique genetic makeup in the European diversity landscape due to the existence of African footprints from recent admixture events, especially of North African components (> 20%). The underrepresentation of non-Europeans in genetic studies and the sizable North African ancestry, which is nearly absent from all existing catalogs of worldwide genetic diversity, justify the need to develop CIRdb, a population-specific reference catalog of natural genetic variation in the Canary Islanders. Based on array genotyping of the selected unrelated donors and comparisons against available datasets from European, sub-Saharan, and North African populations, we illustrate the intermediate genetic differentiation of Canary Islanders between Europeans and North Africans and the existence of within-population differences that are likely driven by genetic isolation. Here we describe the overall design and the methods that are being implemented to further develop CIRdb. This resource will help to strengthen the implementation of Precision Medicine in this population by contributing to increase the diversity in genetic studies. Among others, this will translate into improved ability to fine map disease genes and simplify the identification of causal variants and estimate the prevalence of unattended Mendelian diseases.

Admixture Mapping of Sepsis in European Individuals With African Ancestries

Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction. Although the ancestral genetic background is a relevant factor for sepsis susceptibility, there is a lack of studies using the genetic singularities of a recently admixed population to identify loci involved in sepsis susceptibility. Here we aimed to discover new sepsis loci by completing the first admixture mapping study of sepsis in Canary Islanders, leveraging their distinctive genetic makeup as a mixture of Europeans and African ancestries. We used a case-control approach and inferred local ancestry blocks from genome-wide data from 113,414 polymorphisms genotyped in 343 patients with sepsis and 410 unrelated controls, all ascertained for grandparental origin in the Canary Islands (Spain). Deviations in local ancestries between cases and controls were tested using logistic regressions, followed by fine-mapping analyses based on imputed genotypes, in silico functional assessments, and gene expression analysis centered on the region of interest. The admixture mapping analysis detected that local European ancestry in a locus spanning 1.2 megabases of chromosome 8p23.1 was associated with sepsis (lowest p = 1.37 × 10⁻⁴; Odds Ratio [OR] = 0.51; 95%CI = 0.40–0.66). Fine-mapping studies prioritized the variant rs13249564 within intron 1 of MFHAS1 gene associated with sepsis (p = 9.94 × 10⁻⁴; OR = 0.65; 95%CI = 0.50–0.84). Functional and gene expression analyses focused on 8p23.1 allowed us to identify alternative genes with possible biological plausibility such as defensins, which are well-known effector molecules of innate immunity. By completing the first admixture mapping study of sepsis, our results revealed a new genetic locus (8p23.1) harboring a number of genes with plausible implications in sepsis susceptibility.

New Insight into the human genetic diversity in North African populations by genotyping of SNPs in DRD3, CSMD1 and NRG1 genes

BACKGROUND

The single nucleotide polymorphisms (SNPs) of the dopamine D3 receptor (DRD3), the CUB and sushi multiple domains 1 (CSMD1) and the neuregulin 1 (NRG1) genes were used to study the genetic diversity and affinity among North African populations and to examine their genetic relationships in worldwide populations.

METHODS

The rs3773678, rs3732783 and rs6280 SNPs of the DRD3 gene located on chromosome 3, the rs10108270 SNP of the CSMD1 gene and the rs383632, rs385396 and rs1462906 SNPs of the NRG1 gene located on chromosome 8 were analysed in 366 individuals from seven North African populations (Libya, Kairouan, Mehdia, Sousse, Kesra, Smar and Kerkennah).

RESULTS

The low values of FST indicated that only 0.27%-1.65% of the genetic variability was due to the differences between the populations. The Kairouan population has the lowest average heterozygosity among the North African populations. Haplotypes composed of the ancestral alleles ACC and ACAT were more frequent in the Kairouan population than in other North African populations. The PCA and the haplotypic analysis showed that the genetic structure of populations in North Africa was closer to that of Europeans, Admixed Americans, South Asians and East Asians. However, analysis of the rs3732783 and rs6280 SNPs revealed that the CT microhaplotype was specific to the North African population.

CONCLUSIONS

The Kairouan population exhibited a relatively low rate of genetic variability. The North African population has undergone significant gene flow but also evolutionary forces that have made it genetically distinct from other populations.

Whole-exome analysis in Tunisian Imazighen and Arabs shows the impact of demography in functional variation

Human populations are genetically affected by their demographic history, which shapes the distribution of their functional genomic variation. However, the genetic impact of recent demography is debated. This issue has been studied in different populations, but never in North Africans, despite their relevant cultural and demographic diversity. In this study we address the question by analyzing new whole-exome sequences from two culturally different Tunisian populations, an isolated Amazigh population and a close non-isolated Arab-speaking population, focusing on the distribution of functional variation. Both populations present clear differences in their variant frequency distribution, in general and for putatively damaging variation. This suggests a relevant effect in the Amazigh population of genetic isolation, drift, and inbreeding, pointing to relaxed purifying selection. We also discover the enrichment in Imazighen of variation associated to specific diseases or phenotypic traits, but the scarce genetic and biomedical data in the region limits further interpretation. Our results show the genomic impact of recent demography and reveal a clear genetic differentiation probably related to culture. These findings highlight the importance of considering cultural and demographic heterogeneity within North Africa when defining population groups, and the need for more data to improve knowledge on the region's health and disease landscape.

Comparison of sequencing data processing pipelines and application to underrepresented African human populations

Background

Population genetic studies of humans make increasing use of high-throughput sequencing in order to capture diversity in an unbiased way. There is an abundance of sequencing technologies, bioinformatic tools and the available genomes are increasing in number. Studies have evaluated and compared some of these technologies and tools, such as the Genome Analysis Toolkit (GATK) and its “Best Practices” bioinformatic pipelines. However, studies often focus on a few genomes of Eurasian origin in order to detect technical issues. We instead surveyed the use of the GATK tools and established a pipeline for processing high coverage full genomes from a diverse set of populations, including Sub-Saharan African groups, in order to reveal challenges from human diversity and stratification.

Results

We surveyed 29 studies using high-throughput sequencing data, and compared their strategies for data pre-processing and variant calling. We found that processing of data is very variable across studies and that the GATK “Best Practices” are seldom followed strictly. We then compared three versions of a GATK pipeline, differing in the inclusion of an indel realignment step and with a modification of the base quality score recalibration step. We applied the pipelines on a diverse set of 28 individuals. We compared the pipelines in terms of count of called variants and overlap of the callsets. We found that the pipelines resulted in similar callsets, in particular after callset filtering. We also ran one of the pipelines on a larger dataset of 179 individuals. We noted that including more individuals at the joint genotyping step resulted in different counts of variants. At the individual level, we observed that the average genome coverage was correlated to the number of variants called.

Conclusions

We conclude that applying the GATK “Best Practices” pipeline, including their recommended reference datasets, to underrepresented populations does not lead to a decrease in the number of called variants compared to alternative pipelines. We recommend to aim for coverage of > 30X if identifying most variants is important, and to work with large sample sizes at the variant calling stage, also for underrepresented individuals and populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04407-x.

Genetic diversity of the North African population revealed by the typing of SNPs in the DRD2/ANKK1 genomic region

The dopamine - related genes, like dopamine D2 receptor (DRD2) gene and ankyrin repeat and kinase domain containing 1 (ANKK1) gene are implicated in neurological functions. Some polymorphisms of the DRD2/ANKK1 locus (TaqIA, TaqIB, TaqID) have been used to study genetic diversity and the evolution of human populations. The present investigation aims to assess the genetic diversity in seven North African populations in order to explore their genetic structure and to compare them to others worldwide populations studied for the same locus. Nine single nucleotide polymorphisms (SNPs) from the DRD2/ANKK1 locus (rs1800497 TaqIA, rs2242592, rs1124492, rs6277, rs6275, rs1079727, rs2002453, rs2234690 and rs1079597 TaqIB) were typed in 366 individuals from seven North African populations: six from Tunisia (Sousse, Smar, Kesra, Kairouan, Mehdia and Kerkennah) and one from Libya. The allelic frequencies of rs2002453 and rs2234690 were higher in the Smar population than in the other North African populations. More, the Smar population showed the lowest average heterozygosity (0.313). The principal component analysis (PCA) showed that the Smar population was clearly separated from others. Furthermore, linkage disequilibrium analysis shown a high linkage disequilibrium in the North African population and essentially in Smar population. Comparison with other world populations has shown that the heterozygosity of North African population was very close to that of the African and European populations. The PCA and the haplotypic analysis suggested the presence of an important Eurasian genetic component for the North African population. These results suggested that the Smar population was isolated from the others North Africans ones by its peculiar genetic structure because of isolation, endogamy and genetic drift. On the other hand, the North African population is characterized by a multi ancestral gene pool from Eurasia and sub-Saharan Africa due to human migration since prehistoric times.

Population history of North Africa based on modern and ancient genomes

Compared with the rest of the African continent, North Africa has provided limited genomic data. Nonetheless, the genetic data available show a complex demographic scenario characterized by extensive admixture and drift. Despite the continuous gene flow from the Middle East, Europe and sub-Saharan Africa, an autochthonous genetic component that dates back to pre-Holocene times is still present in North African groups. The comparison of ancient and modern genomes has evidenced a genetic continuity in the region since Epipaleolithic times. Later population movements, especially the gene flow from the Middle East associated with the Neolithic, have diluted the genetic autochthonous component, creating an east to west gradient. Recent historical movements, such as the Arabization, have also contributed to the genetic landscape observed currently in North Africa and have culturally transformed the region. Genome analyses have not shown evidence of a clear correlation between cultural and genetic diversity in North Africa, as there is no genetic pattern of differentiation between Tamazight (i.e. Berber) and Arab speakers as a whole. Besides the gene flow received from neighboring areas, the analysis of North African genomes has shown that the region has also acted as a source of gene flow since ancient times. As a result of the genetic uniqueness of North African groups and the lack of available data, there is an urgent need for the study of genetic variation in the region and its implications in health and disease.

Recent Common Origin, Reduced Population Size, and Marked Admixture Have Shaped European Roma Genomes

The Roma Diaspora-traditionally known as Gypsies-remains among the least explored population migratory events in historical times. It involved the migration of Roma ancestors out-of-India through the plateaus of Western Asia ultimately reaching Europe. The demographic effects of the Diaspora-bottlenecks, endogamy, and gene flow-might have left marked molecular traces in the Roma genomes. Here, we analyze the whole-genome sequence of 46 Roma individuals pertaining to four migrant groups in six European countries. Our analyses revealed a strong, early founder effect followed by a drastic reduction of ∼44% in effective population size. The Roma common ancestors split from the Punjabi population, from Northwest India, some generations before the Diaspora started, <2,000 years ago. The initial bottleneck and subsequent endogamy are revealed by the occurrence of extensive runs of homozygosity and identity-by-descent segments in all Roma populations. Furthermore, we provide evidence of gene flow from Armenian and Anatolian groups in present-day Roma, although the primary contribution to Roma gene pool comes from non-Roma Europeans, which accounts for >50% of their genomes. The linguistic and historical differentiation of Roma in migrant groups is confirmed by the differential proportion, but not a differential source, of European admixture in the Roma groups, which shows a westward cline. In the present study, we found that despite the strong admixture Roma had in their diaspora, the signature of the initial bottleneck and the subsequent endogamy is still present in Roma genomes.

An integrated personal and population-based Egyptian genome reference

A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

Berbers and Arabs: Tracing the genetic diversity and history of Southern Tunisia through genome wide analysis

OBJECTIVES

Tunisia has been a crossroads for people from Africa, Europe, and the Middle East since prehistoric times. At present, it is inhabited by two main ethnic groups, Arabs and Berbers, and several minorities. This study aims to advance knowledge regarding their genetic structure using new population samplings and a genome-wide approach.

MATERIALS AND METHODS

We investigated genomic variation, estimated ancestry components and dated admixture events in three Berber and two Arab populations from Southern Tunisia, mining a dataset including Middle Eastern, sub-Saharan, and European populations.

RESULTS

Differences in the proportion of North African, Arabian, and European ancestries and the varying impact of admixture and isolation determined significant heterogeneity in the genetic structure of Southern Tunisian populations. Admixture time estimates show a multilayer pattern of admixture events, involving both ethno-linguistic groups, which started around the mid XI century and lasted for nearly five centuries.

DISCUSSION

Our study provides evidence that the relationships between genetic and cultural diversity of old and new inhabitants of North Africa in southern Tunisia follow different patterns. The Berbers seem to have preserved a significant part of their common genomic heritage despite Islamization, Arab cultural influence, and linguistic diversity. Compared to Morocco and Algeria, southern Tunisian Arabs have retained a higher level of Arabian ancestry. This is more evident in the semi-nomad R'Baya, who have kept their original Bedouin lifestyle, than in the population from Douz, who have undergone multiple events of stratification and admixture.