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,441
- views
-
- 407
- downloads
-
- 21
- 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
-
- Cell Biology
- Physics of Living Systems
The regulation of cell physiology depends largely upon interactions of functionally distinct proteins and cellular components. These interactions may be transient or long-lived, but often affect protein motion. Measurement of protein dynamics within a cellular environment, particularly while perturbing protein function with small molecules, may enable dissection of key interactions and facilitate drug discovery; however, current approaches are limited by throughput with respect to data acquisition and analysis. As a result, studies using super-resolution imaging are typically drawing conclusions from tens of cells and a few experimental conditions tested. We addressed these limitations by developing a high-throughput single-molecule tracking (htSMT) platform for pharmacologic dissection of protein dynamics in living cells at an unprecedented scale (capable of imaging >106 cells/day and screening >104 compounds). We applied htSMT to measure the cellular dynamics of fluorescently tagged estrogen receptor (ER) and screened a diverse library to identify small molecules that perturbed ER function in real time. With this one experimental modality, we determined the potency, pathway selectivity, target engagement, and mechanism of action for identified hits. Kinetic htSMT experiments were capable of distinguishing between on-target and on-pathway modulators of ER signaling. Integrated pathway analysis recapitulated the network of known ER interaction partners and suggested potentially novel, kinase-mediated regulatory mechanisms. The sensitivity of htSMT revealed a new correlation between ER dynamics and the ability of ER antagonists to suppress cancer cell growth. Therefore, measuring protein motion at scale is a powerful method to investigate dynamic interactions among proteins and may facilitate the identification and characterization of novel therapeutics.
-
- Neuroscience
- Physics of Living Systems
Recent experimental studies showed that electrically coupled neural networks like in mammalian inferior olive nucleus generate synchronized rhythmic activity by the subthreshold sinusoidal-like oscillations of the membrane voltage. Understanding the basic mechanism and its implication of such phenomena in the nervous system bears fundamental importance and requires preemptively the connectome information of a given nervous system. Inspired by these necessities of developing a theoretical and computational model to this end and, however, in the absence of connectome information for the inferior olive nucleus, here we investigated interference phenomena of the subthreshold oscillations in the reference system Caenorhabditis elegans for which the structural anatomical connectome was completely known recently. We evaluated how strongly the sinusoidal wave was transmitted between arbitrary two cells in the model network. The region of cell-pairs that are good at transmitting waves changed according to the wavenumber of the wave, for which we named a wavenumber-dependent transmission map. Also, we unraveled that (1) the transmission of all cell-pairs disappeared beyond a threshold wavenumber, (2) long distance and regular patterned transmission existed in the body-wall muscles part of the model network, and (3) major hub cell-pairs of the transmission were identified for many wavenumber conditions. A theoretical and computational model presented in this study provided fundamental insight for understanding how the multi-path constructive/destructive interference of the subthreshold oscillations propagating on electrically coupled neural networks could generate wavenumber-dependent synchronized rhythmic activity.