Teasing apart the joint effect of demography and natural selection in the birth of a contact zone
Vast population movements induced by recurrent climatic cycles have shaped the genetic structure of plant species. This is especially true in Scandinavia that was repeatedly glaciated. During glacial periods trees were confined to refugia, south and east of the ice sheet, from which they recolonized Scandinavia as the ice melted away. This multi-pronged recolonization led to large contact zones in most species. We leverage large genomic data from 5000 trees to reconstruct the demographic history of Norway spruce (Picea abies) and test for the presence of natural selection during the recolonization process and the establishment of the contact zone. Sweden is today made up of two large genetic clusters, a southern one originating from the Baltics and a Northern one originating from Northern Russia. The contact zone delineating these two clusters closely matches the limit between two major climatic regions. This suggests that natural selection contributed to the establishment and the maintenance of the contact zone. To test this hypothesis we first used Approximate Bayesian Computation; an Isolation-with migration model with genomewide linked selection fits the data better than a purely neutral one. Secondly, we identified loci characterized by both extreme allele frequency differences between geographic regions and association to the variables defining the climatic zones. These loci, many of which are related to phenology, form clusters present on all linkage groups. Altogether, the current genetic structure reflects the joint effect of climatic cycles, recolonization and selection on the establishment of strong local adaptation and con-tact zones.