Advertisement
Research| Volume 104, ISSUE 1, P346-356, January 2021

Factors that affect heat production in lactating Jersey cows

Open ArchivePublished:November 11, 2020DOI:https://doi.org/10.3168/jds.2020-18976

      ABSTRACT

      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.

      Key words

      INTRODUCTION

      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 (
      • Kleiber M.
      The Fire of Life: An Introduction to Animal Energetics.
      ). In lactating Jersey cows, heat production (HP) accounts for 25 to 32% of the gross energy (GE) intake (
      • Drehmel O.R.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Fernando S.C.
      • Miller P.S.
      • Hales K.E.
      • Kononoff P.J.
      The influence of fat and hemicellulose on methane production and energy utilization in lactating Jersey cattle.
      ;
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Energy balance and diurnal variation in methane production as affected by feeding frequency in Jersey cows in late lactation.
      ;
      • Reynolds M.A.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Herrick K.J.
      • Hales K.E.
      • Watson A.K.
      • Kononoff P.J.
      Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming common coproducts.
      ). 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 (
      • NRC
      Nutritional Energetics of Domestic Animals and Glossary of Energy Terms..
      ). 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 (
      • Belyea R.L.
      • Adams M.W.
      Energy and nitrogen utilization of high versus low producing dairy cows.
      ), forage inclusion (
      • Reynolds C.K.
      • Tyrrell H.F.
      • Reynolds P.J.
      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.
      ), and N intake above requirements (
      • Tyrrell H.
      • Moe P.W.
      • Flatt W.P.
      Influence of excess protein intake on energy metabolism of the dairy cow.
      ;
      • Reed K.F.
      • Bonfa H.C.
      • Dijkstra J.
      • Casper D.P.
      • Kebreab E.
      Estimating the energetic cost of feeding excess dietary nitrogen to dairy cows.
      ). Various factors are likely to influence HP and HI, and quantification of these may improve our ability to estimate NEL of diets.
      We recently observed that the
      • NRC
      Nutrient Requirements of Dairy Cattle.
      model underestimated dietary NEL content compared with measured NEL, when NEM was assumed to equal 0.080 Mcal per metabolic BW (BW0.75, MBW;
      • Morris D.L.
      • Brown-Brandl T.M.
      • Hales K.E.
      • Harvatine K.J.
      • Kononoff P.J.
      Effects of high-starch or high-fat diets formulated to be isoenergetic on energy and nitrogen partitioning and utilization in lactating Jersey cows.
      ). Furthermore, it has been shown that the efficiency of converting ME into NEL (denoted as kL) is affected by diet (
      • Moraes L.E.
      • Kebreab E.
      • Strathe A.B.
      • Dijkstra J.
      • France J.
      • Casper D.P.
      • Fadel J.G.
      Multivariate and univariate analysis of energy balance data from lactating dairy cows.
      ;
      • Morris D.L.
      • Judy J.V.
      • Kononoff P.J.
      Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
      ). Historically, the NEL system assumed a constant or near constant kL (
      • Moe P.W.
      • Tyrrell H.F.
      • Flatt W.P.
      Net energy value of feeds for lactation.
      ;
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ). 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 (
      • Moraes L.E.
      • Kebreab E.
      • Strathe A.B.
      • Dijkstra J.
      • France J.
      • Casper D.P.
      • Fadel J.G.
      Multivariate and univariate analysis of energy balance data from lactating dairy cows.
      ). Biologically, kL is likely dependent upon the substrates utilized and products formed (
      • Baldwin R.L.
      Modeling Ruminant Digestion and Metabolism..
      ). For example, HI has been reported to be greater for acetate than propionate or butyrate when these substrates were used for maintenance (
      • Holter J.B.
      • Heald C.W.
      • Colovos N.F.
      Heat increments of steam-volatile fatty acids infused separately and in a mixture into fasting cows.
      ). 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 (
      • Foth A.J.
      • Brown-Brandl T.
      • Hanford K.J.
      • Miller P.S.
      • Garcia Gomez G.
      • Kononoff P.J.
      Energy content of reduced-fat dried distillers grains with solubles for lactating dairy cows.
      ;
      • Drehmel O.R.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Fernando S.C.
      • Miller P.S.
      • Hales K.E.
      • Kononoff P.J.
      The influence of fat and hemicellulose on methane production and energy utilization in lactating Jersey cattle.
      ;
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Energy balance and diurnal variation in methane production as affected by feeding frequency in Jersey cows in late lactation.
      ,
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Harvatine K.J.
      • Miller P.S.
      • Kononoff P.J.
      Increasing the concentration of linolenic acid in diets fed to Jersey cows in late lactation does not affect methane production.
      ,
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Reducing methane production with corn oil and calcium sulfate: Responses on whole-animal energy and nitrogen balance in dairy cattle.
      ;
      • Knoell A.L.
      • Judy J.V.
      • Wilson H.C.
      • Herrick K.J.
      • Fernando S.C.
      • Kononoff P.J.
      Energy utilization in lactating Jersey cows consuming a mixture of distillers dried grains with solubles and straw in replacement of alfalfa hay.
      ;
      • Reynolds M.A.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Herrick K.J.
      • Hales K.E.
      • Watson A.K.
      • Kononoff P.J.
      Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming common coproducts.
      ;
      • Morris D.L.
      • Brown-Brandl T.M.
      • Hales K.E.
      • Harvatine K.J.
      • Kononoff P.J.
      Effects of high-starch or high-fat diets formulated to be isoenergetic on energy and nitrogen partitioning and utilization in lactating Jersey cows.
      ,
      • Morris D.L.
      • Judy J.V.
      • Kononoff P.J.
      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 (
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Harvatine K.J.
      • Miller P.S.
      • Kononoff P.J.
      Increasing the concentration of linolenic acid in diets fed to Jersey cows in late lactation does not affect methane production.
      ,
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Reducing methane production with corn oil and calcium sulfate: Responses on whole-animal energy and nitrogen balance in dairy cattle.
      ), a covariate period was conducted where DE, ME, and nutrient digestibility were not measured. Additionally, during the experimental period of
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Harvatine K.J.
      • Miller P.S.
      • Kononoff P.J.
      Increasing the concentration of linolenic acid in diets fed to Jersey cows in late lactation does not affect methane production.
      , 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
      Item
      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.
      nMeanSDMinimumMedianMaximum
      GE intake, Mcal/d29379.811.141.679.8106.4
      HP, Mcal/d29322.12.8615.322.129.2
      HP, % of GE intake29328.13.7018.827.642.2
      Animal descriptions
       DIM, d29318272.344174410
       Parity2933.00.952.03.05.0
       BW, kg29345344342452568
       Metabolic BW, kg0.7529398.27.279.598.0116.4
       BCS
      On a 1 to 5 scale.
      2723.270.352.003.254.25
       DMI, kg/d29318.22.539.618.225.0
      Milk and components yield
       Milk yield, kg/d29325.05.497.824.643.0
       Milk yield, kg/kg of DMI2931.380.2420.721.372.15
       Fat, kg/d2931.400.2920.481.392.19
       Fat, kg/kg of DMI2930.07680.01270.04850.07660.116
       Protein, kg/d2930.900.1650.330.891.29
       Protein, kg/kg of DMI2930.04930.00580.03320.04900.0701
       Lactose, kg/d2901.200.2850.331.192.10
       Lactose, kg/kg of DMI2900.06580.01300.03390.06560.110
      Digested energy and nutrients
       DE, Mcal/kg of DM2612.880.2192.302.893.77
       dCP, kg/d2772.250.4041.322.213.39
       CP digestibility, %27768.26.349.468.382.8
       dNDF, kg/d2772.880.7220.812.835.06
       NDF digestibility, %27747.36.624.546.965.8
       dSTA, kg/d2774.021.101.524.196.62
       Starch digestibility, %27795.02.793.595.2100.0
       Urinary N excretion, g/d26118857.456.9182394
      Dietary composition, % of DM
       CP28517.70.9515.217.719.5
       NDF28531.93.5125.532.043.0
       Starch28525.03.4816.225.831.1
       Crude fat2814.20.742.24.26.4
      1 HP = heat production; GE = gross energy; NEM = net energy for maintenance [0.100 × metabolic BW (
      • Moraes L.E.
      • Kebreab E.
      • Strathe A.B.
      • Dijkstra J.
      • France J.
      • Casper D.P.
      • Fadel J.G.
      Multivariate and univariate analysis of energy balance data from lactating dairy cows.
      )]; 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.
      2 On a 1 to 5 scale.
      In all 9 experiments, O2 consumption and CO2 and CH4 production were determined using headbox-type indirect calorimeters as described by
      • Foth A.J.
      • Brown-Brandl T.
      • Hanford K.J.
      • Miller P.S.
      • Garcia Gomez G.
      • Kononoff P.J.
      Energy content of reduced-fat dried distillers grains with solubles for lactating dairy cows.
      . 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
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Energy balance and diurnal variation in methane production as affected by feeding frequency in Jersey cows in late lactation.
      , 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
      • Foth A.J.
      • Brown-Brandl T.
      • Hanford K.J.
      • Miller P.S.
      • Garcia Gomez G.
      • Kononoff P.J.
      Energy content of reduced-fat dried distillers grains with solubles for lactating dairy cows.
      , 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 (
      • Foth A.J.
      • Brown-Brandl T.
      • Hanford K.J.
      • Miller P.S.
      • Garcia Gomez G.
      • Kononoff P.J.
      Energy content of reduced-fat dried distillers grains with solubles for lactating dairy cows.
      ;
      • Drehmel O.R.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Fernando S.C.
      • Miller P.S.
      • Hales K.E.
      • Kononoff P.J.
      The influence of fat and hemicellulose on methane production and energy utilization in lactating Jersey cattle.
      ;
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Energy balance and diurnal variation in methane production as affected by feeding frequency in Jersey cows in late lactation.
      ,
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Harvatine K.J.
      • Miller P.S.
      • Kononoff P.J.
      Increasing the concentration of linolenic acid in diets fed to Jersey cows in late lactation does not affect methane production.
      ,
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Reducing methane production with corn oil and calcium sulfate: Responses on whole-animal energy and nitrogen balance in dairy cattle.
      ;
      • Knoell A.L.
      • Judy J.V.
      • Wilson H.C.
      • Herrick K.J.
      • Fernando S.C.
      • Kononoff P.J.
      Energy utilization in lactating Jersey cows consuming a mixture of distillers dried grains with solubles and straw in replacement of alfalfa hay.
      ;
      • Reynolds M.A.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Herrick K.J.
      • Hales K.E.
      • Watson A.K.
      • Kononoff P.J.
      Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming common coproducts.
      ;
      • Morris D.L.
      • Brown-Brandl T.M.
      • Hales K.E.
      • Harvatine K.J.
      • Kononoff P.J.
      Effects of high-starch or high-fat diets formulated to be isoenergetic on energy and nitrogen partitioning and utilization in lactating Jersey cows.
      ,
      • Morris D.L.
      • Judy J.V.
      • Kononoff P.J.
      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
      • Brouwer E.
      Report of sub-committee on constants and factors.
      as follows:
      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
      • Morris D.L.
      • Judy J.V.
      • Kononoff P.J.
      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 milked 3 times per day, and
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Energy balance and diurnal variation in methane production as affected by feeding frequency in Jersey cows in late lactation.
      , 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
      • Foth A.J.
      • Brown-Brandl T.
      • Hanford K.J.
      • Miller P.S.
      • Garcia Gomez G.
      • Kononoff P.J.
      Energy content of reduced-fat dried distillers grains with solubles for lactating dairy cows.
      and
      • Reynolds M.A.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Herrick K.J.
      • Hales K.E.
      • Watson A.K.
      • Kononoff P.J.
      Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming common coproducts.
      , 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 (
      • Foth A.J.
      • Brown-Brandl T.
      • Hanford K.J.
      • Miller P.S.
      • Garcia Gomez G.
      • Kononoff P.J.
      Energy content of reduced-fat dried distillers grains with solubles for lactating dairy cows.
      ;
      • Drehmel O.R.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Fernando S.C.
      • Miller P.S.
      • Hales K.E.
      • Kononoff P.J.
      The influence of fat and hemicellulose on methane production and energy utilization in lactating Jersey cattle.
      ;
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Energy balance and diurnal variation in methane production as affected by feeding frequency in Jersey cows in late lactation.
      ,
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Harvatine K.J.
      • Miller P.S.
      • Kononoff P.J.
      Increasing the concentration of linolenic acid in diets fed to Jersey cows in late lactation does not affect methane production.
      ,
      • Judy J.V.
      • Bachman G.C.
      • Brown-Brandl T.M.
      • Fernando S.C.
      • Hales K.E.
      • Miller P.S.
      • Stowell R.R.
      • Kononoff P.J.
      Reducing methane production with corn oil and calcium sulfate: Responses on whole-animal energy and nitrogen balance in dairy cattle.
      ;
      • Knoell A.L.
      • Judy J.V.
      • Wilson H.C.
      • Herrick K.J.
      • Fernando S.C.
      • Kononoff P.J.
      Energy utilization in lactating Jersey cows consuming a mixture of distillers dried grains with solubles and straw in replacement of alfalfa hay.
      ;
      • Reynolds M.A.
      • Brown-Brandl T.M.
      • Judy J.V.
      • Herrick K.J.
      • Hales K.E.
      • Watson A.K.
      • Kononoff P.J.
      Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming common coproducts.
      ;
      • Morris D.L.
      • Brown-Brandl T.M.
      • Hales K.E.
      • Harvatine K.J.
      • Kononoff P.J.
      Effects of high-starch or high-fat diets formulated to be isoenergetic on energy and nitrogen partitioning and utilization in lactating Jersey cows.
      ,
      • Morris D.L.
      • Judy J.V.
      • Kononoff P.J.
      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 (
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ). 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 (
      • Reynolds C.K.
      • Tyrrell H.F.
      • Reynolds P.J.
      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 thumbnail gr1
      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.
      Models were fit in R (version 3.5.2) with the lmer function (
      • Kuznetsova A.
      • Brockhoff P.B.
      • Christensen R.H.B.
      lmerTest package: Tests in linear mixed effects models: Version 3.0-1.
      ). 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 (
      • Roman-Garcia Y.
      • White R.R.
      • Firkins J.L.
      Meta-analysis of postruminal microbial nitrogen flows in dairy cattle. I. Derivation of equations.
      ), 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
      Items are repeated horizontally and vertically in the same order and the same units. HP = heat production (Mcal); MBW = metabolic BW (kg0.75); DMI (kg/d); MY = milk yield (kg/d); MY/DMI = milk yield/DMI (kg/kg); Fat = milk fat yield (kg/d); fat/DMI = milk fat yield/DMI (kg/kg); pro = milk protein yield (kg/d); pro/DMI = milk protein yield/DMI (kg/kg); lact = milk lactose yield (kg/d); lact/DMI = milk lactose yield/DMI (kg/kg); DE/DM = digested energy (Mcal/kg of DM); dCP = apparent digested CP (kg/d); CP digestibility (%); dNDF = digested NDF (kg/d), NDF digestibility (%); dSTA = digested starch (kg/d); starch digestibility (%); UN = urinary N excretion (g/d).
      ItemVariable
      MBWDMIMYMY/ DMIFatFat/ DMIProPro/ DMILactLact/ DMIDE/ DMdCPCP digestibilitydNDFNDF digestibilitydSTAStarch digestibilityUN
      HP0.26
      P < 0.05.
      0.68
      P < 0.05.
      0.50
      P < 0.05.
      0.090.40
      P < 0.05.
      −0.080.59
      P < 0.05.
      0.12
      P < 0.05.
      0.49
      P < 0.05.
      0.10−0.14
      P < 0.05.
      0.43
      P < 0.05.
      −0.21
      P < 0.05.
      0.38
      P < 0.05.
      0.070.32
      P < 0.05.
      0.15
      P < 0.05.
      0.01
      MBW0.20
      P < 0.05.
      −0.05−0.26
      P < 0.05.
      0.03−0.17
      P < 0.05.
      0.07−0.16
      P < 0.05.
      −0.09−0.28
      P < 0.05.
      0.060.13
      P < 0.05.
      0.04−0.03−0.13
      P < 0.05.
      0.30
      P < 0.05.
      0.000.16
      P < 0.05.
      DMI0.63
      P < 0.05.
      0.010.63
      P < 0.05.
      −0.030.78
      P < 0.05.
      0.080.60
      P < 0.05.
      0.03−0.050.58
      P < 0.05.
      −0.30
      P < 0.05.
      0.55
      P < 0.05.
      0.090.54
      P < 0.05.
      −0.02−0.02
      MY0.77
      P < 0.05.
      0.78
      P < 0.05.
      0.48
      P < 0.05.
      0.86
      P < 0.05.
      0.64
      P < 0.05.
      0.99
      P < 0.05.
      0.78
      P < 0.05.
      −0.050.25
      P < 0.05.
      −0.39
      P < 0.05.
      0.21
      P < 0.05.
      −0.010.47
      P < 0.05.
      0.12
      P < 0.05.
      −0.18
      P < 0.05.
      MY/DMI0.50
      P < 0.05.
      0.66
      P < 0.05.
      0.48
      P < 0.05.
      0.78
      P < 0.05.
      0.79
      P < 0.05.
      0.99
      P < 0.05.
      −0.03−0.13
      P < 0.05.
      −0.25
      P < 0.05.
      −0.16
      P < 0.05.
      −0.090.16
      P < 0.05.
      0.17
      P < 0.05.
      −0.18
      P < 0.05.
      Fat0.76
      P < 0.05.
      0.77
      P < 0.05.
      0.50
      P < 0.05.
      0.75
      P < 0.05.
      0.49
      P < 0.05.
      −0.030.22
      P < 0.05.
      −0.40
      P < 0.05.
      0.16
      P < 0.05.
      −0.110.53
      P < 0.05.
      0.01−0.03
      Fat/DMI0.34
      P < 0.05.
      0.60
      P < 0.05.
      0.47
      P < 0.05.
      0.63
      P < 0.05.
      0.01−0.18
      P < 0.05.
      −0.25
      P < 0.05.
      −0.24
      P < 0.05.
      −0.21
      P < 0.05.
      0.21
      P < 0.05.
      0.020.00
      Pro0.67
      P < 0.05.
      0.83
      P < 0.05.
      0.47
      P < 0.05.
      −0.060.40
      P < 0.05.
      −0.33
      P < 0.05.
      0.28
      P < 0.05.
      −0.010.57
      P < 0.05.
      0.03−0.09
      Pro/DMI0.63
      P < 0.05.
      0.75
      P < 0.05.
      −0.04−0.04−0.17
      P < 0.05.
      −0.18
      P < 0.05.
      −0.13
      P < 0.05.
      0.26
      P < 0.05.
      0.07−0.10
      Lact0.81
      P < 0.05.
      −0.050.25
      P < 0.05.
      −0.36
      P < 0.05.
      0.22
      P < 0.05.
      0.020.42
      P < 0.05.
      0.13
      P < 0.05.
      −0.18
      P < 0.05.
      Lact/DMI−0.03−0.08−0.22
      P < 0.05.
      −0.10−0.040.120.16
      P < 0.05.
      −0.18
      P < 0.05.
      DE/DM0.37
      P < 0.05.
      0.55
      P < 0.05.
      0.18
      P < 0.05.
      0.34
      P < 0.05.
      −0.120.080.08
      dCP0.54
      P < 0.05.
      0.49
      P < 0.05.
      0.27
      P < 0.05.
      0.06−0.040.41
      P < 0.05.
      CP digestibility0.040.30
      P < 0.05.
      −0.44
      P < 0.05.
      −0.020.50
      P < 0.05.
      dNDF0.73
      P < 0.05.
      −0.19
      P < 0.05.
      −0.07−0.06
      NDF digestibility−0.55
      P < 0.05.
      0.03−0.08
      dSTA0.14
      P < 0.05.
      −0.12
      P < 0.05.
      Starch digestibility−0.14
      P < 0.05.
      1 Items are repeated horizontally and vertically in the same order and the same units. HP = heat production (Mcal); MBW = metabolic BW (kg0.75); DMI (kg/d); MY = milk yield (kg/d); MY/DMI = milk yield/DMI (kg/kg); Fat = milk fat yield (kg/d); fat/DMI = milk fat yield/DMI (kg/kg); pro = milk protein yield (kg/d); pro/DMI = milk protein yield/DMI (kg/kg); lact = milk lactose yield (kg/d); lact/DMI = milk lactose yield/DMI (kg/kg); DE/DM = digested energy (Mcal/kg of DM); dCP = apparent digested CP (kg/d); CP digestibility (%); dNDF = digested NDF (kg/d), NDF digestibility (%); dSTA = digested starch (kg/d); starch digestibility (%); UN = urinary N excretion (g/d).
      * P < 0.05.
      To compare among models with the same response variables, concordance correlation coefficients (CCC;
      • Lin L.I.
      A concordance correlation coefficient to evaluate reproducibility.
      ) 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 (
      • White R.R.
      • Hall M.B.
      • Firkins J.L.
      • Kononoff P.J.
      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 (
      • White R.R.
      • Hall M.B.
      • Firkins J.L.
      • Kononoff P.J.
      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 (σc), study (σs), and error (σe), which is more appropriate than the root mean squared error for evaluating mixed effect models (
      • Boerman J.P.
      • Firkins J.L.
      • St-Pierre N.R.
      • Lock A.L.
      Intestinal digestibility of long-chain fatty acids in lactating dairy cows: A meta-analysis and meta-regression.
      ). 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
      • White R.R.
      • Roman-Garcia Y.
      • Firkins J.L.
      • Kononoff P.
      • VandeHaar M.J.
      • Tran H.
      • McGill T.
      • Garnett R.
      • Hanigan M.D.
      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 (
      • Baldwin R.L.
      Modeling Ruminant Digestion and Metabolism..
      ); it is commonly estimated in ruminants using the
      • Brouwer E.
      Report of sub-committee on constants and factors.
      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 (
      • McLean J.
      • Tobin G.
      Animal and Human Calorimetry..
      ). Total HP is an aggregate of heat for basal metabolism, activity, digestion and absorption, fermentation, product formation, thermoregulation, and waste formation and excretion (
      • NRC
      Nutritional Energetics of Domestic Animals and Glossary of Energy Terms..
      ). 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 σe 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 σe = 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 σe. 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)
      ItemINTMILKCOMPINT+MYINT+MILKCOMPDIGINT+DEINT+DIG
      ValueSEValueSEValueSEValueSEValueSEValueSEValueSE
      Model1234567
      Variable
      MBW = BW0.75; dCP = apparent digested CP; dSTA = digested starch; dNDF = digested NDF.
       Intercept−0.198
      Not different from 0 (P > 0.20).
      1.80−2.30
      Not different from 0 (P > 0.20).
      1.90−3.522.16−4.782.250.719
      Not different from 0 (P > 0.20).
      2.025.022.58−2.344.84
       MBW, kg0.750.1220.0180.1650.0180.1290.0180.1320.0180.1470.0200.1140.0200.1110.019
       DMI, kg/d0.5600.0480.5910.480.5770.0470.5980.0510.6000.048
       Milk fat, kg/d1.830.66
       Milk protein, kg/d6.171.13
       Milk yield, kg/ kg of DMI1.540.57
       Milk fat, kg/ kg of DMI16.412.0
       Milk protein, kg/kg of DMI40.223.0
       dCP, kg/d1.630.37
       dSTA, kg/d0.7620.167
       DE, Mcal/kg of DM−1.720.59
       CP digestibility, %−0.04320.027
       NDF digestibility, %−0.06660.018
       Starch digestibility, %0.09100.046
      Fit statistics
      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).
       n293293293293277261277
       BIC1,1611,1781,1591,1481,1571,0451,118
       CCC0.890.880.890.890.860.890.90
       uCCC0.840.810.850.840.790.850.84
      σ^c0.610.600.590.560.680.620.64
      σ^s1.591.821.511.622.021.471.63
      σ^e1.381.441.361.371.521.381.33
      a Not different from 0 (P > 0.20).
      1 MBW = BW0.75; dCP = apparent digested CP; dSTA = digested starch; dNDF = digested NDF.
      2 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).
      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 (
      • Brody S.
      Bioenergetics and Growth..
      ;
      • Kleiber M.
      The Fire of Life: An Introduction to Animal Energetics.
      ), 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
      • NRC
      Nutrient Requirements of Dairy Cattle.
      . However,
      • Moraes L.E.
      • Kebreab E.
      • Strathe A.B.
      • Dijkstra J.
      • France J.
      • Casper D.P.
      • Fadel J.G.
      Multivariate and univariate analysis of energy balance data from lactating dairy cows.
      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 (
      • Dong L.F.
      • Yan T.
      • Ferris C.
      • McDowell D.
      Comparison of maintenance energy requirement and energetic efficiency between lactating Holstein-Friesian and other groups of dairy cows.
      ). 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 (
      • Dong L.F.
      • Yan T.
      • Ferris C.
      • McDowell D.
      Comparison of maintenance energy requirement and energetic efficiency between lactating Holstein-Friesian and other groups of dairy 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 (
      • Purwanto B.P.
      • Abo Y.
      • Sakamoto R.
      • Furumoto F.
      • Yamamoto S.
      Diurnal Patterns of heat-production and heart-rate under thermoneutral conditions in Holstein Friesian cows differing in milk-production.
      ;
      • Reynolds C.K.
      • Tyrrell H.F.
      • Reynolds P.J.
      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.
      ). Up to 50% of total body O2 consumption in ruminants can be attributed to metabolism in the gastrointestinal tract and liver (
      • Seal C.J.
      • Reynolds C.K.
      Nutritional implications of gastrointestinal and liver metabolism in ruminants.
      ). In sheep, increasing DMI increased digestive tract and organ weight as a proportion of BW (
      • McLeod K.R.
      • Baldwin R.L.
      Effects of diet forage:concentrate ratio and metabolizable energy intake on visceral organ growth and in vitro oxidative capacity of gut tissues in sheep.
      ). In growing heifers, increased DMI increased total O2 consumption by the portal-drained viscera and liver (
      • Reynolds C.K.
      • Tyrrell H.F.
      • Reynolds P.J.
      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 (
      • Baldwin R.L.
      Modeling Ruminant Digestion and Metabolism..
      ), 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 (
      • Dado R.G.
      • Mertens D.R.
      • Shook G.E.
      Metabolizable energy and absorbed protein requirements for milk component production.
      ); whereas, conversion of dietary fat or tissue fat into milk fat has an energetic efficiency of 0.94 to 0.97 (
      • Baldwin B.R.
      • Forsberg N.E.
      • Hu C.Y.
      Potential for altering energy partition in the lactating cow.
      ). Given that milk fat is 50% dietary or tissue in origin (
      • Bauman D.E.
      • Griinari J.M.
      Nutritional regulation of milk fat synthesis.
      ), efficiency of milk fat synthesis is 0.83. Synthesis of milk protein is estimated to require approximately 2 Mcal of ME/kg (
      • Dado R.G.
      • Mertens D.R.
      • Shook G.E.
      Metabolizable energy and absorbed protein requirements for milk component production.
      ), 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 (
      • Lemosquet S.
      • Guinard-Flament J.
      • Raggio G.
      • Hurtaud C.
      • Van Milgen J.
      • Lapierre H.
      How does increasing protein supply or glucogenic nutrients modify mammary metabolism in lactating dairy cows?.
      ) to 300% of milk protein yield in lactating dairy goats (
      • Hanigan M.D.
      • France J.
      • Mabjeesh S.J.
      • McNabb W.C.
      • Bequette B.J.
      High rates of mammary tissue protein turnover in lactating goats are energetically costly.
      ), 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 (
      • Holter J.B.
      • Heald C.W.
      • Colovos N.F.
      Heat increments of steam-volatile fatty acids infused separately and in a mixture into fasting cows.
      ), and catabolism of absorbed AA contributes to HI (
      • Reynolds C.
      Splanchnic amino acid metabolism in ruminants.
      ). 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 (
      • NRC
      Nutritional Energetics of Domestic Animals and Glossary of Energy Terms..
      ) 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 (
      • Reynolds C.K.
      • Tyrrell H.F.
      • Reynolds P.J.
      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.
      • Tyrrell H.
      • Moe P.W.
      • Flatt W.P.
      Influence of excess protein intake on energy metabolism of the dairy cow.
      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 (
      • Reed K.F.
      • Bonfa H.C.
      • Dijkstra J.
      • Casper D.P.
      • Kebreab E.
      Estimating the energetic cost of feeding excess dietary nitrogen to dairy cows.
      ). 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 (
      • Kleiber M.
      The Fire of Life: An Introduction to Animal Energetics.
      ).
      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 σe 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 (
      • Reynolds C.K.
      • Tyrrell H.F.
      • Reynolds P.J.
      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.
      ;
      • Cantalapiedra-Hijar G.
      • Lemosquet S.
      • Rodriguez-Lopez J.M.
      • Messad F.
      • Ortigues-Marty I.
      Diets rich in starch increase the posthepatic availability of amino acids in dairy cows fed diets at low and normal protein levels.
      ). 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 (
      • Tyrrell H.
      • Moe P.W.
      • Flatt W.P.
      Influence of excess protein intake on energy metabolism of the dairy cow.
      ;
      • Reed K.F.
      • Bonfa H.C.
      • Dijkstra J.
      • Casper D.P.
      • Kebreab E.
      Estimating the energetic cost of feeding excess dietary nitrogen to dairy cows.
      ). 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 σe 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
      MBW = BW0.75.
      ItemHeat production
      ValueSE
      Variable
      MBW = BW0.75.
       Intercept−5.302.28
       MBW, kg0.750.1200.19
       DMI, kg/d0.6020.52
       Milk fat, kg/kg DMI18.612.5
       Milk protein, kg/kg DMI46.424.1
       Urinary N excretion, kg/d5.322.43
      Fit statistics
      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).
       n261
       BIC1,030
       CCC0.89
       uCCC0.84
      σ^c0.52
      σ^s1.65
      σ^e1.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; σ^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).
      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,
      • Reed K.F.
      • Bonfa H.C.
      • Dijkstra J.
      • Casper D.P.
      • Kebreab E.
      Estimating the energetic cost of feeding excess dietary nitrogen to dairy cows.
      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
      • NRC
      Nutrient Requirements of Dairy Cattle.
      . In the
      • Reed K.F.
      • Bonfa H.C.
      • Dijkstra J.
      • Casper D.P.
      • Kebreab E.
      Estimating the energetic cost of feeding excess dietary nitrogen to dairy cows.
      analysis, HP increased from 6.5 to 7.6 Mcal/kg of excess N when accounting for energy intake, MBW, and milk energy output. Coefficients from the
      • Reed K.F.
      • Bonfa H.C.
      • Dijkstra J.
      • Casper D.P.
      • Kebreab E.
      Estimating the energetic cost of feeding excess dietary nitrogen to dairy cows.
      analysis are larger than ours, which is likely because some urinary N excretion is a function of maintenance and milk production (
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ). 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 (
      • Spek J.W.
      • Dijkstra J.
      • van Duinkerken G.
      • Hendriks W.H.
      • Bannink A.
      Prediction of urinary nitrogen and urinary urea nitrogen excretion by lactating dairy cattle in northwestern Europe and North America: A meta-analysis.
      ), and the energetic cost of ureagenesis from ammonia is estimated to be 3.8 Mcal/kg of N (
      • Martin A.K.
      • Blaxter K.L.
      The energy cost of urea synthesis in sheep.
      ). 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 (
      • Lobley G.E.
      Control of the metabolic fate of amino acids in ruminants: A review.
      ;
      • Firkins J.
      • Reynolds C.
      Whole animal nitrogen balance in cattle.
      ;
      • Reynolds C.
      Splanchnic amino acid metabolism in ruminants.
      ).
      • Reynolds C.
      Splanchnic amino acid metabolism in ruminants.
      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
      • Firkins J.
      • Reynolds C.
      Whole animal nitrogen balance in cattle.
      , 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 (
      • Moorby J.M.
      • Theobald V.J.
      Short communication: The effect of duodenal ammonia infusions on milk production and nitrogen balance of the dairy cow.
      ) or increase liver O2 consumption (
      • Lobley G.E.
      • Connell A.
      • Lomax M.A.
      • Brown D.S.
      • Milne E.
      • Calder A.G.
      • Farningham D.A.
      Hepatic detoxification of ammonia in the ovine liver: Possible consequences for amino acid catabolism.
      ;
      • Firkins J.
      • Reynolds C.
      Whole animal nitrogen balance in cattle.
      ). 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% (
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ) 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 (
      • Spek J.W.
      • Dijkstra J.
      • van Duinkerken G.
      • Hendriks W.H.
      • Bannink A.
      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;
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ); 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 (
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ;
      • Ferraretto L.F.
      • Crump P.M.
      • Shaver R.D.
      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.

      Supplementary Material

      REFERENCES

        • Baldwin B.R.
        • Forsberg N.E.
        • Hu C.Y.
        Potential for altering energy partition in the lactating cow.
        J. Dairy Sci. 1985; 68: 3394-3402
        • Baldwin R.L.
        Modeling Ruminant Digestion and Metabolism..
        Chapman and Hall, London, UK1995
        • Bauman D.E.
        • Griinari J.M.
        Nutritional regulation of milk fat synthesis.
        Annu. Rev. Nutr. 2003; 23 (12626693): 203-227
        • Belyea R.L.
        • Adams M.W.
        Energy and nitrogen utilization of high versus low producing dairy cows.
        J. Dairy Sci. 1990; 73 (2345194): 1023-1030
        • Boerman J.P.
        • Firkins J.L.
        • St-Pierre N.R.
        • Lock A.L.
        Intestinal digestibility of long-chain fatty acids in lactating dairy cows: A meta-analysis and meta-regression.
        J. Dairy Sci. 2015; 98 (26409970): 8889-8903
        • Brody S.
        Bioenergetics and Growth..
        Reinhold Publishing Corporation, New York, NY1945
        • Brouwer E.
        Report of sub-committee on constants and factors.
        in: Proc. the third Symposium of Energy Metabolism, Troon, Scotland. Academic Press, London, UK1965: 441-443
        • Cantalapiedra-Hijar G.
        • Lemosquet S.
        • Rodriguez-Lopez J.M.
        • Messad F.
        • Ortigues-Marty I.
        Diets rich in starch increase the posthepatic availability of amino acids in dairy cows fed diets at low and normal protein levels.
        J. Dairy Sci. 2014; 97 (24931534): 5151-5166
        • Dado R.G.
        • Mertens D.R.
        • Shook G.E.
        Metabolizable energy and absorbed protein requirements for milk component production.
        J. Dairy Sci. 1993; 76 (8326031): 1575-1588
        • Dong L.F.
        • Yan T.
        • Ferris C.
        • McDowell D.
        Comparison of maintenance energy requirement and energetic efficiency between lactating Holstein-Friesian and other groups of dairy cows.
        J. Dairy Sci. 2015; 98 (25483199): 1136-1144
        • Drehmel O.R.
        • Brown-Brandl T.M.
        • Judy J.V.
        • Fernando S.C.
        • Miller P.S.
        • Hales K.E.
        • Kononoff P.J.
        The influence of fat and hemicellulose on methane production and energy utilization in lactating Jersey cattle.
        J. Dairy Sci. 2018; 101 (29908813): 7892-7906
        • Ferraretto L.F.
        • Crump P.M.
        • Shaver R.D.
        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.
        J. Dairy Sci. 2013; 96 (23164230): 533-550
        • Firkins J.
        • Reynolds C.
        Whole animal nitrogen balance in cattle.
        in: Pfeffer E. Hristov A. Nitrogen and Phosphorus Nutrition of Cattle and Environment. CAB International, Cambridge, MA2005: 167-185
        • Foth A.J.
        • Brown-Brandl T.
        • Hanford K.J.
        • Miller P.S.
        • Garcia Gomez G.
        • Kononoff P.J.
        Energy content of reduced-fat dried distillers grains with solubles for lactating dairy cows.
        J. Dairy Sci. 2015; 98 (26233444): 7142-7152
        • Hanigan M.D.
        • France J.
        • Mabjeesh S.J.
        • McNabb W.C.
        • Bequette B.J.
        High rates of mammary tissue protein turnover in lactating goats are energetically costly.
        J. Nutr. 2009; 139 (19403714): 1118-1127
        • Holter J.B.
        • Heald C.W.
        • Colovos N.F.
        Heat increments of steam-volatile fatty acids infused separately and in a mixture into fasting cows.
        J. Dairy Sci. 1970; 53: 1241-1247
        • Judy J.V.
        • Bachman G.C.
        • Brown-Brandl T.M.
        • Fernando S.C.
        • Hales K.E.
        • Harvatine K.J.
        • Miller P.S.
        • Kononoff P.J.
        Increasing the concentration of linolenic acid in diets fed to Jersey cows in late lactation does not affect methane production.
        J. Dairy Sci. 2019; 102 (30612804): 2085-2093
        • Judy J.V.
        • Bachman G.C.
        • Brown-Brandl T.M.
        • Fernando S.C.
        • Hales K.E.
        • Miller P.S.
        • Stowell R.R.
        • Kononoff P.J.
        Energy balance and diurnal variation in methane production as affected by feeding frequency in Jersey cows in late lactation.
        J. Dairy Sci. 2018; 101 (30268600): 10899-10910
        • Judy J.V.
        • Bachman G.C.
        • Brown-Brandl T.M.
        • Fernando S.C.
        • Hales K.E.
        • Miller P.S.
        • Stowell R.R.
        • Kononoff P.J.
        Reducing methane production with corn oil and calcium sulfate: Responses on whole-animal energy and nitrogen balance in dairy cattle.
        J. Dairy Sci. 2019; 102 (30612805): 2054-2067
        • Kleiber M.
        The Fire of Life: An Introduction to Animal Energetics.
        Kreiger Co., Huntington, NY1975
        • Knoell A.L.
        • Judy J.V.
        • Wilson H.C.
        • Herrick K.J.
        • Fernando S.C.
        • Kononoff P.J.
        Energy utilization in lactating Jersey cows consuming a mixture of distillers dried grains with solubles and straw in replacement of alfalfa hay.
        J. Dairy Sci. 2019; 102: T145
        • Kuznetsova A.
        • Brockhoff P.B.
        • Christensen R.H.B.
        lmerTest package: Tests in linear mixed effects models: Version 3.0-1.
        https://cran.r-project.org/web/packages/lmerTest/index.html
        Date: 2017
        Date accessed: November 28, 2019
        • Lemosquet S.
        • Guinard-Flament J.
        • Raggio G.
        • Hurtaud C.
        • Van Milgen J.
        • Lapierre H.
        How does increasing protein supply or glucogenic nutrients modify mammary metabolism in lactating dairy cows?.
        in: Energy and Protein Metabolism and Nutrition. Wageningen Academic Publishers, Wageningen, the Netherlands2010: 175-186
        • Lin L.I.
        A concordance correlation coefficient to evaluate reproducibility.
        Biometrics. 1989; 45 (2720055): 255-268
        • Lobley G.E.
        Control of the metabolic fate of amino acids in ruminants: A review.
        J. Anim. Sci. 1992; 70 (1429303): 3264-3275
        • Lobley G.E.
        • Connell A.
        • Lomax M.A.
        • Brown D.S.
        • Milne E.
        • Calder A.G.
        • Farningham D.A.
        Hepatic detoxification of ammonia in the ovine liver: Possible consequences for amino acid catabolism.
        Br. J. Nutr. 1995; 73 (7626587): 667-685
        • Martin A.K.
        • Blaxter K.L.
        The energy cost of urea synthesis in sheep.
        in: Proc. the 3rd Symposium of Energy Metabolism, Troon, Scotland. Academic Press, London, UK1965: 83-91
        • McLean J.
        • Tobin G.
        Animal and Human Calorimetry..
        Cambridge University Press, Cambridge, UK1987
        • McLeod K.R.
        • Baldwin R.L.
        Effects of diet forage:concentrate ratio and metabolizable energy intake on visceral organ growth and in vitro oxidative capacity of gut tissues in sheep.
        J. Anim. Sci. 2000; 78 (10764085): 760-770
        • Moe P.W.
        • Tyrrell H.F.
        • Flatt W.P.
        Net energy value of feeds for lactation.
        J. Dairy Sci. 1972; 55: 945-958
        • Moorby J.M.
        • Theobald V.J.
        Short communication: The effect of duodenal ammonia infusions on milk production and nitrogen balance of the dairy cow.
        J. Dairy Sci. 1999; 82 (10575611): 2440-2442
        • Moraes L.E.
        • Kebreab E.
        • Strathe A.B.
        • Dijkstra J.
        • France J.
        • Casper D.P.
        • Fadel J.G.
        Multivariate and univariate analysis of energy balance data from lactating dairy cows.
        J. Dairy Sci. 2015; 98 (25892698): 4012-4029
        • Morris D.L.
        • Brown-Brandl T.M.
        • Hales K.E.
        • Harvatine K.J.
        • Kononoff P.J.
        Effects of high-starch or high-fat diets formulated to be isoenergetic on energy and nitrogen partitioning and utilization in lactating Jersey cows.
        J. Dairy Sci. 2020; 103 (32197845): 4378-4389
        • Morris D.L.
        • Judy J.V.
        • Kononoff P.J.
        Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming diets with increasing inclusion of hydrolyzed feather meal.
        J. Dairy Sci. 2020; 103 (32197849): 4206-4217
        • NRC
        Nutritional Energetics of Domestic Animals and Glossary of Energy Terms..
        National Academies Press, Washington, DC1981
        • NRC
        Nutrient Requirements of Dairy Cattle.
        7th rev. ed. National Academies Press, Washington, DC2001
        • Purwanto B.P.
        • Abo Y.
        • Sakamoto R.
        • Furumoto F.
        • Yamamoto S.
        Diurnal Patterns of heat-production and heart-rate under thermoneutral conditions in Holstein Friesian cows differing in milk-production.
        J. Agric. Sci. 1990; 114: 139-142
        • Reed K.F.
        • Bonfa H.C.
        • Dijkstra J.
        • Casper D.P.
        • Kebreab E.
        Estimating the energetic cost of feeding excess dietary nitrogen to dairy cows.
        J. Dairy Sci. 2017; 100 (28711249): 7116-7126
        • Reynolds C.
        Splanchnic amino acid metabolism in ruminants.
        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
        • Reynolds C.K.
        • Tyrrell H.F.
        • Reynolds P.J.
        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.
        J. Nutr. 1991; 121 (2051243): 994-1003
        • Reynolds M.A.
        • Brown-Brandl T.M.
        • Judy J.V.
        • Herrick K.J.
        • Hales K.E.
        • Watson A.K.
        • Kononoff P.J.
        Use of indirect calorimetry to evaluate utilization of energy in lactating Jersey dairy cattle consuming common coproducts.
        J. Dairy Sci. 2019; 102 (30343910): 320-333
        • Roman-Garcia Y.
        • White R.R.
        • Firkins J.L.
        Meta-analysis of postruminal microbial nitrogen flows in dairy cattle. I. Derivation of equations.
        J. Dairy Sci. 2016; 99 (27448861): 7918-7931
        • Seal C.J.
        • Reynolds C.K.
        Nutritional implications of gastrointestinal and liver metabolism in ruminants.
        Nutr. Res. Rev. 1993; 6 (19094308): 185-208
        • Spek J.W.
        • Dijkstra J.
        • van Duinkerken G.
        • Hendriks W.H.
        • Bannink A.
        Prediction of urinary nitrogen and urinary urea nitrogen excretion by lactating dairy cattle in northwestern Europe and North America: A meta-analysis.
        J. Dairy Sci. 2013; 96 (23664347): 4310-4322
        • Tyrrell H.
        • Moe P.W.
        • Flatt W.P.
        Influence of excess protein intake on energy metabolism of the dairy cow.
        in: Proc. the 5th EAAP Symposium on Energy Metabolism of Farm Animals. Juris Verlag, Zurich, Switzerland1970: 69-72
        • White R.R.
        • Hall M.B.
        • Firkins J.L.
        • Kononoff P.J.
        Physically adjusted neutral detergent fiber system for lactating dairy cow rations. I: Deriving equations that identify factors that influence effectiveness of fiber.
        J. Dairy Sci. 2017; 100 (28987580): 9551-9568
        • White R.R.
        • Roman-Garcia Y.
        • Firkins J.L.
        • Kononoff P.
        • VandeHaar M.J.
        • Tran H.
        • McGill T.
        • Garnett R.
        • Hanigan M.D.
        Evaluation of the National Research Council (2001) dairy model and derivation of new prediction equations. 2. Rumen degradable and undegradable protein.
        J. Dairy Sci. 2017; 100 (28259406): 3611-3627