Molecular Footprints of Quaternary Climate Fluctuations in the Circumpolar Tundra Shrub Dwarf Birch

The Arctic tundra biome is undergoing rapid shrub expansion ('shrubification') in response to anthropogenic climate change. During the previous ~2.6 million years, glacial cycles caused substantial shifts in Arctic vegetation, leading to changes in species' distributions, abundance and connectivity, which have left lasting impacts on the genetic structure of modern populations. Examining how shrubs responded to past climate change through genetic data reveals the demographic dynamics that shaped their current diversity and distribution and sheds light on the resilience of Arctic shrubs. Here we test scenarios of Quaternary demographic history of dwarf birch species (Betula nana L. and Betula Glandulosa Michx.) using Single Nucleotide Polymorphism (SNP) markers obtained from RAD sequencing and approximate Bayesian computation. We compare the timings of modelled population events with ice sheet reconstructions and other paleoenvironmental information to untangle the impacts of alternating cold and warm periods on dwarf birch. Our best supported model suggested that the species diverged in the Mid-Pleistocene Transition as glaciations intensified. We found support for a complex history of inter- and intraspecific divergences and gene flow, and secondary contact occurred during both ice sheet expansion and retreat. Our spatiotemporal analysis suggests that the modern genetic structure of dwarf birch results from transitions in climate between glacials and interglacials, with ice sheets acting alternatively as a barrier or an enabler of population mixing. Tundra shrubs may have had more nuanced responses to past climatic changes than phylogeographic analyses have often suggested, with implications for future eco-evolutionary responses to anthropogenic climate change.