Effects of P leistocene climatic fluctuations on the phylogeography, demography and population structure of a high‐elevation snake species, T hermophis baileyi , on the Tibetan Plateau

Abstract

Aim Our aims were to investigate the geographical distribution, population history and demographics of the hot‐spring snake T hermophis baileyi , and to examine how the Pleistocene glacial stages and geomorphological configurations resulting from the uplift of the T ibetan– H imalayan orogen shaped the genetic diversity and evolutionary history of the species. Location T ibetan P lateau. Methods We surveyed potential habitat sites throughout T ibet. Sequence data were obtained for three mitochondrial markers ( ND2 , ND4 and cytb ; 2185 bp) in 181 individuals from 18 populations across the species' distribution range. We conducted population‐genetic ( samova ; baps ), phylogenetic (maximum likelihood; B ayesian), demographic ( T ajima's D ; Fu's F S ; mismatch distribution) and ecological analyses (logistic regression for presence/absence data) to examine the recent distribution, evolutionary history and diversification of the species. Results The geographical range of T . baileyi is a restricted area between the T ranshimalaya and the H imalaya, along the central part of the Y arlung T sangpo suture zone. The divergence times of sequences matched the G uxiang glaciation (300–130 ka) and the last glacial period (70–10 ka). We observed two genetically distinct lineages north of the Yarlung Tsangpo river. These lineages are known from previous work and are largely congruent with a division across the drainage divide of the Nyainqêntanglha rift. Haplotype grouping based on cluster analysis indicated the existence of riverine migration corridors between populations. Main conclusions The evolutionary history of T . baileyi is consistent with a model of range expansion from different refugia during interglacial and post‐glacial times. All lines of evidence suggest that Pleistocene glacial oscillations, triggered by tectonic activity during plateau uplift, have had the strongest influence on intraspecific differentiation. Physical barriers, such as high‐elevation mountains and the drainage system, however, also appear to have affected the population structure and distribution of this snake. Population clustering suggests that there has been long‐distance dispersal along riverine corridors.

Bibliographic Information