Hierarchical Symmetry-Breaking Model for Stem Cell Differentiation

Abstract

Waddington envisioned stem cell differentiation as a marble rolling down a hill, passing through hierarchically branched valleys representing the cell’s temporal state. The terminal valleys at the bottom of the hill indicate the possible committed cells of the multicellular organism. Although originally proposed as a metaphor, Waddington’s hypothesis establishes the fundamental principles for characterizing the differentiation process as a dynamic system: the generated equilibrium points must exhibit hierarchical branching, robustness to perturbations (homeorhesis), and produce the appropriate number of cells for each cell type. This article aims to capture these characteristics using a mathematical model based on two fundamental hypotheses. First, it is assumed that the gene regulatory network consists of hierarchically coupled subnetworks of genes (modules), each modeled as a dynamical system exhibiting supercritical pitchfork or cusp bifurcation. Second, the gene modules are spatiotemporally regulated by feedback mechanisms originating from epigenetic factors. Analytical and numerical results show that the proposed model exhibits self-organized multistability with hierarchical branching. Moreover, these branches of equilibrium points are robust to perturbations, and the number of different cells produced can be determined by the system parameters.