Sequence variation present within the mitochondrial genome was used to investigate genetic diversity within sheep breeds from Asia and Europe. Comparison of 2027 bp of sequence from 121 animals revealed 44 phylogenetically informative nucleotide positions and a single insertion/deletion. A total of 57 haplotypes were observed which formed two distinct clades. Type A haplotypes were found in breeds from Asia (India, Indonesia, Mongolia, and Tibet), while type B haplotypes were observed at the highest frequency in breeds sourced from Europe (nine breeds from Austria, Åland, Finland, Spain, and northwestern Russia). The distribution of haplotypes indicates sheep appear to have the weakest population structure and the highest rate of intercontinental dispersal of any domestic animal reported to date. Only 2.7% of the sequence variation observed was partitioned between continents, which is lower than both goat (approximately 10%) and cattle (approximately 50%). Diagnostic restriction fragment length polymorphism polymerase chain reaction (RFLP-PCR) tests which distinguish type A and B haplotypes were used to test an additional 223 animals from 17 breeds of European and Asian origin. A mixture of the two lineages was found in every breed except Suffolk and the Indian Garole, indicating introgression has played a major part during breed development and subsequent selection.

The goals of the current work were to analyse the population structure of 11 Spanish ovine breeds and to detect genomic regions that may have been targeted by selection. A total of 141 individuals were genotyped with the Infinium 50 K Ovine SNP BeadChip (Illumina). We combined this dataset with Spanish ovine data previously reported by the International Sheep Genomics Consortium (N = 229). Multidimensional scaling and Admixture analyses revealed that Canaria de Pelo and, to a lesser extent, Roja Mallorquina, Latxa and Churra are clearly differentiated populations, while the remaining seven breeds (Ojalada, Castellana, Gallega, Xisqueta, Ripollesa, Rasa Aragonesa and Segureña) share a similar genetic background. Performance of a genome scan with BayeScan and hapFLK allowed us identifying three genomic regions that are consistently detected with both methods i.e. Oar3 (150–154 Mb), Oar6 (4–49 Mb) and Oar13 (68–74 Mb). Neighbor-joining trees based on polymorphisms mapping to these three selective sweeps did not show a clustering of breeds according to their predominant productive specialization (except the local tree based on Oar13 SNPs). Such cryptic signatures of selection have been also found in the bovine genome, posing a considerable challenge to understand the biological consequences of artificial selection.