A Bayesian Sequential Learning Framework to Parameterise Continuum Models of Melanoma Invasion into Human Skin

We present a novel framework to parameterise a mathematical model of cell invasion that describes how a population of melanoma cells invades into human skin tissue. Using simple experimental data extracted from complex experimental images, we estimate three model parameters: (i) the melanoma cell proliferation rate, $$\lambda $$λ; (ii) the melanoma cell diffusivity, D; and (iii) $$\delta $$δ, a constant that determines the rate that melanoma cells degrade the skin tissue. The Bayesian sequential learning framework involves a sequence of increasingly sophisticated experimental data from: (i) a spatially uniform cell proliferation assay; (ii) a two-dimensional circular barrier assay; and (iii) a three-dimensional invasion assay. The Bayesian sequential learning approach leads to well-defined parameter estimates. In contrast, taking a naive approach that attempts to estimate all parameters from a single set of images from the same experiment fails to produce meaningful results. Overall, our approach to inference is simple-to-implement, computationally efficient, and well suited for many cell biology phenomena that can be described by low-dimensional continuum models using ordinary differential equations and partial differential equations. We anticipate that this Bayesian sequential learning framework will be relevant in other biological contexts where it is challenging to extract detailed, quantitative biological measurements from experimental images and so we must rely on using relatively simple measurements from complex images.

[1]  H M Byrne,et al.  An ordinary differential equation model for full thickness wounds and the effects of diabetes. , 2014, Journal of theoretical biology.

[2]  Matthew J. Simpson,et al.  Nonmonotone chemotactic invasion: High-resolution simulations, phase plane analysis and new benchmark problems , 2007, J. Comput. Phys..

[3]  E. T. Gawlinski,et al.  A reaction-diffusion model of cancer invasion. , 1996, Cancer research.

[4]  G J Pettet,et al.  Modelling the action of proteinase and inhibitor in tissue invasion. , 1998, Mathematical biosciences.

[5]  Geoffrey W. Stevens,et al.  Cell migration and proliferation during monolayer formation and wound healing , 2009 .

[6]  Mark M. Tanaka,et al.  Sequential Monte Carlo without likelihoods , 2007, Proceedings of the National Academy of Sciences.

[7]  P. Maini,et al.  The role of acidity in solid tumour growth and invasion. , 2005, Journal of theoretical biology.

[8]  John K Kruschke,et al.  Bayesian data analysis. , 2010, Wiley interdisciplinary reviews. Cognitive science.

[9]  Matthew J Simpson,et al.  Inferring parameters for a lattice-free model of cell migration and proliferation using experimental data , 2017, bioRxiv.

[10]  Alexander G. Fletcher,et al.  A hierarchical Bayesian model for understanding the spatiotemporal dynamics of the intestinal epithelium , 2017, PLoS Comput. Biol..

[11]  L. G. Pillis,et al.  A Comparison and Catalog of Intrinsic Tumor Growth Models , 2013, Bulletin of mathematical biology.

[12]  S. McCue,et al.  A Bayesian Computational Approach to Explore the Optimal Duration of a Cell Proliferation Assay , 2017, Bulletin of Mathematical Biology.

[13]  Matthew J Simpson,et al.  Quantifying the efficacy of first aid treatments for burn injuries using mathematical modelling and in vivo porcine experiments , 2017, Scientific Reports.

[14]  Matthew J. Simpson,et al.  Sensitivity of Edge Detection Methods for Quantifying Cell Migration Assays , 2013, PloS one.

[15]  Kerry A Landman,et al.  Multi-scale modeling of a wound-healing cell migration assay. , 2007, Journal of theoretical biology.

[16]  Matthew J. Simpson,et al.  Quantitative comparison of the spreading and invasion of radial growth phase and metastatic melanoma cells in a three-dimensional human skin equivalent model , 2017, PeerJ.

[17]  David Gavaghan,et al.  Inference-based assessment of parameter identifiability in nonlinear biological models , 2018, Journal of The Royal Society Interface.

[18]  Ruth E. Baker,et al.  Multilevel rejection sampling for approximate Bayesian computation , 2017, Comput. Stat. Data Anal..

[19]  Matthew J Simpson,et al.  Quantifying rates of cell migration and cell proliferation in co-culture barrier assays reveals how skin and melanoma cells interact during melanoma spreading and invasion. , 2017, Journal of theoretical biology.

[20]  D. L. Sean McElwain,et al.  Interpreting scratch assays using pair density dynamics and approximate Bayesian computation , 2014, Open Biology.

[21]  J. Sherratt,et al.  A two parameter family of travelling waves with a singular barrier arising from the modelling of extracellular matrix mediated cellular invasion , 1999 .

[22]  Joe Pitt-Francis,et al.  Bayesian Calibration, Validation and Uncertainty Quantification for Predictive Modelling of Tumour Growth: A Tutorial , 2017, Bulletin of Mathematical Biology.

[23]  V. Quaranta,et al.  Integrative mathematical oncology , 2008, Nature Reviews Cancer.

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

[25]  John Norbury,et al.  Discontinuous travelling wave solutions for certain hyperbolic systems , 2002 .

[26]  Matthew J Simpson,et al.  Optimal Quantification of Contact Inhibition in Cell Populations. , 2017, Biophysical journal.

[27]  D. L. Sean McElwain,et al.  Three-dimensional experiments and individual based simulations show that cell proliferation drives melanoma nest formation in human skin tissue , 2018, bioRxiv.

[28]  D. L. Sean McElwain,et al.  Travelling waves in a wound healing assay , 2004, Appl. Math. Lett..

[29]  Graeme J. Pettet,et al.  Chemotactic Cellular Migration: Smooth and Discontinuous Travelling Wave Solutions , 2003, SIAM J. Appl. Math..

[30]  Matthew J Simpson,et al.  Mathematical models for cell migration with real-time cell cycle dynamics , 2017, bioRxiv.

[31]  Antonio Fasano,et al.  Slow and fast invasion waves in a model of acid-mediated tumour growth. , 2009, Mathematical biosciences.

[32]  A. J. Perumpanani,et al.  Traveling Shock Waves Arising in a Model of Malignant Invasion , 1999, SIAM J. Appl. Math..

[33]  Paul Fearnhead,et al.  Constructing summary statistics for approximate Bayesian computation: semi‐automatic approximate Bayesian computation , 2012 .

[34]  Philip K Maini,et al.  Traveling wave model to interpret a wound-healing cell migration assay for human peritoneal mesothelial cells. , 2004, Tissue engineering.

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

[36]  H. Byrne,et al.  Biphasic behaviour in malignant invasion. , 2006, Mathematical medicine and biology : a journal of the IMA.

[37]  K. A. Landman,et al.  Diffusive and Chemotactic Cellular Migration: Smooth and Discontinuous Traveling Wave Solutions , 2005, SIAM J. Appl. Math..

[38]  M. Meyer-Hermann,et al.  Exploiting the Synergy between Carboplatin and ABT-737 in the Treatment of Ovarian Carcinomas , 2014, PloS one.

[39]  Graeme J. Pettet,et al.  Tactically-driven nonmonotone travelling waves , 2008 .

[40]  OpenStax Introductory Statistics , 2013 .

[41]  David Welch,et al.  Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems , 2009, Journal of The Royal Society Interface.

[42]  Stephen Lynch,et al.  The Image Processing Toolbox , 2014 .

[43]  Ruth E. Baker,et al.  Multiple types of data are required to identify the mechanisms influencing the spatial expansion of melanoma cell colonies , 2013, BMC Systems Biology.

[44]  C. Please,et al.  Experimental characterization and computational modelling of two-dimensional cell spreading for skeletal regeneration , 2007, Journal of The Royal Society Interface.

[45]  Peter Canoll,et al.  Glial progenitor cell recruitment drives aggressive glioma growth: mathematical and experimental modelling , 2012, Journal of The Royal Society Interface.

[46]  J. Sherratt,et al.  Travelling wave solutions to a haptotaxis-dominated model of malignant invasion , 2001 .

[47]  J. Møller Discussion on the paper by Feranhead and Prangle , 2012 .