It takes two to tango – Phylogeography, taxonomy and hybridization in grass snakes and dice snakes (Serpentes: Natricidae:

Using two mitochondrial DNA fragments and 13 microsatellite loci, we examined the phylogeographic structure and taxonomy of two codistributed snake species (Natrix natrix, N. tessellata) in their eastern distribution area, with a focus on Turkey. We found evidence for frequent interspecific hybridization, previously thought to be extremely rare, and for backcrosses. This underscores that closely related sympatric species should be studied together because otherwise the signal of hybridization will be missed. Furthermore, the phylogeographic patterns of the two species show many parallels, suggestive of a shared biogeographic history. In general, the phylogeographies follow the paradigm of southern richness to northern purity, but the dice snake has some additional lineages in the south and east in regions where grass snakes do not occur. For both species, the Balkan Peninsula and the Caucasus region served as glacial refugia, with several mitochondrial lineages occurring in close proximity. Our results show that the mitochondrial divergences in both species match nuclear genomic differentiation. Yet, in the former glacial refugia of grass snakes there are fewer nuclear clusters than mitochondrial lineages, suggesting that Holocene range expansions transformed the glacial hotspots in melting pots where only the mitochondrial lineages persisted, bearing witness of former diversity. On the other hand, the deep mitochondrial divergences in N. tessellata across its entire range indicate that more than one species could be involved, even though lacking microsatellite data outside of Turkey prevent firm conclusions. On the contrary, our microsatellite and mitochondrial data corroborate that N. megalocephala is invalid and not differentiated from sympatric populations of N. natrix. For Cypriot grass snakes, our analyses yielded conflicting results. A critical assessment of the available evidence suggests that N. natrix is a genetically impoverished recent invader on Cyprus and taxonomically not distinct from a subspecies also occurring in western Anatolia and the southern Balkans. Based on combined mitochondrial and nuclear genomic evidence we propose that for grass snakes the following subspecies should be recognized in our study region: (1) Natrix natrix vulgaris Laurenti, 1768, southeastern Central Europe and northern Balkans; (2) Natrix natrix moreoticus (Bedriaga, 1882), southern Balkans, western Anatolia, and Cyprus; and (3) Natrix natrix scutata (Pallas, 1771), eastern Anatolia, Caucasus region, Iran, northeastern distribution range (from eastern Poland and Finland to Kazakhstan and the Lake Baikal region). Thus, Natrix natrix cypriaca (Hecht, 1930) becomes a junior synonym of N. n. moreoticus and Natrix natrix persa (Pallas, 1814) becomes a junior synonym of N. n. scutata. Due to insufficient material, we could not resolve the status of Natrix natrix syriaca (Hecht, 1930) from the Gulf of İskenderun, southeastern Turkey. Vertebrate Zoology 71, 2021, 813–834 | DOI 10.3897/vz.71. e76453 Copyright Marika Asztalos et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Marika Asztalos et al.: Snake tango – Natrix phylogeography, taxonomy and hybridization 814