Molecular phylogeny and systematic of the Schneider's skink Eumeces schneiderii (Daudin, 1802) (Squamata: Scincidae)

Eumeces is a little‐known group of skinks with few phylogenies published. Within its distribution range, five subspecies were identified for Schneider's Skink (Eumeces schneiderii) based on color pattern, pholidosis, intra‐ and inter‐subspecific relationships, but the validity of these taxa has yet to be assessed using molecular markers. Here, a robust phylogeny of the E. schneiderii group is presented based on three molecular markers (Cytb, 16S rRNA, and c‐mos) from 80 samples collected across Anatolia and the Iranian Plateau. Both Maximum Likelihood and Bayesian Inference approaches were used to infer phylogenetic relationships within this group. The results revealed monophyly of E. schneiderii subspecies and Eumeces persicus. These clades also differ in habitat preferences, with the small‐bodied, striped forms generally inhabiting lower elevations and drier environments than the larger, uniform morph E. schneiderii subspecies. Molecular‐clock dating based on Cytb implies that most speciation events took place during the Upper Miocene to Lower Pliocene. A strongly supported basal dichotomy of Eumeces schneiderii zarudnyi + E. persicus in the southern and eastern slopes of the Zagros Mountains, with the E. schneiderii subspecies in the western Zagros, was identified in all individual and concatenated trees. Eumeces schneiderii zarudnyi should be elevated to species level and, therefore, we propose the name Eumeces zarudnyi for this taxon, increasing the number of species in the genus to seven described species. We did sampling across the distribution range of the Eumeces schneiderii complex (2016–2019) and sampled 70 complete individuals of five currently known subspecies. Our dataset comprised 80 sequences of 1560 bp consisting of 16S, Cytb, c‐mos. We performed probabilistic approaches including Maximum Likelihood and Bayesian Inference for the dataset and acquired two major clades (Clade A and Clade B [B1 + B2]). The center of the divergence determined in central Iran for the clade A and Eastern clade B1 and NE Mediterranean Basin as for the Western clade (B2).