|
|
||||||||


* Nutritional Physiology, Department of Animal Science, Iowa State University, Ames 50011
Department of Animal Science and Industry, Kansas State University, Manhattan 66506
Animal Breeding and Genetics, Department of Animal Science, Iowa State University, Ames 50011
3 Corresponding author: dcbeitz{at}iastate.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: dairy cow genetic variation heritability milk fatty acid
Typical milk fat of cows comprises 70% saturated, 25% monounsaturated, and 5% polyunsaturated fatty acids (Grummer, 1991). A concern to consumers of dairy products is the high ratio of saturated to unsaturated fatty acids because of the linkage between intake of saturated fatty acids and various biological markers for cardiovascular disease risk in humans, including elevated blood pressure, insulin resistance, and hyper-lipidemia, particularly of low-density lipoprotein cholesterol (Vessby et al., 2001; Sacks and Katan, 2002; Mensink et al., 2003; Rasmussen et al., 2006). Cow nutrition has been used primarily to alter milk fatty acid composition (Palmquist et al., 1993; Henning et al., 2006). Feeding cows rumen-protected unsaturated lipids decreases the proportion of saturated fatty acids, specifically myristic (C14:0) and palmitic acid (C16:0), in dairy products, which results in decreased serum total and low-density lipoprotein cholesterol in humans consuming the modified dairy products (Noakes et al., 1996; Poppitt et al., 2002).
Genetic selection has been proposed to alter milk fatty acid composition (Palmquist et al., 1993; Soyeurt et al., 2006). Differences in milk fatty acid profiles within and across cattle breeds have been reported (DePeters et al., 1995; Kelsey et al., 2003; Auldist et al., 2004; Soyeurt et al., 2006). Less is known about the genetic variability of milk fatty acids within breeds. Studies that estimated sources of genetic variation of milk fatty acids were conducted in Europe several decades ago, and none of the studies used an animal model (Edwards et al., 1973; Renner and Kosmack, 1974; Karijord et al., 1982; Syrstad et al., 1982). Thus, the objective of this study was to estimate sources of genetic variation of milk fatty acids in US Holstein cows by using a single-trait, linear animal model. We hypothesized that fatty acids that are synthesized de novo in the mammary gland (C6:0 to C14:0 and partly C16:0; Palmquist et al., 1993) have significant genetic variation and, therefore, have the potential to be altered by genetic selection.
Monthly (from August 1993 to July 1994), milk samples (each sample consisted of milk taken from one morning and one evening milking in equal amounts) were collected from each lactating Holstein cow (24 to 144 mo of age, 1 to 10 lactations) at the Iowa State University Dairy Breeding Research Facility (Ankeny, IA). For fatty acid analysis, a representative subset of milk samples was selected (Bobe et al., 2007). The sampled herd comprised 2 genetic lines (high and average milk PTA), and sires are genetically representative of the US population (Bobe et al., 2007). Details about cow management and chemical analyses of milk samples have been described previously (Bobe et al., 2007).
Raw means for yields, concentrations, and relative proportions of milk fatty acids were obtained using PROC MEANS of SAS (version 9.1.3; SAS Institute, 2002). The phenotypic and genetic variation for yields, concentrations, and relative proportions of milk fatty acids were estimated by using a single-trait, mixed linear animal model:
![]() |
where Yijklmno = vector of 592 daily milk records (e.g., milk fatty acid yield expressed in g/d, milk fatty acid concentration expressed in g/L, and relative proportion of a milk fatty acid expressed as weight % of an individual milk fatty acid to total milk fatty acids). The dependent variables were distributed normally. The independent fixed effects were: 3 parity groups (i = 1, 2 and 3, and >3); SC was 7 season of calving groups (j = winter 1992, spring 1993, ..., summer 1994); LS was lactation stage (k = 1, 2, ..., 10, and >10 mo postpartum); SM was year-month of sampling (l = August 1993, September 1993, ..., July 1994); 3 diet groups (m = high, medium, and low lactational performance level); and SCCG was SCC group (n = <250,000, 250,000–1,000,000, and >1,000,000 cells/mL SCC). The matrices Z1 and Z2 were incidence matrices for the random additive genetic (a) and permanent environmental (pe) effects of the cow o, respectively. The random residual effect for each milk record is eijklmno. In addition to animals with records, 53 sires, 132 dams, 92 grandsires, and 206 granddams without production records were included in the relationship matrix as well. All fixed effects and their classifications contributed significantly to the model.
The additive genetic, permanent environmental, and residual variances for yields, concentrations, and relative proportions of milk fatty acids were estimated iteratively via derivative-free REML in MTDFREML (Boldman et al., 1995). The convergence criterion for the variance component estimates was 10–9. Heritabilities were estimated from the ratio of
a2 to the sum of
a2,
pe2, and
e2. Repeatabilities were estimated from the ratio of the sum of
a2 and
pe2 to the sum of
a2,
pe2, and
e2. Standard errors were estimated by the delta method (Dodenhoff et al., 1998). Variances were assumed to be the same for the various stages of lactation.
The estimates of phenotypic variation for yields, concentrations, and relative proportions of individual milk fatty acids and groups of milk fatty acids were large (Table 1
), which is consistent with literature values (Kelsey et al., 2003; Auldist et al., 2004; Soyeurt et al., 2006). The relative proportions of saturated fatty acids ranged from 47.15 to 77.90%, the relative proportions of monounsaturated fatty acids ranged from 18.67 to 48.56%, and the relative proportion of polyunsaturated fatty acids ranged from 1.73 to 4.91%. According to the Wisconsin Milk Marketing Board Milk Fat Roundtable (1988), the ideal nutritional milk fat would contain up to 8% saturated fatty acids and less than 10% polyunsaturated fatty acids, with the remaining fatty acids being monounsaturated fatty acids (ODonnell, 1989). Although the difference between the least saturated milk fat in this study and the ideal milk fat is large, the results indicate that some cows produce milk fat with a more desirable (nutritionally) fatty acid profile.
|
|
For relative proportions of fatty acids, the heritability estimates for individual fatty acids were small to moderate (Table 2
). Low estimates were obtained for C4:0 (h2 = 0.00), C6:0 (h2 = 0.00), C14:0 (h2 = 0.00), C18:1 (h2 = 0.06), and C18:2 (h2 = 0.00), and moderate estimates were obtained for C12:0 (h2 = 0.18), C16:1 (h2 = 0.49), and C18:0 (h2 = 0.24). Edwards et al. (1973) obtained heritability estimates between 0.64 (C12:0) and 0.95 (C18:0) in an Ayrshire twin cow study in the United Kingdom, but those results were subject to overestimation because of the study design. Our estimates are similar in magnitude to those obtained by Renner and Kosmack (1974) and Karijord et al. (1982). Therefore, we conclude that there is additive genetic variability for milk fatty acid composition between US Holstein cows.
The repeatability estimates for the relative proportions of most individual fatty acids are more than twice as high as the heritability estimates and range from 0.07 (C4:0) to 0.49 (C16:1) with most being between 0.35 and 0.49, which is close to those obtained by Soyeurt et al. (2006). Again, Karijord et al. (1982) reported repeatability estimates slightly lower than ours. The coefficients of variation (CV) of individual milk fatty acids ranged from 9.5% (C16:0) and 29.8% (C16:1) with C16:0 being the only fatty acid with a CV below 15% (Table 1
). Other milk fatty acids with high CV were C10:0 (CV = 24.5%) and C12:0 (CV = 25.4%). Given the considerable phenotypic and genetic variation in milk fatty acid composition, our results suggest that it might be feasible to alter the milk fatty acid profile by genetic selection.
There are several potential limitations to our genetic estimates. Multiparous cows were included in the analysis, which could introduce bias because culling based on milk yield might affect milk fatty acid composition. However, we reported, for the same data set, that selection for milk yield had little effect on milk fatty acid composition (Bobe et al., 2007). One has to take into consideration that this study was based on a limited number of cows with a relationship matrix based on 2 generations, and, thus, genetic estimates have large standard errors (Tables 2
), which limits reaching definite conclusions about the genetic potential to alter milk fatty acid composition in the US Holstein population. Genetic estimates for milk fatty acids in US Holstein cows, however, are not available.
In conclusion, this study demonstrated differences in yields, concentrations, and relative proportions of individual milk fatty acids between US Holstein cows. The genetic parameter estimates partially infer genetic variability in yields, concentrations, and relative proportions of individual milk fatty acids, in particular in medium- and long-chain fatty acids (C8:0 to C18:0). Genetic variability in yields and relative proportions of individual milk fatty acids would provide an opportunity to alter milk fatty acid composition by selective breeding for individual or groups of fatty acids.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
2 Present address: Cancer Prevention Fellowship Program, Office of Preventive Oncology, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892; e-mail: bobeg{at}mail.nih.gov. ![]()
Received for publication April 2, 2007. Accepted for publication November 11, 2007.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
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] |
||||
![]() |
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] |
||||
![]() |
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] |
||||
![]() |
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 | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |