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. In this work, a 12-gene network model describing the regulatory interactions driving differentiation into erythrocytes and neutrophils is investigated. We perform nonlinear stability analysis on this high-dimensional system to identify its qualitative behaviors in terms of stability characteristics. The results yield a deeper understanding of the stability structure underlying hematopoietic regulation and offer a foundation for the development of targeted therapeutic strategies.