Fast decisions reflect biases; slow decisions do not.

Decisions are often made by heterogeneous groups of individuals, each with distinct initial biases and access to information of different quality. We show that in groups of independent agents who accumulate evidence the first to decide are those with the strongest initial biases. Their decisions align with their initial bias, regardless of the underlying truth. In contrast, agents who decide last make decisions as if they were initially unbiased and hence make better choices. We obtain asymptotic expressions in the large population limit quantifying how agents' initial inclinations shape early decisions. Our analysis shows how bias, information quality, and decision order interact in nontrivial ways to determine the reliability of decisions in a group.

Pattern Formation and Bistability in a Synthetic Intercellular Genetic Toggle.

Differentiation within multicellular organisms is a complex process that helps to establish spatial patterning and tissue formation within the body. Often, the differentiation of cells is governed by morphogens and intercellular signaling molecules that guide the fate of each cell, frequently using toggle-like regulatory components. Synthetic biologists have long sought to recapitulate patterned differentiation with engineered cellular communities, and various methods for differentiating bacteria have been invented. Here, we couple a synthetic corepressive toggle switch with intercellular signaling pathways to create a "quorum-sensing toggle". We show that this circuit not only exhibits population-wide bistability in a well-mixed liquid environment but also generates patterns of differentiation in colonies grown on agar containing an externally supplied morphogen. If coupled to other metabolic processes, circuits such as the one described here would allow for the engineering of spatially patterned, differentiated bacteria for use in biomaterials and bioelectronics.

Fast decisions reflect biases, slow decisions do not.

Decisions are often made by heterogeneous groups of individuals, each with distinct initial biases and access to information of different quality. We show that in large groups of independent agents who accumulate evidence the first to decide are those with the strongest initial biases. Their decisions align with their initial bias, regardless of the underlying truth. In contrast, agents who decide last make decisions as if they were initially unbiased, and hence make better choices. We obtain asymptotic expressions in the large population limit that quantify how agents' initial inclinations shape early decisions. Our analysis shows how bias, information quality, and decision order interact in non-trivial ways to determine the reliability of decisions in a group.

Inferring delays in partially observed gene regulation processes.

MOTIVATION: Cell function is regulated by gene regulatory networks (GRNs) defined by protein-mediated interaction between constituent genes. Despite advances in experimental techniques, we can still measure only a fraction of the processes that govern GRN dynamics. To infer the properties of GRNs using partial observation, unobserved sequential processes can be replaced with distributed time delays, yielding non-Markovian models. Inference methods based on the resulting model suffer from the curse of dimensionality. RESULTS: We develop a simulation-based Bayesian MCMC method employing an approximate likelihood for the efficient and accurate inference of GRN parameters when only some of their products are observed. We illustrate our approach using a two-step activation model: an activation signal leads to the accumulation of an unobserved regulatory protein, which triggers the expression of observed fluorescent proteins. With prior information about observed fluorescent protein synthesis, our method successfully infers the dynamics of the unobserved regulatory protein. We can estimate the delay and kinetic parameters characterizing target regulation including transcription, translation, and target searching of an unobserved protein from experimental measurements of the products of its target gene. Our method is scalable and can be used to analyze non-Markovian models with hidden components. AVAILABILITY AND IMPLEMENTATION: Our code is implemented in R and is freely available with a simple example data at https://github.com/Mathbiomed/SimMCMC.

Tunable Dynamics in a Multistrain Transcriptional Pulse Generator.

One challenge in synthetic biology is the tuning of regulatory components within gene circuits to elicit a specific behavior. This challenge becomes more difficult in synthetic microbial consortia since each strain's circuit must function at the intracellular level and their combination must operate at the population level. Here we demonstrate that circuit dynamics can be tuned in synthetic consortia through the manipulation of strain fractions within the community. To do this, we construct a microbial consortium comprised of three strains of engineered Escherichia coli that, when cocultured, use homoserine lactone-mediated intercellular signaling to create a multistrain incoherent type-1 feedforward loop (I1-FFL). Like naturally occurring I1-FFL motifs in gene networks, this engineered microbial consortium acts as a pulse generator of gene expression. We demonstrate that the amplitude of the pulse can be easily tuned by adjusting the relative population fractions of the strains. We also develop a mathematical model for the temporal dynamics of the microbial consortium. This model allows us to identify population fractions that produced desired pulse characteristics, predictions that were confirmed for all but extreme fractions. Our work demonstrates that intercellular gene circuits can be effectively tuned simply by adjusting the starting fractions of each strain in the consortium.

Sampling-based Bayesian inference in recurrent circuits of stochastic spiking neurons.

Two facts about cortex are widely accepted: neuronal responses show large spiking variability with near Poisson statistics and cortical circuits feature abundant recurrent connections between neurons. How these spiking and circuit properties combine to support sensory representation and information processing is not well understood. We build a theoretical framework showing that these two ubiquitous features of cortex combine to produce optimal sampling-based Bayesian inference. Recurrent connections store an internal model of the external world, and Poissonian variability of spike responses drives flexible sampling from the posterior stimulus distributions obtained by combining feedforward and recurrent neuronal inputs. We illustrate how this framework for sampling-based inference can be used by cortex to represent latent multivariate stimuli organized either hierarchically or in parallel. A neural signature of such network sampling are internally generated differential correlations whose amplitude is determined by the prior stored in the circuit, which provides an experimentally testable prediction for our framework.

Rapid compensatory plasticity revealed by dynamic correlated activity in monkeys in vivo.

To produce adaptive behavior, neural networks must balance between plasticity and stability. Computational work has demonstrated that network stability requires plasticity mechanisms to be counterbalanced by rapid compensatory processes. However, such processes have yet to be experimentally observed. Here we demonstrate that repeated optogenetic activation of excitatory neurons in monkey visual cortex (area V1) induces a population-wide dynamic reduction in the strength of neuronal interactions over the timescale of minutes during the awake state, but not during rest. This new form of rapid plasticity was observed only in the correlation structure, with firing rates remaining stable across trials. A computational network model operating in the balanced regime confirmed experimental findings and revealed that inhibitory plasticity is responsible for the decrease in correlated activity in response to repeated light stimulation. These results provide the first experimental evidence for rapid homeostatic plasticity that primarily operates during wakefulness, which stabilizes neuronal interactions during strong network co-activation.

Hierarchical Modular Structure of the Drosophila Connectome.

The structure of neural circuitry plays a crucial role in brain function. Previous studies of brain organization generally had to trade off between coarse descriptions at a large scale and fine descriptions on a small scale. Researchers have now reconstructed tens to hundreds of thousands of neurons at synaptic resolution, enabling investigations into the interplay between global, modular organization, and cell type-specific wiring. Analyzing data of this scale, however, presents unique challenges. To address this problem, we applied novel community detection methods to analyze the synapse-level reconstruction of an adult female Drosophila melanogaster brain containing >20,000 neurons and 10 million synapses. Using a machine-learning algorithm, we find the most densely connected communities of neurons by maximizing a generalized modularity density measure. We resolve the community structure at a range of scales, from large (on the order of thousands of neurons) to small (on the order of tens of neurons). We find that the network is organized hierarchically, and larger-scale communities are composed of smaller-scale structures. Our methods identify well-known features of the fly brain, including its sensory pathways. Moreover, focusing on specific brain regions, we are able to identify subnetworks with distinct connectivity types. For example, manual efforts have identified layered structures in the fan-shaped body. Our methods not only automatically recover this layered structure, but also resolve finer connectivity patterns to downstream and upstream areas. We also find a novel modular organization of the superior neuropil, with distinct clusters of upstream and downstream brain regions dividing the neuropil into several pathways. These methods show that the fine-scale, local network reconstruction made possible by modern experimental methods are sufficiently detailed to identify the organization of the brain across scales, and enable novel predictions about the structure and function of its parts.Significance Statement The Hemibrain is a partial connectome of an adult female Drosophila melanogaster brain containing >20,000 neurons and 10 million synapses. Analyzing the structure of a network of this size requires novel and efficient computational tools. We applied a new community detection method to automatically uncover the modular structure in the Hemibrain dataset by maximizing a generalized modularity measure. This allowed us to resolve the community structure of the fly hemibrain at a range of spatial scales revealing a hierarchical organization of the network, where larger-scale modules are composed of smaller-scale structures. The method also allowed us to identify subnetworks with distinct cell and connectivity structures, such as the layered structures in the fan-shaped body, and the modular organization of the superior neuropil. Thus, network analysis methods can be adopted to the connectomes being reconstructed using modern experimental methods to reveal the organization of the brain across scales. This supports the view that such connectomes will allow us to uncover the organizational structure of the brain, which can ultimately lead to a better understanding of its function.

Signaling in microbial communities with open boundaries.

Microbial communities such as swarms or biofilms often form at the interfaces of solid substrates and open fluid flows. At the same time, in laboratory environments these communities are commonly studied using microfluidic devices with media flows and open boundaries. Extracellular signaling within these communities is therefore subject to different constraints than signaling within classic, closed-boundary systems such as developing embryos or tissues, yet is understudied by comparison. Here, we use mathematical modeling to show how advective-diffusive boundary flows and population geometry impact cell-cell signaling in monolayer microbial communities. We reveal conditions where the intercellular signaling lengthscale depends solely on the population geometry and not on diffusion or degradation, as commonly expected. We further demonstrate that diffusive coupling with the boundary flow can produce signal gradients within an isogenic population, even when there is no flow within the population. We use our theory to provide new insights into the signaling mechanisms of published experimental results, and we make several experimentally verifiable predictions. Our research highlights the importance of carefully evaluating boundary dynamics and environmental geometry when modeling microbial cell-cell signaling and informs the study of cell behaviors in both natural and synthetic systems.

Signaling in microbial communities with open boundaries.

Microbial communities such as swarms or biofilms often form at the interfaces of solid substrates and open fluid flows. At the same time, in laboratory environments these communities are commonly studied using microfluidic devices with media flows and open boundaries. Extracellular signaling within these communities is therefore subject to different constraints than signaling within classic, closed-boundary systems such as developing embryos or tissues, yet is understudied by comparison. Here, we use mathematical modeling to show how advective-diffusive boundary flows and population geometry impact cell-cell signaling in monolayer microbial communities. We reveal conditions where the intercellular signaling lengthscale depends solely on the population geometry and not on diffusion or degradation, as commonly expected. We further demonstrate that diffusive coupling with the boundary flow can produce signal gradients within an isogenic population, even when there is no flow within the population. We use our theory to provide new insights into the signaling mechanisms of published experimental results, and we make several experimentally verifiable predictions. Our research highlights the importance of carefully evaluating boundary dynamics and environmental geometry when modeling microbial cell-cell signaling and informs the study of cell behaviors in both natural and synthetic systems.