IF50.5｜基因特异性选择性清除在人类肠道菌群中普遍存在

新统计量iLDS揭示数百个跨宿主传播的适应性片段，饮食与工业化是核心选择压力。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

人类肠道微生物组由高度多样的物种组成，在宿主内部和宿主之间持续进化。已有研究表明，肠道共生菌在健康成人体内可发生快速进化，新突变在数天至数月内即可达到高频率，即使没有抗生素等明显扰动。然而，对于适应性变异如何在宿主间长期传播，以及同源重组在其中的作用，仍知之甚少。肠道微生物组是水平基因转移（HGT）的热点，允许适应性等位基因轻松重组到新的遗传背景上。此前，HGT已被证明在抗生素抗性基因等跨物种传播中发挥关键作用，但HGT在多大程度上促进同一物种菌株间适应性等位基因的传播，尤其是通过同源重组，尚不清楚。

科学问题

卡在哪里： 基于连锁不平衡（LD）的选择扫描统计量在细菌中应用有限，部分原因是肠道共生菌中重组的普遍性和动态性刚刚开始被揭示，且LD统计量可能被人口收缩等非选择性进化力量混淆。此外，如何区分选择性清除与其他导致LD升高的因素，缺乏有效方法。

本文要回答： 开发一种新的选择扫描统计量，能够检测通过迁移和水平基因转移在宿主微生物组间传播的适应性等位基因的清除，并应用该统计量揭示人类肠道共生菌中基因特异性选择性清除的普遍性及其驱动因素。

技术路线

作者主要用了：iLDS统计量、SLiM模拟、MIDAS准相位、UHGG基因组集、Jaccard指数。

作者首先通过正向群体遗传模拟验证假设：选择性清除可导致常见非同义变异间的LD高于同义变异，且该信号在纯化选择与正选择共同作用下出现。随后，他们开发了整合LD得分（iLDS），该统计量结合了局部LD升高和非同义/同义LD差异两个成分，以区分正选择与其他非选择性因素。iLDS在模拟数据和已知存在重组介导选择性清除的艰难梭菌中验证后，被应用于来自全球24个人群的约30种常见肠道共生菌的宏基因组数据，以鉴定选择性清除位点并分析其功能富集和人群分布。

图 1　Fig. 1。图1展示了基因特异性选择性清除的机制和统计量原理。a部分示意携带适应性变异的基因组片段在不同菌株间转移，黄色区域代表重组到不同菌株上的适应性等位基因片段。b部分显示在中性条件下，非同义和同义变异的LD衰减曲线相似，AUC(rN2−rS2)接近零。c部分说明统计量通过比较两条曲线下面积差来量化LD升高。d部分显示在正选择强度超过纯化选择且纯化选择超过漂变时，AUC(rN2−rS2)显著大于零，表明该信号是选择性清除的特异性标志。

核心亮点

亮点 1｜iLDS有效识别已知选择性清除

在模拟中，iLDS对正在进行和近期的强选择性清除具有高检测效力，且在中性条件下假阳性率低于1%。在艰难梭菌的135个分离株中，iLDS识别出6个推定的选择性清除，其中两个与已知的毒力因子tcdB和S层盒重叠，证实了其敏感性。此外，iLDS在幽门螺杆菌和黑腹果蝇中也成功恢复了阳性对照，表明该统计量适用于任何经历重组的物种。

图 2　Fig. 2。图2展示了R. bromii和P. copri中常见变异（MAF≥0.2）的LD衰减模式。两种物种中，非同义变异的LD均显著高于同义变异，且AUC(rN2−rS2)显著大于零。这一基因组范围的模式与纯化选择和正选择共同作用于非同义位点一致，为后续使用iLDS识别特定清除位点提供了基础。

亮点 2｜肠道共生菌中广泛存在选择性清除

对32种肠道菌群物种应用iLDS，共鉴定出155个独特的清除，每个物种中位数为4个。这些清除跨越447个基因，功能多样，但显著富集于碳水化合物转运和代谢基因（COG类别G，调整后单侧P<5×10^-7），特别是糖苷水解酶（EC 3.2.1，调整后单侧P=0.02）。其中，susC/susD淀粉利用系统基因在5个物种中被反复检测到，与先前在宿主内短期进化中的发现一致。

图 3　Fig. 3。图3展示了iLDS在艰难梭菌和R. bromii中的应用。a部分显示艰难梭菌基因组中iLDS的扫描结果，橙色点表示显著窗口，蓝色条标记预测的毒力因子位置。iLDS在tcdB和S层盒等已知毒力因子处出现峰值，验证了其有效性。b和c部分分别显示E. siraeum和R. bromii中mdxEF位点的iLDS信号，表明该区域存在选择性清除。d部分总结了32个物种中每个物种检测到的清除数量，中位数为4个。

亮点 3｜mdxEF位点适应工业化饮食

在Eubacterium siraeum和Ruminococcus bromii中，iLDS检测到mdxEF基因（编码麦芽糊精ABC转运体）的选择性清除。麦芽糊精是超加工食品中常用的合成淀粉衍生物。在R. bromii中，该位点在所有14个工业化人群（欧洲、亚洲、北美）中均被检测到，但在两个非工业化人群（斐济和萨尔瓦多）中缺失，表明该物种可能正在适应工业化饮食中的新淀粉来源。

图 4　Fig. 4。图4展示了工业化与非工业化人群间选择性清除的差异。a部分显示分析人群的地理分布。b部分显示工业化人群之间共享清除的Jaccard指数显著高于与非工业化人群的共享。c部分显示R. bromii中mdxEF位点在所有工业化人群中均被检测到，但在非工业化人群中缺失，突出了该位点对工业化饮食的适应。

亮点 4｜工业化与非工业化人群选择靶点分化

在来自24个人群的16个物种中，共鉴定出309个独特的清除，其中201个为单一人群特有，提示局部适应。工业化人群之间共享清除的频率（Jaccard指数J=0.21）显著高于与非工业化人群的共享（J=0.11，P<10^-4），而非工业化人群之间的共享（J=0.18）与工业化人群无显著差异。此外，32个清除在一种生活方式的人群中频率≥50%而在另一种中缺失，其中24个为工业化特有，8个为非工业化特有。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： iLDS统计量、SLiM模拟、MIDAS准相位、UHGG基因组集、Jaccard指数
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事微生物组进化、群体遗传学或生物信息学方法开发的研究人员。

带走一句

iLDS为细菌群体遗传学提供了检测重组介导选择性清除的新工具，可应用于宏基因组数据，揭示适应性进化的驱动因素。