This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by van Schaik, W.
Right arrow Articles by Abee, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by van Schaik, W.
Right arrow Articles by Abee, T.

 Previous Article  |  Next Article 

Journal of Bacteriology, June 2007, p. 4384-4390, Vol. 189, No. 12
0021-9193/07/$08.00+0     doi:10.1128/JB.00313-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Identification of the {sigma}B Regulon of Bacillus cereus and Conservation of {sigma}B-Regulated Genes in Low-GC-Content Gram-Positive Bacteria{triangledown} ,{dagger}

Willem van Schaik,1,2,{ddagger} Menno van der Voort,1,2 Douwe Molenaar,1,3 Roy Moezelaar,1,4 Willem M. de Vos,1 and Tjakko Abee1,2*

Wageningen Centre for Food Sciences, P.O. Box 557, 6700AN Wageningen, The Netherlands,1 Laboratory of Food Microbiology, Wageningen University and Research Centre, P.O. Box 8129, 6700EV Wageningen, The Netherlands,2 NIZO Food Research, P.O. Box 20, 6710BA Ede, The Netherlands,3 Food Technology Centre, Wageningen University and Research Centre, P.O. Box 17, 6700AA Wageningen, The Netherlands4

Received 2 March 2007/ Accepted 29 March 2007


arrow
ABSTRACT
 
The alternative sigma factor {sigma}B has an important role in the acquisition of stress resistance in many gram-positive bacteria, including the food-borne pathogen Bacillus cereus. Here, we describe the identification of the set of {sigma}B-regulated genes in B. cereus by DNA microarray analysis of the transcriptome upon a mild heat shock. Twenty-four genes could be identified as being {sigma}B dependent as witnessed by (i) significantly lower expression levels of these genes in mutants with a deletion of sigB and rsbY (which encode the alternative sigma factor {sigma}B and a crucial positive regulator of {sigma}B activity, respectively) than in the parental strain B. cereus ATCC 14579 and (ii) increased expression of these genes upon a heat shock. Newly identified {sigma}B-dependent genes in B. cereus include a histidine kinase and two genes that have predicted functions in spore germination. This study shows that the {sigma}B regulon of B. cereus is considerably smaller than that of other gram-positive bacteria. This appears to be in line with phylogenetic analyses where {sigma}B of the B. cereus group was placed close to the ancestral form of {sigma}B in gram-positive bacteria. The data described in this study and previous studies in which the complete {sigma}B regulon of the gram-positive bacteria Bacillus subtilis, Listeria monocytogenes, and Staphylococcus aureus were determined enabled a comparison of the sets of {sigma}B-regulated genes in the different gram-positive bacteria. This showed that only three genes (rsbV, rsbW, and sigB) are conserved in their {sigma}B dependency in all four bacteria, suggesting that the {sigma}B regulon of the different gram-positive bacteria has evolved to perform niche-specific functions.


arrow
INTRODUCTION
 
The ability of bacteria to respond rapidly to changing environmental conditions is a prerequisite for survival in their habitats. This bacterial stress response is triggered by a change in the microorganism's growth conditions. Such a change triggers a cascade of events that will lead to the increased stress resistance of the bacterial cell, most often against not only the stress to which it was exposed but also other stresses, thereby ensuring its survival under a variety of conditions. A common strategy that bacteria use to counter stressful conditions is to activate a specific alternative sigma factor, which leads to the transcription of a set of genes (a so-called regulon), the products of which protect the cell against adverse conditions (17). In several gram-positive bacteria, the alternative sigma factor {sigma}B is the key sigma factor controlling the stress response (17, 26, 36).

In the gram-positive model organism Bacillus subtilis and the human pathogens Listeria monocytogenes and Staphylococcus aureus, numerous studies have shown that the activation of {sigma}B leads to an increased resistance to stress and that in pathogenic Listeria and Staphylococcus species, {sigma}B has a role in virulence (17, 36). In the bacteria that comprise the physiologically diverse Bacillus cereus group, the role of {sigma}B has also been studied. The most prominent members of the B. cereus group are B. anthracis (a mammalian pathogen that causes the disease anthrax), B. thuringiensis (an insect pathogen that is widely used as a biopesticide), and B. cereus, which is a common cause of food-borne infections (15). Recent genomics studies have revealed that these bacteria share a common core set of genes and that their phenotypic characteristics are determined largely by the presence or absence of large plasmids that may encode different virulence factors (30). Throughout the B. cereus group, the genes encoding {sigma}B and the regulators of its activity are conserved. In B. cereus, {sigma}B was shown to be activated under a number of stress conditions but most strongly upon a mild heat shock from 30°C to 42°C. Furthermore, the deletion of the sigB gene affected the heat-adaptive response of vegetative cells and some spore properties, most notably germination characteristics (4, 33). Also, in B. anthracis, {sigma}B was activated upon a heat shock, and the B. anthracis sigB deletion mutant was shown to have slightly lower virulence than the parental strain (8).

The regulation of {sigma}B activity in gram-positive bacteria occurs via a so-called partner-switching mechanism. Under nonstress conditions, {sigma}B is held in an inactive state by sequestration to an anti-sigma factor (RsbW). Under stress, an anti-sigma factor antagonist, RsbV, binds to RsbW, and this leads to the release of {sigma}B from its complex with RsbW. Subsequently, {sigma}B can bind to core RNA polymerase, leading to the transcription of {sigma}B-dependent genes. RsbW is not only the anti-sigma factor for {sigma}B, it also acts as a kinase on RsbV. In its phosphorylated form, RsbV is unable to bind to RsbW and activate {sigma}B. However, under stress conditions, the phosphorylated form of RsbV can be dephosphorylated by the action of a PP2C-type phosphatase, which then leads to the formation of the RsbV-RsbW complex and the activation of {sigma}B (reviewed in references 17, 26, and 36). In the B. cereus group, the PP2C phosphatase RsbY appears to be the only phosphatase responsible for the activation of {sigma}B under stress, as the deletion of the rsbY gene in B. cereus leads to an abolished {sigma}B response under several stress conditions, including a heat shock (34).

The identification of the complete {sigma}B regulons in B. subtilis (9, 24, 27), L. monocytogenes (16), and S. aureus (1, 22) by DNA microarray technology has revealed that a surprisingly low number of genes has an obvious role in the stress response of these organisms. Instead, {sigma}B-dependent specialized metabolic adaptations may contribute to stress resistance. Furthermore, in both L. monocytogenes and S. aureus, pathogenic traits seem to be governed by {sigma}B. These findings suggest that {sigma}B has evolved in response to specific niches that are populated by different bacteria.

Here, we describe the {sigma}B regulon of B. cereus, which was determined by comparing the transcriptome of the parental B. cereus strain ATCC 14579 and its sigB and rsbY deletion mutants upon a mild heat shock. It appears that the {sigma}B regulon of B. cereus is considerably smaller than those of other gram-positive bacteria, which may reflect the previously proposed close-to-ancestral state of {sigma}B in the B. cereus group (7). In other studies, the sets of {sigma}B-regulated genes in B. subtilis (9, 24, 27), L. monocytogenes (16), and S. aureus (1, 22) were determined by DNA microarray technology. This allowed us to compare the {sigma}B regulons of B. cereus, B. subtilis, L. monocytogenes, and S. aureus with each other to determine to what extent {sigma}B-regulated genes are conserved in the different bacteria and to assess how far the {sigma}B dependency of these genes is conserved.


arrow
MATERIALS AND METHODS
 
Strains, growth conditions, and RNA isolation. B. cereus strain ATCC 14579 and its isogenic derivatives, the sigB deletion mutant FM1400 (33) and the rsbY deletion mutant FM1401 (34), were used in this study. Strains were cultured overnight in 100-ml Erlenmeyer flasks containing 20 ml brain heart infusion broth at 30°C with aeration by rotary shaking at 150 rpm. Subsequently, prewarmed brain heart infusion broth (20 ml in 100-ml Erlenmeyer flasks) was inoculated with 100 µl of a single culture grown overnight. These cultures were then incubated at 30°C as described above. One of the three cultures was used to measure the absorbance at 600 nm (A600). When the A600 of this culture reached 0.4 (mid-exponential growth phase), one of the other two cultures was transferred into a shaking water bath at 42°C and incubated for a further 10 min, which was followed by RNA isolation. We previously found that at 10 min after the start of the heat shock, the levels of {sigma}B protein reached maximum levels (34), which leads to high-level expression of {sigma}B-dependent genes at this time point. RNA was isolated immediately from the culture in the mid-exponential growth phase. RNA isolation was performed by transferring the cultures into a 50-ml Falcon tube, and the cultures were spun down at 13,000 x g for 20 s. After decanting the supernatant, the cell pellets were snap-frozen in liquid nitrogen. The time between removal from the incubator and freezing of the cell pellets was approximately 60 s. Within 20 min after freezing the cell pellets, RNAwiz (Ambion, Huntingdon, United Kingdom) was added to the pellets, and RNA was extracted as described previously (33). Residual chromosomal DNA was removed by treating samples with DNA-free (Ambion). Extracted RNA samples were stored in 70% ethanol-83 mM sodium acetate buffer (pH 5.2) at –80°C.

cDNA synthesis and microarray design and hybridization. Before further analysis, the RNA samples were spun down at 13,000 rpm in a microtube centrifuge prechilled at 4°C. The RNA pellet was then washed with 70% ethanol (prechilled at –20°C) and resuspended in RNase-free water. Prior to cDNA synthesis, the quality of the extracted RNA was determined by analysis with an Agilent (Palo Alto, CA) 2100 Bioanalyzer according to the manufacturer's instructions. cDNA with amino-allyl-labeled dUTP (Ambion) was prepared in reverse transcription reactions using Superscript III (Invitrogen, Breda, The Netherlands). Cy3 and Cy5 labeling of the cDNAs was performed with the CyScribe Post-Labeling kit (GE Healthcare, Diegem, Belgium) as previously described (3). Labeled cDNAs were purified using the CyScribe GFX purification kit (GE Healthcare). Final cDNA yields and incorporation of the Cy labels were determined by Nanodrop analysis (Isogen Life Science, IJsselstein, The Netherlands). Custom-made B. cereus microarrays (produced by Agilent) (see below for details) were hybridized with 200 to 300 ng labeled cDNA. The microarray experiments for the comparison of the transcriptomes of the wild-type (WT) strain and the sigB and rsbY mutants were performed with three and two independent biological replicates, respectively. Microarray experiments (using two biological duplicates) were also performed to identify the set of genes that were differentially expressed in B. cereus ATCC 14579 after a heat shock from 30°C to 42°C for 10 min.

The microarrays used in this study were custom-made B. cereus ATCC 14579 microarrays using the 11K platform developed by Agilent Technologies. A total of 10,263 spots represented 5,058 B. cereus open reading frames, meaning that 96.6% of the predicted chromosomal B. cereus ATCC 14579 open reading frames (NCBI accession number NC_004722) are represented on the microarray. Ninety-five percent of the open reading frames for which probes could be designed were represented by two or three nonoverlapping probes on the array. The remaining 5% of the open reading frames were represented by a single oligonucleotide. After hybridization at 60°C for 17 h, the microarrays were washed according to the manufacturer's instructions with 6x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate)-0.005% Triton X-102 at room temperature for 10 min and with prechilled (at 4°C) 0.1x SSC-0.005% Triton X-102 for 5 min. Slides were then scanned using an Agilent microarray scanner (G2565BA), and data were extracted from the scanned microarrays with Feature Extraction software, which includes a Lowess normalization step for the raw data.

Analysis of microarray data. After removal of the data for the different controls printed on the microarray slides, the normalized data for each spot from the microarrays were analyzed for statistical significance using the Web-based VAMPIRE microarray suite (12, 13). A spot was found to be differentially expressed between two samples using the threshold of a false discovery rate smaller than 0.05. Subsequently, the data for the single spots were integrated to obtain expression ratios for an open reading frame. An open reading frame was found to be differentially expressed when all spots representing the open reading frame were significantly differentially expressed between samples. The average expression ratio was determined by calculating the average of the log values of the expression ratios over all experiments. This value (R) was then used to calculate the average expression ratio (10R). Finally, changes of 1.5-fold (for upregulated genes in the parental strain) and 0.67-fold (for downregulated genes in the parental strain) were also introduced as significance limits.

Validation of microarray data. The expression ratios of six genes were determined using quantitative real-time PCR (qPCR) with the oligonucleotides listed in Table 1 on the RNA samples isolated from strains ATCC 14579 and FM1400. cDNA synthesis and qPCR were performed as described previously, with tufA as a reference gene (34). All qPCR reactions were performed in triplicate. Expression ratios between the strains were determined by REST-tool (25).


View this table:
[in this window]
[in a new window]

  		TABLE 1. Oligonucleotides used in this study

Prediction of {sigma}B-dependent promoter sites. The prediction of {sigma}B-dependent promoter sites in B. cereus was performed by using the 500-bp upstream regions of the 26 genes identified in the {sigma}B regulon. Subsequently, the upstream regions were analyzed with AlignACE 3.1 (31) and DBTBS 4.1 (20) with a threshold of 5%. A total of 12 putative {sigma}B-dependent promoter sites were identified by these methods. These promoter sites were aligned with MUSCLE 3.6 (6), and a hidden Markov model (HMM) was constructed with the HMMER package (5). The HMM was used for a first screening for {sigma}B-dependent promoters in the genome sequence of B. cereus ATCC 14579 (14). Hits that were more than 300 nucleotides upstream of the start codon of an open reading frame or that had an E score higher than 1 x 10–3 were discarded. Two of the putative promoter sites identified in the first round of the HMMER search did not meet these standards. Therefore, an iteration of the alignment, the HMM, and the genome search were performed with the 10 remaining promoters. Predicted sites that were more than 300 nucleotides upstream of the start codon of an open reading frame or that had an E score higher than 1 x 10–3 were again discarded. The 10 {sigma}B-dependent promoter sites upstream of genes belonging to the {sigma}B regulon that were identified by this approach were aligned, and this alignment was visualized with WebLogo (2).

Comparison of the {sigma}B regulons of gram-positive bacteria. To determine the conservation of {sigma}B-dependent genes in gram-positive bacteria, the set of genes that is directly regulated by {sigma}B was compiled for B. cereus (this study), B. subtilis (9, 24, 27), L. monocytogenes (16), and S. aureus (1). Direct regulation was defined as genes for which {sigma}B dependency was determined by experimental confirmation (for example, by promoter mapping through primer extension analysis) or by the prediction of a {sigma}B-dependent promoter (for example, by the application of an HMM to identify {sigma}B-dependent promoters) upstream of the gene. To obtain the set of genes directly regulated by {sigma}B in B. subtilis, three different studies in which the {sigma}B regulon was determined were compared, and only genes that were found to be {sigma}B dependent in at least two of the three studies were used in this comparative analysis. For the definition of the genes that were directly regulated by {sigma}B in S. aureus, data from the study described previously by Bischoff et al. (1) were used. Only genes directly regulated by {sigma}B and that had a homolog in S. aureus N315 (18) were included in the analysis. The more limited study of the {sigma}B regulon described previously by Pane-Farre et al. (22) is not further discussed here, but a comparison of the results from this {sigma}B regulon with those of the other bacteria is included in the Table S4 in the supplemental material and is in line with the data described previously by Bischoff et al. (1).

TIGR's Multi-Genome Homology Comparison tool (http://cmr.tigr.org/tigr-scripts/CMR/shared/MakeFrontPages.cgi?page=circular_display), applying a P value of ≤1.0 x 10–5 for significance, was used for comparing the set of {sigma}B-dependent genes of each organism.

Microarray data accession number. Microarray data were submitted to the GEO database under accession number GSE6005.


arrow
RESULTS AND DISCUSSION
 
Validation of transcriptome data by qPCR. Because this is the first report of the use of B. cereus whole-genome microarrays, it was deemed essential to first validate the expression ratios generated by microarray analysis by qPCR. We used microarray data that were obtained in the comparison of the sigB mutant and its parental strain during mid-exponential growth at 30°C and upon a heat shock to 42°C for 10 min. We chose six genes with significantly different expression ratios under at least one growth condition, and the expression ratios found by microarray experiments exhibited a good correlation with the qPCR data (Fig. 1).


Figure 1
View larger version (9K):
[in this window]
[in a new window]

  		FIG. 1. Validation of microarray results by qPCR. Three independent RNA samples were isolated from cultures of B. cereus strain ATCC 14579 (WT) and the sigB deletion mutant during mid-exponential growth at 30°C (ME) and after a heat shock to 42°C for 10 min (HS). The expression ratios between both strains of the indicated genes were measured using microarrays (closed squares) and qPCR (open circles). Small symbols indicate that no significant difference was measured for the expression levels of the target gene between the two strains.

Identification of {sigma}B-dependent genes of B. cereus upon mild heat shock. To identify the set of {sigma}B-regulated genes in B. cereus, the transcriptomes of the sigB deletion mutant FM1400 and the rsbY deletion mutant FM1401 were compared with the transcriptome of the parental strain B. cereus ATCC 14579 after a heat shock to 42°C for 10 min. This mild heat shock strongly activates {sigma}B activity in the parental strain, while this response does not occur in both mutant strains (33, 34). Both the sigB and rsbY deletion mutants were used in these experiments to minimize the influence of mutant-specific, pleiotropic effects, which could give rise to false-positive results in the definition of the B. cereus {sigma}B regulon under this condition.

We observed that upon a heat shock to 42°C for 10 min, 45 genes were significantly upregulated in the parental strain in comparison with the sigB mutant (Fig. 2). Similarly, 38 genes were found to be upregulated after an identical heat shock in the WT strain compared to the rsbY mutant. Expression data for these experiments can be found in Table S1 in the supplemental material. A total of 28 genes from these two data sets were found to be overlapping, which suggests that these genes are candidates for being dependent on {sigma}B for their expression. No genes with expression ratios of ≤0.67 (WT/mutant) were found in either the sigB or the rsbY mutant, so neither {sigma}B nor {sigma}B-dependent proteins appear to have a role in the downregulation of genes in B. cereus upon a heat shock.


Figure 2
View larger version (9K):
[in this window]
[in a new window]

  		FIG. 2. Comparison of differentially expressed genes between B. cereus ATCC 14579 (WT) and the sigB and rsbY deletion mutants. The numbers of genes that are differentially expressed between the parental strain (WT) and the rsbY and sigB mutants upon a heat shock to 42°C for 10 min are indicated in the circles. The black overlapping area indicates that the same genes were identified as being differentially expressed in both the sigB and the rsbY deletion mutants under the indicated conditions. Twenty-four genes were also upregulated upon a heat shock, and these genes form the {sigma}B regulon of B. cereus.

A further defining feature of a {sigma}B-dependent gene is that the gene exhibits upregulated expression upon stress. Therefore, we subsequently performed microarray experiments to define the set of genes that are upregulated upon a heat shock from 30°C to 42°C for 10 min in the parental strain B. cereus ATCC 14579 and compared these data sets (see Table S2 in the supplemental material) with the above-described set of putative {sigma}B-dependent genes (Fig. 2). In total, 681 genes were upregulated upon a heat shock from 30°C to 42°C for 10 min, which corresponds to 12.4% of the total number of genes in B. cereus. This substantial transcriptional reprogramming corresponds to data for the heat shock response in B. subtilis, where over 10% of the total number of genes was induced upon a heat shock (9). The genes that were upregulated upon a heat shock in B. cereus include dnaK, groEL, groES, and the different clp genes, which is in line with data from a previous proteome-based study of the heat shock response of B. cereus (23). Twenty-four of the 28 previously identified genes are also part of the B. cereus heat stress stimulon, and these genes were thus assigned to the {sigma}B regulon of B. cereus (Table 2). Indeed, the expression ratios between mid-exponentially-growing and heat-shocked B. cereus ATCC 14579 (Table 2) are practically identical to the expression ratios when the mutant and parental strains were compared with one another upon a heat shock (Table 2). This reflects both the inactivity of {sigma}B under exponential growth conditions and the absence of a {sigma}B response in the mutants. Four genes (bc1436, bc5148, bc5149, and bc5150), which exhibited higher levels of expression in B. cereus ATCC 14579 than in the sigB and rsbY mutants, were not assigned to the {sigma}B regulon, as they were not upregulated in their expression upon a heat shock in the parental strain.


View this table:
[in this window]
[in a new window]

  		TABLE 2. Identification of {sigma}B-regulated genes in B. cereus ATCC 14579 upon mild heat shocka

Function of {sigma}B-dependent genes in B. cereus. The set of {sigma}B-dependent genes, which were identified by microarray analysis, includes genes that were previously identified in a study using a combined proteomics- and in vitro transcription-based approach (35), but a number of additional genes can now also be assigned to the {sigma}B regulon of B. cereus (Table 2). Note that two genes (bc0999 and bc1012) that were not represented on the DNA microarray are included in Table 2, as they were previously experimentally proven to be {sigma}B dependent (35). The newly identified genes include bc1007, which encodes a predicted CheR-type methyltransferase, and bc1008, a histidine kinase fused to a Chase3-like N-terminal domain (37) and a C-terminal CheY-type response regulator domain, which could be involved in a complex regulatory route in the B. cereus group. A large proportion of genes in the B. cereus {sigma}B regulon appears to have no known function. However, three of these genes have homologs in other bacteria, where they have a role in the stress response. These are bc0998, which has a homolog in B. subtilis (yflT), where it has a role in resistance to ethanol and hyperosmotic conditions (10); bc0999, which is homologous to csbD in B. subtilis and which has a role in hyperosmotic and cold stress (10); and bc1000, which is homologous to the GlsB protein of Enterococcus faecalis, where it has a role in resistance to bile (32). It is interesting that in the region delimited by bc0995 and bc1012 (a 12.4-kb fragment), {sigma}B is involved in the expression of 16 of the 18 open reading frames in this region. This region also includes the genes encoding all known regulators of the {sigma}B activity of B. cereus. Possibly, the clustering of both regulators and important {sigma}B-dependent genes can be advantageous, or it may be a remainder of an early event in the evolution of {sigma}B in gram-positive bacteria.

Another newly identified member of the {sigma}B regulon of B. cereus is bc2108, which encodes the extracytoplasmic function-type sigma factor {sigma}Z. This apparent regulation of the transcription of bc2108 by {sigma}B may very well be indirect, as no candidate {sigma}B promoter site could be identified upstream of bc2108. Two other genes (bc2638 and bc5391), which also appear to be indirectly dependent on {sigma}B for their expression, have predicted roles in the processes of germination and sporulation. The bc2638 gene is predicted to encode the C subunit of a germination receptor. These receptors are important for sensing nutrients and triggering the germination of the bacterial spore (21). The bc2638 gene is not part of one of the seven germination receptor operons of B. cereus (11), and conceivably, the encoded protein can act together with A and B subunits of other germination receptor complexes. The bc5391 gene also has a proposed role in germination. The B. subtilis homolog of this gene (named gerQ or ywdL) is essential for the presence of the cortex-lytic enzyme CwlJ in the spore coat (29). In B. subtilis, GerQ is cross-linked into high-molecular-mass forms by a transglutaminase. This may be important for building a spore coat that contributes to the resistance of the spore to environmental insults (28). Consequently, the germination defect and decreased heat resistance of the spores of the sigB deletion mutant of B. cereus (4) may be explained by the lower expression levels of bc2638 and bc5391 during sporulation.

Prediction of {sigma}B-dependent promoters. A HMMER search was performed using the genes of the {sigma}B regulon of B. cereus to identify {sigma}B-dependent promoters in the B. cereus ATCC 14579 genome (14). The HMMER search identified promoters regulating 21 of the 26 genes belonging to the {sigma}B regulon of B. cereus, suggesting that these genes are under the direct control of {sigma}B (Table 2). Promoters were identified in front of operons coding for 15 open reading frames of the island that are part of the {sigma}B regulon and in front of the bc0863-bc0862 and bc3132-bc3130 operons. Seven out of the 10 predicted promoter sites were experimentally defined previously (33, 35).

By aligning the predicted and experimentally verified {sigma}B-dependent promoter sites in B. cereus, a promoter consensus sequence for {sigma}B in B. cereus could be generated (Fig. 3). When the genome of B. cereus ATCC 14579 was searched with this consensus sequence, we found that the promoters of genes identified in this study were all found to have a low E score and therefore a significant fit to the HMM (see Table S3 in the supplemental material). Indeed, the 10 promoters with the lowest E scores (correlating with optimal fit to the HMM) were all identified in this study. Generally, the promoters with the best fit to the HMM have the highest change values (n-fold) in the microarray analysis, suggesting that differences between expression levels of the genes of the {sigma}B regulon are determined primarily by the recognition of the promoter sequence by {sigma}B. The other promoters are in all likelihood weakly (if at all) dependent on {sigma}B for their expression and/or may be expressed under other growth or stress conditions than those used in this study. Even though we cannot completely rule out that under other stress conditions, more {sigma}B-dependent genes may be identified, the fact that we identified the 10 promoters with the best fit to the HMM by microarray analysis implies that if additional {sigma}B-dependent genes can be identified under other conditions, their dependency on {sigma}B will probably be weak. Additionally, it is noteworthy that in our previous proteomics-based study, only a relatively small number of {sigma}B-dependent genes could be identified, even though we used an inducible overexpression system in B. cereus to obtain artificially high levels of {sigma}B and concomitant strong activation of {sigma}B-dependent genes (35). Both findings imply that the {sigma}B regulon of B. cereus is not much larger than approximately 25 to 30 genes.


Figure 3
View larger version (14K):
[in this window]
[in a new window]

  		FIG. 3. Consensus sequence of the {sigma}B-dependent promoter sites in B. cereus. The promoter sequences indicated in Table 2 were visualized with WebLogo (2). The height of the nucleotide is indicative of its frequency at that position in the promoter sites. Variable spacing between the –35 and –10 binding sites for {sigma}B is signified by the absence of a nucleotide at position 22.

Evolutionary implications of the small {sigma}B regulon of B. cereus. When the findings of this study are compared to the results for other gram-positive bacteria in which transcriptome profiling was used to determine the set of {sigma}B-regulated genes, it is striking that the {sigma}B regulon of B. cereus appears to be considerably smaller than those of other gram-positive bacteria. In L. monocytogenes, 54 genes were upregulated in a {sigma}B-dependent fashion upon stress exposure. However, the total number of {sigma}B-dependent genes in L. monocytogenes may be higher, as a microarray with limited genome coverage was used in that study (16). In both S. aureus and B. subtilis, the number of genes in the {sigma}B regulon is considerably larger and appears to be around 100 (1, 9, 22, 24, 26).

Possibly, the relatively small size of the {sigma}B regulon of B. cereus is a reflection of the way that {sigma}B has evolved throughout gram-positive bacteria. The similarities between the operon structures for sigB and sigF (encoding {sigma}F, the sigma factor that is active in the prespore) have been used to suggest that both operons have evolved from a common ancestor. Phylogenetic analysis indicates that the sigB operon of the B. cereus group most closely resembles the ancestral operon structure, whereas those of other gram-positive bacteria seem to have evolved further (7). Our finding of the relatively small size of the {sigma}B regulon of B. cereus seems to fit this hypothesis. While {sigma}B of B. cereus has a role in the development of stress resistance and is activated upon a number of different stress conditions (33), it does not seem to be very important for the regulation of metabolic pathways or the expression of virulence determinants as in the gram-positive pathogens L. monocytogenes and S. aureus (17). Possibly, the evolution of {sigma}B has progressed in a modular fashion in which functional groups of genes were added to the {sigma}B regulon in the different gram-positive bacteria based on selective pressures encountered in the different environments. In B. cereus, regulatory pathways other than {sigma}B may have evolved to regulate metabolic processes and the expression of virulence determinants, leaving B. cereus {sigma}B solely with the role of countering acute stresses and with a minor role in the process of sporulation.

Conservation of {sigma}B-dependent genes and {sigma}B-dependent gene expression in gram-positive bacteria. In this study, the {sigma}B regulon of B. cereus was determined by DNA microarray analysis. This allowed us to compare the set of {sigma}B-regulated genes in B. cereus with the previously published {sigma}B regulons of B. subtilis, L. monocytogenes, and S. aureus, which encompasses the phylogenetic diversity of gram-positive bacteria that have a bona fide sigB gene (7). We performed this analysis to determine to what extent {sigma}B-regulated genes are conserved among gram-positive bacteria and to assess how far the {sigma}B dependency of these genes is conserved (Fig. 4).


Figure 4
View larger version (18K):
[in this window]
[in a new window]

  		FIG. 4. Conservation of genes from the {sigma}B regulon in B. cereus, B. subtilis, L. monocytogenes, and S. aureus. The number of genes that are directly regulated by {sigma}B in a given organism is shown in the large circle. The conservation of {sigma}B-dependent genes from this organism in three other gram-positive bacteria (B. c., B. cereus; B. s., B. subtilis; L. m., L. monocytogenes; S. a., S. aureus) was determined with TIGR's Multi-Genome Homology Comparison tool (http://cmr.tigr.org/tigr-scripts/CMR/shared/MakeFrontPages.cgi?page=circular_display) with a cutoff for significance of a P value of ≤1.0 x 10–5. The numbers of conserved genes of the {sigma}B regulon are indicated in the arrows toward the small circles. The numbers in parentheses are the numbers of genes from the {sigma}B regulon of the query organism that have a homolog in the target organism. The second number in the arrow indicates the number of genes from this group, which are dependent on {sigma}B for their expression in both organisms. For comparative purposes, the conservation of all protein-coding genes between the genomes is indicated at the bottom of the figure.

This comparison showed that genes from the {sigma}B regulon are generally conserved among the four gram-positive bacteria to the same extent as all the genes of the genome. For example, of the 81 {sigma}B-dependent genes of B. subtilis, 58 genes (71.6%) have a homolog in B. cereus. When all 4,100 genes of B. subtilis are considered, 72.3% have a homolog in B. cereus. It should be noted, however, that {sigma}B-dependent genes that are strongly induced upon stress in B. subtilis, such as ctc, gsiB, and gspA, have no homologs in B. cereus. This already suggests that {sigma}B-dependent genes that are important for the stress response in one organism are not necessarily conserved in other organisms, even when these organisms are relatively closely related, like B. subtilis and B. cereus. It is also noteworthy that the genes of the {sigma}B regulon of L. monocytogenes are relatively poorly conserved in other bacteria. The peculiar lifestyle of L. monocytogenes as a gastrointestinal intracellular pathogen may have led to the divergent evolution of the {sigma}B regulon in this organism.

The most important conclusion that can be drawn from Fig. 4 is that even though {sigma}B-dependent genes from one organism are generally well conserved in other organisms, their {sigma}B dependency generally is not. This is clearly illustrated by comparing the {sigma}B regulons of the phylogenetically relatively closely related bacteria B. subtilis and B. cereus. Of the 81 {sigma}B-dependent genes of B. subtilis, 58 have a homolog in B. cereus. However, only eight of these {sigma}B-dependent B. subtilis genes are also {sigma}B dependent in B. cereus. The complete data set for the genes that have conserved their {sigma}B-dependent expression in the compared bacteria discussed here can be found in Table S4 in the supplemental material.

Only three genes are conserved in their {sigma}B dependency in all bacteria discussed here. These are the sigB gene itself and the genes encoding RsbV and RsbW, the anti-anti-sigma factor and anti-sigma factor of {sigma}B, respectively. The bc0999 gene from B. cereus also has {sigma}B-dependent homologs in B. subtilis (ywmG or csbD) and S. aureus (sa0772). It is also present in L. monocytogenes (lmo2518), but the {sigma}B dependency of this gene has not been described, possibly because this gene was not represented on the microarray used for the determination of the L. monocytogenes {sigma}B regulon (16). However, a putative {sigma}B-dependent promoter upstream of lmo2518 can be identified (data not shown), which strongly suggests that lmo2518 is dependent on {sigma}B for its expression. In that case, it would be the fourth gene that is {sigma}B dependent in all four different gram-positive bacteria, and this would be an indication that csbD-like genes may have an important role in the {sigma}B-mediated stress response. Interestingly, csbD-like genes are widely distributed in eubacteria, strongly suggesting that they perform an important, but hitherto-unrecognized, function in eubacteria.

The poorly conserved control of stress response genes by {sigma}B in the different gram-positive bacteria is in line with data reported in a recent article by Lozada-Chávez et al. (19). That study predicted that bacterial regulatory networks are extremely flexible during evolution and that transcription factors appear to be primarily responsible for the plasticity of transcriptional regulatory networks. The evolution of {sigma}B and its regulon in gram-positive bacteria appears to be an experimentally determined proof of this prediction. Even though {sigma}B's main function appears to be to control the stress response in different organisms, it appears that it does so by using a unique strategy in each gram-positive bacterium through controlling the expression of different output genes.


arrow
ACKNOWLEDGMENTS
 
We thank Jos Boekhorst and Mark W. de Been for their assistance in setting up tools for the analysis of the B. cereus microarrays.


arrow
FOOTNOTES
 
* Corresponding author. Mailing address: Laboratory of Food Microbiology, Wageningen University, Bomenweg 2, 6703HD Wageningen, The Netherlands. Phone: 31-317-484981. Fax: 31-317-484978. E-mail: tjakko.abee{at}wur.nl Back

{triangledown} Published ahead of print on 6 April 2007. Back

{dagger} Supplemental material for this article may be found at http://jb.asm.org/. Back

{ddagger} Present address: Toxines et Pathogénie Bacterienne, Institut Pasteur, Paris, France. Back


arrow
REFERENCES
 
    1
  1. Bischoff, M., P. Dunman, J. Kormanec, D. Macapagal, E. Murphy, W. Mounts, B. Berger-Bachi, and S. Projan. 2004. Microarray-based analysis of the Staphylococcus aureus {sigma}B regulon. J. Bacteriol. 186:4085-4099.[Abstract/Free Full Text]
    2. 2
  2. Crooks, G. E., G. Hon, J. M. Chandonia, and S. E. Brenner. 2004. WebLogo: a sequence logo generator. Genome Res. 14:1188-1190.[Abstract/Free Full Text]
    3. 3
  3. den Hengst, C. D., S. A. F. T. van Hijum, J. M. W. Geurts, A. Nauta, J. Kok, and O. P. Kuipers. 2005. The Lactococcus lactis CodY regulon: identification of a conserved cis-regulatory element. J. Biol. Chem. 280:34332-34342.[Abstract/Free Full Text]
    4. 4
  4. de Vries, Y. P., L. M. Hornstra, R. D. Atmadja, W. van Schaik, W. M. de Vos, and T. Abee. 2005. Deletion of sigB in Bacillus cereus affects spore properties. FEMS Microbiol. Lett. 252:169-173.[CrossRef][Medline]
    5. 5
  5. Durbin, R., S. R. Eddy, A. Krogh, and G. Mitchison. 1998. Biological sequence analysis: probabilistic models of proteins and nucleic acids. Cambridge University Press, Cambridge, United Kingdom.
    6. 6
  6. Edgar, R. C. 2004. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5:1-19.[Free Full Text]
    7. 7
  7. Ferreira, A., M. Gray, M. Wiedmann, and K. J. Boor. 2003. Comparative genomic analysis of the sigB operon in Listeria monocytogenes and in other gram-positive bacteria. Curr. Microbiol. 48:39-46.
    8. 8
  8. Fouet, A., O. Namy, and G. Lambert. 2000. Characterization of the operon encoding the alternative {sigma}B factor from Bacillus anthracis and its role in virulence. J. Bacteriol. 182:5036-5045.[Abstract/Free Full Text]
    9. 9
  9. Helmann, J. D., M. F. Wu, P. A. Kobel, F. J. Gamo, M. Wilson, M. M. Morshedi, M. Navre, and C. Paddon. 2001. Global transcriptional response of Bacillus subtilis to heat shock. J. Bacteriol. 183:7318-7328.[Abstract/Free Full Text]
    10. 10
  10. Hoper, D., U. Volker, and M. Hecker. 2005. Comprehensive characterization of the contribution of individual SigB-dependent general stress genes to stress resistance of Bacillus subtilis. J. Bacteriol. 187:2810-2826.[Abstract/Free Full Text]
    11. 11
  11. Hornstra, L. M., Y. P. de Vries, M. H. J. Wells-Bennik, W. M. de Vos, and T. Abee. 2006. Characterization of germination receptors of Bacillus cereus ATCC 14579. Appl. Environ. Microbiol. 72:44-53.[Abstract/Free Full Text]
    12. 12
  12. Hsiao, A., T. Ideker, J. M. Olefsky, and S. Subramaniam. 2005. VAMPIRE microarray suite: a Web-based platform for the interpretation of gene expression data. Nucleic Acids Res. 33:W627-W632.[Abstract/Free Full Text]
    13. 13
  13. Hsiao, A., D. S. Worrall, J. M. Olefsky, and S. Subramaniam. 2004. Variance-modeled posterior inference of microarray data: detecting gene-expression changes in 3T3-L1 adipocytes. Bioinformatics 20:3108-3127.[Abstract/Free Full Text]
    14. 14
  14. Ivanova, N., A. Sorokin, I. Anderson, N. Galleron, B. Candelon, V. Kapatral, A. Bhattacharyya, G. Reznik, N. Mikhailova, A. Lapidus, L. Chu, M. Mazur, E. Goltsman, N. Larsen, M. D'Souza, T. Walunas, Y. Grechkin, G. Pusch, R. Haselkorn, M. Fonstein, E. S. Dusko, R. Overbeek, and N. Kyrpides. 2003. Genome sequence of Bacillus cereus and comparative analysis with Bacillus anthracis. Nature 423:87-91.[CrossRef][Medline]
    15. 15
  15. Jensen, G. B., B. M. Hansen, J. Eilenberg, and J. Mahillon. 2003. The hidden lifestyles of Bacillus cereus and relatives. Environ. Microbiol. 5:631-640.[CrossRef][Medline]
    16. 16
  16. Kazmierczak, M. J., S. C. Mithoe, K. J. Boor, and M. Wiedmann. 2003. Listeria monocytogenes {sigma}B regulates stress response and virulence functions. J. Bacteriol. 185:5722-5734.[Abstract/Free Full Text]
    17. 17
  17. Kazmierczak, M. J., M. Wiedmann, and K. J. Boor. 2005. Alternative sigma factors and their roles in bacterial virulence. Microbiol. Mol. Biol. Rev. 69:527-543.[Abstract/Free Full Text]
    18. 18
  18. Kuroda, M., et al. 2001. Whole genome sequencing of methicillin-resistant Staphylococcus aureus. Lancet 357:1225-1240.[CrossRef][Medline]
    19. 19
  19. Lozada-Chávez, I., S. C. Janga, and J. Collado-Vides. 2006. Bacterial regulatory networks are extremely flexible in evolution. Nucleic Acids Res. 34:3434-3445.[Abstract/Free Full Text]
    20. 20
  20. Makita, Y., M. Nakao, N. Ogasawara, and K. Nakai. 2004. DBTBS: database of transcriptional regulation in Bacillus subtilis and its contribution to comparative genomics. Nucleic Acids Res. 32:D75-D77.[Abstract/Free Full Text]
    21. 21
  21. Moir, A., B. M. Corfe, and J. Behravan. 2002. Spore germination. Cell. Mol. Life Sci. 59:403-409.[CrossRef][Medline]
    22. 22
  22. Pane-Farre, J., B. Jonas, K. Forstner, S. Engelmann, and M. Hecker. 2006. The {sigma}B regulon in Staphylococcus aureus and its regulation. Int. J. Med. Microbiol. 296:237-258.[CrossRef][Medline]
    23. 23
  23. Periago, P. M., W. van Schaik, T. Abee, and J. A. Wouters. 2002. Identification of proteins involved in the heat stress response of Bacillus cereus ATCC 14579. Appl. Environ. Microbiol. 68:3486-3495.[Abstract/Free Full Text]
    24. 24
  24. Petersohn, A., M. Brigulla, S. Haas, J. D. Hoheisel, U. Volker, and M. Hecker. 2001. Global analysis of the general stress response of Bacillus subtilis. J. Bacteriol. 183:5617-5631.[Abstract/Free Full Text]
    25. 25
  25. Pfaffl, M. W., G. W. Horgan, and L. Dempfle. 2002. Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 30:e36.[Abstract/Free Full Text]
    26. 26
  26. Price, C. W. 2002. General stress response, p. 369-384. In A. L. Sonenshein, J. A. Hoch, and R. Losick (ed.), Bacillus subtilis and its closest relatives: from genes to cells. ASM Press, Washington, DC.
    27. 27
  27. Price, C. W., P. Fawcett, H. Ceremonie, N. Su, C. K. Murphy, and P. Youngman. 2001. Genome-wide analysis of the general stress response in Bacillus subtilis. Mol. Microbiol. 41:757-774.[CrossRef][Medline]
    28. 28
  28. Ragkousi, K., and P. Setlow. 2004. Transglutaminase-mediated cross-linking of GerQ in the coats of Bacillus subtilis spores. J. Bacteriol. 186:5567-5575.[Abstract/Free Full Text]
    29. 29
  29. Ragkousi, K., P. Eichenberger, C. van Ooij, and P. Setlow. 2003. Identification of a new gene essential for spore germination of Bacillus subtilis spores with Ca2+-dipicolinate. J. Bacteriol. 185:2315-2319.[Abstract/Free Full Text]
    30. 30
  30. Rasko, D. A., M. R. Altherr, C. S. Han, and J. Ravel. 2005. Genomics of the Bacillus cereus group of organisms. FEMS Microbiol. Rev. 29:303-329.[CrossRef][Medline]
    31. 31
  31. Roth, F. P., J. D. Hughes, P. W. Estep, and G. M. Church. 1998. Finding DNA regulatory motifs within unaligned noncoding sequences clustered by whole-genome mRNA quantitation. Nat. Biotechnol. 16:939-945.[CrossRef][Medline]
    32. 32
  32. Teng, F., E. C. Nannini, and B. E. Murray. 2005. Importance of gls24 in virulence and stress response of Enterococcus faecalis and use of the Gls24 protein as a possible immunotherapy target. J. Infect. Dis. 191:472-480.[CrossRef][Medline]
    33. 33
  33. van Schaik, W., M. H. Tempelaars, J. A. Wouters, W. M. De Vos, and T. Abee. 2004. The alternative sigma factor {sigma}B of Bacillus cereus: response to stress and role in heat adaptation. J. Bacteriol. 186:316-325.[Abstract/Free Full Text]
    34. 34
  34. van Schaik, W., M. H. Tempelaars, M. H. Zwietering, W. M. de Vos, and T. Abee. 2005. Analysis of the role of RsbV, RsbW, and RsbY in regulating {sigma}B activity in Bacillus cereus. J. Bacteriol. 187:5846-5851.[Abstract/Free Full Text]
    35. 35
  35. van Schaik, W., M. H. Zwietering, W. M. de Vos, and T. Abee. 2004. Identification of {sigma}B-dependent genes in Bacillus cereus by proteome and in vitro transcription analysis. J. Bacteriol. 186:4100-4109.[Abstract/Free Full Text]
    36. 36
  36. van Schaik, W., and T. Abee. 2005. The role of {sigma}B in the stress response of gram-positive bacteria—targets for food preservation and safety. Curr. Opin. Biotechnol. 16:218-224.[CrossRef][Medline]
    37. 37
  37. Zhulin, I. B., A. N. Nikolskaya, and M. Y. Galperin. 2003. Common extracellular sensory domains in transmembrane receptors for diverse signal transduction pathways in Bacteria and Archaea. J. Bacteriol. 185:285-294.[Abstract/Free Full Text]

Journal of Bacteriology, June 2007, p. 4384-4390, Vol. 189, No. 12
0021-9193/07/$08.00+0     doi:10.1128/JB.00313-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.



This article has been cited by other articles:

  • Wang, S.-W., Chen, C.-Y., Tseng, J. T., Liang, S.-H., Chen, S.-C., Hsieh, C., Chen, Y.-h., Chen, C.-C. (2009). orf4 of the Bacillus cereus sigB Gene Cluster Encodes a General Stress-Inducible Dps-Like Bacterioferritin. J. Bacteriol. 191: 4522-4533 [Abstract] [Full Text]  
  • den Besten, H. M. W., Mols, M., Moezelaar, R., Zwietering, M. H., Abee, T. (2009). Phenotypic and Transcriptomic Analyses of Mildly and Severely Salt-Stressed Bacillus cereus ATCC 14579 Cells. Appl. Environ. Microbiol. 75: 4111-4119 [Abstract] [Full Text]  
  • Emerson, J. E., Stabler, R. A., Wren, B. W., Fairweather, N. F. (2008). Microarray analysis of the transcriptional responses of Clostridium difficile to environmental and antibiotic stress. J Med Microbiol 57: 757-764 [Abstract] [Full Text]  

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by van Schaik, W.
Right arrow Articles by Abee, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by van Schaik, W.
Right arrow Articles by Abee, T.

Copyright © 2009 by the American Society for Microbiology.   For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»