Epithelial–mesenchymal transition (EMT) and its reverse mesenchymal–epithelial transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGF β need not be irreversible, but have relatively longer time scales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model to account for the chromatin-mediated epigenetic silencing of the miR-200 family by ZEB family, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGF β withdrawal. Our results suggest that epigenetic modifiers can govern the extent and time scale of EMT reversibility and advise caution against labelling phenotypic changes seen in long-term EMT induction as ‘irreversible’.

Phenotypic heterogeneity is a hallmark of aggressive cancer behaviour and a clinical challenge. Despite much characterisation of this heterogeneity at a multi-omics level in many cancers, we have a limited understanding of how this heterogeneity emerges spontaneously in an isogenic cell population. Some longitudinal observations of dynamics in epithelial-mesenchymal heterogeneity, a canonical example of phenotypic heterogeneity, have offered us opportunities to quantify the rates of phenotypic switching that may drive such heterogeneity. Here, we offer a mathematical modeling framework that explains the salient features of population dynamics noted in PMC42-LA cells: (a) predominance of EpCAMhigh subpopulation, (b) re-establishment of parental distributions from the EpCAMhigh and EpCAMlow subpopulations, and (c) enhanced heterogeneity in clonal populations established from individual cells. Our framework proposes that fluctuations or noise in content duplication and partitioning of SNAIL—an EMT-inducing transcription factor—during cell division can explain spontaneous phenotypic switching and consequent dynamic heterogeneity in PMC42-LA cells observed experimentally at both single-cell and bulk level analysis. Together, we propose that asymmetric cell division can be a potential mechanism for phenotypic heterogeneity.

Intratumoral heterogeneity can exist along multiple axes: Cancer stem cells (CSCs)/non-CSCs, drug-sensitive/drug-tolerant states, and a spectrum of epithelial–hybrid–mesenchymal phenotypes. Further, these diverse cell-states can switch reversibly among one another, thereby posing a major challenge to therapeutic efficacy. Therefore, understanding the origins of phenotypic plasticity and heterogeneity remains an active area of investigation. While genomic components (mutations, chromosomal instability) driving heterogeneity have been well-studied, recent reports highlight the role of non-genetic mechanisms in enabling both phenotypic plasticity and heterogeneity. Here, we discuss various processes underlying phenotypic plasticity such as stochastic gene expression, chromatin reprogramming, asymmetric cell division and the presence of multiple stable gene expression patterns (‘attractors’). These processes can facilitate a dynamically evolving cell population such that a subpopulation of (drug-tolerant) cells can survive lethal drug exposure and recapitulate population heterogeneity on drug withdrawal, leading to relapse. These drug-tolerant cells can be both pre-existing and also induced by the drug itself through cell-state reprogramming. The dynamics of cell-state transitions both in absence and presence of the drug can be quantified through mathematical models. Such a dynamical systems approach to elucidating patterns of intratumoral heterogeneity by integrating longitudinal experimental data with mathematical models can help design effective combinatorial and/or sequential therapies for better clinical outcomes.

Epithelial-Mesenchymal Transition (EMT) and its reverse Mesenchymal-Epithelial Transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGFβ need not be irreversible, but have relatively longer timescales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model incorporating the epigenetic silencing of miR-200 by ZEB, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGFβ withdrawal. Our results suggest that epigenetic modifiers can govern the extent and timescale of EMT reversibility, and advise caution against labelling phenotypic changes seen in long-term EMT induction as 'irreversible'.

Phenotypic heterogeneity is a hallmark of aggressive cancer behaviour and a clinical challenge. Despite much characterisation of this heterogeneity at a multi-omics level in many cancers, we have a limited understanding of how this heterogeneity emerges spontaneously in an isogenic cell population. Some longitudinal observations of dynamics in epithelial-mesenchymal heterogeneity, a canonical example of phenotypic heterogeneity, have offered us opportunities to quantify the rates of phenotypic switching that may drive such heterogeneity. Here, we offer a mathematical modeling framework that explains the salient features of population dynamics noted in PMC42-LA cells: a) predominance of EpCAMhigh subpopulation, b) re-establishment of parental distributions from the EpCAMhigh and EpCAMlow subpopulations, and c) enhanced heterogeneity in clonal populations established from individual cells. Our framework proposes that fluctuations or noise in content duplication and partitioning of SNAIL - an EMT-inducing transcription factor - during cell division can explain spontaneous phenotypic switching and consequent dynamic heterogeneity in PMC42-LA cells observed experimentally at both single-cell and bulk level analysis. Together, we propose that asymmetric cell division can be a potential mechanism for phenotypic heterogeneity.