Question: Enrichment analysis of Differentially expressed genes
gravatar for Ambika
2.6 years ago by
United States
Ambika30 wrote:

Hello everyone,

I used STAR-HTseq-DESeq2 pipleline to analyse my RNAseq data and now I have list of differentially expressed genes. I want to do the enrichment analysis of the list of genes obtained from DEseq2. I tried to use DAVID but it didn't work for me since I am working with non-model organism (fungi). Could you please suggest any other programs that are good for enrichment analysis with non model organisms.

Thank you, Ambika

enrichment analysis rna-seq • 1.0k views
ADD COMMENTlink written 2.6 years ago by Ambika30

Is there a published genome for your organism? Is it annotated?

edit: just as a side note, HTSeq is not needed here, because STAR can output counts with the parameter --quantMode GeneCounts.

ADD REPLYlink modified 2.6 years ago • written 2.6 years ago by h.mon31k

No it is not published and not annotated

ADD REPLYlink written 2.6 years ago by Ambika30
gravatar for h.mon
2.6 years ago by
h.mon31k wrote:

Are you still working with the same fungi as here: Extract exon Intron reads from RNA seq data ? You have an assembled draft genome, a gtf with predicted genes, and mapped your RNAseq to this draft genome? Assuming so, you will have to annotate this genome somehow. A quick and dirty annotation would just be blasting / diamonding the Augustus predicted genes against UniProt, to associate each gene with a UniProt identifier. Then you could use DAVID.

Another quick and dirty option would be to map your RNAseq to a close, annotated genome - in case one exists. Tthis way you can leverage the existing annotation to perform enrichment analysis.

ADD COMMENTlink written 2.6 years ago by h.mon31k

Yes you are right I have an assembled draft genome. Right now I am trying to blast my draft genome using Blast2Go and looks like it will take forever. I will try the other idea you said using UniProt. Thank you h.mon for helping.

ADD REPLYlink written 2.6 years ago by Ambika30

Hi, h.mon, I was looking on how can I associate each identified gene with UniProt identifier, but I am lost. Can you suggest a way I can do that.

ADD REPLYlink written 2.6 years ago by Ambika30
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1408 users visited in the last hour