In Vivo Microrheology Reveals Local Elastic and Plastic Responses Inside 3D Bacterial Biofilms.

Bacterial biofilms are highly abundant 3D living materials capable of performing complex biomechanical and biochemical functions, including programmable growth, self-repair, filtration, and bioproduction. Methods to measure internal mechanical properties of biofilms in vivo with spatial resolution on the cellular scale have been lacking. Here, thousands of cells are tracked inside living 3D biofilms of the bacterium Vibrio cholerae during and after the application of shear stress, for a wide range of stress amplitudes, periods, and biofilm sizes, which revealed anisotropic elastic and plastic responses of both cell displacements and cell reorientations. Using cellular tracking to infer parameters of a general mechanical model, spatially-resolved measurements of the elastic modulus inside the biofilm are obtained, which correlate with the spatial distribution of the polysaccharides within the biofilm matrix. The noninvasive microrheology and force-inference approach introduced here provides a general framework for studying mechanical properties with high spatial resolution in living materials.

Simultaneous spatiotemporal transcriptomics and microscopy of Bacillus subtilis swarm development reveal cooperation across generations.

Development of microbial communities is a complex multiscale phenomenon with wide-ranging biomedical and ecological implications. How biological and physical processes determine emergent spatial structures in microbial communities remains poorly understood due to a lack of simultaneous measurements of gene expression and cellular behaviour in space and time. Here we combined live-cell microscopy with a robotic arm for spatiotemporal sampling, which enabled us to simultaneously acquire phenotypic imaging data and spatiotemporal transcriptomes during Bacillus subtilis swarm development. Quantitative characterization of the spatiotemporal gene expression patterns revealed correlations with cellular and collective properties, and phenotypic subpopulations. By integrating these data with spatiotemporal metabolome measurements, we discovered a spatiotemporal cross-feeding mechanism fuelling swarm development: during their migration, earlier generations deposit metabolites which are consumed by later generations that swarm across the same location. These results highlight the importance of spatiotemporal effects during the emergence of phenotypic subpopulations and their interactions in bacterial communities.

Topological packing statistics of living and nonliving matter.

Complex disordered matter is of central importance to a wide range of disciplines, from bacterial colonies and embryonic tissues in biology to foams and granular media in materials science to stellar configurations in astrophysics. Because of the vast differences in composition and scale, comparing structural features across such disparate systems remains challenging. Here, by using the statistical properties of Delaunay tessellations, we introduce a mathematical framework for measuring topological distances between general three-dimensional point clouds. The resulting system-agnostic metric reveals subtle structural differences between bacterial biofilms as well as between zebrafish brain regions, and it recovers temporal ordering of embryonic development. We apply the metric to construct a universal topological atlas encompassing bacterial biofilms, snowflake yeast, plant shoots, zebrafish brain matter, organoids, and embryonic tissues as well as foams, colloidal packings, glassy materials, and stellar configurations. Living systems localize within a bounded island-like region of the atlas, reflecting that biological growth mechanisms result in characteristic topological properties.

Schrödinger Dynamics and Berry Phase of Undulatory Locomotion.

Spectral mode representations play an essential role in various areas of physics, from quantum mechanics to fluid turbulence, but they are not yet extensively used to characterize and describe the behavioral dynamics of living systems. Here, we show that mode-based linear models inferred from experimental live-imaging data can provide an accurate low-dimensional description of undulatory locomotion in worms, centipedes, robots, and snakes. By incorporating physical symmetries and known biological constraints into the dynamical model, we find that the shape dynamics are generically governed by Schrödinger equations in mode space. The eigenstates of the effective biophysical Hamiltonians and their adiabatic variations enable the efficient classification and differentiation of locomotion behaviors in natural, simulated, and robotic organisms using Grassmann distances and Berry phases. While our analysis focuses on a widely studied class of biophysical locomotion phenomena, the underlying approach generalizes to other physical or living systems that permit a mode representation subject to geometric shape constraints.

Proliferating active matter.

The fascinating patterns of collective motion created by autonomously driven particles have fuelled active-matter research for over two decades. So far, theoretical active-matter research has often focused on systems with a fixed number of particles. This constraint imposes strict limitations on what behaviours can and cannot emerge. However, a hallmark of life is the breaking of local cell number conservation by replication and death. Birth and death processes must be taken into account, for example, to predict the growth and evolution of a microbial biofilm, the expansion of a tumour, or the development from a fertilized egg into an embryo and beyond. In this Perspective, we argue that unique features emerge in these systems because proliferation represents a distinct form of activity: not only do the proliferating entities consume and dissipate energy, they also inject biomass and degrees of freedom capable of further self-proliferation, leading to myriad dynamic scenarios. Despite this complexity, a growing number of studies document common collective phenomena in various proliferating soft-matter systems. This generality leads us to propose proliferation as another direction of active-matter physics, worthy of a dedicated search for new dynamical universality classes. Conceptual challenges abound, from identifying control parameters and understanding large fluctuations and nonlinear feedback mechanisms to exploring the dynamics and limits of information flow in self-replicating systems. We believe that, by extending the rich conceptual framework developed for conventional active matter to proliferating active matter, researchers can have a profound impact on quantitative biology and reveal fascinating emergent physics along the way.

Ultrafast reversible self-assembly of living tangled matter.

Tangled active filaments are ubiquitous in nature, from chromosomal DNA and cilia carpets to root networks and worm collectives. How activity and elasticity facilitate collective topological transformations in living tangled matter is not well understood. We studied California blackworms (Lumbriculus variegatus), which slowly form tangles in minutes but can untangle in milliseconds. Combining ultrasound imaging, theoretical analysis, and simulations, we developed and validated a mechanistic model that explains how the kinematics of individual active filaments determines their emergent collective topological dynamics. The model reveals that resonantly alternating helical waves enable both tangle formation and ultrafast untangling. By identifying generic dynamical principles of topological self-transformations, our results can provide guidance for designing classes of topologically tunable active materials.

Learning hydrodynamic equations for active matter from particle simulations and experiments.

Recent advances in high-resolution imaging techniques and particle-based simulation methods have enabled the precise microscopic characterization of collective dynamics in various biological and engineered active matter systems. In parallel, data-driven algorithms for learning interpretable continuum models have shown promising potential for the recovery of underlying partial differential equations (PDEs) from continuum simulation data. By contrast, learning macroscopic hydrodynamic equations for active matter directly from experiments or particle simulations remains a major challenge, especially when continuum models are not known a priori or analytic coarse graining fails, as often is the case for nondilute and heterogeneous systems. Here, we present a framework that leverages spectral basis representations and sparse regression algorithms to discover PDE models from microscopic simulation and experimental data, while incorporating the relevant physical symmetries. We illustrate the practical potential through a range of applications, from a chiral active particle model mimicking nonidentical swimming cells to recent microroller experiments and schooling fish. In all these cases, our scheme learns hydrodynamic equations that reproduce the self-organized collective dynamics observed in the simulations and experiments. This inference framework makes it possible to measure a large number of hydrodynamic parameters in parallel and directly from video data.

Shared biophysical mechanisms determine early biofilm architecture development across different bacterial species.

Bacterial biofilms are among the most abundant multicellular structures on Earth and play essential roles in a wide range of ecological, medical, and industrial processes. However, general principles that govern the emergence of biofilm architecture across different species remain unknown. Here, we combine experiments, simulations, and statistical analysis to identify shared biophysical mechanisms that determine early biofilm architecture development at the single-cell level, for the species Vibrio cholerae, Escherichia coli, Salmonella enterica, and Pseudomonas aeruginosa grown as microcolonies in flow chambers. Our data-driven analysis reveals that despite the many molecular differences between these species, the biofilm architecture differences can be described by only 2 control parameters: cellular aspect ratio and cell density. Further experiments using single-species mutants for which the cell aspect ratio and the cell density are systematically varied, and mechanistic simulations show that tuning these 2 control parameters reproduces biofilm architectures of different species. Altogether, our results show that biofilm microcolony architecture is determined by mechanical cell-cell interactions, which are conserved across different species.

4-bit adhesion logic enables universal multicellular interface patterning.

Multicellular systems, from bacterial biofilms to human organs, form interfaces (or boundaries) between different cell collectives to spatially organize versatile functions1,2. The evolution of sufficiently descriptive genetic toolkits probably triggered the explosion of complex multicellular life and patterning3,4. Synthetic biology aims to engineer multicellular systems for practical applications and to serve as a build-to-understand methodology for natural systems5-8. However, our ability to engineer multicellular interface patterns2,9 is still very limited, as synthetic cell-cell adhesion toolkits and suitable patterning algorithms are underdeveloped5,7,10-13. Here we introduce a synthetic cell-cell adhesin logic with swarming bacteria and establish the precise engineering, predictive modelling and algorithmic programming of multicellular interface patterns. We demonstrate interface generation through a swarming adhesion mechanism, quantitative control over interface geometry and adhesion-mediated analogues of developmental organizers and morphogen fields. Using tiling and four-colour-mapping concepts, we identify algorithms for creating universal target patterns. This synthetic 4-bit adhesion logic advances practical applications such as human-readable molecular diagnostics, spatial fluid control on biological surfaces and programmable self-growing materials5-8,14. Notably, a minimal set of just four adhesins represents 4 bits of information that suffice to program universal tessellation patterns, implying a low critical threshold for the evolution and engineering of complex multicellular systems3,5.

Nematic bits and universal logic gates.

Liquid crystals (LCs) can host robust topological defect structures that essentially determine their optical and elastic properties. Although recent experimental progress enables precise control over nematic LC defects, their practical potential for information storage and processing has yet to be explored. Here, we introduce the concept of nematic bits (nbits) by exploiting a quaternionic mapping from LC defects to the Poincaré-Bloch sphere. Through theory and simulations, we demonstrate how single-nbit operations can be implemented using electric fields, to construct LC analogs of Pauli, Hadamard, and other elementary logic gates. Using nematoelastic interactions, we show how four-nbit configurations can realize universal classical NOR and NAND gates. Last, we demonstrate the implementation of generalized logical functions that take values on the Poincaré-Bloch sphere. These results open a route toward the implementation of classical digital and nonclassical continuous computation strategies in topological soft matter systems.

Odd dynamics of living chiral crystals.

Active crystals are highly ordered structures that emerge from the self-organization of motile objects, and have been widely studied in synthetic1,2 and bacterial3,4 active matter. Whether persistent  crystalline order can emerge  in groups of autonomously developing multicellular organisms is currently unknown. Here we show that swimming starfish embryos spontaneously assemble into chiral crystals that span thousands of spinning organisms and persist for tens of hours. Combining experiments, theory and simulations, we demonstrate that the formation, dynamics and dissolution of these living crystals are controlled by the hydrodynamic properties and the natural development of embryos. Remarkably, living chiral crystals exhibit self-sustained chiral oscillations as well as various unconventional deformation response behaviours recently predicted for odd elastic materials5,6. Our results provide direct experimental evidence for how non-reciprocal interactions between autonomous multicellular components may facilitate non-equilibrium phases of chiral active matter.

Learning developmental mode dynamics from single-cell trajectories.

Embryogenesis is a multiscale process during which developmental symmetry breaking transitions give rise to complex multicellular organisms. Recent advances in high-resolution live-cell microscopy provide unprecedented insights into the collective cell dynamics at various stages of embryonic development. This rapid experimental progress poses the theoretical challenge of translating high-dimensional imaging data into predictive low-dimensional models that capture the essential ordering principles governing developmental cell migration in complex geometries. Here, we combine mode decomposition ideas that have proved successful in condensed matter physics and turbulence theory with recent advances in sparse dynamical systems inference to realize a computational framework for learning quantitative continuum models from single-cell imaging data. Considering pan-embryo cell migration during early gastrulation in zebrafish as a widely studied example, we show how cell trajectory data on a curved surface can be coarse-grained and compressed with suitable harmonic basis functions. The resulting low-dimensional representation of the collective cell dynamics enables a compact characterization of developmental symmetry breaking and the direct inference of an interpretable hydrodynamic model, which reveals similarities between pan-embryo cell migration and active Brownian particle dynamics on curved surfaces. Due to its generic conceptual foundation, we expect that mode-based model learning can help advance the quantitative biophysical understanding of a wide range of developmental structure formation processes.

Estimating Entropy Production from Waiting Time Distributions.

Living systems operate far from thermal equilibrium by converting the chemical potential of ATP into mechanical work to achieve growth, replication, or locomotion. Given time series observations of intra-, inter-, or multicellular processes, a key challenge is to detect nonequilibrium behavior and quantify the rate of free energy consumption. Obtaining reliable bounds on energy consumption and entropy production directly from experimental data remains difficult in practice, as many degrees of freedom typically are hidden to the observer, so that the accessible coarse-grained dynamics may not obviously violate detailed balance. Here, we introduce a novel method for bounding the entropy production of physical and living systems which uses only the waiting time statistics of hidden Markov processes and, hence, can be directly applied to experimental data. By determining a universal limiting curve, we infer entropy production bounds from experimental data for gene regulatory networks, mammalian behavioral dynamics, and numerous other biological processes. Further considering the asymptotic limit of increasingly precise biological timers, we estimate the necessary entropic cost of heartbeat regulation in humans, dogs, and mice.

Emergence and melting of active vortex crystals.

Melting of two-dimensional (2D) equilibrium crystals is a complex phenomenon characterized by the sequential loss of positional and orientational order. In contrast to passive systems, active crystals can self-assemble and melt into an active fluid by virtue of their intrinsic motility and inherent non-equilibrium stresses. Currently, the non-equilibrium physics of active crystallization and melting processes is not well understood. Here, we establish the emergence and investigate the melting of self-organized vortex crystals in 2D active fluids using a generalized Toner-Tu theory. Performing extensive hydrodynamic simulations, we find rich transition scenarios. On small domains, we identify a hysteretic transition as well as a transition featuring temporal coexistence of active vortex lattices and active turbulence, both of which can be controlled by self-propulsion and active stresses. On large domains, an active vortex crystal with solid order forms within the parameter range corresponding to active vortex lattices. The melting of this crystal proceeds through an intermediate hexatic phase. Generally, these results highlight the differences and similarities between crystalline phases in active fluids and their equilibrium counterparts.

Active topolectrical circuits.

The transfer of topological concepts from the quantum world to classical mechanical and electronic systems has opened fundamentally different approaches to protected information transmission and wave guidance. A particularly promising emergent technology is based on recently discovered topolectrical circuits that achieve robust electric signal transduction by mimicking edge currents in quantum Hall systems. In parallel, modern active matter research has shown how autonomous units driven by internal energy reservoirs can spontaneously self-organize into collective coherent dynamics. Here, we unify key ideas from these two previously disparate fields to develop design principles for active topolectrical circuits (ATCs) that can self-excite topologically protected global signal patterns. Realizing autonomous active units through nonlinear Chua diode circuits, we theoretically predict and experimentally confirm the emergence of self-organized protected edge oscillations in one- and two-dimensional ATCs. The close agreement between theory, simulations, and experiments implies that nonlinear ATCs provide a robust and versatile platform for developing high-dimensional autonomous electrical circuits with topologically protected functionalities.

Emergent order in hydrodynamic spin lattices.

Macroscale analogues1-3 of microscopic spin systems offer direct insights into fundamental physical principles, thereby advancing our understanding of synchronization phenomena4 and informing the design of novel classes of chiral metamaterials5-7. Here we introduce hydrodynamic spin lattices (HSLs) of 'walking' droplets as a class of active spin systems with particle-wave coupling. HSLs reveal various non-equilibrium symmetry-breaking phenomena, including transitions from antiferromagnetic to ferromagnetic order that can be controlled by varying the lattice geometry and system rotation8. Theoretical predictions based on a generalized Kuramoto model4 derived from first principles rationalize our experimental observations, establishing HSLs as a versatile platform for exploring active phase oscillator dynamics. The tunability of HSLs suggests exciting directions for future research, from active spin-wave dynamics to hydrodynamic analogue computation and droplet-based topological insulators.

Topological braiding and virtual particles on the cell membrane.

Braiding of topological structures in complex matter fields provides a robust framework for encoding and processing information, and it has been extensively studied in the context of topological quantum computation. In living systems, topological defects are crucial for the localization and organization of biochemical signaling waves, but their braiding dynamics remain unexplored. Here, we show that the spiral wave cores, which organize the Rho-GTP protein signaling dynamics and force generation on the membrane of starfish egg cells, undergo spontaneous braiding dynamics. Experimentally measured world line braiding exponents and topological entropy correlate with cellular activity and agree with predictions from a generic field theory. Our analysis further reveals the creation and annihilation of virtual quasi-particle excitations during defect scattering events, suggesting phenomenological parallels between quantum and living matter.

Combinatorial patterns of graded RhoA activation and uniform F-actin depletion promote tissue curvature.

During development, gene expression regulates cell mechanics and shape to sculpt tissues. Epithelial folding proceeds through distinct cell shape changes that occur simultaneously in different regions of a tissue. Here, using quantitative imaging in Drosophila melanogaster, we investigate how patterned cell shape changes promote tissue bending during early embryogenesis. We find that the transcription factors Twist and Snail combinatorially regulate a multicellular pattern of lateral F-actin density that differs from the previously described Myosin-2 gradient. This F-actin pattern correlates with whether cells apically constrict, stretch or maintain their shape. We show that the Myosin-2 gradient and F-actin depletion do not depend on force transmission, suggesting that transcriptional activity is required to create these patterns. The Myosin-2 gradient width results from a gradient in RhoA activation that is refined through the balance between RhoGEF2 and the RhoGAP C-GAP. Our experimental results and simulations of a 3D elastic shell model show that tuning gradient width regulates tissue curvature.

Improved bounds on entropy production in living systems.

Living systems maintain or increase local order by working against the second law of thermodynamics. Thermodynamic consistency is restored as they consume free energy, thereby increasing the net entropy of their environment. Recently introduced estimators for the entropy production rate have provided major insights into the efficiency of important cellular processes. In experiments, however, many degrees of freedom typically remain hidden to the observer, and, in these cases, existing methods are not optimal. Here, by reformulating the problem within an optimization framework, we are able to infer improved bounds on the rate of entropy production from partial measurements of biological systems. Our approach yields provably optimal estimates given certain measurable transition statistics. In contrast to prevailing methods, the improved estimator reveals nonzero entropy production rates even when nonequilibrium processes appear time symmetric and therefore may pretend to obey detailed balance. We demonstrate the broad applicability of this framework by providing improved bounds on the energy consumption rates in a diverse range of biological systems including bacterial flagella motors, growing microtubules, and calcium oscillations within human embryonic kidney cells.