Abstract

Escherichia coli is commonly found in intestine of human, and any changes in their adaptation or evolution may affect the human body. The relationship between E. coli and food additives is less studied as compared to antibiotics. E. coli within our human gut are consistently interacting with the food additives; thus, it is important to investigate this relationship. In this paper, we observed the evolution of E. coli cultured in different concentration of food additives (sodium chloride, benzoic acid, and monosodium glutamate), singly or in combination, over 70 passages. Adaptability over time was estimated by generation time and cell density at stationary phase. Polymerase chain reaction (PCR)/restriction fragments length polymorphism (RFLP) using 3 primers and restriction endonucleases, each was used to characterize adaptation/evolution at genomic level. The amplification and digestion profiles were tabulated and analyzed by Nei-Li dissimilarity index. Our results demonstrate that E. coli in every treatment had adapted over 465 generations. The types of stress were discovered to be different even though different concentrations of same additives were used. However, RFLP shows a convergence of genetic distances, suggesting the presence of global stress response. In addition, monosodium glutamate may be a nutrient source and support acid resistance in E. coli.

1. Introduction

A change in the genetic code, also known as mutation, is the primary source of genetic variation which gives rise to diversity within a population. When accumulated over generations, these genetic variations may improve the adaptability; hence, the survival of organisms in different environmental conditions [1, 2]. This may in turn induce or preferentially select for further advantageous changes for better adaptation within the environment [35].

Although mutations-conferring advantageous traits have been observed in animals such as lizards [6] and fish [7], it is difficult to study these effects in a laboratory setting due to space and time constraints. For example, it took 36 years for the lizards to show distinct features [6]. On the other hand, bacteria has a number of advantages-fast generation time [2], ability, to fossilize, and resurrection of historical generations [8]. Escherichia coli, a common intestinal bacterium, has been used in a long-term evolutionary experiment spanning more than 2 decades [3, 911]. A number of stress adaptation studies had demonstrated that the growth phases may impact E. coli adaptation. Nair and Finkel [12] suggested that a nonspecific DNA binding protein, dps, may confer multiple stress tolerance at stationary phase, which concur with Jolivet-Gougeon et al. [13]. In addition, the genome of several strains of E. coli had been sequenced, representing a reliable source of genetic knowledge.

In terms of the effects of chemical treatments, bacterial resistance and tolerance to antibiotics are well established and the mechanisms have been widely studied [1417]. In contrast, mechanisms of insusceptibility to nonantibiotic agents, such as food preservatives and antiseptics which might include tolerance or resistance, are less well understood. For example, citric acid inhibits the growth of proteolytic strains of Clostridium botulinum [18], sodium chloride can inhibit the growth of many bacteria such as Listeria monocytogenes [19], Ochrobactrum anthropi [20], and Lactobacillus plantarum [21] by lowering the water activity [22], and fatty acid such as formic [23], propionic [14], and acetic acid [24, 25] are also capable of inhibiting bacterial growth.

As an intestinal bacterium, E. coli comes into contact with the food and chemicals that we consume. A study which treated pigs with ampicillin, a common antibiotic, demonstrated a significant increase in the occurrence of ampicillin-resistant E. coli from 6% to more than 90% after a course of 7 days [26]. It has been suggested that incomplete absorption in the large intestine led to the presence of subtherapeutic doses of antibiotics in the faeces, resulting in evolutionary pressure for intestinal bacteria such as E. coli towards antibiotic resistance [27, 28].

Although the interactions between antibiotics and bacteria have been well studied [16, 26, 29, 30], the interactions between food additives and bacteria remain elusive. These food additives may have an impact on the evolution of intestinal flora. It had been demonstrated that benzoic acid [31] and sodium chloride [32, 33] can affect E. coli physiology. However, the effect at the genome level is not known.

This paper has two aims. Firstly, we aim to examine the adaptability of E. coli ATCC 8739 (a sequenced strain) in a long-term culture environment in the presence of benzoic acid, sodium chloride, and monosodium glutamate (MSG, a common taste enhancer in Asian cooking), singly and in combination. Two concentrations of each additive were used to evaluate the effects of concentration in the adaptability of E. coli. Generation time across passages is used as an estimation of adaptation where decreased generation time across passages in an additive demonstrated that the cells are growing faster compared to an earlier passage. Thus, a decrease in generation time across passages indicates that the cells are adapting to the additive, whereas an increase in generation time suggests stress. The rate of decrease is an estimator of the rate of adaptation. Secondly, we aim to estimate the genomic effects of these adaptations using Nei and Li distance to estimate the genetic distances between the samples after polymerase chain reaction (PCR)/restriction fragments length polymorphism (RFLP).

We hypothesized that, E. coli ATCC 8739 is able to adapt to the food additives; thereby, demonstrating decrease in generation time across passages. Generation time analysis demonstrated E. coli ATCC 8739 is able to adapt to the additives over extended culture, and DNA fingerprinting suggests that benzoic acid, sodium chloride, and monosodium glutamate are exerting evolutionary pressure on the bacterium.

2. Materials and Methods

2.1. Main Culture Experiment

Lysophilised Escherichia coli ATCC 8739 strain (Reference Passage 4 from ATCC) was revived on nutrient agar plate and incubated at 37°C before inoculating into 8 different treatments supplementation in 10 mL Nutrient Broth. These 8 treatments consist of 4 sets of additives with 2 different concentrations each as shown in Table 1.

The subculturing and generation time estimation protocol were as follows (Figure 1). Subculturing was performed by transferring 1% (100 μL) of the previous culture on every Monday, Wednesday, and Friday to the next passage in order for adaptation to occur at the stationary phase of growth. Optical density (OD) readings were taken before the next subculture at 600 nm wavelength to estimate the number of generations within the current passage and to also determine the number of cells that are being inoculated into the new passage. Generation time was measured on every 3rd passage. Glycerol stocks for each treatment were made from 1% of the culture for every 12th passage after culturing on MacConkey agar.

2.2. Swap Experiment

The swap experiment was done fortnightly (6-7 passages interval), involving the transfer of Escherichia coli cells cultured in different treatments to other treatments for the measurement of generation time. Four types of swaps were carried out, whereby the cells were inoculated into the new treatment in a 100 times dilution. The first set of swap involves the inoculation of basal medium- (L SALT) treated cells into the six nonsalt treatments. An example would be inoculating cells grown in L SALT into H MSG treatment. For the second set, cells cultured in high and low concentrations of each treatment were swapped for all treatments. For example, cells growing in H MSG were inoculated into the L MSG media and vice versa. In the third set, cells of high concentration treatments (H MSG, H BA, and H SALT) were each inoculated into the H COMB treatment. The last set is similar to the previous set except that cells of the low concentration were swapped. Cells from low concentration treatments (L MSG, L BA, and L SALT) were each inoculated into L COMB media. OD600 readings were recorded down at intervals, and generation times were calculated for each interval.

2.3. Genomic DNA Extraction

Treatment cultures from every 12th passage interval were used for genomic DNA extraction using the phenol-chloroform method of DNA extraction for Gram-negative bacteria [34]. The DNA pellet was air-dried and dissolved to 100 ng/μL in pH 8.0 Tris/HCl buffer and stored at −20°C.

2.4. Polymerase Chain Reaction

Each reaction consisted of 50 μL of mixture prepared using 200 ng of DNA template in 10pmoles of DNTPs, 50 pmoles of primer, 1 unit of Taq polymerase, and 1X standard buffer (with 1.5 mM of MgCl2) provided by the supplier (New England Biolabs, Inc.). Primer 5, CgCgCTggC; Primer 6, gCTggCggC; Primer 7, CAggCggCg were used separately. The PCR reaction was carried out (Hybaid Limited, PCR express) under the cycling condition of initial denaturation at 95°C for 10 minutes; 35 cycles of amplification at 95°C for 1 minute, 27°C for 1 minute, 72°C for 3 minutes, followed by a final extension at 72°C for 10 minutes before gel electrophoresis in 2% (w/v) agarose gel with 1X GelRed. The primers used were generated by a previously described method [35] using the following rules: (1) the same primer must be suitable as forward and reverse primers, (2) each primer must be between 6 to 15 bases, (3) the predicted amplicon size must be between 300 to 3100 bases in order for resolution in 2% (w/v) agarose gel, and (4) each primer should be predicted to yield 3 or 4 amplicons.

2.5. Restriction Fragments Length Polymorphism

11 μL of PCR product was digested with 1 unit of restriction endonuclease (TaqI, HinfI, or MspI), in a reaction mixture consisting of 1X restriction digestion buffer and 100 ng/μL acetylated BSA made to a total volume of 20 μL with distilled water. HinfI and MspI reaction mixtures were incubated at 37°C, while the TaqI reaction mixture was incubated at 65°C. All reaction mixtures were incubated for 16 hours before analysis in 2% (w/v) agarose gel with 1X GelRed.

2.6. Data Analysis

Cell density was calculated from OD600 readings using the correction suggested by Sezonov et al. [36]. Briefly, the cell density is directly proportional to OD600 readings when OD600 reading is below or equal to 0.3, at which the cell density is equivalent to 5 × 107 cells per milliliter. If OD600 reading is above 0.3, the cell density is estimated by the equation of Cell Density = 52137400 × In (OD600 reading) + 118718650. Generation time for all experiments was calculated from difference in cell density at intervals between 120 and 300 minutes after the inoculation of cells into fresh media, and the geometric mean was calculated. Changes in generation time across passages were tested using t-test for regression coefficient [37]. The migration distance of the bands of PCR and RFLP of different treatments within the same passage was tabulated and a Nei-Li dissimilarity index (DI) [38], where the maximum value of 1 is obtained when there are no common bands when comparing between the 2 treatments, while a minimum of 0 will be obtained when the 2 treatments have exactly the same bands [35]. The correlation coefficient (CC) value between DIs across passages statistically tested against the CC value of 0.95 (~1) using the 𝑍 -test for two correlation coefficients [37] where the 𝑃 value of more than 0.05 would indicate that the null hypothesis (CC is equal to 0.95) is not rejected.

3. Results

3.1. Generation Time

Analysis of the generation times showed that all eight treatments over the passages displayed different rates of decreasing generation times as shown in Table 2. The steepest decline in generation time occurs in H COMB treatment where approximately 2.02 minutes reduction in generation time per passage over 70 passages was observed, followed by L MSG (1.87 minutes), L BA (1.39 minutes), L SALT (1.24 minutes), L COMB (1.22 minutes), H BA (1.15 minutes), H SALT (1.12 minutes), and finally H MSG (0.906 minutes). The regression intercept may be used to estimate the generation time of the cells in each treatment media for the first passage which is indicative of the level of initial stress on the cells. On this basis, the treatment exerting the highest level of stress on the cells would be H COMB, followed by L MSG, H BA, L BA, L COMB, L SALT, H SALT, and H MSG.

3.2. Swap Experiment

The linear regression of the generation time across passages demonstrated that the gradients of the equations are not equal to zero which indicates that the generation times are not constant for the six swaps. Although there are changes in the general trend of generation time across the passages, the 𝑃 values calculated for the six swaps were more than 0.05 which is not significant: L SALT cells to H MSG media, 0.475509; L SALT to L MSG media is 0.421721; L SALT cells to H BA media is 0.250415, L SALT cells to L BA media is 0.4660235; L SALT cells to H COMB media is 0.484887; L SALT cells to L COMB is 0.443381.

The generation times trend of the four swaps (MSG, BA, Salt, and Combination) from low-treatment to high-treatment over 12 swaps change in a decreasing manner (Table 3). With reference to the regression equations, the linear regression gradient of low treatment cells into high treatment media for combination treatment is the steepest followed by that for BA, MSG, and salt treatments. At swap count between eight and nine, the generation time is almost the same for MSG, salt, and combination, but the generation time for BA is still distantly higher.

The generation time trends of the four swaps (MSG, BA, salt, and combination), from high treatment to low treatment over 12 swaps, changed in a decreasing manner. With reference to the regression equations, the linear regression gradient of high-treatment cells into low-treatment media for MSG treatment was the steepest followed by that for BA, combination, and salt treatments. The generation time of high-treatment cells into low-treatment media for MSG treatment was almost the same as for salt treatment between swap count four and five, and as for combination treatment between swap count eight and nine. The generation time for BA is consistently-lowest among the treatments.

The generation time of swapping high-concentration treated cells into high-combination medium showed similar trends. At the 8th swap, H BA suddenly increased in generation time, due to the retarding growth. The OD600 reading ranges from 0.069 to 0.084 over a period of three days. After the 10th swap, generation time for all treatments remain, constant at 200 minutes till the 12th swap.

The generation time of swapping low-concentration treated cells in to low-combination medium showed similar trends. At the 2nd swap, L Salt suddenly increased in generation time, due to the slow growth. The OD600 reading ranges from 0.077 to 0.545 over a period of 3 days. After the 4th swap, generation time for all treatment followed a similar trend till the 12th swap.

3.3. PCR/RFLP

Electrophoresis agarose gels of the PCR and RFLP products for the eight treatments were used to study the differences between the genome of the E. coli cells of the treatments across the passages. Nei-Li dissimilarity index (DI), which had been shown to be suitable for RFLP [39], was utilised to mathematically calculate the dissimilarity between pair-wise comparisons of the treatments.

The dissimilarity index of the 28 comparisons showed a trend of convergence from PCR/RFLP number 4 onwards (Figure 2).

This trend is further elaborated with the estimation of the maximum and minimum mean values (Figure 3) which shows converging linear regression line across the 6 PCR/RFLP.

Six resulting effects (Table 4) obtained from the comparisons were analysed. The similarity among the six resulting effects is that each type of effects had two originating comparisons. Therefore by plotting the two comparisons against each other and testing for significance, we can deduce whether the genomic differences in each of the two comparisons are actually a consequent effect from the resulting effects.

All resulting effects were not statistically significant except for 10BA + SALT. This suggests that the PCR/RFLP-inferred genetic distance between H MSG and H COMB, and H BA and H SALT varied independently (not correlated).

4. Discussion

In this paper, we present one of the first comprehensive investigations of the effect of E. coli cells’ adaptations to a variety of food additives using a long-term culture approach. Our results suggest that cells grown under different stress condition are able to adapt to the environment which can be observed by decreased generation time and genetic variations.

4.1. Nutrient Broth Does Not Prime Cells for Growth in Other Treatments

Since E. coli cells were grown in NB with the various supplementations of treatments, it was important that any changes to the cells were a direct result of the treatments rather than from the NB. The generation time trend for E. coli cells from L SALT (NB) inoculated into six different media (Table 3) showed that nutrient broth did not appear to impact on adaptability as none of the regression gradients were statistically different from a gradient of zero suggesting that the general generation time trend remained almost the same; therefore, nutrient broth (L SALT media) was unlikely the cause of any adaptations observed.

4.2. Cells Adapt to Their Individual Treatments

The low-concentration treated cells were observed to be adapting to their environment as seen from the decreasing generation time. The cells are dividing at a faster rate with increasing passages suggesting lowered stress level in later passages comparing to early passages (Table 3). However, as the concentration of additives is the same throughout the passages, decreased generation time suggests that the cells are adapting to the stress. Low-treatment cells inoculated into high-treatment media also showed decrease in generation time across the swaps. Growth rate has been used as a measure for adaptation to a stressed environment in previous studies where Chen and Shakhnovich [40] had demonstrated increase in growth rate of 35 bacterial species upon adaptation to thermal stress. In addition, Zhu and Yang [41] also demonstrated increase in growth rate; thus, decrease in generation time, when Clostridium tyrobutyricum adapts to the presence of butyric acid. Our results suggest that the low-treatment cells had gradually adapted to its own individual treatment before the swap, causing it to be less stressed when swapped into high-treatment media.

The generation times of cells from low single treatments to L COMB were observed to increase gradually across passages (Table 3). This suggests that low MSG, BA, or Salt treatments were not stressful enough to induce significant adaptations such that when they were placed into the L COMB treatment which now contains additional stress inducers, the cells could not cope. This may suggest that the cells may be gradually optimized to grow in a specific treatment; thus, increasingly specialized to their specific environment. Similar cases had been reported in other evolutionary studies using E. coli [4244]. The effects of L COMB did not increase the adaptability but instead decreases it as seen from the generation time analysis (Table 2). This suggests that L COMB may be less stressful compared to the individual treatments. It has been suggested that the presence of MSG counteracts the effects of drop in pH caused by BA [14]. This is achieved by increasing the resistance of E. coli cells against the lowered pH, which will otherwise kill the cells. This suggests that the effects of L MSG and L BA cancel each other out, leaving only L SALT which is further supported by the similarity between the adaptability of L SALT (−1.24 minutes per generation, Table 2) and L COMB (−1.22 minutes per generation, Table 2).

E. coli cultured in H COMB treatment had the greatest decline in the generation time over 70 passages. Since higher stress level may force the cells to adapt quickly in order to survive [45, 46], suggesting that the E. coli cells in H COMB treatment experienced the highest level of stress among the eight treatments (Table 2) in contrast to H MSG which induced the lowest decrease in generation time. This may suggest that the presence of glutamate in MSG may be aiding the growth of cells [47] as glutamate can serve as an additional source of nutrient for the cells. Thus, H MSG may cause the least amount of stress but instead led to better growth resulting in a lower rate of adaptation; hence, the lowest decrease in generation time.

When high-concentration treated cells were swapped into low-treatment media containing the same type of stress, reduced generation time was observed. This suggests that these cells have adapted to its individual treatment. Although both high and low concentrations appear to result in some adaptations as measured by generation time, the rate of adaptation differs. The general decline in generation times of the cells from low-treatment to high-treatment media is steeper than that of the reverse. High-concentration treated cells inoculated into low-concentration treatment media appeared more stressed. This is surprising as low-concentration treated cells inoculated into high-treatment media should experience more stress than high-concentration treated cells inoculated into low-treatment media [48], provided that the type of stress is similar. A possible explanation for this is that the type of stress may differ, even though low and high treatment contained the same type of additives and differing only in concentration. This may be explained by the catabolism of stress-induced molecules [49]. In Takatsu [49], the level of cardiolipin, a salt stress marker of Staphylococcus aureus, took longer to return to basal level upon reculture in basal media after stressed in 5% NaCl than in 10% NaCl. This may suggest that function of stress-induced molecules may be concentration dependent. It may be plausible that high-concentration treated cells may induce stress-induced molecules which may add to the level of metabolic stress at a lower concentration. However, the adaptive nature of different concentrations requires further studies.

The swap from individual high-concentration treatments (H MSG, H BA, and H SALT) to H COMB showed decreased generation time. This may suggest that pretreatment of cells in a stressed environment may condition them to adapt to another stress environment which has been demonstrated by other adaptation studies using E. coli [5053]. However, the swap from individual low-concentration treatments (L MSG, L BA, and L SALT) to L COMB showed the opposite trend with increased generation time. This further corroborates that the adaptive nature of different concentrations may be different; thus, requires further studies.

4.3. Cells from Different Treatments Become Genetically Similar

Our results from PCR/RFLP showed a converging trend in DI indicating that the E. coli from all treatments are getting similar (Figure 3) suggesting that they mutate in a similar manner. This suggests that they may evolve the same type of stress mechanism and DNA repair. However, all the treatments originate from the same bacterial clone, suggesting that the initial stress adaptation may involve mutation as it has been suggested that hypermutation is a feature in initial stress adaptations [54].

It is known that E. coli exposed to stresses would respond to counteract the effects. Tucker et al. [55] has shown that E. coli in nitric oxide (NO) will reduce the NO to nitrous oxide under anaerobic conditions which is harmless to the cell. In another study, E. coli produces methylglyoxal to counteract toxic electrophiles [56].

While the cells from all treatments may have experienced different types and levels of stresses, it is likely that the cells might have adapted and activated similar stress-responsive mechanism by evolving similarly. This interpretation is supported by a number of studies suggesting the presence of global stress response in E. coli [5759]. In addition, Cebrián et al. [50] found that adaptation to pH stress may protect Staphylococcus aureus against oxidative stress by hydrogen peroxide, suggesting that adaptation to particular stress may confer tolerance to other stresses.

In this paper, the similar response and mutations to the E. coli of the treatments could not be determined for the whole genome. Primers 5, 6, and 7 amplified a random sample amounting to 0.37% of the whole genome which is a limitation of this study. However, several studies had demonstrated that PCR-based DNA fingerprinting using a small number of primers is suitable to examine genetic diversity in E. coli [60], Candida dubliniensis [61], and mackerel [62]. Another limitation is that DNA fingerprinting was performed on the entire population and not on isolated colonies. Hence, this paper can only imply on areas of the genome that were amplified, and analysis was approaching genetic similarity at a population scale. This method had been used other studies examining metagenomics in environmental bacterial samples [63, 64] and human myopia [65]. The genes responsible for stress-handling mechanism may also not be present in the amplified regions of the genome.

It is unlikely that the genetic distance of E. coli in the eight treatments reaches zero as a population, suggesting that the declining trend is likely to taper off. In addition, spontaneous mutation may prohibit them from being identical. On the other hand, previous studies in E. coli [44] and Herminiimonas arsenicoxydans [66] demonstrated the presence of ecological specialization. Our results showed that the generation time decreased over passages, suggesting the possibility of ecological specialization. This may indicate the presence of both global stress response and ecological specialization in E. coli. Global stress response allows for adaptation to new stress environments, but extended stress may lead to ecological specialization. Hence, it may be hypothesized in future studies that continued culture may lead to ecological specialization which may be seen as a divergence in the genetic distance.

Statistical analysis of the selected comparisons indicates that the effects of all the treatments were insignificant except 10BA + SALT (Table 4). Statistical tests for MSG and BA effect suggested that different gels provided constant results; thus, suggesting reliability in our study.

Statistical tests suggest that MSG and BA, and MSG and SALT are likely to interact with each other (Table 4). However, 10BA + SALT do not appear to have an interacting effect. This suggests that high-combination media contains 10MSG + 10BA + SALT and the interacting effects of MSG and BA, and MSG and salt. Results from swap analysis indicated that low-salt cells when swapped to high-combination media showed an increase in generation time. This suggests that the high stress in high combination media results in difficulties for the cells from L SALT to grow and caused an increase in doubling time. This might be due to the additional combined stress produced by MSG interacting with BA and SALT. The presence of the additional interacting stress of the combination treatment can also be observed from the analysis of generation time where the stress level of H COMB is much higher than the three individual high-concentration treatments.

However, some interacting effects may not be hindering the growth of the bacteria. BA kills bacteria by lowering the pH of the media, whereas MSG has the effects of pH resistance on the cells. Hence, the presence of MSG may aid the growth of the E. coli living in low-pH environment [14] caused by the presence of BA. On the other hand, combined effects from MSG + S could be harmful to the cells as salt may increase the high sodium content provided by MSG in media. The high-sodium environment actually changes the environment to be even more selective.

5. Conclusion and Future Work

This paper had demonstrated that E. coli is able to adapt to food additives over an extended period of time by observing the decreased in generation time over a period of 70 passages. This may have implications in using sublethal doses of bacteriocidal agents such as disinfectants and preservatives. Our results suggest the likelihood of both global stress responses and ecological specialization. Hence, it may be hypothesized that increasing passages may demonstrate a shift towards ecological specialization.

Acknowledgments

The authors wish to thank Phil Au (CSIRO, Australia) for his constructive review on the initial draft of this paper. This project is sponsored by Singapore Polytechnic (Account no. 11-27801-45-2672).