Hi,
I have a long list of Gene Ontology(GO) term like below and I would like to put this list in a on-line free available webpage and then receive the genes which categorize in each GO term ID.
Thanks
GO:0007099 GO:0043516 GO:0045665 GO:0051651 GO:0007017 GO:0051235 GO:0007059 GO:0046777 GO:0051640 GO:0000278 GO:0051301 GO:0048285 GO:0007049
3 answers
You can do this in BioMart on Ensembl.
There's a help video on using BioMart here:
The tool is here: http://www.ensembl.org/biomart/martview/c7904f8b46d31846e072964b3ca0955a
You can just filter by GO accessions, then get gene names as attributes.
I assume that you're filtering by GO:0044702, then selecting GO term accession as an attribute?
What's happening here is that the GO terms that you get out of BioMart are the most specific ones available. GO terms are hierarchical, so GO:0044702 is a very general term, and GO:0001541 is a more specific term which is a subset of GO:0044702. So everything with the GO term GO:0001541 will pass the filter for GO:0044702, but only GO:0001541, or other more specific GO terms, such as GO:0001554 and GO:0007276, will be displayed.
Hi, since this thread is live again. Feel like I must also contribute with some updated/easy solution. (Note: Sorry, this is command-line solution)
You can retrieve complete list of GO terms with corresponding genes as plan text table format using GeneSCF's prepare_database module.
Download all GO terms and list of genes from human:
./prepare_database -db=GO_all -org=goa_human
If you like to retrieve genes only for you GO term list, grep your terms against GeneSCF downloaded table.
grep -fF PATH/YOUR_TERMS_oneTermPerLine.list PATH/geneSCF-master-source-v1.1-p2/class/lib/db/goa_human/GO_all_sym.txt > PATH/YOUR_TERMS_withGenes.table.txt
library(GO.db)
Term("GO:0006281")
GO:0006281 "DNA repair"
While your contribution to an old thread for completeness is valued, this barely constitutes an answer in its present form. This is presumably some R code based on the use of library(), but it isn't at all clear. As a minimum, please format the code (if that is what it is) with proper code highlighting (the 101010 button).
Log in to answer this question.
