Warlock: an automated computational workflow for simulating spatially structured tumour evolution

Abstract

A primary goal of modern cancer research is to characterize tumour growth and evolution, to improve clinical forecasting and individualized treatment. Agent-based models support this endeavour but existing models either oversimplify spatial structure or are mathematically intractable. Here we present warlock, an open-source automated computational workflow for fast, efficient simulation of intratumour population genetics in any of a diverse set of spatial structures. Warlock encapsulates a deme-based oncology model (demon), designed to bridge the divide between agent-based simulations and analytical population genetics models, such as the spatial Moran process. Model output can be readily compared to multi-region and single-cell sequencing data for model selection or biological parameter inference. An interface for High Performance Computing permits hundreds of simulations to be run in parallel. We discuss prior applications of this workflow to investigating human cancer evolution.