The increased complexity of synthetic microbial biocircuits highlights the need for distributed cell functionality due to concomitant increases in metabolic and regulatory burdens imposed on single-strain topologies. Distributed systems, however, introduce additional challenges since consortium composition and spatiotemporal dynamics of constituent strains must be robustly controlled to achieve desired circuit behaviors. Here, we address these challenges with a modeling-based investigation of emergent spatiotemporal population dynamics using cell-length control in monolayer, two-strain bacterial consortia. We demonstrate that with dynamic control of a strain's division length, nematic cell alignment in close-packed monolayers can be destabilized. We find that this destabilization confers an emergent, competitive advantage to smaller-length strains-but by mechanisms that differ depending on the spatial patterns of the population. We used complementary modeling approaches to elucidate underlying mechanisms: an agent-based model to simulate detailed mechanical and signaling interactions between the competing strains, and a reductive, stochastic lattice model to represent cell-cell interactions with a single rotational parameter. Our modeling suggests that spatial strain-fraction oscillations can be generated when cell-length control is coupled to quorum-sensing signaling in negative feedback topologies. Our research employs novel methods of population control and points the way to programming strain fraction dynamics in consortial synthetic biology.

MOTIVATION: Advances in experimental and imaging techniques have allowed for unprecedented insights into the dynamical processes within individual cells. However, many facets of intracellular dynamics remain hidden, or can be measured only indirectly. This makes it challenging to reconstruct the regulatory networks that govern the biochemical processes underlying various cell functions. Current estimation techniques for inferring reaction rates frequently rely on marginalization over unobserved processes and states. Even in simple systems this approach can be computationally challenging, and can lead to large uncertainties and lack of robustness in parameter estimates. Therefore we will require alternative approaches to efficiently uncover the interactions in complex biochemical networks. RESULTS: We propose a Bayesian inference framework based on replacing uninteresting or unobserved reactions with time delays. Although the resulting models are non-Markovian, recent results on stochastic systems with random delays allow us to rigorously obtain expressions for the likelihoods of model parameters. In turn, this allows us to extend MCMC methods to efficiently estimate reaction rates, and delay distribution parameters, from single-cell assays. We illustrate the advantages, and potential pitfalls, of the approach using a birth-death model with both synthetic and experimental data, and show that we can robustly infer model parameters using a relatively small number of measurements. We demonstrate how to do so even when only the relative molecule count within the cell is measured, as in the case of fluorescence microscopy. AVAILABILITY AND IMPLEMENTATION: Accompanying code in R is available at https://github.com/cbskust/DDE_BD. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Modulation of gene network activity allows cells to respond to changes in environmental conditions. For example, the galactose utilization network in Saccharomyces cerevisiae is activated by the presence of galactose but repressed by glucose. If both sugars are present, the yeast will first metabolize glucose, depleting it from the extracellular environment. Upon depletion of glucose, the genes encoding galactose metabolic proteins will activate. Here, we show that the rate at which glucose levels are depleted determines the timing and variability of galactose gene activation. Paradoxically, we find that Gal1p, an enzyme needed for galactose metabolism, accumulates more quickly if glucose is depleted slowly rather than taken away quickly. Furthermore, the variability of induction times in individual cells depends non-monotonically on the rate of glucose depletion and exhibits a minimum at intermediate depletion rates. Our mathematical modeling suggests that the dynamics of the metabolic transition from glucose to galactose are responsible for the variability in galactose gene activation. These findings demonstrate that environmental dynamics can determine the phenotypic outcome at both the single-cell and population levels.

The creation of protein from DNA is a dynamic process consisting of numerous reactions, such as transcription, translation and protein folding. Each of these reactions is further comprised of hundreds or thousands of sub-steps that must be completed before a protein is fully mature. Consequently, the time it takes to create a single protein depends on the number of steps in the reaction chain and the nature of each step. One way to account for these reactions in models of gene regulatory networks is to incorporate dynamical delay. However, the stochastic nature of the reactions necessary to produce protein leads to a waiting time that is randomly distributed. Here, we use queueing theory to examine the effects of such distributed delay on the propagation of information through transcriptionally regulated genetic networks. In an analytically tractable model we find that increasing the randomness in protein production delay can increase signaling speed in transcriptional networks. The effect is confirmed in stochastic simulations, and we demonstrate its impact in several common transcriptional motifs. In particular, we show that in feedforward loops signaling time and magnitude are significantly affected by distributed delay. In addition, delay has previously been shown to cause stable oscillations in circuits with negative feedback. We show that the period and the amplitude of the oscillations monotonically decrease as the variability of the delay time increases.