《Frontiers in Microbiology》：Genomic and functional analysis of stress-responsive prophages in Lactobacillus helveticus
编辑推荐：
为了解决传统硬质奶酪生产中噬菌体在自然乳清发酵剂(NWS)中的功能角色尚不明确的问题，本研究从特伦蒂诺格拉纳奶酪NWS中分离了六株温和噬菌体。通过结合基因组学、宿主范围分析和应激诱导实验，研究人员系统探究了热激、营养限制、乳酸及氧化应激等乳品加工相关胁迫对前噬菌体诱导动态的影响。结果表明，不同胁迫条件可引发噬菌体特异性、剂量依赖性的诱导响应，其中热激是最稳定的前噬菌体激活诱导因素。该研究揭示了一种从部分基因组复制到裂解周期的非二元诱导模型，表明前噬菌体可作为NWS稳定性的情境依赖性调节因子，而非直接的破坏者。这为理解传统发酵乳制品中微生物群落的长期稳定性与功能鲁棒性提供了新见解。
在传统奶酪的世界里，尤其是意大利北部特伦蒂诺格拉纳（Trentingrana）这样受原产地名称保护（PDO）的硬质格拉纳奶酪，其独特的风味和质地深深植根于其使用的自然乳清发酵剂（Natural Whey Starter, NWS）。这些发酵剂并非单一、明确的菌种，而是由包括瑞士乳杆菌（Lactobacillus helveticus）在内的多种嗜热乳酸菌（LAB）构成的复杂微生物群落，通过“回加法”（back-slopping，即用上一批成功的发酵产物作为下一批的接种物）世代相传。这种实践塑造了稳定的微生物生态系统，但也引入了不稳定的“隐形居民”——噬菌体（Bacteriophages, phages）。在乳品工业中，噬菌体攻击常被视为“破坏者”，可能导致发酵失败。然而，越来越多的研究发现，在NWS这样的复杂群落中，噬菌体（尤其是能整合进细菌基因组、处于潜伏状态的温和噬菌体或前噬菌体）与宿主细菌的互动可能更为微妙和持久，甚至可能扮演着生态调节者的角色。一个关键但尚未被充分理解的问题是：在奶酪制作过程中反复出现的非生物胁迫（如温度变化、酸积累、氧化应激等）是否会触发这些“沉睡”的前噬菌体？它们是被统一激活，还是各自拥有独特的“应激响应程序”？这对发酵剂的长期稳定性和奶酪质量又意味着什么？为了回答这些问题，一项发表于《Frontiers in Microbiology》的研究，以特伦蒂诺格拉纳NWS为天然模型，对其中六株瑞士乳杆菌来源的温和噬菌体进行了基因组和功能层面的深入探究。
为开展这项研究，作者运用了几个关键技术方法。首先，从特伦蒂诺格拉纳联盟五个乳品厂收集的NWS样本中分离出392株瑞士乳杆菌，并从中通过丝裂霉素C（MC）诱导，分离出六株温和噬菌体（φCR191, φP185, φCV244, φS16, φS193, φT280）。其次，利用透射电子显微镜（TEM）观察噬菌体形态，并通过噬菌斑裂解试验测定其宿主范围。第三，采用牛津纳米孔技术对噬菌体进行基因组测序，利用一系列生物信息学工具（如PHANOTATE, PHAROKKA, VIBRANT, vConTACT2, ViPTree等）进行基因组注释、比较基因组学、基因共享网络和系统发育分析，特别关注了调控模块的多样性。最后，通过定量PCR（qPCR）靶向噬菌体标志基因，结合检测细胞外噬菌体颗粒，研究了在热激（65°C）、营养限制（培养基稀释）、乳酸和过氧化氢胁迫下，前噬菌体诱导的动力学。
研究结果
3.1 噬菌体宿主范围
六株噬菌体对67株测试的瑞士乳杆菌菌株显示出不同的裂解特异性。其中φCV244和φP185能感染一个以上的菌株，而φCR191、φS16和φS193除了其原始宿主外，仅能裂解另一株菌。这揭示了噬菌体与宿主间的高度特异性，但也存在有限的交叉感染。
3.2 噬菌体形态
透射电子显微镜图像显示，六株噬菌体均具有头尾结构，属于有尾噬菌体纲（Caudoviricetes）。其中五株（φCR191, φP185, φS16, φS193, φT280）呈现肌病毒样（myovirus-like）形态（具有可收缩的尾），而φCV244则呈现长尾病毒样（siphovirus-like）形态（具有不可收缩的长尾）。
3.3 噬菌体基因组特征与比较
六株噬菌体基因组大小在31.5至39.5 kbp之间，GC含量在35.8%至39.3%之间，平均预测开放阅读框（ORF）约72个。BLASTn分析显示它们与已知的瑞士乳杆菌噬菌体φAQ113（NC_019782.1）高度相似（95.29%至98.62%）。基因组注释揭示了典型的噬菌体周期相关基因（如结构蛋白、DNA包装、宿主裂解、溶原/整合、DNA修复与复制），同时含有大量功能未知的假设蛋白。所有噬菌体均被鉴定为温和噬菌体，含有整合酶等溶原相关基因。ViPTree蛋白组树分析进一步确认了它们与φAQ113的聚类关系，而φCV244则与噬菌体phig1e（NC_004305.1）更为接近。
3.3.1 特伦蒂诺格拉纳前噬菌体调控模块的多样化
比较基因组学分析显示，尽管核心的基因组骨架（结构、DNA包装模块）相对保守，但调控溶原-裂解周期转换的基因区域在基因组成和组织上存在显著变异。这些噬菌体编码了不同的Cro/C1样阻遏蛋白、抗阻遏蛋白（如KilAC样、BRO相关）以及辅助DNA结合蛋白（如ArpU样、Rha样、IrrE样）的组合。系统发育分析表明，这些调控蛋白并不严格遵循全基因组相似性的聚类关系，而是与不同乳制品相关噬菌体的同源蛋白交织在一起。这表明调控模块是密切相关的噬菌体之间功能多样化的一个主要轴心，可能对应了它们对胁迫信号的不同感知和响应策略。
3.4 前噬菌体DNA复制研究
通过qPCR检测胁迫处理2小时后噬菌体DNA拷贝数的变化，结合细菌生长曲线监测和4小时后上清液中噬菌体颗粒的检测，系统评估了前噬菌体的诱导动态。研究发现，胁迫暴露引发了噬菌体特异性和剂量依赖性的反应。总体而言，热激（65°C）是最一致的前噬菌体激活诱导因素，在多个噬菌体中引起了DNA拷贝数的显著增加（如φCV244的MHP基因拷贝数增加约1.5 log2FC）。丝裂霉素C（MC）作为阳性对照也有效诱导了所有噬菌体。相比之下，乳酸和氧化应激（H2O2）产生了阈值依赖性或抑制性效应，特别是在较高浓度下。例如，0.2 mM H2O2显著降低了φCV244的DNA拷贝数。营养限制单独作用的影响有限，暗示其可能扮演“许可”而非“触发”角色。值得注意的是，经常观察到噬菌体基因组拷贝数的早期增加与后续病毒颗粒释放之间的“解耦”现象，即DNA复制增加并不总是导致大量噬菌体产生。生长曲线分析显示，热激和MC处理对宿主细菌生长产生了最显著的抑制。
研究结论与讨论
本研究综合表征了来自特伦蒂诺格拉纳NWS的六株瑞士乳杆菌温和噬菌体。这些噬菌体具有狭窄的宿主范围，与参考噬菌体φAQ113基因组高度相似，属于有尾噬菌体纲，但形态（肌病毒样与长尾病毒样）和关键的调控模块构成存在差异。最重要的发现在于，它们对奶酪制作相关胁迫的响应是异质性的，而非统一的。热激是最有效的诱导因素，而其他胁迫则可能根据噬菌体种类和胁迫强度的不同，产生诱导、抑制或无显著影响。
这些结果支持了一个关于前噬菌体诱导的“非二元模型”。传统的观点认为诱导是从溶原状态直接切换到完全裂解周期的“开关”。而本研究揭示的图景更为复杂，涵盖了一系列中间状态——从部分基因组复制、不完全的裂解周期激活，到完全的生产性裂解循环。这些过程强烈受宿主生理状态的制约。例如，氧化应激可能通过损害宿主细胞，间接抑制了前噬菌体的成功激活和扩增。
从生态和应用角度，这项研究具有深远意义。它表明，在特伦蒂诺格拉纳NWS这样的复杂、长期维持的发酵剂中，前噬菌体可能并非简单的“定时炸弹”或“破坏者”，而是充当了“情境依赖性的调节者”。它们对胁迫的异质性响应，使得微生物群落能够在面对环境波动时，通过调节特定细菌亚群的数量和动态，来维持整体的稳定性和功能鲁棒性。这种应激响应的前噬菌体动力学很可能在回加法过程中影响着微生物种群的更替和恢复力，最终为传统发酵乳制品的长期稳定性和独特品质形成做出了贡献。未来的研究需要进一步在群落水平上验证这些相互作用，并探索如何利用这些知识来更好地管理和优化传统发酵食品的生产。
