If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Heat production (HP) represents a major energy cost in lactating dairy cows. Better understanding of factors that affect HP will improve our understanding of energy metabolism. Our objective was to derive models to explain variation in HP of lactating Jersey cows. Individual animal-period data from 9 studies (n = 293) were used. The data set included cows with a wide range (min to max) in days in milk (44–410) and milk yield (7.8–43.0 kg/d). Diets included corn silage as the predominate forage source, but diets varied (min to max on DM basis) in crude protein (CP; 15.2–19.5%), neutral detergent fiber (NDF; 35.5–43.0%), starch (16.2–31.1%), and crude fat (2.2 to 6.4%) contents. Average HP was (mean ± standard deviation) 22.1 ± 2.86 Mcal/d, or 28.1 ± 3.70% of gross energy intake. Eight models were fit to explain variation in HP: (1) dry matter intake (DMI; INT); (2) milk fat, protein, and lactose yield (MILKCOMP); (3) INT and milk yield (INT+MY); (4) INT and MILKCOMP/DMI (INT+MILKCOMP); (5) mass of digested NDF, CP, and starch (DIG); (6) INT and digested energy (INT+DE); (7) INT and NDF, CP, and starch digestibility (INT+DIG); or (8) INT+MILKCOMP model plus urinary N excretion (INT+MILKCOMP+UN). For all HP models, metabolic body weight was included. All models were derived via a backward elimination approach and included the random effects of study, cow, and period within block within study. The INT models adequately explained variation in HP with a nonrandom effect–adjusted concordance correlation coefficient of 0.84. Similar adjusted concordance correlation coefficients (0.79–0.85) were observed for other HP models. The HP associated with milk protein yield and supply of digestible protein was greater than other milk production and nutrient digestibility variables. The HP associated with urinary N excretion was 5.32. Overall, HP can be adequately predicted from metabolic body weight and DMI. Milk component yield, nutrient digestibility, or urinary N excretion explained similar variation as DMI. Coefficients for milk protein and protein digestion suggest that digestion and metabolism of protein and synthesis of milk protein contribute substantially to HP of a dairy cow.
Heat in animals is generated from the catabolism of organic compounds. According to the laws of thermodynamics, heat generated from reactions in the body is equal to the difference in chemical energy between all substrates and products in a reaction (
). Heat production can be categorized into NEM—the basal energy expenditure of an animal—and heat increment (HI)—the HP that is associated with digestion and metabolism of food (
). Because HI represents a loss of energy and is the difference between ME and NEL, changes in HI may affect net milk production and tissue energy accretion of a dairy cow. Energetic efficiency, HP, and HI are affected by plane of production (
Effects of diet forage-to-concentrate ratio and intake on energy metabolism in growing beef heifers: Whole body energy and nitrogen balance and visceral heat production.
Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
). A recent analysis of indirect calorimetry data from the 1960s through the 1990s reported that kL ranged from 0.60 to 0.70 and was associated with heart rate, milk yield, and dietary crude fat (
). Therefore, the use of a static kL value may misrepresent underlying understanding of energy transformations. Consequently, the estimation of dietary NEL may improve with models that predict HP by accounting for sources of energy substrates and end products formed. The objective of this study was to derive models that explain variation in HP of lactating Jersey cows. We hypothesized that HP could be estimated from MBW and DMI, and that the addition of variables that represent the nature of nutrients contained in feed consumed would further explain variation by better representing the biological basis of energy transformations.
MATERIALS AND METHODS
Data Collection
Data from 9 research experiments conducted at the University of Nebraska–Lincoln's Dairy Metabolism Unit were used (
Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
). Across these 9 experiments, data from 293 animal-periods were collected on 54 Jersey cows. However, digestible energy (DE), and nutrient digestibility were not measured on all observations (n = 261–277 for the data set with DE and nutrient digestibility). In 2 studies (
, nutrient digestibility was measured, but, of the energy fractions variables, only HP was measured. Descriptive statistics for variables are listed in Table 1.
Table 1Descriptive statistics of the data used for the analysis to estimate heat production of lactating Jersey cows
HP = heat production; GE = gross energy; NEM = net energy for maintenance [0.100 × metabolic BW (Moraes et al., 2015)]; DE = digestible energy; dCP = apparent digested CP; DE CP = digested energy from CP; dNDF = digested NDF; DE NDF = digested energy from NDF; dSTA = digested starch; DE STA = digested energy from starch.
)]; DE = digestible energy; dCP = apparent digested CP; DE CP = digested energy from CP; dNDF = digested NDF; DE NDF = digested energy from NDF; dSTA = digested starch; DE STA = digested energy from starch.
. Cows were housed in a temperature-controlled room (20°C) in tiestalls equipped with rubber mats. In all studies except for 1, cows were adapted to experimental diets for 21 or 28 d. In
, a 10-d adaptation period was used; however, in this study, experimental diets remained the same throughout the experiment, and only feeding frequency (once or twice per day) changed. In all experiments except for 1, a single group of cows were used. In
, cows were divided into 2 groups (termed block in the statistical model). All studies used a 23-h gas collection period that occurred after adaptation. All gas data were adjusted to a 24-h period. Gas data were collected for 1 or 2 d. For the data set where gas was collected for 2 d (n = 155), daily variation in HP was low (CV = 4.6%). Therefore, given our main objectives and the conditions of the current experiment, in which cattle were housed a climate-controlled tiestall facility, we believed that 1 d of gas collection was adequate to estimate HP for this study. Gas collection for 1 d may not be adequate when environmental conditions can vary or feed delivery is not held constant. Cows were fed within the headboxes and allowed free access to water via a water bowl. Details for the specific methods used to quantify gas consumption and production can be found in individual publications (
Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
). Heat production was estimated through calculation of observed O2 consumption and CO2 production with correction for CH4 production and urinary N excretion according to
HP (Mcal) = 0.003866 × O2 (L) + 0.001200 × CO2 (L) − 0.000518 × CH4 (L) − 0.001431 × urinary N (g).
[1]
In all 9 experiments, milk production and composition were quantified for 4 consecutive days. Cows were milked twice per day and fed once per day in all experiments except in
Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
, where cows were fed once or twice per day. In all studies, fecal and urine output was quantified by total collection for 4 consecutive days. Chemical composition of feeds, refusals, and feces was determined by Cumberland Valley Analytical Services Inc. (Waynesboro, PA) except for
, where analyses were completed at the Ruminant Nutrition Laboratory of the University of Nebraska–Lincoln. Methods were similar across location and are described in individual publications (
Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
). Crossover experimental designs were used in all 9 experiments.
Model Derivation Procedure
Response variable was HP (Mcal/d). Models were fit in a hierarchical structure (Figure 1). The first level included individual models using the variable categories of DMI, milk, and digestibility. In the second level, milk or digestibility variables were combined with DMI, where milk yield was used as the milk variable and DE was used as the digestibility variable, and each variable was expressed as a function of DMI to minimize collinearity. In the third level, the constituents of milk yield (fat, protein, and lactose) and DE (digested CP, NDF, and starch) were used with all variables expressed as a function of DMI. To calculate DE from digested nutrients, we multiplied digested CP, NDF, and starch (kg) by their corresponding enthalpies of 5.65, 4.2, and 4.23 Mcal/kg, respectively (
). This resulted in 7 models for each response variable. Specifically, models are as follows, with reference names for each model listed parenthetically in bold: (1) DMI (INT); (2) milk fat, protein, and lactose yield (MILKCOMP); (3) DMI and milk yield (INT+MY); (4) DMI and MILKCOMP variables as a function of DMI (INT+MILKCOMP); (5) daily intake of digested NDF (dNDF), CP (dCP), and starch (dSTA; DIG); (6) DMI and DE (INT+DE); and (7) DMI and NDF, CP, and starch digestibility (INT+DIG). An eighth model was fit to determine HP associated with urinary N excretion that accounted for DMI and milk fat, protein, and lactose yields as a function of DMI (INT+MILKCOMP+UN). The MBW was included in all models to estimate HP. As discussed previously, although forage inclusion (specifically forage NDF) can affect HP (
Effects of diet forage-to-concentrate ratio and intake on energy metabolism in growing beef heifers: Whole body energy and nitrogen balance and visceral heat production.
), it was not evaluated in the current experiment because in most studies the forage NDF was similar across treatments; consequently, differences between treatments was small (mean ± SD; 22.4 ± 2.28%). Initially, models with DMI included BCS and DIM; however, neither BCS nor DIM were significant, thus neither were included in derivation of final models. Fatty acid digestibility was only measured in 3 studies (n = 87), and its effects were not further explored.
Figure 1Illustration of 7 models that were fit for heat production. Models for heat production included metabolic BW. dCP = digested CP; dNDF = digested NDF; dSTA = digested starch; DE = digestible energy.
). To account for the variation associated with using individual animal observations, all models included the random effects of study (n = 9), cow (n = 55), and period within block within study (n = 36). If parameter estimates were close to 0 and nonsignificant (P > 0.15), the corresponding explanatory variables were removed from the model.
Model Evaluation
Simple correlation coefficients among all variables were assessed (Table 2). These correlation coefficients were used to assess potential multicollinearity among explanatory variables and used to aid in the explanation of observed parameters. In final models, multicollinearity among variables was assessed by variance inflation factor (
), which was less than 5 for all variables, and in most cases less than 2. This indicated that the inflation in variance due to multicollinearity was minimal in the models presented.
Table 2Simple correlation coefficients among heat production, animal description, milk yield and components, and digestibility variables included in model development
) were calculated from predicted values with (CCC) or without (uCCC) adjustment for study, cow, and period within block within study. Concordance correlation coefficients assess both accuracy and precision of models. Because the CCC value includes the random intercept terms for each random effect, they typically produce more favorable values compared with uCCC (
Physically adjusted neutral detergent fiber system for lactating dairy cow rations. I: Deriving equations that identify factors that influence effectiveness of fiber.
). When mixed-effects models are used in an uninformed setting, the random effects will not be known; therefore, removal of the random effects, although not a statistically accurate depiction of model fit, is conducted to provide information on model fit when using the model for future use in naïve settings (
Physically adjusted neutral detergent fiber system for lactating dairy cow rations. I: Deriving equations that identify factors that influence effectiveness of fiber.
). Furthermore, for each model, we reported the estimated standard deviation for cow
, study
, and error
, which is more appropriate than the root mean squared error for evaluating mixed effect models (
). The units for all
values were the same as the response variable (Mcal/d).
The objective of the current study was to explain variation in HP rather than to develop a model that would definitively predict HP. We explored generating a prediction error of models using a Monte Carlo cross-validation as described by
Evaluation of the National Research Council (2001) dairy model and derivation of new prediction equations. 2. Rumen degradable and undegradable protein.
; however, the current data set only contained 9 studies that differ greatly by DIM, DMI, milk production, and dietary chemical composition. This resulted in prediction errors, slope, and intercept biases that were affected by which studies were selected by the Monte Carlo simulation. Therefore, results from the cross-validation may be misleading and are not reported herein.
RESULTS AND DISCUSSION
The number of observations for each variable tested were between 261 and 293 (Table 1). Average HP was 22.1 ± 2.86 Mcal/d. The data set included cows with a wide range in DIM (44–410), DMI (9.6–25.0 kg/d), and milk yield (7.8–43.0 kg/d).
Residual versus predicted plots for all HP models are reported in Supplemental Figure S1 (https://doi.org/10.3168/jds.2020-18976). Although a slope bias was observed (P < 0.10) for some models, the magnitude was less than 2% of error variance in all cases, and was thus not a concern.
Factors that Affect HP
Heat arises from the inefficiency in the conversion of a substrate into a product and from the complete oxidation of a substrate (
equation. Negative coefficients are included in the Brouwer equation for urinary N excretion because protein oxidation leads to urinary urea excretion, and thus protein oxidation is not complete. Therefore, O2 and CO2 associated with this incomplete oxidation must be accounted for when estimating HP (
). Total HP is an aggregate of heat for basal metabolism, activity, digestion and absorption, fermentation, product formation, thermoregulation, and waste formation and excretion (
). We hypothesized that MBW and DMI could predict HP. Furthermore, we predicted that the addition of milk component yield, digestibility of nutrients, and urinary N excretion would explain more variation in HP by decreasing
and increasing CCC. Specifically, we expected the following: (1) MBW to account for the NEM component of HP; (2) DMI to indirectly account for the digestion and absorption, fermentation, and product formation component of HP; (3) milk fat, protein, and lactose yield to account for heat of production formation; (4) digested NDF, CP, and starch to account for heat of digestion and absorption, and fermentation; and (5) urinary N excretion to account for heat associated the with waste formation and excretion, and protein catabolism.
In agreement with our hypothesis, INT explained most of the variation in HP (uCCC = 0.84, and
= 1.38; Table 3). However, counter to our hypothesis, the addition of nutrient digestibility and milk component yield did not improve our ability to explain variation in HP based on similar uCCC and
. A lack of improvement in model performance with inclusion of nutrient digestibility and milk component yield may have occurred because most of the variation in HP in the current study was explained by MBW and DMI. Furthermore, no difference in model performance would suggest that similar variation was explained by milk component yield or digested nutrients compared with DMI, which is not surprising given the correlation between most of these variables (r > 0.55; P < 0.05; Table 2). The relationship between HP and other response variables may improve our understanding of biological sources of variation in HP.
Table 3Parameter estimates and fit statistics for models to predict heat production (Mcal/d) of Jersey cows from metabolic BW (MBW) plus DMI (INT); milk fat, protein, and lactose yield (MILKCOMP); INT plus milk yield (INT+MY); INT plus MILKCOMP variables as a function of DMI (INT+MILKCOMP); digested NDF, CP, and starch (DIG); INT plus digested energy (INT+DE); or INT plus NDF, CP, and starch digestibility (INT+DIG)
n = number of observations; CCC = concordance correlation coefficient; BIC = Bayesian information criterion; uCCC = concordance correlation coefficient without random effects adjustment;
σ^c = square root of the estimated variance associated with cow (Mcal);
σ^s = square root of the estimated variance associated with study (Mcal);
σ^e = square root of the estimated variance associated with error (Mcal).
2 n = number of observations; CCC = concordance correlation coefficient; BIC = Bayesian information criterion; uCCC = concordance correlation coefficient without random effects adjustment;
= square root of the estimated variance associated with cow (Mcal);
= square root of the estimated variance associated with study (Mcal);
= square root of the estimated variance associated with error (Mcal).
Given the other variables included in the models, increasing MBW increased HP on average by 0.111 to 0.165 Mcal per unit of MBW (Table 3). Because basal metabolism is a function of MBW (
), MBW was expected to account for a large fraction of the variation in HP. The MILKCOMP and DIG models, which did not include DMI, resulted in MBW coefficients of 0.165 and 0.147. These values are likely larger than actual NEM, and thus demonstrate the importance of accounting for DMI when explaining variation in HP. The MILKCOMP and DIG models aside, our average MBW coefficient was 0.122 Mcal, which is 53% greater than the current NEM value used by
recently reported that NEM has increased over time from 0.074 to 0.124 for measurements taken during 1963 to 1973, and 1984 to 1995, respectively. From 935 observations from European calorimetry studies published between 1992 to 2010, NEM averaged 0.105 Mcal/MBW (
). Most of the historical energetics data were collected on Holsteins cows; however, energy utilization and maintenance energy requirements did not differ between Holsteins and Norwegian, Jersey × Holstein, or Norwegian × Holstein cows (
On average, increasing DMI increased HP by 0.560 to 0.600 Mcal/kg. Increased DMI typically corresponded with increased HP, which has been previously observed (
Effects of diet forage-to-concentrate ratio and intake on energy metabolism in growing beef heifers: Whole body energy and nitrogen balance and visceral heat production.
Effects of diet forage:concentrate ratio and metabolizable energy intake on visceral organ growth and in vitro oxidative capacity of gut tissues in sheep.
Effects of diet forage-to-concentrate ratio and intake on energy metabolism in growing beef heifers: Whole body energy and nitrogen balance and visceral heat production.
). Additionally, DMI was correlated with milk production (r = 0.63; Table 2); therefore, a strong relationship between DMI and HP was not surprising.
Synthesis of milk constituents is a major contributor to the HI of lactating cows. Although empirical models do not directly measure the underlying mechanistic nature of biology (
), interpretation of parameter estimates may be useful in determining the relationship between response and explanatory variables. When accounting for MBW and DMI, HP increased on average by 1.54 ± 0.57 Mcal per unit increase in milk yield over DMI (kg/kg of DMI); therefore, we were interested in partitioning the HP associated with MY into HP due to fat, protein, or lactose synthesis. For the MILKCOMP and INT+MILKCOMP models to estimate HP, the parameter estimate for milk lactose yield was not different from 0 (P > 0.27), and was consequently removed from the final model. Although synthesis of milk lactose results in HP, milk lactose was highly correlated with milk fat (r = 0.75) and milk protein (r = 0.83; Table 2); therefore, changes in milk fat and protein yield likely explained the variation in HP that arose from changes in milk lactose yield. The parameter estimate for milk protein was 2-fold greater than that of milk fat when expressed on a kg/d (6.17 ± 1.13 vs. 1.83 ± 0.656) or kg/kg of DMI basis (40.2 ± 23.0 vs. 16.4 ± 12.0; Table 3). This difference may be explained by a greater energetic efficiency for synthesizing milk fat compared with milk protein. De novo synthesis of milk fat is about 0.70 efficient (
), and thus the theoretical maximum efficiency of milk protein synthesis is 0.74. However, because of protein turnover, actual protein synthesis in mammary glands will always exceed milk protein yield. Estimated rates of protein turnover range from 140% of milk protein yield in lactating dairy cows (
), with the prediction depending on which labeled AA was used for quantification. When accounting for the fact that actual protein synthesis is 40 to 200% greater than milk protein yield, energetic efficiency of milk protein synthesis ranges from 0.48 to 0.67. Consequently, milk fat production will lead to a lower HI compared with milk protein synthesis.
Increasing the supply of digested nutrients was expected to increase HP because of the underlying correlation with DMI. However, digestion of NDF and starch gives rise to different VFA profiles, HI associated with intermediary metabolism differs by VFA (
in: Sejrsen K. Hvelplund T. Nielsen M.O. Ruminant Physiology: Digestion, Metabolism and Impact of Nutrition on Gene Expression, Immunology and Stress. Wageningen Academic Publishers,
Wageningen, the Netherlands2006: 225-248
). Therefore, we expected to see differences in HP associated with digested NDF, CP, and starch. The parameter estimates for dNDF (kg/d) in the DIG HP model were not different from 0 (0.146 ± 0.242 Mcal/kg; P = 0.55), and thus were removed from the final model. Because heat of fermentation is a component of HP (
) and NDF is primarily digested via ruminal fermentation, we were surprised that the effect of dNDF on HP was small. However, dNDF was correlated with dCP (r = 0.49) and dSTA (r = −0.19; Table 2); therefore, variation in HP associated with dNDF was likely explained by dCP and dSTA. Furthermore, in the data set used in the current study, change in dietary NDF content was primarily driven by changes in byproduct inclusion (r = 0.83 for dietary NDF and nonforage NDF). It is likely that nonforage NDF does not stimulate a similar increase in HP as does increasing forage inclusion (
Effects of diet forage-to-concentrate ratio and intake on energy metabolism in growing beef heifers: Whole body energy and nitrogen balance and visceral heat production.
). Dietary NDF content, rather than digestible NDF, may be a better variable for explaining variation in HP. Therefore, we tested the effects of dietary NDF content in place of NDF digestibility. When dietary NDF content replaced dNDF in the DIG model, HP increased by 0.221 ± 0.0587 Mcal per percentage unit increase in dietary NDF on average (data not shown). Additionally, for the DIG HP model, both dCP and dSTA were associated with HP when MBW was included (Table 3). For the DIG HP model, the coefficient for dCP was almost 2-fold that of dSTA (1.34 ± 0.372 vs. 0.762 ± 0.167 Mcal/kg; Table 3). An increase in HP with increased dCP was expected.
reported that feeding CP in excess of requirements resulted in decreased NEL, which may be attributable to an increase in HP associated with metabolizing excess CP and transforming DE into urinary energy (
). The HP associated with excess N may originate from ureagenesis or AA catabolism (see discussion on the effects of urinary N excretion on HP and HI). Most dietary starch is digested to VFA in the rumen and will subsequently be used as an energy source for the body or be converted into products. Our results suggest that a greater quantity of HP is associated with each kilogram of dCP compared with dSTA, which supports research conducted over 100 yr ago by Rubner, who observed that the HP associated with oxidation of protein is greater than fat or carbohydrates (
Additional digestibility models were fit that included DE and digestibility of CP, NDF, and starch. We expected INT+DE and INT+DIG to improve model performance compared with INT; however, as previously discussed, these models resulted in similar uCCC and
values. For the INT+DE model to estimate HP, increasing DE (Mcal/kg of DM) decreased HP (−1.72 ± 0.592 Mcal). In the current database, DE content was positively correlated (r = 0.25) with dietary crude fat and negatively correlated (r = −0.29) with dietary NDF content (data not shown). As discussed previously, increased dietary fat and decreased NDF should result in a diet that will lead to less HP. For the INT+DIG HP model, decreased HP was observed with increased digestibility of CP and NDF (Mcal/%; −0.0432 ± 0.27 and −0.0666 ± 0.018); as starch digestibility increased, HP increased on average by 0.0910 ± 0.046 Mcal/%. Starch digestibility was positively correlated with milk and lactose yield (Table 2), which will lead to increased HP. Negative coefficients for CP and NDF digestibility may have occurred because CP digestibility were negatively correlated (r < −0.25) with DMI and milk component yields, and NDF digestibility was negatively correlated (r = −0.13) with MBW. As NDF digestibility increased, energy expenditure associated with digestion, and HP from mass of ruminal tissue may have decreased because of changes in the source of NDF (
Effects of diet forage-to-concentrate ratio and intake on energy metabolism in growing beef heifers: Whole body energy and nitrogen balance and visceral heat production.
). Additionally, at a constant DMI, increasing NDF and CP digestibility may increase the partitioning of energy toward milk and away from heat, leading to decreased HP.
The Effects of Urinary N Excretion on HP and Increment
Feeding excess protein to dairy cows increases HP and decreases energy balance (
). Urinary N excretion can serve as an indicator of AA oxidation, and thus excess protein; therefore we hypothesized that the inclusion of urinary N excretion could improve our ability to explain variation in HP and HI. Compared with other models for HP and HI, the models that included urinary N excretion had similar uCCC (0.83) and
values (1.40 Mcal; Table 4). Overall, partitioning of the variation in HP was not improved with the inclusion of urinary N excretion.
Table 4Parameter estimates and fit statistics for models to predict heat production (Mcal/d) of Jersey cows from metabolic BW (MBW; heat production only); DMI; milk fat, protein, and lactose yield; and urinary N excretion
n = number of observations; BIC = Bayesian information criterion; CCC = concordance correlation coefficient; uCCC = concordance correlation coefficient without random effects adjustment;
σ^c = square root of the estimated variance associated with cow (Mcal);
σ^s = square root of the estimated variance associated with study (Mcal);
σ^s = square root of the estimated variance associated with error (Mcal).
n
261
BIC
1,030
CCC
0.89
uCCC
0.84
0.52
1.65
1.39
1 MBW = BW0.75.
2 n = number of observations; BIC = Bayesian information criterion; CCC = concordance correlation coefficient; uCCC = concordance correlation coefficient without random effects adjustment;
= square root of the estimated variance associated with cow (Mcal);
= square root of the estimated variance associated with study (Mcal);
= square root of the estimated variance associated with error (Mcal).
A secondary objective of determining the relationship between HP and urinary N excretion was to quantify the contribution of protein catabolism to heat energy. On average, when accounting for MBW, DMI, milk fat yield, and milk protein yield, HP increased by 5.32 ± 2.43 Mcal/kg of urinary N excretion. In an analysis of a historical indirect calorimetry data set,
determined that HP was associated with excess N intake, which was defined as the digested N supply minus N requirements for maintenance, milk production, and pregnancy as defined by
). Consequently, metabolic fecal N was subtracted out as a component of fecal N excretion and maintenance N requirements, thus underestimating excess N. Both will lead to a greater coefficient than the method used in the current study. In lactating dairy cows, the primary form of urinary N is urea (
Prediction of urinary nitrogen and urinary urea nitrogen excretion by lactating dairy cattle in northwestern Europe and North America: A meta-analysis.
). Ureagenesis occurs as a result of catabolism of AA by rumen microbes, catabolism of AA that are supplied in excess of requirements, and turnover of body protein stores (
in: Sejrsen K. Hvelplund T. Nielsen M.O. Ruminant Physiology: Digestion, Metabolism and Impact of Nutrition on Gene Expression, Immunology and Stress. Wageningen Academic Publishers,
Wageningen, the Netherlands2006: 225-248
in: Sejrsen K. Hvelplund T. Nielsen M.O. Ruminant Physiology: Digestion, Metabolism and Impact of Nutrition on Gene Expression, Immunology and Stress. Wageningen Academic Publishers,
Wageningen, the Netherlands2006: 225-248
suggests that ureagenesis may not be the primary contributor to increased HP with excess protein intake; rather, catabolism of AA leads to increased HP with increase CP intake. As described by
, because ruminants are evolutionarily adapted to be dependent on urea recycling, a large energetic cost to ureagenesis would be constrictive to urea recycling. Supporting evidence includes that increasing ammonia supply or feeding urea does not reduce milk yield (
). The HP that was associated with UN in the current experiment on average was 1.0 Mcal/d (1.3% of GE intake). In the current study, urinary N excretion explained variation in HP that was associated with the catabolism of protein.
Due to the amount of heat associated with urinary N excretion, the question arises: will increasing dietary CP above requirements to replace carbohydrates decrease NEL supply? Assuming a true digestibility of dietary CP of 85% (
) and that all dietary CP is excreted as urinary N, the HP associated with urinary N excretion is 0.72 Mcal/kg of CP. Additionally, N from excess CP is primarily excreted as urea in urinary N (
Prediction of urinary nitrogen and urinary urea nitrogen excretion by lactating dairy cattle in northwestern Europe and North America: A meta-analysis.
), which has an enthalpy of 5.4 Mcal/kg of N or 0.86 Mcal/kg of CP. Thus, the sum of energy associated with increased HP and urinary energy excretion from feeding CP as an energy source is 1.6 Mcal/kg. This is similar to the difference in enthalpy between CP and carbohydrates (5.6 vs. 4.2 Mcal/kg;
); therefore, differences in energy supply between feeding CP or carbohydrates is dependent on digestibility of the CP or carbohydrate source. Average true digestibility of CP is slightly less than starch (85 vs. 92%), but true digestibility of CP is much greater than NDF (
Effect of cereal grain type and corn grain harvesting and processing methods on intake, digestion, and milk production by dairy cows through a meta-analysis.
). Thus, replacing dietary NDF with CP may improve NEL supply.
CONCLUSIONS
From 293 individual animal observations on lactating Jersey cows, HP averaged 22.1 ± 2.86 Mcal/d or 28.1 ± 3.70% of GE intake. Our results suggest that variation in HP is explained by MBW and DMI. Counter to our hypothesis, the inclusion of milk component yield and nutrient digestibility did not improve our ability to explain variation in HP and HI. Parameter estimates from these models agree with other biological and biochemical estimates of energy transaction in lactating cows. For example, HP associated with milk protein synthesis was approximately 2-fold that of milk fat synthesis. Additionally, for every kilogram of urinary N excreted, HP increased on average by 5.32 Mcal.
ACKNOWLEDGMENTS
The authors thank the University of Nebraska Dairy Metabolism (Lincoln, NE) staff and students for care of the experimental animals and assistance with collections. Gratitude is also extended to the American Jersey Cattle Associations AJCC Research Foundation (Reynoldsburg, OH) for funding. The authors have not stated any conflicts of interest.
Effect of cereal grain type and corn grain harvesting and processing methods on intake, digestion, and milk production by dairy cows through a meta-analysis.
Effects of diet forage:concentrate ratio and metabolizable energy intake on visceral organ growth and in vitro oxidative capacity of gut tissues in sheep.
Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
in: Sejrsen K. Hvelplund T. Nielsen M.O. Ruminant Physiology: Digestion, Metabolism and Impact of Nutrition on Gene Expression, Immunology and Stress. Wageningen Academic Publishers,
Wageningen, the Netherlands2006: 225-248
Effects of diet forage-to-concentrate ratio and intake on energy metabolism in growing beef heifers: Whole body energy and nitrogen balance and visceral heat production.
Prediction of urinary nitrogen and urinary urea nitrogen excretion by lactating dairy cattle in northwestern Europe and North America: A meta-analysis.
Physically adjusted neutral detergent fiber system for lactating dairy cow rations. I: Deriving equations that identify factors that influence effectiveness of fiber.
Evaluation of the National Research Council (2001) dairy model and derivation of new prediction equations. 2. Rumen degradable and undegradable protein.