《Tropical Plant Pathology》：Comprehensive pangenomic, phylogenetic, and genome comparison analyses demonstrate that the causal agent of castor bean (Ricinus communis L.) bacterial blight is a pathovar of the Xanthomonas euvesicatoria species
编辑推荐：
为厘清蓖麻细菌性疫病病原（原命名Xanthomonas axonopodis pv. ricini）的分类学疑点，本研究对其模式菌株ICMP 5747进行了基因组测序，并综合运用多基因序列分析（MLSA）、系统发育基因组学、泛基因组学及基因组序列比较分析。结果表明，该病原实为Xanthomonas euvesicatoria的一个致病变种，并提议将其与X. axonopodis pv. commiphoreae一并转移至X. euvesicatoria物种下。此研究为明晰黄单胞菌属的系统发育关系及植物病害病原学确定提供了关键依据。
蓖麻，这种以生产蓖麻油而闻名于医药、化妆品及生物柴油等工业领域的重要经济作物，一直受到一种名为细菌性叶疫病（Bacterial Leaf Blight, BLB）的严重威胁。自1928年在日本首次被描述以来，导致该病害的细菌病原历经数次分类学上的更名，从最初的Bacterium ricinicola，到Xanthomonas campestris pv. ricini，再到目前最广为接受的名称Xanthomonas axonopodis pv. ricini。然而，其分类学地位始终模糊不清。近年来的多项多基因序列分析（Multilocus Sequence Analysis, MLSA）研究发现，被标记为此病原的某些菌株并未与X. axonopodis聚集，而是与Xanthomonas euvesicatoria的菌株关系更近。更令人警觉的是，一项2024年的研究揭示，其模式菌株ICMP 5747不仅能在桉树（Eucalyptus）上致病，而且在系统发育分析中也与X. euvesicatoria聚为一类，这对其宿主专一性和分类归属都提出了根本性质疑。准确鉴定病原的分类学身份，是理解其致病机理、宿主范围、制定有效防控策略乃至进行可靠检疫检测的基石。为了彻底澄清这一悬而未决的科学问题，一个国际研究团队在《Tropical Plant Pathology》期刊上发表了一项研究，旨在通过最先进的基因组学手段，为这位“身份不明”的植物病原正名。
为了达成研究目标，作者团队采用了多层次、互补的分析策略。首先，他们对目标菌株——X. axonopodis pv. ricini的模式菌株ICMP 5747进行了Illumina NovaSeq 6000平台测序，组装获得了一个5.1 Mb的高质量草图基因组。随后，他们开展了四项核心分析：1）基于四个看家基因（dnaK, fyuA, gyrB, rpoD）部分序列的多基因序列分析（MLSA），分别采用贝叶斯推断和最大似然法构建系统发育树；2）系统发育基因组学分析，利用Cognac软件包从41个黄单胞菌（包含Xylella fastidiosa作为外群）基因组中鉴定出368个保守同源基因，并以此构建最大似然树；3）泛基因组学分析，使用Roary比较ICMP 5747与X. euvesicatoria、X. axonopodis多个菌株的基因组成；4）基因组比较，计算平均核苷酸一致性（Average Nucleotide Identity, ANI）和数字化DNA-DNA杂交（digital DNA-DNA Hybridization, dDDH）值，并使用Mauve进行基因组共线性比对，以评估菌株间的基因组相似性和结构保守性。
研究结果
ICMP 5747模式菌株基因组序列的特征
组装获得的ICMP 5747基因组大小为5.1 Mb，GC含量为65.0%，完整度达99.7%。注释预测了4,201个编码序列和65个RNA基因。功能分析显示，其基因组编码了与细菌基础代谢、致病性（如III型和VI型分泌系统）及抗性相关的众多基因。
MLSA将ICMP 5747菌株与X. euvesicatoria模式菌株归为一组
基于dnaK, fyuA, gyrB, rpoD四个基因部分序列构建的系统发育树（贝叶斯推断和最大似然法结果一致）清晰显示，菌株ICMP 5747与X. euvesicatoria（包括其致病变种X. euvesicatoria pv. perforans的模式菌株）聚集在一个高支持率的进化支中。而X. axonopodis则位于另一个与之明显分歧的进化支。这初步表明ICMP 5747在系统发育上更接近于X. euvesicatoria。
系统发育基因组学将菌株ICMP 5747划归X. euvesicatoria物种
为了获得更可靠的证据，研究基于368个在所有分析菌株中保守的同源基因构建了系统发育基因组树。结果再次证实，ICMP 5747菌株位于X. euvesicatoria进化支内部，而X. axonopodis模式菌株则形成一个独立的单系群，与ICMP 5747明显分开。这从全基因组水平强有力地支持了其分类学归属。
菌株ICMP 5747的基因组序列与X. euvesicatoria更为相似
基因组相似性指数计算提供了定量证据。ICMP 5747与X. euvesicatoria致病变种（如pv. euvesicatoria和pv. perforans）的ANI值高于97%（覆盖度>86%），dDDH值高于77%，均超过了公认的物种界定阈值（ANI >95-96%， dDDH >70%）。相反，其与X. axonopodis模式菌株的ANI值仅为92.41%，dDDH值仅为52.0%，远低于物种阈值。这从基因组序列相似性角度确证了其与X. euvesicatoria属于同一物种。
X. axonopodis致病变种ricini的基因组与X. euvesicatoria模式菌株更为相似
进一步的泛基因组学分析显示，ICMP 5747的基因组在基因组成上与X. euvesicatoria菌株聚为一簇，而与X. axonopodis菌株分离。通过Mauve进行的基因组共线性比对发现，ICMP 5747的基因组与X. euvesicatoria模式菌株ATCC 11633的基因组之间，存在更多、更大的局部共线性区块，基因组重排（倒位、易位）更少，表现出更高的同线性。这从基因组结构和组织层面证明了其与X. euvesicatoria的紧密关系。
意外发现：X. axonopodis pv. commiphoreae也属于X. euvesicatoria
在研究过程中，一个意外的发现是，另一个致病变种X. axonopodis pv. commiphoreae的模式菌株LMG 26789，在泛基因组学分析中也与X. euvesicatoria菌株聚类。计算其ANI和dDDH值，与X. euvesicatoria菌株的相似性也远高于与其他黄单胞菌物种的相似性，且基因组共线性分析也得出了与ICMP 5747一致的结论。这表明LMG 26789同样被错误地归类于X. axonopodis之下。
结论与讨论
本研究通过整合多基因序列分析、系统发育基因组学、泛基因组学和全基因组比较等多种前沿生物信息学方法，提供了确凿无疑的证据，表明长期以来被认为是蓖麻细菌性叶疫病病原的Xanthomonas axonopodis pv. ricini，其正确分类学位置应为Xanthomonas euvesicatoria。基于同样坚实的证据，研究同时提出，Xanthomonas axonopodis pv. commiphoreae也应被重新划归为Xanthomonas euvesicatoria的一个致病变种。
因此，作者正式提议进行如下分类学转移：将Xanthomonas axonopodis pv. ricini转移为Xanthomonas euvesicatoria pv. ricini comb. nov.（新组合），将Xanthomonas axonopodis pv. commiphoreae转移为Xanthomonas euvesicatoria pv. commiphoreae comb. nov.（新组合），并对X. euvesicatoria的物种描述进行了相应修订。
这项研究的意义重大而深远。首先，它澄清了一个长期存在的分类学疑点，有助于建立更准确、反映真实进化关系的黄单胞菌属系统发育框架。其次，正确的分类是病害管理的基础。将病原准确归入X. euvesicatoria，提示其宿主范围可能不仅限于蓖麻（研究已证实其可侵染桉树），这为评估其潜在的生态风险和制定跨作物的防控策略提供了关键信息。例如，广泛分布于热带地区的野生或入侵蓖麻植株，可能作为病原的“蓄水池”，威胁商业化种植的桉树林及其他可能感病的经济作物。此外，该研究为未来开发针对X. euvesicatoria pv. ricini的特异性分子检测方法（如基于病原特异性序列的PCR检测）奠定了基因组学基础，对于病害监测、种苗检疫和早期预警至关重要。最后，该研究也再次表明，尽管近年来黄单胞菌属的分类已有不少修订，但仍有一些具有经济重要性的植物病原分类地位不清，需要运用基因组学工具进行持续而细致的梳理，以服务于植物病理学和农业生产的实际需求。
