《Frontiers in Plant Science》：Salt stress shapes rhizosphere bacterial communities and is associated with salt tolerance in tiger nut (Cyperus esculentus L.)
编辑推荐：
背景：尽管土壤盐渍化严重制约作物生产，但植物可通过与根际微生物的互作增强其耐盐性。油莎豆（Cyperus esculentus L.）是一种具有较强抗逆性的粮油兼用型植物，然而其耐盐机制及其与根际微生物的关联尚不明确。
方法：本研究开展盆栽试验，设置六个NaC
背景：尽管土壤盐渍化严重制约作物生产，但植物可通过与根际微生物的互作增强其耐盐性。油莎豆（Cyperus esculentus L.）是一种具有较强抗逆性的粮油兼用型植物，然而其耐盐机制及其与根际微生物的关联尚不明确。
方法：本研究开展盆栽试验，设置六个NaCl盐浓度梯度（0、0.36、0.72、1.08、1.44和1.80 M）。通过生理生化分析、16S rRNA测序、PICRUSt2功能预测及偏最小二乘结构方程模型（PLS-SEM），研究人员探究了盐胁迫梯度下油莎豆生理响应与根际细菌群落演替之间的耦合关系。
结果：结果表明，盐胁迫显著提高了油莎豆可溶性糖（SS）、脯氨酸（Pro）和可溶性蛋白（SP）含量，增强了超氧化物歧化酶（SOD）和过氧化氢酶（CAT）活性（P < 0.05），并与生物量分配向块茎发育方向转变相关。盐度是驱动根际细菌群落定向重塑的首要因素。变形菌门（Proteobacteria）和厚壁菌门（Firmicutes）的相对丰度显著增加，且Planococcus和Paracoccus显著富集。方差分解分析（VPA）表明，植物渗透调节因子解释了根际细菌群落结构变异的最大独立比例。相关性分析结果显示，预测的与耐盐性相关的细菌功能通路与土壤碳、氮循环功能以及植物抗逆性状和块茎干重显著正相关（P < 0.05）。PLS-SEM分析结果表明，植物渗透调节物质和预测的细菌代谢功能在连接盐胁迫与细菌功能、以及潜在耐盐根际细菌与植物抗氧化系统的通路中分别表现出统计学中介关联。
结论：综上所述，本研究提示油莎豆在盐胁迫下植物渗透调节与根际细菌富集之间存在潜在协同关联。特别是，植物渗透调节物质的积累与具有报道耐盐潜力的根际细菌属的富集相关，且这些富集细菌与植物抗逆性状密切相关，可能通过其耐盐代谢途径发挥作用。这些发现为植物渗透调节与具有报道抗逆潜力的根际细菌属富集之间的关联提供了相关性证据，该关联可能与油莎豆的耐盐性有关。然而，仍需根际分泌物分析或微生物接种实验的直接证据来证实这些关联。
油莎豆耐盐性与根际微生物协同适应机制研究解读
一、研究背景与科学问题
土壤盐渍化是全球农业生产面临的最严峻环境挑战之一，目前全球盐渍化土壤面积已超过13亿公顷，占全球陆地面积的10.7%，每年因盐胁迫废弃的农田超过500万公顷。预计到2050年，全球近50%的耕地将受到土壤盐渍化影响。盐分在土壤中的积累不仅导致土壤理化性质恶化、酶活性降低和肥力下降，还会通过离子毒害和渗透胁迫抑制植物生长发育。尽管植物的耐盐适应性常被归因于遗传分化，但根际微生物作为植物的"第二基因组"，在植物适应逆境环境中发挥着关键作用。植物促生根际细菌（PGPR）可通过调节离子稳态、分泌植物激素、促进细胞渗透物质积累、提高抗氧化活性和增强养分吸收等多种途径促进植物的耐盐碱性。
油莎豆（Cyperus esculentus L.）是莎草科莎草属一年生草本植物，具有根系发达、抗逆性强、生物产量高等特点，对我国盐碱地的开发利用具有重要战略意义。已有研究表明油莎豆耐盐性存在品种间差异，且其根际细菌群落组成受盐胁迫显著影响。然而，现有研究多局限于静态生理指标测定和单一促生菌株的筛选验证，缺乏对油莎豆在连续盐度梯度下根际微生物群落动态组装规律、介导根际微生物分化的核心驱动因素，以及不同盐度下特征根际微生物和核心菌群对油莎豆耐盐性状贡献的系统分析。
基于此，本研究提出三个科学假说：（1）盐度梯度变化影响油莎豆叶片中脯氨酸和可溶性糖等渗透调节物质的积累；（2）渗透调节物质是驱动油莎豆根际微生物群落结构差异的关键因素；（3）盐胁迫可引起油莎豆根际微生物的代谢重编程，且这些重编程的代谢通路与油莎豆的抗逆性状显著相关。
二、研究设计与技术方法
研究人员在内蒙古农牧业科学院试验站（40°45′ N，111°40′ E，海拔1040 m）开展盆栽试验。供试油莎豆品种为"西班牙1号"，栽培于装有13 kg风干沙土的盆中（盆径27.5 cm，底径22 cm，盆高31 cm）。试验设置6个NaCl浓度处理：0 M（CK）、0.36 M（CC1）、0.72 M（CC2）、1.08 M（CC3）、1.44 M（CC4）和1.80 M（CC5），每个处理3个重复，共18盆。每盆施加2 L对应浓度NaCl溶液，折算后土壤最终盐浓度为0、3.24、6.48、9.72、12.96和16.2 g kg-1，覆盖内蒙古盐渍化土壤的主要盐度范围。盐溶液分10次等量（200 mL）分步施加以避免渗透冲击。
主要技术方法包括：（1）植物生理生化测定：采用紫外分光光度法测定CAT活性，NBT光化学还原法测定SOD活性，愈创木酚比色法测定POD活性，硫代巴比妥酸（TBA）比色法测定MDA含量，蒽酮比色法测定SS含量，考马斯亮蓝染色法测定SP含量，茚三酮比色法测定Pro含量；（2）土壤指标测定：土壤过氧化氢酶（S_CAT）、脲酶（URE）、碱性磷酸酶（ALP）活性及微生物生物量碳（MBC）、氮（MBN）、磷（MBP）；（3）微生物组分析：采用DNeasy PowerSoil? Pro Kit提取根际土壤总DNA，扩增细菌16S rRNA基因V3-V4区（引物338F/806R），在Illumina MiSeq PE300平台高通量测序，共获得1,115,672条原始读段，质控后保留965,427条，按97%相似度聚类获得5,969个OTU，稀释后保留1,326个OTU，测序覆盖率99.54%；（4）生物信息学分析：使用PICRUSt2预测KEGG同源基因（KO）丰度，并将KO归入渗透胁迫（KO_OS）、抗氧化胁迫（KO_AS）、离子转运（KO_IT）和氮代谢（KO_NIr）功能类别；采用FAPROTAX预测碳、氮、硫循环功能；构建偏最小二乘结构方程模型（PLS-SEM）验证中介效应。
三、研究结果
3.1 油莎豆对不同盐胁迫水平的响应特征
盐胁迫显著抑制油莎豆营养生长，而中度盐胁迫与生物量分配向块茎发育方向转变相关。与对照（CK）相比，各盐处理下株高（PH）和根干重（RDW）均有不同程度降低，其中CC5处理下株高最低（21.0 cm），较CK显著降低50.44%。块茎干重（TDW）在盐胁迫下不同程度增加，CC3处理下达到最大值0.45 g，较CK显著增加41.80%。渗透调节方面，脯氨酸、可溶性糖和可溶性蛋白含量均在CC3处理下达到最高，分别为57.55 μg·g-1 FW、64.04 mg·g-1 FW和101.94 mg·g-1 FW，较CK分别增加40.46%、9.20%和40.63%。抗氧化防御方面，SOD和CAT活性随盐浓度梯度呈先升后降趋势，均在CC3处理下达到最大值203.00 U·g-1 FW和141.83 U·g-1 FW。MDA含量在低盐处理（CC1、CC2）下与CK无显著差异，但在高盐处理（CC4、CC5）下显著升高。SPAD值在高盐处理（CC5）下显著降低。上述结果表明油莎豆通过积累渗透调节物质和增强抗氧化酶活性抵抗一定程度的盐胁迫损伤，其中1.08 M NaCl（CC3）是诱导油莎豆盐胁迫生理响应的最适浓度。
3.2 油莎豆根际土壤酶活性和微生物生物量对盐胁迫梯度的响应
盐胁迫显著影响油莎豆根际土壤酶活性和微生物生物量。S_CAT活性随盐浓度增加持续升高，较CK增加9.56%-16.85%。URE和ALP活性在盐胁迫下较CK显著降低，最大降幅分别达17.27%和31.95%。MBC、MBN和MBP在盐胁迫下均呈显著下降趋势。结果表明盐胁迫对土壤酶活性和微生物生物量具有差异化调控效应，反映了土壤微生物群落应对盐胁迫的结构和功能分化适应策略。
3.3 不同盐胁迫水平下油莎豆根际细菌群落特征
盐胁迫显著改变油莎豆根际细菌群落的多样性和组成。α多样性方面，Shannon和Sobs指数随盐浓度增加呈先降后升趋势，在CC3处理下达到最低。β多样性方面，基于OTU水平的PCoA分析显示处理组沿第一主坐标轴明显分离，PERMANOVA分析证实不同盐处理间微生物群落结构差异显著（R = 0.8494，P = 0.001）。门水平上，变形菌门（Proteobacteria）、放线菌门（Actinobacteria）和厚壁菌门（Firmicutes）为所有处理下的优势菌门。随盐浓度增加，放线菌门相对丰度呈下降趋势，厚壁菌门先增后减，在CC3处理下达到最大值（为CK的4.88倍）。属水平上，Planococcus、Paracoccus和Allorhizobium–Neorhizobium–Pararhizobium–Rhizobium在盐胁迫下均显著增加，在CC3处理下增幅最大，分别为CK的104.70、4.99和72.35倍。Bacillus相对丰度同样在CC3处理下达到峰值（为CK的2.73倍）。LEfSe分析显示，CK组主要富集放线菌门的Arthrobacter；CC1和CC2主要富集变形菌门的Sphingomonas和根瘤菌科类群；CC3主要富集拟杆菌门的黄杆菌科、厚壁菌门的Bacillus和Planococcus以及变形菌门的Paracoccus；CC4以变形菌门的Halomonas为主；CC5以放线菌门的丙酸杆菌目和norank_f_JG30-KF-CM45为特征。
3.4 驱动油莎豆根际细菌群落组装的关键因素
冗余分析（RDA）结果显示，所选土壤和植物因子共同解释了根际微生物群落结构变异的80.38%（RDA1：66.26%，RDA2：14.12%）。方差分解分析（VPA）表明，所有测定变量对微生物群落变异的约束解释总量为63.14%，其中植物生理生化（P）的独立贡献为28.97%，土壤性状（S）的独立贡献为5.03%。进一步将植物性状分为抗氧化系统（AS，包括CAT、POD、SOD和MDA）和渗透调节系统（OS，包括SP、SS和Pro）后分析发现，OS性状的独立贡献最大，达17.75%，其次为AS，土壤性状独立贡献最小（5.03%）。相关性分析显示，经FDR校正后，Planococcus和Paracoccus与植物SS（r = 0.9071，r = 0.8700）、SP（r = 0.8766，r = 0.8384）、Pro（r = 0.8533，r = 0.8543）、SOD（r = 0.6591，r = 0.8526）和CAT（r = 0.8225，r = 0.8638）均呈极显著正相关（q < 0.01）。Paracoccus与S_CAT显著正相关（r = 0.6201），与ALP、MBC、MBN和MBP显著负相关。Sphingomonas和norank_f_Gemmatimonadaceae与ALP、MBC、MBN和MBP显著正相关，与S_CAT极显著负相关。上述结果表明植物渗透调节物质是盐胁迫下根际微生物群落重塑的首要关联因素。
3.5 盐胁迫下根际细菌功能与植物抗逆性及生物量的相关性分析
FAPROTAX功能预测显示，盐胁迫显著改变根际细菌群落碳、氮、硫循环相关功能的预测相对丰度。化学异养和好氧化学异养是所有处理下的优势碳循环功能。硝化、反硝化等氮循环功能在低中盐处理（CC1-CC3）下显著富集，在高盐处理（CC4、CC5）下逐渐降低。PICRUSt2预测结果表明，KO_NIr（如亚硝酸还原酶）、KO_AS（如谷胱甘肽过氧化物酶）、KO_IT（如Na+转运NADH：泛醌氧化还原酶）和KO_OS（如吡咯啉-5-羧酸还原酶）相关功能在CC1-CC3处理下持续较高，在CC3处理下达到最高后下降。Mantel检验和相关性分析显示，Planococcus和Paracoccus与KO_NIr、KO_IT、KO_AS和KO_OS均呈极显著正相关（Mantel相关系数r = 0.33-0.78）。KO_IT、KO_AS和KO_OS与RFW、SP、SS、CAT和Pro经FDR校正后均呈高度显著正相关（r > 0.64，q < 0.01）。PLS-SEM分析显示，盐胁迫对植物渗透调节和抗氧化系统的直接效应分别为0.537和0.775（P < 0.01）；植物渗透调节物质对根际优势潜在耐盐菌的直接效应为0.857；优势潜在耐盐菌对耐盐代谢通路的直接效应为0.969。中介效应检验表明，植物渗透调节系统-优势菌-细菌耐盐功能通路的链式中介效应为0.831（95%置信区间：0.746-0.909），具有统计学显著性。细菌耐盐功能通路对植物抗氧化系统的直接效应为0.399（P < 0.01），优势耐盐菌-细菌耐盐功能-植物抗氧化系统通路的链式中介效应为0.387（95%置信区间：0.164-0.754），同样具有统计学显著性。
四、讨论与结论
讨论部分指出，油莎豆通过积累渗透调节物质、增强抗氧化防御系统以及将生物量分配转向块茎发育来适应盐胁迫环境，其中CC3（1.08 M NaCl）是诱导生理响应的最适浓度。盐胁迫对土壤酶活性和微生物生物量具有差异化影响，S_CAT活性持续升高可能是土壤微生物消除活性氧的适应性调整，而URE和ALP活性降低及微生物生物量减少反映了微生物应对渗透胁迫的能量资源重新分配。盐胁迫下厚壁菌门和变形菌门显著增加，Planococcus、Paracoccus和Bacillus等具有耐盐潜力的功能菌群富集，这些微生物具有渗透物合成、抗氧化防御和离子稳态调节的遗传潜力。植物渗透调节物质对根际微生物群落结构变异的独立贡献最大，Planococcus和Paracoccus与SS、SP和Pro高度正相关，且二者在CC3处理下达到最大相对丰度，与油莎豆渗透调节物质积累峰值完全吻合，提示植物渗透调节与这些类群富集之间存在潜在关联。预测的细菌耐盐代谢通路与土壤碳氮循环功能及植物抗逆性状显著正相关，PLS-SEM进一步揭示了植物渗透调节物质和细菌代谢功能分别在盐胁迫-细菌功能通路和耐盐菌-植物抗氧化系统通路中的统计中介作用。
研究结论：油莎豆具有一定的耐盐胁迫能力。盐胁迫条件下，植物将生物量分配转向块茎发育，伴随渗透调节物质的显著积累和抗氧化防御谱的变化，这可能与耐盐性相关。在众多处理中，1.08 M NaCl引起油莎豆最显著的生理变化。盐度与油莎豆根际细菌群落变化密切相关，根际细菌α多样性随盐度梯度呈U型变化，Planococcus和Paracoccus等潜在耐盐菌属显著富集。植物渗透调节物质如可溶性糖（SS）和脯氨酸（Pro）与潜在耐盐菌群富集显著相关。预测的与盐胁迫相关的细菌代谢通路与土壤养分循环功能和植物抗逆性状存在协变关系。PLS-SEM分析表明，植物渗透调节因子在潜在耐盐菌群结构与盐胁迫之间表现出统计学中介关系，潜在细菌耐盐功能在耐盐根际细菌与植物抗氧化系统之间也表现出统计学中介关系。总体而言，该数据集表明盐胁迫条件下植物生理特性、根际细菌组成和细菌耐盐代谢之间存在协变关系。油莎豆的耐盐性可能涉及潜在的植物-微生物关联，这为改良盐胁迫土壤中的作物产量提供了有价值的参考。
