噬菌体转录因子 Gp53 如何抑制铜绿假单胞菌分裂？

Gp53 是 ΦPA1019 的早期高表达调控蛋白

论文从污水样品中分离出铜绿假单胞菌噬菌体 ΦPA1019，其基因组为 94,085 bp，含 132 个开放阅读框。研究者依据基因组注释挑选了 12 个候选早期基因，用 RT-qPCR 在感染后 0、5、30、60 分钟检测转录水平，并以晚期衣壳基因 gp5 作为内参。gp46、gp47、gp53、gp60、gp63 和 gp111 在感染后 5 分钟显著上调，其中 gp53 表达最高，相对对照基因升高 4.4 倍。

序列分析显示 gp53 编码一个 126 个氨基酸的 RHH 家族 DNA 结合蛋白，与噬菌体 P22 的 Arc 阻遏蛋白和铜绿假单胞菌自身的转录调控因子 AmrZ 同属该家族。这一结构线索提示 Gp53 可能像它们一样直接结合 DNA 并调控转录，而不是通过蛋白-蛋白相互作用干扰宿主。

需要区分的是，这些早期基因的筛选与表达量测定来自噬菌体感染过程，而后续功能实验主要依赖把 gp53 克隆到 IPTG 诱导的质粒上过表达。作者用 RT-qPCR 比较了两种条件：0.5 mM IPTG 诱导下的 gp53 转录水平约为感染后 5 分钟时的 2 倍，二者处于同一数量级，因此作者认为该过表达系统适合做功能表征。这一比较支持但不等同于“Gp53 在自然感染中产生完全相同的效应”。

过表达 Gp53 同时压低生长、运动与生物膜

把 gp53 在铜绿假单胞菌 PAO1 中过表达后，多个表型同时改变。生物膜生物量约为野生型 PAO1 的 1/3，即降低约 3 倍；群集运动（swarming）和泳动（swimming）都明显减弱。作为对照，同样过表达 gp63 对这些表型没有显著影响，说明这些变化不是过表达质粒或诱导条件本身造成的。

生长曲线显示 gp53 过表达菌株的增殖速率明显低于野生型，琼脂点样实验也显示生长受抑。诱导 6 小时后，激光共聚焦显微镜观察到细胞显著伸长，呈丝状结构。生物膜、运动性和丝状化分别对应不同的宿主过程，因此这些结果本身只能说明 Gp53 影响范围广，不能单独证明它们由同一分子事件引起。

这些表型实验的统计单位是三次独立生物学重复（n = 3），生物膜以 560 nm 吸光度定量，运动性以运动圈直径衡量。它们描述的是实验室 LB 培养基、特定诱导浓度下的表现，而不是感染过程中的细菌状态。

Z 环定位异常把表型连到分裂装置

铜绿假单胞菌的分裂依赖 Z 环，Z 环主要由 FtsZ 及其相关蛋白（如 ZapA）组成。由于直接给 FtsZ 加荧光标签可能干扰其功能，研究者改用被 FtsZ 招募到分裂位点的 ZapA-GFP 来观察 Z 环位置。在野生型 PAO1 中，ZapA-GFP 荧光主要定位在细胞中部，对应正常的 Z 环组装；而在 gp53 过表达菌株中，信号弥散或定位错误，说明 Z 环的空间组织被破坏。

这一实验把“细胞变长”这一形态观察与具体的分裂装置缺陷联系起来：细胞伸长是分裂受阻的结果，而 Z 环定位异常提供了可能的直接原因。不过，荧光定位反映的是 ZapA 的分布，论文据此推断 Z 环组装异常，而非直接测量 FtsZ 的聚合状态。

转录组显示 Gp53 是全局性调控因子

RNA-seq 比较了携带空载体 pMMB67EH 与 pMM-gp53 的 PAO1 菌株，以 |log2 fold change| ≥ 1 且校正 P 值 < 0.05 为标准，共鉴定出 3,342 个差异表达基因，其中 1,642 个上调、1,700 个下调。KEGG 富集分析显示差异基因集中在 20 条通路，包括双组分系统（140 个基因）、ABC 转运（126 个）、生物膜形成（72 个）、群体感应（62 个）和鞭毛组装（31 个）。

在具体基因层面，与 Z 环形成相关的 ftsZ、ftsA、zapA 表达下降，而分裂抑制因子 sulA 和 minC/D/E 表达上升；生物膜合成相关基因（pslB/C、amrZ、algC、cdrB）和运动相关基因（flgA/B/C、fleQ、motC/D、pilA/U）下调；毒力相关基因 lasR/rhlR、exsC、phzA1/A2、rhlA/B 也下降。研究者用 RT-qPCR 验证了 9 个代表性基因，以 rpsL 为内参，结果与 RNA-seq 一致。

这些数字说明变化的规模，但不能说明每个基因都被 Gp53 直接结合。转录组反映的是过表达条件下的表达差异，其中既包含直接调控，也包含下游继发效应；区分二者需要独立的 DNA 结合数据。

DAP-seq 与 EMSA 锁定直接靶启动子

为找出 Gp53 直接结合的位点，研究者做了全基因组 DAP-seq：Gp53 特异结合 PAO1 基因组中的 78 个 DNA 区域，其中 62 个（79.5%）位于基因启动子的核心调控区，76.92% 的结合峰集中在转录起始位点上游 −200 至 +200 bp 范围内。靶基因包括生物膜相关基因 amrZ、pslA、pelA，运动相关基因 flgA、fleQ，毒力基因 phzA2、wzz2，c-di-GMP 代谢酶基因 mucR、roeA，以及分泌系统基因 tssA1、hsiA2。

体外的电泳迁移率变动分析（EMSA）进一步验证了直接结合：纯化的 Gp53 蛋白与 fleQ、pslA、amrZ 启动子探针呈剂量依赖性结合，以 rpoC 作为非靶标对照；对 ftsZ 全长启动子及其两个推定结合位点也观察到特异性结合，而对 minC 启动子没有结合。

DAP-seq 给出的是全基因组范围内的结合位点集合，EMSA 给出的是体外特定探针的结合证据，两者回答的问题不同但相互补充：前者说明 Gp53 在基因组上可能作用于哪些位置，后者说明 Gp53 蛋白能否直接结合这些序列。它们共同支持 Gp53 作为转录抑制因子直接作用于 ftsZ 等启动子，但不直接测量结合后的转录抑制强度。

补回 ftsZ 能逆转分裂与 Z 环缺陷

论文的关键因果证据来自补救实验：在 gp53 过表达背景下同时表达 ftsZ，Z 环定位异常和细胞分裂缺陷得到恢复。摘要将这一点表述为 Gp53 介导的 ftsZ 抑制破坏 Z 环形成并抑制细胞分裂，而 ftsZ 表达可恢复这些缺陷。

这一比较把 ftsZ 从众多差异基因中区分出来：转录组里 ftsZ、ftsA、zapA 同时下降，但只有补回 ftsZ 被报告能逆转分裂相关表型，因此 ftsZ 抑制在这条通路中具有特殊地位。需要注意的是，补救实验针对的是分裂与 Z 环表型；论文并未报告补回 ftsZ 是否同样恢复生物膜、运动性或毒力基因表达，因此不能把 ftsZ 的作用外推到全部表型。

结构建模与突变定位 DNA 结合残基

研究者用 AlphaFold Server 预测了 Gp53-DNA 复合物的三维结构，并用 PyMOL 分析极性接触与氢键，据此标注参与 DNA 结合的关键氨基酸残基。同时，他们通过重叠延伸 PCR 构建 gp53 定点突变体，将突变基因分别克隆到 pET-28a-SUMO 用于蛋白表达、克隆到 pMMB67EH 用于表型分析，以检验这些残基对 DNA 结合能力的影响。

这部分工作的定位是结构-功能对应关系：它指出哪些残基可能参与 Gp53 与 DNA 的相互作用，并用突变影响 DNA 结合能力来支持这一判断。结构来自计算预测而非实验解析，因此残基层面的结论应视为模型指导下的假设加突变验证，而不是高分辨率结构测定结果。

这项研究最值得关注的是直接靶向而非广谱效应

把各层证据放在一起，这篇论文最有价值的判断不是“Gp53 影响很多基因”，而是它把一条可追踪的链条接了起来：Gp53 直接结合 ftsZ 启动子，ftsZ 转录下降，Z 环定位异常，细胞分裂受阻，补回 ftsZ 可逆转分裂缺陷。3,342 个差异基因和 78 个结合区域说明作用范围广，但真正被实验串成因果链的是 ftsZ 这一支。

因此读者在评估后续工作时，应优先区分两类主张：Gp53 直接抑制的启动子（ftsZ、fleQ、pslA、amrZ 等，有 DAP-seq 与 EMSA 支持）与仅由转录组相关的变化（如 lasR/rhlR、phzA1/A2 等毒力基因的下调，论文未逐一证明直接结合）。前者可作为机制起点，后者更适合作为待验证的下游效应。

这一优先级也界定了结论的适用边界：功能证据主要来自 PAO1 在 LB 培养基中的质粒过表达，诱导水平约为感染后 5 分钟的 2 倍；结构证据来自计算预测。要把 Gp53 视为感染过程中协调抑制分裂与毒力的通用机制，还需要在自然感染条件下验证这些直接靶点的抑制程度。

本文引用的文献