整合转录组学与生理学解析甜瓜抗细菌性果斑病早期防御的关键基因与调控网络

《Horticulturae》：Integrated Transcriptomic and Physiological Analyses Reveal Key Genes and Regulatory Network for Early-Stage Defense Against Bacterial Fruit Blotch in Melon Dan Zhou, Jianli Shang, Lanying Lu, Jianfei Song, Nannan Li, Shuangwu Ma and Na Li

本研究通过整合转录组、生理与激素谱学分析，系统揭示了抗/感甜瓜种质在细菌性果斑病（BFB）早期（接种后12-24小时）的防御机制差异。研究表明，抗性植株能快速启动以水杨酸（SA）信号为主导、核心免疫通路（如植物-病原互作、MAPK通路）激活、苯丙烷代谢上调的协同防御网络，而感病植株则反应延迟、失调，表现为SA被抑制、茉莉酸（JA）激增及多通路广泛但非协同的激活，最终导致生理紊乱。通过比较分析，鉴定了9个早期响应及21个维持抗性的关键基因，为解析甜瓜抗细菌病害的分子与生理机制提供了整合理论框架。

　　

1. 引言

作为葫芦科重要成员，甜瓜是全球重要的经济作物。细菌性果斑病（BFB）由西瓜嗜酸菌（Acidovorax citrulli, Aac）引起，是危害甜瓜产业最具破坏性的细菌性病害。尽管已有抗BFB甜瓜品种的筛选鉴定研究，但目前尚无商业化抗病品种。因此，挖掘抗性种质并深入研究其抗性机制，对减轻BFB危害至关重要。

2. 材料与方法

本研究以高抗甜瓜种质ZT145和高感种质ZT146为材料。在幼苗一叶一心期，用浓度为1×10⁶CFU·mL^-1的Aac菌悬液喷雾接种叶片，对照组喷施无菌水。分别在接种后0、12、24小时（hpi）采集叶片样本，用于转录组测序、防御相关酶活性（PAL、SOD、CAT、PPO）测定以及植物激素（SA、JA、ABA、IAA等）含量分析。RNA测序采用Illumina平台PE150策略，使用HISAT2将Clean reads比对到甜瓜参考基因组（Cucumis melo DHL92 v4.0），利用DESeq2和edgeR进行差异表达基因（DEGs）筛选（|log₂FC| ≥ 1, FDR < 0.05），并进行GO和KEGG富集分析。通过qRT-PCR对9个候选基因进行表达验证。

3. 结果

3.1. 病害表型与生理生化指标

接种后，感病种质ZT146在12 hpi即出现水渍状病斑，24 hpi病情迅速扩展，72 hpi植株接近死亡。而抗病种质ZT145直到24 hpi才出现少量小病斑。

激素分析显示，抗病植株的整体SA水平显著高于感病植株。在12-24 hpi，抗病植株的SA被病原菌特异性诱导（接种组 > 对照组），而感病植株的SA生物合成似乎被抑制（接种组 < 对照组）。JA水平在感病植株中持续上升，24 hpi接种组激增；而抗病植株的JA含量在12-24 hpi呈下降趋势。IAA含量变化在两种基因型中呈相反调控方向。ABA响应模式复杂。

抗病植株的CAT、SOD和PAL活性均显著高于感病植株，而PPO活性较低。

3.2. 转录组测序与差异表达基因分析

对24个样本进行测序，数据质量良好，样本重复性好。与0 h相比，抗病种质ZT145在12 hpi和24 hpi分别鉴定到4647和3993个DEGs，感病种质ZT146则分别有5459和7699个DEGs。比较接种与未接种样本，抗病植株在12 hpi和24 hpi分别有109和969个DEGs，而感病植株分别有458和6184个DEGs。

3.3. KEGG与GO富集分析

KEGG分析表明，在12 hpi，抗病植株显著上调了苯丙烷生物合成、植物-病原互作、MAPK信号通路-植物等关键通路。到24 hpi，富集通路还包括类黄酮生物合成、芪类化合物生物合成等。感病植株在12 hpi富集了苯丙烷生物合成、苯丙氨酸代谢、亚油酸代谢等通路，到24 hpi，植物-病原互作等病害相关通路才被富集。

GO分析显示，抗病植株在12和24 hpi均持续上调了“响应氧化应激”、“肉桂酸生物合成过程”、“先天免疫应答”等防御相关过程，同时光合作用相关术语被下调。感病植株早期（12 hpi）有苯丙烷代谢术语上调，但到24 hpi减弱，光合作用相关过程出现失调。

3.4. 抗病相关关键基因挖掘

通过比较分析，鉴定了44个在所有四个接种处理组中共同的DEGs（核心响应基因）。进一步聚焦抗病植株在12 hpi的特异性响应，通过严格筛选得到了9个关键候选基因，它们富集于植物-病原互作、苯丙烷生物合成、MAPK信号、植物激素信号转导和磷脂酰肌醇信号系统五个关键通路。这9个基因是：MELO3C001988、MELO3C003421、MELO3C003632、MELO3C003889、MELO3C005427、MELO3C014624、MELO3C015216、MELO3C019528、MELO3C029753。

到24 hpi，在抗病植株中鉴定出152个特异性DEGs，其中重点关注了与植物-病原互作、MAPK信号通路、植物激素信号转导和苯丙烷生物合成这四个抗病相关通路相关的21个基因，它们可能在感染后期维持抗性中发挥作用。

3.5. R基因、MAPK和MLO基因分析

在12 hpi，无论抗感基因型，所分析的104个基因（包括70个R基因、20个MAPK基因和14个MLO基因）均未检测到显著的诱导差异表达。到24 hpi，仅在抗病植株中检测到6个R基因和1个MAPK基因被诱导上调。而在感病植株中，则出现了19个R基因、10个MAPK基因和5个MLO基因的差异表达，涵盖上调和下调模式。

3.6. qRT-PCR验证

随机选择的9个基因的qRT-PCR验证结果显示，其表达变化趋势与RNA-Seq结果一致，证实了转录组数据的可靠性。

4. 讨论

比较转录组分析揭示了抗病（ZT145）和感病（ZT146）甜瓜植株应对BFB的防御机制存在根本差异。抗病植株ZT145早在12 hpi就启动了高效、特异的防御程序。其SA水平被病原菌特异性诱导，同时核心通路如“植物-病原互作”、“MAPK信号”和“苯丙烷生物合成”早期即被富集。值得注意的是，在12 hpi，尽管KEGG/GO分析显示“植物-病原互作”和“MAPK信号通路”已被激活，但所分析的104个R、MAPK、MLO基因均无显著差异表达。这一看似矛盾的现象可能源于早期植物免疫信号传导（PTI/ETI）主要依赖于翻译后修饰，而观察到的转录组变化可能是这些早期信号事件导致的转录后果。到24 hpi，仅少数R和MAPK基因在抗病植株中被特异性诱导，反映了一种高度精确有序的防御反应。

抗病植株持续较高的PAL活性为转录水平上“苯丙烷生物合成”等通路的持续富集提供了生化基础。同时，较高的CAT和SOD活性表明抗病植株拥有更强大的抗氧化缓冲系统。激素数据显示抗病植株内源激素网络处于积极的再平衡状态，这与其资源重分配的生理基础相一致，表明植株主动实施了“生长-防御权衡”策略。

相比之下，感病植株ZT146的反应表现出延迟且协调性差的特点。其SA合成在接种后被抑制，可能导致早期防御信号转导失败。JA水平大幅且持续上升，GO分析也显示其早期反应以“茉莉酸介导的损伤响应”及相关负调控通路为主。在24 hpi，感病植株中“先天免疫应答的负调控”和“蛋白酶体”通路的富集，以及众多通路广泛的上调和下调，表明植株产生了一种广泛但协调性差的响应。尽管在24 hpi检测到“植物-病原互作”通路的上调，但晚于抗病植株，且防御响应分散在多个通路，可能是一种不足以有效遏制病原菌的广义胁迫反应。

5. 结论

本研究阐明，甜瓜对BFB的抗性不依赖于防御相关基因的绝对数量，而取决于植株在感染早期建立和维持协同防御网络的能力。该网络以SA主导的信号、强大的代谢和抗氧化支持以及整合的信号转导为特征。抗病植株利用此网络实现了防御资源的战略性调控，而感病植株由于该网络关键节点的失效，表现出延迟且协调性差的响应，最终不足以遏制病原菌增殖。研究鉴定的9个早期响应基因及21个维持抗性相关基因为后续功能验证和抗病育种提供了重要靶点。