香蕉种质资源抗枯萎病相关的根际细菌群落结构与功能特征

《Biology》：Rhizosphere Bacterial Community Structure and Functional Characteristics Associated with Fusarium Wilt Resistance in Banana Germplasms

　　

根际细菌群落是植物健康和抗病性的重要决定因素。为探究与香蕉枯萎病抗性相关的细菌群落特征，研究人员采用全长16S rRNA扩增子测序技术，对7个具有对比抗性水平的香蕉种质资源的根际细菌群落进行了表征。分析检测到种质资源间存在显著的基因型依赖性细菌群落组成差异。α多样性分析表明，细菌丰富度和多样性与枯萎病抗性并不一致相关。相比之下，高抗种质资源中富集了潜在有益分类群，包括芽孢杆菌属（*Bacillus*）和硝化螺旋菌属（*Nitrospira*），而Chujaibacter属和不动杆菌属（*Acinetobacter*）在中抗和低抗种质资源中丰度更高。β多样性分析进一步揭示了不同种质资源间群落结构的差异。功能预测显示，高抗种质资源相关的细菌群落具有更强的碳水化合物代谢、次生代谢物生物合成和抗生素生物合成潜力。共现网络分析揭示了不同抗性水平种质资源间细菌互作模式的差异。高抗种质资源中多个关键细菌分类群之间存在强关联，而中抗和低抗种质资源则表现出不同的正负关联模式。综合来看，这些结果表明香蕉枯萎病抗性并非与细菌多样性增加相关，而是与潜在有益细菌分类群的富集、增强的细菌功能潜力以及独特的细菌互作模式相关。

**论文解读：香蕉种质资源抗枯萎病相关的根际细菌群落结构与功能特征**

**研究背景与问题**
香蕉枯萎病（Fusarium wilt，由尖孢镰刀菌古巴专化型*Fusarium oxysporum* f. sp. *cubense*，Foc引起）是全球香蕉生产中最具破坏性的土传病害之一，尤其是热带4号生理小种（Tropical Race 4，Foc TR4）的扩散持续威胁亚洲、非洲、中东及南美洲的主要产区。由于多数栽培香蕉品种抗性有限，现有化学防治、土壤熏蒸、轮作等策略效果不稳定，且病原菌厚垣孢子可在土壤中长期存活，导致病害长期防控困难。近年来，根际微生物介导的生物防治因通过生态过程抑制病原菌而受到关注——健康土壤中常富集有益微生物类群，而感病土壤微生物网络趋于简化。植物基因型可选择性招募特定根际微生物，抗病品种往往富集具有拮抗或有益功能的类群。然而，目前对香蕉不同抗性种质资源根际细菌群落的系统比较仍有限，尤其在群落水平上细菌群落组成、功能潜力与宿主抗性之间的关系尚未被全面评估。

**研究内容与意义**
研究人员选取7个具有不同枯萎病抗性水平的香蕉种质资源（均为*Musa acuminata* AAA组卡文迪许亚群），包括低抗（Dianjiao No. 1, DJ1H）、中抗（Hongyan No. 5, HY5H；Dianjiao No. 3, DJ3H）和高抗（Dianjiao No. 4, DJ4H；Nantianhuang, NTH；Baodaojiao, BDJ；Dajiao, DJ），利用全长16S rRNA基因扩增子测序（PacBio平台）表征其根际细菌群落，比较细菌多样性、群落组成、功能潜力及共现网络特征。研究发现，抗性差异与细菌多样性无一致关联，但高抗种质资源显著富集潜在有益类群（如*Bacillus*、*Nitrospira*），并具有更强的碳水化合物代谢、次生代谢物及抗生素生物合成功能潜力，同时细菌互作网络呈现更紧密的合作关联模式。该研究揭示了香蕉枯萎病抗性并非单纯依赖细菌物种数量，而是与根际细菌群落组成、功能潜力及互作模式的整体差异相关，为基于微生物组管理的枯萎病可持续防控策略提供了理论基础。论文发表在《Biology》。

**主要关键技术方法**
研究基于云南省红河热带农业科学研究所（河口县）的香蕉种质资源圃，采集7个种质资源（各4个重复，共28份）健康植株抽穗期根际土壤。采用全长16S rRNA基因（引物F: AGRGTTTGATYNTGGCTCAG, R: TASGGHTACCTTGTTASGACTT）PacBio Sequel II平台单分子实时测序，获得环形一致性测序（CCS）序列。以97%相似度聚类为操作分类单元（OTU），基于SILVA 138数据库进行物种注释。使用QIIME2进行α多样性分析（ANOVA及Tukey检验），β多样性基于Bray-Curtis距离进行主成分分析（PCA）和非度量多维尺度分析（NMDS）及ANOSIM检验。差异分类群通过LEfSe分析（LDA>4.0, p<0.05）识别。功能预测采用PICRUSt2基于KEGG和COG数据库。共现网络基于Spearman秩相关（|r|>0.5, p<0.01）构建。

**研究结果**
**3.1 16S rRNA测序质量评估与OTU分析**
总计获得28个样品的高质量CCS序列（50,922–64,919条/样品），序列长度1447–1456 bp，稀释曲线趋于饱和，表明测序深度充足。所有种质资源共享1491个核心OTU，以变形菌门（Proteobacteria）为主，其次为厚壁菌门（Firmicutes）、拟杆菌门（Bacteroidota）、放线菌门（Actinobacteriota）和硝化螺旋菌门（Nitrospirota）；属水平上*Nitrospira*、*Bacillus*、*Sphingomonas*等为优势核心类群。高抗种质资源DJ拥有最多独有OTU（1283个），低抗种质资源DJ1H次之（917个），中抗种质资源DJ3H最少（456个）。

**3.2 根际细菌群落组成分析**
共鉴定37门、92纲、307目、681科、1589属、4283种。门水平优势菌门为变形菌门、放线菌门、酸杆菌门（Acidobacteriota）、拟杆菌门和厚壁菌门。属水平上，相对丰度差异显著：Chujaibacter在中抗种质资源DJ3H中丰度达17.31–18.76%，在高抗种质资源DJ4H中为9.11–10.14%，其他种质资源低于0.03%；*Acinetobacter*在中抗种质资源HY5H中富集（12.63–13.92%）；*Sphingomonas*在低抗种质资源DJ1H中最高（5.08–5.41%）；*Bacillus*在高抗种质资源DJ和BDJ中丰度较高（3.93–4.61%和2.15–2.29%），中抗和低抗种质资源中仅0.05–0.07%。

**3.3 α多样性分析**
ACE、Chao1、Shannon、Simpson及PD whole tree指数在不同种质资源间存在显著差异（ANOVA, p<0.05）。HY5H和BDJ具有较高丰富度；DJ和BDJ显示较高多样性（Shannon指数DJ最高为11.46–11.58）；DJ3H和HY5H多样性较低。但α多样性指标与枯萎病抗性无一致关联。

**3.4 β多样性分析**
PCA和NMDS（stress=0.084）显示不同抗性种质资源根际细菌群落结构明显分离。ANOSIM（R=0.945, p=0.001）证实种质资源间差异显著大于种内重复，表明群落结构存在基因型依赖性差异。

**3.5 差异细菌分类群与关键根际细菌类群分析**
LEfSe分析识别出种质资源特异性生物标志物：低抗种质资源DJ1H富集*Burkholderia–Caballeronia–Paraburkholderia*、*Sphingomonas*、*Rhodanobacter*；中抗种质资源DJ3H富集Chitinophagaceae、Acidobacteriales等，HY5H富集*Pseudomonas*、*Acinetobacter*；高抗种质资源DJ富集*Bacillus*及Gemmatimonadaceae，NTH富集*Nitrospira*、Myxococcaceae，BDJ富集Rhizobiales、Xanthobacteraceae。各高抗种质资源均有独特的富集谱。

**3.6 功能预测分析**
PICRUSt2预测显示，KEGG Level 2层次上，高抗种质资源DJ、BDJ、NTH的碳水化合物代谢功能丰度较高（0.0874、0.0861、0.0859）；DJ和NTH在翻译、复制与修复通路丰度较高。KEGG Level 3层次上，高抗种质资源NTH、BDJ、DJ的次生代谢物生物合成（0.0771、0.0765、0.0770）和抗生素生物合成（0.0574、0.0570、0.0576）功能丰度较高；DJ4H和BDJ在群体感应（quorum sensing）功能丰度较高。低抗种质资源DJ1H的次生代谢物和抗生素生物合成功能丰度与部分高抗种质资源相当。

**3.7 共现网络分析**
属水平共现网络显示，低抗种质资源（DJ1H）网络以*Acidibacter*、*Rhodanobacter*、Chujaibacter等为核心，存在复杂的正负相关；中抗种质资源（HY5H、DJ3H）中*Pseudomonas*、*Bacillus*、*Acidibacter*等形成多种关联，Chujaibacter与多个类群呈负相关；高抗种质资源（DJ、BDJ、NTH）中Terrimonas、Chujaibacter、*Acidibacter*、MND1、*Pedomicrobium*、*Gemmatimonas*等形成多个强正相关合作关联，如MND1与*Candidatus_Omnitrophus*正相关，*Gemmatimonas*与未分类的Acidobacteriales正相关。高抗种质资源中关键类群呈现更紧密的合作结构。

**总结讨论与结论**
讨论部分指出，根际细菌群落组成和结构受宿主基因型强烈影响，不同抗性种质资源存在独特的细菌群落配置。细菌多样性本身与抗性无一致关联，进一步支持抗性取决于特定有益类群富集、整体群落结构及互作模式的观点。高抗种质资源中*Bacillus*等有益菌的富集，以及碳水化合物代谢、次生代谢物和抗生素生物合成路径的增强，可能共同贡献于抗性。低抗种质资源虽也富集某些有益菌（如*Burkholderia–Caballeronia–Paraburkholderia*），但共现网络模式不同，提示抗性差异不只依赖单一类群存在与否。研究局限性包括：土壤类型单一（砂壤土），可能限制结论的普适性；许多关键分类群的生态功能尚不明确，需要结合培养组学、宏基因组学和功能验证实验进一步阐明。

**结论部分翻译**：本研究旨在探究香蕉种质资源中与枯萎病抗性相关的根际细菌群落特征。具有不同枯萎病抗性水平的香蕉种质资源拥有独特的根际细菌群落，这些群落在群落组成、预测功能谱和细菌关联模式上存在差异。种质资源间观察到的变异表明，宿主基因型在塑造根际细菌群落中发挥重要作用。结果表明，香蕉枯萎病抗性并非单纯与细菌多样性相关，而是与根际细菌群落组织的更广泛差异相关。这些发现为未来研究疾病相关细菌分类群的生态作用及其在基于微生物组的香蕉枯萎病管理中的潜在应用奠定了基础。