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

Einar Bjarki Gunnarsson,Jasmine Foo,Kevin Leder

DOI: https://doi.org/10.1016/j.jtbi.2023.111497

2023-06-14

Abstract: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.

Quantitative Methods,Populations and Evolution

What problem does this paper attempt to address?

The paper primarily aims to address the issue of estimating dynamic parameters in the processes of cancer cell proliferation and phenotypic switching. Specifically, the research focuses on the significant role played by non-genetic (epigenetic) mechanisms in cancer evolution. These mechanisms can lead to dynamic switching of cells between two or more states, and cells in different states respond differently to drug treatments. To understand how cancer evolves over time and how it responds to treatment, it is necessary to quantitatively describe key parameters such as cell proliferation rate, cell death rate, and phenotypic switching rate.

The paper proposes a statistical inference framework to estimate the aforementioned parameters using common cell line experimental data, where cells are separated into different phenotypes and expanded in culture dishes. This framework considers the stochastic dynamics of cell division, cell death, and phenotypic switching, and provides likelihood-based confidence intervals to assess the uncertainty of model parameters. The input data can be the proportion of cells in each state or the number of cells, collected at single or multiple time points.

Through theoretical analysis and numerical simulations, the study finds that:
- When using cell proportion data, only the phenotypic switching rate may be accurately estimated;
- Whereas using cell number data, not only can the net division rate of each phenotype be accurately estimated, but also the state-dependent cell division rate and cell death rate.

Finally, the authors validate the proposed framework using publicly available datasets.

In summary, the core issue that this paper attempts to solve is the development of a statistical method to accurately estimate the key dynamic parameters in the processes of cancer cell proliferation and phenotypic switching, which is crucial for understanding cancer development and guiding new treatment strategies.