Y-Chromosome SNPs Reveal Paternal Lineages and Candidate Circadian-Regulatory Variation in Indian YaksY-Chromosome SNPs Reveal Paternal Lineages and Candidate Circadian-Regulatory Variation in Indian Yaks
63
Authors
-
MAHESH SHIVANAND DIGE
National Bureau of Animal Genetic Resources
Abstract
Yak (Bos grunniens), an economically and culturally important bovine species, is widely regarded as the lifeline of the Himalayan region. However, Indian yak populations remain comparatively understudied, with limited genomic resources to support conservation and genetic improvement. The mammalian Y chromosome, owing to its paternal inheritance, reduced recombination, and distinct evolutionary dynamics, provides a powerful framework for resolving male lineages and male-linked adaptation. In this study, we profiled Y-chromosome–specific variation in four yak populations—Arunachali, Himachali, Ladakhi, and Chinese Jinchuan—using whole-genome resequencing. We identified 274,828, 243,143, 283,774, and 194,228 Y-linked SNPs, respectively. Across populations, SNPs were most abundant in intergenic regions, followed by intronic regions; Jinchuan yak showed fewer SNPs, consistent with reduced diversity expected under intensive selective breeding. The transition/transversion (Ts/Tv) ratios were 1.57, 1.67, 1.55, and 1.99 for Arunachali, Himachali, Ladakhi, and Jinchuan yaks, respectively. Functional annotation using SnpEff indicated that most variants in Indian yaks occur in non-coding regions, while gene-based analyses highlighted key loci including ASMT, ASMTL, and SRY. Notably, ASMT and ASMTL participate in the melatonin biosynthesis pathway, implicating Y-linked (or PAR-adjacent) variation as candidate contributors to circadian regulation, a trait potentially relevant to adaptation under extreme high-altitude photoperiods. Finally, LD pruning retained only ~20% of SNPs, demonstrating that standard pruning parameters can substantially reshape Y-chromosome datasets and should be selected carefully to avoid loss of informative loci. Overall, this first Y-chromosome SNP resource for Indian yaks establishes a foundation for paternal haplotyping and population-genetic studies and motivates functional investigations into Y-linked contributions to chronobiological and high-altitude adaptation.
Downloads
References
Abramovs, N., Brass, A., & Tassabehji, M. (2020). Hardy-Weinberg equilibrium in the large-scale genomic sequencing era. Frontiers in Genetics, 11, 210. https://doi.org/10.3389/fgene.2020.00210
Andrews, S. (2010). FastQC: a quality control tool for high-throughput sequence data.
Bachtrog, D. (2013). Y-chromosome evolution: emerging insights into processes of Y-chromosome degeneration. Nature Reviews Genetics, 14(2), 113–124. https://doi.org/10.1038/nrg3366.
Bolger, A.M., Lohse, M., & Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 30(15), 2114–2120.
Brandli, L., Handley, L.J., Vogel, P., & Perrin, N. (2005). Evolutionary history of the greater white-toothed shrew (Crocidura russula) inferred from analysis of mtDNA, Y, and X chromosome markers. Molecular Phylogenetics and Evolution, 37, 832–844.
Chang, T.C., Yang, Y., Retzel, E.F., & Liu, W.S. (2013). Male-specific region of the bovine Y chromosome is gene-rich with high transcriptomic activity in testis development. Proceedings of the National Academy of Sciences, 110(30), 12373-12378.
Cingolani, P., Platts, A., Wang, L.L., Coon, M., Nguyen, T., Wang, L., Land, S.J., Lu, X., Ruden, D.M. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly, 6(2), 80–92. https://doi.org/10.4161/fly.19695.
Cortez, D., Marin, R., Toledo-Flores, D., Froidevaux, L., Liechti, A., Waters, P.D., Grützner, F., & Kaessmann, H. (2014). Origins and functional evolution of Y chromosomes across mammals. Nature, 508(7497), 488–493. https://doi.org/10.1038/nature13151.
Danecek, P., Bonfield, J.K., Liddle, J., Marshall, J., Ohan, V., Pollard, M.O., Whitwham, A., Keane, T., McCarthy, S.A., Davies, R.M., & Li, H. (2021). Twelve years of SAMtools and BCFtools. Gigascience, 10(2), giab008. https://doi.org/10.1093/gigascience/giab008.
Edwards, C.J., Gaillard, C., Bradley, D.G., & MacHugh, D.E. (2000). Y-specific microsatellite polymorphisms in a range of bovid species. Animal Genetics, 31, 127–130.
Goli, R. C., Mahar, K., Chishi, K. G., Koloi, S., Gunturu, T., Diwakar, V., ... & Kanaka, K. K. (2025). X-linked signatures of positive selection and adaptive genetic plasticity shaping high-altitude resilience in Indian Changthangi sheep. Small Ruminant Research, 247, 107492.
Guo, A. (1983). A comparative study on chromosomes of oxen (Bos taurus), yak (Bos grunniens), and their hybrids. Journal of Genetics and Genomics, 10, 137–143.
Hallast, P., Ebert, P., Loftus, M., Yilmaz, F., Audano, P.A., Logsdon, G.A., Bonder, M.J., Zhou, W., Höps, W., Kim, K., Li, C., Hoyt, S.J., Dishuck, P.C., Porubsky, D., Tsetsos, F., Kwon, J.Y., Zhu, Q., Munson, K.M., Hasenfeld, P., Harvey, W.T., Lewis, A.P., Kordosky, J., Hoekzema, K.; Human Genome Structural Variation Consortium (HGSVC); O'Neill, R.J., Korbel, J.O., Tyler-Smith, C., Eichler, E.E., Shi, X., Beck, C.R., Marschall, T., Konkel, M.K., Lee, C. (2023). Assembly of 43 human Y chromosomes reveals extensive complexity and variation. Nature, 621, 355-364. https://doi.org/10.1038/s41586-023-06425-6.
Han, J.L., Ochieng, J.W., Rege, J.E.O., & Hanotte, O. (2000). Low level of cattle introgression in yak populations from Bhutan and China: evidence from Y-specific microsatellites and mitochondrial DNA markers. In Proceedings of the Third International Congress on Yak (pp. 190–196). International Livestock Research Institute, Nairobi.
Hellborg, L., & Ellegren, H. (2004). Low levels of nucleotide diversity in mammalian Y chromosomes. Molecular Biology and Evolution, 21, 158–163.
Hughes, J.F., Skaletsky, H., Pyntikova, T., Minx, P.J., Graves, T., Rozen, S., Wilson, R.K., & Page, D.C. (2005). Conservation of Y-linked genes during human evolution revealed by comparative sequencing in chimpanzee. Nature, 437, 100–103.
Hurst, L.D., Atlan, A., & Bengtsson, B.O. (1996). Genetic conflicts. Quarterly Review of Biology, 71, 317–364.
Kanaka, K. K., Prabhu, I. G., Goli, R. C., Sruthi, S., Rathi, P., Mahar, K., ... & Sukhija, N. (2025). Hallmarks of selection located in scaffold-level polar bear genome sheds light on genomic rhythmicity. Biological Rhythm Research, 1-14.
Kanaka, K.K., Sukhija, N., Goli, R.C., Singh, S., Ganguly, I., Dixit, S.P., Dash, A., & Malik, A.A. (2023). On the concepts and measures of diversity in the genomics era. Current Plant Biology, 33, 100278. https://doi.org/10.1016/j.cpb.2023.100278.
Kumar, A., Dige, M., Niranjan, S.K., Ahlawat, S., Arora, R., Kour, A., & Vijh, R.K. (2024). Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations. Animal Biotechnology, 35(1), 2282723. https://doi.org/10.1080/10495398.2023.2282723.
Lachance, J., & Tishkoff, S.A. (2013). SNP ascertainment bias in population genetic analyses: why it is important, and how to correct it. Bioessays, 35(9), 780–786.
Lahn, B.T., & Page, D.C. (1999). Four evolutionary strata on the human X chromosome. Science, 286(5441), 964–967. https://doi.org/10.1126/science.286.5441.964.
Laird, N. M., & Lange, C. (2011). The Fundamentals of Modern Statistical Genetics. Springer Science & Business Media.
Lan, D., Xiong, X., Mipam, T.D., Fu, C., Li, Q., Ai, Y., Hou, D., Chai, Z., Zhong, J., & Li, J. (2018). Genetic diversity, molecular phylogeny, and selection evidence of Jinchuan yak revealed by whole-genome resequencing. G3: Genes, Genomes, Genetics, 8, 945–952. https://doi.org/10.1534/g3.118.300572.
Lawson, L.J., & Hewitt, G.M. (2002). Comparison of substitution rates in ZFX and ZFY introns of sheep and goat-related species supports the hypothesis of male-biased mutation rates. Journal of Molecular Evolution, 54, 54–61.
Li, H., & Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 25(14), 1754–1760.
Li, R., Chen, S., Li, C., Xiao, H., Costa, V., Bhuiyan, M. S. A., Baig, M., & Beja-Pereira, A. (2022). Whole-genome analysis deciphers population structure and genetic introgression among bovine species. Frontiers in Genetics, 13, 847492. https://doi.org/10.3389/fgene.2022.847492
Li, R., Wang, S.Q., Xu, S.Y., Huang, J.P., Wang, F.Q., Ma, Z.J., Dang, R.H., Lan, X.Y., Chen, H., & Lei, C.Z. (2014). Novel Y-chromosome polymorphisms in Chinese domestic yak. Animal Genetics, 45, 449–452.
Li, Y., Wu, D. D., Boyko, A. R., Wang, G. D., Wu, S. F., Irwin, D. M., & Zhang, Y. P. (2014). Population variation reveals high-altitude adaptation of Tibetan mastiffs. Molecular Biology and Evolution, 31(5), 1200–1205.
Li, Z., Löytynoja, A., Fraimout, A., & Merilä, J. (2019). Effects of marker type and filtering criteria on QST-FST comparisons. Royal Society Open Science, 6, 190666. https://doi.org/10.1098/rsos.190666
Liu, W., Huang, Z., Zhang, Y., Zhang, S., Cui, Z., Liu, W., Li, L., Xia, J., Zou, Y., & Qi, Z. (2023). ASMT determines gut microbiota and increases neurobehavioral adaptability to exercise in female mice. Communications Biology, 6(1), 1126. https://doi.org/10.1038/s42003-023-05520-8.
Luo, J., Wei, X., Liu, W., Chen, S., Ahmed, Z., Sun, W., Lei, C., & Ma, Z. (2022). Paternal genetic diversity, differentiation, and phylogeny of three white yak breeds/populations in China. Scientific Reports, 12(1), 19331. https://doi.org/10.1038/s41598-022-23453-w.
Ma, T., Tao, J., Yang, M., He, C., Tian, X., Zhang, X., Zhang, J., Deng, S., Feng, J., Zhang, Z., Wang, J., Ji, P., Song, Y., He, P., Han, H., Fu, J., Lian, Z., & Liu, G. (2017). An AANAT/ASMT transgenic animal model constructed with CRISPR/Cas9 system serving as the mammary gland bioreactor to produce melatonin-enriched milk in sheep. Journal of Pineal Research, 63(1). https://doi.org/10.1111/jpi.12406.
Ma, Z.J., Chen, S.M., Sun, Y.G., Xi, Y.L., Li, R.Z., Xu, J.T., & Lei, C.Z. (2015). Y-STR INRA189 polymorphisms in Chinese yak breeds. Genetics and Molecular Research, 14, 18859–18862.
Ma, Z.J., Xia, X.T., Chen, S.M., Zhao, X.C., Zeng, L.L., Xie, Y.L., Chao, S.Y., Xu, J.T., Sun, Y.G., Li, R.Z., Guanque, Z.X., Han, J.L., & Lei, C.Z. (2018). Identification and diversity of Y-chromosome haplotypes in Qinghai yak populations. Animal Genetics, 49(6), 618-622. https://doi.org/10.1111/age.12723.
Mahar, K., Goli, R. C., Chishi, K. G., Ganguly, I., Dixit, S. P., Singh, S., ... & Kareningappa, K. K. (2024). Runs of homozygosity Decipher genetic diversity in cattle breed dwelling in the colder regions of the world. Cytogenetic and Genome Research, 164(3-4), 154-164.
Mahar, K., Gurao, A., Kumar, A., et al. (2025). Genomic insights into high-altitude adaptation and evolutionary dynamics of Indian yaks in the Trans-Himalayan region. Conservation Genetics, 26, 49–62. https://doi.org/10.1007/s10592-024-01650-6.
Mahar, K., Gurao, A., Kumar, A., Pratap Singh, L., Chitkara, M., Gowane, G.R., Ahlawat, S., Niranjan, S.K., Pundir, R.K., Kataria, R.S., & Dige, M.S. (2024). Genomic inbreeding analysis reveals resilience and genetic diversity in Indian yak populations. Gene, 928, 148787. https://doi.org/10.1016/j.gene.2024.148787.
Mipam, T. D., Wen, Y., Fu, C., Li, S., Zhao, H., Ai, Y., Li, L., Zhang, L., & Zou, D. 2012. Maternal phylogeny of a newly-found yak population in China. International Journal of Molecular Sciences 13(9), 11455-11470. https://doi.org/10.3390/ijms130911455
Nei, M., & Li, W.-H. (1973). Linkage disequilibrium in subdivided populations. Genetics, 75(3), 213–219.
Nguyen, T.T., Genini, S., Menetrey, F., Malek, M., Voegeli, P., Goe, M.R., & Stranzinger, G. (2005). Application of bovine microsatellite markers for genetic diversity analysis of Swiss yak (Poephagus grunniens). Animal Genetics, 36, 484–489.
Popescu, C.P. (1969). Idiograms of yak (Bos grunniens), cattle (Bos taurus), and their hybrid. Genetics Selection Evolution, 1, 207. https://doi.org/10.1186/1297-9686-1-3-207.
Qi, X.B., Han, J.L., Rege, J.E.O., & Hanotte, O. (2002). Y-chromosome-specific microsatellite polymorphisms in Chinese yak. In Proceedings of the 7th World Congress on Genetics Applied to Livestock Production (pp. 509–512). Institut National de la Recherche Agronomique, Montpellier, France.
Qiu, Q., Zhang, G., Ma, T., Qian, W., Wang, J., Ye, Z., Cao, C., Hu, Q., Kim, J., Larkin, D.M., Auvil, L., Capitanu, B., Ma, J., Lewin, H.A., Qian, X., Lang, Y., Zhou, R., Wang, L., Wang, K., Xia, J., Liao, S., Pan, S., Lu, X., Hou, H., Wang, Y., Zang, X., Yin, Y., Ma, H., Zhang, J., Wang, Z., Zhang, Y., Zhang, D., Yonezawa, T., Hasegawa, M., Zhong, Y., Liu, W., Zhang, Y., Huang, Z., Zhang, S., Long, R., Yang, H., Wang, J., Lenstra, J.A., Cooper, D.N., Wu, Y., Wang, J., Shi, P., Wang, J., & Liu, J. (2012). The yak genome and adaptation to life at high altitude. Nature Genetics, 44(8), 946-949. https://doi.org/10.1038/ng.2343.
Quintana-Murci, L., & Fellous, M. (2001). The human Y chromosome: the biological role of a “functional wasteland”. Journal of Biomedicine and Biotechnology, 1(1), 18.
Rashid, I., Baisvar, V.S., Singh, M., Kumar, R., Srivastava, P., Kushwaha, B., & Pathak, A.K. (2020). Isolation and characterization of hypoxia-inducible heme oxygenase 1 (HMOX1) gene in Labeo rohita. Genomics, 112(3), 2327-333. https://doi.org/10.1016/j.ygeno.2020.01.004.
Roldan, E.R., & Gomendio, M. (1999). The Y chromosome as a battleground for sexual selection. Trends in Ecology & Evolution, 14, 58–62.
Sarkar, M., & Prakash, B. S. (2005). Circadian variations in plasma concentrations of melatonin and prolactin during breeding and non-breeding seasons in yak (Poephagus grunniens L.). Animal Reproduction Science, 90(1-2), 149-162.
Sivalingam, J., Vineeth, M.R., Surya, T., Singh, K., Dixit, S.P., Niranjan, S.K., ... & Ravikumar, D. (2020). Genomic divergence reveals unique populations among Indian yaks. Scientific Reports, 10(1), 3636.
Slatkin, M. (2008). Linkage disequilibrium—understanding the evolutionary past and mapping the medical future. Nature Reviews Genetics, 9(6), 477–485.
Sukhija, N., Goli, R.C., Sukhija, M., et al. (2024). Evolutionary stamps for adaptation traced in Cervus nippon genome using reduced representation sequencing. Conservation Genetics Resources, 16, 135–146. https://doi.org/10.1007/s12686-023-01339-x.
Tian, Q., Tong, P., Chen, G., Deng, M., Cai, T., Tian, R., Zhang, Z., Xia, K., & Hu, Z. (2023). GLRA2 gene mutations cause high myopia in humans and mice. Journal of Medical Genetics, 60(2), 193-203. https://doi.org/10.1136/jmedgenet-2022-108425.
Visser, C., Lashmar, S. F., Van Marle-Köster, E., Poli, M. A., & Allain, D. (2016). Genetic diversity and population structure in South African, French, and Argentinian Angora goats from genome-wide SNP data. PLOS ONE, 11(5), e0154353. https://doi.org/10.1371/journal.pone.0154353
Wang, X., Pei, J., Bao, P., Cao, M., Guo, S., Song, R., Song, W., Liang, C., Yan, P., & Guo, X. (2021). Mitogenomic diversity and phylogeny analysis of yak (Bos grunniens). BMC Genomics, 22(1), 325. https://doi.org/10.1186/s12864-021-07650-x.
Wang, Y., Cai, H., Luo, X., Ai, Y., Jiang, M., & Wen, Y. (2020). Insight into unique somitogenesis of yak (Bos grunniens) with one additional thoracic vertebra. BMC Genomics, 21, 1–15.
Wu, H., Cui, X., Guan, S., Li, G., Yao, Y., Wu, H., Zhang, J., Zhang, X., Yu, T., Li, Y., Lian, Z., Zhang, L., & Liu, G. (2022). The improved milk quality and enhanced anti-inflammatory effect in acetylserotonin-O-methyltransferase (ASMT) overexpressed goats: An association with the elevated endogenous melatonin production. Molecules, 27(2), 572. https://doi.org/10.3390/molecules27020572.
Yao, S., Liu, Y., Liu, X., & Liu, G. (2022). Effects of SNPs in AANAT and ASMT genes on milk and peripheral blood melatonin concentrations in Holstein cows (Bos taurus). Genes, 13(7), 1196.
Zhang, G.W., Guan, J.Q., Luo, Z.G., Zhang, W.X., Wang, L., Luo, X.L., & Zuo, F.Y. (2016). A tremendous expansion of TSPY copy number in crossbred bulls (Bos taurus × Bos grunniens). Journal of Animal Science, 94(4), 1398-1407.
Zhang, R., Wu, Y., Zhao, M., Liu, C., Zhou, L., Shen, S., Liao, S., Yang, K., Li, Q., & Wan, H. (2009). Role of HIF-1α in the regulation of ACE and ACE2 expression in hypoxic human pulmonary artery smooth muscle cells. American Journal of Physiology - Lung Cellular and Molecular Physiology, 297(4), L631-640. https://doi.org/10.1152/ajplung.90415.2008.
Published
Issue
Section
License
Copyright (c) 2026 The Indian Journal of Animal Sciences
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.
The copyright of the articles published in The Indian Journal of Animal Sciences is vested with the Indian Council of Agricultural Research, which reserves the right to enter into any agreement with any organization in India or abroad, for reprography, photocopying, storage and dissemination of information. The Council has no objection to using the material, provided the information is not being utilized for commercial purposes and wherever the information is being used, proper credit is given to ICAR.
How to Cite
