J. Dairy Sci. 89:3690-3695
© American Dairy Science Association, 2006.
Estimating Fatty Acid Content in Cow Milk Using Mid-Infrared Spectrometry
H. Soyeurt*,
,1,2,
P. Dardenne
,
F. Dehareng
,
G. Lognay*,2,
D. Veselko
,
M. Marlier*,2,
C. Bertozzi#,
P. Mayeres*,#,2 and
N. Gengler*,||,2
* Gembloux Agricultural University, B-5030 Gembloux, Belgium
Fonds pour la Formation à la Recherche dans lIndustrie et lAgriculture (FRIA), B-1000 Brussels, Belgium
Walloon Agricultural Research Centre, Quality Department, B-5030 Gembloux, Belgium
Milk Committee, B-4651 Battice, Belgium
# Walloon Breeders Association, B-5530 Ciney, Belgium
|| National Fund for Scientific Research, B-1000 Brussels, Belgium
1 Corresponding author: soyeurt.h{at}fsagx.ac.be
 |
ABSTRACT
|
|---|
Interest in the fatty acid composition of dairy products is increasing; however, the measurement of fatty acids requires using gas-liquid chromatography. Although this method is suitable, it involves a time-consuming procedure, expensive reagents, and qualified staff. By comparison, the mid-infrared (MIR) spectrometry method could be a good alternative for assessing the fatty acid profile of dairy products. The objective of this study was to explore the calibration of MIR spectrometry for estimating fatty acid concentrations in milk and milk fat. Estimated concentrations in milk fat were less reliable than those for the same fatty acids in milk. Results also showed that when the fatty acid concentrations in milk increased, the efficiency of the infrared analysis method in predicting these values simultaneously increased. Selected prediction equations must have a high cross-validation coefficient of determination, a high ratio of standard error of cross-validation to standard deviation, and good repeatability of chromatographic data. Results from this study showed that the calibration equations predicting 12:0, 14:0, 16:0, 16:1cis-9, 18:1, and saturated and monounsaturated fatty acids in milk could be used. Thus, with its potential for use in regular milk recording, this infrared analysis method offers the possibility of assessing and improving the quality of milk produced. Indeed, it enables the fatty acid composition in milk to be estimated for each cow and the estimates to be used as indicator traits to determine the genetic values of underlying fatty acid concentrations. The knowledge of these genetic values would open up opportunities for animal selection aimed at improving the nutritional quality of cow milk.
Key Words: milk fatty acid mid-infrared quality
 |
INTRODUCTION
|
|---|
Cow milk fat typically contains 70% saturated fatty acids (SAT), 25% monounsaturated fatty acids (MONO), and 5% polyunsaturated fatty acids (POLY; Grummer, 1991). A milk lipid composition more favorable to human health would be about 30% SAT (Pascal, 1996), 60% MONO, and 10% POLY (Hayes and Kosla, 1992). The fatty acid profile of cow milk is therefore far from optimal. However, the observed variations in SAT, MONO, and POLY suggest that the milk fat composition could be modified by various means, such as through feeding and genetics (Palmquist et al., 1993), and could come closer to the optimal profile. Several researchers have been focusing on ways to improve the nutritional quality of bovine milk fat by feed supplementation (e.g., Demeyer and Doreau, 1999; Chilliard et al., 2000). However, all the fatty acids in a specific class (SAT, MONO, or POLY) do not have the same effects on human health. In the case of SAT, although myristic acid is known for its negative effects on cardiovascular diseases, stearic acid does not seem to have this effect (Hu et al., 1999). Similarly, in POLY, the n-6 fatty acids appear to have negative effects on human health because of their prevalence in Western nutrition. Indeed, the current ratio n-6:n-3 is estimated to be 15:1 to 20:1 (Simopoulos, 2003). It is therefore important to check the global fatty acid profile in milk if one wants to assess the nutritional quality of bovine milk fat. As stated earlier, influencing the nutritional quality of milk fat has been the topic of several recent research papers (e.g., Demeyer and Doreau, 1999; Chilliard et al., 2000), some of which have concentrated on milk fat enriched in n-3, in line with the current interest in dairy products that are enriched in n-3. However, there has been less focus on assessing the fatty acid content, because this requires chromatographic analysis. Although this method is suitable (e.g., Dorey et al., 1988; Collomb and Bühler, 2000), it is time-consuming and requires skilled staff.
Mid-infrared (MIR) spectrometry is an alternative to gas chromatography, with advantages such as a very high throughput (up to 500 samples/h; FOSS, 2005), ease of use, and availability. The infrared spectrum is caused by the absorptions of electromagnetic radiation at frequencies that are correlated to the vibrations of specific chemical bonds within a molecule (Coates, 2000). The spectrum therefore illustrates these absorptions at different wavenumbers (cm1) for a specific chemical composition (Smith, 1996). Mid-infrared spectrometry (400 to 4,000 cm1) is particularly interesting because it is very highly sensitive to the chemical environment, as the fundamental absorptions of molecular vibrations occur in this region (Belton, 1997). Mid-infrared spectrometry can be used to estimate various traits quantitatively based on calibration equations. The purpose of our research was to develop the predicted equations necessary for measuring the fatty acid content in milk and milk fat using MIR spectrometry.
 |
MATERIALS AND METHODS
|
|---|
Sampling and Recording Spectra Files
Milk samples were taken from cows in 7 herds selected according to the following criteria: their participation in the milk recording system in Wallonia, the observed variation in the percentage of milk fat and the number of breeds in the herds. An 80-mL sample of milk taken during routine milk recording was divided into 2 parts (60 + 20 mL). The samples were collected for all the cows milked in the herds on a given test day. Following standard procedures (International Committee for Animal Recording, 2004), the samples represented 50% morning milk and 50% evening milk. The 20-mL sample was then analyzed using MIR spectrometry (MilkoScan FT6000; FOSS, 2005) following the normal milk recording procedure (International Committee for Animal Recording, 2004). The MilkoScan FT6000 works within the MIR region from 1,000 to 5,000 cm1 and uses an interferometer. From the resulting interferogram, MIR spectra are generated by means of fast Fourier transformations (FOSS, 2005). The spectra files generated were recorded in a database. The second sample (60 mL) was frozen at 26 ± 2°C. In this procedure, 600 samples were taken between April and June 2005 from 275 cows from 6 breeds (Dual Purpose Belgian Blue, Holstein-Friesian, Jersey, Normande, Montbeliarde, and Red and White). Not all the farms were tested 3 times, and some cows were dried off or had calved during the study. The milk fat percentage of the collected samples ranged between 2.97 and 7.73 g/dL of milk. This large variation indicated that a good calibration could be made.
Reference Values
Using principal component analysis based on the spectral variability, 49 samples were chosen and used from the 600 samples collected. The milk fat was extracted according to ISO Standard 14156:2001 (International Organization for Standardization, 2001). These milk fat samples were analyzed using gas chromatography, based on a method derived from Collomb and Bühler (2000). The gas chromatograph (model 6890; Agilent Technologies, Inc., Palo Alto, CA) was equipped with a CPSil-57 CB capillary column (Varian, Inc., Palo Alto, CA), with a length of 50 m, an internal diameter of 0.25 mm, and a film thickness of 0.20 µm. The retention gap was a "methyl deactivated nonpolar" (Varian, Inc.) with a length of 20 cm and an internal diameter of 0.53 mm. The conditions for the chromatographic analyses were as follows: carrier gas, helium; average velocity, 35 cm/s; cold on-column injector; flame ionization detector at 265°C; and a temperature program from 40°C (2 min) to 150°C (at 30°C/min), then 150 to 250°C (at 2°C/min). The volume injected was 0.5 µL. To measure the fatty acid concentrations, the response factors used were the same as those described by Collomb and Bühler (2000) because the experimental conditions were similar.
For each sample, 2 groups of reference values were generated (g of fatty acid/100 g of fat, and g of fatty acid/dL of milk) and were recorded in the database. The repeatability and accuracy of this method were estimated from several reference butter samples produced by BIPEA (Bureau InterProfessionnel dEtude Analytique, http://www.bipea.org).
Calibration Equations
From the chromatographic and spectral data, a specific program for multivariate calibration (WINISI III; http://www.winisi.com/) was used to compute the calibration equations using partial least squares regression (PLS). Partial least squares regression has 2 important advantages over multiple linear regression or regression on principal components. First, like principal components regression, it uses all the spectral data for the calibration (Frank et al., 1984). Second, in one step, the PLS method compresses the data (Martens and Jensen, 1982) and maximizes the variability of the dependent variable (Martens and Naes, 1987). The PLS method is therefore considered more efficient for calibration than the regression on principal components or multiple linear regression (Prévot, 2004). The number of factors used in the equation was determined by cross-validation, which was also needed to estimate its robustness. Although various pretreatments were studied, none was used on the spectral data before the calibration. Finally, 41 prediction equations were elaborated to estimate the fatty acid profile in milk and milk fat. To assess the efficiency of the calibration equations, various statistical parameters were estimated and analyzed: mean, standard deviation, standard error of calibration (SEC), calibration coefficient of determination, standard error of cross-validation (SECV), and cross-validation coefficient of determination (RCV2). The ratio of SECV to standard deviation (RPD) was also calculated (Williams and Norris, 2001) to assess the efficiency of the calibration.
 |
RESULTS AND DISCUSSION
|
|---|
Reference Values
Figure 1
presents a typical gas chromatogram recorded under the conditions described earlier. Because of the poor resolution between 18:1 isomers, they were not studied individually. Thus, the concentration indicated for 18:1 in this study is the result of the sum of the different 18:1 isomer contents.
The accuracy of the chromatographic method was assessed by comparing the mean values obtained with the reference values set by BIPEA (Table 1
). For all fatty acids, with the exception of 4:0 and 18:1cis-9, the mean results obtained were in the variation limits of the reference values. The divergence for 4:0 with the values set by BIPEA might be explained by the more cautious procedure for volatile fatty acids followed during the experiment: The Sovirel tubes (VWR International, Leuven, Belgium) were closed by a hermetic cap between all manipulations. The divergence for 18:1cis-9 might be explained by the poor chromatographic resolution between this fatty acid and other 18:1 isomers.
Table 1
shows the coefficients of variation for the butter samples analyzed by BIPEA. In this study, repeatability for a fatty acid was considered acceptable if its coefficient of variation was less than or equal to 5%. With good repeatability obtained for most of the fatty acids studied, 2 decimal places were used in Table 1
. In this study, 6:0 and 8:0 showed poor repeatability, which was attributed to some losses owing to their relative volatility. Nevertheless, when pooled in the group of SAT, the repeatability was acceptable. Generally, good repeatability was observed for fatty acids with longer carbon chains. Overall, except for short-chain fatty acids, the accuracy of the chromatographic method used was acceptable.
Calibration Equations
The applied PLS analysis resulted in equations with approximately 10 factors combining more than 500 values in each equation.
The potential for estimating the fatty acid composition in milk using MIR spectrometry might be explained by the absorptions of electromagnetic radiation at frequencies that are correlated to the vibrations of specific chemical bonds within molecules (Coates, 2000). This explanation is easy with a simple matrix as a mix of 2 different components, but the milk matrix is very complex. The spectrum is therefore the result of successive interactions due to the chemical bonds from all the constituents (fatty acids, proteins, lactose, etc.). By comparing the milk spectrum and specific fatty acid spectrum, the principal MIR regions that were implicated in estimating the fatty acid profile were located between 1,736 and 1,805 cm1 and between 2,823 and 3,016 cm1. The implication of the first region is logical because Coates (2000) indicated that 1,745 cm1 is the frequency correlated with the vibration of the fatty acid carbonyl group.
Table 2
shows the estimated statistical parameters for each calibration equation. These show that the correlations for predicting fatty acid concentrations in milk were better than those for predicting the same fatty acid concentrations in milk fat. This might be explained by a different dispersion of values obtained for concentrations of fatty acids in milk or milk fat. Indeed, 2 milk samples can have the same fat profile but different percentages of fat in the milk. The profile values expressed in grams per deciliter of milk were autocorrelated more closely than those expressed in grams per 100 grams of fat. The reference values used to establish the predicted equation for fat in milk came from the predicted values obtained using the Milkoscan FT6000. This explains the high result of R2CV obtained for this calibration equation.
View this table:
[in this window]
[in a new window]
|
Table 2. Estimated statistical parameters for each calibration equation that estimated the concentrations of fatty acid in milk (g/dL of milk) and in milk fat (g/100 g of fat)1
|
|
The SEC parameter underestimates the mean square error of the model because the residual variance is not taken into account. Therefore, SEC neglects the variance of regression coefficients and, in the context of calibration model validation, SECV is preferred (Prévot, 2004). To estimate the efficiency of the calibration, the RPD was calculated (Table 2
, Figure 2
). For potential use, high R2CV and high RPD parameters would be required. Generally, if the value of R2CV was high, the value of RPD was high (Table 2
).

View larger version (9K):
[in this window]
[in a new window]
|
Figure 2. Variation trend of RPD, the ratio of the standard error of cross-validation (SECV) to the standard deviation (SD) as a function of the concentration of fatty acids in milk.
|
|
Based on the results in Table 2
, our study also showed that there was a second-order polynomial relationship (y = 0.0374x2 + 0.7301x + 1.3446) between the fatty acid concentrations in milk and the value of RPD (r2 = 0.75). Generally, therefore, if the concentration of fatty acid is high, the potential for predicting the concentration using MIR spectrometry analysis is also high.
To verify whether the predicted concentrations of fatty acids obtained by the calibration equations were due to real absorbance of these fatty acids or only to the correlations between the total fat content and fatty acids, the correlations between total fat and the studied fatty acids were calculated (Table 3
). If the calibration correlations (RCV) were not due to real absorbances specific to fatty acids, these correlations would not be higher than the correlations between total fat and fatty acids. Thus, the predicted concentrations for these fatty acids resulted more from a real infrared prediction than the correlation with the fat. It is interesting to observe that the differences between the correlation with fat and RCV for short-chain fatty acids (4:0, 6:0, 8:0) were higher than the others. Therefore, specific spectral information can be extracted by the PLS model independently from the correlation with total fat.
View this table:
[in this window]
[in a new window]
|
Table 3. Correlations between the percentage of milk fat and different concentrations of studied fatty acids in milk1
|
|
 |
CONCLUSIONS
|
|---|
The estimation of fatty acid concentrations in milk and in milk fat using MIR spectrometry seems feasible. Results from this study showed that the calibration equations predicting 10:0, 12:0, 14:0, 16:0, 16:1cis-9, 18:1, 18:2cis-9,cis-12, SAT, and MONO in milk could be used. Although some fatty acids present in low concentrations in milk could not be predicted accurately (e.g., n-3 and 14:1cis-9), MIR spectrometry could predict most fatty acids (e.g., 14:0, 16:0, 18:0, 18:1, SAT, and MONO). Indeed, most associated coefficients of determination were clearly lower than 1, but they were significantly different from zero. Therefore, and because of the speed of analysis and the present application of the methodology in routine milk testing, MIR spectrometry is an important alternative in the dairy sector for providing indications of the fatty acid profiles in cow milk. This option seems particularly useful for application in milk recording schemes. The composition of milk fat for each milk sample from an animal could be estimated using the calibration equations established in this study. The estimates based on MIR could be used as indicator traits for real underlying fatty acid concentrations, potentially in a multitrait setting using the appropriate selection index theory. Knowledge of these genetic values would open up opportunities for animal selection aimed at improving the nutritional quality of cow milk.
 |
ACKNOWLEDGEMENTS
|
|---|
Hélène Soyeurt acknowledges the support of the FRIA through a grant scholarship. Nicolas Gengler, a Research Associate of the National Fund for Scientific Research (Brussels, Belgium), acknowledges the Funds support. The authors wish to thank Danny Trisman for his laboratory work and to acknowledge the technical support provided by the Walloon Breeding Association (AWE), the Walloon Milk Committee, and the Walloon Agricultural Research Centre. The partial financial support provided by the Walloon Regional Ministry of Agriculture (Ministère de la Région Wallonne, Direction Générale de lAgriculture, Namur, Belgium) is also acknowledged.
 |
FOOTNOTES
|
|---|
2 H. Soyeurt, P. Mayeres, and N. Gengler are with the Animal Science Unit, G. Lognay is with the Analytical Chemistry Unit, and M. Marlier is with the General and Organic Chemistry Unit of Gembloux Agricultural University. 
Received for publication December 23, 2005.
Accepted for publication March 14, 2006.
 |
REFERENCES
|
|---|
Belton, P. S. 1997. Spectroscopic approaches to the measurement of food quality. Pure Appl. Chem. 69:4750.
Chilliard, Y., A. Ferlay, R. M. Mansbridge, and M. Doreau. 2000. Ruminant plasticity: Nutritional control of saturated, polyunsaturated, trans and conjugated fatty acids. INRA Ann. Zootechnol. 49:181205.
Coates, J. 2000. Interpretation of infrared spectre, a practical approach. Pages 1081510837 in Encyclopedia of Analytical Chemistry. R. A. Meyers, ed. John Wiley & Sons, New York, NY.
Collomb, M., and T. Bühler. 2000. Analyse de la composition en acides gras de la graisse de lait. Mitt. Lebensm. Hyg. 91:306332.
Demeyer, D., and M. Doreau. 1999. Pourquoi et comment modifier les lipides du lait et de la viande bovine? Cah. Nutr. Diet. 34:301308.
Dorey, F., D. Brodin, J.-F. Le Querler, and S. Kudzal-Savoie. 1988. Analyse des acides gras du beurre par chromatographie en phase gazeuse couplée avec la spectrométrie de masse. IAA (June):437441.
FOSS. 2005. MilkoScanTM FT6000. http://www.foss.dk/c/p/solutions/products/showprodfamily.asp?prodfamilypkid=81 Accessed Nov. 17, 2005.
Frank, I. E., J. Feikema, N. Constantine, and B. R. Kowalski. 1984. Prediction of product quality from spectral data using the partial least-squares method. J. Chem. Inf. Comput. Sci. 24:2024.
Grummer, R. R. 1991. Effect of feed on the composition of milk fat. J. Dairy Sci. 74:32443257.[Abstract]
Hayes, K. C., and D. R. Khosla. 1992. Dietary fatty acid thresholds and cholesterolemia. FASEB J. 6:26002607.[Abstract]
Hu, F. B., M. J. Stampfer, J. E. Manson, A. Ascherio, G. A. Colditz, F. E. Speizer, C. H. Hennekens, and W. C. Willet. 1999. Dietary saturated fat and their food sources in relation to the risk of coronary heart disease in women. Am. J. Clin. Nutr. 70:10011008.[Abstract/Free Full Text]
International Committee for Animal Recording (ICAR). 2004. Section 2. ICAR rules, standards and guidelines for dairy production recording. http://www.icar.org/docs/Rules%20and%20regulations/Guidelines/Guidelines_2005_final_low_resolution.pdf Accessed Nov. 17, 2005.
International Organization for Standardization (ISO). 2001. Lait et produits laitiersMéthodes dextraction des lipides et des composés liposolubles. ISO 14156:2001. FIL-IDF 172:2001. ISO, Geneva, Switzerland.
Martens, H., and S. A. Jensen. 1982. Partial least squares regression: A new two-stage NIR calibration method. Pages 607647 in Proc. 7th World Cereal and Bread Congr. J. Holas and J. Kratochvil, ed. Elsevier, Amsterdam.
Martens, H., and T. Naes. 1987. Multivariate Calibration by Data Compression. John Wiley & Sons, Chichester, UK. 420 pp.
Palmquist, D. L., A. D. Beaulieu, and D. M. Barbano. 1993. Feed and animal factors influencing milk fat composition. J. Dairy Sci. 76:17531771.[Abstract]
Pascal, G. 1996. Les apports quotidiens recommandés en lipids et en acides gras. OCL 3:205210.
Prévot, H. 2004. Comparaison de méthodes statistiques et neuronales pour létablissement déquations de calibrage en spectrométrie de réflexion diffuse dans le proche infrarouge. Ph.D. Thesis, Gembloux Agricultural University, Belgium. 382 pp.
Simopoulos, A. P. 2003. Importance of the ratio of omega-6/omega-3 essential fatty acids: Evolutionary aspects. World Rev. Nutr. Diet. 92:122.[Medline]
Smith, B. C. 1996. Fundamentals of Fourier Transform Infrared Spectroscopy. CRC Press, Boca Raton, FL.
William, P., and K. Norris. 2001. Near-Infrared Technology in the Agricultural and Food Industries. American Association of Cereal Chemists, St. Paul, MN.
This article has been cited by other articles:

|
 |

|
 |
 
H. Soyeurt, P. Dardenne, F. Dehareng, C. Bastin, and N. Gengler
Genetic Parameters of Saturated and Monounsaturated Fatty Acid Content and the Ratio of Saturated to Unsaturated Fatty Acids in Bovine Milk
J Dairy Sci,
September 1, 2008;
91(9):
3611 - 3626.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
H. Soyeurt, F. Dehareng, P. Mayeres, C. Bertozzi, and N. Gengler
Variation of {Delta}9-Desaturase Activity in Dairy Cattle
J Dairy Sci,
August 1, 2008;
91(8):
3211 - 3224.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
W. M. Stoop, J. A. M. van Arendonk, J. M. L. Heck, H. J. F. van Valenberg, and H. Bovenhuis
Genetic Parameters for Major Milk Fatty Acids and Milk Production Traits of Dutch Holstein-Friesians
J Dairy Sci,
January 1, 2008;
91(1):
385 - 394.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
H. Soyeurt, A. Gillon, S. Vanderick, P. Mayeres, C. Bertozzi, and N. Gengler
Estimation of Heritability and Genetic Correlations for the Major Fatty Acids in Bovine Milk
J Dairy Sci,
September 1, 2007;
90(9):
4435 - 4442.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
H. Soyeurt, F. G. Colinet, V. M.-R. Arnould, P. Dardenne, C. Bertozzi, R. Renaville, D. Portetelle, and N. Gengler
Genetic Variability of Lactoferrin Content Estimated by Mid-Infrared Spectrometry in Bovine Milk
J Dairy Sci,
September 1, 2007;
90(9):
4443 - 4450.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. Moioli, G. Contarini, A. Avalli, G. Catillo, L. Orru, G. De Matteis, G. Masoero, and F. Napolitano
Short Communication: Effect of Stearoyl-Coenzyme A Desaturase Polymorphism on Fatty Acid Composition of Milk
J Dairy Sci,
July 1, 2007;
90(7):
3553 - 3558.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
H. Soyeurt, P. Dardenne, A. Gillon, C. Croquet, S. Vanderick, P. Mayeres, C. Bertozzi, and N. Gengler
Variation in Fatty Acid Contents of Milk and Milk Fat Within and Across Breeds
J Dairy Sci,
December 1, 2006;
89(12):
4858 - 4865.
[Abstract]
[Full Text]
[PDF]
|
 |
|