Experimental Site
The study was conducted at the Horticulture and Agroforestry Research Farm (HARF) near New Franklin, Missouri, USA (39°01'"N, 92°75'64"W). This area is comprised of silt loam, well-drained, loess soils that roll along the upland ridgetops, and low order streams converge into the Missouri river. Mean annual precipitation during the experimental period was 1060 mm (with approximately 64% of the precipitation occurred between April to September). Mean temperature during July was 25.6°C and mean temperature during January was − 2.06°C. About 2000 m2 area was used for planting elderberry plants in four plots. Each plot had 4 rows of 25 m in length. Details regarding soil fertility, bed preparation, nutrient analysis before bed preparation, and planting details (row-to-row distance and plant-to-plant distance) are available in previous publication (McGowan et al., 2018).
Plant fodder harvest and ensiling
The study was conducted in 2 × 3 factorial design structure in a completely randomized design. Fodder was harvested at three time points after the harvest of berries in August (15 Sept, 6 Oct, and 30 Nov of 2023), and ensiled at immediately after cutting and chopping (38–47% DM) or after cutting, wilting, and chopping (47–54% DM). The target packing density was 240 kg DM/m3. Each treatment combination was replicated three times.
The elderberry canes (whole plant) were coppiced at 12.7 cm above the ground level. Harvested canes were passed through a commercial mulch chipper (Wallenstein BX42) to approximately 5 to 7.6-cm size. Chopped material was passed through a mesh to screen pieces greater than 7.6 cm in length. The screened material was divided into two, with one portion containing material that was ensiled immediately and other portion wilted on a tarp in open air and sunlight for approximately 3 hours before ensiling to target 55% DM. Target DM was estimated using quick method of microwave drying and measured with a balance scale (A&D Engineering, Model EK-30KL balance, Ann Arbor, MI, USA). Pre-ensiled subsamples were collected to determine DM and nutritional parameters prior to placing into the mini-silo for all treatment combinations, but not each replicate.
Eighteen laboratory polyvinyl chloride (PCV) silos were constructed (10.2 cm diameter × 29.2 cm height) with rubber end caps and metal brackets to secure the caps, which allowed for 2.4 L of elderberry material as described in (Nieman & Conway-Anderson, 2022). Silos were packed, sealed, and placed in dry indoor maintained at approximately 10°C during the fermentation period. After 60 days, the ensiled silos were opened, contents were removed, weighted, and samples were individually packed in Ziplock plastic bags, and frozen at -20°C until nutrient and fermentation analysis.
Determination of fruit yield
Determination of fruit yield
The field layout is described in detail in a previous publication (McGowan et al., 2018). Four plots of elderberry plants were established, with each plot further divided into four rows, each 25 m in length. In three plots, the elderberry plant was coppiced on three different dates (15 Sep, 6 Oct, and 30 Nov), while the fourth (control) plot was coppiced on the following year on March. The elderberry plants were spaced at 3 m apart between rows and 2.4 m apart within rows. To determine yield, berries were harvested on 6 August, 2024 when the fruit was fully set and ~ approximately 50% ripened. This should occur prior to any evidence of pest pressure that might artificially reduce yield data. Each block of silage treatment consists of 4 rows. From each row, the number of cymes were counted and snipped from each row as close as possible without leaving extra stem and weighed in bulk. Clippers were used to snip cymes, and they were weighed using a bucket and portable scale.
Nutritive Analysis
Pre-ensiled samples (150 g) were taken from each batch of wilted and unwilted elderberry fodder on each date for nutrient analysis, and another 150 g of pre-ensiled samples from each treatment combination were kept separately in a forced air oven (Blue M Electric, POM-1406F, Watertown, WI, USA) at 60°C for 48 hours to calculate DM content. The remaining 150 g of material was shipped to Dairy One Forage Laboratory (Ithaca, NY, USA) for nutrient analysis. Major nutrient analysis included CP, ADF, NDF, TDN, NEL, NEM, NFC, NEG, RFV, and minerals (Ca, P, Mg, K, Na, Fe, Zn, Cu, Mn, Mo, S, and Ash).
Crude protein was analyzed using the combustion technique with CN928 Carbon/Nitrogen determinator (AOAC, 2003), ADF and NDF were analyzed using the Ankom200 Fiber analyzer (ANKOM, 2020), and non-fiber carbohydrates were calculated with the formula: 100 - (CP + NDF + ether extract + ash). Similarly, TDN was calculated using the summative energy equation by Dairy One, while net energy for lactation, NEG, and NEM were also determined by Dairy One (NRC, 2001). The relative feed value was calculated by multiplying digestible dry matter by dry matter intake, and then dividing by 1.29 (Moore & Undersander, 2002). For mineral analysis, samples were digested in 50 mL MARS Xpress vessels (CEM Corporation, Matthews, NC, USA) using a MARS 6 Microwave Digestion System (CEM Corporation Matthews, NC USA), and mineral concentrations were measured using a Thermo iCAP Pro XP ICP (Thermo Fisher Scientific Inc. Waltham, MA, USA).
Fermentation Analysis
Post-ensiled fodder materials were sent to Rock River Laboratory (Watertown, WI, USA) to determine fermentation parameters 2,3-butanediol, pH, acetic acid, ammonia, butyric acid,1,2-Prop, ethanol, formic acid, lactic acid, and succinic acid. Ensiled samples were mixed with deionized water for pH determination. Volatile fatty acids (VFA) were extracted at a ratio of 1:10 (silage sample to deionized water), centrifuged, and the supernatant solution was combined with calcium hydroxide and copper sulfate. This mixture was centrifuged a second time, and the supernatant was analyzed using high-performance liquid chromatography equipped with a reverse-phase ion exclusion column and a refractive index detector (Waters Corporation, Milford, MA). The same supernatant collected during VFA production was used for ammonia-nitrogen analysis with the Skalar San +++ Segmented Flow SA 500 Analyzer (Skalar, Breda, The Netherlands) (Krom, 1980).
Statistical Analyses
Nutritional and fermentation parameters of silage were analyzed using the R statistical software program (R version 3.6.2) to examine the impact of unwilted or wilted, harvesting dates, and the interaction of wilting status and harvesting dates with a 2×3 factorial treatment in ANOVA. A Tukey test was performed to identify significant differences among treatments. Individual treatment means were summarized using Least-Squares Means (LSMEANS), and the level of significance was determined at P < 0.05, with tendencies toward significance between 0.10 and 0.05. Also, the measurement of elderberry fruit yield was compared among early, mid, and late sampling season plots, including control, using analysis of variance (ANOVA) in R statistical software.