Dynamo: Mapping Vector Field of Single Cells
Inclusive model of expression dynamics with metabolic labeling based scRNA-seq / multiomics, vector field reconstruction, potential landscape mapping, differential geometry analyses, and most probably paths / in silico perturbation predictions.
Installation - Ten minutes to dynamo - Tutorials - API - Citation - Theory
Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo, which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo’s power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.
Highlights of dynamo
- Robust and accurate estimation of RNA velocities for regular scRNA-seq datasets:
Three methods for the velocity estimations (including the new negative binomial distribution based approach)
Improved kernels for transition matrix calculation and velocity projection
Strategies to correct RNA velocity vectors (when your RNA velocity direction is problematic)
- Inclusive modeling of time-resolved metabolic labeling based scRNA-seq:
Explicitly model RNA metabolic labeling, in conjunction with RNA bursting, transcription, splicing and degradation
Comprehensive RNA kinetic rate estimation for one-shot, pulse, chase and mixture metabolic labeling experiments
- Move beyond RNA velocity to continuous vector field function for functional and predictive analyses of cell fate transitions:
Dynamical systems approaches to identify stable cell types (fixed points), boundaries of cell states (separatrices), etc
Calculate RNA acceleration (reveals early drivers), curvature (reveals master regulators of fate decision points), divergence (stability of cell states) and RNA Jacobian (cell-state dependent regulatory networks)
Various downstream differential geometry analyses to rank critical regulators/effectors, and visualize regulatory networks at key fate decision points
- Non-trivial vector field prediction of cell fate transitions:
Least action path approach to predict the optimal paths and transcriptomic factors of cell fate reprogrammings
In silico perturbation to predict the gene-wise perturbation effects and cell fate diversion after genetic perturbations
News
5/30/2023: dynamo 1.3.0 released!
3/1/2023: We welcome @Sichao25 to join the dynamo develop team!
1/28/2023: We welcome @Ukyeon to join the dynamo develop team!
15/12/2022: Thanks for @elfofmaxwell and @MukundhMurthy’s contribution. dynamo 1.2.0 released
11/11/2022: the continuing development of dynamo and the Aristotle ecosystem will be supported by CZI. See here
4/14/2022: dynamo 1.1.0 released!
3/14/2022: Since today dynamo has its own logo! Here the arrow represents the RNA velocity vector field, while the helix the RNA molecule and the colored dots RNA metabolic labels (4sU labeling). See readthedocs
2/15/2022: primers and tutorials on least action paths and in silico perturbation are released.
2/1/2022: after 3.5+ years of perseverance, our dynamo paper is finally online in Cell today!
Discussion
Please use github issue tracker to report coding related issues of dynamo. For community discussion of novel usage cases, analysis tips and biological interpretations of dynamo, please join our public slack workspace: dynamo-discussion (Only a working email address is required from the slack side).
Contribution
If you want to contribute to the development of dynamo, please check out CONTRIBUTION instruction: Contribution