Co-optimization of safety, quality and legislation: opening Pandora’s box?

The probability of bacterial spores surviving a thermal process: The 12D myth and other issues with its quantitative assessment

The concepts of "D-value," "thermal death time" and "commercial sterility," innovative and useful at their inception, are based on untenable assumptions, notably that the log-linear isothermal inactivation model has universal applicability, that extrapolation over several orders of magnitude below the detection level is permissible, and that total microbial inactivation is theoretically impossible. Almost all commonly observed inactivation patterns, the log-linear is just a special case, can be described by both deterministic and fully stochastic models, examples of which are given. Unlike the deterministic, the stochastic models predict either complete elimination of the targeted cells or spores in realistic finite time, or residual survival. In most cases, the published survival data do not contain enough information to establish which actually happens. The microbial safety of thermally processed foods can be compromised not only by under-processing but also by a variety of mishaps whose occurrence probabilities are unrelated to the inactivation kinetics. Moreover, the available sampling plans to detect microbial contamination in sterilized containers through incubation alone are insensitive to levels of potential safety concerns. In principle, many of these issues could be resolved by developing new dramatically improved detection methods and/or verifiable methods to predict very low levels of microbial survival.

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.

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.