From transcriptional landscapes to the identification of biomarkers for robustness
© Abee et al; licensee BioMed Central Ltd. 2011
Published: 30 August 2011
Skip to main content
© Abee et al; licensee BioMed Central Ltd. 2011
Published: 30 August 2011
The ability of microorganisms to adapt to changing environments and gain cell robustness, challenges the prediction of their history-dependent behaviour. Using our model organism Bacillus cereus, a notorious Gram-positive food spoilage and pathogenic spore-forming bacterium, a strategy will be described that allows for identification of biomarkers for robustness. First an overview will be presented of its two-component systems that generally include a transmembrane sensor histidine kinase and its cognate response regulator, allowing rapid and robust responses to fluctuations in the environment. The role of the multisensor hybrid kinase RsbK and the PP2C-type phosphatase RsbY system in activation of the general stress sigma factor σB is highlighted. An extensive comparative analysis of transcriptional landscapes derived from B. cereus exposed to mild stress conditions such as heat, acid, salt and oxidative stress, revealed that, amongst others σB regulated genes were induced in most conditions tested. The information derived from the transcriptome data was subsequently implemented in a framework for identifying and selecting cellular biomarkers at their mRNA, protein and/or activity level, for mild stressinduced microbial robustness towards lethal stresses. Exposure of unstressed and mild stress-adapted cells to subsequent lethal stress conditions (heat, acid and oxidative stress) allowed for quantification of the robustness advantage provided by mild stress pretreatment using the plate-count method. The induction levels of the selected candidate-biomarkers, σB protein, catalase activity and transcripts of certain proteases upon mild stress treatment, were significantly correlated to mild stress-induced enhanced robustness towards lethal thermal, oxidative and acid stresses, and were therefore suitable to predict these adaptive traits. Cellular biomarkers that are quantitatively correlated to adaptive behavior will facilitate our ability to predict the impact of adaptive behavior on cell robustness and will allow to control and/or exploit these adaptive traits. Extrapolation to other species and genera is discussed such as avenues towards mechanism-based design of microbial fitness and robustness.
Bacillus cereus is a member of the genus Bacillus, which is a highly heterogeneous group of bacteria that includes species with large variations in phenotypes, nutritional requirements and other physiological and metabolic characteristics. Within the genus Bacillus, B. cereus and its closest relatives form a highly homogeneous subdivision, which has been termed the “B. cereus group”. This group comprises the species B. cereus, B. anthracis, B. thuringiensis, B. mycoides, B. pseudomycoides and B. weihenstephanensis. B. anthracis is the etiological agent of the acute and fatal disease anthrax in mammals. It is the most notorious member of the B. cereus group, especially since it was used in bioterrorist attacks in the USA in 2001 .
Members of the B. cereus group are ubiquitously present in soil and can adapt to a wide range of environmental conditions [2–4]. The ecological niches and life cycles of B. anthracis and B. thuringiensis may be more specialized than those of B. cereus. B. anthracis spores are ingested by herbivores and germinate within the host to produce vegetative cells, which multiply and produce virulence factors, ultimately killing the host. Upon death, large numbers of B. anthracis cells are released into the soil, where they sporulate upon contact with air, thus completing the B. anthracis life cycle . The life cycle of B. thuringiensis may be comparable to that of B. anthracis. Upon death of its insect host, B. thuringiensis is released into the soil where it is a ubiquitous inhabitant and where it can germinate and grow under favorable conditions. . The niches and life cycle of pathogenic B. cereus remain more obscure. B. cereus spores can be found in many types of soils, sediments, dust and plants [3, 4, 6]. For this reason, B. cereus is considered to be mainly a soil inhabitant and, indeed, B. cereus can germinate, grow and sporulate in soil, thus demonstrating a saprophytic life cycle . However, genome analyses have indicated that B. cereus specializes in protein metabolism, pointing towards a symbiotic or parasitic life cycle . This is further supported by the fact that B. cereus has been isolated from guts of soil-dwelling arthropods  and from stool samples of healthy humans . Furthermore, B. cereus has been isolated from surface waters, a niche that may allow this bacterium to easily enter food-processing lines . Many types of food have been associated with B. cereus food-borne disease, including meats, vegetables, puddings, milk, rice, pasta and noodles [3, 11].
B. cereus can cause two distinct types of food-borne disease: the diarrhoeal and the emetic type. Although both types are generally mild and self-containing, more serious and even lethal cases have been reported [12–15]. The diarrhoeal disease is often associated with protein-rich foods, such as meat, vegetables, puddings and milk products. The diarrhoeal disease is thought to be caused by vegetative cells (ingested as viable cells or spores) that produce enterotoxins in the small intestine. Typical symptoms include abdominal pains, watery diarrhoea, nausea and vomiting. The incubation time generally ranges between 8-16 hours after ingestion, while the symptoms normally last for 12-24 hours. However, longer incubation times have been observed and the symptoms can last for up to several days . The emetic disease is often associated with starch-rich foods, such as fried and cooked rice, pasta and noodles. It was first identified in 1974 when B. cereus was linked to several outbreaks caused by eating cooked rice in the United Kingdom in the early 1970s . The emetic disease is caused by the B. cereus emetic toxin, cereulide, which is produced in foods before ingestion. Symptoms mainly include nausea and vomiting, which occur between 30 minutes to 6 hours after ingestion and which generally last for 6-24 hours .
Because of the medical and economic relevance of B. cereus group members, a large number of genomes has been sequenced. Starting with the first publications in 2003 [8, 18], the list of sequenced members now includes five B. anthracis, eleven B. cereus, three B. thuringiensis, one B. weihenstephanensis, and an additional 89 sequences are in progress or available with high coverage draft sequence (Dec 2010; http://www.ncbi.nlm.nih.gov/genomes/Iproks.cgi). These sequences provide an extensive source for the identification of signaling systems and key factors in adaptive responses towards environmental and food-relevant stresses.
The ubiquitous presence and high adaptive capabilities of B. cereus and its closest relatives may be explained by the fact that these species harbor dedicated signaling pathways that allow rapid and robust responses to fluctuations in their environment. In bacteria, including archaebacteria, the dominant signaling pathways for monitoring environmental cues are the so-called two-component systems (TCS), which include a transmembrane sensor histidine kinase (HK) and its cognate response regulator (RR). The mode of signal transduction by TCSs involves a phosphotransfer reaction between a conserved histidine and aspartate residue located in the HK phosphotransferase and RR receiver domain, respectively. RRs generally function as transcription factors that, upon phosphorylation, bind to specific sites on the DNA to alter the expression of genes involved in adaptive responses . These systems are almost ubiquitously present in bacteria  and analyses of whole-genome and metagenome sequences have revealed that TCS protein domains are the second most numerous Pfam domains in bacteria, exceeded only by ABC-type transporter domains [21, 22]. TCSs are known to monitor a wide variety of conditions, such as nutrient deprivation, cold/heat shock and the presence of antimicrobial compounds [23–25], and they can control virtually all types of microbial behaviour, including motility and chemotaxis, sporulation, biofilm formation and quorum sensing [26–29].
Recent studies have identified the complex, multisensor hybrid HK RsbK and the PP2C-type phosphatase-containing RR RsbY (Fig. 1), to be involved in controlling σB activity [32, 33]. In several Gram-positive bacteria, including the model Gram-positive Bacillus subtilis, the alternative sigma factor σB plays a prominent role in redirecting gene expression under stress conditions . In B. cereus, activation of σB is also induced under a variety of stress conditions, most notably heat stress, and B. cereus sigB and rsbK deletion strains displayed an impaired heat adaptive response [32, 35]. In B. anthracis, σB also plays a role in heat stress and a B. anthracis sigB deletion strain was slightly affected in its virulence .
Further experiments in combination with in silico analysis [32, 37] revealed that σB activation in B. cereus significantly differs from known σB activation pathways in other low-GC Gram-positives. In B. subtilis, but most likely also in many of its closest relatives, a massive protein complex called the stressosome integrates multiple environmental stress signals to orchestrate dephosphorylation of RsbV and the eventual release and activation of σB[38, 39]. In addition, B. subtilis even uses a second pathway, which monitors energy stress, to regulate σB activity . Both of the above pathways are absent in members of the B. cereus group, where RsbY is the sole PP2C-type phosphatase for the control of σB activity . RsbY most likely receives its input signals from the hybrid HK RsbK, either indirectly via an extended phosphorelay or directly. The domain architecture of the unique sensory protein RsbK suggests that in members of the B. cereus group, environmental and intracellular stress signaling routes are combined into one single protein. This activation mechanism may reflect unique niches in which B. cereus group members can reside .
The sodium hypochlorite landscape consists of 37 genes in which there is a clear overrepresentation of genes predicted to have reductase activity (BC1835, BC1844, BC2194, BC3024,) and arsenate resistance (BC3152-BC3154). In addition several transport systems involved in the transport of metal ions (copper, cadmium, zinc) are activated (BC3731, BC3732, BC0595, BC0596). Typically, we also observed the activation of the transcriptional regulator spxA (BC3402) and three MarR-type transcriptional regulators (BC3025, BC4474, BC5038,). These observations extend those of Ceragioli et al. , who performed a comparative analysis with cellular responses to three other disinfectants tested, i.e., peracetic acid, benzalkoniumchloride and hydrogen peroxide, and concluded that exposure to sodium hypochlorite induced genes involved in metabolism of sulfur and sulfur-containing amino acids. These results correlated with the excessive oxidation of sulfhydryl groups observed in sodium hypochlorite-stressed cells. For the genes found in the acid cluster it is less clear what their role could be in the response to the stress; apart from three genes coding for proteins with a putative function (BC2087, BC2858, BC4640), some single genes were found to play role in quite diverse processes such as hydrolysis (BC3414), phospholipid metabolism (lipase; BC4333), iron uptake (BC5381), drug resistance (BC0860) and peptide breakdown (BC0789). Notably, expression activation of genes with putative functions in lipid turnover, iron metabolism and drug resistance have previously been reported for sorbic acid-stressed cells of B. subtilis.
In contrast to the clusters that show clear response to only a single stress, other clusters were found that displayed a more diverse response. One of these clusters contained a large fraction of the σB regulon, including sigB, rsbW, an anti-sigma factor antagonist (BC1002), BC099, and katE, and several involved in stress response and DNA repair (BC0062, BC0103). Another cluster included four genes of the CtsR-regulon (BC0099-0102); genes known to code for proteins involved in the repair of protein damage . Other stress response genes were present in clusters of genes active under many different stress conditions, e.g., hrcA, spxA, groES/EL, grpE, clpX and BC0377 (putative alkyl hydroperoxide reductase) (data not shown). This pointed to a significant role of these components in general stress adaptation and, as expected, they have known canonical roles in stress adaptation including regulation of stress responses, removing reactive oxygen species, and protein repair and protein quality maintenance [2, 51, 52]. Although the transcriptome response of the genes in the clusters seems to be very broad, the genes do not respond to all stresses. Especially for the sigB cluster, there is a clear distinction between stress conditions under which gene expression is present (benzalkonium chloride, peracetic acid, sodium hypochlorite, heat and salt) and conditions where gene expression was not observed at a high level (cold and hydrogen peroxide). Upon acid stress, the response is even variable between different conditions (depending on the type of acid and the time of exposure). It will be of interest now to find the common denominator among the stresses that trigger σB as this may point to the exact stimulus/stimuli to which RsbK responds.
In order to evaluate whether catalase activity could also predict the robustness enhancement towards multiple lethal stresses upon pretreatment to mild oxidative stress, the conditions where catalase activity functioned as biomarker, namely lethal oxidative stress and lethal heat stress were combined. The obtained correlation was tested for significance. Because the correlation remained significant also after combining these responses, the results highlight that catalase activity could also act as biomarker for mild stress induced robustness towards multiple lethal stresses . Mild oxidative stress treatment provided cross-protection towards multiple lethal stresses, e.g. lethal oxidative, heat and acid stress. The induction of catalase activity and activation of other stress-related responses upon mild oxidative treatment may have contributed to these multiple cross-protective effects. Recently, it was demonstrated that aerobic lethal heat and acid stress exposure imposed an oxidative stress burden in B. cereus cells [55, 56], and likely, mild oxidative stress pre-treatment provides a survival advantage for the cells.
The predictive potential of the candidate-biomarkers was highly influenced by the functional cell levels at which the candidate-biomarkers were measured. The katA transcript of the main vegetative catalase and the catE transcript were both not suitable to predict the robustness status of mild oxidative stress adapted cells towards the three lethal stresses tested . This underlined the importance to evaluate the predictive potential of cellular indicators and different functional cell levels. The better predictive potential of enzyme activity compared to transcript levels might reflect the more transient nature of gene expression than of enzyme activity level. None of the other nine selected candidate-biomarkers could predict the robustness advantage of mild oxidative stress-treated cells, but they emerged as biomarkers for cells that were adapted to mild heat, acid and/or salt. The SigB protein proved to be a suitable biomarker for mild stress-induced robustness towards lethal heat stress, whereas the proteases transcripts clpC and clpP could predict the mild stress-induced robustness advantage towards lethal acid stress .
Catalases have known crucial roles in adaptive stress responses and are widely conserved among bacteria and eukarya [52, 57]. The significant role of oxidative stress defense mechanisms in stress adaptation suggests that the predictive potential of oxidative stress related candidate-biomarkers may extend beyond the species B. cereus. Moreover, the other selected candidate-biomarkers that have known functions in regulation of stress adaptation and maintenance of protein quality proved to be suitable to predict the robustness advantage of B. cereus for various stress adaptive traits . Regulators of general stress responses and cellular systems involved in oxidative stress defense and maintenance of protein quality are among the most consistently induced components in microbial stress responses. Induction of these stress responses upon stress treatment has been demonstrated for a wide range of other micro-organisms, including yeast , spoilage organisms  and, recently, also for bifidobacteria  which are applied in functional foods. Induction of robustness must be averted for food-borne pathogens, but can give organisms used in functional foods a survival advantage. The canonical role of the evaluated stress response systems may point to possible predictive potential also for micro-organisms other than B. cereus. Variations in conservation of these components and their regulatory networks underline the need for profound validation of the predictive quality of promising candidate-biomarkers because that might be species or even strain specific. This is clearly illustrated by the regulation of katA in B. cereus as compared to that in other species. B. cereus katA seems to be negatively controlled by the regulator PerR, just as is the case in L. monocytogenes, S. aureus, and B. subtilis, . However, katA of B. cereus also seems to be (weakly) positively controlled by σB under stress conditions : a regulatory module/connection that has also been found in S. aureus[39, 61], but does not seem to occur in e.g. B. subtilis and L. monocytogenes. In turn, the activation mechanisms of σB in B. cereus differs from that in S. aureus, B. subtilis, and L. monocytogenes i.e., σB is activated via RsbK-RsbY in B. cereus (Fig. 1) [32, 37], via RsbU (and putative upstream components) in S. aureus, and via RsbRSTU in B. subtilis and L. monocytogenes. In addition, B. subtilis uses the RsbPQ pathway to control σB activity .
Our systematic and quantitative approach to search for promising candidate-biomarkers and to evaluate their predictive potential provides perspectives to search for biomarkers in other micro-organisms. The search for biomarkers for stress adaptive behavior will contribute to a better understanding of stress adaptation mechanisms and can serve to detect and predict stress adaptive traits at an early stage.
In this review, we have integrated three different research strategies that have led to the identification of biomarkers important for robustness in Bacillus cereus. Firstly, a general overview of the TCS arsenal of B. cereus was obtained using comparative genomics. From this global analysis, a potential key TCS protein (RsbK) could be linked to the earlier defined RsbY-SigB stress-responsive network and further molecular work indeed confirmed the role of RsbK in controlling SigB activity and its regulon. Secondly, a meta-transcriptome data analysis again showed, amongst others, the importance of the RsbKY-controlled SigB regulon in the general B. cereus stress response. Thirdly, it was shown that SigB protein levels, protease transcript levels and catalase activity could predict the robustness level of stress-adapted B. cereus cells towards lethal stresses.
The drive to use more mechanism-based approaches for evaluating robustness of micro-organisms urged the search for biomarkers to early detect and predict stress adaptive behavior. Prediction of phenotypic behavior using cellular indicators is a key area of research in the field of food-borne pathogens, spoilage organisms, and organisms used in functional food applications [60, 63–65]. Quantitatively correlating microbial responses at molecular and phenotypic levels can provide mechanistic understanding of stress adaptive behavior and leads for identifying cellular indicators for bacterial performance. Next to the general stress response regulator sigmaB, cellular components such as catalase activity and proteases were demonstrated to be suitable to predict the robustness level of stress-adapted B. cereus cells towards lethal stresses. The predictive quality of the biomarkers was influenced by the functional cell level at which the biomarker was measured, i.e., transcript, protein or enzyme activity level. The level at which the candidate-biomarkers could be employed as biomarker was demonstrated to be stress-dependent. Not only protein levels and enzyme activity but also transcript levels showed high predictive quality. However, the transient nature of expression of transcripts might point to a better predictive potential for the corresponding biomarker if it is assessed at protein and/or enzyme activity level as shown for sigmaB levels and catalase activity. This underlined that evaluation of predictive potential at different functional cell levels is crucial to select robust biomarkers.
Obviously, the different research strategies that were employed may reveal or have already revealed other potential biomarkers for B. cereus robustness behaviour, besides the ones mentioned above. For instance, when considering the large arsenal of B. cereus TCS proteins, it is noteworthy to mention that of the 51 different HK-RR pairs found in the B. cereus group, around 50% was found to be homologous to HK-RR pairs with a known function in at least some form of (general) stress responsive behavior in other species [30, 31]. In addition, the biological role for around 20% of the B. cereus HK-RR pairs remains completely unknown and it is plausible that at least some of these modules are involved in general stress responsive behavior of B. cereus. Finally, 14 different “orphan” HKs are likely to function in sporulation initiation in the B. cereus group (Fig. 1) and these may turn out to provide excellent biomarkers for another type of B. cereus robustness behavior, namely the formation of dormant endospores.
Further investigation along these lines by varying stress intensities and testing other micro-organisms will reveal how vigorous these quantified correlations are and is required to understand the limits of reliability of the predictive potential of biomarkers. Recently, the response of Bifidobacterium breve, a Gram-positive bacterium isolated from the nursling stool of a breast-fed infant, to several stresses (heat, osmotic, solvent, oxidative) was assessed using transcriptomics. Integration of the results with in silico analysis, allowed the formulation of a model for an interacting regulatory network for stress response in this bifidobacterium , involving HspR, ClgR and HrcA. The latter regulator, and other stress-induced genes in all tested conditions belonging to oxidative stress response and protective proteins such as molecular chaperones and proteases were also generally up-regulated in stressed B. cereus suggesting that at least partial overlap in putative biomarkers for robustness may exist. This is further supported by transcriptomic and phenotypic analysis of stress adaptive responses in a range of gram-positive bacteria including lactic acid bacteria [66, 67], Bacillus spp. [68, 69], and food-borne human pathogens such as L. monocytogenes[70, 71], that revealed activation of (homologues of ) the proposed canonical biomarkers under a variety of stress conditions.
Studies could be extended with assessment of biomarkers for host-microbe interactions including virulence . Life cycle transcriptomics and comparative genomics have been used recently to unravel successive and coordinated global transcriptional changes during infection of amongst others L. monocytogenes, and pointed to previously unknown mechanisms in bacteria with a crucial role for σB-mediated activation of virulence genes in the host intestinal lumen .
Robust biomarkers will be of indisputable significance to better understand the mechanisms of stress adaptation, and will complement an empirical approach to evaluate stress robustness of micro-organisms. Insights obtained may on the one hand aid in the design and development of more efficient (novel combination) treatments for control of spoilage organisms and food-borne pathogens and on the other hand provide tools to enhance robustness of organisms used in functional food applications.
The authors gratefully the acknowledge the fruitful and efficient cooperation with other members of the TI Food and Nutrition project on Food Preservation and Safety, and the support of colleagues at the different affiliations.
This article has been published as part of Microbial Cell Factories Volume 10 Supplement 1, 2011: Proceedings of the 10th Symposium on Lactic Acid Bacterium. The full contents of the supplement are available online at http://www.microbialcellfactories.com/supplements/10/S1.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.