r/bioinformatics • • 16h ago

science question Is Metascape reliable for non-mammalian model organisms like C. elegans?

Post image
6 Upvotes

I ask because when I used it for my dataset from C. elegans, the homologene gene symbols and descriptions are obviously wrong. A lot of C. elegans genes do not have names, and instead use sequence names based of the BAC and YAC positions they were cloned onto. Examples are C30F2.4 and R08F11.4. When I ran such genes through analysis, the homolog instead showed completely different genes: C30F2.4 - C30F2.6, and R08F11.4 - fpn-1.3 (which has the sequence name R08F11.6). What's more, some genes with proper names that are well-annotated, like zipt-2.3, wasn't recognized at all. zipt-2.3 is a member of ZIP family, a family of zinc transporters, but the line for that gene is blank. This resulted in problems in enrichment analyses downstream.

I initially used WormEnrichr for my dataset, but the libraries were outdated (GO term 2018) and redundant GO terms weren't grouped. This made it very hard to make figures. Other tools like WormCat didn't have my zinc transporters annotated, either. Enrichment tool from Wormbase also trimmed the zinc transport GO terms. Under "zinc transport" on Wormbase, there are a list of genes associated with the term. If you paste this list to the Wormbase enrichment tool, you won't get the same GO term back in the results.

Right now I'm completely baffled. What am I supposed to use for my enrichment analyses, and what am I supposed to do with the results I get back from Metascape?