木薯中介导对Xpm抗性的2-氧戊二酸依赖性双加氧酶的整合基因组与转录组分析

《Scientia Horticulturae》：Integrated genomic and transcriptomic analysis of 2-oxoglutarate-dependent dioxygenases mediating resistance to Xpm in cassava

　　

木薯细菌性枯萎病（cassava bacterial blight，CBB）由木薯黄单胞菌致病变种（Xanthomonas phaseoli pv. manihotis，Xpm）引起，至今仍是制约全球木薯产量的主要因素之一，因此亟需鉴定协调宿主免疫应答的基因组组分。2-氧戊二酸依赖性双加氧酶超家族（2-oxoglutarate-dependent dioxygenase superfamily，2OGD）因其在模式植物中调控特化代谢物与防御信号方面的作用而受到广泛关注，是提升木薯抗病性的潜在而尚未被充分开发的遗传资源。本研究中，RNA测序（RNA-seq）显示，木薯接种Xpm后差异表达基因（differentially expressed genes，DEGs）主要富集于丝裂原活化蛋白激酶（mitogen-activated protein kinase，MAPK）信号通路和植物-病原体互作通路。值得注意的是，在显著富集的防御相关酶家族中，2OGD超家族构成了一个主要代谢枢纽。全基因组鉴定共解析出138个Me2OGD基因家族成员，系统发育上将其归类为DOXA和DOXC两个亚家族。进化分析表明，这些基因在保持结构保守性的同时，也发生了谱系特异性分化，这一过程可能受到纯化选择驱动，以适应热带病原压力。研究进一步采用加权基因共表达网络分析（weighted gene co-expression network analysis，WGCNA）锁定了表现出感染阶段特异性表达模式的“khaki1”模块。在该网络中，Me2OGD61被鉴定为中心枢纽基因，与27个核心防御调控因子共表达，并在木薯中表现出强烈的早期诱导。整合表达谱分析与实时荧光定量逆转录PCR（quantitative real-time reverse transcription PCR，qRT-PCR）共同证实其具有关键调控潜力。最后，通过病毒诱导基因沉默（virus-induced gene silencing，VIGS）进行功能验证表明，沉默Me2OGD61会显著提高木薯对CBB的感病性，从而证实该基因为抗性正调控因子。本研究为Me2OGD基因家族提供了系统的基因组学框架，并鉴定出Me2OGD61这一可用于培育抗CBB木薯品种的高价值分子靶标。

该论文发表于《Scientia Horticulturae》，围绕木薯细菌性枯萎病的抗性分子基础展开，核心目标是从2-氧戊二酸依赖性双加氧酶（2OGD）超家族中挖掘参与木薯抗Xpm的关键基因。木薯是热带和亚热带地区重要的粮食与淀粉作物，具有耐旱、耐瘠薄等突出农业性状，但其生产长期受木薯细菌性枯萎病严重威胁。Xpm侵染后可导致水渍状斑、萎蔫与坏死，并造成显著减产。由于抗病品种培育是防控该病最经济有效的路径，因此解析宿主在病原侵染下的防御调控网络、筛选可用于分子育种的抗病基因，具有明确的理论与应用价值。已有研究表明，2OGD家族广泛参与植物激素代谢、次生代谢、防御信号转导以及表观遗传调控，但其在木薯中的组成、演化和抗病机制尚不清晰，这正构成了本研究开展的直接背景。

研究人员首先基于Xpm接种后的时序转录组数据，系统刻画木薯在病原胁迫下的全局转录重编程过程。结果显示，木薯在感染早期即发生显著转录响应，差异表达基因主要集中于植物-病原体互作、MAPK信号通路和植物激素信号转导等免疫相关途径，同时苯丙烷生物合成、黄酮生物合成、角质/栓皮质/蜡质生物合成等次生代谢途径也同步增强，提示木薯在病原识别后不仅迅速启动信号传导，也同时进行代谢流重定向，以支撑防御反应。在多类防御相关酶家族中，2OGD家族表现出突出的富集程度和较为一致的早期诱导特征，因此被优先选为深入解析对象。

在方法上，研究主要采用以下技术路线：其一，利用来源于NCBI BioProject PRJNA881631的木薯叶片Xpm接种时序RNA-seq数据，对CK、3 h、6 h、1 d、3 d和6 d等阶段进行主成分分析、差异表达分析、层次聚类和KEGG富集分析；其二，结合BLASTP与隐马尔可夫模型（HMM）在木薯参考基因组中全基因组鉴定2OGD家族成员，并开展染色体定位、复制类型、共线性、Ka/Ks、系统发育、保守基序、基因结构和顺式作用元件分析；其三，利用WGCNA筛选与感染阶段显著相关的共表达模块及枢纽基因，结合公开组织表达谱和qRT-PCR进行表达验证；其四，以SC8品种为材料，采用VIGS对候选基因Me2OGD61进行功能验证，并通过病斑面积评估其对抗病性的影响。

在研究结果方面，论文首先在“Global transcriptomic reprogramming and temporal dynamics under Xpm infection”部分表明，Xpm侵染可驱动木薯产生明显的时间依赖性转录组重塑。主成分分析显示各生物学重复聚类良好，说明数据可靠；6 d感染样本沿主成分轴偏移最显著，提示晚期感染阶段转录重编程最为剧烈。差异表达分析进一步显示，不同时间点均出现大量上调和下调基因，而表达聚类揭示部分基因在3 h或6 h达到峰值，代表早期应答模块；另一些基因则在6 d显著升高，代表晚期应答模块。

在“Functional transition from immune signaling to metabolic flux redirection”部分，研究显示木薯对Xpm的反应具有由免疫信号快速激活到代谢通量重分配的层级特征。早期差异表达基因显著富集于MAPK信号、植物-病原体互作和植物激素信号通路，说明宿主迅速感知病原并启动防御信号。与此同时，多条次生代谢通路增强，表明防御相关代谢重编程同步发生。进一步按酶学功能归类差异表达基因后发现，2OGD与细胞色素P450、UDP-糖基转移酶、酰基转移酶和过氧化物酶共同构成主要防御酶家族，其中2OGD家族因其与早期免疫信号高度同步而受到重点关注。

在“Identification of 2OGD gene family members in cassava”及后续基因组分析部分，研究人员共鉴定出138个高可信度Me2OGD成员，分布于木薯18条染色体上，且大多编码酸性蛋白，预测主要定位于细胞质。复制事件分析显示，片段复制/全基因组复制（WGD/segmental duplication）是Me2OGD家族扩张的主导力量，串联复制也有一定贡献；所有同源基因对Ka/Ks均小于1，说明该家族总体受到强烈纯化选择。系统发育分析将Me2OGD划分为DOXA和DOXC两个亚家族，其中绝大多数成员属于DOXC，并进一步分为12个功能分支，涉及赤霉素（GA）代谢、黄酮代谢、生物碱合成及其他特化代谢过程。这一结果说明木薯2OGD家族在保守的催化框架下发生了功能分化，为参与多样化防御代谢提供了基础。与拟南芥、番茄、水稻和玉米的共线性比较进一步表明，该家族在不同植物间具有较强保守性，同时也反映出双子叶与单子叶谱系间的演化差异。

在“Conserved motif and gene structure analysis of Me2OGDs”与“Cis-elements analysis in 138 Me2OGDs”部分，研究显示Me2OGD家族在保守基序和外显子-内含子结构上兼具共性与差异。多个基序在全家族中高度保守，提示其构成Fe²⁺结合与α-酮戊二酸辅因子配位所需的催化核心；而部分基序仅分布于特定分支，暗示与信号转导或次生代谢特化功能相关。启动子分析则发现，多数Me2OGD基因含有逆境响应、激素响应和病原防御相关顺式元件，如W-box、TC-rich repeats以及水杨酸（SA）和茉莉酸（JA）响应元件，这从调控层面支持其参与病原胁迫应答。

在“Protein-protein interaction network of Me2OGD proteins”与“Identification of Me2OGDs in the hub module of the co-expression network”部分，研究进一步从网络层面筛选关键成员。基于拟南芥同源信息构建的蛋白互作网络提示，Me2OGD家族可能参与植物激素代谢与免疫反应。WGCNA进一步识别出与Xpm感染早期显著正相关的khaki1模块，该模块含有最多的Me2OGD成员。通过模块内连接度排序构建枢纽子网络后，Me2OGD61与27个防御相关基因共表达，并表现出显著的早期诱导，因此被确定为重点候选基因。该筛选并非仅依据差异表达倍数，而是综合考虑了感染时序响应、模块归属及网络中心性。

在“Tissue-specific expression profiling of 138 Me2OGDs”与“Expression profiling of 138 Me2OGDs under xpmCHN11 infection and qRT-PCR validation”部分，研究利用公开转录组数据和实验验证进一步描绘了Me2OGD家族的表达图谱。不同成员在叶、贮藏根、根尖分生组织等组织中表现出明显的组织偏好性，表明其在生长发育与组织分化中承担不同功能。病原诱导表达分析显示，部分Me2OGD在6 h显著上调，反映其参与早期防御激活；另一些则在6 d增强，提示其可能作用于后期防御阶段。对8个代表性基因进行qRT-PCR验证后，Me2OGD61等基因的表达趋势与RNA-seq结果一致，增强了转录组结论的可信度。

在“Me2OGD61 positively regulates cassava resistance against bacterial blight”部分，研究通过VIGS完成了关键功能验证。Me2OGD61被沉默后，植株在形态上无明显异常，但接种XpmCHN11后叶片病斑面积显著大于阴性对照，说明该基因被抑制会降低木薯对病原的抗性。由此可见，Me2OGD61在木薯抗CBB过程中发挥正调控作用，是连接2OGD家族表达响应与实际抗病表型的关键证据。

讨论部分指出，本研究将转录组学、比较基因组学、系统发育分析和功能验证相结合，系统揭示了Me2OGD家族在木薯抗CBB中的潜在作用。研究认为，Me2OGD家族，尤其是DOXC亚家族，可能主要通过调控激素稳态和次生代谢物积累参与病害防御。Me2OGD61在系统发育上与已知参与JA代谢的JOX类成员接近，又在感染早期快速诱导并处于防御相关共表达网络中心，说明其可能是木薯在生长与防御权衡中维持JA稳态的重要节点。该研究不仅扩展了木薯功能基因组学中对2OGD家族的认识，也为解析热带作物病害抗性机制提供了新的分子线索。

研究结论可译为：综上，本研究对木薯2OGD基因家族进行了系统的基因组与功能分析。RNA-seq分析表明，2OGD在差异表达基因中占较高比例。研究共鉴定出138个Me2OGD家族成员，并证明其扩张主要由片段复制驱动，其结构特征支撑了功能分化。整合转录组分析与全基因组鉴定结果显示，特定Me2OGD成员，尤其是Me2OGD61，在Xpm侵染过程中受到动态调控。通过VIGS进行的功能验证明确表明，Me2OGD61是木薯细菌性枯萎病抗性的正调控因子。

总体而言，该论文的重要意义在于：一方面建立了木薯Me2OGD家族较为完整的基因组学资源框架；另一方面从大规模家族成员中成功锁定并验证了抗病关键基因Me2OGD61，为木薯抗CBB分子育种提供了可操作的靶标，也为后续解析2OGD介导的激素-代谢-免疫耦联机制奠定了基础。