This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to find differentially expressed genes after mapping with Stampy?

Hi! I'm new in using Stampy and DESeq2. The problem is I have mapped RNA-seq data to the reference transcriptome and got sam files. I have treated and untreated samples. Now I want to find differentially expressed genes. I want to use DESeq2 for it. But I cannot figure out how to import my files. Could somebody give some hints?

stampy deseq2

2 answers

Take a look at these links and get a hang of it. You can also see this. Another thing is you have to convert sam to bam and for each samples use featurecounts to extract the raw read counts and then take this matrix of rew read counts for all the sample and feed them in DESeq2 and run the differential expression analysis ( obviously you have to make your custom design and also choose the normalization method that best fits for your data).

sam ->bam-> featureCounts on all bam files->obtain countmatrix data file -> DESeq2 -> Perform DEA

Thank you very much! That was very useful!

Hi katerina192,

When you have sam file, firstly you should get the count data. For a beginner, please make use of Trinity pipeline which has in built scripts to count the data from sam/bam files and find the diff expressed genes using DESeq2 and other tools. It has extensive manual which is quite easy to follow as a beginner. there is mailing list where you could post and get help immediately.

https://github.com/trinityrnaseq/trinityrnaseq/wiki/Post-Transcriptome-Assembly-Downstream-Analyses

https://github.com/trinityrnaseq/trinityrnaseq/wiki/Trinity-Transcript-Quantification

https://github.com/trinityrnaseq/trinityrnaseq/wiki/Trinity-Differential-Expression

Log in to answer this question.