Dietary starch levels modulate production performance and whole-body nutrient metabolism in dairy cows
Usman Arshad, Martin Huser, Mario A. Barrientos-Blanco, Md Zakirul Islam, Xiaoqi Ma, Rong Peng, Kai Wang, Mutian Niu

TL;DR
Increasing dietary starch in dairy cows improves feed efficiency and energy use by shifting metabolism toward carbohydrates.
Contribution
This study shows that higher starch diets enhance energy utilization and production in lactating cows without reducing dry matter intake.
Findings
Higher starch diets increased milk yield and energy intake without affecting dry matter intake.
Starch increased carbohydrate oxidation and metabolic heat production while reducing fat oxidation.
Milk true protein yield increased, and milk urea nitrogen decreased with higher starch.
Abstract
Summary: The objective of this study was to evaluate the effects of increasing dietary starch concentrations on production performance and whole-body nutrient metabolism in lactating dairy cows. Twelve multiparous mid-lactation Holstein cows were assigned to a replicated 3 × 3 Latin square design with three 21-day periods, each consisting of 14 days of dietary adaptation followed by 7 days of sampling. Cows were fed low-, medium-, or high-starch diets, and dry matter intake (DMI) and milk production were recorded, whereas gas exchange was measured using a head chamber (GreenFeed system) to estimate net carbohydrate (COXnet) and fat oxidation (FOXnet), and metabolic heat production (mHP). Increasing dietary starch enhanced production performance and feed efficiency without affecting DMI, and shifted whole-body substrate oxidation toward carbohydrates, resulting in increased COXnet,…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsRuminant Nutrition and Digestive Physiology · Effects of Environmental Stressors on Livestock · Reproductive Physiology in Livestock
Optimizing nutrient utilization in lactating dairy cows is essential for improving production efficiency and reducing the environmental footprint of dairy farming. Among dietary macronutrients, starch is a major determinant of ruminal fermentation patterns, systemic energy metabolism, and milk synthesis in ruminants (Huntington, 1997; Boerman et al., 2015). As a rapidly fermentable carbohydrate in the rumen, it is plausible to speculate that dietary starch may influence the balance between net carbohydrate oxidation (COX_net_) and net fat oxidation (FOX_net_), with potential implications for the partitioning of energy toward milk synthesis and maintenance. Investigating these systemic metabolic effects requires assessment of gas exchange and nutrient oxidation in animals.
Indirect calorimetry via respiration chambers is considered the “gold standard” for assessing whole-animal energy metabolism by quantifying O_2_ consumption, CO_2_ and CH_4_ emissions, and metabolic heat production (mHP; Troy et al., 2016). These measurements enable detailed insights into nutrient oxidation and energy partitioning. Technological advances such as GreenFeed (GF) system, a head chamber technique, have provided noninvasive, real-time assessments of gas exchange, allowing metabolic measurements under more natural and less restrictive feeding environments (Hristov et al., 2015). Such measurements may support investigations into how dietary interventions influence systemic substrate oxidation, mHP, and nutrient partitioning in lactating cows. The present experiment aimed to evaluate the effects of increasing dietary starch concentrations on whole-animal nutrient metabolism using GF in lactating dairy cows. We hypothesized that increasing dietary starch content would alter mHP and substrate oxidation patterns, ultimately affecting nutrient partitioning relevant to lactational performance and feed efficiency in dairy cows.
All experimental processes involving live animals were approved by the Cantonal Veterinary Office of Kanton Zürich (ZH150/2024). The experiment was conducted with dairy cows at the AgroVet-Strickhof Farm (Lindau, Switzerland). Twelve multiparous lactating cows (mean ± SD: lactation number 3.50 ± 1.17; 89 ± 13 DIM; 36.8 ± 3.5 kg/d milk yield; 726 ± 48 kg BW) were assigned to a replicated 3 × 3 Latin square design with 21-d periods. Cows were blocked into squares (n = 4) based on their 7-d average milk yield measured before the start of the experiment to ensure balanced production levels across treatments. Each period included 14 d of dietary adaptation followed by 7 d of sample collection. Treatments consisted of low-starch (LS; 9.66% of DM), medium-starch (MS; 20.6% of DM), and high-starch (HS; 29.1% of DM) diets. The rationale for using a wide range of dietary starch concentrations was to create a broad physiological gradient that would allow us to capture metabolic adaptations across different levels of starch availability. The ingredients and chemical composition of the diets are presented in Table 1. The total duration of the experiment was 63 d. During the first 14 d of each experimental period, the cows were housed in a freestall barn with wheat straw bedding and equipped with a feeding trough (Waagen Döhrn GmbH, Germany) and a Calan-gate system (American Calan Inc., Northwood, NH) for individual feeding. The diet was fed ad libitum at 110% of the previous day's intake, and the TMR was mixed before each morning feeding. Cows were milked twice daily, at 0530 and 1530 h, and their milk yield was recorded by the system. Cows were fed twice daily at 0700 (60%) and 1600 (40%) h. On d 15, cows were moved to a tiestall barn with individual access to drinking water and individual feeding plates (Mettler Toledo GmbH, Greifensee, Switzerland), and samples were collected. Samples of individual feedstuffs and TMR were collected in each sampling week and processed for DM and ingredients composition. Milk yield and DMI were monitored daily throughout the sample collection periods. A total of 4 milk samples per period were collected for composition analysis: the evening of d 5, the morning and evening of d 6, and the morning of d 7 of the sampling week. Milk samples were preserved and analyzed for composition. Milk urea N was calculated as urea concentration (mg/dL) × 0.466 and milk true protein was calculated as milk CP (%) − MUN (mg/dL) × 6.38/1,000. Yield of ECM was calculated according to NASEM (2021): [(12.82 × fat yield) + (7.13 × true protein yield) + (0.323 × milk yield)], and feed efficiency was calculated as the ratio of ECM yield to DMI. Body weight was recorded at the start and end of each period of the experiment.Table 1. Ingredient and nutrient composition (mean ± SD1) of diets fed to mid-lactating Holstein cowsItemTreatment diet2LSMSHSIngredient,3 % of DM Grass silage23.737.135.1 Corn silage16.84.779.50 Alfalfa haylage36.916.44.53 UFA-2497.387.9611.3 UFA-1737.758.959.96 Mineralized concentrate7.387.103.62 Ground corn, dried—17.726.0Diet composition, DM basis NEL,4 Mcal/kg of DM1.601.621.65 DM, %57.4 ± 1.0755.1 ± 1.0154.4 ± 1.04 CP, %15.9 ± 0.4914.7 ± 0.6414.7 ± 0.99 NDF, %47.9 ± 3.2238.3 ± 1.2632.9 ± 0.95 Forage NDF,4 %42.533.827.2 ADF, %31.0 ± 1.1525.8 ± 0.9521.1 ± 2.10 Starch, %9.66 ± 0.0720.6 ± 1.2129.1 ± 1.71 Fatty acids,4 %1.832.202.531Means of 3 composite samples analyzed.2Mid-lactating Holstein cows were fed diets containing low starch (LS), medium starch (MS), or high starch (HS).3UFA-249 and UFA-173 (UFA, Sursee, Switzerland), and mineralized concentrate (no. 50109, Getreidesammelstelle, Thalheim, Switzerland; % DM basis) contained DM = 86.9, 88.0, 90.0; NDF = 9.00, 15.9, 2.60; ADF = 24.4, 7.40, 1.00; CP = 39.0, 27.3, 34.16; starch = 17.0, 32.4; 17.0; crude fat = 5.50, 8.50, 3.19; and ash = 6.50, 7.88, 7.89, respectively.4Calculated using the NASEM Dairy (NASEM, 2021) according to the chemical composition of the dietary ingredients and adjusted for observed DMI for each treatment.
Gas exchange data, including CH_4_, CO_2_, and O_2_, were measured with the GF system (C-Lock Technology Inc., Rapid City, SD) in a tiestall barn. Gas samples were collected on the last day of sampling week in each period at 1000, 1600, 2200, and 0400 h. The system was equipped with an automated collection unit to capture gas samples with a default sample flow rate of 1 L/min operated by a rotameter (10A6130 Glass Tube Purge and Low Flow Meters, Dwyer Instruments Inc.). Each gas collection lasted 5 min, with the GF system collecting background barn air for 2 min between cows to flush the system. During each collection, the cow received a maximum of 8 drops of bait feed, each weighing 30 g, in the form of pellets (UFA-118, UFA, Sursee, Switzerland). Gas exchange data were collected only when proximity sensors confirmed that the cow's head was inside the GF head chamber, and gases were monitored in real time using the Control Feed mobile application (C-Lock Technology Inc., Rapid City, SD) as described previously (Barrientos-Blanco et al., 2025). All parameters related to metabolic gas exchange and nutrient oxidation were calculated based on the average values of gas measurements collected across 4 times during the sampling week of each period. The concentrations of CH_4_ and CO_2_ measured using GF were recorded in grams per day and subsequently converted to liters per day using their respective molecular weights and gas constants. Specifically, CO_2_ (g/d) and CH_4_ (g/d) were converted to volume (L/d) using the molecular weight of each gas and the molar volume at standard temperature and pressure (22.4 L/mol), assuming ideal gas behavior. Fermentative CO_2_ (fCO_2_) was calculated according to Chwalibog et al. (1996) as
whereas metabolic CO_2_ (mCO_2_) was determined by subtracting fCO_2_ from total CO_2_ production as reported (Chwalibog et al., 1996; Kennedy and Kuhla, 2023):
The respiratory quotient (RQ) was calculated as the ratio of metabolic CO_2_ to O_2_ consumption over a 24-h period as reported (Kennedy and Kuhla, 2023):
where the COX_net_ and FOX_net_ were calculated using stoichiometric equations previously described by Kennedy and Kuhla (2023):
Heat production (HP) was estimated based on observed O_2_ consumption and total CO_2_ production, corrected for CH_4_ production and urinary N, using Brouwer's (1965) equation:
and urinary N was calculated according to the Huhtanen et al. (2015) using following equation:
Mean ± SD values of CH_4_ (492 vs. 515 vs. 440 ± 59.4 g/d), CO_2_ (14,403 vs. 15,394 vs. 15,780 ± 1,243 g/d), O_2_ (9,850 vs. 10,601 vs. 10,726 ± 847 g/d), and urinary N (222 vs. 203 vs. 194 ± 15.8 g/d) were observed or used in LS-, MS-, and HS-fed cows, respectively. All metabolic outputs such as COX_net_, FOX_net_, and HP were further normalized to metabolic BW (mBW), defined as BW^0.75^, and expressed both in absolute values and per kilogram of mBW.
A priori power analysis was performed using the pwr4exp package in R (version 4.5.1, R Core Team) following Wang et al. (2025). Three dietary starch levels (LS, MS, HS) were evaluated using preplanned orthogonal polynomial contrasts for linear and quadratic trends. Expected treatment means of milk yield were 34.8, 36.8, and 38.8 kg/d, and error variance (σ^2^ = 9.99 kg^2^/d^2^) was derived from the SEM reported by Boerman et al. (2015). In a replicated Latin square with 4 squares (12 cows total, 4 per treatment, 36 experimental units), the power to detect the linear contrast at α = 0.05 was 0.83. Normality of residuals and homogeneity of variance were assessed for each continuous dependent variable after fitting the statistical model, and none of the response variables violated the normality assumption. Data were analyzed by ANOVA with linear mixed model using the MIXED procedure of SAS (version 9.4, SAS Institute Inc.) using the following model:
where Y_ijkl_ is the dependent variable, μ is the overall mean, T_i_ is the fixed effect of treatment (i = 1 to 3), P_j_ is the fixed effect of period (j = 1 to 3), S_k_ is the fixed effect of square (k = 1 to 4), T_i_ × P_j_ is the interaction between T_i_ and P_j_, T_i_ × S_k_ is the interaction between T_i_ and S_k_, P_j_ × S_k_ is the interaction between P_j_ and S_k_, C(S)lk is the random effect of cow (l = 1 to 3) nested within square (k = 1 to 4), and e_ijkl_ is the random residual. In all mixed models, the Kenward-Roger method was used to approximate the denominator df for the F tests. The covariance structure was modeled based on fit statistics using the smallest corrected Akaike information criterion (AIC). The AR(1) covariance structure produced the lowest AIC and was used for all variables. In all analyses, single df orthogonal polynomial contrasts included to determine linear and quadratic effects of treatments on the response analyzed. The interactive matrix language procedure of SAS (PROC IML) was used to generate the orthogonal coefficients to adjust for the unequal spacing of starch in treatments. The resulting coefficients for LS, MS, and HS were, respectively, −0.735, 0.059, and 0.676 for the linear contrast, and 0.356, −0.814, and 0.458 for the quadratic contrast. Data are reported as LSM and SEM. Evidence of statistical significance against the null hypothesis was considered at P ≤ 0.05, and tendency was considered at 0.05 < P ≤ 0.10.
All 12 cows completed the 63-d experimental period and contributed data for statistical analysis. Data on production performance and feed efficiency are presented in Table 2. Increasing dietary starch concentration linearly (P ≤ 0.03) increased net energy intake (38.9 vs. 40.3 vs. 40.7 ± 0.87 Mcal/d) and milk yield (33.5 vs. 36.6 vs. 38.9 ± 1.24 kg/d), and tended (P = 0.07) to linearly increase ECM yield (39.1 vs. 43.5 vs. 43.1 ± 1.50 kg/d), without a concurrent change in BW, BW change, or DMI, resulting a tendency for improved feed efficiency (ECM/DMI; 1.62 vs. 1.77 vs. 1.78 ± 0.06 kg/kg; P = 0.06) in dairy cows. Milk fat concentration tended to decrease linearly (P = 0.08), from 4.45% in LS to 4.27% in HS. True protein concentration responded quadratically (P < 0.01), with higher values for MS (3.64%) and HS (3.61%) compared with LS (3.44%). Lactose concentration increased linearly (P = 0.02), from 4.92% in LS to 4.98% in HS. Fat and lactose yields were unaffected by dietary starch level, whereas true protein yield increased linearly (P = 0.05), with values of 1.17, 1.37, and 1.34 ± 0.06 kg/d for LS, MS, and HS, respectively. Milk urea nitrogen decreased linearly (P < 0.01) from 18.1 mg/dL in LS to 10.4 ± 0.56 mg/dL in HS-fed cows. Results on whole-body nutrient metabolism are summarized in Table 3. Respiratory quotient showed a tendency (P = 0.06) for a quadratic pattern, as cows fed LS and MS had similar RQ (LS = 0.90 vs. MS = 0.90 ± 0.01), whereas cows fed HS had an RQ of 0.94 ± 0.01, reflecting a metabolic shift toward greater carbohydrate oxidation. Net carbohydrate oxidation rose linearly (P < 0.01), increasing from 5,432 to 7,276 ± 366 g/d, and tended (P = 0.09) toward a quadratic effect when adjusted for mBW, rising from 39.0 to 52.0 ± 2.6 g/kg mBW for LS and HS, respectively. In contrast, FOX_net_ exhibited a quadratic response (P = 0.03), peaking in MS-fed cows and declining in HS-fed cows, with values of 742, 906, and 450 ± 150 g/d for LS, MS, and HS, respectively. A similar pattern was observed when FOX_net_ was normalized to mBW. Heat production increased linearly (P < 0.01), rising from 34.8 to 38.1 ± 0.73 Mcal/d, and showed a similar linear increase when expressed per unit of mBW for LS and HS, respectively.Table 2. Effect of dietary starch level on DMI, productive performance, and feed efficiency in mid-lactating Holstein cowsItemTreatment1SEMPolynomial contrast2LSMSHSLinearQuadraticBW, kg7267337249.50.900.13BW change,3 kg/d0.70−0.820.390.870.700.17DMI, kg/d24.324.924.70.530.400.35Net energy intake, Mcal/d38.940.340.70.870.030.52Concentration, % Fat4.454.474.270.110.080.17 True protein3.443.643.610.06<0.01<0.01 Lactose4.924.924.980.050.020.07Yield, kg/d Milk33.536.638.91.24<0.010.83 Fat1.521.681.600.080.510.30 True protein1.171.371.340.060.050.19 Lactose1.691.861.860.070.150.44 ECM39.143.543.11.500.070.27MUN, mg/dL18.112.910.40.56<0.010.07Feed efficiency, kg/kg ECM:DMI1.621.771.780.060.060.431Mid-lactating Holstein cows were fed diets containing low starch (LS), medium starch (MS), or high starch (HS).2Linear = linear effect of dietary starch; quadratic = quadratic effect of dietary starch.3Body weight change was calculated during the last day of the sampling week in each period.Table 3. Effect of dietary starch level on nutrient partitioning in mid-lactating Holstein cowsItemTreatment1SEMPolynomial contrast2LSMSHSLinearQuadraticRespiratory quotient0.900.900.940.010.020.06COX_net_, g/d5,4325,8487,276366<0.010.11COX_net_, g/kg mBW39.041.652.02.6<0.010.09FOX_net_, g/d7429064501500.110.03FOX_net_, g/kg mBW5.286.423.291.070.130.03HP, Mcal/d34.837.438.10.73<0.010.20HP, Mcal/kg mBW0.250.260.270.01<0.010.361Mid-lactating Holstein cows were fed diets containing low starch (LS), medium starch (MS), or high starch (HS).2Linear = linear effect of dietary starch; quadratic = quadratic effect of dietary starch.
Increasing dietary starch concentration in mid-lactation dairy cows resulted in coordinated shifts in whole-body energy metabolism and productive performance, demonstrating that graded starch supplementation can effectively modulate nutrient partitioning to support lactation. Despite similar BW and DMI across treatments, starch enrichment increased net energy intake, which was reflected in greater milk and ECM yield. The increase in ECM highlights the ability of HS diets to direct nutrients toward productive functions, indicating improved feed efficiency, and this observation is consistent with the findings of Boerman et al. (2015). Results from this study further suggest that the HS diet may have selected for a ruminal microbial community enriched in propionate-producing bacteria (Li et al., 2022). The greater propionate production would increase its absorption and delivery to the liver, where propionate serves as the primary substrate for hepatic gluconeogenesis. The resulting increase in glucose availability supports greater lactose synthesis, which is the key osmotic driver of milk yield. This enhanced glucose supply also contributes to a broader “carbohydrate sufficiency” signal transmitted through gut-liver, gut-pancreas, and gut-peripheral tissue axes. Consistent with this, propionate-derived glucose in response to increased dietary starch may stimulate insulin secretion (Zhang et al., 2015; Piccioli-Cappelli et al., 2022), and elevated insulin provides a systemic anabolic signal that attenuates lipid mobilization and redirects conserved and newly synthesized energy toward the mammary gland (Peters and Elliot, 1984). In this study, DMI was not affected by the treatments, which is consistent with observations by Allen et al. (2009), who reported that feeding cows higher amounts of starch increased ruminal propionate concentrations, leading to smaller meal sizes and reduced DMI. In contrast, the LS diet contained substantially more forage-derived NDF, which can slow digesta passage and increase rumen fill (Salfer et al., 2018), potentially limiting intake. Notably, increasing dietary starch leads to reduced CH_4_ emissions (Juckem et al., 2025), largely because propionate serves as a hydrogen sink. Thus, it is plausible that improvement in energetic efficiency primarily reflects shifts in rumen fermentation and substrate metabolism, with CH_4_ reduction occurring in parallel rather than directly driving milk synthesis. The effect of starch on protein metabolism was also evident, though secondary to the overarching energetic response. Higher starch intake was associated with increased milk protein concentration and yield, reflecting enhanced ruminal microbial protein synthesis supported by improved synchronization of energy and N availability. The observed decline in MUN with increasing dietary starch further supports reduced ammonia losses and more efficient N capture (Al-Dehneh et al., 1997), indicating coordinated partitioning of both energy and N toward productive functions.
In this study, whole-body energetics measurements provide mechanistic insight into the observed production responses. The linear increase in COX_net_ with increasing starch indicates a direct shift toward glucose-driven energy metabolism. In contrast, FOX_net_ exhibited a quadratic response, peaking at intermediate starch intake and declining at the highest starch level, suggesting progressive displacement of lipid catabolism as carbohydrate availability increased. These patterns are consistent with the proposed endocrine and metabolic signaling framework in which insulin and propionate act as key regulators of substrate selection and nutrient partitioning (Peters and Elliot, 1984; Allen et al., 2009; Zhang et al., 2015). Consistent with these changes, the increased RQ observed with the HS diet likely reflects a preferential oxidation of carbohydrates over lipids to support the energetic demands of lactogenesis. Heat production, both in absolute terms and relative to mBW, increased linearly with dietary starch, paralleling the rise in ECM. Although increasing dietary starch has been shown to increase BW gain in early- and late-lactation cows (Piccioli-Cappelli et al., 2022), we did not observe effects on BW or BW change in this study. However, it is plausible to indicate that the increases in COX_net_, RQ, and HP indicate that the additional starch may have been utilized through enhanced carbohydrate oxidation to support lactational metabolism rather than body tissue retention. Overall, these findings demonstrate that increasing dietary starch enhances efficiency by promoting carbohydrate oxidation and directing energy toward milk synthesis in dairy cows.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Al-Dehneh A.A.Kennedy A.D.Mc Allister T.A.Incorporation of recycled urea—N into ruminal bacteria flowing to the small intestine of dairy cows fed a high-grain or high-forage diet Anim. Feed Sci. Technol.68199732733310.1016/S 0377-8401(97)00059-X · doi ↗
- 2Allen M.S.Bradford B.J.Oba M.Board-invited review: The hepatic oxidation theory of the control of feed intake and its application to ruminants J. Anim. Sci.8720093317333410.2527/jas.2009-17791964850019648500 · doi ↗ · pubmed ↗
- 3Barrientos-Blanco M.A.Arshad U.Giannoukos S.Islam M.Z.Kunz C.Peng R.Räisänen S.E.Zenobi R.Niu M.A sampling method for differentiating breath and ruminal exhaled volatile organic compounds in dairy cows using methane as a marker JDS Commun.6202543844310.3168/jdsc.2024-07324045815540458155 PMC 12126823 · doi ↗ · pubmed ↗
- 4Boerman J.P.Potts S.B.Vande Haar M.J.Allen M.S.Lock A.L.Milk production responses to a change in dietary starch concentration vary by production level in dairy cattle J. Dairy Sci.9820154698470610.3168/jds.2014-89992598107525981075 · doi ↗ · pubmed ↗
- 5Brouwer E.Report of sub-committee on constants and factors Proc. Third Symposium of Energy Metabolism, Troon, Scotland 1965 Academic Press London, UK 441443
- 6Chwalibog A.Jensen K.Thorbek G.Oxidation of nutrients in bull calves treated with β-adrenergic agonists Arch. Tierernahr.49199625526110.1080/1745039960938188889883128988312 · doi ↗ · pubmed ↗
- 7Hristov A.N.Giallongo O.J.Frederick F.Weeks T.Zimmerman H.Harper P.R.Hristova M.T.Zimmerman R.S.Branco A.F.The use of an automated system (Green Feed) to monitor enteric methane and carbon dioxide emissions from ruminant animals J. Vis. Exp.1032015 e 529042638388610.3791/52904 PMC 469258726383886 · doi ↗ · pubmed ↗
- 8Huhtanen P.Cabezas-Garcia S.Krizsan S.J.Shingfield K.J.Evaluation of between-cow variation in milk urea and rumen ammonia nitrogen concentrations and the association with nitrogen utilization and diet digestibility in lactating cows J. Dairy Sci.9820153182319610.3168/jds.2014-82152577106025771060 · doi ↗ · pubmed ↗
