Pleistocene origin and population history of a neoendemic alpine butterfly

Alpine environments underwent dramatic transformation during glacial–interglacial cycles, with the consequence that geographical, ecological and demographic changes of alpine populations provided the opportunity for formation of neoendemic species. Several biogeographical models have been proposed to account for the unique history of alpine populations, with different expectations of genetic divergence and speciation. The expanding alpine archipelago model proposes that alpine populations expand spatially and demographically during glacial events, dispersing between mountain ranges. Under this model, alpine populations are unlikely to diverge in isolation due to substantial interpopulation gene flow. In contrast, the alpine archipelago refuge model proposes that gene flow during glacial phases is limited and populations expand demographically during interglacial phases, increasing genetic isolation and the likelihood of speciation. We assess these models by reconstructing the evolutionary history of Colias behrii, a morphologically and ecologically distinct alpine butterfly restricted to the California Sierra Nevada. C. behrii exhibits very low genetic diversity at mitochondrial and nuclear loci, limited population structure and evidence of population expansion. C. behrii and Rocky Mountain C. meadii share identical mitochondrial haplotypes, while in contrast, nuclear data indicate common ancestry between C. behrii and Cascades Range Colias pelidne. The conflict in gene genealogies may be a result of recent expansion in North American Colias, but an isolation with migration analysis indicates that genetic patterns in C. behrii might result from differential introgression following hybridization. Based on the timing of population expansion and gene flow between mountain ranges, the expanding alpine archipelago model is supported in C. behrii.