# Genetic homogenization of indigenous sheep breeds in Northwest Africa

Ibrahim Belabdi, Abdessamad Ouhrouch, Mohamed Lafri, Semir Bechir Suheil Gaouar, Elena Ciani, Ahmed Redha Benali, Hakim Ould Ouelhadj, Abdelmajid Haddioui, Francois Pompanon, Véronique V. Blanquet, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Ibrahim Belabdi, Abdessamad Ouhrouch, Mohamed Lafri, Semir Bechir Suheil Gaouar, Elena Ciani, et al.. Genetic homogenization of indigenous sheep breeds in Northwest Africa. Scientific Reports, 2019, 9 (1), pp.1-13. 10.1038/s41598-019-44137-y. hal-02627209

HAL Id: hal-02627209

https://hal.inrae.fr/hal-02627209

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# Genetic homogenization of indigenous sheep breeds in Northwest Africa

Received: 18 April 2018

Accepted: 7 May 2019

Published online: 28 May 2019

Ibrahim Belabdi<fnref n="1" />,², Abdessamad Ouhrouch³,⁷, Mohamed Lafri<fnref n="1" />,², Semir Bechir Suheil Gaouar⁴, Elena Ciani⁵, Ahmed Redha Benali<fnref n="1" />,², Hakim Ould Ouelhadj⁶, Abdelmajid Haddioui⁷, François Pompanon⁹, Véronique Blanquet⁹, Dominique Taurisson-Mouret¹⁰, Sahraoui Harkat<fnref n="1" />,², Johannes A. Lenstra<fnref n="1" />,¹¹, Badr Benjelloun³,⁸ & Anne Da Silva⁹

Northwest-African sheep represent an ideal case-study for assessing the potential impact of genetic homogenization as a threat to the future of traditional breeds that are adapted to local conditions. We studied ten Algerian and Moroccan breeds of sheep, including three transboundary breeds, distributed over a large part of the Maghreb region, which represents a geographically and historically coherent unit. Our analysis of the dataset that involved carrying out Genome-wide SNP genotyping, revealed a high level of homogenization (ADMIXTURE, NetView, fineSTRUCTURE and IBD segments analyses), in such a way that some breeds from different origins appeared genetically undistinguished: by grouping the eight most admixed populations, we obtained a mean global FST value of 0.0024. The sPCA analysis revealed that the major part of Morocco and the Northern part of Algeria were affected by the phenomenon, including most of the breeds considered. Unsupervised cross-breeding with the popular Ouled-Djellal breed was identified as a proximate cause of this homogenization. The issue of transboundary breeds was investigated, and the Hamra breed in particular was examined via ROH fragments analysis. Genetic diversity was considered in the light of historical archives and anthropological works. All of these elements taken together suggest that homogenization as a factor affecting the Maghrebin sheep stock, has been particularly significant over the last few decades, although this process probably started much earlier. In particular, we have identified the policies set by the French administration during the colonial period of the region's history as a causal factor that probably contributed significantly to this process. The genetic homogenization that we have observed calls into question the integrity of the farm animal genomic resources represented by these local breeds, whose conservation is of critical importance to the future of the livestock sector.

Local breeds are defined as those that have been in a country "for a sufficient time to be genetically adapted to one or more traditional production systems or environments in the country"<fnref n="1" />. Due to these adaptations, local breeds represent valuable genetic resources, especially in the context of climate change². However, a large proportion of them are at serious risk of extinction³. One of the main threats is uncontrolled cross-breeding (i.e., breeding not carried out within the framework of selection plans) with more productive local or imported breeds being favoured by farmers, who face increasing economic pressure. These practices, assessed by FAO to establish

<footnote n="1">¹Science Veterinary Institute, University of Blida, BP 270, Blida, 09000, Algeria. ²Laboratory of Biotechnology related to Animal Reproduction (LBRA), University of Blida, BP 270, Blida, 09000, Algeria. ³National Institute of Agronomic Research (INRA Maroc), Regional Centre of Agronomic Research, Beni-Mellal, Morocco. ⁴Department of Biology, Aboubakr Belkaid Tlemcen University, laboratory of Physiopathologie et biochimie de la Nutrition (PpBioNut), Tlemcen, Algeria. ⁵Department of Biosciences, Biotechnologies and Biopharmaceutics, University of Bari, Bari, Italy. ⁶Institut technique des Elevages Saida, Saida, Algeria. ⁷Laboratoire de Biotechnologies et Valorisation des Ressources Phytogénétiques (LBVRP), Université Sultan Moulay Slimane, BéniMellal, Morocco. ⁸Univ. Grenoble Alpes, Univ. Savoie Mont-Blanc, CNRS, LECA, 38000, Grenoble, France. ⁹Univ. Limoges, INRA, PEREINE EA7500, USC1061 GAMAA, F-87000, Limoges, France. ¹⁰CNRS, UMR 5815, Dynamiques du droit, Université de Montpellier, Montpellier, France. ¹¹Utrecht University, Faculty of Veterinary Medicine, Utrecht, Netherlands. Ibrahim Belabdi, Abdessamad Ouhrouch, Badr Benjelloun and Anne Da Silva contributed equally. Correspondence and requests for materials should be addressed to A.D.S. (email: anne.blondeau@unilim.fr)</footnote>
risk status, lead to genetic dilution, *i.e.* overwhelmingly high numbers of non-local genes lowering the relative frequency of native genes in the population. They favor genetic homogenization via the loss of rare or specific variants, and induce disruption of co-adapted gene complexes, which jeopardizes the integrity and viability of the local populations⁴,⁵.

The stock of North African domestic sheep (*Ovis aries*) provides an ideal case-study for assessing the problem of genetic homogenization among local breeds and its impact on farm animal genetic resources. This stock is 7000 years old⁶, and contains a remarkable diversity of populations, which have been maintained under traditional farming systems over millennia⁷. These systems have enabled the emergence of breeds locally adapted. In the Maghreb, the diversity of territories (mountains, steppes, desert, oasis, etc.) and climates, in combination with the history of the region's people, have shaped a rich livestock heritage.

In recent decades it has been observed that in Algeria there is a strong preference among farmers for the Ouled-Djellal breed, due to its larger conformation, and this has led to the intensification of anarchic crossing practices between local breeds. The breeders hope through such practices to increase their productivity⁸⁻¹⁰. This results in a strong admixture between Ouled-Djellal, Rembi, Taàdmit, and Berber breeds¹¹,¹². For the Berber breed, a primitive breed that is thought to be the most ancient breed of the Maghreb¹³, absorption by an Arab breed-type had already begun in 1900. Indeed, in an observation dating from that time, Couput¹⁴ mentions that only isolated populations, generally located in the mountains, remained preserved. Today, DNA markers do not detect any differentiation between Berber and Ouled-Djellal, Rembi or Taàdmit breeds, while the Sidaoun breed (raised by the Tuareg people), the D'Man breed (mostly encountered in oases), the Tazegzawth breed (from the Kabyle mountains) and the Hamra breed (called Beni-Guil in Morocco) all appear to be more genetically preserved¹¹,¹².

Our study of the genetic diversity of Maghrebin sheep relies on breeds from Morocco and Algeria, two countries that represent a substantial part of the Maghreb area. The stabilization of national borders in the Maghreb region started in the 18th century. For this reason an analysis of the processes shaping these livestock, whose origins date as far back as 7000 years BP, would be better and more coherently approached under a broader geographic scale (*i.e.* the diversity of sheep breeds in the Maghreb should be considered at a scale above the national level). The aims of the study were to assess the genetic structuration of the Maghrebin sheep stock. Particular attention was paid to the case of the transboundary (*i.e.* Hamra/Beni-Guil, D'Man and Ouled-Djellal) breeds, in order to see how their history and genetic characteristics can provide insights regarding the evolution of the Maghrebin sheep stock.

## Results

A total of 180 sheep (93 from Algeria and 87 from Morocco) from ten breeds (three trans-boundary breeds, five Algerian and two Moroccan breeds, Table S1) were analyzed through a panel of 36,376 SNPs after having merged the Moroccan and Algerian data.

The average $F_{IS}$ value was rather low (mean = -0.044, s.d. = 0.06), with 129 individuals (*i.e.* 71.7%) showing negative values. A total of 22 individuals showed a pronounced deficit in heterozygosity with $F_{IS}$ values > 0.1; this deficit was particularly apparent among the Barbarine (BRBA) and Moroccan D'Man (DMNM) breeds, with 40% and 43.3% of individuals displaying such values, respectively.

Calculation of pair-wise $F_{ST}$ (Table 1) was used to explore the genetic relationships among breeds. Two groups were detected: (i) the first group was composed of Ouled-Jellal subpopulations from Algeria and Morocco (OLDM, OLDA), D'Man from Morocco (DMNM), Hamra from Algeria (HAMA), Beni-Guil (BIGM, also called Hamra in Algeria), Berber (BERA), Sardi (SRDM), Timahdite (TMHM), and Rembi (RMBA); this group had a mean global $F_{ST}$ value of 0.0024 [0.0021-0.0026]$_{95\%}$ and showed very low values of pair-wise $F_{ST}$ (<0.006), with nine of these values not significantly different from zero (see confidence intervals). (ii) in the second group, which included the five Algerian remaining breeds (Sidaoun, D'Man, Hamra from pilot farms, Barbarine and Tazegzawth), the pair-wise values of $F_{ST}$ ranged from 0.022 (Barbarine, BRBA/Rembi, RMBA) to 0.113 (Algerian preserved Hamra, HAMAP/Tazegzawth, TZGA).

Considering the case of transboundary breeds, we found that: (i) Beni-Guil (BIGM, called Hamra in Algeria) was not differentiated from Hamra (HAMA) sampled in private Algerian farms ($F_{ST}$ value not significantly different from zero), but was clearly differentiated from the Algerian Hamra population preserved in pilot farms (HAMAP) ($F_{ST}$ = 0.051); (ii) D'Man populations from Morocco and Algeria (DMNM and DMNA) showed pair-wise $F_{ST}$ of 0.026; (iii) Ouled-Djellal populations from Morocco and Algeria (OLDM and OLDA) showed pair-wise $F_{ST}$ close to zero ($F_{ST}$ = 0.002). The NeighborNet graph based on the $F_{ST}$ genetic distances (Fig. S1) allows visualization of these results.

The ADMIXTURE analysis revealed the same pattern, whereby the sheep were characterized by two groups (Fig. 1): In the “non-homogeneous” group, the Sidaoun (SDNA) breed was differentiated from all the other breeds for $K=2$ (*i.e.* the lowest cross-validation (CV) error) and showed genetic overlap with Algerian D'Man (DMNA). Increasing $K$ from 3 to 6, assigned clusters to the preserved Hamra population (HAMAP), Tazegzawth (TZGA), Barbarine (BRBA) and Algerian D'Man (DMNA). The second, or “homogeneous”, group was characterized by an important phenomenon of admixture, and included all of the remaining breeds. At $K=7$, a few D'Man individuals from Morocco (DMNM) showed genetic distinctness. For $K > 7$, intra-breed structuration was also apparent for the SDNA, HAMAP, BRBA and DMNA breeds.

NetView was used to represent the genetic structure of the dataset along a gradient of $k$ values (from 2 to 100), allowing for the investigation of a range of structures, from fine scale population structures, revealed by a small value of $k$, to large scale population structures, revealed by a large value of $k$ (see Fig. 2). As individuals showing $F_{IS}$ value > 0.1 were removed from the analysis (see)¹⁵, we were led to remove the Barbarine (BRBA) breed, which was then only represented by three individuals. The simulated data network was based on a minimum spanning tree. At fine scale (Fig. 2a, $k=15$) the topology of the network showed the existence of four homogeneous clusters
<table><caption>**Table 1.** Pair-wise F<sub>ST</sub> between Algerian and Moroccan breeds. *In bold values with two zeros after the decimal point. BERA = Berber; BRBA = Barbarine; DMNA = D'Man from Algeria; DMNM = D'Man from Morocco; HAMA = Hamra; BIGM = Beni-Guil; HAMAP = Hamra from pilot farms; OLDA = Ouled-Djellal from Algeria; OLDM = Ouled-Djellal from Morocco; RMBA = Rembi; SDNA = Sidaoun; SRDM = Sardi; TMHM = Timahdite; TZGA = Tazegzawth.</caption><thead><tr><th></th><th>SDNA</th><th>BERA</th><th>RMBA</th><th>BRBA</th><th>DMNA</th><th>OLDA</th><th>HAMA</th><th>HAMAP</th><th>TZGA</th><th>BIGM</th><th>SRDM</th><th>OLDM</th><th>DMNM</th></tr></thead><tbody><tr><td>SDNA</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BERA</td><td>[0.036-0.041]<br>0.038</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>RMBA</td><td>[0.036-0.040]<br>0.038</td><td>[-0.000-0.003]<br>0.001</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BRBA</td><td>[0.059-0.062]<br>0.060</td><td>[0.021-0.026]<br>0.023</td><td>[0.021-0.024]<br>0.022</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>DMNA</td><td>[0.041-0.044]<br>0.042</td><td>[0.032-0.036]<br>0.034</td><td>[0.032-0.037]<br>0.035</td><td>[0.052-0.058]<br>0.055</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>OLDA</td><td>[0.037-0.039]<br>0.038</td><td>[0.002-0.006]<br>0.004</td><td>[-0.002-0.002]<br>0.000</td><td>[0.022-0.025]<br>0.023</td><td>[0.034-0.037]<br>0.035</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HAMA</td><td>[0.032-0.036]<br>0.034</td><td>[0.001-0.004]<br>0.002</td><td>[0.001-0.003]<br>0.002</td><td>[0.022-0.026]<br>0.024</td><td>[0.031-0.034]<br>0.032</td><td>[0.001-0.002]<br>0.001</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HAMAP</td><td>[0.077-0.082]<br>0.080</td><td>[0.049-0.053]<br>0.051</td><td>[0.050-0.054]<br>0.052</td><td>[0.071-0.075]<br>0.073</td><td>[0.080-0.085]<br>0.082</td><td>[0.050-0.054]<br>0.052</td><td>[0.043-0.046]<br>0.045</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>TZGA</td><td>[0.096-0.100]<br>0.098</td><td>[0.064-0.068]<br>0.066</td><td>[0.064-0.068]<br>0.066</td><td>[0.085-0.091]<br>0.087</td><td>[0.097-0.102]<br>0.098</td><td>[0.064-0.068]<br>0.066</td><td>[0.060-0.064]<br>0.062</td><td>[0.111-0.118]<br>0.113</td><td></td><td></td><td></td><td></td></tr><tr><td>BIGM</td><td>[0.032-0.036]<br>0.034</td><td>[-0.001-0.002]<br>0.001</td><td>[-0.001-0.003]<br>0.002</td><td>[0.021-0.025]<br>0.022</td><td>[0.030-0.035]<br>0.032</td><td>[0.001-0.004]<br>0.002</td><td>[-0.002-0.001]<br>-0.000</td><td>[0.048-0.053]<br>0.051</td><td>[0.064-0.069]<br>0.067</td><td></td><td></td><td></td></tr><tr><td>SRDM</td><td>[0.032-0.035]<br>0.034</td><td>[0.003-0.005]<br>0.004</td><td>[0.004-0.006]<br>0.005</td><td>[0.027-0.031]<br>0.029</td><td>[0.031-0.033]<br>0.029</td><td>[0.005-0.007]<br>0.006</td><td>[0.002-0.004]<br>0.003</td><td>[0.048-0.052]<br>0.050</td><td>[0.064-0.067]<br>0.065</td><td>[0.001-0.003]<br>0.002</td><td></td><td></td></tr><tr><td>OLDM</td><td>[0.034-0.037]<br>0.036</td><td>[0.000-0.004]<br>0.002</td><td>[-0.000-0.002]<br>0.001</td><td>[0.021-0.025]<br>0.022</td><td>[0.030-0.033]<br>0.031</td><td>[0.000-0.003]<br>0.002</td><td>[-0.000-0.002]<br>0.001</td><td>[0.048-0.052]<br>0.050</td><td>[0.062-0.067]<br>0.065</td><td>[-0.002-0.001]<br>0.000</td><td>[0.003-0.005]<br>0.004</td><td></td></tr><tr><td>DMNM</td><td>[0.026-0.029]<br>0.028</td><td>[0.001-0.003]<br>0.002</td><td>[0.001-0.004]<br>0.002</td><td>[0.022-0.025]<br>0.023</td><td>[0.025-0.028]<br>0.026</td><td>[0.002-0.004]<br>0.003</td><td>[0.002-0.003]<br>0.002</td><td>[0.047-0.051]<br>0.049</td><td>[0.061-0.065]<br>0.063</td><td>[-0.002-0.000]<br>-0.001</td><td>[0.003-0.005]<br>0.004</td><td>[0.000-0.003]<br>0.001</td></tr><tr><td>TMHM</td><td>[0.033-0.036]<br>0.035</td><td>[0.002-0.005]<br>0.003</td><td>[0.003-0.006]<br>0.004</td><td>[0.026-0.029]<br>0.027</td><td>[0.031-0.033]<br>0.032</td><td>[0.004-0.006]<br>0.005</td><td>[0.002-0.003]<br>0.002</td><td>[0.045-0.049]<br>0.047</td><td>[0.062-0.065]<br>0.063</td><td>[0.000-0.002]<br>0.001</td><td>[0.002-0.003]<br>0.002</td><td>[0.002-0.004]<br>0.003</td><td>[0.002-0.003]<br>0.002</td></tr></tbody></table>

corresponding to Sidaoun (SDNA), Algerian D'Man (DMNA), Tazegzawth (TZGA) and Algerian preserved Hamra (HAMAP). It was interesting to note that HAMAP was closely linked to one HAMA individual (Algerian Hamra sampled in private farm). Most of the other individuals were highly interconnected. At k = 25 (Fig. 2b), Algerian D'Man (DMNA) was connected (i) with the central and admixed core, via a Moroccan D'Man (DMNM) individual, (ii) and on another side, with the Sidaoun (SDNA) cluster. At k = 75 (Fig. 2c), clear connections appeared between Algerian D'Man (DMNA) and Sidaoun (SDNA); moreover only Tazegzawth (TZGA) individuals remained unlinked to the central core. The connection of TZGA to the central core was only observed at k = 100 (Fig. 2d) and via a Berber (BERA) individual (i.e. a breed reared in Kabyle mountains as TZGA). The tree based on Allele-Sharing Distance (ASD) distance between individuals (Fig. S2) also showed genetic relatedness between a number of individuals belonging to Sardi and D'Man breeds (i.e. Moroccan breeds).

The characteristics of ROH, were used as source of information about demographic history. Given the limited size of some samples, the purpose of the ROH and IBD segment analyses (see next paragraph) was not to study each breed according to the values obtained, but to assess the extent to which strong patterns, i.e. clear splits between groups of breeds, were also highlighted by these methods. Indeed, plotting the number of ROHs found for each individual genome against ROH total size (Fig. 3) revealed a clear separation between two groups: in the first group (Fig. 3a) of eight breeds (i.e. a group containing the DMNA breed plus all breeds, except DMNM and SRDM, that were identified by the previous analyses of F<sub>ST</sub>, Admixture and NetView as belonging to the “homo- geneous” group), most individuals carried few and short ROHs. The 50 K SNP panel did not detect any ROH in Beni-Guil (BIGM, called Hamra in Algeria), Hamra sampled in private Algerian farms (HAMA), Algerian Ouled-Djellal population (OLDA) and Rembi (RMBA). Such profiles are characteristic of genetically diverse and most probably admixed populations<sup>16,17</sup>. The second group (Fig. 3b), included Barbarine (BRBA), Tazegzawth (TZGA), the preserved Algerian Hamra population (HAMAP), Sidaoun (SDNA) (i.e. all of the breeds identified by the previous analyses, F<sub>ST</sub>, Admixture and NetView, as belonging to the “non-homogeneous” group, except
<figure><img src="image_3.png" /></figure>

**Figure 1.** Bayesian clustering performed with ADMIXTURE software on Moroccan and Algerian sheep. K = number of clusters; HAMAP = Hamra form pilot farms; HAMA = Hamra; BIGM = Beni-Guil; SDNA = Sidaoun; BRBA = Barbarine; BERA = Berber; OLDA = Ouled-Djellal from Algeria; OLDM = Ouled-Djellal from Morocco; DMNA = D'Man from Algeria; DMNM = D'Man from Morocco; SRDM = Sardi; TMHM = Timahdite; RMBA = Rembi; TZGA = Tazegzawth.

<figure><img src="image_4.png" /></figure>

**Figure 2.** Mutual nearest-neighbour graphs obtained from NetView and considering the following k values: k = 15 (a), k = 25 (b), k = 75 (c) and k = 100 (d). Color shades represent different sheep breeds. HAMAP = Hamra form pilot farms; HAMA = Hamra; SDNA = Sidaoun; BERA = Berber; OLDM = Ouled-Djellal from Morocco; DMNA = D'Man from Algeria; DMNM = D'Man from Morocco; SRDM = Sardi; TMHM = Timahdite; TZGA = Tazegzawth.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Runs of homozygosity (ROH) identified for Moroccan and Algerian sheep breeds. The number of ROHs found for each individual genome (y-axis) is plotted against ROH total size (i.e. the number of Kb covered by ROH in each genome, x-axis); (a) breeds in which ROH length was close to zero, (b) others. HAMAP = Hamra form pilot farms; HAMA = Hamra; BIGM = Beni-Guil; SDNA = Sidaoun; BRBA = Barbarine; BERA = Berber; OLDA = Ouled-Djellal from Algeria; OLDM = Ouled-Djellal from Morocco; DMNA = D'Man from Algeria; DMNM = D'Man from Morocco; SRDM = Sardi; TMHM = Timahdite; RMBA = Rembi; TZGA = Tazegzawth.</figcaption></figure>

DMNA), and also Sardi (SRDM), and the Moroccan D'Man population (DMNM). This second group was characterized by a higher number of ROHs. Tazegzawth (TZGA) showed several long ROH fragments, a pattern characteristic of recent inbreeding, resulting in $F_{ROH} > 0.1$ for four individuals out of the six considered (Table S2). Similarly, long ROH fragments were detected in Barbarine and $F_{ROH}$ showed very high values (i.e. 0.22 and 0.34) for two of the five individuals. The Hamra subpopulation, preserved in Algerian pilot farms (HAMAP), received special attention through the analysis of the ROH distribution according to their size (Fig. 4). For this breed, ROH showed a rather uniform distribution for the short and intermediate length categories. Such a pattern is characteristic of an isolated population that is based on a source population of quite large size, whereas a different pattern showing an overrepresentation of short ROHs would have been the signature of a strong population bottleneck$^{18,19}$. Moreover, the low abundance of long ROHs suggests that management practices have been used to avoid inbreeding for HAMAP, which was confirmed by the $F_{ROH}$ values (Table S2). We further noted recent inbreeding for the Moroccan D'Man population (DMNM), indeed ROH size distribution was relatively uniform for classes between 5 and 25 Mb whereas a peak revealed very high values for extremely long ROHs ($>30$ Mb). The $F_{ROH}$ consideration indicated that 43% of individuals showed a high value (i.e. $F_{ROH} > 0.1$) (Table S2).

IBD segments analysis (Table 2) revealed that all the breeds except Tazegzawth (TAZ), the preserved Hamra population (HAMAP) and the Algerian D'Man population (DMNA), were characterized by absence or quasi absence of IBD segments. It is worth noting that Sidaoun (SDNA) also contained almost no IBD segments. For comparison purposes, similar analyses were conducted on Italian local breeds (i.e. breeds standardized during the
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Classification of ROHs in seven categories (x-axis) according to size (from 1 to 5 Mb to more than 30 Mb) and mean sum of ROHs (y-axis, measured in kilobases) within each ROH category and averaged per breed. We analyzed sheep breed from Algeria and Morocco. Results were not presented for breed with sample size ≤6. HAMAP = Hamra form pilot farms; HAMA = Hamra; SDNA = Sidaoun; OLDM = Ouled-Djellal from Morocco; DMNM = D’Man from Morocco; SRDM = Sardi; TMHM = Timahdite.</figcaption></figure>

<table><caption>**Table 2.** Number and lengths of IBD segments. Nb., number of individuals considered (i.e. individuals showing F<sub>IS</sub> values > 0.1); SD, Standard Deviation; *values given as an indication, since the limited size of the samples biases the average estimates.</caption><thead><tr><th>Breeds</th><th>Nb. of animals</th><th>Nb. of IBD segments</th><th>Nb. of IBD segment/ Nb. animals</th><th>% IBD segment intra-breed</th><th>% IBD segment inter-breed</th><th>Mean length (SD) Mb</th></tr></thead><tbody><tr><td>Hamra from Algerian pilot farms</td><td>10</td><td>456</td><td>45.6</td><td>100</td><td>0</td><td>14.11 (11.14)</td></tr><tr><td>Tazegzawth</td><td>5</td><td>245</td><td>49*</td><td>100</td><td>0</td><td>14.92 (9.87)*</td></tr><tr><td>D’Man from Algeria</td><td>5</td><td>19</td><td>3.8*</td><td>100</td><td>0</td><td>14.06 (8.50)*</td></tr><tr><td>Ouled-Djellal from Algeria</td><td>6</td><td>1</td><td>0.17*</td><td>100</td><td>0</td><td>4.99 (—)*</td></tr><tr><td>Sidaoun</td><td>36</td><td>4</td><td>0.11</td><td>100</td><td>0</td><td>9.91 (5.19)</td></tr><tr><td>Sardi</td><td>24</td><td>1</td><td>0.04</td><td>100</td><td>0</td><td>8.45 (—)</td></tr><tr><td>Timahdite</td><td>16</td><td>7</td><td>0.44</td><td>29</td><td>71</td><td>7.15 (1.78)</td></tr><tr><td>Berber</td><td>6</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Beni-Guil</td><td>6</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Barbarine</td><td>3</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>D’Man from Morocco</td><td>17</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Hamra</td><td>10</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Ouled-Djellal from Morocco</td><td>8</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Rembi</td><td>6</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr></tbody></table>

20th centuries) and South-West Asian breeds (see Supplementary Text 1 and 2) highlighting a clear relationship between the occurrence of admixture and a decrease, in number and length, of IBD segments.

The co-ancestry heatmap (Fig. 5) obtained with CHROMOPAINTER/fineSTRUCTURE presents the number of shared haplotypes (“chunks”) between Algerian and Moroccan individuals. As recommended by Van Dorp *et al.*²⁰ individuals showing high values of F<sub>IS</sub> (>0.1) were excluded from the analysis. The individuals belonging to the Sidaoun (SDNA) and Algerian D’Man (DMNA) breeds clustered and were separated from the other populations on the basis of shared haplotypes. The Tazegzawth (TZGA) and Algerian preserved Hamra (HAMAP) breeds (and one HAMA sheep, *i.e.* the same individual who was identified in the Netview analysis, Fig. 2) were also clearly individualized as animals showing the highest number of genomic shared chunks in the dataset, whereas the remaining breeds appeared clearly admixed. When considering Sidaoun (SDNA) in more detail, it turns out that the number of chunks shared within the breed was significant but that these shared haplotypes were also characterized by their relatively small size (matrix of chunk length not shown), which accounts for the fact that the previous IBD analysis did not detect them (Table 2).

In the spatial PCA analyses using Delaunay triangulation as connection networks, the consideration of eigenvalues suggested the possible occurrence of a spatial pattern, as the three first positive scores were distinguished from other eigenvalues (Fig. 6b). The global Monte Carlo test (10,000 iterations) indicated significant global spatial structure (p-value = 0.009). The Mantel test showed no significant correspondence between geographic and genetic distances (p-value = 0.07). Indeed, while the north constituted a largely homogeneous block, individuals separated by more than 2000 km appeared very close from a genetic point of view. The eigenvectors of the three first global scores plotted according to the geographical coordinates indicated a clear homogeneity of the
<figure><img src="image_7.png" /><figcaption>**Figure 5.** FineSTRUCTURE clustering for Algerian and Moroccan sheep breeds. The color of each bin in the matrix indicates the number of “genomic chunks” copied from a donor (column) to a recipient individual (row). HAMAP = Hamra form pilot farms; HAMA = Hamra; BIGM = Beni-Guil; SDNA = Sidaoun; BRBA = Barbarine; BERA = Berber; OLDA = Ouled-Djellal from Algeria; OLDM = Ouled-Djellal from Morocco; DMNA = D’Man from Algeria; DMNM = D’Man from Morocco; SRDM = Sardi; TMHM = Timahdite; RMBA = Rembi; TZGA = Tazegzawth.</figcaption></figure>

Northern stock (i.e., Moroccan and Northern-Algeria sheep) with the Hamra preserved in pilot farms (HAMAP) and Tazegzawth (TZGA) as exceptions (Fig. 6a). The first coordinate highlighted the contrast of Sidaoun (SDNA) and Barbarine (BRBA) with the other breeds (Fig. 6c) and the second coordinate highlighted the separate position of the Tazegzawth (TZGA) (Fig. 6d).

Discussion

This study provides the opportunity to assess genetic diversity in several traditional sheep breeds in the Maghreb, a region which constitutes a coherent unit from geographical and historical perspective.

All of the analyses revealed a high level of genetic homogenization for most breeds, irrespective of their origin. More specifically, the breeds affected by this phenomenon were Ouled-Djellal (populations from Morocco and Algeria), Timahdite, Sardi, Hamra/Beni-Guil (populations from Morocco and Algerian private farms, but not those from Algerian state farms), D’Man (the population from Morocco but not the population from Algeria), Rembi and Berber. We observed a strong spatial pattern with two major areas identified. The first one was the Northern area, corresponding to the major part of Morocco and the Northern part of Algeria delimited to the South by steppe. This area appeared homogeneous for all the breeds previously mentioned except for two differentiated breeds, i.e. the Hamra breed preserved in pilot farms and the Tazegzawth breed from the Kabyle Mountains. The second area identified was the Southern area, including breeds from the Sahara and Oriental Great Erg regions, i.e. the Sidaoun and Barbarine breeds. This second area appeared clearly differentiated from the Northern one.

Chromopainter analysis showed that breeds in the homogeneous group, clustered together and shared few genetic chunks (haplotypes). The characteristics of ROHs, used as sources of information about the demographic history of the breeds¹⁸,¹⁹,²¹ gave rise to consistent conclusions. Indeed, most of the breeds mentioned above were characterized by a very limited number of ROHs with very small sizes. In particular, four of them, Rembi, Ouled-Djellal from Algeria and Beni-Guil/Hamra from Morocco and private Algerian farms, were characterized by the quasi absence of ROH, clearly indicating admixture¹⁶,¹⁷. Moreover, these breeds were characterized by extremely few or no IBD segments, what was shown to also be a signature of genetic homogenization (see Supplementary Text 1 and 2). Indeed, it seems likely that strong and continuous admixture, in small ruminants
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Spatial PCA (sPCA) Analysis of Algerian and Moroccan sheep breeds (a); eigenvalues for each global and local axis (positive and negative eigenvalues indicate global and local structures, respectively) (b); first global score of sPCA (c); second global score of sPCA (d). HAMAP = Hamra form pilot farms; SDNA = Sidaoun; BRBA = Barbarine; DMNA = D'Man from Algeria.</figcaption></figure>

(characterized by relatively short generation time) breaks up contiguous genomic blocks resulting in very small haplotypes.

Among the breeds identified as being genetically differentiated from other breeds, our analysis appeared to show that the Tazegzawth breed, reared by Kabyle people in geographically isolated conditions, is threatened by inbreeding (see ROH analysis). These results were not surprising given that the breed accounts for less than 0.02% (i.e. a few hundred heads) of the present-day Algerian sheep stock. In contrast, the Sidaoun breed, reared by Touareg people, did not appear to be affected by inbreeding. This breed was genetically differentiated from the other studied breeds (as shown by e.g. sPCA, Admixture, NetView or FST analyses), but would have experienced a significant level of admixture with breeds/populations not considered in this study. Indeed, the Chromopainter and IBD segments analyses reveal a very particular profile for this breed with quite a large number of intra-breed shared segments, but with such segments characterized as being relatively small in size. Further investigations will be needed, for example focusing on the genetic relationship between this breed and the Malian sheep stock, which is thought to be the cradle of the Sidaoun breed. Another southern breed, the Barbarine, a transboundary breed related to the Tunisian Barbarine (i.e. the main breed in Tunisia), showed clear genetic distinctness. This fat-tail breed originates from West-Asia, which may explain its genetic distance from all other breeds considered in this study. Analyses of ROH (distribution and FROH) raised concerns about possible inbreeding problems, however, only examination of a larger sample could lead to robust conclusions. Such conclusions would not be surprising if we consider the size reductions suffered by the Algerian Barbarine population in recent decades. Indeed, whilst the fat tail of the Algerian Barbarine sheep means that the animal is remarkably well adapted to arid environments, this same trait hinders the animal's marketability²², its slaughter quality being considered as diminished because of it (the amount of caudal fat can be 2 to 4 kg, which is not marketable). The Algerian Barbarine breed has therefore lost the interest of breeders, who have abandoned it in favor of the Ouled-Djellal breed²³. The effect of this preference among Algerian breeders is so significant that the analysis of information contained in the DAGRIS database (http://dagris.ilri.cgiar.org/) now indicates that the Barbarine breed is in danger of being "absorbed by Ouled-Djellal" in Algeria.

Focusing on transboundary breeds, it is interesting to note that the genetic specificity of Algerian D'Man, sampled from state farms that have been tasked with breed preservation, was not observed in Moroccan samples, which were collected in private farms. In the Moroccan population, many very long stretches of ROH revealed strong levels of recent inbreeding for a number of individuals. This breed, which is highly adapted to oasis environments and is encountered in Morocco and Algeria, has a probable Moroccan origin. Indeed, its Algerian name is "mouton Filali" suggesting that the cradle of the breed could well be Tafilalet (in Morocco). The French agronomist G. Toutain²⁴, lamented, even as far back as 1979, the lack of due attention given to this breed, which he qualified as a "genetic treasure" shaped by traditional practices over centuries. In particular, when making
those comments Toutain regretted that the Merino breed was being considered as good candidate to introduce in the oases. An investigation of archive records might provide a way to trace the history of this breed and thereby understand the process that led to the breed's present-day genetic dilution.

In our study another transboundary breed, Hamra/Beni-Guil, was sampled from different sources, with samples taken both in Morocco and in Algeria, and in different types of farms (pilot farms that have been tasked with preserving the breed, and private farms). The study of this breed through these different populations has been particularly instructive. Among the three sample sites, only sheep reared in the Algerian pilot farms appeared not to be genetically diluted (with also one individual sampled in an Algerian private farm). The main state farm in charge of the preservation of the Hamra breed (Algeria), gathered its initial herd at the end of the 1980s. Since then, a few carefully selected progenitors from private farms have occasionally been introduced into the initial core, over the years. Moreover, a strict breeding program was devised to prevent inbreeding. All the genetic results were in accordance with the management practices reported for the pilot farm as they showed (i) a pattern of short to intermediate ROH sizes, which is a pattern specific to that of a population that is isolated but that is nonetheless based on a source population of significant size¹⁸,¹⁹ (i.e. not a pattern that indicates a strong bottleneck effect), (ii) a low abundance of long ROHs suggesting management practices avoiding inbreeding. Taking all of these elements into account, we can assume that Hamra's dilution has probably intensified over the past 30 years. Indeed the pilot state farms were established at that time, yet they have visibly allowed for the preservation of at least part of the genetic specificity of the breed, even if the founding effect experienced when the farm was created undoubtedly led to a reduction in genetic diversity. This hypothesis is largely supported by various anthropological studies. These have shown that over the last few decades there has been a significant intensification of crossing practices between Ouled-Djellal and this breed in private farms: Bechchari et al.²⁵ described an "invasion phenomenon" characterized by a growing number of Ouled-Djellal in Beni-Guil herds. According to interviews conducted by Bechchari et al.²⁵, it seems that the expansion of the Ouled-Djellal breed is mainly the result of the sedentarization and intensification of agriculture, i.e., breeders, who are under increasing pressure to improve productivity, favor Ouled-Djellal for its market value. These results were confirmed by Brisebarre²⁶ who, moreover, reported massive inflows of Ouled-Djellal originating from Algeria. A figure of 200,000 to 500,000 heads of Ouled-Djellal (depending on the interviewed source) was estimated have been illegally entering Morocco each year²⁶.

If we consult the historical archives, we learn that it has long been recognised that the Hamra breed is an ancient one, belonging to a stock that was passed down within a tribe from one generation to the next, and that this practice likely originated with the Zénète people. The breed's geographic isolation is also related to the specific geological features of the Oran plateau where it lives²⁷. The strong local adaptation of this breed, with respect to climate, vegetation and disease resistance, was emphasised as far back as 1891, in a report commissioned by the General Governor of Algeria²⁸, which concluded that introduction of other breeds to this territory and the crosses induced would jeopardize the survival capacity of the Hamra breed.

Moreover, it may be noted that this breed, called “Le petit Oranais” by the settlers, was particularly appreciated in France for its organoleptic qualities and conformation. This breed of sheep was therefore particularly affected by the sheep export activities in Algeria, with the animals transported live by boat. At the end of the 19th century, it is estimated that this trade, from Algeria to France²⁹, involved more than one million sheep heads per year (all breeds considered). Today, the main pilot farm responsible for the breed's preservation keeps 600 sheep or 1% of the total number of Hamra heads in Algeria. Since 1980 the Moroccan government has implemented countrywide a so-called “sheep strategy”, under which most of the country has been divided into cradles of the main breeds. In these areas, breed phenotypic standardization was highly favored and measures were taken to limit gene flows³⁰. Our study did not identify any preserved areas from a genetic point of view. However, it should be stated that our sampling design in Morocco did not focus on the cradle of breeds, contrary to the Algerian sampling.

If we refer to the historical archives, it appears that the colonial period in the region's history probably contributed significantly to the expansion of the Ouled-Djellal breed. Indeed, at the end of the 19th century, studies commissioned by the colonial government concluded that the “white breed” of the steppes was the breed on which all efforts should focus³¹,³². The aim of the settlers was to drastically increase the number of sheep by promoting this breed through various means (steppe development and irrigation works, creation of “model” farms, etc.), with guidelines explicitly recommending the elimination of breeds considered useless, of poor quality, or judged “abnormal” with regard to their conformation (the Berber and the Barbarine breeds were the ones that were mainly affected by such guidelines).

The diversity pattern of the Maghrebin sheep appears very similar to that revealed by Ouchene-Khelifi et al.³³ in their study looking at the Maghrebin goat in the same area. Indeed, the results of this goat study also showed very weak genetic structuration among Moroccan breeds (i.e. Black, Draa, Nord) and the most widespread of the Algerian breeds (i.e. Arabia). Despite the breeds being from different geographical origins (i.e. Black from Morocco is Berber type and Arabia from Algeria is Sahelian type), they appeared quite indistinct from each other from a genetic perspective (i.e. mean FST was of 0.005 [0.004–0.006]IC95% considering the four breeds together). Two breeds were distinguished from this main homogeneous group: the Algerian ancestral Kabyle and M'Zabite breeds, which have both been reared in isolated conditions, due to social and political factors, since ancient times. In order to assess how the patterns observed in Maghrebin small ruminants were influenced by migration processes, outward from centers of domestication, we examined the results provided by Colli et al.³⁴ for goats. Considering the 12 African countries for which at least two local breeds were genotyped by country (i.e. a total of 37 breeds, distributed from north to south across the continent), it can be seen that of the 703 pair-wise FST calculated, only 12 showed values lower than 0.01 and none of them had a value lower than 0.006. These results postulate for a clear structuring of African goat breeds, whatever the migration routes considered, with occasional cases of more intense gene flows, generally occurring between two breeds. The data currently available do not allow for a study looking specifically at breeds of sheep in Africa. However, Gaouar et al.¹² observed that, out of a total of 4460 pair-wise FST values, and considering Algerian and worldwide breeds (from Kijas et al.³⁵) all together,
instances of particularly low pair-wise FST values were only found for Algerian breeds and Australian Merino populations (i.e. strains of Merino that are not breeds as such). Thus, to our knowledge, as things currently stand, it seems unlikely that the genetic homogenization identified is an effect that results from ancestral migration processes (given that the FST values are outside the norm and affect no less than seven breeds, in a geographical area covering more than several hundred thousand km²).

Maghrebin breeds are representative of traditional breeds more generally, that, worldwide, are a cornerstone of the rural economies on which smallholder farmers and herders rely. These Maghrebin breeds were not standardized during the 20th century, unlike European breeds, and therefore they continue to show intra-breed phenotypic diversity. They could be qualified as landrace breeds, representing an initial stage of breed development, with a genetic specificity mainly resulting from the combination of founder effects, isolation and natural selection. As emphasized by Sponenberg et al.³⁶ isolation, resulting from geographic factors, but also from social practices (e.g. tribe organization) and political contexts, represents a key factor for the development and the preservation of the genetic integrity of landrace breeds. In the Maghreb, such isolation has been disturbed, inducing the genetic homogenization observed in most Maghrebin sheep breeds. In recent decades, motorized transport has certainly played an important role in increasing long-distance gene flow.

In conclusion, our study highlights a strong homogenization of the Maghreb sheep stock affecting the major part of Morocco and most of Northern Algeria. We still lack information regarding the scale at which this phenomenon of homogenization might be occurring worldwide. However, the observation of such a phenomenon clearly calls into question the integrity of farm animal genomic resources represented by these local breeds, whose conservation is of critical importance for the development of livestock breeding in the context of global changes².

## Material and Methods

**Ethics statement.** Samples were collected by veterinarians during routine blood sampling, for medical care or follow up, such that no ethical authorization was required. All the samples and data processed in our study were obtained with the breeders and breeding organizations' consent.

**Breeds.** We studied ten Algerian and Moroccan breeds representing most of the local breeds for this area. Three transboundary breeds were sampled in both countries: Ouled-Djellal, Hamra (Algerian name)/Beni-Guil (Moroccan name) and D'Man. Most of the breeds were “wooled thin-tailed sheep”, but the D'Man breed belongs to the “mixed hair-wool sheep” group, Sidaoun to the “hairy sheep” group and Barbarine was the only “fat-tailed” breed in our study. Algerian breeds were described by Gaouar et al.¹¹,¹². Supplementary Table 1 provides essential information about the breeds studied.

The geographical distribution of the breeds based on their cradle of origin, as described in the 1990s by Chellig³⁷ for Algeria and by Boujenane³⁸ for Morocco, is presented in Fig. S3. This distribution reflects a pattern described by French settlers at the beginning of the 19th century, responsible for drawing up an inventory of the country's wealth. In addition, the consultation of several historical archives¹⁴,²⁴,²⁷,²⁸,³¹,³²,³⁹⁻⁴¹ made it possible to record essential information about these landrace breeds, which have been extensively described by zootechnicians and agronomists, including: (i) phenotypic traits specific to the different breeds (neckline shape, coat colour, wool quality, etc.); (ii) detailed descriptions of the areas occupied by the breeds (vegetation, topography, etc.), and of the geographical factors implying a more or less extensive isolation of the populations; (iii) history and breeding practices of the tribes that have maintained the traditional breeds; (iv) crossbreeding experiments undertaken by the settlers with imported breeds but also between local breeds. A general observation was the strong consistency of the genetic material within the indigenous breeds as evidenced by the heritability of the characteristic traits.

**Sampling design.** The whole dataset included 46 Algerian samples from Gaouar et al.¹², 87 Moroccan samples provided by the NextGen project (Grant Agreement no. 244356), plus 47 Algerian samples genotyped specially for this study (i.e. 14 additional samples of the Hamra breed and 33 additional samples of the Sidaoun breed). If we consider the case of the Hamra breed: previously available genotypes¹² (n = 6) were only representative of one Algerian pilot farm (Aïn El Hadjar). With the further sampling, 3 samples from two other Algerian state farms and 11 samples from private farms were added to the dataset. The sampling design was developed in such a way that it was possible to study the trajectories of different populations of the same breed. The sheep were sampled, as much as possible, from different flocks in order to limit relatedness among individuals and optimize representativeness of breeds. In cases where several samples were obtained from the same farm, the choice was based on the pedigree documentation (if available) and/or genealogical information provided by the breeder. For a few breeds only few samples were available, but this did not affect the robustness of the results, as almost all analytical methods used were based on the consideration of individuals (the very few calculations based on average values were reported only as being of indicative value in the case of small samples). Moreover, we used the Weir & Cockerham FST estimator for quantifying differentiation between breeds, since in contrast to the Reynolds' distance it is not significantly biased by a small sample size provided that, as in our case, the number of markers used is large⁴²,⁴³.

**Genotyping and SNP quality control.** Blood samples were cryo-preserved until DNA extraction and analysis were carried out. Genomic DNA was purified from whole blood by protease K digestion and a salting-out procedure⁴⁴. All Algerian animals were genotyped for 54,241 SNPs, using the Illumina® Ovine SNP50 BeadChip (Illumina, Inc.), following standard operating procedures recommended by the manufacturer. Genotyping was performed by the Van Haeringen Laboratorium (Wageningen, The Netherlands). Moroccan SNP data were extracted from the WGS variation using the Ovine SNP50 BeadChip coordinates of SNPs on the OAR v3.1
reference genome assembly using Vcftools⁴⁵. High concordance between 12× re-sequencing data and the Ovine SNP50 BeadChip was demonstrated in sheep by Benjelloun et al.⁴⁶. Moroccan and Algerian datasets were merged, and SNPs and animals were pruned with PLINK v1.07⁴⁷ using the following filtering out parameters: (i) SNP call rate ≤ 97%, (ii) SNP minor allele frequency (MAF) ≤ 1%, (iii) animals displaying ≥ 10% of missing genotypes.

**Data analysis.** Genetic diversity, runs of homozygosity and segments of identity-by-descent (IBD). Inbreeding coefficient ($F_{IS}$) per individual was assessed in each breed using PLINK. Runs of homozygosity (ROHs) were identified with PLINK using a sliding windows approach. The fraction of 50 SNP-windows that were almost completely homozygous (i.e. allowing: one heterozygous SNP, five missing genotypes and a maximum gap of 1000 kb) was calculated following recommendations of Manunza et al.¹⁹. Indeed, it has been shown that this approach (i.e. considering windows of SNPs that are not completely homozygous) is more powerful in detecting truly autozygous segments⁴⁸. To reduce the occurrence of spurious runs of homozygosity (ROHs), the minimum number of SNPs required to constitute a ROH was estimated using the method proposed by Lencz et al.⁴⁹. The inbreeding coefficient based on ROH ($F_{ROH}$)⁵⁰ was calculated by dividing total ROH length per individual by total autosomal SNP coverage (2.44 Gb).

BEAGLE 4.1⁵¹ was used for detection of IBD segments. The ibdtrim parameter was set to 40. Segments with a LOD score < 4 and a length shorter than 0.5 cM were excluded. Individuals showing $F_{IS}$ values > 0.1 were excluded from analyses.

*Genetic structure.* The extent of population subdivision was examined by calculating pair-wise $F_{ST}$ values of Weir and Cockerham⁵² and the associated 95% confidence intervals using Genetic Data Analysis (GDA) software⁵³.

We used the network-based approach implemented in NetView v.1.1¹⁵,⁵⁴ available in R. This visualization tool allows for the analysis of complex genetic structure via the construction of population networks through mutual k-nearest neighbours thresholds applied to genome-wide SNPs. PLINK was used to construct the initial ASD matrix (–distance-matrix, i.e. 1 - Identity by Similarity). As recommended by Neuditschko et al.¹⁵, the dataset was first subjected to quality control (i.e. minor allele frequency (<0.01) and significant deviation from Hardy–Weinberg equilibrium (P < 0.001)). We explored the data in a range of k (i.e. parameter determining the number of mutual nearest neighbours) from 2 to 100 that was determined following recommendation of Neuditschko et al.¹⁵ with the “selection plot” function.

A tree based on ASD distance between individuals and a NeighborNet graph based on $F_{ST}$ genetic distance were constructed using Splitstree⁵⁵.

ADMIXTURE software⁵⁶ was used to investigate the relationship between sheep breeds. ADMIXTURE was run for $K=2$ through $K=14$, and ten independent runs were performed for each value of $K$. The program CLUMPAK⁵⁷ available at http://clumpak.tau.ac.il, was used to analyze the multiple independent runs at a single $K$ and visualize the results.

The approach of fineSTRUCTURE⁵⁸ based on haplotypes was also used. The dataset was filtered and phased using SHAPEIT ver. 2⁵⁹. We used CHROMOPAINTER⁵⁸ to analyze the painted dataset in order to identify homogenous clusters. Visualization of the posterior distribution of clusters was then performed using the tree-building algorithm of fineSTRUCTURE.

Spatial analysis of Principal Components (sPCA) was performed with the R⁶⁰ package ADEGENET⁶¹, using the Delaunay triangulation as connection network⁶². Monte Carlo tests were used to check the statistical significance of spatial structures (global and/or local spatial structure) for 10,000 iterations. The sPCA results were visualized by plotting the samples according to their geographic coordinates and coloring them according to their respective scores along the third first sPCA components.

**Data Availability**

Genotypes obtained with the Illumina Ovine SNP50 BeadChip (genotype files and markers files by chromosomes) for Sidaoun and Hamra breeds in a Plink format are available from the Dryad Digital Repository, https://doi.org/10.5061/dryad.24p1k82.

**References**

1. FAO. Characterization and value addition to local breeds and their products in the Near East and North Africa – Regional Workshop, Rabat, Morocco, 19–21 November 2012. Animal Production and Health Report No. 3. Rome (2014).

2. Bruford, M. W. et al. Prospects and challenges for the conservation of farm animal genomic resources, 2015–2025. *Front Genet* **6**, 314 (2015).

3. FAO. The Second Report on the State of the World's Animal Genetic Resources for Food and Agriculture, Ed. by Scherf, B. D & Pilling, D. FAO Commission on Genetic Resources for Food and Agriculture Assessments. Rome (2015).

4. Shrestha, J. N. B. Conserving domestic animal diversity among composite populations. *Small Ruminant Research* **56**, 3–20 (2005).

5. Scherf, B., Rischkowsky, B., Pilling, D. & Hoffmann, I. The state of the world's animal genetic resources. 8th Wrld Congr. Gen. Appl. to Livest. Prod., Belo Horizonte, M. G., Brazil, p. 13–18 (2006).

6. Hassan, F. A., Blench, R. M. & MacDonald, K. C. Climate and cattle in North Africa: a first approximation. The origins and development of African livestock: archaeology, genetics, linguistics and ethnography, London UCL Press, p. 61–86 (2000).

7. Taberlet, P. et al. Are cattle, sheep, and goats endangered species? *Mol. Ecol* **17**, 275–284 (2008).

8. MATET (Ministère de l'Aménagement du Territoire, de l'Environnement et du Tourisme). Quatrième rapport national sur la mise en œuvre de la convention sur la diversité biologique au niveau national. Alger, Algérie (2009).

9. Madani, T., Yakhlef, H. & Abbache, N. Les races bovines, ovines, caprines et camelines. In: Abdelguerfi, A. & S.A. Ramdane. Evaluation des besoins en matière de renforcement des capacités nécessaires à la conservation et l'utilisation durable de la biodiversité importante pour l'agriculture en Algérie, Alger (2003).

10. Djaout, A. et al. Etat de la biodiversité des «races» ovines Algériennes. *Genetics and Biodiversity. Journal* **1**, 1–16 (2017).

11. Gaouar, S. B. S. et al. Admixture and Local Breed Marginalization Threaten Algerian Sheep Diversity. *Plos One* **10**, e0122667 (2015).

12. Gaouar, S. B. S. et al. Genome-wide analysis highlights genetic dilution in Algerian sheep. *Heredity (Edinb)* **118**, 293–301 (2016).

13. Mason, I. L. A Dictionary Of Livestock Breeds (1969).
14. Couput, M. Espèce Ovine, Laine et Industrie Lainière. Ed. Pour L'exposition Universelle De 1900, Section Algérie de M. Couput, directeur du Service des bergeries (1900).

15. Neuditschko, M., Khatkar, M. S. & Raadsma, H. W. NetView: a high-definition network-visualization approach to detect fine-scale population structures from genome-wide patterns of variation. *PLoS ONE* 7, e48375 (2012).

16. Ai, H., Huang, L. & Ren, J. Genetic Diversity, Linkage Disequilibrium and Selection Signatures in Chinese and Western Pigs Revealed by Genome-Wide SNP Markers. *PLoS ONE* 8, e56001 (2013).

17. Upadhyay, M. R. et al. Genetic origin, admixture and population history of aurochs (*Bosprimigenius*) and primitive European cattle. *Heredity* **118**, 169–176 (2017).

18. Purfield, D. C., Berry, D. P., McParland, S. & Bradley, D. G. Runs of homozygosity and population history in cattle. *BMC Genetics* **13**, 70 (2012).

19. Manunza, A. et al. A genome-wide perspective about the diversity and demographic history of seven Spanish goat breeds. *Genet Sel Evol* **48** (2016).

20. Van Dorp, L. et al. Evidence for a Common Origin of Blacksmiths and Cultivators in the Ethiopian Ari within the Last 4500 Years: Lessons for Clustering-Based Inference. *PLoS Genet.* **11**, e1005397 (2015).

21. Bosse, M. et al. Regions of Homozygosity in the Porcine Genome: Consequence of Demography and the Recombination Landscape. *PLoS Genet* **8** (2012).

22. Ben Gara, A. Définition des objectifs de la sélection des ovins de race Barbarine en Tunisie. In: Gabiña D. (ed.) Analysis and definition of the objectives in genetic improvement programs in sheep and goats. An economic approach to increase their profitability. Zaragoza: CIHEAM, p. 111–116 (2000).

23. Chekkal, F. et al. Guide de caractérisation phénomotypique des races ovines de l'Algérie. Ed. CRSTRA (2015).

24. Toutain, G. De la Recherche au Développement - éléments d'Agronomie Saharienne, Ed. INRA, 279 p (1979).

25. Bechchari, A., El Aich, A., Mahyou, H., Baghdad, M. & Bendaou, M. Analyse de l'évolution du système pastoral du Maroc oriental. *Revue d'Elevage et de Médecine Vétérinaire des pays Propicaux* **67**, 4 (2015).

26. Brisebarre, A. Races ovines, systèmes d'élevage et représentations des éleveurs, in Bonte P., Elloumi M., Guillaume H. & Mahdi M. (dir.), Développement rural, environnement et enjeux territoriaux. Regards croisés Oriental marocain et Sud-Est tunisien. Tunis, Cérès Ed., p. 63–78 (2009).

27. Sagne, J. L'Algérie pastorale, ses origines, sa formation, son passé, son présent, son avenir. Préface de Pierre Jore d'Arces. Ed.: Alger, Imprimerie Fontana (1950).

28. Mathieu, M. & Trabut, L. Les hauts plateaux Oranais, rapport de mission. *Alger* (1891).

29. Lebaudy, G. De la steppe aux alpages. La transhumance des moutons africains. In Berbères de rives en Rèves, Abbaye de Daoulas, Finistère, p. 118–128 (2008).

30. MARA. Plan Moutonnier. Rabat, Maroc, Direction de l'Élevage (1980).

31. Cambon, J. Le pays du mouton: des conditions d'existence des troupeaux sur les hauts-plateaux et dans le sud de l'Algérie. Ouvrage publié par ordre de J. Cambon, Gouveneur Général de l'Algérie. Alger, typographie Giralt (1893).

32. Viger, A. Etude sur la question ovine en Algérie. Imprimerie typographique et lithographique G. Mont-Louis, Clermont-Ferrand (1892).

33. Ouchene-Khelifi, N.-A. et al. Genetic homogeneity of North-African goats. *PLoS ONE* **13**, e0202196 (2018).

34. Colli, L. et al. Genome-wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post-domestication migration routes. *Genet Sel Evol* **50**, 58 (2018).

35. Kijas, J. W. et al. Genome-Wide Analysis of the World's Sheep Breeds Reveals High Levels of Historic Mixture and Strong Recent Selection. *PLOS Biology* **10**, e1001258 (2012).

36. Sponenberg, D. P. & Bixby, D. E. Managing Breeds for a Secure Future: Strategies for Breeders and Breed Associations. ALBC (2007).

37. Chellig, R. Les races ovines Algériennes. Office des Publications Universitaires. Alger, p. 1–80 (1992).

38. Boujenane, I., Aït Bihi, N. & Bou Kallouche, A. Amélioration génétique des ovins au Maroc. In: Gabiña D. (Ed.). Strategies for sheep and goat breeding. Zaragoza: CIHEAM. *Cahiers Options Méditerranéennes* **11**, 109–119 (1995).

39. Sanson, A. Les moutons. Histoire naturelle et zootechnie. Quatrième édition. Ouvrage orné de 56 gravures. Paris, Librairie agricole de la maison rustique (1885).

40. Velu, H. Deux ans d'expérimentation à la ferme de Aïn Djemaa. Les Journées Marocaines du Mouton. Casablanca, 13–15 octobre (1934).

41. Couvreur, G. La vie pastorale dans le Haut Atlas Central. *Revue de Géographie du Maroc* **13**, 347 (1968).

42. Nazareno, A. G., Bemmels, J. B., Dick, C. W. & Lohmann, L. G. Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species. *Molecular Ecology Resources* **17**, 1136–1147 (2017).

43. Willing, E.-M., Dreyer, C. & van Oosterhout, C. Estimates of genetic differentiation measured by FST do not necessarily require large sample sizes when using many SNP markers. *PLoS ONE* 7, e42649 (2012).

44. Miller, S. A., Dykes, D. D. & Polesky, H. F. A simple salting out procedure for extracting DNA from human nucleated cells. *Nucleic Acids Res* **16**, 1215 (1988).

45. Danecek, P. et al. The variant call format and VCFtools. *Bioinformatics* **27**, 2156–2158 (2011).

46. Benjelloun, B. Diversité des génomes et adaptation locale des petits ruminants d'un pays méditerranéen: le Maroc. Thèse de Doctorat. Biodiversité, Ecologie, Environnement. Université Grenoble Alpes. Grenoble, France, 207p (2015).

47. Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. *Am. J. Hum. Genet.* **81**, 559–575 (2007).

48. Howrigan, D., Simonson, M. & Keller, M. Detecting autozygosity through runs of homozygosity: a comparison of three autozygosity detection algorithms. *BMC Genomics* **12**, 460 (2011).

49. Lencz, T. et al. Runs of homozygosity reveal highly penetrant recessive loci in schizophrenia. *Proc. Natl. Acad. Sci. USA* **104**, 19942–19947 (2007).

50. McQuillan, R. et al. Runs of homozygosity in European populations. *Am J Hum Genet.* **83**, 359–372 (2008).

51. Browning, S. R. & Browning, B. L. Rapid and Accurate Haplotype Phasing and Missing-Data Inference for Whole-Genome Association Studies By Use of Localized Haplotype Clustering. *Am J Hum Genet.* **81**, 1084–1097 (2007).

52. Weir, B. S. & Cockerham, C. C. Estimating F-Statistics for the Analysis of Population-Structure. *Evolution* **38**, 1358–1370 (1984).

53. Lewis, P. O. & Zaykin, D. Genetic data analysis. Computer program for the analysis of allelic data. Version 1.0. Lewis Labs, Univ. of Connecticut, Storrs, C. T., (1999).

54. Steinig, E. J., Neuditschko, M., Khatkar, M. S., Raadsma, H. W. & Zenger, K. R. Netview p: a network visualization tool to unravel complex population structure using genome-wide SNPs. *Mol. Ecol. Resour.* **16**, 216–227 (2016).

55. Huson, D. H. SplitsTree: analyzing and visualizing evolutionary data. *Bioinformatics* **14**, 68–73 (1998).

56. Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. *Genome Research* **19**, 1655–1664 (2009).

57. Kopelman, N. M., Mayzel, J., Jakobsson, M., Rosenberg, N. A. & Mayrose, I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. *Mol. Ecol. Resour.* **15**, 1179–1191 (2015).

58. Lawson, D. J., Hellenthal, G., Myers, S. & Falush, D. Inference of Population Structure using Dense Haplotype Data. *PLoS Genet* **8**, e1002453 (2012).

59. Delaneau, O., Marchini, J. & Zagury, J.-F. A linear complexity phasing method for thousands of genomes. *Nat. Methods* **9**, 179–181 (2012).
## Acknowledgements

The authors acknowledge the local farmers for their involvement in data collection, as well as D. Benlaichi and Y. Baci for their valuable assistance in the field. The assistance of ITELV and pilot farms, as well as discussions with A.M. Brisebarre, G. Lebaudy, P. Vidal, A. Bechchari and M. Ibnelbachyr were very helpful and appreciated. We also thank A. Djaout, A. Laoun, R. El-Bouyahiaoui and A. Boubekeur for sharing their field experiences with us. The calculations presented in this article were performed on the CALI calculator of the University of Limoges, funded by the Limousin region, the European Union, the XLIM institutes, IPAM, GEIST, and the University of Limoges. This work was supported by the European Union 7th framework project NEXTGEN (Grant Agreement no. 244356), as regards Moroccan data.

## Author Contributions

A.D.S. conceived the study, performed the analyses and drafted the manuscript. B.B., M.L. and F.P. conceived the study and contributed to the writing of the manuscript. I.B., A.O. and B.B. contributed to the analyses and performed sampling. S.B.S.G., A.R.B., S.H. and H.O.O. performed sampling in Algeria. J.A.L. contributed to the analyses and the text. V.B. and A.H. contributed to the conception of the study. E.C. contributed to the writing of the manuscript. D.T.-M. carried out all the analysis work from the historical archives. All authors reviewed and approved the final manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-44137-y.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019

<footnote n="60">60. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/ (2013).</footnote>

<footnote n="61">61. Jombart, T. adegenet: a R package for the multivariate analysis of genetic markers. *Bioinformatics* **24**, 1403–1405 (2008).</footnote>

<footnote n="62">62. Upton, G. J. & Fingleton, B. Spatial data analysis by example, volume 1: Point pattern and quantitative data. Wiley, Toronto Singapore, Brisbane, New York, Chichester (1985).</footnote>