单核细胞增生李斯特菌(Listeria monocytogenes)简称单增李斯特菌,能在广泛的盐浓度、温度和pH值条件下生存,是全球重要食源性致病菌[1-2],易感人群为老年人、新生儿、孕妇和免疫功能低下的患者,致死率高达20%~30%[3]。李斯特菌病主要通过食用受污染的食品传播给人类,例如肉类、蔬菜、鱼类、加工食品、即食食品和乳制品等食品[4],食源性比例高达99%[5]。2011年,美国28个州的147人因食用污染的哈密瓜感染李斯特菌病,其中33名患者死亡[6];2018年,南非因猪肉香肠污染,导致1 034例病例和204例死亡,预估损失近3亿美元[7]。在中国2008—2017年部分临床病例(759例)报告中,有128人(18%)死亡,其中新生儿93人,占73%[8]。
近年来,随着居民饮食习惯的改变、冷链物流的成熟以及即食食品的消费需求增加,我国居民发生李斯特菌病的风险呈上升态势[9]。因此,我国亟需开展相关食品的单增李斯物菌的风险评估工作。已有研究发现,从不同环境中分离的致病菌在生长、毒力、生物被膜形成和耐药性等表型方面存在广泛的差异性,即微生物的异质性[10-12]。研究者们[13-15]探究了不同来源和不同基因型的副溶血性弧菌在纯培养、食品基质和人工模拟胃液环境下的生长异质性,结果表明致病基因型(tlh+/tdh+/trh-)菌株在纯培养和模拟胃液环境下表现出更大的生长异质性。微生物的生长异质性直接影响风险评估“危害特性描述”和“暴露评估”过程,影响剂量反应模型结果的可靠性,给风险评估工作带来巨大的挑战[16]。
根据菌体抗原H和鞭毛抗原O的特异性不同,L.monocytogenes可分为13个血清型(1/2a、1/2b、1/2c、3a、3b、3c、4a、4ab、4b、4c、4d、4e和7),其中1/2a、1/2b、1/2c和4b为主要血清型,占食品和临床分离株的95%以上[2,17]。1/2a和1/2b型菌株在食品中检出较高,感染后可导致腹泻、发烧、呕吐等症状[18]。在李斯特菌病的流行暴发中,4b型菌株为优势致病血清型,患者感染4b型菌株时的死亡率为26%,而感染1/2a或1/2b型菌株的患者死亡率为16%[19]。在我国,食品和临床病例中最常见的是1/2b型菌株[20]。
因此,为了构建基于“中国国情”的异质性数据库以减少异质性对风险评估的影响,本文选取了来自我国不同地区的40株L.monocytogenes,通过Biosceen C全自动微生物生长曲线分析仪测定不同来源及血清型的L.monocytogenes在37、25、10 ℃下的生长数据,建立生长预测模型拟合得到以μmax为主的生长参数,分析不同来源和血清型的菌株之间的差异,丰富具有“中国特色”的L.monocytogenes生长异质性数据库,为开展L.monocytogenes精准定量风险评估提供数据支撑。
40株L.monocytogenes由国家食品安全风险评估中心鉴定保存,具体菌株信息如表1所示。
表1 四十株单增李斯特菌菌株信息
Table 1 Information of 40 L.monocytogenes strains
编号血清型来源编号血清型来源LM14b羊杂LM214b患者LM21/2a羊肉LM221/2b患者LM31/2c羊肉LM231/2a羊肉LM41/2b患者LM241/2c素凉菜LM54b素凉菜LM251/2b患者LM61/2b患者LM261/2b患者LM71/2b患者LM271/2b案板拭子LM81/2b熟肉片LM281/2b患者LM91/2a熟羊肉LM291/2b患者LM101/2b患者LM304b患者LM111/2c羊杂LM311/2c羊杂LM124b患者LM321/2b患者LM131/2b患者LM331/2c熟肉LM141/2a羊杂LM344b熟肉LM151/2b熟肉LM351/2c熟食摊位LM161/2b熟羊杂LM361/2b患者LM171/2a患者LM371/2a熟肉LM181/2b患者LM381/2a熟羊肉LM191/2c患者LM391/2b熟肉LM201/2b患者LM401/2b患者
BHI肉汤、PALCAM培养基基础、PALCAM选择性添加剂,北京陆桥技术有限责任公司。
Bioscreen FP-1100C型全自动微生物生长曲线分析仪,芬兰Oy Growth Curves Ab公司;ZQZY-70B振荡培养箱,上海知楚仪器有限公司;GHP-9270型隔水式恒温培养箱,上海一恒科技有限公司。
1.3.1 不同温度条件下L.monocytogenes的生长实验
将L.monocytogenes从-80 ℃保存的甘油管中活化,划线至PALCAM平板,37 ℃培养48 h。挑取单菌落接种至9 mL BHI肉汤中,振荡培养,37 ℃,150 r/min,培养5~8 h,再次接种至BHI肉汤中,37 ℃振荡培养,菌株生长到达指数中期,以此为初始接种液。
参考LIU等[13]的方法,在Bioscreen C配套的100微孔板中进行梯度稀释。吸取初始接种液20 μL接种至第一孔180 μL BHI肉汤中,以此为第一稀释梯度。用移液器吹吸混匀后进行逐级稀释,依次梯度稀释到第5孔(~104 CFU/mL)BHI肉汤中,按相同方式做3个平行。Bioscreen C设置吸光度为600 nm的宽波段,温度分别为37、2、和10 ℃,时间分别为24、48、96 h,测定生长曲线,每组试验重复3次。
1.3.2 生长参数拟合
采用修正Gomperz生长模型[21],得到菌株的延滞期(λ)、最大比生长速率(μmax)、最大最小生长密度之差(ymax),计算拟合决定系数(R2)、变异系数(coefficient of variation,CV)。调整优化后的修正Gomperz生长模型表达式如公式(1)所示:
(1)
式中:A为最初的细菌接种量OD值;C为最大菌种量与初始菌种量OD值之差;μmax为最大比增长速率,OD/h;t为时间,h;λ为延滞期,h。模型评价采用R2,决定系数R2值在0~1,且越靠近1即表明预测模型的参考价值越高。
CV的计算如公式(2)所示:
变异系数
(2)
1.3.3 数据处理与统计分析
生长曲线拟合分析采用Origin Pro 2021软件,显著性差异比较采用SPSS statistics 26.0软件进行单因素方差分析(ANOVA),P<0.05,差异显著,具有统计学意义。
40株L.monocytogenes在37、25、10 ℃下的生长数据均可用修正Gomperz生长模型进行拟合,拟合决定系数R2在0.972~0.998,表明所选模型拟合效果良好。37、25、10 ℃条件下40株L.monocytogenes的μmax值分别为0.206~0.332、0.167~0.232、0.036~0.066 OD/h,结果见图1。
图1 四十株不同来源L.monocytogenes在37、25和
10 ℃下的μmax
Fig.1 Specific growth rate μmax of 40 L.monocytogenes strains from different sources at 37、25, and 10 ℃
由图1可知,随着温度降低,L.monocytogenes菌株的μmax显著减小。以37 ℃的最佳生长条件为参照,LM1~LM40菌株在37 ℃条件下的平均μmax值呈现阶梯状增长的趋势。与37 ℃条件相比,25、10 ℃条件下菌株的生长是随机的,没有固定的趋势,这与LIU等[13]的研究中的现象相同。25 ℃条件下,LM28菌株的μmax值最大,LM23菌株的μmax值最小,但在10 ℃条件下,LM23菌株的μmax值最大,LM27菌株的μmax值最小。
采用修正Gomperz生长模型对40株L.monocytogenes在不同温度下的生长数据进行拟合,得到40株不同来源的L.monocytogenes在37、25、10 ℃条件下的延滞期(λ)、最大比生长速率(μmax),计算CV,结果见表2。
表2 四十株L.monocytogenes在37、25、10 ℃下的动力学生长参数
Table 2 Kinetic growth parameters of 40 L.monocytogenes strains at 37、25, and 10 ℃
T/℃λ(延滞期)/hμmax(最大比生长速率)/(OD/h)CV/%1052.376±5.856c0.053±0.006a12.3032514.884±1.380b0.201±0.017b8.566379.895±0.756a0.288±0.030c10.528
注:同列不同字母表示存在显著性差异(P<0.05)(下同)。
由表2可知,40株L.monocytogenes的平均λ值随着温度的升高而减少,平均μmax值随着温度的升高而增大(P<0.05)。菌株的最大比生长速率的CV值代表了L.monocytogenes菌株的生长异质性。随着温度的升高,L.monocytogenes菌株的CV值先降低后增加。在25 ℃条件下,L.monocytogenes菌株的CV值最小,此时L.monocytogenes菌株对环境变化的响应小。另外,在37、10 ℃的生长环境中进行菌株选育,有助于减小菌株生长异质性对风险评估的影响。
统计临床和食品分离菌株在不同温度下的平均μmax值进行显著性差异分析,结果见图2。在37、25、10 ℃下,临床菌株的平均μmax值均大于食品分离菌株。在25 ℃下,临床和食品分离菌株的μmax值分别为(0.212±0.014)、(0.191±0.013) OD/h,存在显著性差异,其他条件下差异不显著。
图2 37、25、10 ℃下40株临床和食品分离L.monocytogenes
菌株之间的箱线图
Fig.2 Box plot between 40 clinical and food isolates of L.monocytogenes strains at 37, 25, and 10 ℃
注:*表示与临床株相比差异显著(P<0.05)。
如图3所示,随着温度升高,临床和食品分离菌株的CV值均是先减少后增大,在25 ℃时,菌株的生长异质性最小,与表2显示的整体40株菌株在不同温度下的CV值的趋势相同。临床菌株在37、25、10 ℃下的CV值均小于食品分离菌株,这表明临床菌株的生长异质性小于食品分离株。已有研究发现,与食品分离菌株相比,临床菌株往往具有更强的毒力和致病力,更好地适应宿主肠道环境[22-23]。37 ℃条件下,强致病力的临床菌株表现出较小的生长异质性,可能导致其更好地定植于肠腔,更有效地侵入宿主组织,增加感染人员健康风险。
图3 不同来源的L.monocytogenes菌株在37、25、10 ℃下的CV值
Fig.3 The CV values of L.monocytogenes strains from different sources at 37, 25, and 10 ℃
不同血清型L.monocytogenes菌株在不同温度下的μmax值呈现出不同程度的差异性,结果见表3。在不同温度条件下,不同血清型菌株之间的μmax值没有明显规律。在25 ℃下,1/2b、1/2c和4b型菌株的μmax值显著大于1/2a型菌株(P<0.05)。在10、37 ℃下,不同血清型的菌株之间的μmax值并没有显著性差异。
表3 不同血清型的L.monocytogenes菌株在37、25、10 ℃下的μmax均值 单位:OD/h
Table 3 Mean μmax values of L.monocytogenes strains of different serotypes at 37, 25, and 10 ℃
血清型10 ℃25 ℃37 ℃1/2a(n=7)0.057±0.0070.181±0.009a0.281±0.0391/2b(n=20)0.051±0.0060.207±0.016b0.289±0.0291/2c(n=7)0.055±0.0070.199±0.015b0.282±0.0454b(n=6)0.054±0.0040.208±0.014b0.274±0.045
不同血清型的L.monocytogenes菌株的致病力有所不同。绝大多数暴发或散发的侵袭型李斯特菌病病例由4b型引起,而肠胃炎型病例多由1/2a和1/2b引起,其中1/2b和4b型属于谱系I,1/2a和1/2c型属于谱系Ⅱ[24-25],而谱系I的菌株的毒力强于谱系Ⅱ内的菌株[22]。在25 ℃下,致病力更强的1/2b和4b型菌株的μmax值显著大于1/2a型菌株;谱系I内不同血清型菌株之间在不同温度下的μmax值无明显差异,谱系Ⅱ的1/2c型菌株在25 ℃时的μmax值显著大于1/2a型菌株。
不同血清型菌株在不同温度条件下的μmax的CV值如图4所示。随着温度的升高,不同血清型菌株的生长异质性都先减小后增大,在25 ℃时,菌株的生长异质性最小,与表2显示的整体40株菌株在不同温度下的CV值的趋势相同。除1/2b型菌株外,其他血清型菌株的CV在37 ℃条件下达到最大,表现出最大的生长异质性;在25 ℃时,1/2a的CV值最低;在10 ℃时,4b型菌株的生长异质性最小,表现出对较低温环境的良好适应。
图4 血清型对L.monocytogenes在37、25、10 ℃下的生长异质性的影响
Fig.4 The influence of serotype on the growth heterogeneity of L.monocytogenes at 37, 25, and 10 ℃
已有研究表明,L.monocytogenes之间的生长异质性并不能明显地与菌株来源或血清型联系起来,只有在少数情况下,临床和食物分离株之间或不同血清型之间有一些趋势[26-28]。LIANOU等[28]研究发现4b型分离菌株作为一个群体的平均生长速率与代表其他血清型的所有菌株的平均生长速率没有显著差异。等[29]探究了21株L.monocytogenes在热灭活、生长和存活方面的差异,结果表明,L.monocytogenes菌株之间在35 ℃条件下的生长速率无显著差异,而致病力强的菌株在5 ℃时表现出更大的生长速率。本研究也仅在25 ℃条件下观察到临床菌株和1/2b和4b菌株明显的生长优势。因此,L.monocytogenes在最适条件下的生长差异并不显著,而致病力强的菌株能够更好地适应复杂恶劣的环境。另外,由于不同来源菌株之间的生长差异,我国食品分离的单增李斯特菌的优势序列类型为ST9和ST8,临床的为ST87和ST1[30-31],与国外有较大差异,直接使用基于国外菌株建立的参数与模型进行相关评估和标准制定工作,存在较大的不确定。
本研究发现,37、10 ℃下,不同来源、不同致病力的L.monocytogenes之间的μmax无显著性差异,这表明在人体环境下,温度对不同致病力的L.monocytogenes的生长速率无显著影响。在25 ℃条件(常温环境)下,1/2b型菌株具有显著的生长优势,其在我国的临床和食品样本的高检出率可能与此有关[20]。另外,谱系I内的1/2b和4b菌株在25 ℃条件时表现出更佳的生长状况,这表明在常温条件下,致病力强的1/2b和4b型L.monocytogenes菌株在食品中的污染水平更高。感染李斯特菌病的关键因素之一是食物带菌量,被污染食物可能因达不到感染量而不引起疾病[3]。因此,污染水平更高的1/2b和4b型L.monocytogenes菌株进入人体后,更容易穿过消化道、血胎和血脑屏障,从而增加婴儿死亡、脑膜炎和败血症等疾病风险[32]。1/2b和4b菌株在常温下的生长优势可能是导致其具有强致病力的原因之一。
微生物生长异质性是为了适应复杂多变的生存环境,从而表现出生长行为的差异,是微生物变异和多样性的基础。L.monocytogenes作为一种高度多样化的物种,宽泛的生存环境决定了L.monocytogenes菌株在生长、毒力等表型存在不同程度的异质性。正是由于这种菌株间的生长异质性,导致剂量-反应模型结果的不准确性,直接影响了风险评估“危害特性描述”和“暴露评估”过程,从而夸大或降低了食源性致病菌的真实风险,威胁我国人民的健康安全。然而这种菌株之间的差异源于其固有属性,无法降低或消除。因此,只有将生长异质性作为微生物风险评估中筛选特征良好菌株或混合菌株的科学依据,并纳入精准剂量-反应模型及风险表征的考量范畴,才能更加科学、准确地描述食源性致病菌在食品中的存在情况,客观、精准地评估微生物对不同人群的健康风险,保障人民的生命安全。
综上,本研究对40株不同来源和不同血清型的单增李斯特菌在不同温度(10、25、37 ℃)下的生长动力学特征进行了研究。研究结果表明,不同温度下临床菌株的平均μmax值大于食品分离株,临床菌株的生长变异系数小于食品分离株,表现出较小的温度生长异质性;同时在25 ℃时,强致病力的1/2b和4b型菌株表现出了更显著的生长优势;4b型菌株对较低温度具有更好的适应性。本研究结果将为我国L.monocytogenes精准风险评估以及实验中的菌株选择提供数据支撑,为科学地进行食源性致病菌疾病防控提供依据。
[1] GANDHI M, CHIKINDAS M L.Listeria:A foodborne pathogen that knows how to survive[J].International Journal of Food Microbiology, 2007, 113(1):1-15.
[2] 马彦宁, 赵悦, 郭云昌, 等.即食食品中单核细胞增生李斯特菌的血清学分型和毒力基因分析[J].中国食品卫生杂志, 2017, 29(1):14-18.
MA Y N, ZHAO Y, GUO Y C, et al.Serological typing and virulence-associated genes analysis of Listeria monocytogenes isolated from ready-to-eat foods[J].Chinese Journal of Food Hygiene, 2017, 29(1):14-18.
[3] BUCHANAN R L, GORRIS L G M, HAYMAN M M, et al.A review of Listeria monocytogenes:An update on outbreaks, virulence, dose-response, ecology, and risk assessments[J].Food Control, 2017, 75:1-13.
[4] DUZE S T, MARIMANI M, PATEL M.Tolerance of Listeria monocytogenes to biocides used in food processing environments[J].Food Microbiology, 2021, 97:103758.
[5] SCALLAN E, GRIFFIN P M, ANGULO F J, et al.Foodborne illness acquired in the United States:Unspecified agents[J].Emerging Infectious Diseases, 2011, 17(1):16-22.
[6] MCCOLLUM J T, CRONQUIST A B, SILK B J, et al.Multistate outbreak of listeriosis associated with cantaloupe[J].The New England Journal of Medicine, 2013, 369(10):944-953.
[7] SALAMA P J, BEN EMBAREK P K, BAGARIA J, et al.Learning from Listeria:Safer food for all[J].The Lancet, 2018, 391(10137):2305-2306.
[8] CHEN S S, MENG F Z, SUN X W, et al.Epidemiology of human listeriosis in China during 2008-2017[J].Foodborne Pathogens and Disease, 2020, 17(2):119-125.
[9] 闫军. 北京市售食品中单增李斯特菌致病力及比较基因组学研究[D].大庆:黑龙江八一农垦大学, 2018.
YAN J.Study on pathogenicity and comparative genomics of Listeria monocytogenes in food sold in Beijing[D].Daqing:Heilongjiang Bayi Agricultural University, 2018.
[10] HASEGAWA A, HARA-KUDO Y, OGATA K, et al.Differences in the stress tolerances of Vibrio parahaemolyticus strains due to their source and harboring of virulence genes[J].Journal of Food Protection, 2013, 76(8):1456-1462.
[11] 郭丹凤, 张昭寰, 肖莉莉, 等.不同耐药性致病性副溶血性弧菌的生长特性比较研究[J].食品工业科技, 2014, 35(19):137-141.
GUO D F, ZHANG Z H, XIAO L L, et al.Comparison of growth characteristic of different antibiotic resistant pathogenic Vibrio parahaemolyticus[J].Science and Technology of Food Industry, 2014, 35(19):137-141.
[12] LI L, WONG H C, NONG W Y, et al.Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus[J].BMC Genomics, 2014, 15(1):1135.
[13] LIU B X, LIU H Q, PAN Y J, et al.Comparison of the effects of environmental parameters on the growth variability of Vibrio parahaemolyticus coupled with strain sources and genotypes analyses[J].Frontiers in Microbiology, 2016, 7:994.
[14] 刘冰宣. 水产品中副溶血性弧菌生长异质性的研究[D].上海:上海海洋大学, 2017.
LIU B X.Study on growth heterogeneity of Vibrio parahaemolyticus in aquatic products[D].Shanghai:Shanghai Ocean University, 2017.
[15] WANG Y M, ZHAO Y, PAN Y J, et al.Comparison on the growth variability of Vibrio parahaemolyticus coupled with strain sources and genotypes analyses in simulated gastric digestion fluids[J].Frontiers in Microbiology, 2020, 11:212.
[16] 俞文英, 张昭寰, 刘海泉, 等.食源性致病菌菌株多相异质性对微生物风险评估的影响研究进展[J].食品科学, 2019, 40(13):296-303.
YU W Y, ZHANG Z H, LIU H Q, et al.Effect of strain heterogeneity of food-borne pathogens on microbiological risk assessment:A review[J].Food Science, 2019, 40(13):296-303.
[17] JAMALI H, THONG K L.Genotypic characterization and antimicrobial resistance of Listeria monocytogenes from ready-to-eat foods[J].Food Control, 2014, 44:1-6.
[18] 周继福. 单增李斯特菌不同家系、血清型、耐药性菌株代谢组学分析研究[D].南京:南京财经大学, 2021.
ZHOU J F.Metabonomics analysis of Listeria monocytogenes strains from different families, serotypes and drug resistance[D].Nanjing:Nanjing University of Finance and Economics, 2021.
[19] GERNER-SMIDT P, ETHELBERG S, SCHIELLERUP P, et al.Invasive listeriosis in Denmark 1994-2003:A review of 299 cases with special emphasis on risk factors for mortality[J].Clinical Microbiology and Infection, 2005, 11(8):618-624.
[20] WANG Y, JIAO Y, LAN R T, et al.Characterization of Listeria monocytogenes isolated from human Listeriosis cases in China[J].Emerging Microbes &Infections, 2015, 4(8):e50.
[21] JUNEJA V K, VALENZUELA MELENDRES M, HUANG L H, et al.Modeling the effect of temperature on growth of Salmonella in chicken[J].Food Microbiology, 2007, 24(4):328-335.
[22] MAURY M M, TSAI Y H, CHARLIER C, et al.Uncovering Listeria monocytogenes hypervirulence by harnessing its biodiversity[J].Nature Genetics, 2016, 48(3):308-313.
[23] MAURY M M, BRACQ-DIEYE H, HUANG L, et al.Hypervirulent Listeria monocytogenes clones’ adaption to mammalian gut accounts for their association with dairy products[J].Nature Communications, 2019, 10(1):2488.
[24] HASEBE R, NAKAO R, OHNUMA A, et al.Listeria monocytogenes serotype 4b strains replicate in monocytes/macrophages more than the other serotypes[J].The Journal of Veterinary Medical Science, 2017, 79(6):962-969.
[25] PALUMBO J D, BORUCKI M K, MANDRELL R E, et al.Serotyping of Listeria monocytogenes by enzyme-linked immunosorbent assay and identification of mixed-serotype cultures by colony immunoblotting[J].Journal of Clinical Microbiology, 2003, 41(2):564-571.
[26] BUNCIC S, AVERY S M, ROCOURT J, et al.Can food-related environmental factors induce different behaviour in two key serovars, 4b and 1/2a, of Listeria monocytogenes?[J].International Journal of Food Microbiology, 2001, 65(3):201-212.
[27] VIALETTE M, PINON A, CHASSEIGNAUX E, et al.Growths kinetics comparison of clinical and seafood Listeria monocytogenes isolates in acid and osmotic environment[J].International Journal of Food Microbiology, 2003, 82(2):121-131.
[28] LIANOU A, STOPFORTH J D, YOON Y, et al.Growth and stress resistance variation in culture broth among Listeria monocytogenes strains of various serotypes and origins[J].Journal of Food Protection, 2006, 69(11):2640-2647.
[29] DE A J, WHITING R C.Thermal inactivation, growth, and survival studies of Listeria monocytogenes strains belonging to three distinct genotypic lineages[J].Journal of Food Protection, 2003, 66(9):1611-1617.
[30] LI W W, BAI L, FU P, et al.The epidemiology of Listeria monocytogenes in China[J].Foodborne Pathogens and Disease, 2018, 15(8):459-466.
[31] LU B H, YANG J W, GAO C Y, et al.Listeriosis cases and genetic diversity of their L.monocytogenes isolates in China, 2008-2019[J].Frontiers in Cellular and Infection Microbiology, 2021, 11:608352.
[32] VZQUEZ-BOLAND J A, WAGNER M, SCORTTI M.Why are some Listeria monocytogenes genotypes more likely to cause invasive (brain, placental) infection?[J].mBio, 2020, 11(6):e0312.