冷卻牛肉中沙門氏菌生長預(yù)測模型的建立和驗(yàn)證*

鄭婷，董鵬程，王仁歡，張一敏，梁榮蓉，毛衍偉，朱立賢，羅欣

(山東農(nóng)業(yè)大學(xué)食品科學(xué)與工程學(xué)院，山東 泰安，271018)

沙門氏菌(Salmonella)屬腸桿菌科，革蘭氏陰性，兼性厭氧菌，包括鼠傷寒沙門氏菌、豬霍亂沙門氏菌和腸炎沙門氏菌等［1］。在世界范圍內(nèi)，危害公共健康的致病菌中沙門氏菌占據(jù)重要地位［2-3］。據(jù)估計(jì)，全球每年約有2億～13億人感染沙門菌，造成約300萬人死亡［4］;美國每年因沙門氏菌感染患病約有140萬人，約400人因此而喪生［5］;在我國，細(xì)菌性食物中毒中有70%～80%是由沙門氏菌引起的，其中90%以上的沙門氏菌感染都是由肉類動(dòng)物性食品引起［6］，尤其是畜肉類及其肉制品最易感染沙門氏菌［7］。

冷卻牛肉正以滋味鮮美、口感細(xì)膩及營養(yǎng)價(jià)值高等優(yōu)點(diǎn)，代表了生鮮牛肉的生產(chǎn)和消費(fèi)的發(fā)展方向［8］。但其營養(yǎng)豐富，水分活度高，極易受微生物的污染。在正常情況下，冷卻肉始終處于0～4℃低溫環(huán)境中，使沙門氏菌的生長受到了抑制。因此持續(xù)低溫流通是維持肉品質(zhì)并保證微生物安全的最重要的方法。但我國冷鏈系統(tǒng)不完善，在其貯藏、運(yùn)輸和流通過程中出現(xiàn)控制不當(dāng)，從而使溫度經(jīng)常異常波動(dòng)，微生物會(huì)迅速繁殖，從而影響牛肉產(chǎn)品的安全，對人類健康和生命安全造成潛在威脅。

傳統(tǒng)的微生物檢測費(fèi)時(shí)費(fèi)力，而且結(jié)果往往具有滯后性，起不到預(yù)測的作用［9］。而應(yīng)用微生物預(yù)測模型可以快速有效地對重要微生物的生長、存活和死亡進(jìn)行預(yù)測，從而確保食品在生產(chǎn)、運(yùn)輸、貯存等過程中的安全和穩(wěn)定，是一種防止病原微生物對食品污染的有效的預(yù)警工具［10］。目前，已經(jīng)建立了萵苣中、新鮮豬肉濾汁中和碎牛肉中沙門氏菌的生長動(dòng)力學(xué)研究［11-13］，但是有關(guān)冷卻牛肉中沙門氏菌生長動(dòng)力學(xué)模型方面的研究未見報(bào)道。因此，本實(shí)驗(yàn)建立了不同恒溫條件下冷卻牛肉中沙門氏菌的生長模型，預(yù)測其在不同貯藏溫度下的生長動(dòng)態(tài)，并驗(yàn)證了模型的有效性，以期為在冷卻牛肉生產(chǎn)中更好地控制該菌提供有效的手段，以避免肉制品中沙門氏菌的污染和減少由于其污染而造成的食品安全問題。

1 材料和方法

1.1 材料

1.1.1 菌株

本實(shí)驗(yàn)采用4株不同血清型的沙門氏菌混合菌株，分別是鼠傷寒沙門氏菌(Salmonella typhimurium ATCC14028)，腸炎性沙門氏菌(S.enteritidis)，阿貢沙門氏菌(S.agona)和德爾卑沙門氏菌(S.derby)，除第1株為標(biāo)準(zhǔn)菌株，其他3株均為山東農(nóng)業(yè)大學(xué)畜產(chǎn)實(shí)驗(yàn)室保存菌株。

1.1.2 主要試劑和儀器

Hektoen Enteric Agar HE瓊脂，購自北京陸橋技術(shù)有限責(zé)任公司;胰蛋白胨大豆肉湯培養(yǎng)基(TSB)，購自北京陸橋技術(shù)有限責(zé)任公司;BagFilter-400拍打袋，購自法國Interscience公司;Bag Mixer-400W均質(zhì)器，購自法國Interscience公司;HP-9272恒溫培養(yǎng)箱，購自上海一恒科學(xué)儀器有限公司;Gilson P200移液器，購自濟(jì)南莊盟生物技術(shù)有限公司;G154D立式自動(dòng)壓力蒸汽滅菌器，購自廈門致微儀器有限公司;IKA VORTEX1圓周振蕩器，購自廣州艾卡儀器設(shè)備有限公司;Healsafe1200生物安全柜，購自法國力申科學(xué)儀器有限公司。

1.2 樣品處理及接種

1.2.1 樣品的準(zhǔn)備

從山東省某肉牛屠宰公司無菌操作取背最長肌，冷藏(0～4℃)條件下3 h內(nèi)帶回實(shí)驗(yàn)室，紫外燈照射30 min去除外周及筋腱后，垂直于肌纖維方向切成25 cm2大小，厚約2 cm，重約25 g的肉塊。將每份肉樣裝入均質(zhì)袋內(nèi)并密封，冷藏備用。

1.2.2 菌懸液的制備

無菌條件下，將4株不同血清型的沙門氏菌分別平板劃線后挑取單菌落，分別接入TSB液體培養(yǎng)基中，37℃培養(yǎng)24 h，得到生長至穩(wěn)定期的菌液，最終等量混合后備用。

1.2.3 接種

用無菌生理鹽水將活化后的菌液進(jìn)行梯度稀釋，制成3～4 lg(CFU/mL)的菌懸液，把肉塊懸浮于其中15 min進(jìn)行接種，置于無菌的倒扣鐵網(wǎng)上低溫下靜置2 h以去掉多余水分，然后分裝于無菌拍打袋中。

1.3 沙門氏菌生長曲線的測定

將接種后的肉樣置于 0、4、7、10、15 和 20 ℃恒溫培養(yǎng)箱中貯藏，每隔適當(dāng)時(shí)間取出樣品進(jìn)行微生物計(jì)數(shù)。在無菌條件下，每隔相應(yīng)時(shí)間取不同溫度下的每組樣品(3個(gè))相應(yīng)肉樣25 g牛肉置于均質(zhì)拍打袋中，加入225 mL 0.1%無菌蛋白胨水，以最大速度(9 h/s)充分拍打2 min。然后用無菌蛋白胨水對濾液進(jìn)行10倍的梯度稀釋，取3個(gè)適宜梯度的稀釋液1 mL涂布于HE平板，37℃需氧培養(yǎng)24 h，計(jì)數(shù)結(jié)果進(jìn)行對數(shù)轉(zhuǎn)換。

1.4 數(shù)據(jù)處理

實(shí)驗(yàn)數(shù)據(jù)采用SAS 9.0統(tǒng)計(jì)軟件對選取的適當(dāng)模型進(jìn)行擬合和回歸。

1.5 數(shù)學(xué)模型的建立和驗(yàn)證

1.5.1 沙門氏菌初級(jí)模型的建立

分別將在0、4、7、10、15 和20 ℃貯藏條件下接種的冷卻牛肉中得到的沙門氏菌生長的實(shí)驗(yàn)數(shù)據(jù)擬合方程，用修正的 Gompertz方程［14-15］或 Linear方程［16］描述不同溫度下的生長動(dòng)態(tài)。

(1)式中:t為貯藏時(shí)間，h;N(t)為t時(shí)的菌數(shù)，CFU/g;Nmax和N0分別為最大和初始菌數(shù)，CFU/g;μmax微生物生長的最大比生長速率，h-1;λ為微生物生長的延滯時(shí)間，h。

(2)式中，t為貯藏時(shí)間，h;D表示降低一個(gè)對數(shù)級(jí)所需要的時(shí)間，h。

1.5.2 溫度對沙門氏菌生長影響的動(dòng)力學(xué)模型(二級(jí)模型)的建立

溫度對沙門氏菌生長影響的動(dòng)力學(xué)模型用平方根(Ratkowsy)方程描述，是 Ratkowsky等［17］根據(jù)微生物在0～20℃溫度條件下，最大比生長速率(μmax)和遲滯期(λ)倒數(shù)的平方根與溫度之間存在線性關(guān)系，提出的一個(gè)數(shù)學(xué)模型，方程式如下:

式中T是培養(yǎng)溫度，℃;Tmin是最低生長溫度，指微生物沒有代謝活動(dòng)時(shí)的溫度，即在此溫度下最大比生長速率為0，b是方程的常數(shù)。

1.5.3 沙門氏菌生長模型的驗(yàn)證

(1)恒定溫度的驗(yàn)證。應(yīng)用建立的生長動(dòng)力學(xué)模型求得9℃和12℃貯藏時(shí)的預(yù)測值，與冷卻牛肉中9℃和12℃貯藏實(shí)驗(yàn)中實(shí)際的沙門氏菌生長數(shù)值進(jìn)行比較，采用偏差度(Bias factor，Bf)和準(zhǔn)確度(Accuracy factor，Af)來評(píng)價(jià)已經(jīng)建立的沙門氏菌生長動(dòng)力學(xué)模型的可靠性［18］。

式中N實(shí)測是實(shí)驗(yàn)實(shí)際測得的微生物數(shù)量;N預(yù)測是應(yīng)用生長動(dòng)力學(xué)模型得到的與N實(shí)測同一時(shí)間的微生物數(shù)量;n是實(shí)驗(yàn)次數(shù)。

(2)波動(dòng)溫度的驗(yàn)證。根據(jù)實(shí)際流通過程，按4℃24 h—10℃8 h—20℃1 h—4℃197 h的貯藏條件進(jìn)行波動(dòng)溫度條件下的驗(yàn)證實(shí)驗(yàn)，每隔適當(dāng)時(shí)間取出樣品進(jìn)行微生物計(jì)數(shù)。利用EGT原理［19-21］的得到的下面的方程結(jié)合恒定溫度的參數(shù)值，建立波動(dòng)溫度的沙門氏菌生長模型，方程式如下:

同樣，將波動(dòng)溫度下沙門氏菌生長預(yù)測值與實(shí)驗(yàn)值進(jìn)行比較，通過Bf和Af來評(píng)價(jià)所建模型在波動(dòng)溫度下的可靠性。

2 結(jié)果與分析

2.1 牛肉中沙門氏菌的生長動(dòng)力學(xué)模型的建立

根據(jù)沙門氏菌在恒定溫度下的計(jì)數(shù)結(jié)果，擬合Linear方程或修正的Gompertz方程，繪制了在不同溫度下的生長曲線，如圖1所示。

圖1 在0、4、7、10、15和20 ℃貯藏溫度下牛肉利用Modified Gompertz和Linear方程擬合的沙門氏菌的生長曲線Fig.1 Growth data and fitted Modified Gompertz curves and Linear curves of Salmonella on beef stored at 0，4，7，10，15 and 20 ℃

從圖1可以看出，除0℃下沙門氏菌的生長呈線性外，其他溫度下的生長曲線均呈良好的S型，并能由修正的Gompertz方程較好地?cái)M合。根據(jù)沙門氏菌生長的實(shí)驗(yàn)值，由修正的Gompertz方程或Linear方程回歸得到沙門氏菌在不同溫度下的生長動(dòng)力學(xué)模型:

Linear方程:0℃ lg N(t)=2.799 2-0.002 71×t

修正的Gompertz方程:

4℃lg N(t)=1.942 9+1.293 1×exp(-exp(0.009 96×2.718/1.293 1× (104.2-t)+1))

7℃lg N(t)=1.863 9+3.383×exp(-exp(0.019 5×2.718/3.383 0× (50.723 0-t)+1))

10℃lg N(t)=2.014 2+4.043 1×exp(-exp(0.023 7×2.718/4.043 1× (23.101 6-t)+1))

15℃lg N(t)=2.118 1+6.101 4×exp(-exp(0.143 4×2.718/6.101 4× (8.675 6-t)+1))

20℃lg N(t)=2.030 0+6.783 7×exp(-exp(0.230 2×2.718/6.783 7× (4.044 2-t)+1))

2.1.2 沙門氏菌生長動(dòng)力學(xué)模型的參數(shù)

表1是沙門氏菌在不同溫度下的生長動(dòng)力學(xué)模型的參數(shù)。

由表2可知，溫度對沙門氏菌的生長有重要影響。在0～20℃范圍內(nèi)，沙門氏菌的最大比生長速率(μmax)和遲滯期(λ)隨溫度不同而發(fā)生顯著變化，即沙門氏菌的生長速率和遲滯期與溫度密切相關(guān)。隨著溫度的逐漸升高，沙門氏菌的最大比生長速率(μmax)和最大菌落數(shù)(Nmax)不斷增大，而遲滯期(λ)不斷減小?？梢钥闯?，溫度是影響冷卻牛肉中沙門氏菌的重要環(huán)境因素，不管是生長速率還是遲滯期對溫度都有很強(qiáng)的依賴性。初始菌數(shù)(N0)比較穩(wěn)定，是因?yàn)楸緦?shí)驗(yàn)采用同一的接種過程導(dǎo)致的。用Linear線性方程來擬合0℃貯藏過程中得到的沙門氏菌生長實(shí)驗(yàn)值，R2為0.909;用修正的Gompertz方程來擬合其他溫度下的沙門氏菌生長實(shí)驗(yàn)值，R2在0.992到0.999之間，表明不同溫度下沙門氏菌的生長動(dòng)力學(xué)模型擬合度較好，能準(zhǔn)確可靠地預(yù)測冷卻牛肉中沙門氏菌的生長狀態(tài)。

表1 牛肉在0、4、7、10、15和20℃貯藏中的初始的沙門氏菌數(shù)、最大菌數(shù)、最大比生長速率、遲滯期、RMSE和R2Table 1 Initial Salmonella number(N0)，maximum cell number(Nmax)，Maximum specific growth rate(μmax)、lag phase(λ)，Mean Square Error，and R2modelled by Modified Gompertz(Linear)of beef stored at 0，4，7，10，15 and 20 ℃

2.2 溫度對沙門氏菌生長動(dòng)力學(xué)參數(shù)的影響(二級(jí)模型的建立與驗(yàn)證)

2.2.1 二級(jí)模型的建立

由上述分析可知，Gompertz方程能很好地預(yù)測在4～20℃溫度下沙門氏菌的生長，但卻無法描述溫度的變化對沙門氏菌生長的影響，而在冷卻牛肉的實(shí)際加工處理和貯藏過程中，冷藏的溫度往往不斷發(fā)生變化。因此對于溫度對微生物生長的影響選用平方根模型(Ratkowsy方程)描述。圖2表示用Ratkowsy方程描述的溫度與最大比生長速率(μmax)和溫度與延滯時(shí)間(λ)的關(guān)系。

圖2 溫度與最大比生長速率及與延滯期的關(guān)系Fig.2 Relationship between temperature and maximum specific growth rate，also with lag phase

由圖2可看出，溫度與μmax及λ間均呈現(xiàn)良好的線性關(guān)系，溫度與μmax及λ間的平方根方程分別為:兩方程的R2值分別為0.962和0.986，表明在0～20℃范圍內(nèi)溫度與及的線性關(guān)系是可信的。

2.2.2 二級(jí)模型的可靠性評(píng)價(jià)

對沙門氏菌平方根模型的可靠性進(jìn)行評(píng)價(jià)，沙門氏菌最大比生長速率平方根的殘差見圖3。

圖3 沙門氏菌生長速率平方根的殘差分析圖Fig.3 Residual graph(μobserved1/2-μpredicted1/2):comparison of predicted value and observed value

由圖3可知，沙門氏菌平方根模型的殘差值均落在±0.045范圍內(nèi)，殘差絕對值小于0.05，因此所對應(yīng)觀測值都比較可信。因此應(yīng)用平方根模型能夠很好的描述不同溫度對沙門氏菌生長的影響［22］。

表2是根據(jù)溫度與λ及μmax平方根模型得到的實(shí)驗(yàn)值和預(yù)測值，從中可以看出，沙門氏菌生長速率平方根的殘差和遲滯期倒數(shù)的平方根的殘差值均在±0.05范圍內(nèi)，因此可認(rèn)為該模型在α=0.05的水平上是可信的。

表2 溫度與延滯時(shí)間(λ)及最大比生長速率(μmax)平方根模型的殘差值Table 2 Residuals of square-root model of temperature with lag time and maximal specific growth rate

2.3 沙門氏菌初級(jí)模型的驗(yàn)證

2.3.1 恒定溫度下的驗(yàn)證

表3是冷卻牛肉在9℃和12℃溫度下沙門氏菌的實(shí)驗(yàn)值與利用模型計(jì)算得到的預(yù)測值比較，得到預(yù)測方程的偏差度(Bf)和準(zhǔn)確度(Af)。由此可知，預(yù)測值上下波動(dòng)的幅度和誤差均在10%左右，表明本研究中建立的預(yù)測模型能有效預(yù)測牛肉在不同恒定溫度條件下沙門氏菌的生長情況。

2.3.2 波動(dòng)溫度下的驗(yàn)證

根據(jù)EGT的假設(shè)原理，利用修正的Gompertz方程擬合出波動(dòng)溫度下沙門氏菌的生長模型并繪制了生長曲線，如圖4所示。

從圖4可以看出，在設(shè)計(jì)的波動(dòng)溫度條件下，沙門氏菌的預(yù)測生長曲線與實(shí)測生長曲線走勢相同，生長曲線呈“S”型，表明評(píng)價(jià)模型擬合效果較好。由表3可知，波動(dòng)溫度下的準(zhǔn)確度Af和偏差度Bf分別為1.147和0.889，預(yù)測值和實(shí)測值之間的差異即Af和預(yù)測值上下波動(dòng)的幅度Bf均在15%以內(nèi)，表明該數(shù)學(xué)模型也能有效預(yù)測波動(dòng)溫度下沙門氏菌的生長動(dòng)態(tài)。

表3 冷卻牛肉在9℃、12℃和波動(dòng)溫度下沙門氏菌生長預(yù)測值的偏差度(Bf)和準(zhǔn)確度(Af)Table 3 Bias factor，and accuracy factor of Salmonella at 9℃，12℃ and fluctuating temperature

圖4 波動(dòng)溫度下冷卻牛肉中沙門氏菌的實(shí)測值與預(yù)測值曲線Fig.4 Predicted and observed curves of TPC grown on chilled beef during fluctuated temperatures

3 討論

影響沙門氏菌生長的主要因素是外在的溫度變化，因此本文主要研究溫度的變化對該菌生長的影響。并且冷卻肉主要是在冷鏈狀態(tài)進(jìn)行運(yùn)輸和銷售，在20℃以上產(chǎn)品的流通期很短［23］，因此本實(shí)驗(yàn)選取的是實(shí)際中可能出現(xiàn)的溫度變化范圍0～20℃，此溫度范圍內(nèi)建立的沙門氏菌在冷卻牛肉中的生長預(yù)測模型具有實(shí)際意義。很多關(guān)于微生物的預(yù)測模型是在液體培養(yǎng)基中進(jìn)行研究［24-26］，因?yàn)樵谂囵B(yǎng)基上比在真實(shí)食品上更容易控制不同的影響因子，能在短時(shí)間內(nèi)獲得建模所需的大量實(shí)驗(yàn)數(shù)據(jù)，但未考慮到不同的肌肉組織對微生物生長給他影響，簡單的液體培養(yǎng)基很難對肉品進(jìn)行真實(shí)的描述，因而預(yù)測結(jié)果往往要快于實(shí)際的生長速度，有時(shí)誤差甚至高達(dá)40%［27］，而通過向原料食品上接種特定微生物來獲得數(shù)據(jù)而建立模型，雖然方法繁瑣，但所建模型能有效預(yù)測特定微生物在實(shí)際食品上的生長。因此，本實(shí)驗(yàn)將沙門氏菌接種到冷卻牛肉表面建立了沙門氏菌的生長動(dòng)力學(xué)模型，并驗(yàn)證了該模型的可靠性。

國內(nèi)外很多對于沙門氏菌的生長模型的研究采用的是單一菌種的人工污染［28］，但菌種的種類對沙門氏菌有一定影響［29］。因此采用單一菌種得到的模型，很難做到準(zhǔn)確評(píng)價(jià)肉品中微生物的真實(shí)存在狀態(tài)。此外，美國農(nóng)業(yè)部食品安全檢驗(yàn)局建議使用混合菌株或不同血清型的混合菌株進(jìn)行模擬接種，國外在進(jìn)行食品體系中生長/失活或抑菌實(shí)驗(yàn)時(shí)大多是通過接種混合菌株［30-32］。對此，本實(shí)驗(yàn)選取了4個(gè)不同血清型的菌株混合液進(jìn)行冷卻牛肉接種，以便更好地模擬沙門氏菌實(shí)際污染情況。

沙門氏菌的生長動(dòng)力學(xué)模型很多文獻(xiàn)［11-13］已有報(bào)道，多數(shù)文獻(xiàn)選用修正的Gompertz模型作為初級(jí)模型，并聯(lián)用平方根模型來描述不同溫度下沙門氏菌的生長動(dòng)態(tài)。選用修正的Gompertz模型聯(lián)用平方根方程可獲得一些重要數(shù)據(jù)，如最大比生長速率(μmax)、遲滯期(λ)及所需溫度下的函數(shù)方程，因此，本實(shí)驗(yàn)選擇修正的Gompertz模型和Linear作為初級(jí)模型，Ratkowsy平方根方程作為二級(jí)模型。應(yīng)用修正的Gompertz模型建立了4～20℃不同恒定溫度條件下冷卻牛肉中沙門氏菌生長的動(dòng)力學(xué)模型，相關(guān)系數(shù)R2均大于0.992，RMSE均小于0.12，表明該方程能很好地描述不同溫度下沙門氏菌在冷卻牛肉中的生長動(dòng)態(tài);溫度對μ和λ的影響采用平方根模型進(jìn)行描述，即，其 R2分別為 0.962 和0.942，呈現(xiàn)良好的線性關(guān)系，且殘差值在0上下浮動(dòng)。用9℃和12℃和波動(dòng)溫度下沙門氏菌生長預(yù)測值對模型進(jìn)行驗(yàn)證，得出偏差度(Bf)為 0.912、0.991和 0.889，準(zhǔn)確度(Af)為1.118、1.019和1.147，說明該模型能有效預(yù)測牛肉中沙門氏菌在4～20℃范圍內(nèi)的生長情況。

綜上所述，本研究所建立的沙門氏菌生長動(dòng)力學(xué)模型可以快速可靠地實(shí)時(shí)預(yù)測沙門氏菌在牛肉中的生長動(dòng)態(tài)，可以為冷卻牛肉中沙門氏菌污染時(shí)的預(yù)測和監(jiān)控提供有效的工具，以提高產(chǎn)品的安全性。

［1］ Min L，Malladi S，Hurd H S，et al.Salmonella spp.in lymph nodes of fed and cull cattle:Relative assessment of risk to ground beef［J］.Food Control，2015，50(2):423-434.

［2］ Rhoades J R，Duffy G，Koutsoumanis K.Prevalence and concentration of verocytotoxig-enic Escherichia coli，Salmonella enterica and Listeria monocytogenes in the beef production chain:A review［J］.Food Microbiology，2009，26(4):357-376.

［3］ Portal D A，Karras D J.Surveillance for foodborne disease outbreak in United States，2009-2010［J］.Morbidity and Mortality Weekly Report(MMWR)，2013，62(3):91-93.

［4］ YANG B W，QIAO L P，CUI S H.Serotyping，antimicrobial susceptibility and pulse field gel electrophoresis analysis of Salmonella isolates from retail foods in Henan Province，China［J］.Food Control，2013，32(1):228-235.

［5］ Van T H，Nguyen N K，Smooker P M，et al.The antibiotic resistance characteristics of non-typhoidal Salmonella enterica isolated from food-producing animals，retail meat and humans in South East Asia［J］.International Journal of Food Microbiology，2012，154(3):98-106.

［6］ 王軍，鄭增忍，王晶鈺.動(dòng)物源性食品中沙門氏菌的風(fēng)險(xiǎn)評(píng)估［J］.中國動(dòng)物檢疫，2007，24(4):23-25.

［7］ EFSA-ECDC.The European Union summary report on trends and sources of zoonoses，zoonotic agents and foodborne outbreaks in 2010［J］.EFSA Journal，2012，17(10):2 597.

［8］ 李飛燕，張一敏，羅欣，等.貯藏過程中冷卻牛肉中微生物模型的建立［J］.食品與發(fā)酵工業(yè)，2011，37(7):220-224.

［9］ McMeekin T，Olley J，Ratkowsky D，et al.Predictive microbiology theory and application:Is it all about rates?［J］.Food Control，2013，29(2):290-299.

［10］ Pin C，Avenda?o G，Cosciani C E，et al.Modelling Salmonella concentration through the pork supply chain by considering growth and survival in fluctuating conditions of temperature，pH and aw［J］.International Journal of Food Microbiology，2011，145(1):96-102.

［11］ Sant'Ana A S，F(xiàn)ranco B，Schaffner D W.Modeling the growth rate and lag time of different strains of Salmonella enterica and Listeria monocytogenes in ready-to-eat lettuce［J］.Food Microbiology，2012，30(1):267-273.

［12］ 高麗娟，甕紹蘇，劉清珺，等.鮮豬肉中沙門氏菌預(yù)測模型的建立［J］.食品研究與開發(fā)，2011，32(12):143-146.

［13］ Juneja V K，Melendres M V，Huang L，et al.Mathematical modeling of growth of Salmonella in raw ground beef under isothermal conditions from 10 to 45 ℃［J］.International Journal of Food Microbiology，2009，131(2):106-111.

［14］ 李苗云，張秋會(huì)，趙改名，等.冷卻豬肉中單增李斯特菌生長動(dòng)力學(xué)模型的建立與評(píng)價(jià)［J］.食品科學(xué)，2009，30(21):234-236.

［15］ Belda-Galbis C M，Pina-Pérez M C，Espinosa J，et al.Use of the modified Gompertz equation to assess the Stevia rebaudiana Bertoniantilisterial kinetics［J］.Food Microbiology，38(7):56-61.

［16］ Rajkowski K T.Thermal inactivation of Escherichia coli O157:H7 and Salmonella on catfish and tilapia［J］.Food Microbiology，2012，30(2):427-431.

［17］ Ratkowsky D A，Olley J，McMeekin T A.Relationship between temperature and growth rate of bacterial cultures［J］.Journal of Bacteriology，1982，149(1):1-5.

［18］ Mccconnel J A，Schaffner D W.Validation of mathematical models for Salmonella growth in raw ground beef under dynamic temperature conditions representing loss of refrigeration［J］.Journal of Food Protection，2014，77(7):1 110-1 115.

［19］ Nychas G E，Skandamis P N，Tassou C C.Meat spoilage during distribution［J］.Meat Science，2008，78(1):77-89.

［20］ HUANG L.Estimation of growth of Clostridium perfringens in cooked beef under fluctuating temperature conditions［J］.Food Microbiology，2003，20(2):549-559.

［21］ McMeekin T A，Baranyi J，Bowman J，et al.Information systems in food safety management［J］.International Journal of Food Microbiology，2006，112(3):181-194.

［22］ Robins T P，Ocio M J，Kaloti A.The effect of growth environment on the lag phase of Listeria monocytogenes［J］.International Journal of Food Microbiology，1998，44(1):83-92.

［23］ 張媛，李金寒，陳興，等.超市保障零售環(huán)節(jié)冷鏈肉類安全的研究［J］.食品工業(yè)，2014，35(4):168-170.

［24］ Lanciotti R，Sinigaglia M，Gardini F，et al.Growth/no growth interfaces of Bacillus cereus，Staphylococcus aureus and Salmonella enteritidis in model systems based on water activity，pH，temperature and ethanol concentration［J］.Food Microbiology，2001，18(1)，659-668.

［25］ Basti A A，Razavilar V.Growth response and modeling of the effects of selected factors on the time-to-detection and probability of growth initiation of Salmonella typhimurium［J］.Food Microbiology，2004，21(4):431-438.

［26］ 丁甜，董慶利，王璐，等.單增李斯特菌在營養(yǎng)肉湯中最大比生長速率的預(yù)測模型［J］.華中農(nóng)業(yè)學(xué)報(bào)，2010，29(4):522-526.

［27］ Velugoti P R，Bohra L K，Juneja V J，et al.Dynamic model for predicting growth of Salmonella spp.in ground sterile pork［J］.Food Microbiology，2011，28(7):796-803.

［28］ Wen X，Dickson J S.Survival of Campylobacter jejuni and Salmonella entericatyphimurium in vacuum-packed，moisture enhanced pork［J］.Journal of Food Protection，2012，75(3):576-580.

［29］ Margas E，Meneses N，Beatrice C P，et al.Survival and death kinetics of Salmonella strains at low relative humidity，attached to stainless steel surfaces［J］.International Journal of Food Microbiology，2014，187(18):33-40.

［30］ Juneja V K，Barron U G，Butler F，et al.Predictive thermal inactivation model for the combined effect of temperature，cinnamaldehyde and carvacrol on starvation stressed multiple Salmonella serotypes in ground chicken［J］.International Journal of Food Microbiology，2013，165(2):184-199.

［31］ M?ller C O，Aabo S，Christensen B B，et al.Effect of natural microbiota on growth of Salmonella spp.in fresh porka predictive microbiology approach［J］.Food Microbiology，2013，34(2):284-295.

［32］ HUANG L.Determination of thermal inactivation kinetics of Listeria monocytogenes in chicken meats by isothermal and dynamic methods［J］.Food Control，2013，33(2):484-488.