  Predictive modeling of the evolutionary dynamics of cancer is a challenge
issue in computational cancer biology. In this paper, we propose a general
mathematical model framework for the evolutionary dynamics of cancer with
plasticity and heterogeneity in cancer cells. Cancer is a group of diseases
involving abnormal cell growth, during which abnormal regulations in stem cell
regeneration are essential for the dynamics of cancer development. In general,
the dynamics of stem cell regeneration can be simplified as a $\mathrm{G_0}$
phase cell cycle model, which lead to a delay differentiation equation. When
cell heterogeneity and plasticity are considered, we establish a
differential-integral equation based on the random transition of epigenetic
states of stem cells during cell division. The proposed model highlights cell
heterogeneity and plasticity, and connects the heterogeneity with cell-to-cell
variance in cellular behaviors, e.g. proliferation, apoptosis, and
differentiation/senescence, and can be extended to include gene
mutation-induced tumor development. Hybrid computations models are developed
based on the mathematical model framework, and are applied to the process of
inflammation-induced tumorigenesis and tumor relapse after CAR-T therapy.
Finally, we give rise to several mathematical problems related to the proposed
differential-integral equation. Answers to these problems are crucial for the
understanding of the evolutionary dynamics of cancer.
