Search

ABSTRACT The physiology and ecology of particle-associated marine bacteria are of growing interest, but our knowledge of their aggregation behavior and mechanisms controlling their association with particles remains limited. We have found that a particle-associated isolate, Alteromonas sp. ALT199 strain 4B03, and the related type-strain A. macleodii 27126 both form large (>500 μm) aggregates while growing in rich medium. A non-clumping variant (NCV) of 4B03 spontaneously arose in the lab, and whole-genome sequencing revealed a partial deletion in the gene encoding UDP-glucose-4-epimerase (galEΔ308–324). In 27126, a knock-out of galE (ΔgalE::kmr) resulted in a loss of aggregation, mimicking the NCV. Microscopic analysis shows that both 4B03 and 27126 rapidly form large aggregates, whereas their respective galE mutants remain primarily as single planktonic cells or clusters of a few cells. Strains 4B03 and 27126 also form aggregates with chitin particles, but their galE mutants do not. Alcian Blue staining shows that 4B03 and 27126 produce large transparent exopolymer particles (TEP), but their galE mutants are deficient in this regard. This study demonstrates the capabilities of cell-cell aggregation, aggregation of chitin particles, and production of TEP in strains of Alteromonas, a widespread particle-associated genus of heterotrophic marine bacteria. A genetic requirement for galE is evident for each of the above capabilities, expanding the known breadth of requirement for this gene in biofilm-related processes.IMPORTANCEHeterotrophic marine bacteria have a central role in the global carbon cycle. Well-known for releasing CO2 by decomposition and respiration, they may also contribute to particulate organic matter (POM) aggregation, which can promote CO2 sequestration via the formation of marine snow. We find that two members of the prevalent particle-associated genus Alteromonas can form aggregates comprising cells alone or cells and chitin particles, indicating their ability to drive POM aggregation. In line with their multivalent aggregation capability, both strains produce TEP, an excreted polysaccharide central to POM aggregation in the ocean. We demonstrate a genetic requirement for galE in aggregation and large TEP formation, building our mechanistic understanding of these aggregative capabilities. These findings point toward a role for heterotrophic bacteria in POM aggregation in the ocean and support broader efforts to understand bacterial controls on the global carbon cycle based on microbial activities, community structure, and meta-omic profiling.

Abstract Quantifying the contribution of individual molecular components to complex cellular processes is a grand challenge in systems biology. Here we establish a general theoretical framework (Functional Decomposition of Metabolism, FDM) to quantify the contribution of every metabolic reaction to metabolic functions, e.g. the synthesis of biomass building blocks. FDM allowed for a detailed quantification of the energy and biosynthesis budget for growing Escherichia coli cells. Surprisingly, the ATP generated during the biosynthesis of building blocks from glucose almost balances the demand from protein synthesis, the largest energy expenditure known for growing cells. This leaves the bulk of the energy generated by fermentation and respiration unaccounted for, thus challenging the common notion that energy is a key growth-limiting resource. Moreover, FDM together with proteomics enables the quantification of enzymes contributing towards each metabolic function, allowing for a first-principle formulation of a coarse-grained model of global protein allocation based on the structure of the metabolic network.

Abstract Metabolic cross-feeding plays vital roles in promoting ecological diversity. While some microbes depend on exchanges of essential nutrients for growth, the forces driving the extensive cross-feeding needed to support the coexistence of free-living microbes are poorly understood. Here we characterize bacterial physiology under self-acidification and establish that extensive excretion of key metabolites following growth arrest provides a collaborative, inter-species mechanism of stress resistance. This collaboration occurs not only between species isolated from the same community, but also between unrelated species with complementary (glycolytic vs. gluconeogenic) modes of metabolism. Cultures of such communities progress through distinct phases of growth-dilution cycles, comprising of exponential growth, acidification-triggered growth arrest, collaborative deacidification, and growth recovery, with each phase involving different combinations of physiological states of individual species. Our findings challenge the steady-state view of ecosystems commonly portrayed in ecological models, offering an alternative dynamical view based on growth advantages of complementary species in different phases.

ABSTRACT While Vibrio splendidus is best known as an opportunistic pathogen in oysters, Vibrio splendidus strain 1A01 was first identified as an early colonizer of synthetic chitin particles incubated in seawater. To gain a better understanding of its metabolism, a genome-scale metabolic model (GSMM) of V. splendidus 1A01 was reconstructed. GSMMs enable us to simulate all metabolic reactions in a bacterial cell using flux balance analysis. A draft model was built using an automated pipeline from BioCyc. Manual curation was then performed based on experimental data, in part by gap-filling metabolic pathways and tailoring the model’s biomass reaction to V. splendidus 1A01. The challenges of building a metabolic model for a marine microorganism like V. splendidus 1A01 are described. IMPORTANCE A genome-scale metabolic model of V. splendidus 1A01 was reconstructed in this work. We offer solutions to the technical problems associated with model reconstruction for a marine bacterial strain like V. splendidus 1A01, which arise largely from the high salt concentration found in both seawater and culture media that simulate seawater.

Micro-scale microbial community dynamics can substantially alter the fate of sinking particulates in the ocean thus playing a key role in setting the vertical flux of particulate carbon in the ocean.

Abstract Accurate measurements of cellular protein concentrations are invaluable to quantitative studies of gene expression and physiology in living cells. Here, we developed a versatile mass spectrometric workflow based on data‐independent acquisition proteomics (DIA/SWATH) together with a novel protein inference algorithm (xTop). We used this workflow to accurately quantify absolute protein abundances in Escherichiacoli for > 2,000 proteins over > 60 growth conditions, including nutrient limitations, non‐metabolic stresses, and non‐planktonic states. The resulting high‐quality dataset of protein mass fractions allowed us to characterize proteome responses from a coarse (groups of related proteins) to a fine (individual) protein level. Hereby, a plethora of novel biological findings could be elucidated, including the generic upregulation of low‐abundant proteins under various metabolic limitations, the non‐specificity of catabolic enzymes upregulated under carbon limitation, the lack of large‐scale proteome reallocation under stress compared to nutrient limitations, as well as surprising strain‐dependent effects important for biofilm formation. These results present valuable resources for the systems biology community and can be used for future multi‐omics studies of gene regulation and metabolic control in E. coli.

Abstract To respond to fluctuating conditions, microbes typically need to synthesize novel proteins. As this synthesis relies on sufficient biosynthetic precursors, microbes must devise effective response strategies to manage depleting precursors. To better understand these strategies, we investigate the active response of Escherichia coli to changes in nutrient conditions, connecting transient gene expression to growth phenotypes. By synthetically modifying gene expression during changing conditions, we show how the competition by genes for the limited protein synthesis capacity constrains cellular response. Despite this constraint cells substantially express genes that are not required, trapping them in states where precursor levels are low and the genes needed to replenish the precursors are outcompeted. Contrary to common modeling assumptions, our findings highlight that cells do not optimize growth under changing environments but rather exhibit hardwired response strategies that may have evolved to promote fitness in their native environment. The constraint and the suboptimality of the cellular response uncovered provide a conceptual framework relevant for many research applications, from the prediction of evolution to the improvement of gene circuits in biotechnology.

Much recent progress has been made to understand the impact of proteome allocation on bacterial growth; much less is known about the relationship between the abundances of the enzymes and their substrates, which jointly determine metabolic fluxes. Here, we report a correlation between the concentrations of enzymes and their substrates in Escherichia coli. We suggest this relationship to be a consequence of optimal resource allocation, subject to an overall constraint on the biomass density: For a cellular reaction network composed of effectively irreversible reactions, maximal reaction flux is achieved when the dry mass allocated to each substrate is equal to the dry mass of the unsaturated (or “free”) enzymes waiting to consume it. Calculations based on this optimality principle successfully predict the quantitative relationship between the observed enzyme and metabolite abundances, parameterized only by molecular masses and enzyme–substrate dissociation constants (Km). The corresponding organizing principle provides a fundamental rationale for cellular investment into different types of molecules, which may aid in the design of more efficient synthetic cellular systems. This study shows that in E. coli, the cellular mass of each metabolite approximately equals the combined mass of the free enzymes waiting to consume it; this simple relationship arises from the optimal utilization of cellular dry mass, and quantitatively describes available experimental data.

Abstract Microorganisms adjust metabolic activity to cope with diverse environments. While many studies have provided insights into how individual pathways are regulated, the mechanisms that give rise to coordinated metabolic responses are poorly understood. Here, we identify the regulatory mechanisms that coordinate catabolism and anabolism in Escherichia coli. Integrating protein, metabolite, and flux changes in genetically implemented catabolic or anabolic limitations, we show that combined global and local mechanisms coordinate the response to metabolic limitations. To allocate proteomic resources between catabolism and anabolism, E. coli uses a simple global gene regulatory program. Surprisingly, this program is largely implemented by a single transcription factor, Crp, which directly activates the expression of catabolic enzymes and indirectly reduces the expression of anabolic enzymes by passively sequestering cellular resources needed for their synthesis. However, metabolic fluxes are not controlled by this regulatory program alone; instead, fluxes are adjusted mostly through passive changes in the local metabolite concentrations. These mechanisms constitute a simple but effective global regulatory program that coarsely partitions resources between different parts of metabolism while ensuring robust coordination of individual metabolic reactions.

Fast-growing bacteria produce many proteins in excess of what seems optimal for exponential growth. Here, the authors present a mathematical model and experimental evidence supporting that this overexpression serves as a strategic reserve to quickly meet demand upon sudden improvement in growth conditions.

The physical interactions of growing bacterial cells with each other and with their surroundings significantly affect the structure and dynamics of biofilms. Here a 3D agent-based model is formulated to describe the establishment of simple bacterial colonies expanding by the physical force of their growth. With a single set of parameters, the model captures key dynamical features of colony growth by non-motile, non EPS-producing E. coli cells on hard agar. The model, supported by experiment on colony growth in different types and concentrations of nutrients, suggests that radial colony expansion is not limited by nutrients as commonly believed, but by mechanical forces. Nutrient penetration instead governs vertical colony growth, through thin layers of vertically oriented cells lifting up their ancestors from the bottom. Overall, the model provides a versatile platform to investigate the influences of metabolic and environmental factors on the growth and morphology of bacterial colonies.

ABSTRACT In nature, bacteria frequently experience many adverse conditions, including heat, oxidation, acidity, and hyperosmolarity, which all tend to slow down if not outright stop cell growth. Previous work on bacterial stress mainly focused on understanding gene regulatory responses. Much less is known about how stresses compromise protein synthesis, which is the major driver of cell growth. Here, we quantitatively characterize the translational capacity of Escherichia coli cells growing exponentially under hyperosmotic stress. We found that hyperosmotic stress affects bacterial protein synthesis through reduction of the translational elongation rate, which is largely compensated for by an increase in the cellular ribosome content compared with nutrient limitation at a similar growth rate. The slowdown of translational elongation is attributed to a reduction in the rate of binding of tRNA ternary complexes to the ribosomes. IMPORTANCE Hyperosmotic stress is a common stress condition confronted by E. coli during infection of the urinary tract. It can significantly compromise the bacterial growth rate. Protein translation capacity is a critical component of bacterial growth. In this study, we find for the first time that hyperosmotic stress causes substantial slowdown in bacterial ribosome translation elongation. The slowdown of translation elongation originates from a reduced binding rate of tRNA ternary complex to the ribosomes.

New experimental results on bacterial growth inspire a novel top-down approach to study cell metabolism, combining mass balance and proteomic constraints to extend and complement Flux Balance Analysis. We introduce here Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic costs associated to growth are accounted for in an effective way through a single additional genome-wide constraint. Its roots lie in the experimentally observed pattern of proteome allocation for metabolic functions, allowing to bridge regulation and metabolism in a transparent way under the principle of growth-rate maximization. We provide a simple method to solve CAFBA efficiently and propose an "ensemble averaging" procedure to account for unknown protein costs. Applying this approach to modeling E. coli metabolism, we find that, as the growth rate increases, CAFBA solutions cross over from respiratory, growth-yield maximizing states (preferred at slow growth) to fermentative states with carbon overflow (preferred at fast growth). In addition, CAFBA allows for quantitatively accurate predictions on the rate of acetate excretion and growth yield based on only 3 parameters determined by empirical growth laws.

Abstract Understanding how the homeostasis of cellular size and composition is accomplished by different organisms is an outstanding challenge in biology. For exponentially growing Escherichia coli cells, it is long known that the size of cells exhibits a strong positive relation with their growth rates in different nutrient conditions. Here, we characterized cell sizes in a set of orthogonal growth limitations. We report that cell size and mass exhibit positive or negative dependences with growth rate depending on the growth limitation applied. In particular, synthesizing large amounts of “useless” proteins led to an inversion of the canonical, positive relation, with slow growing cells enlarged 7‐ to 8‐fold compared to cells growing at similar rates under nutrient limitation. Strikingly, this increase in cell size was accompanied by a 3‐ to 4‐fold increase in cellular DNA content at slow growth, reaching up to an amount equivalent to ~8 chromosomes per cell. Despite drastic changes in cell mass and macromolecular composition, cellular dry mass density remained constant. Our findings reveal an important role of protein synthesis in cell division control.

Abstract Bacteria must constantly adapt their growth to changes in nutrient availability; yet despite large‐scale changes in protein expression associated with sensing, adaptation, and processing different environmental nutrients, simple growth laws connect the ribosome abundance and the growth rate. Here, we investigate the origin of these growth laws by analyzing the features of ribosomal regulation that coordinate proteome‐wide expression changes with cell growth in a variety of nutrient conditions in the model organism Escherichia coli. We identify supply‐driven feedforward activation of ribosomal protein synthesis as the key regulatory motif maximizing amino acid flux, and autonomously guiding a cell to achieve optimal growth in different environments. The growth laws emerge naturally from the robust regulatory strategy underlying growth rate control, irrespective of the details of the molecular implementation. The study highlights the interplay between phenomenological modeling and molecular mechanisms in uncovering fundamental operating constraints, with implications for endogenous and synthetic design of microorganisms.

Different codons encoding the same amino acid are not used equally in protein-coding sequences. In bacteria, there is a bias towards codons with high translation rates. This bias is most pronounced in highly expressed proteins, but a recent study of synthetic GFP-coding sequences did not find a correlation between codon usage and GFP expression, suggesting that such correlation in natural sequences is not a simple property of translational mechanisms. Here, we investigate the effect of evolutionary forces on codon usage. The relation between codon bias and protein abundance is quantitatively analyzed based on the hypothesis that codon bias evolved to ensure the efficient usage of ribosomes, a precious commodity for fast growing cells. An explicit fitness landscape is formulated based on bacterial growth laws to relate protein abundance and ribosomal load. The model leads to a quantitative relation between codon bias and protein abundance, which accounts for a substantial part of the observed bias for E. coli. Moreover, by providing an evolutionary link, the ribosome load model resolves the apparent conflict between the observed relation of protein abundance and codon bias in natural sequences and the lack of such dependence in a synthetic gfp library. Finally, we show that the relation between codon usage and protein abundance can be used to predict protein abundance from genomic sequence data alone without adjustable parameters.

Cells employ a myriad of signaling circuits to detect environmental signals and drive specific gene expression responses. A common motif in these circuits is inducible auto-activation: a transcription factor that activates its own transcription upon activation by a ligand or by post-transcriptional modification. Examples range from the two-component signaling systems in bacteria and plants to the genetic circuits of animal viruses such as HIV. We here present a theoretical study of such circuits, based on analytical calculations, numerical computations, and simulation. Our results reveal several surprising characteristics. They show that auto-activation can drastically enhance the sensitivity of the circuit's response to input signals: even without molecular cooperativity, an ultra-sensitive threshold response can be obtained. However, the increased sensitivity comes at a cost: auto-activation tends to severely slow down the speed of induction, a stochastic effect that was strongly underestimated by earlier deterministic models. This slow-induction effect again requires no molecular cooperativity and is intimately related to the bimodality recently observed in non-cooperative auto-activation circuits. These phenomena pose strong constraints on the use of auto-activation in signaling networks. To achieve both a high sensitivity and a rapid induction, an inducible auto-activation circuit is predicted to acquire low cooperativity and low fold-induction. Examples from Escherichia coli's two-component signaling systems support these predictions.

Predictive understanding of the myriads of signal transduction pathways in a cell is an outstanding challenge of systems biology. Such pathways are primarily mediated by specific but transient protein-protein interactions, which are difficult to study experimentally. In this study, we dissect the specificity of protein-protein interactions governing two-component signaling (TCS) systems ubiquitously used in bacteria. Exploiting the large number of sequenced bacterial genomes and an operon structure which packages many pairs of interacting TCS proteins together, we developed a computational approach to extract a molecular interaction code capturing the preferences of a small but critical number of directly interacting residue pairs. This code is found to reflect physical interaction mechanisms, with the strongest signal coming from charged amino acids. It is used to predict the specificity of TCS interaction: Our results compare favorably to most available experimental results, including the prediction of 7 (out of 8 known) interaction partners of orphan signaling proteins in Caulobacter crescentus. Surveying among the available bacterial genomes, our results suggest 15∼25% of the TCS proteins could participate in out-of-operon "crosstalks". Additionally, we predict clusters of crosstalking candidates, expanding from the anecdotally known examples in model organisms. The tools and results presented here can be used to guide experimental studies towards a system-level understanding of two-component signaling.

An increasing number of small RNAs (sRNAs) have been shown to regulate critical pathways in prokaryotes and eukaryotes. In bacteria, regulation by trans-encoded sRNAs is predominantly found in the coordination of intricate stress responses. The mechanisms by which sRNAs modulate expression of its targets are diverse. In common to most is the possibility that interference with the translation of mRNA targets may also alter the abundance of functional sRNAs. Aiming to understand the unique role played by sRNAs in gene regulation, we studied examples from two distinct classes of bacterial sRNAs in Escherichia coli using a quantitative approach combining experiment and theory. Our results demonstrate that sRNA provides a novel mode of gene regulation, with characteristics distinct from those of protein-mediated gene regulation. These include a threshold-linear response with a tunable threshold, a robust noise resistance characteristic, and a built-in capability for hierarchical cross-talk. Knowledge of these special features of sRNA-mediated regulation may be crucial toward understanding the subtle functions that sRNAs can play in coordinating various stress-relief pathways. Our results may also help guide the design of synthetic genetic circuits that have properties difficult to attain with protein regulators alone.