Advertisement
Research| Volume 102, ISSUE 10, P9468-9480, October 2019

Reducing milking frequency from thrice to twice daily in early lactation improves the metabolic status of high-yielding dairy cows with only minor effects on yields

  • U. Moallem
    Correspondence
    Corresponding author
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel
    Search for articles by this author
  • H. Kamer
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel

    Department of Animal Science, the Robert H. Smith Faculty of Agriculture, Food and Environment, the Hebrew University of Jerusalem, Rehovot 76100, Israel
    Search for articles by this author
  • A. Hod
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel

    Department of Animal Science, the Robert H. Smith Faculty of Agriculture, Food and Environment, the Hebrew University of Jerusalem, Rehovot 76100, Israel
    Search for articles by this author
  • L. Lifshitz
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel
    Search for articles by this author
  • G. Kra
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel
    Search for articles by this author
  • S. Jacoby
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel
    Search for articles by this author
  • Y. Portnick
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel
    Search for articles by this author
  • M. Zachut
    Affiliations
    Department of Ruminant Science, Institute of Animal Sciences, The Volcani Center, 68 HaMaccabim Road, Rishon LeZion 7505101, Israel
    Search for articles by this author
Open ArchivePublished:August 07, 2019DOI:https://doi.org/10.3168/jds.2019-16674

      ABSTRACT

      Reducing milk production during early lactation might be of interest to improve the energy balance (EB) of high-yielding dairy cows. Therefore, the objective of this study was to determine how reducing the milking frequency (MF) of high-yielding dairy cows from thrice to twice a day during the first 30 d in milk (DIM) affects yields, intake, efficiency, metabolic status, and carryover effects. To this end, 42 multiparous cows were divided into 2 groups according to their previous lactation performance, parity, and body weight. The control cows were milked 3 times a day (3ML) and the treated cows were milked twice a day (2ML) until 30 DIM and then both groups were milked 3 times a day. Milk samples were taken twice a week from 2 or 3 consecutive milkings until 45 DIM for analysis of milk solids, and both groups were followed until 100 DIM to determine the carryover effects of MF until 30 DIM. Individual dry matter intake (DMI), milk yield, and body weight were recorded daily. Blood samples were taken 3 times weekly from 14 d prepartum until 45 DIM. Milk yield during the first 30 DIM was 8.6% higher (49.3 and 45.4 kg/d, respectively), milk fat percentage was lower (3.96 and 4.27%, respectively), and the yields of all milk solids were higher in the 3ML cows than in the 2ML cows. Dry matter intake and 4% fat-corrected milk were similar between groups. The EB during the first 30 DIM was lower in the 3ML cows than in the 2ML cows, and milk yield, but not 4% fat-corrected milk yield, per unit of DMI was higher in the 3ML cows. No differences were observed between groups from 31 to 100 DIM in milk yield (∼56.3 kg/d for both groups), milk solids yield, DMI, or milk/DMI; however, fat percentage was lower and EB was higher in the 3ML cows. Blood glucose concentrations between 0 and 30 DIM were lower and β-hydroxybutyrate concentrations were higher in the 3ML cows than in the 2ML cows, but nonesterified fatty acids concentrations were lower, which may be attributed to the lower clearance frequency of nonesterified fatty acids from the blood stream in the 2ML cows. A lower proportion of the 3ML cows (10%) ovulated ≤15 DIM compared with the 2ML cows (40%), with no beneficial effects on preovulatory follicle characteristics. Reducing the MF from thrice to twice a day during the first 30 DIM improved EB and metabolic status, with only minor effects on production.

      Key words

      INTRODUCTION

      Negative energy balance (NEB), which occurs in early lactation in high-yielding dairy cows, is associated with an increase in the incidence of metabolic disorders (
      • Grummer R.R.
      Etiology of lipid-related metabolic disorders in periparturient dairy cows.
      ;
      • Esposito G.
      • Irons C.P.
      • Webb E.C.
      • Chapwanya A.
      Interactions between negative energy balance, metabolic diseases, uterine health and immune response in transition dairy cows.
      ) and reduced fertility (
      • Lucy M.C.
      Reproductive loss in high-producing dairy cattle: Where will it end?.
      ;
      • Butler W.R.
      Energy balance relationships with follicular development, ovulation and fertility in postpartum dairy cows.
      ). Hence, improving energy balance (EB) during early lactation is of great interest. There are 3 main approaches to improving EB: (1) increasing energy input, which includes enhancing energy consumption and, to some extent, improving the metabolism; (2) decreasing energy output, which means reducing the total energy produced in milk (or body mass); and (3) optimizing body tissue reserves at calving (BCS) to avoid massive body tissue mobilization postpartum (PP;
      • Morrow D.A.
      Fat cow syndrome.
      ). Several nutritional manipulations to alleviate NEB during early lactation by increasing energy input have been tested with mixed results. The most common nutritional strategy that has been intensely investigated is to increase the fat content of the cow's fresh diet to enhance energy density and possibly energy consumption. However, this strategy was not effective in reducing adipose body mobilization and the concentrations of plasma nonesterified fatty acids (NEFA;
      • Chilliard Y.
      Dietary fat and adipose tissue metabolism in ruminants, pigs, and rodents: A review.
      ;
      • Beam S.W.
      • Butler W.R.
      Energy balance, metabolic hormones, and early postpartum follicular development in dairy cows fed prilled lipid.
      ;
      • Bertics S.J.
      • Grummer R.R.
      Effects of fat and methionine hydroxy analog on prevention or alleviation of fatty liver induced by feed restriction.
      ). This was mainly because fat supplementation stimulated milk production (
      • Sklan D.
      • Moallem U.
      • Folman Y.
      Effects of feeding calcium soaps of fatty acids on production and reproduction responses in high production lactating cows.
      ;
      • Chilliard Y.
      Dietary fat and adipose tissue metabolism in ruminants, pigs, and rodents: A review.
      ) or sometimes decreased DMI (
      • Chilliard Y.
      Dietary fat and adipose tissue metabolism in ruminants, pigs, and rodents: A review.
      ;
      • Beam S.W.
      • Butler W.R.
      Energy balance, metabolic hormones, and early postpartum follicular development in dairy cows fed prilled lipid.
      ;
      • Moallem U.
      • Katz M.M.
      • Arieli A.
      • Lehrer H.
      Effects of peripartum propylene glycol or fats differing in fatty acid profiles on feed intake, production, and plasma metabolites in dairy cows.
      ) and consequently was not effective in improving EB. In a review paper,
      • Leroy J.L.
      • Sturmey R.G.
      • Van Hoeck V.
      • De. Bie J.
      • McKeegan P.J.
      • Bols P.E.
      Dietary fat supplementation and the consequences for oocyte and embryo quality: Hype or significant benefit for dairy cow reproduction?.
      concluded that dietary fat supplementation has no added value in alleviating NEB during the transition period in high-producing dairy cows. Another strategy tested was to increase the glucogenic precursor's consumption during the transition period. In a study by
      • Moallem U.
      • Katz M.M.
      • Arieli A.
      • Lehrer H.
      Effects of peripartum propylene glycol or fats differing in fatty acid profiles on feed intake, production, and plasma metabolites in dairy cows.
      , cows were supplemented with propylene glycol from 21 d prepartum through 21 DIM, but no improvement in EB was achieved. In a study by
      • Butler S.T.
      • Pelton S.H.
      • Butler W.R.
      Energy balance, metabolic status, and the first PP ovarian follicle wave in cows administered propylene glycol.
      , cows were drenched daily with propylene glycol from d 10 prepartum until d 25 PP; the metabolic status of the cows improved but with no beneficial effect on the reproductive system. Other efforts attempted to alleviate EB during the transition period, such as reducing milk energy output by supplementing the isomers of CLA, which are potent inhibitors of mammary lipid synthesis; however, this approach achieved very limited results (
      • Bernal-Santos G.
      • Perfield II, J.W.
      • Barbano D.M.
      • Bauman D.E.
      • Overton T.R.
      Production responses of dairy cows to dietary supplementation with conjugated linoleic acid (CLA) during the transition period and early lactation.
      ;
      • Shingfield K.J.
      • Beever D.E.
      • Reynolds C.K.
      • Gulati S.K.
      • Humphries D.J.
      • Lupoli B.
      • Hervás G.
      • Griinari J.M.
      Effect of rumen protected conjugated linoleic acid on energy metabolism of dairy cows during early to mid-lactation.
      ;
      • Castañeda-Gutiérrez E.
      • Overton T.R.
      • Butler W.R.
      • Bauman D.E.
      Dietary supplements of two doses of calcium salts of conjugated linoleic acid during the transition period and early lactation.
      ).
      The relative failure to improve EB through nutritional manipulation led to efforts to develop other management strategies to achieve this goal. Shortening the dry period has been proposed as a management strategy to improve the EB of dairy cows in early lactation (
      • Grummer R.R.
      • Rastani R.R.
      Why reevaluate dry period length?.
      ). Cows with continuous milking or a short dry period produced 12 to 25% less milk in the subsequent lactation; however, these cows also had less severe negative EB PP compared with cows with a traditional dry period length (∼60 d;
      • Rastani R.R.
      • Grummer R.R.
      • Bertics S.J.
      • Gümen A.
      • Wiltbank M.C.
      • Mashek D.G.
      • Schwab M.C.
      Reducing dry period length to simplify feeding transition cows: Milk production, energy balance, and metabolic profiles.
      ). This was also demonstrated by reduced PP plasma NEFA and liver triglyceride (TG) concentrations (
      • Rastani R.R.
      • Grummer R.R.
      • Bertics S.J.
      • Gümen A.
      • Wiltbank M.C.
      • Mashek D.G.
      • Schwab M.C.
      Reducing dry period length to simplify feeding transition cows: Milk production, energy balance, and metabolic profiles.
      ), with a beneficial effect on the reproductive system (
      • Gümen A.
      • Rastani R.R.
      • Grummer R.R.
      • Wiltbank M.C.
      Reduced dry periods and varying prepartum diets alter postpartum ovulation and reproductive measures.
      ). The relatively great loss in milk yield and other relevant aspects, such as the effect of continuous milking or a short dry period on animal welfare, preclude adopting this management strategy as a viable practice in the field.
      Increasing milking frequency (MF) is well established as a management tool to increase milk production in dairy cows, which has a central effect on EB.
      • Erdman R.A.
      • Varner M.
      Fixed yield responses to increased milking frequency.
      found a fixed milk yield response, with an increase of 3.5 kg/d, when MF was increased from twice to thrice a day regardless of the yield for cows milked twice daily. Dairy cows in Israel are commonly milked 3 times a day. We therefore hypothesized that in high-yielding dairy cows the magnitude of the decrease in milk yield when reducing MF from thrice to twice a day in early lactation would be small but might be effective in alleviating EB. Hence, the objective of the present study was to examine how reducing the MF of multiparous high-yielding dairy cows from thrice to twice daily during the first 30 DIM affected milk and milk solids yields, feed intake, efficiency, EB, and metabolic status as well as the reproductive system. In addition, we followed the cows until 100 DIM to determine any potential carryover effects of reduced MF during the early PP period.

      MATERIALS AND METHODS

      The experimental protocol and procedures were approved by the Volcani Center Animal Care Committee (IL 637/16). The study was conducted at the experimental dairy farm of the Volcani Center in Rishon LeZion, Israel. Forty-two high-yielding multiparous Israeli-Holstein pregnant cows were stratified into 2 treatment groups at 256 d of pregnancy, blocked on milk production and fat content during the first 120 d of the previous lactation, parity, BW at drying-off, and the expected time of the approaching parturition. The cows were group housed at 256 d of pregnancy in a shaded loose pen that was equipped with a real-time electronic individual feeding system and fed a dry-cow diet. Each feeding station included an individual identification system (ID tag; S.A.E. Afikim, Kibbutz Afikim, Israel) that allowed each cow to enter a specific feeding station and automatically recorded each meal. The calving period continued for 2 mo. Immediately after calving, the cows were moved to the side of the same barn housing the lactating cows. The cows were divided into 2 subgroups: the control cows were milked thrice a day (3ML; at 0050, 1300, and 2000 h) and the treatment cows were milked twice a day (2ML; at 0700 and 1900 h) until 30 DIM. From 30 DIM, the 2ML cows joined the 3ML cow group, and all cows were then milked 3 times a day. The cows were followed until 100 DIM to determine the carryover effects. Milk samples were taken twice a week (on Monday and Thursday) from d 5 after calving until 30 DIM from 2 consecutive milkings for the 2ML cows or 3 consecutive milkings for the 3ML cows. From 30 to 45 DIM, milk samples were taken twice a week and then monthly from 3 consecutive milkings until 100 DIM. The milk samples were analyzed for milk fat, protein, lactose, and urea by infrared analysis (standard IDF 141C:2000;

      IDF (International Dairy Federation). 2000. IDF Standard 141c:2000 (Determination of milk fat, protein and lactose content—Guidance on the operation of mid-infrared instruments. International Dairy Federation, Brussels, Belgium.).

      ) at the laboratories of the Israeli Cattle Breeders' Association (Caesarea, Israel). Somatic cell counts were also determined in the same laboratory.
      Cows were automatically weighed 3 times daily, after each milking, with a walk-on electronic scale (S.A.E. Afikim). We also determined the average BW at wk 1, 5, and 15 for each cow using daily records and then calculated the individual gain between these time points. The cows were also equipped with collar-mounted tags (HR-Tags; SCR Engineers Ltd., Hadarim, Netanya, Israel) that monitored and transmitted their rumination time. Rumination data were recorded automatically and stored in 2-h blocks. The cows were also equipped with a PlusTM sensor (S.A.E. Afikim), which monitored their activity (steps/h) and lying time (min). The BCS of cows (scale of 1 to 5;
      • Edmonson A. J
      • Lean I.J.
      • Weaver L.D.
      • Farver T.
      • Webster G.
      A body condition scoring chart for Holstein dairy cows.
      ) was determined weekly by 1 technician from 21 d prepartum to 100 PP.
      Postpartum, the cows were fed a standard Israeli milking cow ration. The composition and content of the diet are presented in Table 1. The diet was offered once daily at 1000 h ad libitum to about 5% orts. The individual amounts offered and the daily leftovers were recorded daily to calculate individual feed intake. The diet was sampled weekly and dried at 60°C for 48 h for DM determination.
      Table 1Ingredients and chemical composition of the diet
      ItemValue
      Ingredient, % of DM
       Corn, ground20.2
       Barley, rolled3.3
       Wheat grain, rolled4.8
       Rapeseed meal7.1
       Whole cottonseed2.1
       Gluten feed14.8
       Wheat silage10.8
       Corn silage14.2
       Oat hay9.0
       Dried distillers grains8.9
       Mother liquor
      Dairy industry byproduct.
      0.6
       Calcium salt of fatty acids1.0
       Urea0.3
       Bicarbonate1.0
       Limestone0.4
       NaCl1.4
       Vitamins and minerals
      Containing (per kg of DM) 20,000,000 IU of vitamin A; 2,000,000 IU of vitamin D; 15,000 IU of vitamin E; 6,000 mg of Mn; 6,000 mg of Zn; 2,000 mg of Fe; 1,500 mg of Cu; 120 mg of I; 50 mg of Se; 20 mg of Co.
      0.1
      Chemical composition
       NEL,
      Based on NRC (1989) values of individual ingredients.
      Mcal/kg of DM
      1.78
       CP, %16.5
       NDF, %32.2
       Forage NDF, %17.5
       Total fat (ether extract), %4.9
       Ca, %0.009
       P, %0.005
      1 Dairy industry byproduct.
      2 Containing (per kg of DM) 20,000,000 IU of vitamin A; 2,000,000 IU of vitamin D; 15,000 IU of vitamin E; 6,000 mg of Mn; 6,000 mg of Zn; 2,000 mg of Fe; 1,500 mg of Cu; 120 mg of I; 50 mg of Se; 20 mg of Co.
      3 Based on
      • National Research Council
      Nutrient Requirements of Dairy Cattle.
      values of individual ingredients.
      Blood samples were taken 3 times weekly (Sunday, Tuesday, and Thursday) from 14 d before the expected calving date until 45 DIM. The blood samples were collected by coccygeal venipuncture after the morning milking into vacuum tubes containing lithium heparin (Becton Dickinson, Cowley, UK), and tubes were immediately placed in ice. Plasma was separated by centrifugation at 1,000 × g for 15 min at 4°C, divided into 4 tubes, and stored at −80°C pending analysis. Concentrations of glucose, NEFA, BHB, TG, and aspartate aminotransferase (AST) were determined in blood prepartum and PP, and progesterone blood concentrations were determined from 7 to 45 DIM to determine the cyclicity resumption PP.

      Energy Calculations

      Energy-corrected milk and EB were calculated according to
      • National Research Council
      Nutrient Requirements of Dairy Cattle.
      as follows:
      NEC (Mcal/d) = NEL/kg of DM × DMI


      NEM (Mcal/d) = BW0.75 × 0.08


      ECM (Mcal/d) = milk (kg) × [(0.0929 × fat %) + (0.0547 × protein %) + (0.0395 × lactose %)]


      EB (Mcal/d) = NEC − (NEM + ECM),


      where NEC = net energy consumed and ECM = net energy output in milk. The NEL values for feedstuffs used in the formulated diets were calculated using
      • National Research Council
      Nutrient Requirements of Dairy Cattle.
      values with a few adjustments.

      Ovarian Ultrasonography

      From d 9 to 11 PP, the ovarian follicular activity of 8 cows from each group was examined twice a week using a MylbFive linear array ultrasonography transducer (Pie Medical, Maastricht, the Netherlands) until the day of the first ovulation, followed by the corpus luteum existence, or 34 DIM. At each monitoring, the follicles on both ovaries were counted according to the following criteria: 3- to 5-mm follicles (class 1), 6- to 9-mm follicles (class 2), and ≥10-mm follicles (class 3). Follicular cysts were defined as nonovulatory follicles that were ≥25 mm in diameter and that persisted for ≥10 d.
      In addition, blood samples from all cows from 7 d PP and later were analyzed for progesterone concentrations to determine the time of the first ovulation PP. The first ovulations afterward were estimated to occur 3 d before the first PP increase in progesterone ≥1 ng/mL.

      Follicular Fluid Aspiration

      Follicular fluid (FF) aspiration from preovulatory follicles was performed between 60 and 80 DIM in 8 cows from each group that was randomly selected to assess the effect of a different MF on the characteristics of preovulatory follicles. Ovaries of cows were monitored with a Aquila 5-MHz linear array ultrasonography transducer (Pie Medical) at 50 DIM, and cows that were confirmed as cycling by the presence of a corpus luteum (CL) were treated intramuscularly with a PGF analog (cloprostenol sodium, Estroplan; Parnell Living Science, Alexandria, NSW, Australia) at a dosage of 500 µg/cow. The cows were examined for signs of estrus, and at 14 to 15 d after a visible behavioral estrus, they received another PGF injection followed by FF aspiration after 48 h. Cows that were not confirmed as cycling by CL presence at 50 DIM were treated with a GnRH analog consisting of 200 µg of gonadorelin acetate (Gonabreed; Parnell Living Science). Ten days after the GnRH injection, the cows were monitored with ultrasound. After the presence of CL was confirmed they were injected with a PGF analog, and 14 to 15 d after a visible behavioral estrus they received another PGF injection, followed by FF aspiration after 48 h.
      Forty-eight hours after the second PGF administration, we performed FF aspiration as described by
      • Moallem U.
      • Folman Y.
      • Bor A.
      • Arav A.
      • Sklan D.
      Effect of calcium soaps of fatty acids and administration of somatotropin on milk production, preovulatory follicular development, and plasma and follicular fluid lipid composition in high yielding dairy cows.
      . Briefly, the cows were sedated with an intramuscular administration of 1 mL of 2% Sedaxylan (xylazine base, 20 mg/mL; Eurovet Animal Health, AE Bladel, the Netherlands) and were given a local anesthesia in the form of 5 mL of 2% lidocaine HCl (200 mg/10 mL; Esracain 2%; Rafa Laboratories, Jerusalem, Israel) injected epidurally between the last sacral and the first caudal vertebrae in the sacrococcygeal space. Ovaries were monitored, and follicles ≥7 mm in diameter were aspirated separately. The FF was aspirated with the aid of an ultrasound scanner (Pie Medical) connected to a 7.5-MHz vaginal sector transducer equipped with a needle guide and an MP86 suction pump (Biometra, Goettingen, Germany) set at a flow rate of 25 to 30 mL/min. The needles used were 18 gauge and were changed between follicles. After collection, the FF was centrifuged at 3,000 × g for 15 min at 4°C, separating the sediment, including granulosa cells, which was separated from the fluids, and both fractions were frozen at −32°C pending analysis. Jugular venous blood samples were collected in lithium heparin-coated vacuum tubes (Becton Dickinson, Cowley, UK) on the day of the second PGF treatment and 48 h later (i.e., in parallel with the aspiration of the preovulatory follicles) to verify that the cows were in the expected stage of the estrous cycle.

      Chemical Analysis

      The dried diet samples were analyzed for their content of N (
      • AOAC
      Official Methods of Analysis.
      , method 984.13; Kjeltec 8400 Analyzer, Foss Analytical AB, Höganäs, Sweden), NDF, and ADF (equipment from Ankom Technology, Fairport, NY) using α-amylase and sodium sulfite for NDF according to
      • Van Soest P.J.
      • Robertson J.B.
      • Lewis B.A.
      Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition.
      ; Ca (
      • AOAC
      Official Methods of Analysis.
      , method 935.13); P (
      • AOAC
      Official Methods of Analysis.
      , method 964.06); and ether extracts (
      • AOAC
      Official Methods of Analysis.
      , method 996.06). Samples were dried at 550°C for 3 h for ash determination.
      • National Research Council
      Nutrient Requirements of Dairy Cattle.
      NEL values of feedstuffs were used.
      The concentrations of NEFA in plasma were determined with an NEFA C Test Kit (Wako Chemicals GmbH, Neuss, Germany). The intra- and interassay coefficients of variation for the NEFA assay were 5.9 and 6.1%, respectively.
      The plasma BHB concentration was determined with a Ranbut D-3-Hydroxybutyrate kit (Randox, Crumlin, UK) in which a reaction of 3-hydroxybutyrate and dehydrogenase generated a UV emission that was correlated with the sample's BHB concentration. The samples were examined at 340 nm with an optical density reader, and the results were calibrated against known BHB concentrations. The intra- and interassay coefficients of variation for the BHB assay were 1.3 and 1.6%, respectively. Glucose, TG, and AST concentrations in plasma were analyzed using the Cobas C111 Autoanalyser (Roche Holding GmbH, Grenzach-Wyhlen, Germany) using appropriate reagents.
      The FF of all aspirated preovulatory follicles (n = 40) was assessed for the relevant concentrations of progesterone (P4) and estrogen (E2; Progest-RIA and Ester-CTRIA, respectively; Cisbio Bioassays, Parc Marcel Boiteux, France) as well as androstenedione (A4; Immunotech s.r.o, Prague, Czech Republic). All FF samples were diluted 100-, 500-, or 30-fold for the P4, E2, and A4 assays, respectively, so that the assay results matched the detection ranges.
      The minimal detectable concentrations of P4, E2, and A4 were 0.2, 20, and 0.1 ng/mL, respectively. The intra-assay and interassay coefficients of variation for the P4, E2, and A4 assays were 9.1, 4.0, and 6.0% and 8.4, 3.4, and 4.3%, respectively. We considered follicles with E2:P4 ratios >1 to be E2 active (
      • Ireland J.J.
      • Roche F.
      Development of antral follicles in cattle after prostaglandin-induced luteolysis: Changes in serum hormones, steroids in follicles fluids and gonadotropin receptors.
      ) and subjected them to further analysis and evaluation. The P4 concentration in blood was determined using the same kit (Progest-RIA; Cisbio Bioassays) with a P4 detection threshold of 0.2 pg/mL.

      Statistical Analysis

      Production, EB, and efficiency calculations were analyzed as 3 periods—0 to 30 DIM (period 1), 31 to 100 DIM (period 2), and 0 to 100 DIM—for determining the effects of the whole period on production and efficiency. Continuous variables (milk, milk solids, DMI, and efficiency calculations) were analyzed as repeated measurements with PROC MIXED of SAS version 9.2 (
      • SAS Institute
      SAS User's Guide: Statistics. Version 9.2.
      ). When relevant, variables were analyzed using the specific data of the previous lactation data as covariates.
      The model used was as follows:
      Yijkl = µ + Ti + Lj + C(T × L)ijk + Pl + TPm + DIMijklmn + Eijklmno,


      where Yijkl = dependent variable; µ = overall mean; Ti = treatment effect (i = 1 to 2); Lj = parity (j = 2 or >2); C(T × L)ijk = cow k nested in treatment i and parity j; Pl = period effect; TPm = treatment × period; DIMijklmn = day with milk as a continuous variable; and Eijklmno = random residual. The treatment × DIM interaction was tested and found to be not significant and therefore was excluded from the model. For the 0 to 100 DIM analysis, the effects of period and treatment × period were removed from the model.
      Metabolite concentrations in blood were analyzed for 2 periods (0–30 DIM and 31–50 DIM) for determining the carryover effects using the above model. Body weight gain was analyzed using PROC GLM of SAS version 9.2 (
      • SAS Institute
      SAS User's Guide: Statistics. Version 9.2.
      ). Distributions of ovulation onset were analyzed using PROC FREQ of SAS version 9.2 (
      • SAS Institute
      SAS User's Guide: Statistics. Version 9.2.
      . Least squares means and adjusted SEM are presented in Tables 2, 3, 4, 5, and 6; P < 0.05 was accepted as significant unless otherwise stated.
      Table 2Milk and milk solids yields from 0 to 30 and 31 to 100 DIM for high-yielding dairy cows milked thrice daily (3ML) or twice daily (2ML) from calving until 30 DIM and then milked 3 times daily
      Item0–30 DIM31–100 DIMSEMP-value
      3ML2ML3ML2MLMF
      Milking frequency until 30 DIM.
      PeriodMF × period
      Milk, kg/d49.3
      Within rows, means with different superscript letters are statistically different.
      45.4
      Within rows, means with different superscript letters are statistically different.
      56.4
      Within rows, means with different superscript letters are statistically different.
      56.2
      Within rows, means with different superscript letters are statistically different.
      0.3>0.001>0.001>0.001
      Fat, %3.96
      Within rows, means with different superscript letters are statistically different.
      4.27
      Within rows, means with different superscript letters are statistically different.
      3.87
      Within rows, means with different superscript letters are statistically different.
      4.09
      Within rows, means with different superscript letters are statistically different.
      0.07>0.0010.170.39
      Protein, %3.283.273.203.190.320.750.070.91
      Lactose, %4.93
      Within rows, means with different superscript letters are statistically different.
      4.87
      Within rows, means with different superscript letters are statistically different.
      5.03
      Within rows, means with different superscript letters are statistically different.
      4.94
      Within rows, means with different superscript letters are statistically different.
      0.030.0070.040.64
      Fat, kg/d2.06
      Within rows, means with different superscript letters are statistically different.
      1.92
      Within rows, means with different superscript letters are statistically different.
      2.04
      Within rows, means with different superscript letters are statistically different.
      2.15
      Within rows, means with different superscript letters are statistically different.
      0.040.580.080.001
      Protein, kg/d1.72
      Within rows, means with different superscript letters are statistically different.
      1.49
      Within rows, means with different superscript letters are statistically different.
      1.72
      Within rows, means with different superscript letters are statistically different.
      1.68
      Within rows, means with different superscript letters are statistically different.
      0.02<0.0010.007<0.001
      Lactose, kg/d2.62
      Within rows, means with different superscript letters are statistically different.
      2.26
      Within rows, means with different superscript letters are statistically different.
      2.73
      Within rows, means with different superscript letters are statistically different.
      2.67
      Within rows, means with different superscript letters are statistically different.
      0.05<0.0010.0010.001
      4% FCM, kg/d46.5
      Within rows, means with different superscript letters are statistically different.
      45.8
      Within rows, means with different superscript letters are statistically different.
      52.6
      Within rows, means with different superscript letters are statistically different.
      55.3
      Within rows, means with different superscript letters are statistically different.
      0.40.002<0.001<0.001
      ECM, kg/d36.6
      Within rows, means with different superscript letters are statistically different.
      34.3
      Within rows, means with different superscript letters are statistically different.
      40.0
      Within rows, means with different superscript letters are statistically different.
      40.8
      Within rows, means with different superscript letters are statistically different.
      0.30.001<0.001<0.001
      MUN, mg/dL15.8
      Within rows, means with different superscript letters are statistically different.
      17.2
      Within rows, means with different superscript letters are statistically different.
      17.0
      Within rows, means with different superscript letters are statistically different.
      16.4
      Within rows, means with different superscript letters are statistically different.
      0.50.340.800.02
      a–d Within rows, means with different superscript letters are statistically different.
      1 Milking frequency until 30 DIM.
      Table 3Dry matter intake, energy balance, and efficiency calculations from 0 to 30 and 31 to 100 DIM for high-yielding dairy cows milked thrice daily (3ML) or twice daily (2ML) from calving until 30 DIM and then milked 3 times daily
      Item0–30 DIM31–100 DIMSEMP-value
      3ML2ML3ML2MLMF
      Milking frequency until 30 DIM.
      PeriodMF × period
      DMI, kg/d27.827.931.531.00.20.22>0.0010.07
      Energy balance, Mcal/d−0.7
      Within rows, means with different superscript letters are statistically different.
      1.3
      Within rows, means with different superscript letters are statistically different.
      7.2
      Within rows, means with different superscript letters are statistically different.
      4.8
      Within rows, means with different superscript letters are statistically different.
      0.30.46<0.001<0.001
      Milk/DMI,
      Milk yield (kg/d) divided by DMI (kg/d).
      kg/kg
      1.83
      Within rows, means with different superscript letters are statistically different.
      1.68
      Within rows, means with different superscript letters are statistically different.
      1.83
      Within rows, means with different superscript letters are statistically different.
      1.85
      Within rows, means with different superscript letters are statistically different.
      0.01<0.001<0.001<0.001
      4% FCM/DMI,
      4% FCM (kg/d) divided by DMI (kg/d).
      kg/kg
      1.71
      Within rows, means with different superscript letters are statistically different.
      1.70
      Within rows, means with different superscript letters are statistically different.
      1.70
      Within rows, means with different superscript letters are statistically different.
      1.81
      Within rows, means with different superscript letters are statistically different.
      0.010.0010.005<0.001
      ECM/DMI,
      4% ECM (Mcal/d) divided by DMI (kg/d).
      Mcal/kg
      1.36
      Within rows, means with different superscript letters are statistically different.
      1.29
      Within rows, means with different superscript letters are statistically different.
      1.29
      Within rows, means with different superscript letters are statistically different.
      1.35
      Within rows, means with different superscript letters are statistically different.
      0.010.850.98<0.001
      Rumination time, min/d488.7
      Within rows, means with different superscript letters are statistically different.
      479.7
      Within rows, means with different superscript letters are statistically different.
      519.7
      Within rows, means with different superscript letters are statistically different.
      529.7
      Within rows, means with different superscript letters are statistically different.
      4.80.90<0.0010.02
      Lying time, min/d612.1
      Within rows, means with different superscript letters are statistically different.
      628.2
      Within rows, means with different superscript letters are statistically different.
      575.5
      Within rows, means with different superscript letters are statistically different.
      614.3
      Within rows, means with different superscript letters are statistically different.
      6.0<0.0010.0050.03
      a–d Within rows, means with different superscript letters are statistically different.
      1 Milking frequency until 30 DIM.
      2 Milk yield (kg/d) divided by DMI (kg/d).
      3 4% FCM (kg/d) divided by DMI (kg/d).
      4 4% ECM (Mcal/d) divided by DMI (kg/d).
      Table 4Milk and milk solids yields and efficiency calculations from 0 to 100 DIM for high-yielding dairy cows milked thrice daily (3ML) or twice daily (2ML) from calving until 30 DIM and then milked 3 times daily
      ItemTreatmentSEMP-value
      3ML2ML
      Milk, kg/d51.550.60.80.44
      Fat, %3.914.210.060.001
      Protein, %3.303.290.050.96
      Lactose, %4.964.890.020.04
      Fat, kg/d2.042.020.030.64
      Protein, kg/d1.701.570.020.001
      Lactose, kg/d2.562.380.060.03
      4% FCM, kg/d47.949.70.980.20
      ECM, kg/d37.437.50.590.89
      DMI, kg/d29.829.90.30.92
      Energy balance, Mcal/d4.593.840.390.19
      BW gain from 0 to 100 DIM, kg−27.9−14.95.90.12
      Milk/DMI,
      Milk yield (kg/d) divided by DMI (kg/d).
      kg/kg
      1.811.770.020.20
      4% FCM/DMI,
      4% FCM (kg/d) divided by DMI (kg/d).
      kg/kg
      1.691.760.020.01
      ECM/DMI,
      4% ECM (Mcal/d) divided by DMI (kg/d).
      Mcal/kg
      1.301.330.010.22
      Rumination time, min/d501.6507.46.60.54
      Lying time, min/d578.0616.05.40.001
      1 Milk yield (kg/d) divided by DMI (kg/d).
      2 4% FCM (kg/d) divided by DMI (kg/d).
      3 4% ECM (Mcal/d) divided by DMI (kg/d).
      Table 5Blood metabolite concentrations from 0 to 30 and 31 to 50 DIM for high-yielding dairy cows milked thrice daily (3ML) or twice daily (2ML) from calving until 30 DIM and then milked 3 times daily
      Item
      NEFA = nonesterified fatty acids; TG = triglyceride; AST = aspartate aminotransferase.
      0–30 DIM31–50 DIM
      TreatmentSEMP-valueTreatmentSEMP-value
      3ML2ML3ML2ML
      Glucose, mg/dL58.663.80.900.00162.464.81.00.11
      BHB, mg/dL7.216.280.250.027.37.30.340.90
      NEFA, μEq/L327.0439.122.60.002248.3290.718.20.11
      TG, mg/dL7.588.00.160.128.68.40.210.52
      AST, U/L75.974.02.20.5368.871.82.20.35
      1 NEFA = nonesterified fatty acids; TG = triglyceride; AST = aspartate aminotransferase.
      Table 6Concentrations of progesterone, androstenedione, and estradiol in estrogen-active follicles aspirated from high-yielding dairy cows milked thrice daily (3ML) or twice daily (2ML) from calving until 30 DIM and then milked 3 times daily
      Preovulatory follicles were aspirated 48 h after prostaglandin injection from 8 cows from each group.
      ItemTreatmentSEMP-value
      3ML2ML
      Cows, no.88
      Follicles aspirated, no.2119
      Atretic follicles, no.811
      Estrogen-active follicles, no.139
      Diameter of aspirated follicles, mm12.911.40.880.82
      Volume of aspirated follicles, mL1.551.610.300.93
      Progesterone, ng/mL49.455.114.10.78
      Androstenedione, ng/mL55.970.517.50.56
      Estradiol, ng/mL1,509.41,948.5265.00.26
      Estradiol:progesterone ratio56.945.612.50.53
      Total estradiol content, μg2.573.781.10.46
      Total progesterone content, μg0.140.100.070.70
      Total androstenedione content, μg0.080.160.050.23
      1 Preovulatory follicles were aspirated 48 h after prostaglandin injection from 8 cows from each group.

      RESULTS

      Milk and milk solid yields are presented in Table 2. Milk yield was 8.6% greater from 0 to 30 DIM in the 3ML cows than in the 2ML cows (P < 0.001: Figure 1A) but was similar between groups from 31 to 100 DIM (P = 0.51). The effects of period and treatment × period for milk yield were significant (P < 0.001). Fat content was 0.31 percentage units lower from 0 to 30 DIM in the 3ML cows than in the 2ML cows (P = 0.001), whereas it was 0.22 percentage units lower from 31 to 100 DIM in the 3ML cows (P = 0.03). The effect of period and treatment × period for fat content was not significant. Protein content in milk was not different between groups in both periods, and the effect of period and treatment × period was not significant. Lactose percentage until 30 DIM tended to be higher in the 3ML cows (P = 0.07) and was significantly higher than in the 2ML cows from 31 to 100 DIM (P = 0.04). Fat yield (kg/d) until 30 DIM was higher in the 3ML cows (P = 0.003), but from 31 to 100 DIM fat yield tended to be lower in the 3ML cows (P = 0.07). The effect of period on fat yield tended to be higher (P = 0.08), and the effect of treatment × period was significantly higher (P = 0.001). Protein and lactose yields were higher in the 3ML cows from 0 to 30 DIM, but from 31 to 100 DIM they were similar.
      Figure thumbnail gr1
      Figure 1(A) Milk yield, (B) 4% FCM yield, and (C) ECM yield from calving until 100 DIM for high-yielding dairy cows milked 3 times daily (3ML; ▴) or twice daily (2ML; ○) until 30 DIM (vertical arrow) and then milked 3 times daily. Values are LSM, and error bars represent SE.
      The 4% FCM yield from 0 to 30 DIM was similar between groups but was 5.1% higher from 31 to 100 DIM in the 2ML cows (P < 0.001; Figure 1B). The effects of period and treatment × period for 4% FCM yield were significant (P < 0.001). The ECM yield from 0 to 30 DIM was higher in the 3ML group (P < 0.001; Figure 1C) but from 31 to 100 DIM was higher in the 2ML cows (P = 0.001). The MUN was 8.9% higher in the 2ML cows during the first 30 DIM (P = 0.002) but was similar between groups from 31 to 100 DIM.
      The DMI from 0 to 30 DIM was similar between groups (Table 3; Figure 2A) but was higher in the 3ML cows between 31 and 100 DIM (P = 0.009). The effect of period on DMI was significant (P < 0.001), and the treatment × period interaction tended to be higher (P = 0.07). The average EB was higher from 0 to 30 DIM in the 2ML cows than in the 3ML cows (P < 0.001; Figure 2B), but from 31 to 100 DIM the EB was higher in the 3ML cows than in the 2ML cows (P < 0.001). The effects of period and treatment × period on EB were significant (P < 0.001). No differences were observed between groups in BCS measurements (data not presented). The milk yield per unit of DMI (milk/DMI) until 30 DIM was higher in the 3ML cows (P < 0.001) and from 31 to 100 DIM it tended to be higher in the 3ML cows (P = 0.1). The effects of period and treatment × period interaction on milk/DMI ratio were significant (P < 0.001). The 4% FCM yield per unit of DMI (FCM/DMI) was not different between groups until 30 DIM, but from 31 to 100 DIM it was higher in the 2ML cows (P < 0.001). The effect of period and treatment × period interaction on FCM/DMI ratio was significant (P < 0.005). The ECM yield (Mcal/d) per unit of DMI (ECM/DMI) until 30 DIM was higher in the 3ML cows (P < 0.001), but from 31 to 100 DIM it was higher in the 2ML cows (P < 0.001). The effect of period on ECM/DMI was not significant, but the treatment × period interaction was significant (P < 0.001). A similar result was found for the conversion ratio between the energy intake and ECM (data not shown). The daily rumination time until 30 DIM was similar between groups, but from 31 to 100 DIM it was higher in the 2ML cows (P = 0.03). The effect of period on rumination time and the treatment × period interaction were significant (P = 0.02 and P < 0.001, respectively). The daily resting time until 30 DIM tended to be higher in the 2ML cows (P = 0.06), and from 31 to 100 DIM it was higher in the 2ML cows (P < 0.001). The effects of period and treatment × period on resting time were significant (P = 0.005 and P = 0.03, respectively).
      Figure thumbnail gr2
      Figure 2(A) Dry matter intake, (B) energy balance, and (C) BW from calving until 100 DIM for high-yielding dairy cows milked 3 times daily (3ML; ▴) or twice daily (2ML; ○) until 30 DIM (vertical arrow) and then milked 3 times daily. Values are LSM, and error bars represent SE.
      We also analyzed the main production and efficiency results from calving until 100 DIM to determine the overall effect on production (Table 4). Milk yield was not different between groups, but fat percentage was 0.3 percentage units higher in the 2ML cows than in the 3ML cows (P < 0.001). The protein percentage was not different and the lactose content was higher in the 3ML cows than in the 2ML cows (P = 0.04). Fat yield was not different, but protein and lactose yields were higher in the 3ML cows than in the 2ML cows. The 4% FCM and ECM yields as well as DMI and EB were not different between groups for 100 DIM. The BW loss from calving to 100 DIM was 27.9 kg in the 3ML cows compared with 14.9 kg in the 2ML cows (P = 0.12). The BW gain from 31 to 100 DIM was 17.4 kg in the 2ML cows compared with 6.4 kg in the 3ML cows (P = 0.03). The milk and ECM yields per unit of DMI were not different between groups, but 4% FCM yield per unit of DMI was 4.1% higher in the 2ML cows than in the 3ML cows (P = 0.01). On average, the lying time of the 2ML cows was 38 min/d greater than that of the 3ML cows (P < 0.001), with no differences in rumination time.

      Blood Metabolites

      As presented in Table 5, the plasma glucose concentrations from calving until 30 DIM were 8.9% higher in the 2ML cows than in the 3ML cows (P = 0.001; Figure 3A) and tended to be higher between 31 and 50 DIM. The plasma BHB concentrations were 12.9% lower in the 2ML cows than in the 3ML cows from 0 to 30 DIM (P = 0.02; Table 5; Figure 3B), with no differences found later in lactation. The NEFA concentrations from calving until 30 DIM were 34% higher in the 2ML cows than in the 3ML cows (P = 0.002; Figure 3C) and tended to be higher between 31 and 50 DIM. The TG concentrations tended to be higher in the 2ML cows between 0 and 30 DIM, and the AST concentrations were not affected by MF.
      Figure thumbnail gr3
      Figure 3Concentrations of (A) glucose, (B) BHB, and (C) nonesterified fatty acids (NEFA) in blood from 14 d precalving until 45 DIM for high-yielding dairy cows milked 3 times daily (3ML; ▴) or twice daily (2ML; ○) until 30 DIM (vertical arrow) and then milked 3 times daily. Values are LSM, and error bars represent SE.

      Reproduction System

      The PP ovarian follicular dynamics until 34 DIM are presented in Figure 4. As shown, there was no effect of MF on follicular dynamics until 34 DIM. Two out of 8 cows (25%) in the 3ML group developed follicular cysts, whereas no cow in the 2ML group did.
      Figure thumbnail gr4
      Figure 4Ovarian follicular development from 10 until 32 DIM of high-yielding dairy cows milked 3 times daily (3ML; ▴) or twice daily (2ML; ○) until 30 DIM and then milked 3 times daily. The follicles were divided into 3 categories: (A) 3 to 5 mm (small), (B) 6 to 9 mm (medium), and (C) ≥10 mm (large). Values are LSM, and error bars represent SE.
      The time at first ovulation was determined by P4 concentrations in plasma; it was found that 10% of cows from the 3ML group ovulated ≤15 DIM compared with 40% of the 2ML cows (P < 0.03). In addition, the first ovulation occurred only at ≥30 DIM in 20% of the 3ML cows compared with 5% of the 2ML cows (P < 0.1).
      The size, hormone concentrations, and the content of preovulatory follicles are presented in Table 6. No differences were observed in the size, volume, and E2, P4, or A4 concentrations or the content in preovulatory follicles that were aspirated from 60 to 80 DIM. Although the number of cows used for analyzing the fertility data was small, the fertility performance was higher in the 2ML cows than in the 3ML cows. The average days from calving to first detected estrus was 68.0 and 69.4 d in the 2ML and 3ML cows, respectively (not significant). The conception rates from the first AI were 50 and 44.4% and from the second AI were 50 and 40% in the 2ML and 3ML cows, respectively (not significant). In addition, average days open was 9 d shorter in the 2ML cows than in the 3ML cows (100 and 109 d, respectively; not significant).

      DISCUSSION

      Milking frequency has a great effect on lactation performance in dairy cows, especially at early lactation, and improving EB during early lactation is of great interest. Therefore, in the present study we aimed to improve the metabolic status of high-yielding dairy cows by reducing the MF during the first 30 DIM and to determine the carryover effects. In general, we found a higher milk production in cows milked thrice daily but similar 4% FCM yield until 30 DIM, with no negative carryover effects on yields until 100 DIM in cows milked twice a day until 30 DIM. The metabolic status of the cows milked twice a day was improved compared with those milked 3 times a day, and the effects of reducing MF on yield were negligible when considering the yields and efficiency from 0 to 100 DIM. An improvement in EB was observed in the 2ML cows until 30 DIM, with some indications of the beneficial effects on the reproductive system.

      MF Effects on Yields

      As mentioned above,
      • Erdman R.A.
      • Varner M.
      Fixed yield responses to increased milking frequency.
      summarized data from 19 reports that examined the effect of increasing the MF from twice to thrice daily and found a fixed milk yield response of 3.5 kg/d regardless of the yield for those milked twice daily. Interestingly, the results of our study confirm this finding of a fixed milk yield, with cows producing double yields, compared with those reviewed in
      • Erdman R.A.
      • Varner M.
      Fixed yield responses to increased milking frequency.
      ; ∼46 vs. ∼22 kg/d, respectively). Furthermore, whereas the increment increase in milk yield from 2× to 3× in our study was 8.6%, in
      • Erdman R.A.
      • Varner M.
      Fixed yield responses to increased milking frequency.
      it was as high as 18.3%. This half magnitude in the response in the modern high-yielding cow compared with the low-yielding cow is very interesting with regard to the mechanism underlying the influence of MF on yield and regarding the economic impact of MF. Several studies that examined how increasing MF at early lactation affects yields found positive carryover effects in which milk production remained elevated for a period after the MF treatment was terminated (
      • Bar-Peled U.
      • Maltz E.
      • Bruckental I.
      • Folman Y.
      • Kali Y.
      • Gacitua H.
      • Lehrer A.R.
      • Knight C.H.
      • Robinzon B.
      • Voet H.
      Relationship between frequent milking or suckling in early lactation and milk production of high producing dairy cows.
      ;
      • Hale S.A.
      • Capuco A.V.
      • Erdman R.A.
      Milk yield and mammary growth effects due to increased milking frequency during early lactation.
      ). In our present study, we did not find any negative carryover effects on yields of reducing MF to twice a day during the first 30 DIM (treatment × period was significant), and the 2ML cows closed the gap with the 3ML cows within 10 d after shifting to 3 milkings a day.
      A sophisticated study by
      • Wall E.H.
      • McFadden T.B.
      The milk yield response to frequent milking in early lactation of dairy cows is locally regulated.
      in which a unilateral frequent milking (2× vs. 4×) was conducted concluded that the immediate and carryover effects of frequent milking during early lactation on milk production are regulated locally in the mammary gland. In
      • Hale S.A.
      • Capuco A.V.
      • Erdman R.A.
      Milk yield and mammary growth effects due to increased milking frequency during early lactation.
      , cows were milked 2× or 4× during the first 3 wk of lactation. A positive carryover effect was found for the remainder of lactation without a significant difference in proliferation. From the aforementioned studies, we can conclude that the effect of frequent milking during early lactation is regulated locally in the mammary gland and that the mammary cell population dynamics are not significantly altered by changes in MF. In most of the studies that examined the effects of increased MF the differences in MF between treatments was at least 2 or sometimes even 3 milkings a day, whereas in our study the difference in MF between groups was only 1 milking (2× vs. 3×). It is possible that the changes in cell turnover in our study were small so that no negative carryover effects occurred for the 2ML cows, and within 10 d they produced the same high amount as the 3ML cows. In addition, we could not find in the literature a tentative explanation for the fixed milk yield responses of 3.0 to 4.0 kg/d upon increasing MF from 2 to 3 milkings a day. Intramammary pressure was suggested as a potential regulator of mammary gland blood flow and milk secretion (
      • Wall E.H.
      • McFadden T.B.
      How the mammary gland adapts to changes in milking frequency.
      ). Taking into account the local regulation, it can be assumed that the intramammary pressure of the 2ML cows during the first 30 DIM was higher than in those milked 3 times a day. However, it is possible that no damage occurred to the secretory cells; therefore, the 2ML cows produced amounts of milk similar to the 3ML cows within 10 d, although this speculation should be elucidated in further studies.

      Effect on Milk Solids

      The fat percentage in milk was 7.8 and 5.7% higher in the 2ML cows than in the 3ML cows from 0 to 30 and 31 to 100 DIM, respectively. The decrease in milk fat content with increasing MF was found in many previous studies (
      • Bar-Peled U.
      • Maltz E.
      • Bruckental I.
      • Folman Y.
      • Kali Y.
      • Gacitua H.
      • Lehrer A.R.
      • Knight C.H.
      • Robinzon B.
      • Voet H.
      Relationship between frequent milking or suckling in early lactation and milk production of high producing dairy cows.
      ;
      • Erdman R.A.
      • Varner M.
      Fixed yield responses to increased milking frequency.
      ;
      • Hale S.A.
      • Capuco A.V.
      • Erdman R.A.
      Milk yield and mammary growth effects due to increased milking frequency during early lactation.
      ;
      • Soberon F.
      • Ryan C.M.
      • Nydam D.V.
      • Galton D.M.
      • Overton T.R.
      The effects of increased milking frequency during early lactation on milk yield and milk composition on commercial dairy farms.
      ) but not by
      • Wall E.H.
      • McFadden T.B.
      The milk yield response to frequent milking in early lactation of dairy cows is locally regulated.
      . The trend of higher milk fat content in the 2ML cows after shifting to 3× milking (31–100 DIM) was also found by
      • Hale S.A.
      • Capuco A.V.
      • Erdman R.A.
      Milk yield and mammary growth effects due to increased milking frequency during early lactation.
      . However, fat yield was higher in the 3ML cows during the first 30 DIM, a trend that was also found by
      • Alex A.P.
      • Collier J.L.
      • Hadsell D.L.
      • Collier R.J.
      Milk yield differences between 1× and 4× milking are associated with changes in mammary mitochondrial number and milk protein gene expression, but not mammary cell apoptosis or SOCS gene expression.
      but not by
      • Hale S.A.
      • Capuco A.V.
      • Erdman R.A.
      Milk yield and mammary growth effects due to increased milking frequency during early lactation.
      with increasing MF. The higher fat percentage in the 2ML cows could be explained by the 8.6% lower milk yield in these cows compared with the 3ML cows. However, the higher milk fat content in the 2ML cows from 31 to 100 DIM occurred even though the milk yields of both groups were similar, which demonstrates a positive carryover effect (treatment × period was not significant) of the lower MF until 30 DIM on fat synthesis later in lactation.
      No differences in the protein percentage in milk were observed between groups during the treatment period (until 30 DIM) or from 31 to 100 DIM. This was similar to the results of
      • Hale S.A.
      • Capuco A.V.
      • Erdman R.A.
      Milk yield and mammary growth effects due to increased milking frequency during early lactation.
      and
      • Wall E.H.
      • McFadden T.B.
      The milk yield response to frequent milking in early lactation of dairy cows is locally regulated.
      ; however, in a study with an increased MF on the mean half-udder,
      • Alex A.P.
      • Collier J.L.
      • Hadsell D.L.
      • Collier R.J.
      Milk yield differences between 1× and 4× milking are associated with changes in mammary mitochondrial number and milk protein gene expression, but not mammary cell apoptosis or SOCS gene expression.
      found a higher protein percentage in the increased (4×) MF compared with 1× milking. Because of the higher milk yields, the protein yields (kg/d) were higher in the 3ML cows during the first 30 DIM but not from 31 to 100 DIM.
      Despite the differences in milk, 4% FCM was similar between groups during the treatment period, and from 31 to 100 DIM it was even higher (5.1%) in the 2ML cows than in the 3ML cows. The energy output in milk was higher in the 3ML cows during the first 30 DIM, but from 31 to 100 DIM it was higher in the 2ML cows. Unfortunately, the FCM and ECM yields were not reported in most studies that investigated the effect of MF on yields; these yields are crucial parameters for evaluating whether the differences in milk production lie only in the milk volume or in the production potential of the cows in terms of energy output. According to our results, the differences between the 2ML and 3ML cows were mainly in the milk volume, and to a lesser extent in the total production potential, which includes the solid content in milk. This may have important economic considerations because in most countries the payment for milk includes a combination of the milk volume and the solids content. For example, according to the Israeli formula for milk payment, the average yields from 0 to 100 DIM were numerically higher in the 2ML cows than in the 3ML cows (51.9 vs. 51.3 kg/d, respectively).
      An analysis of the whole study period (0–100 DIM) revealed no difference between groups regarding milk, 4% FCM, and ECM yields and protein percentage; however, the fat percentage and yields were higher in the 2ML cows, and the protein and lactose yields were higher in the 3ML cows. In addition, the efficiency in FCM production (FCM/DMI) was higher in the 2ML cows than in the 3ML cows. These results indicated minor effects on production when reducing the MF to twice a day during the first 30 DIM. Moreover, the economic benefits of reducing MF become noteworthy when taking into consideration the additional savings regarding manpower, electricity, and the physical deterioration of facilities when milking twice daily at early lactation compared with thrice daily.

      Metabolic Status of the Cows

      The main goal of this study was to improve the metabolic status of cows during early lactation by reducing the MF. Indeed, the glucose plasma concentrations were higher and the BHB concentrations were lower in the 2ML cows during the first period (0–30 DIM), which is in agreement with the positive EB of the 2ML cows compared with the 3ML cows. Increasing the MF at early lactation in
      • Bar-Peled U.
      • Aharoni Y.
      • Robinzon B.
      • Bruckental I.
      • Lehrer R.
      • Maltz E.
      • Knight C.
      • Kali J.
      • Folman Y.
      The effect of enhanced milk yield of dairy cows by frequent milking or suckling on intake and digestibility of the diet.
      from 3 to 6 times a day increased milk yield by 22% but decreased glucose concentrations and increased NEFA plasma concentrations. In a study by
      • Rastani R.R.
      • Grummer R.R.
      • Bertics S.J.
      • Gümen A.
      • Wiltbank M.C.
      • Mashek D.G.
      • Schwab M.C.
      Reducing dry period length to simplify feeding transition cows: Milk production, energy balance, and metabolic profiles.
      , eliminating the dry period decreased milk yield by 20% but improved EB and the metabolic status of the cows compared with cows with a traditional dry period length. However, in our present study, the metabolic status of the cows was improved with a smaller reduction in milk yields, with no negative carryover effects. Additionally, the 2ML cows rested (lying time) on average 38 min/d more than the 3ML cows.
      Interestingly, although the glucose and BHB concentrations in blood indicated an improvement in the metabolic status of the 2ML cows, NEFA concentrations were higher in this group, which may indicate greater body mass mobilization in the 2ML cows. However, the BW loss from calving to wk 5 was not different between groups. The blood concentration of a metabolite is the equilibrium between the release and clearance from the blood stream. Because a great part of the released NEFA in the high-yielding dairy cow is transported to the mammary gland for milk fat synthesis, it can be assumed that the higher NEFA in the 2ML cows occurs because of the lesser clearance compared with the 3ML cows.
      The EB from 31 to 100 DIM was greater in the 3ML cows than in the 2ML cows. This was primarily because the increase in ECM yield from the first to the second period in the 2ML cows was 19% compared with 9.3% in the 3ML cows, whereas the increase in intake was almost similar (11.2% in the 2ML cows vs. 13.3% in the 3ML cows).
      • Capuco A.V.
      • Wood D.L.
      • Baldwin R.
      • Mcleod K.
      • Paape M.J.
      Mammary cell number, proliferation, and apoptosis during a bovine lactation: Relation to milk production and effect of bST.
      found that during early lactation, the milk secretion, which increased to peak levels, was due to the continued differentiation of the mammary secretory cells and the increased secretory activity per cell; however, no net mammary growth was observed. It can be assumed that the increase in milk yield from milking 2 times to milking 3 times a day after 30 DIM in the 2ML group was mainly due to an increase in the secretory capacity of the cells rather than to an increase in the mammary gland mass. The intake of the 2ML cows was similar to that of the 3ML cows, and the curve shape (Figure 2A) was identical between groups, which means that perhaps the capacity of the digestive tract was at the maximum level in both groups and therefore the EB of the 2ML cows was lower from 31 to 100 DIM. Both groups were in positive EB during the second period of the study.

      Reproductive System

      Many studies have been published investigating the association among yields, EB, and fertility (
      • Butler W.R.
      • Smith R.D.
      Interrelationships between energy balance and postpartum reproductive function in dairy cattle.
      ;
      • Heravi Moussavi A.R.
      • Gilbert R.O.
      • Overton T.R.
      • Bauman D.E.
      • Butler W.R.
      Effects of feeding fish meal and n-3 fatty acids on ovarian and uterine responses in early lactating dairy cows.
      ;
      • Cheong S.H.
      • Sá Fliho O.G.
      • Absalón-Medina V.A.
      • Pelton S.H.
      • Butler W.R.
      • Gilbert R.O.
      Metabolic and endocrine differences between dairy cows that do or do not ovulate first postpartum dominant follicles.
      ). One of the main effects of negative EB is delayed resumption of ovulation and cyclicity PP. Early resumption of ovarian cyclicity PP has been associated with improved reproductive performance (
      • Galvão K.N.
      • Frajblat M.
      • Butler W.R.
      • Brittin S.B.
      • Guard C.L.
      • Gilbert R.O.
      Effect of early postpartum ovulation on fertility in dairy cows.
      ). Indeed, cows ovulating before 21 d PP had better subsequent reproductive performance than those ovulating for the first time between 21 and 49 d PP (
      • Galvão K.N.
      • Frajblat M.
      • Butler W.R.
      • Brittin S.B.
      • Guard C.L.
      • Gilbert R.O.
      Effect of early postpartum ovulation on fertility in dairy cows.
      ). In our present study, more cows in the 2ML group ovulated before 15 DIM compared with the 3ML cows, which is in line with the improved EB in the former cows. However, although there is evidence of earlier ovulation in the 2ML cows, no differences were observed in the characteristics of preovulatory follicles aspirated from 60 to 80 DIM. This could be explained by the EB of the 3ML cows, which was higher than that of the 2ML cows from 31 DIM and later, which could have compensated for the earlier NEB in this group.

      CONCLUSIONS

      Reducing the MF from thrice to twice daily from calving until 30 DIM had beneficial effects on EB and the metabolic status of high-yielding dairy cows and only minor effects on production. Furthermore, the results of the present study confirm the fixed milk yield increase with increasing MF from twice to thrice reported by
      • Erdman R.A.
      • Varner M.
      Fixed yield responses to increased milking frequency.
      . Although the milk yield level was doubled, the 3 to 4 kg/d difference was maintained, which means that the proportion of the decrease with reducing MF is diminished in high-yielding cows. No negative carryover effects were observed, and the 2ML cows closed the gap in yields within 10 d after shifting to 3 milkings. There was also an indication of earlier PP ovulation without positive effects on preovulatory follicles aspirated from 60 to 80 DIM. It can be concluded that reducing the MF to twice a day during the first 30 DIM and then increasing it to 3 times daily may benefit the metabolic status of high-yielding dairy cows, with negligible effects on milk production.

      ACKNOWLEDGMENTS

      This research was financially supported by the Israeli Dairy Board, Yehud, Israel (grant number 362-0449-16). The authors thank the entire staff of the experimental dairy farm at the Volcani Center, Rishon LeZion, Israel, for their assistance with animal care.

      REFERENCES

        • Alex A.P.
        • Collier J.L.
        • Hadsell D.L.
        • Collier R.J.
        Milk yield differences between 1× and 4× milking are associated with changes in mammary mitochondrial number and milk protein gene expression, but not mammary cell apoptosis or SOCS gene expression.
        J. Dairy Sci. 2015; 98 (25981061): 4439-4448
        • AOAC
        Official Methods of Analysis.
        Association of Official Analytical Chemists, Arlington, VA1990
        • Bar-Peled U.
        • Aharoni Y.
        • Robinzon B.
        • Bruckental I.
        • Lehrer R.
        • Maltz E.
        • Knight C.
        • Kali J.
        • Folman Y.
        The effect of enhanced milk yield of dairy cows by frequent milking or suckling on intake and digestibility of the diet.
        J. Dairy Sci. 1998; 81 (9621246): 1420-1427
        • Bar-Peled U.
        • Maltz E.
        • Bruckental I.
        • Folman Y.
        • Kali Y.
        • Gacitua H.
        • Lehrer A.R.
        • Knight C.H.
        • Robinzon B.
        • Voet H.
        Relationship between frequent milking or suckling in early lactation and milk production of high producing dairy cows.
        J. Dairy Sci. 1995; 78 (8675755): 2726-2736
        • Beam S.W.
        • Butler W.R.
        Energy balance, metabolic hormones, and early postpartum follicular development in dairy cows fed prilled lipid.
        J. Dairy Sci. 1998; 81 (9493087): 121-131
        • Bernal-Santos G.
        • Perfield II, J.W.
        • Barbano D.M.
        • Bauman D.E.
        • Overton T.R.
        Production responses of dairy cows to dietary supplementation with conjugated linoleic acid (CLA) during the transition period and early lactation.
        J. Dairy Sci. 2003; 86 (14594242): 3218-3228
        • Bertics S.J.
        • Grummer R.R.
        Effects of fat and methionine hydroxy analog on prevention or alleviation of fatty liver induced by feed restriction.
        J. Dairy Sci. 1999; 82 (10629820): 2731-2736
        • Butler S.T.
        • Pelton S.H.
        • Butler W.R.
        Energy balance, metabolic status, and the first PP ovarian follicle wave in cows administered propylene glycol.
        J. Dairy Sci. 2006; 89 (16840609): 2938-2951
        • Butler W.R.
        Energy balance relationships with follicular development, ovulation and fertility in postpartum dairy cows.
        Livest. Prod. Sci. 2003; 83: 211-218
        • Butler W.R.
        • Smith R.D.
        Interrelationships between energy balance and postpartum reproductive function in dairy cattle.
        J. Dairy Sci. 1989; 72 (2654227): 767-783
        • Capuco A.V.
        • Wood D.L.
        • Baldwin R.
        • Mcleod K.
        • Paape M.J.
        Mammary cell number, proliferation, and apoptosis during a bovine lactation: Relation to milk production and effect of bST.
        J. Dairy Sci. 2001; 84 (11699449): 2177-2187
        • Castañeda-Gutiérrez E.
        • Overton T.R.
        • Butler W.R.
        • Bauman D.E.
        Dietary supplements of two doses of calcium salts of conjugated linoleic acid during the transition period and early lactation.
        J. Dairy Sci. 2005; 88 (15738242): 1078-1089
        • Cheong S.H.
        • Sá Fliho O.G.
        • Absalón-Medina V.A.
        • Pelton S.H.
        • Butler W.R.
        • Gilbert R.O.
        Metabolic and endocrine differences between dairy cows that do or do not ovulate first postpartum dominant follicles.
        Biol. Reprod. 2016; 94 (26632612): 18
        • Chilliard Y.
        Dietary fat and adipose tissue metabolism in ruminants, pigs, and rodents: A review.
        J. Dairy Sci. 1993; 76 (8132894): 3897-3931
        • Edmonson A. J
        • Lean I.J.
        • Weaver L.D.
        • Farver T.
        • Webster G.
        A body condition scoring chart for Holstein dairy cows.
        J. Dairy Sci. 1989; 72: 68-78
        • Erdman R.A.
        • Varner M.
        Fixed yield responses to increased milking frequency.
        J. Dairy Sci. 1995; 78 (7622729): 1199-1203
        • Esposito G.
        • Irons C.P.
        • Webb E.C.
        • Chapwanya A.
        Interactions between negative energy balance, metabolic diseases, uterine health and immune response in transition dairy cows.
        Anim. Reprod. Sci. 2014; 144 (24378117): 60-71
        • Galvão K.N.
        • Frajblat M.
        • Butler W.R.
        • Brittin S.B.
        • Guard C.L.
        • Gilbert R.O.
        Effect of early postpartum ovulation on fertility in dairy cows.
        Reprod. Domest. Anim. 2010; 45 (19845883): e207-e211
        • Grummer R.R.
        Etiology of lipid-related metabolic disorders in periparturient dairy cows.
        J. Dairy Sci. 1993; 76 (8132893): 3882-3896
        • Grummer R.R.
        • Rastani R.R.
        Why reevaluate dry period length?.
        J. Dairy Sci. 2004; 87: E77-E85
        • Gümen A.
        • Rastani R.R.
        • Grummer R.R.
        • Wiltbank M.C.
        Reduced dry periods and varying prepartum diets alter postpartum ovulation and reproductive measures.
        J. Dairy Sci. 2005; 88 (15956303): 2401-2411
        • Hale S.A.
        • Capuco A.V.
        • Erdman R.A.
        Milk yield and mammary growth effects due to increased milking frequency during early lactation.
        J. Dairy Sci. 2003; 86 (12836942): 2061-2071
        • Heravi Moussavi A.R.
        • Gilbert R.O.
        • Overton T.R.
        • Bauman D.E.
        • Butler W.R.
        Effects of feeding fish meal and n-3 fatty acids on ovarian and uterine responses in early lactating dairy cows.
        J. Dairy Sci. 2007; 90 (17183083): 145-154
      1. IDF (International Dairy Federation). 2000. IDF Standard 141c:2000 (Determination of milk fat, protein and lactose content—Guidance on the operation of mid-infrared instruments. International Dairy Federation, Brussels, Belgium.).

        • Ireland J.J.
        • Roche F.
        Development of antral follicles in cattle after prostaglandin-induced luteolysis: Changes in serum hormones, steroids in follicles fluids and gonadotropin receptors.
        Endocrinology. 1982; 111 (6291909): 2077-2086
        • Leroy J.L.
        • Sturmey R.G.
        • Van Hoeck V.
        • De. Bie J.
        • McKeegan P.J.
        • Bols P.E.
        Dietary fat supplementation and the consequences for oocyte and embryo quality: Hype or significant benefit for dairy cow reproduction?.
        Reprod. Domest. Anim. 2014; 49 (24697981): 353-361
        • Lucy M.C.
        Reproductive loss in high-producing dairy cattle: Where will it end?.
        J. Dairy Sci. 2001; 84 (11417685): 1277-1293
        • Moallem U.
        • Folman Y.
        • Bor A.
        • Arav A.
        • Sklan D.
        Effect of calcium soaps of fatty acids and administration of somatotropin on milk production, preovulatory follicular development, and plasma and follicular fluid lipid composition in high yielding dairy cows.
        J. Dairy Sci. 1999; 82 (10575602): 2358-2368
        • Moallem U.
        • Katz M.M.
        • Arieli A.
        • Lehrer H.
        Effects of peripartum propylene glycol or fats differing in fatty acid profiles on feed intake, production, and plasma metabolites in dairy cows.
        J. Dairy Sci. 2007; 90 (17638995): 3846-3856
        • Morrow D.A.
        Fat cow syndrome.
        J. Dairy Sci. 1976; 59 (965540): 1625-1629
        • National Research Council
        Nutrient Requirements of Dairy Cattle.
        6th ed. Natl. Acad. Sci., Washington, DC1989
        • National Research Council
        Nutrient Requirements of Dairy Cattle.
        7th rev. ed. National Academy Press, Washington, DC2001
        • Rastani R.R.
        • Grummer R.R.
        • Bertics S.J.
        • Gümen A.
        • Wiltbank M.C.
        • Mashek D.G.
        • Schwab M.C.
        Reducing dry period length to simplify feeding transition cows: Milk production, energy balance, and metabolic profiles.
        J. Dairy Sci. 2005; 88 (15738235): 1004-1014
        • SAS Institute
        SAS User's Guide: Statistics. Version 9.2.
        SAS Institute Inc., Cary, NC2002
        • Shingfield K.J.
        • Beever D.E.
        • Reynolds C.K.
        • Gulati S.K.
        • Humphries D.J.
        • Lupoli B.
        • Hervás G.
        • Griinari J.M.
        Effect of rumen protected conjugated linoleic acid on energy metabolism of dairy cows during early to mid-lactation.
        J. Dairy Sci. 2004; 87 (Abstr.): 307
        • Sklan D.
        • Moallem U.
        • Folman Y.
        Effects of feeding calcium soaps of fatty acids on production and reproduction responses in high production lactating cows.
        J. Dairy Sci. 1991; 74 (2045560): 510-517
        • Soberon F.
        • Ryan C.M.
        • Nydam D.V.
        • Galton D.M.
        • Overton T.R.
        The effects of increased milking frequency during early lactation on milk yield and milk composition on commercial dairy farms.
        J. Dairy Sci. 2011; 94 (21854913): 4398-4405
        • Van Soest P.J.
        • Robertson J.B.
        • Lewis B.A.
        Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition.
        J. Dairy Sci. 1991; 74: 3583-3597
        • Wall E.H.
        • McFadden T.B.
        The milk yield response to frequent milking in early lactation of dairy cows is locally regulated.
        J. Dairy Sci. 2007; 90 (17235148): 716-720
        • Wall E.H.
        • McFadden T.B.
        How the mammary gland adapts to changes in milking frequency.
        J. Anim. Sci. 2012; 90 (22205668): 1695-1707