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.
Keywords: gastric digestion process; multilevel analysis; quantitative microbial risk assessment; uncertainty; variability.