DESeq2-->GO term heat map
1
0
Entering edit mode
shaunnak • 0
@shaunnak-19108
Last seen 2.9 years ago

Hello,

I have RNA-seq data that I would like to visualize via a GOterm heat map between time points. I have successfully been able to run DESeq2 on my transcriptome, but the output file with logpvalue/pvalue is my entire transcriptome, not compartmentalized by time point (0, 6, 12, 24, 48hr). I was able to make a heat map from this via time point but again, I have A LOT of genes so wanted to group by GOterm.

I have examined GO enrichment with programs like GOMWU, but again this examines the entire transcriptome, not each time point to compare between samples. Is there a way to use DESeq2 results and merge with GO terms to create a heat map?

I'm hoping my wording for this questions isn't too disorganized. Any suggestions on packages or programs that might help with my goal? I'm fairly new to all of this so apologies in advance.

deseq2 goterm • 436 views
ADD COMMENT
0
Entering edit mode

There are choices to how you define the sample's performance in each GO term. I have some ideas for the way to do it. Try the package rgsepd which was designed for this purpose. 

ADD REPLY
0
Entering edit mode

You might also try restating your goals. Saying you want to 'visualize via a GOterm heat map between time points' doesn't really explain exactly what you are after. One could assume that you want to find GO terms, for which the underlying set of genes appears to be perturbed in your experiment, and given those GO terms might want to make either individual heatmaps (by GO term), or a largish heatmap containing the genes in the top N GO terms (maybe separated into groups?).

But if you don't say exactly what you are after, you leave people to either guess at what you want or ignore your question entirely.

If I am right about what you want, things like goseq will get you halfway there (e.g., give you the GO terms that appear to be perturbed), and then you can make heatmaps using the genes in each term you care about.

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 4 days ago
United States

DESeq2 doesn't have any gene set visualization features, maybe try posting to Biostars?

ADD COMMENT

Login before adding your answer.

Traffic: 269 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