Mathematical model of a gene regulatory network reconciles effects of genetic perturbations on hematopoietic stem cell emergence.

Interlinked gene regulatory networks (GRNs) are vital for the spatial and temporal control of gene expression during development. The hematopoietic transcription factors (TFs) Scl, Gata2 and Fli1 form one such densely connected GRN which acts as a master regulator of embryonic hematopoiesis. This triad has been shown to direct the specification of the hemogenic endothelium and emergence of hematopoietic stem cells (HSCs) in response to Notch1 and Bmp4-Smad signaling. Here we employ previously published data to construct a mathematical model of this GRN network and use this model to systematically investigate the network dynamical properties. Our model uses a statistical-thermodynamic framework to describe the combinatorial regulation of gene expression and reconciles, mechanistically, several previously published but unexplained results from different genetic perturbation experiments. In particular, our results demonstrate how the interactions of Runx1, an essential hematopoietic TF, with components of the Bmp4 signaling pathway allow it to affect triad activation and acts as a key regulator of HSC emergence. We also explain why heterozygous deletion of this essential TF, Runx1, speeds up the network dynamics leading to accelerated HSC emergence. Taken together our results demonstrate that the triad, a master-level controller of definitive hematopoiesis, is an irreversible bistable switch whose dynamical properties are modulated by Runx1 and components of the Bmp4 signaling pathway.