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
1 answer
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.
Log in to answer this question.
-
Curated model organism gene functional annotations
written by ian.will •
-
GO and KEGG enrichment analysis for non-model organism
written by Shakunthala Natarajan •
-
GO term annotation for Callosobruchus
written by sherafzalk769 •
-
GO Enrichment Pichia pastoris
written by Adishi •
-
gene ID conversion
written by Adishi •
-
How to combine p-values of differentially expressed genes from different timepoints
written by AP •
-
Gene enrichment analysis
written by AP •
-
Extract exon Intron reads from RNA seq data
written by AP •
-
GO enrichment analysis. Non-model organism.
written by Frogs •
-
GO enrichment in non-model organism
written by neethav •
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.No it is not published and not annotated