Unraveling the genetic diversity of Belgian Milk Sheep using medium-density SNP genotypes

Genomic characterisation and diversity assessment of eight endangered Belgian sheep breeds

Assessing the genetic diversity of local breeds is essential for conserving these unique breeds, which may possess unique traits. This study provides the genomic characterisation of eight indigenous sheep breeds in Belgium based on pedigree and single nucleotide polymorphism (SNP) analysis. A total of 687 sheep were genotyped and were subjected to a rigorous quality control, resulting in a set of 45 978 autosomal SNPs. Pedigree analysis showed breed-average inbreeding estimates between 3.3% and 11.3%. The genomic analysis included an assessment of runs of homozygosity (ROH) to examine the genomic inbreeding coefficient, with breed-average inbreeding coefficients estimated between 4.1% and 8.5%. Runs of homozygosity islands were identified in six of the eight breeds studied, with some exhibiting an incidence of up to 58%. Interestingly, several ROH islands overlapped with other breeds included in this study, as well as with international sheep breeds. Pedigree-based effective population sizes were estimated below 100 for all breeds, whereas genomic-based effective population sizes were below 24, indicating that these eight local sheep breeds are endangered. Principal component analysis, admixture analyses, and Fst computations were used to study the population structure and genetic differences. A neighbour-joining tree using 95 international sheep breeds positioned the eight local breeds in the group of milksheep, Texel sheep and the Scandinavian breeds. Additionally, the investigation of paternal oY1 genotypes revealed diverse lineage origins within the Belgian sheep population. This study refines and deepens our knowledge about the local sheep breeds in Belgium, thereby improving their management and conservation. Moreover, as these breeds are linked to other international breeds, these insights are significant for the global scientific community.

Association between single-nucleotide polymorphism in PIK3CD gene and litter size in Small Tail Han sheep

The p110δ isoform of the PI3K catalytic subunit (encoded by the PIK3CD gene) is a key component of the PI3K pathway for follicle growth in mammalian ovarian granulosa cells. Nevertheless, little is known about the association of its polymorphisms with ovine litter size. In this study, the distribution of different genotypes of two SNPs in the PIK3CD gene was calculated in more than forty sheep breeds, and the associations between SNPs and litter size in Small Tail Han (STH) sheep were also analyzed. Besides, the mRNA expression of the PIK3CD gene was also detected in some reproduction-related tissues. The results showed that the "A" allele frequency was higher in rs412889931 (g.41926327G > A) in a typical polytocous sheep breed (p < 0.01). The association's analysis showed rs412889931 was correlated with ovine fecundity as assessed by three parity litter sizes (p < 0.05). Finally, we found the expression of PIK3CD in the ovary had significant differences in different fecundity sheep breeds, indicating that SNP may regulate the litter size by influencing the PIK3CD gene expression. The present results demonstrated that rs412889931 could be used in the marker-assisted selection of the litter size in sheep breeding.

Genetic Diversity and Trends of Ancestral and New Inbreeding in German Sheep Breeds by Pedigree Data

In Germany, many autochthonous sheep breeds have developed, adapted to mountain, heath, moorland, or other marginal sites, but breeds imported from other countries have also contributed to the domestic breeds, particularly improving wool and meat quality. Selective breeding and the intense use of rams may risk losing genetic diversity and increasing rates of inbreeding. On the other hand, breeds with a low number of founder animals and only regional popularity may not leave their endangered status, as the number of breeders interested in the breed is limited. The objective of the present study was to determine demographic measures of genetic diversity and recent as well as ancestral trends of inbreeding in all autochthonous German sheep breeds and sheep of all breeding directions, including wool, meat, and milk. We used pedigree data from 1,435,562 sheep of 35 different breeds and a reference population of 981,093 sheep, born from 2010 to 2020. The mean number of equivalent generations, founders, effective founders, effective ancestors, and effective founder genomes were 5.77, 1669, 123.2, 63.5, and 33.0, respectively. Genetic drift accounted for 69% of the loss of genetic diversity, while loss due to unequal founder contributions was 31%. The mean inbreeding coefficient, individual rate of inbreeding (∆F_i), and realized effective population size across breeds were 0.031, 0.0074, and 91.4, respectively, with a significantly decreasing trend in ∆F_i in 11/35 breeds. New inbreeding, according to Kalinowski, contributed to 71.8% of individual inbreeding, but ancestral inbreeding coefficients showed an increasing trend in all breeds. In conclusion, in our study, all but one of the mountain-stone sheep breeds and the country sheep breed Wald were the most vulnerable populations, with N_e < 50. The next most endangered breeds are exotic, country, and heath breeds, with average N_e of 66, 83, and 89, respectively. The wool, meat, and milk breeds showed the highest genetic diversity, with average N_e of 158, 120, and 111, respectively. The results of our study should help strengthen conservation program efforts for the most endangered sheep breeds and maintain a high genetic diversity in all sheep breeds.

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 F_ST were higher in Wallachian sheep (F_ST = 0.14) than in Sumava sheep (F_ST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the N_e 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 (F_HOM) 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.

Analysis of Family Structure and Paternity Test of Tan Sheep in Yanchi Area, China

Tan sheep is a special breed of locally protected sheep in China, one of the best quality meat sheep in the world. Due to the unclear pedigree of the rams on the Ningxia Tan sheep breeding farm, we investigated 74 rams in the field and explored a new method for family division. Genomic DNA was extracted from the blood of breeding rams. Using Plink software, GCTA tools and R language, we analyzed the genetic structure, kinship, and inbreeding coefficient of the breeding sheep, which revealed the genetic relationship between the individuals. The results showed that there was no obvious clustering phenomenon in the PCA, and the genetic background of the samples was similar. The G matrix and IBS distance matrix indicated that most individuals were far away from each other. Paternity testing identified 24 pairs of unknown parent-child pairs, and all the Tan sheep could be divided into 12 families, which provided a reference for sheep breeding. The average inbreeding coefficient based on the ROH of this population was 0.049, so there was a low degree of inbreeding and the rams in the field were able to maintain high genetic diversity. Overall, we explored a more accurate method through paternity and kinship analysis; it provides a scientific basis for pedigree construction, which has an important application value for Tan sheep breeding.

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 (F_IS, F_ROH). The standardized F_ST (d_i) and cross-population extended haplotype homozygosity (XPEHH) methods were used to detect signatures of selection. We observed the lowest heterozygosity (H_O = 0.351) and effective population size (Ne₅ = 25, Ne₅₀ = 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. F_ROH estimates (0.023–0.077) were highest for the Balusha population, although evidence of decreased inbreeding was observed in F_IS 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.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.

Genomics Confirm an Alarming Status of the Genetic Diversity of Belgian Red and Belgian White Red Cattle

Simple Summary

Genetic diversity of livestock is vitally important as it enables the adaptation of future populations to changing environments. Therefore, preserving a sufficiently large genetic diversity is key. However, for many local cattle populations, little is known about their genetic diversity such as inbreeding level, effective size etc. We studied the genetic diversity of two local Belgian red cattle populations (Belgian Red and Belgian White Red cattle) using state-of-the-art genomic techniques. These tools assessed diversity at the population and individual level, and allowed the positioning of these two breeds in an international context of 52 other (European) cattle breeds. Accordingly, we contribute to the general knowledge of European red cattle, and more specifically we help the breeders, breed organization and the government to manage the genetic diversity of both breeds.

Abstract

Genetic diversity is increasingly important for researchers and society. Small and local populations deserve more attention especially, as they may harbor important characteristics. Moreover, small populations are at greater risk and their genetic management is often more challenging. Likewise, European red cattle populations are threatened, as they are outcompeted by more specialized cattle breeds. In this study, we investigate the genetic diversity of two local Belgian red cattle breeds: Belgian Red and Belgian White Red cattle. A total of 270 animals were genotyped via medium density SNP arrays. Genetic diversity was assessed using runs of homozygosity screening, effective population size estimation and F_st analyses. Genomic inbreeding coefficients based on runs of homozygosity were estimated at 7.0% for Belgian Red and 6.1% for Belgian White Red cattle, and both populations had a low effective population size (68 and 86, respectively). PCA, F_st and admixture analyses revealed the relationship to 52 other international breeds, where they were closest related to some Belgian, French, Scandinavian and Dutch breeds. Moreover, F_st analyses revealed for Belgian Red cattle a signature of selection on BTA6, adjacent to the KIT gene. This study gains important knowledge on the genetic diversity of these two small local red cattle breeds, and will aid in their (genetic) management.

A comprehensive genome-wide scan detects genomic regions related to local adaptation and climate resilience in Mediterranean domestic sheep

Background

The management of farm animal genetic resources and the adaptation of animals to climate change will probably have major effects on the long-term sustainability of the livestock sector. Genomic data harbour useful relevant information that needs to be harnessed for effectively managing genetic resources. In this paper, we report the genome characterization of the highly productive Mediterranean Chios dairy sheep and focus on genetic diversity measures related with local adaptation and selection and the genetic architecture of animal resilience to weather fluctuations as a novel adaptative trait linked to climate change.

Results

We detected runs of homozygosity (ROH) and heterozygosity (ROHet) that revealed multiple highly homozygous and heterozygous hotspots across the Chios sheep genome. A particularly highly homozygous region was identified on chromosome 13 as a candidate of directional genetic selection associated with milk traits, which includes annotated genes that were previously shown to be linked to local adaptation to harsh environmental conditions. Favourable heterozygosity related with a potentially protective role against livestock diseases and enhanced overall fitness was revealed in heterozygous-rich regions on sheep chromosomes 3, 10, 13 and 19. Furthermore, genomic analyses were conducted on sheep resilience phenotypes that display changes in milk production in response to weather variation. Sheep resilience to heat stress was a significantly heritable trait (h² = 0.26) and genetically antagonistic to milk production. Genome-wide association and regional heritability mapping analyses revealed novel genomic markers and regions on chromosome 5 that were significantly associated with sheep resilience to climate change. Subsequently, an annotation analysis detected a set of genes on chromosome 5 that were associated with olfactory receptor complexes that could participate in heat stress mitigation through changes in respiration rate and respiratory evaporation. Other genes were grouped in previously reported biological processes relevant to livestock heat dissipation, including stress and immune response.

Conclusions

Our results may contribute to the optimal management of sheep genetic resources and inform modern selective breeding programmes that aim at mitigating future environmental challenges towards sustainable farming, while better balancing animal adaptation and productivity. Our results are directly relevant to the studied breed and the respective environmental conditions; however, the methodology may be extended to other livestock species of interest.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00682-7.

An Overview of the Use of Genotyping Techniques for Assessing Genetic Diversity in Local Farm Animal Breeds

Simple Summary

The number of local farm animal breeds is declining worldwide. However, these breeds have different degrees of genetic diversity. Measuring genetic diversity is important for the development of conservation strategies and, therefore, various genomic analysis techniques are available. The aim of the present work was to shed light on the use of these techniques in diversity studies of local breeds. In summary, a total of 133 worldwide studies that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that over time, almost all available genomic techniques were used and various diversity parameters were calculated. Therefore, the present results provide a comprehensive overview of the application of these techniques in the field of local breeds. This can provide helpful insights into the advancement of the conservation of breeds with high genetic diversity.

Abstract

Globally, many local farm animal breeds are threatened with extinction. However, these breeds contribute to the high amount of genetic diversity required to combat unforeseen future challenges of livestock production systems. To assess genetic diversity, various genotyping techniques have been developed. Based on the respective genomic information, different parameters, e.g., heterozygosity, allele frequencies and inbreeding coefficient, can be measured in order to reveal genetic diversity between and within breeds. The aim of the present work was to shed light on the use of genotyping techniques in the field of local farm animal breeds. Therefore, a total of 133 studies across the world that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that diversity of cattle was most often investigated with microsatellite use as the main technique. Furthermore, a large variety of diversity parameters that were calculated with different programs were identified. For 15% of the included studies, the used genotypes are publicly available, and, in 6%, phenotypes were recorded. In conclusion, the present results provide a comprehensive overview of the application of genotyping techniques in the field of local breeds. This can provide helpful insights to advance the conservation of breeds.

Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area

Under a climate change perspective, the genetic make-up of local livestock breeds showing adaptive traits should be explored and preserved as a priority. We used genotype data from the ovine 50 k Illumina BeadChip for assessing breed autozygosity based on runs of homozygosity (ROH) and fine-scale genetic structure and for detecting genomic regions under selection in 63 Tunisian sheep samples. The average genomic inbreeding coefficients based on ROH were estimated at 0.017, 0.021, and 0.024 for Barbarine (BAR, n = 26), Noire de Thibar (NDT, n = 23), and Queue fine de l'Ouest (QFO, n = 14) breeds, respectively. The genomic relationships among individuals based on identity by state (IBS) distance matrix highlighted a recent introgression of QFO into the BAR and a genetic differentiation of NDT samples, possibly explained by past introgression of European gene pools. Genome-wide scan for ROH across breeds and within the BAR sample set identified an outstanding signal on chromosome 13 (46.58-49.61 Mbp). These results were confirmed using F_ST index, differentiating fat vs. thin-tailed individuals. Candidate genes under selection pressure (CDS2, PROKR1, and BMP2) were associated to lipid storage and probably preferentially selected in fat-tailed BAR animals. Our findings suggest paying more attention to preserve the genetic integrity and adaptive alleles of local sheep 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 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{4}$$\end{document}f4 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.

How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species

Background

PLINK is probably the most used program for analyzing SNP genotypes and runs of homozygosity (ROH), both in human and in animal populations. The last decade, ROH analyses have become the state-of-the-art method for inbreeding assessment. In PLINK, the --homozyg function is used to perform ROH analyses and relies on several input settings. These settings can have a large impact on the outcome and default values are not always appropriate for medium density SNP array data. Guidelines for a robust and uniform ROH analysis in PLINK using medium density data are lacking, albeit these guidelines are vital for comparing different ROH studies. In this study, 8 populations of different livestock and pet species are used to demonstrate the importance of PLINK input settings. Moreover, the effects of pruning SNPs for low minor allele frequencies and linkage disequilibrium on ROH detection are shown.

Results

We introduce the genome coverage parameter to appropriately estimate F_ROH and to check the validity of ROH analyses. The effect of pruning for linkage disequilibrium and low minor allele frequencies on ROH analyses is highly population dependent and such pruning may result in missed ROH. PLINK’s minimal density requirement is crucial for medium density genotypes and if set too low, genome coverage of the ROH analysis is limited. Finally, we provide recommendations for the maximal gap, scanning window length and threshold settings.

Conclusions

In this study, we present guidelines for an adequate and robust ROH analysis in PLINK on medium density SNP data. Furthermore, we advise to report parameter settings in publications, and to validate them prior to analysis. Moreover, we encourage authors to report genome coverage to reflect the ROH analysis’ validity. Implementing these guidelines will substantially improve the overall quality and uniformity of ROH analyses.