Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Novel Systems Biology Techniques

Species Deletions from Microbiome Consortia Reveal Key Metabolic Interactions between Gut Microbes

Natalia Gutiérrez, Daniel Garrido
John F. Rawls, Editor
Natalia Gutiérrez
aDepartment of Chemical and Bioprocess Engineering, School of Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel Garrido
aDepartment of Chemical and Bioprocess Engineering, School of Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daniel Garrido
John F. Rawls
Duke University School of Medicine
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00185-19
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

The gut microbiome is a complex microbial community that plays a key role in human health. Diet is an important factor dictating gut microbiome composition. This is mediated by multiple microbe-microbe interactions that result in the fermentation of nondigestible carbohydrates and the production of short-chain fatty acids. Certain species play key metabolic roles in the microbiome, and their disappearance could result in dysbiosis. In this work, a synthetic consortium of 14 gut microbes was studied during the utilization of prebiotic inulin in batch bioreactors. Fermentations were repeated leaving one species out every time, in order to evaluate the impact of their elimination on the system. Substrate consumption, microbial composition, and metabolite production were determined. Single deletions never resulted in a complete loss of bacterial growth or inulin consumption, suggesting functional redundancy. Deletions of Bacteroides dorei and Lachnoclostridium clostridioforme resulted in lower biomass and higher residual inulin. The absence of B. dorei impacted the abundance of the other 10 species negatively. Lachnoclostridium symbiosum, a butyrate producer, appeared to be the most sensitive species to deletions, being stimulated by the presence of Escherichia coli, Bifidobacterium adolescentis, B. dorei, and Lactobacillus plantarum. Conversely, bioreactors without these species did not show butyrate production. L. clostridioforme was observed to be essential for propionate production, and B. dorei for lactate production. Our analysis identified specific members that were essential for the function of the consortium. In conclusion, species deletions from microbial consortia could be a useful approach to identify relevant interactions between microorganisms and defining metabolic roles in the gut microbiome.

IMPORTANCE Gut microbes associate, compete for, and specialize in specific metabolic tasks. These interactions are dictated by the cross-feeding of degradation or fermentation products. However, the individual contribution of microbes to the function of the gut microbiome is difficult to evaluate. It is essential to understand the complexity of microbial interactions and how the presence or absence of specific microorganisms affects the stability and functioning of the gut microbiome. The experimental approach of this study could be used for identifying keystone species, in addition to redundant functions and conditions that contribute to community stability. Redundancy is an important feature of the microbiome, and its reduction could be useful for the design of microbial consortia with desired metabolic properties enhancing the tasks of the keystone species.

INTRODUCTION

The gut microbiome is a complex microbial community characterized by a high cell density, reaching numbers similar to somatic cells (1). Dominant phyla are Firmicutes and Bacteroidetes, which could reach up to 90% of the total bacteria (2). Other phyla with minor but significant representation are Proteobacteria, Actinobacteria, and Verrucomicrobia (2).

The impact of dietary substrates on the gut microbiome is one of the most remarkable and best studied (3, 4). Certain components not absorbed in the small intestine reaching the colon, such as complex carbohydrates and partially degraded proteins, are accessed by gut microbes. Some of these substrates are considered prebiotics and are of great interest in order to modify the gut microbiome for promoting health (5, 6). Fermentation of these macromolecules results in the production of short-chain fatty acids (SCFA), such as acetate, propionate, and butyrate (7). These molecules display several physiological effects on the host (8, 9), such as influencing the stability and physiology of gut environment and serving as an energy source for colonocytes (10). These SCFA represent up to 90% of total acids and have an approximate molar ratio ranging from 75:15:10 to 40:40:20, respectively (11, 12). Among several roles, butyrate helps to maintain the integrity of the mucosa, protecting against cellular inflammation and promoting the removal of dysfunctional cells (13, 14).

The influence of diet on gut microbiome composition and diversity is probably mediated by metabolic interactions between gut microbes (15, 16). These interactions are the basis of complex ecological networks, characterized by cooperative and competitive relationships (17). These could in turn modulate the activity and stability of the gut microbiome. Gut microbes play distinct metabolic roles specializing in the degradation of complex polysaccharides, the fermentation of simple monomers, or the production of essential metabolic intermediates (16, 18). Cross-feeding is a strong force guiding microbiome composition, where certain intermediates produced or released by one microbe are utilized by another. These metabolites derive from the degradation of complex macromolecules, resulting in simple monomers, or from the production of fermentation products such as SCFA or amino acids (19).

The metabolic functions assigned to the gut microbiome are remarkably more conserved than microbial composition, suggesting redundant functionalities (20). However, microbial diversity in the gut microbiome has been shown in several studies to be an important parameter associated with health. In the gut microbiome, certain species have been shown to display a large impact on the structure and function of the community (21). These keystone species could produce unique metabolites connecting two populations or providing essential metabolites for the host. Therefore, they have a crucial role in maintaining the organization of communities due to their functional capabilities and their biotic interactions with other members of the community (22, 23).

Conditions resulting in the loss of keystone species could lead to a dysbiotic state that impairs the integrity of the gut ecosystem (24, 25). Dysbiosis has been linked to certain disorders such as atopic dermatitis and obesity (26) and inflammatory bowel diseases (IBD) (27). Studies have associated IBD with a decrease in the abundance of butyrate-producing species such as Faecalibacterium prausnitzii and Roseburia sp. (10, 28). Interestingly, their reincorporation into the microbiome restores the production of butyrate and consequently promotes the recovery of gut homeostasis (28). This highlights their pivotal importance in the stability and function of the gut microbiome.

To understand the role of the gut microbiome in health and even target the microbiome as a therapeutic target, it is important to define better the metabolic roles performed by gut microbes and their metabolic interactions (29). Recent studies have shown that paired cocultures of gut microbes could explain the behavior of more-complex communities (17, 30). However, it is not well known how individual species could affect the stability and performance of larger microbiome consortia. In this study, we evaluated a species deletion approach to synthetic microbiome consortia, in order to evaluate the impact of single species on the microbial composition and the metabolic function of the microbial community.

RESULTS

General properties of the microbial consortium.The setup of the study is shown in Fig. 1 We included 14 microorganisms frequently found in the gut microbiome (31–33), which belong to the Firmicutes (5 species), Bacteroidetes (7 species), Proteobacteria, and Actinobacteria (Fig. 1). All microorganisms had their genomes sequenced (Table 1). Genome comparison with other strains in the same species indicated that the selected microorganisms were in general representative of their species, with average nucleotide identity values higher than 97% compared to other strains (see data sets at https://figshare.com/s/a3c67977ccf6fe7292eb). Some species had previous evidence for inulin utilization or are endowed with fructofuranosidases (Table 1). The encoded potential for SCFA production was deduced from their genomes (Table 1), indicating that the consortium was able to produce lactate, acetate, propionate, or butyrate. Two species had the potential to produce butyrate: Lachnoclostridium symbiosum and Lachnoclostridium clostridioforme.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Microorganisms in the synthetic consortium and setup of the study. Distance tree indicates genetic distances between 16S rRNA sequences of the microorganisms. The colored circles indicate the samples analyzed for the quantification of variables in the bottom right. Blue and yellow circles indicate OD and inulin consumption determinations, respectively, and green circles indicate microbial composition analysis. The orange circle indicates total SCFA quantification at the end of each experiment.

View this table:
  • View inline
  • View popup
TABLE 1

Description of the strains used in this study and genomic features important for the study

Reproducibility of the culture system.We first evaluated the reproducibility of the culture system, culturing the 14 microorganisms in a batch bioreactor using inulin as carbon source. Microorganisms were inoculated at similar initial concentrations. Culturing the consortium with all species (designated All) showed a coefficient of variation (CV) of 3.3% in optical density (OD) and only 3.4% in the final inulin amount left (see Fig. S1A at https://figshare.com/s/a3c67977ccf6fe7292eb). The relative abundances of each species showed an average CV of 18.6% in the final composition of the replicates (see Fig. S1B at https://figshare.com/s/a3c67977ccf6fe7292eb). The initial composition between the bioreactors showed an average CV of 8.1% (see Fig. S1C at https://figshare.com/s/a3c67977ccf6fe7292eb). These results indicate that the bioreactor is a reproducible culture system for the purpose of evaluating the impact of species deletions.

Bacterial growth.Microbial growth, growth rate, and inulin utilization were determined for each bioreactor after every species deletion (Fig. 2). Every deletion experiment was performed once. The exponential phase started between times 4 and 6 h, and all bioreactors reached stationary phase between 12 and 14 h (Fig. 2B; see also Fig. S2 at https://figshare.com/s/a3c67977ccf6fe7292eb). Escherichia coli and Bacteroides thetaiotaomicron deletions resulted in the highest biomass and growth rates obtained, indicating that these microorganisms contribute negatively to these parameters in the whole consortium (Fig. 2). In contrast, when leaving Bacteroides dorei and Lachnoclostridium clostridioforme out of the consortium, we observed less biomass and a lower growth rate, suggesting they are important for the growth of the consortium. Their deletion could therefore free carrying capacity of the system. Considering all dropouts, the maximum deviation between each bioreactor was an OD of 1 relative to the average (Fig. 2B). The maximum OD variation between deletion experiments was not very high; it went from 4.23 (B. dorei deletion) to 7.36 (E. coli deletion; see Fig. S2 at https://figshare.com/s/a3c67977ccf6fe7292eb). In aggregate, these results indicate that deletion of individual species does not dramatically affect microbial consortia in terms of total growth. However, deletion of certain species resulted in certain alterations increasing or decreasing the biomass of the system.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Growth (maximum OD and growth rate [μ]) and inulin consumption in deletion experiments. (A) Comparative relationship between growth parameters of each bioreactor (see key). Each parameter was normalized to the maximum value observed between experiments (100%). Deleted species are denoted after each triangle: Ec, E. coli; Ba, B. adolescentis; Bd, B. dorei; Bv, B. vulgatus; Bf, B. finegoldii; Bfr, B. fragilis; Bc, B. cellulosilyticus; Bo, B. ovatus; Bt, B. thetaiotaomicron; Lp, L. plantarum; Rg, R. gnavus; Fp, F. plautii; Ls, L. symbiosum; Lc, L. clostridioforme. (B) Average of growth (blue) and inulin consumption (orange) values of all bioreactors. Shaded areas represent the standard deviation range of each time point.

Inulin consumption.Substrate consumption showed variations of up to 40% between bioreactors (Fig. 2B). An almost complete inulin consumption was observed after deleting Bacteroides fragilis (see Fig. S2 at https://figshare.com/s/a3c67977ccf6fe7292eb). Deletion of this microorganism resulted in slower growth compared to all other replicates (Fig. 2A). These results indicate that B. fragilis in this system interferes with inulin consumption. The opposite was observed in L. clostridioforme, B. thetaiotaomicron, and B. dorei deletion bioreactors, where the absence of these microorganisms resulted in an important amount of residual inulin at the end of the fermentation (see Fig. S2 at https://figshare.com/s/a3c67977ccf6fe7292eb). This indicates that either directly or indirectly, they are critical in the consumption of the prebiotic. In summary, this analysis identified certain microorganisms important for inulin utilization in the consortium and also one member that negatively impacted substrate consumption. Interestingly, none of the members of the consortium appeared to be irreplaceable regarding inulin consumption, since all bioreactors displayed fast growth and utilized at least 40% of the inulin. This suggests functional redundancy in the metabolism of the consortium (Table 1).

Microbial composition of bioreactors.The microbial composition of each consortium was determined by quantitative PCR (qPCR) (Fig. 3). Cell copy numbers were in general correlated with OD values (see Fig. S3 at https://figshare.com/s/a3c67977ccf6fe7292eb). In general, the dominant species were Ruminococcus gnavus, E. coli, B. dorei, L. symbiosum, and Lactobacillus plantarum. L. plantarum showed a considerable increase in the absence of B. dorei and L. symbiosum, reaching 50% of the total abundance and suggesting a competitive relationship between these microbes (Fig. 3D and N). Bifidobacterium adolescentis was not particularly dominant in any bioreactor, which was unexpected considering that it is a prominent inulin utilizer (34, 35). One remarkable exception was the E. coli deletion bioreactor, where B. adolescentis increased its relative abundance to more than 60% (Fig. 3B).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Relative bacterial abundances in every bioreactor during growth in inulin at 0, 10, 20, and 30 h. Cell copy numbers were determined by qPCR and were expressed as a proportion of the total number of cell copies in each time point in each bioreactor. Colors at the bottom indicate the relative abundance of each microorganism. Bar plots show the relative microbial composition of the complete consortium (A) or after the deletion of the microorganism indicated in the center at the top (B to O). See Fig. 2A legend for definitions of organism abbreviations.

To quantify the effect of the absence of one microorganism on the abundance of other species in the consortium, an abundance ratio was calculated for each deleted microorganism relative to the abundance of the same strain in the All bioreactor (Fig. 4). Deletions of B. adolescentis and Bacteroides finegoldii from the consortium did not alter the abundance of any other species. While Bacteroides vulgatus and B. finegoldii were mostly unaffected by any dropout, B. adolescentis and Flavonifractor plauti showed reduced abundances in the absence of almost any other member. Finally, deletion of B. dorei from the consortium caused the most deleterious effect on the abundance of 10 microorganisms (Fig. 4), indicating a key role of B. dorei in this consortium.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Effect of the absence of a strain on the abundance of the remaining microorganisms in the consortium. Colored cells indicate the abundance ratio of each bacterial species (top) between a deletion experiment (right) and the abundance of the same microorganism in the All bioreactor. A ratio of 1 is represented in white. Orange boxes indicate that the strain increased its abundance in the bioreactor indicated to the right with respect to the control, and blue boxes indicate a decrease in abundance with respect to the control (All bioreactor). See Fig. 2A legend for definitions of organism abbreviations.

Certain species benefited from the absence of others. The abundance of L. symbiosum was stimulated in the absence of eight species of the consortium but also negatively impacted by the absence of four others (E. coli, B. adolescentis, B. dorei, and L. plantarum). It could be hypothesized that the latter four are important for L. symbiosum growth since their removal impacted L. symbiosum abundance negatively. This also indicates that L. symbiosum is the most sensitive member relative to the presence or absence of other microorganisms in this consortium.

Production of metabolites.The concentrations of lactate, acetate, propionate, and butyrate were determined at the end of each run and normalized to the maximum value observed (Fig. 5; see also Fig. S4 at https://figshare.com/s/a3c67977ccf6fe7292eb). The production of SCFA in the All consortium was 83.9 mM lactate, 21.5 mM acetate, 5.9 mM propionate, and no butyrate.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Production of SCFA in every bioreactor, quantified by HPLC. The concentration of each acid (key at the bottom) was normalized relative to the maximum value of each SCFA determined (100%). SCFA were quantified at 30 h of each experiment from bacterial supernatants. See Fig. 2A legend for definitions of organism abbreviations.

The maximum production of lactate was reached in the B. vulgatus deletion bioreactor. Deletion of this microorganism also allowed the bioreactor to achieve the highest total SCFA concentration (see Fig. S4 at https://figshare.com/s/a3c67977ccf6fe7292eb). In contrast, deletion of B. dorei resulted in the minimal lactate concentration observed, concomitant with the maximum production of acetate observed. This indicates that B. dorei is essential for lactate production, and when it is absent, the consortium switches to high acetate production.

The maximum propionate concentration was observed in the B. finegoldii deletion bioreactor (16.1 mM), and the concomitant loss of six species did not result in the production of propionate (Fig. 5). Three of these were propionate-producing Bacteroides species (B. dorei, B. fragilis, and Bacteroides cellulosilyticus). Three others (E. coli, L. symbiosum, and L. clostridioforme) could contribute indirectly to propionate production.

Finally, the higher concentration of butyrate was observed in R. gnavus and B. vulgatus deletion bioreactors (10.5 mM [Fig. 5]; see also Fig. S4 at https://figshare.com/s/a3c67977ccf6fe7292eb). This suggests that their presence has a negative impact on the production of the acid. In contrast, deletions of L. symbiosum, E. coli, B. adolescentis, B. dorei, and L. plantarum resulted in a lack of butyrate. A linear regression analysis of butyrate production and microbial abundance indicated that the major contribution to the production of this SCFA was attributed to L. symbiosum (Fig. 6). However, the abundance of other species such as F. plautii, L. clostridioforme, and L. plantarum also showed a positive correlation with the quantity of butyrate produced (Fig. 6).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Percent contribution of each species to the butyrate production. The estimated contributions represent a linear relationship between the abundance of the species and butyrate production in each bioreactor. See Fig. 2A legend for definitions of organism abbreviations.

Correlation of variables.To integrate these observations, a multivariate analysis was performed with all the variables described previously. A principal-component analysis showed four distinct clusters between the species deleted in each bioreactor (Fig. 7; see also Fig. S5 at https://figshare.com/s/a3c67977ccf6fe7292eb). These patterns could indicate functional similarities between these species. The cluster containing the All bioreactors included species which did not show much deviation from this control (B. adolescentis, B. finegoldii, L. plantarum, and Bacteroides ovatus), in addition to E. coli. A second cluster contained mostly Bacteroides species (B. vulgatus, B. cellulosilyticus, and B. thetaiotaomicron). Another cluster contained mostly Firmicutes species, which were important for SCFA production (F. plautii and L. clostridioforme) or controlled inulin utilization (B. fragilis). A fourth cluster was obtained with B. dorei and L. symbiosum, where their deletions from the system resulted in comparable results. Similar observations were also obtained when considering all data points (see Fig. S5 at https://figshare.com/s/a3c67977ccf6fe7292eb), where all bioreactors clustered closely at the beginning of the study, and most of the above associations were also found, including additional data points.

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

Principal-component analysis and clustering of consortia. Data analysis included as variables microbial species abundances, SCFA production, inulin consumption, OD values, and growth rate at the final point (t = 30 h). Colors represent the clusters formed by the deletion of each microorganism. See Fig. 2A legend for definitions of organism abbreviations.

While the absence of B. dorei altered several parameters of the community, L. symbiosum abundance was essential for butyrate production (see Fig. S6 at https://figshare.com/s/a3c67977ccf6fe7292eb) and sensitive to the absence of several microorganisms. Finally, the abundance of L. clostridioforme was positively correlated with propionate production, that of L. symbiosum was positively correlated with butyrate generation, and that of L. plantarum was positively correlated with acetate production (see Fig. S6 at https://figshare.com/s/a3c67977ccf6fe7292eb).

We finally built an interaction network derived from a Spearman correlation matrix (Fig. 8). This network shows only statistically significant relationships between the variables (P < 0.05). Several negative correlations between microorganisms and parameters were found (Fig. 8). As expected, the butyrate production was correlated with L. symbiosum. Interestingly, L. symbiosum presence was positively correlated with only two species, E. coli and B. dorei. The presence of these microorganisms also formed a network of positive interactions with L. clostridioforme, B. cellulosilyticus, B. ovatus, L. plantarum, and R. gnavus. The abundances of B. adolescentis, B. fragilis, and B. thetaiotaomicron and acetate production were not significantly correlated with any other variable. After applying a correction of P values by multiple testing, only the interaction between L. symbiosum and butyrate and the interaction between E. coli and B. dorei remained significant.

FIG 8
  • Open in new tab
  • Download powerpoint
FIG 8

Interaction network from the correlation matrix between species abundances, SCFA production, and growth variables (inulin consumption, OD, and μ). These variables were obtained at the end of each run and represented by black dots. Red lines represent negative interactions, and green lines indicate positive interactions. Distance between nodes shows the intensity of the correlation between variables (the closer the dots, the stronger the interaction). See Fig. 2A legend for definitions of organism abbreviations.

Finally, Fig. 9 summarizes some of the findings of this study. A strong positive interaction was observed between B. dorei and L. symbiosum, which appeared essential for the function of the system converting inulin into SCFA, especially butyrate. L. clostridioforme, E. coli, and L. plantarum were species that with a minor contribution appeared supporting the function of B. dorei and L. symbiosum. Certain members of the consortium appeared to be unimportant (B. fragilis, B. thetaiotaomicron, and F. plautii), and a few appeared to contribute negatively to certain parameters (B. finegoldii, B. vulgatus, and R. gnavus).

FIG 9
  • Open in new tab
  • Download powerpoint
FIG 9

Working model of the study summarizing major results. Arrows indicate positive relationships. Dashed lines indicate negative relationships, and lines with points indicate a relationship with positive or negative status unknown. Thick lines indicate that these observations were statistically significant. Orange lines show the relationship between inulin and species that consume it, and blue lines show the relationship between species and the production of SCFA. See Fig. 2A legend for definitions of organism abbreviations.

DISCUSSION

Microbiome complexity and ecological interactions.There is great interest in understanding the forces dictating gut microbiome composition and how the loss of the stability in the gut microbiome contributes to disease. This represents a great challenge, considering that the gut microbiome is composed of hundreds of species and multiple microbe-microbe and host-microbe interactions (36). One unifying principle is that this microbial community is highly responsive to dietary substrates arriving at the colon. This indicates that the chemical nature of dietary substrates being accessed by gut microbes dictates microbial metabolism and microbial interactions, resulting in complex cross-feeding networks. In this study, we analyzed the impact of species deletions on the composition and metabolic function of a synthetic consortium composed of 14 gut microbes during consumption of inulin. Every deletion experiment was performed once, considering that batch bioreactors showed to be reproducible. However, coculture experiments have shown that gut microbe growth could be very sensitive to the initial concentrations of microorganisms (37). Therefore, the lack of replicates is a limitation of this work.

Species deletion is a common approach in ecological theory of food webs (38, 39). Sometimes, the disappearance of one species in a biological community has detrimental effects on the whole community, causing further extinctions of species and impacting ecosystem performance and fitness (40). Interestingly, here we never found a species deletion causing a catastrophic effect in the community: total growth in the bioreactors always reached OD values higher than 4 and inulin utilization was always higher than 50% of the initial concentration. One possible explanation is that medium composition is directly responsible for the failure to observe larger metabolic dependencies between gut microbes. Considering that the substrate was readily consumed, another explanation is functional redundancy in this consortium, where the role of one species was readily replaced by another, leading to comparable final states. Redundancy is a crucial property of the adult gut microbiome, which has been directly associated with microbiome resilience and stability against perturbations (41).

In contrast to these similarities, microbial compositions and metabolic activities were markedly different after deletions of single microorganisms (Fig. 5; see also Fig. S4 at https://figshare.com/s/a3c67977ccf6fe7292eb). These observations are important since it is actually the metabolic activity of the gut microbiome, rather than the composition or taxonomy, which has the most profound effect on host physiology (42). B. dorei was the microbe that when absent caused the most substantial change in the consortium: 10 species showed a decrease in their relative abundance when B. dorei was missing. Concomitantly, total biomass was the lowest of all bioreactors, and a clear switch in SCFA was observed. Therefore, B. dorei could be classified as a keystone species in this consortium (40).

Interestingly, metabolic profiles did not necessarily reflect the metabolism of the dominant microbe in each deletion experiment. The deletion of L. plantarum resulted in low acetate concentrations, also showing a positive correlation between the abundance of L. plantarum and acetate level. This is intriguing since L. plantarum does not produce acetate under these conditions. L. plantarum has a facultative heterofermentative metabolism, producing large amounts of lactate under anaerobic conditions (43). In addition, we did not observe a decrease in lactate concentration in its absence, indicating that other microorganisms contribute to the lactate pool.

In fecal samples, lactate is not typically detected since it is a cross-feeding metabolite that is readily fermented by other gut microbes (44, 45). In feces from subjects with inflammatory bowel diseases, lactate could reach up to 100 mM (46). In this study, the lowest lactate production was observed by deleting B. dorei. Bacteroides species are well known for producing acetate and propionate as end products. However, certain studies indicate that they could also release lactate to the medium, for example in the absence of organic nitrogen sources and especially fermentable carbohydrates such as fructans (47, 48). Therefore, we hypothesize that B. dorei is the major lactate producer in the consortium. Interestingly, B. dorei deletion resulted in the highest acetate concentration obtained and the concomitant dominance of L. plantarum (Fig. 3). This observation could be explained by L. plantarum switching to heterofermentative metabolism or by another dominant microbe in this consortium, such as E. coli, being responsible for the high acetate observed.

Mechanistic evidence for metabolic interactions.Inulin is a well-studied prebiotic that has been generally associated with the stimulation of Lactobacillus and Bifidobacterium (35, 49, 50). The impact of this fructan in the gut microbiome is regarded to be more complex now, stimulating other species such as Anaerostipes and decreasing Bilophila (5). In our consortium, at least five species are well known for inulin utilization (B. adolescentis, B. thetaiotaomicron, L. plantarum, B. dorei, and B. ovatus), and several others are endowed with fructofuranosidases (L. symbiosum, L. clostridioforme, B. cellulosilyticus, B. finegoldii, and R. gnavus) (Table 1) (51). Among these, L. plantarum and B. adolescentis displayed a limited representation in the synthetic consortium, and their growth was usually repressed in the presence of other microbes. Their abundance increased only after prominent members in the consortium were deleted. B. dorei or L. symbiosum deletion increased L. plantarum growth, and E. coli deletion promoted B. adolescentis. Our working model shows that B. dorei is the dominant microbe accessing inulin and releasing smaller fructans and SCFA, which could support E. coli and L. symbiosum growth. How this combination is able to outcompete L. plantarum or B. adolescentis in vitro is intriguing and could be explained by higher growth rates on the substrate or the release of inhibitory molecules.

Our consortium included seven Bacteroides species, including B. dorei, discussed above. Bacteroides species are important members in the microbiome and are endowed with a wide array of polysaccharide utilization machineries. They are mostly generalists able to explore a broad set of nutrients, including proteins (52). Depending on dietary fibers being fed, Bacteroides species could assume distinct metabolic roles, where some species act as primary fermenters accessing complex substrates and releasing smaller degradation products to other species. For example, it was shown that B. ovatus, while growing in inulin, releases smaller fructans which are utilized by other Bacteroides. B. ovatus also supports the growth of B. vulgatus in gnotobiotic mice (53). These two microorganisms showed little representation in this study, which could be explained by a higher inulin-degrading capability of other members of the consortium, such as B. dorei. Certain Bacteroides species are also sensitive to low pH (54), which could limit their growth in our setup (pH 5.5, simulating the conditions of the proximal colon).

Lachnoclostridium symbiosum was a prominent member of this consortium producing butyrate. Although the production of butyrate in L. symbiosum occurs preferably via the 4-aminobutyrate/succinate pathway (55), it can also be synthesized from acetate and lactate (56–58). We hypothesize that acetate is being cross-fed to L. symbiosum, resulting in butyrate production. While this was not demonstrated here, additional evidence is the presence of homologs of an acetate permease in its genome. Another study suggested that L. symbiosum could produce butyrate from lactate (59). Clostridium butyricum is able to use both lactate and acetate for butyrate production (60). L. symbiosum was the most sensitive strain with regard to the absence or presence of other microbes. Deletions of E. coli, B. adolescentis, B. dorei, and L. plantarum resulted in limited growth of L. symbiosum, suggesting they are important for its growth. A recent report studying the pairwise interactions showed that in a 12-species consortium, Faecalibacterium prausnitzii, a keystone species and butyrate producer, was also the most sensitive species, receiving the most positive interactions from other microbes (17). We have also previously reported that L. symbiosum, in the presence of B. dorei, produced high concentrations of butyrate which were shown to reduce inflammation mediated by tumor necrosis factor alpha (TNF-α) in a cellular model (48). These positive interactions, which are dependent on diet, could be further exploited in order to enrich the gut microbiome for butyrate-producing species.

Conclusions.In this study, we analyzed the impact of the species deletions in a synthetic consortium of gut microbes. This approach was useful to identify significant microbial interactions between microorganisms but also specific microbes playing key metabolic roles in the community. We also observed emergent interactions which cannot be detected using paired cocultures. We observed a strong correlation between E. coli, B. dorei, and L. symbiosum, and we defined essential species for butyrate production. In addition to identifying key roles in the gut microbiome, this approach could be useful for the design of microbial consortia with desired metabolic properties.

MATERIALS AND METHODS

Strains and culture media.Strains used in this study are shown in Table 1, and were obtained from BEI Resources, the ATCC, or the UC Davis Culture Collection. For routine experiments, microorganisms were cultured in their respective culture media described in Table 1. Luria-Bertani medium (LB; Becton, Dickinson, Franklin Lakes, NJ) was used directly, while reinforced clostridium medium (RCM; Becton, Dickinson, Franklin Lakes, NJ) was supplemented with 0.5 g/liter of l-cysteine (Loba Chemie, India). de Man, Rogosa, and Sharpe medium (MRS; Becton, Dickinson, Franklin Lakes, NJ) was supplemented with 0.5 g/liter of l-cysteine, except for L. plantarum. Brain heart infusion (BHI; Becton, Dickinson, Franklin Lakes, NJ) was supplemented with 0.5 g/liter of l-cysteine and 0.01 g/liter of hemin. All incubations were performed at 37°C for 24 to 48 h in an anaerobic jar (Anaerocult; Merck, Darmstadt, Germany) with anaerobic packs (GasPak EM; Becton, Dickinson, Franklin Lakes, NJ).

Genomic analysis.Genomes were obtained and analyzed using the Integrated Microbial Gene (IMG) (51). Pairwise average nucleotide identity (ANI) (61) was calculated for genomes used in this study, compared to available finished genomes in the same species. A genetic distance tree was obtained based on the 16S rRNA sequence of each microorganism in IMG. Fructofuranosidase presence in select genomes was estimated using the EC number 3.2.1.26. Inulin utilization and acetate and lactate production were determined from the literature (Table 1) (16). Propionate and butyrate production was confirmed by finding the respective IMG pathways in select genomes. Finally, sporulation was confirmed by the presence of Spo0A homologs and a full set of sporulation genes in select genomes.

Experimental design and bioreactor operation.Sixteen batch experiments were carried out in a 250-ml bioreactor connected to a MyControl system (Mini-bio; Applikon Biotechnology, Netherlands). Two of them were inoculated with the complete consortium (All) run in duplicate. The remaining experiments replicated the same conditions as the control with all strains but leaving one species out; therefore, they were inoculated with 13 strains. The bioreactors were named with Δ followed by an abbreviation of the microorganism left out in each case (Table 1). Microorganisms were cultured using the optimized formulation mZMB (62), with a fixed pH of 5.5. Before being inoculated in each experiment, microorganisms were grown individually in mZMB supplemented with lactose (20 g/liter) as a carbon source under anaerobic conditions for 48 h at 37°C. After 48 h, the OD was measured, and the volume required of each culture was considered to use an initial OD630 of 1 for each microorganism. Cultures were centrifuged at 3,000 rpm for 5 min. The supernatant was discarded, and the pellet was resuspended in the fixed group of mZMB and used to inoculate the reactor. For each bioreactor, mZMB was supplemented with inulin (20 g/liter) (Piping Rock, Ronkonkoma, NY). Tryptone (Becton, Dickinson, Franklin Lakes, NJ) and l-cysteine were used at 34.2 g/liter and 1 g/liter, respectively, and autoclaved directly in the bioreactor in 70 ml of distilled water. The remaining components of mZMB were sterilized using 0.22-μm filters and incorporated into the bioreactor after being autoclaved (inulin and 0.005 g/liter hemin). The bioreactor was inoculated at an initial OD630 of 1 for each microorganism. Foam level was controlled by injecting 200 μl of silicone antifoam (polydimethylsiloxane) with the inoculum. To generate an anaerobic environment inside the bioreactor, nitrogen (99.99% purity grade) was injected at the beginning of the fermentation and remained anaerobic during each experiment. The temperature was set at 37°C, and stirring was set at 90 rpm. pH was maintained constant throughout the fermentation at 5.5, using an automatic injection of 3 M NaOH and 3 M HCl. Samples were taken every 2 h up to 30 h and immediately centrifuged at 10,000 × g for 2 min. Pellets and supernatants were stored at −20°C. Pellets were subsequently used for DNA extraction and relative abundance determination, and supernatants were used for quantification of inulin consumption and production of SCFA.

Substrate consumption.Inulin utilization was quantified using a previously adapted phenol-sulfuric acid method (62, 63). Supernatants from the bioreactors were diluted 1:200 (vol/vol). Inulin standards between 12 and 0.08 μg were prepared. The experiments were performed in 96-well microplates containing 50 μl of a cold-diluted sample (4°C), 150 μl of concentrated sulfuric acid (98% H2SO4), and 30 μl of 5% phenol. Microplates were heated at 90°C for 5 min and cooled on ice for 5 min. Absorbance at 490 nm was determined in a Tecan Infinite M200 Pro plate reader (Tecan Trading AG, Grödig, Austria). Measurements were performed in replicates.

Quantification of SCFA by HPLC.Acetic, lactic, propionic, and butyric acid were quantified at the end of each run in the bioreactors by high-performance liquid chromatography (HPLC) according to a protocol previously reported (64). An ion exchange column of organic acids and carbohydrates, Aminex HPX-87H (Bio-Rad, Hercules, CA), was used in a Lachrom L-700 HPLC system (Hitachi, Japan) at 35°C with a flow rate of 0.45 ml/min of 5 mM H2SO4. As standard, concentrations from 30 g/liter to 0.155 g/liter of each acid were used.

Determination of relative bacterial abundances.Total DNA extraction from cell pellets was performed by an adapted phenol-chloroform-isoamyl protocol (62). DNA was quantified using a NanoQuant plate in a Tecan Infinite M200 Pro microplate reader and diluted to 10 ng/μl. Bacterial relative and absolute abundances were determined by qPCR, using a set of species-specific primers based on unique genes present in each microorganism (see Table S1 at https://figshare.com/s/a3c67977ccf6fe7292eb). qPCRs and concentration measurements were performed as in the work of Medina et al. (65), with annealing temperatures of 62°C. One exception was R. gnavus qPCR, which was performed at 58°C. Reactions were performed in triplicates, and threshold cycle (CT) values were converted into genome copy numbers per milliliter as in the work of Medina et al. (65).

Data analysis.The contribution of microorganisms to the butyrate production was calculated as a linear regression of the abundances (variables) to the production butyrate (response). A single contribution factor was defined by species. The contribution factor was multiplied by the average of the abundance of the respective bacterium determining its contribution. For the multivariate analysis, principal-component analysis and hierarchical clustering were performed with the ClustVis tool (66). An interaction network was built from the Spearman correlation matrix using the ‘corr’ command in Matlab and showing the statistically significant correlations (P value < 0.05). Later, we used the bonf_holm function in Matlab (Bonferroni-Holm correction) for multiple testing.

Data availability.Raw data and calculations of average nucleotide identity between microorganisms in this study and other strains in the same species are available at https://figshare.com/s/a3c67977ccf6fe7292eb.

ACKNOWLEDGMENTS

We thank Veronica Ortuzar, Rodrigo Santibañez, and Pamela Thomson for their support in this study.

This work was funded by FONDEF IDEA ID16i10045.

FOOTNOTES

    • Received March 12, 2019.
    • Accepted June 26, 2019.
  • Copyright © 2019 Gutiérrez and Garrido.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Sender R,
    2. Fuchs S,
    3. Milo R
    . 2016. Revised estimates for the number of human and bacteria cells in the body. PLoS Biol 14:e1002533-14. doi:10.1371/journal.pbio.1002533.
    OpenUrlCrossRef
  2. 2.↵
    1. Zhu L,
    2. Qin S,
    3. Zhai S,
    4. Gao Y,
    5. Li L
    . 2017. Inulin with different degrees of polymerization modulates composition of intestinal microbiota in mice. FEMS Microbiol Lett 364:fnx075. doi:10.1093/femsle/fnx075.
    OpenUrlCrossRef
  3. 3.↵
    1. Kuo S
    . 2013. The interplay between fiber and the intestinal microbiome in the inflammatory response. Adv Nutr 4:16–28. doi:10.3945/an.112.003046.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Lahner E,
    2. Annibale B
    . 2017. Dietary fiber and gut microbiota modulation—a possible switch from disease to health for subjects with diverticular disease? JSM Nutr Disord 1:1001.
    OpenUrl
  5. 5.↵
    1. Makki K,
    2. Deehan EC,
    3. Walter J,
    4. Bäckhed F
    . 2018. The impact of dietary fiber on gut microbiota in host health and disease. Cell Host Microbe 23:705–715. doi:10.1016/j.chom.2018.05.012.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Singh RK,
    2. Chang HW,
    3. Yan D,
    4. Lee KM,
    5. Ucmak D,
    6. Wong K,
    7. Abrouk M,
    8. Farahnik B,
    9. Nakamura M,
    10. Zhu TH,
    11. Bhutani T,
    12. Liao W
    . 2017. Influence of diet on the gut microbiome and implications for human health. J Transl Med 15:1–17. doi:10.1186/s12967-017-1175-y.
    OpenUrlCrossRef
  7. 7.↵
    1. Conlon MA,
    2. Bird AR
    . 2014. The impact of diet and lifestyle on gut microbiota and human health. Nutrients 7:17–44. doi:10.3390/nu7010017.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Sun M,
    2. Wu W,
    3. Liu Z,
    4. Cong Y
    . 2017. Microbiota metabolite short chain fatty acids, GPCR, and inflammatory bowel diseases. J Gastroenterol 52:1–8. doi:10.1007/s00535-016-1242-9.
    OpenUrlCrossRef
  9. 9.↵
    1. Li M,
    2. van Esch B,
    3. Wagenaar GTM,
    4. Garssen J,
    5. Folkerts G,
    6. Henricks P
    . 2018. Pro- and anti-inflammatory effects of short chain fatty acids on immune and endothelial cells. Eur J Pharmacol 831:52–59. doi:10.1016/j.ejphar.2018.05.003.
    OpenUrlCrossRef
  10. 10.↵
    1. Morrison DJ,
    2. Preston T
    . 2016. Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut Microbes 7:189–200. doi:10.1080/19490976.2015.1134082.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Bergman EN
    . 1990. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol Rev 70:567–590. doi:10.1152/physrev.1990.70.2.567.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    1. Macfarlane S,
    2. Macfarlane GT
    . 2003. Regulation of short-chain fatty acid production. Proc Nutr Soc 62:67–72. doi:10.1079/PNS2002207.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Chung WSF,
    2. Meijerink M,
    3. Zeuner B,
    4. Holck J,
    5. Louis P,
    6. Meyer AS,
    7. Wells JM,
    8. Flint HJ,
    9. Duncan SH
    . 2017. Prebiotic potential of pectin and pectic oligosaccharides to promote anti-inflammatory commensal bacteria in the human colon. FEMS Microbiol Ecol 93:fix127. doi:10.1093/femsec/fix127.
    OpenUrlCrossRef
  14. 14.↵
    1. Fung KYC,
    2. Cosgrove L,
    3. Lockett T,
    4. Head R,
    5. Topping DL
    . 2012. A review of the potential mechanisms for the lowering of colorectal oncogenesis by butyrate. Br J Nutr 108:820–831. doi:10.1017/S0007114512001948.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Sonnenburg JL,
    2. Bäckhed F
    . 2016. Diet-microbiota interactions as moderators of human metabolism. Nature 535:56–64. doi:10.1038/nature18846.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Sung J,
    2. Kim S,
    3. Cabatbat JJT,
    4. Jang S,
    5. Jin YS,
    6. Jung GY,
    7. Chia N,
    8. Kim PJ
    . 2017. Global metabolic interaction network of the human gut microbiota for context-specific community-scale analysis. Nat Commun 8:15393. doi:10.1038/ncomms15393.
    OpenUrlCrossRef
  17. 17.↵
    1. Venturelli OS,
    2. Carr AC,
    3. Fisher G,
    4. Hsu RH,
    5. Lau R,
    6. Bowen BP,
    7. Hromada S,
    8. Northen T,
    9. Arkin AP
    . 2018. Deciphering microbial interactions in synthetic human gut microbiome communities. Mol Syst Biol 14:e8157. doi:10.15252/msb.20178157.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Faust K,
    2. Raes J
    . 2012. Microbial interactions: from networks to models. Nat Rev Microbiol 10:538–550. doi:10.1038/nrmicro2832.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Deehan EC,
    2. Duar RM,
    3. Armet AM,
    4. Perez-Muñoz ME,
    5. Jin M,
    6. Walter J
    . 2017. Modulation of the gastrointestinal microbiome with nondigestible fermentable carbohydrates to improve human health. Microbiol Spectr 5(5):BAD-0019-2017. doi:10.1128/microbiolspec.BAD-0019-2017.
    OpenUrlCrossRef
  20. 20.↵
    1. Bauer E,
    2. Thiele I
    . 2018. From network analysis to functional metabolic modeling of the human gut microbiota. mSystems 3:e00209-17. doi:10.1128/mSystems.00209-17.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Gibson TE,
    2. Bashan A,
    3. Cao HT,
    4. Weiss ST,
    5. Liu YY
    . 2016. On the origins and control of community types in the human microbiome. PLoS Comput Biol 12:e1004688-21. doi:10.1371/journal.pcbi.1004688.
    OpenUrlCrossRef
  22. 22.↵
    1. Trosvik P,
    2. de Muinck EJ
    . 2015. Ecology of bacteria in the human gastrointestinal tract—identification of keystone and foundation taxa. Microbiome 3:44. doi:10.1186/s40168-015-0107-4.
    OpenUrlCrossRef
  23. 23.↵
    1. Shetty SA,
    2. Hugenholtz F,
    3. Lahti L,
    4. Smidt H,
    5. de Vos WM
    . 2017. Intestinal microbiome landscaping: insight in community assemblage and implications for microbial modulation strategies. FEMS Microbiol Rev 41:182–199. doi:10.1093/femsre/fuw045.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Jordán F
    . 2009. Keystone species and food webs. Philos Trans R Soc Lond B Biol Sci 364:1733–1741. doi:10.1098/rstb.2008.0335.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Hooks KB,
    2. O’Malley MA
    . 2017. Dysbiosis and its discontents. mBio 8:e01492-17. doi:10.1128/mBio.01492-17.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Davis CD
    . 2016. The gut microbiome and its role in obesity. Nutr Today 51:167–174. doi:10.1097/NT.0000000000000167.
    OpenUrlCrossRef
  27. 27.↵
    1. Nishida A,
    2. Inoue R,
    3. Inatomi O,
    4. Bamba S,
    5. Naito Y,
    6. Andoh A
    . 2018. Gut microbiota in the pathogenesis of inflammatory bowel disease. Clin J Gastroenterol 11:1–10. doi:10.1007/s12328-017-0813-5.
    OpenUrlCrossRef
  28. 28.↵
    1. Geirnaert A,
    2. Calatayud M,
    3. Grootaert C,
    4. Laukens D,
    5. Devriese S,
    6. Smagghe G,
    7. De Vos M,
    8. Boon N,
    9. Van De Wiele T
    . 2017. Butyrate-producing bacteria supplemented in vitro to Crohn’s disease patient microbiota increased butyrate production and enhanced intestinal epithelial barrier integrity. Sci Rep 7:11450. doi:10.1038/s41598-017-11734-8.
    OpenUrlCrossRef
  29. 29.↵
    1. Jordán F,
    2. Lauria M,
    3. Scotti M,
    4. Nguyen T-P,
    5. Praveen P,
    6. Morine M,
    7. Priami C
    . 2015. Diversity of key players in the microbial ecosystems of the human body. Sci Rep 5:15920. doi:10.1038/srep15920.
    OpenUrlCrossRef
  30. 30.↵
    1. Pinto F,
    2. Medina DA,
    3. Pérez-Correa JR,
    4. Garrido D
    . 2017. Modeling metabolic interactions in a consortium of the infant gut microbiome. Front Microbiol 8:2507. doi:10.3389/fmicb.2017.02507.
    OpenUrlCrossRef
  31. 31.↵
    1. Kraal L,
    2. Abubucker S,
    3. Kota K,
    4. Fischbach MA,
    5. Mitreva M
    . 2014. The prevalence of species and strains in the human microbiome: a resource for experimental efforts. PLoS One 9:e97279. doi:10.1371/journal.pone.0097279.
    OpenUrlCrossRef
  32. 32.↵
    1. Yamada T,
    2. Le Paslier D,
    3. Pons N,
    4. Linneberg A,
    5. Guarner F,
    6. Bork P,
    7. Parkhill J,
    8. Qin J,
    9. Zhou Y,
    10. Yu C,
    11. Kristiansen K,
    12. Xu J,
    13. Sicheritz-Ponten T,
    14. Tap J,
    15. Levenez F,
    16. Nielsen T,
    17. Bertalan M,
    18. Li Y,
    19. Zhu H,
    20. Lepage P,
    21. Li S,
    22. Liang H,
    23. Weissenbach J,
    24. Li D,
    25. Li J,
    26. Li S,
    27. Batto J-M,
    28. Brunak S,
    29. Wang J,
    30. Jian M,
    31. Li R,
    32. Nielsen HB,
    33. Manichanh C,
    34. Qin N,
    35. Wang J,
    36. Zhang X,
    37. Cao J,
    38. Zheng H,
    39. Xie Y,
    40. Pedersen O,
    41. Pelletier E,
    42. Turner K,
    43. Yang H,
    44. Renault P,
    45. Raes J,
    46. Li S,
    47. Burgdorf KS,
    48. Wang B,
    49. Doré J,
    50. Ehrlich SD,
    51. Arumugam M,
    52. Hansen T,
    53. Mende DR
    . 2010. A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464:59–65. doi:10.1038/nature08821.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Forster SC,
    2. Kumar N,
    3. Anonye BO,
    4. Almeida A,
    5. Viciani E,
    6. Stares MD,
    7. Dunn M,
    8. Mkandawire TT,
    9. Zhu A,
    10. Shao Y,
    11. Pike LJ,
    12. Louie T,
    13. Browne HP,
    14. Mitchell AL,
    15. Neville BA,
    16. Finn RD,
    17. Lawley TD
    . 2019. A human gut bacterial genome and culture collection for improved metagenomic analyses. Nat Biotechnol 37:186–192. doi:10.1038/s41587-018-0009-7.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Salazar N,
    2. Dewulf EM,
    3. Neyrinck AM,
    4. Bindels LB,
    5. Cani PD,
    6. Mahillon J,
    7. de Vos WM,
    8. Thissen JP,
    9. Gueimonde M,
    10. de los Reyes-Gavilán CG,
    11. Delzenne NM
    . 2015. Inulin-type fructans modulate intestinal Bifidobacterium species populations and decrease fecal short-chain fatty acids in obese women. Clin Nutr 34:501–507. doi:10.1016/j.clnu.2014.06.001.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Vandeputte D,
    2. Falony G,
    3. Vieira-Silva S,
    4. Wang J,
    5. Sailer M,
    6. Theis S,
    7. Verbeke K,
    8. Raes J
    . 2017. Prebiotic inulin-type fructans induce specific changes in the human gut microbiota. Gut 66:1968–1974. doi:10.1136/gutjnl-2016-313271.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Proal A,
    2. Lindseth IA,
    3. Marshall TG
    . 2017. Microbe-microbe and host-microbe interactions drive microbiome dysbiosis and inflammatory processes. Discov Med 23:51–60.
    OpenUrl
  37. 37.↵
    1. D’hoe K,
    2. Vet S,
    3. Faust K,
    4. Moens F,
    5. Falony G,
    6. Gonze D,
    7. Lloréns-Rico V,
    8. Gelens L,
    9. Danckaert J,
    10. De Vuyst L,
    11. Raes J
    . 2018. Integrated culturing, modeling and transcriptomics uncovers complex interactions and emergent behavior in a three-species synthetic gut community. Elife 7:e37090. doi:10.7554/eLife.37090.
    OpenUrlCrossRef
  38. 38.↵
    1. Quince C,
    2. Higgs PG,
    3. McKane AJ
    . 2005. Deleting species from model food webs. Oikos 110:283–296. doi:10.1111/j.0030-1299.2005.13493.x.
    OpenUrlCrossRefWeb of Science
  39. 39.↵
    1. Symstad AJ,
    2. Tilman D,
    3. Willson J,
    4. Knops J
    . 1998. Species loss and ecosystem functioning: effects of species identity and community composition. Oikos 81:389–397. doi:10.2307/3547058.
    OpenUrlCrossRefWeb of Science
  40. 40.↵
    1. Paine RT
    . 1969. A note on trophic complexity and community stability. Am Nat 103:91–93. doi:10.1086/282586.
    OpenUrlCrossRefWeb of Science
  41. 41.↵
    1. Lozupone CA,
    2. Stombaugh JI,
    3. Gordon JI,
    4. Jansson JK,
    5. Knight R
    . 2012. Diversity, stability and resilience of the human gut microbiota. Nature 489:220–230. doi:10.1038/nature11550.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Krishnan S,
    2. Alden N,
    3. Lee K
    . 2015. Pathways and functions of gut microbiota metabolism impacting host physiology. Curr Opin Biotechnol 36:137–145. doi:10.1016/j.copbio.2015.08.015.
    OpenUrlCrossRef
  43. 43.↵
    1. Quatravaux S,
    2. Remize F,
    3. Bryckaert E,
    4. Colavizza D,
    5. Guzzo J
    . 2006. Examination of Lactobacillus plantarum lactate metabolism side effects in relation to the modulation of aeration parameters. J Appl Microbiol 101:903–912. doi:10.1111/j.1365-2672.2006.02955.x.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Bourriaud C,
    2. Robins RJ,
    3. Martin L,
    4. Kozlowski F,
    5. Tenailleau E,
    6. Cherbut C,
    7. Michel C
    . 2005. Lactate is mainly fermented to butyrate by human intestinal microfloras but inter-individual variation is evident. J Appl Microbiol 99:201–212. doi:10.1111/j.1365-2672.2005.02605.x.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Takahashi K,
    2. Terashima H,
    3. Kohno K,
    4. Ohkohchi N
    . 2013. A stand-alone synbiotic treatment for the prevention of D-lactic acidosis in short bowel syndrome. Int Surg 98:110–113. doi:10.9738/CC169.
    OpenUrlCrossRef
  46. 46.↵
    1. Mayeur C,
    2. Gratadoux JJ,
    3. Bridonneau C,
    4. Chegdani F,
    5. Larroque B,
    6. Kapel N,
    7. Corcos O,
    8. Thomas M,
    9. Joly F
    . 2013. Faecal D/L lactate ratio is a metabolic signature of microbiota imbalance in patients with short bowel syndrome. PLoS One 8:e54335. doi:10.1371/journal.pone.0054335.
    OpenUrlCrossRef
  47. 47.↵
    1. Ríos-Covián D,
    2. Ruas-Madiedo P,
    3. Margolles A,
    4. Gueimonde M,
    5. de Los Reyes-Gavilán CG,
    6. Salazar N
    . 2016. Intestinal short chain fatty acids and their link with diet and human health. Front Microbiol 7:185. doi:10.3389/fmicb.2016.00185.
    OpenUrlCrossRef
  48. 48.↵
    1. Thomson P,
    2. Medina DA,
    3. Ortúzar V,
    4. Gotteland M,
    5. Garrido D
    . 2018. Anti-inflammatory effect of microbial consortia during the utilization of dietary polysaccharides. Food Res Int 109:14–23. doi:10.1016/j.foodres.2018.04.008.
    OpenUrlCrossRef
  49. 49.↵
    1. Rossi M,
    2. Corradini C,
    3. Amaretti A,
    4. Nicolini M,
    5. Pompei A,
    6. Zanoni S,
    7. Matteuzzi D
    . 2005. Fermentation of fructooligosaccharides and inulin by bifidobacteria: a comparative study of pure and fecal cultures. Appl Environ Microbiol 71:6150–6158. doi:10.1128/AEM.71.10.6150-6158.2005.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. Takemura N,
    2. Ozawa K,
    3. Kimura N,
    4. Watanabe J,
    5. Sonoyama K
    . 2010. Inulin-type fructans stimulated the growth of exogenously administered Lactobacillus plantarum No. 14 in the mouse gastrointestinal tract. Biosci Biotechnol Biochem 74:375–381. doi:10.1271/bbb.90794.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Chen IMA,
    2. Markowitz VM,
    3. Chu K,
    4. Palaniappan K,
    5. Szeto E,
    6. Pillay M,
    7. Ratner A,
    8. Huang J,
    9. Andersen E,
    10. Huntemann M,
    11. Varghese N,
    12. Hadjithomas M,
    13. Tennessen K,
    14. Nielsen T,
    15. Ivanova NN,
    16. Kyrpides NC
    . 2017. IMG/M: integrated genome and metagenome comparative data analysis system. Nucleic Acids Res 45:D507–D516. doi:10.1093/nar/gkw929.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Vieira-Silva S,
    2. Falony G,
    3. Darzi Y,
    4. Lima-Mendez G,
    5. Garcia Yunta R,
    6. Okuda S,
    7. Vandeputte D,
    8. Valles-Colomer M,
    9. Hildebrand F,
    10. Chaffron S,
    11. Raes J
    . 2016. Species-function relationships shape ecological properties of the human gut microbiome. Nat Microbiol 1:16088. doi:10.1038/nmicrobiol.2016.88.
    OpenUrlCrossRef
  53. 53.↵
    1. Rakoff-Nahoum S,
    2. Coyne MJ,
    3. Comstock LE
    . 2014. An ecological network of polysaccharide utilization among human intestinal symbionts. Curr Biol 24:40–49. doi:10.1016/j.cub.2013.10.077.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Walker AW,
    2. Duncan SH,
    3. McWilliam Leitch EC,
    4. Child MW,
    5. Flint HJ
    . 2005. pH and peptide supply can radically alter bacterial populations and short-chain fatty acid ratios within microbial communities from the human colon. Appl Environ Microbiol 71:3692–3700. doi:10.1128/AEM.71.7.3692-3700.2005.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Vital M,
    2. Howe A,
    3. Tiedje J
    . 2014. Revealing the bacterial synthesis pathways by analyzing (meta) genomic data. mBio 5:e00889-14. doi:10.1128/mBio.00889-14.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. Van Den Abbeele P,
    2. Belzer C,
    3. Goossens M,
    4. Kleerebezem M,
    5. De Vos WM,
    6. Thas O,
    7. De Weirdt R,
    8. Kerckhof FM,
    9. Van De Wiele T
    . 2013. Butyrate-producing Clostridium cluster XIVa species specifically colonize mucins in an in vitro gut model. ISME J 7:949–961. doi:10.1038/ismej.2012.158.
    OpenUrlCrossRef
  57. 57.↵
    1. Duncan SH,
    2. Louis P,
    3. Flint HJ
    . 2004. Lactate-utilizing bacteria, isolated from human feces, that produce butyrate as a major fermentation product. Appl Environ Microbiol 70:5810–5817. doi:10.1128/AEM.70.10.5810-5817.2004.
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    1. Morrison DJ,
    2. Mackay WG,
    3. Edwards CA,
    4. Preston T,
    5. Dodson B,
    6. Weaver LT
    . 2006. Butyrate production from oligofructose fermentation by the human faecal flora: what is the contribution of extracellular acetate and lactate? Br J Nutr 96:570–577.
    OpenUrlPubMed
  59. 59.↵
    1. Lozupone C,
    2. Faust K,
    3. Raes J,
    4. Faith JJ,
    5. Frank DN,
    6. Zaneveld J,
    7. Gordon JI,
    8. Knight R
    . 2012. Identifying genomic and metabolic features that can underlie early successional and opportunistic lifestyles of human gut symbionts. Genome Res 22:1974–1984. doi:10.1101/gr.138198.112.
    OpenUrlAbstract/FREE Full Text
  60. 60.↵
    1. Detman A,
    2. Mielecki D,
    3. Chojnacka A,
    4. Salamon A,
    5. Błaszczyk MK,
    6. Sikora A
    . 2019. Cell factories converting lactate and acetate to butyrate: Clostridium butyricum and microbial communities from dark fermentation bioreactors. Microb Cell Fact 18:36. doi:10.1186/s12934-019-1085-1.
    OpenUrlCrossRef
  61. 61.↵
    1. Richter M,
    2. Rosselló-Móra R
    . 2009. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A 106:19126–19131. doi:10.1073/pnas.0906412106.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Medina DA,
    2. Pinto F,
    3. Ovalle A,
    4. Thomson P,
    5. Garrido D
    . 2017. Prebiotics mediate microbial interactions in a consortium of the infant gut microbiome. Int J Mol Sci 18:E2095. doi:10.3390/ijms18102095.
    OpenUrlCrossRef
  63. 63.↵
    1. Tuomivaara ST,
    2. Yaoi K,
    3. O’Neill MA,
    4. York WS
    . 2015. Generation and structural validation of a library of diverse xyloglucan-derived oligosaccharides, including an update on xyloglucan nomenclature. Carbohydr Res 402:56–66. doi:10.1016/j.carres.2014.06.031.
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    1. Mendoza SN,
    2. Cañón PM,
    3. Contreras Á,
    4. Ribbeck M,
    5. Agosín E
    . 2017. Genome-scale reconstruction of the metabolic network in Oenococcus oeni to assess wine malolactic fermentation. Front Microbiol 8:534. doi:10.3389/fmicb.2017.00534.
    OpenUrlCrossRef
  65. 65.↵
    1. Medina DA,
    2. Pedreros JP,
    3. Turiel D,
    4. Quezada N,
    5. Pimentel F,
    6. Escalona A,
    7. Garrido D
    . 2017. Distinct patterns in the gut microbiota after surgical or medical therapy in obese patients. PeerJ 2017:e3443. doi:10.7717/peerj.3443.
    OpenUrlCrossRef
  66. 66.↵
    1. Metsalu T,
    2. Vilo J
    . 2015. ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res 43:W566–W570. doi:10.1093/nar/gkv468.
    OpenUrlCrossRefPubMed
  67. 67.
    1. Reichardt N,
    2. Duncan SH,
    3. Young P,
    4. Belenguer A,
    5. McWilliam Leitch C,
    6. Scott KP,
    7. Flint HJ,
    8. Louis P
    . 2014. Phylogenetic distribution of three pathways for propionate production within the human gut microbiota. ISME J 8:1323–1335. doi:10.1038/ismej.2014.14.
    OpenUrlCrossRefPubMed
  68. 68.
    1. Carlier J-P,
    2. Bedora-Faure M,
    3. K’ouas G,
    4. Alauzet C,
    5. Mory F
    . 2010. Proposal to unify Clostridium orbiscindens Winter et al. 1991 and Eubacterium plautii (Seguin 1928) Hofstad and Aasjord 1982, with description of Flavonifractor plautii gen. nov., comb. nov., and reassignment of Bacteroides capillosus to Pseudoflavonifractor capillosus gen. nov., comb. nov. Int J Syst Evol Microbiol 60:585–590. doi:10.1099/ijs.0.016725-0.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.
    1. Alexander CJ,
    2. Citron DM,
    3. Brazier JS,
    4. Goldstein E
    . 1995. Identification and antimicrobial resistance patterns of clinical isolates of Clostridium clostridioforme, Clostridium innocuum, and Clostridium ramosum compared with those of clinical isolates of Clostridium perfringens. J Clin Microbiol 33:3209–3215.
    OpenUrlAbstract/FREE Full Text
  70. 70.
    1. Rios-Covian D,
    2. Salazar N,
    3. Gueimonde M,
    4. de Los Reyes-Gavilan CG
    . 2017. Shaping the metabolism of intestinal Bacteroides population through diet to improve human health. Front Microbiol 8:376. doi:10.3389/fmicb.2017.00376.
    OpenUrlCrossRef
  71. 71.
    1. Robert C,
    2. Chassard C,
    3. Lawson PA,
    4. Bernalier-Donadille A
    . 2007. Bacteroides cellulosilyticus sp. nov., a cellulolytic bacterium from the human gut microbial community. Int J Syst Evol Microbiol 57:1516–1520. doi:10.1099/ijs.0.64998-0.
    OpenUrlCrossRefPubMed
  72. 72.
    1. Adamberg S,
    2. Tomson K,
    3. Vija H,
    4. Puurand M,
    5. Kabanova N,
    6. Visnapuu T,
    7. Jaμgi E,
    8. Alamae T,
    9. Adamberg K
    . 2014. Degradation of fructans and production of propionic acid by Bacteroides thetaiotaomicron are enhanced by the shortage of amino acids. Front Nutr 1:21. doi:10.3389/fnut.2014.00021.
    OpenUrlCrossRefPubMed
  73. 73.
    1. Liu H,
    2. Kang J,
    3. Qi Q,
    4. Chen G
    . 2011. Production of lactate in Escherichia coli by redox regulation genetically and physiologically. Appl Biochem Biotechnol 164:162–169. doi:10.1007/s12010-010-9123-9.
    OpenUrlCrossRefPubMed
  74. 74.
    1. Ramirez-Farias C,
    2. Slezak K,
    3. Fuller Z,
    4. Duncan A,
    5. Holtrop G,
    6. Louis P
    . 2009. Effect of inulin on the human gut microbiota: stimulation of Bifidobacterium adolescentis and Faecalibacterium prausnitzii. Br J Nutr 101:541–550. doi:10.1017/S0007114508019880.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Species Deletions from Microbiome Consortia Reveal Key Metabolic Interactions between Gut Microbes
Natalia Gutiérrez, Daniel Garrido
mSystems Jul 2019, 4 (4) e00185-19; DOI: 10.1128/mSystems.00185-19

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Species Deletions from Microbiome Consortia Reveal Key Metabolic Interactions between Gut Microbes
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Species Deletions from Microbiome Consortia Reveal Key Metabolic Interactions between Gut Microbes
Natalia Gutiérrez, Daniel Garrido
mSystems Jul 2019, 4 (4) e00185-19; DOI: 10.1128/mSystems.00185-19
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

bioreactor
keystone species
microbiome
butyrate
metabolic interaction

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077