Next Article in Journal
Optimisation of Stingless Bee Honey Nanoemulsions Using Response Surface Methodology
Next Article in Special Issue
Dynamics of Microbial Inactivation and Acrylamide Production in High-Temperature Heat Treatments
Previous Article in Journal
Pasta as a Source of Minerals in the Diets of Poles; Effect of Culinary Processing of Pasta on the Content of Minerals
Previous Article in Special Issue
The Inclusion of the Food Microstructural Influence in Predictive Microbiology: State-of-the-Art
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction in the Dynamics and Spoilage of Shewanella putrefaciens in Bigeye Tuna (Thunnus obesus) by Gas Sensors Stored at Different Refrigeration Temperatures

1
College of Food Science & Technology, Shanghai Ocean University, Shanghai 201306, China
2
Shanghai Professional Technology Service Platform on Cold Chain Equipment Performance and Energy Saving Evaluation, Shanghai 201306, China
3
National Experimental Teaching Demonstration Center for Food Science and Engineering, Shanghai 201306, China
4
Shanghai Engineering Research Center of Aquatic Product Processing & Preservation, Shanghai 201306, China
*
Author to whom correspondence should be addressed.
Submission received: 27 July 2021 / Revised: 31 August 2021 / Accepted: 7 September 2021 / Published: 9 September 2021
(This article belongs to the Special Issue Dynamic Modelling and Simulation of Food Systems)

Abstract

:
Shewanella putrefaciens have a faster growth rate and strong spoilage potential at low temperatures for aquatic products. This study developed a nondestructive method for predicting the kinetic growth and spoilage of S. putrefaciens in bigeye tuna during cold storage at 4, 7 and 10 °C by electronic nose. According to the responses of electronic nose sensor P30/2, the fitted primary kinetic models (Gompertz and logistic models) and secondary model (square root function model) were able to better simulate the dynamic growth of S. putrefaciens, with high R2 and low RMSE values in the range of 0.96–0.99 and 0.021–0.061, respectively. A partial least squares (PLS) regression model based on both electronic nose sensor response values and electrical conductivity (EC) values predicted spoilage of S. putrefaciens in bigeye tuna more accurately than the PLS model based on sensor signal values only. In addition, SPME/GC-MS analysis suggested that 1-octen-3-ol, 2-nonanone, 2-heptanone, dimethyl disulfide and methylamine, N, N-dimethyl- are the key VOCs of tuna inoculated with S. putrefaciens.

Graphical Abstract

1. Introduction

Bigeye tuna (Thunnus obesus) is a widely distributed and commercially important fish, favored by consumers because of its good taste and abundant nutrition [1]. However, bigeye tuna is an extremely perishable fish because of microbial spoilage and certain biochemical reactions during processing and storage. Its superior nutritional value and delicious meat make it important to preserve bigeye tuna [2]. Some methods have been used for the preservation of tuna, such as gas packaging, cool store, freezing processing, cryopreservation, etc. Refrigeration is currently an effective storage method used to slow down fish deterioration [3]. The main factor contributing to seafood spoilage during the refrigeration process is the activity of microorganisms. Many studies have reported that the specific spoilage organisms in refrigerated seafood were Shewanella spp., Pseudomonas spp., Aeromonas spp., and Acinetobacter spp. [3,4]. S. putrefaciens is the main spoilage microorganism of seafood in low-temperature storage, which has the potential for decomposing proteins and trimethyl-amine-N-oxide (TMAO) into ammonia, trimethylamine (TMA), and H2S, producing a fishy odor [5]. Shewanella and Pseudomonas species isolated from spoiled tuna were considered as potential main contributors to spoilage in tuna during refrigerated storage [6]. Many studies have reported that the growth of S. putrefaciens may cause tuna to deteriorate during refrigeration [7,8]. Therefore, evaluating the spoilage potential of S. putrefaciens is crucial in the spoilage control of tuna at low temperatures.
Spoilage influences shelf life, marketing options, and safety of the product; no one buys spoiled foods nor should spoiled foods be on the market. Researchers usually use many conventional methods to identify fish spoilage, including sensory evaluation techniques, chemical methods [9], and microbiological methods [10]. Although conventional microbiological techniques are economical and simple to perform, these methods are time-consuming and cannot be continuously monitored in real time. The odor is an important indicator to evaluate the freshness of fish. Odors such as amines, ammonia, trimethylamine, and volatile sulfides are produced in marine fish during spoilage [11]. These volatile compounds can be potential indicators of spoilage in marine fish. Gas chromatography/mass spectrometry (GC-MS) has become a standard instrument for quantitative analysis of volatile substances in laboratories [12]. However, it is expensive, time-consuming, and unsuitable for large-scale detection. Moreover, the electronic nose (E-nose), a gas sensor array technology, has become an effective tool for predicting fish spoilage [13]. It can mimic the human olfactory system with sensitive sensors that interact with multiple odors to generate different electrical signals [14]. The E-nose has many advantages in predicting food spoilage, such as portability, non-destructive samples, low cost, short time consumption, and high sensitivity. For example, Semeano et al. [15] developed a system based on a gas sensing gel material coupled with an optical E-nose to detect tilapia deterioration, and the system predicted microbial growth well. In addition to the rapid detection of fish spoilage by the E-nose, fish spoilage can also be predicted using a simple physical index—electrical conductivity (EC). The decomposition of tissues and the outflow of electrolytes during the fish spoilage because of the catabolic activity of microorganisms and the oxidation of the fish body eventually lead to the rise of EC [16]. Heising et al. [17] found that EC values of aqueous solutions of volatile compounds produced by cod were positively correlated with freshness. Other researchers have also found a strong correlation between fish conductivity or electrical impedance and fish spoilage, and EC may be a good predictor of fish spoilage [18,19]. However, few studies have involved the use of electronic noses to predict the dynamic growth of spoilage bacteria in seafood and to predict the spoilage of seafood inoculated with spoilage organisms.
The dynamic growth of S. putrefaciens and the spoilage potential in aquatic products are of great importance. It is necessary to find a quick and easy method to study the spoilage of marine fish contaminated by S. putrefaciens. At present, there are few studies to predict the spoilage as well as the dynamic growth of specific spoilage bacteria in marine fish at different refrigeration temperatures. For this purpose, sterile tuna blocks were inoculated with S. putrefaciens and the changes in the total number of S. putrefaciens (TNS), total volatile basic nitrogen, EC, and volatiles at different refrigeration temperatures were investigated. The sensor P30/2 of the E-nose was selected to simulate the primary and secondary dynamic growth of S. putrefaciens in tuna. A partial least squares (PLS) regression model based on both E-nose sensor response values and EC values was used to predict the spoilage potential of S. putrefaciens in bigeye tuna. Finally, some key volatile organic compounds (VOCs) of tuna inoculated with S. putrefaciens were identified by GC/MS, and the correlation between VOCs and gas sensor signal values was analyzed.

2. Materials and Methods

2.1. Bacterial Strains and Cultural Conditions

S. putrefaciens was isolated and identified from spoiled bigeye tuna (Zhejiang Fenghui Ocean Fishing Company Ltd., Zhoushan, Zhejiang, China) and was identified based on 16S rRNA gene sequences, compared in GenBank using the BLAST function. Spoiled bigeye tuna was evaluated by trained panelists from the College of Food Science and Technology, Shanghai Ocean University, based on odor, color, and appearance. Bacteria were stored in tryptone soya broth (TSB) containing 25% glycerin at −80 °C. Before use, S. putrefaciens was precultured in brain–heart perfusion infusion (BHI) at 30 °C for 18 h and then cultured in TSB until the maximal concentration (108 CFU/mL).

2.2. Sample Preparation and Inoculation

The back muscle blocks of 20 kg tuna were purchased from Zhejiang Fenghui Ocean Fishing Company Ltd., Zhoushan, Zhejiang, China, and divided into rectangular blocks of about 30 g. Three replications were taken for TNS, pH, TVB-N, TMA, and VOC measurements, with ten replications for EC and e-nose measurement. Then, the blocks were sterilized by soaking in 0.5% (v/v) formalin solution for 40 s and washed in sterile water 2 times. Each sterile block was immersed in a bacterial suspension for 30 s of S. putrefaciens inoculation to achieve an inoculum level of 3.0–4.0 log CFU/g. Non-inoculated blocks immersed in sterile normal saline (0.85% NaCl) were used as the control check (CK) group. All samples were packed in a clean tray in a sterile environment and stored at 4, 7, and 10 °C.

2.3. Physicochemical Analysis

The physicochemical analysis included the measurement of pH, EC, total volatile basic nitrogen (TVB-N), and trimethylamine (TMA) values.
The EC of tuna blocks was measured using the method described by Yao et al. [19]. Briefly, tuna flesh (10 g) was homogenized and stirred for 30 min in 100 mL of distilled water. The mixture was filtered, and the EC of the filtrate was measured using an EC meter (Mettler Toledo FE20/EL20, Shanghai, China).
The pH value was determined by the method of [20]. The sample treatment was consistent with EC measurement and the pH of the filtrate was measured using a digital pH meter (Cyberscan Model 510; Eutech Instruments Pvt. Ltd., Singapore).
Total volatile basic nitrogen (TVB-N) was performed according to the method of [21]. Five grams of minced tuna flesh was accurately weighed. The TVB-N value was measured by an Automatic Kjeldahl Apparatus (KjeltecTM8400; FOSS Quality Assurance Co., Ltd., Copenhagen, Denmark).
TMA content was determined by Colorimetric Picric Acid Method [22]. Briefly, fish samples and trichloroacetic acid (TCA) were homogenized and mixed. After centrifugation, the supernatant was mixed with formaldehyde, saturated potassium carbonate solution, and toluene. The toluene layer solution and picric acid were mixed thoroughly, and absorbance readings were taken at 410 nm.
Measurements of pH, EC, TVB-N, and TMA were taken every 2 days for 12 days for samples stored at 4 °C, every 2 days for 10 days for samples stored at 7 °C, and every 1 day for 6 days for samples stored at 10 °C.

2.4. Microbiological Analysis and Growth Curve Fitting

The total number of S. putrefaciens (TNS) was determined by a basic method as described by Qian et al. [23] using iron agar. Briefly, 25 g of tuna flesh were put in 225 mL of sterilized saline water (NaCl, 0.85%, w/v) and homogenized for 2 min. Then, 0.1 mL of the dilution was spread on iron agar (IA) plates after serial dilution and incubated at 30 °C for 48 h. Black colonies were enumerated in IA plates. Plate counting agar was used to count the total viable count (TVC) in the CK group at 30 °C for 48 h. The measurement cycle of TNS and TVC was the same as the above physicochemical indexes.

2.4.1. Primary Models

The primary models, namely, Gompertz and logistic models, were used to simulate the growth curves of S. putrefaciens in tuna. They are represented by the following equations according to Gibson et al. [24]:
N = N 0 + ( N m a x N 0 ) × exp ( exp ( μ m a x × e N m a x N 0 × ( λ t ) + 1 ) ) ,
N = N 0 + ( N m a x N 0 ) / ( 1 + exp ( μ m a x × ( λ t ) )
where N is the cell concentration (log CFU/g) at time t, N0 and Nmax represent the initial and maximum cell numbers (log CFU/g) of S. putrefaciens, respectively. λ is the lag time (h), t is real time, and μmax represents the maximum growth rate (per h).

2.4.2. Secondary Models

To describe the temperature effect on μmax and λ, a second model (square root model) was used as follows:
1 / λ = a 1 × ( T T m i n 1 ) ,
μ m a x = a 2 × ( T T m i n 2 )
where a1 and a2 are regression coefficients; T is the real temperature in °C; Tmin1 and Tmin2 are the theoretical minimum growth limits in °C.

2.5. E-Nose Analysis

Detection of the volatile compounds of the tuna flesh was performed by an electronic nose (E-nose, Fox 4000 Alpha-MOS, France). The 18 sensors are designed as follows: LY2/LG “Chlorine, fluoride, nitrogen oxide”, LY2/G “Ammonia, amines, carbon oxides”, LY2/AA “Alcohol”, LY2/gCTL “Hydrogen sulphide”, LY2/gCT “Propane and butane”, T30/1 “Polar compounds, hydrogen chloride”, P10/1 “Hydrocarbon, ammonia, chlorine”, P10/2 “Methane and ethane”, P40/1 “Fluoride and chlorine”, T70/2 “Toluene and xylene”, PA/2 “Ethanol, ammonia, amines”, P30/1 “Hydrocarbon”, P40/2 “Chlorine and fluoride”, P30/2 “Hydrogen sulfide and ketones” T40/2 “chlorine”, T40/1 “Fluoride”, TA/2 “Alcohol”. A total of 2.0 g of flesh was placed in a 10 mL vial for 10 min at 50 °C to generate balanced headspace samples. The gas flow rate was 2.5 mL/min, and the sensor cleaning time was 8 min. Then, the sensor response of the E-nose was determined as G/G0 (G0 and G represent the conductivity of the sensor exposed to the zero gas and sample gas). The measurement cycle of the E-nose of the sample was the same as above.

2.6. Headspace Solid Phase Microextraction Gas Chromatography/Mass Spectrometry (SPME-GC/MS) Analysis

According to the method of Li et al. [25] with minor modifications, 2 g of the minced sample was placed into a 20 mL glass vial and equilibrated at 40 °C for 20 min. The SPME extraction fiber was exposed to headspace for 30 min. Gas chromatography coupled with mass spectrometry (GC-MS) was used to measure the volatiles in bigeye tuna. The carrier gas was helium (high purity 99.999%), with a constant flow rate of 1 mL/min. The oven temperature program was as follows: initial temperature of 40 °C for 5 min, 5 °C/min to 120 °C, then 10 °C/min to 250 °C, and held for 5 min. Next, the volatiles were transferred to the MS system, MS source and quadrupole: 230 and 150 °C, respectively. Mass spectra were obtained within the mass range of 20–400 m/z at 70 eV. The VOCs of samples stored at 4 and 10 °C were measured. The samples at 4 °C were measured on days 4, 8 and 12, and the samples at 10 °C were measured on days 2, 4 and 6.

2.7. Statistical Analysis

The measurement experiments of pH, TVB-N, TMA, TNS, and SPME-GC/MS were repeated three times, and the EC and E-nose experiments were repeated ten times. The data of VOCs were expressed as the mean ± standard deviation. The growth kinetic model of S. putrefaciens was fitted using MATLAB 2017b (Math Works Inc., Natick, MA, USA). Pearson correlation analysis was used to evaluate the correlation between sensor response values and TNS to select appropriate sensors for predicting the dynamic growth of S. putrefaciens. Duncan’s test and Pearson correlation coefficient were performed using SPSS 19.0 (SPSS Inc., Chicago, IL, USA). PLS regression was used to predict the spoilage potential of S. putrefaciens in tuna including TVB-N, TMA, and TNS. It is well known that the predictive performance of the calibration model cannot be determined merely by the internal validation but should also be externally validated based on predictions for samples not included in the calibration test. Data measured by the E-nose and electric conductivity meter were randomly divided into a training set (developing fitted models) and a testing set (validating models) in the ratio of 7:3. The training set of response values of the E-nose sensors was used as the data for building models, and the testing set was used as independent data to verify the accuracy of models. The data of pH, TVB-N, TMA, TNS, and SPME-GC/MS were not divided into training and validation sets. This was mainly because the data obtained from the E-nose and EC meter were used as source data of rapid detection to predict spoilage indicators (including TVB-N, TMA, TNS, and SPME-GC/MS) in the samples. Therefore, the spoilage indicators were used as a predicted object without being divided into training and validation sets (each indicator was considered as a whole). In addition, TNS was also one of the predicted indicators and was not divided into training and validation sets for the growth curve fitting. The model accuracy was evaluated by determination coefficient (R2) and root-mean-square error (RMSE). RMSE was calculated as follows:
R M S E = 1 n i = 1 n ( y m y p ) 2
where ym and yp are measured and predicted values.

3. Results and Discussion

3.1. TVB-N and TMA

Changes in the TVB-N content of bigeye tuna blocks are shown in Figure 1A. The TVB-N values increased in inoculated tuna blocks throughout the storage period at different temperatures. TVB-N values of all groups showed a slow increase followed by a rapid increase. This result was consistent with Li et al. [26], who reported that when S. putrefaciens was inoculated into blunt snout bream flesh stored at 4 °C for 14 days, TVB-N values of samples were slow on the first 8 days, but increased rapidly on the last 4 days.
The TMA value increased gradually during storage at 4, 7, and 10 °C (Figure 1B). Similar to the variation pattern of TVB-N values, the TMA values of inoculated tuna samples stored at 10 °C for 6 days were higher than those stored at 4 and 7 °C. This result indicated that temperature is a vital factor for microbial activity. The possible reason for the significant difference between the results for the CK and inoculated group may be that the production of TMA was promoted by S. putrefaciens in the inoculated samples, and the type of bacteria determined the ability to produce TMA. The TMA of samples increased exponentially with storage time, which was in agreement with [27], who reported that the change in TMA of yellowfin tuna fitted an exponential growth during the refrigeration.

3.2. pH and EC

The changes in the pH of aquatic products were closely associated with a series of chemical reactions caused by endogenous enzymes and microorganisms [28]. Changes in pH values in tuna samples at different temperatures are presented in Figure 1C. The initial pH value was 6.21 with a decreasing and then an increasing trend. A decrease in pH of the samples was due to the generation of lactic acid and the release of inorganic phosphate by decomposition of ATP [1], while the increase in the pH value was related to the release of alkaline decomposition products, such as histamine and TMA produced by spoilage microorganisms [29]. The pH of the sterile fish blocks (CK) stored at 4 °C fluctuated around 6.21, while those stored at 7 and 10 °C increased slowly. This may be due to the growth and metabolism of residual microorganisms in the CK group, and a similar phenomenon was observed in TVB-N and TMA.
EC can be used to detect meat quality from the efflux and excessive breakdown of body fluids from fish tissue during storage. Variations in EC value during cold storage are presented in Figure 1D. Initially, the EC in samples was 1057 μS/cm, and the EC of each group of samples increased significantly during storage. The change rate of EC values was higher for samples stored at higher temperatures. The significant increase in EC may be due to the autolysis of tuna cells after death and the decomposition of various nutrients in the cells into ions and small molecule metabolites with electrical conductivity under the action of enzymes and microorganisms, which enhances the electrical conductivity of the cell leachate [30]. Similar to the TVB-N and TMA curves, EC values increased slowly in the CK groups stored at 4, 7, and 10 °C.

3.3. Results of the E-Nose Analysis

As shown in Figure 2, odor maps were obtained with the E-nose from samples stored at 4, 7 and 10 °C. Response values from each sensor represent the average of 10 measurements. The E-nose responses of bigeye tuna samples stored at different storage temperatures showed a similar trend. Furthermore, sensors T70/2, PA/2, P30/1, and P30/2, which were sensitive to aromatic compounds, amines, hydrocarbons, and hydrogen sulfide, increased significantly at all storage temperatures. However, the signal values of sensors PA/2 and P30/2 increased more remarkably at 7 and 10 °C than those at 4 °C, indicating that storage temperature was an essential factor affecting the production of some compounds of tuna in storage. This finding was consistent with other studies that found temperature to be an important factor influencing the production of metamorphic substances in fish during refrigeration [31]. Changes in sensor signal values over time at different temperatures may be due to an increase in volatile gas concentrations and the production of new gas species [32], which may be related to the growth temperature and the number of S. putrefaciens.
The corresponding values from the sensor arrays were explored to determine the differences in volatilization patterns of tuna during different cold storage periods using principal component analysis (PCA). To test whether the electronic nose could distinguish samples with different storage times, PCA was employed to investigate the feasibility of distinguishing tuna inoculated with S. putrefaciens sampled at different times and temperatures based on E-nose signals. As shown in Figure 3, the differences in tuna with different storage times can be represented using two main principal components (PCs), which accounted for 88.95% (4 °C), 93.78% (7 °C), and 80.8% (10 °C) of the total variance in PC1 and 5.22%, 4.71%, and 14.38% in PC2. Regarding the samples stored at 4 °C, the data points of groups 0, 2 and 4 d were placed in the first cluster, whereas the other groups were separated into another three clusters (6, 8 and 10–12 d). The data points of groups at 0, 2, and 4 d were similar, indicating that the odor profiles of the early contaminated samples were relatively similar. But the samples at 6, 8, 10, 12 d were clearly classified by PCA. For the samples stored at 7 °C, groups at 0 and 2 d were mixed in a cluster, and groups at 4 and 6 d had a clear right downshift along the ordinate (PC 1), located into the second cluster. Groups at 8 h were separated into another cluster along the abscissa (PC2), located away from other clusters. In 10 °C groups, the data of groups at 0–2 d located into the first cluster, and another two clusters contained 3–5 d and 6 d, respectively. For samples stored at 7 and 10 °C, PCA analysis indicated that the data points of the samples at the beginning of storage (0 and 2 d for 7 °C, 0 and 1 d for 10 °C) overlapped and were initially difficult to distinguish but could be distinguished at later time points. The samples stored at 4, 7, and 10 °C showed that the E-nose based on PCA analysis poorly identified very early contaminated samples but had good overall resolution for tuna inoculated with S. putrefaciens. Therefore, E-nose could be employed as a promising approach to realize the prediction of dynamics and spoilage of Shewanella putrefaciens.

3.4. Dynamic Growth of S. putrefaciens in Tuna

The changes in the viable count of bigeye tuna are shown in Figure 1E. Gompertz and logistic models were fitted to the dynamic growth of S. putrefaciens. As shown in Table 1, the high R2 and low RMSE values indicated a good fit of models, with ranges of 0.985–0.999 and 0.0654–0.283, respectively. The lag time of the lag phase (λ) and the maximum specific growth rate (μmax) were two vital parameters for predicting the growth of microorganisms [33]. λ is especially vital to monitor food microorganisms and ensure food safety [34]. In our work, it was clear that the value of λ for S. putrefaciens decreased with increasing temperature, while μmaxe was the opposite (Table 1). We observed that the storage temperature had a significant effect on the growth of S. putrefaciens, with slow growth rates at 4 and 7 °C, while a significant growth was observed at 10 °C based on λe and μmaxe values. The growth of S. putrefaciens in tuna was in agreement with previous studies [35]. In addition, the initial colony count (N0) of S. putrefaciens in tuna ranged from 2.3 to 3.3 log CFU/g; the value of the maximum colony count (Nmax) peaked between 8.5 and 10.4 log CFU/g, according to two growth models.
However, the primary models cannot estimate the effect of temperature on the growth of S. putrefaciens in tuna, but the secondary model can evaluate it. Therefore, the square root model was used to describe the relationship between the growth parameters ( μ m a x e and λ e 1 ) and the storage temperature for the microbe. As shown in Figure 4, the results showed a strong linear correlation between the kinetic parameters and the storage temperature, with R2 values higher than 0.98. Our study also predicted the minimum growth temperature (Tmin) of S. putrefaciens in tuna in the range of −8.5 to −4.6 °C based on a secondary model, which was slightly higher than −11.4 °C predicted by [36]. This may be due to the difference in fish samples and handling.

3.5. Modeling the Kinetics of S. putrefaciens in Tuna with E-Nose Sensors

In this study, Pearson correlation analysis was used to determine the correlation between the signal values of sensors and the number of S. putrefaciens colonies in tuna. Finally, sensor P30/2 was selected (data not shown). In addition, the relatively high response values of sensor P30/2 and the considerable variation with storage time indicated its sensitivity to tuna samples during storage. Therefore, in the same way, an attempt was made to simulate the growth of S. putrefaciens by fitting the response of the sensor with Gompertz and logistic functions. Sensor P30/2 responses were fitted via Gompertz and logistic models to simulate S. putrefaciens growth according to the training set, and the validation set was used to verify the quality of the prediction models. The parameters of the generated mathematical equations are shown in Table 2, and the λe and μmaxe of CFU were derived from Table 1. The sensor fitted both models well, with high Rc2 and low RMSEc values, in a range of 0.971–0.994 and 0.0208–0.0472, respectively. Validated with the testing set, the fitting models were credible, with similar high Rp2 and low RMSEp values of 0.963–0.987 and 0.0301–0.0613, respectively.
In addition, the P/30 λs and μmaxs values obtained by the sensor based on the logistic model fit were closer to the results obtained from the actual growth model compared to the Gompertz model (Table 2). It indicated that the logistic model was more suitable for predicting of the growth of S. putrefaciens by the gas sensor in our study, which was different from the result reported by Gu et al. [34]. This difference may be due to differences in the gas sensors selected and strains. Kinetic parameters (λs and μmaxs) are of special interest in predicting microbiology and are of high practical value in monitoring food quality and safety [37]. A high correlation coefficient (r) was obtained by comparing the growth fit models generated by sensor P/30 with the models in Table 1, which indicated that the response changes of sensor P/30 to the sample were similar to the growth of S. putrefaction. Microbial kinetic models according to microbial counting methods often have difficulty in obtaining λs and μmaxs of microbial growth because of long training times [38]. Our study showed that a microbial odor response sensor may be used to simulate the dynamics of S. putrefaciens in tuna, but this requires more experimental verification.
A secondary model of the gas sensor P30/2 was also established and showed a good fit with R2 in the range of 0.965–0.998, and the linear relationship based on logistic equation was slightly stronger than that of the Gompertz equation for λs and μmaxs (Figure 5). Most studies predicted the number of microorganisms based on the chemometric method using multiple E-nose sensors [31,39,40], but rarely predicted the dynamic growth of spoilage microorganisms in aquatic products by a single sensor. In our study, primary and secondary kinetic models of S. putrefaciens in tuna were fitted by a single sensor, and the kinetic parameters were obtained.

3.6. Prediction of the Spoilage of S. putrefaciens in Tuna

To predict the spoilage level of inoculated bigeye tuna, the PLS algorithm was used to evaluate the correlation between the E-nose responses and spoilage indicators (TVB-N, TMA, and TNS) of the tuna. In addition, the changes in tissues caused by fish spoilage can be reflected by electrical conductivity [41,42], and the measurement method of EC is relatively simple. Therefore, the sensor responding values of E-nose and EC values were also combined to predict the spoilage of tuna. For PLS regression modeling, E-nose and EC measurement values from each group were randomly separated into two sets: seven samples used as the calibration set for and the remaining three samples as the prediction set. The leave-one-out cross-validation method was used to improve the accuracy of the PLS model.
Scatter graphs of the count of S. putrefaciens in tuna stored at 4 °C based on the PLS with and without EC values are given in Figure 6. The R2 and RMSE values between the predicted and experimental values are shown in Table 3. It is well known that the key to evaluating the quality of a predictive model is not only internal validation (calibration) but also external validation of samples not included in the calibration set [43]. In the two sets, the predictions of the PLS models with and without EC values performed well, and the PLS model with EC values was better than those without EC values, except for the prediction of TMA values for tuna stored at 10 °C (but the difference was not significant). The reason for this result was that the accuracy of the PLS prediction model improved by adding the index (EC) that had a great correlation with the prediction object [44].

3.7. Volatile Compounds in Tuna Samples According to HS-SPME/GC-MS

In this work, a total of 30 VOCs were detected in the control and inoculated groups (Table 4). The changing trend of VOCs at 10 °C was greater than that at 4 °C, and various VOCs in this study were previously reported to be products of protein or lipid oxidation metabolism [45]. To exclude the effect of oxidation of the fish itself on VOCs, this study also measured the VOCs produced by non-inoculated tuna blocks during storage.
The increasing tendency of 1-Penten-3-ol, 1-Octen-3-ol, and 2-Hexen-1-ol, (Z)- was observed in tuna during storage, and their increase was associated with auto-oxidative derivatization of polyunsaturated fatty acids [46]. Ethanol was present only in the early stages of tuna storage (day 0, day 4–4, and day 2–10) and was not detected as the storage period increased, which was similar to the results of Liu et al. [10]. This may be related to the fact that the metabolic process of Shewanella does not produce ethanol but can use it [4]. In addition, hexanal, heptanal, octanal, nonanal, and propanal were detected at the early stage of storage, and the above substances were confirmed to be produced by fat oxidation and had a fishy odor. The content of hexanal increased with storage time and storage temperature, which may indicate that the fat in tuna was oxidized. The changes of aldehydes reflect the degree of oxidation of polyunsaturated fatty acids such as linoleic acid in bigeye tuna, which can be used as a basis for judging the freshness of tuna. Ketones including 2-nonanone, 2-undecanone, and 2-heptanone were significantly higher in inoculated tuna compared to the control group. These ketones may originate from lipolysis and dehydrogenation by spoilage bacteria [47], which explains the low ketone content of fresh fish samples (day 0). Some hydrocarbons were also detected in this study, which were mainly derived from the decomposition of alkoxy radicals of fatty acids. Various hydrocarbons were present in the volatiles of crustaceans and fish, but they had a high threshold and made little contribution to the flavor of bigeye tuna [48]. Methylamine, N, N-dimethyl-(so-called trimethylamine) were detected only after 8–12 days at 4 °C and 4–6 days at 10 °C, which may be because of higher concentrations of TMA in the inoculated tuna during the late storage period. Dimethyl disulfide was detected in the inoculated tuna at the late storage period. This compound derived from the methionine catabolism that was produced by S. putrefaciens [49].

3.8. Relationship between E-Nose Results and VOCs

In this work, characteristic VOCs, including alcohols, aldehydes, ketones, amines, and sulfide compounds, played a significant role in distinguishing tuna samples infected with S. putrefaciens by GC-MS at different temperatures. Sensor P30/2 was sensitive to hydrogen sulfide and ketone. Comparing the VOCs with the sensor P30/2, the presence of alcohols, ketones, amines, and sulfide compounds had a significant impact on the sensor response. The significantly increased response values of the P30/2 were consistently correlated with the increased concentrations of ketones and sulfides. In addition, the relationship between sensor responses and VOCs was investigated by Pearson correlation coefficients (Figure 7). Correlation coefficient results indicated that the sensor responses were positively correlated with 1-octen-3-ol, 2-nonanone, 2-heptanone, dimethyl disulfide, and methylamine, N, N-dimethyl-, which revealed that the P30/2 was sensitive to representative VOCs of tuna inoculated with S. putrefaciens. In particular, dimethyl disulfide, and methylamine, N, N-dimethyl-, as the characteristic volatile compounds of S. putrefaciens [10], contained high content, which contributed significantly to the high signals of P30/2 to inoculated tuna. Therefore, the selected sensor could be used to discriminate tuna infected with S. putrefaciens through the specific response to characteristic VOCs.

4. Conclusions

This work demonstrated that the growth of S. putrefaciens in tuna samples stored at 4, 7, and 10 °C was consistent with two primary kinetic models (Gompertz and logistic), with high R2 and low RMSE values in the range of 0.985–0.999 and 0.0654–0.283, respectively, as well as a secondary kinetic model with high R2 values in the range of 0.9859–0.9958. The selected sensor P30/2 accurately predicted the dynamic growth of S. putrefaction. In addition, the secondary model was used to characterize the relationship between the storage temperature of the samples and the growth kinetic parameters of S. putrefaciens. The secondary model fitted with the sensor P30/2 accurately estimated the influence of temperature on the kinetic parameters of S. putrefaciens and the minimum growth temperature range of S. putrefaciens. The PLS model based on the E-nose response values with the EC values was more accurate than the model without the EC values in predicting the spoilage of tuna inoculated with S. putrefaciens. Based on the GC-MS analysis, several alcohols, ketones, amines, and sulfide compounds, especially 1-octen-3-ol, 2-nonanone, 2-heptanone, dimethyl disulfide, and methylamine, N, N-dimethyl- were determined as characteristic VOCs in tuna infected with S. putrefaciens stored at 4 and 10 °C. These results revealed that the E-nose can have a wide range of applications for predicting the growth of spoilage microorganisms and performing a quantitative analysis of spoilage in tuna.

Author Contributions

Z.Y.: Conceptualization, Methodology, Software, Investigation, Writing. J.X.: Validation, Formal analysis, Writing—review and editing, Examination, Funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Key R&D Program of China (2019YFD0901604), the Key Project of the Science and Technology Commission of Shanghai Municipality (19DZ1207503), and the Shanghai Professional Technology Service Platform on Cold Chain Equipment Performance and Energy Saving Evaluation (19DZ2284000). All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Authors can confirm that all relevant data are included in the article.

Acknowledgments

We would like to thank the anonymous reviewers for their helpful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaewprachu, P.; Osako, K.; Benjakul, S.; Suthiluk, P.; Rawdkuen, S. Shelf life extension for Bluefin tuna slices (Thunnus thynnus) wrapped with myofibrillar protein film incorporated with catechin-Kradon extract. Food Control 2017, 79, 333–343. [Google Scholar] [CrossRef]
  2. Rossi, S.; Lee, C.; Ellis, P.; Pivarnik, L. Biogenic amines formation in bigeye tuna steaks and whole skipjack tuna. J. Food Sci. 2002, 67, 2056–2060. [Google Scholar] [CrossRef]
  3. Ruiz-Capillas, C.; Moral, A. Sensory and biochemical aspects of quality of whole bigeye tuna (Thunnus obesus) during bulk storage in controlled atmospheres. Food Chem. 2005, 89, 347–354. [Google Scholar] [CrossRef]
  4. Parlapani, F.F.; Mallouchos, A.; Haroutounian, S.A.; Boziaris, I.S. Microbiological spoilage and investigation of volatile profile during storage of sea bream fillets under various conditions. Int. J. Food Microbiol. 2014, 189, 153–163. [Google Scholar] [CrossRef] [PubMed]
  5. Wright, M.H.; Shalom, J.; Matthews, B.; Greene, A.C.; Cock, I.E. Terminalia ferdinandiana Exell: Extracts inhibit Shewanella spp. growth and prevent fish spoilage. Food Microbiol. 2019, 78, 114–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lee, B.-H.; Wu, S.-C.; Shen, T.-L.; Hsu, Y.-Y.; Chen, C.-H.; Hsu, W.-H. The applications of Lactobacillus plantarum-derived extracellular vesicles as a novel natural antibacterial agent for improving quality and safety in tuna fish. Food Chem. 2021, 340, 128104. [Google Scholar] [CrossRef]
  7. Grande-Tovar, C.D.; Serio, A.; Delgado-Ospina, J.; Paparella, A.; Rossi, C.; Chaves-López, C. Chitosan films incorporated with Thymus capitatus essential oil: Mechanical properties and antimicrobial activity against degradative bacterial species isolated from tuna (Thunnus sp.) and swordfish (Xiphias gladius). J. Food Sci. Technol. 2018, 55, 4256–4265. [Google Scholar] [CrossRef] [PubMed]
  8. Serio, A.; Fusella, G.C.; Lopez, C.C.; Sacchetti, G.; Paparella, A. A survey on bacteria isolated as hydrogen sulfide-producers from marine fish. Food Control 2014, 39, 111–118. [Google Scholar] [CrossRef]
  9. Itoh, D.; Koyachi, E.; Yokokawa, M.; Murata, Y.; Murata, M.; Suzuki, H. Microdevice for on-site fish freshness checking based onk-value measurement. Anal. Chem. 2013, 85, 10962–10968. [Google Scholar] [CrossRef]
  10. Liu, X.; Huang, Z.; Jia, S.; Zhang, J.; Li, K.; Luo, Y. The roles of bacteria in the biochemical changes of chill-stored bighead carp (Aristichthys nobilis): Proteins degradation, biogenic amines accumulation, volatiles production, and nucleotides catabolism. Food Chem. 2018, 255, 174–181. [Google Scholar] [CrossRef]
  11. Odeyemi, O.A.; Burke, C.M.; Bolch, C.C.; Stanley, R. Seafood spoilage microbiota and associated volatile organic compounds at different storage temperatures and packaging conditions. Int. J. Food Microbiol. 2018, 280, 87–99. [Google Scholar] [CrossRef]
  12. Wang, X.-Y.; Xie, J. Study on the volatile organic compounds and its correlation with water dynamics of bigeye tuna (Thunnus obesus) during cold storage. Molecules 2019, 24, 3119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wu, L.; Pu, H.; Sun, D.-W. Novel techniques for evaluating freshness quality attributes of fish: A review of recent developments. Trends Food Sci. Technol. 2019, 83, 259–273. [Google Scholar] [CrossRef]
  14. Rasekh, M.; Karami, H.; Wilson, A.; Gancarz, M. Classification and identification of essential oils from herbs and fruits based on a MOS electronic-nose technology. Chemosensors 2021, 9, 142. [Google Scholar] [CrossRef]
  15. Semeano, A.T.; Maffei, D.F.; Palma, S.; Li, R.W.; Franco, B.D.; Roque, A.C.; Gruber, J. Tilapia fish microbial spoilage monitored by a single optical gas sensor. Food Control 2018, 89, 72–76. [Google Scholar] [CrossRef] [PubMed]
  16. Luan, L.; Fu, S.; Yuan, C.; Ishimura, G.; Chen, S.; Chen, J.; Hu, Y. Combined effect of superchilling and tea polyphenols on the preservation quality of hairtail (Trichiurus haumela). Int. J. Food Prop. 2017, 20, S992–S1001. [Google Scholar] [CrossRef] [Green Version]
  17. Heising, J.; Bartels, P.; van Boekel, M.; Dekker, M. Non-destructive sensing of the freshness of packed cod fish using conductivity and pH electrodes. J. Food Eng. 2014, 124, 80–85. [Google Scholar] [CrossRef]
  18. Sun, J.; Zhang, R.; Zhang, Y.; Liang, Q.; Li, G.; Yang, N.; Xu, P.; Guo, J. Classifying fish freshness according to the relationship between EIS parameters and spoilage stages. J. Food Eng. 2018, 219, 101–110. [Google Scholar] [CrossRef]
  19. Yao, L.; Luo, Y.; Sun, Y.; Shen, H. Establishment of kinetic models based on electrical conductivity and freshness indictors for the forecasting of crucian carp (Carassius carassius) freshness. J. Food Eng. 2011, 107, 147–151. [Google Scholar] [CrossRef]
  20. Song, Y.; Liu, L.; Shen, H.; You, J.; Luo, Y. Effect of sodium alginate-based edible coating containing different anti-oxidants on quality and shelf life of refrigerated bream (Megalobrama amblycephala). Food Control 2011, 22, 608–615. [Google Scholar] [CrossRef]
  21. Wang, X.-Y.; Xie, J. Growth kinetics and spoilage potential of co-culturing acinetobacter johnsonii and pseudomonas fluorescens from bigeye tuna (Thunnus obesus) during refrigerated storage. Curr. Microbiol. 2020, 77, 1637–1646. [Google Scholar] [CrossRef]
  22. Dyer, W.J. Amines in Fish Muscle: I. Colorimetric determination of trimethylamine as the picrate salt. J. Fish. Res. Board Can. 1945, 6d, 351–358. [Google Scholar] [CrossRef]
  23. Qian, Y.-F.; Xie, J.; Yang, S.-P.; Wu, W.-H.; Xiong, Q.; Gao, Z.-L. In vivo study of spoilage bacteria on polyphenoloxidase activity and melanosis of modified atmosphere packaged Pacific white shrimp. Food Chem. 2014, 155, 126–131. [Google Scholar] [CrossRef]
  24. Gibson, A.M.; Bratchell, N.; Roberts, T.A. The effect of sodium chloride and temperature on the rate and extent of growth of Clostridium botulinum type A in pasteurized pork slurry. J. Appl. Bacteriol. 1987, 62, 479–490. [Google Scholar] [CrossRef]
  25. Li, D.; Zhang, J.; Song, S.; Feng, L.; Luo, Y. Influence of heat processing on the volatile organic compounds and microbial diversity of salted and vacuum-packaged silver carp (Hypophthalmichthys molitrix) fillets during storage. Food Microbiol. 2018, 72, 73–81. [Google Scholar] [CrossRef] [PubMed]
  26. Li, Y.; Jia, S.; Hong, H.; Zhang, L.; Zhuang, S.; Sun, X.; Liu, X.; Luo, Y. Assessment of bacterial contributions to the biochemical changes of chill-stored blunt snout bream (Megalobrama amblycephala) fillets: Protein degradation and volatile organic compounds accumulation. Food Microbiol. 2020, 91, 103495. [Google Scholar] [CrossRef]
  27. Sofra, C.; Tsironi, T.; Taoukis, P.S. Modeling the effect of pre-treatment with nisin enriched osmotic solution on the shelf life of chilled vacuum packed tuna. J. Food Eng. 2018, 216, 125–131. [Google Scholar] [CrossRef]
  28. Lyu, F.; Gao, F.; Ding, Y. Effects of gamma radiation combined with cinnamon oil on qualities of smoked salmon slices inoculated with Shewanella putrefaciens. Food Sci. Nutr. 2018, 6, 806–813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Feng, X.; Ng, V.K.; Mikš-Krajnik, M.; Yang, H. Effects of fish gelatin and tea polyphenol coating on the spoilage and degradation of myofibril in fish fillet during cold storage. Food Bioprocess Technol. 2017, 10, 89–102. [Google Scholar] [CrossRef]
  30. Shen, S.; Jiang, Y.; Liu, X.; Luo, Y.; Gao, L. Quality assessment of rainbow trout (Oncorhynchus mykiss) fillets during super chilling and chilled storage. J. Food Sci. Technol. 2015, 52, 5204–5211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Shi, C.; Yang, X.; Han, S.; Fan, B.; Zhao, Z.; Wu, X.; Qian, J. Nondestructive prediction of tilapia fillet freshness during storage at different temperatures by integrating an electronic nose and tongue with radial basis function neural networks. Food Bioprocess Technol. 2018, 11, 1840–1852. [Google Scholar] [CrossRef]
  32. Sun, J.; Wang, Q.J.; Huang, J.; Hou, Y.D.; Chen, Y.F.; Su, X.R. Influence of heating temperature on the development of volatile compounds in bigeye tuna meat (Thunnus obesus) as assessed by E-nose and SPME-GC/MS. Int. Food Res. J. 2013, 20, 3077–3083. [Google Scholar]
  33. Medveďová, A.; Šipošová, P.; Mančušková, T.; Valík, L. The effect of salt and temperature on the growth of fresco culture. Fermentation 2018, 5, 2. [Google Scholar] [CrossRef] [Green Version]
  34. Gu, X.; Feng, L.; Zhu, J.; Li, Y.; Tu, K.; Dong, Q.; Pan, L. Application of gas sensors for modelling the dynamic growth of Pseudomonas in pork stored at different temperatures. Meat Sci. 2021, 171, 108282. [Google Scholar] [CrossRef] [PubMed]
  35. Yang, Z.-Q.; Tao, X.-Y.; Zhang, H.; Rao, S.-Q.; Gao, L.; Pan, Z.-M.; Jiao, X.-A. Isolation and characterization of virulent phages infecting Shewanella baltica and Shewanella putrefaciens, and their application for biopreservation of chilled channel catfish (Ictalurus punctatus). Int. J. Food Microbiol. 2019, 292, 107–117. [Google Scholar] [CrossRef] [PubMed]
  36. Taoukis, P.S.; Koutsoumanis, K.; Nychas, G.J.E. Use of time-temperature integrators and predictive modelling for shelf life control of chilled fish under dynamic storage conditions. Int. J. Food Microbiol. 1999, 53, 21–31. [Google Scholar] [CrossRef] [Green Version]
  37. Tango, C.N.; Khan, I.; Park, Y.S.; Oh, D.H. Growth of Staphylococcus aureus in cooked ready-to-eat ground fish as affected by inoculum size and potassium sorbate as food preservative. LWT Food Sci. Technol. 2016, 71, 400–408. [Google Scholar] [CrossRef]
  38. Ding, T.; Liao, X.-Y.; Dong, Q.-L.; Xuan, X.-T.; Chen, S.-G.; Ye, X.-Q.; Liu, D.-H. Predictive modeling of microbial single cells: A review. Crit. Rev. Food Sci. Nutr. 2018, 58, 711–725. [Google Scholar] [CrossRef] [PubMed]
  39. Wang, H.; Hu, Z.; Long, F.; Guo, C.; Yuan, Y.; Yue, T. Early detection of Zygosaccharomyces rouxii—Spawned spoilage in apple juice by electronic nose combined with chemometrics. Int. J. Food Microbiol. 2016, 217, 68–78. [Google Scholar] [CrossRef] [PubMed]
  40. Feng, H.; Zhang, M.; Liu, P.; Liu, Y.; Zhang, X. Evaluation of IoT-enabled monitoring and electronic nose spoilage detection for salmon freshness during cold storage. Foods 2020, 9, 1579. [Google Scholar] [CrossRef] [PubMed]
  41. Zhu, S.; Luo, Y.; Hong, H.; Feng, L.; Shen, H. Correlation between electrical conductivity of the gutted fish body and the quality of bighead carp (Aristichthys nobilis) heads stored at 0 and 3 °C. Food Bioprocess Technol. 2013, 6, 3068–3075. [Google Scholar] [CrossRef]
  42. Shi, C.; Cui, J.; Luo, Y.; Zhu, S.; Zhou, Z. Post-mortem changes of silver carp (Hypophthalmichthys Molitrix) stored at 0 °C assessed by electrical conductivity. Int. J. Food Prop. 2015, 18, 415–425. [Google Scholar] [CrossRef] [Green Version]
  43. Huang, X.-C.; Yuan, Y.-H.; Wang, X.-Y.; Jiang, F.-H.; Yue, T.-L. Application of electronic nose in tandem with chemometric analysis for detection of Alicyclobacillus acidoterrestris-spawned spoilage in apple juice beverage. Food Bioprocess Technol. 2015, 8, 1295–1304. [Google Scholar] [CrossRef]
  44. Prabhakar, P.K.; Vatsa, S.; Srivastav, P.P.; Pathak, S.S. A comprehensive review on freshness of fish and assessment: Analytical methods and recent innovations. Food Res. Int. 2020, 133, 109157. [Google Scholar] [CrossRef] [PubMed]
  45. Edirisinghe, R.K.; Graffham, A.J.; Taylor, S.J. Characterisation of the volatiles of yellowfin tuna (Thunnus albacares) during storage by solid phase microextraction and GC–MS and their relationship to fish quality parameters. Int. J. Food Sci. Technol. 2007, 42, 1139–1147. [Google Scholar] [CrossRef]
  46. Zhang, Q.; Ding, Y.; Gu, S.; Zhu, S.; Zhou, X.; Ding, Y. Identification of changes in volatile compounds in dry-cured fish during storage using HS-GC-IMS. Food Res. Int. 2020, 137, 109339. [Google Scholar] [CrossRef] [PubMed]
  47. Ercolini, D.; Russo, F.; Nasi, A.; Ferranti, P.; Villani, F. Mesophilic and psychrotrophic bacteria from meat and their spoilage potential in vitro and in beef. Appl. Environ. Microbiol. 2009, 75, 1990–2001. [Google Scholar] [CrossRef] [Green Version]
  48. Zhang, Y.; Ma, X.; Dai, Z. Comparison of nonvolatile and volatile compounds in raw, cooked, and canned yellowfin tuna (Thunnus albacores). J. Food Process. Preserv. 2019, 43, e14111. [Google Scholar] [CrossRef]
  49. Wang, H.; Liu, X.; Zhang, Y.; Lu, H.; Xu, Q.; Shi, C.; Luo, Y. Spoilage potential of three different bacteria isolated from spoiled grass carp (Ctenopharyngodon idellus) fillets during storage at 4 °C. LWT Food Sci. Technol. 2017, 81, 10–17. [Google Scholar] [CrossRef]
Figure 1. Changes in total volatile base nitrogen (TVB-N) (A), trimethylamine (TMA) (B), pH (C), Electrical conductivity (EC) (D), and total number of S. putrefaciens (TNS) of inoculation groups and total viable count (TVC) of control check (CK) groups (E) in bigeye tuna stored at different temperatures (each point is the mean value of three determinations). CK4, CK7, and CK10 are CK groups stored at 4, 7, and 10 °C. a–g in the same column with different superscripts are significantly different (p  <  0.05). The arrow indicates that the viable count was below 2.0 log CFU/g.
Figure 1. Changes in total volatile base nitrogen (TVB-N) (A), trimethylamine (TMA) (B), pH (C), Electrical conductivity (EC) (D), and total number of S. putrefaciens (TNS) of inoculation groups and total viable count (TVC) of control check (CK) groups (E) in bigeye tuna stored at different temperatures (each point is the mean value of three determinations). CK4, CK7, and CK10 are CK groups stored at 4, 7, and 10 °C. a–g in the same column with different superscripts are significantly different (p  <  0.05). The arrow indicates that the viable count was below 2.0 log CFU/g.
Foods 10 02132 g001
Figure 2. Average responses of 18 sensors in tuna samples inoculated with S. Putrefaciens at different temperatures during storage.
Figure 2. Average responses of 18 sensors in tuna samples inoculated with S. Putrefaciens at different temperatures during storage.
Foods 10 02132 g002
Figure 3. PCA score plots based on electronic nose measurements of tuna inoculated with S. Putrefaciens with different incubation times.
Figure 3. PCA score plots based on electronic nose measurements of tuna inoculated with S. Putrefaciens with different incubation times.
Foods 10 02132 g003
Figure 4. The secondary growth models of S. putrefaciens based on modified Gompertz (A,C) and logistic model (B,D).
Figure 4. The secondary growth models of S. putrefaciens based on modified Gompertz (A,C) and logistic model (B,D).
Foods 10 02132 g004
Figure 5. The secondary growth model of S. putrefaciens based on the response values of selected sensors P30/2. (A,B): Square root model based on Gompertz model; (C,D): Square root model based on logistic model.
Figure 5. The secondary growth model of S. putrefaciens based on the response values of selected sensors P30/2. (A,B): Square root model based on Gompertz model; (C,D): Square root model based on logistic model.
Foods 10 02132 g005
Figure 6. Reference measured data versus predicted data from the PLS models for prediction of the total number of S. putrefaciens (TNS) in tuna stored at 4 °C. (A,B): PLS models for calibration and validation set without EC values; (C,D): PLS models for calibration and validation set with EC values.
Figure 6. Reference measured data versus predicted data from the PLS models for prediction of the total number of S. putrefaciens (TNS) in tuna stored at 4 °C. (A,B): PLS models for calibration and validation set without EC values; (C,D): PLS models for calibration and validation set with EC values.
Foods 10 02132 g006
Figure 7. Correlations between sensor responses and GC-MS data (A): 4 °C; (B): 10 °C.
Figure 7. Correlations between sensor responses and GC-MS data (A): 4 °C; (B): 10 °C.
Foods 10 02132 g007
Table 1. The primary growth models of S. putrefaciens in tuna at different temperatures based on modified Gompertz and logistic equation by CFU.
Table 1. The primary growth models of S. putrefaciens in tuna at different temperatures based on modified Gompertz and logistic equation by CFU.
Fitted
Models
T/°CEquationsλe (h)μmaxe (h−1)NoNmaxR2RMSE
Gompertz4f(x) = 2.276 + 6.222 × exp(−exp(0.07551/6.222 × (68.8 − x) + 1))68.80.027782.2768.4980.9970.101
7f(x) = 3.292 + 7.149 × exp(−exp(0.1375/7.149 × (42.67 − x) + 1))42.670.050593.29210.4410.9940.1387
10f(x) = 3.252 + 6.331 × exp(−exp(0.2/6.331 × (24.36 − x) + 1))24.360.073583.2529.5830.9860.283
Logistic4f(x) = 3.053+5.783/(1 + exp(0.02028 × (166.6 − x)))166.60.020283.0538.8360.9980.0763
7f(x) = 3.092 + 6.657/(1 + exp(0.03285 × (105.9 − x)))105.90.032853.0929.7490.9990.0654
10f(x) = 2.83 + 6.48/(1 + exp(0.04478 × (63.6 − x)))63.60.044782.839.310.9850.167
Table 2. Parameters of dynamic growth models of S. putrefaciens in tuna stored at different temperatures based on P30/2 by modified Gompertz and logistic equations.
Table 2. Parameters of dynamic growth models of S. putrefaciens in tuna stored at different temperatures based on P30/2 by modified Gompertz and logistic equations.
ModelT/°CEquationsTraining SetTesting SetSensorCFU r
Rc2RMSEcRp2RMSEpcλs (h)μmaxs (h−1)λe (h)μmaxe (h−1)
Gompertz4f(x) = 0.280 + 0.401 × exp(−exp(0.0102/0.4023 × (90.05 − x) + 1))0.9780.02610.9680.030790.050.00375368.80.027780.986
7f(x) = 0.282 + 0.6125 × exp(−exp(0.0229/0.6125 × (52.27 − x) + 1))0.9940.02080.9850.033452.270.00842542.670.050590.976
10f(x) = 0.2872 + 0.509 × exp(−exp(0.03126/0.509 × (40.51 − x) + 1))0.9830.03410.9870.032440.510.0115024.360.073580.986
Logistic4f(x) = 0.2598 + 0.4307/(1 + exp(0.03033 × (156.50 − x)))0.9710.03150.9650.0329156.500.03033166.60.020280.996
7f(x) = 0.2352 + 0.658/(1 + exp(0.03812 × (95.24 − x)))0.9770.04720.9630.061395.240.03812105.90.032850.995
10f(x) = 0.2629 + 0.5553/(1 + exp(0.04864 × (60.36 − x)))0.9780.03910.9870.030160.360.0486463.60.044780.999
Table 3. Calibration and validation results in tuna stored at different temperatures based on the PLS model with and without EC values.
Table 3. Calibration and validation results in tuna stored at different temperatures based on the PLS model with and without EC values.
IndicatorsTemperatures (°C)PLS without EC ValuesPLS with EC Values
Calibration SetValidation SetCalibration SetValidation Set
Rc2RMSEcRv2RMSEvRc2RMSEcRv2RMSEv
TVB-N (mg/100 g)40.97130.93270.97750.86730.98120.75190.99020.6448
TMA (mg/100 g)0.990.19320.98620.24250.99050.1890.9880.2237
TNS (log10 CFU/mL)0.98080.26790.97020.35890.98680.2220.98430.2653
TVB-N (mg/100 g)70.98740.64610.98630.69810.99250.49660.99190.5053
TMA (mg/100 g)0.99560.14180.98510.3220.99620.13260.99660.1292
TNS (log10 CFU/mL)0.99560.15940.99230.2250.99570.15650.9950.1742
TVB-N (mg/100 g)100.99320.54110.98960.67550.99580.44850.99630.4214
TMA (mg/100 g)0.98970.2760.98710.34780.98760.30200.98260.3837
TNS (log10 CFU/mL)0.98570.27350.9690.41850.99630.13960.98640.2705
Table 4. Volatile organic compounds (VOCs) in tuna at 4 and 10 °C. Day 4–4, day 8–4, and day 12–4: inoculated tuna stored for 4, 8, and 12 days at 4 °C, respectively. Day 2–10, day 4–10, and day 6–10: inoculated tuna stored for 2, 4, and 6 days at 10 °C, respectively. CK 12–4 and CK 6–10: non-inoculated tuna stored for 12 and 6 days at 4 and 10 °C, respectively.
Table 4. Volatile organic compounds (VOCs) in tuna at 4 and 10 °C. Day 4–4, day 8–4, and day 12–4: inoculated tuna stored for 4, 8, and 12 days at 4 °C, respectively. Day 2–10, day 4–10, and day 6–10: inoculated tuna stored for 2, 4, and 6 days at 10 °C, respectively. CK 12–4 and CK 6–10: non-inoculated tuna stored for 12 and 6 days at 4 and 10 °C, respectively.
VOCsRelative Concentration (Area 10−6)
Day 0Day 4–4Day 8–4Day 12–4CK 12–4Day 2–10Day 4–10Day 6–10CK 6–10
Alcohols
1-Penten-3-ol3.01 ± 1.133.56 ± 0.563.89 ± 1.364.65 ± 1.984.54 ± 0.646.89 ± 0.988.99 ± 2.369.29 ± 1.997.98 ± 2.36
1-Octen-3-ol7.28 ± 1.8621.69 ± 0.9820.79 ± 3.2546.13 ± 6.5620.36 ± 2.5618.61 ± 0.5635.79 ± 2.1166.98 ± 5.2435.98 ± 4.65
Ethanol9.65 ± 1.231.23 ± 0.21NDNDND2.23 ± 0.56NDNDND
1-Hexanol12.56 ± 0.336.22 ± 0.652.12 ± 0.32ND1.36 ± 0.053.56 ± 0.65NDND2.98 ± 0.06
(3-Methyl-oxiran-2-yl)-methanolND0.66 ± 0.230.68 ± 0.241.65 ± 0.660.12 ± 0.04ND2.38 ± 0.563.01 ± 0.210.22 ± 0.08
2-Hexen-1-ol, (Z)-ND0.54 ± 0.041.89 ± 0.035.65 ± 0.69NDND2.32 ± 0.059.28 ± 0.32ND
2-Nonen-1-ol0.4 ± 0.021.25 ± 0.110.65 ± 0.010.33 ± 0.020.43 ± 0.030.54 ± 0.061.56 ± 0.020.51 ± 0.030.36 ± 0.02
Aldehydes
Hexanal7.86 ± 0.6519.23 ± 1.8933.21 ± 3.1264.93 ± 5.3613.54 ± 2.3622.52 ± 1.3535.26 ± 3.2271.26 ± 5.6221.22 ± 2.56
Heptanal4.22 ± 0.124.12 ± 0.639.52 ± 1.325.24 ± 0.973.22 ± 0.868.79 ± 1.221.23 ± 0.063.54 ± 0.071.98 ± 0.21
Nonanal2.21 ± 0.082.28 ± 0.214.60 ± 0.118.21 ± 0.993.21 ± 0.0718.47 ± 2.3612.11 ± 2.6515.89 ± 3.0611.21 ± 2.10
Propanal2.35 ± 0.095.61 ± 0.467.22 ± 1.236.11 ± 0.122.12 ± 0.335.31 ± 0.224.33 ± 1.116.85 ± 0.645.11 ± 0.21
2-Decenal, (E)ND3.12 ± 0.182.47 ± 0.211.56 ± 0.06NDND0.43 ± 0.052.55 ± 0.650.35 ± 0.05
2-Nonanal, (E)-ND0.22 ± 0.033.56 ± 0.244.54 ± 0.683.72 ± 0.210.31 ± 0.055.33 ± 0.985.99 ± 1.333.56 ± 0.09
2-Dodecanal, (E)-NDND1.263 ± 0.052.112 ± 0.324.6 ± 0.235.68 ± 1.224.22 ± 0.093.68 ± 1.212.23 ± 0.21
Decanal0.65 ± 0.021.36 ± 0.05NDND2.55 ± 0.651.55 ± 0.962.36 ± 0.121.32 ± 0.081.86 ± 0.09
Octanal3.86 ± 0.331.69 ± 0.350.98 ± 0.049.98 ± 2.652.65 ± 0.221.98 ± 0.215.65 ± 0.693.23 ± 0.935.36 ± 0.97
4-Heptanal, (Z)-ND0.09 ± 0.010.98 ± 0.096.22 ± 1.230.28 ± 0.020.18 ± 0.042.32 ± 0.261.69 ± 0.660.56 ± 0.08
Ketones
2,3-OctanedioneNDND0.85 ± 0.151.7 ± 0.56ND0.08 ± 0.011.23 ± 1.232.56 ± 1.56ND
2,3-Pentanedione0.04 ± 0.010.06 ± 0.020.03 ± 0.011.56 ± 0.051.78 ± 0.582.53 ± 0.041.03 ± 0.054.22 ± 1.863.29 ± 0.12
2-NonanoneNDND0.06 ± 0.020.12 ± 0.04NDND0.08 ± 0.050.25 ± 0.08ND
2-UndecanoneNDND0.75 ± 0.022.67 ± 0.660.96 ± 0.06ND1.19 ± 0.284.55 ± 1.321.89 ± 0.64
2-HeptanoneNDND0.12 ± 0.020.89 ± 0.130.06 ± 0.01ND1.35 ± 0.612.46 ± 0.050.09 ± 0.01
Hydrocarbons
HeptacosaneNDND1.23 ± 0.130.85 ± 0.042.35 ± 0.122.79 ± 0.130.46 ± 0.020.77 ± 0.071.23 ± 0.21
Pentadecane2.23 ± 0.148.63 ± 1.364.62 ± 0.883.33 ± 0.213.26 ± 1.029.59 ± 1.334.35 ± 0.982.23 ± 0.291.9 ± 0.06
Tetradecane4.23 ± 0.322.51 ± 0.110.56 ± 0.02ND1.03 ± 0.051.54 ± 0.21NDND0.81 ± 0.06
Others
Ethyl acetate5.03 ± 0.097.54 ± 1.323.21 ± 0.351.22 ± 0.352.13 ± 0.863.23 ± 0.780.77 ± 0.122.27 ± 0.651.22 ± 0.05
Methylamine, N, N-dimethyl-NDND3.31 ± 0.098.02 ± 1.04ND3.28 ± 0.3910.65 ± 2.6521.23 ± 4.65ND
Methoxy-phenyl-oximeNDND0.64 ± 0.060.46 ± 0.030.11 ± 0.02ND0.83 ± 0.042.22 ± 0.250.23 ± 0.03
Naphthalene0.07 ± 0.010.14 ± 0.020.08 ± 0.010.07 ± 0.020.12 ± 0.030.09 ± 0.010.16 ± 0.040.06 ± 0.010.16 ± 0.04
Dimethyl disulfideNDNDND0.85 ± 0.06NDNDND1.32 ± 0.31ND
ND means not detected.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yi, Z.; Xie, J. Prediction in the Dynamics and Spoilage of Shewanella putrefaciens in Bigeye Tuna (Thunnus obesus) by Gas Sensors Stored at Different Refrigeration Temperatures. Foods 2021, 10, 2132. https://0-doi-org.brum.beds.ac.uk/10.3390/foods10092132

AMA Style

Yi Z, Xie J. Prediction in the Dynamics and Spoilage of Shewanella putrefaciens in Bigeye Tuna (Thunnus obesus) by Gas Sensors Stored at Different Refrigeration Temperatures. Foods. 2021; 10(9):2132. https://0-doi-org.brum.beds.ac.uk/10.3390/foods10092132

Chicago/Turabian Style

Yi, Zhengkai, and Jing Xie. 2021. "Prediction in the Dynamics and Spoilage of Shewanella putrefaciens in Bigeye Tuna (Thunnus obesus) by Gas Sensors Stored at Different Refrigeration Temperatures" Foods 10, no. 9: 2132. https://0-doi-org.brum.beds.ac.uk/10.3390/foods10092132

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop