Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds

Whole-genome resequencing reveals the population structure and domestication processes of endemic endangered goose breeds (Anser cygnoides)

In recent years, the dwindling population of these endangered geese has hindered our understanding of their phenotypic variations and the genes associated with important traits. To investigate the population structure and genetic diversity of this breed, the whole-genome data of 90 individuals from a conservation farm were obtained using the Illumina 6000 paired-end platform. The research results indicate that each locally endangered goose variety has formed a monophyletic population. The Baizi (BZ), Lingxian White (LX), and Xupu (XP) geese exhibiting higher genetic diversity than the other goose breeds. Tree-Mix analysis revealed the presence of five gene flows events between goose populations, with Yangjiang (YJ) geese consistently exhibiting significant genetic distance from the other breeds. Under strong pressures from the natural environment and artificial selection, whole-genome selective scanning revealed 394 overlapping genes. Gene Ontology (GO) enrichment analysis of the putative candidate genes (PCGs) revealed significant enrichment of 20 terms (P < 0.05). Similarly, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed significant enrichment of PCGs in 23 terms (P < 0.05). Examination of overlapping genes identified through at least two selection methods revealed a set of genes associated with key traits, including growth and development (CCND1, DES, CCNO, SMC5, and NUBP1), immunity (ABCA2, ABCC8, UHRF2, and ABCA1), and body aging (KAT6B). Our findings provide insights into the genetic basis of endangered geese at the whole-genome level, laying the foundation for future molecular research on genetic variation and phenotypic changes. In summary, our results provide invaluable resources for delineating the uniqueness of endangered goose breeds.

Genome-Wide Runs of Homozygosity Reveal Inbreeding Levels and Trait-Associated Candidate Genes in Diverse Sheep Breeds

BACKGROUND

Quantifying and controlling the inbreeding level in livestock populations is crucial for the long-term sustainability of animal husbandry. However, the extent of inbreeding has not been fully understood in sheep populations on a global scale.

METHODS

Here, we analyzed high-depth genomes of 210 sheep from 20 worldwide breeds to identify the pattern and distribution of genome-wide runs of homozygosity (ROH) and detect candidate selected genes in ROH islands for agronomic and phenotypic traits.

RESULTS

Leveraging whole-genome sequencing data, we found a large number of short ROH (e.g., <1.0 Mb) in all breeds and observed the overall higher values of ROH statistics and inbreeding coefficient in European breeds than in Asian breeds and Dorper sheep. We identified some well-known candidate genes (e.g., CAMK4, HOXA gene family, ALOX12, FGF11, and MTOR) and 40 novel genes (e.g., KLHL1, FGFRL1, WDR62, GDF6, KHDRBS2, and PAX1) that are functionally associated with fecundity, body size, and wool-related traits in sheep. Based on the candidate genes, we revealed different genetic bases for the fecundity traits of European and Asian sheep.

CONCLUSIONS

This study improves the resolution of ROH detection and provides new insights into genomic inbreeding and trait architecture in sheep as well as useful markers for future breeding practice.

Identification of selection signatures and genetic diversity in the sheep

In the study, data obtained from OvineSNP50K SNP chips using the Illumina® iScan platform for Eşme sheep were used. The integrated haplotype score (iHS) and runs of homozygosity (ROH) statistical approaches were used to identify selection signatures. Using the iHS analysis, it was discovered that there are 10 genomic regions and 51 genes on ovine chromosomes 1, 9, 11, and 12 that are under selection. Three genomic regions and 97 genes on ovine chromosomes 6 and 11 were found to be under selection using the ROH analysis. Candidate genes associated with economic and ecological traits were detected using both approaches. Among the genetic diversity parameters considered in this study, the minor allele frequency (MAF), the genetic distance between individuals (D), as well as observed (Ho) and expected heterozygosities (He) values were 0.300, 0.309, 0.388, and 0.390, respectively. The obtained Ho, He and D values indicate a moderate level of genetic diversity. The ratio of polymorphic SNPs (PN) was 0.947, and the average values of FROH and FHOM were 0.030 and 0.029, respectively. Considering the PN value obtained in the study, it is evident that the SNPs in the population exhibit a high level of polymorphism at 94.7%. While the FROH value obtained indicates high genetic diversity among the individuals in the present study, the FHOM value suggests that the population is predominantly composed of heterozygous individuals. As a result, evidence indicating genetic advancements have been made for target traits in breeding programs within the population. Additionally, candidate genes suitable for future molecular marker-supported breeding programs have been identified. In addition, a better understanding of the genetic structure and production potential of the population has been achieved. Findings have shown that Eşme sheep are a breed with high meat production potential and strong adaptation abilities.

Genome-Wide Patterns of Homozygosity and Heterozygosity and Candidate Genes in Greek Insular and Mainland Native Goats

BACKGROUND

Runs of homozygosity (ROHs) and heterozygosity (ROHets) serve for the identification of genomic regions as candidates of selection, local adaptation, and population history.

METHODS

The present study aimed to comprehensively explore the ROH and ROHet patterns and hotspots in Greek native dairy goats, Eghoria and Skopelos, genotyped with the Illumina Goat SNP50 BeadChip. SNP and functional enrichment analyses were conducted to further characterize hotspots and the candidate genes located within these genomic regions. Genetic relationships between and within breeds and inbreeding coefficients were also evaluated.

RESULTS

Clear genetic differentiation and diversified management practices were depicted between the two native populations. The ROH and ROHet average genome coverage for Skopelos (65.35 and 35 Mb) and Eghoria (47.64 and 43 Mb) indicated differences in mainland and insular goats, with Skopelos showing more long ROH fragments, reflecting its geographic isolation and small population size. An ROH hotspot (CHR12: 43.59-44.61 Mb) detected in the Skopelos population has been also reported across European goats and co-localizes with a selection signal detected in the Egyptian Barki goats and sheep adapted to hot-arid conditions. A novel ROH hotspot (CHR18: 60.12-61.81 Mb), shared among the Greek breeds, harbors candidate genes enriched in biosynthesis, metabolism, and immune response. Two well-conserved ROHet islands were detected in Greek goats on chromosomes 1 and 18, with genes participating in development and embryogenesis. The Eghoria population showed the highest number of ROHet islands, potentially reflecting its adaptability to diverse environments.

CONCLUSIONS

These findings offer new insights into the environmental adaptation and artificial selection in Greek goats and could be utilized in future breeding strategies for sustainable goat farming.

Genetic Diversity Analysis and Identification of Candidate Genes for Growth Traits in Chengkou Mountain Chicken

Growth traits constitute critical factors in the breeding program of broiler chickens. The Chengkou mountain chicken A-lineage (CMC-A) represents a breed specifically bred for meat production. To further explore the growth performance of the CMC-A population, this study conducted whole-genome sequencing on 464 CMC-A roosters to systematically evaluate their genetic diversity. Additionally, runs of homozygosity (ROH) islands and genome-wide association studies (GWASs) were employed to identify the loci and functional genes influencing the growth traits in Chengkou mountain chickens. The results revealed a high level of genetic diversity and low levels of inbreeding in Chengkou mountain chickens. Several genes associated with stress resistance, muscle growth, and fat deposition were pinpointed through ROH island identification. Moreover, 52 SNP loci were detected, along with 71 candidate genes. These findings enhance our understanding of the genetic architecture underlying the growth traits in Chengkou mountain chickens and provide a theoretical foundation for subsequent breeding endeavors.

Analysis of genetic evolutionary differences among four Tibetan pig populations in China

Tibetan pigs are a locally bred domestic pig breed originating from the Tibetan Plateau in China. They can be categorized into four distinct groups based on their geographical locations: Sichuan Tibetan pigs, Tibetan pigs from Tibet, Yunnan Tibetan pigs, and Gansu Tibetan pigs. This study aimed to explore population diversity, genetic structure and selection signals among Tibetan pigs in four Chinese national nature reserves. The results show that there is different observed heterozygosity among Tibetan pig populations (0.1957-0.1978). Ratio of runs of homozygosity (Froh) calculation of four Tibetan pig populations by runs of homozygosity (ROH) revealed the presence of inbreeding within the population (0.0336-0.0378). Analysis of the genetic structure demonstrated distinct population stratification among the four Tibetan pig populations, with each showing relatively independent evolutionary directions. Furthermore, Five methods (FST, Piratio, ROD, Tajima's D, XP-CLR) were used to artificially select evolutionary trajectories. The results mainly involved processes such as DNA repair, immune regulation, muscle fat deposition and adaptation to hypoxia. In conclusion, this study enhances our understanding of the genetic characteristics of Tibetan pig populations and provides a theoretical reference for the conservation of resources across different populations of Tibetan pigs.

Genome-Wide Association Studies and Runs of Homozygosity Reveals Genetic Markers Associated with Reproductive Performance in Korean Duroc, Landrace, and Yorkshire Breeds

BACKGROUND

Reproductive performance is critical in the pig industry, and improved sow performance could lead to increased economic benefits. GWAS and ROH analyses based on SNP array data were conducted to identify the breed-specific genetic architecture underlying the variation in NBA and TNB.

METHODS

A total of 7488 breeding pigs with phenotypic data from 1586 Duroc, 2256 Landrace, and 3646 Yorkshire breeds, along with 76,756 SNP markers from Korean grand-grand-parent (GGP) breeding farms, were used.

RESULTS

In the Duroc breeds, SNPs on SSC 9 and 17 were found to be associated with the SIDT2 and TGM2 genes, respectively. In the Landrace breed, PPP1R9A, LMTK2, and GTF2H3 on SSCs 9, 3, and 14, respectively, were associated with both TNB and NBA. With the Yorkshire breed genome, GRID1, DLGAP2, ZZEF1, PARG, RNF17, and NDUFAF5 in SSCs 14, 15, 12, 14, 11, and 17, respectively, were associated with NBA and TNB traits. These genes have distinct functions, ranging from synaptic transmission and cytoskeletal organization to DNA repair and cellular energy production. In the Duroc breed, six genes identified in the ROH islands were associated with various biological pathways, molecular functions, and cellular components. NT5DC1 was associated with metaphyseal chondrodysplasia, CRTAC1 with ion binding, CFAP43 with spermatogenic failure, CASC3 with intracellular mRNA localization, ERC2 with cellular component organization, and FOCAD with Focadhesin. In the Landrace and Yorkshire breeds, PDE6D was associated with GTPase inhibitor activity.

CONCLUSIONS

Through GWAS and ROH analyses, we identified breed-specific SNP markers associated with NBA and TNB in three breed genotypes, providing insights for improving reproductive performance efficiency and contributing to future breeding strategies.

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.

Runs of homozygosity analysis for selection signatures in the Yellow Korean native chicken

OBJECTIVE

Yellow Korean native chicken (KNC-Y) is one of the five pure Korean indigenous chicken breeds that were restored through a government project in 1992. KNC-Y is recognized for its superior egg production performance compared to other KNC lines. In this study, we performed runs of homozygosity (ROH) analysis to discover selection signatures associated with egg production traits in the KNC-Y population.

METHODS

A total of 675 DNA samples from KNC-Y were genotyped to generate single nucleotide polymorphism (SNP) data using custom 60K Affymetrix SNP chips. ROH analysis was performed using PLINK software, with predefined parameters set for the analysis. The threshold of ROH island was defined as the top 1% frequency of SNPs withing the ROH among the population.

RESULTS

In the KNC-Y population, a total of 29,958 runs of homozygosity (ROH) fragments were identified. The average total length of ROH was 120.84 Mb, with each ROH fragment having an average length of 2.71 Mb. The calculated ROH-based inbreeding coefficient (FROH) was 0.13. Furthermore, we revealed the presence of ROH islands on chromosomes 1, 2, 4, 5, 7, 8, and 11. Within the identified regions, a total of 111 genes were annotated, and among them were genes related to economic traits, including PRMT3, ANO5, HDAC4, LSS, PLA2G4A, and PTGS2. Most of the overlapping quantitative trait locus regions with ROH islands were found to be associated with production traits.

CONCLUSION

This study conducted a comprehensive analysis of ROH in the KNC-Y population. Notably, among the findings, the PTGS2 gene is believed to play a crucial role in influencing the laying performance of KNC-Y.

Genetic Diversity and Selection Signal Analysis of Hu Sheep Based on SNP50K BeadChip

This research is designed to examine the genetic diversity and kinship among Hu sheep, as well as to discover genes associated with crucial economic traits. A selection of 50 unrelated adult male Hu sheep underwent genotyping with the SNP50K BeadChip. Seven indicators of genetic diversity were assessed based on high-quality SNP data: effective population size (Ne), polymorphic information content (PIC), polymorphic marker ratio (PN), expected heterozygosity (He), observed heterozygosity (Ho), effective number of alleles, and minor allele frequency (MAF). Plink software was employed to compute the IBS genetic distance matrix and detect runs of homozygosity (ROHs), while the G matrix and principal component analysis were performed using GCTA software. Selective sweep analysis was carried out using ROH, Pi, and Tajima's D methodologies. This study identified a total of 64,734 SNPs, of which 56,522 SNPs remained for downstream analysis after quality control. The population displayed relatively high genetic diversity. The 50 Hu sheep were ultimately grouped into 12 distinct families, with families 6, 8, and 10 having the highest numbers of individuals, each consisting of 6 sheep. Furthermore, a total of 294 ROHs were detected, with the majority having lengths between 1 and 5 Mb, and the inbreeding coefficient FROH was 0.01. In addition, 41, 440, and 994 candidate genes were identified by ROH, Pi, and Tajima's D methods, respectively, with 3 genes overlapping (BMPR1B, KCNIP4, and FAM13A). These results offer valuable insights for future Hu sheep breeding, genetic assessment, and population management.

Genomic and Conventional Inbreeding Coefficient Estimation Using Different Estimator Models in Korean Duroc, Landrace, and Yorkshire Breeds Using 70K Porcine SNP BeadChip

The purpose of this study was to estimate the homozygosity distribution and compute genomic and conventional inbreeding coefficients in three genetically diverse pig breed populations. The genomic and pedigree data of Duroc (1586), Landrace (2256), and Yorkshire (3646) were analyzed. We estimated and compared various genomic and pedigree inbreeding coefficients using different models and approaches. A total of 709,384 ROH segments in Duroc, 816,898 in Landrace, and 1,401,781 in Yorkshire, with average lengths of 53.59 Mb, 56.21 Mb, and 53.46 Mb, respectively, were identified. Relatively, the Yorkshire breed had the shortest ROH segments, whereas the Landrace breed had the longest mean ROH segments. Sus scrofa chromosome 1 (SSC1) had the highest chromosomal coverage by ROH across all breeds. Across breeds, an absolute correlation (1.0) was seen between FROH total and FROH1-2Mb, showing that short ROH were the primary contributors to overall FROH values. The overall association between genomic and conventional inbreeding was weak, with values ranging from 0.058 to 0.140. In contrast, total genomic inbreeding (FROH) and ROH classes showed a strong association, ranging from 0.663 to 1.00, across the genotypes. The results of genomic and conventional inbreeding estimates improve our understanding of the genetic diversity among genotypes.

Genome-Wide Association Studies of Growth Trait Heterosis in Crossbred Meat Rabbits

The application of heterosis can not only effectively improve the disease resistance and meat quality of livestock, but also significantly enhance the reproduction and growth of livestock and poultry. We conducted genome-wide association studies using data from F2 crossbred meat rabbits to screen out candidate genes with significant dominant effects associated with economic trait variation. High-throughput sequencing technology was used to obtain SNPs covering the whole genome to evaluate the homozygosity of the population genome, and analyze the number, length, frequency, and distribution of ROHs in the population. Candidate genes related to economic traits of meat rabbits were searched based on high-frequency ROH regions. After quality control filtering of genotype data, 380 F2 crossbred rabbits were identified with 78,579 SNPs and 42,018 ROHs on the autosomes. The fitting of the Logistic growth curve model showed that 49-day-old rabbits were a growth inflection point. Then, through genome-wide association studies, 10 SNP loci and seven growth trait candidate genes were found to be significantly related to body weight in meat rabbits at 84 days of age. In addition, we revealed the functional roles and locations of 20 candidate genes in the high-frequency ROH region associated with economic traits in meat rabbits. This study identified potential genes associated with growth and development in the high-frequency ROH region of meat rabbits. In this study, the identified candidate genes can be used as molecular markers for assisted selection in meat rabbits. At the same time, the inbreeding situation based on ROH assessment can provide reference for breeding and breeding preservation of meat rabbits.

Genomic Characterization of Local Croatian Sheep Breeds-Effective Population Size, Inbreeding & Signatures of Selection

The Istrian (IS) and the Pag sheep (PS) are local Croatian breeds which provide significant income for the regional economy and have a cultural and traditional importance for the inhabitants. The aim of this study was to estimate some important population specific genetic parameters in IS (N = 1293) and PS (N = 2637) based on genome wide SNPs. Estimates of linkage disequilibrium effective population size (Ne) evidenced more genetic variability in PS (Ne = 838) compared to IS (Ne = 197), regardless of historical time (both recent and ancient genetic variability). The discrepancy in the recent genetic variability between these breeds was additionally confirmed by the estimates of genomic inbreeding (FROH), which was estimated to be notably higher in IS (FROH>2 = 0.062) than in PS (FROH>2 = 0.029). The average FROH2-4, FROH4-8, FROH8-16, and FROH>16 were 0.26, 1.65, 2.14, and 3.72 for IS and 0.22, 0.61, 0.75, and 1.58 for PS, thus evidencing a high contribution of recent inbreeding in the overall inbreeding. One ROH island with > 30% of SNP incidence in ROHs was detected in IS (OAR6; 34,253,440-38,238,124 bp) while there was no ROH islands detected in PS. Seven genes (CCSER1, HERC3, LCORL, NAP1L5, PKD2, PYURF, and SPP1) involved in growth, feed intake, milk production, immune responses, and resistance were associated with the found autozygosity. The results of this study represent the first comprehensive insight into genomic variability of these two Croatian local sheep breeds and will serve as a baseline for setting up the most promising strategy of genomic Optimum Contribution Selection.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

Selection of Male Donors in Local Chicken Breeds to Implement the Italian Semen Cryobank: Variability in Semen Quality, Freezability and Fertility

Native breed conservation is an important component of poultry biodiversity. The aim of this work is to describe different steps that lead to donor selection for the implementation of the Italian Semen Cryobank of Autochthonous Chicken and Turkey Breeds. The variability within and between breeds was evaluated, and the stored semen reproductive capacity was in vivo tested using artificial insemination. Semen from Bionda Piemontese, Bianca di Saluzzo and Pepoi roosters was collected and processed. Concentration, volume, sperm membrane integrity, total motile sperm, progressive motile sperm and kinetic parameters were analyzed; sperm parameters accounting for bird variability were used to select male donors. Fresh semen quality parameters measured in donor ejaculates showed significant differences between breeds; no differences were found after cryopreservation. Variability in the fertilizing ability of cryopreserved semen was found within a breed (5-16%) and between birds within a breed (BP = 3-7%; BS = 7-31%; PP = 6-22%); only sperm quality parameters measured in fresh ejaculates, not frozen/thawed, may be associated with in vivo fertility results. In conclusion, sperm concentration and progressive motility were successfully used as selection parameters to identify chicken male donors with improved sperm quality for sperm cryobanking. However, new reliable sperm markers to predict cryopreserved semen's fertilizing ability are required.

Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle

Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

Distribution of Runs of Homozygosity and Their Relationship with Candidate Genes for Productivity in Kazakh Meat-Wool Sheep Breed

Increasing the fertility of sheep remains one of the crucial issues of modern sheep breeding. The Kazakh meat-wool sheep is an excellent breed with high meat and wool productivity and well adapted to harsh conditions. Nowadays, runs of homozygosity (ROHs) are considered a suitable approach for studying the genetic characteristics of farm animals. The aims of the study were to analyze the distribution of ROHs, describe autozygosity, and detect genomic regions with high ROH islands. In this study, we genotyped a total of 281 Kazakh meat-wool sheep using the Illumina iScan® system (EquipNet, Canton, MA, USA) via Ovine SNP50 BeadChip array. As a results, a total of 15,069 ROHs were found in the three Kazakh meat-wool sheep populations. The mean number of ROH per animal across populations varied from 40.3 (POP1) to 42.2 (POP2) in the category 1+ Mb. Furthermore, the number of ROH per animal in ROH1-2 Mb were much higher than ROH2-4 Mb and ROH8-16 Mb in the three sheep populations. Most of individuals had small number of ROH>16 Mb. The highest and lowest genomic inbreeding coefficient values were observed in POP2 and POP3, respectively. The estimated FROH presented the impact that recent inbreeding has had in all sheep populations. Furthermore, a set of interesting candidate genes (BMP2, BMPR2, BMPRIB, CLOCK, KDM2B, TIAM1, TASP1, MYBPC1, MYOM1, and CACNA2D1), which are related to the productive traits, were found. Collectively, these findings will contribute to the breeding and conservation strategies of the Kazakh meat-wool sheep breed.

A Comprehensive Genomic Analysis of Chinese Indigenous Ningxiang Pigs: Genomic Breed Compositions, Runs of Homozygosity, and Beyond

Ningxiang pigs are a renowned indigenous pig breed in China, known for their meat quality, disease resistance, and environmental adaptability. In recent decades, consumer demand for meats from indigenous breeds has grown significantly, fueling the selection and crossbreeding of Ningxiang pigs (NXP). The latter has raised concerns about the conservation and sustainable use of Ningxiang pigs as an important genetic resource. To address these concerns, we conducted a comprehensive genomic study using 2242 geographically identified Ningxiang pigs. The estimated genomic breed composition (GBC) suggested 2077 pigs as purebred Ningxiang pigs based on a ≥94% NXP-GBC cut-off. The remaining 165 pigs were claimed to be crosses, including those between Duroc and Ningxiang pigs and between Ningxiang and Shaziling pigs, and non-Ningxiang pigs. Runs of homozygosity (ROH) were identified in the 2077 purebred Ningxiang pigs. The number and length of ROH varied between individuals, with an average of 32.14 ROH per animal and an average total length of 202.4 Mb per animal. Short ROH (1-5 Mb) was the most abundant, representing 66.5% of all ROH and 32.6% of total ROH coverage. The genomic inbreeding estimate was low (0.089) in purebred Ningxiang pigs compared to imported western pig breeds. Nine ROH islands were identified, pinpointing candidate genes and QTLs associated with economic traits of interest, such as reproduction, carcass and growth traits, lipid metabolism, and fat deposition. Further investigation of these ROH islands and candidate genes is anticipated to better understand the genomics of Ningxiang pigs.

A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds

BACKGROUND

To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones.

RESULTS

The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses.

CONCLUSIONS

To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

The coefficients of inbreeding revealed by ROH study among inbred individuals belonging to each type of the first cousin marriage: A preliminary report from North India

BACKGROUND

Genome-wide runs of homozygosity (ROH) are appropriate to estimate genomic inbreeding, determine population history, unravel the genetic architecture of complex traits and disorders.

OBJECTIVE

The study sought to investigate and compare the actual proportion of homozygosity or autozygosity in the genomes of progeny of four subtypes of first cousin mating in humans, using both pedigree and genomic measures for autosomes and sex chromosomes.

METHODS

For this purpose, Illumina Global Screening Array-24 v1.0 BeadChip followed by cyto-ROH analysis through Illumina Genome Studio was used to characterise the homozygosity in five participants from North Indian state (Uttar Pradesh). PLINK v.1.9 software was used to estimate the genomic inbreeding coefficients viz. ROH-based inbreeding estimate (F_ROH) and homozygous loci-based inbreeding estimate (F_HOM).

RESULTS

A total of 133 ROH segments were detected with maximum number and genomic coverage in Matrilateral Parallel (MP) type and minimum in outbred individual. ROH pattern revealed that MP type has a higher degree of homozygosity than other subtypes. The comparison of F_ROH, F_HOM, and pedigree-based inbreeding estimate (F_PED) showed some difference in theoretical and realised proportion of homozygosity for sex-chromosomal loci but not for autosome for each type of consanguinity.

CONCLUSIONS

This is the very first study to compare and estimate the pattern of homozygosity among the kindreds of first cousin unions. However, a greater number of individuals from each type of marriage is required for statistical inference of no difference between theoretical and realized homozygosity among different degrees of inbreeding prevalent in humans worldwide.

Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep

BACKGROUND

The detection of selective traits in different populations can not only reveal current mechanisms of artificial selection for breeding, but also provide new insights into phenotypic variation in new varieties and the search for genes associated with important traits. Panou sheep is a cultivated breed of Tibetan sheep in China with stable genetic performance, consistent appearance and fast growth and development after decades of artificial selection and cultivation. Due to long-term adaptation to the high altitude, cold and hypoxic environment in the plateau area, they may have formed a unique gene pool that is different from other Tibetan sheep breeds. To explore the genetic resources of Panou sheep, we used next-generation sequencing technology for the first time to investigate the genome-wide population structure, genetic diversity, and candidate signatures of positive selection in Panou sheep.

RESULTS

Comparative genomic analysis with the closely related species Oula sheep (a native breed of Tibetan sheep in China) was used to screen the population selection signal of Panou sheep. Principal component analysis and neighbor joining tree showed that Panou sheep and Oula sheep had differences in population differentiation. Furthermore, analyses of population structure, they came from the same ancestor, and when K = 2, the two populations could be distinguished. Panou sheep exhibit genetic diversity comparable to Oula sheep, as shown by observed heterozygosity, expected heterozygosity and runs of homozygosity. Genome-wide scanning using the Fst and π ratio methods revealed a list of potentially selected related genes in Panou sheep compared to Oula sheep, including histone deacetylase 9 (HDAC9), protein tyrosine kinase 2 (PTK2), microphthalmia-related transcription factor (MITF), vesicular amine transporter 1 (VAT1), trichohyalin-like 1 (TCHHL1), amine oxidase, copper containing 3 (AOC3), interferon-inducible protein 35 (IFI35).

CONCLUSIONS

The results suggest that traits related to growth and development and plateau adaptation may be selection targets for the domestication and breeding improvement of Tibetan sheep. This study provides the fundamental footprints for Panou sheep breeding and management.

Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (F_PED ), or runs of homozygosity (ROH; F_ROH ) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of F_ROH . The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.

High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context

Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise F_ST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and F_ST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.

Genomic Tools for the Characterization of Local Animal Genetic Resources: Application in Mascaruna Goat

Italy contains a large number of local goat populations, some of which do not have a recognized genetic structure. The "Mascaruna" is a goat population reared for milk production in Sicily. In this study, a total of 72 individuals were genotyped with the Illumina Goat_IGGC_65K_v2 BeadChip with the aim to characterize the genetic diversity, population structure and relatedness with another 31 Italian goat populations. The results displayed a moderate level of genetic variability for Mascaruna, in concordance with the estimated values for Italian goats. Runs of homozygosity islands are linked to genes involved in milk production, immune response and local adaptation. Population structure analyses separated Mascaruna from the other goat populations, indicating a clear genetic differentiation. Although they are not conclusive, our current results represent a starting point for the creation of monitoring and conservation plans. Additional analyses and a wider sampling would contribute to refine and validate these results. Finally, our study describing the diversity and structure of Mascaruna confirms the usefulness of applied genomic analyses as valid tools for the study of the local uncharacterized genetic resources.

Genetic Structure of the Ca Rater Mallorquí Dog Breed Inferred by Microsatellite Markers

Ca Rater Mallorquí is a dog breed from the Island of Mallorca (Spain) traditionally used as a hunting and ratting dog to prevent disease spread and economic losses related to rodent activities on farms. However, the census data shows a population decline that should be addressed by implementing a conservation program. The first step to implementing a conservation plan is knowing the genetic situation of the Ca Rater Mallorquí population. Therefore, we aimed to genetically characterise the breed in our study. We analysed 33 microsatellites recommended by the International Society of Animal Genetics (ISAG) in 77 samples. Data were obtained from 13 samples of Balearic, Spanish, and international dog breeds to study the genetic diversity among breeds. The population did not significantly deviate from the Hardy-Weinberg equilibrium with heterozygosity (Ho) of 0.655 and expected heterozygosity (He) of 0.685. The Wright's fixation indices, the Factorial Correspondence Analysis (FCA), a dendrogram representing Reynolds genetic distance between populations, and the pairwise F_ST values establish the Ca Rater Mallorquí as an independent breed distinct from the Balearic, Spanish, and international breeds.

How Geography and Climate Shaped the Genomic Diversity of Italian Local Cattle and Sheep Breeds

Simple Summary

In this paper, we study the inter-relationships among geography, climate, and genetics in Italian local cattle and sheep breeds. In terms of genetic diversity, geography (latitude and longitude) appears to play a larger role in sheep (26.4%) than that in cattle (13.8%). Once geography is accounted for, 10.1% of cattle genomic diversity and 13.3% of that of sheep are attributable to climatic effects. Stronger geographic effects in sheep can be related to a combination of higher predomestication genetic variability together with biological and productive specializations. The climate alone seems to have had less impact on the current genetic diversity in both species even if climate and geography are greatly confounded. Results confirm that both species are the result of complex evolutionary histories triggered by interactions between human needs and environmental conditions.

Abstract

Understanding the relationships among geography, climate, and genetics is increasingly important for animal farming and breeding. In this study, we examine these inter-relationships in the context of local cattle and sheep breeds distributed along the Italian territory. To this aim, we used redundancy analysis on genomic data from previous projects combined with geographical coordinates and corresponding climatic data. The effect of geographic factors (latitude and longitude) was more important in sheep (26.4%) than that in cattle (13.8%). Once geography had been partialled out of analysis, 10.1% of cattle genomic diversity and 13.3% of that of sheep could be ascribed to climatic effects. Stronger geographic effects in sheep can be related to a combination of higher pre-domestication genetic variability together with biological and productive specificities. Climate alone seems to have had less impact on current genetic diversity in both species, even if climate and geography are greatly confounded. Results confirm that both species are the result of complex evolutionary histories triggered by interactions between human needs and environmental conditions.

Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China

Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of F_ROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.

Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems

Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from −0.009 to 0.029. A low (−0.12–0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (&lt;4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5–1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.

High-Density Genomic Characterization of Native Croatian Sheep Breeds

A recent comprehensive genomic analysis based on 50K SNP profiles has shown that the regional Balkan sheep populations have considerable genetic overlap but are distinctly different from surrounding breeds. All eight Croatian sheep breeds were represented by a small number of individuals per breed. Here, we genotyped 220 individuals representing the native Croatian sheep breeds (Istrian Sheep, Krk Island Sheep, Cres Island Sheep, Rab Island Sheep, Lika Pramenka, Pag Island Sheep, Dalmatian Pramenka, Dubrovnik Sheep) and mouflon using the Ovine Infinium® HD SNP BeadChip (606,006 SNPs). In addition, we included publicly available Balkan Pramenka and other Mediterranean sheep breeds. Our analyses revealed the complex population structure of Croatian sheep breeds and their origin and geographic barriers (island versus mainland). Migration patterns confirmed the historical establishment of breeds and the pathways of gene flow. Inbreeding coefficients (FROH&gt;2 Mb) between sheep populations ranged from 0.025 to 0.070, with lower inbreeding coefficients observed in Dalmatian Pramenka and Pag Island Sheep and higher inbreeding in Dubrovnik sheep. The estimated effective population size ranged from 61 to 1039 for Krk Island Sheep and Dalmatian Pramenka, respectively. Higher inbreeding levels and lower effective population size indicate the need for improved conservation management to maintain genetic diversity in some breeds. Our results will contribute to breeding and conservation strategies of native Croatian sheep breeds.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data

Background

Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (F_ROH) and homozygosity (F_HOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein.

Results

The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10–20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. F_ROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, F_ROH and F_HOM values were highly correlated (0.952–0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits.

Conclusion

Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08583-9.

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.

Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern

Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome‐wide single‐nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model‐based clustering and Neighbour‐Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (H_e) and contemporary effective population size (cNe) were reported in Siciliano (H_e = 0.323, cNe = 397), and the lowest were reported in Purosangue Orientale Siciliano (H_e = 0.277, cNe = 10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome‐wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.

Genomic analysis of 10 years of artificial selection in community‐based breeding programs in two Ethiopian indigenous sheep breeds

In recent times, community-based breeding programs (CBBPs) have been advocated as the best strategy for genetic improvement of local breeds in smallholder farms in developing countries. Since 2009, CBBPs have been implemented for Ethiopian Bonga and Menz sheep to improve growth rates resulting in significant genetic gains in 6-month weights. With the hypothesis that selection could be impacting their genomes, we systematically screened for possible genome changes in the two breeds by analyzing 600K BeadChip genotype data of 151 individuals (with the highest breeding values for 6-month weights) from CBBP flocks against 98 individuals from non-CBBP flocks. We observed no differences in genetic diversity and demographic dynamics between CBBP and non-CBBP flocks. Selection signature analysis employing ROH, logistic regression genome-wide association study , F_ST , XP-EHH and iHS revealed 5 (Bonga) and 11 (Menz) overlapping regions under selection, that co-localized with QTLs for production (body size/weight, growth, milk yield), meat/milk quality, and health/parasite resistance, suggesting that the decade-long selection has likely started to impact their genomes. However, genome-wide genetic differentiation between the CBBP and non-CBBP flocks is not yet clearly evident.

Sequencing Reveals Population Structure and Selection Signatures for Reproductive Traits in Yunnan Semi-Fine Wool Sheep (Ovis aries)

Yunnan semi-fine wool sheep are among the most important cultivated sheep breeds in China. However, their population structure, genetic characteristics and traits of interest are poorly studied. In this study, we systematically studied the population characteristics and selection signatures of 40 Yunnan semi-fine wool sheep using SNPs obtained from whole-genome resequencing data. A total of 1393 Gb of clean data were acquired. The mapping rate against the reference genome was 91.23% on average (86.01%–92.26%), and the average sequence depth was 9.51X. After filtering, 28,593,198 SNPs and 4,725,259 indels with high quality were obtained. The heterozygosity rate, inbreeding coefficient and effective population size of the sheep were calculated to preliminarily explore their genetic characteristics. The average heterozygosity rate was 0.264, the average inbreeding coefficient was 0.0099, and the effective population size estimated from the heterozygote excess (HE) was 242.9. Based on the Tajima’s D and integrated haplotype score (iHS) approaches, 562 windows and 11,356 core SNPs showed selection signatures in the Yunnan semi-fine wool sheep population. After genome annotation and gene enrichment analysis, we found traces of early domestication in sensory organs, behavioural activity and the nervous system as well as adaptive changes in reproductive and wool traits under selection in this population. Some selected genes related to litter size, including FSHR, BMPR1B and OXT, were identified as being under selection. Specific missense mutations of the FSHR gene that differed from the reference genome were also identified in the population, and we found some SNP variations that may affect litter size. Our findings provide a theoretical basis for the conservation and utilization of Yunnan semi-fine wool sheep. Furthermore, our results reveal some changes common to sheep after domestication and provide a new opportunity to investigate the genetic variation influencing fecundity within a population evolving under artificial selection.

Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo

Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (F_ROH) ranged from 0.0201 to 0.0371. Pearson correlations between F_ROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.

Identification of Copy Number Variations and Genetic Diversity in Italian Insular Sheep Breeds

Copy number variants (CNVs) are one of the major contributors to genetic diversity and phenotypic variation in livestock. The aim of this work is to identify CNVs and perform, for the first time, a CNV-based population genetics analysis with five Italian sheep breeds (Barbaresca, Comisana, Pinzirita, Sarda, and Valle del Belìce). We identified 10,207 CNVs with an average length of 1.81 Mb. The breeds showed similar mean numbers of CNVs, ranging from 20 (Sarda) to 27 (Comisana). A total of 365 CNV regions (CNVRs) were determined. The length of the CNVRs varied among breeds from 2.4 Mb to 124.1 Mb. The highest number of shared CNVRs was between Comisana and Pinzirita, and only one CNVR was shared among all breeds. Our results indicated that segregating CNVs expresses a certain degree of diversity across all breeds. Despite the low/moderate genetic differentiation among breeds, the different approaches used to disclose the genetic relationship showed that the five breeds tend to cluster in distinct groups, similar to the previous studies based on single-nucleotide polymorphism markers. Gene enrichment was described for the 37 CNVRs selected, considering the top 10%. Out of 181 total genes, 67 were uncharacterized loci. Gene Ontology analysis showed that several of these genes are involved in lipid metabolism, immune response, and the olfactory pathway. Our results corroborated previous studies and showed that CNVs represent valuable molecular resources for providing useful information for separating the population and could be further used to explore the function and evolutionary aspect of sheep genome.

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.

Detection and Visualization of Heterozygosity-Rich Regions and Runs of Homozygosity in Worldwide Sheep Populations

In this study, we chose 17 worldwide sheep populations of eight breeds, which were intensively selected for different purposes (meat, milk, or wool), or locally-adapted breeds, in order to identify and characterize factors impacting the detection of runs of homozygosity (ROH) and heterozygosity-rich regions (HRRs) in sheep. We also applied a business intelligence (BI) tool to integrate and visualize outputs from complementary analyses. We observed a prevalence of short ROH, and a clear distinction between the ROH profiles across populations. The visualizations showed a fragmentation of medium and long ROH segments. Furthermore, we tested different scenarios for the detection of HRR and evaluated the impact of the detection parameters used. Our findings suggest that HRRs are small and frequent in the sheep genome; however, further studies with higher density SNP chips and different detection methods are suggested for future research. We also defined ROH and HRR islands and identified common regions across the populations, where genes related to a variety of traits were reported, such as body size, muscle development, and brain functions. These results indicate that such regions are associated with many traits, and thus were under selective pressure in sheep breeds raised for different purposes. Interestingly, many candidate genes detected within the HRR islands were associated with brain integrity. We also observed a strong association of high linkage disequilibrium pattern with ROH compared with HRR, despite the fact that many regions in linkage disequilibrium were not located in ROH regions.

Biodiversity of Russian Local Sheep Breeds Based on Pattern of Runs of Homozygosity

Russian sheep breeds traditionally raised in specific environments are valuable parts of sociocultural heritage and economic component of the regions. However, the import of commercial breeds negatively influences the population sizes of local sheep populations and might lead to biodiversity loss. Estimation of the runs of homozygosity (ROH) in local sheep genomes is an informative tool to address their current genetic state. In this work, we aimed to address the ROH distribution and to estimate genome inbreeding based on SNP data to evaluate genetic diversity in Russian local sheep breeds. Materials for this study included SNP-genotypes from twenty-seven Russian local sheep breeds which were generated using the Illumina OvineSNP50 BeadChip (n = 391) or the Illumina Ovine Infinium HD BeadChip (n = 315). A consecutive runs method was used to calculate ROH which were estimated for each animal and then categorized in the ROH length classes. The ROH were found in all breeds. The mean ROH length varied from 86 to 280 Mb, while the ROH number ranged from 37 to 123. The genomic inbreeding coefficient varied from 0.033 to 0.106. Our findings provide evidence of low to moderate genomic inbreeding in major local sheep populations.

Grid search approach to discriminate between old and recent inbreeding using phenotypic, pedigree and genomic information

Background

Although inbreeding caused by the mating of animals related through a recent common ancestor is expected to have more harmful effects on phenotypes than ancient inbreeding (old inbreeding), estimating these effects requires a clear definition of recent (new) and ancient (old) inbreeding. Several methods have been proposed to classify inbreeding using pedigree and genomic data. Unfortunately, these methods are largely based on heuristic criteria such as the number of generations from a common ancestor or length of runs of homozygosity (ROH) segments. To mitigate these deficiencies, this study aimed to develop a method to classify pedigree and genomic inbreeding into recent and ancient classes based on a grid search algorithm driven by the assumption that new inbreeding tends to have a more pronounced detrimental effect on traits. The proposed method was tested using a cattle population characterized by a deep pedigree.

Results

Effects of recent and ancient inbreeding were assessed on four growth traits (birth, weaning and yearling weights and average daily gain). Thresholds to classify inbreeding into recent and ancient classes were trait-specific and varied across traits and sources of information. Using pedigree information, inbreeding generated in the last 10 to 11 generations was considered as recent. When genomic information (ROH) was used, thresholds ranged between four to seven generations, indicating, in part, the ability of ROH segments to characterize the harmful effects of inbreeding in shorter periods of time. Nevertheless, using the proposed classification method, the discrimination between new and old inbreeding was less robust when ROH segments were used compared to pedigree. Using several model comparison criteria, the proposed approach was generally better than existing methods. Recent inbreeding appeared to be more harmful across the growth traits analyzed. However, both new and old inbreeding were found to be associated with decreased yearling weight and average daily gain.

Conclusions

The proposed method provided a more objective quantitative approach for the classification of inbreeding. The proposed method detected a clear divergence in the effects of old and recent inbreeding using pedigree data and it was superior to existing methods for all analyzed traits. Using ROH data, the discrimination between old and recent inbreeding was less clear and the proposed method was superior to existing approaches for two out of the four analyzed traits. Deleterious effects of recent inbreeding were detected sooner (fewer generations) using genomic information than pedigree. Difference in the results using genomic and pedigree information could be due to the dissimilarity in the number of generations to a common ancestor. Additionally, the uncertainty associated with the identification of ROH segments and associated inbreeding could have an effect on the results. Potential biases in the estimation of inbreeding effects may occur when new and old inbreeding are discriminated based on arbitrary thresholds. To minimize the impact of inbreeding, mating designs should take the different inbreeding origins into consideration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07872-z.

Genome-wide insights on gastrointestinal nematode resistance in autochthonous Tunisian sheep

Gastrointestinal nematode (GIN) infections have negative impacts on animal health, welfare and production. Information from molecular studies can highlight the underlying genetic mechanisms that enhance host resistance to GIN. However, such information often lacks for traditionally managed indigenous livestock. Here, we analysed 600 K single nucleotide polymorphism genotypes of GIN infected and non-infected traditionally managed autochthonous Tunisian sheep grazing communal natural pastures. Population structure analysis did not find genetic differentiation that is consistent with infection status. However, by contrasting the infected versus non-infected cohorts using ROH, LR-GWAS, F_ST and XP-EHH, we identified 35 candidate regions that overlapped between at least two methods. Nineteen regions harboured QTLs for parasite resistance, immune capacity and disease susceptibility and, ten regions harboured QTLs for production (growth) and meat and carcass (fatness and anatomy) traits. The analysis also revealed candidate regions spanning genes enhancing innate immune defence (SLC22A4, SLC22A5, IL-4, IL-13), intestinal wound healing/repair (IL-4, VIL1, CXCR1, CXCR2) and GIN expulsion (IL-4, IL-13). Our results suggest that traditionally managed indigenous sheep have evolved multiple strategies that evoke and enhance GIN resistance and developmental stability. They confirm the importance of obtaining information from indigenous sheep to investigate genomic regions of functional significance in understanding the architecture of GIN resistance.

Genome-wide scan for runs of homozygosity identifies candidate genes in Wannan Black pigs

Objective

Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that can reveal inbreeding levels, selection pressure, and mating schemes. In this study, ROHs were evaluated in Wannan Black pigs to assess the inbreeding levels and the genome regions with high ROH frequency.

Methods

In a previous study, we obtained 501.52 GB of raw data from resequencing (10×) of the genome and identified 21,316,754 single-nucleotide variants in 20 Wannan Black pig samples. We investigated the number, length, and frequency of ROH using resequencing data to characterize the homozygosity in Wannan Black pigs and identified genomic regions with high ROH frequencies.

Results

In this work, 1,813 ROHs (837 ROHs in 100 to 500 kb, 449 ROHs in 500 to 1,000 kb, 527 ROHs in >1,000 kb) were identified in all samples, and the average genomic inbreeding coefficient (F_ROH) in Wannan Black pigs was 0.5234. Sixty-one regions on chromosomes 2, 3, 7, 8, 13, 15, and 16 harbored ROH islands. In total, 105 genes were identified in 42 ROH islands, among which some genes were related to production traits.

Conclusion

This is the first study to identify ROH across the genome of Wannan Black pigs, the Chinese native breed of the Anhui province. Overall, Wannan Black pigs have high levels of inbreeding due to the influence of ancient and recent inbreeding due to the genome. These findings are a reliable resource for future studies and contribute to save and use the germplasm resources of Wannan Black pigs.

The pattern of runs of homozygosity and genomic inbreeding in world-wide sheep populations

Genome-wide pattern of runs of homozygosity (ROH) across ovine genome can provide a useful resource for studying diversity and demography history in sheep. We analyzed 50 k SNPs chip data of 2536 animals to identify pattern, distribution and level of ROHs in 68 global sheep populations. A total of 60,301 ROHs were detected in all breeds. The majority of the detected ROHs were <16 Mb and the average total number of ROHs per individual was 23.8 ± 13.8. The ROHs greater than 1 Mb covered on average 8.2% of the sheep autosomes, 1% of which was related to the ROHs with 1-4 Mb of length. The mean sum of ROH length in two-thirds of the populations was less than 250 Mb ranging from 21.7 to near 570 Mb. The level of genomic inbreeding was relatively low. The average of the inbreeding coefficients based on ROH (F_ROH) was 0.09 ± 0.05. It was rising in a stepwise manner with distance from Southwest Asia and maximum values were detected in North European breeds. A total of 465 ROH hotspots were detected in 25 different autosomes which partially surrounding 257 Refseq genes across the genome. Most of the detected genes were related to growth, body weight, meat production and quality, wool production and pigmentation. In conclusion, our analysis showed that the sheep genome, compared with other livestock species such as cattle and pig, displays low levels of homozygosity and appropriate genetic diversity for selection response and genetic merit gain.

Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia

The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.

Genome-wide analysis reveals the patterns of genetic diversity and population structure of 8 Italian local chicken breeds

The aim of this study was to conduct a genome-wide comparative analysis of 8 local Italian chicken breeds (Ermellinata di Rovigo, Millefiori di Lonigo [PML], Polverara Bianca, Polverara Nera, Padovana, Pepoi [PPP], Robusta Lionata, and Robusta Maculata), all under a conservation plan, to understand their genetic diversity and population structure. A total of 152 animals were analyzed using the Affymetrix Axiom 600 K Chicken Genotyping Array. The levels of genetic diversity were highest and lowest in PML and PPP, respectively. The results of genomic inbreeding based on runs of homozygosity (ROH; FROH) showed marked differences among breeds and ranged from 0.161 (PML) to 0.478 (PPP). Furthermore, in all breeds, short ROH (<4 Mb in length) were more frequent than long segments. Patterns of genetic differentiation, model-based clustering, and neighbor networks showed that most breeds formed nonoverlapping clusters and were clearly separate populations. The 2 Polverara breeds shared a similar genetic background and showed the lowest genetic differentiation in comparison with purebred lines; the local populations showed separated groups. PPP and PML were closer to the group of the purebred broiler lines (BRSA, BRSB, BRDA, and BRDB). Six genomic regions are presented as hotspots of autozygosity among the Italian chicken breeds, with candidate genes involved in multiple morphological phenotypes as breast muscle, muscle dry matter content, and body weight. This study is the first exhaustive genome-wide analysis of the diversity of these Italian local chickens from Veneto region. We conclude that breeds have conserved authentic genetic patterns. The results are of significant importance because they will help design and implement conservation strategies. In fact, the conservation of these breeds may also have positive impacts on the local economy, niche traditional markets, and offering a source of high-quality products to consumers. In this context, genomic information may play a crucial role in the management of local breeds.

A selection signatures study among Middle Eastern and European sheep breeds

Selection, both natural and artificial, leaves patterns on the genome during domestication of animals and leads to changes in allele frequencies among populations. Detecting genomic regions influenced by selection in livestock may assist in understanding the processes involved in genome evolution and discovering genomic regions related to traits of economic and ecological interests. In the current study, genetic diversity analyses were conducted on 34,206 quality-filtered SNP positions from 450 individuals in 15 sheep breeds, including six indigenous breeds from the Middle East, namely Iranian Balouchi, Afshari, Moghani, Qezel, Karakas and Norduz, and nine breeds from Europe, namely East Friesian Sheep, Ile de France, Mourerous, Romane, Swiss Mirror, Spaelsau, Suffolk, Comisana and Engadine Red Sheep. The SNP genotype data generated by the Illumina OvineSNP50 Genotyping BeadChip array were used in this analysis. We applied two complementary statistical analyses, F_ST (fixation index) and xp-EHH (cross-population extended haplotype homozygosity), to detect selection signatures in Middle Eastern and European sheep populations. F_ST and xp-EHH detected 629 and 256 genes indicating signatures of selection, respectively. Genomic regions identified using F_ST and xp-EHH contained the CIDEA, HHATL, MGST1, FADS1, RTL1 and DGKG genes, which were reported earlier to influence a number of economic traits. Both F_ST and xp-EHH approaches identified 60 shared genes as the signatures of selection, including four candidate genes (NT5E, ADA2, C8A and C8B) that were enriched for two significant Gene Ontology (GO) terms associated with the adenosine metabolic procedure. Knowledge about the candidate genomic regions under selective pressure in sheep breeds may facilitate identification of the underlying genes and enhance our understanding on these genes role in local adaptation.