Comparative proteomics and secretomics revealed virulence, and coresistance-related factors in non O1/O139 Vibrio cholerae recovered from 16 species of consumable aquatic animals
比较蛋白质组学与分泌组学揭示16种食用水产品来源非O1/O139群霍乱弧菌的毒力与共耐药相关因子
来源：Journal of Proteomics 251 (2022) 104408
1.论文摘要核心内容
霍乱弧菌是可引发人类霍乱大流行的革兰氏阴性菌，广泛分布于全球水生环境中，鉴定水产品中霍乱弧菌的风险因子对保障食品安全至关重要。本研究对分离自16种食用水产品的20株非O1/O139群霍乱弧菌，完成了毒力相关基因检测、抗生素与重金属交叉耐药表型分析及基因组指纹图谱分析；借助二维凝胶电泳（2D-GE）和液相色谱-串联质谱（LC-MS/MS）技术，获得了不同基因型和表型菌株的分泌组与蛋白质组图谱。比较分泌组学分析发现，20株菌株间存在4个共同胞外蛋白和45个差异胞外蛋白，其中包含13个毒力相关蛋白、8个耐药相关蛋白；共鉴定到21972个胞内蛋白，比较蛋白质组学分析显示菌株间存在215个共同胞内蛋白和913个差异胞内蛋白，其中包含22个毒力相关蛋白、8个耐药相关蛋白。此外，鱼类和贝类来源的霍乱弧菌菌株呈现出显著不同的分泌组与蛋白质组特征，同时还在菌株中发现了大量功能未知的新蛋白和菌株特异性蛋白。本研究首次解析了水产品来源非O1/O139群霍乱弧菌的蛋白质组与分泌组，鉴定了一批毒力与共耐药相关因子，填补了该类菌株致病与耐药机制的研究空白，也为全球范围内这一重要水传播病原体的新型诊断靶点与疫苗研发提供了理论支撑。
2.关键词（中文）
霍乱弧菌、分泌组、蛋白质组、毒力、耐药性、水生动物
3.研究目的
① 明确16种食用水产品来源非O1/O139群霍乱弧菌的遗传多样性、毒力基因谱、抗生素耐药与重金属耐受表型，评估其食品安全与公共卫生风险；
② 首次系统解析水产品来源非O1/O139群霍乱弧菌的分泌组与全蛋白质组图谱，鉴定菌株间的共同蛋白与差异蛋白，填补该类菌株蛋白组学研究的空白；
③ 筛选并鉴定非O1/O139群霍乱弧菌中潜在的毒力相关蛋白、抗生素与重金属共耐药相关蛋白，揭示其不携带霍乱核心毒素却仍具致病潜力的分子基础；
④ 对比鱼类与贝类来源菌株的分泌组、蛋白质组差异，解析霍乱弧菌对不同水生宿主的适应性进化特征；
⑤ 发掘菌株特异性蛋白与功能未知的新型蛋白，为霍乱弧菌的新型检测靶点、抑菌剂研发及疫苗设计提供候选分子。
4.研究思路
本研究以16种常见食用水产品中分离的20株非O1/O139群霍乱弧菌为研究对象，围绕“风险表型鉴定-遗传背景解析-蛋白组全景分析-功能因子挖掘”的核心逻辑开展系统性研究，具体研究思路如下：
第一步，样本采集与菌株分离鉴定：从上海两大水产市场采集16种食用水产品（12种贝类、4种鱼类），通过选择性培养、生化试验、特异性PCR（lolB、16S rRNA基因）完成霍乱弧菌的分离与菌种鉴定；
第二步，菌株风险表型系统分析：通过PCR检测菌株的核心毒力基因（ctxAB、tcpA）与附属毒力相关基因，通过药敏试验测定10种抗生素的耐药表型，通过耐受性试验测定8种重金属的耐受表型，明确菌株的毒力潜力与共耐药特征；
第三步，菌株遗传多样性解析：通过ERIC-PCR基因组指纹图谱完成菌株基因分型与聚类分析，通过Illumina平台完成19株菌株的全基因组草图测序与组装注释，基于同源基因序列构建系统发育树，明确菌株的遗传进化关系与多样性特征；
第四步，分泌组与蛋白质组学分析：通过Bioscreen全自动生长曲线分析仪测定菌株生长曲线，收集对数生长期菌株的胞外分泌蛋白与胞内总蛋白；优化2D-GE电泳条件完成菌株分泌组的分离与图谱构建，通过LC-MS/MS技术完成2D-GE差异蛋白点的鉴定，同时通过LC-MS/MS完成菌株全胞内蛋白的定性定量分析；
第五步，功能蛋白系统筛选与验证：对鉴定的胞外与胞内蛋白进行功能注释与GO分类，筛选毒力相关、抗生素与重金属耐药相关蛋白，对比鱼类与贝类来源菌株的蛋白组差异，发掘菌株特异性蛋白与功能未知新蛋白；通过qRT-PCR技术在转录水平验证代表性蛋白编码基因的表达水平，验证蛋白组学结果的可靠性；
最后，整合菌株的表型、遗传与蛋白组学数据，明确水产品来源非O1/O139群霍乱弧菌的致病与耐药相关核心因子，评估其食品安全风险，形成完整的研究结论。
5.研究亮点
① 首次完成了16种食用水产品来源非O1/O139群霍乱弧菌的分泌组与全蛋白质组系统性解析，填补了该类环境源菌株蛋白组学研究的空白，突破了过往研究多聚焦于临床O1/O139流行株的局限；
② 明确了我国市售水产品中非O1/O139群霍乱弧菌的风险特征：菌株均不携带霍乱大流行相关的ctxAB、tcpA等核心毒素基因，但普遍携带rtxABCD、hlyA、tlh、hapA等附属毒力相关基因，76.2%的菌株呈现多重耐药表型，同时存在普遍的重金属耐受性，揭示了该类菌株的潜在公共卫生风险；
③ 系统鉴定了该类菌株中35个毒力相关蛋白（13个分泌型、22个胞内型）、16个耐药相关蛋白（8个分泌型、8个胞内型），阐明了其不携带核心霍乱毒素却仍具致病潜力、以及抗生素-重金属共耐药的分子基础，为该菌的防控提供了全新的分子靶点；
④ 首次发现鱼类与贝类来源的霍乱弧菌菌株存在显著的分泌组与蛋白质组特异性差异，鉴定了宿主来源特异性的蛋白标志物，为解析霍乱弧菌在不同水生宿主中的适应性进化机制提供了关键的蛋白质组学证据；
⑤ 发掘了大量菌株特异性蛋白与181个功能未知的新型蛋白，拓展了霍乱弧菌的基因功能研究边界，同时建立了一套适用于水产品源弧菌的分泌组-蛋白质组联合分析技术体系，为其他食源性致病弧菌的风险因子挖掘提供了可复制的研究范式；
⑥ 结合基因组、表型组与蛋白组多组学数据，完整揭示了水产品源非O1/O139群霍乱弧菌的遗传多样性、表型异质性与分子特征，为水产品中霍乱弧菌的风险监测、检测方法开发与安全防控提供了全面的理论支撑。
6.可延伸的方向
① 对本研究发现的功能未知新型蛋白、菌株特异性蛋白开展体内外功能验证，通过基因敲除与回补试验明确其在霍乱弧菌致病、耐药、环境适应中的生物学功能，完善非O1/O139群霍乱弧菌的毒力调控网络；
② 深入解析鱼类与贝类来源菌株蛋白组差异的分子机制，结合比较基因组学分析，揭示霍乱弧菌对不同水生宿主的适应性进化规律，明确其在水产品中定殖、存活的关键分子机制；
③ 探究该类菌株抗生素与重金属共耐药的分子机制，明确耐药相关蛋白的功能，解析共选择效应的遗传基础，评估水产养殖环境中抗生素与重金属滥用对霍乱弧菌耐药性进化的驱动作用；
④ 基于鉴定的毒力与耐药相关核心蛋白，筛选特异性的检测靶点，开发水产品中非O1/O139群霍乱弧菌的快速检测方法；同时筛选靶向关键毒力/耐药蛋白的小分子抑制剂，研发新型抑菌剂用于水产品中该菌的防控；
⑤ 通过细胞模型与动物感染模型，验证本研究中分离菌株的致病力，明确不同毒力相关蛋白在菌株感染宿主过程中的具体作用，量化评估水产品源非O1/O139群霍乱弧菌的人体健康风险；
⑥ 探究该类菌株在水产品冷链储存、加工、杀菌等环节中的蛋白质组表达变化，明确其环境胁迫响应的关键蛋白与调控机制，为水产品加工过程中该菌的灭活控制提供理论依据；
⑦ 开展非O1/O139群与O1/O139群霍乱弧菌的比较蛋白质组学分析，解析两类菌株的毒力进化规律，揭示环境源非O1/O139群菌株向致病流行株进化的潜在风险与分子特征。
7.测量的数据、研究意义及对应图表
① 20株霍乱弧菌的毒力基因谱、抗生素耐药表型、重金属耐受表型数据：包含菌株的宿主来源、10种毒力相关基因携带情况、10种抗生素的耐药谱、8种重金属的耐受谱，数据来自Table 1。研究意义：明确了水产品来源非O1/O139群霍乱弧菌的核心毒力基因特征，证实菌株均不携带霍乱大流行相关的ctxAB、tcpA等核心毒素基因，但普遍携带rtxABCD、hlyA等附属毒力相关基因；同时揭示了菌株高多重耐药率与重金属普遍耐受的表型特征，为菌株的食品安全风险评估提供了最基础的表型数据。
② 菌株的基因组指纹图谱与系统发育分析数据：包含20株菌株的ERIC-PCR基因组指纹图谱、UPGMA聚类结果，以及基于19株菌株基因组草图与参考菌株序列构建的系统发育树，数据来自Fig.1（A为ERIC-PCR基因组指纹图谱，B为全基因组同源基因构建的系统发育树）。研究意义：证实了水产品来源非O1/O139群霍乱弧菌具有丰富的遗传多样性，20株菌株分为20个ERIC基因型与8个聚类簇；系统发育分析显示所有水产品源菌株均聚为独立分支，与临床流行株遗传距离较远，明确了养殖环境对菌株遗传进化的塑造作用，为菌株的溯源分析提供了关键的遗传背景依据。
③ 菌株的分泌组2D-GE图谱数据：包含20株霍乱弧菌与标准菌株GIM1.449的胞外蛋白二维凝胶电泳图谱、各菌株分泌蛋白数量统计、共同蛋白与差异蛋白的分布定位，数据来自Fig.2（A、B为20株菌株的2D-GE分泌组图谱）。研究意义：首次获得了16种水产品来源非O1/O139群霍乱弧菌的分泌组全景图谱，明确了不同菌株胞外蛋白分泌的显著异质性，鉴定出4个所有菌株共有的核心分泌蛋白与45个差异分泌蛋白，为后续毒力、耐药相关分泌蛋白的鉴定提供了基础图谱与筛选依据。
④ 分泌组与蛋白质组中毒力相关蛋白的鉴定数据：包含通过LC-MS/MS鉴定的13个胞外毒力相关蛋白、22个胞内毒力相关蛋白的参考序列编号、基因名、分子量、等电点、序列覆盖度与功能注释，数据来自Table 2。研究意义：系统筛选并鉴定了非O1/O139群霍乱弧菌的关键毒力相关蛋白，涵盖黏附侵袭、鞭毛运动、免疫逃逸、毒素分泌、群体感应调控等多个致病相关过程，揭示了该类菌株不携带核心霍乱毒素却仍具潜在致病性的分子基础，填补了该类菌株分泌型与胞内毒力因子的研究空白。
⑤ 分泌组与蛋白质组中耐药相关蛋白的鉴定数据：包含通过LC-MS/MS鉴定的8个胞外耐药相关蛋白、14个胞内耐药相关蛋白的参考序列编号、基因名、分子量、等电点、序列覆盖度与功能注释，数据来自Table 3。研究意义：明确了该类菌株中与抗生素耐药、重金属耐受相关的核心功能蛋白，包括ABC转运蛋白、β-内酰胺酶、金属离子结合与调控蛋白、外排系统相关蛋白等，揭示了菌株多重耐药与抗生素-重金属共耐受的分子机制，为耐药性防控提供了潜在的分子靶点。
⑥ 菌株共有胞内蛋白的GO功能分类数据：包含215个所有菌株共有的胞内蛋白在生物过程、细胞组分、分子功能三大GO类别的功能注释与丰度分布，数据来自Fig.3。研究意义：明确了非O1/O139群霍乱弧菌核心生命活动相关的保守蛋白功能谱，证实代谢过程、细胞过程、催化活性是菌株最核心的保守生物学功能，为理解该菌的基础生理特征、核心代谢通路提供了蛋白质组层面的系统依据。
⑦ 菌株差异胞内蛋白的GO功能分类数据：包含913个菌株间差异胞内蛋白的GO功能注释、各类别丰度分布，以及181个无同源匹配的功能未知新蛋白统计，数据来自Fig.4。研究意义：揭示了不同菌株蛋白表达的显著异质性，明确了差异蛋白主要参与的生物学过程与分子功能，同时发掘了大量功能未知的新型蛋白，为解析菌株的表型多样性、宿主适应性差异提供了分子依据，也拓展了霍乱弧菌的基因功能研究边界。
⑧ 菌株分泌组共同与差异胞外蛋白的详细鉴定数据：包含14个共同胞外蛋白、45个差异胞外蛋白的LC-MS/MS详细鉴定结果、序列覆盖度、功能注释，数据来自补充表Table S3（共同胞外蛋白）、Table S4（差异胞外蛋白）。研究意义：完整解析了该类菌株胞外蛋白的组成与功能，明确了核心分泌蛋白主要参与细菌代谢、物质转运、鞭毛合成等过程，为理解其环境适应、宿主互作机制提供了全面的分泌蛋白数据集。
⑨ 鱼类与贝类来源菌株特异性胞内蛋白的鉴定数据：包含鱼类来源菌株特有153个差异胞内蛋白、贝类来源菌株特有48个差异胞内蛋白的鉴定结果与功能注释，数据来自补充表Table S5。研究意义：首次揭示了不同水生宿主来源霍乱弧菌的蛋白表达特异性，鉴定了宿主来源特异性的蛋白标志物，为解析霍乱弧菌在鱼类和贝类中的定殖机制、宿主适应性进化提供了关键的蛋白质组证据。
⑩ 菌株基因组测序与组装注释数据：包含19株菌株的基因组草图基本信息，如测序深度、基因组大小、预测基因数量、蛋白编码基因数量，数据来自补充表Table S2。研究意义：获得了菌株的全基因组序列信息，为系统发育分析、功能基因注释、蛋白序列比对提供了完整的基因组参考，也为蛋白质组学的蛋白鉴定结果提供了基因组层面的支撑。
⑪ 蛋白编码基因的转录水平验证数据：通过qRT-PCR测定的代表性毒力、耐药相关蛋白编码基因的转录表达水平，数据来自补充图Fig.S1。研究意义：在转录水平验证了蛋白质组与分泌组的鉴定结果，证实了蛋白表达丰度与对应基因转录水平的一致性，为所鉴定蛋白的功能可靠性提供了转录层面的佐证。
8.研究结论
① 本研究首次系统解析了16种食用水产品来源非O1/O139群霍乱弧菌的蛋白质组与分泌组，所有菌株均不携带霍乱大流行相关的ctxAB、tcpA、ace、zot等核心毒素基因，但普遍携带高丰度的附属毒力相关基因，其中rtxABCD携带率90.5%-100%、hlyA携带率95.2%、tlh携带率90.5%、hapA携带率100%，证实该类菌株仍具备潜在致病风险。
② 水产品来源非O1/O139群霍乱弧菌呈现严重的耐药风险，76.2%的菌株具有多重耐药（MDR）表型，氨苄西林耐药率最高（76.2%），其次为大观霉素（57.1%）、链霉素与利福平（均为52.4%）；同时菌株对重金属普遍存在耐受性，其中Zn²⁺耐受率47.6%、Ni²⁺耐受率42.9%、Cu²⁺与Hg²⁺耐受率均为28.6%，证实菌株存在抗生素与重金属共耐受特征。
③ ERIC-PCR基因组指纹图谱与全基因组系统发育分析证实，20株水产品来源菌株具有丰富的遗传多样性，分为20个独立的ERIC基因型，且所有环境源菌株均聚为独立的系统发育分支，与临床O1/O139流行株遗传距离较远，表明水产养殖环境塑造了菌株独特的遗传进化特征。
④ 比较分泌组学分析显示，20株菌株间存在4个核心共同胞外蛋白与45个差异胞外蛋白，其中鉴定出13个毒力相关蛋白、8个耐药相关蛋白；比较蛋白质组学共鉴定到21972个胞内蛋白，菌株间存在215个共同胞内蛋白与913个差异胞内蛋白，其中包含22个毒力相关蛋白、8个耐药相关蛋白，这些蛋白是菌株致病与耐药表型的核心分子基础。
⑤ 鱼类与贝类来源的霍乱弧菌菌株呈现显著不同的分泌组与蛋白质组特征，鉴定出宿主来源特异性的表达蛋白，揭示了霍乱弧菌对不同水生宿主的适应性蛋白表达差异。
⑥ 研究在菌株中发现了大量菌株特异性蛋白与181个功能未知的新型蛋白，为霍乱弧菌的基因功能研究、新型检测靶点与疫苗研发提供了全新的候选分子。
⑦ 本研究结果填补了水产品来源非O1/O139群霍乱弧菌致病与耐药机制的研究空白，全面评估了该类菌株的食品安全风险，为全球范围内霍乱弧菌的防控、新型诊断方法与疫苗研发提供了重要的理论支撑与数据基础。
9.芬兰Bioscreen仪器测量的微生物生长曲线数据的研究意义详细解读
本研究中使用的Bioscreen全自动生长曲线分析仪，是芬兰OY Growth Curves AB公司研发的Bioscreen C全自动微生物生长曲线分析系统，该仪器是微生物生长动力学研究的金标准设备，可实现高通量、全自动、无干扰的细菌生长动态连续监测，论文中通过该仪器完成了霍乱弧菌菌株的生长曲线测定，其核心研究意义分为以下六大层面：
第一，为分泌组学分析提供了标准化的取样基准，是分泌蛋白提取实验可靠性的核心前提。分泌组学分析对细菌的生长时期要求极高，不同生长阶段的菌株蛋白分泌谱会发生剧烈变化，只有处于同一生长时期（对数生长期）的菌株，其分泌组数据才具备组间可比性。Bioscreen仪器可通过高频率的OD值连续监测，精准绘制每一株霍乱弧菌的生长曲线，明确其延滞期、对数生长期、稳定期的具体时间节点，从而确保所有20株菌株的胞外蛋白均在完全一致的对数生长期中期提取，彻底消除了生长时期差异带来的分泌蛋白谱系统误差，保证了不同菌株间分泌组数据的可比性、重复性与准确性，是整个分泌组学研究的基础质控核心。
第二，验证了不同菌株的基础生长活力，排除了“蛋白表达差异、表型异质性是由菌株基础生长缺陷导致”的可能性。本研究中不同菌株的耐药表型、蛋白分泌能力、毒力相关蛋白表达存在显著差异，而这些差异可能由菌株本身的生长繁殖能力缺陷导致。通过Bioscreen仪器绘制的完整生长曲线，可精准对比20株菌株的生长速率、最大生物量、代时等核心生长动力学参数，证实菌株间的表型与蛋白组差异并非源于基础生长活力的差异，而是菌株自身的基因与蛋白表达特异性调控的结果，为所有表型与蛋白组学结论的可靠性提供了最基础的阴性对照与逻辑支撑。
第三，为整个研究的所有生物学实验提供了统一的菌体制备标准化依据，确保所有实验组的初始生理状态完全一致。本研究涉及药敏试验、重金属耐受试验、ERIC-PCR、基因组测序、蛋白提取、qRT-PCR验证等大量平行实验，而细菌的初始生长状态、活菌浓度、生理时期会直接影响所有实验的结果稳定性。Bioscreen仪器绘制的生长曲线，可明确每一株菌株的生长周期特征，确定其对数生长期的OD600值与活菌数的对应关系，为所有实验提供了标准化的菌体收集时间、初始接种浓度与菌液制备标准，彻底消除了不同菌株间生理状态差异带来的实验误差，大幅提升了整个研究所有实验数据的重复性、组间可比性与结论的严谨性。
第四，为菌株的表型差异分析提供了基础生长动力学数据支撑。本研究中不同菌株的抗生素耐药性、重金属耐受性、蛋白分泌能力存在显著差异，而这些表型与菌株的生长速率、环境适应能力直接相关。Bioscreen仪器可高通量完成不同培养条件下的菌株生长曲线测定，其获得的生长动力学数据，可用于分析菌株生长特征与耐药表型、蛋白分泌能力的相关性，为解析“生长特征-耐药表型-蛋白表达”的关联机制提供了基础数据，也为后续研究菌株在水产养殖环境、水产品基质中的存活与定殖能力提供了生长动力学参考。
第五，为后续的蛋白功能验证实验提供了成熟的技术方法与基础数据参照。本研究发现了大量毒力、耐药相关蛋白与功能未知的新蛋白，后续需通过基因敲除、过表达等手段开展功能验证。Bioscreen仪器建立的霍乱弧菌生长曲线标准化测定方法，可快速对比野生型与基因突变株的生长动力学差异，明确目标基因对菌株生长、环境适应、胁迫响应的影响，是后续基因功能研究的核心表型检测技术，本研究中获得的野生型菌株生长曲线数据，也为后续功能验证实验提供了可靠的野生型参照基准。
第六，与基因组、蛋白组学数据形成完整的证据链，完善了菌株的多组学特征解析。本研究结合了基因组、表型组与蛋白组多组学数据，而Bioscreen仪器获得的生长曲线数据，是菌株表型组的核心基础数据。该数据与基因组的功能基因注释、蛋白组的代谢通路蛋白表达数据相结合，可完整解析菌株的基因型-蛋白表达-生长表型的关联机制，比如通过生长曲线与蛋白组数据联合分析，明确代谢相关差异蛋白对菌株生长速率、环境适应能力的调控作用，让整个研究的结论形成完整、闭环的逻辑证据链，大幅提升了研究结果的系统性与说服力。