paga icon indicating copy to clipboard operation
paga copied to clipboard

RNA velocity and fig3 in the Genome Bio paper

Open DHGK opened this issue 5 years ago • 4 comments

Hi PAGA community,

In the PAGA genome bio paper fig3 c,d, it shows how users could integrate RNA velocity information with PAGA to generated a directed PAGA graph. I really think that this is an amazing feature and wish to try out!

However, when going through fig3 example notebook, I feel the example stops right after fig3 b without showing vignettes for c,d. I'm wondering if anyone has tried to incorporate RNA velocity info and would like to share how you did it?

Thanks!

DHGK avatar Jun 16 '19 20:06 DHGK

Getting the figure based on velocyto involved quite a bit of book keeping of objects etc.

These days, if you use @VolkerBergen's scvelo, it should be much easier to get results like this. @VolkerBergen adapted the code to match his package: Do you have an example that we could link to also from here, so that people can find it? The Dentate Gyrus, for instance?

falexwolf avatar Jun 17 '19 14:06 falexwolf

@DHGK you can follow this example to obtain velocities, and compute directed PAGA via

scv.tl.velocity_graph(adata)

sc.tl.paga(adata)
sc.tl.paga(adata, use_rna_velocity=True)

sc.pl.paga(adata, transitions='transitions_confidence')

VolkerBergen avatar Jun 17 '19 21:06 VolkerBergen

@falexwolf and @VolkerBergen ,

Thank you so much for your suggestions. I'll try out the scvelo package to see what I get.

DHGK avatar Jun 19 '19 00:06 DHGK

An example of applying directed PAGA can be found at https://scvelo.readthedocs.io/VelocityBasics.html.

VolkerBergen avatar Jun 21 '20 12:06 VolkerBergen