“…Theoretically, cell-fate transitions and switch-like behavior has been modeled in terms of gene-regulatory networks (GRNs) and an associated system of ordinary differential equations (ODEs) that describe the dynamic changes in TF concentrations. 4 , 7 , 8 , 9 In most cases however, the full underlying GRNs are unknown, and the associated ODEs typically include many unknown parameters whose precise values may strongly affect the resulting dynamics, rendering this ODE-modeling approach impractical. In the face of these difficulties, single-cell RNA-Seq (scRNA-Seq) data 10 and other single-cell omic data types 11 offer the unprecedented opportunity to study the underlying bifurcation dynamics from a more data-driven perspective, as exemplified by the development of numerous lineage-trajectory inference algorithms 12 , 13 and statistical methods to detect TFs controlling the cell-fate commitment process.…”