When South meets North: a joint contact zone coinciding with environmental gradients in three boreal tree species

Post-glacial recolonization of Fennoscandia created secondary contact zones in many species, offering opportunities to study how gene flow and selection contribute to their establishment and maintenance. Here, we analyse genomic data from three boreal tree species--Picea abies, Betula pendula, and Pinus sylvestris--sampled along a latitudinal gradient in Sweden. Despite differences in colonization timing and dispersal ecology, all three species exhibit north-south genetic structuring aligned with environmental gradients. Most notably, the two main genetic clusters within each species overlap in a shared contact zone, corresponding to the climatic transition between Swedens two major environmental zones. The extent and structure of the contact zone differ among species: P. abies shows stronger genetic structure and moderate gene flow, B. pendula exhibits intermediate differentiation and gene flow, and P. sylvestris displays the weakest structure with stronger gene flow. All three species also show genomic signatures of local adaptation, with distinct underlying architectures. In P. abies, adaptive loci are broadly distributed across the genome, while, strikingly, they are mostly found within an inversion on chromosome 1 in B. pendula. In P. sylvestris, local adaptation likely relies on subtle allele frequency shifts across many loci with weak signals. These patterns align with theoretical expectations for polygenic local adaptation under varying migration regimes. Our comparative approach demonstrates how gene flow and selection jointly shape genomic landscapes in shared environments and contributes to understanding local adaptation in forest trees, with implications for predicting species responses to climate change.