Abstract
In the last decade, taxonomic studies have drastically increased the number of species known to inhabit the Arabian deserts. While ongoing phylogenetic studies continue to identify new species and high levels of intraspecific genetic diversity, few studies have yet explored the biogeographic patterns in this arid region using an integrative approach.
In the present work, we apply different phylogenetic methods to infer relationships within the Palearctic naked-toed geckos. We specifically address for the first time the taxonomy and biogeography of Bunopus spatalurus Anderson, 1901, from Arabia using multilocus concatenated and species tree phylogenies, haplotype networks and morphology.
We also use species distribution modelling and phylogeographic interpolation to explore the phylogeographic structure of Bunopus spatalurus hajarensis in the Hajar Mountains and the roles of climatic stability and possible biogeographic barriers on lineage occurrence and contact zones in this arid mountain endemism hot spot. According to the inferred topology recovered using concatenated and species tree methods, the genus 'Bunopus' is polyphyletic.
Bunopus tuberculatus and B. blanfordii form a highly supported clade closely related to Crossobamon orientalis, while the two subspecies of 'Bunopus' spatalurus branch together as an independent highly supported clade that diverged during the Miocene according to our estimations. Within B. s. hajarensis, three geographically structured clades can be recognized that according to our estimations diverged during the Late Miocene to Pliocene.