A mathematical framework for human neutrophil state transitions inferred from single-cell RNA sequence data

Neutrophils, the most abundant immune cells in the human circulation, play a central role in the innate immune system. While neutrophil heterogeneity is a topic of increasing research interest, few efforts have been made to model the dynamics of neutrophil population subsets. We develop a mathematical model to describe the dynamics that characterizes the states and transitions involved in the maturation of human neutrophils. We use single-cell gene expression data to identify five clusters of healthy human neutrophils, and pseudo-time analysis to inform model structure. We find that precursor neutrophils transition into immature neutrophils, which then either transition to an interferon-responsive state or continue to mature through two further states. The key model parameters are the transition rates (the inverse of a transition rate is the mean waiting time in one state before transitioning to another). In this framework, the transition from the precursor to immature state (mean time less than an hour) is more rapid than subsequent transitions (mean times more than 12 hours). Approximately a quarter of neutrophils are estimated to follow the interferon-responsive path; the remainder continue along the standard maturation pathway. We use Bayesian inference to describe the variation, between individuals, in the fraction of cells within each cluster.