益生枯草芽孢杆菌MA139增殖培养基的优化.doc

腿肄虿蚀衿艿薅虿羁肂薁蚈膃莈蒇蚇袃膀莃蚆羅莆蚁蚆肈腿薇蚅膀莄蒃螄袀膇荿螃羂莂芅螂肄膅蚄螁袄蒁薀螁羆芄蒆螀聿葿莂蝿膁节蚀螈袁肅薆袇羃芀蒂袆肅肃莈袅螅芈莄袅羇肁蚃袄聿莇蕿袃膂腿蒅袂袁莅莁袁羄膈蚀羀肆莃薆罿膈膆蒂罿袈莂蒈薅肀芄莄薄膃蒀蚂薃袂芃薈薂羅蒈蒄薂肇芁莀蚁腿肄虿蚀衿艿薅虿羁肂薁蚈膃莈蒇蚇袃膀莃蚆羅莆蚁蚆肈腿薇蚅膀莄蒃螄袀膇荿螃羂莂芅螂肄膅蚄螁袄蒁薀螁羆芄蒆螀聿葿莂蝿膁节蚀螈袁肅薆袇羃芀蒂袆肅肃莈袅螅芈莄袅羇肁蚃袄聿莇蕿袃膂腿蒅袂袁莅莁袁羄膈蚀羀肆莃薆罿膈膆蒂罿袈莂蒈薅肀芄莄薄膃蒀蚂薃袂芃薈薂羅蒈蒄薂肇芁莀蚁腿肄虿蚀衿艿薅虿羁肂薁蚈膃莈蒇蚇袃膀莃蚆羅莆蚁蚆肈腿薇蚅膀莄蒃螄袀膇荿螃羂莂芅螂肄膅蚄螁袄蒁薀螁羆芄蒆螀聿葿莂蝿膁节蚀螈袁肅薆袇羃芀蒂袆肅肃莈袅螅芈莄袅羇肁蚃袄聿莇蕿袃膂腿蒅袂袁莅莁袁羄膈蚀羀肆莃薆罿膈膆蒂罿袈莂蒈薅肀芄莄薄膃蒀蚂薃袂芃薈薂羅蒈蒄薂肇芁莀蚁腿肄虿蚀衿艿薅虿羁肂薁蚈膃莈蒇蚇袃膀莃蚆羅莆蚁蚆肈腿薇蚅膀莄蒃螄袀膇荿螃羂莂芅螂肄膅蚄螁袄蒁薀螁羆芄蒆螀聿葿莂蝿膁节蚀螈袁肅薆袇羃芀蒂袆肅肃莈袅螅芈莄袅羇肁蚃袄聿莇蕿袃膂腿蒅袂袁莅莁袁羄膈蚀羀肆莃薆罿膈膆蒂罿袈莂蒈薅肀芄莄薄膃蒀蚂 益生枯草芽孢杆菌MA139增殖培养基的优化郭小华1 陆文清1,2 邓萍1 黄德仕(1.中国农业大学农业部饲料工业中心，北京100193；2.北京晟育达生物科技有限公司，北京100193)摘 要 为得到益生芽抱杆菌Bacillus subtilis MA139产芽抱的最佳的培养基，以摇瓶发酵的方法，用P1ackett-Burman设计从l1种原料中筛选出4种对芽抱产量有显著影响的因素，即玉米粉、大豆粉、蛋白胨和MnSO4H2O，然后针对这4个主要因素，用最陡爬坡试验及中心组合设计优化产芽抱的最佳培养基。结果证明，当培养基的配方为：玉米粉3.17g/L、大豆粉5.80g/L、蛋白胨3.62g/L、MnSO4H2O1.06g/L、葡萄糖5g/L、尿素3g/L、MgSO47H2O1.5g/L和KH2PO43g/L时，MAI39发酵36h细菌总数可以从8.32108cfu/ml提高到3.10109cfu/mL，芽抱率达到96%。试验表明，通过统计优化培养条件可有效提高B.subtilis MA139产芽孢的得率。关键词:枯草芽孢杆菌；发酵；培养基优化；试验设计抗生素添加剂在动物日粮中广泛使用，但其弊端和危害日益体现出来1，而抗生素替代品的微生态制剂在畜牧业已得到应用2。微生态制剂常用的细菌是乳酸菌，而芽孢杆菌作为动物肠道的外源菌群，也有广泛应用，它能够改善动物胃肠道的微生态平衡，促进动物生长，提高饲料利用率3-6。优良的微生态制剂必须有足够的活菌数量7。因此，为了尽量提高产品芽孢的数量，需要对其产芽孢的培养基进行系统优化，以降低生产成本，提高芽孢得率8。芽孢杆菌液体深层发酵培养基的优化一般采用单因子试验以及结合均匀设计和正交设计的方法9-10。目前趋向采用统计软件构建高效试验设计来进行微生物培养基的优化，通过对试验结果进行数学模拟和优化，可以简化试验步骤，提高准确性11-14。本试验旨在以本试验室分离筛选得到的一株枯草芽孢杆菌MA139为出发菌株，优化产芽孢的最适培养基。1、材料与方法1.1、菌株枯草芽孢杆菌MA139，本实验室自行分离和保存。该菌株经中科院微生物所鉴定，其16S rRNA基因序列在GenBank的登录号为DQ415893。该菌株具有饲用益生菌的应用潜力。1.2、培养基1)种子及斜面保存培养基(SC培养基)成分为葡萄糖5g/L；蛋白胨5g/L；酵母粉1g/L；牛肉膏3g/L；MgSO47H2O 0.5g/L；MnSO4H2O 0.005g/L；pH7.0。固体培养基则在液体培养基中添加1.8%的琼脂作为凝固剂。2)发酵培养基分别按照Plackett-Burman设计、最陡爬坡试验设计和中心组合设计的方法来配制。1.3、试剂和原料玉米、大豆为市售，粉碎后过2O目筛备用；试验中其他试剂均为分析纯。1.4、培养条件1)种子培养条件 250 mL三角瓶装种子培养基5OmL，接斜面活化后的种子一环，37，225r/min旋转振荡培养20h。2)发酵培养条件250mL三角瓶装发酵培养基5OmL，按照4的接种量接种后在37，225r/min旋转振荡培养36h。1.5、细菌总数和芽孢总数的测定平板稀释法计数活菌数。芽孢计数则将稀释液放入8O水浴中处理20min后培养计数。1.6、培养基优化试验设计1)Plackett-Burman设计在优化初期利用Plackett-Burman设计法15，以葡萄糖(X1)、可溶性淀粉(X2)、玉米粉(X3)、大豆粉(X4)、(NH4)2SO4(X5)、尿素(X6)、蛋白胨(X7)、MgSO47H2O(X8)、KH2PO4(X9)、CaC12(X10)、MnSO4H2O(X11)等为试验因子，对这11种原料进行全面考察，并将其编码为-1和+1的高低二阶层，其所代表的质量浓度值，见表1。选用”n=12的Plackett-Burman设计，并加上3个中心点来筛选和探讨11个试验因子对细菌总数和芽孢总数的影响，每个试验结果报告为2个重复试验的平均值，见表2。表1 Plackett-Burman试验设计中各参数的编码值及其质量浓度 g/L独立变量编码水平-1+1葡萄糖(X1)515可溶性淀粉(X2)515玉米粉(X3)515大豆粉(X4)1030(NH4)2S04(X5)26尿素(X6)39蛋白胨(X7)515MgS047H2O(X8)0.51.5KH2P04(X9)13CaCl2(X10)26MnS04H2O(X11)0.050.52)最陡爬坡试验设计响应面拟合方程只有在考察的紧接邻域里才充分近似真实结果，故只有在最先逼近最大目标产物产量区域后才能建立有效的响应面拟合方程16。以试验值变化的梯度方向为爬坡方向，根据各因素效应值的大小确定变化步长，从而快速逼近产量最大值。根据Plackett-Burman设计的结果，来确定是否进行最陡爬坡试验。3)中心组合设计根据Plackett-Burman设计筛选出的试验因子和最陡爬坡试验确定的逼近浓度，以JMP5.0软件设计中的DOE程序进行中心组合试验，拟合数据得到一个描述因变量(细菌总数)与自变量(培养基组分)关系的二阶经验模型11-17：y=bo+biXi +bijXiXj+biiX2ii；式中：y为预测响应值，即细菌总数(1gNtot(cfu/mL)；b为回归系数；Xi为自变量的编码水平。用JMP5.0软件对试验数据进行回归拟合，并对拟合方程作显著性检验和方差分析。2、结果与讨论2.1、影响芽孢杆菌增殖的重要因素按照Plackett-Burman试验设计，在发酵36h时取样测定发酵液当中的细菌总数(total number，Ntot)和芽孢总数(spore number，Ns)。Plackett-Burman试验设计结果见表2并对其进行一阶模式回归分析(表3)。细菌总数和芽孢总数2种结果回归分析的检定系数分别为0.96和0.94，说明回归分析能够确切地描述试验数据。根据分析结果，发现对B.subtilis MA139发酵起主要作用有4种原料，即玉米粉、大豆粉、蛋白胨和MnSO4H2O，并对11种原料成分归纳调整，作为下一步试验分析的依据。质量浓度调整说明见表4。表2 Plackeet-Burman试验设计结果(n=12)编号变量水平响应值X1X2X3X4X5X6X7X8X9X10X11lgNtot/(cfu/ml)lgNs/(cfu/ml)1-11-111-11-1-11-18.698.482-1-1-1-11-111-1-119.379.353000000000009.018.9941-1-111-11-11-1-18.418.405-111-11-1-111-1-19.319.296111111111118.438.327-11-1-1-111-11-119.259.24811-11-1-1-11-1-119.079.069-1-111-1111-1-1-18.238.1810000000000009.059.0511111-1-1-11-1-11-18.308.26121-1-1-1-11-1111-19.159.1513-1-111-1-1-1-11119.068.9414000000000009.029.01151-11-111-1-1-1-118.798.74 表3 Plackeet-Burman试验设计回归分析结果 cfu/L截距参数细菌总数Ntot芽孢总数Ns参数估计tP参数估计tPintercept8.876203.0440.0008.831163.2940.000X1-0.147-3.0010.058-0.129-2.1360.122X20.0030.0680.950-0.009-0.1520.889X3-0.152-3.1030.053*-0.163-2.6880.075*X4-0.190-3.8880.030*-0.221-3.6520.035*X5-0.005-0.1020.925-0.021-0.3450.753X6-0.082-1.6710.193-0.099-1.6400.200X7-0.173-3.5470.038*-0.159-2.6330.078*X80.0881.8070.1680.1081.7880.173X90.0971.9780.1420.1061.7500.178X10-0.005-0.1020.925-0.034-0.5650.612X110.1573.2050.049*0.1582.6050.080*P=0.082，RSq=0.96 P=0.125，RSq=0.94 注：*为P0.1，*0.05。表4 Plackeet-Burman试验设计各因子结果的质量浓度分析试验因子质量浓度（g/L）说明调整前调整后葡萄糖(X1)5-155对细菌总数和芽孢总数为负相关，对细菌总数影响达到显著水平，而对芽孢总数的影响不显著，为了减少试验因子，维持在低水平可溶性淀粉(X2)5-15删除无显著影响，可信度大于20%，删除此因子玉米粉(X3)5-151-5负显著，降低筛选浓度大豆粉(X4)10-301-10负极显著，降低筛选浓度(NH4)2S04(X5)2-6删除无显著影响，可信度大于20%，删除此因子尿素(X6)3-93不显著，负相关，可信度小于20%，为减少试验因子，维持在低水平蛋白胨(X7)5-151-5负显著，降低筛选浓度MgS047H2O(X8)0.5-1.51.5不显著，可信度小于20%，为减少试验因子，维持在高水平KH2P04(X9)1-33不显著，正相关，P值相对较小，维持在高水平CaCl2(X10)2-6删除不显著，可信度大于20%，删除此因子MnS04H2O(X11)0.05-0.50.5-1.5显著，正相关，提高筛选浓度2.2、最陡爬坡试验选取对 B.subtilis MA139发酵 4种主要原料。对Plackett-Burman试验设计中得到的结果进行t检验，比较中心点水平的平均响应值与试验点处的平均响应值间的差异显著性。试验结果显示，中心点水平与试验点间在a=005水平上不显著(P值分别是0.460和0.388)，这就表明本试验设计中选取的中心点水平还没有接近最大的响应区域内，因此还必须进一步寻找中心点值，也就是最陡爬坡试验。以回归系数最大的大豆粉(X4)为基准，以每次减少浓度1%为基本步长()。处理5中的细菌总数和芽孢总数达到了最高，这一结果表明，处理5的玉米粉(X3)、大豆粉(X4)、蛋白胨(X7)、MnSO4H20(X11)这4种物质浓度已在最优点附近，可以此浓度为中心点采用中心组合设计来对培养基作进一步的优化(表5)。而且在Plackett-Burman设计及最陡爬坡试验中，细菌总数与芽孢总数的分析结果一致，故在以下的中心组合设计仅以发酵液的细菌总数作为考察目标。表5 最陡爬坡试验设计及其结果处理序号培养基成分/（g/L）结果玉米粉（X3）大豆粉（X4）蛋白胨（X7）MnSO4H20（X11）lgNtot/(cfu/ml)lgNs/(cfu/ml)15.0010.005.000.509.339.2824.529.004.520.619.379.3034.088.004.080.729.139.1243.647.003.640.839.349.2953.206.003.200.949.439.4062.765.002.761.059.399.3472.324.002.321.169.109.0181.883.001.881.279.179.1591.442.001.441.389.068.95101.001.001.001.508.918.842.3、中心组合设计以最陡爬坡试验中处理5的玉米粉(X3)、大豆粉(X4)、蛋白胨(X7)、MnSO4H2O(X11)的浓度为中心点进行优化，选择中心试验点为3，星号臂长=1.547。各自变量水平见表6，试验设计及结果见表7。分析结果表明(表8)，模型极显著(P0.01)，并且有很好的确定系数，RSq=0.93。这表明B.subtilis MA139细菌总数93%的可靠性可由模型来预测。并可以通过回归分析结果得出试验模型：y=9.434-0.003X3-0.025X4+0.028X7-0.013X11-0.006X3X4-0.02X3X7-0.009X4X7-0.012X3X11-0.026X4X11-0.03X7Xll-0.068X3X3-0.013X4X4-0.043X7X7-0.025X11X11 自变量前的数值为回归模型对应项的 P值。根据JMP模型的分析结果以及二次回归方程可以得出，此模型具有最大响应值，极值分析结果见表9。根据结果换算成4个主因素的实际浓度值，即玉米粉3.17g/L，大豆粉5.80g/L，蛋白胨3.62g/L，MnSO4H2O1.06g/L。采用优化出的培养基，芽孢杆菌最终细菌浓度预测值为2.84109cfu/mL。为了检验模型预测的准确性，采用优化的培养基重复试验3次，发酵液当中平均细菌数为 2.91l09cfu/mL，这表明实验值和实际值之间具有良好的拟合性，优化模型可靠。表6 中心组合设计各因子及其编码值 g/L独立变量编码水平-1.547-1011.547玉米粉（X3）1.652.203.204.204.75大豆粉（X4）4.455.006.007.007.55蛋白胨（X7）1.652.203.204.204.75MnSO4H20（X11）0.550.690.941.191.32表7 中心组合设计及其结果编号变量水平lgNs/(cfu/ml)编号变量水平lgNs/(cfu/ml)X3X4X7X11X3X4X7X111-1-1119.3615-111-19.392-1-11-19.3616-11-119.2731.5470009.28171-1-1-19.244-1.5470009.271811-1-19.3150001.5479.341900-1.54709.306-11119.23200-1.547009.467111-19.3021-1-1-119.298-1-1-1-19.2022001.54709.37900009.422300009.441000009.4324-11-1-19.20111-1119.292501.547009.341211119.2026000-1.5479.41131-1-119.34271-11-19.361411-119.20表8 中心组合设计(CCD)回归分析结果回归系数参数估计标准差tP回归系数参数估计标准差tPIntercept9.4340.015638.8100.000Intercept9.4340.015638.8100.000X3-0.0030.007-0.3910.702X3X11-0.0120.007-1.6350.128X4-0.0250.007-3.8600.002X4X11-0.0260.007-3.4910.004X70.0280.0074.2820.001X7X11-0.0300.007-4.0670.002X11-0.0130.007-1.9610.074X3X3-0.0680.009-7.6720.000X3X4-0.0060.007-0.8550.409X4X4-0.0130.009-1.5220.154X3X7-0.0200.007-2.6150.023X7X7-0.0430.009-4.8910.000X4X7-0.0090.007-1.2570.233X11X11-0.0250.009-2.8050.016表9 响应面极值分析临界值估计值细菌总数Nt/(cfu/ml)稳定点模型X3X4X7X11-0.03-1.20.420.129.4542.84