Original Research ARTICLE
New Insights Into the Response of Metabolome of Escherichia coli O157:H7 to Ohmic Heating
- 1Beijing Advanced Innovation Center for Food Nutrition and Human Health, College of Food Science and Nutritional Engineering, China Agricultural University, Beijing, China
- 2Beijing Higher Institution Engineering Research Center of Animal Product, China Agricultural University, Beijing, China
The objective of this study was to investigate the effects of ohmic heating and water bath heating (WB) on the metabolome of Escherichia coli O157:H7 cells at the same inactivation levels. Compared to low voltage long time ohmic heating (5 V/cm, 8.50 min, LVLT) and WB (5.50 min), the high voltage short time ohmic heating (10 V/cm, 1.75 min, HVST) had much shorter heating time. Compared to the samples of control (CT), there were a total of 213 differential metabolites identified, among them, 73, 78, and 62 were presented in HVST, LVLT, and WB samples, revealing a stronger metabolomic response of E. coli cells to HVST and LVLT than WB. KEGG enrichment analysis indicated that the significantly enriched pathways were biosynthesis and metabolism of amino acids (alanine, arginine, aspartate, and glutamate, etc.), followed by aminoacyl-tRNA biosynthesis among the three treatments. This is the first metabolomic study of E. coli cells in response to ohmic heating and presents an important step toward understanding the mechanism of ohmic heating on microbial inactivation, and can serve as a theoretical basis for better application of ohmic heating in food products.
Thermal treatments, as conventional technology in food processing, are used for the pasteurization, sterilization, dehydration, evaporation, and blanching of foods. Generally, heat energy is generated externally and then transferred into the internal of food by conduction or convection in conventional thermal treatment methods. These methods are time consuming due to slow heat transfer through the product, particularly for larger diameter products, and may lead to overcooked surface and quality deterioration (Jaeger et al., 2016; Kanjanapongkul, 2017). Therefore, there is growing interest in alternative thermal treatment methods, which can avoid these shortcomings (Zell et al., 2010). Ohmic heating in particular is just one of these methods, where the heat is generated directly in the food when electric current passes through conductive food (Hradecky et al., 2017). Compared to conventional thermal treatment methods, the heat of ohmic heating is generated from the internal of the food, therefore it can prevent the surface of the solid food or particles from becoming overheated and preserve sensory attributes of food with a shorter heating time (Kanjanapongkul, 2017).
Besides the superior processing characteristics of ohmic heating, the inactivation effect on microorganisms also attracts many researchers’ interest. Up to now, ohmic heating has been widely applied to inactivate vegetative cells and spores in various food, such as Salmonella in buffalo milk (Kumar et al., 2014), Escherichia coli O157:H7, Salmonella Typhimurium, and Listeria monocytogenes in orange and tomato juice (Sagong et al., 2011) and in salsa (Kim and Kang, 2017), Listeria innocua in meat (Zell et al., 2010), Alicyclobacillus acidoterrestris spores in apple juice (Kim et al., 2017), Bacillus licheniformis spores in cloudberry jam (Pereira et al., 2007), and Bacillus cereus spores in doenjang (Ryang et al., 2016). Most studies indicated that ohmic heating had a comparable or even better inactivation effect on microorganisms than conventional thermal treatment methods, and the possibility of non-thermal inactivation effect of ohmic heating was also referred (Pereira et al., 2007; Somavat et al., 2012; Tian X.J. et al., 2017). However, few studies have successfully proved the existence of the specific non-thermal effects and the effects of ohmic heating on intracellular content on the molecular level.
In recent years, transcriptomic, proteomic, and metabolomic collectively referred to as functional genomics techniques, are becoming increasingly important in the life sciences. These techniques can give new insights and a better understanding of the biological function of a cell or organism (Wu and Li, 2013). Metabolites (MW < 1000) are the building blocks of DNA, RNA, proteins, and lipids and play important roles in cell metabolism, signaling, and regulation for all organisms (Mousavi et al., 2016), their composition or contents can directly reflect the phenotype changes of living organisms (Tian J. et al., 2017). The metabolomics can generate the metabolic profiles of living systems at a specified time and specific environmental conditions. Recently, microbial metabolomics has received much attention due to its potential applications in a wide range of research areas, and several studies have shown metabolite changes when microorganisms were exposed to limited nutrient (Jozefczuk and Al, 2010; Lu et al., 2016), heat (Jozefczuk and Al, 2010), cold (Jozefczuk and Al, 2010; Alreshidi et al., 2015), bacteriostatic agent (Mousavi et al., 2016), toxic (Planchon et al., 2017), and oxidative stress (Jozefczuk and Al, 2010).
Untargeted metabolomics have the capacity to implicate previously unexplored biochemical pathways in a particular biological condition, because it can simultaneously detect as many metabolites as possible to maximize the opportunity of identifying compounds (Yanes et al., 2011). The purpose of this study was to investigate the metabolomic response of E. coli O157:H7 cells exposure to ohmic heating using untargeted metabolomic method. Furthermore, conventional water bath heating was performed in order to compare the metabolomic response of E. coli to ohmic heating.
Materials and Methods
Bacterial Strain and Culture Conditions
E. coli O157:H7 (NCTC 12900) was used as experimental strain. Cell cultures were obtained by the same method as our previous study (Tian et al., 2018b). When cultures reached the mid-exponential growth phase (4 h, OD600 = 0.314) they were centrifuged at 5000 × g, 4°C for 10 min, washed twice with phosphate buffered saline (pH 7.2 ± 0.1, electrical conductivity 1.23 ± 0.09 S/m, 0.09 M NaCl, 0.03 M Na2HPO4, and 0.01 M NaH2PO4, PBS), and then re-suspended in the same PBS as above at final concentration approximately 2 × 108 colony-forming units per milliliter (CFU/mL) before treatment.
Water Bath Heating
Water bath heating (WB) was carried out by the same thermostatic water bath equipment, heating cell, and heating method as our previous study (Tian et al., 2018b). In brief, the heating cell was filled with approximately 160 mL E. coli cell suspension and was heated at 80°C water bath accompanied by constant shakes (15 times/min) in order to achieve uniform heating. The same inactivation levels among WB and ohmic heating protocols for metabolomic analysis were designed to reduce less than 2 log CFU/mL of E. coli cell on tryptone soy agar (TSA, pH 7.3 ± 0.2, HB0177, Qingdao Hope Bio-Technology Co., Ltd). Finally, it required 5.50 min for WB to achieve this inactivation level. After heating, the samples were cooled to 4°C in an ice-water mixture immediately, then the samples were pelleted by centrifugation at 5000 × g for 12 min (4°C), quenched using liquid nitrogen immediately, and stored at -80°C until metabolomic analysis. Each treatment was performed in six biological replicates, and samples without any treatment were served as control (CT).
The ohmic heating equipment (frequency 50 Hz) used was the same as our previous study (Dai et al., 2013), and the experiment was carried out according to the method of Tian et al. (2018b). Approximately 160 mL E. coli cell suspension was assigned to the heating cell and was treated by 10 V/cm and 1.75 min (high voltage short time, HVST), and 5 V/cm and 8.50 min (low voltage long time, LVLT), respectively. After heating, samples were treated for subsequent analysis using the same method as WB samples.
Measurement of Inactivated and Injured E. coli Cells
The inactivation of E. coli cells was measured by counting plates according to the methods described in our previous study (Tian et al., 2018b). TSA, thin agar layer (TAL), and selective medium Improved-MacConkey sorbitol agar (IMSA, pH 7.0 ± 0.2, 02–328, BeiJing AoBoXing Bio-Technology Co., Ltd) were used to assess the dead, alive, and sublethal status of E. coli cells after different treatments. Samples were diluted serially in sterile 0.85% physiological saline, and an appropriate dilution of 100 μL was spread on the three plates. The plates were incubated at 37°C for 24–48 h before enumeration, and each dilution was performed in duplicate. The sublethal ratio was calculated according to the following equation (1) (Bi et al., 2015):
Where CFU/mLIMSA represented the colony counts on the IMSA; CFU/mLTAL represented the colony counts on the TAL.
Untargeted Metabolomic Analysis
Sample Extraction and Preparation
The cell precipitate was added 20 μL L-2-chlorphenylalanine solution (0.3 mg/mL, dissolved in methanol) as internal standard, and then the mixture was transferred to the glass vial with pre-cooled methanol/water (4:1, v/v). The mixture was added 200 μL chloroform and resuspended. The suspension was kept on ice and the cells were disrupted by an ultrasonic cell disruptor (500 W, 6 min, 6 s on, 4 s off) (JY98-IIIN, NingBo XinZhi Bio-Technology Co., Ltd). Then metabolites inside the cells were released by an ultrasonic cleaner for 20 min (TYHD-600, Beijing TianYou HengDa Bio-Technology Co., Ltd). After ultrasound, the suspension was centrifuged at 10,000 × g, 4°C for 15 min, and the supernatant was collected. One mL of supernatant was added into 1.5 mL centrifuge tube (twice, 0.5 mL each time) and was volatilized. Dried extracts were dissolved in 250 μL methanol/water (7:3, v/v) followed by vortex 30 s, ultrasonic 2 min. Next, the suspension was centrifuged at 10,000 × g, 4°C for 15 min, and 180 μL of supernatant was transferred to the sample bottle with lining tube for HPLC-MS/MS analysis.
Metabolites of E. coli cells with different treatments were analyzed by instrument platform of ultra-high-performance liquid chromatograph-dual pressure linear well-electrostatic field orbital well tandem mass spectrometer (UHPLC-LTQ Orbitrap, Thermo Fisher Scientific, United States). An acquity BEH C18 column (100 mm × 2.1 mm, i.d., 1.7 μm; Waters, Milford, CT, United States) was used. The mixed mobile phase contained formic acid-aqueous solution (0.1%, v/v, A) and acetonitrile solution containing 0.1% formic acid (v/v, B), and gradient elution was as follows: 5–25% B over 0–1.5 min, 25–100% B over 1.5–10 min, holding at 100% B over 10–13 min, 100–5% B over 13–13.5 min, and holding at 5% B over 13.5–14.5 min. The column was maintained at 45°C. Injection volume was 3.00 μL and flow rate was 0.40 mL/min.
An electrospray ionization (ESI) source in either positive or negative ion mode was used to acquire mass spectra profiles. The electrospray capillary voltage, input voltage, and collision voltage were 3.0 kV, 40 V, and 30 eV, respectively. The capillary and ion source temperature were all set at 350°C, with a carrier gas flow rate of 45 L/h. The acquired mass data was collected from 50 to 1000 m/z with the resolution of 30,000.
Quality control (QC) sample was prepared by mixing all samples equivalently to be a pooled sample, and then analyzed using the same method with the analytic samples. The QC samples were injected at regular intervals (every 10 samples) throughout the analytical run to provide a set of data from which repeatability could be assessed.
The metabolomics processing software progenesis QI (Waters Corporation,Milford, CT, United States) was used for baseline filtering, peak identification, integration, retention time correction, peak alignment, and normalization of the data sets. The retention time, mass ratio, and peak intensity was obtained.
The positive and negative data were combined to get a combine data set which was imported into SIMCA-P+ 14.0 software package (Umetrics, Umeå, Sweden). Firstly, the unsupervised principle component analysis (PCA) was carried out to visualize the overall distribution of the samples and the stability of the whole analysis process. Then, the supervised (orthogonal) partial least-squares-discriminant analysis [(O) PLS-DA] was performed to find the inter-group differential metabolites. Variable importance in projection (VIP) ranked the overall contribution of each variable to the (O) PLS-DA model, and those variables with VIP > 1 were considered relevant for group discrimination. In this study, the default 7-round cross-validation was applied with 1/seventh of the samples being excluded from the mathematical model in each round, in order to guard against overfitting. R2 and Q2 values were used to evaluate the accuracy and predictive ability of the models.
Multidimensional analysis of (O) PLS-DA and single-dimensional analysis (student t-test) were used to screen the inter-group differential metabolites (VIP > 1, P < 0.05). The metabolites were identified using human metabolome database1 and METLIN database2. Then differential metabolites were annotated by KEGG database, including participating pathway and functional classification.
One-way variance (ANOVA) analysis was carried out by SPSS 21 software (IBM, United States), and results were considered to be statistically significant at P < 0.05. All experiments were performed in sextuplicate, and the values were given as means ± standard deviation of six replicates.
Results and Discussion
Inactivation and Sublethal Injury of E. coli Cells
In order to obtain appropriate treatment conditions for the metabolomic analysis, E. coli cells were HVST, LVLT and WB-treated for 1.75, 8.50, and 5.50 min with final temperature of 57.90, 58.60, and 59.50°C, respectively (Table 1). Under these treatment conditions, cell counts ranged from 6.81 to 6.83 log CFU/mL on TSA, ranged from 6.40 to 6.48 log CFU/mL on TAL, and ranged from 5.57 to 5.60 log CFU/mL on IMSA. There was no significant difference of logarithmic reduction on the same plate among HVST, LVLT, and WB-treated samples with similar final temperature (P < 0.05). The initial E. coli populations were 8.24, 8.16, and 8.07 log CFU/mL on TSA, TAL, and IMSA, respectively. There should be no significant difference of colony counts on TSA, TAL, and IMSA if no sublethally injured cell was induced by the three treatments (Chueca et al., 2015; Tian et al., 2018a). However, as shown in Table 1, the counts of E. coli cells on TSA and TAL were higher than that of IMSA after HVST, LVLT, and WB treatments, which indicated the existence of sublethal injury cells after the three treatments. Based on the difference of plate counts on TSA before and after treatments, there were more than 90% of the E. coli cells killed. According to Equation (1), more than 80% cells was sublethally injured after the three treatments (Table 1). At this inactivation degree, there were enough living cells and a large proportion of sublethally injured cells could respond to heat treatments. Because these treatments provided appropriate conditions for metabolomic analysis of E. coli cells.
In this study, the heating time for HVST was shorter than that of LVLT with the same inactivation levels of E. coli cells, which suggested that a higher voltage gradient could cause a comparable inactivation effect of E. coli cells at the similar final temperature with a shorter heating time. Our results were similar to those reported by other researchers. Kim et al. (2018) reported that increasing treatment voltage gradients (9.43–12.14 Vrms/cm) was an effective way to inactivate E. coli O157:H7, S. Typhimurium, and L. monocytogenes by continuous-type pulsed ohmic heating in buffered peptone water and tomato juice at final temperature of 80°C. When 30, 40, and 50 V/cm voltage gradients were used to inactivate Alicyclobacillus acidoterrestris spores in orange juice, Baysal and Icıer (2010) found that the higher voltage gradient had a more effective inactivation effect on A. acidoterrestris spores at final temperature of 70°C. Lee et al. (2012) also suggested that the most effective treatment voltage was 40 V/cm compared to 30 and 35 V/cm for inactivating E. coli, S. Typhimurium, and L. monocytogenes in orange juice and tomato juice at final temperature of 75.5°C.
Multivariate Statistical Analysis
Metabolome-based class separation was presented in the PCA score plot (Figure 1), the HVST, LVLT, and WB-treated samples were separated from the CT samples. (O) PLS-DA (Figure 2) models were developed for comparison of HVST vs. CT, LVLT vs. CT, and WB vs. CT treated samples. The models displayed good descriptive and predictive abilities, expressed as follows: R2(Y) = 0.992, and Q2= 0.937 in HVST vs. CT; R2(Y) = 0.997, and Q2 = 0.931 in LVLT vs. CT; R2(Y) = 0.991, and Q2 = 0.973 in WB vs. CT (data was not shown).
Figure 1. Principal component analysis (PCA) of metabolites from CT, HVST, LVLT, and WB-treated E. coli O157:H7 cells. t was the first principal component, and t was the second principal component.
Figure 2. OPLS-DA modeling of metabolites from E. coli O157:H7 cells. S-plots (A–C) showed the analysis of datasets for HVST vs. CT, LVLT vs. CT, and WB vs. CT treated samples, respectively. t was the first principal component (predicted principal component), and to was the second principal component (orthogonal principal component). (D–F) Permutation testing confirmed the good quality of each model.
Statistical analysis indicated that after the three treatments, a total of 213 metabolites demonstrated the existence of metabolomic differences with VIP > 1 and P < 0.05, as compared with CT samples (Table 2). Of the 213 differentially expressed metabolites, 73, 78, and 62 metabolites belonged to HVST vs. CT, LVLT vs. CT, and WB vs. CT treated samples, respectively. These metabolites were distributed into 6 functional groups: lipid metabolism (20 metabolites), amino acid metabolism (30 metabolites), nucleotide metabolism (32 metabolites), energy metabolism (8 metabolites), carbohydrate metabolism (8 metabolites), and others (8 metabolites), and some metabolites might belong to more than one category. Moreover, there were 46 overlapping compounds among HVST, LVLT, and WB treatments.
Table 2. List of differential metabolites from HVST vs. CT, LVLT vs. CT, and WB vs. CT treated E. coli O157:H7 cells.
Classification Analysis of Differential Metabolites
Lipid-based metabolism is vital to many biochemistry reactions and related to many biological functions, especially essential to the formation of cell membrane. There were 14, 12, and 14 lipid-metabolism-related differential compounds after HVST, LVLT, and WB treatments, respectively, and most of them were up-regulated, expect for sphinganine and lysoPE (0:0/14:1(9Z)) (POS mode). Among them, only arachidic acid, cPA (18:0/0:0), D-Glycerate 3-phosphate, and myristic acid belong to saturated fatty acids, all the others belong to unsaturated fatty acids. Lysophosphatidylcholine (lysoPC), derived from the hydrolysis of phosphatidylcholines (PC) by phospholipase A2, is of great importance to the cell and participates in many physiological functions (Liu et al., 2013). An increase of lysoPC, such as lysoPC [14:1(9Z)] and lysoPC (15:0) indicated a disturbance of phospholipid catabolism in E. coli cells. Lysophosphatidylethanolamine (lysoPE), a constituent of cell membranes, derived from the hydrolysis of phosphatidylethanolamines (PE), which is catalyzed by phospholipase A2 (Tepperman and Soper, 1999). The up-regulation of lysoPE (0:0/14:0), lysoPE [0:0/14:1(9Z)], lysoPE (0:0/15:0), and lysoPE (0:0/16:0) indicated the changes of phospholipid metabolism or cellular damage. Sphinganine involved in the pathway of sphingolipid metabolism, was up-regulated after WB, down-regulated after HVST, but was not affected by LVLT. This indicated that WB promoted sphingolipid metabolism, whereas HVST inhibited its metabolism, which meant that HVST and WB exerted greater damage to sphingolipid metabolism than LVLT. As indispensable components of cell membranes, sphingolipid might be among the first cell component to encounter extracellular stresses (Liu et al., 2013). Fatty acids, especially unsaturated fatty acids, are known to induce decrease in cell respiratory activity, membrane fluidity, and coagulation of cytoplasmic materials, and eventually lead to cell lysis followed by leakage of macromolecules (Mousavi et al., 2016), and play important roles in the environmental stress. The cell membrane should retain its structural integrity as much as possible to antagonize the heat shock stress (Tian J. et al., 2017). Tian J. et al. (2017) reported that the decreased oleic acid content of Saccharomyces cerevisiae might be a self-protection mechanism of ethanol-adapted strains to maintain membrane integrity through decreasing membrane fluidity. In this study, the increased fatty acid content, especially unsaturated fatty acid might mean that the membrane integrity was damaged and the membrane fluidity of E. coli cells increased. The increase of lipid metabolites could mainly attribute to the heat of the treatments, which might result in remodeling the composition and structure of the cell membrane. This result was consistent with previous study that microorganisms could manifest increased resistance to environmental stress and control strategies after sublethal injury (Liu et al., 2018).
Amino Acid Metabolism
In general, stress could reduce membrane fluidity and accelerate the synthesis of some proteins (Zhang and Rock, 2008). The metabolomic analysis revealed that amino acid metabolism was strongly induced in E. coli cells by all the three treatments. There were 24, 24, and 20 differential compounds involved in amino acid metabolism after HVST, LVLT, and WB treatments, respectively. Among them, there were 15 amino acids overlapping the three treatments, and showed consistent trend of content change. However, in contrast to lipid-based metabolism, most of the differential amino acids were down-regulated, only indole and ketoleucine were up-regulated after LVLT. The observed increase in levels of indole and ketoleucine could attributed to protein denaturing and inhibition of protein synthesis, which might be caused by a halt in the synthesis of essential enzymes (Mousavi et al., 2016). Accordingly, the decreased levels of most amino acids might indicate a weakening in enzymatic activity related to protein degradation or a strengthening in enzymatic activity related to protein synthesis. For instance, glutamate is the direct ammonia assimilation product by glutamate dehydrogenase with high external ammonia concentrations, which subsequently serves as a primary precursor in multiple biosynthesis pathways, and it is usually synthesized more in actively growing cells (Lu et al., 2016). Glutamate was down-regulated after the three treatments, indicating that cell activity was reduced. The reductions of amino acids might reflect their consumption to produce new essential proteins or to repair damaged or misfolded proteins, in order to facilitate acclimation to the changing environment. The microorganisms are prone to rapidly changing when the environment conditions change, such as shifts in temperature, osmotic pressure, pH, or nutrient availability, and they have developed many strategies to cope with such unfavorable conditions. Among these strategies, acquisition of thermotolerance is mainly controlled by the activation and regulation of heat stress-related genes involved in the synthesis of specific compounds that protect the microorganism from thermal stress, which involves in induction of several proteins including stress proteins and chaperones (Paul et al., 2012; Dong et al., 2017). Therefore, in order to survive from the three treatments, E. coli cells might synthesize more molecular chaperones to regulate metabolism, which resulted in reduction of the most amino acids. Another reason for more down-regulation of amino acids might due to the ATP deficiency in response to heat shock from the three treatments, where synthesis of amino acids was required ATP as the energy source (Li et al., 2015).
As precursors of DNA and RNA, nucleotides participate in cell signaling and regulate many metabolic pathways, and play a vital role in stress response (Liu et al., 2013). Most purines and pyrimidines are present in the cell as nucleotides, and they are involved in the biosynthesis of genetic information carriers (DNA and RNA) or suppliers of energy (ATP and GTP) (Hu et al., 2017). The changes in purine and pyrimidine metabolism can suggest the increase of DNA damage and cell turnover (Zhou et al., 2017), where the nucleotide biosynthesis are direct indicators of DNA replication, cell division, and growth status, revealing a pronounced effect on cell proliferation (Bhat et al., 2015). The up-regulation and down-regulation changes of metabolites involved in nucleotide metabolism did not behave like lipid metabolism or amino acid metabolism, although more than half of the metabolites showed significant increase after the three treatments. Adenine was up-regulated after the three treatments, but its derivatives (adenosine, cAMP, ADP, and ATP) changed differently among the three treatments. Adenosine was down-regulated after HVST and WB, and was up-regulated after LVLT; cAMP was up-regulated after HVST and WB, and was down-regulated after LVLT; as one of the five kind nucleotides of DNA or RNA synthesis, the change of cAMP meant that the DNA or RNA synthesis was suppressed by HVST and WB, but was promoted by LVLT. In the meantime, ADP and ATP were only up-regulated after LVLT, the reason might be that the more energy was required in DNA synthesis during LVLT. Cytosine and its derivatives (5′-CMP and CTP) were all up-regulated after LVLT, this change was similar to adenines. Hypoxanthine was down-regulated after the three treatments, and the derivative inosine (IMP) could be converted to AMP and GMP, its up-regulation could promote generation of AMP and GMP during LVLT. In this study, guanine and its derivatives (cGMP and dGMP) were up-regulated after the three treatments, which meant that the IMP was mainly converted to AMP, and this result provided evidence that DNA or RNA synthesis were promoted. Uracil was down-regulated after HVST and WB, but was up-regulated after LVLT, which indicated that RNA synthesis was disturbed by the three treatments. However, UMP was down-regulated and UDP was up-regulated after LVLT, which meant that mRNA synthesis was promoted by LVLT; this result was consistent with the down-regulation of most amino acids, where some proteins were synthesized to resist stress. Previous study also proved that genes of Streptococcus agalactiae involved in purine metabolism were significantly up-regulated at 40°C than 30°C in the study of S. agalactiae transcriptomic analysis (Mereghetti et al., 2008).
Several changes were observed in the levels of the metabolites that were involved in carbohydrate metabolism. Specifically, significant changes of metabolites (citric acid, malic acid, succinic acid, and lactic acid) related to tricarboxylic acid cycle (TCA cycle) were observed. Citric acid was up-regulated after HVST and LVLT, which was synthesized from oxaloacetic acid; but malic acid, as the precursor of oxaloacetic acid, was down-regulated after HVST and WB, and up-regulated after LVLT. Succinic acid was down-regulated after the three treatments, but lactic acid was only down-regulated after WB. The changes of metabolites involved in TCA cycle might be due to the energy requirement during synthesis of proteins and nucleotides. Additionally, citric acid is reported as a powerful chelator and it may play a role in managing concentrations of cations such as Ca2+ for survival, and there is evidence suggesting that Ca2+ is involved in the regulation of cell division and gene expression in response to external stimulation in prokaryotes (Alreshidi et al., 2015). Citric acid was up-regulated after HVST and LVLT, and the reason might be that the combination of citric acid and cations was damaged by the electric current during HVST and LVLT. α-ketoisovaleric acid, a branched-chain organic acid, served as a precursor in leucine and valine synthesis, the up-regulation after LVLT might be required by L-Valine synthesis (Li et al., 2017).
Energy is required in basic metabolism, which includes syntheses of proteins, DNA, and RNA (Li et al., 2015). There were 8 metabolites involved in energy metabolism, 7 of them were presented after LVLT, but only 3 and 4 metabolites were presented after HVST and WB. The up-regulation of energy storage compounds (ATP and ADP) during LVLT indicated that they might be used to offset the negative effects of heat or the prolonged electric current stimulation from LVLT, thereby maintaining basic cellular reaction rates (Lu et al., 2016). Similar study reported that temperature variation perturbed the metabolic status of S. agalactiae, including energy metabolism processes, synthesis of proteins, contents of nucleotides, selective utilization of carbon sources, and some cellular materials (Hu et al., 2017). Another study on the metabolomic response of E. coli exposed to titanium dioxide nanoparticles also indicated that metabolites related to energy and growth were up-regulated (Planchon et al., 2017). The higher metabolomic changes of E. coli involved in energy metabolism might also be responsible for the increased heat tolerance of cells, also be partly responsible for the electric current from LVLT.
Enrichment Analysis of the Differential Metabolites
The KEGG pathway enrichment analysis was performed by Fisher’s exact test, and those with P < 0.05 were considered significant pathways. This analysis could provide some additional clues about the complex identified metabolites. As shown in Figure 3, the statistical data revealed that differentially expressed metabolites were enriched to 28, 29, and 25 pathways from HVST vs. CT, LVLT vs. CT, and WB vs. CT treated E. coli cells, respectively. The top 5 pathways of enrichment ratio were aminoacyl-tRNA biosynthesis, 2-oxocarboxylic acid metabolism, alanine, aspartate and glutamate metabolism, arginine biosynthesis, beta-Alanine metabolism response to HVST; alanine, aspartate and glutamate metabolism, aminoacyl-tRNA biosynthesis, arginine biosynthesis, biosynthesis of amino acids, valine, leucine, and isoleucine biosynthesis response to LVLT; and aminoacyl-tRNA biosynthesis, biofilm formation, arginine biosynthesis, alanine, aspartate, and glutamate metabolism, and pantothenate and CoA biosynthesis response to WB, respectively. These results suggested that the most significantly changed metabolites mainly affect biosynthesis and metabolism of amino acid (alanine, arginine, aspartate, and glutamate, etc.) followed by aminoacyl-tRNA biosynthesis among the three treatments.
Figure 3. KEGG pathway enrichment analysis of differentially expressed metabolites from E. coli O157:H7 cells. (A–C) Showed the analysis for HVST vs. CT, LVLT vs. CT, and WB vs. CT treated samples, respectively. ∗ represented P < 0.05, ∗∗ represented P < 0.01, and ∗∗∗ represented P < 0.001.
In summary, in order to obtain similar inactivation levels of E. coli cells by HVST, LVLT, and WB, the required time for HVST (1.75 min) was shorter than LVLT (8.50 min) and WB (5.50 min). The major functional group of metabolites that displayed up-regulation after the three treatments were metabolites involved in lipid metabolism, while a down-regulation was metabolites involved in amino acid metabolism. On the whole, a stronger metabolomic response was caused by HVST and LVLT compared with WB, indicating that electric current might target partial metabolites during ohmic heating. This study provided a detailed description of overall metabolic responses of E. coli cells to ohmic heating, which would facilitate the understandings of ohmic heating on microbial inactivation on the molecular level. In addition, the results described here could provide a theoretical basis for ohmic heating on microbial inactivation in food products, and further facilitate the application of ohmic heating in food industry.
RD was the fund manager of the grants received from the National Key R&D Program of China (2016YFD040040302) and National Natural Science Foundation of China (No. 31271894), and directed and supervised the whole experimental and writing process. XT performed the research plan, experimental process, data analysis, and manuscript writing. QY participated in part of the experimental process and conducted part of the data analysis. DY and FJ participated in part of the experimental process. LS participated in part of the data analysis. XL participated in part of the research plan and provided valuable advice. ZL helped in the metabolomic analysis of microbes. TH helped in the draft revision of the manuscript.
This work was supported by grants from the National Key R&D Program of China (2016YFD040040302) and grants from the National Natural Science Foundation of China (No. 31271894).
Conflict of Interest Statement
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.
Alreshidi, M. M., Dunstan, R. H., Macdonald, M. M., Smith, N. D., Gottfries, J., and Roberts, T. K. (2015). Metabolomic and proteomic responses of Staphylococcus aureus to prolonged cold stress. J. Proteomics 121, 44–55. doi: 10.1016/j.jprot.2015.03.010
Baysal, A. H., and Icıer, F. (2010). Inactivation kinetics of Alicyclobacillus acidoterrestris spores in orange juice by ohmic heating: effects of voltage gradient and temperature on inactivation. J. Food Protect. 73, 299–304. doi: 10.4315/0362-028X-73.2.299
Bhat, S. V., Booth, S. C., Vantomme, E. A., Afroj, S., Yost, C. K., and Dahms, T. E. (2015). Oxidative stress and metabolic perturbations in Escherichia coli exposed to sublethal levels of 2,4-dichlorophenoxyacetic acid. Chemosphere 135, 453–461. doi: 10.1016/j.chemosphere.2014.12.035
Bi, X. F., Wang, Y. T., Zhao, F., Sun, Z. J., Hu, X. S., and Liao, X. J. (2015). Sublethal injury and recovery of Escherichia coli O157:H7 by high pressure carbon dioxide. Food Control 50, 705–713. doi: 10.1016/j.foodcont.2014.10.014
Chueca, B., Pagán, R., and García-Gonzalo, D. (2015). Transcriptomic analysis of Escherichia coli MG1655 cells exposed to pulsed electric fields. Innov. Food Sci. Emerg. Technol. 29, 78–86. doi: 10.1016/j.ifset.2014.09.003
Dai, Y., Miao, J., Yuan, S. Z., Liu, Y., Li, X. M., and Dai, R. T. (2013). Colour and sarcoplasmic protein evaluation of pork following water bath and ohmic cooking. Meat Sci. 93, 898–905. doi: 10.1016/j.meatsci.2012.11.044
Dong, Z., Chen, Z., Wang, H., Tian, K., Jin, P., Liu, X., et al. (2017). Tandem mass tag-based quantitative proteomics analyses reveal the response of Bacillus licheniformis to high growth temperatures. Ann. Microbiol. 67, 501–510. doi: 10.1007/s13213-017-1279-x
Hradecky, J., Kludska, E., Belkova, B., Wagner, M., and Hajslova, J. (2017). Ohmic heating: a promising technology to reduce furan formation in sterilized vegetable and vegetable/meat baby foods. Innov. Food Sci. Emerg. Technol. 43, 1–6. doi: 10.1016/j.ifset.2017.07.018
Hu, W. T., Guo, W. L., Meng, A. Y., Sun, Y., Wang, S. F., Xie, Z. Y., et al. (2017). A metabolomic investigation into the effects of temperature on Streptococcus agalactiae from Nile tilapia (Oreochromis niloticus) based on UPLC–MS/MS. Vet. Microbiol. 210, 174–182. doi: 10.1016/j.vetmic.2017.09.012
Jaeger, H., Roth, A., Toepfl, S., Holzhauser, T., Engel, K. H., and Knorr, D. (2016). Opinion on the use of ohmic heating for the treatment of foods. Trends Food Sci. Technol. 55, 84–97. doi: 10.1016/j.tifs.2016.07.007
Kim, N. H., Ryang, J. H., Lee, B. S., Kim, C. T., and Rhee, M. S. (2017). Continuous ohmic heating of commercially processed apple juice using five sequential electric fields results in rapid inactivation of Alicyclobacillus acidoterrestris spores. Int. J. Food Microbiol. 246, 80–84. doi: 10.1016/j.ijfoodmicro.2017.01.002
Kim, S. S., and Kang, D. H. (2017). Combination treatment of ohmic heating with various essential oil components for inactivation of food-borne pathogens in buffered peptone water and salsa. Food Control 80, 29–36. doi: 10.1016/j.foodcont.2017.04.001
Kim, S. S., Park, S. H., and Kang, D. H. (2018). Application of continuous-type pulsed ohmic heating system for inactivation of foodborne pathogens in buffered peptone water and tomato juice. LWT 93, 316–322. doi: 10.1016/j.fm.2016.10.021
Lee, S. Y., Sagong, H. G., Ryu, S., and Kang, D. H. (2012). Effect of continuous ohmic heating to inactivate Escherichia coli O157:H7, Salmonella typhimurium and Listeria monocytogenes in orange juice and tomato juice. J. Appl. Microbiol. 112, 723–731. doi: 10.1111/j.1365-2672.2012.05247.x
Li, R., Sakir, H. G., Li, J., Shin, H., Du, G., Chen, J., et al. (2017). Rational molecular engineering of L-amino acid deaminase for production of α-ketoisovaleric acid from L-valine by Escherichia coli. RSC Adv. 7, 6615–6621. doi: 10.1039/C6RA26972A
Li, S. S., Hu, X., Zhao, H., Li, Y. X., Zhang, L., Gong, L. J., et al. (2015). Quantitative analysis of cellular proteome alterations of Pseudomonas putida to naphthalene-induced stress. Biotechnol. Lett. 37, 1645–1654. doi: 10.1007/s10529-015-1828-y
Liu, Q., Wu, J., Zhi, Y. L., Lai, S., Lee, N., and Yang, H. (2018). Metabolite profiling of Listeria innocua for unravelling the inactivation mechanism of electrolysed water by nuclear magnetic resonance spectroscopy. Int. J. Food Microbiol. 271, 24–32. doi: 10.1016/j.ijfoodmicro.2018.02.014
Liu, Y., Lin, Z. B., Tan, G. G., Chu, Z. Y., Lou, Z. Y., Zhang, J. P., et al. (2013). Metabonomic studies on potential plasma biomarkers in rats exposed to ionizing radiation and the protective effects of hong shan capsule. Metabolomics 9, 1082–1095. doi: 10.1007/s11306-013-0529-6
Lu, H., Ulanov, A. V., Nobu, M., and Liu, W. T. (2016). Global metabolomic responses of Nitrosomonas europaea 19718 to cold stress and altered ammonia feeding patterns. Appl. Microbiol. Biotechnol. 100, 1843–1852. doi: 10.1007/s00253-015-7095-y
Mereghetti, L., Sitkiewicz, I., Green, N. M., and Musser, J. M. (2008). Remodeling of the Streptococcus agalactiae transcriptome in response to growth temperature. PLoS One 3:e2785. doi: 10.1371/journal.pone.0002785
Mousavi, F., Bojko, B., Bessonneau, V., and Pawliszyn, J. (2016). Cinnamaldehyde characterization as an antibacterial agent toward E. coli metabolic profile using 96-blade solid-phase microextraction coupled to liquid chromatography–mass spectrometry. J. Proteome Res. 15, 963–975. doi: 10.1021/acs.jproteome.5b00992
Paul, S. C., Jain, P., Mitra, J., Dutta, S., Bhattacharya, P., Bal, B., et al. (2012). Induction of Cr(VI) reduction activity in an Anoxybacillus strain under heat stress: a biochemical and proteomic study. FEMS Microbiol. Lett. 331, 70–80. doi: 10.1111/j.1574-6968.2012.02555.x
Pereira, R., Martins, J., Mateus, C., Teixeira, J. A., and Vicente, A. A. (2007). Death kinetics of Escherichia coli in goat milk and Bacillus licheniformis in cloudberry jam treated by ohmic heating. Chem. Pap. 61, 121–126. doi: 10.2478/s11696-007-0008-5
Planchon, M., Léger, T., Spalla, O., Huber, G., and Ferrari, R. (2017). Metabolomic and proteomic investigations of impacts of titanium dioxide nanoparticles on Escherichia coli. PLoS One 12:e0178437. doi: 10.1371/journal.pone.0178437
Ryang, J. H., Kim, N. H., Lee, B. S., Kim, C. T., and Rhee, M. S. (2016). Destruction of Bacillus cereus spores in a thick soy bean paste (doenjang) by continuous ohmic heating with five sequential electrodes. Lett. Appl. Microbiol. 63, 66–73. doi: 10.1111/lam.12588
Sagong, H. G., Park, S. H., Choi, Y. J., Ryu, S., and Kang, D. H. (2011). Inactivation of Escherichia coli O157:H7, Salmonella Typhimurium, and Listeria monocytogenes in orange and tomato juice using ohmic heating. J. Food Protect. 74, 899–904. doi: 10.4315/0362-028X.JFP-10-552
Somavat, R., Hmh, M., Chung, Y. K., Yousef, A. E., and Sastry, S. K. (2012). Accelerated inactivation of Geobacillus stearothermophilus spores by ohmic heating. J. Food Eng. 108, 69–76. doi: 10.1016/j.jfoodeng.2011.07.028
Tian, J., Zhang, S. X., and Li, H. (2017). Changes in intracellular metabolism underlying the adaptation of Saccharomyces cerevisiae strains to ethanol stress. Ann. Microbiol. 67, 195–202. doi: 10.1007/s13213-016-1251-1
Tian, X. J., Wu, W., Yu, Q. Q., Hou, M., Gao, F., Li, X. M., et al. (2017). Bacterial diversity analysis of pork longissimus lumborum following long term ohmic cooking and water bath cooking by amplicon sequencing of 16s rRNA gene. Meat Sci. 123, 97–104. doi: 10.1016/j.meatsci.2016.09.007
Tian, X. J., Yu, Q. Q., Shao, L. L., Li, X. M., and Dai, R. T. (2018a). Sublethal injury and recovery of Escherichia coli O157:H7 after ohmic heating. Food Control 94, 85–92. doi: 10.1016/j.foodcont.2018.06.028
Tian, X. J., Yu, Q. Q., Wu, W., Li, X. M., and Dai, R. T. (2018b). Comparative proteomic analysis of Escherichia coli O157:H7 following ohmic and water bath heating by capillary-HPLC-MS/MS. Int. J. Food Microbiol. 285, 42–49. doi: 10.1016/j.ijfoodmicro.2018.06.005
Wu, Y., and Li, L. (2013). Development of isotope labeling liquid chromatography-mass spectrometry for metabolic profiling of bacterial cells and its application for bacterial differentiation. Anal. Chem. 85, 5755–5763. doi: 10.1021/ac400330z
Zell, M., Lyng, J. G., Cronin, D. A., and Morgan, D. J. (2010). Ohmic cooking of whole beef muscle evaluation of the impact of a novel rapid ohmic cooking method on product quality. Meat Sci. 86, 258–263. doi: 10.1016/j.meatsci.2010.04.007
Keywords: sublethal injury, untargeted metabolomic analysis, HPLC-MS/MS, lipid metabolism, amino acid metabolism
Citation: Tian X, Yu Q, Yao D, Shao L, Liang Z, Jia F, Li X, Hui T and Dai R (2018) New Insights Into the Response of Metabolome of Escherichia coli O157:H7 to Ohmic Heating. Front. Microbiol. 9:2936. doi: 10.3389/fmicb.2018.02936
Received: 17 August 2018; Accepted: 15 November 2018;
Published: 06 December 2018.
Edited by:Learn-Han Lee, Monash University Malaysia, Malaysia
Reviewed by:Ilkin Yucel Sengun, Ege University, Turkey
Brandon Luedtke, University of Nebraska at Kearney, United States
Copyright © 2018 Tian, Yu, Yao, Shao, Liang, Jia, Li, Hui and Dai. 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.
*Correspondence: Ruitong Dai, firstname.lastname@example.org