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
ADD COMMENT
3
Entering edit mode
ADD REPLY
1
Entering edit mode
ADD REPLY
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.

ADD REPLY
1
Entering edit mode
5 weeks ago
Badran ▴ 50

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).

ADD COMMENT

Login before adding your answer.

Traffic: 1912 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6