Source-linked AI summary
From Microscopic Damage to Macroscopic Games: A Dimensionality Reduction of Stem Cell Homeostasis
Jiguang Yu, Louis Shuo Wang, Shihan Ban
TL;DR
The paper addresses the gap between microscopic damage-structured models and macroscopic game-theoretic descriptions of tissue homeostasis. It derives an exact low-dimensional closure under uniform mortality and shows that homeostasis forms an induced Nash equilibrium with equalized per-capita growth rates. Applied to murine intestinal crypts, the framework models dedifferentiation-driven regeneration and supports aggregate-data tests of its assumptions and predictions.
Problem
Existing models usually treat microscopic cellular wear or macroscopic population games separately, limiting rigorous links between damage dynamics and tissue-level stability.
Method
The paper derives an exact closure of a damage-structured PDE into a planar ODE under uniform TD mortality and interprets the resulting dynamics as an induced evolutionary game.
Results
The framework identifies tissue homeostasis with Nash equilibrium and applies its Ratio and Equalization Laws to dedifferentiation-driven murine intestinal regeneration.
Takeaways & Limitations
Aggregate stem and TD measurements can test the closed balance laws and estimate hard-to-measure dedifferentiation fluxes.
Takeaways & Limitations
Exact two-dimensional closure depends on uniform mortality; damage-dependent mortality makes mortality flux depend on the distribution of TD damage and requires additional information.
Abstract
from arXiv · showhide
Tissues must maintain macroscopic homeostasis despite the continuous microscopic accumulation of cellular damage. Theoretical models of this process often suffer from a disconnect between microscopic biophysics and macroscopic phenomenological games. Here, we bridge this gap by deriving an exact dimensionality reduction of a physiologically structured partial differential equation (PDE) into a low-dimensional dynamical system. Under the condition of uniform mortality, we mathematically demonstrate that tissue homeostasis operates as an induced Nash equilibrium, where the per-capita net growth rates of stem and differentiated phenotypes perfectly equalize. This reduction yields closed-form algebraic rules, the Ratio and Equalization Laws, that map continuous microscopic state dynamics to measurable macroscopic observables. To demonstrate the biological utility of this framework, we present a concrete, falsifiable case study of the murine intestinal crypt. By modeling crypt regeneration following irradiation-induced stem cell depletion, our framework successfully recovers the experimentally observed reliance on progenitor dedifferentiation. Furthermore, the model generates explicit, testable predictions, enabling the in vivo estimation of hard-to-measure lineage plasticity rates directly from aggregate static cell counts. This work provides a rigorous, predictive mathematical foundation for understanding how fast-renewing tissues filter microscopic noise to sustain macroscopic regenerative capacity.
Author summary
The paper bridges microscopic cell damage and macroscopic population games by translating complex damage dynamics into a simpler framework. It identifies homeostasis with balanced stem and mature-cell growth and applies the framework to intestinal regeneration.
- The framework bridges models of microscopic cellular wear and macroscopic population games, which are rarely treated together.
- Homeostasis is represented as a Nash equilibrium in which stem and mature cells have no competitive growth advantage.
- Application to mouse intestine regeneration predicts how tissues heal after severe injury.
- Aggregate cell-population snapshots can be used to calculate otherwise hidden tissue healing rates.
1 Introduction
The introduction frames a scale gap between damage-structured PDEs and phenomenological game theory, then proposes an exact closure under uniform mortality. It develops homeostatic game laws, thresholds, and falsifiable intestinal-crypt predictions.
- Rapidly renewing tissues must maintain macroscopic stability while cells accumulate damage and experience stress.
- Structured PDEs capture continuous damage and division but are infinite-dimensional and difficult to convert into simple tissue-level laws.
- Game-theoretic payoff functions are often phenomenological, leaving unclear whether tissue stability is biological or merely mathematical.
- Contributions: Under uniform TD mortality, the PDE total-mass balances close exactly into a planar ODE, making the induced game an emergent property of lineage dynamics.
- Contributions: Homeostasis is mathematically equivalent to a Nash equilibrium when stem and TD phenotypes equalize their per-capita net growth rates.
- Contributions: The intestinal-crypt case study models dedifferentiation-driven regeneration and estimates dedifferentiation flux from static aggregate flow-cytometry observables.
- Falsifiable readouts: The framework is falsifiable because measured totals must satisfy closed balance laws, while perturbations can test the predicted extinction–growth threshold.
- Falsifiable readouts: At steady state, the Ratio and Equalization Laws connect stem-to-TD abundance, mortality, feedback, and dedifferentiation rates to aggregate measurements.
2 Damage-structured PDE model
The model describes stem and TD populations structured by accumulating damage, with division, dedifferentiation, mortality, and feedback-regulated fate probabilities. It uses transport and characteristic methods to formulate the PDE and establish global dynamics.
- Stem and TD densities depend on time and a scalar damage coordinate, with total masses obtained by integrating over damage.
- Damage increases through transport, while division produces stem, TD, or mixed daughters with deterministic inherited damage fractions.
- TD cells dedifferentiate into the stem compartment, and their mortality may depend on damage through δ(x).
- Feedback regulation: Hill-type feedback maps regulate self-renewal and differentiation probabilities as tissue burden changes.
- Mortality assumptions: The model distinguishes general bounded mortality, used for well-posedness and balance laws, from uniform mortality, used for exact closure.
- Mortality assumptions: Uniform mortality is a structural identifiability condition rather than a claim that real tissues have damage-independent death.
- Well-posedness: Characteristic trajectories and integrating factors support a mild formulation, while Lipschitz feedback, positivity, and mass bounds yield a unique global solution.
3 Balance laws and exact reduction
The PDE balance laws close exactly to a two-dimensional ODE when terminally differentiated-cell mortality is uniform. Variable mortality instead couples total dynamics to the damage distribution, making finite-dimensional closure generally approximate and empirically testable.
- Balance laws: Theorem 3.1 establishes unique global nonnegative PDE solutions and balance laws for the total masses.The balance laws hold for almost every t > 0 under the stated assumptions.
- Failure of closure: The balance laws are not closed in general because mortality depends on the full differentiated-cell distribution W(t, ·), not only its total mass.Additional weighted moments are then required, and finite-dimensional closures are typically approximate.
- Exact reduction: Uniform mortality δ(x) ≡δ converts the mortality flux into δ ¯W(t), yielding an exact planar ODE for total stem and differentiated-cell masses.The nonnegative quadrant is forward invariant and solutions are global.
- Interpretation and testing: Uniform mortality is a structural identifiability condition for exact closure, not a literal claim that differentiated-cell death is damage-independent in vivo.Exact closure requires the mortality flux to be expressible as a function of ¯W(t) alone.
- Interpretation and testing: Aggregate trajectories that violate the closed balance laws beyond discretization error empirically falsify the uniform-mortality closure for that regime.The figure contrasts total-mass-dependent mortality with distribution-shape-dependent mortality.
4 Replicator mapping and game-theoretic interpretation
Under uniform mortality, the exact reduced dynamics induce a state-dependent replicator game whose payoffs are phenotype-specific per-capita net growth rates. Homeostasis corresponds to payoff equalization, while the Ratio and Equalization Laws turn this condition into experimentally testable relations among aggregate observables and effective rates.
- Replicator mapping: The exact two-compartment closure induces replicator dynamics with stem and differentiated phenotypes as strategies and per-capita net growth rates as state-dependent payoffs.The game-theoretic structure is derived from PDE lineage bookkeeping rather than postulated.
- Replicator mapping: The sign of FP −FW determines selection on stemness: FP > FW increases stem frequency s, whereas otherwise s decreases.This follows from the replicator equation ˙s = s(1 −s)(FP −FW ).
- Nash/homeostasis equivalence: At homeostasis, payoff equalization gives stem and differentiated phenotypes identical per-capita net growth, so neither has an aggregated selective advantage.This is the model’s induced Nash-equilibrium interpretation.
- Nash/homeostasis equivalence: An interior equilibrium occurs if and only if FP = FW = 0, while FP = FW alone identifies a replicator rest point for the stem frequency.The additional zero common growth-rate condition distinguishes a fixed point from frequency stationarity.
- Closed-form laws and observables: The Ratio and Equalization Laws express equilibrium through algebraic identities among measurable compartment sizes and effective division, death, and dedifferentiation rates.Steady-state sizes can be measured by cell counting or flow cytometry, while rates can be estimated using proliferation, death, and lineage-tracing assays.
- Closed-form laws and observables: Once δ and λP(¯W∗) are estimated, the steady-state ratio and dedifferentiation rate required to offset differentiation losses become directly testable predictions.The framework therefore connects aggregate measurements to effective lineage plasticity rates.
5 Dynamics: thresholds, uniqueness, and global stability
Under uniform differentiated-cell mortality, the PDE balance laws reduce exactly to a planar system whose dynamics feature an extinction–growth threshold, conditions for at most one interior equilibrium, and global convergence without oscillations. Uniform gain scaling changes total abundance while preserving the stem-to-TD composition.
- Dimensionality reduction: Uniform TD death closes the PDE balance laws exactly to a planar ODE system on total stem and differentiated-cell populations.The closed quadrant is forward invariant and solutions are global.
- Extinction–growth threshold: The threshold ∆pcrit = −ˆλR/δ separates locally stable extinction from a growth regime in which the origin is a saddle with a one-dimensional unstable manifold.At equality, the determinant vanishes and one eigenvalue is zero.
- Uniqueness: The Ratio law defines a ratio curve, and substituting it into the Equalization law produces a scalar root condition for interior equilibria.Under the stated slope condition, the residual is strictly decreasing and has at most one root.
- Global stability: Negative divergence excludes nontrivial periodic orbits, and every bounded interior trajectory converges to the unique interior equilibrium when it exists.The global conclusion is stronger than a local eigenvalue calculation.
- Scaling invariance: Uniform amplification of Hill feedback gains preserves the equilibrium stem-to-TD ratio while rescaling total abundances by 1/A.The scaling symmetry separates tissue composition from the level of homeostasis.
6 Numerical methods and reproducible verification
The numerical workflow pairs conservative PDE discretization with certified checks against the exact two-compartment ODE closure and analytical predictions. Refinement, domain, residual, equilibrium, phase-plane, replicator, bifurcation, and gain-scaling tests support the reported behavior while exposing the uniform-mortality scope.
- Verification framework: The verification framework separates numerical equation-solving accuracy from certified checks of theorem-predicted qualitative behavior.Deterministic, parameter-locked scripts emit the reported quantities and map analytical results to computational tests.
- Numerical schemes: Conservative finite-volume remapping preserves the mass identity required for balance-law residuals and PDE-versus-ODE closure checks.The method uses exact overlap fractions rather than pointwise interpolation, while negative-value corrections are monitored through induced mass defects.
- Verification tests: At T = 5, balance-law residuals are O(10^-3), decrease under grid refinement, and remain independent of xmax ∈{10, 15, 20}.The residual magnitudes confirm the balance laws to the accuracy of the spatial discretization.
- Closure verification: With xmax = 20, PDE-versus-ODE total-mass discrepancies decrease at the expected first-order rate under grid refinement, with zero positivity-fix mass defect.The remaining discrepancy is attributed to the O(∆x) upwind scheme rather than closure failure.
- Analytical predictions: Eight trajectories converge monotonically to the unique interior equilibrium, while the replicator field has a stable Nash frequency and the divergence is strictly negative throughout the sampled domain.The gain-scaling test additionally confirms equilibrium mass rescaling by 1/A and invariant composition, with ratio errors ≤3 × 10^-15.
- Parameter scans: The extinction-growth scan resolves the threshold ∆pcrit = −0.6, and increasing δ shifts the equilibrium toward smaller total mass.The threshold separates the extinction regime from the positive equilibrium branch in the admissible parameter window.
7 Case study: murine intestinal crypt homeostasis and regeneration
The calibrated two-compartment reduction represents the murine crypt with aggregate stem and differentiated compartments, reproducing the homeostatic target and modeling post-irradiation recovery through progenitor dedifferentiation. It yields testable predictions while explicitly limiting interpretation because transit-amplifying cells are lumped into the effective differentiated compartment and recovery trajectories are semi-quantitative.
- Calibration: The calibration targets 14 stem cells and 278 differentiated cells, corresponding to a total crypt cellularity of approximately 292 cells.The 14-cell stem target is consistent with approximately 14 ± 2 Lgr5hi stem cells reported for murine small-intestinal crypts.
- Model setup: The crypt is modeled by identifying P with Lgr5+ CBC stem cells and W with all remaining crypt-resident differentiated cells, including transit-amplifying progenitors.This is an aggregate coarse-graining rather than a one-to-one transcriptional cell-state model.
- Calibration: 4.9643 day^-1 is the calibrated aggregate effective production rate from the Ratio law, not the literal single-cell Lgr5+ CBC cycling rate.The reduction incorporates the transit-amplifying cascade into the effective differentiated compartment.
- Homeostatic equilibrium: The calibrated equilibrium is (P̄*, W̄*) = (14.0000, 278.0000), with equalization and ratio residuals below 10^-17.The reduced system reproduces the prescribed aggregate homeostatic state to floating-point precision.
- Homeostatic equilibrium: The equilibrium stem fraction is 4.79%, where stem and differentiated phenotypes have equal per-capita net growth rates.Thus, feedback-regulated lineage dynamics maintain the stem phenotype without intrinsic competitive dominance.
- Regeneration: After irradiation-like stem depletion, recovery is initially driven by differentiated-to-stem return flux, matching experimentally documented progenitor dedifferentiation.The model predicts a transient positive stem-versus-differentiated payoff difference that restores the homeostatic composition.
- Predictions and limitations: The case study generates experimentally testable predictions linking epithelial loss, stem-to-differentiated ratios, and lineage plasticity rates.The framework is presented as a proof-of-principle rather than a formal statistical fit to one harmonized cohort.
- Predictions and limitations: The deliberate lumping of transit-amplifying cells prevents direct comparison between the effective production rate and the CBC single-cell cycling rate.A three-compartment stem/TA/TD model would map more literally to crypt biology but would abandon the exact two-dimensional closure.
8 Discussion
The discussion interprets the Ratio and Equalization laws as experimentally measurable balance statements and clarifies the scope of the exact reduction. The framework requires damage-independent mortality for exact closure, while damage-dependent mortality may still permit homeostatic stability but requires richer approximations.
- Measurable laws: The Ratio law links the steady-state stem-to-differentiated ratio to differentiated-cell clearance and the effective stem-cycle rate.The clearance rate can be measured with labeled-cohort survival or apoptosis assays, while the effective cycle rate can be estimated from EdU/BrdU incorporation.
- Measurable laws: The Equalization law provides a balance statement for the differentiated-to-stem return flux required to maintain steady tissue composition.Lineage tracing or fate mapping can estimate the return rate under steady state or controlled micro-injury.
- Experimental implications: Increasing effective differentiated-to-stem return should shift the steady ratio along the algebraic laws.This provides a perturbation logic for testing plasticity experimentally.
- Scope of the reduction: The exact evolutionary-game derivation relies on damage-independent mortality, δ(x) ≡ δ.This assumption corresponds biologically to removal driven primarily by extrinsic stochastic factors such as mechanical shedding or random differentiation-linked loss.
- Scope of the reduction: When mortality depends on damage, exact moment closure fails, although homeostatic stability may not necessarily be destroyed.The reduced system is therefore a transparent theoretical benchmark, while future moment-closure approximations could extend the framework.
A.1 Notation and auxiliary bounds
The notation section imposes boundedness and Lipschitz regularity on the feedback maps, providing global bounds used in the model’s analysis.
- Regularity assumptions: The feedback maps are bounded and Lipschitz on [0, ∞) under Assumption 2.1.This regularity condition supports the stated analytical framework.
- Lipschitz control: The complementary feedback map p_3 is Lipschitz with constant L_p3 := L_p1 + L_p2.This follows from p_3(·) = 1 − p_1(·) − p_2(·).
- Global bounds: The model bounds p_i between 0 and p̂_i, p_3 between 0 and 1, and the rates λ_P and λ_R by their hatted maxima.These are the global bounds labeled (A1).
A.2 Kernel mass identities
The auxiliary kernel identities establish the mass-preserving properties of the nonlocal birth operators used in the balance-law reduction.
- Mass preservation: The reduction relies on nonlocal birth operators preserving L1-mass in each deterministic inheritance branch.This is the key structural ingredient for deriving total-mass balance laws.
- Change of variables: Lemma A.1 supplies the change-of-variables identity for α ∈ (0, 1) and integrable f.The substitution y = x/α with dx = α dy underlies the kernel calculations.
- Mass formulas: Lemma A.2 obtains the birth-operator mass formulas by applying Lemma A.1 term by term and using p_1 + p_2 + p_3 = 1.The partition-of-probability identity closes the aggregate mass accounting.
A.3 Vanishing of transport boundary flux under homogeneous inflow
The lemma formalizes cancellation of transport boundary terms when integrating over x ≥0 under homogeneous inflow. The proof uses compactly supported test functions, integrability at infinity, and zero inflow at x = 0.
- Boundary flux cancellation: Under homogeneous inflow u(t, 0) = 0, integrating the transport equation over x ≥0 produces no boundary injection at x = 0.The result applies to weak solutions with v > 0 and integrable source term g.
- Conclusion: The resulting identity holds for almost every t in (0, T) in the distributional sense used to derive balance laws.
- Proof strategy: A cutoff test function ψR equal to one on [0, R] and supported in [0, R + 1] localizes the weak formulation before taking R →∞.The cutoff derivative contribution vanishes by dominated convergence because u(t, ·) ∈L1.