来源：Water Regime Determines Dose-Dependent Growth and Biomass Allocation Responses of Soybean Seedlings to Bacillus subtilis and Bacillus amyloliquefaciens.（Plants (Basel)）| 编译：DNA Lab Space
导读
大豆苗期遭遇水分亏缺会限制叶片扩展和根—茎发育，而微生物生物刺激素的作用可能仅在特定水分背景下显现。本研究在受控环境托盘细胞实验中，以4种水分制度×11个菌剂剂量（264个实验单元）的完全随机设计，揭示了大豆幼苗对枯草芽孢杆菌+淀粉芽孢杆菌混合菌悬液的响应具有强烈的水分依赖性。研究发现，严重干旱（W0）下0.05 mL低剂量可使根干重增加116.7%，而充分灌溉条件下各剂量均未稳定超越对照。多变量分析清晰分离出"温暖低扩展"的干旱表型与"凉爽高活力"的充分灌溉表型。
研究背景
大豆（Glycine max (L.) Merr.）苗期建成依赖于基质、根系与扩展中茎叶之间的快速水分连通性。当供水量低于蒸散需求时，膨压下降、气孔调节、胁迫信号传导和碳分配改变会在生物量库完全耗竭之前约束细胞扩展。在大豆中，渐进性水分限制引发热学、结构、代谢和光合特性的协同变化，而田间响应的幅度因基因型和生产环境而异。这一早期阶段在农艺上至关重要，因为叶片部署受限和根系探索不足会放大基质水分可用性的短期波动，并延迟冠层闭合。
水分限制条件下的管理策略日益涉及根际微生物组的调控。植物相关微生物可影响养分 mobilization、根系发育、植物激素平衡、渗透调节、氧化还原稳态和胁迫后恢复。然而，这些效应是有条件的而非普适的——干旱驱动的微生物组重构、微生物对根系生长的维持、微生物群—内皮层对矿物质稳态的控制、宿主偏好以及胁迫特异性促生效应均表明，最终结果取决于植物基因型、微生物身份和环境历史。
研究方法
4.1 植物材料与受控环境条件
大豆（Glycine max (L.) Merr.）品种BRS 1056 IPRO的种子购自巴西隆德里纳Embrapa大豆研究中心。本研究为受控环境托盘细胞实验，非田间试验。幼苗在25 ± 2 °C、相对湿度60 ± 10%、16小时光周期、光合光子通量密度约500 µmol photons m⁻² s⁻¹的条件下生长。每个填充基质的托盘细胞相当于一个小型盆钵，构成一个实验单元。
在水分制度实施前，所有细胞在相同供水条件下管理以标准化建植过程。在共同建植期间，所有幼苗每周接受相同体积（5 mL）的完全Hoagland营养液，以限制处理间的营养失衡。未施用慢生根瘤菌；有意操控的根际因子为混合芽孢杆菌处理。所选制剂包含两种微生物的共培养混合物。
4.2 实验设计、细菌预调理与水分制度
实验采用完全随机的4 × 11因子设计，包含4个水分制度、11个菌剂剂量和每组合6个生物学重复（共264个独立细胞）。菌剂为含有枯草芽孢杆菌+淀粉芽孢杆菌、浓度为5 × 10⁸ CFU mL⁻¹的未稀释悬浮液。来源记录将微生物鉴定至物种水平，但未包含菌株编号；因此，本研究结论支持的是该混合物种制剂的效应，而非菌株特异性表现。单次根际定向施用0、0.05、0.1、0.2、0.4、0.6、0.8、1.0、1.2、1.75或2.0 mL cell⁻¹。这些处理被称为菌剂剂量，因为悬浮液浓度恒定而施用量变化。0 mL对照接受不含细菌菌剂的水。
菌剂在水力限制之前施用，以确保所有处理的初始建植条件一致，并使悬浮液在干旱发展前能够接触根际。这一顺序代表预防性根际预调理，不包括种子处理、干旱建立后的施用或重复施用。基质田间持水量在饱和并自由排水后通过重量法校准。在随后的15天限制期内，W100细胞维持在100%操作田间持水量；W50和W25分别接受恢复W100所需水量的50%和25%；W0不接受额外水分。因此，W50和W25表示补充比例，而非50%和25%体积含水量的固定测量值。W0作为严重端元处理：植物带着基质中储存的水分进入限制期并逐渐消耗之，而非在初始干燥基质中建植。托盘定期重新摆放以减少位置效应。
4.3 生长、生物量分配与热成像测定
限制期结束时，测量同一植株的根长（RL）、茎长（SL）、茎直径和总叶面积（LA）。离体小叶经扫描后用ImageJ 1.54 s版（美国国立卫生研究院，Bethesda, MD, USA）分析。图像使用每次扫描中包含的参考标尺进行空间校准。植株分为根、茎和叶，组织在70 °C鼓风干燥箱中干燥至恒重。根干重（RDM）、茎干重（SDM）、叶干重（LDM）和总干重（TDM）以g plant⁻¹表示。总长（TL）、干重分数和比叶面积（SLA）按照既定的叶片构建和生物量分配定义计算，公式如下：
TDM = RDM + SDM + LDM，TL = RL + SL
RMF = RDM/TDM，SMF = SDM/TDM，LMF = LDM/TDM，SLA = LA/LDM
SLA以cm² g⁻¹表示，RMF、SMF和LMF以g g⁻¹表示。热辐射图像使用FLIR EDGE Pro热像仪（FLIR Systems, Danderyd, Sweden；7.5–13.5 µm光谱波段）从固定距离垂直于冠层采集完整植株图像。叶片发射率设为0.98，遵循先前应用于渐进性水分亏缺大豆的方案。热图像采集和解释遵循既定的植物胁迫表型鉴定建议。
4.4 统计分析
数据以均值±标准误（SEs）报告，每个水分×剂量组合含6个独立生物学重复。水分制度（W100、W50、W25和W0）、菌剂剂量（0、0.05、0.1、0.2、0.4、0.6、0.8、1.0、1.2、1.75和2.0 mL）及W × B作为III型双因素方差分析的固定效应。剂量在主因子分析中作为分类变量处理，仅在图示中按顺序排列。残差正态性和异方差性分别通过Shapiro–Wilk检验和Breusch–Pagan检验考察。III型检验对分类因子使用和为零对比，并进行了异方差稳健检验（HC3）作为敏感性分析：
Var̂ᴴᶜ³β̂ = (X′X)⁻¹X′diag(eᵢ²/(1−hᵢᵢ)²)X(X′X)⁻¹
其中Var̂ᴴᶜ³β̂为估计回归系数的HC3异方差一致性协方差矩阵；β̂为估计模型系数向量；X为模型矩阵；X′为X的转置；(X′X)⁻¹为叉积矩阵的逆；eᵢ为第i个观测的普通最小二乘残差；hᵢᵢ为第i个观测的杠杆值，对应帽子矩阵H = X(X′X)⁻¹X′的第i个对角元素；diag(·)表示对角矩阵，其第i个元素为eᵢ²/(1−hᵢᵢ)²。因子(1−hᵢᵢ)⁻²增加了高杠杆观测残差的贡献，提供异方差稳健标准误、置信区间和Wald检验。
所有配对估计均值比较采用Tukey调整，而各水分制度内每个非零剂量与对应0 mL对照的计划对比采用Holm多重性调整。图中代码ns、、和分别表示不显著（p ≥ 0.05）、p < 0.05、p < 0.01和p < 0.001。模型整体显著性、百分比变化、调整后制度内对比和综合排名见补充图S1–S11。残差诊断、效应估计、置信区间和多重性调整p值被联合解释，而非仅以p值作为推断依据。
综合评分并非额外的假设检验。它将10个测量响应压缩为一个描述性平均z-score，仅用于在各水分制度内对处理组合进行排序；这正是摘要、第2.6节和第5节所报告排名的计算方法。对于处理组合j和性状k，处理细胞均值在44个水分×剂量组合中进行标准化：
zⱼₖ = (x̄ⱼₖ − μₖ)/sₖ
其中分子为处理细胞均值与性状k跨细胞均值之差，分母为该性状44个处理细胞均值的标准差。纳入10个性状：叶面温度、根长、茎长、茎直径、叶面积、总长、根干重、茎干重、叶干重和总干重。脚本定义的综合平均z-score计算如下：
Iⱼ = (1/10) Σₖ₌₁¹⁰ zⱼₖ
所有十个标准化性状权重相等，叶面温度按其观测方向保留，如分析脚本中所实现。因此，较高的分数表示十个标准化响应的算术平均值较高；它本身并非通用的生理效益评分。分数分别在W0、W25、W50和W100内按降序排列。
研究结果
2.1 水分制度主导热学与结构响应
因子方差分析显示，水分制度（W）、菌剂剂量（B）及W × B交互作用显著影响叶面温度、根长、茎长、茎直径、叶面积和总长（图1B–G）。对于这些性状，W和W × B在p < 0.001水平显著；剂量效应在p < 0.001水平显著，总长除外（p < 0.01）。可见光和热辐射图像与模型结果一致：W0幼苗冠层覆盖较少、表面较温暖，而W100幼苗叶片较大、表面较凉爽（图1A）。
跨菌剂剂量平均，W0较W100使叶面温度升高8.8%（22.94对21.09 °C），根长减少44.3%，茎长减少25.4%，茎直径减少38.2%，叶面积减少72.3%，总长减少35.6%（图1B–G）。W50和W25通常处于中间位置，确认处理产生了有序的供水梯度，而非简单的充分灌溉对干旱对比。叶面积是响应最敏感的结构性状，其次是根长和茎直径（图S1–S11）。
显著的交互作用反映了各水分制度内不同的剂量响应曲线（图1B–G）。在W0下，最大平均根长和总长出现在最低剂量，尤其是0.05 mL，而增加剂量并未导致渐进性增长。在W25下，茎长和茎直径在0.2 mL附近最有利，而总长在0.8 mL附近达到峰值。在W50下，若干中高剂量改善了特定性状，包括1.2 mL下的茎长和2.0 mL下的叶面积。在W100下，未处理基线已较高，有利响应具有性状特异性而非在所有指标上一致（图S9）。
多重性调整后的制度内比较进一步说明了这种条件性。与相应未处理对照相比，茎长在W25的0.2 mL下增加28.8%（p = 0.0116），在W50的1.2 mL下增加58.3%（p < 0.001）。叶面积在W50的2.0 mL下增加89.8%（p = 0.0098），在W100的0.2 mL下增加119.4%（p < 0.001）。相比之下，若干视觉上呈阳性的W0响应在多重性调整后未显著变化，表明严重水分限制增加了重复间的离散度，降低了一些剂量对比的分辨率。这些性状的完整因子显著性图谱、百分比变化和多重性调整后剂量对对照估计值见补充图S1、S2、S5和S9–S11。
【数据】叶面温度：W0较W100升高8.8%（22.94 vs 21.09 °C）；根长减少44.3%；茎长减少25.4%；茎直径减少38.2%；叶面积减少72.3%；总长减少35.6%；茎长：W25×0.2 mL增加28.8%（p=0.0116），W50×1.2 mL增加58.3%（p<0.001）；叶面积：W50×2.0 mL增加89.8%（p=0.0098），W100×0.2 mL增加119.4%（p<0.001）
2.2 生物量响应取决于器官、剂量与水分制度
水分制度显著影响根、茎、叶和总干重，菌剂剂量显著影响全部四个生物量性状（图2）。W × B交互作用对根干重、叶干重和总干重显著；茎干重的交互作用不显著（p = 0.0552）。因此，水分可用性设定了整体生物量上限，而细菌响应因器官和水分制度而异。
与W100相比，W0使根干重减少23.2%，茎干重减少24.6%，总干重减少20.1%（图2A,B,D）。叶干重仅减少9.5%，远低于叶面积72.3%的降幅。这一差异表明，严重水分限制对叶片构建和扩展的影响远强于对叶片生物量库的消耗。
W0下最清晰的阳性生物量响应出现在最低非零剂量。与未处理W0对照相比，0.05 mL使根干重增加116.7%（p < 0.001），总干重增加39.3%（p = 0.00368）（图2A,D）。在W50下，1.2 mL使茎干重增加37.6%（p = 0.0121），而最大平均叶干重出现在2.0 mL。对于茎干重，经典III型方差分析的W × B交互作用不显著（p = 0.0552），但HC3敏感性分析显著（p = 0.0147）；因此对该器官特异性交互作用持谨慎解释。在W100下，没有任何单一剂量在全部四个生物量区室中一致超过未处理对照。生物量的完整百分比变化和调整后对比可视化见补充图S3和S6。
总体而言，生物量数据不支持单一的推荐剂量。低剂量与W0下根干重和总干重最相关，而中高剂量在有一定水分可用时对特定器官有利。这些是处理特异性的幼苗响应，但没有证据表明增加菌剂体积必然增加生物量，尽管生物量分配确实受到影响（图2和图S10）。
【数据】根干重：W0较W100减少23.2%；茎干重减少24.6%；总干重减少20.1%；叶干重仅减少9.5%；W0×0.05 mL：根干重增加116.7%（p<0.001），总干重增加39.3%（p=0.00368）；W50×1.2 mL：茎干重增加37.6%（p=0.0121）；茎干重W×B交互：p=0.0552（经典ANOVA）/p=0.0147（HC3）
2.3 水分限制改变生物量分配与叶片构建
水分制度、菌剂剂量及其交互作用显著影响根质量分数（RMF）、茎质量分数（SMF）、叶质量分数（LMF）和比叶面积（SLA）（图3）。W × B项对RMF、LMF和SLA在p < 0.001水平显著，对SMF在p = 0.003水平显著。这些响应表明处理效应涉及比例性生物量分配以及绝对植株大小。
跨剂量平均，W0较W100使SLA降低70.4%（50.19对169.78 cm² g⁻¹），但使LMF增加17.3%（图3C,D）。SLA急剧降低、叶面积大幅减小而叶干重损失相对温和的组合表明，W0幼苗每单位叶生物量产生的表面积更少。RMF和SMF在W0–W100均值对比中变化较小，但两个分数在个别水分制度内仍对菌剂剂量有响应。
各分配性状的剂量排名不同。在W0下，0.05 mL产生最高RMF，而1.0 mL产生最高观测SMF和LMF均值。SLA在W50和W100下均在1.75 mL时最高，分别达到113.0和247.9 cm² g⁻¹（图3）。因此，菌剂可根据水分可用性将根系投资、地上部分配和叶片扩展效率引向不同方向。相应的百分比变化和调整后对比可视化见补充图S4、S7和S11。
【数据】SLA：W0较W100降低70.4%（50.19 vs 169.78 cm² g⁻¹）；LMF：W0较W100增加17.3%；SLA峰值：W50×1.75 mL=113.0 cm² g⁻¹，W100×1.75 mL=247.9 cm² g⁻¹；W×B：RMF/LMF/SLA p<0.001，SMF p=0.003
2.4 热状态与协调生长的模块共变
在264个生物学重复中，叶面温度与大多数扩展性状呈负相关，包括叶面积（r = −0.712）、茎直径（r = −0.703）、SLA（r = −0.689）、总长（r = −0.680）、茎长（r = −0.613）和根长（r = −0.549）（图4）。因此，较温暖的幼苗通常较小，单位叶质量的叶表面积较少。
正相关定义了互补的性状模块。总长与根长（r = 0.916）和茎长（r = 0.710）强相关，叶面积与SLA（r = 0.891）强相关，根干重与RMF（r = 0.897）和总干重（r = 0.823）相关（图4）。这些关系为多变量分析提供了依据，但也表明若干变量代表了全株大小和分配的协调、部分冗余方面。
【数据】叶面温度与叶面积相关系数：r=−0.712；与茎直径：r=−0.703；与SLA：r=−0.689；与总长：r=−0.680；与茎长：r=−0.613；与根长：r=−0.549；总长与根长：r=0.916；总长与茎长：r=0.710；叶面积与SLA：r=0.891；根干重与RMF：r=0.897；根干重与总干重：r=0.823
2.5 层次聚类解析处理水平表型
对44个水分×剂量处理的层次聚类根据标准化多变量谱将组合分离开来（图5）。W100处理产生相对较高的叶面积、SLA、总长、根长、茎长和茎直径，而W0处理产生相对较高的温度（蓝色）和较低的扩展（红色）。W25和W50占据中间区域，尽管个别剂量在树状图中的位置有所移动（图5）。z标准化矩阵包含每个处理×性状组合的一个数值。
【数据】聚类方法：Euclidean距离+完全链接聚合；标准化：z-score跨44个处理均值；矩阵：44×14
2.6 主成分分析
对所有264个独立托盘细胞观测进行主成分分析。14个输入变量为叶面温度、根长、茎长、茎直径、叶面积、总长、根干重、茎干重、叶干重、总干重、RMF、SMF、LMF和SLA。变量在分析前进行中心化和单位方差缩放。PC1和PC2分别解释总方差的39.8%和20.1%。W0处理在PC1负方向形成紧密集群，具有较高温度和较低扩展特征；W100处理在正方向，具有较低温度和较高活力特征。95%多元t分布数据椭圆进一步确认了四个水分制度在多变量空间中的分离。
【数据】PC1解释方差：39.8%；PC2解释方差：20.1%；观测数：264；椭圆：95%多元t分布数据椭圆
讨论与解读
3.1 水分限制约束蒸腾冷却与器官扩展
第一个假设得到支持：水分制度对幼苗表现施加了主导性约束。W0植株相对温暖，在叶面积、根长、茎直径、总长和SLA方面受限最大，而W25和W50通常处于中间位置。这种协调响应与涉及膨压丧失、水力限制、气孔调节、胁迫信号传导和生长—防御权衡的一般干旱生理学一致。根系水力学、结构可塑性和表型框架解释了为何响应同时表现在伸长、叶片部署和生物量分配上。重要的是，叶面积的降幅远大于总干重或叶干重的降幅——生物量的部署方式在生物量库被比例性耗竭之前就已改变。叶面温度作为整合性、非破坏性指标反映了这一水分限制综合征，而非任何单一生理过程的直接测量。
3.2 细菌响应具有剂量依赖性和条件性
第二个假设同样得到支持：细菌响应为非线性且因水分制度而异。最引人注目的效应出现在W0下0.05 mL低剂量，根干重增加116.7%，总干重增加39.3%。然而，在W100下，没有任何剂量能一致超越未处理对照。这一模式与"微生物促生效应取决于环境背景"的文献一致——在资源充足条件下，宿主可能已经接近其生长潜力，额外的微生物投入不再产生正回报；而在胁迫条件下，低剂量微生物可能通过激素调节、渗透保护或营养动员等机制产生显著增益。综合z-score排名显示，最佳剂量随水分可用性变化：W0为0.05 mL，W25为0.2 mL，W50为1.2 mL，W100为0.2 mL。
编译者解读
本研究的设计亮点在于将剂量梯度（11个水平）与水分梯度（4个水平）正交组合，系统揭示了"微生物促生效应依赖环境水分背景"这一核心结论。从合成生物学应用角度看，该工作提醒我们：菌剂开发不能仅追求高剂量或广谱适用，而应针对特定水分生态位优化施用策略。低剂量在极端干旱下的显著促生效应尤其值得关注——这可能意味着微生物信号触发而非营养供给是主要机制。局限在于未提供菌株编号，且为受控环境单时间点测定，田间验证和菌株层面解析将是后续方向。
参考来源
da Silva JPA, Crusiol LGT, Foloni JSS, Farias JRB, Chicati ML. Water Regime Determines Dose-Dependent Growth and Biomass Allocation Responses of Soybean Seedlings to Bacillus subtilis and Bacillus amyloliquefaciens. Plants (Basel). PMID: 42588771.
PMID: 42588771