Wild Yak Research Articles

Identifying low-density, ancestry-informative SNP markers through whole genome resequencing in Indian, Chinese, and wild yak

The current investigation was undertaken to elucidate the population-stratifying and ancestry-informative markers in Indian, Chinese, and wild yak populations using whole genome resequencing (WGS) analysis while employing various selection strategies (Delta, Pairwise Wright’s Fixation Index-FST, and Informativeness of Assignment) and marker densities (5–25 thousand). The study used WGS data on 105 individuals from three separate yak cohorts i.e., Indian yak (n = 29), Chinese yak (n = 61), and wild yak (n = 15). Variant calling in the GATK program with strict quality control resulted in 1,002,970 high-quality and independent (LD-pruned) SNP markers across the yak autosomes. Analysis was undertaken in toolbox for ranking and evaluation of SNPs (TRES) program wherein three different criteria i.e., Delta, Pairwise Wright’s Fixation Index-FST, and Informativeness of Assignment were employed to identify population-stratifying and ancestry-informative markers across various datasets. The top-ranked 5,000 (5K), 10,000 (10K), 15,000 (15K), 20,000 (20K), and 25,000 (25K) SNPs were identified from each dataset while their composition and performance was assessed using different criteria. The average genomic breed clustering of Indian, Chinese, and wild yak cohorts with full density dataset (105 individuals with 1,002,970 markers) was 81.74%, 80.02%, and 83.62%, respectively. Informativeness of Assignment criterion with 10K density emerged as the best combination for three yak cohorts with 86.94%, 96.46%, and 98.20% clustering for Indian, Chinese, and wild yak, respectively. There was an average increase of 7.56%, 22.72%, and 30.35% in genomic breed clustering scores of Indian, Chinese, and wild yak cohorts over the estimates of the original dataset. The selected markers showed overlap multiple protein-coding genes within a 10 kb window including ADGRB3, ANK1, CACNG7, CALN1, CHCHD2, CREBBP, GLI3, KHDRBS2, and OSBPL10. This is the first report ever on elucidating low-density SNP marker sets with population-stratifying and ancestry-informative properties in three yak groups using WGS data. The results gain significance for application of genomic selection using cost-effective low-density SNP panels in global yak species.

PSXIII-11 Large-scale profiling of non-conventional peptides in several tissues of Yak (Bos mutus) by peptidogenomic pipeline

PSXIII-11 Large-scale profiling of non-conventional peptides in several tissues of Yak (<i>Bos mutus</i>) by peptidogenomic pipeline

Comparative Analysis of Gut Bacterial Diversity in Wild and Domestic Yaks on the Qinghai-Tibetan Plateau.

The gut microbiota is a diverse and complex population, and it has a key role in the host's health and adaptability to the environment. The present study investigated the fecal bacterial community of wild grazing (WG) and domestic grazing (DG) yaks on natural grazing pastures, analyzing the gut microbiota using 16S rRNA sequencing to assess bacterial diversity. A total of 48 yak fecal samples were selected from two different grazing habitats. The DG group had more crude proteins and non-fiber carbohydrates. The WG group had more OM, insoluble dietary fiber such as NDF, ADF, ether extract, and TC. There were 165 and 142 unique operational taxonomic units (OTUs) in the WG and DG groups, respectively. Shannon index analysis revealed a higher bacterial diversity in the WG group than in the DG group. At the phylum level, Firmicutes were the dominant bacterial taxa in both groups. The relative abundance of Firmicutes in the WG group was higher than in the DG group. At the family level, the WG group had a significantly higher abundance of Ruminococcaceae (p < 0.001) and Rikenellaceae (p < 0.001) than the DG group. The abundances of Alloprevotella and Succinivibrio were more pronounced in the DG group than in the WG group at the genus level. This study presents a novel understanding of the bacterial communities of ruminants and their potential applications for livestock production.

Maternal Genetic Diversity, Differentiation and Phylogeny of Wild Yak and Four Domestic Yak Breeds in Qinghai, China Inferred from Mitochondrial Cytb Variations

Maternal Genetic Diversity, Differentiation and Phylogeny of Wild Yak and Four Domestic Yak Breeds in Qinghai, China Inferred from Mitochondrial Cytb Variations

Cauldrons of Bronze Age nomads reveals 2700 year old yak milk and the deep antiquity of food preparation techniques

Cauldrons, vessels that are simultaneously common and enigmatic, offer insights into past cultural and social traditions. While assumed to possess a special function, what these cauldrons contained is still largely mysterious. These vessels, such as those made from bronze or copper alloys, function as reservoirs for ancient organics through the antibacterial qualities provided by the metal surfaces. Here we show, through protein analysis, that cauldrons from the Final Bronze Age (ca. 2700 BP) were primarily used to collect blood from ruminants, primarily caprines, likely for the production of sausages in a manner similar to contemporary practices in Mongolia’s rural countryside. Our findings present a different function from the recent findings of cooked meat in copper-alloy vessels from the northern Caucasus 2000 years earlier, exposing the diversity in food preparation techniques. Our secondary findings of bovine milk within the cauldron, including peptides specific to Bos mutus, pushes back their regional domestication into the Bronze Age.

Tibetan terrestrial and aquatic ecosystems collapsed with cryosphere loss inferred from sedimentary ancient metagenomics.

Glacier and permafrost shrinkage and land-use intensification threaten mountain wildlife and affect nature conservation strategies. Here, we present paleometagenomic records of terrestrial and aquatic taxa from the southeastern Tibetan Plateau covering the last 18,000 years to help understand the complex alpine ecosystem dynamics. We infer that steppe-meadow became woodland at 14 ka (cal BP) controlled by cryosphere loss, further driving a herbivore change from wild yak to deer. These findings weaken the hypothesis of top-down control by large herbivores in the terrestrial ecosystem. We find a turnover in the aquatic communities at 14 ka, transitioning from glacier-related (blue-green) algae to abundant nonglacier-preferring picocyanobacteria, macrophytes, fish, and otters. There is no evidence for substantial effects of livestock herding in either ecosystem. Using network analysis, we assess the stress-gradient hypothesis and reveal that root hemiparasitic and cushion plants are keystone taxa. With ongoing cryosphere loss, the protection of their habitats is likely to be of conservation benefit on the Tibetan Plateau.

Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens)

BackgroundThere is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity.ResultsA total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks.ConclusionsYaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

Nuclear sequences of mitochondrial origin in domestic yak

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs). Here, we analysed, for the first time, NUMTs in the domestic yak genome. We obtained 499 alignment matches covering 340.2 kbp of the yak nuclear genome. After a merging step, we identified 167 NUMT regions with a total length of ~ 503 kbp, representing 0.02% of the nuclear genome. We discovered copies of all mitochondrial regions and found that most NUMT regions are intergenic or intronic and mostly untranscribed. 98 different NUMT regions from domestic yak showed high homology with cow and/or wild yak genomes, suggesting selection or hybridization between domestic/wild yak and cow. To rule out the possibility that the identified NUMTs could be artifacts of the domestic yak genome assembly, we validated experimentally five NUMT regions by PCR amplification. As NUMT regions show high similarity to the mitochondrial genome can potentially pose a risk to domestic yak DNA mitochondrial studies, special care is therefore needed to select primers for PCR amplification of mitochondrial DNA sequences.

Meeting your ancestors – Sticks, stones, and discord in Earth’s outposts

An unusual problem at the extreme edges of Earth involves marginalized people and associated livelihoods when their domestic stock is lost to wild ancestors. Animosity to wildlife results due to economic costs, personal injury, and death. Despite IUCN recognition of biodiversity at a policy front, justice and equity for humans have not been well integrated in these realms. In particular, Tibetans, Uyghurs, Changpas, Nepalis, Bhutanese, and Brokpa, as well as Mongolians in the Gobi Desert deal with these realities when the wild progenitors of their livestock attempt to appropriate genetically-related domestic females to biologically propagate. Among the highly endangered ancestors are the likes of wild camels, wild yaks, banteng, and gaur. Beyond Asia, some similarity of issues still thwart African, and South and North American pastoralists and involve species from caribou and bison to wild sheep and goats, and guanacos. Biologically reality has trumped our evolutionary innovation to benefit ourselves through domestication because of the powerful mutual attraction between descendent domestics and their wild ancestors. Conservation biologists and on-the-ground practitioners must engage more with local pastoralists to tackle these complex and infrequently described conflicts on how best to implement protective policies. Sticks and stones are neither the answer for pastoral safety nor for avoiding introgression among these iconic mammals.

The spatio-temporal expression analysis of parathyroid hormone like hormone gene provides a new insight for bone growth of the antler tip tissue in sika deer.

Parathyroid hormone like hormone (PTHLH), as an essential factor for bone growth, is involved in a variety of physiological processes. The aim of this study was to explore the role of PTHLH gene in the growth of antlers. The coding sequence (CDS) of PTHLH gene cDNA was obtained by cloning in sika deer (Cervus nippon), and the bioinformatics was analyzed. The quantitative real-time polymerase chain reaction (qRT-PCR) was used to analyze the differences expression of PTHLH mRNA in different tissues of the antler tip at different growth periods (early period, EP; middle period, MP; late period, LP). The CDS of PTHLH gene was 534 bp in length and encoded 177 amino acids. Predictive analysis results revealed that the PTHLH protein was a hydrophilic protein without transmembrane structure, with its secondary structure consisting mainly of random coil. The PTHLH protein of sika deer had the identity of 98.31%, 96.82%, 96.05%, and 94.92% with Cervus canadensis, Bos mutus, Oryx dammah and Budorcas taxicolor, which were highly conserved among the artiodactyls. The qRT-PCR results showed that PTHLH mRNA had a unique spatio-temporal expression pattern in antlers. In the dermis, precartilage, and cartilage tissues, the expression of PTHLH mRNA was extremely significantly higher in MP than in EP, LP (p<0.01). In the mesenchyme tissue, the expression of PTHLH mRNA in MP was significantly higher than that of EP (p<0.05), but extremely significantly lower than that of LP (p<0.01). The expression of PTHLH mRNA in antler tip tissues at all growth periods had approximately the same trend, that is, from distal to basal, it was first downregulated from the dermis to the mesenchyme and then continuously up-regulated to the cartilage tissue. PTHLH gene may promote the rapid growth of antler mainly through its extensive regulatory effect on the antler tip tissue.

Mapping conservation priorities for wild yak (Bos mutus) habitats on the Tibetan Plateau, China

The wild yak (Bos mutus) is a cold-tolerant herbivore native to the Tibetan Plateau and has been categorized as vulnerable by the International Union for Conservation of Nature and Natural Resources. Low population densities within currently fragmented habitats and unclear landscape conservation priorities warrant attention. Herein, we employed the maximum entropy (MaxEnt) model using over 900 wild yak occurrence records to model wild yak habitat suitability. Our analysis revealed unprotected wild yak landscapes covering 30.79 % of the habitat area, indicating a conservation gap between protected areas (PAs) and wild yak habitats. To protect metapopulation dynamics and mitigate high risks of poaching, habitat degradation and fragmentation, resource competition, and degenerated genetic characterization of wild yaks in fragmented and degraded habitat, we identified eight habitat patches as landscape conservation units (LCUs) and 14 linkages among the LCUs, enhancing the connectivity between LCUs to decrease negative effects of genetic threats. A centrality analysis demonstrated that Changtang, Arjinshan, and Hoh Xil national nature reserves and their linkages are all critical for the maintenance of habitat connectivity. Here, we suggest that habitat- and LCU-specific conservation strategies should be highlighted during the establishment of PAs and transboundary cooperation. Ultimately, our results can assist conservationists and land managers in comprehending wild yak distribution, movement, and habitat requirements, as well as for the development of effective protection strategies. Furthermore, the combined modeling method (MaxEnt–Zonation–InVEST) could be utilized as a component for identifying conservation priorities and linkages between core patches for species and assessing the efficiency of PAs, core habitats, and corridors in achieving conservation goals. Our study can provide a framework in identifying priority conservation and connectivity between habitat patches to facilitate effectively conservation and genetic resilience for endangered species in fragmented habitats.

Influence of environmental factors on physiological state of domestic yak (Bos mutus) (on the example of the Steppe Research Station “Orenburg Tarpania”)

Increased interest in the treatment and prevention of diseases of large herbivorous animals is associated with the processes of their acclimatisation during introduction and reintroduction, modelling behaviour in the conditions of existence in a wild natural environment, with minimal human participation. Since 2014, acclimatisation of domestic yak (Bos mutus Przewalski, 1883) in the steppe zone in the east of the Russian Plain, at the Steppe Research Station “Orenburg Tarpania” has been carried out. “Orenburg Tarpania” is a project of the Steppe Institute of the Ural Branch of the Russian Academy of Sciences to return large ungulates to their natural habitat in order to maintain optimal functioning of natural ecosystems. The ideology of the project corresponds to the concept of ecological rewilding. The animals are kept in a semi-voluntary state. The domestic yak is chosen as a species adapted to existence in extreme conditions with low temperature, extreme solar radiation in the southern latitudes, and relatively arid conditions. The conditions of yak keeping in the steppe zone differ from its natural habitat. During the period of observation, two reasons for the disturbance of the physical condition of animals have been recorded. Firstly, these are injuries received in the process of life activity. Injuries are often complicated by wolfarthiasis. The second disturbance of physical condition of domestic yak is eye damage in some individuals, which occurs due to incomplete adaptation of some animals to high wind activity, typical for the steppe region. The main cause of conjunctivitis and keratitis is dust raised by strong winds. A set of measures has been developed and adapted to provide treatment of animals in conditions close to natural habitat and semivoluntary confinement.

Evolutionary origin of genomic structural variations in domestic yaks

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

Genome Variation Map of Domestic Qinghai-Tibet Plateau Yaks by SLAF-Seq Reveals Genetic Footprint during Artificial Selection.

The yak (Bos grunniens) was domesticated in the high-altitude QTP. Research about their genetic diversity and population structure is limited. In this study, we resequenced the genome of 494 domestic yaks using Specific-Locus Amplified Fragment Sequencing (SLAF-seq). The survey was conducted on six populations sampled from isolated locations in China in order to analyze their structure and genetic diversity. These six domestic populations were clearly grouped into two independent clusters, with Jinchuan, Changtai, and Jiulong showing a tight genetic relationship with the wild yak. Nerve development pathways were enriched with GO enrichment analysis of 334 domesticated genes. Major genomic regions associated with the differentiation of domestic yaks were detected. These findings provide preliminary information on the yak genome variability, useful to understand the genomic characteristics of different populations in QTP.

Comparative analysis of gut fungal composition and structure of the yaks under different feeding models.

The yaks that inhabit the Tibetan plateau are a rare breed that is closely related to local economic development and human civilization. This ancient breed may have evolved a unique gut microbiota due to the hypoxic high-altitude environment. The gut microbiota is susceptible to external factors, but research regarding the effects of different feeding models on the gut fungal community in yaks remains scarce. In this study, we compared and analyzed the composition and variability of the gut fungal community among wild yaks (WYG), house-feeding domestic yaks (HFG), and grazing domestic yaks (GYG). The results revealed that Basidiomycota and Ascomycota were the most preponderant phyla in the gut fungal community, regardless of feeding models. Although the types of dominant fungal phyla did not change, their abundances did. Intergroup analysis of fungal diversity showed that the Shannon and Simpson indices of WYG and GYG were significantly higher than those of HFG. Fungal taxonomic analysis showed that there were 20 genera (Sclerostagonospora and Didymella) that were significantly different between WYG and GYG, and 16 genera (Thelebolus and Cystobasidium) that were significantly different between the WYG and HFG. Furthermore, the proportions of 14 genera (Claussenomyces and Papiliotrema) significantly decreased, whereas the proportions of eight genera (Stropharia and Lichtheimia) significantly increased in HFG as compared to GYG. Taken together, this study indicated that the gut fungal composition and structure differ significantly between yaks raised in different breeding groups.

Evolution of immunogenetic components encoding ultralong CDR H3.

The genomes of most vertebrates contain many V, D, and J gene segments within their Ig loci to construct highly variable CDR3 sequences through combinatorial diversity. This nucleotide variability translates into an antibody population containing extensive paratope diversity. Cattle have relatively few functional VDJ gene segments, requiring innovative approaches for generating diversity like the use of ultralong-encoding IGHV and IGHD gene segments that yield dramatically elongated CDR H3. Unique knob and stalk microdomains create protracted paratopes, where the antigen-binding knob sits atop a long stalk, allowing the antibody to bind both surface and recessed antigen epitopes. We examined genomes of twelve species of Bovidae to determine when ultralong-encoding IGHV and IGHD gene segments evolved. We located the 8-bp duplication encoding the unique TTVHQ motif in ultralong IGHV segments in six Bovid species (cattle, zebu, wild yak, domestic yak, American bison, and domestic gayal), but we did not find evidence of the duplication in species beyond the Bos and Bison genera. Additionally, we analyzed mRNA from bison spleen and identified a rich repertoire of expressed ultralong CDR H3 antibody mRNA, suggesting that bison use ultralong IGHV transcripts in their host defense. We found ultralong-encoding IGHD gene segments in all the same species except domestic yak, but again not beyond the Bos and Bison clade. Thus, the duplication event leading to this ultralong-encoding IGHV gene segment and the emergence of the ultralong-encoding IGHD gene segment appears to have evolved in a common ancestor of the Bos and Bison genera 5-10 million years ago.

"The Yak"-A remarkable animal living in a harsh environment: An overview of its feeding, growth, production performance, and contribution to food security.

Yaks play an important role in the livelihood of the people of the Qinghai-Tibet Plateau (QTP) and contribute significantly to the economy of the different countries in the region. Yaks are commonly raised at high altitudes of ~ 3,000-5,400 m above sea level. They provide many important products, namely, milk, meat, fur, and manure, as well as social status, etc. Yaks were domesticated from wild yaks and are present in the remote mountains of the QTP region. In the summer season, when a higher quantity of pasture is available in the mountain region, yaks use their long tongues to graze the pasture and spend ~ 30-80% of their daytime grazing. The remaining time is spent walking, resting, and doing other activities. In the winter season, due to heavy snowfall in the mountains, pasture is scarce, and yaks face feeding issues due to pasture scarcity. Hence, the normal body weight of yaks is affected and growth retardation occurs, which consequently affects their production performance. In this review article, we have discussed the domestication of yaks, the feeding pattern of yaks, the difference between the normal and growth-retarded yaks, and also their microbial community and their influences. In addition, blood biochemistry, the compositions of the yaks' milk and meat, and reproduction are reported herein. Evidence suggested that yaks play an important role in the daily life of the people living on the QTP, who consume milk, meat, fur, use manure for fuel and land fertilizer purposes, and use the animals for transportation. Yaks' close association with the people's well-being and livelihood has been significant.

Climatic factors and human disturbance influence ungulate species distribution on the Qinghai-Tibet Plateau

Due to human activities and climate change, the habitats of ungulate species on the Qinghai-Tibet Plateau have been greatly affected in recent decades. In this study, the distribution patterns of 19 ungulate species on the Qinghai-Tibet Plateau were identified based on MaxEnt model in the past (1960–1990) and current periods (2000–2015). Then the changes of their habitat distribution and the species richness in different periods were compared. Finally, the Zonation model was used to identify the key protected areas of ungulate species. The results show that the MaxEnt model can well predict the distribution of ungulate species. Most ungulate species are mainly distributed in the southeast of the Qinghai-Tibet Plateau. The distance to lakes and precipitation are the main factors affecting the distribution of most ungulate species. The habitats originally located in the southeast of the Qinghai-Tibet Plateau have mainly extended northeastwards, while the habitats originally located in the northwest has been largely lost. The changes in the habitats of ungulate species in the southeast and northwest are diametrically opposite. The key protected areas identified by Zonation model are mainly located in the southeast of the Qinghai-Tibet Plateau. The existing nature reserves can effectively protect the suitable habitats of the Tibetan antelope, Tibetan wild ass and wild yak. This research can provide scientific basis for coordinating the contradiction between development and protection and promoting the biodiversity conservation on the Qinghai-Tibet Plateau.

The Biological Characteristics and Differential Expression Patterns of TSSK1B Gene in Yak and Its Infertile Hybrid Offspring.

This study aimed to investigate the spatially and temporally expressed patterns and biological characteristics of TSSK1B in male yaks and explore the potential correlation between TSSK1B and male sterility of the yak hybrid offspring (termed cattle-yak). First, the coding sequence (CDS) of TSSK1B was cloned by RT-PCR, and bioinformatics analysis was conducted with relevant software. Quantitative real-time PCR (RT-qPCR) was employed to detect the expression profile of TSSK1B in various tissues of male adult yaks, the spatiotemporal expression of TSSK1B in different stages of yak testes, and the differential expression of TSSK1B between yak and cattle-yak testes. The cellular localization of TSSK1B was determined by immunohistochemistry (IHC). Furthermore, the methylation status of the TSSK1B promoter region was analyzed by bisulfite-sequencing PCR (BSP). The results showed that TSSK1B was 1235 bp long, including 1104 bp of the CDS region, which encoded 367 amino acids. It was a conserved gene sharing the highest homology with Bos mutus (99.67%). In addition, the bioinformatics analysis revealed that TSSK1B was an unstable hydrophilic protein mainly containing the alpha helix of 34.06% and a random coil of 44.41%, with a transmembrane structure of 29 amino acids long. The RT-qPCR results demonstrated that TSSK1B was specifically expressed in yak testes compared with that in other tissues and especially highly expressed in adult yak testes. On the contrary, TSSK1B was hardly expressed in the testis of adult cattle-yak. IHC confirmed that TSSK1B protein was more strongly expressed in the testes of adult yaks than in their fetal and juvenile counterparts. Interestingly, nearly no expression was observed in the testes of cattle-yak compared with the corresponding testes of yak. Bisulfite-sequencing PCR (BSP) revealed that the methylated CpG sites in the TSSK1B promoter region of cattle-yak was significantly higher than that in the yak. Taken together, this study revealed that TSSK1B was specifically expressed in yak testes and highly expressed upon sexual maturity. Moreover, the rare expression in cattle-yak may be related to the hypermethylation of the promoter region, thereby providing a basis for further studies on the regulatory mechanism of TSSK1B in male cattle-yak sterility.

Whole-genome resequencing reveals genetic diversity, differentiation, and selection signatures of yak breeds/populations in Qinghai, China.

The Qinghai Province of China is located in the northeast region of the Qinghai-Tibetan Plateau (QTP) and carries abundant yak genetic resources. Previous investigations of archaeological records, mitochondrial DNA, and Y chromosomal markers have suggested that Qinghai was the major center of yak domestication. In the present study, we examined the genomic diversity, differentiation, and selection signatures of 113 Qinghai yak, including 42 newly sequenced Qinghai yak and 71 publicly available individuals, from nine yak breeds/populations (wild, Datong, Huanhu, Xueduo, Yushu, Qilian, Geermu, Tongde, and Huzhu white) using high-depth whole-genome resequencing data. We observed that most of Qinghai yak breeds/populations have abundant genomic diversity based on four genomic parameters (nucleotide diversity, inbreeding coefficients, linkage disequilibrium decay, and runs of homozygosity). Population genetic structure analysis showed that Qinghai yak have two lineages with two ancestral origins and that nine yak breeds/populations are clustered into three distinct groups of wild yak, Geermu yak, and seven other domestic yak breeds/populations. F ST values showed moderate genetic differentiation between wild yak, Geermu yak, and the other Qinghai yak breeds/populations. Positive selection signals were detected in candidate genes associated with disease resistance (CDK2AP2, PLEC, and CYB5B), heat stress (NFAT5, HSF1, and SLC25A48), pigmentation (MCAM, RNF26, and BOP1), vision (C1QTNF5, MFRP, and TAX1BP3), milk quality (OPLAH and GRINA), neurodevelopment (SUSD4, INSYN1, and PPP1CA), and meat quality (ZRANB1), using the integrated PI, composite likelihood ratio (CLR), and F ST methods. These findings offer new insights into the genetic mechanisms underlying target traits in yak and provide important information for understanding the genomic characteristics of yak breeds/populations in Qinghai.