Effects of Brewers’ Spent Grain Protein Hydrolysates on Gas Production, Ruminal Fermentation Characteristics, Microbial Protein Synthesis and Microbial Community in an Artificial Rumen Fed a High Grain Diet

DOI: https://doi.org/10.21203/rs.3.rs-32671/v1

Abstract

Background: Brewers’ spent grain (BSG) contains 20 ~ 29% protein with high amount of glutamine, proline and hydrophobic and non-polar amino acid residues, making it an ideal material for producing value-added products like bioactive peptides. Protein were extracted from BSG, hydrolyzed with 1% alcalase and flavourzyme, and the generated protein hydrolysates (AlcH and FlaH) showed antioxidant activities. This study was conducted to evaluate effects of AlcH and FlaH on gas production, fermentation characteristics, nutrient disappearance, microbial protein synthesis and microbial community using an artificial rumen system (RUSITEC) fed a high-grain diet.

Results: Supplementation of FlaH decreased (P < 0.01) disappearances of dry matter (DM), organic matter (OM), crude protein (CP) and starch, without affecting fibre disappearances; while AlcH had no effect on nutrient disappearance. Neither AlcH nor FlaH affected gas production and VFA profiles, except they enhanced (P < 0.01) NH3-N but decreased (P < 0.01) H2 production. Supplementation of FlaH decreased (P < 0.01) percentage of CH4 in total gas and dCH4 in dissolved gas. Addition of monensin reduced (P < 0.01) nutrient disappearances, improved fermentation efficiency and reduced CH4 and H2 emission. Total microbial nitrogen production decreased (P < 0.05) but proportion of feed particle associated (FPA) bacteria increased with FlaH and monensin. Numbers of OTUs or Shannon diversity indices of FPA microbial community were unaffected by AlcH and FlaH; whereas both indices were reduced (P < 0.05) by monensin. Taxonomic analysis revealed no effect of AlcH and FlaH on the relative abundance (RA) of bacteria at phylum level; monensin reduced (P < 0.05) the RA of Firmicutes and Bacteroidetes and enhanced Proteobacteria. The FlaH enhanced (P < 0.05) the RA of genus Prevotella, reduced Selenomonas, Shuttleworthia, Bifidobacterium and Dialister as compared to control; monensin reduced (P < 0.05) RA of genus Prevotella but enhaced Succinivibrio.

Conclusions: Inclusion of FlaH in high-grain diet may potentially protect CP and starch from ruminally degradation, without adversely affecting fibre degradation and VFA profiles. The FlaH also showed promising effects on reducing CH4 production by suppressing H2 generation. Protein enzymatic hydrolysates from BSG using flavourzyme showed potential application to high value-added bio-products.

Background

Brewers’ spent grain (BSG) is the most abundant by-products generated in brewing industry, accounting for 85% of the total by-products generated [1]. Currently, the main application of BSG is as animal feed, especially for ruminants, as it contains up to 29% protein and up to 60% fiber on dry matter (DM) basis [2]. Inclusion of low-level wet BSG to conventional finishing diet can maintain growth performance and meat quality of finishing beef cattle [3]. However, there are some potential disadvantages when applying BSG in high amount in ruminant diets, such as reducing feed intake and increasing incidence of liver abscess [4], increasing the releasing of enteric pathogenic bacteria E. coli O157:H7 [5, 6], increasing secretion of N, P and S to the environment [7]. Furthermore, the utilization of BSG in animal production were limited due to its high moisture and fermentable sugar contents, which make it very unstable and easily spoiled, thus not only limited its use generally to animal operations near breweries due to high energy costs of transportation, but also made the storage of wet BSG very challenging because of microbial activity [1]. Although drying is considered the most effective method for lowering storage and transportation costs, drying itself is an energetically expensive process [1, 8]. Therefore, due to the high yield but limited usage of BSG, there is an urgent need for developing new strategies to increase the value of BSG to turn it waste into treasure.

Nowadays, there are increasing interests in the study of by-products from the food processing industry and animal production for extraction of value-added bioactive compounds, which can possibly be used as natural antimicrobials or antioxidants in addition to nutritional use [9]. Barley is the primary grain source used in the brewing industry, and the BSG derived from barley has rich contents of phenolic compounds and hordein and glutelin proteins. Emerging evidence indicate that the phenolic compounds in diet exhibit the ability of anti-inflammatory and antioxidant activities [8, 10]. Meanwhile, the unique structural features of barley proteins, high amount of glutamine and hydrophobic and non-polar residues, offer the possibility of producing bioactive peptides with strong antioxidant and antimicrobial properties upon proteolytic degradation using commercially available proteinase or alkali solutions [11, 12]. Bioactive peptides have a broad range of functions: antibacterial, anti-inflammation, antioxidant and immune function enhancing. Both glutelin and hordein protein hydrolysates have been proved having antioxidant activities [11, 13]. Protein hydrolysates produced from BSG not only exert antioxidant activity under conventional chemical assay, but also exert antioxidant and immunomodulatory effects to cell lines (U937, Jurkat T cell, Caco-2 and HepG2) under oxidative stress in vitro [12, 14]. Consequently, the BSG, a low cost and available in large amount by-product, is an ideal material that can be used for producing bioactive compounds. Furthermore, our previous study proved that the BSG residue resulting from protease aided protein removal can be potentially used as a viable fibre source for ruminant feeding [15], which will make full use of the BSG residue and add more value.

Antibiotics, widely used in north American beef industry, is questioned by the consumers because of potential antibiotic residue accumulation in animal products, and will be banned as growth promotant due to the increasing risk of antimicrobial resistance. A lot of efforts were tried to develop alternatives to in-feed antibiotics using natural source materials. Thus, we hypothesis that inclusion of BSG protein hydrolysates may exert antimicrobial activity in rumen environment, and protein hydrolysates of BSG could be used a functional feed additive to substitute antibiotics in ruminant diets. The objective of this study was to determine if BSG protein hydrolysates prepared with two different proteolytic enzymes (alcalase and flavourzyme) had the ability to modulate rumen fermentation, reduce methane production and alter microbial community in rumen simulation technique (RUSITEC) fed a high-grain diet.

Methods

Preparation of BSG protein hydrolysates

The BSG protein hydrolysates used in the current study were prepared following alkali extraction, acid precipitation and enzyme hydrolysis as previously reported [11, 13]. In brief, BSG was milled through a 0.5-mm screen, solubilized in 0.1 M NaOH at 20% (w/v) rate, followed by continuous stirring at 350 rpm at 50oC for 2 h for protein extraction. Then the supernatant, collected by centrifuging at 8,000 × g for 15 min at 20oC, was passed through an ultrafiltration membrane to obtain the protein fraction > 1 kDa, which was furtherly precipitated by acid (adjusting pH to 3.5) and centrifuged at 8,000 × g for 15 min at 20oC, freeze dried. The obtained protein was dispersed in deionized water to reach a 5% (w/v) solution, hydrolyzed by 1% alcalase and flavourzyme at their optimum pH and temperatures (pH 8.0, 55oC and pH 6.6, 50oC, respectively). At the end of hydrolysis, all the hydrolysate solutions were adjusted to pH 7.0, heated at 95oC for 5 min to inactive the enzyme, and centrifuged at 8,000 × g for 30 min to separate solubilized peptides and amino acids from the non-soluble substrates. The obtained hydrolysates with alcalase and flavourzyme were referred to AlcH and FlaH, respectively. The protein content and sample degree of protein hydrolysis were determined as described [11, 13].

Experimental design and treatments

The experiment was a completely randomized block design with four treatments assigned to sixteen fermentation vessels in two units of RUSITEC apparatus. Treatments were control (no hydrolysates, no antibiotics), 1% AlcH, 1% FlaH (diet dry matter; DM), and antibiotics (33 mg monensin + 11 mg tylosin/kg diet DM; positive control). The substrate was a high concentrate diet that contained 10% barley silage, 87% dry-rolled barley grain, and 3% vitamin and mineral supplement (DM basis; Table 1). The experiment lasted 15 d in duration, including 8 d for adaptation and followed by 7 d for sampling and data collection. The basal diet was prepared as total mixed ration, and grounded through a 4-mm sieve (Arthur Thomas Co., Philadelphia, PA, USA). Approximately 10 g (DM) of diet was weighed into nylon bags (10 × 20 cm; pore size of 50 µm, Ankom Technology Corp., Macedon, NY, USA). The protein hydrolysates and antibiotics were added to bags at the desired amount and manually mixed.

Table 1

Antioxidant activity of brewers’ spent grain (BSG) protein hydrolysates generated using alcalase and flavourzyme

Item

BSG

AlcH1

FlaH2

SEM

P value

BSG protein extraction, % of dry matter

31.7

70.0

70.0

Solubility, %

37.5b

89.1a

86.5a

2.35

0.01

Degree of hydrolysis, %

34.7a

11.9b

5.37

0.01

1AlcH = Alcalase hydrolysates;
2FlaH = Flavourzyme hydrolysates

Inoculum donor

Three ruminally fistulated Angus cross cows (average, 768 ± 95.1 kg BW) fed a high grain diet containing 8.2% barley silage, 89.2% dry rolled barley grain, and 2.6% vitamin and mineral supplement (DM basis) were used as rumen inoculum donor. Two hours after morning feeding, solid and liquid contents were collected from four locations within the rumen of each cow via rumen cannula. The contents were immediately filtered through four layers of cheesecloth, pooled over in equal amount (4 L per cow), and recorded for pH before adding into fermenters. The experimental protocols were reviewed and approved by the Lethbridge Research and Development Centre Animal Care Committee, and the cows were handled in accordance with the guidelines of the Canadian Council on Animal Care [16].

Experimental procedure

The RUSITEC procedure were carried out using two units of RUSITEC apparatus, equipped with eight 920-mL anaerobic fermenters of each unit, as described previously [17, 18]. In brief, to initiate the fermentation, each fermenter was filled with 200 mL of McDougall’s buffer [19] and 700 mL of prepared rumen inoculum, and then one bag containing 20 g of prepared solid rumen digesta and one bag containing 10 g of experimental diet (DM basis) were placed in each fermenter at 0900 h on d 0. Fermenters were placed in a circulating water bath at 39 °C for incubation. After 24 h of incubation, the bag containing solid rumen digesta was replaced with a bag containing experimental diet. Thereafter, a bag was replaced daily, so that each bag remained in each fermenter for 48 h. During the daily feed-bag exchange, the fermenter was flushed with N2 gas to maintain anaerobic condition in the fermenters. The artificial saliva was continuously infused into fermenters using a peristaltic pump (Model ISM 932D, Ismatec, Index Health and science GmbH, Wertheim, Germany) at a dilution rate of 2.9%/h. Effluent and fermentation gasses of each fermenter were collected, respectively, into a 2 L Erlenmeyer flask and a reusable 2 L gas-tight collection bags (CurityR; Conviden Ltd., Mansfield, MA, USA), and they were recorded at the feed bag exchange.

Nutrient disappearance

Disappearance of DM, OM, CP, acid detergent fiber (ADF), neutral detergent fiber (NDF) and starch was measured with 48-h incubated feed bags from d 9 to 13 of the sampling period. Bags were withdrawn from fermenters and washed manually under running cold water until the water was clear, and dried at 55 °C for 48 h (AOAC, 2005; method 930.15) to determine DM disappearance. Thereafter, the bag residues were pooled over 5 d by fermenter, and ground through a 1-mm sieve, for DM, OM, NDF and ADF analysis. A portion of ground sample was furtherly ball ground using a ball mill (Mixer Mill MM2000; Retsch, Haan, Germany) for total N and starch analyze. Disappearance of DM, OM, CP, NDF, ADF and starch was calculated as the differences between the amount of input and the amount remaining in residues.

Gas production and dissolved gases

Fermentation gasses were measured every 24 h using a gas meter (Model DM3A, Alexander-Wright, London, England, UK) from d 9 to 15 of the sampling period. The gas (20 mL) was sampled from each bag once daily for determining gas profiles. On d 14 and 15, fermentation liquids were sampled in duplicates (35 mL each) from each fermenter for measuring concentration of dissolved H2 (dH2) and dissolved CH4 (dCH4). The sampling procedure was carried out using two syringes as reported [20]. In brief, a 50 mL syringe that filled with 35 mL of fermentation liquids connected to a 20 mL syringe that filled with 5 mL of N2 through a T tube with a valve. The N2 was injected from the small syringe into the large syringe and the apparatus was vigorously shaken by hand. The entire gas phase was then transferred from the large syringe into the small one to determine the gas volume. Finally, the small syringe was removed from T tube and 6 mL of gas was sampled for both gas and dissolved gas analyses using a Varian 4900 Gas Chromatograph (Agilent Technologies Canada Inc., Mississauga, Ontario, Canada).

Fermentation parameters

The pH of fermentation fluid of each fermenter was measured daily using a pH meter (Orion model 260A, Fisher Scientific, Toronto, ON, Canada) at the time of feed-bag exchange. From d 9 to 13 of sampling period, effluent (5 mL) was sampled and preserved with 1 mL of 25% metaphosphoric acid for VFA analysis; another 5 mL of effluent were preserved with 1 mL of H2SO4 (1% vol/vol) for NH3-N analysis. All samples were well mixed and frozen at -20 °C until analyzed. The production (mmol/d) of total VFA, individual VFA and NH3-N were calibrated based on daily effluent volume.

Microbial protein synthesis

Bacteria in the fermenters were labeled using 15N. From d 7 to 15, 0.3 g/L (NH4)2SO4 in McDougall’s buffer was replaced with 0.3 g/L 15N-enriched (NH4)2SO4 (Sigma Chemical Co., St. Louis, MO, USA; minimum 15N enrichment 10.01 atom%). From d 9 to 15, the effluents were preserved by adding 3 mL of a sodium azide solution (20%; wt/vol) to each effluent flask. On d 14 and 15, daily volume of effluent from each fermenter was recorded and 35 mL of effluent was sampled and centrifuged at 20,000 × g, 4 °C for 30 min to isolate liquid-associated bacteria (LAB). The obtained pellets were washed with deionized water and centrifuged three times (20,000 × g, 30 min, 4 °C), suspended in distilled water, freezing and lyophilisation for determination of N and 15N.

Feed particle-associated (FPA) and feed particle-bound (FPB) bacterial fractions were measured from 48-h feed residues. After 48-h of incubation, feed bags were squeezed to expel excess liquids, placed individually in a plastic bag with 20 mL of McDougall’s buffer, and processed for 1 min using a Stomacher 400 Laboratory Blender (Seward Medical Ltd, London, UK). Then, the liquid from each feed bag was squeezed into a 50 mL centrifuge tube, the feed residues were washed twice with 10 mL of McDougall’s buffer in each wash, and the washed buffer was pooled over with the squeezed liquid to obtain the FPA bacterial fraction. The washed feed residues were considered as the FPB bacterial fraction. The obtained FPA samples were centrifuged at 500 × g, 4 °C for 10 min to remove large feed particles, the supernatant was centrifuged (20,000 × g, 30 min, 4 °C) to isolate bacterial pellets, which was furtherly processed as aforementioned. Washed feed residues (FPB fraction) were dried at 55 °C for 48 h, weighed for DM determination, and ball grinded (MM400; Retsch Inc., Newtown, PA, USA) for N and 15N analysis. Total microbial protein synthesis was estimated as the sum of LAB, FPA and FPB.

Microbial community

Microbial community of FPA samples were assessed through high-throughput sequencing. Total DNA was extracted from FPA samples (30 mg) using a QIAamp Fast DNA stool mini kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. The quality and quantity of extracted DNA was measured using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), then the extracted DNA was stored at -20 °C until sequencing. The V4 hypervariable region of the archaeal and bacterial 16S rRNA gene was amplified using the modified 515-F and 806-R primers, with the PCR conditions and sequencing steps carried out as previously described by [21]. Briefly, the 16S rRNA gene amplicons were generated using a two-step PCR, and then the amplicons were subjected to Illumina paired-end library preparation, cluster generation, and sequenced on an Illumina MiSeq instrument (Illumina, Inc., San Diego, CA, USA).

The obtained 16S rRNA gene sequencing raw data were processed using QIME2 [22] and the R-package DADA2 (Version 1.4) denoise method as described. In brief, after removing primer sequences and truncating both the forward and reverse reads at 225 bp, quality control was done for the reads using the QIME2, with chimeric sequences identified and removed. Then, the richness (number of OTUs) and diversity (Shannon index) were calculated and non-metric multidimensional scaling (NMDS) was performed based on Bray-Curtis similarity distances using R packages vegan (Version 2.4.4; [23]) and phyloseq (Version 1.20.0; [24]). Fold change of ruminal bacterial at genus level with a threshold of 5% was analyzed using R-package Deseq2 [25].

Chemical analysis

The chemical analyze of the feeds and feed residues were conducted in duplicate, and analysis was repeated when the CV for the replicate analysis was more than 5%. Analytical DM was measured by oven drying at 135 °C for 2 h (AOAC, 2005; method 930.15) [26] and ash content was analyzed by combustion samples at 550 °C for 5 h (AOAC, 2005; method 942.05). The NDF and ADF were determined on a VELP Fiber Digestion System (VELP Scientifica, Burlington, Ontario, Canada) using the method of Van Soest et al. [27] and AOAC (2005; method 973.18), respectively. Total N of feed and residue samples and 15N of LAB, FPA and FPB samples were analyzed using combustion analyzer (NA 2100, Carlo Erba Instruments, Milan, Italy), and CP was calculated as total N × 6.25. Starch was determined by enzymatic hydrolysis of α-linked glucose polymers as reported previously [28]. Concentration of VFA and NH3-N in the effluent was determined using a gas chromatograph (model 5890, Hewlett-Packard Lab, Palo Alto, CA, USA).

Statistical analysis

Data were analyzed according to a completely randomized block design using the MIXED procedure of SAS (Version 16.0.0, SAS Inst. Inc., Cary, NC, USA), with treatment considered as fixed effect, day of sampling as repeated measures, and the fermenter and RUSITEC apparatus as random effects. For the repeated measures, various covariance structures were tested with the final structure chosen based on the lowest Akaike’s information criteria value. Results are reported as least squares means, which were compared using the Tukey correction for multiple comparisons. Significance among treatments was declared at P ≤ 0.05 and a trend at 0.05 < P ≤ 0.10 unless otherwise stated.

Results & Discussion

BSG protein hydrolysates of AlcH and FlaH

The protein content of BSG protein extract was 70%, which was much higher than that of the original BSG substrate (Table 1). The present value was comparable to a recent study that combination of enzyme and ultrasonication, resulting in 69.8% of BSG protein content [29]; whereas, it was greater than pre-treated BSG with carbohydrases followed by direct hydrolysis using proteolytic enzymes (63.1%) [30]. The Degree of hydrolysis (DH) of the protein extracted from BSG was greater (P < 0.01) when it was treated with alcalase (34.7%) than with flavourzyme (11.9%) and both protein hydrolysates showed good solubility in water. The DH expresses the number of peptide bonds cleaved as a percentage of total number of peptide bonds in the substrate; the higher DH can be expected to have lower molecular weight peptides and vice versa. Bamdad et al. [11] showed good antioxidant activities of barley hordein and glutelin enzymatic hydrolysates. Especially the flavourzyme and alcalase hydrolysates had superior DPPH (1,1-diphenyl-2-picryl hydrazyl) free radical scavenging activity, Fe2+-chelating ability and superoxide radical scavenging capacity. The peptides high in hydrophobic and non-polar amino acid residues were most likely responsible for the antioxidant effects [11, 13]. As the hydrophobic and non-polar amino acids are mostly remained in BSG after brewing processing [11], it was hypothesized that the BSG hydrolysates may have beneficial effects in rumen feeds based on their antioxidant activities.

Effects of AlcH and FlaH on nutrient disappearance and fermentation characteristics

Supplementation of FlaH reduced (P < 0.01) the disappearance of DM, OM, CP and starch as compared to the control and AlcH, which had no difference in the nutrient disappearances (Table 2). The reduction of DM disappearance by FlaH consists with the results obtained in our previous study using batch culture technique. We found that the DM disappearance of barley grain linearly decreased from 76.4, 70.7 to 63.8% as the inclusion rate of FlaH increased from 0, 5 to 10% (unpublished data). The reduction in the disappearances of DM and OM by FlaH was primarily resulted from the decreased of CP and starch. It appeared that the effect of AlcH and FlaH on the nutrient disappearances was not associated with their antioxidant activities measured in the present study. In fact, if we admit that the increased the DPPH scavenging or Fe-chelating activity adversely affected the nutrient disappearance as observed with FlaH, the nutrient disappearance with AlcH would be less than that with FlaH because of the greater DPPH scavenging or Fe-chelating activity of AlcH than FlaH. These results suggest that the hydrolysate with chemically determined antioxidant activity may have different response under rumen fermentation condition. Furthermore, the fermentation under current Rusitec condition would be less in oxidative stress compared with in vivo where the oxygen can frequently enter the rumen with feeds and rumination. We also suggest that the hydrolysates may have more like antimicrobial than antioxidant activity in the fermentation. The different bioactivities of AlcH and FlaH in rumen fermentation would be expected, as the activity of protein hydrolysates depends on the specificity of the protease used, hydrolysis conditions and degree of hydrolysis [13]. The flavourzyme is an endo- and exopeptidase enzyme mixture, and mainly produce small peptides and free amino acids, while the alcalase is an endo-protease and mainly generate small- and medium-sized peptides [13]. Furtherly, about 6% of protein in the total protein input was hydrolysate origin as a result of adding FlaH, the decreased CP disappearance with FlaH could be due to the resistance of hydrolysate to microbial degradation. The resistance is related to the amino acid sequences and structures, especially peptides with Gly-Gly, Pro-X or X-Pro residues at the N-terminus, or with modification of N-terminal amino groups [31, 32]. Thus, it is predicted that FlaH generated by flavourzyme is more resistant to rumen microbial digestion, and possess strong bioactivity in modulating rumen fermentation, including reduce the activity of proteolytic and starch utilization microbes. The disappearance of NDF and ADF were not affected by supplementation of either AlcH or FlaH, suggesting that both protein hydrolysates had limited effect on the activity of fibrolytic microbes.

Table 2

Ingredient and chemical composition of experimental diets

Item

Contents

Ingredient, %

 

Barley silage1

10.0

Barley grain2, ground

87.0

Supplement3

3.0

Chemical composition, % of DM

       

DM

91.6

OM

95.2

CP

11.8

NDF

20.0

ADF

8.0

Starch

50.3

1 Composition (DM basis): 31.8% DM, 94.1% OM, 42.9% NDF, 26.7% ADF, 16.4% starch, and 9.6% CP.
2 Composition (DM basis): 90.2% DM, 98.4% OM, 14.9% NDF, 4.1% ADF, 55.9% starch and 12.4% CP.
3 Supplied per kilogram of dietary DM: 565 g barley grain, 100 g canola meal, 250 g calcium carbonate, 25 g molasses, 30 g salt, 20 g urea, 0.66 g vitamin E 500 and 10 g premix. The premix in the supplement contained per kilogram of dietary DM: 15 mg of Cu, 65 of mg Zn, 28 mg of Mn, 0.7 mg of I, 0.2 mg of Co, 0.3 mg of Se, 6,000 IU of vitamin A, 600 IU of vitamin D, and 47 IU of vitamin E.

The pH in the fermenters remained constant throughout the whole experimental period, with slightly higher (P < 0.01) pH with Ant than other treatments (average, 5.76; Table 3). Production of VFA (mmol/d) and individual VFA molar proportion as well as acetate to propionate ratio were not affected with adding AlcH or FlaH except for lower proportion of butyrate with FlaH than control and AlcH. The similar VFA production is consistent with biologically minor difference in the OM disappearance despite of statistical difference between FlaH and control or AlcH. However, AlcH and FlaH had greater (P < 0.01) NH3-N production than control. The greater NH3-N production with FlaH is not expected because of decreased CP degradability by FlaH. It was probably because they provided certain amount of peptides and free amino acid to fulfil partial nitrogen requirements of ruminal microbes. It was reported that ruminal microbes preferred to use peptides or amino acids as a source of nitrogen or as a source of energy, thus lead to accumulation of NH3-N [33]. We speculate that the supplementation of AlcH and FlaH can promote the deaminative activity of high ammonia producing bacteria. Ammonia is mainly produced by the low activity species, but proliferation of the high ammonia producing bacteria will be a problem if certain diets were fed [31]. It was reported that a variety of ruminal bacteria produced ammonia from protein hydrolysates, with strains of Bacteroides ruminicola, Megasphaera elsdenii, and Selenomonas ruminantium being the most active [34].

Table 3

Effect of brewers’ spent grain (BSG) protein hydrolysates or antibiotics on nutrient disappearances and fermentation characteristics in RUSITEC

Item

Treatments1

SEM

P value

Con

AlcH

FlaH

Ant

Nutrient disappearance, %

DM

76.5a

76.1ab

75.3b

72.7c

0.90

0.01

OM

77.7a

77.6a

76.5b

74.1c

0.87

0.01

CP

75.1a

75.4a

73.4b

69.7c

1.12

0.01

NDF

40.3a

38.7a

38.9a

31.8b

1.32

0.01

ADF

29.2a

26.8a

27.2a

21.5b

1.33

0.01

Starch

90.7a

90.8a

89.5b

87.8c

0.92

0.01

Fermentation characteristics

           

pH

5.74b

5.78b

5.77b

5.85a

0.02

0.01

NH3-N, mmol/d

3.66b

4.05a

4.19a

2.57c

0.12

0.01

Total VFA, mmol/d

54.80a

54.51a

53.63a

48.88b

0.86

0.01

Individual VFA, % of total VFA

         

Acetate (A)

30.97a

31.12a

31.67a

29.77b

0.38

0.03

Propionate (P)

39.33b

39.39b

39.42b

42.98a

0.39

0.01

Butyrate

20.37a

19.98a

19.20b

16.76c

0.45

0.01

BCVFA2

2.21a

2.19a

2.19a

1.56b

0.03

0.01

Valerate

4.25b

4.33b

4.78b

7.12a

0.41

0.01

Caproate

2.64a

2.77a

2.80a

1.57b

0.17

0.01

A:P

0.79a

0.79a

0.80a

0.69b

0.01

0.01

a, b, c Least square means within a row with different superscripts differ (P < 0.05).
1 Con = Control, no antioxidant peptide and no antibiotics; AlcH = Alcalase hydrolysates, 10 mg per gram of TMR (DM basis); FlaH = Flavourzyme hydrolysates, 10 mg per gram of TMR (DM basis); Ant = antibiotics, 0.8 mg monensin + 1 mg tylosin per gram of TMR (DM basis).
2 Branched-chain volatile fatty acids (isobutyrate + isovalerate).

The supplementation of antibiotics (P < 0.01) decreased disappearance of DM, OM, CP, NDF, ADF and starch compared with other treatments. Therefore, lower (P < 0.01) production of total VFA and higher (P < 0.01) fermenter pH was observed in Ant treatment. Furthermore, Ant altered the individual VFA proportion as compared to the other treatments, with lower (P < 0.01) molar proportion of acetate, butyrate, BCVFA and caproate, and higher (P < 0.01) molar proportion of propionate and valerate, thus led to lower (P < 0.01) acetate to propionate ratio. Our results confirm the general known mode of action of monensin in the rumen [35]. Additionally, supplementation of monensin and tylosin reduced (P < 0.01) the production of ammonia, which was in alignment with previous reports that ionophores could decrease ammonia production by suppressing high ammonia producing microbial population, as well as the peptidelytic and deaminative activity of the bacteria that grow in the presence of ionophores [35]. Although the FlaH had similar effects on reducing CP and starch disappearance as Ant did, fibre disappearance and fermentation pattern were not altered. It is likely that mode of action in the rumen is different between monensin and FlaH. The addition of FlaH protected protein and starch from rumen fermentation, which can let more protein enter the small intestine for digestion and alleviate risk of rumen acidosis.

Effects of AlcH and FlaH on gas and dissolved gas production

H2 was produced during the oxidizing of reduced electron carries (such as NADH and ferredoxin) by membrane-bound hydrogenases of some rumen microbes, which plays an important role in maintaining the oxidation-reduction homeostasis to guarantee rumen anaerobic fermentation [36]. However, the activity of membrane-bound hydrogenase will be quickly inhibited by high H2 pressure, while the methanogenic Archaea can efficiently use H2 to reduce CO2 to keep a low ruminal H2 pressure, but at the price of producing CH4 [37]. Therefore, either reducing the production of H2 or finding alternative sinks for H2 are theoretically feasible to reduce CH4 emission. It is now clear that ionophore antibiotics, like monensin, can decrease CH4 production via inhibiting growth of H2 producing bacteria, without causing side effects to succinate- and propionate-producing bacteria [35]. In the current study, fermentation liquid and gas were analyzed for measuring dissolved gas (dH2 and dCH4) and gas (H2 and CH4) production, respectively. Production of dissolved gas and total gas (L/d) as well as the production of dH2 were not affected by treatments, whereas the production of dCH4 (% of dGas, mg/d or mg/g digested DM), and CH4 (% of gas) were lower (P < 0.01) with FlaH than with control and AlcH (Table 4). Supplementation of either AlcH or FlaH considerably decreased (P < 0.01) the production of H2, expressed as µg/d or µg/g digested DM. Adding monensin and tylosin consistently decreased (P < 0.05) the production of CH4, H2 and dCH4 compared with other treatments. The reduced production of H2 by adding FlaH (-49%) is of interest, and it consistent with the decrease in dCH4 production and in the proportion of CH4 in total gas. However, the molar proportion of propionate was unchanged by adding FlaH despite of the decreased in production of CH4 and H2. The decrease in OM disappearance with FlaH may be unable to produce sufficient H2 to increase propionate production. These results suggest that the FlaH may especially target at methanogenesis. The decreased production of H2 with AlcH without altering the production of CH4 and propionate is not clear. Anyway, the OM disappearance was not affected by AlcH. The present results confirmed the effect of monensin on reducing CH4 and H2 production and increasing propionate production [35]. Although the FlaH showed reduction of CH4 and H2 as well as the reduction in nutrient disappearance, the magnitude of the reduction by FlaH was much less than by monensin, and in particularly, the FlaH did not increase propionate production and alter fermentation pattern. Therefore, it suggests lower activity of FlaH than monensin and different mode of action between the two additives. Future work also needs to be carried out to determine relationships between their physicochemical and techno-functional properties.

Table 4

Effects of brewers’ spent grain (BSG) protein hydrolysates or antibiotics on gas and dissolved gas production in RUSITEC

Item

Treatments1

SEM

P value

Con

AlcH

FlaH

Ant

Dissolved Gas2

           

Total dGas, L/d

0.95

0.96

0.96

0.94

0.02

0.17

dCH4, % of dGas

0.72a

0.68a

0.63b

0.53c

0.04

0.01

dCH4, mg/d

4.42a

4.21a

3.86b

3.17c

0.23

0.01

dCH4, mg/g DM digested

0.59a

0.58a

0.53b

0.45c

0.03

0.01

dH2, µg/d

0.48

0.63

0.58

0.54

0.20

0.38

dH2, µg/g DM digested

0.06

0.08

0.08

0.07

0.01

0.39

Gas Production3

           

Total Gas, L/d

1.68

1.62

1.59

1.58

0.20

0.87

CH4, % of Gas

2.47a

2.55a

2.11b

1.32c

0.12

0.01

CH4, mg/d

27.49a

27.15a

24.37a

14.92b

5.38

0.01

CH4, mg/g DM digested

3.80a

3.77a

3.49a

2.19b

0.72

0.01

H2, µg/d

105.1a

65.3b

54.3b

18.0c

11.01

0.01

H2, µg/g DM digested

14.6a

9.0b

7.6b

2.6c

1.58

0.01

a, b, c Least square means within a row with different superscripts differ (P < 0.05);
1 Con = Control, no antioxidant peptide and no antibiotics; AlcH = Alcalase hydrolysates, 10 mg per gram of TMR (DM basis); FlaH = Flavourzyme hydrolysates, 10 mg per gram of TMR (DM basis); Ant = antibiotics, 0.8 mg monensin + 1 mg tylosin per gram of TMR (DM basis);
2 Samples from d 14 to 15;
3 Samples from d 9 to 13.

Effects of AlcH and FlaH on microbial protein synthesis and microbial community

Production of total microbial N (LAB + FPB + FPA) and LAB was less (P < 0.05) with FlaH than control and AlcH, whereas the production of FPA was greater (P < 0.01) with AlcH and FlaH without difference in FPB among treatments (Table 5). Supplementation of monensin and tylosin produced greater (P < 0.01) FPA but less LAB and total microbial N compared to control and AlcH without differing with FlaH except for LAB, which was less (P < 0.01) with Ant than FlaH. Microbial N production efficiency did not differ among treatments. Although the total microbial N production reduced with adding FlaH or Ant, proportion of attached microbial biomass (FPA + FPB) in the total biomass was greater for FlaH (30.1%) or Ant (34.4%) than control (24.4%). Our results suggested that adding FlaH or Ant benefits to microbial colonization to feed particle. As the supplementation of AlcH and FlaH to basal diet increased the microbial population in the FPA fraction, which was interesting because the microbial attachment is essential to the feed digestion. Thus, the FPA samples were selected for high-throughput sequencing to assess the microbial community. Neither numbers of OTUs (Fig. 1A) nor Shannon diversity indices (Fig. 1B) of FPA microbial community were affected by supplementing AlcH and FlaH. Whereas, both indices were reduced (P < 0.05) by Ant. Similarly, the results of NMDS analysis indicated that there was no clustering of the FPA microbial community in control, FlaH and AlcH treatments, but the Ant treatment had more clustered community structure (Fig. 2). The study on the effects of BSG protein hydrolysates on ruminal microbiome is lacking. A study applied a blend of plant essential oils with antioxidant activity demonstrated no effect on ruminal microbiome [38]. Meanwhile, our results that monensin greatly reduced the diversity of FPA microbial community were in agreement with in vitro [39] or in vivo [38, 40] studies.

Table 5

Effect of brewers’ spent grain (BSG) protein hydrolysates or antibiotics on microbial N synthesis and microbial community of FPA in RUSITEC

Item

Treatments1

SEM

P value

Con

AlcH

FlaH

Ant

Production of microbial N2, mg/d

           

LAB3

60.7a

60.7a

53.3b

48.2c

1.94

0.05

FPB4

7.0

6.2

6.4

7.0

1.04

0.51

FPA5

13.7c

16.6b

17.6ab

19.4a

1.82

0.01

Total

80.3a

82.5a

76.2b

73.5b

1.83

0.05

Efficiency of microbial protein6

7.0

7.1

6.9

6.9

0.26

0.45

a, b, c Least square means within a row with different superscripts differ (P < 0.05).
1 Con = Control, no antioxidant peptide and no antibiotics; AlcH = Alcalase hydrolysates, 10 mg per gram of TMR (DM basis); FlaH = Flavourzyme hydrolysates, 10 mg per gram of TMR (DM basis); Ant = antibiotics, 0.8 mg monensin + 1 mg tylosin per gram of TMR (DM basis).
2 Samples from d 14 to 15.
3 Liquid associate bacteria.
4 Feed particle-bound bacteria.
5 Feed particle-associated bacteria.
6 Efficiency of microbial protein, mg microbial N production/g OM fermented

Analysis of the taxonomic composition revealed that AlcH and FlaH had no affect on the relative abundance (RA) of bacteria at phylum level (Fig. 3A). Whereas, the addition of antibiotics reduced (P < 0.05) the RA of Firmicutes and Bacteroidetes, and enhanced (P < 0.05) the RA of Proteobacteria at phylum level. For all treatments, phyla of Firmicutes, Bacteroidetes and Proteobacteria were the most predominant phyla, accounting for approximately 90% of the microbial biomass in FPA, which was consistent with the sequencing results of liquid associated microbial community of in vitro [38] and in vivo studies [39, 40]. This suggested that phyla of Firmicutes, Bacteroidetes and Proteobacteria were predominant in both liquid and particle associated ruminal microbiome. While at the genus level, Prevotella, Succinivibrio, Selenomonas, Shuttleworthia, Schwartzia, Bifidobacterium, Lactobacillus, Dialister and Anaerobiospirillum were among the 10 most abundant genera in FPA microbial community (Fig. 3B) and should be considered as the “core bacteria” associated with feed particle. The FlaH enhanced (P < 0.05) the RA of Prevotella, but reduced the RA of Selenomonas, Shuttleworthia, Bifidobacterium and Dialister as compared to the control. The AlcH also had less (P < 0.05) RA of Schwartzia and Bifidobacterium than control. Notably, Ant reduced (P < 0.05) the RA of Prevotella, Selenomonas, Shuttleworthia, Bifidobacterium, Lactobacillus and Dialister, and increased RA of Succinivibrio compared with control. Further Log2 fold change analysis found that several genera from phyla Actinobacteria, Bacteroidetes, Firmicutes and Proteobacteria had 5% more change (enhanced or reduced) than control (Fig. 4). Furthermore, genera of Shuttleworthia, Selenomonas, Lactobacillus and Schwartzia from phyla Firmicutes were especially susceptible to FlaH, AlcH and Ant supplementation. These results indicated that FlaH and AlcH had less antibacterial effects as compared to Ant. As F. succinogenes, R. albus, R. flavefaciens, P. ruminicola, E. cellulosolvens, and E. ruminantium are major rumen fibrolytic bacteria [41], the present results indicated that application of BSG protein hydrolysates had no obvious detrimental effect on fibrolytic bacteria. This was supported by no treatment effect of FlaH and AlcH on the disappearances of NDF and ADF.

More susceptible to FlaH and AlcH for genera of Shuttleworthia, Selenomonas, Lactobacillus and Schwartzia from phyla Firmicutes and genera of Bifidobacterium from phyla Actinobacteria, suggested that the BSG protein hydrolysates with antioxidant activity were more potent to Gram-positive bacteria than to Gram-negative bacteria. The underlined mechanism is unknown, but it seemed that this difference might be attributed to the lack of the protective outer membrane in Gram-positive bacteria. As essential oils with antioxidant activity, had a stronger antibacterial activity than other essential oils [42], we speculated that the antibacterial effects of BSG protein hydrolysates might be due to their antioxidant activities. The antioxidant activity of protein hydrolysates relies on the enzymes and methods used during BSG protein hydrolysis. Similarly, differences in effectiveness of monensin on ruminal Gram-positive and Gram-negative bacteria were also observed in the current study, and confirm the monensin sensitivity to Gram-positive bacteria [38].

Rumen microbial community was usually dominated by Prevotella at genus level when fed high grain diets [43]. Members of Prevotella are able to utilise various nutrients such as starch, proteins and non-cellulosic polysaccharides [44, 45], and to convert lactate into propionate to prevent accumulation of lactate [40]. Moreover, as Prevotella are also H2-consuming bacteria in addition to Selenomonas [46], the reduced H2 with FlaH and AlcH was likely due to the increased relative abundance of genus Prevotella. Meanwhile, members of Prevotella are known as the predominant proteolytic bacteria with a great diversity of extracellular proteolytic activities in the rumen [31, 47], thus, the decreased disappearance of CP by supplementing monensin could be explained by decreased relative abundance of Prevotella. Currently, the probiotics such as direct-fed microbials are often developed from the genera Lactobacillus and Bifidobacterium [48], and the reduced relative abundance of Lactobacillus and Bifidobacterium with Ant suggested a disadvantage of applying monensin in the diet. The genus Schwartzia are asaccharolytic and can ferment succinic acid to produce propionic acid [49]. It was reported that bacteria affiliated with Schwartzia were negatively correlated with methane emissions [50], which was contrary to our results that monensin reduced the RA of Schwartzia, as well as methane emission.

There is limited information on the genus Shuttleworthia, and its function in the rumen is almost unknown [43]. Recent studies found that genus Shuttleworthia was digesta-adherent rumen bacteria in dairy and beef cattle [51], and was a starch and sugar fermenter [52]. It was reported that supplementation of phytogenic compounds decreased the RA of Shuttleworthia in dry cow [52]. In the current study, the decreased Shuttleworthia by BSG hydrolysates or Ant may explain the decreased starch digestibility with FlaH and Ant. A recent study showed negative correlation of Shuttleworthia with lactate and NH3-N concentration [53], thus explained the greater NH3-N production with BSG hydrolysates.

Bacteria from genus Selenomonas, besides of degrading starch and cellulose, play critical role in maintaining normal rumen fermentation through converting lactate and succinate into propionate and reducing lactate accumulation [5456] and maintaining low H2 concentration as H2-consumer [57, 58]. As fumarate and nitrate reducers, Selenomonas were proven a significant H2 sinks in sheep [59]. Ruminants fed high starch diets often result in great rumen dH2 concentration, and would promote the growth of S. ruminantium for disposing electrons derived from fermentation [57]. Therefore, in the present study, the reduced H2 production by applying BSG protein hydrolysates or monensin is consistent with the declined relative abundance of Selenomonas.

Members of genus Succinivibrio can ferment both starch and cellulose into succinate, and succinate then fermented by Selenomonas and other bacteria into propionate via succinate pathway, which is the primary propionate producing pathway in rumen [56]. Therefore, the abundances of members of Succinivibrio in the rumen has been reported to be positively associated with feed efficiency of ruminants [60, 61]. The increased RA of Succinivibrio with monensin in the present study, was correlated with the increased propionate proportion and the increased fermenter pH. Our results agreed with previous report that monensin can lead to an increase in abundance of succinate producers [38], and the report that Succinivibrio is positively correlated to the ruminal pH [55, 62]. The increased pH with Ant would be partly due to the increased RA of genus Dialister, which were reported to play a role in altering the buffering capacity of rumen fluid [63]. Furtherly, although the effects of monensin on methanogens were not determined, the application of monensin reduced the H2 production and CH4 emission. Our results supported that the decreased methane production with Ant most likely ascribed to decrease in nutrient availability for methanogenesis by acting on other rumen microorganisms instead of a direct effect of monensin on the methanogens [38].

Conclusion

Inclusion of FlaH with antioxidant activity to high grain diet can effectively protect CP and starch from ruminal degradation, without affecting fibre degradation, fermentation pH, and VFA profiles. Supplementation of FlaH also showed promising effects on reduced CH4 production by suppressing H2 generation. Whereas, the different effectiveness of FlaH and AlcH on altering in vitro rumen fermentation and gas production is likely due to the different pepetide compositions generated with different hydralysate enzymes. Our results suggested that FlaH had similar but mild effects to monensin. Enzymatically hydrolysis protein extracts from BSG with flavourzyme is a feasible way of processing brewing by-product, which can benefit both the brewing industry and animal production, as well as environmental profit by reducing wastes and greenhouse gas emission.

Abbreviations

A:P, Ratio of acetate to propionate; ADF, Acid detergent fibre; AlcH, Protein hydrolysates generated under Alcalase hydrolysis; BCVFA, Branched-chain volatile fatty acids (isobutyrate + isovalerate); BSG, Brewers’ spent grain; CP, Crude protein; DM, Dry matter; FlaH, Protein hydrolysates generated under Flavourzyme hydrolysis; FPA, Feed particle associated bacteria; FPB, Feed particle-bound bacteria; LAP, Liquid associated bacteria. NDF: Neutral detergent fibre; RA, Relative abundance; TMR, Total mixed ration; VFA, Volatile fatty acid; OM, Organic matter

Declarations

Acknowledgements

The authors thank the Lethbridge Research and Development Centre Metabolism barn staff for their care and management of the animals and Alastair Furtado and Darrell Vedres for their technical assistance.

Authors’ contributions

WZY and LYC designed the study and revised manuscript, RA performed BSG protein hydrolysates generation and antioxidant activity assay and revised manuscript; AMS and XMZ helped on samples collection and analysis; LJ and DYN helped for microbial data analysis and revised manuscript; TR performed this study, analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Funding

This work was financially supported by Agriculture and Agri-Food Canada (AAFC) Growing forward program (GF2#1542).

Availability of data and materials

The data analyzed during the current study are available from the corresponding author on reasonable request.

Ethics approval

All the procedures for the treatment and care of experimental cattle were approved by the Animal Care and Use Committee of Lethbridge Research and Development Centre and followed the guidelines for the Canadian Council on Animal Care (2009).

Declaration of Competing Interest

The authors declare that they do not have any conflict of interest.

Consent for publication

Not applicable.

References

  1. Mussatto SI, Dragone G, Roberto IC. Brewers' spent grain: generation, characteristics and potential applications. J Cereal Sci. 2006;43:1-14.
  2. Aliyu S, Bala M. Brewer’s spent grain: A review of its potentials and applications. Afr J Biotechnol. 2011;10:324-331.
  3. Parmenter RT, Richard JW, James DM. Case study: effects of inclusion of wet brewers’ grains on the growth performance, carcass characteristics, and meat quality of finishing cattle. Prof Anim Sci. 2018;34:505-512.
  4. Yang WZ, Li YL, McAllister TA, McKinnon JJ, Beauchemin KA. Wheat distillers’ grains in feedlot cattle diets: feeding behavior, growth performance, carcass characteristics, and blood metabolites. J Anim Sci. 2012;90:1301-1310.
  5. Dewell GA, Ransom JR, Dewell RD, McCurdy K, Gardner IA, Hill AE, Sofos JN, Belk KE, Smith GC, Salman MD. Prevalence of and risk factors for Escherichia coli O157 in market-ready beef cattle from 12 U.S. feedlots. Foodborne Pathog Dis. 2005;2:70-76.
  6. Yang HE, Yang WZ, McKinnon JJ, Alexander TW, Li YL, McAllister TA. Survival of Escherichia coli O157:H7 in ruminal or fecal contents incubated with corn or wheat dried distillers’ grains with solubles. Can J Microbiol. 2010;56:890-895.
  7. Hao XY, Benke MB, Gibb DJ, Stronks A, Travis G, McAllister TA. Effects of dried distillers’ grains with solubles (wheat-based) in feedlot cattle diets on feces and manure composition. J Environ Qual, 2009;38:1709-1718.
  8. Ikram S, Huang LY, Zhang HJ, Wang J, Yin M. Composition and nutrient value proposition of brewers spent grain. J Food Sci. 2017;82:2232-2242.
  9. Hou YQ, Wu ZL, Dai ZL, Wang GH, Wu GY. Protein hydrolysates in animal nutrition: industrial production, bioactive peptides, and functional significance. J Anim Sci Biotechnol. 2017;8:24.
  10. McCarthy AL, O'Callaghan YC, Piggott CO, FitzGerald RJ, O'Brien NM. Brewers’ spent grain; bioactivity of phenolic component, its role of animal nutrition and potential for incorporation in functional foods: a review. P Nutr Soc. 2013;72:117-125.
  11. Bamdad F, Wu JP, Chen LY. Effects of enzymatic hydrolysis on molecular structure and antioxidant activity of barley hordein. J Cereal Sci. 2011;54:20-28.
  12. Vieira EF, da Silva DD, Carmo H, Ferreira IMPLVO. Protective ability against oxidative stress of brewers' spent grain protein hydrolysates. Food Chem. 2017;228:602-609.
  13. Xia YC, Bamdad F, Gänzle M, Chen LY. Fractionation and characterization of antioxidant peptides derived from barley glutelin by enzymatic hydrolysis. Food Chem. 2012;134:1509-1518.
  14. McCarthy AL, Connolly A, Piggott CO, O'Brien NM, FitzGerald RJ, O'Callaghan YC. In vitro antioxidant and anti-inflammatory effects of brewers' spent grain protein rich isolate and its associated hydrolysates. Food Res Int. 2013;50:205-212.
  15. Shen YZ, Abeynayake R, Sun X, Ran T, Li JG, Chen LY, Yang WZ. Feed nutritional value of brewers' spent grain residue resulting from protease aided protein removal. J Anim Sci Biotechnol. 2019;10:78.
  16. Canadian Council on Animal Care (CCAC). Guide to the care and use of farm animals in research, teaching and testing. CCAC Ottawa, ON, 2009.
  17. Czerkawski JW, Breckenridge G. Design and development of a long-term rumen simulation technique (Rusitec). Br J Nutr. 1977;38:371-84.
  18. Saleem AM, Ribeiro Jr GO, Yang WZ, Ran T, Beauchemin KA, McGeough EJ, Ominski KH, Okine EK, McAllister TA. Effect of engineered biocarbon on rumen fermentation, microbial protein synthesis, and methane production in an artificial rumen (RUSITEC) fed a high forage diet. J Anim Sci. 2018;96:3121-3130.
  19. McDougall EI. Studies on ruminant saliva. 1. The composition and output of sheep’s saliva. Biochem J. 1948;43:99-109.
  20. Wang M, Sun XZ, Janssen PH, Tang SX, Tan ZL. Responses of methane production and fermentation pathways to the increased dissolved hydrogen concentration generated by eight substrates in in vitro ruminal cultures. Anim Feed Sci Tech. 2014;194:1-11.
  21. Walters W, Hyde ER, Berg-Lyons D, Ackermann G, Humphrey G, Parada A, et al. Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems, 2015;1:e00009-15.
  22. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852-857.
  23. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlin D, et al. Vegan: community ecology package. R package Version 2.4-4. 2017.
  24. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.
  25. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
  26. AOAC, Official Methods of Analysis. Vol. I. 18th ed. AOAC, Gaithersburg, MD. 2005.
  27. Van Soest PJ, Robertson JB, Lewis BA. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J Dairy Sci, 1991;74:3583-97.
  28. Rode LM, Yang WZ, Beauchemin KA. Fibrolytic enzyme supplements for dairy cows in early lactation. J Dairy Sci. 1999;82:2121-2126.
  29. Yu DJ, Sun YW, Wang WJ, O’Keefe SF, Neilson AP, Feng H, et al. Recovery of protein hydrolysates from brewer’s spent grain using enzyme and ultrasonication. I J Food Sci Tech. 2020;55:357-368.
  30. Connolly A, Cermeño M, Crowley D, O'Callaghan Y, O'Brien NM, Fitz RJ. Characterisation of the in vitro bioactive properties of alkaline and enzyme extracted brewers' spent grain protein hydrolysates. Food Res Int. 2019;121:524-532.
  31. Wallace RJ, Ruminal microbial metabolism of peptides and amino acids. J Nutr. 1996;126:1326S-1334S.
  32. Witt MW, Newbold CJ, Wallace RJ. Influence of dietary acetylated peptides on fermentation and peptidase activities in the sheep rumen. J Appl Microbiol. 1998;84:847-851.
  33. Bach A, Calsamiglia S, Stern MD. Nitrogen metabolism in the rumen. J Dairy Sci. 2005;88:E9-21.
  34. Bladen HA, Bryant MP, Doetsch RN. A study of bacterial species from the rumen which produce ammonia from protein hydrolyzate. Appl Microbiol. 1961;9:175-180.
  35. Russell JB, Strobel HJ. Effect of ionophores on ruminal fermentation. Appl Environ Microbiol. 1989;55:1-6.
  36. Bodas R, Prieto N, García-González R, Andrés S, Giráldez FJ, López S. Manipulation of rumen fermentation and methane production with plant secondary metabolites. Anim Feed Sci Tech. 2012;176:78-93.
  37. Russell JB. Rumen microbiology and its role in ruminant nutrition. Publ. Co., Ithaca, NY, USA., 2002.
  38. Schären M, Drong C, Kiri K, Riede S, Gardener M, Meyer U, et al. Differential effects of monensin and a blend of essential oils on rumen microbiota composition of transition dairy cows. J Dairy Sci. 2017;100:2765-2783.
  39. Shen JS, Liu Z, Yu ZT, Zhu WY. Monensin and nisin affect rumen fermentation and microbiota differently in vitro. Front Microbiol. 2017;8:1111.
  40. Thomas M, Webb M, Ghimire S, Blair A, Olson K, Fenske GJ, et al. Metagenomic characterization of the effect of feed additives on the gut microbiome and antibiotic resistome of feedlot cattle. Sci Rep. 2017;7:12257.
  41. Stewart CS, Flint HJ, Bryant MP. The rumen bacteria. In: Hobson PN, Stewart CS (eds) The rumen microbial ecosystem. Springer Netherlands, Dordrecht. 1997;10-72.
  42. Patra AK, Yu ZT. Essential oils affect populations of some rumen bacteria in vitro as revealed by microarray (RumenBactArray) analysis. Front Microbiol. 2015;6:297.
  43. Zhang RY, Ye HM, Liu JH, Mao SY. High-grain diets altered rumen fermentation and epithelial bacterial community and resulted in rumen epithelial injuries of goats. Appl Microbiol Biotechnol, 2017;101:6981-6992.
  44. Accetto T, Avgustin G. The diverse and extensive plant polysaccharide degradative apparatuses of the rumen and hindgut Prevotella species: A factor in their ubiquity? Syst Appl Microbiol. 2019;42:107-116.
  45. Purushe J, Fouts DE, Morrison M, White BA, Mackie RI, Coutinho PM, et al., Comparative genome analysis of Prevotella ruminicola and Prevotella bryantii: insights into their environmental niche. Microb Ecol. 2010;60:721-729.
  46. Denman SE, Fernandez GM, Shinkai T, Mitsumori M, McSweeney CS. Metagenomic analysis of the rumen microbial community following inhibition of methane formation by a halogenated methane analog. Front Microbiol. 2015;6:1087.
  47. Griswold KE, White BA, Mackie RI. Diversity of extracellular proteolytic activities among Prevotella species from the rumen. Curr Microbiol. 1999;39:187-194.
  48. McAllister TA, Beauchemin KA, Alazzeh AY, Baah J, Teather RM, Stanford K. Review: The use of direct fed microbials to mitigate pathogens and enhance production in cattle. Can J Anim Sci. 2011;91:193-211.
  49. van Gylswyk NO, Hippe H, Rainey FA. Schwartzia succinivorans gen. nov., sp. nov., another ruminal bacterium utilizing succinate as the sole energy source. Int J of Sys Bacteriol. 1997;47:155-159.
  50. Cunha CS, Veloso CM, Marcondes MI, Mantovani HC, Tomich TR, Pereira LGR, et al., Assessing the impact of rumen microbial communities on methane emissions and production traits in Holstein cows in a tropical climate. Syst Appl Microbiol. 2017;40:492-499.
  51. Noel SJ, Attwood GT, Rakonjac J, Moon CD, Waghorn GC, Janssen PH. Seasonal changes in the digesta-adherent rumen bacterial communities of dairy cattle grazing pasture. PLoS One. 2017;12:e0173819.
  52. Neubauer V, Petri R, Humer E, Kröger I, Mann E, Reisinger N, et al. High-grain diets supplemented with phytogenic compounds or autolyzed yeast modulate ruminal bacterial community and fermentation in dry cows. J Dairy Sci. 2018;101:2335-2349.
  53. Xia GL, Sun J, Fan YT, Zhao FF, Ahmed G, Jin YQ, et al. β-sitosterol attenuates high grain diet-induced inflammatory stress and modifies rumen fermentation and microbiota in sheep. Animals. 2020;10:171.
  54. Wang XX, Li XB, Zhao CX, Hu P, Chen H, Liu ZX et al., Correlation between composition of the bacterial community and concentration of volatile fatty acids in the rumen during the transition period and ketosis in dairy cows. Appl Environ Microbiol. 2012;78:2386-2392.
  55. Pan XH, Xue FG, Nan XM, Tang ZW, Wang K, Beckers Y, et al. Illumina sequencing approach to characterize thiamine metabolism related bacteria and the impacts of thiamine supplementation on ruminal microbiota in dairy cows fed high-grain diets. Front Microbiol. 2017;8:1818.
  56. Lan W, Yang CL. Ruminal methane production: Associated microorganisms and the potential of applying hydrogen-utilizing bacteria for mitigation. Sci Total Environ, 2019;654:1270-1283.
  57. Wang M, Wang R, Liu M, Beauchemin KA, Sun XZ, Tang SX, et al. Dietary starch and rhubarb supplement increase ruminal dissolved hydrogen without altering rumen fermentation and methane emissions in goats. Animal. 2019;13:975-982.
  58. Belanche A, Doreau M, Edwards JE, Moorby JM, Pinloche E, Newbold CJ. Shifts in the rumen microbiota due to the type of carbohydrate and level of protein ingested by dairy cattle are associated with changes in rumen fermentation. J Nutr. 2012;142:1684-1692.
  59. Greening C, Geier R, Wang C, Woods LC, Morales SE, McDonald MJ, et al. Diverse hydrogen production and consumption pathways influence methane production in ruminants. ISME J. 2019;13:2617-2632.
  60. Hernandez-Sanabria E, Goonewardene LA, Wang ZQ, Durunna ON, Moore SS, Guan LL. Impact of feed efficiency and diet on adaptive variations in the bacterial community in the rumen fluid of cattle. Appl Environ Microbiol. 2012;78:1203-1214.
  61. Elolimy AA, Arroyo JM, Batistel F, Iakiviak MA, Loor JJ. Association of residual feed intake with abundance of ruminal bacteria and biopolymer hydrolyzing enzyme activities during the peripartal period and early lactation in Holstein dairy cows. J Anim Sci Biotechnol. 2018;9:43.
  62. Wetzels SU, Mann E, Metzler-Zebeli BU, Pourazad P, Qumar M, Klevenhusen F, Pinior B, et al. Epimural indicator phylotypes of transiently-induced subacute ruminal acidosis in dairy cattle. Front Microbiol. 2016;7:274.
  63. Myer PR, Smith TPL, Wells JE, Kuehn LA, Freetly HC. Rumen microbiome from steers differing in feed efficiency. PLoS One, 2015;10:e0129174.