|
|
||||||||
,1


* Gembloux Agricultural University, Animal Science Unit, B-5030 Gembloux, Belgium
Fonds pour la Formation à la Recherche dans lIndustrie et lAgriculture (FRIA), B-1000, Brussels, Belgium
Walloon Agricultural Research Center, Quality Department, B-5030 Gembloux, Belgium
National Fund for Scientific Research, B-1000 Brussels, Belgium
¶ Walloon Breeders Association, B-5530 Ciney, Belgium
1 Corresponding author: soyeurt.h{at}fsagx.ac.be
| ABSTRACT |
|---|
|
|
|---|
Key Words: milk quality breed difference fatty acid fat
| INTRODUCTION |
|---|
|
|
|---|
Average cow milk fat contains 70% saturated fatty acids (SAT), 25% monounsaturated fatty acids (MONO), and 5% polyunsaturated fatty acids (POLY) (Grummer, 1991). A combination of milk lipids more favorable to human health would be around 30% SAT (Pascal, 1996), 60% MONO, and 10% POLY (Hayes and Khosla, 1992). The fatty acid profile of cows milk is therefore far from optimal. However, the observed variations in SAT, MONO, and POLY suggest that the milk fat composition can be modified by various means (e.g., feeding and genetics) to be brought closer to the optimal profile (Palmquist et al., 1993). Many results are available describing the influence of feeding (e.g., Demeyer and Doreau, 1999; Chilliard et al., 2000); however, many fewer are available on breed and individual genetic differences, even though variation in the overall milk fat among breeds is well known (e.g., Interbull, 2005). Also, not all the fatty acids in a specific class (SAT, MONO, or POLY) have the same effect on human health. In the case of SAT, although myristic acid is known for its negative effect 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 overabundance in the Western diet. The current ratio of n-6:n-3 is estimated to be 1520 to 1 in the human Western diet (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.
An alternative way to study the nutritional quality of milk fat is to analyze the variations in
9-desaturase activity. This enzyme is important to the fatty acid composition in milk because it is responsible for the majority of MONO and the totality of conjugated linoleic acids (CLA) in milk. The
9-desaturase activity can be studied by analyzing the ratio of product to substrate (14:1 cis-9 to 14:0, 16:1 cis-9 to 16:0, 18:1 cis-9 to 18:0, and 18:1 trans-11 to CLA). Through this methodology, Lock and Garnsworthy (2003) studied the seasonal variation in the activity of this enzyme and Peterson et al. (2002) showed individual animal differences in
9-desaturase activity.
Studying the genetics of the nutritional quality of milk requires a large amount of data. Many previous studies have used chromatographic analyses to estimate the fatty acid content in milk fat (e.g., Bobe et al., 1994; DePeters et al., 1995). This method is efficient (e.g., Dorey et al., 1988; Collomb and Bühler, 2000) but requires a lengthy analysis, expensive reagents, and highly skilled staff. Therefore, these studies are generally restricted in the number of animals and samples available (e.g., 234 cows for Bobe et al., 1994; 90 cows for DePeters et al., 1995). Mid-infrared (MIR) spectrometry is a faster method of estimating the different components of milk (up to 500 samples/h; Foss, 2005). This technology is currently integrated into milk recording to measure different components such as the percentage of fat and protein in milk. Unfortunately, it is used only for major components because adapted calibration equations are not available. A recent study by Soyeurt et al. (2006) provided the first calibration equations for estimating the contents of different fatty acids in milk.
Similar to its current use in routine milk recording, MIR spectrometry is an interesting alternative for providing indications of the fatty acid content in cows milk for the dairy sector (Soyeurt et al., 2006). The predicted concentrations of fatty acids based on MIR spectrometry could be used as indicator traits for the actual underlying fatty acid concentrations. The objective of this project was to study predicted results for the fatty acid contents in dairy cattle milk and for
9-desaturase activity based on MIR spectrometry of milk and milk fat.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Samples were collected from April to June 2005 from 6 breeds: Dual Purpose Belgian Blue (DPB), Holstein-Friesian (HOL), Jersey (JER), Montbeliarde (MON), Normande (NOR), and non-Holstein Meuse-Rhine-Yssel type Red and White breeds (RED). This time period covered no grazing feeding for the samples taken in April, partial grazing feeding for the samples taken in May, and total grazing feeding for the samples taken in June. Because of technical issues, the number of test days was not constant for all herds. Also, some cows were dried off or calved during this experiment. Breed composition was determined out of the known pedigrees of the animals. A certain percentage of genes were of unknown origin. Therefore, in the analysis, these genes were treated as though provided by another breed. Only 2 NOR cows were included in the animal population (Table 1
). Also, because of technical problems with spectra file acquisition, only limited fatty acid data were obtained from NOR, and we decided that these data were insufficient to study the variation in fatty acid content for this breed. The presence of the HOL genes (20 to 83.33% of average breed composition of a given herd) on all farms allowed us to study them at the same time and theoretically to separate environmental, breed, and animal effects (Table 1
).
|
Differences in Fatty Acid Profiles
Predicted contents for the different fatty acids in milk or milk fat were analyzed using the following single-trait mixed model:
![]() |
where y is the vector of observations (fat, fatty acid content in milk or milk fat); ß is the vector of fixed effects (herd x test day, parity number, stage of lactation, regressions on breed composition because many animals were crossbred; when the fatty acid content in milk fat was analyzed, a regression for the milk fat percentage was added to prevent bias attributable to variation in the fat percentage in milk); u is the vector of random repetition on animal effect; X and Z are incidence matrices; and e is the vector of random residual effects.
Animals were considered unrelated because tests using the relationship matrix and separated genetic and permanent environmental effects did not converge. Computations were done with PROC MIXED (SAS Institute, 1999), and variance components were estimated by REML. Results for regression on breed composition were reported in comparison with the reference HOL breed. To allow an optimal comparison among studied breeds and traits, results are reported in standardized values. The standard deviations used were based on the total variances given in Tables 2
, 3
, 4
, and 5
.
|
|
|
|
As mentioned,
9-desaturase activity is important for explaining the variation in fatty acid composition because it has been implicated in the production of the majority of MONO and all CLA present in milk. With the MIR predictions for 14:1 cis-9, 16:1 cis-9, 18:1, 14:0, 16:0, and 18:0, the product-to-substrate ratios were computed and analyzed with the same mixed model as that envisaged for the study of fatty acids in milk.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
9-desaturase activity among breeds, respectively. In these 3 tables, differences are presented with their associated P-values, which allowed us to test the overall significance level of breed differences and the significance of comparisons between a given breed and HOL. Because the differences may not be easy to visualize, Figures 1
|
|
|
9-desaturase activity is theoretically the ratio of 14:1 cis-9 to 14:0, because 14:1 cis-9 is derived only through desaturation (Lock and Garnsworthy, 2003), the best predictor of
9 activity for this study was the ratio of 16:1 cis-9 to 16:0. This is because of the low value for the cross-validation coefficient of determination obtained from the calibration equation that predicted the 14:1 cis-9 concentration in milk. The ratio of 16:1 cis-9 to 16:0 would therefore be preferred (Soyeurt et al., 2006). As explained before, because
9-desaturase has been implicated in the production of fatty acids, its activity was studied, because it might partially explain the observed differences in the fatty acid profiles of milk fat.
Globally, the effect of breed was often significant (P < 0.05). Tables 3
, 4
, and 5
show generally low P-values for every trait studied except for 16:1 cis-9 (P = 0.074), 18:1 (P = 0.099) in milk, SAT in fat (P = 0.334), and the ratios of 18:1 to 18:0 (P = 0.092) and 14:1 cis-9 to 14:0 (P = 0.319). In particular, Table 5
shows that breed effect had a significant impact on
9-desaturase activity for the 16:1 cis-9 to 16:0 ratio (P = 0.044), the best predictor used in our study. Kelsey et al. (2003) also showed that
9-desaturase activity differed between breeds (HOL and Brown Swiss), but this effect was minor. This minor impact might be explained by the low number of breeds studied by Kelsey et al. (2003).
As shown in Table 3
, JER showed the greatest difference in milk fatty acid content compared with HOL. The differences between JER and HOL in the fatty acid profiles of milk were generally significant except for 16:1 cis-9 (P = 0.230). In addition, the JER milk contained a SAT fraction of higher nutritional quality than that in HOL milk (Figure 1
), a result which has previously been mentioned by Palmquist and Beaulieu (1992), DePeters et al. (1995), Hermansen and Lund (1990), and Beaulieu and Palmquist (1995). Although the milk of JER cows contained a higher percentage of fat than that of HOL cows (Figure 1
), there were no significant differences between these 2 breeds concerning the proportion of SAT in milk fat after correcting for fat percentage differences (Table 4
; P = 0.460). This observation might be explained by the similar
9-desaturase activity observed for these 2 breeds (Table 5
and Figure 3
). The SAT fraction in JER milk fat was slightly higher, and the unsaturated fraction was therefore a bit lower, than in HOL (Figure 2
), but the difference was not significant (P = 0.460). This result was in agreement with Palmquist and Beaulieu (1992), who indicated that 18:1 was 15% lower in JER than in HOL milk fat, as well as with DePeters et al. (1995) and Beaulieu and Palmquist (1995).
With regard to DPB, Tables 3
and 4
show that the fatty acid profiles linked to milk or fat between this breed and HOL were slightly different; however, except for 16:0 (P = 0.033), these differences were not significant. Although the DPB had the lowest concentrations of milk fat, SAT, and MONO in milk, the concentrations of 18:1 and 18:2 cis-9,cis-12 observed for this breed were not the lowest across the studied breeds (Figure 1
). This confirms an unsaturated fraction in DPB milk fat that was more important than in the milk fat of the other studied breeds. This is supported by the fact that the milk fat of DPB contained the lowest proportion of SAT (Figure 2
). One has to consider these results in the light of a correction for milk fat percentage; therefore, the SAT contents were comparable even though the fat percentages were different. These observations might be explained by the activity of
9-desaturase for this breed (Figure 3
). Indeed, for all 3 ratios, DPB had higher
9-desaturase activity compared with all the other breeds. Even though these results were not significant, they may assist in understanding the findings from this study (Table 5
).
Similarly, the milk fat of RED and MON contained the highest proportion of SAT (Figure 2
). In particular, the MON milk fat contained a significantly higher proportion of 12:0 and 14:0 than the milk fat from HOL cows (Table 4
; respectively, P = 0.02 and P = 0.05). This might be explained by a significantly lower
9 activity for these 2 breeds compared with HOL (Figure 3
). However, Lawless et al. (1999) indicated that the MON milk fat had a low fraction of SAT and that the concentrations of 12:0 and 14:0 were equal to those observed in HOL milk fat. This divergence might be explained by the fact that the results obtained by Lawless et al. (1999) were raw phenotypic values and therefore might have been much more influenced by other factors, for example, by feeding practices and stage of lactation. In the current study, breeds were compared on the same farms, and variation among individual animals was corrected for different fixed effects.
Differences in Fatty Acid Profiles of Milk and Milk Fat Within Breeds
For all the studied components (proportions of SAT, 12:0, and so on in milk and in milk fat), the estimated values for repeatability were high (Table 2
) and ranged from to 41.33 to 60.91%. This could suggest a reasonably moderate heritability for each studied component, similar to what is known for fat (Interbull, 2005). The first study, reported by Edwards et al. (1973) for estimating the heritability of each fatty acid, indicated extremely high heritabilities, between 0.64 for 12:0 and 0.98 for 12:1 and 14:1. These values can be considered as high overestimations because a biased hypothesis was used to compute the heritability. The environmental variance was estimated from the variance component within monozygous pairs. The variance components within dizygous pairs represented the environmental variance and half the genetic variance. In spite of these overestimated values, this study was the first to show high heritability for each fatty acid in milk fat. The repeatability of milk fat components was previously shown by Karijord et al. (1982). For the same contents, the values obtained by Karijord et al. (1982) ranged from 0.10 to 0.33. The differences in values found between these 2 studies might be explained by the data analyzed and the models used. The fixed effects used by Karijord et al. (1982) were only the effect of stage of lactation and the effect of month. The heritabilities estimated by Bobe et al. (1994) were 0.18 for 18:0 and 0.39 for 14:0. A similar conclusion was drawn by Renner and Kosmack (1974a), with 0.26 for short-chain, 0.20 for medium-chain, and 0.30 for long-chain fatty acids. Because the heritability is also a part of the repeatability, our results suggest that the heritability of each fatty acid could be reasonably high.
For
9-desaturase, the estimated repeatability for the ratios of 14:1 cis-9 to 14:0, 16:1 cis-9 to 16:0, and 18:1 to 18:0 were, respectively, 43.12, 41.78, and 45.79%. These values might also suggest a reasonably moderate heritability. Lock and Garnsworthy (2003) previously reported the variation in
9 activity within breeds. Similarly, Royal and Garnsworthy (2005) reported heritability values of 0.30 and 0.19 for the ratios of 14:1 to (14:0 + 14:1) and 18:1 to (18:0 + 18:1), respectively. Lock et al. (2005) and Peterson et al. (2002) confirmed that nutritional and individual variation in
9-desaturase activity were the major factors contributing to the concentration of rumenic acid in milk.
Individual Variation
For the concentrations of fatty acids in milk, the mean scale of variation attributable to individual effects ranged from 1.30 to +2.80 standard deviation units, and the mean scale of variation attributable to breed effects ranged from 0.40 to +1.20 standard deviation units (Table 2
). The range of individual values was higher than the variation across selected breeds. The same conclusion applied to the concentrations of fatty acids in milk fat (Table 2
).
Correlations
Table 6
shows the phenotypic correlations between the studied components in milk. The values shown are Pearson correlation coefficients estimated among the individual effects. In general, high positive correlations were observed, with several exceptions. The correlations between 16:1 cis-9 and 12:0 and between 16:1 cis-9 and 14:0 were low and negative. Low correlations between MONO and 12:0, MONO and 14:0, 14:0 and 18:1, and 14:0 and 18:2 cis-9,cis-12 were also observed. These findings are important because they suggest that animal selection based on the fatty acid profile might be possible. Bobe et al. (1994) obtained a phenotypic correlation of 0.69 between 14:0 and 18:0. This value is in agreement with the value obtained in this study (0.80). Also, Karijord et al. (1982) observed lower correlations between milk components and fat than the values obtained in this study. This might be explained by differences in the data and the models used. In addition, Renner and Kosmack (1974b) obtained high genetic correlations, with fat yield equal to 1.00 for short-chain, 0.48 for medium-chain, and 0.98 for long-chain POLY. The same conclusion applies to the concentrations of fatty acids in milk fat. Indeed, high phenotypic correlations (results not shown) were found in the current study for the concentrations of fatty acids in milk fat.
|
| CONCLUSIONS |
|---|
|
|
|---|
Although the aforementioned results are interesting, they also suggest the need for further research, especially into the estimation of genetic parameters for milk fat components. An objective of future research should be to confirm the possibility of a differentiated quantitative selection of fatty acids to modify the nutritional quality of milk and milk fat. The final objective should be to produce healthier milk products based on genetic selection for modified global fatty acid profiles.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication December 23, 2005. Accepted for publication July 4, 2006.
| REFERENCES |
|---|
|
|
|---|
9-desaturase activity in dairy cows. Livest. Prod. Sci. 79:4759.
9-desaturase enzyme activity in dairy cows. Page 52 in Proc. Br. Soc. Anim. Sci., York, UK. Br. Soc. Anim. Sci., Penicuik, UK.This article has been cited by other articles:
![]() |
A. Schennink, J. M. L. Heck, H. Bovenhuis, M. H. P. W. Visker, H. J. F. van Valenberg, and J. A. M. van Arendonk Milk Fatty Acid Unsaturation: Genetic Parameters and Effects of Stearoyl-CoA Desaturase (SCD1) and Acyl CoA: Diacylglycerol Acyltransferase 1 (DGAT1) J Dairy Sci, May 1, 2008; 91(5): 2135 - 2143. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Bobe, J. A. Minick Bormann, G. L. Lindberg, A. E. Freeman, and D. C. Beitz Short Communication: Estimates of Genetic Variation of Milk Fatty Acids in US Holstein Cows J Dairy Sci, March 1, 2008; 91(3): 1209 - 1213. [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] |
||||
![]() |
P. J. Moate, W. Chalupa, R. C. Boston, and I. J. Lean Milk Fatty Acids. I. Variation in the Concentration of Individual Fatty Acids in Bovine Milk J Dairy Sci, October 1, 2007; 90(10): 4730 - 4739. [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] |
||||
![]() |
M. Mele, G. Conte, B. Castiglioni, S. Chessa, N. P. P. Macciotta, A. Serra, A. Buccioni, G. Pagnacco, and P. Secchiari Stearoyl-Coenzyme A Desaturase Gene Polymorphism and Milk Fatty Acid Composition in Italian Holsteins J Dairy Sci, September 1, 2007; 90(9): 4458 - 4465. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Bobe, G. L. Lindberg, A. E. Freeman, and D. C. Beitz Short Communication: Composition of Milk Protein and Milk Fatty Acids Is Stable for Cows Differing in Genetic Merit for Milk Production J Dairy Sci, August 1, 2007; 90(8): 3955 - 3960. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |