JDS
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J. Dairy Sci. 2008. 91:385-394. doi:10.3168/jds.2007-0181
© 2008 American Dairy Science Association ®

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Stoop, W. M.
Right arrow Articles by Bovenhuis, H.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Stoop, W. M.
Right arrow Articles by Bovenhuis, H.

Genetic Parameters for Major Milk Fatty Acids and Milk Production Traits of Dutch Holstein-Friesians

W. M. Stoop*,1, J. A. M. van Arendonk*, J. M. L. Heck{dagger}, H. J. F. van Valenberg{dagger} and H. Bovenhuis*

* Animal Breeding and Genomics Centre, Wageningen University, PO Box 338, 6700 AH Wageningen, the Netherlands
{dagger} Dairy Science and Technology, Wageningen University, PO Box 8129, 6700 EV Wageningen, the Netherlands

1 Corresponding author: marianne.stoop{at}wur.nl


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
The objective of this study was to estimate genetic parameters for major milk fatty acids and milk production traits. One morning milk sample was collected from 1,918 Holstein-Friesian heifers located in 398 commercial herds in the Netherlands. Each sample was analyzed for total percentages of fat and protein, and for detailed fatty acid percentages (computed as fatty acid weight as a proportion of total fat weight). Intraherd heritabilities were high for C4:0 to C16:0, ranging from 0.42 for C4:0 to 0.71 for C10:0. Saturated and unsaturated C18 fatty acids had intraherd heritability estimates of approximately 0.25, except for C18:2 cis-9, trans-11, which was 0.42. Standard errors of the heritabilities were between 0.07 and 0.12. Genetic correlations were high and positive among C4:0 to C14:0, as well as among unsaturated C18, but correlations of C4:0 to C14:0 with unsaturated C18 were generally weak. The genetic correlation of C16:0 with fat percentage was positive (0.65), implying that selection for fat percentage should result in a correlated increase of C16:0, whereas unsaturated C18 fatty acids decreased with increasing fat percentage (–0.74). Milk fat composition can be changed by means of selective breeding, which offers opportunities to meet consumer demands regarding health and technological aspects.

Key Words: fatty acid • milk composition • genetic parameter • dairy cattle


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Milk fat contains many nutrients necessary for humans, including fat-soluble vitamins, energy, and bio-active lipids (German and Dillard, 2006). During the last decennia, however, milk fat consumption has become negatively associated with human health (Bitman et al., 1995; Jensen, 2002). Cow milk fat contains a relatively low proportion of unsaturated fatty acids (UFA) and a relative high proportion of LDL-cholesterol-increasing fatty acids (FA), mainly C14:0 and C16:0 (German and Dillard, 2006). Both have been associated with increased levels of cholesterol and increased risk of cardiovascular disease (Maijala, 2000; Jensen, 2002). In the past, some studies have questioned whether changing milk fat composition could result in a significant improvement for human health (Gibson, 1991; Maijala, 1995). In recent years, however, market attention has increasingly focused on improving the health aspects of (dairy) products. Possibilities of changing milk fat composition by genetically altering proportions of FA are therefore of interest. To study these possibilities, genetic parameters for FA have to be estimated.

Soyeurt et al. (2006a), in studying 600 milk samples of 275 animals from 5 breeds, found differences in FA proportions among breeds of cattle. Lawless et al. (1999) studied approximately 25 animals per breed and found a relatively high amount of C16:0 and a slightly lower amount of C18:0 in Holstein-Friesians, compared with Normande and Montbeliarde breeds.

Renner and Kosmack (1974b) studied 243 cows originating from 10 sires and found evidence for the existence of within-breed genetic variation in milk fat composition. Karijord et al. (1982) and Famula et al. (1995) also estimated within-breed genetic variation for milk fat composition. The study by Karijord et al. (1982) consisted of approximately 3,000 animals and 7,000 samples, and Famula et al. (1995) studied 523 animals and one sample per animal. Information on within-breed genetic variation in milk fat composition, however, is still limited. The aim of this study was to estimate genetic parameters for major milk FA and to study the relationship of major milk FA with milk production traits in the Dutch Holstein-Friesian population.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Animals and Milk Samples
This study was part of the Milk Genomics Initiative, which focuses on genetic aspects of milk composition. The Milk Genomics Initiative was designed to have approximately 2,000 cows descending from a number of selected bulls; 50 young bulls were aimed to have 20 daughters each, and 5 proven bulls were aimed to have 200 daughters each. The pedigree was the leading reason for the choice of farms; all the farms in the database of the NRS (Arnhem, the Netherlands) with at least one proven bull heifer and one test bull heifer were invited to participate in the study. When 2,000 animals were assigned to the study, no further farms were selected.

To study milk fat composition, data were available from 1,918 cows on 398 commercial herds in the Netherlands. Of those cows, 843 cows were sired by the 50 young bulls, 888 cows by the 5 proven bulls, and 187 cows by 46 other proven bulls, to have at least 3 cows per farm. Each cow was more than 87.5% Holstein-Friesian, and was between 63 and 263 DIM of the first lactation.

One milk sample of 500 mL per cow was collected between February and March 2005. Cows were milked twice daily, but only the morning milk was collected for the study to ensure the quality of the samples. Milk was cooled to 4°C within 3 h after sampling and transported to the laboratory the same morning. Sample bottles contained sodium azide (0.03% wt/wt) for conservation.

Analysis of Milk Samples
Milk fat (butter) was extracted from approximately 400 mL of milk, keeping the remaining 100 mL for other analyses. Fatty acid methyl esters were prepared from milk fat as described in ISO Standard 15884 (International Organization for Standardization-International Dairy Federation, 2002b). The methyl esters were analyzed by gas chromatography according to the 100% FA methyl ester method (International Organization for Standardization-International Dairy Federation, 2002a) with a 100-m polar column (Varian Fame Select CP 7420, Varian Inc., Palo Alto, CA) at the laboratory of the Netherlands Controlling Authority for Milk and Milk Products (Leusden, the Netherlands). The FA were identified and quantified by comparing the methyl ester chromatograms of the milk fat samples with the chromatograms of pure FA methyl ester standards, and were measured as the weight proportion of total fat weight. The chromatograms resulted in approximately 130 measurable FA peaks, of which approximately one-third could be identified. Of these, 16 major FA were used for the estimation of genetic parameters in the present study: the even-numbered FA C4:0 to C18:0; 5 identified C18:1 isomers; C18:2 cis-9,12; C18:3 cis-9,12,15; and conjugated linoleic acid (CLA) cis-9, trans-11. These 16 FA comprised 89% of the total fat. In addition to the individual FA, a number of FA groups were defined based on their potential effect on human health (German and Dillard, 2006; Table 1Go): a "neutral" group (C6–12) containing C6:0, C8:0, C10:0, and C12:0; a "negative" group (C14–16) containing C14:0 and C16:0; and a "positive" group (C18u) containing all unsaturated C18 that were part of the data set. Data were analyzed as weight proportion. To calculate the ratio of SFA:UFA, 11 additional monounsaturated FA and odd-chain FA were included.


View this table:
[in this window]
[in a new window]

 
Table 1. Trait definition: Groups of fatty acids
 
Percentages of fat and protein were determined from a 10-mL milk subsample by infrared spectroscopy by using a Fourier-transformed interferogram (MilkoScan FT 6000, Foss Electric, Hillerød, Denmark) at the certified laboratory of the Milk Control Station (Zutphen, the Netherlands).

In total, 1,918 records were analyzed for fat and protein percentages and fat composition. The NRS supplied the corresponding morning test-day milk yield of the samples; 1,783 NRS records matched our samples. Data for milk yield, fat yield, and protein yield was missing for 135 cows.

Statistical Analyses
Variance components and genetic parameters were estimated by using an Animal model in ASReml (Gilmour et al., 2002):


Formula[1]

where yijklmn is the dependent variable corresponding to the observation of animal n in herd m, with scode l, age at first calving j in season k, and DIM i; µ is the general mean; dimi is the covariate for DIM, modeled with a Wilmink curve (Wilmink, 1987); afcj is the covariate describing the effect of age at first calving; seasonk is the fixed effect with 3 classes for season of calving, summer (June to August 2004), autumn (September to November 2004), and winter (December 2004 to February 2005); scodel is the fixed effect accounting for possible differences between the groups of proven bull daughters and young bull daughters; herdm is the random effect defining groups of animals sampled in the same herd; An is the random additive genetic effect of animal n; and eijklmn is the random residual effect.

Relationships between individuals were accounted for, and the total pedigree consisted of 26,300 animals and was supplied by NRS (Arnhem, the Netherlands). Heritabilities were estimated by using univariate analyses. We defined two heritabilities:


Formula 2[2]

where {sigma}A2is the additive genetic variation, and {sigma}e2 is the residual variation. Heringstad et al. (2006) referred to this heritability as the intraherd heritability. The intraherd heritability is the parameter that is required to predict selection responses of alternative breeding strategies. The ratio of the additive genetic over the total phenotypic variance was defined as


Formula 3[3]

where {sigma}herd2 is the herd test-day variation.

The proportion of variance attributable to herd reflects the relative importance of herd effects such as feed, hygiene, and management, and could be estimated as


Formula 4[4]

To compare the relative importance of genetic and herd effects, the ratio {sigma}A2/{sigma}herd2 was calculated.

Genetic correlations were estimated by using bivariate analyses and model [1] as


Formula 5[5]

where {sigma}A1,A22 is the additive genetic covariance between trait 1 and trait 2, and {sigma}A12 and {sigma}A22 are the additive genetic variance of traits 1 and 2.

A principal component analysis PROC VARCLUS, in combination with PROC TREE in SAS (SAS 9.1, SAS Institute, 1999) was used to graphically visualize genetic correlations among major FA.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Mean and Coefficient of Variation
Mean, coefficient of variation, and 5 and 95% quartiles for individual FA, groups of FA, and milk production traits are shown in Table 2Go. The results showed considerable variation in milk fat composition among cows. Short-chain FA (C4:0 to C12:0) averaged 15% of the total fat, medium-chain FA (C14:0 and C16:0) averaged 44%, and C18:0 averaged 8%. The FA C18:1 cis-9, the most prominent among the UFA, averaged 18%, and the other UFA averaged 3.5%. The total FA in Table 2Go made up approximately 89% of the total fat; the remaining 11% consisted of a large number of FA present in small amounts. These FA included unsaturated medium-chain FA, odd-chain FA, branched-chain FA, and several long-chain FA.


View this table:
[in this window]
[in a new window]

 
Table 2. Number (n), mean,1 SD, CV, and 5 and 95% quartile for fatty acids, groups of fatty acids,2 and milk production traits, measured on a test-day morning milk sample of 1,918 cows in first lactation
 
The ratio of saturated FA (SFA) to UFA averaged 2.8, indicating that approximately 74% of the fat was saturated. This number slightly overestimates the true ratio of SFA:UFA, because the ratio was based on only 27 major FA, whereas trace amounts of some FA, mainly UFA, were not taken into account.

The coefficients of variation for individual FA ranged from 7% for C6:0 to 28% for C18:2. A low coefficient of variation (approximately 10%) was found for most saturated FA (C4:0 to C18:0), and a higher coefficient of variation (approximately 25%) was found for most unsaturated C18 FA. The highest coefficient of variation was found for long-chain polyunsaturated FA: 28% for C18:2, and 27% for C18:3. The coefficient of variation for the entire C18u group was only 11%.

Heritability and Proportion of Variance Attributable to Herd
Intraherd heritabilities for individual FA, groups of FA, and milk production traits are shown in Table 3Go. High intraherd heritabilities of 0.42 to 0.71 were found for C4:0 to C16:0. Intraherd heritabilities of 0.22 to 0.35 were found for both saturated and unsaturated C18 FA, but for CLA cis-9, trans-11, heritability was 0.42. Intraherd heritabilities for the groups C6–12 (0.67) and C18u (0.26) were in line with the results for individual FA. For C14–16, intraherd heritability was rather low (0.16) compared with intraherd heritabilities of 0.59 for C14:0 and 0.43 for C16:0. Standard errors of intraherd heritability estimates were between 0.07 and 0.12.


View this table:
[in this window]
[in a new window]

 
Table 3. Genetic, residual, and herd variation, heritabilities,1 proportion of variance explained by herd,2 and ratio of genetic over herd variance for individual fatty acids, groups of fatty acids,3 and milk production traits, measured on a test-day morning milk sample of 1,918 cows in first lactation
 
Results for the proportion of variance explained by herd are also shown in Table 3Go. Differences among herds most likely represent differences in feeding regimens among herds, but the effects of other management factors also cannot be excluded. Variance attributable to herd was lower for saturated FA (C4:0 to C18:0, approximately 0.20), than for unsaturated C18 FA (approximately 0.50). This difference, however, was not found when comparing different groups of FA: for C6–12, herd explained 27% of the variation, whereas for C18u, herd explained 31%.

The ratio of genetic variance to variance attributable to herd ({sigma}A2/{sigma}herd2; Table 3Go) gives the relative importance of genetic vs. herd effects. For C4:0 to C18:0, genetic effects were generally larger than herd effects, whereas for unsaturated C18, herd effects were larger than genetic effects.

Genetic Correlations Among Individual FA
Genetic correlations among individual FA are shown in Table 4Go. The C4:0 had a moderate negative correlation with most other FA. The C6:0 to C14:0 FA were positively correlated (0.34 to 0.96), with a weak correlation of 0.08 for C6:0 with C14:0. The C16:0 showed negative correlations with all studied FA except for C4:0, C6:0, and C18:0. The C18:0 also showed negative correlations with most other FA, but the correlations were weak. The unsaturated C18 FA were positively correlated (0.25 to 0.99), with a weak correlation of 0.12 for C18:1 cis-9 with C18:1 trans-11. A clustering tree to visualize the genetic correlations among FA is shown in Figure 1Go. Eight clusters explained more than 90% of the variance. The figure shows the clustering of C6:0 to C14:0 in one group and the clustering of unsaturated C18 in another group. This reflects the high correlations within these groups. The C4:0 and C16:0 did not cluster in one of these 2 main groups.


View this table:
[in this window]
[in a new window]

 
Table 4. Genetic correlations1 between individual fatty acids, measured on a test-day morning milk sample of 1,918 cows in first lactation
 

Figure 1
View larger version (10K):
[in this window]
[in a new window]

 
Figure 1. Cluster tree based on principal component analysis of genetic correlations among individual fatty acids. Eight clusters explained more than 90% of the variation (dotted line).

 
Genetic Correlations Among Groups of FA
Genetic correlations among the groups of FA [C4:0, C6–12, C18:0, C18u, CLA cis-9, trans-11, and the ratio of SFA:UFA] are shown in Table 5Go. The C4:0 shows moderate correlations with other FA, as does C6–12. Both C14–16 and C18:0 showed negative correlations with (unsaturated) C18 and CLA cis-9, trans-11, and a positive correlation with the ratio of SFA:UFA. The positive correlation was expected, because C14–16 and C18:0 are the main saturated fractions, so a relative increase will also increase the ratio. The C18u and CLA cis-9, trans-11 were highly correlated (0.71), and both showed a strong negative correlation with the ratio of SFA:UFA, because they are the main unsaturated fractions.


View this table:
[in this window]
[in a new window]

 
Table 5. Genetic correlations1 between 4:0, 18:0, conjugated linoleic acid (CLA) cis-9, trans-11, ratio of SFA:UFA, and groups of fatty acids,2 measured on a test-day morning milk sample of 1,918 cows in first lactation
 
Genetic Correlations of FA with Production Traits
Genetic correlations of individual FA and groups of FA with milk production traits are shown in Table 6Go. Fat percentage showed a positive correlation with C16:0 (0.65), and negative correlations with C14:0 (–0.43) and all unsaturated C18 FA (–0.43 to –0.78). Results were similar for the correlations of fat percentage with groups of FA: fat percentage showed a positive correlation with C14–16 (0.65) and a strong negative correlation with C18u (–0.74).


View this table:
[in this window]
[in a new window]

 
Table 6. Genetic correlations1 between individual fatty acids, groups of fatty acids,2 and milk production traits, measured on a test-day morning milk sample of 1,918 cows in first lactation
 
Protein percentage and fat yield also had positive correlations with C16:0 and C14–16, and negative correlations with all unsaturated C18 FA. For protein yield and milk yield, however, results were opposite: a negative correlation with C14–16 and moderate to strong positive correlations with all unsaturated C18 FA.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
The results of the current study show that there is considerable genetic variation for proportions of FA, with genetic variation being high for C4:0 to C16:0 and moderate for C18 FA. The moderate coefficient of variation in combination with a moderate to high intraherd heritability indicates that FA proportions can be changed by genetic selection. The genetic correlations suggest that the individual FA react differently to selection. Selection on fat percentage seems to have little effect on C6–12, whereas it is expected to lead to an increase in C14–16 and a decrease in C18u.

Samples
In the current study, one morning milk sample per cow was analyzed. Only morning milk was used to ensure a quick transport from the farm to the laboratory. Usually, however, milk production data are analyzed by using a mixed evening-morning sample. This difference could have affected the obtained results. A study by Morris et al. (2005), using an interval of 8 h:16 h, showed high correlations of approximately 0.6 to 0.9 for proportions of FA between the evening and morning samples.

Chromatograms of the milk fat samples showed a total of 133 different peaks. About one-third of these peaks could be identified. The current study focused on the identified major FA, constituting approximately 90% of total fat, which were even-chain FA (C4:0 to C18:0); 5 monounsaturated C18:1; C18:2 cis-9,12; C18:3 cis-9,12,15; and the CLA isomer cis-9, trans-11. Because FA were measured as weight proportion, SFA and UFA were dependent variables and were analyzed as the ratio of SFA:UFA.

Heritabilities
In the current study, the effect of herd was modeled as a random effect, because our interest was also in variation attributable to herd effects. Most studies model herd as a fixed effect. Modeling herd as fixed or random effect has, among others, consequences for the total phenotypic variance, and consequently for the heritability estimates. We also performed an analysis in which we took herd as a fixed effect. The heritabilities that were obtained from these analyses (results not shown) were very similar to what was termed the intraherd heritability (hIH 2 ). This makes the intraherd heritability reported in the current study comparable to the heritability reported by other studies that modeled herd as a fixed effect.

Renner and Kosmack (1974b) and Karijord et al. (1982) estimated heritabilities for short-chain FA between 0.13 and 0.26, for medium-chain FA between 0.06 and 0.11, and for unsaturated C18 at 0.04 (Renner and Kosmack, 1974b), which were lower than estimates found in the current study. Both studies (Renner and Kosmack, 1974b; Karijord et al., 1982) also found low heritabilities for milk yield (0.36 and 0.09) and fat percentage (0.28 and 0.09) compared with the current study. The milk production trait heritabilities found in this study were, however, in line with more recent literature (e.g., Vos and Groen, 1998; Ikonen et al., 1999; Ojala et al., 2004).

Short-Chain vs. Long-Chain FA
Intraherd heritability for C4:0 to C14:0 (approximately 0.60) was higher than for unsaturated C18 (approximately 0.25), which is in agreement with the findings of Renner and Kosmack (1974b) and Karijord et al. (1982). The proportion of variance attributable to herd is smaller for C4:0 to C14:0 (approximately 0.25) than for unsaturated C18 FA (approximately 0.50). There can be many reasons for variance attributable to herd (Jensen, 2002); feed differences among farms, but also other management factors, might play a role. Possible reasons for this could be that short-chain FA are synthesized de novo by the cow. Long-chain FA, however, originate predominantly from dietary FA, and because plant material contains mainly long-chain FA, differences in diet affect long-chain FA more than short-chain FA. The difference caused by herd effects, however, was less clear in the FA groups, where herd consistently explained approximately 30% of the variance. The difference between the FA was further evident from genetic correlations between FA. Strong positive genetic correlations were reported for both C6:0 to C14:0 and unsaturated C18 FA. To the contrary, genetic correlations of C6:0 to C14:0 with unsaturated C18 FA were weak, as were genetic correlations of C4:0, C16:0, and C18:0 with the other FA. Few papers have studied genetic correlations between FA. Genetic correlations found in the current study were comparable to the results reported by Karijord et al. (1982).

Explanation for the grouping of the FA can be found in the biological pathways of synthesis. The FA C4:0 (butyric acid) was negatively correlated with almost all other FA. It is formed partly in the rumen by bacterial processes, together with acetate and propionate, and is a precursor for most other short- and medium-chain FA. Increased de novo synthesis will possibly convert more C4, so less C4 is present in milk, hence the negative correlation. The C6:0 to C16:0 FA are synthesized de novo in a FA cycle starting from C2 and C4 (Bobe et al., 1999). The C16:0, however, is partly synthesized de novo and is partly excreted from blood, which might explain the correlations found for this FA. The unsaturated C18 FA originate mainly from dietary FA, and their proportions are highly dependent on rumen biohydrogenation and on {Delta}9-desaturase enzymatic activity in the mammary gland (MacGibbon and Taylor, 2006).

Because of the synthesis pathways, short- and medium-chain FA are expected to be under stronger genetic control than the long-chain FA. This is also reflected in the higher heritability estimates and the smaller influence of herd for short- and medium-chain FA, compared with long-chain FA. Genetic selection is therefore likely to be more effective for short- and medium-chain FA.

Correlation Between Fat Percentage and FA
Average milk fat percentage in the Netherlands has increased from 3.7% in 1950 to 4.4% in 2005 (NRS, 2006). A positive genetic correlation was found in the present study between fat percentage and C16:0, and a negative genetic correlation was found between fat percentage and unsaturated C18 FA. As a result of the increase in fat percentage, a correlated increase in C16:0 and a decrease in unsaturated C18 FA were expected. In 1974, Renner and Kosmack (1974a) found a fat percentage of 4%, with 25.5% C16:0 and 31.1% unsaturated C18 FA. In the current study, the fat percentage was 4.36%, with 32.6% C16:0 and 21.6% unsaturated C18 FA. Thus, comparing results from this study with those from Renner and Kosmack (1974a), fat percentage has increased, with a strong increase in the proportion of C16:0 and an even stronger decrease in the proportions of unsaturated C18 FA. Although there may be many reasons for these differences, such as breed, season, lactation stage, time of day, or feed, the change in milk fat composition could have been a correlated response to selection for fat, favoring C16:0 rather than increasing all FA simultaneously.

Selection for Fat Composition
In the past, some studies questioned the need to change the milk fat composition for mainly 2 reasons: 1) because large changes in milk fat composition would be required to substantially decrease risks to human health (Maijala, 1995); and 2) because changes that are positive for one product might be detrimental for other products. The latter reason might imply that multiple breeding goals are needed and that the entire production chain would have to be adapted, for example by separately collecting milk for different end uses (Gibson, 1991). In recent years, however, there has been increased awareness among consumers and increased market demands for differentiated, specific, and healthy dairy products, even if these require large adaptations within the current production structure. Even if individual products have small effects because they are only a small part of the total diet, ultimately the sum of all improved products might result in an overall beneficial health effect, for example, by decreasing the intake of C14:0 and C16:0 FA and increasing the intake of unsaturated C18 FA. However, before a breeding program can focus on milk fat composition, several issues need to be resolved. Among others, measuring FA proportions by gas chromatography analysis is too expensive to use in routine analysis, and cheaper methods have to become available, for example, using infrared spectroscopy (Soyeurt et al., 2006b). Farmers in many countries are currently paid based on yields of fat and protein, which tend to increase rather than decrease the amounts of C14–16, so FA-based payment should possibly be introduced. Finally, the direction of selection might depend on the purpose of the milk product, because changes that are favorable for one product might be unfavorable for others. In addition to these practical considerations, it would seem worthwhile to study repeated samples from all cows to estimate parameters more accurately and to study physiological changes when selecting for or against certain FA. This study provides the necessary first data to evaluate the possibilities of improving milk fat composition by selective breeding.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Results have shown considerable genetic variation for all milk FA, with C4:0 to C16:0 having higher intraherd heritabilities (around 0.60) than C18 FA (around 0.25). High genetic correlations exist within the groups of short- and long-chain FA, which coincide with the origin of the FA and the biological pathways of synthesis. Selection for fat may lead to an increased proportion of C16:0 and a decreased proportion of unsaturated C18 FA because of the high genetic correlations. This study shows that it is possible to change the milk FA composition by genetic selection.


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
This study was part of the Milk Genomics Initiative, funded by Wageningen University, NZO (Dutch Dairy Organization, Zoetermeer, the Netherlands), breeding company HG (Arnhem, the Netherlands), and technology foundation STW (Utrecht, the Netherlands). The authors thank the owners of the herds, the Milk Control Station (Zutphen, the Netherlands), and NRS (Arnhem, the Netherlands) for their help in collecting the data.

Received for publication March 9, 2007. Accepted for publication September 7, 2007.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 


Bitman, J., D. L. Wood, R. H. Miller, J. C. Wilk, and E. D. Moore. 1995. Comparison of lipid composition of milk from half-Danish Jersey cows and United States Jersey cows. J. Dairy Sci. 78:655–658.[Abstract]

Bobe, G., D. C. Beitz, A. E. Freeman, and G. L. Lindberg. 1999. Associations among individual proteins and fatty acids in bovine milk as determined by correlations and factor analysis. J. Dairy Sci. 66:523–536.

Famula, T. R., J. F. Medrano, E. J. DePeters, and S. L. Berry. 1995. Estimation of heritability and genetic correlations among fatty acid components of milk in Holstein cows. J. Dairy Sci. 78(Suppl.1):194. (Abstr.)[Abstract]

German, J. B., and C. J. Dillard. 2006. Composition, structure and absorption of milk lipids: A source of energy, fat-soluble nutrients and bio-active molecules. Crit. Rev. Food Sci. Nutr. 46:57–92.[CrossRef][Medline]

Gibson, J. P. 1991. The potential for genetic change in milk fat composition. J. Dairy Sci. 74:3258–3266.[Abstract/Free Full Text]

Gilmour, A. R., B. J. Gogel, B. R. Cullis, S. J. Welham, and R. Thompson. 2002. ASReml User Guide Release 1.0. VSN International Ltd., Hemel Hempstead, UK.

Heringstad, B., D. Gianola, Y. M. Chang, J. Odegard, and G. Klemetsdal. 2006. Genetic associations between clinical mastitis and somatic cell score in early first-lactation cows. J. Dairy Sci. 89:2236–2244.[Abstract/Free Full Text]

Ikonen, T., K. Ahlfors, R. Kempe, M. Ojala, and O. Ruottinen. 1999. Genetic parameters for the milk coagulation properties and prevalence of noncoagulating milk in Finnish Dairy cows. J. Dairy Sci. 82:205–214.[Abstract]

International Organization for Standardization-International Dairy Federation. 2002a. Milk fat—Determination of the fatty acid composition by gas-liquid chromatography. ISO 15885-IDF 184. International Dairy Federation, Brussels, Belgium.

International Organization for Standardization-International Dairy Federation. 2002b. Milk fat—Preparation of fatty acid methyl esters. ISO 15884-IDF 182. International Dairy Federation, Brussels, Belgium.

Jensen, R. G. 2002. Invited review: The composition of bovine milk lipids: January 1995 to December 2000. J. Dairy Sci. 85:295–350.[Abstract]

Karijord, O., N. Standal, and O. Syrstad. 1982. Sources of variation in composition of milk fat. Z. Tierz. Zcht. Biol. 99:81–93.

Lawless, F., C. Stanton, P. L’Escop, R. Devery, P. Dillon, and J. J. Murphy. 1999. Influence of breed on bovine milk cis-9, trans-11-conjugated linoleic acid content. Livest. Prod. Sci. 62:43–49.[CrossRef]

MacGibbon, A. K. H., and M. W. Taylor. 2006. Composition and Structure of Bovine Milk Lipids. Pages 1–42 in Advanced Dairy Chemistry. 3rd ed. Vol. 2: Lipids. P. F. Fox and P. L. H. McSweeney, ed. Springer, New York, NY.

Maijala, K. 1995. Current possibilities of breeding to influence milk composition from the viewpoint of human nutrition and health. Page 23 in Proc. NJF/NMR—Seminar No. 252. Milk in Nutrition: Effects of Production and Processing Factors. Turku, Finland. S. Mantere-Alhonen and K. Maijala, University Printing House, Helsinki, Finland.

Maijala, K. 2000. Cow milk and human development and well-being. Livest. Prod. Sci. 65:1–18.[CrossRef]

Morris, C. A., S. M. Hickey, and M. P. Agnew. 2005. Repeatability of fatty acid composition in pm:am milk samples from cows. Proc. N. Z. Soc. Anim. Prod. 65:297–299.

NRS. 2006. Year Statistics 2005, NRS, Arnhem, the Netherlands.

Ojala, M., T. Seppanen, A. Tyriseva, and T. Ikonen. 2004. Associations of milk protein polymorphism and milk production traits and udder health traits in Finnish Ayrshire cows. EAAP 2004, GPh1.7. Eur. Assoc. Anim. Prod., Bled, Slovenia.

Renner, E., and U. Kosmack. 1974a. Genetic aspects concerning the composition of fatty acids in milk fat. 1. Delimitation to non-genetic factors. Zuchtungskunde 46:90–99.

Renner, E., and U. Kosmack. 1974b. Genetic aspects concerning fatty acid composition of milk fat. 2. Fatty acid pattern of milk from progeny groups. Zuchtungskunde 46:217–226.

SAS Institute. 1999. SAS OnlineDoc, Version 8. SAS Inst. Inc., Cary, NC.

Soyeurt, H., P. Dardenne, F. Dehareng, G. Lognay, D. Veselko, M. Marlier, C. Bertozzi, P. Mayeres, and N. Gengler. 2006b. Estimating fatty acid content in cow milk using mid-infrared spectrometry. J. Dairy Sci. 89:3690–3695.[Abstract/Free Full Text]

Soyeurt, H., P. Dardenne, A. Gillon, C. Croquet, S. Vanderick, P. Mayeres, C. Bertozzi, and N. Gengler. 2006a. Variation in fatty acid contents of milk and milk-fat within and across breeds. J. Dairy Sci. 89:4858–4865.[Abstract/Free Full Text]

Vos, H., and A. F. Groen. 1998. Altering milk protein/fat ratio: Results of a selection experiment in dairy cattle. Livest. Prod. Sci 53:49–55.[CrossRef]

Wilmink, J. B. M. 1987. Adjustment of test-day milk, fat and protein yield for age, season and days-in-milk. Livest. Prod. Sci. 16:335–348.[CrossRef]


This article has been cited by other articles:


Home page
J DAIRY SCIHome page
R. M. Demeter, G. C. B. Schopen, A. G. J. M. O. Lansink, M. P. M. Meuwissen, and J. A. M. van Arendonk
Effects of milk fat composition, DGAT1, and SCD1 on fertility traits in Dutch Holstein cattle
J Dairy Sci, November 1, 2009; 92(11): 5720 - 5729.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
K. J. E. van Hulzen, R. C. Sprong, R. van der Meer, and J. A. M. van Arendonk
Genetic and nongenetic variation in concentration of selenium, calcium, potassium, zinc, magnesium, and phosphorus in milk of Dutch Holstein-Friesian cows
J Dairy Sci, November 1, 2009; 92(11): 5754 - 5759.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
J. M. L. Heck, H. J. F. van Valenberg, J. Dijkstra, and A. C. M. van Hooijdonk
Seasonal variation in the Dutch bovine raw milk composition
J Dairy Sci, October 1, 2009; 92(10): 4745 - 4755.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
M. Glantz, H. Lindmark Mansson, H. Stalhammar, L.-O. Barstrom, M. Frojelin, A. Knutsson, C. Teluk, and M. Paulsson
Effects of animal selection on milk composition and processability
J Dairy Sci, September 1, 2009; 92(9): 4589 - 4603.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
W. M. Stoop, A. Schennink, M. H. P. W. Visker, E. Mullaart, J. A. M. van Arendonk, and H. Bovenhuis
Genome-wide scan for bovine milk-fat composition. I. Quantitative trait loci for short- and medium-chain fatty acids
J Dairy Sci, September 1, 2009; 92(9): 4664 - 4675.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
A. Schennink, W. M. Stoop, M. H. P. W. Visker, J. J. van der Poel, H. Bovenhuis, and J. A. M. van Arendonk
Short communication: Genome-wide scan for bovine milk-fat composition. II. Quantitative trait loci for long-chain fatty acids
J Dairy Sci, September 1, 2009; 92(9): 4676 - 4682.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
L. F. De La Fuente, E. Barbosa, J. A. Carriedo, C. Gonzalo, R. Arenas, J. M. Fresno, and F. San Primitivo
Factors influencing variation of fatty acid content in ovine milk
J Dairy Sci, August 1, 2009; 92(8): 3791 - 3799.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
W. M. Stoop, H. Bovenhuis, J. M. L. Heck, and J. A. M. van Arendonk
Effect of lactation stage and energy status on milk fat composition of Holstein-Friesian cows
J Dairy Sci, April 1, 2009; 92(4): 1469 - 1478.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
G. C. B. Schopen, J. M. L. Heck, H. Bovenhuis, M. H. P. W. Visker, H. J. F. van Valenberg, and J. A. M. van Arendonk
Genetic parameters for major milk proteins in Dutch Holstein-Friesians
J Dairy Sci, March 1, 2009; 92(3): 1182 - 1191.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
P. M. Kgwatalala, E. M. Ibeagha-Awemu, A. F. Mustafa, and X. Zhao
Influence of stearoyl-coenzyme A desaturase 1 genotype and stage of lactation on fatty acid composition of Canadian Jersey cows
J Dairy Sci, March 1, 2009; 92(3): 1220 - 1228.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
M. Mele, R. Dal Zotto, M. Cassandro, G. Conte, A. Serra, A. Buccioni, G. Bittante, and P. Secchiari
Genetic parameters for conjugated linoleic acid, selected milk fatty acids, and milk fatty acid unsaturation of Italian Holstein-Friesian cows
J Dairy Sci, January 1, 2009; 92(1): 392 - 400.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
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]


Home page
J DAIRY SCIHome page
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]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Stoop, W. M.
Right arrow Articles by Bovenhuis, H.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Stoop, W. M.
Right arrow Articles by Bovenhuis, H.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS