Linear viscoelastic properties of the vertex model for epithelial tissues

i-Rheo-optical assay: Measuring the viscoelastic properties of multicellular spheroids

This study introduces a novel mechanobiology assay, named "i-Rheo-optical assay", that integrates rheology with optical microscopy for analysing the viscoelastic properties of multicellular spheroids. These spheroids serve as three-dimensional models resembling tissue structures. The innovative technique enables real-time observation and quantification of morphological responses to applied stress using a cost-effective microscope coverslip for constant compression force application. By bridging a knowledge gap in biophysical research, which has predominantly focused on the elastic properties while only minimally exploring the viscoelastic nature in multicellular systems, the i-Rheo-optical assay emerges as an effective tool. It facilitates the measurement of broadband viscoelastic compressional moduli in spheroids, here derived from cancer (PANC-1) and non-tumoral (NIH/3T3) cell lines during compression tests. This approach plays a crucial role in elucidating the mechanical properties of spheroids and holds potential for identifying biomarkers to discriminate between healthy tissues and their pathological counterparts. Offering comprehensive insights into the biomechanical behaviour of biological systems, i-Rheo-optical assay marks a significant advancement in tissue engineering, cancer research, and therapeutic development.

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.

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.

Couple stresses and discrete potentials in the vertex model of cellular monolayers

The vertex model is widely used to simulate the mechanical properties of confluent epithelia and other multicellular tissues. This inherently discrete framework allows a Cauchy stress to be attributed to each cell, and its symmetric component has been widely reported, at least for planar monolayers. Here, we consider the stress attributed to the neighbourhood of each tricellular junction, evaluating in particular its leading-order antisymmetric component and the associated couple stresses, which characterise the degree to which individual cells experience (and resist) in-plane bending deformations. We develop discrete potential theory for localised monolayers having disordered internal structure and use this to derive the analogues of Airy and Mindlin stress functions. These scalar potentials typically have broad-banded spectra, highlighting the contributions of small-scale defects and boundary layers to global stress patterns. An affine approximation attributes couple stresses to pressure differences between cells sharing a trijunction, but simulations indicate an additional role for non-affine deformations.

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.

Interplay between substrate rigidity and tissue fluidity regulates cell monolayer spreading

Coordinated and cooperative motion of cells is essential for embryonic development, tissue morphogenesis, wound healing and cancer invasion. A predictive understanding of the emergent mechanical behaviors in collective cell motion is challenging due to the complex interplay between cell-cell interactions, cell-matrix adhesions and active cell behaviors. To overcome this challenge, we develop a predictive cellular vertex model that can delineate the relative roles of substrate rigidity, tissue mechanics and active cell properties on the movement of cell collectives. We apply the model to the specific case of collective motion in cell aggregates as they spread into a two-dimensional cell monolayer adherent to a soft elastic matrix. Consistent with recent experiments, we find that substrate stiffness regulates the driving forces for the spreading of cellular monolayer, which can be pressure-driven or crawling-based depending on substrate rigidity. On soft substrates, cell monolayer spreading is driven by an active pressure due to the influx of cells coming from the aggregate, whereas on stiff substrates, cell spreading is driven primarily by active crawling forces. Our model predicts that cooperation of cell crawling and tissue pressure drives faster spreading, while the spreading rate is sensitive to the mechanical properties of the tissue. We find that solid tissues spread faster on stiff substrates, with spreading rate increasing with tissue tension. By contrast, the spreading of fluid tissues is independent of substrate stiffness and is slower than solid tissues. We compare our theoretical results with experimental results on traction force generation and spreading kinetics of cell monolayers, and provide new predictions on the role of tissue fluidity and substrate rigidity on collective cell motion.

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.

Mapping cell cortex rheology to tissue rheology and vice versa

The mechanics of biological tissues mainly proceeds from the cell cortex rheology. A direct, explicit link between cortex rheology and tissue rheology remains lacking, yet would be instrumental in understanding how modulations of cortical mechanics may impact tissue mechanical behavior. Using an ordered geometry built on 3D hexagonal, incompressible cells, we build a mapping relating the cortical rheology to the monolayer tissue rheology. Our approach shows that the tissue low-frequency elastic modulus is proportional to the rest tension of the cortex, as expected from the physics of liquid foams as well as of tensegrity structures. A fractional visco-contractile cortex rheology is predicted to yield a high-frequency fractional visco-elastic monolayer rheology, where such a fractional behavior has been recently observed experimentally at each scale separately. In particular cases, the mapping may be inverted, allowing to derive from a given tissue rheology the underlying cortex rheology. Interestingly, applying the same approach to a 2D hexagonal tiling fails, which suggests that the 2D character of planar cell cortex-based models may be unsuitable to account for realistic monolayer rheologies. We provide quantitative predictions, amenable to experimental tests through standard perturbation assays of cortex constituents, and hope to foster new, challenging mechanical experiments on cell monolayers.

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.