环境科学  2023, Vol. 44 Issue (2): 1104-1119   PDF    
基于宏基因组学揭示咸水滴灌对棉田土壤微生物的影响
杜思垚, 陈静, 刘佳炜, 郭晓雯, 闵伟     
石河子大学农学院, 新疆生产建设兵团绿洲生态农业重点实验室, 石河子 832003
摘要: 咸水灌溉已成为缓解干旱区淡水短缺的重要手段, 但长期咸水灌溉会造成土壤盐分积累, 影响土壤微生物群落结构, 进而影响土壤养分转化.通过宏基因组学的手段探究长期咸水滴灌对棉田土壤微生物群落结构的影响, 试验中灌溉水盐度(ECw)设2个处理:0.35 dS·m-1和8.04 dS·m-1(分别用FW和SW表示), 施氮量分别为0 kg·hm-2和360 kg·hm-2(分别用N0和N360表示).结果表明, 咸水灌溉提高土壤含水量、盐分、有机碳和全氮含量, 降低土壤pH和速效钾含量, 氮肥施用增加土壤有机碳、盐分和全氮含量, 降低土壤含水量、pH和速效钾含量.各处理土壤的优势菌门为:变形菌门、放线菌门、酸杆菌门、绿弯菌门和芽单胞菌门.咸水灌溉显著提高放线菌门、绿弯菌门、芽单胞菌门和厚壁菌门的相对丰度, 显著降低变形菌门、酸杆菌门、蓝细菌和硝化螺旋菌门的相对丰度.氮肥施用显著提高绿弯菌门和硝化螺旋菌门的相对丰度, 显著降低酸杆菌门、芽单胞菌门、浮霉菌门、蓝细菌和疣微菌门的相对丰度.LEfSe分析表明, 咸水灌溉对土壤微生物群落潜在生物的标志物数量无明显影响, 氮肥施用降低土壤微生物群落潜在生物的标志物数量.相关性网络图表明, 20个菌属间相关性程度不同, 其中有44个正相关关系和48个负相关关系, 网络图中核心物种为类诺卡式菌属、StreptomycesPyrinomonasCandidatus_Solibacter和慢生根瘤菌属.咸水灌溉提高反硝化基因nirKnirSnasBnorC的相对丰度, 降低硝化基因amoBamoCnxrA的相对丰度, 氮肥施用提高硝化基因amoAamoBamoChaonxrA的相对丰度, 降低反硝化基因narBnapAnasAnosZ的相对丰度.咸水灌溉会对土壤理化性质产生不利影响, SWC、EC1∶5和BD是影响土壤微生物群落结构和功能基因的主要驱动因子, 土壤微生物通过调节物种组成来适应土壤盐胁迫.
关键词: 宏基因组学      咸水滴灌      氮肥      土壤理化性质      微生物群落结构     
Revealing the Effect of Saline Water Drip Irrigation on Soil Microorganisms in Cotton Fields Based on Metagenomics
DU Si-yao , CHEN Jing , LIU Jia-wei , GUO Xiao-wen , MIN Wei     
Key Laboratory of Oasis Ecological Agriculture Corps, College of Agriculture, Shihezi University, Shihezi 832003, China
Abstract: Saline water irrigation has become an important means to alleviate the shortage of freshwater in arid areas. However, long-term saline water irrigation can cause soil salinity accumulation, affect soil microbial community structure, and then affect soil nutrient transformation. In this study, we used metagenomics to investigate the effects of long-term saline water drip irrigation on soil microbial community structure in a cotton field. In the experiment, the salinity of irrigation water (ECw) was set to two treatments: 0.35 dS·m-1 and 8.04 dS·m-1 (denoted as FW and SW, respectively), and the nitrogen application rates were 0 kg·hm-2and 360 kg·hm-2 (denoted as N0 and N360, respectively). The results showed that saline water irrigation increased soil water content, salinity, organic carbon, and total nitrogen content and decreased soil pH and available potassium content. Nitrogen fertilizer application increased soil organic carbon, salinity, and total nitrogen content and decreased soil water content, pH, and available potassium content. The dominant bacterial phyla in each treatment were: Proteobacteria, Actinobacteria, Acidobacteria, Chloroflexi, and Gemmatimonadetes. Saline water irrigation significantly increased the relative abundances of Actinobacteria, Chloroflexi, Gemmatimonadetes, and Firmicutes but significantly decreased the relative abundances of Proteobacteria, Acidobacteria, Cyanobacteria, and Nitrospira. Nitrogen fertilizer application significantly increased the relative abundances of Chloroflexi and Nitrospira but significantly decreased the relative abundances of Acidobacteria, Gemmatimonadetes, Planctomycetes, Cyanobacteria, and Verrucomicrobia. LEfSe analysis showed that saline water irrigation had no significant effect on the number of potential biomarkers, and nitrogen fertilizer application decreased the number of potential biomarkers in soil microbial communities. The correlation network diagram showed that the 20 genera had different degrees of correlation, including 44 positive correlations and 48 negative correlations. The core species in the network diagram were Nocardioides, Streptomyces, Pyrinomonas, Candidatus_Solibacter, and Bradyrhizobium spp. Saline water irrigation increased the relative abundances of the denitrification genes nirK, nirS, nasB, and norC and decreased the relative abundances of the nitrification genes amoB, amoC, and nxrA, whereas nitrogen fertilizer application increased the relative abundances of the nitrification genes amoA, amoB, amoC, hao, and nxrA and decreased the relative abundances of the denitrifying genes narB, napA, nasA, and nosZ. Saline water irrigation could adversely affect soil physicochemical properties; SWC, EC1∶5, and BD were the main driving factors affecting soil microbial community structure and function genes; and soil microorganisms adapted to soil salt stress by regulating species composition.
Key words: metagenomics      saline water irrigation      nitrogen fertilizer      soil physicochemical properties      microbial community structure     

干旱区淡水资源短缺是限制农业可持续发展的关键问题[1, 2].新疆属于典型的干旱区, 农业用水主要依赖于日益减少的地表水和储量丰富但具有一定含盐量的浅层地下水.因此, 使用咸水灌溉已成为缓解该地区淡水短缺的重要措施[3].干旱区石灰性土壤母质普遍含盐, 在节水滴灌条件下, 盐分并不能被淋洗出土体, 若长期应用咸水灌溉势必造成土壤盐分累积, 增加土壤盐渍化的风险[4].此外, 咸水中的盐分离子Na+和Cl-等进入土壤后会破坏土壤结构, 增加土壤盐分, 降低土壤团聚体稳定性, 导致土壤通气透水性变差, 影响养分转化, 降低土壤肥力[5~8].土壤微生物群落在维持土壤肥力和土壤生态平衡中发挥着重要作用[9].土壤环境的变化势必会造成土壤微生物活性改变[10, 11].诸多学者均一致认为盐分对土壤微生物活性、多样性和群落结构有不利影响[12~14], 但也有研究发现盐分能够促进土壤微生物量增加、提高细菌的丰富度及多样性[15, 16].盐分虽被普遍认为是抑制微生物活性的主要非生物因素之一, 但是部分耐受极端环境的微生物在盐胁迫环境下反而被底物诱导, 表现出较高活性[17].因此咸水灌溉对干旱区农田的土壤理化性质、土壤微生物群落等产生的影响值得引起高度关注, 如何调控咸水灌溉给土壤带来的不良影响成为目前的研究热点.

氮肥一般是作物生长的限制性因素, 在盐渍化土壤中合理施用氮肥可有效降低盐分对作物生长的危害.长期施用氮肥造成土壤养分差异明显, 不仅能显著增加土壤肥力、改善土壤容重、水分含量和团聚体含量等物理化学性质, 同时对土壤微生物多样性、活性和群落组成与结构等也会产生显著影响[18~20].有研究表明, 增加氮肥施用量, 会导致土壤pH值下降, 显著降低土壤微生物生物量[21].适宜的氮肥用量能够有效调控作物根际土壤菌群数量和结构, 保障作物根际土壤微生态的稳定性, 显著提高作物根际土壤细菌和放线菌的数量, 抑制真菌数量[22].也有研究表明长期施肥与不施肥土壤的微生物群落结构差异主要表现在真菌、革兰氏阴性菌和革兰氏阳性菌[23].施入土壤中的氮肥会被微生物进行分解和转化, 其中一部分作为自身营养物质, 大部分会留在土壤中供作物吸收利用, 但是长期咸水灌溉和氮肥施用对土壤微环境的影响尚未明晰.近年来, 随着分子生物学技术的发展, 宏基因组是在微生物群落研究中应用最为广泛的组学, 能够深入剖析微生物群落结构, 对微生物群落结构的研究更为全面.

因此, 本研究针对干旱区咸水灌溉的现实问题, 通过多年田间定位试验, 采用宏基因组学的研究方法, 探究和阐明长期咸水灌溉和氮肥施用对土壤微生物群落结构的影响, 揭示咸水灌溉条件下影响土壤微生物群落的内在机制, 以期为干旱区咸水资源的合理利用和农业的可持续发展提供一定的理论依据.

1 材料与方法 1.1 试验点

试验点位于石河子大学农学院试验站(44°18′N, 86°02′E), 供试土壤为石灰性冲积土(即冲积的灰色沙漠土壤), 该地区为温带干旱地区, 属大陆性气候, 年平均降水量为180~270 mm, 年平均蒸发量为1 000~1 600 mm, 年平均气温为6.5~7.2℃.土壤表层(0~20 cm)基础理化性状为:容重1.33 g·cm-3、pH值(1∶2.5) 7.48、电导率(EC1∶5) 0.13 dS·m-1ω(有机质) 16.8 g·kg-1ω(全氮) 1.08 g·kg-1ω(速效磷) 25.9 mg·kg-1ω(速效钾) 253 mg·kg-1.

1.2 试验设计

本试验于2009年开始, 采用两因素随机区组设计, 其中灌溉水盐度(ECw)设2个处理:0.35 dS·m-1和8.04 dS·m-1(分别用FW和SW表示), 在每个灌溉水盐度下, 施氮量设2个水平:0 kg·hm-2和360 kg·hm-2(分别用N0和N360表示).其中淡水来源于当地深层地下水, 咸水通过在淡水中加入NaCl和CaCl2(质量比1∶1)获得.每个处理重复3次, 共12个试验小区, 每个小区面积为25 m2.

1.3 农艺措施

棉花于每年4月中旬播种, 9月中旬采收.棉花种植采用覆膜栽培, 每膜3管6行, 行距配置为(60+10)cm, 播种密度22.2万株·hm-2.采用干播湿出法, 2020年于4月15日播种, 为保证棉花出苗, 播种后每个处理滴出苗水30 mm.整个棉花生长期间共灌水9次, 于6月中旬开始至8月下旬结束, 灌水周期7~10 d, 灌溉定额450 mm.试验中氮肥使用尿素, 在棉花生长期间分6次随水滴施.磷肥(P2O5)和钾肥(K2O)为基肥, 用量分别为105 kg·hm-2和75 kg·hm-2, 在播种前期一次性施用, 其他管理措施参照当地大田生产.

1.4 土壤样品的采集与测定 1.4.1 样品采集

土壤样品于2020年7月25日采集, 每个处理随机选取3个采样点, 采集深度为0~20 cm.将采集的土样混合均匀并去除其中的杂质、细根等.一部分土样风干粉碎后用于土壤养分的测定, 另一部分土样储存在-80℃下用于宏基因组测序分析.

1.4.2 样品测定

理化性质测定:土壤容重采用环刀法测定; 土壤含水量采用烘干法测定; 土壤盐度采用电导率仪(DDS-308A)测定; 土壤pH值采用pH计(pHS-3C)测定; 土壤有机碳采用TOC仪(multi N/C 2100S, Analytic jena)测定; 土壤全氮采用凯式定氮仪测定; 土壤速效磷采用碳酸氢钠浸提法测定; 土壤速效钾采用火焰光度计法测定.

DNA提取和宏基因组测序:使用DNeasy PowerSoil Pro Kit试剂盒(美国)按照操作说明称取5.0 g新鲜土壤提取总微生物基因组DNA样本, 并储存在-20℃为进一步评估使用.分别使用NanoDrop ND2000分光光度计(美国)和biowest agArose琼脂糖凝胶电泳(西班牙)测量提取DNA的数量和质量.使用NanoDrop2000检测DNA纯度, 使用TBS-380检测DNA浓度.使用超声波破碎仪Covaris M220将DNA打断为约400 bp的片段, 建库试剂盒使用NEXTFLEX Rapid DNA-Seq Kit(美国).通过美吉生物技术有限公司(中国上海)的PE150策略, 使用Illumina HiSeq X-ten平台(美国Illumina)对每个文库进行测序.

原始序列读取经过质量控制, 获得质量过滤后的读数, 以便进一步分析.使用软件fastp(https://github.com/OpenGene/fastp, version 0.20.0)对reads 3′端和5′端的adapter序列进行质量剪切; 使用软件fastp去除剪切后长度小于50 bp、平均碱基质量值低于20以及含N碱基的reads, 保留高质量的pair-end reads和single-end reads; 使用基于succinct de Bruijn graphs原理的拼接软件MEGAHIT(https://github.com/voutcn/megahit, version 1.1.2)对优化序列进行拼接组装.在拼接结果中筛选≥300 bp的contigs作为最终的组装结果.使用MetaGene (http://metagene.cb.k.u-tokyo.ac.jp/)对拼接结果中的contigs进行ORF预测.选择核酸长度≥100 bp的基因, 并将其翻译为氨基酸序列.构建非冗余基因集用CD-HIT(http://www.bioinformatics.org/cd-hit/, version 4.6.1)对所有样品预测出来的基因序列进行聚类(参数为: 90% identity, 90% coverage), 每类取最长的基因作为代表序列.

1.5 数据分析

采用SPSS统计软件包(version SPSS.21.0)对土壤理化性质进行单因素方差分析, Tukey法用于不同处理间显著性的多重比较(P < 0.05).利用MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/)软件制作热图, 对不同样本进行聚类, 并确定微生物群落组成和氮转化功能关键基因.利用软件(https://cloudtutu.com/)绘制相关性热图, 通过计算两种因子之间的相关性(Spearman系数), 将获得的数值矩阵通过热图直观展示, 通过颜色变化反映二维矩阵或表格中的数据信息, 颜色深浅表示数据值的大小.利用分析软件LEfSe (http://huttenhower.sph.harvard.edu/galaxy/root?tool_id=lefse_upload) 根据分类学组成对样本按照不同的分组条件进行线性判别分析(LDA), 找出对样本划分产生显著性差异影响的物种/功能.利用软件Python中的Networkx, 绘制物种分布网络分析图, 用来展示样本与物种之间的分布情况, 通过对不同样本间的物种/功能丰度信息进行相关性分析, 可以获得物种/功能在环境样本中的共存在关系, 突显出样本之间的相似性和差异性.采用R语言进行RDA统计分析和作图.

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

灌溉水盐度、施氮量及两者交互作用均显著影响土壤理化性质(表 1).总体上, 咸水灌溉增加土壤SWC、EC1∶5、SOC和TN, 降低土壤pH和AK; 氮肥施用增加土壤SOC、EC1∶5和TN, 降低土壤SWC、pH和AK.交互作用的影响表现为:在淡水处理下, 施氮显著增加EC1∶5, 而在咸水处理下, 施氮显著降低EC1∶5; 在淡水处理下, 施氮显著降低pH, 而在咸水处理下, 施氮无显著影响; 无论是淡水处理还是咸水处理, 氮肥施用均降低土壤SWC、AP和AK, 但是降低的幅度不同, 在淡水处理下, FWN360处理的土壤SWC、AP和AK分别较FWN0处理降低12.63%、33.61%和8.81%, 在咸水处理下, SWN360处理的土壤SWC、AP和AK分别较SWN0处理降低2.55%, 24.67%和22.22%; SOC和TN的趋势与其相反, 氮肥施用均增加土壤SOC和TN, 但增加的幅度不同, 在淡水处理下, FWN360处理的土壤SOC和TN分别较FWN0处理增加22.73%和35.50%, 在咸水处理下, SWN360处理的土壤SOC和TN分别较SWN0处理增加7.46%和11.39%.

表 1 不同处理对土壤理化性质的影响1) Table 1 Soil physicochemical properties under different treatments

2.2 土壤微生物群落结构

对NCBI数据库中每个样本的Scaftigs序列进行BLASTN比较, 获得物种分类的注释度(表 2).各处理微生物以细菌为主, 占96.06%~97.01%, 其次为古细菌和真核生物, 分别占2.55%~3.44%和0.20%~0.23%.咸水灌溉降低细菌、真核生物和未分类生物的丰度, 增加古菌和病毒的丰度; 氮肥施用降低细菌、真核生物和未分类生物的丰度, 增加古菌和病毒的丰度.

表 2 不同处理下的土壤微生物注释度1)/% Table 2 Annotation of soil microbial community under different treatments/%

2.3 土壤微生物群落多样性

土壤微生物门和属水平的维恩图如图 1所示.在门水平上, FWN0、FWN360、SWN0和SWN360处理总OTUs个数为132, 共有OTUs数为129, 占总OTUs的97.72%, 4个处理均无特异的OTUs[图 1(a)].4个处理中变形菌门(Proteobacteria)、放线菌门(Actinobacteria)、酸杆菌门(Acidobacteria)、绿弯菌门(Chloroflexi)和芽单胞菌门(Gemmatimonadetes)是共有OTUs中的优势菌门.FWN0和SWN0处理特有OTUs数为2个, 为壶菌门(Chytridiomycota)和Candidatus_Atribacteria, FWN0和FWN360特有OTUs数为2个, 为Candidatus_Daviesbacteria和Candidatus_Woesebacteria, FWN0和SWN360特有OTUs数1个, 为球囊菌门(Glomeromycota), FWN360和SWN0处理特有OTUs数2个, 为纤维杆菌门(Fibrobacteres)和硅藻门(Bacillariophyta), SWN0、SWN360和FWN360处理特有OTUs数2个, 为Candidatus_Nomurabacteria和Candidatus_Parcubacteria.

图 1 不同处理土壤微生物群落维恩图 Fig. 1 Venn diagram for soil microbial community under different treatments

在属水平上, FWN0、FWN360、SWN0和SWN36处理总OTUs个数为3 172, 共有OTUs数为2 831, 占总OTUs的89.25%[图 1(b)].FWN0处理具有特异OTUs个数为34, 占总OTUs的1.07%, FWN360处理具有特异OTUs个数为20, 占0.63%, SWN0处理具有特异OTUs个数为17, 占0.53%, SWN360处理具有特异OTUs个数为17, 占0.53%, 说明咸水灌溉和氮肥施用均能够降低物种数量, 改变土壤微生物群落多样性.芽单胞菌属(Gemmatimonas)、Gemmatirosa、土壤红杆菌属(Solirubrobacter)、链霉菌属(Streptomyces)、PyrinomonasConexibacte是4个处理共有OTUs中的优势菌属.MalpighiaVicugnaOstreaPhytolacca等是FWN0处理特有的菌属, BlomiaCucurbitaMesocyclopsSaussurea等是FWN360处理特有的菌属, PetaloniaCyclotellaAntheraeaAcartia等是SWN0处理特有的菌属, EyprepocnemisAhnfeltiopsisChlorosarcinaSyntrichia等是SWN360处理特有的菌属.

2.4 土壤微生物门水平群落组成

门水平上的微生物群落组成情况如图 2所示.其优势菌门为:变形菌门(Proteobacteria)(28.22%, 平均相对丰度, 下同)、放线菌门(Actinobacteria)(22.99%)、酸杆菌门(Acidobacteria)(9.67%)、绿弯菌门(Chloroflexi)(7.83%)和芽单胞菌门(Gemmatimonadetes)(6.91%), 其相对丰度均大于5.00%, 平均相对丰度占总序列的75.61% (75.07%~76.49%).其次为浮霉菌门(Planctomycetes)、厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、蓝细菌(Cyanobacteria)和奇古菌门(Thaumarchaeota), 平均相对丰度均大于1.00%, 其它细菌门类的相对丰度合计为11.43% (10.96%~11.92%).

图 2 不同处理门水平前10种微生物的相对丰度 Fig. 2 Relative abundances of the top ten phyla under by different treatments

咸水灌溉与氮肥施用均显著影响土壤微生物门水平群落结构(图 3).在淡水灌溉下, 与FWN0处理相比, FWN360处理显著增加绿弯菌门(Chloroflexi)、厚壁菌门(Firmicutes)、硝化螺旋菌门(Nitrospirae)、异常球菌-栖热菌门(Deinococcus-Thermus)和Armatimonadetes的相对丰度, 显著降低变形菌门(Proteobacteria)、酸杆菌门(Acidobacteria)、疣微菌门(Verrucomicrobia)和广古菌门(Euryarchaeota)的相对丰度.在咸水灌溉下, 与SWN0处理相比, SWN360显著增加变形菌门(Proteobacteria)、拟杆菌门(Bacteroidetes)、奇古菌门(Thaumarchaeota)和硝化螺旋菌门(Nitrospirae)的相对丰度, 显著降低酸杆菌门(Acidobacteria)、芽单胞菌门(Gemmatimonadetes)、浮霉菌门(Planctomycetes)、厚壁菌门(Firmicutes)、Candidatus_Tectomicrobia、异常球菌-栖热菌门(Deinococcus-Thermus)和candidate_division_NC10的相对丰度.

图 3 不同处理门水平聚类分析热图 Fig. 3 Heat map of dominant phyla under different treatments

在不施氮肥条件下, 与FWN0处理相比, SWN0处理显著降低变形菌门(Proteobacteria)、酸杆菌门(Acidobacteria)、拟杆菌门(Bacteroidetes)、蓝细菌(Cyanobacteria)、硝化螺旋菌门(Nitrospirae)和疣微菌门(Verrucomicrobia)的相对丰度, 显著增加放线菌门(Actinobacteria)、绿弯菌门(Chloroflexi)、芽单胞菌门(Gemmatimonadetes)、厚壁菌门(Firmicutes)、Candidatus_Tectomicrobia、异常球菌-栖热菌门(Deinococcus-Thermus)和Armatimonadetes的相对丰度.在施氮肥条件下, 与FWN360处理相比, SWN360显著降低酸杆菌门(Acidobacteria)、蓝细菌(Cyanobacteria)、硝化螺旋菌门(Nitrospirae)、疣微菌门(Verrucomicrobia)和candidate_division_NC10的相对丰度, 显著增加放线菌门(Actinobacteria)、绿弯菌门(Chloroflexi)、厚壁菌门(Firmicutes)、异常球菌-栖热菌门(Deinococcus-Thermus)和Armatimonadetes的相对丰度.

在门分类水平下选取相对丰度前20的物种对微生物群落组成与土壤理化性质进行相关性分析(图 4).变形菌门(Proteobacteria)、酸杆菌门(Acidobacteria)、蓝细菌(Cyanobacteria)和疣微菌门(Verrucomicrobia)与土壤SOC和TN呈显著负相关, 与土壤pH呈显著正相关, 酸杆菌门(Acidobacteria)、蓝细菌(Cyanobacteria)和疣微菌门(Verrucomicrobia)与土壤AK呈显著正相关.放线菌门(Actinobacteria)和绿弯菌门(Chloroflexi)与土壤SOC和TN呈显著正相关, 绿弯菌门(Chloroflexi)与土壤AK和pH呈显著负相关.芽单胞菌门(Gemmatimonadetes)和厚壁菌门(Firmicutes)与土壤AP呈显著负相关, 与土壤TN呈显著正相关.

红色表示正相关, 蓝色表示负相关, 圆形面积大小表示相关性大小;P≥0.05无标记, *表示P < 0.05, **表示P < 0.01, ***表示P < 0.001 图 4 门水平前20种微生物与土壤理化性质的相关性分析 Fig. 4 Analysis of correlation between bacterial community of the top 20 phyla and soil physical and chemical properties

2.5 土壤微生物属水平群落组成

属水平上的微生物群落组成情况如图 5所示.通过序列比对获得各处理土壤样品中微生物群落相对丰度较高的前10个菌属, 其中10个菌属的平均相对丰度均>1.00%, 占样品总序列的26.99%(25.18%~28.09%).优势菌属有芽单胞菌属(Gemmatimonas)(3.65%, 平均相对丰度, 下同)、Gemmatirosa(3.18%)、土壤红杆菌属(Solirubrobacter)(3.02%)、链霉菌属(Streptomyces) (2.02%)、Pyrinomonas(2.00%)和Conexibacter(1.73%).

图 5 不同处理属水平前10种微生物的相对丰度 Fig. 5 Relative abundances of the top ten genus under different treatments

咸水灌溉与氮肥施用均显著影响土壤微生物属水平群落结构(图 6).在淡水灌溉下, 与FWN0处理相比, FWN360处理显著增加Streptomyces和硝化螺旋菌属(Nitrospira)的相对丰度, 显著降低土壤红杆菌属(Solirubrobacter)和红色杆菌属(Rubrobacter)的相对丰度.在咸水灌溉下, 与SWN0处理相比, SWN360处理显著增加亚硝化螺菌属(Nitrososphaera)、Sphingomonas和硝化螺旋菌属(Nitrospira)的相对丰度, 显著降低芽单胞菌属(Gemmatimonas)、Candidatus_Entotheonella和土壤红杆菌属(Solirubrobacter)的相对丰度.

图 6 不同处理属水平聚类分析热图 Fig. 6 Heat map of dominant genus under different treatments

在不施氮肥条件下, 与FWN0处理相比, SWN0处理显著降低PyrinomonasSphingomonas、硝化螺菌属(Nitrospira)、慢生根瘤菌属(Bradyrhizobium)和土壤红杆菌属(Solirubrobacter)的相对丰度, 显著增加GemmatirosaStreptomyces、类诺卡氏菌属(Nocardioides)、Candidatus_Entotheonella、红色杆菌属(Rubrobacter)、中慢生根瘤菌属(Mesorhizobium)和分支杆菌属(Mycobacterium)的相对丰度.在施氮肥条件下, 与FWN360处理相比, SWN360处理显著降低芽单胞菌属(Gemmatimonas)、Pyrinomonas、硝化螺菌属(Nitrospira)、慢生根瘤菌属(Bradyrhizobium)和土壤红杆菌属(Solirubrobacter)的相对丰度, 显著增加类诺卡氏菌属(Nocardioides)和红色杆菌属(Rubrobacter)的相对丰度.

在属分类水平下选取总丰度前20的物种对微生物群落组成与土壤理化性质进行相关性分析(图 7).芽单胞菌属(Gemmatimonas)和Gemmatirosa与土壤TN呈显著正相关, 与土壤AP呈显著负相关.土壤红杆菌属(Solirubrobacter)与土壤AP呈显著负相关.链霉菌属(Streptomyces)和类诺卡氏菌属(Nocardioides)与土壤SOC和TN呈显著正相关, 与土壤pH呈显著负相关.Pyrinomonas和硝化螺菌属(Nitrospira)与土壤AP呈显著正相关, Pyrinomonas与土壤TN呈显著负相关.

红色表示正相关, 蓝色表示负相关, 圆形面积大小表示相关性大小; P≥0.05无标记, *表示P < 0.05, **表示P < 0.01, ***表示P < 0.001 图 7 属水平前20种微生物与土壤理化性质的相关性分析 Fig. 7 Analysis of correlation between bacterial community of the top 20 genera and soil physical and chemical properties

2.6 微生物群落的LEfSe差异分析和关联分析

使用LEfSe (LDA>3.5, P < 0.05)进行组间比较分析, 分析不同处理土壤微生物群落在相对丰度上有显著差异的物种(图 8).通过LEfSe分析, 细菌共有46个显著差异种群, 其中FWN0处理15个、FWN360处理8个、SWN0处理9个、SWN360处理14个.FWN0处理的差异物种, 门水平上:酸杆菌门(Acidobacteria)、变形菌门(Proteobacteria)、疣微菌门(Verrucomicrobia)和拟杆菌门(Bacteroidetes); 纲水平上:α-变形菌纲(α-Proteobacteria)、β-变形菌纲(β-Proteobacteria)、δ-变形菌纲(δ-Proteobacteria)和鞘脂杆菌纲(Sphingobacteria); 目水平上:伯克氏菌目(Burkholderiales)和鞘脂杆菌目(Sphingobacteriales).FWN360处理的差异物种, 纲水平上:Blastocatellia和硝化螺旋菌纲(Nitrospira); 目水平上:硝化螺旋菌目(Nitrospirales); 属水平上:吡喃单胞菌属(Pyrinomonas)和硝化螺旋菌属(Nitrospira); 科水平上:硝化螺旋菌科(Nitrospiraceae).SWN0处理的差异物种, 门水平上:放线菌门(Actinobacteria)和芽单胞菌门(Gemmatimonadetes); 纲水平上:放线菌纲(Actinobacteria)和芽单胞菌纲(Gemmatimonadetes); 目水平上:芽单胞菌目(Gemmatimonadales); 科水平上:芽单胞菌科(Gemmatimonadaceae)和红螺菌科(Rhodospirillacea); 属水平上:芽单胞菌属(Gemmatimonas)和Gemmatirosa.SWN360处理的差异物种, 门水平上:绿弯菌门(Chloroflexi); 纲水平上:α-变形菌纲(α-Proteobacteria)、γ-变形菌纲(γ-Proteobacteria)和黄杆菌纲(Flavobacteriia); 目水平上:微球菌目(Micrococcales)、硝化螺旋菌目(Nitrospirales)和黄杆菌目(Flavobacteriales); 科水平上:硝化螺旋菌科(Nitrospiraceae)、黄杆菌科(Flavobacteriaceae)和类诺卡氏菌科(Nocardioidaceae).以上结果说明, 氮肥施用会降低群落差异物种数, 在未施氮肥条件下, 咸水灌溉会降低群落差异物种数, 但在施用氮肥条件下, 咸水灌溉会提高群落差异物种数.

不同颜色节点表示在相应组中显著富集并对组间差异有显著影响的微生物群; 浅黄色节点表示在不同组间没有显著差异或对组间差异没有显著影响的微生物群 图 8 土壤微生物群落的LEfSe差异分析 Fig. 8 LEfSe analysis of soil microbial communities

基于属水平的相关性网络分析发现(图 9), 网络中传递性(transitivity)为0.694 6, 网络直径Diameter为3, 节点间平均最短路径长度为1.6, 共有184个有效节点(degree), 说明网络是连通的, 表明4个处理的土壤微生物之间存在相互作用.在所有节点中, 放线菌门(Actinobacteria)、酸杆菌门(Acidobacteria)和芽单胞菌门(Gemmatimonadetes)占比最多, 分别占33.70%、21.74%和17.39%, 说明放线菌门(Actinobacteria)、酸杆菌门(Acidobacteria)和芽单胞菌门(Gemmatimonadetes)与其他物种的相关性较密切.20个菌属间相关性程度不同, 其中有44个正相关关系和48个负相关关系.基于节点中心性(Degree Centrality), 该值最高的5个属为类诺卡式菌属(Nocardioides)、StreptomycesPyrinomonasCandidatus_Solibacter和慢生根瘤菌属(Bradyrhizobium), 表明这些物种在网络中属于核心物种.

1.unclassified_p__Candidatus_Rokubacteria, 2.Nitrososphaera, 3.unclassified_p__Chloroflexi, 4.Candidatus_Solibacter, 5.unclassified_d__bacteria, 6.unclassified_c__β-Proteobacteria, 7.Solirubrobacter, 8.Nocardioides, 9.unclassified_p__Acidobacteria, 10.unclassified_c__Actinobacteria, 11.Gemmatirosa, 12.Bradyrhizobium, 13.Streptomyces, 14.Conexibacter, 15.Sphingomonas, 16.Pyrinomonas, 17.Candidatus_Entotheonella, 18.Gemmatimonas, 19. unclassified_o__Solirubrobacterales, 20. Nitrospira; 网络中的节点表示所属门类, 每个节点的大小与特定属的相对丰度成正比; 节之间的红线表示正相关, 绿线表示负相关, 线的粗细表示相关的大小 图 9 土壤微生物属水平相关性网络图分析 Fig. 9 CorrNetwork analysis of genus level communities

2.7 微生物功能代谢

基于不同处理中微生物的EggNOG注释的功能分类显示, 有50%的功能与代谢相关(图 10).在代谢类别中, 氨基酸转运和代谢(amino acid transport and metabolism)是最丰富的类别, 其次是能量生成和转化(energy production and conversion), 以及碳水化合物的转运和代谢(carbohydrate transport and metabolism).碳水化合物的转运和代谢(carbohydrate transport and metabolism)在FWN0、FWN360、SWN0和SWN360处理中分别占24%、25%、25%和25%, 信号转导机制(signal transduction mechanisms)在FWN0、FWN360、SWN0和SWN360处理中分别占26%、25%、25%和24%, 复制、重组和修复(replication, recombination, and repair)在FWN0、FWN360、SWN0和SWN360处理中分别占23%、24%、26%和27%, 转录(transcription)在FWN0、FWN360、SWN0和SWN360处理中分别占24%、25%、26%和25%.咸水灌溉和氮肥施用提高碳水化合物转运和代谢及复制、重组和修复的相对丰度, 但降低信号转导机制的相对丰度.咸水灌溉提高转录的相对丰度, 氮肥施用对转录无影响.

图 10 样本功能Circos图分析 Fig. 10 Analysis of sample function Circos diagram

咸水灌溉与氮肥施用均显著影响土壤碳氮代谢的功能基因(图 11).在淡水灌溉条件下, 与FWN0处理相比, FWN360显著提高硝化基因(amoAamoBamoChaonxrA)、反硝化基因(nirKnirSnorC)和碳代谢基因(mcrApmoA)的相对丰度, 降低反硝化(napAnarB)的相对丰度.在咸水灌溉条件下, 与SWN0处理相比, SWN360处理显著提高N代谢相关酶硝化基因(nxrAamoAamoBamoChao)和碳代谢基因(pmoA)的相对丰度, 显著降低反硝化基因(napAnasAnorCnosZ)的相对丰度.

1. nasB, 2. narB, 3. napA, 4. nasA, 5. nirK, 6. nirS, 7. hao, 8. nxrB, 9. norC, 10. nosZ, 11. amoA, 12. amoB, 13. amoC, 14. nxrA, 15. mcrA, 16. pmoA 图 11 功能基因相对丰度 Fig. 11 Relative abundance of function genes

在不施氮肥条件下, 与FWN0处理相比, SWN0处理显著提高反硝化基因(nirKnirSnasBnorC)和硝化基因(nxrA)的相对丰度, 降低硝化基因(amoBamoCnxrA)的相对丰度.在施氮肥条件下, 与FWN360处理相比, SWN360提高反硝化基因(nirKnirSnasB)的相对丰度, 降低反硝化(napAnosZ)、硝化(amoAamoC)和碳代谢(pmoAmcrA)的相对丰度.

2.8 RDA分析

微生物门水平群落与土壤理化性质间的关系如图 12所示.微生物门水平群落和土壤理化性质的RDA分析结果显示, 轴1和轴2解释群落结构总变异的85.95%和6.41%. FWN0和FWN360与SWN0和SWN360在轴1分开, 说明咸水灌溉能够显著改变土壤微生物门水平群落结构, FWN0与FWN360处理并未在轴2分开, 但SWN0和SWN360在轴2分开, 说明在咸水灌溉下氮肥施用能显著改变土壤微生物门水平群落结构.微生物门水平群落与SWC(解释度3.24%, P=0.010)、BD(解释度2.49%, P=0.004)和EC1∶5(解释度0.27%, P=0.003)、pH(解释度35.51%, P=0.010)、SOC(解释度51.94%, P=0.010)和AP(解释度14.25%, P=0.038)显著相关.

1. Acidobacteria, 2. unclassified_d__Bacteria, 3. Planctomycetes, 4.Gemmatimonadetes, 5.Actinobacteria, 6.Cyanobacteria, 7.Firmicutes, 8.Proteobacteria, 9.Bacteroidetes, 10.Chloroflexi; 箭头显示具有统计意义的环境参数, 其中箭头的方向表示相关程度最高的方向, 箭头的长度表示相关的大小, 蓝色箭头表示土壤理化性质, 红色箭头表示土壤微生物 图 12 土壤微生物门水平群落结构与土壤理化性质间RDA分析 Fig. 12 RDA analysis of the correlation among soil microbial phylum communities with soil physical and chemical properties

微生物属水平群落与土壤理化性质间的关系如图 13所示.微生物属水平群落和土壤理化性质的RDA分析结果显示, 轴1和轴2解释群落结构总变异的90.76%和4.88%, FWN0和FWN360与SWN0和SWN360在轴1分开, FWN0和SWN0与FWN360和SWN360在轴2分开, 说明咸水灌溉和氮肥施用均能显著改变土壤微生物属水平群落结构.微生物属水平群落与SWC(解释度7.52%, P=0.016)、BD(解释度6.46%, P=0.004)、EC1∶5(解释度2.83%, P=0.004)和pH(解释度30.04%, P=0.034)显著相关.

1. Pyrinomonas, 2. unclassified_p__Acidobacteria, 3. unclassified_c__β-Proteobacteria, 4. unclassified_d__Bacteria, 5. unclassified_p__Chloroflexi, 6. Streptomyces, 7. Solirubrobacter, 8. Gemmatirosa, 9. Conexibacter, 10. Gemmatimonas; 箭头显示具有统计意义的环境参数, 其中箭头的方向表示相关程度最高的方向, 箭头的长度表示相关的大小, 蓝色箭头表示土壤理化性质, 红色箭头表示土壤微生物 图 13 土壤微生物属水平群落结构与土壤理化性质间RDA分析 Fig. 13 RDA analysis of the correlation among soil microbial genus communities with soil physical and chemical properties

功能基因与土壤理化性质间的关系如图 14所示.功能基因和土壤理化性质的RDA分析结果显示, 轴1和轴2解释群落结构总变异的73.75%和15.82%, FWN0和FWN360与SWN0和SWN360在轴1分开, FWN0和SWN0与FWN360和SWN360在轴2分开, 说明咸水灌溉和氮肥施用均能显著改变微生物功能基因.功能基因与SWC(解释度41.66%, P=0.001)、BD(解释度38.25%, P=0.007)和EC1∶5(解释度66.64%, P=0.001)和AP(解释度14.25%, P=0.038)显著相关, 其中nirKnxrAnirSnasBnorCamoApmoAhaoamoBamoC与EC1∶5呈显著正相关, narBnasAnapAnosZmcrA与EC1∶5呈显著负相关.

1.nirK, 2.nxrA, 3.nirS, 4.nasB, 5.norC, 6.amoA, 7.pmoA, 8.hao, 9.amoB, 10.amoC, 11.nxrB, 12.narB, 13.nasA, 14.napA, 15.nosZ, 16.mcrA箭头显示具有统计意义的环境参数, 其中箭头的方向表示相关程度最高的方向, 箭头的长度表示相关的大小, 蓝色箭头表示土壤理化性质, 红色箭头表示功能基因 图 14 功能基因与土壤理化性质间RDA分析 Fig. 14 RDA analysis of the correlation among function genes with soil physical and chemical properties

3 讨论

淡水资源不足是全世界面临的严峻问题, 由于干旱区淡水资源匮乏, 使用咸水灌溉已成为维持干旱区农业可持续发展的必要手段.然而, 长期咸水灌溉会增加土壤盐分, 导致土壤盐渍化, 对土壤理化性质和养分转化产生不良影响.本研究发现, 咸水灌溉会增加土壤盐分、含水量和有机碳含量, 但降低土壤pH, 这是因为咸水灌溉带来的盐分会降低土壤水势, 导致作物受到盐分胁迫, 影响根系吸收水分, 从而使土壤含水量增加[24], 而土壤pH的下降, 主要是因为盐渍化土壤中强酸性离子(如NO3-、SO42-和Cl-)的积累[25].咸水灌溉会提高土壤有机碳含量, 因为持续的咸水灌溉后, 形成的盐渍化环境使土壤微生物活动减弱, 影响微生物矿化, 从而抑制有机质分解, 经过多年积累, 有机质含量增加[26].本研究还发现, 施氮会降低土壤含水量和pH, 因为氮肥会促进作物生长, 增加作物水分吸收和蒸腾作用[27], 从而提高土壤含水量, 土壤pH的降低, 原因可能是化学氮肥添加后土壤中尿素的氨化和硝化均会释放大量质子(H+), 从而降低土壤pH值[28].氮肥施用提高土壤有机碳和全氮含量, 因为施氮为微生物增殖提供了大量氮源, 也在一定程度上对土壤中微生物生物量的增加产生刺激作用, 使得土壤有机质含量进一步提高[29].本研究中氮肥施用显著降低土壤速效钾和速效磷含量, 原因是氮素对磷素有一定的协同作用, 土壤养分之间相互影响, 作物氮素的吸收会促进作物对速效磷和速效钾的利用, 从而降低土壤中速效磷和速效钾含量[30, 31].

干旱区淡水资源短缺, 咸水灌溉对土壤生态环境的影响正受到越来越多的关注.土壤微生物作为土壤养分循环的主要驱动力, 是评价土壤质量的重要指标, 在土壤养分循环中起着至关重要的作用, 包括有机质分解和土壤团聚体形成[32].盐分对土壤微生物群落活性、多样性和结构存在显著影响, 本研究发现咸水灌溉能够提高细菌和古菌的丰度, 原因可能是部分微生物对高盐环境具有一定适应性[33].咸水灌溉可以显著改变微生物群落结构, 原因可能是盐分影响土壤水分的可利用性, 导致土壤微生物活性和群落结构发生改变[34].氮肥也对土壤微生物有重要影响, 有研究发现氮肥施入土壤后很快发生固定, 被微生物同化供其吸收利用[35], 也有研究表明施氮会提高土壤微生物活性和数量[36, 37].本研究发现, 氮肥施用能够改变土壤微生物群落结构, 因为氮肥能够通过改变根系分泌物来间接影响土壤微生物群落结构[38], 施用氮肥还可以提高古菌丰度, 因为氮肥可为土壤古菌提供必需的养分[39].

盐分也会影响土壤微生物群落丰度和结构.本研究中不同处理的土壤微生物优势门类为变形菌门、放线菌门、酸杆菌门、绿弯菌门和芽单胞菌门, 这与Guo等[40]研究结果较为一致.咸水灌溉显著提高放线菌门、芽单胞菌门和厚壁菌门的丰度, 但降低硝化螺旋菌门、硝化螺菌属和疣微菌门的丰度.放线菌门是盐碱土壤中丰富的嗜盐细菌代表, 具有较好的耐盐性, 属于耐盐细菌[41], 因此放线菌对盐分有较好的适应性.厚壁菌门具有较好的耐盐性, 不易受盐分影响, 故盐分对其丰度有促进作用, 这与先前的研究结果相一致[42], 许多作物致病菌属于厚壁菌门[43], 说明咸水灌溉可能会增加土壤致病细菌, 进而对作物健康产生威胁.芽单胞菌在化学转化过程中具有很活跃的作用, 尤其是在高盐土壤中[44], 故咸水灌溉下芽单胞菌门的相对丰度高于淡水灌溉处理.硝化螺旋菌门和硝化螺菌属丰度的降低是因为硝化螺旋菌对盐分比较敏感, 在高盐环境下显著减少[45].疣微菌门是参与土壤氮循环的微生物类群[46], 咸水灌溉降低疣微菌门丰度, 说明咸水灌溉可能会对土壤氮循环产生不利影响.绿弯菌门营养型广泛, 不仅能通过光能和化能自养, 而且会利用各类元素, 推动元素的耦合循环[47], 本研究发现, 施氮能够显著提高绿弯菌门的丰度.氮肥施用还能显著提高硝化螺旋菌门、硝化螺菌属和亚硝化螺菌属的相对丰度, 这是因为施氮会提高土壤氮含量和有机质含量, 增加可供土壤微生物利用的碳源和氮源, 为细菌生物量的增加提供必需养分[39].在本研究中, 芽单胞菌属、类诺卡式菌属、红杆菌属、中慢生根瘤菌属和分支杆菌属为咸水灌溉下丰度较高的菌属, 其中芽单胞菌属归于芽单胞菌门, 可广泛分布在各种自然环境中[48], 可在恶劣环境下生存, 耐受高盐分环境[49], 中慢生根瘤菌属耐受的pH范围为4.0~10.0[50], 红杆菌属归于放线菌门, 有较强抗性, 故可适应高盐分环境.以上结果表明, 长期咸水灌溉和氮肥施用改变土壤微生物群落结构组成, 但由于对盐分和氮肥的耐受力不同, 不同微生物对盐分和氮肥的响应表现不一.

土壤微生物主要通过其功能基因的表达来调节土壤氮的转化, 盐分会导致土壤细菌群落功能发生变化[51].本研究发现, 咸水灌溉会促进反硝化基因nirKnirSnasBnorC的相对丰度, 抑制硝化基因amoBamoCnxrA的相对丰度, 这与硝化和反硝化菌群有关, 咸水中的Cl-会抑制硝化酶和硝化细菌的活性, 而耐盐反硝化菌群会在高盐度条件下进化[52].有研究表明, 盐分显著抑制反硝化基因nirKnirSnosZ的丰度[53], 这与本研究所得到的结论相反, 这是因为不同环境下的反硝化细菌种群不同, 不同反硝化细菌对盐分的耐受度不同.氮肥也会影响土壤微生物功能基因, 本研究发现, 氮肥施用会促进硝化基因amoAamoBamoChaonxrA的相对丰度, 抑制反硝化基因narBnapAnasAnosZ的相对丰度, 这是因为适量氮肥有助于增加土壤微生物种群数量, 提高硝化细菌活性, 促进硝化作用[54].

土壤微生物中的细菌对土壤物质分解、转化和循环等有较大影响[55].土壤微生物群落与土壤理化性质密切相关, 本研究通过分析土壤微生物群落与土壤理化性质间的相关关系, 发现绿弯菌门与土壤有机碳和全氮含量显著正相关, 而与土壤速效磷含量显著负相关, 绿弯菌门与光合作用有关, 能在一定程度上促进作物生长发育[47].本研究发现, 厚壁菌门与土壤全氮显著正相关, 因此土壤全氮对厚壁菌群的丰度有促进作用, 这与以往的研究结果相一致[56]; 酸杆菌门在酸性环境及寡营养条件适合生存[57], 可作为土壤贫瘠的评价指标, 有研究表明酸杆菌门与土壤pH显著负相关[58], 这与本研究结果不一致, 本研究中酸杆菌门与土壤pH显著正相关, 酸杆菌门虽具有嗜酸的特点但也能够在中性和碱性环境中大量存活[59]; 芽单胞菌属与土壤全氮显著正相关, 与土壤速效磷显著负相关, 其适宜生长在中性环境条件下, 可以增加作物抗性, 抑制病原菌侵害[60]; 放线菌门与土壤有机碳和全氮呈显著正相关, 放线菌门在有机物分解和养分循环中发挥重要作用, 这与郑佳华等[61]的研究结果相一致; 有研究表明放线菌与速效钾呈极显著正相关, 细菌与有效磷、全钾、碱解氮呈极显著正相关[62], 本研究发现酸杆菌门(Acidobacteria)、蓝细菌(Cyanobacteria)和疣微菌门(Verrucomicrobia)与土壤速效钾呈显著正相关.RDA分析表明, 土壤盐分、容重和含水量是影响土壤微生物群落结构和功能基因的主导因子, 其中盐分对土壤微生物群落结构和功能基因影响最为显著.土壤微生物对盐分变化十分敏感, 当土壤盐分含量升高时, 盐胁迫直接改变了土壤微生物的生存环境, 导致土壤微生物受到渗透胁迫, 从而降低土壤微生物种群数量和活性[63].

4 结论

(1) 长期咸水灌溉和氮肥施用显著改变土壤理化性质, 咸水灌溉提高土壤含水量、盐分、有机碳和全氮含量, 降低土壤pH和速效钾含量, 氮肥施用提高土壤有机碳和全氮含量, 降低土壤含水量、pH和速效钾含量.

(2) 咸水灌溉和氮肥施用显著改变土壤微生物群落组成, 咸水灌溉显著降低变形菌门、酸杆菌门、蓝细菌和硝化螺旋菌门的相对丰度, 显著增加放线菌门、绿弯菌门、芽单胞菌门和厚壁菌门的相对丰度.氮肥施用显著增加绿弯菌门和硝化螺旋菌门的相对丰度, 显著降低酸杆菌门、芽单胞菌门、浮霉菌门、蓝细菌、和疣微菌门的相对丰度.

(3) 咸水灌溉促进反硝化基因nirKnirSnasBnorC的相对丰度, 抑制硝化基因amoBamoCnxrA的相对丰度, 氮肥施用促进硝化基因amoAamoBamoChaonxrA的相对丰度, 抑制反硝化基因narBnapAnasAnosZ的相对丰度.

(4) 长期咸水灌溉显著改变微生物群落结构, 氮肥施用降低土壤微生物群落的潜在生物标志物的数量, 土壤盐分、容重和含水量是影响土壤微生物群落结构和功能基因的主导因子.

参考文献
[1] Ma L J, Guo H J, Min W. Nitrous oxide emission and denitrifier bacteria communities in calcareous soil as affected by drip irrigation with saline water[J]. Applied Soil Ecology, 2019, 143: 222-235. DOI:10.1016/j.apsoil.2019.08.001
[2] Rosa L, Chiarelli D D, Rulli M C, et al. Global agricultural economic water scarcity[J]. Science Advances, 2020, 6(18). DOI:10.1126/sciadv.aaz6031
[3] Hassan D F, Jafaar A A, Mohamm R J. Effect of irrigation water salinity and tillage systems on some physical soil properties[J]. Iraqi Journal of Agricultural Sciences, 2019, 50(1): 42-47.
[4] Zhou G W, Zhang W, Ma L J, et al. Effects of saline water irrigation and N application rate on NH3 volatilization and N use efficiency in a drip-irrigated cotton field[J]. Water, Air, & Soil Pollution, 2016, 227(4). DOI:10.1007/s11270-016-2806-2
[5] 冯棣, 张俊鹏, 孙池涛, 等. 长期咸水灌溉对土壤理化性质和土壤酶活性的影响[J]. 水土保持学报, 2014, 28(3): 171-176.
Feng D, Zhang J P, Sun C T, et al. Effects of long-term irrigation with saline water on soil physical-chemical properties and activities of soil enzyme[J]. Journal of Soil and Water Conservation, 2014, 28(3): 171-176.
[6] 李丹, 万书勤, 康跃虎, 等. 滨海盐碱地微咸水滴灌水盐调控对番茄生长及品质的影响[J]. 灌溉排水学报, 2020, 39(7): 39-50.
Li D, Wan S Q, Kang Y H, et al. Effects of water-salt regulation on tomato growth and quality under drip irrigation with brackish water in coastal saline-alkali soil[J]. Journal of Irrigation and Drainage, 2020, 39(7): 39-50.
[7] 刘雪艳, 丁邦新, 白云岗, 等. 微咸水膜下滴灌对土壤盐分及棉花产量的影响[J]. 干旱区研究, 2020, 37(2): 410-417.
Liu X Y, Ding B X, Bai Y G, et al. Effects of drip irrigation under a brackish water film with respect to the soil salinity and cotton yield[J]. Arid Zone Research, 2020, 37(2): 410-417.
[8] 吴雨晴, 郑春莲, 孙景生, 等. 长期咸水灌溉对棉田土壤水稳性团聚体的影响[J]. 灌溉排水学报, 2020, 39(9): 58-64, 107.
Wu Y Q, Zheng C L, Sun J S, et al. The effects of long-term saline water irrigation on stability of soil aggregates in a cotton field[J]. Journal of Irrigation and Drainage, 2020, 39(9): 58-64, 107.
[9] Jacobsen C S, Hjelmsø M H. Agricultural soils, pesticides and microbial diversity[J]. Current Opinion in Biotechnology, 2014, 27: 15-20. DOI:10.1016/j.copbio.2013.09.003
[10] Elgharably A, Marschner P. Microbial activity and biomass and N and P availability in a saline sandy loam amended with inorganic N and lupin residues[J]. European Journal of Soil Biology, 2011, 47(5): 310-315. DOI:10.1016/j.ejsobi.2011.07.005
[11] 靳正忠, 雷加强, 李生宇, 等. 咸水滴灌林地风沙土中微生物群落代谢特征[J]. 中国沙漠, 2014, 34(2): 363-370.
Jin Z Z, Lei J Q, Li S Y, et al. Characteristics of sandy soil microbial metabolisms in the forests drip-irrigated with saline water[J]. Journal of Desert Research, 2014, 34(2): 363-370.
[12] 陈绕生, 薛林宝. 园林植物土壤酶活性和微生物活性对盐碱胁迫的响应[J]. 北方园艺, 2021(4): 79-86.
Chen R S, Xue L B. Effects of saline-alkali stress on soil enzyme activities and microorganisms of garden plant[J]. Northern Horticulture, 2021(4): 79-86.
[13] 靳正忠, 雷加强, 徐新文, 等. 沙土微生物多样性与土壤肥力质量的咸水滴灌效应[J]. 生态学报, 2014, 34(13): 3720-3727.
Jin Z Z, Lei J Q, Xu X W, et al. Effect of the saline water irrigation on soil microbial diversity and fertility quality in the Tarim desert highway shelter forest land[J]. Acta Ecologica Sinica, 2014, 34(13): 3720-3727.
[14] 梁悦萍, 李科江, 张俊鹏, 等. 咸水灌溉棉田休耕期土壤胞外酶活性和微生物多样性研究[J]. 农业环境科学学报, 2018, 37(4): 732-740.
Liang Y P, Li K J, Zhang J P, et al. Extracellular enzyme activity and microbial diversity of saline-irrigated cotton field soil during the fallow period[J]. Journal of Agro-Environment Science, 2018, 37(4): 732-740.
[15] Wong V N L, Dalal R C, Greene R S B. Salinity and sodicity effects on respiration and microbial biomass of soil[J]. Biology and Fertility of Soils, 2008, 44(7): 943-953. DOI:10.1007/s00374-008-0279-1
[16] 周永学, 陈静, 李远, 等. 棉秆还田对咸水滴灌棉田土壤酶活性和细菌群落结构多样性的影响[J]. 环境科学, 2022, 43(4): 2192-2203.
Zhou Y X, Chen J, Li Y, et al. Effects of cotton stalk returning on soil enzyme activity and bacterial community structure diversity in cotton field with long-term saline water irrigation[J]. Environmental Science, 2022, 43(4): 2192-2203.
[17] 贺纪正, 张丽梅. 氨氧化微生物生态学与氮循环研究进展[J]. 生态学报, 2009, 29(1): 406-415.
He J Z, Zhang L M. Advances in ammonia-oxidizing microorganisms and global nitrogen cycle[J]. Acta Ecologica Sinica, 2009, 29(1): 406-415.
[18] 代迪, 字洪标, 杨有芳, 等. 高寒草甸土壤微生物功能多样性对氮肥添加的响应[J]. 水土保持研究, 2018, 25(2): 137-144.
Dai D, Zi H B, Yang Y F, et al. Responses of soil microbial functional diversity to nitrogen addition in an alpine meadow on northwestern plateau of Sichuan province[J]. Research of Soil and Water Conservation, 2018, 25(2): 137-144.
[19] 王伟华, 刘毅, 唐海明, 等. 长期施肥对稻田土壤微生物量、群落结构和活性的影响[J]. 环境科学, 2018, 39(1): 430-437.
Wang W H, Liu Y, Tang H M, et al. Effects of long-term fertilization regimes on microbial biomass, community structure and activity in a paddy soil[J]. Environmental Science, 2018, 39(1): 430-437.
[20] 徐润宏, 谭梅, 刘泽华, 等. 高寒湿地土壤微生物区系组成对氮添加的响应[J]. 生态科学, 2022, 41(1): 120-128.
Xu R H, Tan M, Liu Z H, et al. Response of microbial flora to nitrogen addition in alpine wetlands[J]. Ecological Science, 2022, 41(1): 120-128.
[21] 秦玮玺, 斯贵才, 雷天柱, 等. 氮肥添加对土壤微生物生物量及酶活性的影响[J]. 江苏农业科学, 2021, 49(1): 170-175.
Qin W X, Si G C, Lei T Z, et al. Influences of nitrogen fertilizer addition on soil microbial biomass and enzyme activities[J]. Jiangsu Agricultural Sciences, 2021, 49(1): 170-175.
[22] 崔亮, Kristiina L, 刘冠求, 等. 控施氮肥对马铃薯根际土壤微生物及酶活性的影响[J]. 辽宁农业科学, 2022(1): 25-29.
Cui L, Kristiina L, Liu G Q, et al. Effects of drip nitrogen control on microbial and enzyme activitiy in soil of potato[J]. Liaoning Agricultural Sciences, 2022(1): 25-29.
[23] 徐学池, 苏以荣, 王桂红, 等. 秸秆还田配施氮肥对喀斯特农田微生物群落及有机碳矿化的影响[J]. 环境科学, 2019, 40(6): 2912-2919.
Xu X C, Su Y R, Wang G H, et al. Straw returning plus nitrogen fertilizer affects the soil microbial community and organic carbon mineralization in karst farmland[J]. Environmental Science, 2019, 40(6): 2912-2919.
[24] Yuan H Z, Ge T D, Zhou P, et al. Soil microbial biomass and bacterial and fungal community structures responses to long-term fertilization in paddy soils[J]. Journal of Soils and Sediments, 2013, 13(5): 877-886.
[25] 范庆锋, 张玉龙, 陈重, 等. 保护地土壤盐分积累及其离子组成对土壤pH值的影响[J]. 干旱地区农业研究, 2009, 27(1): 16-20.
Fan Q F, Zhang Y L, Chen Z, et al. Effects of soil salinity accumulating and ion constitution on pH in the soil of protected field[J]. Agricultural Research in the Arid Areas, 2009, 27(1): 16-20.
[26] Kida M, Tomotsune M, Iimura Y, et al. High salinity leads to accumulation of soil organic carbon in mangrove soil[J]. Chemosphere, 2017, 177: 51-55.
[27] Min W, Guo H J, Zhang W, et al. Irrigation water salinity and N fertilization: effects on ammonia oxidizer abundance, enzyme activity and cotton growth in a drip irrigated cotton field[J]. Journal of Integrative Agriculture, 2016, 15(5): 1121-1131.
[28] Zhou M H, Zhu B, Brüggemann N, et al. Sustaining crop productivity while reducing environmental nitrogen losses in the subtropical wheat-maize cropping systems: a comprehensive case study of nitrogen cycling and balance[J]. Agriculture, Ecosystems & Environment, 2016, 231: 1-14.
[29] 张欣, 任海燕, 康静, 等. 增温和施氮对内蒙古荒漠草原土壤理化性质的影响[J]. 中国草地学报, 2021, 43(6): 17-24.
Zhang X, Ren H Y, Kang J, et al. Effects of warming and nitrogen addition on soil physical and chemical properties in desert steppe of Inner Mongolia[J]. Chinese Journal of Grassland, 2021, 43(6): 17-24.
[30] 冯爱青, 张民, 李成亮, 等. 控释氮肥对土壤酶活性与土壤养分利用的影响[J]. 水土保持学报, 2014, 28(3): 177-184.
Feng A Q, Zhang M, Li C L, et al. Effects of controlled release nitrogen fertilizer on soil enzyme activities and soil nutrient utilization[J]. Journal of Soil and Water Conservation, 2014, 28(3): 177-184.
[31] 胡诚, 乐群芬, 孙斌, 等. 氮肥用量与施用方法对土壤养分与小麦生长的影响[J]. 西南农业学报, 2017, 30(9): 2017-2023.
Hu C, Yue Q F, Sun B, et al. Effects of nitrogen fertilizer amount and application method on soil nutrient and wheat growth[J]. Southwest China Journal of Agricultural Sciences, 2017, 30(9): 2017-2023.
[32] Amini S, Ghadiri H, Chen C R, et al. Salt-affected soils, reclamation, carbon dynamics, and biochar: a review[J]. Journal of Soils and Sediments, 2016, 16(3): 939-953.
[33] Yang H, Hu J X, Long X H, et al. Salinity altered root distribution and increased diversity of bacterial communities in the rhizosphere soil of Jerusalem artichoke[J]. Scientific Reports, 2016, 6. DOI:10.1038/srep20687
[34] Yuan Y G, Brunel C, Van Kleunen M, et al. Salinity-induced changes in the rhizosphere microbiome improve salt tolerance of Hibiscus hamabo[J]. Plant and Soil, 2019, 443(1-2): 525-537.
[35] 陈珊, 丁咸庆, 祝贞科, 等. 秸秆还田对外源氮在土壤中转化及其微生物响应的影响[J]. 环境科学, 2017, 38(4): 1613-1621.
Chen S, Ding X Q, Zhu Z K, et al. Effect of straw application on the dynamics of exogenous nitrogen and microbial activity in paddy soil[J]. Environmental Science, 2017, 38(4): 1613-1621.
[36] 田飞飞, 纪鸿飞, 王乐云, 等. 施肥类型和水热变化对农田土壤氮素矿化及可溶性有机氮动态变化的影响[J]. 环境科学, 2018, 39(10): 4717-4726.
Tian F F, Ji H F, Wang L Y, et al. Effects of various combinations of fertilizer, soil moisture, and temperature on nitrogen mineralization and soluble organic nitrogen in agricultural soil[J]. Environmental Science, 2018, 39(10): 4717-4726.
[37] 吴春晓, 高小峰, 闫本帅, 等. 长期施肥对黄土高原梯田土壤养分特征和微生物资源限制的影响[J]. 环境科学, 2022, 43(1): 521-529.
Wu C X, Gao X F, Yan B S, et al. Effects of long-term fertilization on soil nutrient characteristics and microbial resource restrictions in a terrace on the loess plateau[J]. Environmental Science, 2022, 43(1): 521-529.
[38] 高宛莉, 刘恩科, 李来福, 等. 玉米栽培模式与施肥方式对土壤微生物和土壤肥力的影响[J]. 中国农业大学学报, 2014, 19(2): 108-117.
Gao W L, Liu E K, Li L F, et al. Effect of corn cultivation mode and fertilization on soil microbial biomass and soil fertility[J]. Journal of China Agricultural University, 2014, 19(2): 108-117.
[39] 刘春华, 吴东梅, 刘雨晖, 等. 氮沉降对米槠天然林土壤有机碳及微生物群落结构的影响[J]. 林业科学研究, 2021, 34(2): 42-49.
Liu C H, Wu D M, Liu Y H, et al. Effects of nitrogen deposition on soil organic carbon and soil microbial communities in a natural castanopsis carlesii forest[J]. Forest Research, 2021, 34(2): 42-49.
[40] Guo H N, Shi X D, Ma L J, et al. Long-term irrigation with saline water decreases soil nutrients, diversity of bacterial communities, and cotton yields in a gray desert soil in China[J]. Polish Journal of Environmental Studies, 2020, 29(6): 4077-4088.
[41] Samaei M R, Mortazavi S B, Bakhshi B, et al. Isolation, genetic identification, and degradation characteristics of n-hexadecane degrading bacteria from tropical areas in Iran[J]. Fresenius Environmental Bulletin, 2013, 22(4): 1304-1312.
[42] 赵娇, 谢慧君, 张建. 黄河三角洲盐碱土根际微环境的微生物多样性及理化性质分析[J]. 环境科学, 2020, 41(3): 1449-1455.
Zhao J, Xie H J, Zhang J. Microbial diversity and physicochemical properties of rhizosphere microenvironment in saline-alkali soils of the Yellow River Delta[J]. Environmental Science, 2020, 41(3): 1449-1455.
[43] Ehrich S, Behrens D, Lebedeva E, et al. A new obligately chemolithoautotrophic, nitrite-oxidizing bacterium, Nitrospira moscoviensis sp. nov. and its phylogenetic relationship[J]. Archives of Microbiology, 1995, 164(1): 16-23.
[44] Zhang H, Sekiguchi Y, Hanada S, et al. Gemmatimonas aurantiaca gen. nov., sp. nov., a gram-negative, aerobic, polyphosphate-accumulating micro-organism, the first cultured representative of the new bacterial phylum Gemmatimonadetes phyl. nov[J]. International Journal of Systematic and Evolutionary Microbiology, 2003, 53(4): 1155-1163.
[45] Zheng W, Xue D M, Li X Z, et al. The responses and adaptations of microbial communities to salinity in farmland soils: a molecular ecological network analysis[J]. Applied Soil Ecology, 2017, 120: 239-246.
[46] Veselá A B, Franc M, Pelantová H, et al. Hydrolysis of benzonitrile herbicides by soil actinobacteria and metabolite toxicity[J]. Biodegradation, 2010, 21(5): 761-770.
[47] 鲜文东, 张潇橦, 李文均. 绿弯菌的研究现状及展望[J]. 微生物学报, 2020, 60(9): 1801-1820.
Xian W D, Zhang X T, Li W J. Research status and prospect on bacterial phylum Chloroflexi[J]. Acta Microbiologica Sinica, 2020, 60(9): 1801-1820.
[48] DeBruyn J M, Nixon L T, Fawaz M N, et al. Global biogeography and quantitative seasonal dynamics of Gemmatimonadetes in soil[J]. Applied and Environmental Microbiology, 2011, 77(17): 6295-6300.
[49] Zeng Y H, Selyanin V, Lukeš M, et al. Characterization of the microaerophilic, bacteriochlorophyll a-containing bacterium Gemmatimonas phototrophica sp. nov., and emended descriptions of the genus Gemmatimonas and Gemmatimonas aurantiaca[J]. International Journal of Systematic and Evolutionary Microbiology, 2015, 65(8): 2410-2419.
[50] Kwon S W, Park J Y, Kim J S, et al. Phylogenetic analysis of the genera Bradyrhizobium, Mesorhizobium, Rhizobium and Sinorhizobium on the basis of 16S rRNA gene and internally transcribed spacer region sequences[J]. International Journal of Systematic and Evolutionary Microbiology, 2005, 55(1): 263-270.
[51] Chen L J, Li C S, Feng Q, et al. Shifts in soil microbial metabolic activities and community structures along a salinity gradient of irrigation water in a typical arid region of China[J]. Science of the Total Environment, 2017, 598: 64-70.
[52] Chambers L G, Osborne T Z, Reddy K R. Effect of salinity-altering pulsing events on soil organic carbon loss along an intertidal wetland gradient: a laboratory experiment[J]. Biogeochemistry, 2013, 115(1): 363-383.
[53] Wang H T, Gilbert J A, Zhu Y G, et al. Salinity is a key factor driving the nitrogen cycling in the mangrove sediment[J]. Science of the Total Environment, 2018, 631-632: 1342-1349.
[54] Kastl E M, Schloter-Hai B, Buegger F, et al. Impact of fertilization on the abundance of nitrifiers and denitrifiers at the root-soil interface of plants with different uptake strategies for nitrogen[J]. Biology and Fertility of Soils, 2015, 51(1): 57-64.
[55] Ngumbi E, Kloepper J. Bacterial-mediated drought tolerance: current and future prospects[J]. Applied Soil Ecology, 2016, 105: 109-125.
[56] 姚兰, 胡立煌, 张焕朝, 等. 黄山土壤细菌群落和酶活性海拔分布特征[J]. 环境科学, 2019, 40(2): 859-868.
Yao L, Hu L H, Zhang H C, et al. Elevational distribution characteristics of soil bacterial community and enzyme activities in mount Huangshan[J]. Environmental Science, 2019, 40(2): 859-868.
[57] 丁凯, 张毓婷, 张俊红, 等. 不同密度杉木林对林下植被和土壤微生物群落结构的影响[J]. 植物生态学报, 2021, 45(1): 62-73.
Ding K, Zhang Y T, Zhang J H, et al. Effects of Chinese fir plantations with different densities on understory vegetation and soil microbial community structure[J]. Chinese Journal of Plant Ecology, 2021, 45(1): 62-73.
[58] 林耀奔, 杨建辉, 叶艳妹. 盐碱地不同土地利用方式下土壤细菌群落结构多样性差异分析[J]. 环境科学学报, 2019, 39(4): 1266-1273.
Lin Y B, Yang J H, Ye Y M. Analysis on diversity of soil bacterial community under different land use patterns in saline-alkali land[J]. Acta Scientiae Circumstantiae, 2019, 39(4): 1266-1273.
[59] 何玉实, 何彤慧, 冯艳琼, 等. 鄂尔多斯台地盐沼滩涂湿地土壤细菌群落结构及特征[J]. 生态学报, 2022, 42(8): 3345-3355.
He Y S, He T H, Feng Y Q, et al. Characteristics and distribution of soil bacterial of salt marsh tidal wetland in Ordos Platform[J]. Acta Ecologica Sinica, 2022, 42(8): 3345-3355.
[60] 李茜, 徐瑞蔓, 陈迪, 等. 不同栽培年限人参根际土壤细菌群落与土壤理化性质和酶活性的相关性[J]. 植物营养与肥料学报, 2022, 28(2): 313-324.
Li X, Xu R M, Chen D, et al. Correlation of bacterial community with soil physicochemical properties and enzyme activities in rhizosphere soil under different cultivation years of ginseng (Panax ginseng C. A. Mey.)[J]. Journal of Plant Nutrition and Fertilizers, 2022, 28(2): 313-324.
[61] 郑佳华, 赵萌莉, 王琪, 等. 利用方式对大针茅草原土壤微生物群落结构及多样性的影响[J]. 生态学报, 2022, 42(12).
Zheng J H, Zhao M L, Wang Q, et al. Effects of management regime on soil microbial community structure and diversity of Stipa grandis grassland[J]. Acta Ecologica Sinica, 2022, 42(12). DOI:10.5846/stxb202103290818
[62] 韩文辉, 赵颖, 熊莲. 污灌区土壤微生物数量及其与土壤因子的关系[J]. 环境与可持续发展, 2017, 42(4): 230-233.
Han W H, Zhao Y, Xiong L. Soil microbial biomass and its relationship with soil factors in the sewage irrigation area[J]. Environment and Sustainable Development, 2017, 42(4): 230-233.
[63] 李新, 焦燕, 代钢, 等. 内蒙古河套灌区不同盐碱程度的土壤细菌群落多样性[J]. 中国环境科学, 2016, 36(1): 249-260.
Li X, Jiao Y, Dai G, et al. Soil bacterial community diversity under different degrees of saline-alkaline in the Hetao Area of Inner Mongolia[J]. China Environmental Science, 2016, 36(1): 249-260.