Integrative taxonomy provides evidence for the species status of the Ibero-Maghrebian grass snakeNatrix astreptophora

The Endangered Sardinian Grass Snake: Distribution Update, Bioclimatic Niche Modelling, Dorsal Pattern Characterisation, and Literature Review

The Sardinian grass snake, Natrix helvetica cetti, is an endangered endemic snake subspecies with a restricted and highly fragmented geographic distribution. Information on its ecology and detailed geographic distribution are scarce and may negatively impact on its conservation status. Therefore, a literature review on its taxonomy, morphology, ecology, and conservation is presented here. Moreover, field records from the authors, citizen science and the existing literature provide an updated geographic distribution highlighting its presence within 13 new and 7 historic 10 × 10 km cells. Bioclimatic niche modelling was then applied to explore patterns of habitat suitability and phenotypic variation within N. h. cetti. The geographic distribution of the species was found to be positively correlated with altitude and precipitation values, whereas temperature showed a negative correlation. Taken together, these outcomes may explain the snake's presence, particularly in eastern Sardinia. In addition, analysis of distribution overlap with the competing viperine snake (N. maura) and the urodeles as possible overlooked trophic resources (Speleomantes spp. and Euproctus platycephalus) showed overlaps of 66% and 79%, respectively. Finally, geographical or bioclimatic correlations did not explain phenotypic variation patterns observed in this highly polymorphic taxon. Perspectives on future research to investigate N. h. cetti's decline and support effective conservation measures are discussed.

Genomics reveals broad hybridization in deeply divergent Palearctic grass and water snakes (Natrix spp.)

Understanding speciation is one of the cornerstones of biological diversity research. Currently, speciation is often understood as a continuous process of divergence that continues until genetic or other incompatibilities minimize or prevent interbreeding. The Palearctic snake genus Natrix is an ideal group to study speciation, as it comprises taxa representing distinct stages of the speciation process, ranging from widely interbreeding parapatric taxa through parapatric species with very limited gene flow in narrow hybrid zones to widely sympatric species. To understand the evolution of reproductive isolation through time, we have sequenced the genomes of all five species within this genus and two additional subspecies. We used both long-read and short-read methods to sequence and de-novo-assemble two high-quality genomes (Natrix h. helvetica, Natrix n. natrix) to their 1.7 Gb length with a contig N50 of 4.6 Mbp and 1.5 Mbp, respectively, and used these as references to assemble the remaining short-read-based genomes. Our phylogenomic analyses yielded a well-supported dated phylogeny and evidence for a surprisingly complex history of interspecific gene flow, including between widely sympatric species. Furthermore, evidence for gene flow was also found for currently allopatric species pairs. Genetic exchange among these well-defined, distinct, and several million-year-old reptile species emphasizes that speciation and maintenance of species distinctness can occur despite continued genetic exchange.

Erythrism in the Eastern Grass Snake, Natrix natrix (Linnaeus, 1758)

We describe the unusual case of erythrism in the Eastern Grass Snake, Natrix natrix. This colour morph is very rare and has not been reported in the literature before. Despite having observed thousands of N. natrix in the field, we personally detected this morph in only three individuals originating from Slovakia, Romania, and mainland Greece, while photos of a fourth individual from a Greek island were provided to us. In addition, a recent study with a large data set from citizen scientists was unable to produce a single reddish Eastern Grass Snake. Such colouration is likewise uncommon in the western members of Grass Snakes (N. helvetica, N. astreptophora), with two examples provided herein. Because the potential biological importance of erythristic colouration is unclear, we encourage other field herpetologists and naturalists to publish their observations of reddish Grass Snakes in the printed literature.

Temperature and Prey Species Richness Drive the Broad-Scale Distribution of a Generalist Predator

The ongoing climate change and the unprecedented rate of biodiversity loss render the need to accurately project future species distributional patterns more critical than ever. Mounting evidence suggests that not only abiotic factors, but also biotic interactions drive broad-scale distributional patterns. Here, we explored the effect of predator-prey interaction on the predator distribution, using as target species the widespread and generalist grass snake (Natrix natrix). We used ensemble Species Distribution Modeling (SDM) to build a model only with abiotic variables (abiotic model) and a biotic one including prey species richness. Then we projected the future grass snake distribution using a modest emission scenario assuming an unhindered and no dispersal scenario. The two models performed equally well, with temperature and prey species richness emerging as the top drivers of species distribution in the abiotic and biotic models, respectively. In the future, a severe range contraction is anticipated in the case of no dispersal, a likely possibility as reptiles are poor dispersers. If the species can disperse freely, an improbable scenario due to habitat loss and fragmentation, it will lose part of its contemporary distribution, but it will expand northwards.

New Evidence on the Distribution of the Highly Endangered Natrix natrix cypriaca and Implications for Its Conservation

The critically endangered Cyprus grass snake (Natrix natrix cypriaca) has been studied for the past 25 years. Although the need for a complete survey on the distribution of its population, outside the strict boundaries of the known mountainous localities, had been stressed, such an effort has not been conducted to date. In this study, we used a rapid survey approach to investigate possible sightings of the species upstream of the known distribution in the Troodos Mountains. We are presenting evidence from 13 sightings of the species that expand the previously known distribution in Cyprus by nearly 70%. This is the first time that new localities for the critically endangered N. n. cypriaca have been reported since the rediscovery of the species in 1992 and the extensive work that followed. Almost all new localities were discovered outside of the currently known species distribution, while individuals were found, for the first time, within the Natura 2000 site of Dasos Machaira (CY2000004) with indications of healthy populations in the area. New localities were recorded within watersheds of the Pediaios and Gialias Rivers, the first sightings outside of Serrachis watershed, indicating an even broader distribution of the species in the Troodos region. The importance of sparse springs within systems of ephemeral streams is highlighted as an exceptional niche for the population in the region. We are discussing the importance of our findings for the conservation of the species and propose targeted conservation actions that will highly improve connectivity of the populations in the region. We underline the necessity of expanding the research on this critically endangered subspecies, based on current findings, so as to fully comprehend its ecology and distribution range in the region.

How often do they do it? An in-depth analysis of the hybrid zone of two grass snake species (Natrix astreptophora and Natrix helvetica)

We examined the contact zone of two parapatric species of grass snake (Natrix astreptophora and Natrix helvetica) in southern France. To this end, we used comprehensive sampling, analysed mtDNA sequences and microsatellite loci, and built Species Distribution Models for current and past climatic conditions. The contact zone had established by the mid-Holocene during range expansions from glacial refuges in the Iberian Peninsula (N. astreptophora) and southern or western France (N. helvetica). The contact zone represents a narrow bimodal hybrid zone, with steep genetic transition from one taxon to the other and rare hybridization, supporting species status for N. astreptophora and N. helvetica. Our results suggest that the steepness of the clines is a more robust tool for species delimitation than cline width. In addition, we discovered in western France, beyond the hybrid zone, a remote population of N. helvetica with genetic signatures of hybridization with N. astreptophora, most likely the result of human-mediated long-distance dispersal. For N. helvetica, we identified a southern and a northern population cluster, connected by broad-scale gene flow in a unimodal hybrid zone running across France. This pattern either reflects genetic divergence caused by allopatry in two microrefuges and subsequent secondary contact or introgression of foreign alleles into the southern cluster.

High pathogen prevalence in an amphibian and reptile assemblage at a site with risk factors for dispersal in Galicia, Spain

Ranaviruses are agents of disease, mortality and population declines in ectothermic vertebrates and emergences have been repeatedly linked to human activities. Ranaviruses in the common midwife toad ranavirus lineage are emerging in Europe. They are known to be severe multi-host pathogens of amphibians and can also cause disease in reptiles. Recurrent outbreaks of ranavirus disease and mortality affecting three species have occurred at a small reservoir in north-west Spain but no data were available on occurrence of the pathogen in the other amphibian and reptile species present or at adjacent sites. We sampled nine species of amphibians and reptiles at the reservoir and nearby sites and screened for ranavirus presence using molecular methods. Our results show infection with ranavirus in all nine species, including first reports for Hyla molleri, Pelophylax perezi, Rana iberica, and Podarcis bocagei. We detected ranavirus in all four local sites and confirmed mass mortality incidents involving Lissotriton boscai and Triturus marmoratus were ongoing. The reservoir regularly hosts water sports tournaments and the risks of ranavirus dispersal through the translocation of contaminated equipment are discussed.

Phylogeny of the Eurasian Wren Nannus troglodytes (Aves: Passeriformes: Troglodytidae) reveals deep and complex diversification patterns of Ibero-Maghrebian and Cyrenaican populations

The Mediterranean Basin represents a Global Biodiversity Hotspot where many organisms show high inter- and intraspecific differentiation. Extant phylogeographic patterns of terrestrial circum-Mediterranean faunas were mainly shaped through Pleistocene range shifts and range fragmentations due to retreat into different glacial refugia. Thus, several extant Mediterranean bird species have diversified by surviving glaciations in different hospitable refugia and subsequently expanded their distribution ranges during the Holocene. Such a scenario was also suggested for the Eurasian Wren (Nannus troglodytes) despite the lack of genetic data for most Mediterranean subspecies. Our phylogenetic multi-locus analysis comprised 18 out of 28 currently accepted subspecies of N. troglodytes, including all but one subspecies which are present in the Mediterranean Basin. The resulting phylogenetic reconstruction dated the onset of the entire Holarctic radiation of three Nannus species to the early Pleistocene. In the Eurasian Wren, two North African subspecies represented separate basal lineages from the Maghreb (N. t. kabylorum) and from the Libyan Cyrenaica (N. t. juniperi), being only distantly related to other Mediterranean populations. Although N. troglodytes appeared to be paraphyletic with respect to the Nearctic Winter Wren (N. hiemalis), respective nodes did not receive strong statistical support. In contrast, paraphyly of the Ibero-Maghrebian taxon N. t. kabylorum was strongly supported. Southern Iberian populations of N. t. kabylorum did not clade with Maghrebian populations of the same subspecies but formed a sister clade to a highly diverse European clade (including nominate N. t. troglodytes and eight further taxa). In accordance with a pattern also found in other birds, Eurasian populations were split into a western clade (Europe, Caucasus) and an eastern clade (Central Asia, Sino-Himalayas, East Asia). This complex phylogeographic pattern revealed cryptic diversification in N. troglodytes, especially in the Iberio-Maghrebian region.

Distribution and hybridisation of barred and common grass snakes (Natrix helvetica, N. natrix) in Baden-Württemberg, South-western Germany

The distribution and hybridisation zone of the two grass snake species occurring in the German state of Baden-Württemberg are described, based on genetic data from maternally inherited mitochondrial DNA (mtDNA, up to 1983 bp) and biparentally inherited microsatellite DNA (13 loci). In agreement with previously published morphological evidence, the barred grass snake (Natrix helvetica) occurs in the Upper Rhine Valley and the Black Forest, while the common grass snake (N. natrix, ‘yellow lineage’) is distributed across the remaining, more eastern parts of Baden-Württemberg. Cline analyses across two transects running through the region of Karlsruhe and the Black Forest indicate that the hybrid zone is similarly narrow here as in the previously characterised stretch near Lake Constance. With respect to nuclear DNA, the Black Forest constitutes no impediment to gene flow in comparison with lowland regions (Karlsruhe, Lake Constance). However, on the eastern slope of the Black Forest, the abrupt replacement of mtDNA of N. helvetica by that of N. natrix indicates male-mediated gene flow and that the Black Forest represents a dispersal barrier for female grass snakes.

Large-scale molecular phylogeny, morphology, divergence-time estimation, and the fossil record of advanced caenophidian snakes (Squamata: Serpentes)

Caenophidian snakes include the file snake genus Acrochordus and advanced colubroidean snakes that radiated mainly during the Neogene. Although caenophidian snakes are a well-supported clade, their inferred affinities, based either on molecular or morphological data, remain poorly known or controversial. Here, we provide an expanded molecular phylogenetic analysis of Caenophidia and use three non-parametric measures of support-Shimodaira-Hasegawa-Like test (SHL), Felsentein (FBP) and transfer (TBE) bootstrap measures-to evaluate the robustness of each clade in the molecular tree. That very different alternative support values are common suggests that results based on only one support value should be viewed with caution. Using a scheme to combine support values, we find 20.9% of the 1265 clades comprising the inferred caenophidian tree are unambiguously supported by both SHL and FBP values, while almost 37% are unsupported or ambiguously supported, revealing the substantial extent of phylogenetic problems within Caenophidia. Combined FBP/TBE support values show similar results, while SHL/TBE result in slightly higher combined values. We consider key morphological attributes of colubroidean cranial, vertebral and hemipenial anatomy and provide additional morphological evidence supporting the clades Colubroides, Colubriformes, and Endoglyptodonta. We review and revise the relevant caenophidian fossil record and provide a time-calibrated tree derived from our molecular data to discuss the main cladogenetic events that resulted in present-day patterns of caenophidian diversification. Our results suggest that all extant families of Colubroidea and Elapoidea composing the present-day endoglyptodont fauna originated rapidly within the early Oligocene-between approximately 33 and 28 Mya-following the major terrestrial faunal turnover known as the "Grande Coupure" and associated with the overall climate shift at the Eocene-Oligocene boundary. Our results further suggest that the caenophidian radiation originated within the Caenozoic, with the divergence between Colubroides and Acrochordidae occurring in the early Eocene, at ~ 56 Mya.

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.

Complex hybridization patterns in European pond turtles (Emys orbicularis) in the Pyrenean Region

Hybrid zones are natural laboratories allowing insights in genetic processes like lineage diversification, speciation and introgression. Using large sampling, 15 microsatellite loci and a mitochondrial marker, we examined the Pyrenean contact zone of three pond turtle taxa (Emys orbicularis orbicularis, E. o. galloitalica, E. o. occidentalis). The Pyrenees are a biogeographically important region separating many lineages endemic to the Iberian Peninsula from their Western European counterparts. We discovered limited admixture, reflecting a complex biogeographic scenario. Simulations using Approximate Bayesian Computing supported that E. o. orbicularis invaded the Iberian Peninsula in the Holocene, circumventing the Pyrenees along the Mediterranean coast, and hybridized in the northern peninsula with the local coastal subspecies galloitalica, and to a lesser extent, with occidentalis. While E. o. occidentalis, and in particular E. o. orbicularis, expanded their ranges considerably during Holocene warming, E. o. galloitalica remained largely confined to its former Iberian refuge. Admixture among the three taxa is surprisingly low, and a future taxonomic investigation that includes the unstudied subspecies of E. orbicularis from North Africa, Eastern Europe and Asia has to determine whether their current status properly reflects their evolutionary divergence or whether certain taxa should be regarded as full species.

Extra-Mediterranean glacial refuges in barred and common grass snakes (Natrix helvetica, N. natrix)

Extra-Mediterranean glacial refugia of thermophilic biota, in particular in northern latitudes, are controversial. In the present study we provide genetic evidence for extra-Mediterranean refugia in two species of grass snake. The refuge of a widely distributed western European lineage of the barred grass snake (Natrix helvetica) was most likely located in southern France, outside the classical refuges in the southern European peninsulas. One genetic lineage of the common grass snake (N. natrix), distributed in Scandinavia, Central Europe and the Balkan Peninsula, had two distinct glacial refuges. We show that one was located in the southern Balkan Peninsula. However, Central Europe and Scandinavia were not colonized from there, but from a second refuge in Central Europe. This refuge was located in between the northern ice sheet and the Alpine glaciers of the last glaciation and most likely in a permafrost region. Another co-distributed genetic lineage of N. natrix, now massively hybridizing with the aforementioned lineage, survived the last glaciation in a structured refuge in the southern Balkan Peninsula, according to the idea of 'refugia-within-refugia'. It reached Central Europe only very recently. This study reports for the first time the glacial survival of a thermophilic egg-laying reptile species in Central Europe.

Ecological preferences of Hyla molleri in the colonisation of arboreal refuges in a human-shaped wetland

Hyla molleri is well adapted to arboreal microhabitats, which are used among breeding seasons. This species is common in wetlands across the Iberian Peninsula and is therefore vulnerable to the loss and degradation of these ecosystems. Due to its secretive habits, the knowledge about the ecology of H. molleri, outside the breeding season, is still scarce. Using artificial refuges as a proxy to natural refuges, we studied how H. molleri uses arboreal microhabitats near reproductive areas and which environmental drivers influence refuge colonisation. We installed 70 PVC pipe refuges in isolated trees and tree patches. Pairs of pipes were installed at two different heights and monitored fortnightly, for four consecutive days, for one year. Each sampling day, we registered the values of variables related with seasonality, microhabitat, dominant plant species, weather and site fidelity. We recorded 2234 individual colonization events by H. molleri, including 516 multiple colonization events, with a maximum of nine individuals in a single refuge. Refuges that were colonized before were more likely to be colonized again. Additionally, colonization was lower in spring and summer and higher in tree clusters than in isolated trees. We found no significant differences on the height or temperature of colonized versus non-colonized refuges. Our results highlight the importance of adequate arboreal microhabitats and the need for the conservation of terrestrial habitats around breeding areas. We also show that artificial refuges can be useful for H. molleri and similar species, namely for habitat improvement or the implementation of citizen-science and monitoring programs.

Hybridization patterns in two contact zones of grass snakes reveal a new Central European snake species

Recent studies found major conflicts between traditional taxonomy and genetic differentiation of grass snakes and identified previously unknown secondary contact zones. Until now, little is known about gene flow across these contact zones. Using two mitochondrial markers and 13 microsatellite loci, we examined two contact zones. One, largely corresponding to the Rhine region, involves the western subspecies Natrix natrix helvetica and the eastern subspecies N. n. natrix, whereas in the other, more easterly, contact zone two lineages meet that are currently identified with N. n. natrix and N. n. persa. This second contact zone runs across Central Europe to the southern Balkans. Our analyses reveal that the western contact zone is narrow, with parapatrically distributed mitochondrial lineages and limited, largely unidirectional nuclear gene flow. In contrast, the eastern contact zone is very wide, with massive nuclear admixture and broadly overlapping mitochondrial lineages. In combination with additional lines of evidence (morphology, phylogeny, divergence times), we conclude that these differences reflect different stages in the speciation process and that Natrix helvetica should be regarded as a distinct species. We suggest a nomenclatural framework for presently recognized grass snake taxa and highlight the need for reconciling the conflicts between genetics and taxonomy.