Bingjun DONG,Jing CHE,Li DING,et al.Testing Hypotheses of Pleistocene Population History Using Coalescent Simulations: Refugial Isolation and Secondary Contact in Pseudepidalea raddei (Amphibia: Bufonidae)[J].Asian Herpetological Research(AHR),2012,3(2):103-113.[doi:10.3724/SP.J.1245.2012.00103]
Click Copy

Testing Hypotheses of Pleistocene Population History Using Coalescent Simulations: Refugial Isolation and Secondary Contact in Pseudepidalea raddei (Amphibia: Bufonidae)
Share To:

Asian Herpetological Research[ISSN:2095-0357/CN:51-1735/Q]

2012 VoI.3 No.2
Research Field:
Original Article
Publishing date:


Testing Hypotheses of Pleistocene Population History Using Coalescent Simulations: Refugial Isolation and Secondary Contact in Pseudepidalea raddei (Amphibia: Bufonidae)
Bingjun DONG1 2 Jing CHE3 Li DING4 Song HUANG5 Robert W. MURPHY3 6 Ermi ZHAO1 4* and Yaping ZHANG3<
1 Key Laboratory for Bio-resources and Eco-environment (Ministry of Education), College of Life Sciences, Sichuan University, Chengdu 610064, Sichuan, China
2 College of Chemistry and Life Sciences, Shenyang Normal University, Shenyang 110034, Liaoning, China
3 State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, Yunnan, China
4 Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu 610064, Sichuan, China
5 Department of Biology, Huangshan University, Huangshan 245021, Anhui, China
6 Centre for Biodiversity and Conservation Biology, Department of Natural History, Royal Ontario Museum, 100 Queen’s Park, Toronto, ON, M5S 2C6, Canada
Mongolian toad Pleistocene refugia phylogeography biogeography
The impact of the Quaternary glaciation on eastern China’s local fauna and flora is a topic of considerable interest. We use mitochondrial DNA (mtDNA) sequences and coalescent simulations to test two general biogeographic hypotheses related to the effect of the Pleistocene climatic fluctuations on a widespread, eastern Chinese amphibian, Pseudepidalea raddei. Genealogical reconstructions are made and they detect major western and eastern lineages, which overlap in northwestern China, and possibly indicate the secondary contact of the populations that had entered the region from separate glacial refugia. Coalescent tests rejected alternative hypotheses of fragmentation of either a widespread ancestor or panmixia. The tests instead supported the hypothesis of geographic isolation and a remarkable dispersal pattern in one of the lineages. Though the Pleistocene climatic events are known to have affected the historical distributions and intra-specific divergence of Chinese squamates, coalescent and non-coalescent demographic analyses indicated that the toad P. raddei was not adversely affected by glacial cycling. Presumably, an increase in the amount of climatically mild habitats in East Asia is due to the development of monsoons since the Mid-late Pleistocene is responsible for the relatively mild effects.


Bauert R., Kalin M., Baltisberger M., Edwards P. J. 1998. No genetic variation detected within isolated relict populations of Saxifraga cernua in the Alps using RAPD markers. Mol Ecol, 7: 1559–1927
Bohonak A. J. 2002. IBD (Isolation by Distance): A program for analyses of isolation by distance. J Hered, 93: 153–154
Charlesworth D., Bartolome C., Schierup M. H., Mable B. K. 2003. Haplotype structure of the stigmatic self-incompatibility gene in natural populations of Arabidopsis lyrata. Mol Biol Evol, 20: 1741–1753
Ding L., Gan X. N., He S. P., Zhao E. M. 2011. A phylogeographic, demographic and historical analysis of the short-tailed pit viper (Gloydius brevicaudus): Evidence for early divergence and late expansion during the Pleistocene. Mol Ecol, 20: 1905–1922
Dong Z. M., Zhou J., Chang Z. J. 2004. Karyotype, C-banding and Ag-NORs of Bufo raddei from Henan. Chin J Zool, 39: 29–35 (In Chinese)
Drummond A. J., Rambaut A., Shapiro B., Pybus O. G. 2005. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol, 22: 1185–1192
Drummond A. J., Rambaut A. 2006. BEAST v 1.4. Available from
Drummond A. J., Ho S. Y. W., Phillips M. J., Rambaut A. 2006. Relaxed phylogenetics and dating with confidence. PLoS Biol, 4: e88
Excoffier L., Smouse P. E., Quattro J. M. 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics, 131: 479–491
Excoffier L., Laval G., Schneider S. 2005. Arlequin version 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online, 1: 47–50
Fontanella F. M., Feldman C. R., Siddall M. E. 2008. Phylogeography of Diadophis punctatus: Extensive lineage diversity and repeated patterns of historical demography in a transcontinental snake. Mol Phylogenet Evol, 46: 1049–1070
Frost D. R., Grant T., Faivovich J., Bain R. H., Haas A., Haddad C. F. B., de Sá, R. O., Channing A., Wilkinson M., Donnellan S. C., Raxworthy C. J., Campbell J. A., Blotto B. L., Moler P., Drewes R. C., Nussbaum R. A., Lynch J. D., Green M. D., Wheeler W. C. 2006. The amphibian tree of life. Bull Am Mus Nat Hist, 297: 1–370
Fu Y. X. 1997. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics, 147: 915–925
Goebel A. M., Donnelly J. M., Atz M. E. 1999. PCR primers and amplification methods of 12S ribosomal DNA, control region, and an overview of PCR primers which have amplified DNA in amphibians successfully. Mol Phylogenet Evol, 11: 163–199
Guiher T. J., Burbrink F. T. 2008. Demographic and phylogeographic histories of two venomous North American snakes of the genus Agkistrodon. Mol Phylogenet Evol, 48: 543–553
Harpending H. 1994. Signature of ancient population growth in a low resolution mitochondrial mismatch distribution. Humn Biol, 66: 131–137
Hewitt G. M. 2000. The genetic legacy of the Quaternary ice ages. Nature, 405: 907–913
Hey J., Nielsen R. 2004. Multilocus methods for estimating population sizes, migration rates, and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persmilis. Genetics, 167: 747–760
Hey J., Nielsen R. 2007. Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics. PNAS, 104: 2785–2790
Ho S. Y. W., Larson G. 2006. Molecular clocks: When times are a-changin'. Trends Genet. 22: 79–83
Huang S., He S. P., Peng Z. G., Zhao K., Zhao E. M. 2007. Molecular phylogeography of endangered sharp-snouted pitviper (Deinagkistrodon acutus; Reptilia, Viperidae) in Mainland China. Mol Phylogenet Evol, 44: 942–952
Hu Y. L., Wu X. B., Jiang Z. G., Peng Y., Su X., Cao S. Y. 2007. Population genetics and phylogeography of Bufo gargarizans in China. Biochem Genet, 45: 697–711
Huelsenbeck J. P., Ronquist F. R. 2001. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics, 17: 754–755
Huelsenbeck J. P., Rannala B. 2004. Frequentist properties of Bayesian posterior probabilities of phylogenetic trees under simple and complex substitution models. Syst Biol, 53: 904–913
Hull J., Girman D. J. 2005. Effects of Holocene climate change on the historical demography of migrating sharp-shinned hawks (Accipiter striatus velox) in North America. Mol Ecol, 14: 159–170
Kimura M. 1983. The neutral theory of molecular evolution. Cambridge, Massachusetts: Cambridge Univ Press
Knowles L. L. 2001. Did the Pleistocene glaciations promote divergence? Tests of explicit refugial models in montane grasshoppers. Mol Ecol, 10: 691–701
Knowles L. L., Maddison W. P. 2002. Statistical phylogeography. Mol Ecol, 11: 2623–2635
Kuchta S. R., Tan A. M. 2005. Isolation by distance and post-glacial range expansion in the rough-skinned newt, Taricha granulose. Mol Ecol, 14: 225–244
Kvist L., Viiri K., Dias P. C., Rytk?nen S., Orell M. 2004. Glacial history and colonization of Europe by the blue tit Parus caeruleus. J Avian Biol, 35: 352–359
Li S. H., Yeung C., Feinstein J., Han L. X., Le M. H., Wang C. X., Ding P. 2009. Sailing through the Late Pleistocene: Unusual historical demography of an East Asian endemic, the Chinese hwamei (Leucodioptron canorum canorum), during the last glacial period. Mol Ecol, 18: 622–633
Liu D. S., Li Z. G. 1996. Geography of Asia. Beijing: Commercial Press
Maddison W. P., Maddison D. R. 2004. Mesquite: A modular system for evolutionary analysis. Ver. 1.01. Available from
Maddison W. P., Maddison D. R. 2008. Mesquite: A modular system for evolutionary analysis. Ver. 2.5. Available from
Mantel N. 1967. The detection of disease clustering and a generalized regression approach. Cancer Res, 27: 209–220
Milá B., Girman D. J., Kimura M., Smith T. B. 2000. Genetic evidence for the effect of a postglacial population expansion on the phylogeography of a North American songbird. Proc Royal Soc B: Biol Sci, 267: 1033–1040
Milá B., Smith T. B., Wayne R. K. 2006. Postglacial population expansion drives the evolution of long-distance migration in a songbird. Evolution, 60: 2403–2409
Milá B., McCormack J. E., Casta?eda G., Wayne R. K., Smith T. B. 2007. Recent postglacial range expansion drives the rapid diversification of a songbird lineage in the genus Junco. Proc Royal Soci B: Biol Sci, 274: 2653–2660
Nei M., Tajima F. 1983. Maximum likelihood estimation of the number of nucleotide substitutions from restriction-site data. Genetics, 105: 207–217
Nielsen R., Wakeley J. 2001. Distinguishing migration from isolation: A Markov Chain Monte Carlo approach. Genetics, 158: 885–896
Posada D., Crandall K. A. 1998. MODELTEST: Testing the model of DNA substitution. Bioinformatics, 14: 817–818
Rambaut A., Drummon A. 2003. Tracer: A program for analysing results from Bayesian MCMC programs such as BEAST & MrBayes. Oxford, U.K. Available from
Rogers A. R., Harpending H. 1992. Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol, 9: 552–569
Rozas J., Rozas R. 1999. DnaSP Ver 3: An integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics, 15: 174–175
Rozas J., Delbarrio J., Messeguer X., Rozas R. 2003. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics, 19: 2496–2497
Sambrook J., Fritsch E. F., Maniatis T. 1989. Molecular cloning. A Laboratory Manual. New York: Cold Spring Harbor Laboratory Press
Schmitt T. 2007. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Front Zool, 4: 11
Shang K. G., Li W. J., Jiang H. Y., Zeng Y. L. 1982. The karyotype of Bufo raddei Strauch I. The cultivation of blood leucocytes of Bufo raddei and an analysis of its karyotype. Acta Genet Sin, 9: 298–302 (In Chinese)
Shi Y. F., Zheng B. X., Li S. J., Ye B. S. 1995. Studies on altitude and climatic environment in the middle and east parts of Tibetan Plateau during Quaternary maximum glaciation. J Glacial Geocryol, 17: 97–112 (In Chinese)
Slatkin M., Maddison W. P. 1989. A cladistic measure of gene flow inferred from the phylogenies of alleles. Genetics, 123: 603–613
Slatkin M., Hudson R. R. 1991. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics, 129: 555–562
Song G., Qu Y. H., Yin Z. H., Li S. H., Liu N. F., Lei F. M. 2009. Phylogeography of the Alcippe morrisonia (Aves: Timaliidae): Long population history beyond late Pleistocene glaciations. BMC Evol Biol, 9: 143–153
St?ck M., Moritz C., Hickerson M., Frynta D., Dujsebayeva T., Eremchenko V., Macey J. R., Papenfuss T. J., Wake D. B. 2006. Evolution of mitochondrial relationships and biogeography of Palearctic green toads (Bufo viridis subgroup) with insights intheir genomic plasticity. Mol Phylogenet Evol, 41: 663–689
Swofford D. L. 2002. PAUP*: Phylogenetic analysis using parsimony and other methods, Ver. 4.0b10, Computer program distributed by Sinauer Associates, Inc., Sunderland, MA
Tajima F. 1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 123: 585–595
Thompson J. D., Gibson T. J., Plewniak F., Jeanmougin F., Higgins D. G. 1997. The CLUSTAL X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res, 25: 4876–4882
Wang Y. J., Wang X. F., Wang X. Z., Li J., Wang Z. S., Chen W. Y. 1997. Chromosome homology between Bufo bufo gargarizans and Bufo raddei. Zool Res, 18: 415–419 (In Chinese)
Wu G. L., Pan B. T., Guan Q. Y., Gao H. S. 2002. The maximum glaciation and desert expansion in China during MIS16. J Glaciol Geocryol, 24: 544–549
Yang H. Y. 1982. The karyotype of the toad Bufo raddei Strauch. Acta Zool Sin, 28: 28–34 (In Chinese)
Yu D., Zhao W. G. 2007. Breeding ecology of Bufo raddei in Harbin region. Sichuan J Zool, 26: 409–411 (In Chinese)
Zhang R. Z. 1999. Zoogeography of China. Beijing: Science Press (In Chinese)
Zhang R. Z. 2002. Geological events and mammalian distribution in China. Acta Zool Sin, 48: 141–153 (In Chinese)
Zhang H., Yan J., Zhang G. Q., Zhou K. Y. 2008. Phylogeography and demographic history of Chinese black-spotted frog populations (Pelophylax nigromaculata): Evidence for independent refugia expansion and secondary contact. BMC Evol Biol, 8: 1–16
Zhao E. M., Adler K. 1993. Herpetology of China. Oxford: Society for the Study of Amphibians and Reptiles
Zhou L. Z., Song J. Y. 1997. A study of feeding ecology of Bufo raddei. Chin J Ecol, 16: 29–34 (In Chinese)


Last Update: 2016-03-15