DOI: 10.3724/SP.J.1245.2014.00213

Asian Herpetological Research 2014/5:4 PP.213-227

Population Dynamics Following the Last Glacial Maximum in Two Sympatric Lizards in Northern China

Phylogeographic studies of Eremias lizards(Lacertidae) in East Asia have been limited, and the impact of major climatic events on their population dynamics remains poorly known. This study aimed to investigate population histories and refugia during the Last Glacial Maximum of two sympatric Eremias lizards(E. argus and E. brenchleyi) inhabiting northern China. We sequenced partial mitochondrial DNA from the ND4 gene for 128 individuals of E. argus from nine localities, and 46 individuals of E. brenchleyi from five localities. Forty-four ND4 haplotypes were determined from E. argus samples, and 33 from E. brenchleyi samples. Population expansion events began about 0.0044 Ma in E. argus, and 0.031 Ma in E. brenchleyi. The demographic history of E. brenchleyi indicates a long-lasting population decline since the most recent common ancestor, while that of E. argus indicates a continuous population growth. Among-population structure was significant in both species, and there were multiple refugia across their range. Intermittent gene flow occurred among expanded populations across multiple refugia during warmer phases of the glacial period, and this may explain why the effective population size has remained relatively stable in E. brenchleyi and grown in E. argus.

Key words:Lacertidae,Eremias lizards,mitochondrial DNA,historical demography,multiple refugia,Last Glacial Maximum

ReleaseDate:2015-03-27 09:32:49

Arnold E. N., Arribas O., Carranza S. 2007. Systematics of the Palaearctic and Oriental lizard tribe lacertini (Squamata:Lacertidae:Lacertinae), with descriptions of eight new genera. Zootaxa, 1430: 1-86

Avise J. C. 2000. Phylogeography: the history and formation of species. Cambridge: Harvard University Press

Bai W. N., Liao W. J., Zhang D. Y. 2010. Nuclear and chloroplast DNA phylogeography reveal two refuge areas with asymmetrical gene flow in a temperate walnut tree from East Asia. New Phytol, 188: 892-901

Bandelt H. J., Forster P., Röhl A. 1999. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol, 16: 37-48

Busack S. D., Lawson R., Arjo W. M. 2005. Mitochondrial DNA, allozymes, morphology and historical biogeography in the Podarcis vaucheri (Lacertidae) species complex. Amphibia-Reptilia, 26: 239-256

Camargo A., Sinervo B., Sites J. W. 2010. Lizards as model organisms for linking phylogeographic and speciation studies. Mol Ecol, 19: 3243-3488

Castoe T. A., Spencer C. L., Parkinson C. L. 2007. Phylogeographic structure and historical demography of the western diamondback rattlesnake (Crotalus atrox): a perspective on North American desert biogeography. Mol Phylogenet Evol, 42: 193-212

Chen B. H. 1991. Lacertidae. In Chen B. H. (Ed.), The Amphibian and Reptilian Fauna of Anhui. Hefei: Anhui Science and Technology Publishing House, 219-230

Clark A. M., Bowen B. W., Branch L. C. 1999. Effects of natural habitat fragmentation on an endemic scrub lizard (Sceloporus woodi): an historical perspective based on a mitochondrial DNA gene genealogy. Mol. Ecol., 8: 1093-1104

Ding G. H., Fu T. B., Ji X. 2012. Caudal autotomy does not increase locomotor costs in the oriental leaf-toed gecko Hemidactylus bowringii. Asian Herpetol Res, 3: 141-146

Drummond A. J., Rambaut A. 2007. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol, 7: 2-14

Dupanloup I., Schneider S., Excoffier L. 2002. A simulated annealing approach to define the genetic structure of populations. Mol Ecol, 11: 2571-2581

Elderkin C. L., Christian A. D., Vaughn C. C., Metcalfe-Smith J. L., Berg D. J. 2007. Population genetics of the freshwater mussel, Amblema plicata (Say 1817) (Bivalvia: Unionidae): Evidence of high dispersal and post-glacial colonization. Conserv Genet, 8: 355-372

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

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

Forcina G., Panayides P., Guerrini M., Nardi F., Gupta B. K., Mori E., Al-Sheikhly O. F., Mansoori J., Khaliq I., Rank D. N., Parasharya B. M., Khan A. A., Hadjigerou P., Barbanera F. 2012. Molecular evolution of the Asian francolins (Francolinus, Galliformes): a modern reappraisal of a classic study in speciation. Mol Phylogenet Evol, 65: 523-534

Fu J. Z. 1998.Toward the phylogeny of the family Lacertidae: implications from mitochondrial DNA 12S and 16S gene sequences (Reptilia: Squamata). Mol Phylogenet Evol, 9: 118-130

Fu Y. X. 1997. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics, 147: 915-925

Hall T. A. 1999. BIOEDIT: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser, 41: 95-98

Hewitt G. M. 1996. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc, 58: 247-293

Hewitt G. M. 1999. Post-glacial re-colonization of European biota. Biol J Linn Soc, 68: 87-112

Hewitt G. M. 2000. The genetic legacy of the quaternary ice ages. Nature, 405: 907-913

Hewitt G. M. 2004. Genetic consequences of climatic oscillations in the Quatemary. Philos Trans R Soc B, 359: 183-195

Hipsley C. A., Himmelmann L., Metzler D., Müller J. 2009. Integration of Bayesian molecular clock methods and fossil-based soft bounds reveals early Cenozoic origin of African lacertid lizards. BMC Evol Biol, 9: 151

Imbrie J., Hays J. D., Martinson D. G., McIntyre A., Mix A. C., Morley J. J. 1984. The orbital theory of Pleistocene climate: support from a revised chronology of the marine δ18O record. In Berger A. L., Imbrie Hays J. H., Kukla G., Saltzman B. (Eds.), Milankovitch and Climate, Part I. Reidel: Dordrecht, 269-305

Jin Y. T., Brown R. P., Liu N. F. 2008. Cladogenesis and phylogeography of the lizard Phrynocephalus vlangalii (Agamidae) on the Tibetan Plateau. Mol Ecol, 17: 1971-1982

Ju L. X., Wang H. J., Jiang D. B. 2007. Simulation of the Last Glacial Maximum climate over East Asia with a regional climate model nested in a general circulation model. Paleogeogr Paleoclimatol Paleoecol, 248: 376-390

Kimura M. 1981. Estimation of evolutionary distances between homologous nucleotide sequences. Proc Natl Acad Sci USA, 78: 454-458

Li J. J., Fang X. M., Ma H. Z., Zhu J. J., Pan B. T., Chen H. L. 1996. Geomorphological evolution of Late Cenozoic river from upper Yellow River and the uplift of Qinghai-Tibetan plateau on climatic changes. Sci China D, 26: 316-322

Li S. H., Yeung C. K. L., Feinstein J., Han L. X., Le M. H., Wang C. X. 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

Li Y. X., Xue X. X., Liu H. J. 2004. Fossil lizards of Qinling Mountains. Vertebr Palasiatic, 42: 171-176

Lin L. H., Ji X., Diong C. H., Du Y., Lin C. X. 2010. Phylogeography and population structure of the Reevese's butterfly lizard (Leiolepis reevesii) inferred from mitochondrial DNA sequences. Mol Phylogenet Evol, 56: 601-607

Lin Z. H., Ji X. 2005. Partial tail loss has no severe effects on energy stores and locomotor performance in a lacertid lizard, Takydromus septentrionalis. J Comp Physiol B, 175: 567-573

Lin Z. H., Qu Y. F., Ji X. 2006. Energetic and locomotor costs of tail loss in the Chinese skink, Eumeces chinensis. Comp Biochem Physiol A, 143: 508-513

Malhotra A., Dawson K., Guo P., Thorpe R. S. 2011. Phylogenetic structure and species boundaries in the mountain pitviper Ovophis monticola (Serpentes: Viperidae: Crotalinae) in Asia. Mol Phylogenet Evol, 59: 444-457

Martrat B., Grimalt J. O., Lopez-Martinez C., Cacho I., Sierro F. J., Flores J. A. 2004. Abrupt temperature changes in the Western Mediterranean over the past 250000 years. Science, 306: 1762-1765

Mayer W., Pavlicev M. 2007. The phylogeny of the family Lacertidae (Reptilia) based on nuclear DNA sequences: Convergent adaptations to arid habitats within the subfamily Eremiainae. Mol Phylogenet Evol, 44: 1155-1163

Mori E., Sforzi A., Di Febbraro M. 2013. From the Apennines to the Alps: recent range expansion of the crested porcupine Hystrix cristata L., 1758 (Mammalia: Rodentia: Hystricidae) in Italy. Ital J Zool, 80: 469-480

Newton A. C., Allnot T. R., Gillies A. C. M., Lowe A. J., Ennos R. A. 1999. Molecular phylogeography, intraspecific variation and conservation of tree species. Trends Ecol Evol, 14: 140-145

Nuñez J. J., Woood N. K., Rabanal F. E., Fontanella F. M., Sites J. W. 2011. Amphibian phylogeography in the Antipodes: refugia and postglacial colonization explain mitochondrial haplotype distribution in the Patagonian frog Eupsophus calcaratus (Cycloramphidae). Mol Phylogenet Evol, 58: 343-352

Nylander J. A. A. 2004. MRMODELTEST version 2.1. Computer program distributed by the author. Uppsala University, Uppsala.

Palumbi S. R. 1996. Nucleic acids II: the polymerase chain reaction. In Hillis D. M., Moritz C., Mable B. K. (Eds.), Moleclar Systematics, 2nd Edition. Massachusetts: Sinauer Associates, 205-247

Parsons T. J., Muniec D. S., Sullivan K, Woodyatt N, Alliston-Greiner R., Wilson M. R., Berry D. L., Holland K. A., Weedn V. W., Gill P., Holland M. M. 1997. A high observed substitution rate in the human mitochondrial DNA control region. Nat Genet, 15: 363-368

Petit R. J., Duminil J., Fineschi S., Hampe A., Salvini D., Vendramin G. G. 2005. Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations. Mol Ecol, 14: 689-701

Qian H., Ricklefs R. E. 2000. Large-scale processes and the Asian bias in temperate plant species diversity. Nature, 407: 180-182

Rambaut A., Drummond A. J. 2007. Tracer, version 1.4.

Rogers A. R., Fraley A. E., Bamshad M. J., Watkins W. S., Jorde L. B. 1996. Mitochondrial mismatch analysis is insensitive to the mutational process. Mol Biol Evol, 13: 895-902

Ronquist F., Huelsenbeck J. P. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics, 19: 1572-1574

Ruiz-González A., Madeira M. J., Randi E., Abramov A. V., Davoli F., Gómez-Moliner B. J. 2013. Phylogeography of the forest-dwelling European pine marten (Martes martes): new insights into cryptic northern glacial refugia. Biol J Linn Soc, 109: 118

Slatkin M. 1996. In defense of founder-flush theories of speciation. Am Nat, 147: 493-505

Stenson A. G., Thorpe R. S. 2003. Phylogeny, paraphyly and ecological adaptation of the colour and pattern in the Anolis roquet complex on Martinique. Mol Ecol, 12: 117-132

Sun Y. Y., Yang J, Ji X. 2009. Many-lined sun skinks (Mabuya multifasciata) do not compensate for the costs of tail loss by increasing feeding rate or digestive efficiency. J Exp Zool A, 311: 125-133

Tamura K., Stecher G., Peterson D., Filipski A., Kumar S. 2013. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol, 30: 2725-2729

Tian B., Liu R. R., Wang L.Y., Qiu Q., Chen K.M., Liu J. Q. 2009. Phylogeographic analyses suggest that a deciduous species (Ostryopsis davidiana Decne., Betulaceae) survived in northern China during the Last Glacial Maximum. J Biogeogr, 36: 2148-2155

Urquhart J., Wang Y. Z., Fu J. Z. 2009. Historical vicariance and male-mediated gene flow in the toad-headed lizards Phrynocephalus przewalskii. Mol Ecol, 18: 3714-3729

Weaver A. J., Eby M, Fanning A. F., Wiebe E. C. 1998. Simulated influence of carbon dioxide, orbital forcing and ice sheets on the climate of the Last Glacial Maximum. Nature, 394: 847-853

Wu C., Zhang X. Q., Ma Y. H. 1999. The Taihang and Yan mountains rose mainly in Quaternary. China Earthq Sci, 17: 17

Zamudio K. R., Savage W. K. 2003. Historical isolation, range expansion, and secondary contact of two highly divergent mitochondrial lineages in spotted salamanders (Ambystoma maculatum). Evolution, 57: 1631-1652

Zhang Z. Y., Yu Q. W., Zhang K. X., Gu Y. S., Xiang S. Y. 2003. Geomorphological evolution of Quaternary river from upper Yellow River and geomorphological evolution investigation for 1: 250000 scale geological mapping in Qinghai-Tibet plateau. J China Univ Geosci, 28: 621-633

Zhao K. T. 1999. Lacertidae. In Zhao E. M., Zhao K. T., Zhou K. Y. (Eds.), Fauna Sinica, Reptilia (Squamata: Lacertilia), Vol. 2. Beijing: Science Press, 219-242

Zhao Q., Liu H. X., Luo L. G., Ji X. 2011. Comparative population genetics and phylogeography of two lacertid lizards (Eremias argus and E. brenchleyi) from China. Mol Phylogenet Evol, 58: 478-491

Zhao Q., Wang Z., Liu L. L., Zhao W. G., Ji X. 2008. Selected body temperature, surface activity and food intake in tailed versus tailless Mongolian racerunners Eremias argus from three populations. Acta Zool Sinica, 54: 60-66

Zink R. M., Groth J. G., Vázquez-Miranda H., Barrowclough G. F. 2013. Phylogeography of the California Gnatcatcher (Polioptila californica) using multilocus DNA sequences and ecological niche modeling: implications for conservation. Auk, 130: 449-458