东北黑土微生物群落对长期施肥及作物的响应

丁建莉，姜 昕，关大伟，马鸣超，赵百锁，周宝库，曹凤明，李 力，李 俊



东北黑土微生物群落对长期施肥及作物的响应

丁建莉1,2，姜 昕1,2，关大伟1,2，马鸣超1,2，赵百锁2，周宝库3，曹凤明1,2，李 力1,2，李 俊1,2

（1中国农业科学院农业资源与农业区划研究所，北京 100081；2农业部微生物产品质量安全风险评估实验室，北京 100081；3黑龙江省农业科学院土壤肥料与环境资源研究所，哈尔滨 150086）

【目的】为表征长达35年轮作与施肥条件下东北黑土微生物群落特征，解析长期施肥及作物对土壤微生物丰度和群落结构的影响，探讨东北黑土微生物群落变化与施肥及不同作物间的相互关系，为进一步改良耕作制度与施肥方式提供依据。【方法】依托黑龙江省农业科学院长期定位试验站，选取玉米和大豆两种作物季的4个不同施肥处理：不施肥处理（CK）、有机肥处理（M）、无机肥处理（NPK）和有机肥配施无机肥处理（MNPK）的耕层土壤为研究对象，处理前加字母m表示玉米季样品，加字母s表示大豆季样品。借助Illumina Miseq高通量测序平台和Real-time PCR技术，以16S rRNA基因为分子标靶，研究作物与施肥方式对黑土中微生物群落结构和丰度的影响，分析群落变化与土壤化学性质的相关性。【结果】玉米季土壤16S rRNA基因拷贝数（6.32×108—8.83×108/ng DNA）比大豆季的低（0.96×109—2.30×109/ng DNA）；玉米季土壤微生物多样性（ACE指数为3 674.58—4 034.84）也低于大豆季（ACE指数为4 167.47—4 887.36)；玉米季的细菌丰度以Acidobacteria（24.47%—27.90%）最高，而大豆季丰度最高的是Proteobacteria（27.78%—34.40%），Bacteroidetes和Actinobacteria丰度在两季作物中差异明显。同一作物季的有机肥无机肥配施处理的16SrRNA基因拷贝数高于无机肥处理，且有机肥配施无机肥处理的微生物α多样性指数也比无机肥处理的高（sMNPK的Chao1指数比sNPK高出11.89%）；不同施肥处理之间群落组成存在差异，Alphaproteobacteria在sMNPK和sNPK处理的相对丰度分别比sCK增加3.31%、5.24%；Gammaproteobacteria在sMNPK和sNPK处理均比sCK处理增加1.72%和1.2%，二者相对丰度变化大。相关性分析显示，16SrRNA基因拷贝数与土壤硝态氮和速效钾正相关；微生物菌落多样性指数与土壤全氮、硝态氮、铵态氮、有效磷和速效钾等化学性质密切相关。【结论】东北黑土不同作物季和不同施肥均会影响土壤微生物丰度、α多样性和群落结构。有机肥无机肥配施能够有效改变微生物群落结构，提高微生物的丰度和多样性，并提高土壤pH，减缓土壤酸化。

施肥；作物；微生物群落结构；土壤化学性质；高通量测序；黑土

0 引言

【研究意义】中国东北黑土带是世界著名的三大黑土带之一，土壤有机质含量高、养分丰富，是重要的商品粮基地，对国家粮食供应和安全起着举足轻重的作用，但长期过量不合理的施肥和耕作，导致了黑土农田质量日益退化[1]。土壤微生物作为农田生态系统的重要组成部分，能够对长期肥料累积作用产生敏感而快速的响应，是评价土壤质量、土壤肥力和作物生产力的重要指标[2]。因此，开展对黑土中微生物群落结构和丰度的研究，能为进一步揭示土壤微生物与生态环境间的响应机制，以及改进施肥和耕作制度提供依据。【前人研究进展】研究发现，施肥会改变土壤微生物学特性，长期施用化肥能显著提高土壤微生物生物量，并引起群落结构变化[3]，长期施氮肥显著影响固氮菌、氨氧化菌等氮循环菌群[4]，磷肥的长期施用能控制一些特异性细菌及真菌的密度，但对褐灰钙土土壤微生物群落总体的丰富度没有明显影响[5]。施用不同有机肥和氮肥均会引起土壤微生物组成的差异[6]，有机肥能够显著改变微生物群落结构及多样性[7]，提高微生物生物量和代谢活性[8]、以及土壤酶活性[9]。粪肥的长期施用会使作物根系分泌物和土壤中可溶性有机碳含量增加，提高土壤中微生物群落多样性[10]，显著降低土壤真菌与细菌的比率，改变细菌和真菌群落结构[11]，长期施用绿肥促进了某些细菌的生长，从而降低了作物枯萎病的发病率[12]。土壤微生物不仅受施肥方式的影响，也受作物类型的影响[13]，蔡艳等[14]报道不施肥条件下粮豆轮作有利于提高细菌多样性。此外，目前的高通量测序技术可对环境中所有微生物群落进行研究，获得数据量大，并更真实地揭示微生物群落的复杂性和多样性，极大促进了环境中不可培养微生物以及痕量菌的更深入研究[15-16]。Illumina Miseq高通量测序可实现多样品的可变区同时测序，前人应用此方法发现东北黑土长期施用氮肥导致土壤微生物多样性降低，群落结构改变[17]；秦杰等[18]通过此技术初步分析了东北黑土长期施用不同无机肥条件下，土壤微生物群落结构特征和组成差异变化，并取得一定进展。【本研究切入点】由于多年连续施用大量无机肥，尤其是氮肥，造成东北黑土土壤pH持续下降，土壤有机质含量下降，土壤中微生物已由“细菌型”向“真菌型”转变[19]。已有的报道显示有机肥能够改善土壤质量、改变土壤微生物组成[20-21]，但其如何影响土壤质量及响应机制并不清楚。同时除常规施肥等农业措施外，不同作物及种植制度也会对土壤微生物群落结构产生影响，但目前这方面报道较少。【拟解决的关键问题】本研究以黑龙江省农业科学院35年定位试验站的连续两种作物种植季的4个不同施肥处理耕层土壤为研究对象，以16S rRNA基因为分子标靶，借助Illumina MiSeq高通量测序平台和定量PCR技术，探讨长期施肥和作物对东北黑土微生物群落结构和丰度的影响，为农田优化施肥和耕作制度，提高土壤可持续生产力提供依据。

1 材料与方法

1.1 试验地概况

试验在黑龙江省农业科学院长期定位试验站进行，该试验站位于黑龙江省哈尔滨市（45°40’ N，126°35’ E），属于松花江二级阶地，地势平坦，海拔151 m，年均气温3.5℃，无霜期135 d，年降水量533 mm。自1980年起，试验按小麦-大豆-玉米顺序轮作，到2014年和2015年分别为大豆茬（第35个生长季，简称大豆季）和玉米茬（第36个生长季，简称玉米季）。本试验设计4个处理：不施肥处理（CK）；有机肥处理（M）；施氮肥、磷肥和钾肥的无机肥处理（NPK）；有机肥配施无机肥处理（MNPK），3次重复。处理名称前加字母m表示玉米季的样品，如，mCK表示玉米季的不施肥处理；加字母s表示大豆季的样品，如sCK表示大豆季的不施肥处理。氮、磷、钾肥分别为尿素、重过磷酸钙、硫酸钾。氮磷钾肥大豆期施用量为75 kg Nžhm-2，150 kg P2O5žhm-2，75 kg K2Ožhm-2；小麦期和玉米期施用量为150 kg Nžhm-2，75 kg P2O5žhm-2，75 kg K2O žhm-2。有机肥为马粪，施用量为18 600 kgžhm-2。

1.2 样品采集

分别于2014年9月和2015年9月，取作物收获季5—25 cm的耕层土壤，每个重复小区随机取10个点，剔除杂质后混匀作为该处理一个重复的土壤样品。土壤保存于封口塑料袋中，一部分温室风干研磨并过0.2 mm筛，用于测定土壤理化性质；另一部分样品保存于-80℃冰箱，待微生物群落分析。

1.3 试验方法

1.3.1 土壤化学性质测定 铵态氮测定采用KCl浸提-靛酚蓝比色法[22]。参照鲁如坤[23]方法测定以下土壤化学性质。其中，土壤pH测定采用酸度计法（土﹕水=1﹕1）；土壤有机质测定采用浓硫酸-重铬酸钾消煮-硫酸亚铁滴定法；土壤全氮的测定采用半微量凯氏定氮法；土壤有效磷的测定采用NaHCO3浸提-钼锑抗比色法；土壤速效钾的测定采用乙酸铵浸提-原子吸收火焰光度法。

1.3.2 土壤DNA的提取和16S rRNA基因的高通量测序 土壤总DNA使用美国MOBIO公司的Power Max Soil DNA Isolation Kit试剂盒，每个样品称取0.25 g土壤，依照试剂盒说明提取土壤总DNA。所提取的土壤总DNA 的浓度和纯度用核酸定量仪（NanoDrop ND-1000）检测。DNA送至北京奥维森基因科技有限公司，应用Illumina MiSeq 平台对细菌和古菌16S rRNA基因的V4（515F-806R）进行测序[24]。测序原始序列上传至NCBI的SRA数据库（SRP062512）。

1.3.3 16S rRNA基因丰度测定 采用SYBR Green定量PCR法测定16S rRNA基因，反应在ABI 7500 Real-time PCR（ABI，USA）仪器上进行。反应体系为FastFire qPCR PreMix （TIANGEN，China）10 μL，ROX Reference Dye 0.4 μL，1μL DNA 模板，10 n mol·L-1引物（515F-806R），补加ddH2O至20 μL。以含有16S rRNA基因的重组 pGEMR-T 载体为标准质粒，参照文献[25]制备质粒和进行后续工作。

1.4 数据分析处理

通过Illumina MiSeq平台进行Paired-end测序，下机数据经过QIIME（v1.8.0）软件过滤、拼接、去除嵌合体[17]，再聚类为用于物种分类的OUT（Operational Taxonomic Units），OTU相似性设置为97%。对比Greengenes数据库，得到每个OTU 对应的物种分类信息。利用Mothur软件（version 1.31.2）进行α多样分析（包括Shannon、ACE和Chao1等3个指数）。基于Weighted Unifrace距离[26]，使用R（v3.1.1）软件包的pheatmap进行聚类分析[27]。采用SPSS 19.1软件，单因素ANOVA分析和Pearson相关性分析，用Turkey 显著差异法分析处理间的差异显著性；同时应用两因素方差分析（Two-way ANOVA）说明作物效应和施肥效应的显著性。

2 结果

2.1 不同作物和施肥处理对土壤化学性质的影响

两种作物季的不同施肥处理对土壤化学性质影响见表1，施肥和作物对土壤化学性质影响见表2。单因素分析结果表明，与不施肥相比，单施无机肥降低了土壤pH，施有机肥提高了pH，大豆季的sNPK处理的土壤pH比sCK降低0.95，有机肥配施无机肥处理（sMNPK）的pH要高出单施无机肥处理（sNPK）0.4。有机肥配施无机肥和单施无机肥处理极显著增加土壤中全氮、硝态氮、有效磷和有机质等的含量（＜0.05），其中玉米季的mMNPK处理的有效磷含量较mCK处理提高60.81 mg·kg-1，mNPK处理也较mCK处理提高57.81mg·kg-1。总体上，施肥对不同处理的铵态氮含量影响不显著，但两个作物季间的铵态氮含量差异显著；玉米季的铵态氮含量明显低于大豆季的，sNPK处理的铵态氮含量为35.80 mg·kg-1，而mNPK的仅为11.26 mg·kg-1，说明作物季的变化是引起土壤化学性质变化的重要因素之一。同时，通过施肥和作物对土壤化学性质影响的方差分析，可知施肥和作物均显著影响土壤化学性质。

表1 两种作物季不同施肥处理土壤化学性质

表中数据为平均值±标准差；同列数值后不同字母表示处理间差异达显著水平(<0.05)。下同

Data are`x ±; Values followed by different letters in a row are significant among treatments at the 0.05 level. The same as below

表2 施肥和作物对土壤化学性质与微生物多样性指数影响的方差分析

2.2 不同作物和施肥处理对黑土中16S rRNA基因丰度的影响

供试土壤中的16S rRNA基因丰度测定结果如图1所示，玉米季土壤中16S rRNA基因拷贝数为6.32× 108—8.83×108/ng DNA，大豆季拷贝数为0.96×109—2.30×109/ng DNA，表明2种作物季存在显著差异。在玉米季中，mMNPK处理的拷贝数最高，为8.83×108/ng DNA，比mCK和mNPK处理分别提高了39.72%、36.63%；另外，mM处理也较mCK和mNPK处理分别提高23.73%、20.86%。在2种作物季中均为有机肥无机肥配施的拷贝数最高，无机肥处理和不施肥处理间没有显著性差异。16S rRNA基因拷贝数与土壤化学性质的Pearson相关性分析（表3），结果表明，基因拷贝数与土壤硝态氮和速效钾极显著正相关，而与铵态氮极显著负相关。

小写字母表示玉米季不同处理间差异达显著水平；大写字母表示大豆季不同处理间差异达显著水平（P<0.05）

表3 16S rRNA基因拷贝数、α多样性指数与土壤化学性质Pearson相关性分析

**，0.01水平具有显著相关性；*，0.05水平具有显著相关性 **, correlation is significant at the 0.01 level; *, correlation is significant at the 0.05 level

表4 两种作物季不同施肥土壤中微生物α多样性指数分析

2.3 不同作物和施肥处理对黑土中微生物α多样性的影响

两种作物季不同施肥处理的微生物 α多样性指数均值统计情况见表4。无论是多样性指数Shannon，还是丰富度指数ACE与Chao1均是大豆季高于玉米季，两种作物季具有显著性差异。ACE指数玉米季为3 674.58—4 034.84，而大豆季在4 167.47—4 887.36。通过分析可见作物明显影响土壤微生物α多样性。

大豆季4个处理α多样性指数比较，单施有机肥处理（sM）均最高；有机肥配施无机肥（sMNPK）的ACE指数和Chao1指数高于无机肥处理（sNPK），更接近于不施肥处理（sCK）。可见长期施肥会引起土壤微生物α多样性的改变。总体上，有机肥配施无机肥处理α多样性比无机肥处理高。

通过施肥和作物对土壤微生物多样性指数影响的方差分析（表2），可知施肥和作物对多样性指数均有显著影响。进一步分析α多样性指数与土壤化学性质Pearson相关性（表3），从表中数据可知α多样性指数与土壤全氮、硝态氮、铵态氮、有效磷和速效钾等具有相关性。Shannon指数与全氮含量具有极显著负相关性，与铵态氮和速效钾含量具有极显著正相关性。

2.4 不同作物和施肥处理对黑土中微生物组成的影响

经高通量测序分析可知，玉米季有88.53%— 91.37%为细菌，古菌比例为7.53%—11.62%，另外的1.09%—1.17%为不可确定；大豆季有93.91%—96.71%为细菌，古菌仅占2.70%—4.10%，有0.32%—2.18%不能确定。且两种作物季的古菌中均约99%以上为Crenarchaeota，细菌群落对不同作物季的响应更为显著。两种作物季的不同施肥土壤中微生物门水平组成见图2，玉米季细菌主要的5个菌门分别为Acidobacteria（24.47%—27.90%）、Proteobacteria（18.60%—21.43%）、Verrucomicrobia（12.23%—16.93%）、Bacteroidetes（7.16%—10.49%）和Gemmatimonadetes （4.07%—6.14%）。大豆季细菌主要的5个菌门分别为Proteobacteria（27.78%—34.40%）、Acidobacteria（16.39%—20.94%）、Gemmatimonadetes（8.85%—11.03%）、Actinobacteria（8.03%—10.34%）和Verrucomicrobia（7.80%—8.37%）。两个作物季的Bacteroidetes和Actinobacteria 相对丰度差异明显。

表5中列出了大豆玉米2种作物季不同施肥处理土壤中微生物纲水平的组成。从表中可知，Acidobacteria在两种作物季中出现较大变化；Chloracidobacteria在玉米季相对丰度为12.01%—18.55%，而大豆季为4.25%—9.37%，玉米季的比大豆季所占的比例提高了7.76%—9.18%。另外，Chloracidobacteria在同一作物季的施肥处理相对丰度要低于不施肥处理，如mM（16.02%）、mMNPK（13.48%）和mNPK（12.01%）的相对丰度均小于mCK（18.55%）的相对丰度。通过微生物组成与土壤化学性质Pearson相关性分析（表6），可知Acidobacteria相对丰度与土壤pH极显著正相关；与铵态氮、有效磷和速效钾含量极显著负相关。

供试土壤中的Proteobacteria在大豆玉米两种作物季中变化也较大（表5），其中，Betaproteobacteria和Alphaproteobacteria的相对丰度在玉米季分别为5.04%—6.38%和3.74%—5.21%；而在大豆季两者分别为9.94%—11.79%和11.05%—16.9%。Alphaproteobacteria的相对丰度在玉米季各施肥处理没有显著性差异，而在大豆季sMNPK（14.97%）和sNPK（16.90%）之间存在显著差异，sCK（11.66%）和sM（11.05%）之间没有显著性差异。Gammaproteobacteria的相对丰度在sMNPK和sNPK处理比不施肥处理增加1.72%和1.2%。另外，Proteobacteria相对丰度与土壤化学性质Pearson相关性分析（表6），Proteobacteria相对丰度与pH极显著负相关；与硝态氮、铵态氮、有效磷和速效钾含量极显著正相关。

图2 两种作物季不同施肥土样中微生物门水平的组成（相对丰度＞1%）

表5 两种作物季不同施肥土壤中微生物纲水平的组成（8组中至少有1组的相对丰度＞1%）

表6 微生物组成与土壤化学性质Pearson相关性分析

从表5中数据显示，Actinobacteria门中的Actinobacteria纲在玉米季的相对丰度为0.81%—1.16%，大豆季的提高到3.61%—3.54%，相同作物季的不同施肥处理间没有显著性差异。Actinobacteria相对丰度与全氮显著负相关；与铵态氮和速效钾含量极显著正相关。Bacteroidetes门中的Saprospirae的相对丰度，在玉米季3个施肥处理间没有差异，但与不施肥处理存在显著差异；而在大豆季的有机肥处理和有机无机配施处理要显著高于不施肥与无机肥处理。Bacteroidetes相对丰度与铵态氮和速效钾含量极显著负相关。

2.5 不同作物季和施肥处理的聚类分析

通过聚类分析，具有相似β多样性的供试样品聚类在一起，本研究各处理的结果（图3）显示主要聚类成两大分支，每一分支由同一作物季的各个施肥处理聚类组成。玉米季的mNPK处理的几个平行先行相聚，再依次与mMNPK、mM和mCK聚类；大豆季则是sCK与sM先聚类，而后依次与sMNPK和sNPK聚类。

图3 两种作物季不同施肥处理的β多样性分析

3 讨论

3.1 不同作物和施肥处理影响黑土中16S rRNA基因丰度的变化

本研究表明，大豆季黑土中16S rRNA基因丰度明显高于玉米季，大豆季的有机肥配施无机肥处理16S rRNA基因的拷贝数比玉米季的提高了1.60倍。因此可知，种植不同作物是影响土壤16S rRNA基因丰度的一个重要因素，种植大豆提高了土壤16S rRNA基因丰度。另外，Han等[28]利用平板计数的研究也发现豆科植物能够明显增加可培养微生物数量。其原因可能是因为种植作物根系分泌物的不同，而引起根系周围及其土壤中细菌等微生物种群结构的变化，大豆是豆科作物，其根瘤内部具有固氮能力的根瘤菌，能为土壤微生物提供更多的氮源；同时大豆根系分泌物中脂肪酸较多，能活化土壤中难溶的养分，提供更多营养物质，更有利于土壤中微生物的生长繁殖[29]。

Zhou等[17]研究发现长期施用氮肥会显著降低土壤中16S rRNA基因的拷贝数。而Wessén等[30]研究得出施用无机肥硝酸钙和硫酸铵会提高16S rRNA基因的拷贝数。我们的研究表明，无机肥处理并没有降低基因拷贝数，施用有机肥，尤其是有机肥配施无机肥能显著提高基因拷贝数。且土壤中16S rRNA基因拷贝数与土壤化学性质有明显的相关性（与硝态氮、速效钾含量显著正相关，而与铵态氮含量负相关）。前面分析已知土壤化学性质同时受施肥和作物的影响，进一步表明黑土中微生物的丰度的变化是施肥和作物共同作用的结果。

3.2 不同作物和施肥处理对黑土中微生物多样性的影响

本研究发现在同一处理中，大豆季的土壤微生物多样性明显高于玉米季，如有机肥配施无机肥处理的Chao1指数大豆季比玉米季提高了16.60%。此外，Wu等[16]研究长期定位试验站的红壤水稻土，结果表明有机肥配施无机肥处理的土壤微生物多样性与无机肥处理之间没有显著性差异；而Tian等[31]发现大量施用有机粪肥较无机肥降低土壤微生物多样性。我们的研究与Doan等[32]的发现一致，即有机肥无机肥配施的微生物多样性高于无机肥。这可能与土壤类型和施肥的时间有关系，Tian是短期定位实验（3年），而本研究是基于35年以上的长期定位实验的结果。

本文的聚类分析表明，群落结构按照作物季分成明显的两大簇，说明其存在较大差异，同时，同一簇的4个不同施肥处理相互比较，得知不同施肥方式也引起了群落结构间的差异。与有机肥无机肥配施相比，施无机肥对黑土土壤微生物群落结构改变更大。因此，长期施肥导致的东北黑土中微生物群落的变化跟作物种类、气候和降水等因素相关，但各因素作用的程度如何还不得而知。同时，根据土壤微生物多样性指数与化学指标的相关性分析，可知多样性指数与土壤全氮、硝态氮、铵态氮、有效磷和速效钾等化学性质密切相关。

3.3 不同作物和施肥处理对黑土中微生物组成的影响

本研究表明在大豆和玉米季中的土壤微生物组成存在较大差异。玉米季的古菌所占比例（7.53%— 11.62%）高于大豆季（仅占2.70%—4.10%）；玉米季Acidobacteria是细菌中相对丰度最高的菌门，而大豆季相对丰度最高的菌门是Proteobacteria；玉米季的Bacteroidetes相对丰度较高而在大豆季较低，大豆季的Actinobacteria相对丰度较高而在玉米季较低。Bacteroidetes和Actinobacteria丰度在两季作物中差异明显，作物影响微生物组成，这个结果与以前的报道一致[13, 33]。且同一作物不同施肥处理之间也存在差异，如Alphaproteobacteria和Gammaproteobacteria相对丰度变化较大，施肥处理均比不施肥处理的相对丰度提高，对有机肥配施无机肥和无机肥均比较敏感。Actinobacteria的相对丰度不同施肥处理之间没有明显差异。Chloracidobacteria在不施肥处理具有较高相对丰度，而对有机肥配施无机肥、有机肥和无机肥均表现不稳定性。以上分析得知，作物和施肥均改变了土壤微生物的群落组成，组成的改变是否会引起土壤微生物整体代谢功能的变化，还有待于更深入的研究。

根据土壤微生物组成（门水平）与化学指标的相关性分析，可知微生物组成与土壤pH、全氮、硝态氮、铵态氮、有效磷和速效钾等化学性质密切相关。土壤质量的改变是缓慢的、长期积累的效应，土壤微生物随着土壤性质的改变而改变。而本论文没有关于施肥和作物引起微生物数量、多样性等变化的贡献量进行量化分析，以及施肥与作物交互作用机理等的研究，这些问题还需进一步研究探讨。

4 结论

本研究表明，施肥和不同作物均对东北黑土土壤微生物丰度、α多样性和群落结构产生了影响。连续施用无机肥会导致土壤酸化，而有机肥无机肥配施能有效地提高土壤pH，缓解酸化。且土壤化学性质与土壤微生物群落变化具有相关性。有机肥配施无机肥显著增加了黑土中微生物的16S rRNA基因拷贝数，而无机肥对其影响不显著；大豆种植季的16S rRNA基因拷贝数显著高于玉米种植季。有机肥配施无机肥的α多样性指数高于无机肥；大豆种植季的α多样性指数较玉米种植季的高。

References

[1] Yin C, Fan F, Song A, Cui P, Li T, Liang Y. Denitrification potential under different fertilization regimes is closely coupled with changes in the denitrifying community in a black soil., 2015, 99: 5719-5729.

[2] KrashevskaV, Klarner B, Widyastuti R, Maraun M, Scheu S. Impact of tropical lowland rainforest conversion into rubber and oil palm plantations on soil microbial communities., 2015, 51(6): 697-705.

[3] 乔洁,毕利东, 张卫建, 沈仁芳, 张斌, 胡锋, 刘艳丽. 长期施用化肥对红壤性水稻土中微生物生物量、活性及群落结构的影响. 土壤,. 2007, 39(5): 772-776.

Qiao J, Bi L D, Zhang W J, Shen R F, Zhang B, Hu F, Liu Y L. Effects of long-term chemical fertilization on soil microbial biomass, activity and community in paddy soil in red soil region of China., 2007, 39(5): 772-776. (in Chinese )

[4] Jorquera M A, Martinez O A, Marileo L G, Acuña J J, Saggar S, Mora M L. Effect of nitrogen and phosphorus fertilization on the composition of rhizobacterial communities of two Chilean Andisol pastures., 2014, 30, 99-107.

[5] Beauregard M S, Hamel C, Atul N, Arnaud M S. Long-term phosphorus fertilization impacts soil fungal and bacterial diversity but not AM fungal community in alfalfa., 2010, 59(2): 379-389.

[6] Börjesson G, Menichetti L, Kirchmann H, Kätterer T. Soil microbial community structure affected by 53 years of nitrogen fertilisation and different organic amendments., 2012, 48(3): 245-257.

[7] Berthrong T, Buckley D H, Drinkwater L E. Agricultural management and labile carbon additions affect soil microbial community structure and interact with carbon and nitrogen cycling., 2013, 66(1): 158-170.

[8] Wei D, Yang Q, Zhang J Z, Wang S, Chen X, Zhang X, Li W Q. Bacterial communities structure and diversity in a black soil as affected by long-term fertilization., 2008, 18 (5): 582-592.

[9] Carlson J, Saxena j, Basta N, Hundal L, Busalacchi D, Dick R P. Application of organic amendments to restore degraded soil: Effects on soil microbial properties., 2015, 187(3): 1-15.

[10] Peacock A D, Mullen M D, Ringelberg D B, Tyler D D, Hedrick D B, Gale P M, White D C. Soil microbial community responses to dairy manure or ammonium nitrate applications., 2001, 33: 1011-1019.

[11] 魏巍, 许艳丽, 朱琳, 韩晓增, Li S. 长期施肥对黑土农田土壤微生物群落的影响. 土壤学报, 2013, 50: 372-380.

Wei W, Xu Y L, Zhu L, Han X Z, Li S. Effectof long-term fertilization on soil microbial communities in farmland of black soil., 2013, 50: 372-380. (in Chinese )

[12] Liu Y, Shi G, Mao LCheng G, Jiang S, Ma X, An L, Du G, Johnson N C, Feng HDirect and indirect influences of 8 yr of nitrogen and phosphorus fertilization on Glomeromycota in an alpine meadow ecosystem., 2012, 194: 523-535.

[13] Zhao J, Ni T, Li Y, Xiong W, Ran W, Shen B, Shen Q, Zhang R. Responses of bacterial communities in arable soils in a rice-wheat cropping system to different fertilizer regimes and sampling times., 2014, 9(1): e85301.

[14] 蔡艳, 郝明德, 张丽琼, 臧逸飞, 何晓雁. 应用454测序技术分析种植制度对黑垆土微生物多样性的影响. 作物学报, , 2015, 41(2): 339-346.

Cai Y, Hao M D, Zhang L Q, Zang Y F, He X Y. Effect of cropping systems on microbial diversity in black loessial soil tested by 454 sequencing technology., 2015, 41(2): 339-346. (in Chinese )

[15] Kiely P D, Call D F, Yates M D, Regan J M, Logan B E. Anodic biofilms in microbial fuel cells harbor low numbers of higher-power-producing bacteria than abundant genera., 2010, 88(1): 371-380.

[16] Shokralla S, Spall J L, Gibson J F, Hajibabaei M. Next-generation sequencing technologies for environmental DNA research., 2012, 21(8): 1794-1805.

[17] Zhou J, Guan D, Zhou B, Zhao B, Ma M, Qin J, Jiang X, Chen S, Cao F, Shen D, Li J. Influence of 34-years of fertilization on bacterial communities in an intensively cultivated black soil in Northeast China., 2015, 90: 42-51.

[18] 秦杰, 姜昕, 周晶, 马鸣超, 关大伟, 周宝库, 赵百锁, 杜秉海, 李俊. 长期不同施肥黑土细菌和古菌群落结构及主效影响因子分析. 植物营养与肥料学报, 2015, 21(6): 1590-1598.

Qin J, Jiang X, Zhou J, Ma M C, Guan D W, Zhou B K, Zhao B S, Du B H, Li J. Characteristics and driving factors of soil bacterial and archaeal communities under long-term fertilization regimes in black soil., 2015, 21(6): 1590-1598. (in Chinese)

[19] Zhou J, Jiang X, Zhou B, Ma M, Guan D, Li J, Chen S, Cao F, Shen D, Qin J. Thirty four years of nitrogen fertilization decreases fungal diversity and alters fungal community composition in black soil in Northeast China., 2016, 95: 135-143.

[20] Wu M, Qin H, Chen Z, Wu J, Wei W. Effect of long-term fertilization on bacterial composition in rice paddy soil., 2011, 47(4): 397-405.

[21] Ding X, Han X, Zhang X, Qiao Y, Liang Y. Continuous manuring combined with chemical fertilizer affects soil microbial residues in a Mollisol., 2013, 49(4): 387-393.

[22] Mulvaney R L. Nitrogen-inorganic forms//Bigham J M.Soil Science Society of America, ASA, Madison , Wisconsin, USA, 1996: 1123-1200.

[23] 鲁如坤. 土壤农业化学分析方法. 北京: 中国农业科技出版社, 1999, 431-472.

Lu R K.. Beijing: China Agricultural Science and Technology Press, 1999: 431-472. (in Chinese)

[24] Peiffer J A, Spor A, Koren O, Jin Z, Tringe S G, Dangl J L, Buckler E S, Ley R E. Diversity and heritability of the maize rhizosphere microbiome under field conditions., 2013, 110(16): 6548-6553.

[25] Chinnadurai C, Gopalaswamy G, Balachandar D. Long term effects of nutrient management regimes on abundance of bacterial genes and soil biochemical processes for fertility sustainability in a semi-arid tropical Alfisol., 2014, 232: 563-572.

[26] Chen J, Bittinger K, Charson E S, Hoffmann C, Lewis J, Wu G D, Collman R G, Bushman F D, Li H. Associating microbiome composition with environmental covariates using generalized UniFrac distances., 2012, 28(16): 2106-2113.

[27] Yu X, Wu C, Fu Y, Brookes P C, Lu S. Three-dimensional pore structure and carbon distribution of macroaggregates in biochar- amended soil., 2016, 67(1): 109-120.

[28] Han X, Wang R, Liu J, Wang M, Zhou J, Guo W.Effects of vegetationtypes on soil microbial community composition and catabolicdiversity assessed by polyphasic methods in North China., 2007, 19(10): 1228-1234.

[29] 李晓慧, 韩晓增, 王树起, 王守宇, 彦君. 土壤微生物对不同作物根系活动的响应. 农业系统科学与综合研究, 2010, 26(2): 156-159.

Li X H, Han X Z, Wang S Q, Wang S Y, Yan J. Soil microorganism amount and population as affected by different crop roots., 2010, 26(2): 156-159. (in Chinese)

[30] WessénE, Hallin S, Philippot L. Differential responses of bacterial and archaeal groups at high taxonomical ranks to soil management., 2010, 42(10): 1759-1765.

[31] Tian W, Wang L, Li Y, Zhuang K, Li G, Zhang J, Xiao X, Xi Y. Responses of microbial activity, abundance, and community in wheat soil after three years of heavy fertilization with manure-based compost and inorganic nitrogen., 2015, 213: 219-227.

[32] Doan T T, Bouvier C, Bettarel Y, Bouvier T, Henry-des-Tureaux T, Janeau J L, Lamballe P, Van Nguyen B, Jouquet P. Influence of buffalo manure, compost, vermicompost and biochar amendments on bacterial and viral communities in soil and adjacent aquatic systems., 2014, 73: 78-86.

[33] Zhao J, Ni T, Li J, Lu Q, Fan f Z, Huang Q, Zhang R, Li R, Shen B, Shen Q. Effects of organic–inorganic compound fertilizer with reduced chemical fertilizer application on crop yields, soil biological activity and bacterial community structure in a rice–wheat cropping system., 2016, 99: 1-12.

（责任编辑 杨鑫浩）

Responses of Micropopulation in Black Soil of Northeast China to Long-Term Fertilization and Crops

DING Jian-li1,2, JIANG Xin1,2, GUAN Da-wei1,2, MA Ming-chao1,2, ZHAO Bai-suo2, ZHOU Bao-ku3, CAO Feng-ming1,2, LI Li1,2, LI Jun1,2

(1Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081;2Laboratory of Quality＆Safety Risk Assessment for Microbial Products, Ministry of Agriculture, Beijing 100081;3Institute of Soil Fertility and Environmental Sources, Heilongjiang Academy of Agricultural Sciences, Harbin 150086)

【Objective】This experiment was carried out to investigate the characteristics of micropopulation in northeast China black soil under more 35 years of fertilization of two crop seasons. The effects of inorganic fertilizer and manure on the abundance and structure of micropopulation were identified. Black soil micropopulation responses to the interaction of long-term fertilization and crops would be revealed. This study will provide evidences for further enhancing fertilization and tillage method.【Method】Based on a long-term fertilization experiment carried out in Heilongjiang Academy of Agricultural Sciences, four different fertilization treatment samples in soybean and maize growing seasons were selected including samples without fertilizer (CK); manure (M); inorganic fertilizer (NPK); and inorganic fertilizers with manure (MNPK). Letter “m” before treatment represented maize season sample, and letter “s” represented soybean season sample. The Illumina Miseq sequencing and quantitative PCR of 16S rRNA gene were applied to analyze the effects of fertilization and cropson micropopulation in black soil. Correlation analysis was carried out between micropopulation and the soil properties.【Result】The 16S rRNA gene copy numbers in maize growing season (6.32×108-8.83×108/ng DNA) was lower than that in soybean growing season (0.96×109-2.30×109/ng DNA). Alpha diversity in maize growing season (ACE index was between 3 674.58-4 034.84) was lower than that of soybean season (ACE index was between 4 167.47-4 887.36), too. The top phylum was Acidobacteria (24.47%-27.90%) in maize growing season, but it was Proteobacteria (27.78%-34.40%) in soybean growing season. The relative abundances of Bacteroidetes and Actinobacteria were significantly different between two crop growing seasons. The 16S rRNA gene copy numbers in treatment of inorganic fertilizers with manure was greater than that of inorganic fertilizers. Alpha diversity in treatment of inorganic fertilizers with manure was higher than that of inorganic fertilizers (Chao1 index of sMNPK compared with that of sNPK was more 11.89%). The composition of micropopulation in different fertilization treatments of one crop growing season was different. The relative abundances of Alphaproteobacteria in sMNPK and sNPK compared with that of sCK were more 3.31% and 5.24%, Gammaproteobacteria in sMNPK and sNPK were higher 1.72% and 1.20% than that in sCK, and were sensitive to fertilizers. Correlation analysis showed that 16S rRNA gene copy number was positively correlated with soil NO3--N and available K, the diversity index and soil total N, NO3--N, NH4+-N, available P and available K were closely related. 【Conclusion】Results of the research demonstrated that of long-term different fertilizations and different crop growing seasons have effects on microbial richness, α diversity and community structure. Inorganic and organic fertilizers improved the soil pH, slowed down soil acidification, changed microbial structures, increased microbial richness and diversity as well as the metabolic activity of micropopulation.

fertilizer; crop; microbial community structure; soil chemical properties; high-throughput sequencing; black soil

2016-05-13；接受日期：2016-09-06

国家自然科学基金（41573066）、国家“973”计划（2015CB150506）、国家“863”计划（2013AA102802-04）、农业微生物产品质量安全风险评估专项任务（GJFP2016013）、国家现代农业产业技术体系建设专项（CARS-04）

丁建莉，E-mail：lndjl@126.com。通信作者姜昕，E-mail：jiangxin@caas.cn。通信作者李俊，E-mail：lijun01@caas.cn