A spatial model predicts that dispersal and cell turnover limit intratumour heterogeneity

Most cancers in humans are large, measuring centimetres in diameter, and composed of many billions of cells. An equivalent mass of normal cells would be highly heterogeneous as a result of the mutations that occur during each cell division. What is remarkable about cancers is that virtually every neoplastic cell within a large tumour often contains the same core set of genetic alterations, with heterogeneity confined to mutations that emerge late during tumour growth. How such alterations expand within the spatially constrained three-dimensional architecture of a tumour, and come to dominate a large, pre-existing lesion, has been unclear. Here we describe a model for tumour evolution that shows how short-range dispersal and cell turnover can account for rapid cell mixing inside the tumour. We show that even a small selective advantage of a single cell within a large tumour allows the descendants of that cell to replace the precursor mass in a clinically relevant time frame. We also demonstrate that the same mechanisms can be responsible for the rapid onset of resistance to chemotherapy. Our model not only provides insights into spatial and temporal aspects of tumour growth, but also suggests that targeting short-range cellular migratory activity could have marked effects on tumour growth rates.

[1]  Michael Thomaschewski,et al.  RGB marking facilitates multicolor clonal cell tracking , 2011, Nature Medicine.

[2]  G. Booth,et al.  BacSim, a simulator for individual-based modelling of bacterial colony growth. , 1998, Microbiology.

[3]  Hans G Othmer,et al.  A Hybrid Model of Tumor–Stromal Interactions in Breast Cancer , 2013, Bulletin of Mathematical Biology.

[4]  K. Foster,et al.  A Quantitative Test of Population Genetics Using Spatiogenetic Patterns in Bacterial Colonies , 2011, The American Naturalist.

[5]  E. Albano,et al.  Interfacial properties in a discrete model for tumor growth , 2013 .

[6]  Frank Jülicher,et al.  Fluidization of tissues by cell division and apoptosis , 2010, Proceedings of the National Academy of Sciences.

[7]  M. Loeffler,et al.  Modeling the effect of deregulated proliferation and apoptosis on the growth dynamics of epithelial cell populations in vitro. , 2005, Biophysical journal.

[8]  Nico Karssemeijer,et al.  Ultrasonics, Ferroelectrics, and Frequency Control , 2011 .

[9]  Triantafyllos Stylianopoulos,et al.  Combining two strategies to improve perfusion and drug delivery in solid tumors , 2013, Proceedings of the National Academy of Sciences.

[10]  N. Komarova Spatial Stochastic Models for Cancer Initiation and Progression , 2006, Bulletin of mathematical biology.

[11]  W. O’Brien,et al.  Young's modulus measurements of soft tissues with application to elasticity imaging , 1996, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[12]  S. Baylin,et al.  Absence of E-Cadherin Expression Distinguishes Noncohesive from Cohesive Pancreatic Cancer , 2008, Clinical Cancer Research.

[13]  D. Drasdo,et al.  Buckling instabilities of one-layered growing tissues. , 2000, Physical review letters.

[14]  Giovanni Parmigiani,et al.  Prevalence of somatic alterations in the colorectal cancer cell genome , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Dominik Wodarz,et al.  Drug resistance in cancer: principles of emergence and prevention. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[16]  R. Arceci Intratumor Heterogeneity and Branched Evolution Revealed by Multiregion Sequencing , 2012 .

[17]  Abbas Samani,et al.  Measuring the elastic modulus of ex vivo small tissue samples. , 2003, Physics in medicine and biology.

[18]  Philip Greulich,et al.  Differentiation imbalance in single oesophageal progenitor cells causes clonal immortalization and field change , 2014, Nature Cell Biology.

[19]  C. Sander,et al.  Tumor Genetic Analyses of Patients with Metastatic Renal Cell Carcinoma and Extended Benefit from mTOR Inhibitor Therapy , 2014, Clinical Cancer Research.

[20]  J H Goldie,et al.  The genetic origin of drug resistance in neoplasms: implications for systemic therapy. , 1984, Cancer research.

[21]  Cynthia A. Reinhart-King,et al.  Physical biology in cancer. 5. The rocky road of metastasis: the role of cytoskeletal mechanics in cell migratory response to 3D matrix topography. , 2014, American journal of physiology. Cell physiology.

[22]  K. Foster,et al.  A Quantitative Test of Population Genetics Using Spatiogenetic Patterns in Bacterial Colonies , 2011, The American Naturalist.

[23]  Cass T. Miller,et al.  A multiphase model for three-dimensional tumor growth , 2013, New journal of physics.

[24]  N L Komarova,et al.  Spatial invasion dynamics on random and unstructured meshes: implications for heterogeneous tumor populations. , 2014, Journal of theoretical biology.

[25]  H Honda,et al.  Establishment of epidermal cell columns in mammalian skin: computer simulation. , 1979, Journal of theoretical biology.

[26]  H. Ohtsuki,et al.  Accumulation of driver and passenger mutations during tumor progression , 2009, Proceedings of the National Academy of Sciences.

[27]  R. Weinberg,et al.  The Biology of Cancer , 2006 .

[28]  C. Maley,et al.  Spatial structure increases the waiting time for cancer , 2011, New journal of physics.

[29]  S. Torquato,et al.  Diversity of dynamics and morphologies of invasive solid tumors. , 2012, AIP advances.

[30]  Z. Stachura,et al.  Elasticity of normal and cancerous human bladder cells studied by scanning force microscopy , 1999, European Biophysics Journal.

[31]  J. Rieu,et al.  Fine Tuning of Tissues' Viscosity and Surface Tension through Contractility Suggests a New Role for α-Catenin , 2013, PloS one.

[32]  Martin A. Nowak,et al.  Genetic Progression and the Waiting Time to Cancer , 2007, PLoS Comput. Biol..

[33]  Krishnendu Chatterjee,et al.  Evolutionary dynamics of cancer in response to targeted combination therapy , 2013, eLife.

[34]  S Torquato,et al.  Emergence of a Subpopulation in a Computational Model of Tumor Growth , 2022 .

[35]  Natalia L Komarova,et al.  Selection in spatial stochastic models of cancer: Migration as a key modulator of fitness , 2010, Biology Direct.

[36]  A. Anderson,et al.  A hybrid mathematical model of solid tumour invasion: the importance of cell adhesion , 2005 .

[37]  J H Goldie,et al.  A stochastic model for the origin and treatment of tumors containing drug-resistant cells. , 1986, Bulletin of mathematical biology.

[38]  E Schöll,et al.  Comparing the growth kinetics of cell populations in two and three dimensions. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[39]  R. Solé,et al.  Metapopulation dynamics and spatial heterogeneity in cancer , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[40]  D. Nelson,et al.  Survival probabilities at spherical frontiers. , 2014, Theoretical population biology.

[41]  Françoise Brochard-Wyart,et al.  Soft Matter Models of Developing Tissues and Tumors , 2012, Science.

[42]  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.

[43]  Alexander R A Anderson,et al.  Quantifying the Role of Angiogenesis in Malignant Progression of Gliomas: in Silico Modeling Integrates Imaging and Histology Nih Public Access Author Manuscript Introduction , 2011 .

[44]  R. Kalluri,et al.  Corrigendum: PGC-1α mediates mitochondrial biogenesis and oxidative phosphorylation in cancer cells to promote metastasis , 2014, Nature Cell Biology.

[45]  J. Goldie,et al.  A model for the resistance of tumor cells to cancer chemotherapeutic agents , 1983 .

[46]  P Hahnfeldt,et al.  Migration rules: tumours are conglomerates of self-metastases , 2009, British Journal of Cancer.

[47]  D. Hartl,et al.  Principles of population genetics , 1981 .

[48]  Melanie J. I. Müller,et al.  Selective sweeps in growing microbial colonies. , 2012, Physical biology.

[49]  Guillaume Charras,et al.  Physical influences of the extracellular environment on cell migration , 2014, Nature Reviews Molecular Cell Biology.

[50]  Luca Toschi,et al.  Preexistence and clonal selection of MET amplification in EGFR mutant NSCLC. , 2010, Cancer cell.

[51]  M. Nowak,et al.  Spatial evolution of tumors with successive driver mutations. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[52]  Dominik Wodarz,et al.  Tumor growth dynamics: insights into evolutionary processes. , 2013, Trends in ecology & evolution.

[53]  N. Komarova Spatial interactions and cooperation can change the speed of evolution of complex phenotypes , 2014, Proceedings of the National Academy of Sciences.

[54]  Ignacio Ramis-Conde,et al.  Multi-scale modelling of cancer cell intravasation: the role of cadherins in metastasis , 2009, Physical biology.

[55]  Z. Szallasi,et al.  Clonal status of actionable driver events and the timing of mutational processes in cancer evolution , 2015, Science Translational Medicine.

[56]  Alexander A. Alemi,et al.  Mechanical Properties of Growing Melanocytic Nevi and the Progression to Melanoma , 2014, PloS one.

[57]  Cristian Picioreanu,et al.  iDynoMiCS: next-generation individual-based modelling of biofilms. , 2011, Environmental microbiology.

[58]  D. Drasdo,et al.  A single-cell-based model of tumor growth in vitro: monolayers and spheroids , 2005, Physical biology.

[59]  M. Nowak,et al.  Distant Metastasis Occurs Late during the Genetic Evolution of Pancreatic Cancer , 2010, Nature.

[60]  K. Kinzler,et al.  Cancer Genome Landscapes , 2013, Science.

[61]  Richard Durrett,et al.  Spatial Moran Models II. Tumor growth and progression , 2012 .

[62]  Natalia L. Komarova,et al.  Accelerated crossing of fitness valleys through division of labor and cheating in asexual populations , 2012, Scientific Reports.

[63]  M. Eden A Two-dimensional Growth Process , 1961 .

[64]  F. Michor,et al.  Evolution of acquired resistance to anti-cancer therapy. , 2014, Journal of theoretical biology.

[65]  Zhong-yu Wang,et al.  New Insight into Ki67 Expression at the Invasive Front in Breast Cancer , 2013, PloS one.

[66]  A. Feinberg,et al.  Genome-scale epigenetic reprogramming during epithelial to mesenchymal transition , 2011, Nature Structural &Molecular Biology.

[67]  K. Leder,et al.  Multifocality and recurrence risk: a quantitative model of field cancerization. , 2013, Journal of theoretical biology.

[68]  R. Durrett,et al.  Spatial Moran models, II: cancer initiation in spatially structured tissue , 2016, Journal of mathematical biology.

[69]  V. P. Collins,et al.  Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics , 2013, Proceedings of the National Academy of Sciences.

[70]  Sridhar Ramaswamy,et al.  Circulating Tumor Cell Clusters Are Oligoclonal Precursors of Breast Cancer Metastasis , 2014, Cell.

[71]  K. Burridge,et al.  The on-off relationship of Rho and Rac during integrin-mediated adhesion and cell migration , 2014, Small GTPases.

[72]  Jelena Pjesivac-Grbovic,et al.  A multiscale model for avascular tumor growth. , 2005, Biophysical journal.

[73]  J. Troge,et al.  Tumour evolution inferred by single-cell sequencing , 2011, Nature.

[74]  B. Waclaw,et al.  Mutational pathway determines whether drug gradients accelerate evolution of drug-resistant cells. , 2012, Physical review letters.

[75]  S. McDougall,et al.  Mathematical modelling of dynamic adaptive tumour-induced angiogenesis: clinical implications and therapeutic targeting strategies. , 2006, Journal of theoretical biology.

[76]  Martin A. Nowak,et al.  Comparative lesion sequencing provides insights into tumor evolution , 2008, Proceedings of the National Academy of Sciences.

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

[78]  J. Bergh,et al.  Clinically used breast cancer markers such as estrogen receptor, progesterone receptor, and human epidermal growth factor receptor 2 are unstable throughout tumor progression. , 2012, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[79]  Jacques Prost,et al.  Undulation instability of epithelial tissues. , 2011, Physical review letters.

[80]  S. McWeeney,et al.  Cancer stem cell tumor model reveals invasive morphology and increased phenotypical heterogeneity. , 2010, Cancer research.

[81]  Johannes G. Reiter,et al.  The molecular evolution of acquired resistance to targeted EGFR blockade in colorectal cancers , 2012, Nature.

[82]  N. Komarova Spatial stochastic models of cancer: fitness, migration, invasion. , 2013, Mathematical biosciences and engineering : MBE.

[83]  S Torquato,et al.  Simulated brain tumor growth dynamics using a three-dimensional cellular automaton. , 2000, Journal of theoretical biology.

[84]  S. Spencer,et al.  An ordinary differential equation model for the multistep transformation to cancer. , 2004, Journal of theoretical biology.

[85]  Michael J. Wiser,et al.  Long-Term Dynamics of Adaptation in Asexual Populations , 2013, Science.

[86]  Philip Gerlee,et al.  The Impact of Phenotypic Switching on Glioblastoma Growth and Invasion , 2012, PLoS Comput. Biol..

[87]  H. Othmer,et al.  A HYBRID MODEL FOR TUMOR SPHEROID GROWTH IN VITRO I: THEORETICAL DEVELOPMENT AND EARLY RESULTS , 2007 .

[88]  T. Vincent,et al.  An evolutionary model of carcinogenesis. , 2003, Cancer research.

[89]  E. Somfai,et al.  Reproduction-time statistics and segregation patterns in growing populations. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[90]  S. Hoehme,et al.  On the Role of Physics in the Growth and Pattern Formation of Multi-Cellular Systems: What can we Learn from Individual-Cell Based Models? , 2007 .

[91]  Neil F. Johnson,et al.  Model for in vivo progression of tumors based on co-evolving cell population and vasculature , 2011, Scientific reports.

[92]  Alan Hall,et al.  Rho family GTPases. , 2012, Biochemical Society transactions.

[93]  A. Frampton,et al.  Gene of the month: E-cadherin (CDH1) , 2013, Journal of Clinical Pathology.

[94]  Martin A Nowak,et al.  Timing and heterogeneity of mutations associated with drug resistance in metastatic cancers , 2014, Proceedings of the National Academy of Sciences.

[95]  Jan Poleszczuk,et al.  Evolution and Phenotypic Selection of Cancer Stem Cells , 2015, PLoS Comput. Biol..

[96]  S. Torquato Toward an Ising model of cancer and beyond , 2010, Physical biology.

[97]  Richard Durrett,et al.  Evolution of Resistance and Progression to Disease during Clonal Expansion of Cancer , 2009 .

[98]  T. C. White,et al.  Increased mRNA levels of ERG16, CDR, and MDR1 correlate with increases in azole resistance in Candida albicans isolates from a patient infected with human immunodeficiency virus , 1997, Antimicrobial agents and chemotherapy.

[99]  Toby C. Cornish,et al.  A Single Institution's 26-Year Experience With Nonfunctional Pancreatic Neuroendocrine Tumors: A Validation of Current Staging Systems and a New Prognostic Nomogram , 2014, Annals of surgery.

[100]  K. Swanson,et al.  Modeling Tumor-Associated Edema in Gliomas during Anti-Angiogenic Therapy and Its Impact on Imageable Tumor , 2013, Front. Oncol..

[101]  I. Fidler,et al.  AACR centennial series: the biology of cancer metastasis: historical perspective. , 2010, Cancer research.

[102]  Deena R. Schmidt,et al.  A waiting time problem arising from the study of multi-stage carcinogenesis , 2007, 0707.2057.

[103]  P. Maini,et al.  Mathematical modeling of cell population dynamics in the colonic crypt and in colorectal cancer , 2007, Proceedings of the National Academy of Sciences.

[104]  Martin A Nowak,et al.  Dynamics of targeted cancer therapy. , 2012, Trends in molecular medicine.

[105]  Krishnendu Chatterjee,et al.  The effect of one additional driver mutation on tumor progression , 2012, Evolutionary applications.