Cancer cell population growth kinetics at low densities deviate from the exponential growth model and suggest an Allee effect

Models of cancer cell population expansion assume exponential growth kinetics at low cell densities, with deviations from exponential growth only at higher densities due to limited resources such as space and nutrients. However, recent pre-clinical and clinical observations of tumor initiation or recurrence indicate the presence of tumor growth kinetics in which growth rates scale with cell numbers. These observations are analogous to the cooperative behavior of species in an ecosystem described by the ecological principle of the Allee effect. In preclinical and clinical models however, tumor growth data is limited by the lower limit of detection (i.e. a measurable lesion) and confounding variables, such as tumor microenvironment and immune responses may cause and mask deviations from exponential growth models. In this work, we present alternative growth models to investigate the presence of an Allee effect in cancer cells seeded at low cell densities in a controlled in vitro setting. We propose a stochastic modeling framework to consider the small number of cells in this low-density regime and use the moment approach for stochastic parameter estimation to calibrate the stochastic growth trajectories. We validate the framework on simulated data and apply this approach to longitudinal cell proliferation data of BT-474 luminal B breast cancer cells. We find that cell population growth kinetics are best described by a model structure that considers the Allee effect, in that the birth rate of tumor cells depends on cell number. This indicates a potentially critical role of cooperative behavior among tumor cells at low cell densities with relevance to early stage growth patterns of emerging tumors and relapse. Author Summary The growth kinetics of cancer cells at very low cell densities are of utmost clinical importance as the ability of a small number of newly transformed or surviving cells to grow exponentially and thus, to “take off” underlies tumor formation and relapse after treatment. Mathematical models of stochastic tumor cell growth typically assume a stochastic birth-death process of cells impacted by limited nutrients and space when cells reach high density, resulting in the widely accepted logistic growth model. Here we present an in-depth investigation of alternate growth models adopted from ecology to describe potential deviations from a simple cell autonomous birth-death model at low cell densities. We show that our stochastic modeling framework is robust and can be used to identify the underlying structure of stochastic growth trajectories from both simulated and experimental data taken from a controlled in vitro setting in which we can capture data from the relevant low cell density regime. This work suggests that the assumption of cell autonomous proliferation via a constant exponential growth rate at low cell densities may not be appropriate for all cancer cell growth dynamics. Consideration of cooperative behavior amongst tumor cells in this regime is critical for elucidating strategies for controlling tumor cell growth.

[1]  N. Ramankutty,et al.  Science–graphic art partnerships to increase research impact , 2019, Communications Biology.

[2]  Kang Yao,et al.  Quorum Sensing: A Prospective Therapeutic Target for Bacterial Diseases , 2019, BioMed research international.

[3]  Gregory J. Kimmel,et al.  Neighborhood size-effects shape growing population dynamics in evolutionary public goods games , 2019, Communications Biology.

[4]  P. Newton,et al.  Cellular interactions constrain tumor growth , 2019, Proceedings of the National Academy of Sciences.

[5]  George D. Cresswell,et al.  Evolutionary dynamics of residual disease in human glioblastoma , 2018, Annals of oncology : official journal of the European Society for Medical Oncology.

[6]  Douglas A. Lauffenburger,et al.  Analysis of Single-Cell RNA-Seq Identifies Cell-Cell Communication Associated with Tumor Characteristics , 2018, Cell reports.

[7]  Jan Poleszczuk,et al.  The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity , 2018, International journal of molecular sciences.

[8]  A. Schäffer,et al.  The evolution of single cell-derived colorectal cancer cell lines is dominated by the continued selection of tumor-specific genomic imbalances, despite random chromosomal instability , 2018, Carcinogenesis.

[9]  Jan Hasenauer,et al.  A Hierarchical, Data-Driven Approach to Modeling Single-Cell Populations Predicts Latent Causes of Cell-To-Cell Variability. , 2018, Cell systems.

[10]  K. Pienta,et al.  Cancer as a Social Dysfunction—Why Cancer Research Needs New Thinking , 2018, Molecular Cancer Research.

[11]  Thomas E. Yankeelov,et al.  Precision Medicine with Imprecise Therapy: Computational Modeling for Chemotherapy in Breast Cancer , 2018, Translational oncology.

[12]  Yohsuke Yagawa,et al.  Systematic Screening of Chemokines to Identify Candidates to Model and Create Ectopic Lymph Node Structures for Cancer Immunotherapy , 2017, Scientific Reports.

[13]  Andras Czirok,et al.  The role of Allee effect in modelling post resection recurrence of glioblastoma , 2017, PLoS Comput. Biol..

[14]  P. Altrock,et al.  Extinction rates in tumour public goods games , 2017, Journal of The Royal Society Interface.

[15]  David Basanta,et al.  Fibroblasts and Alectinib switch the evolutionary games played by non-small cell lung cancer , 2017, Nature Ecology & Evolution.

[16]  Sui Huang,et al.  Extracting Intercellular Signaling Network of Cancer Tissues using Ligand-Receptor Expression Patterns from Whole-tumor and Single-cell Transcriptomes , 2017, Scientific Reports.

[17]  Rafael Meza,et al.  A Systematic Approach to Determining the Identifiability of Multistage Carcinogenesis Models , 2017, Risk analysis : an official publication of the Society for Risk Analysis.

[18]  Mohit Kumar Jolly,et al.  Epithelial/mesenchymal plasticity: how have quantitative mathematical models helped improve our understanding? , 2017, Molecular oncology.

[19]  S. Hanash,et al.  Epithelial-mesenchymal plasticity: How have quantitative mathematical models helped improve our understanding? , 2017, bioRxiv.

[20]  K. Murao,et al.  A method for evaluating the performance of computer-aided detection of pulmonary nodules in lung cancer CT screening: detection limit for nodule size and density. , 2017, The British journal of radiology.

[21]  I. Bose,et al.  Allee dynamics: growth, extinction and range expansion , 2017, bioRxiv.

[22]  Yukyung Jun,et al.  Rapid emergence and mechanisms of resistance by U87 glioblastoma cells to doxorubicin in an in vitro tumor microfluidic ecology , 2016, Proceedings of the National Academy of Sciences.

[23]  T E Yankeelov,et al.  Selection, calibration, and validation of models of tumor growth. , 2016, Mathematical models & methods in applied sciences : M3AS.

[24]  Owen L. Petchey,et al.  The effects of demographic stochasticity and parameter uncertainty on predicting the establishment of introduced species , 2016, Ecology and evolution.

[25]  Joel s. Brown,et al.  Ecological paradigms to understand the dynamics of metastasis. , 2016, Cancer letters.

[26]  G. Genin,et al.  Multi-scale Modeling in Clinical Oncology: Opportunities and Barriers to Success , 2016, Annals of Biomedical Engineering.

[27]  Fabian J. Theis,et al.  Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion , 2016, PLoS Comput. Biol..

[28]  M. Gottesman,et al.  Mathematical Modeling Reveals That Changes to Local Cell Density Dynamically Modulate Baseline Variations in Cell Growth and Drug Response. , 2016, Cancer research.

[29]  J. Lowengrub,et al.  Feedback Regulation in a Cancer Stem Cell Model can Cause an Allee Effect , 2016, Bulletin of mathematical biology.

[30]  L. Sewalt,et al.  Influences of Allee effects in the spreading of malignant tumours. , 2016, Journal of theoretical biology.

[31]  Roland R. Regoes,et al.  Investigating the Consequences of Interference between Multiple CD8+ T Cell Escape Mutations in Early HIV Infection , 2016, PLoS Comput. Biol..

[32]  Edda Klipp,et al.  A Thermodynamic Model of Monovalent Cation Homeostasis in the Yeast Saccharomyces cerevisiae , 2016, PLoS Comput. Biol..

[33]  Paul K. Newton,et al.  An Evolutionary Model of Tumor Cell Kinetics and the Emergence of Molecular Heterogeneity Driving Gompertzian Growth , 2015, SIAM Rev..

[34]  A. Martinez,et al.  Models for Allee effect based on physical principles. , 2015, Journal of theoretical biology.

[35]  D. Sargent,et al.  Evaluating Continuous Tumor Measurement-Based Metrics as Phase II Endpoints for Predicting Overall Survival. , 2015, Journal of the National Cancer Institute.

[36]  Andreas Deutsch,et al.  An Emerging Allee Effect Is Critical for Tumor Initiation and Persistence , 2015, PLoS Comput. Biol..

[37]  A. Pisco,et al.  Non-genetic cancer cell plasticity and therapy-induced stemness in tumour relapse: ‘What does not kill me strengthens me' , 2015, British Journal of Cancer.

[38]  Peter Kraft,et al.  Identification of Novel Genetic Markers of Breast Cancer Survival , 2015, Journal of the National Cancer Institute.

[39]  R. Marschalek,et al.  Optimized Sleeping Beauty transposons rapidly generate stable transgenic cell lines. , 2015, Biotechnology journal.

[40]  K. Pienta,et al.  Ecology meets cancer biology: The cancer swamp promotes the lethal cancer phenotype , 2015, Oncotarget.

[41]  M. Sherman,et al.  The case for conducting a randomized clinical trial to assess the efficacy of a single dose of prophylactic HPV vaccines among adolescents. , 2015, Journal of the National Cancer Institute.

[42]  M. Archetti,et al.  Heterogeneity for IGF-II production maintained by public goods dynamics in neuroendocrine pancreatic cancer , 2015, Proceedings of the National Academy of Sciences.

[43]  Hong Qian,et al.  Nonequilibrium Population Dynamics of Phenotype Conversion of Cancer Cells , 2014, PloS one.

[44]  T. Blackburn,et al.  Quantifying invasion risk: the relationship between establishment probability and founding population size , 2014 .

[45]  Seth Sullivant,et al.  Identifiability Results for Several Classes of Linear Compartment Models , 2014, Bulletin of Mathematical Biology.

[46]  Xuguo Zhou,et al.  Validation of Reference Housekeeping Genes for Gene Expression Studies in Western Corn Rootworm (Diabrotica virgifera virgifera) , 2014, PloS one.

[47]  K. Polyak,et al.  Non-cell autonomous tumor-growth driving supports sub-clonal heterogeneity , 2014, Nature.

[48]  John M. L. Ebos,et al.  Classical Mathematical Models for Description and Prediction of Experimental Tumor Growth , 2014, PLoS Comput. Biol..

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

[50]  S. Gestl,et al.  Tumor cell heterogeneity maintained by cooperating subclones in Wnt-driven mammary cancers , 2014, Nature.

[51]  F. De Filippis,et al.  A Selected Core Microbiome Drives the Early Stages of Three Popular Italian Cheese Manufactures , 2014, PloS one.

[52]  Meike J. Wittmann,et al.  Genetic Diversity in Introduced Populations with an Allee Effect , 2013, Genetics.

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

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

[55]  Donald E Ingber,et al.  Epoxyeicosanoids stimulate multiorgan metastasis and tumor dormancy escape in mice. , 2012, The Journal of clinical investigation.

[56]  E. Lander,et al.  Theory Stochastic State Transitions Give Rise to Phenotypic Equilibrium in Populations of Cancer Cells , 2011 .

[57]  K. Pienta,et al.  Modeling invasion of metastasizing cancer cells to bone marrow utilizing ecological principles , 2011, Theoretical Biology and Medical Modelling.

[58]  Christian P Robert,et al.  Lack of confidence in approximate Bayesian computation model choice , 2011, Proceedings of the National Academy of Sciences.

[59]  E. Lander,et al.  Stochastic State Transitions Give Rise to Phenotypic Equilibrium in Populations of Cancer Cells , 2011, Cell.

[60]  Sophia Hsin-Jung Li,et al.  Paracrine and Autocrine Signals Induce and Maintain Mesenchymal and Stem Cell States in the Breast , 2011, Cell.

[61]  Gordon Broderick,et al.  Using an agent-based model to analyze the dynamic communication network of the immune response , 2011, Theoretical Biology and Medical Modelling.

[62]  Ursula Klingmüller,et al.  Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood , 2009, Bioinform..

[63]  Boris Jerchow,et al.  Molecular evolution of a novel hyperactive Sleeping Beauty transposase enables robust stable gene transfer in vertebrates , 2009, Nature Genetics.

[64]  Hannah H. Chang,et al.  Non-genetic heterogeneity — a mutation-independent driving force for the somatic evolution of tumours , 2009, Nature Reviews Genetics.

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

[66]  S. Ranade,et al.  Stem cell transcriptome profiling via massive-scale mRNA sequencing , 2008, Nature Methods.

[67]  J. Gascoigne,et al.  Allee Effects in Ecology and Conservation , 2008 .

[68]  M. Berens,et al.  Autocrine factors that sustain glioma invasion and paracrine biology in the brain microenvironment. , 2007, Journal of the National Cancer Institute.

[69]  K. Pienta,et al.  Evolution of cooperation among tumor cells , 2006, Proceedings of the National Academy of Sciences.

[70]  E. Wagenmakers,et al.  AIC model selection using Akaike weights , 2004, Psychonomic bulletin & review.

[71]  D. Balding,et al.  Approximate Bayesian computation in population genetics. , 2002, Genetics.

[72]  David Eisenberg,et al.  Bioinformatic identification of potential autocrine signaling loops in cancers from gene expression profiles , 2001, Nature Genetics.

[73]  D. Gillespie The chemical Langevin equation , 2000 .

[74]  A. Raftery Bayes Factors and BIC , 1999 .

[75]  R. Gatenby,et al.  Population ecology issues in tumor growth. , 1991, Cancer research.

[76]  D. Raghavan,et al.  Clonal analysis of a bladder cancer cell line: an experimental model of tumour heterogeneity. , 1990, British Journal of Cancer.

[77]  L. Norton A Gompertzian model of human breast cancer growth. , 1988, Cancer research.

[78]  M. Retsky,et al.  A stochastic numerical model of breast cancer growth that simulates clinical data. , 1984, Cancer research.

[79]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[80]  E. Pacheco A review of models for cancer chemotherapy based on Optimal Control Internship Report INESC ID – Technical report 16 / 2016 , 2017 .

[81]  B. Houchmandzadeh Extracting moments from Master Equations , 2009 .

[82]  J. Neyman,et al.  Stochastic models of population dynamics. , 1959, Science.

[83]  C. Winsor,et al.  The Gompertz Curve as a Growth Curve. , 1932, Proceedings of the National Academy of Sciences of the United States of America.