Substrate specialisation drives an unexpectedly diverse radiation in barking geckos (Ptenopus: Gekkonidae).

Barking geckos (genus Ptenopus) are terrestrial, burrowing lizards endemic to southern Africa, currently with three recognised species. Two species are range-restricted (P. kochi and P. carpi) and display clear differences in substrate preference (soft sand vs. hard gravel). The third and most widespread species, P. garrulus, occurs on a variety of substrates of differing hardness, across potential geographic barriers, and over a steep climatic gradient. Variations in morphology and advertisement calls indicates that P. garrulus may be a species complex. Two subspecies of P. garrulus are currently recognised: P. g. maculatus and P. g. garrulus. To investigate species boundaries, we produced the first comprehensive phylogeny for the genus. We used a novel application of multiple regression on matrices models to assess multiple environmental drivers of diversification, as contrasted to isolation by distance. We show that P. kochi, P. carpi, and P. g. garrulus are valid species, but that P. g. maculatus is a paraphyletic complex of five previously unrecognised taxa. Specialisation onto different substrates was likely the main driver of divergence, with parapatric occurrence of two to four clades occurring at each of the three substrate transition zones identified a priori. The region encompasses diverse bioclimatic regions and potential geographic barriers, and these likely played a role in some divergence events.

Diversifying on the Ark: multiple new endemic lineages of dwarf geckos from the Western Ghats provide insights into the systematics and biogeography of South Asian Cnemaspis (Reptilia: Squamata).

The highly speciose gekkonid genus Cnemaspis Strauch, 1887 is polyphyletic, with three distantly related and geographically isolated clades from Africa, South Asia (SA), and Southeast Asia. At present, there are 85 known species within SA Cnemaspis, although the number continues to increase rapidly with focused surveys and rigorous taxonomic work. Recent studies have provided valuable insights into the diversity and evolutionary history of SA Cnemaspis; however, most of these studies lack sufficient sampling in the Western Ghats (WG), where the genus has its greatest diversity. We addressed this research gap by conducting extensive sampling across the WG and re-examining museum specimens, thus providing a systematic account of various extant Cnemaspis species along with their distribution and natural history. We described 12 new species and a southern WG endemic clade of SA Cnemaspis. Ten of the newly described species are endemic to the forests of the southern WG. We also identified 10 well-supported subclades that can be separated across morphological, geographic, and phylogenetic axes. A time-calibrated phylogeny and ancestral area reconstructions confirmed the Paleocene origin of SA Cnemaspis in the WG and provide insights into its evolutionary history and biogeography. The discovery of multiple endemic and deeply divergent lineages further highlights the evolutionary significance of the WG for lizards.

Keratin intermediate filament chains in the European common wall lizard (Podarcis muralis) and a potential keratin filament crosslinker.

On the basis of sequence homology with mammalian α-keratins, and on the criteria that the coiled-coil segments and central linker in the rod domain of these molecules must have conserved lengths if they are to assemble into viable intermediate filaments, a total of 28 Type I and Type II keratin intermediate filament chains (KIF) have been identified from the genome of the European common wall lizard (Podarcis muralis). Using the same criteria this number may be compared to 33 found here in the green anole lizard (Anole carolinensis) and 25 in the tuatara (Sphenodon punctatus). The Type I and Type II KIF genes in the wall lizard fall in clusters on chromosomes 13 and 2 respectively. Although some differences occur in the terminal domains in the KIF chains of the two lizards and tuatara, the similarities between key indicator residues - cysteine, glycine and proline - are significant. The terminal domains of the KIF chains in the wall lizard also contain sequence repeats commonly based on glycine and large apolar residues and would permit the fine tuning of physical properties when incorporated within the intermediate filaments. The H1 domain in the Type II chain is conserved across the lizards, tuatara and mammals, and has been related to its role in assembly at the 2-4 molecule level. A KIF-like chain (K80) with an extensive tail domain comprised of multiple tandem repeats has been identified as having a potential filament-crosslinking role.

Ecological Opportunity from Innovation, not Islands, Drove the Anole Lizard Adaptive Radiation.

Islands are thought to facilitate adaptive radiation by providing release from competition and predation. Anole lizards are considered a classic example of this phenomenon: different ecological specialists ("ecomorphs") evolved in the Caribbean Greater Antilles (Cuba, Hispaniola, Jamaica, and Puerto Rico), resulting in convergent assemblages that are not observed in mainland Latin America. Yet, the role of islands in facilitating adaptive radiation is more often implied than directly tested, leaving uncertain the role of biogeography in stimulating diversification. Here, we assess the proposed "island effect" on anole diversification using Bayesian phylogenetic comparative methods that explicitly incorporate rate heterogeneity across the tree and demonstrate two cases of would be false positives. We discovered that rates of speciation and morphological evolution of island and mainland anoles are equivalent, implying that islands provide no special context for exceptionally rapid diversification. Likewise, rates of evolution were equivalent between island anoles that arose via in situ versus dispersal-based mechanisms, and we found no evidence for island-specific rates of speciation or morphological evolution. Nonetheless, the origin of Anolis is characterized by a speciation pulse that slowed over time-a classic signature of waning ecological opportunity. Our findings cast doubt on the notion that islands catalyzed the anole adaptive radiation and instead point to a key innovation, adhesive toe pads, which facilitated the exploitation of many arboreal niches sparsely utilized by other iguanian lizards. The selective pressures responsible for arboreal niche diversification differ between islands and the mainland, but the tempo of diversification driven by these discordant processes is indistinguishable. [Anolis; Caribbean; key innovation; morphological evolution; RevBayes; speciation.].

A new species of Asian gracile skink (Scincidae: Lygosominae: emSubdoluseps/em) from the rain-shadow belts of Nilgiri hills, Western Ghats, India.

We describe a new species of Asian gracile skink from the dry leeward slopes of the Nilgiri hills, Tamil Nadu state, India which forms a part of the eastern, rain shadow escarpment of the Western Ghats in peninsular India. The new species, Subdoluseps nilgiriensis sp. nov., is characterized by: slender, small-sized body (47-67 mm); sandy brown above, with each scale tipped with black; a thick black lateral band from snout to tail; a distinct white labial streak; dirty white venter, with throat having mild black striations; 28-29 midbody scale rows; 71-74 mid ventral scales; 66-69 paravertebral scales. The new species is described based on external morphological characters, genetic data and geographical isolation. Based on two mitochondrial DNA genes, we show that the new species shares a sister relationship with Subdoluseps pruthi (Sharma, 1977) which is found in parts of the Eastern Ghats in peninsular India. The discovery of this new population raises two novel scenarios. Firstly, it renders the genus Subdoluseps evolutionarily polyphyletic with respect to the Indian species included in this genus. Secondly, it falsifies the notion that S. pruthi group skinks are restricted to the Eastern Ghats. Our results further indicate that the dry zone of peninsular India has unrealized skink diversity that needs to be further explored.

A taxonomic conundrum: Characterizing a cryptic radiation of Asian gracile skinks (Squamata: Scincidae: Riopa) in Myanmar.

Recognizing species-level diversity is important for studying evolutionary patterns across biological disciplines and is critical for conservation efforts. However, challenges remain in delimiting species-level diversity, especially in cryptic radiations where species are genetically divergent but show little morphological differentiation. Using multilocus molecular data, phylogenetic analyses, species delimitation analyses, and morphological data, we examine lineage diversification in a cryptic radiation of Riopa skinks in Myanmar. Four species of Riopa skinks are currently recognized from Myanmar based on morphological traits, but the boundaries between three of these species, R. anguina, R. lineolata, and R. popae, are not well-defined. We find high levels of genetic diversity within these three species. Our analyses suggest that they may comprise as many as 12 independently evolving lineages, highlighting the extent to which species diversity in the region is underestimated. However, quantitative trait data suggest that these lineages have not differentiated morphologically, possibly indicating that this cryptic radiation represents non-adaptive evolution, although additional data is needed to corroborate this.

Cryptic diversity and non-adaptive radiation of montane New Guinea skinks (Papuascincus; Scincidae).

New Guinea, the world's largest and highest tropical island, has a rich but poorly known biota. Papuascincus is a genus of skinks endemic to New Guinea's mountain regions, comprising two wide-ranging species and two species known only from their type series. The phylogeny of the genus has never been examined and the relationships among its species - as well as between it and closely related taxa - are hitherto unknown. We performed the first large-scale molecular-phylogenetic study of Papuascincus, including sampling across the genus' range in Papua New Guinea. We sequenced three mitochondrial and two nuclear markers from 65 specimens of Papuascincus and reconstructed their phylogenetic relationships. We also performed species-delimitation analyses, estimated divergence times and ancestral biogeography, and examined body-size evolution within the genus. Papuascincus was strongly supported as monophyletic. It began radiating during the mid-Miocene in the area now comprising the Central Cordillera of New Guinea, then dispersed eastward colonising the Papuan Peninsula. We found evidence of extensive cryptic diversity within the genus, with between nine and 20 supported genetic lineages. These were estimated using three methods of species delimitation and predominantly occur in allopatry. Distribution and body-size divergence patterns indicated that character displacement in size took place during the evolutionary history of Papuascincus. We conclude that the genus requires comprehensive taxonomic revision and likely represents a species-rich lineage of montane skinks.

Parachute geckos free fall into synonymy: Gekko phylogeny, and a new subgeneric classification, inferred from thousands of ultraconserved elements.

Recent phylogenetic studies of gekkonid lizards have revealed unexpected, widespread paraphyly and polyphyly among genera, unclear generic boundaries, and a tendency towards the nesting of taxa exhibiting specialized, apomorphic morphologies within geographically widespread "generalist" clades. This is especially true in Australasia, where monophyly of Gekko proper has been questioned with respect to phenotypically ornate flap-legged geckos of the genus Luperosaurus, the Philippine false geckos of the genus Pseudogekko, and even the elaborately "derived" parachute geckos of the genus Ptychozoon. Here we employ sequence capture targeting 5060 ultraconserved elements (UCEs) to infer phylogenomic relationships among 42 representative ingroup gekkonine lizard taxa. We analyze multiple datasets of varying degrees of completeness (10, 50, 75, 95, and 100 percent complete with 4715, 4051, 3376, 2366, and 772 UCEs, respectively) using concatenated maximum likelihood and multispecies coalescent methods. Our sampling scheme addresses four persistent systematic questions in this group: (1) Are Luperosaurus and Ptychozoon monophyletic, and are any of these named species truly nested within Gekko? (2) Are prior phylogenetic estimates of Sulawesi's L. iskandari as the sister taxon to Melanesian G. vittatus supported by our genome-scale dataset? (3) Is the high-elevation L. gulat of Palawan Island correctly placed within Gekko? (4) And, finally, where do the enigmatic taxa P. rhacophorus and L. browni fall in a higher-level gekkonid phylogeny? We resolve these issues; confirm with strong support some previously inferred findings (placement of Ptychozoon taxa within Gekko; the sister taxon relationship between L. iskandari and G. vittatus); resolve the systematic position of unplaced taxa (L. gulat, and L. browni); and transfer L. iskandari, L. gulat, L. browni, and all members of the genus Ptychozoon to the genus Gekko. Our unexpected and novel systematic inference of the placement of Ptychozoon rhacophorus suggests that this species is not grouped with Ptychozoon or even Luperosaurus (as previously expected) but may, in fact, be most closely related to several Indochinese species of Gekko. With our resolved and strongly supported phylogeny, we present a new classification emphasizing the most inclusive, original generic name (Gekko) for these ~60 taxa, arranged into seven subgenera.

Repeated evolution of terrestrial lineages in a continental lizard radiation.

The "early-burst" model of adaptive radiation predicts an early increase in phenotypic disparity concurrent with lineage diversification. Although most studies report a lack of this coupled pattern, the underlying processes are not identified. The continental radiation of Hemidactylus geckos from Peninsular India includes morphologically diverse species that occupy various microhabitats. This radiation began diversifying ~36 Mya with an early increase in lineage diversification. Here, we test the "early-burst" hypothesis by investigating the presence of ecomorphs and examining the pattern of morphological diversification in a phylogenetic framework. Two ecomorphs-terrestrial and scansorial species-that vary significantly in body size and toepad size were identified. Unlike the prediction of the "early-burst" model, we find that disparity in toepad morphology accumulated more recently ~14 Mya and fit the Ornstein-Ulhenbeck model. Ancestral state reconstruction of the two ecomorphs demonstrates that terrestrial lineages evolved independently at least five times from scansorial ancestors, with the earliest diversification in terrestrial lineages 19-12 Mya. Our study demonstrates a delayed increase in morphological disparity as a result of the evolution of terrestrial ecomorphs. The diversification of terrestrial lineages is concurrent with the establishment of open habitat and grasslands in Peninsular India, suggesting that the appearance of this novel resource led to the adaptive diversification.

Evolution of cranial shape in a continental-scale evolutionary radiation of Australian lizards.

A defining character of adaptive radiations is the evolution of a diversity of morphological forms that are associated with the use of different habitats, following the invasion of vacant niches. Island adaptive radiations have been thoroughly investigated but continental scale radiations are more poorly understood. Here, we use 52 species of Australian agamid lizards and their Asian relatives as a model group, and employ three-dimensional geometric morphometrics to characterize cranial morphology and investigate whether variation in cranial shape reflects patterns expected from the ecological process of adaptive radiation. Phylogenetic affinity, evolutionary allometry, and ecological life habit all play major roles in the evolution of cranial shape in the sampled lizards. We find a significant association between cranial shapes and life habit. Our results are in line with the expectations of an adaptive radiation, and this is the first time detailed geometric morphometric analyses have been used to understand the selective forces that drove an adaptive radiation at a continental scale.

Male Meiotic Recombination in the Steppe Agama, Trapelus sanguinolentus (Agamidae, Iguania, Reptilia).

Meiotic recombination rates and patterns of crossover distributions along the chromosomes vary considerably even between closely related species. The adaptive significance of these differences is still unclear due to the paucity of empirical data. Most data on recombination come from mammalian species, while other vertebrate clades are poorly explored. Using immunolocalization of the protein of the lateral element of the synaptonemal complex (SYCP3) and the mismatch-repair protein MLH1, which marks mature recombination nodules, we analyzed recombination rates and crossover distribution in meiotic prophase chromosomes of the steppe agama (Trapelus sanguinolentus, Agamidae, Acrodonta, Iguania) and compared them with data obtained for the genus Anolis (Dactyloidae, Pleurodonta, Iguania). We found that, despite a smaller genome size, the total SC length and the MLH1 focus number per cell are much higher in the agama than in the anoles. The distributions of the MLH1 foci in the agama are multimodal in larger chromosomes and bimodal in smaller chromosomes without a significant centromere effect, resembling the patterns known for birds. A possible relationship between karyotype remodeling and the evolution of recombination in Iguania is discussed.

Does Population Structure Predict the Rate of Speciation? A Comparative Test across Australia's Most Diverse Vertebrate Radiation.

Population divergence is the first step in allopatric speciation, as has long been recognized in both theoretical models of speciation and empirical explorations of natural systems. All else being equal, lineages with substantial population differentiation should form new species more quickly than lineages that maintain range-wide genetic cohesion through high levels of gene flow. However, there have been few direct tests of the extent to which population differentiation predicts speciation rates as measured on phylogenetic trees. Here, we explicitly test the links between organismal traits, population-level processes, and phylogenetic speciation rates across a diverse clade of Australian lizards that shows remarkable variation in speciation rate. Using genome-wide double digest restriction site-associated DNA data from 892 individuals, we generated a comparative data set on isolation by distance and population differentiation across 104 putative species-level lineages (operational taxonomic units). We find that species show substantial variation in the extent of population differentiation, and this variation is predicted by organismal traits that are thought to be proxies for dispersal and deme size. However, variation in population structure does not predict variation in speciation rate. Our results suggest that population differentiation is not the rate-limiting step in species formation and that other ecological and historical factors are primary determinants of speciation rates at macroevolutionary scales.

Integrative species delimitation in practice: Revealing cryptic lineages within the short-nosed skink Plestiodon brevirostris (Squamata: Scincidae).

Integrative taxonomy has been generally considered as a goal in systematics for more than a decade. Here, we employed environmental, molecular, and morphological data to evaluate the species boundaries within the short-nosed skink Plestiodon brevirostris from south-central Mexico, one member of the morphologically conservative P. brevirostris group. Our molecular dataset includes one mitochondrial and two nuclear loci. The mitochondrial fragment includes the full length of the gene coding for the NADH dehydrogenase subunit 1 protein, a segment of the gene coding for 16S ribosomal RNA, and flanking tRNAs. The nuclear dataset includes fragments of the genes coding for the megakaryoblastic leukemia 1 and RNA fingerprint 35 proteins. We employed phylogenetic reconstruction, analyses of population structure and morphological variation, and species delimitation methods (including the integration of the three kinds of data in a unified probabilistic framework) to evaluate species limits. Our results suggest that P. brevirostris represents four distinct species. The information provided by each kind of data allowed us to discern between alternative explanations for the observed patterns of geographic structure. Two of the newly recognized lineages are poorly differentiated morphologically but apparently differ in environmental preferences and are allopatric. Additionally, one lineage is microendemic and parapatric with respect to another one. Moreover, our phylogenetic analyses suggest that other taxa within the P. brevirostris group may represent species complexes. We discuss our results in the context of integrative species delimitation.

Phylogenomics of a rapid radiation: the Australian rainbow skinks.

BACKGROUND: The application of target capture with next-generation sequencing now enables phylogenomic analyses of rapidly radiating clades of species. But such analyses are complicated by extensive incomplete lineage sorting, demanding the use of methods that consider this process explicitly, such as the multispecies coalescent (MSC) model. However, the MSC makes strong assumptions about divergence history and population structure, and when using the full Bayesian implementation, current computational limits mean that relatively few loci and samples can be analysed for even modest sized radiations. We explore these issues through analyses of an extensive (> 1000 loci) dataset for the Australian rainbow skinks. This group consists of 3 genera and 41 described species, which likely diversified rapidly in Australia during the mid-late Miocene to occupy rainforest, woodland, and rocky habitats with corresponding diversity of morphology and breeding colouration. Previous phylogenetic analyses of this group have revealed short inter-nodes and high discordance among loci, limiting the resolution of inferred trees. A further complication is that many species have deep phylogeographic structure - this poses the question of how to sample individuals within species for analyses using the MSC. RESULTS: Phylogenies obtained using concatenation and summary coalescent species tree approaches to the full dataset are well resolved with generally consistent topology, including for previously intractable relationships near the base of the clade. As expected, branch lengths at the tips are substantially overestimated using concatenation. Comparisons of different strategies for sampling haplotypes for full Bayesian MSC analyses (for one clade and using smaller sets of loci) revealed, unexpectedly, that combining haplotypes across divergent phylogeographic lineages yielded consistent species trees. CONCLUSIONS: This study of more than 1000 loci provides a strongly-supported estimate of the phylogeny of the Australian rainbow skinks, which will inform future research on the evolution and taxonomy of this group. Our analyses suggest that species tree estimation with the MSC can be quite robust to violation of the assumption that the individuals representing a taxon are sampled from a panmictic population.

Origins and biogeography of the Anolis crassulus subgroup (Squamata: Dactyloidae) in the highlands of Nuclear Central America.

BACKGROUND: Recent studies have begun to reveal the complex evolutionary and biogeographic histories of mainland anoles in Central America, but the origins and relationships of many taxa remain poorly understood. One such group is the Anolis (Norops) crassulus species subgroup, which contains ten morphologically similar highland taxa, the majority of which have restricted distributions. The nominal taxon A. crassulus has a disjunct distribution from Chiapas, Mexico, through Guatemala, in the highlands of El Salvador, and in the Chortís Highlands of Honduras. We test the relationships of these species using multiple mitochondrial and nuclear loci in concatenated and multispecies coalescent frameworks, in an effort to both resolve long-standing taxonomic confusion and present new insights into the evolution and biogeography of these taxa. RESULTS: Sequences of multiple mitochondrial and nuclear loci were generated for eight of the ten species of the Anolis crassulus species subgroup. We analyzed phylogenetic relationships and estimated divergence times and ancestral ranges of the subgroup, recovering a monophyletic subgroup within Anolis. Within the nominal taxon Anolis crassulus, we recovered multiple genetically distinct lineages corresponding to allopatric populations, and show that the Chortís Highland lineage split from the others over 13 MYA. Additionally, distinct mitochondrial lineages are present within the taxa A. heteropholidotus and A. morazani, and importantly, samples of A. crassulus and A. sminthus previously used in major anole phylogenetic analyses are not recovered as conspecific with those taxa. We infer a Chortís Highland origin for the ancestor of this subgroup, and estimate cladogenesis of this subgroup began approximately 22 MYA. CONCLUSIONS: Our results provide new insights into the evolution, biogeography, and timing of diversification of the Anolis crassulus species subgroup. The disjunctly distributed Anolis crassulus sensu lato represents several morphologically conserved, molecularly distinct anoles, and several other species in the subgroup contain multiple isolated lineages.

Cryptic diversity in Ptyodactylus (Reptilia: Gekkonidae) from the northern Hajar Mountains of Oman and the United Arab Emirates uncovered by an integrative taxonomic approach.

The Hajar Mountains of south-eastern Arabia form an isolated massif surrounded by the sea to the east and by a large desert to the west. As a result of their old geological origin, geographical isolation, complex topography and local climate, these mountains provide an important refuge for endemic and relict species of plants and animals. With 19 species restricted to the Hajar Mountains, reptiles are the vertebrate group with the highest level of endemicity, becoming an excellent model for understanding the patterns and processes that generate and shape diversity in this arid mountain range. The geckos of the Ptyodactylus hasselquistii species complex are the largest geckos in Arabia and are found widely distributed across the Arabian Mountains, constituting a very important component of the reptile mountain fauna. Preliminary analyses suggested that their diversity in the Hajar Mountains may be higher than expected and that their systematics should be revised. In order to tackle these questions, we inferred a nearly complete calibrated phylogeny of the genus Ptyodactylus to identify the origin of the Hajar Mountains lineages using information from two mitochondrial and four nuclear genes. Genetic variability within the Hajar Mountains was further investigated using 68 specimens of Ptyodactylus from 46 localities distributed across the entire mountain range and sequenced for the same genes as above. The molecular phylogenies and morphological analyses as well as niche comparisons indicate the presence of two very old sister cryptic species living in allopatry: one restricted to the extreme northern Hajar Mountains and described as a new species herein; the other distributed across the rest of the Hajar Mountains that can be confidently assigned to the species P. orlovi. Similar to recent findings in the geckos of the genus Asaccus, the results of the present study uncover more hidden diversity in the northern Hajar Mountains and stress once again the importance of this unique mountain range as a hot spot of biodiversity and a priority focal point for reptile conservation in Arabia.

Identification and characterization of tyrosine kinases in anole lizard indicate the conserved tyrosine kinase repertoire in vertebrates.

The tyrosine kinases (TKs) play principal roles in regulation of multicellular aspects of the organism and are implicated in many cancer types and congenital disorders. The anole lizard has recently been introduced as a model organism for laboratory-based studies of organismal function and field studies of ecology and evolution. However, the TK family of anole lizard has not been systematically identified and characterized yet. In this study, we identified 82 TK-encoding genes in the anole lizard genome and classified them into 28 subfamilies through phylogenetic analysis, with no member from ROS and STYK1 subfamilies identified. Although TK domain sequences and domain organization in each subfamily were conserved, the total number of TKs in different species was much variable. In addition, extensive evolutionary analysis in metazoans indicated that TK repertoire in vertebrates tends to be remarkably stable. Phylogenetic analysis of Eph subfamily indicated that the divergence of EphA and EphB occurred prior to the whole genome duplication (WGD) but after the split of Urochordates and vertebrates. Moreover, the expression pattern analysis of lizard TK genes among 9 different tissues showed that 14 TK genes exhibited tissue-specific expression and 6 TK genes were widely expressed. Comparative analysis of TK expression suggested that the tissue specifically expressed genes showed different expression pattern but the widely expressed genes showed similar pattern between anole lizard and human. These results may provide insights into the evolutionary diversification of animal TK genes and would aid future studies on TK protein regulation of key growth and developmental processes.

Molecular systematics, species delimitation and diversification patterns of the Phyllodactylus lanei complex (Gekkota: Phyllodactylidae) in Mexico.

The description of cryptic gecko species worldwide has revealed both that many putative species are, in fact, conformed by a complex of morphologically conserved species that are genetically distinct and highly divergent, and that gecko species diversity could be underestimated. The taxonomy and species delimitation of geckos belonging to the genus Phyllodactylus is still controversial, 16 of which are distributed in Mexico and 13 are endemic. Although the large morphological variation shown by the Phyllodactylus species from Mexico has been amply documented, little is known about their genetic diversity and evolutionary relationships, and much less regarding cryptic speciation. Here, we included the most comprehensive sampling of populations and species of the Phyllodactylus lanei complex distributed in Mexico, and applied an analytical approach that included probabilistic phylogenetic analyses, jointly with species delimitation methods and Bayesian putative species validation analysis. Our results suggest the existence of 10 lineages within the complex, supporting the existence of cryptic species, and in great contrast with the current taxonomic proposal that includes only four subspecies. The most recent common ancestor (MRCA) for the P. lanei clade originated on the Early Eocene (∼54Mya), along the southern coasts of Mexico, followed by the highest diversification of the complex MRCA during the Eocene (34-56Mya). Lineages subsequently dispersed and diversified towards the northwest, and the diversification process ended with the most recent lineages inhabiting two islands on the coasts of Nayarit (Miocene; 5.5-23Mya). Our results highlight three vicariant events associated with the evolution of the lineages, two of them intimately related to the formation of the Sierra Madre del Sur and the Transmexican Volcanic Belt mountain ranges, main geographic barriers that isolated and facilitated the divergence and speciation in this group of geckos. Finally, we propose that there are 10 species in the P. lanei complex, from which four represent taxonomic changes and six are new species and require a formal description. We acknowledge that more analyses, including a detailed evaluation of morphological characters and use of more unlinked nuclear loci with enough variability, are needed to further support their taxonomic description.

Pre-Quaternary divergence and subsequent radiation explain longitudinal patterns of genetic and morphological variation in the striped skink, Heremites vittatus.

BACKGROUND: Many animal and plant species in the Middle East and northern Africa have a predominantly longitudinal distribution, extending from Iran and Turkey along the eastern Mediterranean coast into northern Africa. These species are potentially characterized by longitudinal patterns of biological diversity, but little is known about the underlying biogeographic mechanisms and evolutionary timescales. We examined these questions in the striped skink, Heremites vittatus, one such species with a roughly longitudinal distribution across the Middle East and northern Africa, by analyzing range-wide patterns of mitochondrial DNA (mtDNA) sequence and multi-trait morphological variation. RESULTS: The striped skink exhibits a basic longitudinal organization of mtDNA diversity, with three major mitochondrial lineages inhabiting northern Africa, the eastern Mediterranean coast, and Turkey/Iran. Remarkably, these lineages are of pre-Quaternary origin, and are characterized by p-distances of 9-10%. In addition, within each of these lineages a more recent Quaternary genetic diversification was observed, as evidenced by deep subclades and high haplotype diversity especially in the Turkish/Iranian and eastern Mediterranean lineages. Consistent with the genetic variation, our morphological analysis revealed that the majority of morphological traits show significant mean differences between specimens from northern Africa, the eastern Mediterranean coast, and Turkey/Iran, suggesting lineage-specific trait evolution. In addition, a subset of traits exhibits clinal variation along the eastern Mediterranean coast, potentially indicating selection gradients at the geographic transition from northern Africa to Anatolia. The existence of allopatric, morphologically and genetically divergent lineages suggests that Heremites vittatus might represent a complex with several taxa. CONCLUSIONS: Our work demonstrates that early divergence events in the Pliocene, likely driven by both climatic and geological factors, established the longitudinal patterns and distribution of Heremites vittatus. Subsequent radiation during the Pleistocene generated the genetic and morphological diversity observed today. Our study provides further evidence that longitudinal diversity patterns and species distributions in the Middle East and northern Africa were shaped by complex evolutionary processes, involving the region's intricate geological history, climatic oscillations, and the presence of the Sahara.

Developmental dynamics of ecomorphological convergence in a transcontinental lizard radiation.

Phenotypic convergence has confounded evolutionary biologists for centuries, explained as adaptations to shared selective pressures, or alternatively, the result of limited developmental pathways. We tested the relative roles of adaptation and constraint in generating convergent cranial morphologies across a large lizard radiation, the Lacertidae, whose members inhabit diverse environments throughout the Old World and display high amounts of homoplasy associated with ecological niche. Using 3D X-ray computed tomography, we quantified cranial shape variation associated with ontogeny, allometry, and ecology, covering all lacertid genera and one-third of species diversity. Landmark-based geometric morphometrics showed that cranial shape varied significantly among biomes, with substantial convergence among arid-dwelling lineages. Comparisons of species cranial growth trajectories between biomes revealed that allometric postdisplacement, as evidenced by decreased elevation of a constant ontogenetic slope, drives the convergent paedomorphic appearance of independent arid-dwelling forms. We hypothesize that observed heterochronic changes reflect temporal compression of ancestral life history in response to extreme environments, with associated phenotypes occurring as by-products of adaptive shifts in reproductive investment. Although allometry has long been considered a developmental constraint, our results demonstrate that allometric flexibility during early ontogeny produces convergent ecomorphologies over vast temporal and spatial scales, thus dramatically obscuring underlying phylogenetic signals.