Background: In the Neolithic, domestic sheep migrated into Europe and subsequently spread in westerly and northwesterly directions. Reconstruction of these migrations and subsequent genetic events requires a more detailed characterization of the current phylogeographic differentiation. Results: We collected 50 K single nucleotide polymorphism (SNP) profiles of Balkan sheep that are currently found near the major Neolithic point of entry into Europe, and combined these data with published genotypes from southwest-Asian, Mediterranean, central-European and north-European sheep and from Asian and European mouflons. We detected clines, ancestral components and admixture by using variants of common analysis tools: geography-informative supervised principal component analysis (PCA), breed-specific admixture analysis, across-breed f 4 profiles and phylogenetic analysis of regional pools of breeds. The regional Balkan sheep populations exhibit considerable genetic overlap, but are clearly distinct from the breeds in surrounding regions. The Asian mouflon did not influence the differentiation of the European domestic sheep and is only distantly related to present-day sheep, including those from Iran where the mouflons were sampled. We demonstrate the occurrence, from southeast to northwest Europe, of a continuously increasing ancestral component of up to 20% contributed by the European mouflon, which is assumed to descend from the original Neolithic domesticates. The overall patterns indicate that the Balkan region and Italy served as post-domestication migration hubs, from which wool sheep reached Spain and north Italy with subsequent migrations northwards. The documented dispersal of Tarentine wool sheep during the Roman period may have been part of this process. Our results also reproduce the documented 18th century admixture of Spanish Merino sheep into several central-European breeds.

No abstract

For the purpose of the morphometric characterization and differentiation of local sheep breeds that belong to the group of breeds called Pramenka or Zackel, two Slovenian (Bela Krajina and Istrian Pramenka) and four Montenegrin sheep breeds (Bardoka, Sjenička, Pivska Pramenka, and Zeta Žuja) were studied. The variation of morphometric measures and nine morphometric indices were analysed. Principal component analysis (PCA) was applied in order to provide an easier description of body size and shape. Regarding body size, the Sjenička breed was one of the largest breeds (body weight 76.4 kg, wither height 72.7 cm, chest circumference 100.3 cm), while Zeta Žuja had the smallest body size (37.1 kg, 64.8, and 81.9 cm). On the other hand, Slovenian Istrian Pramenka had the largest body length, chest depth, chest width, and rump width among all included breeds (79.4, 33.6, 22.7, and 21.2 cm). Bela Krajina, Istrian Pramenka, and the Sjenička breed, according to the index of body frame (IBF) value (107–114), have a rectangular body frame, Bardoka and Pivska Pramenka have a square body frame (99.3–100), and Zeta Žuja has a short body frame (91.8). The PCA of all morphometric parameters extracted three components accounting for 96.6 % of the cumulative variance. An unweighted pair–group method with arithmetic mean (UPGMA) cluster analysis by Euclidian distance shows diversity among the studied breeds, through it grouped Pivska Pramenka with Sjenička and Istrian with Bela Krajina Pramenka in two clusters, while Bardoka and Zeta Žuja were clustered separately.

Preservation of genetic diversity is one of the most pressing challenges in the planetary boundaries concept. Within this context, we focused on genetic diversity in a native, unselected and highly admixed domesticated metapopulation. A set of 1,828 individuals from 60 different cattle breeds was analysed using a medium density SNP chip. Among these breeds, 14 Buša strains formed a metapopulation represented by 350 individuals, while the remaining 46 breeds represented the global cattle population. Genetic analyses showed that the scarcely selected and less differentiated Buša metapopulation contributed a substantial proportion (52.6%) of the neutral allelic diversity to this global taurine population. Consequently, there is an urgent need for synchronized maintenance of this highly fragmented domestic metapopulation, which is distributed over several countries without sophisticated infrastructure and highly endangered by continuous replacement crossing as part of the global genetic homogenization process. This study collected and evaluated samples, data and genomewide information and developed genome-assisted cross-border conservation concepts. To detect and maintain genetic integrity of the metapopulation strains, we designed and applied a composite test that combines six metrics based on additive genetic relationships, a nearest neighbour graph and the distribution of semiprivate alleles. Each metric provides distinct information components about past admixture events and offers an objective and powerful tool for the detection of admixed outliers. The here developed conservation methods and presented experiences could easily be adapted to comparable conservation programmes of domesticated or other metapopulations bred and kept in captivity or under some other sort of human control.

The aim of this study was to obtain unbiased estimates of the diversity parameters, the population history, and the degree of admixture in Cika cattle which represents the local admixed breeds at risk of extinction undergoing challenging conservation programs. Genetic analyses were performed on the genome-wide Single Nucleotide Polymorphism (SNP) Illumina Bovine SNP50 array data of 76 Cika animals and 531 animals from 14 reference populations. To obtain unbiased estimates we used short haplotypes spanning four markers instead of single SNPs to avoid an ascertainment bias of the BovineSNP50 array. Genome-wide haplotypes combined with partial pedigree and type trait classification show the potential to improve identification of purebred animals with a low degree of admixture. Phylogenetic analyses demonstrated unique genetic identity of Cika animals. Genetic distance matrix presented by rooted Neighbour-Net suggested long and broad phylogenetic connection between Cika and Pinzgauer. Unsupervised clustering performed by the admixture analysis and two-dimensional presentation of the genetic distances between individuals also suggest Cika is a distinct breed despite being similar in appearance to Pinzgauer. Animals identified as the most purebred could be used as a nucleus for a recovery of the native genetic background in the current admixed population. The results show that local well-adapted strains, which have never been intensively managed and differentiated into specific breeds, exhibit large haplotype diversity. They suggest a conservation and recovery approach that does not rely exclusively on the search for the original native genetic background but rather on the identification and removal of common introgressed haplotypes would be more powerful. Successful implementation of such an approach should be based on combining phenotype, pedigree, and genome-wide haplotype data of the breed of interest and a spectrum of reference breeds which potentially have had direct or indirect historical contribution to the genetic makeup of the breed of interest.