The effect of ensiled paulownia leaves in a high-forage diet on ruminal fermentation, methane production, fatty acid composition, and milk production performance of dairy cows

Haihao Huang, Dorota Lechniak, Malgorzata Szumacher-Strabel, Amlan Kumar Patra, Martyna Koz?owska,,4, Pawel Kolodziejski, Min Gao, Sylwester ?lusarczyk, Daniel Petri? and Adam Cieslak,8*

Abstract Background: The use of industrial by-products rich in bioactive compounds as animal feeds can reduce greenhouse gas production. Paulownia leaves silage (PLS) was supplemented to dairy cows’ diet and evaluated in vitro (Exp. 1; Rusitec) and in vivo (Exp. 2, cannulated lactating dairy cows and Exp. 3, non-cannulated lactating dairy cows). The study investigated the PLS effect on ruminal fermentation, microbial populations, methane production and concentration, dry matter intake (DMI), and fatty acid (FA) proportions in ruminal fluid and milk.Results: Several variables of the ruminal fluid were changed in response to the inclusion of PLS. In Exp. 1, the pH increased linearly and quadratically, whereas ammonia and total volatile fatty acid (VFA) concentrations increased linearly and cubically. A linear, quadratic, and cubical decrease in methane concentration was observed with increasing dose of the PLS. Exp. 2 revealed an increase in ruminal pH and ammonia concentrations, but no changes in total VFA concentration. Inclusion of PLS increased ruminal propionate (at 3 h and 6 h after feeding), isovalerate, and valerate concentrations. Addition of PLS also affected several populations of the analyzed microorganisms. The abundances of protozoa and bacteria were increased, whereas the abundance of archaea were decreased by PLS. Methane production decreased by 11% and 14% in PLS-fed cows compared to the control in Exp. 2 and 3, respectively. Exp. 3 revealed a reduction in the milk protein and lactose yield in the PLS-fed cows, but no effect on DMI and energy corrected milk yield. Also, the PLS diet affected the ruminal biohydrogenation process with an increased proportions of C18:3 cis-9 cis-12 cis-15, conjugated linoleic acid, C18:1 trans-11 FA, polyunsaturated fatty acids (PUFA), and reduced n6/n3 ratio and saturated fatty acids (SFA) proportion in milk. The relative transcript abundances of the 5 of 6 analyzed genes regulating FA metabolism increased.Conclusions: The dietary PLS replacing the alfalfa silage at 60 g/kg diet can reduce the methane emission and improve milk quality with greater proportions of PUFA, including conjugated linoleic acid, and C18:1 trans-11 along with reduction of SFA.

Keywords: Dairy cow, Fatty acid composition in milk, Methane emission, Paulownia leaves

Introduction

Forage is one of the prerequisites to meet nutritional requirements of ruminants. Besides the forage like alfalfa silage, some by-products like pomegranate silage can be effectively utilized to protect the environment [1]. Trees of thePaulowniagenus are known for their rapid growth [2] and thus an abundant source of wood [3]. A vast biomass of leaves becomes available as a by-product of wood processing. Leaves are rich in protein (average 175 g/kg DM), and hence can be used as an alternative source of forage for ruminants [4]. The nutritional value of Paulownia has been investigated in several studies [5–7]. Paulownia leaves are rich in bioactive compounds (BAC) such as phenolic acids and flavonoids as well as in fatty acids (FA) [4]. Besides, ensilaged paulownia leaves contain a higher concentration of BAC compared to its fresh material [4]. Our previous in vitro experiment showed that inclusion of either fresh or ensiled paulownia leaves lowered methane production by reducing the number of methanogens while improving the basic ruminal fermentation characteristics [4, 8]. Considering a high phenolic content of paulownia leaves, they can be a valuable component of dairy cow diet with an additional profit attributed to the improvement of the milk quality by altering its FA profile [4]. We hypothesized that paulownia leaves may beneficially alter ruminal fermentation variables such as pH or methane production, ammonia and volatile fatty acid concentration due to presence of phenolic compounds in PLS. Besides, the inclusion of PLS into the diet of dairy cows may affect milk composition by altering its FA profile by modifying ruminal biohydrogenation. The improvement of milk quality may result from an increase in n-3 PUFA concentration and thus a reduced n-6/n-3 ratio. However, the effect of dietary use of paulownia leaves silage (PLS) on milk production, methane mitigation, and ruminal fermentation has not been studied in dairy cows. To verify the hypothesis, we utilized the ensiled leaves ofPaulownia tomentosa×Paulownia fortuneihybrid plants in an experimental design, which included in vitro and in vivo conditions. We aimed at investigating whether PLS affects (i) methane production and microbial population in vitro as well as (ii) ruminal fermentation parameters, methane production, and milk production and composition in vivo in dairy cows.

Material and methods

Diets and supplements

For the preparation of PLS, paulownia leaves (along with some twigs) were lopped from the trees (three-yearold plantations) at the end of May when the trees were shaped. Leaves were wilted for 8 h and were ensiled using the biological additive (Agricol Sil, Microferm, UK) in 60 kg plastic drums (diameter 40 cm and height 61 cm) according to Huang et al. [4]. The alfalfa was harvested at the second cut at the beginning of June. The alfalfa fresh materials after wilting approximately to 35% of dry matter were chopped to a particle length circa 2.5 cm and ensiled using the same biological additive as paulownia silage. One gram of the additive contained 1011colony-forming units ofLactobacillus plantarumDSMZ 16,627 andPediococcus acidilacticiNCIMB 30,005, as well as an enzyme-producing strain,Lactobacillus paracaseiNCIMB 3015. The closed drums were stored for 8-9 weeks in the case of both plants. The chemical composition of PLS and AS is presented in Table 1. The following fermentation parameters of PLS and AS were observed: pH, 4.75 vs. 4.50; NH3-N, 79.4 vs. 32.3 g/kg total nitrogen; lactic acid, 21.2 vs. 34.9 g/kg DM; acetic acid, 7.04 vs. 19.1 g/kg DM; propionic acid, 1.29 vs. 1.09 g/kg DM; butyric acid, 1.19 vs. 2.43 g/kg DM, respectively.

In vitro experiment (Exp. 1)

The in vitro experiment was conducted using the Rusitec system equipped with four fermenters of 1 L volume each, following the procedures described by Szczechowiak et al. [9]. Rumen fluid and solid digesta for the in vitro experiment were collected 3 h before the morning feeding from four rumen-cannulated (Bar Diamond, Parma, Idaho, USA) multiparous Polish Holstein–Friesian dairy cows (630 ± 25 kg body weight) at their 3rdmonth of lactation. The following diets were tested: control diet (CON) and three PLS diets. The CON diet contained the following forages: corn silage (388 g/kg DM), alfalfa silage (82 g/kg DM), and meadow grass silage (91 g/kg DM; Table 1). The PLS diets contained paulownia silage that replaced alfalfa silage at 25%, 50%, and 75%, which corresponded to the PLS content in the diet at the level of 20, 40, and 60 g/kg DM (Table 1). Ruminal fluid donor cows received the CON diet twice a day ad libitum (Table 1). To meet the nutrient requirements of mid-lactation dairy cows (600 kg body weight, 110 d in milk, and 34 kg/d milk production, 4% of fat content in milk), the in vivo diets were prepared using the FeedExpert software (Rovecom, Hoogeveen, Netherlands). The in vitro experiment was designed in a completely randomized design comprising four diets and three replicates. From d 6 to 10 of each run, fermentation fluid samples were collected under anaerobic conditions from each vessel 3 h before feeding time. The collected fluid samples were analyzed for pH, volatile fatty acids (VFA), ammonia concentration, protozoa, bacteria and methanogen counts. For fatty acids analysis, samples were collected directly from the effluent vessels duringbag replacements. Before feeding time (once a day), the fermentation gas was collected in a gas-tight collection bag (Tecobag 81; Tesseraux Container, Bürstadt, Germany) for the methane concentration measurement. The DM degradability was determined by analyzing feed residues in pre-feeding nylon bag samples for the last 5 d of each run (d 6 to 10).

Table 1 Ingredients and chemical composition of experimental diets (n = 4) used in Rusitec system and in vivo experimetsa and chemical composition of paulownia leaves silage (PLS; n = 4) and alfalfa silage (AS; n = 4)

In vivo experiment using rumen?cannulated dairy cows (Exp. 2)

In Exp. 2, four multiparous cannulated Polish Holstein–Friesian dairy cows (625 ± 20 kg body weight; 4-5thmonth of lactation) were assigned to two dietary treatments (CON vs. PLS60) with two cows in each treatment in a replicated 2 × 2 crossover design. Based on the results from Exp. 1 (mainly pH, total VFA concentration, methane production, and methanogens population), the CON and PLS diets containing the higher level of PLS, i.e., 60 g/kg DM (PLS60) were implemented in Exp. 2. Cows were fed two times a day. Each period lasted for 36 d, with a 21-d adaptation (where cows were also introduced and accustomed to the conditions of the respiration chambers; temperature 12 to 22 °C and the humidity of 50% to 70%) and a 15-d sampling period (5 d of rumen fluid collection and 10 d for gases collections). The cows were housed in tie stalls with rubber mats with individual feeding and had free access to water and salt blocks during the adaptation and sampling (without the period in which the cows were kept in the respiration chambers). The ruminal fluid was collected from each cannulated cow from three locations (top, bottom, and middle) of the midventral sac of the rumen before morning feeding (0 h) and 3 h and 6 h after morning feeding [10]. Rumen samples (about 400 g/animal) were filtered through a two-layer cheesecloth and analyzed for pH value, ammonia and VFA concentrations, and FA profile. For protozoa counting, about 100 g/animal of the rumen content was mixed with an equal amount of 8% formaldehyde solution (w/w), strained through a twolayer cheesecloth into 10 mL polypropylene tube with a screw cup, and stored at 8 °C in a refrigerator until analysis. Quantification of total bacteria and methanogens was carried out only on rumen fluid sampled 3 h after morning feeding. For microbial analysis, rumen content (300 g) was strained through a two-layer cheesecloth into 100 mL polypropylene box, mixed, transferred into two cryotubes of 4.5 mL and frozen in liquid nitrogen. Samples were stored at ? 80 °C until further analyses.

Feed intake, feed residue, and amount of feces were recorded daily from individual cows kept in respiratory chamber (SPA System, Wroclaw, Poland) during the sampling period (d 27 to 36) to determine the total-tract degradability coefficients. The feed and feces subsamples (about 5% wt/wt) were stored at -20 °C for DM,organic matter (OM), neutral detergent fiber (NDF), and CP analyses. The total tract nutrient degradability was calculated as [(nutrient intake – nutrient in feces)/nutrient intake] × 100.

In vivo experiment using commercial dairy cows (Exp. 3)

In Exp. 3, 16 multiparous lactating Polish Holstein–Friesian dairy cows [600 ± 30.4 kg body weight, 2.4 ± 0.45 parity, 160 ± 32 d in milk, and 33 ± 2.1 kg/d milk production; (mean ± SD)] were assigned to two dietary treatments (CON vs. PLS60) with eight cows in each treatment in a replicated 2 × 2 crossover design. Each period consisted of a 21-d adaptation period followed by a 5-d sampling period with a total of 26 d. Cows were randomly assigned to one of the two dietary groups (n= 8): control (CON) and experimental diet (PLS 60) and kept separately in a dedicated area of a barn. The two groups (8 CON and 8 PLS60) had separated access (gate with access control) to a computer-controlled feeder station (De Laval, type FP 204, Tumba, Sweden). The concentrate supplement, contained (as g/kg of DM in concentrate) OM (910), aNDFom (240), CP (17.5), and EE (31), was supplied in both groups at 2.75 kg/d/cow (5 times a day). The rest of the diet was offered twice a day (at 06:00 and 18:00 h) as a partial mixed ration (PMR) in individual feeding boxes located on the feeding table. Due to the limitations of running the experiment under production conditions, the cows were supervised by two workers for 18 h, who ensured that each cow consumed feed from an appropriate feeding box (CON or PLS60). For the remaining 6 h, the feeding boxes were closed. The control PMR contained (g/kg of DM): corn silage (441); alfalfa silage (93), meadow grass silage (103); beet pulp (118), brewer’s grain (108); rapeseed meal (122); mineral and vitamin premix (16) whereas experimental PMR contained (g/kg of DM): corn silage (438); alfalfa silage (29), paulownia leaves silage (68), meadow grass silage (102); beet pulp (117), brewer’s grain (108); rapeseed meal (122); mineral and vitamin premix (16). The nutritive value of the two diets (CON and PLS60) was the same. Cows had access to clean water ad libitum. The dry matter intake was measured daily for the last 5 d (d 22 to 26) of the experiment by weighing the individual amounts of feeds offered and leftovers in the feeding boxes. Feces were individually collected from each cow directly after defecation and the floor was kept clean. Subsamples (about 5% wt/wt) were collected and stored at -20 °C until analysis. Total tract nutrient digestibility was calculated using the method described for Exp. 2. Cows were milked twice a day at 5:30 and 17:30 h in a herringbone milking parlor with 8 milking units, spending a maximum of 1 h/d outside the pen for milking. Milk samples were collected from all cows at each milking during the sampling period (d 22 to 26) based on the proportion of morning and evening yield. The morning milk samples were stored at 4 °C until the evening samples were taken. Subsequently, milk samples were prepared in three equal parts: first part was used to analyze milk basic constituents, second part was stored at -20 °C for FA analysis, and the third part was quickly stored in liquid nitrogen for gene expression analysis.

Sample analysis

Feed and feces samples were analyzed according to AOAC methods [11] for DM (method no. 934.01), ash (method no. 942.05), CP (Kjel-Foss Automatic 16,210 analyzer, Foss Electric, Hiller?d, Denmark; method no. 976.05), ether extract (EE; Soxhlet System HT analyzer; Foss Electric, Hiller?d, Denmark; method no. 973.18) and neutral detergent fiber (aNDF) with amylase and sodium sulfite and expressed without residual ash (Fibertech 1020 Analyzer; Foss., Analytical AB, H?gan?s, Sweden; method according Van Soest et al. [12]. Organic matter was obtained by subtracting ash from DM. Concentration of total phenolic compounds in the diets (Table 1) was calculated by utilizing data for ensiled paulownia leaves published by Huang et al. [4]. The same batch of ensiled paulownia leaves was used in the present study. The silage pH was determined using a pH meter (Elmetron, Type CP-104, Zabrze, Poland). Ammonia was estimated according to the Nessler method [13]. Lactic, acetic, propionic, and butyric acids were determined using a High Performance Liquid Chromatography (Waters 2690, Santa Clara, CA, USA) equipped with Waters 2487 Dualλdetector and Aminex HPX-87H column (300 mm × 7.8 mm, Bio-Rad, Warsaw, Poland). The quantitative and qualitative evaluations of individual peaks were made using the external standard method prepared by mixing individual fatty acids purchased from Supelco (Poznan, Poland) and analysed with the Millennium 2001 software (version 2.15, Waters Corporation, Manchester, England).

The pH of ruminal fluid was measured immediately after sample collection using a pH meter (Elmetron, Type CP-104, Zabrze, Poland). The ammonia concentration was determined using the colorimetric Nessler method and the VFA was analyzed using gas chromatography (GC Varian CP 3380, Sugarland, TX, USA) as described earlier [9, 13]. In the first two experiments, the content of DM degradability (g/kg) was determined by the difference between the initial feed substrate weight and residue weight after incubation.

Protozoa counts in the fermented fluid were carried out under a light microscope (Primo Star 5, Zeiss, Jena, Germany) using an appropriate volume (10 μL for Ophryoscolecidae and 100 μL for Isotrichidae). In Exp. 2, the protozoan genera and species were identified according to size and shape of cells, skeletal plates (if present), macronucleus, and arrangement of ciliature [14]. The methanogens and total bacteria were quantified by fluorescence in situ hybridization (FISH), following the procedure described previously [15] with some modification. Briefly, 6 mL of sterile paraformaldehyde-phosphate buff-ered saline (PBS; pH 7.2) was added to 2 mL of ruminal fluid in a stomacher bag. The mixture was homogenized in a stomacher (Interscience, Saint-Nom-la-Breteche, France) for 2 min. Afterwards, 2 mL of the homogenized mixture was transferred to the eppendorf tubes and fixed with 4% (w/v) sterile paraformaldehyde-PBS for 3 h at 48oC. Mild sonication (two times for 30 s using Hielscher Ultrasonics, Teltow, Germany) was done to avoid the formation of clusters and to optimize homogenization of the samples [16]. The prepared samples were pipetted onto 0.22 μm polycarbonate filters (Frisenette K02BP02500) and vacuumed (Vaccum KNF Vacuport-Neuberg). After vacuuming, the filters were transferred onto cellulose disks for dehydration in ethanol series (50%, 80%, and 90%, 3 min each). For each sample, a series of identical filters were prepared to allow the determination of optimal hybridization. Hybridizations were carried out in 50 μL of hybridization buffer (0.9 mol/L NaCl; 20 mmol/L Tris/HCl, pH 7.2; 0.01% SDS) containing oligonucleotide probes for all methanogens (S-D-Arch-0915-a-A-20) and two order-specific probes (S-O-Mmic-1200-a-A-21 Methanomicrobiales and S-F-Mbac-0310-a-A-22 Methanobacteriales) [17]. After hybridization, the filters were washed with a washing buffer (20 mmol/L Tris/HCl, pH 7.2; 0.01% SDS; 5 mmol/L EDTA) for 20 min at 48 °C. The filters were rinsed gently in distilled water, air-dried, and mounted on object glasses with VectaShield (Vector laboratories nr. H-1000) anti-fading agent containing DAPI (4?,6-diamidino-2-phenylindole). To distinguish the total count of bacteria (DAPI) from methanogens in the ruminal fluid, filters were maintained at 4 °C for one h in the dark until visualization using an Axio Imager M2 microscope (Carl Zeiss Iberia, Madrid, Spain). Besides bacteria counts, relative changes in the population of seven selected species (Ruminococcus flavefaciens, Fibrobacter succinogenes, Streptococcus bovis, Butyrivibrio proteoclasticus, Ruminococcus albus, Butyrivibrio fibrisolvens, Megasphaera elsdenii)and two genera (Prevotellaspp.,Lactobacillusspp.) of ruminal bacteria were determined by quantitative real-time PCR (qPCR). For this purpose, total DNA was extracted from the ruminal fluid using QIAamp DNA Stool mini kit (Qiagen GmbH, Hilden, Germany) according to Szczechowiak et al. [9]. Sequences of primers specific to bacterial genera and species are given in a Supplementary data (Table S1). The specificity of primers (Table S1) was confirmed in the GenBank Database using the BLAST program. The quantitative analysis of particular bacteria was performed with a known starting concentration of bacterial DNA (25 ng/μL) using the QuantStudio 12 Flex PCR system (Life Technologies, Thermo Fisher Scientific, Waltham, MA, USA) [9]. The Power SYBR Green PCR Master Mix (Thermo Fisher Scientific, Waltham, MA, USA) was used for PCR amplification. The reaction mixture (final volume of 10 μL) contained 4 μL of the 2 × Mastermix, 25 ng of template DNA and 0.5 mol/L of each primer. Amplification involved one cycle at 95 °C for 10 min for initial denaturation, 45 cycles of 95 °C for 15 s followed by annealing at temperatures depending upon the individual bacteria, and then primer extension at 60 °C for 62 s. The fluorescent product was monitored in the last step of each cycle. To determine the amplicon specificity, melting analysis was performed after a single amplification (0.1 °C × s?1increment from 65 °C to 95 °C with fluorescence collection at 0.1 °C intervals). Additionally, the size of amplicons was verified by gel electrophoresis. The size of total bacterial populations was referred to the calculated copy number concentrations of the 16S rRNA (rrs) gene [18]. The absolute abundance of bacterial DNA was expressed as a number ofrrsgene copies/mL of ruminal sample.

Table 2 The effect of paulownia leaves silage (PLS) on in vitro ruminal fermentation and methane production (n = 4) (Exp.1)

In the first experiment, the methane concentration was determined using gas chromatograph (SRI PeakSimple 310; Alltech, PA, USA) fitted with Carboxen 1000 column (Supelco, Bellefonte, USA) and thermal conductivity detector according to Szumacher-Strabel et al. [19]. In the in vivo experiments, the methane and carbon dioxide concentration were measured using two separate NDIR (nondispersive infrared spectroscopy) systems (one system per gas) operating in the near-infrared spectrum (detector 1210 Gfx Servomex 4100, Servomex, Crowborough, UK). In Exp. 2, two respiration chambers were used for monitoring of methane production. Briefly, two open-circuit respiration chambers (W × L × H: 300 cm × 400 cm × 220 cm; SPA System, Ltd., Wroclaw, Poland) were used to measure CH4and CO2over 10 d (d 27 to 36). During the sampling period (from d-27 to 36 of the experiment), individual cows were transferred into a respiratory chamber by daily rotation in order to determine the direct CH4emission for 23.5 h consecutively. The time of milking (approximately 30 min) that accompanied the morning and evening feedings was not included in the gas emission calculations. Finally, each cow was monitored for 5 d. During the adaptation period (from d 11 to d 20) to decrease stress, cows were getting accustomed to the respiratory chambers. Cows were restrained within the chambers by a neck yoke on a dedicated platform (180 cm × 126 cm) covered with a rubber mat and had free access to fresh water and salt blocks. Emission of CH4and CO2was determined using two NDIR analyzers operating in the near-infrared spectrum (SERVOMEX 4100, SERVOMEX Ltd, UK, detector 1210 Gfx). Measurements were taken every 2-s interval. Two measuring channels were used: the concentration of CO2in the range of 0–2.5% (0–48,450 mg/m3) and the CH4concentration in the range of 0–1000 ppm (0–706 mg/m3). Samples were collected and ducted to the analyzer via a polyethylene tube with an 8 mm diameter. The sampling rate was 0.6 L/min. Before starting the experiment, analyzers were calibrated using calibration gasses: nitrogen N 5.0 (99,999 vol % purity) and 1210 ppm CH4in nitrogen. The analyzer was equipped with 0.17 L cuvette with 540 mm optical track length for CH4and 0.012 L cuvette with 154 mm optical track length for CO2. In Exp. 3, methane and carbon dioxide concentrations were measured by the same type of detectors during the feeding of concentrate in the feeder station for the last 4 d of the experiment [20]. Air samples were continuously collected from the feed bins in the feeder station at 15 L/min via an 8-mm diameter polyethylene tube and connected to detectors.

The FAs concentrations in collected samples (feed, ruminal fluid, and milk) were determined by gas chromatograph (456-GC, Bruker, USA) with fused-silica capillary column (100 m × 0.25 mm; overlaid with 0.25 μm Agilent HP; Chrompack CP7420; Agilent Technologies, Santa Clara, CA, USA) and flame ionization detector [9].

Daily milk yields were recorded using a milk meter (WB Ezi-Test Meter 33 kg; True-Test, Manukau, New Zealand). The milk composition was measured by infrared analysis (MilkoScan 255 A/S N, FossElectric, Hiller?d, Denmark). For gene expression in milk somatic cells, total RNA was isolated from 10 mL milk samples frozen in liquid nitrogen using previously published procedure [9]. Briefly, milk samples after thawing at 4 °C were centrifuged in 15 mL tubes for 10 min at 3000 ×g. The supernatant was discarded, and the pellet was dissolved in 1 mL TriPure reagent (Roche) and incubated 5 min. After incubation, 200 μL of chloroform were added and shaken vigorously for 30 s. After 10 min incubation at room temperature, the sample was centrifuged for 15 min at 12,000 ×g. The clear phase was transferred to a new tube, 0.5 mL isopropanol was added and incubated for 10 min at room temperature. The next steps involved RNA precipitation with a 75% ethanol and drying on a 40 °C thermoblock. The RNA was resuspended in DEPC treated water following spectrophotometric measurement (Nanodrop c2000, Thermo Scientific, USA) of its concentration and purity. A reverse transcription reaction (RT) was performed using a Transcriptor First Strand cDNA Synthesis Kit (Roche) according to the manufacturer’s protocol. Each sample contained equal concentrations of RNA. The RNA mix composed of RNA (300 ng), random hexameters, oligodT (60 μmol/L and 2.5 mmol/L respectively) and water was incubated at 65 °C for 10 min. Next, reverse transcriptase, RNAse inhibitor, dNTP and buffer were mixed and added to the RNA mix to a final volume of 20 μL. The RT conditions were as follows: 25 °C for 5 min, followed by 42 °C for 45 min and 85 °C for 5 min. Resultant cDNA was stored at ? 20 °C until further analyses. The mRNA expressions of six genes encoding enzymes regulating FA metabolism [acetyl-CoA carboxylase 1 (ACACA), fatty acid synthase (FASN), lipoprotein lipase (LPL), stearoyl-CoA desaturase (SCD), fatty acid desaturase 1 (FADS1) and fatty acid elongase 5 (ELOVL5)] were measured in milk somatic cells using previously published primer pairs [21, 22]. Each gene was analyzed in technical duplicates using a LightCycler 480 instrument (Roche Diagnostics, Basel, Switzerland) and a LightCycler 480 Sybr Green I Master reagent (Roche Diagnostics, Basel, Switzerland).

Statistical analysis

The data of the Exp. 1 (Rusitec) were analyzed using a mixed model procedure (PROC MIXED) of SAS (university edition, version 9.4; SAS Institute, Cary, NC, USA). The dietary treatment was considered as the fixed effect, experimental run as the random effect, and the day (6 to 10 d) as the repeated factor. The linear, quadratic and cubic contrasts were used to determine the effect of PLS dose. In Exp. 2 (cannulated cows), ruminal fermentation and FA data were analyzed using PROC MIXED (ver. 9.4, SAS Institute Inc., Cary, NC) for a crossover design with a model containing group (dietary treatment sequence), period, and treatment as main effects, sampling time as repeated measures, and cow as a random effect. The model for the bacteria, methanogens and degradability analyses contained cow, group, period, and treatment as a main effect. In Exp. 3 (productive dairy cows), the data were subjected to analysis of variance, considering the crossover design, testing the effect of treatment (CON and PLS60), group (dietary treatment sequence), period as fixed effects and cows within a group as a random effect. The analysis of milk production, component yields, and composition was performed on the mean values of the milk variables obtained from two sampling points per day (morning and evening milking). The analysis of FA proportion in the milk as well as expression of six genes was conducted on the mean values of pooled samples from morning and evening milking. The results were tested with an independentt-test where the means of both groups were compared through PROC TTEST procedure. The results were considered significant when theP-values were lesser than 0.05. All values are shown as the means with pooled standard errors of means.

Results

Exp. 1 (Rusitec study)

Ruminal pH increased linearly and quadratically (P< 0.01) with increasing concentrations of PLS with the greatest pH at the highest PLS inclusion (Table 2). The inclusion of PLS in diets increased the ammonia concentration linearly and cubically (P< 0.01) and total VFA concentration linearly (P< 0.01). Molar proportion of acetate and acetate to propionate ratio decreased (P< 0.01) linearly; however molar proportions of other VFA (except isovalerate) increased linearly with increasing levels of PLS (P< 0.01). Degradability of nutrients was not affected by the PLS inclusion in diets. The total gas production increased (P< 0.01) linearly when PLS was included into the diet. Daily methane concentration decreased (P< 0.05) linearly, quadratically, and cubically with increasing doses of PLS in diets. Methane concentration per unit of degraded DM decreased linearly and quadratically (P< 0.05) as the PLS dose increased. The inclusion of PLS into the diet resulted in changes (P< 0.05) of most microbial populations in a linear pattern (Table 3). Total protozoa and Ophryoscolecidae populations decreased linearly (P< 0.01) whereas Isotrichidae counts increased linearly (P< 0.05) with increasing PLS levels. Total archaea, Methanobacteriales and Methanomicrobiales populations were decreased (P< 0.01) linearly by the higher replacement of alfalfa silage with PLS in the diet. The increasing PLS supplementation caused linear increase ofRuminococcus flavefaciens, Fibrobacter succinogenes, Streptococcus bovis,Prevotellaspp.,Butyrivibrio fibrisolvens,andMegasphaera elsdenii(P≤ 0.05) abundances. Besides,M. elsdeniishowed quadratic and cubical responses (P< 0.05).

In vivo experiment (using cannulated cows)

Inclusion of PLS resulted in an increase (P< 0.01) in pH and ammonia concentration in the rumen. Besides, pH values were post-feeding time dependent (P< 0.01) whereas ammonia concentration showed treatment × time interaction (Table 4). The inclusion of PLS decreased the molar proportion of acetate (P< 0.01), but increased the molar proportion of propionate (3 h and 6 h after feeding;P< 0.01), isovalerate (P< 0.05) and valerate (P< 0.01). Molar proportions of most individual VFA showed time-dependent variations (P< 0.05). A/P ratio was lower in the PLS group compared to the control group in 3 h and 6 h after morning feeding (P< 0.01). Ammonia concentration was affected (P< 0.01) by time × treatment interaction, which increased in the PLS diet compared with the control diet (P< 0.01) at 3 and 6 h, but not at 0 h.

Table 3 The effect of paulownia leaves silage (PLS) on in vitro ruminal microbial population (n = 4) (Exp.1)

Table 4 The effect on replacing alfalfa silage with paulownia leaves silage (60 g/kg) on ruminal fermentation characteristics measured in rumen-cannulated cows (n = 4) (Exp. 2)

Total Isotrichidae,Dasytricha ruminantium,Isotricha prostomaandIsotricha intestinalispopulations increased by feeding the PLS diet. The abundances were affected by treatment × time interaction (P< 0.01), except forIsotricha intestinalis.Ostracodinium gracileandPolyplastronmultivesiculatum(P< 0.01) also were affected by treatment × time interaction. The first three protozoa mentioned were increased by PLS at 3 h and 6 h, but were similar at 0 h; whereasOstracodinium gracilepopulation was greater at 3 h and 6 h, andP. multivesiculatumwas greater at 0, 3, and 6 h in the PLS diet than in the control diet. The inclusion of PLS in the diet decreased total protozoa, Ophryoscolecidae andEntodiniumspp. counts (P< 0.01).

Feeding PLS to cannulated cows increased (P< 0.05) the populations of all bacterial species examined in this study, exceptRuminococcus flavefaciens, Ruminococcus albusandLactobacillusspp. (Table 5). The decreases in abundances of total archaea, Methanobacteriales and Methanomicrobiales were noted in the experimental group (P< 0.01). The inclusion of PLS did not affect degradability of DM, OM, and NDF, but a lower crude protein degradability was observed (P< 0.01) in PLS diet. The EE degradability increased due to PLS feeding (P< 0.05). Lowered (P< 0.01) methane production (g/d) and yield (g/kg DM intake) in PLS diet were noted compared to the CON diet; however, PLS diet did not alter CO2production in the rumen-cannulated cows.

In the Exp. 2 (cannulated cows), the ruminal fluid of PLS cows was characterized by altered proportions of selected FAs. The level of C8:0, C10:0, C12:0, C16:0, C16:1, C17:1 and C18:0 decreased (P< 0.05; Table 6) whereas the level of C14:1, C:15:0, C:15:1, and C17:0 increased (P< 0.01). Regarding the proportion of C8:0, C12:0, C14:1, C:15:0, C:15:1, C16:0, and C16:1, the treatment × time interaction was also observed (P< 0.05). The proportion of C18:1trans-10, C18:1trans-11, C18:2cis-9,cis-12, C18:2cis-9,trans-11, and C18:2trans-10,cis-12 increased in both treatment (P< 0.01) and time-dependent (P< 0.01) manners in PLS treatments. The PLS diet tended to decrease C18:1cis-9 (P< 0.05), but increase C18:3n-6 and C18:3cis-9,cis-12,cis-15 (P< 0.01) in treatment-dependent. Sum of unsaturated fatty acids (UFA), monounsaturated fatty acids (MUFA), PUFA, n-6 FA and sum of n-3 FAincreased (P< 0.01) in the experimental group, which resulted in lower sum of SFA (P< 0.01). Time-dependent variation (P< 0.01) was observed in all of the above parameters, in addition to the sums from MUFA. The inclusion of PLS increased the sum oftransC18:1, sum of medium-chain FA, n-6/n-3 FA ratio and PUFA/SFA ratio (P≤ 0.05). The treatment × time interaction was noticed for the following FAs: C18:1trans-10, C18:2cis-9,cis-12, C18:2cis-9,trans-11, and C18:2trans-10,cis-12 and sum of other FA, SFA, UFA, PUFA, n-6 FA,transC18:1, n-6/n-3 FA ratio, and PUFA/SFA ratio (P< 0.05).

Table 5 The effect on replacing alfalfa silage with paulownia leaves silage (60 g/kg) on bacteria, methanogens, methane (CH4) production and digestibility measured in rumen-cannulated cows (n = 4) (Exp. 2)

In vivo experiment (using commercial dairy cows)

The inclusion of PLS in the diet of dairy cows affected milk composition and ruminal methane concentration (Table 7). The PLS diets decreased protein and lactose yield (P< 0.05), however, it did not affect milk yield, energy corrected milk (ECM) and fat yield. The fat content was not affected by replacing alfalfa silage with PLS, but PLS diet decreased protein and lactose content (P< 0.05). Milk urea concentration was increased (P< 0.05) by PLS feeding. Methane concentration in the exhaled gas decreased (P< 0.01) about 14% in the PLS group compared to the CON group.

The PLS diet increased proportions of C15:0, C16:1, C18:2cis-9trans-11, C18:3cis-9cis-12cis-15 and C20:4 n-6 (P< 0.01) (Table 8). The proportions of C18:1trans-10, C18:1trans-11, and the sum oftrans-C18:1 was decreased (P< 0.05) by PLS feeding to dairy cows. The PLS diet decreased the total SFA (P< 0.05) proportion, but increased the total UFA (P< 0.05) and PUFA (P< 0.05) proportions. The PUFA/SFA ratio was higher (P< 0.05) in the PLS group than in the CON group, but the n-6/n-3 FA ratio was lower (P< 0.01) in the PLS group than in the CON group. Supplementation of PLS increased desaturase indices of C14:1, C16:1, and rumenic acid/(vaccenic acid + rumenic acid) (P< 0.05). Furthermore, the PLS affected the relative transcript abundances of five out of six analyzed genes. Replacing alfalfa silage with PLS resulted in increased mRNA expressions of all genes (P< 0.01) exceptACACAgene (Fig. 1).

Table 7 The effect on replacing alfalfa silage with paulownia leaves silage (60 g/kg) on milk production performance and methane concentration of commercial dairy cows (n = 16) (Exp. 3)

Table 8 The effect of replacing alfalfa silage with paulownia leaves silage (60 g/kg) on milk fatty acids (FA) composition (g/100 g FA) and desaturation (DI) of milk of dairy cows (n = 16) (Exp. 3)

Fig. 1 The effect of replacing diet with paulownia leaves silage (PLS: 60 g/kg) on expression of acetyl-CoA carboxylase 1 (ACACA), fatty acid synthase (FASN), stearoyl-CoA desaturase (SCD), fatty acid elongase (ELOVL), fatty acid desaturase 1 (FADS1), and lipoprotein lipase (LPL) genes in the milk of lactating cows (ns–not significant). (Exp. 3)

Discussion

The present study aimed at assessing the effect of replacing alfalfa silage with PLS in the dairy cow diet on ruminal fermentation processes such as VFA profile, methanogenesis, and biohydrogenation. We have demonstrated a beneficial effect of new forages PLS on ruminal fermentation processes and milk quality without a negative impact on milk production performance.

Dynamics of microorganism populations (bacteria, methanogens, protozoa) and ruminal fermentation characteristics (pH, ammonia, VFA, methanogenesis)

The dietary use of paulownia leaves has been studied for pigs, small ruminants, rabbits, and birds [6, 7] but no study has been conducted in high-producing dairy cows. The present study investigated ruminal fermentation, microbiota, milk production performance, and methane production variables. We recently published promising results of PLS on the ruminal environment [4]. The advantage of PLS over the fresh paulownia leaves resulted from high nutritive value (comparable to alfalfa) as well as high content of bioactive phenolic compounds (BAC), including phenolic acid (47 g/kg DM) and flavonoid (13 g/kg DM) contents [4]. The influence of PLS on ruminal fermentation characteristics was likely due to the high content of BAC, such as phenolic acids and flavonoids. The present in vitro and in vivo study is a direct continuation of previous work [4]. The intake of phenolic acid and flavonoid in case of the highest PLS dose (60 g/kg DM) was 65 g/d/cow and 18 g/d/cow, respectively. Dietary BAC can exert modulatory or even antimicrobial effects on ruminal microbial populations and thus the basic ruminal characteristics [24, 25]. Therefore, higher ruminal pH (P< 0.01) of PLS diets can be partly explained by BAC activity, which decreased (P< 0.01) the acetate-to-propionate ratio due to an increase of the numbers of lactate-consuming and propionate-producing bacteria preventing ruminal pH reduction [26]. In the current study, although population ofStreptococcus bovisincreased in the PLS groups (P< 0.01), its growth may not cause lactic acid accumulation in the ruminal fluid and thus pH reduction [27]. Increased pH was probably linked to a higher abundance ofMegasphaera elsdenii(P< 0.01) in the PLS diets as noted both in the present in vitro and in vivo experiments.M. elsdeniiutilizes lactic acid as the main source of energy and thus can protect other ruminal microorganisms from the negative impact of low pH resulting from lactic acid accumulation [28]. High pH positively affects the adhesion of fibrolytic bacteria to the feed particles and thus their degradation [29]. This phenomenon wasconfirmed in the present study which demonstrated an increase (P< 0.01) in the population of some bacteria species such asButyrivibrio fibrisolvensandPrevotellaspp. in response to PLS supplementation. Zhan et al. [30] showed that flavonoids (from 20 to 100 mg/kg body weight) increasedB. fibrisolvenspopulation and feed digestion through changes in the population of ruminal microorganisms. In the present experiment on cannulated cows, flavonoids were supplemented at the rate of 29 mg/kg of body weight, which corresponded to the lower range of Zhan et al. [30] study.

Besides changes in pH, PLS diet increased ammonia concentration (P< 0.01), both in vitro and in vivo. This observation can be explained by higher ruminal CP degradation compared to the control with non-synchronization of available energy and ammonia in the ruminal fluid [4, 31]. On the other hand, the CP total-tract degradability was reduced (P< 0.01) in the PLS cows. Such difference may result from different degradation of protein or distinct solubility of feed protein from alfalfa and paulownia leaves silages [4]. Changes in the observed ammonia concentration may also be associated with alterations in the number of protozoa and bacteria. Although lower protozoa can decrease ammonia concentration due to reduced recycling of engulfed bacterial protein by protozoa, decreased protozoal number may increase bacterial activity or growth in the rumen. In the current study, Ophryoscolecidae decreased (P< 0.01), whereas Isotrichidae protozoa significantly increased (P< 0.01) in the PLS diet. A previous study indicated that Ophryoscolecidae protozoa have higher bacterivory activity [32]. Thus, a decrease in the number of Ophryoscolecidae could increase selected bacterial populations. In the current study, numbers of proteolytic bacteria such asPrevotellaspp.(P< 0.01) andB. fibrisolvens(P< 0.05) in PLS diets increased considerably, both in vitro and in vivo.Prevotellaspp. in the rumen are the predominant proteolytic bacteria with diverse and broad range of peptidase activities and represent 20% to 60% of the bacterial abundance [33]. Therefore, it seems that the PLS diet increased proteolytic activity in the rumen, resulting in greater ammonia concentration in the rumen.

Inhibition of ruminal methanogenesis is usually associated with increase in propionate concentration due to competition for hydrogen. High levels of metabolic hydrogen in the rumen may shift fermentation towards propionate byPrevotellaspp.that increased (P< 0.01) in the PLS diets. RuminalPrevotellaspecies use different pathways for propionate production utilizing metabolic hydrogen via succinate or acrylate pathways for fermentation of sugars and lactate [34]. Moreover, Seradj et al. [35] suggested that elevated propionate can result from flavonoid supplementation (0.2 mg/g DM), which reduced methane production (P< 0.01) and population of methanogenic archaea (P< 0.01), and also increased (P< 0.01) population ofM. elsdeniiwhich has been confirmed in the present study.

In our previous in vitro study, paulownia leaves with high total polyphenols content (31 or 35 g/kg DM) lowered ruminal methane production and methanogen populations without affecting substrate degradability and volatile fatty acid concentrations [8]. Similar results were also noted in the current experiments. Besides, the addition of flavonoid-rich plants may impair the growth of some protozoa-associated methanogens [36]. It has been estimated that protozoa-associated methanogens are responsible for up to 37% of methane production [37]. In the present study, results obtained from rumencannulated dairy cows verified those from the in vitro experiment. This included reduced methane emissions (P< 0.01), increased pH (P< 0.01) and ammonia concentration (P< 0.01), elevated propionate and butyrate concentrations (P< 0.05), decreased methanogens (P< 0.01) and Ophryoscolecidae numbers (P< 0.01), and increased Isotrichidae population (P< 0.01). Differences in protozoa behavior in the ruminal fluid from rumen-cannulated cows (Exp. 2) may suggest their distinct reaction to the BAC of PLS origin.

Another effect of BAC x microorganisms’ interaction is the increase (P< 0.01) in the population ofFibrobacter succinogenes. The increase inF. succinogenesandB. fibrisolvenspopulations might be due to the presence of phenolic acids in PLS. Some phenolic acids such as hydroxycinnamic acid, syringic acid, and p-hydroxybenzoic acid were found to stimulate specific bacterial populations in a dose-dependent manner [38]. Though high concentrations of polyphenols are toxic to ruminal bacteria, a low concentration of polyphenol extract was also found to stimulate theF. succinogenesin the rumen [39]. According to Mitsumori et al. [34], elevated numbers ofF. succinogenesknown as non-H2producers may mitigate methane production through H2restriction. Zhang et al. [40] reported that increased propionate and less acetate negatively affect the abundance ofF. succinogenesandB. fibrisolvens. We did not observe such relations regardingF. succinogenesandB. fibrisolvenspopulation but we noticed decreased acetate (P< 0.01) and increased propionate concentration (P< 0.01) in ruminal fluid in in vitro and in vivo studies.

Dynamics in bacteria populations, rumen (FA) and milk parameters (basic composition, FA, transcript expression)

The present study revealed a tendency of reduction in milk production and a decrease in protein and lactose yield (P< 0.05) in response to PLS treatment. Khorsandi et al. [1] also reported decreased milk protein and lactose yield after using 120 g/kg DM of pomegranate byproduct silage, rich in phenolic compounds (49 g/kg DM). The inclusion of forages rich in phenolic compounds can decrease CP digestibility, as was observed in the current studies, and therefore may negatively affect CP ruminal and post-ruminal digestion. Moreover, in our previous study, the rapidly soluble fraction of protein degradation was higher in PLS than in alfalfa silage, which indicates an imbalance of available energy deficit accompanying increased ammonia concentration in the ruminal fluid [4]. Higher ammonia levels in the ruminal fluid may interact with milk parameters and result in increased urea content (P< 0.05) in milk when the diet was supplemented with PLE [41]. Considering this, it may be important to balance the dairy cow rations for easily fermentable nonstructural carbohydrates such as concentrates.

With regard to the fatty acids, an increase in C15:0 proportion (P< 0.01) in the ruminal fluid and milk suggests higher activity ofF. succinogenes.These microorganisms are responsible for synthesis of some FAs, including odd chain (mainly C15:0) and branched-chain FAs [42]. The PLS diet also modulated biohydrogenation of long chain FA such as C18:2cis-9,cis-12. According to Szczechowiak et al. [9], BAC derived from lingonberry (Vaccinium vitis-idaea) affected the biohydrogenation process with an increase in CLA isomer (C18:2cis-9,trans-11) and vaccenic acid (C18:1trans-11) as intermediates, which was also observed in the current study. The amount of C18:2trans-10,cis-12 in the ruminal fluid may be influenced by the activity ofM. elsdeniias shown by Kim et al. [43]. The amount of this conjugated linoleic acid isomer (CLA), and also α-linolenic acid is usually low in milk [44]. However, we did not observe similar changes in the milk.

In case of two FAs, markers of the biohydrogenation process such as C18:1trans-11 and C18:2cis-9trans-11, parallel changes were noticed in ruminal fluid and in milk. Synthesis of these two isomers is dependent on the activity ofB. fibrisolvens, one of the dominant group A bacteria in ruminal biohydrogenation responsible for transforming UFA into SFA [39, 45]. Despite an increase in abundance ofB. fibrisolvensin ruminal fluid of the PLS group, no elevation in SFA in favor of UFA was noticed in the current study. This may be due to the activity of Δ9-desaturase affecting proportion of the C18:1trans-11 [20]. On the other hand, BAC may reduce Δ9-desaturase activity, which may lead to the reduction of the C18:1trans-11 in the milk [46]. Another factor affecting concentration of the C18:1trans-11 is the activity of theB. proteoclasticus,a B bacteria group capable of efficiently biohydrogenating PUFA to SFA [47].B. proteoclasticusis recognized as an important bacterium that converts C18:1trans-11 to C18:0 and is one of the most sensitive ruminal bacteria to changes in dietary PUFA [48, 49]. Despite an increase in the population ofB. proteoclasticusin response to PLS, the proportion of C18:1trans-11 did not reduce. Elevated PUFA content was also not related to the reduction of this bacterial species. Considering the complexity of ruminal processes, some other factors cannot be ruled out.

Milk quality in terms of FA composition also depends on the activity of genes being expressed in the mammary gland. We investigated the transcript expression of six genes controlling FA metabolism. The relative transcript abundance of 5 genes increased in response to PLS. The fatty acid synthase (FASN), a multifunctional protein, is responsible for de novo biosynthesis of longchain PUFAs [50]. An increase in theFASNtranscript level was not accompanied by an elevation in long-chain FAs although the total and n-3 PUFA content increased, thus improving milk quality. Another factor affecting the elevation of the mentioned PUFAs may be the higher transcript level of theELOVL5gene. This gene encodes a protein responsible for extending long-chain PUFA [51, 52]. According to those authors, an increase in mRNA level ofELOVL5gene elevated both n-6 and n-3 PUFA content in milk. Lack of theELOVL5effect on the n-6 PUFA can be linked to increased transcript abundance ofFADS1,a gene limiting the n-3/n-6 rate of PUFA synthesis [53]. In our study, a major increase was observed for n-3 FAs (e.g., C18:3cis-9,cis-12,cis-15) whereas the proportion of only C20:4 n-6 increased without changing the total n-6 PUFA. We did not observe a limitation of the C18:2trans-10,cis-12 proportion although the increased abundance of theLPLmRNA was noticed. Unlike in the present study, Bryszak et al. [54] described the decreased level of theLPLtranscript accompanied by the reduced proportion of this FA. Milk quality is also influenced by MUFA content, especially C18:1cis-9, which is influenced by theSCDgene-regulating de novo synthesis of endogenous FAs [55]. Despite a higher transcript content ofSCDgene in the PLS group, only an increasing trend of C18:1cis-9 proportion was observed for PLS diet.

Conclusions

Under the conditions of the present study, PLS (60 g/kg DM) reduced methanogenesis, beneficially modulated ruminal fermentation and biohydrogenation processes without a negative impact on milk production performance of lactating dairy cows. Dietary PLS also improved milk FA profile including greater proportions of total UFA, PUFA, conjugated linoleic acid, and C18:1trans-11 along with reduction of n6/n3 ratio. The only negative impact of PLS was an increased ruminal ammonia concentration affecting milk urea content. For this reason, a better energy and protein synergy in PLS-containing diets would be required.

Abbreviations

aNDF: Ash free NDF; A/P ratio: Acetate/Propionate ratio; BAC: Biologically active compounds; BH: Biohydrogenation; CH4: Methane emission; CO2: Carbon dioxide emission; CLA: Conjugated linoleic acid; CP: Crude protein; D?: Desaturation ? at –n; DI: Desaturation index; DM: Dry matter; DMD: Dry matter digestibility; DMI: Dry matter intake; DNA: Deoxyribonucleic acids; EE: Ether extract; EI: Elongase index; ELOVL5: Fatty acid elongase 5; FA: Fatty acids; FADS1: Fatty acid desaturase 1; FASN: Fatty acid synthase; FISH: Fluorescence in situ hybridization; IVDMD: In vitro dry matter digestibility; LA:α-Linoleic acids; LCFA: Long chain fatty acids; LT: Longissimus thoracis; LNA:α-Linolenic acids; LPL: Lipoprotein lipase; MCFA: Medium chain fatty acids; mRNA: Messenger-RNA; MUFA: Monounsaturated fatty acids; NDFD: Neutral detergent fibre digestibility; NH3: Ammonia; OM: Organic matter; PBS: Phosphatebuffered saline; paulownia leaves silage; PLS60: Paulownia leaaves silege; Diets containing 60 g/kg DM of PLS; PUFA: Polyunsaturated fatty acids; qPCR: Quantitative PCR; RA: Rumenic acid; RT: Transcription reaction; Rusitec: Rumen simulation technique; SCD: Stearoyl-CoAdesaturase; SFA: Saturated fatty acids; SPE: Solid phase extraction; TI: Thrombogenicity; UFA: Unsaturated fatty acids; VA: Vaccenic acid; VFA: Volatile fatty acids.

Supplementary Information

The online version contains supplementary material available at https://doi.org/10. 1186/s40104-022-00745-9.

Additional file 1: Table S1.The sequences of primers specific to the analyzed bacteria species.

Acknowledgements

This study was supported by a grant from the National Science Center, Poland (Grant No. 2016/23/B/NZ9/03427). The publication was co-financed within the framework of the Polish Ministry of Science and Higher Education’s program: “Regional Initiative Excellence” in the years 2019–2022 (No. 005/RID/2018/19) “financing amount 12 000 000,00 PLN”. The authors are grateful to Magda Bryszak and Sylwia Inglot for their technical assistance.

Authors’ contributions

HH designed the study protocol, provided laboratory analysis, interpreted the data and wrote first version of the manuscript. DL reviewed and improved the manuscript. MSz improved the manuscript. AKP interpreted the data and reviewed the manuscript. MK provided laboratory analysis of methanogens. PK provided laboratory analysis of bacteria. MG provided statistical analyses. PK performed the laboratory analysis. SS provided laboratory analysis. DP provided laboratory analysis of rumen fluid. AC designed the study protocol, interpreted the data, wrote the manuscript, and finally reviewed the manuscript. All authors read and approved the final manuscript.

Availability of data and materials

Data available by the Author by request (adam.cieslak@up.poznan.pl).

Declarations

Ethics approval and consent to participate

The present study consisted of three consecutive experiments: one in vitro (Rusitec system) and two in vivo (rumen cannulated cows in lactation and milking cows under commercial production system). All experimental procedures were performed according to the National Ethical Commission for Animal Research guidelines (Ministry of Science and Higher Education, Poland). The study was approved by the Local Ethics Commission (decision no. 14/2019).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1Department of Animal Nutrition, Poznań University of Life Sciences, Wo?yńska 33, 60-637 Poznań, Poland.2Department of Genetics and Animal Breeding, Poznań University of Life Sciences, Wolynska 33, Poznań 60-637, Poland.3Department of Animal Nutrition, West Bengal University of Animal and Fishery Sciences, 37 K. B. Sarani, Kolkata, India.4Institute of Genetics and Animal Biotechnology, Polish Academy of Sciences, Post?pu 36A, 05-552, Magdalenka, Warsaw, Poland.5Department of Animal Physiology, Biochemistry and Biostructure, Poznan University of Life Sciences, Wolynska 35, Poznan 60-637, Poland.6Department of Pharmaceutical Biology and Botanic Garden of Medicinal Plants, Wroc?aw Medical University, Wroc?aw 50-556, Poland.7Institute of Animal Physiology, Centre of Biosciences of Slovak Academy of Sciences, ?oltésovej 4-6, 040-01 Ko?ice, Slovak Republic.8Department of Animal Nutrition, Faculty of Veterinary Medicine and Animal Science, Poznań University of Life Sciences, Poznań 60-637, Poland.

Received: 7 March 2022 Accepted: 9 June 2022