RESEARCH ARTICLE

Influence of dietary organic trace minerals on enteric methane emissions and rumen microbiota of heat-stressed dairy steers

A-Rang Son1,#https://orcid.org/0000-0002-3896-0950, Mahfuzul Islam1,2,#https://orcid.org/0000-0001-6010-140X, Seon-Ho Kim1https://orcid.org/0000-0002-9350-1853, Sung-Sill Lee3https://orcid.org/0000-0002-4621-4333, Sang-Suk Lee1,*https://orcid.org/0000-0003-1540-7041
Author Information & Copyright
1Ruminant Nutrition and Anaerobe Laboratory, Department of Animal Science and Technology, Sunchon National University, Suncheon 57922, Korea
2Department of Microbiology and Parasitology, Sher-e-Bangla Agricultural University, Dhaka 1207, Bangladesh
3Institute of Agriculture and Life Science and University-Centered Laboratory, Gyeongsang National University, Jinju 52828, Korea

# These authors contributed equally to this work.

*Corresponding author: Sang-Suk Lee, Ruminant Nutrition and Anaerobe Laboratory, Department of Animal Science and Technology, Sunchon National University, Suncheon 57922, Korea. Tel: +82-61-750-3237, E-mail: rumen@scnu.ac.kr

© Copyright 2023 Korean Society of Animal Science and Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Jul 28, 2022; Revised: Nov 01, 2022; Accepted: Nov 02, 2022

Published Online: Jan 31, 2023

Abstract

Ruminants are the main contributors to methane (CH4), a greenhouse gas emitted by livestock, which leads to global warming. In addition, animals experience heat stress (HS) when exposed to high ambient temperatures. Organic trace minerals are commonly used to prevent the adverse effects of HS in ruminants; however, little is known about the role of these minerals in reducing enteric methane emissions. Hence, this study aimed to investigate the influence of dietary organic trace minerals on rumen fermentation characteristics, enteric methane emissions, and the composition of rumen bacteria and methanogens in heat-stressed dairy steers. Holstein (n=3) and Jersey (n=3) steers were kept separately within a 3×3 Latin square design, and the animals were exposed to HS conditions (Temperature-Humidity Index [THI], 82.79 ± 1.10). For each experiment, the treatments included a Control (Con) consisting of only basal total mixed rations (TMR), National Research Council (NRC) recommended mineral supplementation group (NM; TMR + [Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm]/kg dry matter), and higher concentration of mineral supplementation group (HM; basal TMR + [Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm]/kg dry matter). Higher concentrations of trace mineral supplementation had no influence on methane emissions and rumen bacterial and methanogen communities regardless of breed (p > 0.05). Holstein steers had higher ruminal pH and lower total volatile fatty acid (VFA) concentrations than Jersey steers (p < 0.05). Methane production (g/d) and yield (g/kg dry matter intake) were higher in Jersey steers than in Holstein steers (p < 0.05). The relative abundances of Methanosarcina and Methanobrevibacter olleyae were significantly higher in Holstein steers than in Jersey steers (p < 0.05). Overall, dietary organic trace minerals have no influence on enteric methane emissions in heat-stressed dairy steers; however, breed can influence it through selective alteration of the rumen methanogen community.

Keywords: Dietary minerals; Enteric methane; Heat stress; Holstein and Jersey steers; Rumen methanogens

INTRODUCTION

Enteric CH4 is considered not only an indicator of gross energy losses, but also a potent greenhouse gas with a global warming potential that is 28 times higher than that of carbon dioxide (CO2) [1,2]. Enteric CH4 is produced by methanogenic archaea through methanogenesis following one of three pathways: hydrogenotrophic, methylotrophic, or acetoclastic [35]. In light of this, several dietary mitigation strategies have been conducted to reduce enteric CH4 emissions from ruminants [5,6]; however, little is known about the role of mineral supplementations in achieving this. Li et al. [7] reported decreased CH4 emissions and methanogen populations in lactating cows following dietary mineral salt supplementation. They explained that dietary mineral salt reduced CH4 emission through the reduction of methanogen phenotypes and the A : P ratio, while increasing H+ ion utilization in propionate production. However, that study did not focus on the effects of minerals on heat-stressed dairy cattle or steers.

Ruminant production is greatly hampered by adverse environmental conditions, especially heat stress (HS) [8-10]. HS reduces dry matter intake (DMI), which subsequently hampers energy and protein metabolism, leading to increased metabolic disorder, mineral imbalance, and several other health problems [11-14]. The sensitivity and response to HS varies among breeds [15,16]. Holstein steers have been reported to exhibit a significant reduction in DMI under HS [16]. Additionally, studies have found that both Jersey steers and Jersey dairy cows are less sensitive to HS than Holstein steers [15,16]. Therefore, both Holstein and Jersey steers were considered in this study. Furthermore, the influence of HS on rumen microbial alteration and enteric methane (CH4) emissions in Holstein and Jersey steers has already been reported [16]. During HS, cattle require more supplementation with trace minerals in their diet to prevent the adverse effects of HS [17]. Some trace minerals, such as zinc (Zn), copper (Cu), and selenium (Se), have been used in diets to minimize the adverse effects of HS in ruminants because of their antioxidant effects [14,1823]. However, organic trace minerals have a more beneficial effect than inorganic minerals on ruminants [18]. Therefore, organic trace minerals were considered in this study. Although the National Research Council (NRC) [24] has recommended the dose of these minerals under normal conditions, the optimal dose required to overcome the adverse effects of HS during summer has not yet been established. Hence, during HS, we supplied minerals at a level of 70% based on the maximum tolerable concentration recommended by the NRC [24]. It has been previously reported that different trace minerals have toxic effects on methanogens. Hernández-Sánchez et al. [25] reported that Cu decreases CH4 production because it is toxic to some rumen methanogens. Liu et al. [26] further reported that mineral supplementation reduces enteric methane emissions by altering the rumen microbiome. Dietary Se improves rumen fermentation by altering the rumen microbiome; however, it also has antimicrobial effects that may affect rumen methanogen diversity [27]. Therefore, this study evaluated the effects of supplementation with organic forms of Zn, Cu, and Se minerals on enteric CH4 emissions, rumen fermentation characteristics, rumen bacteria, and methanogens in heat-stressed Holstein and Jersey steers.

MATERIALS AND METHODS

The study was conducted at the Sunchon National University (SCNU) animal farm and the ruminant nutrition and anaerobic laboratory in the Department of Animal Science and Technology, SCNU, Suncheon, Korea. All animals used in this study and all experimental protocols were reviewed and approved by the Institutional Animal Care and Use Committee (approval number: SCNU-IACUC-2020-06).

Animals, experimental design, and diet

Three non-cannulated Holstein steers (710.33 ± 43.02 kg; approximately 30 months old) and three Jersey steers (559.67 ± 32.72 kg; approximately 30 months old), were kept in two separate areas with a 3×3 Latin square design. For both experiments, the diets included a Control group (Con) fed the basal total mixed rations (TMR) without mineral supplementation, an NRC recommended concentration of mineral supplementation group (NM) fed TMR with (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg dry matter (DM), and a higher concentration of mineral supplementation group (HM) fed the basal TMR with (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM. Organic form of Se (Yeast-Selenium; X-SEL 3000TM, Algebra Bio, New South Wales, Australia), Zn-glycinate (BASF SE, Ludwigshafen, Germany), and Cu-glycinate (BASF SE) were used to supplement Se, Zn, and Cu, respectively.

The duration of the feeding experiment for each period was 20 d. This comprised diet adaptation and reaming for the first 14 days and CH4 measurement and rumen fluid sampling for the following 6 days. A 7 d washing period was maintained between each Latin square. The ingredients and chemical composition of the basal TMR are presented in Table 1. Individual stalls with feeding and water facilities were available for the experimental steers. The Con group was fed only basal TMR once a day at 09:00 h with a refusal rate of 5%–10%, whereas the respective concentrations of minerals were mixed well with the basal TMR and fed to the NM and HM groups. The DMI was measured as the difference between the feed offered and feed refused. Basal TMR was collected twice (on days 7 and 14) during the experiment. Dry matter content was estimated using a hot-air oven at 65°C for 72 h [28]. Proximate analysis of TMR was performed using standard methods [29]. The protocols of Van Soest et al. [30] and Van Soest [31] were used to determine the neutral detergent fiber (NDF) and acid detergent fiber (ADF) content, respectively. This study was conducted under the condition of HS with an average temperature-humidity index (THI) of 82.79 ± 1.10, which was calculated as THI = (0.8 × ambient temperature) + [% relative humidity/100 × (ambient temperature − 14.4)] + 46.4 [32].

Table 1. Chemical composition of total mixed ration (TMR)
Ingredients Compositions (% of DM)
Corn gluten feed 8.40
Soybean 6.24
Beet pulp 4.20
Wheat bran 3.15
Corn flakes 2.21
Molasses 1.04
Rice wine residue 5.25
Brewer’s grain residue 21.01
Annual ryegrass straw 27.29
Orchard grass straw 21.01
Limestone 0.10
Sodium bicarbonate 0.01
Salt 0.09
Total 100.00
Chemical composition (DM basis)
 DM (fresh basis) 58.98%
 CP 13.55%
 Crude fiber 21.92%
 Crude fat 3.02%
 Ash 9.21%
 Calcium 1.22%
 Phosphorus 0.47%
 NDF 48.00%
 ADF 25.36%
 Zinc 77.35 ppm
 Copper 17.31 ppm
 Selenium 0.05 ppm

DM, dry matter; CP, crude protein; NDF, neutral detergent fiber; ADF, acid detergent fiber.

Download Excel Table
Measurement of enteric CH4 emission

Enteric CH4 emissions were measured using an automated head chamber system (AHCS) or GreenFeed (GF) unit (C-Lock, Rapid City, SD, USA), as described by Hristov et al. [33] with slight modifications. Briefly, all steers were trained to familiarize themselves with the GF unit before the experiment started to avoid any sort of psychological stress. To measure CH4 emissions, each steer visited the GF at eight different time points within three consecutive days in each measurement period. The 0 h or before feeding at 9 am, 9 h after feeding (6 pm), and 18 h after feeding (3 am) time points were considered for CH4 measurement on day 1, while 3 h after feeding (12 pm), 12 h after feeding (9 pm), and 21 h after feeding (6 am) were considered on day 2, and 6 h after feeding (3 pm) and 15 h after feeding (12 am) were considered on day 3. The GF unit was installed in one corner of a large pen, and at each measurement time point, all steers were successively moved to this pen from their individual stalls. To attract the animals to the GF unit and ensure a proper head-down position within the hood at the time of measurement, molasses-coated concentrated pellets (250–300 g/visit) were used. The number of pellets ingested per steer per day was excluded from DMI calculation. All relevant data of animal entry and exit times to the GF unit, standard gas calibration information, CO2 recovery time point, and amount of gas release data were sent to C-Lock. The calculated CH4 production (g/d) data were obtained using a web-based data management system, and the CH4 yield (g/kg DMI) was calculated.

Sample collection and the processing and recording of rectal temperature

In each period, rumen fluid was collected from each steer on the last day of the feeding trial, before the morning meal, using a stomach tube. The first 300 mL of rumen fluid was discarded to avoid contamination of rumen fluid by saliva. Ruminal pH was measured immediately after collection using a pH meter (SevenCompactTM pH/Ion meter S220, Mettler Toledo, Greifensee, Switzerland). Along the sides, three aliquots were prepared separately from each rumen fluid and transported to the laboratory in the presence of dry ice. These samples were stored at -80°C for further analysis of ammonia nitrogen (NH3-N), volatile fatty acids (VFA), and rumen microbiota. The rectal temperature (RT) of the steers was also recorded at approximately 12 pm on the same day of sampling using a digital thermometer (WPT-1, CAS, Yangju, Korea).

Analyses of ruminal NH3-N and volatile fatty acid concentrations

NH3-N concentration was measured using a Libra S22 spectrophotometer (CB40FJ, Biochrom, Cambourne, UK) according to the protocol described by Chaney and Marbach [34]. VFA concentration was measured using high-performance liquid chromatography (HPLC; Agilent Technologies 1200 series, Agilent Technologies, Waldbronn, Germany) according to the protocol described by Han et al. [35]. To perform HPLC, a UV detector (set at 210 nm and 220 nm), METACARB87H column (Varian, Palo Alto, CA, USA), and buffered solvent (0.0085 N H2SO4; at a flow rate of 0.6 mL/min) were used.

DNA extraction and metataxonomic analysis

For DNA extraction and subsequent metataxonomic analysis of rumen microbiota, all rumen fluid samples were sent to Macrogen (Seoul, Korea). Briefly, DNA from rumen fluid was extracted using a DNeasy PowerSoil Kit (Qiagen, Hilden, Germany), according to the manufacturer’s protocol [36]. The quality and quantity of DNA were assessed using PicoGreen and NanoDrop, respectively. To prepare amplicon libraries of each sample for both bacteria and archaea, two separate sequence runs were performed with two different primer sets. In order to prepare the amplicon library of each sample for bacteria, the Illumina 16S Metagenomic Sequencing Library protocols were used which was performed using two-step polymerase chain reaction (PCR) amplification of the 16S rRNA genes with the primers Bakt_341F (5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3′) and Bakt_805R (5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGAC AGGACTACHVGGGTATCTAATCC-3′) of the V3–V4 region at an annealing temperature of 55°C [37]. For archaea, 787-F (5′-TCGTCGGCAGCGTCAGAT GTGTATAAGAGACAGATTAGATACCCSBGTAGTCC-3′) and 1059-R (5′- GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGCCATGCACCWCCTCT -3′) primer sets of–V5-V6 were used, with an annealing temperature of 63°C. AMPure beads (Agencourt Bioscience, Beverly, MA, USA) were used to purify the products of the first and second PCR. The individual amplicon libraries were normalized after quantification using PicoGreen. They were then size verified using a TapeStation DNA ScreenTape D1000 (Agilent Technologies), pooled at an equimolar ratio, and sequenced on a MiSeq system (Illumina, San Diego, CA, USA) using a 2 × 300 bp kit. After sequencing, Illumina MiSeq raw data were classified by sample using an index sequence, and a paired-end FASTQ file was generated for each sample. The sequencing adapter sequence and F/R primer sequence of the target gene region were removed using Cutadapt (v3.2) [38].

For error correction of the amplicon sequencing process, the DADA2 (v1.18.0) [39] package of the R (v4.0.3) program was used. For paired-end reads of bacterial sequences, the forward sequence (Read1) and reverse sequence (Read2) were cut to 250 bp and 200 bp, respectively, and sequences with expected errors of two or more were excluded. However, 200 bp and 150 bp were considered for the same archaeal sequences. An error model for each batch was then established to remove noise from each sample. After assembling the paired-end sequence corrected for sequencing error into one sequence, the chimera sequence was removed using the DADA2 consensus method to form amplicon sequence variants (ASVs). In addition, for the comparative analysis of the microbial community, the QIIME (v1.9) [40] program was used for normalization by applying subsampling based on the number of reads of the sample with the minimum number of reads among all samples.

For each ASVs sequence, BLAST+ (v2.9.0) [41] was performed in the Reference DB (NCBI 16S Microbial DB for bacteria and NCBI NT DB for archaea), and taxonomic information for the organism of the subject with the highest similarity was assigned. At this time, if the query coverage of the best-hit matching the DB is less than 85% or the identity of the matched area is less than 85%, the taxonomy information is not allocated. A comparative analysis of various microbial communities was performed using QIIME with the above ASVs abundance and taxonomic information. The Shannon index and inverse Simpson index were obtained to check the species diversity and uniformity of the microbial community in the sample, and the alpha diversity information was confirmed using the rarefaction curve and Chao1 value. Based on the weighted and unweighted UniFrac distances, beta diversity between samples (information on microbial community diversity among samples in the comparison group) was obtained, and the relationship between samples was visualized using Principal Coordinates Analysis (PCoA) [40].

Statistical analysis

All data for DMI, CH4 emissions, rumen fermentation, and rumen microbiome were analyzed using the mixed procedure of SAS. Here, we considered breed and trace minerals as factors. Then, we tested whether there were any breed differences using a general linear model along with Duncan’s multiple range test. All analyses were performed using SAS (version 9.4, SAS Institute, Cary, NC, USA) [42]. Statistical significance was set at p < 0.05.

RESULTS

Dry matter intake, enteric CH4 emission, and rumen fermentation characteristics

The DMI and enteric CH4 emissions of Holstein and Jersey steers with varying mineral supplementation levels are presented in Table 2. In both breeds, DMI (kg/d) did not differ significantly among the treatment groups (p > 0.05). CH4 production (g/d) and CH4 yield (g/kg DMI) were significantly higher in Jersey steers than in Holstein steers (p < 0.05); however, a numerical decrease in CH4 production and yield was observed with increasing concentrations of mineral supplementation in both breeds (p > 0.05). The highest RT was recorded in the Con group, followed by the NM and HM groups in both breeds; however, the differences were not statistically significant (p > 0.05). The rumen fermentation characteristics of Holstein and Jersey steers supplemented with different levels of minerals are presented in Table 3. Ruminal pH was significantly higher in Holstein steers than that in Jersey steers (p < 0.05); however, a similar pH was observed among treatment groups (p > 0.05). The NH3-N concentration (mg/dL) was similar between breeds (p > 0.05); however, it was influenced by trace minerals and the interaction between breeds and trace minerals (p < 0.05). In Holstein steers, the Con group had the lowest NH3-N concentration, whereas the HM group had the highest NH3-N concentration when compared with Jersey steer groups (p < 0.05). The total VFA concentration (mmol/L) was higher in Holstein steers than in Jersey steers (p < 0.05); however, no differences were observed in the concentrations of total VFA or molar proportions of propionate and butyrate among different mineral supplemented groups in both breeds (p > 0.05). Although trace mineral supplementation had a significant influence on the molar proportion of acetate, regardless of breed (p < 0.05), no significant differences were observed between groups in each breed (p > 0.05).

Table 2. DMI and enteric methane emission of Holstein and Jersey steers with different levels of mineral supplementation
Parameters Holstein Jersey SEM Mixed p-value
Con NM HM Con NM HM B T B × T
DMI (kg/d) 12.74 13.32 12.90 11.02 11.29 11.37 1.172 0.125 0.947 0.982
CH4 production (g/d) 170.53 157.22 147.75 219.08 184.79 189.63 18.238 0.032 0.379 0.867
CH4 yield (g/kg DMI) 13.58 12.66 11.79 20.23 16.37 17.25 2.399 0.023 0.557 0.840
RT (°C) 39.33 39.30 39.10 39.13 39.10 39.03 0.115 0.185 0.452 0.853

Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM; B, breed effect; T, trace mineral supplementation effect; B × T, interaction effect between breed and trace mineral supplementation.

DMI, dry matter intake; CH4, methane; RT, rectal temperature; TMR, total mixed rations; NRC, National Research Council; DM, dry matter.

Download Excel Table
Table 3. Rumen fermentation characteristics of Holstein and Jersey steers with different levels of mineral supplementation
Parameters Holstein Jersey SEM Mixed p-value
Con NM HM Con NM HM Breed Treat B × T
pH 6.98 7.18 7.11 6.85 6.94 6.93 0.091 0.044 0.340 0.865
NH3-N (mg/dL) 4.18c 5.01a 4.82b 4.25y 4.42y 5.48x 0.035 0.228 < .0001 < .0001
Total VFA (mmol/L) 66.95 68.78 67.09 86.89 89.32 89.34 2.516 < .0001 0.781 0.924
Acetate (%) 57.52 61.37 56.20 56.94 57.81 55.18 0.928 0.135 0.049 0.466
Propionate (%) 26.23 22.79 24.06 25.85 23.53 23.89 1.359 0.967 0.315 0.945
Butyrate (%) 16.25 15.84 19.74 17.21 18.66 20.93 1.471 0.305 0.178 0.857
A:P 2.23 2.69 2.34 2.21 2.46 2.34 0.129 0.584 0.213 0.765

Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM; B, breed effect; T, trace mineral supplementation effect; B × T, interaction effect between breed and trace mineral supplementation.

a-c Means with different superscripts in a row differ significantly among different treatment groups in Holstein steers while

x-z Means with different superscripts in a row differ significantly among different treatment groups in Jersey steers (p < 0.05).

A:P, acetate: propionate; NH3-N, ammonia-nitrogen; VFA, volatile fatty acid; TMR, total mixed rations; NRC, National Research Council; DM, dry matter.

Download Excel Table
Species richness, diversity, and composition of rumen microbiota

A total of 531,527 bacterial and 1,326,280 archaeal quality-filtered sequence reads were obtained from sequencing 18 rumen fluid samples. The average sequences obtained from each sample for bacteria and archaea were greater than 29,500 and 73,600, respectively. Good’s coverage was greater than 99% for each sample. The Holstein steers showed tentatively higher ASV and Chao1 richness estimates (p = 0.051 and 0,052, respectively) and significantly higher Shannon and inverse Simpson diversity indices compared to Jersey steers, regardless of trace mineral supplementation (p < 0.05) (Table 4). Rarefaction measures for rumen bacteria are presented in Fig. 1. The beta diversity data were not significant between breeds or among the trace mineral supplementation groups (Fig. 2). However, the above-mentioned parameters were not influenced by trace mineral supplementation regardless of breed or the interaction between breeds and trace mineral supplementation (p > 0.05) (Table 4). The relative abundance of most bacterial phyla was not influenced by breed, trace mineral supplementation, or the interaction between them (p > 0.05). Jersey steers had a higher relative abundance of the phylum Candidatus Melainabacteria than Holstein steers (p < 0.05); however, no differences were observed among the different trace mineral supplementation groups in either breed (p > 0.05). In Holstein steers, at the phylum level, Bacteroidetes (accounting for 58.82% to 74.78%) and Firmicutes (19.07% to 36.79%) were the two major bacterial taxa among all treatment groups (Fig. 3, Table 5). Likewise, Bacteroidetes (which ranged from 66.21% to 69.63%) and Firmicutes (23.83% to 26.62%) were the top two bacterial phyla among the treated Jersey steers (Fig. 3). At the genus level, most of the more abundant bacterial genera were not influenced by breed, trace mineral supplementation, or the interaction between them (p > 0.05); however, the genus Capnocytophaga was more abundant in Jersey steers than in Holstein steers, regardless of trace mineral supplementation (p < 0.05) (Fig. 4; Table 6).

jast-65-1-132-g1
Fig. 1. Rarefaction measures of the rumen bacterial community of Holstein and Jersey steers with different levels of mineral supplementation. Hol, Holstein; Jer, Jersey.
Download Original Figure
jast-65-1-132-g2
Fig. 2. Principal coordinates analysis plot based on weighted (left) and unweighted-unifrac distance (right) showing the rumen bacterial community of Holstein and Jersey steers with different levels of mineral supplementation. PCoA, principal coordinates analysis; PC, principal component; Hol, Holstein; Jer, Jersey.
Download Original Figure
jast-65-1-132-g3
Fig. 3. Rumen bacterial abundance at the phylum level in Holstein and Jersey steers with different levels of mineral supplementation. Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM. TMR, total mixed rations; DM, dry matter.
Download Original Figure
jast-65-1-132-g4
Fig. 4. Rumen bacterial abundance at the genus level in Holstein and Jersey steers with different levels of mineral supplementation. Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM, and HM: TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM. TMR, total mixed rations; DM, dry matter.
Download Original Figure
Table 4. Species richness and diversity of rumen bacteria in Holstein and Jersey steers with different levels of mineral supplementation
Parameters Holstein Jersey SEM Mixed p-value
CON NM HM CON NM HM B T B × T
ASVs 1,329.00 1,385.33 1,387.00 1,279.33 1,272.00 1,189.00 56.581 0.051 0.836 0.564
Chao1 1,332.60 1,392.55 1,391.04 1,284.44 1,275.88 1,190.71 57.306 0.052 0.823 0.562
Shannon 8.91 9.05 9.17 8.74 8.86 8.67 0.105 0.019 0.605 0.395
Inverse Simpson 0.995 0.995 0.996 0.992 0.993 0.992 0.001 0.023 0.781 0.737
Good’s coverage 1.000 1.000 1.000 1.000 1.000 1.000 0.000 0.846 0.533 0.741

Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM; B, breed effect; T, trace mineral supplementation effect; B × T, interaction effect between breed and trace mineral supplementation.

ASV, amplicon sequence variant; TMR, total mixed rations; NRC, National Research Council; DM, dry matter.

Download Excel Table
Table 5. Rumen bacterial abundance at the phylum level in Holstein and Jersey steers with different levels of mineral supplementation
Phylum Holstein Jersey SEM p-value
Con NM HM Con NM HM Breed Treat B × T
Bacteroidetes 74.78 58.82 65.79 69.63 66.21 67.82 3.866 0.6871 0.1106 0.3591
Firmicutes 19.07 36.79 28.08 23.83 26.62 24.66 4.048 0.4277 0.1058 0.2732
Candidatus melainabacteria 1.24 0.40 1.09 1.33 1.54 1.72 0.266 0.0351 0.401 0.2889
Proteobacteria 1.20 0.56 0.97 1.01 1.03 1.25 0.298 0.4969 0.5659 0.5974
Tenericutes 1.19 1.42 1.44 1.60 1.70 1.95 0.398 0.2782 0.7965 0.9654
Spirochaetes 1.13 0.97 1.44 1.35 1.40 1.40 0.271 0.3912 0.6989 0.7138
Fibrobacteres 0.68 0.63 0.68 0.72 0.89 0.56 0.171 0.6841 0.7338 0.5889

Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM.

TMR, total mixed rations; NRC, National Research Council; DM, dry matter.

Download Excel Table
Table 6. Rumen bacterial abundance at the genus level in Holstein and Jersey steers with different levels of mineral supplementation
Genus Holstein Jersey SEM p-value
Con NM HM Con NM HM Breed Treat B × T
Prevotella 37.76 33.67 35.35 38.82 39.15 39.34 3.870 0.374 0.922 0.8908
Paraprevotella 9.27 7.91 8.48 9.87 5.60 6.99 1.415 0.3936 0.1956 0.6106
Capnocytophaga 1.36 1.69 1.85 3.95 2.74 4.39 0.656 0.0119 0.5814 0.6069
Paludibacter 5.44 6.58 5.66 9.07 8.53 4.43 1.277 0.2423 0.2059 0.2695
Ethanoligenens 3.06 3.11 2.53 1.73 2.20 2.88 0.491 0.1922 0.8429 0.3303
Holdemania 0.38 0.83 0.94 1.21 4.62 4.87 1.196 0.0846 0.4755 0.6513
Succiniclasticum 3.12 2.15 1.98 4.29 1.72 1.96 0.691 0.7285 0.0867 0.6167
Parabacteroides 2.16 1.69 0.80 1.08 1.41 1.02 0.676 0.5632 0.6196 0.711
Barnesiella 1.64 0.56 2.32 1.60 1.09 1.25 0.466 0.6788 0.2364 0.3924
Ruminococcus 1.86 4.13 2.19 1.53 1.86 2.12 0.607 0.1549 0.222 0.281
Lentimicrobium 1.14 2.89 2.75 1.38 1.05 0.97 0.669 0.1192 0.6604 0.3846
Olivibacter 1.54 0.27 0.40 0.84 0.49 0.17 0.335 0.5676 0.1741 0.6565

Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM.

TMR, total mixed rations; NRC, National Research Council; DM, dry matter.

Download Excel Table

Methanobrevibacter (ranging from 57.30% to 72.99% in Holstein steers and from 58.17% to 82.04% in Jersey steers) was the most abundant methanogen among all treatment groups for both breeds (Fig. 5A; Table 7). The relative abundance of Methanobrevibacter was not influenced by breed and trace mineral supplementation (p > 0.05) but was tentatively influenced by the interaction between breed and trace mineral supplementation (p = 0.065). Methanomassiliicoccus (15.74% to 34.57% in Holstein and 14.20% to 37.75% in Jersey steers) was the second most abundant methanogen among all treatment groups for both breeds. The relative abundance of Methanomassiliicoccus was not influenced by breed (p > 0.05); however, it was tentatively influenced by trace mineral supplementation (p = 0.061) and the interaction between breed and trace mineral supplementation (p = 0.072). Jersey steers showed a tendency toward decreasing Methanobrevibacter abundance and increasing Methanomassiliicoccus abundance with higher concentrations of trace mineral supplements (p = 0.072 and 0.086, respectively). However, in Holstein steers, only the lowest abundance of Methanobrevibacter and highest abundance of Methanomassiliicoccus were recorded in the HM group (p > 0.05). Among the remaining genera, Methanosarcina and Methanobacterium were more abundant, while Methanomicrobium was less abundant in the Holstein steers than in Jersey steers, regardless of the trace mineral supplementation (p < 0.05); the relative abundance of these methanogens was not influenced by trace mineral supplementation or the interaction between breed and trace mineral supplementation (p > 0.05). At the species level, the methanogens Methanobrevibacter thaueri, Mbr. olleyae, Mbr. millerae, and Methanomassiliicoccus luminyensis were the top four methanogen species among all the treatment groups in both breeds (Fig. 5B, Table 8). Among them, Mbr. olleyae were more abundant in Holstein steers than in Jersey steers, regardless of trace mineral supplementation (p < 0.05). Relative abundance of Mma. luminyensis was tentatively influenced by trace mineral supplementation and the interaction between breed and trace mineral supplementation (p = 0.061 and 0.072, respectively). A tendency toward an increasing pattern of Mma. luminyensis abundance was observed in Jersey steers with increasing mineral concentration (p = 0.086). Among the remaining methanogens, the Methanosarcina mazei, Mbr. oralis, and Methanobacterium aggregans were more abundant in Holstein steers, whereas the relative abundances of Mbr. boviskoreani, and Methanomicrobium mobile were greater in Jersey steers, regardless of trace mineral supplementation (p < 0.05).

jast-65-1-132-g5
Fig. 5. Rumen methanogen abundance at the genus (a) and species (b) levels in Holstein and Jersey steers with different levels of mineral supplementation. Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM. B, T, and B × T indicate significant (p < 0.05) difference in relative abundance between breeds, trace mineral supplementation, and the interaction between breed and trace mineral supplementation, respectively while b, t, and b × t indicate tentatively significant (0.05 < p < 0.1) difference in relative abundance between breed, trace mineral supplementation and the interaction between breed and trace mineral supplementation, respectively. TMR, total mixed rations; DM, dry matter.
Download Original Figure
Table 7. Rumen archaeal abundance at the genus level in Holstein and Jersey steers with different levels of mineral supplementation
Methanogen Genera Holstein Jersey SEM p-value
Con NM HM Con NM HM Breed Treat B × T
Methanobrevibacter 60.88 72.99 57.30 82.04 62.90 58.17 5.843 0.4348 0.1005 0.0651
Methanomassiliicoccus 27.53 15.74 34.57 14.20 33.20 37.75 5.668 0.6271 0.0608 0.0715
Methanosarcina 7.28 7.22 4.80 0.07 0.05 0.00 0.733 < .0001 0.5314 0.5615
Methanobacterium 1.21 1.38 1.10 0.33 0.45 0.45 0.354 0.0348 0.9264 0.9372
Methanosphaera 2.01 2.24 1.68 1.99 1.31 1.31 0.381 0.2366 0.5218 0.5841
Methanimicrococcus 0.67 0.40 0.33 0.61 0.43 0.38 0.172 0.9709 0.2971 0.9475
Methanomicrobium 0.37 0.01 0.20 0.72 1.41 1.87 0.452 0.0373 0.7141 0.5183

Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM.

TMR, total mixed rations; NRC, National Research Council; DM, dry matter.

Download Excel Table
Table 8. Rumen archaeal abundance at the species level in Holstein and Jersey steers with different levels of mineral supplementation
Methanogen species Holstein Jersey SEM p-value
Con NM HM Con NM HM Breed Treat B × T
Methanobrevibacter thaueri 23.69 16.93 18.35 30.15 26.02 26.57 5.594 0.137 0.645 0.976
Methanomassiliicoccus luminyensis 27.53 15.74 34.57 14.20 33.20 37.75 5.668 0.627 0.061 0.072
Methanobrevibacter olleyae 25.32 27.68 23.28 19.43 13.46 14.09 3.495 0.015 0.687 0.619
Methanobrevibacter millerae 9.21 25.18 12.66 24.91 19.54 13.46 5.160 0.508 0.386 0.280
Methanosarcina mazei 7.28 7.22 4.80 0.07 0.05 0.00 0.733 < .001 0.531 0.562
Methanobrevibacter smithii 1.37 2.29 1.92 4.82 2.52 2.57 1.055 0.146 0.734 0.339
Methanobacterium aggregans 1.08 1.02 0.94 0.31 0.45 0.45 0.294 0.041 0.991 0.907
Methanosphaera stadtmanae 2.01 2.24 1.68 1.99 1.31 1.31 0.381 0.237 0.522 0.584
Methanobrevibacter oralis 1.29 0.91 1.10 1.02 0.76 0.45 0.154 0.027 0.100 0.356
Methanimicrococcus blatticola 0.67 0.40 0.33 0.61 0.43 0.38 0.172 0.971 0.297 0.948
Methanomicrobium mobile 0.37 0.01 0.20 0.72 1.41 1.87 0.452 0.037 0.714 0.518
Methanobrevibacter boviskoreani 0.00 0.00 0.00 1.70 0.60 1.02 0.315 0.017 0.543 0.543

Con, only TMR (without mineral supplementation); NM, TMR + NRC recommended concentration of mineral supplementation (Se 0.1 ppm + Zn 30 ppm + Cu 10 ppm)/kg DM; HM, TMR + higher than recommended concentration of mineral supplementation (Se 3.5 ppm + Zn 350 ppm + Cu 28 ppm)/kg DM.

TMR, total mixed rations; NRC, National Research Council; DM, dry matter.

Download Excel Table

DISCUSSION

Dietary supplementation, particularly mineral supplementation, is one of the most important strategies used to reduce the adverse effects of HS in ruminants. However, determining the minerals to use, as well as their most effective concentrations, is a major challenge. Therefore, this study focused on using higher than the recommended concentrations of dietary mineral supplements (Zn, Cu, and Se) during HS. RT is considered a physiological parameter of HS. Although the RT did not vary significantly among different treatment groups in both breeds, the value confirmed that the animals were in HS, which was also observed in a study by Joo et al. [43]. HS reduces DMI, which affects ruminant performance [11,13]. However, similar DMI was observed among the different treatment groups in both breeds, indicating that, in this study, dietary trace minerals had no influence on DMI. These non-significant findings among the different treatment groups were expected because the present study was conducted under HS conditions, with high THI (82.79 ± 1.10); consequently, these findings cannot be compared with the findings expected under normal conditions.

VFA production is negatively correlated with ruminal pH [44]. In this study, the ruminal pH and total VFA were not influenced by trace mineral supplementation; however, lower ruminal pH was recorded in Jersey steers than in Holstein steers, which is in agreement with the findings of Islam et al. [45]. This might be due to the higher total VFA production by Jersey steers compared to that by Holstein steers in this study. Rumen NH3-N concentration depends on the status of dietary protein breakdown in the rumen, rumen microbial utilization, and ruminal epithelial absorption [46,47]. In this study, breed had no influence on the NH3-N concentration. Although the mineral-supplemented groups had higher NH3-N concentrations than the Con group in both breeds, the concentrations were within normal ranges.

The rumen microbiome helps break down the feed substrate in the rumen and improve animal performance. Bacteroidetes and Firmicutes are the major bacterial phyla in different breeds, including Holstein and Jersey breeds [16,27,4852]. Likewise, in the present study, Bacteroidetes, followed by Firmicutes, were the two most abundant bacterial phyla among all the treatment groups in both breeds. Furthermore, a previous study reported that HS alters the rumen microbiota in Holstein and Jersey steers [16]. However, the similar relative abundance of major bacterial phyla and genera among the treatment groups for both breeds suggests that higher concentrations of mineral supplementation did not alter rumen bacterial community composition in the present study. The higher relative abundance of Capnocytophaga in Jersey steers suggests their preferential growth in Jersey steers compared with Holstein steers, which was supported by the findings of Islam et al. [45].

Methanogenesis is the process of CH4 production by methanogens in the rumen via two different pathways: hydrogenotrophic and methylotrophic [6]. Methanobrevibacter is the major archaeal genus involved in the hydrogenotrophic pathway, whereas Methanomassiliicoccus is involved in the methylotrophic pathway [3]. In this study, while the differences were not significant, Methanobrevibacter was the most abundant archaeal genus among all treatment groups in Holstein (collectively around 90%) and Jersey steers (collectively > 95%), followed by Methanomassiliicoccus. Methanosarcina is another identical methanogen that can produce CH4 via three different pathways, namely the hydrogenotrophic, methylotrophic, and acetoclastic pathways [3]. In this study, a significantly higher relative abundance of Methanosarcina was observed in Holstein steers than in Jersey steers, suggesting their preferential growth in the rumen of Holstein steers. The higher relative abundance of the methanogen species Mbr. olleyae was reported in Holstein steers than in Jersey steers. Similarly, King et al. [53] reported a higher relative abundance of Mbr. olleyae in Holstein cows which produced lower CH4 emissions than Jersey cows. They further reported that CH4 production was negatively correlated with the RO group containing Mbr. ruminantium and Mbr. olleyae and was positively correlated with the SGMT group (consisting of Mbr. smithi, Mbr. gottschalkii, Mbr. millerae, and Mbr. thaurei). The significantly lower CH4 production and yield observed in Holstein steers compared with that in Jersey steers in the present study further corroborates these findings. Dietary mineral salt supplementation has been reported to decrease CH4 emissions and methanogen population in lactating cows [7]. However, trace mineral supplementation had no influence on CH4 production, yield, and methanogen abundance regardless of breed, which may be due to the short-term supplementation of dietary trace minerals in this study.

CONCLUSION

Supplementation with high concentrations of dietary organic trace minerals (selenium, zinc, and copper) did not alter enteric CH4 emissions or the methanogenic community. However, Holstein steers emitted low enteric CH4, with a higher relative abundance of Mbr. olleyae than Jersey steers.

Competing interests

No potential conflict of interest relevant to this article was reported.

Funding sources

This research was supported by the Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ015039032022), Rural Development Administration, Korea.

Acknowledgements

Not applicable.

Availability of data and material

Upon reasonable request, the datasets of this study can be obtained from the corresponding author.

Authors’ contributions

Conceptualization: Son AR, Islam M, Kim SH, Lee Sang-Suk.

Data curation: Son AR, Islam M.

Formal analysis: Son AR, Kim SH.

Methodology: Son AR, Islam M, Kim SH, Lee Sung-Sill, Lee Sang-Suk.

Software: Son AR, Islam M.

Validation: Son AR, Islam M, Kim SH, Lee Sang-Suk.

Investigation: Lee Sung-Sill, Lee Sang-Suk.

Writing - original draft: Son AR.

Writing - review & editing: Son AR, Islam M, Kim SH, Lee Sung-Sill, Lee Sang-Suk.

Ethics approval and consent to participate

All animals used in this research were approved by the Sunchon National University (SCNU) Institutional Animal Care and Use Committee (SCNU-IACUC; approval number: SCNU-IACUC-2020-06).

REFERENCES

1.

Johnson KA, Johnson DE. Methane emissions from cattle. J Anim Sci. 1995; 73:2483-92

2.

Appuhamy JADRN, France J, Kebreab E. Models for predicting enteric methane emissions from dairy cows in North America, Europe, and Australia and New Zealand. Glob Chang Biol. 2016; 22:3039-56

3.

Lambie SC, Kelly WJ, Leahy SC, Li D, Reilly K, McAllister TA, et al. The complete genome sequence of the rumen methanogen Methanosarcina barkeri CM1. Stand Genomic Sci. 2015; 10:57

4.

Ellis JL, Dijkstra J, Kebreab E, Bannink A, Odongo NE, McBride BW, et al. Aspects of rumen microbiology central to mechanistic modelling of methane production in cattle. J Agric Sci. 2008; 146:213-33

5.

Patra A, Park T, Kim M, Yu Z. Rumen methanogens and mitigation of methane emission by anti-methanogenic compounds and substances. J Anim Sci Biotechnol. 2017; 8:13

6.

Islam M, Lee SS. Advanced estimation and mitigation strategies: a cumulative approach to enteric methane abatement from ruminants. J Anim Sci Technol. 2019; 61:122-37

7.

Li X, Liu C, Chen Y, Shi R, Cheng Z, Dong H. Effects of mineral salt supplement on enteric methane emissions, ruminal fermentation and methanogen community of lactating cows. Anim Sci J. 2017; 88:1049-57

8.

Silanikove N. Effects of heat stress on the welfare of extensively managed domestic ruminants. Livest Prod Sci. 2000; 67:1-18

9.

St-Pierre NR, Cobanov B, Schnitkey G. Economic losses from heat stress by US livestock industries. J Dairy Sci. 2003; 86:E52-77

10.

Song X, Luo J, Fu D, Zhao X, Bunlue K, Xu Z, et al. Traditional Chinese medicine prescriptions enhance growth performance of heat stressed beef cattle by relieving heat stress responses and increasing apparent nutrient digestibility. Asian-Australas J Anim Sci. 2014; 27:1513-20

11.

Wheelock JB, Rhoads RP, VanBaale MJ, Sanders SR, Baumgard LH. Effects of heat stress on energetic metabolism in lactating Holstein cows. J Dairy Sci. 2010; 93:644-55

12.

Delfino LJB, de Souza BB, da Silva RMN, Silva WW. Effect of heat stress on the eritrogram of ruminants. Agropecu Cient Semiárido. 2012; 8:1-7

13.

Bernabucci U, Biffani S, Buggiotti L, Vitali A, Lacetera N, Nardone A. The effects of heat stress in Italian Holstein dairy cattle. J Dairy Sci. 2014; 97:471-86

14.

Patra AK, Kar I. Heat stress on microbiota composition, barrier integrity, and nutrient transport in gut, production performance, and its amelioration in farm animals. J Anim Sci Technol. 2021; 63:211-47

15.

Kim DH, Kim MH, Kim SB, Son JK, Lee JH, Joo SS, et al. Differential dynamics of the ruminal microbiome of Jersey cows in a heat stress environment. Animals. 2020; 10:1127

16.

Islam M, Kim SH, Son AR, Ramos SC, Jeong CD, Yu Z, et al. Seasonal influence on rumen microbiota, rumen fermentation, and enteric methane emissions of Holstein and Jersey steers under the same total mixed ration. Animals. 2021; 11:1184

17.

Blair R. Nutrition and feeding of organic cattle. Wallingford, OX: CAB International. 2011

18.

Surai PF, Kochish II, Fisinin VI, Juniper DT. Revisiting oxidative stress and the use of organic selenium in dairy cow nutrition. Animals. 2019; 9:462

19.

Surai PF. Selenium in poultry nutrition and health. Wageningen: Wageningen Academic. 2018

20.

Parashuramulu S, Nagalakshmi D, Rao DS, Kumar MK, Swain PS. Effect of zinc supplementation on antioxidant status and immune response in buffalo calves. Anim Nutr Feed Technol. 2015; 15:179-88

21.

Pappas AC, Zoidis E, Chadio SE. Maternal selenium and developmental programming. Antioxidants. 2019; 8:145

22.

Kloubert V, Rink L. Zinc as a micronutrient and its preventive role of oxidative damage in cells. Food Funct. 2015; 6:3195-204

23.

Cortinhas CS, Botaro BG, Sucupira MCA, Renno FP, Santos MV. Antioxidant enzymes and somatic cell count in dairy cows fed with organic source of zinc, copper and selenium. Livest Sci. 2010; 127:84-7

24.

National Research Council. Nutrient requirements of beef cattle. 7th rev. ed Washington, DC: The National Academies Press. 2000

25.

Hernández-Sánchez D, Cervantes-Gómez D, Ramírez-Bribiesca JE, Cobos-Peralta M, Pinto-Ruiz R, Astigarraga L, et al. The influence of copper levels on in vitro ruminal fermentation, bacterial growth and methane production. J Sci Food Agric. 2019; 99:1073-7

26.

Liu C, Li XH, Chen YX, Cheng ZH, Duan QH, Meng QH, et al. Age-related response of rumen microbiota to mineral salt and effects of their interactions on enteric methane emissions in cattle. Microb Ecol. 2017; 73:590-601

27.

Hendawy AO, Sugimura S, Sato K, Mansour MM, Abd El-Aziz AH, Samir H, et al. Effects of selenium supplementation on rumen microbiota, rumen fermentation, and apparent nutrient digestibility of ruminant animals: a review. Fermentation. 2022; 8:4

28.

Bharanidharan R, Arokiyaraj S, Kim EB, Lee CH, Woo YW, Na Y, et al. Ruminal methane emissions, metabolic, and microbial profile of Holstein steers fed forage and concentrate, separately or as a total mixed ration. PLOS ONE. 2018; 13e0202446

29.

AOAC [Association of Official Analytical Chemists] International. Official methods of analysis of the AOAC International. 18th ed Gaithersburg, MD: AOAC International. 2005

30.

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

31.

Van Soest PJ. Collaborative study of acid-detergent fiber and lignin. J Assoc Off Anal Chem. 1973; 56:781-4

32.

Davis MS, Mader TL, Holt SM, Parkhurst AM. Strategies to reduce feedlot cattle heat stress: effects on tympanic temperature. J Anim Sci. 2003; 81:649-61

33.

Hristov AN, Oh J, Giallongo F, Frederick T, Weeks H, Zimmerman PR, et al. The use of an automated system (GreenFeed) to monitor enteric methane and carbon dioxide emissions from ruminant animals. J Vis Exp. 2015; 103e52904

34.

Chaney AL, Marbach EP. Modified reagents for determination of urea and ammonia. Clin Chem. 1962; 8:130-2

35.

Han SK, Kim SH, Shin HS. UASB treatment of wastewater with VFA and alcohol generated during hydrogen fermentation of food waste. Process Biochem. 2005; 40:2897-905

36.

Claassen S, du Toit E, Kaba M, Moodley C, Zar HJ, Nicol MP. A comparison of the efficiency of five different commercial DNA extraction kits for extraction of DNA from faecal samples. J Microbiol Methods. 2013; 94:103-10

37.

Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013; 41e1

38.

Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011; 17:10-2

39.

Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from illumina amplicon data. Nat Methods. 2016; 13:581-3

40.

Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010; 7:335-6

41.

Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10:421

42.

SAS [Statistical Analysis Software]. Statistical Analysis Software. User’s guide statistics version 9.4. Cary. NC: SAS Institute. 2013

43.

Joo SS, Lee SJ, Park DS, Kim DH, Gu BH, Park YJ, et al. Changes in blood metabolites and immune cells in Holstein and Jersey dairy cows by heat stress. Animals. 2021; 11:974

44.

Aikman PC, Henning PH, Humphries DJ, Horn CH. Rumen pH and fermentation characteristics in dairy cows supplemented with Megasphaera elsdenii NCIMB 41125 in early lactation. J Dairy Sci. 2011; 94:2840-9

45.

Islam M, Kim SH, Ramos SC, Mamuad LL, Son AR, Yu Z, et al. Holstein and Jersey steers differ in rumen microbiota and enteric methane emissions even fed the same total mixed ration. Front Microbiol. 2021; 12:601061

46.

Tamminga S. Protein degradation in the forestomachs of ruminants. J Anim Sci. 1979; 49:1615-30

47.

Wallace RJ. Ruminai microbial metabolism of peptides and amino acids. J Nutr. 1996; 126:1326S-34S

48.

Xue M, Sun H, Wu X, Guan LL, Liu J. Assessment of rumen microbiota from a large dairy cattle cohort reveals the pan and core bacteriomes contributing to varied phenotypes. Appl Environ Microbiol. 2018; 84:e00970-18

49.

Difford GF, Plichta DR, Løvendahl P, Lassen J, Noel SJ, Højberg O, et al. Host genetics and the rumen microbiome jointly associate with methane emissions in dairy cows. PLOS Genet. 2018; 14e1007580

50.

Kim M, Yu Z. Variations in 16S rRNA-based microbiome profiling between pyrosequencing runs and between pyrosequencing facilities. J Microbiol. 2014; 52:355-65

51.

Petri RM, Schwaiger T, Penner GB, Beauchemin KA, Forster RJ, McKinnon JJ, et al. Characterization of the core rumen microbiome in cattle during transition from forage to concentrate as well as during and after an acidotic challenge. PLOS ONE. 2013; 8e83424

52.

Jami E, Mizrahi I. Composition and similarity of bovine rumen microbiota across individual animals. PLOS ONE. 2012; 7e33306

53.

King EE, Smith RP, St-Pierre B, Wright ADG. Differences in the rumen methanogen populations of lactating Jersey and Holstein dairy cows under the same diet regimen. Appl Environ Microbiol. 2011; 77:5682-7