A system-level model reveals that transcriptional stochasticity is required for hematopoietic stem cell differentiation

HSCs differentiation has been difficult to study experimentally due to the high number of components and interactions involved, as well as the impact of diverse physiological conditions. From a 200-node network, that was grounded on experimental data, we derived a 21-node regulatory network by collapsing linear pathways and retaining the functional feedback loops. This regulatory network core integrates key nodes and interactions underlying HSCs differentiation, including transcription factors, metabolic, and redox signaling pathways. We used Boolean, continuous, and stochastic dynamic models to simulate the hypoxic conditions of the HSCs niche, as well as the patterns and temporal sequences of HSCs transitions and differentiation. Our findings indicate that HSCs differentiation is a plastic process in which cell fates can transdifferentiate among themselves. Additionally, we found that cell heterogeneity is fundamental for HSCs differentiation. Lastly, we found that oxygen activates ROS production, inhibiting quiescence and promoting growth and differentiation pathways of HSCs.