Distributed under Creative Commons Cc-by 4.0 an Integrative Systematic Revision and Biogeography of Rhynchocalamus Snakes (reptilia, Colubridae) with a Description of a New Species from Israel

Background. The colubrid snakes of the genus Rhynchocalamus are seldom studied and knowledge of their ecology and life history is scarce. Three species of Rhynchocalamus are currently recognized, R. satunini (from Turkey eastwards to Iran), R. arabicus (Yemen and Oman), and R. melanocephalus (from the Sinai Peninsula northwards to Turkey). All are slender, secretive, mainly nocturnal and rare fossorial snakes. This comprehensive study is the first to sample all known Rhynchocalamus species in order to review the intra-generic phylogenetic relationships and historical biogeography of the genus. Methods. We revised the systematics of Rhynchocalamus using an integrative approach and evaluated its phylogeography. The phylogenetic position within the Colubridae and the phylogenetic relationships within the genus were inferred using 29 individuals belonging to the three known species, with additional sampling of two other closely-related genera, Muhtarophis and Lytorhynchus. We analysed three mitochondrial (12S, 16S, cytb) and one nuclear (c-mos) gene fragments. Phylogenetic trees were reconstructed using maximum likelihood and Bayesian inference methods; the latter method also used to provide the first time-calibrated molecular phylogeny of the genus. We generated a nuclear network and carried out a topology test and species delimitation analysis. Morphological comparisons were used to differentiate among species and to describe a new species from Israel. The studied material was comprised of 108 alcohol-preserved specimens, 15 photographs, and data from the literature for the examination of 17 mensural, 14 meristic, and two categorical characters. Results. The molecular results support Rhynchocalamus as monophyletic, and as having split from its sister genus Lytorhynchus during the Late Oligocene. The three recognized species of Rhynchocalamus comprise four independently evolving groups. The molecular results reveal that the genus began to diverge during the Middle Miocene. We revealed that the best-studied species, R. melanocephalus, is paraphyletic. A population, formally ascribed to this species, from the Negev Mountain area in southern Israel is phylogenetically closer to R. arabicus from Oman than to the northern


INTRODUCTION
Taxonomy today often relies on molecular data for further support and information. Such data are usually preferred over morphology for the reconstruction of evolutionary relationships among organisms (Hebert et al., 2003;Tautz et al., 2003;Blaxter, 2004;Vogler & Monaghan, 2007;Padial et al., 2010). The increasing use and availability of molecular data has led to the development of new methods to study systematics (Sites & Marshall, 2003), and has proven to be an invaluable tool for evaluating the evolutionary relationships between both closely and distantly related species. Recent studies of Middle Eastern snakes have used molecular data to elucidate the inter-and intra-specific relationships among taxa, revealing high levels of genetic differentiation and cryptic diversity that do not accord with the current taxonomy. Such studies have also provided insights into the historical biogeography of the taxa and the processes that triggered their diversification (e.g., Lenk et al., 2001;Utiger et al., 2002;Nagy et al., 2003;Nagy et al., 2004;Pook et al., 2009;Stümpel & Joger, 2009;Kornilios et al., 2013). However, the biodiversity of snakes in the Middle East remains unclear, as systematic and biogeographic data for several genera are still lacking. One such example is that of the colubrid genus Rhynchocalamus Günther, 1864. Rhynchocalamus snakes are secretive, non-venomous and occasionally found near human habitations. They are poorly known and information regarding their natural history is scarce. These are small-sized, slender, fossorial aglyphous snakes that are mostly nocturnal but can also be found active during the day (Gasperetti, 1988;Disi et al., 2001;Baha El Din, 1994;Baha El Din, 2006;Avci et al., 2007;Avci et al., 2008). Morphologically they are characterized by a thin cylindrical body and short tail, small head indistinct from the neck, an enlarged rostral shield wedged between the internasals scales, and divided anal plate and subcaudal scales (Gasperetti, 1988;Disi et al., 2001;Baha El Din, 2006). Rhynchocalamus snakes prefer humid areas with little vegetation and are found in mountainous areas in both Mediterranean and Irano-Turanian ecozones (Fig. 1), on heavy soils but not on sand. In southern Israel, Jordan, and Egypt (i.e., the Sinai Peninsula), they Localities of the material analysed in this study, including type localities (star), samples used for the genetic analyses (circle), specimen vouchers for the morphological examinations (diamond) and photographic material (square). Numbers correspond to samples listed in Table S1 and colours to specimens in Fig. 2 and Figs. S1-S4. Taxon names correspond to changes proposed in this paper. Spatial data modified from various sources (IUCN-http://www.iucnredlist.org; Bar & Haimovitch, 2013;Sindaco, Venchi & Grieco, 2013). are also known from arid and stony steppes, sparsely vegetated rocky slopes and wadis (Disi et al., 2001;Baha El Din, 2006;Amr & Disi, 2011;Bar & Haimovitch, 2013;Werner, 2016).

.continued)
Age estimates based on external calibration points of Hemorrhois and Hierophis subgroups (see 'Materials and Methods') are indicated near the relevant nodes and include the mean and, between parentheses, the 95% highest posterior densities (HPD) confidence interval. (B) Haplotype network of the c-mos nuclear marker. Circle size is proportional to the number of alleles. Sample codes and colours correlate to specimens in Table S1 and in Fig. 1 and Figs. S1-S4. Taxon names correspond to changes proposed in this paper.
Phylogenetic data on this reclusive genus were, until recently, based on morphology alone-it was classified as an Oriental genus with most of the known species occurring in southern Asia, and at times, it was assigned to the South-east Asian genus Oligodon (Boulenger, 1894;Werner, 1905;Chernov, 1937;Bodenheimer, 1944;Haas, 1952;Darevsky, 1970;Amr & Disi, 2011). As no DNA sequences of this genus were available, the squamate phylogeny by Pyron, Burbrink & Wiens (2013) did not include Rhynchocalamus, leaving its phylogenetic position unverified. Recently, two studies provided genetic data on the genus (Avci et al., 2015;Šmíd et al., 2015) helping to confirm the phylogenetic affinities. Rhynchocalamus was found to be closely related to the genus Lytorhynchus within the Western Palearctic colubrid clade (Avci et al., 2015;Šmíd et al., 2015), contradicting the previous hypothesis of an Oriental origin (this relationship, however, was not resolved in the latest squamate phylogeny by Tonini et al., 2016). Although both Avci et al. (2015) and Šmíd et al. (2015) included only four specimens each, the incorporation of molecular data revealed interesting results, such as the first sampling and phylogenetic position of R. arabicus (Šmíd et al., 2015) and the elevation of R. satunini to species level (Avci et al., 2015). Additional taxonomic changes based on the phylogenetic analyses resulted in the classification of Rhynchocalamus barani Olgun et al., 2007, a Turkish endemic, within the new genus Muhtarophis Avci et al., 2015, distantly related to Rhynchocalamus. The current status of the recognized species within Rhynchocalamus, their relationships and distribution, remain relatively unclear, as no study has sampled all known species from the entire distribution range of the genus. In this work, we explore the phylogenetic relationships within Rhynchocalamus by means of a broad sampling coupled with a morphological revision. Using an integrative taxonomic approach (Dayrat, 2005), we seek to produce the most complete phylogeny of Rhynchocalamus to date, in order to clarify its systematics, describe a new species from Israel, and elucidate its biogeographical and evolutionary history.

Taxon sampling, DNA extraction and amplification
In order to resolve the phylogenetic relationships within Rhynchocalamus, 29 individuals belonging to all three recognized species were used in the molecular study, including eight sequences from Avci et al. (2015) and Šmíd et al. (2015) that were retrieved from GenBank. In order to evaluate the phylogenetic position of Rhynchocalamus within the Western Palearctic colubrid clade, we used 40 sequences of different members retrieved from GenBank. Previous studies have identified Rhynchocalamus to be closely related to Lytorhynchus (Šmíd et al., 2015) and the newly-described genus Muhtarophis (Avci et al., 2015). We therefore additionally sequenced seven Lytorhynchus and six Muhtarophis individuals. A list of all Rhynchocalamus individuals and the members of the Western Palearctic colubrid clade included in the molecular analyses, with their localities and GenBank accession codes, is presented in Table S1. Localities of Rhynchocalamus individuals are shown in Fig. 1.
Genomic DNA was extracted from ethanol-preserved tissue samples using the Speed-Tools Tissue DNA Extraction kit (Biotools, Madrid, Spain). Individuals were sequenced for the following markers: three mitochondrial gene fragments, ribosomal 12S rRNA and 16S rRNA (12S and 16S, respectively), cytochrome b (cytb), and the nuclear gene fragment oocyte maturation factor MOS (c-mos). Gene fragments were amplified and sequenced for both strands using published primers (as described in detail by Šmíd et al. (2015); when the cytb long fragment failed to amplify, we used shorter sequences with the Gludg and Cytb2 primers; Kocher et al., 1989;Palumbi, 1996, respectively).

Sequence analysis, phylogenetic analyses and hypothesis testing
Chromatographs were checked manually, assembled and edited using Geneious v.7.1.9 (Biomatter Ltd.). For the c-mos gene fragment, heterozygous positions were identified and coded according to the IUPAC ambiguity codes in both alleles. Coding gene fragments (cytb, c-mos) were translated into amino acids and no stop codons were observed, suggesting that the sequences were all functional and were trimmed to start at the first codon position. DNA sequences were aligned for each gene independently using the online version of MAFFT v.7 (Katoh & Standley, 2013) with default parameters (Auto strategy, Gap opening penalty: 1.53, Offset value: 0.0). For the 12S and 16S ribosomal fragments we applied the Q-INS-i strategy, in which information on the secondary structure of the RNA is considered. To remove poorly aligned positions of the non-protein-coding 12S and 16S we used Gblocks (Castresana, 2000) with low stringency options (Talavera & Castresana, 2007). Inter and intra-specific uncorrected p-distances and the number of variable (V ) and parsimony informative (Pi) sites were calculated in MEGA v.7 (Kumar, Stecher & Tamura, 2016) independently for each gene fragment.
We analysed the Rhynchocalamus data using four datasets assembled using different species-sets (A, B and C) and different DNA sequences: (i) Species-set A was assembled with the aim of resolving the phylogenetic position of Rhynchocalamus within Colubrinae and to obtain dates for some relevant cladogenetic events. Dataset 1 of concatenated mtDNA and nDNA, comprised 82 specimens corresponding to 46 different taxa. Dataset 2, for calibration only, consisted of the concatenated mtDNA and nDNA with a single representative of each mPTP clade of Rhynchocalamus (see below for information on the mPTP species delimitation analysis). It comprised 48 specimens corresponding to 46 different taxa; (ii) species-set B was assembled with the aim of resolving the phylogenetic relationships within Rhynchocalamus. Dataset 3 of concatenated mtDNA and nDNA included 36 specimens corresponding to six different taxa (Lytorhynchus was used to root the tree); (iii) species-set C was assembled with the aim of evaluating the relationships and species boundaries within Rhynchocalamus. Dataset 4 of mtDNA haplotypes only included 28 specimens.
Best-fit partitioning schemes and models of molecular evolution were selected with PartitionFinder v.1.1.1 (Lanfear et al., 2012) using the following parameters: branchlengths (linked); models of evolution (beast); model selection (AIC); data blocks (12S and 16S each as a single partition, cytb and c-mos first and second codon positions of each gene as one partition, and the third codon positions of each gene as another); search scheme (all). A summary of the analyses performed with each dataset, including the different partitions and models is shown in Table S2.
Phylogenetic analyses were performed using maximum likelihood (ML) and Bayesian inference (BI) methods. In both ML and BI analyses, alignment gaps were treated as missing data and the nuclear gene sequences were not phased. Partitions and models for each dataset, with priors, specifications and parameters for each analysis are specified in Table S2. Maximum likelihood trees were estimated in RAxML v.7.4.2 (Stamatakis, 2006) as implemented in raxmlGUI v.1.3 (Silvestro & Michalak, 2012). All ML analyses were performed with a GTR + G model of sequence evolution and 100 replicates. Each inference was initiated with a random starting tree and nodal support was assessed with 1,000 bootstrap pseudoreplicates (Felsenstein, 1985). Bayesian analyses were performed with BEAST v.1.8.2 (Drummond et al., 2012). All BEAST analyses were carried out in CIPRES Science Gateway (Miller, Pfeiffer & Schwartz, 2010), and the .xml file was manually modified to ''Ambiguities = true'' for the nuclear partition (c-mos) to account for variability in the heterozygote positions, rather than treating them as missing data. For all analyses implemented in BEAST, the convergence of runs was assessed by the effective sample size (ESS) values of parameters (>200) using TRACER v.1. 6 (Rambaut et al., 2014). LogCombiner and TreeAnnotator (both available in BEAST package) were used to infer the ultrametric tree after discarding 10% as burn-in. Nodes were considered strongly supported if they received ML bootstrap values ≥ 70% and posterior probability (pp) support values ≥ 0.95 (Wilcox et al., 2002;Huelsenbeck & Rannala, 2004).
To identify divergent lineages within Rhynchocalamus, putative species boundaries were tested using the multi-rate Poisson Tree Processes (mPTP; Kapli et al., 2016) model, using the webserver (http://mptp.h-its.org/). This is an improved method of the previously published species delimitation method PTP (Zhang et al., 2013). As this analysis relies on single locus data, we reconstructed a ML haplotype concatenated mitochondrial phylogenetic tree as specified above for dataset 4.
A nuclear network was constructed for the nuclear gene c-mos. SeqPHASE (Flot, 2010) was used to convert the input files and the software PHASE v.2.1.1 to resolve phased haplotypes (Stephens, Smith & Donnelly, 2001;Stephens & Scheet, 2005). Default settings of PHASE were used, except for phase probabilities, which were set as 0.6. The phased nuclear sequences were used to generate a median-joining network using NETWORKS v.5 (Bandelt, Forster & Röhl, 1999).

Divergence time estimates
Unfortunately, no calibration data for Rhynchocalamus are currently known, precluding the use of internal calibration points and preventing a direct estimation of the time in our phylogeny. Therefore, for a temporal framework we used calibrations of other members of the Western Palearctic colubrid clade. We used two calibration points previously used by several authors (Nagy et al., 2003;Nagy et al., 2004;Wüster et al., 2008;Pook et al., 2009): (i) Hemorrhois divergence between the eastern (H. ravergieri and H. nummifer) and western (H. algirus and H. hippocrepis) subgroups occurred after the contact of Africa-Arabia with Eurasia, 16-18 million years ago (Mya;Nagy et al., 2003); we applied a Normal distribution, mean 18, stdev 2 (95% confidence interval of 14.7-21.3 Mya); (ii) Hierophis subgroup divergence, including Eirenis, 18 Mya according to fossil data (Ivanov, 2002); we applied a Normal distribution, mean 18, stdev 1 (95% confidence interval of 16.4-19.6 Mya). The dataset for the estimation of divergence times in BEAST v.1.8.2 was dataset 2 and comprised the members of the Western Palearctic colubrid clade (species-set A; Table S1) and one representative of each independent mPTP entity of Rhynchocalamus (the nuclear genes unphased; Table S1). Partitions, models, priors and parameters are specified in Table S2.

Morphological material and museum acronyms
The material studied was obtained from the following institutions: [ Material for the morphological revision of Rhynchocalamus comprised 108 alcoholpreserved specimens. In addition, 15 photographed voucher specimens were taxonomically identified and examined (used for meristic, categorical and colouration characters, but not measured), including the lectotype (MHNG1246.77) and previous holotype (BMNH1946.1.3.29) of R. melanocephalus; the holotype of R. satunini (ZISP9343) and the holotype of R. arabicus (FMNH18219, including the pictures published by Šmíd et al., 2015). Localities of the voucher specimens and of the photographic material are presented in Fig. 1. Data from the literature were used to complete missing morphological traits from photographed voucher specimens and to enlarge the sample size for the mensural and meristic comparisons (21 specimens;Nikolsky, 1899;Schmidt, 1933;Reed & Marx, 1959;Darevsky, 1970;Franzen & Bischoff, 1995;Avci et al., 2007;Avci et al., 2008;Šmíd et al., 2015). A list of all examined specimens including their museum accession codes and localities is presented in Table S3.

Morphological characters
Characters for the morphological analyses were selected based on previous taxonomic studies of the genus Rhynchocalamus (Darevsky, 1970;Gasperetti, 1988;Franzen & Bischoff, 1995;Olgun et al., 2007;Avci et al., 2008) and on personal observations. The following mensural characters were taken by the first author (KT; two specimens from the ZMHRU collection were examined by BG) on the right side of each specimen (if bilateral) using Helios callipers with an accuracy to the nearest 0.01 mm and, where necessary, under a stereomicroscope: snout-vent length (SVL), measured from tip of snout to vent; tail length (TL), from vent to tip of tail; pileus length (pilL), measured from tip of snout to posterior margin of parietals; Parietal length (PL); Parietal width (PW); frontal length (FL); frontal width (FW); rostral length (RL); rostral height (RH); rostral width (RW); prefrontal suture length (PFL); internasal suture length (IntNL); eye diameter (EYE); black pattern length (BPL; in R. melanocephalus only), from the head to last black dorsal; distance between nostrils (InD); anterior inframaxillars length (AimL); posterior inframaxillars length (PimL).
In addition to the morphometric continuous variables, the following meristic (pholidotic) variables were collected by the same person (KT) using a dissecting microscope: number of preoculars (PreO); number of postoculars (PostO); number of temporal scales (TS); number of post-temporal scales (PTS); number of loreal scales (LS); number of ventrals (VS); number of subcaudal scales (SCS); number of upper labial scales (UL); number of lower labial scales (LL); number of black dorsal scales (BDS; number of black scales at mid-body from between the parietals to the end of colour pattern; in R. melanocephalus only); number of lower labials in contact with anterior inframaxillar (InfLC); number of gular scales in a row between posterior inframaxillars (GSI); number of gular scales in a row between posterior inframaxillars and 1st ventral (GS); number of dorsal and temporal scales surrounding the margin of parietals (DST).
Categorical characters describing the degree of size and presence of different scales (not modified in preserved specimens) were: size of the 3rd and 4th upper labial scales (34UL; equal/large) and shape of internasal scales (IntN; triangle/trapezoid).

Statistical analyses
Statistical analyses were used to investigate if there are differences in the mensural or meristic characters between R. melanocephalus and the new species described herein. The 17 mensural and 14 meristic characters were analysed independently, and the two categorical characters were directly used in the description of the new species (see systematic account below).
All the mensural variables were log 10 -transformed before the analyses, and the different datasets were tested for normality using the Shaphiro-Wilk's test and homogeneity of variances using the Leven's test; if normality or homogeneity were not present, we used a permutation test. All the statistical analyses were performed in IBM SPSS Statistics 23 (IBM Corp. Armonk, NY, USA). Morphological differences were tested for the meristic variables which presented differences within the assemblage using a one-way ANOVA (for VS, SCS, BDS, DST) and Fisher's exact probability test (for PostO, LL, InfLC). The mensural traits were tested using one-way ANCOVA (SVL as a covariate for size correction; adult specimens only). As a result of the presence of a single female in the new species, sexual size comparisons could only be tested for R. melanocephalus which has an adequate sample size of both males and females, and was tested for each variable as described above.

Species distribution models
We analysed species distribution models (SDM) using Maxent v.3.3 (Phillips, Anderson & Schapire, 2006) to assess which environmental variables shape Rhynchocalamus distribution and whether the species' ranges could potentially overlap. As a presence-only model, Maxent does not require absence data, which are nearly impossible to obtain for these secretive snakes, and it performs well even with datasets of small sample sizes (Phillips & Dudík, 2008), as is the case with the new species from Israel. Specimen localities were used as the input data. In total, we used 83 records of R. melanocephalus, 15 of R. satunini, and six of the new species; R. arabicus with its two known localities was not included. All three models were developed on the same spatial extent that is shown in Fig. 1. The spatial extent was defined as a 250 km buffer around a polygon encompassing all the localities. The models were based on 19 present-day bioclimatic variables (WorldClim database v.1.4; Hijmans et al., 2005), global land cover data (European Space Agency; http://due.esrin.esa.int/page_globcover.php), and global soil type (FAO; http://www.fao.org/), all at a resolution of 30 arc-seconds. We used ENMTools (Warren, Glor & Turelli, 2010) to calculate Pearson's correlation coefficient to measure correlation between the climatic variables and retained only variables correlated less than 0.75: altitude, mean diurnal temperature range (BIO2), isothermality (BIO3), temperature seasonality (BIO4), mean temperature of wettest quarter (BIO8), mean temperature of driest quarter (BIO9), mean temperature of coldest quarter (BIO11), precipitation seasonality (BIO15), precipitation of wettest quarter (BIO16), precipitation of driest quarter (BIO17). In Maxent, we used a maximum of 5,000 iterations, 10 replicate runs and 25% of the data was used as training samples. Other settings were left at default. We reclassified the continuous models into binary presence-absence maps using the maximum training sensitivity plus specificity threshold (MTSS). The area under the receiver operating characteristics curve (AUC) was taken as a measure of overall model accuracy.

Zoobank registration and collection of specimens
The electronic version of this article in Portable Document Format (PDF) will represent a published work according to the International Commission on Zoological Nomenclature (ICZN), and hence the new names contained in the electronic version are effectively published under that Code from the electronic edition alone. This published work and the nomenclatural acts it contains have been registered in ZooBank, the online registration system for the ICZN. The ZooBank LSIDs (Life Science Identifiers) can be resolved and the associated information viewed through any standard web browser by appending the LSID to the prefix http://zoobank.org/. The LSIDs for this publication is: urn:lsid:zoobank.org:pub:1E72A585-C0C3-4835-9466-25E8A7C9FADB. The online version of this work is archived and available from the following digital repositories: PubMed Central and CLOCKSS.
The authors received an ethical permission (Ege University Animal Experiments Ethics Committee, 2010#43) and special permission (2014#51946) for field studies from the Republic of Turkey, Ministry of Forestry and Water Affairs. Genetic and morphological analyses of specimens from other localities were based on museum collections and no specimens were collected for this work. All efforts were made to minimize animal suffering.

Molecular analyses
The phylogenetic results within the Western Palearctic colubrid clade (dataset 1; Fig. S1), for 82 sequences of 46 taxa, agree with the phylogenetic placement of the new genus Muhtarophis as distinct from Rhynchocalamus, although its phylogenetic position within the clade remains unresolved (no support in the ML analysis, but strong support in the BI analysis). Wallaceophis is sister to a clade comprising Rhynchocalamus and Lytorhynchus, though with no support in either ML or BI analyses. Both Rhynchocalamus and Lytorhynchus are recovered as monophyletic sister genera with strong support.
Rhynchocalamus ( Fig. 2A) is divided into four separate lineages. The monophyly of the distinct lineage of R. satunini from Turkey and Iran is strongly supported. Within this species, the three specimens from Turkey form an inner group, while the specimen from Iran is relatively distant. This species is more genetically diverse than the other species (12S: 1%; 16S: 0.4%; cytb: 1.8%; Table S1). The sole specimen of R. arabicus from Oman is genetically distinct. The results reveal that R. melanocephalus is paraphyletic, as the southern Israeli lineage from the Negev Mountain in Israel (the new species described herein, see Systematic account section below) is phylogenetically more closely related to R. arabicus from Oman than to the geographically adjacent populations of R. melanocephalus from the Negev region in Israel and northwards to Turkey (Fig. 1). In addition, within the R. melanocephalus lineage, one sample from Mt. Hermon in Israel (HUJ.R20967) is phylogenetically distinct from the others.
The c-mos nuclear network (Fig. 2B) reveals a pattern of incomplete lineage sorting. Allele sharing is present between the newly-described species from the Negev Mountain in Israel and the Negev population of R. melanocephalus. This is based, however, on only two samples from the Negev Mountain population (the other two samples from this population failed to amplify). Allele sharing also occurs between R. arabicus from Oman and R. satunini from Iran and Turkey.
The mPTP species delimitation approach recognized six distinct entities within Rhynchocalamus (Fig. S2): two discrete entities within R. satunini (the specimen from Iran as distinct from the three specimens from Turkey); a distinct R. arabicus; and three entities within R. melanocephalus-the southern lineage from the Negev Mountain in Israel is distinct (the new species described below), recovered as sister species to R. arabicus, separated from the rest of the R. melanocephalus specimens, and a sister relationship between the sole sample from Mt. Hermon and the rest of the samples from the Negev area in Israel northwards to Turkey.
We performed a topology test in order to better understand the relationships recovered in our analyses, by forcing the monophyly of R. melanocephalus and the new species described herein. The results of this test (AU: p < 0.0001; SH: p = 0) reject the hypothesis that R. melanocephalus and the new species described herein form a clade.
Divergence time estimates for the Western Palearctic colubrid clade are presented in Fig. S3, and near the relevant nodes for Rhynchocalamus in Fig. 2A. We selected six representatives of Rhynchocalamus to be used in the dating analysis according to the mPTP species delimitation method (Figs. S2-S3; Table S1). Our results, based on dataset 2, indicate that Rhynchocalamus split from Lytorhynchus around 26.4 Mya (95% HPD: 21.6-31.9 Mya). Rhynchocalamus started diverging with the split of R. satunini during the Middle Miocene ca. 15.5 Mya (95% HPD: 12-19.4 Mya). The divergence of R. melanocephalus is estimated to have occurred around 11.3 Mya (95% HPD: 8.7-14.4 Mya), with further radiation during the Pliocene, ca. 3.5 Mya (95% HPD: 1.7-5.8 Mya). The separation between R. arabicus and the new species described herein from the Negev Mountain in Israel appears to have occurred approximately 9.5 Mya (95% HPD: 7-12.5 Mya).

Morphological analyses
The morphological database comprised 118 specimens of R. melanocephalus (six of which belong to the new species described herein, see below), 15 specimens of R. satunini and two specimens of R. arabicus. Descriptive statistics for all 33 variables included in the analysis are presented in Table 1.
The morphological characteristics of the four lineages of Rhynchocalamus, comprising 17 mensural, 14 meristic and two categorical traits (Table 1) are as follows: cylindrical body with 15 smooth dorsal scales rows; a small head indistinct from the neck; reduced maxillary dentition with 6-8 maxillary teeth, the posterior ones being long; a rostral shield enlarged, pointed backwards and wedged between the internasals; nostril situated in an undivided elongated nasal scale; small eyes with round pupils; loreal scale usually present in all species (absent in 8.5% of R. melanocephalus); 3rd and 4th upper labial scales usually larger than the 1st and 2nd (equal in size in 19% of R. melanocephalus) and are in contact with the eye; one Table 1 Descriptive statistics for all variables examined for males and females of Rhynchocalamus. Two populations of R. melanocephalus are presented (Southern population from the Negev region in Israel; Northern population from the Mediterranean ecoregion). Mean ± standard deviation and the range (min-max) are given. Mensural variables were taken from adult specimens only and are presented in millimetres. Bilateral meristic characters are presentedt with the left/right sides. Abbreviations of variables as explained in 'Materials and Methods'. Taxon names correspond to changes proposed in this paper.  preocular scale; one postocular scale in R. satunini, R. arabicus and R. melanocephalus (the latter comprises three specimens with two postocular scales on either side), two postocular scales in the new species described herein; one temporal scale, and between one or two post-temporals (mainly two in R. melanocephalus and one in R. satunini, R. arabicus and the new species); seven upper labial scales in R. satunini, six in the others; the shape of the 1st upper labial scale varies between square and trapezoid (34% and 66%, respectively), and that of the 2nd scale is mainly square; mostly seven lower labial scales in R. melanocephalus (11% with eight scales), eight in the others, including in the new species described below; mostly four lower labial scales in contact with the anterior inframaxillars in R. satunini and the new species described below, mostly three in R. melanocephalus (74%); single gular scale in contact with anterior inframaxillars, situated between the posterior inframaxillars; posterior inframaxillars are separated by 1-2 gular scales (22% and 78%, respectively); 2-8 gular rows separate the posterior inframaxillars from the 1st ventral (3-5, 97%); 5-16 black dorsal scales at mid-body in R. melanocephalus, 10-16 in the new species described below; number of dorsal and temporal scales surrounding the margin of parietals is 8-15; 164-240 ventrals; anal and subcaudal scales divided; 29-81 subcaudal scales.

Statistical analyses of morphological data
The results of the sexual dimorphism analysis of R. melanocephalus using one-way ANOVA showed that there are no significant differences in body-size between males and females (SVL; P = 0.06). The one-way ANCOVA results (with SVL as a covariate for size correction) for the remaining mensural characters showed that significant sexual size dimorphism is present for the following 13 characters, for most the males have larger sizes (P < 0.047 for all variables; Table S4): pileus length (pilL), rostral length (RL), rostral height (RH), rostral width (RW), frontal length (FL), frontal width (FW), parietal length (PL), frontal width (PW), eye diameter (EYE), black pattern length (BPL), distance between nostrils (InD), anterior inframaxillars length (AimL), and posterior inframaxillars length (PimL). The remaining three mensural characters did not show significant differences between the sexes (P > 0.17 for all three variables). The results of the one-way ANOVA and the Fisher's exact probability test of seven meristic characters (VS, SCS, BDS, DST and PostO, LL, InfLC, respectively) showed no sexual dimorphism (P > 0.07 and 0.17 for all variables, respectively; Table S4). As a result of the presence of sexual size dimorphism within R. melanocephalus, and due to the presence of only one female of the new species described herein, differences between the two species were only tested among male specimens for each variable as described above. Statistical tests did not detect significant differences in most of the mensural and meristic traits between the new species and R. melanocephalus (P > 0.068 for all characters; Table S4). For the mensural traits significant differences were found for the following four characters (Table S4): (1) Rostral height (RH; P = 0.003), indicating that the new species has a lower rostral shield than R. melanocephalus (1.15-1.66 mm vs. 1.22-2.53 mm, respectively); (2) Parietal width (PW; P = 0.035), indicating that the new species has a narrower head width than R. melanocephalus (3.73-4.71 mm vs. 3.26-6.35 mm, respectively); (3) Distance between nostrils (InD; P < 0.0001), indicating that the new species has a shorter distance than present in R. melanocephalus (1.32-2.11 mm vs. 1.59-3.24 mm, respectively); (4) Black pattern length (BPL; P = 0.001), indicating that the new species has relatively longer dorsal black colour pattern than R. melanocephalus (16.32-22.95 mm vs. 10.77-23.25 mm, respectively). The latter finding, in contrast to the other results, is not significant when the new species is compared to the Negev population of R. melanocephalus (Localities 6-12 in Fig. 1; P = 0.307; 16.32-22.95 mm vs. 14.78-23.24 mm, respectively). For the meristic variables, two characters showed significant differences between the two assemblages (Table S4): number of ventrals (VS; P = 0.015), indicating that the new species has a higher number of ventrals than R. melanocephalus (198-229 vs. 178-234, respectively); and the number of black dorsal scales (BDS; P < 0.0001), corresponding to the mensural trait of black pattern length (BPL), which was also significant, indicating that the new species has more back dorsal scales than R. melanocephalus (10-16 vs. 5-14, respectively). As with the BPL character, this finding is not significant when the new species is compared to the Negev population of R. melanocephalus (Localities 6-12 in Fig. 1; P = 0.872; 10-16 vs. 11-16, respectively). Three meristic characters (PostO, LL, InfLC) showed significant results between the two species for each variable (P < 0.006; Table S4), indicating that the new species has a significantly higher number of postocular scales (two vs. one in R. melanocephalus), higher number of lower labials (eight vs. mostly seven in R. melanocephalus), and a higher number of lower labials in contact with anterior inframaxillar (four vs. mostly three in R. melanocephalus).

Potential distribution of Rhynchocalamus
The distribution models of R. melanocephalus, R. satunini and the new species from southern Israel had either excellent or good predictive accuracy (following Araújo et al., 2005) with AUC averaged over ten replicate runs being 0.977 ± 0.015, 0.835 ± 0.071 and 0.999 ± 0.001, respectively. The potential distributions predicted for the presence of the three species show overlapping distributions (Fig. S4). The range of R. satunini stretches from Iran westwards across the southern Caucasus to southern Turkey where it overlaps with R. melanocephalus which spans from southern Turkey along the eastern Mediterranean coast to the Sinai Peninsula. In its southernmost range R. melanocephalus overlaps with the new species which is restricted mostly to southern Israel. The two latter species have a potential presence in the northern Sinai Peninsula. The main environmental variables that contributed most to the predicted distributions were BIO4 (59.2%), soil type (16.5%), BIO16 (15.2%) for R. melanocephalus; soil type (35.9%), BIO3 (29%), BIO4 (16.3%) for R. satunini; and soil type (51.8%), BIO4 (24.6%), BIO15 (20.4%) for the new species.

Rhynchocalamus melanocephalus
Holotype. HUJ.R21704, Adult male, collected from road no. 40 near Nafha Prison, Negev Mountain, Israel, 30.7317N 34.7709E WGS84, 700 m above sea level (a.s.l.) on the 21st of June 2008 by Gal Vine (Fig. 3).  (Fig. 3) were collected in the Negev Mountain area, Israel, and were included in the molecular analysis (Table S1). Etymology. The specific epithet, ''dayanae,'' is named in honour of Professor Tamar Dayan, director of the Steinhardt Museum of Natural History at Tel Aviv University and curator of the Terrestrial Vertebrate Collection. This naming of the new species constitutes a special recognition of Professor Dayan by two of her former students (KT and SM) to acknowledge her immense contribution to the conservation of Israeli fauna, and her efforts in establishing the National Natural History Museum at Tel Aviv University, and in promoting taxonomy, conservation and ecology studies in Israel.   Differential diagnosis. Rhynchocalamus dayanae sp. nov. differs from the other species of the genus in its head scalation and its head and dorsal colour pattern (Table 1;  . It differs from its phylogenetically closely-related taxon, R. arabicus from Yemen and Oman, in its lower number of ventrals (188-229 vs. 240 in R. arabicus) and subcaudals (54-62 vs. 71-81 in R. arabicus), and the presence of two postocular scales (one in R. arabicus). The overall black colouration of R. arabicus is an additional differentiation from the solely black head and first dorsal scales of R. dayanae sp. nov. It further differs by a genetic distance of 5.9%, 4.2%, and 10% in the mitochondrial 12S, 16S, and cytb genes, respectively (Table S1). Rhynchocalamus dayanae sp. nov. differs from the geographically distant species, R. satunini from Turkey eastwards to Iran, in the lower number of upper labials (6 vs. 7 in R. satunini), and the presence of two postocular scales (one in R. satunini). The head colouration of R. satunini is not uniform black, but features two black patches on a whitish background around the prefrontal scales and the parietals, and a black band around the neck that does not reach the ventrals (the rostral shield, nasals, loreals, upper and lower labials, temporal, and ventrals are whitish). It further differs by a genetic distance of 7.3%, 5.4%, and 11.9% in the mitochondrial 12S, 16S, and cytb genes, respectively (Table S1). Rhynchocalamus dayanae sp. nov. differs from R. melanocephalus (including its Negev population) by the combination of a higher number of lower labials, eight (100%) vs. seven (90%); higher number of lower labials in contact with the anterior inframaxillars, four (100%) vs. three (74%); the presence of two postocular scales (84% vs. 3%); a higher number of ventrals (198-229 vs. 178-234); lower rostral shield (1.15-1.66 mm vs. 1.22-2.53 mm); a narrower head width (3.73-4.71 mm vs. 3.26-6.35 mm); a narrower distance between nostrils (1.32-2.11 mm vs. 1.59-3.24 mm). Rhynchocalamus dayanae sp. nov. differs from the northern, Mediterranean populations of R. melanocephalus (from central Israel northwards) by the larger extent of the black dorsal pattern (BDS, 13 ± 2.1 (10-16) vs. 7.8 ± 1.8 (5-14); BPL,) mm, respectively) and black gulars (absent in the northern R. melanocephalus). Rhynchocalamus dayanae sp. nov. further differs from all populations of R. melanocephalus by a genetic distance of 6.7%, 4%, and 10.2% in the mitochondrial 12S, 16S, and cytb genes, respectively (Table S1).

Diagnosis. A new species of
Description of the holotype. Adult male (HUJ.R21704) (Fig. 3). Snout-vent length 259.8 mm, tail length 59.2 mm. Six to eight maxillary teeth, the posterior teeth are long and strong, broad at the base with an impression but without a longitudinal groove, palatine teeth absent, mandibular teeth slightly longer anteriorly than posteriorly. Body cylindrical and slender, almost consistently thick from the neck to the tail base. Head small and relatively narrow (Pileus length 7.05 mm; Pileus length/parietal width = 1.9), indistinct from the neck, with pointed oblique shape at the anterior side. Head scales, from the rostral to the posterior margin of parietals, including the temporals and post-temporals, not keeled. Rostral scale enlarged, extending backwards, intruding between the internasals, with a length of 1.42 mm, height of 1.15 mm and width of 1.79 mm. Rostral bordered by two upper labials, two nasals and two internasals. Nostrils situated on undivided nasal scales with a distance of 1.32 mm between them. Internasals of triangular shape. The internasal suture (length 0.41 mm) is almost half the length of the prefrontal suture (length 0.78 mm). A square-shaped loreal scale at either side in contact with 2nd and 3rd upper labials. The eyes are small with circular pupil of 1.31 mm in diameter. One preocular scale on both sides. Two postocular scales on the left side, and one on the right side. Frontal scale situated between the supraoculars with a length of 2.66 mm and width of 2.02 mm. Six squarish upper labials, 3rd and 4th in contact with eye. Two large parietal scales with a length of 3.11 mm and width of 3.73 mm. One temporal and one post-temporal scale on each side. Eight lower labials. Four pairs of lower labials in contact with anterior inframaxillars on each side. One gular scale positioned between posterior inframaxillars and in contact with anterior inframaxillars. Anterior and posterior inframaxillars lengths are 2.04 mm and 1.84 mm, respectively. Eleven dorsal and temporal scales surrounding the posterior margin of the parietals. Fifteen dorsal scale rows at mid-body. Two-hundred and one ventrals and 61/61 + 1 subcaudals, including a conical scale at the tail tip.
Colouration in alcohol of the head, neck and first dorsal scales (12 dorsal scales at midbody from between the parietal to end of pattern) is glossy black, including the upper and lower labials. The gulars are partially black with whitish background and the first ventrals are laterally black, while the middle ventral side is whitish. The ground colour of the dorsum is grey-yellowish-brown without maculation, while the ventral side is lighter, white.
Variation. Dentition features of the paratypes are similar to that of the holotype. Paratypes are similar to the holotype in most morphological characters (Table 2 Habitat and ecology. Specimens of R. dayanae sp. nov. were mostly found as road kill. The snake's secretive lifestyle, its remote habitat and its presence, as far as we know, only in a strict nature reserve, hinder a full evaluation and detailed knowledge of its ecological preferences. Its biology, and that of the other Rhynchocalamus species, is poorly known. Rhynchocalamus dayanae sp. nov. is a reclusive, fossorial, ground-dwelling snake. It is mainly nocturnal, and has been observed active only at night on the ground and occasionally inside human households (TAU.R17093 was killed inside a house in the town of Mitzpe Ramon). Observations and collections of these snakes were carried out during the evenings and nights, at elevations between 700 and 1,000 m a.s.l. (G Vine, pers. comm., 2014). The species inhabits arid and rocky or stony steppes and sparsely vegetated areas, but is not found on sand ( Fig. 6; Bar & Haimovitch, 2013;Werner, 2016). The diet consists predominantly of small insects such as ants, crickets and locusts (recorded for R. melanocephalus; Disi et al., 2001).
Distribution. Rhynchocalamus dayanae sp. nov. comprises to date only six specimens, collected from the Negev Mountain area in southern Israel, from the town of Mitzpe Ramon south-west up to the Egyptian border (Fig. 1). Most of the specimens of this newly-described species were observed and collected as road kills at the sides of main roads no. 40 (north-south direction) and no. 171 (from Mitzpe Ramon westwards to the Egyptian border). It probably also occurs on the Egyptian side of the border, but no specimens of R. dayanae sp. nov. from Egypt are known, although some specimens of R. melanocephalus are known from the Sinai Peninsula ( Fig. 1; HUJ.R8856, HUJ.R8885, TAU.R12494); at the moment, therefore, it is considered endemic to Israel. discovery may stem from the reclusive nature of these snakes, and because more specimens from southern Israel have been collected only recently. Subsequently, there are a mere six specimens known within the newly-discovered species in natural history collections, only four of which could be successfully used for genetic analyses (Tables S1 and S3). This new species from the Negev Mountain region in Israel is morphologically and genetically distinct, with a limited distribution from around the town of Mitzpe Ramon westward to the Egyptian border (Fig. 1). The specimens examined in this study reveal R. dayanae sp. nov. and R. melanocephalus to be allopatric species, separated by a distance of 15 km (between specimens HUJ.R22054, R22055 in locations 6-7 in Fig. 1 and HUJ.R21704, location 5). However, due to the low number of specimens from the area and of samples available, this assessment may change with the accumulation of more data. The molecular results reveal that R. dayanae sp. nov. is phylogenetically closely related to R. arabicus from Oman rather than to the geographically adjacent populations of R. melanocephalus in Israel northwards to Turkey. This relationship is surprising, as the two former species are located at a distance of about 2,500 km from one another (Fig. 1), whereas the closest collecting localities of R. melanocephalus and R. dayanae are just 15 km apart. As Rhynchocalamus snakes are rare and observations of them are scarce (specimens were not observed in several surveys in Arabia; e.g., Gasperetti, 1988;Schätti & Gasperetti, 1994;Schätti & Desvoignes, 1999;Van der Kooij, 2001;Gardner, 2013), it is plausible that a connection exist between populations, especially in desert habitats (as suggested in Gasperetti, 1988).
Specimens from the Sinai Peninsula (Table S3; HUJ.R8856; HUJ.R8885; TAU.R12494; Baha El Din, 2006) present a morphological combination of meristic characters similar both to R. dayanae sp. nov. (eight lower labials, 67%; four lower labials in contact with the anterior inframaxillars, 100%) and to R. melanocephalus (one post-ocular scale; two posttemporal scales; 13 dorsal and temporal scales surrounding the parietals). The Sinai specimens' colouration is consistent with that of R. melanocephalus (gular scales and the upper and lower labials white). At this point, due to the low sample size, we cannot determine the taxonomic classification of the Sinai specimens and therefore cannot reject the presence of R. dayanae sp. nov. in the Sinai Peninsula. Because these specimens were kept for decades in formalin, the possibility of a genetic study was not possible.
The distribution of R. melanocephalus in Israel is not well defined, as it is mostly known from the Ecotone region in the northern Negev and northwards, and with small marginal population in the Negev desert (Arbel, 1984;Werner, 1988;Bouskila & Amitai, 2001;Bar & Haimovitch, 2013;Werner, 2016). An early Israeli reptile guide (Barash & Hoofien, 1956) noted that the southern desert populations of R. melanocephalus from the area of Sde Boker in the Negev displayed some morphological differences from the other, northern, specimens of the species, and may represent a new subspecies. We found that the black pattern of the head and the dorsal scales of the Negev populations in Israel (HUJ.R3653, HUJ.R8921, HUJ.R21833, R21834, HUJ.R22054, R22055, R22056, R22057, R22058; locations 6-12 in Fig. 1), adjacent to Sde Boker, do indeed differ from the Mediterranean populations and those from the desert areas in the Sinai Peninsula and Jordan. However, these differences are the sole indicators for differentiation, as their mensural and meristic traits are within the intraspecific variation of R. melanocephalus (Tables S1 and  S3), and they are also genetically closely related (Fig. 2). An interesting finding arising from the phylogenetic analyses is that of the differentiation of the sole sample from Mt. Hermon in Israel (HUJ.R20967), which is otherwise morphologically very similar to other R. melanocephalus in terms of scalation and measurements. As no additional specimens from proximate localities in Lebanon or Syria are available, we are unable to account for this differentiation with the data at hand. In addition, sexual dimorphism is noted within R. melanocephalus. Females, in general, have a larger body size and a higher number of ventrals; males have relatively larger head dimensions and a higher count of sub-caudal scales (see also Werner, 2016).
The distinct lineage of R. satunini from Turkey and Iran, and the morphological revision, support its discrete specific status (as recommended by Reed & Marx, 1959;Baran, 1977;Avci et al., 2015;Šmíd et al., 2015). This molecular and morphological assessment contradicts Franzen & Bischoff (1995) contention that this species shows no morphological differences from R. melanocephalus. That said, the lack of sufficient sampling in Iran and the southern Caucasus, along with the species' disjunct distribution ( Fig. 1 and Fig. S4), and the relatively high genetic diversity within the species, suggest that additional lineages may well be awaiting description.
The enigmatic R. arabicus, one of the rarest snakes in the world, represents its own lineage, based on the sole sample available for the phylogenetic analyses. This snake is known only from one specimen (the holotype) as the recently sampled individual was not collected (Šmíd et al., 2015). Ecological data, such as habitat preferences or biology, are also lacking for this snake. The newly-sampled specimen from Oman was observed during the night, resting between stones in a vegetated wadi with permanent water (Šmíd et al., 2015). This observation accords with the nocturnal behaviour of other Rhynchocalamus species.

Biogeographic assessment of Rhynchocalamus
The estimated divergence times in this study present the first calibrated phylogeny of Rhynchocalamus to date. The estimated times indicate a very old, Late-Oligocene, divergence between Lytorhynchus and Rhynchocalamus at ca. 26 Mya ( Fig. 2A). Lytorhynchus is comprised of six species (Uetz & Hošek, 2016) and is mostly distributed in south-western Asia with a single species also found in Africa (i.e. L. diadema; Sindaco, Venchi & Grieco, 2013). Lytorhynchus is sympatric with Rhynchocalamus in south-western Iran and the Middle East, though they predominantly occur in different habitats (i.e., Lytorhynchus is a sand dweller whereas Rhynchocalamus occurs on heavier soils; Khan, 1996;Disi et al., 2001;Bar & Haimovitch, 2013). The close phylogenetic relationship between the two mostly south-west Asian Lytorhynchus and Rhynchocalamus may help to elucidate the biogeography of the latter genus.
The Late-Oligocene divergence between Lytorhynchus and Rhynchocalamus may have resulted from the changing landscapes and habitat fragmentation around the contact region of the Paratethys and Mediterranean seas with the Indian Ocean in south-western Asia (i.e., today's south-eastern Turkey and northern Syria). This region is where the collision of the African-Arabian plate with the Eurasian landmass took place during the Early Miocene, 16-18 Mya (Rögl, 1999;. The landscape of this region underwent continual changes during the Oligocene and Miocene due to several environmental changes, such as the tectonic motions of the Arabian and African plates resulting in the closure or opening of marine seaways, sea level fluctuations, the uplift of mountain ridges and the creation of regional fault systems (Rögl, 1999;Griffin, 2002;Inwood et al., 2009). The presence of two Levantine species, R. melanocephalus and R. dayanae sp. nov., and another in Arabia, R. arabicus, may allude to the ancestral range of Rhynchocalamus in the Levant (Fig. 1).
Radiation within Rhynchocalamus began with the separation of the northern R. satunini from Iran and Turkey during the Middle Miocene, approximately 15 Mya ( Fig. 2A). This diversification is likely to have been the result of a vicariance event resulting from the severing of the Gomphotherium land bridge in south-eastern Turkey, as was also suggested for other reptilian taxa such as the lacertid genus Apathya (Kapli et al., 2013). This Levantine land bridge, resulting from the African-Arabian-Eurasian tectonic connection, linked Africa and Eurasia ca. 16-18 Mya (Rögl, 1999;Inwood et al., 2009). The Gomphotherium land bridge had major biogeographical implications, enabling faunal dispersal from Africa to Eurasia and vice versa (Tchernov, 1992;Wüster et al., 2007;Pook et al., 2009;Zhou et al., 2012). The marine connection between the Mediterranean and Paratethys seas to the Indian Ocean, severing this land bridge, is hypothesized to have been re-established in the Middle Miocene, around 15-16 Mya (Rögl, 1999;Popov et al., 2004;. The re-opening of the marine connection and the disconnection of Eurasia from Africa-Arabia may have been responsible for the separation of the northern R. satunini from the ancestral, southern, population of the other Rhynchocalamus taxa. The Middle Miocene diversification of R. arabicus, R. dayanae sp. nov. and R. melanocephalus is difficult to interpret due to the paucity of known localities and samples in the Arabian Peninsula connecting R. arabicus to the other two species. The divergence between R. arabicus and R. dayanae sp. nov., ca. 9 Mya, may be related to several environmental events in the Arabian region. It has been hypothesized that continuous Middle-Late Miocene tectonic motions of the Arabian plate caused geological instability in north-western Arabia, leading to the creation of the Aqaba-Levant transform and periodic volcanic activity (Bosworth, Huchon & McClay, 2005). In addition, a temporal land connection of halite deposits existed between Africa and Arabia in the Red Sea, which later became submerged with the expansion and refill of the Red Sea (∼14-10 Mya; Richardson & Arthur, 1988;Rögl, 1999;Bosworth, Huchon & McClay, 2005). These events might have caused habitat fragmentation and facilitated radiation within Rhynchocalamus in the southern Levant and Arabia. The Middle Miocene global climate change resulted in an aridification process of the mid-latitude continental regions, triggering the expansion of sand areas in Arabia and North Africa (Flower & Kennett, 1994;Le Houérou, 1997;Edgell, 2006). The formation and progression of sandy habitats probably promoted vicariance and/or isolation of hard-substrate specific taxa, such as Rhynchocalamus, restricting the snakes to their current distribution, similar to what has been suggested for other reptiles, such as the geckos of the genus Stenodactylus (Metallinou et al., 2012), the Arid clade of