The complex ecosystem in non small cell lung cancer invasion

Many tumors are characterized by genetic instability, producing an assortment of genetic variants of tumor cells called subclones. These tumors and their surrounding environments form complex multi-cellular ecosystems, where subclones compete for resources and cooperate to perform multiple tasks, including cancer invasion. Our recent empirical studies revealed existence of such distinct phenotypes of cancer cells, leaders and followers, in lung cancer. These two cellular subclones exchange a complex array of extracellular signals demonstrating a symbiotic relationship at the cellular level. Here, we develop a computational model of the microenvironment of the lung cancer ecosystem to explore how the interactions between subclones can advance or inhibit invasion. We found that, due to the complexity of the ecosystem, invasion may have very different dynamics characterized by the different levels of aggressiveness. By altering the signaling environment, we could alter the ecological relationship between the cell types and the overall ecosystem development. Competition between leader and follower cell populations (defined by the limited amount of resources), positive feedback within the leader cell population (controlled by the focal adhesion kinase and fibronectin signaling), and impact of the follower cells to the leaders (represented by yet undetermined proliferation signal) all had major effects on the outcome of the collective dynamics. Specifically, our analysis revealed a class of tumors (defined by the strengths of fibronectin signaling and competition) that are particularly sensitive to manipulations of the signaling environment. These tumors can undergo irreversible changes to the tumor ecosystem that outlast the manipulations of feedbacks and have a profound impact on invasive potential. Our study predicts a complex division of labor between cancer cell subclones and suggests new treatment strategies targeting signaling within the tumor ecosystem.

[1]  M. Bazhenov,et al.  Delayed Onset of Symptoms Through Feedback Interference in Chronic Cancers. , 2016, Convergent science physical oncology.

[2]  R. Macarthur.  Species packing and competitive equilibrium for many species. , 1970, Theoretical population biology.

[3]  J. Konen,et al.  Targeting adhesion signaling in KRAS, LKB1 mutant lung adenocarcinoma. , 2017, JCI insight.

[4]  J. Stachowicz.  Mutualism, Facilitation, and the Structure of Ecological Communities , 2001 .

[5]  D. Hanahan,et al.  Hallmarks of Cancer: The Next Generation , 2011, Cell.

[6]  Ignacio A. Rodriguez-Brenes,et al.  Evolutionary dynamics of feedback escape and the development of stem-cell–driven cancers , 2011, Proceedings of the National Academy of Sciences.

[7]  D. Camidge,et al.  Safety, pharmacokinetic, and pharmacodynamic phase I dose-escalation trial of PF-00562271, an inhibitor of focal adhesion kinase, in advanced solid tumors. , 2012, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[8]  E. Ruoslahti,et al.  The Walter Herbert Lecture. Control of cell motility and tumour invasion by extracellular matrix interactions. , 1992, British Journal of Cancer.

[9]  Holger Gerhardt,et al.  VEGF and endothelial guidance in angiogenic sprouting , 2008, Organogenesis.

[10]  P. Friedl.  Prespecification and plasticity: shifting mechanisms of cell migration. , 2004, Current opinion in cell biology.

[11]  H. Frieboes,et al.  Nonlinear modelling of cancer: bridging the gap between cells and tumours , 2010, Nonlinearity.

[12]  Liz Y. Han,et al.  Adrenergic modulation of focal adhesion kinase protects human ovarian cancer cells from anoikis. , 2010, The Journal of clinical investigation.

[13]  Jeff Gore,et al.  Turning ecology and evolution against cancer , 2014, Nature Reviews Cancer.

[14]  Yi Jiang,et al.  A cell-based model exhibiting branching and anastomosis during tumor-induced angiogenesis. , 2007, Biophysical journal.

[15]  S. Smale.  On the differential equations of species in competition , 1976, Journal of mathematical biology.

[16]  N. McGranahan,et al.  The causes and consequences of genetic heterogeneity in cancer evolution , 2013, Nature.

[17]  Inhibition of focal adhesion kinase (FAK) activity prevents anchorage-independent ovarian carcinoma cell growth and tumor progression , 2013, Clinical & Experimental Metastasis.

[18]  J. Berlin,et al.  Bevacizumab plus irinotecan, fluorouracil, and leucovorin for metastatic colorectal cancer. , 2004, The New England journal of medicine.

[19]  G. Hardin.  The competitive exclusion principle. , 1960, Science.

[20]  G. Mcmahon,et al.  PND-1186 FAK inhibitor selectively promotes tumor cell apoptosis in three-dimensional environments , 2010, Cancer biology & therapy.

[21]  Steven J. M. Jones,et al.  Comprehensive molecular portraits of human breast tumors , 2012, Nature.

[22]  Gernot Schaller,et al.  Multicellular tumor spheroid in an off-lattice Voronoi-Delaunay cell model. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[23]  J. Glazier,et al.  3D Multi-Cell Simulation of Tumor Growth and Angiogenesis , 2009, PloS one.

[24]  David Basanta,et al.  Exploiting ecological principles to better understand cancer progression and treatment , 2013, Interface Focus.

[25]  Yi Jiang,et al.  A three-dimensional collagen-fiber network model of the extracellular matrix for the simulation of the mechanical behaviors and micro structures , 2017, Computer methods in biomechanics and biomedical engineering.

[26]  J. Salk.  Clonal evolution in cancer , 2010 .

[27]  Alissa M. Weaver,et al.  Tumor Morphology and Phenotypic Evolution Driven by Selective Pressure from the Microenvironment , 2006, Cell.

[28]  Gavin Thurston,et al.  Inhibition of vascular endothelial growth factor (VEGF) signaling in cancer causes loss of endothelial fenestrations, regression of tumor vessels, and appearance of basement membrane ghosts. , 2004, The American journal of pathology.

[29]  M. Huston.  A General Hypothesis of Species Diversity , 1979, The American Naturalist.

[30]  D. Hayes.  Bevacizumab treatment for solid tumors: boon or bust? , 2011, JAMA.

[31]  Qing Nie,et al.  CORRIGENDUM: Nonlinear Growth Kinetics of Breast Cancer Stem Cells: Implications for Cancer Stem Cell Targeted Therapy , 2013, Scientific Reports.

[32]  Matthew D. Wessel,et al.  Antitumor activity and pharmacology of a selective focal adhesion kinase inhibitor, PF-562,271. , 2008, Cancer research.

[33]  Steven J. M. Jones,et al.  Comprehensive molecular portraits of human breast tumours , 2013 .

[34]  Robert Axelrod,et al.  Ecological therapy for cancer: defining tumors using an ecosystem paradigm suggests new opportunities for novel cancer treatments. , 2008, Translational oncology.

[35]  Anil K Sood,et al.  Gene alterations identified by expression profiling in tumor-associated endothelial cells from invasive ovarian carcinoma. , 2007, Cancer research.

[36]  P. Friedl,et al.  Migration of coordinated cell clusters in mesenchymal and epithelial cancer explants in vitro. , 1995, Cancer research.

[37]  Yi Jiang,et al.  Topography of Extracellular Matrix Mediates Vascular Morphogenesis and Migration Speeds in Angiogenesis , 2009, PLoS Comput. Biol..

[38]  R. Seydel.  Practical Bifurcation and Stability Analysis , 1994 .

[39]  Janna Paulsson,et al.  Prognostic relevance of cancer-associated fibroblasts in human cancer. , 2014, Seminars in cancer biology.

[40]  Eugene S. Kim,et al.  Potent VEGF blockade causes regression of coopted vessels in a model of neuroblastoma , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[41]  H. Dvorak,et al.  Concordant release of glycolysis proteins into the plasma preceding a diagnosis of ER+ breast cancer. , 2012, Cancer research.

[42]  J. Saltz,et al.  Image-guided genomics of phenotypically heterogeneous populations reveals vascular signalling during symbiotic collective cancer invasion , 2017, Nature Communications.

[43]  Kevin J. Cheung,et al.  Mammary collective cell migration involves transient loss of epithelial features and individual cell migration within the epithelium , 2012, Journal of Cell Science.

[44]  M. Gerlinger,et al.  How Darwinian models inform therapeutic failure initiated by clonal heterogeneity in cancer medicine , 2010, British Journal of Cancer.

[45]  Jie Huang,et al.  Role of focal adhesion kinase in regulating YB-1-mediated paclitaxel resistance in ovarian cancer. , 2013, Journal of the National Cancer Institute.

[46]  Peter Chesson,et al.  MacArthur's consumer-resource model , 1990 .

[47]  Natalia L. Komarova,et al.  Dynamics of Cancer: Mathematical Foundations of Oncology , 2014 .

[48]  D. Schlaepfer,et al.  FAK in cancer: mechanistic findings and clinical applications , 2014, Nature Reviews Cancer.

[49]  Mark Alber,et al.  Continuous macroscopic limit of a discrete stochastic model for interaction of living cells. , 2007, Physical review letters.

[50]  Ma Dong,et al.  Bevacizumab plus Irinotecan,Fluorouracil,and Leucovorin for Metastatic Colorectal Cancer , 2006 .

[51]  Kathleen H. Keeler,et al.  The Ecology of Mutualism , 1982 .

[52]  J. Folkman.  Tumor angiogenesis: therapeutic implications. , 1971, The New England journal of medicine.

[53]  Benjamin J. Raphael,et al.  Integrated Genomic Analyses of Ovarian Carcinoma , 2011, Nature.

[54]  Q. Nie,et al.  Nonlinear Growth Kinetics of Breast Cancer Stem Cells: Implications for Cancer Stem Cell Targeted Therapy , 2013, Scientific Reports.

[55]  E. Dejana,et al.  Inhibition of endothelial FAK activity prevents tumor metastasis by enhancing barrier function , 2014, The Journal of cell biology.

[56]  Y. Kuznetsov,et al.  New features of the software MatCont for bifurcation analysis of dynamical systems , 2008 .