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

Abstract

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.

About the Authors:

Qi Liu

Contributed equally to this work with: Qi Liu, Peng Chen, Kai He

Affiliations State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, China, Graduate University of Chinese Academy of Sciences, Beijing, China

Peng Chen

Contributed equally to this work with: Qi Liu, Peng Chen, Kai He

Affiliations State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, China, Graduate University of Chinese Academy of Sciences, Beijing, China

Kai He

Contributed equally to this work with: Qi Liu, Peng Chen, Kai He

Affiliations State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, China, Graduate University of Chinese Academy of Sciences, Beijing, China

C. William Kilpatrick

Affiliation: Department of Biology, University of Vermont, Burlington, Vermont, United States of America

Shao-Ying Liu

Affiliation: Sichuan Academy of Forestry, Chengdu, China

Fa-Hong Yu

* E-mail: [email protected] (FHY); [email protected] (XLJ)

Affiliation: Interdisciplinary Center for Biotechnology Research, University of Florida, Gainesville, Florida, United States of America

Xue-Long Jiang

* E-mail: [email protected] (FHY); [email protected] (XLJ)

Affiliation: State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, China

Introduction

Apodemus species are among the most common small rodents inhabiting woodlands and forests of the Palaearctic and Oriental Region [1], [2], [3], [4]. The genus has been subdivided into four subgenera, Apodemus, Sylvaemus, Alsomys and Karstomys [5] and comprises 20–22 extant species [2], [4], [6].

There are 12 extant Apodemus species in the Oriental Region [2], but only four forms are reported from the Eastern Trans-Himalayas [4], including A. peninsulae, A. latronum, A. chevrieri and the A. draco complex. The Apodemus draco complex include A. draco, A. ilex, and A. orestes and are distributed in mountain areas in China, Myanmar, and India [2]. These taxa have been treated as three valid species, subspecies of a single species, or synonyms of Apodemus draco in different taxonomic revisions [2], [3], [5], [7], [8]. The fossil records and phylogenetic analysis suggested an initial radiation of Apodemus in East Asia into a Japanese endemic (A. argenteus), a Nepalese endemic (A. gurkha) and the ancestral lineage of the remaining eastern Asian species (subgenus Apodemus) after a two-step radiation process associated with the recent tectonic movements that occurred 5–7 Mya (million years ago) and 2–3 Mya, respectively [9], [10], [11]. The phylogeography of Apodemus in the Far East of Asia showed extensive isolations within Apodemus species and could be linked to the presence of many biogeographic barriers such as mountains, rivers, seas, and deserts [1], similar to many other species [9], [12], [13], [14].

The Hengduan Mountains (i.e. the mountains of Southwest China) have the most complex river systems and a profoundly complex and dynamic geological history. The uplifting of the Himalayas and the Qinghai-Tibet Plateau and the successive alternation of glaciation and interglaciation in the Pliocene-Pleistocene contributed to the formation of natural geographical barriers and habitat heterogeneity [15]. This made it not only an important center of relic survival but also decisive evolutionary localities exist. Hence it is an excellent model system for biogeographic studies [16]. Previous studies of Apodemus considered this region to be the Pleistocene refugium or the radiation center for the East Asian Apodemus species [1], [8], [17], [18]. However, there has not been sufficient evidence from morphometric and molecular studies to describe the phylogenetic relationship between A. draco and A. ilex, especially the biogeography in the south of the Hengduan Mountains and the Yunnan-Guizhou Plateau.

In this study, the complete cytochrome b gene sequences (cyt-b) were determined from 203 samples of A. draco/ilex that were collected from southwest China, including the Hengduan Mountains and the Yunnan-Guizhou Plateau. Using phylogenetic and phylogeographical approaches, we examined the different revisions regarding the taxonomic status of A. ilex, the effect of the complex geological structures in the Himalayan regions on the genetic diversity of A. ilex, as well as the hypotheses of the biogeographic patterns and colonization history of A. ilex. In addition, a Bayesian method with a “relaxed” clock model [19] was applied to co-estimate the phylogenetic relationships and divergence times of Apodemus.

Methods

Ethics Statement

All animal samples were obtained following the regulations for the implementation of China on the protection of terrestrial wild animals (State Council Decree [1992] No. 13) and approved by Wildlife Protection Office, Yunnan and Sichuan Provincial Forestry Departments, China as well as the Ethics Committee of Kunming Institute of Zoology, Chinese Academy of Sciences, China.

Sampling and sequencing

A total of 203 samples of Apodemus draco/ilex were collected from 51 localities in the southern Hengduan-Mountains and the Yunnan-Guizhou Plateau in China (Table S1, Figure 1). Specimens were identified based on the description of Thomas [20] and Barrett-Hamilton [21].

[Figure omitted. See PDF.]

Figure 1. Samples of A. ilex used in this study.

Numbers are corresponding to those in Table 1 and presented as pie-charts. Slice size proportional to the frequency of the subclades occurring in the site.

https://doi.org/10.1371/journal.pone.0031453.g001

Total genomic DNAs were extracted from muscle or liver preserved in 95% ethanol at −20°C by using the phenol/proteinase K/sodium dodecyl sulphate method [22]. Mitochondrial cyt-b sequences (1,140 bp) were amplified with the universal primers of L14724 and H15915 [23]. The 50 µl PCR reaction contained 5 µl of 10X PCR buffer, 2 µl of 2 mM dNTP mixture, 2 µl of bovine serum albumin (1 mg/ml), 1 µl of 10 mM of each primer, 2.5 µl of 25 mM MgCl₂, 1.25 U rTaq DNA Polymerase (Takara, Dalian, China) and approximately 100 ng total genomic DNA as template, and DNase/RNase free water diluted to a final volume of 50 µl. A touchdown PCR protocol [24] was used to prevent non-specific amplification, including an initial denaturation at 94°C for 10 min, 10 cycles of denaturation at 94°C for 60 s, annealing at 52.5°C but decreased by 0.5°C per cycle for 60 s, extension at 72°C for 60 s, followed by additional 20 cycles of denaturation at 92°C for 60 s, annealing at 47.5°C for 60 s, extension at 72°C for 60 s, and a final extension at 72°C for 10 min. PCR products were purified with UNIQ-10 spin column DNA gel extraction kit (Shengong, Shanghai, China) and sequenced from both directions with the same PCR primers in an automated DNA sequencer (ABI PRISM 3730) by using Big Dye terminator v3.1 in Tiangen Biotech CO., LTD (Beijing, China).

In addition, 34 cyt-b sequences from GenBank, including 4 sequences of A. ilex, 7 sequences of A. draco, and 23 sequences of A. agrarius, A. alpicola, A. chevrieri, A. flavicollis, A. latronum, A. mystacinus, A. peninsulae, A. semotus, A. speciosus, A. sylvaticus A. uralensis, Mus musculus, and Rattus rattus were included in analyses (Table S2).

Phylogenetic and molecular divergence analysis

The DNA sequences were edited with Seqman and EditSeq (DNASTAR, Lasergene v7.1) and aligned with ClustalX v1.83 [25]. Genetic distance was calculated with MEGA v4.0 [26] with the Kimura 2-parameter (K2P) model [27].

We apply Bayesian inference (BI) and maximum likelihood (ML) to reconstruct phylogenetic relationships. We used RAxML v7.2.8 [28], [29] for ML analyses on the CIPRES Science Gateway v3.1 ([30] http://www.phylo.org). The data set was partitioned by codon. We used the GTRGAMMA model for each partition as recommended by RAxML and selected the novel rapid Bootstrapping algorithm [29] and ran 450 bootstrap replicates.

We used BEAST v1.5.4 for simultaneous Bayesian phylogenetics analysis and “relaxed” molecular dating estimation [31]. The data set was partitioned by codon and the GTR+I+Γ model was selected as the best evolutionary substitution model by likelihood-ratio test in MrModeltest v2.3 [32], [33]. We chose *BEAST model as our tree prior in BEAST analyses [34]. The analyses consisted of a random generated starting tree, uncorrelated lognormal relaxed molecular clock, and the program's default prior distributions of model parameters. Each analysis consisted of 40 million generations, sampled every 1,000 generations. The analyses were repeated four times and convergence was assessed using Tracer v1.5 [31]. Posterior probabilities (PP)≥0.95 are considered to be strongly or significantly supported [35].

Three divergence dates were used as the calibration points in molecular dating estimation and were treated as lognormal distributions [36] in analysis under a relaxed molecular clock model. This is the most appropriate method to use paleontological information [36] and overcomes the problem of rate variation across different timescale [37], [38]. The first calibration date was derived from A. dominans, the common ancestor of A. flavicollis and A. sylvaticus lived about 4.2–4.9 million years ago in European Mammal Neogene (MN) unit 14 [39]. In analysis, the mean sampled age was 4.2 Mya and the older 95% credible interval (CI) was 4.9 Mya (standard deviation = 0.1). The second calibration date was based on the oldest fossil of Apodemus that lived about 11.0–9.88 Mya in MN11 [40]. Estimation of divergence time was calculated with 9.88 Mya (offset = 9.0, mean = 0.88) as the mean sampled age and 11.0 Mya (standard deviation = 0.62) as the older 95% CI. The third calibration date was 12.3–11.0 Mya, the divergence time of Mus and Rattus [41]. The mean sampled age was 11.0 Mya (offset = 10.0, mean = 1.0) and the older 95% CI was 12.3 Mya (standard deviation = 0.63).

Genetic diversity, structure, and population dynamics

The genetic distances among different evolutionary lineages and transversion/transition (Tv/Ts) ratios were estimated with MEGA v4.0 [26]. The number of haplotypes (N), nucleotide diversity (pi), and haplotype diversity (Hd) were calculated with DnaSP v5.10 [42]. The minimum-spanning networks were derived from Network v4.5 using the median-joining approach [43]. Because the Tv/Ts ratio was 5.33 in A. ilex, we set the weight of transversion to 5 and the weight of transition to 1 in analysis. Optimal phylogenetic trees were further calculated with the calculation options of maximum parsimony [44], [45].

Population subdivision was estimated using the hierarchical analysis of molecular variance (AMOVA) and Genetic differences among populations were calculated by pairwise Fst test using Arlequin v3.5 [46]. AMOVA were performed with 10,000 permutations with 3 different grouping options, which were grouped based either on the mtDNA clades/subclades identified in BEAST analysis or on their geographical distributions (see result). Pairwise Fst values were also performed with 10,000 permutations with a significant level of 0.05. An Fst value≥0.25 indicated that the gene flow was limited between two populations [47]. The historical population dynamics were analyzed by mismatch distribution analyses (MDA) [48], Fu's test [49] and Ramos-Onsins and Rozas's R2 test [50]. The behavior of Fs is better for large population sizes, whereas R2 is better for small sample sizes [50]. MDA and the Fu's test was performed using Arlequin v3.5 [46] with 10,000 permutations on each subclade or subgroup if neutrality holds statistical significance. The raggedness index (rg) and sum of squared deviations (SSD) between observed and expected mismatch distributions were estimated simultaneously with MDA. Under population expansion model, the rg and SSD were expected to have lower value [51]. R2 test was also performed using DnaSP v5.10 [42] and significant level was estimated by coalescent simulations with 10,000 replicates. If sudden expansion model was not rejected, the expansion time was calculated with the equation τ = 2ut [48], where t was measured in generations (one year per generation for A. ilex), τ was calculated simultaneously with MDA, u was calculated using formula u = 2 µk, in which 2 µ was the mutation rate per nucleotide and k was the number of nucleotides (1,140 bp).

Results

Phylogenetic analysis and molecular divergence time

Because A. draco and A. ilex are morphologically indistinguishable, the sequences of the samples were identified based on the pairwise comparison with the sequences from the specimens collected at or near their type localities. The type locality for A. draco is located at Kuatum, Fujian, China, while the type locality for A. ilex is at Salween - Mekong divide (28°20′N) [20]. In analysis, the sequences determined from the topotype specimen of A. draco by Liu et al. [8] (Accession number : AY389009) and the specimen from near the type locality of A. ilex from Mt. Meili, China (28°23.8′N) were serve as the reference sequences for A. draco and A. ilex, respectively. Of the 203 sequences generated in this study, 6 were identified as A. draco and 197 as A. ilex (Tables S1, S2). Haplotype analysis of 201 cyt-b sequences of A. ilex, including 4 sequences from GenBank, identified 134 haplotypes. The new identified haplotypes were submitted to GenBank (Accession numbers: JF503102–JF503107 (A. draco) and JF503109–JF503198, JF503200–JF503228, JF503230–JF503240 (A. ilex)).

The phylogeny estimated by RAxML and BEAST were congruent with each other and the topologies were overall highly supported. Thus, only the Bayesian trees were given and both Bayesian posterior probabilities and ML bootstrap support values (BS) were represented (Figures 2, 3). All populations of A. ilex and A. draco formed strongly supported (BS≥94, PP = 1.0) reciprocal monophyletic groups. The sister relationship between A. draco and A. ilex was also supported (BS = 90, PP = 0.87), with 9.0% of a K2P distance. All 201 A. ilex samples were grouped into the Eastern (E) and Western (W) clades, each containing two subclades: E1/E2 and W1/W2 (Figures 2, 3). All clades and subclades were significantly supported by BI analyses (PP = 1.0) and at least moderately supported by ML analyses (BS≥62). The K2P distances between clades and subclades were: E/W = 3.1%, E1/E2 = 1.9% and W1/W2 = 1.6%.

[Figure omitted. See PDF.]

Figure 2. Chronogram of Apodemus based on cyt-b sequences.

Branch lengths represent time; Node bars indicate the 95% CI for the clade age; An asterisk indicates node for calibration; The t_x above the nodes refer to median ages and 95% CI for each node in Table 2; Numbers below the nodes are Bayesian posterior probabilities (PP) and ML bootstrap (BS) values. A ‘-’ indicates the value is lower than 0.5 (PP) or 50 (BS).

https://doi.org/10.1371/journal.pone.0031453.g002

[Figure omitted. See PDF.]

Figure 3. Bayesian phylogenetic analysis of A. draco and A. ilex based on cyt-b sequences.

Branch lengths represent substitution per site and numbers at each node represent the Bayesian posterior probabilities and ML bootstrap values. A ‘-’ indicates the value is lower than 0.5 (PP) or 50 (BS).

https://doi.org/10.1371/journal.pone.0031453.g003

Table 1 presents the divergence times based on the Bayesian relaxed molecular dating estimation. A. ilex and A. draco diverged from their common ancestor at about 2.25 Mya (95% CI = 1.69–2.82). The earliest split within A. draco occurred about 1.15 Mya (95% CI = 0.84–1.46), much earlier than the split of the clades E and W of A. ilex at about 0.62 Mya (95% CI = 0.44–0.84). The divergence times of subclades E1/E2 and W1/W2 were at about 0.33 Mya (95% CI = 0.23–0.45) and 0.32 Mya (95% CI = 0.22–0.45), respectively.

[Figure omitted. See PDF.]

Table 1. Divergence information within and between groups of Apodemus.

https://doi.org/10.1371/journal.pone.0031453.t001

Genetic diversity and structure of A. ilex

Phylogenetics analysis of 201 cyt-b sequences of A. ilex detected 965 conserved sites (84.6% of all sites) and 175 variable sites (15.4% of all sites). The K2P distances between haplotypes of A. ilex ranged from 0.0% to 4.1% (average 2.1%). The overall haplotype diversity (Hd) and nucleotide diversity (Pi) were 0.993 and 0.021, respectively. The pairwise Fst estimation among populations ranged widely from 0.00 to 1.00. Most populations are strongly differentiated from each other (Fst>0.25) indicating restricted gene flow. High levels of gene flow are more often observed within geographically close populations (e.g. population 16 and 18, Fst<0.001; Table S3).

In clade E, the haplotypes in both E1 and E2 are widely distributed in the east and some areas of the west of the Mekong River (populations 1, 14, 34, 36, 39 containing E1 haplotypes and population 1 containing E2 haplotypes; Figure 1). In clade W, the haplotypes in W1 were mainly distributed in the west of the Salween River and two localities east of the Mekong River (populations 10 and 26); while the haplotypes in W2 were distributed mainly in the southern part of the Yunnan-Guizhou Plateau and west of the Mekong River. In addition, sympatric distribution of different maternal lineages were observed in several localities (i.e. population 1, 2, 8, 10, 11, 16, 19, 34, 36, 39; Figure 1)

Further geographical structure was examined with AMOVA using three grouping options, including (1) the populations grouped by the subclades E1, E2, W1 and W2; (2) the populations grouped by geographical distributions, namely, the individuals from the west of the Salween River as group 1, the individuals distributed between the Salween and Mekong rivers as group 2, and the individuals from the east of the Mekong River as group 3; and (3) the populations grouped with the same way as (2) except for populations 1 and 14 which were included in group 3 (Figure 1). The results of AMOVA showed significant genetic structures at all hierarchical levels (P<0.001) and the largest proportion of variances were always found among groups (Table 2). In size order, the variances among groups were the second grouping option (43.30%)<the third grouping option (54.12%)<the first grouping option (55.17%), with the corresponding increased Φst values from 0.715, 0.737 to 0.822, respectively.

[Figure omitted. See PDF.]

Table 2. Results of AMOVA based on different grouping options.

https://doi.org/10.1371/journal.pone.0031453.t002

The network analyses generated eleven most parsimony trees that were similar to the gene tree inferred from the BEAST analyses, including four subclades (E1, E2, E3 and E4; Figure 4). The clade E1 has the most complex structure and can be further divided into 4 subgroups (E1a, E1b, E1c, and E1d). E1a consisted of the haplotypes from Caojian (population 14), E1b consisted of the haplotypes from Mt. Haba, Mt. Yulong, Lushi and Mt. Bangma (populations 3 and 6, 10, 24 and 36, respectively), E1c mainly included the haplotypes from the Mt. Wuliang and Mt. Ailao in central Yunnan (populations 15–22), and E1d contained the haplotypes from the northwest of Yunnan (Table S1). Star-like structures were found in E1c, E1d, W1 and W2 (Figure 4), though the original haplotypes were not found in W1 and W2. These structures are evidence of population expansion [52].

[Figure omitted. See PDF.]

Figure 4. The median-joining network of A. ilex based on cyt-b sequences.

The circle size is proportional to the haplotype frequency and the branch length is proportional to the number of mutations.

https://doi.org/10.1371/journal.pone.0031453.g004

Population historical demography of A. ilex

The analysis of MDA suggested that all four groups showed multimodal distributions (Figure S1) but with small SSD and rg values (Table 3). Moreover, Fu's and R2 test showed the large populations have significant negative Fs (i.e. E1, W1, W2) when small populations have significant small R2 (i.e. E2) (Table 3). Further analyses of the subgroups E1c and E1d revealed smooth unimodal mismatch distributions and significant negative Fu's Fs values (Table 3, Figure S1). The insignificant values of P_SSD and P_rg (P>0.05) in E1c, E1d and W2 confirmed the expansions within these group/subgroups; while the significant P_SSD/P_rg (P<0.05) of W1 might be due to insufficient sample size [53]. Group E2 had multimodal distribution and insignificant negative Fu's Fs value, but the significant small R2 value and the insignificant values of P_SSD and P_rg implied the possibility of expansion in congruent with the network topologies. Because the divergence time between group E and W was at about 0.62 Mya, the population expansions of E1c, E1d and W1 were approximately at 0.064, 0.038 and 0.039 Mya, respectively (Table 3).

[Figure omitted. See PDF.]

Table 3. Neutrality test and Mismatch distribution analyses of A. ilex.

https://doi.org/10.1371/journal.pone.0031453.t003

Discussion

Taxonomic implication of Apodemus ilex

The taxonomic status of A. draco, A. ilex and A. orestes have been controversial for a long time. Apodemus ilex was named based on specimens collected from the Salween-Mekong divide (28°20″N), China [20] but was treated as a synonym of A. orestes [7] or A. draco [5], [54]. Musser et al. [6] included orestes within A. draco, but, after comparing the cranial characteristics between A. orestes and A. draco that are distributed in Wuliang Mountain, China, Jiang and Wang considered A. orestes as a valid species [3]. Patterns of genetic variations observed in the complex of A. draco based on cyt-b genes suggested that A. orestes was a subspecies of A. draco, and A. ilex, which is distributed in the Yunnan-Guizhou Plateau, is a valid species [8]. However, Musser et al. [2] still considered A. ilex and A. orestes as synonyms of A. draco.

With inclusion of 214 samples of A. draco/ilex that were widespread in the southern Hengduan Mountains and Yunnan-Guizhou Plateau, two major phylogroups were identified within the A. draco complex, one representing A. draco that consists of the specimens from eastern and western China, including all specimens from the western Sichuan Plateau, and another representing A. ilex that contains the samples mainly from the southern Hengduan Mountains and the western Yunnan Plateau. The average genetic distance between A. draco and A. ilex was 0.09 (K2P). The molecular dating estimation suggested that the divergence between A. draco and A. ilex was at about 2.25 Mya, earlier than the split of A. alpicola and A. flavicollis or A. agrarius and A. chevrieri at 2.20–1.56 Mya (Table 1, Figure 2). These results support the recognition of A. ilex is a valid species under genetic and phylogenetic species concepts [55], [56]. The taxonomic status of A. orestes will be discussed elsewhere (Chen et al. in preparation).

Phylogeographic structure in A. ilex and topography of mountains and rivers

The Hengduan Mountains have long been recognized as a refugial area for animals [57], [58]. Previous analyses either focused on the northern Hengduan Mountain [59], [60], [61] or treated this area as one refugium Only a few studies have addressed the effect of the extremely complex topography of the southern Hengduan Mountains and the Yunnan Plateau [62], [63]. Our research revealed the significant internal genetic structure within the mountains which is relevant to the “microrefugia” [64] or “refugia within refugia” [65], [66]. These concepts are usually used to explain the phylogenetic structure in the refugia such as the Iberian Peninsula or disjunctive populations surviving in isolated microhabitats. The extremely complex topography [67], climate [68] and habitats [69] in the mountains as well as the mid-high elevation distributed pattern of A. ilex could have lead to the geographically isolation of A.ilex among different mountain areas and the subsequently restricted gene flow, which are respond for the strong geographic structure and the high pairwise Fst values [70].

The minimum-spanning network and AMOVA analyses indicated the geographic structure of A. ilex was also shaped by the Mekong and Salween river systems. When the two rivers were setup as the genetic barriers in the AMOVA analyses, the variances existed mainly among populations in different regions (Φ_ST = 0.715) (Table 2). This result is congruent with the recent proposed hypothesis that deep river valleys may have acted as barriers to Apodemus species [1] as well as to other animals [71]. However, another montane mammal, the Yunnan hare (Lepus comus) has no phylogeographic pattern in the same area. Neither the area of low-elevation nor the river systems seem to be the barrier for the hares [72]. On the other hand, the paleo-drainage systems have facilitated dispersal of a frog species [62]. The discrepancy may due to the different habitats, dispersal abilities and colonization histories. The Nanorana yunnanensis is a semi-aquatic anuran living in cold montane streams [62]; the L. comus is much larger than A. ilex and may have colonized this area only recently [72].

Effect of glaciation

The biogeographic histories of the montane inhabitants are usually affected by Pleistocene glacial cycles [73], [74], [75]. Generally, animals respond to climate change in two different ways [76]. First is by changing their distribution. The montane animals usually have larger distribution during glacial periods (but see [77], [78]) when they spread to lowland [79], [80], and the glacial and interglacial climate fluctuation can result in population isolation and reconnection [73], [74], [81], [82], [83], [84]. Second is by adapting to new environments [75], [76], [81]. Apparently, A. ilex occupies the same habitats as its relatives (i.e. A. draco and A. semotus) and hasn't adapted to a new environment, thus it had to shift distribution responding to climate change.

With a Bayesian method under a “relaxed” clock model, A. ilex diverged from A. draco at around 2.25 Mya. Therefore, the ancestor of A. ilex probably expanded southward from the northern Hengduan Mountains during global cooling in the period 2.7–2.5 Mya [85]. After that, divergence of A. draco/ilex may be attributed to the accelerated uplift of the Qinghai-Tibet Plateau and the resulting geomorphic changes of the plateau and the surrounding areas [86] as well as to turnover of vegetation and habitats [87]. The divergence between the eastern and western populations of A. ilex (∼0.62 Mya) was within the Yulong glaciation (0.73–0.5 Mya) [88] and the simultaneous divergences of the subclades E1/E2 (0.33 Mya) and W1/W2 (0.32 Mya) were consistent with the Lijiang glaciation (0.31–0.13 Mya) [89]. Because the calibration points we used are old (4 Mya to 12 Mya), these divergence time should be taken with caution. However, if the divergence time are “true”, they provide evidence that “nunatak refugia” existed in these mountains. Nunataks are refugia in mountain ranges above the glaciers and are snow-free during glaciations [78]. Indeed, there are evidence of valley glaciers in the Hengduan Mountains [89], [90]. Accordingly, A. ilex probably widely colonized the southern Hengduan Mountains and the Yunnan Plateau before the Yulong glaciation. The development of valley glaciers during the Yulong and Lijiang glaciation resulted in geographical isolation but the populations survived in multiple in situ refugia [77], [78] which result in divergences of the clades/subclades.

Sympatric distributions of different lineages are observed in several localities (Figure 1) which may be attributed to population expansion and the resulting secondary contact during climate fluctuations as well as to the complicated geological history of the drainage system [91]. Both the neutral test and MDA analyses provided evidence of expansion in each clade/subclade. The expansions of E1c, E1d and W1were roughly from 0.089 Mya to 0.023 Mya, corresponding to the last glacial period since 0.11 Mya [92]. Clade W1 and W2 are well segregated by the Salween rivers with only a few exception. On the other hand, without major physical barriers in the east of the Mekong River, populations can expand much more easily and colonize new habitats, resulting in secondary contact of E1 and E2.

Conclusion and Perspectives

The Hengduan Mountains are the most important refugial region in China. Other studies have regarded the mountains as a single refugium, our research has, however, revealed significant internal genetic structure which suggests that the “microrefugia” or “refugia within refugia” models are more relevant.

Our finding suggests that both the low-elevation areas and deep river valleys are strong geographic barriers for A.ilex. However, for aquatic animals in this area, the drainage system is more likely to facilitate dispersal rather than prohibit it [62]. Thus, it seems the drainage system did play a role in shaping geographic patterns, but in different ways for different animals. Furthermore, the evolution of the drainage system may have led to a more complex geographic pattern.

Paleoclimatic change has also shaped genetic structure. The glacial-interglacial cycles not only resulted in inter- and intraspecific divergence, but also led to population expansion and secondary contact.

Our study has shed light on the biodiversity in this area. However, because of the complex topography of the mountains, complicated geological history of the drainage system, Pleistocene climate fluctuation and habitat turnover, it is far from clear how the high endemic biodiversity came into existence. It would be necessary to use comparative phylogeographic approaches of animals distributed in different habitats and with different dispersal abilities to examine how the topography, geographic events and climate change have shaped the biodiversity in the mountains of Southwest China [74], [93], [94].

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

MDA and Fu's Fs test for four subclades of A. ilex .

https://doi.org/10.1371/journal.pone.0031453.s001

(TIF)

Table S1.

Sampling information and genetic variability of A. draco/ilex used in this study.

https://doi.org/10.1371/journal.pone.0031453.s002

(DOC)

Table S2.

Information of outgroups used in this study.

https://doi.org/10.1371/journal.pone.0031453.s003

(DOC)

Table S3.

Pairwise FST values for all pair of populations.

https://doi.org/10.1371/journal.pone.0031453.s004

(XLS)

Acknowledgments

We are grateful to Dr. John R. Stewart and an anonymous reviewer for their useful criticism and suggestions for improvement. We also thank Professor Yong-Gang Yao at Kunming Institute of Zoology, Chinese Academy of Sciences (KIZ, CAS) for network analyses, Dong-Ru Zhang and Min-Sheng Peng at KIZ, CAS for demographic analyses, Joseph Orkin at Washington University at St. Louis for improving the writing, and Kai Zhang at University of Oxford for the early modification.

Author Contributions

Conceived and designed the experiments: KH FHY XLJ. Performed the experiments: QL PC KH. Analyzed the data: QL PC KH. Contributed reagents/materials/analysis tools: SYL. Wrote the paper: QL PC KH CWK FHY XLJ.

Citation: Liu Q, Chen P, He K, Kilpatrick CW, Liu S-Y, Yu F-H, et al. (2012) Phylogeographic Study of Apodemus ilex (Rodentia: Muridae) in Southwest China. PLoS ONE 7(2): e31453. https://doi.org/10.1371/journal.pone.0031453

References

1. Sakka H, Quere J, Kartavtseva I, Pavlenko M, Chelomina G, et al. (2010) Comparative phylogeography of four Apodemus species (Mammalia: Rodentia) in the Asian Far East: evidence of Quaternary climatic changes in their genetic structure. Biological Journal of the Linnean Society 100: 797–821.H. SakkaJ. QuereI. KartavtsevaM. PavlenkoG. Chelomina2010Comparative phylogeography of four Apodemus species (Mammalia: Rodentia) in the Asian Far East: evidence of Quaternary climatic changes in their genetic structure.Biological Journal of the Linnean Society100797821

2. Musser GG, Carleton MD, Wilson DE, Reeder DM (2005) Mammal species of the world: a taxonomic and geographic reference.14.GG MusserMD CarletonDE WilsonDM Reeder2005Mammal species of the world: a taxonomic and geographic reference.14

3. Jiang XL, Wang YX (2000) The field mice (Apodemus) in wuliang mountain with a discussion of A. Orestes. Zoological Research 21: 473–477.XL JiangYX Wang2000The field mice (Apodemus) in wuliang mountain with a discussion of A. Orestes.Zoological Research21473477

4. Corbet GB, Hill JE (1992) The Mammals of the Indomalayan Region: a Systematic Review. Oxford University Press, Oxford. GB CorbetJE Hill1992The Mammals of the Indomalayan Region: a Systematic ReviewOxford University Press, Oxford

5. Corbet GB (1978) 314 p.GB Corbet1978The mammals of the Palaearctic region: a taxonomic review: British Museum (Natural History), London. 314 The mammals of the Palaearctic region: a taxonomic review: British Museum (Natural History), London.

6. Musser G, Carleton M, Wilson D, Reeder D (1993) Mammal species of the world: a taxonomic and geographic reference. Family Muridae, 2nd ed. Smithsonian Institution Press, Washington, DC. pp. 501–756.G. MusserM. CarletonD. WilsonD. Reeder1993Mammal species of the world: a taxonomic and geographic reference. Family Muridae, 2nd edSmithsonian Institution Press, Washington, DC501756

7. Allen GM (1938) The mammals of China and Mongolia (Natural History of Central Asia, vol. 9, Part 2). New York: Ann Mus Nat Hist 26: 1–620.GM Allen1938The mammals of China and Mongolia (Natural History of Central Asia, vol. 9, Part 2).New York: Ann Mus Nat Hist261620

8. Liu XM, Wei FW, Li M, Jiang XL, Feng ZJ, et al. (2004) Molecular phylogeny and taxonomy of wood mice (genus Apodemus Kaup, 1829) based on complete mtDNA cytochrome b sequences, with emphasis on Chinese species. Molecular Phylogenetics and Evolution 33: 1–15.XM LiuFW WeiM. LiXL JiangZJ Feng2004Molecular phylogeny and taxonomy of wood mice (genus Apodemus Kaup, 1829) based on complete mtDNA cytochrome b sequences, with emphasis on Chinese species.Molecular Phylogenetics and Evolution33115

9. Michaux JR, Magnanou E, Paradis E, Nieberding C, Libois R (2003) Mitochondrial phylogeography of the Woodmouse (Apodemus sylvaticus) in the Western Palearctic region. Molecular Ecology 12: 685–697.JR MichauxE. MagnanouE. ParadisC. NieberdingR. Libois2003Mitochondrial phylogeography of the Woodmouse (Apodemus sylvaticus) in the Western Palearctic region.Molecular Ecology12685697

10. Suzuki H, Filippucci MG, Chelomina GN, Sato JJ, Serizawa K, et al. (2008) A biogeographic view of Apodemus in Asia and Europe inferred from nuclear and mitochondrial gene sequences. Biochemical Genetics 46: 329–346.H. SuzukiMG FilippucciGN ChelominaJJ SatoK. Serizawa2008A biogeographic view of Apodemus in Asia and Europe inferred from nuclear and mitochondrial gene sequences.Biochemical Genetics46329346

11. Suzuki H, Sato JJ, Tsuchiya K, Luo J, Zhang YP, et al. (2003) Molecular phylogeny of wood mice (Apodemus, Muridae) in East Asia. Biological Journal of the Linnean Society 80: 469–481.H. SuzukiJJ SatoK. TsuchiyaJ. LuoYP Zhang2003Molecular phylogeny of wood mice (Apodemus, Muridae) in East Asia.Biological Journal of the Linnean Society80469481

12. Avise J (2000) Phylogeography: the history and formation of species. Cambridge, MA, USA: Harvard University Press. 447 p.J. Avise2000Phylogeography: the history and formation of species.Cambridge, MA, USAHarvard University Press447

13. Riddle B, Hafner D, Alexander L, Jaeger J (2000) Cryptic vicariance in the historical assembly of a Baja California Peninsular Desert biota. Proceedings of the National Academy of Sciences of the United States of America 97: 14438–14443.B. RiddleD. HafnerL. AlexanderJ. Jaeger2000Cryptic vicariance in the historical assembly of a Baja California Peninsular Desert biota.Proceedings of the National Academy of Sciences of the United States of America971443814443

14. Serizawa K, Suzuki H, Iwasa MA, Tsuchiya K, Pavlenko MV, et al. (2002) A spatial aspect on mitochondrial DNA genealogy in Apodemus peninsulae from East Asia. Biochemical Genetics 40: 149–161.K. SerizawaH. SuzukiMA IwasaK. TsuchiyaMV Pavlenko2002A spatial aspect on mitochondrial DNA genealogy in Apodemus peninsulae from East Asia.Biochemical Genetics40149161

15. Wu FY, Huang BC, Ye K, Fang AM (2008) Collapsed Himalayan-Tibetan orogen and the rising. Tibetan Plateau Acta Petrologica Sinica 24: 1–30.FY WuBC HuangK. YeAM Fang2008Collapsed Himalayan-Tibetan orogen and the rising.Tibetan Plateau Acta Petrologica Sinica24130

16. Axelrod D, Al Shehbaz I, Raven P (1998) History of the modern flora of China.43–55.D. AxelrodI. Al ShehbazP. Raven1998History of the modern flora of China.4355

17. Gu XD, Liu SY, Wang YZ, Wu H (2009) Development and characterization of eleven polymorphic microsatellite loci from South China field mouse (Apodemus draco). Conservation Genetics 10: 1961–1963.XD GuSY LiuYZ WangH. Wu2009Development and characterization of eleven polymorphic microsatellite loci from South China field mouse (Apodemus draco).Conservation Genetics1019611963

18. Xia WP (1984) A study on Chinese Apodemus with a discussion of its relations to Japanese species. Acta Theriologica Sinica 2: 93–98.WP Xia1984A study on Chinese Apodemus with a discussion of its relations to Japanese species.Acta Theriologica Sinica29398

19. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A (2006) Relaxed phylogenetics and dating with confidence. Plos Biology 4: 699–710.AJ DrummondSYW HoMJ PhillipsA. Rambaut2006Relaxed phylogenetics and dating with confidence.Plos Biology4699710

20. Thomas O (1922) On mammals from the Yunnan Highlands. Annals and Magazine of Natural History 10: 391–406.O. Thomas1922On mammals from the Yunnan Highlands.Annals and Magazine of Natural History10391406

21. Barrett-Hamilton GEH (1900) On geographical and individual variation in Mus sylvaticus and its allies387–428.GEH Barrett-Hamilton1900On geographical and individual variation in Mus sylvaticus and its allies387428

22. Sambrook J, Fritsch EF, Maniatis T (1989) Molecular Cloning: A Laboratory Manual, third edition. Cold Spring Harbor Laboratory Press Cold Spring Harbor, New York. J. SambrookEF FritschT. Maniatis1989Molecular Cloning: A Laboratory Manual, third editionCold Spring Harbor Laboratory Press Cold Spring Harbor, New York

23. Irwin DM, Kocher TD, Wilson AC (1991) Evolution of the Cytochrome-B Gene of Mammals. Journal of Molecular Evolution 32: 128–144.DM IrwinTD KocherAC Wilson1991Evolution of the Cytochrome-B Gene of Mammals.Journal of Molecular Evolution32128144

24. Don RH, Cox PT, Wainwright BJ, Baker K, Mattick JS (1991) Touchdown Pcr to Circumvent Spurious Priming during Gene Amplification. Nucleic Acids Research 19: 4008–4008.RH DonPT CoxBJ WainwrightK. BakerJS Mattick1991Touchdown Pcr to Circumvent Spurious Priming during Gene Amplification.Nucleic Acids Research1940084008

25. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research 25: 4876–4882.JD ThompsonTJ GibsonF. PlewniakF. JeanmouginDG Higgins1997The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.Nucleic Acids Research2548764882

26. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Molecular Biology and Evolution 24: 1596–1599.K. TamuraJ. DudleyM. NeiS. Kumar2007MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0.Molecular Biology and Evolution2415961599

27. Kimura M (1980) A Simple Method for Estimating Evolutionary Rates of Base Substitutions through Comparative Studies of Nucleotide-Sequences. Journal of Molecular Evolution 16: 111–120.M. Kimura1980A Simple Method for Estimating Evolutionary Rates of Base Substitutions through Comparative Studies of Nucleotide-Sequences.Journal of Molecular Evolution16111120

28. Stamatakis A (2006) RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690.A. Stamatakis2006RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models.Bioinformatics2226882690

29. Stamatakis A, Hoover P, Rougemont J (2008) A rapid bootstrap algorithm for the RAxML web servers. Systematic Biology 57: 758–771.A. StamatakisP. HooverJ. Rougemont2008A rapid bootstrap algorithm for the RAxML web servers.Systematic Biology57758771

30. Miller MA, Pfeiffer W, Schwartz T (2010) Creating the CIPRES Science Gateway for inference of large phylogenetic trees. IEEE 1–8.MA MillerW. PfeifferT. Schwartz2010Creating the CIPRES Science Gateway for inference of large phylogenetic trees.IEEE18

31. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology 7: 1471–2148.AJ DrummondA. Rambaut2007BEAST: Bayesian evolutionary analysis by sampling trees.BMC Evolutionary Biology714712148

32. Brandley MC, Schmitz A, Reeder TW (2005) Partitioned Bayesian analyses, partition choice, and the phylogenetic relationships of scincid lizards. Systematic Biology 54: 373–390.MC BrandleyA. SchmitzTW Reeder2005Partitioned Bayesian analyses, partition choice, and the phylogenetic relationships of scincid lizards.Systematic Biology54373390

33. Nylander J (2004) MrModeltest v2. Program distributed by the author. Evolutionary Biology Centre, Uppsala University. J. Nylander2004MrModeltest v2. Program distributed by the authorEvolutionary Biology Centre, Uppsala University

34. Heled J, Drummond AJ (2010) Bayesian Inference of Species Trees from Multilocus Data. Molecular Biology and Evolution 27: 570–580.J. HeledAJ Drummond2010Bayesian Inference of Species Trees from Multilocus Data.Molecular Biology and Evolution27570580

35. Huelsenbeck JP, Rannala B (2004) Frequentist properties of Bayesian posterior probabilities of phylogenetic trees under simple and complex substitution models. Systematic Biology 53: 904–913.JP HuelsenbeckB. Rannala2004Frequentist properties of Bayesian posterior probabilities of phylogenetic trees under simple and complex substitution models.Systematic Biology53904913

36. Ho SYM (2007) Calibrating molecular estimates of substitution rates and divergence times in birds. Journal of Avian Biology 38: 409–414.SYM Ho2007Calibrating molecular estimates of substitution rates and divergence times in birds.Journal of Avian Biology38409414

37. Ho SYW, Phillips MJ, Cooper A, Drummond AJ (2005) Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Molecular Biology and Evolution 22: 1561–1568.SYW HoMJ PhillipsA. CooperAJ Drummond2005Time dependency of molecular rate estimates and systematic overestimation of recent divergence times.Molecular Biology and Evolution2215611568

38. Ho SYW, Shapiro B, Phillips MJ, Cooper A, Drummond AJ (2007) Evidence for time dependency of molecular rate estimates. Systematic Biology 56: 515–522.SYW HoB. ShapiroMJ PhillipsA. CooperAJ Drummond2007Evidence for time dependency of molecular rate estimates.Systematic Biology56515522

39. Martín-Suárez E, Oms O, Freudenthal M, Agustí J, Parés JM (1998) Continental Mio-Pliocene transition in the Granada Basin. Lethaia 31: 161–166.E. Martín-SuárezO. OmsM. FreudenthalJ. AgustíJM Parés1998Continental Mio-Pliocene transition in the Granada Basin.Lethaia31161166

40. Nesin V, Nadachowski A (2001) Late Miocene and Pliocene small mammal faunas (Insectivora, Lagomorpha, Rodentia) of Southeastern Europe. Acta Zoologica Cracoviensia 44: 107–135.V. NesinA. Nadachowski2001Late Miocene and Pliocene small mammal faunas (Insectivora, Lagomorpha, Rodentia) of Southeastern Europe.Acta Zoologica Cracoviensia44107135

41. Benton MJ, Donoghue PCJ (2007) Paleontological evidence to date the tree of life. Molecular Biology and Evolution 24: 26–53.MJ BentonPCJ Donoghue2007Paleontological evidence to date the tree of life.Molecular Biology and Evolution242653

42. Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452.P. LibradoJ. Rozas2009DnaSP v5: a software for comprehensive analysis of DNA polymorphism data.Bioinformatics2514511452

43. Bandelt HJ, Forster P, Rohl A (1999) Median-joining networks for inferring intraspecific phylogenies. Molecular Biology and Evolution 16: 37–48.HJ BandeltP. ForsterA. Rohl1999Median-joining networks for inferring intraspecific phylogenies.Molecular Biology and Evolution163748

44. Polzin T, Daneshmand SV (2002) Extending reduction techniques for the Steiner tree problem. Algorithms-Esa 2002, Proceedings 2461: 795–807.T. PolzinSV Daneshmand2002Extending reduction techniques for the Steiner tree problem.Algorithms-Esa 2002, Proceedings2461795807

45. Polzin T, Daneshmand SV (2003) On Steiner trees and minimum spanning trees in hypergraphs. Operations Research Letters 31: 12–20.T. PolzinSV Daneshmand2003On Steiner trees and minimum spanning trees in hypergraphs.Operations Research Letters311220

46. Excoffier L, Laval G, Schneider S (2005) Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evolutionary Bioinformatics 1: 47–50.L. ExcoffierG. LavalS. Schneider2005Arlequin (version 3.0): An integrated software package for population genetics data analysis.Evolutionary Bioinformatics14750

47. Wright S (1978) Evolution and the Genetics of Population. Variability within and among natural populations. University of Chicago Press, Chicago. S. Wright1978Evolution and the Genetics of Population. Variability within and among natural populationsUniversity of Chicago Press, Chicago

48. Rogers A, Harpending H (1992) Population growth makes waves in the distribution of pairwise genetic differences. Molecular Biology and Evolution 9: 552–569.A. RogersH. Harpending1992Population growth makes waves in the distribution of pairwise genetic differences.Molecular Biology and Evolution9552569

49. Fu Y (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147: 915.Y. Fu1997Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection.Genetics147915

50. Ramos-Onsins SE, Rozas J (2002) Statistical properties of new neutrality tests against population growth. Molecular Biology and Evolution 19: 2092–2100.SE Ramos-OnsinsJ. Rozas2002Statistical properties of new neutrality tests against population growth.Molecular Biology and Evolution1920922100

51. Harpending H (1994) Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Human biology; an international record of research 66: 591.H. Harpending1994Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution.Human biology; an international record of research66591

52. Forster P (2004) Ice Ages and the mitochondrial DNA chronology of human dispersals: a review. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences 359: 255.P. Forster2004Ice Ages and the mitochondrial DNA chronology of human dispersals: a review.Philosophical Transactions of the Royal Society of London Series B: Biological Sciences359255

53. Yao YG, Kong QP, Bandelt HJ, Kivisild T, Zhang YP (2002) Phylogeographic differentiation of mitochondrial DNA in Han Chinese. American Journal of Human Genetics 70: 635–651.YG YaoQP KongHJ BandeltT. KivisildYP Zhang2002Phylogeographic differentiation of mitochondrial DNA in Han Chinese.American Journal of Human Genetics70635651

54. Musser GG, Brothers EM, Carleton MD, Hutterer R (1996) Taxonomy and distributional records of Oriental and European Apodemus, with a review of the Apodemus-Sylvaemus problem. Bonner Zoologische Beiträge 46: 143–190.GG MusserEM BrothersMD CarletonR. Hutterer1996Taxonomy and distributional records of Oriental and European Apodemus, with a review of the Apodemus-Sylvaemus problem.Bonner Zoologische Beiträge46143190

55. Baker RJ, Bradley RD (2006) Speciation in mammals and the genetic species concept. Journal of Mammalogy 87: 643–662.RJ BakerRD Bradley2006Speciation in mammals and the genetic species concept.Journal of Mammalogy87643662

56. Nixon KC, Wheeler QD (1990) An Amplification of the Phylogenetic Species Concept. Cladistics-the International Journal of the Willi Hennig Society 6: 211–223.KC NixonQD Wheeler1990An Amplification of the Phylogenetic Species Concept.Cladistics-the International Journal of the Willi Hennig Society6211223

57. Zhang RZ (2002) Geological events and mammalian distribution in China. Acta Zoologica Sinica 48: 141–153.RZ Zhang2002Geological events and mammalian distribution in China.Acta Zoologica Sinica48141153

58. Zhang RZ (1999) Zoogeography of China. Beijing: Science Press. RZ Zhang1999Zoogeography of ChinaBeijingScience Press

59. Chen W, Li Y, Liu Y, Liu S, Yue B (2011) Complex topographic configuration in the Hengduan Mountains shaped the phylogeographic structure of Chinese white-bellied rats. Journal of Zoology 284: 215–223.W. ChenY. LiY. LiuS. LiuB. Yue2011Complex topographic configuration in the Hengduan Mountains shaped the phylogeographic structure of Chinese white-bellied rats.Journal of Zoology284215223

60. Chen W, Liu S, Liu Y, Hao H, Zeng B, et al. (2010) Phylogeography of the large white-bellied rat Niviventer excelsior suggests the influence of Pleistocene Glaciations in the Hengduan Mountains. Zoological Science 27: 487–493.W. ChenS. LiuY. LiuH. HaoB. Zeng2010Phylogeography of the large white-bellied rat Niviventer excelsior suggests the influence of Pleistocene Glaciations in the Hengduan Mountains.Zoological Science27487493

61. Fan Z, Liu S, Liu Y, Zhang X, Yue B (2011) How Quaternary geologic and climatic events in the southeastern margin of the Tibetan Plateau influence the genetic structure of small mammals: inferences from phylogeography of two rodents, Neodon irene and Apodemus latronum. Genetica 139: 339–351.Z. FanS. LiuY. LiuX. ZhangB. Yue2011How Quaternary geologic and climatic events in the southeastern margin of the Tibetan Plateau influence the genetic structure of small mammals: inferences from phylogeography of two rodents, Neodon irene and Apodemus latronum.Genetica139339351

62. Zhang DR, Chen MY, Murphy RW, Che J, Pang JF, et al. (2010) Genealogy and palaeodrainage basins in Yunnan Province: phylogeography of the Yunnan spiny frog, Nanorana yunnanensis (Dicroglossidae). Molecular Ecology 19: 3406–3420.DR ZhangMY ChenRW MurphyJ. CheJF Pang2010Genealogy and palaeodrainage basins in Yunnan Province: phylogeography of the Yunnan spiny frog, Nanorana yunnanensis (Dicroglossidae).Molecular Ecology1934063420

63. Zhang M, Rao D, Yang J, Yu G, Wilkinson JA (2010) Molecular phylogeography and population structure of a mid-elevation montane frog Leptobrachium ailaonicum in a fragmented habitat of southwest China. Molecular Phylogenetics and Evolution 54: 47–58.M. ZhangD. RaoJ. YangG. YuJA Wilkinson2010Molecular phylogeography and population structure of a mid-elevation montane frog Leptobrachium ailaonicum in a fragmented habitat of southwest China.Molecular Phylogenetics and Evolution544758

64. Rull V (2009) Microrefugia. Journal of Biogeography 36: 481–484.V. Rull2009Microrefugia.Journal of Biogeography36481484

65. Grill A, Amori G, Aloise G, Lisi I, Tosi G, et al. (2009) Molecular phylogeography of European Sciurus vulgaris: refuge within refugia? Molecular Ecology 18: 2687–2699.A. GrillG. AmoriG. AloiseI. LisiG. Tosi2009Molecular phylogeography of European Sciurus vulgaris: refuge within refugia?Molecular Ecology1826872699

66. Gomez A, Lunt DH (2007) Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. Phylogeography in southern European refugia: evolutionary perspectives on the origins and conservation of European biodiversity. Dordrecht: Kluwer Academic Publishers. pp. 155–188.A. GomezDH Lunt2007Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. Phylogeography in southern European refugia: evolutionary perspectives on the origins and conservation of European biodiversityDordrechtKluwer Academic Publishers155188

67. Li BY (1989) Geomorphologic regionalization of the Hengduan Mountainous region. Journal of Mountain Research 7: 13–20.BY Li1989Geomorphologic regionalization of the Hengduan Mountainous region.Journal of Mountain Research71320

68. Zhao E, Yang D (1997) Amphibians and reptiles of the Hengduan Mountain region;. In: Sun HL, editor. Beijing: Sicence Press. E. ZhaoD. Yang1997Amphibians and reptiles of the Hengduan Mountain region;HL SunBeijingSicence Press

69. Yao YH, Zhang BP, Han F, Pang Y (2010) Spatial Pattern and Exposure Effect of Altitudinal Belts in the Hengduan Mountains. Journal of Mountain Science 28: 11–20.YH YaoBP ZhangF. HanY. Pang2010Spatial Pattern and Exposure Effect of Altitudinal Belts in the Hengduan Mountains.Journal of Mountain Science281120

70. McCormack J, Huang H, Knowle L (2009) Sky islands. Encyclopedia of Islands 4: 841–843.J. McCormackH. HuangL. Knowle2009Sky islands.Encyclopedia of Islands4841843

71. Song G, Qu Y, Yin Z, Li S, Liu N, et al. (2009) Phylogeography of the Alcippe morrisonia (Aves: Timaliidae): long population history beyond late Pleistocene glaciations. BMC Evolutionary Biology 9: 143.G. SongY. QuZ. YinS. LiN. Liu2009Phylogeography of the Alcippe morrisonia (Aves: Timaliidae): long population history beyond late Pleistocene glaciations.BMC Evolutionary Biology9143

72. Wu CH, Li HP, Wang YX, Zhang YP (2000) Low genetic variation of the Yunnan hare (Lepus comus G. Allen 1927) as revealed by mitochondrial cytochrome b gene sequences. Biochemical Genetics 38: 147–153.CH WuHP LiYX WangYP Zhang2000Low genetic variation of the Yunnan hare (Lepus comus G. Allen 1927) as revealed by mitochondrial cytochrome b gene sequences.Biochemical Genetics38147153

73. Hewitt G (2004) Genetic consequences of climatic oscillations in the Quaternary. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences 359: 183–195.G. Hewitt2004Genetic consequences of climatic oscillations in the Quaternary.Philosophical Transactions of the Royal Society of London Series B: Biological Sciences359183195

74. Hewitt G (2000) The genetic legacy of the Quaternary ice ages. Nature 405: 907–913.G. Hewitt2000The genetic legacy of the Quaternary ice ages.Nature405907913

75. Hewitt G (1996) Some genetic consequences of ice ages, and their role in divergence and speciation. Biological Journal of the Linnean Society (United Kingdom) 58: 247–276.G. Hewitt1996Some genetic consequences of ice ages, and their role in divergence and speciation.Biological Journal of the Linnean Society (United Kingdom)58247276

76. Stewart JR (2009) The evolutionary consequence of the individualistic response to climate change. Journal of Evolutionary Biology 22: 2363–2375.JR Stewart2009The evolutionary consequence of the individualistic response to climate change.Journal of Evolutionary Biology2223632375

77. Bettin O, Cornejo C, Edwards P, Holderegger R (2007) Phylogeography of the high alpine plant Senecio halleri (Asteraceae) in the European Alps: in situ glacial survival with postglacial stepwise dispersal into peripheral areas. Molecular Ecology 16: 2517–2524.O. BettinC. CornejoP. EdwardsR. Holderegger2007Phylogeography of the high alpine plant Senecio halleri (Asteraceae) in the European Alps: in situ glacial survival with postglacial stepwise dispersal into peripheral areas.Molecular Ecology1625172524

78. Holderegger R, Thiel-Egenter C (2009) A discussion of different types of glacial refugia used in mountain biogeography and phylogeography. Journal of Biogeography 36: 476–480.R. HoldereggerC. Thiel-Egenter2009A discussion of different types of glacial refugia used in mountain biogeography and phylogeography.Journal of Biogeography36476480

79. Stewart JR (2003) Comment on “Buffered Tree Population Changes in a Quaternary Refugium: Evolutionary Implications”. Science 299: 825.JR Stewart2003Comment on “Buffered Tree Population Changes in a Quaternary Refugium: Evolutionary Implications”.Science299825

80. Assefa A, Ehrich D, Taberlet P, Nemomissa S, Brochmann C (2007) Pleistocene colonization of afro-alpine ‘sky islands’ by the arctic-alpine Arabis alpina. Heredity 99: 133–142.A. AssefaD. EhrichP. TaberletS. NemomissaC. Brochmann2007Pleistocene colonization of afro-alpine ‘sky islands’ by the arctic-alpine Arabis alpina.Heredity99133142

81. Stewart JR, Lister AM, Barnes I, Dalén L (2010) Refugia revisited: individualistic responses of species in space and time. Proceedings of the Royal Society B: Biological Sciences 277: 661.JR StewartAM ListerI. BarnesL. Dalén2010Refugia revisited: individualistic responses of species in space and time.Proceedings of the Royal Society B: Biological Sciences277661

82. Hewitt G (1999) Post-glacial re-colonization of European biota. Biological Journal of the Linnean Society 68: 87–112.G. Hewitt1999Post-glacial re-colonization of European biota.Biological Journal of the Linnean Society6887112

83. Yuan SL, Lin LK, Oshida T (2006) Phylogeography of the mole-shrew (Anourosorex yamashinai) in Taiwan: implications of interglacial refugia in a high-elevation small mammal. Molecular Ecology 15: 2119–2130.SL YuanLK LinT. Oshida2006Phylogeography of the mole-shrew (Anourosorex yamashinai) in Taiwan: implications of interglacial refugia in a high-elevation small mammal.Molecular Ecology1521192130

84. McCormack JE, Bowen BS, Smith TB (2008) Integrating paleoecology and genetics of bird populations in two sky island archipelagos. BMC Biology 6: 28.JE McCormackBS BowenTB Smith2008Integrating paleoecology and genetics of bird populations in two sky island archipelagos.BMC Biology628

85. Yamada K, Tanaka Y, Irizuki T (2005) Paleoceanographic shifts and global events recorded in late Pliocene shallow marine deposits (2.80–2.55 Ma) of the Sea of Japan. Palaeogeography, Palaeoclimatology, Palaeoecology 220: 255–271.K. YamadaY. TanakaT. Irizuki2005Paleoceanographic shifts and global events recorded in late Pliocene shallow marine deposits (2.80–2.55 Ma) of the Sea of Japan.Palaeogeography, Palaeoclimatology, Palaeoecology220255271

86. An Z, Wang S, Wu X, Chen M, Sun D, et al. (1999) Eolian evidence from the Chinese Loess Plateau: The onset of the Late Cenozoic great glaciation in the Northern Hemisphere and Qinghai-Xizang Plateau uplift forcing. Science in China Series D: Earth Sciences 42: 258–271.Z. AnS. WangX. WuM. ChenD. Sun1999Eolian evidence from the Chinese Loess Plateau: The onset of the Late Cenozoic great glaciation in the Northern Hemisphere and Qinghai-Xizang Plateau uplift forcing.Science in China Series D: Earth Sciences42258271

87. Cerling TE, Harris JM, MacFadden BJ, Leakey MG, Quade J, et al. (1997) Global vegetation change through the Miocene/Pliocene boundary. Nature 389: 153–158.TE CerlingJM HarrisBJ MacFaddenMG LeakeyJ. Quade1997Global vegetation change through the Miocene/Pliocene boundary.Nature389153158

88. Wu Y, Cui Z, Liu G, Ge D, Yin J, et al. (2001) Quaternary geomorphological evolution of the Kunlun Pass area and uplift of the Qinghai-Xizang (Tibet) Plateau. Geomorphology 36: 203–216.Y. WuZ. CuiG. LiuD. GeJ. Yin2001Quaternary geomorphological evolution of the Kunlun Pass area and uplift of the Qinghai-Xizang (Tibet) Plateau.Geomorphology36203216

89. Zhao XT, Zhang YS, Qu Y, Guo CB (2007) Pleistocene glaciations along the western foot of the Yulong Mountains and their relationhip with the formation and development of the Jinsha River. Quaternary sciences 27: 35–44.XT ZhaoYS ZhangY. QuCB Guo2007Pleistocene glaciations along the western foot of the Yulong Mountains and their relationhip with the formation and development of the Jinsha River.Quaternary sciences273544

90. Zhao X, Qu Y, Zhang Y, Hu D, Guo C (2007) Discovery of Shigu paleolake in the Lijiang area, northwestern Yunnan, China and its significance for the development of the modern Jinsha River valley. Geological Bulletin of China 26: 960–969.X. ZhaoY. QuY. ZhangD. HuC. Guo2007Discovery of Shigu paleolake in the Lijiang area, northwestern Yunnan, China and its significance for the development of the modern Jinsha River valley.Geological Bulletin of China26960969

91. Clark M, Schoenbohm L, Royden L, Whipple K, Burchfiel B, et al. (2004) Surface uplift, tectonics, and erosion of eastern Tibet from large-scale drainage patterns. Tectonics 23: 1006–1029.M. ClarkL. SchoenbohmL. RoydenK. WhippleB. Burchfiel2004Surface uplift, tectonics, and erosion of eastern Tibet from large-scale drainage patterns.Tectonics2310061029

92. Thompson LG, Yao T, Davis ME, Henderson KA, MosleyThompson E, et al. (1997) Tropical climate instability: The last glacial cycle from a Qinghai-Tibetan ice core. Science 276: 1821–1825.LG ThompsonT. YaoME DavisKA HendersonE. MosleyThompson1997Tropical climate instability: The last glacial cycle from a Qinghai-Tibetan ice core.Science27618211825

93. Taberlet P, Fumagalli L, Wust-Saucy AG, Cosson JF (1998) Comparative phylogeography and postglacial colonization routes in Europe. Molecular Ecology 7: 453–464.P. TaberletL. FumagalliAG Wust-SaucyJF Cosson1998Comparative phylogeography and postglacial colonization routes in Europe.Molecular Ecology7453464

94. Qu Y, Lei F, Zhang R, Lu X (2010) Comparative phylogeography of five avian species: implications for Pleistocene evolutionary history in the Qinghai-Tibetan plateau. Molecular Ecology 19: 338–351.Y. QuF. LeiR. ZhangX. Lu2010Comparative phylogeography of five avian species: implications for Pleistocene evolutionary history in the Qinghai-Tibetan plateau.Molecular Ecology19338351