BackgroundEnteric Escherichia coli survives the highly acidic environment of the stomach through multiple acid resistance (AR) mechanisms. The most effective system, AR2, decarboxylates externally-derived glutamate to remove cytoplasmic protons and excrete GABA. The first described system, AR1, does not require an external amino acid. Its mechanism has not been determined. The regulation of the multiple AR systems and their coordination with broader cellular metabolism has not been fully explored.ResultsWe utilized a combination of ChIP-Seq and gene expression analysis to experimentally map the regulatory interactions of four TFs: nac, ntrC, ompR, and csiR. Our data identified all previously in vivo confirmed direct interactions and revealed several others previously inferred from gene expression data. Our data demonstrate that nac and csiR directly modulate AR, and leads to a regulatory network model in which all four TFs participate in coordinating acid resistance, glutamate metabolism, and nitrogen metabolism. This model predicts a novel mechanism for AR1 by which the decarboxylation enzymes of AR2 are used with internally derived glutamate. This hypothesis makes several testable predictions that we confirmed experimentally.ConclusionsOur data suggest that the regulatory network underlying AR is complex and deeply interconnected with the regulation of GABA and glutamate metabolism, nitrogen metabolism. These connections underlie and experimentally validated model of AR1 in which the decarboxylation enzymes of AR2 are used with internally derived glutamate.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-016-0376-y) contains supplementary material, which is available to authorized users.

Microbes face a trade-off between being metabolically independent and relying on neighboring organisms for the supply of some essential metabolites. This balance of conflicting strategies affects microbial community structure and dynamics, with important implications for microbiome research and synthetic ecology. A “gedanken” (thought) experiment to investigate this trade-off would involve monitoring the rise of mutual dependence as the number of metabolic reactions allowed in an organism is increasingly constrained. The expectation is that below a certain number of reactions, no individual organism would be able to grow in isolation and cross-feeding partnerships and division of labor would emerge. We implemented this idealized experiment usingin silicogenome-scale models. In particular, we used mixed-integer linear programming to identify trade-off solutions in communities ofEscherichia colistrains. The strategies that we found revealed a large space of opportunities in nuanced and nonintuitive metabolic division of labor, including, for example, splitting the tricarboxylic acid (TCA) cycle into two separate halves. The systematic computation of possible solutions in division of labor for 1-, 2-, and 3-strain consortia resulted in a rich and complex landscape. This landscape displayed a nonlinear boundary, indicating that the loss of an intracellular reaction was not necessarily compensated for by a single imported metabolite. Different regions in this landscape were associated with specific solutions and patterns of exchanged metabolites. Our approach also predicts the existence of regions in this landscape where independent bacteria are viable but are outcompeted by cross-feeding pairs, providing a possible incentive for the rise of division of labor.IMPORTANCEUnderstanding how microbes assemble into communities is a fundamental open issue in biology, relevant to human health, metabolic engineering, and environmental sustainability. A possible mechanism for interactions of microbes is through cross-feeding, i.e., the exchange of small molecules. These metabolic exchanges may allow different microbes to specialize in distinct tasks and evolve division of labor. To systematically explore the space of possible strategies for division of labor, we applied advanced optimization algorithms to computational models of cellular metabolism. Specifically, we searched for communities able to survive under constraints (such as a limited number of reactions) that would not be sustainable by individual species. We found that predicted consortia partition metabolic pathways in ways that would be difficult to identify manually, possibly providing a competitive advantage over individual organisms. In addition to helping understand diversity in natural microbial communities, our approach could assist in the design of synthetic consortia.

Microbes face a tradeoff between being metabolically independent and relying on neighboring organisms for the supply of some essential metabolites. This balance of conflicting strategies affects microbial community structure and dynamics, with important implications for microbiome research and synthetic ecology. A "gedanken experiment" to investigate this tradeoff would involve monitoring the rise of mutual dependence as the number of metabolic reactions allowed in an organism is increasingly constrained. The expectation is that below a certain number of reactions, no individual organism would be able to grow in isolation, and cross-feeding partnerships and division of labor would emerge. We implemented this idealized experiment using in silico genome-scale models. In particular, we used mixed integer linear programming to identify tradeoff solutions in communities of Escherichia coli strains. The strategies we found reveal a large space of nuanced and nonintuitive metabolic division of labor opportunities, including, for example, splitting the TCA cycle into two separate halves. The systematic computation of possible division of labor solutions for 1-, 2-, and 3-strain consortia resulted in a rich and complex landscape. This landscape displays a nonlinear boundary, indicating that the loss of an intracellular reaction is not necessarily compensated by a single imported metabolite. Different regions in this landscape are associated with specific solutions and patterns of exchanged metabolites. Our approach also predicts the existence of regions in this landscape where independent bacteria are viable, but outcompeted by cross-feeding pairs, providing a possible incentive for the rise of division of labor.

Given an array of phenotypes (e.g., yield across strains and conditions), one can ask how to best choose subsets of conditions that are informative about the whole dataset, enabling efficient system identification and providing a basis vector in phenotype space. Here we introduce a mixed integer linear programming approach to choose explanatory and response variables for a phenotypic matrix. We applied the algorithm to a set of fitness measurements for 462 yeast strains under 38 carbon sources, and to the growth phenotypes of 65 marine bacteria on 11 media. The algorithm identifies environments that can be used as features to predict growth under other conditions, providing biologically interpretable metabolic axes for strain discrimination. Our approach could be used to reduce the number of experiments needed to identify a strain or to map its metabolic capabilities. The generality of the algorithm makes it appropriate for addressing subset selection problems in areas beyond biology.