美国中西部马铃薯生产系统土壤微生物群落与环境梯度的区域分异关系

《Environmental Microbiology Reports》：Regionally Divergent Relationships Between Microbial Communities and Environmental Gradients in Potato Production Systems

　　

将土壤微生物组数据用作土壤健康指标在农业研究中仍具挑战性。微生物组数据复杂，表现出强烈的时空变异及对环境因子的多样化响应。因此，微生物群落的测定结果往往难以解释，且不能提供关于土壤健康状况的可信推断。本研究评估了美国中西部上部三个不同马铃薯种植区土壤微生物群落的分布，并比较了其与土壤化学和管理因素在各区域间的关系。共采集并分析1270份土壤样品，测定微生物特性，包括细菌和真菌丰度、分类学多样性（taxonomic diversity）和功能多样性（functional diversity）及组成。研究人员发现，微生物群落特性在空间尺度上存在显著变异。环境因子与纲（class）水平群落多样性的关联强于更细分类水平。土壤pH、土壤有机质（soil organic matter, SOM）、轮作制度和熏蒸与微生物特性的不同方面相关，但大多数关系具有区域特异性。结果表明，土壤健康的微生物组分包括群落特性的多个方面，制定促进土壤健康的实践需要区域特异性建议。

研究背景与问题：土壤微生物是土壤健康的重要组成部分，土壤健康涵盖物理、化学和生物属性。生物组分因难以建立可反映微生物介导的多维生态系统功能的指标而更具挑战。尽管微生物群落组成是微生物介导土壤过程的核心，但将其作为土壤健康指标或据此进行可信推断并不直接，原因包括群落结构的地理差异强、结构—功能关系复杂。微生物组数据信息量大且在多个分类水平上保留信息，过去多关注单一分类水平；环境因子与不同分类水平群落结构的关系强度尚需比较。农业土壤中微生物丰度、多样性和组成在不同空间尺度上的变异幅度认识仍有限。环境和管理因子可能随地理背景变化，尺度对生态关系推断和管理策略制定尤为重要。美国中西部上部马铃薯生产系统跨三个不同区域，土壤化学、轮作和熏蒸等管理实践在区域间及区域内变异，为研究土壤微生物空间变异及其环境关系提供机会。研究目标是评估该区域马铃薯种植系统土壤微生物群落的空间分布，比较不同分类水平上微生物分布的驱动因子，并比较区域间及区域内微生物丰度和多样性的决定因素；据此提出三个问题并形成假设：田块尺度变异最大、土壤和管理因子在较宽分类水平解释更多变异、环境因子对微生物群落的影响存在区域差异。
研究人员开展了多区域、多季节的田间调查与分子生态分析。共采集1270份土壤样品，来自美国Minnesota和Wisconsin三区域218块商业马铃薯田，分属24个农场（泥炭区40块、壤土区40块、砂土区138块），覆盖2018—2021三个生长季；每块田设2—3个采样点，分别于秋季（11月）和次年夏季（7月）采集0—15 cm土壤，形成预种植期和生长季样品。研究测定土壤化学和管理历史，并分析细菌和真菌丰度、分类学及功能多样性和组成。结论是三个种植区拥有不同微生物群落，群落特性对环境和管理因子的响应多为区域特异，土壤健康促进实践需区域特异性建议。重要意义在于揭示了微生物指标作为土壤健康指标时需考虑空间尺度、分类水平和区域背景，避免“一刀切”管理。
关键技术方法（≤250字）：样本队列来源为美国中西部上部三区域218块商业马铃薯田1270份土样（2018—2021年，24个农场；泥炭、壤土、砂土三区）。采用定量PCR测定细菌16S rRNA和真菌ITS丰度；对细菌16S rRNA V3—V4区和真菌ITS2区进行高通量扩增子测序，经DADA2和Qiime 2处理，分别以SILVA 138和UNITE数据库注释；以逆Simpson指数（Inverse Simpson index）度量α多样性，以Bray—Curtis距离的PCoA和PERMANOVA分析群落组成；用PICRUSt2预测细菌KO和通路及真菌EC功能；用线性混合效应模型分解区域、农场、田块尺度方差，用方差分解分析季节、空间、土壤化学和管理贡献，用线性混合效应回归结合Type III ANOVA评估区域依赖关系。
研究结果：3.1 土壤微生物性质与土壤化学的描述统计。通过区域和季节比较及排序分析，三区域土壤化学性质不同：泥炭区pH最高、SOM最高，壤土区pH最低且P、K较高，砂土区SOM较低；细菌和真菌丰度、分类学及功能多样性在区域和季节间存在差异（region、region×season，p < 0.05）。细菌丰度和纲水平分类学多样性呈泥炭>壤土>砂土，细菌功能多样性呈相反趋势；真菌丰度和功能多样性呈泥炭=壤土>砂土，纲水平分类学多样性在泥炭最低；属水平细菌和真菌多样性区域差异最小。夏季微生物丰度、分类学多样性和功能多样性一般低于秋季，但砂土微生物丰度、泥炭真菌分类学多样性及壤土和砂土微生物功能多样性存在例外。PCoA和PERMANOVA显示区域和季节均结构化细菌和真菌群落；泥炭区细菌分类学和功能组成与壤土、砂土明显分离，真菌分类学组成三区分离，真菌功能组成分离较弱但显著；季节差异一般小于区域差异（如细菌分类学组成r²=0.01 vs. 0.16）。
（1）从区域到农场到田块，土壤微生物丰度、分类学剖面和功能剖面在哪个空间尺度变异最大？通过线性混合效应模型和PERMANOVA进行空间尺度方差分解，结果表明微生物丰度和多样性（细菌和真菌丰度、分类学和功能多样性）以及真菌分类学和功能组成中，田块解释的方差最多，农场次之，区域最少；细菌分类学和功能组成呈不同模式，区域解释最多，田块次之、农场再次；这些趋势在各分类水平上基本一致。
（2）季节、空间位置、土壤化学和管理如何在不同分类尺度上贡献微生物群落变异？通过方差分解分析，土壤化学对纲水平微生物分类学多样性的解释比例高于更细分类水平；单独建模时，土壤化学分别解释细菌和真菌纲水平分类学多样性边际方差的41%和32%，而目水平为5%和6%、科水平为14%和7%、属水平为7%和2%、ASV水平为2%和6%。土壤化学与季节合计解释细菌和真菌纲水平分类学多样性方差的33%和23%；在更细水平，目水平为7%和13%、科水平为15%和12%、属水平为7%和8%、ASV水平为2%和8%。控制季节、化学和轮作后，空间因素单独解释细菌和真菌纲水平分类学多样性剩余方差的19%和13%，目水平29%和17%、科水平19%和27%、属水平18%和26%、ASV水平19%和21%。对分类学组成，土壤化学在各分类水平均解释显著方差；单独建模时，细菌和真菌纲水平分别为43%和15%，目水平43%和12%，科水平44%和14%，属水平42%和13%，ASV水平24%和13%；土壤化学与季节合计解释纲水平30%和14%，目水平30%和13%，科水平31%和13%，属水平30%和12%，ASV水平17%和12%。空间因素独立解释细菌和真菌分类学组成纲水平18%和19%、目水平19%和23%、科水平19%和24%、属水平19%和25%、ASV水平25%和27%。土壤化学通常对细菌属性的解释比例高于真菌；与季节合计解释细菌丰度38%、功能多样性23%、功能组成25%及分类学组成30%—17%，而真菌相应为9%、7%、10%及14%—12%。空间因素通常对真菌属性解释比例高于细菌，但粗分类水平分类学多样性例外。季节单独解释的变异很小，对细菌和真菌丰度贡献均<1%，对真菌分类学和功能多样性及组成略高于细菌。空间、土壤化学和管理之间存在显著重叠，可能与三者相关有关。
（3）土壤和管理因子对生长季土壤微生物丰度和多样性的影响是一致的还是区域特异的？通过线性混合效应回归和Type III ANOVA，许多土壤化学与生长季微生物属性的关系呈区域依赖。土壤pH影响细菌丰度（p=0.005）、真菌丰度（p=0.081）和真菌分类学多样性（p=0.004），对细菌分类学多样性无总体影响（p=0.981），但上述效应均随区域变化（pH×region交互p=0.006、0.005、0.047和0.001）。pH升高与壤土和泥炭区细菌、真菌丰度升高相关，砂土区无此关系；砂土和壤土区pH升高与细菌分类学多样性升高、真菌分类学多样性降低相关，泥炭区细菌和真菌分类学多样性均随pH升高而降低；pH对细菌和真菌功能多样性影响很小。SOM影响细菌丰度（p=0.009）和分类学多样性（p=0.032），对真菌丰度、真菌分类学多样性及细菌和真菌功能多样性影响较小；SOM对细菌丰度和分类学多样性的效应存在区域交互（均p<0.001）：砂土和壤土区SOM与细菌丰度和分类学多样性正相关，泥炭区无此关系。P和K基本不影响微生物丰度和多样性，仅K×region影响真菌功能多样性（p=0.049）。轮作长度影响细菌丰度（p=0.009），且该效应区域间一致（region×rotation length，p=0.151），轮作间隔较长田块细菌丰度较高；轮作长度不影响真菌丰度及细菌和真菌分类学或功能多样性。轮作多样性不影响细菌和真菌丰度或功能多样性，但影响细菌和真菌分类学多样性（p=0.049、0.039）；轮作多样性一般提高真菌分类学多样性，对细菌分类学多样性的影响随区域变化（region×rotation diversity，p=0.033），壤土区提高、砂土和泥炭区不明显。熏蒸主要影响细菌分类学多样性（p=0.001），不影响真菌分类学多样性、细菌和真菌丰度或功能多样性；熏蒸对细菌分类学多样性的影响区域依赖（region×fumigation，p=0.009），预种植熏蒸降低壤土区生长季细菌多样性，砂土区无此效应；泥炭区调查田块未实施熏蒸。
讨论部分总结：空间尺度上，田块是微生物丰度和多样性变异的最重要来源，支持田块尺度变异最大的假设；但“田块”同时包含田内真实生物变异和实验室技术变异，需技术重复和更大田内样本量区分。细菌组成在区域尺度解释较多，与Klasek et al. (2024)多大陆马铃薯田研究中共有ASV仅占14%的发现一致。土壤化学对纲水平分类学多样性解释较多，可能因宽分类群具有共享形态、生理或生活史性状；环境信号在分类学组成上跨分类水平较一致，与假设相反，提示生态性状在窄分类水平也较大程度保留（Philippot et al. 2010；Bay et al. 2020；Lu et al. 2016）。方差分解提示细菌更受确定性过程影响，真菌更受空间/扩散限制等随机过程影响（Chase and Martiny 2018；Naveed et al. 2016；Zhang et al. 2020）。季节对真菌分类学和功能多样性及组成较细菌更重要，但单独解释方差很小；秋季与夏季细菌群落相似，提示秋季样本可能用于预测后续生长季群落。大量变异未被测定变量解释，涉及未测土壤理化生物性质、管理历史、气候和随机过程，也反映DNA提取、PCR偏倚和测序误差；空间、轮作和土壤化学解释存在重叠，区分自相关需分层采样设计（Baker et al. 2009）。区域依赖方面，pH是主导因子但影响随区域变化，可能与pH范围有关；跨区呈单峰关系。pH对丰度的影响在砂土区不同，可能受低SOM共限制。SOM对细菌影响强于真菌，泥炭区可能饱和且SOM与pH负相关。轮作长度提高细菌生物量，轮作多样性提高多样性，壤土区最明显；砂土区轮作长度更重要，泥炭区响应弱。熏蒸主要降低壤土区细菌多样性，SOM可能调节其影响。分类学多样性与功能多样性解耦，提示功能冗余；低SOM砂土功能多样性仍高，可能与功能特化和窄生态位有关；但功能预测受数据库不完整限制，需谨慎解释（Douglas et al. 2019；Ramond et al. 2025）。
研究结论（Conclusions）翻译：总之，三个种植区拥有不同微生物群落，对土壤化学梯度和农场管理实践常以独特方式响应。泥炭区群落以细菌高丰度和广尺度分类学多样性、真菌相应较低丰度和多样性为特征；砂土区是商业马铃薯生产主要区域，细菌和真菌丰度均低，广分类尺度细菌分类学多样性低，壤土区居间。区域间土壤群落对化学性质和管理实践的响应性也不同。砂土区为寡营养、资源受限环境，支持对SOM轻微增加有响应的细菌群落，但其响应受限，对轮作多样性或熏蒸等管理实践响应很小；壤土区对轮作和熏蒸更敏感，可能因较高有机质提供资源基础，使细菌种群对减少扰动或提高植物多样性呈正响应；泥炭土细菌丰度和广尺度分类学多样性在所有样品中均高，因此对管理不响应，原因相反。土壤是复杂系统，反映化学、物理和生物属性，均可能对从数千年到数月尺度的历史因素以独特方式响应；因此，试图定义一套规定的“健康土壤”指标可能困难，健康土壤可能因空间区域、土壤类型或生产系统而异。本研究调查美国中西部上部马铃薯生产田微生物群落跨空间、时间、化学和管理梯度的变异，发现不同微生物组分（真菌与细菌、丰度与多样性）响应不同土壤性质和管理因子，突出即使在同一生产系统定义单一实践促进“健康”微生物群落的困难；多数土壤化学、管理与土壤微生物群落关系具有区域依赖性，强度和方向在三区域间不同，进一步突出需制定区域特异性土壤健康促进实践建议；“一刀切”土壤管理可能无效甚至产生反效果。