Population dynamics of epithelial-mesenchymal heterogeneity in cancer cells

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.