基因祖先揭示了真核生物发生过程中的多样化微生物关联

Gene ancestries reveal diverse microbial associations during eukaryogenesis

作者信息Moisès Bernabeu, Saioa Manzano-Morales, Marina Marcet-Houben, Toni Gabaldón
PMID42271057
期刊Nature
发布时间2026-08
DOI10.1038/s41586-026-10639-9

实验完整度

研究基于计算分析,包括多数据集系统发育重建、负对照和功能推断,但缺乏体外或体内实验验证。

主要模型

最后真核共同祖先(LECA)蛋白质组

重点核对

三个部分重叠的eTOLDB数据集(eTOLDBA、eTOLDBB、eTOLDBC) LECA标准(至少5个物种、3或5个超群、两个主要茎) 系统发育分配阈值(UFBS>75%,分类学bootstrap>95%,供体确定性>95%) 分支长度比率方法用于推断获取时间

摘要

真核生物的起源仍然是生物学中的一个核心谜团。持续的争论一致认为α变形菌与Asgard古菌之间的共生关系起着关键作用。然而,其他潜在细菌伙伴的性质、时间和贡献以及病毒相互作用的角色仍然存在争议。为了解决这些问题,我们使用了先进的系统发育基因组学方法和涵盖细胞生命和病毒已知多样性的全面数据集。我们的分析提供了对最后真核共同祖先(LECA)蛋白质组的修订重建,并追踪了每个蛋白质家族的系统发育起源。我们发现了来自多样化细菌供体的多波水平基因转移的有力证据,其中一些可能先于线粒体内共生。我们推断出了主要供体的可能特征及其对LECA的功能贡献。我们的发现支持水平基因转移对LECA之前祖先蛋白质组的塑造作用,并暗示核质病毒门(Nucleocytoviricota)病毒起到了促进作用。总之,我们的结果表明,古代真核生物可能在复杂的微生物生态系统中通过一系列多样化的关联而起源,这些关联留下了水平转移基因的足迹。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
LECA中的蛋白质家族是否能追踪到除α变形菌和Asgard古菌之外的其他祖先,并且这些贡献是否与病毒有关?
核心机制
多波水平基因转移,主要供体包括浮霉菌门、粘细菌门和核质病毒门,可能塑造了LECA蛋白质组。
主要证据
基于三个数据集重建的LECA蛋白质组,系统发育分析显示约4.5%的获得为病毒介导,其中74%来自核质病毒门;负对照和分支长度分布支持非α变形菌来源。
研究意义
支持真核生物发生是多步骤过程,涉及与多种微生物的连续相互作用和基因转移,病毒作为介导者,挑战简单的二元共生模型。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

LECA蛋白质组重建

重建LECA的蛋白质组,以识别可能存在于祖先中的基因家族。

构建了三个100蛋白质组数据集(eTOLDB),聚类为直系同源群(OG),根据物种分布定义LECA-OG,过滤污染后,通过profile搜索和系统发育重建细化为mLECA-OGs,并注释功能。

2

系统发育起源推断

确定每个基因家族的祖先来源,区分创新、获取和未知。

对每个mLECA-OG重建最大似然树,使用严格阈值分配姐妹群,识别直接获取和病毒介导获取,并应用负对照验证非α变形菌来源。

3

获取时间推断

估计不同供体基因获取的相对时间,以识别转移波。

使用分支长度比率方法计算归一化茎长,通过贝叶斯框架推断每个供体的获取模式,并比较分布峰值。

4

功能与代谢推断

推断每个供体在LECA中的功能贡献和代谢特征。

基于KEGG注释和COG分类,计算每个供体对LECA代谢模块的贡献,并重建供体基因组中的代谢途径。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

序列比对和聚类
系统发育重建
代谢重建

产品清单

实验环节名称品牌货号
BUSCO--v5.5
seqkit--stats v2.8.2
SEG----
MMseqs2--v2
CheckM----
diamond--v2.1.9
PhylomeDB流水线----
KofamScan--v1.3.0

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
eTOLDB数据集构建
物种选择、序列过滤标准(BUSCO完整性、低复杂度区域比例)
阅读提示:Methods: eTOLDB construction
LECA-OG重建
LECA标准(至少5个物种、3或5个超群、两个茎)、聚类参数(MMseqs2阈值、OrthoFinder参数)
阅读提示:Methods: LECA-OG reconstruction
系统发育分配
阈值(UFBS>75%、分类学bootstrap>95%、供体确定性>95%)
阅读提示:Methods: 'Stress test' for identification of main non-eukaryotic contributors
获取时间推断
分支长度比率方法、gamma分布拟合参数
阅读提示:Methods: Inferring the acquisition waves