东印度绿核菌（Ustilaginoidea virens）菌株NRRI-FSM-1的全基因组比较解析基因组演化与遗传可塑性

《Frontiers in Fungal Biology》：Comparative genomics of Eastern-Indian Ustilaginoidea virens strain NRRI-FSM-1 at whole-genome level unravels genome evolution and genetic plasticity

　　

摘要 水稻假黑穗病（Rice False Smut, RFS）由绿核菌（Ustilaginoidea virens，有性型：Villosiclava virens）引起，已成为全球水稻生产的重要威胁，导致产量下降、品质降低及市场价值受损。尽管该病害于19世纪70年代首次在印度报道，但其基因组资源仍十分有限，制约了对病原菌多样性认知和有效防控策略的研发。研究人员对来自东印度的绿核菌分离株NRRI-FSM-1进行了高质量全基因组测序与分析，并结合其他六株绿核菌进行比较基因组研究，以探讨基因组多样性、结构变异及潜在致病相关特征。组装后的NRRI-FSM-1基因组大小为36.3 Mb，包含985个支架（scaffold），N50长度为5,781,932 bp。共鉴定出328,782个变异位点，包括302,430个单核苷酸多态性（Single Nucleotide Polymorphism, SNP）、13,224个插入和13,128个缺失；同时发现5,977个简单序列重复（Simple Sequence Repeat, SSR）和9,257个蛋白编码基因，这是目前已报道的假黑穗病菌基因组中预测基因数最多的。比较基因组分析显示六株菌之间存在显著的基因组多样性，包括候选效应蛋白组差异、基因含量变化以及群体结构的地理分化，其中印度菌株间的多样性尤为显著，表明该病原菌在不同地理区域间存在广泛的基因组变异。该研究拓展了印度及全球绿核菌的病原基因组资源，揭示了其基因组演化和遗传可塑性机制，为未来的病原菌监测、致病机理研究和分子育种抗病策略提供了基础数据。

水稻是全球最重要的粮食作物之一，然而生物胁迫尤其是真菌病害严重威胁其产量与品质。水稻假黑穗病（RFS）由绿核菌（Ustilaginoidea virens）引起，自19世纪70年代在印度首次记录以来，近二十年在全球范围内的流行程度明显上升。该病害不仅造成产量损失，还因产生真菌毒素（如ustiloxins和ustilaginoidins）而带来食品安全风险。这些毒素已在动物模型中显示出肝毒性及其他不良反应，进一步凸显了病害防控的重要性。然而，现有绿核菌的基因组资源多集中在中国和日本菌株，印度及其他地区的基因组数据相对匮乏，限制了对其种群结构、遗传多样性及致病机制的全面认识。在此背景下，研究人员对来自东印度的强毒力菌株NRRI-FSM-1进行了全基因组测序，并与其他已发表的绿核菌基因组进行比较分析，旨在填补区域基因组资源的空白，揭示基因组演化与遗传可塑性的特征，从而为抗病育种和病害管理提供科学依据。本研究发表于《Frontiers in Fungal Biology》。

在技术方法方面，研究人员采用Illumina高通量测序平台对NRRI-FSM-1菌株进行全基因组测序，获得约3.97 Gb的高质量双端读段（paired-end reads），覆盖深度约为221×。组装流程包括质量控制、de novo拼接及染色体坐标分配，并利用BUSCO评估组装完整性。比较基因组分析涵盖六个绿核菌菌株（NRRI-FSM-1、UV-GVT、iJS62、IPU010、P1、UV-8b），辅以两个外群基因组（Magnaporthe oryzae和Ustilago maydis）。研究人员进行了泛基因组构建、存在-缺失变异（Presence-Absence Variation, PAV）分析、共线性和结构变异检测、重复序列注释、基因预测与功能注释、系统发育分析、次级代谢产物基因簇预测、碳水化合物活性酶（Carbohydrate-Active Enzymes, CAZyme）家族分析及分泌蛋白与候选效应蛋白识别。所有分析均采用标准化生物信息学工具与数据库，以确保结果的可比性与可重复性。

研究结果分为多个模块。在基因组组装与质量评估部分，研究人员获得了高质量的NRRI-FSM-1基因组序列，其N50超过5.7 Mb，BUSCO完整性达97.2%，显示了较高的组装质量。与其他菌株相比，NRRI-FSM-1拥有最多的预测蛋白编码基因（9,257个），暗示其可能具备更广泛的代谢与致病潜能。泛基因组与菌株特异性分析显示，六个菌株的泛基因组总大小为43 Mb，NRRI-FSM-1的存在-缺失变异数量最多，且菌株特异性基因密度最高（2.10个/Mb），提示其在基因组水平上存在显著的适应性演化迹象。共线性分析表明NRRI-FSM-1与参考基因组UV-8b之间具有高度保守的基因顺序，皮尔逊相关系数超过0.99，反映出较强的基因组结构稳定性。小变异检测共鉴定出逾32万个SNP与InDel，其中第1号染色体的变异率最高（185/Mb），第6号最低（69/Mb），反映了不同染色体区域的演化速率差异。结构基因组特征分析发现NRRI-FSM-1拥有最多的SSR（5,977个），并以三核苷酸基序为主；转座元件占总基因组的26.52%，高于多数其他菌株，其中长末端重复（Long Terminal Repeat, LTR）类反转录转座子最为丰富。比较系统发育研究显示所有绿核菌菌株聚为一支，外群物种M. oryzae与其亲缘关系较近，而U. maydis则较远，符合其分类地位。次级代谢产物分析表明NRRI-FSM-1拥有最多的次级代谢产物合成基因簇（28个），包括其他印度菌株UV-GVT所缺乏的NAPAA类基因簇，这可能与毒素合成潜力有关。CAZyme家族分析发现NRRI-FSM-1含有273个CAZyme家族，GT59为其特有，而GH24缺失，这可能影响其对宿主细胞壁的降解模式。候选效应蛋白分析显示，NRRI-FSM-1的分泌组中约有38.65%的蛋白被预测为候选效应蛋白，其中质外体效应蛋白占比最高，且包含多个已知致病相关因子（如UvHrip1和Cut1）。

在讨论部分，研究人员指出，NRRI-FSM-1的基因组展现了丰富的结构变异与遗传可塑性，这种多样性可能是区域性适应与致病力分化的基础。泛基因组中存在大量附属基因，菌株特异性基因与重复序列的扩增可能促进了环境适应与宿主侵染能力的提升。次级代谢产物与CAZyme的差异分布，暗示不同菌株在毒素合成与细胞壁降解策略上的分化，这可能在病害流行与防控中起到关键作用。此外，候选效应蛋白库的组成差异，为理解绿核菌与水稻的互作机制提供了新的视角。

结论

东印度绿核菌菌株NRRI-FSM-1的全基因组及比较基因组分析揭示了一个高度动态的基因组，其特征包括显著的基因组多样性、广泛的结构变异及丰富的遗传可塑性。附属基因组的扩张、存在-缺失变异的增加以及菌株特异性基因的积累，均指向该病原菌在区域种群层面的适应性演化。SNP、InDel、SSR及转座元件的变异进一步印证了其持续演化的趋势。功能分析显示，NRRI-FSM-1在候选效应蛋白、次级代谢产物基因簇及CAZyme家族方面均表现出增强的潜能，这可能赋予其更强的宿主侵染、毒素合成及致病能力。总体而言，本研究不仅拓展了绿核菌的基因组资源，还为解析其演化动态、致病机理及制定有效的病害管理策略提供了重要的科学基础。