《Agronomy》：Integrated Sphingolipidomic and Transcriptomic Analysis Reveals Resistance Mechanisms to Bacterial Wilt in Peanut Varieties with Contrasting Tolerance
由细菌Ralstonia solanacearum引起的青枯病是一种毁灭性病害，限制了花生的生产。区分抗病与感病反应的分子机制仍在不断探索中；然而，关于脂质代谢与转录重编程之间相互作用的数据有限。本研究对具有不同抗性的花生品种——抗病品种中花6号（ZH6）和感病品种中花12号（ZH12）——在发病前、发病初期和发病晚期的根部进行了整合鞘脂组学和转录组学分析。抗病品种ZH6表现出特定葡萄糖神经酰胺（GluCer t18:1/h25:0和GluCer t18:0/h23:0）和肌醇磷酸神经酰胺（IPC t18:0/h24:0）的早期积累。代谢趋于稳定，在早期和晚期之间未观察到显著的脂质变化。这一反应得到了脂肪酸α-羟化酶（0538LJ）、中性/碱性神经酰胺酶（IX12GR、KZ47MP.1）和ABA信号传导组分（ABA受体、SRK2A-like）转录激活的支持。相比之下，感病品种ZH12表现出进行性鞘氨醇耗竭以及极长链脂肪酸神经酰胺d18:0/16:0和固醇的延迟积累。KCS4在ZH12中上调，为VLCFA-神经酰胺的积累提供了转录基础。激素信号传导的差异显而易见：ZH6表现出早期的ABA信号传导，随后进入转录静止状态，而ZH12在病理性神经酰胺积累后表现出延迟的ACS4/ACS11过度诱导。这伴随着生长素耗竭（GH3.9上调）、细胞分裂素磷酸中继传递失败（HPt6抑制）以及SA信号传导的解偶联（PR-1诱导而无TGA激活）。因此，抗性由早期协调的鞘脂重塑和ABA信号传导所定义，最终导致稳态稳定。相比之下，感病性表现为延迟的神经酰胺积累和多激素信号传导抑制。这一多组学框架提供了详细的脂质组学和转录组学特征，可用于鉴定候选基因和脂质生物标志物，以进行抗青枯病花生品种的标记辅助育种。
整合鞘脂组学与转录组学揭示花生青枯病抗性机制的研究解读
一、研究背景与科学问题
花生（*Arachis hypogaea* L.）是全球重要的油料和经济作物，预计2026/27年度全球产量将达到5151.9万吨，种植面积约3000万公顷。中国贡献了全球近三分之一的产量，年均产量超过1800万吨。然而，花生产业面临着病害抗性不足、机械化水平参差不齐以及传统产区种植面积萎缩等多重挑战。由土传细菌青枯菌（*Ralstonia solanacearum*）引起的青枯病是危害花生产量最严重的细菌性病害之一，该病原菌在十大植物病原细菌中科学和经济重要性排名第二，可侵染超过450种植物，包括花生、烟草、香蕉和番茄等。在中国，超过10%的花生种植面积受青枯病影响，主要集中在中部和南部地区，造成10%–30%的产量损失，严重时可达50%–100%，全球每年因此造成的损失超过10亿美元。
青枯菌通过根部伤口、根尖和侧根萌发点侵入植物，定殖于根皮层和维管系统，特别是木质部导管，在其中繁殖并产生胞外多糖阻碍水分运输，最终导致植株萎蔫死亡。该病原菌寄主范围广、可在土壤中长期存活，使得防治极为困难。传统管理策略如轮作、土壤改良和化学防治效果有限，而生物防治仅取得部分成功，因此选育和利用抗病品种被认为是最经济有效的策略。中国已培育出多个高抗青枯病花生品种，包括中花6号、鄂花6号和豫花14号等。遗传研究已在多条染色体上定位到青枯病抗性相关数量性状位点（QTLs），并鉴定出多个候选抗性基因。然而，抗性和感病反应的分子机制仍未被完全阐明，特定代谢途径（尤其是脂质代谢途径）在病害进展中的作用尚未被系统表征。
脂质在植物防御中发挥核心作用：脂氧合酶途径产生包括茉莉酸（JA）在内的防御相关oxylipins；膜脂组成在感染过程中被重塑；极长链脂肪酸参与防御信号传导。鞘脂类作为脂质中的重要类别，既是结构膜组分，又是生物活性信号分子。植物鞘脂分为四类：游离长链碱基（LCBs）、神经酰胺、糖基神经酰胺和糖基肌醇磷酸神经酰胺。这些鞘脂调节膜流动性和通透性以及膜结合蛋白活性，还参与程序性细胞死亡、病原体诱导的超敏反应、脱落酸（ABA）依赖性气孔关闭、寄主-病原互作和低温信号转导等信号转导途径。
二、研究目的与方案
鉴于鞘脂在植物免疫中的重要作用以及花生青枯病研究中脂质代谢与转录调控相互作用的资料匮乏，本研究旨在通过整合鞘脂组学和转录组学分析，系统比较抗病品种中花6号（ZH6）和感病品种中花12号（ZH12）在青枯病不同进展阶段的分子动态变化，探究花生青枯病抗性是否与特定鞘脂种类的早期持续积累以及鞘脂代谢基因转录的协调激活相关，鉴定与抗性和感病性相关的关键鞘脂种类，并描绘鞘脂和脂肪酸代谢及植物-病原互作途径相关基因的转录组动态。
研究人员在四川省南充市农业科学院青枯病病圃（具有明确青枯菌 infestation 历史）种植两个花生品种，于2024年5月17日播种，每品种30粒，株距20厘米。病圃为天然发病，未进行人工接种。分别在发病前（顶部叶片尚未萎蔫，ZH6样本记为HR1，ZH12样本记为HS1）、发病初期（顶部叶片开始萎蔫，记为HR2/HS2）和发病晚期（全株萎蔫，记为HR3/HS3）三个时间点采集根部组织。每个品种和时间点组合设三个生物学重复，共18个样本（2品种×3阶段×3重复）。所有根部组织在指定时间点立即采集，液氮冷冻后保存于?80°C。
三、主要关键技术方法
本研究采用了以下主要技术方法：一、基于液相色谱-串联质谱（LC-MS/MS）的靶向鞘脂组学分析：使用岛津Nexera 20AD-HPLC耦合SCIEX QTRAP 6500 PLUS质谱仪，在正离子模式下对鞘脂进行定量分析，游离固醇和固醇酯则采用大气压化学电离（APCI）模式分析。内标包括C17-神经酰胺（d18:1/17:0）、C12-鞘氨醇（d12:0）和C17-葡萄糖神经酰胺（d18:1/17:0）。二、基于Illumina平台的转录组测序（RNA-seq）：提取总RNA后构建文库，在Illumina平台测序，平均每样本产生约1.1亿条原始读段，过滤后获得约1.08亿条有效读段，Q30值均超过96.9%。有效读段比对到花生参考基因组Tifrunner.gnm2.ann1.4K0L，比对率62.2%–96.5%。三、实时荧光定量PCR（qRT-PCR）验证：选取多个候选基因设计特异性引物，以AhActin7为内参基因，在AI7500仪器上使用Vazyme Biotech qPCR SYBR Green Master Mix进行验证，每个样本设三个生物学重复，实验重复三次。四、生物信息学分析：使用fastp进行质控、HISAT2比对、StringTie组装、RSEM定量，DESeq2进行差异表达分析，NR、GO、KEGG、eggnog、SWISS-Prot和Pfam数据库进行功能注释。脂质组学数据使用R软件（v 4.3.3）进行主成分分析（PCA）、Welch's ANOVA和Games-Howell事后检验。五、样本来源为四川省南充市高坪区花园村南充市农业科学院青枯病病圃。
四、主要研究结果
4.1 抗感品种鞘脂组学比较分析揭示不同的脂质重塑动态
比较脂质组学分析在ZH6和ZH12根部三个采样阶段共鉴定出13种脂质类型（103种化合物），包括PhytoCer、PhytoCer-OHFA和Phyto-GluCer等。主成分分析显示样本主要按处理和样本类型分组，PC1和PC2分别解释31.3%和15.7%的变异。大多数鞘脂类别（t-S1P、S1P、phytoSph、Sph、Cer和PhytoCer）在发病前和发病初期HR中的含量高于HS，但到发病晚期，部分类别（t-S1P、PhytoCer和PhytoCer-OHFA）在HS中的含量反而高于HR。甘油酯中，GluCer和Phyto-GluCer在HR1中低于HS1，但在HR2和HR3中分别高于HS2和HS3。IPC和GIPC在所有阶段HR中均持续较高。固醇酯在HS中所有阶段均高于HR，而游离固醇则呈现相反模式。
在品种内部比较中，感病品种ZH12在发病初期（HS2 vs. HS1）有3种脂质显著变化：鞘氨醇d18:1下降44.6%，葡萄糖神经酰胺d18:1/h24:0增加46.8%，总植物葡萄糖神经酰胺增加25.0%。到发病晚期（HS3 vs. HS1），有10种脂质显著变化：鞘氨醇d18:1进一步降至发病前的49.1%；神经酰胺d18:0/16:0惊人地增加11.98倍；总神经酰胺增加2.25倍；多种葡萄糖神经酰胺增加（d18:0/h16:0增加2.73倍，t18:1/h25:0增加1.68倍，t18:1/h26:0增加1.38倍）；谷甾醇增加1.54倍，豆固醇增加1.46倍，总固醇增加1.34倍。发病初期至晚期（HS3 vs. HS2）仅植物神经酰胺t18:1/18:0显著下降28.1%。这些变化表明感病品种ZH12在晚期积累神经酰胺（特别是增加了12倍的d18:0/16:0）和固醇，而鞘氨醇d18:1水平随时间进行性下降。
抗病品种ZH6在发病初期（HR2 vs. HR1）有4种脂质显著增加：葡萄糖神经酰胺t18:1/h25:1（2.73倍）、葡萄糖神经酰胺t18:1/h25:0（2.53倍）、葡萄糖神经酰胺t18:0/h23:0（2.46倍）和肌醇磷酸神经酰胺t18:0/h24:0（2.32倍）。而发病初期至晚期（HR3 vs. HR2）未检测到显著脂质变化（所有p > 0.05）。从发病前到晚期（HR3 vs. HR1），两种脂质保持显著升高：葡萄糖神经酰胺t18:1/h25:0增加2.26倍，葡萄糖神经酰胺t18:0/h28:0增加2.39倍。这表明ZH6表现出特定葡萄糖神经酰胺的早期持续增加，发病初期后无显著变化，反映了代谢稳定化。
品种间比较显示，发病前（HS1 vs. HR1）仅一种脂质有显著差异：感病品种中葡萄糖神经酰胺t18:1/h24:0浓度是抗病品种的1.69倍。发病初期（HS2 vs. HR2），感病品种有4种植物神经酰胺显著低于抗病品种：t18:1/24:0（0.54倍）、t18:0/26:0（0.54倍）、t18:0/h22:1（0.59倍）和t18:0/h24:1（0.49倍）。到发病晚期（HS3 vs. HR3），两品种间无显著脂质差异。品种间最大差异出现在发病初期，此时感病品种的4种植物神经酰胺显著偏低。
4.2 鞘脂代谢相关基因的差异表达分析
在感病品种ZH12中，PT6J7Z基因在HS2 vs. HS1中显著上调（log2FC 9.10），在HS3 vs. HS2中下调。两个可能参与复杂糖基肌醇磷酸神经酰胺（GIPCs）极性头基修饰的α-半乳糖苷酶基因（52QCLV.2和CM8MKE）在HS2 vs. HS1中下调，并在HS3 vs. HS1中进一步受抑制。感病品种中其他鞘脂相关基因没有显著变化，表明在青枯病进展期间该途径的转录激活有限。
在抗病品种ZH6中，鞘脂代谢和防御相关修饰基因广泛转录激活。脂肪酸α-羟化酶（0538LJ）——负责在极长链脂肪酸（VLCFAs）C-2位引入羟基以合成稳定的植物神经酰胺和GIPCs——在HR2 vs. HR1中上调（24.44 log2FC）和HR3 vs. HR1中上调（24.10 log2FC），这与脂质组学中HR植物神经酰胺丰度较高一致。中性/碱性神经酰胺酶IX12GR和KZ47MP.1（水解神经酰胺产生信号活性LCBs和游离脂肪酸）在HR2 vs. HR1（分别为2.49和2.53 log2FC）和HR3 vs. HR1（分别为1.86和1.61 log2FC）中上调。相比之下，另一种中性/碱性神经酰胺酶MG688V在所有阶段感病品种中始终低于抗病品种（HS3 vs. HR3为2.49 log2FC，HS2 vs. HR2为2.16 log2FC，HS1 vs. HR1为2.99 log2FC），表明ZH12中神经酰胺分解代谢的基础水平受到抑制。
值得注意的是，含saposin B结构域的蛋白18JCSN和RCXU7Z（通常作为膜破坏性抗微生物防御肽或靶向植物液泡的结构组分）在HR2 vs. HR1（分别为1.32和1.49 log2FC）和HR3 vs. HR1（分别为1.34和1.29 log2FC）中上调。非溶酶体葡萄糖神经酰胺酶（Y7DSHF）——负责将葡萄糖神经酰胺转化回神经酰胺——在HR2 vs. HR1（1.48 log2FC）和HR3 vs. HR1（2.42 log2FC）中上调。β-己糖胺酶2（Y5DXIN.1）也在HR2 vs. HR1和HR3 vs. HR1中上调。
品种间比较显示，KZ47MP.1在HR1中较HS1显著上调（2.27 log2FC），表明抗病品种在强烈感染前具有较高的基础神经酰胺酶表达。发病初期（HS2 vs. HR2），大多数鞘脂相关DEG在感病品种中较抗病品种显著下调，包括saposin B蛋白（18JCSN、RCXU7Z）、神经酰胺酶（IX12GR、MG688V）、β-己糖胺酶2（Y5DXIN.1）和非溶酶体葡萄糖神经酰胺酶（Y7DSHF）。到发病晚期（HS3 vs. HR3），仅IX12GR神经酰胺酶在感病品种中保持显著下调。
4.3 脂肪酸生物合成、延伸和降解相关基因的差异表达分析
在感病品种ZH12中，多个长链酰基-CoA合成酶（0DV40D、D2PEMM、H4YX61和R34KW1）及一个根质体被膜蛋白（7FN9BJ）在感染后早、晚期较发病前显著上调。7FN9BJ表现出持续上升的诱导模式（HS2 vs. HS1为9.18 log2FC，HS3 vs. HS1为22.87 log2FC）。在脂肪酸延伸方面，3-酮脂酰-CoA合酶4（KCS4，70Q13X）——脂肪酸延伸超过16或18个碳的限速酶——在HS2 vs. HS1（8.89 log2FC）和HS3 vs. HS1（9.43 log2FC）中显著持续上调，并在晚期感染中保持升高（HS3 vs. HS2为0.58 log2FC）。这一显著诱导为感病品种晚期脂质组学中VLCFA神经酰胺（C22:0、C24:0、C24:1和C26:0）的过度积累提供了直接的转录依据。相反，KCS12（Q5DVPF）在HS3 vs. HS1中下调（-1.64 log2FC），表明发病期间存在高度异构体特异性的调控。分解代谢方面，烯酰-CoA水合酶（CY1TAH）和醛脱氢酶（CR7105）在晚期受抑制，而过氧化物酶体酰基-CoA氧化酶（ACX，79B2B8）在早期（HS2 vs. HS1）短暂上调但未能持续表达。
相比之下，ZH6在脂肪酸代谢网络中表现出更广泛更持续的转录反应。多种关键酶在早期和晚期相对于发病前持续上调，包括多个长链酰基-CoA合成酶（D2PEMM、H4YX61、QUY0YV）、酰基-ACP硫酯酶（TC8Y09、V6JT5K.1）、酰基-ACP脱饱和酶（YQG5EZ、ZL0G7D）和3-羟基酰基-ACP脱水酶（ZLL7AW）。酰基-ACP脱饱和酶的强烈上调支持向不饱和脂肪酸生产的代谢通量增加。重要的是，KCS4和KCS12在抗病品种中未因病原菌侵染而上调，KCS4反而在HR2 vs. HR1（-1.39 log2FC）和HR3 vs. HR1（-1.20 log2FC）中略微下调，这与HR脂质谱中不含VLCFA神经酰胺积累的现象相符。脂肪酸降解方面，抗病品种同时上调过氧化物酶体β-氧化关键组分，包括长链酰基-CoA合成酶、ACXs（79B2B8和ZVY5TW）和烯酰-CoA水合酶（CY1TAH），表明脂肪酸分解代谢高度活跃。
品种间比较显示，抗病品种在感染前就具有多种脂肪酸生物合成基因（特别是酰基-ACP硫酯酶和脱水酶）的较高基础表达。值得注意的是，KCS4转录本丰度在HR1中显著高于HS1（HR1 vs. HS1为-9.21 log2FC），表明抗病品种在遇到病原体前具有极高的KCS4组成型基础水平。感染晚期，抗病品种较感病品种表现出脂肪酸激活基因（01GDC8.2、B459M5和G1SIKM）的更高表达。
4.4 α-亚麻酸和亚油酸代谢相关基因的差异表达分析
在感病品种ZH12中，发病早期（HS2 vs. HS1）诱导了两个α/β-水解酶超家族蛋白（41Y8SN/AW6SSY和Y6BJG8/Y8803N），而磷脂酶A2家族蛋白（96ZPFV/N1SVVP）、12-氧植物二烯酸还原酶2（OPR2，L5Q3AX/QLS1U0）和第二个磷脂酶A2（MS236J/XYC619）表达降低，表明茉莉酸前体释放受到短暂抑制。发病晚期与发病前比较（HS3 vs. HS1）鉴定出23个差异表达基因，下调基因包括亚油酸13S-脂氧合酶（13-LOX，3BDP96/5I9T0U/Y6Q89P）和α-双加氧酶（H8DQQQ/0S2PTK）；上调基因包括OPR2、4-香豆酸：CoA连接酶2、ACX1、酰基-CoA合成酶5、乙醇脱氢酶1、丙二烯氧化物环化酶3、α-双加氧酶、AMP依赖性合成酶/连接酶、异戊酰-CoA脱氢酶、13-LOX、patatin样磷脂酶、过氧化物酶体3-酮脂酰-CoA硫解酶3和锌结合乙醇脱氢酶。
在抗病品种ZH6中，早期（HR2 vs. HR1）出现13-LOX、α/β-水解酶、ACX2、AMP依赖性合成酶/连接酶、4-香豆酸：CoA连接酶2、过氧化物酶体硫解酶和种子亚油酸9S-脂氧合酶（9-LOX）的上调；同一下调基因包括13-LOX、磷脂酶A2、丙二烯氧化物合酶和其他9-LOX旁系同源物。上调和下调旁系同源物的共存表明脂氧合酶通量存在异构体特异性分配。
品种间比较显示，发病前（HS1 vs. HR1），ZH12中α-双加氧酶、OPR2和3-羟基酰基-CoA脱氢酶的表达低于ZH6。发病初期（HS2 vs. HR2），ZH12中三个α-双加氧酶的表达高于ZH6。发病晚期（HS3 vs. HR3），ZH12中多个OPR2旁系同源物、锌结合乙醇脱氢酶和S-腺苷-L-甲硫氨酸甲基转移酶的表达较高，而ZH6中一组不同的基因（包括上述甲基转移酶和一个未知蛋白）表达较高，表明晚期oxylipin加工途径存在差异。
4.5 植物-病原互作通路相关基因的差异表达分析
在感病品种ZH12中，感染早期（HS2 vs. HS1）下调环核苷酸门控通道CNGC15（4LCW9K.2、51I078、4X6KYU.2），表明早期钙内流受抑制。同时上调PR-1家族成员、HSP90.1、免疫受体（CC-NBS-LRR、TIR-NBS-LRR和LRR-NB-ARC）和RBOHB。发病晚期（HS3 vs. HS1），上调转录本显著富集于信号和执行网络，包括CDPK1/20、WRKY转录因子、RBOHB、木瓜蛋白酶样半胱氨酸蛋白酶和衰老相关羧酸酯酶101样。时序追踪（HS3 vs. HS2）显示晚期防御反应显著，但这一转录重编程发生太晚，无法阻止以Cer d18:0/16:0增加11.98倍为特征的广泛病理性脂质崩溃。
在抗病品种ZH6中，早期（HR2 vs. HR1）表现出加速的高度协调的防御程序。强烈上调的转录本包括PR-1变体、木瓜蛋白酶、HSP89.1、NLR受体、EDS蛋白和RBOHA/B。这一主动反应恰好与结构葡萄糖神经酰胺（GluCer t18:1/h25:0增加2.53倍，GluCer t18:0/h23:0增加2.46倍）和肌醇磷酸神经酰胺核心（IPC t18:0/h24:0增加2.32倍）的丰度增加同步，表明膜脂微区获得了免疫支持性强化。晚期感染（HR3 vs. HR2），ZH6选择性下调KCS12且缺乏KCS4诱导，从而避免细胞毒性神经酰胺积累，在早期和晚期感染之间实现绝对的脂质组学静止。
品种间比较证实，抗性可能受预激和代谢控制支配。HS1 vs. HR1中，感病品种表现核心免疫调节因子、钙通道和脂质组分（如WRKY、CC-NBS-LRR、KCS4、CML41-like和CNGCs）的抑制，表明ZH6具有更高的基础防御准备水平。发病初期（HS2 vs. HR2），感病品种钙信号组分（WRKY、CML41-like和CNGCs）普遍缺乏，这直接对应其结构植物神经酰胺库的严重耗竭（PhytoCer t18:1/24:0 = 0.54，PhytoCer t18:0/26:0 = 0.54，PhytoCer t18:0/h22:1 = 0.59和PhytoCer t18:0/h24:1 = 0.49），将鞘脂缺乏与早期信号传导受损联系起来。发病晚期（HS3 vs. HR3），虽然脂质丰度总体趋同，但内部遗传程序仍高度分化。
4.6 MAPK和植物激素信号通路相关基因的差异表达分析
在感病品种ZH12中，发病早期（HS2 vs. HS1）上调转录本包括蛋白激酶超家族成员、PR-1蛋白、几丁质酶家族成员、MAPK底物MKS1-like、RBOHB、乙烯传感器和1-氨基环丙烷-1-羧酸合酶11（ACS11）。到发病晚期（HS3 vs. HS1；HS3 vs. HS2），上调转录本包括WRKY和bZIP因子、V类几丁质酶、PR-4家族成员、PP2C调节因子、MKS1、RBOHB、ACS4、ACS11和bHLH因子。ACS4和ACS11这种不断增强的晚期诱导强调了强大但严重延迟的乙烯驱动反应——仅在不可逆的病理性脂质崩溃完全建立后才动员。
在抗病品种ZH6中，感染早期（HR2 vs. HR1）部署了主动的ABA介导防御程序。上调转录本显著包括ABA受体、SnRK2家族激酶SRK2A-like、RBOHA/B/F、bZIP因子、ACS11和PR-1蛋白，而ABA信号的关键负调节因子（PP2Cs）和乙烯周转因子（EIN3结合F-box蛋白1）显著受抑制。早期ABA信号激活在早期和晚期之间（HR3 vs. HR2）转变为稳定转录状态，仅一个PP2C变异体被差异调控，与ZH6中观察到的脂质组学完全静止平行。
品种间比较显示不同的基础和感染后激素变化在病害反应中起作用。发病前（HS1 vs. HR1），WRKY/bZIP转录因子和丝氨酸/苏氨酸蛋白激酶1在抗病品种中表达较高；感病品种则激活CML20、P型ATP酶铜/重金属外排转运蛋白和组织酸激酶1。发病初期（HS2 vs. HR2），转录组差异较窄，与HS2中结构植物神经酰胺屏障的严重品种间差异形成对比。发病晚期（HS3 vs. HR3），感病品种通过ACS4、ACS11和PR-4几丁质酶诱导表现乙烯信号，而抗病品种通过维持WRKY转录因子、受体样蛋白激酶2和PP2C负调节因子的相对较高转录水平来下调不必要的应激级联反应。
qRT-PCR分析显示所选基因的表达趋势与转录组测序数据集一致，证实RNA测序数据高度可靠且具有生物学有效性。
五、讨论和结论总结
在讨论部分，研究人员指出，抗病和感病品种在鞘脂重塑动力学上表现出显著差异。抗病品种ZH6对青枯菌的抗性可能受高度协调的根部鞘脂组的时空丰度调控——感染触发葡萄糖神经酰胺和IPC的立即积累并维持其含量，表明抗病基因型可能进行根部膜的结构性强化。脂肪酸α-羟化酶（0538LJ）的上调使α-羟基化复合鞘脂的合成成为可能，从而稳定脂双层抵抗细菌降解，直接支持了观察到的植物神经酰胺和IPC增加。中性/碱性神经酰胺酶IX12GR和KZ47MP.1的同时上调表明具有信号能力的LCBs的活跃产生。而感病品种ZH12表现为不同的鞘脂谱：晚期信号活性长链碱基鞘氨醇（d18:1）降至发病前水平的49.1%，同时伴随神经酰胺（d18:0/16:0）及游离谷甾醇和豆固醇的积累。神经酰胺酶（MG688V）在感病背景下所有阶段持续受抑制，提示ZH12中神经酰胺分解代谢受损。品种间比较显示，发病初期感病品种四种关键植物神经酰胺（t18:1/24:0，t18:0/26:0，t18:0/h22:1，t18:0/h24:1）丰度低，而极长链植物神经酰胺和复合鞘脂可与固醇互作组织紧密堆积的质膜微结构域，这些纳米结构域物理聚集免疫传递所需的模式识别受体和CNGC通道。saposin B防御蛋白（18JCSN、RCXU7Z）在ZH6中的选择性高表达表明其持续准备进行鞘脂介导的根部防御。
脂肪酸代谢的转录调控可能决定极长链神经酰胺谱和病害结果。感病品种ZH12中KCS4在发病初期和晚期上调，可能驱动含VLCFA神经酰胺（C22:0、C24:0、C24:1、C26:0）的丰度。而ZH6发病前组织中KCS4表达高于ZH12，表明其具有一致的基线脂质组学特征，可满足正常生理VLCFA库需求而无需感染触发的基因表达峰值，从而成功防止毒性简单VLCFA-神经酰胺的积累。抗病品种还执行了更广泛更持续的辅助脂质代谢网络激活，如酰基-ACP脱饱和酶（YQG5EZ、ZL0G7D）上调支持不饱和脂肪酸通量增加，可能调节膜流动性并优化整合免疫受体蛋白的操作构象。ZH6中过氧化物酶体β-氧化级联的诱导（ACXs和烯酰-CoA水合酶在早期和晚期均上调）可能驱动代谢能量产生和活性oxylipin防御信号传导。
激素信号在抗性和感病品种间存在分歧。ZH6中感染触发ABA受体表达变化、PP2C途径负调节因子抑制和SRK2A-like激酶诱导，与ABA作为青枯病抗性正调节因子的新兴作用一致。ZH6在初始ABA激活后维持的转录表达与该品种脂质组学稳定一致——一旦ABA介导的防御机制被激活，系统无需持续的能量密集型免疫级联激活即可稳定。而感病品种ZH12表现出严重延迟和脱偶联的激素轨迹——在细胞毒性神经酰胺积累已造成膜损伤后，乙烯生物合成关键基因ACS4和ACS11晚期上调，表明乙烯在ZH12中可能不再作为防御触发因子，而是作为适应不良的应激反应，加速组织衰老、细胞自溶和维管塌陷。感病性还因多激素信号的更广泛转录变化而加剧：生长素途径中GH3.9生长素酰胺合成酶上调暗示活性生长素信号耗竭；HPt6抑制和APRR2样阻遏物积累表明细胞分裂素信号所需的多步双组分磷酸中继系统崩溃；BRI1激酶在所有发病阶段上调表明活性油菜素甾醇信号传导受抑制；TGA家族成员在晚期ZH12中受抑制表明SA依赖性基因表达向保护性和结构性防御转化失败。这些详细激素信号表达特征表明，花生根部的感病性不仅是防御激活失败，更是激素信号网络的主动适应不良性重编程。
研究结论明确指出，抗病品种ZH6的抗性特征为特定葡萄糖神经酰胺和肌醇磷酸神经酰胺的早期积累、脂肪酸去饱和与过氧化物酶体β-氧化的协调激活以及ABA信号传导的启动，这些过程在晚期感染阶段趋于稳定。相比之下，感病品种ZH12的感病性涉及进行性鞘氨醇耗竭、极长链神经酰胺和固醇的延迟过度积累以及激素信号传导的时间脱节——ABA激活缺失，延迟的乙烯反应无法阻止代谢崩溃。感病背景还表现出生长素耗竭、细胞分裂素磷酸中继系统崩溃、油菜素甾醇抑制和SA信号解偶联的转录证据，表明感病性是激素网络主动重编程的结果而非被动防御失败。青枯菌感染花生的结果取决于鞘脂代谢、脂肪酸延伸和激素信号通路的时间、协调和整合。晚期脂质谱的趋同尽管转录程序存在差异，支持将感病性解释为一种替代性的适应不良代谢轨迹，而非简单的抗性机制缺陷。未来工作应聚焦于候选基因的功能验证，如KCS4、α-羟化酶（0538LJ）和神经酰胺酶IX12GR与MG688V，通过过表达或CRISPR介导敲除在花生或模式系统中实现；需要直接定量感染不同阶段的oxylipins和植物激素以补充转录推断；使用化学抑制剂（如FB1）检测神经酰胺积累是否驱动或伴随细胞死亡将阐明因果关系；诊断性脂质和候选基因生物标志物向育种流程的转化需在多样化花生种质和田间条件下验证。
