Journal of Dairy Science
Volume 90, Issue 7 , Pages 3091-3109, July 2007

Characterization of Flavor and Texture Development Within Large (291kg) Blocks of Cheddar Cheese1

  • M.E. Carunchia Whetstine

      Affiliations

    • Department of Food Science, Southeast Dairy Foods Research Center, North Carolina State University, Raleigh 27695
  • ,
  • P.J. Luck

      Affiliations

    • Department of Food Science, Southeast Dairy Foods Research Center, North Carolina State University, Raleigh 27695
  • ,
  • M.A. Drake

      Affiliations

    • Department of Food Science, Southeast Dairy Foods Research Center, North Carolina State University, Raleigh 27695
    • Corresponding Author InformationCorresponding author.
  • ,
  • E.A. Foegeding

      Affiliations

    • Department of Food Science, Southeast Dairy Foods Research Center, North Carolina State University, Raleigh 27695
  • ,
  • P.D. Gerard

      Affiliations

    • Experimental Statistics Unit, Mississippi State University, Mississippi State 39762
  • ,
  • D.M. Barbano

      Affiliations

    • Department of Food Science, Northeast Dairy Foods Research Center, Cornell University, Ithaca, NY 14853

Received 13 November 2006; accepted 13 March 2007.

Article Outline

Abstract 

Cheddar cheese is a natural product that has a variable flavor and texture profile. Many companies produce 291-kg blocks of Cheddar cheese, which are subsequently cut and shipped, or stored and subsequently cut. Previous research has shown that compositional differences exist within 291-kg blocks and that these differences may influence flavor and texture development. The objectives of this study were to systematically characterize flavor and texture differences within 291-kg blocks. On 2 different occasions, a 291-kg block was manufactured at each of 4 manufacturing facilities. After 7 d, the 291-kg blocks were sliced into sixteen 18-kg sample portions using a predetermined diagram, and each portion was labeled appropriately (outer corner, inner corner, etc.) and stored at 7°C. Cheese from different locations within the 291-kg blocks was evaluated at 1, 4, 8, and 12 mo. At each time point, two 18-kg portions representing an inside and outside location with the 291-kg block cross-section (from inside to outside) were sampled. The moisture content was lower in the inner than outer locations within the 291-kg blocks. Protein hydrolysis was higher in the inner location and inner locations developed aged Cheddar flavors sulfur, nutty, and brothy more rapidly than the outer locations. However, plant-to-plant differences in aging were often larger than differences caused by block location. These differences were due to differences in cheese manufacturing practices among plants. Dynamic headspace results for flavor volatiles were consistent with descriptive sensory flavor results, documenting differences between inner and outer locations within 291-kg blocks. The inner locations were more fracturable and the outer locations were more cohesive and had more residual in the mouth. Inner locations had greater fracture strain than outer locations. Documenting the differences in aging of 291-kg blocks of Cheddar cheese is important in understanding how to make a consistent high-quality Cheddar cheese.

Key words: 291-kilogram block, flavor development, texture development

 

Back to Article Outline

Introduction 

Cheddar cheese is a natural product that has a variable flavor and texture profile. Many companies produce 291-kg blocks of Cheddar cheese. Making these large blocks of Cheddar rather than the more traditional 18-kg blocks reduces labor and handling costs (Mesa-Dishington et al., 1987). However, previous studies have demonstrated that compositional differences exist within 291-kg blocks (Reinbold and Ernstrom, 1988; Reinbold et al., 1992; Barbano, 2001). The process used during pressing and draining the 291-kg blocks may influence the characteristics of the moisture gradient, but generally the gradient from the exterior to the interior is larger in the bottom half of a 291-kg block (Barbano, 2001). These differences in composition may influence flavor and texture development.

In a 291-kg block, the temperature at the outer surfaces of the block will decrease more rapidly than the temperature inside the block. This temperature gradient has an impact on moisture migration (Olabi and Barbano, 2002). Previous studies (Reinbold and Ernstrom, 1988; Reinbold et al., 1992; Barbano, 2001) reported a difference in the moisture content between the inner and outer portions of the 291-kg cheese block. As the temperature of the cheese near the surface of the 291-kg block decreases, the moisture content increases because moisture migrates from the warmer cheese in the center of the block to the colder cheese at the surfaces. Mobility of the moisture appears to be directly related to the ability of CN to hold or release water, which is primarily a function of pH and temperature (Lawrence et al., 2004). Caseins are relatively hydrophobic (Swaisgood, 1992), and lower temperatures cause hydrophobic proteins to favor protein-water interactions (outer location), whereas higher temperatures favor protein-protein interactions (inner location; Lehninger, 1970). Olabi and Barbano (2002) reported that a temperature gradient within cheese caused moisture to move from an area of warm temperature (interior) to an area of cooler temperature (exterior) both with and against the force of gravity in cheeses cooled with a temperature gradient from 27 to 3°C. The result was that the outer portion of the block was higher in moisture than the inner portion (Barbano, 2001; Olabi and Barbano, 2002).

Previous research has also documented a pH gradient initially before cooling of 5.2 at the exterior and 5.38 in the interior to a pH of 5.05 exterior to 4.95 interior in 291-kg cheese blocks (Reinbold et al. 1992). Reinbold et al. (1992) found that during aging, the pH at the interior of the block decreased, whereas there was less change in pH at the exterior of the 291-kg block. Decreasing pH can increase syneresis, which will decrease the moisture content (Pastorino et al., 2003). Proteolysis within the cheese can have dramatic effects on the flavor and texture development and may be different among locations with 291-kg blocks.

There has been no systematic study characterizing flavor and texture development within 291-kg blocks. The objectives of this study were to systematically characterize flavor and texture differences across 291-kg blocks of Cheddar cheese. Instrumental and rheological methods were used as well as descriptive sensory analyses of flavor and texture.

Back to Article Outline

Materials and Methods 

Cheeses 

Blocks of pasteurized milk Cheddar cheeses (291kg) were produced at 4 commercial manufacturing facilities (2 milled-curd, 2 stirred-curd facilities). One block was sampled at each plant, and this sampling was replicated on a block manufactured at a later date in each plant. Stirred-curd facilities were located on the West coast and in the Midwestern United States. Milled-curd facilities were located in the Midwest and Northeastern United States. Blocks at each facility were cooled using forced-air coolers and kept separate from other blocks to ensure even cooling rates through all 6 faces of each block. Blocks from plant 1 were cooled at 5°C for 14 d, whereas blocks from plants 2 to 4 were cooled at 1°C for 3 to 5 d. All plants used fermentation-produced chymosin. The rennet use rates were as follows: Plant 1 used 23.1mL, plants 2 and 3 used 37.5mL, and plant 4 used 21mL per 454kg of standardized milk. Only plant 1 did not add enzymes to increase the rate of flavor development during aging. The temperature of the curd when the 291-kg blocks were filled was lower (28 to 29°C) for plants 1 and 4, whereas it was higher (32°C) for plants 2 and 3. Plants 1 and 3 did not fortify their milk, plant 2 used only a low level of fortification, and plant 4 used a very high level of fortification with nonfat milk solids and cream. When cooling the 291-kg blocks after pressing, plant 1 used forced air at 5°C, whereas the other 3 plants used 1°C.

After cooling (7 d), each 291-kg block was sliced into sixteen 18-kg portions using a predetermined diagram (Figure 1) and each portion was labeled appropriately (outer corner, inner corner, etc.). Eight 18-kg portions from the upper half of the 291-kg block were vacuum sealed and shipped by refrigerated carrier to North Carolina State University. Upon receipt, the 18-kg portions from the 291-kg blocks were examined for damage and then stored in the dark at 7°C.

  • View full-size image.
  • Figure 1. 

    Diagram of 18-kg sample portions cut from a 291-kg block. The upper- and outermost top corner of 18-kg portion C represents an outer location (i.e., the point of the arrow facing into the upper, outside corner of the block) and the lower- and bottom-most corner of 18-kg portion D represents an inner location (i.e., the point of the arrow facing into the center of the block). Blocks E to H are mirrored in the top part of the block. Blocks A, C, E, and G are outer blocks; blocks B, D, F, and H are inner blocks.

The 18-kg portions were trimmed and cut at 1, 4, 8, and 12 mo using steel wires. At each time point, 2 of the 18-kg portions were sampled such that the exterior to interior locations of the upper half of the 291-kg block were represented (locations A and B, C and D, E and F, G and H). Approximately 2kg of cheese representing an inner location (B, D, F, and H in Figure 1) were cut from each 18-kg portion and 2-kg samples of the outer location blocks (A, C, E, and G in Figure 1) were cut from each 18-kg portion for instrumental and sensory analysis. The rest of the 18-kg portions were discarded. The 2-kg samples from inner portions represented adjacent pieces from the interior of the 291-kg block and the outer 2-kg samples represented the exterior top corners of the 291-kg block. The lower half of the 291-kg block was not sampled.

Proximate Analysis 

The pH, fat, moisture, protein, and salt were measured for all (inner and outer) locations in duplicate at each time point using standard methods. Cheese pH was measured using a Xerolyt combination electrode (model HA405; Mettler Toledo, Columbus, OH) and an Accumet pH meter (model AR 25, Fisher Scientific, Pittsburgh, PA) after tempering the cheese to 23°C. Fat content was determined using the Babcock method (Marshall, 1992; method 15.8.A). Cheese moisture was determined gravimetrically in 2g of cheese in a forced-air oven at 100°C for 24h (AOAC, 2000; methods 33.2.44, 990.20) Salt content was determined using the Volhard method (Marshall, 1992; method 15.5.B). The Kjeldahl method was used to determine the total nitrogen content of cheese (Lynch et al., 2002). Crude protein was calculated by multiplying total nitrogen by 6.38. Analyses were performed in duplicate.

Proteolysis 

The o-phthaldialdehyde (OPA) assay, as outlined by Church et al. (1983), was used to determine the degree of protein hydrolysis over time during cheese aging. A 1:20 dilution of grated cheese and water was prepared by taking 1g of cheese and adding it to 19g of deionized water and vortexing. Five grams was removed and 10mL of 0.75 N TCA was added. The samples stood for 10min and were filtered through a Whatman #2 filter paper (Whatman Ltd., Maidstone, UK) to remove particles and obtain the material soluble in the TCA solution. Premade OPA reagent (Pierce, Rockford, IL; 2mL) was added to a 200-μL sample and allowed to sit for exactly 2min to allow the OPA reagent to react with primary amino groups. The absorbance was then measured at 340nm using a Shimadzu UV-260 spectrophotometer (Shimadzu, Columbia, MD). The following formulas were used to determine the degree of protein hydrolysis:

where n is the number of bonds cleaved, ɛ is 6,000 m/cm, M is the molar concentration of protein in the cheese slurry (0.5mM), F is the dilution factor in the assay procedure (0.0017), and ΔA340nm is the absorbance at 340nm of a sample compared with a reagent blank prepared in water,
where 194.2 is the weighted average of the total number of peptide bonds for CN and whey proteins in cheese. Analyses were conducted in duplicate.

At 12 mo of age, a second method (Kjeldahl soluble nitrogen) that has been more commonly used to determine proteolysis in cheese was used to cross-check the results of the OPA method. Total nitrogen and nitrogen soluble in pH 4.6 acetate buffer and 12% TCA were determined in duplicate as described by using the Kjeldahl total nitrogen method for milk (AOAC, 2000) with the sample preparation and size as described by Bynum and Barbano (1985). The pH 4.6 acetate buffer and 12% TCA-soluble nitrogen were expressed as a percentage of total nitrogen.

Dynamic Headspace Analysis GC-MS 

Cheese was frozen and grated using a hand grater. A 50-μL quantity of internal standard (50μL of 2-methyl-3-heptanone and 50μL of 2-methyl pentanoic acid in 5mL of methanol) was added to 50g of grated cheese. The cheese was then kneaded and thoroughly mixed to distribute the internal standard evenly. This reformed cheese sample was refrigerated at 4°C for overnight equilibration and then frozen at −80°C for at least 24h. Cheese was then grated again, and 10g of cheese was added to 20g of water and thoroughly mixed using a hand homogenizer. A 5-g quantity of cheese slurry was loaded into a needle sparger 25-mL purge-and-trap vial (Tekmar, Vernon, British Columbia, Canada) along with 2g of salt. The vial was equilibrated at room temperature for 30min and then placed on a CDS 6000 purge-and-trap apparatus (CDS, Oxford, PA) and purged with nitrogen at 40 mL/min. The initial purge volume was 800mL (32min). Concentrated volatile compounds were desorbed from the trap and then transferred by a heated transfer line to the gas chromatograph and desorbed onto a nonpolar DB-5MS column (30m length×0.25mm i.d.×0.25μm df; J & W Scientific, Folsom, CA) on an HP5890 Series II GC/HP 5972 mass selective detector (Hewlett-Packard, Co., Palo Alto, CA) The oven temperature was programmed at −20 to 60°C at a rate of 4°C/min with a 6-min hold at −20°C, and then 60°C to 220°C at a rate of 6°C/min with a final hold of 5min. Mass selective detector conditions were as follows: capillary direct interface temperature, 280°C; ionization energy, 70eV; mass range, 33 to 330 atomic mass units; electron multiplier voltage (Atune+200V); scan rate, 5 scans/s. Triplicate analyses were performed on each sample. Based on MS results, concentrations of 2,3-butanedione (diacetyl), 2/3-methyl butanal, 3-hydroxy-2-butanone (acetoin), methyl butanoate, ethyl butanoate, hexanal, methyl hexanoate, ethyl hexanoate, acetic acid, and butanoic acid were calculated using an external standard curve.

For positive identifications, retention indices and mass spectra were compared with those of authentic standard compounds analyzed under identical conditions. Tentative identifications were based on comparing mass spectra of unknown compounds with those in the 1992 National Institute of Standards and Technology (Gaithersburg, MD) mass spectral database or on matching the retention index values against those of authentic standards. For the calculation of retention indices, an n-alkane series was used (Van den Dool and Kratz, 1963).

Quantification of Volatile Compounds 

Response factors of selected compounds were calculated by direct addition of known amounts of standards to odor-free water prior to dynamic headspace analysis GC-MS. Response factors were determined using a 5-point standard curve (r2>0.96) on a DB-5 column using GC-MS. With these response factors, the selected compounds were quantified using the response factor and the area ratio of compound to internal standards. 2,3-Butandione (diacetyl), 2/3-methyl butanal, 3-hydroxy-2-butanone (acetoin), methyl butanoate, ethyl butanoate, hexanal, methyl hexanoate, ethyl hexanoate, acetic acid, and butanoic acid were selected for quantification because these compounds were consistently observed in the samples and have previously been shown to affect the flavor in Cheddar cheese (Christensen and Reineccius, 1995; Milo and Reineccius, 1997; Zehentbauer and Reineccius, 2002; Singh et al., 2003; Carunchia Whetstine et al., 2005). All standards were obtained from Aldrich Chemical Company (St. Louis, MO).

Sensory Evaluation of Cheeses 

At each time point (1, 4, 8, and 12 mo), cheeses were sampled for sensory analysis. For sampling, the outer 1cm of each 18-kg portion was trimmed to eliminate flavors caused by packaging or exposure. Flavor and texture were evaluated on different days in different sessions.

Flavor Analysis 

A trained (>100h each) sensory panel (n = 14) evaluated the cheeses using the flavor lexicon developed for Cheddar cheese (Drake et al., 2001). Cheese was presented in 2×2cm cubes and placed into 4-oz. (120mL) soufflé cups with 3-digit codes. Panelists were trained for 100h on flavor, aroma, and feeling factors using the Spectrum method (Meilgaard et al., 1999). The 15-point numerical Spectrum intensity scale was used to mark panelist responses. On this universal intensity scale, which can be applied to all products and flavor intensities, most Cheddar cheese flavors fall between 0 and 5 (Drake et al., 2001, 2005). During evaluation, panelists had free access to water and unsalted crackers. Four cheeses were evaluated per session. Cheeses were evaluated in duplicate by each panelist.

Texture Analysis 

A sensory panel (n = 14) with 50h of training evaluated the cheeses using a previously published texture lexicon (Brown et al., 2003). The scaling technique used for texture analysis was product specific, with a 15-point anchored and referenced scale. The scale was anchored on the right with “very” and on the left with “not.” On this scale, which is specific only to cheeses, texture values fall between 0 and 15. Panelists were provided with cheese references (Brown et al., 2003) during evaluations to minimize variability. At each session, no more than 6 samples were evaluated. Each cheese was cut into 1.27-cm3 cubes and presented at room temperature in lidded soufflé cups (to minimize moisture loss) with 3-digit random codes. At each session, panelists had free access to spring water and unsalted crackers as well as appropriate references. Each cheese was evaluated in duplicate by each panelist.

Instrumental Texture 

Large Strain Analysis 

A torsional method was used to determine fracture properties of cored cheese samples. Cheese was cored from each 18-kg portion using a cylindrical metal corer (minimum diameter, 10mm). The cheese was then cut to a length of 28.7mm. Plastic disks (Gel Consultants, Raleigh, NC) were glued to the ends of the cylinder using cyanoacrylate glue (Loctite 100, Loctite Corporation, Rocky Hill, CT) to enable the samples to be mounted to the grinding and twisting apparatuses. The cylinders were shaped into a capstan shape having a minimum diameter of 10mm using a precision grinding machine (Gel Consultants). Samples were twisted using a Haake 550 viscotester (Gebruder Haake GmbH, Karlsruhe, Germany) fitted with a fabricated apparatus that enabled torsional measurement (Troung and Daubert, 2001). Samples were twisted at a rate of 4.5rpm at 25°C. Samples were twisted at a strain rate of 0.45s−1 until fracture. Six to 8 capstans were produced and measured for each sample. After grinding, sample geometry was measured for calculations. True shear stress (σt) and true shear strain (γt-true) were calculated at fracture based on the method of Diehl et al. (1979) as described by Brown et al. (2003).

Creep and Recovery 

Creep and recovery tests were performed using a Stress Tech controlled stress rheometer (ATS Rheosystems, Bordentown, NJ) fitted with 20-mm-diameter smooth parallel plates. Temperature was maintained at 25°C using an induction heating device. Samples were sliced into approximately 2-mm-thick disks. Cyanoacrylate glue (Loctite 100, Loctite Corporation) was placed on the bottom and top of the cheese disk to attach it firmly to the plates of the instrument. After attaching the cheese to the stationary bottom plate, the top plate was lowered onto the sample until a normal force of 1N was reached. The sample was then trimmed to fit the plate size and a thin film of synthetic lubricant (Superlube, Loctite Corporation) was applied to prevent moisture loss. Tests were conducted at stresses of 400Pa. This stress was chosen based on a previously conducted stress sweep at 25°C. The creep portion of each test consisted of a stress application for 600s; the stress was removed and strain recovery was measured for an additional 1,200s. Three creep and recovery tests were done for each sample.

Relaxation time was determined by the time required for the delayed strain to reach 63.2% of its final value (Steffe, 1996). Compliance (J), defined as the ratio of strain to stress, was monitored as a function of time. Maximum compliance (Jmax) was the peak compliance reached by the material before the constant stress was removed. Percentage of creep recovery gives an indication of the degree of elasticity in the material and was calculated using the following relationship:

where Jmax is the maximum creep compliance and Jr is compliance after recovery.

Statistical Analysis 

Chemical and sensory data were analyzed using the GLM procedure of SAS (version 8.2, 2001; SAS Institute, Cary, NC). A split-plot model was used with cheese plant (n = 4), replicate (n = 2), and location within 291-kg block (n = 2) as the whole-plot category variables and the linear and quadratic form of time of aging and the interactions of time of cheese aging with plant, replicate, and location as subplot variables, with time as a continuous variable. Replicate was not nested within plant because other analytical factors were included in the replicate variation that were independent of plant. Statistical significance of the individual terms in the model was determined when the F-test for the model was ≤0.05. The interaction of plant×location×replicate was used as the error term to determine whether the effect of plant, replicate, location, or their interaction was significant. The full model error was used to test the significance of the subplot terms. Because time of cheese aging was treated as a continuous variable in the ANOVA model, the linear and quadratic terms for time would be correlated. Distortion of the ANOVA by multicolinearity of these terms in the model was minimized by centering the time of cheese aging data using a mathematical transformation (Glantz and Slinker, 2001). The time variable was transformed as follows: time = month of aging[(last monthfirst month)/2]. This transformation made the data set orthogonal with respect to time of aging. The full model with all interaction terms in the whole and the subplot was run for each parameter (e.g., bitterness, moisture, etc.). After running the full model, terms that were not significant were removed in a stepwise fashion, starting with nonsignificant terms in the subplot with the lowest type III sum of squares. Least squares means were calculated and reported for plant and location within 291-kg block. Analysis of the data with this model allowed us to determine the impact of plant and location (inner vs. outer) within 291-kg block.

Of the 4 cheese plants in the study, 2 used only the stirred-curd method and 2 used only the milled-curd method of Cheddar cheese manufacture. Because no manufacturing plants used both methods, it was not possible to clearly separate the impact of curd type from plant of manufacture because plant and curd type were not independent. Therefore, all statistical analyses were done with manufacturing plant in the model without curd type.

Back to Article Outline

Results 

Cheese Composition 

Cheese manufacturing plant explained the largest portion of variation in the moisture content of the cheese (Table 1), and the cheese from plant 1 had higher moisture than cheeses from the other plants (Table 2). A moisture difference was observed between portions from the exterior and interior (i.e., location within block; Table 1) of the 291-kg block (Table 2). The magnitude of the differences in moisture between inner and outer locations within the 291-kg blocks in our study (Table 2) were much smaller than those reported in previous studies (Reinbold and Ernstrom, 1988; Reinbold et al., 1992; Barbano, 2001).

Table 1. Sums of squares (type III) and probability values (in parentheses) for the ANOVA analysis of the impact of plant, location within 291-kg block (inner vs. outer), and aging time (1, 4, 8, and 12 mo) on moisture, pH, fat, fat dry basis, protein, and salt content of portions cut from 291-kg blocks of Cheddar cheeses during 12 mo at 7°C
FactordfMoisture,%pHFat, % (wet basis)Fat, % (dry basis)Protein, %Salt, %Protein hydrolysis, %
Plant317*
(<0.01)
0.19*
(<0.01)
13.6*
(<0.01)
13.5*
(<0.01)
19.3*
(<0.01)
0.78*
(<0.01)
366*
(<0.01)
Replicate10.01
(0.96)
0.88*
(<0.01)
2.0*
(0.04)
3.6*
(0.02)
0.81*
(0.01)
0.25*
(<0.01)
2*
(<0.01)
Location19*
(<0.01)
<0.01
(0.88)
6.8*
(<0.01)
0.6
(0.39)
0.36
(0.10)
<0.01
(0.95)
76*
(<0.01)
Plant×replicate313*
(<0.01)
0.24*
(<0.01)
6.6*
(<0.01)
14.3*
(<0.01)
11.7*
(<0.01)
0.45*
(<0.01)
11*
(<0.01)
Plant×location3NS1NSNSNSNSNS18*
(0.01)
Plant×replicate×location47.2
(0.08)
0.07
(0.18)
1.5
(0.52)
1.5
(0.73)
2.8*
(0.01)
0.03
(0.74)
4.7
(0.44)
Time10.55
(0.43)
0.11*
(<0.01)
3.0*
(0.01)
2.9*
(0.04)
2.3*
(<0.01)
0.24*
(<0.01)
905*
(<0.01)
Time×plant3NSNSNSNSNSNS18*
(<0.01)
Time×location1NSNSNSNSNSNS13*
(<0.01)
Time×time1NSNSNSNSNSNS65*
(<0.01)
Time×time×replicate1NSNSNSNSNSNSNS
Error21111021.55483161.5132
R2 0.280.500.370.290.690.530.92

1NS = nonsignificant terms removed from the model.

2Degrees of freedom of the error term for percentage of protein hydrolysis = 104.

*P<0.05.

Table 2. Least squares mean values of moisture, pH, fat, fat dry basis, protein, salt content, and percentage of protein hydrolysis of cheese cut from 291-kg blocks of Cheddar cheeses averaged across 12 mo of aging at 7°C
VariableMoisture, %pHFat, % (wet basis)Fat, % (dry basis)Protein,%Salt,%Protein hydrolysis,1%
Plant
Plant 136.5a5.23a35.5c55.9bc23.4b1.73c6.40b
Plant 235.6b5.20a35.7bc55.6c23.6a1.93a9.31a
Plant 335.8b5.24a36.4a56.5a22.7c1.76bc6.57b
Plant 435.6b5.14b36.0b56.0b22.8c1.79b4.48c
Location within 291-kg block2
Inner35.6b5.20a36.1a56.1a23.2a1.80a7.47a
Outer36.1a5.20a35.7b55.9a23.1a1.80a5.91b

a–cMeans in a column (within plant or inner vs. outer location) followed by different letters are different (P<0.05).

1Percentage of protein hydrolysis measured by the o-phthaldialdehyde test.

2Inner = interior of 291-kg block; outer = exterior of 291-kg block.

Temperature-induced moisture migration has been studied extensively. As the temperature of the cheese block cools, the moisture content increases, causing a moisture gradient between the cooler outer locations and warmer inner locations. Barbano (2001) found that the outer portions of 291-kg cheese blocks of reduced-fat Cheddar cheese always displayed higher moisture content (about 5%) than the inner portions and that moisture migration occurred, moving up, down, and laterally within the block. However, the size of the moisture gradient was shown to be much smaller in the top half of a 291-kg block than in the bottom half (Barbano, 2001). Olabi and Barbano (2002) confirmed these findings and determined that moisture migrated from a warmer to cooler temperature both with and against gravity. In our study, the inner locations with the 291-kg blocks did not cool as quickly as the outer locations; consequently, the inner locations had a lower moisture content than the outer locations, and this is consistent with previous reports (Reinbold and Ernstrom, 1988; Reinbold et al., 1992; Barbano, 2001). As expected, there was no change (i.e., the linear and quadratic term for time was not significant; Table 1) in the total moisture content over the aging period. No plant×location within block interaction effect (P>0.05) was detected (Table 1), confirming that the temperature-driven moisture migration within the 291-kg blocks was happening in the same way in all 4 plants, regardless of differences from plant to plant in moisture level or cheese manufacturing practice.

The salt content (Table 2) of the cheese differed from plant to plant, and the plant of manufacture explained most of the variation (i.e., highest type III sum of squares in Table 1) in salt content. However, no difference (P>0.05) in salt content between the inner and outer locations (Tables 1 and 2) was detected.

The cheese pH was lower in plant 4 than in the other 3 plants (Table 2). However, the differences in pH between plant 4 and the other 3 plants were not large; the model explained only 50% of the variation (Table 1). Most of that variation was explained by replicate, and very little of the variation was explained by location within block, despite significant differences in moisture caused by location within the 291-kg block. Cheese pH decreased (P<0.05) with time from an average of 5.23 to 5.15. The downward trend in cheese pH with time for all 4 plants was similar (e.g., no interaction between plant and time). A previous study (Reinbold et al., 1992) found that a pH gradient was present initially in 291-kg blocks, but after cooling for 240h, the pH had equilibrated between the inner and outer portions. Because an initial (prior to cutting the 291-kg blocks) pH measurement was not taken in our study, it is not known whether there were large differences in the pH of the inner and outer locations initially.

The fat content of cheese on a wet basis differed from plant to plant (Tables 1 and 2); this would be driven by plant-to-plant differences in the fat-to-CN ratio in the milk used for cheese making and differences in fat loss in the whey during cheese making. Differences in fat content (on a wet basis only) between the inner and outer locations (inner>outer) were consistently observed throughout aging (Tables 3 and 4). Because fat within the structure of cheese does not move and moisture did move with time, it was not surprising that fat on a wet basis was lower in the outer portions of the cheese (Table 2), where the moisture content was high. On average, the fat content on a wet basis increased (P<0.05) with time but the changes in grand mean values (mean of all cheeses) with time were small and variable (35.6, 36.2, 35.6, and 36.2% at 1, 4, 8, and 12 mo, respectively). Although there was a slight but significant trend for an increase in fat content with time, the variation from month to month was large, and this probably represents inconsistency in composition within the 4 quarters of the upper half of the 291-kg block that could not be distinguished from a time effect in the design of our study. These differences in fat content attributable to location were not detected when fat was expressed on a dry basis (Tables 1 and 2), confirming that differences in wet fat content were due to a dilution effect caused by moisture migration during cooling. The protein content of the cheese differed from plant to plant, but this was driven by differences in moisture content of the cheese and protein and fat content differences from plant to plant in the milk used to make the cheese (Tables 1 and 2).

Table 3. Least squares mean values of moisture, protein, pH 4.6-soluble nitrogen, and 12% TCA-soluble nitrogen, expressed as a percentage of total nitrogen content, of cheese cut from 291-kg blocks of Cheddar cheeses after aging for 12 mo at 7°C
VariableMoisture, %Protein, %pH 4.6-soluble nitrogen, %12% TCA-soluble nitrogen, %
Plant
Plant 136.51a23.22ab25.77b17.57c
Plant 235.04b23.72a28.90a21.51a
Plant 335.32b22.60c28.68a19.08b
Plant 435.73ab22.74bc24.23b14.78d
Location within 291-kg block1
Inner35.32b23.16a27.63a18.65a
Outer35.99a22.98a26.16b17.82b

a–dMeans in a column (within plant or inner vs. outer location within 291-kg block) followed by different letters are different (P<0.05).

1Inner = interior of 291-kg block; outer = exterior of 291-kg block.

Table 4. Sums of squares (type III) and probability values (in parentheses) for the ANOVA analysis of the impact of plant, location within 291-kg block (inner vs. outer), and aging time (1, 4, 8, and 12 mo) on sensory flavor1 of cheese cut from 291-kg blocks of Cheddar cheeses after aging for 12 mo at 7°C
FactordfSensory flavor attribute
WheyDiacetylMilk fat/lactoneSulfurBrothyNuttySourSaltySweetUmami
Plant3186*
(<0.01)
26*
(<0.01)
10*
(<0.02)
37*
(<0.01)
216*
(<0.01)
111*
(<0.01)
3*
(<0.01)
18*
(<0.01)
22*
(<0.01)
28*
(<0.01)
Replicate138*
(<0.01)
1
(0.30)
0.01
(0.87)
44*
(<0.01)
32*
(<0.01)
2*
(0.01)
0.01
(0.87)
0.2
(0.48)
11*
(<0.01)
22*
(<0.01)
Location112*
(<0.01)
4*
(<0.01)
2*
(<0.01)
4*
(0.02)
9*
(0.05)
11*
(<0.01)
0.1
(0.41)
1*
(0.05)
4*
(<0.01)
5*
(<0.01)
Panelist18292*
(<0.01)
107*
(<0.01)
140*
(<0.01)
105*
(<0.01)
206*
(<0.01)
132*
(<0.01)
74*
(<0.01)
38*
(<0.01)
93*
(<0.01)
141*
(<0.01)
Plant×replicate3NS2NSNS5*
(<0.01)
NS4*
(0.03)
2*
(0.02)
2*
(0.02)
NS3*
(0.02)
Location×replicate13*
(0.05)
2*
(<0.01)
2*
(<0.01)
NS2
(0.17)
NSNSNSNSNS
Plant×location3NS2*
(0.02)
NSNSNSNSNSNSNSNS
Plant×location×replicate98
(0.38)
3*
(0.04)
6*
(0.05)
3
(0.40)
14 (0.23)1
(0.89)
1
(0.80)
2
(0.17)
3
(0.50)
1
(0.58)
Time1725*
(<0.01)
71*
(<0.01)
9*
(<0.01)
496*
(<0.01)
875*
(<0.01)
148*
(<0.01)
2*
(<0.01)
13*
(<0.01)
131*
(<0.01)
292*
(<0.01)
Time×plant331*
(<0.01)
5*
(<0.01)
5*
(<0.01)
4*
(0.02)
10*
(0.04)
14*
(<0.01)
3*
(<0.01)
3*
(<0.01)
2*
(0.05)
NS
Time×location1NSNSNSNSNS2*
(0.02)
NSNS2*
(<0.01)
NS
Time×replicate118*
(<0.01)
3*
(<0.01)
NSNSNSNS6*
(<0.01)
5*
(<0.01)
6*
(<0.01)
1*
(0.04)
Time×plant×location3NS3*
(0.02)
NSNSNSNSNSNSNSNS
Time×plant×replicate3NSNSNSNSNS9*
(<0.01)
5*
(<0.01)
NSNSNS
Time×time14*
(0.02)
7*
(<0.01)
3*
(<0.01)
3*
(<0.01)
43*
(<0.01)
10*
(<0.01)
3*
(<0.01)
NSNSNS
Time×time×plant322*
(<0.01)
NSNS3
(0.09)
NS8*
(<0.01)
2
(0.06)
NSNSNS
Time×time×replicate133*
(<0.01)
2*
(<0.01)
NS37*
(<0.01)
52*
(<0.01)
7*
(<0.01)
NSNSNSNS
Time×time×location1NSNSNSNSNSNS1*
(0.05)
NSNSNS
Error3 1,2153625016401,842636337340429395
R2 0.570.420.270.570.490.480.250.180.380.56

1Fruity, FFA, catty, and bitter tastes were not detected in cheeses. Cooked flavor is not shown because there were no significant differences among the cheeses.

2NS = nonsignificant terms removed from the model.

3Degrees of freedom in the error term varied from 1,484 to 1,497, depending on the nonsignificant terms removed from the model for each sensory flavor attribute.

*P<0.05.

Degree of Protein Hydrolysis 

There were large differences in proteolysis during aging among the different plants (Tables 1 and 2), and the statistical model explained a large amount (r2 = 0.92) of the variation in the data. Cheese from plant 2 had more than twice as much proteolysis, averaged across 12 mo, as plant 4 (Table 2) when measured by the OPA test. More proteolysis occurred in the inner locations than in the outer locations within the 291-kg blocks. There was also a significant increase in the percentage of protein hydrolysis for all plants with aging time (Figure 2). There was a significant linear and quadratic component (Table 1) to this age-dependent increase, and there was a time×plant and time×location within block interaction. The cheeses with lower moisture (inner location) had a higher degree of hydrolysis. This was related to the temperature gradient during cooling, with lower moisture inner locations remaining at a higher temperature for a longer time than the outer locations.

  • View full-size image.
  • Figure 2. 

    Change in the percentage of protein hydrolysis of inner and outer locations within 291-kg blocks of Cheddar cheese during 12 mo of aging at 7°C (♦ = inner location; ■ = outer location).

The pH 4.6 acetate-soluble nitrogen and 12% TCA-soluble nitrogen content of the cheeses expressed as a percentage of total nitrogen are shown in Table 3. These analyses were conducted in a different laboratory from the laboratory where the OPA testing was done. The differences among plants and treatments for moisture and protein tests on the 12-mo samples analyzed at Cornell (Table 3) were consistent with those in the North Carolina State laboratory (Table 2). Both the pH 4.6 and 12% TCA-soluble nitrogen were different among plants and locations within block. There was more proteolysis in the inner than the outer location across all plants, but the influence of the plant on overall proteolysis was much larger than location within block (Table 3). During the first 30 to 60 d of aging, most of the proteolysis is done by the residual chymosin in the cheese (Grappin et al., 1985). Cheeses from plants 2 and 3 had a higher rate of rennet addition, a higher cheese temperature when the cheese was filled into the forms, and also the highest soluble nitrogen content at 12 mo. The cheese in the internal portion of the 291-kg block stays warmer longer and this may allow more degradation of intact CN by chymosin. The TCA-soluble nitrogen is more selective for lower molecular weight proteolysis products and reflects the completeness or depth of proteolysis caused by the action of proteases and peptides of lysed starter and nonstarter bacteria and their extracellular enzymes, which is sometimes called secondary proteolysis (Rank et al., 1985). The pH 4.6 and 12% TCA-soluble nitrogen were lowest in cheese from plant 4, and this may be due to the combination of a low level of rennet use, a high level of milk fortification, and a lower cheese temperature when the 290-kg blocks were filled.

Flavor Analysis 

Descriptive Sensory Analysis of Flavor 

Cheeses were characterized by high intensities of diacetyl, cooked, whey, and milk fat/lactone flavors at 1 mo. These flavors are commonly found in young Cheddars (Drake et al., 2001). There were significant manufacturing plant effects for all sensory attributes listed in Table 4, and manufacturing site clearly played a key role in the flavor profiles of cheese. No influence (P>0.05) of the factors in our study on cooked flavor was detected (Table 5). Whey and diacetyl flavors decreased over time for all treatments (Figures 3 and 4), as did milk fat flavor (data not shown). There was an interaction (P<0.05) of plant by both the linear and quadratic terms for time (Table 4) for whey flavor; this can be seen in Figure 3, with whey flavor not decreasing as rapidly in cheese from plant 4. Plant-to-plant differences in whey flavor (Figure 3) may be related to differences in the milk standardization and fortification strategies used in the different cheese plants prior to cheese making. Those plants using a higher level of milk solids fortification (e.g., plant 4) may maintain a higher whey flavor intensity longer. Whey flavor intensity was slightly higher in the outer locations of 291-kg blocks (Table 5), which is consistent with this being a zone of higher moisture content.

Table 5. Least squares mean values of sensory flavor of cheese cut from 291-kg blocks of Cheddar cheeses during 12 mo of aging at 7°C by plant and inner vs. outer location1
VariableSensory flavor attribute
CookedWheyDiacetylMilk fat/lactoneSulfurBrothyNuttySourSaltySweetUmami
Plant
Plant 12.38a1.92b0.46b2.52b1.33b1.84c0.50c3.37ab3.66c1.79c1.38c
Plant 22.37a1.46d0.31c2.68a1.43a2.37a1.23a3.42a3.93a2.06a1.62a
Plant 32.43a1.73c0.35c2.55b1.44a2.00b0.76b3.33b3.74b1.90b1.49b
Plant 42.46a2.59a0.66a2.70a0.96c1.31d0.35d3.42a3.70bc1.75c1.25d
Location within 291-kg block2
Inner2.37a1.84b0.40b2.57a1.34a1.95a0.79a3.40a3.78a1.92a1.49a
Outer2.40a2.02a0.50a2.65a1.24b1.80b0.63b3.36a3.73a1.83b1.38b

a–dMeans in a column (within plant, milled vs. stirred curd, or inner vs. outer location) followed by different letters are different (P<0.05).

1Fruity, FFA, catty, and bitter taste were not detected in cheeses.

2Inner = interior of 291-kg block; outer = exterior of 291-kg block.

  • View full-size image.
  • Figure 3. 

    Change in whey flavor from a 291-kg block of Cheddar cheese aged 12 mo at 7°C. Data were averaged across block location because this parameter had no significant effect for this attribute (♦ = plant 1; ■ = plant 2; ▴ = plant 3;×= plant 4).

  • View full-size image.
  • Figure 4. 

    Change in diacetyl flavor from a 291-kg block of Cheddar cheese aged 12 mo at 7°C. Data were averaged across block location (♦ = plant 1; ■ = plant 2; ▴ = plant 3;×= plant 4).

Diacetyl flavor intensity was influenced by location within block (Table 4), with higher diacetyl flavor in the outer locations within the 291-kg blocks (Table 5). There was a significant plant×time interaction (Table 4); this can be seen in Figure 4, with the difference in diacetyl flavor among plants being larger initially and then becoming more similar among plants by 12 mo of age. These are young, undeveloped flavors and are commonly found in Cheddar cheese <4 mo old. With aging time, the intensities of these flavors (i.e., whey, diacetyl, and milk fat/lactone) decreased, consistent with previous studies (Drake et al., 2001; Rehman et al., 2003).

Fruity, FFA, or catty flavors were not detected in any of the cheeses. There was an effect (P<0.05) of the interaction of the linear and quadratic terms for time×plant on nutty flavor (Table 4). Nutty flavor was more intense in the inner than the outer location within 291-kg blocks (Table 5). High nutty flavor intensities are usually observed only in Cheddar cheeses>8 mo old (Avsar et al., 2004). The inner locations (Table 5) had higher intensities of aged flavors (nutty, brothy, and sulfur) and the outer locations had higher intensities of young, undeveloped flavors (whey and diacetyl).

There was a significant effect (P<0.05) of plant and time of aging (Table 4) on all the basic tastes, with the rate of increase with aging time differing among plants (data not shown). Bitter taste was not detected in these cheeses. The inner locations were sweeter than the outer locations (Table 5). There were also significant differences in umami intensity between the inner and outer locations within 291-kg blocks. The inner locations had higher umami intensities than the outer locations. Umami is typically higher in aged Cheddar cheeses (Drake et al., 2001), and umami was positively correlated with the aged flavors sulfur, brothy, and nutty (r2>0.90) in this study (data not shown). There was no difference in salty flavor due to location within the 291-kg block, which is consistent with the fact that no differences in salt content were detected among block locations (Tables 2 and 5). In general, even though there were differences in basic tastes among the cheeses, the differences may not be of practical significance, because differences among the means were very small and variation caused by panelists effects (Table 4) accounted for most of the explained variation in the basic tastes.

Instrumental Volatile Analysis 

Ten volatile compounds were quantified. Two different groups of flavor compounds were clearly identified: compounds that contribute mainly to young, undeveloped flavors (cooked, whey, diacetyl, and milk fat/lactone flavors) and compounds that contribute to aged, developed flavors (sulfur, brothy, and nutty flavors; Tables 6 and 7). 2,3-Butanedione, 3-hydroxy-2-butanone, and acetic acid all contribute to young flavors (Singh et al., 2003), whereas esters and aldehydes contribute more to aged flavors (Singh et al., 2003). In contrast to sensory perception of flavor, production facility did not a have strong direct influence on the level of volatile compounds, but manufacturing plant was significant in various interactions terms in combination with location within block and time for many of the volatile compounds (Table 6). 2,3-Butanedione (diacetyl), 3-hydroxy-2-butanone (acetoin), acetic acid, and butanoic acid were detected at high concentrations at 1 mo, after which the concentrations decreased (Table 6) sharply with time (Figure 5), as confirmed by the fact that the quadratic term for time was significant for 2,3-butanedione (diacetyl), 3-hydroxy-2-butanone (acetoin), and acetic acid (Table 6). The decrease in butanoic acid with time was more linear (Table 6 and Figure 5). 2,3-Butanedione (diacetyl) and 3-hydroxy-2-butanone (acetoin) likely contributed to the buttery/diacetyl flavors (Singh et al., 2003) present in the cheeses after 1 mo of aging. The acetoin (3-hydroxy-2-butanone) content of the cheeses was highly variable and differed greatly among plants at 1 mo (range from 342 to 2,358μg/kg) but by 8 and 12 mo the levels among plants were virtually identical (Table 6, time×time and time×plant interactions). Acetic acid concentrations in cheese were highly variable and differences were mostly due to complex interaction effects among time, plant, and location within block (Table 6). Both acetic and butanoic acids have previously been identified in mild Cheddar cheese (Milo and Reineccius, 1997). Variation in butanoic acid content of the cheeses was mostly explained by interactions of time, plant, and block location (Table 6), with butanoic acid decreasing with time of aging (Figure 5). All these compounds (2,3-butanedione, 3-hydroxy-2-butanone, acetic acid, and butanoic acid) were present above the sensory threshold at 1 mo. Milo and Reineccius (1997) reported these compounds to be key aroma-active compounds in mild Cheddar cheese.

Table 6. Sums of the squares (type III) and probability values (in parentheses) for the ANOVA analysis of the impact of plant, location within 291-kg block (inner vs. outer), and aging time (1, 4, 8, and 12 mo) on the concentration of aroma-active compounds identified in cheese cut from 291-kg blocks of Cheddar cheeses during 12 mo of aging at 7°C
FactordfAroma-active compound
2,3-Butanedione (diacetyl)3-Methyl butanalAcetic acid3-Hydroxy-2-butanone (acetoin)Methyl butanoateEthyl butanoateButanoic acidHexanalMethyl hexanoateEthyl hexanoate
Plant31.1×107
(0.82)
50
(0.97)
1.1×1011
(0.99)
3.2×106
(0.63)
20
(0.50)
1.2×108
(<0.10)
1.4×107
(0.25)
1.21×105
(0.79)
3.4×106
(0.99)
6.0×104
(0.48)
Replicate14.5×106
(0.56)
71
(0.61)
2.6×1011
(0.69)
2.7×105
(0.71)
6
(0.41)
1.3×108
(0.01)
1.6×106
(0.46)
1.45×105*
(0.29)
6.4×107
(0.41)
1.7×105*
(0.02)
Location11.4×107
(0.31)
18
(0.79)
6.4×1010
(0.84)
8.0×105
(0.52)
55
(0.06)
2.1×107
(0.23)
2.26×106
(0.36)
1.3×103
(0.92)
9.0×104
(0.97)
4.0×103
(0.90)
Plant×replicate3NS1753*
(<0.01)
3.5×1012
(<0.01)
NS388*
(<0.01)
8.0×107*
(0.03)
1.5×107*
(0.02)
NS3.3×108*NS
Plant×location35.3×107*
(0.05)
740*
(<0.01)
6.9×1012*
(<0.01)
1.0×107*
(<0.01)
203*
(<0.01)
NS1.3×107*
(0.04)
2.7×105*
(0.01)
1.9×108*
(<0.01)
NS
Replicate×location12.2×107
(0.07)
179*
(0.01)
NSNS75*
(<0.01)
NSNS1.3×105*
(0.02)
6.4×107*
(<0.01)
NS
Plant×replicate×location3 to 107.0×107
(0.11)
660*
(<0.01)
5.7×1012*
(<0.01)
1.2×107*
(<0.01)
20
(0.42)
8.6×107
(0.19)
9.4×106
(0.20)
6.5×105*
(<0.01)
2.1×108*
(<0.01)
2.3×105
(0.59)
Time13.0×108*
(<0.01)
644*
(<0.01)
1.3×1013*
(<0.01)
3.7×107*
(<0.01)
1,365*
(<0.01)
NS4.0×107*
(<0.01)
NS3.0×108*
(<0.01)
1.3×105*
(0.03)
Time×plant3NS1,264*
(<0.01)
6.5×1012*
(<0.01)
1.1×107*
(<0.01)
704*
(<0.01)
NS1.5×107*
(0.02)
NS3.5×108*
(<0.01)
NS
Time×location1NS426*
(<0.01)
2.9×1012*
(<0.01)
NS73*
(<0.01)
NSNSNS8.3×107*
(<0.01)
NS
Time×replicate11×108*
(<0.01)
590*
(<0.01)
7.6×1012*
(<0.01)
NS363*
(<0.01)
NS2.0×107*
(<0.01)
1.5×105*
(0.03)
6.7×107*
(<0.01)
NS
Time×plant×location3 or 4NS1,188*
(<0 01)
8.1×1012*
(<0.01)
1.8×107*
(<0.01)
308*
(<0.01)
NS1.9×107*
(0.02)
NS4.1×108*
(<0.01)
NS
Time×plant×replicate3NS1,054*
(<0.01)
4.3×1012*
(<0.01)
NS519*
(<0.01)
NS1.8×107*
(<0.01)
NS3.3×108*
(<0.01)
NS
Time×time11.8×108*
(<0.01)
82
(0.08)
4.7×1012*
(<0.01)
9.2×106*
(<0.01)
387*
(<0.01)
8.3×107*
(<0.01)
NSNS7.2×107*
(<0.01)
6.3×105*
(<0.01)
Time×time×plant3NS451*
(<0.01)
2.9×1012*
(<0.01)
1.1×107*
(<0.01)
290*
(<0.01)
7.7×107*
(0.03)
NS2.3×105*
(0.03)
1.0×108*
(<0.01)
NS
Time×time×location1NS222*
(<0.01)
1.8×1012*
(<0.01)
NSNSNSNSNS4.1×107*
(<0.01)
NS
Time×time×replicate1 or 25.7×107*
(<0.01)
83
(0.09)
5.0×1012*
(<0.01)
NS101*
(<0.01)
8.3×107*
(<0.01)
1.6×107*
(<0.01)
2.4×105*
(<0.01)
NSNS
Error2 131126124131130127133130124144
R2 0.500.750.760.640.870.260.500.350.770.25

1NS = not significant term removed from the model.

2Degrees of freedom in the error term varied from 124 to 144, depending on the number of nonsignificant terms removed from the model for each aroma-active compound.

*P<0.05.

Table 7. Least squares mean values of volatile concentrations (μg/kg) of aroma-active compounds identified in cheese cut from 291-kg blocks of Cheddar cheeses during 12 mo of aging at 7°C
VariableAroma-active compound
2,3-Butanedione (diacetyl)3-Methyl butanalAcetic acid3-Hydroxy-2 butanone (acetoin)Methyl butanoateEthyl butanoateButanoic acidHexanalMethyl hexanoateEthyl hexanoate
Retention index on DB-5 MS column6686867007187237768508539961,000
Threshold,1μg/kg91.536,0002,900314.7860181.790
Plant
Plant 1860a1.6a28,000a393a1.8b1,862a74a158a336a132a
Plant 21,519a1.3a620,000a355a1.6b379a579a116a99a89a
Plant 31,151a6.1a120,000a291a5.1a716a437a140a3,055a128a
Plant 41,520a1.3a160,000a686a1.3b366a938a129a933a94a
Location within 291-kg block2
Inner946a3.6a110,000a358a3.0a1,220a378a132a1,776a109a
Outer1,579a1.5a350,000a505a1.9b442a636a139a436a112a

a,bMeans in a column (within plant or inner vs. outer location) followed by different letters are different (P<0.05).

1Thresholds reported orthonasally in water by Rychlik et al. (1998).

2Inner = interior of 291-kg block; outer = exterior of 291-kg block.

  • View full-size image.
  • Figure 5. 

    Change in volatile concentration from a 291-kg block of Cheddar cheese aged 12 mo at 7°C. There were no significant differences in concentration between inner and outer locations (P>0.05; ♦ = acetic acid; ■ = 2,3-butandione (diacetyl); ▴ = 3-hydroxy-2-butanone (acetoin);×= butanoic acid).

Prior to 4 mo of aging, 3-methyl butanal was not detected (data not shown). 3-Methyl butanal is one of the compounds that causes a nutty flavor in Cheddar cheese (Avsar et al., 2004). Nutty flavor was higher in the inner location within 291-kg blocks (Table 5), but 3-methyl butanal concentrations were not detected to be higher (P>0.05) in this location within the 291-kg blocks (Tables 6 and 7).

Prior to 8 mo of aging, the esters methyl hexanoate and methyl butanoate were not detected in the cheeses. The concentration of methyl hexanoate and methyl butanoate increased rapidly between 8 and 12 mo of aging, as seen by the significant effects of linear and quadratic terms for time (Table 6). Methyl butanoate concentration was higher in the inner than outer locations within the 291-kg blocks (Table 7), whereas no differences in methyl hexanoate (Table 7) could be detected because of high variation in the data. The esters were likely formed from short- and medium-chain FFA (Singh et al., 2003) and are also not usually present in high concentrations in Cheddars <8 mo old. No differences were detected in the concentrations of ethyl butanoate or ethyl hexanoate. The compounds that contributed to aged flavor (esters and aldehydes) were present in the highest concentrations after 12 mo of aging (data not shown). Because of the high degree of variation in the instrumental volatile analysis data, it was difficult to detect any significant differences in these compounds attributable to location within block. However, the significant time-dependent changes in young and aged cheese flavor compounds (Table 6) generally supported the time-dependent changes in descriptive sensory flavor analysis results (Table 4).

Texture Analysis 

Descriptive Sensory Analysis 

When sensory texture attributes of the cheeses were plotted as a function of time of aging, most attributes increased slightly in value from 1 to 4 mo of aging and then remained the same after 4 mo (data not shown). Two exceptions to this were hand springiness and hand rate of recovery. Both these attributes demonstrated large decreases (P<0.05) in value with age that were predominately linear in character for all treatments (Table 8). The inner locations within 291-kg blocks had higher hand firmness (Tables 8 and 9) than outer locations, which is consistent with the lower moisture content of the inner location. There was a significant time×block location interaction (Table 8), with the magnitude of the difference in hand firmness between the inner and outer locations decreasing with time of aging. The outer locations were more cohesive, had a smoother mass, and had more residual mouth coating than inner locations with 291-kg blocks (Table 9). Initially, there was no difference in cohesiveness between the inner and outer locations, but cohesiveness for the inner and outer locations changed with aging (time×location within block interaction), with the outer locations developing more cohesiveness (Figure 6). Fracturability changed very little with time of aging. No differences in mouth firmness attributable to plant, location within block, or time were detected. After 8 mo, all cheeses were more cohesive, were more adhesive, and had a smoother mouth coating (data not shown) because of the significant impact of time (Table 8).

Table 8. Sums of squares (type III) and probability values (in parentheses) for the ANOVA analysis of the impact of plant, location within 291-kg block (inner vs. outer), and aging time (1, 4, 8, and 12 mo) on sensory texture of cheese cut from 291-kg blocks of Cheddar cheeses during 12 mo of aging at 7°C
FactordfSensory texture attribute
Hand firmnessHand springinessHand rate of recoveryMouth firmnessFracturabilityBreakdownCohesivenessAdhesivenessMouth coatingSmoothness of mass
Plant316
(0.12)
420*
(0.01)
333*
(<0.01)
NS115
(0.77)
622
(0.73)
35
(0.14)
154*
(<0.01)
78*
(<0.01)
134*
(0.02)
Replicate18
(0.07)
899*
(<0.01)
495*
(<0.01)
NS50
(0.10)
641
(0.27)
34*
(0.03)
0.1
(0.89)
0.1
(0.88)
5
(0.27)
Location1109*
(<0.01)
107*
(0.03)
78*
(0.03)
NS29
(0.18)
638
(0.27)
26*
(0.05)
2
(0.46)
29*
(<0.01)
34*
(0.04)
Panelist9302*
(<0.01)
713*
(<0.01)
303*
(0.03)
NS153*
(<0.01)
2,730
(0.61)
156*
(<0.01)
202*
(<0.01)
152*
(<0.01)
146*
(<0.01)
Plant×replicate3121*
(<0.01)
454*
(<0.01)
364*
(<0.01)
NSNSNS159*
(<0.01)
23*
(0.02)
70*
(<0.01)
74*
(<0.01)
Plant×location346*
(<0.01)
48
(0.09)
NSNS10
(0.07)
NSNSNSNS15
(0.06)
Plant×replicate×location4 to 106
(0.56)
41
(0.22)
58
(0.08)
NS94*
(<0.01)
4,700
(0.25)
32
(0.13)
23
(0.18)
10
(0.39)
8
(0.24)
Time1NS2,830*
(<0.01)
1,481*
(<0.01)
NS46*
(<0.01)
NS506*
(<0.01)
647*
(<0.01)
136*
(<0.01)
263*
(<0.01)
Time×plant325*
(<0.01)
NSNSNS37*
(<0.01)
NS36*
(<0.01)
17*
(0.05)
NS22
(0.01)
Time×location116*
(<0.01)
NSNSNSNSNS15*
(0.02)
21*
(<0.01)
NSNS
Time×replicate1NS663*
(<0.01)
529*
(<0.01)
NSNSNS151*
(<0.01)
253*
(<0.01)
110*
(<0.01)
192*
(<0.01)
Time×plant×location3NSNSNSNSNSNS27*
(0.02)
19*
(0.04)
20*
(<0.01)
20*
(0.04)
Time×plant×replicate3 to 638*
(<0.01)
NS100.0*
(<0.01)
NS25*
(<0.01)
NS35*
(<0.01)
33*
(<0.01)
14*
(0.02)
35*
(<0.01)
Time×time123*
(<0.01)
237*
(<0.01)
122*
(<0.01)
NSNSNS255*
(<0.01)
373*
(<0.01)
70*
(<0.01)
186*
(<0.01)
Time×time×plant351*
(<0.01)
NSNSNS10
(0.06)
NS28*
(0.02)
NSNSNS
Time×time×replicate1NS235*
(<0.01)
91*
(<0.01)
NS6*
(0.04)
NSNS51*
(<0.01)
14*
(<0.01)
45
(<0.01)
Error2 1,8757,1404,961668,8301,364373,6942,7722,2121,3611,947
R2 0.300.490.460.020.280.020.390.390.390.41

1NS = nonsignificant terms removed from the model.

2Degrees of freedom in the error term varied from 975 to 993, depending on the number of nonsignificant terms removed from the model for each sensory texture attribute.

*P<0.05.

Table 9. Least squares mean values of sensory texture of cheese cut from 291-kg blocks of Cheddar cheeses during 12 mo of aging at 7°C
VariableSensory texture attribute
Hand firmnessHand springinessHand rate of recoveryMouth firmnessFracturabilityBreakdownCohesivenessAdhesivenessSmoothness of massMouth coating
Plant
Plant 110.2a7.0a5.8ab8.8a5.7a9.6a9.4a8.5b8.6bc9.7b
Plant 210.1a5.6b4.5c11.85a6.0a9.9a9.7a9.1a8.8b9.7b
Plant 39.7a7.0a5.5b8.1a5.2a10.2a10.3a9.2a9.4a10.3a
Plant 410.3a7.3a6.0a8.8a5.7a11.7a9.0a8.3b8.4c9.5b
Location within 291-kg block1
Inner10.4a6.4b5.2b10.5a5.8a9.6a9.5b8.7a8.6b9.6b
Outer9.8b7.1a5.7a8.3a5.5a11.1a9.8a8.8a9.0a10.0a

a–cMeans in a column (within plant or inner vs. outer location) followed by different letters are different (P<0.05).

1Inner = interior of 291-kg block; outer = exterior of 291-kg block.

  • View full-size image.
  • Figure 6. 

    Change in cohesiveness measured by sensory texture in inner and outer locations of Cheddar cheese removed from a 291-kg block aged 12 mo at 7°C (♦ = inner location; ■ = outer location).

Rheological Analysis 

Large Strain Analysis 

There was no effect of production facility on true shear stress, strain, or modulus (data not shown). Fracture stress did not change with time of aging, whereas strain decreased from 0.55 to 0.23 and modulus increased from 69 to 216 kPa with time. In this study, as modulus increased, hand springiness and hand rate of recovery decreased. Previous research has shown that the fracture modulus is negatively correlated with hand springiness and hand rate of recovery for young cheeses (Brown et al., 2003). Location within the 291-kg block had an impact on fracture strain, with the least squares means for inner and outer locations being 0.43 and 0.32, respectively. Fracture strain decreased with time of aging and was higher for the outer location (0.43) than the inner location (0.32) within the 291-kg block. Our data are consistent with those of Brown et al. (2003), who found that fracture strain was highly correlated with firmness and fracturability and similar trends were observed (r2>0.8).

Creep and Recovery Analysis 

No significant effects (P>0.05) were observed among plants (data not shown) on measures of creep compliance or recovery. Significant time-dependent decreases from 1 to 12 mo of aging in percentage of creep recovery (from 49.2 to 34.8%), relaxation time from 128 to 79s, Jmax (from 0.0035 to 0.0023Pa−1), and γrec (from 0.0025 to 0.0013Pa−1). This is expected because during aging, there is an increase in proteolysis, a decrease in hand springiness and hand rate of recovery, and an increase in fracturability. This indicates that the cheese network is breaking down to some extent (as expected during aging). As the cheese network breaks down, the elastic elements of the cheese decrease, and the recovery declines as well. The percentage of creep recovery and Jr are indicators of the visco-elastic properties of cheese and the ability of the cheese to recover from a deformation (Brown et al., 2003). Therefore, during aging as the cheese networks break down, the elastic elements decrease and the cheese becomes more viscous, which is reflected in the decrease of Jr and creep recovery. There was an effect (P<0.05) of location within block, with the values being higher for the inner location for percentage of creep recovery (44 vs. 40%), Jr (5.26×10−6 vs. 3.67×10−6 Pa−1), and Jmax (0.0033 vs. 0.0024Pa−1).

Back to Article Outline

Discussion 

Cheese Manufacturing Conditions 

As would be expected, some important cheese manufacturing procedures and conditions differed from one commercial cheese plant to another (listed in the materials and methods section), and that may be why the term for plant in the ANOVA model explained a large amount of the variation in flavor scores (Table 4). Those conditions that might have influenced flavor and texture development are rennet type and use rate, addition of enzymes to increase flavor development during aging, temperature of the curd going into the 291-kg mold, fortification of milk for cheese making with nonfat milk solids, and cooling conditions after pressing. Those differences in conditions specific to the production facility would be expected to influence the overall flavor and texture development during aging.

Flavor Development 

The main flavor reactions in Cheddar cheese are lipolysis and proteolysis. Flavor development in Cheddar can be greatly affected by the conditions during cheese making, the aging conditions, and the initial composition of the cheese as well as how these composition parameters change throughout aging (Banks et al., 1995; Chen et al., 1998). When looking solely at moisture content, it was unexpected that the inner locations within 291-kg blocks generally had more intense aged Cheddar cheese flavors given that the moisture content of these locations within 291-kg blocks was lower than at the outer locations (Table 5). The salt-to-moisture ratio can also influence proteolysis, but there was very little difference in the salt-to-moisture ratio in inner vs. outer locations within block in our study. During the first few days after manufacture, the salt-to-moisture ratio is the main influence controlling the water activity of the cheese (Lawrence et al., 2004). The water activity influences bacterial growth and enzyme activity, which are the main drivers for proteolysis (Lawrence et al., 2004), with all other factors being equal. Differences in fat content were noted (inner>outer), but the differences were very small (Table 2). However, even though there were differences in the fat content between inner and outer locations, proteolysis was likely the main reaction controlling flavor development in these Cheddar cheeses during 12 mo of aging. Proteolysis is crucial in the formation of flavor compounds. The catabolism of AA results in many different aroma-active compounds that have previously been shown to be important in Cheddar cheese flavor (Christensen and Reineccius, 1995; Milo and Reineccius, 1997; Singh et al., 2003). About half of all flavor compounds in cheese result from the breakdown of AA, especially Met and Leu (Yvon and Rijnen, 2001).

As mentioned, the inner location within the 291-kg blocks displayed higher percentages of protein hydrolysis, pH 4.6 acetate buffer, and TCA-soluble nitrogen, even though these cheeses had lower moisture. This clearly affected aged Cheddar flavor formation in the inner and outer locations (Tables 4 and 5) within the 291-kg blocks in this study. The outer locations not only had higher intensities of the young, undeveloped flavors cooked and diacetyl, but also had lower intensities of the aged, developed flavors sulfur and brothy. These results suggest that aged Cheddar flavor developed more rapidly in the inner locations (more protein hydrolysis and less moisture) than in the outer locations (less protein hydrolysis and more moisture). A direct example of this was nutty flavor formation. This flavor is formed by the aldehydes 2/3-methyl butanal and 3-methyl propanal, which are formed from the Strecker degradation of Leu (Singh et al., 2003; Avsar et al., 2004). Prior to 4 mo of aging, 3-methyl butanal was not detected. It takes several months for this compound to be formed at concentrations above sensory thresholds, which is why nutty flavor is not commonly observed in Cheddar cheeses <8 mo old (Avsar et al., 2004). The presence of this compound was correlated with nutty flavor in this study as well (r2>0.65, P<0.05).

What caused more proteolysis and flavor development in the inner locations of the 291-kg blocks? The combination of a higher rennet use rate and higher temperature of the cheese going into the 291-kg blocks (i.e., longer time to cool) for plants 2 and 3 caused more proteolysis during aging (Table 1, plant×time interaction; Figure 2) and more aged cheese flavor development (Tables 4 and 5). The lowest proteolysis in plant 4 was caused by the combination of the lowest rennet use rate per unit of milk, the high level of milk fortification (i.e., even less rennet per unit of CN), the lower temperature of curd going into the hoop, and the low cooling temperature. Cheese produced in plant 2 had the highest level of proteolysis (Tables 2 and 3) and also had the highest nutty and brothy flavors (Table 5). Plant 4 had the highest intensity of young cheese flavors (Table 5), and this was probably related to the high level of fortification of the milk and low rennet use rate. The time it takes for the inner portion to cool is much longer compared with the outer portion of the block (Reinbold et al., 1992). Depending on the cooling conditions, it can take as long as 12 to 14 d for the interior portion of 291-kg blocks to cool completely, depending on the material used for the cheese hoop (Reinbold et al., 1992). Temperature has a large influence on the rate of proteolysis, and it is likely that chymosin produced more primary proteolysis in the inner portion of the block because of the higher temperature for a longer time than in the outer portion of the 291-kg blocks. This is supported by the higher pH 4.6 acetate buffer soluble nitrogen (which is increased mostly by chymosin activity) observed in the inner location within the 291-kg blocks (Table 3). Also, the higher temperature in the interior portion of the block may favor the growth of a different mixture of nonstarter lactic acid bacteria compared with exterior locations, and this could influence peptidase activity (higher 12% TCA-soluble nitrogen) and the development of other volatile flavor compounds.

Texture Development 

Although proteolysis was likely the driver of differences in flavor formation among different locations within the 291-kg block of cheese, other factors appeared to influence texture. The pH of the cheese influences CN solubility (Lawrence et al., 2004). At pH 5.35, CN hydration is maximized (Lawrence et al., 2004), which would result in a softer texture. However, in our study, the pH was very similar for inner and outer locations within 291-kg blocks (Table 2), and differences in pH did not appear to influence texture. In our study, moisture was likely the main influence in texture differences between locations within 291-kg blocks. The relative ratios of fat, protein, and moisture affect the texture (Hort and Le Grys, 2001). Increased moisture contributed to a softer texture (outer locations), whereas lower moisture contributed to more brittle cheeses (inner locations). These properties were evident by both sensory and rheological measurements. Higher moisture can weaken the protein network because the volume fraction of the protein is lower when more water is present (Lawrence et al., 2004). Increased brittleness is consistent with the fracture strain being higher (P<0.05) in the outer location than in the inner location (0.43 vs. 0.32).

From the proximate analysis results, one can deduce that the moisture gradient had an impact on texture. The inner locations within 291-kg blocks were consistently more firm and fracturable. This result was not expected because the inner locations within 291-kg blocks had a higher degree of protein hydrolysis, which should have made the texture less firm than the outer locations, which had less proteolysis (Lucey et al., 2003). Extensive proteolysis may severely break down the CN structure, resulting in a less firm texture (Lucey et al., 2003); however, in all the cheeses in the current study, hand firmness was almost constant during 12 mo of aging, whereas there were large increases in proteolysis. Therefore, proteolysis may be less important with respect to changes in firmness than for some other texture attributes, such as hand springiness and hand rate of recovery, which decrease greatly with increasing proteolysis.

The higher moisture content of the outer locations within 291-kg blocks also influenced the sensory texture attributes. A higher moisture content is associated with a less rigid cheese matrix (Lucey et al., 2003). The outer locations within 291-kg blocks were more springy and cohesive, and had a more smooth mass in the mouth and mouth coating than the inner locations (Table 9). The inner locations were firmer (hand only; Table 9). Firmness, mouth fracturability, and hand fracturability were highly correlated with each other and were negatively correlated with the breakdown terms (breakdown, cohesiveness, adhesiveness). These findings were in agreement with those observed by Brown et al. (2003) for Mozzarella cheese using the same texture lexicon.

Back to Article Outline

Conclusions 

Large differences were observed in both the cheese flavor and volatile profiles among cheese plants. These differences among plants were caused in part by differences in rennet use rate, time-temperature relationships during cooling, and the level of milk fortification. Differences also were observed between the inner and outer locations within the 291-kg blocks. The inner locations developed aged flavors more rapidly and intensely than the outer locations. A moisture gradient was present (outer>inner), and this undoubtedly influenced flavor and texture. However, the degree of protein hydrolysis was greater in the inner locations than the outer locations. The inner locations within the 291-kg blocks developed more aged flavor. It is likely that the higher temperature of the cheese in the inner location within the 291-kg block for a longer time during cooling caused more aged cheese flavor to develop in the inner location of the blocks. There were also differences in the inner and outer locations for sensory texture. The inner locations were firmer and the outer locations were more cohesive in the mouth. The moisture gradient appeared to have more impact on texture than on flavor in the inner vs. outer locations within the 291-kg blocks. The differences in flavor and texture development observed between inner and outer locations from different 291-kg blocks could be very important when determining end product use. Industrially, cheese can have many applications, such as food ingredients, cheese shreds or slices, or further aging prior to shredding, slicing, or direct consumption. The inner locations, with a firmer texture and more flavor, could be used in some applications, whereas the softer, less flavorful outer blocks may be better suited to other applications. However, the logistics of sorting cheese for different uses based on location within a 291-kg block is probably not practical in a high-speed industrial cheese-cutting operation; therefore, it would be better to develop a technological solution that could achieve more uniform composition and flavor within the 291-kg blocks.

Back to Article Outline

Acknowledgments 

Funding provided in part by the California Dairy Research Foundation (Davis, CA) and Dairy Management, Inc. (Rosemont, IL). Our sincere thanks to the companies that participated in this study. Manuscript FSR0708-2007 of the Department of Food Science, North Carolina State University. The use of trade names in the publication does not imply endorsement by these organizations or criticisms of ones not mentioned.

Back to Article Outline

Supplementary data 

Interpretive summary.

Back to Article Outline

References 

  1. AOAC (Association of Official Analytical Chemists). Official Methods of Analysis. 17th ed.. Gaithersburg, MD: AOAC Int.; 2000;
  2. Avsar Y, Karagul-Yüceer Y, Drake MA, Singh T, Yoon Y, Cadwallader KR. Characterization of nutty flavor in Cheddar cheese. J. Dairy Sci. 2004;87:1999–2010
  3. Banks JM, Brechany EY, Christie WW. Cheddar cheese flavor and chemical indices: Changes during maturation. In:  Cheese E,  Malin L,  Tunkick MH editor. Chemistry of Structure-Function Relationships. New York, NY: Plenum Press; 1995;p. 99–120
  4. Barbano DM. Moisture nonuniformity and sampling errors in large Cheddar cheese blocks. J. AOAC Int. 2001;84:613–619
  5. Brown JA, Foegeding EA, Daubert CR, Drake MA, Gumpertz M. Relationships among rheological and sensorial properties of young cheeses. J. Dairy Sci. 2003;86:3054–3067
  6. Bynum DG, Barbano DM. Whole milk reverse osmosis retentates for Cheddar cheese manufacture: Chemical changes during aging. J. Dairy Sci. 1985;68:1–10
  7. Carunchia Whetstine ME, Cadwallader KR, Drake MA. Characterization of aroma compounds responsible for rosy/floral flavor in Cheddar cheese. J. Agric. Food Chem. 2005;53:3126–3132
  8. Chen M, Irudayaraj J, McMahon DJ. Examination of full fat and reduced fat Cheddar cheese during ripening by Fourier transform infrared spectroscopy. J. Dairy Sci. 1998;81:2791–2797
  9. Christensen K, Reineccius G. Aroma extract dilution analysis of aged Cheddar cheese. J. Food Sci. 1995;60:218–220
  10. Church FC, Swaisgood HC, Porter DH, Catignani GL. Spectrophotometric assay using o-phthaldialdehyde for determination of proteolysis in milk and isolated milk proteins. J. Dairy Sci. 1983;66:1219–1227
  11. Diehl KC, Hamann DD, Whitfield JK. Structural failure in selected raw fruits and vegetables. J. Texture Stud. 1979;10:371–400
  12. Drake MA, Keziah MD, Gerard PD, Delahunty CM, Sheehan C, Turnbull RP, et al. Comparison of cross-cultural differences between lexicons for descriptive analysis of Cheddar cheese flavor in Ireland, New Zealand, and the United States. Int. Dairy J. 2005;15:473–483
  13. Drake MA, McIngvale S, Gerard P, Cadwallader KR, Civille GV. Development of a descriptive language for Cheddar cheese. J. Food Sci. 2001;66:1422–1427
  14. Glantz SA, Slinker BK. Multicolinearity and what to do about it. Primer of Applied Regression and Analysis of Variance. 2nd ed.. New York, NY: McGraw-Hill; 2001;Pages 185–187
  15. Grappin R, Rank TC, Olson NF. Primary proteolysis of cheese during ripening: A review. J. Dairy Sci. 1985;68:531–540
  16. Hort J, Le Grys G. Developments in the textural and rheological properties of UK Cheddar cheese during ripening. Int. Dairy J. 2001;11:475–481
  17. Lawrence RC, Gilles J, Creamer LK, Crow VL, Heap HA, Honore CG, et al. Cheddar cheese and related dry-salted cheese varieties. In:  Fox PFP,  McSweeney LH,  Cogan TM,  Guinee TP editor. Cheese Chemistry, Physics, and Microbiology. 3rd ed.. London, UK: Elsevier; 2004;p. 71–102
  18. Lehninger AL. Biochemistry: The Molecular Basis of Cell Structure and Function. New York, NY: Worth Publishers; 1970;
  19. Lucey JA, Johnson ME, Horne DS. Invited review: Perspectives on the basis of the rheology and texture properties of cheese. J. Dairy Sci. 2003;86:2725–2743
  20. Lynch JM, Barbano DM, Fleming JR. Determination of the total nitrogen content of hard, semihard and processed cheese by the Kjeldahl method: Collaborative study. J. AOAC. 2002;85:445–455
  21. In:  Marshall RT editors. Standard Methods for the Examination of Dairy Products. 16th ed.. Washington, DC: Am. Publ. Health Assoc.; 1992;
  22. Meilgaard MC, Civille GV, Carr BT. Selection and training of panel members. Sensory Evaluation Techniques. 3rd ed.. Boca Raton, FL: CRC Press; 1999;Pages 174–176
  23. Mesa-Dishington JK, Barbano DM, Aplin RD. Cheddar cheese manufacturing costs: Economies of size and effects of different technologies. In: Agric. Res. Bull.. Cornell University, Ithaca, NY. 1987;p. 1–8587-3
  24. Milo C, Reineccius G. Identification and quantification of potent odorants in regular-fat and low-fat mild Cheddar cheese. J. Agric. Food Chem. 1997;45:3590–3594
  25. Olabi A, Barbano DM. Temperature-induced moisture migration in reduced-fat Cheddar cheese. J. Dairy Sci. 2002;85:2114–2121
  26. Pastorino AJ, Hansen CL, McMahon DJ. Effect of pH on the chemical composition and structure-function relationships of Cheddar cheese. J. Dairy Sci. 2003;86:2751–2760
  27. Rank TC, Grappin R, Olson NF. Secondary proteolysis of cheese during ripening: A review. J. Dairy Sci. 1985;68:801–805
  28. Rehman SU, Farkye N, Drake MA. Reduced-fat Cheddar cheese from a mixture of cream and liquid milk protein concentrate. Int. J. Dairy Technol. 2003;56:94–98
  29. Reinbold RS, Ernstrom CA. Effect of nonuniform cooling on moisture, salt, and pH distribution in 290-kg blocks of stirred curd Cheddar cheese. J. Dairy Sci. 1988;71:1499–1506
  30. Reinbold RS, Ernstrom CA, Hansen CL. Temperature, pH, and moisture profiles during cooling of 290-kg stirred curd Cheddar cheese blocks. J. Dairy Sci. 1992;75:2071–2082
  31. Singh T, Drake MA, Cadwallader KR. Flavor of Cheddar cheese: A chemical and sensory perspective. Compr. Rev. Food Sci. Food Saf. 2003;2:139–162
  32. Steffe JF. Rheological Methods in Food Process Engineering. 2nd ed.. East Lansing, MI: Freeman Press; 1996;
  33. Swaisgood HE. Chemistry of the caseins. In:  Fox PF editors. Advanced Dairy Chemistry: Proteins. Vol. 1:Essex, UK: Elsevier Applied Science; 1992;p. 63–110
  34. Troung VD, Daubert CR. Textural characterization of cheeses using vane rheometry and torsion analysis. J. Food Sci. 2001;66:716–721
  35. Van den Dool H, Kratz P. A generalization of the retention index system including linear programmed gas liquid partition chromatography. J. Chromatogr. 1963;11:463–471
  36. Yvon M, Rijnen L. Cheese flavour formation by amino acid catabolism. Int. Dairy J. 2001;11:185–201
  37. Zehentbauer G, Reineccius G. Determination of key aroma components of Cheddar cheese using dynamic headspace dilution assay. Flavour Fragr. J. 2002;17:300–305

PII: S0022-0302(07)71757-5

doi:10.3168/jds.2006-755

Journal of Dairy Science
Volume 90, Issue 7 , Pages 3091-3109, July 2007