Marius Lange Profile picture
Jan 13, 2022 17 tweets 10 min read Read on X
Finally published in @naturemethods: CellRank for directed #SingleCell fate mapping, our take on studying cellular dynamics using RNA velocity. Fantastic collab between labs of @dana_peer @sloan_kettering and @fabian_theis @HelmholtzMunich, read at go.nature.com/3Fm17lt. 1/17
Motivation: scRNA-seq is great to study cellular heterogeneity, but its application to continuous processes is limited by the fact that cells are destroyed upon sequencing. Hence, we obtain many genes across many cells, but only through static snapshots. 2/17
Trajectory inference algorithms have been developed to reconstruct dynamics from static snapshots. They order cells according to similarity and even detect branches. But what about the direction? Often unclear, esp. in perturbed settings like regeneration, reprogramming etc. 3/17
A powerful solution is RNA velocity @GioeleLaManno which uses the ratio of spliced/unspliced counts to infer gene up/down reg. Led by @VolkerBergen, we generalized the concept with scVelo.org, however, how do to interpret the high-dim., noisy vector fields? 4/17
We developed CellRank.org for this task, which reconstructs dynamics based on RNA velocity and exp. similarity in high dimensions. Our aims: robustness to noise, stochastic formulation, highly scalable and able to capture gradual nature of fate commitment. 5/17
Algorithm: set up a Markov chain to describe cell-state changes. KNN graph construction so that transitions are restricted to nearest neighbors and compute correlation of velocity vectors with nearest neighbors -> T. Propagate velocity uncert. into transitions (not shown). 6/17
We teamed up with Bernhard Reuter to coarse-grain T by adapting GPCCA to the single-cell context. This gives us macrostates and transitions among them -> identify initial, intermediate and terminal states. Available in CellRank through pyGPCCA: bit.ly/pygpcca. 7/17
How likely is each cell to reach each terminal state? CellRank finds out by computing smooth fate probabilities. To infer putative decision driver genes, we correlate these against gene expression. 8/17
Using fate probabilities, CellRank charts smooth, trajectory specific gene expression trends in any precomputed pseudotime. Of course, you can use CellRank's detected initial state to root that pseudotime. 9/17
Application: we started with validation on @morris_lab in-vitro reprogramming data which has CellTagging lineage barcodes. Using these as ground truth, we show that CellRank successfully predicts reprogramming outcome. 10/17
Moving beyond in-vitro, we applied CellRank to pancreas dev. at E15.5 where we correctly predict initial & terminal states, comp. fate probs and recover exp. trends of the main drivers. Collab @bakhti_mostafa from @LickertHeiko lab. 11/17
To demonstrate robustness, we varied key method parameters and re-computed fate probabilities. For each parameter and trajectory, we correlated fate probabilities for different parameter value pairs. 12/17
Benchmarking: How do competing methods perform on the pancreas data, with and without RNA velocity information? CellRank is the only method to correctly predict all of initial & terminal states, fate probs and driver gene exp. trends. 13/17
Scalability is key to us - thanks to Michal Klein's implementation magic, CellRank runs on 100k cells in <2 min while being extremely memory efficient. 14/17
Finally, we demonstrate that CellRank generalizes beyond normal development on lung regeneration data where it predicted a novel dedifferentiation trajectory, collab @schniering_j and @Meshal_Ansari from @SchillerLab. 15/17
If this dedifferentiation trajectory exists, we would expect to find novel intermediate cell states between Goblet and Basal cells after injury but not in control mice; that's exactly what @schniering_j found in immunofluorescence stainings. 16/17
Wrap up: CellRank.org is a scalable, easy to use framework to compute directed cell-state trajectories. Massive thanks to co-authors I couldn't mention above, @VolkerBergen & @SettyM, and thanks for support from @BayerFoundation, @jherzstiftung, @MunichQbm. 17/17

• • •

Missing some Tweet in this thread? You can try to force a refresh
 

Keep Current with Marius Lange

Marius Lange Profile picture

Stay in touch and get notified when new unrolls are available from this author!

Read all threads

This Thread may be Removed Anytime!

PDF

Twitter may remove this content at anytime! Save it as PDF for later use!

Try unrolling a thread yourself!

how to unroll video
  1. Follow @ThreadReaderApp to mention us!

  2. From a Twitter thread mention us with a keyword "unroll"
@threadreaderapp unroll

Practice here first or read more on our help page!

More from @MariusLange8

Apr 17, 2023
How can we infer precise differentiation trajectories for complex biological systems? We're addressing this question with moslin, our new algorithm combining state & lineage information to link cells across time @mor_nitzan @fabian_theis 🧵 1/18 --> biorxiv.org/content/10.110… Image
This work grew out of my 2021 lab exchange with the Nitzan lab @mor_nitzan . I teamed up with these fantastic people to co-lead the project: @zoe_piran from Mor's lab, Michal Klein from Apple and Bastiaan Spanjaard from Jan Philipp Junkers lab @MDC_Berlin @BIMSB_MDC 2/18
Cells can be linked across time-points using gene expression (GEX) similarity (WOT @geoffschieb); however, this does not necessarily limply lineage relationships. To make the mapping more reliable, we sought to incorporate lineage information from single-cell lineage tracing 3/18 Image
Read 19 tweets

Did Thread Reader help you today?

Support us! We are indie developers!


This site is made by just two indie developers on a laptop doing marketing, support and development! Read more about the story.

Become a Premium Member ($3/month or $30/year) and get exclusive features!

Become Premium

Don't want to be a Premium member but still want to support us?

Make a small donation by buying us coffee ($5) or help with server cost ($10)

Donate via Paypal

Or Donate anonymously using crypto!

Ethereum

0xfe58350B80634f60Fa6Dc149a72b4DFbc17D341E copy

Bitcoin

3ATGMxNzCUFzxpMCHL5sWSt4DVtS8UqXpi copy

Thank you for your support!

Follow Us!

:(