# Short-Chain and Unsaturated Fatty Acids Increase Sequentially From the Lag Phase During Cold Growth of Bacillus cereus

Marina Français, Romain Bott, Claire Dargaignaratz, Christian Giniès, Frédéric Carlin, Véronique Broussolle, Christophe Nguyen-Thé

<figure><img src="image_1.png" /></figure>

## To cite this version:

Marina Français, Romain Bott, Claire Dargaignaratz, Christian Giniès, Frédéric Carlin, et al.. Short-Chain and Unsaturated Fatty Acids Increase Sequentially From the Lag Phase During Cold Growth of Bacillus cereus. Frontiers in Microbiology, 2021, 12, pp.694757. 10.3389/fmicb.2021.694757. hal-03313562

HAL Id: hal-03313562

https://hal.inrae.fr/hal-03313562

Submitted on 4 Aug 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Short-Chain and Unsaturated Fatty Acids Increase Sequentially From the Lag Phase During Cold Growth of *Bacillus cereus*

Marina Français, Romain Bott, Claire Dargaignaratz, Christian Giniès, Frédéric Carlin, Véronique Broussolle and Christophe Nguyen-Thé*

INRAE, Avignon Université, UMR SQPOV, Avignon, France

Fatty acids of two mesophilic and one psychrotrophic strains of the foodborne pathogen *Bacillus cereus* were analyzed by gas chromatography coupled to mass spectrometry during growth at cold (10 and 12°C) vs. optimal (30°C) temperatures and during the whole growth process (6–7 sampling times) from lag to stationary phase. In all these strains, a sequential change of fatty acids during cold growth was observed. Fatty acids were modified as soon as the end of lag, with an increase of the short-chain fatty acids (less than 15 carbons), particularly i13. These short-chain fatty acids then reached a maximum at the beginning of growth and eventually decreased to their initial level, suggesting their importance as a rapid cold adaptation mechanism for *B. cereus*. In a second step, an increase in $\Delta^{5,10}$ di-saturated fatty acids and in monounsaturated fatty acids in $\Delta^{5}$ position, at the expense of unsaturation in $\Delta^{10}$, started during exponential phase and continued until the end of stationary phase, suggesting a role in growth consolidation and survival at cold temperatures. Among these unsaturated fatty acids, those produced by unsaturation of n16 increased in the three strains, whereas other unsaturated fatty acids increased in some strains only. This study highlights the importance of kinetic analysis of fatty acids during cold adaptation.

**Keywords:** fatty acids, *Bacillus cereus*, cold adaptation, lag phase, kinetic, short-chain fatty acids, unsaturated fatty acids

**OPEN ACCESS**

**Edited by:**
Christian Sohlenkamp,
National Autonomous University
of Mexico, Mexico

**Reviewed by:**
Yu Sang,
Oregon Health and Science
University, United States
Marta Palusinska-Szysz,
Marie Curie-Sklodowska University,
Poland

*Correspondence:
Christophe Nguyen-Thé
christophe.nguyen-the@inrae.fr

**Specialty section:**
This article was submitted to
Microbial Physiology and Metabolism,
a section of the journal
Frontiers in Microbiology

**Received:** 13 April 2021

**Accepted:** 28 June 2021

**Published:** 22 July 2021

**Citation:**
François M, Bott R, Dargaignaratz C, Giniès C, Carlin F, Broussolle V and Nguyen-Thé C (2021) Short-Chain and Unsaturated Fatty Acids Increase Sequentially From the Lag Phase During Cold Growth of *Bacillus cereus*. *Front. Microbiol.* 12:694757.
doi: 10.3389/fmicb.2021.694757

## INTRODUCTION

*Bacillus cereus sensu lato* (*B. cereus sl*), a major cause of foodborne outbreaks in Europe (EFSA-BIOHAZ, 2016), is composed of closely related species (*B. cereus sensu stricto*, *Bacillus thuringiensis*, *Bacillus cytotoxicus*, *Bacillus weihenstephanensis*, *Bacillus wiedmanii*, *Bacillus toyonensis*, *Bacillus mycoides*, and *Bacillus pseudomycoides*) able to grow over a wide range of temperatures and in various environments (Guinebretiere et al., 2008). *B. cereus sl* is divided into seven phylogenetic groups, including the species mentioned earlier, according to their growth temperature range, from psychrotrophic (able to grow at 4°C) to nearly thermophilic groups (with a maximum growth temperature of 55°C). Adaptation to cold temperatures and the formation of highly resistant spores (Luu-Thi et al., 2014) make *B. cereus* an important hazard for heat-processed and refrigerated foods (EFSA, 2005). Inappropriate storage practices with exposure at abuse temperature or improper cooling, regularly observed at consumers' home or in catering, for instance, are also opportunities for multiplication of mesophilic strains of *B. cereus sl* (ANSES, 2017).
B. cereus cold adaptation relies on multiple mechanisms (Brillard and Broussolle, 2012), such as glucose metabolism (Chung et al., 1976), over-synthesis of specific proteins such as cold acclimation proteins, including RNA-helicases (Broussolle et al., 2010; Pandiani et al., 2010, 2011), cold-shock proteins, DNA gyrase (Beckering et al., 2002), or activation of the CasKR two-component system (Diomande et al., 2014). A crucial mechanism for bacterial growth at cold temperatures is also change in fatty acid (FA) composition of the membrane to reduce its melting point, maintain fluidity, and hence exchanges with the extracellular environment (Haque and Russell, 2004; Russell, 2008; de Sarrau et al., 2012). The most common FAs of bacteria, including B. cereus, are straight-chain FAs, branched-chain FAs (BCFAs) in iso or anteiso position, saturated FAs (SFAs), and unsaturated FAs (UFAs) (Diomande et al., 2015c). In most bacteria, the melting point of membrane FAs at cold temperatures is lowered by the synthesis of short-chain FAs (SCFAs) (Kaneda, 1972), increase in UFAs relative abundance, and increase in anteiso/iso BCFAs ratio (Kaneda, 1977; Freese et al., 2008). Mesophilic B. cereus growth at relatively cold temperatures (12 or 15°C) was characterized by an increase in UFA and SCFA relative abundance (Haque and Russell, 2004; Brillard et al., 2010; de Sarrau et al., 2012; Chazarreta Cifre et al., 2013; Diomande et al., 2015a). However, the kinetic of these FA changes in B. cereus during growth is not known. In other bacterial species, a few studies investigated FA changes during the stationary phase (Annous et al., 1999) and during the whole growth at optimal and high temperatures (Lopes et al., 2019a,b), but none considered the whole cold adaptation and cold growth process, from lag to late stationary.

In particular, little is known on what happens during the early phase of growth at cold temperatures, when the bacteria has to adapt to its new growth environment. This hampers an accurate and robust estimation of growth and risk in refrigerated foods (Bertrand, 2019). This is particularly true for FA composition, which has so far been studied at a single phase of cold growth, either exponential or stationary. A dynamic view of FA modifications during cold temperatures adaptation is therefore lacking.

This work aims to examine FAs' relative abundance during the whole growth process of B. cereus, from lag to stationary phase, to detect critical stages in changes of FA profiles, considering different B. cereus strains presenting *de facto* contrasted capacity to grow at cold temperatures. Consequently, we monitored the FA composition of two mesophilic and one psychrotrophic B. cereus strains, all along with growth, at 10, 12, and 30°C.

# MATERIALS AND METHODS

## Bacterial Strains and Media

B. cereus sl strains American Type Culture Collection (ATCC) 14579T and ATCC 10876 (both mesophilic strains, from B. cereus phylogenetic group IV, able to grow between 10 and 45°C) and MM3 (a psychrotrophic strain, from phylogenetic group II, growing between 7 and 40°C and recently affiliated to B. wiedmanii; Miller et al., 2016) were used. Strains were stored

at -80°C in a 30% glycerol solution. For experimental purpose, B. cereus cells were grown in mAOAC made of synthetic broth AOAC (Wright and Mundy M334 broth, HiMedia Laboratories) at pH 7.0, sterilized by autoclaving at 121°C for 20 min and supplemented with a filter-sterilized glucose solution to a final concentration of 6.9 mM (Diomande et al., 2015a; Français et al., 2019).

## Growth Experiments

One purified colony of a 24-h culture was grown overnight at 30°C in mAOAC broth for each tested B. cereus strain. An aliquot of each of these cultures was diluted into 100 ml of mAOAC medium to obtain an A600 of 0.1. The inoculated medium was incubated at 30°C, under shaking at 200 rpm until an A600 of 0.5 and was used as inoculum for growth experiments.

Growth experiments were done at 30°C for warm temperatures and at 10 and 12°C for cold temperatures. Thirty degrees Celsius represented conditions close to the optimum for both the mesophilic and psychrotrophic strains (Diomande et al., 2015b). Twelve degrees Celsius represented a cold temperature that permitted a substantial growth of the mesophilic strains in the AOAC medium (Diomande et al., 2015a). Ten degrees Celsius was used for the psychrotrophic strain to have the same interval from the minimum growth temperature as for the mesophilic strains.

For each strain, the inoculum was 10-fold diluted in 500 ml of mAOAC. Cultures were then incubated in a thermostatically controlled incubator at the tested temperature (30, 12, or 10°C) under shaking at 200 rpm. Colony-forming unit (CFU) counts and A600 were regularly measured to monitor growth. CFU counts were obtained by plating serial dilutions on Luria Bertani (Biokar) agar plates incubated overnight at 30°C. Lag times were determined by fitting the growth curves to the Baranyi model using the Combase online tool.<fnref n="1" /> Additional growth was done in microtiter plates in a controlled temperature, automated optical density reader (Flx-Xenius XMA, Safas, Monaco). Each experiment was performed in triplicates with independently prepared bacterial cultures. Volumes of 40 ml of bacterial cultures were collected at regular time points to obtain samples at all phases of bacterial growth. Sampling times at cold (warm) temperatures were designated by Cn (Wn) with $0 \le n \le 7$.

## Fatty Acid Analysis

Samples were successively centrifuged at 10,000 × g for 10 and 5 min. Supernatants were discarded, and pellets were washed with 30 ml of saline buffer followed by two centrifugations at 7,000 × g for 5 min. Transesterification of lipids to produce fatty acid methyl esters (FAMEs) and gas chromatography coupled to mass spectrometry analysis were done as previously described (de Sarrau et al., 2012; Ginies et al., 2016). FAME extracts were injected into a gas chromatography coupled to mass spectrometry system (Shimadzu QP 2010) equipped with an UBWax capillary column 30 m × 0.25 mm × 0.5 µm (Interchim, Montluçon, France). Each FAME was identified by its retention time, using the equivalent chain length method and

<footnote n="1">¹https://www.combase.cc/index.php/en/</footnote>
its characteristic ions (Brillard et al., 2010; de Sarrau et al., 2012; Diomande et al., 2015b). FA composition was determined by internal normalization and calculated by the following equation:

$$ \% FAME = \frac{\text{Area FAME}}{\text{Total Area FAMEs}} $$

## Statistical Analysis

A principal component analysis (PCA) was done for each strain with XLSTAT software (Addinsoft) using FAs' relative abundance as quantitative variables and growth steps at each temperature as observations. An analysis of variance (ANOVA) was done for each strain considering growth temperatures (two levels, cold and warm) and growth step (seven levels) as factors, FAs relative abundance as quantitative variables, and their combinations. Comparison of mean abundance in each FA between each growth temperature-growth step combination was performed by Tukey's honest significant difference test at 5% level.

## RESULTS

### Growth of the Three B. cereus Strains at Cold and Warm Temperatures and Selection of Samples Analyzed for Fatty Acids

At 30°C, the three strains ATCC 10876, ATCC 14579T, and MM3 had already increased in number as soon as 0.5 h, without noticeable lag time (Figure 1). For all strains, the first three sampling times for FA analysis at 0.5, 1, and 2 h corresponded to a phase of rapid growth, the fourth sampling time at 4 h corresponded to the end of rapid growth, the fifth to seventh sampling points (6, 7, and 24 h) to stationary phase.

At 12°C, MM3 was the first strain to initiate growth, then ATCC 10876 and ATCC 14579T (Figure 1 for CFU and Supplementary Figure 1A for OD600). At 10°C, the growth of MM3 seemed to start immediately (lag time lower than 1 h, Figure 1). Lag time at 12°C could be estimated to 6 h for ATCC 10876 and to 18 h for ATCC 14579T from CFU data. No estimation could be done for the lag time of MM3 at 12°C. Additional growth experiments with a high density of sampling times in microtiter plates with an automated optical reader confirmed this strain order ranking of increasing lag time at 12°C, from MM3 to ATCC 14579T (Supplementary Figure 1B).

MM3 differed from the two other strains by a decline in numbers during the stationary phase, e.g., after a maximum reached at 24 h at 12°C. At 10°C, MM3 reached a maximum at 48 h, similarly to ATCC 14579T and ATCC 10876 at 12°C. This is consistent with the fact that MM3 has a minimum growth temperature lower by 2°C compared with ATCC 10876 and ATCC 14579T (Carlin et al., 2013). For this reason, in the present work, cold temperatures considered were 10°C for MM3 and 12°C for ATCC 10876 and ATCC 14579T. During growth at 12°C, the FA composition of MM3 changed according to similar patterns to 10°C (Supplementary Figure 2) and will not be presented in the subsequent sections.

## Global Fatty Acids Composition and Changes During Growth

Twenty-five FAs were identified and were the same for the three strains. Their chain length ranged from 12 to 18 carbons and consisted of linear, branched-chain iso, branched-chain anteiso, monosaturated, and di-saturated FAs. Considering all strains, growth temperatures, and times, the 10 most abundant FAs were nine FAs between 15 and 17 carbons long and the branched-chain i13 FA (Table 1). They contained linear, iso branched-chain, monounsaturated fatty in positions 5 or 10 and one di-UFA in positions 5 and 10. The same most abundant FAs (except i16) were also the ones with the highest range of variation observed over the complete experimental design (i.e., all temperatures and times tested).

Changes in the relative abundance of FAs were observed during growth at both cold and warm temperatures, as illustrated in Figure 2 for strain ATCC 10876 and for the two other strains, MM3 and ATCC 14579T, in Supplementary Figure 2. In Figure 2 and Supplementary Figure 2, FAs' relative abundance is presented as an increase or decrease from the abundance in the inoculum. The FAs' composition of the inoculum for the three strains is available in Supplementary Figure 3. For both cold and warm growth experiments, the inoculum was a mid-exponential culture at a warm temperature. Changes observed during growth at warm temperatures are therefore linked to growth phases only, whereas those occurring during growth at cold temperatures are linked to both cold adaptation and growth phases.

### Both Growth Temperature and Growth Phase Modified Fatty Acid Composition of B. cereus

The PCA of ATCC 10876 samples is shown in Figure 3. Samples grown at cold (12°C) and warm (30°C) temperatures were clearly separated along the first axis F1, which accounted for 57.2% of the variance. Time zero samples from warm and cold experiments clustered together with warm samples as both cold and warm growth experiments started from warm grown inocula. The second axis F2 (18.2%) of the variance discriminated samples according to the growth time. Samples at cold early growth separated from samples at late growth (growth times C1–C3 vs. C4–C6, corresponding to, respectively, the 3–24-h and 48–96-h growth periods). At warm temperature, growth time 4 (corresponding to 4 h) was the furthest from growth time 1 along axis F2, whereas the latest growth time (7, corresponding to 24 h) was close to time zero.

The variables (FAs' relative abundances) contributing the most to the first axis in the direction of the cold grown samples were Δ5 mono-UFAs (particularly n16:1Δ5), the Δ5,10 di-UFA, and the BCFAs i14 and a15. In the opposite direction, on the side of warm grown, samples were located Δ10 mono-UFAs and the long BCFAs i17, a17, and i16. The short BCFAs i13 and a13 contributed to both axes on the side of early and cold growth samples, whereas n16, i12, and i15 contributed to both axes on the side of early and warm growth samples (Figure 3).

In PCA of the three strains together, all strain samples were distributed according to growth temperatures and growth
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Growth of the three strains of *B. cereus* ATCC 10876 (green squares), ATCC 14579<sup>T</sup> (blue circles), and MM3 (red triangles) at warm (upper panel) and cold (lower panels) temperatures. Warm temperature was 30°C for all strains. Cold temperatures were 12°C (close symbols) and 10°C for MM3 only (open symbols). Arrows indicate growth times at which cultures were sampled for fatty acid analysis (black: 0, 0.5, 1, 2, 4, 6, 7, and 24 h at 30°C, 0, 3, 7, 24, 48, 72, and 96 h at 10 and 12°C) and for MM3 at 10°C only (red: 5 h). At cold temperatures, growth during the first 7 h is detailed in the lower right-hand panel. Error bars represent standard deviation (n = 3).</figcaption></figure>

periods (Figure 4) as shown for strain ATCC 10876 in Figure 3. ATCC 10876 and ATCC 14579T samples grown at warm temperatures were distributed together, whereas those from MM3 stood clearly apart. At cold and late growth, ATCC 14579T samples were located at more extreme positions of both axes than ATCC 10876.

Few Fatty Acids Changed Similarly
During Growth Among the Three Strains
of *B. cereus*

For each strain, an ANOVA was done for each of the 25 FAs'
relative abundances as quantitative variables, considering the two
factors, growth temperature and growth time.

For the three strains, growth temperature, growth time, and
their interaction had a significant effect (*p* < 0.05) for most FAs,
and for all FAs, at least one factor and/or the interaction had a
significant effect (Supplementary Table 1).

For each FA, Tukey's honest significant difference test at 5% level was used to compare all pairs of mean relative abundance for each growth temperature-growth time combination. This resulted in a classification of eight patterns (named ColdA, ColdB, ColdC, WarmA, WarmB, WarmC, ND, and NC) of changes in each FA abundance with temperature and time (Table 2). Some FAs were significantly more abundant at cold temperatures than at warm temperatures, at least at some stages of growth (patterns ColdA and ColdB in Table 2). FAs from ColdA significantly increased and then decreased, peaking in the middle of the growth experiment, while those from ColdB
<table><caption>**TABLE 1** | Fatty acids identified in *B. cereus* ordered by mean maximal relative abundance (left) and range of variation observed over all the strains (ATCC 10876, ATCC 14579<sup>T</sup>, and MM3), growth temperatures, and times considered in study.</caption><thead><tr><th></th><th>Mean* maximal abundance (%)</th><th>Range of variation**</th></tr></thead><tbody><tr><td>n12</td><td>0.6</td><td>n12 0.5</td></tr><tr><td>n17</td><td>0.6</td><td>i18 0.6</td></tr><tr><td>i18</td><td>0.6</td><td>n17 0.6</td></tr><tr><td>n15</td><td>1.5</td><td>n15 1.5</td></tr><tr><td>i12</td><td>1.8</td><td>i12 1.5</td></tr><tr><td>n18:1</td><td>2.2</td><td>n18:1 2.1</td></tr><tr><td>i16:1 Δ<sup>10</sup></td><td>4.9</td><td>i16:1 Δ<sup>10</sup> 4.7</td></tr><tr><td>a13</td><td>5.8</td><td>a13 4.8</td></tr><tr><td>n14</td><td>6.6</td><td>a15 5.2</td></tr><tr><td>n18</td><td>6.6</td><td>i14 5.2</td></tr><tr><td>a17</td><td>6.9</td><td>n14 5.9</td></tr><tr><td>a17:1 Δ<sup>10</sup></td><td>7.4</td><td>a17 6.0</td></tr><tr><td>i17:1 Δ<sup>10</sup></td><td>7.6</td><td>n18 6.1</td></tr><tr><td>i14</td><td>8.0</td><td>a17:1 Δ<sup>10</sup> 7.3</td></tr><tr><td>a15</td><td>8.7</td><td>i16 7.3</td></tr><tr><td><b>i16:1 Δ<sup>5</sup></b></td><td><b>8.8</b></td><td><b>i17:1 Δ<sup>10</sup> 7.4</b></td></tr><tr><td><b>i16</b></td><td><b>9.1</b></td><td><b>i16:1 Δ<sup>5</sup> 8.1</b></td></tr><tr><td><b>n16:1 Δ<sup>5</sup></b></td><td><b>10.2</b></td><td><b>i13 8.4</b></td></tr><tr><td><b>i13</b></td><td><b>11.4</b></td><td><b>n16:1 Δ<sup>5</sup> 10.1</b></td></tr><tr><td><b>n16:2 Δ<sup>5,10</sup></b></td><td><b>12.6</b></td><td><b>n16 12.5</b></td></tr><tr><td><b>i17</b></td><td><b>13.5</b></td><td><b>n16:2 Δ<sup>5,10</sup> 12.6</b></td></tr><tr><td><b>n16</b></td><td><b>14.0</b></td><td><b>i17 13.3</b></td></tr><tr><td><b>n16:1 Δ<sup>10</sup></b></td><td><b>17.8</b></td><td><b>n16:1 Δ<sup>10</sup> 14.4</b></td></tr><tr><td><b>i17:1 Δ<sup>5</sup></b></td><td><b>18.1</b></td><td><b>i17:1 Δ<sup>5</sup> 15.0</b></td></tr><tr><td><b>i15</b></td><td><b>31.8</b></td><td><b>i15 22.1</b></td></tr></tbody></table>

significantly increased during growth and reached a plateau. ColdC pattern was present only in MM3 and corresponded to FA with a slight but significant decrease at warm, vs. no changes at cold, temperature. Some FAs were associated with growth at warm temperature (WarmA, WarmB, and WarmC in Table 2) because they significantly increased during growth at warm temperature (WarmB and WarmC) and/or significantly decreased during growth at cold temperatures (WarmA and WarmC). Other FAs changed similarly at both warm and cold temperatures (Pattern ND) or showed no significant changes during growth or punctual changes with no clear trend (NC).

Type ColdA corresponded to saturated SCFAs (<15 carbons), whereas ColdB-C corresponded to other SCFAs, mono-Δ⁵ UFAs, and the diΔ⁵,¹⁰ UFA. Most FAs were strain-specific, except i13 for ColdA and n16:1 Δ⁵, n16:2 Δ⁵,¹⁰, and n18:1 for ColdB (Table 2), which had the same pattern for all strains. These FAs were among the most abundant except n18-1 (Table 1). The warm temperature-associated FAs (WarmA and WarmB) were saturated FAs longer than 15 carbons (with n14 in ATCC 10876 being one exception) and Δ¹⁰ UFAs. Only n16 increased during growth at a warm temperature for all strains (categories WarmB).

Two FAs, i16:1 Δ¹⁰ and i17:1 Δ¹⁰ decreased at cold temperatures for all strains but did not have the same patterns of changes at a warm temperature for all strains. The FA i15 decreased during growth at both warm and cold temperatures for the three strains.

Strain MM3 differed from the two other strains by fewer FAs associated with warm temperature (n16:1 Δ¹⁰ and i17 with an NC pattern for MM3, instead of a WarmB or WarmC patterns) and by a smaller number of Δ⁵ mono-UFAs increasing at cold temperature (n16:1 Δ⁵ in ColdB only for MM3 against n16:1 Δ⁵, i16:1 Δ⁵, and i17:1 Δ⁵ for the two other strains).

## Changes in Fatty Acids Occurred at Different Growth Periods

Among FAs with similar patterns of changes among the three strains (Table 2), i13 increased rapidly during growth at cold temperatures, with a maximum relative abundance at growth step C2 (7 h). The FA n16:2 Δ⁵,¹⁰ increased later, at step C3 (24 h), reaching a maximum at growth step C4 (48 h) (Figure 5 for strain ATCC 10876 and Supplementary Figure 4 for other strains). At warm temperature, i16:1 Δ¹⁰ increased at step W4 (4 h) and decreased at the end of growth, whereas n16 increased at the very last step (24 h) (Figure 5 for strain ATCC 10876 and Supplementary Figure 4). It should be stressed that all replicate experiments were run independently, so some differences in the mean composition of the inoculum between cold and warm growth experiments were observed, although not significant, as shown in Figure 5. The same comment applied to Figure 6 in the next section.

## Changes of Fatty Acids Categories During Growth at Cold and Warm Temperatures

The relative abundance of FA categories was calculated, UFAs, FAs unsaturated in Δ⁵ and in Δ¹⁰ positions, FAs with less than 15 carbons (<C15), anteiso and iso branched FAs, and the ratio of Δ⁵-Δ¹⁰ and the ratio of anteiso to iso FAs. The two factors, growth temperatures and growth times, and their interaction had a significant effect (p < 0.05 according to the ANOVA) on all these variables except growth temperature for the “anteiso-to-iso” ratio (Supplementary Table 1).

Figure 6 presents the changes of these FA categories or their ratios for strain ATCC 10876. Total UFAs increased during growth at cold and warm temperatures, although the highest level reached at a cold temperature was higher than at a warm temperature. In addition, UFAs decreased at the last growth time that was analyzed at a warm temperature, i.e., at 24 h, but increased and reached a plateau at cold temperature. The position of the unsaturations strongly marked the differences between cold and warm temperatures. The ratio of Δ⁵-Δ¹⁰ increased sharply during growth at a cold temperature but remained constant at a warm temperature. The proportion of FAs with a short carbon chain, containing less than 15 carbons, increased to a peak at C3 (24 h) at a cold temperature and then decreased to the initial level, whereas it symmetrically decreased and increased at a warm temperature.

<footnote>*Mean of three replicate samples.</footnote>

<footnote>**Difference between mean maximal abundance and mean minimal abundance among tested samples.</footnote>

<footnote>Ten more abundant fatty acids or with the highest range are indicated in bold.</footnote>
<figure><img src="image_5.png" /></figure>

**FIGURE 2** | Changes in relative abundance (%) of fatty acids during growth of *B. cereus* strain ATCC 10876 at cold temperature (12°C, upper panel) and warm temperature (30°C, lower panel), relatively to the inoculum. Bars represent mean of three independent experiments and present for each growth time fatty acid relative abundance minus fatty acid relative abundance of inoculum. Lines over or below bars are standard deviations. For each fatty acid, growth times were, from left to right, 3, 7, 24, 48, 72, and 96 h at 12°C and 0.5, 1, 2, 4, 6, 7, and 24 h at 30°C. Colors ranged from dark blue for first sampling time to red for sixth (brown for the seventh at 30°C) sampling time.

These two features of cold growth occurred at different growth times. SCFAs started to increase as soon as growth time C1 (3 h) and decrease after C3 (24 h), whereas the Δ⁵/Δ¹⁰ ratio remained constant until C2 (7 h), started to increase at C3 (24 h), and reached a plateau at C5 (72 h).

The ratio of anteiso to iso BCFAs increased during growth, to similar levels, at both cold and warm temperatures and then declined at the end of growth at both temperatures. The total relative abundance of BCFAs (anteiso plus iso) tended to decrease during growth at cold temperatures, resulting in a higher relative abundance in samples grown at warm temperatures until a decline at the last growth phase at a warm temperature.

Even if most FA considered individually had a strain-specific pattern of change during growth (**Table 2**), categories of FAs showed similar patterns among the three strains. Results presented in **Figure 6** for ATCC 10876 were similar for ATCC 14579T and MM3, however, with some variations in the extent and kinetics of changes at cold temperatures (**Table 3** and **Supplementary Figure 5**). In particular, at step C0 (inoculum),

the Δ⁵/Δ¹⁰ ratio was around 1 in MM3 strain against 0.3–0.4 in the strains ATCC 10876 and ATCC 14579T.

DISCUSSION

Changes in FA profiles occurred during growth at all temperatures, even at warm temperatures (30°C). As the inoculums were grown at 30°C in all experiments, changes at warm temperatures are reasonably due to the growth phase and not to adaptation to a new temperature. A few previous studies have shown FA changes in bacterial species other than *B. cereus*, during the stationary phase for *Pediococcus* sp. at 28°C (*Annous et al.*, 1999) or during the whole growth for *Bacillus licheniformis* at 37 and 50°C (*Lopes et al.*, 2019a,b). Globally, no significant changes at 30°C were observed during active growth, i.e., before 4 h (step W4) compared with the inoculum (W0). The inoculum was grown until the mid-exponential phase, and cells presumably did not need to change their FAs. Changes were observed at 4 h
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Principal component analysis of *B. cereus* ATCC 10876 samples grown at warm (W) and cold (C) temperatures. Samples were grown for 0, 3, 7, 24, 48, 72, and 96 h at cold (corresponding to samples C0 to C6, blue points) and 0, 0.5, 1, 2, 4, 6, 7, and 24 h at warm temperatures (corresponding to samples W0 to W7, orange points). Inoculum at both cold (CO) and warm (WO) growth experiments were grown at a warm temperature. Distribution of samples (upper panel) and position of variables (lower panel). Variables were relative abundances of fatty acids.</figcaption></figure>

(step W4; **Figures 5, 6** and **Table 3**), at a time corresponding to the end of growth entry in the stationary phase (**Figure 1**). These changes were characterized by increases in UFAs (notably some i16:1 Δ¹⁰), in anteiso/ratio, and by a decrease in SCFAs (**Figures 5, 6**). They presumably reflect membrane modifications linked to the shift from growth to survival during the stationary phase. However, the oldest cells tested (step W7 corresponding to 24 h at 30°C) were characterized by a fall of UFAs, of the anteiso-to-iso ratio, of BCFAs, and an increase of SCFA (**Figures 5, 6**). As a result, FA profiles of cells at step W7 were closer to the inoculum ones (W0), as indicated by the PCA (**Figure 3**). The increase in the relative abundance of n16, a straight chain and saturated FA, at step W7 (**Figure 5**), might result from the

reduced synthesis of UFA and BCFAs in old cells. The increase of UFAs during growth at warm temperatures is surprising, as this is usually associated with cold growth.

UFAs have a lower melting point than their corresponding saturated FAs. An increase in the relative abundance of UFAs is a well-known mechanism to increase membrane fluidity at cold temperatures (Diomande et al., 2015c; de Carvalho and Caramujo, 2018). Previous studies reported higher monounsaturated proportions at cold than at optimal temperatures in *B. cereus* (Kaneda, 1972; Haque and Russell, 2004; de Sarrau et al., 2012; Chazarreta Cifre et al., 2013). In our work, mono-UFAs increased during growth at both cold and optimal temperatures, but their composition was markedly
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Distribution of samples in the principal component analysis of samples from three *B. cereus* strains MM3 (italic red), ATCC 14579T (blue), and ATCC 10876 (bold green). Codes for growth temperatures and growth time are same as in **Figure 3**.</figcaption></figure>

<table><caption>**TABLE 2** | Classification of fatty acids according to their patterns of changes during growth at cold and warm temperatures for three *B. cereus* strains ATCC 10876, ATCC 14579T, and MM3.</caption><thead><tr><th>Patterns of changes</th><th>Changes at cold temperature</th><th>Changes at warm temperature</th><th>ATCC 10867-12°C</th><th>ATCC 14579T-12°C</th><th>MM3-10°C</th></tr></thead><tbody><tr><td>ColdA</td><td>Increase followed by a decrease</td><td>Stable or decrease</td><td>i13, a13, a15, n12</td><td>i12, n12, i13, a13, a15</td><td>i12, i13, n14</td></tr><tr><td>ColdB</td><td>Increase to a plateau or steady increase</td><td>Stable or decrease</td><td>i14, n15, i16:1 Δ<sup>5</sup>, <b>n16:1 Δ<sup>5</sup>, n16:2 Δ<sup>5,10</sup></b>, i17:1 Δ<sup>5</sup>, <b>n18:1</b></td><td>i14, i16:1 Δ<sup>5</sup>, <b>n16:1 Δ<sup>5</sup></b>, <b>n16:2 Δ<sup>5,10</sup></b>, i17:1 Δ<sup>5</sup>, <b>n18:1</b></td><td>a13, i14, a15, <b>n16:1 Δ<sup>5</sup></b>, <b>n16:2 Δ<sup>5,10</sup></b>, <b>n18:1</b></td></tr><tr><td>ColdC</td><td>Stable</td><td>Decrease</td><td></td><td></td><td>i16:1 Δ<sup>5</sup></td></tr><tr><td>WarmA</td><td>Decrease</td><td>Stable</td><td>i16, i17, a17</td><td>i16:1 Δ<sup>10</sup>, n16:1 Δ<sup>10**</sup>, n18, a17:1 Δ<sup>10</sup></td><td>i16, i16:1 Δ<sup>10**</sup>, i17:1 Δ<sup>10</sup></td></tr><tr><td>WarmB</td><td>Stable</td><td>Increase</td><td>n14, <b>n16*</b>, n16:1 Δ<sup>10*</sup>, n17</td><td>i16, <b>n16</b>, i17, i18</td><td><b>n16</b>, a17, a17:1 Δ<sup>10</sup></td></tr><tr><td>WarmC</td><td>Decrease</td><td>Increase</td><td>i16:1 Δ<sup>10</sup>, i17:1 Δ<sup>10</sup>, a17:1 Δ<sup>10</sup>, i18</td><td>i17:1 Δ<sup>10</sup>, a17</td><td></td></tr><tr><td>ND</td><td>Decrease</td><td>Decrease</td><td>i12, i15</td><td>i15</td><td>i15</td></tr><tr><td>NC</td><td>No clear change</td><td>No clear change</td><td>n18</td><td>n14, n15, n17</td><td>n12, n15, n16:1 Δ<sup>10</sup>, i17, i17:1 Δ<sup>5</sup>, n17, i18, n18</td></tr></tbody></table>

Fatty acids with same patterns of changes for three strains are indicated in bold.

*Early and transient increase at cold temperature, relative abundance similar at start and end of cold growth.

**Early and transient increase at cold temperature, relative abundance lower at the end than at the start of cold growth.

***Decrease at the very end of the growth period at warm temperature.

ND: No difference between behavior at cold and warm temperatures. NC: No clear changes at cold and warm temperatures.

different, with a shift from Δ¹⁰ to Δ⁵ during growth at cold temperatures. An increase of Δ⁵ UFAs in *B. cereus* at 25°C compared with 37°C was previously observed (Chazarreta Cifre et al., 2013). However, Δ⁵ UFAs represented a minority of the UFAs at both temperatures, in contrast to our study done at much lower cold temperatures. To our knowledge, the shift between the proportions of unsaturation positions among mono-UFAs as a

mechanism of cold adaptation in bacteria has not been previously reported. The melting point of mono-UFAs with unsaturation in an odd position (such as Δ⁵) is lower than that in an even position (such as Δ¹⁰) (Gunstone and Ismail, 1967; Knothe and Dunn, 2009). For cis-isomers of linear mono-UFAs of n-C18, the melting point difference between odd and even unsaturation is in the range of 10–15°C (Gunstone and Ismail, 1967;
<figure><img src="image_8.png" /></figure>

**FIGURE 5** | Changes of relative abundance at cold (12°C, C0–6, blue bars) and warm (30°C, W0–7, orange bars) temperatures of fatty acids with representative patterns of changes (see **Table 2**): i13 (A), n16:2 Δ⁵,¹⁰ (B), i16:1 Δ¹⁰ (C), and n16 (D). Results presented are from *B. cereus* strain ATCC 10876. Growth times C0–6 and W0–7 correspond to, respectively, 0, 3, 7, 24, 48, 72, and 96 h and 0, 0.5, 1, 2, 4, 6, 7, and 24 h. Results are mean of three independent experiments. In each panel, bars sharing same letter are not significantly different according to Tukey's honest significant difference test at the 5% level.

Knothe and Dunn, 2009). In comparison, the introduction of cis-monounsaturation in the n-C18 SFA decreases the melting point by 50–60°C (Jalal et al., 1982). The UFAs are diverse in *B. cereus*, including n-C16, i16, i17, and a17, and impacts of unsaturation positions on melting points of these FA are not documented yet. However, we assume that balance between Δ¹⁰ and Δ⁵ mono-UFAs allows a fine adjustment of membrane fluidity in relation to growth temperature, with a dominance of Δ¹⁰ unsaturations at warm temperature, less impacting membrane fluidity. A shift from Δ¹⁰ to Δ⁵ during growth at cold temperatures was observed for the three strains studied, stressing its importance for *B. cereus* cold adaptation. In addition, strain MM3, which is psychrotrophic, already had a higher Δ⁵−Δ¹⁰ ratio than the two mesophilic strains ATCC 10876 and ATCC 14579T before growing at cold temperatures, which might contribute to its better cold adaptation. In the same line, a previous study also showed that mesophilic strains of *B. cereus* have a Δ⁵−Δ¹⁰ ratio lower than psychrotrophic ones when grown at optimal temperatures (Diomande et al., 2015b). Taken together, these results indicate the importance of the Δ⁵−Δ¹⁰ ratio in cold adaptation of *B. cereus*.

The increase in the di-UFA n16:2 Δ⁵,¹⁰ at cold temperatures for the three strains is consistent with previous findings on strain ATCC 14579T (Chazarreta Cifre et al., 2013). *B. cereus*

possesses two FA desaturases, DesA and DesB, both produced at cold temperatures and acting, respectively, on Δ⁵ and Δ¹⁰ positions (Chazarreta Cifre et al., 2013). In strain ATCC 14579T, Kristoffersen et al. (2012) found that at 30°C, desB was expressed at a much higher level than desA. At cold temperature, it was shown in the same strain that desA expression increased markedly (Diomande et al., 2016). These previous studies are consistent with our results showing a higher proportion of Δ¹⁰ unsaturation at a warm temperature and a higher proportion of Δ⁵ unsaturation at cold temperatures. In a mutant strain of ATCC 14579T deleted for a sensor-regulator system and impaired for growth at cold temperature, desA was not induced at cold temperature, and Δ⁵ unsaturations did not increase at cold temperature (Diomande et al., 2015a). This further confirms the importance of DesA and Δ⁵ unsaturations in cold adaptation of *B. cereus*. However, the desaturase DesB presumably plays an important role together with DesA to both produce di-UFAs and adjust the Δ⁵−Δ¹⁰ ratio of mono-UFAs in the function of growth temperature. The only Δ¹⁰ FA with an early and transient increase at cold temperatures was n16 Δ¹⁰ (**Table 2** and **Figure 2**). It then decreased concomitantly with the increase of n16:2 Δ⁵,¹⁰, indicating that DesA might desaturate the product of DesB action on n16. Diverse Δ⁵ mono-UFAs were observed over the three strains. However, among these mono-UFAs, only
<figure><img src="image_9.png" /></figure>

**FIGURE 6** | Changes for ATCC 10876 strain of *B. cereus* at cold (12°C, C0–6, blue bars) and warm (30°C, W0–7, orange bars) temperatures of relative abundance of total unsaturated fatty acids (A), of fatty acids with less than 15 carbons (B), of relative abundance of total branched-chain fatty acids (C), of ratio Δ⁵ to Δ¹⁰ monounsaturated fatty acids (D), and of ratio anteiso to iso branched-chain fatty acids (E). Growth times C0–6 and W0–7 correspond to, respectively, 0, 3, 7, 24, 48, 72, and 96 h and 0, 0.5, 1, 2, 4, 6, 7, and 24 h. Results are mean of three independent experiments. In each panel, bars sharing same letter are not significantly different according to Tukey's honest significant difference test at the 5% level.

<table><caption>**TABLE 3** | Growth time at which changes in fatty acid categories were reached for *B. cereus* strains ATCC 10876, ATCC 14579T, and MM3 at cold and warm temperatures.</caption><thead><tr><th></th><th>ATCC 10876-12°C</th><th>ATCC 14579T-12°C</th><th>MM3-10°C</th></tr></thead><tbody><tr><td><strong>Cold temperature</strong></td><td></td><td></td><td></td></tr><tr><td>Maximum for short-chain (&lt; C15) fatty acids</td><td>7 h</td><td>24 h</td><td>24 h</td></tr><tr><td>Maximum for unsaturated fatty acids (UFA)</td><td>48 h</td><td>96 h</td><td>48 h</td></tr><tr><td>Maximum for Δ<sup>5</sup> to Δ<sup>10</sup> ratio</td><td>72 h</td><td>96 h</td><td>96 h</td></tr><tr><td>Maximum for anteiso to iso ratio</td><td>24 h</td><td>48 h</td><td>72 h</td></tr><tr><td><strong>Warm temperature</strong></td><td><strong>ATCC 10876-30°C</strong></td><td><strong>ATCC 14579T-30°C</strong></td><td><strong>MM3-30°C</strong></td></tr><tr><td>Maximum for unsaturated fatty acids (UFA)</td><td>4 h</td><td>ns</td><td>6 h</td></tr><tr><td>Minimum for short-chain (&lt;C15) fatty acids</td><td>4 h</td><td>ns</td><td>4 h</td></tr><tr><td>Δ<sup>5</sup>-Δ<sup>10</sup> ratio</td><td>No change</td><td>No change</td><td>No change</td></tr><tr><td>Maximum for anteiso to iso ratio</td><td>4 h</td><td>4 h</td><td>4 h</td></tr></tbody></table>

ns: Changes not significant.

n16:1 Δ⁵ consistently increased at cold temperatures for the three strains. This suggests that n16 could be the main substrate for *B. cereus* desaturase at cold temperatures, with some variations among strains concerning the other FAs desaturated.

Both n16:2 Δ⁵,¹⁰ and Δ⁵ mono-UFAs started to increase between 7 and 24 h at cold temperatures (steps C2 and C3) for ATCC 10876 (Figures 5, 6) and between 24 and 48 h for ATCC 14579T and MM3 (steps C3–C4 for ATCC 14579T
and C4–C5 for MM3, **Supplementary Figures 4, 5**). This
 corresponded to the exponential phase for the three strains
 (Figure 1). These FAs continued to increase or remained stable
during the stationary phase. This suggests their importance for
the growth and survival of *B. cereus* in cold conditions. The
marked increase of Δ⁵ unsaturation during the exponential phase
is consistent with previous results showing that *desA* was highly
expressed during an exponential phase in ATCC 14579T strain
(Diomande et al., 2016).

Short-chain (<15C) FAs changed in opposite ways during
growth at cold and optimal temperatures. They decreased at
warm temperatures and increased at cold temperatures. At cold,
they started to increase very early, between 3 and 7 h (steps C1-
C2) for ATCC 10876 (Figure 6) and between 7 and 24 h for
ATCC 14579T and MM3 (steps C2–C3 for 14579T and steps C3–
C4 for MM3, Supplementary Figure 5). This corresponded to a
phase where growth was not yet detectable for ATCC 10876 (i.e.,
lag phase) and has only just started (i.e., early growth) for the two
other strains (Figure 1). Then, the SCFAs decreased after 24 h for
ATCC 10876 and MM3 strains (steps C3 for 10876 and C4 for
MM3) and after 48 h for ATCC 14579T (step C4), corresponding
in the three strains to the end of growth. SCFAs presumably
have a role in the early steps of cold growth but not during the
stationary phase. Although very little is known on what happens
during the lag and early phase of bacterial growth (Bertrand,
2019), our work shows that an increase in SCFAs occurred
during this early phase of *B. cereus* cold temperature adaptation.
Among SCFAs, i13 increased particularly early during growth
(Figure 5 and Supplementary Figure 4) and for the three strains
studied. Previous studies showed the role of lipase in the increase
of i13 at cold temperatures, as deletion of this lipase gene
prevented the i13 increase and reduced *B. cereus* ATCC 14579T
growth at cold temperature (Brillard et al., 2010). Whether this
lipase could modify preexisting FAs and contribute to a rapid
reduction in the size of FAs during growth at cold temperatures
is not known yet.

In conclusion, the kinetics view of FAs composition of three
*B. cereus* strains during growth at cold and optimal temperature

revealed a succession of two main changes during cold growth.
First, cold adaptation starts with an increase in SCFAs before
or during growth initiation and until the exponential phase,
followed by an increase in di-unsaturated and in Δ⁵ UFAs
at the expense of Δ¹⁰ UFAs during exponential phase and
stationary phase.

DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in online
repositories. The names of the repository/repositories and
accession number(s) can be found below: Data are deposited in
the open public repository "data INRAE" (https://data.inra.fr/)
and accessible at: https://doi.org/10.15454/KCSN4H.

AUTHOR CONTRIBUTIONS

MF, VB, FC, CD, and CN-T contributed to the conception and
the design of the study. MF, CD, and RB contributed to the
acquisition of the data. MF, CD, RB, CG, and CN-T analyzed
the data. MF, CG, and CN-T performed the statistical analysis.
CN-T provided acquisition and project administration. MF wrote
the first draft of the manuscript. All authors contributed to its
revision, read and approved the submitted version.

FUNDING

This work was funded by the Agence Nationale de la Recherche
(France) under the ANR-15-CE21-0011 Opticold project.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found
online at: https://www.frontiersin.org/articles/10.3389/fmicb.
2021.694757/full#supplementary-material

REFERENCES

Annous, B. A., Kozempel, M. F., and Kurantz, M. J. (1999). Changes in membrane fatty acid composition of *Pediococcus* sp. strain NRRL B-2354 in response to growth conditions and its effect on thermal resistance. *Appl. Environ. Microbiol.* 65, 2857–2862. doi: 10.1128/aem.65.7.2857-2862.1999

ANSES. (2017). “Étude individuelle nationale des consommations alimentaires 3 (INCA 3)," in *Avis de l'ANSES - Rapport D'expertise Collective*, ed. E. Scientifique (Maisons-Alfort: ANSES).

Beckering, C. L., Steil, L., Weber, M. H. W., Volker, U., and Marahiel, M. A. (2002). Genomewide transcriptional analysis of the cold shock response in *Bacillus subtilis*. *J. Bacteriol.* 184, 6395–6402. doi: 10.1128/jb.184.22.6395-6402.2002

Bertrand, R. L. (2019). Lag phase - a dynamic, organized, adaptive, and evolvable period that prepares bacteria for cell division. *J. Bacteriol.* 201:e00697. doi: 10.1128/jb.00697-18

Brillard, J., and Broussolle, V. (2012). "Mechanisms involved in low-temperature adaptation in *Bacillus cereus*," in *Stress Response in Microbiology*, ed. J. M. Requena (Wymondham: Caister Academic Press), 125–145.

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. doi: 10. 1128/aem.02348-09

Broussolle, V., Pandiani, F., Haddad, N., Michaud, C., Carlin, F., Nguyen-the, C., et al. (2010). Insertional mutagenesis reveals genes involved in *Bacillus cereus* ATCC 14579 growth at low temperature. *FEMS Microbiol. Lett.* 306, 177–183. doi: 10.1111/j.1574-6968.2010.01953.x

Carlin, F., Albagnac, C., Rida, A., Guinebretiere, M. H., Couvert, O., and Nguyen-
the, C. (2013). Variation of cardinal growth parameters and growth limits
according to phylogenetic affiliation in the *Bacillus cereus* Group. consequences
for risk assessment. *Food Microbiol.* 33, 69–76. doi: 10.1016/j.fm.2012.
08.014

Chazarreta Cifre, L., Alemany, M., de Mendoza, D., and 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. doi: 10.1128/aem.01761-13

Chung, B. H., Cannon, R. Y., and Smith, R. C. (1976). Influence of growth temperature on glucose-metabolism of a psychrotrophic strain of *Bacillus cereus*. *Appl. Environ. Microbiol.* 31, 39–45. doi: 10.1128/aem.31.1.39-45.1976

de Carvalho, C., and Caramujo, M. J. (2018). The various roles of fatty acids. *Molecules* 23:2583. doi: 10.3390/molecules23102583
de Sarrau, B., Clavel, T., Clerte, C., Carlin, F., Ginies, C., and Nguyen-the, C. (2012). Influence of anaerobiosis and low temperature on *Bacillus cereus* growth, metabolism, and membrane properties. *Appl. Environ. Microbiol.* **78**, 1715–1723. doi: 10.1128/aem.06410-11

Diomande, S. E., Chamot, S., Antolinos, V., Vasai, F., Guinebretiere, M. H., Bornard, I., et al. (2014). The CasK/R two-component system is required for the growth of mesophilic and psychrotolerant *Bacillus cereus* strains at low temperatures. *Appl. Environ. Microbiol.* **80**, 2493–2503. doi: 10.1128/aem.00090-14

Diomande, S. E., Christophe, N. T., Abee, T., Tempelaars, M. H., Broussolle, V., and Brillard, J. (2015a). Involvement of the CasK/R two-component system in optimal unsaturation of the *Bacillus cereus* fatty acids during low-temperature growth. *Int. J. Food Microbiol.* **213**, 110–117. doi: 10.1016/j.ijfoodmicro.2015.04.043

Diomande, S. E., Doublet, B., Vasai, F., Guinebretiere, M. H., Broussolle, V., and Brillard, J. (2016). Expression of the genes encoding the CasK/R two-component system and the DesA desaturase during *Bacillus cereus* cold adaptation. *FEMS Microbiol. Lett.* **363**:fnw174. doi: 10.1093/femsle/fnw174

Diomande, S. E., Guinebretiere, M. H., De Sarrau, B., Nguyen-the, C., Broussolle, V., and Brillard, J. (2015b). Fatty acid profiles and desaturase-encoding genes are different in thermo- and psychrotolerant strains of the *Bacillus cereus* Group. *BMC Res. Notes* **8**:329. doi: 10.1186/s13104-015-1288-4

Diomande S. E., Nguyen-The, C., Guinbretiere, M. H., Broussolle, V., and Brillard, J. (2015c). Role of fatty acids in *Bacillus* environmental adaptation. *Front. Microbiol.* **6**:813. doi: 10.3389/fmicb.2015.00813

EFSA. (2005). Opinion of the scientific panel on biological hazards on *Bacillus cereus* and other *Bacillus spp* in foodstuffs. *EFSA J.* **175**, 1–48. doi: 10.1128/EFSA-BIOHAZ. (2016). Risks for public health related to the presence of *Bacillus cereus* and other *Bacillus spp*. including *Bacillus thuringiensis* in foodstuffs. *EFSA J.* **14**:4524. doi: 10.2903/j.efsa.2016.4524

Français, M., Carlin, F., Broussolle, V., and Nguyen-Thé, C. (2019). *Bacillus cereus* cshA is expressed during the lag phase of growth and serves as a potential marker of early adaptation to low temperature and pH. *Appl. Environ. Microbiol.* **85**:e00486. doi: 10.1128/aem.00486-19

Freese, E., Sass, H., Rutters, H., Schiedjewski, R., and Rullkotter, J. (2008). Variable temperature-related changes in fatty acid composition of bacterial isolates from German Wadden sea sediments representing different bacterial phyla. *Org. Geochem.* **39**, 1427–1438. doi: 10.1016/j.orggeochem.2008.06.005

Ginies, C., Brillard, J., and Nguyen-The, C. (2016). Identification of fatty acids in *Bacillus cereus*. *J. Vis. Exp.* **118**:e54960. doi: 10.3791/54960

Guinebretiere, M. H., Thompson, F. L., 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

Gunstone, F. D., and Ismail, I. A. (1967). Fatty acids part 14 the conversion of the cis octadecenoic acids to their trans isomers. *Chem. Phys. Lipids* **1**, 264–269. doi: 10.1016/0009-3084(67)90033-3

Haque, M. A., and Russell, N. J. (2004). Strains of *Bacillus cereus* vary in the phenotypic adaptation of their membrane lipid composition in response to low water activity, reduced temperature and growth in rice starch. *Microbiology* **150**, 1397–1404. doi: 10.1099/mic.0.26767-0

Jalal, I. M., Zografi, G., Rakshit, A. K., and Gunstone, F. D. (1982). Thermal analysis of fatty acids. *Chem. Phys. Lipids* **31**, 395–404. doi: 10.1016/0009-3084(82)90074-3

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. doi: 10.1016/0005-2760(72)90097-5

Kaneda, T. (1977). Fatty-acids of genus *Bacillus* - example of branched-chain preference. *Bacteriol. Rev.* **41**, 391–418. doi: 10.1128/mmbr.41.2.391-418.1977

Knothe, G., and Dunn, R. O. (2009). A comprehensive evaluation of the melting points of fatty acids and esters determined by differential scanning calorimetry. *J. Am. Oil Chem. Soc.* **86**, 843–856. doi: 10.1007/s11746-009-1423-2

Kristoffersen, S. M., Haase, C., Weil, M. R., Passalacqua, K. D., Niazi, F., Hutchison, S. K., et al. (2012). Global mRNA decay analysis at single nucleotide resolution reveals segmental and positional degradation patterns in a Gram-positive bacterium. *Genome Biol.* **13**:R30. doi: 10.1186/gb-2012-13-4-r30

Lopes, C., Barbosa, J., Maciel, E., da Costa, E., Alves, E., Domingues, P., et al. (2019a). Lipidomic signature of *Bacillus licheniformis* 189 during the different growth phases unravelled by high-resolution liquid chromatography-mass spectrometry. *Arch. Biochem. Biophys.* **663**, 83–94. doi: 10.1016/j.abb.2018.12.024

Lopes, C., Barbosa, J., Maciel, E., da Costa, E., Alves, E., Ricardo, F., et al. (2019b). Decoding the fatty acid profile of *Bacillus licheniformis* 189 and its adaptation to different growth conditions to investigate possible biotechnological applications. *Lipids* **54**, 245–253. doi: 10.1002/lipd.12142

Luu-Thi, H., Khadka, D. B., and Michiels, C. W. (2014). Thermal inactivation parameters of spores from different phylogenetic groups of *Bacillus cereus*. *Int. J. Food Microbiol.* **189**, 183–188. doi: 10.1016/j.ijfoodmicro.2014.07.027

Miller, R. A., Beno, S. M., Kent, D. J., Carroll, L. M., Martin, N. H., Boor, K. J., et al. (2016). *Bacillus wiedmannii* sp nov., a psychrotolerant and cytotoxic *Bacillus cereus* group species isolated from dairy foods and dairy environments. *Int. J. Syst. Evol. Microbiol.* **66**, 4744–4753. doi: 10.1099/ijsem.0.001421

Pandiani, F., Brillard, J., Bornard, I., Michaud, C., Chamot, S., Nguyen-the, C., et al. (2010). Differential involvement of the five RNA helicases in adaptation of *Bacillus cereus* ATCC 14579 to low growth temperatures. *Appl. Environ. Microbiol.* **76**, 6692–6697. doi: 10.1128/aem.00782-10

Pandiani, F., Chamot, S., Brillard, J., Carlin, F., Nguyen-The, C., and Broussolle, V. (2011). Role of the five RNA helicases in the adaptive response of *Bacillus cereus* ATCC 14579 cells to temperature, pH, and oxidative stresses. *Appl. Environ. Microbiol.* **77**, 5604–5609. doi: 10.1128/aem.02974-10

Russel, N. J. (2008). "Membrane components and cold sensing," in *Psychrophiles: From Biodiversity to Biotechnology*, eds R. Margesin, F. Schinner, J. C. Marx, and C. Gerday (Berlin: Springer), 177–190. doi: 10.1007/978-3-540-74335-4_11

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Français, Bott, Dargaignaratz, Giniès, Carlin, Broussolle and Nguyen-Thé. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.