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

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.


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 H 2 S, 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

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.

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.

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.

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]: where N is the cell concentration (log CFU/g) at time t, N 0 and N max 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).

Secondary Models
To describe the temperature effect on µ max and λ, a second model (square root model) was used as follows: where a 1 and a 2 are regression coefficients; T is the real temperature in • C; T min1 and T min2 are the theoretical minimum growth limits in • C.

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.

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 (R 2 ) and root-mean-square error (RMSE). RMSE was calculated as follows: where y m and y p are measured and predicted values.

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.

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.

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. perature 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 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. 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.

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 R 2 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

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 R 2 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 (N 0 ) of S. putrefaciens in tuna ranged from 2.3 to 3.3 log CFU/g; the value of the maximum colony count (N max ) 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 ( √ µ maxe 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 R 2 values higher than 0.98. Our study also predicted the minimum growth temperature (T min ) 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.

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 Rc 2 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 Rp 2 and low RMSEp values of 0.963-0.987 and 0.0301-0.0613, respectively.

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 R c 2 and low RMSE c 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 R p 2 and low RMSE p 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 R 2 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.

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.

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 R 2 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].  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].

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]. 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.

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.

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 R 2 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 R 2 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, 2heptanone, 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.

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 R 2 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 R 2 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.