Phylogenomics reveal species limits and inter-relationships in the narrow-range endemic lycian salamanders

Salamanders of the genus Lyciasalamandra are represented by as many as 20 narrow-range endemic taxa inhabiting the Mediterranean coast of Turkey and a handful of Aegean Islands. Despite recent molecular phylogenetic studies, the genus is rife with uncertainty about the number of contained species and their phylogenetic relationships, both of which can interfere with needed conservation actions. To test species limits and infer interrelationships we generated as many as 113,176 RAD loci containing 229,427 single nucleotide polymorphisms (SNPs), for 110 specimens of Lyciasalamandra representing 19 of the 20 described taxa. Through a conservative species delimitation approach, we found support for eight species in the genus which broadly agree with currently described species-level diversity. We then use multiple coalescent-based species tree methods to resolve relationships in this relatively old, synchronous species radiation. We recommend synonymization of the largely over-split subspecific taxa, and the elevation of L. luschani finikensis to full species status as L. finikensis. Our hope is that this revised taxonomic framework provides a stable foundation for conservation management in these fragile, microendemic taxa.