有益微生物可以促进植物生长，而有害微生物可能抑制植物生长。过去很多研究聚焦于单一微生物与宿主植物的互作，但植物如何响应整个土壤微生物群落，尤其是在不同磷供应条件下如何改变生长和根系转录组，目前仍了解不足。
这项研究关注一个核心问题：土壤微生物对玉米生长的影响是否取决于土壤磷水平。
作者利用一个长期定位试验中的土壤，选择了两种不同磷肥管理背景：
然后将玉米种植在两类土壤中，并设置：
作者主要分析了：
- 土壤微生物对玉米生长的作用取决于磷输入水平。土壤微生物并不总是促进植物生长。在不同长期磷肥背景下，它们可能促进或抑制玉米生长。
- 土壤微生物调控玉米根系磷饥饿响应基因。微生物可以促进或抑制不同的 PSR genes，说明根系对磷缺乏的响应受到土壤微生物群落调控。
- 土壤微生物诱导防御相关代谢过程。无论外源磷水平如何，土壤微生物都能诱导一些防御相关代谢过程。这说明玉米根系在感知微生物群落时，会同时启动营养响应和防御响应。
- 高磷长期施肥背景下，微生物降低 PSR 和防御响应强度，并抑制玉米生长。在长期高磷输入土壤中，微生物群落与磷肥历史共同作用，使玉米根系的磷饥饿响应和防御响应较弱，最终表现为生长受抑。
- 低磷背景下，微生物更强地诱导 PSR 和防御响应，并促进磷吸收和玉米生长。在低磷土壤中，土壤微生物诱导更强的磷饥饿响应和防御相关响应，使玉米能够更有效吸收磷，并促进生长。
这项研究表明，土壤微生物对玉米生长的影响具有明显的土壤磷依赖性。在低磷条件下，微生物可以帮助玉米增强磷饥饿响应、提高磷吸收并促进生长；但在长期高磷输入条件下，微生物可能改变根系转录响应并抑制玉米生长。
这篇文章的重要意义在于，它把植物营养响应和微生物诱导防御响应联系起来，说明植物对土壤微生物的响应不是单纯的“促生”或“抑制”，而是取决于养分背景。
1. 试验设计和土壤来源
盆栽实验在中国农业大学温室中进行。土壤来自中国农业大学上庄试验站的长期定位试验。该长期试验始于 2009 年，设置两个磷肥输入水平：
| |
|---|
| |
| 高磷，150 kg P ha⁻¹ yr⁻¹，以过磷酸钙施入 |
作者取 0–30 cm 土层土壤，并在盆栽前补充不同水平的 KH₂PO₄：
同时，所有土壤都补充相同水平的氮和钾，均为 200 mg kg⁻¹，分别使用 Ca(NO₃)₂·4H₂O 和 K₂SO₄。
2. 微生物处理
作者将 LP 和 HP 土壤各分成两部分：
灭菌处理使用 25 kGy γ 射线辐照 3 天。
最终形成 4 个处理：
每个处理有 7 个生物学重复，每个重复是一盆玉米，每盆一株玉米幼苗。玉米材料为自交系 B73。
3. 玉米种植和采样
玉米种子先经过处理：
每盆装 6 kg 干土。播种后用蒸馏水浇灌，并用无菌塑料袋覆盖，减少外源灰尘污染。
培养 40 天，即六叶期后采样。
采样内容包括：
- 根系，用于 AMF 侵染率、根长、根密度、转录组，以及根内细菌、真菌和 AMF 群落分析；
每组中有 3 株玉米用于根际微生物和根系相关分析。
4. 生物量、根系性状和离子组测定
作者将玉米根系分为：
这些根系使用 Epson 4880 扫描仪扫描，并用 WinRhizo Pro Vision 5.0a 分析根系形态。
测定内容包括：
地上部和根系样品在 65°C 烘干后称重，并研磨过 1 mm 筛。约 0.3 g 样品用浓硝酸和 H₂O₂ 微波消解，然后使用 ICP-OES 测定离子组。
离子组用于分析土壤微生物和磷水平是否改变玉米对矿质元素的吸收。
5. 细菌、真菌和 AMF 群落分析
作者从根系和土壤中提取总 DNA，使用 FastDNA SPIN Kit。
不同微生物类群使用不同扩增体系：
| | |
|---|
| | |
| | |
| | 第一轮 GeoA2/AML2，第二轮 NS31/AMDGR |
测序平台为 Illumina MiSeq，由北京诺禾致源完成。
生物信息分析流程包括：
- FLASH 合并 paired-end reads；
此外，作者还使用：
- FAPROTAX
- FUNGuild
- Bray–Curtis 距离和 PCoA 分析群落结构。
6. AMF 根系侵染率测定
作者选取代表性侧根，用台盼蓝染色后在显微镜下观察。
测定指标主要包括：
该指标用于判断长期磷输入和土壤微生物是否改变玉米与丛枝菌根真菌的共生关系。
7. 玉米根系转录组分析
作者提取玉米根系 RNA，并构建 RNA-seq 文库。
主要流程包括：
- 使用 RNeasy Plus Universal Mini Kit 提取 RNA；
- 使用 Agilent 2100 Bioanalyzer 检测 RNA 质量；
- 使用 NEBNext Ultra RNA Library Prep Kit 构建文库；
- Illumina HiSeq 平台测序，产生 150 bp paired-end reads。
每个实验处理设置 3 个生物学重复用于转录组分析。
8. 转录组数据处理和差异基因分析
RNA-seq 原始数据先进行质控，去除：
clean reads 比对到玉米 B73 参考基因组 RefGen_v3，使用 TopHat v2.0.12。随后用 HTSeq v0.6.1 统计基因 reads，并计算 FPKM。
差异表达基因使用 DESeq2 分析，筛选标准为：
作者重点比较的是：自然微生物土 M 与灭菌土 S 在不同磷水平 LP 和 HP 下引起的根系基因表达差异。
9. 功能注释、代谢通路和互作网络分析
作者使用 MapMan 将差异表达基因映射到代谢通路中，用于分析哪些生物学过程受微生物和磷水平调控。
同时，作者使用 STRING v12 构建显著富集基因的蛋白互作网络，连接置信度阈值为 > 0.4。
此外，作者还使用 Mantel test 分析差异基因表达与离子组之间的关系，使用 R 中的 linkET 包完成。
10. 这篇文章的方法逻辑
这篇文章的方法链条可以概括为：
长期低磷和高磷定位土壤→ 设置自然微生物土和灭菌土→ 种植玉米 B73→ 比较玉米生物量、根系结构和离子组→ 分析根际、根内细菌、真菌和 AMF 群落→ 测定 AMF 侵染率→ 进行根系 RNA-seq→ 解析微生物如何在不同磷背景下调控磷饥饿响应、免疫防御和元素吸收相关基因→ 解释低磷和高磷条件下微生物对玉米生长的相反作用
