Abstract:Abstract Hematopoietic malignancies arise from driver alterations acquired during specific, often intermediate and transient, differentiation states. When adequately sampled, “snapshot” single-cell measurements can capture dynamic biological processes, including transient states. While there are effective computational modeling solutions to order and capture the relationships among these cell states (often projecting onto a lower dimensional manifold with directionality), methods to infer the causal gene regulatory mechanisms driving cell state transitions remain limited. Recently, a mathematical framework for modeling cell dynamics from snapshot data using a physics-based drift-diffusion equation was proposed. This drift-diffusion equation framework describes the dynamics of cell distributions with respect to time, wherein stereotypical driving forces (e.g., hematopoietic development) are captured by the drift term, and biological stochasticity is attributed to the diffusion term. The existing modeling solutions within this framework are forced to assume a fixed diffusion parameter across cell states, ascribing all observed dynamics to drift alone. Here, we expand on current drift-diffusion models to learn both the drift and diffusion terms that describe every observed cell state. These models are constructed from stochastic neural differential equations, continuous deep neural networks that approximate the theoretically complex differential equations underlying gene regulatory networks in hematopoiesis. Lineage-barcoded, multi-time point single-cell RNA sequencing experiments can approximate ground truth cellular dynamics. Using such data, we demonstrate the predictive accuracy of our proposed approach through benchmarking against the state-of-the-art methods. Through our drift-diffusion model, we can now distinguish, for the first time, between the deterministic and stochastic contributions to cell fate decision making. In doing so, we identify genes whose expression is associated with cellular drift and indeed show that they correspond to well-known markers of specific cell types; moreover, we identify genes that are associated with cell states that exhibit notable diffusion. Finally, we demonstrate the potential of such models to capture the essence of the dynamic processes through prediction on out-of-distribution data (referred to as “transfer learning” in the machine learning field). For example, when the model is trained on in vitro hematopoiesis data, our model is able to accurately predict cell trajectories in the corresponding in vivo mouse model. We anticipate that such insights will be useful in identifying key dynamic and time-dependent regulatory processes in normal development and cancer. Citation Format: Michael E. Vinyard, Anders W. Rasmussen, Ruitong Li, Luca Pinello, Gad Getz. Modeling single-cell dynamics using stochastic generative models based on neural differential equations. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5371.

Inferring Population Dynamics from Single-Cell RNA-sequencing Time Series Data.

Beyond pseudotime: Following T-cell maturation in single-cell RNAseq time series

Concepts and Limitations for Learning Developmental Trajectories from Single Cell Genomics.

Generalizing RNA Velocity to Transient Cell States Through Dynamical Modeling

Reconstructing developmental trajectories using latent dynamical systems and time-resolved transcriptomics

CellRank for Directed Single-Cell Fate Mapping

Inference of Differentiation Time for Single Cell Transcriptomes Using Cell Population Reference Data

Modeling single cell trajectory using forward-backward stochastic differential equations

Resolving cellular dynamics using single-cell temporal transcriptomics

Reconstructing gene regulatory dynamics from high-dimensional single-cell snapshot data

Decoding the Regulatory Network of Early Blood Development from Single-Cell Gene Expression Measurements.

Mapping transcriptomic vector fields of single cells

Inferring extrinsic factor-dependent single-cell transcriptome dynamics using a deep generative model

Oscillatory differentiation dynamics fundamentally restricts the resolution of pseudotime reconstruction algorithms

Distance Covariance Entropy Reveals Primed States and Bifurcation Dynamics in Single-Cell RNA-Seq Data

Gene trajectory inference for single-cell data by optimal transport metrics

Unraveling the timeline of gene expression: A pseudotemporal trajectory analysis of single-cell RNA sequencing data

Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology

Abstract 5371: Modeling single-cell dynamics using stochastic generative models based on neural differential equations

Towards a mathematical framework for modelling cell fate dynamics