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.

The Silk roads: phylogeography of Central Asian dice snakes (Serpentes: Natricidae) shaped by rivers in deserts and mountain valleys

Influenced by rapid changes in climate and landscape features since the Miocene, widely distributed species provide suitable models to study the environmental impact on their evolution and current genetic diversity. The dice snake Natrix tessellata, widely distributed in the Western Palearctic is one such species. We aimed to resolve a detailed phylogeography of N. tessellata with a focus on the Central Asian clade with 4 and the Anatolia clade with 3 mitochondrial lineages, trace their origin, and correlate the environmental changes that affected their distribution through time. The expected time of divergence of both clades began at 3.7 Mya in the Pliocene, reaching lineage differentiation approximately 1 million years later. The genetic diversity in both clades is rich, suggesting different ancestral areas, glacial refugia, demographic changes, and colonization routes. The Caspian lineage is the most widespread lineage in Central Asia, distributed around the Caspian Sea and reaching the foothills of the Hindu Kush Mountains in Afghanistan, and Eastern European lowlands in the west. Its distribution is limited by deserts, mountains, and cold steppe environments. Similarly, Kazakhstan and Uzbekistan lineages followed the Amu Darya and the Syr Darya water systems in Central Asia, with ranges delimited by the large Kyzylkum and Karakum deserts. On the western side, there are several lineages within the Anatolia clade that converged in the central part of the peninsula with 2 being endemic to Western Asia. The distribution of both main clades was affected by expansion from their Pleistocene glacial refugia around the Caspian Sea and in the valleys of Central Asia as well as by environmental changes, mostly through aridification.

The first comprehensive data on the distribution of reptiles within the Southern Bug eco-corridor, Ukraine

This study provides distribution data for 12 species of reptiles in the Southern Bug eco-corridor located within the steppe zone, Mykolaiv Oblast (province), Ukraine. We compiled 935 records from published literature (324), from public databases (68), and from 12 field surveys we made in 2016–2020 (543). All records were mapped on to a grid of 10×10 km, representing the 294 cells of the studied area. We present new distributional data for Podarcis tauricus. For all the other species, our records add new localities that fall within areas where these reptiles were previously reported. Species richness and Shannon’s H’ index show that herpetofauna diversity was higher in the areas along the Southern Bug River. The maximum number of records within a single cell was 118 (for 10 species) in this oblast. Also the high number of identified reptile chorotypes (nine) within the Mykolaiv Oblast revealed that southern Ukraine is an important zoogeographical territory in Europe. These data provide the basis for future biogeographical and ecological studies and conservation priorities.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in > 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Landscape Genomics of a Widely Distributed Snake, Dolichophis caspius (Gmelin, 1789) across Eastern Europe and Western Asia

Across the distribution of the Caspian whipsnake (Dolichophis caspius), populations have become increasingly disconnected due to habitat alteration. To understand population dynamics and this widespread but locally endangered snake's adaptive potential, we investigated population structure, admixture, and effective migration patterns. We took a landscape-genomic approach to identify selected genotypes associated with environmental variables relevant to D. caspius. With double-digest restriction-site associated DNA (ddRAD) sequencing of 53 samples resulting in 17,518 single nucleotide polymorphisms (SNPs), we identified 8 clusters within D. caspius reflecting complex evolutionary patterns of the species. Estimated Effective Migration Surfaces (EEMS) revealed higher-than-average gene flow in most of the Balkan Peninsula and lower-than-average gene flow along the middle section of the Danube River. Landscape genomic analysis identified 751 selected genotypes correlated with 7 climatic variables. Isothermality correlated with the highest number of selected genotypes (478) located in 41 genes, followed by annual range (127) and annual mean temperature (87). We conclude that environmental variables, especially the day-to-night temperature oscillation in comparison to the summer-to-winter oscillation, may have an important role in the distribution and adaptation of D. caspius.