Comparative phylogeography of four Apodemus species (Mammalia: Rodentia) in the Asian Far East: evidence of Quaternary climatic changes in their genetic structure

Current Status of Trypanosoma grosi and Babesia microti in Small Mammals in the Republic of Korea

Small mammals, such as rodents and shrews, are natural reservoir hosts of zoonotic diseases, including parasitic protozoa. To assess the risk of rodent-borne parasitic protozoa in the Republic of Korea (ROK), this study investigated the status of parasitic protozoa, namely Trypanosoma, Babesia, and Theileria, in small mammals. In total, 331 blood samples from small mammals were analyzed for parasites using PCR and sequenced. Samples were positive for Trypanosoma grosi (23.9%; n = 79) and Babesia microti (10%; n = 33) but not Theileria. Small mammals from Seogwipo-si showed the highest infection rate of T. grosi (48.4%), while the highest B. microti infection rate was observed in those from Gangneung-si (25.6%). Sequence data revealed T. grosi to be of the AKHA strain. Phylogenetic analysis of B. microti revealed the US and Kobe genotypes. B. microti US-type-infected small mammals were detected throughout the country, but the Kobe type was only detected in Seogwipo-si. To our knowledge, this is the first nationwide survey that confirmed T. grosi and B. microti infections at the species level in small mammals in the ROK and identified the Kobe type of B. microti. These results provide valuable information for further molecular epidemiological studies on these parasites.

Spatio-temporal genetic structure of the striped field mouse (Apodemus agrarius) populations inhabiting national parks in South Korea: Implications for conservation and management of protected areas

Population or habitat connectivity is a key component in maintaining species and community-level regional biodiversity as well as intraspecific genetic diversity. Ongoing human activities cause habitat destruction and fragmentation, which exacerbate the connectivity due to restricted animal movements across local habitats, eventually resulting in the loss of biodiversity. The Baekdudaegan Mountain Range (BMR) on the Korean Peninsula represents “biodiversity hotspots” and eight of the 22 Korean national parks are located within the BMR. Given the striped field mouse (Apodemus agrarius) is the most common and ecologically important small mammals in these protected areas, the population genetic assessment of this species will allow for identifying “genetic diversity hotspots” and also “genetic barriers” that may hinder gene flow, and will therefore inform on effective conservation and management efforts for the national park habitats. We collected samples from hair, tail, or buccal swabs for 252 A. agrarius individuals in 2015 and 2019. By using mitochondrial DNA cytochrome b (cyt b) sequences and nine microsatellite loci, we determined levels of genetic diversity, genetic differentiation, and gene flow among eight national park populations of A. agrarius along the BMR. We found high levels of genetic diversity but the occurrences of inbreeding for all the nine samples analyzed. Our results also indicated that there was detectable temporal genetic variation between the 2015 and 2019 populations in the Jirisan National Park, which is probably due to a short-term decline in genetic diversity caused by reduced population sizes. We also found a well-admixed shared gene pool among the national park populations. However, a significant positive correlation between geographic and genetic distances was detected only in mtDNA but not microsatellites, which might be attributed to different dispersal patterns between sexes. There was a genetic barrier to animal movements around the Woraksan National Park areas. The poor habitat connectivity surrounding these areas can be improved by establishing an ecological corridor. Our findings of the presence of genetic barriers in some protected areas provide insights into the conservation and management efforts to improve the population or habitat connectivity among the national parks.

Phylogeography of the striped field mouse (Apodemus agrarius Pallas, 1771) in light of new data from central part of Northern Eurasia

A phylogeographic analysis of A. agrarius based on the complete mtDNA cytochrome b and control region sequences has been performed using data obtained for the first time for the species from large regions of the central part of Northern Eurasia (23 localities of Altai, Western Siberia, and the Urals). The obtained results have demonstrated a complex intraspecific differentiation of A. agrarius, which has manifested not only in the isolation of the isles populations in Southeast Asia (Jeju and Taiwan), but also in the genetic heterogeneity of mainland populations, which has reflected the history of the modern intraspecific genetic diversity formation against the background of changing physiographic conditions of Eurasia in the Quaternary. The divergence of genetic lineages has taken place apparently simultaneously (in mid-Pleistocene) on the territory of the Eastern part of the modern disjunctive range, where all the identified lineages are present today. The demographic history and possible evolutionary scenarios for A. agrarius in the Western part of the range have been considered. TMRC reconstructions have shown that the lifetime of the common ancestor of the lineage that expanded in the Western Palearctic is about 17.7 [95% HPD 13.2–22.5] kyr. This suggests that the transcontinental expansion of A. agrarius is a relatively recent event that has occurred after the LGM.

Distinct body-size responses to warming climate in three rodent species

In mammals, body-size responses to warming climates are diverse, and the mechanisms underlying these different responses have been little investigated. Using temporal and spatial datasets of three rodent species distributed across different climatic zones in China, we investigated temporal and spatial trends of body size (length and mass), identified the critical drivers of these trends, and inferred the potential causes underlying the distinct body-size responses to the critical drivers. We found that body mass of all species remained stable over time and across space. Body length, however, increased in one species over time and in two species across space. Generally, body-length variation was predicted best by minimum ambient temperature. Moreover, in two species, body length changed linearly with temperature differences between ancestral and colonization areas. These distinct temperature-length patterns may jointly be caused by species-specific temperature sensitivities and experienced magnitudes of warming. We hypothesize that species or populations distributed across distinct temperature gradients evolved different intrinsic temperature sensitivities, which affect how their body sizes respond to warming climates. Our results suggest that size trends associated with climate change should be explored at higher temporal and spatial resolutions, and include clades of species with similar distributions.

Phylogeographic study of the Bufo gargarizans species complex, with emphasis on Northeast Asia

We conduct a phylogeographic and population genetic study of the Asiatic toad (Bufo gargarizans) to understand its evolutionary history, and the influence of geology and climate. A total of 292 individuals from 94 locations were genotyped for two mitochondrial loci (cytb, ND2) and five nuclear introns (Sox9-2, Rho-3, CCNB2-3, UCH-2, and DBI-2), and we performed a suite of phylogenetic, population genetic, and divergence dating analyses. The phylogenetic trees constructed using mitochondrial loci inferred B. gargarizans being divided into two major groups: China mainland and Northeast Asia (Northeast China, Russia, and Korean Peninsula). As with previous studies of this species, we recover population genetic structure not tied to geographic region. Additionally, we discover a new genetic clade restricted to Northeast Asia that points towards the Korean Peninsula being a glacial refugium during the Pleistocene. The weak phylogeographic pattern of B. gargarizans is likely the result of multiple biological, anthropogenic, and historical factors – robust dispersal abilities as a consequence of physiological adaptations, human translocation, geologic activity, and glacial cycles of the Pleistocene. We highlight the complex geologic and climatic history of Northeast Asia and encourage further research to understand its impact on the biodiversity in the region.

Late Pleistocene Expansion of Small Murid Rodents across the Palearctic in Relation to the Past Environmental Changes

We investigated the evolutionary history of the striped field mouse to identify factors that initiated its past demographic changes and to shed light on the causes of its current genetic structure and trans-Eurasian distribution. We sequenced mitochondrial cyt b from 184 individuals, obtained from 35 sites in central Europe and eastern Mongolia. We compared genetic analyses with previously published historical distribution models and data on environmental and climatic changes. The past demographic changes displayed similar population trends in the case of recently expanded clades C1 and C3, with the glacial (MIS 3–4) expansion and postglacial bottleneck preceding the recent expansion initiated in the late Holocene and were related to environmental changes during the upper Pleistocene and Holocene. The past demographic trends of the eastern Asian clade C3 were correlated with changes in sea level and the formation of new land bridges formed by the exposed sea shelf during the glaciations. These data were supported by reconstructed historical distribution models. The results of our genetic analyses, supported by the reconstruction of the historical spatial distributions of the distinct clades, confirm that over time the local populations mixed as a consequence of environmental and climatic changes resulting from cyclical glaciation and the interglacial period during the Pleistocene.

Ring distribution patterns-diversification or speciation? Comparative phylogeography of two small mammals in the mountains surrounding the Sichuan Basin

Studying the genetic differentiation in a unique geographical area contributes to understanding the process of speciation. Here, we explore the spatial genetic structure and underlying formation mechanism of two congeneric small mammal species (Apodemus draco and A. chevrieri), which are mainly distributed in the mountains surrounding the lowland Sichuan Basin, southwest China. We applied a set of comparative phylogeographical analyses to determine their genetic diversification patterns, combining mitochondrial (Cytb and COI) and nuclear (microsatellite loci) markers, with dense sampling throughout the range (411 A. draco from 21 sites and 191 A. chevrieri from 22 sites). Moreover, we performed three complementary statistical methods to investigate the correlation between genotype and geographical and environmental components, and predicted the potential suitable distributional range under the present and historical climate conditions. Our results suggest that both species have experienced allopatric differentiation and admixture in historical periods, resulting in a ring-shape diversification, under the barrier effect of the Sichuan Basin. We infer that the tectonic events of the Qinghai-Tibetan Plateau and climatic oscillations during the Quaternary played an important role on the genetic divergence of the two species by providing environmental heterogeneity and geographical variation. Our study reveals a case of two sympatric small mammals following a ring-shaped diversification pattern and provides insight into the process of differentiation.

Contrasting phylogeographic histories between the continent and islands of East Asia: Massive mitochondrial introgression and long-term isolation of hares (Lagomorpha: Lepus)

Hares of the genus Lepus are distributed worldwide, and introgressive hybridization is thought to be pervasive among species, leading to reticulate evolution and taxonomic confusion. Here, we performed phylogeographic analyses of the following species of hare across East Asia: L. timidus, L. mandshuricus, L. coreanus, and L. brachyurus collected from far-eastern Russia, South Korea, and Japan. Nucleotide sequences of one mitochondrial DNA and eight nuclear gene loci were examined, adding sequences of hares in China from databases. All nuclear DNA analyses supported the clear separation of three phylogroups: L. timidus, L. brachyurus, and the L. mandshuricus complex containing L. coreanus. On the other hand, massive mitochondrial introgression from two L. timidus lineages to the L. mandshuricus complex was suggested in continental East Asia. The northern population of the L. mandshuricus complex was mainly associated with introgression from the continental lineage of L. timidus, possibly since the last glacial period, whereas the southern population of the L. mandshuricus complex experienced introgression from another L. timidus lineage related to the Hokkaido population, possibly before the last glacial period. In contrast to continental hares, no evidence of introgression was found in L. brachyurus in the Japanese Archipelago, which showed the oldest divergence amongst East Asian hare lineages. Our findings suggest that glacial-interglacial climate changes in the circum-Japan Sea region promoted distribution shifts and introgressive hybridization among continental hare species, while the geographic structure of the region contributed to long-term isolation of hares on the islands, preventing inter-species gene flow.

Biogeography of Korea's top predator, the yellow-throated Marten: evolutionary history and population dynamics

BACKGROUND

Peninsulas often harvest high genetic diversity through repeated southward migrations of species during glacial maxima. Studies addressing within-species evolutionary responses to climate fluctuations in northeast Asia are limited compared to other regions of the world, and more so in the Korean Peninsula. In this study, we conducted the first population-level study of the yellow-throated marten, Martes flavigula, from the Korean Peninsula, Russian, Taiwanese and Chinese localities in a biogeographic framework using mitochondrial (cyt-b, nd2, cr) and nuclear gene sequencing (ghr).

RESULTS

Bayesian analyses revealed a rather young population, with a split from the most recent common ancestor at around 125 kya. Martes flavigula likely colonized the Korean Peninsula from Mainland China through the Yellow Sea twice, ca. 60 kya and 20 kya. Korean martens diversified during the Late Pleistocene with at least two dispersal events out of Korea, towards the southwest to Taiwan (ca. 80 kya) and towards the North into Russia and eastern China; most likely after the Last Glacial Maxima (ca. 20 kya). We argue that the lack of population structure and mixed populations is possibly a consequence of the high dispersal capability of the species. The Bayesian skyline plot revealed a population decline within the last 5000 years, suggesting potential negative biotic and anthropogenic effects in the area. We find that local populations are not genetically differentiated, therefore no perceptible population structure within Korea was found.

CONCLUSIONS

The topography and geography of the Korean Peninsula has played a pivotal role in its colonization. Connections between the Korean Peninsula and the Mainland through sea-level drops of the Yellow Sea at times of glacial maxima and the high dispersal capability of M. flavigula adds to the lack of geographical structure in this species and the paraphyly of Korean lineages.

Identification of novel haplotypes and interpretation of gene flow of mitochondrial DNA control region of Eurasian otter (Lutra lutra) for the effective conservation

The number and distribution of Eurasian otters have declined during twentieth century due to human activity and water pollution. The global conservation status of Eurasian otter is presently ‘Near Threatened (NT)’ and strictly protected by being listed on the international legislation and conventions. A number of studies using the mitochondrial DNA (mtDNA) control region (CR) have been conducted in order to effectively apply conservation and reintroduction programs, especially in Europe. However, aside from Europe, there have been few studies concerning genetic diversity and phylogeny of Eurasian otters. Therefore, in this study, we sequenced partial mtDNA CR sequences (232 bp) from five South Korean Eurasian otters and analyzed 27 otters originating from parts of northeast Asia (South Korea, China, Japan and Russia (Sakhalin)), and Europe. Out of 232 bp partial mtDNA CR sequences, 13 polymorphic sites (5.6%) were identified and 4 novel mtDNA CR haplotypes (Lut16–19) were discovered from 12 Eurasian otters originating from northeast Asian region. In this study, a comprehensive analysis of genetic diversity and population structure of Eurasian otter between Europe and northeast Asia continents were conducted. Of these, different past demographic histories in Pleistocene period might have largely impacted the genetic structure of each population differently. In addition, low degree of gene flow, isolation by distance (IBD) pattern from geographically wide distanced dataset and analysis of molecular variance (AMOVA) also represented distinct genetic characteristics of Eurasian otter between Europe and northeast Asia.

Mysteries of host switching: Diversification and host specificity in rodent-coccidia associations

Recent studies show that host switching is much more frequent than originally believed and constitutes an important driver in evolution of host-parasite associations. However, its frequency and ecological mechanisms at the population level have been rarely investigated. We address this issue by analyzing phylogeny and population genetics of an extensive sample, from a broad geographic area, for commonly occurring parasites of the genus Eimeria within the abundant rodent genera Apodemus, Microtus and Myodes, using two molecular markers. At the most basal level, we demonstrate polyphyletic arrangement, i.e. multiple origin, of the rodent-specific clusters within the Eimeria phylogeny, and strong genetic/phylogenetic structure within these lineages determined at least partially by specificities to different host groups. However, a novel and the most important observation is a repeated occurrence of host switches among closely related genetic lineages which may become rapidly fixed. Within the studied model, this phenomenon applies particularly to the switches between the eimerians from Apodemus flavicollis/Apodemus sylvaticus and Apodemus agrarius groups. We show that genetic differentiation and isolation between A. flavicollis/A. sylvaticus and A. agrarius faunas is a secondary recent event and does not reflect host-parasite coevolutionary history. Rather, it provides an example of rapid ecology-based differentiation in the parasite population.

Genetic diversity of Bartonella genotypes found in the striped field mouse (Apodemus agrarius) in Central Europe

We investigated the diversity of Bartonella in Apodemus agrarius, an important rodent of peri-domestic habitats, which has spread into Europe in the past 1000 years. Spleen samples of 344 A. agrarius from Eastern Slovakia were screened for the presence of Bartonella spp. using 16S-23S rRNA internal transcribed spacer region and bacteria were detected in 9% of rodents. Based on sequencing of three housekeeping genes (gltA, rpoB and groEL) Bartonella genotypes were ascribed to the species typical for mice and voles: B. grahamii, B. taylorii and B. birtlesii. However, the study also confirmed presence of genotypes belonging to the B. clarridgeiae/B. rochalimae clade, and the B. elizabethae/B. tribocorum clade, which are not commonly found in woodland rodents. In addition, a potential recombination event between these two genotypes was noted, which highlights an important role of A. agrarius in shaping Bartonella diversity and evolution.

Continental Refugium in the Mongolian Plateau during Quaternary Glacial Oscillations: Phylogeography and Niche Modelling of the Endemic Desert Hamster, Phodopus roborovskii

The Mongolian Plateau (MP), which is situated in the interior of Asia and possesses a typical continental climate, experienced harsh climatic conditions during the Quaternary glacial fluctuations. Although these events likely had huge impacts on the local animal populations, the current effects have hardly been explored. To investigate whether the MP supported a refugium along an oceanic-continental gradient (ROCG), and whether this refugium was glacial or interglacial, we investigated the demographic and phylogeographic history of an endemic mammal species, the desert hamster Phodopus roborovskii. We reconstructed the demographic variation, the phylogeographic diffusion, and modelled the potential habitat during historical periods. The genetic diversity in the MP was the highest among all the localities, and the MP was a suitable habitat throughout the modelled historical periods. A phylogeographic diffusion analysis emphasized the importance of the MP as the centre of origin, preservation and spread for P. roborovskii. The homogeneous landscape provided the opportunity for a wide gene flow, which resulted in low resolution of the phylogenetic relationships. Moreover, P. roborovskii was favoured by the interglacial condition, with both its demographical and geographical ranges expanded within the interglacial periods. The range variation from the Last Glacial Maximum to the current condition reflects a distinct longitudinal shift, while both ranges largely contracted from that of the Last Interglacial. Our results support that the MP served as a refugium and spread centre for P. roborovskii during the Quaternary climate fluctuations. The interglacial expansion and the longitudinal shifts highlighted the important effects of precipitations on the distribution range of species adapted to arid and semi-arid during glacial oscillations.

Genetic diversity and genetic structure of the striped field mouse Apodemus agrarius coreae (Muridae, Rodentia) in Korea

The aim of this study was to investigate the genetic diversity and genetic structure of the striped field mouse Apodemus agrarius coreae in Korea. The Korean A. a. coreae is characterized by high levels of haplotype diversity (Hd=0.967) and low levels of nucleotide diversity (π=0.00683). Haplogroup 1 is well separated from the haplotypes of the neighboring regions of the Korean Peninsula, while the other haplogroups are closely related to those from the Russian Far East. Thus, further investigations are required to confirm the validity of the subspecies status of A. a. coreae by implementing additional morphological characters as well as genetic data from the populations present in the Korean Peninsula and its neighboring countries. Haplogroup 1 includes most Korean haplotypes and forms a star-like haplotype network structure, which reveals relatively low levels of sequence divergence and high frequency of unique mutations (only few mutations are shared in most of the haplotype nodes). The results indicate that the haplotypes of Haplogroup 1 might have experienced population expansion since their migration into Korea, which was further corroborated with negative results of neutrality tests for Korean population of A. a. coreae.

Ancient colonization and within-island vicariance revealed by mitochondrial DNA phylogeography of the mountain hare (Lepus timidus) in Hokkaido, Japan

We examined the phylogenetic status and history of the mountain hare Lepus timidus in and around Hokkaido using mitochondrial cytochrome b (cyt b) sequences from 158 samples from Hokkaido and 14 from Sakhalin, as well as four samples from the Korean hare, L. coreanus. The phylogenetic analysis of the cyt b sequences generated in this study and obtained from DNA databases showed the clear genetic specificity of the Hokkaido lineage as a clade. The Hokkaido lineage was estimated to have diverged from the other conspecific and L. coreanus lineages 0.46 and 0.30 million years ago (Mya), respectively. These results suggest that the common ancestor of the mitochondrial lineage in Hokkaido and Korea inhabited Far East Asia before colonization by the present continental lineages of L. timidus, including the Sakhalin population. We estimated the time of the most recent common ancestor of the Hokkaido population to be 0.17 Mya, and found two distinct haplogroups within the island. One group had greater genetic diversity (mean number of pairwise differences: π = 0.0188 ± 0.0108) and appears to have expanded from the west to the entire island of Hokkaido. The other had lower genetic diversity (π = 0.0038 ± 0.0037) and its distribution was concentrated in the east. These contrasting west/east trends indicate that the Hokkaido population was fragmented in the past, and then subsequently expanded. Our study suggests that Hokkaido was an important refugium for boreal species in the far eastern region, and allowed the formation of various population genetic structures within the island.

Phylogeography of the South China field mouse (Apodemus draco) on the southeastern Tibetan Plateau reveals high genetic diversity and glacial refugia

The southeastern margin of the Tibetan Plateau (SEMTP) is a particularly interesting region due to its topographic complexity and unique geologic history, but phylogeographic studies that focus on this region are rare. In this study, we investigated the phylogeography of the South China field mouse, Apodemus draco, in order to assess the role of geologic and climatic events on the Tibetan Plateau in shaping its genetic structure. We sequenced mitochondrial cytochrome b (cyt b) sequences in 103 individuals from 47 sampling sites. In addition, 23 cyt b sequences were collected from GenBank for analyses. Phylogenetic, demographic and landscape genetic methods were conducted. Seventy-six cyt b haplotypes were found and the genetic diversity was extremely high (π = 0.0368; h = 0.989). Five major evolutionary clades, based on geographic locations, were identified. Demographic analyses implied subclade 1A and subclade 1B experienced population expansions at about 0.052-0.013 Mya and 0.014-0.004 Mya, respectively. The divergence time analysis showed that the split between clade 1 and clade 2 occurred 0.26 Mya, which fell into the extensive glacial period (EGP, 0.5-0.17 Mya). The divergence times of other main clades (2.20-0.55 Mya) were congruent with the periods of the Qingzang Movement (3.6-1.7 Mya) and the Kun-Huang Movement (1.2-0.6 Mya), which were known as the most intense uplift events in the Tibetan Plateau. Our study supported the hypothesis that the SEMTP was a large late Pleistocene refugium, and further inferred that the Gongga Mountain Region and Hongya County were glacial refugia for A. draco in clade 1. We hypothesize that the evolutionary history of A. draco in the SEMTP primarily occurred in two stages. First, an initial divergence would have been shaped by uplift events of the Tibetan Plateau. Then, major glaciations in the Pleistocene added complexity to its demographic history and genetic structure.

The complete mitochondrial genome of the striped field mouse, Apodemus agrarius (Rodentia, Murinae) from Korea

The striped field mouse, Apodemus agrarius, is the most common rodents in temperate forests of the Korean peninsula. We have determined the complete mitochondrial genome (JN629047) of A. agrarius from Korea. Apodemus agrarius mitogenome is 16,263 bp in length and has a base composition of 33.9% A, 29.9% T, 23.6% C, and 12.5% G. The total length of the 13 protein-coding genes is 11,398 bp long.

Phylogeographic study of Apodemus ilex (Rodentia: Muridae) in Southwest China

BACKGROUND

The Mountains of southwest China have complex river systems and a profoundly complex topography and are among the most important biodiversity hotspots in the world. However, only a few studies have shed light on how the mountains and river valleys promote genetic diversity. Apodemus ilex is a fine model for investigating this subject.

METHODOLOGY/PRINCIPAL FINDINGS

To assess the genetic diversity and biogeographic patterns of Apodemus ilex, the complete cytochrome b gene sequences (1,140 bp) were determined from 203 samples of A. draco/ilex that were collected from southwest China. The results obtained suggested that A. ilex and A. draco are sistergroups and diverged from each other approximately 2.25 million years ago. A. ilex could be divided into Eastern and Western phylogroups, each containing two sub-groups and being widespread in different geographical regions of the southern Hengduan Mountains and the western Yunnan - Guizhou Plateau. The population expansions of A. ilex were roughly from 0.089 Mya to 0.023 Mya.

CONCLUSIONS

Our result suggested that A. ilex is a valid species rather than synonym of A. draco. As a middle-high elevation inhabitant, the phylogenetic pattern of A. ilex was strongly related to the complex geographical structures in southwest China, particularly the existence of deep river valley systems, such as the Mekong and Salween rivers. Also, it appears that the evolutionary history of A. ilex, such as lineage divergences and population expansions were strongly affected by climate fluctuation in the Late Pleistocene.

A mitogenome of the Chevrier's field mouse (Apodemus chevrieri) and genetic variations inferred from the cytochrome b gene

The Chevrier's field mouse (Apodemus chevrieri) is an endemic species to China and is an important pest in agriculture and human diseases. In this study, the complete mitochondrial genome of this species was sequenced and its size was 16,298 bases (accession no.: HQ896683). The mitogenome structure was similar compared with other reported rodent mitochondrial genomes and includes 13 protein-coding genes, 2 rRNA genes (12S rRNA and 16S rRNA), 22 tRNA genes, and 1 control region. This was the first complete mitogenome sequenced in genus Apodemus. The phylogenetic analyses based on the sequences of 12 heavy-strand protein-coding genes demonstrated that A. chevrieri clustered together with genus Mus. Additionally, extremely high haplotype and nucleotide diversities (h=0.978, π=2.6%) were observed based on 44 mitochondrial cytochrome b (cyt b) gene sequences. This suggests adaptive divergence of this species to a variety of living habitats and potential refuges in the eastern margin of the Hengduan Mountains during the Quaternary ice ages. No population expansions or genetic bottlenecks were observed in demographic analyses. The phylogenetic analysis of cyt b sequences and haplotypes revealed a genetic differentiation between north and south populations. The divergence between north clade and south clade occurred probably in the middle Pleistocene 1.1815 million years ago (Mya) (95% highest posterior density 2.3189-0.2737 Mya), which was congruent with the periods of the most tense uplift events in the Tibetan Plateau.