scVelo on Seurat umap
1
1
Entering edit mode
17 months ago
V ▴ 320

Hello,

Using scVelo, after you've computed velocities etc and it gives you your velocity plot. Is there a way to overlay the velocity plots onto the embedding of the UMAP you obtain from your Seurat analysis?

Thanks

Seurat scVelo RNA • 3.4k views
3
Entering edit mode
1
Entering edit mode
0
Entering edit mode

Thank you. Is the best way to output from Seurat the ‘asLoom’ command? I’ve tried various ways and getting different errors that’s why I resulted to re-analysing from scratch and then wanted to embed it.

1
Entering edit mode
5 weeks ago

There is also this guide from Seurat, which shows you how to do the pre-processing in R and velocity analysis in python using Jupyter notebooks (also possible using Rstudio).

The key point here is that you convert the Seurat object to an h5seurat object, which is then converted to an h5ad object; this h5ad object is what is then loaded into python to project RNA velocity on Seurat's UMAP embeddings. DON'T try to convert a Seurat object using as.loom() because the function will inevitably fail since the Seurat object does not contain the necessary information to create the loom file (the loom file is generated with information about spliced and unspliced counts, which the Seurat object does not normally contain).