The colonization and divergence patterns of Brandt’s vole (Lasiopodomys brandtii) populations reveal evidence of genetic surfing
© The Author(s). 2017
Received: 9 March 2017
Accepted: 7 June 2017
Published: 21 June 2017
The colonial habit of Brandt’s vole (Lasiopodomys brandtii) differs from that of most other species of the genus Microtus. The demographic history of this species and the patterns shaping its current genetic structure remain unknown. Here, we explored patterns of genetic differentiation and infered the demographic history of Brandt’s vole populations through analyses of nuclear microsatellite and D-loop sequences.
Phylogenetic analyses divided the sampled populations into three main clusters, which represent the southeastern, northeastern and western parts of the total range in Mongolia and China. Molecular data revealed an ancestral area located in the southeast of the extant range, in the Xilinguole District, Inner Mongolia, China, from where Brandt’s vole populations began expanding. A gene flow analysis suggested that the most likely colonization route was from the ancestral area and was followed by subsequent northeastward and westward range expansions. We identified decreases in genetic diversity with increasing distance from the founder population within the newly occupied regions (northeastern and western regions), clinal patterns in the allele frequencies, alleles that were rare in the original area that have become common in the newly occupied regions, and higher genetic differentiation in the expanded range compared with the original one.
Our results indicate that L. brandtii most likely originated from the southeastern part of its current geographic range, and subsequently colonized into the northeastern and western parts by expansion. The genetic patterns among the derived populations and with respect to the original population are consistent with that expected under genetic surfing models, which indicated that genetic drift, rather than gene flow, is the predominant factor underlying the genetic structure of expanding Brandt’s vole populations.
The evolutionary history of most species generally is characterized by an interplay among environmental changes, episodes of range expansion and migration, population admixtures, and local extinctions [1, 2]. These and other factors that play roles in the past leave detectable signatures in the genetic structures of modern populations . For a number of Arvicolinae rodents in the Northern Hemisphere such population genetic and phylogeographic studies have been reported [3–8].
Brandt’s vole (Lasiopodomys brandtii) (Radde, 1861) is a steppe-dwelling rodent species that is currently distributed from the central parts of Inner Mongolia, through the central and eastern of Republic of Mongolia, and to the southern borders of Mongolia in Trans-Baikalia, Russia [9, 10]. Brandt’s vole forms its own lineage among arvicolids rodents . Over the past decades, the evolutionary origins of Brandt’s vole have been studied based on paleontological, cytological and nuclear DNA data [10–15]. However, prior to this study, no direct molecular evidence ragarding the geographical origins of the ancestral populations or their colonization and differentiation patterns has been reported.
Brandt’s vole is well adapted to the colonization of patchy, ephemeral habitats [16, 17], where the populations exhibit seasonal and multi-annual fluctuations in abundance [17–19]. The species distributed discontinuously across the Mongolian plateau . However, the habitats are rather homogeneous  with respect to plant community attributes, particularly species composition, vegetation height and degree of plant cover . Brandt’s vole prefers degraded grasslands , which are often interspersed with less favorable habitats characterized by livestock grazing [16, 21, 23]. The populations of Brandt’s vole show declines when the grass is shorter than approximately 5 cm and sparse (less than 40% cover), and do not persist well in habitats with dense coverage (more than 80%) by tall grasses (more than 17 to 20 cm) [17, 21, 24]. The environmental attributes affecting Brandt’s vole populations generally are related to climate and human activities [25, 26].
Climatic changes that alter the availability of suitable habitat have can trigger range shifts and local extinctions of populations for most organisms . Such effects have also been observed on the Mongolian plateau, where Brandt’s vole occurs. There cold and dry climate during the Pleistocenens led to increased coverage by tundra steppe and then dry steppe habitats [28, 29]. The forest-steppe boundary of the plateau was situated further south during the early-mid Holocene stage (10–5 kyr BP), which caused by the enhanced aridity occurred during the period . Currently, it is believed that climate change and overgrazing are driving major factors that result in the degradation and desertification of Mongolian plateau pastureland [30, 31]. The shifts in the distribution of Brandt’s vole populations are expected to result from species migration and/or adaptation to the environment. Biotope evolutionary changes, such as soil development , hydrological changes , desertification , and landscape evolution  in the Mongolia plateau inhabited by extant and extinct Brandt’s voles should impact the population dynamics. Such changes should also interfere with the genetic composition of the population. Through fine-scale geographic sampling in well-defined ecological and historical contexts, it is possible to detect cases of spatially varying selection that involve subtle shifts in the allele frequency among locally adapted populations [35–38].
Here we used microsatellite data and mtDNA D-loop sequences to describe the phylogeography and genetic diversity of Brandt’s vole populations over space and time, and to infer the population demography history. We reconstructed the likely geographic origin of the species and past colonization routes. Genetic surfing is a phenomenon associated with species displaying rapid population growth, a patchy distribution, and rapid population growth [27, 39]. As part of the genetic surfing phenomenon, the frequency of alleles arising on the leading edges of the wave of range expansion can increase to high levels due to genetic drift; in other words, once rare alleles might become predominant in populations colonizing the new territory . The phenomenon has attained some attention because under such a scenario, the relative position of an individual and the alleles it carries more strongly affect the fate of alleles than selection or standard genetic drift occurring in a Fisher-Wright population or stable populations of small effective size . Thus, during our study we paid particular attention to the possibility of empirically demonstrating the phenomenon in Brandt’s vole.
DNA extraction, microsatellite genotyping and D-loop sequencing
Genomic DNA was extracted using the proteinase-K digestion and phenol-chloroform extraction method as described by Sambrook et al. . All individuals were genotyped at 12 microsatellite loci [42, 43]. Primer sets are provided in Additional file 1: Table S2. DNA amplification was performed in a thermocycler (Model 9700, Applied Biosystems, USA). The 15-μL reaction mixtures contained ~50 ng of DNA, 1× PCR buffer, 0.2 μM of each primer, 1.0–2.5 mM MgCl2, 0.2 mM of each nucleotide, 1 U of Taq polymerase (TaKaRa Bio Company, Dalian, China), and ddH2O. The amplification program consisted of an initial denaturation for 10 min at 95 °C followed by 35 cycles of 40 s at 94 °C for denaturation, 40 s for annealing, and 40 s at 72 °C for extension and a final extension for 10 min at 72 °C. The PCR products were analyzed using an ABI PRISM3730 DNA Sequencer (Applied Biosystems). Genotypes were scored using GeneMarker 1.7 (Applied Biosystems). Two observers identified genotypes independently. If the genotypes of the same vole and locus differed, the genotyping process was repeated until a consensus was reached. Allelic dropout and the presence of null alleles were screened using Micro-Checker version 2.2 . Ninety-eight randomly selected samples (12% of all samples) were genotyped again, and the results were analyzed using the Gimlet program  to estimate the rates of different types of errors.
A 734-bp fragment of the D-loop sequence was amplified by PCR from 746 individuals using the primer pairs (5′-ACCATCAACACCCAAAGC-3′ and 3′-GTACTTGATACCCTCTCC-5′) designed with Primer 5  based on the sequence of Lasiopodomys mandarinus with GenBank accession number NC_025283.1. The PCR was performed in 25-μL reaction mixtures containing ~25 ng of genomic DNA, 1× PCR buffer, 0.2 μM of each primer, 1.5 mM MgCl2, 0.2 mM of each nucleotide, 1 U of Taq polymerase (TaKaRa Bio Company, Dalian, China), and ddH2O. The reaction was optimized and programmed over 35 cycles using a GeneAmp PCR system 9700 thermocycler (Applied Biosystems) with the following temperature profile: 5 min at 94 °C for denaturation; followed by 10 cycles of 30 s at 94 °C for denaturation, 30 s at 60 °C for annealing, and 70 s at 72 °C for extension; 15 cycles of 30 s at 94 °C, 30 s at 55 °C, 70 s at 72 °C; and 10 cycles of 94 °C for 30 s, 50 °C for 30 s, 72 °C for 70 s; and a final extension step for 7 min at 72 °C. All PCR products were sequenced in a single direction using an automated ABI 3730 DNA sequencer (Sangon Biotech Ltd). The sequences were aligned to the corresponding D-loop regions of L. mandarinus using the programs Lasergene 7.2  and CLUSTALX  and were checked and edited manually using Lasergene 7.2. The D-loop haplotype sequences from Brandt’s voles were deposited in GenBank (GenBank nos: KY354521-KY354550).
Genetic diversity and phylogenetic analyses
Hardy-Weinberg equilibrium (HWE) tests based on the excess of heterozygosity were performed for all loci using GENEPOP v4.5 . Each test was run for 1000 dememorization steps, followed by 100 batches of 1000 steps each. P-values were adjusted using the Bonferroni correction . We used GENALEX 6.5  to estimate the number of alleles (N a ), allele frequencies, private and rare alleles (allelic frequencies under 0.05), and observed (H o ) and expected heterozygosity (H e ) for each sampling site and for cluster inferred from phylogenetic analyses (see below). Pairwise population differences were also explored with the use of Wright’s F ST values. The allelic richness (A r ) was calculated for each site population based on a minimum sample size of 14 using the rarefaction procedure implemented in HP-Rare .
We used Poptree2  to construct a neighbor-joining (NJ) tree for all Brandt’s vole populations based on the pairwise genetic distances (Da) generated from microsatellite genotypes, and we used a Bayesian MCMC assignment method implemented in STRUCTURE 2.3  to assign individuals to inferred clusters based on multilocus genotypes. The admixture ancestry and the correlated allele frequencies model were implemented in STRUCTURE, with K (the number of clusters) set to a value between 1 and 23, and 20 independent iterations consisting of 100,000 burn-in steps and an iteration length of 100,000 were performed for each value of K. The results subsequently were processed by STRUCTURE HARVESTER (http://taylor0.biology.ucla.edu/structureHarvester/)  to delineate the most likely level of population subdivision (the appropriate K), which was identified using the maximal values of a statistic ΔK based on the rate of change in the log probability of the data between successive K-values . The output from STRUCTURE was visualized using the DISTRUCT program (http://www.stanford.edu/group/rosenberglab/distructDownload.html). We also performed AMOVA analysis with ARLEQUIN  to test for the genetic differentiation within and among clusers.
We estimated the phylogenetic relationships among haplotypes using the Bayesian method implemented in BEAST v2.1.3 . We consider three rodent species to be the most adequate outgroups of L. brandtii: L. mandarinus (GenBank nos. NC_025283.1, JX014233.1 and KF819832.1) always form the sister group with L. brandtii, and it can be rather closely related but clearly outside of L. brandtii in the phylogenetic relationship [3, 60]. Microtus oeconomus (GenBank nos. HM135921, AJ616853 and HM135928) and Microtus arvalis (GenBank no. KP013595) have more distant in the molecular phylogeny of Microtus species [3, 60, 61]. Paleontological and molecular phylogenetic analysis suggested that Lasiopodomys might be part of the sister clade to Microtus [10, 62]. We used the HKY substitution model as selected by Akaike’s Information Criterion  in MODELTEST  and a strict clock model generated by Bayes factors . The Yule model was selected as tree prior model; MCMC chains were run for 50 million iterations with parameters sampled every 1000 iterations; Trees were combined in TreeAnnotator v2.1.3  with 50% burn-in values. Tracer v1.5.0  was used to check for the convergence of Markov chains and to ensure sufficient sampling. Expected sample sizes (ESSs) for all parameters were greater than 100. The maximum clade credibility tree was visualized using Figtree v1.4.2 . We constructed a Median-Joining (MJ) network for D-loop haplotypes using Network5.0 . We performed AMOVA to test for the differentiation within and among cluster and haplogroup using ARLEQUIN . The geographic distributions of the haplotypes were visualized on a map using ArcGIS v10.1 .
Divergence times separating haplogroups were estimated using the Bayesian method implemented in BEAST v2.1.3. To conduct the dating runs, we manipulated a Calibrate Yule tree prior as the tree model and a strict clock model. The divergence time between L.mandarinus and L.brandtii (0.5–0.95 Ma)  was assumed to be the calibration time. MCMC chains were run for 50 million iterations with parameters sampled every 5000 generations and 10% burn-in values. Tracer 1.5 was used to assess ESSs for all paremeters and to verify that the posterior distribution of the divergence calibration matches the prior distribution. Trees were combined in TreeAnnotator v2.1.3 and node ages were visualized using Figtree v1.4.2.
Inferring ancestral areas and modeling colonization routes
We used three clusters (as obtained from D-loop) for ancestral area reconstruction. They were partitioned in the western (W), southeastern (SE) and northeastern (NE) based on the D-loop phylogenetic analysis (Fig. 3c; Table 2). We performed a Statistical-Dispersal-Vicariance-Analysis (S-DIVA) in a Bayesian framework implemented in RASP v3.2 [71, 72], and Dispersal-Extinction-Cladogenesis process (DEC) in a likelihood framework using Lagrange [73, 74] to reconstruct the likely geographic history of Brandt’s vole. We ran S-DIVA on 10,000 trees from MCMC output to account for phylogeny uncertainties. We opted for not constraining the dispersal probabilities to avoid over-parameterization for the DEC analysis. The number of maximum areas was set as 2 for both two analyses. The input tree for the analyses was generated from BEAST (Fig. 2). For Lagrange, a Python script was created using the online Lagrange configurator.
To explore the dispersal process and historical gene flows, we applied MIGRATE-N v.3.6.1 [75, 76] to the microsatellite data for both two distinct clusters (eastern and western distribution) (Additional file 2: Fig. S1) and three clusters (southeastern, northeastern and western distribution) supported by the phylogenetic results (Fig. 1; Additional file 2: Fig. S1; Table 2). MIGRATE-N, implements coalescent-based MCMC simulations, was applied to estimate Θ, M and the marginal likelihood of the specified gene flow models . The population size Θ is defined as 4N e μ, where N e is the effective population size and μ is the mutation rate of loci per generation, and the migration rate M is equivalent to m/μ, where m is the immigration rate per generation . Preliminary runs were performed using the full migration model to determine the optimal parameters. We performed runs through Bayesian inference using the following parameters: slice sampling, an exponential prior distribution for Θ (min: 0, mean: 35, max: 70), an exponential prior distribution for migration (min: 0, mean: 750, max: 1500), static heating with temperatures of 1.00, 1.50, 3.00 and 106; potential occurrence of swapping among chains at every step; and 104 burn-in steps followed 5 × 107 parameter samplings recorded at intervals of 103. We randomly selected 20 individuals per population for our MIGRATE-N analyses due to computational demands and the evidence that more than 20 individuals do not increase the accuracy of parameter estimations . All final reported runs met the convergence criteria for ESSs greater than 104 and showed good agreement between the mean and median estimates for all parameters. The posterior distribution plot of each parameter was also used to visually test for convergence. The marginal likelihood and probability of each model was approximated using a Bezier-corrected thermodynamic integration [76, 78].
Testing for range expansion
We used DNASP 5.0  to conduct neutrality tests [79–81] on D-loop sequence data. Past population size changes were inferred by computing a raggedness index (r) obtained from a site mismatch analysis [82, 83] and by comparing the observed distribution to the distribution under the null model of constant population size.
Effective population size fluctuations over time were inferred using a Bayesian Skyline plot method (BSP) implemented in the program package BEAST v2.1.3. The GTR model was run for the entire population and for each D-loop haplogroup separately. Substitution rate for the D-loop sequences estimated for Meriones was similar to Microtus . We used the substitution rate of 1.27 × 10−7 substitutions per site per year in D-loop estimated for Meriones meridianus (95% confidence interval = 1.2 × 10−7 to 1.33 × 10−7/site/year)  to scale the time axis of BSPs. The chains were run for 200 (entire samples) and 40 (haplogroups) million iterations respectively, of which the first 10% were discarded. Model parameters were sampled every 1000 iterations. The posterior samples were combined using LogCombiner . Skyline plots were generated using Tracer v1.5.0. ESSs for all parameters were greater than 100.
The isolation by distance was tested for all site populations in Brandt’s vole through a regression of pairwise F ST/ (1- F ST) values on the lnGD (geographic distance among populations)  using SPSS v.20 . The pairwise F ST valuas were calculated based on the microsatellite data using ARLEQUIN 3.5 . We also performed IBD tests for different clusters as described above (Fig. 1b).
In a pairwise manner for each of the three clusters, we iteratively calculated the average F ST between samples obtained from one site of each cluster and all samples from all sites that are part of another cluster. According to Graciá et al. , those samples (and the sampling site they represent) that have the lowest average F ST and the highest allelic richness (A r ) was regarded as those representing the likely starting point from which the species began colonizing, and samples (and the sampling site they represent) with the lowest average F ST value was considered as those representing probable newly established derived populations occupying the new territory, or arrival sites. Range expansion is expected to cause a decrease in intrademe heterozygosity with an increase in F ST as a function of the distance from the inferred starting point of colonization [88–91]. We regressed the H e , A r and allele frequencies of each sampling site against the distance to the inferred starting or arrival sites within each cluster as inferred as described above from microsatellite data (Fig. 1b). Due to departures from IBD for the southeastern, northeastern and western clusters (Figs. 1b and 6), we conducted these analyses using the corresponding F ST values instead of the geographic distance, as in such a case F ST might more accurately reflect the dispersal distance between sites. Statistical analyses were performed using SPSS v.20.
Genetic diversity of Brandt’s vole populations
No significant linkage disequilibrium was detected amongst any of the 12 pairwise microsatellite loci combinations. The genotyping results are listed in Additional file 1: Table S1, and characteristics of the microsatellite markers are summarized in Additional file 1: Table S2. The MICRO-CHECKER results suggested no notable scoring errors (the error rate of most loci was less than 0.05) due to stuttering or allele dropout. The number of alleles for each population across loci ranged from 3.83 (HW) to 8.17 (BY), with an average of 6.33. The mean observed heterozygosity (H o ) and expected heterozygosity (H e ) were 0.68 (0.54 for BL to 0.75 for AL) and 0.66 (0.55 for SH to 0.75 for AL), respectively. The mean allelic richness (A r ) was 4.66 (3.39 for SH to 5.61 for AL) (Additional file 1: Table S3). Significant deviations from the Hardy-Weinberg expectation were found for 11 loci in some populations (P < 0.05) (Additional file 1: Table S4).
The analysis of D-loop sequences revealed 30 haplotypes based on 29 segregating sites (including insertions and deletions) observed in the 746 individuals sequenced. Of the 30 haplotypes, 14 were singletons, six were shared among individuals within geographical locations, and 10 were shared among different geographical locations (Additional file 1: Table S5). The global haplotype diversity (Hd) was 0.73, with a range from 0 for MD and TS to 0.83 for EH (Additional file 1: Table S6).
Polymorphism and demographic statistics inferred from D-loop data for haplogroups in Brandt’s vole
Genetic structure of Brandt’s vole populations
The partitioning of the entire data set using the Bayesian method implemented in the STRUCTURE revealed the presence of two or three genetic clusters, with a maximum log likelihood of posterior probability [lnP(X/K) = −33,067] and maximum (ΔK = 226.2) at K = 2 (Additional file 2: Fig. S1). The sample sites at AL, BE, HH, BG, SH, TS and DR belong to a western cluster, in western Mongolia. The other locations grouped as eastern cluster were predominantly found in Inner Mongolia, China (Fig. 1b & c). With the exception that DR was partitioned into the eastern, the phylogenetic relationships were broadly consistent with NJ tree constructed from the microsatellite pairwise genetic distance (Da) matrix (Fig. 1a). At K = 3 (Additional file 2: Fig. S1), BT, HL, TM, DR and BO emerged as a distinct cluster in the northeastern region, and the other locations were divided into two clusters that geographically corresponded to the western (including MM, CG, HW, AL, BE, HH, BG, SH and TS) and southeastern distributions (including BL, AQ, QL, BX, MD, BY, EH, WL and TG) (Fig. 1b & d).
AMOVA analysis of microsatellite and D-loop data in Brandt’s vole populations
Percentage of variation
Percentage of variation
Among populations within clusters
Among populations within clusters
2 clusters (DR in the western)
Among populations within clusters /groups
The Bayes tree (Fig. 2) and MJ network (Fig. 3a) constructed from the D-loop data revealed a phylogeographic pattern where related haplogroups generally were distributed following our two-cluster model for D-loop sequences (Fig. 3b). Haplogroup 1(Hg1) consisted of 18 haplotypes associated primarily with the eastern region, whereas Hg2 associated mainly with the western region (Fig. 3b; Additional file 1: Table S5). Five haplotypes (H1, H20, H22, H23 and H30) in Hg2 were deviated from the pattern and mainly occurred in the eastern. However, except for H1 that was frequent (found in 27 samples) all others are singletons. All clades of L. brandtii in the Bayes tree give weight to a monophyletic species. Estimates of posterior probability in the Bayes tree (Fig. 2) support the monophyly of L. brandtii, consisting of haplogroups 1 and 2.
AMOVAs using D-loop data revealed that 52.46% (P < 0.0001) of the total variation was distributed among populations for all samples; 75.99% (P < 0.0001) of the total variation was distributed among the two haplogroups; and 59.7% (P < 0.0001) of the total variation was distributed between the eastern and western clusters, resulting in maximization of genetic differences (Fig. 3c; Table 2). Under the three clusters scenario a 50.86% (P < 0.001) of the total variation was distributed among clusters (Table 2). Divergence time estimation showed that Hg1 and Hg2 split occurred at about 0.053Mya (Million years age) with a 95% highest posterior density (HPD) of 0.017–0.088. The basal differentiation of Hg1 clade occurred at 0.041 Mya and 0.035 Mya for Hg2 (Fig. 2).
Reconstruction of ancestral area and migration route
Based on the estimates of divergence time, Hg1(eastern) appears to have diverged earlier than Hg2 (Fig. 2). The ancestral areas inferred using S-DIVA and DEC were in broad agreement. The basal node 65 from S-DIVA displayed two possible ancestral ranges, SE (southeastern) and NE+ SE (northeastern + southeastern), and the probability of the two ancestral ranges were 63.5% and 20.9%, respectively (Fig. 3d; Additional file 1: Table S8). DEC yielded the same two ancestral ranges with S-DIVA, and the probability of these ranges at node 65 were 0.76 and 0.12 for SE and NE + SE, respectively (Additional file 1: Table S8). These results provided support for the ancestors of Brandt’s vole occurring in the southeastern portion of its distribution (Fig. 3c).
Marginal log-likelihoods and model probabilities for 15 migration models (Fig. 4) among two and three clusters using microsatellite data in Brandt’s vole
Number of samplings
5 × 107
5 × 107
Range expansion analysis of Brandt’s vole
Analysis of microsatellite data divided Brandt’s voles population into three clusters from southeastern, northeastern and western of its current distribution (Fig. 1b). Using this population structure and patterns of migration and population expansion we posited the hypothesis that the conditions were in place for gene surfing to have occurred in this species.
Allele frequencies showing spatial clines in southeast, northeast and west clusters
spatial correlation (R2)
spatial correlation (R2)
spatial correlation (R2)
spatial correlation (R2)
Brandt’s vole is a species of Lasiopodomys that was recently separated from the genus Microtus . The species has some peculiar biological and ecological characteristics, i.e., seasonal and multi-annual population fluctuations [17, 18] and the discrete and homogeneous habitat environment. Our work described here is an attempt to re-construct the geographic origin of the species in current distribution and subsequent events that led to the present-day distribution of the species. In addition, during this study the phenomenon of genetic surfing emerged as a possible explanation for the peculiar patterns of genetic differentiation, and the success of genetic lineages.
Ancestral area of Brandt’s vole
Phylogeographic and STRUCTURE analysis of microsatellite data and D-loop sequences revealed that Brandt’s vole currently is comprised out of two to three genetically differentiated populations (Figs. 2 and 3b). To refer to the spatial structuring of the species we refer to these as eastern and western clusters (two symbols in Figs. 1b and 3c), or if we consider are more refined structuring we refer to these as southeastern, northeastern and western clusters (two symbols in Figs. 1b and 3c). We used the three clusters in the analyses of origination and colonization route, the southeastern cluster presented the maximum probability ancestry (P = 63.5% for S-DIVA and P = 0.76 for DEC) of the most primitive ancestral haplotype inferred by RASP and Lagrange (Fig. 3d; Additional file 3: Table S9). These results strongly supported the southeastern region (specifically Xilinhaote District, Inner Mongolia of China) as being the ancestral location of Brandt’s vole prior to its expansion (Fig. 3c). This coincides with the pleistocene fossil evidence observed in northern China .
Range expansion of Brandt’s vole
Neutral tests (Table 1), mismatch analysis (Fig. 5) and MJ network (Fig. 3b; Additional file 1: Table S5) of D-loop sequences indicating a past rapid expansion of the population of Brandt’s vole from a few founders [83, 89, 94, 95]. The diversity values (Table 1), mismatch distributions (Fig. 5b & c) and BSP (Fig. 5e & f) revealed that Hg2 expanded later than Hg1, with the temporal scale and spatial scope essentially agreeing with the optimal colonization route, which was from the ancestral area (the southeastern) to the western through the northeastern (Fig. 4). These findings supported a scenario where a small number of founders disperse and initiate colonization followed by a rapid population expansion . Pleistocene fossils of the species further indicate that its previous distribution was wider than the current one . Thus, the molecular data and the fossil evidence remain to be reconciled, the fact that the past distribution was wider than the present distribution implies that the populations spread from their initial location.
Range expansion can be viewed as a series of successive founder events  that result in a decrease in intrademe heterozygosity with increasing distance from the ancestral population to other demes within a region [88–91]. The IBD for all Brandt’s vole populations showed that the population expansion generally satisfied the one-dimensional stepping-stone model (R2 = 0.149, P < 0.0001) . However, none of the regressions for the southeastern, northeastern and western were significant, which suggests the possible existence of geographic barriers or a greatly heterogeneous environment. Therefore, we used pairwise F ST values instead of the geographical distance to test the range expansion of Brandt’s vole. The results that the H e values significantly decrease with increasing F ST values to the probable arrival (BT and AL) or starting (BX and BO) sites for the species (Fig. 8) supported the hypothesis of range expansion.
The fossil evidence indicates that Lasiopodomys was present in the late Early Pleistocene . We inferred the expansion time of L.brandtii happed earlier than 0.041 Mya (about middle-late Pleistocene), which is for the species’ diversion into Hg1 clade (Fig. 2). In this period, the climate changed towards cooler and more arid, causing land degeneration and desertification [28, 29], big area of meadow-steppes degenerated gradually into the dry steppes and semi-deserts in Mongolia [3, 10, 96], where were favourable habitats for Brandt’s voles. That might have promoted the dispersal and increase in the species populations over time [17, 24]. Towards the Holocene, the climate gradually became mild, resulting expanding forest in the west and central Mongolia  and the Transbaikal plains, Russia . As a result, the suitable habitats of Brandt’s vole significantly decreased. The changes would destroy some preferred habitats of Brandt’s vole and cause fragmentation and reduction in their distribution , which possibly contributed to the differences between the populations from the original and colonized regions.
Genetic surfing in Brandt’s vole population
A steady reduction in genetic diversity with increasing F ST to the starting site centroid (BX and BO) or the probable arrival site centroid (BT and AL) within each of the three regions (Fig. 8) .
The alleles present on the wave front of an expansion could increase in frequency and reach very high proportions and even fixations far away from their original areas [39, 98, 99], and this change is particularly observed in the frequencies of rare alleles . In the surveyed Brandt’s vole populations, some rare alleles in both the southeastern and the northeastern have become common in the corresponding expansion ranges, namely the northeastern and western (Additional file 3: Table S9). Additionally, more clinal patterns of allele frequencies involving different loci are found in the colonized areas compared with the original areas (Table 4), although this difference is not significant.
The spatial differentiation is stronger in the more recently established range than in the original one , and this finding is also supported by the mean pairwise F ST values within each region calculated from microsatellite data, which were found to equal 0.042, 0.044 and 0.059 for the western, southeastern and northeastern, respectively (Additional file 1: Table S7).
Genetic drift acting on the advancing front could cause the extant spatial genetic patterns of Brandt’s vole during its expansion within the distribution. Some individuals, particularly those in an outbreak population, possibly disperse outward from the original habitat due to intraspecific competition. These voles were likely the pioneers in the new habitats. Such dispersal can lead to large gene frequency changes and can determine the genetic diversity of colonies propagated by aggressors , as was typically observed in the populations from BX in the southeastern to BT in the northeastern, and from BO in the northeastern to AL in the western (Fig. 1b). Consequently, this distribution would lead to increased global genetic differentiation and strengthen the presence of migration among populations [91, 99, 101, 102]. These patterns are also supported by the results of genetic structure analysis (Fig. 1; Table 2) and the most likely migration models (Fig. 5; Table 3) inferred from microsatellite data.
The spatial genetic patterns of Brandt’s vole populations are unlikely to have been produced by two important alternative scenarios: adaptive events generated by selective sweeps [27, 103] and introgression [104, 105]. For example, clines in allele frequency are often attributed to geographic variations in selection intensity . Although Brandt’s vole distribution was composed of many patchy homogeneous habitats , the geographic variation in selection intensity was low. Genetic patterns such as gradients of introgression [104, 105] or bidirectional introgression close to the introduction area of an invasive species are often interpreted as footprints of selection . Introgression generally occurs in inter-divergent geographical populations, and the divergence level of the two/three clusters in Brandt’s vole population did not support the taxonomy of subspecies (Fig. 1a).
A combination of demographic properties including population size, growth rate and migration rate, determines the success of mutations in populations under range expansion [39, 106, 107]. Our results illustrate the existence of the surfing phenomenon during Brandt’s vole range expansion, but the demographic properties of the successful colonized populations, such as the populations around BX and BO sites, require further investigation. Otherwise, the heterogeneous landscape might play an important role in determining the fate of mutations, such as their locations [108, 109]. It is undoubtedly more complex than corroborating the effects of these processes in natural populations.
Hypothesis of colonization pattern in Brandt’s vole
In a patch structure of an environment, the habitat types selected by individuals should be a combined result of the intrinsic “quality” of the habitat type and the population density (intensity of competition) in the habitat [110, 111]. Population outbreaks of Brandt’s voles occur irregularly, with an interval of five to seven years [17, 18]. Based on repetitions of this pattern of spatial change throughout history, an outbreak population with opportunistic demographic properties, i.e., suitable population size, growth rate and migration rate, inhabiting the surroundings of the BX site located in the current southeastern distribution was expanding spatially. Some “excess” individuals arrived and successfully settled in a new habitat in the area surrounding the BT site, which was never previously occupied by the species and that is equal in quality to the original habitat. These settlers acted as a seed for the population colonizing the new territories located in the current northeastern distribution. Similar processes occurred during the colonization from BO site in the northeastern to AL site in the western (Fig. 1c) and perhaps further westward until Brandt’s vole population covered the whole distribution area. Such expansion processes can lead to the spatial patterns of genetic structure coinciding with the traits of the genetic surfing phenomenon [39, 101, 112]. Our results from genetic data match the predictions made based on the genetic surfing theory.
We hypothesized that Brandt’s voles originated from the extant southeastern distribution. Under the integrative effect of various factors, such as temperature, rainfall, vegetation and landscape traits, this species colonized from the original area surrounding BX in Xilinguole region, propagating by wave within the region to reach the probable arrival site BT in the northeastern and expanding thereafter within the area; the subsequent expansion of Brandt’s voles population to the western at the arrival site AL from the starting site BO in the northeastern followed the same surfing pattern. Landscape heterogeneity, a low density of individuals in the front of the wave and a low dispersal capacity promote the strong genetic structuration of the three clusters. Brandt’s vole thus formed the current patterns, which are characterized by isolated, patchy, unstable habitats with some genetic characteristics.
We thank Yi Jin, Wurenqimuge, Subuda (Inner Mongolia Agricultural University, China), Enkhbold Nanj and Nyamdorj Batsaikhan (Mongolian State University of Agriculture, Mongolia) for providing assistance with the field trap investigations; Dr. Marcel Holyoak (University of California, Davis, U.S.A.) and Dr. Ying Song (Institute of Plant Protection, CAAC) for providing comments on an earlier version of this manuscript.
Avaliability of data and materials
D-loop haplotype sequences from Brandt’s voles were submitted to GenBank®: Accession nos. KY354521-KY354550. Additional analyses are also provided as supporting information were deposited in GenBank.
This work was supported by International Science & Technology Cooperation Program of China (2014DFG31760) and the National Natural Science Foundation of China (31,371,955; 30,800,727).
DW and KL conceived the experiments. DW, DZS and XRW contributed samples. KL, SMZ and DW conducted the experiments. KL and DW analysed the data. KL, DW and M.H.K wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Permission to trap animals was approved by the Animal Care & Welfare Committee of China Agricultural University (CAU20080823–1).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hewitt GM. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc. 1996;58(3):247–76.View ArticleGoogle Scholar
- Hewitt GM. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13.View ArticlePubMedGoogle Scholar
- Petrova TV, Zakharov ES, Samiya R, Abramson NI. Phylogeography of the narrow-headed vole Lasiopodomys (Stenocranius) gregalis (Cricetidae, Rodentia) inferred from mitochondrial cytochrome b sequences: an echo of Pleistocene prosperity. J Zool Syst Evol Res. 2015;53(2):97–108.View ArticleGoogle Scholar
- Brunhoff C, Galbreath KE, Fedorov VB, Cook JA, Jaarola M. Holarctic phylogeography of the root vole (Microtus oeconomus): implications for late Quaternary biogeography of high latitudes. Mol Ecol. 2003;12(4):957–68.View ArticlePubMedGoogle Scholar
- Runck AM, Cook JA. Postglacial expansion of the southern red-backed vole (Clethrionomys gapperi) in North America. Mol Ecol. 2005;14(5):1445–56.View ArticlePubMedGoogle Scholar
- Heckel G, Burri R, Fink S, Desmet JF, Excoffier L. Genetic structure and colonization processes in European populations of the common vole. Microtus arvalis Evolution. 2005;59(10):2231–42.View ArticlePubMedGoogle Scholar
- Abramson NI, Petrova TV, Dokuchaev NE, Obolenskaya EV, Lissovsky AA. Phylogeography of the grey red-backed vole Craseomys rufocanus (Rodentia: Cricetidae) across the distribution range inferred from nonrecombining molecular markers. Russian Journal of Theriology. 2012;11(2):137–56.View ArticleGoogle Scholar
- Prost S, Guralnick RP, Waltari E, Fedorov VB, Kuzmina E, Smirnov N, et al. Losing ground: past history and future fate of Arctic small mammals in a changing climate. Glob Chang Biol. 2013;19(6):1854–64.View ArticlePubMedGoogle Scholar
- Shenbrot GI, Krasnov BR. Atlas of the geographic distribution of the arvicoline rodents of the world (Rodentia, Muridae: Arvicolinae). Moscow: Pensoft Press; 2005.Google Scholar
- Alexeeva N, Erbajeva M, Khenzykhenova F. Lasiopodomys brandtii In Pleistocene of Transbaiklia and adjacent territories: distribution area, evolutionary development in context of global and regional events. Quat Int. 2015;355:11–7.View ArticleGoogle Scholar
- Lyapunova EA, Mirokhanov YM. Description of the Field Vole Chromosome Sets (Stenocranius, Lasiopodomus, Blanfordimus, Microtus). Materialy II Vses. sov. po mlekopitayushchim,1969; 134–138.Google Scholar
- Chaline J. Les rongeurs du pléistocène moyen et supérieur de France (systématique- biostratigraphie- paléoclimatologie). Cahiers de paléontologie. 1972;Google Scholar
- Erbajeva MA, Sotnikova MV, Shevtchenko VK. New Eopleistocene locality of mammal fauna in Transbaikalia. Paleontological Foundation of the Anthropogene Stratigraphy, Moscow. 1977;103Google Scholar
- Agadzhanyan AK, Yatsenko VN. Phylogenetic interrelationships in voles of northern Eurasia. Archives of Zoological Museum of Moscow State University. 1984;22:135–90. (in Russian)Google Scholar
- Abramson NI, Golenishchev FN, Kostygov AY, Tesakov AS. Taxonomic interpretation of the molecular-genetic cladogram of tribe Microtini (Arvicolinae, Rodentia) based on nuclear genes. In: Materials of the International Meeting “Theriofauna of Russia and Adjacent Areas”. 2011; 1e4, Moscow, 7 (in Russian).Google Scholar
- Zhong WQ, Zhou QQ, Sun CL. The basic characteristics of the rodent pests on the pasture in Inner Mongolia and the ecological strategies of controlling. Acta Theriologica Sinica. 1985;5:241–9. (in Chinese)Google Scholar
- Zhang ZB, Pech R, Davis S, Shi DZ, Wan XR, Zhong WQ. Extrinsic and intrinsic factors determine the eruptive dynamics of Brandt’s voles. Microtus brandti in Inner Mongolia, China. Oikos. 2003a;100:299–310.View ArticleGoogle Scholar
- Li ZL, Liu TC. Population cycle of Brandt’s vole (Microtus brandti) and rodent community succession in Xilingele-Meng. Zoology Research. 1999;20:284–7. (in Chinese)Google Scholar
- Wang WG, Wu HS, Chen SK, Gao HB. Rodent and damages. Control of grassland pests in Hulunbeir, China, China Agricultural Press. 2005:98–101. (in Chinese)Google Scholar
- Jiang YE, Enkhbold N, Batsaikhan N, Wang D, Yi J. Wurenqimuge, Shan LY, du GL, Shi D. Study on the habitat characteristics of Brandt’s vole. Acta Agrestia Sinica. 2012;20(2):179–82. (in Chinese)Google Scholar
- Zhong WQ, Wang MJ, Wan XR. Ecological management of Brandt’s vole (Microtus brandti) in Inner Mongolia, China. In: Singleton GR, Hinds LA, Leirs H, Zhang Z, editors. Ecologically-based management of rodent pests. Canberra: ACIAR Monograph. Australian Centre for International Agricultural Research; 1999. p. 199–214.Google Scholar
- Zhong WQ, Zhou QQ, Wang GH, Sun CL, Zhou PL, Liu WZ, et al. The design for the ecological management of Brandt’s vole pest and its application. Acta Theriologica Sinica. 1991;11:204–12. (in Chinese)Google Scholar
- Luo TS, Hao SS, Liang ZA, Niu DF, Zao HC. Some biological observations of the Brandt’s vole (Microtus brandti Radde) control in Hulun-Beier grassland. Acta Zool Sin. 1975;21:5–61. (in Chinese)Google Scholar
- Zhang ZB, Zhong WQ, Fan NC. Rodent problems and management in the grasslands of China. ACIAR Monograph series. 2003b;96:316–9.Google Scholar
- An CB, Chen FH, Barton L. Holocene environmental changes in Mongolia: a review. Glob Planet Chang. 2008;63(4):283–9.View ArticleGoogle Scholar
- Lehmkuhl F, Hilgers A, Fries S, Hülle D, Schlütz F, Shumilovskikh L, et al. Holocene geomorphological processes and soil development as indicator for environmental change around Karakorum, upper Orkhon Valley (Central Mongolia). Catena. 2011;87(1):31–44.View ArticleGoogle Scholar
- Excoffier L, Foll M, Petit RJ. Genetic consequences of range expansions. Annu Rev Ecol Evol Syst. 2009;40:481–501.View ArticleGoogle Scholar
- Li JH, Pan HW, Wang G. Degradation causes of typical steppe in Inner Mongolia. Pratacultural Scientia. 2004;21:49–51. (In Chinese)Google Scholar
- Li SY, Li XB, Wang DD. Prediction of grassland degradation in Xilinhaote of Inner Mongolia based on Markov process model. Chinese Journal of Ecology. 2007;26:78–82. (In Chinese)Google Scholar
- Natsagdorj L, Dulamsuren J. Some aspect of assessment of the dzud phenomena. Papers in Institute of Meteorology and Hydrology. 2001;23(3):3–18.Google Scholar
- Hoshino B, Ganzorig S, Sawamukai M, Kawashima K, Baba K, Kai K, Nurtazin S. The impact of land cover change on patterns of zoogeomorphological influence: Case study of zoogeomorphic activity of Microtus brandti and its role in degradation of Mongolian steppe. In 2014 IEEE Geoscience and Remote Sensing Symposium. 2014; 3518–3521.Google Scholar
- Zhang C, Zhang W, Feng Z, Mischke S, Gao X, Gao D, et al. Holocene hydrological and climatic change on the northern Mongolian plateau based on multi-proxy records from Lake gun Nuur. Palaeogeogr Palaeoclimatol Palaeoecol. 2012;323:75–86.View ArticleGoogle Scholar
- Liu H, Xu L, Cui H. Holocene history of desertification along the woodland-steppe border in northern China. Quat Res. 2002;57(2):259–70.View ArticleGoogle Scholar
- Feng ZD. Gobi dynamics in the northern Mongolian plateau during the past 20,000+ yr: preliminary results. Quat Int. 2001;76:77–83.View ArticleGoogle Scholar
- Hancock AM, Witonsky DB, Gordon AS, Eshel G, Pritchard JK, Coop G, et al. Adaptations to climate in candidate genes for common metabolic disorders. PLoS Genetic. 2008;4(2):32.View ArticleGoogle Scholar
- Coop G, Pickrell JK, Novembre J, Kudaravalli S, Li J, Absher D, et al. The role of geography in human adaptation. PLoS Genetic. 2009;5(6):e1000500.View ArticleGoogle Scholar
- Novembre J, Di Rienzo A. Spatial patterns of variation due to natural selection in humans. Nat Rev Genet. 2009;10(11):745–55.View ArticlePubMed CentralPubMedGoogle Scholar
- Storz JF, Wheat CW. Integrating evolutionary and functional approaches to infer adaptation at specific loci. Evolution. 2010;64(9):2489–509.View ArticlePubMed CentralPubMedGoogle Scholar
- Klopfstein S, Currat M, Excoffier L. The fate of mutations surfing on the wave of a range expansion. Mol Biol Evol. 2006;23(3):482–90.View ArticlePubMedGoogle Scholar
- Wood MD, Slade NA. Comparison of ear-tagging and toe-clipping in prairie voles, Microtus Ochrogaster[J]. J Mammal. 1990;71(2):252–5.View ArticleGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual [M]. Cold spring harbor laboratory press. 1989;Google Scholar
- Wang D, Shi DZ. Isolation and characterization of polymorphic microsatellite loci from Brandt’s voles (Lasiopodomys brandtii). Mol Ecol Notes. 2007;7:671–3.View ArticleGoogle Scholar
- Yue L, Wang D, Huang B, Liu X. Characterization of nine novel microsatellite markers from Brandt’s vole (Lasiopodomys brandtii). Mol Ecol Resour. 2009;9:1194–6.View ArticlePubMedGoogle Scholar
- Van OC, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.View ArticleGoogle Scholar
- Valière N. GIMLET: a computer program for analysing genetic individual identification data. Mol Ecol Notes. 2002;2(3):377–9.Google Scholar
- Lalitha S. Primer premier 5. Biotech Software & Internet Report: The Computer Software Journal for Scient. 2000;1(6):270–2.View ArticleGoogle Scholar
- Burland TG. DNASTAR’s Lasergene sequence analysis software. Bioinformatics Methods and Protocols. 1999;132:71–91.View ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.View ArticlePubMed CentralPubMedGoogle Scholar
- Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.View ArticleGoogle Scholar
- Rice WR. Analyzing tables of statistical tests. Evolution. 1989;43:223–5.View ArticlePubMedGoogle Scholar
- Peakall R, Smouse PE. GenAlEx 6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.View ArticleGoogle Scholar
- Kalinowski ST. The computer program STRUCTURE does not reliably identify the main genetic clusters within species: simulations and implications for human population structure. Heredity. 2011;106:625–32.View ArticlePubMedGoogle Scholar
- Takezaki N, Nei M, Tamura K. POPTREE2: software for constructing population trees from allele frequency data and computing other population statistics with windows interface. Mol Biol Evol. 2010;27(4):747–52.View ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMed CentralPubMedGoogle Scholar
- Earl D.A., vonHoldt B.M. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. 2012; 4: 359–361.Google Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.View ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S. Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005;1:47–50.Google Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.View ArticlePubMedGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.View ArticlePubMed CentralPubMedGoogle Scholar
- MartíNková N, Moravec J. Multilocus phylogeny of arvicoline voles (Arvicolini, Rodentia) shows small tree terrace size[J]. Folia Zool. 2012;61(3/4):254.Google Scholar
- Lemskaya NA, Romanenko SA, Golenishchev FN, Rubtsova NV, Sablina OV, Serdukova NA, et al. Chromosomal evolution of Arvicolinae (Cricetidae, Rodentia). III. Karyotype relationships of ten Microtus species[J]. Chromosom Res. 2010;18(4):459–71.View ArticleGoogle Scholar
- Bužan EV, Kryštufek B, HÄNFLING B, Hutchinson WF. Mitochondrial phylogeny of Arvicolinae using comprehensive taxonomic sampling yields new insights[J]. Biol J Linn Soc. 2008;94(4):825–35.View ArticleGoogle Scholar
- Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19:716–23.View ArticleGoogle Scholar
- Posada D, Crandall KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14:817–8.View ArticlePubMedGoogle Scholar
- Kass RE, Raftery AE. Bayes factors[J]. J Am Stat Assoc. 1995;90(430):773–95.View ArticleGoogle Scholar
- Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolution Biology. 2007;7:214.View ArticleGoogle Scholar
- Rambaut A, Drummond AJ. Tracer v1.5, Available: http://beast.bio.ed.ac.uk/Tracer. Accessed 2014 Jan 30. 2007.Google Scholar
- Rambaut A. FigTree v1.3.1, Available from: http://tree.bio.ed.ac.uk/software/figtree. 2009.
- Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.View ArticlePubMedGoogle Scholar
- ESRI (Environmental Systems Resource Institute). ArcMap 10.1. ESRI, Redlands, California. 2012.Google Scholar
- Yu Y, Harris AJ, He X. S-DIVA (statistical dispersal-Vicariance analysis): a tool for inferring biogeographic histories. Mol Phylogenet Evol. 2010;56:848–50.View ArticlePubMedGoogle Scholar
- Yu Y, Harris AJ, He XJ. RASP (Reconstruct Ancestral State in Phylogenies) 2.1 beta. Available: http://mnh.scu.edu.cn/soft/blog/RASP. Access: 2014 Jan 30. 2013.
- Ree R, Moore BR, Webb CO, Donoghue MT. A likelihood framework for inferring the evolution of geographic range on phylogenetic trees. Evolution. 2005;59:2299–311.View ArticlePubMedGoogle Scholar
- Ree RH, Smith SA. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol. 2008;57:4–14.View ArticlePubMedGoogle Scholar
- Beerli P. Comparison of Bayesian and maximum likelihood inference of population genetic parameters. Bioinformatics. 2006;22:341–5.View ArticlePubMedGoogle Scholar
- Beerli P, Palczewski M. Unified framework to evaluate panmixia and migration direction among multiple sampling locations. Genetics. 2010;185:313–26.View ArticlePubMed CentralPubMedGoogle Scholar
- Kuhner MK. LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006;22:768–70.View ArticlePubMedGoogle Scholar
- Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. Springer Science & Business Media. 2003;Google Scholar
- Fu Y, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709.PubMed CentralPubMedGoogle Scholar
- Fu Y. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.PubMed CentralPubMedGoogle Scholar
- Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMed CentralPubMedGoogle Scholar
- Slatkin M, Hudson RR. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991;129:555–62.PubMed CentralPubMedGoogle Scholar
- Rogers AR, Harpending H. Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992;9:552–69.PubMedGoogle Scholar
- Wang Y, Zhao LM, Fang FJ, Liao JC, Liu NF. Intraspecific molecular phylogeny and phylogeography of the Meriones meridianus (Rodentia: Cricetidae) complex in northern China reflect the processes of desertification and the Tianshan Mountains uplift [J]. Biol J Linn Soc. 2013;110(2):362–83.View ArticleGoogle Scholar
- Rousset F. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997;145:1219–28.PubMed CentralPubMedGoogle Scholar
- SPSS. SPSS advanced statistics 20.0. SPSS: Chicago, IL; 2011.Google Scholar
- Graciá E, Botella F, Anadón JD, Edelaar P, Harris DJ, Giménez A. Surfing in tortoises? Empirical signs of genetic structuring owing to range expansion. Biol Lett. 2013;9(3):20121091.View ArticlePubMed CentralPubMedGoogle Scholar
- Austerlitz F, Jung-Muller B, Godelle B, Gouyon PH. Evolution of coalescence times, genetic diversity and structure during colonization. Theor Popul Biol. 1997;51(2):148–64.View ArticleGoogle Scholar
- Wegmann D, Currat M, Excoffier L. Molecular diversity after a range expansion in heterogeneous environments. Genetics. 2006;174(4):2009–20.View ArticlePubMed CentralPubMedGoogle Scholar
- DeGiorgio M, Degnan JH, Rosenberg NA. Coalescencetime distributions in a serial founder model of human evolutionary history. Genetics. 2011;189:579–93.View ArticlePubMed CentralPubMedGoogle Scholar
- Slatkin M, Excoffier L. Serial founder effects during range expansion: a spatial analog of genetic drift. Genetics. 2012;191:171–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Felauer T, Schlütz F, Murad W, Mischke S, Lehmkuhl F. Late Quaternary climate and landscape evolution in arid Central Asia: a multiproxy study of lake archive Bayan Tohomin Nuur¢, Gobi desert, southern Mongolia. J Asian Earth Sci. 2012;48:125–35.View ArticleGoogle Scholar
- Luo ZX, Chen W, Gao W. Fauna Sinica Mammalia (Vol 6, Rodentia). Part III: Cricetidae. Science Press, Beijing; 2000. (in Chinese)Google Scholar
- Li WH. Distribution of nucleotide differences between two randomly chosen cistrons in a finite population. Genetics. 1977;85:331–7.PubMed CentralPubMedGoogle Scholar
- Ray N, Currat M, Excoffier L. Intra-deme molecular diversity in spatially expanding populations. Mol Biol Evol. 2003;20(1):76–86.View ArticlePubMedGoogle Scholar
- Peck JA, Khosbayar P, Fowell SJ, Pearce RB, Ariunbileg S, Hansen BC, et al. Mid to late Holocene climate change in north central Mongolia as recorded in the sediments of Lake Telmen[J]. Palaeogeogr Palaeoclimatol Palaeoecol. 2002;183(1):135–53.View ArticleGoogle Scholar
- Wu ZJ, Li YM. Effects of habitat fragmentation on survival of animal populations. Acta Ecol Sin. 2003;23:2423–35. (In Chinese)Google Scholar
- Travis JM, Münkemüller T, Burton OJ, Best A, Dytham C, Johst K. Deleterious mutations can surf to high densities on the wave front of an expanding population. Mol Biol Evol. 2007;24(10):2334–43.View ArticlePubMedGoogle Scholar
- Hallatschek O, Nelson DR. Life at the front of an expanding population. Evolution. 2009;64:193–206.View ArticlePubMedGoogle Scholar
- Hofer T, Ray N, Wegmann D, Excoffier L. Large allele frequency differences between human continental groups are more likely to have occurred by drift during range expansions than by selection. Ann Hum Genet. 2009;73(1):95–108.View ArticlePubMedGoogle Scholar
- Excoffier L, Ray N. Surfing during population expansions promotes genetic revolutions and structuration. Trends Ecol Evol. 2008;23(7):347–51.View ArticlePubMedGoogle Scholar
- Hallatschek O, Hersen P, Ramanathan S, Nelson DR. Genetic drift at expanding frontiers promotes gene segregation. Proc Natl Acad Sci. 2007;104:19926–30.View ArticlePubMed CentralPubMedGoogle Scholar
- Currat M, Excoffier L, Maddison W, Otto SP, Ray N. Comment on “Ongoing adaptive evolution of ASPM, a brain size determinant in Homo sapiens” and “Microcephalin, a gene regulating brain size, continues to evolve adaptively in humans.” Science 2006; 313: 172; author reply 72.Google Scholar
- Rendine S, Piazza A, Cavalli-Sforza LL. Simulation and separation by principal components of multiple demic expansions in Europe. Am Nat. 1986;128:681–706.View ArticleGoogle Scholar
- Currat M, Ruedi M, Petit RJ, Excoffier L. The hidden side of invasions: massive introgression by local genes. Evolution. 2008;62:1908–20.PubMedGoogle Scholar
- Kimura M. On the probability of fixation of mutant genes in a population. Genetics. 1962;47:713–9.PubMed CentralPubMedGoogle Scholar
- Maruyama T. On the fixation probability of mutant genes in a subdivided population. Genet Res. 1970;15(2):221–5.View ArticlePubMedGoogle Scholar
- Burton OJ, Travis JM. The frequency of fitness peak shifts is increased at expanding range margins due to mutation surfing. Genetics. 2008a;179:941–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Burton OJ, Travis JM. Landscape structure and boundary effects determine the fate of mutations occurring during range expansions. Heredity. 2008b;101:329–40.View ArticlePubMedGoogle Scholar
- Fretwell SD, Lucas HL. On territorial behavior and other factors influencing habitat distribution in birds. Acta Biotheor. 1969;19:16–36.View ArticleGoogle Scholar
- Fretwell SD. Populations in a seasonal environment (No. 5). Princeton University Press. 1972.Google Scholar
- Edmonds CA, Lillie AS, Cavalli-Sforza LL. Mutations arising in the wave front of an expanding population. Proc Natl Acad Sci U S A. 2004;101(4):975–9.View ArticlePubMed CentralPubMedGoogle Scholar