Collective cell migration: a physics perspective

Learning dynamical models of single and collective cell migration: a review

Single and collective cell migration are fundamental processes critical for physiological phenomena ranging from embryonic development and immune response to wound healing and cancer metastasis. To understand cell migration from a physical perspective, a broad variety of models for the underlying physical mechanisms that govern cell motility have been developed. A key challenge in the development of such models is how to connect them to experimental observations, which often exhibit complex stochastic behaviours. In this review, we discuss recent advances in data-driven theoretical approaches that directly connect with experimental data to infer dynamical models of stochastic cell migration. Leveraging advances in nanofabrication, image analysis, and tracking technology, experimental studies now provide unprecedented large datasets on cellular dynamics. In parallel, theoretical efforts have been directed towards integrating such datasets into physical models from the single cell to the tissue scale with the aim of conceptualising the emergent behaviour of cells. We first review how this inference problem has been addressed in both freely migrating and confined cells. Next, we discuss why these dynamics typically take the form of underdamped stochastic equations of motion, and how such equations can be inferred from data. We then review applications of data-driven inference and machine learning approaches to heterogeneity in cell behaviour, subcellular degrees of freedom, and to the collective dynamics of multicellular systems. Across these applications, we emphasise how data-driven methods can be integrated with physical active matter models of migrating cells, and help reveal how underlying molecular mechanisms control cell behaviour. Together, these data-driven approaches are a promising avenue for building physical models of cell migration directly from experimental data, and for providing conceptual links between different length-scales of description.

PIEZO1 regulates leader cell formation and cellular coordination during collective keratinocyte migration

The collective migration of keratinocytes during wound healing requires both the generation and transmission of mechanical forces for individual cellular locomotion and the coordination of movement across cells. Leader cells along the wound edge transmit mechanical and biochemical cues to ensuing follower cells, ensuring their coordinated direction of migration across multiple cells. Despite the observed importance of mechanical cues in leader cell formation and in controlling coordinated directionality of cell migration, the underlying biophysical mechanisms remain elusive. The mechanically-activated ion channel PIEZO1 was recently identified to play an inhibitory role during the reepithelialization of wounds. Here, through an integrative experimental and mathematical modeling approach, we elucidate PIEZO1's contributions to collective migration. Time-lapse microscopy reveals that PIEZO1 activity inhibits leader cell formation at the wound edge. To probe the relationship between PIEZO1 activity, leader cell formation and inhibition of reepithelialization, we developed an integrative 2D continuum model of wound closure that links observations at the single cell and collective cell migration scales. Through numerical simulations and subsequent experimental validation, we found that coordinated directionality plays a key role during wound closure and is inhibited by upregulated PIEZO1 activity. We propose that PIEZO1-mediated retraction suppresses leader cell formation which inhibits coordinated directionality between cells during collective migration.

ECM-derived biomaterials for regulating tissue multicellularity and maturation

Recent breakthroughs in developing human-relevant organotypic models led to the building of highly resemblant tissue constructs that hold immense potential for transplantation, drug screening, and disease modeling. Despite the progress in fine-tuning stem cell multilineage differentiation in highly controlled spatiotemporal conditions and hosting microenvironments, 3D models still experience naive and incomplete morphogenesis. In particular, existing systems and induction protocols fail to maintain stem cell long-term potency, induce high tissue-level multicellularity, or drive the maturity of stem cell-derived 3D models to levels seen in their in vivo counterparts. In this review, we highlight the use of extracellular matrix (ECM)-derived biomaterials in providing stem cell niche-mimicking microenvironment capable of preserving stem cell long-term potency and inducing spatial and region-specific differentiation. We also examine the maturation of different 3D models, including organoids, encapsulated in ECM biomaterials and provide looking-forward perspectives on employing ECM biomaterials in building more innovative, transplantable, and functional organs.

Residual cells and nutrient availability guide wound healing in bacterial biofilms

Biofilms are multicellular heterogeneous bacterial communities characterized by social-like division of labor, and remarkable robustness with respect to external stresses. Increasingly often an analogy between biofilms and arguably more complex eukaryotic tissues is being drawn. One illustrative example of where this analogy can be practically useful is the process of wound healing. While it has been extensively studied in eukaryotic tissues, the mechanism of wound healing in biofilms is virtually unexplored. Combining experiments in Bacillus subtilis bacteria, a model organism for biofilm formation, and a lattice-based theoretical model of biofilm growth, we studied how biofilms recover after macroscopic damage. We suggest that nutrient gradients and the abundance of proliferating cells are key factors augmenting wound closure. Accordingly, in the model, cell quiescence, nutrient fluxes, and biomass represented by cells and self-secreted extracellular matrix are necessary to qualitatively recapitulate the experimental results for damage repair. One of the surprising experimental findings is that residual cells, persisting in a damaged area after removal of a part of the biofilm, prominently affect the healing process. Taken together, our results outline the important roles of nutrient gradients and residual cells on biomass regrowth on macroscopic scales of the whole biofilm. The proposed combined experiment-simulation framework opens the way to further investigate the possible relation between wound healing, cell signaling and cell phenotype alternation in the local microenvironment of the wound.

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.

Polarization and motility of one-dimensional multi-cellular trains

Collective cell migration, whereby cells adhere to form multi-cellular clusters that move as a single entity, play an important role in numerous biological processes, such as during development and cancer progression. Recent experimental work focused on migration of one-dimensional cellular clusters, confined to move along adhesive lanes, as a simple geometry in which to systematically study this complex system. One-dimensional migration also arises in the body when cells migrate along blood vessels, axonal projections, and narrow cavities between tissues. We explore here the modes of one-dimensional migration of cellular clusters ("trains") by implementing cell-cell interactions in a model of cell migration that contains a mechanism for spontaneous cell polarization. We go beyond simple phenomenological models of the cells as self-propelled particles by having the internal polarization of each cell depend on its interactions with the neighboring cells that directly affect the actin polymerization activity at the cell's leading edges. Both contact inhibition of locomotion and cryptic lamellipodia interactions between neighboring cells are introduced. We find that this model predicts multiple motility modes of the cell trains, which can have several different speeds for the same polarization pattern. Compared to experimental data, we find that Madin-Darby canine kidney cells are poised along the transition region where contact inhibition of locomotion and cryptic lamellipodia roughly balance each other, where collective migration speed is most sensitive to the values of the cell-cell interaction strength.

Soliton approximation in continuum models of leader-follower behavior

Complex biological processes involve collective behavior of entities (bacteria, cells, animals) over many length and time scales and can be described by discrete models that track individuals or by continuum models involving densities and fields. We consider hybrid stochastic agent-based models of branching morphogenesis and angiogenesis (new blood vessel creation from preexisting vasculature), which treat cells as individuals that are guided by underlying continuous chemical and/or mechanical fields. In these descriptions, leader (tip) cells emerge from existing branches and follower (stalk) cells build the new sprout in their wake. Vessel branching and fusion (anastomosis) occur as a result of tip and stalk cell dynamics. Coarse graining these hybrid models in appropriate limits produces continuum partial differential equations (PDEs) for endothelial cell densities that are more analytically tractable. While these models differ in nonlinearity, they produce similar equations at leading order when chemotaxis is dominant. We analyze this leading order system in a simple quasi-one-dimensional geometry and show that the numerical solution of the leading order PDE is well described by a soliton wave that evolves from vessel to source. This wave is an attractor for intermediate times until it arrives at the hypoxic region releasing the growth factor. The mathematical techniques used here thus identify common features of discrete and continuum approaches and provide insight into general biological mechanisms governing their collective dynamics.

Coarse-graining the vertex model and its response to shear

Tissue dynamics and collective cell motion are crucial biological processes. Their biological machinery is mostly known, and simulation models such as the active vertex model exist and yield reasonable agreement with experimental observations such as tissue fluidization or fingering. However, a good and well-founded continuum description for tissues remains to be developed. In this work, we derive a macroscopic description for a two-dimensional cell monolayer by coarse-graining the vertex model through the Poisson bracket approach. We obtain equations for cell density, velocity, and the cellular shape tensor. We then study the homogeneous steady states, their stability (which coincides with thermodynamic stability), and especially their behavior under an externally applied shear. Our results contribute to elucidate the interplay between flow and cellular shape. The obtained macroscopic equations present a good starting point for adding cell motion, morphogenetic, and other biologically relevant processes.

A planar polarized MYO6-DOCK7-RAC1 axis promotes tissue fluidification in mammary epithelia

Tissue fluidification and collective motility are pivotal in regulating embryonic morphogenesis, wound healing, and tumor metastasis. These processes frequently require that each cell constituent of a tissue coordinates its migration activity and directed motion through the oriented extension of lamellipodium cell protrusions, promoted by RAC1 activity. While the upstream RAC1 regulators in individual migratory cells or leader cells during invasion or wound healing are well characterized, how RAC1 is controlled in follower cells remains unknown. Here, we identify a MYO6-DOCK7 axis essential for spatially restricting RAC1 activity in a planar polarized fashion in model tissue monolayers. The MYO6-DOCK7 axis specifically controls the extension of cryptic lamellipodia required to drive tissue fluidification and cooperative-mode motion in otherwise solid and static carcinoma cell collectives.

Microfluidic platform for the reproduction of hypoxic vascular microenvironments

Vascular endothelial cells (ECs) respond to mechanical stimuli caused by blood flow to maintain vascular homeostasis. Although the oxygen level in vascular microenvironment is lower than the atmospheric one, the cellular dynamics of ECs under hypoxic and flow exposure are not fully understood. Here, we describe a microfluidic platform for the reproduction hypoxic vascular microenvironments. Simultaneous application of hypoxic stress and fluid shear stress to the cultured cells was achieved by integrating a microfluidic device and a flow channel that adjusted the initial oxygen concentration in a cell culture medium. An EC monolayer was then formed on the media channel in the device, and the ECs were observed after exposure to hypoxic and flow conditions. The migration velocity of the ECs immediately increased after flow exposure, especially in the direction opposite to the flow direction, and gradually decreased, resulting in the lowest value under the hypoxic and flow exposure condition. The ECs after 6-h simultaneous exposure to hypoxic stress and fluid shear stress were generally aligned and elongated in the flow direction, with enhanced VE-cadherin expression and actin filament assembly. Thus, the developed microfluidic platform is useful for investigating the dynamics of ECs in vascular microenvironments.

3D printing of conch-like scaffolds for guiding cell migration and directional bone growth

Regeneration of severe bone defects remains an enormous challenge in clinic. Developing regenerative scaffolds to directionally guide bone growth is a potential strategy to overcome this hurdle. Conch, an interesting creature widely spreading in ocean, has tough spiral shell that can continuously grow along the spiral direction. Herein, inspired by the physiological features of conches, a conch-like (CL) scaffold based on β-TCP bioceramic material was successfully prepared for guiding directional bone growth via digital light processing (DLP)-based 3D printing. Benefiting from the spiral structure, the CL scaffolds significantly improved cell adhesion, proliferation and osteogenic differentiation in vitro compared to the conventional 3D scaffolds. Particularly, the spiral structure in the scaffolds could efficiently induce cells to migrate from the bottom to the top of the scaffolds, which was like “cells climbing stairs”. Furthermore, the capability of guiding directional bone growth for the CL scaffolds was demonstrated by a special half-embedded femoral defects model in rabbits. The new bone tissue could consecutively grow into the protruded part of the scaffolds along the spiral cavities. This work provides a promising strategy to construct biomimetic biomaterials for guiding directional bone tissue growth, which offers a new treatment concept for severe bone defects, and even limb regeneration.

•

A conch-like scaffold was firstly developed for guiding directional bone growth.

•

The CL scaffolds efficiently induced cells “climbing stairs”- like-migrating.

•

The CL scaffolds showed improved bioactivities benefited from the spiral structure.

•

This work provided a new treatment concept for severe bone defects.

Liquid-hexatic-solid phases in active and passive Brownian particles determined by stochastic birth and death events

We study the effects of stochastic birth and death processes on the structural phases of systems of active and passive Brownian particles subject to volume exclusion. The total number of particles in the system is a fluctuating quantity, determined by the birth and death parameters and on the activity of the particles. As the birth and death parameters are varied, we find liquid, hexatic, and solid phases. For passive particles, these phases are found to be spatially homogeneous. For active particles, motility-induced phase separation (coexisting hexatic and liquid phases) occurs for large activity and sufficiently small birth rates. We also observe a reentrant transition to the hexatic phase when the birth rate is increased. This results from a balance of an increasing number of particles filling the system, and a larger number of defects resulting from the birth and death dynamics.

Unified multiscale theory of cellular mechanical adaptations to substrate stiffness

Rigidity of the extracellular matrix markedly regulates many cellular processes. However, how cells detect and respond to matrix rigidity remains incompletely understood. Here, we propose a unified two-dimensional multiscale framework accounting for the chemomechanical feedback to explore the interrelated cellular mechanosensing, polarization, and migration, which constitute the dynamic cascade in cellular response to matrix stiffness but are often modeled separately in previous theories. By combining integrin dynamics and intracellular force transduction, we show that substrate stiffness can act as a switch to activate or deactivate cell polarization. Our theory quantitatively reproduces rich stiffness-dependent cellular dynamics, including spreading, polarity selection, migration pattern, durotaxis, and even negative durotaxis, reported in a wide spectrum of cell types, and reconciles some inconsistent experimental observations. We find that a specific bipolarized mode can determine the optimal substrate stiffness, which enables the fastest cell migration rather than the largest traction forces that cells apply on the substrate. We identify that such a mechanical adaptation stems from the force balance across the whole cell. These findings could yield universal insights into various stiffness-mediated cellular processes within the context of tissue morphogenesis, wound healing, and cancer invasion.

Modelling the modulation of cortical Up-Down state switching by astrocytes

Up-Down synchronization in neuronal networks refers to spontaneous switches between periods of high collective firing activity (Up state) and periods of silence (Down state). Recent experimental reports have shown that astrocytes can control the emergence of such Up-Down regimes in neural networks, although the molecular or cellular mechanisms that are involved are still uncertain. Here we propose neural network models made of three populations of cells: excitatory neurons, inhibitory neurons and astrocytes, interconnected by synaptic and gliotransmission events, to explore how astrocytes can control this phenomenon. The presence of astrocytes in the models is indeed observed to promote the emergence of Up-Down regimes with realistic characteristics. Our models show that the difference of signalling timescales between astrocytes and neurons (seconds versus milliseconds) can induce a regime where the frequency of gliotransmission events released by the astrocytes does not synchronize with the Up and Down phases of the neurons, but remains essentially stable. However, these gliotransmission events are found to change the localization of the bifurcations in the parameter space so that with the addition of astrocytes, the network enters a bistability region of the dynamics that corresponds to Up-Down synchronization. Taken together, our work provides a theoretical framework to test scenarios and hypotheses on the modulation of Up-Down dynamics by gliotransmission from astrocytes.

Neural networks in many brain regions can display synchronized activities. During the so-called “Up-Down” synchronization regimes for instance, the whole local population of neurons switches in a spontaneous and synchronized fashion between phases of high activity (Up states) and phases of low activity (Down states). The mechanisms responsible for this behaviour are still not well understood, but recent experimental reports have suggested that another type of brain cells, the astrocytes, at least partly control these oscillations. Astrocytes are increasingly believed to play a role in the propagation of signals between neurons, via their connections to neuronal synapses, but how this mechanism could control Up-Down regimes is not understood. To address this issue we present here simple mathematical models of neuronal networks that incorporate astrocytes in addition to neurons according to various levels of description. Using bifurcation analysis and numerical simulations we explore how astrocytes control Up-Down synchronization of the neuronal networks. In particular, astrocytes in the model are found to change the localization of the bifurcation points in the parameter space, so that the neurons enter the region of Up-Down regime when astrocytes are present. We also give some theoretical predictions that can be tested experimentally to test the validity of our models.

Interplay between cell height variations and planar pulsations in epithelial monolayers

Biological tissues change their shapes through collective interactions of cells. This coordination sets length and time scales for dynamics where precision is essential, in particular during morphogenetic events. However, how these scales emerge remains unclear. Here, we address this question using the pulsatile domains observed in confluent epithelial MDCK monolayers where cells exhibit synchronous contraction and extension cycles of [Formula: see text] h duration and [Formula: see text] length scale. We report that the monolayer thickness changes gradually in space and time by more than twofold in order to counterbalance the contraction and extension of the incompressible cytoplasm. We recapitulate these pulsatile dynamics using a continuum model and show that incorporation of cell stiffness dependent height variations is critical both for generating temporal pulsations and establishing the domain size. We propose that this feedback between height and mechanics could be important in coordinating the length scales of tissue dynamics.

Role of Delta-Notch signalling molecules on cell-cell adhesion in determining heterogeneous chemical and cell morphological patterning

Cell mechanics and motility are responsible for collective motion of cells that result in overall deformation of epithelial tissues. On the other hand, contact-dependent cell-cell signalling is responsible for generating a large variety of intricate, self-organized, spatial patterns of the signalling molecules. Moreover, it is becoming increasingly clear that the combined mechanochemical patterns of cell shape/size and signalling molecules in the tissues, for example, in cancerous and sensory epithelium, are governed by mechanochemical coupling between chemical signalling and cell mechanics. However, a clear quantitative picture of how these two aspects of tissue dynamics, i.e., signalling and mechanics, lead to pattern and form is still emerging. Although, a number of recent experiments demonstrate that cell mechanics, cell motility, and cell-cell signalling are tightly coupled in many morphogenetic processes, relatively few modeling efforts have focused on an integrated approach. We extend the vertex model of an epithelial monolayer to account for contact-dependent signalling between adjacent cells and between non-adjacent neighbors through long protrusional contacts with a feedback mechanism wherein the adhesive strength between adjacent cells is controlled by the expression of the signalling molecules in those cells. Local changes in cell-cell adhesion lead to changes in cell shape and size, which in turn drives changes in the levels of signalling molecules. Our simulations show that even this elementary two-way coupling of chemical signalling and cell mechanics is capable of giving rise to a rich variety of mechanochemical patterns in epithelial tissues. In particular, under certain parametric conditions, bimodal distributions in cell size and shape are obtained, which resemble experimental observations in cancerous and sensory tissues.

On the origins of order

A cardinal feature common to embryonic development and tissue reorganization, as well as to wound healing and cancer cell invasion, is collective cellular migration. During collective migratory events the phenomena of cell jamming and unjamming are increasingly recognized, and underlying mechanical, genomic, transcriptional, and signaling events are increasingly coming to light. In this brief perspective I propose a synthesis that brings together in a new way two key concepts. On the one hand, it has been suggested that the unjammed phase of the cellular collective evolved under a selective pressure favoring fluid-like migratory dynamics as would be required so as to accommodate episodes of tissue evolution, development, plasticity, and repair. Being dynamic, such an unjammed migratory phase is expected to be energetically expensive compared with the jammed non-migratory phase, which is presumed to have evolved under a selective pressure favoring a solid-like homeostatic regime that, by comparison, is energetically economical and mechanically stable. On the other hand, well before the discovery of cell jamming and unjamming Kauffman proposed the general biological principle that living systems exist in a solid regime near the edge of chaos, and that natural selection achieves and sustains such a poised state. Here I propose that, in certain systems at least, this poised solid-like state as predicted in the abstract by Kauffman is realized in the particular by the jammed regime just at the brink of unjamming.

The role of cell geometry and cell-cell communication in gradient sensing

Cells can measure shallow gradients of external signals to initiate and accomplish a migration or a morphogenetic process. Recently, starting from mathematical models like the local-excitation global-inhibition (LEGI) model and with the support of empirical evidence, it has been proposed that cellular communication improves the measurement of an external gradient. However, the mathematical models that have been used have over-simplified geometries (e.g., they are uni-dimensional) or assumptions about cellular communication, which limit the possibility to analyze the gradient sensing ability of more complex cellular systems. Here, we generalize the existing models to study the effects on gradient sensing of cell number, geometry and of long- versus short-range cellular communication in 2D systems representing epithelial tissues. We find that increasing the cell number can be detrimental for gradient sensing when the communication is weak and limited to nearest neighbour cells, while it is beneficial when there is long-range communication. We also find that, with long-range communication, the gradient sensing ability improves for tissues with more disordered geometries; on the other hand, an ordered structure with mostly hexagonal cells is advantageous with nearest neighbour communication. Our results considerably extend the current models of gradient sensing by epithelial tissues, making a step further toward predicting the mechanism of communication and its putative mediator in many biological processes.

Impact of contact inhibition on collective cell migration and proliferation

Contact inhibition limits migration and proliferation of cells in cell colonies. We consider a multiphase field model to investigate the growth dynamics of a cell colony, composed of proliferating cells. The model takes into account the mechanism of contact inhibition of proliferation by local mechanical interactions. We compare nonmigrating and migrating cells, in order to provide a quantitative characterization of the dynamics and analyze the velocity of the colony boundary for both cases. Additionally, we measure single cell velocities, number of neighbor distributions, as well as the influence of stress and age on positions of the cells and with respect to each other.

p53 directs leader cell behavior, migration, and clearance during epithelial repair

Epithelial cells migrate across wounds to repair injured tissue. Leader cells at the front of migrating sheets often drive this process. However, it is unclear how leaders emerge from an apparently homogeneous epithelial cell population. We characterized leaders emerging from epithelial monolayers in cell culture and found that they activated the stress sensor p53, which was sufficient to initiate leader cell behavior. p53 activated the cell cycle inhibitor p21WAF1/CIP1, which in turn induced leader behavior through inhibition of cyclin-dependent kinase activity. p53 also induced crowding hypersensitivity in leader cells such that, upon epithelial closure, they were eliminated by cell competition. Thus, mechanically induced p53 directs emergence of a transient population of leader cells that drive migration and ensures their clearance upon epithelial repair.

Sculpting tissues by phase transitions

Biological systems display a rich phenomenology of states that resemble the physical states of matter - solid, liquid and gas. These phases result from the interactions between the microscopic constituent components - the cells - that manifest in macroscopic properties such as fluidity, rigidity and resistance to changes in shape and volume. Looked at from such a perspective, phase transitions from a rigid to a flowing state or vice versa define much of what happens in many biological processes especially during early development and diseases such as cancer. Additionally, collectively moving confluent cells can also lead to kinematic phase transitions in biological systems similar to multi-particle systems where the particles can interact and show sub-populations characterised by specific velocities. In this Perspective we discuss the similarities and limitations of the analogy between biological and inert physical systems both from theoretical perspective as well as experimental evidence in biological systems. In understanding such transitions, it is crucial to acknowledge that the macroscopic properties of biological materials and their modifications result from the complex interplay between the microscopic properties of cells including growth or death, neighbour interactions and secretion of matrix, phenomena unique to biological systems. Detecting phase transitions in vivo is technically difficult. We present emerging approaches that address this challenge and may guide our understanding of the organization and macroscopic behaviour of biological tissues.

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.

Disentangling cadherin-mediated cell-cell interactions in collective cancer cell migration

Cell dispersion from a confined area is fundamental in a number of biological processes, including cancer metastasis. To date, a quantitative understanding of the interplay of single-cell motility, cell proliferation, and intercellular contacts remains elusive. In particular, the role of E- and N-cadherin junctions, central components of intercellular contacts, is still controversial. Combining theoretical modeling with in vitro observations, we investigate the collective spreading behavior of colonies of human cancer cells (T24). The spreading of these colonies is driven by stochastic single-cell migration with frequent transient cell-cell contacts. We find that inhibition of E- and N-cadherin junctions decreases colony spreading and average spreading velocities, without affecting the strength of correlations in spreading velocities of neighboring cells. Based on a biophysical simulation model for cell migration, we show that the behavioral changes upon disruption of these junctions can be explained by reduced repulsive excluded volume interactions between cells. This suggests that in cancer cell migration, cadherin-based intercellular contacts sharpen cell boundaries leading to repulsive rather than cohesive interactions between cells, thereby promoting efficient cell spreading during collective migration.

Multiphase field models for collective cell migration

Confluent cell monolayers and epithelia tissues show remarkable patterns and correlations in structural arrangements and actively driven collective flows. We simulate these properties using multiphase field models. The models are based on cell deformations and cell-cell interactions and we investigate the influence of microscopic details to incorporate active forces on emerging phenomena. We compare four different approaches, one in which the activity is determined by a random orientation, one where the activity is related to the deformation of the cells, and two models with subcellular details to resolve the mechanochemical interactions underlying cell migration. The models are compared with respect to generic features, such as coordination number distribution, cell shape variability, emerging nematic properties, as well as vorticity correlations and flow patterns in large confluent monolayers and confinements. All results are compared with experimental data for a large variety of cell cultures. The appearing qualitative differences of the models show the importance of microscopic details and provide a route towards predictive simulations of patterns and correlations in cell colonies.

Are cell jamming and unjamming essential in tissue development?

The last decade has seen a surge of evidence supporting the existence of the transition of the multicellular tissue from a collective material phase that is regarded as being jammed to a collective material phase that is regarded as being unjammed. The jammed phase is solid-like and effectively 'frozen', and therefore is associated with tissue homeostasis, rigidity, and mechanical stability. The unjammed phase, by contrast, is fluid-like and effectively 'melted', and therefore is associated with mechanical fluidity, plasticity and malleability that are required in dynamic multicellular processes that sculpt organ microstructure. Such multicellular sculpturing, for example, occurs during embryogenesis, growth and remodeling. Although unjamming and jamming events in the multicellular collective are reminiscent of those that occur in the inert granular collective, such as grain in a hopper that can flow or clog, the analogy is instructive but limited, and the implications for cell biology remain unclear. Here we ask, are the cellular jamming transition and its inverse --the unjamming transition-- mere epiphenomena? That is, are they dispensable downstream events that accompany but neither cause nor quench these core multicellular processes? Drawing from selected examples in developmental biology, here we suggest the hypothesis that, to the contrary, the graded departure from a jammed phase enables controlled degrees of malleability as might be required in developmental dynamics. We further suggest that the coordinated approach to a jammed phase progressively slows those dynamics and ultimately enables long-term mechanical stability as might be required in the mature homeostatic multicellular tissue.

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.

Permeation Instabilities in Active Polar Gels

We present a theory of active, permeating, polar gels, based on a two-fluid model. An active relative force between the gel components creates a steady-state current. We analyze its stability, while considering two polar coupling terms to the relative current: a permeation-deformation term, which describes network deformation by the solvent flow, and a permeation-alignment term, which describes the alignment of the polarization field by the network deformation and flow. Novel instability mechanisms emerge at finite wave vectors, suggesting the formation of periodic domains and mesophases. Our results can be used to determine the physical conditions required for various types of multicellular migration across tissues.

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.

Endocytosis in the context-dependent regulation of individual and collective cell properties

Endocytosis allows cells to transport particles and molecules across the plasma membrane. In addition, it is involved in the termination of signalling through receptor downmodulation and degradation. This traditional outlook has been substantially modified in recent years by discoveries that endocytosis and subsequent trafficking routes have a profound impact on the positive regulation and propagation of signals, being key for the spatiotemporal regulation of signal transmission in cells. Accordingly, endocytosis and membrane trafficking regulate virtually every aspect of cell physiology and are frequently subverted in pathological conditions. Two key aspects of endocytic control over signalling are coming into focus: context-dependency and long-range effects. First, endocytic-regulated outputs are not stereotyped but heavily dependent on the cell-specific regulation of endocytic networks. Second, endocytic regulation has an impact not only on individual cells but also on the behaviour of cellular collectives. Herein, we will discuss recent advancements in these areas, highlighting how endocytic trafficking impacts complex cell properties, including cell polarity and collective cell migration, and the relevance of these mechanisms to disease, in particular cancer.

Decoding leader cells in collective cancer invasion

Collective cancer invasion with leader-follower organization is increasingly recognized as a predominant mechanism in the metastatic cascade. Leader cells support cancer invasion by creating invasion tracks, sensing environmental cues and coordinating with follower cells biochemically and biomechanically. With the latest developments in experimental and computational models and analysis techniques, the range of specific traits and features of leader cells reported in the literature is rapidly expanding. Yet, despite their importance, there is no consensus on how leader cells arise or their essential characteristics. In this Perspective, we propose a framework for defining the essential aspects of leader cells and provide a unifying perspective on the varying cellular and molecular programmes that are adopted by each leader cell subtype to accomplish their functions. This Perspective can lead to more effective strategies to interdict a major contributor to metastatic capability.

Photothermal Agarose Microfabrication Technology for Collective Cell Migration Analysis

Agarose photothermal microfabrication technology is one of the micropatterning techniques that has the advantage of simple and flexible real-time fabrication even during the cultivation of cells. To examine the ability and limitation of the agarose microstructures, we investigated the collective epithelial cell migration behavior in two-dimensional agarose confined structures. Agarose microchannels from 10 to 211 micrometer width were fabricated with a spot heating of a focused 1480 nm wavelength infrared laser to the thin agarose layer coated on the cultivation dish after the cells occupied the reservoir. The collective cell migration velocity maintained constant regardless of their extension distance, whereas the width dependency of those velocities was maximized around 30 micrometer width and decreased both in the narrower and wider microchannels. The single-cell tracking revealed that the decrease of velocity in the narrower width was caused by the apparent increase of aspect ratio of cell shape (up to 8.9). In contrast, the decrease in the wider channels was mainly caused by the increase of the random walk-like behavior of component cells. The results confirmed the advantages of this method: (1) flexible fabrication without any pre-designing, (2) modification even during cultivation, and (3) the cells were confined in the agarose geometry.

Theoretical Constraints Imposed by Gradient Detection and Dispersal on Microbial Size in Astrobiological Environments

The capacity to sense gradients efficiently and acquire information about the ambient environment confers many advantages such as facilitating movement toward nutrient sources or away from toxic chemicals. The amplified dispersal evinced by organisms endowed with motility is possibly beneficial in related contexts. Hence, the connections between information acquisition, motility, and microbial size are explored from an explicitly astrobiological standpoint. By using prior theoretical models, the constraints on organism size imposed by gradient detection and motility are elucidated in the form of simple heuristic scaling relations. It is argued that environments such as alkaline hydrothermal vents, which are distinguished by the presence of steep gradients, might be conducive to the existence of "small" microbes (with radii of ≳0.1 μm) in principle, when only the above two factors are considered; other biological functions (e.g., metabolism and genetic exchange) could, however, regulate the lower bound on microbial size and elevate it. The derived expressions are potentially applicable to a diverse array of settings, including those entailing solvents other than water; for example, the lakes and seas of Titan. The article concludes with a brief exposition of how this formalism may be of practical and theoretical value to astrobiology.

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. The patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

Fluctuations can induce local nematic order and extensile stress in monolayers of motile cells

Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix.

Dynamically Re-Organized Collagen Fiber Bundles Transmit Mechanical Signals and Induce Strongly Correlated Cell Migration and Self-Organization

Correlated cell migration in fibrous extracellular matrix (ECM) is important in many biological processes. During migration, cells can remodel the ECM, leading to the formation of mesoscale structures such as fiber bundles. However, how such mesoscale structures regulate correlated single-cells migration remains to be elucidated. Here, using a quasi-3D in vitro model, we investigate how collagen fiber bundles are dynamically re-organized and guide cell migration. By combining laser ablation technique with 3D tracking and active-particle simulations, we definitively show that only the re-organized fiber bundles that carry significant tensile forces can guide strongly correlated cell migration, providing for the first time a direct experimental evidence supporting that matrix-transmitted long-range forces can regulate cell migration and self-organization. This force regulation mechanism can provide new insights for studies on cellular dynamics, fabrication or selection of biomedical materials in tissue repairing, and many other biomedical applications.

Statistical mechanics of active Ornstein-Uhlenbeck particles

We study the statistical properties of active Ornstein-Uhlenbeck particles (AOUPs). In this simplest of models, the Gaussian white noise of overdamped Brownian colloids is replaced by a Gaussian colored noise. This suffices to grant this system the hallmark properties of active matter, while still allowing for analytical progress. We study in detail the steady-state distribution of AOUPs in the small persistence time limit and for spatially varying activity. At the collective level, we show AOUPs to experience motility-induced phase separation both in the presence of pairwise forces or due to quorum-sensing interactions. We characterize both the instability mechanism leading to phase separation and the resulting phase coexistence. We probe how, in the stationary state, AOUPs depart from their thermal equilibrium limit by investigating the emergence of ratchet currents and entropy production. In the small persistence time limit, we show how fluctuation-dissipation relations are recovered. Finally, we discuss how the emerging properties of AOUPs can be characterized from the dynamics of their collective modes.

Study of Wound Healing Dynamics by Single Pseudo-Particle Tracking in Phase Contrast Images Acquired in Time-Lapse

Cellular contacts modify the way cells migrate in a cohesive group with respect to a free single cell. The resulting motion is persistent and correlated, with cells’ velocities self-aligning in time. The presence of a dense agglomerate of cells makes the application of single particle tracking techniques to define cells dynamics difficult, especially in the case of phase contrast images. Here, we propose an original pipeline for the analysis of phase contrast images of the wound healing scratch assay acquired in time-lapse, with the aim of extracting single particle trajectories describing the dynamics of the wound closure. In such an approach, the membrane of the cells at the border of the wound is taken as a unicum, i.e., the wound edge, and the dynamics is described by the stochastic motion of an ensemble of points on such a membrane, i.e., pseudo-particles. For each single frame, the pipeline of analysis includes: first, a texture classification for separating the background from the cells and for identifying the wound edge; second, the computation of the coordinates of the ensemble of pseudo-particles, chosen to be uniformly distributed along the length of the wound edge. We show the results of this method applied to a glioma cell line (T98G) performing a wound healing scratch assay without external stimuli. We discuss the efficiency of the method to assess cell motility and possible applications to other experimental layouts, such as single cell motion. The pipeline is developed in the Python language and is available upon request.

Role of cell polarity dynamics and motility in pattern formation due to contact-dependent signalling

A key challenge in biology is to understand how spatio-temporal patterns and structures arise during the development of an organism. An initial aggregate of spatially uniform cells develops and forms the differentiated structures of a fully developed organism. On the one hand, contact-dependent cell-cell signalling is responsible for generating a large number of complex, self-organized, spatial patterns in the distribution of the signalling molecules. On the other hand, the motility of cells coupled with their polarity can independently lead to collective motion patterns that depend on mechanical parameters influencing tissue deformation, such as cellular elasticity, cell-cell adhesion and active forces generated by actin and myosin dynamics. Although modelling efforts have, thus far, treated cell motility and cell-cell signalling separately, experiments in recent years suggest that these processes could be tightly coupled. Hence, in this paper, we study how the dynamics of cell polarity and migration influence the spatiotemporal patterning of signalling molecules. Such signalling interactions can occur only between cells that are in physical contact, either directly at the junctions of adjacent cells or through cellular protrusional contacts. We present a vertex model which accounts for contact-dependent signalling between adjacent cells and between non-adjacent neighbours through long protrusional contacts that occur along the orientation of cell polarization. We observe a rich variety of spatiotemporal patterns of signalling molecules that is influenced by polarity dynamics of the cells, relative strengths of adjacent and non-adjacent signalling interactions, range of polarized interaction, signalling activation threshold, relative time scales of signalling and polarity orientation, and cell motility. Though our results are developed in the context of Delta-Notch signalling, they are sufficiently general and can be extended to other contact dependent morpho-mechanical dynamics.

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.

Quantifying Material Properties of Cell Monolayers by Analyzing Integer Topological Defects

In developing organisms, internal cellular processes generate mechanical stresses at the tissue scale. The resulting deformations depend on the material properties of the tissue, which can exhibit long-ranged orientational order and topological defects. It remains a challenge to determine these properties on the time scales relevant for developmental processes. Here, we build on the physics of liquid crystals to determine material parameters of cell monolayers. Specifically, we use a hydrodynamic description to characterize the stationary states of compressible active polar fluids around defects. We illustrate our approach by analyzing monolayers of C2C12 cells in small circular confinements, where they form a single topological defect with integer charge. We find that such monolayers exert compressive stresses at the defect centers, where localized cell differentiation and formation of three-dimensional shapes is observed.

Integer topological defects of cell monolayers: Mechanics and flows

Monolayers of anisotropic cells exhibit long-ranged orientational order and topological defects. During the development of organisms, orientational order often influences morphogenetic events. However, the linkage between the mechanics of cell monolayers and topological defects remains largely unexplored. This holds specifically at the timescales relevant for tissue morphogenesis. Here, we build on the physics of liquid crystals to determine material parameters of cell monolayers. In particular, we use a hydrodynamical description of an active polar fluid to study the steady-state mechanical patterns at integer topological defects. Our description includes three distinct sources of activity: traction forces accounting for cell-substrate interactions as well as anisotropic and isotropic active nematic stresses accounting for cell-cell interactions. We apply our approach to C2C12 cell monolayers in small circular confinements, which form isolated aster or spiral topological defects. By analyzing the velocity and orientational order fields in spirals as well as the forces and cell number density fields in asters, we determine mechanical parameters of C2C12 cell monolayers. Our work shows how topological defects can be used to fully characterize the mechanical properties of biological active matter.

Dynamic Self-Organization of Idealized Migrating Cells by Contact Communication

This Letter investigates what forms of cellular dynamic self-organization are caused through intercellular contact communication based on a theoretical model in which migrating cells perform contact following and contact inhibition and attraction of locomotion. Tuning those strengths causes varieties of dynamic patterns. This further includes a novel form of collective migration, snakelike dynamic assembly. Scrutinizing this pattern reveals that cells in this state can accurately respond to an external directional cue but have no spontaneous global polar order.

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.

In primary airway epithelial cells, the unjamming transition is distinct from the epithelial-to-mesenchymal transition

The epithelial-to-mesenchymal transition (EMT) and the unjamming transition (UJT) each comprises a gateway to cellular migration, plasticity and remodeling, but the extent to which these core programs are distinct, overlapping, or identical has remained undefined. Here, we triggered partial EMT (pEMT) or UJT in differentiated primary human bronchial epithelial cells. After triggering UJT, cell-cell junctions, apico-basal polarity, and barrier function remain intact, cells elongate and align into cooperative migratory packs, and mesenchymal markers of EMT remain unapparent. After triggering pEMT these and other metrics of UJT versus pEMT diverge. A computational model attributes effects of pEMT mainly to diminished junctional tension but attributes those of UJT mainly to augmented cellular propulsion. Through the actions of UJT and pEMT working independently, sequentially, or interactively, those tissues that are subject to development, injury, or disease become endowed with rich mechanisms for cellular migration, plasticity, self-repair, and regeneration.

Universal Thinning of Liquid Filaments under Dominant Surface Forces

Theory and numerical simulations of the thinning of liquid threads at high superficial concentration of surfactants suggest the existence of an asymptotic regime where surface tension balances surface viscous stresses, leading to an exponential thinning with an e-fold time F(Θ)(3μ_{s}+κ_{s})/σ, where μ_{s} and κ_{s} are the surface shear and dilatational viscosity coefficients, σ is the interfacial tension, Θ=κ_{s}/μ_{s}, and F(Θ) is a universal function. The potential use of this phenomenon to measure the surface viscosity coefficients is discussed.

Pattern Formation and Defect Ordering in Active Chiral Nematics

Many biological systems display intriguing chiral patterns and dynamics. Here, we present an active nematic theory accounting for individual spin to explore the collective handedness in chiral rod-shaped aggregations. We show that coordinated individual spin and motility can engender a vortex-array pattern with chirality and drive ordering of topological defects. During this chiral process, the stationary trefoil-like defects self-organize into a periodic, hexagon-dominated polygonal network, which segregates persistently rotating cometlike defects in pairs within each polygon, leading to a translation symmetry at the global scale while a broken reflection symmetry at the local scale. Such defect ordering agrees exactly with the Voronoi tiling of two-dimensional space and the emergence of the hexagonal symmetry is deciphered in analogy with topological charge neutralization. We calculate energy barriers to the topological transition of the defect ordering and explain the existing metastable states with nonhexagonal polygons. Our findings shed light on the chiral morphodynamics in life processes and also suggest a potential route towards tuning self-organization in active materials.

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.