Effect of different buffering capacities on in vitro rumen fermentation parameters
The buffering capacity of BC0.9% and BC0.5% were significantly greatest (P < 0.05) after 24 h incubation compared to BC0.3% and CON (Table 1). Both BC0.9% and BC0.5% exhibited significantly highest (P < 0.05) buffering capacity value of 106.00 meq/L, hence had a similar effect on in vitro after 24 h. The ruminal pH obtained from BC0.9% and BC0.5% showed similar effects and were consistently higher (P < 0.05) than the other treatments throughout the incubation period. In gas production, BC0.9%, BC0.5%, and BC0.3% had significantly higher (P < 0.05) gas produced than CON and showed similar effects at 3 and 12 h incubation. Ruminal NH3-N concentration was significantly higher (P < 0.05) in BC0.9% and BC0.5%, thus, it also had a similar effect on this parameter. However, at 6 and 12 h, no effect observed on treatments except that BC0.5% tended to increase (P = 0.073) NH3-N concentration followed by BC0.9% and the rest treatments.
Significantly higher concentrations (P < 0.05) of acetate were observed in BC0.9% at 12 h; however, BC0.5% and BC0.3% obtained the highest value (P < 0.05) after 24 h (Table 2). Propionate and butyrate concentrations were both highest (P < 0.05) in BC0.3% and BC0.9% at 6 h. Subsequently, distinct effects of BC0.3%, BC0.5%, and BC0.9% were observed at 24 h which had significantly higher (P < 0.05) propionate concentrations than CON. A similar pattern was noticeable with butyrate at 12 h such that BC0.3%, BC0.5%, and BC0.9% obtained the highest concentration (P < 0.05) compared with CON. During this period, a similar effect can be seen between the 3 treatments; however, no significant effect was observed after 24 h. Total volatile fatty acid contents were greater (P < 0.05) in BC0.3%, BC0.5% and BC0.9% at 12 h but had a slight change after 24 h. At this time point, treatments BC0.3% and BC0.5% were highest (P < 0.05) compared to BC0.9% and CON. Furthermore, there were no treatment effects on acetate to propionate ratio after 24 h incubation. Consequently, increasing the concentration of buffering capacity showed linear effects (P < 0.05) on pH, total gas production, NH3-N, and at some certain time point of individual VFA.
Table 1 Effect of different buffering capacity concentrations on in vitro rumen fermentation parameters at 3, 6, 12 and 24 h
Parameters
|
Time (h)
|
Treatment e
|
SEM
|
P-value
|
CON
|
BC0.3%
|
BC0.5%
|
BC0.9%
|
All
|
Linear
|
Buffering capacity (meq/L)
|
3
|
76.44d,z
|
83.89c,y
|
87.56b,x
|
90.11a
|
0.398
|
<0.001
|
<0.001
|
6
|
76.45c,z
|
84.17b,y
|
92.89a,x
|
92.89a
|
0.246
|
<0.001
|
<0.001
|
12
|
85.33c,y
|
87.22c,y
|
96.67b,x
|
99.22a
|
0.469
|
<0.001
|
<0.001
|
24
|
100.22c,z
|
102.56b,y
|
106.00a,x
|
106.00a
|
0.462
|
<0.001
|
<0.001
|
pH
|
3
|
6.00c
|
6.03bc
|
6.08ab
|
6.13a
|
0.013
|
0.012
|
0.052
|
6
|
5.68c,z
|
5.80b,y
|
5.91a,x
|
5.93a
|
0.014
|
<0.001
|
<0.001
|
12
|
5.42c,z
|
5.54b,y
|
5.60a,x
|
5.62a
|
0.014
|
<0.001
|
<0.001
|
24
|
5.14b,y
|
5.15b,y
|
5.22a,x
|
5.24a
|
0.009
|
0.001
|
0.003
|
Total gas (mL)
|
3
|
74.67b,y
|
82.00a,x
|
81.67a,x
|
82.00a
|
1.287
|
0.016
|
0.004
|
6
|
122.33
|
124.33
|
123.00
|
123.67
|
1.353
|
0.777
|
0.488
|
12
|
169.00b,y
|
187.00a,x
|
179.00a,x
|
182.67a
|
2.492
|
0.007
|
0.010
|
24
|
251.67
|
269.00
|
256.00
|
256.67
|
3.877
|
0.078
|
0.195
|
NH3-N (mg/dL)
|
3
|
11.26c,y
|
12.24b,x
|
13.05a,x
|
13.39a
|
0.205
|
0.001
|
0.003
|
6
|
13.35
|
13.60
|
16.09
|
14.74
|
0.569
|
0.178
|
0.095
|
12
|
13.90y
|
14.95xy
|
16.23x
|
15.84
|
0.482
|
0.073
|
0.010
|
24
|
19.90c
|
20.39bc
|
21.35ab
|
22.43a
|
0.378
|
0.012
|
0.052
|
a-d Means with different superscripts in a row differ significantly (P < 0.05)
e CON (no buffer added); BC0.3% (0.3% buffer); BC0.5% (0.5% buffer); BC0.9% (0.9% buffer)
x,y,z Means within a row indicate linear effect among CON, BC0.3%, and BC0.5% (P < 0.05)
Table 2 Volatile fatty acid production during in vitro rumen fermentation incubated at 3, 6, 12, and 24 h
Parameters
|
Time (h)
|
Treatment e
|
SEM
|
P-value
|
CON
|
BC0.3%
|
BC0.5%
|
BC0.9%
|
All
|
Linear
|
Acetate (mmol/L)
|
3
|
100.13
|
94.81
|
93.41
|
90.36
|
2.505
|
0.265
|
0.237
|
6
|
100.03
|
100.31
|
99.79
|
100.39
|
0.592
|
0.901
|
0.775
|
12
|
102.52b
|
104.34b
|
103.36b
|
107.76a
|
0.291
|
0.001
|
0.412
|
24
|
103.71c,y
|
112.05ab,xy
|
117.27a,x
|
108.47bc
|
2.182
|
0.018
|
0.009
|
Propionate (mmol/L)
|
3
|
27.44
|
26.51
|
24.02
|
25.59
|
1.395
|
0.598
|
0.280
|
6
|
27.81b
|
28.67a
|
27.66b
|
28.37a
|
0.085
|
0.001
|
0.421
|
12
|
32.19
|
32.65
|
30.59
|
33.85
|
0.756
|
0.227
|
0.350
|
24
|
35.24b,y
|
40.32a,xy
|
42.31a,x
|
38.09ab
|
1.203
|
0.045
|
0.020
|
Butyrate (mmol/L)
|
3
|
17.60
|
20.15
|
19.60
|
19.42
|
0.501
|
0.088
|
0.098
|
6
|
22.23a,x
|
22.78a,x
|
20.95b,y
|
22.81a
|
0.276
|
0.009
|
0.006
|
12
|
27.46b,y
|
40.52a,x
|
39.69a,x
|
40.08a
|
1.660
|
0.007
|
<0.001
|
24
|
53.85
|
54.58
|
56.43
|
54.94
|
0.783
|
0.327
|
0.125
|
Total VFA (mmol/L)
|
3
|
145.16
|
141.47
|
137.02
|
135.37
|
3.381
|
0.430
|
0.285
|
6
|
150.07
|
151.76
|
148.40
|
151.57
|
0.803
|
0.077
|
0.196
|
12
|
162.18b,y
|
177.51a,x
|
173.63a,x
|
181.69a
|
2.424
|
0.005
|
0.005
|
24
|
192.80c,y
|
206.95ab,xy
|
216.01a,x
|
201.51bc
|
3.447
|
0.013
|
0.008
|
A:P ratio
|
3
|
3.74
|
3.57
|
3.89
|
3.53
|
0.134
|
0.618
|
0.670
|
6
|
3.60ab
|
3.50c
|
3.61a
|
3.54bc
|
0.018
|
0.014
|
0.719
|
12
|
3.19
|
3.20
|
3.41
|
3.18
|
0.084
|
0.465
|
0.266
|
24
|
2.94
|
2.78
|
2.78
|
2.85
|
0.069
|
0.434
|
0.178
|
a-c Means with different superscripts in a row differ significantly (P < 0.05)
e CON (no buffer added); BC0.3% (0.3% buffer); BC0.5% (0.5% buffer); BC0.9% (0.9% buffer)
x,y Means within a row indicate linear effect among CON, BC0.3%, and BC0.5% (P < 0.05)
Effect of different buffering capacities on rumen fermentation characteristics in Hanwoo steers
The effect of different buffering capacity concentrations on rumen fermentation characteristics of Hanwoo steers in four treatments are presented in Table 3. Average pH had no significant effects among CON and treatments. However, BC0.3%, BC0.5%, and BC0.9% had significantly higher (P < 0.05) buffering capacity value than CON, and showed linearly significant effect (P < 0.05). Ammonia-nitrogen, acetate to propionate ratio, individual and total VFA concentrations of rumen fluid from steers under all treatments were not significant and showed similar effects after the in vivo experiment.
Table 3 Effect of different buffering capacity concentrations on rumen fermentation characteristics in Hanwoo steers
Parameters
|
Treatment c
|
SEM
|
P-value
|
CON
|
BC0.3%
|
BC0.5%
|
BC0.9%
|
All
|
Linear
|
Average pH
|
6.28
|
6.56
|
6.77
|
6.48
|
0.146
|
0.226
|
0.528
|
Buffering capacity (meq/L)
|
83.47b,y
|
99.53a,x
|
101.00a,x
|
94.80a
|
3.058
|
0.015
|
0.004
|
NH3-N (mg/dL)
|
4.84
|
4.26
|
5.40
|
4.59
|
1.260
|
0.945
|
0.794
|
Total VFA (mmol/L)
|
47.24
|
45.29
|
54.16
|
55.87
|
5.234
|
0.632
|
0.563
|
Acetate (mmol/L)
|
29.67
|
27.28
|
30.70
|
35.56
|
2.421
|
0.532
|
0.890
|
Propionate (mmol/L)
|
11.30
|
9.26
|
14.37
|
12.75
|
2.190
|
0.412
|
0.380
|
Butyrate (mmol/L)
|
6.28
|
8.76
|
9.08
|
7.56
|
1.556
|
0.652
|
0.299
|
A:P ratio
|
2.81
|
3.10
|
2.13
|
2.91
|
0.451
|
0.431
|
0.370
|
a-b Means with different superscripts in a row differ significantly (P < 0.05)
c CON (no buffer added); BC0.3% (0.3% buffer); BC0.5% (0.5% buffer); BC0.9% (0.9% buffer)
x,y Means within a row indicate linear effect among CON, BC0.3%, and BC0.5% (P < 0.05)
Bacterial diversity of the rumen contents of Hanwoo steers
The boxplot representation of alpha diversity indices is shown in Figure 1. Alpha diversity indices are composite indices that reflect abundance and consistency. Chao1 which reflect the OTU abundance in the samples showed that BC0.9% was the highest among treatments followed by BC0.5% and the rest of the treatments (Figure 1a). Shannon index which reflects the diversity of the OTU in samples presented BC0.9% as the most diverse among treatments and BC0.3% being the least (Figure 1b). Moreover, Figure 1c showed the boxplot of OTUs of observed species from the samples. The number of OTUs in BC0.9% was higher followed by BC0.5% and the rest of the treatments. The diversity index is used to analyze the temporal and spatial changes in species composition which reflects whether bacterial communities between groups have differences. Our results showed that the rumen bacterial composition of BC0.5% and BC0.9% had overall higher alpha diversity than other treatment groups, although no significant difference was observed after statistical analysis.
Effect of treatments on bacterial community composition of Hanwoo steers rumen contents
Bacterial taxonomic compositions at the phylum, genera, and species level are shown in Figure 2. Results at the phylum level revealed that 15 bacterial phyla were identified in the rumen digesta samples of Hanwoo steers (Figure 2a). The majority of the sequences obtained from all treatments belonged to Bacteroidetes followed by Firmicutes. It was noticeable that BC0.3% had the highest abundance of Bacteroidetes (71.90%) and lowest Firmicutes (22.13%). On the contrary, BC0.9% had the lowest abundance of Bacteroidetes (54.19%) among treatments; however, BC0.5% had the highest Firmicutes (33.84%) relative abundance. Furthermore, Proteobacteria was dominant upon incorporating BC0.9% in the diet resulting in declining its abundance from that of BC0.3% and BC0.5%. Spirochaetes increased its abundance when animals received supplementation of BC0.9%, whereas there was a sudden decrease once the steers received BC0.5% in their diet. At the genus level (Figure 2b), Prevotella was the predominant among genera in all treatments followed by Bacteroides (10.36%) and Ruminococcus (8.93%) in BC0.3% and BC0.5%, respectively. Furthermore, Succiniclasticum was dominant among genera after supplementing BC0.5%; however, a decreasing abundance of Paludibacter was noticeable as BC0.5% is incorporated into the diet, which had a reverse effect as did CON. Species-level analyses revealed that Prevotella ruminicola predominated the treatments CON, BC0.3%, BC0.5% and BC0.9% with the relative abundance of 24.85%, 32.16%, 26.73%, and 23.17%, respectively (Figure 2c). The comparison of single species analyzed through statistical analysis showed a significant effect of the treatments only in the case of Prevotella brevis. This species was more abundant (P = 0.015) in the CON and as steers received a diet supplemented with BC0.9%, BC0.5% and BC0.3% its abundance decreased. Owing to the BC0.5% supplemented in the diet, a decreasing abundance of Paludibacter propionicigenes was observed; however, it increased in CON. Incorporation of BC0.5% increased the microbial population of Ruminococcus bromii and Succiniclasticum ruminis. Moreover, the smaller percentage of BC0.3% resulted in a higher abundance of Bacteroides massiliensis which led to a sudden decrease in its population as the concentration of treatments increases. Supplementing buffers of different buffering capacity concentration may affect the rumen microbiota through the relative abundance of bacterial species.
The core, shared and unique bacterial community of observed species of the rumen microbiome after treatment of buffer agents with varying level of buffering capacity is presented in Figure 3 as Venn diagram. A total of 211 (59.6%) observed species can be found across all the samples (core), 79 (22.32%) for shared by 2 or 3 samples, and 64 (18.08%) are specific and are distributed to the four samples.
The normalized data presented in Figure 4 shows the clustering based on the similarity of relative abundance between representative families of OTUs (row), and treatments (column). The analysis divided the representative families into two major clusters distinguishing families which represents low relative abundance on all treatments (upper cluster in red), and families that have varying relative abundance between treatments (lower cluster, colored from peach to blue). On the cluster presenting varying abundance between treatments, two sub-clusters were also distinguishable; (1) families which represent variation from very low (red) to average (peach) abundance, and (2) cluster representing families which have average to high (blue) abundance. The cluster in the bottom of the heatmap (labelled) contains the families that represent the above average relative abundance. On this cluster, family Prevotellaceae had branched out because it presents the highest abundance with very small variations between treatments (p = 0.092). Family Ruminicoccaceae were also found in all treatments, but varying relative abundance was observed, with BC0.5% presenting the highest. Families Acidaminococcaceae and Lachnospiraceae were significantly highest (P < 0.05) in BC0.3% and BC0.5%, respectively. Also, the Unclassified Clostridiales had significantly highest (P < 0.05) relative abundance in BC0.5%. A certain unclassified family under order Bacteroidales also showed major abundance especially in BC0.3%, while families Vibrionaceae and Spirochaetaceae were highest in BC0.9%. Meanwhile, the normalized data of relative abundance of representative families of observed OTUs have grouped the control sample together with BC0.3% in a single cluster, while BC0.5% and BC0.9% are on their own cluster, which could indicate the effect of varying buffering capacity of buffer agents.
The comparison of the bacterial communities by principal coordinate analysis (PCoA) is presented in Figure 5. The PCoA plots showed close similarity within and between treatments and control, whereas those under BC0.9% and BC0.3% groups showed dispersed distribution of bacterial communities. The PCoA plot showed dissimilarity of bacterial community and revealed a distinct structure between high buffering capacity and low buffering capacity.
Monitoring of acidosis
The changes in the 24 h mean ruminal pH monitored for 30 d is presented in Table 4. During this period of the feeding challenges, mean pH values were >5.8. Minimum pH was lowest in CON, whereas it was highest in BC0.5%. Additionally, BC0.9% had a low minimum pH value second to that of CON. It was noticeable that BC0.3% and BC0.5% had higher minimum and mean pH values compared to BC0.9% and CON. Obtained results indicated that the duration of time where pH was <5.8 and 5.8 < 6.0 was longer in CON followed by BC0.9% and BC0.3%. Meanwhile, BC0.3% also exhibited good results in the duration of time where pH was approximately 6.0 and above; however, BC0.5% had even better effects and did not show any signs of acidosis in the rumen. Based on the data gathered, BC0.5% stabilized the pH of rumen preventing it from becoming acidotic.
Table 4. Changes in the 24 h mean ruminal pH, duration of time where pH was <5.8 monitored for 30 d challenge diet
Item
|
Treatment a
|
CON
|
BC0.3%
|
BC0.5%
|
BC0.9%
|
24 h mean ruminal pH
|
|
|
|
|
Minimum
|
5.28
|
5.72
|
5.83
|
5.40
|
Mean
|
6.36
|
6.56
|
6.47
|
6.22
|
Duration of ruminal pH
|
|
|
|
|
pH <5.8, min/d
|
66.29
|
1.45
|
0.00
|
13.06
|
pH 5.8 < 6.0, min/d
|
90.97
|
21.77
|
10.16
|
106.45
|
pH 6.0 and above, min/d
|
1214.52
|
1296.77
|
1202.42
|
999.68
|
a CON (no buffer added); BC0.3% (0.3% buffer); BC0.5% (0.5% buffer); BC0.9% (0.9% buffer)