henhaobio
Sign in

NatureEvolutionary BiologyPeer-reviewed journal article

Dominant contribution of Asgard archaea to eukaryogenesis

Victor TobiassoniD, Jacob LuoiD, Yuri I. WolfiD, Eugene V. KooniniD

Abstract

The origin of eukaryotes is one of the key problems in evolutionary biology 1,2. The demonstration that the Last Eukaryotic Common Ancestor (LECA) already contained the mitochondrion, an endosymbiotic organelle derived from an alphaproteobacterium, and the discovery of Asgard archaea, the closest archaeal relatives of eukaryotes 3-7, inform and constrain evolutionary scenarios of eukaryogenesis 8. We undertook a comprehensive analysis of the origins of core eukaryotic genes tracing to the LECA within a rigorous statistical framework centered around evolutionary hypotheses testing using constrained phylogenetic trees. The results reveal dominant contributions of Asgard archaea to the origin of most of the conserved eukaryotic functional systems and pathways. A limited contribution from Alphaproteobacteria was identified, primarily relating to the energy transformation systems and Fe-S cluster biogenesis, whereas ancestry from other bacterial phyla was scattered across the eukaryotic functional landscape, without clear consistent trends. These findings suggest a model of eukaryogenesis in which key features of eukaryotic cell organization evolved in the Asgard lineage leading to the LECA, followed by the capture of the Alphaproteobacterial endosymbiont, and augmented by numerous but sporadic horizontal acquisition of genes from other bacteria both before and after endosymbiosis.

In-depth summary

Claude Opus 5.5 · xhigh effort · 7 Nov 2026 · from the full text of version 3 (14 Jul 2025) · CC0 1.0

Constrained-tree hypothesis testing on 4,290 core eukaryotic protein clusters attributes 50% of likelihood weight to Asgard archaea; Alphaproteobacterial ancestry is largely confined to oxidative phosphorylation and iron-sulfur cluster biogenesis.

Earlier phylogenomic studies found that bacterial genes outnumbered archaeal ones among the core eukaryotic genes traceable to the last eukaryotic common ancestor (LECA), with archaeal genes mostly confined to information processing. Using a larger genome collection and likelihood tests of competing sister clades, this work extends Asgard ancestry to most eukaryotic functional systems, including much of core metabolism, and finds no sign of another symbiotic partner besides the mitochondrion.

Key results

  1. ShownAsgard archaea are the largest inferred source of conserved eukaryotic protein families likely present in the LECA. The result holds across clustering schemes and taxonomies.

    Across 13,500 Eukaryotic-Prokaryotic Orthologous Clusters (EPOCs) covering 3950 KEGG Orthology groups (KOGs), Asgard had the highest average expected likelihood weight (aELW). In the core set of 4290 EPOCs (2100 KOGs) left after removing EPOCs with maximum ELW below 0.4, Asgard held 50% of all likelihood weight. The Asgard lead persisted with tree-annealed clusters and under GTDB taxonomy. · Phylogenomic analysis and global evolutionary…; Cluster annealing; Taxonomy remapping and EPOC construction under…; Figure legends

  2. ShownAsgard ancestry covers almost all core information processing and reaches many other systems, including N-linked glycosylation, glycosylphosphatidylinositol (GPI) anchor synthesis, nucleocytoplasmic transport, sphingolipid metabolism and the mevalonate pathway; fatty acid pathways look more mosaic.

    Strong Asgard aELW for ribosome, RNA and DNA polymerases, Ras-like GTPases, ubiquitin-proteasome degradation, DNA repair (nucleotide excision, mismatch), recombination and aminoacyl-tRNA biosynthesis; also for the oligosaccharyltransferase core, 7-subunit GPI-GnT complex and GPI transamidase. Fatty acid synthesis and decay gave widely spread aELW, with Actinomycetota and Alphaproteobacteria signals. · Phylogenomic analysis and global evolutionary…; Broad Asgard contributions to eukaryogenesis

  3. ShownAlphaproteobacterial ancestry is clear but narrow, essentially limited to mitochondrial functions: oxidative phosphorylation, ubiquinone synthesis, mitochondrial translation and 4Fe-4S cluster assembly.

    Alphaproteobacteria were associated with oxidative phosphorylation complexes, ubiquinone synthesis, the 4Fe-4S assembly platforms, IBA57 and ferredoxin-1 and 2. In contrast, the 2Fe-2S scaffold ISCU was mosaic with a weaker Asgard signal, while NFS1, cysteine and methionine biosynthesis and the cytosolic CIA components CIA1 and CIA2B were Asgard-associated. · Limited and specific Alphaproteobacterial…; Phylogenomic analysis and global evolutionary…; Discussion

  4. ShownAbout a third of core families point to bacteria other than Alphaproteobacteria, but these signals are scattered across functions; apart from a narrow Myxococcota nucleotide-metabolism signal, no pathway traces to one bacterial taxon.

    600 of 2,100 KOGs and 1,580 of 4,300 EPOCs had their top association with neither Asgard nor Alphaproteobacteria. Only Alphaproteobacteria (oxidative phosphorylation) and Myxococcota (nicotinate and nucleotide pathways, mostly nucleotidases) beat Asgard aELW in a KEGG pathway with over 20 EPOCs. Of 130 strictly filtered EPOCs (ELW > 0.7), 97 were Asgard and 33 diverse bacteria, with no functional trend. · Absence of functionally consistent contributions…; Figure legends

  5. SuggestedNormalized stem lengths behave poorly as a clock for gene acquisition: shorter Alphaproteobacterial stems (the earlier mitochondria-late signal) appear only among long stems, and co-acquired oxidative phosphorylation genes vary widely.

    Of 5,850 stems (sharp peak near 0.05), Alphaproteobacterial stems below 0.3 (70% of all stems) were longer than Asgard stems, but those above 0.35 were shorter on average; the pattern held in 20 technical replicates of 175 sampled stems. Oxidative phosphorylation stems spanned more than two orders of magnitude, and information processing genes had shorter than average stems. · Ancestral stem lengths and the amount of…; Figure legends

Methods, data and reproducibility

Soft-core eukaryotic and prokaryotic pangenomes were clustered (mmseqs2) into HMM profiles, and HHblits profile-profile hits defined Eukaryotic-Prokaryotic Orthologous Clusters (EPOCs). Per EPOC, IQtree2 constraint trees forcing each of up to the 12 nearest prokaryotic clades as sister to the eukaryotic clade gave expected likelihood weights (ELW), averaged over KEGG-annotated families.

  • Proteins from complete prokaryotic genomes in NCBI GenBank, November 2023 (prok2311) · bacteria and archaea · n = 47,545 genomes; 75 million proteins
  • Asgard archaeal genome assemblies (metagenomes) added to the prokaryotic database · Asgard archaea · n = 63 assemblies; 441150 proteins
  • Eukaryotic proteomes: 72 NCBI GenBank genomes plus EukProt v3 (Euk72Ep), screened for prokaryotic contaminants · eukaryotes · n = 993 species in EukProt v3; 30.3 million sequences
Methods & tools
mmseqs2, HH-suite (HHblits 3.3.0), muscle5 5.1, IQtree2 2.3.5, FAMSA 2.0.1, FastTree 2.1.1, Prodigal v2.6.3, ete3 3.1.3, UMAP 0.5.3, HDBSCAN (scikit-learn 1.3.0), KEGG release 110, GTDB version 220
Software & environment
Custom Python 3.9 scripts; stated versions: Prodigal 2.6.3, FAMSA 2.0.1, FastTree 2.1.1, muscle5 5.1, HHblits 3.3.0, IQtree2 2.3.5, ete3 3.1.3, UMAP 0.5.3, scikit-learn 1.3.0, scipy 1.11.1, pandas 2.0.3, Altair 4.2.2.
Reproducibility gaps
Hardware, run time and memory are not reported, and no container or dependency lock file is mentioned; the UniProt ID-mapping link is printed in a malformed form.
Limitations (6)
  • Prokaryotic genome sampling is biased and limited; nearly all Asgard genomes are incomplete metagenomic assemblies, making the size of the Asgard contribution uncertain. · stated by the authors
  • Sequenced Alphaproteobacteria skew toward symbionts and parasites with reduced genomes, and the mitochondrial ancestor appears to lie outside known Alphaproteobacterial diversity. · stated by the authors
  • Other bacterial phyla are even more under-sampled, so their contributions may be underestimated; for example, only part of a larger Myxococcota signal may have been captured. · stated by the authors
  • The pangenome filter against late HGT could miss contributions from narrow or poorly sampled clades and may not offset extensive ongoing horizontal transfer among eukaryotes. · stated by the authors
  • The LECA gene set was not explicitly reconstructed and annotation relies on KEGG; results would need revisiting if KEGG poorly represents the LECA gene complement. Membrane lipid biosynthesis ancestry remains unresolved. · stated by the authors
  • The headline 50% is a share of averaged ELW, which the authors treat as proportional to sister-clade probability, and it depends on chosen cutoffs (maximum ELW 0.4; ELW 0.99 or above excluded except for Alphaproteobacterial oxidative phosphorylation). · noted in this summary

Read this if you work on eukaryogenesis or Asgard archaea, or want a reusable pipeline and released EPOC data for testing which prokaryotic clade is sister to eukaryotic protein families.

Citations
40
FWCI
63.54
Field percentile
top 1%

Reproducibility ledger

Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.

Sign in to file a reproducibility report.