Correlating cell shape and cellular stress in motile confluent tissues

CP-Net: Instance-aware part segmentation network for biological cell parsing

Instance segmentation of biological cells is important in medical image analysis for identifying and segmenting individual cells, and quantitative measurement of subcellular structures requires further cell-level subcellular part segmentation. Subcellular structure measurements are critical for cell phenotyping and quality analysis. For these purposes, instance-aware part segmentation network is first introduced to distinguish individual cells and segment subcellular structures for each detected cell. This approach is demonstrated on human sperm cells since the World Health Organization has established quantitative standards for sperm quality assessment. Specifically, a novel Cell Parsing Net (CP-Net) is proposed for accurate instance-level cell parsing. An attention-based feature fusion module is designed to alleviate contour misalignments for cells with an irregular shape by using instance masks as spatial cues instead of as strict constraints to differentiate various instances. A coarse-to-fine segmentation module is developed to effectively segment tiny subcellular structures within a cell through hierarchical segmentation from whole to part instead of directly segmenting each cell part. Moreover, a sperm parsing dataset is built including 320 annotated sperm images with five semantic subcellular part labels. Extensive experiments on the collected dataset demonstrate that the proposed CP-Net outperforms state-of-the-art instance-aware part segmentation networks.

Minimal cellular automaton model with heterogeneous cell sizes predicts epithelial colony growth

Regulation of cell proliferation is a crucial aspect of tissue development and homeostasis and plays a major role in morphogenesis, wound healing, and tumor invasion. A phenomenon of such regulation is contact inhibition, which describes the dramatic slowing of proliferation, cell migration and individual cell growth when multiple cells are in contact with each other. While many physiological, molecular and genetic factors are known, the mechanism of contact inhibition is still not fully understood. In particular, the relevance of cellular signaling due to interfacial contact for contact inhibition is still debated. Cellular automata (CA) have been employed in the past as numerically efficient mathematical models to study the dynamics of cell ensembles, but they are not suitable to explore the origins of contact inhibition as such agent-based models assume fixed cell sizes. We develop a minimal, data-driven model to simulate the dynamics of planar cell cultures by extending a probabilistic CA to incorporate size changes of individual cells during growth and cell division. We successfully apply this model to previous in-vitro experiments on contact inhibition in epithelial tissue: After a systematic calibration of the model parameters to measurements of single-cell dynamics, our CA model quantitatively reproduces independent measurements of emergent, culture-wide features, like colony size, cell density and collective cell migration. In particular, the dynamics of the CA model also exhibit the transition from a low-density confluent regime to a stationary postconfluent regime with a rapid decrease in cell size and motion. This implies that the volume exclusion principle, a mechanical constraint which is the only inter-cellular interaction incorporated in the model, paired with a size-dependent proliferation rate is sufficient to generate the observed contact inhibition. We discuss how our approach enables the introduction of effective bio-mechanical interactions in a CA framework for future studies.

Motility driven glassy dynamics in confluent epithelial monolayers

As wounds heal, embryos develop, cancer spreads, or asthma progresses, the cellular monolayer undergoes a glass transition between solid-like jammed and fluid-like flowing states. During some of these processes, the cells undergo an epithelial-to-mesenchymal transition (EMT): they acquire in-plane polarity and become motile. Thus, how motility drives the glassy dynamics in epithelial systems is critical for the EMT process. However, no analytical framework that is indispensable for deeper insights exists. Here, we develop such a theory inspired by a well-known glass theory. One crucial result of this work is that the confluency affects the effective persistence time-scale of active force, described by its rotational diffusivity, Deffr. Deffr differs from the bare rotational diffusivity, Dr, of the motile force due to cell shape dynamics, which acts to rectify the force dynamics: Deffr is equal to Dr when Dr is small and saturates when Dr is large. We test the theoretical prediction of Deffr and how it affects the relaxation dynamics in our simulations of the active Vertex model. This novel effect of Deffr is crucial to understanding the new and previously published simulation data of active glassy dynamics in epithelial monolayers.

Micro/Nanomotor-Driven Intelligent Targeted Delivery Systems: Dynamics Sources and Frontier Applications

Micro/nanomotors represent a promising class of drug delivery carriers capable of converting surrounding chemical or external energy into mechanical power, enabling autonomous movement. Their distinct autonomous propulsive force distinguishes them from other carriers, offering significant potential for enhancing drug penetration across cellular and tissue barriers. A comprehensive understanding of micro/nanomotor dynamics with various power sources is crucial to facilitate their transition from proof-of-concept to clinical application. In this review, micro/nanomotors are categorized into three classes based on their energy sources: endogenously stimulated, exogenously stimulated, and live cell-driven. The review summarizes the mechanisms governing micro/nanomotor movements under these energy sources and explores factors influencing autonomous motion. Furthermore, it discusses methods for controlling micro/nanomotor movement, encompassing aspects related to their structure, composition, and environmental factors. The remarkable propulsive force exhibited by micro/nanomotors makes them valuable for significant biomedical applications, including tumor therapy, bio-detection, bacterial infection therapy, inflammation therapy, gastrointestinal disease therapy, and environmental remediation. Finally, the review addresses the challenges and prospects for the application of micro/nanomotors. Overall, this review emphasizes the transformative potential of micro/nanomotors in overcoming biological barriers and enhancing therapeutic efficacy, highlighting their promising clinical applications across various biomedical fields.

Mechanical imbalance between normal and transformed cells drives epithelial homeostasis through cell competition

Cell competition in epithelial tissue eliminates transformed cells expressing activated oncoproteins to maintain epithelial homeostasis. Although the process is now understood to be of mechanochemical origin, direct mechanical characterization and associated biochemical underpinnings are lacking. Here, we employ tissue-scale stress and compressibility measurements and theoretical modeling to unveil a mechanical imbalance between normal and transformed cells, which drives cell competition. In the mouse intestinal epithelium and epithelial monolayer, transformed cells get compacted during competition. Stress microscopy reveals an emergent compressive stress at the transformed loci leading to this compaction. A cell-based self-propelled Voronoi model predicts that this compressive stress originates from a difference in the collective compressibility of the competing populations. A new collective compressibility measurement technique named gel compression microscopy then elucidates a two-fold higher compressibility of the transformed population than the normal population. Mechanistically, weakened cell-cell adhesions due to reduced junctional abundance of E-cadherin in the transformed cells render them collectively more compressible than normal cells. Taken together, our findings unveil a mechanical basis for epithelial homeostasis against oncogenic transformations with implications in epithelial defense against cancer.

Spatial confinement affects the heterogeneity and interactions between shoaling fish

Living objects are able to consume chemical energy and process information independently from others. However, living objects can coordinate to form ordered groups such as schools of fish. This work considers these complex groups as living materials and presents imaging-based experiments of laboratory schools of fish to understand how activity, which is a non-equilibrium feature, affects the structure and dynamics of a group. We use spatial confinement to control the motion and structure of fish within quasi-2D shoals of fish and use image analysis techniques to make quantitative observations of the structures, their spatial heterogeneity, and their temporal fluctuations. Furthermore, we utilize Monte Carlo simulations to replicate the experimentally observed data which provides insight into the effective interactions between fish and confirms the presence of a confinement-based behavioral preference transition. In addition, unlike in short-range interacting systems, here structural heterogeneity and dynamic activities are positively correlated as a result of complex interplay between spatial arrangement and behavioral dynamics in fish collectives.

Environmentally dependent and independent control of 3D cell shape

How cancer cells determine their shape in response to three-dimensional (3D) geometric and mechanical cues is unclear. We develop an approach to quantify the 3D cell shape of over 60,000 melanoma cells in collagen hydrogels using high-throughput stage-scanning oblique plane microscopy (ssOPM). We identify stereotypic and environmentally dependent changes in shape and protrusivity depending on whether a cell is proximal to a flat and rigid surface or is embedded in a soft environment. Environmental sensitivity metrics calculated for small molecules and gene knockdowns identify interactions between the environment and cellular factors that are important for morphogenesis. We show that the Rho guanine nucleotide exchange factor (RhoGEF) TIAM2 contributes to shape determination in environmentally independent ways but that non-muscle myosin II, microtubules, and the RhoGEF FARP1 regulate shape in ways dependent on the microenvironment. Thus, changes in cancer cell shape in response to 3D geometric and mechanical cues are modulated in both an environmentally dependent and independent fashion.

Mechanically-driven growth and competition in a Voronoi model of tissues

The mechanisms leading cells to acquire a fitness advantage and establish themselves in a population are paramount to understanding the development and growth of cancer. Although there are many works that study separately either the evolutionary dynamics or the mechanics of cancer, little has been done to couple evolutionary dynamics to mechanics. To address this question, we study a confluent model of tissue using a Self-Propelled Voronoi (SPV) model with stochastic growth rates that depend on the mechanical variables of the system. The SPV model is an out-of-equilibrium model of tissue derived from an energy functional that has a jamming/unjamming transition between solid-like and liquid-like states. By considering several scenarios of mutants invading a resident population in both phases, we determine the range of parameters that confer a fitness advantage and show that the preferred area and perimeter are the most relevant ones. We find that the liquid-like state is more resistant to invasion and show that the outcome of the competition can be determined from the simulation of a non-growing mixture. Moreover, a mean-field approximation can accurately predict the fate of a mutation affecting mechanical properties of a cell. Our results can be used to infer evolutionary dynamics from tissue images, understand cancer-suppressing effects of tissue mechanics, and even search for mechanics-based therapies.

Matrix confinement modulates 3D spheroid sorting and burst-like collective migration

While it is known that cells with differential adhesion tend to segregate and preferentially sort, the physical forces governing sorting and invasion in heterogeneous tumors remain poorly understood. To investigate this, we tune matrix confinement, mimicking changes in the stiffness and confinement of the tumor microenvironment, to explore how physical confinement influences individual and collective cell migration in 3D spheroids. High levels of confinement lead to cell sorting while reducing matrix confinement triggers the collective fluidization of cell motion. Cell sorting, which depends on cell-cell adhesion, is crucial to this phenomenon. Burst-like migration does not occur for spheroids that have not undergone sorting, regardless of the degree of matrix confinement. Using computational Self-Propelled Voronoi modeling, we show that spheroid sorting and invasion into the matrix depend on the balance between cell-generated forces and matrix resistance. The findings support a model where matrix confinement modulates 3D spheroid sorting and unjamming in an adhesion-dependent manner, providing insights into the mechanisms of cell sorting and migration in the primary tumor and toward distant metastatic sites. STATEMENT OF SIGNIFICANCE: The mechanical properties of the tumor microenvironment significantly influence cancer cell migration within the primary tumor, yet how these properties affect intercellular interactions in heterogeneous tumors is not well understood. By utilizing calcium and calcium chelators, we dynamically alter collagen-alginate hydrogel stiffness and investigate tumor cell behavior within co-culture spheroids in response to varying degrees of matrix confinement. High confinement is found to trigger cell sorting while reducing confinement for sorted spheroids facilitates collective cell invasion. Notably, without prior sorting, spheroids do not exhibit burst-like migration, regardless of confinement levels. This work establishes that matrix confinement and intercellular adhesion regulate 3D spheroid dynamics, offering insights into cellular organization and migration within the primary tumor.

Matrix confinement modulates 3D spheroid sorting and burst-like collective migration

While it is known that cells with differential adhesion tend to segregate and preferentially sort, the physical forces governing sorting and invasion in heterogeneous tumors remain poorly understood. To investigate this, we tune matrix confinement, mimicking changes in the stiffness and confinement of the tumor microenvironment, to explore how physical confinement influences individual and collective cell migration in 3D spheroids. High levels of confinement lead to cell sorting while reducing matrix confinement triggers the collective fluidization of cell motion. Cell sorting, which depends on cell-cell adhesion, is crucial to this phenomenon. Burst-like migration does not occur for spheroids that have not undergone sorting, regardless of the degree of matrix confinement. Using computational Self-Propelled Voronoi modeling, we show that spheroid sorting and invasion into the matrix depend on the balance between cell-generated forces and matrix resistance. The findings support a model where matrix confinement modulates 3D spheroid sorting and unjamming in an adhesion-dependent manner, providing insights into the mechanisms of cell sorting and migration in the primary tumor and toward distant metastatic sites.

Discontinuous Shear Thickening in Biological Tissue Rheology

During embryonic morphogenesis, tissues undergo dramatic deformations in order to form functional organs. Similarly, in adult animals, living cells and tissues are continually subjected to forces and deformations. Therefore, the success of embryonic development and the proper maintenance of physiological functions rely on the ability of cells to withstand mechanical stresses as well as their ability to flow in a collective manner. During these events, mechanical perturbations can originate from active processes at the single-cell level, competing with external stresses exerted by surrounding tissues and organs. However, the study of tissue mechanics has been somewhat limited to either the response to external forces or to intrinsic ones. In this work, we use an active vertex model of a 2D confluent tissue to study the interplay of external deformations that are applied globally to a tissue with internal active stresses that arise locally at the cellular level due to cell motility. We elucidate, in particular, the way in which this interplay between globally external and locally internal active driving determines the emergent mechanical properties of the tissue as a whole. For a tissue in the vicinity of a solid-fluid jamming or unjamming transition, we uncover a host of fascinating rheological phenomena, including yielding, shear thinning, continuous shear thickening, and discontinuous shear thickening. These model predictions provide a framework for understanding the recently observed nonlinear rheological behaviors in vivo.

Minimal vertex model explains how the amnioserosa avoids fluidization during Drosophila dorsal closure

Dorsal closure is a process that occurs during embryogenesis of Drosophila melanogaster . During dorsal closure, the amnioserosa (AS), a one-cell thick epithelial tissue that fills the dorsal opening, shrinks as the lateral epidermis sheets converge and eventually merge. During this process, the aspect ratio of amnioserosa cells increases markedly. The standard 2-dimensional vertex model, which successfully describes tissue sheet mechanics in multiple contexts, would in this case predict that the tissue should fluidize via cell neighbor changes. Surprisingly, however, the amnioserosa remains an elastic solid with no such events. We here present a minimal extension to the vertex model that explains how the amnioserosa can achieve this unexpected behavior. We show that continuous shrinkage of the preferred cell perimeter and cell perimeter polydispersity lead to the retention of the solid state of the amnioserosa. Our model accurately captures measured cell shape and orientation changes and predicts non-monotonic junction tension that we confirm with laser ablation experiments.

Significance Statement

During embryogenesis, cells in tissues can undergo significant shape changes. Many epithelial tissues fluidize, i.e. cells exchange neighbors, when the average cell aspect ratio increases above a threshold value, consistent with the standard vertex model. During dorsal closure in Drosophila melanogaster , however, the amnioserosa tissue remains solid even as the average cell aspect ratio increases well above threshold. We introduce perimeter polydispersity and allow the preferred cell perimeters, usually held fixed in vertex models, to decrease linearly with time as seen experimentally. With these extensions to the standard vertex model, we capture experimental observations quantitatively. Our results demonstrate that vertex models can describe the behavior of the amnioserosa in dorsal closure by allowing normally fixed parameters to vary with time.

Minimal vertex model explains how the amnioserosa avoids fluidization during Drosophila dorsal closure

Dorsal closure is a process that occurs during embryogenesis of Drosophila melanogaster. During dorsal closure, the amnioserosa (AS), a one-cell thick epithelial tissue that fills the dorsal opening, shrinks as the lateral epidermis sheets converge and eventually merge. During this process, the aspect ratio of amnioserosa cells increases markedly. The standard 2-dimensional vertex model, which successfully describes tissue sheet mechanics in multiple contexts, would in this case predict that the tissue should fluidize via cell neighbor changes. Surprisingly, however, the amnioserosa remains an elastic solid with no such events. We here present a minimal extension to the vertex model that explains how the amnioserosa can achieve this unexpected behavior. We show that continuous shrink-age of the preferred cell perimeter and cell perimeter polydispersity lead to the retention of the solid state of the amnioserosa. Our model accurately captures measured cell shape and orientation changes and predicts non-monotonic junction tension that we confirm with laser ablation experiments.

Bridging the gap between collective motility and epithelial-mesenchymal transitions through the active finite voronoi model

We introduce an active version of the recently proposed finite Voronoi model of epithelial tissue. The resultant Active Finite Voronoi (AFV) model enables the study of both confluent and non-confluent geometries and transitions between them, in the presence of active cells. Our study identifies six distinct phases, characterized by aggregation-segregation, dynamical jamming-unjamming, and epithelial-mesenchymal transitions (EMT), thereby extending the behavior beyond that observed in previously studied vertex-based models. The AFV model with rich phase diagram provides a cohesive framework that unifies the well-observed progression to collective motility via unjamming with the intricate dynamics enabled by EMT. This approach should prove useful for challenges in developmental biology systems as well as the complex context of cancer metastasis. The simulation code is also provided.

Nature-inspired designs for disordered acoustic bandgap materials

We introduce an amorphous mechanical metamaterial inspired by how cells pack in biological tissues. The spatial heterogeneity in the local stiffness of these materials has been recently shown to impact the mechanics of confluent biological tissues and cancer tumor invasion. Here we use this bio-inspired structure as a design template to construct mechanical metamaterials and show that this heterogeneity can give rise to amorphous cellular solids with large, tunable acoustic bandgaps. Unlike acoustic crystals with periodic structures, the bandgaps here are directionally isotropic and robust to defects due to their complete lack of positional order. Possible ways to manipulate bandgaps are explored with a combination of the tissue-level elastic modulus and local stiffness heterogeneity of cells. To further demonstrate the existence of bandgaps, we dynamically perturb the system with an external sinusoidal wave in the perpendicular and horizontal directions. The transmission coefficients are calculated and show valleys that coincide with the location of bandgaps. Experimentally this design should lead to the engineering of self-assembled rigid acoustic structures with full bandgaps that can be controlled via mechanical tuning and promote applications in a broad area from vibration isolations to mechanical waveguides.

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.

Long-term, real-time and label-free live cell image processing and analysis based on a combined algorithm of CellPose and watershed segmentation

Developing a rapid and quantitative method to accurately evaluate the physiological abilities of living cells is critical for tumor control. Many experiments have been conducted in the field of biology in an attempt to measure the proliferation and movement abilities of cells, but existing methods cannot provide real-time and objective data for label-free cells. The quantitative imaging technique, including an automatic segmentation algorithm for individual label-free cells, has been a breakthrough in this regard. In this study, we develop a combined automatic image processing algorithm of CellPose and watershed segmentation for the long-term and real-time imaging of label-free cells. This method shows strong reliability in cell identification regardless of cell densities, allowing us to obtain accurate information about the number and proliferation ability of the target cells. Additionally, our results also suggest that this method is a reliable way to assess real-time data on drug cytotoxicity, cell morphology, and cell movement ability.

Force propagation between epithelial cells depends on active coupling and mechano-structural polarization

Cell-generated forces play a major role in coordinating the large-scale behavior of cell assemblies, in particular during development, wound healing, and cancer. Mechanical signals propagate faster than biochemical signals, but can have similar effects, especially in epithelial tissues with strong cell-cell adhesion. However, a quantitative description of the transmission chain from force generation in a sender cell, force propagation across cell-cell boundaries, and the concomitant response of receiver cells is missing. For a quantitative analysis of this important situation, here we propose a minimal model system of two epithelial cells on an H-pattern ('cell doublet'). After optogenetically activating RhoA, a major regulator of cell contractility, in the sender cell, we measure the mechanical response of the receiver cell by traction force and monolayer stress microscopies. In general, we find that the receiver cells show an active response so that the cell doublet forms a coherent unit. However, force propagation and response of the receiver cell also strongly depend on the mechano-structural polarization in the cell assembly, which is controlled by cell-matrix adhesion to the adhesive micropattern. We find that the response of the receiver cell is stronger when the mechano-structural polarization axis is oriented perpendicular to the direction of force propagation, reminiscent of the Poisson effect in passive materials. We finally show that the same effects are at work in small tissues. Our work demonstrates that cellular organization and active mechanical response of a tissue are key to maintain signal strength and lead to the emergence of elasticity, which means that signals are not dissipated like in a viscous system, but can propagate over large distances.

Entropons as collective excitations in active solids

The vibrational dynamics of solids is described by phonons constituting basic collective excitations in equilibrium crystals. Here, we consider a non-equilibrium active solid, formed by self-propelled particles, which bring the system into a non-equilibrium steady-state. We identify novel vibrational collective excitations of non-equilibrium (active) origin, which coexist with phonons and dominate over them when the system is far from equilibrium. These vibrational excitations are interpreted in the framework of non-equilibrium physics, in particular, stochastic thermodynamics. We call them "entropons" because they are the modes of spectral entropy production (at a given frequency and wave vector). The existence of entropons could be verified in future experiments on dense self-propelled colloidal Janus particles and granular active matter, as well as in living systems, such as dense cell monolayers.

Migration and division in cell monolayers on substrates with topological defects

Collective movement and organization of cell monolayers are important for wound healing and tissue development. Recent experiments highlighted the importance of liquid crystal order within these layers, suggesting that +1 topological defects have a role in organizing tissue morphogenesis. We study fibroblast organization, motion, and proliferation on a substrate with micron-sized ridges that induce +1 and -1 topological defects using simulation and experiment. We model cells as self-propelled deformable ellipses that interact via a Gay-Berne potential. Unlike earlier work on other cell types, we see that density variation near defects is not explained by collective migration. We propose instead that fibroblasts have different division rates depending on their area and aspect ratio. This model captures key features of our previous experiments: the alignment quality worsens at high cell density and, at the center of the +1 defects, cells can adopt either highly anisotropic or primarily isotropic morphologies. Experiments performed with different ridge heights confirm a prediction of this model: Suppressing migration across ridges promotes higher cell density at the +1 defect. Our work enables a mechanism for tissue patterning using topological defects without relying on cell migration.

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.

Coordination of non-professional efferocytosis and actomyosin contractility during epithelial tissue morphogenesis

In developing embryos, specific cell populations are often removed to remodel tissue architecture for organogenesis. During urinary tract development, an epithelial duct called the common nephric duct (CND) gets shortened and eventually eliminated to remodel the entry point of the ureter into the bladder. Here we show that non-professional efferocytosis (the process in which epithelial cells engulf apoptotic bodies) is the main mechanism that contributes to CND shortening. Combining biological metrics and computational modeling, we show that efferocytosis with actomyosin contractility are essential factors that drive the CND shortening without compromising the ureter-bladder structural connection. The disruption of either apoptosis, non-professional efferocytosis, or actomyosin results in contractile tension reduction and deficient CND shortening. Actomyosin activity helps to maintain tissue architecture while non-professional efferocytosis removes cellular volume. Together our results demonstrate that non-professional efferocytosis with actomyosin contractility are important morphogenetic factors controlling CND morphogenesis.

Collective Behaviors of Active Matter Learning from Natural Taxes Across Scales

Taxis orientation is common in microorganisms, and it provides feasible strategies to operate active colloids as small-scale robots. Collective taxes involve numerous units that collectively perform taxis motion, whereby the collective cooperation between individuals enables the group to perform efficiently, adaptively, and robustly. Hence, analyzing and designing collectives is crucial for developing and advancing microswarm toward practical or clinical applications. In this review, natural taxis behaviors are categorized and synthetic microrobotic collectives are discussed as bio-inspired realizations, aiming at closing the gap between taxis strategies of living creatures and those of functional active microswarms. As collective behaviors emerge within a group, the global taxis to external stimuli guides the group to conduct overall tasks, whereas the local taxis between individuals induces synchronization and global patterns. By encoding the local orientations and programming the global stimuli, various paradigms can be introduced for coordinating and controlling such collective microrobots, from the viewpoints of fundamental science and practical applications. Therefore, by discussing the key points and difficulties associated with collective taxes of different paradigms, this review potentially offers insights into mimicking natural collective behaviors and constructing intelligent microrobotic systems for on-demand control and preassigned tasks.

A quantitative biophysical principle to explain the 3D cellular connectivity in curved epithelia

Epithelial cell organization and the mechanical stability of tissues are closely related. In this context, it has been recently shown that packing optimization in bended or folded epithelia is achieved by an energy minimization mechanism that leads to a complex cellular shape: the "scutoid". Here, we focus on the relationship between this shape and the connectivity between cells. We use a combination of computational, experimental, and biophysical approaches to examine how energy drivers affect the three-dimensional (3D) packing of tubular epithelia. We propose an energy-based stochastic model that explains the 3D cellular connectivity. Then, we challenge it by experimentally reducing the cell adhesion. As a result, we observed an increment in the appearance of scutoids that correlated with a decrease in the energy barrier necessary to connect with new cells. We conclude that tubular epithelia satisfy a quantitative biophysical principle that links tissue geometry and energetics with the average cellular connectivity.

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.

Polar Fluctuations Lead to Extensile Nematic Behavior in Confluent Tissues

How can a collection of motile cells, each generating contractile nematic stresses in isolation, become an extensile nematic at the tissue level? Understanding this seemingly contradictory experimental observation, which occurs irrespective of whether the tissue is in the liquid or solid states, is not only crucial to our understanding of diverse biological processes, but is also of fundamental interest to soft matter and many-body physics. Here, we resolve this cellular to tissue level disconnect in the small fluctuation regime by using analytical theories based on hydrodynamic descriptions of confluent tissues, in both liquid and solid states. Specifically, we show that a collection of microscopic constituents with no inherently nematic extensile forces can exhibit active extensile nematic behavior when subject to polar fluctuating forces. We further support our findings by performing cell level simulations of minimal models of confluent tissues.

Cell monolayer deformation microscopy reveals mechanical fragility of cell monolayers following EMT

Tissue and cell mechanics are crucial factors in maintaining homeostasis and in development, with aberrant mechanics contributing to many diseases. During the epithelial-to-mesenchymal transition (EMT), a highly conserved cellular program in organismal development and cancer metastasis, cells gain the ability to detach from their original location and autonomously migrate. While a great deal of biochemical and biophysical changes at the single-cell level have been revealed, how the physical properties of multicellular assemblies change during EMT, and how this may affect disease progression, is unknown. Here we introduce cell monolayer deformation microscopy (CMDM), a new methodology to measure the planar mechanical properties of cell monolayers by locally applying strain and measuring their resistance to deformation. We employ this new method to characterize epithelial multicellular mechanics at early and late stages of EMT, finding the epithelial monolayers to be relatively compliant, ductile, and mechanically homogeneous. By comparison, the transformed mesenchymal monolayers, while much stiffer, were also more brittle, mechanically heterogeneous, displayed more viscoelastic creep, and showed sharp yield points at significantly lower strains. Here, CMDM measurements identify specific biophysical functional states of EMT and offer insight into how cell aggregates fragment under mechanical stress. This mechanical fingerprinting of multicellular assemblies using new quantitative metrics may also offer new diagnostic applications in healthcare to characterize multicellular mechanical changes in disease.

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.

Activation of Topological Defects Induces a Brittle-to-Ductile Transition in Epithelial Monolayers

Epithelial monolayers are subjected to various mechanical forces, such as stretching, shearing, and compression. Thus, its mechanical response to external loadings is essential for its biological functions. However, the mechanism of the fracture failure of the epithelial monolayer remains poorly understood. Here, by introducing a new type of topological transition, i.e., detach transition or T4 transition, we develop a modified cellular vertex model to investigate the rupture of the cell monolayer. Interestingly, we find a brittle-to-ductile transition in epithelial monolayers, which is controlled by the mechanical properties of single cells and cell-cell contacts. We reveal that the external loadings can activate cell rearrangement in ductile cell monolayers. The plastic deformation results from the nucleation and propagation of "pentagon-heptagon defects" in analogy with the topological defects commonly seen in 2D materials. By using a simplified four-cell model, we further demonstrate that the brittle-to-ductile transition is induced by the competition between cell rearrangement and cell detachment. Our work provides a new theoretical framework to study the rupture of living tissues and may have important implications for many other biological processes, such as wound healing and tissue morphogenesis.

Large-scale demixing in a binary mixture of cells with rigidity disparity in biological tissues

Physical demixing on large scales of embryonic cell populations is fundamental to metazoan development, but whether a rigidity disparity alone is sufficient to driving large-scale demixing in a binary mixture of cell tissues is still an open question. To answer this question, we study mixing and demixing in a binary mixture of rigidity disparity cell tissues without heterotypic interactions using the Voronoi-based cellular model. Under suitable system parameters, the solid-like cells in the mixture can aggregate into a large cluster and the large-scale demixing occurs, which addresses that a rigidity disparity alone is sufficient to drive large-scale demixing. Remarkably, there exists an optimal temperature or rigidity disparity at which the binary mixture can be separated to the maximum extent. The necessary condition for the separation of mixtures is that the two types of cells are solid-like and liquid-like, respectively. The observation of robust demixing on large scales suggests that the sorting of progenitor cells may occur very early in the development process before robust heterotypic interfacial tensions are established. Our findings are relevant to understanding the mechanisms that drive cell sorting in confluent tissues.

Exploiting differential effects of actomyosin contractility to control cell sorting among breast cancer cells

In order to gain a greater understanding of the factors that drive spatial organization in multicellular aggregates of cancer cells, we investigate the segregation patterns of 6 breast cell lines of varying degree of mesenchymal character during formation of mixed aggregates. Cell sorting is considered in the context of available adhesion proteins and cellular contractility. It is found that the primary compaction mediator (cadherins or integrins) for a given cell type in isolation plays an important role in compaction speed, which in turn is the major factor dictating preference for interior or exterior position within mixed aggregates. In particular, cadherin-deficient, invasion-competent cells tend to position towards the outside of aggregates, facilitating access to extracellular matrix. Reducing actomyosin contractility is found to have a differential effect on spheroid formation depending on compaction mechanism. Inhibition of contractility has a significant stabilizing effect on cell-cell adhesions in integrin-driven aggregation and a mildly destabilizing effect in cadherin-based aggregation. This differential response is exploited to statically control aggregate organization and dynamically rearrange cells in pre-formed aggregates. Sequestration of invasive cells in the interior of spheroids provides a physical barrier that reduces invasion in three-dimensional culture, revealing a potential strategy for containment of invasive cell types.

Topological floppy modes in models of epithelial tissues

Recent advances in topological mechanics have revealed unusual phenomena such as topologically protected floppy modes and states of self-stress that are exponentially localized at boundaries and interfaces of mechanical networks. In this paper, we explore the topological mechanics of epithelial tissues, where the appearance of these boundary and interface modes could lead to localized soft or stressed spots and play a role in morphogenesis. We consider both a simple vertex model (VM) governed by an effective elastic energy and its generalization to an active tension network (ATN) which incorporates active adaptation of the cytoskeleton. By analyzing spatially periodic lattices at the Maxwell point of mechanical instability, we find topologically polarized phases with exponential localization of floppy modes and states of self-stress in the ATN when cells are allowed to become concave, but not in the VM.

Activity waves and freestanding vortices in populations of subcritical Quincke rollers

Virtually all of the many active matter systems studied so far are made of units (biofilaments, cells, colloidal particles, robots, animals, etc.) that move even when they are alone or isolated. Their collective properties continue to fascinate, and we now understand better how they are unique to the bulk transduction of energy into work. Here we demonstrate that systems in which isolated but potentially active particles do not move can exhibit specific and remarkable collective properties. Combining experiments, theory, and numerical simulations, we show that such subcritical active matter can be realized with Quincke rollers, that is, dielectric colloidal particles immersed in a conducting fluid subjected to a vertical DC electric field. Working below the threshold field value marking the onset of motion for a single colloid, we find fast activity waves, reminiscent of excitable systems, and stable, arbitrarily large self-standing vortices made of thousands of particles moving at the same speed. Our theoretical model accounts for these phenomena and shows how they can arise in the absence of confining boundaries and individual chirality. We argue that our findings imply that a faithful description of the collective properties of Quincke rollers need to consider the fluid surrounding particles.

Mechanical feedback and robustness of apical constrictions in Drosophila embryo ventral furrow formation

Formation of the ventral furrow in the Drosophila embryo relies on the apical constriction of cells in the ventral region to produce bending forces that drive tissue invagination. In our recent paper we observed that apical constrictions during the initial phase of ventral furrow formation produce elongated patterns of cellular constriction chains prior to invagination and argued that these are indicative of tensile stress feedback. Here, we quantitatively analyze the constriction patterns preceding ventral furrow formation and find that they are consistent with the predictions of our active-granular-fluid model of a monolayer of mechanically coupled stress-sensitive constricting particles. Our model shows that tensile feedback causes constriction chains to develop along underlying precursor tensile stress chains that gradually strengthen with subsequent cellular constrictions. As seen in both our model and available optogenetic experiments, this mechanism allows constriction chains to penetrate or circumvent zones of reduced cell contractility, thus increasing the robustness of ventral furrow formation to spatial variation of cell contractility by rescuing cellular constrictions in the disrupted regions.

Embryonic Tissues as Active Foams

The physical state of embryonic tissues emerges from non-equilibrium, collective interactions among constituent cells. Cellular jamming, rigidity transitions and characteristics of glassy dynamics have all been observed in multicellular systems, but it is unclear how cells control these emergent tissue states and transitions, including tissue fluidization. Combining computational and experimental methods, here we show that tissue fluidization in posterior zebrafish tissues is controlled by the stochastic dynamics of tensions at cell-cell contacts. We develop a computational framework that connects cell behavior to embryonic tissue dynamics, accounting for the presence of extracellular spaces, complex cell shapes and cortical tension dynamics. We predict that tissues are maximally rigid at the structural transition between confluent and non-confluent states, with actively-generated tension fluctuations controlling stress relaxation and tissue fluidization. By directly measuring strain and stress relaxation, as well as the dynamics of cell rearrangements, in elongating posterior zebrafish tissues, we show that tension fluctuations drive active cell rearrangements that fluidize the tissue. These results highlight a key role of non-equilibrium tension dynamics in developmental processes.

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.

Recent Advances in Cellular and Molecular Bioengineering for Building and Translation of Biological Systems

In January of 2020, the Biomedical Engineering Society (BMES)- Cellular and Molecular Bioengineering (CMBE) conference was held in Puerto Rico and themed “Vision 2020: Emerging Technologies to Elucidate the Rule of Life.” The annual BME-CMBE conference gathered worldwide leaders and discussed successes and challenges in engineering biological systems and their translation. The goal of this report is to present the research frontiers in this field and provide perspectives on successful engineering and translation towards the clinic. We hope that this report serves as a constructive guide in shaping the future of research and translation of engineered biological systems.

Circular swimming motility and disordered hyperuniform state in an algae system

Active matter comprises individually driven units that convert locally stored energy into mechanical motion. Interactions between driven units lead to a variety of nonequilibrium collective phenomena in active matter. One of such phenomena is anomalously large density fluctuations, which have been observed in both experiments and theories. Here we show that, on the contrary, density fluctuations in active matter can also be greatly suppressed. Our experiments are carried out with marine algae ([Formula: see text]), which swim in circles at the air-liquid interfaces with two different eukaryotic flagella. Cell swimming generates fluid flow that leads to effective repulsions between cells in the far field. The long-range nature of such repulsive interactions suppresses density fluctuations and generates disordered hyperuniform states under a wide range of density conditions. Emergence of hyperuniformity and associated scaling exponent are quantitatively reproduced in a numerical model whose main ingredients are effective hydrodynamic interactions and uncorrelated random cell motion. Our results demonstrate the existence of disordered hyperuniform states in active matter and suggest the possibility of using hydrodynamic flow for self-assembly in active matter.

Spatiotemporal Oscillation in Confined Epithelial Motion upon Fluid-to-Solid Transition

Fluid-to-solid phase transition in multicellular assembly is crucial in many developmental biological processes, such as embryogenesis and morphogenesis. However, biomechanical studies in this area are limited, and little is known about factors governing the transition and how cell behaviors are regulated. Due to different stresses present, cells could behave distinctively depending on the nature of tissue. Here we report a fluid-to-solid transition in geometrically confined multicellular assemblies. Under circular confinement, Madin-Darby canine kidney (MDCK) monolayers undergo spatiotemporally oscillatory motions that are strongly dependent on the confinement size and distance from the periphery of the monolayers. Nanomechanical mapping reveals that epithelial tensional stress and traction forces on the substrate are both dependent on confinement size. The oscillation pattern and cellular nanomechanics profile appear well correlated with stress fiber assembly and cell polarization. These experimental observations imply that the confinement size-dependent surface tension regulates actin fiber assembly, cellular force generation, and cell polarization. Our analyses further suggest a characteristic confinement size (approximates to MDCK's natural correlation length) below which surface tension is sufficiently high and triggers a fluid-to-solid transition of the monolayers. Our findings may shed light on the geometrical and nanomechanical control of tissue morphogenesis and growth.

Poisson-bracket formulation of the dynamics of fluids of deformable particles

Using the Poisson-bracket method, we derive continuum equations for a fluid of deformable particles in two dimensions. Particle shape is quantified in terms of two continuum fields: an anisotropy density field that captures the deformations of individual particles from regular shapes and a shape tensor density field that quantifies both particle elongation and nematic alignment of elongated shapes. We explicitly consider the example of a dense biological tissue as described by the Vertex model energy, where cell shape has been proposed as a structural order parameter for a liquid-solid transition. The hydrodynamic model of biological tissue proposed here captures the coupling of cell shape to flow and provides a starting point for modeling the rheology of dense tissue.

Cell cycle-dependent active stress drives epithelia remodeling

Epithelia have distinct cellular architectures which are established in development, reestablished after wounding, and maintained during tissue homeostasis despite cell turnover and mechanical perturbations. In turn, cell shape also controls tissue function as a regulator of cell differentiation, proliferation, and motility. Here, we investigate cell shape changes in a model epithelial monolayer. After the onset of confluence, cells continue to proliferate and change shape over time, eventually leading to a final architecture characterized by arrested motion and more regular cell shapes. Such monolayer remodeling is robust, with qualitatively similar evolution in cell shape and dynamics observed across disparate perturbations. Here, we quantify differences in monolayer remodeling guided by the active vertex model to identify underlying order parameters controlling epithelial architecture. When monolayers are formed atop an extracellular matrix with varied stiffness, we find the cell density at which motion arrests varies significantly, but the cell shape remains constant, consistent with the onset of tissue rigidity. In contrast, pharmacological perturbations can significantly alter the cell shape at which tissue dynamics are arrested, consistent with varied amounts of active stress within the tissue. Across all experimental conditions, the final cell shape is well correlated to the cell proliferation rate, and cell cycle inhibition immediately arrests cell motility. Finally, we demonstrate cell cycle variation in junctional tension as a source of active stress within the monolayer. Thus, the architecture and mechanics of epithelial tissue can arise from an interplay between cell mechanics and stresses arising from cell cycle dynamics.

Glassy dynamics of a binary Voronoi fluid: a mode-coupling analysis

The binary Voronoi mixture is a fluid model whose interactions are derived from the Voronoi-Laguerre tessellation of the configurations of the system. The resulting interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent coherent and incoherent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter λ ≈ 0.746 comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good qualitative agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob-Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.

Bronchoconstriction: a potential missing link in airway remodelling

In asthma, progressive structural changes of the airway wall are collectively termed airway remodelling. Despite its deleterious effect on lung function, airway remodelling is incompletely understood. As one of the important causes leading to airway remodelling, here we discuss the significance of mechanical forces that are produced in the narrowed airway during asthma exacerbation, as a driving force of airway remodelling. We cover in vitro, ex vivo and in vivo work in this field, and discuss up-to-date literature supporting the idea that bronchoconstriction may be the missing link in a comprehensive understanding of airway remodelling in asthma.

Regionalized tissue fluidization is required for epithelial gap closure during insect gastrulation

Many animal embryos pull and close an epithelial sheet around the ellipsoidal egg surface during a gastrulation process known as epiboly. The ovoidal geometry dictates that the epithelial sheet first expands and subsequently compacts. Moreover, the spreading epithelium is mechanically stressed and this stress needs to be released. Here we show that during extraembryonic tissue (serosa) epiboly in the insect Tribolium castaneum, the non-proliferative serosa becomes regionalized into a solid-like dorsal region with larger non-rearranging cells, and a more fluid-like ventral region surrounding the leading edge with smaller cells undergoing intercalations. Our results suggest that a heterogeneous actomyosin cable contributes to the fluidization of the leading edge by driving sequential eviction and intercalation of individual cells away from the serosa margin. Since this developmental solution utilized during epiboly resembles the mechanism of wound healing, we propose actomyosin cable-driven local tissue fluidization as a conserved morphogenetic module for closure of epithelial gaps.

The mechanics of embryonic tissue spreading over spherical eggs is not fully understood. Here, the authors show that during gastrulation in the red flour beetle, extraembryonic tissue epiboly is facilitated by local actomyosin-mediated fluidization of the tissue at the leading edge.

Mesoscopic dynamic model of epithelial cell division with cell-cell junction effects

Cell division is central for embryonic development, tissue morphogenesis, and tumor growth. Experiments have evidenced that mitotic cell division is manipulated by the intercellular cues such as cell-cell junctions. However, it still remains unclear how these cortical-associated cues mechanically affect the mitotic spindle machinery, which determines the position and orientation of the cell division. In this paper, a mesoscopic dynamic cell division model is established to explore the integrated regulations of cortical polarity, microtubule pulling forces, cell deformability, and internal osmotic pressure. We show that the distributed pulling forces of astral microtubules play a key role in encoding the instructive cortical cues to orient and position the spindle of a dividing cell. The present model can not only predict the spindle orientation and position, but also capture the morphological evolution of cell rounding. The theoretical results agree well with relevant experiments both qualitatively and quantitatively. This work sheds light on the mechanical linkage between cell cortex and mitotic spindle, and holds potential in regulating cell division and sculpting tissue morphology.

Force networks, torque balance and Airy stress in the planar vertex model of a confluent epithelium

The vertex model is a popular framework for modelling tightly packed biological cells, such as confluent epithelia. Cells are described by convex polygons tiling the plane and their equilibrium is found by minimizing a global mechanical energy, with vertex locations treated as degrees of freedom. Drawing on analogies with granular materials, we describe the force network for a localized monolayer and derive the corresponding discrete Airy stress function, expressed for each N-sided cell as N scalars defined over kites covering the cell. We show how a torque balance (commonly overlooked in implementations of the vertex model) requires each internal vertex to lie at the orthocentre of the triangle formed by neighbouring edge centroids. Torque balance also places a geometric constraint on the stress in the neighbourhood of cellular trijunctions, and requires cell edges to be orthogonal to the links of a dual network that connect neighbouring cell centres and thereby triangulate the monolayer. We show how the Airy stress function depends on cell shape when a standard energy functional is adopted, and discuss implications for computational implementations of the model.

Control of cell colony growth by contact inhibition

Contact inhibition is a cell property that limits the migration and proliferation of cells in crowded environments. Here we investigate the growth dynamics of a cell colony composed of migrating and proliferating cells on a substrate using a minimal model that incorporates the mechanisms of contact inhibition of locomotion and proliferation. We find two distinct regimes. At early times, when contact inhibition is weak, the colony grows exponentially in time, fully characterised by the proliferation rate. At long times, the colony boundary moves at a constant speed, determined only by the migration speed of a single cell and independent of the proliferation rate. Further, the model demonstrates how cell-cell alignment speeds up colony growth. Our model illuminates how simple local mechanical interactions give rise to contact inhibition, and from this, how cell colony growth is self-organised and controlled on a local level.