饲料油菜作绿肥对后茬麦田土壤肥力及细菌群落的影响

李文广，杨晓晓，黄春国，薛乃雯，夏清，刘小丽，张晓琪，杨思，杨珍平，高志强

饲料油菜作绿肥对后茬麦田土壤肥力及细菌群落的影响

李文广，杨晓晓，黄春国，薛乃雯，夏清，刘小丽，张晓琪，杨思，杨珍平，高志强

（山西农业大学农学院，山西太谷 030801）

【】探讨黄土高原旱地麦后复种饲料油菜还田对后茬麦田土壤养分、酶活性及细菌群落的影响，为改善农田地力提供依据。【】设置饲料油菜播量（S1：小播量；S2：中播量；S3：大播量）和还田时期（D1：早期还田9月10日；D2：中期还田9月20日；D3：晚期还田9月30日）二因素试验，以常规农户复种夏玉米模式为对照，比较后茬麦田土壤养分及酶活性变化，并采用16S rRNA 基因Illumina MiSeq高通量测序技术和PICRUSt基因预测分析方法，比较后茬麦田耕层土壤细菌群落结构、多样性及代谢功能变化。【】与常规农户模式相比，饲料油菜还田可不同程度提升土壤养分及酶活性，其中有机质含量和蔗糖酶活性的提升范围最大，分别为11.73%—60.5% 和1.4%—94.5%；饲料油菜播量显著影响土壤养分和酶活性，还田时期显著影响土壤有机质含量和蔗糖酶活性，二者间互作效应显著影响碱性磷酸酶活性；播量、还田时期及其互作效应均显著影响细菌群落多样性；大播量晚期还田处理（S3D3）在所有指标中均表现最优。细菌群落组成方面，门水平上10个土壤样品共得到19个类群，其中优势菌群为变形菌门、放线菌门、酸杆菌门和厚壁菌门；聚类结果显示，10个处理可聚为大播量、中播量、小播量和常规农户模式4类，表明不同播量油菜还田明显改变了土壤细菌群落构成。RDA分析表明，土壤养分、酶活性与细菌群落多样性密切相关。PICRUSt功能预测表明，土壤细菌群落具有丰富的代谢功能，其中编码新陈代谢的基因具有明显优势；二级预测功能分类中，氨基酸代谢、碳水化合物代谢和能量代谢是主要代谢功能路径。KEGG直系同源基因簇KO（KEGG orthologous groups）丰度热图结果表明，油菜还田增加了与土壤碳、氮代谢相关的细菌群落。【】合理复种饲料油菜并还田可提高后茬麦田土壤养分含量及酶活性，并有效改善细菌群落多样性，增加土壤有益菌群落数量。夏闲期引入饲料油菜作绿肥对提高黄土高原旱地麦田土壤肥力以及合理农作非常有利。

油菜；还田；绿肥；小麦；土壤肥力；细菌群落；高通量测序；功能预测

0 引言

【研究意义】黄土高原旱地是我国西北地区重要的耕地资源，但其土壤全氮含量和有机质含量相对较低，超过一半的土壤有机质含量低于1%[1]。为提高作物产量，往往过量化学氮肥被使用，最终导致土壤酸化板结，微生物多样性降低[2]。作为化学肥料零增长[3]计划和土壤有机质提升计划的重要替代肥源，还田后的绿肥在腐解过程中能释放大量养分[4]，从而提高土壤肥力，增加作物产量[5-6]，对农业可持续发展有着重要意义。土壤微生物作为农田生态系统的重要组成部分，对土壤生态变化十分敏感，易受气候、温度、土壤湿度及施肥等外界因素影响，是评价施肥和耕作制度改进土壤肥力和作物生产力的重要指标[7-9]。施化肥或绿肥均可改变土壤细菌群落结构及数量，进而导致土壤结构及肥力的改变[10]。饲料油菜作为绿肥优势明显，其植株生物量大，盛花期养分含量较其他绿肥高，其粗长的直根系能对土壤形成穿刺效应，破除犁底层，改善土壤通气状况，提高土壤有氧呼吸效率及酶活性[11-12]。因此，开展将饲料油菜作为绿肥还田提高土壤肥力、改善土壤细菌群落结构及其多样性研究具有重要意义。【前人研究进展】由于旱地小麦种植区年降水量不稳定，降水季节分布不均，过去常采用夏季休耕方式来稳定冬小麦产量[13]，但该法会加重黄土高原地区土壤侵蚀。国内外一些研究均表明[14-15]，夏闲期复种绿肥能有效提高旱地土壤蓄水效率，使小麦增产，可代替夏季休耕。李红燕等[16]研究发现，麦后复种长武怀豆和油菜并还田，均可显著提高黄土高原南部土壤有机质、速效磷含量及磷酸酶、蔗糖酶和脲酶活性。何红霞等[17]研究发现，休闲期种植长武怀豆能有效培肥黄土高原旱地土壤肥力。持续干旱条件可能使旱地小麦种植区土壤微生物生物量逐渐下降，而绿肥还田不仅可以增加土壤蓄水能力，还能有效增加土壤微生物活性[18-19]。LONGA等[20-21]发现，绿肥是土壤微生物多样性的重要来源，绿肥还田能显著改变土壤微生物群落组成及丰富度。【本研究切入点】本试验所在旱地小麦种植区，雨水多集中于麦后休闲期，该时期光热资源充足，但该区土壤肥力偏低，不合理复种会导致土壤肥力持续下降。利用饲料油菜作绿肥能充分利用麦后休闲期的水、热、土资源，提高土壤使用可持续性。目前，对绿肥的研究多集中于品种间比较及其对土壤理化性质的影响，鲜有研究将饲料油菜引入北方旱地麦田，探讨还田管理方式对土壤培肥及合理农作的影响。【拟解决的关键问题】本研究应用16S rRNA 基因Illumina MiSeq 高通量测序技术分析细菌群落及代谢功能变化情况，探讨后茬麦田土壤养分、酶活性变化，以期明确饲料油菜还田对土壤养分、酶活性和细菌群落的影响，确定最佳播量及还田时期，为提高旱地麦田土壤肥力，维护微生物生态系统和旱地农业可持续发展提供科学依据。

1 材料与方法

1.1 研究区域与样地概况

试验于2015年6月至2016年6月在山西省运城市垣曲县长直乡鲁家坡村十倾园（111°43.3′E，35°14.4′N）进行。该区年均气温13.5℃，年均日照时数2 026.2 h，年积温4 900℃，全年无霜期236 d左右，降水集中在5—9月份，年降水量600—800 mm，年蒸发量1 200 mm，是典型旱作雨养农业区，属暖温带半湿润大陆性季风气候。供试土壤质地为中壤土，土壤类型属于褐土性红立黄土，pH为8.0。饲料油菜播前测得0—20 cm耕层土壤基本性质为：有机质11.51 g·kg-1，全氮0.87 g·kg-1，全磷0.81 g·kg-1，碱解氮41.13 mg·kg-1，速效磷14.58 mg·kg-1。

1.2 供试品种

饲料油菜品种：华油杂62号，长势中等、含油量高，抗病抗倒伏能力强，由华中农业大学油菜研究室提供。

小麦品种：烟农21号，具有较强的抗寒、抗旱和抗倒伏能力，由侯马市金色农田有限公司提供。

1.3 试验设计

前茬小麦收获后，复种饲料油菜，考察合理的种植密度；后茬小麦播种前，将饲料油菜还田，考察适宜的还田时间。试验采用二因素裂区设计：主区为播量（S），设置7.5 kg·hm-2（小播量，S1）、15 kg·hm-2（中播量，S2）和22.5 kg·hm-2（大播量，S3）3个播量水平。副区为还田时期（D），设置9月10日（早期还田，D1）、9月20日（中期还田，D2）和9月30日（晚期还田，D3）3个不同还田时期，以常规农户的麦后复种夏玉米模式（玉米收获后秸秆全量还田）作为对照。各处理3个小区重复，共计30个小区，各小区面积为15 m×44 m=660 m2。

2016年6月10日麦收后按试验设计复种饲料油菜，条播，行距20 cm。于9月10日、9月20日、9月30日3次对设定小区分别进行饲料油菜翻压还田，10月9日采用施肥播种一体机种植冬小麦，条播，行距20 cm，种植密度112.5 kg·hm-2，所施氮磷钾复合肥用量分别为纯N 135 kg·hm-2、P2O5165 kg·hm-2、K2O 37.5 kg·hm-2，施肥深度为20—40 cm。常规农户模式的麦后复种夏玉米种植采用机械化种肥同穴播种，肥料用量是氮磷钾复合肥750 kg·hm-2，肥料施于距种子10 cm处。

其他管理措施同一般大田生产。

1.4 土壤样品采集与处理

2017年6月1日小麦收割前采集土壤样品。每个试验小区取3个样方作为区内重复，每个样方内随机选取5个点，使用直径为4 cm的灭菌土钻采集0—20 cm的耕层土壤土样，将采集的5个样点的土样去除根系、动植物残体、石块等杂物后混合均匀迅速装入无菌自封袋，暂存于4℃冰盒，带回实验室后将混合土样分为3份，其中2份风干后分别用于土壤养分及酶活性测定，1份放置于-80℃冰箱保存用于土壤DNA提取。

1.5 测定项目及方法

1.5.1 土壤养分及酶活性测定 采用常规方法[22]测定土壤化学性质：有机质采用重铬酸钾外加热法；碱解氮采用碱解扩散法；速效磷采用0.5 mol·L-1NaHCO3浸提-钼锑抗比色法；速效钾采用1.0 mol·L-1NH4OAc 浸提-火焰光度法测定；脲酶活性采用靛酚比色法测定；碱性磷酸酶活性采用磷酸苯二钠比色法测定；蔗糖酶活性采用3,5﹣二硝基水杨酸比色法测定[23]。

1.5.2 饲料油菜还田生物量 还田当日，选取3个具有代表性的、面积为1 m2的样段，将收获鲜植株地上部分放入烘箱，于105℃杀青，30min后再80℃烘至恒重，记录各样段干物质重。

1.5.3 细菌基因组DNA的提取及测序 采用Fast DNA SPIN Extraction Kits 试剂盒（MP Biomedicals，Santa Ana，CA，美国），提取土壤总DNA。使用NealDel- ND1000分光光度计（Thermo Fisher Scientific，Waltham，MA，USA）和琼脂糖凝胶电泳分别测定DNA浓度和纯度。

对所提取基因组DNA的16S rRNA V4区进行扩增，通用引物序列为：520F（AYTGGGYDTAAA GNG），802R（TACNVGGGTATCTAATCC）。PCR反应条件为：98℃预变性2 min；25个循环；98℃变性15 s；55℃退火30 s；72℃延伸30 s；72℃后延长5 min，于4℃保存。PCR扩增产物使用Agencourt AMPure Beads（Beckman Coulter，Indianapolis，IN）试剂盒进行纯化，使用PicoGreen dsDNA Assay Kit（Invitrogen，Carlsbad，CA，USA）检测试剂盒进行定量。样品委托上海派森诺生物科技有限公司采用Illumina MiSeq高通量测序技术平台进行序列测序和分析。

1.5.4 测序数据优化处理 对原始数据进行质量过滤和双端序列连接，运用QIIME（version 1.9.0,//qiime. org/）进行序列过滤，数据过滤标准为[24]：去除5′端引物错配碱基数＞1的序列；去除含有N（模糊碱基）序列；去除含有连续相同碱基数＞8的序列；去除长度≤150 bp的序列。运用MOTHUR软件（version1.31.2，http://www.mothur.org/）中UCLHIME方法去除嵌合体序列，得到最终用于后续分析的优质序列。

1.5.5 OTU聚类分析及注释 在QIIME中调用UCLUST方法对优质序列按序列相似度97%进行聚类，选取每个类中最长的序列为代表序列。在QIIME中调用BLAST方法对序列数据库进行比对，获得每个OTU代表序列的分类学信息。注释数据库为：Greengenes（Release 13.8，http://greengenes. secondgenome.com/）。之后对OTU进行精简处理，去掉丰度值小于总序列数条数0.001%的OTU，得到后续分析使用的精简后的OTU列表。

1.5.6 PICRUSt功能预测 功能和代谢途径预测采用PICRUSt 软件进行分析[25]，利用QIIME获得的closed OTU-table 与KEGG 数据库进行比对，获得不同细菌群落功能预测信息。

1.6 计算和统计方法

根据OTU列表中的各样品物种丰度情况，应用软件MOTHUR中的summary.Single命令，计算菌群Chao、ACE和Shannon指数[26-27]。使用QIIME软件，对Weighted的UniFrac距离矩阵进行UPGMA聚类分析，并使用R软件进行可视化。采用Excel 2003进行基本计算，并制作各水平群落分类学组成和差异分析图。采用SAS9.4软件的ANOVA过程进行二因素裂区设计方差分析及多重比较（文中多重比较含CK，方差分析不含CK），差异显著性水平为＜0.05；使用CANOCO4.5软件对土壤化学性质和细菌群落结构及多样性进行冗余分析（RDA）。

2 结果

2.1 饲料油菜作绿肥还田对后茬麦田土壤养分及酶活性的影响

通过分析饲料油菜还田对后茬麦田土壤养分及酶活性变化情况（表1）可知。与常规种植方式相比，所有饲料油菜还田处理均改变了后茬麦田各样地土壤养分含量及酶活性，其中，土壤中速效磷、碱解氮和有机质含量的增加范围为3.15%—18.5%、3.33%—30.2%和11.73%—60.5%；脲酶、碱性磷酸酶和蔗糖酶活性的增加范围为3.12%—31.25%、4.84%—25.8%和1.4%—94.5%。各指标中，大播量晚期还田（S3D3）处理表现出明显优势，这可能与饲料油菜还田生物量（表1）及腐解速率有关。对9个饲料油菜还田处理作二因素方差分析，F检验表明，播量对土壤养分及酶活性的影响均达到极显著水平（＜0.01），还田时期仅对有机质及蔗糖酶的影响达到显著水平（＜0.05），播量与还田时期互作则对土壤碱性磷酸酶活性的影响达到显著水平（＜0.05）。土壤养分及酶活性均随播量增加及还田时期推迟而增加，以S3与D3表现最优。

2.2 饲料油菜作绿肥还田对后茬麦田土壤细菌群落多样性的影响

饲料油菜还田提高了土壤细菌群落多样性（表2），播量、还田时期以及二者交互作用对Chao 1、ACE、Shannon指数及OTU的影响均达到显著（＜0.05）或极显著（＜0.01）水平。OTU及多样性指数均随饲料油菜播量增加而增加，S3显著高于S1与S2。随还田时期推后，OTU及多样性指数表现出先降低后增加的趋势。

2.3 饲料油菜作绿肥还田对后茬麦田土壤细菌群落组成的影响及差异分析

2.3.1 门水平的细菌群落结构特征 各处理细菌群落在门分类水平上具有较高的多样性。共得到19个类群，其中相对丰度≥1%的有12个，分别为变形菌门（Proteobacteria，13.1%—28.5%）、放线菌门（Actinobacteria，15.2%—30.7%）、酸杆菌门（Acidobacteria，8.4%—21.4%）、厚壁菌门（Firmicutes，7.5%—21.6%）、芽单胞菌门（Gemmatimonadetes，4.7%—11.7%）、硝化螺旋菌门（Nitrospirae，0.8%—9.2%）、浮霉菌门（Planctomycetes，3.7%—6.1%）、绿弯菌门（Chloroflexi，3.6%—5.6%）、GAL15（0%—6.1%）、疣微菌门（Verrucomicrobia，1.5%—3.2%）、拟杆菌门（Bacteroidetes，0.7%—3.9%）、WS3（0.2%—2.4%），它们构成了土壤样品在门水平上的基本结构。其余7个门的平均占比均小于1%，共占1.9%。

对门水平相对丰度前10的菌群进行方差分析（图1），结果表明，变形菌门相对丰度随播量增加而增加，S3D3显著高于其他样地（＜0.05）；厚壁菌门、硝化螺旋菌门相对丰度随播量增加呈现降低趋势；CK的放线菌门相对丰度显著高于其他样地（＜0.05）；饲料油菜还田显著提高土壤酸杆菌门、芽单胞菌门、硝化螺旋菌门、GAL15和疣微菌门的相对丰度（＜0.05）。

表1 饲料油菜播量（S）及还田时期（D）对后茬麦田土壤养分及酶活性的影响

同列不同小写字母表示在5%水平上差异显著；多重比较含CK，方差分析不含CK；*表示在0.05水平上差异显著；**表示在0.01水平上差异显著。下同

Values followed by different letters are significantly different at 0.05 level; Multiple comparisons contain the value of CK, ANOVA does not contain the value of CK; * denote significantly different at 0.05 levels; **denote significant difference at 0.01 level. The same as below

2.3.2 纲水平的微生物群落分布特征 各处理细菌群落隶属63个纲，相对丰度≥2%的有15个。其中相对丰度≥2%的有芽孢杆菌纲（Bacilli，7.5%—16.9%），酸杆菌纲（Acidobacteria，4.5%—14.7%），α-变形菌纲（Alphaproteobacteria，5.6%—12.4%），放线菌纲（Actinobacteria，5.2%—17.8%），嗜热油菌纲（Thermoleophilia，4%—8.7%），硝化螺旋菌纲（Nitrospira，0.8%—9.2%）和β-变形菌纲（Betaproteobacteria，2.5%—7.8%）；其次为γ-变形菌纲（Gammaproteobacteria，1.5%—6.6%），芽单胞菌纲（Gemmatimonadetes，1.2%—6.2%），芽单孢菌门的Gemm-1纲（0.8%—5%）、GAL15纲（0—7%），δ变形菌纲（Deltaproteobacteria，2.6%—3.8%），浮霉菌门的Planctomycetia纲（2.2%—3.8%），放线菌门的MB-A2-108纲（0.6%—4.1%）和酸杆菌门的Chloracidobacterial纲（1.3%—2.9%）。他们构成了纲水平上的基本结构。其余48纲平均相对丰度共占20%。

表2 饲料油菜播量（S）及还田时期（D）对后茬麦田土壤细菌群落多样性的影响

对纲水平相对丰度前10的细菌类群进行方差分析（图2）发现，饲料油菜还田显著增加了酸杆菌纲、硝化螺旋菌纲和芽单胞菌门下Gemm-1纲的相对丰度（＜0.05），显著降低了放线菌纲、嗜热油菌纲和β-变形菌纲的相对丰度（＜0.05）；α-变形菌纲、β-变形菌纲、γ-变形菌纲和芽单孢菌纲的相对丰度随饲料油菜播量增加而增加，以S3D3处理的相对丰度最高。

2.4 不同样本土壤细菌群落与土壤养分、酶活性之间的相关性

基于OTU层次的UPGMA聚类分析表明（图3），在0.075的相似性水平上，包括CK在内的不同样品可聚为4类，这表明，与CK相比，饲料油菜在不同播量及还田时期互作下还田改变了后茬麦田耕层土壤细菌群落结构。

通过对土壤养分、酶活性与细菌群落结构纲分类水平及多样性间的RDA冗余分析发现（图4），土壤养分和酶活性均与Betaproteobacteria、Gammaproteobacteria、Alphaproteobacteria、Gemmatimonadetes及细菌群落多样性间呈显著正相关；Acidobacteria与碱性磷酸酶、碱解氮、速效磷及蔗糖酶呈正相关；Actinobacteria 与脲酶及有机质呈正相关；Bacilli、Thermoleophilia、Nitrospira和Gemm-1则与土壤化学性质呈负相关。表明土壤养分、酶活性和细菌群落三者间密切相关，在土壤物质能量循环中均发挥着重要作用。饲料油菜作绿肥的合理还田模式有利于促进土壤酶活性提高、改善土壤细菌群落多样性，进而提高土壤肥力。

图1 不同土壤样本相对丰度前10的菌群（门水平）

2.5 不同土壤样本细菌群落PICRUSt功能预测分析

2.5.1 功能基因家族组成 将预测得到的功能基因家族（图5）与KEGG数据库比对发现，第一等级上，所有菌群的基因序列注释到的功能可分为6类，其中前3类分别是：新陈代谢、环境信息加工和遗传信息加工。第二等级上，参与新陈代谢的菌群中，以参与氨基酸代谢、碳水化物代谢和能量代谢的为最大类群；参与环境信息加工的菌群中，参与膜运转类群具有明显优势，其次是参与信号转导的类群；参与遗传信息加工的菌群中，参与遗传信息复制和修复机制的类群具有明显优势。

2.5.2 代谢相关基因分析 图6为所测土壤样品KEGG直系同源基因簇（KEGG orthologous groups）丰度热图，横坐标表示所测土壤样本，纵坐标表示代谢相关基因所对应的酶。分析发现，氨基酸代谢、碳水化合物代谢和能量代谢通路上的相关基因表达差异较为显著；各土壤样品可大致聚为CK和饲料油菜还田（S2D3除外）两类，饲料油菜还田明显提高代谢相关基因丰度，其中S3D3在基因丰度上更具优势。对该处理下的优势基因所对应的酶进行筛选并归类，得到：与脂肪酸合成代谢相关的酰基载体蛋白还原酶（K00059），与能源物质代谢相关的重要中间代谢物乙酰CoA（K00626），具有特异催化蛋白质底物功能（丝氨酸、苏氨酸羟基磷酸化）的蛋白激酶（K08884），具有特异性识别天冬氨酸和谷氨酰胺并与相应tRNA形成氨酰tRNA的氨酰tRNA合成酶（K02433），作用于支链氨基酸（亮氨酸、缬氨酸和异亮氨酸）转运系统中的酶（K01995、K01998、K01999），参与脂肪酸β-氧化的烯酰CoA水合酶（K01692），参与脂肪酸合成与分解的活性代谢中间产物酰基CoA合成酶（K01897），以及参与肽/镍转运系统底物结合和渗透的酶蛋白（K02032、K02033、K02034、K02035）。

图2 不同土壤样本前10个细菌类群的相对丰度（纲水平）

图3 不同土壤样本基于UniFrac加权的细菌群落结构UPGMA聚类分析

Table 3 The UPGMA cluster analysis based on weighted UniFrac for different soil samples bacterial community

图4 土壤养分、酶活性与细菌群落纲分类水平及多样性指数的冗余分析

图5 不同土壤样本细菌群落功能基因KEGG丰度图

图6 不同土壤样本细菌群落KEGG直系同源基因簇（KO）丰度热图

3 讨论

3.1 土壤养分及酶活性变化

油菜作绿肥翻压还田可增加后茬麦田耕层土壤酶活性[16,28]，随绿肥还田量增大，土壤酶活性增加变缓[29]。本试验中，与常规农户模式相比，9个饲料油菜还田处理均增加了后茬麦田耕层土壤养分及酶活性，与前人研究结果一致，不同的是，随播量增加、还田时期推迟土壤养分及酶活性均呈逐渐增加趋势，以大播量晚期还田处理S3D3的土壤养分及酶活性最高。这可能与还田油菜生物量及腐解速率有关。本试验表明，随播量增加，还田生物量明显增加；随还田时期推后，还田生物量虽没有明显变化，但气温明显降低，可能会降低土壤微生物及土壤酶活性，进而导致还田油菜腐解速率变缓。有研究还指出，绿肥腐解速率表现为前期快、后期慢，腐解高峰均出现在还田后最初的1个月内[30]。本研究中，绿肥还田时期为9月10、20、30日，对应平均气温分别为15、11、7℃，因此，随气温降低，晚期还田油菜较前、中期还田油菜腐解速率变慢，有利于养分缓慢释放[31]，进而促进了后茬麦田土壤养分及酶活性提高，尤其与碳代谢相关的土壤有机质和蔗糖酶活性显著提高。

3.2 土壤细菌群落结构及多样性分析

LONGA等[20]研究指出，绿肥还田能增加土壤细菌Chao1、Shannon和Simpson指数。本研究结果表明，饲料油菜作绿肥时，其播量、还田时期及二者的互作效应均显著影响细菌群落多样性，以大播量晚期还田处理的多样性指数较高且高于CK，因此应合理规划播量及还田时期。本研究中，多样性指数在中低播量处理中出现低于CK的情况，可能的原因是该播量下还田油菜生物量低，对土壤微生物及土壤酶活性的促进作用低，向土壤释放养分低（表1）。而大播量晚期还田处理多样性指数较高，一方面缘于还田生物量高，有利于活化土壤微生物、酶活性，释放土壤养分；另一方面受低气温影响，还田油菜腐解速率变慢，其中未腐解油菜可能与后茬小麦构成了混作体系，进而增加了土壤细菌群落多样性及种群数量[31-32]。

本研究进一步表明，饲料油菜还田改善了土壤细菌群落组成且有益细菌丰度显著增加。各样地优势菌门为变形菌门、放线菌门、酸杆菌门、厚壁菌门。前人研究表明，绿肥还田后，土壤中酸杆菌门、芽单胞菌门和硝化螺旋菌门相对丰度增加[33-35]，这与本研究结果一致。变形菌门属需营养类群，其相对丰度在高氮及高有机质水平下通常增加[36-37]。本研究中，随播量增加变形菌门相对丰度逐渐增加，说明油菜还田为变形菌门提供了丰富的有机养分。酸杆菌门作为嗜酸菌具有降解植物残体多聚体、参与铁循环、增强光合作用及参与单碳化合物代谢等功能，其相对丰度也可反映土壤酸性条件[38-39]。本研究中，油菜还田显著增加了酸杆菌门的相对丰度，说明油菜还田可能降低了后茬麦田土壤pH。芽单胞菌门对土氮素循环中反硝化作用具有重要意义[40]。本研究中，芽单胞菌门相对丰度显著增加，这可能与油菜还田后土壤中具有合适的碳氮比有关。硝化螺旋菌门适合高pH环境条件生存[41]，本研究中硝化螺旋菌门相对丰度随着绿肥还田量的增加而逐渐降低，可能也是由于绿肥还田导致土壤pH降低引起的，而这种改变有利于改良黄土高原碱性土壤环境。疣微菌门为寡营养性细菌，其分布广泛，兼性好氧，能够利用各种糖类，其相对丰度显著增加可能是由于油菜还田提供大量的营养物质[42-43]。另外，本研究发现油菜还田降低了放线菌门及具有致病性的厚壁菌门的相对丰度。这与林叶春等研究结果一致[34]。关于上述菌群在油菜还田后可能促进土壤碳氮循环及降低土壤pH值中发挥的作用，将在后续研究中进一步探讨。

RDA结果表明，土壤优势菌群中的α-变形菌纲、β-变形菌纲、γ-变形菌纲、酸杆菌纲和芽单孢菌纲与土壤养分、酶活性相关性较强，说明油菜还田后土壤养分的增加有益于改善土壤有益细菌群落的多样性[44]。而硝化螺旋菌纲、厚壁菌门中的芽孢杆菌纲、嗜热油菌纲则与之呈负相关，这也与上述门水平上细菌群落变化相吻合。

3.3 功能预测

LEBRUN[45]研究表明，麦田根际土壤中存在大量具有相对特定生态功能的微生物，并进行着诸多活跃的与新陈代谢相关的活动。本研究发现，油菜还田后，后茬麦田土壤中编码新陈代谢的相关基因具有明显优势，由此推论，土壤化学性质的改变导致相对特定生态功能微生物改变。

KEGG丰度热图初步得出，S3D3处理的多种酶基因相对丰度具有明显优势。研究表明[46]，酰基载体蛋白还原酶能从头合成脂肪酸，烯酰CoA水合酶催化脂肪酸氧化分解中β氧化循环，酰基CoA合成酶在脂肪酸合成和分解代谢中也扮演重要角色[47]。乙酰CoA是微生物碳代谢的关键分子，是体内能源物质代谢的枢纽性物质[48]。本研究上述酶基因丰度较CK明显增加，表明油菜还田能增加土壤中与碳代谢有关的细菌类群，增强了碳源利用能力。当土壤中氨含量较高时可形成谷氨酰胺或天冬酰胺，通过不同的氨基酸转移酶形成多种氨基酸[49]。氨基酸能促进土壤微生物活动，增加土壤微生物数量，有利于加速有机物矿化，促进营养元素释放[44]。其中，支链氨基酸是一类具有多种生理和代谢作用的必需氨基酸[50]。氨基酸合成需要的碳骨架来源于糖酵解、光合作用、氧化磷酸戊糖途径及三羧酸循环[51]，因此，氨基酸的合成是碳代谢与氮代谢的枢纽。本研究显示，油菜还田增加了土壤中与碳代谢及氮代谢相关细菌，促进土壤碳、氮代谢。此外，镍元素能促进氮素转运，是脲酶和氢化酶活性表达所必需的金属辅基，缺乏会导致植物体内尿素的积累[52]，钛元素对增加叶绿素含量、提高光合效率、提高植物体内多种酶的活性有积极作用[53]。本研究发现，S3D3下与钛、镍元素转运相关的酶基因丰度，有利于增加土壤中微量元素钛和有益元素镍的利用。

4 结论

相比常规农户麦后复种夏玉米的种植方式，黄土高原旱地小麦种植区夏闲期引入饲料油菜作绿肥能提高后茬麦田土壤肥力，有效改善土壤细菌群落结构及多样性，增强土壤细菌碳代谢及氮代谢能力，对旱地农业可持续发展及合理农作具有重要意义。

[1] LIU Z P, SHAO M A, WANG Y Q. Effect of environmental factors on regional soil organic carbon stocks across the Loess Plateau region, China., 2011, 142(3): 184-194.

[2] 张卫峰, 马林, 黄高强, 武良, 陈新平, 张福锁. 中国氮肥发展、贡献和挑战. 中国农业科学, 2013, 46(15): 3161-3171.

ZHANG W F, MA L, HUANG G Q, WU L, CHEN X P, ZHANG F S. The development and contribution of nitrogenous fertilizer in China and challenges faced by the country., 2013, 46(15): 3161-3171. (in Chinese)

[3] 金书秦, 邢晓旭. 农业面源污染的趋势研判、政策评述和对策建议. 中国农业科学, 2018, 51(3): 593-600.

JIN S Q, XING X X.Trend analysis, policy evaluation, and recommendations of agricultural non-point source pollution., 2018, 51(3): 593-600. (in Chinese)

[4] ZHANG D B, YAO P W, ZHAO N, YU C W, CAO W D, GAO Y J.Contribution of green manure legumes to nitrogen dynamics in traditional winter wheat cropping system in the Loess Plateau of China., 2016, 72: 47-55.

[5] 陈洪俊, 黄国勤, 杨滨娟, 王晓维. 冬种绿肥对早稻产量及稻田杂草群落的影响. 中国农业科学, 2014, 47(10): 1976-1984.

CHEN H J, HUANG G Q, YANG B J, WANG X W.Effects of different winter planting-green manure on the grain yield of rice and weed community of paddy field., 2014, 47(10): 1976-1984. (in Chinese)

[6] 杨滨娟, 黄国勤, 王超, 林青, 徐宁. 稻田冬种绿肥对水稻产量和土壤肥力的影响. 中国生态农业学报, 2013, 21(10): 1209-1216.

YANG B J, HAUNG G Q, WANG C, LIN Q, XU N. Effects of winter green manure cultivation on rice yield and soil fertility in paddy field., 2013, 21(10): 1209-1216.(in Chinese)

[7] RAMIREZ K S, LAUBER C L, KNIGHT R, BRADFROD M, FIERER N. Consistent effects of nitrogen fertilization on soil bacterial communities in contrasting systems., 2010, 91(12): 3463-3470.

[8] DING J, JIANG X, GUAN D, ZHAO B, MA M, ZHOU B. Influence of inorganic fertilizer and organic manure application on fungal communities in a long-term field experiment of Chinese mollisols., 2017, 11: 114-122.

[9] 王传杰, 肖婧, 蔡岸冬, 张文菊, 徐明岗. 不同气候与施肥条件下农田土壤微生物生物量特征与容量分析. 中国农业科学, 2017, 50(6): 1067-1075.

WANG C J, XIAO J, CAI A D, ZHANG W J, XU M G.Capacity and characteristics of soil microbial biomass under various climate and fertilization conditions across china croplands., 2017, 50(6): 1067-1075. (in Chinese)

[10] CHANG H L, PARK K D, JUNG K Y, ALI M A, KIM P J. Effect of Chinese milk vetch (l.) as a green manure on rice productivity and methane emission in paddy soil.2015, 138(3): 343-347.

[11] YE X F, LIU H E, LI Z, WAGN Y, WAGN Y Y, WANG H F, LIU G S. Effects of green manure continuous application on soil microbial biomass and enzyme activity., 2014, 37(4): 498-508.

[12] NICOLAISEN M H, RISGAARD P N, REVSBECH N P, REICHARDTW, RAMSING N B. Nitrification–denitrification dynamics and community structure of ammonia oxidizing bacteria in a high yield irrigated philippine rice field., 2004, 49(3): 359-369.

[13] TANAKA D L, LYON D J, MILLER P R, MERRILL S D, MCCONKEY B G. Soil and water conservation advances in the semiarid Northern Great Plains., 2010, 82(2): 132-139.

[14] O'Dea J K, Miller P R, JONES C A. Greening summer fallow with legume green manures: On-farm assessment in north-central Montana., 2013, 68(4): 270-282.

[15] 赵刚, 樊廷录, 李尚中, 张建军, 王勇,党翼, 王磊. 夏休闲期复种油菜对旱地土壤水分和小麦产量的影响. 应用生态学报, 2013, 24(10): 2807-2813.

ZHAO G, FAN T L, LI S Z, ZHANG J J, WANG Y, DANG Y, WANG L. Effects of rape cropping in summer fallow period on dryland soil moisture content and winter wheat yield.2013, 24(10): 2807-2813. (in Chinese)

[16] 李红燕, 胡铁成, 曹群虎, 鱼昌为, 曹卫东, 黄冬琳, 翟丙年, 高亚军. 旱地不同绿肥品种和种植方式提高土壤肥力的效果. 植物营养与肥料学报, 2016, 22(5): 1310-1318.

LI H Y, HU T C, CAO Q H, YU C W, CAO W D, HUANG D L, ZHAI B N, GAO Y J. Effect of improving soil fertility by planting different green manures in different patterns in dryland., 2016, 22(5): 1310-1318.(in Chinese)

[17] 何红霞, 李小涵, 包明, 李超, 马小龙, 何刚, 邱炜红, 王朝辉. 栽培模式对旱地小麦产量和土壤肥力的影响. 应用生态学报, 2019, 30(2): 573-582.

HE H X, LI X H, BAO M, LI C, MA X L, HE G, QIU W H, WANG Z H. Effects of cultivation patterns on wheat yield and soil fertility in dryland., 2019, 30(2): 573-582. (in Chinese)

[18] MOSAVI S B, JAFARZADEH A A, NISHABOURI M R, SH O, FEIZIASL V. Application of rye green manure in wheat rotation system alters soil water content and chemical characteristics under dryland condition in Maragheh., 2009, 12(2): 178-182.

[19] TAO J M, LIU X D, LIANG Y L, NIU J, XIAO Y H, GU Y B, MA L Y, MENG D L, ZHANG Y G, HUANG W K, PENG D L, YING H Q. Maize growth responses to soil microbes and soil properties after fertilization with different green manures., 2017, 101(3): 1289-1299.

[20] LONGA C M O, NICOLA L, ANTONIELLI L, MESCALCHIN E, ZANZOTTI R, TURCO E, PERTOT I. Soil microbiota respond to green manure in organic vineyards., 2017, 123(6): 1547-1560.

[21] LAN M L, GAO M. Influence of different straws returning with landfill on soil microbial community structure under dry and water farming., 2015, 36(11): 4252-4259.

[22] 鲍士旦. 土壤农化分析. 3 版. 北京: 中国农业出版社, 2000.

BAO S D.Beijing: China Agriculture Press, 2000. (in Chinese)

[23] 关松荫. 土壤酶及其研究法. 北京: 农业出版社, 1986.

GUAN S Y.Beijing: Agriculture Press, 1986. (in Chinese)

[24] EDGAR R C. Search and clustering orders of magnitude faster than BLAST., 2010, 26(19): 2460.

[25] LANGILLE M G I, ZANEVELD J, CAPORASO J G, MCDONALD D, KNIGHTS D, REYES J A. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences., 2013, 31(9): 814-821.

[26] SCHLOSS P D, GEVERS D, WESTCOTT S L. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies., 2011, 6 (12): 1-14.

[27] SCHLOSS P D, HANDELSMAN J. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness., 2005, 71(3): 1501-1506.

[28] YE X F, LIU H E, LI Z, WANG Y, WANG Y Y, WANG H F, LIU G S. Effects of green manure continuous application on soil microbial biomass and enzyme activity., 2014, 37(4): 498-508.

[29] 叶协锋, 杨超, 李正, 敬海霞. 绿肥对植烟土壤酶活性及土壤肥力的影响. 植物营养与肥料学报, 2013, 19(2): 445-454.

YIE X F, YANG C, LI Z, GING H X. Effects of green manure in corporation on soil enzyme activities and fertility in tobacco-planting soils., 2013, 19(2): 445-454. (in Chinese)

[30] ZHAO N, ZHAO H B,YU C W, CAO Q H, LI M, CAO W D, GAO Y J. Nutrient releases of leguminous green manures in rainfed lands., 2011, 17(5): 1179-1187.

[31] BAUER J, KIRSCHBAUM M U F, WEIHERMULLER L, HUISMAN J A, HERBST M, VEREECKEN H. Temperature response of wheat decomposition is more complex than the common approaches of most multi-pool models., 2008, 40(11): 2780-2786.

[32] 杨珍平, 郝教敏, 苗果园. 混作在黄土母质生土改良中的应用. 应用与环境生物学报, 2011, 17(3): 388-392.

YANG Z P, HAO J M, MIAO G Y. Utilization of mixture cropping to improve immature loess soil., 2011, 17(3): 388-392. (in Chinese)

[33] 王秀呈. 稻—稻—绿肥长期轮作对水稻土壤及根系细菌群落的影响[D]. 北京: 中国农业科学院, 2015.

WANG X C. The influence of long term rice-rice-green manure rotation on rice bulk soil and root-associated bacteria[D].Beijing: Chinese Academy of Agricultural Sciences, 2015. (in Chinese)

[34] 林叶春, 李雨, 陈伟, 陈懿, 高维常, 和凤梅, 黄化刚, 潘文杰. 绿肥压青对喀斯特地区植烟土壤细菌群落特征的影响. 中国土壤与肥料, 2018(3): 161-167.

LIN Y C, LI Y, CHEN W, CHEN Y, GAO W C, HE F M, HUANG H G, PAN W J.Effects of green manures on the bacterial community characteristics of the rhizosphere soil in flue-cured tobacco., 2018(3): 161-167. (in Chinese)

[35] 包明, 何红霞, 马小龙, 王朝辉, 邱炜红. 化学氮肥与绿肥对麦田土壤细菌多样性和功能的影响. 土壤学报, 2018, 55(3): 734-743.

BAO M, HE H X, MA X L, WANG Z H, QIU W H. Effects of chemical nitrogen fertilizer and green manure on diversity and functions of soil bacteria in wheat field.2018, 55(3): 734-743. (in Chinese)

[36] ELLA W, HALLIN S, PHILIPPOT L. Differential responses of bacterial and archaeal groups at high taxonomical ranks to soil management., 2010, 42(10): 1759-1765.

[37] FIERER N, LAUBER C L, RAMIREZ K S, ZANEVELD J, BRADFORD M A, KNIGHT R. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients., 2012, 6(5): 1007-1017.

[38] 王伏伟, 王晓波, 李金才, 叶爱华, 王妍, 车威, 朱林. 施肥及秸秆还田对砂姜黑土细菌群落的影响. 中国生态农业学报, 2015, 23(10): 1302-1311.

WANG F W, WANG X B, LI J C, YE A H, WANG Y, CHE W, ZHU L. Effects of fertilization and straw incorporation on bacterial communities in lime concretion black soil., 2015, 23(10): 1302-1311. (in Chinese)

[39] 魏志文, 李韵雅, 江威, 廖祥儒. 无锡地区常见树木根际土壤酸杆菌多样. 生态学杂志, 2018, 37(9): 2649-2656.

WEI Z W, LI Y Y, JIANG W, LIAO X R. Diversity of acidobacteria in rhizosphere soils of common trees in Wuxi., 2018, 37(9): 2649-2656. (in Chinese)

[40] SAIT M, DAVIS K E R, JANSSEN P H. Effect of pH on isolation and distribution of members of subdivision 1 of the phylum acidobacteria occurring in soil., 2006, 53(2): 215-224.

[41] 曹彦强, 闫小娟, 罗红燕, 贾仲君, 蒋先军. 不同酸碱性紫色土的硝化活性及微生物群落组成.土壤学报, 2018, 55(1): 194-202.

CAO Y Q, YAN X J, LUO H Y, JIA Z J, JIANG X J. Nitrification activity and microbial community structure in purple soils with different pH., 2018, 55(1): 194-202. (in Chinese)

[42] 钱雅丽, 梁志婷, 曹铨, 杨宪龙, 沈禹颖, 王先之. 陇东旱作果园生草对土壤细菌群落组成的影响. 生态学杂志, 2018, 37(10): 3010-3017.

QIAN Y L, LIANG Z T, CAO Q, YANG X L, SHEN Y Y, WANG X Z. Effects of grass-planting on soil bacterial community composition of apple orchard in Longdong arid region., 2018, 37(10): 3010-3017. (in Chinese)

[43] FIERER N, JACKSON R B. The diversity and biogeography of soil bacterial communities., 2006, 103(3): 626-631.

[44] 张贵云, 吕贝贝, 张丽萍, 刘珍, 范巧兰, 魏明峰, 姚众, 袁嘉玮, 柴跃进. 黄土高原旱地麦田26年免耕覆盖对土壤肥力及原核微生物群落多样性的影响. 中国生态农业学报, 2019(3): 358-368.

ZHANG G Y, LV B B, ZHANG L P, LIU Z, FAN Q L, WEI M F, YAO Z, YUAN J W, CHAI Y J.Effect of long-term no-tillage with stubble on soil fertility and diversity of prokaryotic microbiome in dryland wheat soils on the Loess Plateau, China., 2019(3): 358-368. (in Chinese)

[45] LEBRUN E S, KANG S. A comparison of computationally predicted functional metagenomes and microarray analysis for microbial P cycle genes in a unique basalt-soil forest., 2018, 7: 179.

[46] KHURANA P, GOKHALE R S, MOHANTY D. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles., 2010, 11(1): 57.

[47] 李志杰, 张凯, 翟宇佳, 周强军, 耿运琪, 孙飞. 线虫烯脂酰CoA水合酶的克隆、表达、纯化及初步晶体学分析. 生物物理学报, 2010, 26(1): 37-48.

LI Z J, ZHANG K, ZHAI Y J, ZHOU Q J, GENG Y Q, SUN F.Cloning, expression, purification and preliminary crystallographic analysis ofEnoyl-CoA hydratase., 2010, 26(1): 37-48. (in Chinese)

[48] KRIVORUCHKO A, ZHANG Y, SIEWERS V. Microbial acetyl-CoA metabolism and metabolic engineering., 2015, 28: 28-42.

[49] ANDREWS M, LEA P J, RAVEN J A, LINDSEY K. Can genetic manipulation of plant nitrogen assimilation enzymes result in increased crop yield and greater N-use efficiency an assessment., 2015, 145(1): 25-40.

[50] 吕娜娜, 沈宗专, 王东升, 刘红军, 薛超, 李荣, 沈其荣. 施用氨基酸有机肥对黄瓜产量及土壤生物学性状的影响. 南京农业大学学报, 2018, 41(03): 456-464.

LV N N, SHEN Z Z, WANG D S, LIU H J, XUE C, LI R, SHEN Q R.Effects of amino acid organic fertilizer on cucumber yield and soil biological characters.2018, 41(3): 456-464. (in Chinese)

[51] 武姣娜, 魏晓东, 李霞, 张金飞, 谢寅峰. 植物氮素利用效率的研究进展. 植物生理学报, 2018, 54(9): 1401-1408.

WU J N, WEI X D, LI X, ZHANG J F, XIE Y F. Research progress in nitrogen use efficiency in plants., 2018, 54(9): 1401-1408. (in Chinese)

[52] POLACCO J C, MAZZAFERA P, TEZOTTO T. Opinion – nickel and urease in plants: Still many knowledge gaps.2013, 199/200: 79-90.

[53] 张亚玉, 孙海, 高明, 孙长伟, 汪景宽. 吉林省人参土壤中重金属污染水平及生物有效性研究. 土壤学报, 2011, 48(6): 1306-1313.

ZHANG Y Y, SUN H, GAO M, SUN C W, WANG J K. Pollution level and bioavailability of heavy metals in ginseng soil Jilin Province., 2011, 48(6): 1306-1313. (in Chinese)

Effects of Rapeseed Green Manure on Soil Fertility and Bacterial Community in Dryland Wheat Field

LI WenGuang, YANG XiaoXiao, HUANG ChunGuo, XUE NaiWen, XIA Qing, LIU XiaoLi, ZHANG XiaoQi, YANG Si, YANG ZhenPing, GAO ZhiQiang

(College of Agriculture, Shanxi Agricultural University, Taigu 030801, Shanxi)

【】The experiment was conducted to explore the effects of rapeseed planting as green manure on soil nutrients, enzyme activities and bacterial community in dryland wheat field of Loess Plateau in order to improve farmland fertility. 【】The experiment was arranged as two-factor split-plot design including sowing rates (S1: small amount; S2: medium amount; S3: large amount) as main plots and rapeseed composting dates (D1: September 10 as early period; D2: September 20 as medium period; D3: September 30 as late period) as subplots. Soil samples at replanting maize field at local areas were considered as control to compare the soil nutrient and enzyme activities. The high-throughput sequencing and PICRUSt gene prediction analysis method was used to determine soil community composition and metabolic function.【】Rapeseed manure increased the soil nutrients, enzyme activity and soil bacterial community and soil organic matter and sucrase activity were increased the most, ranging from 11.7% to 60.5%, and 1.4% to 94.5%, respectively, compared with control. Soil nutrients and enzyme activities were significantly affected by rapeseed sowing rates and soil organic matter, and sucrase activity were significantly affected by rapeseed composting dates. Alkaline phosphatase activity and soil bacterial community were significantly influenced with interaction between sowing rate and compositing date. Soil nutrients, enzyme activities and soil bacterial community diversity in S3D3 group were the highest as compared to other treatments. For bacterial community composition, 19 bacterial populations were identified in 10 soil samples at the phylum level, of which the dominant bacterial population were Proteobacteria, Actinobacteria, Acidobacteria and Firmicutes. Cluster analysis showed that ten treatments could be clustered into four categories, including large amount, medium amount, small amount and control. Results indicated that rapeseed composting could significantly altered soil bacterial community composition. RDA analysis showed that there was a positive correlation between soil nutrients, enzyme activities and Acidobacteria, Betaproteobacteria, Gammaproteobacteria, Alphaproteobacteria, Gemmatimonadetes, Actinobacteria and bacterial diversity. According to the PICRUSt analysis, soil bacterial community had rich metabolic functions and genes encoding metabolism were highest. Amino acid metabolism, carbohydrate metabolism, energy metabolism and environmental information processes were key metabolic function pathways in secondary prediction function classification. Heat map analysis of KEGG demonstrated that rapeseed composting increased the bacterial community associated with soil carbon and nitrogen metabolism.【】Sowing amount and rapeseed composting date significantly increased the soil nutrients, enzyme activities of succeeding wheat field and effectively, and improved the composition and diversity of soil bacterial community and enhanced the beneficial bacteria. The introduction of feed rapeseed as green manure during the summer fallow period is beneficial to improve the soil fertility and provide reasonable farming system in the dryland of the Loess Plateau.

rape; returning; green manure; wheat; soil fertility; bacterial community; high-throughput sequencing; functional predicting

2019-03-15；

2019-05-29

国家公益性行业（农业）科研专项（201503120）、国家自然科学基金项目（31101113）、国家现代农业产业技术体系建设专项（CARS-03-01-24）、山西省科技创新团队项目（201605D131041）、作物生态与旱作栽培生理山西省重点实验室项目（201705D111007）

李文广，Tel：18404983422；E-mail：18404983422@163.com。

杨珍平，Tel：13903446521；E-mail：yangzp.2@163.com

（责任编辑 李云霞）