Antibiotic-induced changes in the microbiota disrupt redox dynamics in the gut

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

How host and microbial factors combine to structure gut microbial communities remains incompletely understood. Redox potential is an important environmental feature affected by both host and microbial actions. We assessed how antibiotics, which can impact host and microbial function, change redox state and how this contributes to post-antibiotic succession. We showed gut redox potential increased within hours of an antibiotic dose in mice. Host and microbial functioning changed under treatment, but shifts in redox potentials could be attributed specifically to bacterial suppression in a host-free ex vivo human gut microbiota model. Redox dynamics were linked to blooms of the bacterial family Enterobacteriaceae. Ecological succession to pre-treatment composition was associated with recovery of gut redox, but also required dispersal from unaffected gut communities. As bacterial competition for electron acceptors can be a key ecological factor structuring gut communities, these results support the potential for manipulating gut microbiota through managing bacterial respiration.

https://doi.org/10.7554/eLife.35987.001

eLife digest

The gut is home to a large and diverse community of bacteria and other microbes, known as the gut microbiota. The makeup of this community is important for the health of both the host and its residents. For instance, many gut bacteria help to digest food or keep disease-causing bacteria in check. In return, the host provides them with nutrients. When this balance is disturbed, the host is exposed to risks such as infections. In particular, treatments with antibiotics that kill gut bacteria can lead to side effects like diarrhea, because the gut becomes recolonized with harmful bacteria including Clostridium difficile and Salmonella.

Reese et al. have now investigated what happens to the gut environment after antibiotic treatment and how the gut microbiota recovers. Mice treated with broad-spectrum antibiotics showed an increase in the “redox potential” of their gut environment. Redox potential captures a number of measures of the chemical makeup of an environment, and provides an estimate for how efficiently some bacteria in that environment can grow. Some of the change in redox potential came from the host’s own immune system releasing chemicals as it reacted to the effects of the treatment. However, Reese et al. found that treating gut bacteria in an artificial gut – which has no immune system – also increased the redox potential. This experiment suggests that bacteria actively shape their chemical environment in the gut.

After the treatment, bacteria that thrive under high redox potentials, which include some disease-causing species, recovered first and fastest. This, in turn, helped to bring redox potential back to how it was before the treatment. Although the gut’s chemical environment recovered, some bacterial species were wiped out by the antibiotic treatment. The microbiota only returned to its previous state when the treated mice were housed together with non-treated mice. This was expected because mice that live together commonly exchange microbes, for instance by eating each other’s feces, and the treated mice received new species to replenish their microbiota.

These findings are important because they show that the chemical environment shapes and is shaped by the bacterial communities in the gut. Future research may investigate if altering redox potential in the gut could help to keep the microbiota healthier in infections and diseases of the digestive tract.

https://doi.org/10.7554/eLife.35987.002

Introduction

Mammalian gut microbial communities are likely to be structured by both host- and microbial-associated factors. Extensive research to date has focused on how host factors like diet (David et al., 2014; Claesson et al., 2012), genetics (Goodrich et al., 2014), geography (De Filippo et al., 2010; Yatsunenko et al., 2012; Arumugam et al., 2011), and immune state (Hooper et al., 2012) shape the gut microbiota. Yet, work in free-living microbial systems reveals that bacteria typically play active roles in shaping their own environment (Shi and Norton, 2000; Gobbetti, 1998; Goddard, 2008; Osono, 2005; Rui et al., 2009; Gerbersdorf et al., 2009). Identifying how these drivers interact is necessary both for a more complete understanding of the gut microbiota and for developing rational interventions. If ecological forces prove important, then compositional changes will likely be the result of feedbacks between the microbes and their environment.

Redox potential, a metric of the environmental capacity for reducing chemical reactions (i.e., those involving the gain of electrons) to occur, is a composite measurement of various factors that influence gut microbiota structure (Cowley et al., 2015; Friedman et al., 2017; Dhall et al., 2014). Much of our knowledge for how gut redox potentials are determined involves host-associated pathways (Spees et al., 2013; Rivera-Chávez et al., 2016). Passive diffusion of oxygen from the epithelium increases redox potential and stimulates growth of aerobic microbes (Espey, 2013; Albenberg et al., 2014). Secretion of redox-active immune molecules such as reactive oxygen species or nitrate is known to be a feature of inflammation that imposes oxidative stress on commensal microbes and can be exploited by select pathogens to colonize the intestine (Winter et al., 2010; 2013; Faber et al., 2016; Spees et al., 2013; Rivera-Chávez et al., 2016; Kelly et al., 2015; David et al., 2015). Maintenance of redox homeostasis in host tissue can also have spillover effects on luminal redox state (Circu and Aw, 2011).

Yet, redox potential is likely to also be shaped by microbial metabolism. In free-living microbial communities, variation in available electron acceptors dictates where microbes that employ respiration can thrive; the differential microbial metabolism that follows can produce further changes in electron acceptor availability (Morris and Schmidt, 2013; Chen et al., 2017; Noll et al., 2005; Orcutt et al., 2011). Moreover, microbial metabolism has been proposed as a mechanism for low redox potential states in the lung of cystic fibrosis patients (Cowley et al., 2015).

Here, we investigated the nature of redox potential dynamics under antibiotic treatment to assess the importance of host and microbial processes in structuring gut bacterial communities. Antibiotics directly disturb the microbiota but are also expected to alter host biology related to redox potential. Specifically, antibiotics have been found to increase gut epithelium oxygenation as a result of altered microbial composition and metabolic signaling to the host (Kelly et al., 2015; Rivera-Chávez et al., 2016). Increases in luminal oxygen due to diffusion from the microvasculature supplying the epithelium would lead to a higher redox potential under antibiotic treatment. In addition, host inflammation responses to antibiotic treatment and antibiotic-associated pathogen colonization have been shown to produce electron acceptors and other redox active molecules that cause oxidative stress (Faber et al., 2016; Winter et al., 2010; 2013; Spees et al., 2013). While these pathways have been demonstrated previously, their overall impact on redox potential has not been measured. Furthermore, the contribution of microbial metabolism to gut redox potential has not been tested, although it is known that a wide range of resident gut bacteria can respire aerobically and anaerobically (Ravcheev and Thiele, 2014). Antibiotic-driven bacterial inhibition could increase the availability of electron acceptors and thus serve as another mechanism for antibiotic induced changes in gut redox potential. On the other hand, antibiotic inhibition could limit bacterial production of oxidizing agents thereby resulting in an overall decrease in redox potential under antibiotics.

Changes in redox potential under antibiotics would in turn be expected to yield insight into the forces structuring the composition and function of the microbiota. Elevated redox potential due to a host immune response could restrict the microbiota beyond direct antibiotic mortality as inflammation would introduce additional oxidative stress into microbial ecosystems. By contrast, elevated redox potential due to the accumulation of oxygen or anaerobic electron acceptors would foster the growth of respiring bacteria. Antibiotic disturbance produces reproducible community succession in the gut following treatment—most notably, a transient bloom in Enterobacteriaceae (Antonopoulos et al., 2009; David et al., 2015; Young and Schmidt, 2004; Theriot et al., 2014; Peterfreund et al., 2012; Dethlefsen et al., 2008; Jakobsson et al., 2010; Looft and Allen, 2012), which has been ascribed to an increase in oxygen availability (David et al., 2015). But, the role of redox potential during this successional process has not previously been studied. If redox potential changes due to antibiotic treatment, we would predict that redox state recovery would be necessary for community resilience, that is for compositional recovery to a pre-disturbance state to occur (Shade et al., 2012). Furthermore, we expect that feedback between the community (i.e. the biotic component of the system) and their environment (i.e., the abiotic component) would drive further redox potential changes during the successional period. Such a pattern would highlight the potential of manipulating redox potential to alter community dynamics after disturbance.

Here, we combined in vivo and ex vivo antibiotic studies to isolate the effects of host and microbial pathways on redox state in the mammalian gut. We confirmed that redox potential increased during antibiotic treatment in association with some changes in host immune state. However, multiple lines of evidence in both mice and an artificial human gut model suggested that antibiotic-induced changes in microbial metabolism were sufficient to cause an increase in redox potential. After antibiotic treatments ended, we observed redox recovery within a week. A successional return to conventional community composition only occurred, though, when mice were co-housed and shared gut microbiota, indicating that microbial dispersal is necessary above and beyond environmental recovery for the return of normal microbiota community structure.

Results

Antibiotics caused a significant increase in mouse fecal redox potential within a day of treatment

For five days, we orally gavaged a cohort of conventional mice with a cocktail of antibiotics (ampicillin, vancomycin, metronidazole, and neomycin [Reikvam et al., 2011]) to broadly inhibit gut bacteria. We measured redox potential in freshly voided feces with a microelectrode paired with a reference electrode daily. Within sixteen hours after the first dose of antibiotics, redox potential significantly increased from 37 ± 164 mV at baseline to 227 ± 45 mV (p=0.04 Bonferroni-corrected Mann-Whitney U test; Figure 1A). Throughout treatment, redox potential differed overall between treated and control mice (p=0.005 linear mixed effects model likelihood tests).

Figure 1 with 2 supplements see all

Download asset Open asset

Antibiotic treatment effects on gut redox state and host inflammation.

(A) Redox state measured in freshly voided feces of treated (red) and control (black) mice (n = 9–10 per treatment) differed (p=0.005, linear mixed effects model likelihood tests). Replicate data is presented in Figure 1—figure supplement 2. (B) Gene expression for four inflammation associated genes measured with RT-qPCR of RNA isolated from feces on the final day of antibiotic treatment (n = 6–7) differed from control levels (p<0.05, Bonferroni-corrected one-sample t-tests). (C) Intestinal inflammation measured as fecal concentration of the biomarker lipocalin-2 (n = 9–10 per treatment) differed between treated and control mice (p<0.001 linear mixed effects model likelihood tests). Data are shown as means ± SD. Post-hoc test results for individual time points (* indicates p<0.05 Bonferroni-corrected Mann-Whitney U test) are included for visualization purposes. Antibiotic treatment began after Day 0 measurement; red bars indicate treatment duration.

https://doi.org/10.7554/eLife.35987.003

Evidence antibiotics alter host biology to increase redox potential

We examined first whether antibiotics affected redox potential through direct host effects mediated by an immune response. Expression of three host genes in fecal samples were consistent with the hypothesis that redox shifts were associated with immune activation. We measured increases relative to controls in Nos2, which is linked to reactive nitrogen species levels (Dedon and Tannenbaum, 2004; Winterbourn, 2008); increases relative to controls in Rela, which has been found to contribute to IBD type inflammation in mice via a proinflammatory cytokine response (Waddell et al., 2013); and, decreases relative to controls in Apoa4, which has known anti-inflammatory function (Broedl et al., 2007) (p<0.05, Bonferroni-corrected one-sample t-tests; Figure 1B). Yet, other biomarkers did not associate antibiotic treatment with an immune response. Nfkb1 expression, which is associated with inflammation suppression (Cartwright et al., 2016), increased after antibiotic treatment. Antibiotic treatment was also followed by a small, but significant, decrease in lipocalin-2 levels, which is a protein biomarker of inflammation (Chassaing et al., 2012) (p<0.001, linear mixed effects model likelihood tests; Figure 1C). Taken together, our biomarker assays provided equivocal evidence for intestinal inflammation in antibiotic-treated mice.

Bacterial responses to antibiotics associated with redox potential shift

To test if inhibition of bacterial populations by antibiotics could directly contribute to redox potential dynamics, we used an ex vivo human gut system based on a continuous-flow bioreactor (McDonald et al., 2013), which allowed the propagation of a stable microbial community representative of the human gut microbiota with all major phyla represented (Figure 2—figure supplement 1). Treating this system with the same antibiotic cocktail as used in the mouse study led to an increase in redox potential relative to an untreated control (p=0.005, linear mixed effects model likelihood tests; Figure 2A). Redox potential increased by 59 ± 47 mV within fifteen hours of the first antibiotic dose, mirroring the rapidity with which redox shifts occurred in vivo (Figure 1A). Redox potential in our ex vivo model increased again by another 141 ± 37 mV after a second antibiotic dose (Figure 2—figure supplement 1). Thus, in the absence of direct interactions between antibiotics and a host, antibiotic treatment produced shifts in environmental redox potential experienced by a gut microbial community.

Figure 2 with 3 supplements see all

Download asset Open asset

Microbial response to antibiotics ex vivo and in vivo.

(A) Redox potential from control (black) and antibiotic treated (red) bioreactor vessels culturing human gut microbial communities (n = 3 per treatment) differed (p=0.005, linear mixed effects model likelihood tests). (B) Microbial load, measured as 16S rRNA gene copy number in antibiotic-treated and control mice decreased during the first 24 hr of treatment for antibiotic treated mice (n = 9–10 per treatment; p<0.001, Bonferroni-corrected Mann-Whitney U tests). (C) Metabolites measured with NMR spectrometry of feces (n = 9–10 per treatment; see Figure 2—figure supplement 3 and Supplementary file 1) decreased in antibiotic treated mice during the first 24 hr (p<0.001, Bonferroni-corrected Mann-Whitney U tests). Data are shown as means ± SD. Post-hoc test results for individual time points (* indicates p<0.05 Bonferroni-corrected Mann-Whitney U test) are included for visualization purposes. Antibiotic treatment began after Day 0 measurement; red bars indicate treatment duration.

https://doi.org/10.7554/eLife.35987.006

Fecal redox potential under antibiotics is associated with bacterial metabolism and respiration

To understand how antibiotic effects on bacterial populations could lead to shifts in gut redox potential, we investigated the dynamics and metabolism of microbiota across treatment in mice. We observed depressed levels of bacterial load and metabolic activity that occurred within hours of initial redox potential shifts. Fecal bacterial concentrations decreased significantly within twelve hours of antibiotic treatment (p=0.01, Bonferroni-corrected Mann-Whitney U test; Figure 2B), and remained significantly below controls throughout the five days of treatment (p<0.001, Bonferroni-corrected Mann-Whitney U tests; Figure 2—figure supplement 2). Using NMR-based metabolomics, we found that the short-chain fatty acids propionate and acetate, end products of key microbial metabolic pathways, decreased eight hours after antibiotic treatment (p<0.05, Bonferroni-corrected Mann-Whitney U tests; Figure 2—figure supplement 3, Supplementary file 1). Twenty-one of twenty-eight metabolites measured, including other short-chain fatty acids, amino acids, branched-chain amino acids, and a group of bile acids, later decreased significantly under antibiotic treatment (p<0.05, Bonferroni-corrected Mann-Whitney U tests; Figure 2C, Figure 2—figure supplement 3, Supplementary file 1). The dynamics of twelve of those metabolites was significantly associated with treatment overall (p<0.05, linear mixed effects model likelihood tests; Supplementary file 1). Thus, the timing and persistence of decreased gut load and microbial activity coincided with increases in fecal redox potential.

Next, we measured the concentration of three major electron acceptors used during microbial respiration to assess whether changes in these pathways contributed to the redox potential increase. First, measuring luminal oxygen with a novel in vivo sensor system—a hydrogel sensor with covalently attached oxygen sensitive Pd-porphyrin derivative was inserted rectally then read optically through the skin—we observed a significant increase in luminal oxygen levels the day after antibiotic treatment (3.7 ± 3.6 Torr at baseline to 26.5 ± 26.1 Torr; p=0.02, Bonferroni-corrected Mann-Whitney U test; Figure 3A). Previous work (Kelly et al., 2015; Rivera-Chávez et al., 2016) has shown elevated tissue oxygenation under antibiotics, which could result in higher diffusion and therefore higher luminal oxygen. In addition, reduced aerobic respiration due to antibiotic mortality or stress could also increase gut oxygen levels. However, gut oxygen levels in antibiotic-treated mice returned to control levels by forty-eight hours after the first dose (p>0.05, Bonferroni-corrected Mann-Whitney U test), and antibiotics were not associated with an overall effect on gut oxygen levels (p>0.05 linear mixed effects model likelihood tests). Moreover, among treated mice, gut oxygen levels were not correlated with redox potential (p>0.05, repeated measures correlation; Figure 3—figure supplement 1). These data together suggest that increased luminal oxygenation during antibiotic treatment may contribute to early shifts in redox potential, but were not responsible for persistent redox shifts.

Figure 3 with 3 supplements see all

Download asset Open asset

Electron acceptor levels during antibiotic treatment.

(A) Large intestine luminal oxygen concentration measured in vivo for treated (red) and control (black) mice (n = 9–10 per treatment) did not differ overall during treatment (p>0.05, linear mixed effects model likelihood tests). Replicate data are presented in Figure 3—figure supplement 2. (B) Fecal nitrate concentration (n = 2–9 per group per time point) did differ between treated and control mice during treatment (p=0.02, linear mixed effects model likelihood tests). (C) Fecal fumarate concentration (n = 2–9 per group per time point) also differed (p<0.001, linear mixed effects model likelihood tests). Electron acceptor levels below the detection level are not plotted here. For fumarate measurements where all treated animals were below the detection limit, values were set to 200 nmol for statistical tests. Data are shown as means ± SD. Post-hoc test results for individual time points (* indicates p<0.05 Bonferroni-corrected Mann-Whitney U test) are included for visualization purposes. Antibiotic treatment began after Day 0 measurement; red bars indicate treatment duration.

https://doi.org/10.7554/eLife.35987.010

By contrast, we observed sustained and significant increases in nitrate levels during antibiotic treatment. Nitrate is one of the most widely used electron acceptors in the gut (Ravcheev and Thiele, 2014; Fischer and Lindenmayer, 2007) and its reduction to nitrogen gas (at pH 7) has a potential of +0.75 V making it one of the most potent electron acceptors after oxygen. The increase in Nos2 expression (Figure 1B) led us to hypothesize that nitrate would increase as it is a gene which encodes inducible nitric oxide synthase and whose expression has been found to confer a growth advantage to nitrate respiration competent E. coli strain (Winter et al., 2013). Lower levels of microbial respiration could also lead to an accumulation in nitrate separate from changes in host expression. We measured fecal nitrate levels with the chromotropic acid method (West and Ramachandran, 1966) for a subset of time points when redox varied, and found a significant increase in nitrate (p=0.03, Mann-Whitney U test; Figure 3B) beginning less than 24 hr after the first dose of antibiotics. This increase persisted throughout treatment, and overall antibiotic treatment produced a significant increase in nitrate (p=0.02, linear mixed effects model likelihood tests).

Finally, we used a colorimetric enzyme assay to measure fecal fumarate levels for a subset of time points during antibiotic treatment. Fumarate reductases have been found in over one third of gut microbe genomes, and fumarate is the most common terminal electron acceptor for bacterial anaerobic respiration (Kröger et al., 1992). Unlike nitrate and oxygen, though, much of the fumarate in the gut is likely produced by the microbiota itself (Fischbach and Sonnenburg, 2011; El Aidy et al., 2013). Overall, there was a significant effect of treatment on fumarate levels (p<0.001, linear mixed effects model likelihood tests; Figure 3C). We observed a spike in fumarate levels 4 hr after the first antibiotics dose (p=0.004, Bonferroni-corrected Mann-Whitney U test); however, this increase had dissipated by 8 hr after the first dose and later during treatment, we observed a significant decrease (to below the detection limit) in fumarate levels. Thus, not all electron acceptors increased in abundance during antibiotic treatment.

Community composition and gut environment recover following antibiotic treatment

Ecological succession, the ‘somewhat orderly and predictable’ (Fierer et al., 2010) dynamics of a community after a perturbation, is known to follow antibiotic treatment (Antonopoulos et al., 2009; David et al., 2015). What is less well-known is how abiotic conditions change during post-antibiotic succession and, by extension, how biotic and abiotic dynamics interface during that period. Having shown that both the community and gut environment are not resistant to antibiotic disturbance, we next sought to determine how resilient they were, i.e. how quickly they recovered (Shade et al., 2012; Allison and Martiny, 2008). In this study, we indeed observed consistent biotic changes over time in mice that received antibiotics (Figure 4—figure supplement 1). Gut bacterial load recovered rapidly after treatment across all treated mice, exhibiting no difference relative to control mice by two days (p>0.05, Bonferroni-corrected Mann-Whitney U tests; Figure 2—figure supplement 2). A reproducible recovery of beta-diversity was observed after one week (p>0.05, Bonferroni-corrected Mann-Whitney U test; Figure 4—figure supplement 2).

To characterize the recovery of specific microbial taxa, we clustered fecal bacterial genera into groups according to their dynamical patterns across antibiotic treatment (Figure 4A,B). We focused on the five clusters with an average abundance of at least 1% across our dataset (Supplementary file 2). Amongst these five groups, three differed significantly during treatment (p<0.05, linear mixed effects model likelihood tests; Figure 4A,B). Notably, the two clusters elevated during treatment included many facultative anaerobic taxa (Supplementary file 2). One of these, the cluster primarily composed of members of the Enterobacteriacaeae, remained elevated into the beginning of the recovery period as well (p<0.05, Bonferroni-corrected Mann-Whitney U tests). This enrichment is in line with previous findings that Enterobacteriaceae often bloom following antibiotic treatment (Antonopoulos et al., 2009; David et al., 2015; Young and Schmidt, 2004; Theriot et al., 2014; Peterfreund et al., 2012; Dethlefsen et al., 2008; Jakobsson et al., 2010; Looft and Allen, 2012). That cluster, as well as one composed primarily of Akkermansia and one of a consortium of typical commensal taxa, exhibited significant residual effects of treatment during the recovery period (p<0.05, linear mixed effects model likelihood tests; Figure 4A,B). However, by the end of the one week recovery period, the abundance of all clusters was indistinguishable between treated and control animals (p>0.05, Bonferroni-corrected Mann-Whitney U tests).

Figure 4 with 5 supplements see all

Download asset Open asset

Gut community is resilient to antibiotic disturbance.

(**A–B**) Compositional dynamics of most abundant clusters during and after treatment for control (A) and treated (B) mice (n = 9–10 per treatment). Three groups (Typical, Staph/Strep, Enterobacteriaceae) differed significantly between groups (p<0.05, linear mixed effects model likelihood tests) during treatment and two (Enterobacteriaceae, Akkermansia) during recovery. (C) Redox potential measured in freshly voided feces of singly housed mice during the recovery period (n = 5 per treatment) was lower than treated mice initially (p<0.05, Bonferroni-corrected Mann-Whitney U tests). (D) Repeated measures correlation was significant (p=0.02) between change in Enterobacteriaceae abundance and change in redox potential during the first three days of post-antibiotic recovery (n = 64). (E) Days until there is no significant difference (p>0.05 Bonferroni-corrected Mann Whitney U test) between groups of control and treated animals for redox and Enterobacteriaceae abundance, under either singly- housed or co-housed settings. Red bars indicate treatment duration (**A–B**). Red arrows indicate the last antibiotic dose; day six measurements are more than 24 hr after the last dose (C). Data are shown as means ± SD. Post-hoc test results for individual time points (* indicates p<0.05 Bonferroni-corrected Mann-Whitney U test) are included for visualization purposes.

https://doi.org/10.7554/eLife.35987.014

Abiotic conditions in the mouse gut were also resilient and recovery was predictable. Redox potential returned to control levels within a week after treatment ended (p>0.05, Bonferroni-corrected Mann-Whitney U tests; Figure 4C). Notably, though, we observed the day after antibiotic therapy ceased, fecal redox potentials in treated mice were significantly decreased relative to control mice (-138.3 ± 149.8 mV vs 54.6 ± 211 mV; p=0.006, Bonferroni-corrected Mann-Whitney U test), and it was not until later in recovery that redox potential returned to control levels across most mice. Oxygen and inflammation biomarker levels remained low throughout the recovery period (p>0.05, Bonferroni-corrected Mann-Whitney U tests; Figure 3—figure supplement 2, Figure 1—figure supplement 1, respectively) consistent with their minimal changes during treatment. Furthermore, many microbial metabolic products measured by NMR had recovered two days after the last doses, and all had recovered by the end of the recovery week (Figure 2—figure supplement 3; Supplementary file 1). In concert, a largely reproducible biotic and abiotic recovery took place across antibiotic treated mice in the days following antibiotic treatment indicating the overall resilience of this ecosystem to antibiotic disturbance.

Post-antibiotic succession resulted from interplay of biotic and abiotic recovery

In investigating potential interactions between abiotic and biotic processes in antibiotic recovery, we hypothesized that the Enterobacteriaceae could thrive under high redox potential conditions and directly impact gut redox state. Members of this bacterial family were enriched at the end of treatment in our clustering analysis (Figure 4B). Enterobacteriaceae can also employ many terminal electron acceptors to perform aerobic and anaerobic respiration (Ravcheev and Thiele, 2014) and thus could likely thrive under high redox potential conditions and directly impact gut redox state. We found that the relative abundance of the Enterobacteriaceae was marginally correlated with redox potential (p=0.1, r=,0.18, repeated measures correlation; Figure 4—figure supplement 3) throughout treatment and recovery. Moreover, the change in redox potential during the earlier recovery period was negatively correlated with changes in Enterobacteriaceae absolute abundance (calculated as total 16S rRNA gene copy number multiplied by relative abundance) for both treated animals (p=0.03, r = −0.42, repeated measures correlation) and for all animals (p=0.02, r = −0.34, repeated measures correlation; Figure 4D). Together, these observations suggest that increased respiration by the Enterobacteriaceae after antibiotic treatment may lead redox potentials to decrease below conventional levels. More broadly, our observations support the hypothesis that Enterobacteriaceae dynamics contribute to natural redox potential variation in the gut.

We next investigated if abiotic recovery led to compositional recovery or if dispersal (i.e., migration of microbes) from unimpacted microbial populations was necessary for microbiota to exhibit resilience and return to a pre-treatment state. Previous research has shown that singly housed mice treated with antibiotics will exhibit altered community composition for weeks after treatment (Antonopoulos et al., 2009). We therefore co-housed half of the treated mice with control animals during the recovery period and kept half in single housing. Because mice practice coprophagy (i.e., ingestion of feces), we would expect that co-housing would re-introduce normal commensal microbes to compete with the Enterobacteriaceae and other early successional species thereby reducing their abundance. At the end of the one-week treatment period, we indeed observed differential success in community recovery between the housing groups. Co-housed mice communities had returned to control composition while singly housed mice had not. For co-housed mice, by the end of the recovery period, the gut microbiota of treated animals was equally dissimilar from baseline as it was for control mice indicating that the community had recovered (p>0.05, Bonferroni-corrected Mann-Whitney U test; Figure 4—figure supplement 4). Both fecal redox potential and Enterobacteriaceae levels also returned to conventional levels more quickly in co-housed mice than singly housed mice (Figure 4E). By contrast, singly-housed treated mice continued to be more dissimilar at the end of recovery than singly-housed control mice (p<0.001, Bonferroni-corrected Mann-Whitney U test; Figure 4—figure supplement 4). Such drift is qualitatively consistent with arguments that stochastic variation can influence how communities respond to perturbation, as seen in soil microbial communities (Zhang et al., 2016). More broadly, the differences between singly- and co-housed mice here highlight that environmental recovery is not sufficient to produce biotic recovery—dispersal is also necessary for a resilient community.

Discussion

Here, we found that a key feature of the microbial environment, redox potential, can be rapidly and dramatically altered in the gut under antibiotic treatment. This shift was accompanied by changes in electron acceptor availability. The direction and timing of the redox potential shift was associated with microbial dynamics and could be reproduced with an ex vivo model of the human gut absent host effects. Moreover, we found evidence that redox potential kinetics were driven by shifts in microbial load and metabolism, specifically among facultative anaerobes from the family Enterobacteriaceae. Although there was evidence of increases in certain host pathways of electron acceptor production, we did not find persistent and significant changes in host inflammation. Therefore, we propose microbial respiration of host sourced electron acceptors acts as a major determinant of redox potential in the gut, in addition to the host production itself. Redox potential consequently shaped successional dynamics following the cessation of treatment, but return to control redox conditions was insufficient to produce community recovery. Dispersal was required for community structure to exhibit full resilience to antibiotic disturbance.

Our findings imply that electron acceptors are a critical resource normally competed over by the resident bacteria and, as such, are key factors structuring gut microbial ecology. The ability to respire a broad range of electron acceptors is already known to be crucial for pathogens to colonize the intestine (Winter et al., 2010; 2013; Faber et al., 2016). Yet, ecological research to date on commensal microbiota has primarily focused on substrate availability (Donaldson et al., 2016 and references therein, Pereira and Berry, 2017) or, to a lesser extent, oxygen (Espey, 2013). Our results here support the hypothesis that the availability of other electron acceptors is an important ecological factor driving the structure of resident bacterial populations as well (Jones et al., 2007; 2011). If true, this hypothesis would imply that pathogen colonization resistance by commensal microbes could be mediated by bacterial competition involving respiratory pathways.

Future work will be needed to comprehensively chart the dynamics and chemical landscape across which bacteria in the gut compete for electron acceptors. While any redox active molecule could contribute to redox potential changes (hence testing all possible chemical drivers of redox is infeasible), we focused here on three relevant electron acceptor species (oxygen, nitrate, and fumarate) that derive from independent pathways and may be representative of overall dynamics. Our data lead us to conjecture that oxygen availability transiently increased due to the decreased respiration capacity of an antibiotic-constrained bacterial community. But, as oxygen grew increasingly available, the remaining facultative anaerobes switched to aerobic respiration leaving other host-sourced electron acceptors, such as nitrate, to accumulate. Accrued electron acceptors, especially nitrate, may also have been contributed by changes in host physiology as evidenced by an upregulation of the Nos2 gene, which is linked to reactive nitrogen species levels (Dedon and Tannenbaum, 2004; Winterbourn, 2008). Still, not all electron acceptors accumulated during antibiotic treatment. Notably, our data suggest that ones produced by the microbiota rather than the host, like fumarate, were less available under antibiotic pressure.

Our data also support a model where consumption by specific bacterial taxa shape redox homeostasis in the gut. We observe that facultative anaerobes like the Enterobacteriaceae behave as pioneer taxa whose dynamics inversely track redox dynamics after antibiotic treatment, providing new understanding of the deterministic mechanism driving the consistent bloom of these bacteria following antibiotic treatment (Antonopoulos et al., 2009; David et al., 2015; Young and Schmidt, 2004; Theriot et al., 2014; Peterfreund et al., 2012; Dethlefsen et al., 2008; Jakobsson et al., 2010; Looft and Allen, 2012). Enterobacteriaceae are especially suited to exploiting electron acceptor availability because their diverse repertoire of respiration pathways (Ravcheev and Thiele, 2014; Jones et al., 2007; 2011) allow for fast growth rates under a high redox regime. Enterobacteriaceae and respiration pathway genes more generally are also more prevalent in inflamed mouse guts (Hughes et al., 2017), another environment where heightened redox potential is expected. Under experimental inflammation, though, Enterobacteriaceae expansion can be prevented by treatment with tungstate, a respiration inhibitor (Zhu et al., 2018). This result corroborates our model that members of the Enterobacteriaceae use respiration to bloom under disturbed conditions, and it highlights the potential for manipulating electron acceptor availability or usability to engineer the microbiota.

An increase in resource availability is common following ecological disturbance and can promote the growth of pioneer taxa, which are best suited to quickly respond to a resource spike but which draw down the resource thereby bringing about their own displacement (Connell and Slatyer, 1977; Peet and Christensen, 1980; Tilman, 1985). Here, as in many cases of ecological succession, the Enterobacteriaceae were eventually replaced by more conventional community members that typically grow better under more limiting conditions; in fact, the relationship between change in Enterobacteriaceae abundance and change in redox state is consistent with these taxa facilitating the success of secondary colonizers. However, why redox recovery was ultimately hysteretic (i.e., the gut temporarily became a reducing environment following the end of treatment) remains unclear. It has previously been shown that members of the Enterobacteriaceae will produce reducing agents, such as H₂S, in response to antibiotic treatment (Shatalin et al., 2011), but this theory requires further investigation in a community and host-associated context.

Beyond identifying a role for microbiota in shaping gut redox potentials, our results more broadly illustrate how microbial factors can be important in addition to host factors for determining abiotic conditions in the gut through the balancing of production and consumption of resources. While the host is ultimately the source of all substrates in the gut, microbial action determines the realized environment which shapes community composition and function. Such microbial ecosystem engineering (Jones et al., 1994; Wright and Jones, 2006) is consistent with our understanding of how microbes can shape their environments in free-living systems (Gerbersdorf et al., 2009; Goddard, 2008). The altered resource availabilities that result from microbial action then contribute to structuring the microbial community itself. Such feedbacks may contribute to a number of phenomena in host-associated systems. Common constituents of the gut (i.e., the phyla Firmicutes and Bacteroidetes) could perpetuate similar hospitable environments for themselves regardless of host biology leading to the broad, phylum-level consistency of composition seen between humans regardless of diet, geography, or genetics (David et al., 2014; Claesson et al., 2012; Goodrich et al., 2014; De Filippo et al., 2010; Yatsunenko et al., 2012; Arumugam et al., 2011). Variations in how the gut environment is dictated by microbes may also help explain inter-individual differences in response to interventions in the absence of consistent host effects (e.g., Venkataraman et al., 2016). Microbial metabolic signaling additionally can alter host production of substrates including redox active molecules. Most notably reduced butyrate levels, a change observed here, have been demonstrated to increase gut nitrate levels (Byndloss et al., 2017) and host epithelial oxygenation (Rivera-Chávez et al., 2016).

Conclusion

Together, our findings suggest new ecological models for how antibiotics reshape the gut microbiota and for how redox shifts could be associated with enteric disease. Antibiotics are triumphs of modern medicine that have dramatically reduced infectious disease mortality (Armstrong et al., 1999). But, we are increasingly learning that antibiotics also meaningfully reshape the resident gut microbiota, leaving an imprint that can last for months to years after treatment (Dethlefsen and Relman, 2011; Jakobsson et al., 2010) and predisposing hosts to obesity (Cho et al., 2012), food allergy (Stefka et al., 2014), auto-immune disease (Russell et al., 2012), and increased infection risk (Stecher et al., 2007; Buffie et al., 2012; Wiström et al., 2001). While these drugs reduce levels of susceptible organisms (Keeney et al., 2014), an additional ecological mechanism of action is decreasing microbial competition and allowing primary metabolites (e.g., primary bile acids, sugars) (Ng et al., 2013; Theriot et al., 2014), as well as host-sourced electron acceptors like oxygen and nitrate to accumulate. This concept complements recent discoveries that electron acceptors facilitate antibiotic-associated enteric pathogen colonization (Rivera-Chávez et al., 2016; Winter et al., 2013; 2010). Such increases in electron acceptor availability likely are not unique to antibiotic treatment and could generalize to various enteric disturbances. Indeed, germ-free animals (Phillips et al., 1958; Celesk et al., 1976) and humans suffering from inflammatory diseases (Circu and Aw, 2011) and malnutrition (Million et al., 2016) exhibit increased gut redox potential. Thus, novel treatments for microbial disorders or preventing antibiotic-associated infections may include chemical alterations of redox potential or introduction of competitors for excess electron acceptors. More broadly, we propose adding redox potential to the list of abiotic conditions frequently assayed and manipulated to improve host well-being.

Materials and methods

Animal experiments

Mice

All animal experiments were conducted in accordance with National Institute of Health Guide for the Care and Use of Laboratory Animals using protocols approved by the Duke University Institutional Animal Care and Use Committee. Male C57BL/6 mice (Charles River Laboratories, Wilmington MA) 8–10 weeks of age with a native microbiota were used for all experiments. Mice were kept in a conventional laboratory animal facility at Duke University. Baseline measurements and fecal samples were collected at least twenty-four hours before the first dose and then mice were placed in individual housing with supplementary enrichment.

Antibiotic treatment

Request a detailed protocol

Mice were orally gavaged with either 0.25 ml autoclaved deionized water (control, N = 10) or 0.25 ml of an antibiotic cocktail (treated, N = 10) daily for five days (Figure 5A). The mice were randomly assigned a group and researchers collecting data were blinded to the groupings until after the final dose was administered. The sample size was chosen following a power analysis to allow for β less than 0.1. The antibiotic cocktail consisted of ampicillin (Gold Biotechnology, St. Louis MO) 1 mg/ml, vancomycin (Alfa Aesar, Ward Hill MA) 5 mg/ml, neomycin (EMD Millipore, Burlington MA) 10 mg/ml, and metronidazole (Alfa Aesar) 10 mg/ml (after Reikvam et al., 2011). Fresh antibiotic cocktails were prepared every day. Throughout the experiment freshly voided fecal samples were collected and stored at −80°C for subsequent analysis.

Figure 5

Download asset Open asset

Experimental setup.

(A) Experimental sampling regime. (B) Cohousing treatment setup.

https://doi.org/10.7554/eLife.35987.020

Recovery cohousing

View detailed protocol

One day after the final gavage, mice were randomly assigned to cohousing groups (Figure 5B). Cohousing control cages (N = 8, four each from treated and control) contained single mice that were kept in individual housing throughout the recovery period. Cohousing treatment cages (N = 5) contained a control and an antibiotic-treated mouse from the same litter that were placed together in the control mouse cage and kept together throughout the recovery period. The sample size was determined by the number of mice available at the end of the treatment period. Antibiotic-treated mice in cohoused cages were marked by partial shaving to allow for continued sampling. Throughout recovery, freshly voided fecal samples were collected and stored at −80°C for subsequent analysis. While the antibiotic treatment experiment was performed multiple times, the cohousing during recovery experiment was only performed with a single replicate cohort.

Redox potential measurements

View detailed protocol

We measured redox potential in fresh fecal pellets using a redox electrode with a tip diameter of 500 µm in tandem with an Ag/AgCl reference electrode (Unisense, Aarhus Denmark). The values are given relative to the standard hydrogen electrode (SHE) and were determined by measuring the offset of the reference electrode in saturated quinhydrone buffer solutions (pH 4 and pH 7) with known redox potentials. Daily calibrations were done at room temperature; however, because of potential unmeasured variation in calibration procedures (e.g., fluctuation in ambient temperature), we focused our analyses on differences between control and treated animals at any given time point. The electrode has a detection limit of 0.10 mV (Pang and Zhang, 1998).

Fecal pellets were placed on an agar plate, resting against the reference electrode. The electrode was inserted into the pellet using a micromanipulator and data were collected continuously for three minutes. The data presented are the average across three minutes. Redox data was collected on two cohorts of mice with comparable results. Replicate data is shown in Figure 1—figure supplement 2.

Oxygen measurements

Request a detailed protocol

We measured oxygen concentration in vivo using a hydrogel oxygen microsensor coupled with an optical reader. The microsensors are composed of a biocompatible hydrogel, poly (2-hydroxyethyl methacrylate) (pHEMA) and a near infrared (NIR) oxygen-sensitive palladium-benzoporphyrin molecule (Pd-MABP) (Montero-Baker et al., 2015; Wisniewski et al., 2017; Chien et al., 2017). The microsensor measures oxygen based on the principle of phosphorescence quenching of metalloporphyrins, a well-established technique with excellent sensitivity and specificity to physiologic oxygen (Wilson et al., 2006; Rumsey et al., 1988; Lo et al., 1996; Vinogradov et al., 2003). The pHEMA hydrogel is biocompatible, has good oxygen permeability, excellent mechanical properties, and a long history of use in medical devices (Montheard et al., 1992). The Pd-MABP molecules are covalently attached to the pHEMA hydrogel, ensuring that the sensing chemistry is retained in the hydrogel structure. The miniature sensors (0.75 x 0.75 x 2.5 mm) are soft and tissue-like to minimize stress at the material-tissue interface.

The non-invasive optical reader was manually positioned over the sensor and a LED pulsed illumination light into the skin above the sensor. A photodetector collected emission light emanating from the sensor. The phosphorescent lifetime, a property of the oxygen sensitive dye, was measured, thereby providing a signal unaffected by optical path permutations such as sensor depth, blood flow fluctuations, oxy/deoxyhemoglobin ratio, melanin content and hydration (Montheard et al., 1992). The current oxygen sensing system measures lifetime to within 2 µs or better, which equates to ~1 µM O₂at physiological oxygen conditions. The temperature was assumed to be at 37°C to convert the lifetime measurement to oxygen concentration.

The oxygen sensor was placed via rectal insertion in the distal large intestine (~1 cm from anus) of mice under isofluorane anesthesia. The phosphorescent signal was collected for two minutes using the optical reader and all values were averaged before conversion to oxygen concentration. Upon reawakening from anesthesia, the sensor was passed naturally from the mouse via peristalsis. Oxygen data were collected on two cohorts of mice with comparable results. Replicate data is shown in Figure 3—figure supplement 3.

Nitrate measurements

Request a detailed protocol

We measured nitrate in frozen mouse fecal samples using the NitraVer X Nitrogen-Nitrate Reaction Set (Hach Company, Loveland CO) following the manufacturer’s instructions. Data below the lower detection limit of 0.00036 mg were not included in plots or statistical tests. These measurements were collected for only a single cohort of mice due to availability of samples.

Fumarate measurements

Request a detailed protocol

We measured fumarate in frozen mouse fecal samples using the Fumarate Colorimetric Assay Kit (Biovision, Milpitas CA) following the manufacturer’s instructions. Data points below the lower detection limit of 200 nmol were not shown on plots, but were included in post-hoc statistical analyses as a pseudo-count of 200 nmol when all measures for a treatment group were below the detection limit. These measurements were collected for only a single cohort of mice due to availability of samples.

DNA isolation from mouse fecal samples

Request a detailed protocol

Metagenomic DNA was isolated from frozen fecal samples using the MoBio (now Qiagen, Hilden Germany) PowerSoil DNA extraction kit following the manufacturer’s instructions (David et al., 2014; 2015).

16S rRNA gene sequencing and processing

Request a detailed protocol

We performed 16S rRNA gene amplicon sequencing using custom barcoded primers targeting the V4 region of the gene (Caporaso et al., 2011) and published protocols (Caporaso et al., 2011; 2012; Maurice et al., 2013). Sequencing was conducted on an Illumina (San Diego, CA) MiSeq with paired end 250 bp reads using the V3 kit. All samples with fewer than 5000 reads were discarded. Sequencing data were processed using QIIME (version 1.8) (Caporaso et al., 2010) to produce an OTU table with 97% cluster similarity. Sequencing was performed on samples from two cohorts of mice with comparable results, but all results presented here are from a single cohort matching the abiotic data presented.

Because baseline microbiome is expected to differ between mouse litters over time as well as between vendor (Rosshart et al., 2017; Rausch et al., 2016; Ivanov et al., 2009; Campbell et al., 2012), we examined whether the effects observed were reproducible between the two cohorts sequenced. An overall analysis of beta-diversity showed that there were significant differences between cohorts in their microbiome at the OTU level throughout the experiment (p<0.01, R² = 0.17, PERMANOVA; Figure 4—figure supplement 5). This cohort effect did not significantly interact with the treatment effect, however (p=0.3). Because the Enterobacteriaceae bloom response to antibiotics has been found to be consistent between humans and mice, as well as between mouse vendors (Antonopoulos et al., 2009; David et al., 2015; Young and Schmidt, 2004; Theriot et al., 2014; Peterfreund et al., 2012; Dethlefsen et al., 2008; Jakobsson et al., 2010; Looft and Allen, 2012), we believe that the overall trends observed here are indicative of what may be occurring in other conventional gut microbiota contexts.

Clustering

Request a detailed protocol

Bacterial genera were clustered using a previously described bioinformatics pipeline for identifying taxa with similar dynamics (David et al., 2014; 2015). Rare genera (defined as those observed in five or fewer samples) were also excluded from clustering. Rare genera comprised 194 of 306 total observed genera, but accounted for only 0.019% of reads. All remaining genera (n = 112) were ultimately assigned to a cluster. A clustering threshold of 0.9 was manually chosen to balance model simplicity (i.e., building clusters with an interpretable number of genera) and fidelity (i.e., capturing the unique dynamics of distinct genera). Only clusters with an average abundance of >1% (five clusters altogether) were included in statistical analyses.

qPCR: Primers and conditions

Request a detailed protocol

To estimate total bacterial abundance, PCR was performed on DNA extracted from feces using the following primers: forward, 5’-ACTCCTACGGGAGGCAGCAGT-3’, reverse, 5’-GTATTACCGCGGCTGCTGGCAC-3’ (Bergström et al., 2012). qPCR assays were run using SYBR FAST qPCR Master Mix (KAPA, Wilmington MA) on a 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City CA). Ct values were standardized against a dilution curve of known concentration and then adjusted for the weight of fecal matter extracted. Load measurements were performed on samples from two cohorts of mice with comparable results, but all data included here is from a single cohort.

Inflammation assays

Request a detailed protocol

Fecal lipocalin 2 (Lcn-2) was used as a noninvasive biomarker for host intestinal inflammation. Lcn-2 was quantified by ELISA following Chassaing et al. (Chassaing et al., 2012). Assays were conducted with mouse lipocalin-2/ngal duo set ELISA kit from R and D Systems (Minneapolis, MN) using a BioTek ELx405rs plate washer and a BioTek Synergy HT plate reader (Bio-Rad, Hercules CA) at the UNC Center for Gastrointestinal Biology and Disease Advanced Analytics Core. These measurements were collected for only a single cohort of mice due to availability of samples.

RNA isolation and RT-PCR

Request a detailed protocol

Total RNA was isolated from fecal pellets stored in RNALater (Thermo Fisher, Waltham, MA) using the MoBio PowerMicrobiome RNA Isolation kit with an added phenol chloroform extraction step. ≤15 ng RNA was reverse transcribed using cDNA Prep Reverse Transcription Master Mix (Fluidigm, South San Francisco, CA) following manufacturer’s instructions. Target transcripts were preamplified for 18 cycles and then diluted 10x. RT-PCR was performed using a BioMark (Fluidigm) on a 48 × 48 chip with Taqman Fast Advanced Master Mix (Thermo Fisher). Three ERCC RNA Spike-in Mix (Thermo Fisher) positive controls and a nontarget negative control of nuclease-free water were also run on the chip. Expression levels were quantified for thirty genes (Supplementary file 3). Of these, six genes could be quantified for at least 20 time points, and so only these genes were included in statistical analyses. Expression levels were normalized to mouse Actb (ΔCt) expression for each time point for each mouse and then compared to the average of control mice for that time point (ΔΔCt). Fold change (2^ΔΔCt) is represented. These measurements were collected for only a single cohort of mice due to availability of samples.

NMR metabolomics

Request a detailed protocol

Samples were prepared by mixing fecal pellets with 200 μl phosphate buffered saline (PBS) prepared in 95% D₂O, pH 7.5 and homogenizing for 5 min at 30 Hz. Cell debris was removed by centrifuging the samples for 10 min at 14K x g at 4°C. The supernatant was filtered through a 0.22 μm centrifugal filter. To maximize the metabolite extraction, each pellet from the first centrifugation was mixed with 200 μl PBS in 95% D₂O, pH 7.5 and homogenized for 5 min at 30 Hz. The solutions were filtered through a 0.22 μm centrifugal filter as before. The supernatants were combined and 300 μl was taken and mixed with 300 μl of 99.9% D₂O containing 500 μM 3-(trimethylsilyl)−1-propanesulfonic acid-d₆, sodium salt (DSS-d₆).

NMR spectra were collected on an Avance III 700MHz NMR spectrometer (Bruker Biospin, Billerica, MA, USA) with a TXI triple resonance probe operating at 25C. A 1D NOESY pulse sequence was used with a spectral width of 12 ppm. Each spectrum was digitized with 32768 points over a 3.9 s total acquisition time. The sequence used a 100 ms mixing time and an inter scan relaxation delay of 2.0 s. These measurements were collected for only a single cohort of mice due to availability of samples.

The spectra were processed using Advanced Chemistry Development Spectrus Processor (version 2016.1, Toronto, ON, Canada). Processing parameters included zero filling to 65526 points, 0.3 Hz decaying exponential multiplication. Phase correction was applied using the automated algorithms in the ACD software and a third order polynomial baseline correction was applied to all samples. Quantitative profiling of the metabolites was carried out with the Chenomx NMR Suite (version 8.2, Edmonton, AL, Canada). The chemical shifts and concentrations were referenced using DSS-d₆ added to each sample at a concentration of 250 μM and set to a chemical shift of 0.00 ppm. The Chenomx 700MHz library included acetate, acetoacetate, acetoin, alanine, asparagine, aspartate, B-hydroxybutyrate, butyrate, DSS-d6, ethanol, formate, glucose, glutamine, glutamate, glycine, isoleucine, isopropanol, lactate, leucine, lysine, methionine, propionate, sarcosine, threonine, and valine. Signals in the bile acid regions of the spectrum, from approximately 0.6 to 0.9 ppm, were fit as synthetic Lorentizan functions and each signal given a generic assignment as a bile acid methyl group. Concentration for bile acids are presented as arbitrary concentration units and relative comparisons should be made. The fitting of each spectrum was carried out in batch mode. Manual adjustment of the fitting was carried out in some samples to correct for errors arising from spectral overlap.

Ex vivo experiments

Culture conditions and treatment

Request a detailed protocol

Gut microbiota were cultured from human stool in a flow-through bioreactor system (Infors, Bottmingen Switzerland) following previous work (McDonald et al., 2013). Healthy human subjects (n = 3) that self-reported no use of antibiotics within a month of enrollment provided a single donation of stool. The sample size was calculated using a power analysis to allow for β less than 0.1. This experiment was only performed with a single cohort of donors. Informed consent was obtained from all subjects and the protocol was approved by the Duke Health Institutional Review Board. Subjects collected samples by placing disposable commode specimen containers (Fisher Scientific, Hampton NH) under their toilet seats before bowel movements. Intact stool samples were moved within roughly 15 min into anaerobic conditions after collection. All bioreactor vessels had a total vessel volume of 400 ml and a continuous turnover of 400 ml per day for the media to emulate the 24 hr average passage time of the human gut. The bioreactors communities were subject to a two-day adjustment period before treatment began (Figure 2—figure supplement 1). Vessels were then randomly assigned a treatment condition. Antibiotic treated cultures were subjected to daily single doses of 3.896 ml antibiotic treatments on two consecutive days. Antibiotics were prepared as above for mice. Control and antibiotic vessels had 1 L per minute N₂ gas bubbled in. Redox was measured continuously with Hamilton (Reno NV) EasyFerm Plus ORP K8 225 probes, which has a measuring range of ±2000 mV. Redox probes were calibrated in Hamilton ORP buffer before autoclaving the bioreactor setup at the beginning of each replicate run. Because of differences in the construction and calibration between the in vivo and ex vivo redox probes, direct comparisons between the in vivo and ex vivo redox values should not be performed. Samples were collected daily and frozen immediately at −80°C for subsequent analysis of microbial community composition. Abiotic data were averaged for five minute increments for analysis.

Statistical analysis

Request a detailed protocol

We performed linear mixed effects model analysis to determine the effects of antibiotics on redox potential, oxygen concentration, fecal lipocalin-2 concentration, and Bray-Curtis dissimilarity. As fixed effects, we entered antibiotic treatment and time with an interaction term into the model. We included mouse identity as a random effect. P values were obtained by likelihood ratio tests comparing the full model against a model including only time and mouse identity and were performed with the ‘anova’ function in the ‘lme4’ package. Repeated measures correlations were used to assess correlations where multiple time points from the same mouse were included in the statistical analysis. Repeated measures correlations were calculated with the ‘rmcorr’ function in the ‘rmcorr’ package. These and all other statistical analyses were carried out in R (R core team, version 3.3). All statistical tests performed were non-parametric except where a Shapiro-Wilks test indicated that data were normally distributed, in which case parametric tests were used. All data points were included in analyses and outliers were not treated in any manner.

Data availability

Sequencing data have been deposited in the European Nucleotide Archive under Primary Accession Number (PRJEB26446).

The following data sets were generated

1. Aspen T Reese
2. Lawrence A David
(2018) Antibiotic induced changes in the microbiota disrupt redox dynamics in the gut
Publicly available at the European Nucleotide Archive (accession no: PRJEB26446).

https://www.ebi.ac.uk/ena/data/view/PRJEB26446

References

1. Albenberg L
2. Esipova TV
3. Judge CP
4. Bittinger K
5. Chen J
6. Laughlin A
7. Grunberg S
8. Baldassano RN
9. Lewis JD
10. Li H
11. Thom SR
12. Bushman FD
13. Vinogradov SA
14. Wu GD
(2014) Correlation between intraluminal oxygen gradient and radial partitioning of intestinal Microbiota
Gastroenterology 147:1055–1063.

https://doi.org/10.1053/j.gastro.2014.07.020
- PubMed
- Google Scholar
1. Allison SD
2. Martiny JB
(2008) Colloquium paper: resistance, resilience, and redundancy in microbial communities
PNAS 105 Suppl 1:11512–11519.

https://doi.org/10.1073/pnas.0801925105
- PubMed
- Google Scholar
(2009) Reproducible community dynamics of the gastrointestinal Microbiota following antibiotic perturbation
Infection and Immunity 77:2367–2375.

https://doi.org/10.1128/IAI.01520-08
- PubMed
- Google Scholar
(1999) Trends in infectious disease mortality in the United States during the 20th century
Jama 281:61–66.

https://doi.org/10.1001/jama.281.1.61
- PubMed
- Google Scholar
1. Arumugam M
2. Raes J
3. Pelletier E
4. Le Paslier D
5. Yamada T
6. Mende DR
7. Fernandes GR
8. Tap J
9. Bruls T
10. Batto JM
11. Bertalan M
12. Borruel N
13. Casellas F
14. Fernandez L
15. Gautier L
16. Hansen T
17. Hattori M
18. Hayashi T
19. Kleerebezem M
20. Kurokawa K
21. Leclerc M
22. Levenez F
23. Manichanh C
24. Nielsen HB
25. Nielsen T
26. Pons N
27. Poulain J
28. Qin J
29. Sicheritz-Ponten T
30. Tims S
31. Torrents D
32. Ugarte E
33. Zoetendal EG
34. Wang J
35. Guarner F
36. Pedersen O
37. de Vos WM
38. Brunak S
39. Doré J
40. Antolín M
41. Artiguenave F
42. Blottiere HM
43. Almeida M
44. Brechot C
45. Cara C
46. Chervaux C
47. Cultrone A
48. Delorme C
49. Denariaz G
50. Dervyn R
51. Foerstner KU
52. Friss C
53. van de Guchte M
54. Guedon E
55. Haimet F
56. Huber W
57. van Hylckama-Vlieg J
58. Jamet A
59. Juste C
60. Kaci G
61. Knol J
62. Lakhdari O
63. Layec S
64. Le Roux K
65. Maguin E
66. Mérieux A
67. Melo Minardi R
68. M'rini C
69. Muller J
70. Oozeer R
71. Parkhill J
72. Renault P
73. Rescigno M
74. Sanchez N
75. Sunagawa S
76. Torrejon A
77. Turner K
78. Vandemeulebrouck G
79. Varela E
80. Winogradsky Y
81. Zeller G
82. Weissenbach J
83. Ehrlich SD
84. Bork P
85. MetaHIT Consortium
(2011) Enterotypes of the human gut microbiome
Nature 473:174–180.

https://doi.org/10.1038/nature09944
- PubMed
- Google Scholar
(2012) Introducing GUt low-density array (GULDA): a validated approach for qPCR-based intestinal microbial community analysis
FEMS Microbiology Letters 337:38–47.

https://doi.org/10.1111/1574-6968.12004
- PubMed
- Google Scholar
1. Broedl UC
2. Schachinger V
3. Lingenhel A
4. Lehrke M
5. Stark R
6. Seibold F
7. Göke B
8. Kronenberg F
9. Parhofer KG
10. Konrad-Zerna A
(2007) Apolipoprotein A-IV is an independent predictor of disease activity in patients with inflammatory bowel disease
Inflammatory Bowel Diseases 13:391–397.

https://doi.org/10.1002/ibd.20078
- PubMed
- Google Scholar
1. Buffie CG
2. Jarchum I
3. Equinda M
4. Lipuma L
5. Gobourne A
6. Viale A
7. Ubeda C
8. Xavier J
9. Pamer EG
(2012) Profound alterations of intestinal microbiota following a single dose of clindamycin results in sustained susceptibility to Clostridium difficile-induced colitis
Infection and Immunity 80:62–73.

https://doi.org/10.1128/IAI.05496-11
- PubMed
- Google Scholar
1. Byndloss MX
2. Olsan EE
3. Rivera-Chávez F
4. Tiffany CR
5. Cevallos SA
6. Lokken KL
7. Torres TP
8. Byndloss AJ
9. Faber F
10. Gao Y
11. Litvak Y
12. Lopez CA
13. Xu G
14. Napoli E
15. Giulivi C
16. Tsolis RM
17. Revzin A
18. Lebrilla CB
19. Bäumler AJ
(2017) Microbiota-activated PPAR-γ signaling inhibits dysbiotic Enterobacteriaceae expansion
Science 357:570–575.

https://doi.org/10.1126/science.aam9949
- PubMed
- Google Scholar
1. Campbell JH
2. Foster CM
3. Vishnivetskaya T
4. Campbell AG
5. Yang ZK
6. Wymore A
7. Palumbo AV
8. Chesler EJ
9. Podar M
(2012) Host genetic and environmental effects on mouse intestinal microbiota
The ISME Journal 6:2033–2044.

https://doi.org/10.1038/ismej.2012.54
- PubMed
- Google Scholar
1. Caporaso JG
2. Kuczynski J
3. Stombaugh J
4. Bittinger K
5. Bushman FD
6. Costello EK
7. Fierer N
8. Peña AG
9. Goodrich JK
10. Gordon JI
11. Huttley GA
12. Kelley ST
13. Knights D
14. Koenig JE
15. Ley RE
16. Lozupone CA
17. McDonald D
18. Muegge BD
19. Pirrung M
20. Reeder J
21. Sevinsky JR
22. Turnbaugh PJ
23. Walters WA
24. Widmann J
25. Yatsunenko T
26. Zaneveld J
27. Knight R
(2010) QIIME allows analysis of high-throughput community sequencing data
Nature Methods 7:335–336.

https://doi.org/10.1038/nmeth.f.303
- PubMed
- Google Scholar
1. Caporaso JG
2. Lauber CL
3. Walters WA
4. Berg-Lyons D
5. Huntley J
6. Fierer N
7. Owens SM
8. Betley J
9. Fraser L
10. Bauer M
11. Gormley N
12. Gilbert JA
13. Smith G
14. Knight R
(2012) Ultra-high-throughput microbial community analysis on the illumina HiSeq and MiSeq platforms
The ISME Journal 6:1621–1624.

https://doi.org/10.1038/ismej.2012.8
- PubMed
- Google Scholar
(2011) Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample
PNAS 108:4516–4522.

https://doi.org/10.1073/pnas.1000080107
- PubMed
- Google Scholar
(2016) NFKB1: a suppressor of inflammation, ageing and cancer
The FEBS Journal 283:1812–1822.

https://doi.org/10.1111/febs.13627
- PubMed
- Google Scholar
(1976) The size, pH, and redox potential of the cecum in mice associated with various microbial floras
Experimental Biology and Medicine 151:260–263.

https://doi.org/10.3181/00379727-151-39187
- Google Scholar
(2012) Fecal lipocalin 2, a sensitive and broadly dynamic non-invasive biomarker for intestinal inflammation
PLoS ONE 7:e44328.

https://doi.org/10.1371/journal.pone.0044328
- PubMed
- Google Scholar
1. Chen J
2. Hanke A
3. Tegetmeyer HE
4. Kattelmann I
5. Sharma R
6. Hamann E
7. Hargesheimer T
8. Kraft B
9. Lenk S
10. Geelhoed JS
11. Hettich RL
12. Strous M
(2017) Impacts of chemical gradients on microbial community structure
The ISME Journal 11:920–931.

https://doi.org/10.1038/ismej.2016.175
- PubMed
- Google Scholar
1. Chien JS
2. Mohammed M
3. Eldik H
4. Ibrahim MM
5. Martinez J
6. Nichols SP
7. Wisniewski N
8. Klitzman B
(2017) Injectable Phosphorescence-based oxygen biosensors identify post ischemic reactive hyperoxia
Scientific Reports 7:8255.

https://doi.org/10.1038/s41598-017-08490-0
- PubMed
- Google Scholar
1. Cho I
2. Yamanishi S
3. Cox L
4. Methé BA
5. Zavadil J
6. Li K
7. Gao Z
8. Mahana D
9. Raju K
10. Teitler I
11. Li H
12. Alekseyenko AV
13. Blaser MJ
(2012) Antibiotics in early life alter the murine colonic microbiome and adiposity
Nature 488:621–626.

https://doi.org/10.1038/nature11400
- PubMed
- Google Scholar
1. Circu ML
2. Aw TY
(2011) Redox biology of the intestine
Free Radical Research 45:1245–1266.

https://doi.org/10.3109/10715762.2011.611509
- PubMed
- Google Scholar
1. Claesson MJ
2. Jeffery IB
3. Conde S
4. Power SE
5. O'Connor EM
6. Cusack S
7. Harris HM
8. Coakley M
9. Lakshminarayanan B
10. O'Sullivan O
11. Fitzgerald GF
12. Deane J
13. O'Connor M
14. Harnedy N
15. O'Connor K
16. O'Mahony D
17. van Sinderen D
18. Wallace M
19. Brennan L
20. Stanton C
21. Marchesi JR
22. Fitzgerald AP
23. Shanahan F
24. Hill C
25. Ross RP
26. O'Toole PW
(2012) Gut microbiota composition correlates with diet and health in the elderly
Nature 488:178–184.

https://doi.org/10.1038/nature11319
- PubMed
- Google Scholar
1. Connell JH
2. Slatyer RO
(1977) Mechanisms of succession in natural communities and their role in community stability and organization
The American Naturalist 111:1119–1144.

https://doi.org/10.1086/283241
- Google Scholar
1. Cowley ES
2. Kopf SH
3. LaRiviere A
4. Ziebis W
5. Newman DK
(2015) Pediatric cystic fibrosis sputum can be chemically dynamic, anoxic, and extremely reduced due to hydrogen sulfide formation
mBio 6:e00767-15.

https://doi.org/10.1128/mBio.00767-15
- PubMed
- Google Scholar
1. David LA
2. Maurice CF
3. Carmody RN
4. Gootenberg DB
5. Button JE
6. Wolfe BE
7. Ling AV
8. Devlin AS
9. Varma Y
10. Fischbach MA
11. Biddinger SB
12. Dutton RJ
13. Turnbaugh PJ
(2014) Diet rapidly and reproducibly alters the human gut microbiome
Nature 505:559–563.

https://doi.org/10.1038/nature12820
- PubMed
- Google Scholar
1. David LA
2. Weil A
3. Ryan ET
4. Calderwood SB
5. Harris JB
6. Chowdhury F
7. Begum Y
8. Qadri F
9. LaRocque RC
10. Turnbaugh PJ
(2015) Gut microbial succession follows acute secretory diarrhea in humans
mBio 6:e00381-15.

https://doi.org/10.1128/mBio.00381-15
- PubMed
- Google Scholar
(2010) Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa
PNAS 107:14691–14696.

https://doi.org/10.1073/pnas.1005963107
- PubMed
- Google Scholar
1. Dedon PC
2. Tannenbaum SR
(2004) Reactive nitrogen species in the chemical biology of inflammation
Archives of Biochemistry and Biophysics 423:12–22.

https://doi.org/10.1016/j.abb.2003.12.017
- PubMed
- Google Scholar
(2008) The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing
PLoS Biology 6:e280.

https://doi.org/10.1371/journal.pbio.0060280
- PubMed
- Google Scholar
1. Dethlefsen L
2. Relman DA
(2011) Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation
PNAS 108:4554–4561.

https://doi.org/10.1073/pnas.1000087107
- PubMed
- Google Scholar
1. Dhall S
2. Do D
3. Garcia M
4. Wijesinghe DS
5. Brandon A
6. Kim J
7. Sanchez A
8. Lyubovitsky J
9. Gallagher S
10. Nothnagel EA
11. Chalfant CE
12. Patel RP
13. Schiller N
14. Martins-Green M
(2014) A novel model of chronic wounds: importance of redox imbalance and biofilm-forming Bacteria for establishment of chronicity
PLoS ONE 9:e109848.

https://doi.org/10.1371/journal.pone.0109848
- PubMed
- Google Scholar
(2016) Gut biogeography of the bacterial microbiota
Nature Reviews Microbiology 14:20–32.

https://doi.org/10.1038/nrmicro3552
- PubMed
- Google Scholar
1. El Aidy S
2. Derrien M
3. Merrifield CA
4. Levenez F
5. Doré J
6. Boekschoten MV
7. Dekker J
8. Holmes E
9. Zoetendal EG
10. van Baarlen P
11. Claus SP
12. Kleerebezem M
(2013) Gut bacteria-host metabolic interplay during conventionalisation of the mouse germfree colon
The ISME Journal 7:743–755.

https://doi.org/10.1038/ismej.2012.142
- PubMed
- Google Scholar
1. Espey MG
(2013) Role of oxygen gradients in shaping redox relationships between the human intestine and its microbiota
Free Radical Biology and Medicine 55:130–140.

https://doi.org/10.1016/j.freeradbiomed.2012.10.554
- PubMed
- Google Scholar
1. Faber F
2. Tran L
3. Byndloss MX
4. Lopez CA
5. Velazquez EM
6. Kerrinnes T
7. Nuccio SP
8. Wangdi T
9. Fiehn O
10. Tsolis RM
11. Bäumler AJ
(2016) Host-mediated sugar oxidation promotes post-antibiotic pathogen expansion
Nature 534:697–699.

https://doi.org/10.1038/nature18597
- PubMed
- Google Scholar
(2010) Changes through time: integrating microorganisms into the study of succession
Research in Microbiology 161:635–642.

https://doi.org/10.1016/j.resmic.2010.06.002
- PubMed
- Google Scholar
1. Fischbach MA
2. Sonnenburg JL
(2011) Eating for two: how metabolism establishes interspecies interactions in the gut
Cell Host & Microbe 10:336–347.

https://doi.org/10.1016/j.chom.2011.10.002
- PubMed
- Google Scholar
1. Fischer J
2. Lindenmayer DB
(2007) Landscape modification and habitat fragmentation: a synthesis
Global Ecology and Biogeography 16:265–280.

https://doi.org/10.1111/j.1466-8238.2007.00287.x
- Google Scholar
1. Friedman N
2. Shriker E
3. Gold B
4. Durman T
5. Zarecki R
6. Ruppin E
7. Mizrahi I
(2017) Diet-induced changes of redox potential underlie compositional shifts in the rumen archaeal community
Environmental Microbiology 19:174–184.

https://doi.org/10.1111/1462-2920.13551
- PubMed
- Google Scholar
(2009) Microbial assemblages as ecosystem engineers of sediment stability
Journal of Soils and Sediments 9:640–652.

https://doi.org/10.1007/s11368-009-0142-5
- Google Scholar
1. Gobbetti M
(1998) The sourdough microflora: interactions of lactic acid bacteria and yeasts
Trends in Food Science & Technology 9:267–274.

https://doi.org/10.1016/S0924-2244(98)00053-3
- Google Scholar
1. Goddard MR
(2008) Quantifying the complexities of Saccharomyces cerevisiae's ecosystem engineering via fermentation
Ecology 89:2077–2082.

https://doi.org/10.1890/07-2060.1
- PubMed
- Google Scholar
1. Goodrich JK
2. Waters JL
3. Poole AC
4. Sutter JL
5. Koren O
6. Blekhman R
7. Beaumont M
8. Van Treuren W
9. Knight R
10. Bell JT
11. Spector TD
12. Clark AG
13. Ley RE
(2014) Human genetics shape the gut microbiome
Cell 159:789–799.

https://doi.org/10.1016/j.cell.2014.09.053
- PubMed
- Google Scholar
(2012) Interactions between the Microbiota and the immune system
Science 336:1268–1273.

https://doi.org/10.1126/science.1223490
- PubMed
- Google Scholar
1. Hughes ER
2. Winter MG
3. Duerkop BA
4. Spiga L
5. Furtado de Carvalho T
6. Zhu W
7. Gillis CC
8. Büttner L
9. Smoot MP
10. Behrendt CL
11. Cherry S
12. Santos RL
13. Hooper LV
14. Winter SE
(2017) Microbial respiration and formate oxidation as metabolic signatures of inflammation-associated dysbiosis
Cell Host & Microbe 21:208–219.

https://doi.org/10.1016/j.chom.2017.01.005
- PubMed
- Google Scholar
1. Ivanov II
2. Atarashi K
3. Manel N
4. Brodie EL
5. Shima T
6. Karaoz U
7. Wei D
8. Goldfarb KC
9. Santee CA
10. Lynch SV
11. Tanoue T
12. Imaoka A
13. Itoh K
14. Takeda K
15. Umesaki Y
16. Honda K
17. Littman DR
(2009) Induction of intestinal Th17 cells by segmented filamentous bacteria
Cell 139:485–498.

https://doi.org/10.1016/j.cell.2009.09.033
- PubMed
- Google Scholar
(2010) Short-term antibiotic treatment has differing long-term impacts on the human throat and gut microbiome
PLoS ONE 5:e9836.

https://doi.org/10.1371/journal.pone.0009836
- PubMed
- Google Scholar
(1994) Organisms as ecosystem engineers
Oikos 69:373–386.

https://doi.org/10.2307/3545850
- Google Scholar
1. Jones SA
2. Chowdhury FZ
3. Fabich AJ
4. Anderson A
5. Schreiner DM
6. House AL
7. Autieri SM
8. Leatham MP
9. Lins JJ
10. Jorgensen M
11. Cohen PS
12. Conway T
(2007) Respiration of Escherichia coli in the mouse intestine
Infection and Immunity 75:4891–4899.

https://doi.org/10.1128/IAI.00484-07
- PubMed
- Google Scholar
1. Jones SA
2. Gibson T
3. Maltby RC
4. Chowdhury FZ
5. Stewart V
6. Cohen PS
7. Conway T
(2011) Anaerobic respiration of Escherichia coli in the mouse intestine
Infection and Immunity 79:4218–4226.

https://doi.org/10.1128/IAI.05395-11
- PubMed
- Google Scholar
(2014) Effects of antibiotics on human microbiota and subsequent disease
Annual Review of Microbiology 68:217–235.

https://doi.org/10.1146/annurev-micro-091313-103456
- PubMed
- Google Scholar
1. Kelly CJ
2. Zheng L
3. Campbell EL
4. Saeedi B
5. Scholz CC
6. Bayless AJ
7. Wilson KE
8. Glover LE
9. Kominsky DJ
10. Magnuson A
11. Weir TL
12. Ehrentraut SF
13. Pickel C
14. Kuhn KA
15. Lanis JM
16. Nguyen V
17. Taylor CT
18. Colgan SP
(2015) Crosstalk between microbiota-derived short-chain fatty acids and intestinal epithelial HIF augments tissue barrier function
Cell Host & Microbe 17:662–671.

https://doi.org/10.1016/j.chom.2015.03.005
- PubMed
- Google Scholar
1. Kröger A
2. Geisler V
3. Lemma E
4. Theis F
5. Lenger R
(1992) Bacterial fumarate respiration
Archives of Microbiology 158:311–314.

https://doi.org/10.1007/BF00245358
- Google Scholar
1. Lo LW
2. Koch CJ
3. Wilson DF
(1996) Calibration of oxygen-dependent quenching of the phosphorescence of Pd-meso-tetra (4-Carboxyphenyl) porphine: a phosphor with general application for measuring oxygen concentration in biological systems
Analytical Biochemistry 236:153–160.

https://doi.org/10.1006/abio.1996.0144
- PubMed
- Google Scholar
1. Looft T
2. Allen HK
(2012) Collateral effects of antibiotics on mammalian gut microbiomes
Gut Microbes 3:463–467.

https://doi.org/10.4161/gmic.21288
- PubMed
- Google Scholar
(2013) Xenobiotics shape the physiology and gene expression of the active human gut microbiome
Cell 152:39–50.

https://doi.org/10.1016/j.cell.2012.10.052
- PubMed
- Google Scholar
(2013) Evaluation of microbial community reproducibility, stability and composition in a human distal gut chemostat model
Journal of Microbiological Methods 95:167–174.

https://doi.org/10.1016/j.mimet.2013.08.008
- PubMed
- Google Scholar
1. Million M
2. Tidjani Alou M
3. Khelaifia S
4. Bachar D
5. Lagier JC
6. Dione N
7. Brah S
8. Hugon P
9. Lombard V
10. Armougom F
11. Fromonot J
12. Robert C
13. Michelle C
14. Diallo A
15. Fabre A
16. Guieu R
17. Sokhna C
18. Henrissat B
19. Parola P
20. Raoult D
(2016) Increased gut redox and depletion of anaerobic and methanogenic prokaryotes in severe acute malnutrition
Scientific Reports 6:26051.

https://doi.org/10.1038/srep26051
- PubMed
- Google Scholar
(2015) The First-in-Man "Si Se Puede" Study for the use of micro-oxygen sensors (MOXYs) to determine dynamic relative oxygen indices in the feet of patients with limb-threatening ischemia during endovascular therapy
Journal of Vascular Surgery 61:1501–1510.

https://doi.org/10.1016/j.jvs.2014.12.060
- PubMed
- Google Scholar
(1992) 2-Hydroxyethyl methacrylate (HEMA): chemical properties and applications in biomedical fields
Journal of Macromolecular Science, Part C: Polymer Reviews 32:1–34.

https://doi.org/10.1080/15321799208018377
- Google Scholar
1. Morris RL
2. Schmidt TM
(2013) Shallow breathing: bacterial life at low O(2)
Nature Reviews Microbiology 11:205–212.

https://doi.org/10.1038/nrmicro2970
- PubMed
- Google Scholar
1. Ng KM
2. Ferreyra JA
3. Higginbottom SK
4. Lynch JB
5. Kashyap PC
6. Gopinath S
7. Naidu N
8. Choudhury B
9. Weimer BC
10. Monack DM
11. Sonnenburg JL
(2013) Microbiota-liberated host sugars facilitate post-antibiotic expansion of enteric pathogens
Nature 502:96–99.

https://doi.org/10.1038/nature12503
- PubMed
- Google Scholar
(2005) Succession of bacterial community structure and diversity in a paddy soil oxygen gradient
Environmental Microbiology 7:382–395.

https://doi.org/10.1111/j.1462-2920.2005.00700.x
- PubMed
- Google Scholar
(2011) Microbial ecology of the dark ocean above, at, and below the seafloor
Microbiology and Molecular Biology Reviews 75:361–422.

https://doi.org/10.1128/MMBR.00039-10
- PubMed
- Google Scholar
1. Osono T
(2005) Colonization and succession of fungi during decomposition of Swida controversa leaf litter
Mycologia 97:589–597.

https://doi.org/10.1080/15572536.2006.11832789
- PubMed
- Google Scholar
1. Pang H
2. Zhang TC
(1998) Fabrication of redox potential microelectrodes for studies in vegetated soils or biofilm systems
Environmental Science & Technology 32:3646–3652.

https://doi.org/10.1021/es980024u
- Google Scholar
1. Peet RK
2. Christensen NL
(1980) Succession: a population process
Vegetatio 43:131–140.

https://doi.org/10.1007/BF00121025
- Google Scholar
1. Pereira FC
2. Berry D
(2017) Microbial nutrient niches in the gut
Environmental Microbiology 19:1366–1378.

https://doi.org/10.1111/1462-2920.13659
- PubMed
- Google Scholar
1. Peterfreund GL
2. Vandivier LE
3. Sinha R
4. Marozsan AJ
5. Olson WC
6. Zhu J
7. Bushman FD
(2012) Succession in the gut microbiome following antibiotic and antibody therapies for Clostridium difficile
PLoS ONE 7:e46966.

https://doi.org/10.1371/journal.pone.0046966
- PubMed
- Google Scholar
(1958) Studies on the ameba-bacteria relationship in Amebiasis. II. some concepts on the etiology of the disease
The American Journal of Tropical Medicine and Hygiene 7:392–399.

https://doi.org/10.4269/ajtmh.1958.7.392
- PubMed
- Google Scholar
1. Rausch P
2. Basic M
3. Batra A
4. Bischoff SC
5. Blaut M
6. Clavel T
7. Gläsner J
8. Gopalakrishnan S
9. Grassl GA
10. Günther C
11. Haller D
12. Hirose M
13. Ibrahim S
14. Loh G
15. Mattner J
16. Nagel S
17. Pabst O
18. Schmidt F
19. Siegmund B
20. Strowig T
21. Volynets V
22. Wirtz S
23. Zeissig S
24. Zeissig Y
25. Bleich A
26. Baines JF
(2016) Analysis of factors contributing to variation in the C57BL/6J fecal microbiota across german animal facilities
International Journal of Medical Microbiology 306:343–355.

https://doi.org/10.1016/j.ijmm.2016.03.004
- PubMed
- Google Scholar
1. Ravcheev DA
2. Thiele I
(2014) Systematic genomic analysis reveals the complementary aerobic and anaerobic respiration capacities of the human gut microbiota
Frontiers in Microbiology 5:674.

https://doi.org/10.3389/fmicb.2014.00674
- PubMed
- Google Scholar
1. Reikvam DH
2. Erofeev A
3. Sandvik A
4. Grcic V
5. Jahnsen FL
6. Gaustad P
7. McCoy KD
8. Macpherson AJ
9. Meza-Zepeda LA
10. Johansen FE
(2011) Depletion of murine intestinal microbiota: effects on gut mucosa and epithelial gene expression
PLoS ONE 6:e17996.

https://doi.org/10.1371/journal.pone.0017996
- PubMed
- Google Scholar
1. Rivera-Chávez F
2. Zhang LF
3. Faber F
4. Lopez CA
5. Byndloss MX
6. Olsan EE
7. Xu G
8. Velazquez EM
9. Lebrilla CB
10. Winter SE
11. Bäumler AJ
(2016) Depletion of Butyrate-Producing clostridia from the gut Microbiota drives an aerobic luminal expansion of salmonella
Cell Host & Microbe 19:443–454.

https://doi.org/10.1016/j.chom.2016.03.004
- PubMed
- Google Scholar
(2017) Wild mouse gut microbiota promotes host fitness and improves disease resistance
Cell 171:1015–1028.

https://doi.org/10.1016/j.cell.2017.09.016
- PubMed
- Google Scholar
1. Rui J
2. Peng J
3. Lu Y
(2009) Succession of bacterial populations during plant residue decomposition in rice field soil
Applied and Environmental Microbiology 75:4879–4886.

https://doi.org/10.1128/AEM.00702-09
- PubMed
- Google Scholar
(1988) Imaging of phosphorescence: a novel method for measuring oxygen distribution in perfused tissue
Science 241:1649–1651.

https://doi.org/10.1126/science.3420417
- PubMed
- Google Scholar
1. Russell SL
2. Gold MJ
3. Hartmann M
4. Willing BP
5. Thorson L
6. Wlodarska M
7. Gill N
8. Blanchet MR
9. Mohn WW
10. McNagny KM
11. Finlay BB
(2012) Early life antibiotic-driven changes in microbiota enhance susceptibility to allergic asthma
EMBO Reports 13:440–447.

https://doi.org/10.1038/embor.2012.32
- PubMed
- Google Scholar
1. Shade A
2. Peter H
3. Allison SD
4. Baho DL
5. Berga M
6. Bürgmann H
7. Huber DH
8. Langenheder S
9. Lennon JT
10. Martiny JB
11. Matulich KL
12. Schmidt TM
13. Handelsman J
(2012) Fundamentals of microbial community resistance and resilience
Frontiers in Microbiology 3:417.

https://doi.org/10.3389/fmicb.2012.00417
- PubMed
- Google Scholar
(2011) H2S: a universal defense against antibiotics in bacteria
Science 334:986–990.

https://doi.org/10.1126/science.1209855
- PubMed
- Google Scholar
1. Shi W
2. Norton JM
(2000) Microbial control of nitrate concentrations in an agricultural soil treated with dairy waste compost or ammonium fertilizer
Soil Biology and Biochemistry 32:1453–1457.

https://doi.org/10.1016/S0038-0717(00)00050-X
- Google Scholar
1. Spees AM
2. Wangdi T
3. Lopez CA
4. Kingsbury DD
5. Xavier MN
6. Winter SE
7. Tsolis RM
8. Bäumler AJ
(2013) Streptomycin-induced inflammation enhances Escherichia coli gut colonization through nitrate respiration
mBio 4:e00430-13.

https://doi.org/10.1128/mBio.00430-13
- PubMed
- Google Scholar
1. Stecher B
2. Robbiani R
3. Walker AW
4. Westendorf AM
5. Barthel M
6. Kremer M
7. Chaffron S
8. Macpherson AJ
9. Buer J
10. Parkhill J
11. Dougan G
12. von Mering C
13. Hardt W-D
(2007) Salmonella enterica serovar typhimurium exploits inflammation to compete with the intestinal Microbiota
PLoS Biology 5:e244.

https://doi.org/10.1371/journal.pbio.0050244
- Google Scholar
1. Stefka AT
2. Feehley T
3. Tripathi P
4. Qiu J
5. McCoy K
6. Mazmanian SK
7. Tjota MY
8. Seo GY
9. Cao S
10. Theriault BR
11. Antonopoulos DA
12. Zhou L
13. Chang EB
14. Fu YX
15. Nagler CR
(2014) Commensal bacteria protect against food allergen sensitization
PNAS 111:13145–13150.

https://doi.org/10.1073/pnas.1412008111
- PubMed
- Google Scholar
1. Theriot CM
2. Koenigsknecht MJ
3. Carlson PE
4. Hatton GE
5. Nelson AM
6. Li B
7. Huffnagle GB
8. Z Li J
9. Young VB
(2014) Antibiotic-induced shifts in the mouse gut microbiome and metabolome increase susceptibility to Clostridium difficile infection
Nature Communications 5:3114.

https://doi.org/10.1038/ncomms4114
- PubMed
- Google Scholar
1. Tilman D
(1985) The resource-ratio hypothesis of plant succession
The American Naturalist 125:827–852.

https://doi.org/10.1086/284382
- Google Scholar
(2016) Variable responses of human microbiomes to dietary supplementation with resistant starch
Microbiome 4:33.

https://doi.org/10.1186/s40168-016-0178-x
- PubMed
- Google Scholar
1. Vinogradov SA
2. Grosul P
3. Rozhkov V
4. Dunphy I
5. Shuman L
6. Dugan BW
7. Evans S
8. Wilson DF
(2003) Oxygen distributions in tissue measured by phosphorescence quenching
Advances in Experimental Medicine and Biology 510:181–185.

https://doi.org/10.1007/978-1-4615-0205-0_30
- PubMed
- Google Scholar
1. Waddell A
2. Ahrens R
3. Tsai YT
4. Sherrill JD
5. Denson LA
6. Steinbrecher KA
7. Hogan SP
(2013) Intestinal CCL11 and eosinophilic inflammation is regulated by myeloid cell-specific RelA/p65 in mice
The Journal of Immunology 190:4773–4785.

https://doi.org/10.4049/jimmunol.1200057
- PubMed
- Google Scholar
1. West PW
2. Ramachandran TP
(1966) Spectrophotometric determination of nitrate using chromotropic acid
Analytica Chimica Acta 35:317–324.

https://doi.org/10.1016/S0003-2670(01)81682-2
- Google Scholar
(2006) Imaging oxygen pressure in the rodent retina by phosphorescence lifetime
Advances in Experimental Medicine and Biology 578:119–124.

https://doi.org/10.1007/0-387-29540-2_19
- PubMed
- Google Scholar
1. Winter SE
2. Thiennimitr P
3. Winter MG
4. Butler BP
5. Huseby DL
6. Crawford RW
7. Russell JM
8. Bevins CL
9. Adams LG
10. Tsolis RM
11. Roth JR
12. Bäumler AJ
(2010) Gut inflammation provides a respiratory electron acceptor for Salmonella
Nature 467:426–429.

https://doi.org/10.1038/nature09415
- PubMed
- Google Scholar
1. Winter SE
2. Winter MG
3. Xavier MN
4. Thiennimitr P
5. Poon V
6. Keestra AM
7. Laughlin RC
8. Gomez G
9. Wu J
10. Lawhon SD
11. Popova IE
12. Parikh SJ
13. Adams LG
14. Tsolis RM
15. Stewart VJ
16. Bäumler AJ
(2013) Host-derived nitrate boosts growth of E. coli in the inflamed gut
Science 339:708–711.

https://doi.org/10.1126/science.1232467
- PubMed
- Google Scholar
1. Winterbourn CC
(2008) Reconciling the chemistry and biology of reactive oxygen species
Nature Chemical Biology 4:278–286.

https://doi.org/10.1038/nchembio.85
- PubMed
- Google Scholar
Book
(2017)
Tissue-Integrating Oxygen Sensors: Continuous Tracking of Tissue Hypoxia

In: Halpern H. J, editors. Advances in Experimental Medicine and Biology. Springer International Publishing. pp. 377–383.
- Google Scholar
1. Wiström J
2. Norrby SR
3. Myhre EB
4. Eriksson S
5. Granström G
6. Lagergren L
7. Englund G
8. Nord CE
9. Svenungsson B
(2001) Frequency of antibiotic-associated diarrhoea in 2462 antibiotic-treated hospitalized patients: a prospective study
Journal of Antimicrobial Chemotherapy 47:43–50.

https://doi.org/10.1093/jac/47.1.43
- PubMed
- Google Scholar
1. Wright JP
2. Jones CG
(2006) The concept of organisms as ecosystem engineers ten years on: progress, limitations, and challenges
BioScience 56:203–209.

https://doi.org/10.1641/0006-3568(2006)056[0203:TCOOAE]2.0.CO;2
- Google Scholar
1. Yatsunenko T
2. Rey FE
3. Manary MJ
4. Trehan I
5. Dominguez-Bello MG
6. Contreras M
7. Magris M
8. Hidalgo G
9. Baldassano RN
10. Anokhin AP
11. Heath AC
12. Warner B
13. Reeder J
14. Kuczynski J
15. Caporaso JG
16. Lozupone CA
17. Lauber C
18. Clemente JC
19. Knights D
20. Knight R
21. Gordon JI
(2012) Human gut microbiome viewed across age and geography
Nature 486:222–227.

https://doi.org/10.1038/nature11053
- PubMed
- Google Scholar
1. Young VB
2. Schmidt TM
(2004) Antibiotic-associated diarrhea accompanied by large-scale alterations in the composition of the fecal microbiota
Journal of Clinical Microbiology 42:1203–1206.

https://doi.org/10.1128/JCM.42.3.1203-1206.2004
- PubMed
- Google Scholar
1. Zhang X
2. Johnston ER
3. Liu W
4. Li L
5. Han X
(2016) Environmental changes affect the assembly of soil bacterial community primarily by mediating stochastic processes
Global Change Biology 22:198–207.

https://doi.org/10.1111/gcb.13080
- PubMed
- Google Scholar
1. Zhu W
2. Winter MG
3. Byndloss MX
4. Spiga L
5. Duerkop BA
6. Hughes ER
7. Büttner L
8. de Lima Romão E
9. Behrendt CL
10. Lopez CA
11. Sifuentes-Dominguez L
12. Huff-Hardy K
13. Wilson RP
14. Gillis CC
15. Tükel Ç
16. Koh AY
17. Burstein E
18. Hooper LV
19. Bäumler AJ
20. Winter SE
(2018) Precision editing of the gut microbiota ameliorates colitis
Nature 553:208–211.

https://doi.org/10.1038/nature25172
- PubMed
- Google Scholar

Article and author information

Author details

Aspen T Reese
1. Department of Biology, Duke University, Durham, United States
2. Department of Molecular Genetics and Microbiology, Duke University, Durham, United States
Contribution
Conceptualization, Data curation, Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9004-9470
Eugenia H Cho

Department of Bioengineering, University of Pennsylvania, Philadelphia, United States

Contribution
Investigation, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Bruce Klitzman

Department of Surgery, Duke University Medical Center, Durham, United States

Contribution
Methodology, Writing—review and editing

Competing interests
Duke University has received support from Profusa, Inc. as an NIH subaward and additional corporate support for research. Dr. Klitzman has no equity ownership. He now has a consulting agreement and has received professorial sponsorship through Duke University from Profusa, Inc.
Scott P Nichols

Profusa, Inc., South San Francisco, United States

Contribution
Software, Investigation, Methodology, Writing—review and editing

Competing interests
Is the Director of Pre-Clinical and Advanced Technologies at Profusa, Inc.
Natalie A Wisniewski

Profusa, Inc., South San Francisco, United States

Contribution
Software, Methodology, Writing—review and editing

Competing interests
Is a Co-Founder and Chief Technology Officer of Profusa, Inc.
Max M Villa

Department of Molecular Genetics and Microbiology, Duke University, Durham, United States

Contribution
Validation, Investigation, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Heather K Durand

Department of Molecular Genetics and Microbiology, Duke University, Durham, United States

Contribution
Investigation, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4059-8484
Sharon Jiang

Department of Molecular Genetics and Microbiology, Duke University, Durham, United States

Contribution
Investigation, Writing—review and editing

Competing interests
No competing interests declared
Firas S Midani

Program in Computational Biology and Bioinformatics, Duke University, Durham, United States

Contribution
Formal analysis, Methodology, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2473-7758
Sai N Nimmagadda

Department of Biomedical Engineering, Duke University, Durham, United States

Contribution
Investigation, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-1837-4054
Thomas M O'Connell

Department of Otolaryngology-Head and Neck Surgery, Indiana University School of Medicine, Indianapolis, United States

Contribution
Data curation, Software, Investigation, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Justin P Wright

Department of Biology, Duke University, Durham, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Writing—review and editing

Competing interests
No competing interests declared
Marc A Deshusses

Department of Civil and Environmental Engineering, Duke University, Durham, United States

Contribution
Conceptualization, Writing—review and editing

Competing interests
No competing interests declared
Lawrence A David
1. Department of Molecular Genetics and Microbiology, Duke University, Durham, United States
2. Program in Computational Biology and Bioinformatics, Duke University, Durham, United States
3. Department of Biomedical Engineering, Duke University, Durham, United States
4. Center for Genomic and Computational Biology, Duke University, Durham, United States
Contribution
Conceptualization, Formal analysis, Supervision, Funding acquisition, Project administration, Writing—review and editing

For correspondence
lawrence.david@duke.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-3570-4767

Funding

National Science Foundation

Aspen T Reese
Justin P Wright
Lawrence A David

Hartwell Foundation

Lawrence A David

Alfred P. Sloan Foundation

Lawrence A David

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

Elise Cowley provided advice on experimental procedures and interpretation of redox dynamics. Lionel Watkins aided in the chemostat experiments. The manuscript was improved by comments from Cari Ficken and Robert Dunn. Inflammation assays were performed by Carlton Anderson at the UNC Center for Gastrointestinal Biology and Disease (CGIBD) Advanced Analytics Core funded by NIH grant P30-DK034987. RT-QPCR was carried out by members of the Duke Sequencing and Genomic Technologies core. This work was significantly improved by the comments of two anonymous reviewers and Jon Clardy.

Ethics

Human subjects: Informed consent was obtained from all subjects and the protocol (Pro00049498) was approved by the Duke Health Institutional Review Board.

Animal experimentation: All animal experiments were conducted in accordance with National Institute of Health Guide for the Care and Use of Laboratory Animals using protocols approved by the Duke University Institutional Animal Care & Use Committee (Protocol A140-14-04).

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.