Muscle foods, such as meat and poultry, are described as spoiled if organoleptic changes make them unacceptable to the consumer. These organoleptic characteristics may include changes in appearance (discoloration), the development of off odors, slime formation, changes in taste, or any other characteristic which makes the food undesirable for consumption (
25,
26). While endogenous enzymatic activity within muscle tissue postmortem can contribute to changes during storage (
1,
25,
32,
44), it is generally accepted that detectable organoleptic spoilage is a result of decomposition and the formation of metabolites caused by the growth of microorganisms (
28,
40). The organoleptic changes which take place will also vary according to the species of microflora present, the characteristics of the meat, processing methods, product composition, and the environment in which the food is stored (
25).
The ideal method for the on-line microbiological analysis of meat would be rapid, noninvasive, reagentless, and relatively inexpensive, and these requirements can be met via the application of a spectroscopic approach, in combination with any appropriate data deconvolution strategy based on statistics or machine learning. Such statistical methods include partial least-squares (PLS) regression (
36), while a popular and powerful series of machine learning strategies (
37) are based on methods of evolutionary computing (
4), such as genetic algorithms (GAs) (
4,
19,
24) and genetic programming (GP) (
6,
33). Fourier transform infrared (FT-IR) spectroscopy involves the observation of vibrations of molecules that are excited by an infrared beam, and an infrared absorbance spectrum represents a “fingerprint” which is characteristic of any chemical or biochemical substance (
18,
43). This technique is also very rapid (taking seconds) and has been shown to be a valuable tool for the characterization of axenically cultured bacteria (
22,
34,
38,
47), including single-gene knockout strains (
41).
RESULTS AND DISCUSSION
The comminution of samples in order to accelerate the spoilage process was successful, as the final mean log
10(TVC) of 9.02 (Table
1) was an order of magnitude above the 10
8 CFU g
−1 generally accepted as the point at which organoleptic spoilage becomes readily detectable (
13,
15). Using 10
8 CFU g
−1 as the indicator for postspoilage, the average spoilage time over the series of experiments was 13.6 h (Table
1). The spoilage of the samples within 24 h at room temperature was anticipated, as comminution ruptures cell walls, releasing a source of nutrients; increases the surface-area-to-volume ratio; and distributes bacteria that would normally be restricted to the surface throughout the meat substrate. The initial mean pH range of fresh samples during the three experiments (5.7 to 5.9) was within those described previously in the literature (
13,
15). The use of pH as an indicator of spoilage or remaining shelf life in meats would be insufficient, as the pH fluctuates prior to spoilage, only rising significantly when levels of bacteria reach ∼10
8 CFU g
−1. At this level, sensory spoilage is readily detectable, and this is partly a consequence of the increase in pH and the production of malodorous substances, such as ammonia, dimethylsulfide, and diacetyl, by the catabolic action of the resident microflora (
13,
25).
The 150 FT-IR spectra from experiment 1 are shown in Fig.
1 and illustrate the reproducibility of both HATR as a sampling method and the experimental protocol that was undertaken over a period of 6 weeks. Typical FT-IR spectral data from the 1,750- to 700-cm
−1 wave number range from measurements of meat at the pre- and postspoilage stages are shown in Fig.
2. These spectra are from chicken breast meat carrying ∼7 · 10
6 and ∼2 · 10
9 CFU g
−1, respectively, and are both data rich and not biased to any particular group of chemicals associated with a particular group of metabolites. Furthermore, the spectra are complex and multidimensional in nature, so they do not easily lend themselves to simple visual interpretation; this is compounded by the fact that the data set for all three experiments is substantial, with a total of 450 spectra, each containing 441 wave numbers. However, with the advent of modern machine learning approaches, the opportunity now exists to analyze such complex high-dimensional spectral patterns (
7,
46) and to extract an answer to a question of biological interest with much lower dimensionality, i.e., “What is the bacterial load on the meat surface?”
Therefore, as described above, the supervised-learning method of PLS regression was calibrated and cross validated with the FT-IR spectral data and the known log
10(TVC) values from experiments 1 and 2 (Table
2 shows the details and TVC levels) before being challenged by the independent and “unseen” test set of data from experiment 3. The plots of the estimates versus the known log
10(TVC) (Fig.
3) show that the FT-IR and PLS predictions were virtually indistinguishable from the expected log proportional fit [i.e., log(
y) = log(
x)] and so show that this approach can be used to accurately assess the spoilage status of meat. As can be seen in Fig.
3 and Table
2, the lowest level of spoilage encountered was ∼2 · 10
6 CFU g
−1, and this necessarily restricts the detection limit. The TVC for chicken immediately postslaughter is 10
3 cm
−2, rising to 10
4 to 10
5 cm
−2 after packaging (
26). That PLS gave accurate results at 2 · 10
6 CFU g
−1 suggests that it will be possible to reach lower levels, and this will be the subject of further study with freshly killed chickens. From Fig.
3, it is evident that the spectra obtained by direct FT-IR analysis of meat do contain biochemical information that allows correlation with the spoilage status of the chicken, for data used to produce the PLS model and, more importantly, for data from a completely new experiment. The obvious question that needs to be addressed is that of which biochemical species the FT-IR is measuring that are related to the spoilage status of the chicken.
The Pearson correlation coefficients between the absorbances at each wave number in the FT-IR spectra from experiments 1 and 2 and the log
10(TVC) were calculated and are also plotted in Fig.
2. It can be seen that most peaks from 1,500 to 700 cm
−1 are positively correlated with spoilage, but no single peak appears uniquely dominant; this necessarily means that it is difficult to pinpoint the cause of microbial spoilage to a single (or a small group of) biochemical species using this correlation approach. Therefore, GAs and GPs were evolved to discriminate qualitatively between meat carrying <10
7 and ≥10
7 bacteria (as TVC) per cm
2.
GA-MLR was applied so as to extract subsets of two, three, and five wave numbers that could discriminate between fresh (<107 bacteria/cm2) and spoiled (≥107 bacteria/cm2) chicken. Because the starting population for each GA run was random, 60 GA-MLR runs were performed, and the following subsets were found to be optimal for selecting just two or three wave numbers, respectively: (1,096, 1,227 cm−1) and (1,312, 1,235, 1,088 cm−1). When the algorithm was used to look for five wave numbers, it was found that the degree of discrimination did not improve compared with selecting subsets of three, and no consistent areas of the FT-IR spectra were found to be dominant in the GA expressions; however, vibrations at 1,096 and 1,305 cm−1 were found within the best subsets.
GP analyses (i) using the same 10
7-bacteria/cm
−1 threshold as above and (ii) evolved to predict the log
10(TVC) levels produced trees which could easily discriminate between fresh and spoiled chicken and quantify the level of spoilage, respectively; a typical GP parse tree is shown in Fig.
4. As with the GAs, the initial populations were produced randomly; therefore, 10 separate GPs were evolved. For the threshold GP analysis, the number of times each input (wave number) was used for the 10 evolved populations was calculated and plotted against the wave number of the infrared light (Fig.
5). Figure
5 clearly shows that the dominant area of the spectra for discriminating between fresh (<10
7 bacteria/cm
2) and spoiled (≥10
7 bacteria/cm
2) chicken was 1,088 to 1,096 cm
−1; moreover, these wave numbers were also selected by the GA-MLR method. The functional group vibration in the region 1,088 to 1,096 cm
−1 is ascribable to C-N stretching, most plausibly from amines (
14,
35).
The most intense peaks that appear in fresh meat are the amide I (C⋕O vibration at 1,640 cm
−1) and amide II (N—H deformation at 1,550 cm
−1) bands from proteins and peptides, and from the Pearson correlation coefficients, the amide II band is the only vibration that is negatively correlated with spoilage (Fig.
2). This strongly suggests that the protein content of the meat was decreasing during spoilage. By contrast, the peaks at 1,240 and 1,088 cm
−1, which are both ascribable to C-N stretching from amines from free amino acids, are positively correlated. Indeed, the rule in Fig.
4 shows that spoilage can be ascribed simply to the ratio of 1,096 to 1,683 cm
−1 from these vibrations from amines and amides, respectively.
Plots of the absorbances of these vibrations versus the time for the second experiment are shown in Fig.
6. It is clear that the amide I and II bands are constant, although the amide II band does decrease very slightly after 16 h while the peaks at 1,240 and 1,088 cm
−1 start to increase significantly after 16 h. It is noteworthy that the onset of spoilage, as characterized by a TVC of >10
8 g
−1, was at 17 h, and this was the point at which the absorbance due to free amines started to increase. This was also found to be the case for experiments 1 and 3 (data not shown). These correlations, and the fact that the GAs and GPs both pick the region 1,088 to 1,096 cm
−1 as the most significant area of the FT-IR spectra for the prediction of spoilage of chicken which is attributable to free amino acids, makes it clear that the most significant metabolic process that occurs at spoilage is the start of proteolysis. This is indeed highly likely, since it is known that spoilage in meat is most frequently associated with the postglucose utilization of amino acids by aerobic microorganisms, such as pseudomonads, and the onset of the enzymatic degradation of proteins and peptides, leading to the production of free amino acids (
8,
12,
40).
In conclusion, FT-IR spectroscopy, in combination with appropriate machine learning methods, presents itself as a novel method for the quantitative detection of food spoilage. Using FT-IR, we were able to acquire a metabolic snapshot (
30) and quantify, noninvasively, the microbial loads of food samples accurately and rapidly (within 60 s) directly from the sample surface. We believe that this approach has considerable potential for further development and will aid both the food safety regulatory bodies and the Hazard Analysis Critical Control Point system. In particular, we will conduct future studies testing our method for quantifying the numbers of spoilage organisms on muscle foods at the production, processing, packaging, and storage levels.