Model for breast cancer diversity and spatial heterogeneity

We present and analyze a growth model of an avascular tumor that considers the basic biological principles of proliferation, motility, death and genetic mutations of the cell. From a regulatory network analysis and an analysis of genomic data we identify two sets of genes-a set of six genes and a set of sixteen genes- that are believed to play an important role in the evolution of breast cancer. Considering that cancer cells shape the tissue microenvironment and niches to their competitive advantage, the model assumes that cancer and normal cells compete for essential nutrients and that the rate of the driver mutations depends on nutrient availability. To this end, we propose a coupling between the transport of nutrients and gene mutations dynamics. Gene mutation dynamics are modeled as a Yule-Furry Markovian process, while transport of nutrients is described with a system of reaction-diffusion equations. For each representative tumor we calculate its diversity, represented by the Shannon index, and its spatial heterogeneity, measured by its fractal dimension. These quantities are important in the clinical diagnosis of tumor malignancy. A tumor malignancy diagram, obtained by plotting diversity versus fractal dimension, is calculated for different values of a parameter β, which is related to the occurrence of driver mutations. It is found that when β 1. More importantly, it is found that the results and conclusions are similar when we use the six-gene set versus sixteen-gene set.