Figures
Abstract
To evaluate the conservation status of a species or population it is necessary to gain insight into its ecological requirements, reproduction, genetic population structure, and overall genetic diversity. In our study we examined the genetic diversity of Rhinopithecus brelichi by analyzing microsatellite data and compared them with already existing data derived from mitochondrial DNA, which revealed that R. brelichi exhibits the lowest mitochondrial diversity of all so far studied Rhinopithecus species. In contrast, the genetic diversity of nuclear DNA is high and comparable to other Rhinopithecus species, i.e. the examined microsatellite loci are similarly highly polymorphic as in other species of the genus. An explanation for these differences in mitochondrial and nuclear genetic diversity could be a male biased dispersal. Females most likely stay within their natal band and males migrate between bands, thus mitochondrial DNA will not be exchanged between bands but nuclear DNA via males. A Bayesian Skyline Plot based on mitochondrial DNA sequences shows a strong decrease of the female effective population size (Nef) starting about 3,500 to 4,000 years ago, which concurs with the increasing human population in the area and respective expansion of agriculture. Given that we found no indication for a loss of nuclear DNA diversity in R. brelichi it seems that this factor does not represent the most prominent conservation threat for the long-term survival of the species. Conservation efforts should therefore focus more on immediate threats such as development of tourism and habitat destruction.
Citation: Kolleck J, Yang M, Zinner D, Roos C (2013) Genetic Diversity in Endangered Guizhou Snub-Nosed Monkeys (Rhinopithecus brelichi): Contrasting Results from Microsatellite and Mitochondrial DNA Data. PLoS ONE 8(8): e73647. https://doi.org/10.1371/journal.pone.0073647
Editor: Tom Gilbert, Natural History Museum of Denmark, Denmark
Received: June 3, 2013; Accepted: July 24, 2013; Published: August 29, 2013
Copyright: © 2013 Kolleck et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The Margot Marsh Biodiversity Foundation and the German Primate Center financially supported this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Many animal species are threatened by extinction and for a number of species long-term survival is questionable without a proper conservation and population management [1]–[5]. Population management is crucial for several species, in particular to keep or increase genetic diversity of a population, e.g., to reduce inbreeding (e.g., Delacour's langur [Trachypithecus delacouri] in Vietnam [6]; grey wolf [Canis lupus] in the northern Rocky Mountains, United States [7]; and Siberian jay [Perisoreus infaustus] in Finland [8]). Important prerequisites for conservation and population management are information on population size, population dynamics, genetic diversity, genetic population structure, the mating system, and the current and historical distribution [6], [9]–[15].
Nearly 50% of the World's primate species are threatened (classified as vulnerable, endangered or critically endangered on the IUCN Red List [3]), among them the Guizhou snub-nosed monkey (Rhinopithecus brelichi). This endangered species is endemic to China and it occurs only in a single population of less than 900 individuals in the Fanjingshan National Nature Reserve (FNNR), Guizhou Province of southern China [16]–[20]. Compared to the closely related species R. roxellana and R. bieti, R. brelichi shows a low diversity of mitochondrial DNA (mtDNA) [21]. Since mtDNA is only maternally inherited, variation in this marker reflects only a part of the overall genetic variation within a population. A low diversity of mtDNA compared to related species describes only a low female effective population size (Nef) sometime in the past and does not give feasible information about the total genetic variation of the species. To get a better understanding of the genetic diversity of a population or species, a combined analysis of nuclear DNA (nDNA) and mtDNA diversity is more powerful. Since for R. brelichi only data on mtDNA diversity is available [21], we genotyped wild R. brelichi samples at eight microsatellite loci to gain information on nDNA diversity. Since the population size of R. brelichi is relatively small, we expect to find a similar impoverished nuclear genetic diversity as for mtDNA, which then would make a population genetic management plan more important [22].
Materials and Methods
Ethics statement
Research complied with protocols approved by the Forest Ministry of Guizhou Province, China, and adhered to the legal requirements of China and Germany, where research was conducted. Permits to collect samples were provided by the State Forestry Administration of China and Fanjingshan National Nature Reserve, and the staff of the nature reserve helped to collect samples and supported the study. The research was carried out in compliance with respective animal care regulations and the principles of the American Society of Primatologists and the German Primate Center for the ethical treatment of nonhuman primates. Fecal samples from wild individuals were collected non-invasively after the animals left the site without disturbing, threatening or harming them in their natural behavior.
Field site and sample collection
The FNNR, established in 1986, covers the main peaks of the Wuling Mountains in Guizhou Province of southern China (27°45′10′′–28°02′24′′N, 108°35′03′′–108°48′23′′E, 480 – 2,750 m, Figure 1). Its size is just 38,300 ha (25,000 ha core zone, 13,300 ha buffer zone). About 15,000 people live in the buffer zone and about 1,100 in the core area of the reserve [19]. The vegetation in the area is subtropical, semi-evergreen mixed and evergreen deciduous forest. In the area also coniferous forest is present, but R. brelichi is rarely found there. The species occurs manly in secondary forest, because primary forest is almost completely gone [20].
Collection of samples for genetic analyses was carried out in the gray region around the Yangaoping field station.
The core distribution of the monkeys is in the northern part of FNNR close to the field station Yangaoping (27°58′54.37′′N, 108°54′36.04′′E, 1,627 m, Figure 1). About 400–500 monkeys live around the station [18], [19], [23]. In this area, a total of 146 R. brelichi samples were collected during surveys in 2007 and 2008 [21]. Details about sample collection are described in Yang et al. [21].
Laboratory methods
DNA extraction followed methods outlined in Yang et al. [21]. To prevent contamination, all steps were performed in a sterile workbench (Captair® flow) and after each sample the workbench was cleaned with ethanol and decontaminated with UV-light. Concentration of DNA was measured on a NanoDrop® ND-1000 spectrophotometer (NanoDrop Technologies). All 146 samples were genotyped at eight microsatellite loci (D1S533, D2S1326, D6S264, D6S501, D7S2204, D8S505, D10S1432, and D17S1290) using fluorescent-labelled primers (FAM, TET, or TAMRA) (Table 1). PCR reactions were carried out in a total volume of 30 µl, including 18.6–19.6 µl H2O (HPLC quality), 3 µl 10x buffer (contains 15 mM MgCl2, Biotherm), 4 µl BT buffer (contains 10 mg/ml BSA, 0.7% Triton X-100), 1 µl of each primer (10 µM), 0.2 µl dNTPs (25 mM), 0.2 µl Taq DNA polymerase (Biotherm, 5 u/µl) and 1–2 µl (50–100 ng) DNA. The PCR profile was 2 min at 94°C; then 40 cycles of 30 seconds at 94°C, 30 seconds at varying annealing temperatures (Table 1), 30 seconds at 72°C; followed by 5 min at 72°C. PCR amplification success was checked by visualization of 1 μl of PCR product under UV light after electrophoresis on 1% agarose gels containing ethidium bromide. Concentration of PCR products was estimated by comparison with diluted concentration standards (Fermentas pUC19 DNA 0.5 µg/µl) of 5 ng/μl and 10 ng/μl, respectively. For analysis, 1–5 ng PCR product was taken and analyzed on an ABI 3130xL Genetic Analyser (Applied Biosystems) along with the Genescan™ 400 HD [Rox™] Size Standard and Hi-Di™ Formamid, according to manufactures instructions (8.775 µl Hi-Di™, 0.225 µl standard, and 1 µl PCR product). GeneMapper® v3.7 software (Applied Biosystems) was applied to determine allele sizes and genotypes. To avoid false inferences, the genotyping was repeated multiple times following the approach proposed by Morin et al. [24]. For heterozygous at least three repetitions, and for homozygous three to seven repetitions were performed from independent DNA extractions. Since fecal samples contain mainly degraded and low amounts of nDNA, allelic dropout and null alleles could lead to incorrect genotypes [24]. Thus, Micro-Checker v2.2.3 [25] was applied to test for the presence of possible null alleles, large allelic dropout or scoring errors due to stuttering.
Microsatellite analysis
First, we tested if there is any population structuring in our data set using a Bayesian approach in STRUCTURE v2.3.4 [26]. Program settings were set to the admixture model and a total run length of 1,000,000 iterations, a burnin of 100,000 and values of K (assumed number of populations) from 1–3. GENPOP v4.0.10 was applied to test for linkage disequilibrium, using default settings (demomorization number: 1,000; number of batches: 100; number of iterations per batch: 1,000) [27], [28]. We used GenAlEx v6.5 [29] for calculating allelic range, observed and effective number of alleles, observed (Ho) and expected (He) heterozygosities, and fixation index (Fis) to measure deviation from Hardy-Weinberg Equilibrium (HWE) based on heterozygosity excess, with significance of deviations from zero assessed through bootstrapping with 1,000 replications across loci to generate 95% confidence intervals [30].
A significant negative Fis value indicates heterozygous excess, contrarily, a significantly positive Fis value indicates heterozygous deficit (homozygote excess) in the population. Hamilton [30] suggested that a deficit of heterozygosity provides evidence for consanguineous mating taking place in a population, while heterozygosity excess indicates a reduction of population size.
To detect indication for a recent bottleneck or a reduction of population size in R. brelichi, we used Bottleneck v1.2.02 [31]. We performed the one-tailed Wilcoxon signed rank test to test for significant excess of heterozygotes. Piry et al. [31] suggest that the Wilcoxon test gives most reliable results when less than 20 polymorphic loci are studied. When a shrinking population gets small, the allele number is faster reduced than the expected heterozygosity. Thus, the heterozygosity excess (He > Heq: He based on allele frequencies; Heq: heterozygosity expected at mutation equilibrium, based on the number of alleles) was tested with the Two-Phase Model (TPM) [32], which combines the Infinite Allele Model (IAM) [33] with the Stepwise Mutation Model (SMM) [34]. The true model of mutation for most loci is between the IAM and SMM [35] and since the SMM is more suitable for microsatellite data than the IAM [36], we selected the TPM method with 90% single-step and just 10% multiple-step mutations. A significant heterozygosity excess (He > Heq) indicates a recent bottleneck event [36]. For the TPM, 10,000 iterations were applied. The variance was set to 12 as recommended [31]. We also applied a qualitative graphical method for detecting recent bottlenecks based on a mode-shift distortion of the allele frequency distribution in Bottleneck. Bottlenecks cause alleles at low frequency (<0.1) to become less abundant than alleles in one or more intermediate allele frequency class (e.g., 0.1–0.2) thus shifting the mode of the normally L-shaped frequency distribution into higher frequency classes [35], [37], [38].
Mitochondrial DNA analysis
To infer timing and magnitude of past changes in R. brelichi Nef, we performed a Bayesian Skyline Plot (BSP) in BEAST v1.7.4 [39], [40]. Therefore, we used sequences of the hypervariable region I (HVI) of the mitochondrial control region published by Yang et al. [21]. In a total of 141 individuals (as determined by genotyping), only five haplotypes were found (HQ891105-HQ891109) with a frequency of the most dominant haplotype of 0.723 [21]. The haplotype network presented in Yang et al. [21] and the fact that there are obviously no physical barriers for primates in FNNR suggest that the R. brelichi population at FNNR is not sub-structured. Hence, false inferences from the BSP due to population structure are unlikely [41].
For the BSP, we applied the HKY model [42] as selected by the Bayesian Information Criterion in jModelTest v0.1.1 [43], [44]. The clock model was set to a strict clock with a generation time of 12 years according to a female age of first reproduction of 8–9 years, and a lifespan of about 18 years [45]. We used a substitution rate of 0.1643 substitutions per nucleotide per million years [21] and a coalescent Bayesian Skyline tree prior with 10 groups under a piecewise-constant skyline model. The length of chains was set to 50 million generations with logged parameters every 1,000 generations. Tracer v1.5.0 was applied to check chain convergence and to visualize the BSP.
Results
Microsatellites
Of the 146 collected samples, genotyping showed that they derived from 141 different individuals. We found no evidence for null alleles, large allelic dropout, scoring errors due to stuttering or sub-structuring in our samples. The observed number of alleles (Na), effective number of alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), and the fixation index (Fis) are shown in Table 2. The Na per locus ranges from 3 to 13 with a mean of 6.63±3.34. The mean Ne resulted in 3.84±1.43 ranging between 2.23 and 5.77. The mean He across eight polymorphic loci (0.70±0.12) is similar to the mean Ho (0.71±0.12). All values indicate that all eight microsatellite loci are highly polymorphic in the R. brelichi population (Table 2). The values for Ho and He for all microsatellite loci are not different (p = 0.85). Thus, no deviation from HWE is indicated. Similarly, the overall Fis value for the eight microsatellite loci is not significantly different from zero (p = 0.72) also suggesting that the R. brelichi population does not depart from HWE. However, D1S533 (Fis = 0.14) and D17S1290 (Fis = 0.22) show significant positive Fis values (p<0.05) indicating heterozygosity deficits. A test for linkage disequilibrium revealed no significant linkage pattern.
A comparison of He and Ho values of R. bieti (data taken from Liu et al. [46]) and R. brelichi for six of our eight loci (D1S533, D2S1326, D6S264, D7S2204, D8S505, and D17S1290) revealed no significant differences between both species (R. brelichi: He: 0.72±0.13; Ho: 0.71±0.14; R. bieti: He: 0.78±0.09; Ho: 0.58±0.27; He: t = −1.5, p = 0.19; Ho: t = 0.94, p = 0.39; Table 3).
The TPM implemented in the Bottleneck program does not support a recent bottleneck in the R. brelichi population (p = 0.098). Additionally the allele frequency distribution does not indicate any distortion form a normal L-shape distribution (Figure 2), hence also this graphic test does not support any recent bottleneck of the R. brelichi population.
Bars represent the proportion of alleles found in each allele frequency class. The distribution is L-shaped, as expected for a stable population under mutation-drift equilibrium, thus not indicating a recent bottleneck.
Mitochondrial DNA
To infer timing and magnitude of past changes in R. brelichi Nef, we performed a BSP using mitochondrial HVI region sequence data of the 141 individually assigned samples. The BSP (Figure 3) displays changes in the Nef through time and indicates a reduction in population size between 3,500 and 4,000 years ago.
Calculations are based on 603 bp of the mitochondrial HVI region. Shown are the median (black) and the 95% highest posterior probability density (dashed lines) around the estimate. The arrow indicates the start of a reduction in Nef.
Discussion
As shown by Yang et al. [21], R. brelichi (haplotype diversity [H]: 0.457±0.084; nucleotide diversity [π]: 0.014±0.007) exhibits a low diversity in mtDNA as compared to its congeners R. roxellana (H: 0.845±0.026; π: 0.034±0.017) and R. bieti (H: 0.948±0.006; π: 0.036±0.018). In contrast to the low diversity in mtDNA, diversity in nDNA does not differ from other Rhinopithecus species. All eight microsatellite loci were highly polymorphic as suggested by the high allelic richness (6.63±3.34) and heterozygosity (He: 0.70±0.12; Ho: 0.71±0.12). Similar allelic richness (7.5) and heterozygosity (He: 0.70; Ho: 0.61) were found in R. bieti [47] with a total population of less than 2000 individuals [3] and in R. roxellana (allelic richness: 4.5; He: 0.59; Ho: 0.59) [48] with a total population in Shennongjia Nature Reserve of approx. 500 individuals [49].
Since only one panmixing population of R. brelichi remained, we could not compare the genetic diversity among (sub-)populations to detect a loss of the overall genetic diversity, but we matched our results with respective data from R. bieti [46]. Both the expected and the observed rates of heterozygosity do not differ largely between the two species. Pan et al. [50] indicated that R. roxellana is highly polymorphic at microsatellite loci. Thus, with similar heterozygosity values, a similar high level of genetic diversity as in R. roxellana can be found in R. brelichi implying no indication for a reduction in terms of nDNA diversity in this species.
A heterozygosity deficit was detected at two of the eight loci: D1S533 (Fis = 0.14) and D17S1290 (Fis = 0.22). The deviation from HWE might be caused by null alleles or consanguineous mating taking place in this population [30]. However, the results from Micro-Checker do not show any evidence for null alleles. Since the population is relatively small and mating opportunities are limited, even more so if certain males are able to monopolize a number of females, on one hand consanguineous mating might happen frequently in R. brelichi and reduces their fitness (inbreeding depression) [51]. On the other hand R. brelichi could have a strong natural avoidance for inbreeding, which could lead to a decline of progeny and so to a decrease of population size [52].
The TPM provides no sufficient evidence for a recent bottleneck in the R. brelichi population. With a sample size of more than 15% of the total population we can exclude the weakness of Bottleneck tests caused by limited sample sizes [53]. Similarly false results caused by deviation from HWE are unlikely [35]. In the analysis of allele frequencies we found the most abundant alleles in the lowest frequency class (<0.1). The overall distribution of allele proportions in various frequency classes follows a normal L-shape (Figure 2). This supports the assumption that no bottleneck occurred in the recent history of R. brelichi. The high genetic diversity of nDNA and lack of evidence for a recent bottleneck strongly suggests that the R. brelichi population was in equilibrium during its recent history.
The analysis of mtDNA draws another picture. Although the BSP suggests a stable female effective population size (Nef) in R. brelichi over at least 25,000 to 35,000 years, a strong decline is indicated 3,500 to 4,000 years ago (Figure 3). A similar population trend was found in R. bieti [47] and R. roxellana [54]. According to geographic provenances of fossils and historical records, all snub-nosed monkey species in China were once widespread. Simultaneously with the decline of Nef in the three Chinese snub-nosed monkeys, the human population increased and with it, significant modifications of the environment due to agricultural practices started, thus largely reducing suitable habitat for snub-nosed monkeys in China [55], [56].
The obvious contrast
High nDNA diversity combined with low mtDNA diversity is uncommon and not often observed in mammals. One example is the Scandinavian brown bear (Ursus arctos), which was close to extinction in the 1930s [57]. In the bear population, the nDNA diversity is much larger than the mtDNA diversity [57]. Waits and colleagues [57] suggested as a possible reason for this discrepancy the greater maximum dispersal distances of male bears in contrast to females. This would lead to a localized distribution of mtDNA haplotypes but a broader geographic distribution of nDNA variants. A similar scenario can be envisioned for R. brelichi. In historical times, when the species had a much larger distribution with several demes [56], female philopatry and male migration would contribute to a mtDNA population structure consisting of local populations with certain haplotype assemblages, since males do not transfer their mitochondria. On the nDNA level however, male migration would contribute to gene flow and a population with no or only a weak spatial genetic structure (Figure 4). After the extinction of almost all demes, many mtDNA haplotypes were lost and only those few in the surviving deme remained. In contrast, the remaining population would nevertheless comprise a large proportion of the former nDNA diversity (Figure 4).
In historical times, R. brelichi was distributed over a larger area than today comprising several subpopulations or demes (A–G) with respective different mtDNA haplogroup assemblages. Due to male migration between these demes, nDNA was transferred between them and equalized nDNA diversity among demes, but not so for mtDNA. After partial habitat and population loss in this example only one deme survived (E; dashed circles), containing just a subset of the original mtDNA haplotypes but almost all nDNA diversity. Thus, mtDNA diversity was strongly reduced while nDNA diversity remained comparatively high.
Our results suggest that the overall genetic diversity of the endangered R. brelichi is not reduced to a degree that it might become a direct threat for the long-term survival of the species. In contrast, due to the very limited distribution and population size of R. brelichi, conservation strategies should be focused on direct threats to this species (such as tourist infrastructure development, deforestation, grazing of livestock, using the reserve as source for food and fire wood) and try to remove or limit them. An enlargement of the reserve would be advantageous, but it seems even more important to establish at least one additional, but geographically isolated population [3]. A translocation of some individuals into another habitat might minimize the risk that a stochastic event like a forest fire or a disease could lead to extinction of the species. There could be several sites within the Wuling Mountains that might offer such habitats [3], but thorough surveys and habitat suitability analysis are required to find them.
Acknowledgments
We are grateful to the staff of the Fanjingshan National Nature Reserve in Guizhou province, China, for helping to collect samples and supporting this study. We also thank Fanjingshan National Nature Reserve for permits to collect samples and the State Forestry Administration of China for research permits. Christiane Schwarz is thanked for her excellent technical assistance in the laboratory.
Author Contributions
Conceived and designed the experiments: DZ CR. Performed the experiments: JK MY. Analyzed the data: JK MY DZ CR. Contributed reagents/materials/analysis tools: JK MY DZ CR. Wrote the paper: JK DZ CR.
References
- 1. Ellstrand NC, Elam DR (1993) Population genetic consequences of small population size: Implications for plant conservation. Annu Rev Ecol Evol Syst 24: 217–242.
- 2.
Frankham R, Ballou JD, Briscoe DA (2002) Introduction to conservation genetics. Cambridge: Cambridge University. 642 p.
- 3.
IUCN (2012) IUCN red list of threatened species 2012.2. Available: http://www.iucnredlist.org. Accessed: 05 March 2013.
- 4. Kirkpatrick RC (1995) The natural history and conservation of the snub-nosed monkeys (genus Rhinopithecus). Biol Conserv 72: 363–369.
- 5. Vucetich JA, Waite TA (2003) Spatial patterns of demography and genetic processes across the species’ range: Null hypotheses for landscape conservation genetics. Conserv Genet 4: 639–645.
- 6. Ebenau A, Nadler T, Zinner D, Roos C (2011) Genetic population structure of the critically endangered Delacour's langur (Trachypithecus delacouri) in Van Long Nature Reserve, Vietnam. Vietn J Primatol 1(5): 1–15.
- 7. vonHoldt BM, Stahler DR, Bangs EE, Smith DW, Jimenez MD, et al. (2010) A novel assessment of population structure and gene flow in grey wolf populations of the northern Rocky Mountains of the United States. Mol Ecol 19: 4412–4427.
- 8. Li MH, Merilä J (2010) Extensive linkage disequilibrium in a wild bird population. Heredity 104: 600–610.
- 9. Austin JD, Jelks HL, Tate B, Johnson AR, Jordan F (2011) Population genetic structure and conservation genetics of threatened Okaloosa darters (Etheostoma okaloosae). Conserv Genet 12: 981–989.
- 10. Chaves PB, Alvarenga CS, Possamai CB, Dias LG, Boubli JP, et al. (2011) Genetic diversity and population history of a critically endangered primate, the northern muriqui (Brachyteles hypoxanthus). PLoS ONE 6: e20722.
- 11.
Lande R, Barrowclough GF (1987) Effective population size, genetic variation, and their use in population management. In: Soulé ME, editor. Variable populations for conservation. Cambridge: Cambridge University Press. 87–123.
- 12. Mace GM (2004) The role of taxonomy in species conservation. Philos Trans R Soc London B Biol Sci 359: 711–719.
- 13. Mbora DNM, McPeek MA (2010) Endangered species in small habitat patches can possess high genetic diversity: The case of the Tana river red colobus and mangabey. Conserv Genet 11: 1725–1735.
- 14. Plesnar-Bielak A, Skrzynecka AM, Prokop ZM, Radwan J (2012) Mating system affects population performance and extinction risk under environmental challenge. Proc Biol Sci 279: 4661–4667.
- 15. Schwartz MK, Luikart G, Waples RS (2007) Genetic monitoring as a promising tool for conservation and management. Trends Ecol Evol 22: 25–33.
- 16. Bleisch W (1992) Preliminary comments on the Guizhou snub-nosed monkey (Rhinopithecus brelichi). Asian Primates 1: 4.
- 17.
Jablonski NG (1998) The natural history of the doucs and snub-nosed monkeys. Singapore: World Scientific. 382 p.
- 18. Xiang ZF, Nie SG, Lei XP, Chang ZF, Wei FW, et al. (2009) Current status and conservation of the gray snub-nosed monkey Rhinopithecus brelichi (Colobinae) in Guizhou, China. Biol Conserv 142: 469–476.
- 19.
Yang YQ, Lei XP, Yang CD (2002) Fanjingshan research: Ecology of the wild Guizhou snub-nosed monkey (Rhinopithecus brelichi). Guiyang: Guizhou Science Press. 186 p.
- 20.
Bleisch W, Yongcheng L, Richardson M (2008) Rhinopithecus brelichi. In: IUCN 2012. IUCN red list of threatened species 2012.2. Available: http://www.iucnredlist.org. Accessed 05 March 2013.
- 21. Yang M, Yang Y, Cui D, Fickenscher G, Zinner D, et al. (2012) Population genetic structure of Guizhou snub-nosed monkeys (Rhinopithecus brelichi) as inferred from mitochondrial control region sequences, and comparison with R. roxellana and R. bieti. Am J Phys Anthropol 147: 1–10.
- 22. Frankham R (1996) Relationship of genetic variation to population size in wildlife. Conserv Biol 10: 1500–1508.
- 23.
Sun DY, Li MJ, Zhang WY (1990) A report of investigation on non-human primate resources in Tongren region. In: Zhou ZH, Yang YQ, editors. Research on the Fanjing Mountain. Guiyang: People’s Publishing House. 442–456.
- 24. Morin PA, Chambers KE, Boesch C, Vigilant L (2001) Quantitative polymerase chain reaction analysis of DNA from non-invasive samples for accurate microsatellite genotyping of wild chimpanzees (Pan troglodytes verus). Mol Ecol 10: 1835–1844.
- 25. van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P (2004) Micro-Checker: Software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4: 535–538.
- 26. Pritchard J, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.
- 27. Raymond M, Rousset F (1995) Genepop (version-1.2) – population-genetics software for exact tests and ecumenicism. J Hered 86: 248–249.
- 28. Rousset F (2008) Genepop ‘007: a complete re-implementation of the Genepop software for Windows and Linux. Mol Ecol Resour 8: 103–106.
- 29. Peakall R, Smouse PE (2006) GenAlEx 6: Genetic analysis in Excel: Population genetic software for teaching and research. Mol Ecol Notes 6: 288–295.
- 30.
Hamilton MB (2009) Population genetics. Oxford: Wiley-Blackwell. 407 p.
- 31. Piry S, Luikart G, Conuet JM (1999) Bottleneck: A computer program for detecting recent reduction in the effective population size using allele frequency data. J Hered 90: 502–503.
- 32. Di Rienzo A, Peterson AC, Garza JC, Valades AM, Slatkin M, et al. (1994) Mutational process of simple-sequence repeat loci in human populations. Proc Natl Acad Sci USA 91: 3166–3170.
- 33. Kimura M, Crow JF (1964) The number of alleles that can be maintained in a finite population. Genetics 49: 725–738.
- 34. Ota T, Kimura M (1973) A model of mutation appropriate to estimate the number of electrophoretically detectable alleles in a finite population. Genet Res 22: 201–204.
- 35. Luikart G, Cornuet JM (1998) Empirical evaluation of a test for identifying recently bottlenecked populations from alleles frequencies data. Conserv Biol 12: 228–237.
- 36. Cornuet JM, Luikart G (1996) Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 144: 2001–2014.
- 37. Luikart G, Sherwin WB, Steele BM, Allendorf FW (1998) Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change. Mol Ecol 7: 963–974.
- 38. Luikart G, Allendorf FW, Cornuet JM, Sherwin WB (1998) Distortion of allele frequency distributions provides a test for recent population bottlenecks. J Hered 89: 238–247.
- 39. Drummond AJ, Rambaut A, Shapiro B, Pybus OG (2005) Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 22: 1185–1192.
- 40. Drummond AJ, Suchard MA, Xie D, Rambaut A (2012) Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 29: 1969–1973.
- 41. Heller R, Chikhi L, Siegismund HR (2013) The confounding effect of population structure on Bayesian skyline plot inferences of demographic history. PLoS ONE 8: e62992.
- 42. Hasegawa M, Kishino H, Yano T (1985) Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol 22: 160–174.
- 43. Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696–704.
- 44. Posada D (2008) jModelTest: Phylogenetic model averaging. Mol Biol Evol 25: 1253–1256.
- 45. Yang M, Sun DY, Zinner D, Roos C (2009) Reproductive parameters in Guizhou snub-nosed monkeys (Rhinopithecus brelichi). Am J Primatol 71: 266–270.
- 46. Liu ZJ, Ren BP, Hao YL, Zhang HR, Wei FW, et al. (2008) Identification of 13 human microsatellite markers via cross-species amplification of fecal samples from Rhinopithecus bieti. Int J Primatol 29: 265–272.
- 47. Liu ZJ, Ren BP, Wu RD, Zhao L, Hao YL, et al. (2009) The effect of landscape features on population genetic structure in Yunnan snub-nosed monkeys (Rhinopithecus bieti) implies an anthropogenic genetic discontinuity. Mol Ecol 18: 3831–3846.
- 48. Chang ZF, Liu ZJ, Yang JY, Li M, Vigilant L (2012) Noninvasive genetic assessment of the population trend and sex ratio of the Shennongjia population of Sichuan snub-nosed monkeys (Rhinopithecus roxellana). Chinese Science Bulletin 57: 1135–1141.
- 49. Zhu ZQ (2003) Ecological research of golden monkeys (Rhinopithecus roxellana) in Shennongjia Nature Reserve. Hubei Forestry Science and Technology Z1: 46–52.
- 50. Pan D, Li Y, Hu HX, Meng SJ, Men ZM, et al. (2005) Microsatellite polymorphism of Sichuan golden monkeys. Chinese Science Bulletin 50: 2850–2855.
- 51. Bijlsma R, Bundgaard J, Boerema AC (2000) Does inbreeding affect the extinction risk of small populations: Predictions from Drosophila. J Evol Biol 13: 502–514.
- 52. IG, Taylor SS, Tracy LN, Kokko H, Armstrong DP (2009) Why some species of birds do not avoid inbreeding: Insights from New Zealand robins and saddlebacks. Behav Ecol 20: 575–584.
- 53. Peery MZ, Kirby R, Reid BN, Stoelting R, Doucet-Beer E (2012) Reliability of genetic bottleneck tests for detecting recent population declines. Mol Ecol 21: 3403–3418.
- 54. Lou MF, Liu ZJ, Pan HJ, Zhao L, Li M (2012) Historical geographic dispersal of the golden snub-nosed monkey (Rhinopithecus roxellana) and the influence of climatic oscillations. Am J Primatol 74: 91–101.
- 55. Zhou W, Yu X, Jull AJT, Burr G, Xiao JY, et al. (2004) High-resolution evidence from southern China of an early Holocene optimum and a mid-Holocene dry event during the past 18,000 years. Quatern Res 62: 39–48. Corrigendum (2008) Quatern Res 70: 129–130.
- 56. Li B, Pan R, Oxnard CE (2002) Extinction of snub-nosed monkeys in China during the past 400 years. Int J Primatol 23: 1227–1244.
- 57. Waits L, Taberlet P, Swenson JE, Sandegren F, Franzén R (2000) Nuclear DNA microsatellite analysis of genetic diversity and gene flow in the Scandinavian brown bear (Ursus arctos). Mol Ecol 9: 421–431.