Lateral Gene Transfer from the Dead

In phylogenetic studies, the evolution of molecular sequences is assumed to have taken place along the phylogeny traced by the ancestors of extant species. In the presence of lateral gene transfer (LGT), however, this may not be the case, because the species lineage from which a gene was transferred may have gone extinct or not have been sampled. Because it is not feasible to specify or reconstruct the complete phylogeny of all species, we must describe the evolution of genes outside the represented phylogeny by modelling the speciation dynamics that gave rise to the complete phylogeny. We demonstrate that if the number of sampled species is small compared to the total number of existing species, the overwhelming majority of gene transfers involve speciation to, and evolution along extinct or unsampled lineages. We show that the evolution of genes along extinct or unsampled lineages can to good approximation be treated as those of independently evolving lineages described by a few global parameters. Using this result, we derive an algorithm to calculate the probability of a gene tree and recover the maximum likelihood reconciliation given the phylogeny of the sampled species. Examining 473 near universal gene families from 36 cyanobacteria, we find that nearly a third of transfer events – 28% – appear to have topological signatures of evolution along extinct species, but only approximately 6% of transfers trace their ancestry to before the common ancestor of the sampled cyanobacteria.

💡 Research Summary

The paper challenges the conventional assumption in phylogenetics that molecular sequences evolve solely along the lineage of extant species. In microbial evolution, lateral gene transfer (LGT) frequently moves genes between lineages, and the donor lineage may be extinct or simply unsampled. Because reconstructing the full tree of all species is impossible, the authors propose a probabilistic framework that models the speciation‑extinction dynamics of the entire (known and unknown) species pool. By treating the total number of species as generated by a birth‑death process, they derive analytical expressions for the probability that a gene transfer involves a lineage that is not represented in the sampled phylogeny. Their key theoretical result is that when the proportion of sampled species is small relative to the total diversity—a common situation in microbial studies—the overwhelming majority of LGT events are expected to pass through extinct or unsampled lineages. Moreover, the evolutionary history of genes along these hidden lineages can be approximated as independent lineages described by a handful of global parameters: the speciation rate, extinction rate, transfer rate, and sampling fraction.

Building on this insight, the authors develop a dynamic‑programming algorithm that computes the likelihood of any given gene tree under the model and identifies the maximum‑likelihood reconciliation between the gene tree and the sampled species tree. The algorithm treats each hidden lineage as an independent “ghost” branch, allowing the calculation to remain tractable despite the combinatorial explosion that would arise from explicitly enumerating all possible extinct taxa.

To validate the method, they applied it to 473 near‑universal gene families from 36 cyanobacterial genomes. The analysis revealed that approximately 28 % of inferred transfer events bear the topological signature of having passed through extinct or unsampled lineages. In contrast, only about 6 % of transfers can be traced back to a donor that existed before the most recent common ancestor of the sampled cyanobacteria. This empirical finding aligns with the theoretical prediction that most LGT events are recent and mediated by lineages that are not present in the sampled data set.

The study has several important implications. First, it quantifies the hidden contribution of extinct or unsampled taxa to gene flow, highlighting that ignoring these “ghost” lineages can lead to substantial underestimation of LGT frequency and misinterpretation of evolutionary histories. Second, the model’s reliance on a small set of global parameters makes it computationally feasible for large‑scale metagenomic analyses where only a fraction of the community is sampled. Third, the predominance of recent, hidden‑lineage transfers suggests that microbial communities maintain a dynamic gene pool, facilitating rapid adaptation to environmental changes. Finally, the framework can be extended to other mobile genetic elements such as plasmids, phages, and transposons, offering a versatile tool for evolutionary biologists, ecologists, and public‑health researchers interested in the spread of antimicrobial resistance or virulence factors across microbial networks.