Meta-analysis of factors affecting milk component yields in dairy cattle
© Lee et al.; licensee BioMed Central Ltd. 2014
Received: 25 April 2014
Accepted: 8 May 2014
Published: 5 June 2014
The objectives of this study were thus to identify most significant factors that determine milk component yield (MCY) using a meta-analysis and, if possible, to develop equations to predict MCY using variables that can be easily measured in the field. A literature database was constructed based on the research articles published in the Journal of Dairy Science from Oct., 2007 till May, 2010. The database consisted of a total of 442 observed means for MCY from 118 studies. The candidate factors that determine MCY were those which can be routinely measured in the field (e.g. DMI, BW, dietary forage content, chemical composition of diets). Using a simple linear regression, the best equations for predicting milk fat yield(MFY) and milk protein yield (MPY) were MFY = 0.351 (±0.068) + 0.038 (±0.003) DMI (R2 = 0.27), and MPY = 0.552 (±0.071) + 0.031 (±0.002) DMI - 0.004 (±0.001) FpDM (%, forage as a percentage of dietary DM) (R2 = 0.38), respectively. The best equation for predicting milk fat content (%) explained only 12% of variations in milk fat content, and none of a single variable can explain more than 5% of variations in milk protein content. We concluded that among the tested variables, DMI was the only significant factor that affects MFY and both DMI and FpDM significantly affect MPY. However, predictability of linear equations was relatively low. Further studies are needed to identify other variables that can predict milk component yield more accurately.
The manipulation of milk composition has been of interest to improve the nutritional value of milk and to increase the efficiency in manufacturing and processing of raw milk for dairy products (Jenkins and McGuire,). In addition, the decision of milk price depends on the amounts of milk component yield (MCY) in most countries including Korea, performing of quantitative analysis of factors that influence on MCY is crucially important in dairy industry.
Many studies have tried to identify such factors affecting milk yield (MY) and MCY for past decades. DMI, dietary energy and protein or individual carbohydrate (CHO) and protein fractions might be important factors in controlling MY and milk protein yield (MPY) in dairy cows ((Sutton,, Grummer,, DePeters and Cant,, Palmquist et al.,, Hristov et al.,[6, 7], Jenkins and McGuire,, Huhtanen and Hristov,). Smoler et al. suggested individual CHO fractions might be better predictors of MPY than total CHO. Hristov et al. showed a moderate linear relationship (R2 = 0.47) between DMI and MY using data set published in Journal of Dairy Science (Volumes 1 through 82). According to NRC, dietary CP was not correlated (P > 0.25) with milk protein percent, but was correlated weakly (r = 0.14; P < 0.01) with MPY. Huhtanen and Hristov reported that metabolizable protein (MP) intake was better predictor of MPY compared with CP intake.
In addition, there have been attempts to develop equations to predict MY and MPY of dairy cows (NRC,, Hristov et al.,) through a meta-analysis approach. They, however, included RDP and RUP in the equation, which cannot be easily measured in the field. NRC also presented equations to predict MY and MPY with DMI and CP contents in diets; however, the predictability was insufficient and study effect was not accounted for in developing equations. The objectives of this study were thus to identify most significant factors that determine MCY using meta-analysis based on recent studies conducted from last decade and, if possible, to derive equations to predict MCY with variables that can be easily measured in the field.
Results and discussions
Descriptive statistics for the database used for developing equations in this study
Forage DMI, kg
Milk yield and composition
Milk yield, kg/d
3.5% FCM, kg
4.0% FCM, kg
Fat yield, kg
Protein yield, kg
The best equation for predicting milk fat content (%) explained only 12% of the variations in milk fat content, and none of a single variable could explain more than 5% of the variations in milk protein content. According to the review by Jenkins and McGuire, the most sensitive component of milk to dietary manipulation was fat content, which could be changed over a range of 3 percentage units. Milk protein was more responsive to diet (over a 0.5 percentage unit range) than lactose, but less responsive than fat. Sutton reported that milk fat concentration was affected by the amount of roughage, the forage-to-concentrate ratio, the carbohydrate composition of concentrate mix, lipids, intake, and meal frequency. A reduction in the dietary forage-to-concentrate ratio usually decreases milk fat content although the degree of response varies (Sutton,). Milk fat content was fairly stable until the proportion of forage in the diet on a DM basis falls to about 50%, but with further reductions in the proportion of forage, a decrease in milk fat content occurs (Thomas and Martin,). Smith et al. indicated that the response in milk fat content to dietary supplementation of lipid was highly variable by the amount, physical form, and fatty acid composition of lipid. Sporndly observed no significant correlation between protein content of milk and protein concentration of the diet (r = 0.06), while milk protein yield and dietary protein level were correlated (r = 0.37). Jenkins and McGuire reported that reducing the proportion of forage in the diet increased both protein content and yield. Milk protein content could increase by 0.4 percentage units or more when forage proportion in the diet reduced to 10% or less of the diet DM. In addition, they indicated that low transfer efficiency (25 to 30%) of dietary protein to milk was a major factor accounting for the inability of diet to markedly alter milk protein content. Therefore, large variations or low responses to dietary manipulation resulted in low predictability for milk fat and protein content with dietary factors.
Estimates, standard errors and significance for the random components of the final candidate models for milk component yield (kg/d; n = 442)
P > |z|
The relationship of MCY with diet composition variables and intake of specific nutrients in dairy cows was also reported by previous studies (Sutton,, Grummer,, DePeters and Cant,, Palmquist et al.,, Hristov et al.,[6, 7], Jenkins and McGuire,, Huhtanen and Hristov,).The supplementation of lipids at up to 6 to 8% in the diet DM generally increases milk yields but the response in milk fat content varies widely, so that milk fat yield usually remains unchanged or increases (Macleod and Wood,, Van der Honing et al.,). Rook et al. reported that DMI, dietary NDF concentration, and digestibility of dietary OM were important in predicting MPY. Sporndly observed crude protein intake was correlated positively with milk protein content (r =0.25) and protein yield (r = 0.69). In addition to dietary CP, estimated MP (NRC,), intakes of RDP or RUP were significant predictors for MPY in dairy cows (Smoler et al.,, NRC,, Hristov et al.,). Hristov et al. showed the negative correlation of MPY with fermentable NDF and protein fraction B1 intakes. Diets for lactating dairy cows were formulated to be highly digestible and, in most situations, DMI was strongly related to intake of total digestible nutrients or energy (net energy for lactation or metabolizable energy) intake (Hristov et al.,). In Huhtanen and Hristov the best prediction models for MPY were based on total digestible nutrients (TDN), CP intake and CP degradability. Hristov et al. also suggested MY and MPY in dairy cows can be better predicted based on intake of individual nutrients than DMI alone.
Because genetic traits as well as production performance of dairy cows have been improved, a new development of equations for predicting MCY on the basis of animal parameters and/or feed characteristics was needed. Therefore, the recent experimental data obtained from 10 years before was used to re-evaluate major factors affecting MCY through meta-analysis in this study. We concluded that MFY and MPY can be predicted by DMI and FpDM which can easily be measured in field. However, predictability was relatively low. There was no variable or a combination of variables that were routinely measured and can be used to predict milk composition. More researches are needed to identify other variables that can predict milk component yield and milk composition more accurately.
The database consisted of a total of 442 treatment means from 118 studies published in the Journal of Dairy Science from Oct., 2007 till May, 2010 (volumes 90 through 93). The detailed descriptive statistics of the database was described in Table 1. The average (±SD) body weight (BW) of the cows involved in this study was 638.2 ± 53.41 kg with a minimum of 443 kg and a maximum of 778 kg. The average DMI were 22.53 ± 3.51 kg/d with a minimum of 7.0 and a maximum of 32.8 kg/d. In nutrient composition of diets, average DM, CP, EE, NDF, RDP and RUP were 54.9 ± 9.58%, 17.4 ± 1.89%, 4.4 ± 1.30%, 33.2 ± 4.81%, 10.6 ± 1.17% and 7.4 ± 5.20%, respectively. Average MY was 34.53 ± 6.76 kg/d with 3.57 ± 0.51% of average milk fat content and 3.14 ± 0.37% of average milk protein content. Average MFY and MPY were 1.22 ± 0.26 and 1.06 ± 0.22 kg, respectively.
Model development and statistical analyses
The independent variables used for explaining the variations in MCY (i.e. MPY, MFY, milk protein content in milk, and milk fat content in milk) were DMI (kg/d), BW (kg), DMIpBW (%, dry matter intake as a percentage of BW), FpDM (%, forage as a percentage of dietary DM), NDF (% DM), CP (% DM), starch (% DM), ADF (% DM), CP intake (kg/d), NDF intake (kg/d), starch intake (kg/d), and forage intake (kg/d). Among these variables, four to five predictive variables that sufficiently explained the variations in each dependent variable were selected using step-wise regression.
where y is the vector of observed MCY of which size is N; X the n × p matrix of xi,j; β the fixed effects parameter vector of which size is p; Z the designed N × (s × p) matrix that was blocked diagonally corresponding to each study (ni × p) to account for the random effect of each study; e the unknown vector of independent, identically and normally distributed random errors with mean 0 and variance σ2; N the total number of observations; s the number of studies; ni the number of observations that ith study; p the number of parameters, which is one (intercept) plus the number of variables used in the equation; and u and e are the normally distributed.
Among the acceptable regression models that had a linear combination of significant fixed effect variables, a model that had the lowest value of -2 restricted log likelihood, Akaike information criterion (AIC), the corrected Akaike information criterion (AICc) and Schwarz’s Bayesian criterion (SBC) was selected. The lowest value of those criteria above indicates a better model considering the number of observations, the number of parameters, and the maximum likelihood estimates. In the second phase, the parameters of the variables in the best model for each dependent variable, identified in the first phase, were estimated by fitting the prediction equations to a multiple regression model using GLM procedure of SAS (SAS Institute, Inc., Cary, U.S.A).
This study was financially supported by research fund of Chungnam National University.
- Jenkins TC, McGuire MA: Major advances in nutrition: Impact on milk composition. J Dairy Sci. 2006, 89: 1302-1310. 10.3168/jds.S0022-0302(06)72198-1.View ArticlePubMedGoogle Scholar
- Sutton JD: Altering milk-composition by feeding. J Dairy Sci. 1989, 72: 2801-2814. 10.3168/jds.S0022-0302(89)79426-1.View ArticleGoogle Scholar
- Grummer RR: Effect of feed on the composition of milk-fat. J Dairy Sci. 1991, 74: 3244-3257. 10.3168/jds.S0022-0302(91)78510-X.View ArticlePubMedGoogle Scholar
- DePeters EJ, Cant JP: Nutritional factors influencing the nitrogen composition of bovine milk: A review. J Dairy Sci. 1992, 75: 2043-2070. 10.3168/jds.S0022-0302(92)77964-8.View ArticlePubMedGoogle Scholar
- Palmquist DL, Beaulieu AD, Barbano DM: Feed and animal factors influencing milk-fat composition. J Dairy Sci. 1993, 76: 1753-1771. 10.3168/jds.S0022-0302(93)77508-6.View ArticlePubMedGoogle Scholar
- Hristov AN, Price WJ, Shafii B: A meta-analysis examining the relationship among dietary factors, dry matter intake, and milk and milk protein yield in dairy cows. J Dairy Sci. 2004, 87: 2184-2196. 10.3168/jds.S0022-0302(04)70039-9.View ArticlePubMedGoogle Scholar
- Hristov AN, Price WJ, Shafii B: A meta-analysis on the relationship between intake of nutrients and body weight with milk volume and milk protein yield in dairy cows. J Dairy Sci. 2005, 88: 2860-2869. 10.3168/jds.S0022-0302(05)72967-2.View ArticlePubMedGoogle Scholar
- Huhtanen P, Hristov AN: A meta-analysis of the effects of dietary protein concentration and degradability on milk protein yield and milk N efficiency in dairy cows. J Dairy Sci. 2009, 92: 3222-3232. 10.3168/jds.2008-1352.View ArticlePubMedGoogle Scholar
- Smoler E, Rook AJ, Sutton JD, Beever DE: Prediction of milk protein concentration from elements of the metabolizable protein system. J Dairy Sci. 1998, 81: 1619-1623. 10.3168/jds.S0022-0302(98)75728-5.View ArticlePubMedGoogle Scholar
- Hristov AN, Hristova KA, Price WJ: Relationship between dry matter intake, body weight, and milk yield in dairy cows: A summary of published data. J Dairy Sci. 2000, 83 (Suppl.1): 260-Google Scholar
- NRC: Nutrient Requirements of Dairy Cattle. 2001, Washington, DC, USA: National Academy Press, 7Google Scholar
- SAS, Institute Inc: User's guide: Statistics, Version 9th ed. 2002, Cary, NC: SAS Institute, Inc.Google Scholar
- Thomas PC, Martin PA: The influence of nutrient balance on milk yield and composition. Nutrition and Lactation in the Dairy Cow. 1988, London, England: ButterworthsGoogle Scholar
- Smith NE, Dunkley WL, Franke AA: Effects of feeding protected tallow to dairy cows in early lactation. J Dairy Sci. 1978, 61: 747-756. 10.3168/jds.S0022-0302(78)83643-1.View ArticleGoogle Scholar
- Sporndly E: Effects of diet on milk composition and yield of dairy cows with specical emphasis on milk protein content. Swed J Agric Res. 1989, 19: 99-106.Google Scholar
- Emery RS: Feeding for increased milk protein. J Dairy Sci. 1978, 61: 825-828. 10.3168/jds.S0022-0302(78)83656-X.View ArticleGoogle Scholar
- Macleod GK, Wood AS: Influence of amount and degree of saturation of dietary fat on yield and quality of milk. J Dairy Sci. 1972, 55: 439-445. 10.3168/jds.S0022-0302(72)85513-9.View ArticlePubMedGoogle Scholar
- Van der Honing Y, Weiman BJ, Steg A, Van Donselaar B: The effect of fat supplementation of concentrates on digestion and utilization of energy by productive dairy cows. Neth J Agric Sci. 1981, 29: 79-85.Google Scholar
- Rook AJ, Sutton JD, France J: Prediction of the yields of milk constituents in dairy-cows offered silage ad-libitum and concentrates at a flat rate. Anim Prod. 1992, 54: 313-322. 10.1017/S0003356100020766.View ArticleGoogle Scholar
- Hristov AN, Price WJ, Shafii B: Proceedings of Pacific Northwest Animal Nutrition Conference. An Overview of Dietary Factors Influencing Dry Matter Intake and Milk and Milk Protein Yields In Dairy Cows. 2002, 147-165.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.