[BioC] Pathview with minor species

Luo Weijun luo_weijun at yahoo.com
Wed Jan 29 21:25:56 CET 2014


Someone asked a similar question on using pathview or gage with honey bee. Here is the solution:
http://seqanswers.com/forums/showthread.php?t=35472#6

Like honey bee, B. subtilis (bsu) is not a major species with Bioconductor gene annotation package as those listed in data(bods) under pathview package. similarly you need to map Locus_tags to Entrez gene ID manually. Again, you can download the gene_info data file (for all bacteria) from NCBI ftp site:
ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/GENE_INFO/Archaea_Bacteria/Bacteria.gene_info.gz
Column 2-6 should be (Entrez) GeneID, Symbol, LocusTag, Synonyms, dbXrefs. Then follow the route described in the honey bee case.

HTHs.

--------------------------------------------
On Wed, 1/29/14, David wrote:

 Subject: Pathview with minor species

 Date: Wednesday, January 29, 2014, 11:09 AM
 
 Dear Mr Luo,
 
 In our course we got B. subtilis (bsu) microarray data
 and I want to visualize some expression changes with
 Pathview. But I have a problem. We got Locus_tags like
 BSU00010 and no gene IDs and I don't know whether Pathview
 can use it or not. I really don't understand the part about
 it in the vignette. So I thought of converting them to
 UniProt IDs. But first I wanted to do a test run. So I
 searched for the UniProt IDs of 3 genes of the glycolysis
 from B. subtilis and put them into a data.frame with
 fictional data (dummy).
 
 library(pathview)
 dummy<-data.frame(name=c("P18159","P80860","P13243"),ratio=c(1,0,-1),row.names="name")
 dummy<-as.matrix(dummy)
 
 pv.out <- pathview(gene.data = dummy[, "ratio"],
 pathway.id = "00010",species = "bsu", out.suffix = "dummy",
 kegg.native=T, gene.idtype="UNIPROT")
 
 
 When I tried it this way, an error occured: No proper gene
[[elided Yahoo spam]]
 Afterwards I tested it with 3 genes from human and hsa as
 species and than it worked fine. Nothing I tried to solve
 the problem, worked. I have the latest version of
 Bioconductor and Pathview installed. Does Pathview not work
 with Bacillus subtilis subsp. subtilis 168 (=bsu) even
 though it's a KEGG organism? Or what am I doing wrong?
 
 I look forward to hearing from you.
 
 Sincerely yours,
 David



More information about the Bioconductor mailing list