INTRODUCTION
The
Bacillus cereus group, also known as
B. cereus sensu lato, consists of ubiquitous Gram-positive spore-forming bacteria. This group comprises seven closely related species, including
B. cereus sensu stricto and
B. cytotoxicus, that cause emetic and/or diarrheal syndromes in food poisoning (
1–3). Some rare strains of
B. weihenstephanensis, another species of this bacterial group, can also be emetic (
4). These bacteria cause a number of food safety issues, as they are able to produce heat-resistant spores and thus survive in food even after processes such as cooking or pasteurization. Compared to the incidence of food poisoning caused by
B. cereus sensu stricto, food poisoning caused by
B. cytotoxicus or
B. weihenstephanensis is rare (
1,
3). The risk for consumers is mostly tied to high doses (10
5 to 10
8 CFU) of cells or spores contaminating the ingested food and thus depends on the ability of the bacterium to multiply during the food's shelf life (
2,
5,
6). As some strains of
B. cereus have the ability to grow at temperatures found in the chill chain, identifying the mechanisms involved in low-temperature adaptation may help predict how this bacterium behaves in refrigerated food and enable more accurate risk prediction and better risk prevention.
Bacterial adaptation to low temperature is a complex and multifactorial process involving both the genetic background of the bacteria (
7) and an array of mechanisms (
5). Concerning the genetic background, seven phylogenetic groups (groups I to VII) were defined in
B. cereus sensu lato, and
B. cereus sensu stricto was broadly positioned in these groups (groups II to VI), while
B. weihenstephanensis and
B. cytotoxicus were merged with groups VI and VII, respectively (
7).
B. cereus sensu stricto is thus expected to show a broadly diverse genetic background, in contrast to
B. cytotoxicus or
B. weihenstephanensis. Interestingly, each phylogenetic group (groups I to VII) was assigned a specific range of growth temperatures on the basis of both genetic and phenotypic criteria (
7). These seven phylogenetic groups could also be seen as seven thermotypes, running from a psychrotolerant (cold-tolerant) group to moderately psychrotolerant, mesophilic, and moderately thermotolerant (heat-tolerant) groups.
B. cereus sensu stricto thus spans different thermotypes containing psychrotolerant, mesophilic, or intermediate strains, whereas
B. weihenstephanensis contains only psychrotolerant strains and
B. cytotoxicus contains only moderately thermotolerant strains. It was suggested that in the course of evolution, changes in temperature tolerance limits have fashioned historical patterns of global ecological diversification in
B. cereus sensu lato (
7).
The mechanisms involved in adaptation to low temperature are equally multifactorial, as illustrated by the various mechanisms identified so far in other model bacteria and also described in
B. cereus (
5). For instance, membrane fatty acid (FA) composition adjustments can increase the proportion of low-melting-point FAs (like unsaturated FAs and branched-chain FAs) (
8–10). RNA helicases that enable the RNA unfolding needed for proper translation and/or RNA degradation also play a major role in
B. cereus low-temperature adaptation (
11).
Among the mechanisms that allow organisms to change in response to environmental conditions, two-component systems (2CSs) are signal transduction systems that are almost ubiquitous in bacteria (
12). 2CSs are known to sense a wide range of environmental stressors, enabling cells to elaborate a response by regulating the expression of genes required for adaptation (
13–17). 2CSs basically comprise a histidine kinase (HK) that senses an environmental stimulus (either directly or after interaction with accessory proteins) and a cognate response regulator (RR) that usually functions as a transcriptional regulator. One of the first studies on a 2CS involved in cold adaptation identified DesKR, which is able to sense an increase in membrane thickness in response to a decrease in temperature in
B. subtilis (
18,
19). DesKR consequently regulates the expression of a desaturase gene responsible for maintaining membrane fluidity during
B. subtilis growth at low temperature. More recently, a
Clostridium botulinum 2CS important for cold tolerance was discovered, and the mechanisms involved were investigated (
20,
21). In a Gram-negative bacterium, the CheA/CheY 2CS mutant of
Yersinia pseudotuberculosis was impaired during growth at 3°C (
22). 2CSs of other Gram-negative bacteria have been described to be temperature sensors for bacterial virulence control, such as CorSR in
Pseudomonas (
23) and PhoPQ in
Edwardsiella (
24). Recent research has revealed how a 2CS contributes to the high adaptability of
B. cereus strains that enables these bacteria to persist in processed foods (
25). Some 2CSs were shown to play diverse roles in
B. cereus adaptation (
26–30), but most of the 2CSs found among
B. cereus sensu lato strains have an unknown function (
31), and none has yet been shown to be involved in low-temperature adaptation.
In the course of experiments to better understand the role of the numerous 2CSs with unknown function found among B. cereus sensu lato strains, we mutated one of them (BC_2216-BC_2217) in the mesophilic model strain B. cereus sensu stricto ATCC 14579. We found that this new 2CS is necessary for low-temperature adaptation not only in mesophilic strains but also in two psychrotolerant strains belonging to distinct phylogenetic groups of B. cereus sensu lato.
RESULTS
The BC_2216 and BC_2217 genes encode a 2CS and are overexpressed at low temperature.
According to their annotation, the BC_2216 and BC_2217 genes found in the B. cereus ATCC 14579 mesophilic strain encode, respectively, a putative histidine kinase and a response regulator. The predicted BC_2216 protein indeed contains the classical transmitter domains of histidine kinases: a ATP binding domain (HATPase domain, residues 288 to 372) and a dimerization and phosphoacceptor domain (Pfam:HisKA_3 domain, residues 186 to 247). BC_2216 has no predicted transmembrane domain, but its N-terminal region displays a GAF domain (Pfam 01590). BC_2217, the presumed cognate response regulator of BC_2216, displays a phosphoacceptor site (REC domain, residues 6 to 120) and a DNA-binding motif (helix-turn-helix domain, residues 151 to 207), suggesting that BC_2217 has DNA regulatory ability. The BC_2216 and BC_2217 proteins thus presumably form a 2CS.
The level of expression of the BC_2216 gene was quantified by RT-quantitative PCR during low-temperature (12°C) growth and compared to that at an optimal temperature (37°C). RNA samples were extracted from cells collected at three times during growth kinetics, in mid-exponential phase (OD600 = 0.5), end exponential phase (OD600 = 1.0), or stationary phase (OD600 = 2.5). Results showed 7.2-fold, 10.4-fold, and 9.7-fold overexpression, respectively, during growth at 12°C compared to that during growth at 37°C, suggesting an important role of the BC_2216-BC_2217 2CS under this growth condition.
Mutation of the BC_2216 and BC_2217 genes causes impaired low-temperature growth.
A BC_2216-BC_2217 mutant of the
B. cereus ATCC 14579 strain was constructed by allelic exchange between the two 2CS-encoding genes and a kanamycin resistance cassette (
Table 1), and the phenotype of this mutant was compared to that of its parental strain during growth at various temperatures (
Fig. 1).
At the Tmin of this strain (i.e., 10°C), the wt growth curves varied slightly between experiments. The lag phase of the 3 replicates was about 114 ± 14 h (mean ± standard error of the mean [SEM]), and they reached stationary phase in 196 ± 21 h after inoculation, with a maximal OD of 2.3 ± 0.1. Under the same conditions, the BC_2216-BC_2217 mutant did not show any growth, even after 300 h of incubation.
At 12°C, the mutant was able to grow, but the mean ± SEM lag phases of the wt strains (91.3 ± 6.9 h) and the mutant strains (189.3 ± 6.7 h) from three independent cultures were significantly different (P < 0.05, Student t test), indicating that the mutant was slower to adapt to this cold temperature than the wt strain. The maximum OD was also significantly lower for the mutant than for the wt (3.68 ± 0.10 versus 4.58 ± 0.19, respectively; P < 0.05). Thus, mutant growth was impaired at 12°C.
In contrast, at the optimal temperature (i.e., 37°C), the growth curves of the two strains were similar. The cells reached the stationary phase, with only a slight delay of 1.19 h (P < 0.05) for the mutant compared to the time to reach the stationary phase for the wt strain, but maximal ODs were not significantly different (2.01 ± 0.09 for the mutant versus 2.06 ± 0.01 for the wt; P > 0.05).
Microscopic observations were performed during growth at low and optimal temperatures (
Fig. 2). Both the wt and mutant strains showed slight cell elongation during growth at 12°C compared to the cell size with growth at 37°C (
Fig. 2A). There was no difference in cell size between the wt and mutant strains. Cells of the two strains checked at a higher magnification by TEM showed a similar cell structure at the optimal growth temperature (i.e., 37°C) (
Fig. 2B). In contrast, when incubated at low temperature (12°C), cells of the mutant strain displayed an atypical morphology compared to that of the wt strain, as their cytoplasm displayed a wavy shape, suggesting an irregular thickness of the cell wall.
Given that cold is a major stressor for bacteria, we investigated whether the growth of the BC_2216-BC_2217 mutant was impaired when it encountered other stressful conditions. Growth of both strains was measured under various conditions: under conditions of high temperature, low pH, high pH, and high osmolarity and in the presence of ethanol, ion chelators, or oxidative agents. Under all the tested conditions, the BC_2216-BC_2217 mutant strain showed growth similar to that of the wt (see Table S1 in the supplemental material).
In addition, the phenotype of the mutant was compared to that of its parental strain after growth on various media. Growth on Mossel agar and growth on sheep blood agar showed that the mutant and the wt strains displayed the same lecithinase activity and hemolytic profile, respectively. Tests run with an API 50CH kit to determine the bacterial oxidation of 49 different carbohydrates did not find any differences between the two strains.
Taken together, these data suggest that the major phenotype of this mutant compared to that of the wt was its growth impairment at low temperature. We therefore propose the name casK (where cas represents cold adaptation sensor) for the BC_2216 gene, which encodes a putative histidine kinase, and casR for the BC_2217 gene, which encodes a putative response regulator.
Survival at temperatures below the minimal temperature of growth.
Given the impaired growth at low temperature, we investigated the ability of the ATCC 14579 wt and Δ
casKR strains to survive at temperatures below the minimal temperature of growth (i.e., below 10°C). The two strains were incubated at 4°C in LB, and the CFU were enumerated over time. As shown in
Fig. 3, viable counts of both strains decreased regularly over time during incubation at 4°C, but the viability loss was significantly faster for the
casKR mutant than for the wt strain. A viability loss of 3 log CFU was reached in 103 ± 6 h (mean ± SEM) for the mutant strain, whereas a viability loss of 3 log CFU was reached in 292 ± 46 h for the wt (
P < 0.05, Student
t test). Thus, after 300 h of incubation, the viable counts were significantly higher for the wt (4.7 ± 0.6 log CFU, mean ± SEM) than for the
casKR mutant (2.5 ± 0.3 log CFU) (
P < 0.05).
In order to determine whether long-term storage in LB medium at 37°C could also impair survival in the same way, we ran a control experiment. The results indicated that viability loss was very limited, with a 1.1-log-unit loss for the ATCC 14579 wt versus a 1.3-log-unit loss for the casKR mutant (means of duplicate experiments) after 400 h of incubation at 37°C (see Fig. S2 in the supplemental material). This result indicates that the ability to survive long-term storage in LB involves CasKR only at low temperature.
Complementation of the mutant phenotype.
Complementation of the Δ
casKR strain was performed by introducing on a plasmid the
casKR genes and a large upstream region presumably containing the
casKR promoter (
Table 1). While no growth was observed for the Δ
casKR strain without complementation at the
Tmin (i.e., 10°C), growth ability was restored for the Δ
casKR strain complemented with pHT-casKR (
Fig. 4). Growth was detected on LB agar even at a 10
−7 dilution, similar to what was observed for the ATCC 14579 wt strain (
Fig. 4). This result confirms that the
casKR deletion is genuinely responsible for the cold growth impairment of the mutant strain.
The CasKR paralog is not required for B. cereus low-temperature growth.
The B. cereus ATCC 14579 BC_5412 and BC_5411 genes present high sequence similarity with the casK and casR genes located elsewhere on the chromosome (49.5% and 54.0% identity, respectively). Consequently, the BC_5412-BC_5411 2CS could be considered a paralog of CasKR. We therefore investigated whether BC_5412-BC_5411 could play a role similar to that of CasKR during B. cereus growth at low temperature. A BC_5412-BC_5411 mutant was constructed, and its growth was tested at the optimal temperature and at low temperature. There was no observable difference between the BC_5412-BC_5411 mutant and its parental strain in terms of the kinetics of growth at 37°C (data not shown) or 12°C (see Fig. S3 in the supplemental material). At Tmin (i.e., 10°C), this mutant was still able to grow on LB agar similarly to the wt (see Fig. S3 in the supplemental material).
Role of CasKR among other B. cereus sensu lato strains.
To determine whether the
casKR genes are present among other strains of the various phylogenetic groups of
B. cereus sensu lato, we performed a search for
casKR (BC_2216-BC_2217) orthologs with conserved synteny of gene neighborhoods (see Fig. S4 in the supplemental material). Sequence similarity and E values are reported for each ortholog in Table S5 in the supplemental material. Orthologs of
casKR were found in strains belonging to five of the seven known phylogenetic groups of
B. cereus sensu lato, i.e., groups II to VI (
Table 3; see Fig. S4 and Table S5 in the supplemental material). While they were absent in group VII (a moderately thermotolerant group),
casKR orthologs were present in mesophilic group I, but synteny was not conserved, particularly for downstream genes.
Outside B. cereus sensu lato, tBLASTN analysis revealed that the best similarity was found in B. megaterium (51% identity for casK and 62% for casR), but the synteny was not conserved. In other species, similarity was lower and the synteny was still not conserved. For instance, in B. subtilis strain 168, the best similarity was obtained with yhcYZ, a 2CS of unknown function (48% identity for casK and 47% for casR). These genes, found outside B. cereus sensu lato strains, may be considered putative casKR orthologs, but their function remains to be investigated to clarify this point.
Thus, the
casKR genes are widespread among
B. cereus sensu lato strains. Psychrotolerant and moderately psychrotolerant strains are clustered in three out of seven phylogenetic groups of
B. cereus sensu lato (
7), and they all displayed orthologs of
casKR genes. We thus investigated whether the role of these orthologs was similar to that observed in the mesophilic ATCC 14579 strain (group IV), despite a lower
Tmin for these strains. Orthologs of the
casKR genes were also found among strains belonging to phylogenetic group III, which are considered to have a higher
Tmin (12°C) than the ATCC 14579 strain (10°C). We therefore investigated the role played by
casKR in three strains representative of three additional phylogenetic groups by constructing a
casKR mutant for each strain.
For the mesophilic AH187 strain (phylogenetic group III), the psychrotolerant MM3 strain (group II), and the psychrotolerant Rock 3-28 strain (group V) studied, the
casKR mutants and their parental strains displayed a similar growth ability at 37°C (
Fig. 5). In contrast, all three
casKR mutants showed impaired growth at low temperature. In addition, the mutants of the two psychrotolerant strains did not show any growth at 8°C, a temperature close to their
Tmins (i.e., 8°C for the Rock 3-28 strain and 7°C for the MM3 strain). These results showed that the CasKR 2CS is required for growth at low temperature among mesophilic as well as psychrotolerant
B. cereus sensu lato strains, whatever the value of the
Tmin. In addition, CasKR seems to be particularly efficient when the bacteria are submitted to a temperature close to their
Tmin.
DISCUSSION
In this study, we identified a pair of genes that play a major role in
B. cereus growth and survival at low temperature. CasKR, the newly identified 2CS, is widespread among
B. cereus sensu lato strains. 2CSs have a wide array of functions, and a few 2CSs have already been shown to be involved in the low-temperature response in other bacteria (DesKR in
B. subtilis, CheAY in
Y. pseudotuberculosis, CBO0365-CBO0366 in
C. botulinum) (
18,
20–22). According to the proposed 2CS classification scheme (
45,
46), CasK belongs to the class II family of histidine kinases and CasR belongs to the NarL family of cognate response regulators, which puts CasKR in the same family as DesKR. However, despite their signal transduction function, these two 2CSs have completely different genetic organizations (for instance, the desaturase-encoding genes are located on different loci in
B. cereus, contrary to what is observed in
B. subtilis) and the two HKs display no similarity within their N-terminal sensory domains. Among the various 2CSs present in the ATCC 14579 model strain of
B. cereus, CasKR (the BC_2216-BC_2217 2CS) and its paralog, the BC_5412-BC_5411 2CS, seem to be atypical, as their kinases do not display a transmembrane domain, suggesting that they belong to those rare HKs that have a cytoplasmic location (
47). Despite the sequence similarity between these two paralogs, only CasKR is overexpressed at low temperature (S. Chamot and J. Brillard, unpublished results) and only CasKR is necessary for optimal growth at low temperature. According to SMART analysis, the sensory N-terminal region of CasK contains a GAF domain (
44). These domains are able to bind cyclic nucleotides and are present, for instance, in cyclic GMP-specific phosphodiesterases or in phytochromes (
48). To our knowledge, a link between such a domain and sensing of low temperatures has never been described. This raises the question of the signal perceived by CasK under cold conditions. A 2CS from
Staphylococcus aureus, YhcSR, presents some sequence similarity with CasKR. As YhcS lacks a transmembrane domain, this histidine kinase probably has an intracytoplasmic location, like CasK (
49). This 2CS was shown to be essential for cell viability, and it controls the expression of an ABC transporter that seems to play a role under high-osmolarity conditions (
50). The full viability of our CasKR mutant and its unaffected growth in the presence of high NaCl concentrations compared to the growth of the wt suggest that YhcSR and CasKR play different roles in cell physiology.
Whatever B. cereus sensu lato strains among the several phylogenetic groups that were tested, ΔcasKR mutants of those strains showed impaired growth at low temperature compared to the growth of their parental strains. In ATCC 14579, survival at 4°C was also impaired in the ΔcasKR mutant compared to that of the wt. In contrast, the wt and the mutant showed similar survival at a subfreezing temperature (−20°C) (S. Diomandé and J. Brillard, unpublished data). These results suggest a major role of CasKR when cells are physiologically active (10°C) or when metabolic activity is strongly reduced (at 4°C), but not when cells are frozen.
Another phenotype of the CasKR mutant is the modified cell shape with an irregular cell wall instead of the elongated rod shape regularly observed for the wt cells at low temperature. Strong modifications of the structures of
B. cereus cells have been described before, for instance, in a
cshA RNA helicase mutant with a substantially stronger phenotype (
11). In another study under conditions involving both low temperature and low redox potential, elongation of the cells coupled to a default in cell separation was observed (
35).
Model strains are good tools to identify new mechanisms, but they may significantly differ phenotypically from food-poisoning strains (
51), which may also be true for
B. cereus (
29,
52,
53).
B. cereus sensu lato has the singularity of displaying a wide panel of strains: some are pathogenic and mainly belong to mesophilic or thermotolerant groups III, IV, and VII, whereas others belong to psychrotolerant groups II and V, where fewer pathogenic strains have been so far described, and the most psychrotolerant strains from group VI have been only marginally associated with human pathogenicity (
3,
4). Therefore, the ability of some
B. cereus sensu lato strains both to grow at low temperature and to cause human infections or food poisoning makes it important to identify the mechanisms involved in cold growth of both mesophilic and psychrotolerant strains. The CasKR 2CS has been identified in various strains belonging to phylogenetic groups II, III, IV, V, and VI, and we proved its involvement in low-temperature growth in strains representative of phylogenetic groups II, III, IV, and V, two of which are classified as psychrotolerant and two of which are classified as mesophilic. It is tempting to speculate that CasKR could also play a similar role in psychrotolerant strains belonging to phylogenetic group VI (
B. weihenstephanensis), but this was not investigated here because several attempts to obtain a
casKR mutant in two different strains from this group failed. Such unsuccessful assays could reveal a major role of CasKR in these strains, but clarification of this point will have to be deciphered by further investigations. Excluding the AH187 emetic strain (
54), the strains chosen here were not isolated from a food-poisoning outbreak. However, all these strains have the pathogenic potential to cause food-borne illness, as their genomes contain genes involved in virulence (e.g.,
nhe,
plcR) (
2,
55). Our results illustrate that psychrotolerant
B. cereus sensu lato strains could require CasKR to achieve efficient growth during food storage prior to a possible food-poisoning event.
Thus, CasKR seems to be a dedicated 2CS among
B. cereus sensu lato strains that could participate in the low-temperature adaptation of both mesophilic and psychrotolerant strains, despite their different ranges of growth temperature. Interestingly,
casKR orthologs are absent from the genomes of strains of phylogenetic group VII, i.e., the most thermotolerant strains of
B. cereus sensu lato. Through this coincidence, it could be assumed that the lack of
casKR genes might be linked to the inability of strains from this group to grow at temperatures below 18°C. However, these strains also probably lack many of the tools needed for low-temperature growth, given that they have a fairly smaller genome than other
B. cereus sensu lato strains (
53). Although the chromosomal region of
casKR seems to have been conserved (see Fig. S4 in the supplemental material), these strains are also phylogenetically remote from other
B. cereus sensu lato strains and therefore constitute a distinct species (
B. cytotoxicus) in this group (
1).
Conclusion.
Food safety issues caused by
B. cereus occur because this food-borne bacterium is able to proliferate in food prior to ingestion by the consumer (
2). Despite the use of refrigeration processes to limit bacterial proliferation in food, some psychrotolerant strains of
B. cereus can still grow at these low temperatures. In addition, occasional leaks in the chill chain may also create conditions that allow the growth of mesophilic strains of
B. cereus. Understanding the mechanisms involved in the cold adaptation of the vegetative cells of
B. cereus may help to more accurately estimate the risk of proliferation in food. Efforts to identify such mechanisms have generally been performed in model mesophilic strains (
8,
11,
56). Here, we identified CasKR, a 2CS with a previously unknown function which seems to be a general mechanism among
B. cereus sensu lato strains that could participate in the low-temperature adaptation of both mesophilic and psychrotolerant strains. The detailed mechanism of this 2CS in cold sensing and adaptive response will have to be deciphered by further studies.