The distributions of amphibians, birds and mammals have underpinned global and local conservation priorities, and have been fundamental to our understanding of the determinants of global biodiversity. In contrast, the global distributions of reptiles, representing a third of terrestrial vertebrate diversity, have been unavailable. This prevented the incorporation of reptiles into conservation planning and biased our understanding of the underlying processes governing global vertebrate biodiversity. Here, we present and analyse the global distribution of 10,064 reptile species (99% of extant terrestrial species). We show that richness patterns of the other three tetrapod classes are good spatial surrogates for species richness of all reptiles combined and of snakes, but characterize diversity patterns of lizards and turtles poorly. Hotspots of total and endemic lizard richness overlap very little with those of other taxa. Moreover, existing protected areas, sites of biodiversity significance and global conservation schemes represent birds and mammals better than reptiles. We show that additional conservation actions are needed to effectively protect reptiles, particularly lizards and turtles. Adding reptile knowledge to a global complementarity conservation priority scheme identifies many locations that consequently become important. Notably, investing resources in some of the world’s arid, grassland and savannah habitats might be necessary to represent all terrestrial vertebrates efficiently

A phylogenetic tree for acrodont lizards (Chamaeleonidae and Agamidae) is established based on 1434 bases (1041 informative) of aligned DNA positions from a 1685-1778 base pair region of the mitochondrial genome. Sequences from three protein-coding genes (ND1, ND2, and COI) are combined with sequences from eight intervening tRNA genes for samples of 70 acrodont taxa and two outgroups. Parsimony analysis of nucleotide sequences identifies eight major clades in the Acrodonta. Most agamid lizards are placed into three distinct clades. One clade is composed of all taxa occurring in Australia and New Guinea; Physignathus cocincinus from Southeast Asia is the sister taxon to the Australia-New Guinea clade. A second clade is composed of taxa occurring from Tibet and the Indian Subcontinent east through South and East Asia. A third clade is composed of taxa occurring from Africa east through Arabia and West Asia to Tibet and the Indian Subcontinent. These three clades contain all agamid lizards except Uromastyx, Leiolepis, and Hydrosaurus, which represent three additional clades of the Agamidae. The Chamaeleonidae forms another clade weakly supported as the sister taxon to the Agamidae. All eight clades of the Acrodonta contain members occurring on land masses derived from Gondwanaland. A hypothesis of agamid lizards rafting with Gondwanan plates is examined statistically. This hypothesis suggests that the African/West Asian clade is of African or Indian origin, and the South Asian clade is either of Indian or Southeast Asian origin. The shortest tree suggests a possible African origin for the former and an Indian origin for the latter, but this result is not statistically robust. The Australia-New Guinea clade rafted with the Australia-New Guinea plate and forms the sister group to a Southeast Asian taxon that occurs on plates that broke from northern Australia-New Guinea. Other acrodont taxa are inferred to be associated with the plates of Afro-Arabia and Madagascar (Chameleonidae), India (Uromastyx), or southeast Asia (Hydrosaurus and Leiolepis). Introduction of different biotic elements to Asia by way of separate Gondwanan plates may be a major theme of Asian biogeography. Three historical events may be responsible for the sharp faunal barrier between Southeast Asia and Australia-New Guinea, known as Wallace's line: (1) primary vicariance caused by plate separations; (2) secondary contact of Southeast Asian plates with Eurasia, leading to dispersal from Eurasia into Southeast Asia, and (3) dispersal of the Indian fauna (after collision of that subcontinent) to Southeast Asia. Acrodont lizards show the first and third of these biogeographic patterns and anguid lizards exhibit the second pattern. Modern faunal diversity may be influenced primarily by historical events such as tectonic collisions and land bridge connections, which are expected to promote episodic turnover of continental faunas by introducing new faunal elements into an area. Repeated tectonic collisions may be one of the most important phenomena promoti...

Identifying and dating historical biological events is a fundamental goal of evolutionary biology, and recent analytical advances permit the modeling of factors known to affect both the accuracy and the precision of molecular date estimates. As the use of multilocus data sets becomes increasingly routine, it becomes more important to evaluate the potentially confounding effects of rate heterogeneity both within (e.g., codon positions) and among loci when estimating divergence times. Here, using Plestiodon lizards as a test case, we examine the effects of accommodating rate heterogeneity among data partitions on divergence time estimation. Plestiodon inhabits both East Asia and North America, yet both the geographic origin of the genus and timing of dispersal between the continents have been debated. For each of the eight independently evolving loci and a combined data set, we conduct single model and partitioned analyses. We found that extreme saturation has obscured the underlying rate of evolution in the mitochondrial DNA (mtDNA), resulting in severe underestimation of the rate in this locus. As a result, the age of the crown Plestiodon clade was overestimated by 15-17 Myr by the unpartitioned analysis of the combined loci data. However, the application of partition-specific models to the combined data resulted in ages that were fully congruent with those inferred by the individual nuclear loci. Although partitioning improved divergence date estimates of the mtDNA-only analysis, the ages were nonetheless overestimated, thus indicating an inadequacy of our current models to capture the complex nature of mtDNA evolution in over large time scales. Finally, the statistically incongruent age distributions inferred by the partitioned and unpartitioned analyses of the combined data support mutually exclusive hypotheses of the timing of intercontinental dispersal of Plestiodon from Asia to North America. Analyses that best capture the rate of evolution in the combined data set infer that this exchange occurred via Beringia ∼18.0-30 Ma.