Resolving complex phylogeographic patterns in the Balkan Peninsula using closely related wall-lizard species as a model system.

The Balkan Peninsula constitutes a biodiversity hotspot with high levels of species richness and endemism. The complex geological history of the Balkans in conjunction with the climate evolution are hypothesized as the main drivers generating this biodiversity. We investigated the phylogeography, historical demography, and population structure of closely related wall-lizard species from the Balkan Peninsula and southeastern Europe to better understand diversification processes of species with limited dispersal ability, from Late Miocene to the Holocene. We used several analytical methods integrating genome-wide SNPs (ddRADseq), microsatellites, mitochondrial and nuclear DNA data, as well as species distribution modelling. Phylogenomic analysis resulted in a completely resolved species level phylogeny, population level analyses confirmed the existence of at least two cryptic evolutionary lineages and extensive within species genetic structuring. Divergence time estimations indicated that the Messinian Salinity Crisis played a key role in shaping patterns of species divergence, whereas intraspecific genetic structuring was mainly driven by Pliocene tectonic events and Quaternary climatic oscillations. The present work highlights the effectiveness of utilizing multiple methods and data types coupled with extensive geographic sampling to uncover the evolutionary processes that shaped the species over space and time.