Comparison of different bacteria growth models on chilled Pseudosciaena crocea
-
摘要:
主要研究了养殖大黄鱼0、5℃冷藏过程鱼肉中细菌变化情况,采用修正的Logistic和Gompertz模型拟合细菌生长曲线, 经非线性回归分析求出生长动力学参数和细菌生长预测模型。结果表明,所建立的模型相关系数均大于0.99, 能有效描述菌落总数的动态变化,预测不同贮藏时间内的菌落总数。为验证模型的适用性,把不同时间的菌落总数预测值和实测值比较,依据均方根评价建立的细菌生长预测模型的适用性。0、5℃冷藏过程Gompertz模型均方根分别为0.077和0.100,Logistic模型均方根分别为0.114和0.138,2种模型相比,Gompertz预测结果更为理想。
-
关键词:
- 养殖大黄鱼 /
- 冷藏 /
- Gompertz模型 /
- Logistic模型
Abstract:Bacterial changes of fish flesh in cultured Pseudosciaena crocea stored aerobically at 0, 5℃ were mainly analyzed, bacteria growth curves of chilled cultured P.crocea was fitted using modified Logistic and Gompertz equations.Kinetic parameters of bacterial growth and predictive model were developed with nonlinear estimation equations. The results showed that the correlation coefficient of developed models was over 0.99, indicating predictive models make it possible to describe the dynamic changes and predict the bacterial number stored at different time. Difference between predicted values and observed values were compared with root mean squares (RMS) for validating the goodness of predictive models of the bacterial growth. RMS of the Gompertz model were 0.077 and 0.100, and RMS of the Logistic model were 0.114 and 0.138 at 0, 5℃, respectively. Predicted values of the Gompertz were significant compared with the Logistic model.
-
Keywords:
- cultured Pseudosciaena crocea /
- cold storage /
- Gompertz model /
- Logistic model
-
20世纪80年代开始,澳大利亚、美国、英国、丹麦等国对食品微生物生长、失活、残存和风险评估模型进行了广泛研究[1-3],目的是使用微生物预测模型描述不同物理和化学条件下微生物变化情况,客观评价食品在加工、流通、销售、贮藏过程中各种影响因子对食品安全和品质的影响,运用数学模型对微生物的动态变化进行快速有效预测和评价。近年来,美国、英国、丹麦等国致力于构建微生物特征数据库,开发了Food Micromodel, Pathogen Modeling Program, Seafood Safety and Spoilage Predictor等专家系统,用于食品品质预测和安全评估,实现关键控制点分析,协作管理者进行管理。
食品微生物生长模型通常有Monod, Gompertz, Baranyi, Logistic等模型,其中修正的Logistic和Gompertz模型被广泛采用。Baranyi和Robert[4]综合Logistic模型和Michaeli-Menton模型开发了Baranyi模型,该模型相对来说比较复杂。本研究中养殖大黄鱼(Pseudosciaena crocea)取自福建闽东三都湾养殖区(26°35′~26°55′N),2005年3月海水温度变化幅度为10.5~11.1℃。依据以前对0、5℃恒温冷藏过程中采用感官、挥发性盐基氮(TVBN)、菌落总数(TVC)进行品质变化研究的基础上,本文采用修正的Logistic和Gompertz模型拟合细菌生长曲线,经非线性回归分析,求出预测模型的动力学参数,建立了细菌生长预测模型,并对2种模型的适用性进行了分析。
1. 材料与方法
1.1 试验材料
大黄鱼在福建省闽东三都湾养殖场捕获(2005年3月),立即放入冰水中冷休克。选用大小基本一致的个体(300~400 g·尾-1)。
1.2 贮藏实验
到达实验室后,将鱼腹部朝上装入下有篦子能沥水的塑料盆中,盖上有漏气孔的盖,分别放入高精度低温培养箱(Sanyo MIR 153, 日本)中,控制贮藏温度在0、5±0.1℃。每隔适当时间取出试样鱼进行感官鲜度评价、TVBN、TVC测定,综合确定产品货架期,本次研究未给出和分析感官、TVBN实验数据。
1.3 样品处理
随机抽取2尾试样鱼,先进行生鱼感官评价,然后去鳞去内脏去腮洗净,用干净吸纸擦干。沿脊骨剖切,取半条鱼肉(带鱼皮),用组织捣碎机打碎,用于TVBN和TVC测定;其余半条鱼蒸熟后用于感官评价。
1.4 感官鲜度评价
由6名经过训练的评价员组成感官评价小组,评价生鱼的气味和蒸熟后鱼的气味和味道。采用3分法进行评分,0为最好品质,1为鲜鱼的鲜香味消失,0~1为高品质期,2为明显出现臭味和异味即可接受界限[5]。当半数或以上评价员评价2或以上时,即为货架期终点(感官拒绝点)。
蒸熟时将带头的半条鱼分别用铝箔包好,待锅中水沸腾后,放入锅内的金属篦子上,盖上锅盖蒸20 min,打开锅盖后立即进行感官评价。
1.5 菌落总数计数
称取鱼肉浆10.0 g,加入90 mL 0.1%蛋白胨无菌生理盐水,高速振荡后,以10倍稀释将鱼肉浆稀释,取3个浓度合适的稀释液0.1 mL,涂布于标准琼胶培养基(中国科学院上海昆虫科技开发公司康乐培养基有限公司)平板表面。每个稀释液涂布2个平皿,25℃培养48 h。
1.6 细菌生长曲线拟合
0、5℃大黄鱼贮藏实验得到的细菌增殖动态数据,采用修正的Gompertz方程和Logistic方程[6]描述其生长动态。修正Gompertz方程如下:
$$ \lg N(t)=A+C \times \exp \{-\exp [-B \times(t-M)]\} $$ (1) 式中t为时间(h),N(t)为t时的菌数(lgCFU·g-1),A为初始菌数N0(lgCFU·g-1),C为最大菌数Nmax与N0之差(lgCFU·g-1),M为1/2Nmax时的时间(h),B为M时比生长速率(h-1),最大比生长速率μmax=BC/e。
修正Logistic方程如下:
$$ \lg N(t)=A+C / \exp \{1+[-B \times(t-M)]\} $$ (2) 式中t为时间(h),N(t)为t时的菌数(lgCFU·g-1),A为最小菌数N0(lgCFU·g-1),C为最大菌数Nmax与N0之差(lgCFU·g-1),B为最大比生长速率μmax (h-1),M为1/2Nmax时的时间(h)。
实验数据用Statistica (Release 5.5)统计软件采用最小平方差法进行非线性回归。
1.7 细菌生长模型的可靠性评价
细菌生长动力学模型求得的预测值,与大黄鱼贮藏实验所得的细菌生长的实测值比较,依据均方根[7]评价建立的生长动力学预测模型的可靠性。均方根用下式表示:
$$ \text { RMS }=\left[\frac{\sum_i\left(\lg N_{i, \text { predicted }}-\lg N_{i, \text { obsered }}\right)^2}{n}\right]^{0.5} $$ (3) 式中lgNi, predicted为冷藏实验中菌落总数预测值(lgCFU·g-1),lgNi, observed为菌落总数实测值(lgCFU·g-1),n为测试次数。
2. 结果
2.1 0℃冷藏大黄鱼细菌生长模型
0℃冷藏期间菌落总数实测值和预测值见表 1。开始贮藏4~5 d,细菌生长缓慢,菌落总数小于6.0 lgCFU·g-1,这是由于暖带海域水温高于温带海域,中温菌数量多, 冷藏过程不耐低温,生长受到抑制,甚至死亡。同时一些嗜冷菌逐渐适应低温环境,随着贮藏期的延长,细菌生长加快,进入指数生长期,好冷菌增殖速度逐渐达到高峰,细菌数呈几何级数增加,腐败终点依据感官评分、TVBN、TVC(7.31 lgCFU·g-1)确定产品货架期为409 h。
时间/h
time实测值
observed valuesGompertz预测模型predictive model Logistic预测模型predictive model 预测值
predicted values残存值
residual values预测值
predicted values残存值
residual values0 5.20 5.20 0.00 5.21 -0.01 49 5.10 5.20 -0.10 5.24 -0.14 76 5.25 5.21 0.04 5.29 -0.04 123.5 5.34 5.42 -0.08 5.48 -0.14 167 6.05 5.95 0.10 5.89 0.16 214 6.50 6.55 -0.05 6.52 -0.02 262 6.89 6.96 -0.07 6.99 -0.10 286 7.10 7.08 0.02 7.12 -0.02 334 7.34 7.22 0.12 7.23 0.11 383 7.22 7.29 -0.07 7.27 -0.05 409 7.31 7.31 0.00 7.28 0.03 图 1、方程4和方程5是采用修正Gompertz方程和Logistic方程回归,得到的0℃冷藏大黄鱼中菌落总数变化曲线和方程,生长动力学参数见表 2。
参数
parametersN0
(lgCFU·g-1)Nmax
(lgCFU·g-1)μmax(h-1) M(h) Ns
(lgCFU·g-1)货架期/h
shelf lifeGompertz模型model 5.20 7.33 0.013 170.88 7.31 409 Logistic模型model 5.20 7.26 0.024 195.56 注:Ns为感官终点细菌数
Note: Ns denotes population at the time of organoleptic rejection.$$ \begin{aligned} & \quad\quad \lg N(t)=5.20+2.13 \times \exp \{-\exp [-0.015 \times \\ & (t-170.85)]\} \end{aligned} $$ (4) $$ \begin{aligned} & \quad\quad \lg N(t)=5.20+2.06 /\{1+\exp [-0.024 \times(t- \\ & 195.56)]\} \end{aligned} $$ (5) 2.2 5℃冷藏大黄鱼细菌生长模型
5℃冷藏期间菌落总数实测值和预测值见表 3。贮藏初期,细菌生长缓慢,进入指数生长期,细菌增殖速度快于0℃冷藏大黄鱼中细菌增殖速度(μmax见表 2,表 4),货架期终点菌落总数为7.34 lgCFU·g-1,产品货架期为291 h。
时间/h
time实测值
observed valuesGompertz预测模型predictive model Logistic预测模型predictive model 预测值
predicted values残存值
residual values预测值
predicted values残存值
residual values0 5.20 5.20 -0.00 5.22 -0.02 49 5.41 5.21 0.20 5.29 0.12 76 5.24 5.29 -0.05 5.39 -0.15 123.5 5.79 5.82 -0.04 5.81 -0.02 167 6.60 6.50 0.10 6.48 0.18 216 6.90 7.03 -0.14 7.08 -0.18 264 7.40 7.30 0.10 7.30 0.10 291 7.34 7.37 -0.03 7.34 -0.00 参数
parametersN0
(lgCFU·g-1)Nmax
(lgCFU·g-1)μmax(h-1) M(h) Ns
(lgCFU·g-1)货架期/h shelf life Gompertz模型model 5.20 7.29 0.016 137.17 7.34 291 Logistic模型model 5.20 7.38 0.030 154.99 图 2、方程6和方程7是采用修正Gompertz方程和Logistic方程回归,得到的5℃冷藏大黄鱼中菌落总数变化曲线和方程,生长动力学参数见表 4。
$$ \begin{aligned} & \quad\quad \lg N(t)=5.20+2.29 \times \exp \{-\exp [-0.019 \times \\ &(t-137.15)]\} \end{aligned} $$ (6) $$ \begin{aligned} & \quad\quad \lg N(t)=5.20+2.06 /\{1+\exp [-0.024 \times(t- \\ & 195.56)]\} \end{aligned} $$ (7) 2.3 细菌生长模型评价
修正Gompertz﹑Logistic模型含有4个参数,参数值和RMS值见表 2﹑表 4﹑表 5。从表 5看出,2种模型的相关系数均大于0.99, 表示2种模型均能很好拟合实验数据。0℃贮藏Gompertz﹑Logistic模型的RMS分别为0.077和0.138, 5℃贮藏Gompertz﹑Logistic模型的RMS分别为0.100和0.114, 2种预测模型相比,Gompertz的RMS均较小,其预测结果更为理想。
参数
parameters均方根(RMS)
root mean square相关系数(R)
correlation coefficient0℃ 5℃ 0℃ 5℃ Gompertz模型model 0.077 0.100 0.994 0.993 Logistic模型model 0.138 0.114 0.992 0.995 3. 讨论
(1) Whiting和Buchanan[7]依据预测模型的发展阶段分为菌数增殖变化模型、环境要素模型和专家系统模型。细菌生长曲线通常呈不对称S形,使用Logistic模型定量描述细菌生长时,不能有效表示细菌生长的延滞期,Gibson等[8, 9]修正了Logistic模型和Gompertz模型,能更好的描述细菌生长情况,然而受食品成分构成、外界环境因子影响,数学模型的适用性显现出较大差异。许多研究者[9, 10]对不同数学模型的性能和适用性进行分析与比较,结果显示相同的数学模型针对不同的研究对象预测结果存在差异。即使如此,大家普遍认为Logistic模型和Gompertz模型的预测效果较好,并得以广泛应用,例如修正的Gompertz模型被应用于预测微生物软件程序Food Micromodel和Pathogen Modeling Program等[9],Dalgaard等[11, 12]把修正的Logistic模型用于新鲜和气调包装鱼的微生物品质分析, 建立了3参数和4参数预测模型。
(2) 建立在恒温条件下的微生物生长模型,很难有效预测在实际生产、流通、贮藏、消费过程食品微生物的生长变化情况。在实际过程中温度是随机波动的,无法直接用数学式来描述时间-温度的变化规律,现在多根据实际过程把时间-温度的变化设定为多个短的假设为恒温的时间间隔,然后使用分段的数学模型来描述微生物的生长[13]。本研究建立的数学模型可以快速有效预测大黄鱼0、5℃恒温冷藏过程细菌生长变化情况,有待进一步研究建立符合实际流通过程的波动温度预测模型。
致谢: 福建宁德夏威水产食品有限公司管崇乐、上海水产大学食品学院2004级硕士研究生周彩华等对实验工作给予了帮助和协作,谨此致谢! -
表 1 0℃冷藏大黄鱼菌落总数
Table 1 Total viable counts of P.crocea stored aerobically at 0℃
lgCFU·g-1 时间/h
time实测值
observed valuesGompertz预测模型predictive model Logistic预测模型predictive model 预测值
predicted values残存值
residual values预测值
predicted values残存值
residual values0 5.20 5.20 0.00 5.21 -0.01 49 5.10 5.20 -0.10 5.24 -0.14 76 5.25 5.21 0.04 5.29 -0.04 123.5 5.34 5.42 -0.08 5.48 -0.14 167 6.05 5.95 0.10 5.89 0.16 214 6.50 6.55 -0.05 6.52 -0.02 262 6.89 6.96 -0.07 6.99 -0.10 286 7.10 7.08 0.02 7.12 -0.02 334 7.34 7.22 0.12 7.23 0.11 383 7.22 7.29 -0.07 7.27 -0.05 409 7.31 7.31 0.00 7.28 0.03 表 2 0℃冷藏大黄鱼细菌生长动力学参数
Table 2 Kinetic parameters of bacteria growth on P.crocea stored aerobically at 0℃
参数
parametersN0
(lgCFU·g-1)Nmax
(lgCFU·g-1)μmax(h-1) M(h) Ns
(lgCFU·g-1)货架期/h
shelf lifeGompertz模型model 5.20 7.33 0.013 170.88 7.31 409 Logistic模型model 5.20 7.26 0.024 195.56 注:Ns为感官终点细菌数
Note: Ns denotes population at the time of organoleptic rejection.表 3 5℃冷藏大黄鱼菌落总数
Table 3 Total viable counts of P.crocea stored aerobically at 5℃
lgCFU·g-1 时间/h
time实测值
observed valuesGompertz预测模型predictive model Logistic预测模型predictive model 预测值
predicted values残存值
residual values预测值
predicted values残存值
residual values0 5.20 5.20 -0.00 5.22 -0.02 49 5.41 5.21 0.20 5.29 0.12 76 5.24 5.29 -0.05 5.39 -0.15 123.5 5.79 5.82 -0.04 5.81 -0.02 167 6.60 6.50 0.10 6.48 0.18 216 6.90 7.03 -0.14 7.08 -0.18 264 7.40 7.30 0.10 7.30 0.10 291 7.34 7.37 -0.03 7.34 -0.00 表 4 5℃冷藏大黄鱼细菌生长动力学参数
Table 4 Kinetics parameters of bacteria growth on P.crocea stored aerobically at 5℃
参数
parametersN0
(lgCFU·g-1)Nmax
(lgCFU·g-1)μmax(h-1) M(h) Ns
(lgCFU·g-1)货架期/h shelf life Gompertz模型model 5.20 7.29 0.016 137.17 7.34 291 Logistic模型model 5.20 7.38 0.030 154.99 表 5 非线性回归方程分析
Table 5 Analysis of nonlinear estimation equation
参数
parameters均方根(RMS)
root mean square相关系数(R)
correlation coefficient0℃ 5℃ 0℃ 5℃ Gompertz模型model 0.077 0.100 0.994 0.993 Logistic模型model 0.138 0.114 0.992 0.995 -
[1] EL Marrakchi A, Bouchriti N, Hamama A, et al. Sensory, chemical and microbiological assessment of Moroccan sardin (Sardina pilchadus) stored in ice[J]. J Food Prot, 1990, 53(7): 600-605. doi: 10.4315/0362-028X-53.7.600
[2] Surette M E, Gill T A, Leblanc P J. Biochemical basis of post-mortern nucleotide catabolism in cod (Gadus morhua) and its relationship to spoilage[J]. J Agric Food Chem, 1988, 36(1): 19-22. doi: 10.1021/jf00079a005
[3] Koutsoumanis K, Giannakourou M C, Taoukis P S, et al. Application of shelf life decision system (SLDS) to marine cultured fish quality[J]. Int J Food Microbiol, 2002, 73(2-3): 375-382. doi: 10.1016/S0168-1605(01)00659-6
[4] Baranyi J, Roberts T A. A dynamic approach to predicting bacterial growth in food. [J]. Int J Food Microbiol, 1994, 23(3-4): 277-294. doi: 10.1016/0168-1605(94)90157-0
[5] Taoukis P S, Koutsoumanis K, Nychas G J E. Use of time temperature integrators and predictive modeling for shelf life control of chilled fish under dynamic storage conditions[J]. Int J Food Microbiol, 1999, 53(1): 21-31. doi: 10.1016/S0168-1605(99)00142-7
[6] Zwietering M H, Jongenburger I, Rombouts F M, et al. Modeling of the bacterial growth curve[J]. Appl Environ Microbiol, 1990, 56(6): 1875-1881. doi: 10.1128/aem.56.6.1875-1881.1990
[7] Whitng R C, Buchanan R L. Microbial modeling[J]. Food Tech, 1994, 48(6): 113-120.
[8] Gibson A M, Bratchell N, Robertrs T A. Predicting microbial growth: growth responses of Salmonella in a laboratory medium as affected by pH, sodium chloride and storage temperature[J]. Int J Food Microbiol, 1988, 6(2): 155-178. doi: 10.1016/0168-1605(88)90051-7
[9] Buchanan R C, Whiting R C, Damert C D. When is simple good enough: a comparison of Gompertz, Baranyi, and three-phase linear models for fitting bacterial growth curves[J]. Food Microbiol, 1997, 14(4): 313-326. doi: 10.1006/fmic.1997.0125
[10] Fujikawa Hiroshi, Kai A, Morozumi Satoshi. A new Logistic model for bacterial growth[J]. J Food Hyg Soc Japan, 2003, 44(3): 155-160. doi: 10.3358/shokueishi.44.155
[11] Dalgaard P. Modeling of microbial activity and prediction of shelf life of packed fresh fish[J]. Int J Food Microbiol, 1995, 26(4): 305-317. doi: 10.1016/0168-1605(94)00136-T
[12] Koutsoumanis K, Nychas G J E. Application of a systematic experimental procedure to develop a microbial model for rapid fish shelf life prediction[J]. Int J Food Microbiol, 2000, 60(2-3): 171-184. doi: 10.1016/S0168-1605(00)00309-3
[13] Fu B, Labuza T P. Predictive microbiology for monitoring spoilage of dairy products with time temperature integrators[J]. J Food Sci, 1991, 56(5): 1209-1215.