Epithelial-to-mesenchymal transition proceeds through directional destabilization of multidimensional attractor
Abstract
How a cell changes from one stable phenotype to another one is a fundamental problem in developmental and cell biology. Mathematically a stable phenotype corresponds to a stable attractor in a generally multi-dimensional state space, which needs to be destabilized so the cell relaxes to a new attractor. Two basic mechanisms for destabilizing a stable fixed point, pitchfork and saddle-node bifurcations, have been extensively studied theoretically, however direct experimental investigation at the single cell level remains scarce. Here we performed live cell imaging studies and analyses in the framework of dynamical systems theories on epithelial-to-mesenchymal transition (EMT). While some mechanistic details remain controversial, EMT is a cell phenotypic transition (CPT) process central to development and pathology. Through time-lapse imaging we recorded single cell trajectories of human A549/Vim-RFP cells undergoing EMT induced by different concentrations of exogenous TGF-β in a multi-dimensional cell feature space. The trajectories clustered into two distinct groups, indicating that the transition dynamics proceeds through parallel paths. We then reconstructed the reaction coordinates and the corresponding quasi-potentials from the trajectories. The potentials revealed a plausible mechanism for the emergence of the two paths where the original stable epithelial attractor collides with two saddle points sequentially with increased TGF-β concentration, and relaxes to a new one. Functionally the directional saddle-node bifurcation ensures a CPT proceeds towards a specific cell type, as a mechanistic realization of the canalization idea proposed by Waddington.
Data availability
The computer code are shared on GitHub, so other researchers can run to reproduce Figure 3, 4, and 5. The processed single cell trajectory data are on Dryad
-
A549 VIM-RFP EMT single cell trajectory dataDryad Digital Repository, RRID:SCR_005910.
Article and author information
Author details
Funding
National Institute of Diabetes and Digestive and Kidney Diseases (R01DK119232)
- Jianhua Xing
National Cancer Institute (R37 CA232209)
- Jianhua Xing
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2022, Wang et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,527
- views
-
- 417
- downloads
-
- 23
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Physics of Living Systems
The tumor stroma consists mainly of extracellular matrix, fibroblasts, immune cells, and vasculature. Its structure and functions are altered during malignancy: tumor cells transform fibroblasts into cancer-associated fibroblasts, which exhibit immunosuppressive activities on which growth and metastasis depend. These include exclusion of immune cells from the tumor nest, cancer progression, and inhibition of T-cell-based immunotherapy. To understand these complex interactions, we measure the density of different cell types in the stroma using immunohistochemistry techniques on tumor samples from lung cancer patients. We incorporate these data into a minimal dynamical system, explore the variety of outcomes, and finally establish a spatio-temporal model that explains the cell distribution. We reproduce that cancer-associated fibroblasts act as a barrier to tumor expansion, but also reduce the efficiency of the immune response. Our conclusion is that the final outcome depends on the parameter values for each patient and leads to either tumor invasion, persistence, or eradication as a result of the interplay between cancer cell growth, T-cell cytotoxicity, and fibroblast activity. However, despite the existence of a wide range of scenarios, distinct trajectories, and patterns allow quantitative predictions that may help in the selection of new therapies and personalized protocols.
-
- Neuroscience
- Physics of Living Systems
Neurons generate and propagate electrical pulses called action potentials which annihilate on arrival at the axon terminal. We measure the extracellular electric field generated by propagating and annihilating action potentials and find that on annihilation, action potentials expel a local discharge. The discharge at the axon terminal generates an inhomogeneous electric field that immediately influences target neurons and thus provokes ephaptic coupling. Our measurements are quantitatively verified by a powerful analytical model which reveals excitation and inhibition in target neurons, depending on position and morphology of the source-target arrangement. Our model is in full agreement with experimental findings on ephaptic coupling at the well-studied Basket cell-Purkinje cell synapse. It is able to predict ephaptic coupling for any other synaptic geometry as illustrated by a few examples.