Animals diet
Bulls were fed a basal diet comprised of 70% concentrate containing corn grain offered ad libitum and protein supplement (soybean meal; premix composed of: urea, vitamins and minerals; limestone; yeast and salt) and 30% corn silage for 62 days (Table 1).
Table 1
Ingredients and chemical composition of basal diet (g/kg DM)
Ingredients | Diet |
Corn silage | 275.9 |
Corn grain | 613.2 |
Soybean meal | 51.0 |
Premix1 | 50.5 |
Mineral salt | 4.5 |
Limestone | 4.5 |
Yeast | 0.4 |
Chemical composition | |
Dry matter | 577 |
Crude protein | 132 |
Organic matter | 968 |
Ash | 31.4 |
Ether extract | 40.1 |
Neutral detergent fibre | 288 |
Acid detergent fibre | 117 |
Total digestible nutrients | 790 |
Metabolisable energy (MJ/kg DM) | 11.9 |
Calcium | 6.82 |
Phosphorus | 3.56 |
1Premix: magnesium (57 g/kg), sodium (81 g/kg), sulphur (3.75 g/kg), cobalt (20 mg/kg), copper (500 mg/kg), iodine (25 mg/kg), manganese (1 500 mg/kg), selenium (10 mg/kg), zinc (2 000 mg/kg), vitamin A (400 000 UI/kg), vitamin D3 (50 000 UI/kg), vitamin E (750 UI/kg), ether extract (168 g/kg) and urea (200 g/kg). |
Feeding behavior activities
There were no effects of NA blend addition to bull diets on rumination, feed intake, water intake and idle time (P > 0.05; Table 2).
|
Experimental diets
|
SEM6
|
P – value
|
Activities, min/day
|
CON1
|
NA152
|
NA303
|
NA454
|
NA605
|
Rumination
|
245.0
|
219.5
|
209.5
|
262.0
|
245.0
|
9.911
|
0.550
|
Feed intake
|
343.5
|
349.5
|
344.5
|
305.5
|
337.5
|
9.182
|
0.394
|
Water Ingestion
|
35.0
|
34.5
|
38.0
|
32.0
|
37.0
|
2.451
|
0.932
|
Idle
|
816.5
|
836.5
|
848.0
|
840.5
|
820.5
|
11.392
|
0.883
|
¹CON = control (without natural additives); 2NA15 – addition of 1.5 g/animal/day of natural additives; 3NA30 – addition of 3.0 g/animal/ day of natural additives; 4NA45 – addition of 4.5 g/animal/day of natural additives; 5NA60 – addition of 6.0 g/animal/day of natural additives. Naturals additives contained clove leaf essential oil, castor and cashew functional oils and a commercial blend composed of vanillin, eugenol and thymol; 6Standard error of means; 7Linear effect; 8Quadratic effect.
|
Table 2
Feeding behavior from young bulls with and without natural additive addition to the diet
Animal performance
The initial body weight and final body weight (FBW) were similar for all diets (P > 0.05), nonetheless average daily gain (ADG) of bulls increased linearly (P < 0.01) when the NA blend was added in the diets (Table 3). The addition of NA in the diets had no effect (P > 0.05) on Dry Matter Intake (DMI) (kg/day – 9.9 or kg/100 kg body weight – 2.3%). However, feed efficiency improved linearly (P < 0.04) with the addition of NA to the diets (Table 3). In addition, the HCW (Hot Carcass Weight) and HCD (Hot Carcass Dressing) did not differ between cattle fed with blend of NA (P > 0.05; Table 3).
|
Experimental diets
|
|
P – value
|
Item
|
CON1
|
NA152
|
NA303
|
NA454
|
NA605
|
SEM6
|
L7
|
Q8
|
0% vs blend
|
Initial body weight, kg
|
382.8
|
388.0
|
385.6
|
385.4
|
387.3
|
2.941
|
0.762
|
0.641
|
0.623
|
Final body weight, kg
|
473.0
|
478.7
|
481.4
|
486.9
|
490.0
|
3.942
|
0.131
|
0.322
|
0.267
|
Average daily gain, kg
|
1.43
|
1.44
|
1.52
|
1.61
|
1.63
|
0.031
|
0.013
|
0.047
|
0.145
|
Dry matter intake, kg/d
|
9.85
|
9.80
|
9.83
|
10.12
|
10.09
|
0.144
|
0.300
|
0.521
|
0.706
|
Dry matter intake, %/BW
|
2.30
|
2.26
|
2.27
|
2.32
|
2.33
|
0.024
|
0.542
|
0.670
|
0.909
|
Feed efficiency, kg
|
0.145
|
0.147
|
0.155
|
0.160
|
0.160
|
0.014
|
0.043
|
0.134
|
0.216
|
Hot carcass weight, kg
|
248.1
|
252.0
|
246.6
|
253.9
|
246.1
|
2.521
|
0.900
|
0.879
|
0.816
|
Hot carcass dressing, %
|
52.37
|
52.62
|
51.25
|
52.18
|
51.51
|
0.302
|
0.178
|
0.195
|
0.357
|
¹CON = control (without natural additives); 2NA15 – addition of 1.5 g/animal/day of natural additives; 3NA30 – addition of 3.0 g/animal/ day of natural additives; 4NA45 – addition of 4.5 g/animal/day of natural additives; 5NA60 – addition of 6.0 g/animal/day of natural additives. Naturals additives contained clove leaf essential oil, castor and cashew functional oils and a commercial blend composed of vanillin, eugenol and thymol; 6Standard error of means; 7Linear effect; 8Quadratic effect.
|
Table 3
Animal performance of young bulls with and without natural additive addition to the diet
Ruminal ammonia and volatile fatty acid (VFA).
The addition of a blend of NA affected rumen fermentative characteristics and resultant VFAs produced (Table 4). The major VFAs: acetate, butyrate, isovalerate, and valerate were reduced linearly when animals were fed NA (P < 0.05). Similarly, NA addition in the diets linearly reduced (P < 0.02) the acetate/propionate ratio. NA supplementation of diets resulted in a quadratic effect on propionate and isobutyric acid concentrations (P < 0.05). Furthermore, animals supplied with NA had linear reductions in rumen methane concentration (P < 0.001). Ammonia concentration had a quadratic effect following NA blend supplementation of bull diets (P < 0.001). The ruminal pH was not affected (P > 0.05) by inclusion of NA in the diets (Table 4).
|
Experimental diets
|
|
P – value
|
Item
|
CON1
|
NA152
|
NA303
|
NA454
|
NA605
|
SEM6
|
L7
|
Q8
|
0% vs blend
|
Acetate (mmol/mol)
|
56.15
|
56.16
|
43.64
|
43.98
|
43.74
|
1.31
|
<.0001
|
<.0001
|
<.0001
|
Propionate (mmol/mol)
|
17.45
|
17.00
|
14.44
|
16.37
|
13.69
|
0.73
|
0.350
|
0.054
|
0.682
|
Isobutyric (mmol/mol)
|
0.91
|
1.18
|
0.85
|
0.79
|
0.93
|
0.03
|
<.0001
|
0.038
|
0.623
|
Butyrate (mmol/mol)
|
10.87
|
13.89
|
8.67
|
6.33
|
7.30
|
0.67
|
<.0001
|
0.221
|
0.262
|
Isovaleric (mmol/mol)
|
3.07
|
3.75
|
2.08
|
1.85
|
2.39
|
0.18
|
0.002
|
0.055
|
0.144
|
Valeric (mmol/mol)
|
1.23
|
1.33
|
0.94
|
0.92
|
1.07
|
0.06
|
0.018
|
0.210
|
0.226
|
A/P* ratio
|
3.22
|
3.37
|
3.02
|
2.73
|
3.24
|
0.12
|
0.023
|
0.945
|
0.434
|
Ammonia (mg/dL)
|
21.82
|
5.95
|
5.94
|
3.02
|
4.20
|
1.72
|
0.006
|
<.0001
|
<.0001
|
pH
|
6.91
|
6.95
|
7.05
|
6.95
|
7.07
|
0.06
|
0.270
|
0.968
|
0.326
|
¹CON = control (without natural additives); 2NA15 – addition of 1.5 g/animal/day of natural additives; 3NA30 – addition of 3.0 g/animal/ day of natural additives; 4NA45 – addition of 4.5 g/animal/day of natural additives; 5NA60 – addition of 6.0 g/animal/day of natural additives. Naturals additives contained clove leaf essential oil, castor and cashew functional oils and a commercial blend composed of vanillin, eugenol and thymol; 6Standard error of means; 7Linear effect; 8Quadratic effect. *A/P = acetate/propionate ratio.
|
Table 4
Ruminal parameters of young bulls with and without natural additive addition to the diet
Rumen bacterial diversity and abundance
In our study, the major phyla present in the rumen were Bacteroidetes (47%) and Firmicutes (36%; Figure 1). Bacteroidetes (P < 0.05) were reduced when NA was included in the diet. A quadratic response was seen for Candidatus Saccharibacteria, Chytridiomycota, Elusimicrobia, Eukaryota Unassigned, Fibrobacteres, Firmicutes, Spirochaetes, Synergistetes and Tenericutes (P < 0.05). Source data are included in supplementary material (Table S1).
The families Prevotellaceae (43%) and Ruminococcaceae (20%) were observed as the most abundant across treatments (Figure 2). Significant changes were observed in the families causing quadratic responses in Cardiobacteriaceae, Clostridiales_Family_XIII_Incertae_Sedis, Prevotellaceae, Ruminococcaceae. Our data also showed a decrease in Acidaminococcaceae, Coriobacteriaceae, Defluviitaleaceae, Desulfovibrionaceae, Neisseriaceae, Paenibacillaceae, Peptococcaceae, Porphyromonadaceae and an increase in Christensenellaceae, Bacillaceae, Lactobacillaceae, Ophryoscolecidae, Rikenellaceae, Trichomonadidae (P < 0.05) post NA supplementation of bull diets. Source data are included in supplementary material (Table S2).
The most common rumen bacterial genera across the treatments were Succinivibrio, Succiniclasticum, Marvinbryantia and Prevotella (12%, 11%, 9% and 6%, respectively; Figure 3). A quadratic effect was observed when NA was supplemented into the bull diet with respect to the genera Alistipes, Asteroleplasma, Dorea, Elusimicrobium, Entodinium, Faecalibacterium, Haemophilus, Holdemanella, Paraprevotella, Pseudoscardovia, Pyramidobacter, Roseburia, Ruminobacter, Sphaerochaeta, Subdoligranulum, Syntrophococcus. A decrease in Acetitomaculum, Acidaminococcus, Akkermansia, Alloprevotella, Candidatus_Saccharimonas, Citreitalea, Clostridium, Fretibacterium, Mailhella, Moryella, Phascolarctobacterium, Prevotella, Robinsoniella, Succinimonas, Suttonella, Tetratrichomonas and Weissella and an increase in Anaerostipes, Atopobium, Bacillus, Bavariicoccus, Fibrobacter, Hydrogenoanaerobacterium, Paenibacillus and Sporobacter (P < 0.05) was noted post NA dietary supplementation. Source data are included in supplementary material (Table S3).
Methanogen diversity and abundance
Archaeal abundance was reduced on the whole with the inclusion of NA in the bull diets (P < 0.05; Table 5). The families Methanobacteriaceae and Methanomicrobiaceae (P < 0.05); orders Methanomicrobiales, Methanobacteriales and Methanomassiliicoccales (P < 0.05) and the genera Methanobrevibacter and Methanosphaera, showed a significant decrease with NA supplementation, whilst the genus Methanomicrobium showed a tendency to be present at lower abundance (P = 0.051). Furthermore, on a species level, a decrease in Methanobrevibacter ruminantium, Methanobrevibacter sp D5 and Methanobrevibacter sp G16 was seen following NA supplementation of bull diets (P < 0.05).
Table 5
Archaea diversity and abundances from young bulls with and without natural additive addition to the diet
| Experimental diets | | P – value |
Archaea taxonomy | CON1 | NA152 | NA303 | NA454 | NA605 | SEM6 | L7 | Q8 | 0% vs blend |
Archaea Euryarchaeota | 2.00 | 2.22 | 2.08 | 1.74 | 1.93 | 0.422 | 0.434 | 0.847 | 0.977 |
f_Methanobacteriaceae | 88.53 | 17.29 | 19.61 | 18.40 | 14.56 | 13.800 | 0.956 | 0.918 | < .0001 |
f_Methanomicrobiaceae | 0.27 | 0.02 | 0.03 | 0.00 | 0.06 | 0.060 | 0.826 | 0.739 | 0.002 |
o_Methanomicrobiales | 21.89 | 2.74 | 5.77 | 4.09 | 3.94 | 4.794 | 0.844 | 0.692 | 0.005 |
o_Methanobacteriales | 19.84 | 2.90 | 3.34 | 3.87 | 2.48 | 3.442 | 0.845 | 0.991 | < .0001 |
o_Methanomassiliicoccales | 1.66 | 0.13 | 0.27 | 0.02 | 0.17 | 0.207 | 0.728 | 0.440 | < .0001 |
g_Methanobrevibacter | 211.22 | 42.48 | 36.31 | 40.08 | 43.89 | 14.733 | 0.909 | 0.786 | < .0001 |
g_Methanomicrobium | 0.74 | 0.04 | 0.23 | 0.05 | 0.18 | 0.264 | 0.981 | 0.557 | 0.051 |
g_Methanosphaera | 7.56 | 2.12 | 1.97 | 1.83 | 2.63 | 1.255 | 0.869 | 0.995 | < .0001 |
s_Methanobrevibacter ruminantium | 0.72 | 0.04 | 0.43 | 0.20 | 0.21 | 0.118 | 0.336 | 0.044 | 0.001 |
s_Methanobrevibacter sp D5 | 0.98 | 0.40 | 0.24 | 0.12 | 0.26 | 0.139 | 0.162 | 0.936 | < .0001 |
s_Methanobrevibacter sp G16 | 0.74 | 0.05 | 0.13 | 0.04 | 0.11 | 0.262 | 0.983 | 0.791 | 0.039 |
¹CON = control (without natural additives); 2NA15 – addition of 1.5 g/animal/day of natural additives; 3NA30 – addition of 3.0 g/animal/ day of natural additives; 4NA45 – addition of 4.5 g/animal/day of natural additives; 5NA60 – addition of 6.0 g/animal/day of natural additives. Naturals additives contained clove leaf essential oil, castor and cashew functional oils and a commercial blend composed of vanillin, eugenol and thymol; 6Standard error of means; 7Linear effect; 8Quadratic effect; f_ = family taxonomy, g_ genus taxonomy; o_ = order taxonomy; s_= species taxonomy. |
Gene Network correlations
We observed close to 13,000 functionally annotated genes in total across the experimental samples using shotgun metagenomics and 28 were significantly differentially abundant when the bull diet contained NA (Fig. 4; Table S4). Functional annotation data showed significantly biological responses due to the NA addition whereas mostly related to protection against foreign attack to DNA and DNA maintenance, replication and repair (Restriction endonuclease, type II, XamI; Restriction endonuclease, type II, EcoRV; Host-nuclease inhibitor protein Gam; RecG, N-terminal antiparallel four helix bundle; Type IV secretion system protein TraG/VirD4; Type IV secretion system, VirB10 / TraB / TrbI and Transposase, ISC1217). There were also functional process associated with membrane protection and maintenance (ABC-2 transporter; Conjugal transfer, TrbG/VirB9/CagX and Capsule biosynthesis protein CapC), metabolic role (Lyase, catalytic; Acyl-protein synthetase, LuxE; Phenolic acid decarboxylase, bacterial; Peptidase G2, IMC autoproteolytic cleavage domain; Glycyl radical enzyme, HI0521, predicted; Transposase, ISC1217 and Tetrahydrodipicolinate-N-succinyltransferase, chain A, domain 1), oxidative stress response (Thiol peroxidase conserved site and Aldehyde oxidase/xanthine dehydrogenase, molybdopterin binding), attack protection and resistance (Bacterial virulence protein VirB8; KorB, C-terminal and Siphovirus Gp157), plasmid replication (KorB, C-terminal), and unknown biologic process (Protein of unknown function DUF4244; Protein of unknown function DUF4054; Protein of unknown function DUF4912; Protein of unknown function DUF4294; Protein of unknown function DUF4416; Protein of unknown function DUF3853).
Specifically, the functional annotations Restriction endonuclease, type II, XamI; Lyase, catalytic; Acyl-protein synthetase, LuxE; Host-nuclease inhibitor protein Gam; ABC-2 transporter; Transposase, ISC1217; RecG, N-terminal antiparallel four helix bundle and Protein of unknown function DUF4294 were decreased with NA inclusion in the diet. Furthermore, the annotations that showed an increase post NA inclusion in the diet were: Glycyl radical enzyme, HI0521, predicted; Aldehyde oxidase/xanthine dehydrogenase, molybdopterin binding, Peptidase G2, IMC autoproteolytic cleavage domain; Siphovirus Gp157; Type IV secretion system protein TraG/VirD4; Type IV secretion system, VirB10 / TraB / TrbI; Conjugal transfer, TrbG/VirB9/CagX; KorB, C-terminal and Protein of unknown function DUF4416. Nevertheless, a quadratic response was also noted for: Bacterial virulence protein VirB8; Capsule biosynthesis protein CapC; Phenolic acid decarboxylase, bacterial; Restriction endonuclease, type II, EcoRV; Thiol peroxidase conserved site; Tetrahydrodipicolinate-N-succinyltransferase, chain A, domain 1; Protein of unknown function DUF3853 and Protein of unknown function DUF4912.
The family Succinivibrionaceae had a strong positive correlation (average r = > 0.9) with Tetrahydrodipicolinate-N-succinyltransferase, chain A, domain 1; Type IV secretion system, VirB10 / TraB / TrbI; Phenolic acid decarboxylase, bacterial; Thiol peroxidase conserved site; Type IV secretion system, VirB10 / TraB / TrbI; Bacterial virulence protein VirB8; Conjugal transfer, TrbG/VirB9/CagX; KorB, C-terminal gene abundances. The Paenibacillaceae bacterial family (Phylum Firmicutes) had a positive correlation (r = > 0.9) with Peptidase G2 and Glycyl radical enzyme, HI0521, predicted gene abundance. The Victivallaceae interacted with Protein Function DUF4416 and Capsule Biosynthesis Protein CapC (r = > 0.9). The Glycyl radical enzyme, HI0521, predicted, showed a major correlation with Bacillaceae (r = > 0.9). Prevotellaceae had a negative correlation (r = − 0.8) with Ruminococcaceae, and Methanobacteriaceae also had a negative correlation (r = > − 0.7) with Protein Function DUF4294 and ABC-2 transporter gene abundances. Source data are included in supplementary material (Table S4, Fig. 6).