Question: Retrieving the gene list related to all GO ID
0
gravatar for Sreeraj Thamban
2.3 years ago by
Indian Institute of Science Education and Research
Sreeraj Thamban140 wrote:

Hi Biostars, I used topGO to do GO enrichment analysis for my list of interesting genes. After the run, I used the following script to retrieve significant genes for each GO ID. Can anyone help me to automate this process in such a way that I will get a table with all enriched GO IDs with corresponding significant genes? Thank you

# retrieving the gene list related to a GO ID 
# retrieve genes2GO list from the "expanded" annotation in GOdata
allGO <- genesInTerm(topGOobject)
#Getting all annotated genes for a GO ID
allGO["GO:123456"]
RetrivedGenes <- lapply(allGO,function(x) x[x %in% INT.GENES$V1] ) # where INT.GENES$V1 is my list of interesting genes
# Your significant genes for GO:123456
RetrivedGenes[["GO:123456"]]
rna-seq topgo gene ontology R • 1.5k views
ADD COMMENTlink modified 2.3 years ago by e.rempel770 • written 2.3 years ago by Sreeraj Thamban140
1
gravatar for EagleEye
2.3 years ago by
EagleEye6.4k
Sweden
EagleEye6.4k wrote:

Give a try with GeneSCF (simple to use),

Here is an example to download complete Gene Ontology as simple table with corresponding genes

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by EagleEye6.4k

Hi EagleEye, I checked GeneSCF , but we have been using topGO in our pipeline for long and repeating all the analysis using GeneSCF is not feasible at this point of time. It will be great if you can suggest any modification in the above script. Thank you

ADD REPLYlink written 2.3 years ago by Sreeraj Thamban140
1
gravatar for Chirag Parsania
2.3 years ago by
Chirag Parsania1.5k
University of Macau
Chirag Parsania1.5k wrote:

Hi,

I have automated the package for species of my interest. Use below github link for automated script of topGO. Input format provided in the sample data folder.

https://github.com/cparsania/topGO_Customized

~Chirag.

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by Chirag Parsania1.5k
1
gravatar for e.rempel
2.3 years ago by
e.rempel770
Germany, Heidelberg, COS
e.rempel770 wrote:

Hi,

I encountered the same issue and solved it this way:

results <- runTest(topGOobject, algorithm = "elim", statistic = "fisher") 
results.tab <- GenTable(object = topGOobject, elimFisher = results)  # results.tab is the table with enriched GO
AnnotatedGenes = lapply(results.tab$GO.ID, function(x) as.character(unlist(genesInTerm(object = topGOobject, whichGO = x)))) # list containg genes annotated to significant GO terms
SignificantGenes = lapply(AnnotatedGenes, function(x) intersect(x, INT.GENES$V1)) # where INT.GENES$V1 is your list of interesting genes

HTH

ADD COMMENTlink written 2.3 years ago by e.rempel770
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 961 users visited in the last hour