HI to everybody, i'm using the "gene2pathway" package but i'm not able to obtain something right. i'm using gene2pathway.signaltrans in order to predict a gene's membership to a KEGG signaling pathway and/or pathway component via the contained InterPro domains and in my case for the gene="19713", that represent RET proto oncogene (http://www.uniprot.org/uniprot/P35546). But something is wrong as you can see down there....someone can help me? > pred.comp = gene2pathway.signaltrans("19713", organism = "hsa") Using KEGG information from SOAP service ... Mapping to signal transduction pathway components via KEGG database ... ---> Information found for 0 genes done. 1 genes to predict Retrieving information from InterPro database for organism ' hsa ' via Ensembl ... Checking attributes ... ok Checking filters ... ok Error in gene2pathway:::getInterProDomains(topredict, gene2Domains = gene2Domains, : No mapping Entrez gene ID -> InterPro found! Best regards. -- ----------------------------------------------------- Dr. Alberto Goldoni Bologna, Italy

Hi Alberto, It's not completely clear to me from you message how you got the gene id 19713, but if that is supposed to be an entrez gene, I note that it is the entrez gene for the mouse version of RET. But in your code you specified 'hsa', which is the 3 letter code from KEGG for humans. Marc On 03/29/2010 08:53 AM, Alberto Goldoni wrote: > HI to everybody, > i'm using the "gene2pathway" package but i'm not able to obtain something right. > > i'm using gene2pathway.signaltrans in order to predict a gene's > membership to a KEGG signaling pathway and/or pathway component via > the contained InterPro domains and in my case for the gene="19713", > that represent RET proto oncogene > (http://www.uniprot.org/uniprot/P35546). > But something is wrong as you can see down there....someone can help me? > > >> pred.comp = gene2pathway.signaltrans("19713", organism = "hsa") >> > Using KEGG information from SOAP service ... > Mapping to signal transduction pathway components via KEGG database ... > ---> Information found for 0 genes > done. > 1 genes to predict > Retrieving information from InterPro database for organism ' hsa ' via > Ensembl ... > Checking attributes ... ok > Checking filters ... ok > Error in gene2pathway:::getInterProDomains(topredict, gene2Domains = > gene2Domains, : > No mapping Entrez gene ID -> InterPro found! > > > Best regards. > > >

Dear Alberto the answer is below your email. >________________________________________ >Da: bioconductor-bounces at stat.math.ethz.ch [bioconductor-bounces at stat.math.ethz.ch] per conto di Alberto Goldoni [alberto.goldoni1975 at gmail.com] >Inviato: luned? 29 marzo 2010 17.53 >A: BioC >Oggetto: [BioC] error in gene2pathway results > >HI to everybody, >i'm using the "gene2pathway" package but i'm not able to obtain something right. > >i'm using gene2pathway.signaltrans in order to predict a gene's >membership to a KEGG signaling pathway and/or pathway component via >the contained InterPro domains and in my case for the gene="19713", >that represent RET proto oncogene >(http://www.uniprot.org/uniprot/P35546). >But something is wrong as you can see down there....someone can help me? > >> pred.comp = gene2pathway.signaltrans("19713", organism = "hsa") >Using KEGG information from SOAP service ... >Mapping to signal transduction pathway components via KEGG database ... >---> Information found for 0 genes >done. >1 genes to predict >Retrieving information from InterPro database for organism ' hsa ' via >Ensembl ... >Checking attributes ... ok >Checking filters ... ok >Error in gene2pathway:::getInterProDomains(topredict, gene2Domains = >gene2Domains, : > No mapping Entrez gene ID -> InterPro found! > > >Best regards. > > >-- >----------------------------------------------------- >Dr. Alberto Goldoni >Bologna, Italy > >_______________________________________________ >Bioconductor mailing list >Bioconductor at stat.math.ethz.ch >https://stat.ethz.ch/mailman/listinfo/bioconductor >Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor There is no mistake in the result you are getting. In fact, you are indicating Homo Sapiens (hsa) as the reference organism but you are using a mouse GeneID to indicate the RET gene, and of course this is resulting in no information found. If you try using 5979 (which is the human RET), this is what you get: > pred.comp = gene2pathway.signaltrans("5979", organism = "hsa") Using KEGG information from SOAP service ... Mapping to signal transduction pathway components via KEGG database ... ---> Information found for 0 genes done. 1 genes to predict Retrieving information from InterPro database for organism ' hsa ' via Ensembl ... Checking attributes ... ok Checking filters ... ok done: Information found for 1 out of 1 genes Model prediction possible for 1 genes ...........done Preparing output finished ...and as you can see there are valid results this way. I hope this will be of some help. All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) Visiting address: Experimental Vascular Pathology group, Dept of Pathology - Room5.08, Maastricht University Medical Center, P. Debeijelaan 25, 6229 HX Maastricht E-mail: m.manca at path.unimaas.nl Office telephone: +31(0)433874633 Personal mobile: +31(0)626441205 Twitter: @markomanka ********************************************************************** *********************************************** This email and any files transmitted with it are confide...{{dropped:15}}