Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(10): e72855.
Published online 2013 October 9. doi:  10.1371/journal.pone.0072855
PMCID: PMC3794015

Evaluating the Significance of Paleophylogeographic Species Distribution Models in Reconstructing Quaternary Range-Shifts of Nearctic Chelonians

Michael Hofreiter, Editor


The climatic cycles of the Quaternary, during which global mean annual temperatures have regularly changed by 5–10°C, provide a special opportunity for studying the rate, magnitude, and effects of geographic responses to changing climates. During the Quaternary, high- and mid-latitude species were extirpated from regions that were covered by ice or otherwise became unsuitable, persisting in refugial retreats where the environment was compatible with their tolerances. In this study we combine modern geographic range data, phylogeny, Pleistocene paleoclimatic models, and isotopic records of changes in global mean annual temperature, to produce a temporally continuous model of geographic changes in potential habitat for 59 species of North American turtles over the past 320 Ka (three full glacial-interglacial cycles). These paleophylogeographic models indicate the areas where past climates were compatible with the modern ranges of the species and serve as hypotheses for how their geographic ranges would have changed in response to Quaternary climate cycles. We test these hypotheses against physiological, genetic, taxonomic and fossil evidence, and we then use them to measure the effects of Quaternary climate cycles on species distributions. Patterns of range expansion, contraction, and fragmentation in the models are strongly congruent with (i) phylogeographic differentiation; (ii) morphological variation; (iii) physiological tolerances; and (iv) intraspecific genetic variability. Modern species with significant interspecific differentiation have geographic ranges that strongly fluctuated and repeatedly fragmented throughout the Quaternary. Modern species with low genetic diversity have geographic distributions that were highly variable and at times exceedingly small in the past. Our results reveal the potential for paleophylogeographic models to (i) reconstruct past geographic range modifications, (ii) identify geographic processes that result in genetic bottlenecks; and (iii) predict threats due to anthropogenic climate change in the future.


Quaternary glacial-interglacial cycles have caused dramatic range shifts of species' distributions [1][5]. Many species in high- and mid-latitudes have been extirpated from northern and high altitude areas, often accompanied by range fragmentation and persistence in so-called glacial refugia [6][9]. These geographic processes have led to phylogeographic differentiation at the population, subspecies and species levels in both plants and animals [10][11]. Range fragmentation has led to allopatric genetic differentiation, and range expansions have led to the formation of tension hybrid zones, processes that lead to genetic diversity [8], [12][15]. Large-scale geographic range contractions have also resulted in the loss of genetic diversity through repeated founder effects and population bottlenecks [16]. Species can respond to changing climate in two non-exclusive ways: (i) adapting its climatic niche to fit the new conditions or (ii) changing its geographic range to track the conditions that are compatible with its climatic niche [17][19]. One of the biggest gaps in our understanding of biotic responses to climate change is the relative balance between the two. Many studies have addressed geographic range tracking in relation to past and future climates [20][25], but only a few have attempted to evaluate the rates and magnitudes by which species respond by adapting to climate change (e.g. [26][29]). This area is in urgent need of investigation because it is highly relevant to assessing the impact of rapidly changing climate.

This study integrates the Quaternary fossil record, the oxygen isotopic record of global temperature history, paleoclimatic models of the geographic distribution of past climates, phylogenetic comparative analysis, and species distribution modeling to better understand the tempo and mode of how turtle species have responded to past climate changes. These integrated data are used to produce testable, geographically explicit, temporally continuous models of species ranges through the full range of climate variability of the last three Quaternary climate cycles. We test our models against fossil occurrence data, which are the only positive source of information about where species occurred in the past. The fossil occurrences test whether the past geographic ranges of these turtle species can accurately be reconstructed from paleoclimate, which is one of the fundamental assumptions of our paleophylogeographic modeling approach, indeed to all ecological niche modeling. The models are then used to evaluate the roles that past changes in climate and geography have played in determining physiological, genetic, and phylogeographic patterns. Finally, our results are used to partition the responses of turtles to climate change into phylogenetic (adaptive) and purely geographic components to assess the relative contribution of each mode.

We call our approach paleophylogeographic modeling. It is derived from species distribution modeling (SDM), uses georeferenced occurrences to characterize the range of climate that a species actually inhabits (the realized climate niche) and which has become the method of choice for projecting modern climatic niches onto paleoclimate reconstructions [30][32]. By assuming that the realized niche actually represents its true climatic niche (the fundamental climate niche) [24], the expected geographic range of a species can be determined for any future and past climate model. Our paleophylogeographic approach expands on the SDM approach in several important ways: (i) it uses a temporally continuous model of changing paleoclimate with which to model changes in turtle ranges through the past 350 Ka; (ii) it uses phylogenetic comparative methods to account for evolutionary changes that happened in the climate niche that occurred over the period being modeled; (iii) it uses fossil occurrences to test whether the reconstructed ranges are accurate; (iv) it also uses the fossil occurrences to test the core assumption that the modern realized niche is a good proxy for the fundamental climate niche; and (v) it partitions past responses to climate change into adaptive and purely geographic components so that the rates of climatic adaptation and geographic range tracking can be compared.

The conceptual framework of our paleophylogeographic modeling, like species distribution modeling, is based on the concepts of niches in climate space and their relationship to distributions in geographic space [33][35] (Figure 1). Hutchinson conceptualized niches as volumes within a multidimensional space whose axes are variables important to the life of the species [35]. Grinnell, who first introduced the niche concept, limited the concept to abiotic environmental conditions, especially climate [34]. We adopt Grinnell's strictly climatic view of the niche because we are focusing on the effects of Quaternary climatic changes on species ranges. Thus, in our study, the multidimensional niche space is defined by climatic variables (E-space, where “E” stands for environmental). The distribution of a species in this niche space is closely related to its distribution in geographic space (G-space), and vice versa. The niche concept can be logically subdivided into fundamental and realized niches [36][39]. A fundamental niche is the full range of E-space in which a species can flourish – its range of climatic tolerance. The realized niche is the range of E-space that the species actually occupies, which may be smaller. Similarly, a realized geographic distribution is the range of territory in G-space that it actually occupies. Both biotic and accessibility factors can prevent a species from realizing the full niche space or territory where it could otherwise flourish based on its climate tolerances. In E-space, biotic exclusion occurs when one species excludes another from part of the climate it can tolerate, in G-space when one species excludes another from a geographic area. Accessibility in E-space occurs where climate is actually available for the species, whereas in G-space species accessibility occurs when there are no geographic barriers, such as oceans, mountains, or rivers, to prevent a species from occupying a territory that it finds climatically tolerable.

Figure 1
Relationships between fundamental niches, realized niches, potential niches and geographic distributions as (a) a Venn diagram (after [36], [157]); (b) in climatic niche space (E-space); and (c) in geographic space (G-space).

Accessibility in E-space is an important concept in our study, and it relates to what we call analogue and non-analogue climates [40][41]. Some areas of E-space represent climate that does not actually occur in G-space, which means that there is no “accessible climate” for the species in E-space. For example, consider a mean annual temperature of 50°C, which is far greater than any mean annual temperature found on Earth today and would not have a corresponding point in G-space. Points such as this in E-space are called non-analogue climate points. Some non-analogue climates actually existed in the geological past [40][41]. For example, in the early Eocene mean annual temperature was more than 12°C higher than it is today and there were no areas of the globe where low temperatures dropped below the freezing point. The maximum mean annual temperature of the early Eocene world fell well above the maximum boundary in E-space defined by modern climate and was thus a non-analogue with respect to modern climate. Conversely, today's minimum mean annual temperature falls well below the minimum boundary in E-space defined by early Eocene climate and could thus be described as non-analogue with respect to Eocene climate. Accessibility in E-space is thus something that changes as climate changes. In our study, changes in climate accessibility provide important clues to how closely the realized niche approximates the fundamental niche.

In this paper we make another important distinction about potential niches and geographic range [42]. A potential range is the territory where a species is capable of flourishing and which is accessible to it (Figure 1A,C). In other words, it is the territory where we expect a species to live based on its fundamental niche (climate tolerances) in the absence of biotic barriers but in the presence of geographic barriers. A potential niche is the climatic space that falls within the fundamental niche and which actually exists at the time a species lives (available climate). The potential range of a species changes as climate changes, as can its potential niche change as the range of available climate changes. We make use of the changes in potential niches through the last 320 Ka of Quaternary climate change to assess the extent to which turtle species could live in climates that do not exist in the present day.

We make use of fossils in two ways, first to test whether modern niches are congruent with niches in the past by testing the power of modern realized niches to predict past geographic ranges that are compatible with fossil occurrences once phylogenetic changes and changes in available climate have been taken into account. Previous studies that have used fossils as tests have found that SDMs are usually successful at regional and continental scales in estimating the past ranges of species, the location of glacial refugia, and the extent of postglacial range expansions (e.g. [18], [23], [26][27], [43][45]. Furthermore, SDMs constructed from fossil occurrences and paleoenvironmental data have often successfully predicted modern geographic distributions [46][47]. However, there are notable exceptions in which SDMs grossly failed to agree with fossil data or molecular phylogeographic data [27], [30], [43], [48][49]. Thus, we recommend that paleogeographic range reconstructions always be tested against fossil data if at all possible. Independent lines of evidence such as physiological data, the fossil record and genetic signatures offer important insights into niche dynamics [18], [49][50]. We make full use of these tests in this study.

We also use fossils to study changes in the potential niche by determining whether species occupied non-analogous climate conditions in the past, which in turn helps better understand whether changes in the niche over time are due to changes in biotic factors, like interspecies competition, adaptive changes in the niche, or changes in the climate system itself. Non-analogous climates can easily be identified by plotting modern and paleoclimate data in the same E-space: non-analogue climates are simply those paleoclimate data points that fall outside the boundaries of modern climate. A fossil can fall in one of three places in E-space: (1) inside the species potential niche, (2) outside the available paleoclimate, or (3) inside the available paleoclimate but outside the potential niche (Figure 1B). The first case indicates that the potential niche as estimated from the modern realized niche is compatible with the fossil evidence; the second indicates that either the fossil or the paleoclimate data are erroneous; the third indicates that the species occurred in an available climate outside the range predicted by its modern tolerances, which suggests that its niche has changed over time. These three outcomes help discriminate between changing climate accessibility and evolving fundamental niches.

We used multivariate environmental similarity surfaces (MESS) to measure the positions of fossil occurrences in E-space [51][52]. MESS is the Euclidean distance of any given point in E-space to boundaries defined by a set of reference points. The MESS score is positive when the point lies within the reference boundaries and negative when it lies outside them. Non-analogous climate points are thus those points with negative MESS scores that lie outside the boundaries of modern climate (for more details see below). We also used MESS scores to measure the position of fossil occurrences of turtle species, with respect to the boundaries of their realized and potential niches, to determine whether the niches have changed over the last 320 Ka. The MESS scores thus indicate whether a species has tolerated climatic conditions beyond what it experiences today, they help define the full potential niche of the species, and they help forecast the fate of the species over the next century.

An important feature of our paleophylogeographic approach is that we use phylogeny to estimate the evolutionary changes that have occurred in climate tolerances over a clades' history. Because of evolution, we cannot expect a priori that a species' niche in the past was the same as it is today [23], [30], [49], [53]. The logic of this approach is based on the observation that the climatic niches of modern members of a clade differ from one another, whereas the ancestor of the clade would have had a single climatic niche. Using phylogenetic comparative methods, ancestral niches can be estimated at each node and along each branch of the tree [45], [54][57]. We use the phylogenetic information to adjust each species niche in E-space for where it is along the phylogenetic branch leading back to its common ancestor with other species [23]. The phylogenetic changes can be partitioned from geographic tracking in our paleophylogeographic models [23] and thus help us address which mode of response has been more important in the clade's history.

Turtles and tortoises are a good group for studying the response of species to climate change because: (i) their current distributions are well documented [58][60]; (ii) they have a good Quaternary fossil record (e.g. [61] with references therein); (iii) their phylogenetic relationships are well resolved [62]; and (iv) phylogeographic data based on molecular markers are available for many species (Table 1). Even more importantly, turtles and tortoises are ectotherms that are more directly affected by ambient temperature than mammals and birds [55], [63] and thus are expected to have realized niches that are more directly controlled by climate. Indeed, incubation temperature determines sex in most species ([64] and references therein). Fifty nine of the 61 living North American turtles are included in this study.

Table 1
Comprehensive overview of phylogeographic and morphometric analyses of North American turtle species.

To summarize, in this study we use modern geographic distributions to estimate realized and potential niches for modern turtle species. We use paleophylogeographic modeling to reconstruct the history of their geographic range changes over the last 320 Ka using (i) a temporally continuous set of paleoclimatic models that we constructed from paleoclimate general circulation model data, (ii) the oxygen isotope record of global mean annual temperature over the same period of time, and (iii) the phylogenetic changes in the niche that occurred along each species' branch as estimated from the phylogeny of the entire clade. We use fossil occurrences to test whether the modern realized niches are good predictors of paleogeographic distributions and to measure the frequency with which turtles lived in non-analogue climates in the past. We test whether the paleogeographic ranges are congruent with physiological and phylogeographic data and we determine whether evolutionary adaptation or geographic range tracking has been more important as species responded to the three glacial-interglacial cycles included in our study. We use our results to characterize the full potential niches of species and to identify ones that may be more vulnerable to the 21st Century climate changes that are forecast for the next century.

Materials and Methods

Species occurrences

We analyzed 59 of the 61 chelonian species in North America. Although the phylogenetic relationships of these species are well-resolved and comparatively stable, their taxonomy is controversial, especially at the sub-species level [62], [65][67]. For consistency, we have arbitrarily chosen to follow the most recent taxonomic classification [62], but our modeling depends on the phylogeny rather than on the classification. Occurrences for most of the species were compiled from the World Turtle Database [60]. Adequate sample sizes are crucial for sound models [68][69], so for species with fewer than ten records in the World Turtle Database, we added occurrences from the Global Biodiversity Information Facility (GBIF, All occurrences were checked for accuracy. To avoid the effects of uneven sampling, we resampled the point occurrences with a grid with ten arc-minutes resolution, retaining only one occurrence point per grid cell. In total, 19,240 unique resampled points for the 59 species were retained for further analysis.

Fossil occurrences

We compiled fossil occurrences of extant species for the last 320 Ka from published literature [61], [70][74] and the Paleobiology Database ( We did not include occurrences whose location or taxonomy was uncertain. In total 141 fossil occurrences were found for 18 of the 59 species (Appendix S3 in Material S1). The geological age estimates for these occurrences vary considerably in their precision and uncertainty. Most fossils younger than 50 Ka are radiocarbon dated, which provides precise ages accurate to within a few decades. Older fossils extend beyond the range of radiocarbon and are dated with lithostratigraphy and biostratigraphy. Occurrences between 50 Ka and 125 Ka are assigned to one of three Wisconsinan (last glacial) stages or to the Sangamonian interglacial (last interglacial), which have a precision of about 10 Ka. Older occurrences are assigned to one of two North American Land Mammal Ages (NALMA) and have precisions of tens of thousands of years [75]. The maximum and minimum age of each occurrence is reported in Appendix S3 in Material S1. Our strategy for dealing with the age uncertainty is described below.

Modern climate data and paleoclimate interpolation

Modern climate data with a spatial resolution of 2.5 arc min were obtained from WorldClim [76]. This set of 19 bioclimatic variables characterize climate in terms of means (e.g. mean annual temperature) and extremes (e.g. maximum temperature of the coldest quarter) (see below; modified from [77]). We selected a subset of eight moderately to uncorrelated climate variables whose coefficients of determination were <0.75 to avoid the effects that correlations have on Maxent analyses (see below) [78]: mean diurnal temperature range (BIO2), maximum temperature of warmest month (BIO6), mean temperature of wettest quarter (BIO8), mean temperature of warmest quarter (BIO10), precipitation of wettest quarter (BIO16), precipitation of driest quarter (BIO17), precipitation of warmest quarter (BIO18) and precipitation of coldest quarter (BIO19).

We used a continuous temporal series of interpolated paleoclimate models previously developed for the late Quaternary of North America [23]. These interpolated models were constructed by using oxygen isotope data from deep sea cores [79], which are a proxy for global mean annual temperature, to interpolate between modern climate (an interglacial climate end-member) and a global circulation climate model (GCM) of paleoclimate at the last glacial maximum at 21 Ka (MIROC3.2 GCM, [80]) (a glacial climate end-member). Each end-member climate data set consisted of the same 19 bioclimatic variables described above. All the variables were scaled relative to their mean annual temperature using the continuous record of global mean annual temperature from the oxygen isotope record [23]. The interpolated paleoclimate models thus capture the full range of climate variability over the last 320 Ka, including the warmer climates of the last (Sangamonian) interglacial and the intense variability within glacial cycles. The paleoclimate models were constructed at 4 thousand year intervals, the same temporal resolution as the isotope data. The interpolated models were validated by comparing them with a GCM at 6 kya and a GCM that describes climate between 120 and 140 kya (see Figure 3 and Figure 4 in [23]). The interpolated models agreed well with the direct GCM models, differing no more from the GCMs than two GCMs differ from one another. Paleoclimate interpolation was performed at a resolution of 50 km2 equal-area grid cells to avoid latitudinal bias in sampling density [23]. The sampling grid coordinates are available at

Figure 3
Relationships between pairwise changes in mean annual temperatures (MAT) and pairwise changes in potential distribution range sizes as well as between MAT and pairwise changes in geographic centers of potential distributions in 59 Nearctic chelonians ...
Figure 4
Historic niche dynamics in Nearctic chelonians based on 141 fossils and the five most proximal (SIM5) 50 km2 points in climate space in those cases where species find no analogous climates to their past potential niche.

Species distribution models (SDMs)

Modern SDMs were produced for each species using Maxent 3.3.3e [51], [81][82]. Maxent provides good predictions of modern geographic ranges from point occurrence data [83]. Default settings were used with logistic output ranging from 0 (no probability of occurrence) to 1 (optimal) [82]. For model training, the available climate space was limited to the climate within the Level 2 watersheds where the species is known to occur (see below). 100 SDMs were computed for each species, each trained with a randomly selected 70% of the occurrences and evaluated with the remaining 30% of the occurrences. The models were tested using the AUC approach of Phillips et al. [81] where an AUC score of 0.5 indicates that the estimated range is not better than random and a score of 1.0 indicates perfect estimation. The 100 replicate SDMs were combined to create presence-absence maps using the minimum training presence logistic threshold and the ten-percentile training omission threshold methods as recommended by Liu et al. [84].

Paleophylogeographic species distribution models (PPGMs)

To estimate geographic ranges of the species in the geological past, we used the paleophyogeographic approach described above, which incorporates phylogenetic correction for niche evolution and the continuous series of paleoclimatic models [23]. Potential niches were estimated from rectilinear BIOCLIM envelopes around modern E-space species occurrence points. The BIOCLIM approach is preferred for paleogeographic modeling because (i) it makes fewer assumptions about the constancy of the species-climate relationship than Maxent, despite the fact that it tends to over-predict ranges in the modern world (the over-prediction is because it tends not to identify range limits defined by biotic interactions) [18], [23], [49], [83]; and (ii) because its rectilinear shape has the potential to include non-analogue climates that are part of the fundamental niche, but which are not available in modern climate and thus not part of more restrictive niche envelopes. Importantly, rectilinear BIOCLIM climate envelopes are unbiased by correlations between climate variables, which is critical because the correlations between variables are not expected to be the same in the past as they are today. The full potential niche and one that excludes the outlying 10% of the E-space points were estimated. The suitability of the two estimates was assessed by their ability to correctly estimate the modern geographic distributions. The reduced envelopes were found to be too restrictive for many of the turtle species.

The modern potential niches were projected onto a phylogenetic tree of the 59 species (see below) to estimate the phylogenetic component of niche change (see [23] for complete description). Ancestral niches at the nodes of the tree and along its branches can be reconstructed from a generalized linear model (GLM) regression. The correlation structure we used in the GLM method assumes that the climatic niches evolved according to a Brownian motion at a stochastically constant rate. The regression was used to estimate the niche boundaries of each species at 4 ky intervals along their respective branches. These phylogenetically corrected potential niches were then used to estimate each species' paleogeographic range at each 4 ky interval. We refer to these estimated ranges as paleophylogeographic models (PPGMs). Each PPGM was masked using Level-1 watersheds from Pfafstetter System, which are large-scale river drainages like the Mississippi Basin (Figure 2D), in order to exclude spurious points in the reconstructed ranges that lie beyond the geographic barriers that make them inaccessible to the species. The 4 ky PPGM slices were compiled into an animated map that shows the changing potential range for each species through the last three glacial-interglacial cycles.

Figure 2
Predicted current species richness according to 100% and 90% environmental envelopes (A) as well as historic fluctuations as projected for the last glacial maximum 21 ky BP (B) and the last interglacial (C) according to palaeophylogeographic models of ...

To partition the geographic responses of each species into adaptive and purely geographic components, we subtracted the phylogenetically corrected PPGMs from uncorrected PPGMs [23]. The phylogenetic (adaptive) component is the difference in geographic area and centroid between the two. The purely geographic component is the residual of the uncorrected PPGMs. The rate of geographic response per degree change in mean annual temperature (BIO1) was calculated as a function of total area of the geographic range and position of the range's geographic centroid using simple linear regression [23] and, for visualization, kernel density estimation using kde2d function of the MASS package for R 15.0 [85].


We constructed a time-calibrated molecular phylogeny for the 59 chelonian species using a supermatrix approach. Data from four mitochondrial genes (COI, control region, cytb, ND4) and three nuclear genes (R35, RAG-1, TB29) were downloaded from GenBank. The number of loci varied across the taxa, but even with missing data the concatenated data are still expected to provide a reasonable estimate of the phylogeny [86]. GenBank accession numbers for the sequences are reported in Appendix S1 in Material S1. Each gene was aligned separately using the MUSCLE algorithm (MUSCLE version 3.6, [87]) and checked manually. Ambiguously aligned or randomly similar sequence sections were detected with ALISCORE [88] and removed. The final, concatenated data set contained 5742 bp (including gaps). Parameters for the substitution model (GTR+G) were estimated separately for each gene by partitioning the data set during subsequent analyses. A Yule speciation process was used to estimate the prior probability of the tree. Divergence times were estimated with BEAST 1.6.2 with a relaxed molecular clock model with rates with an uncorrelated lognormal distribution [89][90]. Two fossil Trachemys species, T. idahoensis (4 Ma) and T. inflata (7 Ma) were used as calibration points for the most recent common ancestors (MRCA) of T. scripta/T. gaigeae and of the North American clade/Tropical clade of Trachemys, respectively [91]. Glyptemys valentinensis (14.5 Ma) was used to calibrate the MRCA of the genus Glyptemys; [92], and Kinosternon pojoaque (13.3 Ma) was used to calibrate the MRCA of K. arizonense, K. durangoense and K. hirtipes; [93]) For T. idahoensis, a uniform prior distribution with the lower bound at 4 Ma and the upper bound at 7 Ma (i.e. the lower bound of the older fossil species T. inflata) was used (cf. [92]). We used single-sided hard bounds for the minimum age of divergence and an exponential distribution as a soft bounds for the maximum age. Every 1000th trees were sampled from three independent runs with 107 generations. The initial 5% of each run was discarded as burn-in, resulting in a total of 28,500 trees. Maximum likelihood bootstrap values (1,000 replicates) were calculated using RAxML v.7.2.X [94].

Dispersal limits and accessible climate

Drainage basin boundaries have been shown to be important barriers to turtle dispersal [66], [95]. Geographic data on basin boundaries were obtained from the USGS HYDRO1k dataset ( Level-2 drainage basin units were used as to define each species' available climate space for Maxent training [96] and for calculation of MESS scores (see below). Level-1 units (Figure 2D) were used to mask PPGMs. The masking allows us to incorporate the concept of accessibility into both the niche modeling and the geographic distribution modeling (see discussion in Introduction).

Multivariate environmental similarity surface (MESS) scores

Multivariate environmental similarity surface (MESS) scores [51] were used to identify non-analogue paleoclimates and to describe the position of each species niche within its available climate space. Negative MESS scores identify past occurrences that fall outside the available climate space on one or more of the climate variables (i.e., non-analogue occurrences). Identification of non-analogue occurrences is important for distinguishing opening of potential niche space from evolution of the fundamental niche due to changing climate [97].

A MESS score is assigned to each point in the E-space based on where it falls with respect to the boundaries of a set of reference points. If the point falls within the boundaries, then a positive MESS score is assigned with a range from 0 to 100 depending on how close the point is to the center of the bounded area (a score of 100 indicates that the point lies in the middle of the bounded area). If the point lies outside the bounded area, then a negative MESS score is assigned with a range from 0 to -∞ as the Euclidean distance between it and the nearest boundary divided by the width of the bounded area. Thus if the distance between the point and the nearest boundary is equal to the width of the bounded area, then the MESS score is −100. For more details see Elith et al. [52].

Four kinds of MESS map were computed. A map of MESS scores was produced for each fossil occurrence with respect to (a) the boundaries of the species modern potential niche space and (b) for the boundaries of the combined available climate space of all the species using first the combined set of climate variables and then for each individual climate variable.

We used the MESS scores to deal with the uncertainty of the ages of the older fossils. We found the time slice in which the paleoclimate was the most compatible with the fossil occurrence as the one whose potential niche MESS score was the most positive (see below). Thus, we conservatively assumed that the best estimate of the age of the fossil is the time slice where the paleoclimate is most compatible with its geographic occurrence. Because we interpret negative MESS scores as indicators that potential climate space has been opened, this approach minimizes ad hoc assumptions that the fundamental niche has changed. Note that for two fossils we found two time slices with identical MESS scores so we included both slices as being equally likely.

For some time slices, the PPGMs of some species contracted to nothing, suggesting that the paleoclimate provided no suitable climatic habitat for the species and that it could easily have become extinct. We know with hindsight that the species are, in fact, still living, so for those species we arbitrarily assigned them a geographic range of the 5 equal-area grid cells (250 km2) whose potential niche MESS scores were the least negative (SIM5). We assumed that the SIM5 range was the most parsimonious refugium for the species (these arbitrarily assigned ranges are highlighted in the video montage maps of the PPGMs). We used potential niche MESS scores for the SIM5 as an estimate of how far a species can be pushed outside its modern climate limits and still survive.

We used the potential niche MESS scores computed with all climate variables and single variables for the past distributions relative to their modern potential niche in E-space as a quantitative measure of how much the potential niches have changed over the last 320 ky. We visualized the changes using the function of the sm package [98] in R 15.0.

Alternative ways to measure niche change exist (e.g. [54], [99][101]). These include the amount of overlap between niches in E-space as a parameter that ranging from zero to one (e.g., Schoener's D), where zero means no overlap and one means complete overlap [100][101]. We preferred MESS scores because these alternatives do not provide information about the scale of difference between niches that do not overlap. Furthermore, they cannot be used to compare single point occurrences, e.g. fossils, to a niche distribution.

Comparisons of paleophylogeographic models (PPGM) with molecular phylogeographic data

We compared our paleophylogeographic range maps to published phylogeographic scenarios derived from molecular and morphologic data for the 34 species with such data (the phylogeographic data are summarized in the left margin of the corresponding animated paleophylogeographic maps in Appendix S4 in Material S1). We compared the number of genetically distinct intraspecific phylogroups (based on molecular data) and their spatial distribution (e.g. phylogenetic lineages identified via allele frequencies, haplotypes and/or morphometric data) to the PPGM range maps to determine the concordance of these data with past range fragmentations and contractions. We used a simple binary scoring to summarize whether the molecular phylogeographic data were consistent (1) or inconsistent (X) with our paleophylogeographic maps in the number of genetically distinct phylogroups and the location of refugia. We chose to not statically quantify the geographic overlap between molecularly predicted refugia and our paleophylogeographic reconstructions [102] because of lack of geographic specificity of most phylogeographic data and because the conclusions we draw are sufficiently supported by our simple classification.

Physiological data

Physiological parameters that might restrict the distribution of chelonians were compared to the modern SDM and Quaternary PPGM reconstructions to determine whether the realized niches in E-space have limits defined by the animals' physiology. Data on the following parameters were compiled from 226 field and lab experimental studies: maximum temperature tolerated by the species (critical thermal maximum, CTmax), minimum temperature tolerated by the species (critical thermal minimum. CTmin), and temperature range at which incubation occurs. We found data for 48 of the 59 species (Appendix S5 in Material S1). Maximum temperature of the warmest month (BIO5) was compared to CTmax, the minimum temperature of the coldest month (BIO6) was compared to CTmin, and the mean temperature of the warmest quarter (BIO10) was compared to incubation temperature range. Comparisons were made for each species' modern available climate space and its realized niche.



A single, well-supported tree resulted from the multi-gene Bayesian analysis (Appendix S6 in Material S1). All higher taxa (families, subfamilies and genera) involved were found to be monophyletic, except the genera Sternotherus and Trachemys (see discussion). The time since common ancestry of several genera were estimated to be very shallow (e.g. Chrysemys, Graptemys and Pseudemys) and most species are estimated to have originated in the last 2 my. The genus Apalone is estimated to have diverged from other taxa about 10 Ma.

Modern SDM geographic range models

The Maxent SDM estimates of modern geographic ranges generally agreed well with the actual distributions (AUC values ranged from 0.740 to 0.972 with a mean 0.891 where 1.0 is perfect match and 0.5 is no better than random) (Appendix S7 in Material S1). All but two species, Kinosternon durangoense and Apalone spinifer, had modeled ranges that closely matched the real ranges with AUC>0.8. One of those, K. durangoense, a Mexican endemic, is very poorly studied and only a few occurrence points are known, thus contributing to uncertainty about its range. The other, A. spinifera, a very widespread species that occurs in many climates, has a broad general realized niche that is compatible with a much larger geographic range than it actually inhabits.

The climate variables that placed the strongest limits on the niches and thus the SDM ranges were mean temperature of warmest quarter (BIO10), minimum temperature of coldest month (BIO6), precipitation of driest quarter (BIO17), and precipitation of coldest quarter (BIO19). Most of the eight climate variables that we used for the SDMs were similarly relevant for each species, with standard deviations of there contributions ranging from 7.1 to 15.6% across species (a low standard deviation means that the variable has a similar effect on all the species). Nevertheless, the variation in how much some variables contribute to one species compared to another is high for some variables, indicating that the species differ in the variables that define the limits of their climate niches. For example, precipitation in the coldest quarter, BIO19, whose contribution ranges from 0.1% in Pseudomys alabamensis, a species with a tiny distribution along a small segment of the Gulf Coast, to 76% in Graptemys oculifera, a species that also has a small distribution but one that extends northward into the continental interior in Alabama.

All Maxent results are provided in Appendix S7 in Material S1. Maps showing the SDM geographic range models for each species are provided in Appendix S4 in Material S1.

Paleophylogeographic models of Quaternary geographic ranges (PPGMs)

Paleophylogeographic models (PPGMs) are presented as animated maps in Appendix S4 in Material S1 to show the modeled changing geographic ranges of each of the 59 species through the last three glacial-interglacial cycles (320 Ka). Each map shows the PPGM models of the expansions, contractions, and fragmentation of the paleogeographic range along with the fossil occurrences, molecular phylogeographic data, and map of parts of the continent with non-analogue climates. The paleogeographic range models take into account evolutionary changes in the species' climatic niche as estimated from the phylogenetic tree presented above. Each animated map consists of 80 frames separated by 4 ky intervals. The global oxygen isotope curve is shown to the right of each map with a moving index marker to show the age of each frame. These animated maps are key to understanding our conclusions.

Summary maps showing species diversity at three key times are shown in Figure 2. Modern species diversity calculated from the SDM geographic range models is shown in 2A. Paleo species diversity calculated from the PPGM paleogeographic range models is shown in 2B for the last glacial maximum (20 Ka) and in 2C for the last interglacial maximum (140 Ka). The regions that are today covered by temperate forests, temperate grasslands, deserts, and lake systems were regions where geographic ranges changed the most over the last 320 ky. In contrast, the regions along the Pacific Coast, the mixed mountain highland biomes, and the tropics were the regions where the ranges changed the least (Appendix S4 in Material S1, Figure 2).

There were strong relationships between change in mean annual temperature (PCMAT) and change of the geographic centers of species' ranges (PCGC) (R2 = 0.904, p<0.0001, PCGC = 28.66+73.17*PCMAT) and the total areas of their geographic range (PCRS) (R2 = 0.796, p<0.0001, PCRS = 24.50+2135.68*PCMAT) based on the PPGM models (Figure 3). In other words, the average geographic range of a species expanded or contracted by more than 2,000 km2 for every degree that the mean annual temperature increased or decreased. When the phylogenetic correction was omitted, the relationship between change mean annual temperature and purely climate driven change in the geographic center of the range was about the same (R2 = 0.904, p<0.0001, PCGC = 8.32+78.74*PCMAT) and the relationship with change in area was stronger (R2 = 0.837, p<0.0001, PCRS = −42.09+2323.9*PCMAT) indicating that adaptation has contributed by only a small amount to the responses to climate change over the last 320 Ka. This conclusion was confirmed by modeling the purely phylogenetic component of change, which had no significant relationship between PCMAT and PCGC (R2 = 0.016, p>0.05) or between PCMAT and PCRS (R2 = 0.019, p>0.05). Geographic tracking has been by far the dominant mechanism by which chelonians have responded to changing climates over the last 320 Ka. The pace of phylogenetic diversification of their niches has been far slower than the 100 ky glacial-interglacial cycles of the Quaternary.

The greatest contraction of geographic ranges in PPGMs occurred during glacial periods. In fact, there were two or more 4-ky time slices when there was no suitable climate for 33 species: min = 2 slices for Pseudemys gorzugi; max = 96 slices for Graptemys pearlensis; median = 45 slices for the 33 species (see Appendix S4 in Material S1). Most of the species that had such extreme geographic bottlenecks in our models have very small modern ranges or are found close to the Gulf of Mexico.

Paleogeographic ranges compared to phylogeographic data and fossil occurrences

One of the main goals of this study was to determine whether the potential climatic niches of North American turtles are adequate proxies for their fundamental niche by testing how well paleogeographic ranges derived from them agree with phylogeographic data and fossil occurrences. Our paleophylogeographic range models were usually successful at predicting modern patterns of within-species genetic differentiation. The patterns of range fragmentation during glacial cycles corresponded to the spatial distribution of living subspecies (11 out of 14 cases) and an even higher proportion of the fragmentation patterns corresponded to the modern distribution of genetically distinct phylogeographic lineages (21 of 23 species) (Table 1, Appendix S4 in Material S1). Species with deep genetic splits between subspecies or populations often experienced long periods of climatically driven range fragmentation, in many cases cyclically repeated allopatric separation (e.g., Apalone mutica [103], A. spinifera [104], Kinosternon subrubrum [105]), whereas species whose past ranges were not fragmented by climate cycles show little evidence of genetic differentiation (e.g. Actinemys marmorata [106]; Table 1, Appendix S4 in Material S1). In rare cases, species with marked range fragmentation in the past have little genetic structure today (e.g. Chelydra serpentina [106]). Species whose paleogeographic ranges contracted to nothing or nearly nothing almost always have very low genetic diversity and small geographic ranges today, suggesting that they experienced climate-driven geographic bottlenecks in the past (e.g., Chrysemys dorsalis [107], Glyptemys insculpta [108], Glyptemys muhlenbergii [109], Graptemys gibbonsi [110]; Table 1, Appendix S4 in Material S1). Morphologic data generally corresponded with these results. Intraspecific morphological diversities were smallest in those species with small distributions in the present or past or where large range changes likely caused bottleneck effects.

Most of the fossils fell within their potential niches (87%), indicating strong niche stability (Table 2). E-space plots of the fossils, the available climate, the realized niches, and the potential niches for each species are shown in Appendix S2 in Material S1. Of the 13% that fell outside the modern available niche space, most did so on mean temperature of the warmest quarter (BIO10; 66.7%) and mean diurnal temperature range (BIO2; 9.7%) (Table 3, Figure 4, Appendix S3 in Material S1). In all, fossils were outliers on only seven of the 19 bioclimatic variables (BIO1, 2, 5, 6, 9, 10, 11). We interpreted this to mean that niches were conserved for the other 12 bioclimatic variables, and we interpret the 87% agreement between fossils and potential niches to mean that niche conservatism is strong (after taking into account phylogenetic evolution).

Table 2
Summary statistics for fossil occurrences in Nearctic chelonians in terms of availability and niche position relative to the species' modern niches (for more details see Appendix S3 in Material S1).
Table 3
Multivariate environmental similarity scores (MESS) for fossils that fell outside the potential distribution of their species. Negative scores indicate the distance outside the potential niche as a percentage of the niche's size.

Physiological constraints on species ranges and niches

Physiological data were generally congruent with potential niches. The maximum temperature that a species is able to tolerate (CTmax) was greater than the maximum warmest temperature (BIO5) in the past and present potential niches of 100% of the species for which CTmax was known (27 out of 59 species [111][114]). Chrysemys picta and Actinemys marmorata were the only species where the maximum temperature of their potential niche came to the limits of CTmax (Appendices S2 and S5 in Material S1). This suggests that most turtle species can tolerate higher maximum temperatures than they have been subjected to in the present or past. Incubation temperature range was also generally congruent with past and present potential niches. Eight species with fossils had incubation temperature data. Of those only one had incubation ranges that were incongruent with its potential niche.

Conversely, the minimum temperature that a species can tolerate (CTmin) was strongly incongruent with the potential niche data. CTmin data and fossil data were available for eight out of 59 species, seven of which had potential niche minimum temperatures (BIO6) that were lower in the past or present than CTmin, suggesting either that the species were able to behaviorally insulate themselves from winter colds or that CTmin has changed over time, or that the potential niches in the past were limited by CTmin.


Our paleophylogeographic models illustrate the dramatic range changes that most turtles have experienced during the last three glacial-interglacial cycles (Appendix S4 in Material S1). The geographic ranges of most species fragmented and contracted into one or more southern refugia, often in peninsular Florida, coastal Texas and northern Mexico. Most of the Mexican and Central American species experienced very little change in the geographic centers of their ranges, but the areas of their ranges often contracted. Some local populations in these southern species may have persisted for geologically long time periods in the same place, which is not the case for the northern species. The PPGM projections agree well with molecular phylogeographic patterns (82.4% of 17 species are congruent), genetic subgroupings (91.3% of 23 species are congruent), and subspecies-level differentiation (78.6% of 14 species are congruent; Table 1). Genetic diversities are highest in species for which our PPGM models indicate either stable or large potential distributions, while low genetic diversities were mostly evident in species for which our PPGM maps show significant range changes or severely restricted distributions. Fossil occurrences are generally congruent with the geographic ranges and potential niche reconstructions of the PPGMs, which suggests that climate tolerances of these species have not changed substantially over the last 320 Ka. However, there are exceptions that do indicate some niche changes (Table 2).


Our phylogenetic tree is largely congruent with previous phylogenetic work (for a synopsis see Iverson et al. [115]). Most discrepancies were found within the Kinosternidae and some genera of the Deirochelinae. Regarding the first, the position of Sternotherus depressus leads to a paraphyly of Sternotherus in respect to Kinosternon. This might be an artifact of the sparse and often non-homologous data available for the members of these two genera and is also reflected by generally low support values in this clade. However, there is currently no comprehensive molecular phylogeny of this group available for verification. Although intrageneric relationships of Graptemys and Pseudemys in our results slightly differ to previous hypotheses, the posterior probabilities strongly support our topology (Appendix S6 in Material S1). The major ambiguity is the polyphyly of Trachemys due to T. nebulosa, which is the sister species of Malaclemys according to our results, yet this has poor support.

Paleophylogeographic reconstructions agree with molecular phylogeography and subspecies taxonomy

Our animated paleophylogeographic range maps (Appendix S4 in Material S1) are strongly congruent with biogeographic patterns inferred from population genetic studies and phylogenies. The number and location of past refugia in our models mostly agree with reconstructions derived from population genetic and phylogeographic analyses. Refugia in the PPGMs correspond to deep genetic splits (e.g., subspecies) in most species (c.f., Apalone mutica, Chrysemys picta, Trachemys scripta, Rhinoclemmys pulcherrima, Terrapene nelson, Trachemys nebulosa and Trachemys venusta), though minor geographic entities were sometimes less well resolved via PPGMs (c.f., Macrochelys temminckii and Malaclemys terrapin). Our modeling exercises suggest that several species had no suitable climate space for at least some period during the past three interglacials. This is unlikely to be a methodological artifact resulting from masking the PPGMs by watershed, because most of the Level-1 watersheds are oriented such that a species could track climate for very long north and south distances.

It is instructive to look at the details of species for which there appears to have been no suitable climate during Pleistocene cold phases because the details of each species are different [116]. The freshwater emydid Emydoidea blandingii is one. This species inhabits the Missouri river system as well as rivers along the Upper Mississippi catchment, meaning that the entire Mississippi basin was treated as available space in our models. The species should be easily able to track its niche southward, if suitable climate space had been available. Our SIM5 points (the five 50 km2 points that are climatically most compatible with the potential niche of the species) suggest three main areas where this species may have found refugia that were close to the climate it is known to tolerate: on the great plains of Kansas and Nebraska, in the Midwest south of the Great Lakes, or on the mid-Atlantic seaboard (Appendix S4 in Material S1). A fossil occurrence in the Great Plains confirms the first of these as a possibility, and a genetic study by Mockford et al. [117] is also consistent with this picture. Using mtDNA, these authors argued for two past barriers of gene flow that match the breaks in our PPGM models.

Glyptemys insculpta, a northerly distributed semi-freshwater species is another example. Amato et al. [108] conducted range-wide genetic analyses on this species and detected very weak differentiation based on mitochondrial genes. The authors suggest a combination of the origin of a single glacial refugium, as well as a selective sweep of the mtDNA genome as the most likely explanation for their findings. A selective sweep was supported through a comparative study finding higher differentiation in microsatellites [118] leading to a mismatch in genetic structuring between the two marker systems. Our PPGM models of G. insculpta's geographic history suggest the species had two refugia separated by the Appalachians. Fossil occurrences confirm the presence of a population during the Late Pleistocene – Early Holocene in the east of the Appalachians around the Alabama/Georgia region (cf. Amato et al. [108] and references therein). High differentiation of populations suggested by microsatellite analyses [118] makes the existence of further refugia in the west likely. Increased selective pressure within the mtDNA could lead to a dominant genotype of the eastern refugial clade, leading to a selective sweep also in individuals originated from the western refugium after postglacial expansion and secondary contact.

A similar picture becomes evident when looking on southerly distributed species, but for different reasons. In contrast to northerly distributed species, the Gulf of Mexico represents an absolute barrier for southward movements in species with ranges adjacent to the Gulf coast. Several species – especially those with very small range extents such as some Graptemys species – may have found non analogous climate space during the past within their restricted geographic distributions. Because all species studied here have survived the climatic variations during the past, there are two main explanations: (1) species have undergone shifts in their realized niche which would give insights into the dynamics of their respective fundamental niche, or alternatively (2) microhabitat features not covered by the comparatively coarse grained resolution of our modeling act as microrefugia as proposed by and Rull [120], which remain unconsidered in our framework. The latter, for example, may be an explanation for the discrepancy between the distribution of currently observed genetic structure and PPGM results in Macrochelys temminckii. An analytical comparison of both genetic analyses and PPGMs would shed light into this issue.

The history of intraspecific diversity

One major implication of the paleophylogeographic approach is that we can disentangle the geographic history of species that have experienced loss of genetic information due to genetic drift. Genetic drift (sensu Wright [121]) is a phenomenon that can lead to positive or negative impacts on intraspecific diversity, in which genetic information can get lost and/or recessive information (alleles or haplotypes) can become dominant. Intraspecific (genetic or morphologic) variability differs strongly among species as well as within populations of each species and highly depends on the species' history. In the field of biogeography, we have to distinguish two main processes severely influencing the intraspecific variability: (i) the positive impact of genetic drift may be the consequence of long-lasting geographic isolation and long-term restriction to geographically small retreats, and (ii) a negative impact (like losses of genetic information) in the wake of repeated expansion-retraction dynamics over large geographical space.

When comparing the potential distribution projections from SDMs and PPGMs with molecular data, we find congruencies between range shifts during the past with species' intraspecific diversities in the majority of the analyzed cases. Low genetic diversities are detectable in those species showing the most extensive range contractions. Small and isolated refugia often lead to small population sizes. Individuals from such populations suffer much stronger under environmental stochasticity, population fluctuations and thus repeated population bottlenecks [122]. Genetic data reveal strong population bottlenecks for various turtle and tortoise species from the southern US and Mexico [104], [123][124], except for species where different refugial ranges became evident during modeling (e.g. [104], [106], [125], but see [126]). For example, SDMs obtained for Emydoidea blandingii, Glyptemys insculpta and Glyptemys muhlenbergii highlight strongly restricted and small refugia over a long time period and,molecular analyses reveal a negative impact on intraspecific diversity. In contrast, SDMs indicate large and rather stable refugia for Apalone spinifera, which might had consequences on the survival of population networks and led to the maintenance of high genetic diversities over time, as revealed by molecular data [105]. The second negative effect of genetic drift can be caused by repeated expansion-retraction dynamics of a species' distributional area in the wake of tracking its climatic niche. These range shifts (e.g. the (re)colonization of (new) suitable habitats) may imply severe population bottlenecks in the wake of small founder populations and in consequence might lead to the elimination of alleles [3], [127][128]. Species with restricted dispersal (e.g. stepping-stone wise movement) are more affected by this phenomenon than species showing a phalanx-wise migration mode (see [3]). Interestingly, several species' PPGM ranges dwindle to nothing during cooling events (i.e., glacial periods) in the PPGM framework. This highlights the species-specific niche dynamics or, alternatively, scale effects not covered by our framework such as the presence of microrefugia [119], [120].

Nearctic chelonian thermal biology

Correspondence between physiological variables and the potential niche described by SDMs and PPGMs vary, where CTmax is congruent and CTmin is incongruent. The potential niches of all species compared are in congruence with CTmax. This suggests CTmax is represented accurately in the SDMs and that it is a relatively conserved variable. On the other hand, CTmin is incongruent with the potential niches. The comparison between the minimum temperature of the coldest month (BIO6) with CTmin was intended to explain the extreme surface temperatures during hibernation within the modern potential range. However, many Nearctic chelonians are aquatic and aquatic species usually spend their hibernation period submerged in water at 4°C (e.g. [129]). This hibernation strategy enables the habitation of environments that only temporarily provide suitable conditions. In addition, hatchlings of numerous species are reported to hibernate within the nest and withstand temperatures well below the critical thermal minima reported for adults (e.g. [129]).

Niche dynamics revealed by fossils and SIM5

Analyzing the relative niche position of turtle fossils unexpectedly reveals that the fossils that deviate from the species potential niche do so mostly in mean temperature of the warmest quarter (BIO10). Within the deviant fossils and their associated species, this suggests large niche shifts of up to one niche breadth in multivariate MESS analyses using the comprehensive set of 19 bioclimatic variables, as well as in univariate MESS analyses (Figure 4). Comparing our fossil and SIM5 data with our PPGM ranges, we detected a multidimensional niche shift in high proportions of the species (Figure 2, Tables 2, ,3).3). Notably, the mean temperature of the warmest quarter (BIO 10) shows a much higher plasticity of pivotal incubation temperatures than previously assumed. This result is somewhat unexpected as the majority (42 of 59 species are documented in the literature) of the Nearctic chelonians have temperature-dependent sex determination (TSD), in which ambient temperature during egg incubation in the summer months determines sex, and only a few (four species have documented cases) have genotypic sex determination (GSD) (see Appendix S5 in Material S1). In TSD species, the sex is environmentally determined wherein, depending on the TSD mode, a balanced sex ratio among hatchlings is only possible within a relatively narrow temperature range [64], [130]. Incubation temperatures exceeding this range in either direction, colder or warmer, may cause the development of only one sex and ultimately skewed sex ratios in populations. In previous studies, authors assume that sex-determining mechanisms may strongly affect the survival probability of a species in a geological time frame [131]. Recently, there is a considerable debate about how the impact of human-induced climate changes might affect species with TSD [132][133], and some studies underline that climate change might jeopardize sex ratios and therefore cause severe population decline and extinction [134]. On the other hand, Silber et al. [135] provided evidence that a higher proportion of TSD species compared with GSD species may have survived historic climate fluctuations suggesting some advantages of TSD over GSD, wherein Valenzuela & Adams [136] provided evidence that chelonians might have adapted to prehistoric climate change by shifting from TSD to GSD. These results suggest that the species' niches may be much more flexible than previously assumed.

Revolution in the application of distribution models

The paleophylogeographic approach can be of high relevance in understanding the Quaternary background of species, including past range shifts with subsequent differentiation (e.g. the evolution of Evolutionary Significant Units, [137]) and genetic bottlenecks. All of our comparisons with PPGMs indicate strong congruencies for most species analyzed here for physiological, morphologic and molecular characters. Hereby, the data show strong overlap in past differentiation processes among/within taxa and the evolution of genetic diversity, as well as novel insights into past niche dynamics. These results underline the applicability of PPGMs as a valuable tool to model past range shifts and to picture potential processes standing in concert with long-term isolation and severe range contraction (e.g. population fluctuations and subsequent bottlenecks). Thus, species distribution models have the potential to shift from classical biogeographic range-shift analyses (from past to future) to conservation biology and the characterization of species, and could be a useful tool to predict intraspecific structures (such as differentiation, lineages and diversities) - without expensive molecular analyses. In regard to the 59 turtles and tortoise species, PPGMs are especially of high relevance since the modern distributions for these taxa will turn into non-analogous climatic conditions within the next few decades [59], [138][139], which might ultimately lead to cryptic biodiversity loss as suggested for numerous taxa [140].

Supporting Information

Material S1

Zip file containing Appendices S1 through S7. Appendix S1: Genbank accession numbers of sequences of four mitochondrial genes (COI, control region, cytb, ND4) and three nuclear genes (R35, RAG-1, TB29) used for phylogenetic analyses; Appendix S2: Plots showing the bioclimatic conditions within each species' currently available climate space, its potential niche as well as its realized niche both in terms of species' modern and fossil occurrences; Appendix S3: Fossils of Nearctic chelonians compiled from Holman [61] and supplemented with data obtained through the Paleo Biology Database (available through Negative MESS scores indicate conditions within the species' current realized/potential niche and positive MESS scores indicate bioclimatic conditions outside of the current realized/potential niche; Appendix S4: Videos showing the results of the Maxent and PPGM analyses, i.e. distribution changes through time in comparison to current distributions according to Maxent, non-analogous climate and information on subspecies and phylogenetic lineages; Appendix S5: Comprehensive summary of data on thermal biology of 48 Nearctic chelonian species; Appendix S6: Dated molecular phylogeny for all 59 extant species of Nearctic turtles inferred by Bayesian analyses; Appendix S7: Summary statistics of each 100 Maxent models per species.



We are grateful to U. Bott for valuable help during the preparation of the manuscript and to A. Donath for many discussions how to set up the analyses pipelines. We acknowledge the Paleobiology Database ( for collecting and archiving data and would like to thank J.B. Iverson, J.A. Kimberling, A.R. Kiester, L.E. Hughes and J. Nicolello for making high quality species records available through the World Turtle Database. We are grateful to D. Nogués-Bravo and J. Stewart as well as anonymous referees for providing comments and excellent suggestions to improve an earlier version of this manuscript.

Funding Statement

This work was supported by the German Academic Exchange Service (Deutscher Akademischer Austauschdienst, DAAD; PostDoc fellowship JCH), through a postdoctoral fellowship to AML at the National Institute for Mathematical and Biological Synthesis (NIMBioS), which is sponsored by the National Science Foundation, the U.S. Department of Homeland Security, and the U.S. Department of Agriculture through National Science Foundation (NSF) Award #EF-0832858, with additional support from The University of Tennessee, Knoxville, and through a National Science Foundation Grant EAR-0843935 to PDP. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Bellard C, Bertelsmeier C, Leadley P, Thuiller W, Courchamp F (2012) Impacts of climate change on the future of biodiversity. Ecol Lett 15: 365–377 [PubMed]
2. Graham RW, Lundelius EL Jr, Graham MA, Schroeder EK, Toomey RS III, et al. (1996) Spatial response of mammals to Late Quaternary environmental fluctuations. Science 272: 1601–1606 [PubMed]
3. Hewitt G (1996) Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc 58: 247–276
4. Riddle BR (1998) The historical assembly of continental biotas: late Quaternary range-shifting, areas of endemism, and biogeographic structure in the North American Mammal fauna. Ecogeography 21: 437–446
5. Avise JC (2000) Phylogeography: the History and Formation of Species. Boston: Harvard University Press. 447 p.
6. Hewitt GM (2000) The genetic legacy of the Quaternary ice ages. Nature 405: 907–913 [PubMed]
7. Sommer RS, Zachos FE (2009) Fossil evidence and phylogeography of temperate species: “glacial refugia” and post-glacial recolonization. J Biogeogr 36: 2013–2020
8. Swenson NG, Howard DJ (2005) Clustering of contact zones, hybrid zones, and phylogeographic breaks in North America. Am Nat 166: 581–591 [PubMed]
9. Shafer ABA, Cullingham CI, Cote SD, Coltman DW (2010) Of glaciers and refugia: a decade of study sheds new light on the phylogeography of northwestern North America. Mol Ecol 19: 4589–4621 [PubMed]
10. Comes HP, Kadereit JW (1998) The effect of Quaternary climatic changes on plant distribution and evolution. Trends Plant Sci 3: 432–438
11. Habel JC, Schmitt T, Müller P (2005) The fourth paradigm pattern of postglacial range expansion of European terrestrial species: the phylogeography of the Marbled White butterfly (Satyrinae, Lepidoptera). J Biogeogr 32: 1489–1497
12. Hewitt GM (2001) Speciation, hybrid zones and phylogeography – or seeing genes in space and time. Mol Ecol 10: 537–549 [PubMed]
13. Hewitt GM (2004) Genetic consequences of climatic oscillation in the Quaternary. Phil Trans R Soc B 359: 183–195 [PMC free article] [PubMed]
14. Gómez A, Lunt DH (2007) Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. In: Weiss S, Ferrand N, editors. Phylogeography of southern European refugia. Heidelberg: Springer. pp. 15–188.
15. Hewitt GM (2011) Mediterranean Peninsulas: The Evolution of Hotspots. In: Zachos FE, Habel JC, editors. Biodiversity Hotspots, Distribution and Protection of Conservation Priority Areas. Heidelberg: Springer. pp. 123–147.
16. Hampe A, Petit RJ (2005) Conserving biodiversity under climate change: the rear edge matters. Ecol Lett 8: 461–467 [PubMed]
17. Darwin CR (1859) On the origin of species by means of natural selection, or the preservation of favoured races in the struggle for life (1st Ed.). London: John Murray. 502 p.
18. Svenning J-C, Fløjgaard C, Marske KA, Nógues-Bravo D, Normand S (2011) Applications of species distribution modeling to paleobiology. Quaternary Sci Rev 30: 2930–2947
19. Stewart JR, Lister AM, Barnes I, Dalén L (2010) Refugia revisited: individualistic responses of species in space and time. Proc Royal Soc B-Biol Sci 277: 661–671 [PMC free article] [PubMed]
20. Peterson AT (2011) Ecological niche conservatism: a time-structured review of evidence. J Biogeogr 38: 817–827
21. Wiens JJ (2011) The niche, biogeography and species interactions. Philos T Roy Soc B 366: 2336–2350 [PMC free article] [PubMed]
22. Espíndola A, Pellissier L, Maiorano L, Hordijk W, Guisan A, et al. (2012) Predicting present and future intra-specific genetic structure through niche hindcasting across 24 millennia. Ecol Lett 15: 649–657 [PubMed]
23. Roberts DR, Hamann A (2012) Predicting potential climate change impacts with bioclimate envelope models: a palaeoecological perspective. Glob Ecol Biogeogr 21: 121–133
24. Guisan A, Thuiller W (2005) Predicting species distribution: offering more than simple habitat models. Ecol Lett 8: 993–1009
25. Chen I-C, Hill JK, Ohlemüller R, Roy DB, Thomas CD (2011) Rapid range shifts of species associated with high levels of climate warming. Science 333: 1024–1026 [PubMed]
26. Bradshaw WE, Holzapfel CM (2006) Evolutionary response to rapid climate change. Science 312: 1477–1478 [PubMed]
27. Parmesan C (2006) Ecological and evolutionary responses to recent climate change. Ann Rev Ecol Evol Syst 37: 637–669
28. Vieites DR, Nieto-Román S, Wake DB (2009) Reconstruction of the climate envelopes of salamanders and their evolution through time. Proc Nat Acad Sci USA 109: 19715–19722 [PubMed]
29. Lawing AM, Polly PD (2011) Pleistocene climate, phylogeny, and climate envelope models: an integrative approach to better understand species' response to climate change. PLoS ONE 6: e28554. [PMC free article] [PubMed]
30. Varela S, Lobo JM, Hortal J (2011) Using species distribution models in paleobiogeography: A matter of data, predictors and concepts. Palaeogeogr Palaeoclim Palaeoecol 310: 451–463
31. Guisan A, Zimmermann N (2000) Predictive habitat distribution models in ecology. Ecol Model 135: 147–186
32. Elith J, Leathwick JR (2009) Species distribution models: ecological explanation and prediction across space and time. Annu Rev Ecol Evol Syst 40: 677–697
33. Elton C (1927) Animal Ecology. London: Sedgwick and Jackson. 207 p.
34. Grinnell J (1917) The niche-relationships of the California Thrasher. Auk 34: 427–433
35. Hutchinson GE (1965) The Ecological Theater and the Evolutionary Play. New Haven: Yale University Press. 139 p.
36. Soberón J, Peterson AT (2005) Interpretation of models of fundamental ecological niches and species' distributional areas. Biodiv Inform 2: 1–10
37. Soberón J (2007) Grinnellian and Eltonian niches and geographic distributions of species. Ecol Lett 10: 1115–1123 [PubMed]
38. Godsoe W (2010) I can't define the niche but I know it when I see it: a formal link between statistical theory and the ecological niche. Oikos 119: 53–60
39. Puliam HR (2000) On the relationship between niche and distribution. Ecol Lett 3: 349–361
40. Jackson ST, Overpeck JT (2000) Responses of plant populations and communities to environmental changes of the late Quaternary. Palaeobiology 26: 194–220
41. Veloz SD, Williams JW, Blois JL, He F, Otto-Bliesner B, et al. (2012) No-analogue climates and shifting realized niches during the late quaternary: implications for 21st-century predictions by species distribution models. Glob Change Biol doi:10.1111/j.1365-2486.2011.02635.x
42. Rödder D, Lötters S (2010) Explanative power of variables used in species distribution modelling: an issue of general model transferability or niche shift in the invasive Greenhouse frog (Eleutherodactylus planirostris). Naturwissenschaften 97: 781–796 [PubMed]
43. Waltari E, Hijmans RJ, Peterson AT, Nyári AS, Perkins SL, et al. (2007) Locating pleistocene refugia: comparing phylogeographic and ecological niche model predictions. PLoS ONE 2: e563. [PMC free article] [PubMed]
44. Pearman PB, Randin CF, Broennimann O, Vittoz P, van der Knaap WO, et al. (2008) Prediction of plant species distributions across six millennia. Ecol Lett 11: 357–369 [PubMed]
45. Yesson C, Culham A (2006) Phyloclimatic modeling: combining phylogenetics and bioclimatic modeling. Sys Biol 55: 785–802 [PubMed]
46. Stigall AL, Lieberman BS (2006) Quantitative paleobiogeography: GIS, phylogenetic biogeographical analysis, and conservation insights. J Biogeogr 33: 2051–2060
47. Maguire KC, Stigall AL (2009) Using ecological niche modeling for quantitative biogeographic analysis: a case study of Miocene and Pliocene Equinae in the Great Plains. Paleobiology 35: 587–611
48. Poelchau MF, Hamrick JL (2011) Palaeodistribution modelling does not support disjunct Pleistocene refugia in several Central American plant taxa. J Biogeogr doi:10.1111/j.1365-2699.2011.02648.x
49. Polly PD, Eronen JT (2011) Mammal associations in the Pleistocene of Britain: implications of ecological niche modelling and a method for reconstrucing palaeoclimate. In: Ashton N, Lewis SG, Stringer C, editors. The Ancient Human Occupation of Britain. Amsterdam: Elsevier. pp. 279–304.
50. Stewart JR (2009) The evolutionary consequence of the individualistic response to climate change. J Evol Biol 22: 2363–2375 [PubMed]
51. Elith J, Kearney M, Phillips S (2010) The art of modelling range-shifting species. Method Ecol Evol 1: 330–342
52. Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, et al. (2011) A statistical explanation of MaxEnt for ecologists. Divers Distrib 17: 43–57
53. McMahon SM, Harrison SP, Armbruster WS, Bartlein PJ, Beale CM, et al. (2011) Improving assessment and modelling of climate change impacts on global terrestrial biodiversity. Trends Ecol Evol 26: 249–259 [PubMed]
54. Graham CH, Ron SR, Santos JC, Schneider CJ, Moritz C (2004) Integrating phylogenetics and environmental niche models to explore speciation mechanisms in dendrobatid frogs. Evolution 58: 1781–1793 [PubMed]
55. Hardy CR, Linder HP (2005) Intraspecific variability and timing in ancestral ecology reconstruction: A test case from the Cape flora. Sys Biol 54: 299–316 [PubMed]
56. Evans MEK, Smith SA, Flynn RS, Donoghue MJ (2009) Climate, niche evolution, and diversification of the “bird-cage” evening primroses (Oenothera, sections Anogra and Kleinia). Am Nat 173: 225–240 [PubMed]
57. Stephens PR, Wiens JJ (2009) Bridging the gap between community ecology and historical biogeography: niche conservatism and community structure in emydid turtles. Mol Ecol 18: 4664–4679 [PubMed]
58. Iverson JB (1992) Species richness maps of the freshwater and terrestrial turtles of the world. Smithsonian Herpetological Information Service 88: 1–18
59. Ihlow F, Dambach J, Engler JO, Flecks M, Hartmann T, et al. (2012) On the brink of extinction? How climate change may affect global chelonian species richness and distribution. Glob Change Biol 18: 1520–1530
60. World Turtle Database (2012) Available: Accessed 2011 Aug 26.
61. Holman JA (1995) Pleistocene Amphibians and Reptiles in North America. New York: Oxford University Press. 243 p.
62. van Dijk PP, Iverson JB, Shaffer BH, Bour R, Rhodin AGJ (2011). Conservation Biology of Freshwater Turtles and Tortoises: A Compilation Project of the IUCN/SSC Tortoise and Freshwater Turtle Specialist Group. In: Rhodin AGJ, Pritchard PCH, van Dijk PP, Saumure RA, Buhlmann KA, Iverson JB, Mittermeier RA, editors. Chelonian Research Monographs 5. pp. 165–242.
63. Makarieva AM, Gorshkov VG, Li B-L (2005) Gigantism, temperature, and metabolic rate in terrestrial poikilotherms. Proc Royal Soc B-Biol Sci 272: 2325–2328 [PMC free article] [PubMed]
64. Valenzuela N (2004) Temperature-dependent sex determination. In: Deeming DE, editor. Reptilian incubation. Nottingham: Nottingham Press. pp. 1–16.
65. Fritz U, Havaš P (2007) Checklist of Chelonians of the World. Vertebr Zool 57: 149–368
66. Buhlmann KA, Akre TSB, Iverson JB, Karapatakis D, Mittermeier RA, et al. (2009) A Global Analysis of Tortoise and Freshwater Turtle Distributions with Identification of Priority Conservation Areas. Chelonian Conserv Biol 8: 116–149
67. Rhodin AGJ, van Dijk PP, Iverson J, Shaffer HB (2010) Turtles of the World, 2010 update: Annotated checklist of taxonomy, synonymy, distribution, and conservation status. Chelon Res Monogr 5doi:10.3854/crm.5.000.checklist.v3.2010
68. Hernandez PA, Graham CH, Master LL, Albert DL (2006) The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography 29: 773–785
69. Wisz MS, Hijmans RJ, Li J, Peterson AT, Graham CH, et al. (2008) Effects of sample size on the performance of species distribution models. Divers Distrib 14: 763–773
70. Hibbard CW, Taylor DW (1960) Two late Pleistocene Faunas from Southwestern Kansas. UM Contr Mus Paleontol 16: 1–223
71. Preston RE (1979) Late Pleistocene cold-blooded vertebrate faunas from the mid-continental United States. I. Reptilia; Testudines, Crocodilia. Ann Arbor: The Museum of Paleontology, University of Michigan. 53 p.
72. Holman JA (1971) Herpetofauna of the Sandahl Local Fauna (Pleistocene: Illinoian) of Kansas. UM Contr Mus Paleontol 23: 349–335
73. Holman JA (1987) Climatic significance of a late Illinoian Herpetofauna from Southwestern Kansas. UM Contr Mus Paleontol 27: 129–141
74. Harris AH (1993) Quaternary Vertebrates of New Mexico. New Mexico Mus Nat Hist Sci Bull 2: 179–197
75. Woodburne MO (2004) Late Cretaceous and Cenozoic Mammals of North America. New York: Columbia University Press. 391 p.
76. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 25: 1965–1978
77. Nix HA (1986) A biogeographic analysis of Australian elapid snakes. In: Longmore R, editor. Atlas of elapid snakes of Australia. Australian Flora and Fauna, Ser.7. Canberra: Australian Government Publishing Service. pp. 4–15.
78. Heikkinen RK, Luoto M, Araújo MB, Virkkala R, Thuiller W, et al. (2006) Methods and uncertainties in bioclimatic envelope modeling under climate change. Progr Phys Geogr 30: 751–777
79. Lisiecki LE, Raymo ME (2005) A Pliocene-Pleistocene stack of 57 globally distributed benthic d18O records. Paleoceanogr 20: PA10023
80. Braconnot P, Otto-Bliesner B, Harrison S, Joussaume S, Peterschmitt J-Y, et al. (2007) Results of PMIP2 coupled simulations of the mid-Holocene and Last Glacial Maximum. Part 1: experiments and large-scale features. Clim Past 3: 261–277
81. Phillips SJ, Anderson RP, Schapire RE (2006) Maximum entropy modeling of species geographic distributions. Ecol Model 190: 231–259
82. Phillips SJ, Dudík M (2008) Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography 31: 161–175
83. Nogués-Bravo D (2009) Predicting the past distribution of species climatic niches. Glob Ecol Biogeogr 18: 521–531
84. Liu C, Pam M, Dawson TP, Pearson RG (2005) Selecting thresholds of occurrence in the prediction of species distributions. Ecography 28: 385–393
85. Venables WN, Ripley BD (2002) Modern Applied Statistics with S. Fourth Edition. New York: Springer. 512 p.
86. Thomson RC, Shaffer HB (2010) Rapid progress on the vertebrate tree of life. BMC Biology 8: 19. [PMC free article] [PubMed]
87. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797 [PMC free article] [PubMed]
88. Kück P, Meusemann K, Dambach J, Thormann B, von Reumont BM, et al. (2010) Parametric and non-parametric masking of randomness in sequence alignments can be improved and leads to better resolved trees. Front Zool 7: 10. [PMC free article] [PubMed]
89. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A (2006) Relaxed Phylogenetics and Dating with Confidence. PLoS Biol 4: e88.doi:10.1371/journal.pbio.0040088 [PMC free article] [PubMed]
90. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7: 214. [PMC free article] [PubMed]
91. Fritz U, Stuckas H, Vargas-Ramírez M, Hundsdörfer AK, Maran J, et al. (2012) Molecular phylogeny of Central and South American slider turtles: implications for biogeography and systematics (Testudines: Emydidae: Trachemys). J Zool Syst Evol Res 50: 125–136
92. Holman JA, Fritz U (2001) A new emydine species from the Middle Miocene (Barstovian) of Nebraska, USA with a new generic arrangement for the species of Clemmys sensu McDowell (1964) (Reptilia: Testudines: Emydidae). Zool Abh 51: 331–354
93. Bourque JR (2012) A fossil mud turtle (Testudines, Kinosternidae) from the early middle Miocene (early Barstovian) of New Mexico. J Vertebr Paleontol 32: 836–853
94. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690 [PubMed]
95. Walker D, Avise JC (1998) Principles of phylogeography as illustrated by Freshwater and Terrestrial turtles in the Southeastern United States. Ann Rev Ecol Syst 29: 23–58
96. Phillips SJ, Dudík M, Elith J, Graham CH, Lehmann A, et al. (2009) Sample selection bias and presence-only distribution models: implications for background and pseudo-absence data. Ecol Appl 19: 181–97 [PubMed]
97. Rocchini D, Hortal J, Lengyel S, Lobo JM, Jimenez-Valverde A, et al. (2011) Accounting for uncertainty when mapping species distributions: The need for maps of ignorance. Prog Phys Geogr 35: 211–226
98. Bowman AW, Azzalini A (2010) R package “sm”: nonparametric smoothing methods (version 2.2-4). Available:,
99. Peterson AT, Soberón J, Sánchez-Cordero V (1999) Conservatism of ecological niches in evolutionary time. Science 285: 1265–1267 [PubMed]
100. Warren DL, Glor RE, Turelli M (2008) Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution 62: 2868–2883 [PubMed]
101. Broennimann O, Fitzpatrick MC, Pearman PB, Petitpierre B, Pellissier L, et al. (2012) Measuring ecological niche overlap from occurrence and spatial environmental data. Global Ecology and Biogeography 21: 481–497
102. Carstens BC, Richards CL (2007) Integrating coalescent and ecological niche modeling in comparative phylogeography. Evolution 61: 1439–1454 [PubMed]
103. Weisrock DW, Janzen FJ (2000) Comparative molecular phylogeography of North American Softshell Turtles (Apalone): Implications for regional and wide-scale historical evolutionary forces. Mol Phylogen Evol 13: 152–164 [PubMed]
104. McGaugh SE, Eckermann CM, Janzen FJ (2008) Molecular phylogeography of Apalone spinifera (Reptilia, Trionychidae). Zool Scr 37: 289–304
105. Walker DE, Moler PE, Buhlmann KA, Avise JC (1998) Phylogeographic uniformity in mitochondrial DNA of the snapping turtle (Chelydra serpentina). Anim Conserv 1: 55–60
106. Spinks PQ, Shaffer B (2005) Range-wide molecular analysis of the western pond turtle (Emys marmorata): cryptic variation, isolation by distance, and their conservation implications. Mol Ecol 14: 2047–2064 [PubMed]
107. Starkey DE, Shaffer HB, Burke RL, Forstner MRJ, Iverson JB, et al. (2003) Molecular systematics, phylogeography and the effects of Pleistocene glaciation in the Painted turtle (Chrysemys picta) complex. Evolution 57: 119–128 [PubMed]
108. Amato ML, Brooks RJ, Fu J (2008) A phylogeographic analysis of populations of the wood turtle (Glyptemys insculpta) throughout its range. Mol Ecol 17: 570–581 [PubMed]
109. Rosenbaum PA, Robertson JM, Zamudio KR (2007) Unexpectedly low genetic divergences among populations of the threatened bog turtle (Glyptemys muhlenbergii). Conserv Genet 8: 331–342
110. Ennen JR, Lovich JE, Kreiser BR, Selman W, Qualls CP (2010) Genetic and morphological variation between populations of the Pascagoula Map Turtle (Graptemys gibbonsi) in the Pearl and Pascagoula rivers with description of a new species. Chelonian Conserv Biol 9: 98–113
111. Hutchinson VH, Vinegar A, Kosh RJ (1966) Critical thermal maxima in turtles. Herpetologica 22: 32–41
112. Pritchard PCH, Greenhood W (1968) The sun and the turtle. Int Turtle Tort Soc J 2: 20–25
113. Mahmoud IY (1969) Comparative ecology of the kinosternid turtles of Oklahoma. Southwest Natur 14: 31–66
114. Bury RB (1972) Habits and home range of the pacific pond turtle Clemmys marmorata, in a stream community. PhD thesis. Berkeley: University of California.
115. Iverson JB, Brown RM, Akre TS, Near TJ, Le M, et al. (2007) In Search of the Tree of Life for Turtles. Chelonian Res Monogr 4: 85–106
116. Stewart JR (2008) The progressive effect of the individualist response of species to Quaternary climate change: an analysis of British mammalian faunas. Quaternary Science Reviews 27: 2499–2508
117. Mockford SW, Herman TB, Snyder M, Wright JM (2007) Conservation genetics of Blanding's turtle and its application in the identification of evolutionary significant units. Conserv Genet 8: 209–219
118. Tessier N, Paquette S, Lapointe FJ (2005) Conservation genetics of the wood turtle (Glyptemys insculpta) in Quebec, Canada. Canadian J Zool 83: 765–772
119. Dobrowski SZ (2011) A climatic basis for microrefugia: the influence of terrain on climate. Glob Change Biol 17: 1022–1035
120. Rull V (2009) Microrefugia. J Biogeogr 36: 481–484
121. Wright S (1943) Isolation by Distance. Genetics 28: 114–138 [PubMed]
122. Melbourne BA, Hastings A (2008) Extinction risk depends strongly on factors contributing to stochasticity. Nature 454: 100–103 [PubMed]
123. Serb JM, Phillips CA, Iverson JB (2001) Molecular phylogeny and biogeography of Kinosternon flavescens based on complete mitochondrial control region sequences. Mol Phyl Evol 18: 149–162 [PubMed]
124. Bailey LA, Dixon JR, Hudson R, Forstner MRJ (2008) Minimal genetic structure in the Rio Grande Cooter (Pseudemys gorzugi). Southwest Nat 53: 406–411
125. Hauswaldt JS, Glenn TC (2005) Population genetics of the diamondback terrapin (Malaclemys terrapin). Mol Ecol 14: 724–732 [PubMed]
126. Echelle AA, Hackler JC, Lack JB, Ballard SR, Roman J, et al. (2010) Conservation genetics of the alligator snapping turtle: cytonuclear evidence of range-wide bottleneck effects and unusually pronounced geographic structure. Conserv Genet 11: 1375–1387
127. Reinig WF (1938) Elimination und Selektion. Jena: Gustav Fischer. 154 p.
128. Dalén L, Nyström V, Valdiosera C, Germonpré M, Sablin M, et al. (2007) Ancient DNA reveals lack of postglacial habitat tracking in the arctic fox. PNAS 104: 6726–6729 [PubMed]
129. Ernst CH, Lovich JE (2009) Turtles of the United States and Canada. Second Edition. Baltimore: The Johns Hopkins University Press. 827 p.
130. Crews D, Bergeron JM, Bull JJ, Flores D, Tousignant A, et al. (1994) Temperature-dependent sex determination in reptiles: proximate mechanisms, ultimate outcomes, and practical applications. Dev Genet 15: 297–312 [PubMed]
131. Escobedo-Galván AH, González-Salazar C (2012) Survival and extinction of sex-determining mechanisms in Cretaceous tetrapods. Cretaceous Res 36: 116–118
132. Hulin V, Delmas V, Girondot M, Godfrey MH, Guillon J-M (2009) Temperature-dependent sex determination and global change: are some species at greater risk? Oecologia 160: 493–506 [PubMed]
133. Mitchell NJ, Janzen FJ (2010) Temperature-dependent sex determination and contemporary climate change. Sex Dev 4: 129–140 [PubMed]
134. Walther G-R, Post E, Convey P, Menzel A, Parmesan C, et al. (2002) Ecological responses to recent climate change. Nature 416: 389–395 [PubMed]
135. Silber S, Geisler JH, Bolortsetseg M (2011) Unexpected resilience of species with temperature-dependent sex determination at the Cretaceous-Palaeogene boundary. Biol Lett 7: 295–298 [PMC free article] [PubMed]
136. Valenzuela N, Adams DC (2011) Chromosome number and sex determination coevolve in turtles. Evolution 65: 1808–1813 [PubMed]
137. Mortiz C (1994) Defining ‘Evolutionary Significant Units’ for conservation. Trends Ecol Evol 9: 373–375 [PubMed]
138. Jackson ST, Williams JW (2004) Modern analogs in Quaternary paleoecology: Here today, gone yesterday, gone tomorrow? Annu Rev Earth Planet Sci 32: 495–537
139. Williams JW, Jackson ST (2007) Novel climates, no-analog plant communities, and ecological surprises: past and future. Front Ecol Evol 5: 475–482
140. Bálint M, Domisch S, Engelhardt CHM, Haase P, Lehrian S, et al. (2011) Cryptic biodiversity loss linked to global climate change. Nat Clim Change 1: 1–6doi: 10.1038/nclimate1191
141. Phillips CA, Dimmick WW, Carr JL (1994) Conservation genetics of the common Snapping Turtle (Chelydra serpentina). Conserv Biol 10: 397–405
142. Schwarz A (1956) Geographic variation in the Chicken turtle Deirochelys reticularia Latreille. Fieldiana Zool 34: 461–503
143. Mockford SW, McEachern L, Herman TB, Snyder M, Wright JM (2004) Population genetic structure of a disjunct population of Blanding's turtle (Emydoidea blandingii) in Nova Scotia, Canada. Biol Conserv 123: 373–380
144. Folkerts GW, Mount RH (1969) A new subspecies of the turtle Graptemys nigrinoda Cagle. Copeia 1969: 677–682
145. Smith AD (2008) Intraspecific phylogeography of Graptemys ouachitensis MsC thesis, Ohio University.
146. Walker DE, Moler PE, Buhlmann KA, Avise JC (1998) Phylogeographic patterns in Kinosternon subrubrum and K. baurii based on mitochondrial DNA restriction analyses. Herpetologica 54: 174–184
147. Roman J, Santhuff SD, Moler PE, Bowen BW (1999) Population structure and cryptic evolutionary units in the Alligator Snapping Turtle. Conserv Biol 13: 135–142
148. Lamb T, Avise JC (1992) Molecular and population genetic aspects of mitochondrial DNA variability in the diamondback terrapin, Malaclemys terrapin. J Heredity 83: 262–269
149. Walker D, Burke VJ, Barák I, Avise JC (1995) A comparison of mtDNA restriction sites vs. control region sequences in phylogeographic assessment of the musk turtle (Sternotherus minor). Mol Ecol 4: 365–673 [PubMed]
150. Walker DE, Nelson WS, Buhlmann KA, Avise JC (1997) Mitochondrial DNA phylogeography and subspecies issues in the monotypic freshwater turtle Sternotherus odoratus. Copeia 1: 16–21
151. Dodd K (2002) North American Box Turtles: a natural history. Norman: University of Oklahoma Press. 231p.
152. Butler JM, Dodd CK, Aresco M, Austin JD (2011) Morphological and molecular evidence indicates that the Gulf Coast box turtle (Terrapene carolina major) is not a distinct evolutionary lineage in the Florida Panhandle. Biol J Linn Soc Lond 102: 889–901
153. Howeth JG, McGaugh SE, Hendrickson DA (2008) Contrasting demographic and genetic estimates of dispersal in the endangered Coahuilan box turtle: a contemporary approach to conservation. Mol Ecol 17: 4209–4221 [PubMed]
154. Kuo C-H, Janzen FJ (2004) Genetic effects of a persistent bottleneck on a natural population of ornate box turtles. Conserv Genet 5: 425–437
155. Seidel ME (2002) Taxonomic observations on extant species and subspecies of Slider Turtles, genus Trachemys. J Herpetol 36: 285–292
156. McCord WP, Joseph-Ouni M, Hagen C, Blanck T (2010) Three new subspecies of Trachemys venusta (Testudines: Emydidae) from Honduras, northern Yucatán (Mexico), and pacific coastal Panama. Reptilia 71: 39-39–49
157. Araújo MB, Peterson AT (2012) Uses and misuses of bioclimatic envelope modeling. Ecology 93: 1527–1539 [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science