纤维素分解菌与乳酸菌复合发酵条件的响应曲面优化分析

1.1 主要材料

本研究主要菌种为购自中国农业微生物菌种保藏管理中心(agricultural culture collection of china,ACCC)的植物乳杆菌(Lactobacillus plantarum,ACCC 11016)、戊糖片球菌(Pediococcus acidilactici,ACCC 05481)、枯草芽孢杆菌(bacillus subtilis,ACCC 19374)、黑曲霉(aspergillus niger,ACCC 30134)和绿色木霉(trichoderma viride,ACCC 30595)。购自中国工业微生物菌种保藏管理中心(China Center of Industrial Culture Collection,CICC)的布氏乳杆菌(Lactobacillus buchneri,CICC 20293)。

1.2 试验设计

通过查阅关于添加同、异型发酵乳酸菌(LAB)发酵青贮玉米的相关文献,初步确定两类LAB 的添加量为1×105、3×105 和5×105 cfu·g-1 (即5、5.48、5.70 lg cfu·g-1),其中同型发酵乳酸菌为植物乳杆菌和戊糖片球菌等比例混合(菌种数量比为1∶1)的复合菌,异型发酵LAB为布氏乳杆菌。纤维素分解菌菌群的添加比例为,黑曲霉∶绿色木霉∶枯草芽孢杆菌分别为1∶1∶2、1∶2∶1和2∶1∶1,添加量分别为0.1、0.2和0.3％。然后进行响应曲面分析,优化各菌种比例和添加量。

1.3 指标测定及方法

将各个菌种复壮后,接种至选择性培养基中扩繁,其中LAB采用MRS液体培养基[配方:酪蛋白胨10.0g,牛肉粉8.0g,酵母粉4.0g,葡萄糖20.0g,硫酸镁0.2g,乙酸钠5.0g,柠檬酸二铵2.0g,磷酸氢二钾2.0g,硫酸锰0.05g,吐温801.0g,蒸馏水1L,pH (6.2±0.2)],真菌Aspergillus niger、Trichoderma viride采用PDA 液体培养基(配方:马铃薯提取液1L,葡萄糖20.0g,自然pH),bacillus subtilis采用LB培养基(配方:酵母提取物5.0g,蛋白胨10.0g,氯化钠10.0g,蒸馏水1L,pH7.0±0.2)。然后进行响应曲面优化,设计方法为BBD(Box-Behnken Design),本研究中响应值为LAB和真菌的数量,自变量分别为同型发酵LAB添加量、异型发酵LAB添加量、纤维素分解菌菌群比例和剂量,分别以A、B、C和D代表,每个自变量的3个试验水平(低、中、高)编码分别为-1、0、1(表1)。

然后将各个扩繁后的菌液按照BBD 设计的添加比例和量进行混合,厌氧培养3d后(35℃,pH4.5)测定液体培养基中LAB和真菌的数量,具体采用传统菌种计数法计数，每个梯度做3个重复。取菌落数在30~300的平板作为有效计数板,菌落以肉眼可见为准,每克鲜样(FM)中微生物的数量,即菌落形成单位(cfu),计算方法为:

微生物数量＝菌落数×稀释倍数×1000/涂布吸样量。

1.4 数据统计分析

数据的初步整理和统计采用Microsoft Excel 2011软件。采用DesignExpert软件进行响应曲面优化,设通过最小二乘法拟合的二次多项式模型为:

式中:Y 为预测响应值,即LAB或真菌数量的预测值;xi 和xj 为自变量的编码值,c0 为常数项,ai 为线性系数,bij 为二次项系数,n 为因子数(n ＝4)。按照BBD试验设计的统计学要求对二次多项式中各试验进行回归拟合,拟合后获得回归方程,运用软件中自带功能,选取两个预测值均为最大(Y1 为LAB,Y2 为真菌数量),以此得到最佳的各水平编码值,最终计算求得各个变量的最佳组合。

2.1 模型建立与显著性检验分析

采用BBD设计得到变量代码,按照变量代码进行试验,结果测得LAB最大值为9.96 lg cfu·g-1,其代码为1(A)、1(B)、0(C)、0(D);最小值为7.93 lg cfu·g-1,其代码为-1(A)、0(B)、-1(C)、0(D),总体平均值为8.63 lg cfu·g-1,总体标准差(SD)为0.16,变异系数(CV)为1.89％.真菌测得最大值为5.23lgcfu·g-1,其代码为0(A)、1(B)、0(C)、1(D);最小值为3.92lgcfu·g-1,其代码为1(A)、1(B)、0(C)、0(D),总体平均值为4.49lgcfu·g-1,总体SD为0.27,CV 为6.05％(表2)。

利用DesignExpert软件,将表2的Y1 和Y2 数据按二次多项式模型进行多元回归拟合,获得Y1 和Y2 对编码自变量A、B、C和D的二次多项回归方程(表3):

 Y1 ＝8.56＋0.50A ＋0.50B ＋0.07C ＋0.04D ＋0.09A2＋0.11B2-0.01C2－2.58×10－3D2＋0.21AB ＋0.06AC-0.04AD-0.04BC-0.06BD＋0.08CD;

Y2＝4.38-0.02A -0.05B－8.33×10－3C＋0.32D－3.58×10－3A2 -0.04B2 -0.03C2＋0.35D2 -0.17AB＋0.30AC-0.01AD-0.18BC＋0.18BD-0.04CD。

通过对Y1、Y2 两个模型的方差分析可知(表3),本研究所选用的二次多项式模型显著性为:Y1,F＝17.23,P＜0.01;Y2,F＝2.83,P＝0.03.失拟项均不显著(Y1:P＝0.50,Y2:P＝0.68).决定系数R2 分别为Y1＝0.95,Y2 ＝0.74.矫正决定系数R2adj 分别为Y1＝0.89,Y2＝0.48.显著性检验结果表明(表3),Y1模型中试验因素A 和B对LAB数量的曲面效应影响均极显著(P＜0.01),因素A 和B的交互作用影响显著(P＝0.02),其余因素及各因素间的交互作用影响均不显著(P＞0.05).Y2 模型中仅试验因素D、因素A 和B的交互作用对真菌数量的曲面效应影响显著(P＜0.05),其余各因素及因素间的交互作用影响均不显著(P＞0.05)。

2.2 复合发酵真菌和乳酸菌数量的响应

曲面分析通过对Y1 所绘制的响应曲面图及等高线图可知(图1),A 和B两个因素的交互对Y1 影响最大,当A和B的代码越接近于1,Y1 的值越趋于最大值(9.96 lg cfu·g-1)。A 和C、D 因素的交互作用中,起主要作用的均为因素A.B和C、D因素的交互作用中,起主要作用的均为因素B.C和D两个因素的交互对Y1 的值基本无影响(P ＝0.36)。总体来看,仅A 和B因素的交互作用影响显著(P ＝0.02),其余两因素的交互作用对Y1 的值影响均不大(P＞0.05)(表3)。通过对Y2 所绘制的响应曲面图及等高线图可知(图2),对Y2 的值起主要作用的为因素D,且因素D 和其他因素的交互作用中因素D也也起到主要作用,但总体来看,除A、B因素交互作用外,其余各因素之间的交互作用对Y2 的值影响均不大(P ＞0.05)(表3)。最终利用软件求出Y1 和Y2 均为最大值的代码解(最优解)为A＝1.000,B＝0.853,C＝1.000,D＝0.996,将代码解换算成实际自变量可得,同型发酵LAB添加量为5×105 cfu·g-1(5.70 lg cfu·g-1),异型发酵LAB 添加量为4.7×105 cfu·g-1 (5.67 lg cfu·g-1),纤维素分解菌比例为2∶1∶1,添加量为0.3％。

以往研究很少对多种类菌种的复合发酵进行研究,类似研究多在发酵堆肥中使用。本研究对LAB和纤维素分解菌复合发酵的发酵比例剂量进行了优化,而在实际发酵过程中,由于真菌和细菌的生长环境不同(好氧/厌氧)而在复合培养情况下仅能确定一种生长环境,因此本研究模拟青贮发酵初期的好养环境对复合菌种进行了优化,得到的结果能很好地为今后青贮复合发酵添加剂的开发提供理论方法。此外,由于多种菌剂的复合发酵体系比较复杂,本研究仅对其进行了初步的探索,因此还需要继续开展实践研究。

本研究对复合菌种发酵体系内的LAB数量进行优化,而未对传统概念上认识的pH、乳酸或乙酸含量等发酵产物进行优化。其一方面原因为,LAB的代谢产物较多,产生情况较为复杂,且受发酵环境和底物的影响极大,而无论其代谢产物种类如何,其必然受自身数量的影响。另一方面,由于本研究中探究复合发酵体系,单纯以某种LAB的代谢产物进行优化,并和真菌数量进行响应面分析,不能反映整体情况,且指标上不具有统一性。另外,选取多个LAB的多种代谢产物进行优化,其操作性较差,不易进行响应曲面优化。复合体系中将真菌数量作为另一个优化指标,其主要原因为本研究选取的能产生胞外酶的微生物,其产酶含量与数量呈显著正相关关系,以数量为响应值即能很好的反映产酶的含量。

众所周知,LAB有抑制真菌生长的特性,而本研究需要真菌在青贮玉米发酵初期发挥其产胞外酶的特性,因此对好氧情况下真菌和LAB均达到最大值的发酵条件进行了优化。研究发现,在确定的4个发酵条件中,同、异型发酵LAB的数量以及其交互作用对总体LAB数量影响最大(表3),且当两类LAB数量最多时,总体LAB数量也最多,这与LAB生物学基础理论相同。在发酵体系内加入的纤维素分解菌菌群的比例和剂量两个因素对LAB数量无显著影响(P ＞0.05),也进而证明了LAB数量并未受纤维素分解菌的抑制,这与诸多的研究结果相似,他们也都认为在厌氧条件下真菌不会抑制LAB 的生长。另外,本研究为多种类菌种的复合发酵,因此选用的培养条件不可能达到所有菌种最适的生长环境,且实际生产中青贮玉米的发酵环境也难于控制,所以本研究仅模拟了青贮发酵初期(3d内)的环境进行研究，这可能也是导致LAB生长不受真菌抑制的主要原因。本研究还发现,在确定的4个发酵条件中,影响真菌数量的主要因素为纤维素分解菌的添加剂量,且同、异型发酵LAB添加量的交互作用对真菌也有显著影响(P＜0.05),而纤维素分解菌的添加比例对其无显著影响(表3)(P＞0.05),由此也进一步证实了LAB对真菌的抑制作用。另外,纤维素分解菌的比例和添加量的交互作用中以添加量为主导,其添加量越高,真菌数量越高(图2。这也说明了真菌添加剂量对其数量起到了决定性作用。值得注意的是,当纤维素分解菌的添加比例分别与同/异型发酵LAB添加量交互时,其产生了较大变异(图2),这可能是不同种LAB抑制真菌的种类不同所致，其具体原因还有待进一步分析研究。除此之外,当纤维素分解菌的添加剂量分别与同/异质型LAB 添加量交互时,抑制真菌的效果异性发酵优于同型发酵(图２),这表明异型发酵LAB产抑制真菌生长的乙酸含量高于同型发酵LAB。

Y2 表示真菌数量。

Y2 is the fungi number。