A mathematical model for phenotypic heterogeneity in breast cancer with implications for therapeutic strategies

Understanding and leveraging phenotypic plasticity during metastasis formation

Cancer metastasis is the process of detrimental systemic spread and the primary cause of cancer-related fatalities. Successful metastasis formation requires tumor cells to be proliferative and invasive; however, cells cannot be effective at both tasks simultaneously. Tumor cells compensate for this trade-off by changing their phenotype during metastasis formation through phenotypic plasticity. Given the changing selection pressures and competitive interactions that tumor cells face, it is poorly understood how plasticity shapes the process of metastasis formation. Here, we develop an ecology-inspired mathematical model with phenotypic plasticity and resource competition between phenotypes to address this knowledge gap. We find that phenotypically plastic tumor cell populations attain a stable phenotype equilibrium that maintains tumor cell heterogeneity. Considering treatment types inspired by chemo- and immunotherapy, we highlight that plasticity can protect tumors against interventions. Turning this strength into a weakness, we corroborate current clinical practices to use plasticity as a target for adjuvant therapy. We present a parsimonious view of tumor plasticity-driven metastasis that is quantitative and experimentally testable, and thus potentially improving the mechanistic understanding of metastasis at the cell population level, and its treatment consequences.

The role of tumor microenvironment in drug resistance: emerging technologies to unravel breast cancer heterogeneity

Breast cancer is a highly heterogeneous disease, at both inter- and intra-tumor levels, and this heterogeneity is a crucial determinant of malignant progression and response to treatments. In addition to genetic diversity and plasticity of cancer cells, the tumor microenvironment contributes to tumor heterogeneity shaping the physical and biological surroundings of the tumor. The activity of certain types of immune, endothelial or mesenchymal cells in the microenvironment can change the effectiveness of cancer therapies via a plethora of different mechanisms. Therefore, deciphering the interactions between the distinct cell types, their spatial organization and their specific contribution to tumor growth and drug sensitivity is still a major challenge. Dissecting intra-tumor heterogeneity is currently an urgent need to better define breast cancer biology and to develop therapeutic strategies targeting the microenvironment as helpful tools for combined and personalized treatment. In this review, we analyze the mechanisms by which the tumor microenvironment affects the characteristics of tumor heterogeneity that ultimately result in drug resistance, and we outline state of the art preclinical models and emerging technologies that will be instrumental in unraveling the impact of the tumor microenvironment on resistance to therapies.

Statistical inference of the rates of cell proliferation and phenotypic switching in cancer

Recent evidence suggests that nongenetic (epigenetic) mechanisms play an important role at all stages of cancer evolution. In many cancers, these mechanisms have been observed to induce dynamic switching between two or more cell states, which commonly show differential responses to drug treatments. To understand how these cancers evolve over time, and how they respond to treatment, we need to understand the state-dependent rates of cell proliferation and phenotypic switching. In this work, we propose a rigorous statistical framework for estimating these parameters, using data from commonly performed cell line experiments, where phenotypes are sorted and expanded in culture. The framework explicitly models the stochastic dynamics of cell division, cell death and phenotypic switching, and it provides likelihood-based confidence intervals for the model parameters. The input data can be either the fraction of cells or the number of cells in each state at one or more time points. Through a combination of theoretical analysis and numerical simulations, we show that when cell fraction data is used, the rates of switching may be the only parameters that can be estimated accurately. On the other hand, using cell number data enables accurate estimation of the net division rate for each phenotype, and it can even enable estimation of the state-dependent rates of cell division and cell death. We conclude by applying our framework to a publicly available dataset.

Tipping the Balance: A Criticality Perspective

Cell populations are often characterised by phenotypic heterogeneity in the form of two distinct subpopulations. We consider a model of tumour cells consisting of two subpopulations: non-cancer promoting (NCP) and cancer-promoting (CP). Under steady state conditions, the model has similarities with a well-known model of population genetics which exhibits a purely noise-induced transition from unimodality to bimodality at a critical value of the noise intensity σ2. The noise is associated with the parameter λ representing the system-environment coupling. In the case of the tumour model, λ has a natural interpretation in terms of the tissue microenvironment which has considerable influence on the phenotypic composition of the tumour. Oncogenic transformations give rise to considerable fluctuations in the parameter. We compute the λ−σ2 phase diagram in a stochastic setting, drawing analogies between bifurcations and phase transitions. In the region of bimodality, a transition from a state of balance to a state of dominance, in terms of the competing subpopulations, occurs at λ = 0. Away from this point, the NCP (CP) subpopulation becomes dominant as λ changes towards positive (negative) values. The variance of the steady state probability density function as well as two entropic measures provide characteristic signatures at the transition point.