This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DESeq2 workflow for RNAseq

Hi everyone,

I just made some graphics Pheatmap and MA plot using DESeq2 for my RNAseq differential expression analysis. However now I am really confused on what to do next. Can I make some table and calculate number of expressed genes across samples. If I could compare the data for my replications and treatments through barplots and venndiagrams .I tried to look on so many stuffs which is just making me more confused.

Looking forward for the suggestions and guidelines.

Thank you,

Ambika

rna-seq
• 3,154 views • link
updated by Devon Ryan Devon Ryan • written by AP • AP

1 answer

You can plot the (normalized) counts across your samples for particular genes if you'd like. However, looking for "expressed" genes that change between samples is usually not useful (there's no really useful definition of "expressed"). Normally one takes the statistical results from DESeq2 and then goes from there, with things like GO/pathway enrichment or simply looking at what the various DE genes are.

• 0 views • link
written by Devon Ryan Devon Ryan

Devon Ryan thank you for the clarification. I think I will just be looking at what DE genes are. So in order to plot the counts across my samples how can I write a table for transformed values?

Thank you, Ambika

• 0 views • link
written by AP • AP

Just use the plotCounts() function.

• 0 views • link
written by Devon Ryan Devon Ryan

That would be fine and it would be better to post questions like this as standalone posts.

• 1 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this