# American Institute of Mathematical Sciences

October  2017, 14(5&6): 1379-1397. doi: 10.3934/mbe.2017071

## A mathematical model of stem cell regeneration with epigenetic state transitions

 Zhou Pei-Yuan Center for Applied Mathematics, Tsinghua University, Beijing 100084, China

* Corresponding author: Jinzhi Lei

Received  May 18, 2016 Accepted  January 15, 2017 Published  May 2017

Fund Project: This work is supported by National Natural Science Foundation of China (91430101 and 11272169).

In this paper, we study a mathematical model of stem cell regeneration with epigenetic state transitions. In the model, the heterogeneity of stem cells is considered through the epigenetic state of each cell, and each epigenetic state defines a subpopulation of stem cells. The dynamics of the subpopulations are modeled by a set of ordinary differential equations in which epigenetic state transition in cell division is given by the transition probability. We present analysis for the existence and linear stability of the equilibrium state. As an example, we apply the model to study the dynamics of state transition in breast cancer stem cells.

Citation: Qiaojun Situ, Jinzhi Lei. A mathematical model of stem cell regeneration with epigenetic state transitions. Mathematical Biosciences & Engineering, 2017, 14 (5&6) : 1379-1397. doi: 10.3934/mbe.2017071
##### References:

show all references

##### References:
Model illustration. During stem cell regeneration, cells in the resting phase either enter the proliferating phase with a rate $\beta$, or be removed from the resting pool with a rate $\gamma$. The proliferating cells undergo apoptosis with a probability $\mu$. Each daughter cell generated from mitosis is either a differentiated cell (with a probability $\kappa$) or a stem cell (with a probability $(1-\kappa)$)
Transition dynamics. (A) The cell population dynamics. (B) The percentage of epigenetic-state cells at the equilibrium state. Here, results of $\nu = 0$ (green), $20$ (blue), and $200$ (red) are shown. In simulations, the initial cell population is taken as $N(0) = 300$, and $N(0, X_i)=1, (i=1, \cdots, 300)$
Simulation of cell-state dynamics. Dynamics of cell-state proportion with different initial states (left: $(S, B, L)=(99.9, 0.05, 0.05)$, middle: $(S, B, L)=(0.05, 0.05, 99.9)$, right: $(S, B, L)=(0.05, 99.9, 0.05)$). Markers are data taken from [13]. Parameters are listed in Table 1
Parameter values in the model of cancer cell state transition. Left: the probabilities $\gamma, \kappa, \mu$ for cells of these three states. Right: the transition matrix $p(X, Y), X, Y\in \Omega$
 Parameter S B L S B L $\gamma$ 0.95 0.7 0.65 S 0.58 0.04 0.01 $\kappa$ 0.02 0.03 0 B 0.07 0.47 0 $\mu$ 0.1 0.1 0.1 L 0.35 0.49 0.09
 Parameter S B L S B L $\gamma$ 0.95 0.7 0.65 S 0.58 0.04 0.01 $\kappa$ 0.02 0.03 0 B 0.07 0.47 0 $\mu$ 0.1 0.1 0.1 L 0.35 0.49 0.09
 [1] Ebraheem O. Alzahrani, Muhammad Altaf Khan. Androgen driven evolutionary population dynamics in prostate cancer growth. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020426 [2] A. M. Elaiw, N. H. AlShamrani, A. Abdel-Aty, H. Dutta. Stability analysis of a general HIV dynamics model with multi-stages of infected cells and two routes of infection. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020441 [3] Laurent Di Menza, Virginie Joanne-Fabre. An age group model for the study of a population of trees. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020464 [4] H. M. Srivastava, H. I. Abdel-Gawad, Khaled Mohammed Saad. Oscillatory states and patterns formation in a two-cell cubic autocatalytic reaction-diffusion model subjected to the Dirichlet conditions. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020433 [5] Soniya Singh, Sumit Arora, Manil T. Mohan, Jaydev Dabas. Approximate controllability of second order impulsive systems with state-dependent delay in Banach spaces. Evolution Equations & Control Theory, 2020  doi: 10.3934/eect.2020103 [6] Reza Chaharpashlou, Abdon Atangana, Reza Saadati. On the fuzzy stability results for fractional stochastic Volterra integral equation. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020432 [7] Scipio Cuccagna, Masaya Maeda. A survey on asymptotic stability of ground states of nonlinear Schrödinger equations II. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020450 [8] Weiwei Liu, Jinliang Wang, Yuming Chen. Threshold dynamics of a delayed nonlocal reaction-diffusion cholera model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020316 [9] Manil T. Mohan. First order necessary conditions of optimality for the two dimensional tidal dynamics system. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020045 [10] Cuicui Li, Lin Zhou, Zhidong Teng, Buyu Wen. The threshold dynamics of a discrete-time echinococcosis transmission model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020339 [11] Shao-Xia Qiao, Li-Jun Du. Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, , () : -. doi: 10.3934/era.2020116 [12] Xin-Guang Yang, Lu Li, Xingjie Yan, Ling Ding. The structure and stability of pullback attractors for 3D Brinkman-Forchheimer equation with delay. Electronic Research Archive, 2020, 28 (4) : 1395-1418. doi: 10.3934/era.2020074 [13] Chao Xing, Jiaojiao Pan, Hong Luo. Stability and dynamic transition of a toxin-producing phytoplankton-zooplankton model with additional food. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020275 [14] José Luis López. A quantum approach to Keller-Segel dynamics via a dissipative nonlinear Schrödinger equation. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020376

2018 Impact Factor: 1.313