环境科学  2023, Vol. 44 Issue (2): 1095-1103   PDF    
粉垄耕作对耕地土壤细菌群落多样性及微生物网络结构的影响
夏皖豫1,2, 陈彦云1     
1. 宁夏大学生命科学学院, 银川 750021;
2. 宁夏大学生态环境学院, 银川 750021
摘要: 为了解粉垄耕作对耕地土壤细菌群落多样性和结构分布的影响,以及微生物生态分子网络结构和物种之间的相互作用关系.以宁夏引黄灌区不同耕作方式下耕地土壤为研究对象,设置粉垄耕作(DVRT)和传统耕作(CT)这2种处理,采用Illumina MiSeq技术对不同处理的耕地土壤细菌群落进行了分子生态网络分析.结果表明,0~20 cm土层中,DVRT比CT处理显著提升了总有机碳(TOC)、总氮(TN)、碱解氮(AN)、总磷(TP)、有效磷(AP)和总钾(TK)的含量;20~40 cm土层中,DVRT比CT处理显著提升了TOC、TN、TP、AP、TK和有效钾(AK)的含量.DVRT处理的含水量(WC)显著提高了8.13%~13.30%,pH和电导率(EC)分别降低了4.51%~5.85%和12.5%~13.33%.玉米不同生育期和土层中,DVRT处理的Shannon和Ace指数均增加.优势细菌门为:放线菌门、变形菌门、绿弯菌门、酸杆菌门和芽单胞菌门.细菌群落β多样性结果显示,DVRT和CT处理的土壤细菌群落结构差异明显.冗余分析和偏曼特尔检验结果表明,pH、EC和TP是影响细菌群落结构多样性的关键环境因子.分子生态网络分析表明,DVRT处理的微生物网络由更多功能相关的微生物模块组成,关键微生物的拓扑作用与CT处理不同.综上,粉垄耕作可以提升不同土层的土壤养分和WC,降低pH和EC,改善耕地土壤细菌群落多样性和微生物网络结构并增强耕地土壤中潜在的生态系统功能.
关键词: 粉垄耕作(DVRT)      分子生态网络      细菌群落结构      高通量测序      土壤理化性质     
Effect of Deep Vertical Rotary Tillage on Soil Bacterial Community Diversity and Microbial Network Structure in Cultivated Land
XIA Wan-yu1,2 , CHEN Yan-yun1     
1. School of Life Sciences, Ningxia University, Yinchuan 750021, China;
2. School of Ecology and Environment, Ningxia University, Yinchuan 750021, China
Abstract: The aim of this study was to understand the effect of deep vertical rotary tillage on the diversity and structural distribution of a soil bacterial community in cultivated land, as well as the interactions between microbial ecological molecular network and species. In this study, the cultivated soil under different tillage methods in the Yellow River diversion irrigation area in Ningxia was selected as the research object, and two treatments were set up, namely, deep vertical rotary tillage (DVRT) and conventional tillage (CT). Then, Illumina MiSeq technology was used to analyze the molecular ecological network of soil bacterial community under different tillage methods in the Ningxia Yellow River irrigation area. The results showed that DVRT significantly increased the contents of total organic carbon (TOC), total nitrogen (TN), alkali-hydrolyzed nitrogen (AN), total phosphorus (TP), available phosphorus (AP), and total potassium (TK) in the 0-20 cm soil layer compared with those in the CT treatment. Compared with those in the CT treatment, DVRT significantly increased the contents of TOC, TN, TP, AP, TK, and available potassium (AK) in the 20-40 cm soil layer. DVRT treatment significantly increased water content (WC) by 8.13%-13.30% and decreased pH and electrical conductivity (EC) by 4.51%-5.85% and 12.5%-13.33%, respectively. In different growth stages and soil layers of maize, the Shannon and Ace indices were increased in the DVRT treatments. The dominant bacterial phyla were: Actinobacteria, Proteobacteria, Chloroflexi, Acidobacteriota, and Gemmatimonadota. The results of bacterial community β diversity showed that DVRT and CT had significant differences in bacterial community structure. The redundancy analysis and Partial Mantel test results showed that pH, EC, and TP were the key environmental factors affecting the structural diversity of bacterial communities. Molecular ecological network analysis revealed that DVRT treatment microbial networks consisted of more functionally related microbial modules, and the topological roles of key microorganisms were different from those in the CT treatment. In conclusion, DVRT can enhance soil nutrients and WC in different soil layers, reduce pH and EC, improve soil bacterial community diversity and microbial network structure, and also enhance the potential ecosystem functions in cultivated soil.
Key words: deep vertical rotary tillage(DVRT)      molecular ecological network      bacterial community structure      high-throughput sequencing      soil physicochemical properties     

玉米(Zea mays L.)是禾本科的一年生草本植物, 是人类重要的粮食和饲料作物.目前, 传统耕作(CT)是我国主要的耕作方式, 宁夏引黄灌区用传统耕作种植玉米.然而, 传统耕作会导致犁盘紧凑且封闭, 阻碍深层土层中土壤水、肥料和热量的循环, 会加速土壤侵蚀和土壤退化, 并严重影响生态系统功能[1, 2].为解决上述问题, 农业研究人员研究了多种耕作方式来提高土壤肥力, 保障土壤的可持续性利用性[3~5].粉垄耕作(DVRT)是我国实施的一种新的耕作方法, 它结合了深松、深耕、轮作和垂直耕作的优点[6], 广泛应用于不同气候带的耕地.DVRT由配备的6个垂直螺旋钻头的粉垄机械将土壤垂直粉碎至预期深度, 在不干扰土壤层和生产新的硬质盘的情况下打破压实的犁盘和松散土壤层, 形成更适合作物生长的土壤环境[7, 8].因此, 选择适宜的耕作方式, 对改变耕地土壤硬实度, 打破原有犁底层提升土壤质量有着非常重要的作用.

土壤微生物在维持生物多样性、改善土壤碳固存、提供生态系统功能和在农业生态系统中实现养分循环方面发挥着至关重要的作用[9].土壤微生物群落受耕作方法、气候、施肥和pH值等多种因素影响[10~12].已有研究表明, DVRT对土壤微生物群落结构和多样性产生有益的影响[13, 14].然而, CT的应用可能使土壤肥力下降, 并导致土壤微生物多样性下降[15].目前, 有关DVRT如何影响土壤微生物群落和生态系统功能的研究, 已成为相关科研工作者研究的热点.

分子生态网络分析技术是基于随机矩阵理论来研究生态系统内作用关系的一种方法[16].网络分析经常揭示可能反映社区组织的非随机共变模式.例如, 直接相互作用或共同的生态位, 并为调查难以在微生物群落中评估的生态概念提供解决途径[17], 对分类共现网络模式的分析提供了对关键种群以及生物群落中重要模块成员的新见解[18].分子生态网络分析还可用于探索间作土壤中微生物相互作用或反应的变化[19].有研究已经使用这种方法来检测微生物群落之间的共现关系, 如Sun等[20]对西藏地区对照和放牧条件的土壤微生物网络碳循环和氮循环基因进行分析, 发现放牧条件下的网络关键基因(模块枢纽和连通者)与对照显著不同.胡志娥等[21]研究表明, 长期覆膜促进了细菌和真菌共生网络中物种间的协同作用, 降低了物种间的相互竞争.此外, 陈兆进等[22]基于高通量测序的分子生态网路分析表明, 新乡市重金属污染土壤的关键细菌分别为: 节杆菌属、大理石雕菌、NocardioidesFerruginibacterFlavitalea、硝化螺旋菌和Lysobacter等.然而, 针对宁夏引黄灌区中土壤细菌网络的互作分析还鲜有报道, 影响细菌网络结构的环境因子有待研究.

因此本研究中采用Illumina MiSeq技术和分子生态网络分析, 探讨了DVRT和CT对耕地土壤细菌多样性、群落组成结构及生态系统功能的影响, 以探究土壤环境因子与土壤细菌群落变化规律及其与细菌多样性的关系, 以期为探寻合理的耕作方式, 提升耕地质量, 改善土壤微生物功能提供科学依据.

1 材料与方法 1.1 试验样地概况

田间试验于2020年4~9月在宁夏回族自治区石嘴山市平罗县头扎村(38°55′52″N, 106°39′39″E) 进行, 属温带大陆性气候, 海拔为1 091 m, 年蒸发量为1 755 mm, 平均风速为2~3 m ·s-1, 年无霜期为171 d, 年降水量为184 mm, 年平均气温为8.21℃.试验区土壤为灌淤土. 2020年9月20日试验玉米收获并测产, 粉垄耕作处理玉米产量为9 096 kg ·hm-2, 传统耕作处理玉米产量为6 078 kg ·hm-2.土壤环境背景值:WC(含水率)为12.38%、pH为8.53、EC(电导率)为1.38 mS ·cm-1ω[TOC(总有机碳)]为20.63 g ·kg-1ω[TN(总氮)]为0.86 g ·kg-1ω[AN(碱解氮)]为38.38 mg ·kg-1ω[TP(总磷)]为0.58 g ·kg-1ω[AP(有效磷)]为34.69 mg ·kg-1ω[TK(总钾)]为20.63 g ·kg-1ω[AK(有效钾)]为42.81 mg ·kg-1.

1.2 试验设计

本研究对2种处理进行了评估, 即传统耕作(CT)和粉垄耕作(DVRT).这2种处理在1个完整的区组设计中每个重复4次, 总共8个小区.每块小区面积为67.5 m2(4.5 m×15 m), 四周设有保护行.试验中玉米品种为Dika 5.粉垄耕作于2020年4月11日进行, 所用机器为悬挂式粉垄机.传统耕作采用旋耕机进行作业.施用商业复合肥(N ∶P2O5 ∶K2O≈12% ∶25% ∶5%), 施用量为450 kg ·hm-2.试验田田间管理均一致. 2020年4月11日机械播种, 行距0.5 m, 株距0.25 m. 2020年6月22日第1次灌水, 2020年7月15日第2次灌水.

1.3 土壤样本采集

土样采集时间为苗期(T1, 2020-05-22)、大喇叭口期(T2, 2020-07-15)和成熟期(T3, 2020-09-15), 每个时间点采用五点取样法在小区0~20 cm和20~40 cm土层进行3次重复采样取平均值.将装有样本的无菌袋运回实验室, 过2 mm筛网, 保存于-80℃, 以备后续DNA提取.收获期土壤风干后过2 mm筛网, 用于土壤理化性质的测定.

1.4 DNA抽提和PCR扩增

各土层DVRT和CT一共36个土壤样品进行Illumina MiSeq测序.根据E.Z.N.A.® soil DNA kit (Omega Bio-tek, Norcross, GA, U.S.)说明书进行微生物总DNA抽提, 使用1%的琼脂糖凝胶电泳检测DNA的提取质量, 使用NanoDrop2000测定DNA浓度和纯度; 细菌引物采用338F (5′-ACTCCTAC GGGAGGCAGCAG-3′)和806R (5′-GGACTACHVGG GTWTCTAAT-3′)[23], 扩增程序如下:95℃预变性3 min, 27个循环(95℃变性30 s, 55℃退火30 s, 72℃延伸30 s), 然后72℃稳定延伸10 min, 最后在4℃进行保存(PCR仪:ABI GeneAmp® 9700型).PCR反应体系为:5×TransStart FastPfu缓冲液4 μL, 2.5 mmol ·L-1 dNTPs 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.5 Illumina MiSeq测序和测序数据处理

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

使用fastp[24]软件对原始测序序列进行质控, 使用FLASH[25]软件进行拼接.使用UPARSE软件, 根据97%[26]的相似度对序列进行OTU聚类并剔除嵌合体.利用RDP classifier[27]对每条序列进行物种分类注释, 细菌比对Silva 16S rRNA数据库, 设置比对阈值为70%.

1.6 土壤理化分析

土壤含水量(WC)采用烘干法测定; pH采用水浸提电位法测定; 电导率(EC)的测定采用5 ∶1的水土比浸提, 用DDS-308A电导仪进行测定; 总有机碳(TOC)由TOC仪器(multi N/C 2100S, Analytik Jena)测定; 总氮(TN)由凯氏定氮仪定量; 碱解氮(AN)采用碱解扩散法测定; 总磷(TP)采用酸溶-钼锑抗比色法测定; 有效磷(AP)采用碳酸氢钠浸提-钼锑抗分光光度法测定; 总钾(TK)采用NaOH碱熔-火焰光度计法测定; 有效钾(AK)采用乙酸铵浸提-火焰光度法测定[28].

1.7 数据分析

采用SPSS 18.0分析软件(SPSS Inc., 美国) 进行独立样本t检验.使用Canoco(5.0)用于冗余分析(RDA).使用MENAP (http://129.15.40.240/mena/) 构建生态网络.使用Cytoscape软件对导出数据网络进行可视化.主坐标分析(PCoA)、偏曼特尔检验(Partial Mantel test)、置换多因素方差分析(PERMANOVA), 使用R (Version 4.0.5) vegan制图表.

2 结果与分析 2.1 土壤理化性质

表 1可知, 在0~20 cm土层中, CT处理中WC、TOC、TN、TP、AN、AP和TK含量均显著低于DVRT处理(P < 0.05), 但是pH和EC均显著高于DVRT处理(P < 0.05).在20~40 cm土层, CT处理中WC、TOC、TN、TP、AP和TK含量均显著低于DVRT处理(P < 0.05), 但是AN、pH和EC均显著高于DVRT处理(P < 0.05).

表 1 不同土层土壤理化性质1) Table 1 Soil physicochemical properties of different soil layers

2.2 不同耕作处理对细菌群落的影响

分析了Shannon指数和Ace指数, 以评估在不同耕作方式下耕地土壤细菌群落的α多样性. 如表 2所示, 在T1时期, 各土层DVRT处理的Ace指数显著高于CT处理; 在T2时期中, 0~20 cm土层的Shannon指数高于CT处理, 20~40 cm土层的Shannon指数和Ace指数均显著高于CT处理; 在T3时期, 各土层DVRT处理的Shannon指数和Ace指数均显著高于CT处理.

表 2 不同生育期内土壤细菌群落多样性指数和丰富度估计量1) Table 2 Soil bacterial diversity index and community richness in different growth periods

图 1为丰度排名前10的细菌门, 所有样本的优势门平均(相对丰度>5%)为: 放线菌门(26.53%)、变形菌门(20.32%)、绿弯菌门(17.65%)、酸杆菌门(17.03%)和芽单胞菌门(5.35%), 占微生物序列的86.87%.在0~20 cm土层中, T1时期DVRT的放线菌门和酸杆菌门丰度显著高于CT处理; T2时期DVRT的放线菌门、酸杆菌门和芽单胞菌门丰度显著高于CT处理; T3时期DVRT的放线菌门丰度显著高于CT处理.在20~40 cm土层中, T2时期DVRT的放线菌门和拟杆菌门丰度显著高于CT处理, T1和T3时期DVRT和CT处理前10细菌门丰度不显著.分析不同土层和耕作方式下耕地土壤细菌群落的β多样性.细菌群落的PCoA (图 2)和置换多因素方差分析(表 3)表明, DVRT和CT处理显著影响不同土层和不同时期的土壤细菌群落分布, 并且不同土层下DVRT与CT处理细菌群落分离.总之, 在DVRT和CT处理下, 不同土层和耕作方式中土壤细菌群落的相对丰度和组成存在显著差异.

同生育期的细菌门不同小写字母表示显着差异(t-text) 图 1 不同时期和不同处理下细菌群落门水平相对丰度 Fig. 1 Relative abundance on bacterial community at phylum level in different periods and treatments

图 2 不同处理下细菌群落组成的主坐标分析(PCoA) Fig. 2 Principal co-ordinates analysis (PCoA) of bacterial community composition under different treatments

表 3 基于Bray-Curtis距离的置换多因素方差分析1) Table 3 PERMANOVA based on Bray-Curtis distance

2.3 影响土壤细菌网络结构的环境因子

为了研究土壤环境因子与细菌群落结构之间的关系, 进行了冗余分析(RDA), 揭示了上述环境因子的重要性及其对不同土层土壤细菌群落的影响(图 3).由图 3 (a)可知, 0~20 cm土层中, 环境因子对细菌群落的总解释量为91.41%, RDA1和RDA2的解释量分别为74.53%和16.88%, DVRT处理与TOC、AN和TP相关性较高, CT处理与EC和pH相关性较高.由图 3 (b)可知, 20~40 cm土层中, 环境因子对细菌群落的总解释量为98.52%, RDA1和RDA2的解释量分别为96.12%和2.40%, DVRT处理与AP、TK和TP相关性较高, CT处理与EC和pH相关性较高.使用基于Bray-Curtis距离的偏曼特尔检验来确定影响细菌群落的关键环境因子(表 4).其中, TP、pH和EC这3个指标的Pearson相关系数均在0.5≤ |r |≤1之间, 是高度线性相关的环境因子.进一步证实环境因子(TP、pH和EC)与细菌群落显著相关(P < 0.05). 因此, 本研究的TP、pH和EC是关键的环境因子, 对耕地土壤中的细菌群落有显著影响.

(a):1.Acidobacteriota, 2.Actinobacteriota, 3. Proteobacteria, 4.Nitrospirota, 5.Unclassified, 6.Myxococcota, 7.Bacteroidota, 8.Gemmatimonadota, 9.Verrucomicrobiota, 10. Firmicutes, 11. Latescibacterota, 12.Planctomycetota, 13.Desulfobacterota, 14. Chloroflexi, 15. Methylomirabilota; (b):1.Desulfobacterota, 2.Unclassified, 3.Chloroflexi, 4.Nitrospirota, 5.Verrucomicrobiota, 6.Planctomycetota, 7.Acidobacteriota, 8. Methylomirabilota, 9.Actinobacteriota, 10.Firmicutes, 11.Proteobacteria, 12.Myxococcota, 13.Gemmatimonadota, 14.Bacteroidota, 15.Latescibacterota; 红色箭头连线表示环境因子, 蓝色箭头连线表示细菌门 图 3 16S rRNA基因与土壤特征的冗余分析(RDA) Fig. 3 Redundancy analysis (RDA) of 16S rRNA genes and soil characteristics

表 4 土壤细菌门水平与环境因子偏曼特尔检验结果1) Table 4 Partial Mantel test between soil bacterial community at phylum level and environmental factors

2.4 DVRT和CT处理下土壤细菌群落的分子生态网络分析

为了揭示不同土层中DVRT和CT处理下的耕地土壤细菌群落结构的差异, 进行了从16S rRNA基因数据得出的pMENs分析(图 4), 基于RMT方法构建的4个网络的相关性阈值均为0.94, 并总结了网络的拓扑性质(表 5).基于经验生态网络的平均聚集系数和模块性, 均高于相应的随机网络, 表明DVRT和CT处理网络具有典型小世界的特征和模块化拓扑(表 5).不同土层中的DVRT处理的细菌网络节点数比CT处理的多, 物种间关系更加复杂(表 5).从pMENs中研究了DVRT和CT处理下的土壤中网络模块数.0~20 cm土层中, DVRT的可视化网络[图 4 (a)]和CT的可视化网络[图 4 (b)]处理方法分别具有30个模块和23个模块. 20~40 cm土层中, DVRT的可视化网络[图 4 (c)]和CT的可视化网络[图 4 (d)]处理方法分别具有30个模块和24个模块.本研究构建的4个网络均以正相关关系占优, 表明了各处理的细菌之间均以协同合作关系为主, 竞争关系较弱, 但不同土层DVRT网络规模大于CT, 粉垄耕作使耕地微生物之间的合作关系更强.

大于5个节点的模块在各自的网络中用不同的颜色标记; 红色的链接表示负相关, 蓝色的链接表示正相关 图 4 RMT分析的耕地土壤中细菌群落的系统发育分子生态网络(pMENs) Fig. 4 Phylogenetic molecular ecology networks (pMENs) of bacterial communities in cultivated land soils based on RMT analysis

表 5 微生物分子生态网络特征参数 Table 5 Characteristic parameters of the molecular ecological network of soil microorganisms

利用Zi-Pi图来显示基于DVRT和CT拓扑处理的土壤中的分布(图 5). 根据简化的分类, 网络中的所有节点在Zi=2.5和Pi=0.62的阈值下都分为4个象限.网络枢纽在第一象限, 均无节点.模块枢纽在第二象限, 与连接器的拥有连通性, 可以在本模块内和/或在不同模块之间连接不同的节点.外围节点在第三象限, 拥有大于97%的DVRT和CT处理的节点, 通常连接到本模块中的节点[29].连接器在第四象限, 是土壤细菌群落的关键微生物, 这些模块连接器与本模块中的许多节点高度连接.由图 5(a)可知, 0~20 cm土层中, DVRT处理的连接器为放线菌门(OTU 78).由图 5(b)可知, 20~40 cm土层中, DVRT处理的连接器为酸杆菌门(OTU 4494)和放线菌门(OTU 1314).CT处理均没有连接器.以上结果表明, DVRT处理极大地改变了本试验样地土壤中关键微生物的拓扑作用.

图 5 土壤细菌网络的拓扑角色分布 Fig. 5 Topological role distribution map of bacterial networks

3 讨论

不同的农业耕作方式可以影响土壤理化性质.李浩等[30]研究了粉垄耕作对宿根蔗和新植蔗生长过程中土壤养分含量影响, 得出粉垄耕作可以改善土壤理化性状, 提高速效氮磷含量.哈斯格日乐等[31]研究指出, 粉垄耕作相较于常规耕作降低了土壤电导率.孙美乐等[32]研究指出, 粉垄耕作改良了棉花不同时期的土壤结构, 增加了土壤含水量.以上发现与本研究结果一致, 即粉垄耕作可以降低耕地土壤pH和电导率, 使土壤养分含量得到提升.粉垄耕作可以提供给耕地疏松的土壤环境条件, 使土壤养分的循环增强, 并且土壤孔隙度的增加可以使盐分向下转移, 同时降低了pH值, 改善了耕地土壤盐渍化.然而, 在20~40 cm土层中, 粉垄耕作的碱解氮含量低于传统耕作, 可能的原因是, 粉垄耕作可以促进作物及其根系生长, 加速了作物对碱解氮的吸收与利用, 使这个土层的碱解氮含量降低.

土壤微生物受耕作方式影响[12, 33], 但关于粉垄耕作影响土壤微生态及其土壤环境的研究较少.周佳等[33]研究指出, 粉垄处理下水稻根际土壤细菌群落多样性香农指数显著高于免耕处理.本研究指出粉垄耕作可以增加不同生育期和不同土层的细菌群落结构的Shannon指数、Ace指数和优势菌门的丰度, 并且基于Bray-Curtis距离的PCoA和置换多因素方差分析结果表明, 不同生育期各处理组间散点分离, 这与覃锋燕等[13]和刘洪等[14]研究的结果相似.由此可知, 粉垄耕作改善了耕作层土壤环境, 使土壤细菌群落α多样性和β多样性也发生了改变.进一步采用RDA分析和偏曼特尔检验分析了不同耕作方式下土壤细菌群落与环境因子的关系, 筛选出了TP、EC和pH这3个高度相关的环境因子.已有研究指出, 土壤微生物群落结构和丰度受pH、EC和TP的影响[33~35].说明本研究中, 粉垄耕作通过改变土壤环境因子, 特别是pH、TP和EC, 进而改变了土壤细菌群落多样性.

农业生产中土壤生态系统功能受耕作方式影响, 不同的耕作方法导致细菌群落结构和丰度有所差异, 这可能表明不同耕作方法可能会影响土壤生态系统的功能群落[33, 36, 37].为了更好地揭示土壤细菌的相互作用, 本研究采用功能分子生态网络方法研究土壤微生物群落结构和功能, 评估了不同耕作方式和土层之间耕地土壤生态系统功能的差异.在生态系统中, 微生物以复杂的网络结构共存和相互作用[38].图 4揭示了2种耕作方式土壤细菌群落的不同结构.模块中的高度连接的微生物可能在群落内具有相似的生态特征[39].与传统耕作相比, 粉垄耕作的细菌群落组成分布和多样性更丰富, 并且具有更多功能模块和多样化的生态功能组.在分子生态网络中, 节点数、连接数和平均连通度用来指示网络的规模和复杂程度, 平均连通度越大表明网络节点之间具有更为复杂的联系[40].粉垄耕作的节点数较高, 其网络中物种间有共生关系[41], 表明不同土层中粉垄耕作的土壤细菌网络规模更大, 物种间关系更加复杂.网络中物种间的相互关系可描述为正相关和负相关, 正相关表示物种的生态位相同或具有共生关系, 负相关代表着竞争或捕食关系[42].粉垄耕作的网络正相关连线数高于传统耕作, 表明粉垄耕作加深了土壤微生物之间的合作.由此可推断, 当外界环境发生扰动时, 粉垄耕作细菌网络可以减缓扰动的传递, 使其结构保持稳定, 并且稳定的细菌网络可能改善参与土壤养分循环的细菌群落多样性, 这与胡晓婧等[16]研究结果相似.因此, 粉垄耕作比传统耕作的土壤更具组织性和多样化的生态功能组, 增强了土壤微生物对养分利用的能力, 使耕地微生态得到改善.

Zi-Pi图可以显示出不同耕作处理的拓扑作用(图 5), 也显示了土壤微生物网络中的关键节点, 这些关键节点称为通才[43].不同处理的OTU通才被分配到不同的模块枢纽和连接器, 具有连通性. 0~20 cm土层中, 粉垄耕作处理连接器的通才属于放线菌门. 20~40 cm土层中, 粉垄耕作处理连接器的通才属于放线菌门和酸杆菌门.它们可以在自己的模块内和/或在不同模块之间连接不同的节点, 表明这些通才可以在不同的耕作处理下塑造耕地土壤微生物群落.已有研究表明, 放线菌门可以降解土壤中各种不溶性有机物质以供细胞代谢所需的各种营养, 参与养分循环[44].酸杆菌门可以降解植物残体多聚物, 并参与单碳类化合物代谢[45].总之, 粉垄耕作改变了耕地中的特征性OTU和通才微生物的拓扑作用, 增强了土壤微生物网络的养分循环功能, 这可能是粉垄耕作改善耕地土壤微生物多样性和群落结构的重要因素之一.

本研究利用Illumina MiSeq测序技术结合分子生态网络, 对宁夏引黄灌区不同耕作方式的土壤细菌群落进行研究, 发现粉垄耕作改变了细菌网络结构和关键微生物, 改善了潜在的生态系统功能, 土壤电导率、pH值和全磷含量是影响耕地土壤细菌群落结构的主要环境因子.此外, 本研究仅关注了3个玉米生育期耕作土壤的细菌结构分析.未来应针对粉垄耕作对土壤养分和土壤微生态的变化进行长期定位研究, 还需要考虑对土壤团聚体、有机碳和土壤呼吸等方向进行深入研究, 以便全面了解粉垄耕作改善耕地土壤生态系统功能的机制.

4 结论

(1) 粉垄耕作能有效提高各土层的土壤养分含量, 0~20 cm土层的WC、TOC、TN、TP、AN、AP和TK含量显著提高, 20~40 cm土层的WC、TOC、TN、TP、AP、AK和TK含量显著提高.粉垄耕作能降低各土层的土壤pH和EC.

(2) 采用高通量测序技术表明, 在不同生育期中, 粉垄耕作能够提升不同土层耕地土壤细菌的α多样性指数.在0~20 cm土层中, 粉垄耕作显著增加了苗期放线菌门和酸杆菌门的丰度, 大喇叭口期的放线菌门、酸杆菌门和芽单胞菌门的丰度, 成熟期放线菌门的丰度; 在20~40 cm土层中, 显著增加了大喇叭口期的放线菌门和拟杆菌门的丰度.细菌β多样性分析发现, 土壤细菌群落分布受耕作方式的影响.筛选出EC、pH和TP是影响耕地土壤细菌群落分布的主要环境因子.

(3) 分子生态网络分析表明, 粉垄耕作改变了耕地中的特征性OTU和关键微生物的拓扑作用, 增强了土壤微生物网络的养分循环功能.

参考文献
[1] 谷思玉, 朱玉伟, 郭兴军, 等. 不同耕作方式下黑土物理性状及其对玉米苗期生长的影响[J]. 华北农学报, 2018, 33(4): 226-231.
Gu S Y, Zhu Y W, Guo X J, et al. Effects of different tillage ways on seedling growth of maize and soil physical properties in Mollisol region[J]. Acta Agriculturae Boreali-Sinica, 2018, 33(4): 226-231.
[2] Zhang Y J, Wang S L, Wang H, et al. The effects of rotating conservation tillage with conventional tillage on soil properties and grain yields in winter wheat-spring maize rotations[J]. Agricultural and Forest Meteorology, 2018, 263: 107-117. DOI:10.1016/j.agrformet.2018.08.012
[3] 陈明刚, 向剑明, 刘勇军, 等. 垂直深旋耕配施有机碳肥对植烟土壤pH和物理性状及养分的影响[J]. 作物研究, 2020, 34(4): 349-353.
Chen M G, Xiang J M, Liu Y J, et al. Effects of deep vertical rotary tillage combined with organic carbon fertilizer on pH, physical property and nutrients of tobacco-planting soil[J]. Crop Research, 2020, 34(4): 349-353.
[4] 陈娟, 马忠明, 刘莉莉, 等. 不同耕作方式对土壤有机碳、微生物量及酶活性的影响[J]. 植物营养与肥料学报, 2016, 22(3): 667-675.
Chen J, Ma Z M, Liu L L, et al. Effect of tillage system on soil organic carbon, microbial biomass and enzyme activities[J]. Plant Nutrition and Fertilizer Science, 2016, 22(3): 667-675.
[5] 蒋云峰, 屈明秋, 王月, 等. 不同耕作方式对耕层土壤性质的影响[J]. 吉林师范大学学报(自然科学版), 2016, 37(1): 144-146.
Jiang Y F, Qu M Q, Wang Y, et al. Effects of different tillage methods on the soil properties[J]. Jilin Normal University Journal (Natural Science Edition), 2016, 37(1): 144-146.
[6] Wu F, Zhai L C, Xu P, et al. Effects of deep vertical rotary tillage on the grain yield and resource use efficiency of winter wheat in the Huang-Huai-Hai plain of China[J]. Journal of Integrative Agriculture, 2021, 20(2): 593-605. DOI:10.1016/S2095-3119(20)63405-0
[7] Wei B H, Shen Z Y, Zhou J, et al. Initial exploration on effect of saline-alkali land rebuilding and utilization by Fenlong cultivation[J]. Agricultural Science & Technology, 2017, 18(12): 2396-2400.
[8] Wei B H. Fenlong cultivation-the fourth set of farming methods invented in China[J]. Agricultural Science & Technology, 2017, 18(11): 2045-2048, 2052.
[9] Hinsinger P, Bengough A G, Vetterlein D, et al. Rhizosphere: biophysics, biogeochemistry and ecological relevance[J]. Plant and Soil, 2009, 321(1-2): 117-152. DOI:10.1007/s11104-008-9885-9
[10] Essel E, Xie J H, Deng C C, et al. Bacterial and fungal diversity in rhizosphere and bulk soil under different long-term tillage and cereal/legume rotation[J]. Soil and Tillage Research, 2019, 194. DOI:10.1016/j.still.2019.104302
[11] Jat H S, Choudhary M, Datta A, et al. Temporal changes in soil microbial properties and nutrient dynamics under climate smart agriculture practices[J]. Soil and Tillage Research, 2020, 199. DOI:10.1016/j.still.2020.104595
[12] 张常仁, 杨雅丽, 程全国, 等. 不同耕作模式对东北黑土微生物群落结构和酶活性的影响[J]. 土壤与作物, 2020, 9(4): 335-347.
Zhang C R, Yang Y L, Cheng Q G, et al. Effects of different tillages on soil microbial community structure and enzyme activity in Mollisols of China[J]. Soil and Crop, 2020, 9(4): 335-347.
[13] 覃锋燕, 杨慰贤, 彭晓辉, 等. 粉垄耕作木薯根际与非根际土壤的细菌群落结构多样性差异[J]. 西南农业学报, 2022, 35(4): 729-739.
Qin F Y, Yang W X, Peng X H, et al. Difference in the diversity of bacterial community structure in rhizosphere and non-rhizosphere soil of cassava in Fenlong tillage[J]. Southwest China Journal of Agricultural Sciences, 2022, 35(4): 729-739.
[14] 刘洪, 韦本辉, 党柯柯, 等. 粉垄耕作对甘蔗土壤微生物群落的影响[J]. 热带作物学报, 2022, 43(3): 597-605.
Liu H, Wei B H, Dang K K, et al. Effect of deep vertical rotary tillage on microbial community in sugar cane soil[J]. Chinese Journal of Tropical Crops, 2022, 43(3): 597-605.
[15] Franzluebbers A J, Arshad M A. Soil microbial biomass and mineralizable carbon of water-stable aggregates[J]. Soil Science Society of America Journal, 1997, 61(4): 1090-1097. DOI:10.2136/sssaj1997.03615995006100040015x
[16] 胡晓婧, 刘俊杰, 魏丹, 等. 东北黑土区不同纬度农田土壤真菌分子生态网络比较[J]. 应用生态学报, 2018, 29(11): 3802-3810.
Hu X J, Liu J J, Wei D, et al. Comparison on fungal molecular ecological networks of agricultural soils with different latitudes in the black soil region of Northeast China[J]. Chinese Journal of Applied Ecology, 2018, 29(11): 3802-3810.
[17] Shi S J, Nuccio E E, Shi Z J, et al. The interconnected rhizosphere: high network complexity dominates rhizosphere assemblages[J]. Ecology Letters, 2016, 19(8): 926-936.
[18] Chen Z J, Zheng Y, Ding C Y, et al. Integrated metagenomics and molecular ecological network analysis of bacterial community composition during the phytoremediation of cadmium-contaminated soils by bioenergy crops[J]. Ecotoxicology and Environmental Safety, 2017, 145: 111-118.
[19] Deng Y, Jiang Y H, Yang Y F, et al. Molecular ecological network analyses[J]. BMC Bioinformatics, 2012, 13. DOI:10.1186/1471-2105-13-113
[20] Sun X, Wang S P, Lin Q Y, et al. Molecular ecological network analyses revealing the effects of livestock grazing on soil microbial community in the Tibetan grassland[J]. Microbiology China, 2015, 42(9): 1818-1831.
[21] 胡志娥, 肖谋良, 丁济娜, 等. 长期覆膜条件下农田土壤微生物群落的响应特征[J]. 环境科学, 2022, 43(10): 4745-4754.
Hu Z E, Xiao M L, Ding J N, et al. Response characteristics of soil microbial community under long-term film mulching[J]. Environmental Science, 2022, 43(10): 4745-4754.
[22] 陈兆进, 李英军, 邵洋, 等. 新乡市镉污染土壤细菌群落组成及其对镉固定效果[J]. 环境科学, 2020, 41(6): 2889-2897.
Chen Z J, Li Y J, Shao Y, et al. Bacterial community composition in cadmium-contaminated soils in Xinxiang city and its ability to reduce cadmium bioaccumulation in Pak Choi (Brassica chinensis L.)[J]. Environmental Science, 2020, 41(6): 2889-2897.
[23] 张雅洁, 李珂, 朱浩然, 等. 北海湖微生物群落结构随季节变化特征[J]. 环境科学, 2017, 38(8): 3319-3329.
Zhang Y J, Li K, Zhu H R, et al. Community structure of microorganisms and its seasonal variation in Beihai lake[J]. Environmental Science, 2017, 38(8): 3319-3329.
[24] Chen S F, Zhou Y Q, Chen Y R, et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): i884-i890.
[25] Mago Dč T, Salzberg S L. FLASH: fast length adjustment of short reads to improve genome assemblies[J]. Bioinformatics, 2011, 27(21): 2957-2963.
[26] Edgar R C. UPARSE: highly accurate OTU sequences from microbial amplicon reads[J]. Nature Methods, 2013, 10(10): 996-998.
[27] Wang Q, Garrity G M, Tiedje J M, et al. Naiïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy[J]. Applied and Environmental Microbiology, 2007, 73(16): 5261-5267.
[28] 鲁如坤. 土壤农业化学分析方法[M]. 北京: 中国农业科技出版社, 2000.
[29] Gu S S, Hu Q L, Cheng Y Q, et al. Application of organic fertilizer improves microbial community diversity and alters microbial network structure in tea (Camellia sinensis) plantation soils[J]. Soil and Tillage Research, 2019, 195. DOI:10.1016/j.still.2019.104356
[30] 李浩, 黄金玲, 李志刚, 等. 粉垄耕作提高土壤养分有效性并促进甘蔗维管组织发育和养分吸收[J]. 植物营养与肥料学报, 2021, 27(2): 204-214.
Li H, Huang J L, Li Z G, et al. Fenlong tillage increase soil nutrient availability, and benefit vascular tissue structure and nutrient absorption of sugarcane[J]. Journal of Plant Nutrition and Fertilizers, 2021, 27(2): 204-214.
[31] 哈斯格日乐, 屈忠义, 王凡. 粉垄耕作下施加脱硫石膏和生物炭对盐渍土壤水热盐的影响研究[J]. 节水灌溉, 2019(9): 19-22, 28.
Hasigerile, Qu Z Y, Wang F. Water-heat-salt effects of applying desulphurization gypsum and biochar on saline-alkali soil under smashing ridge tillage[J]. Water Saving Irrigation, 2019(9): 19-22, 28.
[32] 孙美乐, 蔺国仓, 回经涛, 等. 粉垄耕作对新疆盐碱土理化性质及棉花生长影响[J]. 中国土壤与肥料, 2020(6): 58-64.
Sun M L, Lin G C, Hui J T, et al. Effects of smash-ridging cultivation on physical and chemical properties of saline-alkali soil and cotton growth in Xinjiang[J]. Soil and Fertilizer Sciences in China, 2020(6): 58-64.
[33] 周佳, 周灵芝, 劳承英, 等. 短期不同耕作方式对水稻根际土壤细菌群落结构多样性的影响[J]. 南方农业学报, 2020, 51(10): 2401-2411.
Zhou J, Zhou L Z, Lao C Y, et al. Effects of short-term different tillage methods on the diversity of bacterial community in rice rhizosphere soils[J]. Journal of Southern Agriculture, 2020, 51(10): 2401-2411.
[34] Yin C, Fan F L, Song A L, et al. Denitrification potential under different fertilization regimes is closely coupled with changes in the denitrifying community in a black soil[J]. Applied Microbiology and Biotechnology, 2015, 99(13): 5719-5729.
[35] 郑兰香, 杨桂钦, 高礼, 等. 宁夏第三排水沟岸边带土壤微生物多样性及影响因素分析[J]. 干旱区资源与环境, 2021, 35(9): 164-170.
Zheng L X, Yang G Q, Gao L, et al. Analysis of soil microbial diversity and influencing factors in the riparian zone of the third drainage ditch in Ningxia[J]. Journal of Arid Land Resources and Environment, 2021, 35(9): 164-170.
[36] 王凯悦, 张仲富, 王好才, 等. 哈尼梯田传统耕作方式下的水稻土壤微生物多样性及功能特征分析[J]. 西南林业大学学报, 2022, 42(2): 56-66.
Wang K Y, Zhang Z F, Wang H C, et al. Analysis of microbial diversity and functional characteristics of paddy soil under traditional tillage methods in Hani terraced fields[J]. Journal of Southwest Forestry University, 2022, 42(2): 56-66.
[37] 王倩. 保护性耕作下土壤微生物群落组成及功能研究[D]. 长春: 中国科学院东北地理与农业生态研究所, 2021.
Wang Q. Soil microbial community compositions and functions under conservation tillage[D]. Changchun: Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, 2021.
[38] 陈兆进, 林立安, 李英军, 等. 镉胁迫对芒草根际细菌群落结构、共发生网络和功能的影响[J]. 环境科学, 2021, 42(8): 3997-4004.
Chen Z J, Lin L A, Li Y J, et al. Shifts in rhizosphere bacterial community structure, co-occurrence network, and function of miscanthus following cadmium exposure[J]. Environmental Science, 2021, 42(8): 3997-4004.
[39] Kladivko E J. Tillage systems and soil ecology[J]. Soil and Tillage Research, 2001, 61(1-2): 61-76.
[40] 李冰, 李玉双, 魏建兵, 等. 不同土地利用方式对土壤细菌分子生态网络的影响[J]. 环境科学, 2020, 41(3): 1456-1465.
Li B, Li Y S, Wei J B, et al. Effects of different land use typess on the molecular ecological network of soil bacteria[J]. Environmental Science, 2020, 41(3): 1456-1465.
[41] Banerjee S, Kirkby C A, Schmutter D, et al. Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil[J]. Soil Biology and Biochemistry, 2016, 97: 188-198.
[42] Layeghifard M, Hwang D M, Guttman D S. Disentangling interactions in the microbiome: a network perspective[J]. Trends in Microbiology, 2017, 25(3): 217-228.
[43] Jiang Y J, Sun B, Li H X, et al. Aggregate-related changes in network patterns of nematodes and ammonia oxidizers in an acidic soil[J]. Soil Biology and Biochemistry, 2015, 88: 101-109.
[44] 王文晓, 李小伟, 黄文广, 等. 蒙古沙冬青根际土壤细菌群落组成及多样性与生态因子相关性研究[J]. 生态学报, 2020, 40(23): 8660-8671.
Wang W X, Li X W, Huang W G, et al. Correlations between the composition and diversity of bacterial communities and ecological factors in the rhizosphere of Ammopiptanthus mongolicus[J]. Acta Ecologica Sinica, 2020, 40(23): 8660-8671.
[45] 王光华, 刘俊杰, 于镇华, 等. 土壤酸杆菌门细菌生态学研究进展[J]. 生物技术通报, 2016, 32(2): 14-20.
Wang G H, Liu J J, Yu Z H, et al. Research progress of acidobacteria ecology in soils[J]. Biotechnology Bulletin, 2016, 32(2): 14-20.