Acta Herpetologica 13(2): 117-123, 2018 ISSN 1827-9635 (print) © Firenze University Press ISSN 1827-9643 (online) www.fupress.com/ah DOI: 10.13128/Acta_Herpetol-23189 New insights on the phylogenetic position and population genetic structure of the Critically Endangered Karpathos marsh frog Pelophylax cerigensis (Amphibia: Anura: Ranidae) Elisavet A. Toli1,*, Souzanna Siarabi1, Anastasios Bounas1, Panayiotis Pafilis2, Petros Lymberakis3, Kon- stantinos Sotiropoulos1 1 Molecular Ecology & Conservation Genetics Lab, Department of Biological Applications & Technology, University of Ioannina, 45110, Ioannina, Greece. *Corresponding author. E-mail etoli@cc.uoi.gr 2 Section of Zoology and Marine Biology, Department of Biology, National and Kapodistrian University of Athens, 15784, Athens, Greece 3 Natural History Museum of Crete, University of Crete, 71409, Irakleio, Greece Submitted on: 2018, 24th April; revised on: 2018, 9th August; accepted on: 2018, 26th August Editor: Marcello Mezzasalma Abstract. The genus Pelophylax, which currently comprises 26 species, is a well studied group due to its complex his- tory and high diversification, although some phylogenies remain unresolved. Here we assess the phylogenetic position and the population genetic structure of the Critically Endangered Karpathos frog, Pelophylax cerigensis, endemic to Karpathos Island. A total of 42 cytb sequences were examined including specimens from Rhodes Island, and ampli- fied fragment length polymorphisms (AFLPs) were generated to investigate the genetic structure and connectivity of the only two known populations on Karpathos. Molecular analyses reveal two major waterfrog lineages in the eastern Mediterranean: Clade A comprises Pelophylax bedriagae from the Middle East and the island of Cyprus, while clade B includes both P. bedriagae from the east Aegean Sea and P. cerigensis. Specimens from Karpathos and Rhodes consti- tute a single clade, contrasting previous studies, thus indicating the occurrence of Karpathos frog also to the neighbor- ing Rhodes. The AFLP markers revealed low but statistically significant levels of genetic divergence between the two Karpathos’ populations and similar levels of genetic diversity. Our results suggest that the current taxonomy of the species should be re-evaluated. We also strongly recommend the need of conservation actions to maintain the levels of diversity in the declining population of the Karpathos frog. Keywords. Pelophylax cerigensis, insular populations, phylogeny, cytochrome b, genetic structure, AFLPs, conserva- tion. INTRODUCTION The elucidation of phylogeographic history and the study of intra- and inter-population diversity are inextri- cably linked to the paleogeographic events that took place during the Pleistocene. The Balkan Peninsula, which is considered a Pleistocenic refugium, shows high levels of biodiversity and endemism (Sfenthourakis et al., 2001; Michaux et al., 2004; Sotiropoulos et al., 2007) and the long-term favorable environmental conditions allowed the preservation of genetically differentiated populations (Canestrelli et al., 2010) resulting in the complex history of many terrestrial species (Beerli et al., 1996; Douris et al., 1998; Poulakakis et al., 2003). The Genus Pelophylax Fitzinger, 1843, which cur- rently comprises 26 species, is one of the most well stud- ied groups (Lymberakis et al., 2007; Plötner et al., 2012; Vervust et al., 2013; Litvinchuk et al., 2015) due to its 118 Elisavet A. Toli et alii complex history and high diversification in the Balkan Peninsula and east Mediterranean. Isolation of land areas during the Pleistocene and the formation of islands in the Aegean region had a great impact in the genetic differ- entiation of the genus Pelophylax. Currently, researchers suggest long-separated lineages and diverse evolutionary histories of the western Palearctic water frogs (Lymbe- rakis et al., 2007; Akın et al., 2010), although some phy- logenies are still under question and further revision of their status is needed. Among the three major lineages in the eastern Mediterranean region, the Balkan-Anatolia (ridibunda/ bedriagae) lineage which comprises the spe- cies P. cretensis, P. epeiroticus, P. bedriagae, P. cerigensis, P. kurtmuelleri and P. ridibundus, is considered to have emerged from a widely-distributed common ancestor through vicariance events during the climatic changes in the Pliocene (Lymberakis et al., 2007). Pelophylax cerigensis is a medium-sized, mainly insectivore water frog, endemic to Karpathos Island (Valakos et al., 2008; Pafilis et al., 2018). Karpathos is the second largest island of the Dodecanese archipelago (Aegean Sea, Greece), located approximately 47 kilom- eters southwest of Rhodes. The island is characterized by xeric habitats with dry pine and oak forests, Mediter- ranean maquis, and phrygana. According to Beerli et al. (1994), who first described the species, phylogenetic anal- yses based on electrophoretic data cluster the Karpathos and Rhodes populations together. In addition, Akın et al. (2010) found a shared haplotype in both Karpathos and Rhodes, a result further corroborated by Lymberakis et al. (2007). However, the taxonomy of the Rhodes popula- tions remains unresolved and further evidence is needed for the clarification of the species status. The Karpathos frog is regarded as the most threat- ened frog in Europe (Beerli et al., 2009; Temple and Cox, 2009) due to its limited distribution and current popula- tion decline. Habitat loss and degradation of aquatic ter- rains in this dry island are considered the major factors compromising the species survival whereas they could potentially affect the populations’ genetic status. There- fore, studies on population genetic diversity and structure are necessary to design and implement proper manage- ment actions for the conservation of the species. Due to their simplicity and cost-effectiveness as genetic markers, AFLPs have been widely used in popu- lation genetic studies in the past decades. Several stud- ies have examined the genetic structure through AFLP markers using F-statistics (Wright, 1950; Meudt et al., 2007) and quantified the levels of differentiation and con- nectivity between populations of several species. Here, we aim to evaluate the phylogenetic position of the Karpathos frog within the bedriagae lineage as well as to assess the genetic structure and connectivity of the population between the two known localities in the island of Karpathos. MATERIAL AND METHODS Sample collection and laboratory procedures A total of 30 buccal swab samples of P. cerigensis were col- lected during 2016 from the two known breeding sites (Argoni and Nati rivers) on Karpathos (Table A1, Fig. 1). Additionally, three samples of water frogs (presumed P. bedriagae) from Rhodes Island (NHMC80.2.99.40, NHMC80.2.99.41, NHMC80.2.99.42) were used in the phylogenetic study. Total genomic DNA was extracted from buccal swabs using the NucleoSpin Tissue kit (Macherey-Nagel) following the manufacturer’s protocol. A partial sequence of approximately 340 bp from the mito- chondrial cytb gene was successfully amplified for 26 samples (including the three from Rhodes, Table A1), with the universal primers L14841 and H15149 (Kocher et al., 1989). PCR amplifi- cations were carried out in 12.5μl volume reactions containing 0.05 U Taq (Kapa Biosystems), 1.5 mM MgCl2, 0.2 mM dNT- Ps, 0.4 μΜ of each primer, 1X Taq buffer (Kapa Biosystems) and 10-20 ng DNA template, under the following conditions: an initial denaturation step at 94 °C for 5 min, followed by 35 cycles of denaturation at 94 °C for 60 sec, annealing at 47 °C for 60 sec, extension at 72 °C for 60 sec and a final extension for 7 min at 72 °C. PCR products were purified using the Nucle- oSpin ExtractII (Macherey-Nagel) cleanup kit and single strand sequencing was conducted by CEMIA (Cellular & Molecular Immunological Applications, Larisa, Greece). AFLP analysis was performed according to a modified protocol from Vos et al. (1995), for 28 samples of high qual- ity DNA, collected from the two breeding sites on Karpathos, Argoni (18 samples) and Nati (10 samples). Total genomic DNA was digested using two restriction enzymes TaqI and EcoRI (Takara Bio Inc.). Digestion was carried out in a final volume of 20 μl containing 1 mM TaqI buffer, 0.1 mM BSA, 0.5U of EcoRI and 200 ng of genomic DNA for 2 h at 37 °C. After the addition of 0.5U TaqI samples were incubated for two more hours at 65 °C. The ligation was carried out in a final volume of 30 μl con- taining 1X ligase buffer (Takara), 1.8 μM of each adaptor, 1U T4 DNA ligase and the digested DNA, and the reaction was incu- bated overnight at 16 °C. The digested-ligated DNA fragments were diluted 25-fold to be used as templates for the pre-ampli- fication reaction. Each 50 µl pre-selective reaction contained 1X Taq buffer, 1.5 mM MgCl2, 0.2 mM of each dNTP, 0.3 μΜ of each primer (T01P2, T02P2), 0.3 μΜ of Presel ECO primer and 10 μl of the diluted ligation product. The PCR amplifications were carried out using the following profile: an initial denatura- tion step at 94 °C for 120 sec followed by 20 cycles of 30 sec at 94 °C, 60 sec at 56 °C, and 60 sec at 72 °C with a final extension step at 72 °C for 300 sec. The pre-amplification products were diluted 25-fold to be used as template for the selective ampli- fication. The selective amplifications were performed in a total volume 25 μl containing 10 mM Taq buffer (Takara), 3 mM 119Phylogeny and population structure of the Karpathos frog MgCl2, 0.3 mM of each dNTP, 0.2 μM of EcoRI primers, 0.5 μM of 5 selective primers (T101P2, T105P2, T106P2, T204P2 & T205P2) and 5 μl of diluted pre-amplified DNA. Selective amplification was carried out using a touchdown protocol with an initial denaturation step at 94 °C for 120 sec, 30 sec at 94 °C, 60 sec at 65 °C followed by 11 cycles where the annealing temperature was gradually reduced 0.7°C per cycle followed by 23 cycles of 30 sec at 94 °C, 30 sec at 56°C, and 60 sec at 72 °C. The primers used in each pre-selective and selective PCRs are presented in Table A2. Selective products were separated in Fragment Analyzer (Advanced Analytical Technologies Inc.) using the dsDNA 910 Reagent Kit. AFLP patterns were visual- ized and processed using the PROsize 2.0 Software (Advanced Analytical Technologies) simultaneously by two persons to reduce scoring bias. Scoring bias was quantified by calculating Cohen’s kappa coefficient for inter-rater agreement. AFLP pro- files were scored according to presence/absence of peaks. Genetic analyses Mitochondrial sequences were edited by eye in MEGA v.7 (Kumar et al., 2016) and aligned with CLUSTAL W (Thompson et al., 1994). Additionally, 15 sequences of P. bedriagae, which were retrieved from GenBank, and one of P. cretensis that was used as outgroup, were added in the analyses. To visualize the relationships among the detected haplotypes, a median joining (MJ) network (Bandelt et al., 1999) was constructed with the software PopArt (Leigh et al., 2015), excluding the outgroup sequence and setting the parameter ε equal to zero. Addition- ally, a neighbor-joining (NJ) tree (Saitou and Nei, 1987) was constructed according to the Kimura-2 parameter substitution model, implemented in MEGA. The reliability of the nodes was assessed by 50000 bootstrap replications. Between-population uncorrected sequence divergences (p-distance) were estimated using MEGA. Additionally, haplotype and nucleotide diversity values within each recognized clade were calculated with dnaSP v.5 software (Librado and Rojas, 2009). AFLP data were used to estimate Nei’s (1973) gene diver- sity and levels of genetic differentiation between the two popu- lations (Argoni and Nati), using a Bayesian approach with non- uniform prior distribution through 1000 bootstrap replicates implemented in the software AFLP-SURV (Vekemans, 2002) which is based on the methods described by Lynch and Milli- gan (1994). To evaluate the presence of genetic substructure in the population, a Principal Coordinates Analysis (PCoA) was Fig. 1. Sampling localities of P. cerigensis used in the study. Shaded areas correspond to the current distribution of P. bedriagae. 120 Elisavet A. Toli et alii performed in R 3.4.1 (R core team, 2017) to visualize the clus- tering of individuals based on their AFLP band patterns. Indi- viduals with missing data were excluded from the analysis. RESULTS The final dataset included 42 sequences (Table A1) and the analysis of 286 bp of Cytb revealed 12 haplo- types with high levels of haplotypic diversity (Hd= 0.62). We detected 49 variable sites (excluding the outgroup) of which 27 were parsimony-informative. Both the hap- lotype network (Fig. 2) and the NJ phylogenetic tree (Fig. A2) revealed two well-defined lineages, which are in accordance with the geographical origin of the speci- mens. Clade A (5 haplotypes, Hd = 0.89, Pi = 0.014) cor- responds to bedriagae specimens of Syria and Cyprus, while Clade B (7 haplotypes, Hd =0.42, Pi = 0.019) con- sists of the bedriagae specimens from Turkey, the Aegean islands and P. cerigensis. Between clade p-distance was 6%, while p-distances between the various geographical locations ranged from 0.1 to 7.2% (Table 1). Specimens from the islands of Karpathos and Rhodes share one hap- lotype (B-1) while one additional haplotype (differing by one substitution) was found in Rhodes (B-2). In total, 87 AFLP fragments (Fig. A1) were produced from all primer combinations and the mean number of segregating peaks per individual were 59. Cohen’s kap- pa coefficient value was 0.83, suggesting that no scoring bias was present in our dataset. The results are summa- rized in Table 2. The two breeding populations showed low but statistically significant levels of genetic differen- tiation, while levels of gene diversity are similar between them (Table 2). In the PCoA plot, individuals from each breeding population formed two separated groups that widely overlapped (Fig. 3). The first and the second axis accounted for 24% and 15.9% of the total variance, respectively. Fig. 2. Median joining network constructed using the 12 detected haplotypes presenting the two major waterfrog lineages (A, B). Ver- tical lines correspond to the mutational steps observed between the different haplotypes. Table 1. % p-distances (below diagonal) with the corresponding standard error (above diagonal) between the geographical locations of P. bedriagae and P. cerigensis. Karpathos isl. 0.1 1.4 1.3 1.9 1.2 1.2 1.7 1.3 Rhodes isl. 0.1 1.4 1.3 1.9 1.3 1.2 1.8 1.3 Astypalaia isl. 4.6 4.7 0.9 1.7 0.9 0.9 1.5 1 Chios isl. 3.6 3.8 1.8 1.5 0.6 0.6 1.3 0.8 Cyprus 6.5 6.6 4.6 3.7 1.4 1.5 1 1.4 Dadia 3.6 3.8 1.8 0.9 2.7 0.6 1.2 0.8 Lesvos isl. 3.6 3.8 1.8 0.9 3.7 0.9 1.3 0.8 Syria 7.1 7.2 5.2 4.2 2.8 3.3 4.1 1.3 Turkey 4.1 4.2 2.3 1.3 3.2 1.3 1.3 3.9 Table 2. Gene diversity within the two populations of P. cerigensis. Hj: Nei’s gene diversity, S.E.: standard error, P %: level of polymor- phism and fixation index (Fst value). Population Sample size Hj S.E. P % Fst (95% CI) Argoni 18 0.35 0.013 71.2 Nati 10 0.38 0.011 67.8 0.04 (-0.03-0.01) Fig. 3. A two-dimensional PCoA representation of the dissimilari- ties among individuals with 95% confidence ellipses. Points repre- sent sampled individuals from the two breeding sites (Argoni, Nati). 121Phylogeny and population structure of the Karpathos frog DISCUSSION Taxonomic implications Our results reveal two major water frog lineages in the eastern Mediterranean: clade A comprises P. bedria- gae from the Middle East, while clade B includes both P. bedriagae specimens of east Aegean Sea and P. cerigensis specimens. The results suggest homogeneity in respect to mitochondrial DNA for the specimens from Karpathos and Rhodes and reveal one shared haplotype between the Karpathos and Rhodes specimens and a second one in Rhodes Island. The two island populations comprise a single clade, in contrast to previous studies, where Pelophylax specimens of Rhodes clustered with specimens of P. bedriagae from Asia Minor (Lymberakis et al., 2007). This does not seem to be the case in the present study, thus raising an argument concerning the taxonomic validity of the Karpathos frog. However, the results of the above study were underpinned by different genetic mark- ers, validating the need of utilizing both mitochondrial and nuclear markers in each study. The shared haplotype present in both islands in addition to the low levels of genetic divergence, support the presence of P. cerigensis in the neighboring Rhodes, as proposed by Akin et al. (2010). However, according to geological data Karpathos and Rhodes were already completely isolated during the Pleistocene (Beerli et al., 1996). Although it is not clear when the two islands were separated from Asia Minor, geological records propose that Karpathos was isolated earlier while Rhodes was still connected to Anatolia until late Pliocene or early Pleisto- cene (Boger and Dermitzakis, 1985; Beerli et al., 1994). This disagreement with the paleogeographic history should be investigated using nuclear data from several species along with geological information, which could unravel the complex biogeographic history of the Aegean archipelago. The high pairwise p-distances within clade B (3.6-4.7%) corroborate the divergent and complicated history of the genus Pelophylax, which is linked with the events during the Pleistocene. Our results confirm the necessity of revision of the current taxonomy of the species by utilizing both mito- chondrial and nuclear markers. In addition, further anal- ysis of bioacoustic data could possibly shed light on the subject as they are proven to play an important role in anuran taxonomy (Schneider and Sinsch, 1999; Padial et al., 2008). Additionally, the use of more specimens cov- ering the species’ distribution is highly recommended to verify its taxonomic validity. Population structure and implications for conservation According to the population genetic analysis, both Karpathos populations show moderate levels of genetic variation and differentiation. The AFLP markers showed low but statistically significant levels of genetic divergence between the two population groups and similar levels of gene diversity. According to the PCoA analysis, individu- als from Argoni and Nati form two largely overlapping clusters. The observed low inter-group differentiation reflects ongoing gene flow hence dispersal movements between the two localities may occur. In fact, these two localities are small brooks located in the opposite sides of the central mountain ridge and are separated by the main island road. The smallest distance between them is few hundred meters (approx. 500 m), a distance that lies with- in the dispersal capability of the species, thus dispersal movements along an elevational gradient could take place leading in a higher connectivity between them. Further analysis with other genetic markers such as microsatellites could offer an insight using more polymorphic loci, and assess recent events of gene flow (Manel et al., 2005). A primary objective of conservation actions is to maintain levels of diversity and heterozygosity, which are strongly recommended due to the observed popula- tion decline on Karpathos. Currently, the largest threat to populations of P. cerigensis is the loss of suitable breeding sites due to fragmentation and climatic change (Beerli et al., 2009). Small wetlands are unique ecosystems preserv- ing amphibian populations and providing patches of suit- able habitats and the overall pressure of climatic change is also apparent in other species inhabiting such fragile eco- systems. Hence, there should be increased efforts in the conservation of these areas that provide suitable habitats for the Karpathos marsh frog along with future studies in the biology and behavior of the species. ACKNOWLEDGEMENTS Sampling took place according to the Hellenic National Law (Presidential Decree 67/81) and under a special permit (code: 6ΨT04653Π8-ΥΟ3) issued by the Ministry of the Environment. The study was fully funded by the Mohamed bin Zayed Species Conservation Fund. SUPPLEMENTARY MATERIAL Supplementary material associated with this article can be found at manuscript number 23189. 122 Elisavet A. Toli et alii REFERENCES Akın, Ç., Can Bilgin, C., Beerli, P., Westaway, R., Ohst, T., Litvinchuk, S.N., Uzzel, T., Bilgin, M., Hotz, H., Guex, G.D., Plotner, J. (2010): Phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs were determined by geological processes and climate change in the Late Cenozoic. J. Biogeogr. 37: 2111- 2124. Bandelt, H., Forster, P., Röhl, A. (1999): Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16: 37-48. Beerli, P., Hotz, H., Tunner, H.G., Heppich, S., Uzzell, T. (1994): Two new water frog species from the Aegean islands Crete and Karpathos (Amphibia, Salientia, Ranidae). Notulae Naturae 470: 1-9. Beerli, P., Hotz, H., Uzzell, T. (1996): Geologically dated sea barriers calibrate a protein clock for Aegean water frogs. Evolution 50: 1676-1687. Beerli, P., Uzzell, T., Lymberakis, P. (2009): Pelophylax cerigensis. The IUCN Red List of Threatened Species 2009: e. T58567A11787309. http://dx.doi.org/10.2305/ IUCN.UK.2009.RLTS.T58567A11787309.en [accessed on 29 January 2018]. Boger, H., Dermitzakis, M.D. (1985): Neogene paleogeog- raphy in the Central Aegean region. In: VIIIth Inter- national Neogene Congress, vol. 70, Budapest. Canestrelli, D., Aloise, G., Cecchetti, S., Nascetti, G. (2010): Birth of a hotspot of intraspecific genetic diversity: notes from the underground. Mol. Ecol. 19: 5432-5451. Douris, V., Cameron, R.A., Rodakis, G.C., Lecanidou, R. (1998): Mitochondrial phylogeography of the land snail Albinaria in Crete: long‐term geological and short‐term vicariance effects. Evolution 52: 116-125. Evanno, G., Regnaut, S., Goudet, J. (2005): Detecting the number of clusters of individuals using the soft- ware STRUCTURE: a simulation study. Mol. Ecol. 14: 2611-2620. Kocher, T.D., Thomas, W.K., Meyer, A., Edwards, S.V., Pääbo, S., Villablanca, F.X., Wilson, A. C. (1989): Dynamics of mitochondrial DNA evolution in ani- mals: amplification and sequencing with conserved primers. P. Natl. Acad. Sci. 86: 6196-6200. Kumar, S., Stecher, G., Tamura, K. (2016): MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33: 1870-1874. Le, S., Josse, J., Husson, F. (2008): FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw. 25: 1-18.10.18637/jss.v025.i01. Leigh, J.W, Bryant, D. (2015): PopART: Full-feature soft- ware for haplotype network construction. Methods Ecol. Evol. 6: 1110-1116. Librado, P., Rozas, J. (2009): DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451-1452. Litvinchuk, S., Skorinov, D., Rosanov, J. (2015): Natural spontaneous autotriploidy in the genus Pelophylax (Anura: Ranidae). Russ. J. Herpetol. 22: 318-320. Lymberakis, P., Poulakakis, N., Manthalou, G., Tsigeno- poulos, C., Magoulas, A., Mylonas, M. (2007): Mito- chondrial phylogeography of Rana (Pelophylax) pop- ulations in the Eastern Mediterranean region. Mol. Phylogenet. Evol 44: 115-125. Lynch, M., Milligan, B.G. (1994): Analysis of population genetic structure with RAPD markers. Mol. Ecol. 3: 91-99. Manel, S., Gaggiotti, O.E., Waples, R.S. (2005): Assign- ment methods: matching biological questions with appropriate techniques. Trends Ecol. Evol 20: 136-142. Meudt, H.M., Clarke, A.C. (2007): Almost forgotten or latest practice? AFLP applications, analyses and advances. Trends Plant Sci. 12: 106-117. Michaux, J., Libois, R., Paradis, E., Filippucci, M.G. (2004): Phylogeographic history of the yellow-necked fieldmouse (Apodemus flavicollis) in Europe and in the Near and Middle East. Mol. Phylogenet. Evol. 32: 788-798. Nei, M. (1973): Analysis of gene diversity in subdivided populations. P. Natl. Acad. Sci. USA 70: 3321-3323. Padial, J.M., Köhler, J., Munoz, A., De la Riva, I. (2008): Assessing the taxonomic status of tropical frogs through bioacoustics: geographical variation in the advertisement calls in the Eleutherodactylus discoidalis species group (Anura). Zool. J. Linn. Soc. 152: 353-365. Pafilis, P., Kapsalas, G., Lymberakis, P., Protopappas, D., Sotiropoulos, K. (2018): Diet composition of the Kar- pathos marsh frog (Pelophylax cerigensis): what does the most endangered frog in Europe eat? Anim. Bio- div. Conserv. 42. doi: 10.32800/abc.2018.42.0001. Plötner, J., Baier, F., Akın, C., Mazepa, G., Schreiber, R., Beerli, P., Litvinchuk, S. N., Bilgin, C.C., Borkin, L., Uzzell, T. (2012): Genetic data reveal that water frogs of Cyprus (genus Pelophylax) are an endemic species of Messinian origin. Zoosyst. Evol. 88: 261-283. Poulakakis, N., Lymberakis, P., Antoniou, A., Chalkia, D., Zouros, E., Mylonas, M., Valakos, E. (2003): Molecu- lar phylogeny and biogeography of the wall-lizard Podarcis erhardii (Squamata: Lacertidae). Mol. Phylo- genet. Evol. 28: 38-46. Pritchard, J.K., Stephens, M., Donnelly, P. (2000): Infer- ence of population structure using multilocus geno- type data. Genetics 155: 945-959. R Core Team (2017): R: A language and environment for statistical computing. R Foundation for Statistical 123Phylogeny and population structure of the Karpathos frog Computing, Vienna, Austria. URL https://www.Rpro- ject.org/. Saitou, N., Nei, M. (1987): The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4: 406-425. Schneider, H., Sinsch, U. (1999): Taxonomic reassessment of Middle Eastern water frogs: Bioacoustic variation among populations considered as Rana ridibunda , R. bedriagae or R. levantina. J. Zool. Syst. Evol. Research 37: 57-65. Sfenthourakis, S., Legakis, A. (2001): Hotspots of endem- ic terrestrial invertebrates in southern Greece. Biodiv- ers. Conserv. 10: 1387-1417. Sotiropoulos, K., Eleftherakos, K., Džukić, G., Kalezić, M., Legakis, A., Polymeni, R. (2007): Phylogeny and biogeography of the alpine newt Mesotriton alpestris (Salamandridae, Caudata), inferred from mtDNA sequences. Mol. Phylogenet. Evol. 45: 211-226. Temple, H.J., Cox, N.A. (2009): European Red List of Amphibians. Office for Official Publications of the European Communities, Luxembourg. Thompson, J.D., Higgins, D.G., Gibson, T.J. (1994): CLUSTAL W: improving the sensitivity of progres- sive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22: 4673-4680. Valakos, E.D., Pafilis, P., Sotiropoulos, K., Lymperakis, P., Maragou, P. (2008): The amphibians and reptiles of Greece. Edition Chimaira, Frankfurt am Main. Vekemans, X. (2002): AFLP-SURV 1.0: A program for genetic diversity analysis with AFLP (and RAPD) population data. Laboratoire de Génétique et d’Ecologie Végétales, Université Libre de Bruxelles. Vervust, B., Brecko, J., Pafilis, P. (2013): Preliminary data on the differentiation of Milos Island water frogs. J. Biol. Res.-Thessalon. 20: 290-295. Vos, P., Hogers, R., Bleeker, M., Reijans, M., van de Lee,T., Hornes, M., Friters, A., Pot, J., Paleman, J., Kui- per, M., Zabeau, M. (1995): AFLP a new technique for DNA fingerprinting. Nucleic Acids Res. 23: 319-32. Wright, S. (1950): Genetical structure of populations. Nature 166: 247-249. Acta Herpetologica Vol. 13, n. 1 - June 2018 Firenze University Press