Cells, cancer, and rare events: homeostatic metastability in stochastic non-linear dynamics models of skin cell proliferation

A recently proposed single progenitor cell model for skin cell proliferation [Clayton et al., Nature v446, 185 (2007)] is extended to incorporate homeostasis as a fixed point of the dynamics. Unlimited cell proliferation in such a model can be viewed as a paradigm for the onset of cancer. A novel way in which this can arise is if the homeostatic fixed point becomes metastable, so that the cell populations can escape from the homeostatic basin of attraction by a large but rare stochastic fluctuation. Such an event can be viewed as the final step in a multi-stage model of carcinogenesis. This offers a possible explanation for the peculiar epidemiology of lung cancer in ex-smokers.

💡 Research Summary

The paper builds on the single‑progenitor model of skin epidermal turnover introduced by Clayton et al. (Nature 2007) and reformulates it as a stochastic nonlinear dynamical system with an explicit homeostatic fixed point. In the original framework, two progenitor states (A and B) and a differentiated state (C) interconvert with constant rates, leading to a steady‑state balance of cell numbers. The authors replace the constant transition rates with nonlinear functions of the current cell fractions and external cues, thereby allowing multiple fixed points. One of these fixed points corresponds to physiological homeostasis: the fractions of A, B, and C remain constant and the total cell population is on average stable.

Because cell division, differentiation, and death are intrinsically stochastic, the authors embed intrinsic noise using a chemical master equation and simulate trajectories with the Gillespie algorithm. When the noise amplitude is modest, trajectories fluctuate around the homeostatic point and remain confined within its basin of attraction. However, if the homeostatic fixed point becomes only metastable—i.e., it retains local stability but the barrier separating it from an unbounded growth regime is low—rare large fluctuations can push the system over the barrier. Once this occurs, the progenitor pool expands unchecked, differentiation is suppressed, and the system enters an “unlimited proliferation” regime that mathematically mirrors the onset of a tumor.

The probability of such an escape event is analyzed with large‑deviation theory. The escape rate scales as P ≈ exp(−ΔS/ε), where ΔS is the action associated with the optimal fluctuation path and ε quantifies noise strength. Numerical simulations confirm the exponential dependence and show that, although the event is exceedingly rare on short timescales, the cumulative probability over a human lifetime becomes non‑negligible.

By mapping this stochastic escape onto the classic multi‑stage carcinogenesis paradigm, the authors propose a new interpretation of the final stage: instead of a deterministic accumulation of mutations, the decisive step may be a stochastic breach of homeostatic stability. Early mutational hits (e.g., caused by tobacco‑related DNA damage) can subtly alter the nonlinear transition functions, lowering the barrier and rendering the homeostatic point metastable. Subsequent stochastic fluctuations—perhaps amplified by chronic inflammation or other micro‑environmental stresses—can then trigger the catastrophic transition to uncontrolled growth.

The model offers a mechanistic explanation for the persistent incidence of lung cancer among former smokers. Even after cessation, residual genetic lesions keep the epidermal homeostatic fixed point precariously close to instability, allowing a rare but possible stochastic transition that manifests as cancer years later. This perspective reconciles epidemiological data showing a long‑tail risk in ex‑smokers with a biologically plausible dynamical mechanism.

Finally, the authors outline experimental avenues for validation: (i) measuring cell‑state transition rates at single‑cell resolution to infer the shape of the nonlinear functions; (ii) manipulating noise levels in mouse skin (e.g., via controlled inflammation) to observe metastable escape events; and (iii) long‑term cohort studies to quantify the statistical signature of rare stochastic transitions. Together, these approaches could calibrate the model, test its predictions, and potentially inform preventive strategies that aim to reinforce the stability of the homeostatic attractor.