《Frontiers in Agronomy》：Orchard context dominates first-year soil–plant variation during repeated liquid humus application in commercial citrus orchards
编辑推荐：
：腐殖质类产品在多年生作物系统中应用日益广泛，但其在商业条件下的早期土壤-植株响应可能在很大程度上取决于果园背景（orchard context）。研究评估了西班牙南部两个背景迥异的商业甜橙（sweet orange）果园中重复施用液态腐殖质（liquid h
：腐殖质类产品在多年生作物系统中应用日益广泛，但其在商业条件下的早期土壤-植株响应可能在很大程度上取决于果园背景（orchard context）。研究评估了西班牙南部两个背景迥异的商业甜橙（sweet orange）果园中重复施用液态腐殖质（liquid humus）的首个年度周期。在对照与处理商业田间区域各监测30株树，共进行5次滴灌施肥（fertigation）施用（300 L ha-1）。每个果园内的比较涉及一个对照区与一个处理区，缺乏独立的田间区域重复或处理前基线测量。研究人员采用互补的单变量与多变量方法，评估了土壤理化性质、土壤/根际相关代谢组学（metabolomic）图谱、叶片营养、果实品质、采后失重（postharvest weight loss）及冠层标准化产量（canopy-standardized production）。果园背景主导了系统结构，解释了共享代谢物矩阵（shared-metabolite matrix）中25.4%的变异；全联合矩阵（all-union matrices）中更大的Site组分（50.7–61.4%）还纳入了分析特征表征的差异。土壤、代谢组学与叶片图谱在全局上强烈关联（Spearman's ρ = 0.662–0.758），但大多数关系在单个果园内未检测到。处理解释了全局单个代谢物变异中较小但统计显著的比例（共享矩阵中5.4%，全联合单个代谢物矩阵中2.45%），但这些信号在两个果园内均未重现，因此并不表明跨果园背景存在一致的代谢组学响应。其他对照-处理田间区域差异具有果园特异性：Andarax处理树叶片Fe和B浓度较高，且累积采后失重更大；而COARE处理树叶片Na较高，果实可滴定酸度（titratable acidity）较低。两个预设的冠层标准化产量终点均未检测到统计支持的对照-处理差异。这些结果表明，首年田间区域相关差异具有选择性和果园特异性，在强烈的果园级结构背景下显现，且无可重现的果园内代谢组学响应或统计支持的产量差异。
论文解读文章
一、研究背景与问题
腐殖质类产品（humic-based products）作为土壤改良剂和生物刺激素（biostimulant）在农业中应用日益广泛，但其效应因产品组成、施用策略、作物与环境条件而异。田间与荟萃分析证据显示，植物响应存在显著异质性，包括地下部或营养性状发生变化但生物量或产量未同步增加的情况。针对柑橘类作物中腐殖质类产品的证据仍然有限，田间研究报告了施用后矿质营养与果实品质性状的变化，但响应因制剂和施肥策略而异，这限制了对不同产品间的直接外推，并留下关于商业果园条件下早期响应如何在土壤与植株各组分间分布的疑问。商业柑橘园是长寿命、空间异质性强的系统，新投入品叠加在多个生长季积累的土壤与管理遗产之上。地中海柑橘园的比较研究记录了不同管理与灌溉制度下土壤物理、化学和生物性质的显著变异，近期证据表明环境与农艺背景强烈塑造柑橘土壤微生物群落。因此，对新投入品的早期响应需相对于果园背景并在土壤-植株系统的多个组分间进行解读，而非仅通过短期产量判断。非靶向土壤与根际代谢组学为表征与根系输入、微生物代谢及土壤有机质转化相关的下方化学结构提供了补充手段，但这些图谱高度依赖背景。然而，在商业柑橘园中，将根际代谢组学结构与土壤、叶片、果实、采后及产量响应相整合的田间评估仍然稀缺。本研究评估了重复液态腐殖质施用首年在两个背景迥异的商业柑橘园中的土壤-植株变异，旨在量化果园背景与对照-处理田间区域分化对土壤与植株响应矩阵的相对贡献，识别果园与组分特异的田间区域差异，并评估土壤、代谢组学、叶片、果实、采后与产量组分间的关联。
二、研究内容与主要结论
研究在西班牙南部两个商业甜橙园开展，即Almería的Cítricos del Andarax（有机管理，'Navel Powell'品种，Citrus macrophylla砧木）与Córdoba的COARE（常规管理，'Navel Foyos'品种，Citrus aurantium L.砧木）。两果园在气候、土壤、管理与品种上差异显著。每个果园内设对照与处理两个相邻灌溉扇区，处理区经滴灌系统施用商业液态腐殖质产品HUMUS LIQUIDO（Nostoc Biotechnologies, S.L.），2025年4月至8月每月初施用一次，每次60 L ha-1，首年累计300 L ha-1。Andarax监测10株树（对照5株、处理5株），COARE监测20株树（对照10株、处理10株），两果园监测密度均为每公顷10株。研究无独立田间区域重复，亦无处理前基线测量，因此所有统计分析均不作确证性检验。主要结论为：果园背景主导了首年土壤-植株系统结构；对照-处理田间区域差异具有选择性和果园特异性，涉及叶片营养、果实酸度与采后失重等不同组分；未检测到可重现的果园内代谢组学响应；两个预设冠层标准化产量终点均无统计支持的差异。
三、关键技术方法
研究采用多项互补技术。土壤理化分析由经ENAC认可的LABCOLOR实验室（COEXPHAL）完成，包括pH、电导率、饱和百分比、粒径分布、交换性阳离子、饱和浸提液离子、DTPA提取态金属、有机质、总碳酸盐与速效磷等。土壤/根际相关代谢组学分析由NGAlab（ATENS）完成，采用80%甲醇-0.1%甲酸提取后经UHPLC/Q-TOF（Agilent 1290-G6550）非靶向分析，MS-DIAL 4.9进行特征检测与对齐，多阶段注释流程结合MassBank of North America与自建库、ChemDistiller及Retip。叶片营养分析含N、P、Mg、Ca、K、Na、Zn、Fe、Mn、Cu、B及水分含量。果实品质分析含果汁pH、可滴定酸度、可溶性固形物、硬度、果皮厚度与CIELAB颜色。采后失重评估为25°C环境条件下五周每周称重。产量通过人工计数总果数乘以单果均重估算，冠层体积按三轴椭球近似，冠层标准化产量终点为冠层调整产量（kg m-3）与果数密度（fruits m-3）。统计分析采用单变量ANOVA（Benjamini-Hochberg FDR校正）、PCA、PERMANOVA、PERMDISP、RELATE、DistLM与dbRDA（PRIMER 7 + PERMANOVA+）。
四、研究结果
3.1 果园系统间的土壤分化：全球PCA第一轴解释66.7%变异，样本沿第一轴按果园强烈分离。Andarax主要关联盐度、硫酸盐与钠吸附比相关变量，COARE关联交换性阳离子与金属相关变量（包括土壤Cu）。两果园在质地、盐度/钠吸附比、碳酸盐含量、交换性Mg及DTPA提取态Cu和Zn上差异显著。单变量FDR校正后无任何土壤变量在任一果园内显示对照-处理差异。Andarax质地块检测到探索性多变量对比（解释31.1%变异，permutation p = 0.033），但因土壤粒径分布不可能因五次液态腐殖质施用而在单个年度周期内发生实质性变化，且PERMDISP显示多变量离散度存在差异（p = 0.007），该结果被解释为既有空间异质性和/或离散度差异而非腐殖质响应。
3.2 果园特异的叶片、果实、采后与产量对比：Andarax处理树叶片Fe和B浓度显著高于对照（均q = 0.027），组均值分别高约43%和41%；完整叶片营养图谱亦显示对照-处理差异（解释23.4%变异，permutation p = 0.032）。COARE处理树叶片Na较高（q = 0.048，组均值高约58%），但完整叶片图谱无显著差异（permutation p = 0.161），尽管PERMDISP显示离散度差异（p = 0.002）。果实品质全局分析显示Site显著（p = 0.001）与Site × Treatment交互显著（p = 0.012），Treatment不显著（p = 0.543），但Site与交互项因离散度差异需谨慎解读；果园内对比在Andarax（p = 0.052）与COARE（p = 0.061）均未达显著。内部果实品质中，COARE处理果可滴定酸度显著较低（q = 0.015，组均值低约27%），内部品质图谱显示对照-处理差异（permutation p = 0.022）；Andarax对应对比不显著（p = 0.071）。物理品质与颜色图谱在任一果园内均无显著对照-处理分离。采后失重方面，Andarax处理果周失重图谱显著不同（解释49.2%变异，permutation p = 0.017），失重AUC显著较高（q = 0.043，组均值高约15%）；COARE无显著对比（p = 0.576）。两个预设冠层标准化产量终点在任一果园内均无统计支持的对照-处理差异。
3.3 土壤/根际相关代谢组学图谱：Andarax数据集含386个特征，COARE含205个特征。两果园内在单个代谢物、NPClassifier超类与类水平上均无显著对照-处理多变量分离。共享代谢物矩阵（75个代谢物）中，Site解释25.4%变异（pseudo-F = 9.90，permutation p = 0.001），Treatment解释5.4%（pseudo-F = 2.10，permutation p = 0.037），Site × Treatment不显著（p = 0.300）；但该全局Treatment信号未伴随任一果园内的显著对照-处理分离。全联合矩阵中Site解释50.7–61.4%变异，但因包含结构零值（structural zeros），Site组分同时反映生物分化与分析特征表征差异。单个代谢物水平全联合Treatment项解释2.45%变异（permutation p = 0.021），Site × Treatment解释2.10%（permutation p = 0.047），但后者因离散度差异需谨慎解读；超类与类水平Treatment均不显著。单个特征筛查中Andarax 6个、COARE 13个特征显示名义p < 0.05，但FDR校正后无一保留。
3.4 首年评估中的跨块关联：全局RELATE分析显示土壤与代谢组学图谱强关联（中位ρ = 0.758，permutation p = 0.001），土壤与叶片营养关联（ρ = 0.708，p = 0.001），代谢组学与叶片营养关联（中位ρ = 0.662，p = 0.001）。盐度/钠吸附比组分对代谢组学结构提供最强单独解释（54.3–51.0%变异）。但这些全局关联在单个果园内未检测到。内部果实品质与叶片营养（ρ = 0.336）、土壤（ρ = 0.274）及代谢组学（中位ρ = 0.188）呈较弱全局关联。产量与任何其他响应块无显著全局关联；COARE内土壤与产量关联（ρ = 0.317，p = 0.010），金属相关土壤组分解释34.9%产量变异。果实颜色与土壤（ρ = 0.715）、叶片营养（ρ = 0.735）及代谢组学（中位ρ = 0.713）呈强全局关联，但果园内大多不显著。
五、讨论与结论
讨论部分指出，果园背景是首年土壤-植株结构的主要来源，Site代表整合的商业果园背景，包含地理、土壤、管理历史、品种、环境条件与采样日期等不可分割的差异，不能归因于任何单一因素。共享代谢物矩阵提供更保守的跨果园生物学比较基础，而全联合矩阵的Site组分必然结合生物分化与数据集结构效应。全局跨块关联主要代表果园间协调分化而非果园内树间耦合，且因土壤采样先于叶片与果实采样31–33天，应视为首年采样窗口内的整合而非严格同步耦合。Andarax叶片Fe和B差异的幅度与营养背景相关：按已发表柑橘叶片分析指南，Fe浓度约从缺乏类（<35 mg kg-1）移至偏低类（35–59 mg kg-1），B从适宜类上限移至偏高类（101–200 mg kg-1）但未达过量，这些仅作背景基准而非正式营养诊断。COARE叶片Na差异虽统计支持，但绝对值仍远低于高Na参考范围，不应解读为营养益处或Na毒害证据；可滴定酸度降低被解读为成分变化而非明确品质改善。研究结论为：首年田间区域相关差异具有选择性和果园特异性，在强烈果园级结构背景下显现，无可重现的果园内代谢组学响应，亦无统计支持的产量差异。液态腐殖质重复施用的持续多年度项目将继续对同一批标识树进行后续年度周期评估。
