AbstractIt has been generally acknowledged that glacial climates at the time of the Pleistocene altered the patterns of species distributions, prompting latitudinal and altitudinal distribution shifts in several species, including poikilothermic species commonly known for their thermal sensitivity. However, the historical phylogeographic patterns of such species have remained largely unknown. Here, we present the historical biogeographic, phylogenetic, and phylogeographic relationships of the Caucasian pit viper, G. h. caucasicus, based on two mtDNA (cyt b and ND4) and one nDNA (c-mos) genes. This pit viper represents the westernmost member of the Crotalinae subfamily in the Palearctic and occurs in a variety of habitats, from 30 m to 3,000 m above sea level. In Iran, it is distributed on the northern and southern slopes of the Alborz Mountains, rendering it a target for phylogenetic and phylogeographic studies of a terrestrial poikilothermic animal. Our study identified four Iranian lineages of G. h. caucasicus along the northeastern to northwestern slopes of the Alborz Mountains and southern Azerbaijan (Talysh Mountains). Diversification of the Iranian lineages highlights population expansion and subsequent isolation into four plausible refugial areas during the Quaternary paleo-climatic oscillations, confirmed by our molecular dating and historical biogeographic analyses. The results of coalescence-based simulations support the incursion of the species from northeastern Iran to the western end of the Alborz, and then toward Transcaucasia via two directions: northern and southern slopes of the Alborz Mountains. Furthermore, our results clearly implied that G. h. caucasicus should be elevated to species rank and further referred to as G. caucasicus (Nikolsky, 1916).