Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region.

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and âˆ¼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

Hybrid zone analysis confirms cryptic species of banded newt and does not support competitive displacement since secondary contact.

When two putatively cryptic species meet in nature, hybrid zone analysis can be used to estimate the extent of gene flow between them. Two recently recognized cryptic species of banded newt (genus Ommatotriton) are suspected to meet in parapatry in Anatolia, but a formal hybrid zone analysis has never been conducted. We sample populations throughout the range, with a focus on the supposed contact zone, and genotype them for 31 nuclear DNA SNP markers and mtDNA. We determine the degree of genetic admixture, introgression, and niche overlap. We reveal an extremely narrow hybrid zone, suggesting strong selection against hybrids, in line with species status. The hybrid zone does not appear to be positioned at an ecological barrier, and there is significant niche overlap. Therefore, the hybrid zone is best classified as a tension zone, maintained by intrinsic selection against hybrids. While the two banded newt species can evidently backcross, we see negligible introgression and the pattern is symmetric, which we interpret as supporting the fact that the hybrid zone has been practically stationary since its origin (while extensive hybrid zone movement has been suggested in other newt genera in the region). Our study illustrates the use of hybrid zone analysis to test cryptic species status.

Phenotypic variation in Heremites vittatus (Olivier, 1804) (Sauria: Scincidae) from Iran and Turkey.

According to a large morphological dataset of specimens from Turkey to Iran and based on several morphological analyses, the Iranian populations of the skink Heremites vittatus are separated from other populations of this taxon in Turkey. The values of most of morphological characters were higher in the Turkish populations. Morphological variation among populations of H. vittatus (Olivier, 1804) from Turkey and the western slopes of the Zagros Mountains in Kermanshah province in Iran may be the result of different dispersal and vicariance events. Comparison the current study dataset with specimens from Egypt can definite the taxonomic status of Iranian and Turkish populations.

Body size miniaturization in a lineage of colubrid snakes: Implications for cranial anatomy.

As body size strongly determines the biology of an organism at all levels, it can be expected that miniaturization comes with substantial structural and functional constraints. Dwarf snakes of the genus Eirenis are derived from big, surface-dwelling ancestors, considered to be similar to those of the sister genus Dolichophis. To better understand the structural implications of miniaturization on the feeding apparatus in Eirenis, the morphology of the cranial musculoskeletal system of Dolichophis schmidti was compared with that of the miniature Eirenis punctatolineatus and E. persicus using high-resolution µCT data. The gape index was compared between D. schmidti and 14 Eirenis species. Our results show a relatively increased neurocranium size and decreased maximal jaw muscle force in E. persicus, compared with the D. schmidti, and an intermediate situation in E. punctatolineatus. A significant negative allometry in gape index relative to body size is observed across the transition from the Dolichophis to Pediophis and Eirenis subgenera. However, the gape index relative to head size showed a significant negative allometry only across the transition from the Dolichophis to Pseudocyclophis subgenus. In Dolichophis-Eirenis dwarfing lineages, different structural patterns are observed through miniaturization, indicating that overcoming the challenge of miniaturization has achieved via different adaptations.

Slow life-history strategies are associated with negligible actuarial senescence in western Palaearctic salamanders.

Actuarial senescence has been viewed for a long time as an inevitable and uniform process. However, the work on senescence has mainly focused on endotherms with deterministic growth and low regeneration capacity during the adult stage, leading to a strong taxonomic bias in the study of ageing. Recent studies have highlighted that senescence could indeed display highly variable trajectories that correlate with species life-history traits. Slow life histories and indeterminate growth seem to be associated with weak and late senescence. Furthermore, high regenerative abilities could lead to negligible senescence in ectotherms. However, demographic data for species that would allow testing of these hypotheses are scarce. Here, we investigated senescence patterns in 'true salamanders' from the western Palaearctic. Our results showed that salamanders have slow life histories and that they experience negligible senescence. This pattern was consistent at both intra- and interspecific levels, suggesting that the absence of senescence may be a phylogenetically conserved trait. The regenerative capacities of salamanders, in combination with other physiological and developmental features such as an indeterminate growth and a low metabolic rate, probably explain why these small ectotherms have lifespans similar to that of large endotherms and, in contrast with most amniotes, undergo negligible senescence. Our study seriously challenges the idea that senescence is a ubiquitous phenomenon in the tree of life.

The biogeography of Elaphe sauromates (Pallas, 1814), with a description of a new rat snake species.

BACKGROUND: The rat snake genus Elaphe once comprised several dozens of species distributed in temperate through tropical zones of the New and Old World. Based on molecular-genetic analyses in early 2000s, the genus was split into several separate genera, leaving only 15 Palearctic and Oriental species as its members. One of the three species also occurring in Europe is Elaphe sauromates, a robust snake from the Balkans, Anatolia, Caucasus, Ponto-Caspian steppes, and Levant that has been suspected to be composed of two or more genetically diverse populations. Here, we studied the genetic structure and morphological variation of E. sauromates, aiming to better understand its inter-population relationships and biogeography, and subsequently revise its taxonomy. METHODS: We reconstructed the phylogeography and analyzed the genetic structure of E. sauromates populations originating from most of its geographic range using both mitochondrial (COI, ND4) and nuclear (C-MOS, MC1R, PRLR, RAG1) DNA gene fragments. We employed Maximum likelihood and Bayesian inference methods for the phylogenetic tree reconstructions, supplemented with species delimitation methods, analysis of haplotype networks, and calculation of uncorrected p-distances. Morphological variation in 15 metric and 18 meristic characters was studied using parametric univariate tests as well as multivariate general linearized models. In total, we analyzed sequences originating from 63 specimens and morphological data from 95 specimens of E. sauromates sensu lato. RESULTS: The molecular phylogeny identified two clearly divergent sister lineages within E. sauromates, with both forming a lineage sister to E. quatuorlineata. The genetic distance between them (5.80-8.24% in mtDNA) is similar to the distances among several other species of the genus Elaphe. Both lineages are also moderately morphologically differentiated and, while none of the characters are exclusively diagnostic, their combination can be used for confident lineage identification. Here, following the criteria of genetic and evolutionary species concepts, we describe the lineage from eastern Anatolia and parts of the Lesser and Great Caucasus as a new species E. urartica sp. nov. DISCUSSION: Elaphe urartica sp. nov. represents a cryptic species whose ancestors presumably diverged from their common ancestor with E. sauromates around the Miocene-Pliocene boundary. The intraspecific genetic structure indicates that the recent diversity of both species has been predominantly shaped by Pleistocene climatic oscillations, with glacial refugia mainly located in the Balkans, Crimea, and/or Anatolia in E. sauromates and Anatolia and/or the Caucasus in E. urartica sp. nov.

Multilocus phylogeny and coalescent species delimitation in Kotschy's gecko, Mediodactylus kotschyi: Hidden diversity and cryptic species.

Kotschy's Gecko, Mediodactylus kotschyi, is a small gecko native to southeastern Europe and the Levant. It displays great morphological variation with a large number of morphologically recognized subspecies. However, it has been suggested that it constitutes a species complex of several yet unrecognized species. In this study, we used multilocus sequence data (three mitochondrial and three nuclear gene fragments) to estimate the phylogenetic relationships of 174 specimens from 129 sampling localities, covering a substantial part of the distribution range of the species. Our results revealed high genetic diversity of M. kotschyi populations and contributed to our knowledge about the phylogenetic relationships and the estimation of the divergence times between them. Diversification within M. kotschyi began approximately 15 million years ago (Mya) in the Middle Miocene, whereas the diversification within most of the major clades have been occurred in the last 5 Mya. Species delimitation analysis suggests there exists five species within the complex, and we propose to tentatively recognize the following taxa as full species: M. kotschyi (mainland Balkans, most of Aegean islands, and Italy), M. orientalis (Levant, Cyprus, southern Anatolia, and south-eastern Aegean islands), M. danilewskii (Black Sea region and south-western Anatolia), M. bartoni (Crete), and M. oertzeni (southern Dodecanese Islands). This newly recognized diversity underlines the complex biogeographical history of the Eastern Mediterranean region.

A genomic footprint of hybrid zone movement in crested newts.

Speciation typically involves a stage in which species can still exchange genetic material. Interspecific gene flow is facilitated by the hybrid zones that such species establish upon secondary contact. If one member of a hybridizing species pair displaces the other, their hybrid zone would move across the landscape. Although theory predicts that moving hybrid zones quickly stagnate, hybrid zones tracked over one or a few decades do not always follow such a limitation. This suggests that hybrid zones have the potential to traverse considerable distances over extended periods of time. When hybrid zones move, introgression is predicted to result in biased gene flow of selectively neutral alleles, from the receding species into the advancing species. We test for such a genomic footprint of hybrid zone movement in a pair of crested newt species (genus Triturus) for which we have a priori support for westward hybrid zone movement. We perform a multilocus phylogeographical survey and conduct Bayesian clustering analysis, estimation of ancestry and heterozygosity, and geographical cline analysis. In a 600 km wide area east of the present day hybrid zone a genomic footprint constitutes empirical evidence consistent with westward hybrid zone movement. The crested newt case suggests that hybrid zone movement can occur over an extensive span of time and space. Inferring hybrid zone movement provides fundamental insight into historical biogeography and the speciation process, and we anticipate that hybrid zones will prove to be far more mobile than currently appreciated.

A reinvestigation of phylogeny and divergence times of the Ablepharus kitaibelii species complex (Sauria, Scincidae) based on mtDNA and nuDNA genes.

Morphological and DNA data support that the East Mediterranean snake-eyed skink Ablepharus kitaibelii represents a species complex that includes four species A. kitaibelii, A. budaki, A. chernovi, and A. rueppellii, highlighting the need of its taxonomic reevaluation. Here, we used Bayesian and Maximum Likelihood methods to estimate the phylogenetic relationships of all members of the complex based on two mitochondrial (cyt b, 16S rRNA) and two nuclear markers (MC1R, and NKTR) and using Chalcides, Eumeces, and Eutropis as outgroups. The biogeographic history of the complex was also investigated through the application of several phylogeographic (BEAST) and biogeographic (BBM) analyses. Paleogeographic and paleoclimatic data were used to support the inferred phylogeographic patterns. The A. kitaibelli species complex exhibits high genetic diversity, revealing cases of hidden diversity and cases of non-monophyletic species such as A. kitaibelii and A. budaki. Our results indicate that A. pannonicus branches off first and a group that comprises specimens of A. kitaibelli and A. budaki from Kastelorizo Island group (southeast Greece) and southwest Turkey, respectively is differentiated from the rest A. kitaibelli and A. budaki populations and may represent a new species. The estimated divergence times place the origin of the complex in the Middle Miocene (∼16Mya) and the divergence of most currently recognized species in the Late Miocene. The inferred ancestral distribution suggests that the complex originated in Anatolia, supposing that several vicariance and dispersal events that are related with the formation of the Mid-Aegean Trench, the Anatolian Diagonal and the orogenesis of the mountain chains in southern and eastern Anatolia have led to current distribution pattern of A. kitaibelii species complex in the Balkans and Middle East.

Phylogeny and diversification of mountain vipers (Montivipera, Nilson et al., 2001) triggered by multiple Plio-Pleistocene refugia and high-mountain topography in the Near and Middle East.

The Near and Middle East is a hotspot of biodiversity, but the region remains underexplored at the level of genetic biodiversity. Here, we present an extensive molecular phylogeny of the viperid snake genus Montivipera, including all known taxa. Based on nuclear and mitochondrial data, we present novel insights into the phylogeny of the genus and review the status of its constituent species. Maximum likelihood methods revealed a montane origin of Montivipera at 12.3Mya. We then analyzed factors of mountain viper diversity. Our data support substantial changes in effective population size through Plio-Pleistocene periods. We conclude that climatic oscillations were drivers of allopatric speciation, and that mountain systems of the Near and Middle East have strongly influenced the evolution and survival of taxa, because climatic and topographical heterogeneities induced by mountains have played a crucial role as filters for dispersal and as multiple refugia. The wide diversity of montane microhabitats enabled mountain vipers to retain their ecological niche during climatic pessima. In consequence the varied geological and topographical conditions between refugia favoured genetic isolation and created patterns of species richness resulting in the formation of neoendemic taxa. Our data support high concordance between geographic distributions of Montivipera haplotypes with putative plant refugia.

External morphology and osteology of Darevskia rudis (Bedriaga, 1886), with a taxonomic revision of the Pontic and Small-Caucasus populations (Squamata: Lacertidae).

A broad sample of Darevskia rudis from the main part of its range was reviewed with regard to external morphology (discriminant, UPGMA, MST and ANOVA analyses) and osteology. Darevskia bithynica is raised to species rank, with two subspecies: D. b. bithynica and D. b. tristis. The other subspecies are fairly similar (D. r. rudis being the most different). Two singular populations are described as subspecies: D. r. mirabilis ssp. nov. from Kaçkar Mountains, geographically adjoins the otherwise different D. r. bischoffi and D. r. bolkardaghica ssp. nov., which is geographically isolated but that seems to be very closely related to D. r. obscura.

Molecular phylogeny and historical biogeography of the Anatolian lizard Apathya (Squamata, Lacertidae).

Apathya is a lacertid genus occurring mainly in south-east Turkey and its adjacent regions (part of Iran and Iraq). So far two morphological species have been attributed to the genus; A. cappadocica (with five subspecies, A. c.cappadocica, A. c.muhtari, A. c.schmidtlerorum, A. c. urmiana and A. c.wolteri) and A.yassujica. The first species occupies most of the genus' distribution range, while A. yassujica is endemic of the Zagros Mountains. Here, we explored Apathya's taxonomy and investigated the evolutionary history of the species by employing phylogenetic and phylogeographic approaches and using both mitochondrial (mtDNA) and nuclear markers. The phylogenetic relationships and the genetic distances retrieved, revealed that Apathya is a highly variable genus, which parallels its high morphological variation. Such levels of morphological and genetic differentiation often exceed those between species of other Lacertini genera that are already treated as full species, suggesting the necessity for a taxonomic revision of Apathya. The phylogeographical scenario emerging from the genetic data suggests that the present distribution of the genus was determined by a combination of dispersal and vicariance events between Anatolia and Southwest Asia dating back to the Miocene and continuing up to the Pleistocene. Key geological events for the understanding of the phylogeography of the genus are the movement of the Arabian plate that led to the configuration of Middle East (orogenesis of the mountain ranges of Turkey and Iran) and the formation of Anatolian Diagonal.