Active Fingering Instability in Tissue Spreading

Interface morphodynamics in living tissues

Interfaces between distinct tissues or between tissues and environments are common in multicellular organisms. The evolution and stability of these interfaces are essential for tissue development, and their dysfunction can lead to diseases such as cancer. Mounting efforts, either theoretical or experimental, have been devoted to uncovering the morphodynamics of tissue interfaces. Here, we review the recent progress of studies on interface morphodynamics. The regulatory mechanisms governing interface evolution are dissected, with a focus on adhesion, cortical tension, cell activity, extracellular matrix, and microenvironment. We examine the methodologies used to study morphodynamics, emphasizing the characteristics of experimental techniques and theoretical models. Finally, we explore the broader implications of interface morphodynamics in tissue morphogenesis and diseases, offering a comprehensive perspective on this rapidly developing field.

Epithelial cell-cell interactions in an overcrowded environment: jamming or live cell extrusion

Epithelial tissues respond strongly to the mechanical stress caused by collective cell migration and are able to regulate it, which is important for biological processes such as morphogenesis, wound healing, and suppression of the spread of cancer. Compressive, tensional, and shear stress components are produced in cells when epithelial monolayers on substrate matrices are actively or passively wetted or de-wetted. Increased compressive stress on cells leads to enhanced cell-cell interactions by increasing the frequency of change the cell-cell distances, triggering various signalling pathways within the cells. This can ultimately lead either to cell jamming or to the extrusion of live cells. Despite extensive research in this field, it remains unclear how cells decide whether to jam, or to extrude a cell or cells, and how cells can reduce the compressive mechanical stress. Live cell extrusion from the overcrowded regions of the monolayers is associated with the presence of topological defects of cell alignment, induced by an interplay between the cell compressive and shear stress components. These topological defects stimulate cell re-alignment, as a part of the cells' tendency to re-establish an ordered trend of cell migration, by intensifying the glancing interactions in overcrowded regions. In addition to individual cell extrusion, collective cell extrusion has also been documented during monolayer active de-wetting, depending on the cell type, matrix stiffness, and boundary conditions. Cell jamming has been discussed in the context of the cells' contact inhibition of locomotion caused by cell head-on interactions. Since cell-cell interactions play a crucial role in cell rearrangement in an overcrowded environment, this review is focused on physical aspects of these interactions in order to stimulate further biological research in the field.

Segregation of co-cultured multicellular systems: review and modeling consideration

Cell segregation caused by collective cell migration (CCM) is crucial for morphogenesis, functional development of tissue parts, and is an important aspect in other diseases such as cancer and its metastasis process. Efficiency of the cell segregation depends on the interplay between: (1) biochemical processes such as cell signaling and gene expression and (2) physical interactions between cells. Despite extensive research devoted to study the segregation of various co-cultured systems, we still do not understand the role of physical interactions in cell segregation. Cumulative effects of these physical interactions appear in the form of physical parameters such as: (1) tissue surface tension, (2) viscoelasticity caused by CCM, and (3) solid stress accumulated in multicellular systems. These parameters primarily depend on the interplay between the state of cell-cell adhesion contacts and cell contractility. The role of these physical parameters on the segregation efficiency is discussed on model systems such as co-cultured breast cell spheroids consisting of two subpopulations that are in contact. This review study aims to: (1) summarize biological aspects related to cell segregation, mechanical properties of cell collectives, effects along the biointerface between cell subpopulations and (2) describe from a biophysical/mathematical perspective the same biological aspects summarized before. So that overall it can illustrate the complexity of the biological systems that translate into very complex biophysical/mathematical equations. Moreover, by presenting in parallel these two seemingly different parts (biology vs. equations), this review aims to emphasize the need for experiments to estimate the variety of parameters entering the resulting complex biophysical/mathematical models.

Tumor proliferation and invasion are intrinsically coupled and unraveled through tunable spheroid and physics-based models

Proliferation and invasion are two key drivers of tumor growth that are traditionally considered independent multicellular processes. However, these processes are intrinsically coupled through a maximum carrying capacity, i.e., the maximum spatial cell concentration supported by the tumor volume, total cell count, nutrient access, and mechanical properties of the tissue stroma. We explored this coupling of proliferation and invasion through in vitro and in silico methods where we modulated the mechanical properties of the tumor and the surrounding extracellular matrix. E-cadherin expression and stromal collagen concentration were manipulated in a tunable breast cancer spheroid to determine the overall impacts of these tumor variables on net tumor proliferation and continuum invasion. We integrated these results into a mixed-constitutive formulation to computationally delineate the influences of cellular and extracellular adhesion, stiffness, and mechanical properties of the extracellular matrix on net proliferation and continuum invasion. This framework integrates biological in vitro data into concise computational models of invasion and proliferation to provide more detailed physical insights into the coupling of these key tumor processes and tumor growth. STATEMENT OF SIGNIFICANCE: Tumor growth involves expansion into the collagen-rich stroma through intrinsic coupling of proliferation and invasion within the tumor continuum. These processes are regulated by a maximum carrying capacity that is determined by the total cell count, tumor volume, nutrient access, and mechanical properties of the surrounding stroma. The influences of biomechanical parameters (i.e., stiffness, cell elongation, net proliferation rate and cell-ECM friction) on tumor proliferation or invasion cannot be unraveled using experimental methods alone. By pairing a tunable spheroid system with computational modeling, we delineated the interdependencies of each system parameter on tumor proliferation and continuum invasion, and established a concise computational framework for studying tumor mechanobiology.

[The Dynamic Model of the Active-Inactive Cell Interface]

Objective

To explore the morphodynamics of the active-inactive cell monolayer interfaces by using the active liquid crystal model.

Methods

A continuum mechanical model was established based on the active liquid crystal theory and the active-inactive cell monolayer interfaces were established by setting the activity difference of cell monolayers. The theoretical equations were solved numerically by the finite difference and the lattice Boltzmann method.

Results

The active-inactive cell interfaces displayed three typical morphologies, namely, flat interface, wavy interface, and finger-like interface. On the flat interfaces, the cells were oriented perpendicular to the interface, the -1/2 topological defects were clustered in the interfaces, and the interfaces were negatively charged. On the wavy interfaces, cells showed no obvious preference for orientation at the interfaces and the interfaces were neutrally charged. On the finger-like interfaces, cells were tangentially oriented at the interfaces, the +1/2 topological defects were collected at the interfaces, driving the growth of the finger-like structures, and the interfaces were positively charged.

Conclusion

The orientation of the cell alignment at the interface can significantly affect the morphologies of the active-inactive cell monolayer interfaces, which is closely associated with the dynamics of topological defects at the interfaces.

Fingering Instability Accelerates Population Growth of a Proliferating Cell Collective

During the growth of a cell collective, such as proliferating microbial colonies and epithelial tissues, the local cell growth increases the local pressure, which in turn suppresses cell growth. How this pressure-growth coupling affects the growth of a cell collective remains unclear. Here, we answer this question using a continuum model of a cell collective. We find that a fast-growing leading front and a slow-growing interior of the cell collective emerge due to the pressure-dependent growth rate. The leading front can exhibit fingering instability, and we confirm the predicted instability criteria numerically with the leading front explicitly simulated. Intriguingly, we find that fingering instability is not only a consequence of local cell growth but also enhances the entire population's growth rate as positive feedback. Our work unveils the fitness advantage of fingering formation and suggests that the ability to form protrusions can be evolutionarily selected.

Cell jamming-to-unjamming transitions and vice versa in development: Physical aspects

Collective cell migration is essential for a wide range of biological processes such as: morphogenesis, wound healing, and cancer spreading. However, it is well known that migrating epithelial collectives frequently undergo jamming, stay trapped some period of time, and then start migration again. Consequently, only a part of epithelial cells actively contributes to the tissue development. In contrast to epithelial cells, migrating mesenchymal collectives successfully avoid the jamming. It has been confirmed that the epithelial unjamming cannot be treated as the epithelial-to-mesenchymal transition. Some other mechanism is responsible for the epithelial jamming/unjamming. Despite extensive research devoted to study the cell jamming/unjamming, we still do not understand the origin of this phenomenon. The origin is connected to physical factors such as: the cell compressive residual stress accumulation and surface characteristics of migrating (unjamming) and resting (jamming) epithelial clusters which depend primarily on the strength of cell-cell adhesion contacts and cell contractility. The main goal of this theoretical consideration is to clarify these cause-consequence relations.

Instabilities in a growing system of active particles: scalar and vectorial systems

The physics of micron-scale biological colonies usually benefits from different out-of-equilibrium sources. In bacterial colonies and cellular tissues, the growth process is among the important active sources that determine the dynamics. In this article, we study the generic dynamical instabilities associated with the growth phenomena that may arise in both scalar and vectorial systems. In vectorial systems, where the rotational degrees of particles play a role, a phenomenological growth-mediated torque can affect the rotational dynamics of individual particles. We show that such a growth-mediated torque can result in active traveling waves in the bulk of a growing system. In addition to the bulk properties, we analyze the instabilities in the shape of growing interfaces in both scalar and vectorial systems.

The analytical solution to the migration of an epithelial monolayer with a circular spreading front and its implications in the gap closure process

The coordinated behaviors of epithelial cells are widely observed in tissue development, such as re-epithelialization, tumor growth, and morphogenesis. In these processes, cells either migrate collectively or organize themselves into specific structures to serve certain purposes. In this work, we study a spreading epithelial monolayer whose migrating front encloses a circular gap in the monolayer center. Such tissue is usually used to mimic the wound healing process in vitro. We model the epithelial sheet as a layer of active viscous polar fluid. With an axisymmetric assumption, the model can be analytically solved under two special conditions, suggesting two possible spreading modes for the epithelial monolayer. Based on these two sets of analytical solutions, we assess the velocity of the spreading front affected by the gap size, the active intercellular contractility, and the purse-string contraction acting on the spreading edge. Several critical values exist in the model parameters for the initiation of the gap closure process, and the purse-string contraction plays a vital role in governing the gap closure kinetics. Finally, the instability of the morphology of the spreading front was studied. Numerical calculations show how the perturbated velocities and the growth rates vary with respect to different model parameters.

Emergence of Spatial Scales and Macroscopic Tissue Dynamics in Active Epithelial Monolayers

Migrating cells in tissues are often known to exhibit collective swirling movements. In this paper, we develop an active vertex model with polarity dynamics based on contact inhibition of locomotion (CIL). We show that under this dynamics, the cells form steady-state vortices in velocity, polarity, and cell stress with length scales that depend on polarity alignment rate (ζ), self-motility (v0), and cell-cell bond tension (λ). When the ratio λ/v0 becomes larger, the tissue reaches a near jamming state because of the inability of the cells to exchange their neighbors, and the length scale associated with tissue kinematics increases. A deeper examination of this jammed state provides insights into the mechanism of sustained swirl formation under CIL rule that is governed by the feedback between cell polarities and deformations. To gain additional understanding of how active forcing governed by CIL dynamics leads to large-scale tissue dynamics, we systematically coarse-grain cell stress, polarity, and motility and show that the tissue remains polar even on larger length scales. Overall, we explore the origin of swirling patterns during collective cell migration and obtain a connection between cell-level dynamics and large-scale cellular flow patterns observed in epithelial monolayers.

Size-dependent response of cells in epithelial tissue modulated by contractile stress fibers

Although cells with distinct apical areas have been widely observed in epithelial tissues, how the size of cells affects their behavior during tissue deformation and morphogenesis as well as key physical factors modulating such influence remains elusive. Here, we showed that the elongation of cells within the monolayer under anisotropic biaxial stretching increases with their size because the strain released by local cell rearrangement (i.e., T1 transition) is more significant for small cells that possess higher contractility. On the other hand, by incorporating the nucleation, peeling, merging, and breakage dynamics of subcellular stress fibers into classical vertex formulation, we found that stress fibers with orientations predominantly aligned with the main stretching direction will be formed at tricellular junctions, in good agreement with recent experiments. The contractile forces generated by stress fibers help cells to resist imposed stretching, reduce the occurrence of T1 transitions, and, consequently, modulate their size-dependent elongation. Our findings demonstrate that epithelial cells could utilize their size and internal structure to regulate their physical and related biological behaviors. The theoretical framework proposed here can also be extended to investigate the roles of cell geometry and intracellular contraction in processes such as collective cell migration and embryo development.

The role of viscoelasticity in long time cell rearrangement

Cell rearrangement caused by collective cell migration (CCM) during free expansion of epithelial monolayers has become a landmark in our current understanding of fundamental biological processes such as tissue development, regeneration, wound healing or cancer invasion. Cell spreading causes formation of mechanical waves which has a feedback effect on cell rearrangement and can lead to the cell jamming state. The mechanical waves describe oscillatory changes in cell velocity, as well as, the rheological parameters that affect them. The velocity oscillations, obtained at a time scale of hours, are in the form of forward and backward flows. Collision of forward and backward flows can induce an increase in the cell compressive stress accompanied with cell packing density which have a feedback impact on cell mobility, tissue viscoelasticity and alters the tissue stiffness. The tissue stiffness depends on the cell packing density and the active/passive (i.e. migrating/resting) state of single cells and can be used as an indicator of cell jamming state transition. Since cell stiffness can be measured it may directly show in which state the multicellular system is. In this work a review of existing modeling approaches is given along with assortment of published experimental findings, in order to invite experimentalists to test given theoretical considerations in multicellular systems.

A physical model for dynamic assembly of human salivary stem/progenitor microstructures

The in vitro reconstructions of human salivary glands in service of their eventual medical use represent a challenge for tissue engineering. Here, we present a theoretical approach to the dynamical formation of acinar structures from human salivary cells, focusing on observed stick-slip radial expansion as well as possible growth instabilities. Our findings demonstrate the critical importance of basement membrane remodeling in controlling the growth process.

Collective Cell Migration on Collagen-I Networks: The Impact of Matrix Viscoelasticity

Collective cell migration on extracellular matrix (ECM) networks is a key biological process involved in development, tissue homeostasis and diseases such as metastatic cancer. During invasion of epithelial cancers, cell clusters migrate through the surrounding stroma, which is comprised primarily of networks of collagen-I fibers. There is growing evidence that the rheological and topological properties of collagen networks can impact cell behavior and cell migration dynamics. During migration, cells exert mechanical forces on their substrate, resulting in an active remodeling of ECM networks that depends not only on the forces produced, but also on the molecular mechanisms that dictate network rheology. One aspect of collagen network rheology whose role is emerging as a crucial parameter in dictating cell behavior is network viscoelasticity. Dynamic reorganization of ECM networks can induce local changes in network organization and mechanics, which can further feed back on cell migration dynamics and cell-cell rearrangement. A number of studies, including many recent publications, have investigated the mechanisms underlying structural changes to collagen networks in response to mechanical force as well as the role of collagen rheology and topology in regulating cell behavior. In this mini-review, we explore the cause-consequence relationship between collagen network viscoelasticity and cell rearrangements at various spatiotemporal scales. We focus on structural alterations of collagen-I networks during collective cell migration and discuss the main rheological parameters, and in particular the role of viscoelasticity, which can contribute to local matrix stiffening during cell movement and can elicit changes in cell dynamics.

Bridging microscopic cell dynamics to nematohydrodynamics of cell monolayers

It is increasingly being realized that liquid-crystalline features can play an important role in the properties and dynamics of cell monolayers. Here, we present a cell-based model of cell layers, based on the phase-field formulation, that connects cell-cell interactions specified at the single cell level to large-scale nematic and hydrodynamic properties of the tissue. In particular, we present a minimal formulation that reproduces the well-known bend and splay hydrodynamic instabilities of the continuum nemato-hydrodynamic formulation of active matter, together with an analytical description of the instability threshold in terms of activity and elasticity of the cells. Furthermore, we provide a quantitative characterisation and comparison of flows and topological defects for extensile and contractile stress generation mechanisms, and demonstrate activity-induced heterogeneity and spontaneous formation of gaps within a confluent monolayer. Together, these results contribute to bridging the gap between cell-scale dynamics and tissue-scale collective cellular organisation.

Cellular Sensing Governs the Stability of Chemotactic Fronts

In contexts ranging from embryonic development to bacterial ecology, cell populations migrate chemotactically along self-generated chemical gradients, often forming a propagating front. Here, we theoretically show that the stability of such chemotactic fronts to morphological perturbations is determined by limitations in the ability of individual cells to sense and thereby respond to the chemical gradient. Specifically, cells at bulging parts of a front are exposed to a smaller gradient, which slows them down and promotes stability, but they also respond more strongly to the gradient, which speeds them up and promotes instability. We predict that this competition leads to chemotactic fingering when sensing is limited at too low chemical concentrations. Guided by this finding and by experimental data on E. coli chemotaxis, we suggest that the cells' sensory machinery might have evolved to avoid these limitations and ensure stable front propagation. Finally, as sensing of any stimuli is necessarily limited in living and active matter in general, the principle of sensing-induced stability may operate in other types of directed migration such as durotaxis, electrotaxis, and phototaxis.

Active T1 transitions in cellular networks

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

Chemotactic smoothing of collective migration

Collective migration-the directed, coordinated motion of many self-propelled agents-is a fascinating emergent behavior exhibited by active matter with functional implications for biological systems. However, how migration can persist when a population is confronted with perturbations is poorly understood. Here, we address this gap in knowledge through studies of bacteria that migrate via directed motion, or chemotaxis, in response to a self-generated nutrient gradient. We find that bacterial populations autonomously smooth out large-scale perturbations in their overall morphology, enabling the cells to continue to migrate together. This smoothing process arises from spatial variations in the ability of cells to sense and respond to the local nutrient gradient-revealing a population-scale consequence of the manner in which individual cells transduce external signals. Altogether, our work provides insights to predict, and potentially control, the collective migration and morphology of cellular populations and diverse other forms of active matter.

Effect of substrate stiffness on friction in collective cell migration

In collective cell migration, the motion results from forces produced by each cell and transmitted to the neighboring cells and to the substrate. Because inertia is negligible and the migration occurs over long time scales, the cell layer exhibits viscous behavior, where force and motion are connected by an apparent friction that results from the breaking and forming of adhesive bonds at the cell–cell and cell–substrate interfaces. Most theoretical models for collective migration include an apparent friction to connect force and motion, with many models making predictions that depend on the ratio of cell–cell and cell–substrate friction. However, little is known about factors that affect friction, leaving predictions of many theoretical models untested. Here, we considered how substrate stiffness and the number of adhesions affected friction at the cell–substrate interface. The experimental data were interpreted through prior theoretical models, which led to the same conclusion, that increased substrate stiffness increased the number of cell–substrate adhesions and caused increased cell–substrate friction. In turn, the friction affected the collective migration by altering the curvature at the edge of the cell layer. By revealing underlying factors affecting friction and demonstrating how friction perturbs the collective migration, this work provides experimental evidence supporting prior theoretical models and motivates the study of other ways to alter the collective migration by changing friction.

Collective durotaxis of cohesive cell clusters on a stiffness gradient

Many types of motile cells perform durotaxis, namely directed migration following gradients of substrate stiffness. Recent experiments have revealed that cell monolayers can migrate toward stiffer regions even when individual cells do not-a phenomenon known as collective durotaxis. Here, we address the spontaneous motion of finite cohesive cell monolayers on a stiffness gradient. We theoretically analyze a continuum active polar fluid model that has been tested in recent wetting assays of epithelial tissues and includes two types of active forces (cell-substrate traction and cell-cell contractility). The competition between the two active forces determines whether a cell monolayer spreads or contracts. Here, we show that this model generically predicts collective durotaxis, and that it features a variety of dynamical regimes as a result of the interplay between the spreading state and the global propagation, including sequential contraction and spreading of the monolayer as it moves toward higher stiffness. We solve the model exactly in some relevant cases, which provides both physical insights into the mechanisms of tissue durotaxis and spreading as well as a variety of predictions that could guide the design of future experiments.

Collective motility and mechanical waves in cell clusters

Epithelial cell clusters often move collectively on a substrate. Mechanical signals play a major role in organizing this behavior. There are a number of experimental observations in these systems which await a comprehensive explanation. These include: the internal strains are tensile even for clusters that expand by proliferation; the tractions on the substrate are often confined to the edges of the cluster; there can exist density waves within the cluster; and for cells in an annulus, there is a transition between expanding clusters with proliferation and the case where cells fill the annulus and rotate around it. We formulate a mechanical model to examine these effects. We use a molecular clutch picture which allows "stalling"-inhibition of cell contraction by external forces. Stalled cells are passive from a physical point of view and the un-stalled cells are active. By attaching cells to the substrate and to each other, and taking into account contact inhibition of locomotion, we get a simple picture for many of these findings as well as predictions that could be tested.

Fingering instability in spreading epithelial monolayers: roles of cell polarisation, substrate friction and contractile stresses

Collective cell migration plays a crucial role in many developmental processes that underlie morphogenesis, wound healing, or cancer progression. In such coordinated behaviours, cells are organised in coherent structures and actively migrate to serve different biological purposes. In some contexts, namely during epithelial wound healing, it is well known that a migrating free-edge monolayer develops finger-like instabilities, yet the onset is still under debate. Here, by means of theory and numerical simulations, we shed light on the main mechanisms driving the instability process, analysing the linear and nonlinear dynamics of a continuum compressible polar fluid. In particular, we assess the role of cell polarisation, substrate friction, and contractile stresses. Linear theory shows that it is crucial to analyse the perturbation transient dynamics, since we unravel a plethora of crossovers between different exponential growth rates during the linear regime. Numerical simulations suggest that cell-substrate friction could be the mechanism responsible for the formation of complex finger-like structures at the edge, since it triggers secondary fingering instabilities and tip-splitting phenomena. Finally, we obtain a critical contractile stress that depends on cell-substrate friction and the initial-to-nematic length ratio, characterising an active wetting-dewetting transition. In the dewetting scenario, the monolayer retracts and becomes stable without developing finger-like structures.

Force-mediated cellular anisotropy and plasticity dictate the elongation dynamics of embryos

We developed a unified dynamic model to explain how cellular anisotropy and plasticity, induced by alignment and severing/rebundling of actin filaments, dictate the elongation dynamics of Caenorhabditis elegans embryos. It was found that the gradual alignment of F-actins must be synchronized with the development of intracellular forces for the embryo to elongate, which is then further sustained by muscle contraction-triggered plastic deformation of cells. In addition, we showed that preestablished anisotropy is essential for the proper onset of the process while defects in the integrity or bundling kinetics of actin bundles result in abnormal embryo elongation, all in good agreement with experimental observations.

Mechanism underlying dynamic scaling properties observed in the contour of spreading epithelial monolayer

We found evidence of dynamic scaling in the spreading of Madin-Darby canine kidney (MDCK) cell monolayer, which can be characterized by the Hurst exponent α=0.86 and the growth exponent β=0.73, and theoretically and experimentally clarified the mechanism that governs the contour shape dynamics. Dynamic scaling refers to the roughness of the surface scales, both spatially and temporally. During the spreading of the monolayer, it is known that so-called leader cells generate the driving force and lead the other cells. Our time-lapse observations of cell behavior showed that these leader cells appeared at the early stage of the spreading and formed the monolayer protrusion. Informed by these observations, we developed a simple mathematical model that included differences in cell motility, cell-cell adhesion, and random cell movement. The model reproduced the quantitative characteristics obtained from the experiment, such as the spreading speed, the distribution of the increment, and the dynamic scaling law. Analysis of the model equation shows that the model can reproduce different scaling laws from (α=0.5,β=0.25) to (α=0.9,β=0.75), where the exponents α and β are determined by two dimensionless quantities determined by the microscopic cell behavior. From the analytical result, parameter estimation from the experimental results was achieved. The monolayer on the collagen-coated dishes showed a different scaling law, α=0.74,β=0.68, suggesting that cell motility increased ninefold. This result was consistent with the assay of the single-cell motility. Our study demonstrated that the dynamics of the contour of the monolayer were explained by the simple model, and we propose a mechanism that exhibits the dynamic scaling property.

Tracking collective cell motion by topological data analysis

By modifying and calibrating an active vertex model to experiments, we have simulated numerically a confluent cellular monolayer spreading on an empty space and the collision of two monolayers of different cells in an antagonistic migration assay. Cells are subject to inertial forces and to active forces that try to align their velocities with those of neighboring ones. In agreement with experiments in the literature, the spreading test exhibits formation of fingers in the moving interfaces, there appear swirls in the velocity field, and the polar order parameter and the correlation and swirl lengths increase with time. Numerical simulations show that cells inside the tissue have smaller area than those at the interface, which has been observed in recent experiments. In the antagonistic migration assay, a population of fluidlike Ras cells invades a population of wild type solidlike cells having shape parameters above and below the geometric critical value, respectively. Cell mixing or segregation depends on the junction tensions between different cells. We reproduce the experimentally observed antagonistic migration assays by assuming that a fraction of cells favor mixing, the others segregation, and that these cells are randomly distributed in space. To characterize and compare the structure of interfaces between cell types or of interfaces of spreading cellular monolayers in an automatic manner, we apply topological data analysis to experimental data and to results of our numerical simulations. We use time series of data generated by numerical simulations to automatically group, track and classify the advancing interfaces of cellular aggregates by means of bottleneck or Wasserstein distances of persistent homologies. These techniques of topological data analysis are scalable and could be used in studies involving large amounts of data. Besides applications to wound healing and metastatic cancer, these studies are relevant for tissue engineering, biological effects of materials, tissue and organ regeneration.

Reversible elastic phase field approach and application to cell monolayers

Motion and generation of forces by single cells and cell collectives are essential elements of many biological processes, including development, wound healing and cancer cell migration. Quantitative wound healing assays have demonstrated that cell monolayers can be both dynamic and elastic at the same time. However, it is very challenging to model this combination with conventional approaches. Here we introduce an elastic phase field approach that allows us to predict the dynamics of elastic sheets under the action of active stresses and localized forces, e.g. from leader cells. Our method ensures elastic reversibility after release of forces. We demonstrate its potential by studying several paradigmatic situations and geometries relevant for single cells and cell monolayers, including elastic bars, contractile discs and expanding monolayers with leader cells.

Size-dependent patterns of cell proliferation and migration in freely-expanding epithelia

The coordination of cell proliferation and migration in growing tissues is crucial in development and regeneration but remains poorly understood. Here, we find that, while expanding with an edge speed independent of initial conditions, millimeter-scale epithelial monolayers exhibit internal patterns of proliferation and migration that depend not on the current but on the initial tissue size, indicating memory effects. Specifically, the core of large tissues becomes very dense, almost quiescent, and ceases cell-cycle progression. In contrast, initially-smaller tissues develop a local minimum of cell density and a tissue-spanning vortex. To explain vortex formation, we propose an active polar fluid model with a feedback between cell polarization and tissue flow. Taken together, our findings suggest that expanding epithelia decouple their internal and edge regions, which enables robust expansion dynamics despite the presence of size- and history-dependent patterns in the tissue interior.

Criterion for Fingering Instabilities in Colloidal Gels

We sandwich a colloidal gel between two parallel plates and induce a radial flow by lifting the upper plate at a constant velocity. Two distinct scenarios result from such a tensile test: (i) stable flows during which the gel undergoes a tensile deformation without yielding, and (ii) unstable flows characterized by the radial growth of air fingers into the gel. We show that the unstable regime occurs beyond a critical energy input, independent of the gel's macroscopic yield stress. This implies a local fluidization of the gel at the tip of the growing fingers and results in the most unstable wavelength of the patterns exhibiting the characteristic scalings of the classical viscous fingering instability. Our work provides a quantitative criterion for the onset of fingering in colloidal gels based on a local shear-induced yielding in agreement with the delayed failure framework.

Relationship between velocities, tractions, and intercellular stresses in the migrating epithelial monolayer

The relationship between velocities, tractions, and intercellular stresses in the migrating epithelial monolayer are currently unknown. Ten years ago, a method known as monolayer stress microscopy (MSM) was suggested from which intercellular stresses could be computed for a given traction field. The core assumption of MSM is that intercellular stresses within the monolayer obey a linear and passive constitutive law. Examples of these include a Hookean solid (an elastic sheet) or a Newtonian fluid (thin fluid film), which imply a specific relation between the displacements or velocities and the tractions. Due to the lack of independently measured intercellular stresses, a direct validation of the 2D stresses predicted by a linear passive MSM model is presently not possible. An alternative approach, which we give here and denote as the Stokes method, is based on simultaneous measurements of the monolayer velocity field and the cell-substrate tractions. Using the same assumptions as those underlying MSM, namely, a linear and passive constitutive law, the velocity field suffices to compute tractions, from which we can then compare with those measured by traction force microscopy. We find that the calculated tractions and measured tractions are uncorrelated. Since the classical MSM and the Stokes approach both depend on the linear and passive constitutive law, it follows that some serious modification of the underling rheology is needed. One possible modification is the inclusion of an active force. In the special case where this is additive to the linear passive rheology, we have a new relationship between the active force density and the measured velocity (or displacement) field and tractions, which by Newton's laws, must be obeyed.