A mathematical model of syntrophic cocultures in the chemostat

MI-Sim: A MATLAB package for the numerical analysis of microbial ecological interactions

Food-webs and other classes of ecological network motifs, are a means of describing feeding relationships between consumers and producers in an ecosystem. They have application across scales where they differ only in the underlying characteristics of the organisms and substrates describing the system. Mathematical modelling, using mechanistic approaches to describe the dynamic behaviour and properties of the system through sets of ordinary differential equations, has been used extensively in ecology. Models allow simulation of the dynamics of the various motifs and their numerical analysis provides a greater understanding of the interplay between the system components and their intrinsic properties. We have developed the MI-Sim software for use with MATLAB to allow a rigorous and rapid numerical analysis of several common ecological motifs. MI-Sim contains a series of the most commonly used motifs such as cooperation, competition and predation. It does not require detailed knowledge of mathematical analytical techniques and is offered as a single graphical user interface containing all input and output options. The tools available in the current version of MI-Sim include model simulation, steady-state existence and stability analysis, and basin of attraction analysis. The software includes seven ecological interaction motifs and seven growth function models. Unlike other system analysis tools, MI-Sim is designed as a simple and user-friendly tool specific to ecological population type models, allowing for rapid assessment of their dynamical and behavioural properties.

Flexibility of syntrophic enzyme systems in Desulfovibrio species ensures their adaptation capability to environmental changes

The mineralization of organic matter in anoxic environments relies on the cooperative activities of hydrogen producers and consumers obligately linked by interspecies metabolite exchange in syntrophic consortia that may include sulfate reducing species such as Desulfovibrio. To evaluate the metabolic flexibility of syntrophic Desulfovibrio to adapt to naturally fluctuating methanogenic environments, we studied Desulfovibrio alaskensis strain G20 grown in chemostats under respiratory and syntrophic conditions with alternative methanogenic partners, Methanococcus maripaludis and Methanospirillum hungatei, at different growth rates. Comparative whole-genome transcriptional analyses, complemented by G20 mutant strain growth experiments and physiological data, revealed a significant influence of both energy source availability (as controlled by dilution rate) and methanogen on the electron transfer systems, ratios of interspecies electron carriers, energy generating systems, and interspecies physical associations. A total of 68 genes were commonly differentially expressed under syntrophic versus respiratory lifestyle. Under low-energy (low-growth-rate) conditions, strain G20 further had the capacity to adapt to the metabolism of its methanogenic partners, as shown by its differing gene expression of enzymes involved in the direct metabolic interactions (e.g., periplasmic hydrogenases) and the ratio shift in electron carriers used for interspecies metabolite exchange (hydrogen/formate). A putative monomeric [Fe-Fe] hydrogenase and Hmc (high-molecular-weight-cytochrome c3) complex-linked reverse menaquinone (MQ) redox loop become increasingly important for the reoxidation of the lactate-/pyruvate oxidation-derived redox pair, DsrC(red) and Fd(red), relative to the Qmo-MQ-Qrc (quinone-interacting membrane-bound oxidoreductase; quinone-reducing complex) loop. Together, these data underscore the high enzymatic and metabolic adaptive flexibility that likely sustains Desulfovibrio in naturally fluctuating methanogenic environments.

The mathematical analysis of a syntrophic relationship between two microbial species in a chemostat

A mathematical model involving a syntrophic relationship between two populations of bacteria in a continuous culture is proposed. A detailed qualitative analysis is carried out as well as the analysis of the local and global stability of the equilibria. We demonstrate, under general assumptions of monotonicity which are relevant from an applied point of view, the asymptotic stability of the positive equilibrium point which corresponds to the coexistence of the two bacteria. A syntrophic relationship in the anaerobic digestion process is proposed as a real candidate for this model.

Relative importance of trophic group concentrations during anaerobic degradation of volatile fatty acids

Although obligate syntrophic reactions cannot proceed without hydrogenotrophs, it has been unclear from the literature whether potential improvements are achievable with higher concentrations of hydrogenotrophs. In this study, the relative importance of formate-/H(2)-utilizing and acetate-utilizing trophic groups in the anaerobic degradation of butyrate and propionate was assessed by adding various proportions of these enriched cultures to a mixed anaerobic seed inoculum. The improvement resulting from the additional acetate-utilizing cultures was much greater than with formate/H(2) utilizers. Furthermore, formate/H(2) utilizers did not improve propionate utilization significantly, suggesting the importance of optimum utilization of hydrogenotrophic capacity. During most of the volatile fatty acid (VFA) degradation period, the system responded with characteristic hydrogen levels to maintain the Gibbs free energy of oxidation approximately constant for both butyrate (-6 kJ) and propionate (-14 kJ). These free-energy values were independent of methanogenic activity, as well as the volume of the seed inoculum and the VFA concentrations present. By comparing the experimental results with kinetic and mass transfer models, it was postulated that the diffusional transfer of reducing equivalents was the major limiting factor for efficient VFA degradation. Therefore, for optimum utilization of the hydrogenotrophs, low acetate concentrations are vital to enable the system to respond with higher formate/H(2) levels, thus leading to improved transfer of reducing equivalents. Due to the small number of propionate utilizers (and hence their limited surface area) and low bulk liquid concentrations, the additional formate/H(2) utilizers were of minimal use for improving the degradation rate further. The butyrate degradation rates strongly correlated with the cumulative activity of hydrogenotrophs and acetotrophs over the experimental range studied, indicating the need to model obligate syntrophic reactions as a dependent function of methanogenic activity.