A reference genome assembly for Quercus canariensis Willd

This study presents the first chromosome-scale reference genome assembly for *Quercus canariensis*, generated using PacBio HiFi reads to produce two high-quality haplotypes with over 98% BUSCO completeness, thereby providing a foundational resource for future genomic and evolutionary research on this understudied white oak species.

The Gist

Imagine you are trying to understand how a massive, ancient library works. For centuries, scientists have been trying to read the books inside the library of the Algerian Oak (Quercus canariensis), a tree famous for its ability to survive hot, dry summers. But until now, they only had scattered, torn-up pages and rough drafts. They knew the library existed, but they couldn't see the whole story, let alone understand how the books were organized on the shelves.

This paper is the moment someone finally built a complete, high-definition map of that library.

Here is the story of how they did it, explained simply:

1. The Challenge: A Library in the Heat

European forests are getting hotter and drier. Trees are struggling. Scientists want to help by moving "super-trees" from warmer southern regions (like Spain and North Africa, where the Algerian Oak lives) to cooler northern regions. But to do this safely, we need to know exactly why these trees are so tough.

The problem? We didn't have the "instruction manual" (the genome) for this specific tree. We had manuals for its cousins (like the English Oak), but every tree has its own unique quirks. Without the specific manual for the Algerian Oak, we were flying blind.

2. The Tool: A Super-Powerful Camera

To build this manual, the researchers didn't use a standard camera. They used PacBio HiFi sequencing.

• The Analogy: Imagine trying to read a book that has been shredded into tiny confetti. A normal camera might take a picture of a few pieces and guess what the sentence says.
• The HiFi Method: This technology is like a super-camera that takes high-definition photos of long, continuous strips of the shredded paper. It captures the sentences in their original flow, making it much easier to put the puzzle back together without losing words.

They took a sample from a single, healthy oak tree in southern Spain and generated enough data to cover the tree's entire genetic code 39 times over. This "redundancy" is like reading the same book 39 times to make sure you didn't miss a single typo.

3. The Result: Two Perfect Copies

Trees are special because, like humans, they carry two sets of instructions (one from mom, one from dad). The researchers didn't just mash these two sets together into a blurry average; they separated them into two distinct, perfect copies called Haplotype 1 and Haplotype 2.

• The Assembly: They stitched the long strips of DNA together to form 12 giant chromosomes for each copy. Think of these as the 12 main volumes of an encyclopedia.
• The Quality: The map is incredibly complete.
  • 98% Complete: If you imagine the tree's "essential survival toolkit" as a set of 2,326 tools, they found 98% of them in the right place.
  • The "Missing" Pieces: Only a tiny fraction (about 1-3%) of the DNA was too small or messy to place on the main shelves. They put these in a "miscellaneous box" (Chromosome 0), but it's so small it doesn't matter much.
  • Accuracy: The text is so accurate that the error rate is less than 1 mistake in 100 million letters.

4. What's Inside the Manual?

Now that they have the map, they can start reading the chapters:

• The Gene Count: They found roughly 50,000 to 52,000 genes (the "recipes" for making the tree).
• The Junk Drawer (Transposable Elements): About half of the tree's DNA is made up of "jumping genes" or ancient viral leftovers that don't code for proteins but make up the structure of the genome. It's like the filler text and advertisements in a newspaper that take up half the space but are part of the paper's history.
• The Translation: They didn't just list the genes; they translated them into functions. They figured out which genes make enzymes (chemical workers), which ones handle stress, and which ones are unique to this tree.

5. Why Does This Matter?

This paper is the foundation for everything that comes next.

• Climate Change: Now that we have the map, we can hunt for the specific "genes of drought tolerance." We can see exactly which parts of the Algerian Oak's DNA allow it to survive in the scorching sun.
• Forest Management: If we want to plant these trees in France to save forests from heatwaves, we can now check if the local trees can safely "borrow" these tough genes through natural mixing (hybridization).
• Evolution: It helps us understand how this tree is related to its white oak cousins and how they evolved together.

The Bottom Line

Think of this paper as the moment someone finally finished building the Google Maps for the Algerian Oak. Before, we were driving through the forest with a paper map from a different country. Now, we have a satellite view, street names, and traffic reports for this specific tree. It opens the door to saving European forests by using the best tools nature has to offer.

Technical Summary

1. Problem Statement

European forests face increasing threats from climate change, specifically recurrent droughts and heatwaves, leading to the decline of deciduous tree populations, including European white oaks. While "assisted migration" (moving species from southern to northern regions) is a proposed management strategy to enhance adaptive capacity, it relies heavily on comprehensive genomic resources.

• The Gap: Quercus canariensis (Willd.), a sub-Mediterranean oak species with high drought and heat tolerance, is a promising candidate for northward reforestation and adaptive introgression into northern white oak populations. However, despite its ecological importance and potential for hybridization with species like Q. robur and Q. petraea, no complete, chromosome-scale reference genome assembly existed for this species.
• The Need: High-quality genomic resources are essential to clarify evolutionary relationships within the white oak species complex, identify genomic regions underlying climate adaptation, and detect signatures of selection.

2. Methodology

The authors generated a diploid, chromosome-scale reference genome using a combination of long-read sequencing and reference-guided scaffolding.

• Biological Sampling:
  • Primary Sample: A single adult Q. canariensis (~40–50 years old) was sampled in Cádiz, southern Spain, in a humid microhabitat.
  • TE Library Source: An admixed Q. petraea/pubescens individual from Switzerland was used to generate a transposable element (TE) consensus library, as Q. canariensis TEs were not yet characterized.
• Sequencing:
  • Technology: PacBio Revio System (HiFi long reads).
  • Data Volume: 1,820,767 reads yielding 30.7 Gb of data (~39X coverage).
  • Library Prep: High Molecular Weight (HMW) DNA was extracted, sheared to ~20 kb, and size-selected (>10 kb).
• Assembly Pipeline:
  • Assembler: hifiasm (v0.24.0) was used to produce a phased, diploid assembly resolved into two haplotypes.
  • Polishing/Redundancy Removal: purge_dups was used to remove redundant haplotigs.
  • Scaffolding: A reference-guided approach using RagTag (v2.1.0) was employed. The reference used was a newly assembled Q. robur (dhQueRobu3.1) genome derived from public PacBio HiFi and Hi-C data. This allowed the ordering and orientation of Q. canariensis contigs into 12 chromosomes.
  • Workflow: The entire process was managed via the Asm4pg pipeline (v2.0.0).
• Annotation:
  • Structural: The Eugene pipeline was used for gene prediction, trained with Helixer and supported by RNA-seq data and protein homology.
  • Functional: InterProScan, eggNOG-mapper, and the E2P2 pipeline were used for functional annotation and enzyme prediction.
  • Transposable Elements (TEs): A de novo TE library was built from the Swiss Q. petraea/pubescens sample using the REPET package (TEdenovo and TEannot), then used to annotate TEs in the Q. canariensis assembly.

3. Key Contributions & Results

A. Genome Assembly Statistics

The study produced two high-quality, haplotype-resolved assemblies:

• Haplotype 1: 816.0 Mb (845.4 Mb including unplaced sequences).
• Haplotype 2: 804.8 Mb (815.9 Mb including unplaced sequences).
• Contiguity: Both haplotypes were scaffolded into 12 chromosome-scale pseudomolecules (matching the species' 2n=24 karyotype).
  • N50: 20.1 Mb (Hap1) and 16.8 Mb (Hap2).
  • Unplaced Sequences: Only 3.48% (Hap1) and 1.36% (Hap2) of the assembly remained unplaced (assigned to "chr0").
• Completeness & Accuracy:
  • BUSCO (Embryophyta): 98.3% complete (Hap1) and 98.2% complete (Hap2).
  • LAI (LTR Assembly Index): 22.99 (Hap1) and 24.46 (Hap2), indicating excellent reconstruction of repeat-rich regions.
  • Telomeres: 19 and 16 telomeric motifs detected for Hap1 and Hap2, respectively.
  • Consensus Quality (QV): 72.84 for Hap2 (error rate ~5.19 × 10⁻⁸).
  • Heterozygosity: Estimated at 1.9%.

B. Gene Annotation

• Gene Count:
  • Haplotype 1: 65,629 total genes (51,882 protein-coding).
  • Haplotype 2: 55,146 total genes (46,482 protein-coding).
  • Note: When restricted to the 12 chromosomes, counts were 51,343 and 50,615 genes, respectively.
• Functional Annotation:
  • ~95% of predicted proteins received functional annotations via InterProScan, eggNOG, and E2P2.
  • Enzymatic functions were assigned to ~27–28% of proteins.

C. Transposable Elements (TEs)

• TEs account for 54.35% of the genome assembly.
• Identified 3,606 TE families and over 1.1 million TE copies.
• Mean sequence identity across families was 81.57%.

4. Significance

• First Reference Genome: This is the first chromosome-scale, diploid reference genome for Quercus canariensis, filling a critical gap in the genomic resources for the European white oak species complex.
• Climate Adaptation Research: The assembly provides the necessary foundation to investigate the genetic basis of drought and heat tolerance, which are crucial for predicting how these species will respond to climate change.
• Evolutionary Insights: By resolving two distinct haplotypes, the assembly facilitates the study of allelic diversity, introgression patterns, and evolutionary relationships between Q. canariensis and other white oaks (e.g., Q. robur, Q. petraea).
• Conservation & Management: The resource supports "assisted migration" strategies by enabling precise identification of adaptive alleles that could be transferred to northern populations to enhance their resilience.
• Data Availability: All raw data, assemblies, and annotations are publicly available (ENA Project PRJEB110096), serving as a foundational dataset for an ongoing European white oak pangenome project.