The contribution of abortive infection to preventing populations of Lactococcus lactis from succumbing to infections with bacteriophage
流产感染机制在乳酸乳球菌抵御噬菌体群体消亡中的作用研究
论文来源：PLOS ONE，2024年第19卷第4期，文章编号e0298680，DOI：10.1371/journal.pone.0298680
一、摘要内容总结
乳制品发酵常因烈性噬菌体污染造成乳酸乳球菌发酵失效，流产感染（AbiZ）是乳酸菌关键抗噬菌体免疫机制：菌体被侵染后提前裂解，阻断子代噬菌体完整增殖。本研究以两株携带/缺失AbiZ的乳酸乳球菌（NCK4/NCK5、IL6/IL7）、两类工业常见噬菌体p2、P335为试验对象，结合数学动力学模型与连续传代体外培养，解析AbiZ对菌群的保护边界与局限。结果表明：仅当AbiZ阻断效率高于94%且无逃逸突变噬菌体时，才能完全保护菌群；噬菌体易产生衣壳蛋白错义突变形成Abi逃逸株，单一Abi机制不足以长期维持菌群存活，菌体后续进化出包膜抗性突变才是种群长期存续核心保障；不同菌株体系进化速率存在差异，p2体系包膜抗性出现更快，P335体系需更长传代周期才能产生稳定抗性。该研究厘清AbiZ免疫的生态局限性，揭示噬菌体逃逸突变、细菌包膜抗性的协同军备进化规律，为乳制品发酵抗噬菌体菌种选育提供理论依据。
二、关键词
三、研究目的
1. 构建包含宿主菌、原始噬菌体、Abi逃逸噬菌体、抗性菌株的质量作用动力学微分模型，量化Abi阻断效率、突变速率等关键参数对菌群存亡的影响。
2. 短期液体培养试验，对比AbiZ⁺与AbiZ⁻菌株在p2、P335侵染下活菌、噬菌体密度动态变化。
3. 鉴定p2ev、P335ev两类Abi逃逸噬菌体突变位点，明确衣壳蛋白突变是逃逸分子基础。
4. 开展5代连续传代长期进化试验，解析菌群、噬菌体长期共进化规律，区分Abi与包膜抗性各自贡献。
5. 对比两套菌-噬菌体体系进化差异，阐明菌株背景对抗性进化速度的影响。
6. 系统评价AbiZ单一机制的应用局限，完善乳制品抗噬菌体菌种开发理论。
四、研究思路
1. 菌株与噬菌体构建：无AbiZ宿主NCK4、IL6，电转质粒pTRK914获得AbiZ⁺菌株NCK5、IL7；实验室分离获得逃逸噬菌体p2ev、P335，测序定位突变位点。
2. 生长参数测定：Bioscreen高通量仪测定各菌株最大生长速率、突变株适应性代价，拟合模型参数。
3. 7h短期共培养：定时梯度稀释涂布，分别计数活菌CF、总噬菌体、逃逸噬菌体PFU，绘制短期种群动态曲线。
4. 数学模拟：Berkeley Madonna求解微分方程，分无逃逸噬菌体、有逃逸无包膜抗性、同时存在逃逸噬菌体与抗性菌三种情景模拟菌群结局。
5. 24h静态培养验证模拟预测，对比初始与终末菌、噬菌体丰度。
6. 5次连续24h传代长期进化试验，每代计数、分离单菌落检测包膜抗性表型。
7. 噬菌体全基因组测序，对比原始与逃逸株衣壳蛋白序列差异；波动试验测定抗性突变频率。
8. 多组统计学对比模型模拟与实测数据拟合度，总结共进化规律。
五、研究亮点
1. 首次建立兼顾Abi机制、噬菌体逃逸、细菌包膜抗性的多组分动力学数学模型，量化Abi发挥保护作用的临界效率阈值94%。
2. 证实AbiZ仅短期缓冲噬菌体压力，一旦产生衣壳突变逃逸噬菌体，单一流产感染机制无法阻止菌群衰退。
3. 明确噬菌体逃逸核心突变位点：两类逃逸噬菌体均发生主要衣壳蛋白错义突变，揭示Abi识别靶标为噬菌体衣壳结构。
4. 区分两套试验体系进化差异：NCK5+p2体系短时间即可产生高适应性包膜抗性，IL7+P335体系抗性突变进化周期更长，存在过渡敏感亚群。
5. 长期传代证明即便存在抗性菌，噬菌体仍可持续低水平增殖，抗性属于“渗漏型”，无法完全清除噬菌体。
6. 整合数学模拟、高通量生长、短期/长期体外进化、基因组测序多维证据，完整阐明乳酸菌与噬菌体协同军备进化生态机制。
六、可延伸的方向
1. 构建携带多重Abi系统的工程菌株，评估多Abi叠加抑制噬菌体逃逸效果。
2. 模拟乳制品真实发酵基质（牛乳）开展共进化试验，验证模型工业场景适用性。
3. 筛选靶向噬菌体衣壳的天然抑制剂，阻断逃逸突变产生。
4. 扩大噬菌体种类（多血清型）、混合噬菌体体系，完善多噬菌体共存动力学模型。
5. 解析包膜抗性对应的受体蛋白突变位点，改造宿主受体开发永久抗噬菌体工业菌株。
6. 基于动力学模型预测发酵罐噬菌体污染爆发时间，建立车间污染预警模型。
七、测量的数据、配套图表与对应研究价值
1. Bioscreen测得各菌株最大生长速率、突变株适应性代价数据，配套表1模型参数表；精准获取vn、va、突变率等核心建模参数，量化包膜抗性带来的生长适合度损耗，为数学模拟提供真实试验基础。
2. 7h短期共培养活菌、原始/逃逸噬菌体时序丰度数据，配套图1；直观证明AbiZ⁺菌株短期种群保护优势，同时观测逃逸噬菌体随时间持续增殖趋势，验证Abi机制短期局限性。
3. 三组情景数学模拟菌/噬菌体丰度柱状图，配套图3；计算得出Abi起效临界阻断效率q=0.94，清晰区分无逃逸、仅逃逸、逃逸+抗性三种菌群存亡结局。
4. 24h培养初始/终点菌、噬菌体计数数据，配套图4；体外实测匹配模型预测，证实逃逸噬菌体出现后Abi保护作用大幅衰减。
5. 5代连续传代长期动态曲线，配套图5、图6；对比两套体系抗性进化速度差异，说明菌株遗传背景决定菌群长期存续能力。
6. 噬菌体基因组测序突变位点、包膜抗性筛选菌落统计数据，配套补充表格；从分子层面解释噬菌体逃逸机制，证实衣壳蛋白突变是核心逃逸通路。
7. 质量作用动力学微分模型通路示意图，配套图2；可视化宿主、原始噬菌体、逃逸噬菌体、抗性菌株相互转化与侵染逻辑，统一全文理论框架。
八、得出的结论
1. 流产感染AbiZ对乳酸乳球菌的保护存在严格条件：阻断效率≥94且不存在可逃逸突变噬菌体时，才能避免菌群被噬菌体完全清除。
2. p2、P335噬菌体均可通过主要衣壳蛋白错义突变产生Abi逃逸株，该突变会削弱Abi识别功能，使单一Abi机制失效。
3. AbiZ仅能短期延缓菌群消亡，种群长期稳定存活依赖后续进化产生包膜抗性突变，包膜抗性是抵御逃逸噬菌体的核心屏障。
4. 不同菌-噬菌体组合进化特征存在差异：NCK5+p2体系包膜抗性进化快、适合度代价小；IL7+P335体系抗性突变生成周期长，短期传代易出现菌群数量暴跌。
5. 乳酸乳包膜抗性属于渗漏型，即使抗性菌株占主导，噬菌体仍可持续低水平复制，发酵体系无法彻底消除噬菌体污染压力。
6. 仅依靠单一Abi流产感染机制不足以构建稳定抗噬菌体发酵菌种，需结合包膜抗性筛选或多重Abi叠加改造，才能适配乳制品长期发酵生产。
九、芬兰Bioscreen仪器微生物生长曲线数据的研究意义
1. 96孔高通量同步测定野生、Abi⁺、包膜抗性多菌株生长曲线，一次性平行重复消除分批培养温度、溶氧系统误差。
2. 仪器连续自动采集OD600浊度，完整记录菌株全生长周期，精准计算最大比生长速率，量化抗性突变带来的适合度增减。
3. 直接提供vn、va等动力学基础参数，为微分数学模拟提供真实试验校准数据，提升模型预测准确性。
4. 全程恒温统一培养环境，浊度检测重复性远优于人工单点分光计数。
5. 输出标准化生长时序数据，作为菌株适应性差异定量核心支撑图表。
6. 该高通量检测流程可标准化用于抗噬菌体乳酸菌突变株生长适合度批量筛选。