**Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of Squamates in China**

Yong Huang, Xianguang Guo and Yuezhao Wang

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/55719

### **1. Introduction**

[21] Namzalov, B. B. (1996). The mountain forest-steppe of Southern Siberia landscape phenomenon of Central Asia. The Flora and Vegetation of Siberia and the Far East. Lectures in the memory of L.M.Cherepnin. Book of Abstracts of the nd All-Russia

[22] Peshkova, G. A. (1962). Forest-steppe relationship in the Angara region. Trans. East-

[24] Peshkova, G. A. (1985). The Vegetation of Siberia, Prebaikalia and Transbaikalia, No‐

[25] Peshkova, G. A. (2001). Florogenetic Analysis of the Steppe Flora of the Mountains in

[26] Popov, M. G. (1953). On the relationship of forest (taiga) and steppe in Middle Sibe‐

[27] Popov, M. G. (1957). The steppe and rock floras of the Baikal's western shore area.

[29] Sizykh, A. P. Protection of the taiga-steppe ecosystems of the western shore of Lake Baikal // Journal of Design & Nature and Ecodynamics, UK, (2009). , 4, 66-71.

[30] Snytko, V A. Dan'ko L V, Kuz'min S B, Sizykh A Diversity of geosystems of taiga and steppe contact on the Western coast of Lake Baikal. Geography and Natural Re‐

[31] ShiytovS. GTerent'ev M.M., Fomin V.V. Spatiotemporal Dynamics of Forest-Tundra Communities in the Polar Urals // Russian Journal of Ecology, (2005). , 36, 69-75.

[33] The Vegetation (South of East Siberia)(1972). Map 1: 1 500 000 (ed. by A.V.Belov),

[34] Thomas, T. Veblen and Diane C. Lorenz. (1988). Recent vegetation changes along the forest/steppe ecotone of Northern Patagonia. Annals of the Association of American

[35] Tsybzhitov Ts KhTsybikdorzhiev Ts Ts, Tsybzhitov A I. (1999). Soils of Lake Baikal

[36] VoroninV. I Sizykh A. P Oskolkov V. A Voronin, Sheifer E.V. Structural-dynamics organization of the plant communities of the basic ecotones of the Lake Baikal // En‐

[32] The Types of Forests in the South of East Siberia(1980). Novosibirsk: , 236-243.

[28] Savina, L. N. (1986). Taiga Forests of Northern Asia in the Holocen, Novosibirsk.

[23] Peshkova, G. A. (1972). The Steppe Flora of Baikalian Siberia, Moscow.

ria. Bull. of Moscow Soc. of Naturalists, Biol. Dept. 58, 5: 81-85.

Conference, Krasnoyarsk: 215-217., 2.

Sib. Biol. Inst. , 1, 90-99.

Southern Siberia, Novosibirsk.

Trans. Baik. Limnol. st. , 15, 408-426.

sources: 61-68., 2001.

Geographers. , 78, 93-111.

vironments- 2010, Russia, Tomsk, (2010). , 45-46.

GUGK.

basin: 128.

vosibirsk.

168 Current Progress in Biological Research

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 distribu‐ tion, 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 envi‐ ronmental constraints, whereas the latter focuses on taxonomic groups and historical biogeo‐ graphical 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 *et al*. [5], China has a relatively independent geological history. Six primary plates are involved in its tectonic history, namely the North China, Tarim, Yangtze, Cathaysian plates, and parts of the Siberia and Gondwana plates [6-7]. In the Pleistocene, seven collisions and integrations of these plates had united the ancient Siberian and European plates. With approximately 200 tectonic-facies, China has a complex topography, including towering mountains, basins of various sizes, undulating plateaus and hills, and fertile plains. As noted by Meng *et al*. [5], they are assigned to four different terraces in general. The highest terrace, the Qinghai-Tibetan Plateau – the 'roof of the world' – has an average elevation of over 4000 m. Tremendous differences in latitude (a span of more than 50°latitude), longitude (a span of more than 60°longitude), and altitude (a span of more than 88km) create the conditions for

© 2013 Huang et al.; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2013 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

extremely diverse climate and a highly astounding heterogeneous landscape. Thus, a high degree of species richness and endemism in China as might have been expected in terms of squamates diversity. It is conservatively estimated that today there are about 422 squamates species in this country, belonging to 106 genera in 17 families and two suborders (see Appendix S1 in Supplementary Material).

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

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

171

http://dx.doi.org/10.5772/55719

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 *et al*. [11](Figure 3), in order to find areas of congruent

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, topo‐ graphical 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, vegeta‐ tion, water system, farm belt, and so on [11]. More details about the OGUS see these studies.

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 *Fauna Sinica Reptilia, Vol.2 Squamates Lacertilia* [25], *Herpetology of China* [26], *Fauna Sinica Reptilia, Vol.3 Squamates Serpentes* [27], *Snakes of China (2 Volume set)* [28]*, Zootaxa, Zoological Journal of the Linnean Society,Bonn Zoological Bulletin, Asian Herpetological Research*, *Zoological Research, Acta Herpetologica Information, Acta Herpetologica Sinica, Acta Zootaxonomica Sinica, Sichuan Journal of Zoology*, *Chinese Journal of Zoology*, and so on. For the taxonomic revisions and updated information recently, we follow the taxonomy of Fauna

used in many biogeographical studies in recent years [5,22-23].

distributional patterns in China.

*2.2.1. Study area and operational geographical units*

**2.2. Materials and methods**

*2.2.2. Data sources and dataset*

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 biogeograph‐ ical 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 *et al*. [11](Figure 3), in order to find areas of congruent distributional patterns in China.

#### **2.2. Materials and methods**

extremely diverse climate and a highly astounding heterogeneous landscape. Thus, a high degree of species richness and endemism in China as might have been expected in terms of squamates diversity. It is conservatively estimated that today there are about 422 squamates species in this country, belonging to 106 genera in 17 families and two suborders (see Appendix

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 biogeograph‐ ical 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

**2. Compared biogeographical patterns of Chinese squamates based on parsimony analysis of endemicity (PAE) at different natural area units**

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

S1 in Supplementary Material).

170 Current Progress in Biological Research

pattern [13].

**2.1. Introduction**

#### *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, topo‐ graphical 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, vegeta‐ tion, 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 *Fauna Sinica Reptilia, Vol.2 Squamates Lacertilia* [25], *Herpetology of China* [26], *Fauna Sinica Reptilia, Vol.3 Squamates Serpentes* [27], *Snakes of China (2 Volume set)* [28]*, Zootaxa, Zoological Journal of the Linnean Society,Bonn Zoological Bulletin, Asian Herpetological Research*, *Zoological Research, Acta Herpetologica Information, Acta Herpetologica Sinica, Acta Zootaxonomica Sinica, Sichuan Journal of Zoology*, *Chinese Journal of Zoology*, and so on. For the taxonomic revisions and updated information recently, we follow the taxonomy of Fauna Sinica Reptilia [26-27], Snakes of China [28] and the reptile database (http://www.reptiledatabase.org/). To reduce the sources of potential errors records, uncertain distributions and identifications, such as isolated distribution records, we checked again and removed erroneous records. If possible, we tried to overcome these difficulties by adding as many taxa as possible to the analyses. The distributional data with longitudes and latitudes were reviewed from museum databases or literature, and otherwise, we compiled coordinates from Google Earth. Finally, an initial data set of 422 Squamates with species names (see Appendix S1 in Supple‐ mentary Material) represented by a total of 50,749 records, of these, 61.0% records were lizards, while the remaining records were snakes.

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

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

173

http://dx.doi.org/10.5772/55719

**Figure 1.** (a) The division of 28 areas (A01–A28) in China used in the PAE analysis. (b) 50 % majority consensus tree generated by PAE analysis. The terminals correspond to the OGUs shown in (a). Bootstrap values were shown above

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

support value of 64%. The remaining subclades showed ambiguous.

the branches

*2.3.2. Biological areas 54*

#### *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 autapo‐ morphies [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 *et al* [5]. The parsimony analysis obtained twelve most-parsimonious trees [Tree length = 676; Consis‐ tency index (CI) = 0.3609; Homoplasy index (HI) = 0.6391; Retention index (RI) = 0.6415]. Although two clear groups, D and E, were identified by the analysis, they did not support either the delimitation between palaearctic and oriental realms in China or the geographical barrier within the confines of the Qinling Mountains and Huai River, and east of the Hengduan Mountains [5]. The bootstrap value of group D (64%) was almost the same as that of group E (61%). The group D includes area 21 and area 23, which corresponded to the Tibetan Plateau and Taiwan Island in Map (Figure 1).

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.

**Figure 1.** (a) The division of 28 areas (A01–A28) in China used in the PAE analysis. (b) 50 % majority consensus tree generated by PAE analysis. The terminals correspond to the OGUs shown in (a). Bootstrap values were shown above the branches

#### *2.3.2. Biological areas 54*

Sinica Reptilia [26-27], Snakes of China [28] and the reptile database (http://www.reptiledatabase.org/). To reduce the sources of potential errors records, uncertain distributions and identifications, such as isolated distribution records, we checked again and removed erroneous records. If possible, we tried to overcome these difficulties by adding as many taxa as possible to the analyses. The distributional data with longitudes and latitudes were reviewed from museum databases or literature, and otherwise, we compiled coordinates from Google Earth. Finally, an initial data set of 422 Squamates with species names (see Appendix S1 in Supple‐ mentary Material) represented by a total of 50,749 records, of these, 61.0% records were lizards,

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 autapo‐ morphies [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.

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 *et al* [5]. The parsimony analysis obtained twelve most-parsimonious trees [Tree length = 676; Consis‐ tency index (CI) = 0.3609; Homoplasy index (HI) = 0.6391; Retention index (RI) = 0.6415]. Although two clear groups, D and E, were identified by the analysis, they did not support either the delimitation between palaearctic and oriental realms in China or the geographical barrier within the confines of the Qinling Mountains and Huai River, and east of the Hengduan Mountains [5]. The bootstrap value of group D (64%) was almost the same as that of group E (61%). The group D includes area 21 and area 23, which corresponded to the Tibetan Plateau

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

while the remaining records were snakes.

172 Current Progress in Biological Research

*2.2.3. Parsimony Analysis of Endemicity (PAE)*

**2.3. Results**

*2.3.1. Biological areas 28*

and Taiwan Island in Map (Figure 1).

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%.

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 corre‐ sponded 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

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

175

http://dx.doi.org/10.5772/55719

(Emin Valley) formed Clade AZ with a well support proportion of 97%.

patterns may lead to incongruence in the distributions patterns [30].

Past major geological events have played important roles in shaping the biogeographic distribution of extant organisms. PAE originally aimed to find areas of congruent distribu‐ tional 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

**2.4. Conclusion**

**Figure 2.** (a) The division of 54 areas (A01–A54) in China used in the PAE analysis. (b) 50 % majority consensus tree generated by PAE analysis. The terminals correspond to the OGUs shown in (a). Bootstrap values were shown above the branches

#### *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 corre‐ sponded 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**

areas A50 and A51 with 75% bootstrap value. Clade L contained Taiwan Island areas with

**Figure 2.** (a) The division of 54 areas (A01–A54) in China used in the PAE analysis. (b) 50 % majority consensus tree generated by PAE analysis. The terminals correspond to the OGUs shown in (a). Bootstrap values were shown above

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

the branches

*2.3.3. Biological areas 124*

174 Current Progress in Biological Research

similar Squamates faunas, areas A52 and A53; its well-supported value was 100%.

Past major geological events have played important roles in shaping the biogeographic distribution of extant organisms. PAE originally aimed to find areas of congruent distribu‐ tional 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].

in China, such as Southern China geological province versus Southern regions, The Lauren‐ tian/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 biogeo‐

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

177

http://dx.doi.org/10.5772/55719

**3. Distribution patterns in species diversity of lizards in China and their**

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

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

(environmental) hypotheses have been widely discussed and accepted [39].

graphic patterns.

**3.1. Introduction**

**relationships to ecological factors**

**Figure 3.** (a) The division of 124 areas (A01–A124) in China used in the PAE analysis. (b) 50 % majority consensus tree generated by PAE analysis. The terminals correspond to the OGUs shown in (a). Bootstrap values were shown above the branches

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, evolution‐ ary 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 Lauren‐ tian/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 biogeo‐ graphic patterns.

### **3. Distribution patterns in species diversity of lizards in China and their relationships to ecological factors**

#### **3.1. Introduction**

**Figure 3.** (a) The division of 124 areas (A01–A124) in China used in the PAE analysis. (b) 50 % majority consensus tree generated by PAE analysis. The terminals correspond to the OGUs shown in (a). Bootstrap values were shown above

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, evolution‐ ary 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

the branches

176 Current Progress in Biological Research

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 thor‐ ough. 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.

commission error to 40%, this approach added a component of conservatism in predict‐ ing 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 generat‐ ing 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

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

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

**1.** 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

**2.** 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

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

**3.** Contemporary climate hypothesis—eight variables are associated with this hypothesis within each cell, including: mean annual temperature (AT) (data from 1961 to 1990 with

 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

resolution, available at http://www.data.ac.cn/index.asp), mean annual wet-

annual solar radiation (RAD) (data from 1950 to 1980 with 1 km2

resolution, available at http://www.geodata.cn/Portal).

resolution, Data source: Environment and Ecology Scientific Data

resolution, available at

Squamates in China

179

http://dx.doi.org/10.5772/55719

richness, summing the value of overlaid corresponding grid cells.

*3.2.3. Environmental data*

related variables are:

1981 to 1996 with 500 m2

to 2008 with 1 km2

1 km2

with 1 km2

http://www.geodata.cn/Portal).

#### **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' distribu‐ tional 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 predict‐ ing 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 generat‐ ing 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*

hypotheses. Because their taxonomy is well resolved and distributional data are quite thor‐ ough. 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.

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

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' distribu‐ tional potential [48, 55-58], for example in applications to invasive species [57-58], tree species

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 =

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

**3.2. Materials and methods**

178 Current Progress in Biological Research

*3.2.1. Data collection*

deviation = 38.3).

*3.2.2. Ecological niche modeling and species richness*

[59-60], squamate species [48, 61], and so on.

hard; 10 models under hard omission threshold.

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:


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]).

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

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

181

http://dx.doi.org/10.5772/55719

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

**Figure 4.** Number of species per grid cell based on (a) the raw data of museum collections, and (b) the ecological ni‐ che modeling of 151 species. The grid corresponds to the approximate area of China and the area of each cell is 100

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

km2. Blank cells have no specimen based on the major collections [34].

**3.3. Results**

*3.3.1. Patterns of species richness*

using Akaike weights (wi) was used to confront model selection uncertainty [73].

collected coincide with the areas the model indicated as higher species richness) [48].

**4.** 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 *I* coefficient were used to examine the model residuals of the spatial autocorrelation in reducing spatial autocorrela‐ tion [70]. These filters were then used as candidate predictor variables, together with other environmental predictors formed in the full model. In this way, the effects of environmental predictors are evaluated as partial effects, taking spatial factors into account explicitly [69]. The total explanatory power, r2 values, was divided into three parts: a part explained by space, a part explained by environmental variables, and a part of shared explained variance.

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].

**Figure 4.** Number of species per grid cell based on (a) the raw data of museum collections, and (b) the ecological ni‐ che modeling of 151 species. The grid corresponds to the approximate area of China and the area of each cell is 100 km2. Blank cells have no specimen based on the major collections [34].

#### **3.3. Results**

day frequency (WET) (number days with >0.1 mm precipitation per month) and mean

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.

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

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

eigenvector filters in an iterative process, by minimizing both the spatial autocorrelation among residuals and the number of filters used in regression. Moran's *I* coefficient were used to examine the model residuals of the spatial autocorrelation in reducing spatial autocorrela‐ tion [70]. These filters were then used as candidate predictor variables, together with other environmental predictors formed in the full model. In this way, the effects of environmental predictors are evaluated as partial effects, taking spatial factors into account explicitly [69].

a part explained by environmental variables, and a part of shared explained variance.

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

values, was divided into three parts: a part explained by space,

correlated (r<0.80) and there were one variable represented each hypothesis at least.

statistically significant (P<0.05) and had sufficient explanatory power (r2

resolution, Data source: Environment and

>0.02). We selected

annual relative humidity (REH) (data from 1961 to 1990 with 10' resolution [64]). **4.** 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

*3.2.4. Statistical analyses*

180 Current Progress in Biological Research

The total explanatory power, r2

#### *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).

richness, such as historical factors and biotic interactions (e.g., competition, predation, and

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

183

http://dx.doi.org/10.5772/55719

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 Quater‐ nary 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

The genus *Phrynocephalus*, i.e., toad-headed agamas, are a modern group of lizards including about ~37 distinct species [83, but see 84]. They are typically grouped together with Asian Rock Agamas (*Laudakia*), *Bufoniceps*, and *Trapelus* in the subfamily Agaminae. Toad-headed agamas inhabit arid regions from northwestern China to the western side of the Caspian Sea, across the Tibetan Plateau, and southwest Asia to the Arabian Peninsula. They constitute one of the major components of the central Asian desert fauna and are highly adapted to sand dunes and desert environments. The reproductive biology of *Phrynocephalus* is notable in that there exist two reproductive modes: viviparity and oviparity. All six viviparous species are endemic to China and mainly restricted to high elevations in the Tibetan Plateau. Despite the progress that has been made in recent years for the systematics of certain *Phrynocephalus* groups, the large-scale pattern of their evolution in time and space remains open. Here, we review the origin and diversification of toad-headed agamas with their relationships to geological and

Until recently, our knowledge of the phylogenetic relationships and historical biogeography of *Phrynocephalus* was mainly based on anatomical studies [85-86]. Researchers usually agree that the broad distribution of *Phrynocephalus* may be explained by different hypotheses. Ananjeva and Tuniyev [87] speculated that there were two original centers for the *Phrynoce‐ phalus* species: Central Asia (a northern Tethys origin) and Middle Asia (a southern Tethys

**4. Climate and history explain the origin and diversification for toadheaded agamas (genus** *Phrynocephalus***) and racerunner lizards (genus**

parasitism).

*Eremias***)**

patterns [82].

climatic changes.

**4.1. Toad-headed agamas (***Phrynocephalus***)**

#### *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 *I* index was used to evaluate the pattern of spatial autocorrelation in the residuals of the regression. The multiple regression model explained a total of 80.1% variance of lizard richness in China (r2 = 0.801; F = 203.47; P<0.001). The total explanatory power explained by predictors alone was 17.0%, and explained by space alone was 5.4%, and the shared explained variance was 57.7%.

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 toadheaded 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 Quater‐ nary 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***)**

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

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 *I* index was used to evaluate the pattern of spatial autocorrelation in the residuals of the regression. The multiple regression

P<0.001). The total explanatory power explained by predictors alone was 17.0%, and explained

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

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

= 0.801; F = 203.47;

Mongolia-Xinjiang China, west-south Qinghai-Tibet China areas (Figure 4b).

model explained a total of 80.1% variance of lizard richness in China (r2

by space alone was 5.4%, and the shared explained variance was 57.7%.

*3.3.2. Species richness and environmental variables*

182 Current Progress in Biological Research

(>10, [69]) and low values of Akaike weights.

**3.4. Conclusion**

The genus *Phrynocephalus*, i.e., toad-headed agamas, are a modern group of lizards including about ~37 distinct species [83, but see 84]. They are typically grouped together with Asian Rock Agamas (*Laudakia*), *Bufoniceps*, and *Trapelus* in the subfamily Agaminae. Toad-headed agamas inhabit arid regions from northwestern China to the western side of the Caspian Sea, across the Tibetan Plateau, and southwest Asia to the Arabian Peninsula. They constitute one of the major components of the central Asian desert fauna and are highly adapted to sand dunes and desert environments. The reproductive biology of *Phrynocephalus* is notable in that there exist two reproductive modes: viviparity and oviparity. All six viviparous species are endemic to China and mainly restricted to high elevations in the Tibetan Plateau. Despite the progress that has been made in recent years for the systematics of certain *Phrynocephalus* groups, the large-scale pattern of their evolution in time and space remains open. Here, we review the origin and diversification of toad-headed agamas with their relationships to geological and climatic changes.

Until recently, our knowledge of the phylogenetic relationships and historical biogeography of *Phrynocephalus* was mainly based on anatomical studies [85-86]. Researchers usually agree that the broad distribution of *Phrynocephalus* may be explained by different hypotheses. Ananjeva and Tuniyev [87] speculated that there were two original centers for the *Phrynoce‐ phalus* species: Central Asia (a northern Tethys origin) and Middle Asia (a southern Tethys origin). They also inferred that diversification was related to geological events such as alpine orogenesis causing the isolation of valleys and basins and changes in the direction of river courses. Through a phylogenetic analysis of allozymes, Macey *et al.* [88] suggested that *Phrynocephalus* represents an old radiation that has been evolving in response to the Indian collision with Eurasia 35 Ma, and a clade in the former Soviet Union may have diverged as a result of internal Eurasian block movements in Afghanistan caused by the indenting Indian continent. Later, Arnold [86] proposed that *Phrynocephalus* originated in the Arabia-NW India area rather than in Central Asia, and achieved its current distribution by dispersal. However, Wang and Macey [89] initially considered the origin of the viviparous species group as a result of a vicariance event associated with the uplifting of the Tibetan Plateau, which was subse‐ quently supported by Zeng *et al.* [90] and Pang *et al.* [91]. Meanwhile, Pang *et al.* [91] using mitochondrial DNA data, partial sequences of the mitochondrial genes (12S rRNA, 16S rRNA, cyt *b*, ND4-tRNALeu), corroborated the monophyly of the viviparous group for the first time, and hypothesized that the oviparous group achieved their current distribution by dispersal. Pang *et al.* [91]. also provided the first time estimate for the origin of Chinese *Phrynocephalus* at 6.2–3.1 Ma, using average evolutionary rates of rate-constant genes to apply a clock.

both the analysis of Guo and Wang [92] and RAG-1 analysis did not include the basal species *P. interscapularis*, which was included in the mtDNA analysis by Melville *et al.* [93]. From this analysis, Melville *et al.* [93] concluded that including *P. interscapularis*, the age of *Phrynocepha‐ lus* would probably be Late Oligocene or Early Miocene, which coincides with the onset of aridification in Central Asia. Thus, their findings are consistent with a scenario that *Phrynce‐ phalus* is an old radiation, rather than it being a more recent diversification associated with the global deterioration of the Cenozoic climates when aridification intensified in the Late Miocene and into the Plio-Pleistocene. In addition, the deep phylogenetic split between the Chinese and Central Asian lineages of the *P. alpherakii*/*P. guttatus*/*P. versicolor* species group is probably related to the uplift of the Altay and Tien-Shan Mountains [93]. It was also suggested that gradual intensification of orogenetic processes and caused by them progressive aridization in the Middle East and Central Asia served as one of the factors favoring the evolution and initial

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

185

http://dx.doi.org/10.5772/55719

Dunayev *et al.* [95] pursued the evolutionary history of 35 *Prhynocepahlus* species from COI segments (600-750 bp), which is, so far, the phylogenetic study with the most comprehensive taxonomic coverage. This phylogeny shows that the clade of Iranian species (e.g. *P. scutella‐ tus* and *P. maculatus*) is sister to the remaining taxa, which joins Middle-Asian and Central-Asian taxa. Within the latter clade the representatives of the subgenus *Megalochilus* form a monophyletic group (*P. mystaceus*, *P. interscapularis*, *P. sogdianus*) of probably Irano-Turanian origin. The monophyletic Tibeto-Himalayan subgenus *Oreosaura* joins viviparous species (*P. vlangalii*, *P. thebaldi*, *P. forsythia* etc.) and is grouped together with oviparous species of Central Asia. Within Middle-Asian taxa three main species groups can be distinguished: (i) the groups of sun-watcher agamas (*P. helioscopus*, *P. ersicus*), (ii) *P. ocellatus-P. strauchi* group and (iii) *P. guttatus-P. versicolor* group. They further inferred that various bigoegoraphic events might cause diversification of *Phrynocepahlus* in continental Asia: starting with orogenetic processes on the territory of present-day Iranian and Tibetan plateaus, transgressions of Paratethys as major factors for old-splits, and up to aridisation processes and dynamics of major river valleys

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 *et al.* [96] studied the phylogeographic patterns of the Qinghai toad-headed lizard *Phrynocephalus vlangalii* by analysing sequence data from Mitochondrial DNA (mtDNA) sequences (partial ND2, tRNATrp and partial tRNAAla). The Qinghai toad-headed lizard is viviparous and lives at elevations of 2000–4500 m, and is the dominant terrestrial species in the QTP [26]. It is primarily distributed in the Qaidam Basin and surrounding area, with three major mountain chains, the Kunlun, the Arjin and the A'nyemaqen Mountains, in its range. With a thorough sampling across the entire distribution and a traditional phylogenetic analysis, five deeply diverged mtDNA lineages were recovered and geographical distribution of these lineages had almost no overlap [96]. Not only the spatial locations of the mountain chains coincided with the divisions of the lineages, but the times of the mountain uplift were also concordant with the times of mtDNA lineage divergence. Their NCPA analysis further

diversification of the species complex of *P. helioscopus* [94].

as factors determining recent speciation processes.

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 *et al.* [91] by employing partition-specific modeling in a combined DNA analysis to clarify existing gaps in the phylogeny of Chinese *Phrynocephalus*. Using this phylogenetic framework, they inferred the genus' historical biogeography by using weighted ancestral-area analysis and dispersalvicariance analysis in combination with a Bayesian relaxed molecular-clock approach. They drew three major conclusions: (i) the uplift of the QTP played a fundamental role in the diversification of viviparous *Phrynocepalus*; (ii) an evolutionary scenario combining aspects of vicariance and dispersal is necessary to explain the distribution of *Phrynocephalus*; (iii) Chinese *Phrynocephalus* originated at the Middle-Late Miocene boundary (12.78 Ma, 95% CI 17.61–8.25 Ma), and diversified from Late Miocene to Pleistocene from a center of origin in Central Asia, Tarim Basin, and Juggar Basin temperate desert.

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 *Phrynocephalus* diversification and evolution; (ii) intensification of aridity and geologic events in the Plio-Pleistocene and Quaternary glacial cycling probably had a significant influence on intraspecific diversification patterns within *Phrynocephalus*. Melville *et al.*[93] generated age estimates using two data. One method relied on mtDNA data, while the other relied on nuclear DNA data. The mtDNA data suggests that the estimated age of the common ancestor for *Phrynocephalus* was at 28.9 million years ago (Ma) with a 95% credibility interval of 36.2-21.1 Ma, a finding that is consistent with previous estimates [88]. Using the nuclear gene RAG-1excluding *P. interscapularis*, the estimate for the age of the common ancestor was somewhat younger, 15.8 Ma (95% CI 23–11.8 Ma), in accord with estimates produced by Guo and Wang [92]. However, both the analysis of Guo and Wang [92] and RAG-1 analysis did not include the basal species *P. interscapularis*, which was included in the mtDNA analysis by Melville *et al.* [93]. From this analysis, Melville *et al.* [93] concluded that including *P. interscapularis*, the age of *Phrynocepha‐ lus* would probably be Late Oligocene or Early Miocene, which coincides with the onset of aridification in Central Asia. Thus, their findings are consistent with a scenario that *Phrynce‐ phalus* is an old radiation, rather than it being a more recent diversification associated with the global deterioration of the Cenozoic climates when aridification intensified in the Late Miocene and into the Plio-Pleistocene. In addition, the deep phylogenetic split between the Chinese and Central Asian lineages of the *P. alpherakii*/*P. guttatus*/*P. versicolor* species group is probably related to the uplift of the Altay and Tien-Shan Mountains [93]. It was also suggested that gradual intensification of orogenetic processes and caused by them progressive aridization in the Middle East and Central Asia served as one of the factors favoring the evolution and initial diversification of the species complex of *P. helioscopus* [94].

origin). They also inferred that diversification was related to geological events such as alpine orogenesis causing the isolation of valleys and basins and changes in the direction of river courses. Through a phylogenetic analysis of allozymes, Macey *et al.* [88] suggested that *Phrynocephalus* represents an old radiation that has been evolving in response to the Indian collision with Eurasia 35 Ma, and a clade in the former Soviet Union may have diverged as a result of internal Eurasian block movements in Afghanistan caused by the indenting Indian continent. Later, Arnold [86] proposed that *Phrynocephalus* originated in the Arabia-NW India area rather than in Central Asia, and achieved its current distribution by dispersal. However, Wang and Macey [89] initially considered the origin of the viviparous species group as a result of a vicariance event associated with the uplifting of the Tibetan Plateau, which was subse‐ quently supported by Zeng *et al.* [90] and Pang *et al.* [91]. Meanwhile, Pang *et al.* [91] using mitochondrial DNA data, partial sequences of the mitochondrial genes (12S rRNA, 16S rRNA, cyt *b*, ND4-tRNALeu), corroborated the monophyly of the viviparous group for the first time, and hypothesized that the oviparous group achieved their current distribution by dispersal. Pang *et al.* [91]. also provided the first time estimate for the origin of Chinese *Phrynocephalus* at 6.2–3.1 Ma, using average evolutionary rates of rate-constant genes to apply a clock.

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 *et al.* [91] by employing partition-specific modeling in a combined DNA analysis to clarify existing gaps in the phylogeny of Chinese *Phrynocephalus*. Using this phylogenetic framework, they inferred the genus' historical biogeography by using weighted ancestral-area analysis and dispersalvicariance analysis in combination with a Bayesian relaxed molecular-clock approach. They drew three major conclusions: (i) the uplift of the QTP played a fundamental role in the diversification of viviparous *Phrynocepalus*; (ii) an evolutionary scenario combining aspects of vicariance and dispersal is necessary to explain the distribution of *Phrynocephalus*; (iii) Chinese *Phrynocephalus* originated at the Middle-Late Miocene boundary (12.78 Ma, 95% CI 17.61–8.25 Ma), and diversified from Late Miocene to Pleistocene from a center of origin in Central Asia,

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 *Phrynocephalus* diversification and evolution; (ii) intensification of aridity and geologic events in the Plio-Pleistocene and Quaternary glacial cycling probably had a significant influence on intraspecific diversification patterns within *Phrynocephalus*. Melville *et al.*[93] generated age estimates using two data. One method relied on mtDNA data, while the other relied on nuclear DNA data. The mtDNA data suggests that the estimated age of the common ancestor for *Phrynocephalus* was at 28.9 million years ago (Ma) with a 95% credibility interval of 36.2-21.1 Ma, a finding that is consistent with previous estimates [88]. Using the nuclear gene RAG-1excluding *P. interscapularis*, the estimate for the age of the common ancestor was somewhat younger, 15.8 Ma (95% CI 23–11.8 Ma), in accord with estimates produced by Guo and Wang [92]. However,

Tarim Basin, and Juggar Basin temperate desert.

184 Current Progress in Biological Research

Dunayev *et al.* [95] pursued the evolutionary history of 35 *Prhynocepahlus* species from COI segments (600-750 bp), which is, so far, the phylogenetic study with the most comprehensive taxonomic coverage. This phylogeny shows that the clade of Iranian species (e.g. *P. scutella‐ tus* and *P. maculatus*) is sister to the remaining taxa, which joins Middle-Asian and Central-Asian taxa. Within the latter clade the representatives of the subgenus *Megalochilus* form a monophyletic group (*P. mystaceus*, *P. interscapularis*, *P. sogdianus*) of probably Irano-Turanian origin. The monophyletic Tibeto-Himalayan subgenus *Oreosaura* joins viviparous species (*P. vlangalii*, *P. thebaldi*, *P. forsythia* etc.) and is grouped together with oviparous species of Central Asia. Within Middle-Asian taxa three main species groups can be distinguished: (i) the groups of sun-watcher agamas (*P. helioscopus*, *P. ersicus*), (ii) *P. ocellatus-P. strauchi* group and (iii) *P. guttatus-P. versicolor* group. They further inferred that various bigoegoraphic events might cause diversification of *Phrynocepahlus* in continental Asia: starting with orogenetic processes on the territory of present-day Iranian and Tibetan plateaus, transgressions of Paratethys as major factors for old-splits, and up to aridisation processes and dynamics of major river valleys as factors determining recent speciation processes.

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 *et al.* [96] studied the phylogeographic patterns of the Qinghai toad-headed lizard *Phrynocephalus vlangalii* by analysing sequence data from Mitochondrial DNA (mtDNA) sequences (partial ND2, tRNATrp and partial tRNAAla). The Qinghai toad-headed lizard is viviparous and lives at elevations of 2000–4500 m, and is the dominant terrestrial species in the QTP [26]. It is primarily distributed in the Qaidam Basin and surrounding area, with three major mountain chains, the Kunlun, the Arjin and the A'nyemaqen Mountains, in its range. With a thorough sampling across the entire distribution and a traditional phylogenetic analysis, five deeply diverged mtDNA lineages were recovered and geographical distribution of these lineages had almost no overlap [96]. Not only the spatial locations of the mountain chains coincided with the divisions of the lineages, but the times of the mountain uplift were also concordant with the times of mtDNA lineage divergence. Their NCPA analysis further demonstrated that the allopatric fragmentation resulted from historical vicariant events dominated the history of this species. Consequently, Jin *et al.* [96] concluded that the uplifts of the three major mountain chains formed the physical barrier that caused the initial vicariant events. They also inferred that populations from the Qaidam Basin appeared to have under‐ gone major demographic and range expansions in the early Pleistocene (1.6 Ma), consistent with the colonization of areas previously covered by the huge Qaidam palaeolake, which desiccated at the onset of the Pleistocene [96].

**4.2. Racerunner lizards (***Eremias***)**

changes.

The lacertid genus *Eremias* Fitzinger 1834, is considered to comprise approximately 34 species, which inhabit sand, steppe, and desert regions from northern China, Mongolia, Korea, Central and Southwest Asia to Southeastern Europe. The reproductive biology of *Eremias* is notable in that there exist two reproductive modes: viviparity and oviparity. Most are oviparous, whereas the *Eremias multiocellata* complex (comprising 6 subspecies), *E. buechneri*, *E. kokshaaliensis*, *E. yarkandensis*, *E. quadrifrons*, and *E. przewalskii* are viviparous [102]. Despite the progress that has been made in recent years for the systematics of certain *Eremias* groups, the large-scale pattern of their evolution in time and space remains open. Here, we review the origin and diversification of racerunner lizards with their relationships to geological and climatic

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

187

http://dx.doi.org/10.5772/55719

Until recently, our understanding of the phylogenetic relationships and historical biogeogra‐ phy of *Eremias* was mainly based on anatomical studies [103-105]. Boulenger [106] and FitzSimons [107] assigned most of the species now placed in *Pedioplanis* to the subgenus *Mesalina* within the large genus *Eremias*. Szczerbak [108] regarded *Eremias* sensu lato poly‐ phyletic and considered *Eremias* (s. s.) endemic to Asia. Based on morphological characters and geographic distribution, Szczerbak [103] subdivided the inclusive genus *Eremias*(s. l.) into two distinct genera: the genus *Mesalina* as a north African and lowland Southwest Asian clade, and the genus *Eremias* (s. s.), which is endemic to Asia. Furthermore, Szczerbak [103] subdi‐ vided *Eremias* into five distinct subgenera: *Eremias* Fitzinger in Wiegmann, 1834 (group E. velox), *Rhabderemias* Lantz, 1928 (group *Eremias scripta*–*Eremias lineolata*), *Ommateremias* Lantz, 1928 (group *E. arguta*), *Scapteira* Fitzinger in Wiegmann, 1834 (group *Eremias grammica*), and *Pareremias* Szczerbak, 1973 (group *E. multiocellata*). The five subgenera were supported by

Wan *et al.* [109] made the first attempt to elucidate the phylogenetic relationships among the Chinese racerunner lizards on the basis of mitochondrial 16S rRNA data. They found that *E. brenchleyi* and *E. argus* formed a clade as the sister group of *E. multiocellata*. They further inferred that *E. arguta*, *E.grammica* and *E. velox* originated from Central Asia and the rest species in China originated from East Asia. However, these conclusions were tentative due to both limited taxa sampling as well as use of relatively short, partly nondiagnostic, gene fragment. On the other hand, one potential shortcoming of the analyses of Wan *et al.* [109] was that they

Orlova *et al.* [110] studied the phylogeographic pattern of the steppe-runner lizard *Eremias arguta* by analyzing mtDNA cyt *b* gene sequences. The steppe-runner lizard inhabits steppes and semi-deserts of Eastern Europe and Middle Asia from Romania to Western Mongolia and China. They found an old split between *E. arguta uzbekistanica* and all other taxa, probably coming from the area of Ustyurt plateau. Consequently, Orlova *et al.* [110] inferred that this

More recently, a study [111 using combined DNA data sets (3925 bp) from two mitochondrial genes (cyt *b*, 12S rRNA) and one nuclear gene (RAG-1) and comprehensive taxonomic coverage resulted in three well-supported conclusions: (i) the species of the traditional genus *Eremias*

did not incorporate secondary structural constraints of 16S rRNA into analyses.

vicariant event is likely to have been caused by Paratethys regression.

Arnold (1986) on the basis of the hemipenial characters.

Using data from 11 microsatellite DNA loci, Wang *et al.* [97] further supported the hypotheis of Jin *et al.* [96] that the uplift of the Arjin and the A'nyemaqen Mountains caused the initial vicariant events that led to the formation of the diverged lineages within *P. vlangalii*. Later, using a mitochondrial fragment ND4-tRNALeu, Guo *et al.* [98]confirmed that the uplift of the A'nyemaqen Mountains and glaciations since the mid-late Pleistocene, especially during the Kunlun Glaciation, are considered to have promoted the allopartric divergence of *P. vlangalii*. The diversification of *P. putjatia* may be triggered by the tectonic movement in the Huangshui River valley during the C phase of Qingzang Movement. Subsequently, the glacial climate throughout the Pleistocene may have continued to impede the gene flow of *P. putjatia*, eventually resulting in the genetic divergence of *P. putjatia* in the allopatric regions.

Jin and Liu [99] described the phylogeography of a unique endemic agamid lizard, *Phrynoce‐ phalus erythruru*s from the Qiangtang Plateau by analyzing sequence data from a mtDNA segment ND2- tRNAAla. *P. erythrurus* is viviparous and occupies the highest regions of any reptile on earth [26]. They found that (i) this species diverged into two major lineages/ subspecies at 3.7 Ma corresponding to the Northern and Southern Qiangtang Plateau; (ii) the Northern Qiangtang lineage diverged into two subpopulations at 2.8 Ma separated by the Beilu River Region and Wulanwula Mountains [99]. Their NCPA analysis further demonstrated that the allopatric fragmentation and restricted gene flow were the most likely mechanisms of population differentiation. Their results also indicated the presence of at least three refugia since the Hongya glaciation. Consequently, Jin and Liu [99] concluded that the uplift of Tanggula Mountains movement and glaciations since mid-Pliocene have shaped phylogenetic patterns of *P. erythrurus*.

An aridification of the Tarim Basin and adjacent areas since middle Pleistocene has produced significant genetic structuring of the local fauna. Zhang *et al.* [100] compared the phylogeo‐ graphic patterns, population structure and history of *Phrynocephalus axillaris* and *P. forsythii* using a mitochondrial fragment ND4-tRNALeu. They demonstrated that the two species might have experienced different evolutionary history throughout their current distribution. For *P. forsythii*, a vicariant event, as a consequence of geological isolation by the initiation of Quater‐ nary folding at 2.1–1.2 Ma in the eastern Tien-Shan [101], and desert expansion, might have produced the significant divergence between the Tarim and the Yanqi populations. For *P. axillaris*, populations of the Yanqi, Turpan and Hami Basins might have been established through dispersal during demographic expansion. Climatic fluctuations caused alternate expansion and shrinkage of rivers and oases several times, which likely led to habitat frag‐ mentation for both species. Interaction between vicariance, dispersal and habitat fragmenta‐ tion produced the current distribution and genetic diversity.

#### **4.2. Racerunner lizards (***Eremias***)**

demonstrated that the allopatric fragmentation resulted from historical vicariant events dominated the history of this species. Consequently, Jin *et al.* [96] concluded that the uplifts of the three major mountain chains formed the physical barrier that caused the initial vicariant events. They also inferred that populations from the Qaidam Basin appeared to have under‐ gone major demographic and range expansions in the early Pleistocene (1.6 Ma), consistent with the colonization of areas previously covered by the huge Qaidam palaeolake, which

Using data from 11 microsatellite DNA loci, Wang *et al.* [97] further supported the hypotheis of Jin *et al.* [96] that the uplift of the Arjin and the A'nyemaqen Mountains caused the initial vicariant events that led to the formation of the diverged lineages within *P. vlangalii*. Later, using a mitochondrial fragment ND4-tRNALeu, Guo *et al.* [98]confirmed that the uplift of the A'nyemaqen Mountains and glaciations since the mid-late Pleistocene, especially during the Kunlun Glaciation, are considered to have promoted the allopartric divergence of *P. vlangalii*. The diversification of *P. putjatia* may be triggered by the tectonic movement in the Huangshui River valley during the C phase of Qingzang Movement. Subsequently, the glacial climate throughout the Pleistocene may have continued to impede the gene flow of *P. putjatia*,

eventually resulting in the genetic divergence of *P. putjatia* in the allopatric regions.

Jin and Liu [99] described the phylogeography of a unique endemic agamid lizard, *Phrynoce‐ phalus erythruru*s from the Qiangtang Plateau by analyzing sequence data from a mtDNA segment ND2- tRNAAla. *P. erythrurus* is viviparous and occupies the highest regions of any reptile on earth [26]. They found that (i) this species diverged into two major lineages/ subspecies at 3.7 Ma corresponding to the Northern and Southern Qiangtang Plateau; (ii) the Northern Qiangtang lineage diverged into two subpopulations at 2.8 Ma separated by the Beilu River Region and Wulanwula Mountains [99]. Their NCPA analysis further demonstrated that the allopatric fragmentation and restricted gene flow were the most likely mechanisms of population differentiation. Their results also indicated the presence of at least three refugia since the Hongya glaciation. Consequently, Jin and Liu [99] concluded that the uplift of Tanggula Mountains movement and glaciations since mid-Pliocene have shaped phylogenetic

An aridification of the Tarim Basin and adjacent areas since middle Pleistocene has produced significant genetic structuring of the local fauna. Zhang *et al.* [100] compared the phylogeo‐ graphic patterns, population structure and history of *Phrynocephalus axillaris* and *P. forsythii* using a mitochondrial fragment ND4-tRNALeu. They demonstrated that the two species might have experienced different evolutionary history throughout their current distribution. For *P. forsythii*, a vicariant event, as a consequence of geological isolation by the initiation of Quater‐ nary folding at 2.1–1.2 Ma in the eastern Tien-Shan [101], and desert expansion, might have produced the significant divergence between the Tarim and the Yanqi populations. For *P. axillaris*, populations of the Yanqi, Turpan and Hami Basins might have been established through dispersal during demographic expansion. Climatic fluctuations caused alternate expansion and shrinkage of rivers and oases several times, which likely led to habitat frag‐ mentation for both species. Interaction between vicariance, dispersal and habitat fragmenta‐

tion produced the current distribution and genetic diversity.

desiccated at the onset of the Pleistocene [96].

186 Current Progress in Biological Research

patterns of *P. erythrurus*.

The lacertid genus *Eremias* Fitzinger 1834, is considered to comprise approximately 34 species, which inhabit sand, steppe, and desert regions from northern China, Mongolia, Korea, Central and Southwest Asia to Southeastern Europe. The reproductive biology of *Eremias* is notable in that there exist two reproductive modes: viviparity and oviparity. Most are oviparous, whereas the *Eremias multiocellata* complex (comprising 6 subspecies), *E. buechneri*, *E. kokshaaliensis*, *E. yarkandensis*, *E. quadrifrons*, and *E. przewalskii* are viviparous [102]. Despite the progress that has been made in recent years for the systematics of certain *Eremias* groups, the large-scale pattern of their evolution in time and space remains open. Here, we review the origin and diversification of racerunner lizards with their relationships to geological and climatic changes.

Until recently, our understanding of the phylogenetic relationships and historical biogeogra‐ phy of *Eremias* was mainly based on anatomical studies [103-105]. Boulenger [106] and FitzSimons [107] assigned most of the species now placed in *Pedioplanis* to the subgenus *Mesalina* within the large genus *Eremias*. Szczerbak [108] regarded *Eremias* sensu lato poly‐ phyletic and considered *Eremias* (s. s.) endemic to Asia. Based on morphological characters and geographic distribution, Szczerbak [103] subdivided the inclusive genus *Eremias*(s. l.) into two distinct genera: the genus *Mesalina* as a north African and lowland Southwest Asian clade, and the genus *Eremias* (s. s.), which is endemic to Asia. Furthermore, Szczerbak [103] subdi‐ vided *Eremias* into five distinct subgenera: *Eremias* Fitzinger in Wiegmann, 1834 (group E. velox), *Rhabderemias* Lantz, 1928 (group *Eremias scripta*–*Eremias lineolata*), *Ommateremias* Lantz, 1928 (group *E. arguta*), *Scapteira* Fitzinger in Wiegmann, 1834 (group *Eremias grammica*), and *Pareremias* Szczerbak, 1973 (group *E. multiocellata*). The five subgenera were supported by Arnold (1986) on the basis of the hemipenial characters.

Wan *et al.* [109] made the first attempt to elucidate the phylogenetic relationships among the Chinese racerunner lizards on the basis of mitochondrial 16S rRNA data. They found that *E. brenchleyi* and *E. argus* formed a clade as the sister group of *E. multiocellata*. They further inferred that *E. arguta*, *E.grammica* and *E. velox* originated from Central Asia and the rest species in China originated from East Asia. However, these conclusions were tentative due to both limited taxa sampling as well as use of relatively short, partly nondiagnostic, gene fragment. On the other hand, one potential shortcoming of the analyses of Wan *et al.* [109] was that they did not incorporate secondary structural constraints of 16S rRNA into analyses.

Orlova *et al.* [110] studied the phylogeographic pattern of the steppe-runner lizard *Eremias arguta* by analyzing mtDNA cyt *b* gene sequences. The steppe-runner lizard inhabits steppes and semi-deserts of Eastern Europe and Middle Asia from Romania to Western Mongolia and China. They found an old split between *E. arguta uzbekistanica* and all other taxa, probably coming from the area of Ustyurt plateau. Consequently, Orlova *et al.* [110] inferred that this vicariant event is likely to have been caused by Paratethys regression.

More recently, a study [111 using combined DNA data sets (3925 bp) from two mitochondrial genes (cyt *b*, 12S rRNA) and one nuclear gene (RAG-1) and comprehensive taxonomic coverage resulted in three well-supported conclusions: (i) the species of the traditional genus *Eremias* form six clades; (ii) the Iranian plateau is the center of origin of the genus as a whole; (iii) genus *Eremias* should be divided into 4 distinct genera, and at least 14 new species should be described within these genera. Using a molecular clock calibrated with well dated paleological events as calibration, the average evolutionary rate for the mitochondrial genes was estimated as 1.3%-1.6% per million years. Based on the mitochondrial DNA data and using the globe molecular clock method, the origin time for the genus *Eremias* was dated to the early Miocene and the start of intrageneric divergence to the middle Miocene.

**Appendix S1**

**Lacertilia:**

List of the 422 squamates species.

*salvator, Zootoca vivipara.*

**Serpentes**

*Acanthosaura armata, A. lepidogaster, Alsophylax pipiens, A. przewalskii, Asymblepharus alaicus, A. ladacensis, A. sikimmensis, Ateuchosaurus chinensis, Calotes emma, C. medogensis, C. mystaceus, C. versicolor, Cyrtodactylus khasiensis, C. zhaoermii, Cyrtopodion dadunense, C. elongatum, C. medogense, C. stoliczkai, C. tibetanus, Dibamus bogadeki, D. bourreti, Dopasia gracilis, D. hainanensis, D. harti, D. ludovici, Draco blanfordii, D. maculatus, Eremias argus, E. arguta, E. brenchleyi, E. buechneri, E. grammica, E. kokshaaliensis, E. multiocellata, E. przewal‐ skii, E. quadrifrons, E. stummeri, E. velox, E. vermiculata, E. yarkandensis, Eutropis longicauda‐ ta, E. multicarinata, E. multifasciata, Gehyra mutilate, Gekko auriverrucosus, G. chinensis, G. gecko, G. hokouensis, G. japonicas, G. kikuchii, G. melli, G. palmatus, G. reevesii, G. scabridus, G. similignum, G. swinhonis, G. taibaiensis, G. wenxianensis, Goniurosaurus bawanglingensis, G. hainanensis, G. lichtenfelderi, G. luii, G. yingdeensis, Hemidactylus aquilonius, H. bowringii, H. frenatus, H. garnotii, H. platyurus, H. stejnegeri, Hemiphyllodactylus typus, H. yunnanensis, Japalura andersoniana, J. batangensis, J. brevipes, J. brevicauda, J. dymondi, J. fasciata, J. flavi‐ ceps, J. graham, J. luei, J. kumaonensis, J. makii, J. micangshanensis, J. polygonata, J. splendida, J. swinhonis, J. tricarinata, J. varcoae, J. yulongensis, J. yunnanensis, J. zhaoermii, Lacerta agilis, Lamprolepis smaragdina, Laudakia badakhshana, L. himalayana, L. papenfussi, L. sacra, L. stoliczkana, L. tuberculata, L. wui, Leiolepis guttata, L. reevesii, Lepidodactylus lugubris, L. yami, Lygosoma bowringii, L. quadrupes, Mediodactylus russowii, Oriocalotes Paulus, Phrynocephalus axillaris, P. forsythia, P. guinanensis, P. guttatus, P. helioscopus, P. mystaceus, P. przewalskii, P. putjatai, P. roborowskii, P. theobaldi, P. versicolor, P. vlangalii, Physignathus cocincinus, Plestiodon capito, P. chinensis, P. elegans, P. liui, P. marginatus, P. popei, P. quadrilineatus, P. tamdaoen‐ sis, P. tunganus, Pseudocalotes brevipes, P. kakhienensis, P. kingdonwardi, P. microlepis, Ptyctolae‐ mus gularis, Scincella barbouri, S. doriae, S. formosensis, S. huanrenensis, S. modesta, S. monticola, S. potanini, S. przewalskii, S. reevesii, S. schmidti, S. tsinlingensis, Shinisaurus crocodilurus, Sphenomorphus courcyanum, S. incognitus, S. indicus, S. maculatus, S. taiwanensis, Takydro‐ mus amurensis, T. formosanus, T. hsuehshanensis, T. intermedius, T. kuehnei, T. luyeanus, T. sauteri, T. septentrionalis, T. stejnegeri, T. sexlineatus, T. sylvaticus, T. viridipunctatus, T. wolteri, Teratoscincus przewalskii, T. roborowskii, T. scincus, T. toksunicus, Trapelus agilis, T. sanguinolen‐ tus, Tropidophorus berdmorei, T. guangxiensis, T. hainanus, T. sinicus, Varanus bengalensis, V.*

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

189

http://dx.doi.org/10.5772/55719

*Acalyptophis peronei, Achalinus ater, A. formosanus, A. hainanus, A. jinggangensis, A. meiguensis, A. niger, A. rufescens, A. spinalis, Acrochordus granulatus, Ahaetulla prasina, Aipysurus eydouxii, Amphiesma atemporale, A. bitaeniatum, A. boulengeri, A. craspedogaster, A. johannis, A. khasiense, A. metusia, A. miyajimae, A. modestum, A. octolineatum, A. optatum, A. parallelum, A. platyceps, A. popei, A. sauteri, A. stolatu, A. venningi, A. vibakari, Amphiesmoides ornaticeps, Archelaphe bella, Atretium yunnanensis, Azemiops feae, Blythia reticulata, Boiga cyanea, B. guangxiensis, B. kraepelini,*

Another phylogenetic study [112] using mitochondrial 16S rRNA segment results in five major conclusions: (i) monophyly of *Eremias* and a clade comprising *Eremias*, *Acanthodactylus* and *Latastia* are recognized; (ii) monophyly of the subgenus *Pareremias* is corroborated, with *Eremias argus* being the sister taxon to *Eremias brenchleyi*; (iii) the first evidence that viviparous species form a monophyletic group was presented; (iv) *Eremias* probably diversified at about 9.9 million years ago (with the 95% credibility interval ranging from 7.6 to 12 Ma); (v) specifically, the divergence time of the subgenus *Pareremias* was dated to about 6.3 million years ago (with the 95% credibility interval ranging from 5.3 to 8.5 Ma), suggesting that the diversification of this subgenus might be correlated with the evolution of an East Asian monsoon climate triggered by the rapid uplift of the Tibetan Plateau approximately 8 Ma. A recent phylogeo‐ graphic study [113] indicated that the Yellow River and Taihang Mountains may have acted as important barriers to gene flow in *Eremias brenchleyi*.

Rastegar-Pouyani *et al.* [114] explored the genealogical relationships and intraspecific differ‐ entiation of the rapid fringe-toed lizard *Eremias velox* by analyzing mtDNA cyt *b* and 12S rRNA gene sequences. The rapid fringe-toed lizard is widely distributed in the Iranian Plateau and Central Asia. In combination with molecular dating, they inferred that the *E. velox* complex originated on the Iranian Plateau in the Middle Miocene, which is in accordant with the result of Rastegar-Pouyani [115]. Another founding was that the northern Iranian clade diverged first some 11–10 Ma, caused by first uplifting of the Elburz Mountains in the late Miocene, and that the Central Asian lineages split from the northeastern Iranian lineage approximately 6 Ma, most likely as a result of uplifting of the Kopet-Dagh Mountains in the northern margin of the Iranian Plateau.

#### **4.3. Conclusion**

As exemplified with toad-headed agamas (genus *Phrynocephalus*) and racerunner lizards (genus *Eremias*), geomorphic and climatic changes in this big area definitely have remark‐ able influences on the regional and adjacent biogeographic patterns. Future comparative studies between *Eremias, Phrynocephalus* and other groups that distribute at the same area and inhabit similar habitat, can help to elucidate the origination, diversification and dispersion of these lizards, and test the link between speciation, adaptation, and histori‐ cal changes in their biogeography. Specifically, further studies could search for differen‐ ces in the ecological niches of species and clusters in order to assess ecological divergence among the several groups [116].

### **Appendix S1**

List of the 422 squamates species.

#### **Lacertilia:**

form six clades; (ii) the Iranian plateau is the center of origin of the genus as a whole; (iii) genus *Eremias* should be divided into 4 distinct genera, and at least 14 new species should be described within these genera. Using a molecular clock calibrated with well dated paleological events as calibration, the average evolutionary rate for the mitochondrial genes was estimated as 1.3%-1.6% per million years. Based on the mitochondrial DNA data and using the globe molecular clock method, the origin time for the genus *Eremias* was dated to the early Miocene

Another phylogenetic study [112] using mitochondrial 16S rRNA segment results in five major conclusions: (i) monophyly of *Eremias* and a clade comprising *Eremias*, *Acanthodactylus* and *Latastia* are recognized; (ii) monophyly of the subgenus *Pareremias* is corroborated, with *Eremias argus* being the sister taxon to *Eremias brenchleyi*; (iii) the first evidence that viviparous species form a monophyletic group was presented; (iv) *Eremias* probably diversified at about 9.9 million years ago (with the 95% credibility interval ranging from 7.6 to 12 Ma); (v) specifically, the divergence time of the subgenus *Pareremias* was dated to about 6.3 million years ago (with the 95% credibility interval ranging from 5.3 to 8.5 Ma), suggesting that the diversification of this subgenus might be correlated with the evolution of an East Asian monsoon climate triggered by the rapid uplift of the Tibetan Plateau approximately 8 Ma. A recent phylogeo‐ graphic study [113] indicated that the Yellow River and Taihang Mountains may have acted

Rastegar-Pouyani *et al.* [114] explored the genealogical relationships and intraspecific differ‐ entiation of the rapid fringe-toed lizard *Eremias velox* by analyzing mtDNA cyt *b* and 12S rRNA gene sequences. The rapid fringe-toed lizard is widely distributed in the Iranian Plateau and Central Asia. In combination with molecular dating, they inferred that the *E. velox* complex originated on the Iranian Plateau in the Middle Miocene, which is in accordant with the result of Rastegar-Pouyani [115]. Another founding was that the northern Iranian clade diverged first some 11–10 Ma, caused by first uplifting of the Elburz Mountains in the late Miocene, and that the Central Asian lineages split from the northeastern Iranian lineage approximately 6 Ma, most likely as a result of uplifting of the Kopet-Dagh Mountains in the northern margin

As exemplified with toad-headed agamas (genus *Phrynocephalus*) and racerunner lizards (genus *Eremias*), geomorphic and climatic changes in this big area definitely have remark‐ able influences on the regional and adjacent biogeographic patterns. Future comparative studies between *Eremias, Phrynocephalus* and other groups that distribute at the same area and inhabit similar habitat, can help to elucidate the origination, diversification and dispersion of these lizards, and test the link between speciation, adaptation, and histori‐ cal changes in their biogeography. Specifically, further studies could search for differen‐ ces in the ecological niches of species and clusters in order to assess ecological divergence

and the start of intrageneric divergence to the middle Miocene.

as important barriers to gene flow in *Eremias brenchleyi*.

of the Iranian Plateau.

188 Current Progress in Biological Research

among the several groups [116].

**4.3. Conclusion**

*Acanthosaura armata, A. lepidogaster, Alsophylax pipiens, A. przewalskii, Asymblepharus alaicus, A. ladacensis, A. sikimmensis, Ateuchosaurus chinensis, Calotes emma, C. medogensis, C. mystaceus, C. versicolor, Cyrtodactylus khasiensis, C. zhaoermii, Cyrtopodion dadunense, C. elongatum, C. medogense, C. stoliczkai, C. tibetanus, Dibamus bogadeki, D. bourreti, Dopasia gracilis, D. hainanensis, D. harti, D. ludovici, Draco blanfordii, D. maculatus, Eremias argus, E. arguta, E. brenchleyi, E. buechneri, E. grammica, E. kokshaaliensis, E. multiocellata, E. przewal‐ skii, E. quadrifrons, E. stummeri, E. velox, E. vermiculata, E. yarkandensis, Eutropis longicauda‐ ta, E. multicarinata, E. multifasciata, Gehyra mutilate, Gekko auriverrucosus, G. chinensis, G. gecko, G. hokouensis, G. japonicas, G. kikuchii, G. melli, G. palmatus, G. reevesii, G. scabridus, G. similignum, G. swinhonis, G. taibaiensis, G. wenxianensis, Goniurosaurus bawanglingensis, G. hainanensis, G. lichtenfelderi, G. luii, G. yingdeensis, Hemidactylus aquilonius, H. bowringii, H. frenatus, H. garnotii, H. platyurus, H. stejnegeri, Hemiphyllodactylus typus, H. yunnanensis, Japalura andersoniana, J. batangensis, J. brevipes, J. brevicauda, J. dymondi, J. fasciata, J. flavi‐ ceps, J. graham, J. luei, J. kumaonensis, J. makii, J. micangshanensis, J. polygonata, J. splendida, J. swinhonis, J. tricarinata, J. varcoae, J. yulongensis, J. yunnanensis, J. zhaoermii, Lacerta agilis, Lamprolepis smaragdina, Laudakia badakhshana, L. himalayana, L. papenfussi, L. sacra, L. stoliczkana, L. tuberculata, L. wui, Leiolepis guttata, L. reevesii, Lepidodactylus lugubris, L. yami, Lygosoma bowringii, L. quadrupes, Mediodactylus russowii, Oriocalotes Paulus, Phrynocephalus axillaris, P. forsythia, P. guinanensis, P. guttatus, P. helioscopus, P. mystaceus, P. przewalskii, P. putjatai, P. roborowskii, P. theobaldi, P. versicolor, P. vlangalii, Physignathus cocincinus, Plestiodon capito, P. chinensis, P. elegans, P. liui, P. marginatus, P. popei, P. quadrilineatus, P. tamdaoen‐ sis, P. tunganus, Pseudocalotes brevipes, P. kakhienensis, P. kingdonwardi, P. microlepis, Ptyctolae‐ mus gularis, Scincella barbouri, S. doriae, S. formosensis, S. huanrenensis, S. modesta, S. monticola, S. potanini, S. przewalskii, S. reevesii, S. schmidti, S. tsinlingensis, Shinisaurus crocodilurus, Sphenomorphus courcyanum, S. incognitus, S. indicus, S. maculatus, S. taiwanensis, Takydro‐ mus amurensis, T. formosanus, T. hsuehshanensis, T. intermedius, T. kuehnei, T. luyeanus, T. sauteri, T. septentrionalis, T. stejnegeri, T. sexlineatus, T. sylvaticus, T. viridipunctatus, T. wolteri, Teratoscincus przewalskii, T. roborowskii, T. scincus, T. toksunicus, Trapelus agilis, T. sanguinolen‐ tus, Tropidophorus berdmorei, T. guangxiensis, T. hainanus, T. sinicus, Varanus bengalensis, V. salvator, Zootoca vivipara.*

#### **Serpentes**

*Acalyptophis peronei, Achalinus ater, A. formosanus, A. hainanus, A. jinggangensis, A. meiguensis, A. niger, A. rufescens, A. spinalis, Acrochordus granulatus, Ahaetulla prasina, Aipysurus eydouxii, Amphiesma atemporale, A. bitaeniatum, A. boulengeri, A. craspedogaster, A. johannis, A. khasiense, A. metusia, A. miyajimae, A. modestum, A. octolineatum, A. optatum, A. parallelum, A. platyceps, A. popei, A. sauteri, A. stolatu, A. venningi, A. vibakari, Amphiesmoides ornaticeps, Archelaphe bella, Atretium yunnanensis, Azemiops feae, Blythia reticulata, Boiga cyanea, B. guangxiensis, B. kraepelini,* *B. multomaculata, B. nigriceps, Bungarus bungaroides, B. fasciatus, B. multicinctus, Calamaria buchi, C. pavimentata, C. septentrionalis, C. yunnanensis, Calliophis maculiceps, Chitulia inornata, C. ornate, C. torquata, Chrysopelea ornate, Coelognathus radiates, Cyclophiops doriae, C. major, C. multicinctus, C. ruffus, Daboia siamensis, Deinagkistrodon acutus, Dendrelaphis biloreatus, D.hollinrakei, D. ngansonensis, D. pictus, D. subocularis, Dinodon flavozonatum, D. rosozonatum, D. rufozonatum, D. septentrionalis, Disteira stokesii, Elaphe anomala, E. bimaculata, E. carinata, E. davidi, E. dione, E. schrenckii, E. zoigeensis, Emydocephalus ijimae, Enhydris bennettii, E. bocourti, E. chinensis, E.en‐ hydris, E. plumbea, Eryx miliaris, E. tataricus, Euprepiophis mandarinus, E. perlacea, Gloydius blomhoffii, G. brevicaudus, G. halys, G. intermedius, G. lijianlii, G. monticola, G. saxatilis, G. shedaoensis, G. strauchi, G. ussuriensis, Hemorrhois ravergieri, Hierophis spinalis, Hydrophis fasciatus, H. gracilis, H. parviceps, Kerilia jerdonii, Lapemis curtus, L. hardwickii, L. colubrine, L. laticaudata, Leioselasma cyanocincta, L. melanocephala, L. spiralis, Liopeltis frenatus, Lycodon capucinus, L. fasciatus, L. futsingensis, L. gongshan, L. laoensis, L. liuchengchaoi, L. ruhstrati, L. subcinctus, L. synaptor, Macropisthodon rudis, Naja atra, N. kaouthia, Natrix natrix, N. tessellate, Oligodon albocinctus, O. catenatus, O. chinensis, O. cinereus, O. cyclurus, O. eberhardti, O. formo‐ sanus, O. lacroixi, O. lungshenensis, O. melanozonatus, O. multizonatus, O. ningshaanensis, O. ornatus, O. taeniatus, Oocatochus rufodorsatus, Ophiophagus Hannah, Opisthotropis andersonii, O. balteata, O. cheni, O. guangxiensis, O. jacobi, O. kuatunensis, O. lateralis, O. latouchii, O. maculosa, O. maxwelli, Oreocryptophis porphyraceus, Orthriophis hodgsoni, O. moellendorffi, O. taeniurus, Ovophis monticola, O. tonkinensis, O. zayuensis, Paratapinophis praemaxillaris, Pareas boulengeri, P. carinatus, P. chinensis, P. formosensis, P. hamptoni, P. margaritophorus, P. monticola, P. nigriceps, P. stanleyi, Pelamis platura, Plagiopholis blakewayi, P. nuchalis, P. styani, P. unipostocularis, Polyodontognathus caerulescens, Praescutata viperina, Protobothrops cornutus, P. jerdonii, P. kaulbacki, P. mangshanensis, P. maolanensis, P. mucrosquamatus, P. xiangchengensis, Psammody‐ nastes pulverulentus, Psammophis lineolatus, Pseudolaticauda semifasciata, Pseudoxenodon bambusi‐ cola, P. karlschmidti, P. macrops, P. stejnegeri, Ptyas carinata, P. dhumnades, P. korros, P. mucosa, P. nigromarginata, Python bivittatus, Ramphotyphlops albiceps, R. braminus, R. lineatus, Rhabdophis adleri, R. callichroma, R. chrysargos, R. himalayanus, R. leonardi, R. nigrocinctus, R. nuchalis, R. subminiatus, R. swinhonis, R. tigrinus, Rhabdops bicolor, Rhadinophis frenatus, R. prasinus, Rhyn‐ chophis boulengeri, Sibynophis chinensis, S. collaris, Sinomicrurus hatori, S. kelloggi, S. macclellandi, Sinonatrix aequifasciata, S. annularis, S. percarinata, S. yunnanensis, Thermophis baileyi, T. zhaoermii, Trachischium monticola, T. tenuiceps, Trimeresurus albolabris, T. gracilis, T. gramineus, T. gum‐ prechti, T. medoensis, T. sichuanensis, T. stejnegeri, T. tibetanus, T. yunnanensis, Typhlops diardii, T. koshunensis, T. lazelli, Vipera berus, V. renardi, V. sachalinensis, V. ursinii, Xenochrophis flavipunc‐ tatus, X. piscator, X. hainanensis, X. unicolor.*

**Author details**

Chengdu Sichuan, PR China

30/3/399.full.pdf

2003.

Press; 1959.

1365-2699.2007.01843.x/pdf.

, Xianguang Guo2

\*Address all correspondence to: guoxg@cib.ac.cn

ect.com/science/article/pii/S0169534704002745

and Yuezhao Wang2

2 Department of Herpetology, Chengdu Institute of Biology, Chinese Academy of Sciences,

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

191

http://dx.doi.org/10.5772/55719

[1] Wiens JJ, Donoghue MJ. Historical biogeography, ecology, and species richness. Trends in Ecology and Evolution 2004; 19(12): 639-644. http://www.sciencedir‐

[2] Crisci JV, Sala OE, Katinas L, Posadas, P. Bridging historical and ecological ap‐ proaches 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/

[4] Posadas P, Crisci JV, Katinas L. Historical biogeography: A review of its basic con‐ cepts and critical issues. Journal of Arid Environments 2006; 66(3): 389-403. http://

[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.

[6] Hsü KJ, Chen H. Geological Atlas of China: An Application of the Tectonic Facies

[7] Wan T. The Outline of Geotectonics of China. Beijing: Geological Publishing House;

[8] Chinese Academy of Sciences, Nature Division Committee. China Animal Geogra‐ phy Division and China Insect Geography Division (First Draft). Beijing: Science

www.sciencedirect.com/science/article/pii/S0140196306000346

Concept to the Geology of China. New York: Elsevier Science; 1999.

[9] Zhang R. Zoogeography of China. Beijing: Science Press; 1999.

1 Guangxi Botanical Garden of Medicinal Plants, Nanning Guangxi, PR China

Yong Huang1

**References**

#### **Acknowledgements**

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).

### **Author details**

*B. multomaculata, B. nigriceps, Bungarus bungaroides, B. fasciatus, B. multicinctus, Calamaria buchi, C. pavimentata, C. septentrionalis, C. yunnanensis, Calliophis maculiceps, Chitulia inornata, C. ornate, C. torquata, Chrysopelea ornate, Coelognathus radiates, Cyclophiops doriae, C. major, C. multicinctus, C. ruffus, Daboia siamensis, Deinagkistrodon acutus, Dendrelaphis biloreatus, D.hollinrakei, D. ngansonensis, D. pictus, D. subocularis, Dinodon flavozonatum, D. rosozonatum, D. rufozonatum, D. septentrionalis, Disteira stokesii, Elaphe anomala, E. bimaculata, E. carinata, E. davidi, E. dione, E. schrenckii, E. zoigeensis, Emydocephalus ijimae, Enhydris bennettii, E. bocourti, E. chinensis, E.en‐ hydris, E. plumbea, Eryx miliaris, E. tataricus, Euprepiophis mandarinus, E. perlacea, Gloydius blomhoffii, G. brevicaudus, G. halys, G. intermedius, G. lijianlii, G. monticola, G. saxatilis, G. shedaoensis, G. strauchi, G. ussuriensis, Hemorrhois ravergieri, Hierophis spinalis, Hydrophis fasciatus, H. gracilis, H. parviceps, Kerilia jerdonii, Lapemis curtus, L. hardwickii, L. colubrine, L. laticaudata, Leioselasma cyanocincta, L. melanocephala, L. spiralis, Liopeltis frenatus, Lycodon capucinus, L. fasciatus, L. futsingensis, L. gongshan, L. laoensis, L. liuchengchaoi, L. ruhstrati, L. subcinctus, L. synaptor, Macropisthodon rudis, Naja atra, N. kaouthia, Natrix natrix, N. tessellate, Oligodon albocinctus, O. catenatus, O. chinensis, O. cinereus, O. cyclurus, O. eberhardti, O. formo‐ sanus, O. lacroixi, O. lungshenensis, O. melanozonatus, O. multizonatus, O. ningshaanensis, O. ornatus, O. taeniatus, Oocatochus rufodorsatus, Ophiophagus Hannah, Opisthotropis andersonii, O. balteata, O. cheni, O. guangxiensis, O. jacobi, O. kuatunensis, O. lateralis, O. latouchii, O. maculosa, O. maxwelli, Oreocryptophis porphyraceus, Orthriophis hodgsoni, O. moellendorffi, O. taeniurus, Ovophis monticola, O. tonkinensis, O. zayuensis, Paratapinophis praemaxillaris, Pareas boulengeri, P. carinatus, P. chinensis, P. formosensis, P. hamptoni, P. margaritophorus, P. monticola, P. nigriceps, P. stanleyi, Pelamis platura, Plagiopholis blakewayi, P. nuchalis, P. styani, P. unipostocularis, Polyodontognathus caerulescens, Praescutata viperina, Protobothrops cornutus, P. jerdonii, P. kaulbacki, P. mangshanensis, P. maolanensis, P. mucrosquamatus, P. xiangchengensis, Psammody‐ nastes pulverulentus, Psammophis lineolatus, Pseudolaticauda semifasciata, Pseudoxenodon bambusi‐ cola, P. karlschmidti, P. macrops, P. stejnegeri, Ptyas carinata, P. dhumnades, P. korros, P. mucosa, P. nigromarginata, Python bivittatus, Ramphotyphlops albiceps, R. braminus, R. lineatus, Rhabdophis adleri, R. callichroma, R. chrysargos, R. himalayanus, R. leonardi, R. nigrocinctus, R. nuchalis, R. subminiatus, R. swinhonis, R. tigrinus, Rhabdops bicolor, Rhadinophis frenatus, R. prasinus, Rhyn‐ chophis boulengeri, Sibynophis chinensis, S. collaris, Sinomicrurus hatori, S. kelloggi, S. macclellandi, Sinonatrix aequifasciata, S. annularis, S. percarinata, S. yunnanensis, Thermophis baileyi, T. zhaoermii, Trachischium monticola, T. tenuiceps, Trimeresurus albolabris, T. gracilis, T. gramineus, T. gum‐ prechti, T. medoensis, T. sichuanensis, T. stejnegeri, T. tibetanus, T. yunnanensis, Typhlops diardii, T. koshunensis, T. lazelli, Vipera berus, V. renardi, V. sachalinensis, V. ursinii, Xenochrophis flavipunc‐*

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).

*tatus, X. piscator, X. hainanensis, X. unicolor.*

**Acknowledgements**

190 Current Progress in Biological Research

Yong Huang1 , Xianguang Guo2 and Yuezhao Wang2


2 Department of Herpetology, Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu Sichuan, PR China

### **References**


[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.

ces and track analysis. Biological Journal of the Linnean Society 2007; 92(3): 405-417.

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

193

http://dx.doi.org/10.5772/55719

[23] Melo Santos AM, Cavalcanti DR, da Silva JC, Tabarelli M. Biogeographical relation‐ ships 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/

[24] Posadas P. Distributional patterns of vascular plants in Tierra del Fuego: a study ap‐ plying parsimony analysis of endemicity (PAE). Biogeographica 1996; 72: 161-177.

[25] Zhao E, Zhao K, Zhou K. Fauna Sinica, Reptilia, Vol. 2. (Squamates, Lacertilia). Bei‐

[26] Zhao E, Adler K. Herpetology of China. Oxford: Ohio: Society for the Study of Am‐ phibians and Reptiles in cooperation with Chinese Society for the Study of Amphib‐

[27] Zhao E, Huang M, Zong Y. Fauna Sinica Reptilia, Vol.3 Squamates Serpentes. Beijing:

[28] Zhao E. Snake of China. Hefei: Anhui Science and Technology Publishing House;

[29] Swofford DL. PAUP\*. Phylogenetic Analysis Using Parsimony (\*and Other Meth‐ ods). Version 4. Sinauer Associates, Sunderland, MA; 1998. http://www.sinauer.com/

[30] Glasby CJ, Álvarez B. Distribution patterns and biogeographic analysis of Austral Polychaeta (Annelida). Journal of Biogeography 1999; 26(3): 507–533. Http://onlineli‐

[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

[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):

[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/

[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#

brary.wiley.com/doi/10.1046/j.1365-2699.1999.00297.x/pdf

186-194. www.jstor.org/stable/1566107

10.1007%2Fs10531-005-4382-4

Herpetology 2000; 34(1): 54-60. http://www.jstor.org/stable/1565238

http://onlinelibrary.wiley.com/doi/10.1111/j.1095-8312.2007.00844.x/full

full

jing: Science press; 1999.

ians and Reptiles; 1993.

Science press; 1999.

detail.php?id=8060

2006.


ces 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 relation‐ ships 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

[10] MacKinnon J, Meng S, Cheung C, Carey G, Zhu X, Melville D. A Biodiversity Review

[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/

[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.sciencedir‐

[13] Morrone JJ, Escalante T. Parsimony analysis of endemicity (PAE) of Mexican terres‐ trial mammals at different area units: when size matters. Journal of Biogeography 2002; 29(8):1095-1104. http://onlinelibrary.wiley.com/doi/10.1046/j.

[14] Humphries CJ, Parenti LR. Cladistic Biogeography, 2nd edn. Oxford Biogeography

[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;

[17] Cracraft J. Patterns of diversification within continental biotas: hierarchical congru‐ ence among the areas of endemism of Australian vertebrates. Australian Systematic

[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://

[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.

[20] Brooks DR, van Veller MGP. Critique of parsimony analysis of endemicity as a meth‐ od of historical biogeography. Journal of Biogeography 2003; 30(6): 819-825. http://

[21] Nihei SS. Misconceptions about parsimony analysis of endemicity. Journal of Bio‐ geography 2006; 33(12): 2099-2106. http://onlinelibrary.wiley.com/doi/10.1111/j.

[22] Contreras-Medina R, Vega IL, Morrone JJ. Application of parsimony analysis of en‐ demicity to Mexican gymnosperm distributions: grid-cells, biogeographical provin‐

Botany 1991; 4(1): 211-227. http://www.publish.csiro.au/paper/SB9910211.htm

onlinelibrary.wiley.com/doi/10.1046/j.1365-2699.2001.00576.x/full

onlinelibrary.wiley.com/doi/10.1046/j.1365-2699.2003.00848.x/full

of China. Hong Kong: World Wide Fund for Nature International; 1996.

p1355416n7h87706/

192 Current Progress in Biological Research

1365-2699.2002.00753.x/full

1095-8312.2000.tb01265.x/pdf

1365-2699.2006.01619.x/full

ect.com/science/article/pii/S1002007108001457

Series 12, Oxford: Oxford University Press; 1999.

43(3): 438-441. http://www.jstor.org/stable/2413679


[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#

[47] Wright DH. Species–energy theory: An extension of species–area theory. Oikos 1983;

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

195

http://dx.doi.org/10.5772/55719

[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.

[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.

[50] Hortal J, Triantis KA, Meiri S, Thébault E, Sfenthourakis S. Island species richness in‐ creases with habitat diversity. The American Naturalist 2009; 174(6): E205-E217.

[51] McCain CM. Global analysis of reptile elevational diversity. Global Ecology and Bio‐ geography 2010; 19 (4): 541–553. http://onlinelibrary.wiley.com/doi/10.1111/j.

[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://onlineli‐

[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/

[54] Stockwell D, Peters D. The GARP modelling system: problems and solutions to auto‐ mated 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' distribu‐ tions from small numbers of occurrence records: A test case using cryptic geckos in Madagascar. Journal of Biogeography 2007; 34(1): 102-117. http://onlineli‐

[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.

[57] Peterson AT, Robins CR. Using ecological-niche modeling to predict barred owl in‐ vasions 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/

[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://

http://onlinelibrary.wiley.com/doi/10.1111/j.0906-7590.2007.05102.x/full

http://onlinelibrary.wiley.com/doi/10.1111/j.0906-7590.2004.03732.x/full

41(3): 496-506. http://www.jstor.org/stable/3544109

http://www.jstor.org/stable/10.1086/645085

journal/v398/n6728/full/398611a0.html

brary.wiley.com/doi/10.1111/j.0906-7590.2007.05025.x/full

brary.wiley.com/doi/10.1111/j.1365-2699.2006.01594.x/full

www.bioone.org/doi/full/10.1614/P2002-081

1472-4642.2007.00369.x/pdf

1466-8238.2010.00528.x/pdf

full


[47] Wright DH. Species–energy theory: An extension of species–area theory. Oikos 1983; 41(3): 496-506. http://www.jstor.org/stable/3544109

[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):

[36] Gaston KJ. Global patterns in biodiversity. Nature 2000; 405(6783): 220-227. http://

[37] von Humboldt A. Ansichten der Natur mit wissenschaftlichen Erlauterungen. Ger‐

[38] Hawkins BA, Field R, Cornell HV, Currie DJ, Guegan JF, Kaufman DM, Kerr JT, Mit‐ telbach 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.

[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 con‐ trasting 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 rich‐ ness in high-energy regions. Nature 1997; 385(6613): 252-254. http://

[42] Kerr JT, Southwood TRE, Cihlar J. Remotely sensed habitat diversity predicts butter‐ fly species richness and community similarity in Canada. Proceedings of the Nation‐ al Academy of Sciences of the United States of America 2001; 98(20): 11365-11370.

[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

[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://

[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.

[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

http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2311.2003.00551.x/full

2001; 98(8):4534-4539. http://www.pnas.org/content/98/8/4534.full.pdf+html

www.nature.com/nature/journal/v385/n6613/abs/385252a0.html

http://www.pnas.org/content/98/20/11365.full.pdf+html

www.jstor.org/stable/3565855

343-353. http://www.zoores.ac.cn/EN/abstract/abstract3193.shtml#

www.nature.com/nature/journal/v405/n6783/pdf/405220a0.pdf

http://www.esajournals.org/doi/full/10.1890/03-8006

many: J. G. Cotta, Tübingen; 1808.

194 Current Progress in Biological Research


[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.sciencedir‐ ect.com/science/article/pii/S1617138108000484

[71] Burnham KP, Anderson DR. Model Selection and Inference: A Practical Information

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

197

http://dx.doi.org/10.5772/55719

[72] Burnham KP, Anderson DR. Multimodel inference—understanding AIC and BIC in model selection. Sociological Methods and Research 2004; 33(2): 261-304. http://

[73] Powney GD, Grenyer R, Orme CDL, Owens IPF, Meiri S. Hot, dry and different Aus‐ tralian lizard richness is unlike that of mammals, amphibians, and birds. Global Ecol‐ ogy and Biogeography 2010; 19(3): 386-396. http://onlinelibrary.wiley.com/doi/

[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://

[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/

[76] Willis KJ, Whittaker RJ. Species diversity–scale matters. Science 2002; 295(5558):

[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/

[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. Na‐ ture 1997; 386: 788-795. http://www.nature.com/nature/journal/v386/n6627/abs/

[79] Shi Y, Zhao J, Wang J. New Understanding of Quaternary Glaciations in China.

[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

[81] Zhou S, Li J, Zhao J, Wang J, Zheng J. Quaternary Glaciations: Extent and Chronolo‐ gy in China. In: Ehlers J, Gibbard PL, Hughes PD. (eds) Quaternary Glaciations-Ex‐ tent and Chronology, Volume 15: A closer look (Developments in Quaternary Science). Amsterdam: Elsevier Press; 2011. p981-1002. Available from http://

[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.sciencedir‐

www.sciencedirect.com/science/article/pii/B9780444534477000702

onlinelibrary.wiley.com/doi/10.1111/j.1469-7998.2000.tb00605.x/pdf

1245-1248. http://www.sciencemag.org/content/295/5558/1245.full

Theoretic Approach. New York: Springer-Verlag; 1998.

smr.sagepub.com/content/33/2/261.full.pdf+html

10.1111/j.1466-8238.2009.00521.x/full

10.1111/j.1365-2699.1997.tb00049.x/pdf

Shanghai: Shanghai Popular Science Press; 2011.

ect.com/science/article/pii/S1002007109000513

full/416159a.html

386788a0.html

2011; 66(7): 867-884.


[71] Burnham KP, Anderson DR. Model Selection and Inference: A Practical Information Theoretic Approach. New York: Springer-Verlag; 1998.

[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.sciencedir‐

[60] Menon S, Choudhury BI, Khan LM, Peterson AT. Ecological niche modeling and lo‐ cal knowledge predict new populations of Gymnocladus assamicus a critically en‐ dangered tree species. Endangered Species Research 2010; 11: 175-181. http://

[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 spe‐ cies in Madagascar. Nature 2003; 426(6968): 837-841. http://www.nature.com/nature/

[62] Anderson RP, Lew D, Peterson AT. Evaluating predictive models of species' distri‐ butions: criteria for selecting optimal models. Ecological Modelling 2003; 162(3):

[63] Ahn CH, Tateishi R. Development of a global 30-minute grid potential evapotranspi‐ ration data set. Journal of the Japan Society Photogrammetry and Remote Sensing

[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/arti‐

[65] Tabachnick BG, Fidell LS. Using Multivariate Statistics, 4th edn. Boston: Allyn and

[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):

[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.

[69] Rangel TF, Diniz-Filho JAF, Bini LM. SAM: A comprehensive application for spatial analysis in macroecology. Ecography 2010; 33(1): 46-50. http://onlineli‐

[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://

211-232. http://www.sciencedirect.com/science/article/pii/S0304380002003496

ect.com/science/article/pii/S1617138108000484

www.int-res.com/articles/esr2010/11/n011p175.pdf

journal/v426/n6968/full/nature02205.html

1659-1673. http://www.jstor.org/stable/1939924

brary.wiley.com/doi/10.1111/j.1600-0587.2009.06299.x/full

onlinelibrary.wiley.com/doi/10.1046/j.1466-822X.2003.00322.x/full

1994; 33(2): 12-21.

196 Current Progress in Biological Research

Bacon, MA; 2000.

cles/cr2002/21/c021p001.pdf

1600-0587.2002.250508.x/full


[83] Barabanov AV, Ananjeva NB. Catalogue of the available scientific species-group names for lizards of the genus *Phrynocephalus* Kaup, 1825 (Reptilia, Sauria, Agami‐ dae). Zootaxa 2007; 1399: 1-56. http://www.mapress.com/zootaxa/2007f/ z01400p068f.pdf

complex *Phrynocephalus* superspecies *helioscopus* (Pallas 1771) (Reptilia: Agamidae). Russian Journal of Genetics 2011; 47(7): 842-856. http://www.springerlink.com/

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

199

http://dx.doi.org/10.5772/55719

[95] Dunayev EA, Ivanova N, Poyarkov NA, Borisenko A, Dujsebayeva T, Hebert PD. Molecular perspective on the evolution and barcoding of toad-headed agamas (ge‐ nus *Phrynocephalus*, Agamidae) in Middle Asia. Programme and Abstracts of 14th Eu‐ ropean Congress of Herpetology and SHE Ordinary General Meeting, 19-23

[96] Jin Y, Brown RP, Liu N. Cladogenesis and phylogeography of the lizard *Phrynocepha‐ lus 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 contempo‐ rary 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.

[98] Guo X, Liu L, Wang Y. Phylogeography of the *Phrynocephalus vlangalii* species com‐ plex 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-jour‐

[99] Jin Y, Liu N. Phylogeography of *Phrynocephalus erythrurus* from the Qiangtang Pla‐ teau 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://

[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.science‐

[102] Szczerbak NN. Guide to the Reptiles of the Eastern Palearctic. USA: Krieger Publish‐

[103] Szczerbak NN. Yashchurki Palearktiki (*Eremias* lizards of the Palearctic). Axadeimiya

[104] Arnold EN. The hemipenis of lacertid lizards (Reptilia: Lacertidae): structure, varia‐ tion and systematic implications. Journal of Natural History 1986; 20: 1221-1257.

[105] Arnold EN. Towards the phylogeny and biogeography of the Lacertidae: relation‐ ships within an Old-World family of lizards derived from morphology. Bulletin of

http://www.tandfonline.com/doi/abs/10.1080/00222938600770811#preview

Nauk Ukrainskoi USSR Institut Zoologii. Naukova Dumka, Kiev; 1974.

nal.com/index.php?module=case&act=List&ClassID=19

www.sciencedirect.com/science/article/pii/S1055790310004094

direct.com/science/article/pii/S0040195103001379

ing Company, Malabar, Florida; 2003.

content/h245501257785793/

1469-7998.2009.00564.x/full

September 2007, Porto (Portogal); 2007.


complex *Phrynocephalus* superspecies *helioscopus* (Pallas 1771) (Reptilia: Agamidae). Russian Journal of Genetics 2011; 47(7): 842-856. http://www.springerlink.com/ content/h245501257785793/

[83] Barabanov AV, Ananjeva NB. Catalogue of the available scientific species-group names for lizards of the genus *Phrynocephalus* Kaup, 1825 (Reptilia, Sauria, Agami‐ dae). Zootaxa 2007; 1399: 1-56. http://www.mapress.com/zootaxa/2007f/

[84] Uetz P. The Reptile Database, http://www.reptile-database.org (accessed Aug 3,

[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*, Aga‐ midae) based on morphology. Bulletin of the British Museum (Natural History) Zo‐

[87] Ananjeva NB, Tuniyev BS. Historical biogeography of the *Phrynocephalus* species of the USSR. Asiatic Herpetology Research 1992; 4: 76-98. http://www.asiatic-herpeto‐

[88] Macey JR, Ananjeva NB, Zhao E, Wang Y, Papenfuss TJ. An allozyme-based phylo‐ genetic hypothesis for *Phrynocephalus* (Agamidae) and its implications for the histori‐ cal 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

[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 Chi‐

[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-eval‐ uation. Molecular Phylogenetics and Evolution 2007; 45(2): 643-662. http://

[93] Melville J, Hale J, Mantziou G, Ananjeva NB, Milto K, Clemann N. Historical bio‐ geography, 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/

[94] Solovyeva EN, Poyarkov NA, Dunayev EA, Duysebayeva TN, Bannikova AA. Mo‐ lecular differentiation and taxonomy of the sunwatcher toad-headed agama species

First Asian Herpetological Meeting. Beijing: China Forestry Press; 1993.

nese *Phrynocephalus*. Acta Zoologica Sinica 1997; 43(4): 409-418.

www.sciencedirect.com/science/article/pii/S1055790307002175

z01400p068f.pdf

198 Current Progress in Biological Research

ology 1999; 65: 1-13.

S1055790309001754

logical.org/Archive/Volume%2004/04\_12.pdf

2012).


the British Museum (Natural History) Zoology 1989; 55(2): 209-257. http:// biostor.org/reference/107021

prints. Iranian Journal of Animal Biosystematics 2009; 5(1): 33-46. http://

Historical and Ecological Factors Affecting Regional Patterns of Endemism and Species Richness: The Case of

Squamates in China

201

http://dx.doi.org/10.5772/55719

[116] Warren DL, Glor RE, Turelli M. Environmental niche equivalency versus conserva‐ tism: 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

www.sid.ir/en/VEWSSID/J\_pdf/116420090105.pdf


prints. 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 conserva‐ tism: 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

the British Museum (Natural History) Zoology 1989; 55(2): 209-257. http://

[106] Boulenger GA. Monograph of the Lacertidae, vol. II. London: Trustees of the British

[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

[109] Wan L, Sun S, Jin Y, Yan Y, Liu N. Molecular phylogeography of the Chinese lacer‐ tids of the genus *Eremias* (Lacertidae) based on 16S rRNA mitochondrial DNA se‐ quences. Amphibia-Reptilia 2007; 28(1): 33-41.. http://www.ingentaconnect.com/

[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 Eu‐ ropean Congress of Herpetology and SHE Ordinary General Meeting, 19-23 Septem‐

[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. Pro‐ gramme and Abstracts of 15th European Congress of Herpetology and SHE Ordinary

[112] Guo X, Dai X, Chen D, PapenfussTJ, Ananjeva NB, Melnikov DA, Wang Y. Phyloge‐ ny and divergence times of some racerunner lizards (Lacertidae: *Eremias*) inferred from mitochondrial 16S rRNA gene segments. Molecular Phylogenetics and Evolu‐ tion 2011; 61(2): 400-412. http://www.sciencedirect.com/science/article/pii/

[113] Zhao Q, Liu HX, Luo L, Ji X. Comparative population genetics and phylogeography of two lacertid lizards (*Eremias argus* and *E. brenchleyi*) from China. Molecular Phylo‐ genetics and Evolution 2011; 58(3): 478-491. http://www.sciencedirect.com/science/

[114] Rastegar-Pouyani E, Noureini SK, Joger U, Wink M. Molecular phylogeny and intra‐ specific differentiation of the *Eremias velox* complex of the Iranian Plateau and Cen‐ tral Asia (Sauria, Lacertidae). Journal of Zoological Systematics and Evolutionary Research 2012; 50(3): 220-229. http://onlinelibrary.wiley.com/doi/10.1111/j.

[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 finger‐

General Meeting, 28 September-02 October 2009, Kuşadasi, Turkey; 2009.

[107] Fitzinger L. Systema Reptilium. Vienna: Fasciculus Primus; 1843.

content/brill/amre/2007/00000028/00000001/art00005

biostor.org/reference/107021

1971; 2: 48-55.

200 Current Progress in Biological Research

Museum of Natural History; 1921.

ber 2007, Porto (Portogal); 2007.

S1055790311003113

article/pii/S1055790310004987

1439-0469.2012.00662.x/pdf

**Chapter 9**

*In vitro* **Propagation of Critically Endangered Endemic**

*Rhaponticoides mykalea* **(Hub.-Mor.) by Axillary Shoot**

Turkey is one of the richest countries in variability of flora. It has nearly 9000 plant species about 3000 of which are endemic [1]. *Asteraceae*, is represented by 50 species in Turkey with an endemism of nearly 54% [2]. *Rhaponticoides mykalea* (Hub.-Mor.) M.V. Agab. & Greuter which belongs to the *Asteraceae* family, falls within the CR (Critically Endangered) category in the Red Data Book of Turkey [1]. While *R. mykalea* (Hub.-Mor.) was classified under the section *Centaurea* as *Centaurea mykalea* (Hub.-Mor.) before now. Today it has been separated from the section Centaurea [3]. It spreads very scarce in Kuşadası (Aydın), Muğla and Isparta, and faces with the danger of extinction. *R. mykalea* that has very limited number of individuals is under strong anthropogenic pressure such as the gradually increase in ongoing urbanization due to rapid developments of tourism sector, the conversion of natural habitats into human domi‐ nated lands, the over-grazing and collecting capitula of *R. mykalea* by local people for food. The species has already been under the threat of extinction and the situation above will increase the risk of extinction of this species even more [4]. For this reason, local protection measures

Nowadays, the conservation of wild plant genetic resources is very important for preventing a decrease in genetic variability. Conservation of the endemic or threatened plants is carried out using different strategies. *In vitro* culture is an efficient method for *ex situ* conservation of plant diversity [6,7], because many endangered species can be quickly propagated and preserved from a minimum of plant material with low impact on wild populations with this technology [8]. In recent years, there has been an increased interest in *in vitro* techniques that offer powerful tools for germplasm conservation and the mass multiplication of many

and reproduction in any medium, provided the original work is properly cited.

© 2013 Emek and Erdag; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

© 2013 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution,

**Proliferation**

Yelda Emek and Bengi Erdag

http://dx.doi.org/10.5772/55599

**1. Introduction**

Additional information is available at the end of the chapter

and global conservation strategies are necessary [5].
