There is increasing evidence that most parapatric cryptic/sister taxa are reproductively compatible across their areas of contact. Consequently, the biological species concept, which assumes absence of interbreeding, is becoming a not so effective criterion in evolutionary ecology. Nevertheless, the few parapatric sister taxa showing complete reproductive barriers represent interesting models to study speciation processes and the evolution of reproductive isolation. In this study, we examined contact populations in northwestern Italy of two butterfly species, Zerynthia polyxena and Z. cassandra, characterized by different genitalic morphotypes. We studied levels of divergence among 21 populations distributed from Sicily to France using three genetic markers (the mitochondrial COI and ND1 genes and the nuclear wingless gene) and genitalic geometric morphometrics. Moreover, we performed species distribution modelling to estimate different climatic requirements of Z. polyxena and Z. cassandra. We projected climatic data into glacial maximum scenarios in order to verify if and to which extent glacial cycles could have contributed to speciation processes. Genetic and morphometric analyses identified two main groups. All specimens showed a concordant pattern of diversification, including those individuals sampled in the contact area. Haplotype distribution and climatic models showed that during glacial maxima both species experienced a strong range contraction and presumably remained separated into different microrefugia in southern France, in the Italian Peninsula and on the islands of Elba and Sicily. Long term separation was probably favoured by reduced dispersal ability and high phylopatry, while genitalic diversification probably favoured interbreeding avoidance. Conversely, the aposematic wing pattern remained almost identical. We compared our results with those obtained in other species and concluded that Z. polyxena and Z. cassandra represent a valuable model in the study of speciation.