Geometric Frustration and Solid-Solid Transitions in Model 2D Tissue

Regulation of chromatin modifications through coordination of nucleus size and epithelial cell morphology heterogeneity

Cell morphology heterogeneity is pervasive in epithelial collectives, yet the underlying mechanisms driving such heterogeneity and its consequential biological ramifications remain elusive. Here, we observed a consistent correlation between the epithelial cell morphology and nucleus morphology during crowding, revealing a persistent log-normal probability distribution characterizing both cell and nucleus areas across diverse epithelial model systems. We showed that this morphological diversity arises from asymmetric partitioning during cell division. Next, we provide insights into the impact of nucleus morphology on chromatin modifications. We demonstrated that constraining nucleus leads to downregulation of the euchromatic mark H3K9ac and upregulation of the heterochromatic mark H3K27me3. Furthermore, we showed that nucleus size regulates H3K27me3 levels through histone demethylase UTX. These findings highlight the significance of cell morphology heterogeneity as a driver of chromatin state diversity, shaping functional variability within epithelial tissues.

Regulation of Chromatin Modifications through Coordination of Nucleus Size and Epithelial Cell Morphology Heterogeneity

Cell morphology heterogeneity is pervasive in epithelial collectives, yet the underlying mechanisms driving such heterogeneity and its consequential biological ramifications remain elusive. Here, we observed a consistent correlation between the epithelial cell morphology and nucleus morphology during crowding, revealing a persistent log-normal probability distribution characterizing both cell and nucleus areas across diverse epithelial model systems. We further showed that this morphological diversity arises from asymmetric partitioning during cell division. Moreover, we provide insights into the impact of nucleus morphology on chromatin modifications. We demonstrated that constraining nucleus leads to downregulation of the euchromatic mark H3K9ac and upregulation of the heterochromatic mark H3K27me3. Furthermore, we showed that nucleus size regulates H3K27me3 levels through histone demethylase UTX. These findings highlight the significance of cell morphology heterogeneity as a driver of chromatin state diversity, shaping functional variability within epithelial tissues.

Generic Elasticity of Thermal, Underconstrained Systems

Athermal (i.e., zero-temperature) underconstrained systems are typically floppy, but they can be rigidified by the application of external strain, which is theoretically well understood. Here and in the companion paper [C. T. Lee and M. Merkel, Phys. Rev. E 110, 064147 (2024)PRESCM2470-004510.1103/PhysRevE.110.064147], we extend this theory to finite temperatures for a very broad class of underconstrained systems. In the vicinity of the athermal transition point, we derive from first principles expressions for elastic properties such as isotropic tension t and shear modulus G on temperature T, isotropic strain ϵ, and shear strain γ, which we confirm numerically. These expressions contain only three parameters: entropic rigidity κ_{S}, energetic rigidity κ_{E}, and a parameter b_{ϵ} describing the interaction between isotropic and shear strain, which can be determined from the microstructure of the system. Our results imply that in underconstrained systems, entropic and energetic rigidity interact like two springs in series. This also allows for a simple explanation of the previously numerically observed scaling relation t∼G∼T^{1/2} at ϵ=γ=0. Our work unifies the physics of systems as diverse as polymer fibers and networks, membranes, and vertex models for biological tissues.

Finite elasticity of the vertex model and its role in rigidity of curved cellular tissues

Using a mean field approach and simulations, we study the non-linear mechanical response of the vertex model (VM) of biological tissue to compression and dilation. The VM is known to exhibit a transition between solid and fluid-like, or floppy, states driven by geometric incompatibility. Target perimeter and area set a target shape which may not be geometrically achievable, thereby engendering frustration. Previously, an asymmetry in the linear elastic response was identified at the rigidity transition between compression and dilation. Here we show that the asymmetry extends away from the transition point for finite strains. Under finite compression, an initially solid VM can completely relax perimeter tension, resulting in a drop discontinuity in the mechanical response. Conversely, an initially floppy VM under dilation can rigidify and have a higher response. These observations imply that re-scaling of cell area shifts the transition between rigid and floppy states. Based on this insight, we calculate the re-scaling of cell area engendered by intrinsic curvature and write a prediction for the rigidity transition in the presence of curvature. The shift of the rigidity transition in the presence of curvature for the VM provides a new metric for predicting tissue rigidity from image data of curved tissues in a manner analogous to the flat case.

Constitutive model for the rheology of biological tissue

The rheology of biological tissue is key to processes such as embryo development, wound healing, and cancer metastasis. Vertex models of confluent tissue monolayers have uncovered a spontaneous liquid-solid transition tuned by cell shape; and a shear-induced solidification transition of an initially liquidlike tissue. Alongside this jamming/unjamming behavior, biological tissue also displays an inherent viscoelasticity, with a slow time and rate-dependent mechanics. With this motivation, we combine simulations and continuum theory to examine the rheology of the vertex model in nonlinear shear across a full range of shear rates from quastistatic to fast, elucidating its nonlinear stress-strain curves after the inception of shear of finite rate, and its steady state flow curves of stress as a function of strain rate. We formulate a rheological constitutive model that couples cell shape to flow and captures both the tissue solid-liquid transition and its rich linear and nonlinear rheology.

The role of non-affine deformations in the elastic behavior of the cellular vertex model

The vertex model of epithelia describes the apical surface of a tissue as a tiling of polygonal cells, with a mechanical energy governed by deviations in cell shape from preferred, or target, area, A₀, and perimeter, P₀. The model exhibits a rigidity transition driven by geometric incompatibility as tuned by the target shape index, . For with p_*(6) the perimeter of a regular hexagon of unit area, a cell can simultaneously attain both the preferred area and preferred perimeter. As a result, the tissue is in a mechanically soft compatible state, with zero shear and Young's moduli. For p₀ < p_*(6), it is geometrically impossible for any cell to realize the preferred area and perimeter simultaneously, and the tissue is in an incompatible rigid solid state. Using a mean-field approach, we present a complete analytical calculation of the linear elastic moduli of an ordered vertex model. We analyze a relaxation step that includes non-affine deformations, leading to a softer response than previously reported. The origin of the vanishing shear and Young's moduli in the compatible state is the presence of zero-energy deformations of cell shape. The bulk modulus exhibits a jump discontinuity at the transition and can be lower in the rigid state than in the fluid-like state. The Poisson's ratio can become negative which lowers the bulk and Young's moduli. Our work provides a unified treatment of linear elasticity for the vertex model and demonstrates that this linear response is protocol-dependent.

Geometric theory of mechanical screening in two-dimensional solids

Holes in mechanical metamaterials, quasilocalized plastic events in amorphous solids, and bound dislocations in a hexatic matter are different mechanisms of generic stress relaxation in solids. Regardless of the specific mechanism, these and other local stress relaxation modes are quadrupolar in nature, forming the foundation for stress screening in solids, similar to polarization fields in electrostatic media. We propose a geometric theory for stress screening in generalized solids based on this observation. The theory includes a hierarchy of screening modes, each characterized by internal length scales, and is partially analogous to theories of electrostatic screening such as dielectrics and Debye-Hückel theory. Additionally, our formalism suggests that the hexatic phase, traditionally defined by structural properties, can also be defined by mechanical properties and may exist in amorphous materials.

Truchet-tile structure of a topologically aperiodic metal-organic framework

When tiles decorated to lower their symmetry are joined together, they can form aperiodic and labyrinthine patterns. Such Truchet tilings offer an efficient mechanism of visual data storage related to that used in barcodes and QR codes. We show that the crystalline metal-organic framework [OZn₄][1,3-benzenedicarboxylate]₃ (TRUMOF-1) is an atomic-scale realization of a complex three-dimensional Truchet tiling. Its crystal structure consists of a periodically arranged assembly of identical zinc-containing clusters connected uniformly in a well-defined but disordered fashion to give a topologically aperiodic microporous network. We suggest that this unusual structure emerges as a consequence of geometric frustration in the chemical building units from which it is assembled.

Instabilities and Geometry of Growing Tissues

We present a covariant continuum formulation of a generalized two-dimensional vertexlike model of epithelial tissues which describes tissues with different underlying geometries, and allows for an analytical macroscopic description. Using a geometrical approach and out-of-equilibrium statistical mechanics, we calculate both mechanical and dynamical instabilities of a tissue, and their dependences on various variables, including activity, and cell-shape heterogeneity (disorder). We show how both plastic cellular rearrangements and the tissue elastic response depend on the existence of mechanical residual stresses at the cellular level. Even freely growing tissues may exhibit a growth instability depending on the intrinsic proliferation rate. Our main result is an explicit calculation of the cell pressure in a homeostatic state of a confined growing tissue. We show that the homeostatic pressure can be negative and depends on the existence of mechanical residual stresses. This geometric model allows us to sort out elastic and plastic effects in a growing, flowing, tissue.

Anomalous elasticity of a cellular tissue vertex model

Vertex models, such as those used to describe cellular tissue, have an energy controlled by deviations of each cell area and perimeter from target values. The constrained nonlinear relation between area and perimeter leads to new mechanical response. Here we provide a mean-field treatment of a highly simplified model: a uniform network of regular polygons with no topological rearrangements. Since all polygons deform in the same way, we only need to analyze the ground states and the response to deformations of a single polygon (cell). The model exhibits the known transition between a fluid/compatible state, where the cell can accommodate both target area and perimeter, and a rigid/incompatible state. We calculate and measure the mechanical resistance to various deformation protocols and discover that at the onset of rigidity, where a single zero-energy ground state exists, linear elasticity fails to describe the mechanical response to even infinitesimal deformations. In particular, we identify a breakdown of reciprocity expressed via different moduli for compressive and tensile loads, implying nonanalyticity of the energy functional. We give a pictorial representation in configuration space that reveals that the complex elastic response of the vertex model arises from the presence of two distinct sets of reference states (associated with target area and target perimeter). Our results on the critically compatible tissue provide a new route for the design of mechanical metamaterials that violate or extend classical elasticity.

Linear viscoelastic properties of the vertex model for epithelial tissues

Epithelial tissues act as barriers and, therefore, must repair themselves, respond to environmental changes and grow without compromising their integrity. Consequently, they exhibit complex viscoelastic rheological behavior where constituent cells actively tune their mechanical properties to change the overall response of the tissue, e.g., from solid-like to fluid-like. Mesoscopic mechanical properties of epithelia are commonly modeled with the vertex model. While previous studies have predominantly focused on the rheological properties of the vertex model at long time scales, we systematically studied the full dynamic range by applying small oscillatory shear and bulk deformations in both solid-like and fluid-like phases for regular hexagonal and disordered cell configurations. We found that the shear and bulk responses in the fluid and solid phases can be described by standard spring-dashpot viscoelastic models. Furthermore, the solid-fluid transition can be tuned by applying pre-deformation to the system. Our study provides insights into the mechanisms by which epithelia can regulate their rich rheological behavior.

Shear-Driven Solidification and Nonlinear Elasticity in Epithelial Tissues

Biological processes, from morphogenesis to tumor invasion, spontaneously generate shear stresses inside living tissue. The mechanisms that govern the transmission of mechanical forces in epithelia and the collective response of the tissue to bulk shear deformations remain, however, poorly understood. Using a minimal cell-based computational model, we investigate the constitutive relation of confluent tissues under simple shear deformation. We show that an initially undeformed fluidlike tissue acquires finite rigidity above a critical applied strain. This is akin to the shear-driven rigidity observed in other soft matter systems. Interestingly, shear-driven rigidity can be understood by a critical scaling analysis in the vicinity of the second order critical point that governs the liquid-solid transition of the undeformed system. We further show that a solidlike tissue responds linearly only to small strains and but then switches to a nonlinear response at larger stains, with substantial stiffening. Finally, we propose a mean-field formulation for cells under shear that offers a simple physical explanation of shear-driven rigidity and nonlinear response in a tissue.

Energetic rigidity. II. Applications in examples of biological and underconstrained materials

This is the second paper devoted to energetic rigidity, in which we apply our formalism to examples in two dimensions: Underconstrained random regular spring networks, vertex models, and jammed packings of soft particles. Spring networks and vertex models are both highly underconstrained, and first-order constraint counting does not predict their rigidity, but second-order rigidity does. In contrast, spherical jammed packings are overconstrained and thus first-order rigid, meaning that constraint counting is equivalent to energetic rigidity as long as prestresses in the system are sufficiently small. Aspherical jammed packings on the other hand have been shown to be jammed at hypostaticity, which we use to argue for a modified constraint counting for systems that are energetically rigid at quartic order.

Jamming and arrest of cell motion in biological tissues

Collective cell motility is crucial to many biological processes including morphogenesis, wound healing, and cancer invasion. Recently, the biology and biophysics communities have begun to use the term 'cell jamming' to describe the collective arrest of cell motion in tissues. Although this term is widely used, the underlying mechanisms are varied. In this review, we highlight three independent mechanisms that can potentially drive arrest of cell motion - crowding, tension-driven rigidity, and reduction of fluctuations - and propose a framework that connects all three. Because multiple mechanisms may be operating simultaneously, this emphasizes that experiments should strive to identify which mechanism dominates in a given situation. We also discuss how specific cell-scale and molecular-scale biological processes, such as cell-cell and cell-substrate interactions, control aspects of these underlying physical mechanisms.

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.

Peristalsis by pulses of activity

Peristalsis by actively generated waves of muscle contraction is one of the most fundamental ways of producing motion in living systems. We show that peristalsis can be modeled by a train of rectangular-shaped solitary waves of localized activity propagating through otherwise passive matter. Our analysis is based on the Fermi-Pasta-Ulam (FPU) type discrete model accounting for active stresses and we reveal the existence in this problem of a critical regime which we argue to be physiologically advantageous.

Softness, anomalous dynamics, and fractal-like energy landscape in model cell tissues

Epithelial cell tissues have a slow relaxation dynamics resembling that of supercooled liquids. Yet, they also have distinguishing features. These include an extended short-time subdiffusive transient, as observed in some experiments and recent studies of model systems, and a sub-Arrhenius dependence of the relaxation time on temperature, as reported in numerical studies. Here we demonstrate that the anomalous glassy dynamics of epithelial tissues originates from the emergence of a fractal-like energy landscape, particles becoming virtually free to diffuse in specific phase space directions up to a small distance. Furthermore, we clarify that the stiffness of the cells tunes this anomalous behavior, tissues of stiff cells having conventional glassy relaxation dynamics.

Hyperuniformity and density fluctuations at a rigidity transition in a model of biological tissues

The suppression of density fluctuations at different length scales is the hallmark of hyperuniformity. Here, we explore the presence of this hidden order in a manybody interacting model of biological tissue, known to exhibit a transition, or sharp crossover, from a solid to a fluid like phase. We show that the density fluctuations in the rigid phase are only suppressed up to a finite lengthscale. This length scale monotonically increases and grows rapidly as we approach the fluid phase reminiscent to divergent behavior at a critical point, such that the system is effectively hyperuniform in the fluid phase. Furthermore, complementary behavior of the structure factor across the critical point also indicates that hyperuniformity found in the fluid phase is stealthy.

Small-scale demixing in confluent biological tissues

Surface tension governed by differential adhesion can drive fluid particle mixtures to sort into separate regions, i.e., demix. Does the same phenomenon occur in confluent biological tissues? We begin to answer this question for epithelial monolayers with a combination of theory via a vertex model and experiments on keratinocyte monolayers. Vertex models are distinct from particle models in that the interactions between the cells are shape-based, as opposed to distance-dependent. We investigate whether a disparity in cell shape or size alone is sufficient to drive demixing in bidisperse vertex model fluid mixtures. Surprisingly, we observe that both types of bidisperse systems robustly mix on large lengthscales. On the other hand, shape disparity generates slight demixing over a few cell diameters, a phenomenon we term micro-demixing. This result can be understood by examining the differential energy barriers for neighbor exchanges (T1 transitions). Experiments with mixtures of wild-type and E-cadherin-deficient keratinocytes on a substrate are consistent with the predicted phenomenon of micro-demixing, which biology may exploit to create subtle patterning. The robustness of mixing at large scales, however, suggests that despite some differences in cell shape and size, progenitor cells can readily mix throughout a developing tissue until acquiring means of recognizing cells of different types.

Linear and nonlinear mechanical responses can be quite different in models for biological tissues

The fluidity of biological tissues - whether cells can change neighbors and rearrange - is important for their function. In traditional materials, researchers have used linear response functions, such as the shear modulus, to accurately predict whether a material will behave as a fluid. Similarly, in disordered 2D vertex models for confluent biological tissues, the shear modulus becomes zero precisely when the cells can change neighbors and the tissue fluidizes, at a critical value of control parameter s0* = 3.81. However, the ordered ground states of 2D vertex models become linearly unstable at a lower value of control parameter (3.72), suggesting that there may be a decoupling between linear and nonlinear response. We demonstrate that the linear response does not correctly predict the nonlinear behavior in these systems: when the control parameter is between 3.72 and 3.81, cells cannot freely change neighbors even though the shear modulus is zero. These results highlight that the linear response of vertex models should not be expected to generically predict their rheology. We develop a simple geometric ansatz that correctly predicts the nonlinear response, which may serve as a framework for making nonlinear predictions in other vertex-like models.

Supersonic kinks and solitons in active solids

To show that steadily propagating nonlinear waves in active matter can be driven internally, we develop a prototypical model of a topological kink moving with a constant supersonic speed. We use a model of a bi-stable mass-spring (Fermi-Pasta-Ulam) chain capable of generating active stress. In contrast to subsonic kinks in passive bi-stable chains that are necessarily dissipative, the obtained supersonic solutions are purely anti-dissipative. Our numerical experiments point towards the stability of the obtained kink-type solutions and the possibility of propagating kink-anti-kink bundles reminiscent of solitons. We show that even the simplest quasi-continuum approximation of the discrete model captures the most important features of the predicted active phenomena. This article is part of the theme issue 'Modelling of dynamic phenomena and localization in structured media (part 2)'.

Mechanical Heterogeneity in Tissues Promotes Rigidity and Controls Cellular Invasion

We study the influence of cell-level mechanical heterogeneity in epithelial tissues using a vertex-based model. Heterogeneity is introduced into the cell shape index (p_{0}) that tunes the stiffness at a single-cell level. The addition of heterogeneity can always enhance the mechanical rigidity of the epithelial layer by increasing its shear modulus, hence making it more rigid. There is an excellent scaling collapse of our data as a function of a single scaling variable f_{r}, which accounts for the overall fraction of rigid cells. We identify a universal threshold f_{r}^{*} that demarcates fluid versus solid tissues. Furthermore, this rigidity onset is far below the contact percolation threshold of rigid cells. These results give rise to a separation of rigidity and contact percolation processes that leads to distinct types of solid states. We also investigate the influence of heterogeneity on tumor invasion dynamics. There is an overall impedance of invasion as the tissue becomes more rigid. Invasion can also occur in an intermediate heterogeneous solid state that is characterized by significant spatial-temporal intermittency.

A minimal-length approach unifies rigidity in underconstrained materials

We present an approach to understand geometric-incompatibility-induced rigidity in underconstrained materials, including subisostatic 2D spring networks and 2D and 3D vertex models for dense biological tissues. We show that in all these models a geometric criterion, represented by a minimal length [Formula: see text], determines the onset of prestresses and rigidity. This allows us to predict not only the correct scalings for the elastic material properties, but also the precise magnitudes for bulk modulus and shear modulus discontinuities at the rigidity transition as well as the magnitude of the Poynting effect. We also predict from first principles that the ratio of the excess shear modulus to the shear stress should be inversely proportional to the critical strain with a prefactor of 3. We propose that this factor of 3 is a general hallmark of geometrically induced rigidity in underconstrained materials and could be used to distinguish this effect from nonlinear mechanics of single components in experiments. Finally, our results may lay important foundations for ways to estimate [Formula: see text] from measurements of local geometric structure and thus help develop methods to characterize large-scale mechanical properties from imaging data.

A simple landscape of metastable state energies for two-dimensional cellular matter

The mechanical behavior of cellular matter in two dimensions can be inferred from geometric information near its energetic ground state. Here it is shown that the much larger set of all metastable state energies is universally described by a systematic expansion in moments of the joint probability distribution of size (area) and topology (number of neighbors). The approach captures bounds to the entire range of metastable state energies and quantitatively identifies any such state. The resulting energy landscape is invariant across different classes of energy functionals, across simulation techniques, and across system polydispersities. The theory also finds a threshold in tissue adhesion beyond which no metastable states are possible. Mechanical properties of cellular matter in biological and technological applications can thus be identified by visual information only.