Teasing apart the joint effect of demography and natural selection in the birth of a contact zone

AbstractVast 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 genome-wide 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 contact zones.Significance StatementUnderstanding how past climatic events, human actions and evolutionary forces contributed to the present distribution of genetic diversity is crucial to predict their reaction to the current climate crisis. Vast distribution shifts induced by past environmental changes, local ecological processes, natural selection and human transfers contributed to the current distribution of Norway spruce across Northern Europe. Genome-wide polymorphisms from thousands of individuals show that Scandinavia was recolonized after the Last Glacial from both south and north. This two-pronged recolonization established a contact zone between two genetic clusters that matches the limit between two major climate zones. The contact zone is shaped and maintained by natural selection on a large number of loci that form blocks of co-adapted loci spread genome-wide.