Two evolutionary histories in one nucleus: genome remodeling and allelic regulation underlying heterosis in hybrid oil palm

Technical Summary: Two evolutionary histories in one nucleus: genome remodeling and allelic regulation underlying heterosis in hybrid oil palm

Problem Statement
Oil palm (Elaeis) is the world's leading source of vegetable oil, primarily derived from the African oil palm (Elaeis guineensis). While interspecific hybrids between E. guineensis and the American oil palm (E. oleifera) exhibit significant heterosis (hybrid vigor), combining high yield with improved oil quality and disease resistance, the genetic and regulatory mechanisms driving this phenomenon remain unclear. Existing genomic resources rely on separate reference genomes for each parent, which obscures haplotype-specific structural variants and biases allele-specific analyses. Consequently, the precise interplay between parental genome divergence, structural remodeling, and allelic regulation within a single hybrid nucleus has not been resolved.

Methodology
The authors generated a haplotype-resolved, telomere-to-telomere (T2T) reference genome for the interspecific F1 hybrid oil palm 'Reyou 40'.

• Sequencing and Assembly: The study utilized PacBio HiFi and Oxford Nanopore ultra-long reads, integrated with Hi-C chromatin conformation data. The assembly was phased using hifiasm, resulting in two distinct haplotypes: HapG (derived from E. guineensis) and HapO (derived from E. oleifera).
• Genomic Characterization: The authors performed comparative genomics across the palm family (including Cocos nucifera, Areca catechu, etc.) and evolutionary analyses to assess gene family expansion/contraction, whole-genome duplication (WGD) events, and synteny.
• Structural and Repeat Analysis: They analyzed presence-absence variation (PAV), transposable element (TE) dynamics (specifically LTR-retrotransposons), and centromere architecture using CENH3 ChIP-seq.
• Introgression Detection: SubPhaser was employed to identify ancient introgressed regions in HapG by detecting enrichment of E. oleifera-specific k-mers.
• Transcriptomics and 3D Architecture: Hi-C maps and RNA-seq data from five fruit developmental stages were integrated to analyze three-dimensional chromatin organization (TADs and A/B compartments) and allele-specific expression (ASE) patterns.

Key Results

1. Haplotype-Resolved Assembly: The study produced a complete 1.73 Gb T2T assembly for HapG (16 gap-free chromosomes) and a near-T2T 1.84 Gb assembly for HapO (16 chromosomes with 17 gaps). Despite 91.56% genome-wide sequence identity, the haplotypes exhibit distinct repeat landscapes and centromere architectures (e.g., HapG chromosome 10 centromere lacks canonical repeats, instead containing a mosaic of TEs).
2. Complementary Functional Specialization: The two haplotypes show substantial divergence in gene content. HapG contains specific clusters of genes related to lipid metabolism (e.g., GDSL lipases), while HapO is enriched for disease-resistance genes (NB-ARC/LRR types). These PAV genes are largely associated with transposable elements.
3. Evolutionary History and Lipid Expansion: Comparative genomics revealed that while the palm family shares ancient WGD events, the oil palm lineages underwent asymmetric gene-family turnover. E. guineensis (HapG) expanded genes related to fatty-acid synthesis and oil-body formation, whereas E. oleifera (HapO) expanded genes involved in lipid modification and stress responses.
4. Ancient Introgression and Remodeling: Six ancient introgressed regions (~64 Mb) were identified in HapG, originating from E. oleifera. These regions, estimated to be ~4.0 Ma old, were subsequently remodeled by TE proliferation and tandem duplications. Notably, these regions show an enrichment of resistance genes, including an expanded RPM1 tandem cluster on HapG chromosome 2, suggesting that ancient gene flow broadened the defense repertoire of the African lineage.
5. Allelic Expression and Dosage Regulation: Transcriptomic analysis of 26,034 allelic gene pairs revealed that 82.2% maintain balanced expression across developmental stages. The remaining pairs show dynamic or stable biases. HapG-biased genes are enriched for defense and lipid biosynthesis, while HapO-biased genes are linked to hormone-mediated stress and lipid transport.
6. Tandem Duplication Effects: Tandemly duplicated genes exhibited significantly lower expression per copy compared to non-tandem genes and showed reduced allelic bias, suggesting a dosage buffering mechanism where copy-number expansion does not linearly increase transcriptional output.

Significance and Claims
The paper claims to provide the first haplotype-resolved T2T reference for an interspecific oil palm hybrid, offering a direct view of how two divergent evolutionary histories coexist within a single nucleus. The authors posit that heterosis in 'Reyou 40' is not driven by the uniform activation of one parental genome but by the integration of complementary gene repertoires (lipid production from HapG and stress resistance/metabolic flexibility from HapO) regulated by a balance of allelic expression.

The study highlights that:

• Structural Divergence vs. Functional Conservation: While the genomes show extensive structural divergence driven by TE dynamics and PAVs, essential reproductive modules remain conserved, facilitating successful hybridization.
• Regulatory Stability: The hybrid maintains transcriptional homeostasis through broad allelic balance, with locus-specific biases accommodating parental functional differences.
• Breeding Applications: These haplotype-resolved genomic resources and the identified regulatory mechanisms (e.g., dosage buffering in tandem loci, ancient introgression of resistance genes) provide vital targets for molecular breeding to develop high-yielding and resilient oil palm cultivars.

The authors conclude that their multi-omics framework links parental genome divergence, evolutionary remodeling, and regulatory balance, establishing a genomic basis for understanding heterosis in oil palm.