Abstract
Sika deer is indigenous to East Asia. Currently the species is extinct or are nearing the verge of extinction, in numerous places. Reduced population size is sensitive to climate change and has poor adaptation as a result of inbreeding. In this study genome-wide Single nucleotide polymorphisms (SNPs) were identified using double-digest restriction site associated digestion (ddRAD) approach and annotated further to identify genes related to fitness. A total 85,001 SNPs were retained following the quality control. Composite Likelihood Ratio (CLR) statistics found 13 stronger selective sweeps, identified sweeps overlapped with genes related to adaptation, reproduction, cellular functions. Average heterozygosity and nucleotide diversity were 0.27 and 0.29 respectively. Admixture analysis furnished a single cluster indicating no further sub-divergence of Sika deer. These findings would help fabricate appropriate genetic improvement programs by inclusion in the SNP chips, in order to get deeper insights into deer genetics, evolution, domestication and conservation.
Graphical Abstract
Bioinformatic analysis for unravelling genomic variants and selective sweep regions in Cervus nippon
Similar content being viewed by others
References
Abramovs N, Brass A, Tassabehji M (2020) Hardy-Weinberg equilibrium in the large scale genomic sequencing era. Front Genet 11:210
Alberts, B., Johnson, A., Lewis, J., Raff, M., Roberts, K., & Walter, P. (2002). The shape and structure of proteins. In Molecular biology of the cell (4th ed). Garland Science
Alexander DH, Novembre J, Lange K (2009) Fast model-based estimation of ancestry in unrelated individuals. Genome Res 19(9):1655–1664. https://doi.org/10.1101/gr.094052.109
Andrews, S. (2010). FastQC: A quality control tool for high throughput sequence data. http://www.bioinformatics.babraham.ac.uk/projects/fastqc. Accessed 12 Jan 2023
Ba H, Jia B, Wang G, Yang Y, Kedem G, Li C (2017) Genome-wide SNP discovery and analysis of genetic diversity in farmed sika deer (Cervus nippon) in Northeast China using double-digest restriction site-associated DNA sequencing. G3 7(9):3169–3176. https://doi.org/10.1534/g3.117.300082
Bahbahani H, Clifford H, Wragg D, Mbole-Kariuki MN, Van Tassell C, Sonstegard T, Hanotte O (2015) Signatures of positive selection in East African Shorthorn Zebu: a genome-wide single nucleotide polymorphism analysis. Sci Rep 5(1):11729
Barbato M, Orozco-terWengel P, Tapio M, Bruford MW (2015) SNeP: a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data. Front Genet 6:109. https://doi.org/10.3389/fgene.2015.00109
Barczak E, Wolc A, Wójtowski J, Ślósarz P, Szwaczkowski T (2009) Inbreeding and inbreeding depression on body weight in sheep. J Anim Feed Sci 18(1):42–50. https://doi.org/10.22358/jafs/66366/2009
Ben-Jemaa S, Adam G, Boussaha M, Bardou P, Klopp C, Mandonnet N, Naves M (2023) Whole genome sequencing reveals signals of adaptive admixture in Creole cattle. Sci Rep 13(1):12155. https://doi.org/10.1038/s41598-023-38774-7
Brouard JS, Boyle B, Ibeagha-Awemu EM, Bissonnette N (2017) Low-depth genotyping-by-sequencing (GBS) in a bovine population: strategies to maximize the selection of high quality genotypes and the accuracy of imputation. BMC Genet 18(1):32. https://doi.org/10.1186/s12863-017-0501-y
Callus BA, Verhagen AM, Vaux DL (2006) Association of mammalian sterile twenty kinases, Mst1 and Mst2, with hSalvador via C-terminal coiled-coil domains, leads to its stabilization and phosphorylation. FEBS J 273(18):4264–4276. https://doi.org/10.1111/j.1742-4658.2006.05427.x
Catchen JM, Amores A, Hohenlohe P, Cresko W, Postlethwait JH (2011) Stacks: building and genotyping loci de novo from short-read sequences. G3 1(3):171–182. https://doi.org/10.1534/g3.111.000240
Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM (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
DaCosta JM, Sorenson MD (2016) ddRAD-seq phylogenetics based on nucleotide, indel, and presence–absence polymorphisms: analyses of two avian genera with contrasting histories. Mol Phylogenetics Evol 94:122–135. https://doi.org/10.1016/j.ympev.2015.07.026
Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, McVean G, Durbin R, & 1000 Genomes Project Analysis Group (2011) The variant call format and VCFtools. Bioinformatics 27(15):2156–2158. https://doi.org/10.1093/bioinformatics/btr330
DePaoli-Roach AA, Park IK, Cerovsky V, Csortos C, Durbin SD, Kuntz MJ, Sitikov A, Tang PM, Verin A, Zolnierowicz S (1994) Serine/threonine protein phosphatases in the control of cell function. Adv Enzyme Regul 34:199–224. https://doi.org/10.1016/0065-2571(94)90017-5
Eymard-Pierre E, Yamanaka K, Haeussler M, Kress W, Gauthier-Barichard F, Combes P, Cleveland DW, Boespflug-Tanguy O (2006) Novel missense mutation in ALS2 gene results in infantile ascending hereditary spastic paralysis. Ann Neurol 59(6):976–980. https://doi.org/10.1002/ana.20879
Fariello MI, Boitard S, Naya H, SanCristobal M, Servin B (2013) Detecting signatures of selection through haplotype differentiation among hierarchically structured populations. Genetics 193(3):929–941. https://doi.org/10.1534/genetics.112.147231
Frankham R, Ballou SEJD, Briscoe DA, Ballou JD (2002) Introduction to conservation genetics. Cambridge University Press, Cambridge
Fritz G, Kaina B (1999) Phosphorylation of the DNA repair protein APE/REF-1 by CKII affects redox regulation of AP-1. Oncogene 18(4):1033–1040. https://doi.org/10.1038/sj.onc.1202394
George L, Alex R, Sukhija N, Jaglan K, Vohra V, Kumar R, Verma A (2023) Genetic improvement of economic traits in Murrah buffalo using significant SNPs from genome-wide association study. Trop Anim Health Prod (3):199. https://doi.org/10.1007/s11250-023-03606-3
Gianola D, Simianer H, Qanbari S (2010) A two-step method for detecting selection signatures using genetic markers. Genet Res 92(2):141–155. https://doi.org/10.1017/S0016672310000121
Harris RB (2014) IUCN red list of threatened species: Cervus nippon. IUCN Red List of Threatened Species, UK
He Y, Wang ZH, Wang XM (2014) Genetic diversity and population structure of a Sichuan sika deer (Cervus sichuanicus) population in Tiebu Nature Reserve based on microsatellite variation. Dong Wu Xue Yan Jiu—Zool Res 35(6):528–536. https://doi.org/10.13918/j.issn.2095-8137.2014.6.528
Jaglan K, Ravikumar D, Sukhija N, George L, Alex R, Vohra V, Verma A (2023a) Genomic clues of association between clinical mastitis and SNPs identified by ddRAD sequencing in Murrah buffaloes. Anim Biotechnol 1–9. https://doi.org/10.1080/10495398.2023.2165937
Jaglan K, Surati U, Sukhija N, Mohan M, Jayakumar S, Verma A, Niranjan SK (2023b) Restriction site-associated DNA (RAD) sequencing: A novel and cost-effective genotyping approach using next generation sequencing. Cattle Pract 31(9):46–83
Jia B, Wang G, Zheng J, Yang W, Chang S, Zhang J, Liu Y, Li Q, Ge C, Chen G, Liu D, Yang F (2020) Development of novel EST microsatellite markers for genetic diversity analysis and correlation analysis of velvet antler growth characteristics in sika deer. Hereditas 157(1):24. https://doi.org/10.1186/s41065-020-00137-x
Kanaka KK, Sukhija N, Goli RC, Singh S, Ganguly I, Dixit SP, Dash A, Malik AA (2023) On the concepts and measures of diversity in the genomics era. Curr Plant Biol 33:100278. https://doi.org/10.1016/j.cpb.2023.100278
Kim S, Xu X, Hecht A, Boyer TG (2006) Mediator is a transducer of Wnt/beta-catenin signaling. J Biol Chem 281(20):14066–14075. https://doi.org/10.1074/jbc.M602696200
Koike N, Maita H, Taira T, Ariga H, Iguchi-Ariga SM (2000) Identification of heterochromatin protein 1 (HP1) as a phosphorylation target by Pim-1 kinase and the effect of phosphorylation on the transcriptional repression function of HP1(1). FEBS Lett 467(1):17–21. https://doi.org/10.1016/s0014-5793(00)01105-4
Kour A, Niranjan SK, Malayaperumal M, Surati U, Pukhrambam M, Sivalingam J, Kumar A, Sarkar M (2022) Genomic diversity profiling and breed-specific evolutionary signatures of selection in Arunachali yak. Genes. https://doi.org/10.3390/genes13020254
Kumar DR, Devadasan MJ, Surya T, Vineeth MR, Choudhary A, Sivalingam J, Verma A (2020) Genomic diversity and selection sweeps identified in Indian swamp buffaloes reveals it’s uniqueness with riverine buffaloes. Genomics 112(3):2385–2392
Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nat Methods 9(4):357–359. https://doi.org/10.1038/nmeth.1923
Li H (2011) A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27(21):2987–2993. https://doi.org/10.1093/bioinformatics/btr509
Liu H, Ju Y, Tamate H, Wang T, Xing X (2021) Phylogeography of sika deer (Cervus nippon) inferred from mitochondrial cytochrome-b gene and microsatellite DNA. Gene 772:145375. https://doi.org/10.1016/j.gene.2020.145375
Ma Y, Ding X, Qanbari S, Weigend S, Zhang Q, Simianer H (2015) Properties of different selection signature statistics and a new strategy for combining them. Heredity 115(5):426–436. https://doi.org/10.1038/hdy.2015.42
Malik AA, Sharma R, Ahlawat S, Deb R, Negi MS, Tripathi SB (2018) Analysis of genetic relatedness among Indian cattle (Bos indicus) using genotyping-by-sequencing markers. Anim Genet 49(3):242–245. https://doi.org/10.1111/age.12650
Messina G, Atterrato MT, Prozzillo Y, Piacentini L, Losada A, Dimitri P (2017) The human cranio facial development protein 1 (Cfdp1) gene encodes a protein required for the maintenance of higher-order chromatin organization. Sci Rep. https://doi.org/10.1038/srep45022
Mishra DC, Sikka P, Yadav S, Bhati J, Paul SS, Jerome A, Singh I, Nath A, Budhlakoti N, Rao AR, Rai A, Chaturvedi KK (2020) Identification and characterization of trait-specific SNPs using ddRAD sequencing in water buffalo. Genomics 112(5):3571–3578. https://doi.org/10.1016/j.ygeno.2020.04.012
Negishi M, Saraya A, Mochizuki S, Helin K, Koseki H, Iwama A (2010) A novel zinc finger protein Zfp277 mediates transcriptional repression of the Ink4a/arf locus through polycomb repressive complex 1. PloS one 5(8):e12373. https://doi.org/10.1371/journal.pone.0012373
Niu R, Shang HT, Wei H, Huang ZB, Zeng YZ (2001) Genetic analysis of 35 microsatellite loci in 5 lineages of Xishuangbanna miniature pig inbred line. Yi Chuan Xue Bao—Acta Genet Sin 28(6):518–526
Ohtaishi N, Gao Y (1990) A review of the distribution of all species of deer (Tragulidae, Moschidae and Cervidae) in China. Mammal Rev 20(2–3):125–144. https://doi.org/10.1111/j.1365-2907.1990.tb00108.x
Pavlidis P, Živković D, Stamatakis A, Alachiotis N (2013) SweeD: Likelihood-based detection of selective sweeps in thousands of genomes. Mol Biol Evol 30(9):2224–2234. https://doi.org/10.1093/molbev/mst112
Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE (2012) Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PloS one 7(5):e37135. https://doi.org/10.1371/journal.pone.0037135
Salzberg SL (2019) Next-generation genome annotation: we still struggle to get it right. Genome Biol 20(1):92. https://doi.org/10.1186/s13059-019-1715-2
Sivalingam J, Vineeth MR, Surya T, Singh K, Dixit SP, Niranjan SK, Tantia MS, Gupta ID, Ravikumar D (2020) Genomic divergence reveals unique populations among Indian Yaks. Sci Rep 10(1):3636. https://doi.org/10.1038/s41598-020-59887-3
Spee MJ, Albers MJ (1999) Genetic linkage: Concepts and methods. Genetics of cerebrovascular disease. Blackwell Publishing, Hoboken, pp 25–26
Stillwell EE, Zhou J, Joshi HC (2004) Human ninein is a centrosomal autoantigen recognized by CREST patient sera and plays a regulatory role in microtubule nucleation. Cell Cycle 3(7):923–930. https://doi.org/10.4161/cc.3.7.947
Štohlová Putnová L, Štohl R, Ernst M, Svobodová K (2021) A microsatellite genotyping-based genetic study of interspecific hybridization between the red and sika deer in the Western Czech Republic. Anim: Open Access J MDPI 11(6):1701. https://doi.org/10.3390/ani11061701
Sukhija N, Kanaka KK, Rangasai Chandra GR, Kapoor P, Sivalingam J, Verma A, Sharma R, Tripathi SB, Malik AA (2023a) The flight of chicken genomics and allied omics-a mini review. Ecol Genet Genomics. https://doi.org/10.1016/j.egg.2023.100201
Sukhija N, Malik AA, Devadasan JM, Dash A, Bidyalaxmi K, Ravi Kumar D, Kousalaya Devi M, Choudhary A, Kanaka KK, Sharma R, Tripathi SB (2023b) Genome-wide selection signatures address trait specific candidate genes in cattle indigenous to arid regions of India. Anim Biotech 1–5. https://doi.org/10.1080/10495398.2023.2290521
Tensen L (2016) Under what circumstances can wildlife farming benefit species conservation? Global Ecol Conserv 6:286–298. https://doi.org/10.1016/j.gecco.2016.03.007
Teshima KM, Coop G, Przeworski M (2006) How reliable are empirical genomic scans for selective sweeps? Genome Res 16(6):702–712. https://doi.org/10.1101/gr.5105206
Thévenon S, Thuy LT, Ly LV, Maudet F, Bonnet A, Jarne P, Maillard JC (2004) Microsatellite analysis of genetic diversity of the vietnamese sika deer (Cervus nippon pseudaxis). J Hered 95(1):11–18. https://doi.org/10.1093/jhered/esh001
Vineeth MR, Surya T, Sivalingam J, Kumar A, Niranjan SK, Dixit SP, Singh K, Tantia MS, Gupta ID (2020) Genome-wide discovery of SNPs in candidate genes related to production and fertility traits in Sahiwal cattle. Trop Anim Health Prod 52(4):1707–1715. https://doi.org/10.1007/s11250-019-02180-x
Xing X, Ai C, Wang T, Li Y, Liu H, Hu P, Wang G, Liu H, Wang H, Zhang R, Zheng J, Wang X, Wang L, Chang Y, Qian Q, Yu J, Tang L, Wu S, Shao X, Yang F (2022) The first high-quality reference genome of sika deer provides insights for high-tannin adaptation. Genomics Proteomics Bioinform 21:203–215
Ya HE, Zheng-Huan WANG, Xiao-Ming WANG (2014) Genetic diversity and population structure of a Sichuan sika deer (Cervus sichuanicus) population in Tiebu Nature Reserve based on microsatellite variation. Zool Res 35(6):528
Acknowledgements
The authors thank the Director, ICAR-National Dairy Research Institute for the support to carry out this research work. We would also like to extend our gratitude to Hengxing Ba for allowing us to use their data in our research/project. Without their contributions, this work would not have been possible.
Funding
The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.
Author information
Authors and Affiliations
Contributions
NS, MS and KK conceived concept. NS, RG, MS and KK analyzed the data. NS, RG, PB, AD and KJ wrote the main manuscript text. PR, MS and KJ prepared tables and figures. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interest
The authors have no relevant financial or non-financial interests to disclose.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Sukhija, N., Goli, R.C., Sukhija, M. et al. Evolutionary stamps for adaptation traced in Cervus nippon genome using reduced representation sequencing. Conservation Genet Resour 16, 135–146 (2024). https://doi.org/10.1007/s12686-023-01339-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s12686-023-01339-x