Cold-blooded in the Ice Age: "refugia within refugia", inter-and intraspecific biogeographic diversification of European whipsnakes (Squamata, Colubridae, Hierophis).

In this work, we performed a biogeographic analysis with Bayesian binary MCMC (BBM) statistical dispersal-vicariance analysis (S-DIVA) and species distribution models (SDM) on three phylogenetically closely related Mediterranean whipsnakes (Hierophis gemonensis, H. carbonarius, H. viridiflavus), to investigate the pathways of their geographical diversification and locate putative refugial areas in the last glacial maximum (LGM). Our analysis suggests that the diversification processes between the studied species overall followed an east-west route, from eastern Greece to the Iberian Peninsula and continental France, highlighting a significant role of dispersal and vicariance processes at both inter- and intraspecific levels. In particular, the main lineage-splitting events between H. gemonensis, H. carbonarius and H. viridiflavus coincide with two events of vicariance, involving respectively eastern Greece and eastern Italy, and eastern Italy and western Italy, Iberian Peninsula and continental France. SDM analyses highlight the occurrence of multiple putative glacial refugia in the Balkans, Italy and southern France, which represent well the occurrence of distinct "refugia within refugia" in the main Mediterranean peninsulas. Furthermore, our results suggest how these refugia may have played an important role during the Quaternary climatic oscillations in shaping the current haplotype distribution of European whipsnakes.