Genomic analysis of conservation status, population structure, and admixture in local Czech and Slovak dairy goat breeds.

Although dairy goat production, characterized by traditional production on small farms, is an important source of income in the Czech Republic and Slovakia, locally adapted breeds have not been fully consolidated over the last 100 yr due to large fluctuations in population size and inconsistent breeding programs that allowed for different crossbreeding strategies. Our main objective in this study was therefore to assess the conservation status of 4 Czech (Alpine Goat, White Shorthair, Brown Shorthair, and Czech Landrace) and 1 Slovak (Slovak White Shorthair) local goat breeds, to analyze their population structure and admixture, and to estimate their relatedness to several neighboring breeds. Our analyses included 142 goats belonging to 5 local breeds genotyped with the Illumina 50K BeadChip, and 618 previously genotyped animals representing 15 goat breeds from Austria and Switzerland (all analyses based on 46,862 autosomal SNPs and 760 animals). In general, the conservation status of the Czech and Slovak local goat breeds was satisfactory, with the exception of the Brown Shorthair goat, as the analyzed parameters (heterozygosity, haplotype richness, runs of homozygosity-based inbreeding, and effective population size) were mostly above the median of 20 breeds. However, for all 5 Czech and Slovakian breeds, an examination of historical effective population size indicated a substantial decline about 8 to 22 generations ago. In addition, our study revealed that the Czech and Slovakian breeds are not fully consolidated; for instance, White Shorthair and Brown Shorthair were not clearly distinguishable. Considerable admixture, especially in Czech Landrace (effective number of parental clusters = 4.2), and low but numerous migration rates from other Austrian and Swiss breeds were found. These results provide valuable insights for future breeding programs and genetic diversity management of local Czech and Slovak goat breeds.

Genetic Parameters for Limousine Interbeef Genetic Evaluation of Calving Traits.

The aim of this study was to estimate across-country genetic correlations for calving traits (birth weight, calving ease) in the Limousine breed. Correlations were estimated for eight populations (Czech Republic, joint population of Denmark, Finland, and Sweden, France, Great Britain, Ireland, Slovenia, Switzerland, and Estonia). An animal model on raw performance accounting for across-country interactions (AMACI) was used. (Co)variance components were estimated for pairwise combinations of countries. Fixed and random effects were defined by each country according to its national genetic evaluation system. The average across-country genetic correlation for the direct genetic effect was 0.85 for birth weight (0.69-0.96) and 0.75 for calving ease (0.62-0.94). The average correlation for the maternal genetic effect was 0.57 for birth weight and 0.61 for calving ease. After the estimation of genetic parameters, the weighted bending procedure was used to compute the full Interbeef genetic correlation matrix. After bending, direct genetic correlations ranged from 0.62 to 0.84 (with an average of 0.73) for birth weight and from 0.58 to 0.82 (with an average of 0.68) for calving ease.

Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds.

BACKGROUND: The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type. RESULTS: We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis. CONCLUSIONS: Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.

Genetic Parameters for a Weighted Analysis of Survivability in Dairy Cattle.

The genetic parameters for the survival of Holstein cows, analysed in nine consecutive time periods during the first three calving intervals, were estimated. The earlier the animals are culled, the more they are informationally underestimated. This undervaluing can be remedied by using a weighted analysis that balances the amount of information. If the method of estimating breeding values changes, the genetic parameters will also change. The Holstein cattle dataset from 2005 to 2017 used in this study included 1,813,636 survival records from 298,290 cows. The pedigree with three generations of ancestors included 660,476 individuals. Linear repeatability models estimated genetic parameters for overall and functional survivability. Due to weights, heritability increased from 0.013 to 0.057. Repeatability with weights was 0.505. The standard deviations of breeding values were 1.75 and 2.18 without weights and 6.04 and 6.20 with weights. Including weights in the calculation increased the additive variance proportion and the breeding values' reliabilities. We conclude that the main contribution of the weighted method we have presented is to compensate for the lack of records in culled individuals with a positive impact on the reliability of the breeding value.

Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms.

Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of FST were higher in Wallachian sheep (FST = 0.14) than in Sumava sheep (FST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the Ne estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (FHOM) among individuals, from -0.04 to 0.16 in Sumava and from -0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.

Possible effects of the MC4R Asp298Asn polymorphism on pig production traits under ad libitum versus restricted feeding.

The missense mutation Asp298Asn in the melanocortin 4 receptor (MC4R) is associated with daily gain or fatness in pigs (Sus scrofa domesticus). However, to the best of our knowledge, no study has directly compared the effects of the polymorphism between different feeding levels, even though diet plays a vital role in the swine industry. To explore possible differences, data from 439 mostly commercial hybrids fattened ad libitum and 119 commercial hybrids fattened with restricted feed ration were collected. The recorded traits were average daily gain (ADG), feed conversion ratio (FCR), carcass weight (CW), dressing percentage (DP), lean meat content (LM), backfat thickness (BFT), lean cuts weight, and meat quality parameters such as pH, temperature, drip loss, and CIELAB colour space. The general linear model revealed that the overall effect of MC4R was not statistically significant, but significant differences (p < 0.05) were found in ADG, FCR, CW, DP, LM, and BFT. In the ad libitum category, the AA genotype (298Asn/298Asn) tended to be the most favourable for growth-related traits, with the lowest LM, which is consistent with previous findings. In the restricted category, on the other hand, GA heterozygotes (298Asp/298Asn) achieved the best performance in terms of growth, whereas AA homozygotes showed the worst performance. Therefore, these results raise the possibility of an interaction between MC4R and the feeding level.

Genome-Wide Association Study for Body Conformation Traits and Fitness in Czech Holsteins.

The aim of this study was a genome-wide association study (GWAS) on conformation traits using 25,486 genotyped Czech Holsteins, with 35,227 common SNPs for each genotype. Linear trait records were collected between 1995 and 2020. The Interbull information from Multiple Across Country Evaluation (MACE) was included for bulls that mostly had daughter records in a foreign country. When using the Bonferroni correction, the number of SNPs that were either significant or approached the significance threshold was low-dairy capacity composite on BTA4, feet and legs composite BTA21, total score BTA10, stature BTA24, body depth BTA6, angularity BTA20, fore udder attachment BTA10. Without the Bonferroni correction, the total number of significant or near of significance SNPs was 32. The SNPs were localized on BTA1,2,4,5,6,7,8,18,22,25,26,28 for dairy capacity composite, BTA15,21 for feet and legs composite, BTA10 for total score, BTA24 stature, BTA6,23 body depth, BTA20 angularity, BTA2 rump angle, BTA9,10 rear legs rear view, BTA2,19 rear legs side view, BTA10 fore udder attachment, BTA2 udder depth, BTA10 rear udder height, BTA12 central alignment, BTA24 rear teat placement, BTA8,29 rear udder width. The results provide biological information for the improvement of body conformation and fitness in the Holstein population.

Genetic Parameters of Performance and Conformation Traits of 3-Year-Old Warmblood Sport Horses in the Czech Republic.

The aim of this study was to estimate the genetic parameters of a one-day performance test together with the linear type traits of 3-year-old warmblood horses. The study of genetic parameters was based on 5958 tested horses in the period 1998-2021. A total of 22 traits of linear description, three quantitatively measured traits, and one summary mark from the performance test were tested. The model equation included the fixed effect of gender and combination effects of classifier-year of evaluation-place. A single-trait animal model was used for the estimation of heritability and genetic variance, while the two-trait animal model was applied for the estimation of variance and covariance between all traits. The heritability of the overall score of the performance test was 0.25. The range for heritability was between 0.04 and 0.33 for the linear type traits and between 0.46 and 0.57 for the quantitatively measured traits. Genetic correlations were between -0.47 and 0.92. The estimated genetic parameters suggest that the results from the performance test can be incorporated into genetic evaluation in the Czech Republic.

Genealogical analysis of European bison population revealed a growing up population despite very low genetic diversity.

In 1919, the European bison population became extinct in the wild. The rescue of the lowland subspecies and the whole species was achieved mainly thanks to individuals from the Bialowieza Forest (Polish-Belarusian border). There are currently two breeding lines-the lowland (purebred B. b. Bonasus) founded by 7 individuals and the lowland-Caucasian (hybrids of B. b. Bonasus and B. b. caucasicus) founded by 12 individuals. This genealogical study was conducted on 15,071 individuals recorded in the pedigree book between 1881 and 2020. Its objective was to determine the level of genetic variability and inbreeding almost 100 years after the rescue measures were initiated. The completeness of the pedigree of the reference population was 77% in the fifth generation backwards. A maximum of 23 generations can be traced back in the pedigree. The average inbreeding coefficient and the mean average relatedness of the reference population were very high, about 17% and 16% respectively. No significant amount of new inbreeding was discovered. The reference population has lost 9.11% of the total genetic diversity compared to the population of founders. A male of the Caucasian subspecies Kaukasus was discovered among the ancestors of the lowland lineage reference population. The effective population size calculated based on the increase in inbreeding was 23.93 individuals, based on complete generations equivalent it was 16.1 individuals. Wright's F-statistics showed very small differences in genotypic frequencies between individuals within the two lineages in the reference population (FIS = 0.10), between individuals and the total population (FIT = 0.04) and low differentiation between lineages (FST = 0.06). The population of the European bison from the Bialowieza Forest is generally very uniform but still shows good fitness.

Identification of Selection Signals on the X-Chromosome in East Adriatic Sheep: A New Complementary Approach.

Sheep are one of the most important livestock species in Croatia, found mainly in the Mediterranean coastal and mountainous regions along the East Adriatic coast, well adapted to the environment and mostly kept extensively. Our main objective was therefore to map the positive selection of the X-chromosome (18,983 SNPs that passed quality control), since nothing is known about the adaptation genes on this chromosome for any of the breeds from the Balkan cluster. Analyses were performed on a sample of eight native Croatian breeds (101 females and 100 males) representing the East Adriatic metapopulation and on 10 mouflons (five females and males), all sampled in Croatia. Three classical within-population approaches (extreme Runs of Homozygosity islands, integrated Haplotype Score, and number of Segregating Sites by Length) were applied along with our new approach called Haplotype Richness Drop (HRiD), which uses only the information contained in male haplotypes. We have also shown that phylogenetic analyses, such as the Median-joining network, can provide additional information when performed with the selection signals identified by HRiD. Our new approach identifies positive selection signals by searching for genomic regions that exhibit a sudden decline in haplotype richness. In total, we identified 14 positive selection signals, 11 using the classical approach and three using the HRiD approach, all together containing 34 annotated genes. The most reliable selection signal was mapped by all four approaches in the same region, overlapping between 13.17 and 13.60 Mb, and assigned to the CA5B, ZRSR2, AP1S2, and GRPR genes. High repeatability (86%) of results was observed, as 12 identified selection signals were also confirmed in other studies with sheep. HRiD offers an interesting possibility to be used complementary to other approaches or when only males are genotyped, which is often the case in genomic breeding value estimations. These results highlight the importance of the X-chromosome in the adaptive architecture of domestic ruminants, while our novel HRiD approach opens new possibilities for research.