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

History and genetic diversity of African sheep: Contrasting phenotypic and genomic diversity

Domesticated sheep have adapted to contrasting and extreme environments and continue to play important roles in local community-based economies throughout Africa. Here we review the Neolithic migrations of thin-tailed sheep and the later introductions of fat-tailed sheep into eastern Africa. According to contemporary pictorial evidence, the latter occurred in Egypt not before the Ptolemaic period (305-25 BCE). We further describe the more recent history of sheep in Egypt, the Maghreb, west and central Africa, central-east Africa, and southern Africa. We also present a comprehensive molecular survey based on the analysis of 50 K SNP genotypes for 59 African breeds contributed by several laboratories. We propose that gene flow and import of fat-tailed sheep have partially overwritten the diversity profile created by the initial migration. We found a genetic contrast between sheep north and south of the Sahara and a west-east contrast of thin- and fat-tailed sheep. There is no close relationship between African and central and east Asian fat-tailed breeds, whereas we observe within Africa only a modest effect of tail types on breed relationships.

Insight into the current genomic diversity, conservation status and population structure of Tunisian Barbarine sheep breed

Local livestock breeds play a crucial role in global biodiversity, connecting natural and human-influenced environments and contributing significantly to ecosystem services. While commercial breeds dominate industrial systems, local livestock breeds in developing countries, like Barbarine sheep in Tunisia, are vital for food security and community maintenance. The Tunisian Barbarine sheep, known for its adaptability and distinctive fat-tailed morphology, faces challenges due to historical crossbreeding. In this study, the Illumina Ovine SNP50K BeadChip array was used to perform a genome-wide characterization of Tunisian Barbarine sheep to investigate its genetic diversity, the genome structure, and the relationship within the context of Mediterranean breeds. The results show moderate genetic diversity and low inbreeding. Runs of Homozygosity analysis find genomic regions linked to important traits, including fat tail characteristics. Genomic relationship analysis shows proximity to Algerian thin-tailed breeds, suggesting crossbreeding impacts. Admixture analysis reveals unique genetic patterns, emphasizing the Tunisian Barbarine's identity within the Mediterranean context and its closeness to African breeds. Current results represent a starting point for the creation of monitoring and conservation plans. In summary, despite genetic dilution due to crossbreeding, the identification of genomic regions offers crucial insights for conservation. The study confirms the importance of preserving unique genetic characteristics of local breeds, particularly in the face of ongoing crossbreeding practices and environmental challenges. These findings contribute valuable insights for the sustainable management of this unique genetic reservoir, supporting local economies and preserving sheep species biodiversity.

Refining the genomic profiles of North African sheep breeds through meta-analysis of worldwide genomic SNP data

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 FST 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.

Marker-assisted introgression of fecundity mutation into Barbarine sheep breed under low input production system

Fecundity gene introgression has become an economical method to improve sheep prolificacy in developing countries. The FecX^Bar and FecGH mutations are variants of the BMP15 and GDF9 sheep genes, respectively, identified in Tunisian Barbarine sheep prolific line "W-INRAT," created through a prolificacy-based selection program. The first mutation caused increased prolificacy in heterozygous state and sterility in homozygous ewes. The aim of this work was to increase the number of effective carriers by the introgression of fecundity mutation into non-carrier conventional flocks based on a marker-assisted breeding program. The genotyping was carried out to follow up on the segregation of prolificacy mutation. The conventional ewes, raised in state farms of Tunisia and inseminated by "W-INRAT" rams, yielded 100% carriers females at heterozygous state. These females were selected to be inseminated with conventional rams and produced carriers' lambs with 66.7% males and 33.3% females. The prolific males will be kept in breeding centers to disseminate the fecundity mutation in commercial flocks. This approach facilitates the dissemination of fecundity genes and contributes to livelihood improvement in communities raising Barbarine sheep.

Efficiency of genotyping by sequencing in inferring genomic relatedness and molecular insights into fat tail selection in Tunisian sheep

In developing countries, the use of simple and cost-efficient molecular technology is crucial for genetic characterization of local animal resources and better development of conservation strategies. The genotyping by sequencing (GBS) technique, also called restriction enzyme- reduced representational sequencing, is an efficient, cost-effective method for simultaneous discovery and genotyping of many markers. In the present study, we applied a two-enzyme GBS protocol (PstI/MspI) to discover and genotype SNP markers among 197 Tunisian sheep samples. A total of 100 333 bi-allelic SNPs were discovered and genotyped with an SNP call rate of 0.69 and mean sample depth 3.33. The genomic relatedness between 183 samples grouped the samples perfectly to their populations and pointed out a high genetic relatedness of inbred subpopulation reflecting the current adopted reproductive strategies. The genome-wide association study contrasting fat vs. thin-tailed breeds detected 41 significant variants including a peak positioned on OAR20. We identified FOXC1, GMDS, VEGFA, OXCT1, VRTN and BMP2 as the most promising for sheep tail-type trait. The GBS data have been useful to assess the population structure and improve our understanding of the genomic architecture of distinctive characteristics shaped by selection pressure in local sheep breeds. This study successfully investigates a cost-efficient method to discover genotypes, assign populations and understand insights into sheep adaptation to arid area. GBS could be of potential utility in livestock species in developing/emerging countries.

Haplotype structure of MSTN, IGF1, and BMP2 genes in Tunisian goats (Capra hircus) and their association with morphometric traits

The aim of the present study was to evaluate variability and haplotype structure of twenty-eight single nucleotide polymorphisms (SNPs) at myostatin (MSTN), insulin-like growth factor 1 (IGF1), and bone morphogenetic protein 2 (BMP2) genes. Association between the polymorphic SNPs and morphometric traits was performed on a population of 263 Tunisian goats. The SNPs analyzed were all polymorphic (except one), and the three genes had different haplotype structures. Significant association of SNPs at MSTN with head length was highlighted in the Tunisian goats. The variability at IGF1 gene was associated with body length, ear length, tail length, and chest depth. For BMP2, significant association was revealed with chest depth. Significant association was also detected between linkage disequilibrium (LD) block 2 at IGF1 with body length. These findings might play a potential role in gene-assisted programs.

Genomic Population Structure of the Main Historical Genetic Lines of Spanish Merino Sheep

Simple Summary

Historical documentation shows that the Spanish Merino sheep was selected over many centuries due to the quality of wool, following which it was used to originate all other Merino breeds around the world, mainly by crossbreeding with local breeds. Today, the historical genetic lines that originated the Spanish Merino are still preserved in several closed herds in which they have been bred for nearly 200 years, maintaining their original genetic purity. Our study demonstrates, using a genomic approach, the exceptional genetic richness and variability of these lines, which are clearly differentiated from modern Merino breeds, and must therefore be protected to safeguard the large genetic pool they represent.

Abstract

According to historiographical documentation, the Romans first began to select Merino sheep in the Iberian Peninsula during the first century, with the aim of obtaining a breed appreciated for the quality of its wool. This process continued locally during the Middle Ages, when Spanish sheep were protected, and their export to foreign countries was banned. It was during the 16th century when individual Merino sheep were allowed to spread around the world to be used to improve the wool quality of local breeds. However, the wool crisis of the 1960s shifted the selection criteria of the Merino breed towards meat production at the expenses of wool. Consequently, individuals that display the genetic and phenotypic characteristics of those sheep originally bred in the kingdom of Spain in the Middle Ages are extremely difficult to find in commercial herds. In this study, we characterized the genetic basis of 403 individuals from the main historical Spanish Merino genetic lines (Granda, Hidalgo, Lopez-Montenegro, Maeso, Donoso and Egea), which were bred in isolation over the last 200 years, using a genomic approach based on genotyping data from the Axiom™ Ovine 50K SNP Genotyping Array. Our analysis included measuring population structure, genomic differentiation indexes, runs of homozygosity (ROH) patterns, and an analysis of molecular variance (AMOVA). The results showed large genetic differences between the historical lines, even though they belong to the same breed. In addition, ROH analysis showed differences due to increased inbreeding among the ancient generations compared with the modern Merino lines, confirming the breed’s ancestral and closed origin. However, our results also showed a high variability and richness within the Spanish historical Merino lines from a genetic viewpoint. This fact, together with their great ability to produce high-quality wool, suggests that ancestral Merino lines from Spain should be considered a valuable genetic population to be maintained as a resource for the improvement of wool-producing sheep breeds all around the world.

Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern

Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome‐wide single‐nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model‐based clustering and Neighbour‐Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (H_e) and contemporary effective population size (cNe) were reported in Siciliano (H_e = 0.323, cNe = 397), and the lowest were reported in Purosangue Orientale Siciliano (H_e = 0.277, cNe = 10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome‐wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.

Genome variation in tick infestation and cryptic divergence in Tunisian indigenous sheep

Background

Ticks are obligate haematophagous ectoparasites considered second to mosquitos as vectors and reservoirs of multiple pathogens of global concern. Individual variation in tick infestation has been reported in indigenous sheep, but its genetic control remains unknown.

Results

Here, we report 397 genome-wide signatures of selection overlapping 991 genes from the analysis, using ROH, LR-GWAS, XP-EHH, and F_ST, of 600 K SNP genotype data from 165 Tunisian sheep showing high and low levels of tick infestations and piroplasm infections. We consider 45 signatures that are detected by consensus results of at least two methods as high-confidence selection regions. These spanned 104 genes which included immune system function genes, solute carriers and chemokine receptor. One region spanned STX5, that has been associated with tick resistance in cattle, implicating it as a prime candidate in sheep. We also observed RAB6B and TF in a high confidence candidate region that has been associated with growth traits suggesting natural selection is enhancing growth and developmental stability under tick challenge. The analysis also revealed fine-scale genome structure indicative of cryptic divergence in Tunisian sheep.

Conclusions

Our findings provide a genomic reference that can enhance the understanding of the genetic architecture of tick resistance and cryptic divergence in indigenous African sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08321-1.

How have sheep breeds differentiated from each other in Morocco? Genetic structure and geographical distribution patterns

Background

Based on the relatively homogeneous origin of the sheep breeds in Morocco that originate mainly from Iberia, it is highly relevant to address the question of how these very diverse sheep populations differentiated from each other. The Mountains of the High Atlas and Middle Atlas are expected to constitute North–South and West–East geographical barriers, respectively, which could have shaped the history of the differentiation of sheep breeds. The aim of this study was to test this hypothesis by considering the genetic structure and the spatial distribution of five major breeds (Sardi, Timahdite, Beni Guil, Boujaad and D’man) and one minor breed (Blanche de Montagne), by analysing the mtDNA control region, using 30 individuals per breed.

Results

Phylogenetic and network analyses did not indicate any clear separation among the studied breeds and discriminant component principal analysis showed some overlap between them, which indicates a common genetic background. The calculated pairwise F_ST values and Nei’s genetic distances revealed that most breeds showed a moderate genetic differentiation. The lowest and highest degrees of differentiation were retrieved in the Beni Guil and Boujaad breeds, respectively. Analysis of molecular variance (AMOVA) indicated that more than 95% of the genetic diversity occurs within individuals, while between- and within-population variabilities represent only 1.332% and 2.881%, respectively. Isolation-by-distance, spatial Principal Component Analysis (sPCA), and spatial AMOVA analyses evidenced clear examples of geographical structuration among the breeds, both between and within breeds. However, several enigmatic relationships remain, which suggest the occurrence of complex events leading to breed differentiation.

Conclusions

The approaches used here resulted in a convergent view on the hypothetic events that could have led to the progressive differentiation between the Moroccan breeds. The major split seems to be linked to the West–East barrier of the Middle Atlas, whereas the influence of the High Atlas is less obvious and incompletely resolved. The study of additional breeds that have settled near the High Atlas should clarify the relationships between the breeds of the West part of the country, in spite of their small population size.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00679-2.

Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area

Under a climate change perspective, the genetic make-up of local livestock breeds showing adaptive traits should be explored and preserved as a priority. We used genotype data from the ovine 50 k Illumina BeadChip for assessing breed autozygosity based on runs of homozygosity (ROH) and fine-scale genetic structure and for detecting genomic regions under selection in 63 Tunisian sheep samples. The average genomic inbreeding coefficients based on ROH were estimated at 0.017, 0.021, and 0.024 for Barbarine (BAR, n = 26), Noire de Thibar (NDT, n = 23), and Queue fine de l'Ouest (QFO, n = 14) breeds, respectively. The genomic relationships among individuals based on identity by state (IBS) distance matrix highlighted a recent introgression of QFO into the BAR and a genetic differentiation of NDT samples, possibly explained by past introgression of European gene pools. Genome-wide scan for ROH across breeds and within the BAR sample set identified an outstanding signal on chromosome 13 (46.58-49.61 Mbp). These results were confirmed using F_ST index, differentiating fat vs. thin-tailed individuals. Candidate genes under selection pressure (CDS2, PROKR1, and BMP2) were associated to lipid storage and probably preferentially selected in fat-tailed BAR animals. Our findings suggest paying more attention to preserve the genetic integrity and adaptive alleles of local sheep breeds.