Functional enrichment analysis for unique gene IDs
1
0
Entering edit mode
26 days ago
Pegasus ▴ 100

Hey everyone,

I've been working on RNA-seq analysis and ended up with a list of differentially expressed (DE) genes identified using edgeR. However, out of 100 genes, I could only find gene symbols for 20 of them using various tools.

The remaining genes are identified by unique locus_tag IDs, which couldn't be converted to other symbols. Despite this, all genes are associated with GO terms, using pannzer2.

I'm now looking for a method or tool to summarize and perform functional enrichment analysis, including identifying biological pathways and Gene Ontology (GO) enrichment.

Is there a tool that can do this based on either the gene sequences or their linked GO terms, since most of them like goseq, david, goprofile etc requires known gene symboles?

Thanks in advance for any suggestions!

RNA-SEQ • 437 views
ADD COMMENT
1
Entering edit mode

If these genes do not have a symbol then it is utterly unlikely that anyone has looked at their function. I would get their Ensembl IDs and quickly paste it into gprofiler2. Probably nothing will come out.

ADD REPLY
0
Entering edit mode

Thanks for your answer, but as far as I know, gProfiler2 does not support bacteria. Could you please suggest alternative websites or software, where Ensembleare accepted. thanks

ADD REPLY
0
Entering edit mode

Since the locus-tags are unclassified, I used the CDS instead of the locus-tags to retrive their ensembl IDs using ensembl bacteria website/datatbse. The website allows to blast only 30 sequeneces at one time, and the task is to select the best corresponding/score ID for each gene/sequence.

Is there any alterantive approach that can I find the best score ensembl ID for all genes/CDs in the genome, instead of going through them 30 by 30 ?

ATpoint

Thanks in advance for any help

ADD REPLY
1
Entering edit mode
25 days ago

If you are looking at producing GO annotations, we currently only recommend InterProScan, and basically only for novel organisms or where there are no other options. See the GO annotation download guide for full options.

InterProScan does not currently provide an R package, but it should be trivial to call InterProScan from R by passing the command to the "system()" function.

ADD COMMENT

Login before adding your answer.

Traffic: 1167 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6