粉垄耕作对耕地土壤酶活性、微生物群落结构和功能多样性的影响

陈彦云,夏皖豫,赵 辉,曾 明

1 宁夏大学生命科学学院,银川 750021 2 宁夏大学生态环境学院,银川 750021 3 宁夏科技特派员创业指导服务中心,银川 750021

土壤微生物是生态系统的重要组成部分,也是土壤生态系统的核心,直接或间接参与调节土壤养分循环、能量流动、有机质转换、土壤肥力形成、污染物的降解及环境净化等,特别在维持生物多样性、改善土壤碳固存、生态系统功能和自然和农业生态系统中的养分循环方面发挥着重要作用[1—3]。土壤微生物组成的质与量的变化是土壤健康状况的重要敏感指示,土壤微生物多样性可以定义为生命的丰富度,代表着微生物群落的稳定性,也反映土壤生态机制及土壤胁迫对微生物群落的影响,因此对土壤微生物群落结构和功能的研究一直是土壤与植物营养学的研究热点[4—6]。研究表明,土壤微生物群落受到各种因素的影响,包括耕作方式、气候、施肥和pH,其中耕作方式是对土壤微生物的影响尤为关键[7—10]。

玉米(ZeamaysL.)是禾本科玉蜀黍属一年生草本植物,是重要的粮食作物和饲料作物,也是全世界总产量最高的农作物。耕作方式一直是提高土壤肥力和提高作物生产力的最重要的农业实践之一,大量研究表明不同耕作方式对土壤特性和作物产量的影响[11—14]。宁夏引黄灌区的玉米种植大多采用传统耕作。然而,传统耕作导致水土流失加剧,环境污染加剧,土壤退化,犁盘紧凑且封闭,阻碍了深层土层中土壤水、肥和热量的循环影响生态系统功能[15—18]。粉垄耕作是一种集深松、深耕、旋耕、垂直耕作优点于一体的深耕技术[19]。它由配备六个立式螺旋钻的新型强力机器进行,可以在不扰乱土壤层和生产新的硬质犁的情况下打碎压实的犁盘并松动土层[20]。有研究表明,粉垄耕作模式对农田土壤质量、蓄水能力、作物产量、根系活力以及土壤微生物生物量和微生物多样性等方面的有益作用更强[19,21—24]。近年来,粉垄耕作技术在农业生产方面取得了巨大成功。然而,目前尚不清楚粉垄耕作对耕地土壤微生物群落和土壤生态系统功能的影响。

本研究利用高通量 16S rRNA和ITS基因 Illumina 测序和Biolog-ECO方法,系统地研究了粉垄耕作对土壤中微生物群落结构和土壤微生物功能多样性的影响。本研究的目的是：(1)探寻不同耕作方式下耕地土壤酶活性、土壤微生物多样性和群落结构的差异；(2)分析不同耕作方式下土壤微生物群落功能多样性及功能基团；(3)粉垄耕作如何影响土壤微生物结构、酶活性与微生物功能多样性间的相关性,使玉米产量得到提升。本研究为粉垄耕作技术在宁夏引黄灌区作物栽培中的应用发展提供了科学理论依据和农业实践理论支撑。

1 材料与方法

1.1 试验样地

试验于2019年4月—10月在宁夏回族自治区石嘴山市平罗县头闸村(106°63′N、38°95′E)进行,属温带大陆性气候,海拔1091 m,日照充足、昼夜温差大、干旱少雨、年蒸发量1755 mm,平均风速2—3 m/s,无霜期171 d,年降雨量为184 mm,年均气温8.21℃。试验区土壤为黏土,土壤状况为全磷0.86 g/kg,全氮1.10 g/kg,碱解氮76.92 mg/kg,速效磷34.69 mg/kg,速效钾214.03 mg/kg,有机碳11.42 mg/kg,总盐0.67 g/kg,含水量12.38%,pH值为8.5。

1.2 试验处理与土样采集

采用大区试验,试验面积为900 m2(60 m×15 m),设置3个处理分别为传统耕作深度20 cm(CK)、粉垄耕作深度35 cm(FL1)、粉垄耕作深度50 cm(FL2)。玉米出苗后将各处理的样地平均分为4个67.5 m2(15 m×4.5 m)的小区,即为每个处理的4次重复,每个小区之间设1 m宽走道,重复小区间设走道,宽1 m,四周设置保护行,宽3 m。

供试玉米品种为“迪卡5号”,采用吴忠市伊禾农机作业服务有限公司提供的悬挂式粉垄机作业,全面深松垂直旋磨粉碎土壤,粉垄后镇压。传统耕作采用拖拉机犁耙整地。2019年5月3日整地,5月16日播种,10月18日收获。播种密度为82500株/hm2,肥料为控释肥0.12 kg/m2(N-P-K：30- 12- 5),作为底肥一次性施入,各处理田间管理均一致。土壤采集时间为10月18日,以五点取样法在0—20 cm,20—40 cm耕作层取土样,每个土壤3个平行样品。将土壤样品装入无菌袋中,一部分鲜土过2 mm筛后4℃保存,用于土壤酶活性及Biology-ECO测定,另一部分鲜土过1 mm筛-80℃保存,用于土壤微生物群落结构的测定。

1.3 土壤酶活性的测定

脲酶用靛酚蓝比色法测定,以24 h后1 g土壤NH3-N的毫克数表示；土壤碱性磷酸酶采用磷酸苯二钠比色法,以1 h后1g土壤中对硝基苯酚的毫克数表示；转化酶用3,5-二硝基水杨酸比色法,以24 h后1 g土壤葡萄糖的毫克数表示[25]。

1.4 土壤微生物功能多样性的测定

土壤微生物功能多样性用Biolog方法进行测定[26,27]：称取5.0 g土壤样品加入装有45 mL 0.85%无菌生理盐水的三角瓶中,25℃ 150 r/min振荡30 min,冰浴1 min,后静置30 min,获得土壤样品的微生物悬浮液。将土壤悬液梯度稀释为10-3g/mL,在超净工作台中用移液枪将制备好的土壤悬液接种至Biolog微平板的各孔中,每孔150 μL,盖好盖子置于25℃培养箱中暗培养7 d,每隔24 h用酶标仪读取在590 nm与750 nm波长的数值。

Shannon指数(H)用于评估物种丰富度,计算公式如下：

H=-∑(Pi×lnPi)

(1)

式中,Pi为第i孔的相对吸光值与整个平板相对吸光值总和的比率。

McIntosh指数(U)用于评估群落均一度：

(2)

式中,ni为第i孔的相对吸光值。

Simpson指数(D)用于评估常见种优势度的指数：

D=1-∑Pi2

(3)

式中,Pi为第i孔的相对吸光值与整个平板相对吸光值总和的比率。

1.5 DNA 抽提和 PCR 扩增

各土层FL1、FL2和CK一共18个土壤样品进行IlluminaMiSeq测序。根据E.Z.N.A.® soil DNA kit (Omega Bio-tek, Norcross, GA, U.S.)说明书进行微生物总DNA抽提,使用1%的琼脂糖凝胶电泳检测DNA的提取质量,使用NanoDrop2000测定DNA浓度和纯度；细菌引物采用338F (5′-ACTCCTACGGGAGGCA GCAG- 3′)和806R (5′-GGACTACHVGGGTWTCTAAT- 3′),真菌引物采用ITS1F (5′-CTTGGTCATTTAGAGGAA GTAA- 3′)和ITS2R (5′-GCTGCGTTCTTCATCGATGC- 3′),扩增程序如下：95℃预变性3 min,27个循环(95℃ 变性 30s,55℃ 退火30s, 72℃ 延伸30s),然后 72℃ 稳定延伸10min,最后在4℃进行保存(PCR仪：ABI GeneAmp® 9700型)。PCR反应体系为：5×TransStart FastPfu 缓冲液4 μL,2.5mmol·L-1dNTPs 2μL,上游引物(5μmol·L-1)0.8 μL,下游引物(5μmol·L-1)0.8 μL, TransStart FastPfu DNA聚合酶0.4 μL,模板DNA 10 ng,补足至20 μL。

1.6 Illumina Miseq 测序和测序数据处理

将同一样本的PCR产物混合后使用2%琼脂糖凝胶回收PCR产物,利用AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA)进行回收产物纯化,2%琼脂糖凝胶电泳检测,并用QuantusTMFluorometer (Promega, USA)对回收产物进行检测定量。使用NEXTFLEX Rapid DNA-Seq Kit进行建库。利用Illumina公司的Miseq PE300平台进行测序(上海美吉生物医药科技有限公司)。

使用fastp[28]软件对原始测序序列进行质控,使用FLASH[29]软件进行拼接。使用UPARSE软件,根据97%[30,31]的相似度对序列进行OTU聚类并剔除嵌合体。利用RDP classifier[32]对每条序列进行物种分类注释,细菌比对Silva 16S rRNA数据库,真菌比对UNITE ITS数据库,设置比对阈值为70%。功能注释分别基于Tax4Fun和FunGuild用于细菌和真菌。

1.7 数据分析

方差分析(ANOVA)采用SPSS 18.0分析软件(SPSS Inc.,美国)和Tukey检验,在P=0.05显着性水平上进行。结构方程模型(SEM)用于评估假设的潜在因素,使用IBM SPSS Amos 24进行分析。使用Microsoft Excel 2016,Origin 2016和R软件绘图。

2 结果Results

2.1 粉垄耕作对土壤酶活性及玉米产量的影响

表1列出了FL1、FL2和CK处理下不同土壤深度的土壤脲酶、碱性磷酸酶和转化酶活性。在不同处理下,0—20 cm深度的3种酶活性高于20—40 cm深度。0—20 cm土层中,FL1处理的土壤脲酶活性显著高于CK处理、碱性磷酸酶和转化酶活性显著高于FL2和CK处理。20—40 cm土层中,FL1和FL2处理的脲酶活性比CK处理提高了28.80%和4.45%,碱性磷酸酶活性比CK处理提高了30.95%和16.67%,FL1处理土壤转化酶比CK处理提高了25.61%,但各处理间差异均没达到显著水平(P>0.05)。FL1、FL2处理产量分别为8.58 t/hm2和 8.38 t/hm2,与对照处理相比提高了37.9%和34.7%(P<0.05)。

表1 不同土层处理间土壤酶活性Table 1 Soil enzyme activity between treatments in different soil layers

2.2 粉垄耕作对微生物群落组成的影响

本研究中,通过高通量测序技术检测所有土壤样品的微生物,获得了30730个OTUs,包括26668个细菌OTUs和4062个真菌OTUs。对各处理平均丰度大于1的所有土壤细菌和真菌OTUs进行了Flower plot分析。对于土壤细菌OTUs,在0—20 cm和20—40 cm的FL1、FL2和CK处理共有2296个相同OTUs。0—20 cm土层中FL1、FL2和CK处理特有细菌OTUs分别是81、206和72个,20—40 cm土层中FL1、FL2和CK处理特有细菌OTUs分别是135、83和194个(图1)。对于土壤真菌OTUs,在0—20 cm和20—40 cm的FL1、FL2和CK处理共有211个相同OTUs。0—20 cm土层中FL1、FL2和CK处理特有真菌OTUs分别是129、85和102个,20—40 cm土层中FL1、FL2和CK处理特有真菌OTUs分别是67、129和88个(图1)。在门水平上分析了样本中的细菌和真菌分布(图1)。不同土层各样本中的细菌在门水平排名前10的物种组成相同(图1)。不同土层各处理间优势菌门为变形杆菌、放线菌门、绿弯菌门和酸杆菌门,达78.85%—81.83%。0—20 cm土层中的FL1、FL2和CK处理的细菌门差异不显著(P>0.05)。

图1 不同处理下土壤微生物群落组成及花瓣图Fig.1 Soil microbial community composition and flower plot under different treatments细菌相对值：Proteobacteria 变形杆菌；Actinobacteria放线菌门；Chloroflexi 绿弯菌门；Acidobacteria 酸杆菌门；Gemmatimonadetes 芽单胞菌门；Firmicutes 厚壁菌门；Bacteroidetes 拟杆菌门；Planctomycetes 浮霉菌门；Nitrospirae 硝化螺旋菌门；真菌相对值：Ascomycota 子囊菌门； Basidiomycota 担子菌门； Mortierellomycota 被孢霉门； unclassified_k__Fungi 未分类_k__真菌

20—40 cm土层中,FL1处理的Rokubacteria与FL2和CK处理差异显著(P<0.05),FL2的变形杆菌门与FL1处理差异显著(P<0.05)。在真菌在不同土层各样本中,对丰度小于1%的物种归类为others。0—20 cm土层中,FL1处理的unclassified_k__Fungi与FL2处理差异显著(P<0.05)。20—40 cm土层中,FL1处理的子囊菌门与FL2处理差异显著(P<0.05)(图1)。

在这些OTUs的基础上,进一步分析了物种丰富度指数和多样性指数(表2)。总的来说,在0—20 cm土层中,FL2与CK处理的细菌Ace、Chao和香农指数差异显著(P<0.05)；FL1与CK处理的真菌Ace指数、Chao指数和辛普森指数差异显著(P<0.05),FL1与FL2处理间Ace指数、Chao指数、香农指数和辛普森指数差异显著(P<0.05)。20—40 cm土层中,FL1与FL2处理的细菌香农指数差异显著(P<0.05),FL2与CK、FL1处理的细菌辛普森指数差异显著(P<0.05)；FL1与CK处理间真菌辛普森指数差异显著(P<0.05)。

表2 不同土层处理间土壤细菌群落丰富度和群落多样性指数Table 2 Soil bacterial community richness and community diversity index between treatments in different soil layers

2.3 粉垄耕作下土壤微生物群落功能多样性

本研究中以培养96 h获得的数据对土壤微生物功能多样性指数进行分析(表3)。0—20 cm土层中,FL1处理的香农指数与FL2处理相比差异显著(P<0.05),麦金托什指数中,FL1处理与FL2、CK处理相比差异显著(P<0.05)。20—40 cm土层中,FL2处理香农指数与FL1、CK处理相比差异显著(P<0.05),辛普森指数FL2与CK处理之间差异显著(P<0.05)。

表3 不同土层处理间土壤微生物功能多样性Table 3 Soil microbial functional diversity between treatments in different soil layers

对不同土层各处理间土壤微生物群落的碳源利用情况进行分析时发现,从土壤微生物群落利用百分比来看,以利用羧酸类、聚合物类、氨基酸类以及糖类为主,胺类和双亲化合物类的利用较低(图2)。在0—20 cm土层中,FL1、FL2和CK处理利用羧酸类碳源差异显著(P<0.05)。在20—40 cm土层中,FL1与FL2处理对胺类碳源的利用差异显著(P<0.05)(图2)。

图2 不同处理下土壤微生物对6种碳源的相对利用 Fig.2 The relative utilization of soil microorganisms to 6 carbon sources under different treatments不同小写表示不同处理间差异显著(P<0.05)

主成分分析结果表明,PC1和PC2的贡献率分别为27.7%和16.6%。影响PC1的碳源主要有糖类、氨基酸类、羧酸类和聚合物类；影响PC2的碳源主要有羧酸类和糖类。0—20 cm土层中,FL1和CK处理位于PC2正端,受PC2上的碳源影响较大,相似度较高,碳源利用差异较小；FL2处理在第三象限,与FL1和CK距离过大,碳源利用差异较大(图3)。20—40 cm土层中,CK位于第四象限,且FL1、FL2和CK处理相对分离,说明各处理间碳源利用有差异(图3)。

图3 不同处理下土壤微生物群落碳源利用主成分分析 Fig.3 Principal component analysis of carbon source utilization of soil microbial communities under different treatments

2.4 粉垄耕作下土壤细菌和真菌的功能预测

使用Tax4Fun进行潜在的微生物KEGG功能预测,代谢作用和环境信息处理是主要一级代谢通路,平均占比分别为63.67%和17.29%(图4)。在二级代谢通路上总共注释了41个特征,图5显示了前15个细菌的相对功能基团丰度。在0—20 cm和20—40 cm土层中,FL1、FL2和CK处理的主要细菌功能基团的功能注释为：氨基酸代谢、碳水化合物代谢、膜运输、信号传递、萜类化合物和聚酮化合物的代谢、能量代谢(图4)。FL1、FL2和CK处理在0—20 cm土层的辅因子和维生素的代谢差异显著(P<0.05)。在20—40 cm土层中,FL1、FL2的碳水化合物代谢显著高于CK处理(P<0.05),FL1的氨基酸代谢、辅因子和维生素的代谢和膜运输显著高于CK处理(P<0.05),FL2的信号传递、能量代谢显著高于CK处理(P<0.05)。

图4 不同处理下不同土层主要土壤细菌的功能结构(二级代谢通路)Fig.4 Functional structure of main soil bacteria in different soil layers under different treatments (Pathway level 2)

使用FUNGuild对真菌进行功能Guild注释(图5)。按真菌的营养类型主要分为9种功能类型：病原-腐生-共生菌、腐生菌、病菌、腐生-共生菌、病原-腐生-共生菌、病菌-腐生菌、共生菌和病菌-共生菌。在0—20 cm土层中,FL1的病菌-腐生菌、病菌-共生菌、腐生菌和腐生-共生菌显著高于CK(P<0.05),FL2的共生菌和腐生-共生菌显著高于CK(P<0.05)。在20—40 cm土层中,FL1的病菌-腐生菌-共生菌显著高于CK(P<0.05),FL2的病菌-腐生菌和共生菌显著高于CK(P<0.05)(图5)。

图5 不同处理下不同土层主要土壤真菌的功能结构(二级代谢通路)Fig.5 Functional structure of main soil fungi in different soil layers under different treatments (Pathway level 2)

2.5 影响微生物功能多样性的生态因素及相关性分析

结构方程模型(SEM)检测控制因素对土壤微生物功能多样性有直接和间接的影响(图6)。耕作方式的直接效应能解释土壤细菌群落多样性总方差的72%(图6)。土壤细菌群落多样性能够解释土壤真菌群落多样性总方差的11%(图6)。耕作方式和土壤微生物(细菌、真菌)群落多样性能解释土壤酶活性总方差的96% (图6)。耕作方式、土壤微生物群落多样性和土壤酶活性能解释土壤微生物功能多样性总方差的69% (图6)。

图6 耕作方法、土壤酶活性、细菌和真菌多样性指数和土壤微生物功能多样性的结构方程模型(SEM)Fig.6 Structural equation model (SEM) of farming methods, soil enzyme activity, bacterial and fungal diversity index and soil microbial functional diversity箭头的粗细表示标准化路径系数的大小;R2 值代表每个内生变量解释方差的比例；GFI：拟合优度指数；RMSEA：近似的均方根误差;*P < 0.05, **P < 0.01, ***P < 0.001

如表4所示,玉米产量与土壤脲酶、转化酶、土壤微生物香农指数(H)、麦金托什指数(U)和细菌群落香农指数(H′)显著正相关(P<0.05) ; 三种土壤酶(Ure、Inv、Alp)均与土壤微生物辛普森指数(D)、细菌群落多样性指数(Ace、Chao、H′、D′)和真菌群落多样性指数(Ace、Chao)显著正相关(P<0.05)；土壤微生物香农指数(H)与细菌群落多样性指数(Ace、Chao、H′、D′)和真菌群落多样性指数(H′、D′)显著正相关(P<0.05)。

表4 土壤酶活性、微生物群落与微生物功能多样性以及产量相关性分析Table 4 Correlation analysis of soil enzyme activity, microbial community and microbial functional diversity and yield

3 讨论

3.1 粉垄耕作对土壤微生物的影响

土壤微生物群落活动和组成受到农业耕作方式影响[33—35]。粉垄耕作较常规耕作有效提高土壤微生物数量,使土壤中的氨化细菌、固氮菌、无机磷细菌和钾细菌的数量高于旋耕耕作处理和深翻耕作处理[24]。杨博[36]研究发现,粉垄耕作为深层土壤微生物的生存提供了有利条件,0—20 cm和20—40 cm土层土壤微生物总菌落数较传统耕作显著提升。本研究中,细菌的OTUs数量明显高于真菌的OTUs,表明土壤细菌占土壤微生物总数量的比例较大。本研究发现,粉垄耕作与传统耕作相比改变了细菌的多样性和群落组成(图1,表2)。

造成的原因可能是,粉垄耕作提高了耕地土壤耕层碳存储,降低了土壤紧实度,促进了土壤大团聚体的形成[37,38]。土壤团聚结构的动态变化又反馈控制土壤微生物活动和土壤有机质组分的分异作用,使土壤有机质组分-土壤团聚结构-微生物群落结构的动态变化之间存在耦合作用[39],从而改变了细菌群落结构和多样性。此外,由于细菌提供信号分子或将介质化合物转化为真菌消耗的物质促进了真菌的生长,提高了真菌的定殖率,甚至促进真菌的完整生命周期[40—42]。这些结果表明,真菌群落间接地受到了粉垄耕作方式的影响,导致真菌群落结构、物种丰富度指数和多样性指数较传统耕作有所提升(图1,表2)。粉垄耕作提升了耕地土壤的微生物群落多样性,使参与土壤养分转化的功能型微生物增加,为维系耕地土壤微生态环境的稳定发挥了重要作用。

3.2 控制因素对土壤酶活性的影响

研究发现,耕作方式和土壤酶活性显著正相关(图6)。土壤中的氮、磷转化以及碳水化合物的代谢和运输贮藏与脲酶、转化酶和碱性磷酸酶密切相关[43,44]。作物根系可以产生酶进行养分吸收[45],粉垄耕作通过物理方式改善作物土壤环境,使作物根系具有活力,向下伸展[46],释放出更多的酶来加强与土壤微生境养分的互换[24],使粉垄耕作比传统耕作具有更高的土壤酶活性。土壤中的酶不仅来源于作物根系,还来源于土壤微生物的代谢活动[47]。异养型细菌为了自身发展,会释放出相应功能的酶,来参与土壤养分的循环[48]。王丹等[49]研究发现,放线菌门(Actinobacteriota)可以降解土壤中各种不溶性有机物质以供细胞代谢所需的各种营养,与水解酶的活性有显著相关。拟杆菌门(Bacteroidota)是具有溶磷作用的富营养菌,参与养分代谢,与脲酶、转化酶和碱性磷酸酶呈正相关[50]。真菌群落的发展使养分循环加快,能量代谢有关的真菌和参与养分循环土壤酶高度相关[51],也进一步证明土壤酶活性的高低受土壤真菌的影响。传统耕作不利于作物表层根系发育和土壤微生物活动[52],导致土壤中酶活性低。粉垄耕作通过改善耕地土壤物理结构,促进了作物根系的生长和土壤微生物自身发展,并释放出相应的酶,加快了土壤养分的循环,使粉垄耕作下的土壤酶活性高于传统耕作(表1),并且脲酶和转化酶活性与产量显著正相关(表4),粉垄耕作处理的玉米产量得到提升。

3.3 控制因素对功能多样性的影响

已有研究表明,根系的生长可以释放酶和碳源[53],从而影响分解碳源的微生物,有利于土壤微生物多样性的增加[54,55]。粉垄耕作与传统耕作之间微生物群落结构存在差异,粉垄耕作深度35 cm显著提升了0—20 cm土层的功能多样性,粉垄耕作深度50 cm显著提升了20—40 cm土层微生物功能多样性。造成的原因可能是,粉垄耕作可以促进根系生长,使根系释放出更多的碳源,从而提高了微生物功能多样性。有研究表明,患病根系增加了羧酸类和胺类碳源的释放,从而提升了土壤微生物对羧酸类和胺类碳源的利用[56,57]。本研究得出,在0—20 cm土层中粉垄耕作对羧酸类碳源的利用率显著降低；在20—40 cm土层中,粉垄耕作深度35 cm和粉垄耕作深度50 cm相比显著降低了对胺类碳源的利用比例。说明粉垄耕作可以减少作物根系患病的几率,患病根系向周围分泌羧酸类和胺类碳源减弱,使微生物利用碳源的能力发生改变。

土壤微生物细菌和真菌群落多样性可以直接影响微生物功能多样性(图6)。已有研究表明,土壤中微生物群落结构和碳源代谢紧密相关[58]。土壤中拥有代谢、遗传信息处理、有机系统3类功能基因的共养生物(copiotrophs、R-strategists)和寡养生物(oligotrophs、K-strategists)相互作用,促进了细菌群落从土壤中吸收营养物质和溶解铁和小分子等来加速自身生长,使土壤细菌代谢旺盛,提高了细菌功能群落多样性[59—61]。真菌功能预测的腐生营养型为最主要的营养型,这可能与子囊菌门为优势菌门有关。子囊菌门大多为腐生菌,是土壤中重要的分解者[62],可以分解难降解的有机质,在养分循环方面起着重要作用[63]。另外,含有共生营养型的丛枝菌根真菌(AMF)与宿主植物建立共生关系后, 不仅显著影响植物生长, 还能引起根系分泌物的变化[64,65]。AMF与植物共生可以影响土壤细菌群落结构和功能多样性[66]。本研究得出,粉垄耕作下细菌功能预测的一级代谢通路代谢作用和环境信息处理比传统耕作显著增加(图4),真菌的功能预测出的营养类型也比传统耕作显著提高(图5)。原因可能是,粉垄耕作可以增加了细菌的功能基团,细菌增强了养分的摄取能力,功能代谢得到提升。真菌不同营养类型的功能基团增加,提升了自己参与碳源代谢的能力。真菌不仅可以分解出细菌利用的碳源,而且间接促进根系向土壤释放碳源[67]。细菌真菌相互作用,使土壤中代谢底物变得丰富多样,土壤微生物功能多样性提升。此外,粉垄耕作下的病菌-腐生菌、病菌-共生菌、腐生菌、腐生-共生菌和病菌-腐生菌-共生菌等功能基团不同程度下均高于传统耕作处理,可能的原因是,粉垄耕作处理显著改变了真菌群落结构,真菌数量增多,导致相应功能的功能基团也增多。

土壤酶活性作为微生物活性的指示物,与微生物功能多样性显著相关[68—71]。段益莉等[72]研究表明,土壤酶活性的高低不仅可以反映出有机物质的水解程度,还能反映出对碳源利用能力的高低。本研究也发现,土壤酶活性显著影响了微生物功能多样性(图6),表明粉垄耕作可以通过提升酶活性间接的使土壤微生物功能多样性提高,使土壤微生物群落利用养分的能力加强,群落功能结构完善,增加了土壤生态系统的稳定性,提升了玉米的产量。

4 结论

本文研究了粉垄耕作对玉米耕地土壤酶活性、微生物群落结构和功能代谢的影响,研究表明：粉垄耕作可以显著改善土壤细菌、真菌群落结构,增加参与养分代谢的功能基团,提高土壤微生物群落多样性,并增加了玉米产量；粉垄耕作提升了参与养分循环的土壤脲酶、转化酶和碱性磷酸酶的活性；粉垄耕作技术促进了土壤微生物养分代谢的能力,使土壤微生物总体的功能多样性得到提升；从土壤微生物群落多样性、功能多样性、酶活性和玉米产量可以看出,粉垄耕作深度35 cm处理优于粉垄耕作深度50 cm。总的来说,从土壤微生态的角度可以初步得出粉垄耕作技术可以促进了养分循环,使土壤耕地质量得到提升,增加了农作物的产量。本研究为粉垄耕作技术在中国西部干旱半干旱地区作物栽培的应用发展提供了科学依据。

耕作方式对耕地质量的影响是一个周期较长的过程。因此,随着长期粉垄耕作的使用,耕地物理结构和土壤微生态也会发生着变化。本团队仅研究了一年的粉垄耕作不同深度对耕地土壤酶活性、微生物群落结构和功能多样性的影响,未来会在此试验田进行长期工作,今后会对土壤微生态的变化进行长期检测和跟踪研究,并进一步增加关于土壤团聚体、土壤碳含量以及土壤病害等研究内容。