环境科学  2023, Vol. 44 Issue (4): 2338-2347   PDF    
菌渣与化肥配施对稻田土壤微生物群落组成及多样性的影响
耿和田1, 王旭东1, 石思博2, 叶正钱1, 周文晶1     
1. 浙江农林大学环境与资源学院, 浙江省污染土壤修复重点实验室, 杭州 311300;
2. 南京农业大学资源与环境科学学院土壤生态实验室, 南京 210095
摘要: 菌渣是一种独特而丰富的有机物料,与化肥配施不仅能改良土壤质量还可以调控微生物群落.然而,土壤细菌和真菌对菌渣与化肥配施的响应是否一致仍不清楚.在水稻田间长期定位试验条件下,设置化肥水平(C)0%、50%和100%,菌渣相对用量(F)0%、50%和100%各3个水平,共9个处理,对土壤肥力与微生物群落的相关指标进行测定.结果表明,土壤全氮(TN)在C0F100处理中最高,碳氮比(C/N)、全磷(TP)、可溶性碳(DOC)和有效磷(AP)在C100F100处理中最高,土壤有机碳(SOC)、碱解氮(AN)、速效钾(AK)和pH在C50F100处理中最高,较对照分别增加了55.56%、26.18%、26.46%、17.13%、279.54%、85.57%、41.61%、29.33%和4.62%.菌渣与化肥配施后,不同处理土壤细菌和真菌α-多样性存在显著变化,与对照C0F0处理相比较,各处理细菌β-多样性并未发生显著变化,却使真菌β-多样性发生了明显差异,且C50F100处理显著降低了真菌子囊菌门(Ascomycota)和粪壳菌纲(Sordariomycetes)的相对丰度.随机森林预测模型表明,AP和C/N分别是细菌和真菌α-多样性的主要驱动因子,AN、pH、SOC和DOC是细菌β-多样性的主要驱动因子,而AP和DOC是真菌β-多样性的主要驱动因子.相关性分析表明,真菌子囊菌门(Ascomycota)和粪壳菌纲(Sordariomycetes)与SOC、TN、TP、AN、AP、AK和C/N显著负相关.置换多元方差分析表明,菌渣对土壤肥力指标、细菌门及纲水平上优势物种、真菌门及纲水平优势物种、细菌多样性和真菌多样性的变化贡献率分别为46.35%、18.47%、41.57%、23.84%和10.42%,菌渣和化肥的交互效应对细菌和真菌多样性变化的贡献率分别为9.90%和35.00%.综上所述,施用菌渣较化肥在影响土壤肥力指标含量和微生物群落变化方面更具优势.
关键词: 菌渣      化肥      稻田土壤      细菌多样性      真菌多样性      微生物群落     
Effects of Combined Application of Fungal Residue and Chemical Fertilizer on Soil Microbial Community Composition and Diversity in Paddy Soil
GENG He-tian1 , WANG Xu-dong1 , SHI Si-bo2 , YE Zheng-qian1 , ZHOU Wen-jing1     
1. Zhejiang Provincial Key Laboratory of Contaminated Soil Remediation, College of Environmental and Resource Science, Zhejiang Agricultural and Forestry University, Hangzhou 311300, China;
2. Soil Ecology Laboratory, College of Resources and Environmental Sciences, Nanjing Agricultural University, Nanjing 210095, China
Abstract: Fungal residue is a unique abundant organic material undervalued in agricultural production. The application of chemical fertilizer combined with fungal residue can not only improve soil quality but also regulate the microbial community. However, it is unclear whether the response of soil bacteria and fungi to the combined application of fungal residue and chemical fertilizer is consistent. Therefore, a long-term positioning experiment in a rice field was conducted with a total of nine treatments. Chemical fertilizer (C) and fungal residue (F) were applied at 0, 50%, and 100% to evaluate 1 the change in soil fertility properties and microbial community structure and 2 the main driving factors of soil microbial diversity and species composition. The results showed that soil total nitrogen (TN) was highest after treatment C0F100 (55.56% higher than in the control), and the carbon to nitrogen ratio (C/N), total phosphorus (TP), dissolved organic carbon (DOC), and available phosphorus (AP) contents were highest after treatment with C100F100(26.18%, 26.46%, 17.13%, and 279.54% higher than in the control, respectively). The amounts of soil organic carbon (SOC), available nitrogen (AN), available potassium (AK), and pH were highest after treatment with C50F100 (85.57%, 41.61%, 29.33%, and 4.62% higher than in the control, respectively). Following the application of fungal residue with chemical fertilizer, there were significant changes in the α-diversity of bacteria and fungi in each treatment. Compared with that of the control (C0F0), different long-term applications of fungal residue with chemical fertilizer did not significantly change soil bacterial β-diversity but resulted in significant differences in fungal β-diversity, and the relative abundance of soil fungal Ascomycota and Sordariomycetes significantly decreased after the application of C50F100. The random forest prediction model indicated that AP and C/N were the main driving factors of bacterial and fungal α-diversity, respectively, and AN, pH, SOC, and DOC were the main driving factors of bacterial β-diversity, whereas AP and DOC were the main driving factors of fungal β-diversity. Correlation analysis suggested that the relative abundance of soil fungal Ascomycota and Sordariomycetes had a significantly negative correlation with SOC, TN, TP, AN, AP, AK, and C/N. PERMANOVA showed that variation in soil fertility properties, dominant species of soil bacteria at the phylum and class level, and dominant species of soil fungi at the phylum and class level were all best explained by fungal residue (46.35%, 18.47%, and 41.57%, respectively), and variation in bacterial diversity was best explained by fungal residue (23.84%) and to a lesser extent by the interaction between fungal residue and chemical fertilizer (9.90%). In contrast, the variation in fungal diversity was best explained by the interaction between fungal residue and chemical fertilizer (35.00%) and to a lesser extent by fungal residue (10.42%). In conclusion, the application of fungal residue has more advantages than chemical fertilizer in influencing soil fertility properties and microbial community structure changes.
Key words: fungal residue      chemical fertilizer      paddy soil      bacterial diversity      fungal diversity      microbial community     

土壤微生物是土壤的重要组成部分, 是维持土壤健康的关键[1, 2].土壤微生物群落可能对植物生长、凋落物分解、驱动养分循环和碳、氮周转起着关键的作用[3~6].Wei等[7]报道了土壤微生物群落的组成和功能决定了植物是否可以抵抗土传病害.有研究表明, 微生物多样性可以显著影响土壤中碳的周转[8~10], Chen等[11]揭示了微生物Shannon多样性指数可以作为植被变化下SOM动态的可测量指标, Xu等[12]的研究发现在底土中, 当土壤细菌和真菌多样性较高时, SOC分解表现出更强的抗变暖能力.因此, 土壤微生物是土壤质量变化的衡量指标.

由于土壤微生物对环境的敏感性[13, 14], 长期施肥等农艺活动会导致土壤环境因子发生变化, 从而影响微生物群落[15, 16].有研究表明, 施氮改变了土壤微生物群落组成, 使细菌群落从以寡养菌群为主转变为以共养菌群为主[17, 18], Dai等[19]的研究发现施氮会降低细菌多样性, 对于真菌, 施氮会降低担子菌(Basidiomycota)的相对丰度[20].土壤微生物群落对有机肥施用也极为敏感[21], 有机物料和有机肥可以缓解化肥对微生物多样性的负面影响[22, 23].Liu等[24]报道了与单施化肥相比, 有机替代增加了微生物多样性.马龙等[25]的研究发现有机肥替代部分化肥可以提高土壤细菌和真菌丰度, 改变土壤优势菌群的组成.因此, 探究有机物料与化肥配施对土壤肥力指标含量和微生物群落的影响, 对正确揭示土壤肥力指标含量与微生物群落之间的关系具有重要意义.

菌渣是一种独特而丰富的有机物料, 是农业生产中被低估的生物质资源[26], 将其还田能为土壤提供大量的营养物质.我国是世界上最大的食用菌生产国, 但菌渣的利用率仅为33%[27], 不恰当处置可能会造成生态环境污染[28].以往对菌渣还田的研究主要集中在土壤养分、土壤酶和作物产量变化方面[29~31].菌渣还田后土壤微生物群落可能会发生变化, 然而, 菌渣还田后对微生物群落的影响目前尚未报道.为此, 在水稻田间长期定位试验条件下, 设置不同比例的菌渣与化肥配施处理, 利用Illumina MiSeq高通量测序技术测定细菌和真菌群落, 通过置换多元方差分析和随机森林预测模型, 定量评估菌渣对土壤肥力指标和微生物群落的影响, 及其微生物多样性和物种组成的主要驱动因子.

1 材料与方法 1.1 试验场地

本试验于2017年6月至2017年12月在中国浙江省嘉兴市水稻田(北纬30°37′~30°40′, 东经120°39′~120°44′, 平均海拔为4 m)进行.该地区四季分明, 属东亚季风区, 年平均降水量为1 168.6 mm, 主要集中在5~9月, 年平均气温为15.9℃, 无霜期为230 d, 年平均日照为2 017 h.长期定位试验始于2010年, 供试土壤试验前的基本理化性质为:pH值6.80, ω[有机碳(SOC)]25.16 g·kg-1, ω[碱解氮(AN)]94.29 mg·kg-1, ω[有效磷(AP)]37.01 mg·kg-1, ω[速效钾(AK)]127.22 mg·kg-1.土壤耕作方式是休耕轮作.菌渣选用种植黑木耳后的桑枝屑物料经发酵后的产物.主要养分为:ω(SOC)451.8 g·kg-1, ω[全氮(TN)]11.4 g·kg-1, ω[全磷(TP)]1.0 g·kg-1, ω[全钾(TK)]6.0 g·kg-1, 碳氮比(C/N)为39.6.成分由纤维素(28.21%)、半纤维素(20.16%)、木质素(15.11%)、灰分(14.83%)和蛋白质(14.3%)组成.供试水稻品种为甬优1540. 2017年6月15日播种育苗, 7月13日插秧移栽秧苗, 12月13日收获.

1.2 试验设计

本试验设置施用化肥C0、C50和C100, 施用菌渣F0、F50和F100, 采用双因素随机区组排列, 共设9个处理, 每个处理3次重复, 共27个小区, 每个小区面积为20 m2.菌渣施用量为0(0%)、10 t·hm-2(50%)和20 t·hm-2(100%), 化肥施用量为当地常规施肥量的0%、50%和100%.各处理具体如下: 不施化肥和菌渣(C0F0)为对照, 施用菌渣10 t·hm-2(C0F50), 施用菌渣20 t·hm-2(C0F100), 50%的化肥用量(C50F0), 50%的化肥用量配施菌渣10 t·hm-2(C50F50), 50%的化肥用量配施菌渣20 t·hm-2(C50F100), 100%的化肥用量(C100F0), 100%的化肥用量配施菌渣10 t·hm-2(C100F50), 100%的化肥用量配施菌渣20 t·hm-2(C100F100).结合当地单季稻种植制度, 水稻移栽前将土壤进行翻耕, 同时施入不同比例混合的菌渣和化肥(基肥).常规化肥100%施用量具体如下:基肥施碳酸氢铵(含氮量17%)300 kg·hm-2;分蘖初期追施尿素(含氮量46%)150 kg·hm-2;分蘖盛期追施控释复合肥(N ∶P ∶K=22 ∶11 ∶13)225 kg·hm-2, 50%化肥处理则将100%化肥处理减半进行, 0%化肥处理不施肥, 其他田间管理策略遵循传统栽培技术.

1.3 土壤采样与分析

土壤样品于水稻收获后, 从每个小区0~20 cm土层多点采集土壤, 清除石块和植物残留物, 彻底混合取样土壤, 通过2 mm筛网进行筛分, 并分为3部分.第一部分用于土壤可溶性碳(DOC)分析, 第二部分存于-84℃用于土壤DNA提取, 第三部分通风阴干用于土壤理化性质分析.

1.4 测定方法

SOC、TP、TN、pH、AN、AP和AK采用常规土壤农化分析方法测定, DOC按水土比为2 ∶1比例浸提.

土壤DNA提取及测序:使用DNA试剂盒(MP Biomedicals, Santa Ana, CA, USA)从1.0 g新鲜土壤中提取DNA总量.为了生成细菌和真菌高通量测序的扩增子库, 细菌16S rRNA基因用引物F515(5′-GTGCCAGCMGCCGCGGTAA-3′)和R806(5′-GGAC TACHVGGGTWTCTAAT-3′)在V4区进行PCR扩增, 真菌ITS基因用引物ITS1F(5′- CTTGGTCA TTTAGAGGAAGTAA-3′)和ITS2(5′-GCTGCGTTCTT CATCGATGC-3′)在ITS1区进行PCR扩增.PCR扩增完成后, 使用QIAquick PCR纯化试剂盒(OMEGA, USA)对产物进行纯化, 用Qubit2.0荧光仪(Invitrogen, USA)进行定量.然后将纯化后的扩增子以等量浓度汇集, 用NEB Next UltraTM DNA文库预备试剂盒(Annoline Biotechnology Ltd., London, UK)构建文库.用Agilent 2100生物分析仪(Agilent Technologies, Inc., Santa Clara, CA, USA)检测文库的最终质量和浓度, 用KAPA文库定量试剂盒(Kapa Bio, Wilmington, MA, USA)进行测定.文库测序所有准备工作均在Illumina MiSeq平台上进行(Shanghai Majorbio Bio-Pharm Technology Co., Ltd., Shanghai, China).

1.5 统计分析

使用SPSS 26.0(IBM, Chicago, IL, USA)对不同处理的SOC、DOC、TN、TP、AN、AP、AK、pH和微生物α-多样性指数进行单因素方差分析和Duncan显著性检测(统计显著水平为P < 0.05), 以比较不同处理间的差异.利用R4.1.2的“vegan”包计算微生物α-多样性指数和绘制非度量多维尺度图以辨别微生物群落的差异.利用R4.1.2的“rfPermute”包预测SOC、DOC、TN、TP、AN、AP、AK和pH对细菌和真菌α-多样性及β-多样性的重要性, 以确定细菌和真菌α-多样性及β-多样性的主要驱动因子.使用R4.1.2的“vegan”包中的“adonis”函数执行置换多元方差(PERMANOVA)分析, 定量评估菌渣和化肥对土壤肥力指标含量和微生物群落的影响.使用R4.1.2的“pheatmap”包对细菌和真菌的门及纲水平上优势物种与土壤肥力指标进行相关性分析.使用R4.1.2的“ggplot2”包绘制细菌和真菌的门及纲水平上物种组成堆叠图.数据可视化在RStudio上操作, 用Adobde Illustrator CC 2018(Adobe Inc., San Jose, CA, USA)进行图形排版.

2 结果与分析 2.1 菌渣与化肥配施对土壤肥力指标的影响

菌渣与化肥配施后, 各处理土壤肥力指标均有所增加(表 1).TN在C0F100处理中最高, C/N、TP、DOC和AP在C100F100处理中最高, pH、SOC、AN和AK均在C50F100处理中最高, 较对照分别增加了55.56%、26.18%、26.46%、17.13%、279.54%、4.62%、85.57%、41.61%和29.33%.在化肥施用C0、C50和C100水平下, 各处理土壤肥力指标的含量均随菌渣施用量的增加而增加.此外, 菌渣对土壤肥力指标含量变化的贡献率为46.35%[图 1(a)].

表 1 土壤肥力指标1) Table 1 Soil fertility properties

*为P < 0.05, **为P < 0.01, ***为P < 0.001 图 1 基于PERMANOVA值的菌渣和化肥对土壤肥力指标和微生物群落变化的贡献率 Fig. 1 PERMANOVA values showing the contribution rate that fungal residue and chemical fertilizer contributed to the variation in soil fertility properties and microbial community

2.2 菌渣与化肥配施对土壤微生物多样性的影响

菌渣与化肥配施后, 土壤细菌和真菌多样性指数存在显著差异.细菌Simpson和Chao指数分别在C0F100和C0F50处理中最大[图 2(a)2(b)], 真菌Simpson和Chao指数分别在C50F0和C100F0处理中最大[图 2(c)2(d)], 较对照(C0F0)分别增加了64.01%、2.82%、85.46%和50.89%.真菌Simpson和Chao指数在C50F100处理中最小[图 2(c)2(d)], 与对照相比分别降低90.60%和29.37%.另外, 各处理土壤真菌群落与对照明显分离[图 3(b)].菌渣对土壤细菌和真菌多样性变化的贡献率分别为23.84%和10.42%[图 1(b)1(c)], 且菌渣和化肥的交互效应对真菌多样性变化的贡献率为35.00%[图 1(c)].

不同小写字母表示不同处理在P < 0.05水平显著差异 图 2 菌渣与化肥配施对土壤细菌和真菌α-多样性的影响 Fig. 2 Effects of combination of fungal residue with chemical fertilizer on soil bacterial and fungal α-diversity

图 3 菌渣与化肥配施对细菌和真菌β-多样性的影响 Fig. 3 Effects of combination of fungal residue with chemical fertilizer on soil bacterial and fungal β-diversity

2.3 菌渣与化肥配施对土壤微生物物种组成的影响

菌渣与化肥配施后, 各处理土壤细菌门及纲水平上物种组成无明显变化[图 4(a)4(b)].然而, 真菌门及纲水平上物种存在明显差异[图 4(c)4(d)], 真菌子囊菌门(Ascomycota, 40.40% ~70.57%)和粪壳菌纲(Sordariomycetes, 30.17% ~55.03%)分别为最优势菌门和菌纲, 其相对丰度在对照C0F0处理中最高, 在C50F100处理中最低, 且在施用化肥C0和C100水平下随着菌渣用量的增加而减小, C50水平下呈先增大后减小[图 4(c)4(d)].施用菌渣对真菌门及纲水平上优势物种变化的贡献率为41.57%[图 1(e)].

图 4 菌渣与化肥配施对土壤微生物物种组成的影响 Fig. 4 Effects of combination of fungal residue with chemical fertilizer on soil microbial species composition

2.4 土壤微生物多样性与土壤肥力指标的相关性分析

SOC、AP和TN对土壤细菌α-Simpson多样性有极显著贡献[P < 0.01, 图 5(a)], AP对α-Chao多样性有显著贡献[P < 0.05, 图 5(c)].AP和SOC对真菌α-Chao多样性呈显著贡献[P < 0.05, 图 5(d)], C/N对真菌α-Simpson和Chao多样性分别有显著[P < 0.05, 图 5(b)]和极显著贡献[P < 0.01, 图 5(d)].AN、pH、SOC和DOC对细菌β-多样性有显著贡献[P < 0.05, 图 5(e)], AP和DOC对真菌β-多样性有显著贡献[P < 0.05, 图 5(f)].

变量的均方误差增加率用于评估这些预测值的重要性, 均方误差增加率越大, 则意味着预测值越重要;*为P < 0.05, **为P < 0.01 图 5 土壤肥力指标作为土壤微生物α-和β-多样性驱动因子的随机森林预测重要性均值 Fig. 5 Random forest mean predicted importance of the soil fertility properties as drivers for soil microbial α- and β-diversity

2.5 土壤细菌和真菌的优势物种与土壤肥力指标的相关性分析

土壤细菌变形菌门(Proteobacteria)和δ-变形菌纲(δ-Proteobacteria)与DOC、TN、TP和AP显著负相关[图 6(a)], 绿弯菌门(Chloroflexi)和厌氧绳菌纲(Anaerolineae)与SOC、DOC、TN、TP、AN、AP和AK显著正相关[图 6(a)].真菌子囊菌门(Ascomycota)和粪壳菌纲(Sordariomycetes)与SOC、TN、AN、AP和C/N显著负相关[图 6(b)], 接合菌门(Zygomycota)、单子菌门(Basidiomycota)和接合菌纲(norank_Zygomycota)都与SOC、DOC、AN、AP、AK和C/N显著正相关[图 6(b)].

a1.Nitrospirae, a2.Actinobacteriota, a3.Acidobacteriota, a4.Chloroflexi, a5.Proteobacteria, a6.Actinobacteria, a7.α-Proteobacteria, a8.δ-Proteobacteria, a9.Anaerolineae, a10.Acidobacteria, b1.Rozellomycota, b2.Zygomycota, b3. unclassified_Fungi, b4.Basidiomycota, b5.Ascomycota, b6.norank_Zygomycota, b7.Dothideomycetes, b8. unclassified_Fungi, b9.Agaricomycetes, b10.Sordariomycetesc; 红色表示显著正相关, 蓝色表示显著负相关, 灰色表示无显著相关;*为P < 0.05, **为P < 0.01, ***为P < 0.001 图 6 土壤细菌和真菌优势物种与土壤肥力指标的相关性 Fig. 6 Correlation between soil bacterial and fungal dominant species and soil fertility properties

3 讨论 3.1 菌渣与化肥配施对土壤肥力指标的影响

Tarkalson等[32]的研究发现, 施用化肥是导致土壤酸化的主要因素, 有研究表明土壤中添加有机物料不仅能调节pH, 还能提高土壤肥力[33].化肥减量配施20 t·hm-2菌渣能够提高土壤pH, 说明菌渣与化肥配施可以缓解土壤酸化[34]. SOC是维持土壤肥力和微生物群落稳定的重要角色[35], 相比单施化肥, 菌渣与化肥配施后SOC含量显著增加, 这与温广婵等[36]的研究结果一致, 马欣等[37]的研究也发现有机物料配施化肥对农田土壤固碳效果更好.DOC是表征土壤活性碳库的重要指标, 虽然其含量占土壤活性碳库的比例较小, 但却是土壤碳库中最为重要和活跃的部分, 且在土壤物质转化和养分循环中起着关键的作用[38~40].顾春朝等[41]的研究发现, 有机肥配施化肥更有利于土壤DOC的释放, 这与本研究的结果类似.N、P和K被认为是大多数生态系统中主要的营养元素以及养分限制因素[42, 43], 土壤N、P和K含量高低不仅影响着农作物的产量和品质, 而且还是影响微生物群落结构的重要因子[4], 有研究表明AP显著增加会抑制微生物生长[44].菌渣与化肥配施后, 较对照土壤TN、AN、TP、AP和AK都显著增加(P < 0.05, 表 1), 这些结果与诸多研究一致[30, 31, 34, 35], 土壤C/N是土壤质量变化的敏感指标之一, 在本研究中土壤C/N随着菌渣添加量的增加而增加, 土壤C/N高低对作物生长发育和土壤养分平衡具有重要的影响, 且决定着有机氮的矿化以及NH4+的硝化作用[45, 46], 另外, 还会介导微生物对有机质分解过程中养分的释放[45].

3.2 菌渣与化肥配施对土壤细菌和真菌多样性的影响

土壤微生物多样性在可持续农业中发挥着重要作用[5, 6].本研究发现, 菌渣与化肥配施后, 微生物多样性存在显著变化, 细菌α-多样性在施用化肥水平下随着菌渣用量的增加而增加[图 1(b)], 这与Liu等[24]秸秆替代部分化肥增加了细菌多样性类似.相反, 化肥减量配施20 t·hm-2菌渣显著降低了真菌α-多样性[图 1(c)1(d)].细菌和真菌对施用菌渣的反应不同, 这可能是细菌和真菌群落获取营养素资源的偏好和策略不同.Six等[47]的研究表明, 与真菌相比较, 细菌更喜欢容易分解的资源.因此, 菌渣对真菌类群的选择性影响更大, 从而降低了真菌α-多样性.对于细菌α-多样性增加可能是菌渣添加后为细菌提供了大量的营养物质, 有利于细菌的生长.秦红灵等[48]的研究表明土壤AP增加能够显著增加细菌群落多样性, 在本研究中随机森林预测模型表明土壤AP对细菌α-多样性有显著贡献(P < 0.05), 这与姜霓雯等[49]的研究结果类似, 说明土壤AP是细菌α-多样性的主要驱动因子.本研究发现AN、SOC、pH和DOC是细菌β-多样性的主要驱动因子, 与Zheng等[50]的报道类似, 陈桂鲜等[51]的研究也表明, DOC是细菌β-多样性的主要驱动因子, 可能解释为DOC易于被微生物分解利用, 为微生物的活动提供了能量.然而, 在一项关于全球变化因子对土壤微生物多样性和功能的影响报道中, 土壤pH主导着α-多样性的反应[52], 但在本研究中并非如此.化肥减量配施20 t·hm-2菌渣降低真菌α-多样性可能归因于菌渣本身含有丰富的碳源, 高量菌渣施用增加了土壤C/N, 土壤C/N是决定微生物多样性的重要因素[53], N源不足从而导致真菌α-多样性降低; 随机森林模型表明, C/N对真菌α-多样性有显著贡献[P < 0.05, 图 5(c)5(d)], 与Ning等[54]的研究结果一致, 在本研究中AP和DOC是真菌β-多样性的主要驱动因子.综上所述, 菌渣对土壤细菌多样性的变化发挥了重要作用, 可能是通过增加AP的含量影响细菌α-多样性, AN、SOC、pH和DOC驱动了细菌β-多样性变化.然而, 菌渣和化肥的交互效应介导了真菌多样性变化, 且化肥减量配施20 t·hm-2菌渣降低了真菌α-多样性, C/N在其中发挥了重要的作用, AP和DOC主导了真菌β-多样性变化.

3.3 菌渣与化肥配施对细菌和真菌优势物种组成的影响

微生物群落对维持土壤生态系统功能至关重要, 且受到土地利用、耕作和施肥等的影响[55~57].在本研究中与细菌相比, 真菌在门及纲水平上物种变化更明显, 这可能是真菌对环境的敏感性, 菌渣与化肥配施后土壤环境因子发生明显变化, 其丰度受到影响[14].子囊菌门(Ascomycota)是所有处理中最优势的菌门, 这与Chen等[58]和Bei等[59]的研究结果一致.化肥减量配施20 t·hm-2菌渣降低了子囊菌门(Ascomycota)和粪壳菌纲(Sordariomycetes)的相对丰度[图 4(c)4(d)], 但Wu等[60]的研究表明, 化肥减量与有机物料相结合显著增加了子囊菌门(Ascomycota)的相对丰度.子囊菌门(Ascomycota)的相对丰度降低可能归因于两方面:首先, 可能是过高的土壤养分抑制了子囊菌门(Ascomycota)生长[61], 子囊菌门(Ascomycota)与SOC、TN、TP、AN、AP、AK、pH和C/N呈显著负相关能解释这一观点[图 6(b)].其次, 可解释为与菌渣本身有较高的C/N密切相关, 子囊菌门(Ascomycota)是土壤中有机物质的主要分解者, 对腐烂有机基质具有降解能力[62], 子囊菌门(Ascomycota)分解需要大量的N, 可能由于N缺乏抑制了子囊菌门(Ascomycota)的活性, 从而导致子囊菌门(Ascomycota)的相对丰度降低.另外, 本研究发现子囊菌门(Ascomycota)与担子菌门(Basidiomycota)之间可能存在竞争关系, 这与Ye等[63]的研究结果类似.担子菌门(Basidiomycota)与子囊菌门(Ascomycota)一样, 是典型的腐生菌门[64], 具有环境友好和分解有机物质的特点[63], 在降解土壤中高木质素含量的凋落物方面发挥了重要的作用[65].综上所述, 化肥减量配施20 t·hm-2菌渣通过增加土壤养分和C/N抑制了子囊菌门(Ascomycota)和粪壳菌纲(Sordariomycetes)的相对丰度.

4 结论

(1) 与对照相比较, 菌渣与化肥配施显著增加了土壤肥力指标的含量, 菌渣主导了土壤肥力指标含量的变化.

(2) 菌渣对土壤细菌多样性的变化发挥了重要作用, 可能是通过增加AP的含量影响细菌α-多样性, AN、SOC、pH和DOC驱动了细菌β-多样性的变化.然而, 菌渣和化肥的交互效应介导了真菌多样性变化, 且化肥减量配施20 t·hm-2菌渣降低了真菌α-多样性, C/N在其中发挥了重要的作用, AP和DOC主导了真菌β-多样性的变化.

(3) 化肥减量配施20 t·hm-2菌渣降低了真菌子囊菌门(Ascomycota)和粪壳菌纲(Sordariomycetes)的相对丰度, 且施用菌渣较化肥在影响土壤细菌和真菌的门及纲水平上优势物种变化方面更具优势.

参考文献
[1] 朱永官, 彭静静, 韦中, 等. 土壤微生物组与土壤健康[J]. 中国科学: 生命科学, 2021, 51(1): 1-11.
Zhu Y G, Peng J J, Wei Z, et al. Linking the soil microbiome to soil health[J]. Scientia Sinica Vitae, 2021, 51(1): 1-11.
[2] Saleem M, Hu J, Jousset A. More than the sum of its parts: microbiome biodiversity as a driver of plant growth and soil health[J]. Annual Review of Ecology, Evolution, and Systematics, 2019, 50: 145-168. DOI:10.1146/annurev-ecolsys-110617-062605
[3] Fang W S, Song Z X, Tao S, et al. Biochar mitigates the negative effect of chloropicrin fumigation on beneficial soil microorganisms[J]. Science of the Total Environment, 2020, 738. DOI:10.1016/j.scitotenv.2020.139880
[4] 马垒, 李燕, 魏建林, 等. 长期秸秆还田对潮土真菌群落、酶活性和小麦产量的影响[J]. 环境科学, 2022, 43(10): 4755-4764.
Ma L, Li Y, Wei J L, et al. Effects of long-term straw returning on fungal community, enzyme activity and wheat yield in a fluvo-aquic soil[J]. Environmental Science, 2022, 43(10): 4755-4764.
[5] Hallama M, Pekrun C, Lambers H, et al. Hidden miners-the roles of cover crops and soil microorganisms in phosphorus cycling through agroecosystems[J]. Plant and Soil, 2019, 434(1-2): 7-45. DOI:10.1007/s11104-018-3810-7
[6] 孙倩, 吴宏亮, 陈阜, 等. 不同轮作模式下作物根际土壤养分及真菌群落组成特征[J]. 环境科学, 2020, 41(10): 4682-4689.
Sun Q, Wu H L, Chen F, et al. Characteristics of soil nutrients and fungal community composition in crop rhizosphere under different rotation patterns[J]. Environmental Science, 2020, 41(10): 4682-4689.
[7] Wei Z, Gu Y A, Friman V P, et al. Initial soil microbiome composition and functioning predetermine future plant health[J]. Science Advances, 2019, 5(9). DOI:10.1126/sciadv.aaw0759
[8] Domeignoz-Horta L A, Pold G, Liu X J A, et al. Microbial diversity drives carbon use efficiency in a model soil[J]. Nature Communications, 2020, 11(1). DOI:10.1038/s41467-020-17502-z
[9] Juarez S, Nunan N, Duday A C, et al. Soil carbon mineralisation responses to alterations of microbial diversity and soil structure[J]. Biology and Fertility of Soils, 2013, 49(7): 939-948. DOI:10.1007/s00374-013-0784-8
[10] Tardy V, Spor A, Mathieu O, et al. Shifts in microbial diversity through land use intensity as drivers of carbon mineralization in soil[J]. Soil Biology and Biochemistry, 2015, 90: 204-213. DOI:10.1016/j.soilbio.2015.08.010
[11] Chen S T, Feng X, Lin Q M, et al. Pool complexity and molecular diversity shaped topsoil organic matter accumulation following decadal forest restoration in a karst terrain[J]. Soil Biology and Biochemistry, 2022, 166. DOI:10.1016/j.soilbio.2022.108553
[12] Xu M, Li X L, Kuyper T W, et al. High microbial diversity stabilizes the responses of soil organic carbon decomposition to warming in the subsoil on the Tibetan Plateau[J]. Global Change Biology, 2021, 27(10): 2061-2075. DOI:10.1111/gcb.15553
[13] 林耀奔, 叶艳妹, 杨建辉, 等. 土地整治对土壤微生物多样性的影响分析[J]. 环境科学学报, 2019, 39(8): 2644-2653.
Lin Y B, Ye Y M, Yang J H, et al. The effect of land consolidation on soil microbial diversity[J]. Acta Scientiae Circumstantiae, 2019, 39(8): 2644-2653.
[14] Ciccolini V, Bonari E, Pellegrino E. Land-use intensity and soil properties shape the composition of fungal communities in Mediterranean peaty soils drained for agricultural purposes[J]. Biology and Fertility of Soils, 2015, 51(6): 719-731. DOI:10.1007/s00374-015-1013-4
[15] 邹湘, 易博, 张奇春, 等. 长期施肥对稻田土壤微生物群落结构及氮循环功能微生物数量的影响[J]. 植物营养与肥料学报, 2020, 26(12): 2158-2167.
Zou X, Yi B, Zhang Q C, et al. Effects of long-term fertilization on the microbial community structure and the population of N cycle-related functional microorganism in paddy soil[J]. Journal of Plant Nutrition and Fertilizers, 2020, 26(12): 2158-2167.
[16] 李琳, 向丹, 武亚芬, 等. 长期不同施肥方式对日光温室番茄土壤养分和微生物群落结构的影响[J]. 应用生态学报, 2022, 33(2): 415-422.
Li L, Xiang D, Wu Y F, et al. Effects of long-term different fertilization patterns on soil nutrients and microbial community structure of tomato in a solar greenhouse[J]. Chinese Journal of Applied Ecology, 2022, 33(2): 415-422.
[17] Fierer N, Lauber C L, Ramirez K S, et al. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients[J]. The ISME Journal, 2012, 6(5): 1007-1017. DOI:10.1038/ismej.2011.159
[18] Ling N, Chen D M, Guo H, et al. Differential responses of soil bacterial communities to long-term N and P inputs in a semi-arid steppe[J]. Geoderma, 2017, 292: 25-33. DOI:10.1016/j.geoderma.2017.01.013
[19] Dai Z M, Su W Q, Chen H H, et al. Long-term nitrogen fertilization decreases bacterial diversity and favors the growth of Actinobacteria and Proteobacteria in agro-ecosystems across the globe[J]. Global Change Biology, 2018, 24(8): 3452-3461. DOI:10.1111/gcb.14163
[20] Zak D R, Pregitzer K S, Burton A J, et al. Microbial responses to a changing environment: implications for the future functioning of terrestrial ecosystems[J]. Fungal Ecology, 2011, 4(6): 386-395. DOI:10.1016/j.funeco.2011.04.001
[21] Bello A, Wang B, Zhao Y, et al. Composted biochar affects structural dynamics, function and co-occurrence network patterns of fungi community[J]. Science of the Total Environment, 2021, 775. DOI:10.1016/j.scitotenv.2021.145672
[22] Sun R B, Dsouza M, Gilbert J A, et al. Fungal community composition in soils subjected to long-term chemical fertilization is most influenced by the type of organic matter[J]. Environmental Microbiology, 2016, 18(12): 5137-5150. DOI:10.1111/1462-2920.13512
[23] Sun R B, Zhang X X, Guo X S, et al. Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw[J]. Soil Biology and Biochemistry, 2015, 88: 9-18. DOI:10.1016/j.soilbio.2015.05.007
[24] Liu H W, Du X F, Li Y B, et al. Organic substitutions improve soil quality and maize yield through increasing soil microbial diversity[J]. Journal of Cleaner Production, 2022, 347. DOI:10.1016/j.jclepro.2022.131323
[25] 马龙, 高伟, 栾好安, 等. 基于宏基因组学方法分析施肥模式对设施菜田土壤微生物群落的影响[J]. 植物营养与肥料学报, 2021, 27(3): 403-416.
Ma L, Gao W, Luan H A, et al. Soil microbial community characteristics in greenhouse vegetable production under different fertilization patterns based on metagenomic analysis[J]. Journal of Plant Nutrition and Fertilizers, 2021, 27(3): 403-416.
[26] Zhou J M, Wang L T, Wang H M, et al. Effects of different ratios of pig manure to fungus residue on physicochemical parameters during composting[J]. Journal of the Air & Waste Management Association, 2016, 66(5): 499-507.
[27] 兰良程. 中国食用菌产业现状与发展[J]. 中国农学通报, 2009, 25(5): 205-208.
Lan L C. The status quo and development of China's edible mushroom industry[J]. Chinese Agricultural Science Bulletin, 2009, 25(5): 205-208.
[28] Shi S B, Wang X D, Ye Z Q, et al. Effect of the combined application of fungal residue and chemical fertilizers on the mineralization of soil organic carbon in paddy fields[J]. Environmental Science and Pollution Research, 2019, 26(23): 23292-23304. DOI:10.1007/s11356-019-05603-2
[29] Huang C, Deng L J, Gao X S, et al. Effects of fungal residues return on soil enzymatic activities and fertility dynamics in a paddy soil under a rice-wheat rotation in Chengdu Plain[J]. Soil and Tillage Research, 2010, 108(1-2): 16-23. DOI:10.1016/j.still.2010.03.011
[30] 赵自超, 赵时锋, 张宏启, 等. 菌渣还田对设施瓜菜产量、品质和土壤肥力的影响[J]. 中国农学通报, 2021, 37(19): 112-118.
Zhao Z C, Zhao S F, Zhang H Q, et al. Effect of mushroom residue on cucumber and melon yield and quality and soil fertility in greenhouse[J]. Chinese Agricultural Science Bulletin, 2021, 37(19): 112-118.
[31] 滕青, 曾梦凤, 林慧凡, 等. 菌渣还田对生菜生长、土壤养分及酶活性的影响研究[J]. 中国农学通报, 2020, 36(6): 30-36.
Teng Q, Zeng M F, Lin H F, et al. Influences of application of mushroom residue on lettuce growth, soil nutrients and enzyme activities[J]. Chinese Agricultural Science Bulletin, 2020, 36(6): 30-36.
[32] Tarkalson D D, Payero J O, Hergert G W, et al. Acidification of soil in a dry land winter wheat-sorghum/corn-fallow rotation in the semiarid U.S. Great Plains[J]. Plant and Soil, 2006, 283(1-2): 367-379. DOI:10.1007/s11104-006-0027-y
[33] Kögel-Knabner I, Amelung W, Cao Z H, et al. Biogeochemistry of paddy soils[J]. Geoderman, 2010, 157(1-2): 1-14. DOI:10.1016/j.geoderma.2010.03.009
[34] 胡杨勇, 马嘉伟, 叶正钱, 等. 稻耳轮作制度下连续菌渣还田对土壤肥力性状的影响[J]. 水土保持学报, 2013, 27(6): 172-176.
Hu Y Y, Ma J W, Ye Z Q, et al. Effects of continuous application of edible fungus residue on soil fertility properties under rice-edible fungus rotation system[J]. Journal of Soil and Water Conservation, 2013, 27(6): 172-176.
[35] Bastida F, Eldridge D J, García C, et al. Soil microbial diversity-biomass relationships are driven by soil carbon content across global biomes[J]. The ISME Journal, 2021, 15(7): 2081-2091.
[36] 温广蝉, 叶正钱, 王旭东, 等. 菌渣还田对稻田土壤养分动态变化的影响[J]. 水土保持学报, 2012, 26(3): 82-86.
Wen G C, Ye Z Q, Wang X D, et al. Effects of edible fungus residue on dynamic changes of soil nutrients in paddy field[J]. Journal of Soil and Water Conservation, 2012, 26(3): 82-86.
[37] 马欣, 魏亮, 唐美玲, 等. 长期不同施肥对稻田土壤有机碳矿化及激发效应的影响[J]. 环境科学, 2018, 39(12): 5680-5686.
Ma X, Wei L, Tang M L, et al. Effects of varying long-term fertilization on organic carbon mineralization and priming effect of paddy soil[J]. Environmental Science, 2018, 39(12): 5680-5686.
[38] Li Y, Wu J S, Liu S L, et al. Is the C ∶N ∶P stoichiometry in soil and soil microbial biomass related to the landscape and land use in southern subtropical China?[J]. Global Biogeochemical Cycles, 2012, 26(4). DOI:10.1029/2012gb004399
[39] 章晓芳, 郑生猛, 夏银行, 等. 红壤丘陵区土壤有机碳组分对土地利用方式的响应特征[J]. 环境科学, 2020, 41(3): 1466-1473.
Zhang X F, Zheng S M, Xia Y H, et al. Responses of soil organic carbon fractions to land use types in hilly red soil regions, China[J]. Environmental Science, 2020, 41(3): 1466-1473.
[40] 石思博, 王旭东, 叶正钱, 等. 菌渣化肥配施对稻田土壤微生物量碳氮和可溶性碳氮的影响[J]. 生态学报, 2018, 38(23): 8612-8620.
Shi S B, Wang X D, Ye Z Q, et al. Effects of the combination of fungal residue and chemical fertilizer on soil microbial biomass carbon and nitrogen and dissolved organic carbon and nitrogen in paddy soil[J]. Acta Ecologica Sinica, 2018, 38(23): 8612-8620.
[41] 顾春朝, 傅民杰, 刘春海, 等. 施肥类型和种植年限对淹水稻田土壤可溶性有机碳的影响[J]. 广东农业科学, 2015, 42(12): 73-77.
Gu C Z, Fu M J, Liu C H, et al. Effects of different fertilizations and cultivation years on dissolved organic carbon (DOC) in paddy soil under water-logging condition[J]. Guangdong Agricultural Sciences, 2015, 42(12): 73-77.
[42] 宋佳杰, 徐郗阳, 白金泽, 等. 秸秆还田配施化肥对土壤养分及冬小麦产量的影响[J]. 环境科学, 2022, 43(9): 4839-4847.
Song J J, Xu X Y, Bai J Z, et al. Effects of straw returning and fertilizer application on soil nutrients and winter wheat yield[J]. Environmental Science, 2022, 43(9): 4839-4847.
[43] Yuan X, Qin W K, Xu H, et al. Sensitivity of soil carbon dynamics to nitrogen and phosphorus enrichment in an alpine meadow[J]. Soil Biology and Biochemistry, 2020, 150. DOI:10.1016/j.soilbio.2020.107984
[44] Sundqvist M K, Liu Z F, Giesler R, et al. Plant and microbial responses to nitrogen and phosphorus addition across an elevational gradient in subarctic tundra[J]. Ecology, 2014, 95(7): 1819-1835.
[45] 王建林, 钟志明, 王忠红, 等. 青藏高原高寒草原生态系统土壤碳氮比的分布特征[J]. 生态学报, 2014, 34(22): 6678-6691.
Wang J L, Zhong Z M, Wang Z H, et al. Soil C/N distribution characteristics of alpine steppe ecosystem in Qinhai-Tibetan Plateau[J]. Acta Ecologica Sinica, 2014, 34(22): 6678-6691.
[46] Wang N, Xu R K, Li J Y. Amelioration of an acid ultisol by agricultural by-products[J]. Land Degradation & Development, 2011, 22(6): 513-518.
[47] Six J, Frey S D, Thiet R K, et al. Bacterial and fungal contributions to carbon sequestration in agroecosystems[J]. Soil Science Society of America Journal, 2006, 70(2): 555-569.
[48] 秦红灵, 袁红朝, 张慧, 等. 红壤坡地利用方式对土壤细菌群落结构的影响[J]. 土壤学报, 2011, 48(3): 594-602.
Qin H L, Yuan H Z, Zhang H, et al. Soil bacteria community structure in upland red soil in relation to land use pattern[J]. Acta Pedologica Sinica, 2011, 48(3): 594-602.
[49] 姜霓雯, 梁辰飞, 张勇, 等. 土地利用变化后不同种植年限香榧土壤微生物群落的组成及多样性[J]. 环境科学, 2022, 43(1): 530-539.
Jiang N W, Liang C F, Zhang Y, et al. Microbial composition and diversity in soil of Torreya grandis cv. merrillii relative to different cultivation years after land use conversion[J]. Environmental Science, 2022, 43(1): 530-539.
[50] Zheng Q, Hu Y T, Zhang S S, et al. Soil multifunctionality is affected by the soil environment and by microbial community composition and diversity[J]. Soil Biology and Biochemistry, 2019, 136. DOI:10.1016/j.soilbio.2019.107521
[51] 陈桂鲜, 吴传发, 葛体达, 等. 土壤多功能性对微生物多样性降低的响应[J]. 环境科学, 2022, 43(11): 5274-5285.
Chen G X, Wu C F, Ge T D, et al. Response of soil multifunctionality to reduced microbial diversity[J]. Environmental Science, 2022, 43(11): 5274-5285.
[52] Zhou Z H, Wang C K, Luo Y Q. Meta-analysis of the impacts of global change factors on soil microbial diversity and functionality[J]. Nature Communications, 2020, 11(1). DOI:10.1038/s41467-020-16881-7
[53] 杨佩文, 李铭刚, 张庆, 等. 洱海流域稻油轮作农田土壤微生物多样性特征[J]. 中国土壤与肥料, 2021(4): 59-68.
Yang P W, Li M G, Zhang Q, et al. Molecular ecological characteristics of soil microorganisms in rice-rape ratation farmland in Erhai Basin[J]. Soils and Fertilizers Sciences in China, 2021(4): 59-68.
[54] Ning Q, Chen L, Jia Z J, et al. Multiple long-term observations reveal a strategy for soil pH-dependent fertilization and fungal communities in support of agricultural production[J]. Agriculture, Ecosystems & Environment, 2020, 293. DOI:10.1016/j.agee.2020.106837
[55] Malik A A, Puissant J, Buckeridge K M, et al. Land use driven change in soil pH affects microbial carbon cycling processes[J]. Nature Communications, 2018, 9(1). DOI:10.1038/s41467-018-05980-1
[56] Deng F B, Wang H J, Xie H T, et al. Low-disturbance farming regenerates healthy deep soil toward sustainable agriculture-Evidence from long-term no-tillage with stover mulching in Mollisols[J]. Science of the Total Environment, 2022, 825. DOI:10.1016/j.scitotenv.2022.153929
[57] Geisseler D, Scow K M. Long-term effects of mineral fertilizers on soil microorganisms-A review[J]. Soil Biology and Biochemistry, 2014, 75: 54-63.
[58] Chen C, Zhang J N, Lu M, et al. Microbial communities of an arable soil treated for 8 years with organic and inorganic fertilizers[J]. Biology and Fertility of Soils, 2016, 52(4): 455-467.
[59] Bei S K, Zhang Y L, Li T T, et al. Response of the soil microbial community to different fertilizer inputs in a wheat-maize rotation on a calcareous soil[J]. Agriculture, Ecosystems & Environment, 2018, 260: 58-69.
[60] Wu X, Zhang T, Zhao J N, et al. Variation of soil bacterial and fungal communities from fluvo-aquic soil under chemical fertilizer reduction combined with organic materials in north China plain[J]. Journal of Soil Science and Plant Nutrition, 2021, 21(1): 349-363.
[61] Zhou J, Jiang X, Zhou B K, et al. Thirty four years of nitrogen fertilization decreases fungal diversity and alters fungal community composition in black soil in northeast China[J]. Soil Biology and Biochemistry, 2016, 95: 135-143.
[62] Treseder K K, Maltz M R, Hawkins B A, et al. Evolutionary histories of soil fungi are reflected in their large-scale biogeography[J]. Ecology Letters, 2014, 17(9): 1086-1093.
[63] Ye G P, Lin Y X, Luo J F, et al. Responses of soil fungal diversity and community composition to long-term fertilization: field experiment in an acidic Ultisol and literature synthesis[J]. Applied Soil Ecology, 2020, 145. DOI:10.1016/j.apsoil.2019.06.008
[64] Guo J J, Liu W B, Zhu C, et al. Bacterial rather than fungal community composition is associated with microbial activities and nutrient-use efficiencies in a paddy soil with short-term organic amendments[J]. Plant and Soil, 2018, 424(1-2): 335-349.
[65] Liu G Y, Chen L L, Shi X R, et al. Changes in rhizosphere bacterial and fungal community composition with vegetation restoration in planted forests[J]. Land Degradation & Development, 2019, 30(10): 1147-1157.