Epithelial-to-mesenchymal transition proceeds through directional destabilization of multidimensional attractor

  1. Weikang Wang  Is a corresponding author
  2. Dante Poe
  3. Yaxuan Yang
  4. Thomas Hyatt
  5. Jianhua Xing  Is a corresponding author
  1. University of Pittsburgh, United States

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

The following data sets were generated

Article and author information

Author details

  1. Weikang Wang

    Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, United States
    For correspondence
    weikang@pitt.edu
    Competing interests
    The authors declare that no competing interests exist.
  2. Dante Poe

    Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Yaxuan Yang

    Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Thomas Hyatt

    Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Jianhua Xing

    Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, United States
    For correspondence
    xing1@pitt.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3700-8765

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,468
    views
  • 411
    downloads
  • 21
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Weikang Wang
  2. Dante Poe
  3. Yaxuan Yang
  4. Thomas Hyatt
  5. Jianhua Xing
(2022)
Epithelial-to-mesenchymal transition proceeds through directional destabilization of multidimensional attractor
eLife 11:e74866.
https://doi.org/10.7554/eLife.74866

Share this article

https://doi.org/10.7554/eLife.74866

Further reading

    1. Cell Biology
    2. Physics of Living Systems
    Marta Urbanska, Yan Ge ... Jochen Guck
    Research Article

    Cell mechanical properties determine many physiological functions, such as cell fate specification, migration, or circulation through vasculature. Identifying factors that govern the mechanical properties is therefore a subject of great interest. Here, we present a mechanomics approach for establishing links between single-cell mechanical phenotype changes and the genes involved in driving them. We combine mechanical characterization of cells across a variety of mouse and human systems with machine learning-based discriminative network analysis of associated transcriptomic profiles to infer a conserved network module of five genes with putative roles in cell mechanics regulation. We validate in silico that the identified gene markers are universal, trustworthy, and specific to the mechanical phenotype across the studied mouse and human systems, and demonstrate experimentally that a selected target, CAV1, changes the mechanical phenotype of cells accordingly when silenced or overexpressed. Our data-driven approach paves the way toward engineering cell mechanical properties on demand to explore their impact on physiological and pathological cell functions.

    1. Physics of Living Systems
    M Julia Maristany, Anne Aguirre Gonzalez ... Jerelle A Joseph
    Research Article

    Proteins containing prion-like low complexity domains (PLDs) are common drivers of the formation of biomolecular condensates and are prone to misregulation due to amino acid mutations. Here, we exploit the accuracy of our residue-resolution coarse-grained model, Mpipi, to quantify the impact of amino acid mutations on the stability of 140 PLD mutants from six proteins (hnRNPA1, TDP43, FUS, EWSR1, RBM14, and TIA1). Our simulations reveal the existence of scaling laws that quantify the range of change in the critical solution temperature of PLDs as a function of the number and type of amino acid sequence mutations. These rules are consistent with the physicochemical properties of the mutations and extend across the entire family tested, suggesting that scaling laws can be used as tools to predict changes in the stability of PLD condensates. Our work offers a quantitative lens into how the emergent behavior of PLD solutions vary in response to physicochemical changes of single PLD molecules.