1. Introduction
Utilizing homegrown forages with a high protein concentration increases the cost-effectiveness and sustainability of high-producing livestock systems. Previous studies highlighted the significant potential, but inconsistent results of whole-plant soybean silage (WSS) for ruminant nutrition [1] [2] [3]. Nevertheless, whole-plant soybean fermentation has limitations and challenges, such as low water-soluble carbohydrates (WSC) and dry matter (DM) concentrations, as well as high buffering capacity [4]. These characteristics impair silage fermentation by promoting proteolysis, and thereby, the production of poor-quality silage with high pH, ammonia-nitrogen (NH3-N), and butyric acid concentrations [5], leading to decreased ruminal nutrient degradation [6] [7].
Molasses is an additive with the potential to enhance WSS fermentation and ruminal degradation by increasing the availability of fermentable substrate to lactic acid bacteria. The addition of molasses (20 g/kg fresh forage) increased lactic acid concentration and reduced WSS pH [4]. Research literature [6] evaluating increasing doses of molasses (0, 30, 60, and 90 g/kg fresh forage) on WSS, observed a proportional increase of DM and lactic acid concentrations as well as a decrease of pH, NH3-N, and butyric acid concentrations in WSS, suggesting a dose-dependent response. Furthermore, greater ruminal degradation of DM, crude protein (CP), neutral detergent fiber inclusive of residual ash (NDF), and non-fiber carbohydrates (NFC) were observed in beef cattle diets containing WSS with molasses addition [7]. Improvements in fermentation profile reduce losses of neutral detergent solubles, which are more digestible than cell wall constituents [5], increasing ruminal degradability.
Moreover, the phenological stage at harvest influences the nutritive value of WSS mainly due to an increase in the proportion of pods and decrease in the proportion of leaves with later phenological stage [8]. These changes could affect the concentrations of digestible and indigestible fractions of WSS. A linear decline in dry matter intake and milk yield, with no effect on milk fat concentration was observed when 50% of forage DM (corn silage) was replaced with WSS harvested at 50% of full seed stage (R5.5) [3]. These responses were associated with impaired WSS fermentation and increased concentration of indigestible NDF as diets were formulated to be isonitrogenous and isoenergetic. Likewise, there was a decrease of soluble fractions and effective ruminal degradability (ERD) of DM and CP for WSS ensiled at beginning maturity (R7) compared to full seed (R6) stage [9]. A wide range of phenological stages from beginning seed (R5) to R7 has been commonly used to define when to begin harvesting WSS [10]. However, which phenological stage supports the greater nutrient degradation in the rumen is still unknown.
The effect of molasses to stimulate fermentation and increase nutrient availability in different phenological stages is not extensively investigated in the literature. We hypothesized that molasses would improve ruminal degradation of DM, CP, and NDF regardless of the phenological stage at harvest. Thus, the objective of this retrospective study was to evaluate the effects of molasses on ruminal DM, CP, and NDF degradation parameters, ERD, and undegraded neutral detergent fiber (uNDF; NDF basis) concentration of WSS harvested at R5, R6, and R7 phenological stages.
2. Material and Methods
2.1. Crop Establishment, Harvesting, and Ensiling
The first phase of the experiment (original study) [11] was carried out at the Dairy Cattle Research Laboratory of the University of São Paulo (LPBL/USP, Pirassununga, Brazil), and the second phase (current retrospective study) was performed at the Dairy Research Unit and Department of Animal Sciences of University of Florida (DRU/UF, Gainesville, FL, USA).
Soybean (Glycine max cultivar Agroeste® 3610 I PRO, Xanxerê, Brazil) was seeded in an experimental field localized 21˚58’19.8" south and 47˚28’15.9" west, at 620 m of altitude, subtropical humid (Cwa), and 1200 mm of mean annual rainfall. Soybean plants were harvested at three different phenological stages (R5: beginning seed, R6: full seed, and R7: beginning maturity [12] at 90, 113, and 124 days after seeded respectively). For each stage, 30 kg of forage material was harvested from each of five different plots previously determined in the field, hand-mixed to get a composite sample of 150 kg of fresh forage for each stage and was chopped (Cremasco® ECT-4000) targeting 20 mm of the theoretical length of cut. Fresh forage samples of each plot and phenological stage were collected for nutrient characterization.
A dose of 40 g/kg fresh forage of molasses was top-dressed on 2400 g of the composite sample and hand-mixed before placing into experimental silos, while the control group was not treated with molasses. The experimental silos (10 silos/treatment - 60 in total; silos were not used as experimental units in the current retrospective study) were made with polyvinyl chloride pipe with 0.5 m height and 0.1 m i.d. to achieve 650 kg/m3 fresh forage. All silos were sealed with plastic canvas (sailcloth 200 µm; Superlona black/white, Electro Plastic, Varginha, Brazil), taped within 2 h after harvest, and stored for 90 d.
2.2. Chemical Analysis of Whole-Plant Soybean Fresh Forage
For chemical profile characterization, fresh forage samples were individually dried at 55˚C for 72 h in a forced-air oven and ground to pass through a 1-mm sieve. The 1-mm ground samples were analyzed for absolute DM (method 950.15; correction of DM at 105˚C) [13], CP (N × 6.25; Kjeldahl method 984.13) [14], ether extract (EE; method 920.39) [14], NDF, acid detergent fiber (ADF), and lignin (sulfuric acid method) concentrations by sequential analysis [15] using a heat stable α-amylase [16] and sodium sulfite (TE-149 fiber analyzer, Tecnal Laboratory Equipment Inc., Piracicaba, Brazil). Non-fiber carbohydrate (NFC) concentration was calculated as NFC = 1000 − (NDF + CP + EE + ash) [17].
Chemical composition of whole-plant soybean fresh forage (average weight g/kg DM) is in Table 1. As expected, concentration of DM increased from 253 to 353 g/kg DM, as phenological stage progressed. Organic matter concentration ranged from 925 to 937 g/kg DM, as was very similar among phenological stages. Concentration of NDF and ADF decreased from 610 to 541 and 462 to 415 g/kg DM, respectively, as phenological stage progressed from R5 to R7. In addition, CP and EE concentrations increased with phenological stage.
2.3. Chemical Analysis of Whole-Plant Soybean Silage
At opening, the top 10 cm of WSS was discarded from each experimental silo and the remaining sample was homogenized for further analysis. For fermentation profile characterization, 15 g from each experimental silo was mixed with 250 mL of distilled water for 30 s and the silage juice was evaluated for pH and buffering capacity [18] using a pH meter (MB-10, Marte Cientifica®, Santa Rita do Sapucaí, Brazil). Silage juice was centrifuged (2000 × g for 15 min.), and the supernatant aliquot (0.9 mL) was mixed (1:1) with 450 µL sulfuric acid (1 N), for ammonia nitrogen (NH3-N) concentration analysis by colorimetric phenol-hypochlorite method [19]. For lactic acid concentration, a supernatant aliquot (1 mL) was frozen until further analysis [20]. For the concentration of organic acids, a supernatant aliquot (1.6 mL) of silage juice was mixed (4:1) with 0.4 mL ortho-fosforic acid (20%) and frozen until analysis. Briefly, organic acids were determined using a gas chromatograph equipment (GC-2014 Plus, Shimadzu©, Kyoto, Japan) equipped with automatic injector (AOC-20i at 300˚C), capillary column (Stabilwax-DA™: 30 m length, and 0.25 mm i.d.; Restek©,
Table 1. Chemical composition of whole-plant soybean (mean ± standard deviation g/kg DM).
1R5: beginning seed; R6: full seed; R7: beginning maturity; *not detected.
Bellefonte, USA) and flame ionization detector using helium as the carrier gas with a linear speed of 42 cm3/s. The temperature ramp of the column started with a gradient from 40˚C to 120˚C, at a rate of 40˚C/min, followed by ranges of 120˚C to 180˚C, at 10˚C/min, and from 180˚C to 240˚C, at 120˚C/min, holding temperature at 240˚C for 3 min at the end of analysis. Data were analyzed using the GCsolution software (v.2.42.00; Shimadzu Corporation, Kyoto, Japan).
Remaining samples from each experimental silo were dried at 55˚C for 72 h in a forced-air oven and ground to pass through a 2-mm sieve and were composited per treatment (6 samples in total, 1 for each treatment used in the present trial) for the DM, CP, and NDF in situ degradability assays. Composited 2-mm ground samples were ground to pass through a 1-mm sieve and analyzed for DM (method 950.15) [13], CP [N × 6.25; total N was analyzed by the Dumas dry combustion method (method 968.06) [13] using a CHNS analyzer (Vario Micro Cube; Elementar, Hanau, Germany)], and NDF (method 2002.04) [13] in an ANKOM 200 Fiber Analyzer (Ankom Technologies, Macedon, NY).
The addition of molasses slightly altered chemical concentrations of DM, NDF, CP, and lignin throughout the late compared to early phenological stages, whereas greatly changes were observed for fermentation parameters regardless of the phenological stage (Table 2).
2.4. In Situ Assays
Ruminal in situ procedures were conducted at the University of Florida (Gainesville, FL) under a protocol approved by the University of Florida, Institute of
Table 2. Chemical composition and fermentation profile (mean ± standard deviation) of whole-plant soybean silage harvested at different phenological stages, with or without the addition of molasses. (n = 10).
1R5: beginning seed. 2R6: full seed. 3R7: beginning maturity. 4Control: no additive. 5Molasses: 40g/kg fresh forage. 6Dry matter. 7Neutral detergent fiber. 8Crude protein. 9 meq. of alkali required of changing the pH from 4 to 6 per 100 g of dry matter [18].
Food and Agricultural Sciences, and Animal Care Research Committee (protocol #201709849). A completely randomized design with a 3 × 2 factorial arrangement of treatments was used considering phenological stage (R5, R6, and R7) and the addition of molasses (0 and 40 g/kg fresh forage) as main effects and rumen environment (cow) as experimental units. Ruminal in situ measurements were performed using three rumen-cannulated, mid-lactation, multiparous Holstein cows fed ad libitum a diet containing (DM basis) corn silage (382 g/kg), alfalfa hay (40 g/kg), dry ground shelled corn (273 g/kg), soybean meal (145 g/kg), citrus pulp (91 g/kg), minerals and supplements (68 g/kg). Cows were consuming this diet for at least 21 d prior to the experiment.
Dacron polyester cloth bags (R1020, 10 × 20 cm and 50 ± 10 μm pores; Ankom Technology, Macedon, NY, USA) containing 5 g of dried and 2-mm ground composited WSS samples (6 samples in total, 1 for each treatment), yielding a ratio of sample mass to bag area of 16.6 mg/cm2 were incubated in duplicate within each cow. A total of 96 bags were incubated within each cow (6 treatments in duplicated for each of the 8 time-points). Average of duplicated bags of each treatment generated the kinetic parameter estimates on each rumen environment (cow).
Bags were placed in a nylon laundry bag (30 × 40 cm) and incubated in the ventral rumen of each cow in reverse chronological order for 240, 120, 72, 48, 24, 12, 6, and 3 h. Additionally, two empty bags for each time point (16 bags in total/cow) were incubated to correct bag infiltration or losses. After the incubation period, bags were removed from the rumen, soaked in cold water to stop fermentation, and rinsed in a washing machine using the mode rinse and spin cycle, set with room temperature water for a 30-min cycle (Roper RTW4516F*, Whirlpool Corp., Benton Harbor, MI). Each treatment had two bags prepared, soaked, and washed together with the rest of in situ bags (0 h bags). After washing, bags were dried in a forced-air oven at 60˚C for 48 h, and the difference between the initial and final DM amount was considered as DM disappearance. Rumen incubated residues from duplicates within each treatment, timepoint and cow were composited, ground to pass through a 1-mm Wiley mill screen, and analyzed for CP (0, 3, 6, 12, 24, 48, and 72 h) and NDF (0, 3, 6, 12, 24, 48, 72, 120, and 240 h) concentrations as previously described (method 968.06 [16] and method 2002.04 [16], respectively). The NDF residue at 240 h was considered as uNDF concentration (DM basis).
The degradation of DM, CP and, NDF were calculated without correction for microbial protein contamination and calculated by the exponential model according to the Equation (1) [21], and the parameters estimated by PROC NLIN of SAS (version 9.3, SAS Institute Inc., Cary, NC, USA) as follows:
, (1)
where
is degradation (% DM) at time t; A is an intercept representing the proportion of soluble fraction (% DM); B is the insoluble but potentially rumen degradable fraction (% DM);
is the fractional rate of soluble degradation of the fraction B (h); t is the incubation time (h).
The effective ruminal degradability (ERD) of DM, CP, and NDF were calculated according to the Equation (2) [21],
(2)
where A, B, and
are the same parameters as described earlier, and
is a passage rate set to 3 %/h [21].
2.5. Statistical Analysis
Data on fractions A, B, and C,
, ERD,
, and uNDF concentration were analyzed as a completely randomized block design with the three cows designated as blocks. Average of duplicated bags of each treatment that generated the kinetic parameter estimates on each block (cow) were the experimental units. Analysis was performed by a mixed model procedure (PROC MIXED in SAS; SAS Institute Inc., Cary, NC, United States) according to the following model:
,
with
and
; where
is the value of the dependent variable, µ is the overall mean,
is the fixed effect of the phenological stage (i = 1 to 3),
is the fixed effect of molasses (j = 1 and 2),
is the random effect of cow (k = 1 to 3),
is the fixed interaction effect, and
is the residual error. N stands Gaussian distribution,
is the residual variance, and
is the variance associated with the random effect of cow. Degrees of freedom were adjusted based on Kenward-Roger’s method, an approach specifically proposed for small sample setting as a precaution for inflated error variance [22]. Significance was declared at P ≤ 0.05.
3. Results
3.1. Dry Matter Kinetic Parameters
Molasses did not affect (P ≥ 0.15) fraction B or the of DM, but both were greater (P = 0.001) at R7 compared to earlier phenological stages (Table 3). A molasses by phenological stage interaction was detected (P ≤ 0.03) on fractions A and C, and ERD of DM. Addition of molasses increased (P ≤ 0.05) fraction A regardless of phenological stage, but differences among phenological stages were only observed when WSS was not treated with molasses. Molasses addition increased (P ≤ 0.05) ERD while decreasing (P ≤ 0.05) fraction C of DM at R5 and R6, but not R7 phenological stage.
3.2. Crude Protein Kinetic Parameters
Molasses enhanced (P ≤ 0.05) the ERD of CP at R5 and R6, but not at R7 (Table 4). The addition of molasses increased (P = 0.01) fraction A and decreased (P =
Table 3. Effects of molasses and phenological stage on ruminal DM degradation kinetic parameters and effective ruminal degradability of whole-plant soybean silage (% DM, unless otherwise stated).
1A: r a-cMeans in the same row with different superscripts differ (P < 0.05). uminal rapidly degradable fraction; B: ruminal slowly degradable fraction; C: ruminal indigestible fraction; kd: degradation rate constant of fraction B; ERD: effective ruminal degradability at ruminal passage rate of 0.03 h−1 by Ørskov and McDonald [21]. 2Control: no additive; 3Molasses: 40 g/kg as-fed forage. 4R5: beginning seed; 5R6: full seed; 6R7: beginning maturity. 7MOL: molasses, PS: phenological stage, and INT: molasses and phenological stage interaction effects.
Table 4. Effects of molasses and phenological stage on ruminal CP degradation kinetic parameters and effective ruminal degradability of whole-plant soybean silage (% CP, unless otherwise stated).
a-bMeans in the same row with different superscripts differ (P < 0.05). 1A: ruminal rapidly degradable fraction; B: ruminal slowly degradable fraction; C: ruminal indigestible fraction; kd: degradation rate constant of fraction B; ERD: effective ruminal degradability at ruminal passage rate of 0.03 h−1 by Ørskov and McDonald [21]. 2Control: no additive; 3Molasses: 40 g/kg of as-fed forage. 4R5: beginning seed; 5R6: full seed; 6R7: beginning maturity. 7MOL: molasses, PS: phenological stage, and INT: molasses and phenological stage interaction effects.
0.01) fraction B, with no differences (P ≥ 0.43) on fraction C or the of CP. Later phenological stages decreased (P = 0.001) fraction A and increased (P ≤ 0.01) fraction B and kd, with no effect (P = 0.18) on fraction C of CP. There were no interaction effects (P ≥ 0.11) on fractions A, B, C, or kd of CP.
3.3. Neutral Detergent Fiber Kinetic Parameters
The addition of molasses decreased (P = 0.01) fraction C, increased (P = 0.001) ERD, and had no effect (P = 0.51) on the kd of NDF (Table 5). Harvesting WSS at a later phenological stage decreased (P = 0.001) fraction C and kd, increasing (P = 0.001) the ERD of NDF. Moreover, an interaction between the addition of molasses and phenological stage was detected (P ≤ 0.04) on fractions A and B of NDF, and uNDF. Fractions A and B of NDF were greater (P ≤ 0.05) at R7 and
Table 5. Effects of molasses and phenological stage on ruminal NDF degradation kinetic parameters and effective ruminal degradability of whole-plant soybean silage (% NDF, unless otherwise stated).
a-dMeans in the same row with different superscripts differ (P < 0.05). 1A: ruminal rapidly degradable fraction; B: ruminal slowly degradable fraction; C: ruminal indigestible fraction; kd: degradation rate constant of fraction B; ERD: effective ruminal degradability at ruminal passage rate of 0.03 h−1; undegraded NDF at 240 h of ruminal incubation (g/kg DM) by Ørskov and McDonald [21]. 2Control: no additive; 3Molasses: 40 g/kg of as-fed forage. 4R5: beginning seed; 5R6: full seed; 6R7: beginning maturity. 7MOL: molasses, PS: phenological stage, and INT: molasses and phenological stage interaction effects.
R5, respectively, while uNDF was decreased (P ≤ 0.05) at R5 and R6 when silages were treated with molasses.
3.4. Degraded Dry Matter, Neutral Detergent Fiber, and Crude Protein
Degraded DM, NDF, and CP concentrations in WSS at each incubation timepoint are presented in Figure 1.
4. Discussion
A positive effect of phenological stage was observed on ruminal degradation parameters of DM, CP, and NDF. However, the magnitude of these effects was lower in molasses-treated silages. The reason for that is related to the chemical composition of soybean plant in late stages which is closer to ideal for a good fermentation than early stages. While the phenological stage effect reflects morphological changes in the plant structure, such as leaves, pods, stalks, and their chemical compositions [8] [23], the effect of molasses could be related to either a reduction in moisture or a greater concentration of substrate for silage fermentation [4] [6]. Although both DM and WSC concentrations are crucial for obtaining an efficient silage fermentation [5], our findings suggest these were not limiting factors to produce a high ruminal degradable WSS at R7 stage. However, concentration of WSC was not measured in the current study and our speculation throughout this section should be considered cautiously. Greater DM increments in early phenological stages (17.2%, 15.7%, and 7.8% respectively to R5, R6, and R7) associated with the increase of 28.0%, 14.0%, or 13.7% in fraction A of DM at R5, R6, and R7, respectively, due to the addition of molasses suggests the additive was effective at increasing WSC concentration [24].
Additionally, fraction C of DM was decreased in molasses-treated silages at R5 and R6 indicating a dilution effect on indigestible compounds by supplying greater nutrient availability [4]. Changes in ruminal DM and nutrient kinetic parameters related to the chemical composition of the forage source have been reported previously [25] [26]. Furthermore, NDF degradation is closely related to DM disappearance.
At later phenological stages, there is a greater proportion of pods than stalk in soybean plants [8] [23], thereby increasing the proportion of more digestible plant compounds [8] [9]. Prior studies reported an increase in NDF and CP ruminal degradability of WSS when phenological stage progressed from R4 up to R6 [8]. In contrast, a decrease in WSS DM and CP ruminal degradability was reported at R7 in comparison to R6 [9]. In the present study, while the ERD of DM and CP were greater at R5 and R6 in molasses-treated silages, ERD of NDF has increased both with later phenological stage or the addition of molasses.
Greater solubility and rate of degradation of sugars brought into the silage by the addition of molasses are responsible for accelerating bacterial attachment and hence increase the CP soluble fraction regardless of phenological stage [4]. On the other hand, the same effect was observed in early phenological stage where the increase in fraction A of CP is at least in part due to greater NH3-N commonly found in high-moisture legume silage (<300 g/kg DM) [27]. It is worth to highlight, that although NH3-N averages were slightly lower than expected [27], the pattern throughout the treatments are coherent, reflecting the effects of phenological stages and molasses on fermentation profile properly. While greater levels of fraction A of CP improve rumen degradable protein to the microbial population in the rumen as percentage of non-protein nitrogen, fraction C contains proteins associated with lignin [28] that reduce enzymatic access, and portions of fraction B of CP show little, or no input. Thus, we speculate that diets formulated to ruminants containing WSS ensiled in early phenological stage (R5 and R6) or treated with molasses at ensiling require greater concentrations of ruminal degradable carbohydrates source to optimize animal performance [29].
Although high-moisture legume silages, such as R5 and R6 in the present study, are more susceptible to an extensive linkage of CP to cell wall components [30], no treatment effects on fraction C of CP were observed. However, the ERD of CP increased at R5 and R6 in molasses-treated silages, suggesting changes in the chemical composition of whole-plant soybean with later phenological stages [8] [23] and a shift of the CP fractions from rapidly to slowly degradable with greater kd, and lower neutral detergent insoluble nitrogen concentration [31] [32].
The uNDF decreased and ERD of NDF increased in WSS harvested at later phenological stages. This was related to the greater proportion of pods, and hence of the NDF digestible fraction, in later phenological stages. This finding is essential, as the morphological distribution of lignin in legumes is challenging to rumen microbial access and degradation of cell walls [33]. Increasing the potentially digestible fraction would allow for a greater inclusion of WSS in ruminant diets. Alternatively, it would allow ruminants to maintain similar performance when fed diets balanced to have similar NDF concentration [3]. Both, phenological stage and the addition of molasses are practices that can be easily implemented by producers to manipulate the degradation of structural polysaccharides and thus, change the NDF pool size, potentially improving animal performance [34].
5. Conclusion
Although a positive effect of phenological stage on ruminal degradation parameters of DM, CP, and NDF was observed, the magnitude of these effects was less pronounced in molasses-treated silages. Thus, despite the more pronounced effects of molasses that were observed in R5 and R6, its inclusion in whole-plant soybean silage is beneficial regardless of phenological stage. Changes in the proportion of digestible nutrients either by later phenological stage or the addition of molasses are powerful strategies to improve ruminal DM and nutrient degradation of whole-plant soybean silage.
Acknowledgements
The authors acknowledge the University of São Paulo (Pirassununga, Brazil), Dairy Cattle Research Laboratory of the University of São Paulo (Pirassununga, Brazil), and University of Florida (Gainesville, USA) for providing infrastructure, staff, land, and animal support. Lucas Ghedin Ghizzi expresses his acknowledgement to São Paulo Research Foundation (FAPESP, São Paulo, Brazil; grant#2019/02653-8) for providing his scholarship.