How to perform pathway analysis of Candida albicans using pathview?
1
0
Entering edit mode
@danielmolecular-9513
Last seen 8.0 years ago

Hi,

I'm doing a study that analyzes RNA-seq results from mouse and Candida albicans. For each organism, I have six different conditions. For a pathway analysis I used the Pathview package with the multi-state option, and it worked pretty fine with mice data. Bellow the code used for mouse data:

>  pv.out=pathview(gene.data=tab[,1:6], gene.idtype ="ENSEMBL", pathway.id="06016" ,
              species="mmu",out.suffix="Huntington", kegg.native = T,limit=2, res=300)

 

However, when I try to analyze Candida data, I get the following error:

> pv.out=pathview(gene.data=tab[,1:6], gene.idtype ="GENENAME", pathway.id="03440",
+                 species="cal",out.suffix="MAPK", kegg.native = T,limit=2, res=300)
Error in pathview(gene.data = tab[, 1:6], gene.idtype = "GENENAME", pathway.id = "03440",  :
  No proper gene annotation package available!

 

Candida's results were annotated using Gene name as ID from the gff file from the Candida Genome database.  I am aware that KEGG database includes C. albicans data, so I wonder what the problem might be.

 

I really appreciate any kind of help.

Best regards.

Daniel

 

> sessionInfo()
R version 3.2.1 (2015-06-18)
Platform: x86_64-apple-darwin10.8.0 (64-bit)
Running under: OS X 10.7.5 (Lion)

locale:
[1] pt_BR.UTF-8/pt_BR.UTF-8/pt_BR.UTF-8/C/pt_BR.UTF-8/pt_BR.UTF-8

attached base packages:
[1] parallel  stats4    stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] org.Mm.eg.db_3.2.3   pathview_1.10.1      org.Hs.eg.db_3.2.3   RSQLite_1.0.0       
 [5] DBI_0.3.1            AnnotationDbi_1.32.0 IRanges_2.4.1        S4Vectors_0.8.2     
 [9] Biobase_2.30.0       BiocGenerics_0.16.1

loaded via a namespace (and not attached):
 [1] graph_1.48.0      KEGGgraph_1.28.0  XVector_0.10.0    magrittr_1.5      zlibbioc_1.16.0  
 [6] R6_2.1.1          stringr_1.0.0     httr_1.0.0        tools_3.2.1       grid_3.2.1       
[11] png_0.1-7         Rgraphviz_2.14.0  KEGGREST_1.10.0   stringi_1.0-1     Biostrings_2.38.0
[16] XML_3.98-1.3 
pathview • 1.4k views
ADD COMMENT
0
Entering edit mode
Luo Weijun ★ 1.6k
@luo-weijun-1783
Last seen 10 months ago
United States
Daniel, Candida albicans is not a model organism annotated in Bioconductor. In pathview process, any gene ID need to be mapped to the gene ID used in KEGG. The mapping fails without proper annotation package for Candida albicans. That’s why you get error. You can check the corresponding korg entry for Candida albicans: data(korg) korg[korg[,1]=="cal",] output is: kegg.code scientific.name common.name entrez.gnodes "cal" "Candida albicans" "" "0" kegg.geneid ncbi.geneid "CaO19.10002" "3639782" So the KEGG gene ID looks like: "CaO19.10002", and it is not Enrtez Gene ID for sure. I am not sure what type of gene ID that is, you should have better idea on this. you need to map your genename to KEGG gene ID manually, or you can obtain this ID mapping table and use mol.sum() function for the data mapping. HTHs, Weijun
ADD COMMENT

Login before adding your answer.

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