Question: PFAM enrichment analysis of RNAseq data from Non-model organism
gravatar for Yunqing Yu
15 months ago by
Yunqing Yu0
Yunqing Yu0 wrote:

I am looking for a software to do pfam enrichment analysis for my differential expressed genes from RNAseq. I have the differentially expressed gene list, and the annotaiton (GO, pfam, KEGG) of the genes in the genome. I did the GO enrichment analysis using GOseq. However, Goseq cannot do pfam. I wonder if there is a similar software for any annotation term that I want to analyze, particularly pfam.

The species I am working on is not model species. So I am looking for a software that allows me to upload my own annotation file.

Thank you very much.

ADD COMMENTlink modified 12 months ago by Michael Dondrup46k • written 15 months ago by Yunqing Yu0

This might be a bit bold as an answer, but AFAIK and as unlikely as it might seem, such an online-tool does not exist, please prove me wrong though! All online-tools I know work with pre-annotated genes only, it seems as supporting the minority of NOMO researchers is not a prime concern. One could use a standard F-test or Chi-square test online tools but they often offer only limited contingency tables and are not adept to the requirements genome-wide analysis.

Assuming, the term relationships in PFAM can be ignored, I would do a Chi-square test in R locally, taking the genomic distribution of PFAM terms as background distribution. Hope this works out. Or check GSEA as an alternative.

Edit: I think I got the (false?) impression that you were looking for an online tool because you wrote "upload". If any tool would do, then R f-test/or Chi-square would be valid.

ADD REPLYlink modified 12 months ago • written 12 months ago by Michael Dondrup46k

Btw, don't get me wrong, these tools are great. If one can get their hands on PFAM terms assigned to genes in a non-model organism and make gene-sets, testing for enrichment is possibly easy-breezy.

ADD REPLYlink written 12 months ago by Michael Dondrup46k
gravatar for reza.jabal
12 months ago by
New York, USA
reza.jabal330 wrote:

Have you tried Enrichr?

Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, McDermott MG, Monteiro CD, Gundersen GW, Ma'ayan A. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Research. 2016;

ADD COMMENTlink written 12 months ago by reza.jabal330

In my understanding, Enrichr does only work for model organisms. Unfortunately, most online enrichment tools are working for model organisms only. I guess for people working only with well annotated organisms, it is often difficult to imagine the hardships we are going though...

ADD REPLYlink modified 12 months ago • written 12 months ago by Michael Dondrup46k
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: 1324 users visited in the last hour