随着丙午马年的到来，农业生物技术领域对新兴植物病害的防控需求日益紧迫。在山桐子——这种在中国西南、东部及南部广泛分布的重要产油木本植物——规模化种植不断扩大的背景下，一种新的病害威胁悄然浮现。山桐子因其种子含油量高、不饱和脂肪酸比例大，在生物柴油生产和食品工业中具有广阔的开发前景，已成为贵州省特色产业发展的重要经济树种。然而，一种名为香蕉黑孢（Nigrospora musae）ST1的真菌，被首次鉴定为导致山桐子叶斑病的病原体，给其种植带来了新挑战。此类由黑孢属真菌引起的病害，已在芝麻、生姜、杨树等多种作物中频繁发生，严重影响着植物生长和作物产量。为了从根源上理解这个“新敌人”的入侵手段，为未来精准防控提供“作战地图”，Yun-Ze Chen、Yan Chen和Jing Yang等研究人员在《Journal of Fungi》期刊上发表了一项开创性研究，他们首次完成了香蕉黑孢ST1的全基因组测序与深度分析，揭开了其基因组“密码”和致病潜能的神秘面纱。
为了回答“香蕉黑孢ST1如何致病”这一核心问题，研究人员采用了一系列前沿的生物技术和生物信息学分析方法。研究首先从贵州省贵定县采集的山桐子病叶样本中分离得到病原菌株。核心的实验技术包括：利用PacBio Sequel II平台进行高精度全基因组测序，结合第二代测序数据进行纠错和组装，从而获得高质量的基因组草图。随后，研究人员通过整合从头预测、同源蛋白预测和转录组证据预测等多种方法，对基因组进行全面的组件预测，包括重复序列、蛋白质编码基因、非编码RNA和假基因。在功能注释方面，研究使用了GO、KOG、KEGG等多个数据库，并专门针对致病性相关基因，比对了PHI、DFVF、TCDB、CAZy等专业数据库。此外，利用antiSMASH工具预测了次级代谢产物生物合成基因簇，并通过信号肽预测、跨膜区分析及EffectorP工具预测了分泌蛋白和效应子。最后，研究人员还将香蕉黑孢ST1的基因组与其他已公开的黑孢属物种基因组进行了比较基因组学分析，包括平均核苷酸一致性（ANI）计算和全基因组共线性分析，以明确其分类地位和进化关系。
3. 结果
3.1. 基因组特征
通过PacBio测序平台，成功完成香蕉黑孢ST1的全基因组测序和从头组装。最终组装的基因组大小为49,259,803 bp，GC含量为56.23%，由42个无缺口的contig组成。通过二代测序reads回贴分析和BUSCO分析评估，显示基因组组装完整性和质量极高。
3.2. 序列预测
3.2.1. 重复序列、蛋白质编码基因、非编码RNA和假基因预测
基因组中重复序列占5.21%。共预测到12,063个蛋白质编码基因，其中95.39%的基因得到同源预测和转录证据支持。此外，还预测了475个非编码RNA和17个假基因。
3.2.2. 基因簇预测
通过antiSMASH工具，共预测到77个次级代谢产物生物合成基因簇，其中绝大多数（96.10%）与已知基因簇无相似性，表明该菌可能产生大量未知的次级代谢产物。鉴定出3个与已知簇100%相似的基因簇，分别与化合物alternapyrone、dimethylcoprogen和serinocyclin A/B的合成相关。
3.3. 基因组功能注释
利用GO数据库对6597个基因进行了功能分类。KOG注释显示，注释基因最多的类别是“仅具一般功能预测”，而KEGG通路注释显示，注释基因数量最多的通路包括氨基酸生物合成、碳代谢和核糖体等。
3.4. 碳水化合物活性酶
在CAZy数据库中共注释到722个碳水化合物活性酶基因。其中，糖苷水解酶（GHs）数量最多，达302个，显著高于多糖裂解酶（PLs，11个）和碳水化合物酯酶（CEs，134个）。
3.5. 致病性系统分析
3.5.1. PHI分析
在病原-宿主互作数据库（PHI）中，共注释到3600个与致病性相关的基因，占总预测基因的29.84%。其中，与毒力降低相关的基因数量最多（1413个）。
3.5.2. 真菌毒力因子数据库注释
在真菌毒力因子数据库中共注释到2503个基因，占总预测基因的20.75%。
3.5.3. 转运蛋白分类数据库注释
在转运蛋白分类数据库中共注释到125个基因，占总预测基因的1.04%。
3.5.4. 细胞色素P450基因
共预测到585个细胞色素P450基因。
3.6. 蛋白质亚细胞定位分析
预测到1462个具有信号肽的蛋白、2713个跨膜蛋白和1071个分泌蛋白，其中最终预测出197个效应子蛋白。进一步筛选出39个候选致病相关效应子，其功能域分布被可视化分析。
对CAZymes、PHI、DFVF和分泌组进行整合分析，维恩图显示有6个基因是这四类数据的交集。这6个分泌型碳水化合物活性酶基因，包括3个糖苷水解酶基因、2个果胶酸裂解酶基因和1个角质酶基因，它们被认为在植物致病过程中可能扮演重要角色。
3.7. 与其他黑孢属物种的比较基因组学分析
选取多个黑孢属物种基因组进行比较。平均核苷酸一致性分析显示，香蕉黑孢ST1与同物种菌株香蕉黑孢贵州的相似性高达97.73%，而与N. lacticolonia、N. oryzae等其他物种相似性较低（84.43–86.49%），这进一步确认了其分类地位。共线性分析表明，香蕉黑孢ST1与香蕉黑孢贵州的基因组结构保守性最高，而与其他物种间存在广泛的基因组重排，显示出显著的基因组结构差异。
4. 讨论与结论
本研究首次完成了引起山桐子叶斑病的病原真菌香蕉黑孢ST1的高质量全基因组测序与解析。基因组大小为49.26 Mb，预测了超过1.2万个蛋白质编码基因。研究发现，该基因组富含碳水化合物活性酶，尤其是大量糖苷水解酶，表明其具有强大的植物细胞壁降解能力，这是其突破植物物理防御屏障、实现侵染的关键。通过对PHI、DFVF、TCDB等致病相关数据库的注释，研究者系统鉴定了大量与毒力、转运、次生代谢相关的基因。尤为重要的是，整合分析发现了6个关键的分泌型碳水化合物活性酶，包括角质酶、木聚糖酶和果胶酸裂解酶，它们很可能是病原菌入侵宿主的“分子武器”。
次级代谢产物生物合成基因簇的分析揭示了该菌巨大的代谢潜能，其拥有的基因簇数量（77个）在所比较的黑孢属物种中最高。其中绝大多数是新簇，预示着该菌能合成大量未知的化合物，为天然产物发现提供了宝藏。鉴定出的少数已知基因簇与具有植物毒性或细胞毒性的化合物（如alternapyrone）合成相关，直接关联了其致病机制。效应子预测则找到了39个候选蛋白，其中包含多种细胞壁降解酶，进一步细化了其攻击宿主的方式。
比较基因组学分析不仅从ANI值上确认了香蕉黑孢ST1的分类地位，更揭示了黑孢属不同物种间存在显著的基因组结构重排，即使在同一物种的不同菌株间也是如此。这种基因组结构可塑性被认为是真菌适应不同宿主和生态环境的重要进化策略。
综上所述，这项研究首次绘制了香蕉黑孢ST1的完整基因组蓝图，系统揭示了其基因组结构特征、功能组成及与致病性密切相关的关键基因资源。它不仅为深入阐明香蕉黑孢ST1侵染山桐子的分子致病机制奠定了坚实的理论基础，也为未来开发针对该叶斑病的特异性分子靶点、新型杀菌剂或抗病育种策略提供了关键的数据支持和基因候选库，对保障山桐子这一重要经济树种的健康发展具有重要的科学意义和应用价值。
