Journal of Dairy Science
Volume 91, Issue 10 , Pages 3814-3823, October 2008

Nutrition, Metabolism, and Fertility in Dairy Cows: 1. Dietary Energy Source and Ovarian Function

University of Nottingham, School of Biosciences, Sutton Bonington Campus, Loughborough LE12 5RD, United Kingdom

Received 18 January 2008; accepted 19 May 2008.

Article Outline

Abstract 

In previous studies, high plasma insulin was associated with earlier resumption of postpartum estrous cycles in dairy cows. The objective of this experiment was to quantify hormonal and ovarian responses to dietary starch and fat contents. Thirty cows were fed on a standard diet from calving until 40 d in milk (DIM) and then 6 cows were allocated to each of 5 isoenergetic diets containing 231, 183, 159, 135, and 87g of starch and 39, 42, 43, 45, and 48g of fat/kg of dry matter (DM) for diets 1 to 5, respectively, until 70 DIM. Estrus was synchronized at 60 DIM. Between 60 and 70 DIM, energy intake, milk yield, and energy balance were similar among diet groups. Plasma insulin-to-glucagon ratio increased with increasing dietary starch and decreasing dietary fat concentrations, reaching a break point at 159g of starch, 43g of fat/kg of DM (diets 1 to 5: mean 3.86, 3.78, 3.59, 2.98, 2.06±standard error 0.22). Growth hormone, insulin-like growth factor-I, and leptin did not vary among diets. The greatest dietary starch concentration was associated with elevated plasma urea-N (diets 1 to 5: mean 3.69, 3.01, 2.94, 2.95, 2.75,±standard error 0.13mmol/L, respectively) and delayed postovulatory progesterone increase (progesterone at 3 to 5 d postovulation for diets 1 to 5: mean 2.7, 5.9, 4.2, 5.6, 4.3±standard error 0.9ng/mL, respectively). The number of small (<5mm) ovarian follicles was positively related to starch intake (r=0.381) and plasma insulin concentration (r=0.402). It is concluded that to maintain adequate insulin-to-glucagon ratio in cows at the start of the breeding period, dietary starch concentration should be above 160g/kg of DM and dietary fat below 44g/kg of DM, and this should have a positive effect on ovarian function.

Key words: insulin-to-glucagon ratio, starch, fatty acid, ovarian function

 

Back to Article Outline

Introduction 

Large increases in milk yield over the past 30 yr were associated with declining fertility of dairy cows all over the world (Butler, 2003). This decline in fertility was attributed partly to unfavorable genetic correlations between milk yield and reproductive traits and partly to increasing imbalance of nutrients leading to metabolic stress (Pryce et al., 2004). Much emphasis has been placed on the strong association between negative energy balance (NEB) in early lactation and length of the postpartum anovulatory period (Garnsworthy et al., 2008). Prolonged periods of NEB were associated with suppression of pulsatile LH secretion, reduced ovarian responsiveness to LH stimulation, and reduced estradiol secretion by the dominant follicle, all of which influenced ovulation of the dominant follicle (Butler, 2003). Mobilization of body fat during NEB increased plasma concentrations of NEFA and BHBA, both of which were associated with reduced fertility (Garnsworthy et al., 2008).

Negative energy balance resulted in loss of BCS as the cow mobilized body fat reserves to support milk production. Greater BCS loss was associated with delayed first ovulation postpartum and reduced conception rate (Butler, 2003). The magnitude and duration of BCS loss was directly related to BCS at calving, because dairy cows adjust their DMI in early lactation to move toward a biological target BCS at around 12wk postpartum (Garnsworthy and Topps, 1982). A recent review (Garnsworthy, 2007) suggested that biological BCS targets were defined by genetics and have reduced over the past 20 yr; therefore, modern dairy cows are more likely to suffer prolonged NEB.

The most common strategy used to reduce the extent of NEB and BCS loss in early lactation is to increase dietary energy concentration by increasing the starch or fat components of the ration at the expense of forage components. Such changes in carbohydrate and fat supplies have implications for rumen function, milk composition, nutrient partition, and metabolic hormones.

Changes in metabolic hormones are pertinent for dairy cow fertility, because they interact with reproductive hormones that control ovarian function (Webb et al., 2004; Garnsworthy et al., 2008). Selection of dairy cows for high milk production was associated with a longer interval from parturition to first ovulation, high plasma concentrations of GH and BHBA, and low plasma concentrations of glucose and insulin (Gutierrez et al., 2006). In an earlier study, we tested the hypothesis that feeding a high-starch diet to increase circulating insulin concentrations for the first 50 DIM could alleviate the delay in the first ovulation postpartum observed in high genetic merit dairy cows (Gong et al., 2002). The high insulin-inducing diet increased the proportion of cows ovulating within 50 d of calving from 55 to 90% and reduced the interval to first ovulation postpartum from 48 to 34 d. In other studies, we found that developmental competence of oocytes was influenced by high versus low dietary concentrations of starch (Fouladi-Nashta et al., 2005) and fat (Fouladi-Nashta et al., 2007) and that responses were related to changes in metabolic hormones. Although these studies provided proof of principles, comparison of extremes does not allow prediction of responses to intermediate diets.

We have conducted a series of experiments to investigate responses to diets designed to alter metabolic hormones, particularly insulin, in high-yielding dairy cows. The overall objective was to determine hormonal responses to diet at a critical stage of the breeding period (40 to 70 DIM) and to relate these to traits of ovarian function.

The specific objective of this experiment was to investigate responses to isoenergetic diets varying in composition from high starch to high fat in a regression approach. The main hypothesis being tested was that plasma insulin concentration would be influenced by dietary concentrations of starch and fat, which would lead to differences in ovarian function.

Back to Article Outline

Materials and Methods 

Animals, Management, and Diets 

Thirty multiparous Holstein dairy cows were group-housed and individually fed through Calan electronic gates from calving until 70 DIM. Diets were fed as TMR, and daily feed allowances were adjusted for each cow to at least 110% of intake, to ensure ad libitum feeding. Daily feed allowances were offered in 2 portions at 0700 and 1300h each day. Cows were milked twice daily at approximately 0500 and 1530h.

Cows were fed a standard TMR (diet 3; Table 1) from calving until 40 DIM and were assigned to 6 blocks of 5 cows on the basis of calving date. Blocks were matched for milk yield recorded during the first 20 DIM. From 38 DIM, cows were transferred at random within blocks to 5 treatment TMR over a period of 2 d. Treatment TMR were fed from 40 to 70 DIM.

Table 1. Formulation and chemical composition of isoenergetic treatment diets varying from high starch (1) to high fat (5)
Diet
Formulation (g per kg of DM)12345
Grass1 silage334334334334334
Corn2 silage163163163163163
Brewers grains5555555555
Wheat24917714210635
Fatty acid supplement3813151823
Molassed sugar beet pulp3799129159220
Soybean meal7881838588
Rapeseed meal6668697072
Minerals and vitamins466666
Dicalcium phosphate44444
Total1,0001,0001,0001,0001,000
Composition5
DM (g/kg)445445446446447
ME (MJ/kg of DM)12.212.212.212.212.2
NDF (g/kg of DM)301317324332347
Starch (g/kg of DM)23118315913587
Rumen bypass starch (g/kg of DM)8466584932
Sugars (g/kg of DM)58727986100
Fat6 (g/kg of DM)3942434548
CP (g/kg of DM)181181181181181
Effective RDP5 (g/kg of DM)115114113112111
Digestible RUP7 (g/kg of DM)5455565657

1Perennial ryegrass (Lolium perenne).

2Maize (Zea mays).

3Megalac (calcium salts of palm fatty acids); Volac International, Royston, Herts., UK.

4Bibby HiPhos: ABN Ltd., Peterborough, UK: calcium, 18%; phosphorus, 10%; magnesium, 5%; salt, 17%; copper, 2,000mg/kg; manganese, 5,000mg/kg; cobalt, 100mg/kg; zinc, 6,000mg/kg; iodine, 500mg/kg; selenium, 25mg/kg; vitamin A, 400,000IU/kg; vitamin D3, 80,000IU/kg; vitamin E, 1,000mg/kg.

5Calculated from laboratory analysis of components.

6Acid hydrolysis method.

7Calculated from rumen degradability values determined in situ and laboratory analysis of components.

The TMR (Table 1) were formulated to be isoenergetic and isonitrogenous, supplying ME and MP requirements (Thomas, 2004) for 45L/d of milk production. All TMR contained the same proportions of grass silage, corn silage, brewers grains, and concentrates. Two concentrates were formulated: one was designed glucogenic (high starch concentration) and the other ketogenic (high fat and fiber concentrations). Five concentrates were produced by mixing these 2 concentrates in the proportions 1:0, 2:1, 1:1, 1:2, and 0:1. These 5 concentrates were mixed with the forage and brewers grains to produce the 5 TMR.

Feed Analysis 

Feed intake was recorded daily, and diet samples were taken for analysis at weekly intervals. Forages were analyzed by near-infrared spectroscopy using Forage Analysis Assurance Group equations to predict nutrient contents. All other feeds were analyzed using AOAC (1990) methods. Dry matter, CP, NDF, starch, sugars, and fat (oil acid hydrolysis; OAH) were determined directly. Metabolizable energy was calculated as (0.14×NCGD + 0.25×OAH), where NCGD = neutral cellulase plus gamanase digestibility (FSR, 2005). Effective RDP, digestible RUP, and rumen bypass starch were calculated from laboratory-determined CP and starch values, using degradation rates determined by the nylon bag method (Allison and Garnsworthy, 2002) in 4 rumen-fistulated cows.

Estrous Synchronization 

Estrous cycles were synchronized by vaginal insertion of a controlled internal drug release (CIDR) progesterone-releasing device (SmithKline Beecham, Tadworth, Surrey, UK) at 50 DIM, with administration of an injection of 2mL of the PGF analog, cloprostenol (Estrumate; Schering-Plough Animal Health, Welwyn, Garden City, UK), when the CIDR was removed at 60 DIM. Cows were allowed to ovulate without further hormonal intervention.

Ovarian activity was monitored by daily transrectal ultrasound scanning from 61 to 70 DIM, using an Aloka SSD-500 scanner equipped with a 5-MHz linear array transducer (Aloka Co. Ltd., Tokyo, Japan). Day of ovulation (3 to 7 d after CIDR removal) was determined as the day on which a large follicle was no longer present and was confirmed by subsequent progesterone increase. Numbers of small (<5mm) and medium-sized (5 to 10mm) follicles were recorded at each scan, as was the diameter of the largest follicle before ovulation and diameter of the developing corpus luteum after ovulation.

Performance Measurements 

Milk yield was recorded daily from calving until 70 DIM. Milk samples were collected at 30, 35, and 60 to 70 DIM and analyzed for fat, protein, and lactose contents by infrared analysis at the National Milk Records Laboratory, Harrogate, Yorkshire, United Kingdom, using AOAC reference method no. 972.16 (AOAC, 1990). Body weight and BCS (1 to 5 scale; Mulvanny, 1977) were recorded weekly from 30 to 70 DIM at 0900h on Monday of each week.

Blood Sampling and Analysis 

Blood samples were collected by jugular venipuncture at 0900h (2h after morning feeding), on 2 occasions between 30 and 38 DIM for determination of baseline hormones and metabolites. Additional blood samples were collected daily at 0900h from 60 to 70 DIM for determination of daily variation in hormones and metabolites.

Blood samples were analyzed for the following hormones (in each case, the reference is followed by mean assay sensitivity, intraassay CV, and interassay CV): insulin (Adamiak et al., 2005;0.045ng/mL, 4.0%, 8.4%), growth hormone (GH; Gutierrez et al., 2006; 1.2ng/mL, 3.9%, 11.3%), IGF-I (Gutierrez et al., 2006; 0.11ng/mL, 3.8%, 12%), glucagon (Linco Research Inc., St. Charles, MO; 40.5pg/mL, 5.3%, 6.2%), leptin (Adamiak et al., 2005; 0.2ng/mL, 3.9%, 12.2%), estradiol (Mann et al., 1995; 0.27pg/mL, 10%, 9.7%), and progesterone (Mann et al., 1995; 0.3ng/mL, 5.7%, 9.4%).

Blood samples were analyzed for the following metabolites on a Bayer opera autoanalyzer (Bayer UK Ltd., Newbury, UK): albumin (Bayer kit T01 137702), total protein (Bayer kit T01 130102), globulin (total protein minus albumin), urea-N (Bayer kit T01 182356), glucose (Bayer kit T01 183356), BHBA (Randox kit Ranbut RB 1008), NEFA (Waiko kit NEFA-C), magnesium (Bayer kit T01 287802), and inorganic phosphorus (Bayer kit T01 130302); intra- and interassay CV were <5%.

Statistical Analyses 

Metabolizable energy balance was calculated for each cow from BW, milk energy output, and metabolizable energy intake, using UPWin software (AGM Systems, Exeter, UK) and Feed Into Milk equations (Thomas, 2004). Dietary effects on mean performance and blood data during the treatment period (60 to 70 DIM) were examined by analysis of covariance, using Genstat 8 (Lawes Agricultural Trust, Rothamsted, UK), looking at the main effect of dietary treatment group, with observations during the standardization period (30 to 38 DIM) as covariates. When the main effect was significant, differences among treatment means were tested by least significant difference. Linear and quadratic effects of dietary starch concentration were examined. Simple correlation coefficients were determined to examine relationships among variables. The all-subsets regression procedure (Genstat 8, Lawes Agricultural Trust) was used to identify significant multiple relationships among variables, and coefficients for these were determined by multiple linear regression.

Reproductive data, which were not normally distributed, were analyzed by using generalized linear models with a Poisson error distribution and log link function. For progesterone, estradiol, ovulatory follicle diameter, and corpus luteum diameter, day of ovulation was included as a covariate.

Back to Article Outline

Results 

The responses are analyzed and reported as linear and quadratic effects of dietary starch concentration, but they can be interpreted as effects of dietary fat, sugars, or fiber concentrations, all of which changed in a linear fashion across diets.

Intake, Milk Production, and Body Weight 

Intakes of DM, NDF, and ME did not vary among treatment groups (Table 2). Starch intake decreased (P<0.001), whereas intakes of sugars and fat increased (P<0.001), with decreasing dietary starch concentration and increasing sugars and fat concentrations.

Table 2. Least squares means over 60 to 70 DIM, adjusted by covariance for values between 30 and 38 DIM, for intakes, milk production and composition, BW, and BCS in cows (n=6) given diets varying from high starch (1) to high fat (5)
Diet1 Contrast2
Item12345SED3DietLinearQuadratic
Intake
DM (kg/d)23.324.523.224.423.81.830.9270.8460.755
Starch (kg/d)5.4a4.5b3.7c3.3c2.1d0.29<0.001<0.0010.694
Fat (kg/d)1.0a1.1ab1.1ab1.2b1.2b0.090.0430.0070.841
Sugars (kg/d)1.3a1.7b1.8bc2.0cd2.3d0.14<0.001<0.0010.796
NDF (kg/d)7.17.87.58.18.20.600.3400.0580.765
ME (MJ/d)28329628029528722.20.9280.8610.755
Milk production
Milk yield (kg/d)43.743.741.243.644.11.480.1130.8700.121
Milk fat (g/kg)36.2a37.9ab40.4abc43.0c40.7abc2.250.0490.0150.182
(g/d)1,5731,6851,6751,8501,766142.10.3890.1100.513
Milk protein (g/kg)30.029.830.929.228.50.820.0880.0780.157
(g/d)1,3081,2921,2681,2571,25647.60.7580.2150.733
Milk lactose (g/kg)47.447.346.447.047.30.390.0800.5600.061
(g/d)2,0862,0741,8762,0232,08169.20.0680.7360.310
Energy balance (MJ/d)−0.14.43.1−6.20.73.750.0770.4910.869
Weight and condition
BW (kg)6116076136126188.10.7420.2930.485
BW change (kg/d)−0.010.100.08−0.130.050.1440.5150.3430.286
BCS2.62.42.52.42.50.190.6690.7940.199
BCS change0.09−0.15−0.04−0.130.020.1800.6700.7910.195

a–dMeans in the same row without a common superscript differ (P<0.05).

1Diets 1 to 5 contained 231, 183, 159, 135, and 87g of starch and 39, 42, 43, 45, and 48g of fat/kg of DM, respectively.

2Significance of diet effect and linear and quadratic effects of dietary starch concentration.

3SED = standard error of the difference between treatment means.

Milk yield did not vary among treatment groups (Table 2). Milk fat concentration varied among treatment groups (P<0.05) and was linearly related (P<0.05) to dietary starch (negative) and fat (positive) concentrations. Concentrations of milk protein and lactose, and yields of milk fat, protein, and lactose, did not vary among treatment groups. In addition, energy balance, BW, BW change, BCS, and change in BCS did not vary among treatment groups (Table 2).

Metabolic Hormones and Metabolites 

Mean plasma insulin and glucagon concentrations, and their ratio, varied (P<0.001) among treatment groups (Table 3), with significant linear (P<0.001) and quadratic (P<0.05) effects of dietary starch concentration. Effects of dietary starch concentration were positive for insulin, negative for glucagon, and positive for the insulin-to-glucagon ratio. For IGF-I, GH, and leptin, there was no difference between treatment groups and no linear or quadratic effect of dietary starch concentration.

Table 3. Least squares means over 60 to 70 DIM, adjusted by covariance for values between 30 and 38 DIM, for plasma concentrations of metabolic hormones and metabolites in cows (n=6) given diets varying from high starch (1) to high fat (5)
Diet1 Contrasts2
Item12345SED3DietLinearQuadratic
Insulin (ng/mL)0.40a0.39a0.38ab0.34b0.25c0.031<0.001<0.0010.019
Glucagon (pg/mL)102a108a108a115b119b5.60.0370.0020.833
Insulin:glucagon ratio3.86a3.78a3.59a2.98b2.06c0.315<0.001<0.0010.018
IGF-I (ng/mL)110711168810221.30.2560.9210.261
Growth hormone (ng/mL)7.967.267.867.287.011.4740.9570.5480.989
Leptin (ng/mL)1.271.291.431.041.160.1740.2740.3270.609
Albumin (g/L)34.834.437.136.337.41.650.3020.0820.980
Globulin (g/L)36.335.336.435.137.71.790.6160.4950.255
Total protein (g/L)71707372753.010.4460.1660.382
Urea-N (mmol/L)3.69a3.01b2.94b2.95b2.75b0.186<0.001<0.0010.033
BHBA (mmol/L)0.670.700.781.241.150.2800.1360.0270.980
NEFA (mmol/L)0.19a0.18a0.19a0.28b0.32b0.0390.004<0.0010.149
Glucose (mmol/L)3.43.63.63.53.40.150.2800.6720.046
Phosphorus (mmol/L)1.421.511.551.411.510.1200.7480.7290.759
Magnesium (mmol/L)1.000.900.940.990.970.0370.0830.7690.114

a–cMeans in the same row without a common superscript differ (P<0.05).

1Diets 1 to 5 contained 231, 183, 159, 135, and 87g of starch and 39, 42, 43, 45, and 48g of fat/kg DM, respectively.

2Significance of diet effect and linear and quadratic effects of dietary starch concentration.

3SED = standard error of the difference between treatment means.

There was no difference between treatment groups, and no effect of dietary starch concentration, for plasma concentrations of albumin, globulin, total protein, phosphorus, or magnesium (Table 3). Plasma urea-N concentration was greater for group 1 (greatest starch) than for any other group; therefore, the overall dietary effect was significant (P<0.001) and significant positive linear (P<0.001) and quadratic (P<0.05) effects of dietary starch concentration were observed. There was no difference among treatment groups for BHBA, although there was a linear (P<0.05) increase in BHBA with decreasing dietary starch concentration. Plasma NEFA concentrations were lower (P<0.01) in diets 1 to 3 than in diets 4 and 5, and the negative linear effect of dietary starch concentration was significant (P<0.001). There was no overall difference among treatment groups for plasma glucose, but diet groups 1 and 5 had lower glucose concentrations than the other groups, and the quadratic effect of dietary starch concentration was significant (P<0.05).

Ovarian Activity 

Ovulation was detected by ultrasound scanning in 28 of the 30 cows. For 1 cow in diet group 1 and 1 cow in diet group 3, scanning was not possible due to behavioral problems, and ovulation was confirmed by assessment of estradiol and progesterone concentrations. Cows showed a normal pattern of ovulation with a mean (±SE) day of ovulation of 4.3±0.2 (range 3 to 7 d) after CIDR removal.

The maximum number of small (<5mm) follicles during the preovulatory follicular wave was greater in diet groups 1 and 2 (high starch) than in diet groups 4 and 5 (high fat), and the overall treatment effect was significant (P<0.001; Table 4). Similar effects (P<0.01) on the maximum number of small follicles were seen in the postovulatory follicular wave.

Table 4. Treatment means1 relative to a synchronized estrus at around 60 DIM for plasma concentrations of estradiol before ovulation and progesterone from d 3 to 5 postovulation, as well as follicle numbers, ovulatory follicle diameter, and corpus luteum (CL) diameter in cows (n=6) given diets varying from high starch (1) to high fat (5)
Diet2
Item12345SED3P-value4
Progesterone (ng/mL)2.7a5.9b4.2ab5.6b4.3ab1.260.043
Estradiol (pg/mL)1.31.41.11.21.30.260.988
Maximum small (<5mm) follicles preovulation13.0a15.7a11.8ab8.0b8.7b2.03<0.001
Maximum small (<5mm) follicles postovulation14.2ab17.2a11.8b9.5c11.5bc2.140.004
Maximum medium-sized (5 to 10mm) follicles postovulation2.6a3.5a6.0b4.2ab4.8ab1.230.044
Ovulatory follicle diameter (mm)19262423222.70.109
CL diameter5 (mm)26a32b26a29ab28ab3.10.041

a–cMeans in the same row without a common superscript differ (P<0.05).

1Values are back-transformed from generalized linear models with Poisson error distribution and log link function. For progesterone, estradiol, ovulatory follicle diameter, and CL diameter, day of ovulation was included as a covariate.

2Diets 1 to 5 contained 231, 183, 159, 135, and 87g of starch and 39, 42, 43, 45, and 48g of fat/kg of DM, respectively.

3SED = standard error of the difference between treatment means.

4P = probability of main treatment effect (diet) on each variable.

5On d 5 after ovulation.

The maximum number of medium-sized (5 to 10mm) follicles during the postovulatory follicular wave (Table 4) was lower (P<0.05) in diet groups 1 and 2 (high starch and low fat) than in diet group 3; the maximum number of medium-sized follicles for diet groups 4 and 5 (low starch and high fat) was not different from groups 1 to 3. In addition, there was no difference among treatment groups in maximum diameter of the ovulatory follicle (Table 4), but the corpus luteum that developed postovulation was larger for diet group 2 than for diet groups 1 and 3.

After ovulation, cows in diet group 1 (high starch) showed a delayed postovulatory increase in progesterone compared with the other 4 groups (Figure 1). This late progesterone increase resulted in lower mean progesterone concentrations from d 3 to d 5 postovulation for diet group 1 than for diet groups 2 and 4 (Table 4). Plasma concentration of estradiol did not differ among dietary groups on any day of the study, and there was no difference among groups in plasma estradiol measured during the preovulatory phase (d -2 to 0; Table 4).

  • View full-size image.
  • Figure 1. 

    Mean plasma concentrations of progesterone in groups (n=6 per group) of lactating Holstein-Friesian cows fed diets varying from high starch (1) to high fat (5). Cycles were synchronized and data then aligned to day of ovulation. Diets 1 to 5 contained 231, 183, 159, 135, and 87g of starch and 39, 42, 43, 45, and 48g of fat/kg of DM, respectively. aMean progesterone was lower (P=0.039) for diet 1 than for the other diets on d 5.

Relationships Among Variables (Across Individuals) 

Plasma insulin was positively correlated with starch intake (r=0.56), energy balance (r=0.40), milk protein content (r=0.50), and plasma glucose concentration (r=0.45). Insulin was negatively correlated with plasma BHBA (r=−0.39) and plasma NEFA (r=−0.58). Plasma IGF-I was negatively correlated with intakes of DM (r=−0.43) and fat (r=−0.37), and with milk yield (r=−0.38), milk fat concentration (r=−0.40), and milk fat yield (r=−0.53). Plasma leptin was negatively correlated with intakes of DM (r=−0.44) and fat (r=−0.41) and with yields of milk (r=−0.50), fat (r=−0.45), protein (r=−0.39), and lactose (r=−0.48). Insulin-like growth factor-I was positively correlated with plasma leptin (r=0.47). None of the traits measured were correlated with glucagon or GH.

The maximum number of small (<5mm) ovarian follicles during the follicular wave preceding ovulation was positively related to insulin, starch intake, and milk yield and negatively related to milk fat and protein concentrations (Table 5). Apart from milk yield, similar relationships were observed during the first follicular wave postovulation.

Table 5. Effects1 of insulin and production traits on numbers of ovarian follicles <5mm before and after a synchronized ovulation at around 60 DIM
Preovulation Postovulation
TraitConstantSESlopeSERSD2P-valueConstantSESlopeSERSDP-value
Insulin (ng/mL)1.270.2753.3010.7392.10<0.0011.880.2511.9330.6901.860.005
Starch intake (kg/d)1.630.2060.2150.0512.12<0.0011.980.1900.1560.0481.810.001
Milk yield (kg/d)1.450.4820.0230.0112.250.040
Milk fat (g/kg)3.740.488−0.0330.0122.220.0073.450.459−0.0230.0121.870.049
Milk protein (g/kg)5.591.11−0.1070.0382.210.0046.031.050−0.1180.0361.800.007

1Values are regression equations from generalized linear models with Poisson error distribution and log link function, fitted to individual cow means. Each equation predicts log of maximum follicle numbers from the relevant trait; for example, preovulation, log[follicles]=1.45 + (milk yield×0.023); so if milk yield=40kg/d, predicted maximum small follicles=10.7.

2Residual standard deviation.

The best multiple regression equation for predicting plasma insulin concentration included positive effects of starch intake and plasma glucose and a negative effect of milk yield (Table 6). This equation accounted for 53% of variation in plasma insulin. The best multiple regression equation for predicting plasma IGF-I concentration included positive effects of BCS change and BW and a negative effect of DMI. This equation accounted for 44% of variation in plasma IGF-I. The best regression equation for predicting plasma leptin concentration included only a negative effect of milk yield. This equation accounted for 23% of variation in plasma leptin.

Table 6. Best regression equations to predict metabolic hormones from intake, milk production, and blood parameters
ItemConstantSlope 1Slope 2Slope 3RSD1% VAR2
Insulin (ng/mL)
Starch intake (kg/d) + milk yield (kg/d) + plasma glucose (mmol/L)−0.160.046−0.0040.140.05453.4
IGF-I (ng/mL)
Condition score change + BW (kg) + DMI (kg/d)136.155.30.23−7.3427.3143.5
Growth hormone (ng/mL)
BW (kg) + DMI (kg/d)10.7−0.0210.396 2.2510.5
Leptin (ng/mL)
Milk yield (kg/d)2.93−0.039 0.3522.6

1Residual standard deviation.

2Percentage of variation accounted for by regression equation.

Back to Article Outline

Discussion 

The design of the experiment allows examination of responses to isoenergetic diets varying widely in composition across the range found in normal commercial practice. An advantage of providing stepped changes in energy sources from high starch to high fat is that intermediate responses can indicate thresholds. Dry matter intake was not affected by diet; therefore, responses can be attributed to changing nutrients rather than energy intake. At the individual cow level, DMI was correlated with milk yield (r=0.76; P<0.001), indicating that cows were eating to meet their physiological requirements for energy and nutrients.

The main hypothesis was that plasma insulin concentration would be influenced by dietary concentrations of starch and fat. This hypothesis was supported by the differences in insulin observed among treatments and the relationships between insulin and starch concentration (r=0.45), fat concentration (r=−0.37), and starch intake (r=0.56). Insulin plays a central role in metabolism by stimulating utilization of glucose in peripheral tissues such as muscle and adipose tissue and by promoting accumulation of glycogen and lipid reserves. Insulin was increased by diets with high starch content (Reynolds, 2006) and was decreased by diets with high fat content (Choi and Palmquist, 1996), although increased insulin concentrations were found when supplementary fat increased energy intake (Palmquist and Moser, 1981).

The curvilinear response of insulin to dietary starch concentration suggests that insulin secretion was reduced at low starch concentrations and reached a break point at approximately 160g of starch/kg of DM and 43g of fat/kg of DM. The high-merit cows used in our previous study (Gong et al., 2002) had insulin concentrations of 0.21ng/mL on the low-insulin diet and 0.32ng/mL on the high-insulin diet over the first 50 DIM; that difference was similar to the difference between equivalent diets (5 and 1) in the current experiment.

Glucagon plays a complementary role to insulin in metabolism; it stimulates mobilization of glycogen and lipid reserves, increases the rate of gluconeogenesis, and promotes a switch from the use of glucose to the use of fatty acids as a source of energy. Diets rich in carbohydrates decrease the concentration of glucagon in the blood, whereas fasting, or feeding low-carbohydrate diets, has the opposite effect. Although glucagon increased linearly with decreasing dietary starch concentration, the difference between diets did not reach significance until dietary starch concentration was below 160g/kg of DM.

Due to the opposing effects of insulin and glucagon on plasma glucose, insulin-to-glucagon ratio could be more important than absolute levels of each hormone. The ratio started to decrease when starch concentration decreased below 160g/kg of DM. Low insulin-to-glucagon ratio could indicate compromised gluconeogenesis from propionate and amino acids and may explain some of the reproductive effects seen with low-insulin diets (Gong et al., 2002).

In the current experiment, insulin was positively related to the maximum number of small follicles in follicular waves before and after ovulation, in agreement with the observations of Gutierrez et al. (1997). Greater numbers of small follicles were associated with lower FSH concentrations, independently of circulating estradiol, inhibin, or IGF-I concentrations (Burns et al., 2005). Declining FSH, accompanied by increasing LH pulses, aids the differentiation and maturation of dominant follicles during early lactation, thereby increasing the chance of ovulation in response to an LH surge (Webb et al., 2004). The increased number of medium-sized follicles for diet 3, compared with diets 1 and 2, suggests that increasing dietary fat concentration stimulated growth of follicles. Importantly, insulin effects in the current experiment were independent of energy balance, whereas in other studies, low insulin accompanied NEB (Butler, 2003).

As with insulin, nutritionally induced changes in systemic IGF-I concentrations have been linked to ovarian activity (Webb et al., 2004). In the current experiment, IGF-I was negatively related to DMI but did not vary among diets and was not related to plasma insulin or GH concentrations. The endocrine role for IGF-I in controlling ovarian activity cannot be viewed in isolation and must be combined with factors affecting its expression in the liver along with nutritionally induced changes in both the serum IGF binding protein profile and IGF binding protein expression in the dominant follicle (Webb et al., 2004). The results of this experiment provide further support for the suggestion that GH does not appear to be directly involved in the physiological mechanisms underlying the nutritional influence on ovarian function in cattle (Webb et al., 2004), although administration of exogenous GH does influence follicular populations (Gong et al., 1991).

Leptin appears to act as a further signal linking nutritional status with reproductive performance (Keisler et al., 1999) and is related to BCS in lactating cows and level of feeding in nonlactating cows (Armstrong et al., 2003). It appears that both insulin and IGF-I are involved in regulating leptin responses to nutrition. In the current study, the lack of variation in leptin among diets is consistent with the lack of variation in BCS and DMI, although the relationship between leptin and DMI (r=−0.44) across individual cows is consistent with the role of leptin in regulation of feed intake.

Mean progesterone concentration during the post-ovulatory increase was 40% lower in the greatest starch (lowest fat) group than in the other 4 groups. During this time, progesterone plays a critical role in supporting uterine function, and numerous studies have linked low progesterone secretion at this time to reduced embryo development and early embryo loss (Mann and Lamming, 2001). Interestingly, cows on diet 2 numerically had the greatest mean progesterone concentration and also the largest mean corpus luteum diameter, both of which were significantly different from diet 1, suggesting that a critical threshold for fatty acid supply had been met. Thus, high dietary starch (>183g/kg of DM) and low dietary fat (<42g/kg of DM) concentrations might be detrimental to pregnancy rate by reducing progesterone. Further support for this proposition comes from studies of developmental competence of oocytes. Fouladi-Nashta et al. (2005) found that cows fed a high-starch diet (equivalent to diet 2) produced oocytes that had a significantly lower proportion developing to the blastocyst stage after in vitro fertilization than cows fed on a high-fiber diet (equivalent to diet 5). In a later study (Fouladi-Nashta et al., 2007), a high-fat diet produced significantly greater blastocyst yields than a low-fat diet.

Diet composition had no effect on plasma metabolites except BHBA, NEFA. and urea-N. There was a linear increase in BHBA with decreasing dietary starch concentration. Grohn et al. (1983) suggested BHBA levels of <1mmol/L represented healthy cows, whereas levels of 1 to 3mmol/L were mildly ketotic cows. Thus, cows on diets 4 and 5 might have been mildly ketotic, but cows on diets 1 to 3 were healthy.

Elevated plasma NEFA concentrations were associated with NEB in early lactation and were linked to reduced reproductive performance (Butler, 2003; Garnsworthy et al., 2008). There was no difference among treatments in energy balance, which was generally positive; therefore, the positive relationship between NEFA and dietary fat concentration suggests that the greater NEFA values seen for diets 4 and 5 resulted from the low insulin-to-glucagon ratios observed in these groups.

The positive relationship between urea-N and dietary starch concentration is unusual; normally high-starch diets improve rumen ammonia capture and decrease plasma urea-N concentrations (Reynolds, 2006). The inverse effect could be due to the use of molassed sugar beet pulp as a fiber source, which reduced plasma urea concentrations (Kenny et al., 2002). The high moisture absorbency of beet pulp, approximately 20× its weight in water, could trap ammonia in solution within the rumen, improving its utilization for microbial protein synthesis. The increases in corpus luteum diameter and progesterone concentration for diet 2 compared with diet 1 could be linked to increased intakes of starch and sugars, or reduced plasma urea-N, the only variables that were significantly different between these diets. Yet, plasma urea-N concentrations were below the standard range (3.3 to 5mmol/L) reported by Ward et al. (1995) in all but 9 cows. The implications of low plasma urea-N concentrations are unknown. High plasma ammonia and urea levels were associated with elevated progesterone concentrations, reduced dominant follicle growth, and reduced cleavage rate of oocytes from medium-sized follicles (Sinclair et al., 2000), but whether the relative changes in urea-N found in the current study are important in this respect remains to be tested.

Back to Article Outline

Conclusions 

The results of this experiment demonstrate that insulin and glucagon can be changed by diet in modern high-yielding Holstein dairy cows and provide quantitative data that can be used to determine critical thresholds. The critical threshold for maintaining adequate insulin-to-glucagon ratios in cows at the start of the breeding period appears above 160g of starch/kg of DM and below 44g of fat/kg of DM. Although there is less direct evidence from this experiment, associated studies suggest that the greatest starch diet (231g of starch/kg of DM) could cause problems from low progesterone or impaired developmental competence of oocytes. Further evidence is provided to support the role of insulin as a metabolic signal influencing ovarian function, because nutritionally induced changes in insulin were associated with changes in follicle dynamics.

Back to Article Outline

Acknowledgments 

This study was part of Project LK0646 in the LINK Sustainable Livestock Production program, which was funded by The Scottish Executive Environment and Rural Affairs Department (SEERAD, Edinburgh, UK), ABNA Ltd. (Peterborough, UK), BOCM PAULS Ltd. (Ipswich, UK), and Provimi Ltd. (Sint-Stevens-Woluwe, Belgium).The following people served on the Programme Management Committee and made valuable contributions to the design and interpretation of the experiment: J. Newbold (Provimi), M. Marsden (ABNA), S. Richards (Provimi), A. Henderson (BOCM PAULS), P. Thomas (Artilus Ltd., UK), D. Armstrong (Roslin Institute, Edinburg, UK), A. Flint (Univ. Nottingham, UK), D. Garwes (Department for Environment, Food and Rural Affairs), and D. Leaver (Royal Agricultural College). We would like to thank the following people for technical assistance: H. Russell (Univ. Nottingham), J. Gong (Roslin Inst.), G. Baxter (Roslin Inst.), M. Mitchell (Univ. Nottingham), N. Armstrong (Univ. Nottingham), D. Scholey (Univ. Nottingham), M. Hunter (Univ. Nottingham), D. Whitaker (Univ. Edinburgh), and C. Smith (National Milk Records). Statistical advice was provided by J. Craigon (Univ. Nottingham).

Back to Article Outline

Supplementary data 

Interpretive summary.

Back to Article Outline

References 

  1. Adamiak SJ, Mackie K, Watt RG, Webb R, Sinclair KD. Impact of nutrition on oocyte quality: Cumulative effects of body composition and diet leading to hyperinsulinemia in cattle. Biol. Reprod. 2005;73:918–926
  2. Allison RD, Garnsworthy PC. Increasing the digestible undegraded protein intake of lactating dairy cows by feeding fishmeal or a rumen protected vegetable protein blend. Anim. Feed Sci. Technol. 2002;96:69–81
  3. AOAC. Official Methods of Analysis. 15th ed.. Arlington, VA: Association of Official Analytical Chemists; 1990;
  4. Armstrong DG, Gong JG, Webb R. Interactions between nutrition and ovarian activity in cattle: Physiological, cellular and molecular mechanisms. Reprod. Suppl. 2003;61:1–12
  5. Burns DS, Jimenez-Krassel F, Ireland JLH, Knight PG, Ireland JJ. Numbers of antral follicles during follicular waves in cattle: Evidence for high variation among animals, very high repeatability in individuals, and an inverse association with serum follicle-stimulating hormone concentrations. Biol. Reprod. 2005;73:54–62
  6. Butler WR. Energy balance relationships with follicular development, ovulation and fertility in postpartum dairy cows. Livest. Prod. Sci. 2003;83:211–218
  7. Choi B, Palmquist DL. High fat diets increase plasma cholecystokinin and pancreatic polypeptide, and decrease plasma insulin and feed intake in lactating cows. J. Nutr. 1996;126:2913–2919
  8. Fouladi-Nashta AA, Gutierrez CG, Garnsworthy PC, Webb R. Effects of dietary carbohydrate source on oocyte/embryo quality and development in high-yielding, lactating dairy cattle. Biol. Reprod. Special Issue. 2005;2005:135–136
  9. Fouladi-Nashta AA, Gutierrez CG, Gong JG, Garnsworthy PC, Webb R. Impact of dietary fatty acids on oocyte quality and development in lactating dairy cows. Biol. Reprod. 2007;77:9–17
  10. FSR. 2005. The Feeding Stuffs (England) Regulations 2005. Statutory Instrument No. 3281. Her Majesty's Stationery Office, London, UK.
  11. Garnsworthy PC. Body condition score in dairy cows: Targets for production and fertility. In:  Garnsworthy PC,  Wiseman J editor. Recent Advances in Animal Nutrition – 2006. Nottingham, UK: Nottingham University Press; 2007;p. 61–86
  12. Garnsworthy PC, Sinclair KD, Webb R. Integration of physiological mechanisms that influence fertility in dairy cows. Animal. 2008;2:1144–1152
  13. Garnsworthy PC, Topps JH. The effect of body condition of dairy cows at calving on their food intake and performance when given complete diets. Anim. Prod. 1982;35:113–119
  14. Gong JG, Bramley TA, Webb R. The effect of recombinant bovine somatotropin on ovarian function in heifers: Follicular populations and peripheral hormones. Biol. Reprod. 1991;45:941–949
  15. Gong JG, Lee WJ, Garnsworthy PC, Webb R. Effect of dietary-induced increases in circulating insulin concentrations during the early postpartum period on reproductive function in dairy cows. Reproduction. 2002;123:419–427
  16. Grohn Y, Lindberg LA, Bruss ML, Farver TB. Fatty infiltration of the liver in spontaneously ketotoic dairy cows. J. Dairy Sci. 1983;66:2320–2328
  17. Gutierrez CG, Gong JG, Bramley TA, Webb R. Selection on predicted breeding value for milk production delays ovulation independently of changes in follicular development, milk production and body weight. Anim. Reprod. Sci. 2006;95:193–205
  18. Gutierrez CG, Oldham J, Bramley TA, Gong JG, Campbell BK, Webb R. The recruitment of ovarian follicles is enhanced by increased dietary intake in heifers. J. Anim. Sci. 1997;75:1876–1884
  19. Keisler DH, Daniel JA, Morrison CD. The role of leptin in nutritional status and reproductive function. J. Reprod. Fertil. Suppl. 1999;54:425–435
  20. Kenny DA, Boland MP, Diskin MG, Sreenan JM. Effect of rumen degradable protein with or without fermentable carbohydrate supplementation on blood metabolites and embryo survival in cattle. Anim. Sci. 2002;74:529–537
  21. Mann GE, Lamming GE. Relationship between the maternal endocrine environment, early embryo development and the inhibition of the luteolytic mechanism in the cow. Reproduction. 2001;121:175–180
  22. Mann GE, Lamming GE, Fray MD. Plasma oestradiol and progesterone during early pregnancy in the cow and the effects of treatment with buserelin. Anim. Reprod. Sci. 1995;37:121–131
  23. Mulvanny, P. M. 1977. Dairy cow condition scoring. Paper 4468. National Institute for Research in Dairying, Reading, UK.
  24. Palmquist DL, Moser EA. Dietary fat effects on blood insulin, glucose utilization, and milk protein content of lactating cows. J. Dairy Sci. 1981;64:1664–1670
  25. Pryce JE, Royal MD, Garnsworthy PC, Mao IL. Fertility in the high producing dairy cow. Livest. Prod. Sci. 2004;86:125–135
  26. Reynolds CK. Production and metabolic effects of site of starch digestion in dairy cattle. Anim. Feed Sci. Technol. 2006;130:78–94
  27. Sinclair KD, Kuran M, Gebbie FE, Webb R, McEvoy TG. Nitrogen metabolism and fertility in cattle: II. Development of oocytes recovered from heifers offered diets differing in their rate of nitrogen release in the rumen. J. Anim. Sci. 2000;78:2670–2689
  28. Thomas C. Feed Into Milk: A New Applied Feeding System for Dairy Cows. Nottingham, UK: Nottingham University Press; 2004;
  29. Ward WR, Murray RD, White AR, Rees EM. The use of blood biochemistry for determining the nutritional status of dairy cows. In:  Garnsworthy PC,  Wiseman J editor. Recent Advances in Animal Nutrition – 1995. Nottingham, UK: Nottingham University Press; 1995;p. 29–51
  30. Webb R, Garnsworthy PC, Gong JG, Armstrong DG. Control of follicular growth: Local interactions and nutritional influences. J. Anim. Sci. 2004;82:E63–E74

PII: S0022-0302(08)71007-5

doi:10.3168/jds.2008-1031

Journal of Dairy Science
Volume 91, Issue 10 , Pages 3814-3823, October 2008