Analysis of Migration and Adaptive Evolution in Tibetan Sheep Populations

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 (<i>HIF1AN</i> [Hypoxia Inducible Factor 1 Alpha Subunit Inhibitor], <i>HBE1</i> [Hemoglobin Subunit Epsilon 1], <i>HBE2</i> [Hemoglobin Subunit Epsilon 2], <i>TNFAIP3</i> [TNF Alpha Induced Protein 3], <i>RAD50</i> [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 (<i>HIF1AN</i> [Hypoxia Inducible Factor 1 Alpha Subunit Inhibitor], <i>HBE1</i> [Hemoglobin Subunit Epsilon 1], <i>HBE2</i> [Hemoglobin Subunit Epsilon 2], <i>TNFAIP3</i> [TNF Alpha Induced Protein 3], <i>RAD50</i> [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 (<i>TP53</i> [Tumor Protein P53], <i>ATG101</i> [Autophagy Related 101], <i>ATP12A</i> [ATPase H+/K+ Transporting Non-Gastric Alpha2 Subunit], <i>KRT80</i> [Keratin 80], <i>KRT7</i> [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.