Genetic Diversity of Montenegrin Local Sheep Breeds Based on Microsatellite Markers.

The Montenegrin sheep population mostly consists of local breeds and their crossbreeds that are very valuable from their genome preservation point of view. The aim of this study was the investigation of the genetic diversity of seven Montenegrin sheep breeds (Jezeropivska-JP, Sora-SOR, Zetska zuja-ZZ, Bardoka-BAR, Sjenicka-SJ, Ljaba-Lj, and Piperska zuja-PIP) using 18 microsatellite sets of markers. The genotyping was done for 291 samples from seven populations using the multiplex amplification of sequences with polymerase chain reaction (PCR). The parameters of genetic diversity were estimated using several software tools. In total, 243 alleles were found, with a range of 6 to 25 by locus. The mean observed heterozygosity (Ho), polymorphism information content (PIC), and Fis values (fixation index) per marker were 0.728, 0.781, and -0.007, respectively. The mean number of alleles per breed varied from 4.889 in ZZ to 10.056 in JP. The highest Ho was estimated for JP (0.763) and the lowest for ZZ (0.640). The genetic structure showed close relations between SOR and JP, and both of them with SJ, while ZZ, LJ, and PIP were more distanced. This study provides useful indicators for the development of further in-depth studies and the creation of appropriate conservation programs.

Large-scale mitogenome sequencing reveals consecutive expansions of domestic taurine cattle and supports sporadic aurochs introgression.

The contribution of domestic cattle in human societies is enormous, making cattle, along with other essential benefits, the economically most important domestic animal in the world today. To expand existing knowledge on cattle domestication and mitogenome diversity, we performed a comprehensive complete mitogenome analysis of the species (802 sequences, 114 breeds). A large sample was collected in South-east Europe, an important agricultural gateway to Europe during Neolithization and a region rich in cattle biodiversity. We found 1725 polymorphic sites (810 singletons, 853 parsimony-informative sites and 57 indels), 701 unique haplotypes, a haplotype diversity of 0.9995 and a nucleotide diversity of 0.0015. In addition to the dominant T3 and several rare haplogroups (Q, T5, T4, T2 and T1), we have identified maternal line in Austrian Murbodner cattle that possess surviving aurochs' mitochondria haplotype P1 that diverged prior to the Neolithization process. This is convincing evidence for rare female-mediated adaptive introgression of wild aurochs into domesticated cattle in Europe. We revalidated the existing haplogroup classification and provided Bayesian phylogenetic inference with a more precise estimated divergence time than previously available. Occasionally, classification based on partial mitogenomes was not reliable; for example, some individuals with haplogroups P and T5 were not recognized based on D-loop information. Bayesian skyline plot estimates (median) show that the earliest population growth began before domestication in cattle with haplogroup T2, followed by Q (~10.0-9.5 kyBP), whereas cattle with T3 (~7.5 kyBP) and T1 (~3.0-2.5 kyBP) expanded later. Overall, our results support the existence of interactions between aurochs and cattle during domestication and dispersal of cattle in the past, contribute to the conservation of maternal cattle diversity and enable functional analyses of the surviving aurochs P1 mitogenome.

On the origin of European sheep as revealed by the diversity of the Balkan breeds and by optimizing population-genetic analysis tools.

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 [Formula: see text] 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. CONCLUSIONS: Our results contribute to a better understanding of the events that have created the present diversity pattern, which is relevant for the management of the genetic resources represented by the European sheep population.

Differentiation of some Pramenka sheep breeds based on morphometric characteristics.

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, Sjenicka, Pivska Pramenka, and Zeta Zuja) 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 Sjenicka breed was one of the largest breeds (body weight 76.4 kg, wither height 72.7 cm, chest circumference 100.3 cm), while Zeta Zuja 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 Sjenicka 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 Zuja 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 Sjenicka and Istrian with Bela Krajina Pramenka in two clusters, while Bardoka and Zeta Zuja were clustered separately.

Conservation of a domestic metapopulation structured into related and partly admixed strains.

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

Extensive polymorphism of the major histocompatibility complex DRA gene in Balkan donkeys: perspectives on selection and genealogy.

The major histocompatibility complex (MHC) contains genes important for immune response in mammals, and these genes exhibit high polymorphism and diversity. The DRA gene, a member of the MHC class II family, is highly conserved across a large number of mammalian species, but it displays exceptionally rich sequence variations in Equidae members. We analyzed allelic polymorphism of the DRA locus in 248 donkeys sampled across the Balkan Peninsula (Albania, Bulgaria, Croatia, Macedonia, Greece and Montenegro). Five known alleles and two new alleles were identified. The new allele Eqas-DRA*0601 was found to carry a synonymous mutation, and new allele Eqas-DRA*0701, a non-synonymous mutation. We further analyzed the historical selection and allele genealogy at the DRA locus in equids. Signals of positive selection obtained by various tests were ambiguous. A conservative conclusion is that DRA polymorphism occurred relatively recently and that positive selection has been acting on the DRA locus for a relatively brief period.

Conservation priorities of genetic diversity in domesticated metapopulations: a study in taurine cattle breeds.

We estimated neutral diversity of 21 European cattle breeds with 105 microsatellites. Nine of them resembled unselected Balkan Busa strains with diffuse breeding barriers and the 12 others were strongly differentiated, isolated breeds. Because of the impact of neutral genetic diversity on long-term population adaptive capacity, we discuss the long-term outcome of different conservation priorities in a subdivided metapopulation of the investigated cattle breeds. The optimal contribution to a pool of total genetic diversity allocated more than 95% of long-term relevant neutral diversity to virtually unselected strains of the Balkan Busa, while the maximization of total variance preferred inbred breeds. Current artificial selection methods, such as genomic selection sped up and a recovery of underestimated traits becomes quickly impossible. We emphasize that currently neutral and even deleterious alleles might be required for future genotypes in sustainable and efficient livestock breeding and production systems of a 21st century. We provide cumulative evidences that long-term survival relies on genetic complexity and complexity relies on allelic diversity. Our results suggest that virtually unselected, nonuniform strains harbor a crucial proportion of neutral diversity and should be conserved with high global priority. As one example, we suggest a cooperative maintenance of the nondifferentiated, highly fragmented, and fast vanishing metapopulation of Balkan Busa.