Mathematical model of adult stem cell regeneration with cross-talk between genetic and epigenetic regulation
Jinzhi Leia, Simon A. Levinb,1, and Qing Niec
aZhou Pei-Yuan Center for Applied Mathematics, Ministry of Education Key Laboratory of Bioinformatics, Tsinghua University, Beijing 100084, China; bDepartment of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ 08544; and cDepartment of Mathematics, University of California, Irvine, CA 92697
Contributed by Simon A. Levin, January 7, 2014 (sent for review September 4, 2013) Adult stem cells, which exist throughout the body, multiply by cell are continuously cycling (9, 10). Each state is likely associated with division to replenish dying cells or to promote regeneration to repair a unique microenvironment (10, 11). Dormant and homeostatic damaged tissues. To perform these functions during the lifetime of HSCs are anchored in endosteal niches through interactions with organs or tissues, stem cells need to maintain their populations in a number of adhesion molecules expressed by both HSCs and a faithful distribution of their epigenetic states, which are suscepti- niche stromal cells (10, 12). Furthermore, injury-activated HSCs ble to stochastic fluctuations during each cell division, unexpected are located near sinusoidal vessels (the perivascular niche). In injury, and potential genetic mutations that occur during many cell response to the loss of hematopoietic cells, surviving dormant divisions. However, it remains unclear how the three processes of HSCs located in their niches develop into injury-activated HSCs differentiation, proliferation, and apoptosis in regulating stem cells to undergo self-renewing divisions. In the recovery stage, injury- collectively manage these challenging tasks. Here, without consid- activated HSCs either differentiate into multipotential pro- ering molecular details, we propose a genetic optimal control model genitor cells or migrate to their osteoblastic niches to reestablish for adult stem cell regeneration that includes the three fundamental the dormant and homeostatic HSC pools (10, 13). processes, along with cell division and adaptation based on differ- The growth and regeneration of many adult stem cell pools are ential fitnesses of phenotypes. In the model, stem cells with a distri- tightly controlled with feedback regulation at different levels. For bution of epigenetic states are required to maximize expected example, HSC self-renewal and differentiation are regulated by performance after each cell division. We show that heteroge- direct HSC–niche interactions and cytokines secreted from stro- neous proliferation that depends on the epigenetic states of stem mal cells through various feedback signals (9–11). Adult intestinal cells can improve the maintenance of stem cell distributions to stem cells residing in a niche in the crypt are regulated by the create balanced populations. A control strategy during each cell paracrine secretion of growth factors and cytokines from sur- division leads to a feedback mechanism involving heterogeneous rounding mesenchymal cells (14–16). In addition, the mammalian proliferation that can accelerate regeneration with less fluctuation olfactory epithelium, a self-renewing neural tissue, is regulated in the stem cell population. When mutation is allowed, apoptosis through negative feedback signals involving the diffusive mole- evolves to maximize the performance during homeostasis after multiple cell divisions. The overall results highlight the importance cules GDF11 and activin (17). of cross-talk between genetic and epigenetic regulation and the Independent of division modes, symmetric or asymmetric cell performance objectives during homeostasis in shaping a desirable divisions may lead to daughter cells with genetic or epigenetic heterogeneous distribution of stem cells in epigenetic states. states different from the normal states. The enormous functional demands and longevity of stem cells suggest that stem cells, par- fitness function | optimization | robustness | dynamic programming | systems biology Significance
dult stem cells are present in most self-renewing tissues, in- This paper examines how adult stem cells maintain their ability Acluding skin, intestinal epithelium, and hematopoietic systems. to carry out a complex set of tasks, including tissue regeneration Stem cells provide regeneration through proliferation, differentia- and replacement of defective cells. To do so, stem cell pop- tion, and apoptosis; therefore, the accumulation of undesirable ulations must coordinate differentiation, proliferation, and epigenetic changes, which are independent of the genetic instruc- cell death (apoptosis) to maintain an appropriate distribution tions but heritable at each cell division, can lead to the causation or of epigenetic states. Using the tools of applied mathematics, progression of diseases (1, 2). Epigenetic effects such as the sto- and borrowing from the theory of intergenerational transfer chastic partitioning of the distribution of regulatory molecules of resources, this paper shows how control strategies during during cell division may change the capability of the cell to undergo cell division should be chosen to maximize expected perfor- differentiation or proliferation (3), and the accumulation of DNA mance, utilizing cross-talk between genetic and epigenetic regulation and performance criteria during homeostasis. Het- errors (or damages) can result in carcinogenesis (4–6). erogeneous proliferation, a mixed strategy in which not all Many stem cells are heterogeneous in their ability to proliferate, cells have the same proliferation probability, is shown to in- self-renew, and differentiate, and they can reversibly switch be- crease robustness, and hence long-term performance. tween different subtypes under stress conditions. Specifically, he-
matopoietic stem cells (HSCs) (see ref. 7 for a review of HSC Author contributions: J.L., S.A.L., and Q.N. designed research; J.L. and Q.N. performed heterogeneity) have distinguished subtypes (such as lymphoid de- research; J.L., S.A.L., and Q.N. contributed new reagents/analytic tools; J.L. and Q.N. an- ficient, balanced, or myeloid deficient) whose distribution depends alyzed data; and J.L., S.A.L., and Q.N. wrote the paper. on their heterogeneity during the differentiation process (7, 8). The authors declare no conflict of interest. HSCs can reversibly acquire at least three proliferative states: a See Commentary on page 3653. dormant state in which the cells are maintained in the quiescent 1To whom correspondence should be addressed. E-mail: [email protected]. stage of the cell cycle, a homeostatic state in which the cells are This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. occasionally cycling, and an injury-activated state in which the cells 1073/pnas.1324267111/-/DCSupplemental.
E880–E887 | PNAS | Published online February 5, 2014 www.pnas.org/cgi/doi/10.1073/pnas.1324267111 Downloaded by guest on September 25, 2021 ticularly the cells from highly regenerative tissues (e.g., epithelium A PNAS PLUS or blood), may be equipped with effective repair mechanisms to ensure genomic integrity over a lifetime (18). Stem cells often respond differently to genetic or epigenetic errors at different proliferation phases (19). Studies regarding the population re- sponse to DNA damage of HSCs have suggested that the system selects for the least damaged cells, and the competition between different cells is controlled by the level of p53 proteins (20, 21). Highly regenerative adult stem cells (e.g., HSCs) need to possess effective strategies that balance long-term regeneration with pro- tection from mutagenesis (for example, cell proliferation or differ- SEE COMMENTARY entiation may be affected by the DNA damage response) (20, 22). B Previous modeling studies based on the cell population dynamics have indicated that feedback regulation to the proliferation is im- portant to maintain the homeostasis of tissue growth (23–25). The exploration and analysis of models that include transit-amplifying progenitor cells and terminally differentiated cells have suggested that multiple feedback mechanisms at different lineage stages can influence the speed of tissue regeneration for better performance – (17, 26 28). These population dynamic models could include age Fig. 1. Model Illustration. (A)Atthetth cell cycle, cells in the resting phase β structure (29), evolution (27, 30), and stochasticity (30, 31); and these either enter the proliferating phase with the probability of t , or differen- models could also be applied to the regulation of cancer (32). Studies tiate into other cell types with the probability of δt . The proliferating cells μ based on spatial modeling have found that diffusive and regulatory undergo apoptosis with the probability of t . Resting phase cells occasionally molecules involved in feedback mechanisms regulating the differ- migrate to the quiescent phase and vice versa under stress. (B) The perfor- entiation capabilities of the cells are important in maintaining the mance function QðNt ,ft ðxÞÞ quantifies how well the tissue fits to its physio- logical properties. The changes in the tissue state ðNt ,ft ðxÞÞ at each cell cycle stem cell niche and shaping tissue stratification (33, 34). β μ δ are determined by the three quantities f t ðxÞ, t ðxÞ, t ðxÞg chosen to maxi- During the tissue self-renewal process driven by adult stem cells, APPLIED mize the performance at the next cycle to give QðNtþ1,ftþ1ðxÞÞ. An evolu- MATHEMATICS how do stem cells maintain a desirable distribution of epigenetic tionary fitness function at homeostasis, denoted by W, is the limit of states over their lifetimes despite many perturbations or accidental QðNt ,ft ðxÞÞ when t→∞. changes? What are the controlling strategies that enable a cell to maximize its performance at each cell cycle while contributing posi- tively to the entire cell population during tissue growth? Additionally, intrinsic cellular states that may change during cell division. are these control strategies able to guide genetic evolution to achieve Here, only epigenetic states that are significant for cell differen- high tissue performance over a long period? Without considering any tiation, proliferation, or apoptosis are considered. Consequently, molecular details, we present a dynamic programming model that the three processes have dependences on the epigenetic state x:
includes stochastic transitions between cell cycles. The model is de- SYSTEMS BIOLOGY δtðxÞ, βtðxÞ, and μtðxÞ, where the subscript t indicates the tth cell fined by the combination of a performance function at each cell di- vision and a fitness function during tissue homeostasis. We sought cycle (Fig. 1). optimal controlling strategies involving proliferation, differentiation The distribution density of stem cells during the resting phase, N and apoptosis that naturally and collectively emerge from achieving whose total population is denoted as t, with different epigenetic performance objectives as well as optimizing fitness. The model, states x, is characterized by ftðxÞ. ðNt; ftðxÞÞ undergoes a trans- which represents stem cells as a distribution of a state variable, formation from one cell cycle to the next (Fig. 1B): emphasizes the cross-talk between genetic evolution and epigenetic states and their stochastic transitions at each cell cycle. The analysis ðNt; ftðxÞÞ↦ Ntþ1; ftþ1ðxÞ : [1] and computation of the model suggest the existence of several critical R controlling strategies that regulate proliferation and apoptosis and During each cell cycle, Nt ftðxÞδtðxÞdxR cells leave the resting are highlighted by heterogeneous dependence on the epigenetic phase due to differentiation, and Nt ftðxÞβtðxÞdx cells enter states in the feedback regulation. the proliferating phase. Each cell in the proliferating phase ei- ther undergoes apoptosis with a probability of μtðxÞ or produces Results two daughter cells. Hence, the cell population after mitosis is A Model of Stem-Cell Regeneration with Epigenetic Transition. The Z Z model is based on the G0 cell cycle model (35, 36) and a dynamic Nt ¼ Nt − Nt ftðxÞδtðxÞdx − Nt ftðxÞβtðxÞdx programming approach for intergenerational resource transfer þ1 Z (37, 38) together with evolutionary dynamics (39). Stem cells at þ 2 Nt ftðxÞβtðxÞð1 − μtðxÞÞdx cell cycling are classified into resting (G0) or proliferating (G1,S, Z and G phases and mitosis) phases (Fig. 1A) (35). During each 2 N f x β x − μ x − δ x dx : cell cycle, a cell in the proliferating phase either undergoes ap- ¼ t 1 þ tð Þ½ tð Þð1 2 tð ÞÞ tð Þ optosis or divides into two daughter cells; however, a cell in the resting phase either irreversibly differentiates into a terminally The integrals are taken over all possible epigenetic states. In this differentiated cell or returns to the proliferating phase. In some derivation, the reversible transition between the resting phase and tissues, resting phase cells (e.g., HSCs) may undergo a reversible the quiescent phase is regarded as perfectly balanced for an transition to a quiescent phase with preserved self-renewal, equilibrium, which may occur during homeostasis. In this paper, we which results in two distinct cell populations. x only considered the effect of this transition for regeneration in To study the heterogeneity of cell responses, we denote as response to a severe loss of differentiated cells (SI Text, section S3). the epigenetic state of a cell, which, for example, can be the We define the observed proliferation probability as expression levels of one or multiple genes, the number of Z nucleosome modifications of a DNA region, or the positions β f x β x − μ x − δ x dx; [2] of DNA methylation. In short, x represents one or several t;obs ¼ 1 þ tð Þ½ tð Þð1 2 tð ÞÞ tð Þ
Lei et al. PNAS | Published online February 5, 2014 | E881 Downloaded by guest on September 25, 2021 then 1.05 14 A 1 Nt Ntβ : [3] 7 þ1 ¼ t; (x) ) obs 0.5 t 0 β 0 Here t;obs is the ratio of the cell population numbers between 14 B 0 100 200 300 1 two consecutive cell cycles. 7 To account for stochastic effects during the inheritance of epi- 0
genetic states that lead to variability of daughter cells in each cell C Cell population (N division (3, 40, 41), we introduced an inheritance probability pðx; yÞ, percentage of cells (%) 0.6 0.3 D which represents the probability that a daughter cell of state x comes R 0 0.95 from a mother cell of state y.Therefore, pðx; yÞdx ¼ 1foranyy. 0 100 200 300 0 1000 2000 3000 4000 Similarly to the above argument, we obtained (SI Text, section S1) x Cell cycle (t) 1 Fig. 2. Distribution of cells at homeostasis under three different combina- f x f x − δ x β x β μ tþ1ð Þ¼ tð Þð1 ð tð Þþ tð ÞÞÞ tions of the epigenetic regulation. (A) Both ðxÞ and GðxÞ are independent βt; obs of x, and δGðxÞ changes with x.(Inset) The performance function χðxÞ is Z μ β δ shown. (B) GðxÞ is independent of x, and ðxÞ and GðxÞ change with x.(C) Both δ x and β x are independent of x, and μ x changes with x. Shadow þ2 ftðyÞβtðyÞð1 − μtðyÞÞpðx; yÞdy : [4] Gð Þ ð Þ Gð Þ regions ðx < 60Þ represent defective states. (D) Time course of Nt under the three conditions (red, green, and blue for conditions A–C, respectively). (See Eqs. 3 and 4 define a transformation between two cell cycles. SI Text, section S5 for details on simulations.) During the tissue homeostasis, Eq. 3 indicates that the ob- β → t→∞ served proliferation satisfies t;obs 1as . Otherwise, either face of uncertainties in apoptosis μtðxÞ and differentiation uncontrolled growth ðβt; > 1Þ or tissue degeneration ðβt; < 1Þ obs obs δ x occurs. Hence, cell proliferation, differentiation, and apoptosis tð Þ, which leads to solving the corresponding Bellman con- dition (38, 43–45) (i.e., fβtðxÞ; μtðxÞ; δtðxÞg) must be dynamically controlled at each cell cycle, for example, through signal molecules released from Q N ; f x β x ; μ x ; δ x ; [8] EμtðxÞ;δtðxÞ max ð t tð Þj tð Þ tð Þ tð ÞÞ downstream cell lineages (17, 32). This dynamic regulation leads βtðxÞ to a limited distribution at homeostasis, where Eμ x ;δ x is the expectation with respect to apoptosis and fðxÞ¼lim ftðxÞ; [5] tð Þ tð Þ t→∞ differentiation probabilities during cell division. The evolutionary fitness function is defined as the perfor- which describes the stem cell distribution as a function of epi- mance at homeostasis after multiple cell divisions (i.e., t→∞; “ ” genetic states, and is termed tissue epigenetics for short. see also Fig. 1B): One possible control strategy for this type of growth may follow evolution akin to natural selection (42). To model this selection, we W ¼ lim QðNt; ftðxÞÞ: [9] first introduced a tissue performance function Q depending on the t→∞ population of stem cells through a function φ as well as the distri- x While the tissue performance function Q is subject to epigenetic bution of epigenetic states inthetissuethroughacellperformance W function χðxÞ, so that the performance at the tth cell cycle is given by regulation at each cell cycle, the fitness function is genetically Z regulated and dependent on the apoptosis μGðxÞ and the differ- entiation δGðxÞ. Evolution selects μGðxÞ and δGðxÞ through muta- QðNt; ftðxÞÞ ¼ φðNtÞ χðxÞftðxÞdx: [6] tions to maximize the fitness W. The overall model defines a principle of a control strategy that incorporates cross-talk be- The cell performance χðxÞ measures the capability of a cell with tween genetic and epigenetic regulation in stem cell regeneration given epigenetic state x in accomplishing its physiological func- and evolution. tions (see Fig. 2 as an example). A larger value corresponds to better performance. Heterogeneous Apoptosis Can Improve the Maintenance of Tissue We assumed that two layers of regulation occur between two cell Epigenetics. During growth, the accumulation of stochastic mod- cycles: one at the epigenetic level that occurs at each cell division, ifications in epigenetic states may produce defective cells that need and one at the genetic level that is selected by mutations over a long to be effectively repaired or removed. Here, we show that hetero- time scale of many cell divisions. In particular, the probability of geneous apoptosis is advantageous in controlling tissue epigenetics. f x t→∞ 3 proliferation βtðxÞ varies at each cell cycle by epigenetic regulation, First, the epigenetic function ð Þ, when we take in Eqs. 4 while the apoptosis probability μtðxÞ¼μGðxÞþμ^tðxÞ,inwhich and with an assumption of no epigenetic uncertainty in dif- μGðxÞ is the average probability at homeostasis and is selected ferentiation and apoptosis, satisfies the following integral equa- SI Text through genetic mutations over a long time scale and μ^tðxÞ is ran- tion ( , section S2): dom at each cell cycle due to epigenetic modulations. Similarly, the Z ^ differentiation probability takes the form of δtðxÞ¼δGðxÞþδtðxÞ f y β y − μ y p x; y dy ^ 2 ð Þ ð Þð1 Gð ÞÞ ð Þ in which δGðxÞ is the average probability at homeostasis and δtðxÞ fðxÞ¼ ; [10] represents epigenetic uncertainty. With these mechanisms of regu- δGðxÞþβðxÞ lation, the performance Q after cell division depends, through Eqs. 2–4, on the proliferation βtðxÞ as well as the stochasticities in apo- where βðxÞ¼limt→∞βtðxÞ satisfies ptosis μtðxÞ and differentiation δtðxÞ. Thus, we can write the per- Z formance function after cell division as (SI Text,sectionS1) fðxÞ βðxÞð1 − 2μGðxÞÞ − δGðxÞ dx ¼ 0: [11] Q Ntþ1; ftþ1ðxÞ ¼ QðNt; ftðxÞjβtðxÞ; μtðxÞ; δtðxÞÞ: [7] Analysis of a simplified model based on Eqs. 10 and 11 shows During each cell cycle, the proliferation βtðxÞ is controlled to that homogenous apoptosis [i.e., μGðxÞ is independent of x] easily achieve maximum tissue performance after cell division in the leads to abnormal or disease conditions for a tissue (SI Text,
E882 | www.pnas.org/cgi/doi/10.1073/pnas.1324267111 Lei et al. Downloaded by guest on September 25, 2021 β PNAS PLUS section S2). This observation is further confirmed by direct modulated such that t;2 changes at each cell cycle. Biologically, simulations of Eqs. 10 and 11 under the condition in which ap- this assumption corresponds to the situation in which, for exam- optosis probability μGðxÞ is either dependent on or independent ple, certain growth factor receptors are active (or expressed) only of x (Fig. 2). Whenever the apoptosis μGðxÞ is independent of x, in type II but not type I cells; however, the receptors are required most cells accumulate in low-performance states (Fig. 2 A and to respond to external signals to control proliferation. B μ x x β ∂Q=∂β ). In contrast, if Gð Þ is dependent on so that the cells with The probability t;2 (strategy B) is determined by t;2 ¼ 0, β 12 lower performance have a greater probability of apoptosis, only which yields an equation for t;obs similar to Eq. . In particular, a small number of low-performance cells are present during when N is close to the value Np (SI Text, section S3), one has homeostasis (Fig. 2C). These results suggest that heterogeneity A′φ N in apoptosis can improve the maintenance of acceptable tissue 1 N t ð pÞ − − β f − μ δ epigenetics during a long lifespan. N p þ N φ″ N 1 1;G t;1 2 t;1 þ t SEE COMMENTARY β ≈ t p ð pÞ : [15] Furthermore, we find that heterogeneity in the cell perfor- t;2 f − 2μt; mance function ðχðxÞÞ is important for successful natural selec- t;2 2 tion of apoptosis strategies, and epigenetic transition during cell All bar terms are averages over cell epigenetic states, with division is helpful for robust tissue epigenetics during homeo- x ∈ Ω SI a subscript 1 for type I cells ð 1Þ and a subscript 2 for type II stasis with respect to accidental changes in the tissue lifespan ( x ∈ Ω δ μ 13 Text, section S3). Interestingly, despite apparent differences in ð 2Þ. Similarly to t and t in Eq. , these average terms tissue epigenetics, homogeneous or heterogeneous apoptosis incorporate genetic and epigenetic regulation in differentiation and apoptosis. Examples of tissue dynamics based on strategy B yields similar dynamics in the cell population Nt (Fig. 2D), f x are shown in SI Text, section S3. demonstrating the importance of introducing the function tð Þ 15 β for epigenetic states into the model. The cell population model Eq. shows that t;2 is a decreasing function of the cell pop- alone may be insufficient to study the control strategies of stem ulation, resulting in a complex negative feedback regulation with respect to cell populations and the epigenetic states of the tissue cell regeneration. β cells. We note that the heterogeneous proliferation probability t;2 An Optimal Control for Proliferation During Each Cell Cycle Depends for the type II cells also depends on the probability of the type I on Complex Feedback Regulation Involving the Epigenetic States and cells, which suggests that an appropriate selection of the unmodu- lated proliferation β ;G can improve the performance at homeo- the Size of the Total Cell Population. Optimal control at each cell 1 APPLIED SI Text cycle involves identifying the proliferation probability to maxi- stasis in comparison with homogenous cells ( ,sectionS3). MATHEMATICS mize the performance Q in Eq. 8 after cell division. To study the Simple feedback via the size of the cell population (strategy C). Optimal system analytically, we considered two cases based on either ho- controls of proliferation based on our model lead to the negative mogeneous or heterogeneous proliferation. feedback regulation of proliferation through the cell population. Homogeneous proliferation (strategy A). When βtðxÞ is independent of Similar regulatory mechanisms with negative feedback have been the epigenetic state x, meaning that all cells in the tissue are alike explicitly introduced in many stem cell population models (17, in their ability to undergo cell cycle reentry, the optimal pro- 27, 36, 46) by use of a Hill function (strategy C) such as liferation (strategy A) is governed by ∂Q=∂βt ¼ 0, which yields m 1 þ ρðNt=KÞ