Investigation into the role of macrophages heterogeneity on solid tumour aggregations

Quantification of spatial and phenotypic heterogeneity in an agent-based model of tumour-macrophage interactions

We introduce a new spatial statistic, the weighted pair correlation function (wPCF). The wPCF extends the existing pair correlation function (PCF) and cross-PCF to describe spatial relationships between points marked with combinations of discrete and continuous labels. We validate its use through application to a new agent-based model (ABM) which simulates interactions between macrophages and tumour cells. These interactions are influenced by the spatial positions of the cells and by macrophage phenotype, a continuous variable that ranges from anti-tumour to pro-tumour. By varying model parameters that regulate macrophage phenotype, we show that the ABM exhibits behaviours which resemble the 'three Es of cancer immunoediting': Equilibrium, Escape, and Elimination. We use the wPCF to analyse synthetic images generated by the ABM. We show that the wPCF generates a 'human readable' statistical summary of where macrophages with different phenotypes are located relative to both blood vessels and tumour cells. We also define a distinct 'PCF signature' that characterises each of the three Es of immunoediting, by combining wPCF measurements with the cross-PCF describing interactions between vessels and tumour cells. By applying dimension reduction techniques to this signature, we identify its key features and train a support vector machine classifier to distinguish between simulation outputs based on their PCF signature. This proof-of-concept study shows how multiple spatial statistics can be combined to analyse the complex spatial features that the ABM generates, and to partition them into interpretable groups. The intricate spatial features produced by the ABM are similar to those generated by state-of-the-art multiplex imaging techniques which distinguish the spatial distribution and intensity of multiple biomarkers in biological tissue regions. Applying methods such as the wPCF to multiplex imaging data would exploit the continuous variation in biomarker intensities and generate more detailed characterisation of the spatial and phenotypic heterogeneity in tissue samples.

The impact of tumor associated macrophages on tumor biology under the lens of mathematical modelling: A review

In this article, we review the role of mathematical modelling to elucidate the impact of tumor-associated macrophages (TAMs) in tumor progression and therapy design. We first outline the biology of TAMs, and its current application in tumor therapies, and their experimental methods that provide insights into tumor cell-macrophage interactions. We then focus on the mechanistic mathematical models describing the role of macrophages as drug carriers, the impact of macrophage polarized activation on tumor growth, and the role of tumor microenvironment (TME) parameters on the tumor-macrophage interactions. This review aims to identify the synergies between biological and mathematical approaches that allow us to translate knowledge on fundamental TAMs biology in addressing current clinical challenges.

Analysis of the Equilibrium Phase in Immune-Controlled Tumors Provides Hints for Designing Better Strategies for Cancer Treatment

When it comes to improving cancer therapies, one challenge is to identify key biological parameters that prevent immune escape and maintain an equilibrium state characterized by a stable subclinical tumor mass, controlled by the immune cells. Based on a space and size structured partial differential equation model, we developed numerical methods that allow us to predict the shape of the equilibrium at low cost, without running simulations of the initial-boundary value problem. In turn, the computation of the equilibrium state allowed us to apply global sensitivity analysis methods that assess which and how parameters influence the residual tumor mass. This analysis reveals that the elimination rate of tumor cells by immune cells far exceeds the influence of the other parameters on the equilibrium size of the tumor. Moreover, combining parameters that sustain and strengthen the antitumor immune response also proves more efficient at maintaining the tumor in a long-lasting equilibrium state. Applied to the biological parameters that define each type of cancer, such numerical investigations can provide hints for the design and optimization of cancer treatments.

Mathematical investigation into the role of macrophage heterogeneity on the temporal and spatio-temporal dynamics of non-small cell lung cancers

Non Small Cell Lung Cancer (NSCLC) is the most common type of lung cancer, and represents the leading cause of cancer-related deaths worldwide. Experimental studies have shown that these solid cancers are heavily infiltrated with macrophages: anti-tumour M1 macrophages, pro-tumour M2 macrophages, and macrophage subtypes sharing both M1 and M2 properties. In this study we aim to investigate qualitatively the role of macrophages with different functional phenotypes (especially those with mixed phenotypes) on cancer dynamics and the success of different immunotherapies for cancer. To this end, we start with two time-evolving mathematical models for cancer-immune interactions that consider: (i) the effect of the two extreme phenotypes, M1 and M2 cells; (ii) the effect of M1 and M2 cells, as well as a macrophage sub-population with a mixed phenotype (throughout this theoretical study we call these cells "M12 cells"). We compare the dynamics of the two models using computational approaches, paying particular attention to the effect of different anti-cancer immunotherapies that focus on macrophages. Since data available for NSCLC and macrophage interactions are incomplete, we perform a global sensitivity analysis to see the influence of input parameters on model outcomes. Finally, we consider extensions of the previous two models to include also the spatial movement of cells, and investigate the role of macrophages with extreme phenotypes and with mixed phenotypes, on the invasion of cancer cells into the surrounding extracellular matrix (ECM). We use numerical simulations to investigate the macrophages phenotypes at the tumour center versus the invasive margin. Again, we examine the impact of immunotherapies for cancer on the spatial dynamics of cancers and immune cells, and observe a shift in the phenotype of macrophages distributed at the tumour center and invasive margin.

Mathematical investigation of innate immune responses to lung cancer: The role of macrophages with mixed phenotypes

Macrophages' role in the evolution of solid tumours is a well accepted fact, with the M1-like macrophages having an anti-tumour role and the M2-like macrophages having a pro-tumour role. Despite the fact that some clinical studies on lung tumours have emphasised also the presence of macrophages with mixed M1 and M2 phenotypes in addition to macrophages with distinct phenotypes, the majority of studies still use the distinct M1-M2 classification to predict the evolution of tumours and patient survival. In this theoretical study we use a mathematical modelling and computational approach to investigate the role of macrophages with mixed phenotype on growth/control/elimination of lung tumours. We show that tumour control in the presence of M2→M1 re-polarising treatments is mainly the result of macrophages with mixed phenotypes (due to the assumption of short half-life of M1-like macrophages). We also show that the half-life of various macrophage phenotypes (distinct M1 or mixed M1/M2 phenotypes) impacts the outcome of various therapeutic strategies targeting tumour-associated macrophages. All these results suggest the need for a better experimental understanding of the kinetics of macrophages inside solid tumours.