Computational modeling of stem and progenitor cell kinetics identifies plausible hematopoietic lineage hierarchies
Lisa Bast,
Michèle C. Buck,
Judith S. Hecker,
Robert A.J. Oostendorp,
Katharina S. Götze,
Carsten Marr
Affiliations
Lisa Bast
Helmholtz Zentrum München–German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg, Germany; Technical University of Munich, Department of Mathematics, Chair of Mathematical Modeling of Biological Systems, Garching, Germany
Michèle C. Buck
Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Department of Internal Medicine III, Munich, Germany
Judith S. Hecker
Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Department of Internal Medicine III, Munich, Germany
Robert A.J. Oostendorp
Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Department of Internal Medicine III, Munich, Germany
Katharina S. Götze
Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Department of Internal Medicine III, Munich, Germany; German Cancer Consortium (DKTK), Heidelberg, Partner Site Munich, Germany; Corresponding author
Carsten Marr
Helmholtz Zentrum München–German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg, Germany; Technical University of Munich, Department of Mathematics, Chair of Mathematical Modeling of Biological Systems, Garching, Germany; Corresponding author
Summary: Classically, hematopoietic stem cell (HSC) differentiation is assumed to occur via progenitor compartments of decreasing plasticity and increasing maturity in a specific, hierarchical manner. The classical hierarchy has been challenged in the past by alternative differentiation pathways. We abstracted experimental evidence into 10 differentiation hierarchies, each comprising 7 cell type compartments. By fitting ordinary differential equation models with realistic waiting time distributions to time-resolved data of differentiating HSCs from 10 healthy human donors, we identified plausible lineage hierarchies and rejected others. We found that, for most donors, the classical model of hematopoiesis is preferred. Surprisingly, multipotent lymphoid progenitor differentiation into granulocyte-monocyte progenitors is plausible in 90% of samples. An in silico analysis confirmed that, even for strong noise, the classical model can be identified robustly. Our computational approach infers differentiation hierarchies in a personalized fashion and can be used to gain insights into kinetic alterations of diseased hematopoiesis.