The Tunisian Noire de Thibar sheep breed is a composite breed, recently selected to create animals that are uniformly black in order to avoid skin photosensitization after the ingestion of toxic “hypericum perforatum” weeds, which causes a major economic loss to sheep farmers. We assessed genetic differentiation and estimated marker F_ST using genotyping-by-sequencing (GBS) data in black (Noire de Thibar) and related white-coated (Queue fine de l’ouest) sheep breeds to identify signals of artificial selection. The results revealed the selection signatures within candidate genes related to coat color, which are assumed to be indirectly involved in the mechanism of photosensitization in sheep. The identified genes could provide important information for molecular breeding.

Introduction: The development of reproducible tools for the rapid genotyping of thousands of genetic markers (SNPs) has promoted cross border collaboration in the study of sheep genetic diversity on a global scale.Methods: In this study, we collected a comprehensive dataset of 239 African and Eurasian sheep breeds genotyped at 37,638 filtered SNP markers, with the aim of understanding the genetic structure of 22 North African (NA) sheep breeds within a global context.Results and discussion: We revealed asubstantial enrichment of the gene pool between the north and south shores of the Mediterranean Sea, which corroborates the importance of the maritime route in the history of livestock. The genetic structure of North African breeds mirrors the differential composition of genetic backgrounds following the breed history. Indeed, Maghrebin sheep stocks constitute a geographically and historically coherent unit with any breed-level genetic distinctness among them due to considerable gene flow. We detected a broad east–west pattern describing the most important trend in NA fat-tailed populations, exhibited by the genetic closeness of Egyptian and Libyan fat-tailed sheep to Middle Eastern breeds rather than Maghrebin ones. A Bayesian F_ST scan analysis revealed a set of genes with potentially key adaptive roles in lipid metabolism (BMP2, PDGFD VEGFA, TBX15, and WARS2), coat pigmentation (SOX10, PICK1, PDGFRA, MC1R, and MTIF) and horn morphology RXFP2) in Tunisian sheep. The local ancestry method detected a Merino signature in Tunisian Noire de Thibar sheep near the SULF1gene introgressed by Merino’s European breeds. This study will contribute to the general picture of worldwide sheep genetic diversity.