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.

Unsupervised machine learning reveals temporal components of gene expression in HeLa cells following release from cell cycle arrest

Label-free cell cycle analysis for high-throughput imaging flow cytometry

Cell cycle time series gene expression data encoded as cyclic attractors in Hopfield systems

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

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

Cell Cycle and Cell Size Dependent Gene Expression Reveals Distinct Subpopulations at Single-Cell Level

Discovering Statistically Significant Periodic Gene Expression

Reconstructing cell cycle pseudo time-series via single-cell transcriptome data

Discovery of Time-Delayed Gene Regulatory Networks Based on Temporal Gene Expression Profiling

Characterizing Gene Expressions Based on Their Temporal Observations

A phase synchronization clustering algorithm for identifying interesting groups of genes from cell cycle expression data

Supervised learning of short and high-dimensional temporal sequences for life science measurements

A Computational Framework for Prediction and Analysis of Cancer Signaling Dynamics from RNA Sequencing Data—Application to the ErbB Receptor Signaling Pathway

Single-cell Analysis of Transcription Kinetics Across the Cell Cycle

Time series-free rhythm profiling using COFE reveals multi-omic circadian rhythms in human cancers

tigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities

Modeling gene expression cascades during cell state transitions

Transient power-law behaviour following induction distinguishes between competing models of stochastic gene expression

Identification of Human Cell Cycle Phase Markers Based on Single-Cell RNA-Seq Data by Using Machine Learning Methods

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

Transcriptional Landscape Of The Human Cell Cycle