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

Multiple Ecological Niche Modeling Reveals Niche Conservatism and Divergence in East Asian Yew (Taxus)

Understanding ecological niche evolution patterns is crucial for elucidating biogeographic history and guiding biodiversity conservation. Taxus is a Tertiary relict gymnosperm with 11 lineages mainly distributed across East Asia, spanning from tropical to subarctic regions. However, the spatiotemporal dynamics of its ecological niche evolution and the roles of ecological and geographical factors in lineage diversification, remain unclear. Using occurrence records, environmental data, and reconstructed phylogenies, we employed ensemble ecological niche models (eENMs), environmental principle components analysis (PCA-env), and phyloclimatic modeling to analyze niche similarity and evolution among 11 Taxus lineages. Based on reconstructed Bayesian trees and geographical distribution characteristics, we classified the eleven lineages into four clades: Northern (T. cuspidata), Central (T. chinensis, T. qinlingensis, and the Emei type), Western (T. wallichiana, T. florinii, and T. contorta), and Southern (T. calcicola, T. phytonii, T. mairei, and the Huangshan type). Orogenic activities and climate changes in the Tibetan Plateau since the Late Miocene likely facilitated the local adaptation of ancestral populations in Central China, the Hengduan Mountains, and the Yunnan-Guizhou Plateau, driving their expansion and diversification towards the west and south. Key environmental variables, including extreme temperature, temperature and precipitation variability, light, and altitude, were identified as major drivers of current niche divergence. Both niche conservatism and divergence were observed, with early conservatism followed by recent divergence. The Southern clade exhibits high heat and moisture tolerance, suggesting an adaptive shift, while the Central and Western clades retain ancestral drought and cold tolerance, displaying significant phylogenetic niche conservatism (PNC). We recommend prioritizing the conservation of T. qinlingensis, which exhibits the highest PNC level, particularly in the Qinling, Daba, and Taihang Mountains, which are highly degraded and vulnerable to future climate fluctuations.

Whole genome data confirm pervasive gene discordance in the evolutionary history of Coenonympha (Nymphalidae) butterflies

Phylogenetic inference is challenged by genealogical heterogeneity amongst molecular markers. Such discordance is driven predominantly by incomplete lineage sorting (ILS) and interspecific gene flow, and bears attendant consequences for the accurate resolution of species relationships. Understanding the distribution of gene conflict in organismal genomes is, therefore, a key aspect of phylogenetic analysis. In this study, three large phylogenomic datasets (i.e., whole mitogenomes, conserved nuclear protein-coding loci, and genomic windows) are used to probe the extent to which discordance pervades the unresolved phylogeny of Coenonympha (Nymphalidae) butterflies. Gene tree discordance is found to be elevated at multiple historically recalcitrant phylogenetic positions. In particular, species relationships near the crown of Coenonympha and within a rapidly diversifying subclade (the hero group) remain difficult to resolve, suggesting that ILS and gene flow have obscured the evolution of this genus. These findings have implications for the taxonomy of this butterfly group and the study of its diversification history. In addition, this work lends support to a growing body of evidence that gene conflict driven by biological processes stands to confound phylogeny, even when extensive data are used.

Biogeography of the Iranian snakes

The events of the Cenozoic era such as mountain formation caused Iran to become one of the most amazing biodiversity hotspots in the world today. This pioneering study on Iranian snake biogeography integrates historical and ecological analyses. A phylogeographic review traces speciation and dispersal, while cluster analysis with a new snake checklist assesses faunistic similarities within Iran and its surroundings. Jaccard and Sorenson indices generate similarity dendrograms, Indicator Species Analysis pinpoints regional key species, and Endemism index calculates regional endemism rates, enriching our knowledge of Iran's species diversity. Phylogeographic analyses identify four biogeographical corridors for snake ingress into Iran: the Arabian region through southwestern Iran, the Western Asian mountainous transition zone via northwestern Iran, the Turanian region into northeastern Iran, and the Indus River Valley into southeastern and eastern Iran. Dendrogram analysis divides snake fauna into three groups. The first group associates western Zagros and Khuzestan fauna with the Sahara and Arabian regions. The second group links Kopet Dagh and Turkmen Steppe fauna with the Turanian region, and Central Plateau and Baluchistan fauna with the Iranian region. The third group connects northwest highlands, Alborz and Zagros mountains, and Caspian Sea coasts with the Western Asian Mountain transition zone. The study validates broad biogeographic patterns via ecoregional associations and indicator species analysis, providing finer resolution. Species like Platyceps najadum in Caspian Hyrcanian mixed forests exemplify ecoregional alignment, while Zagros and Alborz mountains exhibit unique faunal indicators, indicating species-level divergence. Shared indicators among widespread ecoregions reflect habitat continuity; exclusive indicators emphasize regional distinctiveness. Despite endemic species prevalence, they seldom act as significant indicators due to various factors. Our research confirms the Zagros Mountains, Khuzestan Plain, Alborz Mountains, and Persian Gulf coasts as snake diversity hotspots, marked by higher species richness compared to other Iranian regions.

Turtles and tortoises of Syria: Diversity, distribution, and conservation

Syria is home to a wide range of ecosystems and landscapes that support a variety of animal species, including reptiles. Currently, Syria hosts nine species of chelonians (turtles and tortoises) from six different families, making its family-level diversity remarkable for the Mediterranean and the Middle East. We conducted a comprehensive review of existing literature, database information, and museum records, and incorporated new field data collected between 2019 and 2023 to document the distribution and conservation status of chelonians in Syria. As a result, we present an updated species list, accompanied by relevant information and occurrence maps for each turtle and tortoise species found in the country. Our investigation has yielded a total of 321 distribution records for Syrian chelonians (including nests of marine species), 105 of which are new. In addition, we verified 18 sites previously documented in the literature. Among the most frequently recorded species is Testudo graeca terrestris Forskål, 1775. In addition, for the first time, we have systematically evaluated the records of two marine turtle species that use Syrian coasts as nesting sites. This study is the most comprehensive and up-to-date source of information on the distribution and biogeography of chelonians in Syria, encompassing both historical and current geographic records. In addition, our data provide a zoogeographic overview and a substantial dataset of distribution information that can be used for subsequent analysis of macroecological patterns, distribution patterns, and conservation efforts related to this group of reptiles.

Phylogeography and Ecological Niche Modeling of the Alashan Pit Viper (Gloydius cognatus; Reptilia, Viperidae) in Northwest China and Adjacent Areas

The joint impacts of historical geological events and Quaternary climatic oscillations in Northwest China on species evolution have been examined extensively in plant under a phylogeographic perspective. However, animal phylogeographic analyses in this region are still limited. The Alashan pit viper, Gloydius cognatus, occurs primarily in arid Northwest China and adjacent areas. Based on variation at two mtDNA genes (ND4 and Cytb) in 27 individuals representing 24 populations, the spatial genetic structure and demographic history of G. cognatus were examined across its geographic range. Phylogenetic analyses revealed two well-supported allopatric clades (each with two distinct subclades/lineages), distributed across the southern (Qaidam Basin, Lanzhou Basin, and Zoige Basin [S1]; Loess Plateau [S2]) and northern (Ily Basin [N1]; Junggar Basin and Mongolian Plateau [N2]) regions. AMOVA analysis demonstrated that over 76% of the observed genetic variation was related to these lineage splits, indicating substantial genetic differentiation among the four lineages. A strong pattern of isolation-by-distance across the sampling populations suggested that geographic distance principally shaped the genetic structure. The four lineages diverged by 0.9-2.2% for the concatenated data, which were estimated to have coalesced ~1.17 million years ago (Mya), suggesting that the expansions of the Badain Jaran, Tengger, and Mu Us deserts during the Xixiabangma glaciation likely interrupted gene flow and triggered the observed divergence in the southern and northern regions. Subsequently, the early Pleistocene integration of the Yellow River and associated deserts expansion promoted the differentiation of S1 and S2 lineages (~0.9 Mya). Both mitochondrial evidence and ecological niche modeling (ENM) reject the signature of demographic and range contractions during the LGM for G. cognatus. In addition, ENM predicts that the suitable habitat of G. cognatus will contract in the future. As such, the conservation and management of ESUs should be a priority. Our findings provide the first insights on the lineage diversification and population dynamics of the Alashan pit viper in relation to geological history and Pleistocene climatic oscillations in arid Northwest China.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

A new, rare, small-ranged, and endangered mountain snake of the genus Elaphe from the Southern Levant

The genus Elaphe Fitzinger, 1833 includes 17 species of charismatic, large-sized, non-venomous, Eurasian snakes. In the Western Palearctic, the genus is represented by three species from the Elaphe quatuorlineata group ranging from the Apennine peninsula to Central Asia. The southernmost population of this group is distributed in the mountains of the Southern Levant, with more than 400 km gap to other Elaphe populations. This population has been known to science for only 50 years and is virtually unstudied due to its extreme rarity. We studied these snakes' morphological and genetic variation from the three countries where they are known to occur, i.e., Israel (Hermon, the Israeli-controlled Golan Heights), Lebanon, and Syria. We used nine mitochondrial and nuclear genes, complete mitogenome sequences, and a comprehensive morphological examination including published data, our own field observations, and museum specimens, to study its relationship to other species in the group. The three currently recognized species of the group (E. quatuorlineata, E. sauromates, E. urartica), and the Levant population, form four deeply divergent, strongly supported clades. Three of these clades correspond to the abovementioned species while the Southern Levant clade, which is genetically and morphologically distinct from all named congeners, is described here as a new species, Elaphe druzei sp. nov. The basal divergence of this group is estimated to be the Late Miocene with subsequent radiation from 5.1 to 3.9 Mya. The revealed biogeography of the E. quatuorlineata group supports the importance of the Levant as a major center of endemism and diversity of biota in Eurasia. The new species is large-sized and is one of the rarest snakes in the Western Palearctic. Because of its small mountain distribution range, in an area affected by land use and climate change, the new Elaphe urgently needs strict protection. Despite political issues, we hope this will be based on the cooperation of all countries where the new species occurs.

Local envenoming by the coinsnake or Asian racer, Hemorrhois nummifer and mountain racer or leopard snake, Hemorrhois ravergieri (Serpentes: Colubridae, Colubrinae) in Iran: A reminder of the importance of species identification in the medical management of snakebites

In Iran, there are approximately 4500-6500 snakebites per year, but fortunately only 3-9 of these are fatal. However, in some population centers such as Kashan city (Isfahan Province, central Iran), approximately 80% of snakebites are attributed to "non-venomous" snakes that are often comprised of several species of non-front-fanged snakes (NFFS). NFFS comprise a diverse group that constitute approximately 2900 species belonging to an estimated 15 families. We report here two cases of local envenoming from H. ravergieri, and one from H. nummifer that occurred in Iran. The clinical effects consisted of local erythema, mild pain, transient bleeding and edema. Two victims experienced progressive local edema that distressed the victims. The medical team's unfamiliarity with snakebites contributed to the incorrect clinical management of one victim including the contraindicated, ineffective provision of antivenom. These cases provide further documentation about local envenoming caused by these species, and also emphasize the need for regional medical personnel to receive increased training in order to improve familiarity with the local snake fauna and evidence-based snakebite management.

MaxEnt brings comparable results when the input data are being completed; Model parameterization of four species distribution models

Species distribution models (SDMs) are practical tools to assess the habitat suitability of species with numerous applications in environmental management and conservation planning. The manipulation of the input data to deal with their spatial bias is one of the advantageous methods to enhance the performance of SDMs. However, the development of a model parameterization approach covering different SDMs to achieve well-performing models has rarely been implemented. We integrated input data manipulation and model tuning for four commonly-used SDMs: generalized linear model (GLM), gradient boosted model (GBM), random forest (RF), and maximum entropy (MaxEnt), and compared their predictive performance to model geographically imbalanced-biased data of a rare species complex of mountain vipers. Models were tuned up based on a range of model-specific parameters considering two background selection methods: random and background weighting schemes. The performance of the fine-tuned models was assessed based on recently identified localities of the species. The results indicated that although the fine-tuned version of all models shows great performance in predicting training data (AUC > 0.9 and TSS > 0.5), they produce different results in classifying out-of-bag data. The GBM and RF with higher sensitivity of training data showed more different performances. The GLM, despite having high predictive performance for test data, showed lower specificity. It was only the MaxEnt model that showed high predictive performance and comparable results for identifying test data in both random and background weighting procedures. Our results highlight that while GBM and RF are prone to overfitting training data and GLM over-predict nonsampled areas MaxEnt is capable of producing results that are both predictable (extrapolative) and complex (interpolative). We discuss the assumptions of each model and conclude that MaxEnt could be considered as a practical method to cope with imbalanced-biased data in species distribution modeling approaches.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

First report of a confirmed case of Montivipera latifii (Latifi's viper) envenoming and a literature review of envenoming by Montivipera species

Latifi's viper (Montivipera latifii), also known as Lar Valley or Damavandi viper, is endemic to Iran. It has rarely been recorded, as it occurs in a highly-protected national park. In this first clinical report of a confirmed bite by this species, a teenage girl was bitten on the chin, causing rapidly-progressive swelling of the face and oropharyngeal mucosa. At a local hospital, a misleading history given by the patient's relatives of a wasp sting and inadequate inspection of the bite wound misled the physicians from making the correct diagnosis, resulting in a considerable delay in the administration of antivenom. This allowed the development of partial obstruction of the upper airway causing respiratory distress. After transfer to a tertiary hospital, attempts at endotracheal intubation failed, necessitating tracheostomy, but this was not implemented early enough to prevent her developing respiratory failure and losing consciousness. After she was stabilized, snakebite envenoming was diagnosed by a clinical toxicologist who observed two fang puncture marks on her chin. This was later confirmed when a snake, identified as M. latifii, was discovered at the room where the bite had occurred. Her facial swelling and ecchymosis, attributable to envenoming, were effectively controlled by high-dose antivenom therapy. However, she did not recover consciousness, remaining in a vegetative state. About three weeks after the bite, she died as an indirect result of hypoxic brain damage complicated by septicemia. Prompt diagnosis, relief of upper airway obstruction and timely antivenom therapy might have prevented this tragic fatal outcome.

The legacy of Eastern Mediterranean mountain uplifts: rapid disparity of phylogenetic niche conservatism and divergence in mountain vipers

BACKGROUND

The orogeny of the eastern Mediterranean region has substantially affected ecological speciation patterns, particularly of mountain-dwelling species. Mountain vipers of the genus Montivipera are among the paramount examples of Mediterranean neo-endemism, with restricted ranges in the mountains of Anatolia, the Levant, Caucasus, Alborz, and Zagros. Here we explore the phylogenetic and ecological diversification of Montivipera to reconstruct its ecological niche evolution and biogeographic history. Using 177 sequences of three mitochondrial genes, a dated molecular phylogeny of mountain vipers was reconstructed. Based on 320 occurrence points within the entire range of the genus and six climatic variables, ecological niches were modelled and used to infer ancestral niche occupancy. In addition, the biogeographic history and ancestral states of the species were reconstructed across climate gradients.

RESULTS

Dated phylogenetic reconstruction revealed that the ancestor of mountain vipers split into two major clades at around 12.18 Mya followed by multiple vicariance events due to rapid orogeny. Montivipera colonised coastal regions from a mountain-dwelling ancestor. We detected a highly complex ecological niche evolution of mountain vipers to temperature seasonality, a variable that also showed a strong phylogenetic signal and high contribution in niche occupation.

CONCLUSION

Raising mountain belts in the Eastern Mediterranean region and subsequent remarkable changes in temperature seasonality have led to the formation of important centres of diversification and endemism in this biodiversity hotspot. High rates of niche conservatism, low genetic diversity, and segregation of ranges into the endemic distribution negatively influenced the adaptive capacity of mountain vipers. We suggest that these species should be considered as evolutionary significant units and priority species for conservation in Mediterranean mountain ecosystems.

Venom-Induced Blood Disturbances by Palearctic Viperid Snakes, and Their Relative Neutralization by Antivenoms and Enzyme-Inhibitors

Palearctic vipers are medically significant snakes in the genera Daboia, Macrovipera, Montivipera, and Vipera which occur throughout Europe, Central Asia, Near and Middle East. While the ancestral condition is that of a small-bodied, lowland species, extensive diversification has occurred in body size, and niche specialization. Using 27 venom samples and a panel of in vitro coagulation assays, we evaluated the relative coagulotoxic potency of Palearctic viper venoms and compared their neutralization by three antivenoms (Insoserp Europe, VIPERFAV and ViperaTAb) and two metalloprotease inhibitors (prinomastat and DMPS). We show that variation in morphology parallels variation in the Factor X activating procoagulant toxicity, with the three convergent evolutions of larger body sizes (Daboia genus, Macrovipera genus, and Vipera ammodytes uniquely within the Vipera genus) were each accompanied by a significant increase in procoagulant potency. In contrast, the two convergent evolutions of high altitude specialization (the Montivipera genus and Vipera latastei uniquely within the Vipera genus) were each accompanied by a shift away from procoagulant action, with the Montivipera species being particularly potently anticoagulant. Inoserp Europe and VIPERFAV antivenoms were both effective against a broad range of Vipera species, with Inoserp able to neutralize additional species relative to VIPERFAV, reflective of its more complex antivenom immunization mixture. In contrast, ViperaTAb was extremely potent in neutralizing V. berus but, reflective of this being a monovalent antivenom, it was not effective against other Vipera species. The enzyme inhibitor prinomastat efficiently neutralized the metalloprotease-driven Factor X activation of the procoagulant venoms. In contrast, DMPS (2,3-dimercapto-1-propanesulfonic acid), which as been suggested as another potential treatment option in the absence of antivenom, DMPS failed against all venoms tested. Overall, our results highlight the evolutionary variations within Palearctic vipers and help to inform clinical management of viper envenomation.

Macroevolutionary variation and environmental correlates of scalation traits in Eurasian vipers (Serpentes: Viperinae)

Understanding how phenotypic variation across species is shaped by the combination of shared evolutionary history and environmental factors is key to elucidating the processes that underlie biodiversity. In reptiles, morphological traits have traditionally been used to delimit species and make systematic inferences. Recent studies highlight the possibility that phenotypic variation, particularly in scalation traits, might instead be driven by environmental factors and therefore not reflect the phylogenetic relationships among species. In this study, we combined morphological and ecological data in a macroevolutionary framework, in order to describe the morphological variation across species of Eurasian vipers (Serpentes: Viperinae), investigate the phylogenetic structure of scalation traits and test the contribution of environmental factors in shaping morphological patterns. We found considerable variation in all examined traits, which, in most cases, agreed with the phylogenetic relationships among species, reinforcing their usefulness for taxonomic inferences. Interestingly, however, the number of ventral scales exhibited lower phylogenetic signal and a tight association with environmental factors of geographical ranges, suggesting potential adaptive or developmental sources of variation in the trait. This is the first comparative study of macroevolutionary variation in scalation traits in Eurasian vipers, validating the use of most of them for systematic inferences, but also indicating possible environmental factors that might shape phenotypic variation across species.

Reptile species richness associated to ecological and historical variables in Iran

Spatial gradients of species richness can be shaped by the interplay between historical and ecological factors. They might interact in particularly complex ways in heterogeneous mountainous landscapes with strong climatic and geological contrasts. We mapped the distribution of 171 lizard species to investigate species richness patterns for all species (171), diurnal species (101), and nocturnal species (70) separately. We related species richness with the historical (past climate change, mountain uplifting) and ecological variables (climate, topography and vegetation). We found that assemblages in the Western Zagros Mountains, north eastern and north western parts of Central Iranian Plateau have the highest number of lizard species. Among the investigated variables, annual mean temperature explained the largest variance for all species (10%) and nocturnal species (31%). For diurnal species, temperature change velocity shows strongest explained variance in observed richness pattern (26%). Together, our results reveal that areas with annual temperature of 15–20 °C, which receive 400–600 mm precipitation and experienced moderate level of climate change since the Last Glacial Maximum (LGM) have highest number of species. Documented patterns of our study provide a baseline for understanding the potential effect of ongoing climate change on lizard diversity in Iran.

Thermal melanism explains macroevolutionary variation of dorsal pigmentation in Eurasian vipers

Colouration may endorse thermoregulatory and antipredatory functions in snakes. The thermal melanism hypothesis predicts that dark-coloured individuals are ecologically favoured in cool climates. However, the loss of aposematic and cryptic colourations may imply high predation for melanistic snakes. Here, we used the monophyletic group of Eurasian vipers (subfamily Viperinae) to test whether an increase in the extent of dark area inside the characteristic zigzag dorsal pattern is associated to colder environments. We measured two colouration traits in zigzag-patterned individuals (number of dorsal marks and weighted pigmentation index) and used a phylogenetic comparative approach to explore macroevolutionary patterns of dorsal pigmentation and test whether its extent is associated to ecogeographic characteristics of lineages’ ranges. Phylogenetically-naïve and phylogenetically-informed analyses yielded a significant association between the degree of pigmentation of the zigzag pattern and environmental variables such as solar radiation, elevation and latitude. The degree of pigmentation of the zigzag pattern is highlighted as an adaptive trait that matches range attributes mirroring cold environments irrespective of the phylogeny. These results constitute the first large-scale evidence supporting the thermal melanism hypothesis in snakes, opening new avenues of inquiry for the mechanisms that shape the evolution of colour phenotypes.

Origin, extinction and ancient DNA of a new fossil insular viper: molecular clues of overseas immigration

Viperinae is a subfamily of viperid snakes whose fossil record in the Mediterranean islands is, until now, restricted to 12 palaeontological deposits on seven islands. Revision of the material excavated 30 years ago from the Middle/Late Pleistocene–Holocene deposit of Es Pouàs [Eivissa (= Ibiza), Balearic Islands, western Mediterranean] revealed about 6000 bones of a small-sized viper across different stratigraphic levels. Its morphological characteristics are different enough to known species of Vipera to warrant the description of a new species, but the nearly complete mitochondrial genome obtained from this snake based on a sample dated to 16 130 ± 45 bp, suggested it belonged to a new insular population of Lataste’s viper (Vipera latastei), Vipera latastei ebusitana subsp. nov. Phylogenetic analysis indicates that the dispersal of the ancestors of V. l. ebusitana to Eivissa, most probably from a north-east Iberian population, occurred via overwater colonization < 1.5 Mya, well after the Messinian Salinity Crisis (5.97–5.32 Mya) when land bridges allowed terrestrial colonization of the Balearic Islands by mainland faunas. The morphological differences between V. l. ebusitana and the Iberian populations suggest that it is a new dwarf taxon resulting from insular evolutionary processes, becoming extinct shortly after the first human arrival to this island about 4000 years ago.

The Anatolian glacial refugium and human-mediated colonization: a phylogeographical study of the stone marten (Martes foina) in Turkey

The Anatolian Peninsula, comprising most of modern Turkey, is situated at the junction of Europe, the Middle East and Asia. Together with its complex geomorphological and climatic history, this has given rise to a rich fauna and flora, which exhibits a wide range of historical biogeographical patterns. The stone marten (Martes foina) is a small carnivore that is found across the temperate Palaearctic region including Anatolia, where it is often associated with habitats modified by humans, but few genetic data exist for this species. We sequenced a 1840-bp region of the mitochondrial genome from 97 martens sampled across the peninsula and intron 7 of the nuclear β-fibrinogen gene from 53 of these. Two mitochondrial lineages were recovered, with overlapping eastern and western distributions, but there was no geographical structure for the autosomal marker. Coalescent analyses indicated that both of the lineages originated during the Last Glacial Maximum, one of them within an eastern Anatolian refugium and the other in a western Anatolian or Balkan refugium. The western lineage colonized most of Europe in the Holocene, while the eastern lineage may be endemic to Anatolia, from where it colonized the Iberian Peninsula via human translocation. The presence of at least one refugial stone marten population highlights the importance of Anatolia to the preservation of genetic variation and biodiversity.

Conservation Below the Species Level: Suitable Evolutionarily Significant Units among Mountain Vipers (the Montivipera raddei complex) in Iran

Northern and western mountains of Iran are among the most important biodiversity and endemism hot spots for reptiles in the Middle East. Among herpetofauna, the montivipers represent an emblematic and fragmented endemic group for which estimating their level of genetic differentiation and defining conservation priorities is urgently needed. Here, we present the most comprehensive phylogenetic study on the Montivipera raddei species group comprising all 5 known taxa, among which 3 are endemic to Iran. Based on 2 mitochondrial genes, phylogenetic and phylogeographic analyses revealed 3 major lineages each presenting very contrasting distribution areas. The Iranian montivipers are highly structured in clades showing low genetic diversity and corresponding to high altitude summits. Molecular dating revealed the role of Quaternary paleo-climatic oscillations and altitudinal movements of montivipers in shaping genetic diversity and differentiation of these sky-island taxa. In addition, the best scenario of historical biogeography allowed identifying 3 possible refugial areas in Iran most likely arising by vicariance. Based on our mitochondrial results and pending additional data, we recognize 3 candidate species among the M. raddei complex: M. raddei, Montivipera latifii, and Montivipera kuhrangica that are coherent with their geographical distribution. We propose that the most appropriate evolutionary significant units for conservation of the montivipers are represented by 13 units among which 6 are recognized as high priority. Finally, we suggest some recommendations to the IUCN as well as to the Iranian conservation policies with respect to conservation prioritization.

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.

Extinction risks of a Mediterranean neo-endemism complex of mountain vipers triggered by climate change

Climate change is among the most important drivers of biodiversity decline through shift or shrinkage in suitable habitat of species. Mountain vipers of the genus Montivipera are under extreme risk from climate changes given their evolutionary history and geographic distribution. In this study, we divided all Montivipera species into three phylogenetic-geographic Montivipera clades (PGMC; Bornmuelleri, Raddei and Xanthina) and applied an ensemble ecological niche modelling (ENM) approach under different climatic scenarios to assess changes in projected suitable habitats of these species. Based on the predicted range losses, we assessed the projected extinction risk of the species relative to IUCN Red List Criteria. Our result revealed a strong decline in suitable habitats for all PGMCs (63.8%, 79.3% and 96.8% for Xanthina, Raddei and Bornmuelleri, respectively, by 2070 and under 8.5 RCP scenario) with patterns of altitudinal range shifts in response to projected climate change. We found that the mountains close to the Mediterranean Sea are exposed to the highest threats in the future (84.6 ± 9.1 percent range loss). We also revealed that disjunct populations of Montivipera will be additionally highly isolated and fragmented in the future. We argue that leveraging climate niche projections into the risk assessment provides the opportunity to implement IUCN criteria and better assess forthcoming extinction risks of species.

Cryptic diversity in the smooth snake (Coronella austriaca)

The smooth snake, Coronella austriaca, is a common snake species widespread in the Western Palearctic region. It does not form conspicuous morphological variants and, although several evolutionary lineages have been distinguished based on the analyses of the mitochondrial DNA sequences, only two subspecies with very limited distribution have been traditionally recognized. Here we present an mtDNA phylogeography of the species using geographically extended sampling while incorporating biogeographically important areas that have not been analyzed before, such as Anatolia, Crimea, and Iran. We find that the smooth snake comprises 14 distinct phylogenetic clades with unclear mutual relationships, characterized by complex genetic structure and relatively deep divergences; some of them presumably of Miocene origin. In general, the biogeographic pattern is similar to other Western Palearctic reptiles and illustrates the importance of the main European peninsulas as well as the Anatolian mountains, Caucasus, and Alborz Mts. in Iran for the evolution of the present-day diversity. Considerable genetic structure present in the smooth snake populations within these large areas indicates the existence of several regional Plio-Pleistocene refugia that served as reservoirs for dispersal and population expansions after the glacial periods. The current taxonomy of C. austriaca does not reflect the rich genetic diversity, deep divergences, and overall evolutionary history revealed in our study and requires a thorough revision. This will only be possible with an even higher-resolution sampling and integrative approach, combining analyses of multiple genetic loci with morphology, and possibly other aspects of the smooth snake biology.

Rapid and recent diversification patterns in Anseriformes birds: Inferred from molecular phylogeny and diversification analyses

The Anseriformes is a well-known and widely distributed bird order, with more than 150 species in the world. This paper aims to revise the classification, determine the phylogenetic relationships and diversification patterns in Anseriformes by exploring the Cyt b, ND2, COI genes and the complete mitochondrial genomes (mito-genomes). Molecular phylogeny and genetic distance analyses suggest that the Dendrocygna species should be considered as an independent family, Dendrocygnidae, rather than a member of Anatidae. Molecular timescale analyses suggests that the ancestral diversification occurred during the Early Eocene Climatic Optimum (58 ~ 50 Ma). Furthermore, diversification analyses showed that, after a long period of constant diversification, the median initial speciation rate was accelerated three times, and finally increased to approximately 0.3 sp/My. In the present study, both molecular phylogeny and diversification analyses results support that Anseriformes birds underwent rapid and recent diversification in their evolutionary history, especially in modern ducks, which show extreme diversification during the Plio-Pleistocene (~ 5.3 Ma). Therefore, our study support that the Plio-Pleistocene climate fluctuations are likely to have played a significant role in promoting the recent diversification for Anseriformes.

Combined venom profiling and cytotoxicity screening of the Radde's mountain viper (Montivipera raddei) and Mount Bulgar Viper (Montivipera bulgardaghica) with potent cytotoxicity against human A549 lung carcinoma cells

Here we report the first characterization of the endemic Mount Bulgar Viper (Montivipera bulgardaghica) and Radde's mountain viper (Montivipera raddei) venom by a combined approach using intact mass profiling and bottom-up proteomics. The cytotoxicity screening of crude venom as well as isolated serine proteases revealed a high activity against A549 human lung carcinoma cells. By means of intact mass profiling of native and reduced venom we observed basic and acidic phospholipases type A₂. Moreover, the analysis revealed snake venom metalloproteases, cysteine-rich secretory proteins, disintegrins, snake venom serine proteases, C-type lectins, a vascular endothelial growth factor and an L-amino acid oxidase.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

A new population and subspecies of the critically endangered Anatolian meadow viper Vipera anatolica Eiselt and Baran, 1970 in eastern Antalya province

We report on a new population of Vipera anatolica from the Geyik Mountain Range in eastern Antalya Province, Turkey. It represents only the second known location, and is situated in a valley about 200 km east from the terra typica at Kohu Dağ in western Antalya Province. We compare both populations and, based on marked differences in morphology, habitat, genetics, and its geographically isolated location, we describe the recently discovered population as a new subspecies. Aspects of ecology, threats, and conservation needs are discussed.