Bioinactivation: Software for modelling dynamic microbial inactivation

Optimal experimental design (OED) for the growth rate of microbial populations. Are they really more "optimal" than uniform designs?

Secondary growth models from predictive microbiology can describe how the growth rate of microbial populations varies with environmental conditions. Because these models are built based on time and resource consuming experiments, model-based Optimal Experimental Design (OED) can be of interest to reduce the experimental load. In this study, we identify optimal experimental designs for two common models (full Ratkowsky and Cardinal Parameters Model (CPM)) for a different number of experiments (10-30). Calculations are also done fixing one or more model parameters, observing that this decision strongly affects the layout of the OED. Using in silico experiments, we conclude that OEDs are more informative than conventional (equidistant) designs with the same number of experiments. However, OEDs cluster the experiments near the growth limits (Xmin and Xmax) resulting in impractical designs with aggregated experimental runs ~10 times longer than conventional designs. To mitigate this, we propose a novel optimality criterion (i.e., the objective function) that accounts for the aggregated time. The novel criterion provides a reduction in parameter uncertainty with respect to the conventional design, without an increase in the experimental load. These results underline that an OED is only based on information theory (Fisher information), so the results can be impractical when actual experimental limitations are considered. The study also emphasizes that most OED schemes identify where to measure, but do not give an indication on how many experiments should be made. In this sense, numerical simulations can estimate the parameter uncertainty that would be obtained for a particular experimental design (OED or not). These results and methodologies (available in Open Code) can guide the design of future experiments for the development of secondary growth models.

Shelf life estimation of refrigerated vacuum packed beef accounting for uncertainty

This study estimates the shelf life of vacuum packed beef meat (three muscles: striploin (longissimus thoracis et lumborum, LTL), tenderloin (psoas major, PM) and outside chuck (trapezius thoracis, TT)) at refrigeration temperatures (0 °C-10 °C) based on modelling the growth of two relevant groups of spoilage microorganisms: lactic acid bacteria (LAB) and Enterobacteriaceae. The growth models were developed combining a two-step and a one-step approach. The primary modelling was used to identify the parameters affecting the growth kinetics, guiding the definition of secondary growth models. For LAB, the secondary model included the effect of temperature and initial pH on the specific growth rate. On the other hand, the model for Enterobacteriaceae incorporated the effect of temperature on the specific growth rate and the lag phase; as well as the effect of the initial pH on the specific growth rate, the lag phase and the initial microbial count. We did not observe any significant effect of the type of muscle on the growth kinetics. Once the equations were defined, the models were fitted to the complete dataset using a one-step approach. Model validation was carried out by cross-validation, mitigating the impact of an arbitrary division between training and validation sets. The models were used to estimate the shelf life of the product, based on the maximum admissible microbial concentration (7 log CFU/g for LAB, 5 log CFU/g for Enterobacteriaceae). Although LAB was the dominant microbiota, in several cases, both LAB and Enterobacteriaceae reached the critical concentration practically at the same time. Furthermore, in some scenarios, the end of shelf life would be determined by Enterobacteriaceae, pointing at the potential importance of non-dominant microorganisms for product spoilage. These results can aid in the implementation of effective control measures in the meat processing industry.

Dynamic Thermal Treatments in Green Coconut Water Induce Dynamic Stress Adaptation of Listeria innocua That Increases Its Thermal Resistance

The global coconut water market is projected to grow in the upcoming years, attributed to its numerous health benefits. However, due to its susceptibility to microbial contamination and the limitations of non-thermal decontamination methods, thermal treatments remain the primary approach to ensure the shelf-life stability and the microbiological safety of the product. In this study, the thermal inactivation of Listeria innocua, a Listeria monocytogenes surrogate, was evaluated in coconut water and in tryptone soy broth (TSB) under both isothermal (50-60 °C) and dynamic conditions (from 30 to 60 °C, with temperature increases of 0.5, 1 and 5 °C/min). Mathematical models were used to analyse the inactivation data. The Geeraerd model effectively described the thermal inactivation of L. innocua in both TSB and coconut water under isothermal conditions, with close agreement between experimental data and model fits. Parameter estimates and analysis revealed that acidified TSB is a suitable surrogate medium for studying the thermal inactivation of L. innocua in coconut water, despite minor differences observed in the shoulder length of inactivation curves, likely attributed to the media composition. The models fitted to the data obtained at isothermal conditions fail to predict L. innocua responses under dynamic conditions. This is attributed to the stress acclimation phenomenon that takes place under dynamic conditions, where bacterial cells adapt to initial sub-lethal treatment stages, leading to increased thermal resistance. Fitting the Bigelow model directly to dynamic data with fixed z-values reveals a three-fold increase in D-values with lower heating rates, supporting the role of stress acclimation. The findings of this study aid in designing pasteurization treatments targeting L. innocua in coconut water and enable the establishment of safe, mild heat treatments for refrigerated, high-quality coconut water.

Training in modern statistical methodologies and software tools for the definition and analysis of (stochastic) quantitative microbial risk assessment models with a comparison between the Hungarian and Spanish food supply chains

Human pathogenic Salmonella enterica strains have been infecting people since historical times. The original human pathogens, typhoid Salmonella strains (e.g. S. Typhi) played a huge role in the previous centuries but nowadays in the developed world the number of cases or outbreaks caused by these serotypes deceased due to the development of personal and public hygiene. Nowadays in these regions the animal-borne zoonotic serotypes (e.g. S. Enteritidis) became more important because of their high prevalence in intensive animal husbandry. But these bacteria can also appear in fruits and vegetables. The fellow joined the scientific work of the Polytechnic University of Cartagena, Spain about the safety of plant-based products, where he could gain experience in microbiological laboratory exercises and theoretical calculations of statistics and modelling. The activities in the laboratory were part of the research lines already established at the host institution, being based on the protocols they have already implemented. Nonetheless, the fellow had the opportunity to design his own experiment, do the experimental work required and analysed the data within the context of a qualitative microbiological risk assessment. The main focus was on the heat resistance of two strains of zoonotic Salmonella spp. at different temperatures. Experiments were done using a reference strain and an extremely resistant variant to evaluate this rare phenotype. The experiments were executed using a Mastia thermoresistometer, a device patented by the host institution that provides more control when studying thermal treatments than traditional methods. The data was analysed using the principles of predictive microbiology, using the D-value as an estimate of heat resistance that provides insight into the bacterial behaviour. For this, the fellow used the bioinactivation software, developed within the host group. Through the work and results the fellow learned the principles of quantitative microbiological risk assessment (QMRA) and predictive microbiology, which was the aim for the EU-FORA programme.

Training in tools to develop quantitative microbial risk assessment of ready-to-eat food with a comparison between the Romanian and Spanish food supply chains

The prevention and control of bacterial contamination on ready-to-eat (RTE) fresh produce is an essential task to ensure food safety. Therefore, the development of novel and effective decontamination technologies to ensure microbiological safety of fruits and vegetables has gained considerable attention and new sanitisation methods are needed. The antimicrobial activity of essential oils (EOs) is well documented, but their application in fresh produce remains a challenge due to their hydrophobic nature. Thus, nanoemulsions efficiently contribute to support the use of EOs in foods by enhancing their dispersibility, their contact area and facilitating the introduction into bacterial cells. The combination of these factors ultimately increases their antimicrobial activity. Quantitative microbial risk assessment (QMRA) is gaining more attention as an effective tool to assess and prevent potential risks associated with food-borne pathogens. In this context, the current project aims to study the effectiveness of different washing methods based on nanoemulsified EOs, comparing them against traditional methods, using a QMRA model for Escherichia coli O157:H7 on cherry tomatoes. Different simulations within a stochastic risk assessment model were implemented using the biorisk package for R, aiming to describe microbial behaviour and biological risk along the Romanian and Spanish food supply chains of RTE fresh produce. Nanoemulsions were prepared using oregano and rosemary EOs, each from Romania and Spain. The four nanoemulsions were evaluated as decontamination treatments to control the growth of E. coli O157:H7 on artificially contaminated cherry tomatoes. The decontamination treatments showed encouraging results, comparable to commonly used chlorine solutions. Therefore, oregano and rosemary nanoemulsions are promising and could be a feasible alternative for chlorine solutions in the reduction of microbiological contaminants.

Disentangling the contributions of initial heterogeneities and dynamic stress adaptation to nonlinearities in bacterial survival curves

The deviations from log-linearity that are often observed in bacterial survivor curves can be explained using different arguments, both biological and experimental. In this study, we used Bacillus subtilis as a model organism to demonstrate that the generally accepted vitalistic arguments (initial heterogeneities in the stress resistance of the cells in the population) may fail to describe microbial inactivation in some situations. In this sense, we showed how dynamic stress acclimation during an isothermal treatment provides an alternative explanation for survivor curves with an upwards curvature. We also provided an innovative experimental approach based on preadaptation experiments to evaluate which hypothesis is more suitable for the bacterial response. Furthermore, we used our experimental results to define bounds for the possible stress acclimation that may take place during dynamic treatments, concluding that the magnitude of stress acclimation may be larger for dynamic treatments than for isothermal experiments. We also evaluated the contribution of the SigB general stress response system to heat resistance by comparing the heat survival of wt and the ΔsigB mutant. Both strains survived better in 51, 52.5 and 55 °C when cells were pre-adapted at 48 °C than non-pre-adapted cells. However, ΔsigB was less resistant to heat than wt due to the missing SigB general stress system. Although these conclusions were based on B. subtilis as a model organism, this study can be the first step towards the development of a novel methodology able to estimate dynamic effects using only isothermal experiments. This would improve the models developed within the predictive microbiology community, improving our ability to predict microbial inactivation during industrial treatments, which are most often dynamic.

Training in tools to develop quantitative microbial risk assessment along the food chain of Spanish products

Food safety is a widespread challenge. Every year it is estimated that almost 1 in 10 people in the world fall ill after eating contaminated food resulting in over 400,000 deaths. The risk of outbreaks is higher when consuming ready-to-eat (RTE) products because they are eaten without a further cooking process that could inactivate pathogenic microorganisms. Hence, food processing is essential to increase the safety of RTE products. Microbiological risk assessment (MRA) integrates food science, microbiology and data science to provide a comprehensive understanding of the safety of the food system. MRA provides qualitative and/or quantitative information to decision makers, which might promote the adoption of better food practices. In this contest, this project aims to study and implement tools for quantitative microbial risk assessment (QMRA) of food products along the food chain. A common RTE product (cured ham) from Spain was used as a case study. Following, the exposure assessment model was implemented using mathematical models and statistical software to describe the microbial behaviour along the food chain. The study presents the possibility to identify the risk exposure in different scenarios (e.g. growth during different storage conditions, inactivation induced by traditional or innovative decontamination techniques), showing the flexibility of the predictive tools developed.

Microrisk Lab: An Online Freeware for Predictive Microbiology

Microrisk Lab is an R-based online modeling freeware designed to realize parameter estimation and model simulation in predictive microbiology. A total of 36 peer-reviewed models were integrated for parameter estimation (including primary models of bacterial growth/inactivation under static and nonisothermal conditions, secondary models of specific growth rate, and competition models of two-flora growth) and model simulation (including integrated models of deterministic or stochastic bacterial growth/inactivation under static and nonisothermal conditions) in Microrisk Lab. Each modeling section was designed to provide numerical and graphical results with comprehensive statistical indicators depending on the appropriate data set and/or parameter setting. In this study, six case studies were reproduced in Microrisk Lab and compared in parallel with DMFit, GInaFiT, IPMP 2013/GraphPad Prism, Bioinactivation FE, and @Risk, respectively. The estimated and simulated results demonstrated that the performance of Microrisk Lab was statistically equivalent to that of other existing modeling systems. Microrisk Lab allows for a friendly user experience when modeling microbial behaviors owing to its interactive interfaces, high integration, and interconnectivity. Users can freely access this application at https://microrisklab.shinyapps.io/english/ or https://microrisklab.shinyapps.io/chinese/.

Impact of Heating Rates on Alicyclobacillus acidoterrestris Heat Resistance under Non-Isothermal Treatments and Use of Mathematical Modelling to Optimize Orange Juice Processing

Alicyclobacillus acidoterrestris is a spoilage microorganism responsible for relevant product and economic losses in the beverage and juice industry. Spores of this microorganism can survive industrial heat treatments and cause spoilage during posterior storage. Therefore, an effective design of processing treatments requires an accurate understanding of the heat resistance of this microorganism. Considering that industrial treatments are dynamic; this understanding must include how the heat resistance of the microorganism is affected by the heating rate during the heating and cooling phases. The main objective of this study was to establish the effect of heating rates and complex thermal treatments on the inactivation kinetics of A. acidoterrestris. Isothermal experiments between 90 and 105 °C were carried out in a Mastia thermoresistometer, as well as four different dynamic treatments. Although most of the inactivation takes place during the holding phase, our results indicate the relevance of the heating phase for the effectiveness of the treatment. The thermal resistance of A. acidoterrestris is affected by the heating rate during the heating phase. Specifically, higher heating rates resulted in an increased microbial inactivation with respect to the one predicted based on isothermal experiments. These results provide novel information regarding the heat response of this microorganism, which can be valuable for the design of effective heat treatments to improve product safety and stability. Moreover, it highlights the need to incorporate experimental data based on dynamic treatments in process design, as heating rates can have a very significant effect on the thermal resistance of microorganisms.

To pool or not to pool: That is the question in microbial kinetics

Variation observed in heat inactivation of Salmonella strains (data from Combase) was characterized using multilevel modeling with two case studies. One study concerned repetitions at one temperature, the other concerned isothermal experiments at various temperatures. Multilevel models characterize variation at various levels and handle dependencies in the data. The Weibull model was applied using Bayesian regression. The research question was how parameters varied with experimental conditions and how data can best be analyzed: no pooling (each experiment analyzed separately), complete pooling (all data analyzed together) or partial pooling (connecting the experiments while allowing for variation between experiments). In the first case study, level 1 consisted of the measurements, level 2 of the group of repetitions. While variation in the initial number parameter was low (set by the researchers), the Weibull shape factor varied for each repetition from 0.58-1.44, and the rate parameter from 0.006-0.074 h. With partial pooling variation was much less, with complete pooling variation was strongly underestimated. In the second case study, level 1 consisted of the measurements, level 2 of the group of repetitions per temperature experiment, level 3 of the cluster of various temperature experiments. The research question was how temperature affected the Weibull parameters. Variation in initial numbers was low (set by the researchers), the rate parameter was obviously affected by temperature, the estimate of the shape parameter depended on how the data were analyzed. With partial pooling, and one-step global modeling with a Bigelow-type model for the rate parameter, shape parameter variation was minimal. Model comparison based on prediction capacity of the various models was explored. The probability distribution of calculated decimal reduction times was much narrower using multilevel global modeling compared to the usual single level two-step approach. Multilevel modeling of microbial heat inactivation appears to be a suitable and powerful method to characterize and quantify variation at various levels. It handles possible dependencies in the data, and yields unbiased parameter estimates. The answer on the question "to pool or not to pool" depends on the goal of modeling, but if the goal is prediction, then partial pooling using multilevel modeling is the answer, provided that the experimental data allow that.

On optimum dynamic temperature profiles for thermal inactivation kinetics determination

Determination of inactivation kinetics, associated with thermal processing of foods and obtained from dynamic temperature experiments, requires carefully designed experiments, the primary element being the selection of the appropriate temperature profile along with a carefully planned sampling schedule. In the present work, a number of different dynamic temperature profiles were investigated in terms of their ability to generate accurate kinetic parameters with low confidence intervals (CIs). Although alternative models have been also tested, our work was concentrated on thermal inactivation kinetics that could be described by the classical D-z values. A pair of D and z values was assumed, and for each temperature profile tested, concentration data at different processing times were generated through the appropriate models. Next, an error (up to ±2.5% or ±5%) was introduced on these theoretical values to generate pseudo-experimental data, and the back-calculation of the assumed kinetic parameters by non-linear regression was performed. The accuracy and the 95% CIs of the estimated kinetic parameters were evaluated; joint confidence regions were also constructed to investigate parameters correlation. The effect of temperature profile pattern, level of error, number of experimental points, and reference temperature was assessed. A stepwise increasing and a single triangle-pattern temperature profile were the best profiles among those tested. As a general observation, based on different kinetic models investigated, temperature profiles and sampling intervals that result in concentration versus time diagrams having shapes as suggested by the primary model used when isothermally applied are not considered appropriate for parameter estimation.

Evaluating the fate of Escherichia coli O157:H7 and Salmonella spp. on cucumbers

The occurrence of various foodborne disease outbreaks linked to the consumption of cucumbers worldwide in the last years raised concerns regarding the survival ability of foodborne pathogens on this food matrix. This work aimed at evaluating and quantifying the survival of Escherichia coli O157:H7 and Salmonella spp. on cucumber surfaces. Cucumbers were inoculated with a 5-strain cocktail of each microorganism and kept at 25 °C. The survival ability of two green fluorescent protein (GFP) labelled Salmonella strains inoculated individually on cucumbers was also evaluated. The inoculated areas were swabbed at different time intervals (maximum of 72 h) and cells were enumerated by plate count method (log CFU/cm²). The population of both pathogens decreased significantly on cucumber surfaces over time. E. coli O157:H7 could only be recovered up to 8 h while Salmonella spp. could be detected up to 24 h. The GFP-labelled Salmonella strains showed similar behaviour on cucumbers compared to the evaluated Salmonella cocktail. Survival kinetic parameters were estimated by fitting the Weibull model to the survival data. The data obtained in this study indicate that despite of the rapid decrease on concentrations of both pathogens evaluated on cucumbers surfaces, strategies to avoid their contamination during the supply chain as well as proper cleaning and disinfection protocols must be put forward to mitigate both E. coli O57:H7 and Salmonella on cucumbers and therefore, to decrease the exposure of consumers to microbial hazards and to avoid cross-contamination events during distribution, retail and in domestic environments.

Effect of Physical Structures of Food Matrices on Heat Resistance of Enterococcus faecium NRRL-2356 in Wheat Kernels, Flour and Dough

Nonpathogenic surrogate microorganisms, with a similar or slightly higher thermal resistance of the target pathogens, are usually recommended for validating practical pasteurization processes. The aim of this study was to explore a surrogate microorganism in wheat products by comparing the thermal resistance of three common bacteria in wheat kernels and flour. The most heat-resistant Enterococcus faecium NRRL-2356 rather than Salmonella cocktail and Escherichia coli ATCC 25922 was determined when heating at different temperature-time combinations at a fixed heating rate of 5 °C/min in a heating block system. The most heat-resistant pathogen was selected to investigate the influences of physical structures of food matrices. The results indicated that the heat resistance of E. faecium was influenced by physical structures of food matrices and reduced at wheat kernel structural conditions. The inactivation of E. faecium was better fitted in the Weibull distribution model for wheat dough structural conditions while in first-order kinetics for wheat kernel and flour structural conditions due to the changes of physical structures during heating. A better pasteurization effect could be achieved in wheat kernel structure in this study, which may provide technical support for thermal inactivation of pathogens in wheat-based food processing.

Microbiological and process variability using biological indicators of inactivation (BIIs) based on Bacillus cereus spores of food and fish-based animal by-products to evaluate microwave heating in a pilot plant

Microwave processing can be a valid alternative to conventional heating for different types of products. It enables a more efficient heat transfer in the food matrix, resulting in higher quality products. However, for many food products a uniform temperature distribution is not possible because of heterogeneities in their physical properties and non-uniformtiy in the electric field pattern. Hence, the effectiveness of microwave inactivation treatments is influenced by both intrinsic (differences between cells) and extrinsic variability (non-uniform temperature). Interpreting the results of the process and considering its impact on microbial inactivation is essential to ensure effective and efficient processing. In this work, we quantified the variability in microbial inactivation attained in a microwave pasteurization treatment with a tunnel configuration at pilot-plant scale. The configuration of the equipment makes it impossible to measure the product temperature during treatment. For that reason, variability in microbial counts was measured using Biological Inactivation Indicators (BIIs) based on spherical particles of alginate inoculated with spores of Bacillus spp. The stability of the BIIs and the uncertainty associated to them was assessed using preliminary experiments in a thermoresistometer. Then, they were introduced in the food product to analyse the microbial inactivation in different points of the products during the microwave treatment. Experiments were made in a vegetable soup and a fish-based animal by-product (F-BP). The results show that the variation in the microbial counts was higher than expected based on the biological variability estimated in the thermoresistometer and the uncertainty of the BIIs. This is due to heterogeneities in the temperature field (measured using a thermographic camera), which were higher in the F-BP than in the vegetable soup. Therefore, for the process studied, extrinsic variability was more relevant than intrinsic variability. The methodology presented in this work can be a valid method to evaluate pasteurization treatments of foods processed by heating, providing valuable information of the microbial inactivation achieved. It can contribute to design microwave processes for different types of products and for product optimization.

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.

Effect of moderate thermal treatments on the inactivation of a strain of Listeria monocytogenes and physicochemical properties of soursop pulp

Soursop (Annona muricata L.) is a commercially important tropical fruit, whether fresh or processed as a pasteurized or frozen pulp used to prepare juice, drinks, nectar, ice cream, popsicles, and desserts. Besides preserving quality, another preoccupation in the processing of fruit pulps is product safety. Several studies show the association between pulp processing and the development of various microorganisms; however, few have focused on the association between L. monocytogenes and the pulp of sour fruits. The objective was to evaluate the effect of moderate thermal treatments on the inactivation of L. monocytogenes and the physicochemical properties in soursop pulp in order to determine the best processing conditions that will allow to maintain quality as well as to achieve an adequate level of safety. Thermal inactivation kinetics were obtained for L. monocytogenes inoculated in soursop pulp at five levels of temperature (50, 52.5, 55, 57.5, and 60 ℃) and different exposure times (0-60 min). The survival curves did not suggest a log-linear relationship, and were, consequently, fitted to the modified Gompertz equation. The results indicated that the modified Gompertz equation provided an acceptable goodness of fit. Five-log₁₀ cycles reductions of L. monocytogenes were achieved at 50 ℃/60 min, 52.5 ℃/16 min, 55 ℃/10 min, 57.5 ℃/5 min, and 60 ℃/1.25 min. These 5-log₁₀ treatments applied to the soursop pulp indicated that the soursop pulp showed changes in the color parameters and a decrease in the content of total sugars, reducing sugars, ascorbic acid, total phenols, and pH.

Limonene nanoemulsified with soya lecithin reduces the intensity of non-isothermal treatments for inactivation of Listeria monocytogenes

Consumers' demands for ready-to-eat, fresh-like products are on the rise during the last years. This type of products have minimal processing conditions that can enable the survival and replication of pathogenic microorganisms. Among them, Listeria monocytogenes is of special concern, due to its relatively high mortality rate and its ability to replicate under refrigeration conditions. Previous research works have shown that nanoemulsified essential oils in combination with thermal treatments are effective for inactivating L. monocytogenes. However, previous research works were limited to isothermal conditions, whereas actual processing conditions in industry are dynamic. Under dynamic conditions, microorganism can respond unexpectedly to the thermal stress (e.g. adaptation, acclimation or increased sensitivity). In this work, we assess the combination of nanoemulsified D-limonene with thermal treatments under isothermal and dynamic conditions. The nanoemulsion was prepared following an innovative methodology using soya lecithin, a natural compound as well as the essential oil. Under isothermal heating conditions, the addition of the antimicrobial enables a reduction of the treatment time by a factor of 25. For time-varying treatments, dynamic effects were relevant. Treatments with a high heating rate (20 °C/min) are more effective than those with a slow heating rate (1 °C/min). This investigation demonstrates that the addition of nanoemulsified D-limonene can greatly reduce the intensity of the thermal treatments currently applied in the food industry. Hence, it can improve the product quality without impacting its safety.

The Complex Effect of Food Matrix Fat Content on Thermal Inactivation of Listeria monocytogenes: Case Study in Emulsion and Gelled Emulsion Model Systems

Previous studies on the influence of food matrix fat content on thermal inactivation kinetics of food pathogens have shown contradictory results due to the combined influence of fat content and other factors such as composition. Therefore, thermal inactivation of Listeria monocytogenes at 59, 64, and 69°C was systematically investigated in emulsion and gelled emulsion food model systems with various fat content (1, 5, 10, and 20%), such that the effect of fat content was isolated. Thermal conductivity and rheological properties of the model systems were quantified, as well as the effect of these properties on the thermal load of the model systems. Thermal conductivity was complexly related to fat content, the nature of the food matrix (i.e., viscous or gelled), and temperature. For the emulsions, the consistency index K increased with increasing fat content, while the flow behavior index n followed the opposite trend. For the gelled emulsions, the storage modulus G' was always larger than the loss modulus G″ (i.e., measure of elastic and viscous properties, respectively). The phase angle δ [i.e., arctan (G″/G')] was proportional with fat content, but this relation became more complex at higher temperatures. The thermal load of the model systems was not largely affected by food matrix fat content. Thermal inactivation of L. monocytogenes was investigated by means of the maximum specific inactivation rate k max, log reductions, and sublethal injury (SI). Both for emulsions and gelled emulsions, k max decreased with increasing fat content below approximately 60°C, while a more complex behavior was observed at higher temperatures. In the emulsions, log reductions were considerably lower (i.e., 2-3 log) at 1% fat than in systems with higher fat content. In the gelled emulsions, log reductions generally decreased with increasing fat content. SI decreased with increasing fat content, both in emulsions and gelled emulsions. In conclusion, the inactivation rate (i.e., k max) of L. monocytogenes was affected by a complex relation between food matrix fat content, thermal conductivity, rheological properties, and inactivation temperature. Due to the small scale of the model systems, differences in k max did not directly affect the final log reductions in a similar fashion.

Guidelines for the design of (optimal) isothermal inactivation experiments

Kinetic models are nowadays a basic tool to ensure food safety. Most models used in predictive microbiology have model parameters, whose precision is crucial to provide meaningful predictions. Kinetic parameters are usually estimated based on experimental data, where the experimental design can have a great impact on the precision of the estimates. In this sense, Optimal Experiment Design (OED) applies tools from optimization and information theory to identify the most informative experiment under a set of constrains (e.g. mathematical model, number of samples, etc). In this work, we develop a methodology for the design of optimal isothermal inactivation experiments. We consider the two dimensions of the design space (time and temperature), as well as a temperature-dependent maximum duration of the experiment. Functions for its application have been included in the bioOED R package. We identify design patterns that remain optimum regardless of the number of sampling points for three inactivation models (Bigelow, Mafart and Peleg) and three model microorganisms (Escherichia coli, Salmonella Senftenberg and Bacillus coagulans). Samples at extreme temperatures and close to the maximum duration of the experiment are the most informative. Moreover, the Mafart and Peleg models require some samples at intermediate time points due to the non-linearity of the survivor curve. The impact of the reference temperature on the precision of the parameter estimates is also analysed. Based on numerical simulations we recommend fixing it to the mean of the maximum and minimum temperatures used for the experiments. The article ends with a discussion presenting guidelines for the design of isothermal inactivation experiments. They combine these optimum results based on information theory with several practical limitations related to isothermal inactivation experiments. The application of these guidelines would reduce the experimental burden required to characterize thermal inactivation.

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.

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.

Optimization of E. coli Inactivation by Benzalkonium Chloride Reveals the Importance of Quantifying the Inoculum Effect on Chemical Disinfection

Optimal disinfection protocols are fundamental to minimize bacterial resistance to the compound applied, or cross-resistance to other antimicrobials such as antibiotics. The objective is twofold: guarantee safe levels of pathogens and minimize the excess of disinfectant after a treatment. In this work, the disinfectant dose is optimized based on a mathematical model. The model explains and predicts the interplay between disinfectant and pathogen at different initial microbial densities (inocula) and dose concentrations. The study focuses on the disinfection of Escherichia coli with benzalkonium chloride, the most common quaternary ammonium compound. Interestingly, the specific benzalkonium chloride uptake (mean uptake per cell) decreases exponentially when the inoculum concentration increases. As a consequence, the optimal disinfectant dose increases exponentially with the initial bacterial concentration.

Bioinactivation FE: A free web application for modelling isothermal and dynamic microbial inactivation

Mathematical models developed in predictive microbiology are nowadays an essential tool for food scientists and researchers. However, advanced knowledge of scientific programming and mathematical modelling are often required in order to use them, especially in cases of modelling of dynamic and/or non-linear processes. This may be an obstacle for food scientists without such skills. Scientific software can help making these tools more accessible for scientists lacking a deep mathematical or computing background. Recently, the R package bioinactivation was published, including functions (model fitting and predictions) for modelling microbial inactivation under isothermal or dynamic conditions. It was uploaded to the Comprehensive R Archive Network (CRAN), but users need basic R programming knowledge in order to use it. Therefore, it was accompanied by Bioinactivation SE, a user-friendly web application including selected functions in the software for users without a programming background. In this work, a new web application, Bioinactivation FE, is presented. It is an extension of Bioinactivation SE which includes an interface to every function in the bioinactivation package: model fitting of isothermal and non-isothermal experiments, and generation of survivor curves and prediction intervals. Moreover, it includes several improvements in the user interface based on the users' feedback. The capabilities of the software are demonstrated through two case studies using data published in the scientific literature. In the first case study, the response of Escherichia coli to isothermal and non-isothermal treatments is compared, illustrating the presence of an induced thermal resistance. In the second, the effect of nanoemulsified d-limonene on the thermal resistance of Salmonella Senftenberg is quantified.

Mathematical quantification of the induced stress resistance of microbial populations during non-isothermal stresses

This contribution presents a mathematical model to describe non-isothermal microbial inactivation processes taking into account the acclimation of the microbial cell to thermal stress. The model extends the log-linear inactivation model including a variable and model parameters quantifying the induced thermal resistance. The model has been tested on cells of Escherichia coli against two families of non-isothermal profiles with different constant heating rates. One of the families was composed of monophasic profiles, consisting of a non-isothermal heating stage from 35 to 70°C; the other family was composed of biphasic profiles, consisting of a non-isothermal heating stage followed by a holding period at constant temperature of 57.5°C. Lower heating rates resulted in a higher thermal resistance of the bacterial population. This was reflected in a higher D-value. The parameter estimation was performed in two steps. Firstly, the D and z-values were estimated from the isothermal experiments. Next, the parameters describing the acclimation were estimated using one of the biphasic profiles. This set of parameters was able to describe the remaining experimental data. Finally, a methodology for the construction of diagrams illustrating the magnitude of the induced thermal resistance is presented. The methodology has been illustrated by building it for a biphasic temperature profile with a linear heating phase and a holding phase. This diagram provides a visualization of how the shape of the temperature profile (heating rate and holding temperature) affects the acclimation of the cell to the thermal stress. This diagram can be used for the design of inactivation treatments by industry taking into account the acclimation of the cell to the thermal stress.

A Novel Technique for Sterilization Using a Power Self-Regulated Single-Mode Microwave Cavity

In this paper, a novel technique to achieve precise temperatures in food sterilization has been proposed. An accurate temperature profile is needed in order to reach a commitment between the total removal of pathogens inside the product and the preservation of nutritional and organoleptic characteristics. The minimal variation of the target temperature in the sample by means of a monitoring and control software platform, allowing temperature stabilization over 100 °C, is the main goal of this work. A cylindrical microwave oven, under pressure conditions and continuous control of the microwave supply power as function of the final temperature inside the sample, has been designed and developed with conditions of single-mode resonance. The uniform heating in the product is achieved by means of sample movement and the self-regulated power control using the measured temperature. Finally, for testing the sterilization of food with this technology, specific biological validation based on Bacillus cereus as a biosensor of heat inactivation has been incorporated as a distribution along the sample in the experimental process to measure the colony-forming units (CFUs) for different food samples (laboratory medium, soup, or fish-based animal by-products). The obtained results allow the validation of this new technology for food sterilization with precise control of the microwave system to ensure the uniform elimination of pathogens using high temperatures.

Nanoemulsified D-Limonene Reduces the Heat Resistance of Salmonella Senftenberg over 50 Times

Salmonella Senftenberg is a pathogen agent causative of foodborne disease and it is considered the most heat-resistant serovar within this genus. Food industries use heat treatment and chemical antimicrobials in order to eliminate this microorganism in food, but consumers prefer natural antimicrobials as essential oils and their components. This study evaluates the combined effect of thermal treatments and different concentrations of D-limonene nanoemulsion on the inactivation of Salmonella (S.) Senftenberg. The results showed an important effect of the nanoemulsified D-limonene on the heat resistance of S. Senftenberg. The δ_{50 °C} value was reduced by 85%, 96% and 98% when 0.1, 0.5 and 1 mM of nanoemulsified D-limonene was added to the heating medium. The effect was kept along all the heating temperatures researched and the shape of the survival curves did not change with the addition of the antimicrobial. The results obtained in this research could be very useful for food industries for optimizing or improving heat treatments applied to food.