Collective Cell Migration: Leadership, Invasion and Segregation

A number of biological processes, such as embryo development, cancer metastasis or wound healing, rely on cells moving in concert. The mechanisms leading to the emergence of coordinated motion remain however largely unexplored. Although biomolecular signalling is known to be involved in most occurrences of collective migration, the role of physical and mechanical interactions has only been recently investigated. In this paper, a versatile framework for cell motility is implemented in-silico in order to study the minimal requirements for the coordination of a group of epithelial cells. We find that cell motility and cell-cell mechanical interactions are sufficient to generate a broad array of behaviours commonly observed in vitro and in vivo. Cell streaming, sheet migration and susceptibility to leader cells are examples of behaviours spontaneously emerging from these simple assumptions, which might explain why collective effects are so ubiquitous in nature. This analysis provides also new insights into cancer metastasis and cell sorting, suggesting in particular that collective invasion might result from an emerging coordination in a system where single cells are mechanically unable to invade.

Emerging modes of collective cell migration induced by geometrical constraints.

The role of geometrical confinement on collective cell migration has been recognized but has not been elucidated yet. Here, we show that the geometrical properties of the environment regulate the formation of collective cell migration patterns through cell-cell interactions. Using microfabrication techniques to allow epithelial cell sheets to migrate into strips whose width was varied from one up to several cell diameters, we identified the modes of collective migration in response to geometrical constraints. We observed that a decrease in the width of the strips is accompanied by an overall increase in the speed of the migrating cell sheet. Moreover, large-scale vortices over tens of cell lengths appeared in the wide strips whereas a contraction-elongation type of motion is observed in the narrow strips. Velocity fields and traction force signatures within the cellular population revealed migration modes with alternative pulling and/or pushing mechanisms that depend on extrinsic constraints. Force transmission through intercellular contacts plays a key role in this process because the disruption of cell-cell junctions abolishes directed collective migration and passive cell-cell adhesions tend to move the cells uniformly together independent of the geometry. Altogether, these findings not only demonstrate the existence of patterns of collective cell migration depending on external constraints but also provide a mechanical explanation for how large-scale interactions through cell-cell junctions can feed back to regulate the organization of migrating tissues.

Guidance of collective cell migration by substrate geometry.

Collective behavior refers to the emergence of complex migration patterns over scales larger than those of the individual elements constituting a system. It plays a pivotal role in biological systems in regulating various processes such as gastrulation, morphogenesis and tissue organization. Here, by combining experimental approaches and numerical modeling, we explore the role of cell density ('crowding'), strength of intercellular adhesion ('cohesion') and boundary conditions imposed by extracellular matrix (ECM) proteins ('constraints') in regulating the emergence of collective behavior within epithelial cell sheets. Our results show that the geometrical confinement of cells into well-defined circles induces a persistent, coordinated and synchronized rotation of cells that depends on cell density. The speed of such rotating large-scale movements slows down as the density increases. Furthermore, such collective rotation behavior depends on the size of the micropatterned circles: we observe a rotating motion of the overall cell population in the same direction for sizes of up to 200 μm. The rotating cells move as a solid body, with a uniform angular velocity. Interestingly, this upper limit leads to length scales that are similar to the natural correlation length observed for unconfined epithelial cell sheets. This behavior is strongly altered in cells that present a downregulation of adherens junctions and in cancerous cell types. We anticipate that our system provides a simple and easy approach to investigate collective cell behavior in a well-controlled and systematic manner.

Nonautonomous contact guidance signaling during collective cell migration.

Directed migration of groups of cells is a critical aspect of tissue morphogenesis that ensures proper tissue organization and, consequently, function. Cells moving in groups, unlike single cells, must coordinate their migratory behavior to maintain tissue integrity. During directed migration, cells are guided by a combination of mechanical and chemical cues presented by neighboring cells and the surrounding extracellular matrix. One important class of signals that guide cell migration includes topographic cues. Although the contact guidance response of individual cells to topographic cues has been extensively characterized, little is known about the response of groups of cells to topographic cues, the impact of such cues on cell-cell coordination within groups, and the transmission of nonautonomous contact guidance information between neighboring cells. Here, we explore these phenomena by quantifying the migratory response of confluent monolayers of epithelial and fibroblast cells to contact guidance cues provided by grooved topography. We show that, in both sparse clusters and confluent sheets, individual cells are contact-guided by grooves and show more coordinated behavior on grooved versus flat substrates. Furthermore, we demonstrate both in vitro and in silico that the guidance signal provided by a groove can propagate between neighboring cells in a confluent monolayer, and that the distance over which signal propagation occurs is not significantly influenced by the strength of cell-cell junctions but is an emergent property, similar to cellular streaming, triggered by mechanical exclusion interactions within the collective system.

Role of Vinculin in Regulating Force-Sensitive Dynamics of Adherens Junctions During Collective Cell Migration

Dynamic processes such as morphogenesis and tissue patterning require the precise control of many cellular processes, especially cell migration. Historically, these processes are thought to be mediated by genetic and biochemical signaling pathways. However, recent advances have unraveled a previously unappreciated role of mechanical forces in regulating these homeostatic processes in of multicellular systems. In multicellular systems cells adhere to both deformable extracellular matrix (ECM) and other cells, which are sources of applied forces and means of mechanical support. Cells detect and respond to these mechanical signals through a poorly understood process called mechanotransduction, which can have profound effects on processes such as cell migration. These effects are largely mediated by the sub cellular structures that link cells to the ECM, called focal adhesions (FAs), or cells to other cells, termed adherens junctions (AJs).

Overall this thesis is comprised of my work on identifying a novel force dependent function of vinculin, a protein which resides in both FAs and AJs - in dynamic process of collective migration. Using a collective migration assay as a model for collective cell behavior and a fluorescence resonance energy transfer (FRET) based molecular tension sensor for vinculin I demonstrated a spatial gradient of tension across vinculin in the direction of migration. To define this novel force-dependent role of vinculin in collective migration I took advantage of previously established shRNA based vinculin knock down Marin-Darby Canine Kidney (MDCK) epithelial cells.

The first part of my thesis comprises of my work demonstrating the mechanosensitive role of vinculin at AJ’s in collectively migrating cells. Using vinculin knockdown cells and vinculin mutants, which specifically disrupt vinculin’s ability to bind actin (VinI997A) or disrupt its ability to localize to AJs without affecting its localization at FAs (VinY822F), I establish a role of force across vinculin in E-cadherin internalization and clipping. Furthermore by measuring E-cadherin dynamics using fluorescence recovery after bleaching (FRAP) analysis I show that vinculin inhibition affects the turnover of E-cadherin at AJs. Together these data reveal a novel mechanosensitive role of vinculin in E-cadherin internalization and turnover in a migrating cell layer, which is contrary to the previously identified role of vinculin in potentiating E-cadherin junctions in a static monolayer.

For the last part of my thesis I designed a novel tension sensor to probe tension across N-cadherin (NTS). N-cadherin plays a critical role in cardiomyocytes, vascular smooth muscle cells, neurons and neural crest cells. Similar to E-cadherin, N-cadherin is also believed to bear tension and play a role in mechanotransduction pathways. To identify the role of tension across N-cadherin I designed a novel FRET-based molecular tension sensor for N-cadherin. I tested the ability of NTS to sense molecular tension in vascular smooth muscle cells, cardiomyocytes and cancer cells. Finally in collaboration with the Horwitz lab we have been able to show a role of tension across N-cadherin in synaptogenesis of neurons.

REGULATION OF SHAPE DYNAMICS AND ACTIN POLYMERIZATION DURING COLLECTIVE CELL MIGRATION

This thesis aims to understand how cells coordinate their motion during collective migration. As previously shown, the motion of individually migrating cells is governed by wave-like cell shape dynamics. The mechanisms that regulate these dynamic behaviors in response to extracellular environment remain largely unclear. I applied shape dynamics analysis to Dictyostelium cells migrating in pairs and in multicellular streams and found that wave-like membrane protrusions are highly coupled between touching cells. I further characterized cell motion by using principle component analysis (PCA) to decompose complex cell shape changes into a serial shape change modes, from which I found that streaming cells exhibit localized anterior protrusion, termed front narrowing, to facilitate cell-cell coupling. I next explored cytoskeleton-based mechanisms of cell-cell coupling by measuring the dynamics of actin polymerization. Actin polymerization waves observed in individual cells were significantly suppressed in multicellular streams. Streaming cells exclusively produced F-actin at cell-cell contact regions, especially at cell fronts. I demonstrated that such restricted actin polymerization is associated with cell-cell coupling, as reducing actin polymerization with Latrunculin A leads to the assembly of F-actin at the side of streams, the decrease of front narrowing, and the decoupling of protrusion waves. My studies also suggest that collective migration is guided by cell-surface interactions. I examined the aggregation of Dictyostelim cells under distinct conditions and found that both chemical compositions of surfaces and surface-adhesion defects in cells result in altered collective migration patterns. I also investigated the shape dynamics of cells suspended on PEG-coated surfaces, which showed that coupling of protrusion waves disappears on touching suspended cells. These observations indicate that collective migration requires a balance between cell-cell and cell-surface adhesions. I hypothesized such a balance is reached via the regulation of cytoskeleton. Indeed, I found cells actively regulate cytoskeleton to retain optimal cell-surface adhesions on varying surfaces, and cells lacking the link between actin and surfaces (talin A) could not retain the optimal adhesions. On the other hand, suspended cells exhibited enhanced actin filament assembly on the periphery of cell groups instead of in cell-cell contact regions, which facilitates their aggregation in a clumping fashion.

Models of collective cell behaviour with crowding effects : comparing lattice-based and lattice-free approaches

Individual-based models describing the migration and proliferation of a population of cells frequently restrict the cells to a predefined lattice. An implicit assumption of this type of lattice based model is that a proliferative population will always eventually fill the lattice. Here we develop a new lattice-free individual-based model that incorporates cell-to-cell crowding effects. We also derive approximate mean-field descriptions for the lattice-free model in two special cases motivated by commonly used experimental setups. Lattice-free simulation results are compared to these mean-field descriptions and to a corresponding lattice-based model. Data from a proliferation experiment is used to estimate the parameters for the new model, including the cell proliferation rate, showing that the model fits the data well. An important aspect of the lattice-free model is that the confluent cell density is not predefined, as with lattice-based models, but an emergent model property. As a consequence of the more realistic, irregular configuration of cells in the lattice-free model, the population growth rate is much slower at high cell densities and the population cannot reach the same confluent density as an equivalent lattice-based model.

Sensitivity of edge detection methods for quantifying cell migration assays

Quantitative imaging methods to analyze cell migration assays are not standardized. Here we present a suite of two–dimensional barrier assays describing the collective spreading of an initially–confined population of 3T3 fibroblast cells. To quantify the motility rate we apply two different automatic image detection methods to locate the position of the leading edge of the spreading population after 24, 48 and 72 hours. These results are compared with a manual edge detection method where we systematically vary the detection threshold. Our results indicate that the observed spreading rates are very sensitive to the choice of image analysis tools and we show that a standard measure of cell migration can vary by as much as 25% for the same experimental images depending on the details of the image analysis tools. Our results imply that it is very difficult, if not impossible, to meaningfully compare previously published measures of cell migration since previous results have been obtained using different image analysis techniques and the details of these techniques are not always reported. Using a mathematical model, we provide a physical interpretation of our edge detection results. The physical interpretation is important since edge detection algorithms alone do not specify any physical measure, or physical definition, of the leading edge of the spreading population. Our modeling indicates that variations in the image threshold parameter correspond to a consistent variation in the local cell density. This means that varying the threshold parameter is equivalent to varying the location of the leading edge in the range of approximately 1–5% of the maximum cell density.

Mathematical models for collective cell spreading

Collective cell spreading is frequently observed in development, tissue repair and disease progression. Mathematical modelling used in conjunction with experimental investigation can provide key insights into the mechanisms driving the spread of cell populations. In this study, we investigated how experimental and modelling frameworks can be used to identify several key features underlying collective cell spreading. In particular, we were able to independently quantify the roles of cell motility and cell proliferation in a spreading cell population, and investigate how these roles are influenced by factors such as the initial cell density, type of cell population and the assay geometry.

La petite GTPase Rab11 et ses interacteurs orchestrent la migration cellulaire collective et la cytocinèse chez la Drosophile

Le trafic vésiculaire permet un échange coordonné de molécules entre les différents organites de la cellule et dépend largement des petites GTPases de la famille des Rabs dont le nombre varie entre 27 chez la Drosophile et 70 chez l’Homme. Un des prochains défis consiste donc à élucider les mécanismes cellulaires qui coordonnent l’activité de ces Rabs, laquelle garantit un transport vésiculaire ordonné au sein de la cellule. Les Rabs agissent comme des interrupteurs moléculaires grâce à leur capacité à cycler entre un état actif et inactif. L’activité des Rabs est contrôlée par des protéines régulatrices puis des effecteurs en aval coordonnent leurs différentes fonctions. La petite GTPase Rab11 est essentielle au développement de plusieurs organismes incluant la Drosophile, C. elegans et la souris puisqu’elle se retrouve au cœur de différentes voies de transport. D’ailleurs, le trafic de molécules dépendant de Rab11 est perturbé dans plusieurs pathologies. Malgré son rôle central dans le trafic vésiculaire, la régulation de Rab11 reste peu comprise in vivo. Cette thèse se penche sur les mécanismes moléculaires contrôlant les fonctions de Rab11 et de ses effecteurs lors de la migration cellulaire collective et lors de la cytocinèse. Nous avons identifié Evi5 comme un nouvel acteur clé de la migration cellulaire collective, et nous montrons qu’elle possède une activité Rab11-GAP essentielle pour maintenir les récepteurs de guidance actifs de façon polarisée au front de migration. Nous avons ensuite déterminé que Rab11 régule la communication cellulaire lors de la migration collective par l’entremise de son interaction avec la Moésine. Une question reste toutefois en suspens : sachant que Rab11 compte plus de 13 effecteurs, quels sont les mécanismes assurant la spécificité de l’interaction entre cette GTPase et un effecteur particulier? Une partie de la réponse provient peut-être de nos observations que les membres des Rab11-FIPs de classe I, une famille d’effecteurs de Rab11, interagissent avec les protéines d’échafaudage 14-3-3. Chez la Drosophile, Rip11 est le seul représentant des Rab11-FIPs de classe I et nous montrons que Rip11 aurait des fonctions inattendues durant la cytocinèse qui seraient coordonnées par 14-3-3. Nos recherches permettent de dresser un portrait plus authentique des mécanismes moléculaires régulant les différentes fonctions de Rab11 et de ses effecteurs in vivo.

Mechanistic investigations into interactions between IGF-I and IGFBPs and their impact on facilitating cell migration on vitronectin

Numerous studies have reported links between insulin-like growth factors (IGFs) and the extra-cellular matrix protein vitronectin (VN). We ourselves have reported that IGF-I binds to VN via IGF-binding proteins (IGFBPs) to stimulate HaCaT and MCF-7 cell migration. Here, we detail the functional evaluation of IGFBP-1, -2, -3, -4 and -6 in the presence and absence of IGF-I and VN. The data presented here, combined with our prior data on IGFBP-5, suggest that IGFBP-3, -4 and -5 are the most effective at stimulating cell migration in combination with IGF-I and VN. In addition, we demonstrate that different regions within IGFBP-3 and -4 are critical for complex formation. Furthermore, we examine whether multi-protein complexes of IGF-I and IGFBPs associated with fibronectin and collagen IV are also able to enhance functional biological responses.

Models of collective cell spreading with variable cell aspect ratio : a motivation for degenerate diffusion models

Continuum diffusion models are often used to represent the collective motion of cell populations. Most previous studies have simply used linear diffusion to represent collective cell spreading, while others found that degenerate nonlinear diffusion provides a better match to experimental cell density profiles. In the cell modeling literature there is no guidance available with regard to which approach is more appropriate for representing the spreading of cell populations. Furthermore, there is no knowledge of particular experimental measurements that can be made to distinguish between situations where these two models are appropriate. Here we provide a link between individual-based and continuum models using a multi-scale approach in which we analyze the collective motion of a population of interacting agents in a generalized lattice-based exclusion process. For round agents that occupy a single lattice site, we find that the relevant continuum description of the system is a linear diffusion equation, whereas for elongated rod-shaped agents that occupy L adjacent lattice sites we find that the relevant continuum description is connected to the porous media equation (pme). The exponent in the nonlinear diffusivity function is related to the aspect ratio of the agents. Our work provides a physical connection between modeling collective cell spreading and the use of either the linear diffusion equation or the pme to represent cell density profiles. Results suggest that when using continuum models to represent cell population spreading, we should take care to account for variations in the cell aspect ratio because different aspect ratios lead to different continuum models.

A cell migration device that maintains a defined surface with no cellular damage during wound edge generation

Studying the rate of cell migration provides insight into fundamental cell biology as well as a tool to assess the functionality of synthetic surfaces and soluble environments used in tissue engineering. The traditional tools used to study cell migration include the fence and wound healing assays. In this paper we describe the development of a microchannel based device for the study of cell migration on defined surfaces. We demonstrate that this device provides a superior tool, relative to the previously mentioned assays, for assessing the propagation rate of cell wave fronts. The significant advantage provided by this technology is the ability to maintain a virgin surface prior to the commencement of the cell migration assay. Here, the device is used to assess rates of mouse fibroblasts (NIH 3T3) and human osteosarcoma (SaOS2) cell migration on surfaces functionalized with various extracellular matrix proteins as a demonstration that confining cell migration within a microchannel produces consistent and robust data. The device design enables rapid and simplistic assessment of multiple repeats on a single chip, where surfaces have not been previously exposed to cells or cellular secretions.

Gonadotropin-induced ovarian cancer cell migration and proliferation require extracellular signal-regulated kinase 1/2 activation regulated by calcium and protein kinase Cδ

The gonadotropin hypothesis proposes that elevated serum gonadotropin levels may increase the risk of epithelial ovarian cancer (EOC). We have studied the effect of treating EOC cell lines (OV207 and OVCAR-3) with FSH or LH. Both gonadotropins activated the mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase 1/2 (ERK1/2) pathway and increased cell migration that was inhibited by the MAPK 1 inhibitor PD98059. Both extra- and intracellular calcium ion signalling were implicated in gonadotropin-induced ERK1/2 activation as treatment with either the calcium chelator EGTA or an inhibitor of intracellular calcium release, dantrolene, inhibited gonadotropin-induced ERK1/2 activation. Verapamil was also inhibitory, indicating that gonadotropins activate calcium influx via L-type voltage-dependent calcium channels. The cAMP/protein kinase A (PKA) pathway was not involved in the mediation of gonadotropin action in these cells as gonadotropins did not increase intracellular cAMP formation and inhibition of PKA did not affect gonadotropin-induced phosphorylation of ERK1/2. Activation of ERK1/2 was inhibited by the protein kinase C (PKC) inhibitor GF 109203X as well as by the PKCδ inhibitor rottlerin, and downregulation of PKCδ was inhibited by small interfering RNA (siRNA), highlighting the importance of PKCδ in the gonadotropin signalling cascade. Furthermore, in addition to inhibition by PD98059, gonadotropin-induced ovarian cancer cell migration was also inhibited by verapamil, GF 109203X and rottlerin. Similarly, gonadotropin-induced proliferation was inhibited by PD98059, verapamil, GF 109203X and PKCδ siRNA. Taken together, these results demonstrate that gonadotropins induce both ovarian cancer cell migration and proliferation by activation of ERK1/2 signalling in a calcium- and PKCδ-dependent manner.

Models of collective cell motion for cell populations with different aspect ratio : diffusion, proliferation and travelling waves

Continuum, partial differential equation models are often used to describe the collective motion of cell populations, with various types of motility represented by the choice of diffusion coefficient, and cell proliferation captured by the source terms. Previously, the choice of diffusion coefficient has been largely arbitrary, with the decision to choose a particular linear or nonlinear form generally based on calibration arguments rather than making any physical connection with the underlying individual-level properties of the cell motility mechanism. In this work we provide a new link between individual-level models, which account for important cell properties such as varying cell shape and volume exclusion, and population-level partial differential equation models. We work in an exclusion process framework, considering aligned, elongated cells that may occupy more than one lattice site, in order to represent populations of agents with different sizes. Three different idealizations of the individual-level mechanism are proposed, and these are connected to three different partial differential equations, each with a different diffusion coefficient; one linear, one nonlinear and degenerate and one nonlinear and nondegenerate. We test the ability of these three models to predict the population level response of a cell spreading problem for both proliferative and nonproliferative cases. We also explore the potential of our models to predict long time travelling wave invasion rates and extend our results to two dimensional spreading and invasion. Our results show that each model can accurately predict density data for nonproliferative systems, but that only one does so for proliferative systems. Hence great care must be taken to predict density data for with varying cell shape.