《Frontiers in Microbiology》：Proteomic adaptation of the hypervirulent pandemic Vibrio parahaemolyticus ST36 lineage to salinity and temperature stress
编辑推荐：
摘要
副溶血性弧菌（Vibrio parahaemolyticus）是一种嗜盐革兰氏阴性菌，广泛分布于海洋和河口生态系统中，是全球范围内海产品源性胃肠炎的主要病原体。在其众多遗传谱系中，序列型36（ST36）已成为具有显著环境持久性和全球传播能力的大流行克隆。
摘要
副溶血性弧菌（Vibrio parahaemolyticus）是一种嗜盐革兰氏阴性菌，广泛分布于海洋和河口生态系统中，是全球范围内海产品源性胃肠炎的主要病原体。在其众多遗传谱系中，序列型36（ST36）已成为具有显著环境持久性和全球传播能力的大流行克隆。其致病性与环境条件密切相关，尤其是盐度和温度，二者可调节毒力因子、应激反应蛋白及代谢通路的表达。本研究采用无标记定量（LFQ）蛋白质组学方法，研究了副溶血性弧菌在受控盐度（1%和3%）及温度（20?°C、31?°C和37?°C）条件下的差异蛋白表达谱。研究结果表明，副溶血性弧菌蛋白质组受盐热特定组合的动态非线性调控，而非这些因素的简单叠加效应。温度-盐度交互作用驱动了大多数表达变化（2,144种蛋白），其中温度是主要驱动因素（1,302种蛋白）。蛋白质组变化在31?°C时对盐度变化响应最强，在3%盐度时对温度变化响应最强，各影响约80%的蛋白，但具体响应最终取决于单个蛋白的特性。在功能上，蛋白变化体现了从温和胁迫条件下相对平衡的细胞代谢，向高盐度和高温条件下呼吸重塑、中心碳代谢重组、ATP生成及氧化还原防御的转变。此外，毒力、毒素、运动性、趋化性和防御相关蛋白分布于多个趋势聚类中。总体而言，副溶血性弧菌响应温度和盐度波动的蛋白质组模式揭示了毒力、应激反应和代谢蛋白的协同调控，这支撑了其生态适应性。除生物学见解外，本研究所采用的蛋白质组学工作流程为研究细菌胁迫适应提供了稳健的方法学框架。最终，这些发现加深了研究人员对气候变化下病原体动态的理解，并有助于改善海产品安全与风险管理。未来研究有必要对已鉴定的分子通路进行实验验证，并确认其在适应中的功能作用。
论文解读：副溶血性弧菌ST36谱系响应盐度与温度胁迫的蛋白质组适应性研究
一、研究背景与科学问题
副溶血性弧菌（Vibrio parahaemolyticus）是一种广泛分布于海洋和河口生态系统的嗜盐革兰氏阴性菌，是全球海产品相关急性细菌性胃肠炎的主要病原体之一。近年来，亚洲、美洲和欧洲沿海地区由该菌引起的疫情暴发呈上升趋势，这与海产品贸易全球化和海洋升温密切相关。在众多遗传谱系中，序列型36（ST36）已成为最具临床意义的大流行克隆之一，其最初与北美太平洋西北地区相关，近二十年来地理范围持续扩展，并在北美和欧洲引发了多起海产品相关疫情。ST36谱系具有显著的环境适应能力、多样的海洋生态系统定殖能力以及增强的毒力特征，引起了公共卫生和食品安全领域的高度关注。越来越多的证据表明，气候驱动的海水温度和盐度变化可能促进该谱系的持续存在和地理扩张，因此阐明其环境适应的分子机制至关重要。
以往研究多聚焦于单一胁迫因子、转录组分析或宿主相关条件，缺乏对盐度和温度联合效应下蛋白质组水平的系统表征。本研究首次在环境相关条件下，对副溶血性弧菌ST36谱系响应盐度和温度联合胁迫的蛋白质组进行了整合性表征，旨在揭示温度-盐度交互作用驱动的蛋白质组协同重塑规律，为理解该病原体的环境适应性和致病潜力提供分子层面的证据。
二、研究设计与核心发现
研究人员对一株分离自西班牙西南部瓜达尔基维尔河口沉积物样本的环境ST36菌株，设置了两种盐度（1%和3% NaCl）和三种温度（20?°C、31?°C和37?°C）共六种实验条件，每种条件设置三个独立生物学重复，共18个样本，采用无标记定量（LFQ）蛋白质组学方法进行分析。研究结果表明，温度是影响蛋白质组变化的主导因素，单独影响2,163种蛋白；盐度单独影响1,827种蛋白；而温度-盐度交互作用驱动的蛋白变化达2,144种，占全部变化的82%以上，其中温度是1,302种蛋白变化的主要驱动因素，842种蛋白对盐度更为敏感。盐度效应在31?°C时最强，影响80.0%的蛋白；温度效应在3%盐度时最为显著。主成分分析（PCA）显示不同盐热组合形成各自独特的蛋白质组谱，而非两因素效应的简单叠加。蛋白质组变化呈现非线性特征，不同蛋白聚类表现出截然不同的表达轨迹。
三、关键技术方法
本研究采用的主要技术方法包括：基于无标记定量（LFQ）的蛋白质组学分析流程，具体涉及蛋白质提取与酶解、液相色谱-串联质谱（LC–MS/MS）分析（采用timsTOF Pro质谱仪结合diaPASEF数据非依赖采集模式）、DIA-NN软件进行蛋白质鉴定与定量，以及基于R和Python（Julius AI在线软件）的统计分析与功能富集分析。统计分析采用双因素方差分析（two-way ANOVA）评估盐度、温度主效应及二者交互作用，以Benjamini–Hochberg错误发现率（FDR）校正多重检验，差异表达蛋白（DEP）判定阈值为名义p≤0.01且|log?FC|≥1。样本来源为西班牙瓜达尔基维尔河口Bajo de Guia区域沉积物中分离的环境ST36菌株。
四、研究结果
4.1 副溶血性弧菌对盐度和温度胁迫的差异蛋白质组响应
研究共鉴定出2,620种蛋白（在至少70%样本中被检测到）。温度是主要驱动因素，影响2,163种蛋白；盐度单独影响1,827种蛋白；温度-盐度交互作用影响2,144种蛋白。盐度效应在31?°C时最强，影响80.0%的蛋白；温度效应在3%盐度时最显著。PCA分析显示不同盐热组合形成独特的蛋白质组谱，各条件生物学重复具有良好的重现性。
4.2 蛋白质组表达趋势聚类
蛋白被分为六个聚类。Cluster 5（n=534）在3%盐度/37?°C条件下表达升高；Cluster 1（n=453）在1%盐度下相对丰富，但在3%盐度/37?°C下显著降低；Cluster 4（n=313）和Cluster 2（n=369）均在3%盐度/31?°C达到峰值后于37?°C急剧下降；Cluster 3（n=440）在1%盐度下随温度升高而增加，但在3%盐度/31?°C显著降低；Cluster 6（n=511）呈现复杂轨迹。这些结果表明蛋白质组呈非线性动态调控。
4.3 差异蛋白亚细胞分布
大多数差异蛋白定位于细胞质（36.8%）或其他细胞区室（51.6%），周质蛋白占4.9%，外膜蛋白占4.2%，内膜蛋白占1.9%，胞外蛋白占0.6%。Cluster 4具有最强的膜相关特征，与其富集呼吸和氧化磷酸化功能一致。仅18–20%的蛋白具有数据库支持的亚细胞定位注释。
4.4 代谢和核心细胞维持程序重编程
核心细胞维持蛋白（翻译/核糖体相关蛋白、RNA加工因子、DNA复制/修复蛋白）表现出显著的条件特异性变化。能量代谢相关蛋白（糖酵解、三羧酸循环（TCA）、氧化磷酸化、电子传递、ATP合成）在3%盐度/31?°C达到最高平均相对丰度，随后在3%盐度/37?°C下降。Cluster 4包含18个氧化磷酸化/电子传递候选蛋白，KEGG分析确认氧化磷酸化在该聚类中显著富集（FDR=0.0168）。3%盐度/31?°C代表活跃的呼吸重塑状态，3%盐度/37?°C则反映更受限的应激生存状态。
4.5 铁获取、储存和硫酸盐/亚硫酸盐同化响应
铁相关蛋白呈现异质性响应模式。铁（III）受体VPA1435随温度升高而增加，在3%盐度/37?°C达峰值；铁（III）ABC转运蛋白周质结合蛋白VPA0657和VP249在3%盐度/31?°C低丰度，但在3%盐度/37?°C强烈上调。相比之下，硫酸盐/亚硫酸盐同化蛋白（CysJ、CysI、CysH、CysN、CysD、CysC）表现出高度协同的响应，在1%盐度（20–31?°C）高丰度，但在3%盐度/37?°C被强烈抑制。TonB能量转导和转运组分在温暖条件下上调，而TonB依赖性受体在3%盐度/31?°C下降。
4.6 盐度和温度变化下的运动性改变
结构鞭毛蛋白（flaA、flaB/flaD、flaF、flaG）在3%盐度/20?°C达峰值，但在3%盐度/31?°C显著下降，37?°C部分恢复。鞭毛调控蛋白（VPA0262、FlhF、FlhG、FliM）在1%盐度/37?°C强烈上调。趋化蛋白（CheA、CheX、甲基趋化受体蛋白）在1%盐度/20?°C和3%盐度/31?°C高丰度，但在3%盐度/37?°C下降。结构鞭毛蛋白与其调控因子之间的表达差异表明鞭毛组装和调控在环境波动下并非同步。
4.7 温度适应和应激响应
冷休克蛋白（CspD、CspA、VPA0552）在1%盐度/31?°C达峰值，在所有3%盐度条件下均降低。大多数热休克蛋白（16种中9种，包括ClpB、DnaK、DnaJ、HtpG、GroEL1、GroEL2、GroES1、GroES2、GrpE）在3%盐度/31?°C强烈增加并保持相对升高。噬菌体休克蛋白（PspA–E）呈现非均匀的应激响应模式，PspA和PspE在3%盐度/37?°C上调最强。
4.8 毒力调控和致病因子
ToxR在1%盐度/37?°C显著上调，3%盐度抑制了这种升高。RTX毒素激活子在两种盐度下均从20?°C到31?°C增加，但在3%盐度/37?°C下降。多种溶血素异构体在31?°C达到最高水平，在37?°C下降（主要在3%盐度）。第二组溶血素异构体在1%盐度下37?°C表达增加更为显著。毒力相关蛋白并非在同一条件下达到峰值，表明存在不同的调控通路。
4.9 抗氧化防御、群体感应和分子伴侣活性
抗氧化和氧化还原防御蛋白分为四个条件依赖性趋势组：第一组（谷胱甘肽过氧化物酶、细胞色素c过氧化物酶、超氧化物歧化酶（SOD）、谷胱甘肽S-转移酶（GST）、DsbA样硫氧还蛋白）仅在1%盐度/37?°C增加；第二组（过氧化氢酶、AhpF、谷胱甘肽还原酶、谷氧还蛋白相关蛋白、过氧化物氧还蛋白）在3%盐度/31?°C急剧达峰后下降；第三组（硫氧还蛋白、谷氧还蛋白、谷胱甘肽、二硫键异构化相关蛋白）在1%盐度/31?°C高丰度，3%盐度/31?°C显著下降后恢复；第四组（AhpC、KatG1、KatG2、硫氧还蛋白过氧化物酶/还原酶、SOD、GST、KefG）在高盐度下渐进增加，在3%盐度/37?°C达峰值。群体感应相关蛋白LuxP随温度升高而下降，在高盐度下下调更为显著。
五、讨论与结论
本研究表明，副溶血性弧菌ST36对盐度和温度并发变化呈现复杂、非线性和多因素的蛋白质组响应。在相对温和条件下（低盐度），蛋白聚类以核心细胞维持功能为主；随盐度和温度升高（3%盐度/31?°C），蛋白质组转向活跃的代谢重塑状态，呼吸蛋白、氧化磷酸化组分、中心碳代谢酶、ATP合酶亚基、分子伴侣、热响应蛋白、抗氧化防御和渗透应激候选蛋白增加；在最高胁迫条件（3%盐度/37?°C）下呈现更受限的应激生存特征，核心维持和硫酸盐同化蛋白减少，抗氧化、渗透应激和包膜保护响应增强。结构鞭毛蛋白与调控因子的表达解耦表明调控信号并非线性决定鞭毛组装。毒力因子（ToxR、RTX毒素、溶血素）在1%盐度/37?°C大多上调，在3%盐度下于31?°C达峰值后37?°C下调。高盐度促进分子伴侣、抗氧化蛋白和VI型分泌系统组分的过表达，同时降低冷休克蛋白表达。LuxP在高盐度和高温下的渐进下调可能代表一种有利于运动性和环境探索而非生物膜相关群体行为的适应性策略。
本研究与既往研究的主要区别在于实验设计和生物学范围：既往研究多聚焦于低盐度和胆汁盐双重胁迫下的人类肠道环境适应，或仅关注热胁迫下的代谢和生长变化，而本研究从生态学视角系统绘制了盐度和温度两个海洋生态系统核心因子间的交互作用，整合了代谢调控、群体感应、氧化应激防御和毒力相关蛋白质组重塑。研究存在两方面局限：其一，仅使用单一环境ST36分离株，观察到的蛋白质组响应不能外推至整个ST36谱系，未来需整合不同地理分离株的基因组学、转录组学和蛋白质组学多组学研究；其二，差异表达蛋白的丰度变化不能直接证明功能活性或翻译后调控，需通过表型和功能实验加以验证。
研究结论表明，温度和盐度是副溶血性弧菌ST36生理和致病性的关键环境决定因素，二者的交互作用驱动了代谢、毒力、应激反应和信号通路的协同调控。这些蛋白质组改变反映了该海洋病原体的潜在适应机制，提示盐度和温度并发变化可能直接影响其生存、环境持久性和宿主互作能力，为气候变化背景下病原体动态监测和海产品安全管理提供了重要的分子基础。
