尖孢镰刀菌古巴专化型过氧化物酶的全基因组分析揭示氧化应激反应与感染相关表达

《Journal of Fungi》：Genome-Wide Analysis of Fusarium oxysporum f. sp. cubense Peroxidases Reveals Oxidative-Stress Responses and Infection-Associated Expression

　　

香蕉枯萎病由土传真菌尖孢镰刀菌古巴专化型（Fusarium oxysporum f. sp. cubense，Foc）引起，是一种毁灭性病害，对全球香蕉生产构成重大威胁。研究人员在Foc 4号生理小种（Foc4）中鉴定出24个预测编码依赖血红素的过氧化物酶或过氧化氢酶–过氧化物酶相关蛋白的基因，同时单独编目了9个编码依赖硫醇的过氧化物还原蛋白的基因。系统发育、共线性及Ka/Ks分析表明，该过氧化物酶 repertoire在强纯化选择下高度保守，未发生显著的谱系特异性扩增。基序、结构域、基因结构及启动子分析揭示了亚组特异性特征以及大量推定的胁迫和激素响应顺式调控基序。香蕉侵染及H₂O₂处理期间的表达谱显示了不同的时间和氧化应激响应。基因本体论（Gene Ontology，GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）分析将这些基因的亚群与抗氧化活性、过氧化物代谢、过氧化物酶体功能及应激信号传导相关联。基于尖孢镰刀菌番茄专化型（F. oxysporum f. sp. lycopersici）直系同源物的相互作用网络结合共表达分析，优先筛选出5个可能与致病性相关表达程序有关的过氧化物酶。FoCP是根据其对H₂O₂的快速响应及过氧化氢酶–过氧化物酶注释单独选出的，其在酵母中异源表达后增强了氧化应激耐受性。总体而言，Foc4似乎通过条件特异性调控保守的过氧化物酶 repertoire而非基因家族扩增来适应氧化应激。网络优先筛选出的过氧化物酶为后续功能和致病性研究提供了候选基因。

香蕉枯萎病是香蕉产业最具毁灭性的土传病害之一，其病原菌尖孢镰刀菌古巴专化型4号生理小种（Fusarium oxysporum f. sp. cubense race 4，Foc4）寄主范围广、土壤存活时间长，且缺乏有效化学防治手段和抗性品种。在植物与病原菌互作过程中，活性氧（reactive oxygen species，ROS）既是植物防御信号分子，也是病原菌必须克服的屏障。过氧化物酶作为真菌抗氧化系统的核心酶家族，在ROS解毒、氧化还原信号传导和寄主-病原互作中发挥多重功能。然而，Foc4中过氧化物酶基因家族的系统性组成、进化特征及其在氧化应激适应和侵染过程中的作用机制尚不明确。研究人员对此开展了全基因组水平的鉴定和功能表征研究。

研究人员利用NCBI公开的Foc4菌株II5基因组和侵染香蕉根部的转录组数据，结合自行构建的中国Foc4分离株H₂O₂胁迫RNA-seq数据，通过保守结构域搜索鉴定了Foc4过氧化物酶基因家族成员，并综合运用系统发育分析、共线性分析、Ka/Ks选择压力分析、启动子顺式元件分析、侵染和H₂O₂胁迫下的表达谱分析、GO和KEGG富集分析、蛋白质互作网络（protein-protein interaction，PPI）构建、共表达分析以及酵母异源表达等手段，系统解析了该基因家族的进化保守性、表达调控分化和功能特性。

主要研究结果如下：

（1）Foc4过氧化物酶基因家族的鉴定与基本特征。通过Pfam HMM搜索并结合NCBI CDD和Peroxiscan验证，研究人员在Foc4基因组中鉴定到24个编码依赖血红素的过氧化物酶或过氧化氢酶–过氧化物酶相关蛋白的基因，包含Peroxidase_2家族、DyP型过氧化物酶（EfeB/YfeX-like）、抗坏血酸过氧化物酶（ascorbate peroxidase，APX）、植物过氧化物酶样超家族以及过氧化氢酶相关蛋白（包括KatG）。另外单独鉴定到9个非血红素、依赖硫醇的过氧化物还原蛋白（8个过氧化物氧还蛋白和1个BtuE样过氧化物酶）。这些蛋白在长度和亚细胞定位预测上差异显著，反映了家族的结构多样性。

（2）系统发育和进化分析。利用Foc4、尖孢镰刀菌番茄专化型（Fol）、禾谷镰刀菌（Fusarium graminearum）和稻瘟菌（Magnaporthe oryzae）的过氧化物酶蛋白序列构建最大似然树，将24个Foc4成员划分为4个系统发育组（Group I–IV）。Group I为经典血红素过氧化物酶核心谱系，Group II包含过氧化氢酶相关蛋白，Group III包含植物过氧化物酶样蛋白和APX，Group IV为KatG型过氧化氢酶–过氧化物酶。Foc4与Fol之间存在24对推定直系同源对，Ka/Ks分析显示所有配对Ka/Ks值均低于0.3（均值0.081），表明这些基因受到强纯化选择。四个真菌物种的过氧化物酶基因数量为17–28个，家族大小相对稳定，未发生显著谱系特异性扩增。

（3）结构、染色体定位和共线性分析。保守基序分析揭示了组特异性分布模式，如Group I全部含基序3，Group II全部含基序4和9，Group IV共享基序5、8、10；Group III的基序组成变异较大，提示其结构多样性更高。结构域分析显示同一亚家族成员具有高度相似的结构域架构，但外显子-内含子组织在各基因间差异明显。24个过氧化物酶基因分散分布在10条染色体上，彼此距离≥33 kb，未发现串联重复。共线性分析显示，致病性Foc4与Fol之间过氧化物酶基因所在区域高度保守，而非致病性Fo47菌株则呈现频繁的基因组重排。

（4）启动子顺式元件分析。对24个基因上游2 kb启动子区共鉴定到695个推定调控基序，分为激素响应、生长发育、胁迫响应和代谢调控四类。其中光响应元件最多（214个），MeJA响应（190个）和ABA响应（58个）元件在激素相关元件中占主导，并发现43个低氧相关元件（包括厌氧诱导和缺氧诱导）。该结果提示Foc4过氧化物酶表达可能受多种环境和寄主信号调控，但由于PlantCARE主要为植物数据库，这些基序需实验验证。

（5）侵染过程中的动态表达分析。18个在侵染转录组中可检测的过氧化物酶基因在18、32和56小时侵染时间点呈现三种表达模式：早期高表达簇（6个基因，在32 h达峰）、晚期高表达簇（8个基因，在56 h显著上调）和低表达簇（4个基因，FPKM<10）。这种时间上的转录特化表明不同家族成员在侵染不同阶段可能发挥不同作用。

（6）H₂O₂胁迫下的表达响应及qRT-PCR验证。研究人员构建了1.5 mM H₂O₂处理的中国Foc4分离株转录组数据。在8个呈现明显响应模式的基因中，FoxII5_1335143和FoxII5_16016在0.5 h内快速上调（分别达2.3倍和2.2倍）；FoxII5_19678、FoxII5_10492和FoxII5_6369在2–8 h持续上调，其中FoxII5_19678的FPKM从0 h的3.5升至8 h的15.7；FoxII5_19831和FoxII5_9926被下调。qRT-PCR验证结果与RNA-seq趋势一致。通过RefFinder评估，肌动蛋白（actin）被选为最稳定的内参基因。

（7）FoCP异源表达增强酵母H₂O₂耐受性。FoxII5_1335143（命名为FoCP）编码双功能过氧化氢酶–过氧化物酶，在H₂O₂处理后快速诱导表达。将FoCP全长编码序列克隆至pYES2-NTB载体并在酿酒酵母（Saccharomyces cerevisiae）INVSc1中受GAL1启动子驱动表达。在含3 mM H₂O₂的SG/-Ura平板上，空载体对照在所有稀释度均无可见生长，而两个FoCP表达转化子在10⁰–10^-3稀释度均能生长，表明FoCP赋予酵母H₂O₂耐受性。

（8）GO和KEGG富集分析。18个侵染相关过氧化物酶基因显著富集于过氧化物酶活性、以过氧化物为受体的氧化还原酶活性和抗氧化活性等分子功能条目，以及血红素结合、四吡咯结合和辅因子结合条目。基于9个具有KEGG直系同源注释的基因的KEGG分析显示，显著富集于色氨酸代谢、过氧化物酶体、乙醛酸和二羧酸代谢、MAPK信号通路以及碳代谢和次生代谢物生物合成等通路。KEGG结果主要反映过氧化氢酶和过氧化氢酶–过氧化物酶基因，不应推广至整个家族。

（9）PPI网络和共表达分析。以Fol为参考物种构建STRING PPI网络，鉴定出8个中心节点，其中5个对应的Foc4基因在香蕉侵染中表达：FoxII5_19831、FoxII5_10992、FoxII5_11849、FoxII5_16016和FoxII5_6369。将这5个基因与23个来自PHI数据库的致病性相关基因（涉及代谢酶、信号通路、转录调控、细胞壁功能、分泌蛋白/效应子和毒素等类别）进行Pearson相关分析，发现FoxII5_16016与22个PHI基因强正相关（r≥0.7），FoxII5_6369与全部23个PHI基因强正相关，而FoxII5_19831、FoxII5_10992和FoxII5_11849与多数PHI基因强负相关。这些结果提示不同过氧化物酶成员可能与不同的致病性相关表达程序相关联。

讨论部分指出，Foc4过氧化物酶家族的保守基序分布和强纯化选择表明其核心功能高度保守，但侵染和H₂O₂胁迫下的差异表达模式以及启动子元件的多样性暗示转录调控分化可能比基因家族扩增对功能多样化贡献更大。共线性分析结果符合“双速”基因组模型，即过氧化物酶基因位于保守核心基因组区域，但通过调控变异实现功能多样化。FoCP的异源表达实验支持其氧化应激耐受功能，但未检测酵母中FoCP转录或蛋白水平，且直接酶活性未测定。网络优先筛选的5个过氧化物酶与致病性相关基因的共表达关联属于探索性描述，仍需遗传学验证。

研究结论为：该研究首次对Foc4过氧化物酶基因家族进行了全基因组表征，证明该家族进化保守、受强纯化选择、谱系特异性扩增有限。整合表达、功能富集和网络分析揭示了Foc4过氧化物酶的调控多样性，并鉴定出多个可能与氧化应激响应和致病性相关表达程序有关的核心基因。这些发现提示Foc4可能主要通过保守抗氧化 repertoire的调控分化而非基因家族扩增来应对寄主来源的氧化胁迫，FoCP异源表达实验表明其可增强酵母的氧化应激耐受性。网络分析鉴定出的核心过氧化物酶为后续功能验证和病害防控研究提供了有前景的候选靶标。