1. 引言
竹荪(Dictyophora),是一类营养丰富的大型野生食用菌,属担子菌亚门、腹菌纲、鬼笔目、鬼笔科、竹荪属,长裙竹荪是竹荪的主要栽培种之一 [1]。竹荪具有高蛋白、低脂肪、多氨基酸成分的特点,素有“菌中皇后”的美誉,经常食用对人类疾病的预防具有良好的效果 [2]。竹荪是名贵的食用菌,经济价值高,适合在荫蔽林下栽培,在我国福建、广西、云南等省区的林下栽培效益多有研究及报道 [3]。
目前,国内外有关竹荪栽培的研究主要集中在生长环境、栽培基质和连作障碍等方面的研究。相关研究表明,土壤环境影响竹荪的产量和品质 [4],栽培基质的成份及处理方式影响竹荪的菌丝生长 [5],连作影响竹荪栽培土壤微生物数量和质量,不同连作年限下的土壤酶活性差异显著 [6]。自2007年以来,海南省农垦系统和海南林业科学研究所等单位陆续开展了橡胶林下种植的长裙竹荪的试种和推广工作,增加了胶农收入,提高了橡胶林的综合经济效益 [7]。橡胶林下种植长裙竹荪产业既充分利用橡胶林木加工过程产生的“边脚料”资源,又充分利用橡胶林下“空闲”空间资源,实现生态环境和经济效益的双丰收 [8] [9]。但是,随着橡胶林下长裙竹荪种植年限的延长,相继出现产量降低、品质下降等连作障碍带来的典型现象。我省有关橡胶林下种植长裙竹荪的相关研究主要集中在竹荪间作对橡胶园土壤肥力、土壤微生物区系与群落功能多样性的影响,竹荪栽培种覆土层土壤微生物群落功能多样性等研究 [10] [11]。目前,有关海南橡胶林下种植长裙竹荪连作障碍成因的相关研究未见报道。本研究以橡胶林下不同连作年限长裙竹荪覆土层土壤为研究对象,分析土壤微生物、土壤酶活性及其土壤微生物功能多样性等指标特征,旨在探究海南橡胶林下种植长裙竹荪连作障碍的原因,为橡胶林下种植长裙竹荪连作障碍的解决提供基础,助力海南橡胶林下经济竹荪产业的健康发展。
2. 材料与方法
2.1. 样地概况及土壤样品采集
本研究样地位于海南省林业科学研究院(海南省红树林研究院)定安龙州基地内的橡胶林,该橡胶林林龄为15~20年,正常水肥管理和割胶。样地阳光充足,雨量充沛,年平均气温24℃,年平均降水量1953 mm。
以长裙竹荪“出菇”期作为采样时间,采集同一橡胶林地下长裙竹荪种植第一年、第二年(连作第1年)和第三年(连作第2年)的土壤,根据样地地形地貌,按“S”形随机选取5点,去除土壤表面凋落物,采集覆土层(0 cm~20 cm)土壤样品,除去大的石粒等杂物,混合均匀后装入无菌自封袋,放入恒温采样箱内并及时转运至实验室置4℃冰箱中保存。自样品采集24 h内完成样品的预处理、制备及分析工作。
2.2. 土壤微生物
2.2.1. 样品的预处理
取10 g土壤置于装有90 mL 0.85%无菌生理盐水的三角烧瓶中(内置5~10颗玻璃珠),120 r/min、37℃振荡30 min后备用。
2.2.2. 样品的梯度稀释
取样品悬浊液1 mL于9 mL无菌水中,充分混匀,梯度稀释至适宜的稀释度,取合适的三个梯度稀释液100 ul分别涂布于培养细菌、真菌、放线菌和固氮菌的培养基平板上,于对应的温度下培养,其中细菌35℃,真菌26℃,放线菌28℃),24 h~72 h后计数,采用最可能数法(MPN: most probable number)测定土壤样品中氨化细菌、好氧性纤维素分解细菌、亚硝化细菌、反硝化细菌等土壤中功能微生物的数量。
2.3. 土壤酶活性测定
按《土壤农业化学分析》相关方法分析橡胶林下长裙竹荪不同种植年限的覆土层土壤理化性质和酶活性,其中土壤有机质采用K2Cr2O7-H2SO4氧化外加热法,土壤全氮采用半微量凯氏法,土壤全磷采用HClO4-H2SO4法,土壤全钾采用HF-HClO4法,土壤pH值采用原位pH计测定;土壤蔗糖酶采用3,5–二硝基水杨酸比色法,土壤脲酶采用苯酚钠–次氯酸钠比色法,土壤脱氢酶采用容量法,土壤磷酸酶采用磷酸苯二钠比色法 [12]。
2.4. Biolog-ECO分析
2.4.1. 样品的与处理
另取10 g土壤于装有90 mL 0.85%无菌生理盐水的三角烧瓶中(内置5~10颗玻璃珠),120 r/min、37℃振荡30 min后,冰浴静置3 min~5 min。取静置的上清液梯度稀释至1:100作为样品测试液。
2.4.2. 样品的测定
取制备好的样品测试液加至Biolog-ECO板微孔板96孔中,150 ul/孔。28℃恒温避光培养9 d,每隔24 h用Biolog酶标仪测定样品在590 nm和750 nm的吸光值 [13]。
平均吸光值(AWCD值)来评估土壤微生物群落对碳源利用的整体情况及利用活性,其计算公式为:
。式中:Ci为所测定的31个碳源孔吸光值;R对对照孔的吸光值;n为培养基碳源总数,本研究为31。
Shananon-Wiener指数(H'值)来评估微生物物种的丰富度,Shannon-Wiener指数:
,其中
为有培养基的孔和对照孔的光密度值差与整板总差的比值 [14]。
Simpson指数(D值),又称优势度指数,用来评估物种的优势度。Simpson指数:
。
Pielou指数(E值),又称均匀度指数,是微生物群落实测多样性与最大多样性的比率。Pielou指数:
,式中S为被利用的碳源总数。
2.5. 数据处理及分析
本研究数据采用Excel 2010和SPSS 22.0进行统计分析,通过单因素方差分析(Analysis of Variance, ANOVA)来确定不同连作年限下长裙竹荪覆土层土壤理化性质、土壤酶活性、土壤可培养微生物总数、功能微生物及微生物群落功能多样性指数的显著水平。
3. 结果
3.1. 土壤理化性质及酶活性
3.1.1. 土壤理化性质
Table 1. Main physical and chemical properties of soil in different continuous cropping periods
表1. 不同连作时期土壤主要理化性质
注:同一列中,小写字母不同表示差异显著(P < 0.05),大写字母不同表示差异极显著(P < 0.01)。
从表1可知,橡胶林下种植长裙栽培竹荪土壤属于酸性土壤,pH值随长裙竹荪连作年限的增加而降低。不同连作期,各化学元素的含量有所变化。全氮(TN)和全磷(TP)的含量在种植第一年和种植第二年(连作第1年)有所增加,可能与长裙竹荪生长过程中对培养料的快速分解以及连作过程的养分积累有关。种植第二年(连作第1年)土壤的全N、全P含量与种植第三年(连作第2年)、种植前期均存在显著性差异。与种植前期比较,种植第三年(连作第2年)土壤中全N、全P的含量有所增加,可能与种植第二年(连作第1年)残存在土壤中的各类营养物质有关。种植前期至种植第三年(连作第2年),土壤中钾(K)的含量降低了0.67 g/kg。土壤中有机质含量在不同种植期有所不同,总体呈现出减少的趋势,部分连作年限间差异显著。
3.1.2. 土壤酶活性
从表2可以看出,种植长裙竹荪的橡胶林地土壤蔗糖酶、脲酶含量随长裙竹荪连作年限的增长而下降,种植前期显著高于其他种植期;土壤蔗糖酶含量,除种植第二年(连作第1年)和种植第三年(连作第2年)间差异不显著外,其他种植年限间差异均极显著;土壤脲酶含量,除种植第二年(连作第1年)和种植第1年间差异不显著外,其他种植年限间差异均极显著。
土壤磷酸酶含量随长裙竹荪连作年限的增长持续升高,种植第二年(连作第1年)最高,并且显著高于种植前期和种植第一年;种植第二年(连作第1年)和种植第三年(连作第2年)间差异不显著外,其他种植期间差异均极显著。土壤脱氢酶含量随长裙竹荪连作年限的增长呈现先降低后升高再降低的规律,种植第三年(连作第2年)显著高于其他种植期,各种植期之间差异均极显著。
注:同一列中,小写字母不同表示差异显著(P < 0.05),大写字母不同表示差异极显著(P < 0.01)。
3.2. 土壤中可培养微生物数量及功能微生物
3.2.1. 可培养微生物数量
从表3可以看出,橡胶林下长裙竹荪覆土层土壤中可培养的细菌、真菌和放线菌总数均随长裙竹荪连作年限的增长先升高后降低,并且在种植第二年(连作第1年)达到峰值。橡胶林下长裙竹荪种植第一年和种植第二年(连作第1年)土壤中可培养细菌和放线菌总数均差异不显著,但均极显著高于种植前期和种植第三年(连作第2年);橡胶林下长裙竹荪种植前期与种植第三年(连作第2年)间的土壤中可培养细菌和放线菌总数差异不显著。橡胶林下长裙竹荪种植第二年(连作第1年)与种植第三年(连作第2年)土壤中可培养真菌总数均差异不显著,但均极显著高于种植前期和种植第1年;种植前期与种植第一年土壤中可培养真菌总数差异不显著。
另外,从表3还可以在连作期内,看出橡胶林下种植长裙竹荪覆土层土壤可培养微生物总数(细菌、真菌和放线菌)在连作期内整体呈现先升高后降低,在种植的前2年中,长裙竹荪覆土层土壤可培养微生物以细菌总数最多,放线菌总数次之,真菌总数最少;在种植的第3年,以真菌总数最多,细菌总数次之,放线菌总数最少。
Table 3. Number of culturable microorganisms in soil
表3. 土壤中可培养微生物数量
注:同一列中,小写字母不同表示差异显著(P < 0.05),大写字母不同表示差异极显著(P < 0.01)。
3.2.2. 土壤功能微生物
Table 4. Number of soil functional microorganisms
表4. 土壤功能微生物数量
注:同一列中,小写字母不同表示差异显著(P < 0.05),大写字母不同表示差异极显著(P < 0.01)。
从表4可以看出,橡胶林下种植长裙竹荪的土壤固氮菌、好氧纤维素分解菌、亚硝化细菌和反硝化细菌数量随长裙竹荪连作年限的增长先升高后降低,并且在种植第二年(连作第1年)最高,其中好氧纤维素分解菌和亚硝化细菌数量在种植第二年(连作第1年)极显著最高;氨化细菌数量随长裙竹荪连作年限的增长先降低后升高再降低的规律,在种植第二年(连作第1年)显著最高。
土壤固氮菌和反硝化细菌数量,除种植第一年和种植第二年(连作第1年)间差异不显著外,其他种植期间差异均极显著;土壤氨化细菌数量和好氧纤维素分解细菌数量,除种植前期与种植第三年(连作第2年)间差异不显著外,其他种植期间差异均极显著;土壤亚硝化细菌数量在各种植期差异均极显著;种植第已年和种植第二年(连作第1年)间的土壤反硝化细菌数量差异不显著,但与种植前期和种植第三年(连作第2年)之间差异及显著。
3.3. 土壤微生物群落功能多样性
3.3.1. 土壤微生物的AWCD值变化特征
从图1可以看出,AWCD值随着培养时间的延长而提高,各种植期样地土壤微生物的AWCD值均在培养的前24 h内变化不大,说明该阶段土壤微生物活性较低;在48 h~168 h内快速升高,说明该阶段土壤微生物对不同碳源的利用加快。
Figure 1. Changes of average light absorption value of soil microorganism with time in different planting years of Dictyophora indusiata.
图1. 长裙竹荪不同种植年限土壤微生物平均吸光值随时间的变化
96 h处的平均吸光值(AWCD)为种植第一年 > 种植第三年(连作第2年) > 种植第二年(连作第1年) >种植前期,种植第一年的土壤微生物的平均吸光值显著高于其他种植期(表5)。
Table 5. Average light absorption value of soil microorganisms in different planting years of Dictyophora indusiata
表5. 长裙竹荪不同种植年限土壤微生物的平均吸光值
注:同一列中,小写字母不同表示差异显著(P < 0.05),大写字母不同表示差异极显著(P < 0.01)。
3.3.2. 土壤微生物对碳源的利用情况
从图2可以看出,土壤微生物对各类碳源的利用均随时间的延长而增加,并基本在120 h达到峰值,土壤微生物对六大类碳源的利用情况随长裙竹荪连作年限发生变化。
对碳水化合物类碳源的利用情况:前48 h内,种植第一年 > 种植第三年(连作第2年) > 种植第二年(连作第1年) > 种植前期;48 h~96 h内,种植第一年 > 种植第三年(连作第2年) > 种植前期 > 种植第二年(连作第1年);96 h~168 h内,种植第一年 > 种植前期 > 种植第三年(连作第2年) > 种植第二年(连作第1年),且种植第三年和种植第二年基本一致;168 h后,种植前期 > 种植第一年 > 种植第二年(连作第1年) > 种植第三年(连作第2年),且种植第三年和种植第二年基本一致。
对氨基酸类碳源的利用情况:前96 h内,种植第一年 > 种植第三年(连作第2年) > 种植第二年(连作第1年) > 种植前期;96 h~192 h内,种植第一年 > 种植前期 > 种植第三年(连作第2年) > 种植第二年(连作第1年);192 h后,种植第三年(连作第2年) > 种植前期 > 种植第一年 > 种植第二年(连作第1年)。
对多聚物类碳源的利用情况:前96 h内,种植第一年 > 种植第三年(连作第2年) > 种植第二年(连作第1年) > 种植前期;96 h~192 h内,种植第一年 > 种植第三年(连作第2年) > 种植前期 > 种植第二年(连作第1年);192 h后,种植第三年(连作第2年) > 种植第一年 > 种植前期 > 种植第二年(连作第1年)。
对酚酸类碳源的利用情况:前48 h内,种植第一年 > 种植第三年(连作第2年) > 种植第二年(连作第1年) > 种植前期;48 h~96 h内,种植第一年 > 种植第二年(连作第1年) > 种植第三年(连作第2年) > 种植前期;96 h~120 h内,种植第一年 > 种植第二年(连作第1年) > 种植前期 > 种植第三年(连作第2年);120 h以后,种植前期 > 种植第一年 > 种植第二年(连作第1年) > 种植第三年(连作第2年)。
对胺类碳源的利用情况:前48 h内,种植第一年 > 种植第三年(连作第2年) > 种植第二年(连作第1年) > 种植前期;48 h~72 h内,种植第一年 > 种植前期 > 种植第三年(连作第2年) > 种植第二年(连作第1年);72 h~120 h内,种植第一年 > 种植前期 > 种植第二年(连作第1年) > 种植第三年(连作第2年);120 h以后,种植前期 > 种植第一年 > 种植第二年(连作第1年) > 种植第三年(连作第2年)。
对羧酸类碳源的利用情况:前120 h内,种植第一年>种植第三年(连作第2年) > 种植第二年(连作第1年) > 种植前期;120 h~168 h内,种植第三年(连作第2年) > 种植第一年 > 种植第二年(连作第1年) > 种植前期;168 h后,种植第三年(连作第2年) > 种植第二年(连作第1年) > 种植第一年>种植前期。
Figure 2. Soil microbiological utilization of carbon sources by Dictyophora indusiata in different planting years
图2. 长裙竹荪不同种植年限土壤微生利用碳源情况
3.3.3. 土壤微生物群落功能多样性指数
橡胶林下种植长裙竹荪的土壤微生物群落功能的McIntosh指数(U值)前96 h内总体趋势上呈现种植第1年最高,种植第三年(连作第2年)其次,种植前期和种植第二年(连作第1年)的变化情况基本一致(图3(a))。96 h内的Shannon指数(H'值)、Pielou指数(E值)、Simpson指数(D值)整体趋势均为种植第一年 > 种植第二年(连作第1年) > 种植第三年(连作第2年) > 种植前期,96 h后各种植期的土壤群落微生物功能多样性指数基本相同(图3(b)、图3(c)和图3(d))。
(a) (b) (c) (d)
Figure 3. Functional diversity index characteristics of soil microbial community in different planting years of Dictyophora indusiata ((a): McIntosh index; (b): Shannon index; (c): Pielou index; (d): Simpson index)
图3. 长裙竹荪不同种植年限土壤微生物群落功能多样性指数特征((a):McIntosh指数;(b):Shannonss指数;(c):Pielou均匀度指数;(d):Simpson指数)
Table 6. Functional diversity index of soil microbial community in different planting years of Dictyophora indusiata
表6. 长裙竹荪不同种植年限土壤微生物群落功能多样性指数
注:同一行中,小写字母不同表示差异显著(P < 0.05),大写字母不同表示差异极显著(P < 0.01)。
从表6可以看出,土壤微生物群落功能McIntosh指数(U值)整体随种植年限的延长呈现先降低后升高的变化规律,种植第一年最高,且显著高于种植前期和种植第二年(连作第1年),但与种植第三年(连作第2年)差异不显著。土壤微生物群落功能的Shananon-Wiener指数(H'值)和Pielou指数(E)随种植年限的延长持续下降,各种植年限间差异均不显著;种植前期土壤微生物群落功能的Simpson指数(D值)和Pielou指数(E)均显著低于其他种植期。
4. 讨论与结论
4.1 讨论
众所周知,细菌是土壤微生物类群中数量最多的,其次为真菌和放线菌,它们在土壤有机碳、氮等物质循环和转化过程中起着十分重要的作用。本研究中,细菌、真菌和放线菌总数整体均随长裙竹荪连作年限的增加呈现先升高后降低的变化趋势,且在种植第三年(连作第2年)显著最低。植物连作影响土壤中有机碳、氮的转化,并导致土壤肥力的下降。土壤微生物作为驱动土壤营养成分转化最主要的因子之一,影响橡胶林下长裙竹荪的生长及品质。另外,放线菌作为抗生素的主要来源,其生长代谢产物中可能存在部分影响长裙竹荪生长的物质。吕杨兰等研究发现连作地土壤竹荪菌的拮抗菌–巨大芽孢杆菌所产蛋白代谢产物,能显著抑制竹荪菌的生长 [6]。另外,江惠华等人研究发现有害微生物潜伏对土壤中原有的有益微生物形成拮抗作用, 使有益微生物活性降低,病虫源逐年增长,导致产地土壤环境恶化,竹荪菌丝降解吸收养分能力降低,从而引起产量降低,产品质量变劣,等级下降,效益受损 [15]。本研究中,橡胶林下不同连作年限的长裙竹荪覆土层土壤三大类微生物(细菌、真菌、放线菌)数量发生了的变化,既影响了橡胶林下长裙竹荪种植的“微生境”,又影响了土壤养分,两者共同响长裙竹荪的产量和品质,该结果与刘亚峰、胡元森等人的研究结果一致 [16] [17]。众多研究发现土壤从细菌主导型转变为真菌主导型是土壤肥力下降的一种表现形式,并最终表现为连作障碍。本研究存在种植长裙竹荪第三年(连作第2年)的橡胶林土壤从“细菌型”转变为“真菌型”的类似现象,而土壤有机质和全钾含量随种植年限的增长而下降,由此可见土壤肥力存在下降的现状。
固氮细菌是一类能固定空气中氮素的细菌,氨化细菌是一类分解含氮有机化合物释放氨的细菌,好养性纤维素分解细菌是一类能分解纤维素的细菌,亚硝化细菌是一类将氨氮转化为亚硝酸盐的细菌,反硝化细菌是一类将亚硝酸盐氮、硝酸盐氮还原成气态氮的细菌 [18] [19]。本研究中,固氮细菌和氨化细菌随着长裙竹荪连作年限的增加呈现先升高后降低的变化规律,土壤中固氮细菌和氨化细菌数量的增减影响土壤中氮元素的循环,间接影响长裙竹荪的生长。好氧性纤维素分解细菌可分解橡胶树的枯枝落叶等凋落物,有利于营养物质的循环和供给,为固氮菌的生长提供氮源,同时提高土壤肥力,有助于林下长裙竹荪的生长 [20]。但是,随着连作时间的增长,土壤透水、通气性下降,好氧性纤维素分解细菌的数量会随之下降,如大豆、马铃薯、太子参和烤烟等 [21]。本研究发现种植长裙竹荪第三年(连作第2年)的橡胶林土壤中的好氧性纤维素分解细菌数量急剧下降的现象。
土壤中含有的各种酶是土壤发挥功能的重要组成部分,土壤酶活性大小表征土壤中各种生化反应的强度,是评估土壤中物质循环、能量代谢强度及土壤质量水平的重要指标。土壤酶活性随种植作物种类、作物生长过程及作物种植方式的变化而变化,同时受土壤内外环境的影响 [22]。土壤酶活性间接地表明土壤腐质化的程度和土壤肥力演变的趋势,马云华等人在黄瓜连作的研究中发现随着种植年限的增加,大多数土壤酶活性呈倒“马鞍形” [23]。本研究中,土壤蔗糖酶活性随长裙竹荪连作年限的增加先降低后升高,呈“马鞍形”;土壤脲酶和脱氢酶活性随长裙竹荪连作年限的增加先升高后降低,呈倒“马鞍形”。
平均吸光值(AWCD值)反映了土壤微生物利用碳源的整体能力和微生物活性,可作为土壤微生物活性的评价指标之一。本研究中,种植第一年的长裙竹荪覆土层土壤微生物的总碳源利用能力最强。McIntosh指数,是基于群落物种的多维空间上的Euclidian距离的多样性指数,Shananon-Wiener指数用来评估物种丰富度,Simpson指数用来评估样品中常见物种的优势度,Pielou指数用来评估物种的均一性 [24]。常颖萃等人对竹荪不同生长期土壤微生物动态变化的研究表明,进入出菇期,代谢物积累加快,影响土壤微生物的数量、种类及其多样性 [25]。本研究中,长裙竹荪的不同种植年限的橡胶林下土壤微生物的物种丰富度、优势度和均匀性差异不显著,种植第二年(连作第1年)的McIntosh指数显著低于种植第一年和种植第三年(连作第2年)。土壤物种丰富度指数随种植年限的增加持续下降,说明橡胶林下长裙竹荪的连作影响了土壤物种的丰富度,其原因可能是长裙竹荪自身的生长代谢及其连作共同影响土壤中部分微生物种类的生长,降低了土壤微生物的物种丰富度,导致微生物物种多样性趋向单一。相关研究表明长裙竹荪种植过程中,其自身的代谢产物具有一定的抑菌性 [26] [27] [28],未及时采摘的长裙竹荪自身发生“自溶”现象,长裙竹荪“自溶”物以及释放的化学自毒物质在环境中的不断积累对土壤中的微生物具有一定的“毒副”作用,导致土壤微生态系统失衡,影响长裙竹荪的正常生长与发育,病原微生物富集,进而影响长裙竹荪的产量和品质。
4.2. 结论
橡胶林下长裙竹荪的连作导致土壤肥力下降、影响土壤微生物群落结构和组成并降低微生物群落功能多样性,长裙竹荪“自溶”产生的毒副物质的积累,影响土壤微生物的种类、功能微生物数量、土壤酶活性以及土壤微生物对碳源的利用能力不同和碳源利用的“偏好性”。橡胶林下种植长裙竹荪连作障碍的产生是土壤肥力状态、土壤微生物、长裙竹荪菌体及其代谢产物共同作用的微生态环境引起的,是一个多因素协同调控下相互作用和影响的复杂过程,并最终影响长裙竹荪的产量和品质。因此,要解决橡胶林下长裙竹荪种植的连作障碍问题,可以通过轮作、间作、套作等不同种植方式,维持或改善土壤肥力及酶活性,土壤微生物的活性及多样性,抑制有害土壤微生物的繁衍,提高土壤微生物群落的结构、组成和功能的稳定性,协调橡胶林–长裙竹荪–土壤生态系统的持续平衡与稳定,是海南发展橡胶林下长裙竹荪种植可持续发展的重要保障。
致谢
本研究受海南省高等学校科学研究项目(Hnky2019-35):海南橡胶林下种植长裙竹荪连作障碍成因初探资助。