The Hybrid Automata Library (HAL) is a Java Library developed for use in mathematical oncology modeling. It is made of simple, efficient, generic components that can be used to model complex spatial systems. HAL's components can broadly be classified into: on-and off-lattice agent containers, finite difference diffusion fields, a GUI building system, and additional tools and utilities for computation and data collection. These components are designed to operate independently and are standardized to make them easy to interface with one another. As a demonstration of how modeling can be simplified using our approach, we have included a complete example of a hybrid model (a spatial model with interacting agent-based and PDE components). HAL is a useful asset for researchers who wish to build performant 1D, 2D and 3D hybrid models in Java, while not starting entirely from scratch. It is available on GitHub at https://github.com/MathOnco/HAL under the MIT License. HAL requires the Java JDK version 1.8 or later to compile and run the source code.

Nature Research wishes to improve the reproducibility of the work that we publish. This form provides structure for consistency and transparency in reporting. For further information on Nature Research policies, see our Editorial Policies and the Editorial Policy Checklist.

The biology of the metastatic colonization process remains a poorly understood phenomenon. To improve our knowledge of its dynamics, we conducted a modelling study based on multi-modal data from an orthotopic murine experimental system of metastatic renal cell carcinoma. The standard theory of metastatic colonization usually assumes that secondary tumours, once established at a distant site, grow independently from each other and from the primary tumour. Using a mathematical model that translates this assumption into equations, we challenged this theory against our data that included: 1) dynamics of primary tumour cells in the kidney and metastatic cells in the lungs, retrieved by green fluorescent protein tracking, and 2) magnetic resonance images (MRI) informing on the number and size of macroscopic lesions. Critically, when calibrated on the growth of the primary tumour and total metastatic burden, the predicted theoretical size distributions were not in agreement with the MRI observations. Moreover, tumour expansion only based on proliferation was not able to explain the volume increase of the metastatic lesions. These findings strongly suggested rejection of the standard theory, demonstrating that the time development of the size distribution of metastases could not be explained by independent growth of metastatic foci. This led us to investigate the effect of spatial interactions between merging metastatic tumours on the dynamics of the global metastatic burden. We derived a mathematical model of spatial tumour growth, confronted it with experimental data of single metastatic tumour growth, and used it to provide insights on the dynamics of multiple tumours growing in close vicinity. Together, our results have implications for theories of the metastatic process and suggest that global dynamics of metastasis development is dependent on spatial interactions between metastatic lesions.