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


     


J. Dairy Sci. 2008. 91:4282-4292. doi:10.3168/jds.2007-0960
© 2008 American Dairy Science Association ®

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
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 Guo, J.
Right arrow Articles by Kohn, R. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Guo, J.
Right arrow Articles by Kohn, R. A.

Modeling Nutrient Fluxes and Plasma Ketone Bodies in Periparturient Cows

J. Guo, R. R. Peters and R. A. Kohn1

Department of Animal and Avian Sciences, University of Maryland, College Park 20742

1 Corresponding author: rkohn{at}umd.edu


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
A mechanistic model was developed to study the interrelationship between glucose and lipid metabolism in periparturient cows. The driving variables were dry matter intake, feed composition, calf birth weight, milk production, and milk components. The response variables were body fat content and concentrations of plasma glucose, glycerol, nonesterified fatty acids (NEFA), and total ketone bodies (KB). Fetal growth and milk synthesis were assigned the highest priority for glucose demand in the model. The rate of fat mobilization was expressed as a function of glucose deficiency. The model assumed first-order kinetics for utilization of NEFA and KB. Model prediction errors were 19, 43, 48, and 36% of mean predictions for glucose, glycerol, NEFA, and KB, respectively. A linear bias was observed in KB and glycerol predictions. The model may be useful for understanding and explaining ketosis development.

Key Words: mechanistic model • glucose metabolism • lipid metabolism • periparturient cow


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Ketosis causes economic losses in dairy herds directly by decreasing milk production and indirectly by increasing the risk for other periparturient diseases. Ketosis development results from the complicated interaction between glucose and lipid metabolism. Glucose deficiency may occur in periparturient cows because the demand for glucose is increased greatly by milk synthesis while DMI lags behind the requirements for milk production. Krebs (1966) proposed long ago that glucose deficiency depletes oxaloacetate, and that ketosis development results from the deficiency of oxaloacetate needed for the catabolism of acetyl-CoA. However, experimental results often cannot be explained by Krebs’ theory (Drackley et al., 2001) indicating that other interactions between glucose and lipid metabolism are also involved in ketogenesis. Glycerol from fat mobilization may be an important gluconeogenic precursor as the cow adapts to lactation (Drackley et al., 2001). Therefore, ketosis development in periparturient cows may result from excessive fat mobilization for glycerol to compensate for glucose deficiency.

A metabolic model of energy metabolism has been developed (Baldwin, 1995) and evaluated for the peri-parturient period (McNamara and Baldwin, 2000), but this model did not deal with the issue of ketone body accumulation. The objectives of the current study were 1) to quantitatively evaluate the relationship between nutrient fluxes and blood ketone bodies using quantities of blood metabolites and literature-derived rates of metabolite synthesis and utilization, and 2) to develop a mathematical model that may predict plasma ketone body concentrations from easily obtained data such as milk production and dry matter intake. Such a model could be used to demonstrate the theory behind ketosis, and might be used to predict the risk of its occurrence.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Data
The developmental data set used to parameterize the model was described previously (Guo et al., 2007). Twenty-eight multiparous Holstein cows were fed a nonlactating-cow diet prepartum and a lactation diet postpartum. Half of the cows were used as the control group. A transition diet was fed to the other half (the treatment group) for 17 d before calving date until 14 d postcalving. The nonlactating-cow diet contained 1.54 Mcal of NEL/kg, 10.9% CP, and 53.1% NDF, and the lactation diet 1.77 Mcal of NEL/kg, 16.8% CP, and 29.9% NDF. The transition diet (1.71 Mcal of NEL/kg, 16.8% CP, 35.2% NDF) had a lower energy density compared with the lactation diet.

The Model
A representation of the model is shown in Figure 1Go and Tables 1Go and 2Go. The driving variables were DMI, feed composition, calf birth weight, milk production, and milk components. The response variables were body fat content and the concentrations of plasma glucose, glycerol, NEFA, and KB (sum of BHBA, acetoacetate, and acetone).


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

 
Figure 1. A model of glucose and lipid metabolism in periparturient cows.

 

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

 
Table 1. Principal symbols used in the model
 

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

 
Table 2. Model equations
 
Response Variables.
The plasma glucose concentrations depended on the rate of glucose synthesis from glycerol, the rates of gluconeogenesis from propionate and amino acids, the rate of glucose utilization by peripheral tissues, and the rates of glucose utilization for fetal growth and milk synthesis. The model assumed that glycerol released from fat mobilization was completely utilized for glucose synthesis. Glucose utilized by peripheral tissues was expressed as a function of glucose concentration in plasma and metabolic body weight (Table 2Go).

The glycerol concentrations depended on fat mobilization and the rate of gluconeogenesis from glycerol. Glycerol used for lipogenesis in adipose tissues was not accounted for in the kinetics of glycerol concentrations, but was considered glucose consumption by peripheral tissues in the model. Total body fat (kg) depended on the rate of adipose mobilization (kg/d per kg of body fat) and glucose deficiency (ratio of glucose demand for fetal growth, milk synthesis, and peripheral tissue utilization to glucose supply from feed per day). The model assumed that lipolysis and lipogenesis were in equilibrium before d 7 prepartum and that lipogenesis had completely ceased in adipose tissues after d 7 prepartum. The change in plasma NEFA was expressed as a balance between its utilization and fat mobilization. The model assumed that the rate of NEFA utilization (mmol/d per mmol of NEFA) followed first-order kinetics. The KB concentrations depended on the ketogenesis rates from NEFA and butyrate, and the rate of KB utilization (mmol/d per mmol of KB). The rate of KB utilization was represented as a first-order reaction in the model.

Parameters Derived from Published Literature.
The rate of propionate production from rumen fermentation was calculated from DMI and feed composition (Murphy et al., 1982). More recent equations for predicting VFA production (Dijkstra et al., 1992; Bannink et al., 2000) were not chosen because those equations require more inputs such as rumen volume and fractional passage rates. The rumen adaptation index in the model represented the shift to grain fermentation in rumen from the high-forage diet of the nonlactating period. The propionate-producing microbes were assumed to increase when the transition or lactation diet was fed. The shift depends mainly on 2 factors: substrate availability and the microbial population. The index was derived from the logistic growth equation (France and Thornley, 1984): dW/dt = k x W x S, where W is microbial population at time of t; k is a growth rate constant; S is substrate availability for microbial growth, which after integration and rearrangement gives W = W0 x Wf /[W0 + (Wf – W0)–ut], where W0 = 1; Wf = 100; u = 0.95. The values for W0 and Wf were assigned only to represent the change in propionate-producing microbes. The coefficient u corresponds to the rate of adaptation of rumen microbes to the high-grain diet, and it was assigned the value 0.95, which corresponds to a 10-d adaptation period. The 10-d adaptation period is based on the fact that the potential for digestive upsets from rapid fermentation of starch requires gradual introduction of high concentrate diets over a period of at least 10 d (Coe et al., 1999).

All of the propionate produced from starch in the rumen was assumed to be available for gluconeogenesis in the animal. In reality, the rumen epithelium could metabolize 40 to 60% of propionate (Parker, 1990). Products of metabolism are energy, CO2, lactate, and alanine, the last 2 of which are glucogenic. The liver could take up 90% or more of the absorbed propionate (Kristensen and Harmon, 2004). Although this assumption regarding propionate available for gluconeogenesis would lead to an overestimate, we are likely to have underestimated glucose available from the digestion of bypass starch in the small intestine. At least 90% of starch in small grains is fermentable in the rumen, up to 40% of corn starch could escape rumen fermentation (Orskov, 1986). Quantities of starch disappearing in the small intestine vary with digestibilities ranging from 10 to 96% (Harmon, 1992). The variability of the available information makes it difficult to incorporate the above quantitative information into the model. For model simplicity, the contribution of the escape starch to glucose production was not considered and all the propionate was assumed to be converted into glucose by the liver in the model. The rate of butyrate production was estimated in the same way as the rate of propionate production in the model (Table 2Go).

The rate of gluconeogenesis from protein was estimated by the amount of catabolized protein. Catabolism of 100 g of protein was assumed to give rise to 58 g of glucose (Dukes, 1993). Before parturition, one-third of metabolizable protein was catabolized (NRC, 2001) assuming that all the metabolizable protein came from feed. After parturition some of the catabolized protein came from endogenous sources (Komaragiri et al., 1998). The rate of catabolized protein postpartum was estimated from urinary nitrogen excretion as Urinary N (g/d) = 0.026 x BW (kg) x MUN (mg/dL) (Kohn et al., 2002). The equation is based on renal physiology throughout lactation. However, overestimation of postpartum protein catabolism may result from origination of urinary nitrogen from rumen ammonia. Carbohydrate sources and ratios of carbohydrate to protein fractions in the diets could influence the amount of urea from rumen ammonia. The coefficient (58 g of glucose per 100 g of protein) is also associated with large variation (Dukes, 1993). Estimation of gluconeogenesis from protein could be improved once additional information is available.

The rate of glucose utilization for fetal growth was derived from the energy requirement, which was estimated from calf birth weights according to NRC (2001). In the nonlactating period, about 775 of 2,336 kcal/d of energy required for fetal growth came from glucose and lactate (Bell, 1995). Presumably, one-third of the energy requirement was provided in the form of glucose.

The mammary glucose requirement was estimated from milk yield and lactose concentrations. Lactose production in milk accounted for 55 to 78% of the glucose uptake by the mammary gland (Cant et al., 1993; Mackle et al., 2000). In the model, the efficiency of glucose uptake for lactose synthesis was assumed to be 66.5%, which is the average of 55 and 78%.

Parameters Determined by Best Fit.
The parameters listed in Table 3Go were optimized by a modified Powell algorithm for each individual cow to find the minimum sum of squares for deviations between observed data and model predictions. The model predictions for the response variables were calculated using the DMI, feed composition, calf birth weight, milk yield, and milk components from the development data set as the model input. The initial values for the response variables were adapted from the pretreatment measurements from the developmental data set. In the developmental data, body fat content was not actually measured, but was estimated from body weight and BCS (NRC, 2001).


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

 
Table 3. The mean values of the parameters (SD)1 (n = 28)
 
The model was simulated and fit using Scientist software (Scientist Software, 1995, Salt Lake City, UT). Two criteria (Neal and Thornley, 1983) were applied in fitting parameters to the model. First, the predictions of the model must approximate observed data. Second, biological realism should not be violated.

Statistical Analysis
The mean parameter values for all cows (Table 3Go) in the developmental data were used to calculate model predictions, and compared against observations. A mean bias for model predictions was declared if residual (observed – predicted) values were significantly different from zero. Linear bias for model predictions was evaluated by regression analysis for residuals against model predictions. Root mean square prediction error (RMSPE) was calculated from the following equation (Bibby and Toutenburg, 1977): RMSPE = square root of [{sum} (observed – predicted)2/n]. Statistical significance was declared at P < 0.05.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Parameter values are shown in Table 3Go. The rates of fat mobilization and KB utilization postpartum were significantly greater than those prepartum (P < 0.05 and P < 0.01, respectively).

The predicted values for plasma glucose, glycerol, NEFA, and KB followed a similar pattern to the data observed in the animal trial (Figure 2Go). The agreement between predicted and observed KB was further demonstrated in a ketotic cow as an extreme case. One cow was diagnosed with clinical ketosis and her data were excluded from data analysis in the previous paper (Guo et al., 2007) and from model parameterization in the present paper. At d 2 and 3 postpartum, she was administrated intravenously with 1,000 mL of dextrose (50%). The KB concentration increased from 0.6 to 0.9 mM, and then decreased from 0.9 to 0.25 mM after the dextrose treatment had stopped. The agreement between observed and predicted values for KB concentrations is presented in Figure 3Go.


Figure 2
View larger version (17K):
[in this window]
[in a new window]

 
Figure 2. Time course of blood metabolites during periparturient period. {blacktriangleup} = observed values in the control group; {blacksquare} = observed values in the treatment group; ------- = predicted values for the control group; —— predicted values for the treatment group. The cows in the control group were fed a diet for dry cows before calving and a lactation diet postcalving. The cows in the treatment group were fed a transition diet in the last 17 d of gestation and the first 14 d of lactation.

 

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

 
Figure 3. Time course of total plasma ketone bodies in a ketotic cow during periparturient period. {blacksquare} = observed; —— = model predicted. The cow was diagnosed with clinical ketosis and given 1,000 mL of 50% dextrose at d 2 and 3 postpartum. Data around parturition were missing due to severe weather conditions.

 
Residual analysis for the developmental data was conducted by comparing model predictions with residuals (observed – predicted) for the response variables (Figures 4Go and 5Go). The model overpredicted glycerol concentrations with a mean bias of 0.001 mM (P < 0.05; Table 4Go). No mean bias or linear bias was found for body fat predictions (P > 0.05). Linear biases for glucose, glycerol, NEFA, and KB were observed (P < 0.05). The standard error of model predictions was 19, 43, 48, 36, and 4% of mean predictions for glucose, glycerol, NEFA, KB, and body fat predictions.


Figure 4
View larger version (19K):
[in this window]
[in a new window]

 
Figure 4. Plots of residual values against model predictions. Residual = observed – predicted. Data points with {blacktriangleup} were identified as outliers and were not used by regression analysis. The outliers were detected with 3 regression diagnostic statistics: jackknife residuals, leverages, and Cook’s distance (Kleinbaum et al., 1998). Glucose: y = 0.65–0.18 x x; glycerol: y = 0.006 – 0.313 x x; NEFA: y = 0.007 – 0.054 x x; ketone bodies: y = 0.21 – 0.44 x x.

 

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

 
Figure 5. Plot of residual values against model predicted body fat. Residual = observed – predicted. y = 3.48 – 0.06 x x.

 

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

 
Table 4. Residual analysis for the developmental data set (from 28 cows)
 
According to the model prediction, peripheral tissues consumed more glucose in the last 21 d of gestation compared with the first 21 DIM with a surge around parturition. Glucose consumed by peripheral tissues was greater for the treatment cows from d 15 prepartum to d 10 postpartum compared with the control cows (Figure 6Go). The predicted glucose balances differed in the patterns between the treatment and control groups (Figure 7Go). Before parturition, the treatment cows had a greater glucose balance. However, after parturition the control cows had a greater glucose balance compared with the treatment cows. According to the model prediction, glycerol provided 12 and 17% of the glucose demand in the control and treatment groups, respectively, and the treatment group mobilized more glycerol for gluconeogenesis from adipose tissues after parturition compared with the control (Figure 8Go).


Figure 6
View larger version (13K):
[in this window]
[in a new window]

 
Figure 6. Model prediction for glucose consumed by the peripheral tissues. {blacktriangleup} = control group; {blacksquare} = treatment group. The cows in the control group were fed a diet for dry cows before calving and a lactation diet postcalving. The cows in the treatment group were fed a transition diet in the last 17 d of gestation and the first 14 d of lactation.

 

Figure 7
View larger version (14K):
[in this window]
[in a new window]

 
Figure 7. Estimated glucose balances during the periparturient period. {blacktriangleup} = control group; {blacksquare} = treatment group. The cows in the control group were fed a diet for dry cows before calving and a lactation diet postcalving. The cows in the treatment group were fed a transition diet in the last 17 d of gestation and the first 14 d of lactation. Glucose balances were calculated by propionate from feed, catabolized protein, fetus growth, milk synthesis, and peripheral tissue requirement.

 

Figure 8
View larger version (13K):
[in this window]
[in a new window]

 
Figure 8. Model prediction for the contribution of glycerol from fat mobilization to glucose synthesis in the periparturient cows. {blacktriangleup} = control group; {blacksquare} = treatment group. The cows in the control group were fed a diet for dry cows before calving and a lactation diet post-calving. The cows in the treatment group were fed a transition diet in the last 17 d of gestation and the first 14 d of lactation.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Physiological Basis of the Model
The model assumed that glucose was the limiting nutrient during the periparturient period. In ruminants, glucose supply is met mainly by gluconeogenesis, an energetically inefficient pathway compared with hydrolysis of starch in nonruminant animals. Fetal growth increased dramatically in late gestation (House and Bell, 1993), and oxidation of lactate and glucose account for up to 60% of fetal respiration in lambs (Hay et al., 1983). After parturition, the demand for glucose is increased greatly by milk synthesis. Lactate utilization for gluconeogenesis represents recycling of carbon because some of circulating lactate is formed during catabolism of glucose by peripheral tissues (Drackley et al., 2001). Lactate from partial catabolism of propionate by visceral epithelial tissues was accounted for by the assumption that propionate was completely used for gluconeogenesis in the model. The primary function of body protein is by no means to preserve glucose precursors, although amino acids (except for leucine and lysine) are glucogenic. Between 2 and 5 wk postpartum cows mobilized only 12 kg of body protein (Komaragiri et al., 1998), while part of mobilized protein must be excreted as milk protein. Bauman and Elliot (1983) concluded that over the period up to peak lactation the contribution of mobilized tissue protein to gluconeogenesis is small. Depressed feed intake around parturition (Ingvartsen and Andersen, 2000) could further exaggerate glucose deficiency. Thus, glucose deficiency could play a critical role in the orchestration of the entire metabolism in periparturient cows.

Glucose
Glucose utilization by peripheral tissues is regulated mainly by plasma insulin, tissue responses to insulin (Petterson et al., 1993), and glucose availability. The exact mathematical relationship between glucose utilization and these 3 factors was beyond the scope of this model. An aggregate equation was selected to approximately represent that relationship in the model (Table 2Go). In the equation, the dependence of peripheral consumption on glucose availability was defined as the ratio of plasma glucose concentrations to a reference concentration (3.15 mM). The power ‘Q’ to that ratio was assigned to 3 prepartum and to 4 postpartum, which provided the best fit of the data. The biological meaning of the power ‘Q’, which is not known, is presumably related to endocrine status and tissue sensitivity to hormones. Although the equation of glucose utilization by peripheral tissues did not fully represent the biological mechanism, the model predicted blood glucose concentrations without mean bias. The linear bias, although statistically significant, could not be of biological importance because the maximum prediction error resulting from the linear bias is only 0.31 mM [–0.31 = 0.65 – 0.18 x 5.43 – (–0.02), where 5.43 is the maximum glucose concentration observed in the development data set, –0.02 is the mean bias, and the rest of equation is from the regression analysis in Figure 3Go].

The model predicted that more glucose was consumed by peripheral tissues prepartum compared with the postpartum period, with a surge around parturition. A difference in peripheral glucose utilization between the treatment and control groups was also predicted by the model. Before parturition the difference may result from the greater glucose availability in the treatment groups compared with the control group as reported previously (Guo et al., 2007). Feeding a high-concentrate diet during the late gestation period increased plasma insulin concentrations, and this effect carried over into early lactation (Holcomb et al., 2001). The carryover effect may be responsible for the difference between the treatment and control groups after parturition, as the model predicted. According to the model predictions, the amount of glucose consumed by peripheral tissues ranged from 0.04 to 0.07 mol/d per kg of BW0.75 prepartum, and from 0.04 to 0.05 mol/d per kg of BW0.75 postpartum in the control cows. A turnover rate of glucose in ruminants under basal conditions had been reported by Baldwin (1995) between 0.03 and 0.05 mol/d per kg of BW0.75. Compared with the data by Baldwin (1995), the overprediction may be caused by 2 factors: first, the cows were not under the basal condition; second, during the nonlactating period, glucose may be used to produce glycerol moiety for triglyceride synthesis in adipose tissues.

Glycerol
A mean bias was observed for glycerol concentration predictions; however, the absolute value was only 0.001 mM, which is biologically insignificant relative to the glycerol concentrations under the basal condition. A linear bias was also observed for glycerol predictions. The linear bias was mainly caused by a few of the residuals when glycerol predictions were above the normal range. Blood glycerol can be directly incorporated into galactose for lactose synthesis in the mammary gland (Sunehag et al., 2002). The linear bias for glycerol predictions may be caused by the synthesis of galactose from glycerol, which was not considered in the model. Contribution of glycerol to gluconeogenesis ranged from 12 to 17% as predicted by the current model, which agrees well with the range from 15 to 20% of the glucose demand at 4 d postpartum (Bell, 1995).

Removal of glycerol by the liver (Reynolds et al., 2003) and gluconeogenesis (Greenfield et al., 2000) increases greatly after parturition. The model assumed that glycerol released from fat mobilization was completely utilized for glucose synthesis. Contribution of glycerol to gluconeogenesis depended on the concentrations of plasma glycerol that resulted from fat mobilization. Fat mobilization was determined by glucose deficiency in the model. The model predictions for glycerol concentrations agreed well with the observed values from the developmental data (Figure 2Go).

Body Fat Content
The result of model parameterization showed that the postpartum rate of fat mobilization was greater than the prepartum rate. The difference in the rates of fat mobilization prepartum and postpartum is in agreement with the fact that the change in the endocrine profiles and sensitivities to hormones greatly enhanced fat mobilization in early lactation compared with late gestation (Bell, 1995).

The contribution of glucose precursor from fat mobilization is predicted in Figure 8Go. According to the model, 270 g of glucose could be provided from 2.5 kg of fat when maximum fat mobilization occurred. This amount of glucose is small relative to total supply from feed; however, the importance of this amount should not be inconsequential. In the field, bolus i.v. administration of 500 mL of 50% dextrose solution (about 250 g) is a common therapy for ketosis. Drenching of propylene glycol (250 g/d) is also effective in treating ketosis for postpartum cows. These small doses provide a similar amount of glucose as the amount of glycerol from fat mobilization, but have a consequential effect on reestablishing glycemia and reducing blood ketone body concentrations in ketotic cows. The present model demonstrated that the contribution of glucose precursor from fat mobilization should be considered to explain the dynamics of glucose metabolism in peri-parturient cows after the glucose supply from feed and gluconeogenesis from body protein mobilization has been accounted for. Thus, glycerol from fat mobilization may be an important gluconeogenic precursor as the cow adapts to lactation (Drackley et al., 2001).

The actual amount of fat mobilization could be underestimated in the present model. The observed values for body fat content used for model parameterization were estimated from BW and BCS. Changes in BW after parturition could easily be masked by the mass increase in the gastrointestinal tract and its contents. Variation in body condition would not be detected within a short period by BCS, which is a crude estimate. In addition, calculation of glucose precursors from fat mobilization is also underestimated. According to the model, glycerol from 1 kg of triglyceride provides about 110 g of glucose, which means that an extra 1 kg of milk requires 3 kg more fat to be mobilized. In fact, oxidation of odd- and branched-chain fatty acids produces propionyl-CoA, which could be either completely oxidized or converted into glucose. Quantitative information on metabolic fate of the propionyl-CoA is limited. The total proportions of odd-number fatty acids represented 3.7 to 4.2% of total fatty acids (Elias-Calles et al., 1997), and the proportion of branched-chain fatty acids was about 2% (Duncan and Garton, 1978). Thus, it is difficult to incorporate metabolism of odd- and branched-chain fatty acids into the current model.

NEFA
Endocrine status differs greatly before and after parturition. However, according to the model parameterization, the prepartum rate of NEFA utilization was not significantly different from the postpartum rate, probably because NEFA metabolism in the liver of dairy cows is less responsive to hormonal control than is the metabolism in laboratory species (Cadorniga-Valino et al., 1997), and there is no preferential direction of NEFA to the liver at the expense of other tissue (Drackley et al., 2001). Uptake and oxidation of NEFA by the liver and extrahepatic tissues are directly related to plasma concentrations (Pethick et al., 1983). Thus, the model assumed the rate of NEFA utilization followed first-order kinetics and expressed in term of millimoles per day per millimole of NEFA. In the model, the postpartum rate of NEFA utilization rate was 0.2236 (mmol/d per mmol of NEFA). With this value, the predicted whole-body utilization for the cows at wk 4 postpartum studied by Reynolds et al. (1988) would be 4.4 mol/d: 4,400 mmol/d = 0.2236 (mmol/d per mmol, NEFA utilization rate) x 0.328 (mM, NEFA concentration) x 472 (L/BW0.75, cardiac output) x 6400.75 (kg of BW0.75). Given that the percentage of cardiac output flowing through the liver was 31% (Huntington et al., 1990) and that the measured net uptake of NEFA by liver was 60.8 mmol/h (Reynolds et al., 1988), the NEFA utilization by the whole body would be 4.7 mol/d (4,700 mmol/d = 60.8 mmol/h x 24 h/d ÷ 31%), which agreed well with the current model prediction. The linear bias, although statistically significant, could not be of biological importance as the maximum prediction error resulting from the linear bias was only 0.09 mM [–0.09 = 0.007 – 0.054 x 2.05 – (–0.01)], where 2.05 is the maximum NEFA concentration observed in the development data set, –0.01 is the mean bias, and the rest of equation is from the regression analysis in Figure 3Go.

Ketone Bodies
In the model, the prepartum rate of KB utilization was significantly different from the postpartum rate. The model predicted that the pre- and postpartum rates of KB utilization rates were 0.31 and 0.47 mmol/d per mmol KB respectively, corresponding to 7.5 mol/d in late gestation and to 17.5 mol at d 5 postpartum in the control cows. Extrapolated from the published data in sheep when corrected for BW, the cows in the control group would utilize 4.6 mol of KB per day in late gestation, and 12.5 mol of KB at d 5 postpartum according to Heitmann et al. (1987). The utilization rates of KB in the model were slightly greater compared with the published data, probably because the data by Heitmann et al. (1987) did not include acetone metabolism. The quantitative information on acetone metabolism is extremely limited for dairy cows. In rats, after administration of radiolabeled acetone by stomach tube or by injection, demonstrable amounts of radioactivity were recovered in glycogen, urea, cholesterol, fatty acids, amino acids, and heme, and a substantial amount of radiolabeled carbon was recovered in exhaled carbon dioxide (Price and Rittenberg, 1950).

There is a significant linear bias for KB predictions in the model. The bias may result from acetone metabolism as discussed above. In addition, as KB concentrations increase, an increased proportion of KB is lost in the urine or via breathing. The ratio of KB excreted via urine (mg/d) to blood concentration (mg/dL) is <42 in humans under normal conditions, whereas the ratio is about 56 in an untreated diabetic patient (Nelson et al., 2000). Another possibility may come from peroxisomal oxidation, which is not considered in the model. The peroxisomal pathway, an auxiliary pathway to mitochondrial oxidation, may be induced when hepatocellular influx of NEFA is increased (Drackley et al., 2001).

Although the current model represented glucose and lipid metabolism under normal conditions, the KB profile after therapeutic glucose infusion was successfully simulated in a ketotic cow. Clinical ketosis occurs when there is a failure of the homeostatic mechanisms regulating the glucose and fat metabolism. The therapeutic approach is to reestablish normal homeostasis (Herdt and Emery, 1992). The agreement between model prediction and KB profile after therapeutic treatment further supports the proposed interrelationship between glucose and lipid metabolism simulated in the model. However, only one cow was diagnosed with ketosis as presented in Figure 3Go, and more animals are needed to evaluate the applicability of the model under ketotic conditions in future studies.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Using DMI, feed composition, calf birth weight, milk yield, and milk components as driving variables, the model can explain most of the variations in body fat content, and plasma concentrations of glucose, glycerol, NEFA, and KB during the periparturient period. When using this model to quantify metabolite flows, glucose deficiency was closely related to the rate of fat mobilization. The excessive KB could result from elevated fat mobilization due to low plasma and tissue glucose concentrations. The model should be evaluated with data from independent experiments and compared with other models. In addition, more animal experiments are needed to investigate the possibility that fat is mobilized because of low glucose concentrations in periparturient cows. The model may be useful for understanding and explaining ketosis development.

Received for publication December 18, 2007. Accepted for publication July 6, 2008.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 


Baldwin, R. L. 1995. Modeling Ruminant Digestion and Metabolism. 1st ed. Chapman & Hall, London, UK.

Bannink, A., J. Kogut, J. Kijkstra, J. France, S. Tamminga, and A. M. van Vuuren. 2000. Modelling production and portal appearance of volatile fatty acids in dairy cows. CAB International, Wallingford, UK.

Bauman, D. E., and J. M. Elliot. 1983. Control of nutrient partitioning in lactating ruminants. Elsevier, Amsterdam, the Netherlands.

Bell, A. W. 1995. Regulation of organic nutrient metabolism during transition from late pregnancy to early lactation. J. Anim. Sci. 73:2804–2819.[Abstract]

Bibby, J., and H. Toutenburg. 1977. Prediction and Improved Estimation in Linear Models. Wiley, London, UK.

Cadorniga-Valino, C., R. R. Grummer, L. E. Armentano, S. S. Donkin, and S. J. Bertics. 1997. Effects of fatty acids and hormones on fatty acid metabolism and gluconeogenesis in bovine hepatocytes. J. Dairy Sci. 80:646–656.[Abstract/Free Full Text]

Cant, J. P., E. J. DePeters, and R. L. Baldwin. 1993. Mammary amino acid utilization in dairy cows fed fat and its relationship to milk protein depression. J. Dairy Sci. 76:762–774.[Abstract/Free Full Text]

Coe, M. L., T. G. Nagaraja, Y. D. Sun, N. Wallace, E. G. Towne, K. E. Kemp, and J. P. Hutcheson. 1999. Effect of virginiamycin on ruminal fermentation in cattle during adaptation to a high concentrate diet and during an induced acidosis. J. Anim. Sci. 77:2259–2268.[Abstract/Free Full Text]

Dijkstra, J., H. D. Neal, D. E. Beever, and J. France. 1992. Simulation of nutrient digestion, absorption and outflow in the rumen: Model description. J. Nutr. 122:2239–2256.[Abstract/Free Full Text]

Drackley, J. K., T. R. Overton, and G. N. Douglas. 2001. Adaptations of glucose and long-chain fatty acid metabolism in liver of dairy cows during the periparturient period. J. Dairy Sci. 84(E Suppl.):100–112.[Abstract]

Dukes, H. H. 1993. Physiology of Domestic Animals. Comstock Publishing Associates, Ithaca, NY.

Duncan, W. R. H., and G. A. Garton. 1978. Differences in the proportions of branched-chain fatty acids in subcutaneous triacylglycerols of barley-fed ruminants. Br. J. Nutr. 40:29–33.[CrossRef][Medline]

Elias-Calles, J. A., C. T. Gaskins, J. R. Busboom, J. D. Cronrath, S. K. Duckett, J. J. Reeves, R. W. Wright, K. A. Johnson, and C. J. Moreshead. 1997. Comparison of fatty acid composition and palatability characteristics of beef from various Wagyu sires. WSU Wagyu Proc. 3:63–86.

France, J., and J. H. Thornley. 1984. Growth functions. Pages 75–94 in Mathematical Models in Agriculture: A quantitative approach to problems in agriculture and related sciences. Butterworths, London, UK.

Greenfield, R. B., M. J. Cecava, T. R. Johnson, and S. S. Donkin. 2000. Impact of dietary protein amount and rumen undegradability on intake, peripartum liver triglyceride, plasma metabolites, and milk production in transition dairy cattle. J. Dairy Sci. 83:703–710.[Abstract]

Guo, J., R. R. Peters, and R. A. Kohn. 2007. Effect of transition diet on production performance and metabolism in periparturient dairy cows. J. Dairy Sci. 90:5247–5258.[Abstract/Free Full Text]

Harmon, D. L. 1992. Dietary influences on carbohydrases and small intestinal starch hydrolysis capacity in ruminants. J. Nutr. 122:203–210.[Abstract/Free Full Text]

Hay, W. W. Jr., S. A. Myers, J. W. Sparks, R. B. Wilkening, G. Meschia, and F. C. Battaglia. 1983. Glucose and lactate oxidation rates in the fetal lamb. Proc. Soc. Exp. Biol. Med. 173:553–563.[CrossRef][Medline]

Heitmann, R. N., D. J. Dawes, and S. C. Sensenig. 1987. Hepatic ketogenesis and peripheral ketone body utilization in the ruminant. J. Nutr. 117:1174–1180.[Abstract/Free Full Text]

Herdt, T. H., and R. S. Emery. 1992. Therapy of diseases of ruminant intermediary metabolism. Vet. Clin. North Am. Food Anim. Pract. 8:91–106.[Medline]

Holcomb, C. S., H. H. Van Horn, H. H. Head, M. B. Hall, and C. J. Wilcox. 2001. Effects of prepartum dry matter intake and forage percentage on postpartum performance of lactating dairy cows. J. Dairy Sci. 84:2051–2058.[Abstract]

House, W. A., and A. W. Bell. 1993. Mineral accretion in the fetus and adnexa during late gestation in Holstein cows. J. Dairy Sci. 76:2999–3010.[Abstract/Free Full Text]

Huntington, G. B., J. H. Eisemann, and J. M. Whitt. 1990. Portal blood flow in beef steers: Comparison of techniques and relation to hepatic blood flow, cardiac output and oxygen uptake. J. Anim. Sci. 68:1666–1673.[Abstract]

Ingvartsen, K. L., and J. B. Andersen. 2000. Integration of metabolism and intake regulation: A review focusing on periparturient animals. J. Dairy Sci. 83:1573–1597.[Abstract]

Kleinbaum, D. G., L. L. Kupper, K. E. Muller, and Z. Nizam. 1998. Regression diagnostics. Page 228–233 in Applied Regression Analysis and Multivariable Methods. 3rd ed. Brooks/Cole Publishing Company, Pacific Grove, CA.

Kohn, R. A., K. F. Kalscheur, and E. Russek-Cohen. 2002. Evaluation of models to estimate urinary nitrogen and expected milk urea nitrogen. J. Dairy Sci. 85:227–233.[Abstract]

Komaragiri, M. V., D. P. Casper, and R. A. Erdman. 1998. Factors affecting body tissue mobilization in early lactation dairy cows. 2. Effect of dietary fat on mobilization of body fat and protein. J. Dairy Sci. 81:169–175.[Abstract]

Krebs, H. A. 1966. Bovine ketosis. Vet. Rec. 78:187–192.[Medline]

Kristensen, N. B., and D. L. Harmon. 2004. Splanchnic metabolism of volatile fatty acids absorbed from the washed reticulorumen of steers. J. Anim. Sci. 82:2033–2042.[Abstract/Free Full Text]

Mackle, T. R., D. A. Dwyer, K. L. Ingvartsen, P. Y. Chouinard, D. A. Ross, and D. E. Bauman. 2000. Evaluation of whole blood and plasma in the interorgan supply of free amino acids for the mammary gland of lactating dairy cows. J. Dairy Sci. 83:1300–1309.[Abstract]

McNamara, J. P., and R. L. Baldwin. 2000. Estimation of parameters describing lipid metabolism in lactation: Challenge of existing knowledge described in a model of metabolism. J. Dairy Sci. 83:128–143.[Abstract]

Murphy, M. R., R. L. Baldwin, and L. J. Koong. 1982. Estimation of stoichiometric parameters for rumen fermentation of roughage and concentrate diets. J. Anim. Sci. 55:411–421.[Abstract/Free Full Text]

National Research Council. 2001. Nutrient Requirements of Dairy Cattle. 7th rev. ed. National Academy Press, Washington, DC.

Neal, H., and J. H. Thornley. 1983. The lactation curve in cattle: A mathematical model of the mammary gland. J. Agric. Sci. Camb. 101:389–400.

Nelson, D. L., M. M. Cox, and A. L. Lehninger. 2000. Lehninger Principles of Biochemistry. 3rd ed. Worth Publishers, New York, NY.

Orskov, E. R. 1986. Starch digestion and utilization in ruminants. J. Anim. Sci. 63:1624–1633.[Abstract/Free Full Text]

Parker, D. S. 1990. Manipulation of the functional activity of the gut by dietary and other means (antibiotics/probiotics) in ruminants. J. Nutr. 120:639–648.[Abstract/Free Full Text]

Pethick, D. W., D. B. Lindsay, P. J. Barker, and A. J. Northrop. 1983. The metabolism of circulating non-esterified fatty acids by the whole animal, hind-limb muscle and uterus of pregnant ewes. Br. J. Nutr. 49:129–143.[CrossRef][Medline]

Petterson, J. A., F. R. Dunshea, R. A. Ehrhardt, and A. W. Bell. 1993. Pregnancy and undernutrition alter glucose metabolic responses to insulin in sheep. J. Nutr. 123:1286–1295.[Abstract/Free Full Text]

Price, T. D., and D. Rittenberg. 1950. The metabolism of acetone. I. Gross aspects of catabolism and excretion. J. Biol. Chem. 185:449–459.[Free Full Text]

Reynolds, C. K., P. C. Aikman, B. Lupoli, D. J. Humphries, and D. E. Beever. 2003. Splanchnic metabolism of dairy cows during the transition from late gestation through early lactation. J. Dairy Sci. 86:1201–1217.[Abstract/Free Full Text]

Reynolds, C. K., G. B. Huntington, H. F. Tyrrell, and P. J. Reynolds. 1988. Net metabolism of volatile fatty acids, D-beta-hydroxybutyrate, nonesterified fatty acids, and blood gasses by portal-drained viscera and liver of lactating Holstein cows. J. Dairy Sci. 71:2395–2405.[Abstract/Free Full Text]

Scientist Handbook. 1995. Experimental data fitting; Microsoft windows version. SofMicroMath Inc., Salt Lake City, UT.

Sunehag, A. L., K. Louie, J. L. Bier, S. Tigas, and M. W. Haymond. 2002. Hexoneogenesis in the human breast during lactation. J. Clin. Endocrinol. Metab. 87:297–301.[Abstract/Free Full Text]


This article has been cited by other articles:


Home page
J DAIRY SCIHome page
J. Guo, R. R. Peters, and R. A. Kohn
Evaluation of a Mechanistic Model of Glucose and Lipid Metabolism in Periparturient Cows
J Dairy Sci, November 1, 2008; 91(11): 4293 - 4300.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
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 Guo, J.
Right arrow Articles by Kohn, R. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Guo, J.
Right arrow Articles by Kohn, R. A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS