Performing nonlinear stability analysis for a highly interconnected gene network model

Hematopoiesis, the process of forming all the types of blood cells, is a good model to study cell-fate choice. Inferring hematopoietic gene regulatory networks (GRNs) using gene circuit models is a crucial step in unraveling the complex dynamics governing hematopoiesis. A study has been carried out in which, a 12-gene model was used to model the dynamics of the gene regulatory network that governs differentiation into erythrocytes and neutrophils. We apply nonlinear stability analysis to this highly dimensional gene network model, to figure out qualitative behaviors in terms of stability. This will offer valuable insights into hematopoietic gene regulation and guide the development of targeted therapeutic interventions.