Hysteresis control of epithelial-mesenchymal transition dynamics conveys a distinct program with enhanced metastatic ability

Epithelial-mesenchymal transition (EMT) have been extensively characterized in development and cancer, and its dynamics have been modeled as a non-linear process. However, less is known about how such dynamics may affect its biological impact. Here, we use mathematical modeling and experimental analysis of the TGF-β-induced EMT to reveal a non-linear hysteretic response of E-cadherin repression tightly controlled by the strength of the miR-200s/ZEBs negative feedback loop. Hysteretic EMT conveys memory state, ensures rapid and robust cellular response and enables EMT to persist long after withdrawal of stimuli. Importantly, while both hysteretic and non-hysteretic EMT confer similar morphological changes and invasive potential of cancer cells, only hysteretic EMT enhances lung metastatic colonization efficiency. Cells that undergo hysteretic EMT differentially express subsets of stem cell and extracellular matrix related genes with significant clinical prognosis value. These findings illustrate distinct biological impact of EMT depending on the dynamics of the transition.The epithelial-mesenchymal transition (EMT) is a dynamic process that plays important roles in cancer progression and metastasis. Here, the authors characterize a non-linear hysteretic response of E-cadherin repression during TGFβ-induced EMT that is controlled by the strength of the miR-200s/ZEBs negative feedback loop and enhances metastasis.

[1]  Cathryn M. Gould,et al.  Whole-genome multiparametric screening to identify modulators of epithelial-to-mesenchymal transition. , 2014, Assay and drug development technologies.

[2]  Eshel Ben-Jacob,et al.  Coupling the modules of EMT and stemness: A tunable ‘stemness window’ model , 2015, Oncotarget.

[3]  ειδικούς στόχους,et al.  (2016) , 2018 .

[4]  W. Schiesser.  The Numerical Method of Lines: Integration of Partial Differential Equations , 1991 .

[5]  A. Saltelli,et al.  Making best use of model evaluations to compute sensitivity indices , 2002 .

[6]  Jianhua Xing,et al.  Coupled reversible and irreversible bistable switches underlying TGFβ-induced epithelial to mesenchymal transition. , 2013, Biophysical journal.

[7]  I. Sobola,et al.  Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .

[8]  Herschel Rabitz,et al.  Efficient Implementation of High Dimensional Model Representations , 2001 .

[9]  M. Nieto.  Epithelial Plasticity: A Common Theme in Embryonic and Cancer Cells , 2013, Science.

[10]  V. LeBleu,et al.  EMT Program is Dispensable for Metastasis but Induces Chemoresistance in Pancreatic Cancer , 2015, Nature.

[11]  Kakajan Komurov,et al.  Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes , 2010, Proceedings of the National Academy of Sciences.

[12]  A. García de Herreros,et al.  Cooperation, amplification, and feed-back in epithelial-mesenchymal transition. , 2012, Biochimica et biophysica acta.

[13]  Robert A. Weinberg,et al.  Epithelial-mesenchymal transition: at the crossroads of development and tumor metastasis. , 2008, Developmental cell.

[14]  C. Prieur,et al.  Generalized Hoeffding-Sobol Decomposition for Dependent Variables -Application to Sensitivity Analysis , 2011, 1112.1788.

[15]  G. Goodall,et al.  The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1 , 2008, Nature Cell Biology.

[16]  M. F. Shannon,et al.  A double-negative feedback loop between ZEB1-SIP1 and the microRNA-200 family regulates epithelial-mesenchymal transition. , 2008, Cancer research.

[17]  Jing Yang,et al.  Epithelial–mesenchymal plasticity in carcinoma metastasis , 2013, Genes & development.

[18]  J. Hájek,et al.  Asymptotic Normality of Simple Linear Rank Statistics Under Alternatives II , 1968 .

[19]  C. Heldin,et al.  TGF-beta and the Smad signaling pathway support transcriptomic reprogramming during epithelial-mesenchymal cell transition. , 2005, Molecular biology of the cell.

[20]  Sun-Mi Park,et al.  The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2. , 2008, Genes & development.

[21]  Jingyu Zhang,et al.  TGF-β–induced epithelial-to-mesenchymal transition proceeds through stepwise activation of multiple feedback loops , 2014, Science Signaling.

[22]  R. Gomis,et al.  Epithelial-mesenchymal transition can suppress major attributes of human epithelial tumor-initiating cells. , 2012, The Journal of clinical investigation.

[23]  S. A. Robertson,et al.  NONLINEAR OSCILLATIONS, DYNAMICAL SYSTEMS, AND BIFURCATIONS OF VECTOR FIELDS (Applied Mathematical Sciences, 42) , 1984 .

[24]  p53 regulates epithelial-mesenchymal transition (EMT) and stem cell properties through modulating miRNAs , 2010, Nature Cell Biology.

[25]  Simone Brabletz,et al.  The ZEB/miR‐200 feedback loop—a motor of cellular plasticity in development and cancer? , 2010, EMBO reports.

[26]  P. Howe,et al.  Pleiotropic Roles of Non-Coding RNAs in TGF-β-Mediated Epithelial-Mesenchymal Transition and Their Functions in Tumor Progression , 2017, Cancers.

[27]  B. Cieply,et al.  Suppression of the epithelial-mesenchymal transition by Grainyhead-like-2. , 2012, Cancer research.

[28]  Katherine C. Chen,et al.  Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. , 2003, Current opinion in cell biology.

[29]  J. Lieberman,et al.  miR-200 promotes the mesenchymal to epithelial transition by suppressing multiple members of the Zeb2 and Snail1 transcriptional repressor complexes , 2016, Oncogene.

[30]  Eduard Batlle,et al.  Metastatic stem cells: sources, niches, and vital pathways. , 2014, Cell stem cell.

[31]  G. Hooker.  Generalized Functional ANOVA Diagnostics for High-Dimensional Functions of Dependent Variables , 2007 .

[32]  Yutaka Kawakami,et al.  Cancer metastasis is accelerated through immunosuppression during Snail-induced EMT of cancer cells. , 2009, Cancer cell.

[33]  W. Hoeffding.  A Class of Statistics with Asymptotically Normal Distribution , 1948 .

[34]  R. Weinberg,et al.  Upholding a role for EMT in breast cancer metastasis , 2017, Nature.

[35]  R. Weinberg,et al.  Upholding a role for EMT in pancreatic cancer metastasis , 2017, Nature.

[36]  M. Korpal,et al.  The miR-200 Family Inhibits Epithelial-Mesenchymal Transition and Cancer Cell Migration by Direct Targeting of E-cadherin Transcriptional Repressors ZEB1 and ZEB2* , 2008, Journal of Biological Chemistry.

[37]  M. Nieto,et al.  Metastatic colonization requires the repression of the epithelial-mesenchymal transition inducer Prrx1. , 2012, Cancer cell.

[38]  Eshel Ben-Jacob,et al.  Implications of the Hybrid Epithelial/Mesenchymal Phenotype in Metastasis , 2015, Front. Oncol..

[39]  J. Ferrell.  Self-perpetuating states in signal transduction: positive feedback, double-negative feedback and bistability. , 2002, Current opinion in cell biology.

[40]  J. Settleman,et al.  EMT, cancer stem cells and drug resistance: an emerging axis of evil in the war on cancer , 2010, Oncogene.

[41]  Jarno Tuimala,et al.  Gene expression in TGFbeta-induced epithelial cell differentiation in a three-dimensional intestinal epithelial cell differentiation model , 2006, BMC Genomics.

[42]  Kohei Miyazono,et al.  Differential Regulation of Epithelial and Mesenchymal Markers by δEF1 Proteins in Epithelial–Mesenchymal Transition Induced by TGF-β , 2007 .

[43]  Herschel Rabitz,et al.  General formulation of HDMR component functions with independent and correlated variables , 2011, Journal of Mathematical Chemistry.

[44]  Wenjun Guo,et al.  The Epithelial-Mesenchymal Transition Generates Cells with Properties of Stem Cells , 2008, Cell.

[45]  Eshel Ben-Jacob,et al.  MicroRNA-based regulation of epithelial–hybrid–mesenchymal fate determination , 2013, Proceedings of the National Academy of Sciences.

[46]  Double asymptotics for the chi-square statistic. , 2016, Statistics & probability letters.

[47]  H. Rabitz,et al.  General foundations of high‐dimensional model representations , 1999 .

[48]  Kohei Miyazono,et al.  Differential regulation of epithelial and mesenchymal markers by deltaEF1 proteins in epithelial mesenchymal transition induced by TGF-beta. , 2007, Molecular biology of the cell.

[49]  Jean Paul Thiery,et al.  EMT: 2016 , 2016, Cell.

[50]  H. Rabitz,et al.  Random sampling-high dimensional model representation (RS-HDMR) and orthogonality of its different order component functions. , 2006, The journal of physical chemistry. A.

[51]  Paola Annoni,et al.  Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index , 2010, Comput. Phys. Commun..

[52]  Stephen T. C. Wong,et al.  EMT is not required for lung metastasis but contributes to chemoresistance , 2016 .

[53]  Arghavan Louhghalam,et al.  Variance decomposition and global sensitivity for structural systems , 2010 .

[54]  G. Dontu,et al.  In vitro propagation and transcriptional profiling of human mammary stem/progenitor cells. , 2003, Genes & development.

[55]  M. F. Shannon,et al.  An autocrine TGF-β/ZEB/miR-200 signaling network regulates establishment and maintenance of epithelial-mesenchymal transition , 2011, Molecular biology of the cell.

[56]  Héctor Peinado,et al.  Snail, Zeb and bHLH factors in tumour progression: an alliance against the epithelial phenotype? , 2007, Nature Reviews Cancer.

[57]  T. Brabletz,et al.  A reciprocal repression between ZEB1 and members of the miR-200 family promotes EMT and invasion in cancer cells , 2008, EMBO reports.

[58]  A. Takemura.  Tensor Analysis of ANOVA Decomposition , 1983 .

[59]  J. Massagué,et al.  Transcriptional control by the TGF‐β/Smad signaling system , 2000 .

[60]  G. Pan,et al.  Sequential introduction of reprogramming factors reveals a time-sensitive requirement for individual factors and a sequential EMT–MET mechanism for optimal reprogramming , 2013, Nature Cell Biology.

[61]  C. Pilarsky,et al.  The EMT-activator Zeb1 is a key factor for cell plasticity and promotes metastasis in pancreatic cancer , 2017, Nature Cell Biology.

[62]  Hong Peng,et al.  Interactions between cancer stem cells and their niche govern metastatic colonization , 2011, Nature.

[63]  Elgene Lim,et al.  Open Access Research Article Transcriptome Analyses of Mouse and Human Mammary Cell Subpopulations Reveal Multiple Conserved Genes and Pathways , 2022 .

[64]  Yibin Kang,et al.  Distinctive properties of metastasis-initiating cells , 2016, Genes & development.

[65]  Z. Szallasi,et al.  An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients , 2010, Breast Cancer Research and Treatment.

[66]  R. Fisher.  014: On the "Probable Error" of a Coefficient of Correlation Deduced from a Small Sample. , 1921 .

[67]  C. Iacobuzio-Donahue,et al.  TGF-β Tumor Suppression through a Lethal EMT , 2016, Cell.

[68]  E. Çinlar.  Probability and Stochastics , 2011 .

[69]  Alicia Zhou,et al.  Mesenchyme Forkhead 1 (FOXC2) plays a key role in metastasis and is associated with aggressive basal-like breast cancers , 2007, Proceedings of the National Academy of Sciences.

[70]  A. van Oudenaarden,et al.  Plasticity between Epithelial and Mesenchymal States Unlinks EMT from Metastasis-Enhancing Stem Cell Capacity , 2016, Cell reports.

[71]  R. Weinberg,et al.  Distinct EMT programs control normal mammary stem cells and tumour-initiating cells , 2015, Nature.