Quantifying variability on thermal resistance of Listeria monocytogenes

Effects of intrinsic characteristics of Salmonella enterica strains isolated from foods and humans, and their interaction with food matrices during simulated gastric conditions

The stomach's acidic pH is a crucial barrier against foodborne pathogens such as Salmonella enterica. This study investigated the survival of S. enterica under simulated oral and gastric conditions (SGC; pH 2 for 120 min) as a function of intrinsic pathogen characteristics and food matrix. Fifty-seven S. enterica strains isolated from food and human infections (previously characterized by serotype, virulotype, multi-drug resistance, isolation source, and isolation season) were subjected to SGC using water as a vehicle. Population reduction among the 57 isolates ranged from 2.7 to 4.7 log CFU, revealing that human isolates were inactivated less than food isolates (p = 0.0008). Among food strains, strains isolated during the cold season (food sampled from December to February) displayed the highest reduction (p = 0.00002). Six representatives of the 57 S. enterica strains were selected according to their virulotype and antimicrobial profile. They were further used to evaluate their survival under SGC in four food matrices (water, mango, tomato, and chicken), measuring S. enterica at 30 min intervals. The strains in chicken showed the lowest reduction and inactivation rate (1.42 ± 0.35 log CFU; 0.03 ± 0.005 min-1), followed by tomato (3.75 ± 0.57 log CFU; 0.15 ± 0.02 min-1), water (4.23 ± 0.27 log CFU; 0.17 ± 0.02 min-1), and mango (4.49 ± 0.39 log CFU; 0.17 ± 0.03 min-1). These data suggest that not all S. enterica strains have the same ability to survive under SGC, influencing the probability of arriving into the small intestine.

Comparative Analysis of Growth, Survival, and Virulence Characteristics of Listeria monocytogenes Isolated from Imported Meat

Listeria monocytogenes is an important foodborne pathogen with worldwide prevalence. Understanding the variability in the potential pathogenicity among strains of different subtypes is crucial for risk assessment. In this study, the growth, survival, and virulence characteristics of 16 L. monocytogenes strains isolated from imported meat in China (2018-2020) were investigated. The maximum specific growth rate (μmax) and lag phase (λ) were evaluated using the time-to-detection (TTD) method and the Baranyi model at different temperatures (25, 30, and 37 °C). Survival characteristics were determined by D-values and population reduction after exposure to heat (60, 62.5, and 65 °C) and acid (HCl, pH = 2.5, 3.5, and 4.5). The potential virulence was evaluated via adhesion and invasion to Caco-2 cells, motility, and lethality to Galleria mellonella. The potential pathogenicity was compared among strains of different lineages and subtypes. The results indicate that the lineage I strains exhibited a higher growth rate than the lineage II strains at three growth temperatures, particularly serotype 4b within lineage I. At all temperatures tested, serotypes 1/2a and 1/2b consistently demonstrated higher heat resistance than the other subtypes. No significant differences in the log reduction were observed between the lineage I and lineage II strains at pH 2.5, 3.5, and 4.5. However, the serotype 1/2c strains exhibited significantly low acid resistance at pH 2.5. In terms of virulence, the lineage I strains outperformed the lineage II strains. The invasion rate to Caco-2 cells and lethality to G. mellonella exhibited by the serotype 4b strains were higher than those observed in the other serotypes. This study provides meaningful insights into the growth, survival, and virulence of L. monocytogenes, offering valuable information for understanding the correlation between the pathogenicity and subtypes of L. monocytogenes.

Fate of osmotically adapted and biofilm Listeria monocytogenes cells after exposure to salt, heat, and liquid smoke, mimicking the stresses induced during the processing of hot smoked fish

The study aimed to investigate the response of osmotically adapted and detached biofilm Listeria monocytogenes cells following sequential stresses that occur during the processing of hot smoking, such as heating and smoke application. Thermal resistance of L. monocytogenes was significantly affected by previous osmotic adaptation of the cells. D60oC-values of osmotically adapted L. monocytogenes cells were significantly higher than control cells. The osmotically adapted and subsequently heat-injured cells were more resistant to PALCAM and less resistant to TSAYE with 5.00% NaCl (TSAYE/NaCl) than control cells. Detached biofilm cells were more thermotolerant and less resistant to PALCAM and TSAYE/NaCl than control cells. The sequential effect of smoking against heat-treated (60 °C, 20 min) and osmotically adapted or detached L. monocytogenes biofilm cells was investigated using two liquid smoke extracts (L9 and G6). L9 led to significantly higher reductions (>3.00-Log CFU) compared to G6. The heat-treated, detached biofilm cells revealed resistance to L9, presumably due to metabolic downregulation and physical protection by the extracellular polymeric substances (EPS). These data highlight the potential of the food industry to make informed decisions for using safe heat treatments during hot smoking to effectively inactivate L. monocytogenes and maintain rigorous environmental sanitation practices to control biofilm cells.

Effects of osmotic stress on Listeria monocytogenes ATCC 7644: persistent cells and heat resistance

Persistent bacteria are a microbial subpopulation that, exposed to bactericidal treatment, is killed at a slower rate than the rest of the population they are part of. They can be triggered either following stress or stochastically without external signals. The hallmark of persistent bacteria is the biphasic killing curve, a sign that, within a microbial population, two subpopulations are inactivated at a different rate. Furthermore, when plated into a fresh medium and in the absence of stressors, persistent bacteria typically remain in the lag phase longer before resuming active replication. This study aims to evaluate in vitro whether the formation of persistent cells in a strain of Listeria monocytogenes can be triggered by exposure to osmotic stress and if this phenomenon can increase heat resistance in the bacterial population. In a first experiment, the lag time distribution of a L. monocytogenes strain grown in a 6% NaCl broth was evaluated using the software ScanLag. A stationary phase broth culture was inoculated on agar plates placed on an office scanner inside an incubator at 37°C. The plates were scanned every 20' for 4 days and the acquired images were automatically elaborated with the aid of MatLab software in order to evaluate the appearance times of every single colony. The experiment was also carried out on a control culture obtained by growing the strain in the broth without salt. In a second experiment, the same broth cultures, after proper dilutions to rebalance NaCl concentration, were subjected to a heat treatment at 51°C and the death curves obtained were parameterized using the GinaFit system. Results showed that the lag phase of 31.40% of the salt culture colonies was long enough to suppose the formation of persistent bacteria. Analyses of the thermal survival curves showed that the shoulder and tail model was the one that best represented the inactivation trend of the salt culture, unlike the control culture, whose trend was essentially linear. Results of the present study show how exposure to salt could induce the formation of persistent bacteria in a L. monocytogenes strain. The last raises concerns as persistent cells may not only be undetected with common analytical techniques but they even show a greater heat resistance.

A Comprehensive Review of Variability in the Thermal Resistance (D-Values) of Food-Borne Pathogens-A Challenge for Thermal Validation Trials

The thermal processing of food relies heavily on determining the right time and temperature regime required to inactivate bacterial contaminants to an acceptable limit. To design a thermal processing regime with an accurate time and temperature combination, the D-values of targeted microorganisms are either referred to or estimated. The D-value is the time required at a given temperature to reduce the bacterial population by 90%. The D-value can vary depending on various factors such as the food matrix, the bacterial strain, and the conditions it has previously been exposed to; the intrinsic properties of the food (moisture, water activity, fat content, and pH); the method used to expose the microorganism to the thermal treatment either at the laboratory or commercial scale; the approach used to estimate the number of survivors; and the statistical model used for the analysis of the data. This review focused on Bacillus cereus, Cronobacter sakazakii, Escherichia coli, Listeria monocytogenes, and Clostridium perfringens owing to their pathogenicity and the availability of publications on their thermal resistance. The literature indicates a significant variation in D-values reported for the same strain, and it is concluded that when designing thermal processing regimes, the impact of multiple factors on the D-values of a specific microorganism needs to be considered. Further, owing to the complexity of the interactions involved, the effectiveness of regimes derived laboratory data must be confirmed within industrial food processing settings.

The Most Probable Curve method - A robust approach to estimate kinetic models from low plate count data resulting in reduced uncertainty

A novel method is proposed for fitting microbial inactivation models to data on liquid media: the Most Probable Curve (MPC) method. It is a multilevel model that makes a separation between the "true" microbial concentration according to the model, the "actual" concentration in the media considering chance, and the actual counts on the plate. It is based on the assumptions that stress resistance is homogeneous within a microbial population, and that there is no aggregation of microbial cells. Under these assumptions, the number of colonies in/on a plate follows a Poisson distribution with expected value depending on the proposed kinetic model, the number of dilutions and the plated volume. The novel method is compared against (non)linear regression based on a normal likelihood distribution (traditional method), Poisson regression and gamma-Poisson regression using data on the inactivation of Listeria monocytogenes. The conclusion is that the traditional method has limitations when the data includes plates with low (or zero) cell counts, which can be mitigated using more complex (discrete) likelihoods. However, Poisson regression uses an unrealistic likelihood function, making it unsuitable for survivor curves with several log-reductions. Gamma-Poisson regression uses a more realistic likelihood function, even though it is based mostly on empirical hypotheses. We conclude that the MPC method can be used reliably, especially when the data includes plates with low or zero counts. Furthermore, it generates a more realistic description of uncertainty, integrating the contribution of the plating error and reducing the uncertainty of the primary model parameters. Consequently, although it increases modelling complexity, the MPC method can be of great interest in predictive microbiology, especially in studies focused on variability analysis.

Heterogeneity in single-cell outgrowth of Listeria monocytogenes in half Fraser enrichment broth is affected by strain variability and physiological state

The behaviour of pathogens at the single-cell level can be highly variable and can thus affect the detection efficacy of enrichment-based detection methods. The outgrowth of single cells of three Listeria monocytogenes strains was monitored after fluorescence-activated single-cell sorting in non-selective brain heart infusion (BHI) broth and selective half Fraser enrichment broth (HFB) to quantify outgrowth heterogeneity and its effect on the detection probability. Single-cell heterogeneity was higher in HFB compared to non-selective BHI and heterogeneity increased further when cells were heat-stressed. The increase in heterogeneity was also strain-dependent because the fast-recovering strain Scott A showed less outgrowth heterogeneity than the slower-recovering strains EGDe and H7962. Modelling of the outgrowth kinetics during the primary enrichment demonstrated that starting at low cell concentrations could fail detection of L. monocytogenes at least partly due to cell heterogeneity. This highlights that it is important to take single-cell heterogeneity into account when optimizing enrichment formulations and procedures when L. monocytogenes contamination levels are low.

Strain Variability in Growth and Thermal Inactivation Characteristics of Listeria monocytogenes Strains after Acid Adaptation

ABSTRACT

Given the importance of strain variability to predictive microbiology and risk assessment, this study aimed to quantify the magnitude of strain variability in growth and thermal inactivation kinetics behaviors after acid adaptation. Thirty-three Listeria monocytogenes strains were exposed to acid-adapted tryptic soy broth supplemented with yeast extract (TSBYE; pH 5.5) and non-acid-adapted TSBYE (pH 7.0) for 20 h. Next, the growth parameters of these adapted and nonadapted strains that grew in nonbuffered TSBYE at 25°C were estimated. The tested strains were inactivated at 60°C in nonbuffered broth to obtain the heat resistance parameters. The results revealed that strain variability was present in the growth and thermal inactivation characteristics. The maximum specific growth rate ranged from 0.21 to 0.44 h-1 and from 0.20 to 0.45 h-1 after acid and nonacid adaptation, respectively. The lag times were from 0.69 to 2.56 h and from 0.24 to 3.36 h for acid-adapted and non-acid-adapted cells, respectively. The apparent D-values at 60°C of the pathogen ranged between 0.56 and 3.93 min and between 0.52 and 3.63 min for the presence and absence of acid adaptation condition, respectively. Acid adaptation significantly (P < 0.05) increased the magnitude of strain variability in the thermal inactivation characteristics of L. monocytogenes, with the coefficient of variation increasing to 0.17, whereas acid adaptation did not significantly (P ≥ 0.05) influence the variabilities in the growth parameters of the tested strains. Furthermore, the subsequent growth behaviors of all strains did not exhibit significant (P > 0.05) changes after exposure to acidic broth. However, the thermal resistance of most (25 of 33) of the tested strains increased (P < 0.05) after growing in acid-adapted broth. The relevant data generated in the present study can be used to describe the strain variability in predictive microbiology and to deeply understand the behavioral responses of different strains to acid adaptation.

HIGHLIGHTS

In vitro virulence potential, surface attachment and transcriptional response of sublethally injured Listeria monocytogenes following exposure to peracetic acid

The disinfectant Peracetic acid (PAA) can cause high levels of sublethal injury to L. monocytogenes. This study aims to evaluate phenotypic and transcriptional characteristics concerning surface attachment and virulence potential of sublethally injured L. monocytogenes ScottA and EGDe after exposure to 0.75 ppm PAA for 90 min at 4°C and subsequent incubation in TSBY at 4°C. Results showed that injured L. monocytogenes cells (99% of total population) were able to attach (after 2 and 24h) on stainless steel coupons at 4°C and 20°C. In vitro virulence assays using human intestinal epithelial Caco-2 cells showed that injured L. monocytogenes could invade host cells but could not proliferate intracellularly. In vitro virulence response was strain-dependent; injured ScottA was more invasive than EGDe. Assessment of PAA-injury at the transcriptional level showed upregulation of genes (motB, flaA) involved in flagellum motility and surface attachment. The transcriptional response of L. monocytogenes EGDe and ScottA was different; only injured ScottA demonstrated upregulation of the virulence genes inlA and plcA. Downregulation of the stress-related genes fri and kat, and upregulation of lmo0669 was observed in injured ScottA. The obtained results indicate that sublethally-injured L. monocytogenes cells may retain part of their virulence properties as well as their ability to adhere on food processing surfaces. Transmission to food products and introduction of these cells in the food chain is therefore a plausible scenario that is worth taking into consideration in terms of risk assessment. Importance L. monocytogenes is the causative agent of listeriosis a serious food-borne illness. Antimicrobial practices, such as disinfectants used for the elimination of this pathogen in food industry can produce a sublethally injured population fraction. Injured cells of this pathogen, that may survive an antimicrobial treatment, may pose a food safety-risk. Nevertheless, knowledge regarding how sublethal injury may impact important cellular traits and phenotypic responses of this pathogen is limited. This work suggests that sublethally injured L. monocytogenes cells maintain the virulence and surface attachment potential and highlights the importance of the occurrence of sublethally injured cells regarding food safety.

Heat Inactivation of Listeria monocytogenes on Pecans, Macadamia Nuts, and Sunflower Seeds

This project was undertaken to determine the kinetic parameters of thermal inactivation of Listeria monocytogenes on pecans, macadamia nuts, and sunflower seeds subjected to heat treatments simulating industry processes. Five strains were grown in nonselective medium, mixed, and resuspended before inoculating macadamia nuts, pecans, and sunflower seeds (6 to 9 Log CFU/g). Redried inoculated pecans and macadamia nuts were heated in an oven at a temperature range of 90 to 140°C. Unshelled sunflower seeds were heated in sunflower seed oil. The thermal inactivation was determined by measuring viable cell counts using standard microbiological methods. Average count data were fit to the log-linear model, and thermal-death kinetics were calculated. On pecans, the viable Listeria counts were reduced by 3 and 3.5 Log CFU/g after 40 min at 110°C and 8 min at 140°C, respectively. On macadamia nuts, the L. monocytogenes population was reduced by 5 Log CFU/g after 20 min at 120°C. Unshelled sunflower seeds were subjected to heat treatment via a hot-oil bath. On sunflower seeds, >7 Log CFU/g reductions were observed after 15 min at 120°C. The thermal resistance (D value) for inactivation on pecans at 140°C was 3.1 min and on macadamia nuts at 120°C was 4.4 min. The inactivation of L. monocytogenes was influenced by the kind of nut or seed. These results suggest that L. monocytogenes has a relatively high thermal tolerance. The findings from this study will contribute to the assessment of the effectiveness of heat treatment for control of this pathogen on nuts and seeds. IMPORTANCE Listeria monocytogenes is a major concern for the food industry in ready-to-eat (RTE) foods. In recent years, large-scale recalls have occurred with contaminated sunflower seeds and macadamia nuts that triggered product withdrawals. These events stress the importance of understanding Listeria's ability to survive heat treatments in these low-water activity foods. Nuts and seeds are subjected to a variety of thermal treatments typically referred as roasting. To date, no listeriosis outbreak has been linked to nuts and seeds, but the recent recognition that this pathogen can be detected in commercial products stresses the need for research on thermal treatments. The characterization of heat inactivation kinetics at temperatures typically used during roasting processes will be very beneficial for validation studies. This manuscript reports inactivation rates of L. monocytogenes strains inoculated onto macadamia nuts, sunflower seeds, and pecan halves subjected to temperatures between 90 and 140°C.

Listeria monocytogenes sensitivity to antimicrobial treatments depends on cell origin

In this study we investigated how cell origin could affect the efficacy of an antimicrobial treatment (mild heating combined with terpenoids) in Listeria monocytogenes Scott A, considering cells from: 1. single colony, 2. glycerol stock, 3. cold adapted culture, and 4. fresh culture in stationary phase. After treatment, culturability on BHI medium and viability assessed by flow cytometry were evaluated. Our results showed that the cell origin significantly impacted viability and culturability of L. monocytogenes towards antimicrobial treatment. The mild heat treatment combined or not with terpenoids mainly affected culturability rather than viability, although the culturability of cells from single colony was less impacted. Therefore, to mimic the worst scenario, these latter were selected to contaminate Gorgonzola rind and roast beef slices and we evaluated the ability of L. monocytogenes cells to recover their culturability (on ALOA agar medium) and to growth on the food matrix stored at 4 °C for 7 days. Our results suggest that only Gorgonzola rind allowed a partial recovery of the culturability of cells previously heated in presence or not of terpens. In conclusion, we found a connection between the cell history and sensitivity toward an antimicrobial treatment, underlying the importance to standardize the experimental procedures (starting from the cells to be used in the assay) in the assessment of cell sensitivity to a specific treatment. Finally, our study clearly indicated that VBNC cells can resuscitate under favorable conditions on a food matrix, becoming a threat for consumer’s health.

Evaluation of Strain Variability in Inactivation of Campylobacter jejuni in Simulated Gastric Fluid by Using Hierarchical Bayesian Modeling

This study was conducted to quantitatively evaluate the variability of stress resistance in different strains of Campylobacter jejuni and the uncertainty of such strain variability. We developed Bayesian statistical models with multilevel analysis to quantify variability within a strain, variability between different strains, and the uncertainty associated with these estimates. Furthermore, we measured the inactivation of 11 strains of C. jejuni in simulated gastric fluid with low pH, using the Weibullian survival model. The model was first developed for separate pH conditions and then analyzed over a range of pH levels. We found that the model parameters developed under separate pH conditions exhibited a clear dependence of survival on pH. In addition, the uncertainty of the variability between different strains could be described as the joint distribution of the model parameters. The latter model, including pH dependency, accurately predicted the number of surviving cells in individual as well as multiple strains. In conclusion, variabilities and uncertainties in inactivation could be simultaneously evaluated and interpreted via a probabilistic approach based on Bayesian theory. Such hierarchical Bayesian models could be useful for understanding individual-strain variability in quantitative microbial risk assessment. IMPORTANCE Since microbial strains vary in their growth and inactivation patterns in food materials, it is important to accurately predict these patterns for quantitative microbial risk assessment. However, most previous studies in this area have used highly resistant strains, which could lead to inaccurate predictions. Moreover, variability, including measurement errors and variability within a strain and between different strains, can contribute to predicted individual-level outcomes. Therefore, a multilevel framework is required to resolve these levels of variability and estimate their uncertainties. We developed a Bayesian predictive model for the survival of Campylobacter jejuni under simulated gastric conditions taking into account the variabilities and uncertainties. We demonstrated a high correspondence between predictions from the model and empirical measurements. The modeling procedure proposed in this study recommends a novel framework for predicting pathogen behavior, which can help improve quantitative microbial risk assessment during food production and distribution.

Amino acid substitutions in ribosomal protein RpsU enable switching between high fitness and multiple-stress resistance in Listeria monocytogenes

Microbial population heterogeneity contributes to differences in stress response between individual cells in a population, and can lead to the selection of genetically stable variants with increased stress resistance. We previously provided evidence that the multiple-stress resistant Listeria monocytogenes LO28 variant 15, carries a point mutation in the rpsU gene, resulting in an arginine-proline substitution in ribosomal protein RpsU (RpsU^17Arg-Pro). Here, we investigated the trade-off between general stress sigma factor SigB-mediated stress resistance and fitness in variant 15 using experimental evolution. By selecting for higher fitness in two parallel evolving cultures, we identified two evolved variants: 15EV1 and 15EV2. Whole genome sequencing and SNP analysis showed that both parallel lines mutated in the same codon in rpsU as the original mutation resulting in RpsU^17Pro-His (15EV1) and RpsU^17Pro-Thr (15EV2). Using a combined phenotyping and proteomics approach, we assessed the resistance of the evolved variants to both heat and acid stress, and found that in both lines reversion to WT-like fitness also resulted in WT-like stress sensitivity. Proteome analysis of L. monocytogenes LO28 WT, variant 15, 15EV1, and 15EV2 revealed high level expression of SigB regulon members only in variant 15, whereas protein profiles of both evolved variants were highly similar to that of the LO28 WT. Experiments with constructed RpsU^17Arg-Pro mutants in L. monocytogenes LO28 and EGDe, and RpsU^17Arg-His and RpsU^17Arg-Thr in LO28, confirmed that single amino acid substitutions in RpsU enable switching between multiple-stress resistant and high fitness states in L. monocytogenes.

Guidance on validation of lethal control measures for foodborne pathogens in foods

Food manufacturers are required to obtain scientific and technical evidence that a control measure or combination of control measures is capable of reducing a significant hazard to an acceptable level that does not pose a public health risk under normal conditions of distribution and storage. A validation study provides evidence that a control measure is capable of controlling the identified hazard under a worst-case scenario for process and product parameters tested. It also defines the critical parameters that must be controlled, monitored, and verified during processing. This review document is intended as guidance for the food industry to support appropriate validation studies, and aims to limit methodological discrepancies in validation studies that can occur among food safety professionals, consultants, and third-party laboratories. The document describes product and process factors that are essential when designing a validation study, and gives selection criteria for identifying an appropriate target pathogen or surrogate organism for a food product and process validation. Guidance is provided for approaches to evaluate available microbiological data for the target pathogen or surrogate organism in the product type of interest that can serve as part of the weight of evidence to support a validation study. The document intends to help food manufacturers, processors, and food safety professionals to better understand, plan, and perform validation studies by offering an overview of the choices and key technical elements of a validation plan, the necessary preparations including assembling the validation team and establishing prerequisite programs, and the elements of a validation report.

Quantification of Campylobacter jejuni gene expression after successive stresses mimicking poultry slaughtering steps

Broiler meat is considered as the most important source of the foodborne pathogen Campylobacter jejuni. Exposure to stress conditions encountered during the slaughtering process may induce bacterial adaptation mechanisms, and enhance or decrease pathogen resistance to subsequent stress. This adaptation may result from changes in bacterial gene expression. This study aims to accurately quantify the expression of selected C. jejuni genes after stresses inspired from the poultry slaughtering process. RT-qPCR was used to quantify gene expression of 44 genes in three strains after successive heat and cold stresses. Main results indicated that 26 genes out of 44 were differentially expressed following the successive thermal stresses. Three clusters of genes were differentially expressed according to the strain and the stress condition. Up-regulated genes mainly included genes involved in the heat shock response, whereas down-regulated genes belonged to metabolic pathways (such as lipid, amino-acid metabolisms). However, four genes were similarly overexpressed in the three strains; they might represent indicators of the thermal stress response at the species scale. Advances in the molecular understanding of the stress response of pathogenic bacteria, such as Campylobacter, in real-life processing conditions will make it possible to identify technological levers and better mitigate the microbial risk.

Inactivation and Subsequent Growth Kinetics of Listeria monocytogenes After Various Mild Bactericidal Treatments

This study was carried out to investigate the effects of mild heat, lactic acid, benzalkonium chloride and nisin treatments on the inactivation, sublethal injury, and subsequent growth of Listeria monocytogenes. Results showed that the Bigelow model successfully described the thermal inactivation kinetics, while the Log-linear model with tail consistently offered the most accurate fit to LA, BC, and nisin inactivation curves of cells. Differential plating indicated that percentage of sublethal injury for nisin treated cells was significantly higher than that for the other three treatments. Compared to non-treated cells, significant extension of lag time was observed for all treated cells. The longer exposures to heat treatment contributed to the extended lag time of the survivors. While for LA, BC and nisin treated cells, the longest lag time was not observed at the most severe treatment conditions. The correlation analysis of sublethal injury percentage on the duration of lag time revealed that only heat treatment showed the significant correlation. Overall, the lag time analysis could evaluate a wide range of bacterial injury. Lag time of treated cells was significantly influenced by stress treatments and temperatures of recovery, however, there were not any significant changes in the maximum specific growth rate between treated and non-treated cells under isothermal recovery conditions. The information generated from this study is valuable for utilizing intervention strategies in the elimination or growth inhibition of L. monocytogenes.

Assessing the survival and sublethal injury kinetics of Listeria monocytogenes under different food processing-related stresses

The foodborne pathogen L. monocytogenes can be present in food processing environments where it is exposed to various stressors. These antimicrobial factors, which aim to eliminate the pathogen, can induce sub-lethal injury to the bacterial cells. In the present study, we investigated the efficacy of different treatments (stresses) relevant to food processing and preservation as well as sanitation methods, in generating sub-lethal injury at 4 °C and 20 °C to two L. monocytogenes strains, ScottA and EGDe. Additionally, we evaluated the survival and extent of L. monocytogenes injury after exposure to commonly used disinfectants (peracetic acid and benzalkonium chloride), following habituation in nutrient-deprived, high-salinity medium. Each stress had a different impact on the survival and injury kinetics of L. monocytogenes. The highest injury levels were caused by peracetic acid which, at 4 °C, generated high populations of injured cells without loss of viability. Other injury-inducing stresses were lactic acid and heating. Long-term habituation in nutrient-limited and high salinity medium (4 °C) and subsequent exposure to disinfectants resulted in higher survival and injury in benzalkonium chloride and increased survival, yet with lower injury levels, in peracetic acid at 20 °C. Taken together, these results highlight the potential food safety risk emerging from the occurrence of injured cells by commonly used food processing methods. Consequently, in order to accurately assess the impact of an antimicrobial method, its potential of inducing sublethal injury needs to be considered along with lethality.

Heat Resistance of Listeria monocytogenes in Dairy Matrices Involved in Mozzarella di Bufala Campana PDO Cheese

The presence of Listeria monocytogenes in Mozzarella di Bufala Campana Protected Designation of Origin cheeses may depend on curd stretching conditions and post contaminations before packaging. To avoid cross-contamination, thermal treatment of water, brines and covering liquid may become necessary. The present study aimed to improve knowledge about L. monocytogenes thermal resistance focusing on the influence of some cheese making operations, namely curd stretching and heat treatment of fluids in contact with cheese after molding, in order to improve the safety of the cheese, optimize efficacy and sustainability of the processes. Moreover, the role that cheese curd stretching plays in L. monocytogenes inactivation was discussed. The 12 tested strains showed a very heterogeneous heat resistance that ranged from 7 to less than 1 Log₁₀ Cfu/mL reduction after 8 min at 60°C. D-values (decimal reduction times) and z-values (thermal resistance constant) calculated for the most heat resistant strain among 60 and 70°C were highly affected by the matrix and, in particular, heat resistance noticeably increased in drained cheese curd. As cheese curd stretching is not an isothermal process, to simulate the overall lethal effect of an industrial process a secondary model was built. The lethal effect of the process was estimated around 4 Log₁₀ reductions. The data provided may be useful for fresh pasta filata cheese producers in determining appropriate processing durations and temperatures for producing safe cheeses.

Incorporating strain variability in the design of heat treatments: A stochastic approach and a kinetic approach

For the design of thermal processes, the decimal reduction times (D-values) of target organisms can be used. However, many factors influence the D-value, like inherent organism's characteristics (strain variability), the effect of the history of the cells, as well as product factors and process factors. Strain variability is a very large contributor to the overall variation of the D-value. Hence, the overall reduction of microbial contaminants by a heat treatment is a combination of the occurrence of a strain with a certain heat resistance and its reduction given the prevailing conditions. This reduction can be determined using two approaches: a kinetic analysis based on integral equations or a stochastic approach based on Monte Carlo analysis. In this article, these two approaches are compared using as case studies the inactivation of two microorganisms: Listeria monocytogenes in a pasteurization process and the sporeformer Geobacillus stearothermophilus in a UHT process. Both approaches resulted in similar conclusions, highlighting that the strains with the highest heat resistance are determinant for the overall inactivation, even if the probability of cells having such extreme heat resistance is very low.

Variability in lag duration of Listeria monocytogenes strains in half Fraser enrichment broth after stress affects the detection efficacy using the ISO 11290-1 method

A collection of 23 Listeria monocytogenes strains of clinical and food origin was tested for their ability to recover and grow out in half Fraser enrichment broth following the ISO 11290-1:2017 protocol. Recovery of sub-lethally heat-injured cells in half Fraser broth was compared to reference cells with no stress pre-treatment. The enrichments were followed over time by plate counts and the growth parameters were estimated with the 3-phase model which described the data best. The reference cells without stress pre-treatment showed a short lag duration, which ranged from 1.4 to 2.7 h. However, significant variation in the ability to recover after 60 °C heat stress was observed among the tested strains and resulted in a lag duration from 4.7 to 15.8 h. A subset of strains was also exposed to low-temperature acid stress, and the lag duration showed to be also stress dependent. Scenario analyses and Monte Carlo simulations were carried out using the growth parameters obtained in the enrichments. This demonstrated that when starting with one cell, the detection threshold for efficient transfer of at least one cell to the secondary enrichment step, i.e. 2 log₁₀ CFU/ml, was not reached by 11 of 23 strains tested (48%) after exposure to 60 °C heat stress. Increasing the incubation time from 24 to 26 h and the transfer volume from 0.1 to 1.0 ml can increase the average probability to transfer at least one cell to the secondary enrichment step from 79.9% to 99.0%. When optimizing enrichment procedures, it is crucial to take strain variability into account as this can have a significant impact on the detection efficacy.

Variability in the heat resistance of Listeria monocytogenes under dynamic conditions can be more relevant than that evidenced by isothermal treatments

Heterogeneity in the response of microbial cells to environmental conditions is inherent to every biological system and can be very relevant for food safety, potentially being as important as intrinsic and extrinsic factors. However, previous studies analyzing variability in the microbial response to thermal treatments were limited to data obtained under isothermal conditions, whereas in the reality, environmental conditions are dynamic. In this article we analyse both empirically and through mathematical modelling the variability in the microbial response to thermal treatments under isothermal and dynamic conditions. Heat resistance was studied for four strains of Listeria monocytogenes (Scott A, CECT 4031, CECT 4032 and 12MOB052), in three different matrices (buffered peptone water, pH 7 Mcllvaine buffer and semi-skimmed milk). Under isothermal conditions, between-strain and between-media variability had no impact in the heat resistance, whereas it was very relevant for dynamic conditions. Therefore, the differences observed under dynamic conditions can be attributed to the variability in the ability for developing stress acclimation. The highest acclimation was observed in strain CECT 4031 (10-fold increase of the D-value), while the lowest acclimation was observed in strain CECT 4032 (50% increase of the D-value). Concerning the different media, acclimation was higher in buffered peptone water and semi-skimmed milk than in Mcllvaine buffer of pH 7.0. To the knowledge of the authors, this is the first research work that specifically analyses the variability of microbial adaptation processes that take place under dynamic conditions. It highlights that microbial heat resistance under dynamic conditions are sometimes determined by mechanisms that cannot be observed when cells are treated in isothermal conditions (e.g. acclimation) and can also be affected by variability. Consequently, empirical evidence on variability gathered under isothermal conditions should be extrapolated with care for dynamic conditions.

Multilevel modelling as a tool to include variability and uncertainty in quantitative microbiology and risk assessment. Thermal inactivation of Listeria monocytogenes as proof of concept

Variability is inherent in biology and also substantial for microbial populations. In the context of food safety risk assessment, it refers to differences in the response of different bacterial strains (between-strain variability) and different cells (within-strain variability) to the same condition (e.g. inactivation treatment). However, its quantification based on empirical observations and its incorporation in predictive models is a challenge for both experimental design and (statistical) analysis. In this article we propose the use of multilevel models to quantify (different levels of) variability and uncertainty and include them in the predictions. As proof of concept, we analyse the microbial inactivation of Listeria monocytogenes to thermal treatments including different levels of variability (between-strain and within-strain) and uncertainty. The relationship between the microbial count and time was expressed using a (non-linear) Weibullian model. Moreover, we defined stochastic hypotheses to describe the different types of variation at the level of the kinetic parameters, as well as in the observations (microbial counts). The model parameters (kinetic parameters and variances) are estimated using Bayesian statistics. The multilevel approach was compared against an analogous, single-level model. The multilevel methodology shrinks extreme parameter estimates towards the mean according to uncertainty, thus mitigating overfitting. In addition, this approach enables to easily incorporate different levels of variation (between-strain and/or within-strain variability and/or uncertainty) in the predictions. On the other hand, multilevel (Bayesian) models are more complex to define, implement, analyse and communicate than single-level models. Nevertheless, their ability to incorporate different sources of variability in predictions make them very suitable for Quantitative Microbial Risk Assessment.

Inactivation of Listeria monocytogenes in game meat applying sous vide cooking conditions

This study investigated the effects of sous vide cooking at temperatures between 50 °C and 60 °C on the inactivation kinetics of Listeria (L.) monocytogenes. Nutrient broth and minced game meat (Capreolus capreolus and Sus scrofa) were inoculated with three strains of L. monocytogenes and cooked under sous vide conditions (50, 55 or 60 °C for several hours). Results showed that the decimal reduction values (D-values) were largely dependent on the surrounding matrix. D-values of 125.5, 29.7 and 5.1 min were reached for BHI (brain heart infusion) at 50 °C, 55 °C and 60 °C, respectively. For roe deer, D-values of 49.2, 14.9 and 3.7 min and for wild boar, D-values of 100.2, 23.8 and 4.2 min were reached. It can be concluded that microbiologically safe cooking durations under sous-vide conditions below 60 °C should be considered individually for each meat product due to the dramatic influence of the matrix in comparison to higher temperature conditions.

Heat sensitization of hepatitis A virus and Tulane virus using grape seed extract, gingerol and curcumin

Human noroviruses (HNoV) and hepatitis A virus (HAV) are predominantly linked to foodborne outbreaks worldwide. As cell-culture systems to propagate HNoV in laboratories are not easily available, Tulane virus (TV) is used as a cultivable HNoV surrogate to determine inactivation. Heat-sensitization of HAV and TV by "generally recognized as safe'' (GRAS) substances can potentially reduce their time-temperature inactivation parameters during processing to ensure food safety. Curcumin, gingerol (from ginger), and grape seed extract (GSE) reportedly have anti-inflammatory, immune-modulating and antiviral properties. The objective of this study was to determine and compare the D-values and z-values of HAV and TV at 52-68 °C with or without curcumin (0.015 mg/ml), gingerol (0.1 mg/ml), or GSE (1 mg/ml) in 2-ml glass vials. HAV at ~7 log PFU/ml and TV at ~6 log PFU/ml were diluted in phosphate buffered saline (PBS) and added to two sets of six 2-mL sterile glass vials. One set served as the control and the second set had the three extracts individually added for thermal treatments in a circulating water bath for 0-10 min. The D-values for TV in PBS ranged from 4.55 ± 0.28 to 1.08 ± 0.16 min, and for HAV in PBS ranged from to 9.21 ± 0.24 to 0.67 ± 0.19 min at 52-68 °C. Decreased D-values (52-58 °C) for TV with curcumin ranging from 4.32 ± 0.25 to 0.62 ± 0.17 min, gingerol from 4.09 ± 0.18 to 0.72 ± 0.09 min and GSE from 3.82 ± 0.18 to 0.80 ± 0.07 min, with similar trends for HAV were observed. The linear model showed significant differences (p < 0.05) between the D-values of HAV and TV with and without plant extracts for most tested temperatures. This suggests that GRAS substances can potentially lower temperature and time regimens needed to inactivate HAV and TV.

On the use of in-silico simulations to support experimental design: A case study in microbial inactivation of foods

The mathematical models used in predictive microbiology contain parameters that must be estimated based on experimental data. Due to experimental uncertainty and variability, they cannot be known exactly and must be reported with a measure of uncertainty (usually a standard deviation). In order to increase precision (i.e. reduce the standard deviation), it is usual to add extra sampling points. However, recent studies have shown that precision can also be increased without adding extra sampling points by using Optimal Experiment Design, which applies optimization and information theory to identify the most informative experiment under a set of constraints. Nevertheless, to date, there has been scarce contributions to know a priori whether an experimental design is likely to provide the desired precision in the parameter estimates. In this article, two complementary methodologies to predict the parameter precision for a given experimental design are proposed. Both approaches are based on in silico simulations, so they can be performed before any experimental work. The first one applies Monte Carlo simulations to estimate the standard deviation of the model parameters, whereas the second one applies the properties of the Fisher Information Matrix to estimate the volume of the confidence ellipsoids. The application of these methods to a case study of dynamic microbial inactivation, showing how they can be used to compare experimental designs and assess their precision, is illustrated. The results show that, as expected, the optimal experimental design is more accurate than the uniform design with the same number of data points. Furthermore, it is demonstrated that, for some heating profiles, the uniform design does not ensure that a higher number of sampling points increases precision. Therefore, optimal experimental designs are highly recommended in predictive microbiology.

Mathematical modelling of the stress resistance induced in Listeria monocytogenes during dynamic, mild heat treatments

Modelling of stress acclimation induced by thermal inactivation of Listeria monocytogenes under dynamic conditions is analyzed in this work. A mathematical model that separates the effect of the instantaneous temperature from the one of stress acclimation, was used. The model was trained using isothermal inactivation experiments, and one biphasic dynamic treatment with a heating rate of 1 °C/min and a holding phase of 60 °C. These experiments were performed in laboratory media (Tryptic Soy Broth; TSB). The model parameters estimated through these experiments (D₅₅=12.87±0.82min, z=4.58±0.04°C, a=0.11±0.01min^-1, E=0.50±0.01°C and c=1.23±0.03) were successfully used to predict the microbial inactivation for another seven inactivation profiles, with and without a holding phase. Moreover, similar experiments were performed using milk as heating media, obtaining a good agreement between the model predictions and the empirical observations. The results of this study are compatible with the hypothesis that L. monocytogenes is able to develop a physiological response during dynamic treatments that increases its thermal resistance. Also, that the model used can be used to predict microbial inactivation of this microorganism taking into consideration stress acclimation.

Tail or artefact? Illustration of the impact that uncertainty of the serial dilution and cell enumeration methods has on microbial inactivation

The estimation of the concentration of microorganisms in a sample is crucial for food microbiology. For instance, it is essential for prevalence studies, challenge tests (growth and/or inactivation studies) or microbial risk assessment. The application of serial dilutions followed by viable counts in Petri dishes is probably the most extended experimental methodology for this purpose. However, this enumeration technique is also a source of uncertainty. In this article, the uncertainty of the serial dilution and viable count methodology related to the sampling error is analyzed, as well as the approximation of the microbial concentration by the number of colonies in a Petri dish. We analyze from a theoretical point of view (statistical analysis) the application of the binomial and Poisson models, demonstrating that the Poisson distribution increases the variance when used to model individual serial dilutions. On the other hand, the binomial model produces unbiased results. Therefore, the Poisson distribution is only applicable when it is a good approximation of the binomial distribution, so the use of the latter is recommended. The relevance of this uncertainty is demonstrated by Monte Carlo simulations of a generic microbial inactivation experiment, where the only source of uncertainty/variability considered is the one generated by serial plating and viable cell enumeration. Due to both the uncertainty of the methodology and the omission of zero-count plates because of the log-transformation, the simulated survival curve can have a tail. Therefore, this phenomenon, which is usually attributed to biological variability, can be to some extent an artefact of the experimental design and/or methodology.

Thermal Inactivation of Listeria monocytogenes in Crab Meat

Listeria monocytogenes is an important bacterial pathogen in seafood products, but limited information is currently available on the thermal resistance of relevant isolates in seafood. Thermal inactivation studies were undertaken (i) to provide much needed thermal inactivation data for L. monocytogenes in crab meat and (ii) to investigate whether tryptone soya broth (TSB) is representative of crab meat in thermal inactivation studies involving L. monocytogenes. D-values were obtained for a cocktail of two crab isolates (serotypes 1/2a and 4b) at 50, 55, and 60°C. In crab meat, D-values were 174.4, 28.2, and 1.6 min, respectively. Similar D-values of 176.4, 28.8, and 1.4 min were obtained in TSB. The corresponding z-values were 4.9°C (crab meat) and 4.8°C (TSB), respectively. The conclusions were that (i) current pasteurization conditions (e.g., 70°C for 2 min) would achieve complete destruction of any L. monocytogenes present in crab meat and (ii) TSB could be used as a model matrix for assessing the thermal inactivation of L. monocytogenes in crab meat.

Relevance of the Induced Stress Resistance When Identifying the Critical Microorganism for Microbial Risk Assessment

Decisions regarding microbial risk assessment usually have to be carried out with incomplete information. This is due to the large number of possible scenarios and the lack of specific data for the problem considered. Consequently, risk assessment studies are based on the information obtained with a small number of bacterial cells which are considered the most heat resistant and/or more capable of multiplying during storage. The identification of the most resistant strains is usually based on D and z-values, normally estimated from isothermal experiments. This procedure omits the potential effect that the shape of the dynamic thermal profile applied in industry has on the microbial inactivation. One example of such effects is stress acclimation, which is related to a physiological response of the cells during sub-lethal treatments that increases their resistance. In this article, we use a recently published mathematical model to compare the development of thermal resistance for Escherichia coli K12 MG1655 and E. coli CECT 515 using inactivation data already published for these strains. Based only on the isothermal experiments, E. coli K12 MG1655 would be identified as more resistant to the thermal treatment than the CECT 515 strain in the 50-65°C temperature range. However, we conclude that stress acclimation is strain (and/or media)-dependent; the CECT 515 strain has a higher capacity for developing a stress acclimation than K12 MG1655 (300% increase of the D-value for CECT 515, 50% for K12 MG1655). It, thus, has the potential to be more resistant to the thermal treatment than the K12 MG1655 strain for some conditions allowing acclimation. A methodology is proposed to identify for which conditions this may be the case. After calibrating the model parameters representing acclimation using real experimental data, the applicability of the proposed approach is demonstrated using numerical simulations, showing how the CECT 515 strain can be more resistant for some heating profiles. Consequently, the most resistant bacterial strain to a dynamic heating profile should not be identified based only on isothermal experiments (D- and z-value). The relevance of stress acclimation for the treatment studied should also be evaluated.

Listeria monocytogenes contamination of ready-to-eat foods and the risk for human health in the EU

Food safety criteria for Listeria monocytogenes in ready-to-eat (RTE) foods have been applied from 2006 onwards (Commission Regulation (EC) 2073/2005). Still, human invasive listeriosis was reported to increase over the period 2009-2013 in the European Union and European Economic Area (EU/EEA). Time series analysis for the 2008-2015 period in the EU/EEA indicated an increasing trend of the monthly notified incidence rate of confirmed human invasive listeriosis of the over 75 age groups and female age group between 25 and 44 years old (probably related to pregnancies). A conceptual model was used to identify factors in the food chain as potential drivers for L. monocytogenes contamination of RTE foods and listeriosis. Factors were related to the host (i. population size of the elderly and/or susceptible people; ii. underlying condition rate), the food (iii. L. monocytogenes prevalence in RTE food at retail; iv. L. monocytogenes concentration in RTE food at retail; v. storage conditions after retail; vi. consumption), the national surveillance systems (vii. improved surveillance), and/or the bacterium (viii. virulence). Factors considered likely to be responsible for the increasing trend in cases are the increased population size of the elderly and susceptible population except for the 25-44 female age group. For the increased incidence rates and cases, the likely factor is the increased proportion of susceptible persons in the age groups over 45 years old for both genders. Quantitative modelling suggests that more than 90% of invasive listeriosis is caused by ingestion of RTE food containing > 2,000 colony forming units (CFU)/g, and that one-third of cases are due to growth in the consumer phase. Awareness should be increased among stakeholders, especially in relation to susceptible risk groups. Innovative methodologies including whole genome sequencing (WGS) for strain identification and monitoring of trends are recommended.