Effect of Long-term Biogas Slurry Application on Fungal Community in Rhizosphere Soil at Areca Taro Fields
-
摘要:目的 探究施用沼液对槟榔芋根际土壤真菌群落多样性的影响，为沼液更好的还田利用提供理论依据。方法 采集连续施用沼液6年和未施用沼液的槟榔芋种植地块土壤，通过Illumina Miseq高通量测序平台，分析土壤真菌群落多样性、组成及其与土壤理化性质间的关系，揭示影响真菌群落的主要土壤因子。结果 与未施用沼液土壤比较，连续施用沼液6年的土壤有机质、碱解氮和速效钾含量都有明显提高。连续施用沼液后土壤真菌物种OTUs数量、群落丰富度指数和多样性指数也明显增加。真菌门水平鉴定出子囊菌门（Ascomycota）、担子菌门（Basidiomycota）、被孢霉门（Mortierellomycota）、球囊菌门（Glomeromycota）等类群，其中子囊菌门为土壤中丰富度最高的真菌类群。连续施用沼液明显增加了被孢霉门中被孢霉纲的相对丰度，而球囊菌门中球囊菌纲则大幅减少；连续施用沼液后土壤真菌属水平优势类群也发生明显改变。冗余分析结果表明，引起土壤真菌群落多样性变化最主要的环境因子为碱解氮、速效钾和有机质。结论 长期施用沼液改变了槟榔芋根际土壤真菌所处的微环境，进而导致其群落组成和多样性发生变化。Abstract:Objective Effect of long-term biogas slurry application on the diversity of fungal community in rhizosphere soil at areca taro fields was studied.Method The taro planting plots that had been continuously fertilized with biogas slurry in the past 6 years were sampled to compare with areas without such treatment as control. Using the Illumina Miseq high-throughput sequencing platform, relationship between the fungal community and physicochemical characteristics of the soil was analyzed.Result The contents of organic matters, alkaline nitrogen, and available potassium in the soil were significantly increased by the continuous biogas slurry application. The fungal species OUT, community richness, and diversity at the sites also increased significantly. At phylum level, Ascomycota, Basidiomycota, Mortierellomycota, Glomeromycota, etc. were identified in the soils, with Ascomycota being the most abundant. After 6-year biogas slurry application, the relative abundance of Mortierellomycetes in Mortierellomycota increased significantly, while that of Glomeromycetes in Glomeromycota decreased significantly in the soil. The dominant fungi at genus level were also changed significantly from control. The redundancy analysis indicated that the major factors affecting the change of fungal community diversity in soil included alkaline nitrogen, available potassium, and organic matters.Conclusion Long-term biogas slurry fertilization altered the soil microenvironment leading to changes in the composition and diversity of rhizosphere fungal community.
-
0. 引言
【研究意义】近年来，随着畜禽养殖业的不断发展，畜禽废弃物的产量越来越大，沼气工程为处理畜禽养殖场废弃物提供了有效手段[1]。沼液是厌氧发酵后的副产物，含有丰富的营养元素，多以速效性及水溶性状态存在，易被作物吸收利用，是一种优质的有机肥料[2]。沼液施肥后可以增强土壤肥力，在提高作物产量、改善品质和防病抗逆等方面发挥重要作用[3]。通过沼液还田利用可以有效解决沼液的消纳问题，已成为当前畜禽养殖废弃物资源化利用的重要途径[4]。农业农村部在2018年重点推广的十大类型生态农业模式中，就有4类模式涉及沼液还田利用。沼液进入农田后会引起土壤环境发生系列改变，包括土壤的营养结构、土壤中重金属的形态和含量、土壤酸碱环境等，而土壤微生物能够敏感地反映出土壤环境的变化，是评价土壤质量的主要指标之一[5]。因此，研究沼液施肥对农田土壤微生物群落及微生态环境的影响，对于评价土壤质量演变、制定科学沼液还田方案具有重要的指导意义。【前人研究进展】目前关于沼液还田影响土壤微生物的研究主要集中在土壤微生物群落代谢功能多样性、土壤微生物数量、土壤细菌群落特征等方面。郑学博等[6]对旱地红壤的研究表明，沼液还田显著影响旱地红壤微生物群落的代谢活性和多样性。冯丹妮等[7]对水旱轮作耕层土壤的研究表明，长期沼液还田有利于土壤微生物数量和酶活性的提高，但对土壤微生态平衡及土壤质量存在潜在威胁。朱金山等[8]对稻田土壤的研究表明，不同沼液还田年限下稻田耕层土壤细菌群落组成及多样性均发生了明显变化。【本研究切入点】在众多研究中，有关沼液还田对土壤真菌群落特征影响的研究却较少报道。土壤真菌作为土壤微生物中主要的类群之一，是土壤中有机物的重要分解者，在陆地生态系统养分循环中发挥着重要作用，真菌多样性及其群落结构的变化是评价其所在生态系统健康稳定的关键性指标之一[9-10]。【拟解决的关键问题】为此，本文基于高通量测序技术，对采自连续施用沼液6年和未施用沼液的槟榔芋种植地块土壤真菌多样性及群落结构进行研究，探讨沼液施用对土壤真菌群落特征的影响，以及土壤真菌对由沼液施用引起的土壤理化性质变化的响应，为促进沼液资源化利用，维持畜禽养殖业的可持续发展提供理论依据。
1. 材料与方法
1.1 采样点基本情况
采样点位于建瓯市某养猪场粪污处理沼气工程周边的农田（27°13.30′N；118°12.80′E），该区域属于中亚热带海洋性季风气候，四季分明，雨量充足，春夏多雨，秋冬干燥，常年平均气温约18.80 ℃，降雨量约1 700 mm左右，无霜期约286 d。施肥用沼液来自该养猪场基地沼气工程充分发酵后经氧化塘处理的液体（沼液与氧化塘的水配比为1∶2），布管道至田里，按需实施漫灌施肥。供试沼液基本理化性质为：pH 7.99，化学需氧量（COD）285 mg·L−1，五日生化需氧量（BOD5）263 mg·L−1，氨氮475 mg·L−1，总氮978 mg·L−1，总磷22.50 mg·L−1。采样点土壤质地为砂壤土，主要种植作物为槟榔芋，每年只种1季，3月中下旬种植，11月初收获，冬季闲田。沼液灌溉的田块按照基肥和追肥方式灌入，在槟榔芋种植前，每667 m2灌入沼液量约15 t与田土混匀作基肥，随后每隔15 d进行1次追肥，每次追肥灌入沼液量约6 t，大概8月下旬进行最后1次追肥，整个生长期都不施用其他肥料，其他农艺管理措施均按照当地常规习惯统一实施。
1.2 样品的采集与制备
沼液施肥试验开始于2013年3月，于2018年10月（沼液施肥第6年），在槟榔芋种植地块上进行根际土壤样品的采集。采样时除去表层枯枝落叶等杂物，用取土器沿着植株的根部采集0～20 cm土层的土壤样品，标记为BS，同时采集周边一处地势较高未施用沼液的槟榔芋土壤作为对照样品，标记为CK，这2块槟榔芋田同属于一个农户，种植方式和田间管理方式相同。每块样地设置3份重复，每份重复样品均由5点取样法获得5个小样组成，现场混合均匀后装入无菌塑封袋中，放入带有冰块的保温箱中运回实验室。将取回的土样分成2份，1份用于土壤基本理化性质的测定，另一份存放于–80 ℃冰箱内用于DNA提取。
1.3 土壤理化性质测定
土壤pH值采用酸度计法（土∶水=1∶2.5）测定；有机质含量采用重铬酸钾氧化法测定；碱解氮采用扩散法测定；速效磷采用钼锑抗比色法测定；速效钾采用火焰光度计法测定。每个指标都进行全程同步试剂空白控制，每个样均设平行双样。
1.4 土壤基因组DNA提取
采用PowerSoil DNA Isolation Kit（MoBio Laboratories, Carlsbad，CA，USA）试剂盒，参照说明书，对土壤微生物总DNA进行提取。将提取的DNA用1%琼脂糖凝胶进行电泳，使用Nanodrop 2000检验其质量和浓度。
1.5 PCR扩增与高通量测序
选用引物ITS1-F（CTTGGTCATTTAGAGGAAGTAA）和ITS2-R（TGCGTTCTTCATCGATGC），对真菌ITS基因的ITS1可变区进行扩增，扩增程序为：95 ℃预变性5 min，30个循环（95 ℃变性30 s，55 ℃退火30 s，72 ℃延伸30 s），最后72 ℃延伸10 min。扩增体系总体积为20 μL，各反应组分为：4 μL 5×FastPfu Buffer缓冲液，1 μL引物（5 μM），2 μL dNTPs（2.5 mM），2 μL模板DNA，不足部分由无菌超纯水补足。反应结束后，扩增产物通过2%的琼脂糖凝胶电泳进行检测，采用QIAquick Gel Extraction Kit（QIAGEN, Germany）凝胶回收试剂盒进行纯化。将纯化后的扩增片段进行文库构建，构建好的文库经检测合格后，使用Illumina公司的Miseq PE250平台进行测序（由北京奥维森基因科技有限公司协助完成）。
1.6 数据分析
利用QIIME 1.8软件对测序原始数据进行质量筛查，去除barcode和primer并拼接后得到原始序列，经进一步去除嵌合体及短序列后得到高质量序列。使用UPARSE 10.0.240软件，对有效数据按97%的序列相似度进行归并和OTU（operational taxonomic unit）划分，统计样品中每个OTU的丰度信息。基于OTU的结果，计算Chao1指数、Observed species指数、PD whole tree指数和Shannon指数来进行生物多样性分析。采用RDP Classifier（http://sourceforge.net/projects/rdp-classifier/）方法对照Silva 128（http://www.arb-silva.de）数据库对OTU代表性序列进行物种分类注释，比对阈值设置为70%，并在不同分类水平上统计样本的群落组成。
样品在不同分类水平上的群落组成比较柱形图是通过OrginPro 2015软件得到的。利用R语言的Venn Diagram程序包和gplot程序包分别绘制韦恩图和热图，采用Canoco for Windows 4.5软件对真菌群落和理化因子进行典范对应分析。
2. 结果与分析
2.1 长期施用沼液对土壤理化性质的影响
对槟榔芋种植地块土壤pH分析结果表明：连续施用沼液6年后，土壤中pH有降低的趋势，但整体变化不显著（P＞0.05），所测土样均呈微酸性。进一步对槟榔芋地块土壤样本的主要化学参数分析表明：连续施用沼液6年后，土壤有机质、碱解氮和速效钾含量均有一定程度提高，增幅分别为4.18%、10.92%和48.95%（表1），与未施用沼液土壤相比差异均达到显著性水平（P＜0.05）。
表 1 不同施肥方式土壤理化性质Table 1. Physicochemical characteristics of soil under varied fertilizations样品编号
Samples IDpH 有机质
Organic matter/（g·kg−1）碱解氮
Alkali-hydrolysis N/（mg·kg−1）速效磷
Available P/（mg·kg−1）速效钾
Available K/（mg·kg−1）CK 5.29±0.01a 38.02±0.06b 229±0.28b 31.59±0.25a 190±0.76b BS 5.27±0.0a 39.61±0.01a 254±2.47 a 30.17±0.10a 283±2.09a 注：CK，未施用沼液土壤；BS，沼液施肥6年土壤；同列数据后不同小写字母表示处理间在0.05水平差异显著，表2同。
Note: CK, Soil with no biogas fertilization; BS, Soil with 6-year biogas slurry fertilization; Data with different lowercase letters on same column indicate significant difference among treatments at 0.05 level. Same for Table 2.2.2 长期施用沼液对土壤真菌多样性的影响
2.2.1 土壤真菌OTU水平分析
从未施用沼液（CK）和施用沼液（BS）两组土壤中，分别获得真菌总有效序列108 605条和82 479条，测序覆盖率在99%以上，说明测序数据量合理，能够比较真实的反映所测土壤样品中真菌群落组成情况。在97%相似性水平下，对样品的有效数据进行聚类，共划分出2 202个真菌OTUs。未施用沼液土壤（CK）样品中含有真菌种类1 839个，施用沼液土壤（BS）样品中含有1 844个，其中CK和BS有相同种类1 481个，特有种类分别为358个和363个（图1），表明不同施肥方式下土壤真菌OTUs组成有一定的差异。
![]()
图 1 不同施肥方式土壤真菌OTU分布韦恩图注：CK，未施用沼液土壤；BS，沼液施肥6年土壤，下同。Figure 1. Venn graph on fungi OTU distribution in soil under varied fertilizationsNote: CK, Soil with no biogas fertilization; BS, Soil with 6-year biogas slurry fertilization. Same for the following.2.2.2 土壤真菌多样性分析
基于OTU的结果，通过Chao1、Observed species、PD whole tree和Shannon等4种不同的计算指标来进行样品多样性分析，结果如表2所示。Chao1指数表征微生物物种的丰富度，用以估计群落中的OTU数目；Observed species指数用以估计样品中实际观测到的OTU数目；PD whole tree指数表征系统发育的多样性指数；Shannon指数表征样品中微生物群落的多样性程度。由表2可知，与未施用沼液土壤比较，长期施用沼液土壤Chao1和Observed species指数显著上升（P＜0.05），增幅分别为30.38%和26.36%，即真菌群落的物种丰富度增加；PD whole tree和Shannon指数也呈现上升的趋势，其中PD whole tree指数与未施用沼液土壤相比差异达显著性水平（P＜0.05），即长期施用沼液后，土壤真菌群落的多样性增加。
表 2 不同施肥方式土壤真菌多样性分析Table 2. Fungal diversity in soil under varied fertilizations样品编号 Samples ID Chao1 Observed species PD whole tree Shannon CK 1 623±14.72b 1 252±18.88b 248±2.14b 6.72±0.26a BS 2 116±57.86a 1 582±66.11a 305±6.34a 6.87±0.23a 2.2.3 土壤真菌群落组成分析
将2组土壤样品中得到的OTU，分别在门、纲、目、科、属5个水平上进行物种注释，统计不同组土壤样品中的真菌群落组成情况，图2-图4分别展示了门、纲、属3个不同水平的群落结构比较结果。
![]()
图 2 不同施肥方式土壤真菌群落在门水平的相对丰度Figure 2. Relative abundance of fungal community at phylum level in soil under varied fertilizations![]()
图 3 不同施肥方式土壤真菌群落在纲水平的相对丰度Figure 3. Relative abundance of fungal community at class level in soil under varied fertilizations![]()
图 4 不同施肥方式土壤真菌群落在属水平聚类热图Figure 4. Cluster heatmap of fungal community at genus level in soil under varied fertilizations从门的分类水平上看（图2），2组土壤真菌主要由子囊菌门（Ascomycota）、担子菌门（Basidiomycota）、被孢霉门（Mortierellomycota）、未被分类的真菌门（Fungi_unidentified）、球囊菌门（Glomeromycota）、罗兹菌门（Rozellomycota）等组成，其中优势菌门为子囊菌门、担子菌门、被孢霉门和未被分类的真菌门，它们占整体真菌群落总丰度的92%以上。对不同施肥方式土壤中真菌类群比较发现，长期施用沼液土壤中被孢霉门的相对丰度明显高于未施用沼液土壤，增幅为88.98%；而担子菌门和未被分类的真菌门的相对丰度却低于未施用沼液土壤，分别减少了13.35%和27.14%，可见施用沼液改变了土壤真菌优势菌门的相对丰度。此外，长期施用沼液后，土壤中的一些非优势菌门的相对丰度也发生了明显变化，如球囊菌门和罗兹菌门所占比例呈减少趋势，较未施用沼液土壤分别减少了88.70%和46.62%。
从纲的分类水平上看（图3），2组土壤真菌主要由粪壳菌纲（Sordariomycetes）、未被分类的真菌纲（Fungi unidentified）、被孢霉纲（Mortierellomycetes）、伞菌纲（Agaricomycetes）、座囊菌纲（Dothidemycetes）等组成，其中优势菌纲为粪壳菌纲、未被分类的真菌纲、被孢霉纲和伞菌纲，这些类群占整体真菌群落总丰度的67%以上。与未施用沼液土壤比较，长期施用沼液土壤中丰度显著增加的类群为被孢霉纲，增幅为90.12%，而一些非优势菌群如散囊菌纲（Eurotiomycetes）、球囊菌纲（Glomeromycetes）、古根菌纲（Archaeorhizomycetes）和微球黑粉菌纲（Microbotryomycetes）所占比例显著减少，较未施用沼液土壤分别减少了47.84%、90.14%、99.29%和97.68%。
为了更形象直观的呈现长期施用沼液后土壤真菌群落结构的变化，进一步选取相对丰度较高（＞1%）的真菌属进行聚类分析（图4）。
从图4中可以看出，2组土壤真菌群落组成既有一定相似性，又表现出了一定差异性。未被分类的真菌属Fungi unidentified和隶属子囊菌门的Pseudeurotium属聚为一类，在2组土壤样品中的相对丰度无明显变化。隶属子囊菌门的Penicillium属、Acremonium属和Sordaria属、隶属担子菌门的Psilocybe属、隶属球囊菌门的Pacispora属聚为一类，这些类群在两组土壤中的相对丰度大小为未施用沼液＞施用沼液，为未施用沼液土壤的优势菌属；而隶属子囊菌门的Rhodoveronaea属、Trichoderma属、Westerdykella属和Chaetomium属等、隶属担子菌门的Saitozyma属和Ramaria属、隶属被孢霉门的Mortierella属聚为一类，这些类群在长期施用沼液后比例明显增加，成为施用沼液土壤中的优势菌属。
2.3 土壤真菌群落与土壤理化性质的相关性
为了分析土壤理化因子对真菌群落的影响，在门水平上选取相对丰度大于0.1%的真菌类群与土壤的pH、有机质（OM）、碱解氮（AN）、速效磷（AP）和速效钾（AK）等5个理化因子进行典范对应分析。首先对土壤真菌群落做降趋势对应分析（Detrended correspondence analysis, DCA）。结果显示，4个排序轴中梯度最大长度为0.457，小于3，因此本研究选择基于线性模型的冗余分析（Reundancy analysis, RDA）进行排序（图5）。
![]()
图 5 不同施肥方式土壤真菌冗余分析注：OM，有机质；AN，碱解氮；AP，速效磷；AK，速效钾。Figure 5. Redundancy analysis on fungi in soil under varied fertilizationsNote: OM, Organic matter; AN, Alkali-hydrolysis N; AP, Available P; AK, Available K.从图5可以看出，未施用沼液和施用沼液土壤为2个相对独立的群落，表明它们之间的真菌群落差异较大，说明了施用沼液对土壤真菌群落的影响比较明显。RDA前2个排序轴共解释了98.0%的真菌群落变化，其中第1排序轴解释了84.5%的变异。经forward分析，对真菌群落影响程度排前3位的土壤因子依次为碱解氮、速效钾和有机质，蒙特卡罗检验显示这3个因子对土壤真菌群落影响显著（P＜0.05）。碱解氮、速效钾和有机质存在明显的协同作用，其含量与Mortierellomycota、Entorrhizomycota和Cercozoa等菌群丰度呈极显著或显著正相关，而与Glomeromycota丰度呈极显著负相关。
3. 讨论
施用沼液能够提高土壤养分含量，增加土壤中微生物的数量及多样性指数[11]。本研究中，多样性分析结果表明，长期沼液施肥明显增加了土壤真菌的物种丰富度和群落多样性。这主要是由于沼液的长期施用改善了土壤理化性质，为土壤真菌提供了良好的生长环境及充足的能源，从而刺激了真菌微生物的活性，增加其群落的丰富度和多样性。与本研究相似，郑学博等[6]的研究也表明沼液还田会提高土壤微生物群落功能的多样性。左狄等[12]关于不同深度稻田土壤的研究表明，沼液还田显著增加了各深度土壤中微生物的数量。而朱金山等[8]的研究发现，随着沼液施肥年限的增加，稻田土壤微生物的物种丰富度和群落多样性逐渐降低。冯丹妮等[13]关于水稻油菜轮作土壤的研究也发现，沼液连年施用不利于丰富耕层土壤细菌的群落多样性。这可能与土壤质量状况、不同的沼液成分或不同的浇灌措施等都有密切关系，如Jose等[14]研究发现，土壤在受到高浓度重金属污染后，细菌群落多样性会显著降低。
从群落组成上来看，本研究在2组土壤样品中，真菌鉴定出子囊菌门（Ascomycota）、担子菌门（Basidiomycota）、被孢霉门（Mortierellomycota）、球囊菌门（Glomeromycota）、罗兹菌门（Rozellomycota）等主要菌门，其中子囊菌门为两组土壤中丰富度最高的真菌类群。与众多研究结果一致[15-17]，子囊菌门在土壤生态系统中广泛存在，这可能是因为子囊菌可以产生大量的无性孢子，从而快速生长[18]。另外，本研究所测土壤均呈微酸性，这也就证实了子囊菌门为酸性土壤中的主要优势菌群[19]。同样，2组土壤鉴定出的优势菌纲有粪壳菌纲（Sordariomycetes）、被孢霉纲（Mortierellomycetes）和伞菌纲（Agaricomycetes），其中丰富度最高的优势菌为粪壳菌纲，这与徐雪雪等[20]、季凌飞等[21]的研究结果相类似，即该优势菌普遍存在于土壤中且丰富度最高。
本研究揭示，长期施用沼液改变了土壤真菌群落的相对丰度，提高了被孢霉门的相对丰度，降低了担子菌门、球囊菌门、罗兹菌门等的相对丰度。同样，长期施用沼液提高了被孢霉纲的相对丰度，降低了球囊菌纲的相对丰度。大量研究表明，不同施肥模式引起土壤环境条件的改变，会间接影响土壤微生物群落结构[22-24]。在本研究中，沼液施肥后槟榔芋根际土壤的理化性质也发生了不同程度的变化。RDA分析结果显示，土壤碱解氮、速效钾和有机质含量对土壤真菌群落影响显著，之前研究也报道了这些理化因子对土壤真菌群落结构的影响。如冯翠娥等[25]的研究表明碱解氮含量是影响土壤真菌群落组成最主要的因子之一；武俊男等[26]的研究表明速效钾含量是决定土壤真菌群落组成的关键因素之一；Liu等[27]、Luo等[28]认为有机碳含量是影响土壤真菌多样性的重要因素。本文中，长期施用沼液明显提高了土壤中碱解氮、速效钾和有机质含量，这些土壤因子与被孢霉门丰度呈极显著或显著正相关，与球囊菌门丰度却呈极显著负相关，这也就导致了长期施用沼液后被孢霉门丰度增加，而球囊菌门丰度减少。另外，有研究表明，球囊菌门中分布着大量的丛枝菌根真菌，连续施用沼液可能会降低土壤中丛枝菌根真菌含量，进而降低了球囊菌门的相对丰度[29-30]。
本研究发现，与未施用沼液土壤相比，长期施用沼液后土壤真菌属水平类群的相对丰度也发生明显变化。其中隶属子囊菌门的Rhodoveronaea属、Trichoderma属、Westerdykella属和Chaetomium属等、隶属担子菌门的Saitozyma属和Ramaria属、隶属被孢霉门的Mortierella属丰度明显增加。这主要是因为长期的沼液施用为这些菌群的生长提供了丰富的碳源和氮源，进而促进了其生长繁殖。然而，隶属子囊菌门的Penicillium属、Acremonium属和Sordaria属、隶属担子菌门的Psilocybe属、隶属球囊菌门的Pacispora属，这些类群在长期施用沼液后丰度明显降低，说明这些类群对沼液施肥引发的环境变化较为敏感。另外，有研究表明，Trichoderma、Chaetomium、Penicillium、Acremonium等菌属是目前报道的纤维素降解菌，对纤维素分解具有较强的作用[31]。本文施用沼液有效增加了Trichoderma和Chaetomium菌属的相对丰度，说明土壤施用沼液可能有利于某些纤维素降解菌的生长，这为后续从施用沼液土壤中筛选分离能高效分解纤维素的真菌提供了基础。然而，有关沼液施肥土壤中哪些菌为有益菌，哪些菌为有害菌，这些菌与槟榔芋生长发育之间的关系还需后期菌群功能分类分析进一步验证。
4. 结论
本文采用Illumina平台的MiSeq高通量测序技术，从真菌群落多样性、丰度及组成等方面研究了沼液施肥对槟榔芋种植地块土壤真菌的影响。结果表明，施用沼液明显提高了土壤有机质、碱解氮和速效钾含量，土壤理化性质的改变影响着土壤真菌多样性和群落结构特征。与未施用沼液土壤比较，长期施用沼液后土壤真菌物种OTUs数量、丰富度和群落多样性有明显增加，不同分类水平上各真菌类群所占比例也发生了改变。真菌门水平鉴定出子囊菌门（Ascomycota）、担子菌门（Basidiomycota）、被孢霉门（Mortierellomycota）、球囊菌门（Glomeromycota）等类群，其中子囊菌门为土壤中丰富度最高的真菌类群。与未施用沼液土壤比较，在长期施用沼液的土壤中，被孢霉门中被孢霉纲相对丰度明显增加，而球囊菌门中球囊菌纲丰度明显减少；长期施用沼液后土壤真菌属水平优势类群也发生明显改变。RDA分析结果表明，施用沼液组土壤真菌群落结构与未施用沼液组差异较大，土壤碱解氮、速效钾和有机质含量是引起土壤真菌群落变化的重要环境因素。
-
![]()
图 1 不同施肥方式土壤真菌OTU分布韦恩图
注：CK，未施用沼液土壤；BS，沼液施肥6年土壤，下同。
Figure 1. Venn graph on fungi OTU distribution in soil under varied fertilizations
Note: CK, Soil with no biogas fertilization; BS, Soil with 6-year biogas slurry fertilization. Same for the following.
![]()
图 2 不同施肥方式土壤真菌群落在门水平的相对丰度
Figure 2. Relative abundance of fungal community at phylum level in soil under varied fertilizations
![]()
图 3 不同施肥方式土壤真菌群落在纲水平的相对丰度
Figure 3. Relative abundance of fungal community at class level in soil under varied fertilizations
![]()
图 4 不同施肥方式土壤真菌群落在属水平聚类热图
Figure 4. Cluster heatmap of fungal community at genus level in soil under varied fertilizations
![]()
图 5 不同施肥方式土壤真菌冗余分析
注：OM，有机质；AN，碱解氮；AP，速效磷；AK，速效钾。
Figure 5. Redundancy analysis on fungi in soil under varied fertilizations
Note: OM, Organic matter; AN, Alkali-hydrolysis N; AP, Available P; AK, Available K.
表 1 不同施肥方式土壤理化性质
Table 1 Physicochemical characteristics of soil under varied fertilizations
样品编号
Samples IDpH 有机质
Organic matter/（g·kg−1）碱解氮
Alkali-hydrolysis N/（mg·kg−1）速效磷
Available P/（mg·kg−1）速效钾
Available K/（mg·kg−1）CK 5.29±0.01a 38.02±0.06b 229±0.28b 31.59±0.25a 190±0.76b BS 5.27±0.0a 39.61±0.01a 254±2.47 a 30.17±0.10a 283±2.09a 注：CK，未施用沼液土壤；BS，沼液施肥6年土壤；同列数据后不同小写字母表示处理间在0.05水平差异显著，表2同。
Note: CK, Soil with no biogas fertilization; BS, Soil with 6-year biogas slurry fertilization; Data with different lowercase letters on same column indicate significant difference among treatments at 0.05 level. Same for Table 2.
下载: 导出CSV
表 2 不同施肥方式土壤真菌多样性分析
Table 2 Fungal diversity in soil under varied fertilizations
样品编号 Samples ID Chao1 Observed species PD whole tree Shannon CK 1 623±14.72b 1 252±18.88b 248±2.14b 6.72±0.26a BS 2 116±57.86a 1 582±66.11a 305±6.34a 6.87±0.23a
下载: 导出CSV
-
[1] 张东, 徐甦, 陈斌, 等. 畜禽粪便沼气工程处理技术 [J]. 浙江农业科学, 2012(2):223−227. ZHANG D, XU S, CHEN B, et al. Treatment method of livestock manure via biogas project [J]. Journal of Zhejiang Agricultural Sciences, 2012(2): 223−227.(in Chinese
[2] ZIRKLER D, PETERS A, KAUPENJOHANN M. Elemental composition of biogas residues: Variability and alteration during anaerobic digestion [J]. Biomass and Bioenergy, 2014, 67: 89−98. DOI: 10.1016/j.biombioe.2014.04.021
[3] 王国海, 翟静. 沼液对土壤改良的作用 [J]. 北京农业, 2015 (3): 93. WANG G H, ZHAI J. The effects of biogas slurry on soil improvements[J]. Beijing Agriculture, 2015 (3): 93. (in Chinese).
[4] 曾锦, 徐锐, 梁高飞, 等. 畜禽养殖废弃物资源化利用技术及推广模式研究进展 [J]. 畜牧与饲料科学, 2018, 39(8):56−63. ZENG J, XU R, LIANG G F, et al. Research advances on the resource recycling technology and extension model of livestock and poultry wastes [J]. Animal Husbandry and Feed Science, 2018, 39(8): 56−63.(in Chinese
[5] NANNIPIERI P, ASCHER J, CECCHERINI M T, et al. Microbial diversity and soil functions [J]. European Journal of Soil Science, 2003, 54(4): 655−670. DOI: 10.1046/j.1351-0754.2003.0556.x
[6] 郑学博, 樊剑波, 崔键, 等. 沼液还田对旱地红壤微生物群落代谢与多样性的影响 [J]. 生态学报, 2016, 36(18):5865−5875. ZHENG X B, FAN J B, CUI J, et al. Analysis on metabolic characteristics and functional diversity of soil edaphon communities in upland red soil under biogas slurry application [J]. Acta Ecologica Sinica, 2016, 36(18): 5865−5875.(in Chinese
[7] 冯丹妮, 伍钧, 杨刚, 等. 连续定位施用沼液对水旱轮作耕层土壤微生物区系及酶活性的影响 [J]. 农业环境科学学报, 2014, 33(8):1644−1651. FENG D N, WU J, YANG G, et al. Influence of long-term applications of biogas slurry on microbial community composition and enzymatic activities in surface soil under rice-rape rotation [J]. Journal of Agro-Environment Science, 2014, 33(8): 1644−1651.(in Chinese
[8] 朱金山, 张慧, 马连杰, 等. 不同沼灌年限稻田土壤微生物群落分析 [J]. 环境科学, 2018, 39(5):2400−2411. ZHU J S, ZHANG H, MA L J, et al. Diversity of the microbial community in rice paddy soil with biogas slurry irrigation analyzed by illumina sequencing technology [J]. Environmental Science, 2018, 39(5): 2400−2411.(in Chinese
[9] 张晶, 张惠文, 李新宇, 等. 土壤真菌多样性及分子生态学研究进展 [J]. 应用生态学报, 2004, 15(10):1958−1962. ZHANG J, ZHANG H W, LI X Y, et al. Research advances in soil fungal diversity and molecular ecology [J]. Chinese Journal of Applied Ecology, 2004, 15(10): 1958−1962.(in Chinese
[10] HARRIS J. Soil microbial communities and restoration ecology: Facilitators or followers? [J]. Science, 2009, 325(5940): 573−574. DOI: 10.1126/science.1172975
[11] 张无敌, 尹芳, 徐锐, 等. 沼液对土壤生物学性质的影响 [J]. 湖北农业科学, 2009, 48(10):2403−2407. ZHANG W D, YIN F, XU R, et al. Effect of biogas liquid on biological properties of soil [J]. Hubei Agricultural Sciences, 2009, 48(10): 2403−2407.(in Chinese
[12] 左狄, 吕卫光, 李双喜, 等. 沼液还田对稻田土壤养分与氮循环微生物的影响 [J]. 上海农业学报, 2018, 34(2):55−59. ZUO D, LÜ W G, LI S X, et al. Effects of biogas slurry on soil nutrients and nitrogen cycling microorganism in paddy field [J]. Acta Agriculturae Shanghai, 2018, 34(2): 55−59.(in Chinese
[13] 冯丹妮, 伍钧, 杨虎德. 连续施用沼液水稻油菜轮作耕层土壤的细菌多样性 PCR-DGGE分析 [J]. 甘肃农业科技, 2019(8):41−49. FENG D N, WU J, YANG H D. Bacterial diversity in topsoil of continuous application of biogas slurry rice-rape rotation analysised using PCR-DGGE [J]. Gansu Agricultural Science and Technology, 2019(8): 41−49.(in Chinese
[14] JOSE J, GIRIDHAR R, ANAS A, et al. Heavy metal pollution exerts reduction/adaptation in the diversity and enzyme expression profile of heterotrophic bacteria in Cochin estuary, India [J]. Environmental Pollution, 2011, 159(10): 2775−2780. DOI: 10.1016/j.envpol.2011.05.009
[15] 宋三多. 沼肥施用对成都平原稻麦轮作制作物生长及土壤微生物特性影响 [D]. 雅安: 四川农业大学, 2017. SONG S D. Effects of biogas manure on crop growth and soil microbial properties under rice-wheat rotation system in Chengdu plain [D]. Yaan, China: Sichuan Agricultural University, 2017. (in Chinese).
[16] DING H Y, ALI A, CHENG Z H. Dynamics of a soil fungal community in a three-year green garlic/cucumber crop rotation system in northwest China [J]. Sustainability, 2018, 10(5): 1391. DOI: 10.3390/su10051391
[17] 李敏, 闫伟. 海拔对乌拉山油松根围真菌群落结构的影响 [J]. 菌物学报, 2019, 38(11):1992−2006. LI M, YAN W. Effects of altitude on rhizosphere fungal community structure of Pinus tabulaeformis in Wula Mountain, China [J]. Mycosystema, 2019, 38(11): 1992−2006.(in Chinese
[18] AL-SADI A M, AL-MAZROUI S S, PHILLIPS A J L. Evaluation of culture-based techniques and 454 pyrosequencing for the analysis of fungal diversity in potting media and organic fertilizers [J]. Journal of Applied Microbiology, 2015, 119(2): 500−509. DOI: 10.1111/jam.12854
[19] 张岩, 邓娇娇, 朱文旭, 等. 辽西北风沙区典型人工防护林土壤微生物群落功能特征 [J]. 生态学杂志, 2019, 38(6):1688−1696. ZHANG Y, DENG J J, ZHU W X, et al. Functional characteristics of soil microbial community of the typical artificial shelter forests in the windy-sandy region of northwest Liaoning Province, China [J]. Chinese Journal of Ecology, 2019, 38(6): 1688−1696.(in Chinese
[20] 徐雪雪, 王东, 秦舒浩, 等. 沟垄覆膜连作马铃薯根际土壤真菌多样性分析 [J]. 水土保持学报, 2015, 29(6):301−306, 310. XU X X, WANG D, QIN S H, et al. Diversity analysis of fungal communities in potato continuous cropping soil under different patterns of ridge-furrow film mulching [J]. Journal of Soil and Water Conservation, 2015, 29(6): 301−306, 310.(in Chinese
[21] 季凌飞, 倪康, 马立锋, 等. 不同施肥方式对酸性茶园土壤真菌群落的影响 [J]. 生态学报, 2018, 38(22):8158−8166. JI L F, NI K, MA L F, et al. Effect of different fertilizer regimes on the fungal community of acidic tea-garden soil [J]. Acta Ecologica Sinica, 2018, 38(22): 8158−8166.(in Chinese
[22] LI J, COOPER J, LIN Z, et al. Soil microbial community structure and function are significantly affected by long-term organic and mineral fertilization regimes in the North China Plain [J]. Applied Soil Ecology, 2015, 96: 75−87. DOI: 10.1016/j.apsoil.2015.07.001
[23] HUHE, CHEN X J, HOU F J, et al. Bacterial and fungal community structures in loess plateau grasslands with different grazing intensities [J]. Frontiers in Microbiology, 2017(8): 606. DOI: 10.3389/fmicb.2017.00606
[24] 陈浩, 魏立本, 王亚麒, 等. 烤烟不同种植施肥模式对土壤养分、酶活性及细菌多样性的影响 [J]. 南方农业学报, 2019, 50(5):982−989. CHEN H, WEI L B, WANG Y Q, et al. Effects of different planting and fertilizing modes on soil nutrient, enzyme activity and bacterial diversity of tobacco [J]. Journal of Southern Agriculture, 2019, 50(5): 982−989.(in Chinese
[25] 冯翠娥, 岳思君, 简阿妮, 等. 硒砂瓜连作对土壤真菌群落结构的影响 [J]. 中国生态农业学报(中英文), 2019, 27(4):537−544. FENG C E, YUE S J, JIAN A, et al. The effect of continuous cropping of selenium melon on soil fungal community structure [J]. Chinese Journal of Eco-Agriculture, 2019, 27(4): 537−544.(in Chinese
[26] 武俊男, 刘昱辛, 周雪, 等. 基于Illumina MiSeq测序平台分析长期不同施肥处理对黑土真菌群落的影响 [J]. 微生物学报, 2018, 58(9):1658−1671. JUNNAN W, YUXIN L, XUE Z, et al. Effects of long-term different fertilization on soil fungal communities in black soil based on the Illumina Mi Seq platform [J]. Acta Microbiologica Sinica, 2018, 58(9): 1658−1671.(in Chinese
[27] LIU J J, SUI Y Y, YU Z H, et al. Soil carbon content drives the biogeographical distribution of fungal communities in the black soil zone of northeast China [J]. Soil Biology & Biochemistry, 2015, 83: 29−39.
[28] LUO P Y, HAN X R, WANG Y, et al. Influence of long-term fertilization on soil microbial biomass, dehydrogenase activity, and bacterial and fungal community structure in a brown soil of northeast China [J]. Annals of Microbiology, 2015, 65(1): 533−542. DOI: 10.1007/s13213-014-0889-9
[29] WRIGHT S H A, BERCH S M, BERBEE M L. The effect of fertilization on the below-ground diversity and community composition of ectomycorrhizal fungi associated with western hemlock (Tsuga heterophylla) [J]. Mycorrhiza, 2009, 19(4): 267−276. DOI: 10.1007/s00572-008-0218-x
[30] 王幼珊, 刘润进. 球囊菌门丛枝菌根真菌最新分类系统菌种名录 [J]. 菌物学报, 2017, 36(7):820−850. WANG Y S, LIU R J. a checklist of arbuscular mycorrhizal fungi in the recent taxonomic system of Glomeromycota [J]. Mycosystema, 2017, 36(7): 820−850.(in Chinese
[31] 熊冬梅, 周红丽. 纤维素降解菌群的研究进展 [J]. 酿酒科技, 2011(5):94−97. XIONG D M, ZHOU H L. Research progress in cellulose-degrading bacteria community [J]. Liquor-Making Science & Technology, 2011(5): 94−97.(in Chinese
-
期刊类型引用(2)
1. 曾志浩，袁宗胜，陈雪莹，林鸿艳，黄晓南，刘芳. 3种红树植物根际与非根际土壤细菌群落结构及多样性特征. 福建农业学报. 2022(06): 809-816 . 
本站查看
2. 雷菲，张冬明，吴宇佳，谭皓，吉清妹，潘孝忠. 化肥减量配施有机肥对樱桃番茄产量、品质和微生物群落结构的影响. 江西农业大学学报. 2021(06): 1269-1277 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 1293
- HTML全文浏览量: 375
- PDF下载量: 41
- 被引次数: 5
