Systematic Analysis of c-di-GMP Signaling Mechanisms and Biological Functions in Dickeya zeae EC1
玉米迪基氏菌EC1中c-di-GMP信号传导机制与生物学功能的系统分析
来源：mbio November/December 2020 Volume 11 Issue 6 e02993-20
1.论文摘要核心内容
玉米迪基氏菌（Dickeya zeae）是重要的侵袭性植物病原细菌，会给香蕉和水稻种植造成重大经济损失。前期研究发现，该菌EC1菌株中负责c-di-GMP合成与降解的信号蛋白（环化酶/磷酸二酯酶），在细菌从固着到运动的表型转换中发挥重要作用。为明确这些c-di-GMP信号蛋白之间是否存在协同效应，本研究以EC1为亲本菌株，通过连续框内缺失的方法，分别构建了编码双鸟苷酸环化酶（DGC，负责c-di-GMP合成）和磷酸二酯酶（PDE，负责c-di-GMP降解）的系列突变菌株。结果显示，全部推定DGC基因的完全缺失，会导致细菌运动能力显著增强、生物被膜形成能力完全丧失，但不影响其致病性和毒力因子的产生；相反，全部c-di-GMP PDE基因的缺失，会使细菌完全丧失运动能力，无法侵染水稻种子。通过测定所有突变体的胞内c-di-GMP浓度和游泳运动能力，本研究提出c-di-GMP以浓度依赖的方式，通过多级程序调控细菌的游泳行为，二者的关系可描述为L型回归曲线，这一特征与沙门氏菌、新月柄杆菌等其他已报道的细菌物种存在显著差异。进一步分析鉴定出PDE10355、DGC14945和PDE14950这3个c-di-GMP信号蛋白，它们在影响EC1菌株的全局c-di-GMP库中发挥主导作用。本研究结果凸显了不同细菌物种中c-di-GMP调控回路的复杂性与可塑性。
2.关键词（中文）
玉米迪基氏菌、环二鸟苷酸（c-di-GMP）、生物被膜形成、游泳运动能力、连续框内缺失、固着-运动表型转换、生物被膜、植物病原菌、群体感应、根际定殖微生物
3.研究目的
① 突破前期单基因突变研究的局限性，系统解析玉米迪基氏菌EC1菌株中全部c-di-GMP代谢基因（DGC和PDE）的协同作用与全局调控效应，明确c-di-GMP全局库在细菌生理表型与致病性中的核心作用；
② 揭示c-di-GMP浓度与细菌游泳运动、生物被膜形成等关键表型之间的剂量-效应关系，明确该菌中c-di-GMP的调控模式，填补该病原菌中c-di-GMP全局调控机制的研究空白；
③ 鉴定出调控EC1菌株全局c-di-GMP水平的核心“管家”DGC和PDE基因，明确其在c-di-GMP信号网络中的主导地位；
④ 阐明c-di-GMP调控玉米迪基氏菌致病性的分子机制，明确其调控毒力的关键下游表型，为该病原菌的绿色防控提供新的分子靶点与理论依据；
⑤ 丰富革兰氏阴性菌中c-di-GMP调控模式的多样性认知，为其他细菌中复杂c-di-GMP调控机制的研究提供可复制的方法学参考。
4.研究思路
本研究以水稻基腐病、香蕉细菌性软腐病的致病菌玉米迪基氏菌EC1为研究对象，以系统解析c-di-GMP信号传导的全局调控机制为核心目标，按以下逻辑开展系统性研究：
第一步，生物信息学预分析：对EC1菌株基因组中编码c-di-GMP代谢蛋白的基因进行系统解析，明确12个含GGDEF结构域的DGC、4个含EAL/HD-GYP结构域的PDE，以及3个双功能结构域蛋白的序列特征与结构域组成；
第二步，系列突变体构建：以EC1野生型为亲本，通过连续框内缺失的方法，分别构建DGC全缺失突变体（cdG0，15ΔDGC）和PDE全缺失突变体（7ΔPDE），同时构建系列梯度缺失中间突变体，为后续表型与机制分析提供核心实验材料；
第三步，全局表型系统分析：测定DGC全缺失、PDE全缺失突变体的核心表型，包括游泳/群集运动能力、生物被膜形成能力、胞外酶与玉米素（zeamine）等毒力因子的产生能力，以及对水稻种子的致病性与侵染定殖能力；
第四步，分子与形态学验证：通过qRT-PCR测定运动、生物被膜相关功能基因的转录水平，结合透射电子显微镜（TEM）观察细菌鞭毛的数量与组装形态，从转录和细胞层面验证表型变化的内在机制；
第五步，剂量-效应关系解析：通过LC-MS精准定量系列梯度缺失突变体的胞内c-di-GMP浓度，结合对应的游泳运动能力数据，拟合c-di-GMP浓度与运动表型的回归曲线，明确二者的调控规律；
第六步，核心调控基因鉴定：通过单基因突变体的表型测定、胞内c-di-GMP浓度定量，结合体外HPLC酶活实验，鉴定出调控全局c-di-GMP水平的核心DGC和PDE蛋白，并验证其催化活性；
第七步，下游调控通路初探：鉴定出c-di-GMP的下游受体蛋白YcgR同源物，通过基因回补实验验证其在高c-di-GMP水平下对细菌运动的调控作用，完善信号传导通路；
最后，整合所有实验数据，明确玉米迪基氏菌中c-di-GMP的全局调控机制、物种特异性调控模式，以及其调控致病性的核心途径，形成完整的研究结论。
5.研究亮点
① 首次在玉米迪基氏菌中完成了全部c-di-GMP合成与降解基因的系统性敲除，构建了DGC全缺失的c-di-GMP零突变体（cdG0）和PDE全缺失突变体，突破了前期单基因突变研究的局限性，完整揭示了c-di-GMP全局库对细菌生理与毒力的双向调控效应；
② 发现了玉米迪基氏菌中c-di-GMP对运动表型的独特调控模式，即c-di-GMP浓度与游泳运动能力呈L型回归曲线关系，这与新月柄杆菌的倒U型、沙门氏菌的调控模式存在显著差异，拓展了革兰氏阴性菌中c-di-GMP调控模式的认知，凸显了该信号系统在不同细菌中的复杂性与物种可塑性；
③ 修正了该菌中c-di-GMP调控毒力的传统认知，明确其并非通过调控核心毒力因子玉米素的合成影响致病性，而是主要通过调控鞭毛介导的运动能力，进而影响细菌对水稻种子的侵染与定殖，厘清了c-di-GMP调控该病原菌毒力的核心分子机制；
④ 鉴定出3个主导全局c-di-GMP水平的核心信号蛋白（DGC14945、PDE10355、PDE14950），其中DGC14945为非经典GGDEF结构域的SGDEF型环化酶，通过体内外实验验证了其催化活性，为该病原菌的绿色防控提供了全新的关键分子靶点；
⑤ 建立了“连续框内缺失构建梯度突变体-胞内c-di-GMP绝对定量-表型拟合分析”的研究范式，为其他携带多拷贝DGC/PDE基因的细菌中，c-di-GMP全局调控机制的研究提供了可复制的方法学参考。
6.可延伸的方向
① 深入解析核心DGC/PDE蛋白（DGC14945、PDE10355、PDE14950）的酶活调控机制，明确其感应的环境信号、上游调控因子，以及蛋白结构对催化活性的影响，揭示c-di-GMP信号网络响应环境变化的分子机制；
② 系统鉴定玉米迪基氏菌中c-di-GMP的全部受体蛋白，解析各受体与c-di-GMP结合后的下游调控通路，完整构建该菌中c-di-GMP从信号合成、感知到表型输出的全信号网络；
③ 探究c-di-GMP信号系统与该菌中已报道的AHL群体感应、多胺信号、双组分系统VfmIH、转录因子SlyA/Fis等毒力调控系统的交叉调控关系，构建该病原菌的全局毒力调控网络；
④ 基于鉴定的核心c-di-GMP调控靶点，开发靶向DGC/PDE酶活的小分子抑制剂，评估其对玉米迪基氏菌致病性的抑制效果，为水稻、香蕉细菌性病害的绿色防控提供新型抑菌剂；
⑤ 探究c-di-GMP在玉米迪基氏菌与水稻/香蕉宿主互作过程中的动态变化与调控作用，明确其在根际定殖、宿主侵染、系统定殖全侵染过程中的功能，揭示病原菌-宿主互作的新机制；
⑥ 分析不同地理来源、不同宿主分离的玉米迪基氏菌菌株中c-di-GMP代谢基因的保守性与多样性，比较其调控模式的差异，揭示该病原菌的环境适应性进化机制；
⑦ 解析c-di-GMP对玉米迪基氏菌生物被膜形成、环境胁迫抗性、根际种间互作的调控机制，明确其在田间环境存活、定殖中的作用，为该病害的田间生态防控提供理论依据。
7.测量的数据、研究意义及对应图表
① c-di-GMP代谢酶的结构域特征与系列缺失突变体的基因型数据：通过生物信息学分析明确了19个推定c-di-GMP信号蛋白的结构域组成，以及连续缺失构建的15个DGC梯度突变体、7个PDE梯度突变体的基因型，数据来自Fig.1（A为蛋白结构域结构，B为突变体与对应基因型）。研究意义：明确了玉米迪基氏菌EC1中c-di-GMP代谢基因的序列与结构特征，为后续系统性敲除与功能解析提供了序列基础，系列梯度突变体的构建也为解析c-di-GMP浓度梯度的表型效应提供了核心实验材料。
② DGC全缺失突变体15ΔDGC（cdG0）的核心表型数据：测定了该突变体的游泳运动、群集运动、生物被膜形成能力，对水稻种子萌发的抑制率，以及水稻种子侵染定殖能力，数据来自Fig.2（A为游泳运动表型，B为群集运动表型，C为生物被膜形成定量，D为水稻种子萌发毒力测定，E为水稻种子定殖荧光显微观察）。研究意义：首次明确了c-di-GMP零水平对玉米迪基氏菌核心生理与毒力表型的影响，证实低c-di-GMP水平会显著增强细菌运动能力、完全抑制生物被膜形成，但不影响核心毒力因子的致病性，颠覆了c-di-GMP普遍调控细菌毒力的传统认知，明确了该菌中c-di-GMP的表型调控特异性。
③ 15ΔDGC突变体的功能基因转录水平与鞭毛形态数据：通过qRT-PCR测定了纤维素合成相关bcs操纵子、鞭毛转录分级调控相关基因的转录水平，通过TEM观察了细菌鞭毛数量与组装形态，数据来自Fig.3（A为qRT-PCR定量结果，B为细菌鞭毛TEM图，C为鞭毛数量统计）。研究意义：从转录水平和细胞形态层面，揭示了c-di-GMP零水平下运动增强、生物被膜丧失的分子机制，证实c-di-GMP通过调控鞭毛合成相关基因的转录、鞭毛的组装，以及纤维素合成相关基因的表达，实现对运动与生物被膜表型的调控，为表型变化提供了直接的分子与形态学证据。
④ PDE全缺失突变体7ΔPDE的核心表型数据：测定了该突变体的游泳运动、群集运动、生物被膜形成能力，水稻种子萌发抑制率，以及水稻种子侵染定殖能力，数据来自Fig.4（A为游泳运动表型，B为群集运动表型，C为生物被膜形成定量，D为水稻种子萌发毒力测定，E为水稻种子定殖荧光显微观察）。研究意义：明确了高c-di-GMP水平对细菌表型的调控效应，证实高c-di-GMP水平会完全抑制细菌运动能力、增强生物被膜形成，同时显著减弱细菌对水稻种子的侵染能力与致病性，与cdG0突变体形成表型互补，完整揭示了c-di-GMP全局库的双向调控效应。
⑤ 7ΔPDE突变体的功能基因转录水平与鞭毛形态数据：通过qRT-PCR测定了bcs操纵子、鞭毛合成相关基因的转录水平，通过TEM观察了细菌鞭毛的组装情况，数据来自Fig.5（A为qRT-PCR定量结果，B为细菌鞭毛TEM图，C为鞭毛数量统计）。研究意义：揭示了高c-di-GMP水平下运动丧失、生物被膜增强的分子机制，证实高c-di-GMP会显著抑制鞭毛合成相关基因的转录与鞭毛组装，明确了该菌中c-di-GMP通过转录水平调控鞭毛合成，进而影响运动能力与致病性，这与肠杆菌科其他细菌中c-di-GMP仅调控鞭毛旋转的机制存在显著差异，是该菌的物种特异性特征。
⑥ 系列梯度缺失突变体的游泳运动与胞内c-di-GMP浓度数据：测定了7个PDE梯度缺失突变体、15个DGC梯度缺失突变体的相对游泳运动能力，通过LC-MS测定了各突变体的胞内c-di-GMP相对浓度，数据来自Fig.6（A为PDE突变体游泳运动定量，B为PDE突变体c-di-GMP浓度定量，C为DGC突变体游泳运动定量，D为DGC突变体c-di-GMP浓度定量）。研究意义：明确了单个基因缺失对全局c-di-GMP水平与运动表型的贡献，初步鉴定出PDE14950、PDE10355和DGC14945是调控c-di-GMP水平与运动表型的核心基因，同时为拟合c-di-GMP浓度与运动能力的剂量-效应关系提供了核心定量数据。
⑦ c-di-GMP浓度与游泳运动能力的回归曲线数据：选取7个代表性菌株，拟合了c-di-GMP相对浓度与相对游泳运动能力的关系曲线，明确了二者的L型回归关系，数据来自Fig.7（A为L型回归曲线，B为对应菌株的游泳运动表型图）。研究意义：首次揭示了玉米迪基氏菌中c-di-GMP对运动表型的独特L型调控模式，明确了不同c-di-GMP浓度区间的表型效应，与其他细菌的调控模式形成鲜明对比，丰富了细菌中c-di-GMP信号调控的多样性认知，是本研究的核心创新性发现之一。
⑧ 单基因缺失突变体的游泳运动与c-di-GMP浓度数据：测定了所有单个PDE、DGC基因缺失突变体的相对游泳运动能力与胞内c-di-GMP相对浓度，数据来自Fig.8（A为单个PDE突变体游泳运动定量，B为单个DGC突变体游泳运动定量，C为单个PDE突变体c-di-GMP浓度定量，D为单个DGC突变体c-di-GMP浓度定量）。研究意义：在单基因水平验证了核心调控基因的功能，证实只有DGC14945的单基因缺失会导致c-di-GMP水平显著下降、运动能力增强，PDE14950和PDE10355的单基因缺失会导致c-di-GMP水平显著上升、运动能力下降，进一步明确了这三个基因在全局c-di-GMP调控中的主导地位。
⑨ 核心DGC、PDE蛋白的体外酶活数据：通过HPLC测定了纯化的PDE10355、PDE14950的磷酸二酯酶活性，以及DGC14945的双鸟苷酸环化酶活性，数据来自Fig.9（A-D为PDE酶活HPLC分析，E-H为DGC酶活HPLC分析）。研究意义：在体外生化水平直接验证了核心蛋白的催化活性，证实PDE14950可降解c-di-GMP为pGpG，DGC14945可利用GTP合成c-di-GMP，为这两个蛋白的生物学功能提供了最直接的生化证据，明确了其作为c-di-GMP代谢酶的催化功能。
⑩ c-di-GMP受体蛋白的结构域与功能验证数据：分析了YcgR同源物W909_08750的结构域组成，通过在7ΔPDE背景下敲除该基因，测定了突变体的游泳运动能力，数据来自Fig.10（A为蛋白结构域结构，B、C为游泳运动能力测定）。研究意义：鉴定出玉米迪基氏菌中c-di-GMP调控运动的下游受体蛋白，证实高c-di-GMP水平通过该YcgR同源蛋白抑制细菌运动，揭示了c-di-GMP调控运动表型的下游信号通路，完善了该菌中c-di-GMP的信号传导机制。
⑪ 细菌毒力因子产生能力数据：测定了DGC全缺失、PDE全缺失突变体的胞外降解酶与玉米素产生能力，数据来自Fig.S2。研究意义：证实DGC和PDE全缺失均不影响细菌核心毒力因子的产生，明确了c-di-GMP并非通过调控毒力因子合成影响致病性，为其通过运动能力调控毒力提供了关键佐证。
⑫ 鞭毛与纤维素合成基因单突变体的毒力数据：测定了ΔfliG（鞭毛缺失）、ΔbcsA（纤维素合成缺失）突变体的水稻种子侵染能力与致病性，数据来自Fig.S3。研究意义：直接证明了鞭毛介导的运动能力是该菌侵染宿主的关键，而生物被膜相关的纤维素合成不影响侵染能力，完整证实了c-di-GMP主要通过调控运动能力影响致病性。
8.研究结论
① 玉米迪基氏菌EC1中，c-di-GMP全局库是调控细菌运动-固着表型转换的核心信号：c-di-GMP零水平（DGC全缺失突变体cdG0）会导致细菌鞭毛合成相关基因转录显著上调、鞭毛数量增多，游泳/群集运动能力大幅增强，同时纤维素合成相关基因转录下调，生物被膜形成能力完全丧失；而高c-di-GMP水平（PDE全缺失突变体7ΔPDE）则完全抑制鞭毛合成与细菌运动能力，显著增强生物被膜形成。
② c-di-GMP对玉米迪基氏菌致病性的调控，主要通过调控鞭毛介导的运动能力实现，而非调控核心毒力因子的合成：DGC全缺失和PDE全缺失均不影响细菌胞外降解酶、 phytotoxin zeamine等毒力因子的产生；PDE全缺失导致的致病性减弱，是因为运动能力丧失使细菌无法有效侵染水稻种子，而鞭毛结构基因fliG的单基因缺失可完全模拟这一毒力缺陷，证实运动能力是该菌侵染宿主的关键毒力相关表型。
③ 玉米迪基氏菌中，c-di-GMP浓度与细菌游泳运动能力呈独特的L型回归曲线调控关系：当c-di-GMP浓度降至野生型的25%以下时，细菌游泳运动能力提升至野生型的4倍以上；当c-di-GMP浓度与野生型相当时，运动能力恢复至野生型水平；当c-di-GMP浓度超过野生型的2倍后，运动能力随浓度升高显著下降，直至完全丧失。这一调控模式与新月柄杆菌、沙门氏菌等已报道的细菌存在显著差异，体现了c-di-GMP调控回路的物种特异性。
④ 玉米迪基氏菌EC1的19个推定c-di-GMP信号蛋白中，仅有3个发挥全局调控的主导作用：DGC14945是主要的c-di-GMP合成酶，其单基因缺失即可导致胞内c-di-GMP水平显著下降、运动能力增强；PDE14950和PDE10355是主要的c-di-GMP降解酶，二者的单基因缺失会导致c-di-GMP水平显著上升、运动能力下降；其余DGC/PDE蛋白在实验室培养条件下对全局c-di-GMP水平无显著影响，可能在特定环境信号下发挥功能。
⑤ 玉米迪基氏菌中，高c-di-GMP水平通过PilZ结构域蛋白YcgR同源物（W909_08750）实现对细菌运动的抑制，该蛋白是c-di-GMP调控运动表型的关键下游受体。
⑥ 本研究系统揭示了玉米迪基氏菌中c-di-GMP信号系统的全局调控机制，凸显了不同细菌物种中c-di-GMP调控网络的复杂性与可塑性，为该植物病原菌的绿色防控提供了全新的分子靶点，同时为其他细菌中复杂c-di-GMP调控机制的研究提供了重要的方法学参考。
9.芬兰Bioscreen仪器测量的微生物生长曲线数据的研究意义详细解读
本研究中使用的是芬兰OY Growth Curves AB公司研发的Bioscreen-C全自动微生物生长曲线分析仪，该仪器是微生物生长动力学研究的金标准设备，可实现高通量、全自动、无干扰的细菌生长动态监测，通过连续测定培养体系的OD600值，精准呈现细菌在不同培养条件下的延滞期、指数期、稳定期等全生长周期特征。论文在材料与方法中明确使用该仪器完成了野生型与各突变菌株的生长动力学分析，其核心研究意义分为以下层面：
第一，为所有表型实验提供了基础的生长活力对照，排除了“运动/生物被膜/致病性表型变化是由细菌生长活力缺陷导致”的可能性。c-di-GMP是调控细菌多种生理过程的核心第二信使，DGC或PDE基因的全缺失可能会影响细菌的基础生长活力，进而导致表型出现非特异性变化。通过Bioscreen-C仪器全自动、高频率的生长曲线监测，可精准对比野生型EC1、DGC全缺失突变体、PDE全缺失突变体，以及各单基因缺失突变体的生长速率、最大生物量、生长周期等核心生长参数。若各菌株的生长曲线无显著差异，即可证实后续观察到的运动、生物被膜、致病性等表型变化，是c-di-GMP信号调控的特异性结果，而非细菌基础生长活力受损导致的非特异性效应，为所有表型实验的结论可靠性提供了最基础的对照与前提保障。
第二，直接验证了c-di-GMP全局水平的变化不影响玉米迪基氏菌的基础生长繁殖能力，为“c-di-GMP特异性调控运动、生物被膜与侵染相关表型，而非核心代谢”提供了关键证据。论文中发现DGC全缺失和PDE全缺失均不影响细菌胞外酶、玉米素等毒力因子的产生，而Bioscreen-C测得的生长曲线数据则进一步证实，胞内c-di-GMP水平的极端变化（零水平或极高水平），不会影响细菌在实验室培养条件下的基础生长与繁殖能力。这一结果明确了c-di-GMP在该菌中的调控特异性：其主要功能是调控细菌的表型转换（运动-固着）、宿主侵染等环境适应性行为，而非调控细菌的核心碳氮代谢与基础繁殖，为理解该菌中c-di-GMP的生物学功能提供了关键的基础认知。
第三，为整个研究的所有实验提供了标准化的菌体制备依据，保证了实验组间初始条件的一致性。微生物实验中，细菌的生长时期（对数期、稳定期）会显著影响胞内酶活、c-di-GMP浓度、基因转录水平、宿主侵染能力等实验结果。Bioscreen-C仪器绘制的精准生长曲线，可明确野生型与各突变体的生长周期特征，确定其对数生长期的时间节点、OD600对应的活菌数，为后续的胞内c-di-GMP定量、体外酶活实验、水稻种子侵染实验、qRT-PCR取样、TEM制样等所有实验，提供了标准化的菌体收集时间与初始浓度标准，确保了所有实验组的初始生理状态完全一致，消除了生长时期差异带来的系统实验误差，大幅提升了整个研究中所有实验数据的可靠性、重复性与组间可比性。
第四，为后续探究c-di-GMP在不同环境下的调控功能提供了成熟的技术平台与基础方法。Bioscreen-C仪器支持高通量、多条件的平行生长测定，可快速检测不同培养基、pH、温度、渗透压、植物根际分泌物等条件下，野生型与突变体的生长差异。本研究中通过该仪器完成了基础培养条件下的生长曲线测定，建立了该菌生长动力学的标准化检测方法，为后续研究c-di-GMP信号系统在根际环境、宿主互作、胁迫响应等田间场景下的功能，提供了可直接复用的技术方法与基础生长数据参照，大幅拓展了研究的可延伸性。
第五，与其他表型实验形成完整的证据链，让研究结论的逻辑更严谨。论文中通过运动、生物被膜、毒力等实验，证实了c-di-GMP对这些表型的调控作用，而Bioscreen-C测得的生长曲线数据，则从基础生长层面排除了所有非特异性干扰因素，形成了“c-di-GMP水平变化→不影响基础生长繁殖→特异性调控运动/生物被膜/侵染表型→通过鞭毛运动能力调控宿主致病性”的完整证据链，避免了因基础生长差异导致的结论偏差，是整个研究结论可靠性的重要底层支撑。