|
|
||||||||

,2
* Laboratory for Animal Nutrition and Animal Product Quality, Ghent University, Proefhoevestraat 10, 9090 Melle, Belgium
Animal Nutrition Group, Wageningen University, PO Box 338, 6700AH Wageningen, The Netherlands
Institute of Grassland and Environmental Research, Plas Gogerddan, Aberystwyth, Ceredigion SY23 3EB, UK
Animal Sciences Group, Wageningen University, Edelhertweg 15, 8200 AB Lelystad, The Netherlands
# Department of Animal Nutrition and Husbandry, Ministry of the Flemish Community-Agricultural Research Centre, Scheldeweg 68, 9090 Melle, Belgium
1 Corresponding author: veerle.fievez{at}UGent.be
| ABSTRACT |
|---|
|
|
|---|
Key Words: rumen fermentation volatile fatty acid milk odd- and branched-chain fatty acids
| INTRODUCTION |
|---|
|
|
|---|
The objective of this work was to evaluate the potential of odd- and branched-chain fatty acids in milk to predict rumen proportions of VFA. First, we investigated the relationship between concentrations of odd-and branched-chain fatty acids in milk and rumen proportions of VFA, using data from 10 feeding experiments with rumencannulated dairy cows. Second, we validated the developed equations using independent data of experiments reported in literature.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Experiment 1 (n = 24).
Experimental design and diets are as described by Dewhurst et al. (2003). Briefly, the experiment was according to a 4-period incomplete changeover design, in which 6 cows in the beginning of the lactation were used to test 6 dietary treatments. Each experimental period lasted for 3 wk, of which the first 2 wk were for adaptation. Each cow was offered 4 different diets. Cows received 8 kg/d of a dairy concentrate (Dewhurst et al., 2003), in 3 portions: 3 kg at each milking (0730 and 1600 h) and 2 kg at 1200 h. Cows had ad libitum access to 1 of 6 silages: grass, red clover, white clover, alfalfa, and 50/50 (DM basis) mixtures of grass and red clover and grass and white clover. Each forage treatment comprised a proportional mixture of all cuts taken in the year. Fresh forage was distributed daily at 0900 h. Samples of rumen fluid were taken using automated equipment at 2-h intervals during 1 d (12 samples/cow), acidified, and analyzed for VFA (Dewhurst et al., 2003). Milk samples were taken from 4 consecutive milkings and stored frozen.
Experiment 2 (n = 16).
This experiment was a 4 x 4 Latin square. Four dairy cows (90 ± 34 DIM at the beginning of the experiment) were offered diets varying in forage-to-concentrate ratio. Dietary treatments were based on ad libitum access to ryegrass silage and a dairy concentrate (Moorby et al., 2006; Vlaeminck et al., 2006b) with forage:concentrate ratios of 80:20, 65:35, 50:50, 35:65 on a DM basis (Moorby et al., 2006; Vlaeminck et al., 2006b). Each experimental period lasted for 4 wk of which the first 2 wk were for adaptation. Fresh forage was distributed daily at 0900 h whereas concentrates were distributed twice daily in equal portions at milking (0800 and 1600 h; Moorby et al., 2006; Vlaeminck et al., 2006b). Samples of rumen fluid were taken using automated equipment at 2-h intervals during 1 d (12 samples/cow), acidified, and analyzed for VFA (Dewhurst et al., 2003). Milk samples were taken from 4 consecutive milkings and stored frozen.
Experiment 3 (n = 16).
Experimental procedures were described previously (Hindle et al., 2005). Each experimental period lasted for 4 wk of which the first 2 wk were for adaptation. Four multicannulated dairy cows (80 ± 18 DIM at the beginning of the experiment) received either a control diet, consisting of grass silage (43% of DM), ensiled sugar beet pulp (11% of DM), and a concentrate mixture with 70% dried sugar beet pulp. Dried sugar beet pulp of the concentrate was replaced either by native potato starch, cornmeal, or wheat meal in each of the 3 experimental diets. Diets were distributed twice daily in equal portions as a TMR at 0600 and 1700 h. Rumen fluid was sampled at 4-h intervals during one 24-h period. After each sampling, samples were acidified and stored at 20°C until analysis of VFA (Hindle et al., 2005). Milk samples were taken from 8 consecutive milkings and stored frozen.
Experiment 4 (n = 16).
Experimental procedures were described in detail previously (Bruinenberg et al., 2004). Briefly, 4 rumen-cannulated multiparous Holstein cows (249 ± 76 DIM at the beginning of the experiment) were assigned to a 4 x 4 Latin square experiment. Each experimental period lasted for 3 wk, of which the first 2 wk were for adaptation. Forages differed between the 4 dietary treatments and consisted of different combinations of 3 grassland silages (Bruinenberg et al., 2004). At 0600 and 1600 h, cows received 40 and 60% of the daily DM offered, respectively. Over one 24-h period, samples from rumen fluid were taken at 1600, 1800, 2000, 2200, 2400, 0300, 0600, 0800, 1000, and 1300 h, and analyzed for VFA (Bruinenberg et al., 2004). Milk samples were taken from 4 consecutive milkings and stored frozen.
Experiments 5 and 6 (n = 50).
These two experiments were both 5 x 5 Latin squares. Five dairy cows in early lactation (45 ± 14 DIM at the beginning of the experiment; experiment 5) and late lactation (236 ± 14 DIM at the beginning of the experiment; experiment 6) were offered diets varying in source of forage and concentrate (55/45, DM basis). Dietary treatments were based on ad libitum access to 1 of the 5 TMR: 1) a mixture (50/50, DM basis) of ryegrass silage (CP: 16.6% of DM, sugar: 3.6% of DM, NDF: 50.9% of DM) and corn silage (CP: 7.5% of DM, sugar: 0.6% of DM, starch: 31.4% of DM, NDF: 40.0% of DM) as forage and a mixture (50/50, DM basis) of 2 concentrates either rich in structural (CP: 19.4% of DM, sugar: 10.8% of DM, starch: 1.4% of DM, NDF: 32.0% of DM) or in nonstructural carbohydrates (CP: 19.6% of DM, sugar: 10.5% of DM, starch: 30.3% of DM, NDF: 14.6% of DM); 2) rye-grass silage as forage and the concentrate mixture (50/ 50, DM basis) as in diet 1; 3) corn silage as forage and the concentrate mixture (50/50, DM basis) as in diet 1; 4) a mixture (50/50, DM basis) of ryegrass silage and corn silage as forage and a concentrate rich in structural carbohydrates; 5) a mixture (50/50, DM basis) of rye-grass silage and corn silage as forage and a concentrate rich in nonstructural carbohydrates. Each experimental period lasted for 3 wk, of which the first 2 wk were for adaptation. Samples of rumen fluid were taken frequently during one 12-h interval (at least 15 samples/ cow) based on the cows ingestion pattern, acidified, and analyzed for VFA (Chilibroste et al., 1998). Milk samples were taken from 4 consecutive milkings and stored frozen.
Experiments 7 and 8 (n = 18).
These 2 experiments were both 3 x 3 Latin squares. Three dairy cows (294 ± 148 DIM at the beginning of the experiment) were offered diets with corn and grass silage as forage with a standard dairy concentrate. Corn silage was taken from 2 different varieties varying in starch content. Dietary treatments were based on ad libitum access of forage and 5.1 kg/d of the standard dairy concentrate. Each experimental period lasted for 13 d, of which the first 10 d were for adaptation. Samples of rumen fluid were taken before the morning feeding and 1, 2, 3, 5, and 8 h after the morning feeding during the last 3 consecutive days. Before VFA analysis (De Boever et al., 2005), rumen samples from the 3 d were pooled according to their sampling time (6 samples per animal and per treatment). Milk samples were taken from 4 consecutive milkings, pooled and stored frozen.
Experiments 9 and 10 (n = 8).
These 2 experiments were both 2 x 2 Latin squares. Two dairy cows (278 ± 108 DIM at the beginning of the experiment) were offered diets with corn and grass silage as forage with a standard dairy concentrate. Corn silage was from 2 different varieties, inducing variation in rumen bypass starch content. Dietary treatments were based on ad libitum access of the forage and 5.1 kg/d of a standard dairy concentrate. Each experimental period lasted for 13 d of which the first 10 d were for adaptation. Samples of rumen fluid were taken before the morning feeding and 1, 2, 3, 5, and 8 h after the morning feeding during 3 consecutive days. Before VFA analysis (De Boever et al., 2005), rumen samples of the 3 d were pooled according to their sampling time (6 samples per animal and per treatment). Milk samples were taken from 4 consecutive milkings, pooled, and stored frozen.
Milk Fatty Acid Analysis
In all the experiments, milk samples were extracted, methylated, and analyzed by GLC as described by Vlaeminck et al. (2005). Briefly, in the first step, samples were extracted with ammonium hydroxide solution, ethanol, diethyl ether, and petroleum ether. In the second extraction step, ethanol, diethyl ether, and petroleum ether were used, and in the final extraction step, the solvents used were diethyl ether and petroleum ether. Extracts from the 3 consecutive steps were combined, evaporated, methylated, and analyzed separately for short-chain fatty acids (C4:0 to C10:0) and medium- and long-chain fatty acids (C12:0 to C24:0). Standard curves were used to determine the response factors for milk short-chain fatty acids, taking into account tridecanoic acid (Sigma, Bornem, Belgium) as the internal standard, whereas the other fatty acids were quantified with nonadecanoic acid as the internal standard (Sigma). Methylation and analysis by GLC were described previously (Vlaeminck et al., 2005).
Calculations
Mean 24-h molar proportions of VFA were calculated from the results of the individually analyzed samples, taken in relation to the cows feeding pattern as described before for the different experiments. Because time intervals between samplings were not the same in each experiment, hourly values were calculated by linear interpolation of values for the next and previous sampling. From these 24 values, mean molar proportions of acetate, propionate, and butyrate [mmol/mol (acetate + propionate + butyrate)] for each sampling day were calculated.
Milk odd- and branched-chain fatty acids were expressed as g/100 g of fatty acids and individual odd-and branched-chain fatty acids were used to construct the model. However, because cis-9 C17:1 is a desaturation product of C17:0 in the mammary gland (Fievez et al., 2003a), the sum of C17:0 and cis-9 C17:1 was used.
Independent Data Set to Evaluate Accuracy of Prediction
Qualitative Evaluation.
Because of the limited data available in literature reporting both rumen fermentation pattern and milk odd- and branched-chain fatty acids, we first performed a qualitative evaluation of the regression equations: using milk odd- and branched-chain fatty acids, reported by Kraft et al. (2003), Loor et al. (2005a, c), and Van Nespen et al. (2005), we calculated the associated VFA proportions from the regression equations based on milk odd- and branched-chain fatty acids developed in the current study. Because the VFA proportions were not measured in these studies (Kraft et al., 2003; Loor et al., 2005a,c; Van Nespen et al., 2005), the qualitative evaluation assessed the ability of the equations to predict well-known shifts in the fermentation pattern induced by varying roughage proportions (Kraft et al., 2003) or supplementation of fish (Loor et al., 2005a,c) or vegetable oils (Loor et al., 2005a) and rumen-degradable starch (Van Nespen et al., 2005).
Quantitative Evaluation.
Data from the experiments described by Shingfield et al. (2003, 2005) and Loor et al. (2005b) were used to evaluate the accuracy of equations developed to predict rumen proportions of VFA based on milk odd- and branched-chain fatty acids. To our knowledge, these are the only published articles with information on both rumen fermentation pattern and milk odd- and branched-chain fatty acids in the same experiments, albeit in separate papers. The resulting test data set consisted of 14 observations.
Statistics
All statistical analyses were performed using SPSS 12.0 (SPSS software for Windows, release 12.0, SPSS, Inc., Chicago, IL).
Linear Regression Analyses.
Regression analyses were performed using the linear mixed models procedure with the discrete effect of study included as a random variable (St-Pierre, 2001). All individual odd-and branched-chain fatty acids were introduced in the model and equations were developed by running multiple iterations and removing the least significant effect at each iteration. Variables were kept in the model when P-values were below 0.10. Regression equations were evaluated based on the standard error and coefficients of variation and by regressing residuals on the predicted values. Predicted values were centered by subtracting the mean of all predicted values from each prediction (St-Pierre, 2003). Mean biases were assessed using the intercepts of the regression equations, and the slopes of the regression equations were used to determine the presence of linear biases.
Independent Data Sets to Evaluate Accuracy of Prediction.
The regression equations developed in this study were validated with independent data sets using the mean square prediction error (MSPE; Bibby and Toutenberg, 1977):
![]() |
where n is the number of observations, and yi and
i are the observed and predicted values, respectively. The square root of MSPE is expressed in the same units as the observed values and a comparison of the root MSPE as a percentage of the observed mean provides an indication of the overall error of prediction. The MSPE was further decomposed into 3 components; mean bias, line bias (deviation from the regression slope from one), and the disturbance proportion indicating random deviation, which cannot be accounted for by bias or regression deviation. Further, residuals were regressed on the predicted values as described before.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Model Development
The regression equation showed that milk iso C14:0 and iso C15:0 are positively related to rumen proportions of acetate. In a literature review, Vlaeminck et al. (2006a) described a positive relation between dietary NDF content and milk iso C14:0 and iso C15:0. In addition, studies reporting odd- and branched-chain fatty acids of mixed rumen bacteria show an increase in these fatty acids when forage:concentrate ratio increases (Bas et al., 2003; Vlaeminck et al., 2006b). Higher dietary NDF and an increase in forage:concentrate ratio generally increases the importance of cellulolytic bacteria (Dehority and Orpin, 1988; Weimer et al., 1999). Fatty acid analysis of pure strains of rumen bacteria have shown that some major cellulolytic bacteria are enriched in iso C14:0 and iso C15:0 (see, for example, Ifkovitz and Ragheb 1968; Minato et al., 1988). Because the latter bacteria are mainly acetate producers, it is not surprising that the coefficient of iso C14:0 and iso C15:0 is positive in the acetate-predicting equation. Linear odd-chain fatty acids were positively related with rumen proportions of propionate. During in vitro incubations, Vlaeminck et al. (2004) also found a positive relation between C15:0 and rumen proportions of propionate. This suggests that rumen conditions inducing high propionate production are favorable for bacteria synthesizing linear odd-chain fatty acids. Indeed, amylolytic bacteria (e.g., Ruminobacter amylophilus, Succinivibrio dextrinosolvens) show low levels of branched-chain fatty acids and are relatively enriched in linear odd-chain fatty acids (see, for example, Ifkovitz and Ragheb 1968; Minato et al., 1988). In addition, high rumen proportions of propionate might have stimulated 1) incorporation of propionyl-CoA in bacterial fatty acids (Fulco, 1983), and 2) de novo synthesis of linear odd-chain fatty acids from propionyl-CoA in the mammary gland (Vlaeminck et al., 2006a). These finding might also partially explain the positive relation between linear odd-chain fatty acids and rumen proportions of propionate.
Friggens et al. (1998) derived empirical equations based on observations in sheep to quantify the effect of dietary chemical components on rumen fermentation pattern. In their model development, a residual standard deviation of 11.4, 9.9, and 6.6 mmol/mol was reported for acetate, propionate, and butyrate, respectively. Recently, Bannink et al. (2006) estimated stoichiometric coefficients for VFA production from fermented substrates using literature data from experiments with dairy cows. The root mean square error values in their study were 32.6, 29.4, and 14.3 mmol/ mol for acetate, propionate, and butyrate, respectively. Results from the present study show intermediate SE values (15.3, 14.7, and 9.2 mmol/mol for acetate, propionate, and butyrate, respectively). This suggests that prediction of rumen fermentation pattern based on milk odd- and branched-chain fatty acids might be a promising alternative to equations based on dietary characteristics.
Model Validation
Although the relatively low SE supports the validity of the fitted regression model (Neter et al., 1996), comparison of MSPE values obtained with an independent data set is the best way to evaluate the accuracy of different equations. Due to the limited data available in the literature reporting both rumen fermentation pattern and milk odd- and branched-chain fatty acids, we first did a qualitative evaluation of the regression equations using milk odd- and branched-chain fatty acids reported by Kraft et al. (2003), Loor et al. (2005a, c), and Van Nespen et al. (2005). During this evaluation, prediction equations were solely evaluated on well-documented and expected changes in the rumen fermentation pattern because the latter was not measured in these studies.
Qualitative Evaluation.
Ruminal infusion of fish oil (Loor et al., 2005c) resulted in a predicted decrease of acetate whereas propionate increased. This shift in the rumen fermentation pattern toward an increased propionate production when supplementing fish oil is well documented (see Doreau and Chilliard, 1997; Fievez et al., 2003b). Similarly, diets high in rumen-degradable starch (Van Nespen et al., 2005) are known to decrease acetate and increase propionate proportions, as confirmed by the current predictions. Kraft et al. (2003) reported an increase in iso C14:0 and iso C15:0 in milk from cows of organic farms compared with conventional farms, resulting in a higher predicted acetate proportion and lower propionate. This was probably a reflection of the higher roughage proportion of the diets fed on organic farms (Kraft et al., 2003). In conclusion, results from this first qualitative evaluation show that predictions of rumen proportions of VFA based on milk odd- and branched-chain fatty acids were in line with the expected changes.
Quantitative Evaluation.
Validation of the prediction equations with the independent data set suggests that milk odd- and branched-chain fatty acids show potential to predict the rumen fermentation pattern. Indeed, the root MSPE as a percentage of the observed mean was less than 10%. In an evaluation of existing mechanistic models, Bannink et al. (1997) found that current models predict the molar proportions of individual VFA in the rumen of cattle inaccurately with a root MSPE ranging from 8.9 to 50.0 for molar proportions of acetate, 45.5 to 97.0 for propionate, and 15.0 to 53.2 for butyrate. Only the mechanistic model of Baldwin et al. (1987) showed better root MSPE values for molar proportions of acetate compared with the proposed predictions based on milk odd- and branched-chain fatty acids. Using newly estimated values of stoichiometric coefficients, Bannink et al. (2000) found a root MSPE of 24.1, 21.0, and 12.0 for molar proportions of acetate, propionate, and butyrate, respectively. Similar values were found by Nagorcka et al. (2000) using different stoichiometric coefficients for the 3 microbial groups currently represented in mechanistic models (root MSPE of 18.2, 31.3, and 21.0 for molar proportions of acetate, propionate, and butyrate, respectively). These values were slightly higher than the values observed using milk odd- and branched-chain fatty acids. However, an integrated comparison between the different models (mechanistic models vs. models based on milk odd- and branched-chain fatty acids) is hindered because of the different sets of data used. In addition, validations based on data from a few trials only may result in particularly favorable or unfavorable validations if trials near or far from the average trial used to define the model are used, respectively (Oldick et al., 1999).
The partition of MSPE into 3 components (mean bias, line bias, and random variation) suggests that random variation is the major cause of inaccuracy in the prediction of rumen VFA. A relatively high proportion of MSPE (>10%) was due to mean bias in the prediction of acetate and butyrate. The large amount of carbon exchange between acetate and butyrate (Sutton et al., 2003) might explain the underprediction of acetate and concomitant overprediction of butyrate. The biological consequence of this inaccurate prediction of the mean molar proportions of acetate and butyrate seems limited, because both are lipogenic VFA. Indeed, when molar proportions of lipogenic VFA (i.e., sum of acetate and butyrate) or the lipogenic:glucogenic ratio was calculated from the predictions of individual VFA, the root MSPE expressed as a percentage of the observed mean was 2.8 and 11.2% without significant mean and linear bias, despite the significant mean bias for predictions of acetate and butyrate.
The large deviation from the actual values observed for the linseed oil supplemented diets (Loor et al., 2005b) might indicate that the developed equations are not adequate when dairy cows are fed lipid supplements. Indeed, in none of the dietary treatments used for model development were lipids supplemented. The average fatty acid content of the diets used for model development was 22.4 ± 5.4 g/kg of DM (11.8 to 34.8 g/ kg of DM), a value far below the fatty acid content of the linseed supplemented diets (49 g/kg of DM; Loor et al., 2005b). The higher dietary fatty acid content might have reduced de novo synthesis of bacterial fatty acids or have resulted in a dilution with greater concentrations of even-chain fatty acids in milk. Hence, to improve prediction accuracy, dietary fatty acid content might have to be included in the regression equations. Furthermore, additional studies are needed to identify the effects of other factors, such as lactation stage, parity, and contribution of odd- and branched-chain fatty acids from body fat to milk odd- and branched-chain fatty acids, particularly during periods of negative energy balance.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Received for publication February 20, 2006. Accepted for publication May 22, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Craninx, A. Steen, H. Van Laar, T. Van Nespen, J. Martin-Tereso, B. De Baets, and V. Fievez Effect of Lactation Stage on the Odd- and Branched-Chain Milk Fatty Acids of Dairy Cattle Under Grazing and Indoor Conditions J Dairy Sci, July 1, 2008; 91(7): 2662 - 2677. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |