Thermally Driven Order-Disorder Transition in Two-Dimensional Soft Cellular Systems

Cell Division and Motility Enable Hexatic Order in Biological Tissues

Biological tissues transform between solid- and liquidlike states in many fundamental physiological events. Recent experimental observations further suggest that in two-dimensional epithelial tissues these solid-liquid transformations can happen via intermediate states akin to the intermediate hexatic phases observed in equilibrium two-dimensional melting. The hexatic phase is characterized by quasi-long-range (power-law) orientational order but no translational order, thus endowing some structure to an otherwise structureless fluid. While it has been shown that hexatic order in tissue models can be induced by motility and thermal fluctuations, the role of cell division and apoptosis (birth and death) has remained poorly understood, despite its fundamental biological role. Here we study the effect of cell division and apoptosis on global hexatic order within the framework of the self-propelled Voronoi model of tissue. Although cell division naively destroys order and active motility facilitates deformations, we show that their combined action drives a liquid-hexatic-liquid transformation as the motility increases. The hexatic phase is accessed by the delicate balance of dislocation defect generation from cell division and the active binding of disclination-antidisclination pairs from motility. We formulate a mean-field model to elucidate this competition between cell division and motility and the consequent development of hexatic order.

Hydrodynamics and multiscale order in confluent epithelia

We formulate a hydrodynamic theory of confluent epithelia: i.e. monolayers of epithelial cells adhering to each other without gaps. Taking advantage of recent progresses toward establishing a general hydrodynamic theory of p-atic liquid crystals, we demonstrate that collectively migrating epithelia feature both nematic (i.e. p = 2) and hexatic (i.e. p = 6) orders, with the former being dominant at large and the latter at small length scales. Such a remarkable multiscale liquid crystal order leaves a distinct signature in the system's structure factor, which exhibits two different power-law scaling regimes, reflecting both the hexagonal geometry of small cells clusters and the uniaxial structure of the global cellular flow. We support these analytical predictions with two different cell-resolved models of epithelia - i.e. the self-propelled Voronoi model and the multiphase field model - and highlight how momentum dissipation and noise influence the range of fluctuations at small length scales, thereby affecting the degree of cooperativity between cells. Our construction provides a theoretical framework to conceptualize the recent observation of multiscale order in layers of Madin-Darby canine kidney cells and pave the way for further theoretical developments.

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.

Active cell divisions generate fourfold orientationally ordered phase in living tissue

Morphogenesis, the process through which genes generate form, establishes tissue-scale order as a template for constructing the complex shapes of the body plan. The extensive growth required to build these ordered substrates is fuelled by cell proliferation, which, naively, should destroy order. Understanding how active morphogenetic mechanisms couple cellular and mechanical processes to generate order-rather than annihilate it-remains an outstanding question in animal development. We show that cell divisions are the primary drivers of tissue flow, leading to a fourfold orientationally ordered phase. Waves of anisotropic cell proliferation propagate across the embryo with precise patterning. Defects introduced into the nascent lattice by cell divisions are moved out of the tissue bulk towards the boundary by subsequent divisions. Specific cell proliferation rates and orientations enable cell divisions to organize rather than fluidize the tissue. We observe this using live imaging and tissue cartography to analyse the dynamics of fourfold tissue ordering in the trunk segmental ectoderm of the crustacean Parhyale hawaiensis beginning 72 h after egg lay. The result is a robust, active mechanism for generating global orientational order in a non-equilibrium system that sets the stage for the subsequent development of shape and form.

Micro and Macroscopic Stress-Strain Relations in Disordered Tessellated Networks

We demonstrate that for a rigid and incompressible network in mechanical equilibrium, the microscopic stress and strain follows a simple relation, σ=pE, where σ is the deviatoric stress, E is a mean-field strain tensor, and p is the hydrostatic pressure. This relationship arises as the natural consequence of energy minimization or equivalently, mechanical equilibration. The result suggests not only that the microscopic stress and strain are aligned in the principal directions, but also microscopic deformations are predominantly affine. The relationship holds true regardless of the different (foam or tissue) energy model considered, and directly leads to a simple prediction for the shear modulus, μ=⟨p⟩/2, where ⟨p⟩ is the mean pressure of the tessellation, for general randomized lattices.

Local Yield and Compliance in Active Cell Monolayers

The rheology of biological tissue plays an important role in many processes, from organ formation to cancer invasion. Here, we use a multiphase field model of motile cells to simulate active microrheology within a tissue monolayer. When unperturbed, the tissue exhibits a transition between a solidlike state and a fluidlike state tuned by cell motility and deformability-the ratio of the energetic costs of steric cell-cell repulsion and cell-edge tension. When perturbed, solid tissues exhibit local yield-stress behavior, with a threshold force for the onset of motion of a probe particle that vanishes upon approaching the solid-to-liquid transition. This onset of motion is qualitatively different in the low and high deformability regimes. At high deformability, the tissue is amorphous when solid, it responds compliantly to deformations, and the probe transition to motion is smooth. At low deformability, the monolayer is more ordered translationally and stiffer, and the onset of motion appears discontinuous. Our results suggest that cellular or nanoparticle transport in different types of tissues can be fundamentally different and point to ways in which it can be controlled.

Active T1 transitions in cellular networks

In amorphous solids as in tissues, neighbor exchanges can relax local stresses and allow the material to flow. In this paper, we use an anisotropic vertex model to study T1 rearrangements in polygonal cellular networks. We consider two different physical realizations of the active anisotropic stresses: (i) anisotropic bond tension and (ii) anisotropic cell stress. Interestingly, the two types of active stress lead to patterns of relative orientation of T1 transitions and cell elongation that are different. Our work suggests that these two realizations of anisotropic active stresses can be observed in vivo. We describe and explain these results through the lens of a continuum description of the tissue as an anisotropic active material. We furthermore discuss the energetics of the dynamic tissue and express the energy balance in terms of internal elastic energy, mechanical work, chemical work and heat. This allows us to define active T1 transitions that can perform mechanical work while consuming chemical energy.

Topographic Control of Order in Quasi-2D Granular Phase Transitions

We experimentally investigate the nature of 2D phase transitions in a quasi-2D granular fluid. Using a surface decorated with periodically spaced dimples we observe interfacial tension between coexisting granular liquid and crystal phases. Measurements of the orientational and translational order parameters and associated susceptibilities indicate that the surface topography alters the order of the phase transition from a two-step continuous one to a first-order liquid-solid one. The interplay of boundary inelasticity and geometry, either order promoting or inhibiting, controls whether it is the granular crystal or the granular fluid which makes contact with the edge. This order induced wetting has important consequences, determining how coexisting phases separate spatially.

Emergence and melting of active vortex crystals

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

Large-scale simulations of biological cell sorting driven by differential adhesion follow diffusion-limited domain coalescence regime

Cell sorting, whereby a heterogeneous cell mixture segregates and forms distinct homogeneous tissues, is one of the main collective cell behaviors at work during development. Although differences in interfacial energies are recognized to be a possible driving source for cell sorting, no clear consensus has emerged on the kinetic law of cell sorting driven by differential adhesion. Using a modified Cellular Potts Model algorithm that allows for efficient simulations while preserving the connectivity of cells, we numerically explore cell-sorting dynamics over very large scales in space and time. For a binary mixture of cells surrounded by a medium, increase of domain size follows a power-law with exponent n = 1/4 independently of the mixture ratio, revealing that the kinetics is dominated by the diffusion and coalescence of rounded domains. We compare these results with recent numerical studies on cell sorting, and discuss the importance of algorithmic differences as well as boundary conditions on the observed scaling.

Theory and simulation for equilibrium glassy dynamics in cellular Potts model of confluent biological tissue

Glassy dynamics in a confluent monolayer is indispensable in morphogenesis, wound healing, bronchial asthma, and many others; a detailed theoretical framework for such a system is, therefore, important. Vertex-model (VM) simulations have provided crucial insights into the dynamics of such systems, but their nonequilibrium nature makes theoretical development difficult. The cellular Potts model (CPM) of confluent monolayers provides an alternative model for such systems with a well-defined equilibrium limit. We combine numerical simulations of the CPM and an analytical study based on one of the most successful theories of equilibrium glass, the random first-order transition theory, and develop a comprehensive theoretical framework for a confluent glassy system. We find that the glassy dynamics within the CPM is qualitatively similar to that in the VM. Our study elucidates the crucial role of geometric constraints in bringing about two distinct regimes in the dynamics, as the target perimeter P_{0} is varied. The unusual sub-Arrhenius relaxation results from the distinctive interaction potential arising from the perimeter constraint in such systems. The fragility of the system decreases with increasing P_{0} in the low-P_{0} regime, whereas the dynamics is independent of P_{0} in the other regime. The rigidity transition, found in the VM, is absent within the CPM; this difference seems to come from the nonequilibrium nature of the former. We show that the CPM captures the basic phenomenology of glassy dynamics in a confluent biological system via comparison of our numerical results with existing experiments on different systems.

Thermal shape fluctuations of a two-dimensional compressible droplet

Analysis of thermal capillary waves on the surface of a liquid usually assumes incompressibility of the bulk fluid. However, for droplets or bubbles with submicronic size, or for epithelial cells whose out-of-plane elongation can be modeled by an effective 2D bulk modulus, compressibility of the internal fluid must be taken into account for the characterization of their shape fluctuations. We present a theoretical analysis of the fluctuations of a two-dimensional compressible droplet. Analytical expressions for area, perimeter and energy fluctuations are derived and compared with Cellular Potts Model (CPM) simulations. This comparison shows a very good agreement between theory and simulations, and offers a precise calibration method for CPM simulations.

Solid-Liquid Transition of Deformable and Overlapping Active Particles

Experiments and theory have shown that cell monolayers and epithelial tissues exhibit solid-liquid and glass-liquid transitions. These transitions are biologically relevant to our understanding of embryonic development, wound healing, and cancer. Current models of confluent epithelia have focused on the role of cell shape, with less attention paid to cell extrusion, which is key for maintaining homeostasis in biological tissue. Here, we use a multiphase field model to study the solid-liquid transition in a confluent monolayer of deformable cells. Cell overlap is allowed and provides a way for modeling the precursor for extrusion. When cells overlap rather than deform, we find that the melting transition changes from continuous to first order like, and that there is an intermittent regime close to the transition, where solid and liquid states alternate over time. By studying the dynamics of five- and sevenfold disclinations in the hexagonal lattice formed by the cell centers, we observe that these correlate with spatial fluctuations in the cellular overlap, and that cell extrusion tends to initiate near fivefold disclinations.

Hexatic phase in a model of active biological tissues

In many biological processes, such as wound healing, cell tissues undergo an epithelial-to-mesenchymal transition, which is a transition from a more rigid to a more fluid state. Here, we investigate the solid/fluid transition of cell tissues within the framework of the self-propelled Voronoi model, which accounts for the deformability of the cells, for their many-body interactions, and for their polarized motility. The transition is controlled by two parameters, respectively accounting for the strength of the self-propelling force of the cells, and for the mechanical rigidity of the cells. We find the melting transition to occur via a continuous solid-hexatic transition followed by a continuous hexatic-liquid transition, as in the Kosterlitz, Thouless, Halperin, Nelson, and Young scenario. This finding indicates that the hexatic phase may have an unexpected biological relevance.