Dietary selection of metabolically distinct microorganisms drives hydrogen metabolism in ruminants

The ISME Journal volume 16, pages 2535–2546 (2022)Cite this article

4882 Accesses

7 Citations

21 Altmetric

Metrics details

Ruminants are important for global food security but emit the greenhouse gas methane. Rumen microorganisms break down complex carbohydrates to produce volatile fatty acids and molecular hydrogen. This hydrogen is mainly converted into methane by archaea, but can also be used by hydrogenotrophic acetogenic and respiratory bacteria to produce useful metabolites. A better mechanistic understanding is needed on how dietary carbohydrates influence hydrogen metabolism and methanogenesis. We profiled the composition, metabolic pathways, and activities of rumen microbiota in 24 beef cattle adapted to either fiber-rich or starch-rich diets. The fiber-rich diet selected for fibrolytic bacteria and methanogens resulting in increased fiber utilization, while the starch-rich diet selected for amylolytic bacteria and lactate utilizers, allowing the maintenance of a healthy rumen and decreasing methane production (p < 0.05). Furthermore, the fiber-rich diet enriched for hydrogenotrophic methanogens and acetogens leading to increased electron-bifurcating [FeFe]-hydrogenases, methanogenic [NiFe]- and [Fe]-hydrogenases and acetyl-CoA synthase, with lower dissolved hydrogen (42%, p < 0.001). In contrast, the starch-rich diet enriched for respiratory hydrogenotrophs with greater hydrogen-producing group B [FeFe]-hydrogenases and respiratory group 1d [NiFe]-hydrogenases. Parallel in vitro experiments showed that the fiber-rich selected microbiome enhanced acetate and butyrate production while decreasing methane production (p < 0.05), suggesting that the enriched hydrogenotrophic acetogens converted some hydrogen that would otherwise be used by methanogenesis. These insights into hydrogen metabolism and methanogenesis improve understanding of energy harvesting strategies, healthy rumen maintenance, and methane mitigation in ruminants.

Ruminants have evolved to harbor a complex rumen microbial ecosystem of bacteria, archaea, protozoa, and fungi, which convert complex plant-cell carbohydrates into volatile fatty acids (VFA) and microbial protein usable by the host animal [1, 2]. Molecular hydrogen (H2) is produced during carbohydrate fermentation by anaerobic bacteria, protists, and fungi and is primarily consumed by methane (CH4)-producing archaea in the rumen [3, 4]. This H2 is then consumed by a range of bacteria as an energy source and electron donor, which ensures fermentation remains thermodynamically favorable [5]. Hydrogen cycling is catalyzed by H2-producing and H2-consuming hydrogenases classified into [FeFe]-, [NiFe]-, and [Fe]-hydrogenases based on metal content of their H2-binding sites [6]. Hydrogenases are widespread in bacterial and archaeal genomes, reflecting the metabolic importance of H2 transactions in rumen fermentation [5, 7, 8]. Methane produced from enteric fermentation not only represents a loss of dietary energy but also contributes to global anthropogenic greenhouse gas emissions [9, 10]. Importantly, rumen bacteria also encode hydrogenases and terminal reductases for hydrogenotrophic growth, using electron acceptors such as CO2, fumarate, sulfate, and nitrate. Unlike methanogenesis, some of these alternative hydrogenotrophic pathways produce metabolites usable by the host animal, for example acetate through hydrogenotrophic acetogenesis [5]. It is therefore important to understand how hydrogenotrophic bacteria compete with methanogens for the rumen H2, as this knowledge may facilitate metabolic strategies to enhance energy harvesting with reduced CH4 production in ruminants.

The composition of dietary carbohydrates drives microbial composition and function, including H2 metabolism and methanogenesis [11, 12]. When fed cellulosic plant material, rumen microorganisms hydrolyze and ferment lignocellulose primarily through the acetate pathway releasing H2, resulting in high levels of CH4 production [13]. However, some herbivores such as kangaroos and yak produce less CH4 [14, 15]. Hydrogenotrophic acetogenic bacteria can convert H2 and CO2 into acetate via the Wood–Ljungdahl pathway [5]. Some wild animals in regions where forage quantity or quality is low may have evolved digestive systems with enhanced hydrogenotrophic acetogenesis to promote the efficiency of dietary energy extraction and utilization [16]. It remains to be determined whether diet quality alters the fate of H2 produced, and thus influences the ratio of hydrogenotrophic acetogenesis to methanogenesis and efficiency of energy harvesting in farmed ruminants.

Despite being widely used, starch-rich high-grain diets present risks for ruminant health. High starch intake stimulates rapid VFA and lactate production in non-adapted rumens, leading to a rapid and sometimes pronounced drop in rumen pH [17]. Even in animals adapted to grain, a prolonged acidic rumen, a condition termed subacute rumen acidosis [18], results in reduced feed intake, chronic inflammation, rumen papillae dysfunction, and other diseases [19]. Unlike acetate production, lactate production is not associated with H2 production and facilitates reoxidation of redox cofactors (e.g., NADH to NAD+) [20]. However, given lactate has a lower pKa than VFAs, its accumulation decreases rumen pH and hence rapid conversion of lactate to propionate and other VFAs by lactate-utilizing bacteria is required to maintain a normal rumen function. In order to promote a healthy rumen while keeping CH4 production at low levels, it is important to understand how the conversion of lactate to propionate is controlled.

Altogether, previous studies suggest that rumen microbiota exhibits distinct metabolic strategies to adapt to fiber-rich or starch-rich diets. However, we lack a detailed mechanistic understanding of the microbiota and pathways of H2 metabolism affected by dietary shifts. To address this knowledge gap, we conducted one in vivo and two in vitro experiments, and developed a rumen model of fiber-rich and starch-rich selected microbiomes. We observed that feeding fiber-rich versus starch-rich diets selected for distinct rumen microbiome compositions, capabilities, and activities. Using metagenomic profiling, we found that microbiota resulting from feeding these two diets exhibited distinct carbohydrate and H2 metabolic pathways, with differences in downstream methanogenesis and VFA production pathways. Results also suggest a potential role in the rumen of hydrogenotrophic acetogenesis as a H2 sink other than methanogenesis with the fiber-rich dietary treatment, which might compete for H2 available for methanogenesis and produce extra energy available for the host animal. This study emphasizes that dietary interventions modulate rumen microbiota composition and in turn carbohydrate degradation, H2 metabolism, and methanogenesis.

All animals involved in the experiment were cared for in accordance with the Animal Care and Use Guidelines of the Animal Care Committee, Institute of Subtropical Agriculture, the Chinese Academy of Sciences, Changsha, China, with all animal experimental procedures approved by the committee (approval number ISA-W-201901).

A total of 24 local breed (Xiangxi) beef cattle (initial body weight 147 ± 9.8 kg) were randomly assigned to one of two dietary treatments that lasted for 300 d with three periods (Fig. 1A). The fiber-rich diet was kept constant in all three experimental periods and was formulated to have a forage neutral detergent fiber (NDF) content greater than 40% dry matter (DM) basis by including corn stover and rice straw. Three starch-rich diets with increasing starch content were formulated by gradually replacing one third of corn stover with maize meal (DM basis) in each experimental period, until reaching 90% of concentrate (DM basis) in period 3 (Supplemental Table S1). All the cattle were fed twice daily (7 a.m. and 5 p.m.) and had free access to drinking water.

A Overview of daily fiber and starch intake (kg/d) and rumen pH over three experimental periods; B feed digestibility, dissolved hydrogen (H2), lactate and volatile fatty acids profile in vivo; C feed degradation, methane (CH4) production and volatile fatty acids profile in vitro rumen experiment 1; D inspection of morphological changes in rumen papillae; E ratio of emptied rumen to carcass and papillae height; F q-PCR results of genes related to volatile fatty acids absorption and intracellular pH regulation in ruminal epithelium. OM organic matter, NDF neutral detergent fiber, VFA volatile fatty acid, Ace acetate, Pro propionate, But butyrate; others, valerate, isovalerate and isobutyrate; HMGCL 3-hydroxy-3-methylglutaryl-CoA lyase; HMGCS-1 3-hydroxy-3-methylglutaryl-CoA synthase, isoform 1, HMGCS-2 3-hydroxy-3-methylglutaryl-CoA synthase, isoform 2, NHE-1 Na+/H+ exchanger 1, NHE-2 Na+/H+ exchanger 2, NHE-3 Na+/H+ exchanger 3, MCT-1 monocarboxylate transporter, isoform 1. Data with error bars are expressed as mean ± standard error. *p < 0.05, **p < 0.01, ***p < 0.001, n = 12/group.

Full details are provided in Supplementary Materials and Methods. Briefly, nutrient digestibility was measured between days 93 and 98 of each of the three periods. Rumen contents were collected at 0, 2.5 and 6 h after the morning feeding on days 99 and 100 of each of the three periods. Rumen pH and dissolved hydrogen concentrations were immediately measured after collecting rumen contents, while other subsamples were frozen for DNA extraction and fermentation product measurements. All the cattle were slaughtered at the end of experiment after fasting for 12 h, and rumen contents and epithelial samples were collected.

Detailed information is provided in Supplementary Materials and Methods. The present study included two batch culture ruminal in vitro experiments conducted at final experimental period. In vitro experiment 1 compared ruminal fermentation between fiber-rich and starch-rich diets as incubation substrates, with each substrate being incubated with the rumen inoculum from animals fed the corresponding diet. In vitro experiment 2 compared the activity of fiber-rich and starch-rich selected microbiomes by inoculating rice straw with rumen inoculum from animals fed the fiber-rich and starch-rich diets, respectively. In vitro ruminal fermentation was performed according to the procedure of Wang et al. [21].

All detailed analysis procedures are provided in Supplementary Materials and methods. Feed and feces composition were analyzed using the methods of the Association of Official Analytical Chemists [22], Van Soest et al. [23] and Karthner and Theurer [24]. Dissolved H2 and CH4 were extracted from the liquid phase of rumen contents into the gas phase, measured using a gas chromatograph (GC, Agilent 7890A, Agilent Inc., Palo Alto, CA), and concentration calculated using the equation of Wang et al. [25]. Individual volatile fatty acid (VFA) concentrations were analyzed by GC (Agilent 7890A, Agilent Inc., Palo Alto, CA) [25], while the concentration of lactate was analyzed by high-performance liquid chromatography (HPLC, Agilent LC1290, Agilent Inc., Palo Alto, CA) [26].

DNA was extracted from rumen samples using repeated bead beating and column filtration as described [27]. Total RNA of the ruminal epithelium was extracted from the tissue after removing the genomic DNA. Quantitative real-time PCR (q-PCR) was performed to measure relative gene expression of seven genes in ruminal epithelium (HMGCL, HMGCS-1, HMGCS-2, NHE-1, NHE-2, NHE-3, MCT-1), the expression of which were normalized to the reference genes (GAPDH and β-actin) using a 2–ΔΔCt method [28]. Absolute quantification of specific microbial taxa was performed according to the procedures described by Ma et al. [29].

All bioinformatics software, parameters, and databases involved in the study are present in Supplementary Materials and Methods. Briefly, the V3 and V4 region was used for 16S rRNA gene sequencing, with all amplicon libraries preparation and sequencing performed on a MiSeq platform (Illumina, San Diego, CA, USA). A total of 1,238,536 high-quality reads were generated with an average of 51,606 ± 1737 reads per sample for assignment of amplicon sequence variants (ASVs) by DADA2. Taxonomy was annotated using SILVA (release 138, http://www.arb-silva.de). For shotgun metagenome sequencing, 1 μg of genomic DNA was sheared by a Covaris S220 Focused-ultrasonicator (Woburn, MA USA) and sequencing libraries were prepared with a length of approximately 350 bp (ranging from 300 to 400 bp). Sequencing of metagenomic libraries was performed on the HiSeq X platform (Illumina, San Diego, CA, USA) in paired-end 150 bp (PE150) mode. After quality control, assembly, prediction, and clustering into a non-redundant gene set, the annotations of the gene set were performed with KEGG, CAZy, HydDB, and NCBI-NR, and the gene abundance within total samples were calculated. Microbial genomes were further verified by using Hungate1000 collection, which includes virtually all the bacterial and archaeal species that have been cultivated from the rumen of a diverse group of animals [30]. The abundance of each genome was calculated by metawrap quant_bins module with default parameters [31].

The metabolic data were analyzed by one-way analysis of variance (ANOVA) in SPSS 21.0 software (SPSS Inc., Chicago, IL). When sampling time was included, data were analyzed using a linear mixed model with treatment and treatment by sampling time interaction as fixed effects, and sampling time as a repeated measures variable. The microbial community composition data was analyzed using Wilcoxon rank-sum test in the JMP Pro software (JMP Pro version 13.2.1, SAS Institute Inc. SAS Institute, Cary, NC, USA). The genomic ranks of attributes were evaluated by Correlation, ReliefF, Symmetrical Uncert, and multi-cluster feature selection (MCFS) methods in software of Waikato Environment for Knowledge Analysis (WEKA) (version 3.8.4, Hamilton, New Zealand) [32], and further comprehensively analyzed by RobustRankAggreg (RRA) R package [33]. All p values were adjusted for False Discovery Rate (FDR) using the Benjamini-Hochberg method. Statistical significance was declared at p ≤ 0.05 and tendencies at 0.05 < p ≤ 0.10.

Animals were adapted to a starch-rich diet by gradually increasing dietary concentrate content from 50 to 90% over the three 100-d experimental periods (Fig. 1A). With both diets, rumen structure and epithelial morphology were robust and rumen pH remained above 6.0 (Fig. 1A, D Additional Table S4), indicating a healthy rumen function throughout the entire experiment.

Cattle with high fiber intake exhibited greater fiber digestibility, while cattle adapted to the starch-rich diet had greater organic matter and starch digestibility but lower fiber digestibility in vivo and in the in vitro experiment 1 (Fig. 1B, C, Additional Table S4, p < 0.001). Cattle adapted to the starch-rich diet exhibited greater rates of carbohydrate degradation, leading to increased total VFA concentration in periods 1 and 2 in vivo (Additional Table S4), and in the in vitro experiment 1 (Fig. 1C, p < 0.001). The lower VFA concentration (p < 0.01) observed in vivo at period 3 with the starch-rich diet could have been the result of a greater increase in the rate of VFA absorption than production, as suggested by the increased copy number of some of the genes related to VFA absorption and intracellular pH regulation in the ruminal epithelium (Fig. 1F), namely the limiting enzyme in ketone formation 3-hydroxy-3-methylglutaryl-CoA lyase (HMGCL), 3-hydroxy-3-methylglutaryl-CoA synthase isoform 2 (HMGCS-2), Na+/H+ exchanger 2 (NHE-2), and Na+/H+ exchanger 3 (NHE-3) [34,35,36]. Changes in VFA absorption between diets are also suggested by the alternation of rumen organ structure, with a lower ratio of rumen organ to carcass and shorter rumen papilla in the starch-rich diet (Fig. 1D, E).

Acetate molar percentage and the acetate to propionate molar ratio were higher in cattle adapted to the fiber-rich diet (p < 0.001). Such differences in VFA profiles agree with the in vitro ruminal fermentation results (Fig. 1B, C, Additional Table S4). The starch-rich diet increased (p < 0.001) ruminal dissolved H2 (dH2, 2.5-fold higher) and reduced CH4 production (p = 0.003). Altogether, these results suggest that the two diets differentially influence VFA production, H2 metabolism, and methanogenesis, likely due to the promotion of distinct microbial activities.

16S rRNA gene sequencing was used to examine rumen microbial community composition. Across the 24 samples, 3028 amplicon sequence variants (ASVs) were observed. Microbial composition clearly clustered by diet, based on unweighted UniFrac dissimilarity (Fig. 2A), indicating that carbohydrate type drove marked changes in the microbial community composition. In addition, the random forest was further applied to classify groups and selected the 20 most important bacterial genera ranked by mean decrease accuracy of classification of bacterial communities by diet (Fig. 2B, Additional File S1). Among them, Ruminobacter and Succinivibronaceae UCG-002 were the two most representative genera for the starch-rich diet, whilst Fibrobacter was the most representative genus in case of the fiber-rich diet. Furthermore, correlation network analysis of these representative bacterial genera revealed strong positive correlations within each treatment, and negative correlations between the treatments (p < 0.05, Spearman's | r | > 0.5, Additional Fig. S3).

A principal coordinate analysis (PCoA) profile of ruminal bacterial community based on unweighted UniFrac dissimilarity matrix at the taxa (ASVs) level (PERMANOVA, p = 0.001, R2 = 0.19); B random forest analysis of rumen bacterial community. The y-axis, from top to bottom, displays the genera ranked by their relative importance based on mean decrease accuracy in the classification of groups; C relative abundance of the top five phyla; D relative abundance of genus Fibrobacter, Ruminobacter, and Treponema; E the relationships of feed digestibility and volatile fatty acids with bacterial community at genus level (average percent > 0.5%). No rank means there is no specific taxonomic information at the genus level. Only Spearman's significance levels p < 0.05 are shown, yellow and blue indicate negative and positive correction respectively. OM organic matter; NDF neutral detergent fiber. Significance was tested using independent two-group Wilcoxon rank-sum tests. Data with error bars are expressed as mean ± standard error. **p < 0.01, ***p < 0.001, n = 12/group.

At the phylum level, Bacteroidota (Bacteroidetes) and Fibrobacterota (Fibrobacteres) were significantly enriched with the fiber-rich diet, whereas Proteobacteria and Spirochaetota (Spirochaetes) were more abundant with the starch-rich diet (Fig. 2B, p < 0.001). It should be noted that, due to the specificity of the primer pairs used (targeting the V3–V4 primer), Spirochaetota and archaea are potentially underrepresented in this analysis [37, 38]. At the genus level, Fibrobacter (37.3-fold higher), Christensenellaceae R-7 group (2.6-fold higher), and Bacteroidales RF16 group (4.0-fold higher) were enriched with the fiber-rich diet, while Ruminobacter (12.0-fold higher), Treponema (6.9-fold higher), and Succinivibrionaceae UCG-002 (21.8-fold higher) were enriched with the starch-rich diet (Fig. 2D, Additional Fig. S2, p < 0.001). Four of these genera were also among those with the highest discriminatory power based on the random forest analysis (Fig. 2B). In addition, Fibrobacter, Bacteroidales RF16 group, and Christensenellaceae R-7 group were positively correlated with acetate concentration and neutral detergent fiber (NDF) digestibility, while Ruminobacter, Succinivibrionaceae UCG-002, and Treponema were positively correlated with propionate concentration and starch digestibility (Fig. 2E). Differences in microbial composition are thus associated with distinct substrate preferences and degradation abilities.

After removing the reads assigned to the host, we obtained 503 Gb of paired-end sequencing data, which averaged 21.0 Gb (ranging from 17.4 to 29.1 Gb) per sample. We first investigated the metabolic capacities of the rumen microbiomes by using the Kyoto Encyclopedia of Genes and Genomes (KEGG) [39]. There were 71.4 and 33.9% genes assigned to KEGG Orthology (KO) database and KEGG pathways, respectively. Principal coordinate analysis (PCoA) of all KO genes showed that the fiber-rich and starch-rich diets selected for different metabolic functions (Fig. 3A), with "carbohydrate metabolism pathways" being the most abundant category that was significantly different between diets (Additional Fig. S4).

A PCoA profile of all KO genes (PERMANOVA, P < 0.001, R2 = 0.70); B relative abundance of total of CAZymes genes; C PCoA profile of the GH family (PERMANOVA, P < 0.001, R2 = 0.64); D gene abundance of the top five GH family enzymes, and their phylogenetic distribution enriched per phylum (AA auxiliary activity, CBM carbohydrate-binding module, CE carbohydrate esterase, GH glycoside hydrolase, GT glycosyltransferase, PL polysaccharide lyase); E acetate, butyrate, propionate and methanogenesis pathways of KEGG enzymes expressed as ratios of alignments of fiber-rich versus starch-rich treatment (log2 ratio); and pie charts show the phylogenetic distribution of the pathways enriched in each treatment at phylum. Ace-P acetate production pathway, Pyr-But-P pyruvate to butyrate production pathway, Ace-But-P acetate to butyrate production pathway, Pro-Lac-P propionate (lactate) pathway, Pro-Suc-P propionate (succinate) pathway, Met methanogenesis pathway. All KO genes in enriched pathways were assigned to the identified phylum. Details of all KO genes together with identified genus and pathways are in additional files S4 and S5. Significance was tested using independent two-group Wilcoxon rank-sum tests. Data with error bars are expressed as mean ± standard error. Asterisks denote significant adjusted p values: **p < 0.01, ***p < 0.001, n = 12/group.

We screened for carbohydrate-active enzymes (CAZymes) in the assembled metagenomic contigs (Additional File S2). High starch intake increased the abundance of total CAZymes (p = 0.01), including auxiliary activities (AAs), carbohydrate-binding module (CBMs), and glycosyltransferases (GTs) (Fig. 3B). Furthermore, the rank of the top five carbohydrate-active enzyme classes and top ten assigned phyla or genera differed between the two treatments (Additional Fig. S5). The GHs had the highest relative abundance among the six CAZyme classes, and the abundance of GH subfamilies was significantly different between the two dietary treatments (Fig. 3C). The fiber-rich diet selected for a greater abundance of β-xylosidase GH43 (p = 0.001, 1.3-fold higher), which was assigned to Prevotella, Bacteroides, and Fibrobacter at the genus level (Additional Fig. S6). Abundance of α-amylase GH13 was greater with the starch-rich diet (Fig. 3D, Additional File S3, p < 0.001, 1.4-fold higher) and was assigned to Prevotella, Ruminobacter, and Clostridium (Additional Fig. S6). Thus, both dietary treatments promoted bacterial communities with distinct capacities to digest carbohydrates.

We then analyzed the abundance of genes encoding for enzymes related to the pathways of acetate, propionate, butyrate, and methanogenesis (Additional File S4). The fiber-rich diet enriched for the pathways of acetate production (pta gene), acetate to butyrate production (acs gene), propionate production via succinate as intermediate (MUT and sdhA genes), and methanogenesis (mcrA and mcrB genes), while the starch-rich diet enriched for the acrylate pathway of propionate production (ldh and pct genes) (Fig. 3E, Additional Fig. S7, p < 0.05). Bacteroidota, Firmicutes, and Proteobacteria were the major phyla assigned to the acrylate pathway in the starch-rich diet. With the fiber-rich diet, Fibrobacterota was prominent in fermentative acetate and propionate production, and as expected, Firmicutes was the major phylum for acetate to butyrate production pathway, and Methanobacteriota (Euryarchaeota) was the assigned phylum for methanogenesis (Fig. 3E, Additional File S5).

Given the significant differences observed in VFA profiles, dH2 concentrations, and CH4 production (Fig. 1B, C), we screened for the genes encoding for the catalytic subunits of H2-producing and H2-consuming enzymes in the assembled contigs. Of the 2,686 genes identified, 82%, 17% and 1.2% were annotated as [FeFe]-, [NiFe]- and [Fe]-hydrogenases respectively. Hydrogenases were taxonomically assigned to 149 genera, including Bacteroides, Clostridium, Oscillibacter, Methanobrevibacter, and Ruminococcus (Additional File S6). As the two diets exhibited distinct hydrogenase composition (Additional Fig. S8), we then classified hydrogenases into subgroups. The trimeric group A3 [FeFe]-hydrogenases, which mediate the process of electron confurcation during fermentative carbohydrate degradation leading to H2 production [40], were the most abundant hydrogenases with both treatments (Fig. 4A). The diaphorase subunit (HydB) of these hydrogenases was highly enriched with the fiber-rich diet (p < 0.01, 1.8-fold higher), and was primarily encoded by Firmicutes and Bacteroidota at the phylum level (Additional File S6).

A Hydrogenase genes distributions assigned by phylum; B genes of associated terminal reductases distributions assigned by phylum. Fermentative hydrogenases (group B, A1 and A2 FeFe-hydrogenases), electron-bifurcating hydrogenases (group A3 and A4 FeFe-hydrogenases), energy-converting hydrogenases (bidirectional; group 4a, 4c, 4d, 4e, 4f and 4g NiFe-hydrogenases), methanogenic hydrogenases (Fe-hydrogenases, group 3a, 3c, 4h, 4i and 1k NiFe-hydrogenases), respiratory hydrogenases (group 1a, 1b, 1c, 1d, 1k, 2a), sensory hydrogenases (group C FeFe-hydrogenases). HydB hydrogenase-associated diaphorase. NifH nitrogenase. H2 uptake pathways can be coupled to fumarate reduction (FrdA fumarate reductase), nitrate ammonification (NrfA, ammonia-forming nitrite reductase; NarG, dissimilatory nitrate reductase; NapA, periplasmic nitrate reductase), sulfate and sulfite reduction (AprA adenylylsulfate reductase; AsrA alternative sulfite reductase; DsrA, dissimilatory sulfite reductase), dimethyl sulfoxide and trimethylamine N-oxide reduction (DmsA DMSO and TMAO reductase), reductive acetogenesis (AcsB, acetyl-CoA synthase), aerobic respiration (CydA cytochrome bd oxidase), and methanogenesis (McrA methyl-CoM reductase). No rank means there is no specific taxonomic information at the phylum level. Only genes of average relative abundance > 1 were shown, while other were shown in additional file S6. Significance was tested using independent two-group Wilcoxon rank-sum tests. Data with error bars are expressed as mean ± standard error. **p < 0.01, ***p < 0.001, n = 12/group.

The fiber-rich diet enriched genes for methanogenic hydrogenases (group 3a [NiFe]-hydrogenases and [Fe]-hydrogenases) encoded by Methanobacteriota (Fig. 4A, p < 0.05, 1.5-fold higher). The starch-rich diet enriched genes encoding for fermentative group B [FeFe]-hydrogenases (p < 0.01, 2.2-fold higher), energy-converting group 4e and 4g [NiFe]-hydrogenases, respiratory group 1d [NiFe]-hydrogenases, and sensory group C3 [FeFe]-hydrogenases (Fig. 4A, p < 0.05), which were taxonomically assigned mainly to Firmicutes and Spirochaetota (Additional File S6). These results suggest the starch-rich diet resulted in enhanced H2 flow, which was consistent with greater fermentation and higher rumen dH2 concentration (Fig. 1B).

We further analyzed signature genes that support hydrogenotrophic growth, including methanogenesis, acetogenesis, fumarate reduction, sulfidogenesis, nitrate reduction, and aerobic respiration [5]. The fiber-rich diet selected for signature genes associated with methanogenesis, including methyl-CoM reductase (mcrA, Fig. 4B, p = 0.009, 1.6-fole higher), which were mainly affiliated with Methanobrevibacter (Additional File S6). The starch-rich diet selected for a higher abundance of group 1d [NiFe]-hydrogenases (p < 0.05) encoded by Firmicutes, a membrane-bound type of enzymes known to support hydrogenotrophic respiration, along with the signature gene for sulfate reduction (aprA, p = 0.002, 3.0-fold higher). An unexpected result was that the marker gene for hydrogenotrophic acetogenesis (acetyl-CoA synthase, acsB, p < 0.01) was threefold more abundant with the fiber-rich diet (Fig. 4B). Acetyl-CoA synthase is a reversible enzyme, for example mediating acetate utilization in various sulfate-reducing bacteria [41], but most reads (about 75%) were most closely related to those of acetogenic Firmicutes.

The Hungate1000 collection of assembled genomes, which includes virtually all of the bacterial and archaeal species that have been cultivated from the rumen of diverse ruminant species [30, 42], was used to gain a stronger strain-level understanding of the findings of this study. One million random reads from each sample were extracted and aligned to the Hungate1000 genomes, resulting in an average mapping rate of 20.2% (Additional Fig. S9). We found that 257 and 168 genomes were more enriched by the fiber-rich and starch-rich diets, respectively (Additional Fig. S10). Genomes encoding for GH43 were more enriched with the fiber-rich diet and assigned to Bacteroides, Prevotella, and Fibrobacter, while genomes encoding for GH13 were enriched with the starch-rich diet and assigned to Ruminobacter and Clostridium. The fiber-rich diet selected 8 out of 10 putatively methanogenic genomes, which harbor unique hydrogenases (group 3a, 3c, 4h, 4i [NiFe]-hydrogenases and [Fe]-hydrogenases) and the signature gene mcrA, as well as putatively 5 out of 8 hydrogenotrophic acetogenic bacterial genomes encoding electron-bifurcating hydrogenases (group A3, A4 [FeFe]-hydrogenases) or the signature gene acsB. The starch-rich diet selected for 67 out of 130 genomes encoding fermentative hydrogenases (group A1, A2, and B [FeFe]-hydrogenases), assigned to Lachnospirales and Selenomonadales, and 17 out of 33 genomes encoding for respiratory hydrogenases (group 1 and 2a [NiFe]-hydrogenases) assigned to Selenomonadales. These genomes with terminal sulfate and sulfite reductases were more enriched with the starch-rich diet, and mainly assigned to Lachnospirales (Additional File S7).

We performed an analysis of ranks of attributes based on differential abundance between two treatments. Fibrobacter succinogenes subsp. elongatus strain HM2 and Lachnospiraceae bacterium AC2028, both of which have functions in fiber degradation [30], were the most representative genomes enriched with the fiber-rich diet (p < 0.001). Both Ruminobacter sp. RM87 and Succinimonas amylolytica DSM 2873, known to degrade starch [30], were the most representative genomes with the starch-rich diet (Fig. 5A, p < 0.001, Additional Fig. S11). We validated by q-PCR that 16S rRNA gene copies of F. succinogenes and Ruminobacter amylophilus were the more abundant in the fiber-rich and starch-rich treatments, respectively (Additional Fig. S12, p < 0.001). Enhanced fiber degradation was further verified by in vitro experiment 2, which indicated that inoculating microbiome of fiber-rich diet resulted in greater NDF degradation (Fig. 5B, p = 0.026).

A metabolic features of representative microorganisms. Reads from each sample were aligned to sequenced genomes of cultured rumen microorganisms in the Hungate1000 collection using the burrows-wheeler alignment tool. The hydrogenase function, substrate utilization, and metabolite production of each microorganism based on the known growth characteristics [5, 30] are colored in blue, orange, and green, respectively. The ratios between alignments of fiber-rich/starch-rich samples to each genome are presented, and data were expressed as log2 (ratio). Genomes attributes: fibrolytic bacteria, Fibrobacter succinogenes subsp. elongatus strain HM2 and Lachnospiraceae bacterium AC2028; amylolytic bacteria, Ruminobacter sp. RM87 and Succinimonas amylolytica DSM 2873; acetate producer, Lachnospiraceae bacterium G41 and Fibrobacter succinogenes subsp elongatus strain HM2; butyrate producer, Lachnospiraceae bacterium AC2028 and Butyrivibrio sp. LB2008; lactate utilizer and propionate producer, Megasphaera elsdenii strain J1 and Anaerovibrio lipolyticus LB2005; hydrogenotrophic methanogen, Methanobrevibacter thaueri strains DSM 11995 and sp. YE315; acetogens, Acetitomaculum ruminis DSM 5522. Detailed information of other genomes is shown in additional file S7; B methane (CH4) production, neutral detergent fiber (NDF) degradation, pH and volatile fatty acid (VFA) profile of in vitro rumen experiment 2 with rice straw fermentation by inoculating fiber-rich or starch-rich selected microbiome. Data with error bars are expressed as mean ± standard error. **p < 0.01, ***p < 0.001, n = 12/group.

We selected the two most differential abundant genomes of microorganisms inferred to produce acetate, propionate, or butyrate (Additional File S7). The microbiome of the fiber-rich diet favored acetate and butyrate production while being enriched with Lachnospiraceae bacterium G41 and AC2028, F. succinogenes subsp. elongatus strain HM2, and Butyrivibrio sp. LB2008 (p < 0.01), while the microbiome of the starch-rich diet favored propionate production with lactate as an intermediate along with enrichment of Megasphaera elsdenii strain J1 and Anaerovibrio lipolyticus LB2005 (p < 0.001, Additional Fig. S11). These results were further verified by in vitro experiment 2 by incubating rice straw inoculated with fiber-rich or starch-rich selected microbiomes (Fig. 5B). Accordingly, five of six representative genomes encoded fermentative, electron-bifurcating and energy-converting hydrogenases, while one genome (i.e., Megasphaera elsdenii strain J1) had the capacity to use H2 as a substrate for VFA production (Fig. 5A). These results further upheld that dietary carbohydrate selected for microbiomes with distinct pathways of VFA production and H2 metabolism.

We selected genomes encoding marker genes for hydrogenotrophic methanogenesis (mcrA) and acetogenesis (acsB). The most differentially abundant genomes were Methanobrevibacter thaueri strains DSM 11995 and sp. YE315, as well as the hydrogenotrophic acetogen Acetitomaculum ruminis DSM 5522 [43], which were all enriched with the fiber-rich diet (Fig. 5A and Fig. S11, p < 0.01, Additional File S7). Both M. thaueri strains DSM 11995 and sp. YE315 encode methanogenic hydrogenases to use H2 to reduce CO2 to CH4, while A. ruminis DSM 5522 encodes electron-bifurcating group A3 [FeFe]-hydrogenases predicted to use H2 to reduce CO2 to acetate (Fig. 5A). The acsB gene is also encoded in some butyrate-producing bacteria (e.g., Eubacterium limosum), which produces acetate through the Wood–Ljungdahl pathway and converts it to butyrate [44, 45]. In turn, in vitro experiment 2 where rice straw was the sole substrate indicated that inoculation of the fiber-rich selected microbiome resulted in lower CH4 and greater acetate and butyrate production compared to the starch-rich selected microbiome (Fig. 5B, p < 0.05). This suggests the possibility that the fiber-rich diet may have selected for hydrogenotrophic acetogens, which may use some of the H2 that would otherwise be used for methanogenesis to produce acetate.

The fiber-rich diet selected for fibrolytic bacteria (e.g., Fibrobacter and Ruminococcus) and enriched for genes encoding for xylosidases (e.g., GH43, Fig. 6A). The characteristic microbiome and metabolic strategies promoted by the fiber-rich diet are likely to be important for ruminants living in harsh environmental conditions, allowing them to digest and ferment low-quality plant biomass to produce VFA to be used by the host animal for maintenance, reproduction, and growth. Tibetan ruminants, which are well adapted to digest poor-quality forage diet [15, 46], have been shown to harbor more fibrolytic bacteria than lowland ruminants, especially of Ruminococcus albus (100-fold higher) and F. succinogenes (50-fold higher) [47]. The camel, living in the barren desert steppe with access to the lignified coarsest parts of scarce plants, has likewise greater enrichment of fiber-degrading rumen bacteria compared to domestic ruminants such as cattle and sheep [48, 49]. Furthermore, fibrolytic bacteria can also be enriched in the hindgut of non-ruminants including pigs [50], horses [51], and rabbits [52] when the fiber content of the diet is increased. Proliferation of fibrolytic bacteria for enhanced fiber utilization can be an important metabolic strategy for host animals to adapt to eating lignified diets. However, enrichment in fibrolytic bacteria with the fiber-rich treatment could be associated with increased CH4 production, leading to a reduction in dietary energy use efficiency.

A consolidating analysis on enzymes, microorganisms of Hungate1000 collection, metabolites and pathways recruitment. "+ H2", H2 produced; "- H2", H2 consumed. The fold changes (left: fiber-rich/starch-rich; right: starch-rich/fiber-rich) of enzymes or microorganisms in each group are labeled besides, and presented in additional file S8 in details; B proposed rumen H2 metabolism models. Representative species or genus and corresponding hydrogenases are shown next to the arrow, and extracted from Figs. 1B, C, 4A, 5, and additional file S8. The solid and dotted line (square) indicated trends of enhanced and reduced pathways, respectively. H2 hydrogen. The upward arrows, increased metabolites; the downward arrows, decreased metabolites.

The first important finding of this work is that a fiber-rich diet increased the abundance of the acsB gene, encoding for acetyl-CoA synthase of hydrogenotrophic acetogens (e.g., Acetitomaculum spp.), compared to a starch-rich diet. Hydrogenotrophic acetogens have been isolated from the rumen and foregut microbial ecosystems of cattle and kangaroos fed with high forage diets and from the hindgut of wood-eating termites [14, 43, 53]. Some ruminal hydrogenotrophic acetogens (e.g., Acetitomaculum) have the ability to simultaneously use inorganic (H2/CO2) and organic substrates (e.g., cellobiose and glucose) as energy and carbon sources, whereas others (e.g., Blautia sp. Ser8) do not use H2 in the presence of glucose, producing acetate as the sole reduced end product [16, 54]. One discrepancy with other studies is that we observed greater abundance of Spirochaetota in the starch-rich diet, whereas Spirochaetota have been shown to be stimulated by fiber-rich diets and can be hydrogenotrophic acetogens in other animal hosts [55,56,57]. It is unclear if our findings are compounded by the specificity of the primers used to profile community composition. Nevertheless, the metagenomic and activity-based analyses suggest an increase in the abundance of acetogenic bacteria, in addition to fermentative acetate production, can also contribute to increased acetate production with the fiber-rich diet. As some butyrate-producing bacteria (e.g., Eubacterium limosum) encode the acsB gene for acetyl-CoA synthase, a greater proportion of butyrate produced by fiber-rich selected microbiome can go through the acetate pool.

In turn, this work provides novel insights about the possible role of hydrogenotrophic acetogenesis in the rumen microbial communities selected by the fiber-rich diet. Studies report that methanogens can outcompete for acetogens due to their higher affinity for H2 and lower Gibbs energy change of H2 utilization for methanogenesis compared to hydrogenotrophic acetogenesis [7, 58]. The classical theory holds that stimulation of acetogens may be an effective strategy for redirection of metabolic hydrogen towards acetate production in methanogenesis-inhibited scenarios with increased ruminal H2 concentration [5]. However, we found an increased abundance of the marker gene for hydrogenotrophic acetogenesis (e.g., acsB) with lower dH2 concentration in the fiber-rich diet. Our in vitro experiment 2 with rice straw as the sole substrate indicated that inoculation with fiber-rich selected microbiome paradoxically promoted acetate production with decreased CH4 production in comparison to the starch-rich selected microbiome. These results suggest that a greater community of hydrogenotrophic acetogens in the fiber-rich selected microbiome might perhaps compete for H2 with methanogens. There has been speculation that high-altitude ruminants yak and Tibetan sheep might favor acetate production with lower CH4 formation compared to lowland cattle and sheep [15, 59]. Molecular hydrogen incorporation into hydrogenotrophic acetogenesis is estimated to contribute for about 1% of acetate formed and account for 1 to 2% of total reducing equivalents incorporated in in vitro cultures inoculated with sheep rumen fluid [60]. The stimulation of acetate synthesized by acetogens would not only decrease the emissions of greenhouse gas CH4, but would also help improving energy use efficiency for host animals. Such pathways of energy harvesting may be a competitive strategy for ruminants fed poor-quality diets.

High starch intake generally results in decreased pH, and even subclinical or lactic rumen acidosis [17]. We found that gradually increasing dietary starch up to 90% did not drop the pH below 6.0 or harm the rumen epithelium, although lactate concentration significantly increased. The enrichment of M. elsdenii and A. lipolyticus, which use the acrylate pathway to convert lactate to propionate, is likely important for avoiding acidosis. M. elsdenii reportedly utilizes 60–80% of the lactate produced in the rumen [61], and accounts for 20% of lactate utilizers in the rumen of animals fed high-concentrate diets [62]. A. lipolyticus can also utilize lactate and produce propionate as an end product [63, 64]. In addition to preventing lactate accumulation [17, 65], lactate conversion to propionate incorporates reductant into NADH, thus competing with H2 formation and ultimately methanogenesis [66]. A strong community of bacteria metabolizing lactate to propionate via acrylate can prevent lactate accumulation and help maintaining a healthy rumen.

A second important finding of this work is that, in both fiber-rich and starch-rich selected microbiomes, electron-bifurcating group A3 [FeFe]-hydrogenases were inferred to be the primary mediators of ruminal H2 production. This result coincides with the findings about the importance of H2 formation by electron confurcation in the rumen [67]. Equally important is the differential abundance of hydrogenases and its relationship with methanogenesis promoted by the two diets. The fiber-rich diet resulted in low rumen dH2 concentration and increased abundance of the HydB subunit of electron-bifurcating group A3 [FeFe]-hydrogenases. The starch-rich diet caused higher rumen dH2 concentration and increased the abundance of genes encoding fermentative group B [FeFe]-hydrogenases, as well as energy-converting 4e and 4g [NiFe]-hydrogenases, which couple ferredoxin oxidation to H2 production in anaerobic bacteria, including in Firmicutes and Spirochaetota (Fig. 6A). Elevated rumen dH2 concentration with the starch-rich diet was accompanied with decreased abundance of genes encoding methanogenic [NiFe] and [Fe]-hydrogenases, as well as methyl-CoM reductases, which mediate hydrogenotrophic methanogenesis, which is in line with decreased CH4 production. A mechanism explaining less CH4 formation from concentrates based on faster rate of growth and/or lower maximal growth rate of methanogens resulting in greater H2 concentration, and consequently acetate and H2 formation becoming thermodynamically less favorable, with eventually less H2 being available for CH4 formation, has been proposed before [7]. Our finding of decreased abundance of methanogenic hydrogenases with the starch-rich diet is in line with that proposal, as fewer methanogens and methanogenic enzymes would presumably be present in the rumen with the faster rumen outflow rates and lower rumen pH associated to feeding concentrates. Greater H2 concentrations in turn thermodynamically favors H2-incorporating pathways such a propionate production, which agrees with the increased abundance of respiratory 1d [NiFe]-hydrogenases and putative sensory C3 [FeFe]-hydrogenases affiliated to Firmicutes, and supports hydrogenotrophic respiration (e.g., Selenomonas) [5].

Overall, our results indicate that the type of carbohydrate fed to ruminants alters rumen microbiome composition and function with distinct H2 metabolism, which is closely associated with rumen methanogenesis, and VFA production and profile (Fig. 6B). The fiber-rich diet enriched for fibrolytic bacteria and enhanced fiber utilization, and acetate and H2 production, along with increased methanogens abundance and greater CH4 production and decreased H2 concentrations. Enhanced acetate production by the fiber-rich selected microbiome might partly result from hydrogenotrophic acetogenesis, and, if confirmed, would be a novel metabolic strategy for ruminants to use H2 and CO2 with greater efficiency for the host compared to CH4 production. The starch-rich diet enriched for amylolytic bacteria, and enhanced propionate production through the acrylate pathway with lactate as an intermediate, helping to maintain a healthy rumen and decrease production of H2 available for methanogenesis. High starch intake increased fermentation with decreased methanogens abundance, and caused an increase in rumen dH2 concentration accompanied by greater abundance of hydrogenases involved in hydrogenotrophic respiration. These findings provide novel insights of distinct H2 metabolic pathways in rumen microbiota of cattle to adapt to interventions based on feeding carbohydrates with contrasting fermentability, and help to understand energy harvesting strategies, healthy rumen maintenance, and CH4 mitigation in ruminants. More work is needed to elucidate the adaptation, resilience, and colonization of bacteria with the metabolic flexibility to use H2 in the rumen, and their benefit for host animals.

All data are available in the main text or the in supplementary materials. Raw reads of 16S rRNA gene sequencing of ruminal microbiota are available at National Center for Biotechnology Information (NCBI, project number PRJNA736529). Raw reads of metagenomic sequencing of ruminal microbiota are available at NCBI (project number PRJNA779163).

Flint HJ. Polysaccharide breakdown by anaerobic microorganisms inhabiting the Mammalian gut. Adv Appl Microbiol. 2004;56:89–120.

Article CAS PubMed Google Scholar

Abbas W, Howard JT, Paz HA, Hales KE, Wells JE, Kuehn LA, et al. Influence of host genetics in shaping the rumen bacterial community in beef cattle. Sci Rep. 2020;10:15101.

Article CAS PubMed PubMed Central Google Scholar

Hungate RE. Hydrogen as an intermediate in the rumen fermentation. Arch Mikrobiol. 1967;59:158–64.

Article CAS PubMed Google Scholar

Robinson JA, Tiedje JM. Kinetics of hydrogen consumption by rumen fluid, anaerobic digestor sludge, and sediment. Appl Environ Microbiol. 1982;44:1374–84.

Article CAS PubMed PubMed Central Google Scholar

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–32.

Article CAS PubMed PubMed Central Google Scholar

Vignais PM, Billoud B, Meyer J. Classification and phylogeny of hydrogenases. FEMS Microbiol Rev. 2001;25:455–501.

Article CAS PubMed Google Scholar

Janssen PH. Influence of hydrogen on rumen methane formation and fermentation balances through microbial growth kinetics and fermentation thermodynamics. Anim Feed Sci Technol. 2010;160:1–22.

Article CAS Google Scholar

Wolf PG, Biswas A, Morales SE, Greening C, Gaskins HR. H2 metabolism is widespread and diverse among human colonic microbes. Gut Microbes. 2016;7:235–45.

Article CAS PubMed PubMed Central Google Scholar

Wang R, Bai Z, Chang J, Li Q, Hristov AN, Smith P, et al. China's low-emission pathways toward climate-neutral livestock production for animal-derived foods. Innovation. 2022;3:100220.

CAS PubMed PubMed Central Google Scholar

Wang F, Harindintwali JD, Yuan Z, Wang M, Wang F, Li S, et al. Technologies and perspectives for achieving carbon neutrality. Innovation. 2021;2:100180.

CAS PubMed PubMed Central Google Scholar

Henderson G, Cox F, Ganesh S, Jonker A, Young W, Abecia L, et al. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci Rep. 2015;5:14567.

Article CAS PubMed PubMed Central Google Scholar

Knapp JR, Laur GL, Vadas PA, Weiss WP, Tricarico JM. Invited review: enteric methane in dairy cattle production: quantifying the opportunities and impact of reducing emissions. J Dairy Sci. 2014;97:3231–61.

Article CAS PubMed Google Scholar

Beauchemin KA, McGinn SM. Methane emissions from feedlot cattle fed barley or corn diets. J Dairy Sci. 2005;83:653–61.

CAS Google Scholar

Godwin S, Kang A, Gulino LM, Manefield M, Gutierrez-Zamora ML, Kienzle M, et al. Investigation of the microbial metabolism of carbon dioxide and hydrogen in the kangaroo foregut by stable isotope probing. ISME J. 2014;8:1855–65.

Article CAS PubMed PubMed Central Google Scholar

Zhang Z, Xu D, Wang L, Hao J, Wang J, Zhou X, et al. Convergent evolution of rumen microbiomes in high-altitude mammals. Curr Biol. 2016;26:1873–9.

Article CAS PubMed Google Scholar

Joblin KN. Ruminal acetogens and their potential to lower ruminant methane emissions. Aust J Agric Res. 1999;50:1307–13.

Article Google Scholar

Owens FN, Secrist DS, Hill WJ, Gill DR. Acidosis in cattle: a review. J Anim Sci. 1998;76:275–86.

Article CAS PubMed Google Scholar

Malekkhahi M, Tahmasbi AM, Naserian AA, Danesh-Mesgaran M, Kleen JL, AlZahal O, et al. Effects of supplementation of active dried yeast and malate during sub-acute ruminal acidosis on rumen fermentation, microbial population, selected blood metabolites, and milk production in dairy cows. Anim Feed Sci Technol. 2016;213:29–43.

Article CAS Google Scholar

Kleen JL, Hooijer GA, Rehage J, Noordhuizen JPTM. Subacute Ruminal Acidosis (SARA) in dairy cows. J Vet Med. 2003;50:406–14.

Article CAS Google Scholar

Ungerfeld EM. Metabolic hydrogen flows in rumen fermentation: principles and possibilities of interventions. Front Microbiol. 2020;11:589.

Article PubMed PubMed Central Google Scholar

Wang M, Wang R, Tang SX, Tan ZL, Zhou CS, Han XF, et al. Comparisons of manual and automated incubation systems: effects of venting procedures on in vitro ruminal fermentation. Livest Sci. 2016;184:41–45.

Article Google Scholar

AOAC, Horwitz W. Official methods of analysis of AOAC International. 16th ed. Arlington, VA: Association of Official Analytical Chemists; 1995.

Google Scholar

Soest PJ, Robertson JB, Lewis BA. Symposium: carbohydrate methodology, metabolism, and nutritional implications in dairy cattle. J Dairy Sci. 1991;74:3585–97.

Google Scholar

Kartchner RJ, Theurer B. Comparison of hydrolysis methods used in feed, digesta, and fecal starch analysis. J Agric Food Chem. 1981;29:8–11.

Article CAS Google Scholar

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 Technol. 2014;194:1–11.

Article Google Scholar

Canale A, Valente ME, Ciotti A. Determination of volatile carboxylic acids (C1–C5i) and lactic acid in aqueous acid extracts of silage by high performance liquid chromatography. J Sci Food Agric. 1984;35:1178–82.

Article CAS Google Scholar

Ma ZY, Zhang XM, Wang R, Wang M, Liu T, Tan ZL. Effects of chemical and mechanical lysis on microbial DNA yield, integrity, and downstream amplicon sequencing of rumen bacteria and protozoa. Front Microbiol. 2020;11:581227.

Article PubMed PubMed Central Google Scholar

Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.

Article CAS PubMed Google Scholar

Ma Z, Wang R, Wang M, Zhang X, Mao H, Tan Z. Short communication: variability in fermentation end-products and methanogen communities in different rumen sites of dairy cows. J Dairy Sci. 2018;101:5153–8.

Article CAS PubMed Google Scholar

Seshadri R, Leahy SC, Attwood GT, Teh KH, Lambie SC, Cookson AL, et al. Cultivation and sequencing of rumen microbiome members from the Hungate1000 Collection. Nat Biotechnol. 2018;36:359–67.

Article CAS PubMed PubMed Central Google Scholar

Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP-a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome. 2018;6:158.

Article PubMed PubMed Central Google Scholar

Witten IH, Frank E, Hall MA, Pal CJ. Data mining: practical machine learning tools and techniques. 4th ed. Burlington: Morgan Kaufmann; 2016.

Kolde R, Laur S, Adler P, Vilo J. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. 2012;28:573–80.

Article CAS PubMed PubMed Central Google Scholar

Lane MA, Baldwin RLT, Jesse BW. Developmental changes in ketogenic enzyme gene expression during sheep rumen development. J Anim Sci. 2002;80:1538–44.

Article CAS PubMed Google Scholar

Helenius TO, Misiorek JO, Nystrom JH, Fortelius LE, Habtezion A, Liao J, et al. Keratin 8 absence down-regulates colonocyte HMGCS2 and modulates colonic ketogenesis and energy metabolism. Mol Biol Cell. 2015;26:2298–310.

Article CAS PubMed PubMed Central Google Scholar

Diernaes L, Sehested J, Moller PD, Skadhauge E. Sodium and chloride transport across the rumen epithelium of cattle in vitro: effect of short-chain fatty acids and amiloride. Exp Physiol. 1994;79:755–62.

Article CAS PubMed Google Scholar

Lin L, Xie F, Sun D, Liu J, Zhu W, Mao S. Ruminal microbiome-host crosstalk stimulates the development of the ruminal epithelium in a lamb model. Microbiome. 2019;7:83.

Article PubMed PubMed Central Google Scholar

Mao S, Zhang M, Liu J, Zhu W. Characterising the bacterial microbiota across the gastrointestinal tracts of dairy cattle: membership and potential function. Sci Rep. 2015;5:16116.

Article CAS PubMed PubMed Central Google Scholar

Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

Article CAS PubMed PubMed Central Google Scholar

Wang S, Huang H, Kahnt J, Thauer, Rudolf K. A reversible electron-bifurcating ferredoxin- and NAD-dependent [FeFe]-hydrogenase (HydABC) in Moorella thermoacetica. J Bacteriol. 2013;195:1267–75.

Article CAS PubMed PubMed Central Google Scholar

Adam PS, Borrel G, Gribaldo S. Evolutionary history of carbon monoxide dehydrogenase/acetyl-CoA synthase, one of the oldest enzymatic complexes. Proc Natl Acad Sci USA. 2018;115:E1166–E73.

Article CAS PubMed PubMed Central Google Scholar

Creevey CJ, Kelly WJ, Henderson G, Leahy SC. Determining the culturability of the rumen bacterial microbiome. Micro Biotechnol. 2014;7:467–79.

Article Google Scholar

Greening RC, Leedle JA. Enrichment and isolation of Acetitomaculum ruminis, gen. nov., sp. nov.: acetogenic bacteria from the bovine rumen. Arch Microbiol. 1989;151:399–406.

Article CAS PubMed Google Scholar

Kelly WJ, Henderson G, Pacheco DM, Li D, Reilly K, Naylor GE, et al. The complete genome sequence of Eubacterium limosum SA11, a metabolically versatile rumen acetogen. Stand Genom Sci. 2016;11:26.

Article Google Scholar

Genthner BRS, Davis CL, Bryant MP. Features of rumen and sewage-sludge strains of Eubacterium limosum, a methanol-utilizing and H2-CO2-utilizing species. Appl Environ Microbiol. 1981;42:12–19.

Article CAS PubMed PubMed Central Google Scholar

Ding XZ, Long RJ, Kreuzer M, Mi JD, Yang B. Methane emissions from yak (Bos grunniens) steers grazing or kept indoors and fed diets with varying forage:concentrate ratio during the cold season on the Qinghai-Tibetan Plateau. Anim Feed Sci Technol. 2010;162:91–98.

Article CAS Google Scholar

Huang X. Rumen microbial diversity of ruminants in Qinghai Tibetan Plateau, PhD dissertation. Lanzhou, China: Lanzhou University; 2013.

Google Scholar

Gharechahi J, Zahiri HS, Noghabi KA, Salekdeh GH. In-depth diversity analysis of the bacterial community resident in the camel rumen. Syst Appl Microbiol. 2015;38:67–76.

Article PubMed Google Scholar

Kayouli C, Jouany JP, Demeyer DI, Aliali, Taoueb H, Dardillat C. Comparative-studies on the degradation and mean retention time of solid and liquid-phases in the forestomachs of dromedaries and sheep fed on low-quality roughages from Tunisia. Anim Feed Sci Technol. 1993;40:343–55.

Article CAS Google Scholar

Varel VH, Fryda SJ, Robinson IM. Cellulolytic bacteria from pig large-intestine. Appl Environ Microbiol. 1984;47:219–21.

Article CAS PubMed PubMed Central Google Scholar

Goodson J, Tyznik WJ, Cline JH, Dehority BA. Effects of an abrupt diet change from hay to concentrate on microbial numbers and physical-environment in the cecum of the pony. Appl Environ Microbiol. 1988;54:1946–50.

Article CAS PubMed PubMed Central Google Scholar

Boulahrouf A, Fonty G, Gouet P. Establishment, counts, and identification of the fibrolytic microflora in the digestive-tract of rabbit - influence of feed cellulose content. Curr Microbiol. 1991;22:21–25.

Article Google Scholar

Kane MD, Breznak JA. Acetonema longum gen.nov.sp.nov., an H2/CO2 acetogenic bacterium from the termite, Pterotermes occidentis. Arch Microbiol. 1991;156:91–98.

Article CAS PubMed Google Scholar

Morvan B, Fonty G. Mixotrophy by rumen acetogenic bacteria in the utilization of hydrogen and sugars. Ann Zootech. 1996;45:354–354.

Article Google Scholar

Graber JR, Breznak JA. Physiology and nutrition of Treponema primitia, an H2/CO2-acetogenic spirochete from termite hindguts. Appl Environ Microbiol. 2004;70:1307–14.

Article CAS PubMed PubMed Central Google Scholar

Zhang X, Matson EG, Leadbetter JR. Genes for selenium dependent and independent formate dehydrogenase in the gut microbial communities of three lower, wood-feeding termites and a wood-feeding roach. Environ Microbiol. 2011;13:307–23.

Article CAS PubMed Google Scholar

Dougal K, de la Fuente G, Harris PA, Girdwood SE, Pinloche E, Geor RJ, et al. Characterisation of the faecal bacterial community in adult and elderly horses fed a high fibre, high oil or high starch diet using 454 pyrosequencing. PLoS ONE. 2014;9:e87424.

Article PubMed PubMed Central Google Scholar

Cordruwisch R, Seitz HJ, Conrad R. The capacity of hydrogenotrophic anaerobic-bacteria to compete for traces of hydrogen depends on the redox potential of the terminal electron-acceptor. Arch Microbiol. 1988;149:350–7.

Article CAS Google Scholar

Wang W M, Ungerfeld E, Degen AA, Jing X, Guo W, Zhou J, et al. Ratios of rumen inoculum from tibetan and small-tailed han sheep influenced in vitro fermentation and digestibility. Anim Feed Sci Technol. 2020;267:114562.

Article Google Scholar

Raju P. Homoacetogenesis as an alternative hydrogen sink in the rumen: a thesis presented in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Microbiology and Genetics at Massey University. Palmerston North, 2016.

Counotte GH, Prins RA, Janssen RH, Debie MJ. Role of Megasphaera elsdenii in the fermentation of dl-[2-C]lactate in the rumen of dairy cattle. Appl Environ Microbiol. 1981;42:649–55.

Article CAS PubMed PubMed Central Google Scholar

Mackie RI, Gilchrist FMC, Heath S. An in vivo study of ruminal micro-organisms influencing lactate turnover and its contribution to volatile fatty acid production. J Agric Sci. 1984;103:37–51.

Article CAS Google Scholar

Strmpl C, Jarvis GN. Anaerovibrio: Bergey's manual of systematics of archaea and bacteria. Hoboken: Wiley; 2015. p. 1–5.

Prive F, Kaderbhai NN, Girdwood S, Worgan HJ, Pinloche E, Scollan ND, et al. Identification and characterization of three novel lipases belonging to families II and V from Anaerovibrio lipolyticus 5ST. PLoS ONE. 2013;8:e69076.

Article CAS PubMed PubMed Central Google Scholar

Chen L, Shen Y, Wang C, Ding L, Zhao F, Wang M, et al. Megasphaera elsdenii lactate degradation pattern shifts in rumen acidosis models. Front Microbiol. 2019;10:162.

Article PubMed PubMed Central Google Scholar

Kamke J, Kittelmann S, Soni P, Li Y, Tavendale M, Ganesh S, et al. Rumen metagenome and metatranscriptome analyses of low methane yield sheep reveals a Sharpea-enriched microbiome characterised by lactic acid formation and utilisation. Microbiome. 2016;4:56.

Xie F, Jin W, Si H, Yuan Y, Tao Y, Liu J, et al. An integrated gene catalog and over 10,000 metagenome-assembled genomes from the gastrointestinal microbiome of ruminants. Microbiome. 2021;9:137.

Article CAS PubMed PubMed Central Google Scholar

Download references

We thank Zhi Peng Li, Sheng Yong Mao, and Qiang Qiu for discussions. This work was supported by National Natural Science Foundation of China (Grant Nos. 31922080, 32002204), Strategic Priority Research Program of the Chinese Academy of Sciences (Grant Nos. XDA26040203), Hunan Province Science and Technology Plan (2020NK2066, 2022NK2021), China Agriculture Research System of MOF and MARA, Youth Innovation Promotion Association CAS (Grant No. Y202078), Open Fund of Key Laboratory of Agro-ecological Processes in Subtropical Region Chinese Academy of Sciences (Grant No. ISA2021203). CG is supported by an NHMRC EL2 Fellowship (APP1178715).

These authors contributed equally: Qiu Shuang Li, Rong Wang.

Key Laboratory for Agro-Ecological Processes in Subtropical Region, Institute of Subtropical Agriculture, Chinese Academy of Sciences, Changsha, Hunan, China

Qiu Shuang Li, Rong Wang, Xiu Min Zhang, Jin Zhen Jiao, Zhi Liang Tan & Min Wang

University of Chinese Academy of Sciences, Beijing, China

Qiu Shuang Li, Zhi Liang Tan & Min Wang

College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, China

Zhi Yuan Ma

State Key Laboratory for Conservation and Utilization of Bio-Resources in Yunnan, School of Life Sciences, Yunnan University, Kunming, Yunnan, China

Zhi Gang Zhang

Centro Regional de Investigación Carillanca, Instituto de Investigaciones Agropecuarias (INIA), Temuco, Chile

Emilio M. Ungerfeld

Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China

Kang Le Yi & Bai Zhong Zhang

Xiangxi Cattle Engineering Technology Center of Hunan Province, Huayuan, Hunan, China

Kang Le Yi, Bai Zhong Zhang, Liang Long & Yun Long

Hunan De Nong Animal Husbandry Group Co. Ltd, Huayuan, Hunan, China

Liang Long & Yun Long

Shanghai BIOZERON Biotechnology Company Ltd, Shanghai, China

Ye Tao

Shanghai Institute of Nutrition and Health, Chinese Academy of Sciences, Shanghai, China

Tao Huang

Biomedicine Discovery Institute, Department of Microbiology, Monash University, Clayton, Australia

Chris Greening

You can also search for this author in PubMed Google Scholar

Conceptualization and research design: MW, QSL. Research conduction and data acquisition: RW, QSL, ZGZ, YL. Data analysis: QSL, RW, YT, TH, MW, JZJ, ZYM, CG. Investigation: KLY, BZZ, ZLT, LL, MW. Writing—original draft: QSL, MW, ZGZ, ZLT, CG, EMU. Writing—reviewing & editing: all authors.

Correspondence to Zhi Liang Tan or Min Wang.

The authors declare no competing interests.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Li, Q.S., Wang, R., Ma, Z.Y. et al. Dietary selection of metabolically distinct microorganisms drives hydrogen metabolism in ruminants. ISME J 16, 2535–2546 (2022). https://doi.org/10.1038/s41396-022-01294-9

Download citation

Received: 02 February 2022

Revised: 29 June 2022

Accepted: 11 July 2022

Published: 05 August 2022

Issue Date: November 2022

DOI: https://doi.org/10.1038/s41396-022-01294-9

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative