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

This study reveals that while three boreal tree species in Sweden share a common north-south genetic contact zone aligned with environmental gradients, they exhibit distinct patterns of gene flow and local adaptation architectures, ranging from broad genomic distribution to specific chromosomal inversions, illustrating how migration and selection jointly shape genomic landscapes in shared environments.

The Gist

Imagine Sweden as a giant, natural laboratory where three different types of trees—Norway Spruce (Picea abies), Silver Birch (Betula pendula), and Scots Pine (Pinus sylvestris)—are living out a dramatic story of migration, mixing, and survival.

Thousands of years ago, during the last Ice Age, these trees were forced to hide in "safe houses" (refugia) at the very bottom of Europe or tucked away in the far north. When the ice melted, they marched back up the country to reclaim their homes.

Here is the simple story of what happens when these two marching armies meet in the middle of Sweden, and how the scientists figured out what's going on.

1. The Great Meeting in the Middle

Imagine two armies marching toward each other on a battlefield. One comes from the South, the other from the North. When they meet, they don't just stop; they mix. This meeting point is called a Contact Zone.

In this study, the scientists found that all three tree species have a meeting zone right in the middle of Sweden (roughly between 60° and 63° latitude). Interestingly, this isn't just a random spot on the map. It lines up perfectly with a climate change. North of this line, it's cold and short-season; south of it, it's warmer and longer-season. It's like a natural "weather wall" where the trees have to decide: "Do I belong to the cold team or the warm team?"

2. Three Different Personalities

Even though they are all meeting in the same place, the three tree species handle this mixing very differently. Think of them as three different types of people at a party:

• The Norway Spruce (The Strict Traditionalist):
  This tree is very picky. It has strong "family ties" and doesn't mix its genes easily with the other side. It keeps a clear distinction between its Northern and Southern cousins. It's like a person who only talks to people from their own hometown, even at a big party. Because it's so strict, the "border" between the two groups is sharp and well-defined.
• The Silver Birch (The Diplomat with a Secret):
  This tree is in the middle. It mixes more than the Spruce, but it still keeps some boundaries. However, the scientists found something amazing: the Birch has a genetic "super-chip" (a chromosomal inversion) on one of its chromosomes. Think of this as a special backpack that holds all the important survival tools for the cold climate. Because this backpack is zipped shut (due to the inversion), the trees can't accidentally lose these tools when they mix with the warm-climate trees. It's a clever trick to stay adapted without getting confused.
• The Scots Pine (The Social Butterfly):
  This tree is the most relaxed. It mixes its genes freely with everyone. It's like the person at the party who talks to everyone and forgets who they were originally. Because it mixes so much, there is no sharp line between North and South. Instead, the changes happen very gradually, like a smooth fade in a sunset rather than a hard cut.

3. The Battle of Nature vs. Nurture

The big question was: Why do these lines exist? Is it just because the trees haven't mixed enough yet (history), or is nature actively forcing them to stay separate (selection)?

The answer is Nature (Selection).
The scientists found that the trees are constantly being "tested" by the environment.

• If a Southern tree tries to grow in the North, it might freeze or fail to reproduce.
• If a Northern tree tries to grow in the South, it might grow too fast and get damaged by frost.

So, even though pollen flies everywhere (gene flow), natural selection acts like a bouncer at a club, kicking out the trees that don't fit the local climate. This keeps the two groups distinct despite their constant mixing.

4. The "Polygenic" Puzzle

How do the trees know how to adapt?

• Spruce uses a "scattered" approach: It has thousands of tiny genetic tweaks spread all over its genome, like having a million tiny screws holding a machine together.
• Birch uses a "concentrated" approach: It relies heavily on that special "backpack" (the inversion) on one chromosome to hold all the necessary parts together.
• Pine uses a "subtle" approach: It relies on tiny, almost invisible shifts in many genes. It's like a choir where everyone sings a slightly different note, but together they create a perfect harmony that fits the environment.

Why Does This Matter?

This study is like a crystal ball for the future. As the climate changes and gets warmer, these trees will have to move or adapt again.

• The Spruce might struggle to move fast because it's so picky.
• The Birch might be able to adapt quickly because of its clever genetic "backpack."
• The Pine might be the most flexible, able to shift its traits gradually as the climate changes.

In a nutshell: This paper shows us that even when species live in the same place and mix their genes, the environment acts as a powerful editor, keeping them distinct. By comparing these three trees, we learn that there is no single "right way" to survive; nature uses different strategies (strict boundaries, genetic backpacks, or subtle shifts) to help life survive in a changing world.

Technical Summary

1. Problem Statement and Objectives

Post-glacial recolonization of Fennoscandia following the Last Glacial Maximum (LGM) created secondary contact zones where distinct genetic lineages (originating from southern and northern refugia) met and hybridized. While such zones are well-documented in animals and the Norway spruce (Picea abies), it remains unclear how these zones are maintained in other major boreal tree species with different life histories and dispersal capabilities.

The study addresses three primary objectives:

1. Determine Contact Zone Existence: Investigate whether Betula pendula (silver birch) and Pinus sylvestris (Scots pine) exhibit a secondary contact zone separating northern and southern genetic lineages similar to P. abies.
2. Identify Selection Signatures: Test for genomic signatures of natural selection that maintain genetic divergence despite ongoing gene flow.
3. Compare Adaptive Architectures: Compare the genomic architecture of local adaptation across the three species to understand how varying migration regimes influence the maintenance of adaptive variation.

2. Methodology

The researchers conducted a comparative genomic analysis of three species sampled along a latitudinal gradient in Sweden (60°N to 69°N).

• Data Acquisition:

  • Species: Picea abies (154 trees, 12 populations), Betula pendula (150 trees, 12 populations), and Pinus sylvestris (351 trees, 10 populations).
  • Sequencing: P. abies and B. pendula utilized exome capture data; P. sylvestris used Genotyping-by-Sequencing (GBS).
  • SNP Datasets: Final neutral datasets contained ~415k SNPs (P. abies), ~23k SNPs (B. pendula), and ~66k SNPs (P. sylvestris).

• Population Structure Analysis:

  • PCA & Clustering: Principal Component Analysis (PCA) and ADMIXTURE/TESS3 (K=2) were used to identify genetic clusters and admixture proportions.
  • Isolation-by-Distance (IBD): Genetic distance ($F_{ST}/(1-F_{ST})$) was regressed against geographic distance to estimate neighborhood size ($N_S$) as a proxy for dispersal potential.
  • Migration Surfaces: Effective migration surfaces were estimated using FEEMS to visualize barriers to gene flow.

• Selection and Adaptation Detection:

  • Barrier Detection: The program RIDGE (using Approximate Bayesian Computation) was employed to classify loci as barriers to gene flow, testing 14 demographic/genomic models (e.g., Isolation with Migration vs. Secondary Contact).
  • Genome Scans: Two approaches identified outlier loci:
    1. pcadapt: Individual-based PCA outlier detection.
    2. XTX (Bayenv2): Population-based differentiation corrected for structure.
  • Genotype-Environment Association (GEA): Bayenv2 was used to correlate allele frequencies with 19 bioclimatic variables.
  • Cline Analysis: The hzar R package fitted sigmoid and stepped cline models to allele frequencies to estimate cline center ($c$) and width ($w$). Selection strength ($s$) was estimated using the formula $s \approx N_S / w^2$.

• Functional Annotation: Candidate genes were annotated using GO terms and InterProScan to identify enriched biological processes.

3. Key Results

A. Shared Contact Zone with Species-Specific Structure

• Genetic Clustering: All three species exhibited two distinct genetic clusters (Northern and Southern) separated by a contact zone between 60°N and 63°N.
• Environmental Alignment: This genetic contact zone coincides precisely with the transition between two major climatic zones in Sweden (primarily driven by temperature gradients).
• Differential Gene Flow:
  • P. abies: Showed the strongest genetic structure and the most restricted gene flow. The contact zone is narrow, and differentiation is maintained even outside the immediate contact zone.
  • B. pendula: Exhibited intermediate structure. Gene flow is constrained primarily around the contact zone.
  • P. sylvestris: Displayed the weakest structure and highest gene flow (largest neighborhood size). Genetic differentiation is low, and the contact zone is more diffuse.

B. Genomic Architecture of Local Adaptation

The study revealed distinct adaptive architectures corresponding to the species' migration regimes:

1. Diffuse Architecture (P. abies): Adaptive loci are broadly distributed across the genome. Selection acts on many loci with small effects, consistent with lower gene flow allowing selection to counteract migration across the genome.
2. Concentrated Architecture (B. pendula): Strikingly, the vast majority of adaptive signals (pcadapt, XTX, and GEA candidates) are clustered within a ~9 Mb region on Chromosome 1. Linkage disequilibrium (LD) analysis confirmed this region is a chromosomal inversion. This inversion likely captures adaptive alleles and suppresses recombination, protecting them from gene flow.
3. Transient/Weak Architecture (P. sylvestris): Adaptive signals are weak and scattered. Due to high gene flow, adaptation relies on subtle allele frequency shifts across many loci (polygenic adaptation) rather than strong outliers.

C. Cline Analysis and Selection Strength

• Clinal Patterns: Candidate SNPs were significantly enriched for clinal alleles (sharp frequency transitions) in all species.
• Alignment: In P. abies and B. pendula, cline centers aligned perfectly with the environmental transition zone (60–63°N). P. sylvestris showed more variable cline centers, likely due to high gene flow accelerating cline decay.
• Selection Strength: Estimated selection coefficients ($s$) were generally weak ($10^{-5}$ to $10^{-2}$), supporting a polygenic basis for adaptation. However, P. sylvestris (highest gene flow) required the strongest selection coefficients to maintain clines, while P. abies required the weakest.

D. Functional Insights

• Enriched Gene Ontology (GO) terms across all species included photoperiodism, flowering time, growth regulation, and stress responses, confirming that temperature-related selection drives local adaptation.
• Overlap in specific candidate genes between species was low, likely due to technical limitations (different sequencing strategies) and the polygenic nature of traits (different genes achieving similar phenotypes).

4. Significance and Conclusions

• Validation of Suture Zones: The study confirms that a shared secondary contact zone exists across multiple boreal tree species in Sweden, driven by post-glacial recolonization from distinct refugia and maintained by environmental selection.
• Migration-Selection Balance: The research provides empirical support for theoretical models predicting that the genomic architecture of local adaptation shifts based on migration rates:
  • Low migration $\rightarrow$ Diffuse, genome-wide adaptation.
  • Intermediate migration $\rightarrow$ Concentrated adaptation via inversions.
  • High migration $\rightarrow$ Weak, polygenic adaptation relying on small allele frequency shifts.
• Climate Change Implications: Understanding how these species maintain local adaptation despite gene flow is critical for predicting their resilience to rapid climate change. The findings suggest that while P. abies may rely on specific genomic regions, P. sylvestris possesses a flexible, polygenic capacity to adapt, though it faces higher risks of "migration load" if environmental shifts outpace evolutionary responses.
• Methodological Contribution: The comparative approach demonstrates the power of integrating demographic modeling, cline fitting, and functional genomics to dissect the complex interplay between history, selection, and gene flow in forest trees.