1. Introduction
Biogeography is closely tied to both ecology and phylogenetic biology and its main areas of interest are ecological biogeography, i.e. the study of factors influencing the present distribution, and historical biogeography, i.e. the study of causes that have operated in the past [1]. Ecological and historical biogeography therefore applies different concepts in order to explain the distribution of organisms. The former deals with functional groups of species and environmental constraints, whereas the latter focuses on taxonomic groups and historical biogeographical events [2]. Recently, the division between historical and ecological biogeography has been considered as an obstacle to the progress of biogeography and some authors have stressed the benefits of integrating these two points of view [1,3-4]. In this context, the present work attempts at integrating these two approaches to explore the regional patterns of endemism and species richness of the squamates(lizards and snakes) in China.
As summarized by Meng
It has long been recognized that geological complexity and history usually have a profound influence on the distributions of living organisms. China animal geography division was firstly put forward in 1959 [8]. Now based on the distributions of vertebrates, mainly mammals and birds, China was divided into seven major biogeographical regions and fifty-four provinces [9]. The Palaearctic realm includes North-eastern, Northern, Inner Mongolia-Xinjiang, and Qinghai-Tibetan China biogeographical regions, while the Oriental realm was divided into South-western, Central, and Southern biogeographical regions. The study on biogeographical divisions of China was both intriguing and challenging due to the complex topography and historical processes [9]. The zoogeographic division of China was explored into 30 units [10]. However, China’s territory was divided into 124 basic units on the basis of comprehensive natural factors (including altitude, landform, climate, vegetation, water system, farm belt, and so on) by cluster analysis [11]. Recently, based on the distributional patterns of spiders, seven major biogeographical regions in China were investigated by PAE, not corresponding to any previous studies [5]. These studies played an important role in biodiversity conservation, planning and management in China. However, these divisions need further investigation due to primarily on basis of mammals and birds which are highly adapted to environment diversity and higher locomotion [12]. Thus, the ectothermic animals may act as a better indicator to determinate zoogeographical division than mammals and birds [12]. Although only a few studies on biogeographical patterns have been identified and emphasized their biogeographical complexity, data are still not enough to evaluate and compare directly for other taxa, such as squamates with limited dispersal abilities. It is necessary to compare patterns based on different groups of organisms to better understand their biogeography and infer a general pattern [13].
2. Compared biogeographical patterns of Chinese squamates based on parsimony analysis of endemicity (PAE) at different natural area units
2.1. Introduction
One of important goals of historical biogeography is to investigate convergent biogeographical patterns relying on different taxa [14]. These may assist in identifying priority areas or hotspots for biodiversity conservation, particularly today the issue of global concern biodiversity loss. A historical biogeographical method, parsimony analysis of endemicity(PAE), firstly proposed by Rosen [15] and further elaborated by Morrone [16], provided an insight to generate area cladograms to make inferences on historical patterns. PAE was originally used the most parsimony algorithm to reconstruct relations among sampling localities [15], then previously delimited areas [17] and quadrats [16]. Analogous to cladistic methods in phylogenetic analyses, PAE classifies areas (cf. taxa in cladistics) on the basis of the shared presence of taxa (cf. characters in cladistics) [15]. Using PAE, biogeographical studies can investigate biotic similarities between different geographical regions, give static or ecological interpretations [18], and estimate historical hierarchical congruence in target localities or geographical regions [19]. Although there is ongoing debate about the value of PAE [20-21], the PAE has been widely used in many biogeographical studies in recent years [5,22-23].
PAE has been applied to establish relationships among different areas units, for example localities, quadrats, areas of endemism, continents, islands, and so on [13]. The ideal organisms for using PAE are those with limited dispersal abilities and speciation in vicariant events [13,19]. Although the same size and shape quadrats are not required in PAE and do not affect the analyses [24], the best PAE results were obtained with natural areas instead of quadrats [13]. Here, we used PAE to compare squamates biogeographical patterns at different natural area units with previously delimited biogeographic patterns obtained by Meng and Murphy [5] (Figure1), Zhang [9](Figure 2) and Xie
2.2. Materials and methods
2.2.1. Study area and operational geographical units
For a comparable study area with previously provided optimal results in detecting the biogeographical patterns, we assumed the same operational geographic units (OGUs) as suggested by Meng and Murphy [5] (Figure 1), Zhang [9] (Figure 2) and Xie et al. [11] (Figure 3) defined, respectively. According to very similar climatic, geological characteristics, topographical characteristics and natural barriers to dispersal for neighbouring quadrats, 28 biological and physiographical similar areas were combined and divided as OGUs in PAE [5]. Zoogeographical divisions of China were classified into 54 biogeographical provinces based on the distribution of mainly mammals and birds [9]. The whole China was divided into 124 basic units using comprehensive natural factors, such as altitude, landform, climate, vegetation, water system, farm belt, and so on [11]. More details about the OGUS see these studies.
2.2.2. Data sources and dataset
The Squamates species catalogue and distributions were compiled from the most recent and comprehensive references, including specimens, exhaustive field surveys, monograph, published literature and expert reviews. All date Information from herbarium specimens was mainly obtained through Chengdu Institute of Biology, Chinese Academy of Sciences (CIB/CAS). Additional information was also obtained from HerpNET and a variety of published sources, for example
2.2.3. Parsimony Analysis of Endemicity (PAE)
We used PAE to identify biogeographical patterns and followed the procedure modified by Morrone [16]. A taxon/area matrix for the basic PAE data set was built in which the absence of a species in an area was coded as ‘0’ and presence as ‘1’. A hypothetical area coded ‘0’ as an outgroup to provide a root for the final cladogram [13, 16]. We removed species that occurred in all areas, as well as those in only one area due to phylogenetically uninformative autapomorphies [15]. The taxon/area matrix were imported into PAUP * [29] to find the most parsimonious cladograms with a heuristic search of 1000 replicates and random sequence additions. We estimated relative support for each branch using bootstrapping, with 100 replicates and tree bisection–reconnection (TBR) swapping. All characters were weighted equally and a 50% majority consensus tree of the equally parsimonious trees was generated.
2.3. Results
2.3.1. Biological areas 28
The bootstrap 50% majority-rule consensus tree using 244”characters” in the analysis was shown in Figure 1(b). Our results were substantially different from the results of Meng
For Group E, it was very complex and ten clades were discovered. Subclade E1 as Sothern region (bootstrap 88%) corresponded to the clade C in [11], but major regions such as Central region (C2), Eastern Southern region (C3), Western Southern region (C4), and Central Southern region (C5) in [5] cannot be recognized. In subclade E1, A8 mainly representing Jiangsu was basal, followed by separation of A1 and the remaining areas in subclade E1. Areas A22, A24, A25, A26 and A27 formed a clade with a moderate support of 78%, corresponding to Southern China of zoogeographical division of [9] excluding Taiwan Island. Subclade E2 had a bootstrap support of 98%, including two sub-clades, A12 and A5 + A28. It corresponded to northeastern Tsaidam basin, Loess plateau and Alashan Plateau. Subclade E3 (A10 + A19) corresponded to the steppe and desert of north-western China, including Xinjiang and Inner Mongolia. Its bootstrap support value was 72%. Subclade E4 (A14 + A18) consists of the Xiao Xingan and Changbai mountains, with a bootstrap value of 86%. The Qinghai-Tibetan Plateau included A11+A21, however, subclade E6 (A11) was separated from A21 in the results with a bootstrap support value of 64%. The remaining subclades showed ambiguous.
2.3.2. Biological areas 54
The 1507 most parsimonious cladograms of 930 steps with CI of 0.3172, HI of 0.6828, and RI of 0.6340 were found. The 50% majority-rule consensus tree using 295”characters” (Figure 2, b) showed a basal polytomy, but several clades emerged. Clade F consisted of areas A12, A13 and A14 with a bootstrap value of 57%, which corresponded to the Loess plateau subregion. Clade G had a weakly supported group containing areas A18, A19, A20, A21, A23 and A24 with a bootstrap proportion of 57%, which corresponded to Western desert subregion and Tian Shan mountains subregion of Mognolia-Xijiang China. Clade H had a well-support of 95%, including two sub-clades, A30 + A36 and A37, in response with the Himalaya mountains subregion of Southwestern China. Areas A40, A44, A45, A46, and A47 formed a clade I, including Eastern hills and plains subregion and Western mountains and plateau subregion of Central China, and coastal subregions of Guangdong and Fujian provinces of Southern China; its bootstrap value was 57%. Clade J consisted of areas A48 and A49 with a strongly support value of 100%, corresponding to Southern Yunnan mountains subregion of Southern China. Clade K corresponded to Hainan Island subregion of Southern China, consisting of areas A50 and A51 with 75% bootstrap value. Clade L contained Taiwan Island areas with similar Squamates faunas, areas A52 and A53; its well-supported value was 100%.
2.3.3. Biological areas 124
For OGUs of biological areas 124, the 1000 most parsimonious trees were obtained. The 50%-majority consensus tree using 294 characters (tree length = 1407, CI = 0.2090, RI = 0.7910) of the 1000 trees was shown in Figure 2 (b). Like the Figure 3, a basal polytomy were identified where only one of the branches contained a dicotomy in Figure 3. Clade M is moderately supported by Areas A10 and A14 with 88% bootstrap value. It corresponded to Yinshan mountains hill and Ordos Plateau. Areas (A21, ((A22, A23), A17)) formed the clade N with a week support value of 62%, corresponding to (Mountain front mesa, ((Liaodong Peninsula, Lower Liaohe Plain), Changbai mountains)) in Eastern Northeastern China. Clade O consisted of areas A29 and A32, including Southeast Shanxi plain and Western Henan mountains. Its bootstrap support proportion was 60%. Clade P contained areas (A101, (A37, A99)) with 68% bootstrap value, corresponding to the range of Qilian Mountain, Central Gansu incisive hill and Upper Yellow River incisive mountains. Areas A46 (Qinling mountains) and A47 (Daba–Micang mountains) in Central China formed clade Q with a moderate support value of 81%. Clade R was compose of Areas A58 and A63, corresponding to Xiangjiang valley hill and Honghe catchment montane basin; its support proportion was 80%. Clade S included Areas A61, A62 and A63 with a weekly support value of 51%, corresponding to Dalou Mountain mid-land valley, Miaoling hilly plain and Wujiang and Nanpanjiang catchments mid-land valley. Clade T was the largest clade including Areas A67, A68, A73, A74, A75 and A87 with 60% support value. It corresponded to Southeast Yunnan low-heat plateau, South central Yunnan low-heat valley, West Yunnan montane plain, Southwest Yunnan Plateau wide valley, South Yunnan wide valley and Salween and Lancang Rivers parallel valley. Areas A76 (West Guangdong and south Guangxi coastal mesa plain) and A79 (South Hainan montane hill) formed Clade U with 64% support value. Clade V had the greatest bootstrap value (98%), including areas A80, A81 and A83. It corresponded to Northwest subtropical hilly plain, Central subtropical mountain and East tropical coast in Taiwan Island. Clade W consisted of areas A92, A93, A94 and A95 with a bootstrap support of 73%, corresponding to areas in Himalayas mountains including Kangrigebu south wing mountains, Himalayas south wing mountains, Salween and Lancang Rivers incisive mountains and Brahmaputra Great Turn and upper Salween incisive mountains. Similarly, clade X consisted of areas A96, and A97, corresponding to areas in Himalayas mountains including Brahmaputra valley mountains and Himalayas central mountains; its bootstrap support value was 61%. Clade Y was compose of areas A104 and A105 with a week bootstrap support proportion of 59%. It corresponded to areas in Qinghai–Tibetan Plateau including South Qiangtang Plateau mountains and North Tibet plateau northwestern lake basin mountains. Clade Z included areas (A111, (A115, (A120, A122))), which corresponded to West Hexi Corridor, East Tianshan mountains, Central Tianshan mountains, Tarim Basin. Its bootstrap value was 57%. Areas A118 (Junggar Basin), A121 (Ili Valley), and A119 (Emin Valley) formed Clade AZ with a well support proportion of 97%.
2.4. Conclusion
Past major geological events have played important roles in shaping the biogeographic distribution of extant organisms. PAE originally aimed to find areas of congruent distributional patterns, and the best PAE results were obtained with natural areas (e.g. biogeographical provinces, ecoregions) instead of quadrats by increasing the absolute and relative numbers of synapomorphies [13]. If we compare the results of areas 28, 54 and 124, there seems to exist a trend to result in poor resolution of the resultant area cladograms as the size of OGU decreases. Although the samples within the regions are not uniform at different areas unit, in our area cladograms, biogeographic patterns of squamates distributions appear to have a hierarchical structure and general patterns for areas 28. Based on the above comprehensive squamates distributions patterns at different natural area units, seven major congruent biogeographic regions can be identified in China: Eastern Northern region, Tibetan Plateau region, Xinjiang and Inner Mongolia region, Loess plateau and Alashan Plateau region, Taiwan Island region and Southern region. However, there existed several unresolved areas relationships to one another, the uncertain position of certain areas or incongruence in the ‘character’ distributions. There may be basal problems with data themselves such as squamates distribution incomplete known. Sure, extinction, long-distance dispersal, isolation, or other undetected historical patterns may lead to incongruence in the distributions patterns [30].
It is surely not coincidental that different organisms may share a general distribution pattern. A PAE area cladogram might contain areas related by shared ecologies or similar historical events (biotic divergence and isolation) [15]. Historical hypothesis, in other words, evolutionary history has recently been considered to be a driving force determining squamates regional species pools’ differences in China. Geological complexity and history usually have a profound influence on the distributions of living organisms in China [5], though most of China has never been covered by ice sheets. Squamates species are followed this rule. For example, squamates were similar to spiders distributional patterns broadly corresponded to geological provinces in China, such as Southern China geological province versus Southern regions, The Laurentian/Cathaysian Southern, South-western Margin geological province and Tibetan geological provinc vs Tibetan Plateau biogeographical region [9]. Furthermore, biotic and abiotic conditions are also important factors in determinative of the distribution of squamates species. Some correlations between species richness and reproductive modes with geography and ecological conditions have been reported [31-32]. The interpretations of those relationships have postulated that contemporary factors are the main regulatory force of the distribution of squamates taxa in China [33-35]. Thus, both history and ecology may well be inseparable and have a profound impact on not only the diversity of Squamates taxa but also their biogeographic patterns.
3. Distribution patterns in species diversity of lizards in China and their relationships to ecological factors
3.1. Introduction
Determining the causes of the great biodiversity variation across Earth has long been a major challenge for ecologists and biogeographers [36], ever since biotic diversity contrast between equatorial and polar latitudes was discovered two centuries ago [37]. Among the considerable number of hypotheses that aim to explain species richness patterns [36, 38], many ecological (environmental) hypotheses have been widely discussed and accepted [39].
Three alternative variants of ecological hypothesis, the species-energy, contemporary climate and habitat heterogeneity hypotheses, have received a great deal of attention as the primary determinants of species richness [39-43]. The species-energy hypothesis includes at least two versions, the ambient energy and productive energy hypotheses [38]. The ambient energy hypothesis, widely indicated by temperature or allied measures, argues that species richness was influenced by energy inputs into an area that affects the physiological tolerance of organisms [40, 44]. The productive energy hypothesis claims that animal species richness is limited by energy via food webs rather than by physiological requirements. The energy and water availability (i.e., energy–water dynamics) limits the total available plant productivity, which ultimately moves up the food chains [40, 45-47]. The contemporary climate hypothesis states that species richness correlates with contemporary climate conditions, and putative causal mechanisms are in terms of environmental stability, variability, favorability and harshness [40, 47-48]. The habitat heterogeneity hypothesis is measured either as the number of habitat types or the topographic relief (range in elevation) presented within an area [42, 50]. It assumes that high species richness is found in physically or biologically complex habitats, through higher speciation rates and providing more ecological niches [40-42].
However, the knowledge of the determinants of reptile richness remains insufficiently documented among terrestrial vertebrates [39, 51-82]. It is urgent to understand the drivers of reptile richness patterns due to global warming impact on species distribution and abundance [36, 52-53]. Lizards belong to Reptila and are good model systems to test these alternative hypotheses. Because their taxonomy is well resolved and distributional data are quite thorough. They are ectothermic and sensitive to environmental variables. In this study, we examine the correlation between lizard species richness and various environmental factors across China. Our objectives include (1) mapping distributions of Chinese lizards and describing any patterns, and (2) testing various ecological factors in determining species richness patterns.
3.2. Materials and methods
3.2.1. Data collection
We collected locality data for lizard species which occur in China from a variety of sources as above-mentioned. We excluded coastal grid cells with less than 96% land cover and all islands from the analysis in order to remove the effects of insularity. Finally, we built a database of 151 lizard with species names (see Appendix S1 in Supplementary Material) represented by a total of 3,391 records for unique point localities, with a range of 2–288 (mean = 22.5, standard deviation = 38.3).
3.2.2. Ecological niche modeling and species richness
For each of the 151 species, we used the Genetic Algorithm for Rule-set Prediction (GARP) (for free download see: http://www.nhm.ku.edu/desktopgarp/) [54] for reconstructing species distribution maps. GARP uses an evolutionary computing genetic algorithm to search iteratively for non-random correlations between species presence and environmental variables for localities using several different types of rules (i.e., atomic rules, range rules, negated range rules, logistic regression rules), and then creates ecological niche models for each species’ predicted distribution, as contrasted with environmental characteristics across the overall study area [54]. GARP was found that it did not tend to be more sensitive to sampling bias than Maxent, and GARP is a very useful technique to estimate richness and composition of unsampled areas and have been tested to correctly predict the most of the species’ distributional potential [48, 55-58], for example in applications to invasive species [57-58], tree species [59-60], squamate species [48, 61], and so on.
We included a total of eighteen environmental variables in the model. Variables for details, descriptions, and files for download are described in the following text. We set several optimization parameters while running the software following [48]. The parameters included: 20 runs, 0.001 convergence limit, and 1,000 maximum interactions; rule types: atomic, range, negated range, and logistic regression; best subset active, 5% omission error, 40% commission error, and 67% of points for training; omission measure = extrinsic, and omission threshold = hard; 10 models under hard omission threshold.
The estimation output of DesktopGarp produced in Arc/Info grid maps with ‘zeros’, where the species were not predicted to occur, and ‘ones’, where the species were predicted to occur. The area covered by the coincidence of at least seven out of the 10 models in the best subset selection (optimum models considering omission/commission relationships [62]) were used as the predicted distribution of each species. By doing so and by setting the commission error to 40%, this approach added a component of conservatism in predicting distribution by GARP, which might otherwise extrapolate too much and predict areas that are too far from where the species have previously been collected [48]. After generating such maps using the same criteria for all 151 species, we used ARCGIS software to overlay all species prediction maps into a composite map. This final map was used to create a girded of species richness map at a resolutions of 100 km (approximately equivalent to 1°at the equator) on an Albers Equal-Area Conic projection. Consequently, we used the occurrences of 151 lizard species within 827 grid cells to calculate species richness, summing the value of overlaid corresponding grid cells.
3.2.3. Environmental data
We used eighteen environmental variables. We selected these variables based on previous studies and the four associated hypotheses [40-43, 49]. All environmental variables for assessing hypothesized explanations of species richness were re-projected and re-sampled to the same equal-area cell as the species richness data in ARCGIS. The hypotheses and their related variables are:
Ambient energy—five variables are associated this hypothesis within each cell, including: mean annual potential evapotranspiration (PET) ([63], 30'resolution, available at http://www.grid.unep.ch/data/grid/gnv183.html); mean annual highest temperature (HT), and mean annual lowest temperature (LT) (data from 1961 to 1990 with 1 km2 resolution, available at http://www.data.ac.cn/ index.asp); mean annual sum of effective temperature (≥0℃) (SET0) and mean annual sum of effective temperature (≥10℃) (SET10) (data from 1981 to 1996 with 500 m2 resolution, available at http://www.geodata.cn/Portal).
Productive energy—three variables are used to account for productive energy hypothesis, including: mean annual remotely sensed Normalized Difference Vegetation Index (NDVI), obtained from Advanced Very High-Resolution Radiometer (AVHRR) record of monthly changes in the photosynthetic activity of terrestrial vegetation (data from 1998 to 2008 with 1 km2 resolution, Data source: Environment and Ecology Scientific Data Center of western China, National Natural Science Foundation of China, available at http://westdc.westgis.ac.cn), mean annual actual evapotranspiration (AET)( [60], 30' resolution, available at http://www.grid.unep.ch/data/grid/gnv183.html), and mean annual solar radiation (RAD) (data from 1950 to 1980 with 1 km2 resolution, available at http://www.geodata.cn/Portal).
Contemporary climate hypothesis—eight variables are associated with this hypothesis within each cell, including: mean annual temperature (AT) (data from 1961 to 1990 with 1 km2 resolution, available at http://www.data.ac.cn/index.asp); mean annual sunshine (SUN) (percent of daylength), mean annual diurnal temperature range (DTR) and mean annual frost-day frequency (FF) (data from 1961 to 1990 with 10' resolution [64]); and mean annual wind speed (WIND) (data from 1981 to 1996 with 500 m2 resolution, available at http://www.geodata.cn/Portal); mean annual precipitation (PRE) (data from 1961 to 1990 with 1 km2 resolution, available at http://www.data.ac.cn/index.asp), mean annual wet-day frequency (WET) (number days with >0.1 mm precipitation per month) and mean annual relative humidity (REH) (data from 1961 to 1990 with 10' resolution [64]).
Habitat heterogeneity—the count of 300 m elevation range within each quadrat (ELE) (HYDRO1 k data set for Asia, 1 km2 resolution, available at http://eros.usgs.gov/) and the number of vegetation classes (VEG) (1 km2 resolution, Data source: Environment and Ecology Scientific Data Center of western China, National Natural Science Foundation of China, available at http://westdc.westgis.ac.cn) as indicators of habitat heterogeneity.
3.2.4. Statistical analyses
In order to examine the potential predictors of lizard richness patterns in China, we first tested the relationship between lizard richness and environmental variables using a multiple regression analysis. We did not use all environmental variables employed to run GARP, because including many highly correlated variables in a multiple regression creates several theoretical and statistical problems, especially estimating partial regression coefficients [65]. We selected variables previously identified as affecting species richness and were not highly correlated (r<0.80) and there were one variable represented each hypothesis at least.
We used the eigenvector-based filtering, or spatial eigenvector mapping (SEVM) obtained by Principal Coordinates Neighbour Matrices (PCNM) to account for spatial autocorrelation [66]. Spatial autocorrelation is a potential problem when work with large-scale ecological data and explanatory variables [67-68]. Failure to account for spatial autocorrelation could result in inflating Type I error because model fitting may generate artificially narrow standard errors due to the lack of independence among residuals [67-68]. A truncation distance of 102.33 km, calculated in SAM—Spatial Analysis in Macroecology [69], was used to create the spatial filters. Eigenvector filters were chosen when their influence on species richness was both statistically significant (P<0.05) and had sufficient explanatory power (r2>0.02). We selected eigenvector filters in an iterative process, by minimizing both the spatial autocorrelation among residuals and the number of filters used in regression. Moran’s
To test which hypothesis best explains variation in lizard richness in China, we conducted separate regressions to fit each of the hypothesis, with an addition of mixed models using all variables associated with each hypothesis. The sample-size-corrected Akaike information criterion (AICc) was used to evaluate the goodness of model fit. The model with the lowest AICc score was considered the most parsimonious, therefore optimizing the tradeoff between bias and precision in model construction [71]. The difference between any candidate models and the best model (∆AICc) was used to evaluate the relative model fit when their AICc scores were close. The larger the ∆AICc, the less possible is the fitted model as being the best approximating model in the given models set. In general, Models having ∆AICc ≤ 2 have substantial support (evidence), those in which 4 ≤ ∆AICc ≤7 have considerably less support, and models having ∆AICc≥10 have essentially no support [72]. Model-averaging of estimates using Akaike weights (wi) was used to confront model selection uncertainty [73].
Finally, to make a comparison between the actual available data and the lizard distributions predicted by ecological niche modeling, we mapped species locality points’ data and calculated species richness at 100 km resolution (Figure 4a). This method allowed us to check whether a spatial sampling bias was shown in the final modeling map (i.e., areas that have more species collected coincide with the areas the model indicated as higher species richness) [48].
3.3. Results
3.3.1. Patterns of species richness
Figure 4(b) shows the distribution map which summed all 151 lizard species richness in China. Lizard species richness varied between 1 and 38 species per cell (mean: 13 ± 8 SD) and displayed a consistent pattern that species number increased from higher latitudes to lower ones, and from west to east. The Highest species richness occurs in the Oriental Realm tropics, around the border between southern China and southwestern China, and around the Nanling Mountains, at the border between southern China and central China. Other areas with relative high richness included southwestern China, southern China, northwestern and eastern central China of the Oriental realm (Figure 4b). The raw data map shows a slight sampling bias to the northwest central China and northeast southwestern China (Figure 4a), where the largest herpetological museum CIB/CAS in China are located. However, the niche modeling results are not highly influenced by this bias, since areas corresponding to the highest species richness in China do not overlay completely with the pattern. Furthermore, high richness areas were found by the niche modeling in the relatively poor sampling regions, such as the northwest Mongolia-Xinjiang China, west-south Qinghai-Tibet China areas (Figure 4b).
3.3.2. Species richness and environmental variables
The environmental models for the multiple regression analysis using SEVM with adding eigenvector spatial filters (PCNM), were sufficient to reduce autocorrelation in the residuals (filters data not shown). A spatial correlogram based on Moran’s
Based on the multiple regression analysis taking PCNM spatial filters into account, annual frost-day frequency (FF), elevation range (ELE), the number of vegetation classes (VEG), and wet-day frequency (WET) were the best predictors of species richness. FF was negatively correlated with lizard richness, while ELE, VEG, and WET, respectively, were positively correlated with lizard species richness. Based on model selection approach, the model with the lowest AICc value was the mixed model, which contained all variables related to the different hypotheses. It had an Akaike weight of 1.00. Other models had high ∆AICc values (>10, [69]) and low values of Akaike weights.
3.4. Conclusion
Our results indicate that mechanisms related to different ecological hypotheses might work together to account for lizard richness in China. It is important to consider the influence that environmental factors may have on shaping richness patterns [73]. The frost-day frequency, elevation range, vegetation and wet-day frequency were the most important environmental variable predicting lizard species richness in China. The current alternative hypotheses are not mutually exclusive and may work together and best explain patterns of lizard species richness in China. Based on results of the model selection, our conclusion is in concordant with several previous studies that multiple hypotheses may best account for species richness patterns [48,66,74-75]. Clearly, a variety of factors works synergistically to determine species richness patterns. Our results indicate significant conservation implications, and habitat heterogeneity would be taken into account as an assessment of the threat to endemism from habitat loss in the future investigation. Lizards in China might have experienced large radiations and adapted to dramatic climatic fluctuations after the uplifting of the Tibetan Plateau in Pleistocene. For future studies, it is important to test species richness distribution in Asia at different spatial extent and sample resolution [76], as well as to explore other factors known to affect species richness, such as historical factors and biotic interactions (e.g., competition, predation, and parasitism).
4. Climate and history explain the origin and diversification for toad-headed agamas (genus Phrynocephalus ) and racerunner lizards (genus Eremias )
Dramatic geologic events and climatic shifts are often considered to have significantly influenced the diversification and distributions of organisms. Central Asia has a long history of aridity, with the onset of desertification starting at least 22 million years ago (Ma) [77]. It is believed that the Miocene retreat of the Paratethys Sea, an epicontinental sea stretching over Eurasia 30 Ma, and the uplift of the Tibetan Plateau during the Oligocene–Miocene, played a major role in the shift of the Central Asian climate from oceanic to continental [78], leading to increasing levels of aridity. These climatic changes intensified in the Late Miocene and into the Plio-Pleistocene, as part of the global deterioration of the Cenozoic climates [77]. The Quaternary glaciations in the Qinghai-Tibetan Plateau (QTP) and the bordering mountains were the consequence of a combination between climate and local tectonic uplift. In particular, the Kunlun-Huanghe and Gonghe tectonic uplifts have played very important roles in triggering glaciations in high Asia [ [79-81]. The dramatic geological and climatic histories on the Plateau during the Quaternary have a remarkable influence on regional and adjacent biogeographic patterns [82].
4.1. Toad-headed agamas (Phrynocephalus )
The genus
Until recently, our knowledge of the phylogenetic relationships and historical biogeography of
Most recent analyses use relaxed-clock methods, which allow evolutionary rates to vary among genes and lineages. Guo and Wang [92] re-analyzed the data of Pang
A more recent study [93] using a ~1200 bp region of mitochondrial DNA (ND2-tRNATyr) and a ~1200 bp nuclear gene (RAG-1) and samples across Central Asian agamids results in two well supported conclusions: (i) the onset of aridification in Central Asia during the Late Oligocene, resulting from the retreat of the Paratethys Sea and the intensified uplift of the Tibetan-Himalayan complex, appears to have played an important role in
Dunayev
The topographic variation of the Qinghai-Tibetan Plateau, coupled with cyclical climatic changes in the Pleistocene, and alternating glacial-interglacial periods exerted the greatest influence on the current spatial distribution and genetic structure of viviparous toad-headed lizards. Jin
Using data from 11 microsatellite DNA loci, Wang
Jin and Liu [99] described the phylogeography of a unique endemic agamid lizard,
An aridification of the Tarim Basin and adjacent areas since middle Pleistocene has produced significant genetic structuring of the local fauna. Zhang
4.2. Racerunner lizards (Eremias )
The lacertid genus
Until recently, our understanding of the phylogenetic relationships and historical biogeography of
Wan
Orlova
More recently, a study [111 using combined DNA data sets (3925 bp) from two mitochondrial genes (cyt
Another phylogenetic study [112] using mitochondrial 16S rRNA segment results in five major conclusions: (i) monophyly of
Rastegar-Pouyani
4.3. Conclusion
As exemplified with toad-headed agamas (genus
Appendix S1
List of the 422 squamates species.
Acknowledgments
Support was provided by National Science Foundation of China (30700062, 31272281) and Knowledge Innovation Program of the Chinese Academy of Sciences (KSCX2-EW-Q-6).References
- 1.
Wiens JJ, Donoghue MJ. Historical biogeography, ecology, and species richness. Trends in Ecology and Evolution 2004; 19(12): 639-644. http://www.sciencedirect.com/science/article/pii/S0169534704002745 - 2.
Crisci JV, Sala OE, Katinas L, Posadas, P. Bridging historical and ecological approaches in biogeography. Australian Systematic Botany 2006; 19(1): 1–10. http://www.publish.csiro.au/view/journals/dsp_journal_fulltext.cfm?nid=150&f=SB05006 - 3.
Kent M. Numerical classification and ordination methods in biogeography. Progress in Physical Geography 2006; (3): 399-408. http://ppg.sagepub.com/content/30/3/399.full.pdf - 4.
Posadas P, Crisci JV, Katinas L. Historical biogeography: A review of its basic concepts and critical issues. Journal of Arid Environments 2006; 66(3): 389-403. http://www.sciencedirect.com/science/article/pii/S0140196306000346 - 5.
Meng k, Li S, Murphy RW. Biogeographical patterns of Chinese spiders (Arachnida: Araneae) based on a parsimony analysis of endemicity. Journal of Biogeography 2008; 35(7): 1241-1249. http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2699.2007.01843.x/pdf. - 6.
Hsü KJ, Chen H. Geological Atlas of China: An Application of the Tectonic Facies Concept to the Geology of China. New York: Elsevier Science; 1999. - 7.
Wan T. The Outline of Geotectonics of China. Beijing: Geological Publishing House; 2003. - 8.
Chinese Academy of Sciences, Nature Division Committee. China Animal Geography Division and China Insect Geography Division (First Draft). Beijing: Science Press; 1959. - 9.
Zhang R. Zoogeography of China. Beijing: Science Press; 1999. - 10.
MacKinnon J, Meng S, Cheung C, Carey G, Zhu X, Melville D. A Biodiversity Review of China. Hong Kong: World Wide Fund for Nature International; 1996. - 11.
Xie Y, Mackinnon J, Li D. Study on biogeographical divisions of China. Biodiversity and Conservation 2004; 13(7): 1391-1417. http://www.springerlink.com/content/p1355416n7h87706/ - 12.
Chen L, Song Y, Xu S. The boundary of Palaearctic and Oriental realms in western China. Progress in Natural Science 2008; 18(7): 833-841. http://www.sciencedirect.com/science/article/pii/S1002007108001457 - 13.
Morrone JJ, Escalante T. Parsimony analysis of endemicity (PAE) of Mexican terrestrial mammals at different area units: when size matters. Journal of Biogeography 2002; 29(8):1095-1104. http://onlinelibrary.wiley.com/doi/10.1046/j.1365-2699.2002.00753.x/full - 14.
Humphries CJ, Parenti LR. Cladistic Biogeography, 2nd edn. Oxford Biogeography Series 12, Oxford: Oxford University Press; 1999. - 15.
Rosen BR. From Fossils to Earth History: Applied Historical Biogeography. In: Myers AA, Giller PS. (eds). Analytical Biogeography: An Integrated Approach to the Study of Animal and Plant Distributions. London: Chapman and Hall; 1988. p437-481. - 16.
Morrone JJ. On the identification of areas of endemism. Systematic Biology 1994; 43(3): 438-441. http://www.jstor.org/stable/2413679 - 17.
Cracraft J. Patterns of diversification within continental biotas: hierarchical congruence among the areas of endemism of Australian vertebrates. Australian Systematic Botany 1991; 4(1): 211-227. http://www.publish.csiro.au/paper/SB9910211.htm - 18.
Trejo-Torres JC, Ackerman JD. Biogeography of the Antilles based on a parsimony analysis of orchid distributions. Journal of Biogeography 2001; 28(6): 775-794. http://onlinelibrary.wiley.com/doi/10.1046/j.1365-2699.2001.00576.x/full - 19.
Ron SR. Biogeographic area relationships of lowland Neotropical rainforest based on raw distributions of vertebrate groups. Biological Journal of the Linnean Society 2000; 71(3): 379-402. http://onlinelibrary.wiley.com/doi/10.1111/j.1095-8312.2000.tb01265.x/pdf - 20.
Brooks DR, van Veller MGP. Critique of parsimony analysis of endemicity as a method of historical biogeography. Journal of Biogeography 2003; 30(6): 819-825. http://onlinelibrary.wiley.com/doi/10.1046/j.1365-2699.2003.00848.x/full - 21.
Nihei SS. Misconceptions about parsimony analysis of endemicity. Journal of Biogeography 2006; 33(12): 2099-2106. http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2699.2006.01619.x/full - 22.
Contreras-Medina R, Vega IL, Morrone JJ. Application of parsimony analysis of endemicity to Mexican gymnosperm distributions: grid-cells, biogeographical provinces and track analysis. Biological Journal of the Linnean Society 2007; 92(3): 405-417. http://onlinelibrary.wiley.com/doi/10.1111/j.1095-8312.2007.00844.x/full - 23.
Melo Santos AM, Cavalcanti DR, da Silva JC, Tabarelli M. Biogeographical relationships among tropical forest in north-eastern Brazil. Journal of Biogeography 2007; 34(3): 437-446. http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2699.2006.01604.x/full - 24.
Posadas P. Distributional patterns of vascular plants in Tierra del Fuego: a study applying parsimony analysis of endemicity (PAE). Biogeographica 1996; 72: 161-177. - 25.
Zhao E, Zhao K, Zhou K. Fauna Sinica, Reptilia, Vol. 2. (Squamates, Lacertilia). Beijing: Science press; 1999. - 26.
Zhao E, Adler K. Herpetology of China. Oxford: Ohio: Society for the Study of Amphibians and Reptiles in cooperation with Chinese Society for the Study of Amphibians and Reptiles; 1993. - 27.
Zhao E, Huang M, Zong Y. Fauna Sinica Reptilia, Vol.3 Squamates Serpentes. Beijing: Science press; 1999. - 28.
Zhao E. Snake of China. Hefei: Anhui Science and Technology Publishing House; 2006. - 29.
Swofford DL. PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. Sinauer Associates, Sunderland, MA; 1998. http://www.sinauer.com/detail.php?id=8060 - 30.
Glasby CJ, Álvarez B. Distribution patterns and biogeographic analysis of Austral Polychaeta (Annelida). Journal of Biogeography 1999; 26(3): 507–533. Http://onlinelibrary.wiley.com/doi/10.1046/j.1365-2699.1999.00297.x/pdf - 31.
Ji X, Brana F. Among clutch variation reproductive output and egg size in the wall lizard ( Podarcis muralis ) from a lowland population of Northern Spain. Journal of Herpetology 2000; 34(1): 54-60. http://www.jstor.org/stable/1565238 - 32.
Ji X, Du W. Effects of thermal and hydric environments on incubating eggs and hatchling traits in the cobra, Naja naja atra . Journal of Herpetology 2001; 35(2):186-194. www.jstor.org/stable/1566107 - 33.
Fu C, Wang J, Pu Z, Zhang S, Chen H, Zhao B, Chen J, Wu J. Elevational gradients of diversity for lizards and snakes in the Hengduan Mountains, China. Biodiversity and Conservation 2007; 16(3): 707-726. http://www.springerlink.com/article/10.1007%2Fs10531-005-4382-4 - 34.
Huang Y, Dai Q, Chen Y, Wan H, Li J, Wang Y. Lizard species richness patterns in China and its environmental associations. Biodiversity and Conservation 2011; 20(7):1399-1414. http://www.springerlink.com/ article/10.1007/s10531-011-0033-0# - 35.
Cai B, Huang Y, Chen Y, Hu J, Guo X, Wang Y. Geographic patterns and ecological factors correlates of snake species richness in China. Zoological Research 2012; 33(4): 343-353. http://www.zoores.ac.cn/EN/abstract/abstract3193.shtml# - 36.
Gaston KJ. Global patterns in biodiversity. Nature 2000; 405(6783): 220-227. http://www.nature.com/nature/journal/v405/n6783/pdf/405220a0.pdf - 37.
von Humboldt A. Ansichten der Natur mit wissenschaftlichen Erlauterungen. Germany: J. G. Cotta, Tübingen; 1808. - 38.
Hawkins BA, Field R, Cornell HV, Currie DJ, Guegan JF, Kaufman DM, Kerr JT, Mittelbach GG, Oberdorff T, O’Brien EM, Porter EE, Turner JRG. Energy, water, and broad-scale geographic patterns of species richness. Ecology 2003; 84(12): 3105-3117. http://www.esajournals.org/doi/full/10.1890/03-8006 - 39.
Terribile LC, Olalla-Tárraga MÁ, Morales-Castilla I, Rueda M, Vidanes RM, Rodríguez MÁ , Diniz-Filho JAF. Global richness patterns of venomous snakes reveal contrasting influences of ecology and history in two different clades. Oecologia 2009; 159(3): 617-626. http://www.springerlink.com/article/10.1007%2Fs00442-008-1244-2 - 40.
Currie DJ. Energy and large-scale patterns of animal- and plant-species richness. The American Naturalist 1991; 137(1): 27-49. http://www.jstor.org/stable/2462155 - 41.
Kerr JT, Packer L. Habitat heterogeneity as a determinant of mammal species richness in high-energy regions. Nature 1997; 385(6613): 252-254. http://www.nature.com/nature/journal/v385/n6613/abs/385252a0.html - 42.
Kerr JT, Southwood TRE, Cihlar J. Remotely sensed habitat diversity predicts butterfly species richness and community similarity in Canada. Proceedings of the National Academy of Sciences of the United States of America 2001; 98(20): 11365-11370. http://www.pnas.org/content/98/20/11365.full.pdf+html - 43.
Rahbek C, Graves GR. Multiscale assessment of patterns of avian species richness. Proceedings of the National Academy of Sciences of the United States of America 2001; 98(8):4534-4539. http://www.pnas.org/content/98/8/4534.full.pdf+html - 44.
Turner JRG, Gatehouse CM, Corey CA. Does solar energy control organic diversity? Butterflies, moths and British climate. Oikos 1987; 48(2): 195-205. http://www.jstor.org/stable/3565855 - 45.
Hawkins BA, Porter EE. Water-energy balance and the geographic pattern of species richness of western Palearctic butterflies. Ecological Entomology 2003; 28(6): 678-686. http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2311.2003.00551.x/full - 46.
Mittelbach GG, Steiner CF, Scheiner SM, Gross KL, Reynolds HL, Waide RB, Willig MR, Dodson SI, Gough L. What is the observed relationship between species richness and productivity? Ecology 2001; 82(9): 2381-2396. http://www.esajournals.org/doi/full/10.1890/0012-9658%282001%29082%5B2381%3AWITORB%5D2.0.CO%3B2 - 47.
Wright DH. Species–energy theory: An extension of species–area theory. Oikos 1983; 41(3): 496-506. http://www.jstor.org/stable/3544109 - 48.
Costa GC, Nogueira C, Machado RB, Colli GR. Squamate richness in the Brazilian Cerrado and its environmental-climatic associations. Diversity and Distribution 2007; 13(6): 714-724. http://onlinelibrary.wiley.com/doi/10.1111/j.1472-4642.2007.00369.x/pdf - 49.
Tognelli MF, Kelt DA. Analysis of determinants of mammalian species richness in South America using spatial autoregressive models. Ecography 2004; 27(4): 427-436. http://onlinelibrary.wiley.com/doi/10.1111/j.0906-7590.2004.03732.x/full - 50.
Hortal J, Triantis KA, Meiri S, Thébault E, Sfenthourakis S. Island species richness increases with habitat diversity. The American Naturalist 2009; 174(6): E205-E217. http://www.jstor.org/stable/10.1086/645085 - 51.
McCain CM. Global analysis of reptile elevational diversity. Global Ecology and Biogeography 2010; 19 (4): 541–553. http://onlinelibrary.wiley.com/doi/10.1111/j. 1466-8238.2010.00528.x/pdf - 52.
Qian H, Wang XH, Wang SL, Li YL. Environmental determinants of amphibian and reptile species richness in China. Ecography 2007; 30(4): 471-482. http://onlinelibrary.wiley.com/doi/10.1111/j.0906-7590.2007.05025.x/full - 53.
Pounds JA, Fogden MPL, Campbell JH. Biological response to climate change on a tropical mountain. Nature 1999; 398(6728) : 611-615. http://www.nature.com/nature/journal/v398/n6728/full/398611a0.html - 54.
Stockwell D, Peters D. The GARP modelling system: problems and solutions to automated spatial prediction. International Journal of Geographical Information Science 1999; 13(2): 143-158. http://www.tandfonline.com/doi/abs/10.1080/136588199241391 - 55.
Pearson RG, Raxworthy C, Nakamura M, Peterson AT. Predicting species’ distributions from small numbers of occurrence records: A test case using cryptic geckos in Madagascar. Journal of Biogeography 2007; 34(1): 102-117. http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2699.2006.01594.x/full - 56.
Peterson AT, Papefl M, Eaton M. Transferability and model evaluation in ecological niche modeling: A comparison of GARP and Maxent. Ecography 2007; 30(4): 550-560. http://onlinelibrary.wiley.com/doi/10.1111/j.0906-7590.2007.05102.x/full - 57.
Peterson AT, Robins CR. Using ecological-niche modeling to predict barred owl invasions with implications for spotted owl conservation. Conservation Biology 2003; 17(4): 1161-1165. http://onlinelibrary.wiley.com/doi/10.1046/j.1523-1739.2003.02206.x/full - 58.
Peterson AT, Papeß M, Kluza DA. Predicting the potential invasive distributions of four alien plant species in North America. Weed Science 2003; 51(6): 863-868. http://www.bioone.org/doi/full/10.1614/P2002-081 - 59.
Ferreira de Siqueira M, Durigan G, de Marco Ju´nior P, Peterson AT. Something from nothing: using landscape similarity and ecological niche modeling to find rare plant species. Journal for Nature Conservation 2009; 17(1): 25-32. http://www.sciencedirect.com/science/article/pii/S1617138108000484 - 60.
Menon S, Choudhury BI, Khan LM, Peterson AT. Ecological niche modeling and local knowledge predict new populations of Gymnocladus assamicus a critically endangered tree species. Endangered Species Research 2010; 11: 175-181. http://www.int-res.com/articles/esr2010/11/n011p175.pdf - 61.
Raxworthy C, Martínez-Meyer E, Horning N, Nussbaum R, Schneider G, Ortega-Huerta M, Peterson AT. Predicting distributions of known and unknown reptile species in Madagascar. Nature 2003; 426(6968): 837-841. http://www.nature.com/nature/journal/v426/n6968/full/nature02205.html - 62.
Anderson RP, Lew D, Peterson AT. Evaluating predictive models of species’ distributions: criteria for selecting optimal models. Ecological Modelling 2003; 162(3): 211-232. http://www.sciencedirect.com/science/article/pii/S0304380002003496 - 63.
Ahn CH, Tateishi R. Development of a global 30-minute grid potential evapotranspiration data set. Journal of the Japan Society Photogrammetry and Remote Sensing 1994; 33(2): 12-21. - 64.
New M, Lister D, Hulme M, Makin I. A high-resolution data set of surface climate over global land areas. Climate Resarch 2002; 21: 1-25. http://www.int-res.com/articles/cr2002/21/c021p001.pdf - 65.
Tabachnick BG, Fidell LS. Using Multivariate Statistics, 4th edn. Boston: Allyn and Bacon, MA; 2000. - 66.
Diniz-Filho JAF, Bini LM. Modelling geographical patterns in species richness using eigenvectorbased spatial filters. Global Ecology and Biogeography 2005; 14(2): 177-185. http://onlinelibrary.wiley.com/doi/10.1111/j.1466-822X.2005.00147.x/full - 67.
Legendre P. Spatial autocorrelation—trouble or new paradigm. Ecology 1993; 74(6): 1659-1673. http://www.jstor.org/stable/1939924 - 68.
Legendre P, Dale MRT, Fortin MJ, Gurevitch J, Hohn M, Myers D. The consequences of spatial structure for the design and analysis of ecological field surveys. Ecography 2002; 25(5): 601-615. http://onlinelibrary.wiley.com/doi/10.1034/j.1600-0587.2002.250508.x/full - 69.
Rangel TF, Diniz-Filho JAF, Bini LM. SAM: A comprehensive application for spatial analysis in macroecology. Ecography 2010; 33(1): 46-50. http://onlinelibrary.wiley.com/doi/10.1111/j.1600-0587.2009.06299.x/full - 70.
Diniz-Filho JAF, Bini LM, Hawkins BA. Spatial autocorrelation and red herrings in geographical ecology. Global Ecology and Biogeography 2003; 12(1): 53-64. http://onlinelibrary.wiley.com/doi/10.1046/j.1466-822X.2003.00322.x/full - 71.
Burnham KP, Anderson DR. Model Selection and Inference: A Practical Information Theoretic Approach. New York: Springer-Verlag; 1998. - 72.
Burnham KP, Anderson DR. Multimodel inference—understanding AIC and BIC in model selection. Sociological Methods and Research 2004; 33(2): 261-304. http://smr.sagepub.com/content/33/2/261.full.pdf+html - 73.
Powney GD, Grenyer R, Orme CDL, Owens IPF, Meiri S. Hot, dry and different Australian lizard richness is unlike that of mammals, amphibians, and birds. Global Ecology and Biogeography 2010; 19(3): 386-396. http://onlinelibrary.wiley.com/doi/10.1111/j.1466-8238.2009.00521.x/full - 74.
Andrews P, O’Brien EM. Climate, vegetation, and predictable gradients in mammal species richness in southern Africa. Journal of Zoology 2000; 251(20): 205–231. http://onlinelibrary.wiley.com/doi/10.1111/j.1469-7998.2000.tb00605.x/pdf - 75.
Bohning-Gaese K. Determinants of avian species richness at different spatial scales. Journal of Biogeography 1997; 24(1): 49-60. http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2699.1997.tb00049.x/pdf - 76.
Willis KJ, Whittaker RJ. Species diversity–scale matters. Science 2002; 295(5558): 1245-1248. http://www.sciencemag.org/content/295/5558/1245.full - 77.
Guo Z, Ruddiman WF, Hao Q, Wu H, Qiao Y, Zhu R, Peng S, Wei J, Yuan B, Liu T. Onset of Asian desertification by 22 Myr ago inferred from loess deposits in China. Nature 2002; 416(6877): 159-163. http://www.nature.com/nature/journal/v416/n6877/full/416159a.html - 78.
Ramsetin G, Fluteau F, Besse J, Joussaumes S. Effect of orogeny, plate motion and land-sea distribution on Eurasian climate change over the past 30 million years. Nature 1997; 386: 788-795. http://www.nature.com/nature/journal/v386/n6627/abs/386788a0.html - 79.
Shi Y, Zhao J, Wang J. New Understanding of Quaternary Glaciations in China. Shanghai: Shanghai Popular Science Press; 2011. - 80.
Zhao J, Shi Y, Wang J. Comparison between Quaternary glaciations in China and the marine oxygen isotope stage (MIS): An improved schema. Acta Geographica Sinica 2011; 66(7): 867-884. - 81.
Zhou S, Li J, Zhao J, Wang J, Zheng J. Quaternary Glaciations: Extent and Chronology in China. In: Ehlers J, Gibbard PL, Hughes PD. (eds) Quaternary Glaciations-Extent and Chronology, Volume 15: A closer look (Developments in Quaternary Science). Amsterdam: Elsevier Press; 2011. p981-1002. Available from http://www.sciencedirect.com/science/article/pii/B9780444534477000702 - 82.
Yang S, Dong H, Lei F. Phylogeogrpahy of regional fauna on the Tibetan Plateau: A review. Progress in Natural Science 2009; 19(7): 789–799. http://www.sciencedirect.com/science/article/pii/S1002007109000513 - 83.
Barabanov AV, Ananjeva NB. Catalogue of the available scientific species-group names for lizards of the genus Phrynocephalus Kaup, 1825 (Reptilia, Sauria, Agamidae). Zootaxa 2007; 1399: 1-56. http://www.mapress.com/zootaxa/2007f/z01400p068f.pdf - 84.
Uetz P. The Reptile Database, http://www.reptile-database.org (accessed Aug 3, 2012). - 85.
Whiteman RS. Evolutionary history of the lizard genus Phrynocephalus (Lacertilia, Agamidae). Master of Art thesis. California State University, Fullerton; 1978. - 86.
Arnold EN. Phylogenetic relationships of toad-headed lizards ( Phrynocephalus , Agamidae) based on morphology. Bulletin of the British Museum (Natural History) Zoology 1999; 65: 1-13. - 87.
Ananjeva NB, Tuniyev BS. Historical biogeography of the Phrynocephalus species of the USSR. Asiatic Herpetology Research 1992; 4: 76-98. http://www.asiatic-herpetological.org/Archive/Volume%2004/04_12.pdf - 88.
Macey JR, Ananjeva NB, Zhao E, Wang Y, Papenfuss TJ. An allozyme-based phylogenetic hypothesis for Phrynocephalus (Agamidae) and its implications for the historical biogeography of arid Asia. In: Zhao E, Chen B, Papenfuss TJ. (eds.) Proceedings of the First Asian Herpetological Meeting. Beijing: China Forestry Press; 1993. - 89.
Wang Y, Macey JR. On the ecological-geographic differentiation of Chinese species of the genus Phrynocephalus . In: Zhao E, Chen B, Papenfuss TJ. (eds.) Proceedings of the First Asian Herpetological Meeting. Beijing: China Forestry Press; 1993. - 90.
Zeng X, Wang Y, Liu Z, Fang Z, Wu G, Papenfuss TJ, Macey JR. Karyotypes of nine species in the genus Phrynocephalus , with discussion of karyotypic evolution of ChinesePhrynocephalus . Acta Zoologica Sinica 1997; 43(4): 409-418. - 91.
Pang J, Wang Y, Zhong Y, Hoelzel AR, Papenfuss TJ, Zeng X, Ananjeva NB, Zhang YP. A phylogeny of Chinese species in the genus Phrynocephalus (Agamidae) inferred from mitochondrial DNA sequences. Molecular Phylogenetics and Evolution 2003; 27(3): 398-409. http://www.sciencedirect.com/science/article/pii/S1055790303000198 - 92.
Guo X, Wang Y. Partitioned Bayesian analyses, dispersalvicariance analysis, and the biogeography of Chinese toad-headed lizards (Agamidae: Phrynocephalus ): A re-evaluation. Molecular Phylogenetics and Evolution 2007; 45(2): 643-662. http://www.sciencedirect.com/science/article/pii/S1055790307002175 - 93.
Melville J, Hale J, Mantziou G, Ananjeva NB, Milto K, Clemann N. Historical biogeography, phylogenetic relationships and intraspecific diversity of agamid lizards in the Central Asian deserts of Kazakhstan and Uzbekistan. Molecular Phylogenetics and Evolution 2009; 53(1): 99-112. http://www.sciencedirect.com/science/article/pii/S1055790309001754 - 94.
Solovyeva EN, Poyarkov NA, Dunayev EA, Duysebayeva TN, Bannikova AA. Molecular differentiation and taxonomy of the sunwatcher toad-headed agama species complex Phrynocephalus superspecieshelioscopus (Pallas 1771) (Reptilia: Agamidae). Russian Journal of Genetics 2011; 47(7): 842-856. http://www.springerlink.com/content/h245501257785793/ - 95.
Dunayev EA, Ivanova N, Poyarkov NA, Borisenko A, Dujsebayeva T, Hebert PD. Molecular perspective on the evolution and barcoding of toad-headed agamas (genus Phrynocephalus , Agamidae) in Middle Asia. Programme and Abstracts of 14th European Congress of Herpetology and SHE Ordinary General Meeting, 19-23 September 2007, Porto (Portogal); 2007. - 96.
Jin Y, Brown RP, Liu N. Cladogenesis and phylogeography of the lizard Phrynocephalus vlangalii (Agamidae) on the Tibetan plateau. Molecular Ecology 2008; 17(8): 1971-1982. http://onlinelibrary.wiley.com/doi/10.1111/j.1365-294X.2008.03721.x/pdf. - 97.
Wang Y, Zhan A, Fu J. Testing historical phylogeographic inferences with contemporary gene flow data: population genetic structure of the Qinghai toad-headed lizard. Journal of Zoology 2009; 278(2): 149-156. http://onlinelibrary.wiley.com/doi/10.1111/j.1469-7998.2009.00564.x/full - 98.
Guo X, Liu L, Wang Y. Phylogeography of the Phrynocephalus vlangalii species complex in the upper reaches of the Yellow River inferred from mtDNA ND4-tRNALeu segments. Asian Herpetology Research 2012; 3(1): 52-68. http://www.ahr-journal.com/index.php?module=case&act=List&ClassID=19 - 99.
Jin Y, Liu N. Phylogeography of Phrynocephalus erythrurus from the Qiangtang Plateau of the Tibetan Plateau. Molecular Phylogenetics and Evolution 2009; 54(3): 933-940. http://www.sciencedirect.com/science/article/pii/S1055790309004503 - 100.
Zhang Q, Xia L, He J, Wu Y, Fu J, Yang Q. Comparison of phylogeographic structure and population history of two Phrynocephalus species in the Tarim Basin and adjacent areas. Molecular Phylogenetics and Evolution 2010; 57(3): 1091-1104. http://www.sciencedirect.com/science/article/pii/S1055790310004094 - 101.
Fu B, Lin A, Kano K, Maruyama T, Guo J. Quaternary folding of the eastern Tian Shan, northwest China. Tectonophysics 2003; 369(1-2): 79-101. http://www.sciencedirect.com/science/article/pii/S0040195103001379 - 102.
Szczerbak NN. Guide to the Reptiles of the Eastern Palearctic. USA: Krieger Publishing Company, Malabar, Florida; 2003. - 103.
Szczerbak NN. Yashchurki Palearktiki ( Eremias lizards of the Palearctic). Axadeimiya Nauk Ukrainskoi USSR Institut Zoologii. Naukova Dumka, Kiev; 1974. - 104.
Arnold EN. The hemipenis of lacertid lizards (Reptilia: Lacertidae): structure, variation and systematic implications. Journal of Natural History 1986; 20: 1221-1257. http://www.tandfonline.com/doi/abs/10.1080/00222938600770811#preview - 105.
Arnold EN. Towards the phylogeny and biogeography of the Lacertidae: relationships within an Old-World family of lizards derived from morphology. Bulletin of the British Museum (Natural History) Zoology 1989; 55(2): 209-257. http://biostor.org/reference/107021 - 106.
Boulenger GA. Monograph of the Lacertidae, vol. II. London: Trustees of the British Museum of Natural History; 1921. - 107.
Fitzinger L. Systema Reptilium. Vienna: Fasciculus Primus; 1843. - 108.
Szczerbak NN. Taxonomy of the genus Eremias (Sauria, Reptilia) in connection with the focuses of the desert-steppe fauna development in Paleoarctic. Vestnik Zoologii 1971; 2: 48-55. - 109.
Wan L, Sun S, Jin Y, Yan Y, Liu N. Molecular phylogeography of the Chinese lacertids of the genus Eremias (Lacertidae) based on 16S rRNA mitochondrial DNA sequences. Amphibia-Reptilia 2007; 28(1): 33-41.. http://www.ingentaconnect.com/content/brill/amre/2007/00000028/00000001/art00005 - 110.
Orlova VF, Poyarkov NA, Chirikova M, Dolotovskaya S.I. Preliminary molecular phylogeography of wide-spread steppe-runner lizard – Eremias arguta (Lacertidae) and considerations on its subspecific structure. Programme and Abstracts of 14th European Congress of Herpetology and SHE Ordinary General Meeting, 19-23 September 2007, Porto (Portogal); 2007. - 111.
Rastegar-Pouyani E, Joger U, Wink M. A molecular phylogeny of the genus Eremias (Reptillia, Lacertidae) based on the mitochondrial and nuclear DNA sequences. Programme and Abstracts of 15th European Congress of Herpetology and SHE Ordinary General Meeting, 28 September-02 October 2009, Kuşadasi, Turkey; 2009. - 112.
Guo X, Dai X, Chen D, PapenfussTJ, Ananjeva NB, Melnikov DA, Wang Y. Phylogeny and divergence times of some racerunner lizards (Lacertidae: Eremias ) inferred from mitochondrial 16S rRNA gene segments. Molecular Phylogenetics and Evolution 2011; 61(2): 400-412. http://www.sciencedirect.com/science/article/pii/S1055790311003113 - 113.
Zhao Q, Liu HX, Luo L, Ji X. Comparative population genetics and phylogeography of two lacertid lizards ( Eremias argus andE. brenchleyi ) from China. Molecular Phylogenetics and Evolution 2011; 58(3): 478-491. http://www.sciencedirect.com/science/article/pii/S1055790310004987 - 114.
Rastegar-Pouyani E, Noureini SK, Joger U, Wink M. Molecular phylogeny and intraspecific differentiation of the Eremias velox complex of the Iranian Plateau and Central Asia (Sauria, Lacertidae). Journal of Zoological Systematics and Evolutionary Research 2012; 50(3): 220-229. http://onlinelibrary.wiley.com/doi/10.1111/j.1439-0469.2012.00662.x/pdf - 115.
Rastegar-Pouyani E. A Phylogeny of the Eremias velox complex of the Iranian Plateau and Central Asia (Reptilia, Lacertidae): molecular evidence from ISSR-PCR fingerprints. Iranian Journal of Animal Biosystematics 2009; 5(1): 33-46. http://www.sid.ir/en/VEWSSID/J_pdf/116420090105.pdf - 116.
Warren DL, Glor RE, Turelli M. Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution 2008; 62(11): 2868-2883. http://onlinelibrary.wiley.com/doi/10.1111/j.1558-5646.2008.00482.x/full