Prediction of Bacillus weihenstephanensis acid resistance: The use of gene expression patterns to select potential biomarkers

Exposure to mild stress conditions can activate stress adaptation mechanisms and provide cross-resistance towards otherwise lethal stresses. In this study, an approach was followed to select molecular biomarkers (quantitative gene expressions) to predict induced acid resistance after exposure to various mild stresses, i.e. exposure to sublethal concentrations of salt, acid and hydrogen peroxide during 5 min to 60 min. Gene expression patterns of unstressed and mildly stressed cells of Bacillus weihenstephanensis were correlated to their acid resistance (3D value) which was estimated after exposure to lethal acid conditions. Among the twenty-nine candidate biomarkers, 12 genes showed expression patterns that were correlated either linearly or non-linearly to acid resistance, while for the 17 other genes the correlation remains to be determined. The selected genes represented two types of biomarkers, (i) four direct biomarker genes (lexA, spxA, narL, bkdR) for which expression patterns upon mild stress treatment were linearly correlated to induced acid resistance; and (ii) nine long-acting biomarker genes (spxA, BcerKBAB4_0325, katA, trxB, codY, lacI, BcerKBAB4_1716, BcerKBAB4_2108, relA) which were transiently up-regulated during mild stress exposure and correlated to increased acid resistance over time. Our results highlight that mild stress induced transcripts can be linearly or non-linearly correlated to induced acid resistance and both approaches can be used to find relevant biomarkers. This quantitative and systematic approach opens avenues to select cellular biomarkers that could be incremented in mathematical models to predict microbial behaviour. Keywords: Bacillus; Acid resistance; Gene expression; Modelling; Biomarker; Food safety