This is a Discovery ($\Psi_{\text{Discovery}}$) paper. While it introduces a refined implementation of molecular clock calibration (“cross-bracing”), the primary contribution is the biological inference regarding LUCA’s age, genome size, and metabolic nature. The computational methods serve to characterize a specific biological entity.

What is the motivation?

Understanding the Last Universal Common Ancestor (LUCA) is critical for reconstructing the early evolution of life, yet consensus has been elusive due to disparate data and methods.

• Age Conflicts: Estimates vary widely depending on fossil interpretation and molecular clock calibrations, particularly regarding the “Late Heavy Bombardment” (LHB) constraints.
• Physiological Uncertainty: Debates persist over whether LUCA was a simple “progenote” dependent on geochemistry or a complex prokaryote-grade organism.
• Environmental Context: LUCA is often modeled in isolation, ignoring the ecological interactions that would have shaped its survival and impact on the early Earth system.

What is the novelty here?

The study integrates three advanced computational approaches to provide a holistic reconstruction of LUCA:

• Cross-Braced Dating: It employs a “cross-bracing” strategy in Bayesian molecular clocks, using pre-LUCA gene duplications (paralogues) to constrain the root. This allows the same fossil calibrations to be applied to mirrored nodes, significantly reducing uncertainty.
• Probabilistic Reconciliation: It uses the ALE (Amalgamated Likelihood Estimation) algorithm to reconcile 9,365 gene family trees against the species tree. This explicitly models gene transfer, duplication, and loss, allowing for a much broader reconstruction of the proteome.
• Ecosystem Modeling: The physiological reconstruction is coupled with geochemical modeling to propose that LUCA was a member of a productive, hydrogen-recycling early ecosystem.

What experiments were performed?

• Phylogenomics: Inferred a species tree from 57 single-copy marker genes across 700 diverse prokaryotic genomes (350 Archaea, 350 Bacteria) using maximum likelihood (IQ-TREE 2).
• Molecular Dating: Estimated divergence times using MCMCtree with a partitioned dataset of 5 pre-LUCA paralogue pairs (e.g., ATP synthase, EF-Tu/G). Calibrations included 13 fossil constraints and a “soft” maximum bound based on the Moon-forming impact (4.51 Ga).
• Metabolic Reconstruction: Reconciled 9,365 KEGG ortholog families against the species tree to calculate the posterior probability (PP) of each gene’s presence in LUCA. Metabolic potential was inferred from genes with high PP (typically >0.75).
• Genome Size Prediction: Trained a LOESS regression model on modern prokaryotes to predict LUCA’s genome size based on the inferred number of KEGG families.

What outcomes/conclusions?

• Age: LUCA lived approximately 4.2 Ga (95% CI: 4.09-4.33 Ga), surprisingly soon after the Moon-forming impact (~4.5 Ga).
• Complexity: LUCA was a complex, prokaryote-grade organism with a genome size of ~2.75 Mb (encoding ~2,600 proteins), comparable to modern prokaryotes.
• Physiology:
  • Metabolism: Anaerobic acetogen using a complete Wood-Ljungdahl pathway (WLP) for $CO_2$ fixation and an almost complete TCA cycle. Likely thermophilic (reverse gyrase present, PP = 0.97). The paper found no strong evidence for nitrogenase or nitrogen fixation.
  • Immunity: Possessed 19 Class 1 (Type I and Type III) CRISPR-Cas effector protein families. Cas1 and Cas2 were absent, suggesting an early immune system capable of RNA cleavage and binding but lacking the full CRISPR adaptation machinery.
• Ecology: LUCA likely inhabited one of two major habitats: (1) the deep ocean, where hydrothermal vents and serpentinization provided $H_2$ (supported by reverse gyrase presence, PP = 0.97, consistent with hyperthermophily), or (2) the ocean surface, where atmospheric $H_2$ from volcanism and metamorphism could fuel growth. A shallow hydrothermal vent or hot spring is also considered a possibility. LUCA was part of an established ecosystem whose metabolic by-products would have created niches for other metabolisms, including methanogenesis. If methanogens were also present, the $CH_4$ they produced would have been photochemically recycled to $H_2$ in the atmosphere, boosting biosphere productivity by at least an order of magnitude over abiotic $H_2$ input rates.
• Limitation: The placement of two small-genome lineages (CPR, Candidate Phyla Radiation, and DPANN) remained uncertain. The AU (approximately unbiased) test could not reject either topology (p = 0.517), meaning the data cannot discriminate between the two placements. This phylogenetic uncertainty affects inferences about the early bacterial and archaeal stem lineages.

Reproducibility Details

Data

The study relied on publicly available genomic data and specific subsets of marker genes.

Algorithms

Key computational steps involved sequence processing, tree inference, and probabilistic reconciliation.

• Alignment & Trimming: sequences aligned with MAFFT L-INS-i (v7.407) and trimmed with BMGE (v1.12, BLOSUM30 matrix, entropy 0.5).
• Tree Inference: IQ-TREE 2 (v2.1.2). Species tree: LG+C60+F+G (best-fit by BIC from concatenated 57-marker alignment). Gene family trees for ALE reconciliation (9,365 KEGG families): LG+F+G with 1,000 ultrafast bootstraps.
• Reconciliation: ALE (Amalgamated Likelihood Estimation) program ALEml_undated used to calculate gene presence probabilities, accounting for HGT, duplication, and loss.
• Genome Prediction: LOESS regression (Locally Estimated Scatterplot Smoothing) used to map KEGG family counts to total protein counts/genome size.

Models

The analysis employed sophisticated evolutionary models to handle deep time scales and heterogeneity.

• Substitution Models:
  • Species Tree: LG+C60+F+G (mixture model with 60 profiles, best-fit by BIC).
  • Gene Family Trees (for ALE reconciliation): LG+F+G with 1,000 ultrafast bootstraps.
  • Timetree inference: LG+F+G4 for approximate likelihood calculation (CODEML), as CODEML does not implement the CAT mixture model.
• Molecular Clock:
  • MCMCtree (PAML v4.10.7).
  • Relaxed clock models: GBM (Geometric Brownian Motion) and ILN (Independent Lognormal).
  • Cross-Bracing: Specifically models shared divergence times for duplicated nodes (driver and mirror nodes).

Evaluation

Validation focused on robustness across different topologies and clock models.

Hardware

• Software: PAML v4.10.7 (MCMCtree), IQ-TREE 2, ALE v0.4, HMMER v3.3.2.
• Compute: IQ-TREE runs specified usage of 4 CPUs; MCMCtree approximated likelihood calculation (approx method) to reduce computational cost.

Paper Information

Citation: Moody, E. R. R., Álvarez-Carretero, S., Mahendrarajah, T. A., et al. (2024). The nature of the last universal common ancestor and its impact on the early Earth system. Nature Ecology & Evolution, 8(9), 1654-1666. https://doi.org/10.1038/s41559-024-02461-1

Publication: Nature Ecology & Evolution, Volume 8, Number 9, 2024

@article{moodyTheNatureLast2024,
  title={The nature of the last universal common ancestor and its impact on the early Earth system},
  author={Moody, Edmund R. R. and Álvarez-Carretero, Sandra and Mahendrarajah, Tara A. and Clark, James W. and Betts, Holly C. and Dombrowski, Nina and Szánthó, Lénárd L. and Boyle, Richard A. and Daines, Stuart and Chen, Xi and Lane, Nick and Yang, Ziheng and Shields, Graham A. and Szöllősi, Gergely J. and Spang, Anja and Pisani, Davide and Williams, Tom A. and Lenton, Timothy M. and Donoghue, Philip C. J.},
  journal={Nature Ecology & Evolution},
  volume={8},
  number={9},
  pages={1654--1666},
  year={2024},
  publisher={Nature Publishing Group},
  doi={10.1038/s41559-024-02461-1}
}

Open Access: This article is published under CC BY 4.0 and is freely available at the paper URL above.

Artifacts:

This is a Systematization paper with a strong Position component.

Systematization: It synthesizes decades of molecular sequence data (specifically rRNA) to propose a “formal system of organisms” that replaces previous taxonomies.

Position: It argues that the prevailing “Prokaryote-Eukaryote” and “Five Kingdom” models are “outmoded,” “misleading,” and based on “flawed premises” regarding the organization of life.

What is the motivation?

The authors aim to align formal taxonomy with the “natural system emerging from molecular data”.

The Problem: Existing systems (Whittaker’s 5-Kingdoms) were based on morphology and nutrition, which are insufficient for microbial classification.

The Gap: The “Prokaryote” definition was negative (defined by what they lack, a nucleus), obscuring the fact that “Archaebacteria” are as distinct from “Eubacteria” as they are from Eukaryotes.

The Goal: To establish a taxonomic rank above Kingdom that recognizes the three primary evolutionary lineages.

What is the novelty here?

The core contribution is the formal proposal of the Domain as the highest taxonomic rank. Specific novel definitions include:

1. Three Domains:

   • Bacteria (formerly Eubacteria): Membrane lipids are diacyl glycerol diesters; eubacterial rRNA.
   • Archaea (formerly Archaebacteria): Membrane lipids are isoprenoid glycerol diethers/tetraethers; archaeal rRNA. The term “archaebacteria” is abandoned to emphasize their independence.
   • Eucarya (Eukaryotes): Cells with nuclei; glycerol fatty acyl diester lipids; eukaryotic rRNA.

2. Subdivision of Archaea: The domain is formally split into two kingdoms:

   • Euryarchaeota (methanogens, halophiles, thermoplasms, sulfate-reducing Archaeoglobus, and thermophiles Thermococcus and Pyrococcus).
   • Crenarchaeota (sulfur-dependent extreme thermophiles).

What experiments were performed?

This paper is a synthesis of phylogenetic analysis. It relies on:

• rRNA Sequencing: Comparison of 16S (small subunit) ribosomal RNA sequences. The paper cites over 400 known eubacterial cases of a characteristic structural feature (the 6-nucleotide side bulge at positions 500-545).
• Phylogenetic Tree Reconstruction: Analysis of branching orders and lengths based on rRNA sequence comparisons (citing Woese, 1987).
• Paralogous Gene Rooting: Determining the root of the universal tree by comparing duplicated genes (e.g., elongation factors) that diverged before the three lineages separated.

What outcomes/conclusions?

• Tripartite Division: Life divides into three monophyletic groups. The evolutionary differences among the three domains are more profound than those separating traditional kingdoms such as animals and plants.
• Archaea-Eucarya Sisterhood: The root of the tree separates Bacteria from the other two, making Archaea and Eucarya sister groups.
• Molecular Definition: Phenotypes are replaced by molecular signatures. For example, Bacteria are defined by a 6-nucleotide bulge in the 16S rRNA (positions 500-545), whereas Archaea and Eucarya share a 7-nucleotide bulge.
• “Prokaryote” as Invalid Taxon: The paper explicitly argues that “Prokaryote” is not a valid natural taxon. Because it is defined by what the organisms lack (a nucleus), it groups together two deeply divergent domains (Bacteria and Archaea) by a plesiomorphic character. The term should be abandoned in natural classification.
• Domain Replaces Kingdom: Introducing the Domain rank above Kingdom resolves the issue. A bacterium is no more related to an archaeon than either is to a eukaryote, so all three deserve equivalent top-level status.
• Formal Conclusions (adapted from paper):
  1. Life comprises three primary groupings, the Domains Bacteria, Archaea, and Eucarya.
  2. None of these is ancestral to the others; all descend from a common ancestor.
  3. The Archaea comprise two kingdoms, Euryarchaeota and Crenarchaeota.
  4. Both Bacteria and Eucarya will contain numerous kingdoms; for Eucarya, the paper anticipates preserving Plantae, Animalia, and Fungi while replacing Protista with several kingdoms.
  5. “Prokaryote” has no phylogenetic meaning and should not be used as a formal taxon.

Reception and ongoing debate: At publication, abandoning “prokaryote” was a controversial claim. Most microbiology and cell biology textbooks through the 2000s retained the term, and many introductory curricula continue to use it today. The three-domain framework has since been adopted in modern phylogenetics and comparative genomics, but the transition is not yet universal in pedagogy, and some researchers have proposed alternative deep-tree topologies (e.g., the eocyte hypothesis) that differ from Woese’s original Archaea-Eucarya sisterhood.

Reproducibility Details

Note: As a theoretical systematics paper from 1990, “reproducibility” refers to the data sources and criteria used to construct the taxonomy.

Data

The taxonomy rests on comparative analysis of Ribosomal RNA (rRNA), specifically the small subunit (16S in prokaryotes, 18S in eukaryotes).

Algorithms

The paper relies on rRNA sequence comparisons to generate the universal tree in Figure 1, using phylogenetic methods standard at the time.

• Tree Inference: Branching order/lengths taken from Microbiol. Rev. 51, 221-271 (1987).
• Rooting Strategy: The “Outgroup” method using anciently duplicated genes (paralogs) such as Elongation Factors Tu and G, which diverged prior to the Universal Ancestor.

Models

The “Model” proposed is the Three-Domain System:

1. Domain Bacteria: Rooted independently. Includes Thermotogales, Flavobacteria, Cyanobacteria, Purple bacteria, Gram-positives, Green nonsulfur.
2. Domain Archaea:
   • Kingdom Crenarchaeota: “Ancestral” phenotype (thermophily). Includes Pyrodictium, Thermoproteus.
   • Kingdom Euryarchaeota: “Broad” phenotype. Includes Methanogens, Halophiles, Thermoplasma, Archaeoglobus (sulfate-reducing), and Thermococcus and Pyrococcus (thermophilic).
3. Domain Eucarya: Includes Animals, Ciliates, Plants, Fungi, Flagellates, Microsporidia.

Evaluation

The authors validate the model by demonstrating Molecular Invariants: features present in all members of a domain but absent in others:

Paper Information

Citation: Woese, C. R., Kandler, O., & Wheelis, M. L. (1990). Towards a natural system of organisms: Proposal for the domains Archaea, Bacteria, and Eucarya. Proc. Natl. Acad. Sci. USA, 87(12), 4576-4579. https://doi.org/10.1073/pnas.87.12.4576

Publication: Proc. Natl. Acad. Sci. USA, Volume 87, Number 12, 1990

@article{woeseNaturalSystemOrganisms1990,
  title = {Towards a Natural System of Organisms: Proposal for the Domains {{Archaea}}, {{Bacteria}}, and {{Eucarya}}.},
  shorttitle = {Towards a Natural System of Organisms},
  author = {Woese, C R and Kandler, O and Wheelis, M L},
  year = {1990},
  month = jun,
  journal = {Proceedings of the National Academy of Sciences of the United States of America},
  volume = {87},
  number = {12},
  pages = {4576--4579},
  issn = {0027-8424},
  doi = {10.1073/pnas.87.12.4576}
}

Additional Resources:

• Woese’s 1977 Discovery of Archaea
• NCBI Taxonomy Browser