Imprints of multiple glacial refugia in the Pyrenees revealed by phylogeography and palaeodistribution modelling of an endemic spider

Mediterranean mountain ranges harbour highly endemic biota in islandlike habitats. Their topographic diversity offered the opportunity for mountain species to persist in refugial areas during episodes of major climatic change. We investigate the role of Quaternary climatic oscillations in shaping the demographic history and distribution ranges in the spider Harpactocrates ravastellus, endemic to the Pyrenees. Gene trees and multispecies coalescent analyses on mitochondrial and nuclear DNA sequences unveiled two distinct lineages with a hybrid zone around the northwestern area of the Catalan Pyrenees. The lineages were further supported by morphological differences. Climatic niche-based species distribution models (SDMs) identified two lowland refugia at the western and eastern extremes of the mountain range, which would suggest secondary contact following postglacial expansion of populations from both refugia. Neutrality test and approximate Bayesian computation (ABC) analyses indicated...