Morphology of the Podarcis wall lizards (Squamata: Lacertidae) from the Iberian Peninsula and North Africa: patterns of variation in a putative cryptic species complex

Lack of country-wide systematic herpetology collections in Portugal jeopardizes future research and conservation

Natural History Collections (NHCs) represent the world's largest repositories of long-term biodiversity datasets. Specimen collection and voucher deposition has been the backbone of NHCs since their inception, but recent decades have seen a drastic decline in rates of growth via active collecting. Amphibians and reptiles are amongst the most threatened zoological groups on the planet and are historically underrepresented in most worldwide NHCs. As part of an ongoing project to review the Portuguese zoological collections in the country's NHCs, herpetological data from its three major museums and smaller collections was gathered and used to examine the coverage and representation of the different taxa extant in Portugal. These collections are not taxonomically, geographically, or temporally complete. Approximately 90% of the Portuguese herpetological taxa are represented in the country's NHCs, and around half of the taxa are represented by less than 50 specimens. Geographically, the collections cover less than 30% of the country's territory and almost all of the occurring taxa have less than 10% of their known distribution represented in the collections. A discussion on the implications for science of such incomplete collections and a review of the current status of Portuguese NHCs is presented.

The Intriguing Biogeographic Pattern of the Italian Wall Lizard Podarcis siculus (Squamata: Lacertidae) in the Tuscan Archipelago Reveals the Existence of a New Ancient Insular Clade

The Tuscan Archipelago is one of the most ancient and ecologically heterogeneous island systems in the Mediterranean. The biodiversity of these islands was strongly shaped by the Pliocene and Pleistocene sea regressions and transgression, resulting in different waves of colonization and isolation of species coming from the mainland. The Italian wall lizard, Podarcis siculus, is present on the following islands of the Tuscan Archipelago: Elba, Giglio, Giannutri, Capraia, Montecristo and Cerboli. The species in the area displays a relatively high morphological variability that in the past led to the description of several subspecies. In this study, both the genetic and morphological diversity of P. siculus of the Tuscan Archipelago were investigated. Specifically, the meristic characters and the dorsal pattern were analyzed, while the genetic relationships among these populations were explored with mtDNA and microsatellite nuclear markers to reconstruct the colonization history of the Archipelago. Our results converge in the identification of at least two different waves of colonization in the Archipelago: Elba, and the populations of Cerboli and Montecristo probably originate from historical introductions from mainland Tuscany, while those of Giglio and Capraia are surviving populations of an ancient lineage which colonized the Tuscan Archipelago during the Pliocene and which shares a common ancestry with the P. siculus populations of south-eastern Italy. Giannutri perhaps represents an interesting case of hybridization between the populations from mainland Tuscany and the Giglio-Capraia clade. Based on the high phenotypic and molecular distinctiveness of this ancient clade, these populations should be treated as distinct units deserving conservation and management efforts as well as further investigation to assess their taxonomic status.

Identification of morphologically cryptic species with computer vision models: wall lizards (Squamata: Lacertidae: Podarcis) as a case study

Automated image classification is a thriving field of machine learning, and various successful applications dealing with biological images have recently emerged. In this work, we address the ability of these methods to identify species that are difficult to tell apart by humans due to their morphological similarity. We focus on distinguishing species of wall lizards, namely those belonging to the Podarcis hispanicus species complex, which constitutes a well-known example of cryptic morphological variation. We address two classification experiments: (1) assignment of images of the morphologically relatively distinct P. bocagei and P. lusitanicus; and (2) distinction between the overall more cryptic nine taxa that compose this complex. We used four datasets (two image perspectives and individuals of the two sexes) and three deep-learning models to address each problem. Our results suggest a high ability of the models to identify the correct species, especially when combining predictions from different perspectives and models (accuracy of 95.9% and 97.1% for females and males, respectively, in the two-class case; and of 91.2% to 93.5% for females and males, respectively, in the nine-class case). Overall, these results establish deep-learning models as an important tool for field identification and monitoring of cryptic species complexes, alleviating the burden of expert or genetic identification.

Morphological diversification of Mediterranean anurans: the roles of evolutionary history and climate

Investigation of the ecological and evolutionary mechanisms governing the origin and diversification of species requires integrative approaches that often have to accommodate strong discordance among datasets. A common source of conflict is the combination of morphological and molecular characters with different evolutionary rates. Resolution of these discordances is crucial to assess the relative roles of different processes in generating and maintaining biodiversity. Anuran amphibians provide many examples of morphologically similar, genetically divergent lineages, posing questions about the relative roles of phylogeny and ecological factors in phenotypic evolution. We focused on three circum-Mediterranean anuran genera (Hyla, Alytes and Discoglossus), characterizing morphological and environmental disparity and comparing diversity patterns across biological levels of organization. Using a comparative phylogenetic framework, we tested how shared ancestry and climatic factors come together to shape phenotypic diversity. We found higher morphological differentiation within Hyla and Alytes than in Discoglossus. Body size and limb morphology contributed most to inter- and intraspecific morphological variation in Hyla and Alytes, but there was no strong phylogenetic signal, indicating that shared ancestry does not predict patterns of phenotypic divergence. In contrast, we uncovered a significant association between morphology and climatic descriptors, supporting the hypothesis that morphological disparity between species results from adaptive evolution.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae)

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

Status of the largest extant population of the critically endangered Aeolian lizard Podarcis raffonei (Capo Grosso, Vulcano island)

The Aeolian wall lizard Podarcis raffonei is an island endemic that survives only on three tiny islets, and on the Capo Grosso peninsula of the Vulcano island, thus is among the European vertebrates with the smallest range and one of the most threatened by extinction. This species is declining due to competition and hybridization with the non-native lizard Podarcis siculus, but a regular monitoring program is lacking. Here we assessed the size and status of the Capo Grosso population of P. raffonei on Vulcano. In September 2015 we captured 30 individuals showing the typical brown phenotype of P. raffonei, while one single male showed a green phenotype, apparently intermediate between P. raffonei and the non-native Podarcis siculus. In May 2017, only 47% of 131 individuals showed the typical brown phenotype (P. raffonei-like) and 53% showed the green phenotype (P. siculus-like). Based on N-mixture models and removal sampling the estimated size of the Capo Grosso population was of 800–1300 individuals in 2017, being similar to 2015; available data suggest that the total range of the species could be as small as 2 ha. The frequency of individuals with the typical P. raffonei phenotype dramatically dropped between two samplings with a parallel increase of individuals displaying the green phenotype. Observation on outdoor captive-bred individuals demonstrates plasticity for colouration in P. raffonei individuals from Capo Grosso, with several individuals showing the typical brown pattern in September 2017 and a green pattern in March 2021. Non-exclusive hypotheses, including hybridization with P. siculus and plasticity in colour pattern of P. raffonei, are discussed to explain the phenotypic shifts of the P. raffonei population of Capo Grosso. While genomic evidence is required to reach conclusions and investigate eventual hybridization, it is urgent to undertake a programme for the monitoring and management of this lizard.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.

Local Segregation of Realised Niches in Lizards

Species can occupy different realised niches when sharing the space with other congeneric species or when living in allopatry. Ecological niche models are powerful tools to analyse species niches and their changes over time and space. Analysing how species’ realised niches shift is paramount in ecology. Here, we examine the ecological realised niche of three species of wall lizards in six study areas: three areas where each species occurs alone; and three areas where they occur together in pairs. We compared the species’ realised niches and how they vary depending on species’ coexistence, by quantifying niche overlap between pairs of species or populations with the R package ecospat. For this, we considered three environmental variables (temperature, humidity, and wind speed) recorded at each lizard re-sighting location. Realised niches were very similar when comparing syntopic species occurring in the same study area. However, realised niches differed when comparing conspecific populations across areas. In each of the three areas of syntopy, the less abundant species shift its realised niche. Our study demonstrates that sympatry may shift species’ realised niche.

Dehydration constrains thermoregulation and space use in lizards

Climate change is negatively affecting many species. The increase in mean air temperature is often associated with shifts in distribution, changes in phenology, and local extinctions. Other factors that only partially correlate with air temperature, like water shortage, may also contribute to the negative consequences of climate change. Although the effect of temperature on lizards' ecophysiology is highly studied, many lizards are also at risks of increased water loss and dehydration, which are predicted to increase under climate change. Here we aimed for the first time to explore if lacertid lizards exposed to dehydration thermoregulate less precisely than hydrated lizards and if dehydrated lizards are less active, change the daily pattern of thermoregulation and balance water balance against thermoregulation. We exposed four lizard species with differences in the thermal preference to thermal gradients with or without a source of water. We measured preferred body temperatures, daily pattern of thermoregulation, and the use of space. Dehydration negatively affected thermoregulation in all investigated species. Dehydrated lizards reduced their preferred body temperature and showed a species-specific pattern of hourly change in thermal preference. Furthermore, they more frequently used the colder parts of the gradients and spent more time hidden. Lizards experiencing dehydration may suffer a reduction in survival and fitness because of poor thermoregulation. Similarly, they may spend more time hidden, waiting for more favourable weather conditions. Such inactivity may carry ecological costs especially in those regions that undergo either short or prolonged periods of droughts.

Seasonal and interpopulational phenotypic variation in morphology and sexual signals of Podarcis liolepis lizards

Widespread species often show extensive phenotypic variation due to the contrasting abiotic and biotic factors that shape selective pressures in different environments. In this context, the gradual and predictable patterns of variation in climatic and environmental conditions found in mountain areas offer a great opportunity to explore intraspecific phenotypic variation. For instance, temperature is negatively correlated with altitude and virtually all aspects of the behavior and physiology of ectotherms are sensitive to body temperature. In this work, we tested the hypothesis that morphology, dorsal and ventral coloration and the chemical profile of femoral secretions show interpopulational and seasonal variation in the lacertid lizard (Podarcis liolepis). We compared lizards from three populations inhabiting lowland and highland habitats in the French Pyrenees that were closely related genetically. We found that highland lizards were larger, stockier, had longer heads and more femoral pores and had a darker dorsal coloration than lowland ones. In addition, we detected interpopulational differences in both the abundance and the richness of chemical compounds in the glandular secretions, and we also found seasonal variation in the overall chemical composition. Dorsal and ventral coloration differed seasonally and between populations. Ventral and dorsal brightness were higher in lowland than in highland lizards in the reproductive season whereas the reversed trend was found in the non-reproductive season but only for dorsal brightness. In addition, all lizards had browner dorsal coloration in the non-reproductive season, and lowland lizards were greener in the reproductive season. By integrating information from both visual and chemical systems, our works offers a comprehensive view of how these lizards communicate in a multimodal context.

Where does diversity come from? Linking geographical patterns of morphological, genetic, and environmental variation in wall lizards

BACKGROUND

Understanding how phenotypic variation scales from individuals, through populations, up to species, and how it relates to genetic and environmental factors, is essential for deciphering the evolutionary mechanisms that drive biodiversity. We used two species of Podarcis wall lizards to test whether phenotypic diversity within and divergence across populations follow concordant patterns, and to examine how phenotypic variation responds to genetic and environmental variability across different hierarchical levels of biological organization, in an explicit geographic framework.

RESULTS

We found a general concordance of phenotypic variation across hierarchical levels (i.e. individuals and populations). However, we also found that within-population diversity does not exhibit a coherent geographic structure for most traits, while among-population divergence does, suggesting that different mechanisms may underlie the generation of diversity at these two levels. Furthermore, the association of phenotypic variation with genetic and environmental factors varied extensively between hierarchical levels and across traits, hampering the identification of simple rules to explain what yields diversity.

CONCLUSIONS

Our results in some cases comply with general ecological and evolutionary predictions, but in others they are difficult to explain in the geographic framework used, suggesting that habitat characteristics and other regulatory mechanisms may have a more substantial contribution in shaping phenotypic diversity.

Molecular phylogeny and timing of diversification in Alpine Rhithrogena (Ephemeroptera: Heptageniidae)

BACKGROUND

Larvae of the Holarctic mayfly genus Rhithrogena Eaton, 1881 (Ephemeroptera, Heptageniidae) are a diverse and abundant member of stream and river communities and are routinely used as bio-indicators of water quality. Rhithrogena is well diversified in the European Alps, with a number of locally endemic species, and several cryptic species have been recently detected. While several informal species groups are morphologically well defined, a lack of reliable characters for species identification considerably hampers their study. Their relationships, origin, timing of speciation and mechanisms promoting their diversification in the Alps are unknown.

RESULTS

Here we present a species-level phylogeny of Rhithrogena in Europe using two mitochondrial and three nuclear gene regions. To improve sampling in a genus with many cryptic species, individuals were selected for analysis according to a recent DNA-based taxonomy rather than traditional nomenclature. A coalescent-based species tree and a reconstruction based on a supermatrix approach supported five of the species groups as monophyletic. A molecular clock, mapped on the most resolved phylogeny and calibrated using published mitochondrial evolution rates for insects, suggested an origin of Alpine Rhithrogena in the Oligocene/Miocene boundary. A diversification analysis that included simulation of missing species indicated a constant speciation rate over time, rather than any pronounced periods of rapid speciation. Ancestral state reconstructions provided evidence for downstream diversification in at least two species groups.

CONCLUSIONS

Our species-level analyses of five gene regions provide clearer definitions of species groups within European Rhithrogena. A constant speciation rate over time suggests that the paleoclimatic fluctuations, including the Pleistocene glaciations, did not significantly influence the tempo of diversification of Alpine species. A downstream diversification trend in the hybrida and alpestris species groups supports a previously proposed headwater origin hypothesis for aquatic insects.

Effects of environmental disturbance on phenotypic variation: an integrated assessment of canalization, developmental stability, modularity, and allometry in lizard head shape

When populations experience suboptimal conditions, the mechanisms involved in the regulation of phenotypic variation can be challenged, resulting in increased phenotypic variance. This kind of disturbance can be diagnosed by using morphometric tools to study morphological patterns at different hierarchical levels and evaluate canalization, developmental stability, integration, modularity, and allometry. We assess the effect of urbanization on phenotypic variation in the common wall lizard (Podarcis muralis) by using geometric morphometrics to assess disturbance to head shape development. The head shapes of urban lizards were more variable and less symmetric, suggesting that urban living is more likely to disturb development. Head shape variation was congruent within and across individuals, which indicated that canalization and developmental stability are two related phenomena in these organisms. Furthermore, urban lizards exhibited smaller mean head sizes, divergent size-shape allometries, and increased deviation from within-group allometric lines. This suggests that mechanisms regulating head shape allometry may also be disrupted. The integrated evaluation of several measures of developmental instability at different hierarchical levels, which provided in this case congruent results, can be a powerful methodological guide for future studies, as it enhances the detection of environmental disturbances on phenotypic variation and aids biological interpretation of the results.

Cryptic diversity within the Anatololacerta species complex (Squamata: Lacertidae) in the Anatolian Peninsula: evidence from a multi-locus approach

The rapid development of innovative molecular tools for characterizing biodiversity is leading to an extensive and sometimes unexpected renovation of taxonomic classifications. Particularly, for species having allopatric or parapatric distributions or resulting from recent speciation processes, the absence of clear phenotypic differentiation may hinder the recognition of closely related taxa, while intraspecific polymorphism may be confused with the presence of more than one single species. In the present work, we apply different phylogenetic methods in order to infer relationships within the genus Anatololacerta, and to assess the taxonomy of this morphologically diversified group of lizards endemic to western and southern Anatolia and some neighboring Aegean islands. According to morphology, three species have been recognized (Anatololacerta anatolica, A. oertzeni and A. danfordi) as well as several subspecies, but small variation at immunological markers led some authors to join all the populations into one single taxon, A. danfordi. By selecting both mitochondrial and nuclear informative markers, we tested the effectiveness of classical "gene tree" (i.e. Bayesian Inference) vs. innovative (i.e. coalescent-based) "species tree" methods in resolving the Anatololacerta taxonomic enigma, as a case in point for similar studies on species complexes resulting from non-obvious and cryptic diversification patterns. According to our results, the gene tree method failed in resolving phylogenetic relationships among clades, whereas the multi-locus species tree approach, coupled with species delimitation methods, allowed the identification of four well distinct species. These species probably diversified in different allopatric refugia located in southern and western Anatolia, where isolated populations may have persisted during Pleistocene glacial cycles.

A comparison of multiple methods for estimating parasitemia of hemogregarine hemoparasites (apicomplexa: adeleorina) and its application for studying infection in natural populations

Identifying factors influencing infection patterns among hosts is critical for our understanding of the evolution and impact of parasitism in natural populations. However, the correct estimation of infection parameters depends on the performance of detection and quantification methods. In this study, we designed a quantitative PCR (qPCR) assay targeting the 18 S rRNA gene to estimate prevalence and intensity of Hepatozoon infection and compared its performance with microscopy and PCR. Using qPCR, we also compared various protocols that differ in the biological source and the extraction methods. Our results show that the qPCR approach on DNA extracted from blood samples, regardless of the extraction protocol, provided the most sensitive estimates of Hepatozoon infection parameters; while allowed us to differentiate between mixed infections of Adeleorinid (Hepatozoon) and Eimeriorinid (Schellackia and Lankesterella), based on the analysis of melting curves. We also show that tissue and saline methods can be used as low-cost alternatives in parasitological studies. The next step was to test our qPCR assay in a biological context, and for this purpose we investigated infection patterns between two sympatric lacertid species, which are naturally infected with apicomplexan hemoparasites, such as the genera Schellackia (Eimeriorina) and Hepatozoon (Adeleorina). From a biological standpoint, we found a positive correlation between Hepatozoon intensity of infection and host body size within each host species, being significantly higher in males, and higher in the smaller sized host species. These variations can be associated with a number of host intrinsic factors, like hormonal and immunological traits, that require further investigation. Our findings are relevant as they pinpoint the importance of accounting for methodological issues to better estimate infection in parasitological studies, and illustrate how between-host factors can influence parasite distributions in sympatric natural populations.

Lizards from urban areas are more asymmetric: using fluctuating asymmetry to evaluate environmental disturbance

The increase in human activities that leads to wildlife decline and species extinction poses an urgent need for simple indicators of environmental stress in animal populations. Several studies have suggested that fluctuating asymmetry (FA) can be an easy, direct measure of developmental instability because it is associated to environmental stress and, as such, it can be a useful indicator of population disturbance. We examined three different morphological traits in urban and rural populations of the common wall lizard (Podarcis muralis) to test whether anthropogenic disturbance causes an increase in FA. Compared to rural populations, urban ones showed higher levels of FA in all analyzed traits, thus providing evidence that FA can respond to anthropogenic disturbance. However, we also found significant differences in FA among traits, where femoral pores and subdigital lamellae, traits with a functional relevance, were more stable developmentally compared to supracilliar granules which have no evident function. Unsigned FA [abs(right-left)] exhibited significant, but weak, positive correlations among traits, indicating that developmental noise does not have a uniform effect across characters and thus questioning the view of developmental stability as an organism-wide property. The degree of signed FA (right-left) was more similar between structurally associated traits, possibly as an outcome of morphological integration. In conclusion, our results demonstrate that FA can be a reliable indicator of disturbance provided that it is analyzed on multiple traits simultaneously and examined at the population level.

Cryptic speciation patterns in Iranian rock lizards uncovered by integrative taxonomy

While traditionally species recognition has been based solely on morphological differences either typological or quantitative, several newly developed methods can be used for a more objective and integrative approach on species delimitation. This may be especially relevant when dealing with cryptic species or species complexes, where high overall resemblance between species is coupled with comparatively high morphological variation within populations. Rock lizards, genus Darevskia, are such an example, as many of its members offer few diagnostic morphological features. Herein, we use a combination of genetic, morphological and ecological criteria to delimit cryptic species within two species complexes, D. chlorogaster and D. defilippii, both distributed in northern Iran. Our analyses are based on molecular information from two nuclear and two mitochondrial genes, morphological data (15 morphometric, 16 meristic and four categorical characters) and eleven newly calculated spatial environmental predictors. The phylogeny inferred for Darevskia confirmed monophyly of each species complex, with each of them comprising several highly divergent clades, especially when compared to other congeners. We identified seven candidate species within each complex, of which three and four species were supported by Bayesian species delimitation within D. chlorogaster and D. defilippii, respectively. Trained with genetically determined clades, Ecological Niche Modeling provided additional support for these cryptic species. Especially those within the D. defilippii-complex exhibit well-differentiated niches. Due to overall morphological resemblance, in a first approach PCA with mixed variables only showed the separation between the two complexes. However, MANCOVA and subsequent Discriminant Analysis performed separately for both complexes allowed for distinction of the species when sample size was large enough, namely within the D. chlorogaster-complex. In conclusion, the results support four new species, which are described herein.