The natural and human-mediated expansion of a human-commensal lizard into the fringes of Southeast Asia.

BACKGROUND: Human-commensal species often display deep ancestral genetic structure within their native range and founder-effects and/or evidence of multiple introductions and admixture in newly established areas. We investigated the phylogeography of Eutropis multifasciata, an abundant human-commensal scincid lizard that occurs across Southeast Asia, to determine the extent of its native range and to assess the sources and signatures of human introduction outside of the native range. We sequenced over 350 samples of E. multifasciata for the mitochondrial ND2 gene and reanalyzed a previous RADseq population genetic dataset in a phylogenetic framework. RESULTS: Nuclear and mitochondrial trees are concordant and show that E. multifasciata has retained high levels of genetic structure across Southeast Asia despite being frequently moved by humans. Lineage boundaries in the native range roughly correspond to several major biogeographic barriers, including Wallace's Line and the Isthmus of Kra. Islands at the outer fringe of the range show evidence of founder-effects and multiple introductions. CONCLUSIONS: Most of enormous range of E. multifasciata across Southeast Asia is native and it only displays signs of human-introduction or recent expansion along the eastern and northern fringe of its range. There were at least three events of human-introductions to Taiwan and offshore islands, and several oceanic islands in eastern Indonesia show a similar pattern. In Myanmar and Hainan, there is a founder-effect consistent with post-warming expansion after the last glacial maxima or human introduction.

How the African house gecko (Hemidactylus mabouia) conquered the world.

Alien species are among the greatest threats to biodiversity, but the evolutionary origins of invasiveness remain obscure. We conducted the first range-wide sampling of Hemidactylus mabouia from more than 120 localities across Africa, Madagascar and the Neotropics to understand the evolutionary history of one of the most widely distributed, invasive vertebrates in the world. We used a multi-locus phylogeny, species delimitation, fossil-calibrated timetree, ancestral area reconstruction and species distribution models (SDMs) to determine how many putative species-level lineages are contained within H. mabouia, the timing and tempo of diversification, and the origins of commensality-providing insights into the evolutionary origins of invasiveness. Our analyses suggest 'H. mabouia' originated in the Miocene in the Zambezian biogeographic region and includes as many as 20 putative species-level lineages, of which only Hemidactylus mabouia sensu stricto is invasive and widely distributed, including all Neotropical records. Zambezia is the hotspot for diversity within the group with 14 species in southeastern Zambezia. SDMs suggest that H. mabouia was able to establish in the Neotropics due to habitat suitability, and globalization and the slave trade probably allowed it to cross the Atlantic. Distribution models for the H. mabouia complex overpredict the range of the invasive H. mabouia sensu stricto-highlighting the importance of taxonomy in invasive species management.

Multilocus phylogeny of Bornean Bent-Toed geckos (Gekkonidae: Cyrtodactylus) reveals hidden diversity, taxonomic disarray, and novel biogeographic patterns.

The gekkonid genus Cyrtodactylus is a highly diverse group of lizards (280 + species), which covers an expansive geographic range. Although this genus has been the focus of many taxonomic and molecular systematic studies, species on the Southeast Asian island of Borneo have remained understudied, leading to an unclear evolutionary history with cascading effects on taxonomy and biogeographic inferences. We assembled the most comprehensive multilocus Bornean dataset (one mitochondrial and three nuclear loci) that included 129 novel sequences and representatives from each known Cyrtodactylus species on the island to validate taxonomic status, assess species diversity, and elucidate biogeographic patterns. Our results uncovered a high proportion of cryptic diversity and revealed numerous taxonomic complications, especially within the C. consobrinus, C. malayanus, and C. pubisulcus groups. Comparisons of pairwise genetic distances and a preliminary species delimitation analysis using the Automatic Barcode Gap Discovery (ABGD) method demonstrated that some wide-ranging species on Borneo likely comprise multiple distinct and deeply divergent lineages, each with more restricted distributional ranges. We also tested the prevailing biogeographic hypothesis of a single invasion from Borneo into the Philippines. Our analyses revealed that Philippine taxa were not monophyletic, but were likely derived from multiple separate invasions into the geopolitical areas comprising the Philippines. Although our investigation of Bornean Cyrtodactylus is the most comprehensive to-date, it highlights the need for expanded taxonomic sampling and suggests that our knowledge of the evolutionary history, systematics, and biogeography of Bornean Cyrtodactylus is far from complete.

Optimizing Phylogenomics with Rapidly Evolving Long Exons: Comparison with Anchored Hybrid Enrichment and Ultraconserved Elements.

Marker selection has emerged as an important component of phylogenomic study design due to rising concerns of the effects of gene tree estimation error, model misspecification, and data-type differences. Researchers must balance various trade-offs associated with locus length and evolutionary rate among other factors. The most commonly used reduced representation data sets for phylogenomics are ultraconserved elements (UCEs) and Anchored Hybrid Enrichment (AHE). Here, we introduce Rapidly Evolving Long Exon Capture (RELEC), a new set of loci that targets single exons that are both rapidly evolving (evolutionary rate faster than RAG1) and relatively long in length (>1,500 bp), while at the same time avoiding paralogy issues across amniotes. We compare the RELEC data set to UCEs and AHE in squamate reptiles by aligning and analyzing orthologous sequences from 17 squamate genomes, composed of 10 snakes and 7 lizards. The RELEC data set (179 loci) outperforms AHE and UCEs by maximizing per-locus genetic variation while maintaining presence and orthology across a range of evolutionary scales. RELEC markers show higher phylogenetic informativeness than UCE and AHE loci, and RELEC gene trees show greater similarity to the species tree than AHE or UCE gene trees. Furthermore, with fewer loci, RELEC remains computationally tractable for full Bayesian coalescent species tree analyses. We contrast RELEC to and discuss important aspects of comparable methods, and demonstrate how RELEC may be the most effective set of loci for resolving difficult nodes and rapid radiations. We provide several resources for capturing or extracting RELEC loci from other amniote groups.

Uncovering karst endemism within Borneo: two new Cyrtodactylus species from Sarawak, Malaysia.

The island of Borneo lies within one of the most biodiverse regions in the world. Despite this, its documented gekkonid diversity is not commensurate with other areas of Southeast Asia. The megadiverse genus Cyrtodactylus is especially underrepresented. Limestone-karst ecosystems, in particular, harbor many endemic Cyrtodactylus species, but only one karst-dwelling species is currently recognized from Borneo. This paper adds two additional karst-dwelling Cyrtodactylus species-C. muluensis sp. nov. and C. limajalur sp. nov.-from Sarawak, Malaysia. Cyrtodactylus muluensis sp. nov. is endemic to Gunung Mulu and is distinguished from its congeners by having a precloacal groove, 31-38 ventral scales, a maximum SVL of at least 88 mm, enlarged subcaudals, 19-20 subdigital lamellae, and a banded dorsal body pattern. Cyrtodactylus limajalur sp. nov. is endemic to the Serian region and is distinguished from its congeners by having 33-42 ventral scales, enlarged subcaudals, a precloacal pit, a maximum SVL of at least 94 mm, 5-6 enlarged femoral scales, 19-22 subdigital lamellae, and five distinct bands on the dorsum. Both species are phylogenetically distinct and deeply divergent from all other congeners. The description of two new karst-dwelling species highlights the need to conserve karst habitats and the endemic species they harbor.

A species-level phylogeny of Trachylepis (Scincidae: Mabuyinae) provides insight into their reproductive mode evolution.

Trachylepis (Mabuyinae) includes ∼80 species of fully-limbed skinks found primarily in Africa and Madagascar, but a robust species-level phylogeny for this genus is lacking and this impedes studies on a wide-range of topics from biogeography to character evolution. Trachylepis and its close relatives (which together form the Mabuya group or Mabuyinae) are notable in that they have undergone multiple transitions and remarkable specializations in their reproductive modes. A Trachylepis phylogeny will be particularly useful for investigating reproductive evolution, because it includes species that exhibit oviparity, viviparity, and bimodal parity (species with both oviparous and viviparous populations). We sequenced DNA at four mitochondrial and five nuclear loci for 67 (∼84% of) Trachylepis species to infer a phylogeny for this genus. We performed stochastic character mapping of parity mode under a variety of parity mode transition models to infer ancestral parity mode states and the number and type of parity mode transitions. We recovered a strongly supported phylogeny of Trachylepis that is generally consistent with earlier phylogenetic studies. The best-fit model of reproductive mode evolution supports an oviparous ancestor for Trachylepis, and supports at least three viviparity to oviparity transitions. We compared parity mode evolution under the overall best-fit model (no constraints on parity mode transitions) to the best-fit model among the subset of models that assume viviparity to oviparity transitions are impossible. Our results support a model of reproductive evolution that allows for reversibility from viviparity to oviparity, a process that is not generally accepted. Alternatively, the best-fit model of evolution among the set of models that eliminate reversals from viviparity to oviparity suggests that bimodal reproduction may have persisted for millions of years within multiple lineages.

Phylogenetic and morphological investigation of the Mochlus afer-sundevallii species complex (Squamata: Scincidae) across the arid corridor of sub-Saharan Africa.

The aridification of Africa resulted in the fragmentation of forests and the expansion of an arid corridor stretching from the northeast to southwest portion of sub-Saharan Africa, but the role this corridor has had in species-level diversification of southern African vertebrates is poorly understood. The skink species Mochlus afer and M. sundevallii inhabit wide areas of the arid corridor and are therefore an ideal species pair for studying patterns of genetic and phenotypic diversity associated with this landscape. However, species boundaries between these taxa have been controversial. Using multi-locus molecular and morphological datasets, we investigate diversification patterns of the M. afer-sundevallii Species Complex across the arid corridor. Although analyses of genetic data reveals some genetic structure among geographic populations, results of phylogenetic and morphological analyses provide little support for two distinct evolutionary lineages, suggesting that populations previously referred to as M. afer and M. sundevallii represent a single species, Mochlus sundevallii. Genetic diversity is unequally distributed across the arid corridor, with observed patterns consistent with aridification-facilitated diversification southward across southern Africa. Additional geographic and population-level sampling is necessary before more conclusive inferences can be drawn about the role historical climate transitions have played in skink diversification patterns across southern Africa.

Ancient divergence time estimates in Eutropis rugifera support the existence of Pleistocene barriers on the exposed Sunda Shelf.

Episodic sea level changes that repeatedly exposed and inundated the Sunda Shelf characterize the Pleistocene. Available evidence points to a more xeric central Sunda Shelf during periods of low sea levels, and despite the broad land connections that persisted during this time, some organisms are assumed to have faced barriers to dispersal between land-masses on the Sunda Shelf. Eutropis rugifera is a secretive, forest adapted scincid lizard that ranges across the Sunda Shelf. In this study, we sequenced one mitochondrial (ND2) and four nuclear (BRCA1, BRCA2, RAG1, and MC1R) markers and generated a time-calibrated phylogeny in BEAST to test whether divergence times between Sundaic populations of E. rugifera occurred during Pleistocene sea-level changes, or if they predate the Pleistocene. We find that E. rugifera shows pre-Pleistocene divergences between populations on different Sundaic land-masses. The earliest divergence within E. rugifera separates the Philippine samples from the Sundaic samples approximately 16 Ma; the Philippine populations thus cannot be considered conspecific with Sundaic congeners. Sundaic populations diverged approximately 6 Ma, and populations within Borneo from Sabah and Sarawak separated approximately 4.5 Ma in the early Pliocene, followed by further cladogenesis in Sarawak through the Pleistocene. Divergence of peninsular Malaysian populations from the Mentawai Archipelago occurred approximately 5 Ma. Separation among island populations from the Mentawai Archipelago likely dates to the Pliocene/Pleistocene boundary approximately 3.5 Ma, and our samples from peninsular Malaysia appear to coalesce in the middle Pleistocene, about 1 Ma. Coupled with the monophyly of these populations, these divergence times suggest that despite consistent land-connections between these regions throughout the Pleistocene E. rugifera still faced barriers to dispersal, which may be a result of environmental shifts that accompanied the sea-level changes.

The measure of success: geographic isolation promotes diversification in Pachydactylus geckos.

BACKGROUND: Geckos of the genus Pachydactylus and their close relatives comprise the most species-rich clade of lizards in sub-Saharan Africa. Many explanations have been offered to explain species richness patterns of clades. In the Pachydactylus group, one possible explanation is a history of diversification via geographic isolation. If geographic isolation has played a key role in facilitating diversification, then we expect species in more species-rich subclades to have smaller ranges than species in less diverse subclades. We also expect traits promoting geographic isolation to be correlated with small geographic ranges. In order to test these expectations, we performed phylogenetic analyses and tested for correlations among body size, habitat choice, range sizes, and diversification rates in the Pachydactylus group. RESULTS: Both body size and habitat use are inferred to have shifted multiple times across the phylogeny of the Pachydactylus group, with large size and generalist habitat use being ancestral for the group. Geographic range size is correlated with both of these traits. Small-bodied species have more restricted ranges than large-bodied species, and rock-dwelling species have more restricted ranges than either terrestrial or generalist species. Rock-dwelling and small body size are also associated with higher rates of diversification, and subclades retaining ancestral conditions for these traits are less species rich than subclades in which shifts to small body size and rocky habitat use have occurred. The phylogeny also illustrates inadequacies of the current taxonomy of the group. CONCLUSIONS: The results are consistent with a model in which lineages more likely to become geographically isolated diversify to a greater extent, although some patterns also resemble those expected of an adaptive radiation in which ecological divergence acts as a driver of speciation. Therefore, the Pachydactylus group may represent an intermediate between clades in which radiation is adaptive versus those in which it is non-adaptive.

A single origin of extreme matrotrophy in African mabuyine skinks.

Most mammals and approximately 20% of squamates (lizards and snakes) are viviparous, whereas all crocodilians, birds and turtles are oviparous. Viviparity evolved greater than 100 times in squamates, including multiple times in Mabuyinae (Reptilia: Scincidae), making this group ideal for studying the evolution of nutritional patterns associated with viviparity. Previous studies suggest that extreme matrotrophy, the support of virtually all of embryonic development by maternal nutrients, evolved as many as three times in Mabuyinae: in Neotropical Mabuyinae (63 species), Eumecia (2 species; Africa) and Trachylepis ivensii (Africa). However, no explicit phylogenetic hypotheses exist for understanding the evolution of extreme matrotrophy. Using multilocus DNA data, we inferred a species tree for Mabuyinae that implies that T. ivensii (here assigned to the resurrected genus Lubuya) is sister to Eumecia, suggesting that extreme matrotrophy evolved only once in African mabuyine skinks.

Resolving the higher-order phylogenetic relationships of the circumtropical Mabuya group (Squamata: Scincidae): An out-of-Asia diversification.

Despite an abundance of phylogenetic studies focused on intrageneric relationships of members of the Mabuya group, the intergeneric relationships have remained difficult to resolve. The most-persistent unresolved regions of the phylogeny of the group include: (1) the placement of the Middle-Eastern Trachylepis with respect to the Afro-Malagasy Trachylepis and its taxonomic status; (2) the phylogenetic position of the Cape Verdean Chioninia within the larger Mabuya group; (3) support for the placement of Dasia with respect to the entire group; and (4) the phylogenetic placement of Eutropis novemcarinata with respect to other Eutropis and Dasia. In this study, we include representatives of all these taxa as well as African Eumecia and Neotropical Mabuya. We seek to address these phylogenetic and systematic issues by generating a well-resolved and supported phylogeny for the Mabuya group as a whole that can be used to develop a stable taxonomy and reconstruct the geographic patterns of diversification within the group. To meet these goals, we built a large multi-locus dataset of 11 markers (nine nuclear and two mitochondrial), and performed concatenated and species tree analyses to generate a well-supported phylogeny for the group. Statistical topology tests reject the monophyly of Middle-Eastern Trachylepis with Afro-Malagasy Trachylepis, and to reflect monophyly we place the Middle-Eastern species into a previously described genus, Heremites. Cape-Verdean Chioninia are resolved as the strongly supported sister-group to Afro-Malagasy Trachylepis. Monophyly of the Southeast-Asian genera, Eutropis and Dasia, is not supported, with a clade composed of Dasia+Eutropis novemcarinata more closely related to the rest of the Mabuya group than to the remaining Eutropis. The phylogenetic position of E. novemcarinata renders Eutropis polyphyletic, and we therefore describe and place E. novemcarinata into a new monotypic genus, Toenayar, to preserve monophyly among the genera. In light of these novel findings, we review and discuss the historical biogeography of the entire Mabuya group.

Mitochondrial introgression via ancient hybridization, and systematics of the Australian endemic pygopodid gecko genus Delma.

Of the more than 1500 species of geckos found across six continents, few remain as unfamiliar as the pygopodids - Family Pygopodidae (Gray, 1845). These gekkotans are limited to Australia (44 species) and New Guinea (2 species), but have diverged extensively into the most ecologically diverse limbless radiation save Serpentes. Current phylogenetic understanding of the family has relied almost exclusively on two works, which have produced and synthesized an immense amount of morphological, geographical, and molecular data. However, current interspecific relationships within the largest genus Delma Gray 1831 are based chiefly upon data from two mitochondrial loci (16s, ND2). Here, we reevaluate the interspecific relationships within the genus Delma using two mitochondrial and four nuclear loci (RAG1, MXRA5, MOS, DYNLL1), and identify points of strong conflict between nuclear and mitochondrial genomic data. We address mito-nuclear discordance, and remedy this conflict by recognizing several points of mitochondrial introgression as the result of ancient hybridization events. Owing to the legacy value and intraspecific informativeness, we suggest the continued use of ND2 as a phylogenetic marker. Results identify strong support for species groups, but relationships among these clades, and the placement of several enigmatic taxa remain uncertain. We suggest a more careful review of Delma australis and the 'northwest Australia' clade. Accurately assessing and addressing species richness and relationships within this endemic Australian Gekkotan genus is relevant for understanding patterns of squamate speciation across the region.

A molecular phylogeny of Afromontane dwarf geckos (Lygodactylus) reveals a single radiation and increased species diversity in a South African montane center of endemism.

Afromontane habitats throughout eastern sub-Saharan Africa support remarkable levels of microendemism. However, despite being the subject of decades of research interest, biogeographical patterns of diversification throughout this disjunct montane system still remain largely unknown. We examined the evolutionary relationships of diurnal dwarf geckos (Lygodactylus) from several Afromontane regions throughout southeastern Africa, focusing primarily on two species groups (rex and bonsi groups). Using both mitochondrial and nuclear markers, we generate a molecular phylogeny containing all members of the rex and bonsi groups, to evaluate the monophyly of these groups along with previous biogeographic hypotheses suggesting independent southward invasions into the greater Drakensberg Afromontane center of endemism in northeastern South Africa by each group. Our results provide no support for these taxonomic and biogeographic hypotheses, and instead reveal geographically circumscribed patterns of diversification. One clade is restricted to the highlands of southern Malawi and northern Mozambique and the other to the greater Drakensberg region of northeastern South Africa and Swaziland. Interestingly, L. bernardi from the Nyanga Highlands of eastern Zimbabwe is nested within the primarily savanna-dwelling capensis group. We use Bayesian species delimitation methods to evaluate species limits within the greater Drakensberg clade, which support the elevation of the subspecies of L. ocellatus and L. nigropunctatus, thus bringing the total to eight species within a relatively confined geographic area. These results further highlight the greater Drakensberg Afromontane region as both an important center of endemism, as well as a center of diversification contributing to the accumulation of southern Africa's rich species diversity.

Insights into Himalayan biogeography from geckos: a molecular phylogeny of Cyrtodactylus (Squamata: Gekkonidae).

The India-Asia collision profoundly influenced the climate, topography and biodiversity of Asia, causing the formation of the biodiverse Himalayas. The species-rich gekkonid genus Cyrtodactylus is an ideal clade for exploring the biological impacts of the India-Asia collision, as previous phylogenetic hypotheses suggest basal divergences occurred within the Himalayas and Indo-Burma during the Eocene. To this end, we sampled for Cyrtodactylus across Indian areas of the Himalayas and Indo-Burma Hotspots and used three genes to reconstruct relationships and estimate divergence times. Basal divergences in Cyrtodactylus, Hemidactylus and the Palaearctic naked-toed geckos were simultaneous with or just preceded the start of the India-Asia collision. Diversification within Cyrtodactylus tracks the India-Asia collision and subsequent geological events. A number of geographically concordant clades are resolved within Indo-Burmese Cyrtodactylus. Our study reveals 17 divergent lineages that may represent undescribed species, underscoring the previously undocumented diversity of the region. The importance of rocky habitats for Cyrtodactylus indicates the Indo-Gangetic flood plains and the Garo-Rajmahal Gap are likely to have been important historical barriers for this group.

A phylogenetic analysis of the southern African gecko genus Afroedura Loveridge (Squamata: Gekkonidae), with the description of nine new species from Limpopo and Mpumalanga provinces of South Africa.

A molecular phylogeny of the largely rupicolous geckos of the gekkonid genus Afroedura is presented based on a combination of mitochondrial and nuclear gene sequence data. Previously recognized species groups are only partly recovered, with A. pondolia retrieved as very distantly related to the congeners to which it was previously considered allied. Afroedura hawequensis forms a monotypic group that is sister to all other species, which are allocated to the A. nivaria, A. transvaalica, A. africana, A. multiporis, A. marleyi, and A. langi groups. The taxonomic status of species occurring in the former Transvaal province of South Africa (now Limpopo and Mpumalanga provinces) is reviewed based on morphological and molecular evidence. Afroedura rupestris sp. nov. is described in the A. multiporis group, A. maripi sp. nov., A. pongola sp. nov., and A. rondavelica sp. nov., are described in the A. marleyi group, and A. broadleyi sp. nov., A. granitica sp. nov., A. leoloensis sp. nov., A. pienaari sp. nov., and A. waterbergensis sp. nov. are described in the A. langi group. In addition, A. haackei, A. namaquensis, and A. tirasensis, are all raised to specific status. The description of nine new species highlights Limpopo and Mpumalanga provinces as previously unrecognized centers of diversity for the genus Afroedura. A key to the species of Afroedura is provided. 

A new species of lizard in the genus Caledoniscincus (Reptilia: Scincidae) from southern New Caledonia and a review of Caledoniscincus atropunctatus (Roux).

A new species of skink, Caledoniscincus notialis sp. nov., is described from the ultramafic ranges in southern New Caledonia. It is most similar to, and has previously been referred to Caledoniscincus atropunctatus (Roux), a species with a widespread distribution throughout the Grand Terre and Loyalty Islands. The new species is distinct genetically from all other Caledoniscincus and can be distinguished by a unique pattern of dorsal coloration in males. Its range extends over much of the extensive ultramafic block in the south of the island, including the Goro Plateau and the mountain ranges at the southern edge of the Chaine Centrale north to Mt. Humboldt. It occurs mainly in humid forest habitat, much of which is now present only as isolated fragments in southern New Caledonia. Its preference for a habitat that has undergone a significant reduction in area of extent in a region under pressure from a range of anthropogenic threats suggests this new species is of conservation concern, and could be ranked as Vulnerable under IUCN listing. The genetic relationships of a redefined Caledoniscincus atropunctatus identifies two major subgroups, one located mainly in the northern and central-east regions of Grand Terre and the other in the southern and central-west regions and also including the population on the Loyalty Islands.

A preliminary phylogeny of the Palearctic naked-toed geckos (Reptilia: Squamata: Gekkonidae) with taxonomic implications.

Palearctic naked-toed geckos are a group of gekkonid geckos that range from North Africa to northern India and western China, with their greatest diversity in Iran and Pakistan. Relationships among the constituent genera remain incompletely resolved and the monophyly of key genera remains unverified. Further, competing classifications are in current use and many species have been allocated to different genera by different authors. We used both mitochondrial (ND2) and nuclear genes (RAG1, PDC) to explore relationships among representatives of all but one genus in the group (Rhinogecko), including four genera not previously included in phylogenetic analyses (Asiocolotes, Altigekko, Indogekko, and Siwaligekko). Siwaligekko (and presumably other Tibeto-Himalayan species often referred to Cyrtopodion) are more closely related to tropical Asian Cyrtodactylus than to Palearctic naked-toed geckos. Sampled species of Asiocolotes and Altigekko are sister taxa, but both genera are here considered junior subjective synonyms of Altiphylax. Cyrtopodion sensu lato is non-monophyletic; Mediodactylus and Tenuidactylus, which have variably been considered as subgenera or synonyms of Cyrtopodion are both valid genera. Indogekko is embedded within Cyrtopodion and is here treated as a subgenus. Bunopus and Crossobamon are closely related to one-another, and with Agamura are interdigitated among taxa previously assigned to Cyrtopodion. Our data confirm the previous identification of a Saharo-Arabian Stenodactylus/Tropiocolotes/Pseudoceramodactylus clade and verify that Microgecko and Alsophylax are not members of the main clade of Palearctic naked-toed geckos. Osteological differences between Tropiocolotes and Microgecko, formerly treated as congeneric, are discussed and illustrated. The divergence between Cyrtodactylus and the Palearctic naked-toed clade predates the initial collision of the Indian and Eurasian plates, but deeper divergences within both groups are consistent with mountain building in the Himalayas and adjacent ranges as promoting cladogenic events. Miocene divergences within Tenuidactylus are consistent with vicariant speciation caused by uplift events in the Iranian and Transcaspian regions. Taxonomic implications of our phylogenetic results are discussed and a preliminary allocation of all species of padless Palearctic gekkonids to genus is provided.

Phylogeny of bent-toed geckos (Cyrtodactylus) reveals a west to east pattern of diversification.

The Asian/Pacific genus Cyrtodactylus is the most diverse and among the most widely distributed genera of geckos, and more species are continually being discovered. Major patterns in the evolutionary history of Cyrtodactylus have remained largely unknown because no published study has broadly sampled across the geographic range and morphological diversity of the genus. We assembled a data set including sequences from one mitochondrial and three nuclear loci for 68 Cyrtodactylus and 20 other gekkotan species to infer phylogenetic relationships within the genus and identify major biogeographic patterns. Our results indicate that Cyrtodactylus is monophyletic, but only if the Indian/Sri Lankan species sometimes recognized as Geckoella are included. Basal divergences divide Cyrtodactylus into three well-supported groups: the single species C. tibetanus, a clade of Myanmar/southern Himalayan species, and a large clade including all other Cyrtodactylus plus Geckoella. Within the largest major clade are several well-supported subclades, with separate subclades being most diverse in Thailand, Eastern Indochina, the Sunda region, the Papuan region, and the Philippines, respectively. The phylogenetic results, along with molecular clock and ancestral area analyses, show Cyrtodactylus to have originated in the circum-Himalayan region just after the Cretaceous/Paleogene boundary, with a generally west to east pattern of colonization and diversification progressing through the Cenozoic. Wallacean species are derived from within a Sundaland radiation, the Philippines were colonized from Borneo, and Australia was colonized twice, once via New Guinea and once via the Lesser Sundas. Overall, these results are consistent with past suggestions of a Palearctic origin for Cyrtodactylus, and highlight the key role of geography in diversification of the genus.

Evolution of gliding in Southeast Asian geckos and other vertebrates is temporally congruent with dipterocarp forest development.

Gliding morphologies occur in diverse vertebrate lineages in Southeast Asian rainforests, including three gecko genera, plus frogs, snakes, agamid lizards and squirrels. It has been hypothesized that repeated evolution of gliding is related to the dominance of Asian rainforest tree floras by dipterocarps. For dipterocarps to have influenced the evolution of gliding in Southeast Asian vertebrates, gliding lineages must have Eocene or later origins. However, divergence times are not known for most lineages. To investigate the temporal pattern of Asian gliding vertebrate evolution, we performed phylogenetic and molecular clock analyses. New sequence data for geckos incorporate exemplars of each gliding genus (Cosymbotus, Luperosaurus and Ptychozoon), whereas analyses of other vertebrate lineages use existing sequence data. Stem ages of most gliding vertebrates, including all geckos, cluster in the time period when dipterocarps came to dominate Asian tropical forests. These results demonstrate that a gliding/dipterocarp correlation is temporally viable, and caution against the assumption of early origins for apomorphic taxa.

Repeated origin and loss of adhesive toepads in geckos.

Geckos are well known for their extraordinary clinging abilities and many species easily scale vertical or even inverted surfaces. This ability is enabled by a complex digital adhesive mechanism (adhesive toepads) that employs van der Waals based adhesion, augmented by frictional forces. Numerous morphological traits and behaviors have evolved to facilitate deployment of the adhesive mechanism, maximize adhesive force and enable release from the substrate. The complex digital morphologies that result allow geckos to interact with their environment in a novel fashion quite differently from most other lizards. Details of toepad morphology suggest multiple gains and losses of the adhesive mechanism, but lack of a comprehensive phylogeny has hindered efforts to determine how frequently adhesive toepads have been gained and lost. Here we present a multigene phylogeny of geckos, including 107 of 118 recognized genera, and determine that adhesive toepads have been gained and lost multiple times, and remarkably, with approximately equal frequency. The most likely hypothesis suggests that adhesive toepads evolved 11 times and were lost nine times. The overall external morphology of the toepad is strikingly similar in many lineages in which it is independently derived, but lineage-specific differences are evident, particularly regarding internal anatomy, with unique morphological patterns defining each independent derivation.