硝铵比6:4驱动连作芝麻根际微生物组重塑与防御代谢物富集：基于宏基因组与代谢组的多组学解析

《Microorganisms》：Rhizosphere Soil Microbiome and Metabolome Responses of Sesame to Different Nitrate-to-Ammonium Nitrogen Ratios: A Multi-Omics Approach

　　

理解根际微生物群落和根系相关代谢物如何响应硝态氮与铵态氮比例，对于优化芝麻氮素管理至关重要。研究人员设置三个NO₃^?-N:NH₄⁺-N比例（NB1=10:0、NB2=6:4、NB4=0:10），通过宏基因组测序（metagenomic sequencing）和非靶向代谢组学（untargeted metabolomics）评估其对芝麻农艺性状、产量、微生物群落结构、功能谱和代谢组谱的影响。样品采集自四个区室：健康和患病植株的根际（H0和D0）与非根际土（bulk soil，H1和D1）。结果表明，NB2处理（6:4）获得最高种子产量（923.43 kg·ha^?1），较NB1（10:0）增加7.04%。LEfSe（linear discriminant analysis effect size，线性判别分析效应量）分析鉴定出处理间显著差异类群，其中NB2RH0含有8个差异细菌类群（s_Methylobacterium radiodurans、s_Flavisolibacter ginsenosidimutans、s_Sphingomonas sp_LK11等）和3个差异真菌类群（g_Cercospora、o_Mycosphaerellales和f_Mycosphaerellaceae）。通过KEGG（Kyoto Encyclopedia of Genes and Genomes，京都基因与基因组百科全书）、CAZy（Carbohydrate-Active enZYmes Database，碳水化合物活性酶数据库）、eggNOG（evolutionary genealogy of genes: Non-supervised Orthologous Groups，直系同源基因簇数据库）和GO（Gene Ontology，基因本体）功能注释，揭示处理特异性富集的通路，涉及过氧化物酶体（Peroxisome，map04146）、多种植物次生代谢物生物合成（Biosynthesis of various plant secondary metabolites，map00999）、EC 4.1.1.48、GH78和M00848。PLS-DA（partial least squares discriminant analysis，偏最小二乘判别分析）和单变量分析鉴定出多个与抗病性密切相关的代谢物，包括正离子模式下的Hamaudol、Nodakenin、Psoralenoside、Rhapontigenin和Columbianetin acetate，其在健康植株根际土壤中显著更丰富。负离子模式下，Pyranonigrin A、Pteroic acid和Botryaloic acid acetate在健康植株根际土壤中的丰度显著低于患病植株。总体而言，这些发现表明这些氮形态比例显著塑造了根际微生物群落、功能谱和代谢组谱。本研究为通过氮形态管理改善芝麻种植土壤健康提供了有价值的见解。

研究背景与问题：硝态氮（NO₃^?）与铵态氮（NH₄⁺）比例是土壤微生物群落的重要调控因子，影响生态系统功能和农业生产力。芝麻（Sesamum indicum L.）是重要特色油料作物，连作障碍严重，茎枯病和青枯病等土传病害缺乏有效化学药剂，农户连作导致病害加重。缓解措施包括减施氮肥、生物炭改良和优化硝铵比。已有研究多基于16S rDNA扩增子测序，宏基因组和代谢组层面的机制仍不清楚。研究人员假设最佳硝铵比通过重塑根际微生物群、促使微生物平衡由病原向有益类群转变及改变关键土壤代谢谱来缓解芝麻连作障碍，整合微生物组-代谢组相关而非单一独立因素驱动有效抑病。研究在江西南昌进贤县罗溪镇2024年生长季开展，红壤旱地pH 5.8、有机碳1.8%、阳离子交换量12.4 cmol(+) kg^?1、初始矿质氮32.5 kg ha^?1；采用随机完全区组设计，3个氮形态比例、4次重复、12个小区，总氮75 kg N ha^?1，NB1为100% NO₃^?-N，NB2为60% NO₃^?-N+40% NH₄⁺-N，NB4为100% NH₄⁺-N。芝麻品种Yuzhi No. ND837连续单作两年诱导自然青枯病，第三年成熟期采集健康植株根际（RH0）、健康区非根际土（RH1）、患病植株根际（RD0）和患病区非根际土（RD1）。产量以单因素方差分析（ANOVA）比较，土壤养分、微生物群落和代谢物按12组比较。研究得出NB2产量最高（923.43 kg·ha^?1），较NB1增加7.04%，但农艺性状和产量组分影响较小；氮形态比例显著塑造微生物群落结构、功能谱和代谢组，健康植株根际正离子模式多种抗病相关代谢物更高，负离子模式则相反。该研究为氮形态管理改善土壤健康和缓解连作障碍提供依据。

关键技术方法（不超过250字）：2024年江西进贤县田间试验设置3个NO₃^?-N:NH₄⁺-N比例（10:0、6:4、0:10），总氮75 kg N ha^?1，4重复；芝麻品种Yuzhi No. ND837连续单作两年，第三年成熟期采集健康/患病植株根际与非根际土共12组。采用Illumina PE150宏基因组测序，Kraken2/Bracken进行物种注释，HUMAnN2/UniRef90进行功能注释，CARD注释抗性基因，结合LEfSe、DESeq2、PCoA/NMDS、RDA和Spearman相关分析；非靶向代谢组学采用UPLC-QTOF-MS正负离子模式，XCMS处理，MetaboAnalystR进行PLS-DA、t检验/折叠变化、通路分析和随机森林。样本队列来源为上述田间试验，序列提交NCBI-SRA PRJNA1521067。

研究结果：3.1 芝麻农艺性状和种子产量的影响：田间产量和农艺性状测定及ANOVA显示，NB4结实穗长较NB2增加3.4 cm，产量组分无显著差异；NB2种子产量最高，较NB1增加7.04%。土壤化学性质中仅pH、有效磷（AP）和有效钾（AK）处理间显著，NB1下位置间pH差异显著，NB1RD0较NB1RD1和NB1RH1分别增加5.68%和4.85%；NB1RD1的AP较NB1RH0和NB1RH1分别高64.32%和62.64%。3.2 微生物群落结构的影响：宏基因组物种注释、NMDS和PERMANOVA显示细菌（F=1.481，p=0.001）和真菌（F=1.241，p=0.034）群落受处理显著影响，古菌和病毒无显著。LEfSe显示NB2RD0差异细菌最多（27个），NB2RH0含8个差异细菌（s_Methylobacterium radiodurans、g_Flavisolibacter、s_Flavisolibacter ginsenosidimutans、s_Sphingomonas sp LK11等）和3个差异真菌（g_Cercospora、o_Mycosphaerellales、f_Mycosphaerellaceae）；NB4RH0差异真菌最多（6个）。三因素ANOVA仅少数物种受单一因素显著影响。3.3 细菌群落功能的影响：KEGG、CAZy、eggNOG和GO注释及LEfSe/DESeq2显示，KEGG Level 2有4个处理显著，Level 3中NB2RH0富集map04146和map00999；模块水平NB2RH0含M00345、M00848、M00344和M00052；EC水平NB2RH0含EC 4.1.1.48等；CAZy中NB2RH0含GH78和GH140。差异功能潜力分析共鉴定380个显著项，下降320个、增加60个；NB2RH0与NB2RD0热图明显分离，核糖体等通路富集，PPI呈核心-外围结构。3.4 代谢组谱的影响：UPLC-QTOF-MS正负离子模式PCA、PLS-DA和单变量分析显示，正离子模式鉴定14个显著差异代谢物，单变量分析15个，其中7个与抗病相关，Hamaudol、Nodakenin、Psoralenoside、Rhapontigenin和Columbianetin acetate在健康植株根际更高，4’-Methoxychalcone和Lobaric acid例外；负离子模式鉴定11个显著差异代谢物，单变量分析15个，6个与抗病相关，Pyranonigrin A、Pteroic acid和Botryaloic acid acetate在健康植株根际显著低于患病植株，另3个相反。三因素ANOVA显示Hamaudol受Status和Location×Status影响，Nodakenin受Location、Status及交互影响，Psoralenoside受Location、Status和Location×Status影响；负离子模式三个代谢物均受Status和Location×Status显著影响。3.5 微生物类群和功能通路对环境变量的响应：RDA和Spearman分析显示环境因子显著塑造群落结构（p=0.001）和功能组成（p=0.015），RDA1解释功能变异78.77%、群落结构27.14%；pH和AK为主要驱动因子，AN和AB与RDA2相关；多数细菌属与AK、pH、AP显著正相关，pH、AP、AK与多数代谢功能强正相关，主要富集于能量代谢和分泌系统。

讨论总结：讨论认为6:4比例提高连作芝麻产量，机制涉及缓冲根际酸化、铵态氮适度刺激局部根系增殖及碳氮代谢协同；连续种植条件可能放大其产量优势。关键微生物方面，s_Methylobacterium radiodurans、g_Flavisolibacter和s_Sphingomonas ginsengisol等被视为氮形态驱动生态位选择的潜在指示类群，其中Sphingomonas的促生和抑病机制支持最强；Cercospora富集被解释为真菌功能与生境过滤结果，而非简单病原风险指标。功能通路方面，map04146过氧化物酶体参与活性氧（ROS）稳态和氧化还原防御，map00999涉及香豆素、呋喃香豆素、苯并噁嗪酮等防御性次生代谢物生物合成；EC 4.1.1.48（indole-3-glycerol phosphate synthase，IGPS）连接色氨酸和苯并噁嗪酮防御，GH78与植物细胞壁降解和致病性相关，M00848与aurachin类抗生素生物合成相关。代谢物方面，正离子模式7个抗病相关代谢物构成多层防御模块，香豆素类、查尔酮类和芪类分别参与微生物组调控、直接抑菌和群体感应抑制；健康土中5个更高被解释为主动保护状态，4’-Methoxychalcone和Lobaric acid的例外可能与胁迫诱导、周转或微生物转化有关。负离子模式中Pyranonigrin A与真菌氧化胁迫适应相关，Pteroic acid和Botryaloic acid acetate被视为候选病害相关代谢标志；健康土代谢组偏向微生物支持与抑病，病土积累胁迫适应、真菌持久或失调相关代谢物。整合而言，6:4比例通过优化碳氮平衡和根际pH、重塑微生物群落并促进防御代谢物积累，多层面缓解连作障碍；因未检验交互项，差异应解释为组水平关联。研究结论（翻译）：本研究考察硝铵比对芝麻根际微生物群和代谢物的影响，6:4比例产量最高；氮形态比例明显塑造微生物群落结构、功能谱和代谢网络，正离子模式下若干抗病相关代谢物在健康植株根际富集，负离子模式则相反。这些发现强调氮形态管理在调控根际生态方面比直接促进产量作用更大。未来应通过合成菌群（SynCom）实验和基因敲除方法验证关键差异代谢物和关键类群的功能，并评估优化氮比例在不同田间条件和芝麻品种下对土壤健康和抑病的长期影响。