Non-genetic adaptation by collective migration.

Collective behaviors require coordination of individuals. Thus, a population must adjust its phenotypic distribution to adapt to changing environments. How can a population regulate its phenotypic distribution? One strategy is to utilize specialized networks for gene regulation and maintaining distinct phenotypic subsets. Another involves genetic mutations, which can be augmented by stress-response pathways. Here, we studied how a migrating bacterial population regulates its phenotypic distribution to traverse across diverse environments. We generated isogenic Escherichia coli populations with varying distributions of swimming behaviors and observed their phenotype distributions during migration in liquid and porous environments. Surprisingly, we found that during collective migration, the distributions of swimming phenotypes adapt to the environment without mutations or gene regulation. Instead, adaptation is caused by the dynamic and reversible enrichment of high-performing swimming phenotypes within each environment. This adaptation mechanism is supported by a recent theoretical study, which proposed that the phenotypic composition of a migrating population results from a balance between cell growth generating diversity and collective migration eliminating the phenotypes that are unable to keep up with the migrating group. Furthermore, by examining chemoreceptor abundance distributions during migration towards different attractants, we found that this mechanism acts on multiple chemotaxis-related traits simultaneously. Our findings reveal that collective migration itself can enable cell populations with continuous, multi-dimensional phenotypes to flexibly and rapidly adapt their phenotypic composition to diverse environmental conditions. Significance statement: Conventional cell adaptation mechanisms, like gene regulation and random phenotypic switching, act swiftly but are limited to a few traits, while mutation-driven adaptations unfold slowly. By quantifying phenotypic diversity during bacterial collective migration, we discovered an adaptation mechanism that rapidly and reversibly adjusts multiple traits simultaneously. By dynamically balancing the elimination of phenotypes unable to keep pace with generation of diversity through growth, this process enables populations to tune their phenotypic composition based on the environment, without the need for gene regulation or mutations. Given the prevalence of collective migration in microbes, cancers, and embryonic development, non-genetic adaptation through collective migration may be a universal mechanism for populations to navigate diverse environments, offering insights into broader applications across various fields.

Direct measurement of dynamic attractant gradients reveals breakdown of the Patlak-Keller-Segel chemotaxis model.

Chemotactic bacteria not only navigate chemical gradients, but also shape their environments by consuming and secreting attractants. Investigating how these processes influence the dynamics of bacterial populations has been challenging because of a lack of experimental methods for measuring spatial profiles of chemoattractants in real time. Here, we use a fluorescent sensor for aspartate to directly measure bacterially generated chemoattractant gradients during collective migration. Our measurements show that the standard Patlak-Keller-Segel model for collective chemotactic bacterial migration breaks down at high cell densities. To address this, we propose modifications to the model that consider the impact of cell density on bacterial chemotaxis and attractant consumption. With these changes, the model explains our experimental data across all cell densities, offering insight into chemotactic dynamics. Our findings highlight the significance of considering cell density effects on bacterial behavior, and the potential for fluorescent metabolite sensors to shed light on the complex emergent dynamics of bacterial communities.

Direct measurement of dynamic attractant gradients reveals breakdown of the Patlak-Keller-Segel chemotaxis model.

Chemotactic bacteria not only navigate chemical gradients, but also shape their environments by consuming and secreting attractants. Investigating how these processes influence the dynamics of bacterial populations has been challenging because of a lack of experimental methods for measuring spatial profiles of chemoattractants in real time. Here, we use a fluorescent sensor for aspartate to directly measure bacterially generated chemoattractant gradients during collective migration. Our measurements show that the standard Patlak-Keller-Segel model for collective chemotactic bacterial migration breaks down at high cell densities. To address this, we propose modifications to the model that consider the impact of cell density on bacterial chemotaxis and attractant consumption. With these changes, the model explains our experimental data across all cell densities, offering new insight into chemotactic dynamics. Our findings highlight the significance of considering cell density effects on bacterial behavior, and the potential for fluorescent metabolite sensors to shed light on the complex emergent dynamics of bacterial communities.

Collective behavior and nongenetic inheritance allow bacterial populations to adapt to changing environments.

Collective behaviors require coordination among a group of individuals. As a result, individuals that are too phenotypically different from the rest of the group can be left out, reducing heterogeneity, but increasing coordination. If individuals also reproduce, the offspring can have different phenotypes from their parent(s). This raises the question of how these two opposing processes-loss of diversity by collective behaviors and generation of it through growth and inheritance-dynamically shape the phenotypic composition of an isogenic population. We examine this question theoretically using collective migration of chemotactic bacteria as a model system, where cells of different swimming phenotypes are better suited to navigate in different environments. We find that the differential loss of phenotypes caused by collective migration is environment-dependent. With cell growth, this differential loss enables migrating populations to dynamically adapt their phenotype compositions to the environment, enhancing migration through multiple environments. Which phenotypes are produced upon cell division depends on the level of nongenetic inheritance, and higher inheritance leads to larger composition adaptation and faster migration at steady state. However, this comes at the cost of slower responses to new environments. Due to this trade-off, there is an optimal level of inheritance that maximizes migration speed through changing environments, which enables a diverse population to outperform a nondiverse one. Growing populations might generally leverage the selection-like effects provided by collective behaviors to dynamically shape their own phenotype compositions, without mutations.

Activation-induced substrate engagement in ERK signaling.

The extracellular signal-regulated kinase (ERK) pathway is an essential component of developmental signaling in metazoans. Previous models of pathway activation suggested that dissociation of activated dually phosphorylated ERK (dpERK) from MAPK/ERK kinase (MEK), a kinase that phosphorylates ERK, and other cytoplasmic anchors, is sufficient for allowing ERK interactions with its substrates. Here, we provide evidence for an additional step controlling ERK's access to substrates. Specifically, we demonstrate that interaction of ERK with its substrate Capicua (Cic) is controlled at the level of ERK phosphorylation, whereby Cic binds to dpERK much stronger than to unphosphorylated ERK, both in vitro and in vivo. Mathematical modeling suggests that the differential affinity of Cic for dpERK versus ERK is required for both down-regulation of Cic and stabilizing phosphorylated ERK. Preferential association of Cic with dpERK serves two functions: it prevents unproductive competition of Cic with unphosphorylated ERK and contributes to efficient signal propagation. We propose that high-affinity substrate binding increases the specificity and efficiency of signal transduction through the ERK pathway.

Maximizing the information learned from finite data selects a simple model.

We use the language of uninformative Bayesian prior choice to study the selection of appropriately simple effective models. We advocate for the prior which maximizes the mutual information between parameters and predictions, learning as much as possible from limited data. When many parameters are poorly constrained by the available data, we find that this prior puts weight only on boundaries of the parameter space. Thus, it selects a lower-dimensional effective theory in a principled way, ignoring irrelevant parameter directions. In the limit where there are sufficient data to tightly constrain any number of parameters, this reduces to the Jeffreys prior. However, we argue that this limit is pathological when applied to the hyperribbon parameter manifolds generic in science, because it leads to dramatic dependence on effects invisible to experiment.

The Design Space of the Embryonic Cell Cycle Oscillator.

One of the main tasks in the analysis of models of biomolecular networks is to characterize the domain of the parameter space that corresponds to a specific behavior. Given the large number of parameters in most models, this is no trivial task. We use a model of the embryonic cell cycle to illustrate the approaches that can be used to characterize the domain of parameter space corresponding to limit cycle oscillations, a regime that coordinates periodic entry into and exit from mitosis. Our approach relies on geometric construction of bifurcation sets, numerical continuation, and random sampling of parameters. We delineate the multidimensional oscillatory domain and use it to quantify the robustness of periodic trajectories. Although some of our techniques explore the specific features of the chosen system, the general approach can be extended to other models of the cell cycle engine and other biomolecular networks.

Parallel imaging of Drosophila embryos for quantitative analysis of genetic perturbations of the Ras pathway.

The Ras pathway patterns the poles of the Drosophila embryo by downregulating the levels and activity of a DNA-binding transcriptional repressor Capicua (Cic). We demonstrate that the spatiotemporal pattern of Cic during this signaling event can be harnessed for functional studies of mutations in the Ras pathway in human diseases. Our approach relies on a new microfluidic device that enables parallel imaging of Cic dynamics in dozens of live embryos. We found that although the pattern of Cic in early embryos is complex, it can be accurately approximated by a product of one spatial profile and one time-dependent amplitude. Analysis of these functions of space and time alone reveals the differential effects of mutations within the Ras pathway. Given the highly conserved nature of Ras-dependent control of Cic, our approach provides new opportunities for functional analysis of multiple sequence variants from developmental abnormalities and cancers.

Long-term dynamics of multisite phosphorylation.

Multisite phosphorylation cycles are ubiquitous in cell regulation systems and are studied at multiple levels of complexity, from molecules to organisms, with the ultimate goal of establishing predictive understanding of the effects of genetic and pharmacological perturbations of protein phosphorylation in vivo. Achieving this goal is essentially impossible without mathematical models, which provide a systematic framework for exploring dynamic interactions of multiple network components. Most of the models studied to date do not discriminate between the distinct partially phosphorylated forms and focus on two limiting reaction regimes, distributive and processive, which differ in the number of enzyme-substrate binding events needed for complete phosphorylation or dephosphorylation. Here we use a minimal model of extracellular signal-related kinase regulation to explore the dynamics of a reaction network that includes all essential phosphorylation forms and arbitrary levels of reaction processivity. In addition to bistability, which has been studied extensively in distributive mechanisms, this network can generate periodic oscillations. Both bistability and oscillations can be realized at high levels of reaction processivity. Our work provides a general framework for systematic analysis of dynamics in multisite phosphorylation systems.

A Transport Model for Estimating the Time Course of ERK Activation in the C. elegans Germline.

The Caenorhabditis elegans germline is a well-studied model system for investigating the control of cell fate by signaling pathways. Cell signals at the distal tip of the germline promote cell proliferation; just before the loop, signals couple cell maturation to organism-level nutrient status; at the proximal end of the germline, signals coordinate oocyte maturation and fertilization in the presence of sperm. The latter two events require dual phosphorylation and activation of ERK, the effector molecule of the Ras/MAPK cascade. In C. elegans, ERK is known as MPK-1. At this point, none of today's methods for real-time monitoring of dually phosphorylated MPK-1 are working in the germline. Consequently, quantitative understanding of the MPK-1-dependent processes during germline development is limited. Here, we make a step toward advancing this understanding using a model-based framework that reconstructs the time course of MPK-1 activation from a snapshot of a fixed germline. Our approach builds on a number of recent studies for estimating temporal dynamics from fixed organisms, but takes advantage of the anatomy of the germline to simplify the analysis. Our model predicts that the MPK-1 signal turns on ∼30 h into germ cell progression and peaks ∼7 h later.