Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia

The Genome Landscape of Tibetan Sheep Reveals Adaptive Introgression from Argali and the History of Early Human Settlements on the Qinghai-Tibetan Plateau

Tibetan sheep are the most common and widespread domesticated animals on the Qinghai-Tibetan Plateau (QTP) and have played an essential role in the permanent human occupation of this high-altitude region. However, the precise timing, route, and process of sheep pastoralism in the QTP region remain poorly established, and little is known about the underlying genomic changes that occurred during the process. Here, we investigate the genomic variation in Tibetan sheep using whole-genome sequences, single nucleotide polymorphism arrays, mitochondrial DNA, and Y-chromosomal variants in 986 samples throughout their distribution range. We detect strong signatures of selection in genes involved in the hypoxia and ultraviolet signaling pathways (e.g., HIF-1 pathway and HBB and MITF genes) and in genes associated with morphological traits such as horn size and shape (e.g., RXFP2). We identify clear signals of argali (Ovis ammon) introgression into sympatric Tibetan sheep, covering 5.23-5.79% of their genomes. The introgressed genomic regions are enriched in genes related to oxygen transportation system, sensory perception, and morphological phenotypes, in particular the genes HBB and RXFP2 with strong signs of adaptive introgression. The spatial distribution of genomic diversity and demographic reconstruction of the history of Tibetan sheep show a stepwise pattern of colonization with their initial spread onto the QTP from its northeastern part ∼3,100 years ago, followed by further southwest expansion to the central QTP ∼1,300 years ago. Together with archeological evidence, the date and route reveal the history of human expansions on the QTP by the Tang-Bo Ancient Road during the late Holocene. Our findings contribute to a depth understanding of early pastoralism and the local adaptation of Tibetan sheep as well as the late-Holocene human occupation of the QTP.

A Novel SNP in EIF2AK4 Gene Is Associated with Thermal Tolerance Traits in Chinese Cattle

Simple Summary

China harbors two lineages of cattle (Bos taurus and Bos indicus) that display pronounced geographical distribution differences. Northern Chinese cattle predominantly belong to B. taurus and southern Chinese cattle belong to B. indicus. Both B. taurus and B. indicus contribute to the admixture of cattle in central China. Thermal stress induces oxidative stress and DNA damage in mammals. In general, B. indicus are more resistant to thermal stress than B. taurus. Eukaryotic translation initiation factor 2-alpha kinase 4 (EIF2AK4), which pertains to the family of serine–threonine kinase, is a candidate gene for thermal stress. However, the effects of the bovine EIF2AK4 gene on the thermal tolerance traits of Chinese cattle breeds remain unknown. Our results suggest that a variant of the EIF2AK4 gene is associated with thermal tolerance traits in Chinese cattle.

Abstract

Eukaryotic translation initiation factor 2-alpha kinase 4 (EIF2AK4, also known as GCN2), which pertains to the family of serine–threonine kinase, is involved in oxidative stress and DNA damage repair. A missense single-nucleotide polymorphism (SNP) (NC_037337.1 g.35615224 T > G) in exon 6 of the EIF2AK4 gene which encodes a p.Ile205Ser substitution was observed in the Bovine Genome Variation Database and Selective Signatures (BGVD). The purpose of the current study is to determine the allelic frequency distribution of the locus and analyze its association with thermal tolerance in Chinese indigenous cattle. In our study, the allelic frequency distribution of the missense mutation (NC_037337.1 g.35615224 T > G) in Chinese cattle was analyzed by sequencing 1105 individuals of 37 breeds including 35 Chinese indigenous cattle breeds and two exotic breeds. In particular, association analysis was carried out between the genotypes and three environmental parameters including annual mean temperature (T), relative humidity (RH), and temperature–humidity index (THI). The frequency of the mutant allele G (NC_037337.1 g.35615224 T > G) gradually decreased from the southern cattle groups to the northern cattle groups, whereas the frequency of the wild-type allele T showed an opposite pattern, consistent with the distribution of indicine and taurine cattle in China. In accordance with the association analysis, genotypes were significantly associated with T (P < 0.01), RH (P < 0.01), and THI (P < 0.01), suggesting that the cattle with genotype GG were found in regions with higher T, RH, and THI. Thus, our results suggest that the mutation (NC_037337.1 g.35615224 T > G) of the EIF2AK4 gene is associated with thermal tolerance traits in Chinese cattle.

Exploring genetic diversity and phylogenic relationships of Chinese cattle using gene mtDNA 16S rRNA

The objective of this research was to characterize the genetic diversity and phylogenetic diversity among 12 cattle breeds (10 Chinese breeds and two foreign taurine breeds as controls) utilizing gene mtDNA 16S rRNA. The complete sequences of the mtDNA 16S rRNA genes of the 251 animals were 1570 bp long. The mean percentages of the four nitrogen bases were 37.8 % for adenine (A), 23.7 % for thymine (T), 20.9 % for cytosine (C), and 17.6 % for guanine (G). The mtDNA 16S rRNA gene base percentages had a strong bias towards A  +  T. All detected nucleotide variations in gene mtDNA 16S rRNA were either transitions (62.3 %) or transversions (37.7 %); no indels (insertions and deletions) were found. A total of 40 haplotypes were constructed based on these mutations. A total of 36 haplotypes of these 40 haplotypes were present in 10 Chinese cattle breeds. The haplotype diversity of all Chinese cattle populations was 0.903 ± 0.077 , while the nucleotide diversity was 0.0071 ± 0.0039 . Kimura's two-parameter genetic distances between pairs of the studied 12 breeds ranged from 0.001 to 0.010. The phylogenetic analysis assigned the 10 Chinese breeds to two distinct lineages that likely differed in their percentage of Bos taurus and Bos indicus ancestry.

A First Y-Chromosomal Haplotype Network to Investigate Male-Driven Population Dynamics in Domestic and Wild Bactrian Camels

Polymorphic markers on the male-specific part of the Y chromosome (MSY) provide useful information for tracking male genealogies. While maternal lineages are well studied in Old World camelids using mitochondrial DNA, the lack of a Y-chromosomal reference sequence hampers the analysis of male-driven demographics. Recently, a shotgun assembly of the horse MSY was generated based on short read next generation sequencing data. The haplotype network resulting from single copy MSY variants using the assembly as a reference revealed sufficient resolution to trace individual male lines in this species. In a similar approach we generated a 3.8 Mbp sized assembly of the MSY of Camelus bactrianus. The camel MSY assembly was used as a reference for variant calling using short read data from eight Old World camelid individuals. Based on 596 single nucleotide variants we revealed a Y-phylogenetic network with seven haplotypes. Wild and domestic Bactrian camels were clearly separated into two different haplogroups with an estimated divergence time of 26,999 ± 2,268 years. Unexpectedly, one wild camel clustered into the domestic Bactrian camels' haplogroup. The observation of a domestic paternal lineage within the wild camel population is concerning in view of the importance to conserve the genetic integrity of these highly endangered species in their natural habitat.

Genomic variants identified from whole-genome resequencing of indicine cattle breeds from Pakistan

The primary goal of cattle genomics is the identification of genome-wide polymorphism associated with economically important traits. The bovine genome sequencing project was completed in 2009. Since then, using massively parallel sequencing technologies, a large number of Bos taurus cattle breeds have been resequenced and scanned for genome-wide polymorphisms. As a result, a substantial number of single nucleotide polymorphisms (SNPs) have been discovered across European Bos taurus genomes, whereas extremely less number of SNPs are cataloged for Bos indicus breeds. In this study, we performed whole-genome resequencing, reference-based mapping, functional annotation and gene enrichment analysis of 20 sires representing eleven important Bos indicus (indicine) breeds of Pakistan. The breeds sequenced here include: Sahiwal, Red Sindhi, Tharparkar and Cholistani (tropically adapted dairy and dual purpose breeds), Achai, Bhagnari, Dajal and Lohani (high altitude adapted dual and drought purpose breeds); Dhanni, Hisar Haryana and Gabrali (dairy and light drought purpose breeds). A total of 17.4 billion QC passed reads were produced using BGISEQ-500 next generation sequencing platform to generate 9 to 27-fold genome coverage (average ~16×) for each of the 20 sequenced sires. A total of 67,303,469 SNPs were identified, of which 3,850,365 were found novel and 1,083,842 insertions-deletions (InDels) were detected across the whole sequenced genomes (491,247 novel). Comparative analysis using coding region SNPs revealed a close relationship between the best milking indicine breeds; Red Sindhi and Sahiwal. On the other hand, Bhagnari and Tharparkar being popular for their adaptation to dry and extremely hot climates were found to share the highest number of SNPs. Functional annotation identified a total of 3,194 high-impact (disruptive) SNPs and 745 disruptive InDels (in 275 genes) that may possibly affect economically important dairy and beef traits. Functional enrichment analysis was performed and revealed that high or moderate impact variants in wingless-related integration site (Wnt) and vascular smooth muscle contraction (VSMC) signaling pathways were significantly over-represented in tropically adapted heat tolerant Pakistani-indicine breeds. On the other hand, vascular endothelial growth factor (VEGF) and hypoxia-inducible factor 1 (HIF-1) signaling pathways were found over-represented in highland adapted Pakistani-indicine breeds. Similarly, the ECM-receptor interaction and Jak-STAT signaling pathway were significantly enriched in dairy and beef purpose Pakistani-indicine cattle breeds. The Toll-like receptor signaling pathway was significantly enriched in most of the Pakistani-indicine cattle. Therefore, this study provides baseline data for further research to investigate the molecular mechanisms of major traits and to develop potential genomic markers associated with economically important breeding traits, particularly in indicine cattle.

Deep Genome Resequencing Reveals Artificial and Natural Selection for Visual Deterioration, Plateau Adaptability and High Prolificacy in Chinese Domestic Sheep

Sheep were one of the earliest domesticated animals. Both artificial and natural selection during domestication has resulted in remarkable changes in behavioral, physiological, and morphological phenotypes; however, the genetic mechanisms underpinning these changes remain unclear, particularly for indigenous Chinese sheep. In the present study, we performed pooled whole-genome resequencing of 338 sheep from five breeds representative of indigenous Chinese breeds and compared them to the wild ancestors of domestic sheep (Asian mouflon, Ovis orientalis) for detection of genome-wide selective sweeps. Comparative genomic analysis between domestic sheep and Asian mouflon showed that selected regions were enriched for genes involved in bone morphogenesis, growth regulation, and embryonic and neural development in domestic sheep. Moreover, we identified several vision-associated genes with funtional mutations, such as PDE6B (c.G2994C/p.A982P and c.C2284A/p.L762M mutations), PANK2, and FOXC1/GMSD in all five Chinese native breeds. Breed-specific selected regions were determined including genes such as CYP17 for hypoxia adaptability in Tibetan sheep and DNAJB5 for heat tolerance in Duolang sheep. Our findings provide insights into the genetic mechanisms underlying important phenotypic changes that have occurred during sheep domestication and subsequent selection.

Genome-Wide SNP Discovery in Indigenous Cattle Breeds of South Africa

Single nucleotide polymorphism arrays have created new possibilities for performing genome-wide studies to detect genomic regions harboring sequence variants that affect complex traits. However, the majority of validated SNPs for which allele frequencies have been estimated are limited primarily to European breeds. The objective of this study was to perform SNP discovery in three South African indigenous breeds (Afrikaner, Drakensberger, and Nguni) using whole genome sequencing. DNA was extracted from blood and hair samples, quantified and prepared at 50 ng/μl concentration for sequencing at the Agricultural Research Council Biotechnology Platform using an Illumina HiSeq 2500. The fastq files were used to call the variants using the Genome Analysis Tool Kit. A total of 1,678,360 were identified as novel using Run 6 of 1000 Bull Genomes Project. Annotation of the identified variants classified them into functional categories. Within the coding regions, about 30% of the SNPs were non-synonymous substitutions that encode for alternate amino acids. The study of distribution of SNP across the genome identified regions showing notable differences in the densities of SNPs among the breeds and highlighted many regions of functional significance. Gene ontology terms identified genes such as MLANA, SYT10, and CDC42EP5 that have been associated with coat color in mouse, and ADAMS3, DNAJC3, and PAG5 genes have been associated with fertility in cattle. Further analysis of the variants detected 688 candidate selective sweeps (ZHp Z-scores ≤ -4) across all three breeds, of which 223 regions were assigned as being putative selective sweeps (ZHp scores ≤-5). We also identified 96 regions with extremely low ZHp Z-scores (≤-6) in Afrikaner and Nguni. Genes such as KIT and MITF that have been associated with skin pigmentation in cattle and CACNA1C, which has been associated with biopolar disorder in human, were identified in these regions. This study provides the first analysis of sequence data to discover SNPs in indigenous South African cattle breeds. The information will play an important role in our efforts to understand the genetic history of our cattle and in designing appropriate breed improvement programmes.

Evaluation of Linkage Disequilibrium, Effective Population Size and Haplotype Block Structure in Chinese Cattle

Simple Summary

Evaluation of the population structure and linkage disequilibrium can offer important insights to fully understand the genetic diversity and population history of cattle, which can enable us to appropriately design and implement GWAS and GS in cattle. In this study, we characterized the extent of genome-wide LD and the haplotype block structure, and estimated the persistence of phase of Chinese indigenous cattle with Illumina BovineHD BeadChip. According to our study, 58K, 87K, 95K, 52K, and 52K markers would be necessary for SCHC, NCC, SWC, SIM, and WAG, respectively, in the implementation of GWAS and GS and combining a multipopulation with high persistence of phase is feasible for the implication of genomic selection for Chinese beef cattle.

Abstract

Understanding the linkage disequilibrium (LD) across the genome, haplotype structure, and persistence of phase between breeds can enable us to appropriately design and implement the genome-wide association (GWAS) and genomic selection (GS) in beef cattle. We estimated the extent of genome-wide LD, haplotype block structure, and the persistence of phase in 10 Chinese cattle population using high density BovinHD BeadChip. The overall LD measured by r² between adjacent SNPs were 0.60, 0.67, 0.58, 0.73, and 0.71 for South Chinese cattle (SCHC), North Chinese cattle (NCC), Southwest Chinese cattle (SWC), Simmental (SIM), and Wagyu (WAG). The highest correlation (0.53) for persistence of phase across groups was observed for SCHC vs. SWC at distances of 0–50 kb, while the lowest correlation was 0.13 for SIM vs. SCHC at the same distances. In addition, the estimated current effective population sizes were 27, 14, 31, 34, and 43 for SCHC, NCC, SWC, SIM, and WAG, respectively. Our result showed that 58K, 87K, 95K, 52K, and 52K markers were required for implementation of GWAS and GS in SCHC, NCC, SWC, SIM, and WAG, respectively. Also, our findings suggested that the implication of genomic selection for multipopulation with high persistence of phase is feasible for Chinese cattle.

Whole-Genome Sequencing of Three Native Cattle Breeds Originating From the Northernmost Cattle Farming Regions

Northern Fennoscandia and the Sakha Republic in the Russian Federation represent the northernmost regions on Earth where cattle farming has been traditionally practiced. In this study, we performed whole-genome sequencing to genetically characterize three rare native breeds Eastern Finncattle, Western Finncattle and Yakutian cattle adapted to these northern Eurasian regions. We examined the demographic history, genetic diversity and unfolded loci under natural or artificial selection. On average, we achieved 13.01-fold genome coverage after mapping the sequencing reads on the bovine reference genome (UMD 3.1) and detected a total of 17.45 million single nucleotide polymorphisms (SNPs) and 1.95 million insertions-deletions (indels). We observed that the ancestral species (Bos primigenius) of Eurasian taurine cattle experienced two notable prehistorical declines in effective population size associated with dramatic climate changes. The modern Yakutian cattle exhibited a higher level of within-population variation in terms of number of SNPs and nucleotide diversity than the contemporary European taurine breeds. This result is in contrast to the results of marker-based cattle breed diversity studies, indicating assortment bias in previous analyses. Our results suggest that the effective population size of the ancestral Asiatic taurine cattle may have been higher than that of the European cattle. Alternatively, our findings could indicate the hybrid origins of the Yakutian cattle ancestries and possibly the lack of intensive artificial selection. We identified a number of genomic regions under selection that may have contributed to the adaptation to the northern and subarctic environments, including genes involved in disease resistance, sensory perception, cold adaptation and growth. By characterizing the native breeds, we were able to obtain new information on cattle genomes and on the value of the adapted breeds for the conservation of cattle genetic resources.

Genomic diversity and admixture patterns among six Chinese indigenous cattle breeds in Yunnan

Objective

Yunnan is not only a frontier zone that connects China with South and Southeast Asia, but also represents an admixture zone between taurine (Bos taurus) and zebu (Bos indicus) cattle. The purpose of this study is to understand the level of genomic diversity and the extent of admixture in each Yunnan native cattle breed.

Methods

All 120 individuals were genotyped using Illumina BovineHD BeadChip (777,962 single nucleotide polymorphisms [SNPs]). Quality control and genomic diversity indexes were calculated using PLINK software. The principal component analysis (PCA) was assessed using SMARTPCA program implemented in EIGENSOFT software. The ADMIXTURE software was used to reveal admixture patterns among breeds.

Results

A total of 604,630 SNPs was obtained after quality control procedures. Among six breeds, the highest level of mean heterozygosity was found in Zhaotong cattle from Northeastern Yunnan, whereas the lowest level of heterozygosity was detected in Dehong humped cattle from Western Yunnan. The PCA based on a pruned dataset of 233,788 SNPs clearly separated Dehong humped cattle (supposed to be a pure zebu breed) from other five breeds. The admixture analysis further revealed two clusters (K = 2 with the lowest cross validation error), corresponding to taurine and zebu cattle lineages. All six breeds except for Dehong humped cattle showed different degrees of admixture between taurine and zebu cattle. As expected, Dehong humped cattle showed no signature of taurine cattle influence.

Conclusion

Overall, considerable genomic diversity was found in six Yunnan native cattle breeds except for Dehong humped cattle from Western Yunnan. Dehong humped cattle is a pure zebu breed, while other five breeds had admixed origins with different extents of admixture between taurine and zebu cattle. Such admixture by crossbreeding between zebu and taurine cattle facilitated the spread of zebu cattle from tropical and subtropical regions to other highland regions in Yunnan.

The Transcriptomic Landscape of Yaks Reveals Molecular Pathways for High Altitude Adaptation

Yak is one of the largest native mammalian species at the Himalayas, the highest plateau area in the world with an average elevation of >4,000 m above the sea level. Yak is well adapted to high altitude environment with a set of physiological features for a more efficient blood flow for oxygen delivery under hypobaric hypoxia. Yet, the genetic mechanism underlying its adaptation remains elusive. We conducted a cross-tissue, cross-altitude, and cross-species study to characterize the transcriptomic landscape of domestic yaks. The generated multi-tissue transcriptomic data greatly improved the current yak genome annotation by identifying tens of thousands novel transcripts. We found that among the eight tested tissues (lung, heart, kidney, liver, spleen, muscle, testis, and brain), lung and heart are two key organs showing adaptive transcriptional changes and >90% of the cross-altitude differentially expressed genes in lung display a nonlinear regulation. Pathways related to cell survival and proliferation are enriched, including PI3K-Akt, HIF-1, focal adhesion, and ECM–receptor interaction. These findings, in combination with the comprehensive transcriptome data set, are valuable to understanding the genetic mechanism of hypoxic adaptation in yak.

Maternal genetic and phylogenetic characteristics of domesticated cattle in northwestern China

Northwestern China, an important part of the Silk Road, was the birthplace of the Zhouzu farming culture. The domestication of cattle as an important aspect of farming culture has had a long history in northwestern China. In this study, we assessed the maternal structure and phylogeny of cattle by analyzing the mitochondrial DNA hypervariable segment I (HVS-I) in 698 native cattle from eight areas of northwestern China. The phylogenetic analyses revealed two highly divergent mtDNA clades: clade T, which had four sub-clades (Ta—Td), and clade I. The cattle domesticated from Bos taurus showed a clear dominant distribution pattern in northwestern China. The nucleotide diversity of the Bos indicus clade was lower than that of clades from Bos taurus. In summary, our results suggest that the native cattle of northwestern China were domesticated from two different maternal ancestors, Bos taurus and Bos indicus, which migrated to the central plains of China from the north and south, respectively, with Bos taurus remaining at the edges of the region. The population expansion of the cattle domesticated from Bos taurus occurred in the Longdong region of Gansu Province, and these cattle formed four relatively independent evolutionary branches. Subsequent to this expansion event, Bos indicus migrated from southern to northern China.

Legacies of domestication, trade and herder mobility shape extant male zebu cattle diversity in South Asia and Africa

All tropically adapted humped cattle (Bos indicus or “zebu”), descend from a domestication process that took place >8,000 years ago in South Asia. Here we present an intercontinental survey of Y-chromosome diversity and a comprehensive reconstruction of male-lineage zebu cattle history and diversity patterns. Phylogenetic analysis revealed that all the zebu Y-chromosome haplotypes in our dataset group within three different lineages: Y3_A, the most predominant and cosmopolitan lineage; Y3_B, only observed in West Africa; and Y3_C, predominant in South and Northeast India. The divergence times estimated for these three Zebu-specific lineages predate domestication. Coalescent demographic models support either de novo domestication of genetically divergent paternal lineages or more complex process including gene flow between wild and domestic animals. Our data suggest export of varied zebu lineages from domestication centres through time. The almost exclusive presence of Y3_A haplotypes in East Africa is consistent with recent cattle restocking in this area. The cryptic presence of Y3_B haplotypes in West Africa, found nowhere else, suggests that these haplotypes might represent the oldest zebu lineage introduced to Africa ca. 3,000 B.P. and subsequently replaced in most of the world. The informative ability of Interspersed Multilocus Microsatellites and Y-specific microsatellites to identify genetic structuring in cattle populations is confirmed.

Comprehensive analysis of the mitochondrial DNA diversity in Chinese cattle

Complete mitochondrial DNA D-loop sequences of 1105 individuals were used to assess the diversity of maternal lineages of cattle populations in China. In total, 250 taurine and 88 zebu haplotypes were identified. Five main haplogroups-T1a, T2, T3, T4 and T5-were identified in Bos taurus, whereas Bos indicus harbored two haplogroups-I1 and I2. Our results suggest that the distribution of T1a in Asia was concentrated mainly in the northeast region (northeast China, Korea and Japan); haplogroups T2, T3 and T4 were predominant in Chinese cattle; and T5 was sporadically detected in Mongolian and Pingwu cattle. In contrast to the widespread presence of I1, I2 was distributed only in southwestern China (Yunnan-Guizhou Plateau and the Tibet Autonomous Region) and Xinjiang Uygur Autonomous Region. This is the first time that all five taurine haplogroups and two zebu haplogroups have been found in Mongolian cattle. In addition, eight individuals in Tibetan cattle carried the Bos grunniens mtDNA type. The high mtDNA diversity (H = 0.904 ± 0.008) and the weak genetic structure among the 57 Chinese cattle breeds/populations are consistent with their complex historical background, migration route and ecological environment.

Genetic diversity of Chinese cattle revealed by Y-SNP and Y-STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y-SNPs (UTY19 and ZFY10) and Y-STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y-chromosome haplotypes that represented three haplogroups. Y2-104-158 and Y2-102-158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3-88-156 predominates in southern China. Haplotypes Y2-108-158, Y2-110-158, Y2-112-158 and Y3-92-156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.

Identification and diversity of Y-chromosome haplotypes in Qinghai yak populations

The aim of the present study was to perform a preliminary analysis of the characterization and diversity of Y-chromosome haplotypes/haplogroups in yak of Qinghai Province, China. A total of 322 male yaks from nine populations belonging to three officially recognized breeds (Gaoyuan, Huanhu and Datong) were sampled. Animals were genotyped using six previously reported Y-SNPs present in the SRY, USP9Y, UTY, AMELY and OFD1Y genes and four new Y-SNPs in the OFD1Y gene (g.569A>C, g.578A>C, g.608G>T and g.653G>C) identified in this study. Seven Y-chromosome haplotypes (H1-H7) were identified according to the combination of the 10 Y-SNPs. H1, H2 and H6 were the most common and shared haplotypes across all yak populations/breeds. Private haplotypes H3 and H7 were detected in the Datong breed; H4 in Guoleimude, Qumalai, Qilian, Tianjun and Ganglong populations; and H5 in Qumalai of Gaoyuan breed. Haplotype clustering and network analyses inferred two haplogroups, Y1 and Y2, indicating two divergent lineages of paternal origins of Qinghai yak. The analysis of molecular variance showed a significant difference among individuals (P < 0.0001) with more than 93% of the total genetic variation present within populations, suggesting a weak genetic structure among Qinghai yak populations. The overall Y-haplotype diversity was 0.538 ± 0.028, showing a relatively high diversity in Qinghai yak. The Gaoyuan and Datong breeds had similar haplotype diversities (0.547 ± 0.030 and 0.553 ± 0.083, respectively), which were higher than that of the Huanhu breed (0.441 ± 0.098). Our results support the conservation and sustainable use of unique yak genetic resources in Qinghai.