Theory of defect-mediated morphogenesis

Energetic scaling behavior of patterned epithelium

Cellular monolayers display various degrees of coordinated motion ranging from the small scale of just a few cells to large multi-cellular scales. This collective migration carries important physical cues for creating proper tissue morphology. Previous studies have demonstrated that the energetics of the epithelial monolayer show a linear variation with time in conjunction with an arrest in monolayer motion after confluency. However, little is known about how the energetics of monolayer development are affected by confined geometries. Here, we demonstrate that micropatterned epithelial monolayers display a non-linear change in energetic variables, which coincides with the large-scale coordination of migration. This non-linear scaling behavior was further seen to be associated with the biased alignment of cells and cell-cell adhesion. These findings provide a new understanding of how developing epithelia may be impacted by different conditions in vivo.

Exploring the Synthesis of Self-Organization and Active Motion

Proteins, genetic material, and membranes are fundamental to all known organisms, yet these components alone do not constitute life. Life emerges from the dynamic processes of self-organization, assembly, and active motion, suggesting the existence of similar artificial systems. Against this backdrop, our Perspective explores a variety of chemical phenomena illustrating how nonequilibrium self-organization and micromotors contribute to life-like behavior and functionalities. After explaining key terms, we discuss specific examples including enzymatic motion, diffusiophoretic and bubble-driven self-propulsion, pattern-forming reaction-diffusion systems, self-assembling inorganic aggregates, and hierarchically emergent phenomena. We also provide a roadmap for combining self-organization and active motion and discuss possible outcomes through biological analogs. We suggest that this research direction, deeply rooted in physical chemistry, offers opportunities for further development with broad impacts on related sciences and technologies.

Agent-based modeling of stress anisotropy driven nematic ordering in growing biofilms

Living active collectives have evolved with remarkable self-patterning capabilities to adapt to the physical and biological constraints crucial for their growth and survival. However, the intricate process by which complex multicellular patterns emerge from a single founder cell remains elusive. In this study, we utilize an agent-based model, validated through single-cell microscopy imaging, to track the three-dimensional (3D) morphodynamics of cells within growing bacterial biofilms encased by agarose gels. The confined growth conditions give rise to a spatiotemporally heterogeneous stress landscape within the biofilm. In the core of the biofilm, where high hydrostatic and low shear stresses prevail, cell packing appears disordered. In contrast, near the gel-cell interface, a state of high shear stress and low hydrostatic stress emerges, driving nematic ordering, albeit with a time delay inherent to shear stress relaxation. Strikingly, we observe a robust spatiotemporal correlation between stress anisotropy and nematic ordering within these confined biofilms. This correlation suggests a mechanism whereby stress anisotropy plays a pivotal role in governing the spatial organization of cells. The reciprocity between stress anisotropy and cell ordering in confined biofilms opens new avenues for innovative 3D mechanically guided patterning techniques for living active collectives, which hold significant promise for a wide array of environmental and biomedical applications.

Three-Dimensional Chiral Morphogenesis of Active Fluids

Chirality is an essential nature of biological systems. However, it remains obscure how the handedness at the microscale is translated into chiral morphogenesis at the tissue level. Here, we investigate three-dimensional (3D) tissue morphogenesis using an active fluid theory invoking chirality. We show that the coordination of achiral and chiral stresses, arising from microscopic interactions and energy input of individual cells, can engender the self-organization of 3D papillary and helical structures. The achiral active stress drives the nucleation of asterlike topological defects, which initiate 3D out-of-plane budding, followed by rodlike elongation. The chiral active stress excites vortexlike topological defects, which favor the tip spheroidization and twisting of the elongated rod. These results unravel the chiral morphogenesis observed in our experiments of 3D organoids generated by human embryonic stem cells.

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.

Odd elasticity and topological waves in active surfaces

Odd elasticity describes active elastic systems whose stress-strain relationship is not compatible with a potential energy. As the requirement of energy conservation is lifted from linear elasticity, new antisymmetric (odd) components appear in the elastic tensor. In this work we study the odd elasticity and non-Hermitian wave dynamics of active surfaces, specifically plates of moderate thickness. These odd moduli can endow the vibrational modes of the plate with a nonzero topological invariant known as the first Chern number. Within continuum elastic theory, we show that the Chern number is related to the presence of unidirectional shearing waves that are hosted at the plate's boundary. We show that the existence of these chiral edge waves hinges on a distinctive two-step mechanism. Unlike electronic Chern insulators where the magnetic field at the same time gaps the spectrum and imparts chirality, here the finite thickness of the sample gaps the shear modes, and the odd elasticity makes them chiral.

Self-assembled active actomyosin gels spontaneously curve and wrinkle similar to biological cells and tissues

Living systems adopt a diversity of curved and highly dynamic shapes. These diverse morphologies appear on many length scales, from cells to tissues and organismal scales. The common driving force for these dynamic shape changes are contractile stresses generated by myosin motors in the cell cytoskeleton, that converts chemical energy into mechanical work. A good understanding of how contractile stresses in the cytoskeleton arise into different three-dimensional (3D) shapes and what are the shape selection rules that determine their final configurations is still lacking. To obtain insight into the relevant physical mechanisms, we recreate the actomyosin cytoskeleton in vitro, with precisely controlled composition and initial geometry. A set of actomyosin gel discs, intrinsically identical but of variable initial geometry, dynamically self-organize into a family of 3D shapes, such as domes and wrinkled shapes, without the need for specific preprogramming or additional regulation. Shape deformation is driven by the spontaneous emergence of stress gradients driven by myosin and is encoded in the initial disc radius to thickness aspect ratio, which may indicate shaping scalability. Our results suggest that while the dynamical pathways may depend on the detailed interactions between the different microscopic components within the gel, the final selected shapes obey the general theory of elastic deformations of thin sheets. Altogether, our results emphasize the importance for the emergence of active stress gradients for buckling-driven shape deformations and provide insights on the mechanically induced spontaneous shape transitions in contractile active matter, revealing potential shared mechanisms with living systems across scales.

Axis convergence in C. elegans embryos

Embryos develop in a surrounding that guides key aspects of their development. For example, the anteroposterior (AP) body axis is always aligned with the geometric long axis of the surrounding eggshell in fruit flies and worms. The mechanisms that ensure convergence of the AP axis with the long axis of the eggshell remain unresolved. We investigate axis convergence in early C. elegans development, where the nascent AP axis, when misaligned, actively re-aligns to converge with the long axis of the egg. We identify two physical mechanisms that underlie axis convergence. First, bulk cytoplasmic flows, driven by actomyosin cortical flows, can directly reposition the AP axis. Second, active forces generated within the pseudocleavage furrow, a transient actomyosin structure similar to a contractile ring, can drive a mechanical re-orientation such that it becomes positioned perpendicular to the long axis of the egg. This in turn ensures AP axis convergence. Numerical simulations, together with experiments that either abolish the pseudocleavage furrow or change the shape of the egg, demonstrate that the pseudocleavage-furrow-dependent mechanism is a major driver of axis convergence. We conclude that active force generation within the actomyosin cortical layer drives axis convergence in the early nematode.

Shape-Tension Coupling Produces Nematic Order in an Epithelium Vertex Model

We study the vertex model for epithelial tissue mechanics extended to include coupling between the cell shapes and tensions in cell-cell junctions. This coupling represents an active force which drives the system out of equilibrium and leads to the formation of nematic order interspersed with prominent, long-lived +1 defects. The defects in the nematic ordering are coupled to the shape of the cell tiling, affecting cell areas and coordinations. This intricate interplay between cell shape, size, and coordination provides a possible mechanism by which tissues could spontaneously develop long-range polarity through local mechanical forces without resorting to long-range chemical patterning.

Geometry-structure models for liquid crystal interfaces, drops and membranes: wrinkling, shape selection and dissipative shape evolution

We review our recent contributions to anisotropic soft matter models for liquid crystal interfaces, drops and membranes, emphasizing validations with experimental and biological data, and with related theory and simulation literature. The presentation aims to illustrate and characterize the rich output and future opportunities of using a methodology based on the liquid crystal-membrane shape equation applied to static and dynamic pattern formation phenomena. The geometry of static and kinetic shapes is usually described with dimensional curvatures that co-mingle shape and curvedness. In this review, we systematically show how the application of a novel decoupled shape-curvedness framework to practical and ubiquitous soft matter phenomena, such as the shape of drops and tactoids and bending of evolving membranes, leads to deeper quantitative insights than when using traditional dimensional mean and Gaussian curvatures. The review focuses only on (1) statics of wrinkling and shape selection in liquid crystal interfaces and membranes; (2) kinetics and dissipative dynamics of shape evolution in membranes; and (3) computational methods for shape selection and shape evolution; due to various limitations other important topics are excluded. Finally, the outlook follows a similar structure. The main results include: (1) single and multiple wavelength corrugations in liquid crystal interfaces appear naturally in the presence of surface splay and bend orientation distortions with scaling laws governed by ratios of anchoring-to-isotropic tension energy; adding membrane elasticity to liquid crystal anchoring generates multiple scales wrinkling as in tulips; drops of liquid crystals encapsulates in membranes can adopt, according to the ratios of anchoring/tension/bending, families of shapes as multilobal, tactoidal, and serrated as observed in biological cells. (2) Mapping the liquid crystal director to a membrane unit normal. The dissipative shape evolution model with irreversible thermodynamics for flows dominated by bending rates, yields new insights. The model explains the kinetic stability of cylinders, while spheres and saddles are attractors. The model also adds to the evolving understanding of outer hair cells in the inner ear. (3) Computational soft matter geometry includes solving shape equations, trajectories on energy and orientation landscapes, and shape-curvedness evolutions on entropy production landscape with efficient numerical methods and adaptive approaches.

Hexanematic crossover in epithelial monolayers depends on cell adhesion and cell density

Changes in tissue geometry during developmental processes are associated with collective migration of cells. Recent experimental and numerical results suggest that these changes could leverage on the coexistence of nematic and hexatic orientational order at different length scales. How this multiscale organization is affected by the material properties of the cells and their substrate is presently unknown. In this study, we address these questions in monolayers of Madin-Darby canine kidney cells having various cell densities and molecular repertoires. At small length scales, confluent monolayers are characterized by a prominent hexatic order, independent of the presence of E-cadherin, monolayer density, and underlying substrate stiffness. However, all three properties affect the meso-scale tissue organization. The length scale at which hexatic order transits to nematic order, the "hexanematic" crossover scale, strongly depends on cell-cell adhesions and correlates with monolayer density. Our study demonstrates how epithelial organization is affected by mechanical properties, and provides a robust description of tissue organization during developmental processes.

"Magic Numbers" in Self-Faceting of Alcohol-Doped Emulsion Droplets

Oil-in-water emulsion droplets spontaneously adopt, below some temperature Td , counterintuitive faceted and complex non-spherical shapes while remaining liquid. This transition is driven by a crystalline monolayer formed at the droplets' surface. Here, we show that ppm-level doping of the droplet's bulk by long-chain alcohols allows tuning Td by >50 °C, implying formation of drastically different interfacial structures. Furthermore, "magic" alcohol chain lengths maximize Td . This we show to arise from self-assembly of mixed alcohol:alkane interfacial structures of stacked alkane layers, co-crystallized with hydrogen-bonded alcohol dimers. These structures are accounted for theoretically and resolved by direct cryogenic transmission electron microscopy (cryoTEM), confirming the proposed structures. The discovered tunability of key properties of commonly-used emulsions by minute concentrations of specific bulk additives should benefit these emulsions' technological applicability.

Alcanivorax borkumensis biofilms enhance oil degradation by interfacial tubulation

During the consumption of alkanes, Alcanivorax borkumensis will form a biofilm around an oil droplet, but the role this plays during degradation remains unclear. We identified a shift in biofilm morphology that depends on adaptation to oil consumption: Longer exposure leads to the appearance of dendritic biofilms optimized for oil consumption effected through tubulation of the interface. In situ microfluidic tracking enabled us to correlate tubulation to localized defects in the interfacial cell ordering. We demonstrate control over droplet deformation by using confinement to position defects, inducing dimpling in the droplets. We developed a model that elucidates biofilm morphology, linking tubulation to decreased interfacial tension and increased cell hydrophobicity.

Microdomains and stress distributions in bacterial monolayers on curved interfaces

Monolayers of growing non-motile rod-shaped bacteria act as active nematic materials composed of hard particles rather than the flexible components of other commonly studied active nematics. The organization of these granular monolayers has been studied on flat surfaces but not on curved surfaces, which are known to change the behavior of other active nematics. We use molecular dynamics simulations to track alignment and stress in growing monolayers fixed to curved surfaces, and investigate how these vary with changing surface curvature and cell aspect ratio. We find that the length scale of alignment (measured by average microdomain size) increases with cell aspect ratio and decreases with curvature. Additionally, we find that alignment controls the distribution of extensile stresses in the monolayer by concentrating stress in negative-order regions. These results connect active nematic physics to bacterial monolayers and can be applied to model bacteria growing on droplets, such as oil-degrading marine bacteria.

Tuneable defect-curvature coupling and topological transitions in active shells

Recent experimental observations have suggested that topological defects can facilitate the creation of sharp features in developing embryos. Whereas these observations echo established knowledge about the interplay between geometry and topology in two-dimensional passive liquid crystals, the role of activity has mostly remained unexplored. In this article we focus on deformable shells consisting of either polar or nematic active liquid crystals and demonstrate that activity renders the mechanical coupling between defects and curvature much more involved and versatile than previously thought. Using a combination of linear stability analysis and three-dimensional computational fluid dynamics, we demonstrate that such a coupling can in fact be tuned, depending on the type of liquid crystal order, the specific structure of the defect (i.e. asters or vortices) and the nature of the active forces. In polar systems, this can drive a spectacular transition from spherical to toroidal topology, in the presence of large extensile activity. Our analysis strengthens the idea that defects could serve as topological morphogens and provides a number of predictions that could be tested in in vitro studies, for instance in the context of organoids.

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.

Patterning of morphogenetic anisotropy fields

Orientational order, encoded in anisotropic fields, plays an important role during the development of an organism. A striking example of this is the freshwater polyp Hydra, where topological defects in the muscle fiber orientation have been shown to localize to key features of the body plan. This body plan is organized by morphogen concentration gradients, raising the question how muscle fiber orientation, morphogen gradients and body shape interact. Here, we introduce a minimal model that couples nematic orientational order to the gradient of a morphogen field. We show that on a planar surface, alignment to a radial concentration gradient can induce unbinding of topological defects, as observed during budding and tentacle formation in Hydra, and stabilize aster/vortex-like defects, as observed at a Hydra's mouth. On curved surfaces mimicking the morphologies of Hydra in various stages of development-from spheroid to adult-our model reproduces the experimentally observed reorganization of orientational order. Our results suggest how gradient alignment and curvature effects may work together to control orientational order during development and lay the foundations for future modeling efforts that will include the tissue mechanics that drive shape deformations.

Topological phases and curvature-driven pattern formation in cholesteric shells

We study the phase behaviour of cholesteric liquid crystal shells with different geometries. We compare the cases of tangential anchoring and no anchoring at the surface, focussing on the former case, which leads to a competition between the intrinsic tendency of the cholesteric to twist and the anchoring free energy which suppresses it. We then characterise the topological phases arising close to the isotropic-cholesteric transition. These typically consist of quasi-crystalline or amorphous tessellations of the surface by half-skyrmions, which are stable at lower and larger shell sizes, respectively. For ellipsoidal shells, defects in the tessellation couple to a local curvature, and according to the shell size, they either migrate to the poles or distribute uniformly on the surface. For toroidal shells, the variations in the local curvature of the surface stabilise heterogeneous phases where cholesteric or isotropic patterns coexist with hexagonal lattices of half-skyrmions.

Hydrodynamic Enhancement of p-atic Defect Dynamics

We investigate numerically and analytically the effects of hydrodynamics on the dynamics of topological defects in p-atic liquid crystals, i.e., two-dimensional liquid crystals with p-fold rotational symmetry. Importantly, we find that hydrodynamics fuels a generic passive self-propulsion mechanism for defects of winding number s=(p-1)/p and arbitrary p. Strikingly, we discover that hydrodynamics always accelerates the annihilation dynamics of pairs of ±1/p defects and that, contrary to expectations, this effect increases with p. Our Letter paves the way toward understanding cell intercalation and other remodeling events in epithelial layers.

Alignment rule and geometric confinement lead to stability of a vortex in active flow

Vortices are hallmarks of a wide range of nonequilibrium phenomena in fluids at multiple length scales. In this work, we numerically study the whirling motion of self-propelled soft point particles confined in circular domain, and aim at addressing the stability issue of the coherent vortex structure. By the combination of dynamical and statistical analysis at the individual particle level, we reveal the persistence of the whirling motion resulting from the subtle competition of activity and geometric confinement. In the stable whirling motion, the scenario of the coexistence of the irregular microscopic motions of individual particles and the regular global whirling motion is fundamentally different from the motion of a vortex in passive fluid. Possible orientational order coexisting with the whirling are further explored. This work shows the stability mechanism of vortical dynamics in active media under the alignment rule in confined space and may have implications in creating and harnessing macroscale coherent dynamical states by tuning the confining geometry.

Topological defect-mediated morphodynamics of active-active interfaces

Physical interfaces widely exist in nature and engineering. Although the formation of passive interfaces is well elucidated, the physical principles governing active interfaces remain largely unknown. Here, we combine simulation, theory, and cell-based experiment to investigate the evolution of an active-active interface. We adopt a biphasic framework of active nematic liquid crystals. We find that long-lived topological defects mechanically energized by activity display unanticipated dynamics nearby the interface, where defects perform "U-turns" to keep away from the interface, push the interface to develop local fingers, or penetrate the interface to enter the opposite phase, driving interfacial morphogenesis and cross-interface defect transport. We identify that the emergent interfacial morphodynamics stems from the instability of the interface and is further driven by the activity-dependent defect-interface interactions. Experiments of interacting multicellular monolayers with extensile and contractile differences in cell activity have confirmed our predictions. These findings reveal a crucial role of topological defects in active-active interfaces during, for example, boundary formation and tissue competition that underlie organogenesis and clinically relevant disorders.

Tubulation and dispersion of oil by bacterial growth on droplets

Bacteria on surfaces exhibit collective behaviors, such as active turbulence and active stresses, which result from their motility, growth, and interactions with their local surroundings. However, interfacial deformations on soft surfaces and liquid interfaces caused by active growth, particularly over long time scales, are not well understood. Here, we describe experimental observations on the emergence of tubular structures arising from the growth of rod-shaped bacteria at the interface of oil droplets in water. Using microfluidics and timelapse microscopy, the dimensions and extension rates of individual tubular structures as well as bulk bio-aggregate formation are quantified for hundreds of droplets over 72 hours. Tubular structures are comparable in length to the initial droplet radius and are composed of an outer shell of bacteria that stabilize an inner filament of oil. The oil filament breaks up into smaller microdroplets dispersed within the bacterial shell. This work provides insight into active stresses at deformable interfaces and improves our understanding of microbial oil biodegradation and its potential influence on the transport of droplets in the ocean water column.

The emergence of spontaneous coordinated epithelial rotation on cylindrical curved surfaces

Three-dimensional collective epithelial rotation around a given axis represents a coordinated cellular movement driving tissue morphogenesis and transformation. Questions regarding these behaviors and their relationship with substrate curvatures are intimately linked to spontaneous active matter processes and to vital morphogenetic and embryonic processes. Here, using interdisciplinary approaches, we study the dynamics of epithelial layers lining different cylindrical surfaces. We observe large-scale, persistent, and circumferential rotation in both concavely and convexly curved cylindrical tissues. While epithelia of inverse curvature show an orthogonal switch in actomyosin network orientation and opposite apicobasal polarities, their rotational movements emerge and vary similarly within a common curvature window. We further reveal that this persisting rotation requires stable cell-cell adhesion and Rac-1-dependent cell polarity. Using an active polar gel model, we unveil the different relationships of collective cell polarity and actin alignment with curvatures, which lead to coordinated rotational behavior despite the inverted curvature and cytoskeleton order.

Active Nematic Defects and Epithelial Morphogenesis

Inspired by recent experiments that highlight the role of nematic defects in the morphogenesis of epithelial tissues, we develop a minimal framework to study the dynamics of an active curved surface driven by its nematic texture. Allowing the surface to evolve via relaxational dynamics leads to a theory linking nematic defect dynamics, cellular division rates, and Gaussian curvature. Regions of large positive (negative) curvature and positive (negative) growth are colocalized with the presence of positive (negative) defects. In an ex-vivo setting of cultured murine neural progenitor cells, we show that our framework is consistent with the observed cell accumulation at positive defects and depletion at negative defects. In an in-vivo setting, we show that the defect configuration consisting of a bound +1 defect state, which is stabilized by activity, surrounded by two -1/2 defects can create a stationary ring configuration of tentacles, consistent with observations of a basal marine invertebrate Hydra.

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.

Spontaneous flow created by active topological defects

Topological defects are at the root of the large-scale organization of liquid crystals. In two-dimensional active nematics, two classes of topological defects of charges [Formula: see text] are known to play a major role due to active stresses. Despite this importance, few analytical results have been obtained on the flow-field and active-stress patterns around active topological defects. Using the generic hydrodynamic theory of active systems, we investigate the flow and stress patterns around these topological defects in unbounded, two-dimensional active nematics. Under generic assumptions, we derive analytically the spontaneous velocity and stall force of self-advected defects in the presence of both shear and rotational viscosities. Applying our formalism to the dynamics of monolayers of elongated cells at confluence, we show that the non-conservation of cell number generically increases the self-advection velocity and could provide an explanation for their observed role in cellular extrusion and multilayering. We finally investigate numerically the influence of the Ericksen stress. Our work paves the way to a generic study of the role of topological defects in active nematics, and in particular in monolayers of elongated cells.