Abstract
The Sichuan snub-nosed monkey (Rhinopithecus roxellana) is an endangered species endemic to China, where the smallest population resides in Shennongjia National Park, Hubei Province. In this study, the genetic diversity of the Sichuan snub-nosed monkeys from two areas, Qianjiaping (QJP) and Dalongtan (DLT) in Shennongjia National Park was evaluated using single nucleotide polymorphisms (SNPs) data derived from restriction site-associated DNA sequencing (RAD-seq). A total of 41,260 SNP loci were found in 29 Sichuan snub-nosed monkey individuals. The genome-wide nucleotide diversity (π) of the Shennongjia population was 0.001842. The genetic differentiation (FST) between the QJP and DLT subpopulations was 0.034. The heterozygosity of individuals from QJP was 0.3475 ± 0.03696 and 0.3148 ± 0.03501 for individuals from DLT. Although the DLT and QJP subpopulations did not show significant genetic differences, genetic differentiation between the two subpopulations was confirmed using Bayesian cluster analysis, neighbor-joining trees and principal component analysis. These results suggest that the Shennongjia population of Sichuan snub-nosed monkey has relatively low genetic diversity at the genomic level. The little genetic differentiation noted between the DLT and QJP subpopulations likely due to natural and anthropogenic barriers which may exacerbate loss of genetic diversity of this endangered subpopulation.
Similar content being viewed by others
Data accessibility
Raw sequence data for this study will be available at the NCBI Sequence Read Archive after paper publishing, accession PRJNA471594.
References
Alexander DH, Lange K (2011) Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinform 12(1)
Alexander DH, Novembre J, Lange K (2009) Fast model-based estimation of ancestry in unrelated individuals. Genome Res 19(9):1655–1664
Amores A, Catchen J, Ferrara A, Fontenot O, Postlethwait JH (2011) Genome evolution and meiotic maps by massively parallel DNA sequencing: spotted gar, an outgroup for the teleost genome duplication. Genetics 188:799–808
Andrews KR, Good JM, Miller MR, Luikart G, Hohenlohe PA (2016) Harnessing the power of RADseq for ecological and evolutionary genomics. Nat Rev Genet 17:81–92
Attard CRM, Möller LM, Sasaki M, Hammer MP, Bice CM, Brauer CJ, Carvalho DC, Harris JO, Beheregaray LB (2016) A novel holistic framework for genetic-based captive-breeding and reintroduction programs: holistic framework for reintroductions. Conserv Biol 30:1060–1069
Blanco-Bercial L, Bucklin A (2016) New view of population genetics of zooplankton: RAD-seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus. Mol Ecol 25:1566–1580
Cariou M, Duret L, Charlat S (2013) Is RAD-seq suitable for phylogenetic inference? An in silico assessment and optimization. Ecol Evol 3:846–852
Catchen JM, Amores A, Hohenlohe P, Cresko W, Postlethwait JH (2011) Stacks: building and genotyping loci de novo from short read sequences. G3 1:171–182
Chang Z, Liu Z, Yang J, Li M, Vigilant L (2012a) Noninvasive genetic assessment of the population trend and sex ratio of the Shennongjia population of Sichuan snub-nosed monkeys (Rhinopithecus roxellana). Chin Sci Bull 57:1135–1141
Chang ZF, Luo MF, Liu ZJ, Yang JY, Xiang ZF, Li M, Vigilant L (2012b) Human influence on the population decline and loss of genetic diversity in a small and isolated population of Sichuan snub-nosed monkeys (Rhinopithecus roxellana). Genetica 140:105–114
Chapman CA, Peres CA (2001) Primate conservation in the new millennium: the role of scientists. Evol Anthropol Issues News Rev 10:16–33
Chutimanitsakun Y, Nipper RW, Cuesta-Marcos A, Cistue L, Corey A, Filichkina T et al (2011) Construction and application for QTL analysis of restriction site associated DNA (RAD) linkage map in barley. BMC Genomics 12:4
Clayton D (2012) snpStats: SnpMatrix and XSnpMatrix classes and methods. R Package
Corander J, Majander KK, Cheng L, Merilä J (2013) High degree of cryptic population differentiation in the Baltic Sea herring Clupea harengus. Mol Ecol 22(11):2931–2940
Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, McVean G, Durbin R (2011) The variant call format and VCFtools. Bioinformatics 27:2156–2158
Davey JW, Blaxter ML (2010) RADSeq: next-generation population genetics. Brief Funct Genomics 9:416–423
Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, Blaxter ML (2011) Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet 12:499–510
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ (2011) A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 43:491–498
Eaton DA (2014) PyRAD: assembly of de novo RADseq loci for phylogenetic analyses. Bioinformatics 30:1844–1849
Evans BJ, Zeng K, Esselstyn JA, Charlesworth B, Melnick DJ (2014) Reduced representation genome sequencing suggests low diversity on the sex chromosomes of Tonkean Macaque monkeys. Mol Biol Evol 31:2425–2440
Felsenstein J (2005) PHYLIP (Phylogeny Inference Package) version 3.6. Distributed by the author. Department of Genome. Sciences, University of Washington, Seattle
He L, Zhang Y, Peng H, Li D, Li D (2010) Genetic diversity of Rhinopithecus roxellana in Shennongjia National Nature Reserve as estimated by non-invasive DNA technology. Acta Ecol Sin 30:4340–4350
He W, Zhao S, Liu X, Dong S, Lv J, Liu D, Wang J, Meng Z (2013) ReseqTools: an integrated toolkit for large-scale next generation sequencing based resequencing analysis. Genet Mol Res 12:6275–6283
Hohenlohe PA, Amish SJ, Catchen J, Allendorf FW, Luikart G (2011) Next-generation RAD sequencing identifies thousands of SNPs for assessing hybridization between rainbow and westslope cutthroat trout. Mol Ecol Resour 11:117–122
Ingvarsson PK, Baxter SW, Davey JW, Johnston JS, Shelton AM, Heckel DG, Jiggins CD, Blaxter ML (2011) Linkage mapping and comparative genomics using next-generation RAD sequencing of a non-model organism. PLoS ONE 6:e19315
Ishtiaq F, Prakash V, Green RE, Johnson J (2015) Management implications of genetic studies for ex situ populations of three critically endangered Asian Gyps vultures. Anim Conserv 18:259–270
Ivy JA, Lacy RC (2010) Using molecular methods to improve the genetic management of captive breeding programs for threatened species. In: DeWoody JA, Bickham JW, Michler CH, Nichols KM, Rhodes OE, Woeste KE (eds) Molecular approaches in natural resource conservation and management. Cambridge University Press, Cambridge, pp 267–295
Ivy JA, Miller A, Lacy RC, DeWoody JA (2009) Methods and prospects for using molecular data in captive breeding programs: an empirical example using Parma wallabies (Macropus parma). J Hered 100:441–454
Jeffries DL, Copp GH, Lawson Handley L, Olsén KH, Sayer CD, Hänfling B (2016) Comparing RADseq and microsatellites to infer complex phylogeographic patterns, an empirical perspective in the Crucian carp, Carassius carassius, L. Mol Ecol 25:2997–3018
Khan S, Nabi G, Ullah MW, Yousaf M, Manan S, Siddique R, Hou H (2016) Overview on the role of advance genomics in conservation biology of endangered species. Int J Genomics 2016:1–8
Labuschagne C, Nupen L, Kotzé A, Grobler PJ, Dalton DL (2015) Assessment of microsatellite and SNP markers for parentage assignment in ex situ African Penguin (Spheniscus demersus) populations. Ecol Evol 5:4389–4399
Lew RM, Finger AJ, Baerwald MR, Goodbla A, May B, Meek MH (2015) Using next-generation sequencing to assist a conservation hatchery: a single-nucleotide polymorphism panel for the genetic management of endangered delta smelt. Trans Am Fish Soc 144:767–779
Li Y (2004) The effect of forest clear-cutting on habitat use in Sichuan snub-nosed monkey (Rhinopithecus roxellana) in Shennongjia Nature Reserve, China. Primates 45:69–72
Li H, Durbin R (2010) Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26(5):589–595
Li M, Liu Z, Gou J, Ren B, Pan R, Su Y, Funk SM, Wei F (2007) Phylogeography and population structure of the golden monkeys (Rhinopithecus roxellana): inferred from mitochondrial DNA sequences. Am J Primatol 69:1195–1209
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R (2009) The sequence alignment/map format and SAMtools. Bioinformatics 25:2078–2079
Liu Z, Liu G, Roos C, Wang Z, Xiang Z, Zhu P, Wang B, Ren B, Shi F, Pan H, Li M (2015) Implications of genetics and current protected areas for conservation of 5 endangered primates in China. Conserv Biol 29:1508–1517
Long Y, Richardson M (2008) Rhinopithecus roxellana. The IUCN red list of threatened species 2008: e.T19596A8985735. Available from https://doi.org/10.2305/IUCN.UK.2008.RLTS.T19596A8985735.en
Luo M, Pan H, Liu Z, Li M (2012) Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (Rhinopithecus roxellana). BMC Evol Biol 12:207
Malukiewicz J, Guschanski K, Grativol AD, Oliveira MA, Ruiz-Miranda CR, Stone AC (2016) Application of PE-RADSeq to the study of genomic diversity and divergence of two Brazilian marmoset species (Callithrix jacchus and C. penicillata). Am J Primatol 79:1–12
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297–1303
Nichols KM, Kozfkay CC, Narum SR (2016) Genomic signatures among Oncorhynchus nerka ecotypes to inform conservation and management of endangered Sockeye Salmon. Evol Appl 9:1285–1300
Oates JF (2013) Primate conservation: unmet challenges and the role of the international primatological society. Int J Primatol 34:235–245
Pan D, Hu HX, Meng SJ, Men ZM, Fu YX, Zhang YP (2009) A high polymorphism level in Rhinopithecus roxellana. Int J Primatol 30:337–351
Pan Y, Wang X, Sun G, Li F, Gong X (2016) Application of RAD sequencing for evaluating the genetic diversity of domesticated Panax notoginseng (Araliaceae). PLoS ONE 11:e0166419
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81:559–575
Scally A, Yngvadottir B, Xue Y, Ayub Q, Durbin R, Tyler-Smith C (2013) A genome-wide survey of genetic variation in gorillas using reduced representation sequencing. PLoS ONE 8:e65066
The R Core Team (2016) R: a language and environmental for statistical computing. Reference Index. https://cran.r-project.org/doc/manuals/r-release/R-intro.pdf
Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J, Banks E, Garimella KV, Altshuler D, Gabriel S, DePristo MA (2013) From FastQ data to high-confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinform 43:11.10.1–11.10.33
Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38:1358–1370
Willoughby JR, Ivy JA, Lacy RC, Doyle JM, DeWoody JA (2017) Inbreeding and selection shape genomic diversity in captive populations: implications for the conservation of endangered species. PLoS ONE 12:e0175996
Wragg D, Marti-Marimon M, Basso B, Bidanel JP, Labarthe E, Bouchez O, Conte YL, Vignal A (2016) Whole-genome resequencing of honeybee drones to detect genomic selection in a population managed for royal jelly. Sci Rep 6:27168
Xiang Z, Yu Y, Yang M, Yang J, Niao M, Li M (2011) Does flagship species tourism benefit conservation? A case study of the golden snub-nosed monkey in Shennongjia National Nature Reserve. Chin Sci Bull 56(24):2553–2558
Yang J, Yang W (2018) Study and conservation of Rhinopithecus roxellanae and its habitat in Shennongjia. China Forestry Publishing House, Beijing
Yang J, Liao M, Yu H, Yao H (2008) Protection and research status of golden monkeys in Shennongjia. World Sci Tech R & D 30:418–421
Yang J, Lee SH, Goddard ME, Visscher PM (2011) GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet 88:76–82
Yang W, Zhang B, Xia M, Yao H, Liu J, Yu H, Xu H, Yang J, Xiang Z (2018) The anesthetic effects of ketamine on golden snub-nosed monkeys (Rhinopithecus roxellana). Acta Theriol Sin 38:295–301
Zhang BD, Xue DX, Wang J, Li YL, Liu BJ, Liu JX (2016) Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD-seq for the small yellow croaker (Larimichthys polyactis). Mol Ecol Resour 16:755–768
Zhang Y, Zhong Y, Hong Y, Xue Y, Li D, Zhou C, Liu S (2019) RAD-seq data advance captive-based conservation of wild Bactrian camels (Camelus ferus). Conserv Genet. https://doi.org/10.1007/s10592-019-01173-5
Zhou X, Wang B, Pan Q, Zhang J, Kumar S, Sun X, Liu Z, Pan H, Lin Y, Liu G, Zhan W, Li M, Ren B, Ma X, Ruan H, Cheng C, Wang D, Shi F, Hui Y, Tao Y, Zhang C, Zhu P, Xiang Z, Jiang W, Chang J, Wang H, Cao Z, Jiang Z, Li B, Yang G, Roos C, Garber PA, Bruford MW, Li R, Li Ming (2014) Whole-genome sequencing of the snub-nosed monkey provides insights into folivory and evolutionary history. Nat Genet 46:1303–1310
Zhou X, Meng X, Liu Z, Chang J, Wang B, Li M, Orozco-terWengel P, Tian S, Wen C, Wang Z, Garber PA, Pan H, Ye X, Xiang Z, Bruford MW, Edwards SV, Cao Y, Yu S, Gao L, Cao Z, Liu G, Ren B, Shi F, Peterfi Z, Li D, Li B, Jiang Z, Li J, Gladyshev VN, Li R, Li M (2016) Population genomics reveals low genetic diversity and adaptation to hypoxia in snub-nosed monkeys. Mol Biol Evol 33:2670
Zhou Y, Xue Y, Zhang Y, Yu H, Yang L, Jiang J, Zhang Y (2018) Genetic structure of the golden snub-nosed monkey in Shennongjia National Natural Reserve based on microsatellite DNA markers. Acta Ecol Sin 38:1404–1410
Acknowledgements
This study was funded by the Key Technology Program of Hubei Province (2019AC161) and Shennongjia National Park (2019SNJ003), and National Key Technology R&D Program of China (2013BAD03B02). The authors would like to thank Jingyuan Yang, Wanji Yang and Hui Yao for the fieldwork assistance and Dr. Fei Gao, Mr. Dongming Fang, Dr. Yue Hu, Dr. Qi Hu and Ph.D. student Du Zhang at Shenzhen Agricultural Genome Research Institute, Chinese Academy of Agricultural Sciences for their support of this project. The authors would like to thank Charlotte Hacker for assistance in editing this paper.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Zhang, Y., Zhou, Y., Liu, X. et al. Genetic diversity of the Sichuan snub-nosed monkey (Rhinopithecus roxellana) in Shennongjia National Park, China using RAD-seq analyses. Genetica 147, 327–335 (2019). https://doi.org/10.1007/s10709-019-00073-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10709-019-00073-3