李金春等

摘 要：采用營養(yǎng)肉湯培養(yǎng)基，研究不同溫度、pH值和水分活度對金黃色葡萄球菌生長規(guī)律的影響，并構(gòu)建金黃色葡萄球菌菌落數(shù)的一級生長模型和最大比生長速率、延滯期的二級生長模型，其中一級模型通過Baranyi模型來擬合；二級模型包括單因素二級模型和多因素二級模型，分別通過B?lehrádek平方根模型和二次多項式回歸模型來擬合。結(jié)果表明：金黃色葡萄球菌生長速率隨著溫度、pH值和水分活度的增加而增加，而延滯期隨著溫度、pH值和水分活度的增加而降低，并且溫度與pH值，溫度與水分活度以及pH值與水分活度之間存在明顯的交互作用。一級Baranyi模型、二級平方根模型以及多項式回歸模型擬合性較好（R2>0.85）。通過數(shù)學(xué)檢驗參數(shù)精確因子（Af）、偏差因子（Bf）、均方根誤差（RMSE）對模型進(jìn)行數(shù)學(xué)檢驗，表明構(gòu)建的一級、二級模型適用性、可靠性和擬合性均在可接受范圍內(nèi)，在建模設(shè)定的范圍內(nèi)能預(yù)測金黃色葡萄球菌建模范圍內(nèi)的延滯期和最大比生長速率。

關(guān)鍵詞：金黃色葡萄球菌；營養(yǎng)肉湯；預(yù)測微生物學(xué)；生長模型

Abstract： Effects of temperature， pH， and water activity （aw） on the growth of Staphylococcus aureus in nutrient broth medium were studied， respectively， and the first-order kinetics models of bacterial count （CFU） and the second-order kinetics models of maximum specific growth rate （μmax） and lag time were built subsequently. Results showed that increasing temperature， pH and aw could increase the μmax but decrease the value of lag time. Colony-forming units （CFU） of S. aureus were modeled as a function of T， pH and aw， which revealed a good relationship as evaluated by the Baranyi model. The single-variable second-order models were evaluated by B?lehrádek square-root-type model with correlation coefficients more than 0.85. The multi-variable second-order models， namely polynomial logistic regression equations， were fitted using a stepwise method with correlation coefficient more than 0.85. Then， the goodness of fit of the established models was evaluated by accuracy factor （Af）， bias factor （Bf） and root mean square error （RMSE）. It was found that their applicability， reliability and goodness of fit were within acceptable ranges. The first-order and second-order models provide a useful and accurate approach to predicting the growth of S. aureus within certain limits.

Key words： Staphylococcus aureus； nutrient broth medium； predictive microbiology； growth model

中圖分類號：TS251.1 文獻(xiàn)標(biāo)志碼：A 文章編號：1001-8123（2015）03-0005-06

doi： 10.7506/rlyj1001-8123-201503002

金黃色葡萄球菌（Staphyloccocus aureus）是一種革蘭氏染色陽性球型細(xì)菌，隸屬于葡萄球菌屬，除了常引起皮膚的化膿性炎癥外，其產(chǎn)生的腸毒素可導(dǎo)致食物中毒，是近年來危害食品安全的食源性致病菌[1]。近年來，國內(nèi)外關(guān)于致病菌的生長模型、致死模型均有大量報道，致病菌主要包括沙門氏菌[2-3]、單增李斯特菌[4-6]、大腸桿菌O157∶H7[7-9]、金黃色葡萄球菌等[10-11]，其中對金黃色葡萄球菌的研究多集中于單一因素的一級和二級生長模型，多因素生長模型研究的相對較少[12-14]。利用傳統(tǒng)微生物方法對食品中金黃色葡萄球菌進(jìn)行監(jiān)測費時又費力。通過微生物預(yù)測技術(shù)對食品中的致病菌建立不同條件下的模型，然后根據(jù)模型對致病菌的生長和死亡規(guī)律進(jìn)行分析預(yù)測，可提高效率、節(jié)約人力物力。預(yù)測微生物學(xué)[15-16]是一門運用數(shù)學(xué)模型來描述和預(yù)測給定條件下食品中微生物生長、衰亡規(guī)律的新興技術(shù)，目前已在食品風(fēng)險評估和食品安全控制領(lǐng)域發(fā)揮著重要作用。

Whiting等[17]將模型劃分3 個層次，一級模型描述的是微生物的生物量隨著時間變化的函數(shù)關(guān)系，可得到延滯期、最大比生長速率以及最大菌落數(shù)，主要包括Baranyi、Logistic、Modified-Gompertz等模型[18]；二級模型描述的是一級模型中各參數(shù)（如最大比生長速率、延滯期等）與條件因素（如溫度、pH值、水分活度等）之間的函數(shù)關(guān)系，上世紀(jì)八九十年代由B?lehrádek開創(chuàng)、Ratkowsky等改進(jìn)的平方根模型在擬合不同溫度或pH值下微生物的生長參數(shù)（最大比生長速率和遲滯期）是其中的典型代表[19]。三級模型是以數(shù)學(xué)模型和參數(shù)數(shù)據(jù)庫為基礎(chǔ)編寫的計算機程序或軟件，主要包括PMP、FSP、ComBase等。

本實驗通過研究營養(yǎng)肉湯培養(yǎng)基中金黃色葡萄球菌在pH值、水分活度（aw）和溫度（T）3個因素影響下的生長規(guī)律，建立金黃色葡萄球菌菌落數(shù)的一級生長模型和最大比生長速率（μmax）、延滯期（Lag）的單因素、多因素二級生長模型，通過模型可以獲得金黃色葡萄球菌在建模范圍內(nèi)的μmax和Lag，從而預(yù)測其生長情況，為完善致病微生物預(yù)測技術(shù)提供數(shù)據(jù)支持。

1 材料與方法

1.1 菌種、材料與試劑

金黃色葡萄球菌（Staphylococcus aureus，ATCC 27217），由北京市食品質(zhì)量監(jiān)督檢驗三站提供。菌種采用磁珠法在―80 ℃低溫保藏。

營養(yǎng)肉湯（NB）、胰蛋白胨大豆肉湯（TSB）、平板

計數(shù)瓊脂（PCA） 北京陸橋技術(shù)有限責(zé)任公司；乳酸、氯化鈉 國藥集團化學(xué)試劑有限公司；菌種保藏試劑盒 友康恒業(yè)生物科技有限公司。

1.2 儀器與設(shè)備

BT224S型分析天平 賽多利斯科學(xué)儀器（北京）有限公司；F1-45 Toyo型恒溫培養(yǎng)箱 東洋化學(xué)產(chǎn)業(yè)有限公司；GI54DWS型全自動高壓滅菌器 美國致微（廈門）儀器有限公司；ESCO型二級生物安全柜 上海

金鵬分析儀器有限公司；FSP-625型勻漿機 日本Nihonseiki Kaisha公司；PHT-810型pH測量儀 廣州市毅佰科學(xué)儀器有限公司；Labmaster型水分活度儀 瑞士Novasina公司。

1.3 方法

1.3.1 菌種活化、接種與計數(shù)

菌種活化：將菌種接種到10 mL TSB增菌液中進(jìn)行活化，37 ℃、150 r/min振蕩培養(yǎng)8 h，活化至對數(shù)初期

（約105 CFU/mL）。擴大培養(yǎng)：按0.5%的比例將活化好的菌種接種到100 mL TSB增菌液中進(jìn)行擴大培養(yǎng)，30 ℃靜態(tài)過夜培養(yǎng)至對數(shù)后期（約109 CFU/mL），制得菌種原液。按照GB4789.10—2010《食品微生物檢驗 金黃色葡萄球菌檢驗》[20]進(jìn)行菌落計數(shù)。均做3 個平行，重復(fù)2 次。

1.3.2 單因素試驗

金黃色葡萄球菌的最適培養(yǎng)條件為溫度37 ℃、pH6.0、aw 0.98。選擇初始pH值（4.5、5.0、5.5、6.0、6.5）、aw（0.91、0.93、0.95、0.97、0.99）和培養(yǎng)溫度（15、20、25、30、37 ℃）3 個因素進(jìn)行單因素試驗，考察各個因素對金黃色葡萄球菌生長的影響。其中水分活度通過NaCl溶液來調(diào)節(jié)[21]，營養(yǎng)肉湯培養(yǎng)基中aw與NaCl溶液質(zhì)量濃度符合二次函數(shù)關(guān)系：aw=―1.203x2―0.3484x+0.9872，式中：x為NaCl質(zhì)量濃度/（g/mL）；pH值通過乳酸進(jìn)行調(diào)節(jié)；根據(jù)菌種原液的濃度取適量接種到營養(yǎng)肉湯培養(yǎng)基中，使初始接種量在103 CFU/mL左右，培養(yǎng)時間為24 h，在0、4、6、8、10、12、24 h時取樣檢測金黃色葡萄球菌菌落數(shù)（CFU/mL）。

1.3.3 正交試驗

在單因素試驗的基礎(chǔ)上設(shè)計L9（34）正交試驗來研究pH值、aw和溫度對金黃色葡萄球菌的綜合生長作用。正交試驗因素與水平如表1所示。培養(yǎng)時間為24 h，在0、4、6、8、10、12、24 h時取樣檢測金黃色葡萄球菌菌落數(shù)（CFU/mL）。

1.4 生長模型

1.4.1 一級Baranyi生長模型

1.5 數(shù)據(jù)分析

用Origin 8.0軟件進(jìn)行單因素試驗數(shù)據(jù)的處理、一級生長模型的構(gòu)建以及柱形圖繪制；用DPS 7.05軟件對二級生長模型進(jìn)行數(shù)據(jù)處理和模型的構(gòu)建；用PASW Statistics 18.0軟件進(jìn)行pH值、aw和溫度對金黃色葡萄球菌生長作用的顯著性分析；通過Microsoft Excel 2007軟件進(jìn)行正交試驗以及模型評價因子Af、Bf、RMSE和R2數(shù)據(jù)處理分析。

2 結(jié)果與分析

2.1 單因素試驗

2.1.1 pH值對金黃色葡萄球菌生長的影響

由圖1可知，pH值對金黃色葡萄球菌的μmax影響較大，pH值在4.0～6.0范圍內(nèi)金黃色葡萄球菌的μmax隨著pH值的增加而顯著增加（P<0.05）。而當(dāng)pH值大于6.0時，金黃色葡萄球菌的μmax趨于平穩(wěn)。Stewart等[26]研究了金黃色葡萄球菌的在不同pH值條件下的生長情況，發(fā)現(xiàn)pH值在3.0～7.0范圍內(nèi)隨著pH值的增加，μmax先增加后趨于穩(wěn)定，與本研究相一致。金黃色葡萄球菌的Lag值隨pH值的增高而降低，隨著pH值的增高，對金黃色葡萄球菌的抑制作用減弱，從而使Lag值逐漸減小。

2.1.2 aw對金黃色葡萄球菌生長的影響

金黃色葡萄球菌有高度的耐鹽性，在15g/100mL NaCl溶液中都可以生長[27]。營養(yǎng)肉湯培養(yǎng)基的aw通過NaCl溶液來調(diào)節(jié)，NaCl溶液質(zhì)量濃度越高，aw越低。金黃色葡萄球菌最適生長aw為0.99。由圖2可知，營養(yǎng)肉湯培養(yǎng)基中aw在0.91～0.97范圍內(nèi)，金黃色葡萄球菌μmax隨著aw的增加而顯著增加（P<0.05）；aw大于0.97時，μmax趨于穩(wěn)定。營養(yǎng)肉湯培養(yǎng)基中aw為0.91（NaCl溶液約為0.15 g/mL）時，金黃色葡萄球菌依然能夠生長，說明金黃色葡萄球菌能夠在低aw條件下生長。金黃色葡萄球菌的Lag值隨aw的增高而降低。

2.1.3 溫度對對金黃色葡萄球菌生長的影響

金黃色葡萄球菌的最適生長溫度為37 ℃。溫度低于最適溫度時，金黃色葡萄球菌的生長受到抑制。由圖3可知，在15～37 ℃溫度范圍內(nèi)，隨著溫度的增加金黃色葡萄球菌μmax顯著增加（P<0.05）。金黃色葡萄球菌的Lag值隨溫度的增高而降低，隨著溫度的增高，對金黃色葡萄球菌的抑制作用減弱，從而使Lag逐漸減小。

2.2 正交試驗

由表2可知，各因素不同水平的組合下，金黃色葡萄球菌的μmax均不相同，μmax的變化范圍從0.0604 h-1到0.1957 h-1，分布比較均勻，有利于二級生長模型的建立。μmax最大值出現(xiàn)在第3組試驗，當(dāng)培養(yǎng)溫度25 ℃、aw 0.99、pH 6時，μmax為0.1957 h-1。而最小μmax值出現(xiàn)在第4組試驗，為0.0404 h-1。金黃色葡萄球菌的Lag值在各因素不同水平的組合下分布比較均勻，亦有利于二級生長模型的建立。

2.3 一級Baranyi生長模型

一級模型描述的是微生物的生物量隨著時間變化的函數(shù)關(guān)系，在此基礎(chǔ)上可進(jìn)一步計算出遲滯期、生長速率以及最大菌數(shù)，Baranyi模型是一級模型中的典型代表[28]；Baranyi模型屬于機理模型，是根據(jù)微生物的生長代謝機理推導(dǎo)出的動力學(xué)方程，該模型能夠很好處理延滯生長期，能夠用來擬合完整的“S”形曲線，不僅能描述生長曲線也能描述致死曲線，而且也能用于波動的環(huán)境。pH值、aw與溫度對應(yīng)的Baranyi模型各參數(shù)如表3所示，其中決定系數(shù)R2均在0.98左右，說明Baranyi模型能夠分別擬合不同pH值、aw與溫度條件下金黃色葡萄球菌的生長規(guī)律。

單因素-二級生長模型得到的是單一因素與金黃色葡萄球菌的μmax、Lag之間的定量關(guān)系。溫度、pH值和aw與μmax、Lag之間的定量關(guān)系均符合平方根模型。擬合得到模型方程如表4所示。

2.4.2 多因素多項式回歸模型

多因素-二級生長模型描述的是溫度、pH值以及aw

數(shù)學(xué)模型通常用Af、Bf、RMSE和決定系數(shù)（R2）4 個參數(shù)來評判擬合度的優(yōu)劣。參數(shù)Af和Bf用于衡量模型的有效性，Af值越小表明模型預(yù)測值與實測值越接近，模型越精確；而Bf在0.9～1.05范圍內(nèi)的模型被認(rèn)為是“好的預(yù)測模型”，在0.7～0.9或1.06～1.15范圍的被認(rèn)為“可接受的模型”，在其他范圍內(nèi)則被認(rèn)為是“不可接受的模型[29]。決定系數(shù)R2和RMSE表示模型的可靠度，R2越接近于1，RMSE越小，擬合度越高，模型越可靠。由表5可知，一級生長模型：pH-Baranyi、aw-Baranyi以及T-Baranyi模型的R2均大于0.93，Bf值介于0.9～1.05之間，說明構(gòu)建的一級模型擬合度較高，模型比較可靠。由表6可知，二級生長模型的R2均大于0.85，Bf值除aw-Lag平方根模型（式（11））和pH值、aw、T-Lag多項式回歸模型（式（15））外，其他Bf值均介于0.9～1.05之間，說明構(gòu)建的二級模型擬合度較高，模型較可靠，可用來預(yù)測營養(yǎng)肉湯中金黃色葡萄球菌的生長規(guī)律。

3 結(jié) 論

實驗結(jié)果表明金黃色葡萄球菌生長速率隨著溫度、pH值和水分活度的增加而增加，而延滯期隨著溫度、pH值和水分活度的增加而降低，并且溫度與pH值，溫度與水分活度以及pH值與水分活度之間存在明顯的交互作用。一級Baranyi模型、二級平方根模型以及多項式回歸模型擬合性較好（R2>0.85），在建模設(shè)定的范圍內(nèi)能預(yù)測金黃色葡萄球菌建模范圍內(nèi)的延滯期和最大比生長速率。構(gòu)建模型選取的溫度、pH值和水分活度3 個因素與薩拉米香腸加工過程中發(fā)酵階段的一些重要環(huán)境因素和柵欄因子存在一一對應(yīng)關(guān)系，而各個模型均能夠很好地預(yù)測金黃色葡萄球菌在營養(yǎng)肉湯中的最大比生長速率和延滯期，這為評估薩拉米香腸生產(chǎn)體系風(fēng)險系數(shù)和制定產(chǎn)品質(zhì)量安全控制方案提供重要參考依據(jù)。

參考文獻(xiàn)：

[1] ZHANG Yanjun， ZHANG Junyan， MEI Lingling， et al. Typing and distribution of toxingenes among Staphylococcus aureus isolates from samples of milk[J]. Chinese Journal of Health Laboratory Technology， 2005， 15（6）： 682-683.

[2] THOMAS P， OSCAR I. Development and validation of a tertiary simulation model for predicting the potential growth of Salmonella typhimurium on cooked chicken[J]. International Journal of Food Microbiology， 2002， 76： 177-190.

[3] NIGHTINGALE K K， THIPPAREDDI H， PHEBUS R K， et al. Validation of a traditional italian-style salami manufacturing process for control of Salmonella and Listeria monocytogenes[J]. Journal of Food Protection， 2006， 69（4）： 794-800.

[4] JUNEJA V K， ALTUNTA? E G， AYHAN K， et al. Predictive model for the reduction of heat resistance of Listeria monocytogenes in ground beef by the combined effect of sodium chloride and apple polyphenols[J]. International Journal of Food Microbiology， 2013， 164： 54-59.

[5] LEE Y L， JUNG B S， YOON H J， et al. Predictive model for the growth kinetics of Listeria monocytogenes in raw pork meat as a function of temperature[J]. Food Control， 2014， 44： 16-21.

[6] 李苗云， 張秋會， 趙改名， 等. 冷卻豬肉中單增李斯特氏菌生長動力學(xué)模型的建立與評價[J]. 食品科學(xué)， 2009， 30（21）： 234-236.

[7] de BLACKBURN C W， CURTIS L M， HUMPHESON L， et al. Development of thermal inactivation models for Salmonella enteritidis and Escherichia coli O157：H7 with temperature， pH and NaCl as controlling factors[J]. International Journal of Food Microbiology， 1997， 38： 31-44.

[8] MCKELLAR R C， PER?Z-RODR?GUEZ F， HARRIS L J， et al. Evaluation of different approaches for modeling Escherichia coli O157：H7 survival on field lettuce[J]. International Journal of Food Microbiology， 2014， 184： 74-85.

[9] CHARIMBA G， HUGO C J， HUGO A. The growth， survival and thermal inactivation of Escherichia coli O157：H7 in a traditional South African sausage[J]. Meat Science， 2010， 85： 89-95.

[10] LE MARC Y， VAL?K L， MEDVE ?OV? A. Modelling the effect of the starter culture on the growth of Staphylococcus aureus in milk[J]. International Journal of Food Microbiology， 2009， 129： 306-311.

[11] SCHEERANS C， WICHA S G， MICHAEL J， et al. Concentration-esponse studies and modelling of the pharmacodynamics of linezolid： Staphylococcus aureus versus Enterococcus faecium[J]. International Journal of Antimicrobial Agents， 2014， 4727： 1-7.

[12] P?REZ-RODR?GUEZ F， POSADA-IZQUIERDO G D， VALERO A， et al. Modelling survival kinetics of Staphylococcus aureus and Escherichia coli O157：H7 on stainless steel surfaces soiled with different substrates under static conditions of temperature and relative humidity[J]. Food Microbiology， 2013， 33： 197-204.

[13] 黃和， 李文杰， 雷曉凌， 等. 模擬蟹肉中金黃色葡萄球菌生長預(yù)測模型的研究[J]. 食品科技， 2010， 35（10）： 54-57.

[14] BUENO E， GARC?A P， MART?NEZ B， et al. Phage inactivation of Staphylococcus aureus in fresh and hard-type cheeses[J]. International Journal of Food Microbiology， 2012， 158： 23-27.

[15] 胡潔云， 歐杰， 李柏林. 預(yù)報微生物學(xué)在食品安全評估中的作用[J]. 微生物學(xué)通報， 2009， 36（9）： 1397-1403.

[16] 李柏林， 郭劍飛， 歐杰. 預(yù)測微生物學(xué)數(shù)學(xué)建模的方法構(gòu)建[J]. 食品科學(xué)， 2004， 25（11）： 52-57.

[17] WHITING R C， BUCHANAN R L. Development of a quantitative risk assessment model for Salmonella enteritidis in pasteurized liquid eggs[J]. International Journal of Food Microbiology， 1997， 36（2）： 111-125.

[18] 劉超群， 王宏勛， 侯溫甫. 低溫肉制品微生物控制與預(yù)測模型應(yīng)用研究進(jìn)展[J]. 食品科學(xué)， 2009， 30（21）： 481-484.

[19] 董慶利， 屠康， 李保國， 等. 平方根模型在亞硝酸鹽抑菌預(yù)測中的探討研究[J]. 食品工業(yè)科技， 2009， 30（7）： 161-164.

[20] 全國文獻(xiàn)工作標(biāo)準(zhǔn)化技術(shù)委員會. GB/T 4789.10—2008 食品衛(wèi)生微生物學(xué)檢驗 金黃色葡萄球菌檢驗[S]. 北京： 中國標(biāo)準(zhǔn)出版社， 2008.

[21] RESNIK S L， CHIRIFE J. Proposed theoretical aw values at various temperatures for selected solutions to be used as reference sources in the range of microbial growth[J]. Journal of Food Protection， 1988， 51： 419-423.

[22] BARANYI J， ROBERTS T A. A dynamic approach to predicting bacterial growth in food[J]. International Journal of Food Microbiology， 1994， 23： 277-294.

[23] ZWIETERING M H， ROMBOUTS F M， VAN T RIET K. Comparison of definitions of the lag phase and the exponential phase in bacterial growth[J]. Journal Applied of Bacteriology， 1996， 81（5）： 501-508.

[24] SUTHERLAND J P， BAYLISS A J， BRAXTON D S， et al. Predictive modelling of Escherichia coli O157：H7： inclusion of carbon dioxide as a fourth factor in a pre-existing model[J]. International Journal of Food Microbiology， 1997， 37： 113-120.

[25] ROSS T. Indices for performance evaluation of predictive models in food microbiology[J]. Applied Bacteriology， 1996， 81（5）： 501-508.

[26] STEWART C M， COLE M B， LEGAN J D， et al. Staphylococcus aureus growth boundaries： moving towards mechanistic predictive models based on solute-specific effects[J]. Applied and Environmental Microbiology， 2002， 68： 1864-1871.

[27] JEONG Dowon， HAN Seulhwa， LEE Jonghoon. Safety and technological characterization of Staphylococcus equorum isolates from jeotgal， a Korean high-salt-fermented seafood， for starter development[J]. International Journal of Food Microbiology， 2014， 188： 108-115.

[28] XIONG R， XIE G， EDMONDSON A S， et al. Comparison of the Baranyi model with the modified Gompertz equation for modelling thermal inactivation of Listeria monocytogenes Scott A[J]. Food Microbiology， 1999， 16（3）： 269-279.

[29] 王虎虎， 董洋， 徐幸蓮， 等. 酸化亞氯酸鈉對雞肉致病菌的致死模型研究[J]. 中國食品學(xué)報， 2013， 13（9）： 33-38.