随着农业规模化种植的进程加快,土传病害愈发严重,其具有不易被发现、强传染性、爆发性和摧毁性,治理困难
[1]。农业生产中为追求作物高产,盲目使用和滥用农药、化肥的情况时有发生,致使土壤有机养分降低及微生物群落结构失衡、病虫害频发,最终导致农作物、蔬菜和中药材品质下降
[1-2]。为了解决上述问题,寻找对生态友好的有机肥成为重要途径,已有研究发现蚯蚓肥及其衍生物,不仅可以增加作物生长和产量,还可以可持续地抑制土传病害发生
[3]。2022年蚯蚓堆肥技术被农业农村部列为4种典型农村有机废弃物资源化利用技术模式之一,而华北平原作为重要农牧区,存在着大量农牧废弃物资源,促成该区域蚯蚓养殖技术成熟,蚯蚓养殖规模位居全国前列,成为国内主要优质蚯蚓肥供应地
[4]。因此,华北平原蚯蚓肥应用潜力巨大,蚯蚓堆肥有助于推动该地区生态循环农业的发展。
蚯蚓肥是蚯蚓在降解有机废弃物后的产物,是一种集植物营养、植物保护和生物活性激素为一体的优质长效有机肥和土壤改良剂
[5-6]。施用蚯蚓肥一方面可以通过富含的营养成分和生物活性物质改变土壤理化性状;另一方面,与传统堆肥相比,蚯蚓肥具有出色的生物学特性,能改善并扩大土壤中有益微生物群落的多样性和数量;此外,利用自带的微生物直接参与土壤生物化学活动,同时还能拮抗植物病害
[7-9]。蚯蚓肥因富含大量有益细菌而成为抑制土传病害和调控微生物区系的重要手段
[2-3]。
蚯蚓肥的使用方法存在差异,多数研究倾向于选择使用蚯蚓肥浸提液
[3,10]。蚯蚓肥浸提液作为蚯蚓肥的衍生物,是蚯蚓肥与水混合并经过发酵制成
[11]。蚯蚓肥中的可溶性营养物质和微生物转移到液体溶液中,同样适用于农田生产
[12-14]。研究表明蚯蚓肥浸提液能提高种子活力,促进幼苗生长并增加叶绿素含量
[15]。蚯蚓肥浸提液也能促进微生物的活性和多样性
[13]。Singh等
[16]研究发现,蚯蚓肥浸提液的应用不仅减少了白化病和果实畸形等生理障碍,而且还减少了灰霉病的发生。因此,无论使用固体蚯蚓肥还是蚯蚓肥浸提液,对促进植物的生长和保护都能发挥出巨大的潜力。
目前多认为蚯蚓肥和蚯蚓肥浸提液发挥作用主要是其中含有生物酶和激素等物质
[6,11,17],但有研究认为二者对病原菌的抑制作用本质上是生物因素引起,而不是由化学物质主导
[2]。因此,为探明蚯蚓肥抑菌作用机制,本研究拟以华北平原牛粪蚯蚓堆肥为研究对象,采用微生物高通量测序方法,分析蚯蚓肥细菌、真菌群落结构和其代谢功能,并运用平板对峙方法,以灭菌的蚯蚓肥浸提液为对照,对引起作物腐烂病的茄病镰刀菌进行抑菌试验,以此来验证蚯蚓肥中的微生物是否发挥抑菌作用。以期为实现蚯蚓肥资源高值化利用和减少化学杀菌剂的投入提供依据。
1 材料与方法
1.1 蚯蚓堆肥过程及样品采集
蚯蚓堆肥地点为河北省石家庄市行唐县兴农蚯蚓养殖基地,以牛粪作为主要苗床养殖基质,同时添加适量草炭、碎木屑和碎秸秆并混合均匀,调整碳氮质量比,促进蚯蚓繁殖。所用的蚓种为赤子爱胜蚓(Eisenia foetida),品种是大平2号,该品种个头小、产卵数量多,从卵到成虫需要20~40 d,接种一次就可进行大量繁殖。采用林下覆膜繁殖方法,苗床宽×长×高=1.8 m×10.0 m×0.8 m。通过喷淋措施,保持养殖蚯蚓的苗床湿度控制在40%~60%,温度25 ℃,此时蚯蚓生长迅速,繁殖能力最强。养殖时间4—10月,通过机械将蚯蚓过筛后,剩余基质经过一周左右晾晒灭菌,即得优质蚯蚓肥。供试蚯蚓肥按照m(水)∶m(土)=35%配制。
由于赤子爱胜蚓按照生态型分类属于表层种
[18],即主要消耗表层的有机物。因此,为保证样品的代表性,本研究选择已经完全腐熟的表层蚯蚓肥为研究对象,蚯蚓肥采集时避开表面杂质,去除牛粪、石块以及其他杂物后,用土钻随机采样5次并混合均匀,取样深度约40 cm,将土钻预先经过消毒灭菌处理,将蚯蚓肥过10目土壤筛后,一部分自然晾干,待测定蚯蚓肥pH及过200目筛后测定元素含量;另一部分分装到无菌离心管里,每管3~5 g,密封后立即-80 ℃储存备用。
1.2 蚯蚓肥中养分测定
土壤pH采用电位法测定(
V(水)∶
m(土)=2.5∶1),采用重铬酸钾氧化加热法测定蚯蚓肥中有机质含量(OM)。采用凯氏定氮法和FOSS凯氏8200(FOSS, Hilleroed, Denmark)分析总氮(TN)含量
[19]。蚯蚓肥样品在HNO
3-HClO
4-HF溶液中用高压消解容器消解后,使用电感耦合等离子体-光学发射光谱仪(ICP-OES)测定总磷(TP)、总钾(TK)、总铅(Pb)、总铬(Cr)和总砷(As)
[20]。
1.3 蚯蚓肥微生物高通量测序
1.3.1 DNA提取和PCR扩增
采用CTAB法提取蚯蚓肥样本中DNA,利用1%琼脂糖凝胶电泳检测DNA的纯度和浓度,取适量的蚯蚓肥于离心管中,使用无菌水将蚯蚓肥稀释至1 ng/µL,用稀释后的基因组DNA作为模板,使用高效和高保真酶、Phusion® High-Fidelity PCR Master Mix with GC Buffer及带Barcode的特异引物341F(CCTAYGGGRBGCASCAG)和806R(GGACTACNNGGGTATCTAAT)对V3+V4可变区进行细菌PCR扩增,蚯蚓肥中真菌扩增区域为ITS2区,ITS3-KYO2(GATGAAGAACGYAGY RAA)和ITS4(TCCTCCGCTTATTGATATGC)。
1.3.2 PCR产物的混样与纯化
采用Illumina NovaSeq对测序样本进行测序。根据PCR产物的浓度进行等浓度混样,在充分混匀后,用1×TAE浓度为2%的琼脂糖凝胶电泳纯化PCR产物,对目的条带使用胶回收试剂盒(Qiagen公司, Germany)回收产物。文库的构建通过使用Illumina公司(San Diego, CA, USA)TruSeq
® DNA PCR-Free Sample Preparation Kit建库试剂盒进行,经Qubit定量和文库检测,合格后,使用NovaSeq 6000 PE250进行上机测序。
图1展示了蚯蚓肥中细菌和真菌Alpha多样性稀释曲线。可见:稀释曲线随着测序量的增加逐渐趋于平缓,表明对蚯蚓肥的测序深度已经接近饱和,足以反映当前样本的物种多样性。
1.3.3 测序数据处理
测序完成后,采用QIIME软件进行生物信息学分析。测序数据经过barcode拆分后,运用QIIME2插件中的DADA2插件对所有样品的全部原始序列(input)进行质量控制(filtered),去噪(denoised),拼接(merged),并且去嵌合体(non-chimeric),形成OTU
[21]。选取OTU的代表性序列,与数据库(16S 默认Greengenes Database 13_8版本,ITS默认UNITE数据库8.2版本)进行比对从而获得物种注释信息
[22]。使用PICRUSt2和FAPROTAX对蚯蚓肥细菌群落进行功能预测,利用QIIME获得的colsed OTU Table与KEGG数据库进行比对,获得功能预测信息。通过FUNGuild工具对真菌营养类型及其功能类群分类。使用微科盟生科云平台(
https://www.bioincloud.tech/)进行绘图。
1.4 蚯蚓肥浸提液提取及病原菌培养
蚯蚓肥浸提液制作方法主要参考
[23-24],称取20 g风干粉碎的蚯蚓肥样品于250 mL锥形瓶中,加入100 mL无菌去离子水,密封好后置于常温摇床内振荡72 h。经过4层无菌纱布和真空漏斗过滤使其不含颗粒物质后转移到锥形瓶内,即为本试验用蚯蚓肥浸提液。采用平板对峙的方法,研究蚯蚓肥浸提液对引起作物腐烂的主要病原菌腐皮镰刀菌(
Fusarium solani)的抑制作用
[25]。病原菌在马铃薯葡萄糖琼脂(PDA)上4 ℃保存至使用。浸提液的提取方法如下:在无菌锥形烧瓶(250 mL)中加入蚯蚓肥(10 g)和无菌反渗透水(100 mL)。为了避免污染,摇瓶盖上无菌密封膜,在定轨摇床(ZWY-240, Zhicheng, Shanghai, China)中以180 r/m摇动12 h。浸提液以4200 r/m(Eppendorf AG, Hamburg, Germany)离心10 min,收集上清液并分成两部分。其中一部分在使用前放置在4℃的冰箱中,其他部分立即用高压灭菌锅(MLS-3750; SANYO, Osaka, Japan)在121 ℃下灭菌20 min。用移液枪将灭菌和未灭菌的浸提液(100 μL)加入培养皿中,并用溅射镀膜仪均匀涂布。然后,在病原体菌落的边缘用5 mm直径的打孔机获得菌丝块
[26],并用灭菌镊子将菌丝块放在培养皿的中央。最后,根据菌丝生长情况观察蚯蚓肥浸提液抑菌活性。
2 结果与分析
2.1 蚯蚓肥中养分特点
经过土壤酸碱度、土壤养分和重金属含量测定发现:蚯蚓肥pH 7.98,有机质含量24.4%,全氮16.2%,全磷13.2%,全钾19.4%;总铅含量为21.6 mg/kg,总铬含量为68.5 mg/kg,总砷含量为8.3 mg/kg。说明按照本试验方法所得蚯蚓肥中养分丰富,能够满足作物生长,与NY 525—2012《有机肥料》
[27]重金属含量相比,没有超过规定的标准。
2.2 蚯蚓肥中微生物种群多样性
2.2.1 蚯蚓肥细菌菌群结构特征
根据16S rRNA扩增子测序结果,分析得出蚯蚓肥细菌菌群结构特征见
图2。可见:在门水平上,细菌群落中变形菌门Proteobacteria为相对丰度最大的优势菌门,占33%;群落中放线菌门Actinobacteria(17%)和拟杆菌门Bacteroidetes(13%)的相对丰度>10%,其次是绿弯菌门Chloroflexi(9%)、酸杆菌门Acidobacteria(9%)、厚壁菌门Firmicutes(7%)、芽单胞菌门Gemmatimonadetes(5%)(
图2(a));在目水平上,相对丰度>5%的有3个目,分别为噬纤维菌目Cytophagales(9%)、酸微菌目Acidimicrobiales(7%)、黄色单胞菌目Xanthomonadales(6%),未分类群共占15%(
图2(b));在属水平上平均丰度>1%的有凯氏杆菌属
Kaistobacter(2%)、氨氧化菌属
Candidatus_Nitrososphaera(2%)、乳酸杆菌属
Lactobacillus(1%)3个,未分类群占79%。
2.2.2 蚯蚓肥真菌菌群结构特征
根据ITS扩增子测序结果得出蚯蚓肥真菌菌群结构特征见
图3。可见:在门水平上,真菌群落中相对丰度最大的优势菌门为子囊菌门Ascomycota和被孢霉门Mortierellomycota,分别占75%和24%(
图3(a))。在目水平上,优势目为小子囊菌目Microascales,占44%,其次为被孢霉目Mortierellales和粪壳菌目Sordariales,分别占24%和18%。平均丰度>1%的有2个目,为肉座菌目Hypocreales(3%)和盘菌目Pezizales(1%)。未分类群平均丰度占10%(
图3(b))。在属水平上,优势属为被孢霉属
Mortierella,占24%。其中平均丰度>0.1%的有4个属,分别为镰刀菌属
Fusarium(1%)、盾盘菌属
Scutellinia(0.9%)、假阿利什霉属
Pseudallescheria(0.4%)和奈尼兹皮属
Nannizzia(0.1%)。未分类群平均丰度占72%(
图3(c))。
2.3 蚯蚓肥中微生物群落功能预测
基于KEGG pathway,运用PICRUSt2对蚯蚓肥细菌群落进行代谢功能预测,结果见
图4(a)。可见:在Pathway level 2中,代谢(74%)所占比例最高,其次是遗传信息处理(10%)。蚯蚓肥细菌群落在代谢水平中,氨基酸代谢(15%)、辅因子和维生素的代谢(13%)和碳水化合物代谢(13%)的功能丰度较其他功能最高。为进一步分析细菌群落的功能信息,本研究通过FAPROTAX对蚯蚓肥细菌菌群进行生态功能预测(
图4(b))。共获得63项功能分组,其中与化能异养、好氧化能异养和发酵作用相关基因信息所占比例最高,分别为11%、8%和3%。
图4(c)为所测样本中真菌群落功能注释结果。共统计分析了相对丰度较高的8类预测功能,其中:有4类被标记为腐生营养型、3类为病理营养型、1类为共生营养型;未定义腐生真菌(Undefined Saprotrophs)的相对丰度占比最高,土壤腐生真菌(Soil Saprotroph)和木腐生真菌(Wood Saprotroph)次之,动物病原真菌(Animal pathogen)、昆虫寄生虫(Insect parasite)和植物病原真菌(Plant pathogen)等的相对丰度最低。
2.4 蚯蚓肥浸提液对植物病原菌的抑制作用
本研究对易导致农作物腐烂的病原菌茄病镰刀菌进行了平板对峙试验,结果见
图5。可见:与
图5(a)正常培养的茄病镰刀菌相比,蚯蚓肥浸提液能够有效抑制土传病原真菌菌丝
F. solani的生长(
图5(c)~(f));随着蚯蚓肥浸提液的所含比例的增加,菌丝块生长逐渐受到抑制,200 μL处理下的浸提液抑菌效果最佳(
图5(f)),而灭菌后的蚯蚓肥浸提液抑菌效果(
图5(b))则明显减弱,说明蚯蚓肥浸提液中的抑菌效果与蚯蚓肥中有益微生物活性相关。
3 讨 论
从理化性质来分析,本研究测定的蚯蚓肥pH接近8,这可能与堆肥时间从4—10月,蚯蚓肥中微生物活动丰富、螯合物的形成以及游离氨的释放有关
[28-29],这与赵懿等
[30]以牛粪为饵料,蚯蚓堆肥达到45 d后pH超过8结果一致。对有机肥而言,除了基本的理化性质外,重金属含量是人们最为关注的指标。NY 525—2012《有机肥料》
[27]规定的总铅、总铬和总砷含量不能超过50、150和15 mg/kg,本研究中的蚯蚓肥这3种重金属含量分别为21.6、68.5和8.3 mg/kg,均未超过此标准,说明本研究所用蚯蚓肥比较安全。
已有研究表明蚯蚓堆肥中使用污泥、厨房、植物和动物废弃物等原料会影响蚯蚓肥中细菌群落的丰度和多样性
[31-33]。宋晓贺等
[34]从蚯蚓肥中共检测到34个菌门,其中变形菌门相对丰度占比最大。蔬菜、牛粪和污泥3种蚯蚓肥中丰度最高的同样为变形菌门,分别占48%、31%和36%
[35]。上述细菌群落的高通量测序结果与本研究结果一致,即变形菌门是蚯蚓肥细菌群落中的主要菌门。变形菌门在营养丰富的环境中被刺激生长,并且能够降解复杂的有机化合物
[36],在促植物生长以及土壤中的碳、氮和硫循环中,具有重要作用
[37-39]。对真菌群落进行高通量测序,本研究中子囊菌门相对丰度最大,而牛粪和白葡萄渣为原料进行蚯蚓堆肥
[40]的测序结果同样显示子囊菌门为优势菌门。因此,尽管堆肥原料不一致,但蚯蚓肥细菌和真菌菌群中的优势菌群不会有显著改变。
蚯蚓肥作为一种质优、环境友好的多功能型天然有机肥,在促进作物生长方面发挥着积极作用。Mu等
[41]从蚓肥中筛选出具有广谱抗真菌能力的枯草芽孢杆菌
Bacillus subtilis M29,该菌释放具有挥发性的有机化合物,其能对灰霉病菌
Botrytis cinerea发挥显著的拮抗作用。真菌群落中,毛壳菌属
Chaetomium如球毛壳
Chaetomium globosum和角毛壳
Chaetomium cupreum的一些特定品系能抑制植物病原菌、促进植物生长和诱导植物产生抗性等
[42]。其次,被孢霉属的成员可产生花生四烯酸,是一种植物病原真菌的抑制剂
[43]。此外,蚯蚓肥中的微生物还能促进养分转化,间接调控土壤微生物群落结构。如土壤中的大量磷素约有95%不能被植物直接吸收,可以利用蚯蚓肥中的溶磷菌来研制菌肥以代替传统化肥,提高土壤中可溶性磷的含量,改良土质
[44]。
本研究基于KEGG数据库,应用PICRUSt2预测蚯蚓肥中细菌微生物的代谢功能。代谢组所占比例明显高于其他途径,其次相对丰度较高的是氨基酸代谢,而传统堆肥系统中相对丰度最高的则是碳代谢
[45]。表明细菌的氨基酸代谢群可以增加氨基酸和腐殖酸的含量,氨基酸在堆肥过程中可充当微生物生长和代谢的能量和碳源。此外,碳水化合物代谢可通过降解纤维素和半纤维素提供各种基质,意味着碳相关代谢和有机物质的生物转化具有高活性
[45-46]。本研究对蚯蚓肥真菌群落功能进行预测分析发现,其主要以腐生营养型和病理营养型为主,未定义的腐生菌在序列丰度中占主导地位,该结果与Nguyen等
[47]的研究结果相似。可能是由于蚯蚓肥中微生物种类丰富,本研究中采用扩增子技术挖掘微生物群落,对物种的注释程度存在一定局限性,未来可以继续采用宏基因组技术深入挖掘有益微生物种类及其功能。
病原真菌腐皮镰刀菌(
F. solani)是黄连、黄芪、枸杞、半夏等药材根腐病的主要致病菌之一,目前对其防治手段仍是以化学防治为主
[48]。已有研究表明蚯蚓肥能够显著抑制腐皮镰刀菌
[49]。本研究中,蚯蚓肥浸提液能够对病原真菌的生长起到抑制作用,而其灭菌后的抑制作用减弱,表明可能是蚯蚓肥浸提液中聚集的潜在活性微生物起到了抑制致病菌的效果,结果与Yatoo等
[2]研究发现蚯蚓肥浸提液防治植物病害的本质是生物因素的观点一致。与蚯蚓肥相比,蚯蚓肥浸提液中微生物的多样性不仅没有降低,而且对病原菌有抑制作用的假单胞菌(
Pseudomonads)和梭状芽胞杆菌(
Clostridia)的相对丰度则明显增加
[50]。蚯蚓肥浸提液所含毛霉菌门和变形菌门显示出高比例,而毛霉菌在促进有机物分解方面发挥关键作用
[51],说明蚯蚓肥浸提液中不仅存在防治病害的微生物,而且这类微生物数量较蚯蚓肥中更加丰富。蚯蚓肥浸提液中的矿物质、腐殖酸含量、pH和电导率与蚯蚓堆肥中已知的含量非常相似
[52],蚯蚓肥浸提液代替固体蚯蚓肥,很大程度上减少了所需固体蚯蚓肥的用量
[53]。本研究结果进一步证明了蚯蚓肥中微生物在生物防治方面的开发潜力。
4 结 论
本研究以华北平原牛粪蚯蚓堆肥为研究对象,探究蚯蚓肥中养分、细菌和真菌微生物群落结构及其抑菌功能,主要结论如下:1)蚯蚓肥中富含大量有机质和氮磷钾等,总铅、总铬和总砷的重金属含量远低于现行有机肥标准;2)蚯蚓肥中变形菌门和放线菌门是细菌群落的优势菌门,真菌群落中优势菌门为子囊菌门和被孢霉门;3)蚯蚓肥中细菌群落以氨基酸代谢功能为主,而真菌群落以腐生营养型较多,能够有效促进腐殖质分解。与对照相比,灭菌后蚯蚓肥浸提液抑制茄病镰刀菌的活性减弱,进一步说明蚯蚓肥对病原菌的抑制作用是由微生物本身主导。因此,蚯蚓肥可以作为一种潜在的实现作物养分供应和病害防治双赢的优质有机肥料,在生态农业中具有广阔应用前景。
河北省自然科学基金资助项目(H2022423004)
河北省中医药管理局科研计划项目(2022098)
河北省高等学校科学技术研究项目(QN2023026)
河北中医药大学博士科研基金项目(BSZ2021019)
河北省现代农业产业技术体系创新团队项目(HBCT2023080201)
河北省现代农业产业技术体系创新团队项目(HBCT2023080205)