Owing to its temperature dependence and low vagility, the asp viper (Vipera aspis) is an interesting model species to study the effects of Pleistocene climatic fluctuations on vertebrate genomes. We genotyped 102 specimens from the whole Italian distribution range at three mitochondrial DNA regions (2278 characters, total) and six microsatellite DNA loci (Short Tandem Repeats, STR). The molecular phylogeny was constructed according to Bayesian, Neighbour Joining, Maximum Parsimony and Maximum Likelihood procedures. All methods grouped individuals of the three morphological subspecies (V. a. aspis, V. a. francisciredi, V. a. hugyi) into five different haploclades. Specimens assigned to hugyi clustered in two highly differentiated clades, one being sister group to the complex comprising the second clade of hugyi (i.e., a paraphyletic status), plus two clades of francisciredi. The Bayesian clustering of the STR variability disclosed only two groups, the first including aspis and francisciredi, the second all hugyi. Introgressive hybridization and capture of francisciredi-like lineages in the hugyi mitochondrial genome were suggested to explain the discordance between mitochondrial and nuclear data. The phylogeographic pattern was compatible with population contractions in three glacial refuges. Plausibility of derived hypothesis was checked using coalescence simulations as post hoc tests. Long-term drift and serial founder effects, rather than selection, appeared the main factors affecting the genetic make-up of the Italian asp viper.