Genetic structure of Tunisian sheep breeds as inferred from genome-wide SNP markers

Assessing the status of genetic variability of native sheep breeds could provide important clues for research and policy makers to devise better strategies for the conservation and management of genetic resources. In this study, a genetic investigation of Tunisian sheep breeds using a genome-wide scan of approximately 50,000 SNPs was performed. To reconstruct genetic structure and relationships among four sheep breeds, 40 samples belonging to fat-tailed Barbarine, Queue Fine de l'Ouest, Noire de Thibar and D'Man breeds were genotyped using Illumina Ovine SNP50 BeadChip. Tunisian breeds averaged 96 % polymorphic loci with an expected heterozygosity (He = 0.36). Genetic analysis of relationship between breeds using Bayesian clustering, MDS and Neighbor-Network analysis, and estimation of FST genetic structure, highlighted the genetic differentiation of Noire de Thibar breed from the other local breeds, reflecting the effect of past events of introgression of European gene pool. The Queue Fine de l'Ouest breed showed a genetic heterogeneity and was close to Barbarine and D'Man breeds, as evidenced by MDS and the lowest level of differentiation with Barbarine breed (FST = 1.8 %). The D'Man breed shared a considerable gene flow with the thin-tailed Queue Fine de l'Ouest breed. Possible factors explaining the genetic patterns observed, such as considerable gene flow probably due to anthropogenic activities in the light of population management and conservation programs.