土壤微生物是土壤生态系统中的重要组成部分,直接或间接参与土壤中的众多生化反应,对维持土壤生态平衡、推动土壤养分循环以及保障作物生长发育具有重要意义
[1-2]。土壤微生物对于土壤微环境具有重要影响,直接参与土壤多种养分的转化和土壤理化性质的形成,对土壤养分循环利用和土壤微环境稳定性发挥着重要作用
[3]。调节土壤微生物群落结构是提升土壤健康和作物产量品质的重要手段,其中有机肥和生物肥料的应用对土壤微生物群落结构和功能具有显著影响
[4-5]。有机肥含大量有机质,能够有效增加土壤养分,调节微生物群落结构,增强土壤微生物代谢活动
[6-7]。目前,微生物菌剂在土壤改良及优化土壤性能方面的作用日益凸显
[8]。微生物菌剂不仅可以增加土壤微生物的种群多样性,还能使其中的有益微生物与根际微生物协调共生,维持土壤生态群落稳定,平衡微生物种群
[9-10]。施河丽等
[11]研究发现,将微生物菌剂和黄腐酸联合施用,能显著提高烟草根际土壤中细菌种群的多样性。尽管目前已有较多研究关注有机肥或微生物菌剂单一应用对土壤微生物群落的影响
[12-15],但联合应用对植烟土壤微生物群落的研究相对较少。本课题组前期研究发现,施用芝麻饼肥可有效改善土壤理化性质,通过提高根系生长和香气前体物质的积累提高烟叶品质
[16]。因此,本研究基于芝麻饼肥的应用,利用高通量测序技术,研究不同微生物菌剂配施对植烟土壤微生物及理化特性的影响,以期为烟草生产中芝麻饼肥与适宜微生物菌剂的配合使用提供参考。
1 材料与方法
1.1 试验地概况
试验于2024年在广东省韶关市黄坑镇(114°30′21″ E,25°14′44″ N)进行。供试烤烟品种为‘粤烟97’。试验田前茬作物为玉米,土壤类型为碱性紫色土,土壤pH值 8.26,有机质含量17.12 g/kg,全氮含量0.64 g/kg,全磷含量1.24 g/kg,全钾含量16.96 g/kg,碱解氮含量63.83 mg/kg,有效磷含量6.08 mg/kg,速效钾含量181.48 mg/kg。
1.2 供试微生物菌剂
拜赖青霉菌(北海业盛旺生物科技有限公司):有效活菌数≥1×109/g,含有拜赖青霉分生孢子、发酵代谢产物和固体发酵基质等;枯草芽孢杆菌(广西康绿生物科技有限公司):有效活菌数≥4×1010/g,含有枯草芽孢杆菌等;巨大芽孢杆菌(北海业盛旺生物科技有限公司):有效活菌数≥1×1010/g,含有巨大芽孢杆菌及其代谢产物、适宜载体等;吲哚丁酸钾(湖北清宇化工有限公司):吲哚丁酸钾含量≥98%;复合型微生物菌剂(山西永济市华杰生物科技有限公司):有效活菌数≥2×109/g,含有枯草芽孢杆菌、巨大芽孢杆菌和解淀粉芽孢杆菌,高活性矿源黄腐酸钾≥50%,腐殖酸≥60%,有机质≥65%,有机碳≥45%,氧化钾≥12%。
1.3 试验设计
试验共设置6个处理,分别为芝麻饼肥+清水(CK)、芝麻饼肥+3.0 kg/hm2拜赖青霉菌(T1)、芝麻饼肥+1.5 kg/hm2枯草芽孢杆菌(T2)、芝麻饼肥+7.5 kg/hm2巨大芽孢杆菌(T3)、芝麻饼肥+7.5 kg/hm2吲哚丁酸钾(T4)和芝麻饼肥+3.0 kg/hm2复合型微生物菌剂(T5)。每个处理设置3次重复,共18个小区,小区面积为39.6 m2(11 m×3.6 m)。供试烤烟品种为‘粤烟97’,种植密度为15 000株/hm2,行株距为120 cm×55 cm。微生物菌剂于移栽当天进行灌根,灌根质量浓度均为0.1 mg/L,将微生物菌剂溶于清水,每穴灌溶液500 mL,每隔15 d灌根1次,共灌根4次。采用漂浮育苗,2024年2月19日所有处理均施用750 kg/hm2芝麻饼肥(m(N)∶m(P2O5)∶m(K2O)=3.73∶2.99∶1.50,驻马店市康博汇鑫农业科技有限公司)和330.60 kg/hm2烟草专用复合肥(m(N)∶m(P2O5)∶m(K2O)=13∶8∶15,湖南金叶众望科技有限公司)作为基肥,并在当天进行膜下移栽。基肥追肥比例为5∶5,基肥采用穴施,追肥采用淋施。2024年3月24日所有处理均施用288.45 kg/hm2烟草专用复合肥作为第1次追肥,2024年4月2日所有处理均施用55.95 kg/hm2烟草专用复合肥和139.35 kg/hm2硝酸钾(m(N)∶m(P2O5)∶m(K2O)=12∶0∶45,云南云叶化肥股份有限公司)作为第2次追肥,2024年4月12日所有处理均施用79.65 kg/hm2硝酸钾和205.05 kg/hm2农业用硫酸钾(m(N)∶m(P2O5)∶m(K2O)=0∶0∶52,国投新疆罗布泊钾盐有限责任公司)作为第3次追肥。田间管理按当地优质烟叶生产技术规范进行。
1.4 样品采集与处理
于旺长期(移栽后65 d)随机选取烟垄上2株烟苗正中位置,采集0~20 cm土壤,分为两份,一份自然风干后过筛(孔径0.18 mm),用于测定土壤基本理化性质,另一份在4 °C保存,用于测定土壤酶活性。同时,每个小区选取长势一致的烟株3株,抖根法采集根际土壤,混合后作为一个样品,剔除植物体等杂质后装入无菌管并迅速用液氮速冻后,于-80 °C保存待测,用于DNA高通量测序分析。
1.5 土壤理化指标测定
土壤基本理化性质参照《土壤农化分析》
[17]的方法测定:pH值采用pH计法测量(水土质量比为2.5∶1);有机质含量采用重铬酸钾容量法测定;全氮含量采用全自动凯氏定氮仪Kjeltec
TM 8400(FOSS,丹麦)测定;全磷含量采用HClO
4-H
2SO
4法测定;全钾含量采用氢氧化钠熔融法测定;碱解氮含量采用碱解扩散法测定;有效磷含量采用0.5 mol/L NaHCO
3浸提,钼锑抗比色法测定;速效钾含量采用1 mol/L NH
4OAc浸提,火焰光度法测定。土壤体积质量采用环刀法测定;土壤含水率采用烘干法称重测定。各处理重复测定3次。
土壤脲酶、蔗糖酶、过氧化氢酶和酸性磷酸酶活性使用试剂盒(上海晶抗生物科技有限公司)测定。各处理重复测定3次。
1.6 高通量测序分析
根据E.Z.N.A.
®soil DNA kit (Omega Bio-tek,Norcross,GA,U.S.)说明书进行微生物群落总基因组 DNA 抽提, 1%琼脂糖凝胶电泳检测抽提的基因组DNA质量,NanoDrop2000(美国Thermo Scientific公司)测定DNA 浓度和纯度。使用携带Barcode序列的上游引物338F (5′-ACTCCTACGGGAGGCAGCAG-3′)和下游引物806R (5′-GGACTACHVGGGTWTCTAAT-3′)
[18]对16S rRNA基因V3-V4可变区进行 PCR 扩增,使用2%琼脂糖凝胶回收PCR产物,利用DNA凝胶回收纯化试剂盒(PCR Clean-Up Kit,中国逾华)进行回收产物纯化,并用Qubit 4.0(Thermo Fisher Scientific,USA)对回收产物进行检测定量,使用NEXTFLEX Rapid DNA-Seq Kit对纯化后的PCR产物进行建库,利用Illumina Nextseq2000平台进行测序(上海美吉生物医药科技有限公司)。使用fastp
[19](
https://github.com/OpenGene/fastp,version 0.19.6)软件对双端原始测序序列进行质控,用FLASH
[20](
http://www.cbcb.umd.edu/software/flash,version 1.2.11)软件进行拼接;基于默认参数,使用Qiime2流程
[21]中的DADA2
[22]插件对质控拼接之后的优化序列进行降噪处理,得到ASVs(扩增子序列变体,amplicon sequence variant)。基于Sliva 16S rRNA基因数据库(v138),使用Qiime2中的Naive bayes分类器对ASVs进行物种分类学分析。使用PICRUSt2
[23](version 2.2.0)软件进行16S功能预测分析。
1.7 数据处理
采用Excel 2021、SPSS 26.0、Origin 2021软件进行数据处理和单因素方差分析,用Duncan新复极差法进行多重比较;在美吉生物云平台(
https://cloud.majorbio.com)上进行土壤微生物数据分析。
2 结果与分析
2.1 芝麻饼肥配施不同微生物菌剂对植烟土壤理化性质的影响
2.1.1 植烟土壤养分和物理特性
由
表1可知,芝麻饼肥配施不同微生物菌剂可有效提高植烟土壤有机质、全氮、全钾、碱解氮、有效磷、速效钾含量,并降低土壤体积质量,提高土壤含水率。与CK相比,T3和T5处理对土壤养分提升效果最显著,有机质含量分别显著提高了43.83%和61.66%,全钾含量分别显著提高了20.86%和26.44%,碱解氮含量分别显著提高了19.74%和32.89%,有效磷含量分别显著提高了23.59%和40.40%,速效钾含量分别显著提高了10.65%和12.60%。此外,T2和T5处理能有效改善土壤物理性状,与CK相比,T2和T5处理土壤体积质量分别降低了11.88%和8.91%,而含水率分别提高了40.91%和73.79%。
2.1.2 植烟土壤酶活性
由
表2可知,芝麻饼肥配施微生物菌剂均能显著提高植烟土壤脲酶(除T4处理)、蔗糖酶、过氧化氢酶、酸性磷酸酶活性。与CK相比,T3和T5处理土壤脲酶活性分别显著提高了66.82%和64.67%,T2和T5处理土壤蔗糖酶活性分别显著提高了104.20%和107.87%,T1和T5处理土壤酸性磷酸酶活性分别显著提高了126.67%和120.0%,T2和T5处理土壤过氧化氢酶活性分别显著提高了66.11%和65.11%。
2.2 芝麻饼肥配施不同微生物菌剂对植烟土壤细菌群落的影响
2.2.1 植烟土壤细菌群落 Alpha 多样性
对18个土壤样本(6个处理×3次重复)进行测序,共得到18 371个ASV,归属于49门、136纲、311目、518科、1 079属和2 296种。对不同处理的土壤细菌群落ASV进行Alpha多样性分析,分别用Ace指数(丰富度指数)、Chao指数(丰富度指数)、Shannon指数(多样性指数)、Simpson指数(多样性指数)4个度量指标探究各处理内细菌群落结构特征,结果见
表3。由
表3可知,T2、T3和T5处理的Ace指数、Chao指数与CK相比均显著升高,其中T2处理增幅最大,分别显著提高了15.58%和14.96%。施用微生物菌剂处理的土壤细菌Shannon指数、Simpson指数相较于CK均有所提高,但仅T2处理Shannon指数的差异达到显著水平。
2.2.2 植烟土壤细菌群落组成
由
图1可知,6组样品中注释到共有的细菌属455属,占属总数的42.17%,CK、T1、T2、T3、T4和T5处理的特有细菌属分别为50,33,56,39,45和31属,分别占属总数的4.63%,3.06%,5.19%,3.61%,4.17%和2.87%。
由
图2-A可以看出,各处理土壤细菌群落中,假单胞菌门(Pseudomonadota,相对丰度43.63%~48.97%)和放线菌门(Actinomycetota,相对丰度15.00%~20.91%)为最丰富菌门,其次是芽孢杆菌门(Bacillota,相对丰度8.17%~13.15%)、拟杆菌门(Bacteroidota,相对丰度5.22%~12.15%)。与CK相比,T1、T2、T3和T5处理降低了假单胞菌门、放线菌门的相对丰度,而增加了芽孢杆菌门、拟杆菌门的相对丰度;T4处理降低了放线菌门、拟杆菌门的相对丰度,而增加了假单胞菌门、芽孢杆菌门的相对丰度。
在进行优势物种分析时,相对丰度超过1%即可作为优势物种。由
图2-B可以看出,CK的优势菌属主要有鞘氨醇单胞菌属(
Sphingomonas,相对丰度9.42%,下同)、微枝形杆菌属(
Microvirga,3.72%)、溶杆菌属(
Lysobacter,2.71%)等;T1处理的优势菌属与CK大体一致。相较于CK,T2处理的优势菌属主要有硝化螺旋菌属(
Nitrospira,3.95%)、
MND1菌属(3.46%)、红色杆菌属(
Rubrobacter,3.01%)等,T3处理的优势菌属主要有硝化螺旋菌属(3.91%)、
MND1菌属(2.61%),T4处理的优势菌属主要有鞘氨醇单胞菌属(11.31%)、硝化螺旋菌属(3.42%)、
MND1菌属(2.92%),T5处理的优势菌属主要有微枝形杆菌属(5.56%)、硝化螺旋菌属(2.80%)、
MND1菌属(2.69%)。由此可见,配施微生物菌剂可以增加土壤中优势细菌属的相对丰度。
由
图3可以看出,不同处理土壤优势细菌在属水平上的多样性和丰度差异较大,施用微生物菌剂处理的硝化螺旋菌属、
MND1菌属均高于CK。其中,红色杆菌属、类固醇杆菌属(
Steroidobacter)在T2处理中的占比高于其他处理,鞘氨醇单胞菌属在T4处理中的占比最高,微枝形杆菌属在T5处理中的占比高于其他处理。
选取相对丰度排名前20的细菌属进行分析,结果(
图4)显示,与CK相比,配施不同微生物菌剂均提高了剑菌属(
Ensifer)、
Peribacillus、黄杆菌属(
Flavobacterium)、
Bryobacter、微枝形杆菌属、
Ramlibacter等菌属的相对丰度。
2.2.3 LEfSe 多级物种差异判别分析
LEfSe 是一种用于发现高维生物标识和揭示基因组特征的软件。LEfSe采用线性判别分析(LDA)估算每个组分(物种)丰度对差异效果影响的大小。从门水平到属水平,各个分类单位设定LDA阈值为2。由图
5和
6可知,CK处理放线菌门、溶杆菌目(Lysobacterales)显著富集,T1处理溶杆菌科(Lysobacteraceae)、脱卤球菌纲(Dehalococcoidia)、壤霉菌属(
Agromyces)显著富集,T2处理玫瑰弯菌科(Roseiflexaceae)、酸杆菌纲(Acidobacteriae)显著富集,T3处理黄杆菌目(Flavobacteriales)、黄杆菌属显著富集,T4处理鞘氨醇单胞菌属、芽单胞菌门(Gemmatimonadota)显著富集,T5处理黄色土源菌属(
Flavisolibacter)、房间芽孢杆菌属(
Domibacillus)显著富集。
2.2.4 植烟土壤细菌群落 NMDS分析
基于 Bray-Curtis 距离对 6个处理的土壤细菌群落组成进行非度量多维尺度分析(non-metric multidimensional scaling,NMDS)分析,结果见
图7。由
图7可以看出,T1、T2、T3、T4处理细菌群落组成均与CK、T5处理存在明显差异,且CK与T5处理的组成也有较大差异,说明施用微生物菌剂改变了土壤细菌群落组成;此外,T2和T3处理细菌群落组成较为相似,与T1、T4处理有明显差异,说明不同微生物菌剂的效果有所不同。
2.2.5 功能预测
通过PICRUSt2软件对细菌 ASV信息进行功能预测,COG功能分析结果见
图8。由
图8可以看出,不同处理土壤相对丰度大于5%的功能基因包括翻译核糖体结构和生物发生(translation, ribosomal structure and biogenesis)、转录(transcription)、细胞壁/膜/包膜生物发生(cell wall/membrane/envelope biogenesis)、能量生产与转化(energy production and conversion)、碳水化合物的运输与代谢(carbohydrate transport and metabolism)、氨基酸转运与代谢(amino acid transport and metabolism)、无机物运输和代谢(inorganic ion transport and metabolism),各处理之间关于细菌功能的差别不明显。
2.3 土壤细菌群落组成与土壤理化性质和酶活性的相关分析
由
表4可以看出,土壤理化性质显著影响了土壤酶活性和细菌群落组成。土壤脲酶、蔗糖酶、酸性磷酸酶、过氧化氢酶活性与有机质、碱解氮、有效磷含量呈显著或极显著正相关,土壤pH与蔗糖酶活性呈极显著正相关,土壤体积质量与脲酶、酸性磷酸酶、过氧化氢酶活性呈显著正相关,土壤含水率与蔗糖酶、过氧化氢酶活性呈极显著正相关,而土壤全磷含量与酸性磷酸酶活性呈极显著负相关;土壤有机质含量和脲酶、蔗糖酶、酸性磷酸酶活性与Ace指数、Chao指数、Shannon指数呈显著或极显著正相关,土壤有效磷含量、含水率、过氧化氢酶活性与Ace指数、Chao指数呈显著或极显著正相关,土壤有效磷含量与Simpson指数呈显著正相关,而土壤体积质量与Shannon指数呈显著负相关。
上述结果表明,土壤环境因子变化对土壤酶活性及土壤细菌群落的丰富度和多样性有一定影响,其中土壤有机质、碱解氮和有效磷是影响土壤酶活性的主要环境因子,土壤有机质和有效磷是影响土壤细菌群落丰富度和多样性的主要环境因子,土壤脲酶、蔗糖酶和酸性磷酸酶对土壤细菌群落丰富度和多样性也具有重要影响。
3 讨 论
施用有机肥或微生物菌剂是改善植烟土壤质量和提升烟叶品质的重要手段
[24-25]。本研究结果表明,芝麻饼肥配施不同微生物菌剂均有效提高了植烟土壤有机质、全氮、全钾、碱解氮、有效磷和速效钾含量及土壤含水率,同时降低了土壤体积质量。这可能是因为微生物菌剂中的有效菌分泌多糖等黏性物质,促进土壤团粒结构的形成,增加土壤的通气性、透水性和保水性
[26],从而降低了土壤体积质量。此外,微生物菌剂还能分解土壤中的有机质,将其转化为植物可吸收的无机营养物质(如氮、磷、钾等),促进土壤中难溶性养分的溶解和释放,提高土壤的养分供应能力,最终改善植烟土壤的养分状况
[27]。赵津聪等
[28]研究也表明,微生物菌剂和有机肥配施能提高葡萄根际土壤的有机质、碱解氮、有效磷和速效钾含量。侯会静等
[29]研究发现,功能性微生物进入土壤后可有效活化被固定的磷、钾等养分,从而提升土壤肥力。有机肥为微生物菌剂提供了良好的生长环境和碳源,促进了土壤中微生物数量的增加;而微生物数量的增多又加速了有机质的分解,并推动了土壤难溶性养分的转化过程
[30-31]。
土壤酶作为土壤生态系统中的关键生物催化剂,参与有机质分解和养分转化等重要过程,其活性水平能够反映土壤生物化学反应的强度和方向
[32-33]。本研究发现,施用微生物菌剂(除T4处理)的土壤酶活性均显著提升,且土壤脲酶、蔗糖酶、酸性磷酸酶和过氧化氢酶的活性均与土壤有机质、碱解氮和有效磷含量呈显著或极显著正相关。这可能是由于微生物菌剂的施用促进了土壤和肥料中养分的释放,满足了土壤微生物对养分的需求,为土壤酶的分泌创造了有利条件,从而促进了土壤酶活性的提高
[34]。土壤脲酶活性与土壤氮循环密切相关,是衡量土壤供氮能力的重要指标
[35];蔗糖酶参与高分子有机物的分解,通过水解蔗糖生成葡萄糖和果糖,补充碳源和能量,可反映土壤的熟化程度和肥力水平
[36];磷酸酶的活性直接影响土壤中有机磷的分解转化及其生物有效性,是评价土壤磷素营养状况的关键指标;过氧化氢酶则通过分解过氧化氢,减轻其对生物体的毒害作用,其活性与土壤有机质含量和微生物数量密切相关
[37]。在本试验中,T5(复合型微生物菌剂)处理的土壤酶活性提升最显著,表明复合型微生物菌剂对植烟土壤的供氮能力及营养元素转化具有显著促进作用,这一结果与前人研究结果
[25, 38-40]一致。
土壤微生物是土壤微生态的重要组成部分,对土壤的理化性质、养分循环和肥力水平产生重要影响,进而维持土壤质量并促进植物生长
[12]。微生物菌剂的施用能够调节土壤微生物群落结构,尤其对土壤细菌群落产生显著影响
[41]。本研究发现,芝麻饼肥与枯草芽孢杆菌(T2)配施时,土壤细菌的丰富度和多样性最高。赵诣等
[42]研究也表明,枯草芽孢杆菌的配施能够显著提升土壤细菌群落的丰富度和多样性,并优化其结构。通过对不同处理土壤细菌群落分析发现,施用微生物菌剂后芽孢杆菌门的相对丰度有所增加,芽孢杆菌通过分解有机物质和固氮作用,参与土壤中的碳循环和氮循环,从而为土壤提供养分
[43]。本研究发现,施用微生物菌剂处理(除T1处理外)的优势菌属均有硝化螺旋菌属和
MND1菌属。研究发现硝化螺旋菌属通过促进土壤氮循环来促进植物的生长
[44],
MND1菌属通过改善土壤的养分循环过程间接促进植物生长
[45]。本研究发现,T2、T4和T5处理优势菌属还有红色杆菌属、鞘氨醇单胞菌属和微枝形杆菌属。研究发现红色杆菌属和微枝形杆菌属表现出对高温和干旱环境的普遍适应性
[46],鞘氨醇单胞菌属可以促进植物在干旱等胁迫条件下生长,提高植物抗逆性
[47]。此外,施用微生物菌剂处理的剑菌属、围孢杆菌属、黄杆菌属、布氏杆菌属的相对丰度均有所提高,且不同微生物菌剂处理显著富集的优势细菌种类存在差异,这表明芝麻饼肥与微生物菌剂的配施可以改变植烟土壤的优势细菌种类,增强微生物群落的复杂性和稳定性,芝麻饼肥配施枯草芽孢杆菌(T2)、吲哚丁酸钾(T4)和复合型微生物菌剂(T5)可以增强烟株对干旱环境的适应性。从土壤细菌NMDS分析图可以看出,T1、T2、T3、T4处理细菌群落结构相对于CK和T5处理有明显差异,而T5处理与CK的组成也有较大区别,这说明芝麻饼肥与微生物菌剂的配施能够改变土壤细菌群落结构,且不同微生物菌剂对细菌群落结构的影响不尽相同,这一结果与前人研究结果
[48-49]相一致。然而,芝麻饼肥配施微生物菌剂对植烟土壤细菌的功能类群影响不显著,推测功能类群的生态响应需要较长时间尺度的群落演替,短期试验可能不足以观察到显著的功能分化。
土壤有机质、碱解氮、有效磷和速效钾等关键养分指标是衡量土壤肥力水平的重要参数,对土壤微生物群落的构建和植物生长具有显著影响
[50-51]。本研究中,土壤有机质和有效磷含量的变化显著影响土壤细菌群落的丰富度,而土壤有机质、碱解氮和有效磷含量与土壤酶活性也呈显著正相关,这表明土壤细菌群落的多样性与土壤养分,尤其是速效养分密切相关。此外,土壤的物理特性与细菌群落的丰富度和多样性也表现出显著相关性,这可能与芽孢杆菌门对土壤有机质的分解作用进而改变土壤物理特性有关。土壤酶活性与细菌群落结构的变化也呈现显著相关性,这可能是因为土壤酶活性的增强促进了有机碳的分解和养分的循环,从而为细菌提供了更丰富的底物和更适宜的生存环境,从而丰富了细菌群落
[52]。
4 结 论
在施用芝麻饼肥基础上配施微生物菌剂可有效增加植烟土壤养分含量,改善土壤理化特性,提高土壤酶活性;土壤有机质、有效磷含量变化对土壤酶活性、细菌群落的影响较大。在芝麻饼肥应用基础上配施枯草芽孢杆菌(T2)、复合型微生物菌剂(T5)效果最优,可以有效增加植烟土壤有机质、有效磷含量,提高土壤酶活性,促进土壤细菌群落的繁殖,其主要通过增加芽孢杆菌门、拟杆菌门相对丰度提高土壤细菌丰富度和多样性,使植烟土壤微生物群落结构更加复杂和稳定。
广东中烟工业有限责任公司科研项目(2024440000340013)
广东省烟草专卖局科技计划资助项目(2023440000240044)
广东省基础与应用基础研究省市联合基金项目(2022A1515110353)