Many cellular patterns exhibit a reaction-diffusion component, suggesting that Turing instability may contribute to pattern formation. However, biological gene-regulatory pathways are more complex than simple Turing activator-inhibitor models and generally do not require fine-tuning of parameters as dictated by the Turing conditions. To address these issues, we employ random matrix theory to analyze the Jacobian matrices of larger networks with robust statistical properties. Our analysis reveals that Turing patterns are more likely to occur by chance than previously thought and that the most robust Turing networks have an optimal size, consisting of only a handful of molecular species, thus significantly increasing their identifiability in biological systems. Broadly speaking, this optimal size emerges from a trade-off between the highest stability in small networks and the greatest instability with diffusion in large networks. Furthermore, we find that with multiple immobile nodes, differential diffusion ceases to be important for Turing patterns. Our findings may inform future synthetic biology approaches and provide insights into bridging the gap to complex developmental pathways.

Single-cell organisms and various cell types use a range of motility modes when following a chemical gradient, but it is unclear which mode is best suited for different gradients. Here, we model directional decision-making in chemotactic amoeboid cells as a stimulus-dependent actin recruitment contest. Pseudopods extending from the cell body compete for a finite actin pool to push the cell in their direction until one pseudopod wins and determines the direction of movement. Our minimal model provides a quantitative understanding of the strategies cells use to reach the physical limit of accurate chemotaxis, aligning with data without explicit gradient sensing or cellular memory for persistence. To generalize our model, we employ reinforcement learning optimization to study the effect of pseudopod suppression, a simple but effective cellular algorithm by which cells can suppress possible directions of movement. Different pseudopod-based chemotaxis strategies emerge naturally depending on the environment and its dynamics. For instance, in static gradients, cells can react faster at the cost of pseudopod accuracy, which is particularly useful in noisy, shallow gradients where it paradoxically increases chemotactic accuracy. In contrast, in dynamics gradients, cells form de novo pseudopods. Overall, our work demonstrates mechanical intelligence for high chemotaxis performance with minimal cellular regulation.

Spatial information from cell-surface receptors is crucial for processes that require signal processing and sensing of the environment. Here, we investigate the optimal placement of such receptors through a theoretical model that minimizes uncertainty in gradient estimation. Without requiring a priori knowledge of the physical limits of sensing or biochemical processes, we reproduce the emergence of clusters that closely resemble those observed in real cells. On perfect spherical surfaces, optimally placed receptors spread uniformly. When perturbations break their symmetry, receptors cluster in regions of high curvature, massively reducing estimation uncertainty. This agrees in many scenarios with mechanistic models that minimize elastic preference discrepancies between receptors and cell membranes. We further extend our model to motile receptors responding to cell-shape changes and external fluid flow, demonstrating the biological relevance of our model. Our findings provide a simple and utilitarian explanation for receptor clustering at high-curvature regions when high sensing accuracy is paramount.

SummaryTuring patterns1are well-known self-organising systems that can form spots, stripes, or labyrinths. They represent a major theory of patterning in tissue organisation, due to their remarkable similarity to some natural patterns, such as skin pigmentation in zebrafish2, digit spacing3,4, and many others. The involvement of Turing patterns in biology has been debated because of their stringent fine-tuning requirements, where patterns only occur within a small subset of parameters5,6. This has complicated the engineering of a synthetic gene circuit for Turing patterns from first principles, even though natural genetic Turing networks have been successfully identified4,7. Here, we engineered a synthetic genetic reaction-diffusion system where three nodes interact according to a non-classical Turing network with improved parametric robustness6. The system was optimised inE. coliand reproducibly generated stationary, periodic, concentric stripe patterns in growing colonies. The patterns were successfully reproduced with a partial differential equation model, in a parameter regime obtained by fitting to experimental data. Our synthetic Turing system can contribute to novel nanotechnologies, such as patterned biomaterial deposition8,9, and provide insights into developmental patterning programs10.

Turing patterns are a fundamental concept in developmental biology, describing how homogeneous tissues develop into self-organized spatial patterns. However, the classical Turing mechanism, which relies on linear stability analysis, often fails to capture the complexities of real biological systems, such as multistability, non-linearities, growth, and boundary conditions. Here, we explore the impact of these factors on Turing pattern formation, contrasting linear stability analysis with numerical simulations based on a simple reaction-diffusion model, motivated by synthetic gene-regulatory pathways. We demonstrate how non-linearities introduce multistability, leading to unexpected pattern outcomes not predicted by the traditional Turing theory. The study also examines how growth and realistic boundary conditions influence pattern robustness, revealing that different growth regimes and boundary conditions can disrupt or stabilize pattern formation. Our findings are critical for understanding pattern formation in both natural and synthetic biological systems, providing insights into engineering robust patterns for applications in synthetic biology. Author summary During development, tissues self-organize to go from a single cell to a structured organism. In this process, simple chemical reactions lead to the emergence of the intricate designs we see in nature, like the stripes on a zebra or the labyrinths on a brain cortex. Although multiple theories have been proposed to model this phenomenon, one of the most simple and popular ones was introduced in the 1950s by the mathematician Alan Turing. However, his theory oversimplifies the biological conditions and ignores properties such as non-linearities, boundary effects, or growth in the tissue. In this work, we used a combination of mathematical models and computer simulations to investigate how these real-world factors influence pattern formation. Our findings show that when we account for these realistic effects, the patterns that emerge can be very different from what Turing’s theory would predict. Thus, this work may help us better understand the laws behind pattern formation and could have practical applications in tissue engineering for medical or environmental applications. ### Competing Interest Statement The authors have declared no competing interest.

Turing patterns are self-organizing systems that can form spots, stripes, or labyrinths. Proposed examples in tissue organization include zebrafish pigmentation, digit spacing, and many others. The theory of Turing patterns in biology has been debated because of their stringent fine-tuning requirements, where patterns only occur within a small subset of parameters. This has complicated the engineering of synthetic Turing gene circuits from first principles, although natural genetic Turing networks have been identified. Here, we engineered a synthetic genetic reaction-diffusion system where three nodes interact according to a non-classical Turing network with improved parametric robustness. The system reproducibly generated stationary, periodic, concentric stripe patterns in growing E. coli colonies. A partial differential equation model reproduced the patterns, with a Turing parameter regime obtained by fitting to experimental data. Our synthetic Turing system can contribute to nanotechnologies, such as patterned biomaterial deposition, and provide insights into developmental patterning programs. A record of this paper's transparent peer review process is included in the supplemental information.

Traditional linear stability analysis based on matrix diagonalization is a computationally intensive process for high-dimensional systems of differential equations, posing substantial limitations for the exploration of Turing systems of pattern formation where an additional wave-number parameter needs to be investigated. In this paper, we introduce an efficient and intuitive technique that leverages Gershgorin's theorem to determine upper limits on regions of parameter space and the wave number beyond which Turing instabilities cannot occur. This method offers a streamlined avenue for exploring the phase diagrams of other complex multi-parametric models, such as those found in gene regulatory networks in systems biology. Due to its suitability for the asymptotic limit of infinitely large systems, it predicts the existence of a sweet spot in network size for maximal Jacobian stability.

The diffusion-driven Turing instability is a potential mechanism for spatial pattern formation in numerous biological and chemical systems. However, engineering these patterns and demonstrating that they are produced by this mechanism is challenging. To address this, we aim to solve the inverse problem in artificial and experimental Turing patterns. This task is challenging since patterns are often corrupted by noise and slight changes in initial conditions can lead to different patterns. We used both least squares to explore the problem and physics-informed neural networks to build a noise-robust method. We elucidate the functionality of our network in scenarios mimicking biological noise levels and showcase its application using an experimentally obtained chemical pattern. The findings reveal the significant promise of machine learning in steering the creation of synthetic patterns in bioengineering, thereby advancing our grasp of morphological intricacies within biological systems while acknowledging existing limitations.

The type VI secretion system (T6SS) is an antibacterial weapon that is used by numerous Gram-negative bacteria to gain competitive advantage by injecting toxins into adjacent prey cells. Predicting the outcome of a T6SS-dependent competition is not only reliant on presence-absence of the system but instead involves a multiplicity of factors. Pseudomonas aeruginosa possesses 3 distinct T6SSs and a set of more than 20 toxic effectors with diverse functions including disruption of cell wall integrity, degradation of nucleic acids or metabolic impairment. We generated a comprehensive collection of mutants with various degrees of T6SS activity and/or sensitivity to each individual T6SS toxin. By imaging whole mixed bacterial macrocolonies, we then investigated how these P. aeruginosa strains gain a competitive edge in multiple attacker/prey combinations. We observed that the potency of single T6SS toxin varies significantly from one another as measured by monitoring the community structure, with some toxins acting better in synergy or requiring a higher payload. Remarkably the degree of intermixing between preys and attackers is also key to the competition outcome and is driven by the frequency of contact as well as the ability of the prey to move away from the attacker using type IV pili-dependent twitching motility. Finally, we implemented a computational model to better understand how changes in T6SS firing behaviours or cell-cell contacts lead to population level competitive advantages, thus providing conceptual insight applicable to all types of contact-based competition.

Clustering of membrane-associated molecules is thought to promote interactions with the actomyosin cortex, enabling size-dependent transport by actin flows. Consistent with this model, in the Caenorhabditis elegans zygote, efficient anterior segregation of the polarity protein PAR-3 requires oligomerization. However, through direct assessment of local coupling between motion of PAR proteins and the underlying cortex, we find no links between PAR-3 oligomer size and the degree of coupling. Indeed, both anterior and posterior PAR proteins experience similar advection velocities, at least over short distances. Consequently, differential cortex engagement cannot account for selectivity of PAR protein segregation by cortical flows. Combining experiment and theory, we demonstrate that a key determinant of differential segregation of PAR proteins by cortical flow is the stability of membrane association, which is enhanced by clustering and enables transport across cellular length scales. Thus, modulation of membrane binding dynamics allows cells to achieve selective transport by cortical flows despite widespread coupling between membrane-associated molecules and the cell cortex.

Traditional linear stability analysis based on matrix diagonalization is a computationally intensive $O(n^3)$ process for $n$-dimensional systems of differential equations, posing substantial limitations for the exploration of Turing systems of pattern formation where an additional wave-number parameter needs to be investigated. In this study, we introduce an efficient $O(n)$ technique that leverages Gershgorin's theorem to determine upper limits on regions of parameter space and the wave number beyond which Turing instabilities cannot occur. This method offers a streamlined avenue for exploring the phase diagrams of other complex multiparametric models, such as those found in systems biology.

T cells use sophisticated shape dynamics (morphodynamics) to migrate towards and neutralize infected and cancerous cells. However, there is limited quantitative understanding of the migration process in three-dimensional extracellular matrices (ECMs) and across timescales. Here, we leveraged recent advances in lattice light-sheet microscopy to quantitatively explore the three-dimensional morphodynamics of migrating T cells at high spatio-temporal resolution. We first developed a new shape descriptor based on spherical harmonics, incorporating key polarization information of the uropod. We found that the shape space of T cells is low-dimensional. At the behavioural level, run-and-stop migration modes emerge at approximately 150 s, and we mapped the morphodynamic composition of each mode using multiscale wavelet analysis, finding ‘stereotyped’ motifs. Focusing on the run mode, we found morphodynamics oscillating periodically (every approx. 100 s) that can be broken down into a biphasic process: front-widening with retraction of the uropod, followed by a rearward surface motion and forward extension, where intercalation with the ECM in both of these steps likely facilitates forward motion. Further application of these methods may enable the comparison of T cell migration across different conditions (e.g. differentiation, activation, tissues and drug treatments) and improve the precision of immunotherapeutic development.

Non-equilibrium thermodynamics has long been an area of substantial interest to ecologists because most fundamental biological processes, such as protein synthesis and respiration, are inherently energy-consuming. Microbial communities are a natural system to decipher this mechanistic basis because their interactions in the form of substrate consumption, metabolite production, and cross-feeding can be described explicitly in thermodynamic terms. Previous work has considered how thermodynamic constraints impact competition between pairs of species, but restrained from analysing how this manifests in complex dynamical systems. To address this gap, we develop a thermodynamic microbial community model with fully reversible reaction kinetics, which allows direct consideration of free-energy dissipation. This also allows species to interact via products rather than just substrates, increasing the dynamical complexity, and allowing a more nuanced classification of interaction types to emerge. Using this model, we find that community diversity increases with substrate lability, because greater free-energy availability allows for faster generation of niches. Thus, more niches are generated in the time frame of community establishment, leading to higher final species diversity. We also find that allowing species to make use of near-to-equilibrium reactions increases diversity in a low free-energy regime. In such a regime, two new thermodynamic interaction types that we identify here reach comparable strengths to the conventional (competition and facilitation) types, emphasising the key role that thermodynamics plays in community dynamics. Our results suggest that accounting for realistic thermodynamic constraints is vital for understanding the dynamics of real-world microbial communities.

Most biological processes are driven by non-equilibrium thermodynamics, but despite significant progress in theoretical ecology the constraints this places on ecosystem dynamics has been barely considered. Microbial ecosystems represent a natural place to begin this consideration, as many of the ways they interact, such as metabolite production and cross-feeding, can be described in thermodynamic terms. Previous work considered the impact of thermodynamics such as the rate-yield trade-off on individual species’ competitive ability, but restrained from analyzing complex dynamical systems. To address this gap we developed a thermodynamic microbial consumer-resource model with fully reversible reaction kinetics, which allows direct consideration of free-energy dissipation. Using this model, we show that ecosystem diversity increases with supplied free energy, because greater availability of free energy allows for faster ecosystem development. Thus, when species from the initial community begin to go extinct more possible niches have been formed, facilitating increased diversity. Our model also shows that the inclusion of species utilising near-to-equilibrium reactions increases diversity under conditions of low free-energy supply. At low free energy supply thermodynamic interaction types reach comparable strength to the conventional (competition and facilitation) interactions yielding a more nuanced classification of interactions, and emphasising the key role thermodynamics plays in the dynamics. Though our model is valid for all microbial ecosystems where diversification from an initial substrate occurs it is of particular use when the initial substrate is recalcitrant (low-free energy).

Medicines and agricultural biocides are often discovered using large phenotypic screens across hundreds of compounds, where visible effects of whole organisms are compared to gauge efficacy and possible modes of action. However, such analysis is often limited to human-defined and static features. Here, we introduce a novel framework that can characterize shape changes (morphodynamics) for cell-drug interactions directly from images, and use it to interpret perturbed development of Phakopsora pachyrhizi, the Asian soybean rust crop pathogen. We describe population development over a 2D space of shapes (morphospace) using two models with condition-dependent parameters: a top-down Fokker-Planck model of diffusive development over Waddington-type landscapes, and a bottom-up model of tip growth. We discover a variety of landscapes, describing phenotype transitions during growth, and identify possible perturbations in the tip growth machinery that cause this variation. This demonstrates a widely-applicable integration of unsupervised learning and biophysical modeling.

Multistable non-equilibrium systems are abundant outcomes of nonlinear dynamics with feedback, but still relatively little is known about what determines the stability of the steady states and their switching rates in terms of entropy and entropy production. Here, we will link fluctuation theorems for the entropy production along trajectories with the action obtainable from the Freidlin-Wentzell theorem to elucidate the thermodynamics of switching between states in the large volume limit of multistable systems. We find that the entropy production at steady state plays no role, but the entropy production during switching is key. Steady-state entropy and diffusive noise strength can be neglected in this limit. The relevance to biological, ecological, and climate models is apparent.

Motile cells have developed a variety of migration modes relying on diverse traction-force-generation mechanisms. Before the behavior of intracellular components could be easily imaged, cell movements were mostly classified by different types of cellular shape dynamics. Indeed, even though some types of cells move without any significant change in shape, most cell propulsion mechanisms rely on global or local deformations of the cell surface. In this review, focusing mostly on metazoan cells, we discuss how different types of local and global shape changes underlie distinct migration modes. We then discuss mechanical differences between force-generation mechanisms and finish by speculating on how they may have evolved.

Collective foraging has been shown to benefit organisms in environments where food is patchily distributed, but whether this is true in the case where organisms do not rely on long-range communications to coordinate their collective behaviour has been understudied. To address this question, we use the tractable laboratory model organismCaenorhabditis elegans, where a social strain (npr-1mutant) and a solitary strain (N2) are available for direct comparison of foraging strategies. We first developed an on-lattice minimal model for comparing collective and solitary foraging strategies, finding that social agents benefit from feeding faster and more efficiently simply owing to group formation. Our laboratory foraging experiments withnpr-1and N2 worm populations, however, show an advantage for solitary N2 in all food distribution environments that we tested. We incorporated additional strain-specific behavioural parameters ofnpr-1and N2 worms into our model and computationally identified N2's higher feeding rate to be the key factor underlying its advantage, without which it is possible to recapitulate the advantage of collective foraging in patchy environments. Our work highlights the theoretical advantage of collective foraging owing to group formation alone without long-range interactions and the valuable role of modelling to guide experiments. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.