单细胞成像揭示裂解性噬菌体生活周期的动态异质性

《Cell Systems》：Single-cell imaging of the lytic phage life cycle in bacteria

当一个裂解性噬菌体感染细菌细胞后，它会劫持宿主的资源进行复制，最终裂解细胞以释放新的病毒颗粒。由于噬菌体是专性寄生生物，感染过程的每个阶段都依赖于宿主细胞的生理参数。鉴于遗传学上相同的细菌群体内部存在固有的生理变异性，研究人员探究了噬菌体感染动态如何反映这种异质性。本研究引入了一种延时成像分析方法，用于研究单个T7噬菌体对单个细菌的各个感染步骤的动态过程。这种高通量、时间分辨的方法使研究人员能够同时监测多个细胞中的感染进程，揭示了每个步骤中显著的异质性，并阐明了感染动态与被感染细胞生理状态之间的相关性。对竞争性噬菌体群体的模拟表明，感染动态的异质性会显著影响噬菌体的适应度，凸显其作为噬菌体-细菌相互作用潜在进化驱动力的可能性。

一、研究背景与科学问题

噬菌体（bacteriophages）是感染细菌的病毒，在塑造自然界细菌群落中扮演关键角色，并在医学和生物技术领域作为生物防治剂具有重要前景。其中，裂解性噬菌体（lytic bacteriophages）因其对抗抗生素耐药性感染的潜力而备受关注。由于这些病毒依赖宿主分子机器和前体进行增殖，从感染到裂解的整个过程高度依赖于靶细胞的生理特征，例如宿主代谢和受体密度影响吸附速率，宿主细胞资源和机器影响病毒粒子复制速率及裂解剂的产生。近期单细胞研究揭示，基因相同的细菌细胞间存在显著的生理变异性，这引发了关于这种多样性如何影响噬菌体感染动态，以及感染动态的变异性如何在群体水平上影响噬菌体消除靶细菌的整体效果等问题。

当前研究噬菌体感染步骤动态的大多数方法依赖批量培养（bulk culture）和组学分析，缺乏分析细胞间异质性所需的单细胞分辨率。虽然创新性的改进可通过极高稀释倍数实现单细胞分辨率，但无法将测量的感染参数与宿主细胞生理联系起来。冷冻电镜（cryo-electron microscopy, cryo-EM）虽能高分辨率研究结构方面，但无法提供时间分辨信息以监测单个感染步骤的进展及其相互关系。近期的单细胞研究揭示了温带噬菌体裂解-溶原转换机制的前所未见见解，但很少有研究尝试在单细胞水平分析裂解性噬菌体感染，并理解宿主细胞生理多样性如何影响感染周期。因此，需要一种方法能够：（1）以时间分辨方式追踪单个噬菌体针对单个活细菌发起的每个感染阶段；（2）维持时空均匀的环境以隔离宿主细胞生理内在变化对感染参数的影响；（3）具备足够高的通量以量化靶细菌细胞间的详细分布。

二、研究方法概述

为解决上述挑战，研究人员开发了一套定制化方法。该方法利用微流控平台（microfluidic platform），该平台经过改造（改良型“母机”装置，mother machine device），能够在均匀生长条件下维持分离的目标细胞群体，从而追踪单个细胞被单个噬菌体感染和裂解的感染动态。结合高速扫描时间分辨显微镜（high-speed scanning time-resolved microscopy）以及在模式系统T7噬菌体上使用的荧光标记组合（包括快速成熟的黄色荧光蛋白mVenus NB作为衣壳表达报告子，以及DNA结合染料SYTOX Orange标记噬菌体基因组以观察吸附和注入），研究人员能够跟踪从噬菌体吸附到细胞裂解的单个感染事件。图像数据通过机器学习模型（machine-learning model）处理，该模型使用合成显微照片训练，以量化细胞生理学参数（如大小、形状、分裂时间和生长速率）和感染标志物（如吸附、基因组注入、宿主接管、衣壳表达和裂解）。研究人员还建立了一个数学模型（mathematical model），考虑病毒基因组复制、转录和翻译，以解析蛋白质生产动力学。

三、主要研究结果

1. 单细胞成像揭示T7噬菌体感染周期

研究人员通过荧光显微镜监测T7噬菌体生活周期的不同步骤。通过在T7基因组中引入衣壳基因表达的荧光报告子（gp10A-B，mVenus NB），并用DNA结合染料SYTOX Orange染色，实现了对噬菌体吸附和随后基因组注入的可视化。在低感染复数（multiplicity of infection, MOI）下操作，确保了首次裂解事件源于单个噬菌体感染单个细菌。结果显示，该方法能清晰分辨吸附（t₀）、基因组注入（t₁）、生长停滞（t₂）、衣壳蛋白表达开始（t₃）和细胞裂解（t₅）等关键时间点。

2. 基因组注入动力学呈现多种进入模式

通过定量随时间变化的染色噬菌体DNA荧光信号，研究人员发现基因组注入过程存在显著的异质性。跨多个吸附事件的强度趋势揭示了五种广泛的基因组注入动力学类别，包括缓慢恒定速率进入、经典三步进入以及极高速率进入等模式。T7的平均注入持续时间为4.4分钟，T7^*为5.7分钟，但变异系数（coefficient of variation, CV）高达50%-55%，表明不同感染事件间差异巨大。

3. 宿主细胞关闭与病毒接管的动力学在不同感染事件间高度一致

研究人员通过细胞长度变化确定生长停滞点，并通过荧光报告子监测衣壳蛋白表达。结果表明，生长停滞动力学在不同感染事件间非常稳健，细胞在3-4分钟内从感染前生长速率过渡到完全停止。相比之下，衣壳报告子的表达表现出相当大的变异性，最大强度（I_max）在不同事件间可相差近一个数量级。这种变异性与生长停滞时的细胞大小（L_GA）强相关。从生长停滞到裂解的时间（t₂到t₅，平均9.6分钟，CV=15%）以及从衣壳生产开始到裂解的时间（t₃到t₅，平均7.7分钟，CV=17%）分布非常窄，表明一旦噬菌体接管细胞，后续事件时序高度可重复。

4. 噬菌体蛋白生产动力学与宿主生理变异性紧密耦合

通过建立的数学模型拟合实验数据，研究人员提取了最大翻译速率（λ₂）和噬菌体基因组复制速率（λ₁）。研究发现，λ₂与最大YFP值强相关，并与生长停滞时的细胞大小正相关，表明被感染细胞的核糖体数量是控制噬菌体蛋白生产的关键变量。这解释了为什么在均匀环境下，仅因被感染细胞的生理差异，爆发量（burst size）理论上可变化近一个数量级。此外，λ₁和λ₂之间存在显著的负相关，暗示病毒基因组复制与翻译之间存在负反馈调节机制。

5. 细胞膜的逐渐穿孔先于感染细胞的突然裂解

通过高频成像（100 Hz），研究人员观察到裂解过程分为两个截然不同的阶段：首先是一个较慢的“穿孔”（perforation）阶段，表现为细胞内物质逐渐泄漏，相衬亮度增加，平均持续约6.56秒（CV=50%）；随后是快速的“裂解”（lysis）阶段，即细胞包膜的结构性崩解。穿孔在所有裂解细胞中均被观察到。

6. 裂解时间的异质性主要由感染早期阶段的变异性驱动

通过构建完整的感染周期时间线，研究人员发现裂解时间的总体变异性（CV=21%）与注入结束到裂解的时间间隔（t₁到t₅）的变异性相当。相反，生长停滞到裂解（t₂到t₅）和衣壳生产开始到裂解（t₃到t₅）的持续时间非常稳定。这表明裂解时间变异性的主要来源在于感染初期直至细胞生长停止的阶段。

7. 裂解时间变异性可为噬菌体群体提供适应性优势

通过随机基于主体的模拟（stochastic agent-based simulations）比较不同裂解时间分布的噬菌体竞争，研究人员发现，即使平均裂解时间相同，具有更大标准差（standard deviation）的突变体具有竞争优势；在平均值和标准差相同的情况下，具有负偏态（negative skew）的突变体也具有竞争优势。这表明裂解时间的均值不足以预测噬菌体的竞争能力，分布的高阶矩可以显著改变其适应性优势。

四、总结与讨论

本研究首次在单噬菌体-单细菌分辨率下量化了T7噬菌体裂解周期中各个步骤的动力学。研究揭示，裂解时间的主要变异性来源于感染早期阶段（从病毒DNA进入到细胞生长停滞），而宿主接管后的第二阶段则高度一致。单细胞轨迹分析揭示了五种不同的基因组进入动力学类别。在宿主接管后的生产阶段，尽管事件时序高度稳健，但噬菌体蛋白表达动力学和总产量却存在巨大变异性，这种变异性与宿主细胞大小（反映翻译资源）密切相关，建立了宿主细胞生理与噬菌体爆发量之间的强机制联系。裂解过程包含一个约6.6秒的穿孔期。模拟结果清楚地表明，平均表型值（如平均裂解时间）不足以预测噬菌体群体的适应性优势，分布的高阶矩可产生显著影响，这提出了一个引人入胜但目前探索不足的可能性：噬菌体感染表型的变异性可能受到强烈的选择压力。

研究结论翻译如下：

本研究报告了首个在单噬菌体-单细胞分辨率下量化T7噬菌体裂解周期中各个步骤动力学的工作。该检测方法提供了一种正交于批量分析和结构研究的量化噬菌体-细菌相互作用的方式，能够在精确控制的环境中跨多个感染事件进行动力学测量，同时保持每个事件的个体性，以评估变异性以及噬菌体表型与相应被感染细胞之间的相关性。研究人员发现裂解时间的主要变异性源自感染早期阶段，从病毒DNA进入到细胞生长停滞，而感染过程的第二部分，从宿主接管到细胞裂解，则以显著的稳定性进行。对单细胞轨迹的分析揭示了五种不同的基因组进入动力学类别。在宿主接管后，感染进入一个以快速合成噬菌体蛋白为特征的生产阶段。尽管事件时序高度稳健，但在此阶段观察到表达动力学和总产量的巨大变异性。噬菌体总蛋白产量与生长停滞宿主的大小之间的强相关性表明，这种变异性可能源于被感染细胞间的生理差异。通过数学模型，研究人员发现噬菌体蛋白生产效率强烈依赖于细胞的翻译资源，后者与细胞大小成比例，从而提供了宿主细胞生理与噬菌体爆发量之间的强机制联系。在生产阶段之后，感染以宿主细胞裂解告终。研究人员发现，T7感染细胞的物质损失大约在细胞包膜完全崩解前6.6秒开始，这一阶段被称为穿孔期。感染事件的有序序列，从基因组进入到细胞裂解，在每个步骤都表现出可测量的变异性。准确量化这些相对变异性的来源及其在感染事件间的相关性，不仅对于理解控制噬菌体感染结果的潜在分子机制至关重要，而且可能具有重要的进化意义。模拟结果清楚表明，平均表型值（如平均裂解时间）不足以预测噬菌体群体的适应性优势，分布的高阶矩可产生显著影响。这些发现开辟了一个引人入胜且目前探索不足的可能性：噬菌体感染表型的变异性可能处于强大的选择压力之下，引发了关于进化如何在不同的情景下塑造它的迷人问题。