Question: Importing eggnog annotations into topGO or similar?
gravatar for rmash
8 months ago by
rmash10 wrote:

I'm working with a non-model organism. I took some genes of interest and mapped it to eggnog. I got an annotation file that has a column with sometimes multiple GO terms per row.

What is the best way to look at GO term enrichment?

I can import the relevant columns (ID) and column with GO terms into R.

But do I still need a reference database now that these are mapped?

What's the best way to look like GO term enrichment?

This is what happens when I try to import my custom annation into topGO, where one column are some IDs and the second column is of GOterms

geneID2GO <- readMappings(file = "fc_test.csv", sep="\t") 
Error in `[.data.frame`(a, , 2) : undefined columns selected
ADD COMMENTlink modified 8 months ago by Biostar ♦♦ 20 • written 8 months ago by rmash10

Without seeing the file content, one can only guess at what the problem is although most likely it is a data format issue. Check that your file conforms to what's expected by the readMappings function, i.e. columns are tab-separated and the first one contains gene IDs, the second one a comma-separated list of GO term IDs, e.g.:

068724  GO:0005488, GO:0003774, GO:0001539, GO:0006935, GO:0009288
119608  GO:0005634, GO:0030528, GO:0006355, GO:0045449, GO:0003677, GO:0007275
ADD REPLYlink written 8 months ago by Jean-Karim Heriche21k

Thanks, you were right!

This is more of an issue with EggNOG but do you find that it doesn't assign GO terms to most of the orthologs?

ADD REPLYlink written 3 months ago by rmash10
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: 891 users visited in the last hour