机器学习(Machine-learning)引导的蛋白质组学分析揭示壳聚糖(Chitosan)调控爪哇根结线虫(Meloidogyne javanica)–厚垣轮枝菌(Pochonia chlamydosporia)互作中的寄生过程

《Physiological and Molecular Plant Pathology》：Machine-learning-guided proteomic profiling reveals chitosan-modulated parasitism in the Pochonia chlamydosporia–root-knot nematode interaction

　　

研究人员采用整合蛋白质组学方法结合差异丰度分析、序列特征聚类、系统发育比较及机器学习(Machine-learning)方法，分析了爪哇根结线虫(Meloidogyne javanica, MELJA)卵与厚垣轮枝菌(Pochonia chlamydosporia, Pc)在有无壳聚糖(Chitosan, T8)处理下接种后24 h和96 h的两个时间点(full proteome与secretome)，以解析三方体系分子机制。蛋白质组全谱分析表明Pc呈现刺激依赖性的重组，多元模型识别出早期代谢激活因子与晚期细胞壁相关酶是时间响应关键驱动因子。真菌分泌组(Secretome)由功能趋同性而非进化亲缘关系组织，胞外水解酶和CFEM蛋白主导细胞组分(Cellular Component, CC)预测。壳聚糖与线虫蛋白质组显著分化相关联，抑制卵黄蛋白原(Vitellogenin)、ZP结构域蛋白及Piwi调控因子，同时上调ATP合酶和EF1-α异构体。功能富集提示壳聚糖与线虫细胞骨架强化及染色质活性通路偏移相伴随，而单纯真菌感染诱导基因沉默和代谢应激特征。跨机器学习模型获得13个共识生物标志物(Consensus Biomarkers)。壳聚糖与增强的真菌寄生相伴随，对应Pc保守分泌模块激活及MELJA涉及抗性相关通路的蛋白质组偏移。上述发现为受控体外(in vitro)条件下鉴定线虫敏感性与真菌寄生分子生物标志物提供了机器学习引导框架。

论文解读：机器学习引导的蛋白质组学解析壳聚糖调控厚垣轮枝菌–爪哇根结线虫互作中的寄生机制

研究背景与意义

根结线虫(Meloidogyne spp.)是全球重要的植物寄生线虫(Plant-parasitic nematodes, PPNs)，对茄科及豆科作物造成严重危害，传统杀线虫剂因环境风险受限，亟需可持续防控策略。厚垣轮枝菌(Pochonia chlamydosporia, Pc)是重要的线虫卵寄生性生防真菌(Nematophagous fungus)，可通过附着胞(Appressoria)形成及分泌水解酶（如VCP1、几丁质酶、糖苷水解酶）穿透并寄生根结线虫卵，亦可定殖根际促进植物生长。近期证据表明Pc感染期可重塑细胞壁并局部合成壳聚糖(Chitosan)，外源壳聚糖对多数植物病原菌具抑制作用，但Pc因膜流动性低而对壳聚糖耐受，且壳聚糖可促进Pc碳水化合物代谢、蛋白水解、氧化还原平衡及卵寄生与根部定殖能力。然而Pc遇线虫卵并受壳聚糖调控时的早期分子特征仍不明确。本研究通过高精度蛋白质组学与机器学习框架，在受控体外(in vitro)微宇宙体系中表征Pc及爪哇根结线虫(Meloidogyne javanica, MELJA)卵的时间依赖性蛋白质组响应，以鉴定驱动真菌发育与寄生的蛋白标志物，为壳聚糖协同生防菌的生物防治制剂开发提供分子依据。该论文发表于《Physiological and Molecular Plant Pathology》。

主要关键技术方法

研究人员设置四组处理——Pc单独(对照)、Pc+MELJA卵(MELJA_Pc)、Pc+壳聚糖T8 1 mg/mL(PcT8)、Pc+MELJA卵+壳聚糖(MELJA_Pc_T8)，于接种后24 h与96 h收集样品(n=3)。总蛋白经TCA–脱氧胆酸沉淀、还原烷基化及Trypsin酶解后，进行LC–MS/MS质谱分析，检索UniProtKB中Pc123(UP00052797)及Meloidogyne属(UP000887561)数据库。定量采用蛋白检测采样率(Protein Detection Sampling Rate, PDSR)与平均检出百分比(Average Detection, AD)归一化。差异丰度(Differentially Abundant Proteins, DAPs)以log₂FC±1为阈值筛选。多维统计含UMAP、PERMANOVA/SIMPER、Jaccard/Bray–Curtis距离；监督学习采用支持向量机(Support Vector Machine, SVM)、偏最小二乘判别分析(Partial Least Squares Discriminant Analysis, PLS-DA)提取变量重要性投影(Variable Importance in Projection, VIP)及随机森林(Random Forest, RF, 1500–2000棵树)获取Mean Decrease Accuracy/Gini，跨模型取共识生物标志物；非胞质(Non-cytoplasmic, Non-Cyt)结构域蛋白做氨基酸特征向量化、PCoA及基于非胞质域与全长序列的系统发育树Mantel检验；细胞组分(Cellular Component, CC)用OmicsBox与GO.db注释，并以广义线性混合模型(Generalized Linear Mixed Model, GLMM)＋Tukey/Sidak事后检验验证标志物时间动力学。

研究结果

3.1. Early germination and parasitism stages of Pochonia chlamydosporia 123 under chitosan and nematode egg cues（壳聚糖与线虫卵信号下Pc123早期萌发与寄生阶段）

显微与萌发实验显示Pc123分生孢子在疏水聚苯乙烯表面较亲水玻璃更易萌发，MELJA卵存在可补偿亲水性表面的萌发抑制。壳聚糖(1 mg/mL)显著抑制24 h萌发，但提高卵密度可部分恢复。原生质体(Protoplast)在无细胞壁情况下仍可感知宿主、黏附并穿透卵壳，内部定植早于细胞壁再生；带完整细胞壁的分生孢子侵染速度快于原生质体。表明Pc具宿主感应程序，表面物理信号与卵源化学信号共同调控早期侵染。

3.2. Proteome-wide restructuring and multivariate identification of key drivers in Pc123（Pc123全蛋白质组重构与关键驱动因子多元识别）

UMAP显示处理间与时间点明显分离，壳聚糖处理组96 h相似度高(92.4%)。DAP在24 h以诱导为主，96 h PcRKNT8保留大稳定核心而PcT8出现较多下调。PLS-DA显示24 h高VIP蛋白为ATP合酶α亚基、线粒体外膜孔蛋白等能量相关，96 h为1,3-β-葡聚糖基转移酶等细胞壁相关。RF识别出HSP70-1等为领先驱动因子，共识分析(PLS-DA VIP top60∩RF top30)获8个共识标志物：早期能量相关ATP合酶α(A0A4Q7JSQ1)、甘油激酶(A0A4Q7KBQ4)；晚期β-葡萄糖苷酶1(A0A4Q7K9B3)、葡糖淀粉酶(A0A4Q7K7Z2)等碳水化物活性酶；以及醛脱氢酶、黄素单加氧酶、组氨酸激酶/HSP90样ATP酶。GLMM验证其时间依赖性丰度变化具显著性。

3.3. Functional proteome plasticity of Pc123 under stress（胁迫下Pc123功能蛋白质组可塑性）

GO注释显示24 hpi分子功能(Molecular Function, MF)富集结构分子活性与RNA结合，细胞组分(Cellular Component, CC)中线粒体核糖体/胞质核糖体占达66%；96 hpi转向分解代谢与氨基酸代谢，水解酶活性升高，PcRKNT8尤为明显，线粒体相关CC升至19%。RF按GO类别分类显示24 h MF准确率100%，脂质代谢过程与跨膜转运为最具判别力生物过程(Biological Process, BP)；96 h胞外区域与线粒体为CC最重要判别特征，反映向高能耗寄生外泌状态的转换。

3.4. Structural and functional reconfiguration of the Pc123 secretory proteome（Pc123分泌组结构与功能重配置）

含信号肽(Signal Peptide, SP)或跨膜区(Transmembrane, TM)的非胞质(Non-Cyt)蛋白在PcRKN与PcT8早现增多，PcRKNT8于96 h累积最强。序列特征空间PCoA显示24 h Pc与PcT8富集SP富含蛋白（细胞壁重塑/糖类代谢），PcRKN富集TM关联蛋白。Mantel检验表明在大多数场景下非胞质域定位距离与氨基酸组成特征空间相关性高于全长进化距离，即Pc分泌组由功能定位而非系统发生主导；仅PcRKNT8_96全长进化距离显著相关，暗示联合胁迫招募保守进化子集。RF全局模型识别top20非胞质判别蛋白含胞外水解酶(GH18几丁质酶A0A4Q7JT25、1,3-β-葡聚糖酶A0A4Q7JLT6)、细胞壁蛋白(CWP1–3)、过氧化氢酶等；分区RF显示胞外组top10与全局top20交叠6个（含最高权重Glucan 1,3-β-glucosidase A0A4Q7JLT6及CFEM域蛋白），膜组无重叠（主导为Cupin-type蛋白A0A4Q7K1K2、假定黏附素A0A4Q7K4R2），胞内/胞质组较稳定。

3.5. Chitosan-induced restructuring of the Meloidogyne javanica proteome（壳聚糖诱导M. javanica蛋白质组重构）

MELJA_Pc与MELJA_Pc_T8共有蛋白24 h相似度56.6%，96 h降至19.3%，PERMANOVA证实时间×处理显著。24 h仅卵黄蛋白原(Vitellogenin) A0A915MDA7显著上调；96 h四个EF-1α(A0A915LDY9等)＞6倍上调，两个ZP域蛋白(A0A915LHV2/MC03)下调。RF全局模型将卵黄蛋白原A0A915MDA7/A0A915LIP9列为最重要区分特征。RF、SVM与PLS-DA交集得13个严格共识生物标志物，涵盖ATP合酶β、EF-1α、Piwi域蛋白、ZP域蛋白、卵黄蛋白原。GLMM验证其时间动力学：壳聚糖下Vitellogenin早升晚降，ATP合酶β与EF-1α晚阶段特异升高，ZP域蛋白与Piwi受抑。GO富集示24 h壳聚糖组细胞分裂相关条目略升；96 h对照组ncRNA介导基因沉默与核/线粒体定位蛋白增多，壳聚糖组细胞骨架(17%)与染色体相关(23%)条目升高，RF判别线粒体相关特征最重要。

讨论与结论翻译

讨论指出Pc具从腐生向寄生转换的表型可塑性，可整合表面疏水性与卵源化学线索启动"宿主感应"程序。Pc蛋白质组呈时间层级：24 h为翻译启动/ priming阶段（蛋白合成机器与早期毒力因子上调），96 h转为分解代谢/降解谱（HSP70、ATP合酶及分泌蛋白酶为特征）。壳聚糖兼为胁迫因子与激发子(Elicitor)，促使Pc膜与分泌过程重组，上调糖类活性酶(CAZymes)与胞外水解酶及CFEM域蛋白，与寄生准备态相关联。线虫方面壳聚糖与三类改变相伴随——生殖储备相关(Vitellogenin扰动)、卵壳细胞外基质结构(ZP域蛋白抑制)、代谢应激(ATP合酶β上调、氧化应激标志)，单纯真菌感染诱导标准应激谱，壳聚糖叠加则伴随线虫涉及生殖/结构/代谢通路的更广泛重配置，可能影响寄生发生的生理背景。综上该三方互作中壳聚糖"推"动真菌向分泌与高代谢寄生程序，"拉"动线虫结构生殖通路蛋白质组偏移，为壳聚糖协同生防菌的可持续植病管理提供分子框架。

结论：本研究确立壳聚糖是Pochonia chlamydosporia 123与Meloidogyne javanica三方互作的强力调控因子。壳聚糖逐步将真菌蛋白质组导向以线粒体激活与扩展分泌组为特征的寄生高能程序，同时与线虫涉及生殖与结构通路的蛋白质组改变相伴随，可能改变卵壳微环境以利于寄生。机器学习架构(Random Forest、SVM、PLS-DA)的整合对鉴定稳健生物标志物与功能特征不可或缺，这些计算手段提供了增强寄生响应的预测分子图谱，可为可持续农业中壳聚糖基类杀线虫剂(Bionematicide)的优化提供可扩展框架。