Tracking collective cell motion by topological data analysis (2009.14161v1)

Abstract: 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, spreading tests exhibit finger formation in the moving interfaces, swirls in the velocity field, and the polar order parameter and correlation and swirl lengths increase with time. Cells inside the tissue have smaller area than those at the interface, as observed in recent experiments. In antagonistic migration assays, 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 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 numerical simulations. We use time series of numerical simulation data to automatically group, track and classify advancing interfaces of cellular aggregates by means of bottleneck or Wasserstein distances of persistent homologies. These topological data analysis techniques 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.