Tracking collective cell motion by topological data analysis

• Published in: PLoS computational biology 2020-12, Vol.16 (12), p.e1008407-e1008407
• Main Authors: Bonilla, Luis L, Carpio, Ana, Trenado, Carolina
• Format: Article
• Language: English
• Subjects: Aggregates; Biological activity; Biology and Life Sciences; Cancer; Cell adhesion & migration; Cell migration; Cell Movement; Cells; Colorectal carcinoma; Data Analysis; Earth Sciences; Electronic data processing; Epithelial cells; Experiments; Humans; Interfaces; Mathematical models; Medicine and Health Sciences; Metastases; Metastasis; Methods; Microscopy; Model testing; Models, Biological; Monolayers; Morphogenesis; Motility; Order parameters; Pattern formation; Physical Sciences; Simulation; Swarming; Tissues; Tumors; Velocity; Velocity distribution; Wound Healing
• ISSN: 1553-7358; 1553-734X; 1553-7358
• DOI: 10.1371/journal.pcbi.1008407

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.