Short Report | Open | Published:
Fatty acid profiles and desaturase-encoding genes are different in thermo- and psychrotolerant strains of the Bacillus cereus Group
BMC Research Notesvolume 8, Article number: 329 (2015)
The Bacillus cereus Group consists of closely-related bacteria, including pathogenic or harmless strains, and whose species can be positioned along the seven phylogenetic groups of Guinebretière et al. (I–VII). They exhibit different growth-temperature ranges, through thermotolerant to psychrotolerant thermotypes. Among these, B. cytotoxicus is an atypical thermotolerant and food-poisoning agent affiliated to group VII whose thermotolerance contrasts with the mesophilic and psychrotolerant thermotypes associated to the remaining groups I–VI. To understand the role of fatty acid (FA) composition in these variable thermotypes (i.e. growth behavior vs temperatures), we report specific features differentiating the FA pattern of B. cytotoxicus (group VII) from its counterparts (groups I–VI).
The FA pattern of thermotolerant group VII (B. cytotoxicus) displayed several specific features. Most notably, we identified a high ratio of the branched-chain FAs iso-C15/iso-C13 (i15/i13) and the absence of the unsaturated FA (UFA) C16:1(5) consistent with the absence of ∆5 desaturase DesA. Conversely, phylogenetic groups II–VI were characterized by lower i15/i13 ratios and variable proportions of C16:1(5) depending on thermotype, and presence of the DesA desaturase. In mesophilic group I, thermotype seemed to be related to an atypically high amount of C16:1(10) that may involve ∆10 desaturase DesB.
The levels of i15/i13 ratio, C16:1(5) and C16:1(10) UFAs were related to growth temperature variations recorded between thermotypes and/or phylogenetic groups. These FA are likely to play a role in membrane fluidity and may account for the differences in temperature tolerance observed in B. cereus Group strains.
In bacteria, the fatty acid (FA) composition of the cell membrane varies according to environmental conditions, as it is involved in bacterial adaptation to environmental changes such as temperature, pressure, and O2 availability [1–4]. Among these changes, the effect of temperature on bacterial FA composition is a prominent focus of research. Desaturases responsible for producing unsaturated FAs (UFAs) have been shown to play a role in low-temperature adaptation. The effect of incorporation of low-melting-point UFAs is to maintain membrane fluidity under the physical stress due to cold [5–9]. The number of desaturases varies depending on species, e.g. two desaturases have been identified in Bacillus cereus stricto or ss  whereas B. subtilis counts only one known desaturase . The FA composition of bacterial cells is also known to vary according to species, especially in Bacillus and related genera [12, 13], and has been included among important features for describing new taxa of the aerobic endospore-forming bacteria .
The Bacillus cereus Group (B. cereus sensu lato or sl) includes bacterial strains with a wide range of growth temperatures. These strains can be classified by growth-temperature range, from psychrotrophic to thermotolerant strains . These ranges of growth temperatures fit with the seven major phylogenetic groups (I–VII) established by Guinebretière et al. [15, 16] in the B. cereus Group (see Table 1). This 7-macrogroup classification is the most complete phylogenetic description of the B. cereus Group and is coherent with all MLST, AFLP, MLEE and genomic data produced since 2004 in the literature . It offers a unique setting to investigate the relation between temperature adaptation and hypothetical factors and can be used to resolve problems tied to effective species delimitation in the Group. Indeed, B. cereus sl contains seven closely-related species that, although not all genomospecies, can be clearly positioned by reference to each other in the classification of Guinebretière et al. (Additional file 1). Some of them have been described on the basis of singular phenotypic or pathologic traits such as rhizoidal colonies (B. mycoides), psychrotolerance (B. weihenstephanensis), insecticidal properties (B. thuringiensis) , enterotoxins (B. cereus sensu ss)  and anthrax (B. anthracis) . Only two are true genomospecies: the harmless species B. pseudomycoides (mesophilic group I) and the pathogenic species B. cytotoxicus (thermotolerant group VII) , both of which thus share a homogenous growth-temperature range. B. mycoides/B. weihenstephanensis (psychrotolerant group VI), and B. anthacis (a clonal lineage within the highly mesophilic group III) also have a homogenous growth-temperature range as they belong each to a unique phylogenetic group, whereas B. cereus ss and B. thuringiensis represent highly heterogeneous thermotypes across phylogenetic groups II–VI (Additional file 1) [15, 16].
Bacillus cytotoxicus, though relatively rare, is known as one of the most virulent pathogenic species of B. cereus sl [16, 21]. Its pathogenicity is mainly attributed to the greater expression and cytotoxic activity of the pore-forming cytotoxin K1 (CytK-1) [22–25], a variant of the cytotoxin K found in many strains of B. cereus sl . Bacillus cytotoxicus comprises solely thermotolerant strains [15, 21]. In addition to all these particularities, B. cytotoxicus has been described as displaying a specific FA composition in B. cereus sl . However, this specific FA composition of B. cytotoxicus has actually only been validated with reference to 4 of the 7 phylogenetic groups, and the relation between FA composition and the various thermotypes found in B. cereus sl has not yet been studied.
The aim of this study was to determine to what extent the FA composition of B. cytotoxicus (phylogenetic group VII) differs from that of all the other phylogenetic groups of the B. cereus Group, and to what extent these differences relate to the growth-temperature ranges of the groups. In addition, as desaturases are known to play a role in low-temperature adaptation through FA composition and membrane fluidity, we also investigated the presence of desaturase-encoding genes among the available B. cereus sl genomes and the existence of a putative relation between their presence and the FA composition of strains exhibiting different thermotypes.
The studied strains are listed in Additional file 2 and were representative of all seven phylogenetic groups. Minimal growth temperature (Tmin) and maximal growth temperature (Tmax) were used to determine growth-temperature range (Tmin–Tmax). Their group affiliation and growth-temperature range were determined in previous studies in 2008 and 2013 [15, 21]. As a rule, each phylogenetic group has its own growth-temperature range , as presented in Table 1. In the previous work of 2008, Tmin and Tmax were determined using a standard test described in the Bergey’s Manual , with temperature fluctuation of the incubators being not greater than ±0.2°C for all tested temperatures.
All these strains are referenced in the previous works of 2008 and 2013 [15, 21].These strains have since been conserved in our Laboratory (UMR408 Collection) and the original sources are presented in Additional file 2.
The thermotypes were determined from the growth-temperature ranges (Tmin–Tmax) as the resulting phenotype, and are presented in Table 1 by reference to growth-temperature range.
The FA methyl esters (FAMEs) were extracted using the standardized MIDI protocol (http://www.microbialid.com/PDF/TechNote_101.pdf). Bacterial cells were obtained from culturing at 30°C (at or close to optimal conditions for all seven phylogenetic groups) on trypticase soy broth agar (TSBA, 30 g trypticase soy broth and 15 g Bacto agar; l−1) for 24 h as previously described . After extraction, FAMEs were analyzed by gas chromatography-mass spectrometry (GC–MS) (Shimadzu QP2010 system), as previously described .
In silico analysis
A total of 210 B. cereus sl genomes available in databases at the time of the search were used for this study. Genome affiliations to the phylogenetic groups were established previously , as described in , using panC sequence similarity.
Two desaturases have been described in B. cereus, i.e. DesA and DesB , which are responsible for two different types of unsaturation. The DesA enzyme adds a double bond in the ∆5-position of a saturated FA (SFA) while DesB creates an unsaturation in the ∆10-position of an SFA. The presence of each desaturase-encoding gene was thus investigated among all the available B. cereus sl genomes. The search for desA and desB orthologs with reference to locus_tag BC_2983 and BC_0400, respectively, in the ATCC 14579T genome (i.e. sequence loci described for desaturases in the ATCC 14579 genome) was performed via the Integrated Microbial Genomes (IMG) interface . No other ∆5 or ∆10 desaturase-encoding genes have been reported in the B. cereus Group. First, candidate homologs were identified based on BLASTp similarities with a 1e-2 E-value cutoff and with low-complexity soft masking (-F’m S’) turned on. Second, orthologous relationships between BC_2983 or BC_0400 genes and their respective homologous genes in all other genomes were established through bidirectional best blast hits.
Results and discussion
FA composition in B. cereus sensu lato: general features
The B. cereus Group displays a specific overall FA pattern setting it apart from the other species of the Bacillus genus [21, 30], with short-chain branched FAs (12C and 13C) and a characteristic predominance of iso-C13:0 (Additional file 2). Whatever the phylogenetic group analyzed, three major FAs were identified: iso-C15:0 (i15), iso-C13:0 (i13) and C16:0 (n16). For better visibility, only these three major SFAs and the 7 UFAs (previously linked to cold adaptation [3, 27]) are listed in Table 1. While the n16 SFA did not range widely according to thermotype, the proportion of the other two SFAs (i13 and i15) varied for the most sharply-contrasting thermotypes, i.e. the thermotolerant type (group VII, B. cytotoxicus), the highly mesophilic type (group III), and the highly psychrotolerant type (group VI). The SFA i15 accounted for more than 1/13.3 of total FAs in group VII (B. cytotoxicus) and thus constituted a marker for this group.
i15/i13 ratio as a rough indicator of thermotype
We calculated the i15/i13 ratio, defined as proportion of iso-C15:0 divided by proportion of iso-C13:0, for strains belonging to each phylogenetic group (Fig. 1a). i15/i13 ratio was specifically very high in the thermotolerant strains of group VII (7.13 ± 1.7, Fig. 1a) but comparatively very low in the other groups (I–VI), with some variations: largely above 1 (1.57 ± 0.04) for the highly mesophilic strains of group III, largely below 1 (0.61 ± 0.13) for the highly psychrotolerant strains of group VI, and close to 1 in the remaining strains (groups I, II, V and IV). Taken together, these data strongly suggest that the i15/i13 ratio criterion can roughly discriminate B. cereus sl strains into 4 categories (Table 1): (1) highly psychrotolerant strains (group VI), (2) low psychrotolerant to mesophilic strains (groups II, V, IV, I), (3) highly mesophilic strains (group III), and (4) thermotolerant strains (group VII, B. cytotoxicus).
Presence of ∆5/∆10-desaturases and proportion of C16:1(5) as an accurate indicator of thermotype in the most recent branch of B. cereus sl
UFAs unsaturated at the ∆5 position were present in the most recent phylogenetic groups (II–VI, see Additional file 1) yet absent from thermotolerant group VII (B. cytotoxicus) and mesophilic group I (Table 1), particularly C16:1(5). We therefore searched for the presence of the genes encoding for the desaturases responsible for the synthesis of UFAs unsaturated at the ∆5 or ∆10 location (desA and desB respectively) in this group of bacteria. The results indicated that, contrary to desB, orthologs of desA were not found in all B. cereus sl genomes (Table 1): none of the genomes in groups I and VII contained an ortholog of desA gene. The 2.5–4% of negative genomes in groups III and IV may be due to information missing from draft genomes. However, in groups I (B. pseudomycoides strains) and VII (B. cytotoxicus strains), the ortholog of desA was truly absent, consistent with the near-zero concentration of ∆5 UFAs observed in groups I and VII.
Considering only groups II to VI, C16:1(5) proportion increased with thermotype, from the highly mesophilic to the highly psychrotolerant groups (with no significant difference between the two psychrotolerant groups II and VI) (Fig. 1b). Despite being produced in low amounts, C16:1(5) proportion appeared a good parameter to discriminate psychrotolerant thermotypes (groups II, VI) from other thermotypes and even to discriminate between mesophilic thermotypes (groups III, IV, V), including those that were not discriminated by i15/i13 ratio (groups II, V, IV).
Atypically in group-I strains, UFAs unsaturated at the ∆10 location were produced in higher proportions than in the other groups (Table 1), particularly C16:1(10), offsetting the lack of UFAs unsaturated at the ∆5 location. This difference presumably contributes to membrane fluidity and allows group I (B. pseudomycoides) to exhibit a mesophilic thermotype. In contrast to mesophilic group I (B. pseudomycoides), the thermotolerant group VII (B. cytotoxicus) seems unable to offset this deficiency by producing larger amounts of UFAs unsaturated at the ∆10 location.
FA composition and putative link with cold adaptation in the B. cereus Group
Our study highlighted the relation between i15/i13 ratio and B. cereus Group thermotypes. Another study reported that proportion of i13 was strongly reduced in a B. cereus ATCC 14579 mutant displaying growth impairment at low temperature compared with its parental strain during growth at low temperature . The i15/i13 ratio recalculated from these data  was much higher in the mutant than in the parental strain, emphasizing the putative role of a low i15/i13 ratio for psychrotolerance ability.
The C16:1(5) UFA presumably appeared with ∆5 desaturase DesA in the most recent branch of the phylogeny in B. cereus sl containing groups II–VI (see Additional file 1). Its absence in phylogenetic groups I and VII indicates a link with the whole evolution of the B. cereus Group. Indeed, groups I and VII belong to two other independent branches at the base of the phylogenetic tree (see Additional file 1). Taken together, these results converge towards a differential process of evolution involving ∆5 UFAs for groups II–VI and ∆10 UFAs for group I. As group VII is basal to the evolutionary tree, followed by group I and then the remaining groups, we can posit that the ancestor of the B. cereus Group was devoid of ∆5 UFAs and went through adaptation in a few steps. The first step would involve an increase of ∆10 UFAs through group I. The second step would involve ∆5 UFAs through groups II–VI, with a more efficient adaptation from mesophily to psychrotolerance. This is also consistent with the absence of ∆5 desaturase (DesA) in groups I and VII. Interestingly, the same kind of configuration (absence in groups I and VII) was also observed for the two-component system CasK/R, which was recently described as playing a role in B. cereus-Group cold adaptation . Through this observation, we can also posit that the lack of key genes such as casK/R and desA might be related to the inability of B. cytotoxicus strains to grow at temperatures below 20°C , and that these genes probably took part in a more complex mechanism of adaptation in the most recent branch of the phylogeny.
A link was established between the FA pattern of B. cereus sl strains and ability or inability to grow at low temperature. The FA profile of B. cytotoxicus (group VII) is highly specific compared to that of the phylogenetic groups I–VI and is relatable to its atypical thermotolerance: high i15/i13 ratio, absence of UFAs unsaturated at the ∆5 location (particularly C16:1(5)), absence of a ∆5 desaturase (DesA). In contrast, the presence of ∆5 desaturase DesA and subtle amounts of C16:1(5) seem to be associated with an advanced mechanism of adaptation, resulting in a large panel of thermotypes through groups II–VI (from mesophily to psychrotolerance). Mesophilic strains of Group I seem to exhibit a specific intermediate state of evolution involving a fairly atypical amount of the ∆10-desaturated C16:1(10).
Availability of supporting data
The datasets supporting the results of this article are included in Additional file 2.
- ΔX desaturase:
a desaturase adding a double bond on the Xth carbon–carbon position of the fatty acid, counting from the carboxyl acid end
anteiso-branched fatty acid with x number of carbons
- Cx: y (z):
fatty acid with x number of carbons and y number of unsaturation(s) at location z
fatty acid methyl ester
intermediate between mesophilic and psychrotolerant
iso-branched fatty acid with x number of carbons
saturated fatty acid with x number of carbons
- ss :
- sl :
unsaturated fatty acid
x carbon molecule(s)
Yano Y, Nakayama A, Ishihara K, Saito H (1998) Adaptive changes in membrane lipids of barophilic bacteria in response to changes in growth pressure. Appl Environ Microbiol 64:479–485
Sinensky M (1971) Temperature control of phospholipid biosynthesis in Escherichia coli. J Bacteriol 106:449–455
de Sarrau B, Clavel T, Clerte C, Carlin F, Ginies C, Nguyen-The C (2012) Influence of anaerobiosis and low temperature on Bacillus cereus growth, metabolism, and membrane properties. Appl Environ Microbiol 78:1715–1723
Beranova J, Mansilla MC, de Mendoza D, Elhottova D, Konopasek I (2010) Differences in cold adaptation of Bacillus subtilis under anaerobic and aerobic conditions. J Bacteriol 192:4164–4171
Fulco AJ (1967) The effect of temperature on the formation of delta 5-unsaturated fatty acids by bacilli. Biochim Biophys Acta 144:701–703
Fulco AJ (1972) The biosynthesis of unsaturated fatty acids by Bacilli: III. Uptake and utilization of exogenous palmitate. J Biol Chem 247:3503–3510
Fulco AJ (1969) The biosynthesis of unsaturated fatty acids by Bacilli: I. Temperature induction of the desaturation reaction. J Biol Chem 244:889–895
Kaneda T (1972) Positional preference of fatty acids in phospholipids of Bacillus cereus and its relation to growth temperature. Biochim Biophys Acta 280:297–305
Bredeston LM, Marciano D, Albanesi D, De Mendoza D, Delfino JM (2011) Thermal regulation of membrane lipid fluidity by a two-component system in Bacillus subtilis. Biochem Mol Biol Educ 39:362–366
Chazarreta Cifre L, Alemany M, de Mendoza D, Altabe S (2013) Exploring the biosynthesis of unsaturated fatty acids in Bacillus cereus ATCC 14579 and functional characterization of novel acyl-lipid desaturases. Appl Environ Microbiol. 79:6271–6279
Aguilar PS, Hernandez-Arriaga AM, Cybulski LE, Erazo AC, de Mendoza D (2001) Molecular basis of thermosensing: a two-component signal transduction thermometer in Bacillus subtilis. EMBO J 20:1681–1691
Kampfer P (1994) Limits and possibilities of total fatty acid analysis for classification and identification of Bacillus species. Syst Appl Microbiol 17:86–98
Kampfer P (2002) Whole-cell fatty acid analysis in the systematics of Bacillus and related genera. In: Berkeley R, Heyndrickx M, Logan N, DeVos P (eds) Applications and Systematics of Bacillus and Relatives, pp 271–299
Logan NA, Berge O, Bishop AH, Busse HJ, De Vos P, Fritze D et al (2009) Proposed minimal standards for describing new taxa of aerobic, endospore-forming bacteria. Int J Syst Evol Microbiol 59:2114–2121. doi:10.1099/ijs.0.013649-0
Guinebretière M-H, Thompson FL, Sorokin A, Normand P, Dawyndt P, Ehling-Schulz M et al (2008) Ecological diversification in the Bacillus cereus Group. Environ Microbiol 10:851–865. doi:10.1111/j.1462-2920.2007.01495.x
Guinebretiere MH, Velge P, Couvert O, Carlin F, Debuyser ML, Nguyen-The C (2010) Ability of Bacillus cereus Group strains to cause food poisoning varies according to phylogenetic affiliation (groups I to VII) rather than species affiliation. J Clin Microbiol 48:3388–3391
Tourasse NJ, Helgason E, Klevan A, Sylvestre P, Moya M, Haustant M et al (2011) Extended and global phylogenetic view of the Bacillus cereus group population by combination of MLST, AFLP, and MLEE genotyping data. Food Microbiol 28:236–244. doi:10.1016/j.fm.2010.06.014
Bravo A, Likitvivatanavong S, Gill SS, Soberon M (2011) Bacillus thuringiensis: a story of a successful bioinsecticide. Insect Biochem Mol Biol 41:423–431
Stenfors Arnesen L, Fagerlund A, Granum P (2008) From soil to gut: Bacillus cereus and its food poisoning toxins. FEMS Microbiol Rev 32:579–606
Kalamas AG (2004) Anthrax. Anesthesiol Clin N Am 22:533–540
Guinebretiere MH, Auger S, Galleron N, Contzen M, De Sarrau B, De Buyser ML et al (2013) Bacillus cytotoxicus sp. nov. is a novel thermotolerant species of the Bacillus cereus Group occasionally associated with food poisoning. Int J Syst Evol Microbiol 63:31–40
Lund T, De Buyser M-L, Granum PE (2000) A new cytotoxin from Bacillus cereus that may cause necrotic enteritis. Mol Microbiol 38:254–261. doi:10.1046/j.1365-2958.2000.02147.x
Fagerlund A, Ween O, Lund T, Hardy SP, Granum PE (2004) Genetic and functional analysis of the cytK family of genes in Bacillus cereus. Microbiology 150:2689–2697
Fagerlund A, Brillard J, Furst R, Guinebretiere MH, Granum PE (2007) Toxin production in a rare and genetically remote cluster of strains of the Bacillus cereus group. BMC Microbiol 7:43
Brillard J, Lereclus D (2004) Comparison of cytotoxin cytK promoters from Bacillus cereus strain ATCC 14579 and from a B. cereus food-poisoning strain. Microbiology-Sgm. 150:2699–2705. doi:10.1099/mic.0.27069-0
Claus D, Berkeley RCW (1986) Genus Bacillus Cohn 1872, 174AL. In: Sneath PHA, Mair NS, Sharpe ME, Holt JG (eds) Bergey’s manual of systematic bacteriology. Williams & Wilkins, Baltimore, pp 1105–1139
Brillard J, Jehanno I, Dargaignaratz C, Barbosa I, Ginies C, Carlin F et al (2010) Identification of Bacillus cereus genes specifically expressed during growth at low temperatures. Appl Environ Microbiol 76:2562–2573
Diomandé S, Chamot S, Antolinos V, Vasai F, Guinebretière M-H, Bornard I et al (2014) The CasKR two-component system is required for growth at low temperature of mesophilic and psychrotolerant Bacillus cereus strains. Appl Environ Microbiol. doi:10.1128/aem.00090-14
Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Grechkin Y et al (2012) IMG: the Integrated Microbial Genomes database and comparative analysis system. Nucleic Acids Res. 40:D115–D122
Song Y, Yang R, Guo Z, Zhang M, Wang X, Zhou F (2000) Distinctness of spore and vegetative cellular fatty acid profiles of some aerobic endospore-forming bacilli. J Microbiol Methods 39:225–241
SD organized and drafted the manuscript, analyzed the acquired data, performed the statistical analysis, interpreted the results and contributed to the write-up of the results. MHG conceived the study design, performed FAMEs extraction, contributed to data analysis and interpretation of results, and organized and edited the manuscript. BDS performed GC–MS analysis and extracted the data. CNT and VB were involved in editing and writing up the final version of the manuscript. JB supervised the analysis and interpretation of the outcomes and was involved in writing up the different versions of the manuscript. All the authors read and approved the final manuscript.
The authors thank the INRA-MICA department and Provence-Alpes-Côte d’Azur Regional Council for providing financial support for S.D.’s PhD.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Phylogeny of the B. cereus Group based on rrs and ITS1 gene sequence analysis, and position of species in this phylogeny. This file illustrates the position of the B. cereus sl species over 7 major phylogenetic groups according to the classification scheme given in Guinebretière et al. 2008. It also shows the three main branches of evolution in the B. cereus Group.
Fatty acid (FA) composition of B. cereus sl strains belonging to the 7 phylogenetic groups. This file contains the raw FA composition data of the B. cereus sl strains used to generate Table 1 and Figure 1.