Advertisement
Research| Volume 106, ISSUE 1, P151-167, January 2023

Download started.

Ok

The role of packaging on the flavor of fluid milk

Open AccessPublished:November 07, 2022DOI:https://doi.org/10.3168/jds.2022-22060

      ABSTRACT

      Few studies have addressed the effects of package material in the absence of light on contributions to fluid milk flavor. The objective of this study was to compare the sensory and chemical properties of fluid milk packaged in paperboard cartons, low-density polyethylene, high-density polyethylene (HDPE), polyethylene terephthalate (PET), linear low-density polyethylene (LLDPE), and glass. Pasteurized (high temperature short time, 77°C for 25 s) skim and whole milk were filled (280 mL ± 10 mL) into paperboard cartons, low-density polyethylene, HDPE, PET, LLDPE, and glass (control). Milks were stored at 4°C in the dark and sampled at d 0, 5, 10, and 15. Descriptive analysis was applied to document sensory profiles at each time point, and volatile compounds were extracted and identified by solid-phase microextraction with gas chromatography mass spectrometry and gas chromatography-olfactometry. Tetrad tests with consumers were conducted at d 10. Both skim and whole milks packaged in cartons had noticeable paperboard flavor by d 5 and higher levels of hexanal than skim and whole milks in other package types at d 5. Skim milks packaged in paperboard cartons and LLDPE had distinct refrigerator/stale flavor compared with milks in the other package types, concurrent with increased levels of refrigerator/package-related compounds including styrene, acetophenone and 2-ethyl-1-hexanol. Milks packaged in glass, PET and HDPE were not distinguished by consumers at d 10 post-processing. Package type influences fluid milk flavor, and these effects are greater in skim milk compared with whole milk. Paperboard cartons do not preserve milk freshness, as well as PET, HDPE, or glass, due to flavor migration and scalping. Glass remains an ideal barrier to preserve fluid milk flavor, but in the absence of light, HDPE and PET provide additional benefits while also maintaining fluid milk flavor.

      Key words

      INTRODUCTION

      The main concerns for the dairy industry for pasteurized fluid milk are preserving the quality and safety of the final product (
      • Dainelli D.
      • Gontard N.
      • Spyropoulos D.
      • Zondervan-van den Beuken E.
      • Tobback P.
      Active and intelligent food packaging: legal aspects and safety concerns.
      ), while also preserving the characteristic pleasant flavor of milk (
      • Alvarez V.B.
      Fluid milk and cream products, Chapter 5.
      ). The characteristic flavor active compounds of pasteurized fluid milk are created during thermal processing (
      • Shimamura T.
      • Uke H.
      Maillard reaction in milk—Effect of heat treatment. Chapter 5.
      ). However, the flavor of milk has changed over time due, in part, to changes in packaging used to transport milk to the consumer (
      • Azzara C.D.
      • Campbell L.B.
      Off-flavors of dairy products, Chapter 13.
      ;
      • Dainelli D.
      • Gontard N.
      • Spyropoulos D.
      • Zondervan-van den Beuken E.
      • Tobback P.
      Active and intelligent food packaging: legal aspects and safety concerns.
      ). Light oxidation from transparent containers, as well as scalping and migration from different packaging types, affect the flavor of milk and have led to consumer dissatisfaction (
      • Johnson D.S.
      • Duncan S.E.
      • Bianchi L.M.
      • Chang H.H.
      • Eigel W.N.
      • O'Keefe S.F.
      Packaging modifications for protecting flavor of extended-shelf-life milk from light.
      ;
      • Potts H.L.
      • Amin K.N.
      • Duncan S.E.
      Retail lighting and packaging influence consumer acceptance of fluid milk.
      ;
      • Schiano A.N.
      • Jo Y.
      • Barbano D.M.
      • Drake M.A.
      Does vitamin fortification affect light oxidation in fluid skim milk?.
      ).
      Pasteurized fluid milk is more susceptible to packaging-related off-flavors than many other beverages because of its relatively bland taste. The dairy industry, as such, has focused on light blocking options in regards to packaging. Light-oxidized flavor is a well-established and well-researched off-flavor in retail fluid milk in transparent containers, and is caused by photo oxidation of riboflavin and other naturally occurring photosensitive compounds in milk (
      • Northrop-Clewes C.A.
      • Thurnham D.I.
      The discovery and characterization of riboflavin.
      ;
      • Walsh A.M.
      • Duncan S.E.
      • Potts H.
      • Gallagher D.L.
      Comparing quality and emotional responses as related to acceptability of light-induced oxidation flavor in milk.
      ;
      • Brothersen C.
      • McMahon D.J.
      • Legako J.
      • Martini S.
      Comparison of milk oxidation by exposure to LED and fluorescent light.
      ). Paperboard cartons are generally considered to have good light barrier properties compared with other fluid milk package materials, such as polyethylene terephthalate (PET), high-density polyethylene (HDPE), or glass (
      • Potts H.
      • Duncan S.
      • Johnson M.L.
      • O'Keefe S.F.
      • Marcy J.E.
      • Mallikarjunan K.
      0712 Oxygen barrier and light interference packaging properties for controlling light-induced oxidation in milk.
      ,
      • Potts H.L.
      • Amin K.N.
      • Duncan S.E.
      Retail lighting and packaging influence consumer acceptance of fluid milk.
      ).
      In addition to the light-oxidized flavor, off odors and flavors in milk can be created or enhanced, or both, by flavor active compounds from packaging itself via migration and scalping (
      • Bassette R.
      • Fung D.Y.C.
      • Mantha V.R.
      • Marth E.H.
      Off-flavors in milk.
      ). Food migration is the transfer of unwanted chemical contaminants such as trace amounts of solvents, byproducts, additives, or plastic monomers, or all these contaminants, into a food from the surrounding environment or packaging (
      • Kim-Kang H.
      Volatiles in packaging materials.
      ). Food scalping is the transfer of desirable volatile food flavors and aromas to the packaging through absorption or release of undesirable flavors and aromas into the food from the surrounding packaging (
      • Roland A.M.
      • Hotchkiss J.H.
      Determination of flavor—Polymer interactions by vacuum-microgravimetric method. Chapter 13.
      ). Migration and scalping of flavor active compounds are heavily affected by the volatility and polarity of both the food matrix and package, lipophilicity of the product, physical structure (solid or liquid), and temperature of the product and storage conditions (
      • Sajilata M.G.
      • Savitha K.
      • Singhal R.S.
      • Kanetkar V.R.
      Scalping of flavors in packaged foods.
      ).
      Paperboard cartons and HDPE are the most widely used package types for fluid milk (
      • Brody A.L.
      Packaging milk and milk products, Chapter 20.
      ). Paperboard cartons are widespread because they are malleable and cheaper to produce than many other options (
      • Brody A.L.
      Packaging milk and milk products, Chapter 20.
      ). Paperboard cartons are made from wood chips that have been chemically treated to a pulped consistency. The pulped wood is then poured into molds and left to cure until set in the desired shape. Commonly, paperboard cartons are lined with millimeter-thin layers of low-density polyethylene (LDPE) as a barrier between the paperboard and the product. Low-density polyethylene acts as a fluid barrier because of its resistance to liquids, pH, and alcohol (
      • Nielsen T.J.
      • Jägerstad I.M.
      • Öste R.E.
      Study of factors affecting the absorption of aroma compounds into low-density polyethylene.
      ;
      • Simon M.
      • Hansen A.P.
      Effect of various dairy packaging materials on the shelf life and flavor of pasteurized milk.
      ). Paperboard cartons have potential disadvantages because paperboard is highly permeable to the surrounding environment and can allow volatile compound migration into the container contents (
      • Ziegleder G.
      Volatile and odorous compounds in unprinted paperboard.
      ). Paperboard may also be a direct source of off-flavors because the formation of wood lignin into paperboard is an oxidative process and these compounds can also migrate into product (
      • Ziegleder G.
      Volatile and odorous compounds in unprinted paperboard.
      ). High-density polyethylene, though costlier than paperboard, is commonly used because of its structural stability, moisture resistance, and barrier to the surrounding environment (
      • Crosby N.T.
      Plastic packaging materials, Chapter 1.
      ;
      • Nielsen T.J.
      • Jägerstad I.M.
      • Öste R.E.
      Study of factors affecting the absorption of aroma compounds into low-density polyethylene.
      ). Polyethylene terephthalate is also an alternative for packaging of fluid milk (
      • Brody A.L.
      Packaging milk and milk products, Chapter 20.
      ). Polyethylene terephthalate, when compared with HDPE, is highly resistant to acids, oils, and alcohols, and is generally a good barrier against moisture (
      • Sajilata M.G.
      • Savitha K.
      • Singhal R.S.
      • Kanetkar V.R.
      Scalping of flavors in packaged foods.
      ). Polyethylene terephthalate is not as quite as rigid, durable, or resistant to heat when compared with HDPE, but is cheaper to produce and more malleable, making it easier to manipulate during production and recycling (
      • Sajilata M.G.
      • Savitha K.
      • Singhal R.S.
      • Kanetkar V.R.
      Scalping of flavors in packaged foods.
      ).
      The role of flavors from packaging (scalping or migration) on the sensory quality of fluid milks have not been fully investigated. Previous research has shown that salted and unsalted butter absorbed flavors from packaging and the surrounding environment at 4°C storage (
      • Lozano P.R.
      • Miracle E.R.
      • Krause A.J.
      • Drake M.A.
      • Cadwallader K.R.
      Effect of cold storage and packaging material on the major aroma components of sweet cream butter.
      ).
      • Leong C.M.O.
      • Harte B.R.
      • Partridge J.A.
      • Ott D.B.
      • Downes T.W.
      Off-flavor development in milk packaged in polyethylene-coated paperboard cartons.
      demonstrated that whole, low-fat, and skim milk packaged in polyethylene-coated half-pint paperboard cartons stored at 4°C had package off-flavor by d 1 (
      • Leong C.M.O.
      • Harte B.R.
      • Partridge J.A.
      • Ott D.B.
      • Downes T.W.
      Off-flavor development in milk packaged in polyethylene-coated paperboard cartons.
      ).
      • Leong C.M.O.
      • Harte B.R.
      • Partridge J.A.
      • Ott D.B.
      • Downes T.W.
      Off-flavor development in milk packaged in polyethylene-coated paperboard cartons.
      focused on off-flavors contributed from polyethylene-coated paperboard cartons during fluid milk shelf life and did not conduct volatile compound analysis or compare with other packaging materials.
      • Sipple L.R.
      • Barbano D.M.
      • Drake M.A.
      Maintaining and growing fluid milk consumption by children in school lunch programs.
      evaluated child consumer perception of extrinsic and intrinsic parameters of fluid milk in single-serve half pints of different package materials. They demonstrated that school children preferred skim or 1% unflavored milk packaged in HDPE or PET containers over traditional coated paperboard after approximately 10 d of cold storage. Volatile compound differences, multiple time points, and whole milk were not evaluated. Different package materials for fluid milk have not been directly compared in the absence of light, nor have volatile compounds been determined. The objective of this study was to determine the role of package on the chemical and sensory properties of milk in the absence of light exposure. The absence of light exposure allows control of the study parameters such that the actual role of package material on migration or scalping, or both, can be clearly identified and sourced.

      MATERIALS AND METHODS

      Experimental Overview

      Whole and skim milk were high temperature short time (HTST) processed and filled into the following 6 different packages: amber glass (TraceClean Wide Mouth Packers VWR), LDPE (Wide Mouth, VWR), HDPE (Wide Mouth, VWR), PET (Storage Bottles, Corning), LLDPE (Gusset Bag, Elkay Plastics Co.), or paperboard carton (Evergreen Barrier Paperboard, Evergreen Packaging), and stored at 4°C. Sensory properties (descriptive sensory analysis) and volatile compound analysis [headspace solid-phase microextraction (HS-SPME) GC-MS] were conducted at time 0 and d 5, 10, and 15. Consumer difference testing (tetrad tests) were conducted for whole and skim fluid milk packaged in HDPE versus glass, PET versus glass, and paperboard carton versus glass at d 10. The entire experiment was repeated 3 times.

      Milk Processing

      Raw milk (<10,000 cfu/mL aerobic plate count, <100,000 SCC) was standardized (whole milk had 3.25 ± 0.2% fat and skim milk had <0.10 ± 0.05% fat) and processed at the North Carolina State University Dairy Enterprise System (Raleigh, NC). Skim and whole milk were HTST processed (77°C for 25 s) and filled into the following 6 different 280-mL (± 10 mL) packaging types: amber glass, LDPE, HDPE, PET, LLDPE, or polyethylene-lined paperboard carton. All packaging types were sanitized before filling in a 100-mg/kg (NaOCl) sanitizer/water solution (XY-12, Ecolab). Processing and filling were done with minimal light to reduce light exposure. Packages were filled to the same volume, as follows: 230 mL (± 10 mL) with 50 (± 10 mL) of headspace. Paperboard cartons were filled using an Evergreen EQ-2 gable top filler (Evergreen Packaging) system. All other packaging types were hand filled using a Micro Thermics Ultra-Clean Fill Hood (MicroThermics). Milks were stored at 4°C and covered with black laser cloth to shield from light. Milks were sampled for sensory and instrumental analysis on d 0, 5, 10, 15 post-processing. Total aerobic plate and coliform counts (coliform and aerobic plate count Petrifilms; 3M) were conducted on milks at d 0, 10, 15 to monitor microbial quality.

      Dissolved Oxygen

      The dissolved oxygen of skim and whole milk of each package type was measured using a benchtop dissolved oxygen analyzer (YSI Model 5100) at each time point (0, 5, 10, 15 d). The oxygen probe was placed in the sample, with both ends of the probe submerged, and stirred using the self-stirring probe BOD probe attachment. The reading was then allowed to stabilize for dissolved oxygen and a measurement was taken. All measurements were done at 21°C in triplicate for skim and whole milk on each package type and each experimental replication. Measurements were given in milligrams per liter of oxygen in sample. Calibration for dissolved oxygen was done initially before each time point by placing the probe in a known oxygen environment (250 mL of water in a bottle aerated for 15 min with pure 02, the probe was stabilized for 15 min in the aerated water). The probe was adjusted if readings of known oxygen environment were out of specifications. Barometric pressure was adjusted to current elevation and temperature in the setting before each calibration

      Headspace Oxygen

      The headspace oxygen of skim and whole milk of each package type was measured using a benchtop oxygen analyzer (Quantek Model 905) at each time point (0, 5, 10, 15 d). To prevent oxygen leaking in or out of the packaging, the package was first fitted with a foam rubber septum (Foam Septum). The package was then punctured using a 20/2”/5 leur needle (N720 Hamilton Needle) with a side-slotted opening. The oxygen analyzer was then turned on until a stable oxygen reading was established. Readings were conducted at 21°C. All measurements were done in triplicate for each package type and each experimental replication for skim and whole milk. A calibration was done initially before each time point by filling a heat-sealed plastic bag (size 2” × 2”, thickness 0.0009”, Gusset Bag, Elkay Plastics Co.) with pure N2 through a foam rubber septum, and then measuring the atmosphere in the bag. The reading was adjusted to 0% oxygen (± 0.1%) with the ZERO potentiometer on the back of the analyzer as needed. The atmosphere of the room was then used to check calibration between each sample, and was re-adjusted or re-calibrated if necessary.

      Descriptive Sensory Analysis

      Samples were prepared with the overhead light off to minimize potential lights effects on milk. Sensory profiles of milks from each package type at each time point were evaluated by 7 experienced and trained panelists, using an established sensory language for fluid milk (
      • Lozano P.R.
      • Drake M.A.
      • Benitez D.
      • Cadwallader K.R.
      Instrumental and sensory characterization of heat-induced odorants in aseptically packaged soy milk.
      ;
      • Lee A.P.
      • Barbano D.M.
      • Drake M.A.
      The influence of ultra-pasteurization by indirect heating versus direct steam injection on skim and 2% fat milks.
      ;
      • McCarthy K.S.
      • Lopetcharat K.
      • Drake M.A.
      Milk fat threshold determination and the effect of milk fat content on consumer preference for fluid milk.
      ). Each panelist had a minimum of 90 h of previous experience evaluating the sensory profile of different fat levels of HTST fluid milk. Attributes were scored using a 0- to 15-point scaling consistent with the Spectrum Method (
      • Meilgaard M.
      • Civille G.
      • Carr B.
      Descriptive analysis techniques. Chapter 10.
      ). Each milk was evaluated at each time point (0, 5, 10, 15 d) in duplicate by each panelist. Milks were poured (40 mL) into soufflé cups (65 mL) numbered with randomized 3-digit codes. Panelists evaluated milks from each experimental replication of each fat type (6 milks) in a randomized order in 1 session. A 3-min rest between samples was enforced. Paper ballots were used for data collection.

      Volatile Compound Analysis of Fluid Milk

      HS-SPME-GCMS

      Volatile compounds were extracted from milks at each time point (0, 5, 10, 15 d) by HS-SPME, followed by GC-MS with a ZB-5ms (30 m × 0.25 mm i.d. × 0.25 µm; Phenomenex) and a ZB-Wax (30 m × 0.25 mm i.d. × 0.25 µm) column (Phenomenex) with MS (7820 GC 5975 MSD, Agilent Technologies Inc.). Samples were injected via a CTC Analytics CombiPal Autosampler (CTC Analytics). Extraction methods were modified from
      • Jo Y.
      • Benoist D.M.
      • Barbano D.M.
      • Drake M.A.
      Flavor and flavor chemistry differences among milks processed by high-temperature, short-time pasteurization or ultra-pasteurization.
      . Five grams of milk was measured in triplicate in SPME vials (Microliter Analytical), with 20 uL of internal standard added 81 mg/kg 2-methyl-3-heptanone (Sigma-Aldrich) in methanol and 167 mg/kg 2-methylpentanoic (Sigma-Aldrich) in methanol. Vials were equilibrated at 40°C for 20 min with 3 s of 350 rpm agitation. A single DVB/Carboxen/PDMS 1.0-cm fiber (Supelco) was used for all analysis. The SPME fiber was exposed to the samples for 25 min at depth 3.1 cm. The fiber was retracted and injected at 5.0 cm in the GC inlet for 7.5 min. The GC oven was held at 40°C for 3 min with a gradient of 10°C/min to 90°C with no hold, and the gradient was increased at a rate of 5°C/min to 200°C and held for 10 min, and then increased to 20°C/min 250°C and held for 5 min at completion. Helium was pressurized at 7 psi in the column and maintained at a constant flow rate of 1 mL/min. A standard scan mode was used (30 to 350 m/z) to identify compounds of interest for initial evaluation. Selective ion monitoring mode was subsequently used for compound identification using authentic standards.

      HS-SPME GC-Olfactometry

      The HS-SPME GC-olfactometry was used to identify aroma active compounds. Gas chromatography-olfactometry was conducted on an Agilent 6850 GC (Agilent Technologies Inc.) with a ZB-5ms (30 m × 0.25 mm i.d. × 0.25 µm) column (Phenomenex) or a ZB-Wax plus (30 m × 0.25 mm i.d. × 0.25 µm) column (Phenomenex). Ten milliliters of milk was added to 40-mL vials in duplicate. Vials were equilibrated for 30 min at 40°C and SPME fibers were then exposed to the headspace of the vials for 30 min at 2 cm. Fibers were manually injected at 3-cm depth. The fiber was initially exposed to the inlet for 10 min at 250°C. The GC program was held at 40°C for 3 min and increased gradually to 150°C over 11 min (10°C/min). The temperature was then increased to 250°C over 3.5 min (30°C/min) and then held for 10 min. Carrier gas used was helium at a flow rate of 1.5 mL/min. The splitter was used at a 1:1 ratio between the flame ionization detector and sniffer port. Two experienced sniffers recorded aroma events, retention times, and perceived aroma intensities using a 0- to 5-point post-peak intensity scale.

      Compound Identification and Quantitation

      Compound identification was achieved by MS (
      • NIST (National Institute of Standards and Technology)
      NIST Wiley Registry: NIST Mass Spectral Library.
      ) and odor properties of authentic standards. Retention indices of volatile compounds were calculated (
      • van Den Dool H.
      • Kratz P.D.
      A generalization of the retention index system including linear temperature programmed gas-liquid partition chromatography.
      ) using a C6–C22 alkane series (Sigma-Aldrich) under identical GC olfactometry/MS conditions. External standard curves for 9 selected compounds (styrene, xylene, 2-ethyl-1-hexanol, hexanal, nonanal, p-cymene, acetophenone, toluene, octanal) were generated using skim or whole (HTST) milk packaged in glass as the matrix. Each stock solution was prepared using authentic standards (Sigma-Aldrich) and diluted (methanol, ethanol, or water based on solubility) at separate points for a 5-point standard curve. Each matrix was also spiked with 20 uL of an internal standard (81 mg/kg 2-methyl-3-heptanone in methanol), with the ratio of response calculated along with the external standard in the matrix.

      Direct Solvent Extraction Solvent-Assisted Flavor Extraction

      Direct solvent extraction with solvent-assisted flavor extraction (DSE-SAFE) was chosen to evaluate skim milk stored at 4°C for 15 d from paperboard cartons and glass (control) because of its more sensitive recovery of semi-volatile compounds and higher molecular weight compounds (
      • Engel W.
      • Bahr W.
      • Schieberle P.
      Solvent assisted flavour evaporation—A new and versatile technique for the careful and direct isolation of aroma compounds from complex food matrices.
      ;
      • Lozano P.R.
      • Drake M.A.
      • Benitez D.
      • Cadwallader K.R.
      Instrumental and sensory characterization of heat-induced odorants in aseptically packaged soy milk.
      ;
      • Xie Z.
      • Liu Q.
      • Liang Z.
      • Zhao M.
      • Yu X.
      • Yang D.
      • Xu X.
      The GC/MS analysis of volatile components extracted by different methods from Exocarpium Citri Grandis.
      ). Direct solvent extraction with solvent-assisted flavor extraction was used for the identification and determination of specific compounds related to paperboard cartons that were identified at the limits of quantification by SPME-GCMS.
      Fluid milk was extracted using a modified method of direct solvent extraction, as described by
      • Engel W.
      • Bahr W.
      • Schieberle P.
      Solvent assisted flavour evaporation—A new and versatile technique for the careful and direct isolation of aroma compounds from complex food matrices.
      . Three hundred milliliters of fluid skim milk was separated into 3 equal parts into 3 250-mL Teflon bottles (Thermo Fisher Scientific). One hundred mL of ethyl ether (Sigma-Aldrich) and 40 uL of each internal standard was added to each bottle (2-methy-3-heptanone, 81 mg/kg, neutral/basic fraction, 2-methylpentanoic acid 167 mg/kg acid-fraction; Sigma-Aldrich), and 30 g of NaOH (Sigma-Aldrich) was added to break the emulsion. Each bottle was sealed and then shaken vigorously for 2 min, untightened to let out excess pressure, and then placed for 30 min on a platform mixer (Rotomix type 50800 Thermolyne) at the highest setting. Bottles were then centrifuged at 2,500 × g at 21°C (Thermo Lynx 6000, Thermo Fisher Scientific) for 30 min. The top aqueous layer (ethyl ether) was then removed and dried using anhydrous sodium sulfate (Na2SO4; Sigma-Aldrich) to remove excess water. All samples were concentrated to 100 mL under nitrogen gas.
      Concentrated samples were then distilled using a SAFE apparatus, assembled as previously described (
      • Engel W.
      • Bahr W.
      • Schieberle P.
      Solvent assisted flavour evaporation—A new and versatile technique for the careful and direct isolation of aroma compounds from complex food matrices.
      ;
      • Lozano P.R.
      • Drake M.A.
      • Benitez D.
      • Cadwallader K.R.
      Instrumental and sensory characterization of heat-induced odorants in aseptically packaged soy milk.
      ). A round 1-L flask was attached to the main SAFE head apparatus and submerged in a 40°C water bath to maintain temperature. The SAFE head apparatus was then attached to 2 glass traps, in series, for collection of the sample. Both traps were submerged in liquid nitrogen-filled Nalgene buckets that were constantly filled to maintain temperature. The SAFE Apparatus (Ace Glassware) was then brought to 10−5 Torr vacuum, and attached to the end of the second trap in series. Sample was introduced slowly (over 15 min) via a stopcock into the main SAFE apparatus head into the 1-L flask, volatilizing the sample, and pulling it through the SAFE apparatus and condensing it in the first trap. After all sample was fully introduced, the SAFE apparatus was maintained at constant vacuum (10−5 Torr) and liquid nitrogen was maintained in Nalgene buckets for 2 h for maximum retention. The primary trap with condensed sample was collected and allowed to thaw at room temperature. Sample was then dried with anhydrous sodium sulfate (Na2SO4) in an amber glass jar (250 mL) and then concentrated to 20 mL under a stream of nitrogen gas. The concentrated sample was then transferred to a screw-top glass tube for phase separation.
      To aid in compound detection, all samples were separated into 2 separate fractions, neutral/basic and acidic. The concentrated sample was first washed with 3 mL of 0.5 M sodium bicarbonate (Thermo Fisher Scientific) and shaken for 10 min. The aqueous layer was removed (bottom layer) after each washing and collected. The sodium bicarbonate wash was repeated twice. The sample was then washed with 2 mL of saturated NaCl solution (23% wt/wt; VWR) and shaken for 10 min with the aqueous layer removed (bottom layer) and collected into the same glass tube. The ether phase remaining (upper layer) was the resultant neutral/basic fraction and was then concentrated to 2 mL under a stream of nitrogen gas, dried with anhydrous sodium sulfate (Na2SO4) into a new 2 mL glass vial, and concentrated to 0.5 mL under nitrogen gas.
      The transferred aqueous layer (bottom layer) in the glass tube was acidified through the addition of 6 M hydrochloric acid (Sigma-Aldrich) to a pH of approximately 2.0. Five milliliters of diethyl ether was then added to the tube and shaken for 10 min and the ether phase removed. This process was repeated 3 times. The collected ether phase was the acidic fraction. The acidic fraction was then concentrated to 2 mL under a stream of nitrogen gas, dried with anhydrous sodium sulfate (Na2SO4), and concentrated to 0.5 mL under nitrogen gas.

      Compound Identification for DSE-SAFE

      Three microliters of the acidic and neutral/basic fraction were injected in split-less mode in duplicate on a ZB-5ms (30 m × 0.25 mm i.d. × 0.25 µm) column and a ZB-Wax plus (30 m × 0.25 mm i.d. × 0.25 µm) column installed in an Agilent 7820-GC paired with a 5975-MSD. All injections had a 2.5-min solvent delay to protect the instrument. The GC oven was at held at 40°C for 3 min with a gradient of 10°C/min to 90°C with no hold, the gradient was increased to a rate of 5°C/min to 200°C and held for 10 min, and the gradient was then increased to 20°C/min 250°C and held for 5 min at completion. Selective ion mode and scan mode were used in tandem to identify and quantify compounds of interest. Compound identities were confirmed with authentic standards. Internal standards were used to calculate relative abundance for compounds of interest.

      Consumer Difference Tests (Tetrads)

      Sensory testing was conducted in accordance with the North Carolina State University Institutional Review Board for human subject regulations. Tetrad difference testing (
      • ASTM International
      E3009–15e1 Standard Test Method for Sensory Analysis—Tetrad Test.
      ; E3009) was conducted for both skim and whole milks with 4 packaging types, as follows: paperboard cartons, HDPE, PET, and glass (control) at d 10. Consumers were presented the following 3 separate tetrad tests for whole and skim milk (6 tetrads): HDPE versus glass, PET versus glass, paperboard carton versus glass at d 10 of shelf life. Package types for tetrads were selected based on descriptive analysis and realistic practical applications. High-density polyethylene and PET packages were selected because they make up more than 85% of fluid milk packaging sales (
      • USDA
      Packaged Fluid Milk Sales in Federal Milk Order Markets: By Size and Type of Container and Distribution Method.
      ). Paperboard cartons contribute close to 14% of fluid packaging sales with the majority of that being used for the school lunch programs (
      • USDA
      Packaged Fluid Milk Sales in Federal Milk Order Markets: By Size and Type of Container and Distribution Method.
      ). Consumers (n = 50) were recruited via emails from the North Carolina State University Sensory Service Center database with the qualifier being that they were self-reported regular consumers of fluid milk. During testing, each consumer was instructed both verbally and provided with an instruction guide describing how a tetrad test is performed. They were asked to pair the 4 coded samples based on which ones they believed to be the most similar to each other. Testing was replicated with 2 of the experimental replications for skim and whole milk for all 3 tetrad pairs (n = 100 consumers through each tetrad).
      Milks (60 mL) were evaluated by consumers in 128-mL Styrofoam cups coded with 3-digit codes. Milks were served at 4°C. Each panelist was given 4 coded samples at once and instructed to sample them in the order presented left to right. A 3-min enforced rest occurred between each tetrad during which consumers were requested to rinse their mouths with water and to take a bite of unsalted cracker. Consumers were asked to complete 3 sets of tetrads (1 fat content level) in 1 seating, with the order of presentation of each tetrad test randomized between each consumer. Consumers were compensated at the completion of the test with a $5 gift card. Compusense Cloud (Compusense) was used for data collection.

      Statistical Analysis

      Data analysis was conducted at a 95% confidence (P < 0.05). Analysis software used included XLSTAT version 2019.1.3 (Addinsoft) and SAS version 9.4 (SAS Institute). Two separate analyses were conducted, 1 with whole milk and 1 with skim. Each experiment was designed as a randomized complete block design with a time point by package type factorial arrangement of treatments and evaluated as a mixed model ANOVA (SAS). Principal component analysis was used to visualize differences among volatile compounds among package types within each fluid milkfat content (XLSTAT). Consumer difference data were evaluated using the minimum correct judgement for significance (P < 0.05;
      • ASTM International
      E3009–15e1 Standard Test Method for Sensory Analysis—Tetrad Test.
      ; E3009–15).

      RESULTS AND DISCUSSION

      Milk Microbial Quality

      No coliforms were detected in pasteurized whole or skim milks. All milks had an aerobic plate count of less than 102 cfu/mL initially, and less than 104 cfu/mL at d 15.

      Dissolved and Headspace Oxygen Analysis

      Dissolved oxygen concentration decreased over storage time for all packaging types, regardless of milkfat content (P < 0.05; Tables 1, 2). Dissolved oxygen decreases were consistent with previous studies found for fluid milk (
      • Zygoura P.
      • Moyssiadi T.
      • Badeka A.
      • Kondyli E.
      • Savvaidis I.
      • Kontominas M.G.
      Shelf life of whole pasteurized milk in Greece: Effect of packaging material.
      ;
      • Karatapanis A.E.
      • Badeka A.V.
      • Riganakos K.A.
      • Savvaidis I.N.
      • Kontominas M.G.
      Changes in flavour volatiles of whole pasteurized milk as affected by packaging material and storage time.
      ). Dissolved oxygen for skim and whole milks packaged in LLDPE decreased faster compared with milks in other packages, presumably due to the greater permeability of the LLDPE (P < 0.05; Tables 1, 2, 3). Headspace oxygen (% O2) concentration decreased over storage time for glass, PET, HDPE, and LDPE, regardless of milkfat content (P < 0.05; Tables 4, 5). Headspace oxygen in milks packaged in cartons and LLDPE did not decrease over storage time due to the high permeability of both packaging types (P > 0.05; Tables 3, 4, 5).
      Table 1Concentration of dissolved oxygen (mg/L) in skim milk in different packages across 4 time points (d 0, 5, 10, 15)
      Package type
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Time point
      d 0d 5d 10d 15
      Glass10.31 ± 0.20
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      8.68 ± 0.50
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      5.71 ± 0.66
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      3.90 ± 1.20
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      PET10.59 ± 0.08
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      9.26 ± 0.40
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.77 ± 0.43
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      5.12 ± 1.35
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      HDPE10.66 ± 0.10
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      9.37 ± 0.16
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.70 ± 0.35
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      4.36 ± 0.50
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LDPE10.61 ± 0.13
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      9.00 ± 0.66
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.04 ± 0.67
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      3.90 ± 0.32
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LLDPE10.46 ± 0.35
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.08 ± 0.61
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      4.73 ± 0.96
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      0.96 ± 0.25
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      Carton11.15 ± 0.13
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      10.06 ± 0.26
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      8.31 ± 0.59
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.02 ± 0.99
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      a–d Means in the same column followed by a different superscript are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Table 2Concentration of dissolved oxygen (mg/L) in whole milk in different packages across 4 time points (d 0, 5, 10, 15)
      Package type
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Time point
      d 0d 5d 10d 15
      Glass11.94 ± 0.08
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      10.63 ± 0.37
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.89 ± 0.48
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      5.78 ± 0.36
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      PET12.19 ± 0.11
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      10.09 ± 1.04
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      8.80 ± 1.01
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      4.72 ± 0.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      HDPE12.18 ± 0.16
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      11.04 ± 0.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      8.94 ± 0.56
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      5.37 ± 0.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LDPE12.16 ± 0.14
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      10.57 ± 0.36
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.75 ± 0.65
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      5.09 ± 0.59
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LLDPE11.93 ± 0.07
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      7.57 ± 0.34
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      3.80 ± 0.43
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      1.16 ± 0.18
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      Carton12.10 ± 0.06
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      11.47 ± 0.24
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      10.51 ± 0.26
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      8.38 ± 0.37
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      a–e Means in the same column followed by a different superscript are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Table 3Oxygen transmission rates of packaging types over 24 h
      Oxygen transmission rates were reported from package type manufacturer. Values shown ± SD.
      ItemPackage type
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene.
      GlassPETLDPEHDPELLDPECarton
      Oxygen transmission rate
      Oxygen transmission rate (cc × mm/m2 × 24 h × atm).
      <0.55195721,009801
      Film/packaging wall thickness (mm)1.1 ± 0.11.1 ± 0.11.3 ± 0.21.3 ± 0.20.019 ± 0.010.9 ± 0.2
      1 Oxygen transmission rates were reported from package type manufacturer. Values shown ± SD.
      2 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene.
      3 Oxygen transmission rate (cc × mm/m2 × 24 h × atm).
      Table 4Concentration of headspace oxygen (% O2) in skim milk in different packages across 4 time points (d 0, 5, 10, 15)
      Package type
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Time point
      d 0d 5d 10d 15
      Glass20.2 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.6 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.5 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.4 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      PET20.2 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.4 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      17.7 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      17.5 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      HDPE19.7 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      19.0 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      17.8 ± 0.6
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      17.8 ± 0.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LDPE19.7 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.6 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      17.7 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      17.9 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LLDPE20.8 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.9 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.9 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.9 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      Carton20.6 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.5 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.6 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.7 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      a–d Means in the same column followed by a different superscript are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Table 5Concentration of headspace oxygen (% O2) in whole milk in different packages across 4 time points (d 0, 5, 10, 15)
      Package type
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Time point
      d 0d 5d 10d 15
      Glass20.2 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      19.3 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      19.3 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.8 ± 0.5
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      PET20.1 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.7 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.6 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.3 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      HDPE19.7 ± 0.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      19.0 ± 0.8
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.5 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.5 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LDPE19.7 ± 0.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.8 ± 0.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      19.0 ± 0.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      18.5 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      LLDPE20.9 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.9 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      21.0 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      21.1 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      Carton20.2 ± 0.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.5 ± 0.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.1 ± 0.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      20.3 ± 0.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      a–d Means in the same column followed by a different superscript are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.

      Descriptive Analysis

      Milks (whole and skim) packaged in different packages demonstrated distinct flavor differences initially and with storage time (Tables 6, 7). Skim milks packaged in paperboard cartons and LLDPE had lower cooked and sweet aromatic flavors with storage time and distinct refrigerator/stale flavor when compared with milks in the other packaging types (P < 0.05; Table 6). This flavor difference may be due to the high permeability of paperboard cartons and LLDPE (Table 3). High permeability of packaging has been shown to increase stale flavor and decrease overall flavor acceptability in UHT milks (
      • Wadsworth K.D.
      • Bassette R.
      Effect of oxygen on development of off-flavors in ultrahigh-temperature milk.
      ). Paperboard/cardboard flavor was only detected in milks packaged in cartons, and this flavor was detected within hours of packaging (d 0) and increased with storage time (Table 6). Previous studies have found that paperboard or cardboard off-flavors associated with paperboard cartons are detected by sensory panelists from d 1 to 5 of storage time (
      • Leong C.M.O.
      • Harte B.R.
      • Partridge J.A.
      • Ott D.B.
      • Downes T.W.
      Off-flavor development in milk packaged in polyethylene-coated paperboard cartons.
      ;
      • Karatapanis A.E.
      • Badeka A.V.
      • Riganakos K.A.
      • Savvaidis I.N.
      • Kontominas M.G.
      Changes in flavour volatiles of whole pasteurized milk as affected by packaging material and storage time.
      ). Skim milk packaged in LLDPE had higher refrigerator/stale flavor than skim milk packaged in cartons across storage (P < 0.05; Table 6).
      Table 6Trained panel profiles of skim milks filled into different packages with storage time (d 0, 5, 10, 15)
      TimePackaging type
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene.
      Sensory attribute
      Sensory attribute intensities were scored on a 0- to 15-point universal intensity scale (Meilgaard et al., 2007).
      Overall aromaSweet aromaticCookedPaperboard/ cardboardRefrigerator/ staleSweetSaltyAstringency
      d 0Glass2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
       ND
      ND = not detected.
      ND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND0.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      d 5Glass2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      d 10Glass2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND2.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      d 15Glass2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND0.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      a–d Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene.
      2 Sensory attribute intensities were scored on a 0- to 15-point universal intensity scale (
      • Meilgaard M.
      • Civille G.
      • Carr B.
      Descriptive analysis techniques. Chapter 10.
      ).
      3 ND = not detected.
      Table 7Trained panel profiles of whole milks filled into different packages with storage time (d 0, 5, 10, 15)
      TimePackaging type
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene.
      Sensory attribute
      Sensory attribute intensities were scored on a 0- to 15-point universal intensity scale (Meilgaard et al., 2007).
      Overall aromaSweet aromaticMilkfatCookedPaperboard/ cardboardRefrigerator/ staleSweetSaltyAstringency
      d 0Glass2.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
       ND
      ND = not detected.
      ND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND0.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      d 5Glass2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      d 10Glass2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      0.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      d 15Glass2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      PET2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      HPDE2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.4
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LDPE2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.6
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      NDND2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.7
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.8
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      LLDPE2.3
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND2.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      Carton2.1
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      ND3.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.9
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.0
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1.5
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      2.2
      Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      a–f Means in the same column followed by a different superscript letter are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene.
      2 Sensory attribute intensities were scored on a 0- to 15-point universal intensity scale (
      • Meilgaard M.
      • Civille G.
      • Carr B.
      Descriptive analysis techniques. Chapter 10.
      ).
      3 ND = not detected.
      Whole milks packaged in paperboard carton and LLDPE had lower sweet aromatic flavor across storage time compared with milks filled into glass, PET, HDPE, or LDPE (P < 0.05; Table 7). Similar to skim milk, paperboard/cardboard flavor was only detected in milks packaged in paperboard cartons (Table 7). Unlike skim milk, this flavor was not detected until d 5, but also increased through d 15 (Table 7). This effect has been previously documented in whole milk packaged in paperboard cartons by d 5 (
      • Karatapanis A.E.
      • Badeka A.V.
      • Riganakos K.A.
      • Savvaidis I.N.
      • Kontominas M.G.
      Changes in flavour volatiles of whole pasteurized milk as affected by packaging material and storage time.
      ). Refrigerator/stale flavor was detected in milks packaged in LLDPE at d 0 and in whole milk filled into paperboard carton milk by d 10, and this flavor increased with storage time (Table 7). Though skim and whole milk cannot be directly compared, paperboard/cardboard and refrigerator/stale flavors in milk packaged in paperboard cartons or LLDPE were generally more intense for skim milk at a specific time point compared with whole milk (Table 6, Table 7). This flavor difference between skim and whole milk is consistent with previous studies that off-flavors due to paperboard packaging were more pronounced in skim milk compared with whole milk or low-fat milk (
      • Leong C.M.O.
      • Harte B.R.
      • Partridge J.A.
      • Ott D.B.
      • Downes T.W.
      Off-flavor development in milk packaged in polyethylene-coated paperboard cartons.
      ;
      • Sipple L.R.
      • Barbano D.M.
      • Drake M.A.
      Maintaining and growing fluid milk consumption by children in school lunch programs.
      ). Off-flavor problems in fluid milk due to package-related absorption or migration do increase with a decrease in container size due to the increase in the contact surface area with fluid milk. Comparisons have previously been made among various package sizes such as pint, quart, gallon, and liter (
      • Bradley Jr., R.L.
      Effect of light on alteration of nutritional value and flavor of milk: A review.
      ;
      • Leong C.M.O.
      • Harte B.R.
      • Partridge J.A.
      • Ott D.B.
      • Downes T.W.
      Off-flavor development in milk packaged in polyethylene-coated paperboard cartons.
      ;
      • van Aardt M.
      • Duncan S.E.
      • Marcy J.E.
      • Long T.E.
      • Hackney C.R.
      Effectiveness of poly (ethylene terephthalate) and high-density polyethylene in protection of milk flavor.
      ). Our packaging volume (230 mL) was selected to represent the largest contact area (surface area) or “worst-case” scenario. This volume is also the serving size for school lunch milk in the school lunch program in the United States (8 fl oz, 236 mL) and that milk is also typically served in a paperboard carton (
      • Sipple L.R.
      • Barbano D.M.
      • Drake M.A.
      Maintaining and growing fluid milk consumption by children in school lunch programs.
      ).

      Volatile Compound Analysis

      Skim and whole milks packaged in paperboard cartons or LLDPE had higher concentrations of packaging/refrigerator taint volatile compounds styrene, acetophenone, and 2-ethyl-1-hexanol than milks filled into LDPE, HDPE, PET, or glass (P < 0.05; Tables 8, 9). Styrene concentrations increased with storage across all packaging types except for the glass package. Styrene is a marker for degradation of plastic barriers over storage time (
      • López P.
      • Batlle R.
      • Salafranca J.
      • Nerín C.
      Efficiency of whole and skimmed powdered milk for trapping volatile compounds released from plastic containers in high-temperature applications.
      ), is present in refrigeration environments (
      • Lozano P.R.
      • Miracle E.R.
      • Krause A.J.
      • Drake M.A.
      • Cadwallader K.R.
      Effect of cold storage and packaging material on the major aroma components of sweet cream butter.
      ), and would be expected to migrate in from the environment with permeable containers such as cartons and LLDPE. Ethylbenzene is the main component in the production of styrene, and acetophenone is a byproduct of the oxidation of styrene and ethylbenzene (
      • Azzara C.D.
      • Campbell L.B.
      Off-flavors of dairy products, Chapter 13.
      ;
      • Ziegleder G.
      Volatile and odorous compounds in unprinted paperboard.
      ;
      • Vera P.
      • Canellas E.
      • Nerín C.
      Compounds responsible for off-odors in several samples composed by polypropylene, polyethylene, paper and cardboard used as food packaging materials.
      ). The more permeable barriers of paperboard cartons and LLDPE (Table 3) may be the cause of increased concentrations of styrene and acetophenone, which migrate into the milk from the refrigeration environment. The compound 2-ethyl-1-hexanol does occur naturally in some foods at low concentrations, but, industrially, it is used mainly used as a solvent for household cleaners and as a precursor to plasticizers (
      • Vera P.
      • Canellas E.
      • Nerín C.
      Compounds responsible for off-odors in several samples composed by polypropylene, polyethylene, paper and cardboard used as food packaging materials.
      ). Plasticizers are substances that make materials (in this case plastic) more pliable or flexible (
      • Vera P.
      • Canellas E.
      • Nerín C.
      Compounds responsible for off-odors in several samples composed by polypropylene, polyethylene, paper and cardboard used as food packaging materials.
      ). The compound 2-ethyl-1-hexanol is also known or considered to be an indoor air pollutant coming from the condensation of refrigerators or air-conditioners (
      • Wakayama T.
      • Ito Y.
      • Sakai K.
      • Miyake M.
      • Shibata E.
      • Ohno H.
      • Kamijima M.
      Comprehensive review of 2-ethyl-1-hexanol as an indoor air pollutant.
      ). Styrene is “generally recognized as safe” by the Food and Drug Administration, though it is documented to cause a “chemical/plastic/stale” off-flavor in most food products when at detectable levels (
      • Baner A.L.
      Case study: Styrene monomer migration into dairy products in single serve portion packs. Chapter 14.
      ). Acetophenone is permitted by the Food and Drug Administration for human consumption (21CFR172, 2019) but is known to cause “stale” off-flavors in milk (
      • Arnold R.G.
      • Libbey L.M.
      • Day E.A.
      Identification of components in the stale flavour fraction of sterilized concentrated milk.
      ) and 2-ethyl-1-hexanol has a “chemical/cleaning agent” aroma (
      • Drake M.A.
      • Miracle R.E.
      • Wright J.M.
      Sensory properties of dairy proteins. Chapter 16.
      ).
      Table 8Selected volatile compounds in skim milk packaged in glass, HDPE, PET, LDPE, LLDPE, and paperboard carton at 4°C across 4 time points (d 0, 5, 10, 15)
      Container
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Day of shelf life (time point)Volatile compound
      Volatile compounds were quantified using 6-point standard curves (R2 > 0.95).
      (μg/kg)
      StyreneTolueneP-xyleneP-cymeneHexanalOctanalNonanalLimoneneAcetophenone2-Ethyl-1-hexanol
      Glass0ND
      ND = not detected.
      0.51
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      NDND1.49
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.41
      ND0.57
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.37
      NDND3.09
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.16
      5ND1.11
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      NDND3.26
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.39
      ND1.96
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.30
      0.72
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      ND4.00
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.60
      10ND2.06
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.31
      0.46
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      ND10.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.32
      ND2.42
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.36
      0.98
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      ND4.11
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.05
      15ND1.26
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      0.40
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      ND22.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.99
      0.92
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.49
      3.86
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      0.88
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      0.51
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      4.16
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.35
      PET0ND0.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      NDND1.35
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      ND1.58
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      NDND1.05
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.29
      50.59
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      1.11
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      NDND5.23
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.52
      ND2.42
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      0.75
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      ND4.70
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.36
      10ND0.97
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      0.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      0.41
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      29.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.83
      0.79
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      13.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.12
      0.96
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      ND3.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.16
      151.18
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      1.68
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      0.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      0.42
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      21.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.69
      1.06
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.54
      4.58
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.26
      0.87
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      0.72
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      4.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.24
      HDPE0ND0.97
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      NDND2.5
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.63
      0.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      NDNDND3.72
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.25
      51.00
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      2.87
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.19
      0.89
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      0.85
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      12.8
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.73
      0.93
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.18
      4.28
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.63
      0.48
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      0.81
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      15.8
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.35
      100.99
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      2.42
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.20
      1.15
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      0.96
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      8.48
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.24
      0.87
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.21
      7.12
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.34
      1.30
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      1.07
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.20
      15.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.79
      151.40
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      4.51
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      1.92
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      1.71
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      15.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.29
      1.36
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.39
      4.69
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.13
      1.02
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      1.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      13.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.70
      LDPE0ND1.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      0.85
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      ND1.45
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.23
      ND1.52
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.28
      0.41
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      ND1.24
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.02
      50.77
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      2.53
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.40
      1.28
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.13
      0.47
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      8.89
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.20
      0.61
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      7.73
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      0.70
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.27
      0.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.13
      15.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.01
      101.17
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      1.91
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      1.39
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      0.68
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      1.61
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.32
      1.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.19
      2.05
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.16
      1.24
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      1.09
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.16
      24.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.38
      152.21
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      4.91
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.46
      3.16
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      0.91
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      12.9
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.12
      1.38
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.17
      10.8
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.56
      1.01
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      1.66
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.24
      19.8
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 3.55
      LLDPE01.75
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.39
      0.90
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      NDND1.61
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.62
      ND4.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.41
      0.40
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      ND1.06
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      58.93
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.71
      2.77
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      0.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      1.29
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      9.86
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.43
      0.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      19.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.97
      0.89
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      3.20
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.30
      36.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 4.52
      1012.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.09
      2.06
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.21
      1.38
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.18
      1.90
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      10.9
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.24
      ND9.03
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.89
      2.19
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.22
      5.65
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.35
      59.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 11.36
      1515.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.28
      5.29
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.26
      2.99
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      2.88
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      18.4
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 3.24
      0.52
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      34.6
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.98
      4.48
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      20.0
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.68
      119
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 12.34
      Carton00.74
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      0.53
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      NDND7.17
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.16
      ND0.70
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.38
      NDND2.81
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.54
      511.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.99
      3.42
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.60
      0.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      1.28
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      133
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 7.37
      4.71
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.50
      24.1
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.25
      0.80
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.24
      1.57
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.17
      60.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 5.90
      1013.9
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.09
      3.67
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.25
      1.27
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      1.87
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      115
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 9.25
      2.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.21
      17.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.41
      1.28
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      2.24
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      70.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 5.28
      1518.5
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.00
      4.90
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.43
      2.37
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.41
      2.72
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.16
      159
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 12.9
      2.22
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      28.8
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.48
      3.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.32
      6.54
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.34
      58.3
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.36
      a–m Means in the same column followed by a different superscript are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      2 Volatile compounds were quantified using 6-point standard curves (R2 > 0.95).
      3 ND = not detected.
      Table 9Volatile compounds in whole milk packaged in glass, HDPE, PET, LPDE, LLDPE, and paperboard carton at 4°C across 4 time points (d 0, 5, 10, 15)
      Container
      PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      Day of shelf life (time point)Volatile compound
      Volatile compounds were quantified using 6-point standard curves (R2 > 0.95).
      (μg/kg)
      StyreneTolueneP-xyleneP-cymemeHexanalOctanalNonanalLimoneneAcetophenone2-Ethyl-1-hexanol
      Glass0ND
      ND = not detected.
      0.82
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      NDND0.84
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.13
      ND2.85
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.93
      7.74
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.51
      NDND
      5ND0.98
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.26
      NDND1.30
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.17
      ND1.24
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.47
      3.04
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.70
      ND0.71
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.28
      10ND1.26
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      NDND1.38
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.46
      NDND2.67
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.41
      ND0.46
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      15ND1.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.23
      NDND3.26
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.89
      NDND2.32
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      ND0.41
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      PET0ND0.66
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      NDND2.98
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.30
      ND1.64
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.19
      7.45
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.23
      NDND
      5ND0.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      NDND3.17
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.65
      ND0.66
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.21
      1.69
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.22
      NDND
      10ND0.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      NDND4.11
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.54
      0.51
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.16
      0.75
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.34
      2.08
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      ND0.83
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.17
      15ND0.85
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      NDND2.28
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.06
      NDND1.65
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.13
      NDND
      HDPE0ND0.69
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      ND0.48
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.03
      2.14
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.14
      ND0.80
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      8.44
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.83
      NDND
      5NDNDNDND1.87
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.36
      NDND1.48
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.19
      NDND
      100.69
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      1.04
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.18
      NDND3.68
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.57
      ND0.64
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.34
      2.11
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.32
      ND0.52
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      15ND1.35
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      NDND4.89
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.34
      ND0.57
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      2.01
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      ND1.32
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      LDPE0ND1.67
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      NDND2.92
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.17
      ND1.52
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.16
      8.89
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.32
      NDND
      50.53
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      1.58
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.22
      NDND3.34
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.60
      ND0.60
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      2.56
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.39
      NDND
      100.63
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      3.64
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      0.42
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      ND3.94
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.22
      ND1.00
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.13
      3.01
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.40
      0.46
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      0.45
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      150.47
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.25
      5.09
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.31
      0.95
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.19
      ND2.57
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.18
      ND0.85
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.35
      3.12
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.40
      NDND
      LLDPE02.11
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      1.53
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      ND0.66
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      3.54
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.45
      0.36
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      2.27
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.26
      9.41
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.19
      0.98
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      0.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.19
      53.39
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.46
      1.18
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      0.40
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      0.51
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.08
      2.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.47
      ND2.67
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.49
      2.83
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.29
      2.64
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.40
      1.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      107.45
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.74
      2.05
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.33
      1.59
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.28
      0.94
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      4.09
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.42
      ND3.72
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.55
      3.78
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.80
      6.63
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.83
      8.06
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.22
      1512.7
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.38
      2.41
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.27
      2.16
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.20
      0.94
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.16
      3.54
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.54
      ND7.88
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.63
      4.83
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.58
      8.95
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.61
      7.67
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.60
      Carton00.71
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.02
      0.69
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      NDND6.85
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.39
      ND1.28
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.17
      10.8
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 4.49
      NDND
      50.88
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.01
      0.45
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.04
      NDND7.07
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.23
      ND0.87
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.28
      0.96
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      ND0.46
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.07
      102.66
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.32
      1.40
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.12
      0.60
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      ND13.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 1.77
      ND0.69
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.29
      1.86
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      1.05
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      0.69
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.11
      153.50
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.09
      1.43
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.06
      0.62
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.05
      ND19.2
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 2.19
      ND0.99
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.20
      2.13
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.10
      1.59
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.24
      1.71
      Means in the same column followed by a different superscript are significantly different (P < 0.05).
      ± 0.15
      a–m Means in the same column followed by a different superscript are significantly different (P < 0.05).
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene; LDPE = low-density polyethylene; LLDPE = linear low-density polyethylene. Values shown ± SD.
      2 Volatile compounds were quantified using 6-point standard curves (R2 > 0.95).
      3 ND = not detected.
      Limonene and p-cymene increased in skim milk packaged in LLDPE or paperboard cartons compared with LDPE, HDPE, PET, and glass. A similar increase was documented in whole milk packaged in LLDPE, but not paperboard carton (Tables 8, 9). Limonene is a compound found in fluid milk and is also used as a metric for assessment of environmental air quality (
      • Subramanian P.
      • Reynolds S.
      • Breuer G.
      • Whitten P.
      Analysis of limonene as an indicator for the assessment of environmental air-quality and health symptoms in office buildings.
      ). P-cymene is a degradation compound of limonene and correlates with the increase of limonene (R2 = 0.74, P < 0.05). The increased concentrations of limonene and the further degradation into p-cymene may be due to the higher permeability of LLDPE and paperboard cartons (Tables 3, 8, 9). By d 5 of storage, toluene and p-xylene concentrations were higher in skim milks packaged in HDPE, LDPE, LLDPE, and paperboard cartons compared with skim packaged in PET or glass (Table 8). Toluene and all 3 isomers of xylene (m-xlyene, o-xlyene, p-xlyene) are used as solvents or intermediates in the production of polymers for plastics (
      Agency for Toxic Substances and Disease Registry (ATSDR)
      Toxicological Profile for Toluene.
      ,
      Agency for Toxic Substances and Disease Registry (ATSDR)
      Toxicological Profile for Xylene.
      ). The degradation of plastics may be the cause for increased detection of toluene and p-xylene over time. These volatile compounds are likely contributors to stale/refrigerator flavor in skim and whole milks, consistent with
      • Lozano P.R.
      • Miracle E.R.
      • Krause A.J.
      • Drake M.A.
      • Cadwallader K.R.
      Effect of cold storage and packaging material on the major aroma components of sweet cream butter.
      who associated these compounds with refrigerator/stale flavor in stored butter.
      Skim and whole milks filled into paperboard cartons had increased levels of hexanal compared with milks in other packaging types (P < 0.05; Tables 8, 9). This difference in hexanal concentration may signify the degradation of the surrounding paperboard packaging, which may increase hexanal concentrations in milk across time. Hexanal is a degradation compound formed from the oxidation of linoleic acid during the processing of wood extracts or paper pulp, and has been associated with off-flavor in food products packaged in paperboard and cardboard (
      • Czerny M.
      • Buettner A.
      Odor-active compounds in cardboard.
      ).
      • Ziegleder G.
      Volatile and odorous compounds in unprinted paperboard.
      found that hexanal can be used as a guide for tracking lipid oxidation of paperboard across time. Both of these studies and others have found the association of increased concentrations of aldehydes such as hexanal, heptanal, octanal, nonanal, and decanal to an increase in cardboard/paperboard flavor (
      • Czerny M.
      • Buettner A.
      Odor-active compounds in cardboard.
      ;
      • Whitson M.E.
      • Miracle R.E.
      • Drake M.A.
      Sensory characterization of chemical components responsible for cardboard flavor in whey protein.
      ;
      • Van Caelenberg T.
      • Van Leuven I.
      • Dirinck P.
      An analytical approach for fast odour evaluation of recycled food-grade paperboard materials using HS-SPME-MS-nose technology.
      ). The elevated level of hexanal in milks in paperboard cartons compared with whole and skim milks in the other 5 packaging types may be attributed to hexanal being the largest source of paperboard or cardboard off-flavor (Tables 8, 9). Cardboard flavor in dairy products has also previously been sourced to lipid oxidation compounds including hexanal (
      • Whitson M.E.
      • Miracle R.E.
      • Drake M.A.
      Sensory characterization of chemical components responsible for cardboard flavor in whey protein.
      ). Skim milk packaged in paperboard cartons or LLDPE and whole milk packaged in LLDPE also had higher concentrations of nonanal compared with milks filled into other package types (P < 0.05; Table 8). This oxidation compound may indicate elevated lipid oxidation in milks packaged in these containers, indicative of increased oxygen and greater permeability of these package types compared with the others in this study. Dynamic headspace analysis conducted on commercial milks revealed that milks with higher concentrations of accessible oxygen were correlated with higher concentrations of lipid oxidation compounds (
      • Kim Y.D.
      • Morr C.V.
      Dynamic headspace analysis of light activated flavor in milk.
      ). Increases in lipid oxidation were also documented in ultrapasteurized and HTST milks stored at 4°C for 14 d (
      • Jo Y.
      • Benoist D.M.
      • Barbano D.M.
      • Drake M.A.
      Flavor and flavor chemistry differences among milks processed by high-temperature, short-time pasteurization or ultra-pasteurization.
      ).
      • Karatapanis A.E.
      • Badeka A.V.
      • Riganakos K.A.
      • Savvaidis I.N.
      • Kontominas M.G.
      Changes in flavour volatiles of whole pasteurized milk as affected by packaging material and storage time.
      also documented lipid oxidation in HTST whole milk across 7 d of refrigerated storage (
      • Karatapanis A.E.
      • Badeka A.V.
      • Riganakos K.A.
      • Savvaidis I.N.
      • Kontominas M.G.
      Changes in flavour volatiles of whole pasteurized milk as affected by packaging material and storage time.
      ). Concentrations of volatile compounds were generally lower in whole milk compared with skim milk (Tables 8, 9). This result may signify that packaging-related compounds may be soluble in milkfat and thus less available for headspace extraction, as well as sensory perception, because refrigerator/stale and cardboard flavors by sensory analysis were also less prevalent in whole milks compared with skim milks. In a previous study, fluid milk with higher fat content had lower recoveries of most volatile compounds via headspace extraction compared with skim milk (
      • Jo Y.
      • Benoist D.M.
      • Barbano D.M.
      • Drake M.A.
      Flavor and flavor chemistry differences among milks processed by high-temperature, short-time pasteurization or ultra-pasteurization.
      ).
      Skim milk packaged in paperboard carton or LLDPE were closely aligned across storage time with increasing concentrations of packaging-related volatile compounds styrene, toluene, p-xylene, p-cymene, limonene, acetophenone, 2-ethyl-1-hexanol, and hexanal (Figure 1). Skim milk packaged in LLDPE or paperboard carton after d 5 of storage were distinct from milks filled into glass, LDPE, HDPE, or PET (Figure 1). Whole milks were not as distinct in volatile compound differences compared with skim milks but whole milk packaged in LLDPE was closely aligned across storage time, with increasing concentrations of packaging-related volatile compounds styrene, toluene, p-xylene, p-cymene, acetophenone, 2-ethyl-1-hexanol, and hexanal (Figure 2). Whole milk packaged in LLDPE after d 5 was distinct from milk filled into glass, LDPE, HDPE, and PET (Figure 2).
      • Karatapanis A.E.
      • Badeka A.V.
      • Riganakos K.A.
      • Savvaidis I.N.
      • Kontominas M.G.
      Changes in flavour volatiles of whole pasteurized milk as affected by packaging material and storage time.
      reported slight increases in packaging/storage volatiles in the low part per billion (μg/kg) for toluene and various benzene compounds (ethyl benzene, 1,4-dimethyl benzene, 1,3-dimethyl benzene, 1,3,5-trimethyl benzene) in whole HTST milk across 7 d of refrigerated storage time.
      Figure thumbnail gr1
      Figure 1Principal component biplot of volatile compound concentrations and trained panel sensory attribute intensities of skim milk packaged in glass, high-density polyethylene (HDPE), polyethylene terephthalate (PET), low-density polyethylene (LDPE), linear low-density polyethylene (LLDPE), or paperboard cartons at 4°C across 4 time points (d 0, 5, 10, 15).
      Figure thumbnail gr2
      Figure 2Principal component biplot of volatile compound concentrations and trained panel sensory attribute intensities of whole milk packaged in glass, high-density polyethylene (HDPE), polyethylene terephthalate (PET), low-density polyethylene (LDPE), linear low-density polyethylene (LLDPE), or paperboard cartons at 4°C across 4 time points (d 0, 5, 10, 15).

      Relationship Between Volatile Compounds and Sensory Attributes

      Principal component biplots provide a visualization of the changes in the volatiles and sensory attributes that characterize differences among the milks with storage time, along with correlations between sensory attributes and volatile compounds (Figures 1, 2). All milks decreased in cooked/milky and sweet aromatic flavors with storage time, with milks filled into paperboard cartons or LLDPE having the most distinct storage-related changes (Figure 1). Paperboard/cardboard flavor and the volatile compound hexanal were correlated (R2 = 0.84, P < 0.05; Figure 1). Refrigerator stale flavor was correlated with the volatile compounds styrene (R2 = 0.88, P < 0.05), p-cymene (R2 = 0.73, P < 0.05), acetophenone (R2 = 0.71, P < 0.05), and 2-ethyl-1-hexanol (R2 = 0.82, P < 0.05). Styrene has also been attributed to “refrigerator/stale” flavor in stored butter (
      • Lozano P.R.
      • Miracle E.R.
      • Krause A.J.
      • Drake M.A.
      • Cadwallader K.R.
      Effect of cold storage and packaging material on the major aroma components of sweet cream butter.
      ). Similar to skim milk, whole milk filled into cartons was associated with paperboard/cardboard flavor and astringency, and the volatile compound hexanal (Figure 2). Whole milk packaged in LLDPE and paperboard cartons were also associated with higher overall aroma and refrigerator stale flavor compared with whole milk in other package materials (Figure 2).

      DSE-SAFE

      Skim and whole milk packaged in paperboard cartons had distinct paperboard and refrigerator stale flavors. Significant migration of volatile compounds associated with package or refrigerator taint compounds, or both, was documented at d 15 in whole and skim milk in paperboard cartons (whole and skim milk) by HS-SPME-GCMS (Tables 8, 9). Headspace volatile analysis of skim milk filled in paperboard cartons at d 15 revealed traces of plasticizers and ink-related compounds at detectable limits, but not above base line noise and, thus, not quantifiable (results not shown). Direct solvent extraction was subsequently applied to skim milks packaged in paperboard or glass to facilitate an increased recovery of these specific paperboard-related volatile compounds.
      Taint compounds associated with paperboard packaging that were recovered from skim milks by HS-SPME-GCMS were also recovered by DSE-SAFE, as expected. Relative abundance of package/refrigeration-related compounds in skim milk stored in paperboard cartons versus glass for 15 d by DSE-SAFE were (μg/kg) styrene 283 ± 23.3 versus 1.7 ± 0.7, toluene 930 ± 85.7 versus 8.6 ± 0.5, p-xylene 46 ± 12.4 versus 3.6 ± 2.1, acetophenone 38.7 ± 13.9 versus not detected, and 2-ethyl-1-hexanol 287.3 ± 8.19 versus 25.1 ± 0.90 (P < 0.05). These values were higher in skim milk by SAFE extraction compared with headspace extraction (SPME-GCMS) milk from paperboard carton, as expected, but it is also important to note that these values from SAFE extraction were calculated based on internal standard recovery rather than actual quantitation with standard curves (SPME-GCMS) and, thus, do not represent absolute concentrations. More importantly, as expected, additional compounds were detected in skim milk by DSE-SAFE extraction (but not by SPME-GCMS) including n-propylbenzene (dye solvent), benzophenone (photoiniators, plasticizer), dibutyl phthalate, and diethyl phthalate (plasticizers and softening agents;
      • Frostling H.
      • Hoff A.
      • Jacobsson S.
      • Pfäffli P.
      • Vainiotalo S.
      • Zitting A.
      Analytical, occupational and toxicologic aspects of the degradation products of polypropylene plastics.
      ;
      • Czerny M.
      • Buettner A.
      Odor-active compounds in cardboard.
      ;
      • Leja K.
      • Lewandowicz G.
      Polymer biodegradation and biodegradable polymers—A review.
      ). The compounds n-propylbenzene, benzophenone, dibutyl phthalate, and diethyl phthalate were detected in skim milk packaged in paperboard cartons at 210 ± 28.6, 6 ± 0.3, 56 ± 18, and 300 ± 1.79 parts per billion (μg/kg), respectively, via relative abundance. Comparatively, these compounds were not detected in skim milk packaged in glass. The compound n-propylbenzene is a compound found in plastic packaging for milk (
      • Abrantes S.
      GC-MS identification of styrene and oligomers in polystyrene plastic for milk packaging. Headspace and solution injection technique.
      ), and benzophenone was previously used as a photoinitiator to cure ink on packaging materials (21CFR172:
      • US FDA
      21CFR172: Food additives permitted for direct addition to food for human consumption.
      ; 21CFR177:
      • US FDA
      21CFR177: Indirect food additives: Polymers.
      ). Benzophenone and n-propylbenzene have been previously found in skim milk powders (
      • Shiratsuchi H.
      • Shimoda M.
      • Imayoshi K.
      • Noda K.
      • Osajima Y.
      Volatile flavor compounds in spray-dried skim milk powder.
      ;
      • Sanches-Silva A.
      • Pastorelli S.
      • Cruz J.M.
      • Simoneau C.
      • Castanheira I.
      • Paseiro-Losada P.
      Development of an analytical method for the determination of photoinitiators used for food packaging materials with potential to migrate into milk.
      ), possibly due to these products being packaged in plastic and using photoinitiators to cure the ink on the packages. Diethyl phthalate is classified as a plasticizer by the FDA for specific use in food (21CFR181.27:
      • US FDA
      21CFR181.27: Prior-sanctioned food ingredients.
      ). Dibutyl phthalate and diethyl phthalate have been previously found in milk and other dairy products (
      • Sørensen L.K.
      Determination of phthalates in milk and milk products by liquid chromatography/tandem mass spectrometry.
      ;
      • Fierens T.
      • Van Holderbeke M.
      • Willems H.
      • De Henauw S.
      • Sioen I.
      Transfer of eight phthalates through the milk chain—A case study.
      ). Determination of exact concentrations of these specific taint compounds in skim and whole milks filled into paperboard or other package types, or both, at larger package volumes merits future work.

      Consumer Difference Tests

      Consumers detected differences between skim and whole milk filled into paperboard cartons and glass (P < 0.05; Tables 10, 11). Consistent with minimal to no differences by trained panels, consumers could not detect differences between the PET and glass pairings and HDPE and glass (P > 0.05; Tables 10, 11). These results are consistent with a previous study, where consumers reported similar flavor and acceptability of milks packaged in PET and HDPE (
      • Potts H.L.
      • Amin K.N.
      • Duncan S.E.
      Retail lighting and packaging influence consumer acceptance of fluid milk.
      ). The results are also consistent with our trained descriptive panel and volatile compound results. Skim or whole milks filled into PET, HDPE or glass were similar in overall flavor profile with no paperboard/cardboard or refrigerator stale off-flavors present by d 15 (Tables 6, 7).
      Table 10Consumer tetrad difference tests of skim milk in 3 package types (paperboard, PET, and HDPE) compared with the control (glass) after 10 d at 4°C
      Packaging comparison
      PET = polyethylene terephthalate; HDPE = high-density polyethylene.
      nFluid skim milk
      No. correctSignificant at α = 0.05
      Paperboard carton vs. glass (control)10156Yes
      PET vs. glass (control)10328No
      HDPE vs. glass (control)10336No
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene.
      Table 11Consumer tetrad difference tests of whole milk in 3 package types (paperboard, PET, and HDPE) compared with the control (glass) after 10 d at 4°C
      Packaging comparison
      PET = polyethylene terephthalate; HDPE = high-density polyethylene.
      nFluid whole milk
      No. correctSignificant at α = 0.05
      Paperboard carton vs. glass (control)10447Yes
      PET vs. glass (control)10339No
      HDPE vs. glass (control)10337No
      1 PET = polyethylene terephthalate; HDPE = high-density polyethylene.
      Paperboard cartons are the most widely used packaging type for school meal (breakfast and lunch) programs in the United States, and milks packaged in paperboard cartons showed distinct off-flavors and presence of specific migration volatile compounds. Our findings suggest that industry and policy makers should seek new package alternatives for school meal fluid milk. The consequences of using fluid milk packaging that contributes significant off-flavors, over time, may affect how young children, and those children as adults, perceive fluid milk.

      CONCLUSIONS

      Off-flavors due to packaging and refrigeration environment were detected in fluid skim and whole milk over refrigerated storage time. Off-flavors in milks detected by sensory analysis were correlated with increased specific migration volatile compounds. Milks packaged in paperboard cartons and LLDPE had the highest intensities of off-flavors, due to permeability and migration, with off-flavors present by day 0 in skim milk. Consumer difference tests were consistent with trained panel and volatile compound analysis, suggesting that HTST milks packaged in HDPE, PET, or glass in the absence of light exposure have no discernable sensory differences by day 10 post-processing. In contrast, milks filled into paperboard cartons were differentiated by consumers compared to milks filled into glass.

      ACKNOWLEDGMENTS

      Funding provided in part by the National Dairy Council (Rosemont, IL). The use of tradenames does not imply endorsement nor lack of endorsement of those not mentioned. The authors have not stated any conflicts of interest.

      REFERENCES

        • Abrantes S.
        GC-MS identification of styrene and oligomers in polystyrene plastic for milk packaging. Headspace and solution injection technique.
        J. High Resolut. Chromatogr. 1993; 16: 113-115
        • Agency for Toxic Substances and Disease Registry (ATSDR)
        Toxicological Profile for Toluene.
        U.S. Public Health Service, U.S. Department of Health and Human Services, 2017
        • Agency for Toxic Substances and Disease Registry (ATSDR)
        Toxicological Profile for Xylene.
        U.S. Public Health Service, U.S. Department of Health and Human Services, 2007
        • Alvarez V.B.
        Fluid milk and cream products, Chapter 5.
        in: Clark S. Costello M. Drake M.A. Bodyfelt F. the Sensory Evaluation of Dairy Products. Springer, 2009: 73-134
        • Arnold R.G.
        • Libbey L.M.
        • Day E.A.
        Identification of components in the stale flavour fraction of sterilized concentrated milk.
        J. Food Sci. 1966; 31: 566-573
        • ASTM International
        E3009–15e1 Standard Test Method for Sensory Analysis—Tetrad Test.
        ASTM International, 2015
        • Azzara C.D.
        • Campbell L.B.
        Off-flavors of dairy products, Chapter 13.
        in: Charalambous G. Off-flavors in Foods and Beverages. Elsevier Science Publishing, 1992: 329-375
        • Bassette R.
        • Fung D.Y.C.
        • Mantha V.R.
        • Marth E.H.
        Off-flavors in milk.
        CRC Crit. Rev. Food Sci. Nutr. 1986; 24 (3519085): 1-52
        • Bradley Jr., R.L.
        Effect of light on alteration of nutritional value and flavor of milk: A review.
        J. Food Prot. 1980; 43 (30822864): 314-320
        • Brody A.L.
        Packaging milk and milk products, Chapter 20.
        in: Partridge J. Dairy Processing and Quality Assurance. 1st ed. Wiley-Blackwell, 2009: 443-464
        • Brothersen C.
        • McMahon D.J.
        • Legako J.
        • Martini S.
        Comparison of milk oxidation by exposure to LED and fluorescent light.
        J. Dairy Sci. 2016; 99 (26874418): 2537-2544
        • Crosby N.T.
        Plastic packaging materials, Chapter 1.
        in: Crosby N.T. Food Packaging Materials: Aspects of Analysis and Migration of Contaminants. 1st ed. Applied Science Publishers Ltd, 1981: 1-18
        • Czerny M.
        • Buettner A.
        Odor-active compounds in cardboard.
        J. Agric. Food Chem. 2009; 57: 9979-9984
        • Dainelli D.
        • Gontard N.
        • Spyropoulos D.
        • Zondervan-van den Beuken E.
        • Tobback P.
        Active and intelligent food packaging: legal aspects and safety concerns.
        Trends Food Sci. Technol. 2008; 19: S103-S112
        • Drake M.A.
        • Miracle R.E.
        • Wright J.M.
        Sensory properties of dairy proteins. Chapter 16.
        in: Thompson A. Boland M. Singh H. Milk Proteins: From Expression to Food. 2nd ed. Academic Press, 2014: 473-492
        • Engel W.
        • Bahr W.
        • Schieberle P.
        Solvent assisted flavour evaporation—A new and versatile technique for the careful and direct isolation of aroma compounds from complex food matrices.
        Eur. Food Res. Technol. 1999; 209: 237-241
        • Fierens T.
        • Van Holderbeke M.
        • Willems H.
        • De Henauw S.
        • Sioen I.
        Transfer of eight phthalates through the milk chain—A case study.
        Environ. Int. 2013; 51 (23138015): 1-7
        • Frostling H.
        • Hoff A.
        • Jacobsson S.
        • Pfäffli P.
        • Vainiotalo S.
        • Zitting A.
        Analytical, occupational and toxicologic aspects of the degradation products of polypropylene plastics.
        Scand. J. Work Environ. Health. 1984; 10 (6474110): 163-169
        • Jo Y.
        • Benoist D.M.
        • Barbano D.M.
        • Drake M.A.
        Flavor and flavor chemistry differences among milks processed by high-temperature, short-time pasteurization or ultra-pasteurization.
        J. Dairy Sci. 2018; 101 (29501345): 3812-3828
        • Johnson D.S.
        • Duncan S.E.
        • Bianchi L.M.
        • Chang H.H.
        • Eigel W.N.
        • O'Keefe S.F.
        Packaging modifications for protecting flavor of extended-shelf-life milk from light.
        J. Dairy Sci. 2015; 98 (25660743): 2205-2214
        • Karatapanis A.E.
        • Badeka A.V.
        • Riganakos K.A.
        • Savvaidis I.N.
        • Kontominas M.G.
        Changes in flavour volatiles of whole pasteurized milk as affected by packaging material and storage time.
        Int. Dairy J. 2006; 16: 750-761
        • Kim Y.D.
        • Morr C.V.
        Dynamic headspace analysis of light activated flavor in milk.
        Int. Dairy J. 1996; 6: 185-193
        • Kim-Kang H.
        Volatiles in packaging materials.
        Crit. Rev. Food Sci. Nutr. 1990; 29 (2257079): 255-271
        • Lee A.P.
        • Barbano D.M.
        • Drake M.A.
        The influence of ultra-pasteurization by indirect heating versus direct steam injection on skim and 2% fat milks.
        J. Dairy Sci. 2017; 100 (28088421): 1688-1701
        • Leja K.
        • Lewandowicz G.
        Polymer biodegradation and biodegradable polymers—A review.
        Pol. J. Environ. Stud. 2010; 19: 255-266
        • Leong C.M.O.
        • Harte B.R.
        • Partridge J.A.
        • Ott D.B.
        • Downes T.W.
        Off-flavor development in milk packaged in polyethylene-coated paperboard cartons.
        J. Dairy Sci. 1992; 75: 2105-2111
        • López P.
        • Batlle R.
        • Salafranca J.
        • Nerín C.
        Efficiency of whole and skimmed powdered milk for trapping volatile compounds released from plastic containers in high-temperature applications.
        J. Food Prot. 2008; 71 (18810874): 1889-1897
        • Lozano P.R.
        • Miracle E.R.
        • Krause A.J.
        • Drake M.A.
        • Cadwallader K.R.
        Effect of cold storage and packaging material on the major aroma components of sweet cream butter.
        J. Agric. Food Chem. 2007; 55 (17705437): 7840-7846
        • Lozano P.R.
        • Drake M.A.
        • Benitez D.
        • Cadwallader K.R.
        Instrumental and sensory characterization of heat-induced odorants in aseptically packaged soy milk.
        J. Agric. Food Chem. 2007; 55 (17373812): 3018-3026
        • McCarthy K.S.
        • Lopetcharat K.
        • Drake M.A.
        Milk fat threshold determination and the effect of milk fat content on consumer preference for fluid milk.
        J. Dairy Sci. 2017; 100 (28088417): 1702-1711
        • Meilgaard M.
        • Civille G.
        • Carr B.
        Descriptive analysis techniques. Chapter 10.
        in: Meilgaard M. Civille G. Carr B. Sensory Evaluation Techniques. 4th ed. CRC Press, 2007: 173-186
        • Nielsen T.J.
        • Jägerstad I.M.
        • Öste R.E.
        Study of factors affecting the absorption of aroma compounds into low-density polyethylene.
        J. Sci. Food Agric. 1992; 60: 377-381
        • NIST (National Institute of Standards and Technology)
        NIST Wiley Registry: NIST Mass Spectral Library.
        11th ed. Wiley, 2017
        • Northrop-Clewes C.A.
        • Thurnham D.I.
        The discovery and characterization of riboflavin.
        Ann. Nutr. Metab. 2012; 61 (23183293): 224-230
        • Baner A.L.
        Case study: Styrene monomer migration into dairy products in single serve portion packs. Chapter 14.
        in: Piringer O.G. Baner A.L. Plastic Packaging Materials for Food: Barrier Function, Mass Transport, Quality Assurance and Legislation. Wiley-VCH Verlag GmbH, 2000: 427-443
        • Potts H.
        • Duncan S.
        • Johnson M.L.
        • O'Keefe S.F.
        • Marcy J.E.
        • Mallikarjunan K.
        0712 Oxygen barrier and light interference packaging properties for controlling light-induced oxidation in milk.
        J. Anim. Sci. 2016; 94: 341
        • Potts H.L.
        • Amin K.N.
        • Duncan S.E.
        Retail lighting and packaging influence consumer acceptance of fluid milk.
        J. Dairy Sci. 2017; 100 (27865502): 146-156
        • Roland A.M.
        • Hotchkiss J.H.
        Determination of flavor—Polymer interactions by vacuum-microgravimetric method. Chapter 13.
        in: Risch S.J. Hotchkiss J.H. Food and Packaging Interactions II. American Chemical Society, 1991: 149-160
        • Sajilata M.G.
        • Savitha K.
        • Singhal R.S.
        • Kanetkar V.R.
        Scalping of flavors in packaged foods.
        Compr. Rev. Food Sci. Food Saf. 2007; 6: 17-35
        • Sanches-Silva A.
        • Pastorelli S.
        • Cruz J.M.
        • Simoneau C.
        • Castanheira I.
        • Paseiro-Losada P.
        Development of an analytical method for the determination of photoinitiators used for food packaging materials with potential to migrate into milk.
        J. Dairy Sci. 2008; 91 (18292245): 900-909