Open Access

Stem cell biology is population biology: differentiation of hematopoietic multipotent progenitors to common lymphoid and myeloid progenitors

Theoretical Biology and Medical Modelling201310:5

DOI: 10.1186/1742-4682-10-5

Received: 27 August 2012

Accepted: 22 December 2012

Published: 17 January 2013

Abstract

Abstract

The hematopoietic stem cell (HSC) system is a demand control system, with the demand coming from the organism, since the products of the common myeloid and lymphoid progenitor (CMP, CLP respectively) cells are essential for activity and defense against disease. We show how ideas from population biology (combining population dynamics and evolutionary considerations) can illuminate the feedback control of the HSC system by the fully differentiated products, which has recently been verified experimentally. We develop models for the penultimate differentiation of HSC Multipotent Progenitors (MPPs) into CLP and CMP and introduce two concepts from population biology into stem cell biology. The first concept is the Multipotent Progenitor Commitment Response (MPCR) which is the probability that a multipotent progenitor cell follows a CLP route rather than a CMP route. The second concept is the link between the MPCR and a measure of Darwinian fitness associated with organismal performance and the levels of differentiated lymphoid and myeloid cells. We show that many MPCRs are consistent with homeostasis, but that they will lead to different dynamics of cells and signals following a wound or injury and thus have different consequences for Darwinian fitness. We show how coupling considerations of life history to dynamics of the HSC system and its products allows one to compute the selective pressures on cellular processes. We discuss ways that this framework can be used and extended.

Keywords

Hematopoieitic stem cell Multipotent progenitor Common lymphoid progenitor Common myeloid progenitor Darwinian fitness Natural selection Population dynamics

Introduction

Hematopoiesis (the formation of blood components) is a highly orchestrated and dynamical process. Hematopoietic Stem Cells (HSCs) give rise, through a large array of successively differentiated progeny, to mature blood cells. While progress has been made in understanding the HSC system, particularly at the molecular level [1, 2], Tan et al. (pg 82-83) [3] concluded a recent review on HSCs by identifying critical unanswered questions: 1) Is it possible to manipulate adult stem cells to increase their ability to proliferate in vitro while maintaining stem cell qualities, so that adult stem cells can be used as a sufficient source of tissue for transplants and other therapeutic strategies? 2) What are the intrinsic and extrinsic controls that keep stem cells from differentiating or that direct them along a particular differentiation pathway to form one specialized cell type rather than another? 3) What are the factors responsible for stem cell responses to injury or damage that enable rapid activity and appropriate contribution to tissue repair and regeneration? 4) Can the ‘stress’ signals’ that command facultative stem cells to respond to tissue damage and gain specific regenerative quality be harnessed for therapeutic value?

We will show that these questions can only be fully answered if one considers the connection between the needs of the organism and the HSC system, a demand control system [4] in which the demand comes from the organism via signals to the HSC system. As noted by Metcalf [4] in his classic lectures, appropriate long-term behavior of the HSC system is essential for the health of the organism. Furthermore, HSCs are required to produce the required differentiated cells without depleting the stem cell pool or creating damaged stem cells that lead to cancer. Even in the absence of oxidative damage due to ischemia/reperfusion, blood loss, or infections damage to blood cells and tissues caused by reactive oxygen species is a common process leading to cell apoptosis, and mature blood cells have a steady rate of turnover and are thus constantly replaced [5].

Thus, organismal performance is intimately connected to HSCs and their products [69]. The myeloid products of the HSC system are particularly important for antipredator and foraging behaviors and general immune response and the lymphoid products of the HSC system are essential for specific immune response and neonatal survival. Recent whole organism manipulations have demonstrated the veracity of this proposition (see [10] for myeloid cells; [11] for lymphoid cells). At the same time, when the organism is in a steady state (homeostasis) the HSC system and its products are relatively stable [12]. In addition, the demand is not unlimited; for example, it has been known for a very long time that organismal performance is a peaked function of hematocrit (e.g., [6, 9, 13]).

In this paper, we show that the questions raised by Tan et al. [3] are answered most effectively if one takes an approach to stem cell biology based on population biology, which involves combining the dynamics of cells with considerations of fitness and natural selection [14]. We develop a theoretical framework that complements existing models (Additional file 1: Table S1) to explore the dynamics of the HSC system. Elsewhere ([15]; Figure 1 here) we have developed stochastic and deterministic models of HSC cells and their associated products and applied evolutionary invasion analysis [16] and state dependent life history theory [1720] to show that understanding the dynamics of HSCs and their products requires asking more than whether a stem cell renews, symmetrically differentiates, or asymmetrically differentiates. Understanding the roles of positive and negative feedback is essential for predicting stem cell dynamics. By linking these feedback processes to stochastic population models (which allow uncertainties inherent in the system to be accounted for) we showed how well the overall mean dynamics of the system can be approximated a system of ordinary differential equations. We build on this work here and introduce additional concepts associated with population biology to the biology of stem cells, with the focus on the HSC system. We introduce the Multipotent Progenitor Commitment Response (MPCR) that characterizes the penultimate differentiation of a multipotent progenitor (MPP) to a Common Lymphoid Progenitor (CLP) or a Common Myeloid Progenitor (CMP), i.e. whether they follow a myeloid or lymphoid track. Although we recognize that within the myeloid track there is another decision towards a granulocyte-macrophage progenitor or megakaryocyte-erythrocyte progenitor (see Additional file 1). We show how the fitness (survival and reproduction) of the organism shapes the MPCR, thus providing an approach for modeling the demand control nature of the HSC system.
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Fig1_HTML.jpg
Figure 1

A diagrammatic derivation of Eqns 1 to 6 (details given in Additional file 1 ). a) In the most general case, we consider stem cells (S), a series of Multipotent Progenitor Cells (MPP), a Common Lymphoid Progenitor (CLP) and a Common Myeloid Progenitor (CMP). CLPs give rise to B, NK, and T cells; CMPs give rise to Erythrocytes (E), Granulocytes (G), and Platelets (P). We denote the total numbers of lymphoid and myeloid cells by L and M respectively, rates of differentiation by r ·(with subscript indicating the cell type involved), rates of development of MPP cells by λ ·, feedback from fully differentiated cells on those rates by Φ ·, and rates of cell death by μ ·. The feedback functions have the property that they are 1 when stem cell or fully differentiated cell numbers are low and decline as stem cells or fully differentiated cells increase. Thus, for example, stem cells renew (one stem cell becomes two) at rate r s Φ s (l,m)when the concentrations of lymphoid and myeloid cells are l and m respectively, asymmetrically differentiate (one stem cell becomes two stage-0 progenitors) at rate 2 r p Φ p ( l , m ) Φ s ( l , m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq1_HTML.gif, symmetrically differentiate (one stem cell becomes a stem cell and a stage-0 progenitor) at rate r p Φ s (l,m), and die at rate μ s . Similar interpretations hold for other transitions. The Multipotent Commitment Response (MPCR), denoted by ρ(l,m), is the probability that a MPP in its final stage commits to the lymphoid route. b) To focus on the MPCR, we combine all of the fully differentiated cells into lymphoid and myeloid classes (L and M) and use Michaelis-Menten-like arguments to compress the MPP class into a single stage, assuming that steady states of intermediate stages are rapidly reached, characterized by combination of rate constants Ω N .

Methods

We begin first by describing, in summary here, with details in Additional file 1, the dynamics of the stem cells and their descendants, after which we describe the components of fitness (survival and the reproduction) and their dynamics. We then couple the two together.

Dynamics of stem cells and their descendants

In the Additional file 1, we derive the computational model given below and in Figure 1 we provide a graphical representation of the full mathematical model. We thus consider stem cells, with concentration denoted by [S], Multipotent Progenitor cells, denoted by [MPP], Common Lymphoid Progenitors, denoted by [CLP], Common Myeloid Progenitors, denoted by [CMP], fully differentiated lymphoid cells, denoted by [L](measured in numbers per milliliter), and fully differentiated myeloid cells, denoted by [M](measured in numbers per nanoliter). Central to these dynamics we assume that the stem cell niche can support at most K stem cells and that in absence of all other feedback (described below), the dynamics in the niche follow Gompertzian kinetics (justified in [15]).

The dynamics are described by the following set of coupled ordinary differential equation:
d [ S ] dt = [ S ] · log ( K / [ S ] ) ( r s r p Φ p ( [ L ] , [ M ] ) ) × Φ s ( [ L ] , [ M ] ) μ s [ S ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ1_HTML.gif
(1)
d [ MPP ] dt = [ S ] · log ( K / [ S ] ) ( r s + 2 r p Φ p ( [ L ] , [ M ] ) ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ2_HTML.gif
(2)
× Φ s ( [ L ] , [ M ] ) + ( λ r d , MPP ) Φ p ( [ L ] , [ M ] ) × [ MPP ] μ p [ MPP ] d [ CLP ] dt = r d , MPP Φ p ( L , M ) ρ ( [ L ] , [ M ] ) Ω N [ MPP ] r CLP [ CLP ] μ CLP [ CLP ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ3_HTML.gif
(3)
d [ CMP ] dt = r d , MPP Φ p ( L , M ) ( 1 ρ ( [ L ] , [ M ) ] ) Ω N [ MPP ] r CMP [ CMP ] μ CMP [ CMP ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ4_HTML.gif
(4)
d [ L ] dt = r CLP [ CLP ] + ( r l μ l μ l I V ( t ) > v th ) [ L ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ5_HTML.gif
(5)
d [ M ] dt = r CMP [ CMP ] + ( r m μ m ) [ M ] . https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ6_HTML.gif
(6)
where Φ s ([L],[M]), Φ p ( [ L ] , [ M ] ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq2_HTML.gif and Φ p ([L],[M]) are feedback functions (see below) for the activity of stem cells (s), the asymmetric differentiation of stem cells (p ), and the activity of MPP cells (p), respectively. ρ([L],[M])is the demand control function (see below) that describes the probability of an MPP cell differentiating into a lymphoid or myeloid progenitor. Ω N [MPP] is the survival of MPP cells from initial differentiation through to ultimate differentiation into CLP or CMP cells (see Additional file 1 for its full derivation). μ lcharacterizes the additional mortality when the immune system is activated and I a>b is an indicator function that is 1 if a>band 0 otherwise. All other parameters are defined in Table 1.
Table 1

Variables, parameters, their interpretation, and values

Symbol

Interpretation

Value

t

Non-dimensional time

1-3500

[S]

Concentration of stem cells at time t

Eqn 1

[MPP]

Concentration of Multipotent Progenitor (MPP) cells at time t

Eqn 2

[CLP]

Concentration of Common Lymphoid Progenitor (CLP) cells at time t

Eqn 3

[CMP]

Concentration of Common Myeloid Progenitor (CMP) cells at time t

Eqn 4

[L]

Concentration of fully differentiated Lymphoid (L) cells at time t

Eqn 5

[M]

Concentration of fully differentiated Myeloid (M) cells at time t

Eqn 6

K

Maximum number of stem cells in a niche

10

r s

Maximum rate of stem cell self-renewal

2.5

r p https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq3_HTML.gif

Maximum rate of stem cell asymmetrical division

0.001

Φ p ( [ L ] , [ M ] ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq4_HTML.gif

Feedback control from fully differentiated cells to asymmetric division

Eqn 9

Φ s ([L],[M])

Feedback control from fully differentiated cells to stem cell self-reneval

Eqn 7

Φ p ([L],[M])

Feedback control from fully differentiated cells to symmetric division

Eqn 8

μ s

Rate of stem cell death

0.004

λ

Rate of MPP multiplication

0.25

μ p

Rate of MPP cell death

0.02

Ω N

Combination of intermediate multipotent progenitor rate constants

1.0

r CLP

Rate of division of CLP into fully differentiated lymphoid cells

0.01

μ CLP

Rate of CLP cell death

0.001

r l

Rate of multiplication of lymphoid cells

0.025

μ l

Rate of lymphoid cell death when immune system is not activated

0.028

μ l

Additional rate of lymphoid cell death when immune system is activated

0.01

I a>b

Indicator function for the inequality

=1 if a>b, 0 otherwise

v th

Threshold concentration for pathogens to activate the immune system

0.025

r CMP

Rate of division of CMP into fully differentiated myeloid cells

0.01

μ CMP

Rate of CMP cell death

0.001

r m

Rate of multiplication of myeloid cells

0.0

μ m

Rate of myeloid cell death

0.01

l

Value of [L]

varies

m

Value of [M]

varies

ϕ sl (l)

Feedback control of fully differentiated lymphoid cells on stem cell activity

Eqn 10

ϕ sm (m)

Feedback control of fully differentiated myeloid cells on stem cell activity

Similar to Eqn 10

ϕ pl (l)

Feedback control of fully differentiated lymphoid cells on symmetric renewal

Eqn 11

ϕ pm (m)

Feedback control of fully differentiated myeloid cells on symmetric renewal

Similar to Eqn 10

ϕ p l ( l ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq5_HTML.gif

Feedback control of fully differentiated lymphoid cells on asymmetric renewal

Eqn 11

ϕ p m ( m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq6_HTML.gif

Feedback control of fully differentiated myeloid cells on asymmetric renewal

Similar to Eqn 10

α sl

Feedback parameter in ϕ sl (l)

10

α pl

Feedback parameter in ϕ pl (l)

100

α p l https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq7_HTML.gif

Feedback parameter in ϕ p l ( l ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq8_HTML.gif

20

α sm

Feedback parameter in ϕ sm (m)

0.1

α pm

Feedback parameter in ϕ pm (m)

0.001

α p m https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq9_HTML.gif

Feedback parameter in ϕ p m ( m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq10_HTML.gif

0.2

α

Coefficient in MPP Commitment Response (MPCR)

Varies

γ

Exponent in MPCR

Varies

κ l

Density of lymphoid cells in homeostasis

30

κ m

Density of myeloid cells in homeostasis

30,000

ρ h

Fraction of lymphoid cells in homeostasis

Eqn 15

Δf([M]

Rate of accumulation of fitness when myeloid cell concentration is [M]

Eqn 17

F ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq11_HTML.gif

Fitness accumulated to time t

Eqn 18

ε

Ratio of organismal to cellular time scale

0.05

S(t)

Survival to time t

Eqn 21

μ e ([M])

Total rate of mortality when myeloid cell concentration is [M]

Eqn 20

μ e0

Myeloid independent rate of mortality

0.05

μ e1

Myeloid dependent rate of mortality

5.0

μ i ([V])

Additional rate of mortality when concentration of infectious agents is [V]

Eqn 20

μ i0

Coefficient of [V] in additional mortality

0.02

μ i1

Coefficient of [V]2 in additional mortality

0.002

r v

Replication rate of infectious agents

0.05

c l

Clearance rate of infectious agents by lymphoid cells

0.05

c m

Clearance rate of infectious agents by myeloid cells

0

v 0

Concentration of infectious agents at the start of an infection

1

v th

Concentration of infectious agents below which additional lymphoid mortality does not occur

0.025

(Parameter values are a canoncial fixed set, arbitrarily chosen, to illustrate the general principles of an MPCR).

Feedback control, which requires nonlinear dynamics, is essential for the growth and regeneration of tissues. A recent model of genetic products [21] to characterize the erythroid-myeloid lineage decision, shows how nonlinearities arise. de Graaf et al. [22] showed that the platelet concentration can regulate HSCs, through the concentration of Thrombopoietin (TPO). Other authors have shown that Lkb1, a kinase enzyme best known as a tumor suppressor, provides feedback control of HSCs [2326]. Such work provides the empirical basis for our modeling.

To incorporate feedbacks we extend Lander et al.’s [27] approach: if r denotes a generic reaction rate constant and [χ] the concentration of fully differentiated cells, they model the reaction rate with feedback control is r 1 + ν [ χ ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq12_HTML.gif. We have adapted this framework for the HSC system, with changes. First, two kinds of differentiated cells provide feedback; we use l and m to denote the total concentration of lymphoid and myeloid cells. Second, there is potentially different feedback on the activity of stem cells, the asymmetric differentiation of stem cells, and the activity of MPP cells. Third, the system needs to be active when there is a shortage of either lymphoid or myeloid cells. Thus we set:
Φ s ( l , m ) = max [ ϕ sl ( l ) , ϕ sm ( m ) ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ7_HTML.gif
(7)
Φ p ( l , m ) = max [ ϕ pl ( l ) , ϕ pm ( m ) ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ8_HTML.gif
(8)
Φ p ( l , m ) = max [ ϕ p l ( l ) , ϕ p m ( m ) ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ9_HTML.gif
(9)
where ϕ sl (0) = ϕ sm (0) = 1(representing the feedback of lymphoid and myeloid cells on stem cell activity), etc, and all ϕ ij are decreasing functions of their arguments, as in
ϕ sl ( l ) = 1 1 + α sl · l https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ10_HTML.gif
(10)
ϕ pl ( l ) = 1 1 + α pl · l https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ11_HTML.gif
(11)
ϕ p l ( l ) = 1 1 + α p l · l https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ12_HTML.gif
(12)

where α sl , α pl , α p l https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq13_HTML.gif

are parameters. A similar form is used for the feedback control from myeloid cells. This is the simplest form of feedback between the whole organism and the bone marrow stem cell system; see [28] for alternatives.

In a demand control system the probability of a MPP cell ultimately differentiating into a CLP or CMP cell must depend upon the state of the organism. That is, the current densities of myeloid and lymphoid cells determine the appropriate response. We choose a functional form that is widely used in population biology and similar to Michaelis-Menten enzyme kinetics;
ρ ( l , m ) = α m l γ 1 + α m l γ https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ13_HTML.gif
(13)
where ρ(l,m) represents the demand control function and αand γare parameters that describe this specific asymptotic Michaelis-Menten function. We let κ m :κ l denote the ratio of myeloid to lymphoid cells in homoeostasis. If ρ h denotes the value of ρ(l,m)in homeostasis then on average we have:
ρ h = κ l κ l + κ m https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ14_HTML.gif
(14)
and from Eqn 13 we have
ρ h = α ( κ m / κ l ) γ 1 + α ( κ m / κ l ) γ . https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ15_HTML.gif
(15)
We view the unknowns in this equation as the two parameters α and γ from which we find
α = ρ h 1 ρ h ( κ m / κ l ) γ . https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ16_HTML.gif
(16)
Eqn 16 determines a curve in the γ-αplane and every value on this curve provides the same value of ρ(l,m)=ρ h . However, when out of homeostasis, the value of α(γ)that determines Eqn 16 has a profound effect on the fate of HSC descendants. In Figure 2a we show the relationship between the parameters αand γ of the MPP commitment response when homeostasis corresponds to 1 lymphoid cell per 1000 myeloid cells, a typical ratio for humans. Each point on the lines in these panels correspond to a particular value of the pair (γ,α)consistent with the number of cells at homeostasis. In Figure 2b, we show how ρ(l,m) varies as m/l varies for three values of γ. As γincreases, the MPP commitment response becomes more sensitive to variation in the densities of myeloid and lymphoid cells. These three curves represent just a few possibilities in the infinite space of functional responses corresponding to Eqn 13 and raises the question: how can we predict which response an organism will use (i.e., where on the curve in Figure 2a will a population of organisms sit)? We discuss this below.
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Fig2_HTML.jpg
Figure 2

a) The relationship between the parameters α and γ of the stem cell commitment response when homeostasis corresponds to 1 lymphoid cell per 1000 myeloid cells. b) Different values of γaffect how the MPCR varies with changes in the number of lymphoid and myeloid cells. In the presence of high numbers of myeloid cells, the demand response is to drive the MPPs to make more lymphoid cells.

How an organism goes out of homeostasis depends upon its environment. For example, in an environment when wounds occur frequently, we anticipate the m/lwill be lower than the value in homeostasis while in an environment when infection occurs frequently we anticipate that m/l will be greater than that value in homeostasis.

In this paper, we are interested in ρ(l,m), and in particular how natural selection affects it, in light of the environment of the organism. To do this, we need to couple Eqns 1-13 to organismal fitness, which is what separates our work from all that has come before it. We consider it advisable to step back from modeling a particular situation, but instead consider the more general properties that connect the needs of the organism with the activity of the bone marrow stem cell system.

The components of fitness and their dynamics

The representation of genes in subsequent generations is determined by survival and successful reproduction of the focal organism. Regarding the latter, we assume that the rate at which successful reproduction occurs (Δf([M])) is a function of myeloid cells, justified by the long recognition that organismal performance is a peaked function of hematocrit. Following figure two in [9], we set
Δf ( [ M ] ) = a 0 + a 1 · M + a 2 · M 2 https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ17_HTML.gif
(17)
provided this expression is positive; otherwise we set Δf([M])=0, where a 0=−1.034565, a 1=0.001527, and a 2=−0.0000002864. The peak of Δf([M]) occurs at [M]=2666(Figure 3). Hematocrit is the fraction of myeloid cells in blood; if too great then fitness is impaired as an organism has few lymphoid (immune) cells. If too low then there are insufficient myeloid cells to support the oxygen needs of the organism. Physiologically, the parabolic hematocrit function emerges through the relationship between the oxygen carrying capacity of blood and blood flow. When myeloid cells are low (M→0), blood flow is maximal but oxygen concentration is low. In contrast as M, oxygen concentration is maximal but blood is too viscous to flow. Hence, this proximate physiological explanation and consequently the ultimate evolutionary fitnees cost lead to a peaked (parabolic) hematocrit function. The numerical relationship between reproductive rate and myeloid cell concentration described in equation (17) captures this parabolic shape.
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Fig3_HTML.jpg
Figure 3

We assume that the rate at which successful reproduction accumulates, Δ f ( m ) is a parabolic function of the density of myeloid cells m .

We let F ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq14_HTML.gif denote lifetime fitness accumulated to time t, S ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq15_HTML.gif denote survival to time t, and ε<<1a scaling parameter that relates the organismal and cellular time scales. Then
d F dt = ε · S ( t ) · Δf ( [ M ( t ) ] ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ18_HTML.gif
(18)
To determine survival, we assume that uninfected individuals have a per unit time rate of mortality with myeloid-independent and myeloid-dependent components so that the total rate of mortality is (Figure 4)
μ e ( [ M ] ) = μ e 0 + μ e 1 [ M ] https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ19_HTML.gif
(19)
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Fig4_HTML.jpg
Figure 4

We assume that the rate of mortality declines with increasing numbers of myeloid cells, which has the effect that annual survival increases with increasing densities of myeloid cells; here we artificially hold the myeloid cells constant.

Although we focus on non-fatal diseases here, such diseases can still increase mortality rate, e.g. by reducing the effectiveness of flight responses. Hence we assume that the additional mortality induced by the pathogen is
μ i ( [ V ] ) = μ i 0 [ V ] + μ i 1 [ V ] 2 https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ20_HTML.gif
(20)
so that
d S dt = ε · ( μ e ( [ M ] ) + μ i [ V ] ) S https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ21_HTML.gif
(21)
Finally, we incorporate the dynamics of infectious agents. We assume that in the absence of immune response the growth of the infectious agent is exponential with rate r v and that lymphoid and myeloid cells clear the infection at rate c l and c m respectively. In this model, again for simplicity, we ignore memory in the immune system. Thus, if V(t) denotes the density of pathogens
dV dt = [ r v c l L ( t ) c m M ( t ) ] · V ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ22_HTML.gif
(22)

Integration of Eqns 17-22 forward in time, conditioned on ρ([L],[M]) and a sequence of wound and/or infection events, allows us to compute the lifetime reproductive success associated with that particular MPP commitment response.

Eqns 1-22 are a set of deterministic ordinary differential equations that link the behavior of the stem cell system with the needs of the organism. However, organisms in nature experience wounding and infection in a quasi-random manner. We account for this in the following way. Imagine that there are K w and K i times at which wounds or infections can occur (these values could, of course, be random variables but we treat them as fixed in this paper, only for purposes of simplicity) and then determine a sequence of times T w (k w ),k w =1,2,....,K w and T i (k i ),k i =1,2,....,K i at which either a wound or infection occurs (in principle both could occur at one time). To illustrate the ideas, we assume that when a wound occurs, myeloid cells drop by 40% and that when an infection occurs, the infectious agent increases to the level v 0. These occur instantaneously and we then continue with the solution of the differential equations. For the results shown here, we assume that K w =K i =7 and that the times are uniformly distributed over the interval between day 0 and day 1500.

Results

As introduced above (Figure 2a) a homeostatic ρ h ratio of fully differentiated myeloid to lymphoid cells specifies a curve in the γα plane in which all points on this curve are consistent with ρ h but as illustrated in Figure 2b, different values of γ(and thus α) will lead to different MPP commitment responses when the organism is out of homeostasis. Thus, to understand how natural selection will act on ρ(l,m)requires dealing explicitly with components of fitness and linking them to the dynamics of the HSC system.

To explore this hypothesis, we assume that the rate at which the organism accumulates fitness, Δf(m) is a parabolic function of myeloid cells (Figure 3) and that the rate of mortality is a declining function of myeloid cells (Figure 4) and an increasing function of the density of infectious agents (Eqn 22).

In a ‘deterministic’ or laboratory environment with neither wounding nor infection, organisms still die, so that survival declines with age (Figure 5a) and fitness accumulates but ultimately saturates because of the declining survival (Figure 5b). The strength of selection on ρ(l,m) will depend upon how fitness varies with γ. We show this in Figure 6 fitness as a function of γfor the deterministic case, wounding only, infections only, and wounding and infections. Clearly there is little selection on γ in the laboratory case or of wounding only. Selection does occur when there are infections, with larger values of γ, leading to a more responsive MPCR.
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Fig5_HTML.jpg
Figure 5

Even in a laboratory environment, without wounding or infection, organism do not live forever, so that survival declines with age (panel a) with the consequence that accumulated fitness saturates.

https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Fig6_HTML.jpg
Figure 6

Plotting lifetime accumulated fitness as a function of γ allows us to understand the strength of selection on γ as determined by the environment in which the organism lives.

The alternative to varying γand holding the environment at one stochastic realization is to hold γ constant and consider multiple realizations of the stochastic environment. We show the results of such an approach in Figure 7, for γ=2. Using these, we can compute both a mean and variance for the selection acting on γ.
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Fig7_HTML.jpg
Figure 7

Ten realizations of the model with both wounds and infection, for the case of γ = 2.

Discussion

Here we have developed and analysed a theoretical framework for linking the population biology of the hematopoietic stem cells to the demands of the individual. We have introduced the notion of the MPCR (Multipotent Progenitor Commitment Response) as the response that the describes the penultimate decision of stem cells before commitment to either a myeloid or lymphoid lineage. We use this response to investigate the control dynamics of a hematopoietic stem cell system and show that different values of the ‘shape’ parameters that describe the MPCR give a range of optimal response. Below we discuss the implications of this on the evolutionary dynamics of the HSC system and, more broadly, for developing a theory of stem cell systems based in population biology.

The meaning of a flat fitness surface

The first derivative, F ∂γ https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq16_HTML.gif, of the accumulated fitness with respect to γis a fundamental measure of the strength of natural selection on γ. As the derivative becomes smaller, the strength of natural selection becomes weaker, with the implication that a wider spectrum of values of γwill provide equal values of fitness. But, as we have shown in Figure 2b, differing values of γ will lead to considerably different MPCRs, and thus kinetics of the HSC descendants following an external challenge such as a transplant or perturbation, with the prediction that if one uses animals with little evolutionary history of wounding or infection, a wide range of HSC dynamical responses is expected. For instance, among 44 laboratory mice, Abkowitz et al. [29] observed seven different patterns of donor cell dynamics following hematopoeitic stem cell transplant experiments, suggesting that there is individual heterogeneity in the parameters of the MPCR, as we would predict.

In previous work [15] we showed that the differential equations used here are a good approximation for the mean of underlying stochastic system. Understanding the limitations imposed by stochastic fluctuations on the feedback in our model [30] is an important next step because the comparisons of models and data will require a framework to account for process stochasticity and observation error. Developing appropriate state-space approaches to understand the consequences of these heterogeneities and nonlinearities on stem cell dynamics is clear and obvious future step.

Recently, Huang [31] proposed that a systems biology of stem cells should consider not only the link between observation and process through state-space approaches (e.g. [32]) but also the consequences of nonlinearity and non-genetic heterogeneity in cell systems. Here, to start to understand how nonlinearity effects influence stem cell behavior we have introduced the notion of a MPCR under the control from the needs of the organism. Furthermore, under a paradigm of stem cell heterogeneity, we should anticipate variation at multiple levels within any stem cell system; particularly so with individual variation in the MPCR (γ and α).

How this model can be extended

There are a number of extensions that go beyond the current work for linking population biology and stem cell systems. For instance, experimental validation of the MPCR would require repeated cell count measures of long-term HSCs, short-term HSCs and a range of HSC-derived products. Given such experimental data we envisage that it would be possible to assess goodness of fit between an MPCR model and data (and also characterize unexplained heterogeneties) using computational and statistical methods (Bonsall and Mangel, unpublished).

Furthermore, our framework could be extended to study the consequences of transplant or perturbation effects A perturbation experiment can be modeled by starting the HSC system in its steady state and then reducing the number of lymphoid or myeloid cells and then integrating Eqns 1-22 forward. In addition to predicting the kinetics of fully differentiated cells, we can predict the activity of the stem cells following the perturbation. As described in the Additional file 1, a transplant experiment involves the addition of the dynamics of host cells, which are dying, and allows us to predict the fraction of fully developed lymphoid or myeloid cells. Our approach can also be used, in conjunction with evolutionary invasion analysis as in [15], to predict the outcome of a partial transplant of stem cells, in which stem cells with different kinetic properties are transplanted into an organism with a healthy but not vibrant stem cell system and for which we wish to predict whether the transplanted stem cells will co-exist with the host stem cells, overtake them, or be extinguished by them.

Computing ρ(l,m)from more fundamental principles: the fitness control hypothesis

We began with the assumption that ρ ( l , m ) = α m l γ 1 + α m l γ https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq17_HTML.gif, but there are an infinite number of functions that would have the properties we seek: bounded by 1 and increasing as the the ratio m/lincreases. Given experimental data, it is possible to use Bayesian methods to determine ρ(l,m) (Bonsall and Mangel, unpublished), but we may also ask if ρ(l,m)could be constructed from some more fundamental principles. One way to do approach this woul be through state dependent life history, as implemented by stochastic dynamic programming (SDP) [1720], which we explain now but leave the details for a subsequent paper.

Our results point to organismal fitness being a function of fully differentiated lymphoid and myeloid cells, F ( l , m , t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq18_HTML.gif, which we have modeled (Eqns 17-22). Beginning there, consider the question from the perspective of the organism, rather than from the perspective of the dynamics of cells. Suppose that s denotes the time scale for the organism, so that in the interval s to s + 1 the time Δt>>1elapses in the cellular model. During that elapsed time, we assume that up to P https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq19_HTML.gif fully differentiated cells are produced and the challenge for the organism is to allocate these fully differentiated cells to lymphoid or myeloid cells. That is, if δl and δm denote the number of fully differentiated cells produced during Δt, subject to δl + δm P https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq20_HTML.gif, then we seek the value of P https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq21_HTML.gif and combination of δ l and δ m that maximizes organismal fitness.

To do this, we let F(l,m,s)denote the accumulated lifetime reproduction from time s until time S given that [L(s)]=l,[M(s)]=m and the organism is not infected. Similarly, we let F i (l,m,v,s)denote accumulated lifetime reproduction under the same conditions about lymphoid and myeloid cells and that the density of infectious agents is [V(s)]=v. Since fitness cannot be accumulated after time S, we have the end condition F(l,m,S)=F i (l,m,v,S)=0. Methods of SDP, along with value iteration, allow us to compute the stationary values of these functions for s<<S[15].

Such a method allows us to represent the optimal production of CLP and CMP cells given that the organism is not infected at time s and that [L(s)]=land M[s]=m. Then
P ( l , m , s ) = δ 0 l ( l , m , s ) + δ 0 m ( l , m , s ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ23_HTML.gif
(23)

is the optimal production of MPP cells at time s, given the current state of the organism. Similar arguments generate the optimal production of MPP cells, distributed as CLP or CMP for the case of the organism being infected at time s.

In the stationary state, fitness is only a function of l,mand v; we denote these values by F ¯ ( l , m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq22_HTML.gif and F ¯ i ( l , m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq23_HTML.gif. Similarly, the analog of the quantities in Eqn 23 are δ ¯ 0 l ( l , m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq24_HTML.gif, δ ¯ 0 m ( l , m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq25_HTML.gif and P ¯ ( l , m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_IEq26_HTML.gif. Thus, for an organism that is currently not infected, the fraction of CLP cells produced is
ρ ( l , m ) = δ ¯ 0 l ( l , m ) δ ¯ 0 l ( l , m ) + δ ¯ 0 m ( l , m ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-10-5/MediaObjects/12976_2012_Article_425_Equ24_HTML.gif
(24)

This equation provides an explicit form for ρ(l,m)based directly on fitness of the organism, with a similar one holding for the situation in which the organism is infected.

Connection to empirical studies

Our work complements the rapid recent development of understanding how gene products that regulate HSCs operate [33]. However, translating that understanding into a functional perspective at the organismal level is less well developed. Recent work suggests that epigenetic factors [34] are important in controlling HSC behavior [35], but the scale of those factors is unknown, and evidence suggests that HSC behavior is regulated by circadian oscillations modulated by photic cues [36], as is the behavior of whole organisms [37]. Other the hand, it is less clear how the differentiated T and B cells are regulated in adaptive immunity [38]. DNA damage increases in HSCs with age and such damage is more poorly repaired in older individuals [39, 40]; lymphocytes also age [41]. Rossi et al. [42] argue that the reduction in adaptive immune capacity with age and the increased incidence of myeloproliferative diseases both have origins in changes in the HSC system, presumably due to increased damage of HSCs that could include telomere shortening or other genome maintenance failures, and other damage to DNA or proteins [4346].

Thus, currently much is known about the mechanism of HSC system and its descendants [47, 48] but much less is known about how natural selection has shaped these mechanisms, even though stem cells are units of natural selection [49]. Using classic ecological methods of allometric analysis (e.g., [50]) and compartmental modeling, it has recently been argued [5153] that the fundamental architecture of the HSC system is unchanged across mammalian species (also see [54]), which suggests a long evolutionary history determining the role of HSCs in organismal ecology. Recognizing that evolutionary considerations are essential for understanding every problem in biology [55], our goal here has been to understand how the penultimate route of differentiation towards a common myeloid or lymphoid progenitor cell [5658] can be characterized both dynamically and evolutionarily.

One common way for approaching an understanding of the dynamics of HSCs and their products is through irradiation and transplant experiments (e.g., [59]), which also have important implications for bone marrow transplants [60, 61]. The typical experiment proceeds as follows [61]: a host animal receives a dose of irradiation sufficient to destroy the blood system (about 900 rad) but which leaves the rest of the organ systems intact (the intestinal system requires 1200-12,000 rad for irreversible damage, and the brain more than 12,000). This dose will kill the animal within 2-3 weeks, which gives a sense of the rate of decline of HSC descendants. Cells of a donor are isolated from the bone marrow by flushing sterile medium through the bone marrow and then injected into the host animal. The donor cells first circulate in the blood stream and then repopulate the bone marrow. The success of such transplants, from the perspective of the whole organism is demonstrated when a skin graft from the donor to the host is not rejected. In the period following transplantation, the donor HSCs produce progenitors that then differentiate into myeloid or lymphoid cells. Each differentiation that leads to a commitment of one route or another involves complex signaling [56, 57, 62, 63] and the signals themselves must be shaped by natural selection. We call this a transplant experiment.

A less commonly used approach, but equally important, is a perturbation experiment in which an animal is challenged in a way that reduces its complement of erythrocytes, platelets or granulocytes (commom myeloid progenitor (CMP) descendants) or lymphocytes (common lymphoid progenitor (CLP) descendants) and HSC activity is observed subsequent to the perturbation. For example, Cheshier et al. [10] bled mice on days 3, 6, and 9 and then sacrificed them on day 10 and measured the markers of HSC activity (see [64] for a similar experiment regarding the epidermis). Baldridge et al. [11] showed that quiescent HSCs were activated in response to an infection, another perturbation experiment. The models that we have described here can be used to predict dynamics of stem cells and their descendants for both kinds of experiments.

Conclusions

The use of quantitative models to understand the HSC system can be traced to the classic work of Till et al. [65], who used branching processes to interpret their experimental results on the variation of spleen colonies formed after a transplant experiment but evolutionary considerations rarely appear in models in stem cell biology [15]. The concepts developed for understanding how populations of individuals respond to these abiotic and biotic processes have colloraries at the cellular level. In particular, the size of stem cell systems is determined by the availability of molecular resources, feedbacks from differentiated cells, the effects of the epigenetic environment and the size of the environment capable of supporting stem cells [66]. Finally, population biology forces us to recognize the distinction between typological thinking (that all individuals of the same species are identical with constant characteristics and responses) and population thinking (that variation is real because every individual is unique and individual variation is central in ontogenetic and evolutionary history) [67]. While there are verbal and pictorial models (e.g.,[5, 68, 69]) of how signaling between differentiated cells and HSC may work, we have provided the first general quantitative theory linking how feedbacks through organismal need shapes the performance of the HSC system (cf [70]).

In conclusion, our work raises the critical question of how we connect the MPCR with the vast understanding on how signalling shapes HSC products. For instance, appreciating how the MPCR links to T-cell specification [71] and drives differential levels of the several well characterised precursors (of which lymphoid primed MPPs are one) requires further consideration of the both the key qualitative and quantitative positive and negative feedbacks in these sorts of signalling pathways. Similarly, linking the ideas of demand feedback to erythropoietin production [72], that differs between foetus and adult, requires systematic approach to combining the dynamical drivers with the lifetime fitness of the organism.

While the broad evolutionary ecological of HSC activity such as mounting an immune response to infection are well-known to have a cost on reproductive fitness (e.g., [73]), appreciating the ecological setting suggests that standard laboratory animals (bred for many generations under refined conditions) may not be the most appropriate organisms for the sorts of empirical tests of the theory we propose (also see [74] on metabolic morbidity of laboratory rodents). Studies on the evolutionary ecology of model organisms (see [75] for an approach using zebrafish) should have greater prominence in the understanding the population biology of stem cell responses.

Although raised almost two decades ago [76] this systems biology approach to stem cell dynamics must include evolutionary thinking. Approaching stem cell biology as population biology has much to offer both fields.

Declarations

Acknowledgements

We thank Camilla Forsberg, for many discussions on these topics, in which she endeavored to help us understand how the experiments are done, and Oliver Cinquin, German Encisco, Arthur Lander, John Lowengrub, and Edward Nelson for other very helpful conversations. The work of MM is supported by NSF grant EF-0924195.

Authors’ Affiliations

(1)
Department of Applied Mathematics and Statistics, University of California
(2)
Department of Biology, University of Bergen
(3)
Mathematical Ecology Research Group, Department of Zoology, University of Oxford and St Peter’s College

References

  1. Godin I, Cumano A: Hematopoietic Stem Cell Development. 2006, New York: Kluwer, 178View Article
  2. Kondo M: Hematopoietic Stem Cell Biology. 2010, New York: Humana Press, 264View Article
  3. Tan KY, Kim FS, Wagers AJ, Mayack SR: Hematopoeitic stem cells and somatic stem cells. Hematopoietic Stem Cell Biology. Edited by: Kondo M. 2010, New York: Humana Press, 57-92.View Article
  4. Metcalf D: The Molecular Control of Blood Cells. 1988, Cambridge: Harvard University Press, 192
  5. Ghaffari S: Oxidative stress in regulation of normal and neoplastic hematopoiesis. Antioxid Redox Sign. 2008, 10: 1923-1940. 10.1089/ars.2008.2142.View Article
  6. Guyton AC, Richardson TQ: Effect of hematocrit on venous return. Circ Res. 1961, 9: 157-164. 10.1161/01.RES.9.1.157.View ArticlePubMed
  7. Kanstrup I-L, Ekblom B: Blood volume and hemoglobin concentration as determinants of maximal aerobic power. Med Sci Sport Exer. 1984, 16: 258-262.View Article
  8. Horning M, Trillmich F: Development of hemoglobin, hematocrit, and erythrocyte values in Galapagos fur seals. Mar Mammal Sci. 1997, 13: 100-113. 10.1111/j.1748-7692.1997.tb00614.x.View Article
  9. Schuler B, Arras M, Keller S, Rettich A, Lundby C: Optimal hematorcrit for maximal exercise performance in acute and chronic erythropoietin-treated mice. Proc Natl Acad Sci USA. 2010, 107: 419-423. 10.1073/pnas.0912924107.PubMed CentralView ArticlePubMed
  10. Cheshier SH, Prohaska SS, Weissman IR: The effect of bleeding on hematopoeitic stem cell cycling and self-renewal. Stem Cell Devel. 2007, 16: 707-717. 10.1089/scd.2007.0017.View Article
  11. Baldridge MT, King KY, Boles NC, Weksberg DC, Goodell MA: Quiescent haematopoietic stem cells are activated by IFN-γ in response to chronic infection. Nature. 2010, 465: 793-798. 10.1038/nature09135.PubMed CentralView ArticlePubMed
  12. Rubtsov YP, Niec RE, Josefowicz S, Li L, Darce J: Stability of the regulatory T cell lineage in vivo. Science. 2010, 329: 1667-1671. 10.1126/science.1191996.PubMed CentralView ArticlePubMed
  13. Gaehtgens P, Kreutz F, Albrecht KH: Optimal hematocrit for canine skeletal muscle during rhythmic isotonic exercise. Eur J Appl Physiol. 1979, 41: 27-39. 10.1007/BF00424466.View Article
  14. Schoener TW: The newest synthesis: Understanding the interplay of evolutionary and ecological dynamics. Science. 2011, 331: 426-429. 10.1126/science.1193954.View ArticlePubMed
  15. Mangel M, Bonsall MB: Phenotypic evolutionary models in stem cell biology: replacement, quiescence, and variability. PLoS ONE. 2008, 3 (2): e1591-10.1371/journal.pone.0001591.PubMed CentralView ArticlePubMed
  16. Dercole F, Rinaldi S: Analysis of Evolutionary Processes. The Adaptive Dynamics Approach and Its Applications. 2008, Princeton: Princeton University Press, 352View Article
  17. Mangel M, Clark CW: Dynamic Modeling in Behavioral Ecology. 1988, Princeton: Princeton University Press, 308
  18. Mangel M, Ludwig D: Definition and evaluation of behavioral and developmental programs. Annu Rev Ecol Syst. 1992, 23: 507-536. 10.1146/annurev.es.23.110192.002451.View Article
  19. Houston AI, McNamara JM: Models of adaptive behavior. 1999, Cambridge: Cambridge University Press, 392
  20. Clark CW, Mangel M: Dynamic State Variable Models in Ecology: Methods and Applications. 2000, New York: Oxford University Press, 289
  21. Chickarmane CV, Enver T, Peterson C: Computational modeling of the hematopoietic erythroid-myeloid switch reveals insights into cooperativity, priming, and irreversibility. PLoS Comp Biol. 2009, 5 (1): e1000268-10.1371/journal.pcbi.1000268.View Article
  22. de Graaf CA, Kauppi M, Baldwin T, Hyland CD, Metcalf D: Regulation of hematopoietic stem cells by their mature progeny. Proc Natl Acad Sci USA. 2010, 107: 21689-21694. 10.1073/pnas.1016166108.PubMed CentralView ArticlePubMed
  23. Durand EM, Zon LI: The blood balance. Nature. 2010, 468: 644-645. 10.1038/468644a.View ArticlePubMed
  24. Gan B, Hu J, Jiang S, Liu Y, Sahin E: Lbk1 regulates quiescence and metabolic homoeostasis of haematopoietic stem cells. Nature. 2010, 468: 701-704. 10.1038/nature09595.PubMed CentralView ArticlePubMed
  25. Gurumurthy S, Xie SZ, Alagesan B, Kim J, Yusuf RZ: The Lbk1 metabolic sensor maintains haemotopoietic stem cell survival. Nature. 2010, 468: 659-663. 10.1038/nature09572.PubMed CentralView ArticlePubMed
  26. Nakada D, Saunders TL, Morrison SJ: Lkb1 regulates cell cycle and energy metabolism in haematopoietic stem cells. Nature. 2010, 468: 653-658. 10.1038/nature09571.PubMed CentralView ArticlePubMed
  27. Lander AD, Gokoffski KK, Wan FYM, Nie Q, Calof AL: Cell lineages and the logic of proliferative control. PLoS Biol. 2009, 7: e1000015-PubMed CentralView Article
  28. Wang L, Nie W, Enciso G: Non-essential sites improve phosphorylation switch. Biophys J. 2006, 99: L41-L43.View Article
  29. Abkowitz JL, Golinelli D, Harrison DE, Guttorp P: In vivo kinetics of murine hemopoietic stem cells. Blood. 2000, 96: 3399-3405.PubMed
  30. Lestas I, Vinnicombe G, Paulsson J: Fundamental limits on the suppression of molecular fluctuations. Nature. 2010, 467: 174-178. 10.1038/nature09333.PubMed CentralView ArticlePubMed
  31. Huang S: Systems biology of stem cells: three useful perspectives to help overcome the paradigm of linear pathways. Phil Trans R Soc B. 2011, 366: 2247-2259. 10.1098/rstb.2011.0008.PubMed CentralView ArticlePubMed
  32. Cressie N, Wikle CK: Statistics for Spatio-Temporal Data. 2011, Hoboken: John Wiley and Sons, 588
  33. Forsberg EC, Bhattacharya D, Weissman IL: Hematopoietic stem cells. Expression profiling and beyond. Stem Cell Rev. 2006, 2: 23-30.PubMed
  34. Bossdorf O, Richards CL, Pigliucci M: Epigentics for ecologists. Ecol Letts. 2008, 11: 106-115.
  35. Attema JL, Papthanasiou P, Forsberg EC, Xu J, Smale ST: Epigentic characterization of hematopoietic stem cell differentiation using miniChIP and bisulfite sequencing analysis. Proc Natl Acad Sci USA. 2007, 104: 12371-12376. 10.1073/pnas.0704468104.PubMed CentralView ArticlePubMed
  36. Mendez-Ferrer S, Lucas D, Battista M, Frenette PS: Haematopoietic stem cell release is regulated circadian oscillations. Nature. 2008, 452: 442-448. 10.1038/nature06685.View ArticlePubMed
  37. Thorpe JE, Mangel M, Metcalfe NB, Huntingford FA: Modelling the proximate basis of salmonid life-history variation, with application to Atlantic salmon, Salmo salar L. Evol Ecol. 1998, 12: 581-600. 10.1023/A:1022351814644.View Article
  38. Hawkins ED, Turner ML, Dowling MR, van Gend C, Hodgkin PD: A model of immune regulation as a consequence of randomized lymphocyte division and death times. Proc Natl Acad Sci USA. 2007, 104: 5032-5037. 10.1073/pnas.0700026104.PubMed CentralView ArticlePubMed
  39. Nijnik A, Woodbine L, Marchetti C, Dawson S, Lambe T: DNA repair is limiting for haematopoeitic stem cells during ageing. Nature. 2007, 447: 686-671. 10.1038/nature05875.View ArticlePubMed
  40. Rossi DJ, Bryder D, Seita J, Nussenzweig A, Hoeijmakers J: Deficiencies in DNA damage repair limit the function of haemotopoietic stem cells with age. Nature. 2007, 447: 725-730. 10.1038/nature05862.View ArticlePubMed
  41. Globerson A, Effros RB: Ageing of lymphocytes and lymphocytes in the aged. Immunol Today. 2000, 21: 515-521. 10.1016/S0167-5699(00)01714-X.View ArticlePubMed
  42. Rossi DJ, Bryder D, Weissman IL: Hematopoietic stem cell aging: mechanism and consequence. Exp Gerontol. 2007, 42: 385-390. 10.1016/j.exger.2006.11.019.PubMed CentralView ArticlePubMed
  43. Schlessinger D, Van Zant G: Does functional depletion of stem cells driving aging?. Mech Ageing Dev. 2001, 122: 1537-1553. 10.1016/S0047-6374(01)00299-8.View ArticlePubMed
  44. Aviv A, Levy D, Mangel M: Growth, telomere dynamics and successful and unsuccessful human aging. Mech Ageing Devel. 2003, 124: 829-837. 10.1016/S0047-6374(03)00143-X.View Article
  45. Kamminga LM, de Haan G: Cellular memory and hematopoietic stem cell aging. Stem Cells. 2006, 24: 1143-1149. 10.1634/stemcells.2005-0345.View ArticlePubMed
  46. Ruzankina Y, Asare A, Brown EJ: Replicative stress, stem cells and aging. Mech Ageing Dev. 2008, 129: 460-466. 10.1016/j.mad.2008.03.009.PubMed CentralView ArticlePubMed
  47. Wolpert L, Bedddington R, Jessell T, Lawrence P, Meyerowtiz E: Principles of Development,. 2002, Oxford: Oxford University Press, 569
  48. Passegue E, Wagers AJ, Giuriato S, Anderson WC, Weissman IL: Global analysis of proliferation and cell cycle gene expression in the regulation of hematopoietic stem and progenitor cell fates. J Exp Med. 2005, 202: 1599-1611. 10.1084/jem.20050967.PubMed CentralView ArticlePubMed
  49. Laird DJ, De Tomaso AW, Weissman IL: Stem cells are units of natural selection in a colonial ascidian. Cell. 2005, 123: 1351-1360. 10.1016/j.cell.2005.10.026.View ArticlePubMed
  50. Reiss MJ: The Allometry of Growth and Reproduction. 1989, Cambridge: Cambridge University Press, 200View Article
  51. Dingli D, Pacheco JM: Allometric scaling of the active hematopoietic stem cell pool across mammals. PLoS ONE. 2006, 1: 1-4. 10.1371/journal.pone.0000001.View Article
  52. Dingli D, Traulsen A, Pacheco JM: Compartmental architecture and dynamics of hematopoiesis. PLoS ONE. 2007, 4: 1-4.
  53. Dingli D, Traulsen A, Pacheco JM: Dynamics of haemopoiesis across mammals. Proc R Soc Lond B. 2008, 275: 2389-2392. 10.1098/rspb.2008.0506.View Article
  54. Abkowitz JL, Catlin SN, McCallie MT, Guttorp P: Evidence that the number of hematopoietic stem cells per animas is conserved in mammals. Blood. 2002, 100: 2665-2667. 10.1182/blood-2002-03-0822.View ArticlePubMed
  55. Dobzhansky T: Biology, molecular and organismic. Am Zool. 1964, 4: 443-452.View ArticlePubMed
  56. Ikawa T, Hirose S, Masuda K, Kakugawa K, Satoh R: An essential developmental checkpoint for the production of T cell lineage. Science. 2010, 329: 93-96. 10.1126/science.1188995.View ArticlePubMed
  57. Li L, Leid M, Rothenberg EV: An early T cell lineage commitment checkpoint dependent on the transcription factor Bc/11b. Science. 2010, 329: 89-93. 10.1126/science.1188989.PubMed CentralView ArticlePubMed
  58. Rothenberg EV, Zhang JL, Li L: Multilayered specification of the T-cell lineage fate. Immunol Rev. 2010, 238: 150-168. 10.1111/j.1600-065X.2010.00964.x.PubMed CentralView ArticlePubMed
  59. Forsberg EC, Serwold T, Kogan S, Weissman IL, Passegue E: New evidence supporting megakaryocyte-erythrocyte potential of Flk2/Flt3+ multipotent mematopoietic progenitors. Cell. 2006, 126: 415-426. 10.1016/j.cell.2006.06.037.View ArticlePubMed
  60. Appelbaum FR: Hematopoietic-Cell transplantation at 50. New Engl J Med. 2007, 357: 1472-1475. 10.1056/NEJMp078166.View ArticlePubMed
  61. Zipori D: Biology of Stem Cells and the Molecular Basis of the Stem State. 2009, New York: Humana Press, 304View Article
  62. Rothenberg EV: Negotiation of the T Lineage fate decision by transcription-factor interplay and microenvironmental signals. Immunity. 2007, 26: 690-702. 10.1016/j.immuni.2007.06.005.View ArticlePubMed
  63. Georgescu C, Longabaugh WJR, Scripture-Adams DD, David-Fung E-S, Yui MA: A gene regulatory network armature for T lymphocyte specification. Proc Natl Acad Sci USA. 2008, 105: 20100-20105. 10.1073/pnas.0806501105.PubMed CentralView ArticlePubMed
  64. Kasper M, Jaks V, Are A, Bergström A, Schwäger A: Wounding enhances epidermal tumorigenesis by recruiting hair follicle keratinocytes. Proc Natl Acad Sci USA. 2011, 108: 4099-4104. 10.1073/pnas.1014489108.PubMed CentralView ArticlePubMed
  65. Till JE, McCulloch EA, Siminovitch L: A stochastic model of stem cell proliferation based on the growth of spleen colony-forming cells. Proc Natl Acad Sci USA. 1964, 51: 29-36. 10.1073/pnas.51.1.29.PubMed CentralView ArticlePubMed
  66. Watt FM, Hogan BLM: Out of Eden: Stem cells and their niches. Science. 2000, 287: 1427-1430. 10.1126/science.287.5457.1427.View ArticlePubMed
  67. Mayr E: One Long Argument. Charles Darwin and the Genesis of Modern Evolutionary thought. 1991, Cambridge: Harvard University Press, 224
  68. Forsberg EC, Prohaska SS, Katzman S, Heffner GC, Stuart JM: Differential expression of novel potential regulators in hematopoietic stem cells. PLoS Genet. 2005, 1: 281-294.View Article
  69. Kopp H-G, Avecilla ST, Hooper AT, Rafii S: The bone marrow vascular niche: home of HSC differentiation and mobilization. Physiology. 2005, 20: 349-356. 10.1152/physiol.00025.2005.View ArticlePubMed
  70. Mangel M, Bonsall MB: The evolutionary ecology of stem cells and their niches - the time is now. Oikos. 2007, 116: 1779-1781. 10.1111/j.0030-1299.2007.16248.x.View Article
  71. Rothenberg EV, Moore JE, Yui MA: Launching the T-cell lineage developmental programme. Nat Rev Immunol. 2008, 8: 9-21. 10.1038/nri2232.PubMed CentralView ArticlePubMed
  72. Minamishima YA, Kaelin Jr WG: Reactivation of hepatic EPO synthesis in mice after PHD loss. Science. 2010, 329: 407-10.1126/science.1192811.PubMed CentralView ArticlePubMed
  73. Ahmed AM, Baggott SL, Maingon R, Hurd H: The costs of mounting an immune response are reflected in the reproductive fitness of the mosquito Anopheles gambiae. Oikos. 2002, 97: 371-377. 10.1034/j.1600-0706.2002.970307.x.View Article
  74. Martin B, Ji S, Maudsley S, Mattson MP: “Control” laboratory rodents are metabolically morbid: Why it matters. Proc Natl Acad Sci USA. 2010, 107: 6127-6133. 10.1073/pnas.0912955107.PubMed CentralView ArticlePubMed
  75. Engeszer RE, Patterson LB, Rao AA, Parichy DM: Zebrafish in the wild: A review of natural history and new notes from the field. Zebrafish. 2007, 4: 21-39. 10.1089/zeb.2006.9997.View ArticlePubMed
  76. Novak JP, Necas E: Proliferation-differentiation pathways of murine haematopoiesis: correlation of lineage fluxes. Cell Proliferat. 1994, 27: 597-633. 10.1111/j.1365-2184.1994.tb01377.x.View Article

Copyright

© Mangel and Bonsall; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement