This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Downstream analysis of bam files from scRNA-Seq Data

Hi,

I am new to single cell RNA-Seq, and am trying to find a workflow related question. I have aligned about 1500 cells from two different experiments (but using same experimental procedure). I am wondering how do I go ahead with differential expression analysis. Should I merge all the bam files and then generate a count table? or should I use string tie on each cell (use the hisat-stringtie pipeline). I am a bit confused because the data set is too large. I understand that ultimately I should have an expression matrix (Genes/Expression level/cell). What should be the next step after alignment (i,e creation of individual bam files) ? Please advice.

Thanks

rna-seq single cell rna-seq next-gen
• 2,985 views • link
updated by Devon Ryan Devon Ryan • written by rajeev.vikram • rajeev.vikram

1 answer

Normally one would use featureCounts or something to get the per-gene/per-cell counts. In single-cell studies this tends to involve using UMIs for deduplication, so be sure to account for that if possible. If you happen to know what groups the various cells belong to then you can go ahead with differential expression then. Otherwise you'll need to come up with your own groups (via tSNE or what not) and then do the DE between those groups.

• 0 views • link
written by Devon Ryan Devon Ryan

Thanks Ryan, This makes things a lot more clear.

• 0 views • link
written by rajeev.vikram • rajeev.vikram

Log in to answer this question.

More posts like this