Next Article in Journal
Two New Species of Elaphoidella (Copepoda, Harpacticoida) from Subterranean Waters in Northeast Thailand, with a Record of a Gynandromorphic Specimen and an Up-to-Date Key to Elaphoidella Species from Southeast Asia
Previous Article in Journal
An Assessment of the Implementation of the Convention on International Trade in Endangered Species of Wild Fauna and Flora in Kenya
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 16
Issue 3
10.3390/d16030184
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phylogeography of Dolichophis Populations in the Aegean Region (Squamata: Colubridae) with Taxonomic Remarks

by
Adam Javorčík
1,
Ilias Strachinis
2,
Evanthia Thanou
3,
Panagiotis Kornilios
3,
Aziz Avcı
4,
Nazan Üzüm
4,
Kurtuluş Olgun
4,
Çetin Ilgaz
5,
Yusuf Kumlutaş
5,
Petros Lymberakis
6,
Zoltán T. Nagy
7 and
Daniel Jablonski
1,*
1
Department of Zoology, Comenius University in Bratislava, Ilkovičova 6, Mlynská Dolina, 842 15 Bratislava, Slovakia
2
Department of Genetics, Development and Molecular Biology, School of Biology, Faculty of Natural Sciences, Aristotle University of Thessaloniki, 541 24 Thessaloniki, Greece
3
Department of Biology, University of Patras, 265 04 Patras, Greece
4
Department of Biology, Faculty of Science, Aydın Adnan Menderes University, 090 10 Efeler, Aydın, Turkey
5
Department of Biology, Faculty of Science, Dokuz Eylül University, 353 90 Buca, İzmir, Turkey
6
Natural History Museum of Crete, School of Sciences and Engineering, University of Crete, Knossos Av, 714 09 Irakleio, Greece
7
Independent Researcher, 13158 Berlin, Germany
*
Author to whom correspondence should be addressed.
Diversity 2024, 16(3), 184; https://0-doi-org.brum.beds.ac.uk/10.3390/d16030184
Submission received: 24 January 2024 / Revised: 6 March 2024 / Accepted: 8 March 2024 / Published: 18 March 2024
(This article belongs to the Special Issue Systematics, Phylogenetics, and Phylogeography of Animals in the Mediterranean Region)
Browse Figures
Review Reports Versions Notes

Abstract

:
In this study, we investigate the phylogeographic patterns of Dolichophis species in the Aegean region, aiming to elucidate their genetic diversity and putative historical colonisation routes through mitochondrial and nuclear DNA data. Our findings revealed distinct phylogeographic patterns: D. caspius exhibited a higher level of haplotypes within two shallow mitochondrial lineages, contrasting with D. jugularis, which displayed lower genetic variability in the area. Additionally, we identified evidence showing possible human-mediated historical translocation of D. caspius populations to Karpathos from the Balkans mainland. The mitochondrial variability in D. jugularis remained relatively uniform across southwestern Anatolia and Dodecanese, except for Rhodes Island. The evidence from mitochondrial and nuclear data confirming the previously described morphological differentiation of the Rhodes snakes, and thus the name D. j. zinneri Cattaneo, 2012, described on the island, could be applied to this isolated population. This result addresses the first genetic view on the long-standing taxonomic uncertainties regarding the subspecies status of Rhodes D. jugularis. Our results also raise questions regarding possible historical hybridisations between D. caspius and D. jugularis in the Dodecanese islands, prompting the need for further investigation using extensive field studies and genomic approaches. Ultimately, the Aegean islands, particularly Kos and Rhodes, seem to be important sites for the evolution of these colubrid snakes and their historical dynamics.

1. Introduction

The Aegean region serves as an outdoor evolutionary laboratory within the Mediterranean biodiversity hotspot, where the present-day biota has undergone remarkable development since the Miocene [1]. Over time, owing to unique local conditions, vertebrates within this region have experienced significant microevolutionary processes and speciation events [2,3,4]. These include hybridisations, specific adaptations, habitat contractions due to climate change, subsequent recolonisations, and competitions resulting from encounters with similar species that have colonised the area from various sources. Additionally, the later human influence played a role in shaping the Aegean biota, particularly through the introduction of species in different locations, notably the numerous islands found in the Aegean [1].
The diversity of snake fauna in the Aegean region is rich (19 species in 15 genera) but still includes many gaps in our knowledge as to distribution, morphology, or genetics. Numerous studies have been related to the historical biogeography of the Aegean herpetofauna, including colubrid snakes [5,6,7] that showed contrasting phylogeographic patterns. However, many snake species have never been genetically studied in the area, and their evolutionary history is only hypothesised (e.g., Eirenis Jan, 1863; Hemorrhois Boie, 1826; Platyceps Blyth, 1860).
The two prominent species of the genus Dolichophis Gistel, 1868, within the region are also less explored. The geographical range of this genus spans from Central Europe (Hungary) across the Balkans and the Aegean to Anatolia, Ukraine, the Caucasus region, the Middle East, and Central Asia [8]. Dolichophis caspius (Gmelin, 1789) and D. jugularis (Linnaeus, 1758) meet in the Aegean region and Anatolia [8,9]. The former has a wide distribution encompassing the Balkans and Anatolia’s mainland, along with numerous Aegean islands. The latter reaches the region via southern Anatolia and can be located on several Dodecanese islands (i.e., east Aegean islands close to Anatolia). However, their phylogeography in this area and the implications of their range overlap, particularly on the small Greek island of Kos, have not been thoroughly investigated. Nonetheless, this does not imply a complete lack of information on the genetic patterns of these two snake species in the region, albeit primarily focused on D. caspius [10,11]. The second species has been studied even less across its entire range [12,13]. The intriguing aspect is that the Aegean region serves as a contact zone for these two snakes and includes numerous islands with varying geological periods of isolation, influencing their phylogeography. Moreover, this area harbours two names linked to the genus Dolichophis that have not undergone genetic testing to date, i.e., Coluber caspius eiselti Zinner, 1972 (with NHMW 18618:1 designated as a holotype), and Dolichophis jugularis zinneri Cattaneo, 2012 (holotype ZFMK 92945). Both taxa come from Rhodes Island [14,15,16,17,18], with the former taxon considered invalid (see [19]) due to its description in an unpublished doctoral thesis of Hermann Zinner [20].
Consequently, we aimed to elucidate these snake populations in the Aegean region using genetic data (mitochondrial and nuclear DNA markers) to delineate their local phylogeography. In addition, we tried to provide pertinent insights into the taxonomy of morphologically distinct populations in Rhodes.

2. Material and Methods

2.1. Sampling

In total, we used 35 new samples of D. caspius and 11 samples of D. jugularis from the Aegean islands (Agios Efstratios, Karpathos, Kos, Lesbos, Limnos, Rhodes, Samothraki, Serifos, Tinos) and from the related mainland Balkans and Anatolia (Türkiye). The material was obtained either directly during field surveys or from collections of the Aydın Adnan Menderes University, Türkiye (ADU); Zoology Department, Ege University, Izmir, Türkiye (ZDEU); Natural History Museum of Crete, Greece (NHMC); and the Zoological Museum University of Patras, Greece (ZMUP). Detailed information about sample codes, sampling localities, and GenBank accession numbers can be found in Table 1.

2.2. DNA Extraction, Molecular Markers, PCR Amplification and Sequencing

Tissue samples were preserved in 96% ethanol. DNA was extracted using the E.Z.N.A.® Tissue DNA Kit (Omega Bio-tek, Inc., Norcross, GA, USA). Sequences of one mitochondrial (mtDNA) and three nuclear (nuDNA) genes were targeted: the mitochondrial protein-coding segment cytochrome b (Cyt b), the nuclear protein-coding genes for the brain-derived neurotrophic factor (BDNF), the oocyte maturation factor Mos (C-mos), and recombination activation gene 1 (Rag1). For PCR, we used Red Taq 2X Master Mix 2 mM MgCl. The primers and PCR conditions used for each marker are presented in Table S1. The same primers were also used for sequencing. PCR products were purified with ExoSAP-ITTM PCR Product Cleanup Reagent (USB Europe GmbH, Staufen, Germany, following manufacturer’s protocol). Sequencing was performed by Macrogen Inc. (Amsterdam, Netherlands; http://www.macrogen.com). We performed a BLAST search in GenBank to confirm that the targeted loci and species were amplified. The translation of protein-coding sequences into amino acids was checked using DnaSP 6.00 [21], confirming the absence of stop codons. Then, we combined new sequences with GenBank data (Table 1).
Table 1. The data on Dolichophis caspius and D. jugularis along with their geographic origins, in addition to GenBank sequences utilised in this study.
Table 1. The data on Dolichophis caspius and D. jugularis along with their geographic origins, in addition to GenBank sequences utilised in this study.
SpeciesCode/Voucher Number (Museum)LocalityNECountryGenBank Number
Cyt b/C-mos/BDNF/Rag1		Source
D. jugularis492 (CUHC)Palio Pyli, Kos Island36.845727.1885GreecePP378899/-/-/-This study
D. jugularis493 (CUHC)Tigaki, Kos Island36.888827.2697GreecePP378900/PP378937/PP378942/PP378947This study
D. jugularis12199 (I. Strachinis)Alyki, Kos Island36.877927.1757GreecePP378926/-/-/-This study
D. jugularis12200 (I. Strachinis)Zia, Kos Island36.852227.2316GreecePP378927/-/-/-This study
D. jugularis12204 (I. Strachinis)Gennadi, Rhodes Island36.029827.9236GreecePP378928/PP378940/PP378941/PP378948This study
D. jugularis10568 (ADU)Kızılot, Antalya Province36.686431.6048TürkiyePP378913/PP378938/PP378944/PP378949This study
D. jugularis10608 (ADU)Bafa Lake, Muğla Province37.485927.5435TürkiyePP378918/-/-/-This study
D. jugularis10621 (ADU)Sabuca Village-Koçarlı, Aydın Province37.762327.6765TürkiyePP378920/-/-/-This study
D. jugularis10627 (ADU)Madran Village-Bozdoğan, Aydın Province37.674428.2866TürkiyePP378921/-/-/-This study
D. jugularis10629 (ADU)Topçam-Çine, Aydın Province37.689128.0148TürkiyePP378922/-/-/-This study
D. jugularis10630 (ADU)Güneyköy-Bozdoğan, Aydın Province37.713628.2114TürkiyePP378923/-/-/-This study
D. jugularisMVZ 2302422 km NE Finike, Antalya Province 36.323830.1672TürkiyeAY486917/AY486941/-[13]
D. jugularisZMHRU 2012/50Kocaaliler, Burdur Province37.316630.6833TürkiyeAY486917/-/-[22]
D. jugularisZMHRU 2012/51Kızılseki, Burdur Province37.266630.7500TürkiyeAY486917/-/-[22]
D. jugularis----Jordan-/AY376798/-[12]
D. caspius216 (ZMUP)Apollonia40.644123.4561GreecePP378890/-/-/-This study
D. caspius231 (ZMUP)Farsala39.256422.3159GreecePP378891/-/-/-This study
D. caspius418 (ZMUP)Tinos Island37.574125.1501GreecePP378892/-/-/-This study
D. caspius427 (ZMUP)Tinos Island37.595725.1535GreecePP378893/-/-/-This study
D. caspius430 (ZMUP)Karpathos Island35.523527.1103GreecePP378894/-/-/-This study
D. caspius431 (ZMUP)Karpathos Island35.564127.0993GreecePP378895/-/-/-This study
D. caspius458 (ZMUP)Serifos Island37.130324.4876GreecePP378896/-/-/-This study
D. caspius459 (ZMUP)Serifos Island37.144024.4485GreecePP378897/-/-/-This study
D. caspius489 (ZMUP)Tigaki, Kos Island36.875627.1853GreecePP378898/PP378936/PP378943/PP378946This study
D. caspius763 (ZMUP)Loutros40.877026.0468GreecePP378901/-/-/-This study
D. caspius6146 (CUHC)Geoponiki40.308023.0940GreecePP378902/-/-/-This study
D. caspius6618 (CUHC)Iztok41.990827.5436BulgariaPP378903/-/-/-This study
D. caspius7761/80.3.35.51 (NHMC)Stratoni40.514123.8238GreecePP378904/-/-/-This study
D. caspius9895 (CUHC)Metamorfosi40.561521.3069GreecePP378905/-/-/-This study
D. caspius9908 (CUHC)Prinos39.590021.6188GreecePP378906/-/-/-This study
D. caspius9921 (CUHC)Deskati39.923021.8209GreecePP378907/-/-/-This study
D. caspius9928 (CUHC)Svoronos40.264522.4487GreecePP378908/-/-/-This study
D. caspius9929 (CUHC)Lofos40.252622.4046GreecePP378909/-/-/-This study
D. caspius10551/231/2013 (ZDEU)Çukurca, Domaniç, Kütahya Province39.790829.6895TürkiyePP378910/-/-/-This study
D. caspius10552/179/2014 (ZDEU)İhsaniye, Eyüb, İstanbul Province41.241128.8048TürkiyePP378911/-/-/-This study
D. caspius10557/232/2013 (ZDEU)Yörükakçayır, Eskişehir Province39.743830.3389TürkiyePP378912/-/-/-This study
D. caspius10573 (ADU)Between Çalı and Atlas, Bursa Province40.150628.9042TürkiyePP378914/-/-/-This study
D. caspius10584 (ADU)Kepez, Aydın Province37.861127.8540TürkiyePP378915/PP378939/PP378945/-This study
D. caspius10595 (ADU)Pamukkale, Denizli Province37.929629.1332TürkiyePP378916/-/-/-This study
D. caspius10601 (ADU)Between Simav and Demirci, Kütahya Province39.124928.8580TürkiyePP378917/-/-/-This study
D. caspius10614 (ADU)Kuşadası, Aydın Province37.686627.1604TürkiyePP378919/-/-/-This study
D. caspius12195 (I. Strachinis)Alykes, Lesbos Island39.226026.2429GreecePP378924/-/-/-This study
D. caspius12197 (I. Strachinis)Antimachia, Kos Island36.814427.1086GreecePP378925/-/-/-This study
D. caspius12205 (I. Strachinis)Agios Efstratios Island39.534824.9940GreecePP378929/-/-/-This study
D. caspius12206 (I. Strachinis)Agios Efstratios Island39.539624.9962GreecePP378930/-/-/-This study
D. caspius12207 (I. Strachinis)Kos Island36.766927.0901GreecePP378931/-/-/-This study
D. caspius12208 (I. Strachinis)Samothraki Island40.508725.5579GreecePP378932/-/-/-This study
D. caspius12209 (I. Strachinis)Samothraki Island40.426525.5311GreecePP378933/-/-/-This study
D. caspius12211 (I. Strachinis)Samothraki Island40.443325.4977GreecePP378934/-/-/-This study
D. caspius12214 (I. Strachinis) Repanidi, Limnos Island39.924225.3147GreecePP378935/-/-/-This study
D. caspiusNHMW KCC1Serifos Island37.160824.4810GreeceAY376739/AY376797/-/-[12]
D. caspius-Andros Island37.846724.8961GreeceAY376739/-/-/-[12]
D. caspius-Şile, Istanbul Province41.167029.5923TürkiyeHM210777/-/-/-[10]
D. caspius-Samos Island37.713526.8166GreeceHM210778/-/-/-[10]
D. caspius-Samos Island37.713526.8166GreeceHM210779/-/-/-[10]
D. caspius-Euboea (Evvoia) Island38.534623.8197GreeceHM210780/-/-/-[10]
D. caspius-Thasos Island40.676824.6439GreeceHM210782/-/-/-[10]
D. caspius-Prespes40.742821.1597GreeceHM210782/-/-/-[10]
D. caspius-Veles41.708121.7759North MacedoniaHM210782/-/-/-[10]
D. caspius-Brest43.275921.7287SerbiaHM210782/-/-/-[10]
D. caspius-Zlot44.011921.9880SerbiaHM210782/-/-/-[10]
D. caspius Kolets, Haskovo41.920725.5469BulgariaHM210783/-/-/-[10]
D. caspius-Shumen43.249126.9427BulgariaHM210783/-/-/-[10]
D. caspius-Hagieni, Constanta43.782928.4842RomaniaHM210783/-/-/-[10]
D. caspius-Tekirdağ, Tekirdağ Province40.984527.5020TürkiyeHM210785/-/-/-[10]
D. caspius-Samothraki Island40.448225.5848GreeceHM210787/-/-/-[10]
D. schmidtiZISP 27777Khumlakh39.234747.0514ArmeniaAY376743/AY376801/-/-[12]
D. schmidtiCAS 182953Mt. Dushak, Mary Province37.130060.0200Turkmenistan-/AY486947/-/-[13]
D. andreanusICSTZM.7H.1154Darreh Shahr, Ilam Province33.144447.3828IranMN531565/-/-/-[23]
E. collaris ZISP 27859Arax Fluss, n. Megri38.885246.2589ArmeniaAY376766/-/-/-[12]
Abbreviations: ADU—Aydın Adnan Menderes University, Aydın, Türkiye; CAS—California Academy of Sciences in San Francisco, California, USA; CUHC—Comenius University Herpetological Collection, Bratislava, Slovakia; ICSTZM—International Center for Science and Technology Zoological Museum, Kerman, Iran; MVZ—Museum of Vertebrate Zoology, Berkeley, California, USA; NHMC—Natural History Museum of Crete, Iraklio, Greece; NHMW—Natural History Museum Vienna, Austria; ZDEU—Zoology Department, Ege University, Izmir, Türkiye; ZISP—Zoological Institute of St. Petersburg, Russia; ZMHRU—Zoology Museum of Harran University, Şanlıurfa, Türkiye; ZMUP—Zoological Museum University of Patras, Greece.

2.3. Sequence Alignment and Data Analysis

The sequences were aligned and edited manually using Geneious Prime 2023.1.2 (Biomatters, Auckland, New Zealand). The default Geneious alignment algorithm was used. The phylogenetic relationships of 56 sequences of D. caspius and D. jugularis together with three outgroup sequences [D. andreanus (Werner, 1917), D. schmidti (Nikolsky, 1909) and E. collaris (Menetries, 1832)] of the final alignment were assessed using maximum likelihood (ML) in RaxML 8.0.0 [24]. The best-fit model of sequence evolution (GTR + G with three subsets) was selected using Partitionfinder 2 [25]. The ML clade support was assessed by 1000 bootstrap pseudoreplicates. This mitochondrial dataset also includes sequences from published studies where the same haplotypes, present in different locations, were used (Table 1).
Mitochondrial haplotype networks were examined and drawn using PopArt [26] (http://www.popart.otago.ac.nz) and the implemented parsimony network algorithm of [27], with a 95% connection limit. Independent networks were considered distinct evolutionarily significant units, following Fraser and Bernatchez [28]. The alignment was 1114 bp long in D. caspius, and 1083 bp long in D. jugularis. DnaSP 6.00 [21] was used to estimate the number of haplotypes (h), haplotype diversity (Hd), nucleotide diversity (π), and uncorrected p distances in mitochondrial datasets.
To access the genealogical relationships between selected taxa/populations of Dolichophis in the Aegean region, we separately conducted an allele network analysis on three nuDNA markers, C-mos, BDNF, and Rag1. For comparison, we used previously generated sequences of C-mos AY376797 (D. caspius); AY376798 and AY486941 (both D. jugularis); AY376801 and AY486947 (D. schmidti) from Nagy et al. [12,13]; and our new sequences from Kos Island (sympatric distribution of D. caspius and D. jugularis; samples 489 and 493), Rhodes (D. j. zinneri; sample 12204), and continental Anatolia (samples 10568 and 10584) (due to different quality of resulting sequences, not all of them were equally represented in nuclear datasets; see Table 1). The resulting sequence alignments (567 bp, 661 bp, and 899 bp, respectively) with more than one heterozygous site were resolved in PHASE 2.1.1 [29] for which the input data were prepared in SeqPHASE [30]. PHASE was run under default settings except for the probability threshold, which was set to 0.9. Allele networks of both analysed markers were examined and drawn using PopArt [26] (http://www.popart.otago.ac.nz) and the implemented parsimony network algorithm of TCS with a 95% connection limit [27].

3. Results

In the current study, we present a phylogeography (mitochondrial and nuclear DNA) of two species within the Dolichophis genus from the Aegean region where they show partially sympatric distribution (Figure 1). Our sampling in the Aegean and surrounding areas is comprehensive, totally encompassing 43 mitochondrial sequences for D. caspius and 12 for D. jugularis, representing 50 and 14 localities, respectively (Table 1, Figure 1). Notably, we show the first genetic insights into Rhodes populations, historically associated with the description of two names within the genus. Additionally, for the first time, we conducted mitotyping on populations across various Aegean islands where the genus is present. These islands include Agios Efstratios, Karpathos, Kos, Limnos, Rhodes, Samothraki, Thassos, and Tinos, representing the Dodecanese, Cyclades (Central Aegean), and North Aegean archipelagos of the Aegean Sea. It is worth noting that our sampling did not cover the Sporades and several small islands in the Dodecanese where the genus is known to occur too (Figure 1, Figure 2 and Figure 3). Furthermore, we also provide a first genetic insight into the populations on Kos Island, the only Aegean island where two species of the genus currently coexist.

3.1. Mitochondrial DNA Data

From our analysed dataset, the maximum likelihood (ML) tree revealed four distinct, well-supported (79–100) clades within the Dolichophis genus, aligning with the current species taxonomy, i.e., D. andreanus, D. caspius, D. jugularis, and D. schmidti (Figure 1B). Our examination of the populations of the two focal species, D. caspius and D. jugularis, in the Aegean area and surrounding regions, unveiled low levels of intraspecific diversity. Notably, the genetic distance between these two species in the area was calculated as 12.0%, with average intraspecific distances of 0.76% for D. caspius and 0.26% for D. jugularis.
Regarding their intraspecific diversity, D. caspius exhibited two well-recognised and well-supported but shallow mitochondrial lineages, with a genetic distance of 1.6% between them (refer to Figure 1B). These two D. caspius lineages are roughly distributed on the east and west sides of the Aegean Sea, respectively. The widely distributed lineage (π = 0.20%) represents populations from the west Aegean and continental Balkans, alongside the islands of Thassos, Samothraki, Limnos, and Karpathos, but also Anatolia. The second lineage, the eastern Aegean–Anatolian (π = 0.14%; denoted by * in Figure 1A,B), is situated on the east side of the Aegean and was recorded in the Black Sea coast, east of the Bosporus, in Bursa, along the Büyük Menderes River in western Anatolia, and on Greek islands, i.e., Kos, Samos, Lesbos, and Agios Efstratios.
The mitochondrial clade representing populations of D. jugularis exhibited no significant phylogenetic substructure. Our sampling encompassed populations from continental western Anatolia (the restricted type locality of the species) and two islands, Kos and Rhodes (the island with the type locality of D. j. zinneri). In examining intraspecific mitochondrial diversity among local populations from both the continent and islands, we identified low distances (but see Rhodes and nuclear DNA data) and only at the haplotype level.
Although the number of mitochondrial haplotypes differs between the studied D. caspius and D. jugularis, with 16 haplotypes and seven haplotypes, respectively, the haplotype variability (Hd) has similar values, i.e., 0.88 and 0.86, respectively. However, it is essential to acknowledge that D. caspius is a widely distributed species in the studied area, with a higher number of analysed samples, which could potentially impact the observed level of haplotypes. Nevertheless, within the dataset of D. caspius, we identified six mitochondrial haplogroups differentiated from each other by a range of two to 13 mutation steps (Figure 2). These haplogroups form a west-to-east division pattern, with the yellow haplogroup being widespread from central Greece to the European part of Türkiye, including islands such as Thassos, Samothraki, and surprisingly Karpathos. Two other groups, red and brown, were detected on the islands of the Cyclades and Euboea (Evvoia), also situated in the west part of the Aegean, while the remaining haplogroups (blue, green, orange) were observed in the east part of the Aegean and Anatolia. The two primary lineages of D. caspius converge with their haplotypes near the Bosporus strait, between Limnos and Agios Efstratios islands, as well as near Samos and Anatolia (Figure 2).
Conversely, the haplotype network of D. jugularis demonstrated a clear and straightforward division into green and yellow haplogroups, as defined here (Figure 3). Mitochondrial sequences were grouped into a green group from continental western and southern Anatolia and Kos Island (with a maximum of two mutation steps between them) and a yellow haplogroup that includes the population on Rhodes Island (representing D. j. zinneri), which was distanced by 10 to 12 mutation steps from other haplotypes of D. jugularis (Figure 3). Notably, populations from the green haplogroup represent the region defined as the type locality of D. jugularis (Figure 3).

3.2. Nuclear DNA Data

The network analysis of C-mos and BDNF identified four and five alleles, respectively. In C-mos, alleles corresponded to the different analysed species (D. caspius and D. schmidti alleles showed one and two mutation steps from D. jugularis, respectively), except D. jugularis, where a heterozygotic sequence from Kos Island (493) was identified (Figure 1C). Sample 12204, which represents the Rhodes population of D. jugularis, shares the main C-mos allele with populations of the species from Kos Island (Greece), mainland Türkiye, and even Jordan. Due to the absence of a D. schmidti BDNF sequence, we compared only D. caspius and D. jugularis using this marker. We identified three alleles for D. caspius (two of them unique to the species) and four alleles for D. jugularis (three unique to the species). The two species share one allele representing geographically close populations from Kos (allele 493a), Rhodes (12204), and western Anatolia (Kepez; 10584a), where both species have come in recent or possible historical contact. The Rhodes population is distanced from the mainland population from Kızılot, Antalya, Türkiye (allele 10568), by three mutation steps and by two from the allele 493b from Kos (Figure 1C). The same allele diversity was recorded in Rag1 (six alleles in two species), where D. jugularis individuals from Rhodes (12204) and Kos (493) were heterozygotic and individual from Kızılot (10568) and D. caspius from Kos (489) were homozygotic. Notably, alleles from Rhodes were genetically equally distant to D. caspius (three to four mutation steps) from Kos and to D. jugularis from the same island (Figure 1C).

4. Discussion

4.1. In Situ Evolution and Natural Dispersal

While snakes belonging to the genus Dolichophis are recognised as common and sizable reptiles that play significant ecological roles [31], scientific attention toward their biogeography and diversification remains relatively limited [10,11,13,20]. This is a surprise considering the taxonomic interest surrounding these snakes that were historically ranked among different genera [32]. Several taxonomic names within snakes belonging to the Dolichophis genus have been described [15,20,33,34], hinting at potential speciation events, yet they have not been thoroughly assessed using contemporary molecular taxonomic methods. This study aims to partially address this research gap in the Aegean area where two species form a contact zone.
Generally, the species-level phylogeography presented here does not bring surprising results compared to the paradigm of phylogeographic patterns and colonisations of the biota and snakes in the Aegean region [1,5]. For both studied Dolichophis species, diversification has occurred relatively recently, much after the initial formation of the Aegean Sea 9–12 Mya [1]. Specifically, the low values of the estimated intraspecific p-distances point to Pleistocene diversifications, as already discussed in Nagy et al. [10] for D. caspius. It is well known that the broader Aegean region is mostly inhabited by D. caspius (sister species to D. schmidti; [23]) with a different level of mitochondrial diversity [10]. Although we lack a well-supported biogeographic hypothesis, the origin of the species is expected to be in the Balkans [11] with subsequent colonisations and divergence through space and time. This species was able to colonise most of the current range in the southern Balkans and western Anatolia, including the Aegean islands and further to the east through the steppes of Eurasia [8]. Surprisingly, from the Balkan–Aegean perspective, the species is not present in the Peloponnese peninsula (Figure 2), which is not well explained in the literature but may be connected to high reptile competition in the peninsula or its specific ecological and topographic conditions. In any case, the diversification history of D. caspius is so recent that the Peloponnese was geologically separated from the remaining Balkan mainland when the species reached this specific sea barrier.
On the other hand, due to the distribution, phylogenetic relationship (sister species to D. andreanus; [23]), and genetic diversity, the ancestral area of D. jugularis is probably in Anatolia or the Middle East from where the species colonised the western areas, reaching the Aegean region on the western coast of Anatolia (the north-western border is probably the Büyük Menderes River; see details in Cattaneo [16]) and several Aegean islands. In this area, the species has contact and possible hybrid zones (see our nuclear evidence Figure 1C) with D. caspius that are so far understudied. Similar distribution patterns where Anatolian–Middle Eastern species reach the south-western coast of Anatolia and neighbouring islands are observed in other reptile genera such as Blanus Wagler, 1830; Anatololacerta Arnold, Arribas, Carranza, 2007; or Hemorrhois (Linnaeus, 1758) [8,35,36]. An interesting parallel example of colonisation and contact zone in the eastern Aegean can be observed between the mostly Balkan species Zamenis situla (Linnaeus, 1758) and the Anatolian–Middle Eastern Z. hohenackeri (Strauch, 1873) [37], although the former species is missing on local islands (but see the discussion regarding Rhodes in Cattaneo et al. [18]).
We can observe more details on the intraspecific level of genetic differentiation for both species. Nagy et al. [10] recognise two main mitochondrial lineages for D. caspius (west and east) that are separated by the Aegean Sea and Bosporus. Using a genomic approach, Mahtani-Williams et al. [11] recognised eight clusters, three of which can be associated with the Aegean region or the Balkans in general. Both studies conclude that this genetic diversity is probably a reflection of former glacial refugia. Our data support the previous conclusions, but detailed sampling on both species better highlights the phylogeographic structure related to the historical development of the Aegean region and island isolations. First, we recognised unique haplotypes in the west side of the Aegean Sea, in the Cyclades islands, and Euboea that were not found in other parts of the studied area. Additionally, our results revealed that D. caspius dispersed from the Balkan Peninsula to Anatolia at least twice, since we found a previously undiscovered mitochondrial haplogroup in Anatolia that is nested in the west lineage. This has come into secondary contact with the east lineage in west Türkiye and some of the adjoining islands. We can assume that these two mitochondrial clades introgress, which would explain the population structure patterns based on the genome-wide SNPs of Mahtani-Williams et al. [11]. This pattern of multiple dispersals can be explained by the wider and faster colonisation to the east that was conducted by the most widespread lineage containing red, brown, orange, and yellow haplogroups (Figure 2). On the other hand, the lineage of D. caspius including blue and green haplotypes (Figure 2) colonised areas in western Anatolia and eastern Aegean islands, where it displays certain levels of local genetic diversity (see unique haplotypes for Kos or Agios Efstratios islands; Figure 2).
Unique mitochondrial haplotypes from islands (Kos, Rhodes) were also found in D. jugularis; however, the overall genetic variation in this species is lower in the Aegean region (but see D. j. zinneri below). Nevertheless, these data represent the first broader genetic insight into the phylogeography of D. jugularis. For further conclusions, more genetic data from the rest of the species range are needed.

4.2. Human Mediated Introductions

Unexpectedly, our genetic analysis revealed an intriguing result for the D. caspius populations from Karpathos Island. Mitochondrial sequences from this island exhibited a close genetic affinity with the haplotype observed in the continental Balkans, Thassos, and Limnos islands. This genetic correspondence suggests that the Karpathos population might have an artificial origin, possibly resulting from historical introductions. This scenario mirrors past instances of human-mediated and accidental amphibian and reptile introductions in the Aegean region and the broader Mediterranean area, a phenomenon documented in various snake species such as Eryx jaculus (Linnaeus, 1758), Telescopus fallax Fleischmann, 1831, Hierophis viridiflavus (Lacépède, 1789) (under the name Coluber gemonensis gyaronensis Mertens, 1968), and Zamenis situla (Linnaeus, 1758) [38,39,40,41]. Prior to the utilisation of molecular data, such introductions, particularly concerning snakes, remained largely unnoticed and sometimes led to wrong taxonomic conclusions [42] or uncertain biogeographic hypotheses. Exceptions to this trend are introductions that have been historically well documented [43].
The genus Dolichophis partly exemplifies a well-documented case of introductions. In 1960, populations of this genus were deliberately introduced to Cyprus from continental Anatolia to mitigate rodent populations on the island [44,45]. Although contentious, such introductions represent contemporary approaches to addressing agricultural pest-related issues. Similarly, populations of D. caspius in the Dalmatian islands of Croatia are believed to have an introduced origin [46], an inference substantiated by genomic investigations revealing a human-mediated origin [11]. These instances underscore the complex interplay between human intervention and wildlife distribution, especially in island ecosystems characterised by high levels of endemism, such as Cyprus and Karpathos. This is why we must carefully evaluate the local biota, at least by population genotyping. The phylogeographic pattern of other populations does not suggest further historical introductions, although populations of D. caspius in Anatolia display certain mixed patterns of mitochondrial data that can be considered as natural dispersion in several waves or a combination with human-mediated introduction (see green haplotype at the Black Sea coast; Figure 2).

4.3. Dolichophis jugularis and Island Endemism

Rhodes Island, despite its size and isolation from the Anatolian mainland, is not generally regarded as a sole source of endemism. While the island does harbour certain locally endemic herpetofauna [e.g., Pelophylax cerigensis (Beerli, Hotz, Tunner, Heppich, Uzzell, 1994)], these species are not exclusive to Rhodes and in some cases, they have limited or mixed patterns of distribution there. These endemics have centres of distribution in closely neighbouring islands within the Dodecanese, notably Karpathos. Due to its prolonged isolation within the Aegean region [47], Karpathos is considered a notable example of local endemism, hosting species such as Lyciasalamandra helverseni (Pieper, 1963) or Mediodactylus oertzeni (Boettger, 1888) found throughout the whole island.
Consequently, the geographic distance of Rhodes from the mainland (20 km), its strong chorotype affinity to Anatolian herpetofauna, and the inherent complexities in diagnosis and taxonomy have historically led to uncertainties regarding two names related to the genus Dolichophis originating in Rhodes: Coluber caspius eiselti Zinner, 1972 (specimens NHMW 18618:1, 3; Figure 4 and Figure 5 in this study), and D. jugularis zinneri Cattaneo, 2012 (holotype ZFMK 92945; see [15] (Figure 1 and Figure 2) and Figure 6 in this study). The former is considered invalid due to the unpublished status of the Zinner thesis (see the discussion in [14,19,20]), while the latter’s validity is in doubt due to the aforementioned reasons and the lack of available genetic data from the island. The taxonomic classifications by Zinner [20] and Cattaneo [15] primarily relied on morphological characteristics that can be somewhat plastic, particularly concerning patterns, colouration, and observed behavioural differences in Rhodes populations in comparison to others. Dolichophis j. zinneri stands out for several distinct morphological features compared to its counterparts. Specifically, it exhibits a wider head with a shorter and pointed shape and consistently displays yellow colouration in certain areas such as the supralabials, neck sides, and throat, contrasting with the red salmon hues observed in D. j. jugularis. Additionally, according to Cattaneo [15], this subspecies displays less intense dorsal dark colouration, which tends to darken anterior–posteriorly. Cattaneo [15] also mentions a high proportion of individuals displaying a yellow-bellied phenotype, a trait less common in the nominate subspecies, which typically exhibits red-bellied phenotypes with a lower count of ventral scales (usually fewer than 200). Furthermore, the ventral colouration—either red or yellow but never black—serves as a distinguishing characteristic between D. j. zinneri and other populations. These characteristics can be observed in specimens known under the name C. caspius eiselti (NHMW 18618:1, 3; Figure 4 and Figure 5). According to Cattaneo [15,16], D. j. zinneri is reported to be present on the islands of Halki, Rhodes, Tilos, and Symi, which needs to be confirmed by further genetic screening.
Despite our sequencing efforts being limited to just one specimen from Rhodes, our study provides the first genetic evidence to address this taxonomic challenge. The mitochondrial distinction (although not strong; Figure 1 and Figure 2) and the character of networks of D. jugularis from Rhodes, particularly within nuclear DNA markers, suggest a genetic affinity to the island. The allele difference (alongside the geographic isolation) between the Kızılot population of D. jugularis in western Türkiye (representing a population of nominotypical subspecies; Figure 3) and Rhodes Island indicates that the subspecies status of D. j. zinneri is likely justified. This suggests independent evolution within Rhodes rather than recent colonisation. Our findings may thus help to clarify more than 50 years of uncertainty regarding the distinct phenotype of D. jugularis from Rhodes. Based on the information provided in our results, as well as considering morphology, and noting that a name for the Rhodes population already exists, the designation D. j. zinneri can be applied to differentiate these snakes. Nevertheless, we strongly advocate for additional research and more comprehensive data collection across the Dodecanese to obtain further evidence supporting or refuting subspecies status [48].
Nevertheless, the question regarding the similarity in colour and pattern of some Rhodes specimens resembling D. caspius, which led Zinner to describe the new subspecies under this name, remains open. Although we and existing literature do not anticipate the sympatric presence of both species on the island [14,16,17] as observed in Kos, we cannot exclude the possibility that both species coexisted on Rhodes in the past (Pleistocene) when the island was connected with the mainland due to the sea level being lower than it is today [47]. Hypotheses involving remnants of populations or overseas colonisation cannot be disregarded either. Such interactions could have resulted in historical hybridisations and the fixation of specific alleles leading to phenotypes resembling either species or another (see [14]). Also, the presence of shared BDNF alleles within Kos and the differentiation of this allele may support the presumption of the common evolution of both species.
However, the distinct phenotype differentiation of D. jugularis from Rhodes may also signify an environmentally related phenomenon in the Aegean region. Jablonski et al. [4] discuss that specific environmental and ecological conditions, combined with genetic components and selection, might underlie the high phenotype variability (mostly in colour and pattern) observed in island populations. The Aegean region presents several instances where morphologically based taxonomy does not align with genetic data, such as the Aegean Natrix natrix schweizeri Müller, 1932, or Macrovipera schweizeri (Werner, 1935) [4,49,50]. Additionally, Eastern Mediterranean reptile populations often exhibit exceptions to their common colouration and pattern (see [51,52,53]), a trend also observed in Dolichophis populations [22] (genotyped specimens ZMHRU 2012/50-51; Table 1 and Figure 3) from western Anatolia and the western Aegean [54].
Conclusively, while the Rhodes population appears genetically distinct, its morphological differentiation may also result from the specific ecological conditions on the island. However, without genomic data, we also cannot exclude the possibility of historical hybridisations. In this context, Kos Island and, generally, the Dodecanese represent a great, naturally isolated outdoor laboratory for exploring historical and recent hybridisations between both present Dolichophis species in the Aegean region and their ecological divergence. Extensive field studies, along with DNA barcoding and genomics, are imperative to unravel the evolutionary processes in Dolichophis, not solely in the Aegean region but throughout its range.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/d16030184/s1: Table S1. Primers and PCR conditions used in this study; primer orientation (OR): F = forward, R = reverse [55,56,57,58,59].

Author Contributions

Conceptualisation, D.J.; methodology, all authors; analyses, A.J., E.T., P.K., Z.T.N. and D.J.; investigations, all authors; resources, all authors; data curation, all authors; writing—original draft preparation, D.J.; writing—review and editing, all authors; visualisation, A.J., D.J.; funding acquisition, D.J. and P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the grant of the Scientific Grant Agency of the Slovak Republic VEGA 1/0242/21. P.K. was financed by the funding programme “MEDICUS” of the University of Patras.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The new sequence data presented in this study are openly available in GenBank under the accession numbers PP378890-PP378949, respectively.

Acknowledgments

We would like to thank colleagues and friends who donated tissue samples and Jana Poláková and Martina Lawson for their help in the DNA laboratory. For access to the collection of the Natural History Museum in Vienna and the Zoological Research Museum Alexander Koenig in Bonn, we thank Silke Schweiger, Georg Gassner, and Morris Flecks.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Poulakakis, N.; Kapli, P.; Lymberakis, P.; Trichas, A.; Vardinoyiannis, K.; Sfenthourakis, S.; Mylonas, M. A review of phylogeographic analyses of animal taxa from the Aegean and surrounding regions. J. Zool. Syst. Evol. Res. 2015, 53, 18–32. [Google Scholar] [CrossRef]
  2. Dufresnes, C.; Probonas, N.M.; Strachinis, I. A reassessment of the diversity of green toads (Bufotes) in the circum-Aegean region. Integr. Zool 2021, 16, 420–428. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, W.; Feiner, N.; Pinho, C.; While, G.M.; Kaliontzopoulou, A.; Harris, D.J.; Salvi, D.; Uller, T. Extensive introgression and mosaic genomes of Mediterranean endemic lizards. Nat. Commun. 2021, 12, 2762. [Google Scholar] [CrossRef] [PubMed]
  4. Jablonski, D.; Tzoras, E.; Panagiotopoulos, A.; Asztalos, M.; Fritz, U. Genotyping the phenotypic diversity in Aegean Natrix natrix moreotica (Bedriaga, 1882) (Reptilia, Serpentes, Natricidae). ZooKeys 2023, 1169, 87–94. [Google Scholar] [CrossRef]
  5. Kyriazi, P.; Kornilios, P.; Nagy, Z.T.; Poulakakis, N.; Kumlutaş, Y.; Ilgaz, Ҫ.; Avcı, A.; Göçmen, B.; Lymberakis, P. Comparative phylogeography reveals distinct colonization patterns of Cretan snakes. J. Biogeogr. 2013, 40, 1143–1155. [Google Scholar] [CrossRef]
  6. Thanou, E.; Kornilios, P.; Lymberakis, P.; Leaché, A.D. Genomic and mitochondrial evidence of ancient isolations and extreme introgression in the four-lined snake. Curr. Zool. 2020, 66, 99–111. [Google Scholar] [CrossRef]
  7. Asztalos, M.; Ayaz, D.; Bayrakci, Y.; Afsar, M.; Tok, C.V.; Kindler, C.; Jablonski, D.; Fritz, U. It takes two to tango—Phylogeography, taxonomy and hybridization in grass snakes and dice snakes (Serpentes: Natricidae: Natrix natrix, N. tessellata). Vertebr. Zool. 2021, 71, 813–834. [Google Scholar] [CrossRef]
  8. Sindaco, R.; Venchi, A.; Grieco, C. The Reptiles of the Western Palearctic 2. Annotated checklist and Distributional Atlas of the Snakes of Europe, North Africa, the Middle East and Central Asia, with an Update to the Volume 1; Societas Herpetologica Italica: Latina, Italy, 2013; p. 543. [Google Scholar]
  9. Baran, İ.; Avcı, A.; Kumlutaş, Y.; Olgun, K.; Ilgaz, Ç. Türkiye Amfibi ve Sürüngenleri; Palme Yayınevi: Ankara, Turkey, 2021; p. 230. ISBN 978-605-282-611-619. [Google Scholar]
  10. Nagy, Z.T.; Bellaagh, M.; Wink, M.; Paunović, A.; Korsós, Z. Phylogeography of the Caspian whipsnake in Europe with emphasis on the westernmost populations. Amphibia-Reptilia 2010, 31, 455–461. [Google Scholar] [CrossRef]
  11. Mahtani-Williams, S.; Fulton, W.; Desvars-Larrive, A.; Lado, S.; Elbers, J.P.; Halpern, B.; Herczeg, D.; Babocsay, G.; Lauš, B.; Nagy, Z.T.; et al. Landscape Genomics of a Widely Distributed Snake, Dolichophis caspius (Gmelin, 1789) across Eastern Europe and Western Asia. Genes 2020, 11, 1218. [Google Scholar] [CrossRef]
  12. Nagy, Z.T.; Schmidtler, J.F.; Joger, U.; Wink, M. Systematik der Zwergnattern (Reptilia: Colubridae: Eirenis) und verwandter Gruppen anhand von DNA-Sequenzen und morphologischen Daten. Salamandra 2003, 39, 149–168. [Google Scholar]
  13. Nagy, Z.T.; Lawson, R.; Joger, U.; Wink, M. Molecular systematics of racers, whipsnakes and relatives (Reptilia: Colubridae) using mitochondrial and nuclear markers. J. Zool. Syst. Evol. Res. 2004, 42, 223–233. [Google Scholar] [CrossRef]
  14. Bader, T.; Riegler, C.; Grillitsch, H. The herpetofauna of the Island of Rhodes (Dodecanese, Greece). Herpetozoa 2009, 21, 147–169. [Google Scholar]
  15. Cattaneo, A. Il Colubro gola rossa dell’arcipelago di rodi: Dolichophis jugularis zinneri subsp. nova (Reptilia Serpentes). Nat. Sicil. 2012, 36, 77–103. [Google Scholar]
  16. Cattaneo, A. Morpho-ecology of the black whip snake Dolichophis jugularis (Linnaeus, 1758) of the southeast Aegean area (Dodecanese, SW Turkey) including previously unpublished growth rate data (Reptilia Serpentes). Nat. Sicil. 2018, 42, 3–24. [Google Scholar]
  17. Cattaneo, A.; Cattaneo, C.; Grano, M. Update on the herpetofauna of the Dodecanese Archipelago (Greece). Biodivers. J. 2020, 11, 69–84. [Google Scholar] [CrossRef]
  18. Cattaneo, A.; Cattaneo, C.; Grano, M. The herpetofauna of the Dodecanese Islands. In Life on Islands. 2. Zoological Diversity of the Aegean Archipelago. Studies dedicated to Norma Chapmann; Masseti, M., Ed.; Edizioni Danaus: Palermo, Italy, 2023; pp. 37–58. [Google Scholar]
  19. Szczerbak, N.N.; Böhme, W. Coluber caspius—Kaspische Pfeilnatter oder Springnatter. In Handbuch der Reptilien und Amphibien Europas, Band 3/I., Schlangen (Serpentes) I; Böhme, W., Ed.; Aula-Verlag: Wiesbaden, Germany, 1993; pp. 83–96. [Google Scholar]
  20. Zinner, H. Systematics and Evolution of the Species Group Coluber jugularis Linnaeus, 1758—Coluber caspius Gmelin, 1789 (Reptilia, Serpentes). PhD. Thesis, The Hebrew University of Jerusalem, Jerusalem, Israel, 1972; p. 62. [Google Scholar]
  21. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large datasets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  22. Göçmen, B.; Nagy, Z.T.; Ҫiçek, K.; Akman, B. An unusual juvenile coloration of the whip snake (Linnaeus, 1758) observed in Southwestern Anatolia, Turkey. Herpetol. Notes 2015, 8, 531–533. [Google Scholar]
  23. Rajabizadeh, M.; Pyron, R.A.; Nazarov, R.; Poyarkov, N.A.; Adriaens, D.; Herrel, A. Additions to the phylogeny of colubrine snakes in southwestern Asia, with description of a new genus and species (Serpentes: Colubridae: Colubrinae). PeerJ 2020, 8, C9016. [Google Scholar] [CrossRef]
  24. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post analysis of large phylogenies. Bioinformatics 2019, 30, 1312–1313. [Google Scholar] [CrossRef]
  25. Lanfear, R.; Frandsen, P.B.; Wright, A.M.; Senfeld, T.; Calcott, B. Partitionfinder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 2017, 34, 772–773. [Google Scholar] [CrossRef]
  26. Leigh, J.W.; Bryant, D. Popart: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  27. Clement, M.; Posada, D.; Crandall, K.A. TCS: A computer program to estimate gene genealogies. Mol. Ecol. 2000, 9, 1657–1659. [Google Scholar] [CrossRef] [PubMed]
  28. Fraser, D.J.; Bernatchez, L. Adaptive evolutionary conservation: Towards a unified concept for defining conservation units. Mol. Ecol. 2001, 10, 2741–2752. [Google Scholar] [CrossRef] [PubMed]
  29. Stephens, M.; Smith, N.J.; Donnelly, P.A. New statistical method for haplotype reconstruction from population data. Am. J. Hum. Genet. 2001, 68, 978–989. [Google Scholar] [CrossRef]
  30. Flot, J.F. Seqphase: A web tool for interconverting phase input/output fles and fasta sequence alignments. Mol. Ecol. Resour. 2010, 10, 162–166. [Google Scholar] [CrossRef] [PubMed]
  31. Kaşot, N. Field notes on trophic spectrum of Dolichophis jugularis from Northern Cyprus. Biharean Biol. 2016, 10, 153–155. [Google Scholar]
  32. Schätti, B.; Utiger, U. Hemerophis, a new genus for Zamenis socotrae Günther, and a contribution to the phylogeny of Old World racers, whip snakes, and related genera (Reptilia: Squamata: Colubrinae). Rev. Suisse Zool. 2001, 108, 919–948. [Google Scholar] [CrossRef]
  33. Mertens, R.; Wermuth, H. Die Amphibien und Reptilien Europas (Dritte Liste, nach dem Stand vom 1. Januar 1960); Verlag Waldemar Kramer: Frankfurt am Main, France, 1960; p. XI + 264. [Google Scholar]
  34. Afrasiab, S.R.; Mohammad, M.K.; Hussein, A.M. Description of a new species of the genus Dolichophis Gitstel from the upper Mesopotamian Plain-Iraq (Reptilia-Serpentes-Colubridae). J. Biodivers. Environ. Sci. 2016, 8, 15–19. [Google Scholar]
  35. Sindaco, R.; Kornilios, P.; Sacchi, R.; Lymberakis, P. Taxonomic reassessment of Blanus strauchi (Bedriaga, 1884) (Squamata: Amphisbaenia: Blanidae), with the description of a new species from south-east Anatolia (Turkey). Zootaxa 2014, 3795, 311–326. [Google Scholar] [CrossRef]
  36. Karakasi, D.; Ilgaz, Ҫ.; Kumlutaş, K.; Güçlü, Ö.; Kankılıç, T.; Beşer, N.; Sindaco, R.; Lymberakis, P.; Poulakakis, N. More evidence of cryptic diversity in Anatololacerta species complex Arnold, Arribas and Carranza, 2007 (Squamata: Lacertidae) and re-evaluation of its current taxonomy. Amphibia-Reptilia 2021, 42, 201–216. [Google Scholar] [CrossRef]
  37. Hofmann, S.; Mebert, K.; Schulz, K.D.; Helfenberger, N.; Göçmen, B.; Böhme, W. A new subspecies of Zamenis hohenackeri (Strauch, 1873) (Serpentes: Colubridae) based on morphological and molecular data. Zootaxa 2018, 4471, 137–153. [Google Scholar] [CrossRef] [PubMed]
  38. Warnecke, H. Telescopus fallax (Fleischmann, 1831) auf den ozeanischen Strophaden-Inseln? Salamandra 1988, 24, 16–19. [Google Scholar]
  39. Insacco, G.; Spadola, F.; Russotto, S.; Scaravelli, D. Eryx jaculus (Linnaeus, 1758): A new species for the Italian herpetofauna (Squamata: Erycidae). Acta Herpetol. 2015, 10, 149–153. [Google Scholar] [CrossRef]
  40. Chondropoulos, B.P. A checklist of Greek reptiles. II. The snakes. Herpetozoa 1989, 2, 3–36. [Google Scholar]
  41. Tzoras, E.; Foufopoulos, J.; Volz, M.; Troidl, S.; Troidl, A.; Jablonski, D. Lost in the Cyclades: Genetic affiliation of the Yellow-bellied Toad, Bombina variegata (Anura: Bombinatoridae), from Paros Island, Greece. Salamandra 2023, 59, 92–95. [Google Scholar]
  42. Utiger, U.; Schätti, B. Morphology and phylogenetic relationships of the Cyprus racer, Hierophis cypriensis, and the systematic status of Coluber gemonensis gyarosensis Mertens (Reptilia: Squamata: Colubrinae). Rev. Suisse Zool. 2004, 111, 225–238. [Google Scholar] [CrossRef]
  43. Wettstein, O. Herpetologia aegaea. Sitzungsberichte der Akademie der Wissenschaften mathematisch-naturwissenschaftliche Klasse; Springer: Berlin/Heidelberg, Germany, 1953; Volume 162, pp. 651–833. [Google Scholar]
  44. Schätti, B.; Sigg, H. Die Herpetogauna der Insel Zypern. Teil 1: Die herpetologische Erforschung/Amphibien. Herpetofauna 1989, 11, 9–18. [Google Scholar]
  45. Göçmen, B.; Yildiz, M.Z. Snakes and their relations with humans. Kibr. Bilim (Cyprus Sci.) 2006, 2, 28–31. [Google Scholar]
  46. Kletečki, E.; Lanszki, J.; Trócsányi, B.; Mužinić, J.; Purger, J.J. First record of Dolichophis caspius (Gmelin, 1789), (Reptilia: Colubridae) on the island of Olib, Croatia. Nat. Croat. 2009, 18, 437–442. [Google Scholar]
  47. Beerli, P.; Hotz, H.; Uzzell, T. Geologically dated sea barriers calibrate a protein clock for aegean water frogs. Evolution 1996, 50, 1676–1687. [Google Scholar] [CrossRef]
  48. Dufresnes, C.; Poyarkov, N.; Jablonski, D. Acknowledging more biodiversity without more species. Proc. Natl. Acad. Sci. USA 2023, 120, e2302424120. [Google Scholar] [CrossRef] [PubMed]
  49. Kindler, C.; Böhme, W.; Corti, C.; Gvoždík, V.; Jablonski, D.; Jandzik, D.; Metallinou, M.; Široký, R.; Fritz, U. Mitochondrial phylogeography, contact zones and taxonomy of grass snakes (Natrix natrix, N. megalocephala). Zool. Scr. 2013, 42, 458–472. [Google Scholar] [CrossRef]
  50. Freitas, I.; Ursenbacher, S.; Mebert, K.; Zinenko, O.; Schweiger, S.; Wüster, W.; Brito, J.C.; Crnobrnja-Isailović, J.; Halpern, B.; Fahd, S.; et al. Evaluating taxonomic inflation: Towards evidence-based species delimitation in Eurasian vipers (Serpentes: Viperidae). Amphib. -Reptil. 2020, 41, 285–311. [Google Scholar] [CrossRef]
  51. Kalaentzis, K.; Kazilas, C.; Strachinis, I. Two cases of melanism in the Ring-headed Dwarf Snake Eirenis modestus (Martin, 1838) from Kastellorizo, Greece (Serpentes: Colubridae). Herpetol. Notes 2018, 11, 175–178. [Google Scholar]
  52. Schneider, B. Eine melanistische Schlanknatter, Coluber najadum kalymnensis n. subsp. (Colubridae, Serpentes), von der insel Kalymnos (Dodekanes, Ägäis). Bonn. Zool. Beitr. 1979, 30, 380–384. [Google Scholar]
  53. Kazilas, C.; Kalaentzis, K.; Strachinis, I. A case of piebaldism in Anatolian Worm Lizard Blanus strauchi (Bedriaga, 1884) (Squamata: Blanidae) from Kastellorizo, Greece. Herpetol. Notes 2018, 11, 527–529. [Google Scholar] [CrossRef]
  54. Kalogiannis, S. Cases of melanism in Dolichophis caspius (Gmelin, 1789) (Squamata: Colubridae) from Greece and a new distribution record. Parnass. Arch. 2021, 9, 19–22. [Google Scholar]
  55. Burbrink, F.T.; Lawson, R.; Slowinski, J.B. Mitochondrial DNA phylogeography of the North American rat snake (Elaphe obsoleta): A critique of the subspecies concept. Evolution 2000, 54, 2107–2114. [Google Scholar] [CrossRef]
  56. de Queiroz, A.; Lawson, R.; Lemos-Espinal, J.A. Phylogenetic relationships of North American garter snakes (Thamnophis) based on four mitochondrial genes: How much DNA sequence is enough? Mol. Phylogenet. Evol. 2002, 22, 315–329. [Google Scholar] [CrossRef]
  57. Lawson, R.; Slowinski, J.B.; Crother, B.I.; Burbrink, F.T. Phylogeny of the Colubroidea (Serpentes): New evidence from mitochondrial and nuclear genes. Mol. Phylogenet. Evol. 2005, 37, 581–601. [Google Scholar] [CrossRef]
  58. Vieites, D.R.; Min, M.S.; Wake, D.B. Rapid diversification and dispersal during periods of global warming by plethodontid salamanders. Proc. Nat. Acad. Sci. USA 2007, 104, 19903–19907. [Google Scholar] [CrossRef] [PubMed]
  59. Growth, J.G.; Barrowclough, G.F. Basal Divergences in Birds and the Phylogenetic Utility of the Nuclear RAG-1 Gene. Mol. Phylogenet. Evol. 1999, 12, 115–123. [Google Scholar] [CrossRef] [PubMed]
Diversity 16 00184 g001
Figure 1. The genotyped localities of Dolichophis caspius and D. jugularis in the Aegean region (A) together with the mitochondrial tree (B). Colourations on the map and the tree correspond to taxa. Underlined sequence codes in the mtDNA tree mark those used for nuDNA analysis. The orange shading of Rhodes and surrounding islands marks the range of D. jugularis zinneri according to Cattaneo [16]. The asterisk indicates the distribution of the eastern Aegean–Anatolian lineage of D. caspius. The line in the middle of the Aegean Sea represents the early formation of the Aegean barrier. Nuclear allele networks of the phased sequences of C-mos, BDNF, and Rag1 genes (C) presenting relationships between Rhodes (12204) and selected populations and taxa. Species colours in networks follow those used in the map (A) and tree (B). Circle sizes are proportional to the number of individuals that share a given allele. A small empty circle in networks indicates a missing or hypothetical allele. Different alleles of a single heterozygous sequence is coded as a and b, and without division, they represent homozygous sequences. Refer to Table 1 for locality details and sample codes used in the tree and networks. Inset photographs: Daniel Jablonski (D. caspius from Bulgaria) and Ilias Strachinis (D. jugularis from Rhodes). The map was generated using QGIS 3.28, available at https://qgis.org/.
Figure 1. The genotyped localities of Dolichophis caspius and D. jugularis in the Aegean region (A) together with the mitochondrial tree (B). Colourations on the map and the tree correspond to taxa. Underlined sequence codes in the mtDNA tree mark those used for nuDNA analysis. The orange shading of Rhodes and surrounding islands marks the range of D. jugularis zinneri according to Cattaneo [16]. The asterisk indicates the distribution of the eastern Aegean–Anatolian lineage of D. caspius. The line in the middle of the Aegean Sea represents the early formation of the Aegean barrier. Nuclear allele networks of the phased sequences of C-mos, BDNF, and Rag1 genes (C) presenting relationships between Rhodes (12204) and selected populations and taxa. Species colours in networks follow those used in the map (A) and tree (B). Circle sizes are proportional to the number of individuals that share a given allele. A small empty circle in networks indicates a missing or hypothetical allele. Different alleles of a single heterozygous sequence is coded as a and b, and without division, they represent homozygous sequences. Refer to Table 1 for locality details and sample codes used in the tree and networks. Inset photographs: Daniel Jablonski (D. caspius from Bulgaria) and Ilias Strachinis (D. jugularis from Rhodes). The map was generated using QGIS 3.28, available at https://qgis.org/.
Diversity 16 00184 g001
Diversity 16 00184 g002
Figure 2. The distribution of mitochondrial haplotypes of Dolichophis caspius in the Aegean area (A) and the haplotype network (B) defining haplogroups by colouration. The names in bold of the haplotype network highlight the Aegean islands. Refer to Table 1 for locality details and sample codes used in the network. The line in the middle of the Aegean Sea represents the early formation of the Aegean barrier. The question mark indicates a place of unclear or possible occurrence of the species on Rhodes Island. The distribution range of the species in the studied area is highlighted in orange. The map was generated using QGIS 3.28, available at https://qgis.org/.
Figure 2. The distribution of mitochondrial haplotypes of Dolichophis caspius in the Aegean area (A) and the haplotype network (B) defining haplogroups by colouration. The names in bold of the haplotype network highlight the Aegean islands. Refer to Table 1 for locality details and sample codes used in the network. The line in the middle of the Aegean Sea represents the early formation of the Aegean barrier. The question mark indicates a place of unclear or possible occurrence of the species on Rhodes Island. The distribution range of the species in the studied area is highlighted in orange. The map was generated using QGIS 3.28, available at https://qgis.org/.
Diversity 16 00184 g002
Diversity 16 00184 g003
Figure 3. The distribution of mitochondrial haplotypes of Dolichophis jugularis in the Aegean area (A) and the haplotype network (B) defining haplogroups by colouration. The names in bold of the haplotype network highlight the Aegean islands. Refer to Table 1 for locality details and sample codes used in the network. The question mark indicates a place of possible occurrence of the species in the westernmost Anatolia [8]. The distribution range of the species in the studied area is highlighted in orange. The map was generated using QGIS 3.28, available at https://qgis.org/.
Figure 3. The distribution of mitochondrial haplotypes of Dolichophis jugularis in the Aegean area (A) and the haplotype network (B) defining haplogroups by colouration. The names in bold of the haplotype network highlight the Aegean islands. Refer to Table 1 for locality details and sample codes used in the network. The question mark indicates a place of possible occurrence of the species in the westernmost Anatolia [8]. The distribution range of the species in the studied area is highlighted in orange. The map was generated using QGIS 3.28, available at https://qgis.org/.
Diversity 16 00184 g003
Diversity 16 00184 g004
Figure 4. The adult specimen NHMW 18618:1 (holotype; [20]) representing Coluber caspius eiselti Zinner 1972 from Lindos, Rhodes, Greece. (A,B) Dorsal and ventral view of the body. (C–F) Dorsal, ventral, and lateral view of the head. (G–I) Dorsal and ventral view of the body’s scalation, colouration, and pattern. Photos: Daniel Jablonski.
Figure 4. The adult specimen NHMW 18618:1 (holotype; [20]) representing Coluber caspius eiselti Zinner 1972 from Lindos, Rhodes, Greece. (A,B) Dorsal and ventral view of the body. (C–F) Dorsal, ventral, and lateral view of the head. (G–I) Dorsal and ventral view of the body’s scalation, colouration, and pattern. Photos: Daniel Jablonski.
Diversity 16 00184 g004
Diversity 16 00184 g005
Figure 5. The juvenile specimen NHMW 18618:3 (paratype; [20]) representing Coluber caspius eiselti Zinner 1972 from Lindos, Rhodes, Greece. (A,B) Dorsal and ventral view of the body. (C–F) Dorsal, ventral, and lateral view of the head. (G–I) Dorsal and ventral view of the body’s scalation, colouration, and pattern. Photos: Daniel Jablonski.
Figure 5. The juvenile specimen NHMW 18618:3 (paratype; [20]) representing Coluber caspius eiselti Zinner 1972 from Lindos, Rhodes, Greece. (A,B) Dorsal and ventral view of the body. (C–F) Dorsal, ventral, and lateral view of the head. (G–I) Dorsal and ventral view of the body’s scalation, colouration, and pattern. Photos: Daniel Jablonski.
Diversity 16 00184 g005
Diversity 16 00184 g006
Figure 6. The adult specimen ZFMK 92945 (holotype; [15]) from the type collection of Dolichophis jugularis zinneri Cattaneo, 2012. (A,B) Dorsal and ventral view of the body. (C–F) Dorsal, ventral, and lateral view of the head. (G–I) Dorsal and ventral view of the body’s scalation, colouration, and pattern. Photos: Morris Flecks (A–D) and Daniel Jablonski (E–I).
Figure 6. The adult specimen ZFMK 92945 (holotype; [15]) from the type collection of Dolichophis jugularis zinneri Cattaneo, 2012. (A,B) Dorsal and ventral view of the body. (C–F) Dorsal, ventral, and lateral view of the head. (G–I) Dorsal and ventral view of the body’s scalation, colouration, and pattern. Photos: Morris Flecks (A–D) and Daniel Jablonski (E–I).
Diversity 16 00184 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Javorčík, A.; Strachinis, I.; Thanou, E.; Kornilios, P.; Avcı, A.; Üzüm, N.; Olgun, K.; Ilgaz, Ç.; Kumlutaş, Y.; Lymberakis, P.; et al. Phylogeography of Dolichophis Populations in the Aegean Region (Squamata: Colubridae) with Taxonomic Remarks. Diversity 2024, 16, 184. https://0-doi-org.brum.beds.ac.uk/10.3390/d16030184

AMA Style

Javorčík A, Strachinis I, Thanou E, Kornilios P, Avcı A, Üzüm N, Olgun K, Ilgaz Ç, Kumlutaş Y, Lymberakis P, et al. Phylogeography of Dolichophis Populations in the Aegean Region (Squamata: Colubridae) with Taxonomic Remarks. Diversity. 2024; 16(3):184. https://0-doi-org.brum.beds.ac.uk/10.3390/d16030184

Chicago/Turabian Style

Javorčík, Adam, Ilias Strachinis, Evanthia Thanou, Panagiotis Kornilios, Aziz Avcı, Nazan Üzüm, Kurtuluş Olgun, Çetin Ilgaz, Yusuf Kumlutaş, Petros Lymberakis, and et al. 2024. "Phylogeography of Dolichophis Populations in the Aegean Region (Squamata: Colubridae) with Taxonomic Remarks" Diversity 16, no. 3: 184. https://0-doi-org.brum.beds.ac.uk/10.3390/d16030184

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop