KEGGSOAP: problem with get.pathways.by.genes()
1
0
Entering edit mode
@julien-meunier-3442
Last seen 10.2 years ago
Hi everyone, I am trying to convert entrez id (mouse) to KEGG pathway using the keggsoap package, and having a hard time. get.genes.by.pathway(c("path:mmu04060")) would give me the following result: [1] "mmu:100044793" "mmu:100047704" "mmu:110542" "mmu:11477" "mmu:11479" ... However, the following get.pathways.by.genes(c("mmu:11479")) returns character(0) I also tested get.pathways.by.genes with the IL2 gene (entrez id 16183) involved in several pathways, calling the function with "16183", "mmu:16183" and other variations with the same result. The vignette mentioned they might be "problem on the server side" with get.pathways.by.genes(), and I googled without succes. Any ideas to get this function to work ? Thank you for your time, Julien PS: I known I can use KEGG.db or org.Mm.eg.db, or other locale packages but for some reason I would rather use an online query.
Pathways convert KEGGSOAP Pathways convert KEGGSOAP • 1.4k views
ADD COMMENT
0
Entering edit mode
@saroj-k-mohapatra-3419
Last seen 10.2 years ago
Hi Julien: For some reason, if you use more than one gene, the function does return the pathways. > get.pathways.by.genes(c("mmu:11479")) character(0) > > get.pathways.by.genes(c("mmu:11479", "mmu:11477")) [1] "path:mmu04060" "path:mmu04350" Best, Saroj Julien Meunier wrote: > > > Hi everyone, > > I am trying to convert entrez id (mouse) to KEGG pathway > using the keggsoap package, and having a hard time. > > > get.genes.by.pathway(c("path:mmu04060")) > > would give me the following result: > > [1] "mmu:100044793" "mmu:100047704" "mmu:110542" "mmu:11477" > "mmu:11479" ... > > > > However, the following > > get.pathways.by.genes(c("mmu:11479")) > > returns character(0) > > > I also tested get.pathways.by.genes with the IL2 gene (entrez id 16183) > involved in several pathways, calling the function with "16183", > "mmu:16183" and other variations with the same result. > > The vignette mentioned they might be "problem on the server side" with > get.pathways.by.genes(), and I googled without succes. > > > Any ideas to get this function to work ? > > > Thank you for your time, > Julien > > > PS: I known I can use KEGG.db or org.Mm.eg.db, or other locale packages > but for some reason I would rather use an online query. > > _______________________________________________ > 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 >
ADD COMMENT
0
Entering edit mode
Hi Saroj, Thank you for your reply. Yet, I am still puzzled by the behavior of get.pathway.by.genes() in KEGGSOAP: I have been running this session: > sessionInfo() R version 2.9.0 (2009-04-17) i386-apple-darwin8.11.1 locale: C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] KEGG.db_2.2.11 RSQLite_0.7-1 DBI_0.2-4 [4] AnnotationDbi_1.6.0 Biobase_2.4.1 RCurl_0.94-1 [7] KEGGSOAP_1.18.0 loaded via a namespace (and not attached): [1] SSOAP_0.4-6 XML_2.3-0 and the following code: My.entrez2pathid<-function(refseqid){ require(KEGG.db) stopifnot(is.character(refseqid)) v<-mget(refseqid,KEGGEXTID2PATHID,ifnotfound=NA) return(v) } v1<-My.entrez2pathid(c("16185","11652")) would give this result > v1 $`16185` [1] "mmu04060" "mmu04630" $`11652` [1] "mmu04010" "mmu04012" "mmu04150" "mmu04210" "mmu04370" "mmu04510" [7] "mmu04530" "mmu04620" "mmu04630" "mmu04660" "mmu04662" "mmu04664" [13] "mmu04910" "mmu04920" "mmu05210" "mmu05211" "mmu05212" "mmu05213" [19] "mmu05214" "mmu05215" "mmu05218" "mmu05220" "mmu05221" "mmu05222" [25] "mmu05223" whereas My.ksoap.entrez2pathid<-function(refseqid){ require(KEGGSOAP) stopifnot(is.character(refseqid)) v<-get.pathways.by.genes(refseqid) return(v) } v2<-My.ksoap.entrez2pathid(c("mmu:16185","mmu:11652")) returns > v2 [1] "path:mmu04630" Checking on the KEGG website, the refseq id 11652 is associated with a lot of pathways (those returned in "v1"). Does anybody knows why KEGGSOAp seems to fail to retrieve all the pathways ? Many thanks, Julien Saroj K Mohapatra wrote: > Hi Julien: > > For some reason, if you use more than one gene, the function does return > the pathways. > > > get.pathways.by.genes(c("mmu:11479")) > character(0) > > > > get.pathways.by.genes(c("mmu:11479", "mmu:11477")) > [1] "path:mmu04060" "path:mmu04350" > > Best, > > Saroj > > > Julien Meunier wrote: >> >> >> Hi everyone, >> >> I am trying to convert entrez id (mouse) to KEGG pathway >> using the keggsoap package, and having a hard time. >> >> >> get.genes.by.pathway(c("path:mmu04060")) >> >> would give me the following result: >> >> [1] "mmu:100044793" "mmu:100047704" "mmu:110542" "mmu:11477" >> "mmu:11479" ... >> >> >> >> However, the following >> >> get.pathways.by.genes(c("mmu:11479")) >> >> returns character(0) >> >> >> I also tested get.pathways.by.genes with the IL2 gene (entrez id 16183) >> involved in several pathways, calling the function with "16183", >> "mmu:16183" and other variations with the same result. >> >> The vignette mentioned they might be "problem on the server side" with >> get.pathways.by.genes(), and I googled without succes. >> >> >> Any ideas to get this function to work ? >> >> >> Thank you for your time, >> Julien >> >> >> PS: I known I can use KEGG.db or org.Mm.eg.db, or other locale packages >> but for some reason I would rather use an online query. >> >> _______________________________________________ >> 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 >> >
ADD REPLY
0
Entering edit mode
Perhaps v2 returns the common pathway between both genes? Saroj Julien Meunier wrote: > > Hi Saroj, > > > Thank you for your reply. Yet, I am still puzzled by the behavior of > get.pathway.by.genes() in KEGGSOAP: > > > I have been running this session: > > > sessionInfo() > R version 2.9.0 (2009-04-17) > i386-apple-darwin8.11.1 > > locale: > C > > attached base packages: > [1] stats graphics grDevices utils datasets methods base > > other attached packages: > [1] KEGG.db_2.2.11 RSQLite_0.7-1 DBI_0.2-4 > [4] AnnotationDbi_1.6.0 Biobase_2.4.1 RCurl_0.94-1 > [7] KEGGSOAP_1.18.0 > > loaded via a namespace (and not attached): > [1] SSOAP_0.4-6 XML_2.3-0 > > > > and the following code: > > My.entrez2pathid<-function(refseqid){ > > require(KEGG.db) > stopifnot(is.character(refseqid)) > v<-mget(refseqid,KEGGEXTID2PATHID,ifnotfound=NA) > > return(v) > > } > > v1<-My.entrez2pathid(c("16185","11652")) > > > would give this result > > > > v1 > $`16185` > [1] "mmu04060" "mmu04630" > > $`11652` > [1] "mmu04010" "mmu04012" "mmu04150" "mmu04210" "mmu04370" "mmu04510" > [7] "mmu04530" "mmu04620" "mmu04630" "mmu04660" "mmu04662" "mmu04664" > [13] "mmu04910" "mmu04920" "mmu05210" "mmu05211" "mmu05212" "mmu05213" > [19] "mmu05214" "mmu05215" "mmu05218" "mmu05220" "mmu05221" "mmu05222" > [25] "mmu05223" > > > whereas > > My.ksoap.entrez2pathid<-function(refseqid){ > > require(KEGGSOAP) > stopifnot(is.character(refseqid)) > v<-get.pathways.by.genes(refseqid) > > return(v) > > } > > v2<-My.ksoap.entrez2pathid(c("mmu:16185","mmu:11652")) > > > returns > > > v2 > [1] "path:mmu04630" > > > > Checking on the KEGG website, the refseq id 11652 is associated with > a lot of pathways (those returned in "v1"). > > > Does anybody knows why KEGGSOAp seems to fail to retrieve all the > pathways ? > > > Many thanks, > Julien > > > > > > Saroj K Mohapatra wrote: >> Hi Julien: >> >> For some reason, if you use more than one gene, the function does >> return the pathways. >> >> > get.pathways.by.genes(c("mmu:11479")) >> character(0) >> > >> > get.pathways.by.genes(c("mmu:11479", "mmu:11477")) >> [1] "path:mmu04060" "path:mmu04350" >> >> Best, >> >> Saroj >> >> >> Julien Meunier wrote: >>> >>> >>> Hi everyone, >>> >>> I am trying to convert entrez id (mouse) to KEGG pathway >>> using the keggsoap package, and having a hard time. >>> >>> >>> get.genes.by.pathway(c("path:mmu04060")) >>> >>> would give me the following result: >>> >>> [1] "mmu:100044793" "mmu:100047704" "mmu:110542" "mmu:11477" >>> "mmu:11479" ... >>> >>> >>> >>> However, the following >>> >>> get.pathways.by.genes(c("mmu:11479")) >>> >>> returns character(0) >>> >>> >>> I also tested get.pathways.by.genes with the IL2 gene (entrez id 16183) >>> involved in several pathways, calling the function with "16183", >>> "mmu:16183" and other variations with the same result. >>> >>> The vignette mentioned they might be "problem on the server side" with >>> get.pathways.by.genes(), and I googled without succes. >>> >>> >>> Any ideas to get this function to work ? >>> >>> >>> Thank you for your time, >>> Julien >>> >>> >>> PS: I known I can use KEGG.db or org.Mm.eg.db, or other locale packages >>> but for some reason I would rather use an online query. >>> >>> _______________________________________________ >>> 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 >>> >> >
ADD REPLY

Login before adding your answer.

Traffic: 569 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