Next Article in Journal
Characterisation of the Behavioural Effects of a Thoracic Squeeze in Healthy Newborn Piglets
Previous Article in Journal
Canine Olfaction: Physiology, Behavior, and Possibilities for Practical Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nitrogen Balance of Dairy Cows Divergent for Milk Urea Nitrogen Breeding Values Consuming Either Plantain or Perennial Ryegrass

1
Faculty of Agriculture and Life Sciences, Lincoln University, P.O. Box 85084, Lincoln 7647, New Zealand
2
USDA-ARS, Livestock Nutrient Management Research Unit, 300 Simmons Drive, Unit 10, Bushland, TX 79012, USA
*
Author to whom correspondence should be addressed.
Animals 2021, 11(8), 2464; https://doi.org/10.3390/ani11082464
Submission received: 19 July 2021 / Revised: 19 August 2021 / Accepted: 20 August 2021 / Published: 22 August 2021
(This article belongs to the Section Animal Nutrition)

Abstract

:

Simple Summary

We studied the nitrogen excretion patterns of cows selected for divergent nitrogen excretion consuming either ryegrass or plantain. Both the use of a plantain diet as well as the use of low milk urea nitrogen breeding values were found to reduce the concentration of urinary urea nitrogen per urination event compared to cows with a high milk urea nitrogen breeding value and cows consuming a ryegrass-based diet. These results indicate that both the use of cows with low milk urea nitrogen breeding values and the use of a plantain diet are tools that temperate pastoral dairy production systems can use to reduce nitrogen loses.

Abstract

Inefficient nitrogen (N) use from pastoral dairy production systems has resulted in environmental degradation, as a result of excessive concentrations of urinary N excretion leaching into waterways and N2O emissions from urination events into the atmosphere. The objectives of this study were to measure and evaluate the total N balance of lactating dairy cows selected for milk urea N concentration breeding values (MUNBVs) consuming either a 100% perennial ryegrass (Lolium perenne L.) or 100% plantain (Plantago lanceolata L.) diet. Sixteen multiparous lactating Holstein-Friesian × Jersey cows divergent for MUNBV were housed in metabolism crates for 72 h, where intake and excretions were collected and measured. No effect of MUNBV was detected for total N excretion; however, different excretion characteristics were detected, per urination event. Low MUNBV cows had a 28% reduction in the concentration of urinary urea nitrogen (g/event) compared to high MUNBV cows when consuming a ryegrass diet. Cows consuming plantain regardless of their MUNBV value had a 62% and 48% reduction in urinary urea nitrogen (g/event) compared to high and low MUNBV cows consuming ryegrass, respectively. Cows consuming plantain also partitioned more N into faeces. These results suggest that breeding for low MUNBV cows on ryegrass diets and the use of a plantain diet will reduce urinary urea nitrogen loading rates and therefore estimated nitrate leaching values, thus reducing the environmental impact of pastoral dairy production systems.

1. Introduction

Nitrate leaching is a key environmental concern from temperate pastoral dairy production systems, largely as a result of an oversupply of dietary nitrogen (N) [1]. Dairy cows typically only utilise 30% of dietary N for milk production, with approximately 70% of ingested N being excreted as a waste product [2]. Over 60% of this surplus N is excreted as urinary N (UN) [3,4], which is in an immediately available form for leaching or volatilisation [1]. Approximately 82% of UN is discharged to pastures in temperate pasture-based systems [5,6], with 20–30% of this UN being leached to groundwater as nitrates [7] and 2% transformed to nitrous oxide (N2O) [8]. Excessive concentrations of N in waterways are linked to environmental degradation, such as eutrophication of waterways, while N2O is a potent greenhouse gas (GHG) [9]. This environmental degradation reduces the ability of countries where temperate pastoral dairy production practises are used, such as New Zealand, to meet their commitments under global conventions, such as the Kyoto protocol and Ramsar convention [10], and reduces their ability to meet ethical [11] and sustainable development goals set out by the United Nations [12]. There has also been an association between high concentrations of nitrates in drinking water and several negative human health outcomes, such as increased risk of colorectal cancer [13], methemoglobinemia in infants [14], and neural tube defects [15].
Dietary and management strategies focusing on diverse swards and grazing management have been extensively explored as ways to reduce N losses [4,16,17,18,19], with the use of narrow-leaf plantain (Plantago lanceolata L. (PL)) showing promising results for reducing total N excretions and UN concentrations in dairy cows [20,21]. PL is known to contain secondary bioactive compounds, such as acubin, catalpol, and acteoside [22,23], that have been reported to have antimicrobial and diuretic effects [22,24], which may influence the N partitioning of the animals. The use of animal genetics to reduce environmental impact is being explored. Several studies have proposed the potential of using milk urea nitrogen concentration (MUN) values as a predictor for urinary urea excretions, as urea is known to equilibrate between the body fluids of an animal [25], with a linear relationship between the blood and urinary urea pools being established [26,27,28]. The heritability of MUN has been well documented, with estimates ranging between 0.14 and 0.24 [29,30,31]. Beatson et al. (2019) [31] proposed that breeding for lower MUN may present an opportunity to reduce urinary N excretion. Marshall et al. (2020) [32] demonstrated a relationship between MUN breeding values (MUNBVs) and urinary urea nitrogen (UUN) excretion, showing that in a grazing setting, dairy cows with a low MUNBV had reductions in UUN excretion, thus indicating a potential for animal-based solutions to reduce nitrogen losses.
While the relationship between MUNBV and UUN excretion has been demonstrated in a grazing environment, to our knowledge, this relationship has not been conducted in a controlled environment with all urination, defecation, and milk production events being collected and collated into a total daily N balance. Any interaction between diet and MUNBV is also yet to be reported.
We hypothesized that cows divergent for MUNBVs will have differences in their nitrogen balance, with the effect of diet being additive to MUNBV for reducing UUN excretion. The objectives of this study were, therefore, to capture and collate individual intake, urination, defecation, and milk production events into a total N balance of lactating dairy cows divergent for MUNBV; and secondly, to evaluate the effect of a monoculture diet of PL compared with a monoculture diet of perennial ryegrass (Lolium perenne L. (RG)) diet on relationships between MUNBVs and milk production and excretion parameters.

2. Materials and Methods

This study was conducted at Lincoln University’s Johnstone Memorial Laboratory in Canterbury, New Zealand (43°64′ S, 172°45′ E) during January and February 2020. Sixteen multiparous lactating Holstein-Friesian × Jersey cows (150.4 ± 20.7 days in milk; 506 ± 35 kg live weight; 3.75 ± 0.25 body condition score; 6.7 ± 1.4 years old; 4.6 ± 1.4 lactations) divergent for their estimated MUNBV mg/dL were used. Milk urea nitrogen breeding values were estimated as described by Beatson et al. (2019) [31] and were provided as part of a regular commercial herd test conducted by CRV (New Zealand). Cows were selected based on their level of divergence from the herd’s mean value: the cows used during this study had an MUNBV range from −1.4 to +3.3.
A 2 × 2 factorial design with repeated measures was implemented over 4 runs. Each run consisted of 12 days (9 acclimation days, 3 measurement days in metabolism crates). Each run was conducted concurrently with a 3-day delay so that as the first group of cows left the metabolism crates, the next group entered. During each run, cows were grouped as either high or low for MUNBV (8 high vs. 8 low) and diet type (RG (cv. One50) or PL (cv. Agritonic)). Four cows classified as high for MUNBV and 4 cows classified as low for MUNBV were allocated to each of the 2 dietary treatments (PL or RG), resulting in 1 animal per factorial arrangement (e.g., 1 high PL, 1 high RG, 1 low PL, 1 low RG) per run. The average MUNBV for cows classified as ‘high’ in the PL treatment was 2.19 ± 0.21 and −1.23 ± 0.21 for cows classified as low. Cows in the RG treatment had an average MUNBV of 2.23 ± 0.21 for cows classified as high and −1.09 ± 0.21 for those classified as low. Cows were balanced for production worth, age, and live weight across both dietary treatments and MUNBV groups.
All cows grazed a ryegrass-based sward before the start of the experiment and were all dosed with Rumensin capsules (Elanco, Auckland, New Zealand) 1 week before the start of experimentation to help prevent bloat. A nine-day indoor acclimation period was used to acclimate cows to their allocated diets and feeding bins before cows were housed in metabolism crates for 72 h. The plantain proportion in the diet was increased by 20% each day starting at a baseline level of 20% during the acclimation period, with the remaining proportion of the diet being ryegrass. In this way, cows were on a 100% diet of plantain by day 5. Cows had ad libitum access to fresh water during both the acclimation and measurement periods. During the measurement period, PIUSI K24 (PIUSI, Via pacinotti, Italy) flow rate meters were installed in the troughs in the metabolism crates to measure water intake.

2.1. Herbage Measurements

Fresh herbage was cut daily using a Haldrup (GmbH, Ilshofen, Germany) forage harvester at a height of 3 cm from a ~30 cm height standing sward before being stored in a chiller at 4 °C and fed non-chopped within 24 h of harvest. Twice daily, after morning and afternoon milking (0700 and 1500 h, respectively), refusals from the previous period were removed and weighed before new herbage was weighed and allocated to the cows. A 10% refusal of as fed dry matter (DM) per feeding for each cow was targeted to provide ad libitum access to feed. At each feeding, samples of allocated and refused herbages were taken for analysis of DM content, and chemical and botanical composition.
Botanical composition was determined for both allocated and refused feed by splitting a representative sample into its constituent parts before being oven-dried for 72 h (60 °C). The chemical composition of both allocated and refused feed was determined for each feeding by taking a representative sample from each dietary treatment. Subsamples were lyophilized and ground using a 1-mm sieve in a centrifugal mill (ZM200 Retsch, Haan, Germany) and analysed for quality attributes. Near-infrared spectrophotometry (NIRS, Model: FOSS NIRSystems 5000, Laurel, MD, USA) was used to analyse organic matter (OM), water-soluble carbohydrates (WSCs), neutral and acid detergent fibre (NDF, ADF), crude protein (CP), and DM, OM, and dry OM digestibility (DMD, OMD, DOMD, respectively). The metabolizable energy (ME) content of forage was calculated as per the calculation of the primary Industries Standing Committee (2007) [33] of 0.16 × DOMD. NIRS calibration equations were determined prior to sample analysis for OM by digestibility (R2 = 0.89) [34], WSC (R2 = 0.97) [35], ADF (R2 = 0.97) (AOAC, 1990 [36]; method 973.18), NDF (R2 = 0.98) [37], CP by combustion (R2 = 0.99) (Vario Max CN Analyser Elementar, Langenselbold, Germany), DMD (R2 = 0.97), OMD (R2 = 0.97), and DOMD (R2 = 0.92) [34].

2.2. Animal Measurements

Cows were monitored for the entire duration that they were housed in the metabolism crates by trained technicians who were unaware of each animal’s MUNBV. Cows were fitted with a 3-D-printed attachment constructed of acrylonitrile butadiene styrene printed using an UP-Mini 3-D printer (3D Printing Systems, Auckland, NZ), which fitted over the vulva of the animal and was secured with the use of a biocompatible glue (Loctite 454 was used; Henkel, Düsseldorf, Germany). Attached to the 3-D-printed attachment was approximately 1 m of plastic vinyl sleeve, which ‘hung free’ behind the animal and allowed the flow of a urination event to be diverted into a bucket by a technician. In this manner, the urination pattern was able to be recorded per event with every event being individually captured in a bucket, with the time and weight recorded. A subsample of each urination event was immediately acidified with 0.1 mL of sulfuric acid (concentration ≥95%) from a 1-mL transfer pipette and capped to minimise ammonia (NH3) volatilisation before being frozen. Total collection trays were in place to capture any urine that may have been missed or spilt. These trays were emptied and weighed at both morning and afternoon milkings. Urine from the total collection trays was not used for chemical analysis due to the possibility of faecal contamination. The volume of urine missed from the bucket was calculated to be consistent per event, with the weight of the urine in the tray being divided by the number of events that occurred during these periods and added to the volume of each event to create an adjusted weight. On average, 30 mL were missed per urination event. A commercial enzymatic kinetic technique using a clinical analyser (Randox Daytona; Crumlin, UK) was used to determine UUN concentration. Urine total N content was determined by combustion (Vario Max CN Analyser Elementar, Langenselbold, Germany) from samples that were composited by volume depending on the relative urination event size from each 24-h period that the cows were housed in the metabolism crates.
Nitrogen deposited to pasture was calculated as the N loading rate using Equation (1) of Haynes and Williams (1993) [38] as previously implemented by Marshall et al. (2020) [32]. The loading rate of UUN per urine patch can be estimated as:
UUN   rate   ( kg   UUN   ha 1 ) = Concentration   ( g   UUN   L 1 ) × Vol   ( L ) Area   per   urination   event   ( m 2 ) × 10
The area covered by the urine patch was estimated using the following equation of Romera et al. (2012) [39] (2) assuming a urine column size of 5 mm:
Area   per   urination   event = Average   volume   per   event   × 1000 10,000 Urine   column × 0.1
Purine derivatives were determined from an average weighted urine sample for each 24-h period cows were housed in the metabolism creates. A modified HPLC methodology [40] was used on an Agilent 1100 series chromatograph (Agilent Technologies, Santa Clara, CA, USA) with a prodigy 250 × 4.6 mm, 5 u, ODS, C18 reverse-phase column (phenomenex, Torrance, CA, USA). Total purine excretion was calculated as the sum of allantoin and uric acid [41] multiplied by the total volume of urine per day. Xanthine and hypoxanthine values were negligible and therefore not included in total purine derivatives; these purines are often not detected in cattle urine [42,43]. Total purine absorbed was then calculated using the equation of International Atomic Energy Agency (1997) [41] (3):
y = 0.85 x + ( 0.385   L W 0.75 )
where y = PD excretion and x = purine absorption (mmol/d) and LW 0.75 = metabolic liveweight of the animal (kg).
Microbial CP supply was then estimated using the following Equation (4), where χ = purine absorbed:
M i c r o b i a l N ( g N / d ) = χ ( m m o l / d ) × 70 0.116 × 0.83 × 1000 = 0.727 χ
The N content of purines was assumed to be 70 mg N/mmol, the ratio of purine N to total N in mixed rumen microbes was 11.6:100, and the digestibility of microbial purines was assumed to be 0.83 [41]. Microbial CP was then calculated by multiplying microbial N by 6.25. Microbial synthesis efficiency was calculated as grams of urinary microbial CP per kilogram of truly digestible OM intake [42,44].
Faecal outputs were recorded in the same manner as urine outputs, with each event being collected in a bucket by a technician with the weight and time being recorded. A subsample was taken per event and frozen for analysis. A screen above the urine total collection tray was in place to capture any missed faecal material per event. Technicians were able to collect all faecal material per event as any missed material was collected from the screen. Faecal samples were lyophilized to determine faecal DM before being ground through a 1-mm sieve (ZM200 Retsch, Haan, Germany) and analysed for OM, NDF, ADF, and N percentage by NIRS (Model: FOSS NIRSystems 5000, Laurel, MD, USA). NIRS calibration equations were determined for faeces prior to sample analysis for OM by digestibility (R2 = 0.93) [34], ADF (R2 = 0.99) (AOAC [36], 1990; method 973.18), NDF (R2 = 0.99) [37], N% by combustion (R2 = 0.98) (Vario Max CN Analyser Elementar, Langenselbold, Germany), and DMD (R2 = 0.96) [34].
Cows were milked twice daily at 07:00 and 15:00 h throughout the adaptation and measurement period. A portable milking machine (DeLaval vacuum pump DVP170; DeLaval: Tumba, Sweden) consisting of two separate lines and cisterns was used so that each cow’s milk could be kept separate, with samples collected during the 3 days of the measurement period assessed for yield and quality attributes. Lines were flushed with 20 L of cold water and ‘sucked’ dry between cows to avoid any cross-contamination of milk between cows. Representative 100-mL milk samples were taken per milking and analysed by MilkTest, New Zealand using a CombiFoss machine (Foss Electric, Hillerød, Denmark) for quality attributes (protein, fat, lactose, and somatic cell count) and MUN was measured on skimmed milk by an automated Modular P analyser (Roche/Hitachi; [45]). Milk DM was determined after drying a sample for 5 days at 60 °C.
Rumen samples were collected from each animal by oesophageal tubation on the first day of the adaption period and the last day of the measurement period (days 1 and 12). Volatile fatty acid (VFA) concentrations were determined by gas chromatography [46] using a gas chromatographer (GC-2010, Shimadzu, Kyoto, Japan) fitted with an SGE BP21 30 m × 530 μm × 1.0 μm bore capillary column, VFAs are grouped and reported as glucogenic and non-glucogenic VFAs. NH3 concentration was determined on a rumen fluid sub-sample acidified using sulphuric acid by an enzymatic UV method [47], using an automated clinical analyser (Randox Rx Daytona, Crumlin, UK).

2.3. Statistical Analysis

Data were tested for normality using Shapiro–Wilk’s test to satisfy the assumption of normality for the statistical models. Data that were normally distributed were analysed using a mixed model ANOVA using the ‘lme4′ package, and a type III analysis of variance table using Satterthwaite’s method for determining degrees of freedom was used to determine p values for fixed effects and interactions. Non-normally distributed data were analysed using a generalised linear mixed model ANOVA using the ‘lme4’ package, and a type II analysis of variance table using Wald chi-square method was used for determining p values for fixed effects and interactions. Milk urea nitrogen breeding value group (high or low), diet (RG or PL), and MUNBV group × diet interactions were used as fixed effects for all analysis. Day (1 and 12 (first and last day of the experiment)) was used as a fixed effect for rumen fermentation characteristics, and its interaction terms with diet and MUNBV groups. The day of treatment in the metabolism crate (1, 2, or 3) was nested within run (1, 2, 3, and 4) and run itself was fitted as a random effects for analyses involving faecal and urinary characteristics, water and N balance, digestibility, purine, microbial efficiency, and milk characteristic analysis. Run + date of the trial was fitted as a random effect for intake analysis with individual cow’s milk production breeding value fitted as a covariate. The milk production breeding value was found to have no effect and was thus removed from the final model. Run was included as a random effect for VFA analysis, and the date of the trial was included as a random effect for botanical analysis. Time of day (am or pm) nested within date was used as the random effect for herbage quality analysis. Herbage chemical composition values were analysed for every day animals were housed within metabolism crates and then reported as an average value for the duration of the trial. Statistical significance was determined at p ≤ 0.05 and tendencies are discussed at 0.05 < p ≤ 0.10. All statistical analyses were conducted using R [48].

3. Results

3.1. Herbage Measurements

Ryegrass diets contained 13 ± 6% weed, 4 ± 0.8% dead material, 3 ± 3% reproductive stem, and 80 ± 4% leaf. The PL diet was comprised of 9 ± 5% weed, 2 ± 1% dead material, 17 ± 3% reproductive stem, and 73 ± 5% leaf. The chemical composition of the diets is presented in Table 1. Differences were detected (p ≤ 0.01) between the RG and PL diet for all measured chemical and nutritive value entities, except for DMD% (p = 0.15). Ryegrass had a 45% greater DM% content compared to PL. RG also had a 4% greater OM% content, 11% greater CP% content, 70% greater NDF% content, 9% greater ADF% content, 2% greater OMD% content, and 4% greater ME (MJ/kg DM) content, whilst PL had a 19% greater WSC% content.

3.2. Herbage Intake and Digestibility

No main effects of MUNBV or diet or diet by MUNBV interaction were detected for either herbage DMI or energy (ME, MJ/kg DM/d) intake (p > 0.05), with cows consuming, on average, 15.38 ± 1.13 kg DM/day and an average ME of 181 ± 8.6 MJ/d. Only an effect of diet (p < 0.05) was detected for N intake (g/d), with cows on the RG diet consuming on average 384.7 g N/d, which was a 15% greater N intake than cows consuming PL (333.0 g N/d). Only dietary effects were detected on truly digestible values (Table 2). Cows consuming the RG diet had a greater concentration (p < 0.05) of truly digestible values for DM (6%), OM (10%), N (12%), ADF (35%), and NDF (69%).

3.3. Milk Production

Table 3 presents the milk production data. A tendency for a MUNBV by diet effect (p = 0.06) was seen, with low MUNBV cows on RG having higher milk protein percentages. A tendency (p = 0.06) was only detected for MUNBV on milk protein yield, with low MUNBV cows tending to have a greater yield (kg/d) than high MUNBV cows. No main effect or interaction effects (p > 0.05) were detected for either milk (kg/d) or fat yield and fat percentage. Lactose yield (kg/d) was 9% greater for low MUNBV cows compared with high MUNBV cows (p < 0.01), with no effect detected for diet or the MUNBV by diet interaction (p > 0.05). A tendency (p = 0.09) was only detected for low MUNBV cows to have greater values for milk solids (fat + protein content, kg/d). A main effect of MUNBV (p < 0.01) was detected, with low MUNBV cows having a 7% increase in total solids (total production-water, kg/d) compared with high MUNBV cows. Only diet (p < 0.01) and MUNBV effect (p < 0.01) were detected for daily MUN (mg/dl), with cows on the PL diet having 30% less MUN (mg/dL) compared with the RG diet, whilst low MUNBV cows had a 15% decrease in MUN (mg/dL) compared with high MUNBV cows.

3.4. Rumen

The rumen characteristics are presented in Table 4. Glucogenic VFAs decreased 18% from day 1 to day 12 (p < 0.01) whilst non-glucogenic VFAs increased 5%. During the same period, the acetate to propionate ratio also increased by 26%. No main or interaction effects (p > 0.05) for NH3 or total VFA concentrations were detected.

3.5. Purine Excretion and Microbial Crude Protein Flow and Efficiency

Table 5 presents the purine derivatives, creatinine concentrations, and microbial N dynamics. A main effect of diet (p < 0.01) was detected, with cows fed ryegrass having a 43% greater concentration of creatinine (mmol/L), 75% greater concentration of hippuric acid (mmol/L), 55% greater concentration of allantoin (mmol/L), and 54% greater concentration of total purine derivatives (mmol/L) compared to cows fed PL.
An interaction between MUNBV and diet was detected (p = 0.04) for microbial crude protein flow (g/d), with low MUN cows on the RG diet having an 18% increase in flow compared with all other treatments. Neither diet, MUNBV, or the interaction (p > 0.05) affected the microbial efficiency (g of MCP/kg TDOM).

3.6. Excretory Pattern

Table 6 presents the excretion characteristics for urine and faeces. A MUNBV by diet interaction (p < 0.01) was detected for UUN g/L, resulting in high MUNBV cows consuming the RG diet having a 31% greater concentration than low MUNBV cows on RG. High and low MUNBV cows consuming RG had a 143% and 85% increased UUN g/L concentration compared with both high and low MUNBV cows consuming PL, respectively. Cows consuming PL had no difference in UUN g/L concentrations based on MUNBV (p > 0.05). A MUNBV by diet interaction (p < 0.01) was also detected for UN g/L, resulting in high MUNBV cows consuming the RG diet having a 29% greater concentration than their low MUNBV counterparts consuming the same diet. High MUNBV cows consuming RG and low MUNBV cows consuming RG were found to have a 123% and 73% greater concentration of UN g/L (p < 0.05), respectively, than cows on the PL diet. No difference was detected (p > 0.05) in the UN g/L concentration for cows on the PL diet regardless of MUNBV.
A MUNBV by diet interaction (p < 0.01) was detected for UUN, g/event, with high MUNBV cows consuming RG having a 38% greater concentration per urination event compared with low MUNBV cows on the same diet. A 165% and 92% increase (p < 0.05) in UUN, g/event was detected for high MUNBV and low MUNBV cows (respectively) on the RG diet compared with cows on the PL diet. No difference was detected between MUNBV groups consuming the PL diet (p > 0.05). Figure 1 presents the phenotypical relationship between MUN (mg/dL) and UUN g/event.
A relationship was detected, with g of UUN per event increasing 0.37 g per increase in MUN (mg/dL) across both diets. Diet shifted the intercept, with cows consuming RG having a higher intercept (4.26) compared to cows consuming PL (−0.27). An interaction between MUNBV and diet for UN, g/event (p < 0.01) was also detected. High MUNBV cows had a 39% greater UN excretion per event compared with low MUNBV cows on the RG diet. Both high and low MUNBV cows had a 145% and 77% greater UN excretion per event compared to cows consuming PL, respectively. No difference was detected (p > 0.05) between high and low MUNBV cows on the PL diet for UN excretions per urination event.
A main effect of both MUNBV and diet was detected (p < 0.01) on urine volume per event (L). Cows consuming RG had, on average, a 7% increase in urine volume/event compared with cows on a PL diet, whilst high MUNBV cows had a 7% increase in urine volume/event compared with low MUNBV cows. An interaction was found between MUNBV and diet (p < 0.01) for the number of urination events per day. High and low MUNBV cows consuming PL had no difference in urination events per day; however, on average, cows consuming PL had 3.45 and 6.05 more urination events per day compared to low and high MUNBV cows consuming RG, respectively. An interaction between MUNBV and diet (p < 0.05) was also detected for total daily urine volume (L/d). Low and high MUNBV cows consuming RG had a 22% and 34% reduction in daily urine volume compared to cows consuming PL, respectively. Only dietary effects (p < 0.01) were detected for both UUN, g/d and UN, g/d, with cows consuming RG having a 56% and 36% increase in UUN and UN g/d, respectively, compared to cows consuming PL. A diet by MUNBV interaction was detected for UUN as a percentage of total UN. High MUNBV cows fed RG had a 2% greater content of UUN as a percentage of total UN compared to low MUNBV cows on the same diet. Cows fed PL had a 12% and 10% lower content of UUN as a percentage of UN compared to high and low MUNBV cows fed RG, respectively.
A main effect of MUNBV and diet was detected for urine patch area, with high MUNBV cows having a 7% greater area per urine patch compared to low MUNBV cows. Cows on the RG diet also had a 7% greater area per urine patch compared to cows on the PL diet. An interaction effect was found between diet and MUNBV for UUN loading rate, with high MUNBV cows on RG diet having a 39% greater loading rate than low MUNBV cows on RG. High and low MUNBV cows on RG had a 164% and 90% greater loading rate (respectively) compared to cows on the PL diet, which were not considered different based on MUNBV.
A MUNBV by diet interaction was detected (p < 0.01) for the amount of N (g) per faecal event. Low MUNBV cows offered RG had a 41% increase in N content compared with high MUNBV cows on the same diet. Low MUNBV cows fed RG had a 22% increase whilst high MUNBV cows fed RG had a 14% decrease in N content per faecal event compared to cows of both high and low MUNBV offered PL. An interaction term was also detected for MUNBV by diet (p < 0.01) for OM content per faecal event. Low MUNBV cows offered RG had 45% more OM in faeces than high MUNBV cows on the same diet and were 23% and 11% higher than low and high MUNBV cows offered the PL diet, respectively. An interaction term was also detected for MUNBV by diet (p < 0.01) for ADF content per faecal event, with high MUNBV cows offered the PL diet having 8% more ADF in the faeces than low MUNBV cows offered the same diet and 75% and 23% more than high and low MUNBV cows on the RG diet, respectively. A MUNBV by diet interaction term was also detected (p < 0.01) for NDF content per faecal event, with high MUNBV cows offered PL having 56% more NDF than high MUNBV cows offered RG. High MUNBV cows on RG had an 8% greater NDF content per faecal event than low MUNBV cows irrespective of diet. High MUNBV cows on RG had a 19% lower content of NDF per faecal event compared to low MUNBV cows irrespective of diet. An interaction (p < 0.01) between MUNBV and diet was detected for the DM content (g) per faecal event. Low MUNBV cows offered the RG diet had 45% more DM/event than high MUNBV cows offered the same diet and 31% and 18% more DM/event than low and high MUNBV cows offered the PL diet, respectively. A diet by MUNBV interaction (p = 0.03) was also detected for the number of faecal events per day. Low MUNBV cows offered RG had 29% fewer faecal events than high MUNBV cows offered both diets and low MUNBV cows on the PL diet, all of which did not differ (p > 0.05). Only a diet effect was detected (p < 0.01) for total daily faecal N content (g) and total daily faecal DM content (kg), with cows offered the PL diet having a 15% greater N content and 14% greater total DM content per day compared with cows on the RG diet.

3.7. Water Balance

Table 7 presents the average daily water balance. Main effects of MUNBV (p = 0.03) and diet (p < 0.01) were detected on daily water consumed from the trough. Low MUNBV cows drank 44% more water than high MUNBV cows, whilst cows on the RG diet consumed 6 times more water than cows on the PL diet. A dietary effect was also detected for daily water consumed in feed (p < 0.01), with cows on the PL diet consuming 58% more water in the forage than cows on the RG diet. No effect of MUNBV or interaction between MUNBV and diet were detected for total water consumption per day. However, a tendency (p = 0.06) was detected for diet, with cows consuming PL having larger volumes of water intake than cows consuming RG. An interaction between MUNBV and diet (p < 0.05) resulted in less total daily water excretion via urine from high MUNBV cows consuming RG compared to low MUNBV cows consuming RG, whilst cows consuming PL had the greatest volume of water excreted. An effect of diet (p < 0.01) was detected for water excreted in faeces, with cows on the RG diet having 16% greater volume of water excreted than cows on PL. A tendency (p = 0.09) was detected for MUNBV to have an effect, with high MUNBV cows having larger volumes of water excreted in their faeces. No effect of diet, MUNBV, or their interaction term was detected for the volume of water excreted in milk (p > 0.05). Only a dietary effect was detected (p < 0.01) for total water excreted, with cows on PL having a 10% greater volume of total water excretion than cows on RG diet. No diet, MUNBV, or interaction effects were detected for overall water balance (p > 0.05). A MUNBV (p = 0.02) and diet (p < 0.01) effect was detected for the amount of water absorbed (%) as a proportion of intake. Cows with low MUNBV had 5% more water absorbed as a proportion of intake compared with high MUNBV cows whilst cows on PL had 11% more water absorbed as a proportion of intake than cows on RG.

3.8. Nitrogen Balance and Partitioning

Table 8 presents the N balance and N partitioning. An effect of diet was detected (p < 0.05) for differential N partitioning for milk, urine, and faeces. Cows consuming PL partitioned more N into milk (13%) and faeces (30%) but less into the urine (14%) than cows consuming RG. A diet effect was detected for the total amount of N excreted, with cows consuming PL excreting 9% less N compared to cows consuming RG. No main effects of MUN or diet or an interaction term (p > 0.05) were detected for the N balance.

4. Discussion

We investigated whether cows divergent for MUNBV had a difference in their N balance and if there were any dietary effects or interaction between diet and MUNBV, with the hypothesis that cows divergent for MUNBV would have differences in their N balance and an additive effect of diet and MUNBV on reducing UUN excretions. Our results indicate that there was no difference in the total amount of N excreted based on MUNBV, which is consistent with the recent work of Müller et al. (2021) [49]. There was an effect of diet on total daily N excreted but no interaction with MUNBV. However, our results do indicate that the excretion patterns and characteristics of excreta differ according to MUNBV and diet, which can result in different environmental impacts.
Only dietary effects were detected for the daily N balance, with cows fed the PL diet having a lower amount of total N excreted. This difference can largely be accounted for as the cows on the PL diet had a lower N intake compared with RG, thus resulting in less N excreted. The low content of CP in the RG diet may account for the low N excretion rates from RG compared to the literature [50]. Differential N partitioning from the PL diet was observed, with cows excreting more N through faeces than urine compared with RG cows, which has been reported previously [51,52]. The increased N partitioning to faeces by PL is likely a result of the different plant secondary compounds (aucubin, actoside, catalpol) found in PL compared to RG, which has been reported to influence N partitioning in other studies [51]. A study investigating the use of tannin-containing hay showed an increased partitioning of N into faeces rather than urine [53] similar to what was observed in this study. Increased N partitioning to faeces represents an environmental advantage for PL diets, as N in faeces is less at risk of leaching or volatilising to the environment [54,55]. Different partitioning of N was also observed in the urine, with RG cows having greater concentrations of purine derivatives. The differences in purine derivatives can largely be explained by the different digestibility, energy, and CP content of the diets [56]. The interaction between MUNBV and RG for UUN as a percent UN also indicates a difference in N partitioning within the urine, with high MUNBV cows having a greater level of UUN as a percent of UN when consuming RG. This relationship may warrant further investigation.
A difference in lactose yield was detected, with low MUNBV cows having a greater lactose yield than high MUNBV animals. A similar relationship was reported by Marshall et al. (2020) [32]. Whilst the New Zealand dairy industry does not reward lactose yield in milk payments, there has been growing interest internationally in the refinement of lactose from milk products for use as a component in pharmaceutical products and lactose-fortified food products [57,58]. This could theoretically offer an economic incentive for the use of low MUNBV cows in markets where milk producers pay for lactose content.
Although not different for total daily N excreted, low MUNBV cows consistently had lower concentrations of both UUN and UN when consuming ryegrass, resulting in lower values on a per urination event basis compared with high MUNBV cows. This paired with lower urine volumes per event results in less N being loaded per urine patch for low MUNBV cows consuming ryegrass. The urine patch is considered the “engine room” of N cycling in grazed pastures as it ultimately determines the amount of N deposited to pasture, which can then either be used by plants or lost to the environment [7]. It is possible to estimate the amount of potential nitrate leached based on the urine N loading rates at the urine patch using the empirical model of Di and Cameron (2007) [59] (5):
y = 16.7 + 0.173 x + 0.000071 x 2
where y = total nitrate nitrogen (NO3N) leached (kg N ha 1 ) and x = UUN load onto pasture (kg/ ha 1 ). A 1% difference exists between the estimated leaching rate of high and low MUNBV cows consuming PL per urine patch. Low MUNBV cows consuming RG are estimated to have a 23% reduction in nitrate leaching compared to high MUNBV cows consuming RG. Cows consuming PL regardless of MUNBV are estimated to have a 32% and 47% reduction in nitrate leaching compared to low and high MUNBV cows consuming RG, respectively. These reductions in estimated nitrate leaching from the urine patch indicate the ability for both animal genetics with regards to MUNBV when grazing RG diets and the use of PL to reduce the amount of N lost from leaching events. Reductions in N leaching are considered a key objective for reducing eutrophication [60] and therefore meeting legislative concerns nationally and internationally regarding freshwater quality and for reducing any possible negative human health impacts from high N concentrations in drinking water. Not only will the environmental impact on waterways be reduced through the use of MUNBV on RG diets and the use of PL, but the ability to reduce UUN content in the urine will also help reduce GHG emissions. A lower content of UUN per urine patch would be expected to yield comparatively lower amounts of N2O emissions [61] and thus reduce GHG emissions from low MUNBV cows grazing RG or cows grazing PL.
The tendency for low MUNBV cows and the effect of the PL diet to have more urination events means that N deposited in urine is likely to be spread over a larger spatial area over a longer time period, which could reduce the N loading in any one area. However, on a per hectare basis, this may result in greater overlapping of urine patches. The effectiveness of low MUNBV cows and the PL diet to reduce N leaching should be further assessed on a per hectare basis. The use of PL should also be assessed from an animal welfare perspective in regard to potential collateral long-term effects of diuresis.
The lack of effect from MUNBV on total daily N excretion from urine despite the consistently lower N concentration on RG diets from low MUNBV cows may be explained by the greater daily urine volume from low MUNBV cows on RG. A greater volume of excretion with a lower N concentration from low MUNBV cows on RG can equate to a similar level of N excretion of a smaller volume but a higher N concentration from high MUNBV cows on RG.
A similar relationship as to that seen with urination characteristics and total N balance was observed in faecal characteristics. Although there was no difference in total daily N excreted in faeces by MUNBV, the behavioural characteristics of the faecal events differed. A greater partitioning of N to faeces from low MUNBV cows consuming RG was observed on a per-event basis, which has associated positive environmental benefits as mentioned earlier when discussing the PL effect of N in total daily faeces. The differences caused by the diet faction of the diet by MUNBV interaction can largely be attributed to the different digestibility observed between RG and PL diets; however, no differences were detected in digestibility variables by MUNBV. As these differences are ‘seen’ per event, they might therefore suggest different dynamics in rumen digesta outflow based on MUNBV. Thus, this results in a different temporal pattern of excretion of the non-digestible diet components.
The difference in MCP flow between the MUNBV groups supports the idea that these cows could have had different rumen digesta outflow dynamics. Low MUNBV cows fed RG had a greater MCP flow, which would result in a greater content of amino acids of high biological value supplied to the small intestine and which could have increased animal productivity [44,62] as potentially represented by the tendency for greater milk protein yield from low MUNBV cows on RG. Similar results were detected by Souza et al. (2021) [63], who found that cows considered high for MUN had a suppressed rate of urea transport to the gastrointestinal tract compared to low MUN cows. Souza et al. (2021) [63] went on to hypothesize that these high MUN cows may therefore be more susceptible to rumen degradable protein deficiencies than low MUN cows, which may have also occurred in this study.
The increased MCP flow from the low MUNBV cows paired with greater voluntary water intake is indicative of a difference in the rumen osmotic concentration [64]. Osmolality may be increased due to increased rates of fermentation [65] in low MUNBV cows. The latter may either be as a result of different rumen microbial communities, greater oral processing of ingesta, or both. Cows divergent for genetic parameters have been reported to have different oral processing behaviour of ingesta [66]. If low MUNBV cows have increased oral processing, it could result in smaller particle sizes entering the rumen [67], which would increase rumen fermentation rates and digesta outflow flow, potentially accounting for some of the differences observed based on MUNBV. No differences in total VFA concentration were detected; however, no measurements of total VFA pools or absorption were made, thus further research is required in this area.
If cows did have a difference in rumen digesta outflow, it would also help explain differences in urination behaviour based on MUNBV. Greater rumen digesta outflow rates, especially the liquid phase, represent a greater content of available water in the body’s water pool that can be used in respiration, evaporation, milk, faeces, or urine [68]. Cows are stimulated to urinate when the bladder capacity reaches 50% [68]. If there was a greater continual outflow of liquid from the rumen entering the body’s water pool and eventually the bladder, it would theoretically reach this threshold at a higher frequency, resulting in more urination events per day as per what was seen with low MUNBV cows. Additionally, low MUNBV cows may have a lower urea N concentration threshold, in which case the dilution function of the kidney would be exacerbated and ultimately increasing the liquid flow to the bladder.
A lack of diet effect from PL was seen throughout this experiment for cows differing in MUNBV, e.g., no difference was observed in UUN concentrations between high and low MUNBV despite the fact this relationship has been established previously on an RG diet [32]. This may in part be due to the interaction of the MUNBV model and differing diets. A genetic model was used to predict MUNBV values, and the parameters used are based on cows grazing predominantly RG/white clover pastures in New Zealand. There is the potential that these parameters may be different if these cows were grazing alternative diets, such as PL, compared to the typical RG/white clover diet used in New Zealand.
Figure 1 presents the phenotypical relationship between MUN (mg/dL) and UUN g/event. A clear effect of MUN (mg/dL) was observed, where an increasing MUN (mg/dL) content increases UNN g/event. This represents the physiological relationship that is occurring within these cows, which was unable to be captured based on a BV for the PL diet, indicating that the model’s parameters may not be appropriate for the PL diet as the outcome does not match the data generating process [69]. New Zealand dairy pastoral systems are largely RG based [70], which would allow for the current MUNBV models to be used with sufficient accuracy to determine N losses in the typical herd. However, the model parameters should be re-assessed using the parameters of cows grazing multiple different forage species—a clear opportunity to add functionality to the model—before any legislative or management decisions are made based on cow MUNBV values in production systems that are not based on RG-dominated swards.

5. Conclusions

No effect of MUNBV was detected for total daily urinary N excretions, indicating that cows divergent for MUNBV will be depositing similar amounts of N to pasture via urine on a total daily basis. However, due to the different urination patterns from low MUNBV cows consuming RG, the environmental impact would be expected to differ. Due to the lower concentration of N per event and the subsequent lower urine N loading rate, it is estimated that low MUNBV cows will have less N leaching from urine patches and will excrete more N as faeces when consuming RG pastures compared to high MUNBV cows, thus potentially reducing the environmental impact. Plantain diets also offer a potential opportunity to reduce the environmental impact of pastoral dairy production systems. The PL diet not only resulted in more partitioning of N into faeces, it also reduced the content of N excreted and therefore potentially leached at the urine patch as well as reducing the total daily N excreted compared to the RG diet for both high and low MUNBV cows. Although a limited effect of breeding value for MUN was detected when cows were fed PL, a phenotypical effect was detected with low MUN (mg/dL), resulting in less UUN g/L and therefore reduced environmental impact. The MUNBV model’s parameters must be re-assessed for different forage species, e.g., PL. Once the diet is taken into account in generating MUNBV values, it could be argued that the PL effect may be additive to the positive effect of low MUNBV cows.

Author Contributions

C.J.M.: Conceptualization, methodology, validation, formal analysis, resources, data curation, writing- original draft, writing—review & editing, visualization, supervision, project administration, funding acquisition. M.R.B.; Methodology, validation, formal analysis, investigation, resources, data curation, writing—review & editing. K.G.; Methodology, validation, formal analysis, investigation, resources, data curation, writing—review & editing. G.K.B.: Conceptualization, methodology, writing—review & editing. O.A.-M.: Writing—review & editing. P.G.: Conceptualization, methodology, resources, supervision, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

Thank you to the Lincoln University Open Access Fund for helping fund the publication costs of this paper.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Lincoln University Animal Ethics Committee (AEC 2019-25A).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The first authors thesis is supported by CRV, Agricom, The Lincoln University Pastoral Livestock production lab, Fonterra, and the Lincoln University Doctorial scholarship.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pacheco, D.; Waghorn, G. Dietary nitrogen—Definitions, digestion, excretion and consequences of excess for grazing ruminants. Proc. N. Z. Grassl. Assoc. 2008, 70, 107–116. [Google Scholar] [CrossRef]
  2. Castillo, A.R.; Kebreab, E.; Beever, D.E.; Barbi, J.H.; Sutton, J.D.; Kirby, H.C.; France, J. The effect of protein supplementation on nitrogen utilization in lactating dairy cows fed grass silage diets. J. Anim. Sci. 2001, 79, 247–253. [Google Scholar] [CrossRef]
  3. Kebreab, E.; France, J.; Beever, D.E.; Castillo, A.R. Nitrogen pollution by dairy cows and its mitigation by dietary manipulation. Nutr. Cycl. Agroecosyst. 2001, 60, 275–285. [Google Scholar] [CrossRef]
  4. Gregorini, P.; Beukes, P.C.; Dalley, D.; Romera, A.J. Screening for diets that reduce urinary nitrogen excretion and methane emissions while maintaining or increasing production by dairy cows. Sci. Total Environ. 2016, 551–552, 32–41. [Google Scholar] [CrossRef] [Green Version]
  5. Oudshoorn, F.W.; Kristensen, T.; Nadimi, E. Dairy cow defecation and urination frequency and spatial distribution in relation to time-limited grazing. Livest. Sci. 2008, 113, 62–73. [Google Scholar] [CrossRef]
  6. Clark, C.; Waghorn, G.; Gregorini, P.; Woodward, S.; Clark, D. Diurnal pattern of urinary and faecal nitrogen excretion by dairy cows fed ryegrass pasture twice daily indoors. Adv. Anim. Biosci. 2010, 2, 269. [Google Scholar]
  7. Selbie, D.R.; Buckthought, L.E.; Shepherd, M.A. The challenge of the urine patch for managing nitrogen in grazed pasture systems. Adv. Agron. 2015, 129, 229–292. [Google Scholar]
  8. Eggleston, H.S.; Buendia, L.; Miwa, K.; Ngara, T.; Tanabe, K. IPCC Guidelines for National Greenhouse Gase Inventories; IPCC: Kanagawa, Japan, 2006. [Google Scholar]
  9. Del Grosso, S.J.; Wirth, T.; Ogle, S.; Parton, W.J. Estimating agricultural nitrous oxide emissions. Eos 2008, 89, 529. [Google Scholar] [CrossRef]
  10. Schallenberg, M.; de Winton, M.D.; Verburg, P.; Kelly, D.J.; Hamill, K.D.; Hamilton, D.P. Ecosystem services of lakes. In Ecosystem Services in New Zealand—Conditions and Trends; Dymond, J.R., Ed.; Manaaki Whenua Press: Lincoln, New Zealand, 2013. [Google Scholar]
  11. Gregorini, P.; Maxwell, T. Grazing in future multiscapes—From thoughtscapes to ethical and sustainable foodscapes. Journal 2020, 24, 23–25. [Google Scholar]
  12. United Nations. Sustainable Development Goals. 2015. Available online: https://sustainabledevelopment.un.org/sdgs (accessed on 17 March 2021).
  13. Schullehner, J.; Hansen, B.; Thygesen, M.; Pedersen, C.B.; Sigsgaard, T. Nitrate in drinking water and colorectal cancer risk: A nationwide population-based cohort study. Int. J. Cancer 2018, 143, 73–79. [Google Scholar] [CrossRef]
  14. Johnson, S.F. Methemoglobinemia: Infants at risk. Curr. Probl. Pediatr. Adolesc. Health Care 2019, 49, 57–67. [Google Scholar] [CrossRef]
  15. Brender, J.D.; Olive, J.M.; Felkner, M.; Suarez, L.; Marckwardt, W.; Hendricks, K.A. Dietary nitrites and nitrates, nitrosatable drugs, and neural tube defects. Epidemiology 2004, 15, 330–336. [Google Scholar] [CrossRef]
  16. Thomson, B.C.; Ward, K.; Smith, N.; Gibbs, J.; Muir, P.D. Effect of feeding time on urinary and faecal nitrogen excretion patterns in sheep. N. Z. J. Agric. Res. 2021, 64, 314–319. [Google Scholar] [CrossRef]
  17. Beukes, P.C.; Gregorini, P.; Romera, A.J.; Woodward, S.L.; Khaembah, E.; Chapman, D.F.; Nobilly, F.; Bryant, R.H.; Edwards, G.; Clark, D.A. The potential of diverse pastures to reduce nitrogen leaching on New Zealand dairy farms. Anim. Prod. Sci. 2014, 54, 1971–1979. [Google Scholar] [CrossRef] [Green Version]
  18. Garrett, K.; Beck, M.R.; Marshall, C.J.; Fleming, A.; Logan, C.M.; Maxwell, T.M.R.; Greer, A.W.; Gregorini, P. Functional diversity vs. monotony: The effect of a multiforage diet as opposed to a single forage diet on animal intake, performance, welfare, and urinary nitrogen excretion. J. Anim. Sci. 2021, 99, skab058. [Google Scholar] [CrossRef]
  19. Carmona-Flores, L.; Bionaz, M.; Downing, T.; Sahin, M.; Cheng, L.; Ates, S. Milk production, N partitioning, and Methane emissions in dairy cows grazing mixed or spatially separated simple and diverse pastures. Animals 2020, 10, 1301. [Google Scholar] [CrossRef] [PubMed]
  20. Box, L.A.; Edwards, G.; Bryant, R.H. Milk production and urinary nitrogen excretion of dairy cows grazing plantain in early and late lactation. N. Z. J. Agric. Res. 2017, 60, 470–482. [Google Scholar] [CrossRef]
  21. Bryant, R.H.; Snow, V.O.; Shorten, P.R.; Welten, B.G. Can alternative forages substantially reduce N leaching? Findings from a review and associated modelling. N. Z. J. Agric. Res. 2019, 63, 3–28. [Google Scholar] [CrossRef]
  22. Stewart, A. Plantain (Plantago lanceolata)–A potential pasture species. Proc. N. Z. Grassl. Assoc. 1996, 86, 77–86. [Google Scholar] [CrossRef]
  23. Navarrete, S.; Kemp, P.D.; Pain, S.J.; Back, P.J. Bioactive compounds, aucubin and acteoside, in plantain (Plantago lanceolata L.) and their effect on in vitro rumen fermentation. Anim. Feed Sci. Technol. 2016, 222, 158–167. [Google Scholar] [CrossRef]
  24. O’Connell, C.; Judson, H.; Barrell, G.K. Sustained diuretic effect of plantain when ingested by sheep. N. Z. Soc. Anim. Prod. 2016, 76, 14–17. [Google Scholar]
  25. Butler, W.; Calaman, J.J.; Beam, S.W. Plasma and milk urea nitrogen in relation to pregnancy rate in lactating dairy cattle. J. Anim. Sci. 1996, 74, 858–865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kauffman, A.; St-Pierre, N. The relationship of milk urea nitrogen to urine nitrogen excretion in holstein and jersey cows. J. Dairy Sci. 2001, 84, 2284–2294. [Google Scholar] [CrossRef]
  27. Kohn, R.; Kalscheur, K.; Russek-Cohen, E. Evaluation of models to estimate urinary nitrogen and expected milk urea Nitrogen. J. Dairy Sci. 2002, 85, 227–233. [Google Scholar] [CrossRef]
  28. Kohn, R.A.; Dinneen, M.M.; Russek-Cohen, E. Using blood urea nitrogen to predict nitrogen excretion and efficiency of nitrogen utilization in cattle, sheep, goats, horses, pigs, and rats. J. Anim. Sci. 2005, 83, 879–889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Stoop, W.; Bovenhuis, H.; van Arendonk, J. Genetic parameters for milk urea Nitrogen in relation to milk production traits. J. Dairy Sci. 2007, 90, 1981–1986. [Google Scholar] [CrossRef] [PubMed]
  30. Lopez-Villalobos, N.; Correa-Luna, M.; Burke, J.L.; Sneddon, N.W.; Schutz, M.; Donaghy, D.; Kemp, P. Genetic parameters for milk urea concentration and milk traits in New Zealand grazing dairy cattle. N. Z. J. Anim. Sci. Prod. 2018, 78, 56–61. [Google Scholar]
  31. Beatson, P.R.; Meier, S.; Cullen, N.G.; Eding, H. Genetic variation in milk urea nitrogen concentration of dairy cattle and its implications for reducing urinary nitrogen excretion. Animal 2019, 13, 2164–2171. [Google Scholar] [CrossRef] [Green Version]
  32. Marshall, C.; Beck, M.; Garrett, K.; Barrell, G.; Al-Marashdeh, O.; Gregorini, P. Grazing dairy cows with low milk urea nitrogen breeding values excrete less urinary urea nitrogen. Sci. Total Environ. 2020, 739, 139994. [Google Scholar] [CrossRef] [PubMed]
  33. Primary Industries Standing Committee. Nutrient Requirements of Domesticated Ruminants; CSIRO Publishing: Clayton, Australia, 2007. [Google Scholar] [CrossRef]
  34. Iowerth, D.; Jones, H.; Hayward, M.V. The effect of pepsin pretreatment of herbage on the prediction of dry matter digestibility from solubility in fungal cellulase solutions. J. Sci. Food Agric. 1975, 26, 711–718. [Google Scholar] [CrossRef]
  35. MAFF—Ministry of Agriculture, Fisheries and Food. Carbohydrates, soluble, in herbage. In The Analysis of Agricultural Materials: A Manual of the Analytical Methods Used by the Agricultural Development and Advisory Service; MAFF: London, UK, 1986. [Google Scholar]
  36. AOAC. Official Methods of Analysis, 15th ed.; Association of Official Analytical Chemist: Washington, DC, USA, 1990. [Google Scholar]
  37. Van Soest, P.; Robertson, J.; Lewis, B. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 1991, 74, 3583–3597. [Google Scholar] [CrossRef]
  38. Haynes, R.J.; Williams, P.H. Nutrient cycling and soil fertility in grazed pature ecosystem. Adv. Agron. 1993, 49, 119–199. [Google Scholar]
  39. Romera, A.J.; Levy, G.; Beukes, P.; Clark, D.A.; Glassey, C.B. A urine patch framework to simulate nitrogen leaching on New Zealand dairy farms. Nutr. Cycl. Agroecosyst. 2012, 92, 329–346. [Google Scholar] [CrossRef]
  40. George, S.; Dipu, M.; Mehra, U.; Singh, P.; Verma, A.; Ramgaokar, J. Improved HPLC method for the simultaneous determination of allantoin, uric acid and creatinine in cattle urine. J. Chromatogr. B 2006, 832, 134–137. [Google Scholar] [CrossRef] [PubMed]
  41. International Atomic Energy Agency. Estimation of Rumen Microbial Protein Production From Purine Derivatives In Urine—A Laboratory Manual for FAO/IAEA Co-Ordinated Research Program; IAEA: Vienna, Austria, 1997. [Google Scholar]
  42. Crawford, G.I.; Macdonald, J.C.; Watson, A.K.; Erickson, G.; Klopfenstein, T.J. Diurnal and dietary impacts on estimating microbial protein flow from urinary purine derivative excretion in beef cattle. Transl. Anim. Sci. 2020, 4, txaa140. [Google Scholar] [CrossRef]
  43. Martín-Orúe, S.M.; Balcells, J.; Guada, J.; Fondevila, M. Microbial nitrogen production in growing heifers: Direct measurement of duodenal flow of purine bases versus urinary excretion of purine derivatives as estimation procedures. Anim. Feed. Sci. Technol. 2000, 88, 171–188. [Google Scholar] [CrossRef]
  44. Clark, J.; Klusmeyer, T.; Cameron, M. Microbial protein synthesis and flows of Nitrogen fractions to the duodenum of dairy cows. J. Dairy Sci. 1992, 75, 2304–2323. [Google Scholar] [CrossRef]
  45. Talke, H.; Schubert, G.E. Enzymatische harnstoffbestimmung in blut und serum im optischen test nachwarburg. Klin Wochschr 1965, 43, 174–175. [Google Scholar] [CrossRef]
  46. Playne, M.J. Determination of ethanol, volatile fatty acids, lactic and succinic acids in fermentation liquids by gas chromatography. J. Sci. Food Agric. 1985, 36, 638–644. [Google Scholar] [CrossRef]
  47. Neeley, W.E.; Phillipson, J. Automated enzymatic method for determining ammonia in plasma, with 14-day reagent stability. Clin. Chem. 1988, 34, 1868–1869. [Google Scholar] [CrossRef]
  48. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  49. Müller, C.B.M.; Görs, S.; Derno, M.; Tuchscherer, A.; Wimmers, K.; Zeyner, A.; Kuhla, B. Differences between Holstein dairy cows in renal clearance rate of urea affect milk urea concentration and the relationship between milk urea and urinary Nitrogen excretion. Sci. Total Environ. 2020, 755, 143198. [Google Scholar] [CrossRef]
  50. Pacheco, D.; Lowe, K.; Burke, J.L.; Cosgrove, G.P. Urinary nitrogen excretion from cows at different stage of lactation grazing different ryegrass cultivars during spring or autumn. N. Z. Soc. Anim. Sci. 2009, 69, 196–200. [Google Scholar]
  51. Minnée, E.; Leach, C.; Dalley, D. Substituting a pasture-based diet with plantain (Plantago lanceolata) reduces nitrogen excreted in urine from dairy cows in late lactation. Livest. Sci. 2020, 239, 104093. [Google Scholar] [CrossRef]
  52. Beck, M.R.; Garrett, K.; Thompson, B.R.; Stevens, D.R.; Barrell, G.K.; Gregorini, P. Plantain (Plantago lanceolata) reduces the environmental impact of farmed red deer (Cervus elaphus). Transl. Anim. Sci. 2020, 4, txaa160. [Google Scholar] [CrossRef]
  53. Stewart, E.K.; Beauchemin, K.; Dai, X.; Macadam, J.W.; Christensen, R.G.; Villalba, J.J. Effect of tannin-containing hays on enteric methane emissions and nitrogen partitioning in beef cattle. J. Anim. Sci. 2019, 97, 3286–3299. [Google Scholar] [CrossRef] [PubMed]
  54. Waghorn, G. Beneficial and detrimental effects of dietary condensed tannins for sustainable sheep and goat production—Progress and challenges. Anim. Feed. Sci. Technol. 2008, 147, 116–139. [Google Scholar] [CrossRef]
  55. Powell, J.; Broderick, G.; Grabber, J.; Hymes-Fecht, U. Technical note: Effects of forage protein-binding polyphenols on chemistry of dairy excreta. J. Dairy Sci. 2009, 92, 1765–1769. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, X.B. Excretion of Purine Derivatives by Sheep and Cattle and Its Use for the Estimation of Absorbed Microbial Protein. Ph.D. Thesis, University of Aberdeen, Aberdeen, UK, 1989. [Google Scholar]
  57. Løvendahl, P.; Weisbjerg, M.R. Lactose in milk—How can lactose concentration data be beneficial in management and breeding? ICAR Tech. Ser. 2018, 22, 99–104. [Google Scholar]
  58. Wong, S.Y.; Hartel, R.W. Crystallization in lactose refining—A review. J. Food Sci. 2014, 79, R257–R272. [Google Scholar] [CrossRef] [PubMed]
  59. Di, H.J.; Cameron, K. Nitrate leaching losses and pasture yields as affected by different rates of animal urine nitrogen returns and application of a nitrification inhibitor—A lysimeter study. Nutr. Cycl. Agroecosyst. 2007, 79, 281–290. [Google Scholar] [CrossRef]
  60. Smolders, A.J.P.; Lucassen, E.C.H.E.T.; Bobbink, R.; Roelofs, J.G.M.; Lamers, L.P.M. How nitrate leaching from agricultural lands provokes phosphate eutrophication in groundwater fed wetlands: The sulphur bridge. Biogeochemistry 2009, 98, 1–7. [Google Scholar] [CrossRef] [Green Version]
  61. Chadwick, D.; Cardenas, L.; Dhanoa, M.; Donovan, N.; Misselbrook, T.; Williams, J.; Thorman, R.; McGeough, K.; Watson, C.; Bell, M.; et al. The contribution of cattle urine and dung to nitrous oxide emissions: Quantification of country specific emission factors and implications for national inventories. Sci. Total Environ. 2018, 635, 607–617. [Google Scholar] [CrossRef]
  62. Hackmann, T.J.; Firkins, J.L. Maximizing efficiency of rumen microbial protein production. Front. Microbiol. 2015, 6, 465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Souza, V.; Aguilar, M.; Van Amburgh, M.; Nayananjalie, W.; Hanigan, M. Milk urea nitrogen variation explained by differences in urea transport into the gastrointestinal tract in lactating dairy cows. J. Dairy Sci. 2021, 104, 6715–6726. [Google Scholar] [CrossRef]
  64. Rogers, J.; Marks, B.; Davis, C.; Clark, J. Alteration of rumen fermentation in steers by increasing rumen fluid dilution rate with mineral salts. J. Dairy Sci. 1979, 62, 1599–1605. [Google Scholar] [CrossRef]
  65. Castillo-González, A.; Burrola-Barraza, M.; Domínguez-Viveros, J.; Chávez-Martínez, A. Rumen microorganisms and fermentation. Arch. Med. Veter. 2014, 46, 349–361. [Google Scholar] [CrossRef] [Green Version]
  66. Gregorini, P.; Waghorn, G.; Kuhn-Sherlock, B.; Romera, A.; Macdonald, K. Short communication: Grazing pattern of dairy cows that were selected for divergent residual feed intake as calves. J. Dairy Sci. 2015, 98, 6486–6491. [Google Scholar] [CrossRef]
  67. Beauchemin, K. Invited review: Current perspectives on eating and rumination activity in dairy cows. J. Dairy Sci. 2018, 101, 4762–4784. [Google Scholar] [CrossRef] [Green Version]
  68. Gregorini, P.; Provenza, F.D.; Villalba, J.J.; Beukes, P.; Forbes, M.J. Diurnal patterns of urination and drinking by grazing ruminants: A development in a mechanistic model of a grazing ruminant, MINDY. J. Agric. Sci. 2018, 156, 71–81. [Google Scholar] [CrossRef] [Green Version]
  69. Henley, S.S.; Golden, R.M.; Kashner, T.M. Statistical modeling methods: Challenges and strategies. Biostat. Epidemiol. 2019, 4, 105–139. [Google Scholar] [CrossRef] [Green Version]
  70. Ledgard, S.F.; Wei, S.; Wang, X.; Falconer, S.; Zhang, N.; Zhang, X.; Ma, L. Nitrogen and carbon footprints of dairy farm systems in China and New Zealand, as influenced by productivity, feed sources and mitigations. Agric. Water Manag. 2018, 213, 155–163. [Google Scholar] [CrossRef]
Figure 1. Average daily urinary urea N (UUN) excretion g per event (y-axis) as a function of the average daily milk urea N (MUN, mg/dL) (x-axis) from lactating dairy cows offered either a ryegrass or plantain diet whilst housed in metabolism creates for 72 h. For both diets, a one unit increase in MUN (mg/dL) resulted in a 0.37 ± 0.09 increase in the UUN content (g) per urination event (p < 0.01). Diet shifted the intercept of the regression line to 4.26 ± 1.42 for the ryegrass diet and −0.27 ± 1.42 for the plantain diet (p < 0.01). The regression model R2 = 0.63, with the shaded area representing the 95% confidence interval.
Figure 1. Average daily urinary urea N (UUN) excretion g per event (y-axis) as a function of the average daily milk urea N (MUN, mg/dL) (x-axis) from lactating dairy cows offered either a ryegrass or plantain diet whilst housed in metabolism creates for 72 h. For both diets, a one unit increase in MUN (mg/dL) resulted in a 0.37 ± 0.09 increase in the UUN content (g) per urination event (p < 0.01). Diet shifted the intercept of the regression line to 4.26 ± 1.42 for the ryegrass diet and −0.27 ± 1.42 for the plantain diet (p < 0.01). The regression model R2 = 0.63, with the shaded area representing the 95% confidence interval.
Animals 11 02464 g001
Table 1. Forage chemical composition of plantain and ryegrass diets offered to cows of differing MUNBVs (milk urea N breeding value) housed in metabolism crates.
Table 1. Forage chemical composition of plantain and ryegrass diets offered to cows of differing MUNBVs (milk urea N breeding value) housed in metabolism crates.
Diet p-Value
Item 1PlantainRyegrassRSD 2Diet
DM, %12.818.50.01<0.01
OM, %87.790.90.86<0.01
CP, %13.414.92.000.01
WSC, %19.116.04.390.02
NDF, %27.346.52.45<0.01
ADF, %24.126.21.32<0.01
DMD, %74.075.12.520.15
OMD, %77.679.53.130.04
ME, MJ/kg DM11.111.50.43<0.01
1 DM, Dry Matter; OM, Organic Matter; CP, Crude Protein; WSC, Water-Soluble Carbohydrates; NDF, Neutral Detergent Fibre; ADF, Acid Detergent Fibre; DMD, Dry Matter Digestibility; OMD, Organic Matter Digestibility; ME, Metabolizable Energy. 2 RSD = Residual standard deviation.
Table 2. Effect of diet (ryegrass and plantain) and MUNBV (milk urea N breeding value) classified as either ‘high’ or ‘low’ on true digestible values for dairy cows housed in metabolism crates.
Table 2. Effect of diet (ryegrass and plantain) and MUNBV (milk urea N breeding value) classified as either ‘high’ or ‘low’ on true digestible values for dairy cows housed in metabolism crates.
High MUNBVLow MUNBV p-Value
Item 1PlantainRyegrassPlantainRyegrassRSD 2MUNBVDietMUNBV × Diet
Dry
Matter
7074.569.774.14.450.8<0.010.49
Organic Matter7076.869.576.43.990.72<0.010.5
Nitrogen647264.272.14.790.92<0.010.39
ADF54.974.353.6720.10.32<0.010.33
NDF4678.444.774.80.130.36<0.010.22
1 ADF = Acid Detergent Fibre; NDF = Neutral Detergent Fibre. 2 RSD = Residual standard deviation.
Table 3. Effect of diet (ryegrass and plantain) and MUNBV (milk urea N breeding value) classified as either ‘high’ or ‘low’ on daily milk production characteristics for dairy cows housed in metabolism crates.
Table 3. Effect of diet (ryegrass and plantain) and MUNBV (milk urea N breeding value) classified as either ‘high’ or ‘low’ on daily milk production characteristics for dairy cows housed in metabolism crates.
High MUNBVLow MUNBV p-Value
Item 1PlantainRyegrassPlantainRyegrassRSD 2MUNBVDietMUNBV × Diet
Milk yield, kg18.618.319.3191.730.150.580.74
Fat, %4.544.784.514.750.010.860.190.3
Fat yield, kg0.880.840.840.860.010.290.160.27
Protein, %3.373.513.413.560.260.450.070.06
Protein yield, kg0.620.640.650.670.050.060.370.44
Lactose yield, kg0.930.921.0110.09<0.010.680.76
Milk solids, kg1.431.481.491.530.110.090.170.2
Total Solids, kg2.492.492.662.660.19<0.010.890.38
MUN, mg/dl13.219.411.815.90.17<0.01<0.010.7
1 Milk solids = Fat yield + protein yield; Total solids = total production without water; MUN = milk urea nitrogen. 2 RSD = Residual standard deviation.
Table 4. Effect of diet (ryegrass (RG) and plantain (PL)) and MUNBV (milk urea N breeding value) classified as either ‘high’ or ‘low’ on rumen parameter concentrations for dairy cows housed in metabolism crates. There were no significant interaction terms between day or any other factor.
Table 4. Effect of diet (ryegrass (RG) and plantain (PL)) and MUNBV (milk urea N breeding value) classified as either ‘high’ or ‘low’ on rumen parameter concentrations for dairy cows housed in metabolism crates. There were no significant interaction terms between day or any other factor.
Day 1Day 12
High MUNBVLow MUNBVHigh MUNBVLow MUNBVp-Value
Item 1PLRGPLRGPLRGPLRGRSD 2DayMUNBVDietMUNBV × Diet
Glucogenic21.321.120.920.817.517.417.217.10.08<0.010.530.790.41
Non-Glucogenic76.376.376.676.6808080.480.40.01<0.010.470.960.45
A:P ratio3.313.343.433.474.134.184.324.380.08<0.010.180.690.35
NH37.296.797.166.685.535.245.455.170.440.190.940.750.89
Total VFA73.872.669.568.272.371.167.966.720.790.830.560.870.35
1 A:P = Acetate: Propionate ratio; NH3 = Ammonia; VFA = Volatile Fatty Acids. 2 RSD = Residual standard deviation.
Table 5. Mean values for purine derivatives and creatinine and microbial nitrogen dynamics for dairy cows classified as either ‘high’ or ‘low’ for a milk urea nitrogen breeding value on either a plantain or ryegrass diet housed in metabolism crates.
Table 5. Mean values for purine derivatives and creatinine and microbial nitrogen dynamics for dairy cows classified as either ‘high’ or ‘low’ for a milk urea nitrogen breeding value on either a plantain or ryegrass diet housed in metabolism crates.
Item 1High MUNBVLow MUNBV p-Value
PlantainRyegrassPlantainRyegrassRSD 2MUNBVDietMUNBV × Diet
Purine derivatives (mmol/L)
Creatinine1.952.811.892.690.220.27<0.010.95
Hippuric22.5941.0819.9533.370.350.08<0.010.62
Allantoin4.997.814.87.370.240.42<0.010.65
Uric0.390.430.430.480.030.260.340.28
Total PD5.348.275.268.070.240.52<0.010.72
Microbial N dynamics
MCP flow, g/d1058 b1065 b1041 b1240 a149.230.090.030.04
Microbial
efficiency. g of MCP/kg TDOM
96.590.3107100.721.760.130.360.2
a–b Mean values in the same row with different superscript differ (p < 0.05) for the interaction between diet and MUNBV group. 1 PD = Purine derivatives (sum Allantoin + Uric acid); MCP = Microbial Crude Protein; TDOM = Truly digestible organic matter. 2 RSD = Residual standard deviation.
Table 6. Mean values for urinary and faecal characteristics for dairy cows classified as either ‘high’ or ‘low’ for a milk urea nitrogen breeding value on either a plantain or ryegrass diet housed in metabolism crates.
Table 6. Mean values for urinary and faecal characteristics for dairy cows classified as either ‘high’ or ‘low’ for a milk urea nitrogen breeding value on either a plantain or ryegrass diet housed in metabolism crates.
Item 1High MUNBVLow MUNBV p-Value
PlantainRyegrassPlantainRyegrassRSD 2MUNBVDietMUNBV × Diet
Urination characteristics
UUN, g/L1.69 c4.31 a1.86 c3.29 b0.38<0.01<0.01<0.01
UN, g/L2.50 c5.79 a2.70 c4.49 b0.32<0.01<0.01<0.01
UUN, g/event4.53 c11.90 a4.44 c8.62 b0.47<0.01<0.01<0.01
UN, g/event6.71 c16.20 a6.50 c11.66 b0.04<0.01<0.01<0.01
Volume/event, L2.632.832.472.650.340.01<0.010.16
Events/day15.1 a9.1 c15.2 a11.7 b0.220.06<0.01<0.01
Volume/day, L40.4 a26.4 c39.3 a31.0 b0.230.38<0.010.05
UUN, g/d67.5109.1268.82102.040.230.17<0.010.26
UN, g/d108.1147.83103.76139.850.210.23<0.010.35
UUN, % UN65.7 c75.3 a66.5 c73.5 b4.420.17<0.01<0.01
Urine patch area, m20.530.570.50.530.340.020.010.15
UUN loading rate, kg UUN ha−190.0 c241.1 a92.3 c173.7 b44.45<0.01<0.01<0.01
Faecal characteristics
N, g/event6.65 b5.51 c6.10 b7.75 a0.42<0.010.19<0.01
OM, g/event238.65 b183.10 d216.28 c265.22 a0.39<0.010.05<0.01
ADF, g/event105.47 a60.14 d97.72 b85.66 c0.39<0.01<0.01<0.01
NDF, g/event144.07 a92.23 c133.44 b133.63 b0.4<0.01<0.01<0.01
DM, g/event272.11 b221.35 d243.75 c320.07 a0.4<0.010.12<0.01
Events/day17.65 a17.81 a17.92 a12.71 b0.230.040.040.03
Total Faecal N/d117.66102.5114.69100.240.180.570.010.47
Total daily, kg/d4.74.124.543.990.010.37<0.010.12
a–c Mean values in the same row with different superscript differ (p < 0.05) for the interaction between diet and MUNBV group. 1 UUN, Urinary Urea Nitrogen; UN, Urinary Nitrogen; N, Nitrogen; OM, Organic Matter; ADF, Acid Detergent Fibre; NDF, Neutral Detergent Fibre; DM, Dry Matter. 2 RSD = Residual standard deviation.
Table 7. Average daily water balance for dairy cows on either a ryegrass or plantain diet classified as either ‘high’ or ‘low’ for MUNBV (milk urea N breeding value) housed in metabolism crates.
Table 7. Average daily water balance for dairy cows on either a ryegrass or plantain diet classified as either ‘high’ or ‘low’ for MUNBV (milk urea N breeding value) housed in metabolism crates.
High MUNBVLow MUNBVp-Value
ItemPlantainRyegrassPlantainRyegrassRSD 1MUNBVDietMUNBV × Diet
Water consumed, L/d
Trough0.830.39.435.59.830.03<0.010.21
Forage100.562.894.960.60.120.12<0.010.87
Total101.393.1104.396.113.460.470.060.33
Water excreted, L/d
Urine40.4 a26.4 c 39.3 a31.0 b0.230.38<0.010.05
Faeces29.834.427.532.24.280.09<0.010.88
Milk16.115.716.716.31.560.30.450.76
Total87.477.685.279.20.10.83<0.010.28
Water balance13.915.519.116.98.190.340.360.63
Water
absorbed, %
Prop. of intake70.162.973.566.34.50.02<0.010.29
a–c Mean values in the same row with different superscript differ (p < 0.05) for the interaction between diet and MUNBV group. 1 RSD = Residual standard deviation.
Table 8. Average daily intake and nitrogen balance for cows on either a plantain or ryegrass diet and classified as either ‘high’ or ‘low’ for MUN (milk urea N breeding value) housed in metabolism crates.
Table 8. Average daily intake and nitrogen balance for cows on either a plantain or ryegrass diet and classified as either ‘high’ or ‘low’ for MUN (milk urea N breeding value) housed in metabolism crates.
High MUNBVLow MUNBVp-Value
Item 1PlantainRyegrassPlantainRyegrassRSD 2MUNBVDietMUNBV × Diet
Intake
DMI, kg/d15.515.814.915.30.140.410.570.96
ME, MJ/d1791901721831.920.410.190.95
N, g/d338.5390.1327.6379.256.340.25<0.010.59
Nitrogen excreted
Milk, % of intake30.126.832.9290.20.150.040.69
Urine, % of intake33.138.431.536.40.260.40.020.32
Faeces, % of intake35.427.335.827.50.20.85<0.010.55
Total N excreted, g/d341.4374.2329.5362.836.050.47<0.010.32
N Balance, g/d−2.915.9−1.916.453.50.710.630.96
1 DMI, Dry Matter Intake; ME, Metabolizable Energy; N, Nitrogen. 2 RSD = Residual standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marshall, C.J.; Beck, M.R.; Garrett, K.; Barrell, G.K.; Al-Marashdeh, O.; Gregorini, P. Nitrogen Balance of Dairy Cows Divergent for Milk Urea Nitrogen Breeding Values Consuming Either Plantain or Perennial Ryegrass. Animals 2021, 11, 2464. https://doi.org/10.3390/ani11082464

AMA Style

Marshall CJ, Beck MR, Garrett K, Barrell GK, Al-Marashdeh O, Gregorini P. Nitrogen Balance of Dairy Cows Divergent for Milk Urea Nitrogen Breeding Values Consuming Either Plantain or Perennial Ryegrass. Animals. 2021; 11(8):2464. https://doi.org/10.3390/ani11082464

Chicago/Turabian Style

Marshall, Cameron J., Matthew R. Beck, Konagh Garrett, Graham K. Barrell, Omar Al-Marashdeh, and Pablo Gregorini. 2021. "Nitrogen Balance of Dairy Cows Divergent for Milk Urea Nitrogen Breeding Values Consuming Either Plantain or Perennial Ryegrass" Animals 11, no. 8: 2464. https://doi.org/10.3390/ani11082464

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop