IF50.5｜人肠道温带噬菌体的分离、工程化与生态学

只有18%的预测前噬菌体能被诱导，人类细胞裂解产物可触发诱导。

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

研究背景

人类肠道微生物组中噬菌体多样性极高，但功能研究受限于缺乏可培养的噬菌体-宿主配对。大规模宏基因组学揭示了大量病毒基因组，且高达90%的肠道细菌被预测为溶原菌。然而，这些前噬菌体在何种条件下重新进入裂解周期仍不清楚。前噬菌体的诱导受宿主和噬菌体特异性信号调控，在肠道环境中尤其复杂。培养组学的进步使得许多先前“不可培养”的肠道细菌及其噬菌体得以分离，为实验验证提供了基础。

科学问题

卡在哪里： 尽管计算预测表明肠道细菌中前噬菌体普遍存在，但缺乏系统的实验验证，不清楚有多少预测的前噬菌体真正具有诱导能力，以及哪些环境因素（特别是人类宿主相关因子）能触发诱导。

本文要回答： 系统表征来自人肠道细菌分离株的可诱导前噬菌体，评估诱导效率，并探究人类细胞产物是否作为诱导信号。

技术路线

作者主要用了：细菌纯培养诱导、合成菌群共培养、Caco2细胞共培养、宏基因组测序、qPCR、CRISPR基因敲除。

作者首先对252株人肠道细菌分离株进行前噬菌体计算预测，然后用8种标准诱导条件（如丝裂霉素C、过氧化氢、甜菊糖、饥饿等）处理纯培养物，通过DNA测序检测诱导的前噬菌体。为模拟肠道环境，构建了78株菌的合成菌群，并与Caco2细胞单层共培养，比较有无人类细胞时的诱导情况。进一步用32株菌的纯培养直接暴露于Caco2细胞单层、细胞裂解液或培养基，以区分人类细胞因子的作用。最后通过基因敲除和qPCR验证关键基因功能。

图 1　Fig. 1。图1展示了整体实验设计和诱导结果。a部分示意了三种方法：纯培养诱导、合成菌群与Caco2共培养、纯培养暴露于人类细胞产物。b部分显示252株菌的系统发育树，按门着色。c-f部分显示合成菌群实验中，Caco2共培养诱导的噬菌体物种数（57）显著高于无细胞共培养（29），且群落组成发生改变。g部分显示32株菌纯培养暴露于Caco2细胞产物时，29个前噬菌体被诱导，其中9个为新发现。h-i部分总结纯培养中仅18%的预测前噬菌体被诱导，Bacteroidota诱导率最高。

核心亮点

亮点 1｜仅少数预测前噬菌体可被诱导

在252株细菌中，94%被预测含有高质量前噬菌体区域，但在10种诱导条件下，仅32%的菌株被诱导，18%的高质量前噬菌体预测（134/736）或24%的高质量前噬菌体物种（68/274）对应实验可诱导的前噬菌体。Bacteroidota中诱导率最高（27%的预测），而Pseudomonadota仅12%。合成菌群共培养中检测到36%的前噬菌体物种被诱导，与近期宏基因组研究估计的8-36%一致。这表明肠道中大部分预测前噬菌体可能很少进入裂解周期。

图 2　Fig. 2。图2展示了诱导噬菌体的分类和流行率。a部分为基因共享网络，显示68个可诱导前噬菌体物种与参考数据库的关系，其中Hankyvirus属与Bacteroides噬菌体聚类。b部分显示这些噬菌体在1241个肠道病毒组中的检出率和丰度，LoVEphage最为常见（8%病毒组），而最丰富的Crassvirales基因组在19%病毒组中检出。

亮点 2｜人类细胞裂解产物触发前噬菌体诱导

在合成菌群与Caco2细胞共培养中，35%的噬菌体物种被诱导，显著高于无细胞共培养的17%。进一步用32株菌的纯培养暴露于Caco2细胞单层、细胞裂解液或DMEM培养基，发现29个前噬菌体被诱导，其中25个在细胞裂解液条件下，14个在单层和培养基中。重要的是，9个前噬菌体之前未在标准诱导剂中被检测到，表明人类细胞裂解产物可作为诱导信号。这为炎症性肠病中噬菌体扩增与细胞死亡相关的观察提供了机制联系。

图 3　Fig. 3。图3展示了前噬菌体基因组特征和失活机制。a部分显示诱导、高质量预测和低质量预测前噬菌体的长度分布，呈双峰。b部分显示诱导前噬菌体富集结构基因，而非诱导预测富集辅助基因和未知功能基因。c部分显示前噬菌体dN/dS高于宿主，但诱导与非诱导对之间无差异。d部分显示143个基因dN/dS>1，多数与DGRs相关。e-f部分显示ΦPomma转座酶基因敲除后，在甜菊糖和过氧化氢诱导下qPCR和测序均无诱导信号。

亮点 3｜DGRs在肠道前噬菌体中富集并可能介导宿主切换

在68个可诱导前噬菌体物种中，19%编码多样性生成逆转录元件（DGRs），且主要存在于Bacteroidota噬菌体中（43%）。这些DGRs靶向已知或预测的尾蛋白。同时发现8个Bacteroidota前噬菌体物种在不同细菌物种间活跃复制，其中3个跨属复制，提示DGRs通过尾蛋白多样化扩展宿主范围。此外，4个前噬菌体物种含有第二个可变区，靶向远离逆转录酶盒的基因，邻近防御基因如DNA甲基转移酶，表明DGRs参与噬菌体-宿主军备竞赛。

图 4　Fig. 4。图4展示了诱导条件、多溶原和宿主遗传的影响。a部分热图显示每个分离株在不同条件下诱导的前噬菌体数量，丝裂霉素C诱导最多。b部分显示Bacteroidota中多溶原比例最高。c部分显示前噬菌体数量与诱导条件数正相关。d部分显示5株B. caccae中ΦWilby和ΦPomma的差异诱导。e部分显示ΦPomma在不同菌株中整合位点不同，而ΦWilby位点保守。

亮点 4｜整合/切除基因突变导致前噬菌体失活

比较诱导与非诱导的高相似前噬菌体对（95% ANI），发现非诱导前噬菌体中整合和切除相关基因的dN/dS显著升高（P=0.002），表明这些基因的非功能性突变是前噬菌体失活的途径。通过敲除可诱导前噬菌体Pomma的DNA转座蛋白基因（ΦPomma Δtran），在甜菊糖和过氧化氢诱导下，qPCR显示野生型比突变体分别高3.5和2.6个对数级，测序也证实突变体无诱导，证明该基因对诱导至关重要。

亮点 5｜多溶原现象普遍且影响诱导协调性

在具有可诱导前噬菌体的分离株中，52%为多溶原菌。Bacteroidota中多溶原比例最高（68%），显著高于其他门（29%，P=0.002）。多溶原菌中的前噬菌体平均在2.7种条件下被诱导，而单溶原菌为2.1种（P=0.03），且前噬菌体数量与诱导条件数正相关（τ=0.22, P=0.002），表明多溶原可能促进协同诱导并降低溶原稳定性。

亮点 6｜宿主遗传背景和整合位点影响差异诱导

在5株高度相似（99% ANI）的Bacteroidota caccae分离株中，均携带ΦWilby和ΦPomma两个前噬菌体。qPCR显示在标准培养基中ΦWilby优先诱导（P=0.006），但在过氧化氢处理中无差异。不同分离株间ΦWilby/ΦPomma比值差异显著（ANOVA P=0.012和P=0.0008）。长读长测序发现ΦWilby整合在同一tRNA位点，而转座型ΦPomma在5株菌中位于4个不同基因组位置，提示整合位点差异可能是差异诱导的原因。

生信可带走

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

• 方法栈： 细菌纯培养诱导、合成菌群共培养、Caco2细胞共培养、宏基因组测序、qPCR、CRISPR基因敲除
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事肠道噬菌体、微生物组-宿主互作或培养组学研究的科研人员。

带走一句

计算预测的前噬菌体大多不可诱导，实验验证必不可少；人类细胞裂解产物是重要的诱导信号，值得在机制研究中关注。