Defect-Mediated Morphologies in Growing Cell Colonies

Spatiotemporal control of structure and dynamics in a polar active fluid

We apply optimal control theory to a model of a polar active fluid (the Toner-Tu model), with the objective of driving the system into particular emergent dynamical behaviors or programming switching between states on demand. We use the effective self-propulsion speed as the control parameter (i.e. the means of external actuation). We identify control protocols that achieve outcomes such as relocating asters to targeted positions, forcing propagating solitary waves to reorient to a particular direction, and switching between stationary asters and propagating fronts. We analyze the solutions to identify generic principles for controlling polar active fluids. Our findings have implications for achieving spatiotemporal control of active polar systems in experiments, particularly in vitro cytoskeletal systems. Additionally, this research paves the way for leveraging optimal control methods to engineer the structure and dynamics of active fluids more broadly.

Growing bacterial colonies harness emergent genealogical demixing to regulate organizational entropy

Spatiotemporal organization of individuals within growing bacterial colonies is a key determinant of intraspecific interactions and colony-scale heterogeneities. The evolving cellular distribution, in relation to the genealogical lineage, is thus central to our understanding of bacterial fate across scales. Yet, how bacteria self-organize genealogically as a colony expands has remained unknown. Here, by developing a custom-built label-free algorithm, we track and study the genesis and evolution of emergent self-similar genealogical enclaves, whose dynamics are governed by biological activity. Topological defects at enclave boundaries tune finger-like morphologies of the active interfaces. The Shannon entropy of cell arrangements reduce over time; with faster-dividing cells possessing higher spatial affinity to genealogical relatives, at the cost of a well-mixed, entropically favorable state. Our coarse-grained lattice model demonstrates that genealogical enclaves emerge due to an interplay of division-mediated dispersal, stochasticity of division events, and cell-cell interactions. The study reports so-far hidden emergent self-organizing features arising due to entropic suppression, ultimately modulating intraspecific genealogical distances within bacterial colonies.

Structure and dynamics of motor-driven microtubule bundles

Connecting the large-scale emergent behaviors of active cytoskeletal materials to the microscopic properties of their constituents is a challenge due to a lack of data on the multiscale dynamics and structure of such systems. We approach this problem by studying the impact of depletion attraction on bundles of microtubules and kinesin-14 molecular motors. For all depletant concentrations, kinesin-14 bundles generate comparable extensile dynamics. However, this invariable mesoscopic behavior masks the transition in the microscopic motion of microtubules. Specifically, with increasing attraction, we observe a transition from bi-directional sliding with extension to pure extension with no sliding. Small-angle X-ray scattering shows that the transition in microtubule dynamics is concurrent with a structural rearrangement of microtubules from an open hexagonal to a compressed rectangular lattice. These results demonstrate that bundles of microtubules and molecular motors can display the same mesoscopic extensile behaviors despite having different internal structures and microscopic dynamics. They provide essential information for developing multiscale models of active matter.

Colonization and spreading dynamics of Lactiplantibacillus plantarum spoilage isolates on wet surfaces

The role of lactic acid bacteria, including Lactiplantibacillus plantarum, in food spoilage is well recognized, while the behavior of these non-motile bacteria on wet surfaces, such as those encountered in food processing environments has gained relatively little attention. Here, we observed a fast colony spreading of non-motile L. plantarum spoilage isolates on wet surfaces via passive sliding using solid BHI agar media as a model. We investigated the effect of physical properties of agar hydrogel substrate on the surface spreading of six L. plantarum food isolates FBR1-6 and a model strain WCFS1, using increasing concentrations of agar from 0.25 up to 1.5% (w/v). Our results revealed that L. plantarum strain FBR2 spreads significantly on low agar concentration plates compared to the other strains studied here (with a factor of 50-60 folds higher surface coverage), due to the formation of very soft biofilms with high water content that can float on the surface. The fast-spreading of FBR2 colonies is accompanied by an increased number of cells, elongated cell morphology, and a higher amount of extracellular components. Our finding highlights colonization dynamics and the spreading capacity of non-motile bacteria on surfaces that are relatively wet, thereby revealing an additional hitherto unnoticed parameter for non-motile bacteria that may contribute to contamination of foods by fast surface spreading of these bacteria in food processing environments.

Non-equilibrium structural and dynamic behaviors of active polymers in complex and crowded environments

Active matter systems, which convert internal chemical energy or energy from the environment into directed motion, are ubiquitous in nature and exhibit a range of emerging non-equilibrium behaviors. However, most of the current works on active matter have been devoted to particles, and the study of active polymers has only recently come into the spotlight due to their prevalence within living organisms. The intricate interplay between activity and conformational degrees of freedom gives rise to novel structural and dynamical behaviors of active polymers. Research in active polymers remarkably broadens diverse concepts of polymer physics, such as molecular architecture, dynamics, scaling and so on, which is of significant importance for the development of new polymer materials with unique performance. Furthermore, active polymers are often found in strongly interacting and crowded systems and in complex environments, so that the understanding of this behavior is essential for future developments of novel polymer-based biomaterials. This review thereby focuses on the study of active polymers in complex and crowded environments, and aims to provide insights into the fundamental physics underlying the adaptive and collective behaviors far from equilibrium, as well as the open challenges that the field is currently facing.

Machine learning topological defects in confluent tissues

Active nematics is an emerging paradigm for characterizing biological systems. One aspect of particularly intense focus is the role active nematic defects play in these systems, as they have been found to mediate a growing number of biological processes. Accurately detecting and classifying these defects in biological systems is, therefore, of vital importance to improving our understanding of such processes. While robust methods for defect detection exist for systems of elongated constituents, other systems, such as epithelial layers, are not well suited to such methods. Here, we address this problem by developing a convolutional neural network to detect and classify nematic defects in confluent cell layers. Crucially, our method is readily implementable on experimental images of cell layers and is specifically designed to be suitable for cells that are not rod shaped, which we demonstrate by detecting defects on experimental data using the trained model. We show that our machine learning model outperforms current defect detection techniques and that this manifests itself in our method as requiring less data to accurately capture defect properties. This could drastically improve the accuracy of experimental data interpretation while also reducing costs, advancing the study of nematic defects in biological systems.

Controlling active turbulence by activity patterns

By patterning activity in space, one can control active turbulence. To show this, we use Doi's hydrodynamic equations of a semidilute solution of active rods. A linear stability analysis reveals the resting isotropic fluid to be unstable above an absolute pusher activity. The emergent activity-induced paranematic state displays active turbulence, which we characterize by different quantities including the energy spectrum, which shows the typical power-law decay with exponent -4. Then, we control the active turbulence by a square lattice of circular spots where activity is switched off. In the parameter space lattice constant versus surface-to-surface distance of the spots, we identify different flow states. Most interestingly, for lattice constants below the vorticity correlation length and for spot distances smaller than the nematic coherence length, we observe a multi-lane flow state, where flow lanes with alternating flow directions are separated by a street of vortices. The flow pattern displays pronounced multistability and also appears transiently at the transition to the isotropic active-turbulence state. At larger lattice constants a trapped vortex state is identified with a non-Gaussian vorticity distribution due to the low flow vorticity at the spots. It transitions to conventional active turbulence for increasing spot distance.

[The Dynamic Model of the Active-Inactive Cell Interface]

Objective

To explore the morphodynamics of the active-inactive cell monolayer interfaces by using the active liquid crystal model.

Methods

A continuum mechanical model was established based on the active liquid crystal theory and the active-inactive cell monolayer interfaces were established by setting the activity difference of cell monolayers. The theoretical equations were solved numerically by the finite difference and the lattice Boltzmann method.

Results

The active-inactive cell interfaces displayed three typical morphologies, namely, flat interface, wavy interface, and finger-like interface. On the flat interfaces, the cells were oriented perpendicular to the interface, the -1/2 topological defects were clustered in the interfaces, and the interfaces were negatively charged. On the wavy interfaces, cells showed no obvious preference for orientation at the interfaces and the interfaces were neutrally charged. On the finger-like interfaces, cells were tangentially oriented at the interfaces, the +1/2 topological defects were collected at the interfaces, driving the growth of the finger-like structures, and the interfaces were positively charged.

Conclusion

The orientation of the cell alignment at the interface can significantly affect the morphologies of the active-inactive cell monolayer interfaces, which is closely associated with the dynamics of topological defects at the interfaces.

Biofilms as self-shaping growing nematics

Active nematics are the nonequilibrium analogue of passive liquid crystals. They consist of anisotropic units that consume free energy to drive emergent behaviour. Like liquid crystal molecules in displays, ordering and dynamics in active nematics are sensitive to boundary conditions. However, unlike passive liquid crystals, active nematics have the potential to regulate their boundaries through self-generated stresses. Here, we show how a three-dimensional, living nematic can actively shape itself and its boundary to regulate its internal architecture through growth-induced stresses, using bacterial biofilms confined by a hydrogel as a model system. We show that biofilms exhibit a sharp transition in shape from domes to lenses upon changing environmental stiffness or cell-substrate friction, which is explained by a theoretical model that considers the competition between confinement and interfacial forces. The growth mode defines the progression of the boundary, which in turn determines the trajectories and spatial distribution of cell lineages. We further demonstrate that the evolving boundary and corresponding stress anisotropy define the orientational ordering of cells and the emergence of topological defects in the biofilm interior. Our findings may provide strategies for the development of programmed microbial consortia with emergent material properties.

Instabilities in a growing system of active particles: scalar and vectorial systems

The physics of micron-scale biological colonies usually benefits from different out-of-equilibrium sources. In bacterial colonies and cellular tissues, the growth process is among the important active sources that determine the dynamics. In this article, we study the generic dynamical instabilities associated with the growth phenomena that may arise in both scalar and vectorial systems. In vectorial systems, where the rotational degrees of particles play a role, a phenomenological growth-mediated torque can affect the rotational dynamics of individual particles. We show that such a growth-mediated torque can result in active traveling waves in the bulk of a growing system. In addition to the bulk properties, we analyze the instabilities in the shape of growing interfaces in both scalar and vectorial systems.

Active nematics with deformable particles

The hydrodynamic theory of active nematics has been often used to describe the spatio-temporal dynamics of cell flows and motile topological defects within soft confluent tissues. Those theories, however, often rely on the assumption that tissues consist of cells with a fixed, anisotropic shape and do not resolve dynamical cell shape changes due to flow gradients. In this paper we extend the continuum theory of active nematics to include cell shape deformability. We find that circular cells in tissues must generate sufficient active stress to overcome an elastic barrier to deforming their shape in order to drive tissue-scale flows. Above this threshold the systems enter a dynamical steady-state with regions of elongated cells and strong flows coexisting with quiescent regions of isotropic cells.

Geometrical control of interface patterning underlies active matter invasion

Interaction between active materials and the boundaries of geometrical confinement is key to many emergent phenomena in active systems. For living active matter consisting of animal cells or motile bacteria, the confinement boundary is often a deformable interface, and it has been unclear how activity-induced interface dynamics might lead to morphogenesis and pattern formation. Here, we studied the evolution of bacterial active matter confined by a deformable boundary. We found that an ordered morphological pattern emerged at the interface characterized by periodically spaced interfacial protrusions; behind the interfacial protrusions, bacterial swimmers self-organized into multicellular clusters displaying +1/2 nematic defects. Subsequently, a hierarchical sequence of transitions from interfacial protrusions to creeping branches allowed the bacterial active drop to rapidly invade surrounding space with a striking self-similar branch pattern. We found that this interface patterning is geometrically controlled by the local curvature of the interface, a phenomenon we denote as collective curvature sensing. Using a continuum active model, we revealed that the collective curvature sensing arises from enhanced active stresses near high-curvature regions, with the active length scale setting the characteristic distance between the interfacial protrusions. Our findings reveal a protrusion-to-branch transition as a unique mode of active matter invasion and suggest a strategy to engineer pattern formation of active materials.

Friction-mediated phase transition in confined active nematics

Using a minimal continuum model, we investigate the interplay between circular confinement and substrate friction in active nematics. Upon increasing the friction from low to high, we observe a dynamical phase transition from a circulating flow phase to an anisotropic flow phase in which the flow tends to align perpendicular to the nematic director at the boundary. We demonstrate that both the flow structure and dynamic correlations in the latter phase differ from those of an unconfined, active turbulent system and may be controlled by the prescribed nematic boundary conditions. Our results show that substrate friction and geometric confinement act as valuable control parameters in active nematics.

Mechanical basis and topological routes to cell elimination

Cell layers eliminate unwanted cells through the extrusion process, which underlines healthy versus flawed tissue behaviors. Although several biochemical pathways have been identified, the underlying mechanical basis including the forces involved in cellular extrusion remains largely unexplored. Utilizing a phase-field model of a three-dimensional cell layer, we study the interplay of cell extrusion with cell-cell and cell-substrate interactions in a flat monolayer. Independent tuning of cell-cell versus cell-substrate adhesion forces reveals that extrusion events can be distinctly linked to defects in nematic and hexatic orders associated with cellular arrangements. Specifically, we show that by increasing relative cell-cell adhesion forces the cell monolayer can switch between the collective tendency towards fivefold, hexatic, disclinations relative to half-integer, nematic, defects for extruding a cell. We unify our findings by accessing three-dimensional mechanical stress fields to show that an extrusion event acts as a mechanism to relieve localized stress concentration.

Collective rotational motion of freely expanding T84 epithelial cell colonies

Coordinated rotational motion is an intriguing, yet still elusive mode of collective cell migration, which is relevant in pathological and morphogenetic processes. Most of the studies on this topic have been carried out on epithelial cells plated on micropatterned substrates, where cell motion is confined in regions of well-defined shapes coated with extracellular matrix adhesive proteins. The driver of collective rotation in such conditions has not been clearly elucidated, although it has been speculated that spatial confinement can play an essential role in triggering cell rotation. Here, we study the growth of epithelial cell colonies freely expanding (i.e. with no physical constraints) on the surface of cell culture plates and focus on collective cell rotation in such conditions, a case which has received scarce attention in the literature. One of the main findings of our work is that coordinated cell rotation spontaneously occurs in cell clusters in the free growth regime, thus implying that cell confinement is not necessary to elicit collective rotation as previously suggested. The extent of collective rotation was size and shape dependent: a highly coordinated disc-like rotation was found in small cell clusters with a round shape, while collective rotation was suppressed in large irregular cell clusters generated by merging of different clusters in the course of their growth. The angular motion was persistent in the same direction, although clockwise and anticlockwise rotations were equally likely to occur among different cell clusters. Radial cell velocity was quite low as compared to the angular velocity, in agreement with the free expansion regime where cluster growth is essentially governed by cell proliferation. A clear difference in morphology was observed between cells at the periphery and the ones in the core of the clusters, the former being more elongated and spread out as compared to the latter. Overall, our results, to our knowledge, provide the first quantitative and systematic evidence that coordinated cell rotation does not require a spatial confinement and occurs spontaneously in freely expanding epithelial cell colonies, possibly as a mechanism for the system.

Tilt-induced polar order and topological defects in growing bacterial populations

Rod-shaped bacteria, such as Escherichia coli, commonly live forming mounded colonies. They initially grow two-dimensionally on a surface and finally achieve three-dimensional growth. While it was recently reported that three-dimensional growth is promoted by topological defects of winding number +1/2 in populations of motile bacteria, how cellular alignment plays a role in nonmotile cases is largely unknown. Here, we investigate the relevance of topological defects in colony formation processes of nonmotile E. coli populations, and found that both ±1/2 topological defects contribute to the three-dimensional growth. Analyzing the cell flow in the bottom layer of the colony, we observe that +1/2 defects attract cells and -1/2 defects repel cells, in agreement with previous studies on motile cells, in the initial stage of the colony growth. However, later, cells gradually flow toward -1/2 defects as well, exhibiting a sharp contrast to the existing knowledge. By investigating three-dimensional cell orientations by confocal microscopy, we find that vertical tilting of cells is promoted near the defects. Crucially, this leads to the emergence of a polar order in the otherwise nematic two-dimensional cell orientation. We extend the theory of active nematics by incorporating this polar order and the vertical tilting, which successfully explains the influx toward -1/2 defects in terms of a polarity-induced force. Our work reveals that three-dimensional cell orientations may result in qualitative changes in properties of active nematics, especially those of topological defects, which may be generically relevant in active matter systems driven by cellular growth instead of self-propulsion.

Active nematic defects in compressible and incompressible flows

We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint. Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, while they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a smecticlike pattern, which induce an associated global polar ordering of +1/2 nematic defects. By contrast, divergence-free flows give rise to a local nematic order of the +1/2 defects, consisting of antialigned pairs of neighboring defects, as established in previous studies.

Stress anisotropy in confined populations of growing rods

A central feature of living matter is its ability to grow and multiply. The mechanical activity associated with growth produces both macroscopic flows shaped by confinement, and striking self-organization phenomena, such as orientational order and alignment, which are particularly prominent in populations of rod-shaped bacteria due to their nematic properties. However, how active stresses, passive mechanical interactions and flow-induced effects interact to give rise to the observed global alignment patterns remains elusive. Here, we study in silico colonies of growing rod-shaped particles of different aspect ratios confined in channel-like geometries. A spatially resolved analysis of the stress tensor reveals a strong relationship between near-perfect alignment and an inversion of stress anisotropy for particles with large length-to-width ratios. We show that, in quantitative agreement with an asymptotic theory, strong alignment can lead to a decoupling of active and passive stresses parallel and perpendicular to the direction of growth, respectively. We demonstrate the robustness of these effects in a geometry that provides less restrictive confinement and introduces natural perturbations in alignment. Our results illustrate the complexity arising from the inherent coupling between nematic order and active stresses in growing active matter, which is modulated by geometric and configurational constraints due to confinement.

Active nematic gel with quenched disorder

With quenched disorder, we introduce two-dimensional active nematics suspended in an incompressible fluid. We write the coarse-grained hydrodynamic equations of motion for slow variables, viz. density, orientation, and flow fields. The quenched disorder is introduced such that it interacts with the local orientation at every point with some strength. Disorder strength is tuned from zero to large values. We numerically study the defect dynamics and system kinetics and find that the finite disorder slows the ordering. The presence of fluid induces large fluctuation in the orientation field, further disturbing the ordering. The large fluctuation in the orientation field due to the fluid is so dominant that it reduces the effect of the quenched disorder. We have also found that the disorder effect is almost the same for both the contractile and extensile nature of active stresses in the system. This study can help to understand the impact of quenched disorder on the ordering kinetics of active gels with nematic interaction among the constituent objects.

Activity gradients in two- and three-dimensional active nematics

We numerically investigate how spatial variations of extensile or contractile active stress affect bulk active nematic systems in two and three dimensions. In the absence of defects, activity gradients drive flows which re-orient the nematic director field and thus act as an effective anchoring force. At high activity, defects are created and the system transitions into active turbulence, a chaotic flow state characterized by strong vorticity. We find that in two-dimensional (2D) systems active torques robustly align +1/2 defects parallel to activity gradients, with defect heads pointing towards contractile regions. In three-dimensional (3D) active nematics disclination lines preferentially lie in the plane perpendicular to activity gradients due to active torques acting on line segments. The average orientation of the defect structures in the plane perpendicular to the line tangent depends on the defect type, where wedge-like +1/2 defects align parallel to activity gradients, while twist defects are aligned anti-parallel. Understanding the response of active nematic fluids to activity gradients is an important step towards applying physical theories to biology, where spatial variations of active stress impact morphogenetic processes in developing embryos and affect flows and deformations in growing cell aggregates, such as tumours.

Defect absorption and emission for p-atic liquid crystals on cones

We investigate the ground-state configurations of two-dimensional liquid crystals with p-fold rotational symmetry (p-atics) on fixed curved surfaces. We focus on the intrinsic geometry and show that isothermal coordinates are particularly convenient as they explicitly encode a geometric contribution to the elastic potential. In the special case of a cone with half-angle β, the apex develops an effective topological charge of -χ, where 2πχ=2π(1-sinβ) is the deficit angle of the cone, and a topological defect of charge σ behaves as if it had an effective topological charge Q_{eff}=(σ-σ^{2}/2) when interacting with the apex. The effective charge of the apex leads to defect absorption and emission at the cone apex as the deficit angle of the cone is varied. For total topological defect charge 1, e.g., imposed by tangential boundary conditions at the edge, we find that for a disk the ground-state configuration consists of p defects each of charge +1/p lying equally spaced on a concentric ring of radius d=(p-1/3p-1)^{1/2p}R, where R is the radius of the disk. In the case of a cone with tangential boundary conditions at the base, we find three types of ground-state configurations as a function of cone angle: (i) for sharp cones, all of the +1/p defects are absorbed by the apex; (ii) at intermediate cone angles, some of the +1/p defects are absorbed by the apex and the rest lie equally spaced along a concentric ring on the flank; and (iii) for nearly flat cones, all of the +1/p defects lie equally spaced along a concentric ring on the flank. Here the defect positions and the absorption transitions depend intricately on p and the deficit angle, which we analytically compute. We check these results with numerical simulations for a set of commensurate cone angles and find excellent agreement.

Biotropic liquid crystal phase transformations in cellulose-producing bacterial communities

Biological functionality is often enabled by a fascinating variety of physical phenomena that emerge from orientational order of building blocks, a defining property of nematic liquid crystals that is also pervasive in nature. Out-of-equilibrium, "living" analogs of these technological materials are found in biological embodiments ranging from myelin sheath of neurons to extracellular matrices of bacterial biofilms and cuticles of beetles. However, physical underpinnings behind manifestations of orientational order in biological systems often remain unexplored. For example, while nematiclike birefringent domains of biofilms are found in many bacterial systems, the physics behind their formation is rarely known. Here, using cellulose-synthesizing Acetobacter xylinum bacteria, we reveal how biological activity leads to orientational ordering in fluid and gel analogs of these soft matter systems, both in water and on solid agar, with a topological defect found between the domains. Furthermore, the nutrient feeding direction plays a role like that of rubbing of confining surfaces in conventional liquid crystals, turning polydomain organization within the biofilms into a birefringent monocrystal-like order of both the extracellular matrix and the rod-like bacteria within it. We probe evolution of scalar orientational order parameters of cellulose nanofibers and bacteria associated with fluid-gel and isotropic-nematic transformations, showing how highly ordered active nematic fluids and gels evolve with time during biological-activity-driven, disorder-order transformation. With fluid and soft-gel nematics observed in a certain range of biological activity, this mesophase-exhibiting system is dubbed "biotropic," analogously to thermotropic nematics that exhibit solely orientational order within a temperature range, promising technological and fundamental-science applications.

Self-rotation regulates interface evolution in biphasic active matter through taming defect dynamics

Chirality can endow nonequilibrium active matter with unique features and functions. Here, we explore the chiral dynamics in biphasic active nematics composed of self-rotating units that continuously inject energy and angular momentum at the microscale. We show that the self-rotation of units can regularize the boundaries between two phases, rendering sinusoidal-like interfaces, which allow lateral wave propagation and are characterized by chains of ordered antiferromagnetic cross-interface flow vortices. Through the spontaneous coordination of counter-rotating units across the interfaces, topological defects excited by activity are sorted spatiotemporally, where positive defects are locally trapped at the interfaces but, unexpectedly, are transported laterally in a unidirectional rather than wavy mode, whereas inertial negative defects remain spinning in the bulks. Our findings reveal that individual chirality could be harnessed to modulate interfacial morphodynamics in active systems and suggest a potential approach toward controlling topological defects for programmable microfluidics and logic operations.

Spreading rates of bacterial colonies depend on substrate stiffness and permeability

The ability of bacteria to colonize and grow on different surfaces is an essential process for biofilm development. Here, we report the use of synthetic hydrogels with tunable stiffness and porosity to assess physical effects of the substrate on biofilm development. Using time-lapse microscopy to track the growth of expanding Serratia marcescens colonies, we find that biofilm colony growth can increase with increasing substrate stiffness, unlike what is found on traditional agar substrates. Using traction force microscopy-based techniques, we find that biofilms exert transient stresses correlated over length scales much larger than a single bacterium, and that the magnitude of these forces also increases with increasing substrate stiffness. Our results are consistent with a model of biofilm development in which the interplay between osmotic pressure arising from the biofilm and the poroelastic response of the underlying substrate controls biofilm growth and morphology.

Cellular Sensing Governs the Stability of Chemotactic Fronts

In contexts ranging from embryonic development to bacterial ecology, cell populations migrate chemotactically along self-generated chemical gradients, often forming a propagating front. Here, we theoretically show that the stability of such chemotactic fronts to morphological perturbations is determined by limitations in the ability of individual cells to sense and thereby respond to the chemical gradient. Specifically, cells at bulging parts of a front are exposed to a smaller gradient, which slows them down and promotes stability, but they also respond more strongly to the gradient, which speeds them up and promotes instability. We predict that this competition leads to chemotactic fingering when sensing is limited at too low chemical concentrations. Guided by this finding and by experimental data on E. coli chemotaxis, we suggest that the cells' sensory machinery might have evolved to avoid these limitations and ensure stable front propagation. Finally, as sensing of any stimuli is necessarily limited in living and active matter in general, the principle of sensing-induced stability may operate in other types of directed migration such as durotaxis, electrotaxis, and phototaxis.

Phase separation dynamics in deformable droplets

Phase separation can drive spatial organization of multicomponent mixtures. For instance in developing animal embryos, effective phase separation descriptions have been used to account for the spatial organization of different tissue types. Similarly, separation of different tissue types is also observed in stem cell aggregates, where the emergence of a polar organization can mimic early embryonic axis formation. Here, we describe such aggregates as deformable two-phase fluid droplets, which are suspended in a fluid environment (third phase). Using hybrid finite-volume Lattice-Boltzmann simulations, we numerically explore the out-of-equilibrium routes that can lead to the polar equilibrium state of such a droplet. We focus on the interplay between spinodal decomposition and advection with hydrodynamic flows driven by interface tensions, which we characterize by a Peclet number Pe. Consistent with previous work, for large Pe the coarsening process is generally accelerated. However, for intermediate Pe we observe long-lived, strongly elongated droplets, where both phases form an alternating stripe pattern. We show that these "croissant" states are close to mechanical equilibrium and coarsen only slowly through diffusive fluxes in an Ostwald-ripening-like process. Finally, we show that a surface tension asymmetry between both droplet phases leads to transient, rotationally symmetric states whose resolution leads to flows reminiscent of Marangoni flows. Our work highlights the importance of advection for the phase separation process in finite, deformable systems.

Chemotactic smoothing of collective migration

Collective migration-the directed, coordinated motion of many self-propelled agents-is a fascinating emergent behavior exhibited by active matter with functional implications for biological systems. However, how migration can persist when a population is confronted with perturbations is poorly understood. Here, we address this gap in knowledge through studies of bacteria that migrate via directed motion, or chemotaxis, in response to a self-generated nutrient gradient. We find that bacterial populations autonomously smooth out large-scale perturbations in their overall morphology, enabling the cells to continue to migrate together. This smoothing process arises from spatial variations in the ability of cells to sense and respond to the local nutrient gradient-revealing a population-scale consequence of the manner in which individual cells transduce external signals. Altogether, our work provides insights to predict, and potentially control, the collective migration and morphology of cellular populations and diverse other forms of active matter.

Large-scale orientational order in bacterial colonies during inward growth

During colony growth, complex interactions regulate the bacterial orientation, leading to the formation of large-scale ordered structures, including topological defects, microdomains, and branches. These structures may benefit bacterial strains, providing invasive advantages during colonization. Active matter dynamics of growing colonies drives the emergence of these ordered structures. However, additional biomechanical factors also play a significant role during this process. Here, we show that the velocity profile of growing colonies creates strong radial orientation during inward growth when crowded populations invade a closed area. During this process, growth geometry sets virtual confinement and dictates the velocity profile. Herein, flow-induced alignment and torque balance on the rod-shaped bacteria result in a new stable orientational equilibrium in the radial direction. Our analysis revealed that the dynamics of these radially oriented structures, also known as aster defects, depend on bacterial length and can promote the survival of the longest bacteria around localized nutritional hotspots. The present results indicate a new mechanism underlying structural order and provide mechanistic insights into the dynamics of bacterial growth on complex surfaces.

Polar Fluctuations Lead to Extensile Nematic Behavior in Confluent Tissues

How can a collection of motile cells, each generating contractile nematic stresses in isolation, become an extensile nematic at the tissue level? Understanding this seemingly contradictory experimental observation, which occurs irrespective of whether the tissue is in the liquid or solid states, is not only crucial to our understanding of diverse biological processes, but is also of fundamental interest to soft matter and many-body physics. Here, we resolve this cellular to tissue level disconnect in the small fluctuation regime by using analytical theories based on hydrodynamic descriptions of confluent tissues, in both liquid and solid states. Specifically, we show that a collection of microscopic constituents with no inherently nematic extensile forces can exhibit active extensile nematic behavior when subject to polar fluctuating forces. We further support our findings by performing cell level simulations of minimal models of confluent tissues.

Active Extensile Stress Promotes 3D Director Orientations and Flows

We use numerical simulations and linear stability analysis to study an active nematic layer where the director is allowed to point out of the plane. Our results highlight the difference between extensile and contractile systems. Contractile stress suppresses the flows perpendicular to the layer and favors in-plane orientations of the director. By contrast extensile stress promotes instabilities that can turn the director out of the plane, leaving behind a population of distinct, in-plane regions that continually elongate and divide. This supports extensile forces as a mechanism for the initial stages of layer formation in living systems, and we show that a planar drop with extensile (contractile) activity grows into three dimensions (remains in two dimensions). The results also explain the propensity of disclination lines in three dimensional active nematics to be of twist type in extensile or wedge type in contractile materials.

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.

Mechanical forces drive a reorientation cascade leading to biofilm self-patterning

In growing active matter systems, a large collection of engineered or living autonomous units metabolize free energy and create order at different length scales as they proliferate and migrate collectively. One such example is bacterial biofilms, surface-attached aggregates of bacterial cells embedded in an extracellular matrix that can exhibit community-scale orientational order. However, how bacterial growth coordinates with cell-surface interactions to create distinctive, long-range order during biofilm development remains elusive. Here we report a collective cell reorientation cascade in growing Vibrio cholerae biofilms that leads to a differentially ordered, spatiotemporally coupled core-rim structure reminiscent of a blooming aster. Cell verticalization in the core leads to a pattern of differential growth that drives radial alignment of the cells in the rim, while the growing rim generates compressive stresses that expand the verticalized core. Such self-patterning disappears in nonadherent mutants but can be restored through opto-manipulation of growth. Agent-based simulations and two-phase active nematic modeling jointly reveal the strong interdependence of the driving forces underlying the differential ordering. Our findings offer insight into the developmental processes that shape bacterial communities and provide ways to engineer phenotypes and functions in living active matter.

Physics of liquid crystals in cell biology

The past decade has witnessed a rapid growth in understanding of the pivotal roles of mechanical stresses and physical forces in cell biology. As a result, an integrated view of cell biology is evolving, where genetic and molecular features are scrutinised hand in hand with physical and mechanical characteristics of cells. Physics of liquid crystals has emerged as a burgeoning new frontier in cell biology over the past few years, fuelled by an increasing identification of orientational order and topological defects in cell biology, spanning scales from subcellular filaments to individual cells and multicellular tissues. Here, we provide an account of the most recent findings and developments, together with future promises and challenges in this rapidly evolving interdisciplinary research direction.

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.

Topological defects in the mesothelium suppress ovarian cancer cell clearance

We investigated an in vitro model for mesothelial clearance, wherein ovarian cancer cells invade into a layer of mesothelial cells, resulting in mesothelial retraction combined with cancer cell disaggregation and spreading. Prior to the addition of tumor cells, the mesothelial cells had an elongated morphology, causing them to align with their neighbors into well-ordered domains. Flaws in this alignment, which occur at topological defects, have been associated with altered cell density, motion, and forces. Here, we identified topological defects in the mesothelial layer and showed how they affected local cell density by producing a net flow of cells inward or outward, depending on the defect type. At locations of net inward flow, mesothelial clearance was impeded. Hence, the collective behavior of the mesothelial cells, as governed by the topological defects, affected tumor cell clearance and spreading. Importantly, our findings were consistent across multiple ovarian cancer cell types, suggesting a new physical mechanism that could impact ovarian cancer metastasis.

Morphogenesis and cell ordering in confined bacterial biofilms

Biofilms are aggregates of bacterial cells surrounded by an extracellular matrix. Much progress has been made in studying biofilm growth on solid substrates; however, little is known about the biophysical mechanisms underlying biofilm development in three-dimensional confined environments in which the biofilm-dwelling cells must push against and even damage the surrounding environment to proliferate. Here, combining single-cell imaging, mutagenesis, and rheological measurement, we reveal the key morphogenesis steps of Vibrio cholerae biofilms embedded in hydrogels as they grow by four orders of magnitude from their initial size. We show that the morphodynamics and cell ordering in embedded biofilms are fundamentally different from those of biofilms on flat surfaces. Treating embedded biofilms as inclusions growing in an elastic medium, we quantitatively show that the stiffness contrast between the biofilm and its environment determines biofilm morphology and internal architecture, selecting between spherical biofilms with no cell ordering and oblate ellipsoidal biofilms with high cell ordering. When embedded in stiff gels, cells self-organize into a bipolar structure that resembles the molecular ordering in nematic liquid crystal droplets. In vitro biomechanical analysis shows that cell ordering arises from stress transmission across the biofilm-environment interface, mediated by specific matrix components. Our imaging technique and theoretical approach are generalizable to other biofilm-forming species and potentially to biofilms embedded in mucus or host tissues as during infection. Our results open an avenue to understand how confined cell communities grow by means of a compromise between their inherent developmental program and the mechanical constraints imposed by the environment.

Multicomponent superconducting order parameter in UTe2

An unconventional superconducting state was recently discovered in uranium ditelluride (UTe₂), in which spin-triplet superconductivity emerges from the paramagnetic normal state of a heavy-fermion material. The coexistence of magnetic fluctuations and superconductivity, together with the crystal structure of this material, suggests that a distinctive set of symmetries, magnetic properties, and topology underlie the superconducting state. Here, we report observations of a nonzero polar Kerr effect and of two transitions in the specific heat upon entering the superconducting state, which together suggest that the superconductivity in UTe₂ is characterized by a two-component order parameter that breaks time-reversal symmetry. These data place constraints on the symmetries of the order parameter and inform the discussion on the presence of topological superconductivity in UTe₂.

Active flows and deformable surfaces in development

We review progress in active hydrodynamic descriptions of flowing media on curved and deformable manifolds: the state-of-the-art in continuum descriptions of single-layers of epithelial and/or other tissues during development. First, after a brief overview of activity, flows and hydrodynamic descriptions, we highlight the generic challenge of identifying the dependence on dynamical variables of so-called active kinetic coefficients- active counterparts to dissipative Onsager coefficients. We go on to describe some of the subtleties concerning how curvature and active flows interact, and the issues that arise when surfaces are deformable. We finish with a broad discussion around the utility of such theories in developmental biology. This includes limitations to analytical techniques, challenges associated with numerical integration, fitting-to-data and inference, and potential tools for the future, such as discrete differential geometry.

Topological defects of integer charge in cell monolayers

Many cell types spontaneously order like nematic liquid crystals, and, as such, they form topological defects, which influence the cell organization. While defects with topological charge ±1/2 are common in cell monolayers, defects with charge ±1, which are thought to be relevant in the formation of protrusions in living systems, are more elusive. We use topographical patterns to impose topological charge of ±1 in controlled locations in cell monolayers. We study two types of cells, 3T6 fibroblasts and EpH-4 epithelial cells, and we compare their behavior on such patterns, characterizing the degree of alignment, the cell density near the defects, and their behavior at the defect core. We observe density variation in the 3T6 monolayers near both types of defects over the same length-scale. By choosing appropriate geometrical parameters of our topographical features, we identify a new behavior of 3T6 cells near the defects with topological charge +1, leading to a change in the cells' preferred shape. Our strategy allows a fine control of cell alignment near defects as a platform to study liquid crystalline properties of cells.

Defect Spirograph: Dynamical Behavior of Defects in Spatially Patterned Active Nematics

Topological defects in active liquid crystals can be confined by introducing gradients of activity. Here, we examine the dynamical behavior of two defects confined by a sharp gradient of activity that separates an active circular region and a surrounding passive nematic material. Continuum simulations are used to explain how the interplay among energy injection into the system, hydrodynamic interactions, and frictional forces governs the dynamics of topologically required self-propelling +1/2 defects. Our findings are rationalized in terms of a phase diagram for the dynamical response of defects in terms of activity and frictional damping strength. Different regions of the underlying phase diagram correspond to distinct dynamical modes, namely immobile defects, steady rotation of defects, bouncing defects, bouncing-cruising defects, dancing defects, and multiple defects with irregular dynamics. These dynamic states raise the prospect of generating synchronized defect arrays for microfluidic applications.

Activity pulses induce spontaneous flow reversals in viscoelastic environments

Complex interactions between cellular systems and their surrounding extracellular matrices are emerging as important mechanical regulators of cell functions, such as proliferation, motility and cell death, and such cellular systems are often characterized by pulsating actomyosin activities. Here, using an active gel model, we numerically explore spontaneous flow generation by activity pulses in the presence of a viscoelastic medium. The results show that cross-talk between the activity-induced deformations of the viscoelastic surroundings and the time-dependent response of the active medium to these deformations can lead to the reversal of spontaneously generated active flows. We explain the mechanism behind this phenomenon based on the interaction between the active flow and the viscoelastic medium. We show the importance of relaxation time scales of both the polymers and the active particles and provide a phase space over which such spontaneous flow reversals can be observed. Our results suggest new experiments investigating the role of controlled pulses of activity in living systems ensnared in complex mircoenvironments.

Fluctuations can induce local nematic order and extensile stress in monolayers of motile cells

Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix.

How Physical Interactions Shape Bacterial Biofilms

Biofilms are structured communities formed by a single or multiple microbial species. Within biofilms, bacteria are embedded into extracellular matrix, allowing them to build macroscopic objects. Biofilm structure can respond to environmental changes such as the presence of antibiotics or predators. By adjusting expression levels of surface and extracellular matrix components, bacteria tune cell-to-cell interactions. One major challenge in the field is the fact that these components are very diverse among different species. Deciphering how physical interactions within biofilms are affected by changes in gene expression is a promising approach to obtaining a more unified picture of how bacteria modulate biofilms. This review focuses on recent advances in characterizing attractive and repulsive forces between bacteria in correlation with biofilm structure, dynamics, and spreading. How bacteria control physical interactions to maximize their fitness is an emerging theme.

Robustness of topological defects in discrete domains

Topological defects are singular points in vector fields, important in applications ranging from fingerprint detection to liquid crystals to biomedical imaging. In discretized vector fields, topological defects and their topological charge are identified by finite differences or finite-step paths around the tentative defect. As the topological charge is (half) integer, it cannot depend continuously on each input vector in a discrete domain. Instead, it switches discontinuously when vectors change beyond a certain amount, making the analysis of topological defects error prone in noisy data. We improve existing methods for the identification of topological defects by proposing a robustness measure for (i) the location of a defect, (ii) the existence of topological defects and the total topological charge within a given area, (iii) the annihilation of a defect pair, and (iv) the formation of a defect pair. Based on the proposed robustness measure, we show that topological defects in discrete domains can be identified with optimal trade-off between localization precision and robustness. The proposed robustness measure enables uncertainty quantification for topological defects in noisy discretized nematic fields (orientation fields) and polar fields (vector fields).

Epithelial colonies in vitro elongate through collective effects

Epithelial tissues of the developing embryos elongate by different mechanisms, such as neighbor exchange, cell elongation, and oriented cell division. Since autonomous tissue self-organization is influenced by external cues such as morphogen gradients or neighboring tissues, it is difficult to distinguish intrinsic from directed tissue behavior. The mesoscopic processes leading to the different mechanisms remain elusive. Here, we study the spontaneous elongation behavior of spreading circular epithelial colonies in vitro. By quantifying deformation kinematics at multiple scales, we report that global elongation happens primarily due to cell elongations, and its direction correlates with the anisotropy of the average cell elongation. By imposing an external time-periodic stretch, the axis of this global symmetry breaking can be modified and elongation occurs primarily due to orientated neighbor exchange. These different behaviors are confirmed using a vertex model for collective cell behavior, providing a framework for understanding autonomous tissue elongation and its origins.

Active nematics with quenched disorder

We introduce a two-dimensional active nematic with quenched disorder. We write the coarse-grained hydrodynamic equations of motion for slow variables, viz. density and orientation. Disorder strength is tuned from zero to large values. Results from the numerical solution of equations of motion as well as the calculation of two-point orientation correlation function using linear approximation shows that the ordered steady state follows a disorder dependent crossover from quasi-long-range order to short-range order. Such crossover is due to the pinning of ±1/2 topological defects in the presence of finite disorder, which breaks the system in uncorrelated domains. Finite disorder slows the dynamics of +1/2 defect, and it leads to slower growth dynamics. The two-point correlation functions for the density and orientation fields show good dynamic scaling but no static scaling for the different disorder strengths. Our findings can motivate experimentalists to verify the results and find applications in living and artificial apolar systems in the presence of a quenched disorder.

Role of Friction in Multidefect Ordering

We use continuum simulations to study the impact of friction on the ordering of defects in an active nematic. Even in a frictionless system, +1/2 defects tend to align side by side and orient antiparallel reflecting their propensity to form, and circulate with, flow vortices. Increasing friction enhances the effectiveness of the defect-defect interactions, and defects form dynamically evolving, large-scale, positionally, and orientationally ordered structures, which can be explained as a competition between hexagonal packing, preferred by the -1/2 defects, and rectangular packing, preferred by the +1/2 defects.

Optimal Control of Active Nematics

In this work we present the first systematic framework to sculpt active nematic systems, using optimal control theory and a hydrodynamic model of active nematics. We demonstrate the use of two different control fields, (i) applied vorticity and (ii) activity strength, to shape the dynamics of an extensile active nematic that is confined to a disk. In the absence of control inputs, the system exhibits two attractors, clockwise and counterclockwise circulating states characterized by two co-rotating topological +1/2 defects. We specifically seek spatiotemporal inputs that switch the system from one attractor to the other; we also examine phase-shifting perturbations. We identify control inputs by optimizing a penalty functional with three contributions: total control effort, spatial gradients in the control, and deviations from the desired trajectory. This work demonstrates that optimal control theory can be used to calculate nontrivial inputs capable of restructuring active nematics in a manner that is economical, smooth, and rapid, and therefore will serve as a guide to experimental efforts to control active matter.

Universality class for a nonequilibrium state of matter: A d=4-ε expansion study of Malthusian flocks

We show that "Malthusian flocks"-i.e., coherently moving collections of self-propelled entities (such as living creatures) which are being "born" and "dying" during their motion-belong to a new universality class in spatial dimensions d>2. We calculate the universal exponents and scaling laws of this new universality class to O(ε) in a d=4-ε expansion and find these are different from the "canonical" exponents previously conjectured to hold for "immortal" flocks (i.e., those without birth and death) and shown to hold for incompressible flocks with spatial dimensions in the range of 2<d≤4. We also obtain a universal amplitude ratio relating the damping of transverse and longitudinal velocity and density fluctuations in these systems. Furthermore, we find a universal separatrix in real space (r) between two regions in which the equal-time density correlation 〈δρ(r,t)δρ(0,t)〉 has opposite signs. Our expansion should be quite accurate in d=3, allowing precise quantitative comparisons between our theory, simulations, and experiments.

Moving, Reproducing, and Dying Beyond Flatland: Malthusian Flocks in Dimensions d>2

We show that "Malthusian flocks"-i.e., coherently moving collections of self-propelled entities (such as living creatures) which are being "born" and "dying" during their motion-belong to a new universality class in spatial dimensions d>2. We calculate the universal exponents and scaling laws of this new universality class to O(ε) in an ε=4-d expansion, and find these are different from the "canonical" exponents previously conjectured to hold for "immortal" flocks (i.e., those without birth and death) and shown to hold for incompressible flocks in d>2. Our expansion should be quite accurate in d=3, allowing precise quantitative comparisons between our theory, simulations, and experiments.