Tutorials › Quantitative Systems Pharmacology (QSP) › Stem, Progenitor, and Differentiated Compartments

Cell Population Dynamics · Tutorial 168 of 799

Stem, Progenitor, and Differentiated Compartments

Learn to construct and interpret a hierarchical lineage model by tracking stem-cell renewal, progenitor amplification, terminal differentiation, and cell loss as separate fluxes.

Intermediate 10 min read

What You'll Learn

  • Write balance equations for stem, progenitor, and differentiated cell compartments.
  • Distinguish stem self-renewal from commitment into the progenitor compartment.
  • Interpret progenitor amplification as a source of additional progenitor cells.
  • Track terminal differentiation as a transfer between compartments rather than a loss from the modeled lineage.
  • Check how differentiation fluxes cancel in the total-cell balance.
  • Identify assumptions and common interpretation errors in a three-compartment lineage model.

Why Separate a Lineage Into Compartments?

A total cell count can conceal important changes in population structure. The same total might consist mostly of stem cells at one time and mostly of differentiated cells at another. A useful lineage model therefore distinguishes cell states that have different roles: stem cells maintain the lineage, progenitors expand and supply downstream cells, and differentiated cells represent the terminal output being tracked.

The goal is not to describe every individual division. Instead, a compartment model records the rates at which cells enter, leave, or are added within each state. Keeping these processes distinct makes it possible to ask whether a changing differentiated-cell count reflects altered supply from progenitors, altered loss from the differentiated pool, or both.

Definition: A stem compartment contains cells capable of maintaining the lineage; a progenitor compartment contains cells that expand and contribute to downstream production; and a differentiated compartment contains the terminal cell state represented in the model. A lineage model specifies the cell fluxes between these compartments and any sources or losses within them.

Define the Compartments and Fluxes

Let \(N_S(t)\), \(N_P(t)\), and \(N_D(t)\) denote the numbers of stem, progenitor, and differentiated cells at time \(t\). Consider the unidirectional hierarchy \(S\to P\to D\). In this model, stem cells contribute cells to the progenitor compartment, progenitors contribute differentiated cells, and no cells return upstream.

We distinguish two types of input from the transitions. Self-renewal supplies the stem compartment, helping maintain its population. Progenitor amplification adds progenitor cells within the progenitor compartment. For a simple linear model, represent those effective additions by \(r_S N_S\) and \(r_P N_P\), respectively. Here \(r_S\) and \(r_P\) have units of inverse time. They describe effective net cell additions per source cell per unit time, not counts of individual division events.

The transition flux from stem to progenitor is \(k_{SP}N_S\), and the transition flux from progenitor to differentiated cells is \(k_{PD}N_P\). The coefficients \(k_{SP}\) and \(k_{PD}\) also have units of inverse time. Let \(L_S\), \(L_P\), and \(L_D\) denote the respective cell-loss fluxes, measured in cells per unit time. They are written as fluxes rather than assigned a particular loss mechanism or rate law.

Model assumptions: The three compartments are well-mixed aggregate populations; the modeled lineage moves only from \(S\) to \(P\) to \(D\); the stated inputs and transitions are the only sources or transfers represented; and losses remove cells from the modeled compartments. The linear input and transition expressions apply when their coefficients are constant and the per-cell flux assumptions are appropriate. The loss fluxes may be specified separately, without assuming a particular loss mechanism.

With these definitions, the model is a set of cell-count balances:

$$ \begin{aligned} \frac{dN_S}{dt} &= r_S N_S - k_{SP}N_S - L_S,\\ \frac{dN_P}{dt} &= k_{SP}N_S + r_P N_P - k_{PD}N_P - L_P,\\ \frac{dN_D}{dt} &= k_{PD}N_P - L_D. \end{aligned} $$

Each term has units of cells per unit time. Self-renewal adds cells to \(S\); commitment from \(S\) removes cells from \(S\) and supplies \(P\). Amplification adds cells within \(P\); differentiation removes cells from \(P\) and supplies \(D\). Loss removes cells from the compartment in which it occurs. This accounting makes clear that differentiation is not itself a loss from the whole modeled lineage: it is a transfer between compartments.

Worked Example: Calculate the Three Compartment Rates of Change

Suppose a hypothetical population is observed at an instant with \(N_S=100\) cells, \(N_P=200\) cells, and \(N_D=1{,}000\) cells. Let \(r_S=0.020\) per day, \(k_{SP}=0.010\) per day, \(r_P=0.10\) per day, and \(k_{PD}=0.050\) per day. At that instant, suppose the loss fluxes are \(L_S=1\), \(L_P=8\), and \(L_D=20\) cells per day.

The self-renewal input is \(r_SN_S=(0.020)(100)=2\) cells per day, and the stem-to-progenitor flux is \(k_{SP}N_S=(0.010)(100)=1\) cell per day. Thus,

$$ \frac{dN_S}{dt}=2-1-1=0\ \text{cells/day}. $$

The progenitor amplification input is \(r_PN_P=(0.10)(200)=20\) cells per day. The progenitor-to-differentiated flux is \(k_{PD}N_P=(0.050)(200)=10\) cells per day. Therefore,

$$ \frac{dN_P}{dt}=1+20-10-8=3\ \text{cells/day}, \qquad \frac{dN_D}{dt}=10-20=-10\ \text{cells/day}. $$

At this instant, the stem count has zero net rate of change, the progenitor count is increasing by 3 cells per day, and the differentiated count is decreasing by 10 cells per day. These are instantaneous rates, not predictions that the counts will change linearly at those rates over an arbitrary future interval.

What Happens to the Total Cell Count?

Adding the three balance equations gives a useful check on the model. Each transition appears once as an outflow and once as an inflow: \(k_{SP}N_S\) leaves \(S\) and enters \(P\), while \(k_{PD}N_P\) leaves \(P\) and enters \(D\). The two contributions therefore cancel in the total.

$$ \begin{aligned} \frac{d}{dt}(N_S+N_P+N_D) &=r_SN_S-k_{SP}N_S-L_S\\ &\quad +k_{SP}N_S+r_PN_P-k_{PD}N_P-L_P\\ &\quad +k_{PD}N_P-L_D\\ &=r_SN_S+r_PN_P-(L_S+L_P+L_D). \end{aligned} $$

This identity follows directly from cell accounting: internal transfers redistribute cells but do not change the total number across the three compartments. In the preceding example, total change is \(2+20-(1+8+20)=-7\) cells per day. The same result follows by adding the compartment rates, \(0+3-10=-7\) cells per day.

The cancellation is also a diagnostic. If a proposed set of equations makes a one-way differentiation flux decrease total cell number without an explicit loss, the transition has likely been counted as an outflow without its matching inflow. Conversely, if a transition is added to both compartments as an input, it has not been removed from its source.

Worked Example: How a Change in Differentiation Redistributes Cells

Use the hypothetical counts and loss fluxes from the first example, but compare two possible values of \(k_{PD}\): \(0.050\) and \(0.080\) per day. Hold all other quantities fixed. The progenitor count is \(200\), so the differentiation flux is 10 cells per day in the first case and 16 cells per day in the second.

At \(k_{PD}=0.050\) per day, the relevant rates are

$$ \frac{dN_P}{dt}=1+20-10-8=3,\qquad \frac{dN_D}{dt}=10-20=-10 \quad\text{cells/day}. $$

At \(k_{PD}=0.080\) per day, the relevant rates become

$$ \frac{dN_P}{dt}=1+20-16-8=-3,\qquad \frac{dN_D}{dt}=16-20=-4 \quad\text{cells/day}. $$

The larger transition flux moves more cells out of \(P\) and into \(D\), so the progenitor rate of change falls by 6 cells per day while the differentiated rate rises by 6 cells per day. The total rate is unchanged: it remains \(-7\) cells per day because the other inputs and losses did not change. A higher differentiation flux therefore does not guarantee an increasing differentiated-cell count; here the differentiated count still declines because its loss flux exceeds its incoming flux.

Interpret Fluxes Before Interpreting Counts

The rates \(r_S\) and \(r_P\) represent distinct biological roles in this coarse-grained model. A change in \(r_S\) changes the supply to the stem compartment; a change in \(r_P\) changes the addition of progenitor cells. By contrast, \(k_{SP}\) redirects cells from the stem compartment to the progenitor compartment, and \(k_{PD}\) redirects cells from progenitors to differentiated cells. Changes in those transition coefficients affect both the source and destination compartments, with opposite signs.

A differentiated-cell count is consequently an output of upstream supply and local loss, not a direct readout of any one process. In particular, a low \(N_D\) alone does not establish that differentiation is slow. It could also arise from a small progenitor pool, weak progenitor amplification, or substantial differentiated-cell loss. The model helps organize these explanations, but deciding among them requires information about counts or fluxes beyond a single terminal measurement.

Worked Example: Separate Differentiated-Cell Production From Loss

In a hypothetical system, suppose \(N_P=300\) progenitor cells and \(N_D=500\) differentiated cells are present at an instant. Let \(k_{PD}=0.040\) per day, and let the differentiated-cell loss flux be 15 cells per day. The production flux is

$$ k_{PD}N_P=(0.040)(300)=12\ \text{cells/day}. $$

The differentiated-cell balance is therefore

$$ \frac{dN_D}{dt}=12-15=-3\ \text{cells/day}. $$

The differentiated compartment is declining at 3 cells per day at that instant, even though 12 new differentiated cells per day are being supplied. The count is not declining because production is absent; production is simply smaller than loss in this balance. If the loss flux were instead 9 cells per day, with the same production, the rate would be \(12-9=3\) cells per day. This comparison isolates the balance of incoming and outgoing cells without inferring an unmeasured mechanism from the count alone.

Assumptions, Diagnostics, and Alternative Representations

The three-compartment structure is useful when the main question concerns the balance of stem maintenance, progenitor expansion, and terminal output. Its simplicity depends on assumptions that should be checked against the intended application:

  • Compartment meaning is consistent. Each state must group cells whose modeled roles are sufficiently similar. If a “progenitor” label combines states with very different expansion or differentiation behavior, one effective rate may not describe them well.
  • The lineage is treated as one-way. The equations omit dedifferentiation, reversible switching, and branching into separate fates. If those processes matter, the specified hierarchy is incomplete.
  • Rates and fluxes have appropriate units and signs. Transition and addition coefficients are nonnegative and have units of inverse time; loss fluxes are nonnegative and have units of cells per time. A negative loss flux would represent an unmodeled source, not ordinary loss.
  • The coarse-grained inputs have a defensible interpretation. The \(r_S N_S\) and \(r_P N_P\) terms are effective net additions in their respective compartments. They should not be interpreted automatically as raw division frequencies.

A practical diagnostic is to write a flux table before fitting or simulating the equations: list each source and destination, its units, and whether it adds cells, transfers them, or removes them. Then verify the total-count identity. If observed trajectories cannot be represented without implausible negative fluxes or inconsistent compartment meanings, the model structure or its inputs may be inadequate.

Alternative representations are appropriate for different questions. A maturation chain separates several successive intermediate stages when their residence times or counts matter. A branching-flux model is more suitable when the important feature is a split into distinct downstream fates. A single total-population model is more compact, but it cannot describe how cells are distributed among stem, progenitor, and differentiated states. The three-compartment model is a middle ground: it preserves these functional distinctions without representing every division or every intermediate state.

A common mistake is to treat differentiation as both a loss from the whole lineage and a source of differentiated cells. The incorrect reasoning is that \(k_{PD}N_P\) should be subtracted from the total because it leaves the progenitor pool, while also being added to the differentiated pool. The correction is to distinguish a compartment outflow from a lineage loss: differentiation leaves \(P\) but enters \(D\), so it cancels in the total balance. Only fluxes that leave the modeled lineage reduce its total cell count.

Key takeaway: Build the lineage from explicit fluxes. Self-renewal supplies the stem compartment, progenitor amplification adds progenitors, differentiation transfers cells downstream, and loss removes cells from the modeled pools. This separation clarifies what each rate can—and cannot—explain about observed cell counts.

Check Your Understanding

Use the model and definitions above to answer each question. Solutions follow each problem.

  1. A hypothetical system has \(N_S=80\), \(r_S=0.030\) per day, \(k_{SP}=0.010\) per day, and \(L_S=1\) cell per day. What is the instantaneous rate of change in the stem compartment?
  2. If the progenitor-to-differentiated flux is 14 cells per day, what sign does this flux have in the progenitor and differentiated equations, and what is its contribution to the total-cell rate?
  3. A hypothetical differentiated compartment receives 8 cells per day and loses 11 cells per day. Calculate its rate of change and interpret it in context.
  4. Why can a larger \(k_{PD}\) increase the differentiated-cell production flux without ensuring that the differentiated count increases?

Solutions. (1) Self-renewal supplies \((0.030)(80)=2.4\) cells per day, while stem commitment removes \((0.010)(80)=0.8\) cell per day and loss removes 1 cell per day. Thus \(dN_S/dt=2.4-0.8-1=0.6\) cell per day. (2) The flux contributes \(-14\) cells per day to the progenitor balance and \(+14\) cells per day to the differentiated balance, so its net contribution to the total is zero. (3) The rate is \(8-11=-3\) cells per day; at that instant, the differentiated count is declining because loss exceeds supply by 3 cells per day. (4) The differentiated count depends on the difference between incoming differentiation and differentiated-cell loss. A larger incoming flux can still be smaller than the loss flux, leaving the differentiated compartment declining.

Clinical Biostats