Model-based quantification of metabolic interactions from dynamic microbial-community data

Mark Hanemaaijer; Brett G. Olivier; Wilfred F. M. Röling; Frank J. Bruggeman; Bas Teusink

doi:10.1371/journal.pone.0173183

Abstract

An important challenge in microbial ecology is to infer metabolic-exchange fluxes between growing microbial species from community-level data, concerning species abundances and metabolite concentrations. Here we apply a model-based approach to integrate such experimental data and thereby infer metabolic-exchange fluxes. We designed a synthetic anaerobic co-culture of Clostridium acetobutylicum and Wolinella succinogenes that interact via interspecies hydrogen transfer and applied different environmental conditions for which we expected the metabolic-exchange rates to change. We used stoichiometric models of the metabolism of the two microorganisms that represents our current physiological understanding and found that this understanding - the model - is sufficient to infer the identity and magnitude of the metabolic-exchange fluxes and it suggested unexpected interactions. Where the model could not fit all experimental data, it indicates specific requirement for further physiological studies. We show that the nitrogen source influences the rate of interspecies hydrogen transfer in the co-culture. Additionally, the model can predict the intracellular fluxes and optimal metabolic exchange rates, which can point to engineering strategies. This study therefore offers a realistic illustration of the strengths and weaknesses of model-based integration of heterogenous data that makes inference of metabolic-exchange fluxes possible from community-level experimental data.

Citation: Hanemaaijer M, Olivier BG, Röling WFM, Bruggeman FJ, Teusink B (2017) Model-based quantification of metabolic interactions from dynamic microbial-community data. PLoS ONE 12(3): e0173183. https://doi.org/10.1371/journal.pone.0173183

Editor: Stephen S. Fong, Virginia Commonwealth University, UNITED STATES

Received: September 6, 2016; Accepted: February 16, 2017; Published: March 9, 2017

Copyright: © 2017 Hanemaaijer et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files and on https://sourceforge.net/projects/cbmpy/files/publications/data/2017_Hanemaaijer/.

Funding: This work is financially supported by ZonMW and the Netherlands Genomics Initiative via a Zenith Horizon grant 40-41009-98-10038. BGO is supported by a BE-Basic Foundation grant F08.005.001.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Microbial communities carry out important processes for the planet’s ecosystem, animal health and industrial purposes. Engineering such communities is not straightforward; they are generally composed of many interacting microorganisms, and they are highly dynamic. Developing methods that relate the systemic properties of communities to the underlying metabolic processes and interactions of the community members is a major challenge in microbial ecology. Those interactions drive community behaviour, including its (in-)stability upon environmental perturbations [1]. Metabolic interactions between community members are widespread and generally considered to be the dominating interactions [2–4].

Current approaches focus mostly on correlation and co-occurrence of species members for inference of community-interaction partners [5, 6]. These approaches importantly predict the community-interaction structure, even for large systems. It does not, however, inform us about interaction mechanisms and their importance for species survival and community properties. This can be achieved when we know the metabolic exchange fluxes between microorganisms; then we would capture both mechanism and importance of those fluxes, as those can be linked to intracellular metabolic activities and growth. With this additional information we could rationally design approaches that alter community behaviour upon external perturbations.

Unfortunately, quantifying metabolic exchange fluxes—metabolic interactions—directly from experimental data is rarely possible. Generally, only net fluxes are inferred from dynamic metabolite levels measured at the community level that result from contributions of many species. To determine the individual contributions of those species, we suggested that their metabolic capacities, expressed in terms of quantitative models, should be integrated with experimental data [7]. Here we illustrate that this indeed allows for the identification and quantification of metabolic interactions between microorganisms and that those interactions are dependent on environmental conditions.

Our approach relies on stoichiometric models of metabolism and the linkage of the metabolism of microbial species in the community. Stoichiometric modeling of the metabolism of single microorganisms have been developed in systems biology in the last two decades. Recently, such models are being considered for microbial communities [8–11], but the number of studies that combine such metabolic models with experimental data is still limited. The first study was performed by Stolyar et al. on a methanogenic co-culture [12] and several other co-culture studies followed [13, 14]. Also, purely computational studies investigated the potential interactions in a community [15], designed medium compositions that enforces metabolic interactions [16] or calculated biomass ratios and fluxes under balanced growth conditions of microbial communities [17]. The focus of those studies was mostly the prediction of the community phenotype. However, another application of those models is to use them to infer metabolic exchange fluxes between microorganisms from experimental data [7].

Here we demonstrate that we are able to quantify community interactions by combining experimental data, consisting of biomass abundances and metabolite concentrations, with stoichiometric metabolism models. For this purpose, we designed a synthetic co-culture that serves as a model for anaerobic interspecies hydrogen transfer, which is one of the driving forces in the anaerobic digestion process, and allows for tunable metabolic interactions by changing the environmental conditions. Quantification of the metabolic interactions contributes to identifying how this important process can be improved. Although synthetic communities are not as complex as natural ecosystems, they have shown to be useful for microbial ecologists to examine ecological theories, such as the influence of community evenness on the functionality of a community [18]. Additionally, synthetic communities are used to improve industrial applications, such as bioremediation [19] or bioethanol production [20]. The co-culture we designed consists of the anaerobic bacteria Clostridium acetobutylicum and Wolinella succinogenes. C. acetobutylicum produces H₂ that subsequently is consumed by W. succinogenes to reduce nitrate (NO₃^-) into nitrite (NO₂^-) or ammonium (NH₄⁺). NO₂^- or NH₄⁺ can act as a nitrogen-source for C. acetobutylicum. We varied the nitrogen source to vary the metabolic interactions in the community, which we subsequently inferred, using a model and experimental data. This study illustrates how the integration of a model, representing our knowledge about the metabolism of the two species, with experimental data leads to the quantification of metabolic exchange fluxes.

Materials and methods

Strains and cultivation conditions

Co-culture experiments of Clostridium acetobutylicum DSM-792 with Wolinella succinogenes DSM-1740 were grown in Widdel-based medium containing the following components (per liter): 1 g NaCl, 0.4 g MgCl⋅6H₂O, 0.1 g CaCl⋅2H₂O, 0.5 g KCl, 0.59 g glucose⋅H₂O, 0.15 g cysteine, supplemented with 1 ml of 100x RPMI-1640 vitamin solution (Sigma-Aldrich) and 1 ml of trace elements solution, containing (per liter): 0.5 g EDTA, 3 g MgSO₄⋅7H₂O, 0.5 g MnSO₄⋅H₂O, 1 g NaCl, 0.1 g FeSO₄⋅7H₂O, 0.1 g Co(NO₃)₂⋅6H₂O, 0.1 g CaCl₂, 0.1 g ZnSO₄⋅7H₂O, 10 mg CuSO₄⋅5H₂O, 10 mg AIK(SO₄)₂, 10 mg H₃BO₂, 10 mg Na₂MoO₄⋅2H₂O, 1 mg Na₂SeO₃, 10 mg Na₂WO₄⋅2H₂O, 20 mg NiCl₂⋅6H₂O. A 20 mM Na-phosphate buffer pH 7.0 was used to maintain a constant pH during the growth experiment. The final pH of the medium was 7.0. The medium was kept anaerobic by flushing with a gas mixture of 10% CO₂ and 90% N₂. Resazurine (0.01 mg/l) was used as an indicator for anaerobic conditions. The growth experiment was performed at 37°C in 100 ml septum vials with rubber stoppers containing 50 ml of growth medium. ‘NH₄⁺’ condition contained 0.25 g/L NH₄Cl, ‘NO₃^-’ condition contained 0.85 g/l NaNO₃^-, ‘NH₄⁺ + NO₃^-’ condition contained 0.25 g/l NH₄Cl and 0.85 g/l NaNO₃^- and ‘N₂’ condition contained no N-source, except N₂ in the headspace. Cells were taken from actively growing stock, cultivated in monoculture first in abovementioned medium, but with yeast extract supplemented. The monocultures were added 1:1 (OD/OD) in the media, pre-grown for 70 generations via serial transfer of the co-culture in NH₄⁺ + NO₃^- containing medium to create a stable co-culture. We are aware that mutations could arise during this period which potentially influence the community phenotype, but that is beyond the scope of this study. The co-culture was 1:100 propagated to fresh medium once it was fully grown. The experiment started when the co-culture reached exponential growth and was 1:100 propagated to fresh medium containing different nitrogen sources in triplicate.

Metabolite analysis

Supernatant from the batch cultivation was taken from the septum vials by taking 1 ml of growth medium and filtered through a 0.2 μm polyethersulfone (PES) filter and stored at -20°C until further processing. Samples were analyzed for fermentation products (i.e. formate, succinate, acetate, propionate, lactate, and butyrate) by high performance liquid chromatography (HPLC) with a Phenomenex Rezex ROA 300x7.8mm column at 55°C and occupied with RID 10A and SPD 20A detectors having a flow rate of 0.5ml/min and eluent of 5mM H₂SO₄. NO₃^- was measured according Yang et al. [21] and NO₂^- was measured spectrophotometrically at 520nm by diazotization of sulphanilamide by nitrite. CO₂ and H₂ were measured by taking 100 μl gas sample from the head space with a syringe (Terumo 1 mL syringe) and the gas was immediately analyzed on a Shimadzu GC2010 system equipped with a Carboxen 1010 column (30m⋅0.58mm, 105°C) and a barrier discharge ionization detector (BID). Helium was used as a carrier gas at 28.2 mL/min at a pressure of 37.2 kPa.

Quantification of cells

DNA was extracted with the PowerSoil DNA isolation kit (MO BIO Laboratories, Solana Beach, CA, USA), according the manufacturer protocol, from 2 ml of culture. Species abundance was measured by quantitative PCR using primers for C. acetobutylicum [22] and W. succinogenes (forward primer (5’→3’): CCACACGACACTACCCTCAC and reverse primer (5’→3’): CGGTGCTTACGATTCCCAGT) on a 7300 Real Time PCR System (AB Applied Biosystems, CA, USA). Both primer sets are selective for the target species and the quantification of genes is corrected with the number of genes on the genome. The following PCR program was used for C. acetobutylicum quantification: a start-up step of 50°C for 2 minutes, followed by an initial denaturation step at 94°C for 10 minutes, 40 cycles of 94°C for 30 seconds, 54°C for 30 seconds and 72°C for 40 seconds, plus a final elongation step at 72°C for 5 minutes. For quantification of W. succinogenes we used: a start-up step of 50°C for 2 minutes, followed by an initial denaturation step at 94°C for 10 minutes, 45 cycles of 94°C for 10 seconds, and 60°C for 30 seconds.

Models and software

All simulations were done using CBMPy 0.7.0 (http://cbmpy.sourceforge.net/), a Python-based software package. The genome-scale metabolic model of C. acetobutylicum [13], a coarse-grained metabolic model and a genome-scale metabolic model of W. succinogenes (see supplements), using the SBML level 3 standards [23], were used to perform the simulations (All models and scripts are available on https://sourceforge.net/projects/cbmpy/files/publications/data/2017_Hanemaaijer/). The method to reconstruct the genome-scale metabolic model of W. succinogenes is described in the supplemental material. Dynamic flux balance analysis simulations were based on the method by Zhuang et al. [24], but implemented in Python using the LSODA integrator provided by the SciPy scientific library [25]. The CPLEX-solver (academic license) from IBM^® was used for optimization of the metabolic networks. The method for model reconstruction and co-culture growth simulations are described in the supplements. The Flux Variability Analysis was performed according Mahadevan et al. [26].

Results

Design of co-culture and media composition that gives rise to different metabolic interactions

We designed a synthetic co-culture of C. acetobutylicum and W. succinogenes, serving as a model system for anaerobic interspecies hydrogen transfer. These species are chosen, because they allow for tunable metabolic interactions by changes made to the medium composition, in particular the nitrogen source (Fig 1).

Download:

Fig 1. Source of nitrogen is expected to influence the metabolic interactions between C. acetobutylicum and W. succinogenes.

A shows the expected interactions between C. acetobutylicum and W. succinogenes in the presence of NH₄⁺ and NO₃^-. B shows the interactions in the presence of only NH₄⁺, C in the presence of only NO₃^- and D shows the expected interactions in the presence of only N₂. Although we expect certain interactions in the co-culture during the various environmental conditions, we do not know at what rate the metabolites are exchanged.

https://doi.org/10.1371/journal.pone.0173183.g001

We choose C. acetobutylicum because it produces H₂ from glucose, which can be consumed by W. succinogenes to generate energy for the reduction of NO₃^- to NO₂^- or, further, to NH₄⁺ when the NO₃^- concentration is low. We choose a high NO₃^- concentration in our experimental set-up to prevent this. As an alternative to H₂ W. succinogenes can use formate. Since C. acetobutylicum cannot use NO₃^- as a nitrogen-source it has to fix N₂ gas when NO₃^- is the sole nitrogen source in pure culture [27]. However, in co-culture C. acetobutylicum can use NO₂^- produced by W. succinogenes. Alternatively, when NH₄⁺ is present, both organisms can use it as nitrogen source, but W. succinogenes would still not grow because it lacks an energy source. Varying the nitrogen source also has additional effects. When NH₄⁺ and NO₃^- is present, C. acetobutylicum is independent of W. succinogenes and supplies H₂ to W. succinogenes. C. acetobutylicum becomes dependent on W. succinogenes when only NO₃^- is present. Then W. succinogenes cross feeds NO₂^- to C. acetobutylicum, which, in turn, supplies H₂ that facilitates NO₃^- to NO₂^- reduction by W. succinogenes. Thus we can design co-cultures with uni- and bidirectional metabolic interactions. We studied four different conditions; the associated, expected metabolic interactions are shown in Fig 1.

Experimental data of four growth conditions

The co-culture was grown in four different environmental conditions, in which we varied the nitrogen-source to study its influence on the state of the community, outlined in Fig 1. The media either contained NH₄⁺ and NO₃^-, only NH₄⁺, only NO₃^- or neither. N₂ was the only nitrogen-source in the absence of NH₄⁺ and NO₃^-, with the expectation that only C. acetobutylicum can grow. Across conditions, the biomass abundances and metabolite concentrations were measured in time to study the impact of the nitrogen-source on the co-culture dynamics.

The different conditions influenced the growth and metabolic profiles of the co-culture (Fig 2). The co-culture grew slower in the N₂ condition and led to a 2-fold lower optical density than the NO₃^- condition (Fig 2A). In the N₂ condition, C. acetobutylicum fixes N₂ gas to acquire nitrogen, which is an energetically unfavourable process [28, 29]. The N₂ condition was also the only condition where glucose was not completely consumed at the end of the experiment. A drop in the optical density was observed in two conditions at the end of the experiment, because of induction of sporulation by C. acetobutylicum. We only detected H₂ when NO₃^- was absent (Fig 2B), suggesting that W. succinogenes consumed all the H₂ produced by C. acetobutylicum in the presence of NO₃^- (since no formate was detected). The final amounts of butyrate, acetate, CO₂ (produced by C. acetobutylicum) and reduced NO₃^- (by W. succinogenes) varied between the different conditions (Fig 2C). In the absence of NO₃^-, more butyrate was formed out of glucose, at the expense of acetate. This indicates a bidirectional interaction: in the presence of NO₃^-, when W. succinogenes grows, the H₂ concentration is kept low by W. succinogenes, such that C. acetobutylicum has a different fermentation profile due to alleviation of H₂-inhibition, which is consistent with literature [30]. This is therefore proof of a bidirectional interaction, because C. acetobutylicum behaves as in pure culture in the absence of NO₃^-, with high H₂ and butyrate concentrations. Note that H₂ accumulates in the absence of NO₃^-, because the growth rate of W. succinogenes decreases (Fig 2E). C. acetobutylicum produced 30% more CO₂ from glucose in the ‘N₂’ condition relative to the other conditions.

Download:

Fig 2. Nitrogen-source alters the phenotypic behaviour of the co-culture in batch experiments.

A shows optical density measurements at 600 nm. B is the measured H₂ concentrations. C are fermentation products per consumed glucose at the the end of the cultivation period and is the sum of the production and consumption of the species in the co-culture. The NO₃^- is the amount of NO₃^- reduced. D and E are the quantified gene copies concentration per mL of C. acetobutylicum and W. succinogenes. F shows the ratios of C. acetobutylicum and W. succinogenes based on gene-copy number. This is to show that both strains show balanced growth during batch cultivation. The copy numbers are corrected for the amount of genes on the genome.

https://doi.org/10.1371/journal.pone.0173183.g002

Fig 2F indicates time periods during which a constant biomass ratio of the two species were observed. Only in the presence of NH₄⁺, W. succinogenes was outgrown by C. acetobutylicum after 13.5 hours, because no NO₃^- is available for W. succinogens to convert H₂. In the other cases we found constant biomass ratios, always after 13.5 hours of cultivation, indicating that the two microorganisms grew at the same specific rate. Equal specific growth rates indicate that the growth of the two organisms was coupled. Such growth coupling can either result from bidirectional exchange of growth-supporting metabolites, or a unidirectional exchange from the slower-growing to the faster-growing microorganism. A bidirectional exchange was expected in the presence of only NO₃^-, whereas in the presence of NH₄⁺ and NO₃^- a unidirectional exchange was expected. Coupling occurs because the microorganism that is capable of a higher growth rate becomes limited by the other organism, which cannot supply its resources fast enough.

Which metabolites are exchanged and at what rate across those conditions cannot be deduced directly from the experimental data. For instance, we cannot conclude whether C. acetobutylicum only consumed NH₄⁺, or that NO₂^- was also consumed in the presence of NH₄⁺ and NO₃^-. Also the H₂ production rate of C. acetobutylicum remains unknown, as we could not detect H₂ in two of the four conditions. Even though we can conclude that the species do exchange metabolites with each other, as they grew equally fast, we are not sure which metabolites are being exchanged let alone at what rate. To address those questions we analysed our experimental data with metabolic models of the two species.

Stoichiometric models of co-culture metabolism and growth

We used stoichiometric metabolic models of C. acetobutylicum and W. succinogenes to fit the experimental data for the inference of metabolic interactions between the species during different environmental conditions. The metabolic models only describe the stoichiometry of metabolic reactions, based on genomic and biochemical information, and do therefore not contain any kinetic parameters of all the individual enzymes. In addition, we took the growth rates and the biomass abundances of the microorganisms into account to capture the metabolism and the dynamics of the co-culture, using an existing method called dynamic flux balance analysis (dFBA) [31], which extends the stoichiometric model by adding kinetic information of growth-limiting substrate importers in the form of an irreversible Michaelis-Menten equation (), where v_upt represents the uptake reaction, V_max the maximum rate of the uptake reaction, [S] the substrate concentration and K_m the substrate affinity. At each time step in dFBA the biomass production rate for each organism is independently optimised based on the specific substrate uptake rate. Implementing such dynamic uptake reaction enables the model to respond to changes in external concentrations.

For the simulations we use a genome-scale metabolic model for the metabolism of C. acetobutylicum and for W. succinogenes a coarse-grained model, which lumps metabolic segments into single reactions. We used a coarse-grained model as we expected that we could model our results without the need for full genome-scale reconstruction. We did do the reconstruction to verify that this was indeed the case (see supplemental material). Having a coarse-grained model facilitates computation and interpretation, especially if these methods are applied to larger and more complex microbial ecosystems.

Conversely, the detailed model of C. acetobutylicum has the advantage that it captures changes in its metabolism, as we observed that final acetate and butyrate concentrations, produced by C. acetobutylicum changes when the growth medium was varied. Another advantage is that it describes the metabolic flexibility associated with H₂ production by C. acetobutylicum and leads to oxidization of ferredoxin. Ferredoxin is a conserved moiety whose total concentration does not change (only its redox state), and when it is oxidized by another reaction, it will result in decreased H₂ production. Since C. acetobutylicum contains several reactions that can oxidize or reduce ferredoxin, the H₂ production flux is flexible, influencing the community interactions that we aim to quantify. For instance, NO₂^- uptake decreases the amount of H₂ production, because both reactions require ferredoxin, which results in potential interesting dynamics between the two species, because W. succinogenes provides NO₂^-, but requires H₂ which is produced by C. acetobutylicum. Several genome-scale metabolic models of C. acetobutylicum are available [32, 33]. We decided to use the model created by Salimi and co-workers [13], because this model was already successfully used for dFBA simulations.

We constructed a coarse-grained model of the core metabolism of W. succinogenes, containing eight reactions (S1 Table) that was detailed enough to fit the exchange fluxes to experimental data; the data indicated that C. acetobutylicum shows most metabolic flexibility. The ATP yield of NO₃^- reduction, when H₂ is oxidized, was taken from literature [34]. We fitted the elemental composition of W. succinogenes biomass. It can use acetate as a carbon source [35], but with an unknown biomass yield. Similarly the NH₄⁺ and ATP requirements for the production of one unit of biomass are unclear. We estimated all these values by varying them and used those values that resulted in the best fit with the experimental data. An advantage of working with coarse-grained models is that the model complexity is greatly reduced and the degrees of freedom are limited, whereas in genome-scale metabolic models these can increase rapidly [36–38]. However, we also constructed a draft genome-scale metabolic model of W. succinogenes to compare the results from the simulations with this model with the coarse-grained model (S1 File).

The two metabolic models were coupled through exchange fluxes for H₂, NO₂^- and acetate. The co-culture model was manually fitted to the experimental data, using dFBA, to infer the metabolic exchange fluxes in the microbial community. Manually fitting was done by changing the kinetic parameters (S2 Table), biomass composition of W. succinogenes (S1 Table) and ratio of acetate/butyrate production and NH₄⁺/NO₂^- consumption until the simulation agreed most with the experimental data. ame

The metabolic models can simulate the experimental data

Once the simulations agree with the experimental data, the corresponding metabolic fluxes can be extracted from the simulations, including the exchange fluxes between the two species. This allows us to infer the metabolic interactions between species in a microbial community.

To fit the model with the experiments, we varied the parameters of the models (S2 Table) to find an optimal fit. These parameters influence the rate of consumption and production of the metabolites in the co-cultures and can be different between the various environmental conditions. For instance, the V_max of glucose uptake in C. acetobutylicum changes between the different conditions (S2 Table). Although the parameters can change between the conditions, we did not change the metabolic models itself. The structure of the metabolic models, therefore, remain identical between the different environmental conditions.

For C. acetobutylicum, we had to constrain the ratio of butyrate and acetate production fluxes and the ratio of the NH₄⁺ and NO₂^- consumption to fit the experimental data (S2 Table). Removing those ratios resulted in incorrect simulation of the metabolic profiles, which is the result of apparent suboptimal behaviour of C. acetobutylicum. One reason for this could be that H₂ inhibits certain reactions and results in suboptimal behavior. Several studies showed that H₂ affects the metabolism of Clostridia species [30, 39]. An other reason could be that there is a yield/rate trade-off and a high growth rate would result in a suboptimal biomass yield as observed for other organisms. Flux balance analysis will always find the flux distribution that maximizes this yield, not the rate. We also had to fit V_max and K_m of glucose uptake by C. acetobutylicum and the V_max and K_m of H₂ uptake by W. succinogenes (S2 Table).

The simulation with the metabolic models are consistent with the experimental data, except that W. succinogenes grew in the absence of NO₃^-, which was not predicted with the model (Fig 3). This is most likely caused by carry-over of NO₃^- from the pre-culture and by cell lysis. However, these phenomena were not taken into account in the simulations of the metabolic model and therefore, no growth of W. succinogenes was predicted with the metabolic models (Fig 4B and 4D). Simulations where the genome-scale metabolic model of W. succinogenes was used were similar relative to the simulations with the coarse-grained model of W. succinogenes (S1 File).

Download:

Fig 3. The dFBA simulations agree mostly with the experimental data for the four different cultivation conditions and were used to infer the metabolic fluxes.

The metabolite profile plots contain error bars, but the other subplots not. The biomass ratios are based on the gene-copy data and the carbon balance consisted of the measured metabolites that contained carbon. The remaining missing carbon is assumed to be incorporated into biomass.

https://doi.org/10.1371/journal.pone.0173183.g003

Download:

Fig 4. Growth of W. succinogenes was limited by the growth-rate of C. acetobutylicum, based on the inferred growth-rates.

Balanced growth was observed after 10 hours of cultivation in the NH₄⁺ + NO₃^- (A) and NO₃^- conditions (C). No growth of W. succinogenes was simulated in the NH₄⁺ (B) and N₂ conditions (D).

https://doi.org/10.1371/journal.pone.0173183.g004

Influence of environmental conditions on the metabolic exchange fluxes between C. acetobutylicum and W. succinogenes

Simulations with the fitted model indicates that the growth rates of C. acetobutylicum and W. succinogenes are dynamic during the experiment (Fig 4). The biomass yields on substrate, however, remain constant during the whole experiment (S2 Fig). The dynamic behaviour of the community resulted therefore solely from the interplay between the dynamic concentrations of the environmental nutrients and products and the biomass abundances.

W. succinogenes grew faster than C. acetobutylicum in the first 10 hours in the presence of NH₄⁺ and NO₃^-. This is possible, because C. acetobutylicum is more abundant than W. succinogenes and can provide excess H₂ despite the slower growth of C. acetobutylicum. After 10 hours, both species grew at the same specific rate. This is reflected in the experimental data of the biomass ratios of C. acetobutylicum/W. succinogenes, which decreased in the first hours of the experiment suggesting that W. succinogenes grows faster than C. acetobutylicum. This indicates that C. acetobutylicum determined the growth-rate of the community, because H₂ consumption of W. succinogenes was limited by the H₂ production of C. acetobutylicum after ten hours of cultivation.

Different nitrogen sources lead to changed metabolic interactions between C. acetobutylicum and W. succinogenes (Fig 5). The H₂ production rate of C. acetobutylicum, during cultivation in the presence of NO₃^-, differed between the conditions: 1.10 and 0.94 mol/mol H₂ per consumed glucose, for the ‘NH₄⁺ + NO₃^-’ and ‘NO₃^-’ condition, respectively. H₂ production was reduced by 17% in the ‘NO₃^-’ condition relative to the ‘NH₄⁺ + NO₃^-’ condition. H₂ is associated with oxidization of ferredoxin and this ferredoxin is also oxidized when NO₂^- is reduced to NH₄⁺. Since ferredoxin is a conserved moiety, an increase in the NO₂^- reduction flux causes a lower H₂ production flux, which occurred in the ‘NO₃^-’ condition. Therefore, there is an interesting dynamic between C. acetobutylicum and W. succinogenes when NO₃^- is in the medium, because C. acetobutylicum can use its ferredoxin for the production of H₂, but also for the consumption of NO₂^-, which is produced by W. succinogenes. We show that the nitrogen source can have implications for interspecies hydrogen transfer in a community. An oxidized nitrogen source, in this case NO₂^-, would decrease H₂ transfer and therefore NH₄⁺ is the preferred substrate for maximization of interspecies hydrogen transfer. Contrary to our expectations, a combination of NH₄⁺ and NO₂^- is predicted to be utilized by C. acetobutylicum in the NH₄⁺ + NO₃^- condition. One explanation for C. acetobutylicum’s utilization of NO₂^- is to detoxify the medium. High NO₂^- concentrations inhibit growth and utilizing NO₂^- as a nitrogen source could be a detoxification strategy [40].

Download:

Fig 5. Nitrogen-source also had an impact on the H₂ and the NO₂^- exchange rates.

In A the calculated flux-values for C. acetobutylicum are normalized for the glucose uptake. B shows the calculated flux-values for W. succinogenes and are normalized for the H₂ uptake. Note that H₂ production by C. acetobutylicum is not equal to the specific H₂ uptake rate of W. succinogenes as the biomass abundances should be taken into account, which were not equal during the experiment.

https://doi.org/10.1371/journal.pone.0173183.g005

In contrast to the inferred fluxes of C. acetobutylicum, the inferred fluxes of W. succinogenes were equal between the ‘NH₄⁺ + NO₃^-’ and ‘NO₃^-’ conditions (Fig 5B), suggesting that W. succinogenes was insensitive to the change in environmental conditions. This is primarily due to the limited flexibility of the model of W. succinogenes; however, the model fits with the experimental data, which suggests that the inferred fluxes are realistic. The simulations with the draft genome-scale metabolic model of W. succinogenes showed similar results relative to the coarse-grained model, except that the amount of consumed NO₃^- and produced NO₂^- per H₂ was different (S3 Fig). This is primarily caused by the nitrogen and energy requirements for the production of biomass and is different for both models. We expect that a more accurate biomass function for the genome-scale metabolic model of W. succinogenes would improve the results.

Intracellular fluxes of C. acetobutylicum can be inferred from community-level experimental data

In addition to calculation of the uptake and production fluxes, the genome-scale metabolic model of C. acetobutylicum informs us about intracellular flux values across the different conditions. This gives additional information about the active metabolic pathways. Due to the large number of reactions in the genome-scale metabolic model, the predicted values of the intracellular fluxes in the metabolic network are likely underdetermined, because we have too few experimental data to constrain those fits sufficiently.

We therefore performed a Flux Variability Analysis (FVA) [26], to check whether alternative intracellular flux values of the C. acetobutylicum metabolic network could fit the experimental data equally well. This analysis indicates a small flexibility, only 41 reactions of the 744 reactions had variable values (5.5%) (https://sourceforge.net/projects/cbmpy/files/publications/data/2017_Hanemaaijer/). Those 41 reactions are involved in either malate cycling, proline and nucleotide biosynthesis, and carry very small fluxes relative to the fluxes in C. acetobutylicum’s core-metabolism. Two reactions that did carry a high flux are involved in ferredoxin cycling, but had no impact on the H₂ production flux. There was also no variability identified in other uptake or production fluxes. These results suggest that in our system, the calculated flux distribution within C. acetobutylicum can be robustly predicted from community-level experimental data alone.

Modelling predicts that co-culture behaviour can be enhanced

Besides providing a method for inferring experimental fluxes from experimental data, the metabolic model can also be used to predict the optimal behaviour of the community when not constrained with experimental data. Optimisation informs us about optimal metabolic coupling between the microorganisms when growth of both species are optimised.

Here, we consider the co-culture as optimal when the growth rates of both organisms are maximal. In dFBA this means that we optimise the flux towards biomass of both organisms on the supplied nutrients to the community. We determined the metabolic fluxes (internal and exchange), growth rate and biomass abundances corresponding to the optimal state in the ‘NH₄⁺ + NO₃^-’ and ‘NO₃^-’ condition when we released the constraints on the butyrate/acetate production ratio and the NH₄⁺/NO₂^- consumption ratio of C. acetobutylicum (Fig 6).

Download:

Fig 6. The ‘NO₃^-’ condition behaved more optimal than the ‘NH₄⁺ + NO₃^-’ condition.

The calculated experimental specific uptake and consumption rates (mmol/(gDW⋅h)) for C. acetobutylicum (white bars) and the in-silico optimal rates (gray bars) were compared. All specific uptake and production fluxes were normalized to the specific glucose uptake rate.

https://doi.org/10.1371/journal.pone.0173183.g006

The model predicted that an optimal C. acetobutylicum would have a changed metabolism relative to the experimental data. For instance, acetate rather than butyrate should be produced when C. acetobutylicum grows optimally. This behaviour was seen both in the ‘NH₄⁺ + NO₃^-’ and the ‘NO₃^-’ condition. The explanation is that more ATP can be generated from acetate production than from butyrate production. An optimal C. acetobutylicum would not consume NO₂^- when growing at the ‘NH₄⁺ + NO₃^-’ condition, because NO₂^- reduction by nitrite reductase decreases growth of C. acetobutylicum according to the metabolic model. As explained, the amount of H₂ produced is inversely related with the amount of NO₂^- taken up. As a consequence, 159% more H₂ was produced in the optimal scenario relative to the experiments in the ‘NH₄⁺ + NO₃^-’ condition. A smaller increase in H₂ production (+23%) was shown in the ‘NO₃^-’ condition, suggesting that the interactions between the two species were not optimal in the experiments at those two conditions. However, the experimental growth-rate of C. acetobutylicum was close to its maximal, computed value. For the ‘NH₄⁺ + NO₃^-’ the experimental growth rate was 87% of the optimal value, whereas for the ‘NO₃^-’ this was 94%. These results suggest that the experimental behaviour of the community with imposed bidirectional interactions is closer to the optimal metabolic behaviour than a community with unidirectional interactions.

Discussion

In this study we showed that the combination of a stoichiometric metabolic model with measured community-level concentrations and biomass abundances allows for the inference of metabolic interactions in a microbial community. This approach is complementary with current approaches that are based on correlation or co-occurrence networks, which predicts the community-interaction structure for large complex ecosystems [5, 41, 42]. These networks identify co-occurrence of species, suggesting that certain species in an ecosystem form a positive relationship with each other. However, the mechanisms behind the interactions remain unknown; here we infer them with metabolic models.

Stoichiometric models of metabolism make inference of metabolic exchange fluxes possible because they relate the uptake and production of environmental metabolites, to intracellular metabolic fluxes, organism abundance and growth rate. Such models are therefore integrating heterogeneous data types from which we make inferences about community behaviour. In addition to experimental data of species abundances and metabolite concentrations, molecular data can be used in the construction of stoichiometric metabolic models [43–45]. It is the integrative capability of stoichiometric models that makes them so useful for addressing current challenges in microbial ecology.

Metabolic models have been used to predict the behaviour of industrially relevant co-cultures and the type of interactions required to predict this behaviour [13, 14]. We show that this approach also works for ecologically relevant communities, where interspecies hydrogen transfer is the dominating interaction. Additionally, not only the type of interactions is identified, but we also infer the size of the interactions from the experimental data. With this method we show that the nitrogen source plays an important role in the rate of H₂ exchange in our co-culture. Recently, Embree et al. [46] also inferred metabolic interactions of a microbial community, but optimized each species individually as a steady-state simulation and did not combine all metabolic models into one community. We, on the other hand, simulated the dynamics of the co-culture, because we fitted data from a dynamic system where species could change their metabolic behaviour during the experiment, which is closer to natural environments.

We have shown that our approach works for a co-culture and a next step would be the inference of metabolic interactions in a complex ecosystem. A complex community will have more interactions between the members of the community. In addition to beneficial relationships, competition can also play an important role in microbial communities [1]. The advantage of our dFBA based approach is that these type of interactions can also be inferred. Moreover, our results indicate that it is not even necessary to use full-blown genome-scale stoichiometric models for all members of a complex community as a coarse-grained model of W. succinogenes sufficed for the inference of metabolic interactions, as we recently suggested [7]. More detail can always be added later to refine predictions or resolve discrepancies, depending on the research question and data at hand.

Even though metabolic interactions can be inferred from experimental data, not all data and systems are equally amenable to a modeling approach. Firstly, a community with growing micro-organisms is advantageous. Growing cells have an active metabolism and increase their biomass in time, while interacting with other species. While non-growing species may also have an active metabolism, they do not increase in biomass, causing them to have lower fluxes, which makes it difficult to assign flux-values to such organisms. Secondly, controlled environments allow for more targeted perturbations of the community. Manipulation of the environment results in different dominating processes in the ecosystem [47]. In addition, a basic understanding of the core metabolism of the studied organisms is required. In our system, we had already an idea what type of interactions would take place, but this will not always be possible with complex ecosystems with poorly studied microorganisms. In such cases, metabolic reconstruction based on complete genomes is the best option to get a hint at the basic metabolic capacities. Lastly, a medium that is well-defined and minimal, is advantageous. Metabolites, which are turned over, are more easily identified in such media. A medium which, for instance, contains yeast extract provides all types of vitamins, amino acids and carbon compounds to the community members, and could mask or compromise the metabolic interactions in the ecosystem. With improved analytics, such challenges may be overcome, and we believe in view of the complexity of microbial communities, metabolic models are a promising -and arguable the only- tool for the inference of the metabolic interactions of co-cultures.

To conclude, in our eyes metabolic models are an indispensable tool for the inference of the metabolic interactions, as illustrated with our co-culture study. We have shown that the nitrogen-source quantitatively influences the interspecies hydrogen transfer, which was not feasible with experimental data alone. Moreover, the metabolic models can be used for more purposes than inference of metabolic interactions. They can also be used to infer the intracellular fluxes of species, which informs us about the metabolic adaptations of microorganisms, during different conditions. The models can also be used to investigate whether species in a community grow optimally or whether they are limited by external factors. This is particularly interesting for communities used for industrial processes, such as the anaerobic digestion process, where high biogas yields are desired. Finally, metabolic models could indicate strategies to steer the community in a desired direction, guiding community engineering approaches. Thus metabolic models have great potential –are simply required– to study complex, microbial ecosystems.

Supporting information

S1 Material and Methods.

https://doi.org/10.1371/journal.pone.0173183.s001

(PDF)

S1 Table. A coarse-grained metabolic model of W. succinogenes was detailed enough to fit the experimental data.

All metabolites, except for H⁺ and H₂O, are element-balanced in the metabolic model.

https://doi.org/10.1371/journal.pone.0173183.s002

(PDF)

S2 Table. Parameters used for simulation of the co-culture.

https://doi.org/10.1371/journal.pone.0173183.s003

(PDF)

S1 Fig. The dFBA simulations, with the draft genome-scale metabolic model of W. succinogenes, also agree mostly with the experimental data for the four different cultivation conditions.

The metabolite profile plots contain error bars, but the other subplots not. The biomass ratios are based on the gene-copy data and the carbon balance consisted of the measured metabolites that contained carbon. The remaining missing carbon is assumed to be incorporated into biomass.

https://doi.org/10.1371/journal.pone.0173183.s004

(EPS)

S2 Fig. Product and substrate yields did not change during different growth rates.

The extracted specific uptake and production rates (mmol/(gDW⋅h)) are plotted against the simulated growth rate μ (h-1). All lines are lineair what suggests that the specific uptake and production rates of both species are scaled with the growth. The N₂ condition is not plotted, because the simulated growth rate was constant over the whole simulation.

https://doi.org/10.1371/journal.pone.0173183.s005

(EPS)

S3 Fig. Inference of the metabolic interactions using a draft genome-scale metabolic model of W. succinogenes showed only differences in the NO₂^- production and NO₃^- consumption per utilized H₂.

In A the calculated flux-values for C. acetobutylicum are normalized for the glucose uptake. B shows the calculated flux-values for W. succinogenes and are normalized for the H₂ uptake. Note that H₂ production by C. acetobutylicum is not equal to the specific H₂ uptake rate of W. succinogenes as the biomass abundances should be taken into account, which were not equal during the experiment.

https://doi.org/10.1371/journal.pone.0173183.s006

(EPS)

S1 File. Draft genome-scale metabolic model of W. succinogenes.

The reactions and the corresponding reversibility and bounds are listed in ‘reactions’. The metabolites-id and the corresponding names are listed in ‘metabolites’. All reactions and their corresponding substrates and products are listed in ‘network_react’, where in ‘network_metab’ all reactions per metabolite is listed. All metabolites and their corresponding MIRIAM information is listed in ‘miriam.’

https://doi.org/10.1371/journal.pone.0173183.s007

(XLS)

Acknowledgments

This work is financially supported by ZonMW and the Netherlands Genomics Initiative via a Zenith Horizon grant 40-41009-98-10038. BGO is supported by a BE-Basic Foundation grant F08.005.001.

Author Contributions

Conceptualization: MH WR FB BT.
Formal analysis: MH BO.
Funding acquisition: BT.
Investigation: MH.
Methodology: MH BO.
Project administration: WR BT.
Software: BO.
Supervision: WR BT.
Validation: MH.
Visualization: MH.
Writing – original draft: MH FB.
Writing – review & editing: MH BO FB BT.

References

1. Coyte KZ, Schluter J, Foster KR. The ecology of the microbiome: Networks, competition, and stability. Science. 2015;350(6261):663–666. pmid:26542567
- View Article
- PubMed/NCBI
- Google Scholar
2. De Vuyst L, Neysens P. The sourdough microflora: biodiversity and metabolic interactions. Trends in Food Science & Technology. 2005;16(1):43–56.
- View Article
- Google Scholar
3. Møller S, Sternberg C, Andersen JB, Christensen BB, Ramos JL, Givskov M, et al. In situ gene expression in mixed-culture biofilms: evidence of metabolic interactions between community members. Applied and environmental microbiology. 1998;64(2):721–732. pmid:9464414
- View Article
- PubMed/NCBI
- Google Scholar
4. Stams AJ. Metabolic interactions between anaerobic bacteria in methanogenic environments. Antonie van Leeuwenhoek. 1994;66(1-3):271–294. pmid:7747937
- View Article
- PubMed/NCBI
- Google Scholar
5. Faust K, Sathirapongsasuti JF, Izard J, Segata N, Gevers D, Raes J, et al. Microbial co-occurrence relationships in the human microbiome. PLoS Comput Biol. 2012;8(7):e1002606–e1002606. pmid:22807668
- View Article
- PubMed/NCBI
- Google Scholar
6. Barberán A, Bates ST, Casamayor EO, Fierer N. Using network analysis to explore co-occurrence patterns in soil microbial communities. The ISME journal. 2012;6(2):343–351. pmid:21900968
- View Article
- PubMed/NCBI
- Google Scholar
7. Hanemaaijer M, Röling WF, Olivier BG, Khandelwal RA, Teusink B, Bruggeman FJ. Systems modeling approaches for microbial community studies: from metagenomics to inference of the community structure. Frontiers in microbiology. 2015;6. pmid:25852671
- View Article
- PubMed/NCBI
- Google Scholar
8. Zengler K, Palsson BO. A road map for the development of community systems (CoSy) biology. Nature Reviews Microbiology. 2012;10(5):366–372. pmid:22450377
- View Article
- PubMed/NCBI
- Google Scholar
9. Röling WF, van Bodegom PM. Toward quantitative understanding on microbial community structure and functioning: a modeling-centered approach using degradation of marine oil spills as example. Frontiers in microbiology. 2014;5. pmid:24723922
- View Article
- PubMed/NCBI
- Google Scholar
10. Klitgord N, Segre D. Ecosystems biology of microbial metabolism. Current opinion in biotechnology. 2011;22(4):541–546. pmid:21592777
- View Article
- PubMed/NCBI
- Google Scholar
11. Levy R, Borenstein E. Metabolic modeling of species interaction in the human microbiome elucidates community-level assembly rules. Proceedings of the National Academy of Sciences. 2013;110(31):12804–12809. pmid:23858463
- View Article
- PubMed/NCBI
- Google Scholar
12. Stolyar S, Van Dien S, Hillesland KL, Pinel N, Lie TJ, Leigh JA, et al. Metabolic modeling of a mutualistic microbial community. Molecular systems biology. 2007;3(1). pmid:17353934
- View Article
- PubMed/NCBI
- Google Scholar
13. Salimi F, Zhuang K, Mahadevan R. Genome-scale metabolic modeling of a clostridial co-culture for consolidated bioprocessing. Biotechnology journal. 2010;5(7):726–738. pmid:20665645
- View Article
- PubMed/NCBI
- Google Scholar
14. Hanly TJ, Henson MA. Dynamic flux balance modeling of microbial co-cultures for efficient batch fermentation of glucose and xylose mixtures. Biotechnology and bioengineering. 2011;108(2):376–385. pmid:20882517
- View Article
- PubMed/NCBI
- Google Scholar
15. Heinken A, Sahoo S, Fleming RM, Thiele I. Systems-level characterization of a host-microbe metabolic symbiosis in the mammalian gut. Gut microbes. 2013;4(1):3–15. pmid:23022739
- View Article
- PubMed/NCBI
- Google Scholar
16. Klitgord N, Segre D. Environments that induce synthetic microbial ecosystems. PLoS Comput Biol. 2010;6(11):e1001002. pmid:21124952
- View Article
- PubMed/NCBI
- Google Scholar
17. Khandelwal RA, Olivier BG, Röling WF, Teusink B, Bruggeman FJ. Community Flux Balance Analysis for Microbial Consortia at Balanced Growth. PloS one. 2013;8(5):e64567. pmid:23741341
- View Article
- PubMed/NCBI
- Google Scholar
18. Wittebolle L, Marzorati M, Clement L, Balloi A, Daffonchio D, Heylen K, et al. Initial community evenness favours functionality under selective stress. Nature. 2009;458(7238):623–626. pmid:19270679
- View Article
- PubMed/NCBI
- Google Scholar
19. Pawelczyk S, Abraham WR, Harms H, Müller S. Community-based degradation of 4-chorosalicylate tracked on the single cell level. Journal of microbiological methods. 2008;75(1):117–126. pmid:18593643
- View Article
- PubMed/NCBI
- Google Scholar
20. Patle S, Lal B. Ethanol production from hydrolysed agricultural wastes using mixed culture of Zymomonas mobilis and Candida tropicalis. Biotechnology letters. 2007;29(12):1839–1843. pmid:17657407
- View Article
- PubMed/NCBI
- Google Scholar
21. Yang J, Kim J, Skogley E, Schaff B. A simple spectrophotometric determination of nitrate in water, resin, and soil extracts. Soil Science Society of America Journal. 1998;62(4):1108–1115.
- View Article
- Google Scholar
22. Lee SM, Cho MO, Um Y, Sang BI. Development of Real-Time PCR primer and probe sets for detecting degenerated and non-degenerated forms of the butanol-producing bacterium Clostridium acetobutylicum ATCC 824. Applied biochemistry and biotechnology. 2010;161(1-8):75–83. pmid:19798472
- View Article
- PubMed/NCBI
- Google Scholar
23. Olivier BG, Bergmann FT. The Systems Biology Markup Language (SBML) Level 3 Package: Flux Balance Constraints. Journal of integrative bioinformatics. 2015;12(2):269. pmid:26528567
- View Article
- PubMed/NCBI
- Google Scholar
24. Zhuang K, Izallalen M, Mouser P, Richter H, Risso C, Mahadevan R, et al. Genome-scale dynamic modeling of the competition between Rhodoferax and Geobacter in anoxic subsurface environments. The ISME journal. 2010;5(2):305–316. pmid:20668487
- View Article
- PubMed/NCBI
- Google Scholar
25. Jones E, Oliphant T, Peterson P. {SciPy}: Open source scientific tools for {Python}. 2014;.
26. Mahadevan R, Schilling C. The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metabolic engineering. 2003;5(4):264–276. pmid:14642354
- View Article
- PubMed/NCBI
- Google Scholar
27. Wilson P, Peterson W, Fred E. The relationship between the nitrogen and carbon metabolism of Clostridium acetobutylicum. Journal of bacteriology. 1930;19(4):231. pmid:16559425
- View Article
- PubMed/NCBI
- Google Scholar
28. Rosenblum ED, Wilson P. Fixation of isotopic nitrogen by Clostridium. Journal of bacteriology. 1949;57(4):413. pmid:18120606
- View Article
- PubMed/NCBI
- Google Scholar
29. Chen J, Toth J, Kasap M. Nitrogen-fixation genes and nitrogenase activity in Clostridium acetobutylicum and Clostridium beijerinckii. Journal of Industrial Microbiology and Biotechnology. 2001;27(5):281–286. pmid:11781802
- View Article
- PubMed/NCBI
- Google Scholar
30. Chung KT. Inhibitory effects of H2 on growth of Clostridium cellobioparum. Applied and environmental microbiology. 1976;31(3):342–348. pmid:779644
- View Article
- PubMed/NCBI
- Google Scholar
31. Mahadevan R, Edwards JS, Doyle FJ III. Dynamic Flux Balance Analysis of Diauxic Growth in Escherichia coli. Biophysical journal. 2002;83(3):1331–1340. pmid:12202358
- View Article
- PubMed/NCBI
- Google Scholar
32. Senger RS, Papoutsakis ET. Genome-scale model for Clostridium acetobutylicum: Part I. Metabolic network resolution and analysis. Biotechnology and bioengineering. 2008;101(5):1036. pmid:18767192
- View Article
- PubMed/NCBI
- Google Scholar
33. Lee J, Yun H, Feist AM, Palsson BØ, Lee SY. Genome-scale reconstruction and in silico analysis of the Clostridium acetobutylicum ATCC 824 metabolic network. Applied microbiology and biotechnology. 2008;80(5):849–862. pmid:18758767
- View Article
- PubMed/NCBI
- Google Scholar
34. Bokranz M, Katz J, Schröder I, Roberton AM, Kröger A. Energy metabolism and biosynthesis of Vibrio succinogenes growing with nitrate or nitrite as terminal electron acceptor. Archives of Microbiology. 1983;135(1):36–41.
- View Article
- Google Scholar
35. Schumacher W, Kroneck PM, Pfennig N. Comparative systematic study on “Spirillum” 5175, Campylobacter and Wolinella species. Archives of microbiology. 1992;158(4):287–293.
- View Article
- Google Scholar
36. Reed JL, Famili I, Thiele I, Palsson BO. Towards multidimensional genome annotation. Nature Reviews Genetics. 2006;7(2):130–141. pmid:16418748
- View Article
- PubMed/NCBI
- Google Scholar
37. Teusink B, Wiersma A, Molenaar D, Francke C, de Vos WM, Siezen RJ, et al. Analysis of growth of Lactobacillus plantarum WCFS1 on a complex medium using a genome-scale metabolic model. Journal of Biological Chemistry. 2006;281(52):40041–40048. pmid:17062565
- View Article
- PubMed/NCBI
- Google Scholar
38. Kelk SM, Olivier BG, Stougie L, Bruggeman FJ. Optimal flux spaces of genome-scale stoichiometric models are determined by a few subnetworks. Scientific reports. 2012;2. pmid:22896812
- View Article
- PubMed/NCBI
- Google Scholar
39. Van Andel JG, Zoutberg GR, Crabbendam PM and Breure AM. Glucose fermentation by Clostridium butyricum grown under a self generated gas atmosphere in chemostat culture. Applied microbiology and biotechnology. 1985;23(1):21–26.
- View Article
- Google Scholar
40. Ávila M, Gómez-Torres N, Hernández M, Garde S. Inhibitory activity of reuterin, nisin, lysozyme and nitrite against vegetative cells and spores of dairy-related Clostridium species. International journal of food microbiology. 2014;172:70–75. pmid:24361835
- View Article
- PubMed/NCBI
- Google Scholar
41. Berry D, Widder S. Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Frontiers in microbiology. 2014;5. pmid:24904535
- View Article
- PubMed/NCBI
- Google Scholar
42. Zhou J, Deng Y, Luo F, He Z, Yang Y. Phylogenetic molecular ecological network of soil microbial communities in response to elevated CO2. MBio. 2011;2(4):e00122–11. pmid:21791581
- View Article
- PubMed/NCBI
- Google Scholar
43. Åkesson M, Förster J, Nielsen J. Integration of gene expression data into genome-scale metabolic models. Metabolic engineering. 2004;6(4):285–293. pmid:15491858
- View Article
- PubMed/NCBI
- Google Scholar
44. Maarleveld TR, Boele J, Bruggeman FJ, Teusink B. A data integration and visualization resource for the metabolic network of Synechocystis sp. PCC 6803. Plant physiology. 2014; p. pp–113.
45. Yizhak K, Benyamini T, Liebermeister W, Ruppin E, Shlomi T. Integrating quantitative proteomics and metabolomics with a genome-scale metabolic network model. Bioinformatics. 2010;26(12):i255–i260. pmid:20529914
- View Article
- PubMed/NCBI
- Google Scholar
46. Embree M, Liu JK, Al-Bassam MM, Zengler K. Networks of energetic and metabolic interactions define dynamics in microbial communities. Proceedings of the National Academy of Sciences. 2015;112(50):15450–15455. pmid:26621749
- View Article
- PubMed/NCBI
- Google Scholar
47. Kraft B, Tegetmeyer HE, Sharma R, Klotz MG, Ferdelman TG, Hettich RL, et al. The environmental controls that govern the end product of bacterial nitrate respiration. Science. 2014;345(6197):676–679. pmid:25104387
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Coyte KZ, Schluter J, Foster KR. The ecology of the microbiome: Networks, competition, and stability. Science. 2015;350(6261):663–666. pmid:26542567
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. De Vuyst L, Neysens P. The sourdough microflora: biodiversity and metabolic interactions. Trends in Food Science & Technology. 2005;16(1):43–56.
View Article
Google Scholar

[6] View Article

[7] Google Scholar

[ref3] 3. Møller S, Sternberg C, Andersen JB, Christensen BB, Ramos JL, Givskov M, et al. In situ gene expression in mixed-culture biofilms: evidence of metabolic interactions between community members. Applied and environmental microbiology. 1998;64(2):721–732. pmid:9464414
View Article
PubMed/NCBI
Google Scholar

[9] View Article

[10] PubMed/NCBI

[11] Google Scholar

[ref4] 4. Stams AJ. Metabolic interactions between anaerobic bacteria in methanogenic environments. Antonie van Leeuwenhoek. 1994;66(1-3):271–294. pmid:7747937
View Article
PubMed/NCBI
Google Scholar

[13] View Article

[14] PubMed/NCBI

[15] Google Scholar

[ref5] 5. Faust K, Sathirapongsasuti JF, Izard J, Segata N, Gevers D, Raes J, et al. Microbial co-occurrence relationships in the human microbiome. PLoS Comput Biol. 2012;8(7):e1002606–e1002606. pmid:22807668
View Article
PubMed/NCBI
Google Scholar

[17] View Article

[18] PubMed/NCBI

[19] Google Scholar

[ref6] 6. Barberán A, Bates ST, Casamayor EO, Fierer N. Using network analysis to explore co-occurrence patterns in soil microbial communities. The ISME journal. 2012;6(2):343–351. pmid:21900968
View Article
PubMed/NCBI
Google Scholar

[21] View Article

[22] PubMed/NCBI

[23] Google Scholar

[ref7] 7. Hanemaaijer M, Röling WF, Olivier BG, Khandelwal RA, Teusink B, Bruggeman FJ. Systems modeling approaches for microbial community studies: from metagenomics to inference of the community structure. Frontiers in microbiology. 2015;6. pmid:25852671
View Article
PubMed/NCBI
Google Scholar

[25] View Article

[26] PubMed/NCBI

[27] Google Scholar

[ref8] 8. Zengler K, Palsson BO. A road map for the development of community systems (CoSy) biology. Nature Reviews Microbiology. 2012;10(5):366–372. pmid:22450377
View Article
PubMed/NCBI
Google Scholar

[29] View Article

[30] PubMed/NCBI

[31] Google Scholar

[ref9] 9. Röling WF, van Bodegom PM. Toward quantitative understanding on microbial community structure and functioning: a modeling-centered approach using degradation of marine oil spills as example. Frontiers in microbiology. 2014;5. pmid:24723922
View Article
PubMed/NCBI
Google Scholar

[33] View Article

[34] PubMed/NCBI

[35] Google Scholar

[ref10] 10. Klitgord N, Segre D. Ecosystems biology of microbial metabolism. Current opinion in biotechnology. 2011;22(4):541–546. pmid:21592777
View Article
PubMed/NCBI
Google Scholar

[37] View Article

[38] PubMed/NCBI

[39] Google Scholar

[ref11] 11. Levy R, Borenstein E. Metabolic modeling of species interaction in the human microbiome elucidates community-level assembly rules. Proceedings of the National Academy of Sciences. 2013;110(31):12804–12809. pmid:23858463
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref12] 12. Stolyar S, Van Dien S, Hillesland KL, Pinel N, Lie TJ, Leigh JA, et al. Metabolic modeling of a mutualistic microbial community. Molecular systems biology. 2007;3(1). pmid:17353934
View Article
PubMed/NCBI
Google Scholar

[45] View Article

[46] PubMed/NCBI

[47] Google Scholar

[ref13] 13. Salimi F, Zhuang K, Mahadevan R. Genome-scale metabolic modeling of a clostridial co-culture for consolidated bioprocessing. Biotechnology journal. 2010;5(7):726–738. pmid:20665645
View Article
PubMed/NCBI
Google Scholar

[49] View Article

[50] PubMed/NCBI

[51] Google Scholar

[ref14] 14. Hanly TJ, Henson MA. Dynamic flux balance modeling of microbial co-cultures for efficient batch fermentation of glucose and xylose mixtures. Biotechnology and bioengineering. 2011;108(2):376–385. pmid:20882517
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref15] 15. Heinken A, Sahoo S, Fleming RM, Thiele I. Systems-level characterization of a host-microbe metabolic symbiosis in the mammalian gut. Gut microbes. 2013;4(1):3–15. pmid:23022739
View Article
PubMed/NCBI
Google Scholar

[57] View Article

[58] PubMed/NCBI

[59] Google Scholar

[ref16] 16. Klitgord N, Segre D. Environments that induce synthetic microbial ecosystems. PLoS Comput Biol. 2010;6(11):e1001002. pmid:21124952
View Article
PubMed/NCBI
Google Scholar

[61] View Article

[62] PubMed/NCBI

[63] Google Scholar

[ref17] 17. Khandelwal RA, Olivier BG, Röling WF, Teusink B, Bruggeman FJ. Community Flux Balance Analysis for Microbial Consortia at Balanced Growth. PloS one. 2013;8(5):e64567. pmid:23741341
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref18] 18. Wittebolle L, Marzorati M, Clement L, Balloi A, Daffonchio D, Heylen K, et al. Initial community evenness favours functionality under selective stress. Nature. 2009;458(7238):623–626. pmid:19270679
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref19] 19. Pawelczyk S, Abraham WR, Harms H, Müller S. Community-based degradation of 4-chorosalicylate tracked on the single cell level. Journal of microbiological methods. 2008;75(1):117–126. pmid:18593643
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref20] 20. Patle S, Lal B. Ethanol production from hydrolysed agricultural wastes using mixed culture of Zymomonas mobilis and Candida tropicalis. Biotechnology letters. 2007;29(12):1839–1843. pmid:17657407
View Article
PubMed/NCBI
Google Scholar

[77] View Article

[78] PubMed/NCBI

[79] Google Scholar

[ref21] 21. Yang J, Kim J, Skogley E, Schaff B. A simple spectrophotometric determination of nitrate in water, resin, and soil extracts. Soil Science Society of America Journal. 1998;62(4):1108–1115.
View Article
Google Scholar

[81] View Article

[82] Google Scholar

[ref22] 22. Lee SM, Cho MO, Um Y, Sang BI. Development of Real-Time PCR primer and probe sets for detecting degenerated and non-degenerated forms of the butanol-producing bacterium Clostridium acetobutylicum ATCC 824. Applied biochemistry and biotechnology. 2010;161(1-8):75–83. pmid:19798472
View Article
PubMed/NCBI
Google Scholar

[84] View Article

[85] PubMed/NCBI

[86] Google Scholar

[ref23] 23. Olivier BG, Bergmann FT. The Systems Biology Markup Language (SBML) Level 3 Package: Flux Balance Constraints. Journal of integrative bioinformatics. 2015;12(2):269. pmid:26528567
View Article
PubMed/NCBI
Google Scholar

[88] View Article

[89] PubMed/NCBI

[90] Google Scholar

[ref24] 24. Zhuang K, Izallalen M, Mouser P, Richter H, Risso C, Mahadevan R, et al. Genome-scale dynamic modeling of the competition between Rhodoferax and Geobacter in anoxic subsurface environments. The ISME journal. 2010;5(2):305–316. pmid:20668487
View Article
PubMed/NCBI
Google Scholar

[92] View Article

[93] PubMed/NCBI

[94] Google Scholar

[ref25] 25. Jones E, Oliphant T, Peterson P. {SciPy}: Open source scientific tools for {Python}. 2014;.

[ref26] 26. Mahadevan R, Schilling C. The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metabolic engineering. 2003;5(4):264–276. pmid:14642354
View Article
PubMed/NCBI
Google Scholar

[97] View Article

[98] PubMed/NCBI

[99] Google Scholar

[ref27] 27. Wilson P, Peterson W, Fred E. The relationship between the nitrogen and carbon metabolism of Clostridium acetobutylicum. Journal of bacteriology. 1930;19(4):231. pmid:16559425
View Article
PubMed/NCBI
Google Scholar

[101] View Article

[102] PubMed/NCBI

[103] Google Scholar

[ref28] 28. Rosenblum ED, Wilson P. Fixation of isotopic nitrogen by Clostridium. Journal of bacteriology. 1949;57(4):413. pmid:18120606
View Article
PubMed/NCBI
Google Scholar

[105] View Article

[106] PubMed/NCBI

[107] Google Scholar

[ref29] 29. Chen J, Toth J, Kasap M. Nitrogen-fixation genes and nitrogenase activity in Clostridium acetobutylicum and Clostridium beijerinckii. Journal of Industrial Microbiology and Biotechnology. 2001;27(5):281–286. pmid:11781802
View Article
PubMed/NCBI
Google Scholar

[109] View Article

[110] PubMed/NCBI

[111] Google Scholar

[ref30] 30. Chung KT. Inhibitory effects of H2 on growth of Clostridium cellobioparum. Applied and environmental microbiology. 1976;31(3):342–348. pmid:779644
View Article
PubMed/NCBI
Google Scholar

[113] View Article

[114] PubMed/NCBI

[115] Google Scholar

[ref31] 31. Mahadevan R, Edwards JS, Doyle FJ III. Dynamic Flux Balance Analysis of Diauxic Growth in Escherichia coli. Biophysical journal. 2002;83(3):1331–1340. pmid:12202358
View Article
PubMed/NCBI
Google Scholar

[117] View Article

[118] PubMed/NCBI

[119] Google Scholar

[ref32] 32. Senger RS, Papoutsakis ET. Genome-scale model for Clostridium acetobutylicum: Part I. Metabolic network resolution and analysis. Biotechnology and bioengineering. 2008;101(5):1036. pmid:18767192
View Article
PubMed/NCBI
Google Scholar

[121] View Article

[122] PubMed/NCBI

[123] Google Scholar

[ref33] 33. Lee J, Yun H, Feist AM, Palsson BØ, Lee SY. Genome-scale reconstruction and in silico analysis of the Clostridium acetobutylicum ATCC 824 metabolic network. Applied microbiology and biotechnology. 2008;80(5):849–862. pmid:18758767
View Article
PubMed/NCBI
Google Scholar

[125] View Article

[126] PubMed/NCBI

[127] Google Scholar

[ref34] 34. Bokranz M, Katz J, Schröder I, Roberton AM, Kröger A. Energy metabolism and biosynthesis of Vibrio succinogenes growing with nitrate or nitrite as terminal electron acceptor. Archives of Microbiology. 1983;135(1):36–41.
View Article
Google Scholar

[129] View Article

[130] Google Scholar

[ref35] 35. Schumacher W, Kroneck PM, Pfennig N. Comparative systematic study on “Spirillum” 5175, Campylobacter and Wolinella species. Archives of microbiology. 1992;158(4):287–293.
View Article
Google Scholar

[132] View Article

[133] Google Scholar

[ref36] 36. Reed JL, Famili I, Thiele I, Palsson BO. Towards multidimensional genome annotation. Nature Reviews Genetics. 2006;7(2):130–141. pmid:16418748
View Article
PubMed/NCBI
Google Scholar

[135] View Article

[136] PubMed/NCBI

[137] Google Scholar

[ref37] 37. Teusink B, Wiersma A, Molenaar D, Francke C, de Vos WM, Siezen RJ, et al. Analysis of growth of Lactobacillus plantarum WCFS1 on a complex medium using a genome-scale metabolic model. Journal of Biological Chemistry. 2006;281(52):40041–40048. pmid:17062565
View Article
PubMed/NCBI
Google Scholar

[139] View Article

[140] PubMed/NCBI

[141] Google Scholar

[ref38] 38. Kelk SM, Olivier BG, Stougie L, Bruggeman FJ. Optimal flux spaces of genome-scale stoichiometric models are determined by a few subnetworks. Scientific reports. 2012;2. pmid:22896812
View Article
PubMed/NCBI
Google Scholar

[143] View Article

[144] PubMed/NCBI

[145] Google Scholar

[ref39] 39. Van Andel JG, Zoutberg GR, Crabbendam PM and Breure AM. Glucose fermentation by Clostridium butyricum grown under a self generated gas atmosphere in chemostat culture. Applied microbiology and biotechnology. 1985;23(1):21–26.
View Article
Google Scholar

[147] View Article

[148] Google Scholar

[ref40] 40. Ávila M, Gómez-Torres N, Hernández M, Garde S. Inhibitory activity of reuterin, nisin, lysozyme and nitrite against vegetative cells and spores of dairy-related Clostridium species. International journal of food microbiology. 2014;172:70–75. pmid:24361835
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref41] 41. Berry D, Widder S. Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Frontiers in microbiology. 2014;5. pmid:24904535
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref42] 42. Zhou J, Deng Y, Luo F, He Z, Yang Y. Phylogenetic molecular ecological network of soil microbial communities in response to elevated CO2. MBio. 2011;2(4):e00122–11. pmid:21791581
View Article
PubMed/NCBI
Google Scholar

[158] View Article

[159] PubMed/NCBI

[160] Google Scholar

[ref43] 43. Åkesson M, Förster J, Nielsen J. Integration of gene expression data into genome-scale metabolic models. Metabolic engineering. 2004;6(4):285–293. pmid:15491858
View Article
PubMed/NCBI
Google Scholar

[162] View Article

[163] PubMed/NCBI

[164] Google Scholar

[ref44] 44. Maarleveld TR, Boele J, Bruggeman FJ, Teusink B. A data integration and visualization resource for the metabolic network of Synechocystis sp. PCC 6803. Plant physiology. 2014; p. pp–113.

[ref45] 45. Yizhak K, Benyamini T, Liebermeister W, Ruppin E, Shlomi T. Integrating quantitative proteomics and metabolomics with a genome-scale metabolic network model. Bioinformatics. 2010;26(12):i255–i260. pmid:20529914
View Article
PubMed/NCBI
Google Scholar

[167] View Article

[168] PubMed/NCBI

[169] Google Scholar

[ref46] 46. Embree M, Liu JK, Al-Bassam MM, Zengler K. Networks of energetic and metabolic interactions define dynamics in microbial communities. Proceedings of the National Academy of Sciences. 2015;112(50):15450–15455. pmid:26621749
View Article
PubMed/NCBI
Google Scholar

[171] View Article

[172] PubMed/NCBI

[173] Google Scholar

[ref47] 47. Kraft B, Tegetmeyer HE, Sharma R, Klotz MG, Ferdelman TG, Hettich RL, et al. The environmental controls that govern the end product of bacterial nitrate respiration. Science. 2014;345(6197):676–679. pmid:25104387
View Article
PubMed/NCBI
Google Scholar

[175] View Article

[176] PubMed/NCBI

[177] Google Scholar

Figures

Abstract

Introduction

Materials and methods

Strains and cultivation conditions

Metabolite analysis

Quantification of cells

Models and software

Results

Design of co-culture and media composition that gives rise to different metabolic interactions

Experimental data of four growth conditions

Stoichiometric models of co-culture metabolism and growth

The metabolic models can simulate the experimental data

Influence of environmental conditions on the metabolic exchange fluxes between C. acetobutylicum and W. succinogenes

Intracellular fluxes of C. acetobutylicum can be inferred from community-level experimental data

Modelling predicts that co-culture behaviour can be enhanced

Discussion

Supporting information

S1 Material and Methods.

S1 Table. A coarse-grained metabolic model of W. succinogenes was detailed enough to fit the experimental data.

S2 Table. Parameters used for simulation of the co-culture.

S1 Fig. The dFBA simulations, with the draft genome-scale metabolic model of W. succinogenes, also agree mostly with the experimental data for the four different cultivation conditions.

S2 Fig. Product and substrate yields did not change during different growth rates.

S3 Fig. Inference of the metabolic interactions using a draft genome-scale metabolic model of W. succinogenes showed only differences in the NO2- production and NO3- consumption per utilized H2.

S1 File. Draft genome-scale metabolic model of W. succinogenes.

Acknowledgments

Author Contributions

References

S3 Fig. Inference of the metabolic interactions using a draft genome-scale metabolic model of W. succinogenes showed only differences in the NO₂^- production and NO₃^- consumption per utilized H₂.