Genetic basis of ear length in sheep breeds sampled across the region from the Middle East to the Alps.

Ear length in sheep (Ovis aries) shows a wide range of natural variation, from the absence of an outer ear structure (anotia), to small outer ears (microtia), to regular ear length. Up until now, the underlying genetics of this phenotype has been studied in four sheep breeds from China, Jordan and Italy. These studies revealed a broad range of genes significantly associated with ear length, potentially indicating genetic heterogeneity across breeds or geographic regions. In the current study, we performed genome-wide SNP genotyping and haplotype-based mapping, in a population of 340 individuals, to identify loci influencing ear length variation in additional sheep breeds from Slovenia, Croatia, Cyprus and Greece. Additionally, two previously described candidate variants were also genotyped in our mapping population. The mapping model without candidate variant genotypes revealed only one genome-wide significant signal, which was located next to HMX1 on OAR6. This region was previously described as being associated with ear length variation in the Altay and Awassi sheep breeds. The mapping model including the candidate duplication genotype near HMX1 as a fixed effect explained the phenotypic variance on OAR6 and revealed an additional genome-wide significant locus on OAR13 associated with ear length. Our results, combined with published evidence, suggest that a duplication in the evolutionarily conserved region near HMX1 is the major regulator of ear length in sheep breeds descended from a larger region from Central Asia, to the Middle East, Cyprus, Greece and to the Alps. This distribution suggests an ancient origin of the derived allele.

Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats.

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Genomic regions underlying positive selection in local, Alpine cattle breeds.

We used genome-wide SNP data from 18 local cattle breeds from six countries of the Alpine region to characterize population structure and identify genomic regions underlying positive selection. The geographically close breeds Evolèner, Eringer, Valdostana Pezzata Nera, and Valdostana Castana were found to differ from all other Alpine breeds. In addition, three breeds, Simmental, and Original Braunvieh from Switzerland and Pinzgauer from Austria built three separate clusters. Of the 18 breeds studied, the intra-alpine Swiss breed Evolèner had the highest average inbreeding based on runs of homozygosity (FROH ) and the highest average genomic relationship within the breed. In contrast, Slovenian Cika cattle had the lowest average genomic inbreeding and the lowest average genomic relationship within the breed. We found selection signatures on chromosome 6 near known genes such as KIT and LCORL explaining variation in coat color and body size in cattle. The most prominent selection signatures were similar regardless of marker density and the breeds in the data set. In addition, using available high-density SNP data from 14 of the breeds we identified 47 genome regions as ROH islands. The proportion of homozygous animals was higher in all studied animals of local breeds than in Holstein and Brown Swiss cattle, the two most important commercial breeds in the Alpine region. We report ROH islands near genes related to thermoregulation, coat color, production, and stature. The results of this study serve as a basis for the search for causal variants underlying adaptation to the alpine environment and other specific characteristics selected during the evolution of local Alpine cattle breeds.

How Can We Advance Integrative Biology Research in Animal Science in 21st Century? Experience at University of Ljubljana from 2002 to 2022.

In this perspective analysis, we strive to answer the following question: how can we advance integrative biology research in the 21st century with lessons from animal science? At the University of Ljubljana, Biotechnical Faculty, Department of Animal Science, we share here our three lessons learned in the two decades from 2002 to 2022 that we believe could inform integrative biology, systems science, and animal science scholarship in other countries and geographies. Cultivating multiomics knowledge through a conceptual lens of integrative biology is crucial for life sciences research that can stand the test of diverse biological, clinical, and ecological contexts. Moreover, in an era of the current COVID-19 pandemic, animal nutrition and animal science, and the study of their interactions with human health (and vice versa) through integrative biology approaches hold enormous prospects and significance for systems medicine and ecosystem health.

The Consequences of Mitochondrial T10432C Mutation in Cika Cattle: A "Potential" Model for Leber's Hereditary Optic Neuropathy.

While mitogenome mutations leading to pathological manifestations are rare, more than 200 such mutations have been described in humans. In contrast, pathogenic mitogenome mutations are rare in domestic animals and have not been described at all in cattle. In the small local Slovenian cattle breed Cika, we identified (next-generation sequencing) two cows with the T10432C mitogenome mutation in the ND4L gene, which corresponds to the human T10663C mutation known to cause Leber's hereditary optic neuropathy (LHON). Pedigree analysis revealed that the cows in which the mutation was identified belong to two different maternal lineages with 217 individual cows born between 1997 and 2020. The identified mutation and its maternal inheritance were confirmed by Sanger sequencing across multiple generations, whereas no single analysis revealed evidence of heteroplasmy. A closer clinical examination of one cow with the T10432C mutation revealed exophthalmos, whereas histopathological examination revealed retinal ablations, subretinal oedema, and haemorrhage. The results of these analyses confirm the presence of mitochondrial mutation T10432C with homoplasmic maternal inheritance as well as clinical and histopathological signs similar to LHON in humans. Live animals with the mutation could be used as a suitable animal model that can improve our understanding of the pathogenesis of LHON and other mitochondriopathies.

Capturing Genetic Diversity and Selection Signatures of the Endangered Kosovar Balusha Sheep Breed.

There is a growing concern about the loss of animal genetic resources. The aim of this study was to analyze the genetic diversity and potential peculiarity of the endangered Kosovar sheep breed Balusha. For this purpose, a dataset consisting of medium-density SNP chip genotypes (39,879 SNPs) from 45 Balusha sheep was generated and compared with SNP chip genotypes from 29 individuals of a second Kosovar breed, Bardhoka. Publicly available SNP genotypes from 39 individuals of the relatively closely located sheep breeds Istrian Pramenka and Ruda were additionally included in the analyses. Analysis of heterozygosity, allelic richness and effective population size was used to assess the genetic diversity. Inbreeding was evaluated using two different methods (FIS, FROH). The standardized FST (di) and cross-population extended haplotype homozygosity (XPEHH) methods were used to detect signatures of selection. We observed the lowest heterozygosity (HO = 0.351) and effective population size (Ne5 = 25, Ne50 = 228) for the Balusha breed. The mean allelic richness levels (1.780-1.876) across all analyzed breeds were similar and also comparable with those in worldwide breeds. FROH estimates (0.023-0.077) were highest for the Balusha population, although evidence of decreased inbreeding was observed in FIS results for the Balusha breed. Two Gene Ontology (GO) TERMs were strongly enriched for Balusha, and involved genes belonging to the melanogenesis and T cell receptor signaling pathways, respectively. This could result from selection for the special coat color pattern of Balusha (black head) and resistance to certain infectious diseases. The analyzed diversity parameters highlight the urgency to preserve the local Kosovar Balusha sheep as it is clearly distinguished from other sheep of Southeastern Europe, has the lowest diversity level and may harbor valuable genetic variants, e.g., for resistance to infectious diseases.

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.

Estimation of Genetic Parameters of Type Traits in First Parity Cows of the Autochthonous Cika Cattle in Slovenia.

The aim of this study was to estimate genetic parameters of 26 individual and four composite type traits in first parity Cika cows. An analysis of variance was performed with the generalized linear model procedure of the SAS/STAT statistical package, where the fixed effects of year of recording, cow's age at recording and days after calving as a linear regression were included in the model. The variance components for the direct additive genetic effect and the herd effect in all type traits were estimated using the REML method in the VCE-6 software package. The estimated heritabilities ranged from 0.42 to 0.67 for the measured body frame traits, from 0.36 to 0.80 for the scored autochthonous traits, from 0.11 to 0.61 for the scored body frame traits, and from 0.20 to 0.47 for the scored udder traits. The estimated heritabilities for the composite traits called "autochthonous characteristics", "muscularity", "body frame" and "udder" were 0.55, 0.19, 0.19, and 0.26, respectively. The estimated genetic correlations among the measured body frame traits were positive and high, while the majority of them among the scored body frame traits were low to moderate. The estimated proportions of variance explained by the herd effect for the composite traits "autochthonous characteristics," "muscularity," "body frame" and "udder" were 0.09, 0.28, 0.14, and 0.10, respectively. The estimated heritabilities for the type traits of first parity Cika cows were similar to those reported for other breeds where breeding values have been routinely predicted for a long time. All estimated genetic parameters are already used for breeding value prediction in the Cika cattle population.

Post-genotyping optimization of dataset formation could affect genetic diversity parameters: an example of analyses with alpine goat breeds.

BACKGROUND: Local breeds retained unique genetic variability important for adaptive potential especially in light of challenges related to climate change. Our first objective was to perform, for the first time, a genome-wide diversity characterization using Illumina GoatSNP50 BeadChip of autochthonous Dreznica goat breed from Slovenia, and five and one local breeds from neighboring Austria and Italy, respectively. For optimal conservation and breeding programs of endangered local breeds, it is important to detect past admixture events and strive for preservation of purebred representatives of each breed with low or without admixture. In the second objective, we hence investigated the effect of inclusion or exclusion of outliers from datasets on genetic diversity and population structure parameters. RESULTS: Distinct genetic origin of the Dreznica goat was demonstrated as having closest nodes to Austrian and Italian breeds. A phylogenetic study of these breeds with other goat breeds having SNP data available in the DRYAD repository positioned them in the alpine, European and global context. Swiss breeds clustered with cosmopolitan alpine breeds and were closer to French and Spanish breeds. On the other hand, the Dreznica goat, Austrian and Italian breeds were closer to Turkish breeds. Datasets where outliers were excluded affected estimates of genetic diversity parameters within the breed and increased the pairwise genetic distances between most of the breeds. Alpine breeds, including Dreznica, Austrian and Italian goats analyzed here, still exhibit relatively high levels of genetic variability, homogeneous genetic structure and strong geographical partitioning. CONCLUSIONS: Genetic diversity analyses revealed that the Slovenian Dreznica goat has a distinct genetic identity and is closely related to the neighboring Austrian and Italian alpine breeds. These results expand our knowledge on phylogeny of goat breeds from easternmost part of the European Alps. The here employed outlier test and datasets optimization approaches provided an objective and statistically powerful tool for removal of admixed outliers. Importance of this test in selecting the representatives of each breed is warranted to obtain more objective diversity parameters and phylogenetic analysis. Such parameters are often the basis of breeding and management programs and are therefore important for preserving genetic variability and uniqueness of local rare breeds.

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.

TLR2 polymorphism (rs650082970) is associated with somatic cell count in goat milk.

Pathogens invading the mammary gland are recognized through a range of pattern recognition receptors (PRRs), residing on the plasma membrane of mammary epithelial cells. Toll-like receptor 2 (TLR2) signalling is responsible for recognition of Gram-positive bacteria, which are the most common mastitis-causing pathogens in goats. Somatic cell counts (SCC) in milk are routinely determined in goat dairy flocks and serve as an indicator of milk quality, which is highly correlated to intramammary infections. Recently, a single nucleotide polymorphism of the TLR2 was suggested to be associated with SCC in goat milk. To further test the suggested association, we genotyped 61 Slovenian Alpine goats included in the dataset. The effect of the genotype was analysed using the general linear model (GLM) procedure of SAS/STAT software. We found the TLR2 genotypes significantly (p = 0.0007) associated with milk SCC. Animals with the A/G genotype had significantly (p ≤ 0.05) lower SCC value in milk compared to the G/G genotype. Our data suggest that the A allele is the minor one and is associated with lower milk SCC. In the current study, we provide a validated PCR-RFLP based genotyping assay for the TLR2 SNP (rs650082970) and confirm its association with milk SCC. Further studies to confirm the association on a larger number of animals of different breeds and to explain functional consequences of the polymorphism in relation to SCC are encouraged.

Genetic diversity of local cattle.

The Slovak Pinzgau breed faces the bottleneck effect and the loss of diversity due to unequal use of founders and a significant population decline. Further population size reduction can lead to serious problems. Information obtained here and in other studies from high-throughput genotyping of 179 individuals was used to characterise genetic diversity and differentiation of Slovak Pinzgau, Austrian Pinzgau, Cika and Piedmontese cattle by Bayesian clustering algorithm. A gene flow network for the clusters estimated from admixture results was produced. The low estimate of genetic differentiation (FST) in Pinzgau cattle populations indicated that differentiation among these populations is low, particularly owing to a common historical origin and high gene flow. Changes in the log marginal likelihood indicated Austrian Pinzgau as the most similar breed to Slovak Pinzgau. All populations except the Piedmontese one displayed two ways of gene flow among populations, indicating that Piedmontese cattle was involved in producing of the analysed breeds while these breeds were not involved in creation of Piedmontese. Genetic evaluation represents an important tool in breeding and cattle selection. It is more strategically important than ever to preserve as much of the livestock diversity as possible, to ensure a prompt and proper response to the needs of future generations. Information provided by the fine-scale genetic characterization of this study clearly shows that there is a difference in genetic composition of Slovak and Austrian populations, as well as the Cika and Piedmontese cattle. Despite its population size, the Slovak Pinzgau cattle have a potential to serve as a basic gene reserve of this breed, with European and world-wide importance.

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.

Recovery of native genetic background in admixed populations using haplotypes, phenotypes, and pedigree information--using Cika cattle as a case breed.

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.