[Bioperl-l] Extracting Gene Names Genome Ontology (GO) with Perl

[Chris Mungall]

You could perform the closure locally and then iterate over the  
individual IDs or construct a big disjunctive query to Entrez -  
either way it's not so efficient, especially for less specific nodes  
(distributed queries with ontologies is an interesting challenge).

Soon you'll be able to do the same query over the GO Database / AmiGO  
using a REST API

On Apr 16, 2007, at 12:25 PM, Chris Fields wrote:

> You are correct; it explains why the list is only 120 genes.  The
> only way (currently) to do so would be to perform the closure locally
> somehow (maybe via go-perl or similar).
>
> chris
>
> On Apr 16, 2007, at 1:41 PM, Chris Mungall wrote:
>
>>
>> Unless the Entrez interface has changed since I last looked, the
>> query below for "pyrimidine ribonucleotide biosynthetic process"
>> will NOT perform the transitive closure over the graph; this means
>> genes and gene products annotated to GO:0009174 "pyrimidine
>> ribonucleoside monophosphate biosynthetic process", for example
>>
>> On Apr 16, 2007, at 9:25 AM, Chris Fields wrote:
>>
>>> You can limit EntrezGene searches by Gene Ontology ID using the  
>>> [Gene
>>> Ontology] field in queries.  The following query:
>>>
>>> '9220[Gene Ontology]'
>>>
>>> will give 120 gene IDs.  You can get the same list using the still-
>>> under-development Bio::DB::EUtilities (usual EUtilities caveat: I'm
>>> still working on this):
>>>
>>> my $esearch = Bio::DB::EUtilities->new(-eutil => 'esearch',
>>>                                         -db => 'gene',
>>>                                         -term => '9220[Gene
>>> Ontology]',
>>>                                         -retmax => 300);
>>> $esearch->get_response;
>>> my @ids = $esearch->get_ids;
>>> print join "\n", at ids;
>>>
>>> In my opinion, Sean's idea of using SQL is probably better if you
>>> have tons of searches to do.
>>>
>>> chris
>>>
>>> On Apr 16, 2007, at 9:36 AM, Wijaya Edward wrote:
>>>
>>>>
>>>> Dear all,
>>>>
>>>> Given a GO id, is there a way to extract all
>>>> the related gene names from that id with Perl?
>>>>
>>>> Anybody has experience with that?
>>>> I've looked through GO module in CPAN, but can't seem
>>>> to find any tool that facilitated that searc
>>>>
>>>> Look forward very much for your advice.
>>>>
>>>> --
>>>> Edward WIJAYA
>>>> SINGAPORE
>>>>
>>>> ------------ Institute For Infocomm Research - Disclaimer
>>>> -------------
>>>> This email is confidential and may be privileged.  If you are not
>>>> the intended recipient, please delete it and notify us immediately.
>>>> Please do not copy or use it for any purpose, or disclose its
>>>> contents to any other person. Thank you.
>>>> --------------------------------------------------------
>>>> _______________________________________________
>>>> Bioperl-l mailing list
>>>> Bioperl-l at lists.open-bio.org
>>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>
>>> Christopher Fields
>>> Postdoctoral Researcher
>>> Lab of Dr. Robert Switzer
>>> Dept of Biochemistry
>>> University of Illinois Urbana-Champaign
>>>
>>>
>>>
>>> _______________________________________________
>>> Bioperl-l mailing list
>>> Bioperl-l at lists.open-bio.org
>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>
>>
>
> Christopher Fields
> Postdoctoral Researcher
> Lab of Dr. Robert Switzer
> Dept of Biochemistry
> University of Illinois Urbana-Champaign
>
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>