Morphogenesis driven by nematic defects in active biological networks

Technical Summary: Morphogenesis Driven by Nematic Defects in Active Biological Networks

Problem Statement
The paper addresses the fundamental challenge of modeling cellular morphogenesis, specifically how biological tissues transition from disordered states to ordered structures with specific shapes. While it is established that cytoskeletal filaments can organize into nematic phases with partial orientational order, the precise mechanical role of topological defects within these phases in driving tissue deformation and pattern formation remains unclear. Existing models often rely on linear approximations or pre-defined morphogen gradients, which may fail to capture the complex, non-linear interplay between nematic order, elasticity, and stress-driven growth observed in regenerating organisms like Hydra. The authors aim to develop a continuum framework that treats living tissues as active biological networks (ABNs) where shape changes emerge from the coupling of nematic order, elasticity, and adaptive relaxation, without prescribing the growth pattern a priori.

Methodology
The authors propose a quasi-static continuum framework modeling living tissues as active biological networks composed of cross-linked filament chains capable of forming a nematic phase. The core assumptions and mathematical formulation include:

• Kinematics and Growth: The tissue is modeled as locally incompressible. Morphogenesis is described as a sequence of equilibrium states. Growth is not prescribed but is driven by local mechanical stress. The deformation gradient $\mathbf{F}$ is multiplicatively decomposed into a nematic step-length tensor $\mathbf{L}$, an elastic distortion $\mathbf{F}_e$, a growth tensor $\mathbf{F}_g$, and a relaxation tensor $\mathbf{F}_{rel}$.
• Energetic Formulation: The total free energy $\mathcal{F}$ comprises two main contributions:
  1. Nematic Energy ($\mathcal{F}_n$): Based on Ericksen's continuum theory, this term accounts for spatial variations in the director field $\mathbf{n}$ and the scalar order parameter $s$. Crucially, $s$ is allowed to vary spatially, vanishing at topological defects to regularize singularities and allow for local transitions to an isotropic phase.
  2. Elastic Energy ($\mathcal{F}_e$): A neo-classical strain-energy function extended to include growth and relaxation. It depends on the elastic part of the deformation, capturing the response of the fiber network to distortions.
• Coupling Mechanisms: The model integrates four mechanisms: (a) mechanical stress near topological defects, (b) activity-induced ordering (increasing $s$), (c) stress-driven growth (where new material is added proportional to elastic energy density), and (d) structural relaxation (updating the stress-free configuration to relieve internal stresses).
• Numerical Implementation: The coupled nonlinear problem is solved using a mixed finite element formulation implemented in FEniCSx. To handle the head-tail symmetry of nematic directors and the associated branch cuts for half-integer defects, the director field is reformulated using a symmetric rank-2 tensor field $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$. The simulations employ a Newton method with line search and incremental loading to track the evolution from an initial isotropic spheroid to mature morphologies.

Key Contributions

1. Continuum Framework for ABNs: The paper introduces a mathematical model that couples nematic order, elasticity, and stress-driven growth/relaxation, moving beyond linear shell theories to capture non-linear morphogenetic outcomes.
2. Defect-Driven Morphogenesis: The study demonstrates that the topology of nematic defects (specifically their charge and spatial arrangement) acts as the primary driver for shape evolution. The model shows that defects are not merely singularities but organizing centers that localize stress and direct growth.
3. Reproduction of Complex Morphologies: By varying the initial defect configuration, the model successfully reproduces three distinct morphogenetic outcomes observed in Hydra regeneration:
   • Uniaxial (Head-Foot) Regeneration: Two $+1$ defects at opposite poles drive elongation and the formation of a mouth (apical rupture) and a foot.
   • Bicephalic (Two-Head) Development: Two $+1$ defects in the upper hemisphere, balanced by $-1/2$ defects, lead to the formation of two distinct heads.
   • Tentacle Formation: A configuration with $+1$ defects at tentacle tips and $-1/2$ defects at their bases drives the outgrowth of tentacles.
4. Stress-Driven Growth Mechanism: The simulations explicitly show that growth is localized to regions of high elastic stress (near defects), while relaxation allows for the enlargement of apertures (mouths) where stress can be geometrically relieved, contrasting with persistent stress at the foot or tentacle tips.

Results
Finite element simulations of Hydra regeneration reveal that:

• Stress Localization: Distortions in the nematic director field near topological defects generate high mechanical stresses. These stresses drive the accumulation of new material (growth) in their vicinity.
• Asymmetric Outcomes: In the uniaxial configuration, stress relaxation allows the apical pole (mouth) to rupture and enlarge, while the basal pole (foot) remains closed and structurally stressed.
• Defect Roles: $+1$ defects drive protrusion formation and tissue rupture. $-1/2$ defects act primarily as structural stabilizers or "topological balancers" in early stages, though in the tentacular configuration, stress accumulation at these defects eventually triggers localized growth at the tentacle bases.
• Morphological Evolution: The model captures the transition from a spherical, isotropic tissue fragment to complex, elongated structures with specific features (oral pores, tentacles) solely through the evolution of the nematic order parameter and the resulting stress-driven growth.

Significance and Claims
The authors claim that their work provides a robust theoretical foundation for understanding how topological defects in active nematic networks can drive complex morphogenetic processes without the need for pre-specified growth patterns or external morphogen gradients. By treating the tissue as an active biological network where growth is a direct response to mechanical stress generated by nematic distortions, the model bridges the gap between microscopic cytoskeletal organization and macroscopic tissue shape. The results align with experimental observations in Hydra, suggesting that the "memory" of body axis polarity and the emergence of specific structures are encoded in the topological charge and arrangement of nematic defects. The paper positions this framework as a step toward understanding the intrinsic mechanical drivers of pattern formation in living tissues, distinct from purely biochemical signaling models.