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.

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

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

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

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

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

Inference of Spatiotemporal Effects on Cellular State Transitions from Time-Lapse Microscopy

Pseudocell Tracer-A method for inferring dynamic trajectories using scRNAseq and its application to B cells undergoing immunoglobulin class switch recombination

Resolving cellular dynamics using single-cell temporal transcriptomics

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

A statistical framework for differential pseudotime analysis with multiple single-cell RNA-seq samples

Modeling single cell trajectory using forward-backward stochastic differential equations

T Cell Development: Old Tales Retold By Single-Cell RNA Sequencing.

Trajectory inference from single-cell genomics data with a process time model

Modeling gene expression cascades during cell state transitions

Harnessing spatial transcriptomics to study fate decisions and lineage specification in T cell development

Single-Cell Trajectory Detection Uncovers Progression and Regulatory Coordination in Human B Cell Development

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

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

Oscillatory differentiation dynamics fundamentally restricts the resolution of pseudotime reconstruction algorithms

SCOUT: A New Algorithm for the Inference of Pseudo-Time Trajectory Using Single-Cell Data