Genomic diversity of British native oaks: species differentiation, hybridisation and triploidy

This study utilizes whole-genome sequencing of 418 British oak individuals to reveal distinct ecological niches and geographic distributions for *Quercus robur* and *Q. petraea*, extensive hybridization biased toward *Q. robur* introgression, the presence of rare triploid individuals with superior growth rates, and specific genomic islands of differentiation despite widespread gene flow.

The Gist

Imagine the British countryside as a massive, ancient library. For centuries, people have been trying to sort the books on the shelves into two distinct categories: Pedunculate Oak (Quercus robur) and Sessile Oak (Quercus petraea).

For a long time, librarians (botanists and foresters) were confused. The books looked so similar, and the shelves were so mixed up, that they couldn't tell which was which. Some said, "They're all just one big messy family!" while others insisted, "No, there are two distinct species!"

This new study is like hiring a team of super-sleuths with high-tech DNA scanners to walk through every aisle of this library. They didn't just look at the covers (leaves); they read the entire text of 418 different trees from 60 forests across Britain. Here is what they discovered, translated into everyday language:

1. The Great Geographic Divide (The "North-South" Split)

Think of Britain as a long hallway.

• The Pedunculate Oak is the "Sunshine Oak." It loves the warm, dry, southern, and eastern parts of the country (like London and the Midlands). It prefers soil that isn't too acidic, kind of like a plant that likes a nice, sunny garden bed.
• The Sessile Oak is the "Rainy Mountain Oak." It prefers the cool, wet, rugged, and acidic soils of the north and west (like Wales, Scotland, and Cornwall). It's the tough guy that thrives on steep hillsides.

The Analogy: Imagine a party where the "Sunshine Oak" is hanging out by the pool in the south, and the "Rainy Mountain Oak" is dancing in the rain on the roof in the north. They have their own zones, but they do meet in the middle.

2. The Great Mix-Up (Hybridization)

Here is the twist: These two oaks are terrible at keeping to themselves. They are constantly having "affairs."

• The Result: A huge number of trees are actually hybrids (mixed children of both parents).
• The Direction of the Mix: The mixing isn't fair. It's mostly the "Sunshine Oak" (Q. robur) sneaking its genes into the "Rainy Mountain Oak" (Q. petraea) family. It's like a popular celebrity moving into a small town and slowly changing the local culture, while the locals rarely move to the city.
• The "Resurrection" Theory: Scientists used to think the "Rainy Mountain Oak" was actually just the "Sunshine Oak" that got a few new genes and changed its look to survive the cold. But this study says, "Nope, they are still distinct species." They just mix a lot, but they haven't merged into one big blob.

3. The Genetic "No-Go" Zones

Even though these trees mix their DNA everywhere else, there are specific "fortresses" in their genetic code where they refuse to mix.

• The Analogy: Imagine two rival gangs. They might trade goods and hang out in the neutral zone, but they have a strict rule: "You can never cross the border into our secret clubhouse."
• The Findings: The researchers found these "clubhouses" on specific chromosomes (especially chromosome 2). These areas contain about 2,000 genes that keep the two species distinct. If these genes mixed, the trees might lose their ability to survive in their specific environments.

4. The "Three-Parent" Anomalies (Triploids)

Usually, trees have two sets of chromosomes (one from mom, one from dad), like a pair of shoes. But the researchers found five trees that had three sets (like a pair of shoes plus an extra left shoe).

• The Analogy: These are the "genetic oddballs." They are triploids.
• The Superpower: These five trees weren't just weird; they were super-growth machines. Even when you accounted for the weather and soil, these three-set trees grew significantly faster than the normal two-set trees. They are like the "super-soldiers" of the oak world, potentially useful for planting forests that need to capture carbon quickly.

5. The Growth Race

The study also looked at how fast these trees grew between 1990 and 2019.

• The Winner: The "Sunshine Oak" (Q. robur) grew faster on average.
• The Reason: It wasn't because they are genetically superior; it's because they live in warmer, friendlier places.
• The Twist: When the researchers looked at the environment, they realized the trees were just reacting to their neighborhood. If you put a "Rainy Mountain Oak" in a warm spot, it would grow fast too.

Why Does This Matter?

This study is like getting the first high-definition map of a territory that was previously blurry.

1. For Conservation: We now know exactly where the pure species live and where the hybrids are, helping us protect the unique genetic "fortresses" that keep them distinct.
2. For Forestry: If you want to plant a forest, you now know which oak to plant where. Don't plant the "Sunshine Oak" in the wet, rocky north; it won't thrive.
3. For the Future: Those five "super-growth" triploid trees might hold the secret to planting forests that grow faster and soak up more carbon dioxide to fight climate change.

In a nutshell: British oaks are a complex, intermingled family with a clear north-south split. They mix genes constantly, but they have secret genetic walls that keep them from becoming one single species. And hidden among them are a few "super-trees" that grow like weeds, offering a potential key to a greener future.

Technical Summary

1. Problem Statement

The two native British oak species, Quercus robur (pedunculate oak) and Quercus petraea (sessile oak), are ecologically and economically vital but morphologically difficult to distinguish, leading to frequent misidentification in forestry and ecological surveys. While previous studies have utilized leaf morphology, microsatellites, or limited genomic regions, there has been no comprehensive whole-genome analysis of these species across Great Britain. Key unresolved questions include:

• The precise geographic distribution and niche differentiation of the two species.
• The extent and directionality of hybridization and introgression between them.
• The genomic regions responsible for species maintenance despite gene flow.
• The existence and phenotypic impact of polyploidy (specifically triploidy) in natural populations.
• How genetic factors influence long-term growth rates in the context of environmental change.

2. Methodology

The study employed a large-scale population genomics approach combined with long-term phenotypic data.

• Sample Collection: Researchers resequenced the whole genomes of 418 individuals from 60 natural oak populations across Great Britain. These samples were sourced from the Forest Condition Survey (FCS), a long-term monitoring program (1987–2007) that originally tracked forest health against air pollution.
• Sequencing & Genotyping:
  • DNA was extracted from dried leaf samples and sequenced using Illumina NovaSeq 6000 (150-bp paired-end) to achieve ~20x coverage.
  • Variant Calling: Using the Q. robur reference genome, the pipeline (BWA-MEM, GATK) identified 146,995,863 SNPs and 28,275,691 indels.
  • Filtering: After rigorous quality control (removing transposable elements, filtering for depth, mapping quality, and missingness), a final dataset of 7,112,139 biallelic SNPs (Dataset 1) was generated for population structure analysis.
• Analytical Framework:
  • Ploidy Screening: Inbreeding coefficients ($F$) and allele balance analysis were used to detect triploids (expected 0.33/0.67 allele ratio vs. 0.5/0.5 for diploids).
  • Population Structure: fastSTRUCTURE and PCA were used to assign species and quantify admixture ($q$).
  • Differentiation: FST calculations and LEA (snmf) outlier detection were used to identify genomic regions under selection or maintaining species barriers.
  • Chloroplast Analysis: De novo assembly of chloroplast genomes to construct haplotype networks.
  • Niche Modeling: Principal Component Analysis (PCA) and Beta regression were used to correlate genetic admixture proportions with 19 bioclimatic variables, soil properties (pH, C:N), and topographic features.
  • Growth Analysis: Linear mixed models (LMM) analyzed Diameter at Breast Height (DBH) growth from 1990 to 2019, controlling for site, species, and ploidy.

3. Key Contributions

• First Whole-Genome Survey of British Oaks: This is the first study to utilize high-depth whole-genome sequencing across the entire geographic range of British oaks, moving beyond limited marker sets.
• Discovery of Triploidy: The study successfully identified and characterized rare triploid individuals in natural populations using genomic allele balance, a method applicable to future ploidy screening in trees.
• Quantification of Asymmetric Introgression: The study provides robust evidence of the directionality of gene flow, specifically the bias of Q. robur introgression into Q. petraea.
• Genomic Islands of Differentiation: Identification of specific genomic regions (particularly on Chromosome 2) that resist introgression, offering targets for understanding speciation mechanisms.
• Integration of Genomics and Long-Term Phenotyping: Uniquely links genomic data with 30 years of growth data to disentangle genetic vs. environmental drivers of tree performance.

4. Key Results

A. Species Distribution and Niche Differentiation

• Geographic Pattern: Q. robur is dominant in the south and east of Britain, while Q. petraea is dominant in the north and west. No pure Q. robur samples were found in Scottish sites.
• Niche Preferences:
  • Q. robur prefers warmer, more thermally variable environments with alkaline soils and flatter terrain.
  • Q. petraea prefers cooler, wetter environments with acidic soils, higher precipitation, and more complex, rugged topography.
• Hybridization: Extensive hybridization was found. Of the 418 samples, 61 were classified as hybrids/admixed.
  • Asymmetry: Introgression is heavily biased from Q. robur into Q. petraea.
  • F1 Hybrids: Only 4% of samples were F1 hybrids; the majority were backcrosses, predominantly into the Q. petraea genetic background.

B. Genomic Differentiation

• Outlier SNPs: The study identified 8,238 significantly differentiated SNPs (FDR < 0.01) between the species.
• Genomic Architecture: Differentiation is concentrated in "islands" on Chromosomes 2, 5, 7, and 11, while large portions of the genome (e.g., Chromosomes 3, 9, 12) show free exchange.
• Gene Content: Of the differentiated SNPs, 2,005 are located within 945 predicted genes. Many of these genes overlap with findings from previous European studies, suggesting conserved mechanisms of species differentiation.

C. Triploidy

• Prevalence: Five triploid individuals (1.2% of the sample) were identified.
• Composition: Two were Q. robur, two were Q. petraea, and one was a backcrossed hybrid.
• Growth Impact: Triploids exhibited significantly faster growth rates (16.6% increase in DBH) compared to diploids (12.5%), even after accounting for environmental variables and site effects.

D. Growth Rates and Environment

• Species Growth: Q. robur showed higher mean growth than Q. petraea and hybrids, but this was primarily driven by environmental factors (isothermality, temperature of the coldest quarter, and terrain roughness) rather than intrinsic genetic superiority.
• Environmental Drivers: Growth was positively associated with isothermality and cold-quarter temperatures, and negatively associated with terrain roughness.

E. Chloroplast Haplotypes

• Widespread sharing of chloroplast haplotypes was observed between species, indicating historical and ongoing chloroplast capture via hybridization. Five major haplotypes were shared across both species and hybrids.

5. Significance

• Ecological Insight: The study refutes the "resurrection model" (where Q. petraea spreads via pollen swamping of Q. robur) for the British context. Instead, the distribution suggests direct colonization by Q. petraea into the north/west driven by environmental selection, with Q. robur introgression occurring later.
• Conservation & Forestry: The identification of triploids with superior growth rates suggests potential for using these genotypes in high carbon-capture plantations.
• Species Management: The clear genomic definition of species boundaries and hybrid zones aids in accurate species identification for conservation and timber management, moving beyond unreliable morphological traits.
• Climate Adaptation: The findings highlight how specific environmental gradients (soil pH, temperature seasonality, topography) drive the maintenance of species boundaries in the face of extensive gene flow, providing a model for understanding tree adaptation to climate change.