Skip to main content

H2 generated by fermentation in the human gut microbiome influences metabolism and competitive fitness of gut butyrate producers

Abstract

Background

Hydrogen gas (H2) is a common product of carbohydrate fermentation in the human gut microbiome and its accumulation can modulate fermentation. Concentrations of colonic H2 vary between individuals, raising the possibility that H2 concentration may be an important factor differentiating individual microbiomes and their metabolites. Butyrate-producing bacteria (butyrogens) in the human gut usually produce some combination of butyrate, lactate, formate, acetate, and H2 in branched fermentation pathways to manage reducing power generated during the oxidation of glucose to acetate and carbon dioxide. We predicted that a high concentration of intestinal H2 would favor the production of butyrate, lactate, and formate by the butyrogens at the expense of acetate, H2, and CO2. Regulation of butyrate production in the human gut is of particular interest due to its role as a mediator of colonic health through anti-inflammatory and anti-carcinogenic properties.

Results

For butyrogens that contained a hydrogenase, growth under a high H2 atmosphere or in the presence of the hydrogenase inhibitor CO stimulated production of organic fermentation products that accommodate reducing power generated during glycolysis, specifically butyrate, lactate, and formate. Also as expected, production of fermentation products in cultures of Faecalibacterium prausnitzii strain A2-165, which does not contain a hydrogenase, was unaffected by H2 or CO. In a synthetic gut microbial community, addition of the H2-consuming human gut methanogen Methanobrevibacter smithii decreased butyrate production alongside H2 concentration. Consistent with this observation, M. smithii metabolic activity in a large human cohort was associated with decreased fecal butyrate, but only during consumption of a resistant starch dietary supplement, suggesting the effect may be most prominent when H2 production in the gut is especially high. Addition of M. smithii to the synthetic communities also facilitated the growth of E. rectale, resulting in decreased relative competitive fitness of F. prausnitzii.

Conclusions

H2 is a regulator of fermentation in the human gut microbiome. In particular, high H2 concentration stimulates production of the anti-inflammatory metabolite butyrate. By consuming H2, gut methanogenesis can decrease butyrate production. These shifts in butyrate production may also impact the competitive fitness of butyrate producers in the gut microbiome.

Video Abstract

Background

Hydrogen gas (H2) is a common product of bacterial metabolism in anoxic environments, when electron acceptors for anaerobic respiration are limited. H2 is commonly produced when fermentative microbes use protons as electron acceptors to dispose of reducing power, reducing them to H2 via hydrogenases [1,2,3]. Thermodynamic principles render H2 production less favorable when H2 concentrations are high, impacting the metabolism of H2-producing microbes [4,5,6]. Hydrogenase genes occur in phylogenetically diverse microbes including 71% of the reference genomes in the Human Microbiome Project Gastrointestinal Tract database, suggesting that H2 concentration may be a major factor influencing fermentation in the human gut microbiota [7].

H2 produced during bacterial fermentation in the large intestine can be consumed by other microbes, escape in flatus, or diffuse into the blood stream where it is subsequently released into the lungs and exhaled. The summation of these processes results in a concentration of H2 in intestinal gas ranging from undetectable to over 40% v/v (Supplementary Fig. 1) [8, 9]. Diet is a major determinant of H2 production in the human colon. In particular, fermentable, microbiota-accessible carbohydrates (MACs) [10] largely drive H2 production [11, 12]. Despite the ubiquity of H2 in the environment of the large intestine, concrete information is lacking about how the concentration of hydrogen regulates fermentation of specific gut microbes. An effect of H2 concentration on human gut butyrogens has been predicted [13] and would be of particular significance because of the anti-inflammatory and anticarcinogenic effects of butyrate [14,15,16,17].

The fermentation scheme of typical human gut butyrogens is depicted in Fig. 1A. Carbohydrate substrates (most simply represented by glucose) are first processed via glycolysis. Per glucose, the reactions of glycolysis form two pyruvate and phosphorylate two ADP to form ATP by substrate-level phosphorylation (SLP). Importantly, glucose oxidation to pyruvate reduces two moles of the cofactor NAD+ to NADH [18]. The necessity of regenerating NAD+ from this NADH to maintain redox balance represents both a central constraint on possible fermentations and an opportunity to conserve additional energy [2, 19].

Fig. 1
figure 1

Stoichiometry and thermodynamics of fermentation in human gut butyrogens. A Generic fermentation pathways in human gut butyrogens that can yield combinations of H2, CO2, formate, lactate, acetate, and butyrate (adapted from Louis and Flint, 2017 [18]). Fd ferredoxin, LDH lactate dehydrogenase, PFOR pyruvate:ferredoxin oxidoreductase, PFL pyruvate formate-lyase, Ack acetate kinase, Hyd ferredoxin hydrogenase, Bhbd hydroxybutyryl-CoA dehydrogenase, Etf-Bcd butyryl-CoA dehydrogenase electron-transferring flavoprotein complex, But butyryl-CoA:acetyl CoA transferase, Rnf complex. The division of acetyl-CoA between acetate production (green box) and butyrate production (blue box) is variable; an equal division is shown as a representative case with simple stoichiometry. The division of pyruvate between the PFOR and PFL routes of acetyl-CoA formation is also variable. Stoichiometry is balanced from glucose to either lactate (orange box) or to butyrate and acetate, via either PFOR or PFL. B Illustration of the effect of [H2] on ∆G of glucose fermentation to butyrate and acetate. A range of possible fermentation balances are shown, with the net molar production or consumption of butyrate (but), acetate (ace), and H2 per glucose indicated. Dotted horizontal lines indicate theoretical ∆G thresholds for formation of a given number of ATP, assuming ∆G of ADP phosphorylation as + 70 kJ/mol [19]. Dashed line indicates the least butyrate formation likely to be feasible for human gut butyrogens. Representative physiological conditions were used as indicated. See Supplementary Table 4 for complete reactions and thermochemical parameters

Following glycolysis, pyruvate may be reduced to lactate by lactate dehydrogenase (LDH). This pathway fully reoxidizes the NADH produced in glycolysis, but produces no additional ATP [20]. More commonly, pyruvate is converted into acetyl-CoA and either CO2 or formate. Production of acetyl-CoA and CO2 is catalyzed by pyruvate:ferredoxin oxidoreductase (PFOR) and is coupled with the reduction of ferredoxin (Fd), a small iron–sulfur protein with a lower reduction potential than NAD+ [19, 21, 22]. By contrast, pyruvate cleavage into acetyl-CoA and formate is catalyzed by pyruvate formate-lyase (PFL) and does not generate any additional reduced species [21, 23]. Formate may subsequently be used in anabolic pathways or simply secreted as a fermentation product [23]. In vivo, both the PFOR and PFL pathways of acetyl-CoA formation can be active simultaneously [24].

As with pyruvate, acetyl-CoA in human gut butyrogens can proceed down either of two branched pathways culminating in acetate or butyrate production [2, 13, 18]. In acetate production, the acetyl group is transferred from CoA to phosphate to form acetyl phosphate (acetyl-P). In a reaction catalyzed by acetate kinase (Ack), this phosphate is then transferred to ADP to generate ATP via SLP, releasing acetate. This pathway conserves energy as ATP but does not contribute to reoxidation of NADH or reduced ferredoxin (Fdred) [2, 13, 25]. Butyrate production, by contrast, is an important sink for reducing power. In this pathway, two acetyl-CoA are combined to form acetoacetyl-CoA, which is then reduced to 3-hydroxybutyryl-CoA by hydroxybutyryl-CoA dehydrogenase (Bhbd) in a reaction that reoxidizes one NADH cofactor to NAD+ [26]. The next reaction forms crontonyl-CoA, which is then further reduced to butyryl-CoA by the butyryl-CoA dehydrogenase electron-transferring flavoprotein complex (Bcd-Etf). This electron-bifurcating complex couples crotonyl-CoA reduction to the endergonic reduction of Fdox by NADH in an overall thermodynamically feasible reaction [19, 27]. The resulting Fdred, in addition to that formed by PFOR, is reoxidized by ferredoxin hydrogenase through the reduction of protons to H2 [18, 28]. Sufficient flux through the butyrate production pathway results in an overabundance of Fdred and a shortfall of NADH. Butyrogens can conserve additional energy in this case through anaerobic respiration using the Rnf complex, which couples Fdred oxidation and NAD+ reduction (ferredoxin:NAD+ oxidoreductase) to cation transport across the membrane, synthesizing ATP by a chemiosmotic mechanism [19, 29]. Final release of butyrate is by exchange with free acetate catalyzed by butyryl-CoA:acetyl CoA transferase (But), which results in net acetate consumption when more acetyl-CoA is flows to butyrate production than to acetate production [18, 25, 30].

A salient feature of the branched metabolism of gut butyrogens described above is that increasing butyrate production reduces the reductant available for H2 formation. Conversely, increasing acetate production necessarily entails the formation of more H2 in order to regenerate NAD+ [2, 18]. As a result, thermodynamic equilibrium increasingly favors butyrate over acetate production as ambient H2 concentrations rise (illustrated in Fig. 1B) [6, 13, 31]. This has led to the expectation that higher gut H2 concentration favors butyrate over acetate production [6, 13]. While this result or similar results demonstrating the same principle have been observed in species related to human gut butyrogens [32,33,34], to the best of our knowledge, the effect of ambient H2 concentration on fermentation by the predominant human gut butyrogens Roseburia intestinalis, Eubacterium rectale, and Faecalibacterium prausnitzii [13, 35] has not been directly investigated, with the exception of one study that found no effect of autogenous H2 on cultures of R. intestinalis [36].

H2 concentration in bacterial culture can be reduced by stirring [34, 37], sparging [32], or co-culture with hydrogen-consuming microbes (hydrogenotrophs) [4, 33, 38, 39]. Under these circumstances, hydrogenases can generate more H2 and oxidized end products (e.g., acetate), with correspondingly less production of reduced organic end products (e.g., ethanol, butyrate) [32,33,34, 37, 39]. Three guilds of hydrogenotrophs colonize the human gut: methanogens, sulfate reducers, and reductive acetogens [40,41,42,43,44]. These hydrogenotrophs actively consume H2 and may therefore play an important role in regulating the H2 concentrations to which human gut butyrogens are exposed.

In this study, we investigate the effect of H2 concentration on the profile of fermentation end products of R. intestinalis, E. rectale, and F. prausnitzii. We find evidence that physiologically relevant variations in H2 concentration influence the favored routes of reductant disposal in H2-producing human gut butyrogens, resulting in shifts in the production of fermentation products. Specifically, exposure to high H2 concentrations increases production of butyrate, lactate, and formate at the expense of acetate and, presumably, CO2. These metabolic shifts appear to impact the competitive fitness of certain butyrogens. We propose a model where the profile of fermentation products from these taxa, and metabolically similar fermenters in the human colon, is modulated by local colonic H2 concentration. This, in turn, is a balance of production by fermenters and elimination by hydrogenotrophs. Finally, we report observations from a large human cohort consuming resistant potato starch as a MAC expected to stimulate fermentation and H2 production. Consistent with our model, we found that hydrogenotrophic gut methanogenesis was associated with decreased fecal butyrate during supplement consumption.

Methods

Human cohort

Results from a portion of this study’s human cohort were previously reported by Baxter et al. (2019) [45]. Participants were recruited through Authentic Research Sections of the University of Michigan BIO173 introductory biology course. Subjects were excluded based on self-reported inflammatory bowel syndrome, inflammatory bowel disease, colorectal cancer, and consumption of antibiotics in the past 6 months. De-identified human subject metadata including age and sex is provided in Supplementary Table 1.

Microbial strains and culture

Methanobrevibacter smithii F1 (DSM 2374), Faecalibacterium prausnitizii A2-165 (DSM 17,677), and Roseburia intestinalis L1-82 (DSM 14,610) were obtained from the German Collection of Microorganisms and Cell Cultures GmbH (DSMZ). Ruminococcus bromii VPI 6883 (ATCC 27,255) was obtained from the American Type Culture Collection (ATCC). Eubacterium rectale A1-86 (DSM 17,629), Bacteroides thetaiotaomicron VPI 5482 (DSM 2079), Bacteroides vulgatus Eggerth and Gagnon (ATCC 8482), and Prevotella copri CB7 (DSM 18,205) were obtained from collaborators.

Also included in the synthetic community were the strains Bifidobacterium adolescentis 269–1 and Anaerostipes caccae 127–8-5, which are isolates from fecal samples obtained in the course of the human cohort study. B. adolescentis 269–1 was obtained from a fecal sample serially diluted and plated on Bifidus Selective Medium agar (BSM Agar, Sigma-Aldrich) including the BSM supplement according to the manufacturer’s instructions. Plates were incubated at 37 °C in an anoxic atmosphere of 5% carbon dioxide, between 1.5 and 3.5% H2, and balance N2 in an anaerobic chamber (Coy Laboratory Products Inc., Grass Lake, MI). Bifidobacterium colonies were identified as having a pink center and light brown edge and were restreaked on BSM agar.

A. caccae 127–8-5 was obtained from a fecal sample stored at − 80 °C in an OMNIgene-Gut collection kit tube (DNA Genotek, Ottawa, Ontario, Cat#OMR-200). The fecal sample was serially diluted and plated on SABU agar, a medium containing 2 g/L taurocholate to stimulate spore germination (full list of medium components in Supplementary Table 2). Plates were incubated at 37 °C in the anaerobic chamber described above, and colonies that grew were picked.

The taxonomic identity of the 269–1 and 127–8-5 isolates was determined using Sanger sequencing of the 16S rRNA gene. 16S rRNA was amplified using primers designated 8F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 1492R (5′- GGTTACCTTGTTACGACTT-3′) and sequenced from the 8F primer. These sequencing results have been deposited in Zenodo (https://doi.org/10.5281/zenodo.6643453).

All microbial strains used in this study are available upon request made to the lead contact.

All microbial strains were preserved in frozen stocks at − 80 °C with either 5% DMSO or 20% glycerol as a cryopreservative. To begin cultivation, for all strains except M. smithii F1, a small amount of material was scraped from the frozen stocks and added to 5–10 mL of SAB4 base medium (components in Supplementary Table 2) supplemented with either 4 g/L D-glucose or 2 g/L each of D-glucose and D-fructose in the Coy anaerobic chamber described above. These cultures were incubated at 37 °C and passaged as necessary (no more than four passages, most commonly one or two) to produce mid- or late-exponential phase cultures used to inoculate experimental cultures. For M. smithii F1, frozen stocks were thawed in the anaerobic chamber and transferred using a 1-mL syringe fitted with a 23-gauge needle into a Balch tube sealed with a butyl rubber stopper and aluminum crimp (Chemglass Life Sciences, Vineland, NJ, Cat#CLS-4209) containing 5 mL SAB4 base medium under a headspace of 80% H2 + 20% CO2 mixed gas at 20 psig. These primary cultures were incubated at 37 °C on an orbital shaker at 150 rpm and passaged anaerobically as necessary to produce mid- or late-exponential phase cultures used to inoculate experimental cultures.

Monoculture experiments

Cultures were grown at 37 °C in 10 mL (H2 headspace experiments) or 5 mL (CO headspace experiments) SAB4 base medium (Supplementary Table 2) supplemented with 4 g/L D-glucose and 2.31 g sodium bicarbonate in Balch tubes sealed with butyl rubber stoppers and aluminum crimps (as described in Experimental Model and Subject Details). For shaking cultures, the Balch tubes were placed on their side in an orbital shaker at 150 rpm.

All Balch tubes were prepared with a headspace of 80% N2 + 20% CO2 mixed gas at atmospheric pressure. In experiments involving the addition of H2 to the headspace, all cultures were prepared with a headspace at 3 atm gauge pressure containing the indicated partial pressure of H2 and the balance N2. Gases were added using a custom gas manifold, and a pressure gauge was used to adjust regulators to supply the correct pressure (SSI Technologies Inc., Janesville, WI, Cat#MG-30-A-9 V-R). Ultra-high purity grades of N2 and H2 were used. In experiments involving the addition of carbon monoxide (CO), 2.2 mL of either pure CO or N2 was added using a syringe fitted with a stopcock and needle. All gases used in this study were purchased from Cryogenic Gases Inc., a division of Metro Welding Supply Corp. (Detroit, MI).

Growth curves were obtained by making regular measurements of the OD600 in the culture tubes using a Spec-20 spectrophotometer (Thermo Spectronic Model 333,183). Before each series of measurements, the spectrophotometer was zeroed using a Balch tube containing uninoculated medium from the same batch used in the experiment.

Monoculture experiments under H2 were performed twice for E. rectale and F. prausnitzii and three times for R. intestinalis. One experiment each for E. rectale and R. intestinalis included conditions with ppH2 of 2 and 3 atm in addition to 0 and 1 atm. All conditions in all experiments under H2 had three to five biological replicates. Monoculture experiments with CO were performed once for each butyrogen with four (R. intestinalis, F. prausnitzii) or five (E. rectale) biological replicates.

Synthetic community experiments

Cultures of synthetic community members (excluding M. smithii F1) were grown from stock in SAB4 base medium supplemented with 2 g/L each of D-glucose and D-fructose in the anaerobic chamber described above (and passaged so as to obtain mid- or late-exponential phase cultures of all the microbes simultaneously (as described above in Experimental Model and Subject Details)). Once this was achieved, equal cell numbers of each synthetic community member (estimated using OD600 measurements) were combined to create an inoculation mix, which was used to inoculate Balch tubes for the experimental cultures, which were subsequently sealed. These Balch tubes contained 10 mL of the SAB4 base medium supplemented with D-glucose and D-fructose described above. Since they were inoculated and sealed in the anaerobic chamber, their initial headspace matched that of the anerobic chamber (5% carbon dioxide, between 1.5 and 3.5% H2, balance N2). Each experimental culture was grown for 24 h, then passaged at a 1:100 dilution into another Balch tube for two subsequent 24-h cultures. For shaking cultures, the Balch tubes were placed on their side in an orbital shaker at 150 rpm.

Cultures of M. smithii F1 were grown from stock in Balch tubes and passaged so as to obtain mid- or late-exponential phase cultures at the same time as the other synthetic community members. M. smithii cells were added to the appropriate experimental cultures as an inoculum separate from the inoculation mix described above. Additional inocula of M. smithii were added at each passage of the synthetic community from pure cultures of M. smithii that were maintained in Balch tubes during the course of the experiment. The number of M. smithii cells added in each inoculum was estimated using OD600 and kept consistent.

The synthetic community experiment was performed twice with five biological replicates for each condition each time. Shaking cultures were only included in the second experiment.

Synthetic community relative abundance quantification

One milliliter samples of synthetic community cultures were centrifuged for 2 min at 11,000 g. Genomic DNA was extracted from the pellet using a DNeasy PowerLyzer Microbial Kit (Qiagen, Cat#12,255–50) according to the manufacturer’s instructions. The V4 region of the 16S rRNA gene was amplified and sequenced on the Illumina MiSeq platform using a 2 × 250-bp paired-end kit as described in Kozich et al. (2013) [46].

The resulting 16S amplicon data was analyzed using mothur v1.39.5 [47]. The mothur script and logfile have been deposited in Zenodo (https://doi.org/10.5281/zenodo.6621661). In summary, paired-end reads were merged into contigs, screened for sequencing errors, and aligned to the SILVA v132 reference database [46]. Aligned sequences were pre-clustered at 1 difference, screened for chimeras, and classified using the SILVA v132 reference database. Sequences identified as mitochondria, chloroplasts, or eukaryotes were removed. Sequences were then clustered into 99% OTUs, which reproduced the 9 community members (plus M. smithii) known to be present in the cultures, and a shared file was exported. Microsoft Excel (Microsoft Corporation, Redmond, WA) was used to calculate relative abundances from the shared file, and the results were imported into GraphPad Prism 9 (GraphPad Software, San Diego, CA), where statistical analyses were carried out as described in the figure legends.

Aqueous fermentation product quantification

Samples of 1 mL bacterial culture were centrifuged for 2 min at 11,000 g and the supernatant passed through a 0.22-µm MultiScreenHTS GV 0.22-µm filter plate (Millipore Sigma, Burlington, MA). Similar to the procedure described by Baxter et al. (2019) [45], filtrates were transferred into 100-µl inserts inside 1.5-ml screw cap vials in preparation for analysis by HPLC. Quantification of SCFAs was performed using a Shimadzu HPLC system (Shimadzu Scientific Instruments, Columbia, MD) that included an LC-10AD vp pump A, LC-10AD vp pump B, DGU-14A degasser, CBM-20A communications bus module, SIL-20AC HT autosampler, CTO-10AS vp column oven, RID-10A RID detector, and an Aminex HPX-87H column (Bio-Rad Laboratories, Hercules, CA). We used a mobile phase of 0.01 N H2SO4 at a total flow rate of 0.6 ml per min with the column oven temperature at 50 °C. The sample injection volume was 10 µl, and each sample eluted for 40 min. The concentrations were calculated using standard curves generated for each product from a cocktail of short-chain organic acid standards at concentrations of 40, 20, 10, 5, 2.5, 1, 0.5, 0.25, and 0.1 mM. These standards were run before and after each batch of samples, and standard curves were generated using averaged values. The baseline of the chromatographs was manually corrected to ensure consistency between samples and standards. Samples were analyzed in a randomized order.

Gaseous fermentation product quantification

Gas samples were removed from the headspace of cultures using syringes fitted with stopcocks. Methane content was measured using a Shimadzu GC-2014A greenhouse gas analyzer gas chromatograph (Shimadzu Scientific Instruments, Inc., Columbia, MD) equipped with a flame ionization detector (FID) fed by ultra-high purity H2 and zero-grade air. Ultra-high purity N2 was used as the carrier gas. Sample separation was performed with a 1.0-M Hayesep T 80/100 mesh column, a 4.0-M Hayesep D 80/100 mesh column, and a 0.7-M Shimalite Q 100/180 mesh column. Before each series of measurements, accuracy was checked using a 500-ppm methane standard (Argus-Hazco, Byron Center, MI, Cat#GD40-007-A-221S).

H2 content was measured using a Peak Performer 1 gas chromatograph (Cat#910–105) with a reducing compound photometer (RCP) detector and post-column diluter (Peak Laboratories, Mountain View, CA) calibrated using a 10-ppm H2 standard (GASCO 105L-H2N-10, Cal Gas Direct Incorporate, Huntington Beach, CA). Ultra-high purity N2 was used as the carrier gas. When necessary, samples were diluted in room air using syringes fitted with stopcocks before measurement to reduce the H2 concentration below 100 ppm, which was the upper detection limit.

Total protein quantification

Total protein content was used as an indicator of bacterial biomass; 1 mL samples of microbial cultures at endpoint were centrifuged for 2 min at 11,000 g. The supernatant fraction was stored at − 80 °C for later analysis. The pellet was resuspended in 1.5 mL distilled H2O and sonicated to lyse cells. Sonication was performed on ice using a Branson Digital Sonifier 450 equipped with a 102C converter and microtip, which was placed directly in the bacterial suspension; 35% amplitude was used for a 3-min cycle of 1 s on followed by 14 s off (total 12 s sonication time). Protein concentrations in the resulting lysate and the saved supernatant fraction were determined using a Pierce Coomassie Plus Bradford assay reagent (Thermo Scientific, Cat#23,238) with bovine serum albumin (BSA) standards per the manufacturer’s instructions. The results from the lysate and supernatant were added together to obtain the total protein yield of the culture.

Human cohort study design and sample collection

The study took place during a number of separate semesters over the course of 3 years, from the winter semester of 2016 to the winter semester of 2019. While all supplements consumed consisted of resistant starch from potato (RSP), they varied in source, total dose, and frequency. The supplements consumed were Bob’s Red Mill potato starch (Bob’s Red Mill Natural Foods, Milwaukie, OR) consumed as a 20-g dose twice daily, a 20-g dose mixed with 2.5 g psyllium twice daily, a 40-g dose once daily, or a 40-g dose twice daily; or resistant potato starch from LODAAT Pharmaceuticals (Oak Brook, IL) consumed as a 20-g dose once daily. The supplement and dosage consumed by each subject is documented in Supplementary Table 1.

In each semester, the study followed a 3-week course. During the first week, fecal and breath samples were collected before consumption of RSP. During the second week, RSP consumption began at a half dose and increased to the full dose. During the third week, RSP consumption continued at the full dose while fecal and breath samples were collected.

Human sample analysis

Fecal sample collection, preparation, and quantification of short-chain fatty acid concentration by high-performance liquid chromatography (HPLC) was performed as previously described in Baxter et al. (2019) [45]. Breath samples consisted of 30 mL of end-expiratory breath collected in a 30-mL gastight syringe. Immediately after collection, samples were injected into a QuinTron BreathTracker SC analyzer (QuinTron Instrument Company Inc., Milwaukee, WI, Cat#QTLNRBTGCSC) for analysis. Concentrations of H2, methane, and carbon dioxide gas were measured, and hydrogen and methane measurements were normalized based on an assumed nominal concentration of 3.5% carbon dioxide. The BreathTracker analyzer was calibrated daily using a standard calibration gas containing 150 ppm H2, 75 ppm methane, and 6% carbon dioxide (QuinTron, Cat#QT07500-G). The quantifications of fecal butyrate in each fecal sample and H2 and methane in each breath sample are provided in Supplementary Table 1.

Microbial culture fermentation products

In monocultures of R. intestinalis, E. rectale, and F. prausnitzii, depending on the moles of butyrate formed per glucose fermented, acetate can be either a net product (< 1 mol butyrate per mol glucose) or net substrate (> 1 mol butyrate per mol glucose) of fermentation. Even when it is a net substrate, however, some acetyl-CoA still flows to acetate production and ATP formation by acetate kinase. In order to not obscure this nuance by reporting the net consumption of acetate that was observed in some cultures, results for fermentation products were expressed as the percent of total carbon consumed that was used in the formation of each product, rather than simple carbon recovery. This metric was calculated by first reasoning that since all three of these butyrogens produce butyrate by consuming acetate via the butyryl-CoA:acetate CoA enzyme, each mole of butyrate produced represented a mole of acetate (abundantly available in the SAB4 medium) consumed [30]. Therefore, a molar value of acetate that was theoretically consumed was set equivalent to moles of produced butyrate. Total fermented carbon was then calculated by adding the moles of carbon in the consumed glucose to the moles of carbon in theoretically consumed acetate. Total acetate produced was then calculated by adding the moles of theoretically consumed acetate to the change in acetate measured in the culture at endpoint versus blank medium, which varied from consumption to production depending on strain and condition. Since neither butyrate, formate, nor lactate were present in the medium or expected to be consumed during microbial metabolism, their total quantity produced was simply taken to be their endpoint concentration. Percent total fermented carbon in each substrate was then calculated as the moles of carbon in the produced substrate divided by the total moles of fermented carbon.

In the synthetic communities, the presence of diverse potential metabolic pathways rendered the above approach impractical. Instead, results for each fermentation product were expressed as change in the product in moles (an increase in all but one culture where acetate decreased) divided by the total moles of substrate (glucose and fructose) consumed in the culture.

All statistical analyses of metabolite data obtained from microbial cultures were performed using GraphPad Prism 9 and are described in the figure legends.

Human samples and methanogenesis classification

Concentrations of fecal metabolites obtained from HPLC (described above) were normalized to the wet weight of fecal material. Concentrations of CH4 and H2 in breath samples were quantified as described above. For each fecal metabolite and breath gas, samples with values lying more than three interquartile ranges below the lower quartile or above the upper quartile were excluded from analysis according to the method of Tukey’s Fences [48].

Subjects were classified as methanogenic or non-methanogenic, with separate classifications for the periods before and during supplement consumption. Methanogenic subjects were defined as those with over 4 ppm methane in at least one breath sample. This cutoff was based on a study of responses to consumption of lactulose (a fiber inaccessible to human enzymes but rapidly degraded by the gut microbiota) which suggested that a baseline threshold of 4 ppm above background is best predictive of increased breath methane [49]. We used this threshold because our intent was to identify subjects where methanogenesis was not just present, but a significant component of the gut ecosystem. In the Winter 2016 and Winter 2019 semesters, an elevated baseline concentration of 1 ppm methane was observed across most samples. Since this was likely due to instrument calibration rather than biological activity, this elevated baseline was subtracted before classifying individuals as methanogenic or non-methanogenic. The average concentration of each fecal metabolite and breath gas before and during starch supplement consumption was then calculated for each subject. The average concentration of fecal metabolites and breath gases in methanogenic and non-methanogenic subjects was then compared using two-tailed Student’s t-tests in GraphPad Prism 9.

Results

To test the hypothesis that H2 modulates the production of fermentation products by human gut butyrogens, we studied pure cultures of strains representing abundant butyrogens in the human gastrointestinal tract. Eubacterium rectale A1-86 and Roseburia intestinalis L1-82 were selected to represent the generalized metabolic pathways of butyrogens that could be affected by H2 (Fig. 1). Faecalibacterium prausnitzii A2-165 was chosen as a representative of butyrogens that lack a hydrogenase and are therefore unlikely to be affected by H2.

Replicate cultures of these three butyrogens were grown under a headspace of either H2 or N2 and shaken continuously to equilibrate headspace gases with the culture medium. As predicted, the presence of a H2 headspace shifted the profile of fermentation away from acetate towards more reduced organic acids (e.g., lactate and butyrate) for both H2-producing butyrogens (Fig. 2A, C). The same pattern of fermentation products was recapitulated in the presence of carbon monoxide (CO; Fig. 2B, D), a potent inhibitor of ferredoxin hydrogenase [50]. The profiles of reduced organic acids differed between the hydrogen-producing butyrogens. In cultures of R. intestinalis, reducing power was diverted to butyrate and formate. Lactate production increased very significantly, but remained only a trace product in both conditions (< 0.5% substrate carbon). By contrast, cultures of E. rectale, saw reducing power diverted primarily to lactate, with a smaller diversion to formate and no change in butyrate. Since formate production reduces intracellular Fdred versus the alternative production of CO2, it clearly represents a diversion reducing power away from H2 production via ferredoxin hydrogenase, as does butyrate production (Fig. 1A). Increasing the partial pressure of H2 in the headspace up to 3 atm led to larger shifts in a rough dose–response pattern (Supplementary Fig. 2). Unlike E. rectale and R. intestinalis, F. prausnitzii lacks hydrogenase activity [51]. As expected, its fermentation products were unaffected by the presence of H2 or CO in the headspace (Fig. 2E, F).

Fig. 2
figure 2

Variation in fermentation products in cultures of human gut butyrogens grown under different atmospheres. Endpoint fermentation products in cultures of R. intestinalis (note lactate is shown on smaller scale on right axis) (A–B), E. rectale (C–D), and F. prausnitzii (E–F) grown in shaken cultures with H2, N2, or CO—a potent inhibitor of hydrogenases. Error bars indicate SEM. Statistical significance calculated using two-sided Student’s two-sample t-tests (*p < 0.05; **p < 0.01; ***p < 0.001)

To assess whether H2 has an impact on butyrate production in more complex microbial communities, we assembled a synthetic community of microbes isolated from the human gut. This mixture consisted of representative strains of common butyrate producers (F. prausnitzii, E. rectale, R. intestinalis, and Anaerostipes caccae), two common fiber degraders (Bifidobacterium adolescentis and Ruminococcus bromii), and several members of the abundant gut phylum Bacteroidetes (Bacteroides vulgatus, Bacteroides thetaiotaomicron, and Prevotella copri). We compared butyrate production by this community from equimolar quantities of glucose and fructose in the presence or absence of M. smithii, the predominant H2-consuming methanogen in the human gut [52, 53].

The relative abundances of the constituent microbes, as well as production of CH4, H2, and fermentation products, were monitored over the course of three consecutive subcultures of the synthetic community. The addition of M. smithii resulted in the production of methane (Fig. 3A) and decreased the concentration of H2 and production of butyrate as predicted (Fig. 3B, C). The corresponding increase in acetate, which was observed in monocultures, was only observed in the second subculture of the synthetic community (Supplementary Fig. 3D). The expected shift in acetate production by butyrogens may have been masked by the copious production of acetate by other members of the synthetic community, such as the two Bacteroides species. Lactate and formate were also produced, along with the characteristic Bacteroides fermentation products propionate and succinate (Supplementary Fig. 3E-H). Lactate production was reduced by the addition of M. smithii. Since M. smithii consumes all available formate, its impact on formate production could not be determined.

Fig. 3
figure 3

Influence of methanogenesis on butyrate production products in a synthetic gut community and the human gut microbiota. Methane (A), H2 (B), and butyrate (C) production by a 9-species synthetic human gut community grown with (shaded bars) and without (open bars) the addition of M. smithii. Butyrate was measured in three successive subcultures, with a new inoculum of M. smithii added at each passage. In a human cohort consuming a resistant potato starch supplement, breath and fecal samples were used to determine weekly average breath CH4 (D), breath H2 (E), and fecal butyrate (F) in individuals with and without active gut methanogenesis, defined as at least one breath methane measurement greater than 4 ppm CH4 in the measurement period (MG) or no breath measurements over 4 ppm in the same period (non-MG). Error bars indicate SEM. Statistical significance calculated using two-sided Student’s two-sample t-tests (*p < 0.05; **p < 0.01; ***p < 0.001). One breath H2 measurement shown in the methanogenic group with a value over 60 ppm (data point shown in red) was excluded from statistical analysis

Removal of M. smithii from the synthetic community increased H2 levels and favored the growth of F. prausnitzii (Fig. 4). The increase of F. prausnitzii suggests that accumulation of H2 in the absence of a methanogen forced butyrogens with ferredoxin hydrogenase to shift their fermentation towards less energetically favorable pathways. Consistent with this explanation is the decreased abundance of E. rectale (Fig. 4), which had a slower growth rate and lower yield under higher H2 in monoculture (Supplementary Table 3). In the absence of a methanogen, R. intestinalis exhibited a similar, but less obvious, decrease in abundance (Fig. 4), although its growth rate and yield were not significantly impacted by high H2 in monoculture (Supplementary Table 3).

Fig. 4
figure 4

Influence of methanogenesis on competitive fitness of butyrogens in a synthetic gut community. Fold change in relative abundance of F. prausnitzii, E. rectale, and R. intestinalis in a 9-species synthetic human gut community with the addition of M. smithii compared to the same community without M. smithii. Relative abundance was quantified at the end of three successive 24-h subcultures (first subculture in light gray, second in dark gray, third in black), with a new inoculum of M. smithii added to the appropriate cultures at each passage. For each species, there were five replicate cultures in each condition. Statistical significance calculated from relative abundance values using two-sided Student’s two-sample t-tests (*p < 0.05; **p < 0.01; ***p < 0.001)

Incubation of the synthetic community cultures with vigorous shaking completely abrogated the effect of M. smithii on butyrate production and the relative abundances of butyrate producers (Supplementary Fig. 3C, I). This may be because local accumulation of dissolved H2 is prevented when H2 in the culture medium is rapidly equilibrated with the headspace, preventing H2 consumption by M. smithii from making a difference by reducing this accumulation. Shaking entirely prevented net production of lactate, which also suggests decreased exposure of the butyrogens to high H2 (Supplementary Fig. 3).

To explore whether the influence of M. smithii on butyrate production is relevant in the human gut, we collected breath samples from a human cohort both before and during consumption of resistant starch from potatoes (RSP). Resistant starch is not degraded by human amylases and reaches the gut microbiota undigested, where it can serve as a substrate for fermentation. We previously reported that RSP supplementation in a portion of this cohort increased fecal butyrate overall [45]. In this study, measurement of breath methane was used to assess gut methanogenesis both before and during RSP consumption. During RSP consumption, active gut methanogenesis was associated with lower levels of breath H2 (Fig. 3E) and lower fecal butyrate concentration (Fig. 3F) compared to individuals lacking gut methanogenesis. These findings are consistent with the results from in vitro cultures and suggest that the H2 produced from RSP breakdown in the gut may play an important role in the stimulation of butyrate production by the microbiota. Indeed, individuals with gut methanogenesis did not follow overall trend of increased fecal butyrate during RSP consumption (Supplementary Fig. 4D). Removal of H2 by hydrogenotrophs such as M. smithii appears to shift butyrogen metabolism in vivo as well as in vitro. Interestingly, gut methanogenesis in the same individuals before RSP consumption was not associated with decreased breath H2 or fecal butyrate (Supplementary Fig. 4A–C).

Discussion

The in vitro results reported in this study revealed shifts in fermentation products of human gut butyrogens grown under a headspace containing 1 atm partial pressure of H2. While this quantity of H2 is not found in intestinal gas, the H2 concentration relevant for microbial physiology is not that of the gas above a microbial culture, but rather that of the H2 dissolved in the aqueous phase where microbes dwell [41, 54]. Studies in bioreactors have shown that H2-producing microbial communities experience dissolved H2 concentrations many times greater than equilibrium with the headspace, with reports indicating 3- to 100-fold overconcentrations in various conditions and bioreactor designs [54,55,56,57]. Since H2 in intestinal gas (analogous to bioreactor headspace gas) ranges from < 1% to > 40% (v/v) with a median of approximately 15% (Supplementary Fig. 1) [8, 9], it is likely that dissolved H2 in the human colon ranges above and below that produced by equilibration with 1 atm H2. Therefore, the metabolic shifts induced in human gut butyrogens by the 1 atm H2 headspace used in our in vitro cultures could also occur in vivo. Observations from a human cohort were consistent with this hypothesis (Fig. 3D–F).

Much of the H2 produced in the human gut is consumed in situ by hydrogenotrophic microbes [40, 41]. Accordingly, in this study, we investigated the consumption of H2 by gut methanogens, reasoning that active methanogenesis must at some level lead to a reduction in dissolved H2. The fact that nearly all CH4 in the human gut is produced by the single culturable species, Methanobrevibacter smithii [53, 58], allowed us to use a simple, in vitro synthetic gut community to model the influence of methanogenesis on fermentation in the human gut.

Previous studies of the role of gut methanogenesis have often focused on aspects of human health and sometimes produce conflicting results [44, 59]. Explaining these inconsistencies, and distinguishing between correlation and causation, is difficult without mechanistically founded expectations about the effect of H2 removal on the gut microbiota [59]. In the current study, we hoped to obtain more interpretable results by first studying pure cultures of important gut microbes to test theoretical expectations (Fig. 2). Establishing the effects of H2 concentration in this system allowed us to develop predictions for highly simplified synthetic gut communities in which H2 was modulated by M. smithii (as in the gut) rather than by direct experimental manipulation of the headspace gas (Fig. 3A–C). Finding a methanogen-mediated decrease in butyrate production in this system in turn allowed us to understand the observation of lower fecal butyrate in methanogenic individuals (previously reported by Abell et al. (2009) [60] in a small cohort of eight individuals) as consistent with the predicted effect of gut methanogens rather than simply an intriguing association (Fig. 3D–F). Notably, the increase in lactate production observed in E. rectale under high H2 likely also drives increased fecal butyrate given that lactate in the human colon appears to be rapidly fermented to SCFAs including butyrate [61,62,63]. Certain human gut butyrogens, notably Anaerostipes caccae and Eubacterium hallii appear to specialize in this route of butyrate production when lactate is available, while R. intestinalis, E. rectale, and F. prausnitzii have not been observed to significantly utilize lactate as a substrate [64, 65]. Our synthetic gut community included A. caccae and so provided an in vitro model of the process.

A previous study also reported decreased butyrate and increased acetate production by human gut butyrogens in the prevalent but low-abundnace genus Christensenella in in vitro co-culture with M. smithii [66]. The reported shift in fermentation was similar to our findings in R. intestinalis, indicating that the effects of H2 concentration and methanogenesis we describe are common to other human gut butyrogens beyond the strains we investigated. Another study failed to find any influence of co-culture with M. smithii on R. intestinalis fermentation and found that co-culture with the hydrogenotrophic acetogen Blautia hydrogenotrophica actually increased butyrate production [67]. However, these results were driven by acetate availability, as acetate was not provided in the culture medium and net acetate consumption is required for production of high levels of butyrate (Fig. 1). This likely does not reflect the environment of the human colon, where acetate is abundant [68]. The most direct evidence to date of M. smithii modulating fermentation in vivo does not involve a butyrogen, but rather the commonly studied Bacteroides thetaiotaomicron. A study using gnotobiotic mice showed that M. smithii modulates B. theta fermentation products in vivo, increasing acetate and formate production at the expense of propionate, which the authors interpreted as due to consumption of H2 and/or formate by M. smithii [69].

In the present study, simple in vitro experiments with single species allowed a more specific description of the influence of H2 removal on human gut fermentation beyond the commonly repeated broad description of it as facilitating, enhancing, or improving the efficiency of human gut fermentation on the whole [44, 52, 70,71,72,73]. The principle that H2 removal facilitates H2-producing fermentation in the human gut is well-founded and accounts for the decrease in competitive fitness of the hydrogenogenic butyrogens E. rectale and (marginally) R. intestinalis in the synthetic community experiments reported here (Fig. 4), as well as the impairment of E. rectale growth rate and yield under very high H2 (Supplementary Table 3). However, this perspective obscures the fact that H2 accumulation does not simply shut down fermentation in the human gut, as it does in other well-studied systems such as sewage digesters. There, endergonic oxidation of butyrate and propionate to acetate requires an intimate syntrophic association between fermenters and methanogens [4, 31]. Our findings show that unlike these obligate syntrophs, the predominant human gut butyrogens E. rectale and R. intestinalis [13, 35] can cope with elevated H2 by disposing of reducing equivalents via butyrate and lactate instead. They do suffer some loss of metabolic efficiency, especially in the case of E. rectale which forgoes roughly half of the its ATP formation per glucose with its dramatic shift from butyrate and acetate production to lactate fermentation. However, they continue to grow using “backup” metabolic strategies. They are therefore examples of “facultative syntrophs” [5] for whom H2 accumulation results in a fermentation shift rather than a fermentation arrest. Counterintuitively, high H2 concentration actually stimulates production of the fermentation products butyrate and lactate in these organisms. As predicted by estimates of the Gibbs free energy of a range of fermentation balances (Fig. 1B), exposure to increasing concentrations of H2 shifted fermentation products in a roughly dose–response fashion, showing that the shift is progressive and not governed by a fixed H2 threshold (Supplementary Fig. 2).

The reduced fecal butyrate we report in methanogenic individuals only appears during consumption of an RSP supplement (Fig. 3F) and is not observed in the same individuals before supplement consumption (Supplementary Fig. 4A–C). The most likely explanation of this result is that RSP consumption is necessary in most individuals to stimulate sufficient production of H2 in the colon to change the thermodynamic situation if not efficiently removed. This possibility is supported by higher average H2 during versus before RSP consumption (p = 0.005). Another explanation, not mutually exclusive with the first, is based on the biogeography of methanogens in the human colon. A number of reports indicate that methanogens are more abundant in the distal colon and rectum rather than the proximal colon [43, 74,75,76]. As a refractive substrate, RSP may reach the distal colon in higher quantity than most other substrates in the diet before supplementation. Therefore, RSP fermentation could be more influenced by methanogens than fermentation of substrates that are mostly degraded before reaching the distal colon. Other guilds of human gut hydrogenotrophs—the sulfate reducers and reductive acetogens—may play a greater role in modulating fermentation of these substates. Further work should seek to include these guilds of hydrogenotrophs in our understanding of the role of H2 in modulating fermentation in the gut microbiome.

A final point of consideration is the negative association between active gut methanogenesis and successful stimulation of butyrate by RSP supplementation. While RSP supplementation generally increased fecal butyrate [45], methanogenic individuals showed no increase in fecal butyrate on average (Supplementary Fig. 4D). Although this is a correlative finding, this study provides a theoretical basis for a causal role of methanogenesis in decreasing butyrate production via efficient H2 removal. Given the myriad positive effects of butyrate on colon health [14], consideration should be given to reducing methanogenesis (and perhaps hydrogenotrophy in general) during supplement interventions intended to stimulate butyrate production. An alternative approach would be to administer H2 to stimulate butyrate production directly. A large body of research has studied H2 administration for its apparent antioxidant and anticarcinogenic effects, often via the consumption of water supersaturated with H2 [77,78,79]. Our findings in this study raise the possibility that these treatments may also stimulate butyrate production in the gut microbiota, especially in combination with supplement interventions.

Conclusions

H2 has often been proposed as a regulator of metabolic processes in the human gut microbiota [44], but concrete information is lacking on its specific role in the complex gut ecosystem. In this study, we examined the effect of H2 concentration on one prominent aspect of the human gut microbiota: production of the anti-inflammatory and anti-carcinogenic bacterial metabolite butyrate. Using in vitro approaches, we were able to observe the effect of H2 on three prominent human gut butyrogens: R. intestinalis, E. rectale, and F. prausnitzii. We found that high H2 concentration upregulated butyrate production by R. intestinalis, but not in E. rectale, which instead upregulated lactate production. F. prausnitzii was unaffected by H2. We further found that H2 consumption by the predominant gut methanogen M. smithii was sufficient to alter butyrate production by the H2-regulated butyrogens. Findings from a large human cohort supported a model in which gut H2 concentration, which is a balance between H2 production by fermenting bacteria and H2 consumption by methanogens, influences the total butyrate production by the gut microbiota.

Availability of data and materials

16S V4 amplicon sequences generated from synthetic community cultures grown in the course of this study have been deposited in the NCBI Sequence Read Archive (SRA) as BioProject ID PRJNA956530 (https://www.ncbi.nlm.nih.gov/sra/PRJNA956530). 16S amplicon data was analyzed using mothur v1.39.5 [47]. The mothur script used and logfile generated have been deposited in Zenodo (https://doi.org/10.5281/zenodo.6621661). All other data supporting the conclusions of this study is available in the paper and supplemental materials. In particular, the de-identified human cohort data used in this study is available in Supplemental Table 2. Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

References

  1. Schwartz E, Fritsch J, Friedrich B. H2-metabolizing prokaryotes. In: The Prokaryotes. Berlin: Springer Berlin Heidelberg; 2013. p. 119–99.

  2. Thauer RK, Jungermann K, Decker K. Energy conservation in chemotrophic anaerobic bacteria. Bacteriol Rev. 1977;41:100–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Gottschalk G. Bacterial metabolism. New York: Springer New York; 1986.

  4. Stams AJM. Metabolic interactions between anaerobic bacteria in methanogenic environments. Antonie Van Leeuwenhoek. 1994;66:271–94.

    Article  CAS  PubMed  Google Scholar 

  5. Stams AJM, Plugge CM. Electron transfer in syntrophic communities of anaerobic bacteria and archaea. Nat Rev Microbiol. 2009;7:568–77.

    Article  CAS  PubMed  Google Scholar 

  6. Macfarlane S, Macfarlane GT. Regulation of short-chain fatty acid production. Proc Nutr Soc. 2003;62:67–72.

    Article  CAS  PubMed  Google Scholar 

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

  8. Kirk E. The quantity and composition of human colonic flatus. Gastroenterology. 1949;12:782–94.

    Article  CAS  PubMed  Google Scholar 

  9. Levitt MD. Volume and composition of human intestinal gas determined by means of an intestinal washout technic. N Engl J Med. 1971;284:1394–8.

  10. Sonnenburg ED, Sonnenburg JL. Starving our microbial self: the deleterious consequences of a diet deficient in microbiota-accessible carbohydrates. Cell Metab. 2014;20:779–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Muir JG, Lu ZX, Young GP, Cameron-Smith D, Collier GR, O’Dea K. Resistant starch in the diet increases breath hydrogen and serum acetate in human subjects. Am J Clin Nutr. 1995;61:792–9.

    Article  CAS  PubMed  Google Scholar 

  12. Oku T, Nakamura S. Comparison of digestibility and breath hydrogen gas excretion of fructo-oligosaccharide, galactosyl-sucrose, and isomalto-oligosaccharide in healthy human subjects. Eur J Clin Nutr. 2003;57:1150–6.

    Article  CAS  PubMed  Google Scholar 

  13. Louis P, Flint HJ. Diversity, metabolism and microbial ecology of butyrate-producing bacteria from the human large intestine. FEMS Microbiol Lett. 2009;294:1–8.

    Article  CAS  PubMed  Google Scholar 

  14. Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost FJ, Brummer RJ. Review article: the role of butyrate on colonic function. Aliment Pharmacol Ther. 2008;27:104–19.

  15. Lee C, Kim BG, Kim JH, Chun J, Im JP, Kim JS. Sodium butyrate inhibits the NF-kappa B signaling pathway and histone deacetylation, and attenuates experimental colitis in an IL-10 independent manner. Int Immunopharmacol. 2017;51:47–56.

    Article  CAS  PubMed  Google Scholar 

  16. Weaver GA, Krause JA, Miller TL, Wolin MJ. Short chain fatty acid distributions of enema samples from a sigmoidoscopy population: an association of high acetate and low butyrate ratios with adenomatous polyps and colon cancer. Gut. 1988;29:1539–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. D’Argenio G, Cosenza V, Cave MD, Iovino P, Valle ND, Lombardi G, et al. Butyrate enemas in experimental colitis and protection against large bowel cancer in a rat model. Gastroenterology. 1996;110:1727–34.

    Article  PubMed  Google Scholar 

  18. Louis P, Flint HJ. Formation of propionate and butyrate by the human colonic microbiota. Environ Microbiol. 2017;19:29–41.

    Article  CAS  PubMed  Google Scholar 

  19. Buckel W, Thauer RK. Flavin-based electron bifurcation, ferredoxin, flavodoxin, and anaerobic respiration with protons (Ech) or NAD+ (Rnf) as electron acceptors: a historical review. Front Microbiol. 2018;9:401.

  20. Hillman ET, Kozik AJ, Hooker CA, Burnett JL, Heo Y, Kiesel VA, et al. Comparative genomics of the genus Roseburia reveals divergent biosynthetic pathways that may influence colonic competition among species. Microb Genom. 2020;6:mgen000399.

  21. Ragsdale SW. Pyruvate ferredoxin oxidoreductase and its radical intermediate. Chem Rev. 2003;103:2333–46.

    Article  CAS  PubMed  Google Scholar 

  22. Charon M, Volbeda A, Chabrière E, Pieulle L, Hatchikian EC, Fontecilla-Camps J. Crystal structures of the key anaerobic enzyme pyruvate:ferredoxin oxidoreductase, free and in complex with pyruvate. Nat Struct Biol. 1999;6:182–90.

    Article  PubMed  Google Scholar 

  23. Thauer RK, Kirchniawy FH, Jungermann KA. Properties and function of the pyruvate-formate-lyase reaction in clostridiae. Eur J Biochem. 1972;27:282–90.

    Article  CAS  PubMed  Google Scholar 

  24. Dostal A, Lacroix C, Bircher L, Pham VT, Follador R, Zimmermann MB, et al. Iron modulates butyrate production by a child gut microbiota in vitro. mBio. 2015;6:e01453.

  25. Louis P, Duncan SH, McCrae SI, Millar J, Jackson MS, Flint HJ. Restricted distribution of the butyrate kinase pathway among butyrate-producing bacteria from the human colon. J Bacteriol. 2004;186:2099–106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Bennett GN, Rudolph FB. The central metabolic pathway from acetyl-CoA to butyryl-CoA in Clostridium acetobutylicum. FEMS Microbiol Rev. 1995;17:241–9.

    Article  CAS  Google Scholar 

  27. Li F, Hinderberger J, Seedorf H, Zhang J, Buckel W, Thauer RK. Coupled ferredoxin and crotonyl coenzyme A (CoA) reduction with NADH catalyzed by the butyryl-CoA dehydrogenase/Etf complex from Clostridium kluyveri. J Bacteriol. 2008;190:843–50.

    Article  CAS  PubMed  Google Scholar 

  28. Chen J-S, Mortenson LE. Purification and properties of hydrogenase from Clostridium pasteurianum W5. Biochimica et Biophysica Acta (BBA) - Protein Structure. 1974;371:283–98.

  29. Biegel E, Müller V. Bacterial Na + -translocating ferredoxin:NAD + oxidoreductase. Proc Natl Acad Sci. 2010;107:18138–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Vital M, Howe AC, Tiedje JM, Moran MA. Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data. mBio. 2014;5:e00889–14.

  31. Angenent LT, Karim K, Al-Dahhan MH, Wrenn BA, Domíguez-Espinosa R. Production of bioenergy and biochemicals from industrial and agricultural wastewater. Trends Biotechnol. 2004;22:477–85.

    Article  CAS  PubMed  Google Scholar 

  32. Van Andel JG, Zoutberg GR, Crabbendam PM, Breure AM. Glucose fermentation by Clostridium butyricum grown under a self generated gas atmosphere in chemostat culture. Appl Microbiol Biotechnol. 1985;23:21–6.

    Article  Google Scholar 

  33. Pavlostathis SG, Miller TL, Wolin MJ. Cellulose fermentation by continuous cultures of Ruminococcus albus and Methanobrevibacter smithii. Appl Microbiol Biotechnol. 1990;33:109–16.

    Article  CAS  Google Scholar 

  34. Lamed RJ, Lobos JH, Su TM. Effects of stirring and hydrogen on fermentation products of Clostridium thermocellum. Appl Environ Microbiol. 1988;54:1216–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Duncan SH, Belenguer A, Holtrop G, Johnstone AM, Flint HJ, Lobley GE. Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate-producing bacteria in feces. Appl Environ Microbiol. 2007;73:1073–8.

    Article  CAS  PubMed  Google Scholar 

  36. D’hoe K, Vet S, Faust K, Moens F, Falony G, Gonze D, et al. Integrated culturing, modeling and transcriptomics uncovers complex interactions and emergent behavior in a three-species synthetic gut community. Elife. 2018;7:e37090.

  37. Doremus MG, Linden JC, Moreira AR. Agitation and pressure effects on acetone-butanol fermentation. Biotechnol Bioeng. 1985;27:852–60.

    Article  CAS  PubMed  Google Scholar 

  38. Kafkewitz D, Iannotti EL, Wolin MJ, Bryant MP. An anaerobic chemostat that permits the collection and measurement of fermentation gases. Appl Microbiol. 1973;25:612–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Weimer PJ, Zeikus JG. Fermentation of cellulose and cellobiose by Clostridium thermocellum in the absence of Methanobacterium thermoautotrophicum. Appl Environ Microbiol. 1977;33:289–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Christl SU, Murgatroyd PR, Gibson GR, Cummings JH. Production, metabolism, and excretion of hydrogen in the large intestine. Gastroenterology. 1992;102:1269–77.

    Article  CAS  PubMed  Google Scholar 

  41. Strocchi A, Levitt MD. Factors affecting hydrogen production and consumption by human fecal flora. The critical roles of hydrogen tension and methanogenesis. J Clin Invest. 1992;89:1304–11.

  42. Gibson GR, Macfarlane S, Macfarlane GT. Metabolic interactions involving sulphate-reducing and methanogenic bacteria in the human large intestine. FEMS Microbiol Ecol. 1993;12:117–25.

    Article  CAS  Google Scholar 

  43. Nava GM, Carbonero F, Croix JA, Greenberg E, Gaskins HR. Abundance and diversity of mucosa-associated hydrogenotrophic microbes in the healthy human colon. ISME J. 2012;6:57–70.

    Article  CAS  PubMed  Google Scholar 

  44. Carbonero F, Benefiel AC, Gaskins HR. Contributions of the microbial hydrogen economy to colonic homeostasis. Nat Rev Gastroenterol Hepatol. 2012;9:504–18.

    Article  CAS  PubMed  Google Scholar 

  45. Baxter NT, Schmidt AW, Venkataraman A, Kim KS, Waldron C, Schmidt TM. Dynamics of human gut microbiota and short-chain fatty acids in response to dietary interventions with three fermentable fibers. mBio. 2019;10:e02566.

  46. Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Tukey JW. Exploratory Data Analysis. Reading: Addison-Wesley; 1977.

    Google Scholar 

  49. Gottlieb K, Le C, Wacher V, Sliman J, Cruz C, Porter T, et al. Selection of a cut-off for high- and low-methane producers using a spot-methane breath test: results from a large north American dataset of hydrogen, methane and carbon dioxide measurements in breath. Gastroenterol Rep (Oxf). 2017;doi: 10.1093/gastro/gow048.

  50. Dabrock B, Bahl H, Gottschalk G. Parameters affecting solvent production by Clostridium pasteurianum. Appl Environ Microbiol. 1992;58:1233–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Duncan SH, Hold GL, Harmsen HJM, Stewart CS, Flint HJ. Growth requirements and fermentation products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii gen. nov., comb. nov. Int J Syst Evol Microbiol. 2002;52:2141–6.

  52. Samuel BS, Hansen EE, Manchester JK, Coutinho PM, Henrissat B, Fulton R, et al. Genomic and metabolic adaptations of Methanobrevibacter smithii to the human gut. Proc Natl Acad Sci. 2007;104:10643–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Scanlan PD, Shanahan F, Marchesi JR. Human methanogen diversity and incidence in healthy and diseased colonic groups using mcrA gene analysis. BMC Microbiol. 2008;8:79.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Kraemer JT, Bagley DM. Supersaturation of dissolved H2 and CO2 during fermentative hydrogen production with N2 sparging. Biotechnol Lett. 2006;28:1485–91.

    Article  CAS  PubMed  Google Scholar 

  55. Pauss A, Samson R, Guiot S, Beauchemin C. Continuous measurement of dissolved H2 in an anaerobic reactor using a new hydrogen/air fuel cell detector. Biotechnol Bioeng. 1990;35:492–501.

    Article  CAS  PubMed  Google Scholar 

  56. Pauss A, Guiot SR. Hydrogen monitoring in anaerobic sludge bed reactors at various hydraulic regimes and loading rates. Water Environ Res. 1993;65:276–80.

    Article  CAS  Google Scholar 

  57. Pauss A, Andre G, Perrier M, Guiot SR. Liquid-to-gas mass transfer in anaerobic processes: inevitable transfer limitations of methane and hydrogen in the biomethanation process. Appl Environ Microbiol. 1990;56:1636–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Abu-Ali GS, Mehta RS, Lloyd-Price J, Mallick H, Branck T, Ivey KL, et al. Metatranscriptome of human faecal microbial communities in a cohort of adult men. Nat Microbiol. 2018;3:356–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Levitt MD, Furne JK, Kuskowski M, Ruddy J. Stability of human methanogenic flora over 35 years and a review of insights obtained from breath methane measurements. Clin Gastroenterol Hepatol. 2006;4:123–9.

    Article  CAS  PubMed  Google Scholar 

  60. Abell GCJ, Conlon MA, McOrist AL. Methanogenic archaea in adult human faecal samples are inversely related to butyrate concentration. Microb Ecol Health Dis. 2009;18:154–60.

    Google Scholar 

  61. Bourriaud C, Robins RJ, Martin L, Kozlowski F, Tenailleau E, Cherbut C, et al. Lactate is mainly fermented to butyrate by human intestinal microfloras but inter-individual variation is evident. J Appl Microbiol. 2005;99:201–12.

    Article  CAS  PubMed  Google Scholar 

  62. Belenguer A, Duncan SH, Holtrop G, Anderson SE, Lobley GE, Flint HJ. Impact of pH on lactate formation and utilization by human fecal microbial communities. Appl Environ Microbiol. 2007;73:6526–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Wang SP, Rubio LA, Duncan SH, Donachie GE, Holtrop G, Lo G, et al. Pivotal roles for pH, lactate, and lactate-utilizing bacteria in the stability of a human colonic microbial ecosystem. mSystems. 2020;5:e00645.

  64. Duncan SH, Louis P, Flint HJ. Lactate-utilizing bacteria, isolated from human feces, that produce butyrate as a major fermentation product. Appl Environ Microbiol. 2004;70:5810–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Falony G, Verschaeren A, De Bruycker F, De Preter V, Verbeke K, Leroy F, et al. In vitro kinetics of prebiotic inulin-type fructan fermentation by butyrate-producing colon bacteria: implementation of online gas chromatography for quantitative analysis of carbon dioxide and hydrogen gas production. Appl Environ Microbiol. 2009;75:5884.

  66. Ruaud A, Esquivel-Elizondo S, de la Cuesta-Zuluaga J, Waters JL, Angenent LT, Youngblut ND, et al. Syntrophy via interspecies H2 Transfer between Christensenella and Methanobrevibacter underlies their global cooccurrence in the human gut. mBio. 2020;11:1022.

  67. Chassard C, Bernalier DA. H2 and acetate transfers during xylan fermentation between a butyrate-producing xylanolytic species and hydrogenotrophic microorganisms from the human gut. FEMS Microbiol Lett. 2006;254:116–22.

    Article  CAS  PubMed  Google Scholar 

  68. Cummings JH, Pomare EW, Branch WJ, Naylor CP, Macfarlane GT. Short chain fatty acids in human large intestine, portal, hepatic and venous blood. Gut. 1987;28:1221–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Samuel BS, Gordon JI. A humanized gnotobiotic mouse model of host–archaeal–bacterial mutualism. Proc Natl Acad Sci. 2006;103:10011–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. DiBaise JK, Zhang H, Crowell MD, Krajmalnik-Brown R, Decker GA, Rittmann BE. Gut microbiota and its possible relationship with obesity. Mayo Clin Proc. 2008;83:460–9.

    Article  PubMed  Google Scholar 

  71. Nakamura N, Lin HC, McSweeney CS, Mackie RI, Gaskins HR. Mechanisms of microbial hydrogen disposal in the human colon and implications for health and disease. 2010;1:363–95.

    CAS  Google Scholar 

  72. Smith NW, Shorten PR, Altermann EH, Roy NC, McNabb WC. Hydrogen cross-feeders of the human gastrointestinal tract. Gut Microbes. 2019;10:270–88.

    Article  CAS  PubMed  Google Scholar 

  73. Fischbach MA, Sonnenburg JL. Eating for two: how metabolism establishes interspecies interactions in the gut. Cell Host Microbe. 2011;10:336–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Flourie B, Etanchaud F, Florent C, Pellier P, Bouhnik Y, Rambaud JC. Comparative study of hydrogen and methane production in the human colon using caecal and faecal homogenates. Gut. 1990;31:684–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Pochart P, Lémann F, Flourié B, Pellier P, Goderel I, Rambaud J-C. Pyxigraphic sampling to enumerate methanogens and anaerobes in the right colon of healthy humans. Gastroenterology. 1993;105:1281–5.

    Article  CAS  PubMed  Google Scholar 

  76. Macfarlane GT, Gibson GR, Cummings JH. Comparison of fermentation reactions in different regions of the human colon. J Appl Bacteriol. 1992;72:57–64.

    CAS  PubMed  Google Scholar 

  77. Ohsawa I, Ishikawa M, Takahashi K, Watanabe M, Nishimaki K, Yamagata K, et al. Hydrogen acts as a therapeutic antioxidant by selectively reducing cytotoxic oxygen radicals. Nat Med. 2007;13:688–94.

    Article  CAS  PubMed  Google Scholar 

  78. Hu Y, Wang P, Han K. Hydrogen attenuated inflammation response and oxidative in hypoxic ischemic encephalopathy via Nrf2 mediated the inhibition of NLRP3 and NF-κB. Neuroscience. 2022;485:23–36.

    Article  CAS  PubMed  Google Scholar 

  79. Yang Y, Zhu Y, Xi X. Anti-inflammatory and antitumor action of hydrogen via reactive oxygen species. Oncol Lett. 2018;16:2771.

Download references

Acknowledgements

We would like to thank Kwi Kim for her HPLC analyses and help with Bifidobacterium isolation and Robert Hein for his work organizing human cohort data. We also thank Drs. Clegg Waldron and Ethan Hillman for comments and advice on the text of the manuscript. Figure 1A was created with BioRender.com.

Funding

This work was supported by the US National Institutes of Health (P01HL149633), the University of Michigan Mouse Metabolic Phenotyping Center (U2CDK110768), and the Howard Hughes Medical Institute.

Author information

Authors and Affiliations

Authors

Contributions

A.C. and T.S. designed the experiments and wrote the paper. A.C. performed the in vitro experiments and created figures. A.S. and K.G. collected and organized data from the human cohort. A.C. and K.G. analyzed human cohort data. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Thomas M. Schmidt.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of the University of Michigan Medical School (HUM00094242 and HUM00118951) and was conducted in compliance with the Helsinki Declaration. All human subjects providing samples analyzed in this study provided informed consent for sample collection and the use of data generated from samples in unspecified future research.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Rights and permissions

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Campbell, A., Gdanetz, K., Schmidt, A.W. et al. H2 generated by fermentation in the human gut microbiome influences metabolism and competitive fitness of gut butyrate producers. Microbiome 11, 133 (2023). https://doi.org/10.1186/s40168-023-01565-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-023-01565-3

Keywords