Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Wool fineness affects the quality of wool, and some studies have identified about forty candidate genes that affect sheep wool fineness, but these genes often reveal only a certain proportion of the variation in wool thickness. We further explore additional genes associated with the fineness of sheep wool. Whole-genome resequencing of eight sheep breeds was performed to reveal selection signals associated with wool fineness, including four coarse wool and four fine/semi-fine wool sheep breeds. Multiple methods to reveal selection signals (Fst and θπ Ratio and XP-EHH) were applied for sheep wool fineness traits. In total, 269 and 319 genes were annotated in the fine wool (F vs. C) group and the coarse wool (C vs. F) group, such as LGR4, PIK3CA, and SEMA3C and NFIB, OPHN1, and THADA. In F vs. C, 269 genes were enriched in 15 significant GO Terms (p < 0.05) and 38 significant KEGG Pathways (p < 0.05), such as protein localization to plasma membrane (GO: 0072659) and Inositol phosphate metabolism (oas 00562). In C vs. F, 319 genes were enriched in 21 GO Terms (p < 0.05) and 16 KEGG Pathways (p < 0.05), such as negative regulation of focal adhesion assembly (GO: 0051895) and Axon guidance (oas 04360). Our study has uncovered genomic information pertaining to significant traits in sheep and has identified valuable candidate genes. This will pave the way for subsequent investigations into related traits.

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, F_st & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

In recent years, research has gradually uncovered the mechanisms of animal adaptation to hypoxic conditions in different altitude environments, particularly at the genomic level. However, past genomic studies on high-altitude adaptation have often not delved deeply into the differences between varying altitude levels. This study conducted whole-genome sequencing on 60 Tibetan sheep (Medium Altitude Group (MA): 20 Tao sheep (TS) at 2887 m, High Altitude Group (HA): 20 OuLa sheep (OL) at 3501 m, and Ultra-High Altitude Group (UA): 20 AWang sheep (AW) at 4643 m) from different regions of the Tibetan Plateau in China to assess their responses under varying conditions. Population genetic structure analysis revealed that the three groups are genetically independent, but the TS and OL groups have experienced gene flow with other northern Chinese sheep due to geographical factors. Selection signal analysis identified FGF10, MMP14, SLC25A51, NDUFB8, ALAS1, PRMT1, PRMT5, and HIF1AN as genes associated with ultra-high-altitude hypoxia adaptation, while HMOX2, SEMA4G, SLC16A2, SLC22A17, and BCL2L2 were linked to high-altitude hypoxia adaptation. Functional analysis showed that ultra-high-altitude adaptation genes tend to influence physiological mechanisms directly affecting oxygen uptake, such as lung development, angiogenesis, and red blood cell formation. In contrast, high-altitude adaptation genes are more inclined to regulate mitochondrial DNA replication, iron homeostasis, and calcium signaling pathways to maintain cellular function. Additionally, the functions of shared genes further support the adaptive capacity of Tibetan sheep across a broad geographic range, indicating that these genes offer significant selective advantages in coping with oxygen scarcity. In summary, this study not only reveals the genetic basis of Tibetan sheep adaptation to different altitudinal conditions but also highlights the differences in gene regulation between ultra-high- and high-altitude adaptations. These findings offer new insights into the adaptive evolution of animals in extreme environments and provide a reference for exploring adaptation mechanisms in other species under hypoxic conditions.

Tibetan sheep have developed unique adaptations for survival in the Qinghai–Tibet Plateau environment. However, the functional genes and molecular mechanisms that regulate hypoxia adaptation have not been fully characterized. In this study, based on the whole-genome resequencing data for Tibetan sheep at different altitudes, the population differentiation index (F_ST) and nucleotide diversity ratio (θπ ratio) were evaluated in populations of 20 Oula sheep (3501 m altitude, OL), 20 Zashijia sheep (4369 m altitude, ZSJ), and 20 Awang sheep (4643 m altitude, AW) to reveal candidate loci related to high-altitude hypoxia. We screened 728 and 524 candidate genes in the AW vs. OL and ZSJ vs. OL groups, respectively, of which 134 genes were jointly screened. Candidate genes were mainly enriched in the Ras, melanoma, melanogenesis, VEGF, and PPAR signaling pathways. HIF1AN, PDGFA, PDGFD, ANXA2, SOCS2, NOXA1, WNT7B, MMP14, GNG2, ATF6, PGAM2, PPP3R1, GSTCD, and PPARA may play important roles in the high-altitude adaptation of Tibetan sheep. In conclusion, this study provides valuable insights into the genes and molecular mechanisms that underlie high-altitude hypoxia adaptation in Tibetan sheep.

The genetic basis for Tibetan sheep adaptation to different high-altitude environments remains unknown. This study conducted whole-genome resequencing on 80 Tibetan sheep individuals from four major distribution areas on the Qinghai–Tibet Plateau. Based on the high-quality single-nucleotide polymorphisms (SNPs) obtained, an analysis of population-level genomic selection signals was performed. Population genomic analysis revealed that Tibetan sheep distributed across China originated in northern China but showed evidence of gene flow from South Asian sheep. Between populations from extremely high-altitude and mid-altitude regions, selection analyses identified five strongly positive selected genes (HIF1AN [Hypoxia Inducible Factor 1 Alpha Subunit Inhibitor], HBE1 [Hemoglobin Subunit Epsilon 1], HBE2 [Hemoglobin Subunit Epsilon 2], TNFAIP3 [TNF Alpha Induced Protein 3], RAD50 [RAD50 Double Strand Break Repair Protein]). These genes are associated with adaptation to hypoxia and intense UV radiation in high-altitude environments. Selection analyses between populations from extremely high-altitude and mid-altitude regions identified five strongly selected genes (HIF1AN [Hypoxia Inducible Factor 1 Alpha Subunit Inhibitor], HBE1 [Hemoglobin Subunit Epsilon 1], HBE2 [Hemoglobin Subunit Epsilon 2], TNFAIP3 [TNF Alpha Induced Protein 3], RAD50 [RAD50 Double Strand Break Repair Protein]) associated with hypoxia and intense UV radiation in high-altitude environments. Comparative genomic analyses of populations in cold and arid environments identified several candidate genes related to energy and water homeostasis, as well as hair development (TP53 [Tumor Protein P53], ATG101 [Autophagy Related 101], ATP12A [ATPase H+/K+ Transporting Non-Gastric Alpha2 Subunit], KRT80 [Keratin 80], KRT7 [Keratin 7]). Additionally, Tibetan sheep in the high-altitude arid deserts exhibit stronger adaptive selection for energy homeostasis and water utilization; meanwhile, the HIF-1 [Hypoxia Inducible Factor 1] signaling pathway helps counteract oxidative stress induced by extreme water scarcity in the plateau environment. Our study supports the hypothesis that Tibetan sheep originated in northern China and identifies distinct adaptive features in the Tibetan sheep genome corresponding to their habitats.