Advertisement
Research| Volume 102, ISSUE 12, P10887-10902, December 2019

Field bean inclusion in the diet of early-lactation dairy cows: Effects on performance and nutrient utilization

Open ArchivePublished:September 20, 2019DOI:https://doi.org/10.3168/jds.2019-16513

      ABSTRACT

      The European livestock sector has a significant deficit of high-quality protein feed ingredients. Consequently there is interest in using locally grown protein grain crops to partially or completely replace imported protein feeds in dairy cow rations. Field bean (FB; Vicia faba) has been identified as a locally grown crop with significant potential. The current study was designed to examine the effects of FB on cow performance and nutrient utilization in the diet of early-lactation dairy cows, including high levels of FB (up to 8.4 kg/cow per day). The experiment used 72 dairy cows in a 3-treatment continuous design (from calving until wk 20 of lactation). All cows were given ad libitum access to a mixed ration comprising grass silage and concentrates [45:55 on a dry matter (DM) basis]. Concentrates offered contained either 0, 349, or 698 g of FB/kg of concentrate (treatments FB0, FB-Low, and FB-High, respectively), with FB completely replacing soybean meal, rapeseed meal, maize gluten, and wheat in the concentrate for the FB-High treatment. Following completion of the 20-wk experiment, ration digestibility, nutrient utilization, and methane (CH4) production were measured using 4 cows from each treatment. Neither silage DM intake, total DM intake, nor milk yield were affected by treatment. Cows on FB0 had a higher milk fat content than those on FB-High, and cows on FB0 and FB-Low had higher milk protein contents than did those on FB-High. Field bean inclusion increased the degree of saturation of milk fat produced. Milk fat yield, milk protein yield, and milk fat plus protein yield were higher with FB0 than with either FB-Low or FB-High. Treatment had no effect on the digestibility of DM, organic matter, nitrogen (N), gross energy, or neutral detergent fiber, whereas digestibility of acid detergent fiber was higher with FB0 than with FB-High. Neither the efficiency of gross energy or N utilization, nor any of the CH4 production parameters examined, were affected by treatment. Similarly, none of the fertility or health parameters examined were affected by treatment. The reduction in milk fat observed may have been due to the higher starch content of the FB-High diet, and the reduction in milk protein may have been due to a deficit of methionine in the diet. It is likely that these issues could be overcome by changes in ration formulation, thus allowing FB to be included at the higher range without loss in performance.

      Key words

      INTRODUCTION

      The increasing milk yield potential of dairy cows in many European countries has led to a requirement for more nutrient-dense diets. Increased nutrient density has often been achieved through the adoption of higher-concentrate feed levels, and this in turn has increased the demand for high-quality protein ingredients. However, European agriculture has a significant deficit of high-quality protein feedstuffs (
      • Watson C.
      • Reckling M.
      • Preissel S.
      • Bachinger J.
      • Bergkvist G.
      • Kuhlman T.
      • Lindstrom K.
      • Nemecek T.
      • Topp C.F.E.
      • Vanhatalo A.
      • Zander P.
      • Murphy-Bokern D.
      • Stoddard F.L.
      Grain legume production and use in European agricultural systems.
      ). This is a particular problem in the United Kingdom, driving considerable interest in increased use of locally grown protein ingredients (
      • Wilkins R.J.
      • Jones R.
      Alternative home-grown protein sources for ruminants in the United Kingdom.
      ). One protein crop of particular interest, especially in the cooler, wetter regions of the United Kingdom, is the field bean (FB; Vicia faba).
      Although FB have a moderate CP content (280 g/kg of DM), they have a relatively high starch content (400 g/kg of DM;
      • Ewing W.N.
      The Feeds Directory, Volume 1: Commodity Products Guide.
      ) and, as such, might appear to be ideal in ruminant diets. However, the inclusion of FB in ruminant diets is normally limited due to the perceived risk associated with antinutritional substances, including trypsin inhibitors, tannins, lectins, and protease inhibitors, which can reduce DM intake and animal performance (
      • Dvořák R.
      • Pechová A.
      • Pavlata L.
      • Klejdus B.
      • Kovařčík K.
      • Dostálová J.
      • Culková J.
      • Filípek J.
      • Svajdlenka E.
      • Capkova V.
      Reduction in the content of antinutritional substances in Fava beans (Vicia faba) by different treatments.
      ). In addition, FB contain several phytoestrogens, some of which are known to have negative effects on reproductive performance (
      • Zdunczyk S.
      • Piskula M.
      • Janowski T.
      • Baranski W.
      • Ras M.
      Concentrations of isoflavones in blood plasma of dairy cows with different incidence of silent heat.
      ), which has also contributed to the limited inclusion of FB in dairy cow diets.
      The use of FB in dairy cow diets has been examined in several studies, with the findings sometimes conflicting. For example, although neither DM intake nor milk yield were affected by FB inclusion in studies by
      • Ingalls J.
      • McKirdy J.
      Faba bean as a substitute for soybean meal or rapeseed meal in rations for lactating cows.
      and
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ,
      • Ramin M.
      • Höjer A.
      • Fogelberg F.
      • Hetta M.
      The effects of legume seeds on the lactation performance of dairy cow fed grass silage-based diets.
      observed a reduction in milk yield with FB inclusion. In contrast,
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      found both DM intake and milk yield to be reduced when rapeseed meal was replaced by FB.
      These studies have demonstrated that in some circumstances FB can be included in dairy cow diets at levels between 4.7 and 5.0 kg/d (
      • Ingalls J.
      • McKirdy J.
      Faba bean as a substitute for soybean meal or rapeseed meal in rations for lactating cows.
      ;
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ) with no adverse effects on DM intake and milk production, but higher inclusion levels do not appear to have been examined previously. Given the high starch content of FB, which suggests that it could replace both conventional protein and energy components of the diet, the effects of higher FB inclusion levels need to be examined. In addition, as highlighted by
      • Tufarelli V.
      • Khan R.U.
      • Laudadio V.
      Evaluating the suitability of field beans as a substitute for soybean meal in early lactating dairy cow: Production and metabolic responses.
      , information on the effects of FB inclusion in the diet of dairy cows in early lactation, the period that encompasses the breeding season, is limited at present. This is important, given concerns about the presence of estrogenic compounds in FB, although there is no evidence that dairy cow fertility is reduced when FB are offered. Consequently, the current study was designed to examine the effects on cow performance, nutrient utilization efficiency, and fertility when early-lactation dairy cows were offered diets containing higher levels of FB than examined previously. We hypothesized that these parameters would be unaffected by these higher levels of FB inclusion.

      MATERIALS AND METHODS

      This experiment was conducted at the Agri-Food and Biosciences Institute, Hillsborough, Northern Ireland (NI). All experimental procedures were conducted under an experimental license granted by the Department of Health, Social Services and Public Safety for NI, in accordance with the Animals (Scientific Procedures) Act 1986.

      Animals and Housing

      This continuous-design 20-wk experiment involved 72 Holstein-Friesian dairy cows, 33 primiparous and 39 multiparous; mean parity 2.07 (SD 1.19). During the 3-wk period before calving, cows had ad libitum access to grass silage, while a dry-cow mineral and vitamin mix, along with calcined magnesite (Trouw Nutrition, Cheshire, UK), were mixed with the silage to achieve a target intake of 100 g/cow per day and 50 g/cow per day, respectively. Cows were moved to a straw-bedded maternity pen approximately 24 to 48 h before calving, based on behavioral observations. Following calving (normally within 24 h), cows were transferred to an experimental free-stall barn with solid concrete floors, where they had access to cubicles fitted with rubber mats, which were bedded with sawdust twice daily. The floor was scraped every 3 h, using an automated slurry scraper system.

      Treatments

      The variety of FB (Vicia faba var. Fuego) used in the study was described by
      • Flores F.
      • Hybl M.
      • Knudsen J.C.
      • Marget P.
      • Muel F.
      • Nadal S.
      • Narits L.
      • Raffiot B.
      • Sass O.
      • Solis I.
      • Winkler J.
      Adaptation of spring faba bean types across European climates.
      as having “colored flower, high standing ability, high yield potential, big seeds.” Colored-flower varieties of FB have long been known to have higher tannin contents than white-flowered varieties (
      • Newton S.D.
      • Hill G.D.
      The composition and nutritive value of field beans.
      ), with
      • Abdulla J.M.
      The use of field beans as a feed for broiler chickens. PhD Thesis.
      recently recording the tannin content of the Fuego variety as 6.1 mg/g of DM. This variety was chosen for use in this study because it has been successfully grown by local arable and livestock farmers in NI over several years. The study by
      • Flores F.
      • Hybl M.
      • Knudsen J.C.
      • Marget P.
      • Muel F.
      • Nadal S.
      • Narits L.
      • Raffiot B.
      • Sass O.
      • Solis I.
      • Winkler J.
      Adaptation of spring faba bean types across European climates.
      also demonstrated that this variety could be successfully grown over wide areas of Europe, with this variety having been evaluated previously in a dairy cow study in Lithuania (
      • Kudlinskienė I.
      • Lionikienė J.
      • Klementavičiūtė J.
      • Stanytė G.
      • Gružauskas R.
      • Stankevičius R.
      • Urbšienė D.
      • Bliznikas S.
      • Racevičiūtė-Stupelienė A.
      • Ots M.
      • Kass M.
      Effect of extruded field beans “Fuego” (Vicia Faba) on dairy cow's performance and milk sensory properties.
      ). The beans offered within this study were grown in NI. The crop was sown during spring 2016, harvested on September 28, 2016, and then dried at 80°C for approximately 2 h to achieve a moisture content of approximately 14%.
      Cows were randomly allocated to 1 of 3 treatments at calving (FB0, FB-Low, FB-High), with primiparous and multiparous cows allocated separately. A check was made to ensure that each treatment group remained balanced for parity, BW, and BCS at drying-off, and, in the case of multiparous cows, for previous lactation 305-d milk yield and milk fat and protein contents.
      Treatments differed in concentrate type offered (ingredient composition in Table 1). The concentrate offered with FB0 contained no FB, whereas concentrates offered with FB-Low and FB-High contained 349 and 698 g of FB per kilogram (fresh basis), respectively. The experimental concentrates (in the form of a meal) were mixed with grass silage (forage:concentrate ratio of 45:55 on a DM basis), together with 320 g/cow per day of straw and 50 g/cow per day of a rumen buffer (Acid Buf, Celtic Sea Minerals, Cork, Ireland; buffering capacity approximately 2.5× that of sodium bicarbonate), and offered in the form of a mixed ration. The grass silage component of the diet was produced from a primary growth (approximately 75% of cow feeding days) and a primary regrowth (approximately 25% of cow feeding days) of perennial ryegrass (Lolium perenne)–based swards. The rations were prepared and offered daily at approximately 0900 h (at proportionally 1.07 of the previous day's intake), with uneaten ration removed the following day at approximately 0800 h. Rations were prepared using a feeder wagon (Vari-Cut 12, Redrock, Armagh, NI). The total silage requirement for all 3 treatments was initially mixed for approximately 4 min and then deposited on a clean covered concrete yard. Silage for each individual treatment was then removed from this “pile” and returned to the feeder wagon. The appropriate quantities of concentrate, straw, and buffer were added to the mix, and mixing continued for another 6 min. The rations were then transferred from the wagon to a series of feedboxes mounted on weigh scales, with cows accessing feed in these boxes via an electronic identification system, enabling individual cow intakes to be recorded daily (Controlling and Recording Feed Intake, BioControl, Rakkestad, Norway). Cows had free access to fresh water at all times. In addition, all cows were offered 1.0 kg per day of soy hulls via an in-parlor feeder during milking (0.5 kg at each milking). Cows remained on the experiment until d 140 of lactation.
      Table 1Ingredient composition (g/kg) of concentrates offered during the experiment
      IngredientTreatment
      FB0FB-LowFB-High
      Field beans (FB)0349698
      Soybean meal170850
      Rapeseed meal150750
      Maize gluten feed74370
      Maize meal245178110
      Wheat174870
      Soy hulls (toasted)140139140
      Molaferm
      United Molasses, Belfast, UK.
      252525
      Calcined magnesite666
      Limestone (CaCO3)765
      Dicalcium phosphate0510
      Salt642
      Mineral and vitamin premix
      Devenish Nutrition, Belfast, UK. Contains a minimum of (per kg) vitamin A retinyl acetate (2,500,000 IU), vitamin D3 cholecalciferol (500,000 IU), vitamin E alpha tocopherol (18,750 mg), vitamin B12 (7,500 μg/μg), biotin (500 mg), iodine from calcium iodate anhydrous (1,875 mg), selenium from sodium selenite (150 mg), selenium from organic form of selenized yeast inactivated (25 mg), cobalt from coated granulated cobalt carbonate (40 mg), manganese from manganese(II) oxide (12,500 mg), copper from copper(II) sulfate pentahydrate (5,625 mg), copper from copper(II) chelate of amino acid hydrate (1,875 mg), zinc from zinc oxide (18,750 mg), and zinc from zinc chelate of amino acid hydrate (6,250 mg).
      444
      1 United Molasses, Belfast, UK.
      2 Devenish Nutrition, Belfast, UK. Contains a minimum of (per kg) vitamin A retinyl acetate (2,500,000 IU), vitamin D3 cholecalciferol (500,000 IU), vitamin E alpha tocopherol (18,750 mg), vitamin B12 (7,500 μg/μg), biotin (500 mg), iodine from calcium iodate anhydrous (1,875 mg), selenium from sodium selenite (150 mg), selenium from organic form of selenized yeast inactivated (25 mg), cobalt from coated granulated cobalt carbonate (40 mg), manganese from manganese(II) oxide (12,500 mg), copper from copper(II) sulfate pentahydrate (5,625 mg), copper from copper(II) chelate of amino acid hydrate (1,875 mg), zinc from zinc oxide (18,750 mg), and zinc from zinc chelate of amino acid hydrate (6,250 mg).

      Cow Measurements

      During the experiment cows were milked twice daily (between 0500 and 0700 h, and between 1500 and 1700 h) using a 50-point rotary milking parlor (BouMatic, Madison, WI), with milk yields recorded automatically at each milking and a total daily milk yield for each cow for each 24-h period calculated. Milk samples were taken during 2 consecutive milkings each week, a preservative tablet was added (Broad Spectrum Microtabs II, D and F Control Systems, Advanced Instruments Inc., Norwood, MA), and samples were stored at 4°C until analyzed. Samples were analyzed for fat, protein, and lactose using an infrared milk analyzer (Milkoscan model FT+; Foss Electric, Hillerod, Denmark), and a weighted concentration of each constituent was determined for the 24-h sampling period. A single milk sample was collected for progesterone analysis from each cow on 2 occasions each week (Monday and Thursday p.m.) until 60 d postpartum. Samples were preserved and stored as above, for a maximum of 3 wk, with progesterone concentrations subsequently determined using a competitive ELISA kit (Ridgeway Science Ltd., Gloucestershire, UK). In addition, at 10 and 20 wk post-calving (±1 wk) a milk sample was taken in proportion to yield during 2 successive milkings (a.m. and p.m.), the 2 samples bulked, and the bulked sample frozen at −20°C. These samples were subsequently analyzed for milk FA, as follows: milk fat was extracted from 1.0 mL of homogenized milk using a chloroform methanol extraction method (
      • Bligh E.G.
      • Dyer W.J.
      A rapid method of total lipid extraction and purification.
      ), and FAME prepared. The FA composition was determined using GLC, with an aliquot (1.0 µL) of the FAME extract injected onto a CP Sil88 capillary column (100 m × 0.25-mm internal diameter × 0.2-µm film thickness) in a, Agilent 7800 gas chromatograph (both from Agilent Technologies, Santa Clara, CA), equipped with a temperature programmable injector operated in the split mode and a flame ionization detector. The oven was initially held at 50°C for 4 min, then increased at 8°C/min to 110°C, then at 5°C/min to 170°C (hold time 10 min), and finally at 2°C/min to 225°C (hold time 30 min). Fatty acids were identified by their retention time with reference to commercially available FA standards (37 Supelco FAME mix) and individual standards for those not in the mix (Sigma-Aldrich Co. Ltd., Gillingham, UK), and were quantified using C13 FAME as an internal standard.
      Body weight was recorded twice daily (immediately after each milking) using an automated weighbridge, and a mean weekly BW for each cow was determined. The BCS of each cow was estimated fortnightly according to
      • Edmonson A.J.
      • Lean I.J.
      • Weaver L.D.
      • Farver T.
      • Webster G.
      A body condition scoring chart for Holstein dairy cows.
      by a trained technician. Fecal scores were assessed at 6 and 12 wk post-calving (±1 wk) using a scale of 1 to 5, as described by
      • Ireland-Perry R.L.
      • Stallings C.C.
      Fecal consistency as related to dietary composition in lactating Holstein cows.
      . Blood samples were collected from the tail of each cow before feeding, at 2, 4, 6, 8, 10, and 14 wk (±3 d) post-calving. Samples were collected into evacuated tubes (Becton Dickinson, Oxford, UK), which were coated with either a clot activator or fluoride oxalate. The blood samples were centrifuged (1,800 × g at 17°C for 30 min) to obtain serum (tubes with clot activator) or plasma (fluoride oxalate tubes), which were separated and stored at −20°C until analyzed. Plasma was analyzed for glucose concentrations, whereas serum was analyzed for BHB, nonesterified fatty acids (NEFA), and urea concentrations. Plasma glucose concentrations were determined using the hexokinase method (Roche Diagnostics Ltd., Burgess Hill, UK), and serum biochemistry analysis was carried out on a Randox Imola chemistry analyzer system, using Randox reagent kits (Randox Laboratories, Crumlin, UK).
      A clinical examination of vaginal mucus was conducted at 2, 4, and 6 wk postpartum, as an indicator of immune status. The methodology and scoring system used was as described by
      • Little M.W.
      • O'Connell N.E.
      • Welsh M.D.
      • Mulligan F.J.
      • Ferris C.P.
      Concentrate supplementation of a diet based on medium-quality grass silage for 4 weeks prepartum: Effects on cow performance, health, metabolic status, and immune function.
      , with the combined results expressed on a scale of 0 to 5. In the current study, scores were subsequently grouped into 3 categories for analysis (0, 1, ≥2).

      Nutrient Utilization

      On completion of the 20-wk feeding study, 4 cows from each treatment were selected for use in a nutrient-utilization study, with the cows selected from each treatment group balanced for daily milk yield and BW. Cows were transferred in pairs from the main experimental group into a digestibility unit, with pairs transferred at approximately 0900 h on Monday or Thursday of each week. The first pair transferred comprised cows from treatments FB0 and FB-Low; the second, cows from treatments FB-High and FB0; and the third, cows from FB-Low and FB-High. This pattern was repeated for the remaining 3 pairs. Cows were tied by the neck in individual stalls fitted with a rubber mat, and continued to be offered their experimental rations via a feedbox located at the front of each stall. Rations were offered ad libitum daily at 0900 h, at proportionally 1.07 of the previous d intake, and uneaten ration was removed the following day at 0800 h. Soy hulls were offered (0.5 kg at each milking) via a plastic feed bucket, which was placed within the feedboxes during milking (at 0630 and 1630 h). Cows had access to fresh water at all times via a drinker located within each stall.
      Measurements of nutrient utilization commenced 24 h after cows were placed in the digestibility unit and comprised a 6-d feeding period, followed 48 h later by a 6-d total feces and urine collection period. Following collection of feces and urine on d 3 of the 6-d collection period, each pair of cows was transferred into 2 indirect open-circuit respiration chambers for measurement of CH4 production, with each chamber used twice with each treatment, to remove any possible chamber effects. Cows remained in the chambers for 72 h, with measurements of CH4 and CO2 production recorded during the final 48-h period used in the subsequent analysis. The milking routine, milk sampling, milk analysis, feces and urine collection methodologies, sample management throughout the 6-d collection period, specifications of the respiration chambers, details of measurements of CH4 production, laboratory analysis processes of all feed, milk, urine, and feces, have been described in our previous work (
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ). In the present study, a mean CH4 recovery rate of 98.5% was recorded.

      Feed Analysis

      Samples of grass silage were taken daily throughout the experiment and dried at 85°C for 24 h to determine oven DM (ODM) content. Subsamples of the dried milled silages were taken 3 times weekly and bulked for every 14 d, with the bulked sample analyzed for NDF, ADF, and ash concentrations. A fresh sample of grass silage was taken weekly and analyzed for concentrations of N, ammonia-N, lactic acid, acetic acid, ethanol, propanol, and gross energy (GE), and for pH, whereas ME content was predicted using near-infrared reflectance spectroscopy, as described by
      • Park R.S.
      • Agnew R.E.
      • Gordon F.J.
      • Steen R.W.J.
      The use of near infrared reflectance spectroscopy (NIRS) on undried samples of grass silage to predict chemical composition and digestibility parameters.
      . A sample of each experimental concentrate offered (FB0, FB-Low, and FB-High) was taken weekly, bulked for each 2-wk period, and subsequently dried at 85°C for 24 h to determine ODM content. Additional samples of each experimental concentrate were taken weekly, bulked for each 2-wk period, dried at 60°C for 48 h, and analyzed for N, NDF, ADF, ash, GE, and starch concentrations. These latter samples were further bulked for each 4-wk period and analyzed using the procedure developed by
      • Licitra G.
      • Hernandez T.M.
      • Van Soest P.J.
      Standardization of procedures for nitrogen fractionation of ruminant feeds.
      for ADIN and NDIN. Acid detergent insoluble protein (ADIP) and neutral detergent insoluble protein (NDIP) were calculated as ADIP = 6.25 × ADIN and NIDP = 6.25 × NDIN, respectively.
      Straw and soy hulls offered were sampled weekly, the samples bulked every 4 wk, and the bulked samples dried at 85°C for 24 h to determine ODM content. Dried samples were subsequently analyzed for N, GE, ADF, NDF, and ash content. Laboratory analyses of all feed samples were undertaken as described previously by
      • Purcell P.J.
      • Law R.A.
      • Gordon A.W.
      • McGettrick S.A.
      • Ferris C.P.
      Effect of concentrate feeding method on the performance of dairy cows in early to mid lactation.
      .

      Statistical Analysis

      Two primiparous cows from treatment FB-Low did not complete the study due to lameness and very low milk yield, respectively. Mean weekly data for DMI, milk production, milk composition, milk constituent yields, and BW over the 20-wk experimental period were analyzed using REML repeated measures analysis (autoregressive order 1). The mixed model used included the following terms as fixed effects: Parity (1, 2, 3, >3) + Wk (1–20, as the repeated-measures time factor) + FB inclusion level (FB0, FB-Low, FB-High) + Wk × FB inclusion level, and cow was fitted as a random effect. Locomotion score and BCS data were analyzed using the same model, except that fortnightly data were used as the repeated-measures time factor. Milk FA and blood metabolites were analyzed using a similar model, except that for the former, sampling week (10 and 20 post-calving) was used as the repeated-measure time factor, whereas with the latter, sampling week (2, 4, 6, 8, 10, and 14-post calving) was used as the repeated-measures time factor. Continuous BW, BCS, and fertility data were analyzed by ANOVA, with parity included as a covariate. For each continuous variable, polynomial contrasts of order 1 were fitted to the treatment effects to test the data for linear effects associated with increasing FB level. Binomial fertility and health data were analyzed via generalized linear model regression analysis using the binomial distribution with a logit link function. The model included treatment as a term, with parity included as a covariate in the case of mastitis incidence. Mucus scores of 2 and 3 were combined into a single category (2–3), yielding categories of 0, 1, and 2–3 for analysis. These categorical data were analyzed using a generalized linear mixed model (ordinal logistic regression) with random effects (proportional odds model). The factorial arrangement of treatment (FB0, FB-Low, or FB-High) and sampling week (wk 6 and 12 for feces scores; wk 2, 4, and 6 for mucus scores) were fitted as fixed effects, and cow was fitted as a random effect. Significance was identified using chi-squared. Mean nutrient utilization data over the 6-d measurement period, and mean CH4 production data over the 2-d measurement period, were analyzed using ANOVA. All data were analyzed using GenStat (Release 18.1; VSN International Ltd., Oxford, UK).

      RESULTS

      The silage offered had ODM, CP, and ME contents of 293 g/kg, 145 g/kg of DM, and 11.0 MJ/kg of DM, respectively (Table 2). The 3 experimental concentrates had similar CP contents (mean, 223 g/kg of DM), and the starch contents of FB0, FB-Low, and FB-High concentrates were 291, 313, and 338 g/kg of DM, respectively. However, the NDIP of the concentrates decreased, being 112, 98, and 90 g/kg of DM for FB0, FB-Low, and FB-High, respectively, and the respective values for ADIP were 28.5, 21.6, and 16.6. The straw offered had CP, ADF, NDF, and ash contents of 46 (SD ±11.2), 501 (±9.6), 866 (±20.7), and 57 (±13.2) g/kg of DM, and a GE content of 18.5 (±0.23) MJ/kg of DM. The respective values for the soy hulls offered were 111 (±4.3), 492 (±11.3), 685 (±29.4), 53 (±1.1) g/kg of DM, and 17.5 (±0.06) MJ/kg of DM, with a starch content of 13.9 (±7.84) g/kg of DM.
      Table 2Chemical composition of the grass silage and experimental concentrates (FB0, FB-Low, and FB-High)
      Concentrates FB0, FB-Low, and FB-High contained 0, 349, and 698 g of field beans (FB) per kilogram, respectively.
      offered during the experiment
      ItemGrass silageSDConcentrate
      FB0SDFB-LowSDFB-HighSD
      Oven DM (g/kg)29361.58925.28895.588210.4
      Volatile-corrected oven DM (g/kg)31161.8
      CP (g/kg of DM)14522.12245.822210.12239.2
      Neutral detergent insoluble protein (g/kg of DM)11215.79812.49018.8
      Acid detergent insoluble protein (g/kg of DM)28.52.8321.62.7617.61.06
      Ash (g/kg of DM)876.1685.7644.6594.0
      ADF (g/kg of DM)29423.615921.916129.016424.0
      NDF (g/kg of DM)48042.027942.526741.526631.9
      Gross energy (MJ/kg of DM)18.91.4718.30.1218.10.1318.00.10
      Starch (g/kg of DM)29125.331338.833827.2
      Ammonia N (g/kg of total N)768.0
      pH3.870.216
      Lactic acid (g/kg of DM)11456.8
      Acetic acid (g/kg of DM)2310.6
      Ethanol (g/kg of DM)1410.2
      Propanol (g/kg of DM)0.90.94
      Metabolizable energy (MJ/kg of DM)11.00.5512.7
      Estimated from composition of individual concentrate ingredients.
      12.5
      Estimated from composition of individual concentrate ingredients.
      12.3
      Estimated from composition of individual concentrate ingredients.
      Calcium
      Estimated from composition of individual concentrate ingredients.
      (g/kg of DM)
      6.912.918.5
      Phosphorus
      Estimated from composition of individual concentrate ingredients.
      (g/kg of DM)
      5.55.55.5
      1 Concentrates FB0, FB-Low, and FB-High contained 0, 349, and 698 g of field beans (FB) per kilogram, respectively.
      2 Estimated from composition of individual concentrate ingredients.
      Neither total DMI nor milk yield were affected by treatment (P > 0.05; Table 3). Silage DMI was estimated as 9.7, 9.6, and 9.8 kg/d for FB0, FB-Low, and FB-High, respectively (based on ration formulation), and the respective values for concentrates were 11.8, 11.7, and 12.2 kg/d. Cows on FB0 had higher milk fat content than those on FB-High (P = 0.031), whereas cows on FB0 and FB-Low had higher milk protein content than those on FB-High (P < 0.001). Milk lactose content was not affected by treatment (P > 0.05). Milk fat yield (P < 0.001), milk protein yield (P = 0.035), and milk fat plus protein yield (P = 0.007) were higher with FB0 than with either FB-Low or FB-High. All of these DMI and milk production parameters varied with time post-calving (P < 0.001; see Figure 1 for DMI and Figure 2 for milk yield), but we found no treatment × time interaction for any of these parameters. Milk fat, protein, and lactose concentrations showed linear decrease with FB inclusion (P = 0.012, P < 0.001, and P = 0.046, respectively), as did milk fat yield (P = 0.002), milk protein yield (P = 0.020), and milk fat plus protein yield (P = 0.006).
      Table 3Effects of field bean (FB) inclusion level in dairy cow concentrates on mean DMI and milk production parameters over the 20-wk experimental period, and on mean milk fatty acid (FA) concentrations (mean of 2 sampling occasions)
      ItemTreatmentSEMP-value
      FB0FB-LowFB-HighTreatmentTimeTreatment × timeLinear
      Total DMI (kg/d)21.921.622.10.500.497<0.0010.7220.611
      Milk yield (kg/d)35.733.233.90.890.272<0.0010.7230.539
      Milk fat (g/kg)42.8
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      42.5
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      41.3
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.420.031<0.0010.2970.012
      Milk protein (g/kg)33.8
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      33.6
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      32.2
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.33<0.001<0.0010.460<0.001
      Milk lactose (g/kg)48.348.247.90.1120.113<0.0010.8830.046
      Milk fat yield (kg/d)1.52
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      1.39
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      1.39
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.025<0.001<0.0010.2480.002
      Milk protein yield (kg/d)1.20
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      1.11
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      1.09
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.0260.035<0.0010.7120.020
      Milk fat + protein yield (kg/d)2.71
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      2.49
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      2.47
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.0500.007<0.0010.1960.006
      Milk FA (g/100 g of total FAME identified)
       Total SFA
      Sum of C4:0, C6:0, C8:0, C10:0, C11:0, C12:0, C14:0, C15:0, C16:0, C17:0, C18:0, C20:0, C21:0, C22:0, C23:0, and C24:0.
      73.9
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      75.4
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      75.2
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.420.0280.9080.4480.138
       Total MUFA
      Sum of C14:1, C15:1, C16:1, C17:1, C18:1 cis-9, C18:1 cis-11, C20:1 cis-11, C22:1 cis-13, and C24:1 cis-15.
      22.8
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      21.4
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      21.5
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.370.0160.1970.5510.016
       Total PUFA
      Sum of C18:2, C18:2 cis-9,trans-11, C18:3n-3, C18:3n-6, C20:2 cis-11,14, C20:3n-3, C20:3n-6, C20:4n-6, C20:5n-3, C22:2 cis-13,16, C22:5n-3, and C22:6n-3.
      3.33.23.30.070.495<0.0010.2400.948
       C14:012.912.812.60.180.3280.0600.1440.138
       C16:035.8
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      38.4
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      39.4
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.57<0.0010.0020.571<0.001
       C18:09.8
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      9.2
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      8.4
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.270.002<0.0010.835<0.001
       C18:1 cis-918.2
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      17.0
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      17.1
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.350.0220.9120.5120.020
       C18:21.771.721.810.040.292<0.0010.3870.565
       C18:2 cis-9,trans-11 (CLA)0.54
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.45
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.40
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.019<0.0010.1330.751<0.001
       C18:3n-30.49
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.52
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.56
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.0150.006<0.001<0.0010.002
      a–c Means with the same superscript within a row do not differ significantly (P > 0.05).
      1 Sum of C4:0, C6:0, C8:0, C10:0, C11:0, C12:0, C14:0, C15:0, C16:0, C17:0, C18:0, C20:0, C21:0, C22:0, C23:0, and C24:0.
      2 Sum of C14:1, C15:1, C16:1, C17:1, C18:1 cis-9, C18:1 cis-11, C20:1 cis-11, C22:1 cis-13, and C24:1 cis-15.
      3 Sum of C18:2, C18:2 cis-9,trans-11, C18:3n-3, C18:3n-6, C20:2 cis-11,14, C20:3n-3, C20:3n-6, C20:4n-6, C20:5n-3, C22:2 cis-13,16, C22:5n-3, and C22:6n-3.
      Figure thumbnail gr1
      Figure 1Silage DMI (dashed lines) and total DMI (solid lines) responses of dairy cows over the first 20 wk of lactation, to concentrates containing a range of field bean (FB) inclusion levels (FB0, FB-Low, or FB-High). Silage DMI: SEM = 0.22; treatment P = 0.497; time P < 0.001; treatment × time P = 0.722. Total DMI: SEM = 0.50; treatment P = 0.497; time P < 0.001; treatment × time P = 0.722.
      Figure thumbnail gr2
      Figure 2Milk yield response of dairy cows to concentrates containing a range of field bean (FB) inclusion levels (FB0, FB-Low, or FB-High) over the first 20 wk of lactation. SEM = 0.89; treatment P = 0.272; time P < 0.001; treatment × time P = 0.723.
      The total concentration of SFA in milk fat was lower with FB0 than with either FB-Low or FB-High (P = 0.028), but the reverse was true for the total concentration of MUFA (P = 0.016; Table 3). Treatment had no effect on the total PUFA content of milk fat (P > 0.05). Concentrations of C18:0 (P = 0.002), C18:1 cis-9 (P = 0.022), and C18:2 cis-9,trans-11 (P < 0.001) decreased with increasing FB inclusion; however, concentrations of C16:0 (P < 0.001) and C18:3n-3 (P < 0.006) increased. Concentrations of C14:0 and C18:2 were not affected by treatment (P > 0.05). Although total PUFA concentrations (P < 0.001) and concentrations of C16:0 (P = 0.002), C18:0 (P < 0.001), C18:2 (P < 0.001), and C18:3n-3 (P < 0.001) differed between the 2 sampling periods (time) with the exception of the latter (P < 0.001), we found no treatment × time interactions for any of the FA recorded (P > 0.05). Total MUFA concentrations (P = 0.016) and concentrations of C18:0 (P < 0.001), C18:1 cis-9 (P = 0.020), and CLA (P < 0.001) showed a linear decrease with FB inclusion, but concentrations of C16:0 (P < 0.001) and C18:3n-3 (P = 0.002) showed the reverse trend.
      Although average BCS was unaffected by treatment (P < 0.05), end-of-study BCS tended to increase with increasing levels of FB inclusion (P = 0.054; Table 4), with this overall linear effect significant (P = 0.020). None of the BW parameters examined (average BW, end-of-study BW, nadir BW, BW loss to nadir, and days to nadir BW) were affected by treatment (P > 0.05). Both BCS (P = 0.046) and BW (P < 0.001) varied with time (Figure 3, Figure 4, respectively), but we found no treatment × time interactions for either (P < 0.05). Plasma glucose and BHB concentrations were unaffected by treatment (P > 0.05), although both increased with time post-calving (Figure 5a, P = 0.010; Figure 5b, P = 0.002, respectively). Plasma NEFA concentrations were significantly higher with FB0 than with FB-High (P = 0.022); however, plasma urea concentrations followed the reverse trend (P < 0.001). Plasma NEFA concentrations decreased with time post-calving (Figure 5c, P < 0.001), but plasma urea concentrations increased (Figure 5d, P < 0.001). We detected no treatment × time interactions for any of the blood metabolites examined (P > 0.05). With increasing FB inclusion levels, analysis revealed a tendency for linear decrease in plasma glucose concentrations (P = 0.067), whereas plasma NEFA concentrations showed a linear decrease (P = 0.005). Plasma urea levels increased linearly with FB inclusion (P < 0.001).
      Table 4Effects of field bean (FB) inclusion level in dairy cow concentrates on BCS, BW, and plasma metabolite concentrations
      ItemTreatmentSEMP-value
      FB0FB-LowFB-HighTreatmentTimeTreatment × timeLinear
      Average BCS2.492.532.510.0240.2590.0460.2480.371
      End-of-study BCS2.422.532.570.0470.0540.020
      Average BW (kg)5986005998.00.899<0.0010.7600.646
      End-of-study BW (kg)6146216179.60.8700.779
      Nadir BW (kg)5705655698.10.9020.963
      Loss to nadir BW (kg)2320304.80.2760.270
      Days to nadir BW50354611.20.6080.775
      Plasma metabolites
       Glucose (mmol/L)3.813.733.710.0360.1530.0100.9540.067
       BHB (mmol/L)0.510.510.480.0180.4190.0020.1230.206
       NEFA
      NEFA = nonesterified fatty acids.
      (mEq/L)
      0.43
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.39
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.36
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.0200.022<0.0010.6330.005
       Urea (mmol/L)4.63
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      5.28
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      5.79
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.103<0.001<0.0010.567<0.001
      a–c Means with the same superscript within a row do not differ significantly (P > 0.05).
      1 NEFA = nonesterified fatty acids.
      Figure thumbnail gr3
      Figure 3Body condition score (BCS) response of dairy cows to concentrates containing a range of field bean (FB) inclusion levels (FB0, FB-Low, or FB-High) over the first 20 wk of lactation. SEM = 0.024; treatment P = 0.259; time P < 0.046; treatment × time P = 0.248.
      Figure thumbnail gr4
      Figure 4Body weight response of dairy cows to concentrates containing a range of field bean (FB) inclusion levels (FB0, FB-Low, or FB-High) over the first 20 wk of lactation. SEM = 8.0; treatment P = 0.899; time P < 0.001; treatment × time P = 0.760.
      Figure thumbnail gr5
      Figure 5Effect of field bean (FB) inclusion level in the diet (FB0, FB-Low, or FB-High) on plasma chemical concentrations among dairy cows. (a) β-Hydroxybutyrate: SEM = 0.018; treatment P = 0.419; time P = 0.002; treatment × time P = 0123. (b) Glucose: SEM = 0.036: treatment P = 0.153; time P = 0.010; treatment × time P = 0.954. (c) Nonesterified fatty acids: SEM = 0.020; treatment P = 0.022; time P < 0.001; treatment × time P = 0.633. (d) Urea: SEM = 0.103; treatment P < 0.001; time P < 0.001; treatment × time P = 0.567.
      Fecal scores differed between wk 6 and 12 post-calving (overall SEM = 0.048, P = 0.012), but treatment had no effect on mean fecal scores (probability of having a score of 1, 2, or ≥ 3: 0.08, 0.47, and 0.45 for FB0; 0.05, 0.42, and 0.53 for FB-Low; 0.07, 0.45, and 0.48 for FB-High: overall SEM = 0.061, P = 0.878). We found no interaction between treatment and measurement period (SEM = 0.084, P = 0.783). The digestibility of ADF decreased (P = 0.035: linear, P = 0.013) from FB0 to FB-High, but none of the other digestibility coefficients were affected by FB inclusion level (Table 5). Neither total N intake, digestible N intake, N output in feces, urine, manure, or milk, nor any of the N utilization coefficients examined were affected by treatment (P > 0.05; Table 6). Treatment had no effect (P > 0.05) on GE intake, digestible energy (DE) intake, ME intake, energy output in feces, urine, CH4, heat, or milk, or any of the energy utilization coefficients examined (Table 7). Total CH4 production, CH4 as a proportion of DMI, OM intake, and milk production were unaffected by treatment (P > 0.05; Table 7). Similarly, CH4-E, as a proportion of GE or ME intake, was unaffected by treatment (P > 0.05). There was a linear trend for ME/GE (P = 0.071) to decrease with increasing FB inclusion in the diet.
      Table 5Effects of field bean (FB) inclusion level in dairy cow concentrates on DMI, milk yield, and total ration digestibility coefficients
      ItemTreatmentSEMP-value
      FB0FB-LowFB-HighTreatmentLinear
      Total DMI (kg/d)21.223.420.71.730.5370.830
      Milk yield (kg/d)30.231.031.23.000.9670.811
      Digestibility coefficients (g/g)
       DM0.7840.7750.7710.00940.6090.341
       OM0.8000.7910.7870.00980.6020.339
       Nitrogen0.7160.6940.7080.01860.5170.685
       Gross energy0.7770.7670.7640.00890.5480.303
       ADF0.732
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.671
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.644
      Means with the same superscript within a row do not differ significantly (P > 0.05).
      0.02860.0350.013
       NDF0.6990.6590.6660.02730.3400.256
      a,b Means with the same superscript within a row do not differ significantly (P > 0.05).
      Table 6Effect of field bean (FB) inclusion level in dairy cow concentrates on nitrogen utilization parameters
      ItemTreatmentSEMP-value
      FB0FB-LowFB-HighTreatmentLinear
      N intake and output (g/d)
       Total N intake58264256949.70.5600.851
       Digestible N intake42044440534.70.7290.772
       Feces N16319816417.90.3320.967
       Urine N19918218523.80.8680.868
       Manure N36238034934.430.8170.799
       Milk N13916314016.00.5050.961
      N utilization (g/g)
       Feces N/N intake0.2840.3060.2920.01320.5170.685
       Urine N/N intake0.3390.2870.3320.03160.4770.880
       Manure N/N intake0.6230.5930.6240.02760.6730.983
       Milk N/N intake0.2410.2550.2450.01750.8430.872
       Feces N/manure N0.4640.5160.4700.02850.4030.886
       Urine N/manure N0.5360.4840.5300.02850.4030.886
      Table 7Effects of field bean (FB) inclusion level in dairy cow concentrates on energy utilization parameters and methane production
      Item
      GE = gross energy; DE = digestible energy; CH4-E = methane energy.
      TreatmentSEMP-value
      FB0FB-LowFB-HighTreatmentLinear
      Energy intake and output (MJ/d)
       GE intake39242537531.70.5480.719
       Fecal energy8799888.10.5420.885
       DE intake30532728724.60.5460.610
       Methane energy2829303.00.8900.641
       Urinary energy1313131.00.8890.921
       Milk energy1021051029.60.9640.994
       ME intake26428424321.40.4460.513
       Heat production15115816515.10.8030.519
       Retained energy1121−2318.00.2400.204
      Energy utilization (MJ/MJ)
       DE/GE0.7700.7680.7650.00920.5180.292
       ME/GE0.6750.6670.6500.00880.1690.071
       Heat production/ME0.5640.5730.6790.04950.2430.136
       Milk energy/ME0.3840.3720.4190.01710.1880.182
       Retained energy/ME0.0520.055−0.0970.06210.1950.123
      CH4 production
       CH4 (g/d)51252754954.40.8900.641
       CH4/feed intake or milk yield (g/kg)
      CH4/DMI23.822.826.41.510.2700.258
      CH4/OM intake25.724.528.41.630.2810.274
      CH4/milk yield16.817.217.81.210.8500.850
      CH4/ECM yield15.515.616.70.8230.5060.298
       CH4-E/energy intake (MJ/MJ)
      CH4-E/GE intake0.0710.0690.0810.00670.2410.205
      CH4-E/ME intake0.1060.1040.1240.01030.1500.113
      1 GE = gross energy; DE = digestible energy; CH4-E = methane energy.
      Vaginal mucus scores decreased with week (2, 4, and 6) post-calving (overall SEM = 0.042, P < 0.001), but treatment had no effect on mean mucus scores (probability of having a score of 0, 1, or 2–3: 0.70, 0.13, and 0.17 for FB0; 0.71, 0.13, and 0.16 for FB-Low; 0.67, 0.16, and 0.17 for FB-High; overall SEM = 0.047, P = 0.97). Analysis revealed no interaction between treatment and measurement week (SEM = 0.072, P = 0.476). Treatment had no effect on the proportion of cows showing commencement of luteal activity before d 42 post-calving, days to commencement of luteal activity, or peak progesterone content at commencement of luteal activity (P > 0.05; Table 8). Treatment did not affect conception to first service, conception to first and second service, days to conception, or the proportion of cows pregnant at the end of the breeding season (P > 0.05). The proportion of cows with at least one incident of digestive upset, mastitis, or lameness, as well as mean locomotion score, were unaffected by treatment (Table 8; P > 0.05).
      Table 8Effects of field bean (FB) inclusion level in dairy cow concentrates on fertility parameters and dairy cow health (95% CI in parentheses)
      ItemTreatmentSEMP-value
      FB0FB-LowFB-High
      Pre-d 42
       Proportion of cows showing luteal activity0.75 (0.54–0.88)0.67 (0.44–0.83)0.79 (0.58–0.91)0.632
       Days to commencement of luteal activity22.121.622.12.250.980
       Peak progesterone content at commencement of luteal activity (ng/mL)31.326.932.52.740.353
      Conception to first service (proportion)0.45 (0.27–0.65)0.33 (0.16–0.55)0.30 (0.15–0.51)0.511
      Conception to first and second service (proportion)0.62 (0.42–0.79)0.57 (0.36–0.75)0.50 (0.31–0.69)0.681
      Days to conception84102909.50.412
      Cows pregnant at end of breeding season (proportion)0.91 (0.70–0.98)0.80 (0.57–0.92)0.74 (0.53–0.88)0.307
      Proportion of cows with at least one incidence of illness
       Digestive upset0.17 (0.06–0.44)0.35 (0.17–0.70)0.08 (0.02–0.32)0.117
       Mastitis0.10 (0.03–0.32)0.07 (0.02–0.27)0.08 (0.02–0.29)0.889
       Lameness0.29 (0.14–0.61)0.35 (0.17–0.69)0.37 (0.19–0.72)0.876
      Mean locomotion score2.52.52.47.300.593

      DISCUSSION

      Average daily FB intakes for FB-Low and FB-High were estimated as 4.1 and 8.4 kg/cow, respectively, with the latter substantially higher than FB intake levels used in any previous study. At the maximum inclusion level adopted (FB-High, 698 g of FB per kilogram of concentrate), FB replaced all of the soybean meal, rapeseed meal, maize gluten feed, and wheat in the concentrate, compared with the control treatment (FB0). Although it was not possible to achieve a common starch level with each of the concentrates (291, 309, and 338 g per kilogram of DM for FB0, FB-Low, and FB-High concentrates, respectively), all had a similar CP content, resulting in a total diet CP content of 180 g per kilogram of DM across the 3 diets. This total diet CP content was somewhat higher than planned, yet it remains typical of many grass silage–based diets offered in the United Kingdom.

      Cow Performance

      The potential negative effects of antinutritional factors within FB on DM intake are often cited by nutritionists as a reason for low inclusion levels in dairy cow diets. In addition, preference trials (
      • Hutson G.D.
      • van Mourik S.
      Food preferences of sheep.
      ) have suggested that FB may be unpalatable. However, the current study clearly demonstrated that FB intakes of up to 8.4 kg/cow per day had no negative effects on DM intakes. This is in agreement with most published studies (
      • Ingalls J.
      • McKirdy J.
      Faba bean as a substitute for soybean meal or rapeseed meal in rations for lactating cows.
      ;
      • Tufarelli V.
      • Khan R.U.
      • Laudadio V.
      Evaluating the suitability of field beans as a substitute for soybean meal in early lactating dairy cow: Production and metabolic responses.
      ;
      • Ramin M.
      • Höjer A.
      • Fogelberg F.
      • Hetta M.
      The effects of legume seeds on the lactation performance of dairy cow fed grass silage-based diets.
      ;
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ), albeit the maximum daily FB intake in these studies was 5.0 kg/cow per day, but the diets offered and concentrate feeding strategies adopted also differed considerably from those in the current study. In contrast,
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      found intakes to be reduced when FB replaced rapeseed meal in dairy cow diets. Those authors suggested that this may have been due to a reduction in NDF digestibility, a poorer amino acid profile with FB, and a “pull effect” caused by the higher milk yield with the rapeseed meal treatment. Despite reductions in ADF digestibility with FB inclusion in the current study, this had no effect on DMI.
      In common with DMI, milk yield was unaffected by FB inclusion. However, the reduction in both milk fat and milk protein concentrations with FB inclusion meant that yields of fat, protein, and fat plus protein were lower with both FB-Low and FB-High, compared with the FB0 treatment. Literature evidence concerning the effects of FB inclusion on milk production and milk composition is inconsistent. For example, both
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      and
      • Ramin M.
      • Höjer A.
      • Fogelberg F.
      • Hetta M.
      The effects of legume seeds on the lactation performance of dairy cow fed grass silage-based diets.
      , in studies involving maximum FB intakes of 3.7 and 2.6 kg of DM per day, respectively, observed reductions in milk yield when rapeseed meal was replaced by FB. Milk protein content was reduced in the former study, but milk fat content was unaffected by FB inclusion in either study. In contrast, at an FB intake level of 4.7 kg/d (
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ), neither milk yield nor milk fat content were affected by FB inclusion level, and milk protein content tended to decrease (P = 0.08). In several on-farm studies conducted in Italy (
      • Comellini M.
      • Antonella Volpelli L.
      • Pietro Lo Fiego D.
      • Scipioni R.
      Faba bean in dairy cow diet: Effect on milk production and quality.
      ;
      • Volpelli L.A.
      • Comellini M.
      • Masoero F.
      • Moschini M.
      • Fiego D.P.L.
      • Scipioni R.
      Faba beans (Vicia faba) in dairy cow diet: Effect on milk production and quality.
      ;
      • Tufarelli V.
      • Khan R.U.
      • Laudadio V.
      Evaluating the suitability of field beans as a substitute for soybean meal in early lactating dairy cow: Production and metabolic responses.
      ;
      • Volpelli L.A.
      • Comellini M.
      • Gozzi M.
      • Masoero F.
      • Moschini M.
      Pea (Pisum sativum) and faba beans (Vicia faba) in dairy cow diet: Effect on milk production and quality.
      ), the inclusion of FB in dairy cow diets had no effect on either milk yield or milk fat and protein content, although in a further study (
      • Mordenti A.
      • Merendi F.
      • Fustini M.
      • Formigoni A.
      Effects of different protein plants in cows diet on milk for Parmigiano Reggiano production.
      ) milk yield was reduced but milk fat content was increased with FB inclusion. However, FB inclusion levels in these studies were low, normally less than 2.0 kg/cow per day.
      Energy intake is a key driver of milk protein content, yet the reduction in milk protein content with the FB-Low treatment is unlikely to be driven by energy supply, given the similar intakes across all treatments and the absence of treatment effects on any of the BW or BCS score parameters. Indeed, plasma NEFA concentrations decreased with FB inclusion; however, there was a trend for end-of-study BCS to increase with FB inclusion, suggesting an improved energy balance with FB-High. Rather, it is more likely that this reduction was caused by a deficit of specific amino acids. For example, FB is lower in both lysine and methionine than are soybean meal or rapeseed meal (
      • Ewing W.N.
      The Feeds Directory, Volume 1: Commodity Products Guide.
      ), with both of these amino acids important for milk protein synthesis (
      • Schwab C.G.
      • Satter L.D.
      • Clay B.
      Response of lactating dairy cows to abomasal infusion of amino acids.
      ). When the rations were evaluated using Feed into Milk, the UK dairy cow rationing system (
      • Thomas C.
      Feed into Milk: A New Applied Feeding System for Dairy Cows.
      ), lysine supply relative to requirements were “adequate” for all 3 treatments, whereas for methionine the FB0 ration was classified as “borderline,” with both the FB-Low and FB-High rations classified as “poor.” When
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      replaced rapeseed meal with FB, total concentrations of essential amino acids and plasma concentrations of lysine and methionine were reduced; these authors suggested that this was likely partly responsible for the reduction in milk protein content and milk protein yield that they observed.
      In the current study, FB replaced or partially replaced both maize starch (less degradable than FB starch) and wheat starch (more degradable than FB starch) in the concentrate (
      • Offner A.
      • Bach A.
      • Sauvant D.
      Quantitative review of in situ starch degradation in the rumen.
      ), likely affecting rumen fermentation patterns and the associated reduction in milk fat content with the FB-High treatment. This reduction was in contrast to most previous studies; we propose that with the FB-High treatment, a reduction in rumen acetate production occurred due to the higher-concentrate starch and slightly lower fiber content, and this may have contributed to a reduction in de novo synthesis of milk fat. Possible adverse effects of diet on rumen function are supported in part by the lower ADF digestibility observed with this treatment. In addition, we postulate that altered rumen fermentation and biohydrogenation pathways, leading to inhibition of mammary lipid synthesis by specific fatty acid intermediates, as proposed by
      • Bauman D.E.
      • Griinari J.M.
      Regulation and nutritional manipulation of milk fat: Low milk fat syndrome.
      and reviewed more recently by others (e.g.,
      • Harvatine K.J.
      • Boisclaire Y.R.
      • Bauman D.E.
      Recent advances in the regulation of milk fat synthesis.
      ), may also have contributed to this reduction. In addition, some evidence suggests that dietary methionine promotes milk fat synthesis (
      • Qi H.
      • Meng C.
      • Jin X.
      • Li X.
      • Li P.
      • Gao X.
      Methionine promotes milk protein and fat synthesis and cell proliferation via the SNAT2–P13K signaling pathway in bovine mammary epithelial cells.
      ), and, as such, the likely deficit of methionine in the FB diets in the present study may have contributed to the reduction in milk fat concentrations.
      Effects of FB inclusion on the fatty acid profile of milk do not appear to have been examined previously; we found that FB inclusion in the current study increased the degree of saturation of the milk produced. This appears to have been largely driven by increasing concentrations of C16:0 in milk, with FB particularly high in C16:0 fatty acids (
      • Grela E.R.
      • Gunter K.D.
      Fatty acid composition and tocopherol content of some legume seeds.
      ), with
      • Moate P.J.
      • Chalupa W.
      • Boston R.C.
      • Lean I.J.
      Milk fatty acids II: Prediction of the production of individual fatty acids in bovine milk.
      estimating that 42% of C16:0 is derived directly from the diet. Similarly, decreasing concentrations of C18:1 with FB inclusion may reflect the higher concentration of C18:1 in soybean compared with FB (
      • Grela E.R.
      • Gunter K.D.
      Fatty acid composition and tocopherol content of some legume seeds.
      ).
      It is possible that the reduction in cow performance observed with FB inclusion could be avoided by supplementing the diet with specific limiting amino acids. Nevertheless, cow performance data need to be considered within the context of the total diet CP in the study (180 g/kg of DM), a level not unusual with early-lactation dairy cows offered grass silage–based diets in the United Kingdom and Ireland due to the often variable composition of grass silage on offer. However, it is possible that greater differences between treatments might have been observed with lower-protein diets as, in general, responses to protein source and quality tend to be greater at lower diet protein levels. Nevertheless,
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      observed no interaction between diet protein level (154 vs. 190 g/kg of DM) and milk production when rapeseed meal was partially replaced with FB.

      Nutrient Utilization

      Field beans contain several antinutritional substances, some of which are known to have negative effects on rumen function and digestion ability (
      • Newton S.D.
      • Hill G.D.
      The composition and nutritive value of field beans.
      ;
      • Dixon R.M.
      • Hosking B.J.
      Nutritional value of grain legumes for ruminants.
      ). For example, the FB variety Fuego is a colored-flower variety, which is known to be a tannin-containing type. However, fecal scores—a simple proxy for rumen function—were unaffected by treatment in this study. In addition, with the exception of ADF, no digestibility coefficients for the other parameters examined were affected by FB inclusion. In our previously published work, we found fiber digestibility to be unaffected with an FB intake similar to that used with the FB-Low treatment (
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ), but the reduction in ADF digestibility with the highest FB treatment in the current study is similar to the reduction in NDF digestibility observed by
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      .
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      attributed this to the readily fermentable FB carbohydrate and associated fermentation products, decreasing rumen pH and having a negative effect on cellulolytic bacteria activity, rather than to the presence of antinutritional factors.
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      also found the apparent digestibility of DM, OM, N, and starch to increase with FB inclusion, but they attributed this to the lower proportion of silage in the diets of cows offered FB, rather than to FB inclusion per se. Thus it appears that any effects of FB inclusion on digestibility are likely to be caused by other dietary factors, rather than by the presence of antinutritional factors. This is supported by the findings of
      • Melicharová V.
      • Pechová A.
      • Dvořák R.
      • Pavlata L.
      Performance and metabolism of dairy cows fed bean seeds (Vicia faba) with different levels of anti-nutritional substances.
      , who found cow performance to be unaffected when FB varieties containing either high or low levels of ANF were offered.
      Although our earlier work (
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ) indicated a reduction in the ratio of milk N to N intake with FB inclusion, a similar effect was not observed in the nutrient utilization portion of the current study, in agreement with other previous findings (
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      ;
      • Ramin M.
      • Höjer A.
      • Fogelberg F.
      • Hetta M.
      The effects of legume seeds on the lactation performance of dairy cow fed grass silage-based diets.
      ). Nevertheless, the nutrient utilization phase in the current study involved only 4 cows per treatment and was undertaken at the end of the 20-wk feeding study. When N use efficiency over the entire study period is examined using treatment mean intake and milk production data (with milk N calculated as milk protein/6.38), the ratio of milk N to N intake was 0.30, 0.28, and 0.27 for FB-0, FB-Low, and FB-High, respectively. The decreasing N use efficiency with increasing FB inclusion levels was reflected in an increase in plasma urea N levels, in agreement with
      • Volpelli L.A.
      • Comellini M.
      • Gozzi M.
      • Masoero F.
      • Moschini M.
      Pea (Pisum sativum) and faba beans (Vicia faba) in dairy cow diet: Effect on milk production and quality.
      , the trend (P = 0.078) observed in our previous work (
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ), and the increase in milk urea concentrations observed by
      • Puhakka L.
      • Jaakkola S.
      • Simpura L.
      • Kokkonen T.
      • Vanhatalo A.
      Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
      when the rapeseed component of the diet was replaced by FB. In contrast, other studies (
      • Comellini M.
      • Antonella Volpelli L.
      • Pietro Lo Fiego D.
      • Scipioni R.
      Faba bean in dairy cow diet: Effect on milk production and quality.
      ;
      • Volpelli L.A.
      • Comellini M.
      • Masoero F.
      • Moschini M.
      • Fiego D.P.L.
      • Scipioni R.
      Faba beans (Vicia faba) in dairy cow diet: Effect on milk production and quality.
      ;
      • Tufarelli V.
      • Khan R.U.
      • Laudadio V.
      Evaluating the suitability of field beans as a substitute for soybean meal in early lactating dairy cow: Production and metabolic responses.
      ) have found both blood and milk urea N concentrations to be either reduced or unchanged with FB inclusion. Protein in FB has a higher degradability in the rumen compared with protein from soybean meal or rapeseed meal (
      • Ewing W.N.
      The Feeds Directory, Volume 1: Commodity Products Guide.
      ), and in general blood and milk urea levels increase when FB replaces a less-degradable protein source.
      None of the energy utilization parameters examined were affected by FB inclusion level, which was largely in agreement with our earlier study (
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ), except that these latter authors inexplicably found a quadratic effect of FB inclusion on urine E output. Similarly, none of the CH4 production parameters examined were affected by FB inclusion, in agreement with previous indications (
      • Johnston D.J.
      • Theodoridou K.
      • Ferris C.P.
      The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
      ). However,
      • Ramin M.
      • Höjer A.
      • Fogelberg F.
      • Hetta M.
      The effects of legume seeds on the lactation performance of dairy cow fed grass silage-based diets.
      found that ratios of CH4 to DMI and CH4 to ECM were unaffected when rapeseed meal was replaced by FB, whereas total CH4 production tended to be reduced, which they attributed in part to the higher fat content of rapeseed meal. Despite evidence suggesting that CH4 production may be reduced when ruminants are offered legume forages (
      • Lüscher A.
      • Mueller-Harvey L.
      • Soussana J.F.
      • Rees R.M.
      • Peyraud J.L.
      Potential of legume-based grassland–livestock systems in Europe: A review.
      ), there appears to be little evidence that this is the case with grain legumes, in line with the findings of the current study. Our results demonstrate enteric CH4 emissions to be unaffected by FB inclusion, but the effects of replacing imported protein ingredients with locally grown FB, across the life cycle of milk, needs to be examined.

      Fertility and Health

      Negative energy balance can have a detrimental effect on reproductive performance (
      • Beam S.W.
      • Butler W.R.
      Effects of energy balance on follicular development and first ovulation in postpartum dairy cows.
      ;
      • Butler W.R.
      Energy balance relationships with follicular development, ovulation and fertility in postpartum dairy cows.
      ), but BCS and BW data in the current study suggest that energy balance differed little between treatments. However, the presence of phytoestrogens in many legumes species (for example, daidzein and genistein;
      • Kaufman P.B.
      • Duke J.A.
      • Brielmann H.
      • Boik J.
      • Hoyt J.E.
      A comparative survey of leguminous plants as sources of the isoflavones, genistein and daidzein: Implications for human nutrition and health.
      ), including some FB cultivars, is often highlighted as a cause for concern in relation to dairy cow fertility. Phytoestrogens, compounds that are produced naturally in many plant species, may have adverse effects on ovarian function, including possibly inhibited luteinizing hormone surge, altered development of or lack of dominant follicles, abnormal follicular waves, early embryonic death, and repeat breeding (
      • Mostrom M.
      • Evans T.J.
      Phytoestrogens.
      ). Although these compounds have been shown to suppress the secretion of luteinizing hormone (
      • Zdunczyk S.
      • Piskula M.
      • Janowski T.
      • Baranski W.
      • Ras M.
      Concentrations of isoflavones in blood plasma of dairy cows with different incidence of silent heat.
      ) and to lower progesterone concentrations (
      • Piotrowska K.K.
      • Woclawek-Potocka L.
      • Bah M.M.
      • Piskula M.K.
      • Pilawski W.
      • Bober A.
      • Skarzynski D.J.
      Phytoestrogens and their metabolites inhibit the sensitivity of the bovine corpus luteum to luteotropic factors.
      ), neither the time to onset of cyclicity, the proportion of cows cycling within 42 d of calving, nor the peak progesterone concentrations during the first estrus cycle were affected by diet in the present study. While actual conception rates within the study must be considered within the context of the number of cows on the study, we found no evidence that reproductive outcomes were affected by FB inclusion level. Vaginal mucus scores, assessed as an indicator of immune status, did not differ across treatments, suggesting that FB treatment has no effect on immune status. In addition, none of the health parameters examined were affected by FB inclusion in the diet, which again suggests that immune function was unaffected by treatment.

      CONCLUSIONS

      Including FB (Vicia faba. var. Fuego) in dairy cow diets at up to 8.4 kg/cow per day had no detrimental effects on feed intake, milk production, body tissue reserves, fertility performance, cow health, or diet digestibility. The reduction in milk fat content with the highest FB treatment was likely due to the effect on rumen function of the high starch levels with this treatment, and the reduction in milk protein levels with FB inclusion was likely due to a deficit in specific limiting amino acids. Methane production was unaffected by FB inclusion.

      ACKNOWLEDGMENTS

      This study was funded by the Department of Agriculture, Environment, and Rural Affairs (DAERA) in Belfast, Northern Ireland, and by AgriSearch, Innovation Centre, Hillsborough, Northern Ireland. David Johnston was in receipt of a DAERA Postgraduate Studentship. The authors thank the staff of the Dairy Unit at the Agri-Food and Biosciences Institute (Hillsborough, Northern Ireland) for management of the experimental animals, and the laboratory staff for analysis of samples. The invaluable help provided by Jim Uprichard (Trouw Nutrition, Belfast, UK) in formulating the experimental concentrates was much appreciated.

      REFERENCES

        • Abdulla J.M.
        The use of field beans as a feed for broiler chickens. PhD Thesis.
        Harper Adams University, UK2017
        • Bauman D.E.
        • Griinari J.M.
        Regulation and nutritional manipulation of milk fat: Low milk fat syndrome.
        Livest. Prod. Sci. 2001; 70: 15-29
        • Beam S.W.
        • Butler W.R.
        Effects of energy balance on follicular development and first ovulation in postpartum dairy cows.
        J. Reprod. Fertil. Suppl. 1999; 54 (10692872): 411-424
        • Bligh E.G.
        • Dyer W.J.
        A rapid method of total lipid extraction and purification.
        Can. J. Biochem. Physiol. 1959; 37 (13671378): 911-917
        • Butler W.R.
        Energy balance relationships with follicular development, ovulation and fertility in postpartum dairy cows.
        Livest. Prod. Sci. 2003; 83: 211-218
        • Comellini M.
        • Antonella Volpelli L.
        • Pietro Lo Fiego D.
        • Scipioni R.
        Faba bean in dairy cow diet: Effect on milk production and quality.
        Ital. J. Anim. Sci. 2009; 8: 396-398
        • Dixon R.M.
        • Hosking B.J.
        Nutritional value of grain legumes for ruminants.
        Nutr. Res. Rev. 1992; 5 (19094311): 19-43
        • Dvořák R.
        • Pechová A.
        • Pavlata L.
        • Klejdus B.
        • Kovařčík K.
        • Dostálová J.
        • Culková J.
        • Filípek J.
        • Svajdlenka E.
        • Capkova V.
        Reduction in the content of antinutritional substances in Fava beans (Vicia faba) by different treatments.
        Slov. Vet. Zb. 2006; 43: 174-179
        • 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
        • Ewing W.N.
        The Feeds Directory, Volume 1: Commodity Products Guide.
        Context Publications, Leicestershire, UK1997
        • Flores F.
        • Hybl M.
        • Knudsen J.C.
        • Marget P.
        • Muel F.
        • Nadal S.
        • Narits L.
        • Raffiot B.
        • Sass O.
        • Solis I.
        • Winkler J.
        Adaptation of spring faba bean types across European climates.
        Field Crops Res. 2013; 145: 1-9
        • Grela E.R.
        • Gunter K.D.
        Fatty acid composition and tocopherol content of some legume seeds.
        Anim. Feed Sci. Technol. 1995; 52: 325-331
        • Qi H.
        • Meng C.
        • Jin X.
        • Li X.
        • Li P.
        • Gao X.
        Methionine promotes milk protein and fat synthesis and cell proliferation via the SNAT2–P13K signaling pathway in bovine mammary epithelial cells.
        J. Agric. Food Chem. 2018; 66 (30274521): 11027-11033
        • Harvatine K.J.
        • Boisclaire Y.R.
        • Bauman D.E.
        Recent advances in the regulation of milk fat synthesis.
        Animal. 2009; 3 (22444171): 40-54
        • Hutson G.D.
        • van Mourik S.
        Food preferences of sheep.
        Aust. J. Exp. Agric. 1981; 21: 575-582
        • Ingalls J.
        • McKirdy J.
        Faba bean as a substitute for soybean meal or rapeseed meal in rations for lactating cows.
        Can. J. Anim. Sci. 1974; 54: 87-89
        • Ireland-Perry R.L.
        • Stallings C.C.
        Fecal consistency as related to dietary composition in lactating Holstein cows.
        J. Dairy Sci. 1993; 76 (8387550): 1074-1082
        • Johnston D.J.
        • Theodoridou K.
        • Ferris C.P.
        The impact of field bean inclusion level in dairy cow diets on cow performance and nutrient utilisation.
        Livest. Sci. 2019; 220: 166-172
        • Kaufman P.B.
        • Duke J.A.
        • Brielmann H.
        • Boik J.
        • Hoyt J.E.
        A comparative survey of leguminous plants as sources of the isoflavones, genistein and daidzein: Implications for human nutrition and health.
        J. Altern. Complement. Med. 1997; 3 (9395689): 7-12
        • Kudlinskienė I.
        • Lionikienė J.
        • Klementavičiūtė J.
        • Stanytė G.
        • Gružauskas R.
        • Stankevičius R.
        • Urbšienė D.
        • Bliznikas S.
        • Racevičiūtė-Stupelienė A.
        • Ots M.
        • Kass M.
        Effect of extruded field beans “Fuego” (Vicia Faba) on dairy cow's performance and milk sensory properties.
        Vet Med Zoot. T. 2018; 76
        • Licitra G.
        • Hernandez T.M.
        • Van Soest P.J.
        Standardization of procedures for nitrogen fractionation of ruminant feeds.
        Anim. Feed Sci. Technol. 1996; 57: 347-358
        • Little M.W.
        • O'Connell N.E.
        • Welsh M.D.
        • Mulligan F.J.
        • Ferris C.P.
        Concentrate supplementation of a diet based on medium-quality grass silage for 4 weeks prepartum: Effects on cow performance, health, metabolic status, and immune function.
        J. Dairy Sci. 2017; 100 (28342599): 4457-4474
        • Lüscher A.
        • Mueller-Harvey L.
        • Soussana J.F.
        • Rees R.M.
        • Peyraud J.L.
        Potential of legume-based grassland–livestock systems in Europe: A review.
        Grass Forage Sci. 2014; 69 (26300574): 206-228
        • Melicharová V.
        • Pechová A.
        • Dvořák R.
        • Pavlata L.
        Performance and metabolism of dairy cows fed bean seeds (Vicia faba) with different levels of anti-nutritional substances.
        Acta Vet. Brno. 2009; 78: 57-66
        • Moate P.J.
        • Chalupa W.
        • Boston R.C.
        • Lean I.J.
        Milk fatty acids II: Prediction of the production of individual fatty acids in bovine milk.
        J. Dairy Sci. 2008; 91 (18292274): 1175-1188
        • Mordenti A.
        • Merendi F.
        • Fustini M.
        • Formigoni A.
        Effects of different protein plants in cows diet on milk for Parmigiano Reggiano production.
        Ital. J. Anim. Sci. 2007; 6: 463-465
        • Mostrom M.
        • Evans T.J.
        Phytoestrogens.
        in: Gupta R.C. Veterinary Toxicology: Basic and Clinical Principles. 3rd ed. Academic Press, Cambridge, MA2018: 817-833
        • Newton S.D.
        • Hill G.D.
        The composition and nutritive value of field beans.
        Nutr. Abstr. Rev. Series B. 1983; 53: 99-115
        • Offner A.
        • Bach A.
        • Sauvant D.
        Quantitative review of in situ starch degradation in the rumen.
        Anim. Feed Sci. Technol. 2003; 106: 81-93
        • Park R.S.
        • Agnew R.E.
        • Gordon F.J.
        • Steen R.W.J.
        The use of near infrared reflectance spectroscopy (NIRS) on undried samples of grass silage to predict chemical composition and digestibility parameters.
        Anim. Feed Sci. Technol. 1998; 72: 155-167
        • Piotrowska K.K.
        • Woclawek-Potocka L.
        • Bah M.M.
        • Piskula M.K.
        • Pilawski W.
        • Bober A.
        • Skarzynski D.J.
        Phytoestrogens and their metabolites inhibit the sensitivity of the bovine corpus luteum to luteotropic factors.
        J. Reprod. Dev. 2006; 52 (16276041): 33-41
        • Puhakka L.
        • Jaakkola S.
        • Simpura L.
        • Kokkonen T.
        • Vanhatalo A.
        Effects of replacing rapeseed meal with fava bean at 2 concentrate crude protein levels on feed intake, nutrient digestion, and milk production in cows fed grass silage–based diets.
        J. Dairy Sci. 2016; 99 (27522411): 7993-8006
        • Purcell P.J.
        • Law R.A.
        • Gordon A.W.
        • McGettrick S.A.
        • Ferris C.P.
        Effect of concentrate feeding method on the performance of dairy cows in early to mid lactation.
        J. Dairy Sci. 2016; 99 (26805998): 2811-2824
        • Ramin M.
        • Höjer A.
        • Fogelberg F.
        • Hetta M.
        The effects of legume seeds on the lactation performance of dairy cow fed grass silage-based diets.
        Agric. Food Sci. 2017; 26: 129-137
        • Schwab C.G.
        • Satter L.D.
        • Clay B.
        Response of lactating dairy cows to abomasal infusion of amino acids.
        J. Dairy Sci. 1976; 59 (950397): 1254-1270
        • Thomas C.
        Feed into Milk: A New Applied Feeding System for Dairy Cows.
        Nottingham University Press, Nottingham, UK2004
        • Tufarelli V.
        • Khan R.U.
        • Laudadio V.
        Evaluating the suitability of field beans as a substitute for soybean meal in early lactating dairy cow: Production and metabolic responses.
        Anim. Sci. J. 2012; 83 (22339694): 136-140
        • Volpelli L.A.
        • Comellini M.
        • Gozzi M.
        • Masoero F.
        • Moschini M.
        Pea (Pisum sativum) and faba beans (Vicia faba) in dairy cow diet: Effect on milk production and quality.
        Ital. J. Anim. Sci. 2012; 11: 217-222
        • Volpelli L.A.
        • Comellini M.
        • Masoero F.
        • Moschini M.
        • Fiego D.P.L.
        • Scipioni R.
        Faba beans (Vicia faba) in dairy cow diet: Effect on milk production and quality.
        Ital. J. Anim. Sci. 2010; 9: 138-144
        • Watson C.
        • Reckling M.
        • Preissel S.
        • Bachinger J.
        • Bergkvist G.
        • Kuhlman T.
        • Lindstrom K.
        • Nemecek T.
        • Topp C.F.E.
        • Vanhatalo A.
        • Zander P.
        • Murphy-Bokern D.
        • Stoddard F.L.
        Grain legume production and use in European agricultural systems.
        Adv. Agron. 2017; 144: 235-303
        • Wilkins R.J.
        • Jones R.
        Alternative home-grown protein sources for ruminants in the United Kingdom.
        Anim. Feed Sci. Technol. 2000; 85: 23-32
        • Zdunczyk S.
        • Piskula M.
        • Janowski T.
        • Baranski W.
        • Ras M.
        Concentrations of isoflavones in blood plasma of dairy cows with different incidence of silent heat.
        Bull. Vet. Inst. Pulawy. 2005; 49: 189-191