The epithelial–mesenchymal transition (EMT) and cancer stem cell (CSC) formation are two paramount processes driving tumor progression, therapy resistance, and cancer metastasis. Recent experiments show that cells with varying EMT and CSC phenotypes are spatially segregated in the primary tumor. The underlying mechanisms generating such spatiotemporal dynamics in the tumor microenvironment, however, remain largely unexplored. Here, we show through a mechanism-based dynamical model that the diffusion of EMT-inducing signals such as TGF-β, together with noncell autonomous control of EMT and CSC decision making via the Notch signaling pathway, can explain experimentally observed disparate localization of subsets of CSCs with varying EMT phenotypes in the tumor. Our simulations show that the more mesenchymal CSCs lie at the invasive edge, while the hybrid epithelial/mesenchymal (E/M) CSCs reside in the tumor interior. Further, motivated by the role of Notch-Jagged signaling in mediating EMT and stemness, we investigated the microenvironmental factors that promote Notch-Jagged signaling. We show that many inflammatory cytokines such as IL-6 that can promote Notch-Jagged signaling can (i) stabilize a hybrid E/M phenotype, (ii) increase the likelihood of spatial proximity of hybrid E/M cells, and (iii) expand the fraction of CSCs. To validate the predicted connection between Notch-Jagged signaling and stemness, we knocked down JAG1 in hybrid E/M SUM149 human breast cancer cells in vitro. JAG1 knockdown significantly restricted tumor organoid formation, confirming the key role that Notch-Jagged signaling can play in tumor progression. Together, our integrated computational–experimental framework reveals the underlying principles of spatiotemporal dynamics of EMT and CSCs.

Epithelial–mesenchymal transition (EMT) plays key roles during embryonic development, wound healing and cancer metastasis. Cells in a partial EMT or hybrid epithelial/mesenchymal (E/M) phenotype exhibit collective cell migration, forming clusters of circulating tumour cells—the primary drivers of metastasis. Activation of cell–cell signalling pathways such as Notch fosters a partial or complete EMT, yet the mechanisms enabling cluster formation remain poorly understood. Using an integrated computational–experimental approach, we examine the role of Numb—an inhibitor of Notch intercellular signalling—in mediating EMT and clusters formation. We show via an mathematical model that Numb inhibits a full EMT by stabilizing a hybrid E/M phenotype. Consistent with this observation, knockdown of Numb in stable hybrid E/M cells H1975 results in a full EMT, thereby showing that Numb acts as a brake for a full EMT and thus behaves as a ‘phenotypic stability factor' by modulating Notch-driven EMT. By generalizing the mathematical model to a multi-cell level, Numb is predicted to alter the balance of hybrid E/M versus mesenchymal cells in clusters, potentially resulting in a higher tumour-initiation ability. Finally, Numb correlates with a worse survival in multiple independent lung and ovarian cancer datasets, hence confirming its relationship with increased cancer aggressiveness.

The epithelial-mesenchymal transition (EMT) is a key process implicated in cancer metastasis and therapy resistance. Recent studies have emphasized that cells can undergo partial EMT to attain a hybrid epithelial/mesenchymal (E/M) phenotype – a cornerstone of tumour aggressiveness and poor prognosis. These cells can have enhanced tumour-initiation potential as compared to purely epithelial or mesenchymal ones and can integrate the properties of cell-cell adhesion and motility that facilitates collective cell migration leading to clusters of circulating tumour cells (CTCs) – the prevalent mode of metastasis. Thus, identifying the molecular players that can enable cells to maintain a hybrid E/M phenotype is crucial to curb the metastatic load. Using an integrated computational-experimental approach, we show that the transcription factor NRF2 can prevent a complete EMT and instead stabilize a hybrid E/M phenotype. Knockdown of NRF2 in hybrid E/M non-small cell lung cancer cells H1975 and bladder cancer cells RT4 destabilized a hybrid E/M phenotype and compromised the ability to collectively migrate to close a wound in vitro. Notably, while NRF2 knockout simultaneously downregulated E-cadherin and ZEB-1, overexpression of NRF2 enriched for a hybrid E/M phenotype by simultaneously upregulating both E-cadherin and ZEB-1 in individual RT4 cells. Further, we predict that NRF2 is maximally expressed in hybrid E/M phenotype(s) and demonstrate that this biphasic dynamic arises from the interconnections among NRF2 and the EMT regulatory circuit. Finally, clinical records from multiple datasets suggest a correlation between a hybrid E/M phenotype, high levels of NRF2 and its targets and poor survival, further strengthening the emerging notion that hybrid E/M phenotype(s) may occupy the ‘metastatic sweet spot’.