An introduction to biodbKegg
30 April 2024
Package
biodbKegg 1.10.0
1 Introduction
biodbKegg is a biodb extension package that implements a connector to KEGG Compound database (Kanehisa and Goto 2000).
2 Installation
Install using Bioconductor:
if (!requireNamespace("BiocManager", quietly=TRUE))
install.packages("BiocManager")
BiocManager::install('biodbKegg')3 Initialization
The first step in using biodbKegg, is to create an instance of the biodb
class BiodbMain from the main biodb package. This is done by calling the
constructor of the class:
mybiodb <- biodb::newInst()During this step the configuration is set up, the cache system is initialized and extension packages are loaded.
We will see at the end of this vignette that the biodb instance needs to be
terminated with a call to the terminate() method.
4 Creating a connector to KEGG Compound database
In biodb the connection to a database is handled by a connector instance that you can get from the factory. biodbKegg implements a connector to a remote database. Here is the code to instantiate a connector:
kegg.comp.conn <- mybiodb$getFactory()$createConn('kegg.compound')## Loading required package: biodbKegg5 Accessing entries
To retrieve entries, use:
entries <- kegg.comp.conn$getEntry(c('C00133', 'C00751'))
entries## [[1]]
## Biodb KEGG Compound entry instance C00133.
##
## [[2]]
## Biodb KEGG Compound entry instance C00751.To convert a list of entries into a dataframe, run:
x <- mybiodb$entriesToDataframe(entries, compute=FALSE)
x## accession monoisotopic.mass formula molecular.mass
## 1 C00133 89.0477 C3H7NO2 89.0932
## 2 C00751 410.3913 C30H50 410.7180
## name cas.id ncbi.pubchem.comp.id
## 1 D-Alanine;D-2-Aminopropionic acid;D-Ala 338-69-2 3433
## 2 Squalene;Spinacene;Supraene 111-02-4 4013
## chebi.id
## 1 15570
## 2 15440
## kegg.reaction.id
## 1 R00399;R00401;R01147;R01148;R01149;R01150;R01225;R01344;R02718;R04369;R04611;R05861;R07651;R08850;R09588;R09595;R11965;R12557;R12812;R12863;R12867;R12871;R12873;R12875;R12904
## 2 R02872;R02874;R02875;R02876;R06223;R07322;R07323;R08535;R09712;R10167;R10169;R11401;R12355
## kegg.enzyme.id
## 1 1.4.3.3;1.4.3.19;2.1.2.7;2.3.1.263;2.3.2.14;2.6.1.21;3.1.1.103;3.4.13.22;3.4.17.8;5.1.1.1;6.1.1.13;6.1.2.1;6.3.2.4;6.3.2.16;6.3.2.35
## 2 1.3.1.96;1.14.14.17;1.14.19.-;1.17.8.1;2.1.1.262;2.5.1.21;4.2.1.123;4.2.1.129;5.4.99.17;5.4.99.37
## kegg.pathway.id
## 1 map00470;map00550;map00552;map01100;map01502;map01503;map04742
## 2 map00100;map00909;map00996;map00999;map01060;map01062;map01066;map01070;map01100;map01110
## kegg.compound.id lipidmaps.structure.id
## 1 C00133 <NA>
## 2 C00751 LMPR01060100026 Search for compounds of a certain mass
ids <- kegg.comp.conn$searchForEntries(list(monoisotopic.mass=list(value=64, delta=2.0)), max.results=10)
entries <- mybiodb$getFactory()$getEntry('kegg.compound', ids)7 Add information to a data frame containing KEGG Compound IDs
If you have a data frame containing a column with KEGG Compound IDs, you can add information such as associated KEGG Enzymes, associated KEGG Pathways and KEGG Modules to your data frame, for a specific organism.
For the example we use the list of compound IDs we already have, to construct a data frame:
kegg.comp.ids <- c('C06144', 'C06178', 'C02659')
mydf <- data.frame(kegg.ids=kegg.comp.ids)Using the addInfo() method of KeggCompoundConn class, we add information
about pathways, enzymes and modules for these compounds:
kegg.comp.conn$addInfo(mydf, id.col='kegg.ids', org='mmu')## kegg.ids kegg.enzyme.id
## 1 C06144 4.2.1.27
## 2 C06178 1.4.3.21|2.3.1.74
## 3 C02659 1.14.14.41|2.4.1.63|3.2.1.21
## kegg.reaction.id
## 1 R01611;R01367
## 2 R01853;R02382;R02529|R01613;R07987;R07988
## 3 R10030;R10034;R11597|R03625;R04948;R10037|R00026;R00306;R02558
## kegg.pathway.id
## 1 mmu00650|mmu01100
## 2 mmu00760|mmu01100
## 3 mmu01100
## kegg.pathway.name
## 1 Butanoate metabolism - Mus musculus (house mouse)|Metabolic pathways - Mus musculus (house mouse)
## 2 Nicotinate and nicotinamide metabolism - Mus musculus (house mouse)|Metabolic pathways - Mus musculus (house mouse)
## 3 Metabolic pathways - Mus musculus (house mouse)
## kegg.pathway.pathway.class kegg.module.id
## 1 Metabolism;Carbohydrate metabolism|NA M00027|M00001|M00002
## 2 Metabolism;Metabolism of cofactors and vitamins|NA M00912|M00001|M00002
## 3 <NA> M00001|M00002|M00003
## kegg.module.name
## 1 GABA (gamma-Aminobutyrate) shunt|Glycolysis (Embden-Meyerhof pathway), glucose => pyruvate|Glycolysis, core module involving three-carbon compounds
## 2 NAD biosynthesis, tryptophan => quinolinate => NAD|Glycolysis (Embden-Meyerhof pathway), glucose => pyruvate|Glycolysis, core module involving three-carbon compounds
## 3 Glycolysis (Embden-Meyerhof pathway), glucose => pyruvate|Glycolysis, core module involving three-carbon compounds|Gluconeogenesis, oxaloacetate => fructose-6PNote that, by default, the number of values for each field is limited to 3.
Please see the help page of KeggCompoundConn for more information about
addInfo(), and a description of all parameters.
The list of organisms is available at https://www.genome.jp/kegg/catalog/org_list.html.
8 Getting pathways related to compounds using KEGG databases
In this example we will look for pathways related to specified compounds, count to how many pathways each compound is related, build a pathway graph, and create a decorated pathway graph picture.
For that we will start from a given list of KEGG compound IDs, and explore KEGG to find to which organisms they can be related and try to discover links between them through KEGG pathways.
As an example, we will start from a predefined list of KEGG compound IDs, and focus on one organism, the mouse.
8.3 Build a pathway graph
To build a pathway graph, we need a connector to the KEGG Pathway database:
kegg.path.conn <- mybiodb$getFactory()$getConn('kegg.pathway')Building list of edges and vertices for pathways is done by calling buildPathwayGraph():
kegg.path.conn$buildPathwayGraph(pathways[[1]])## $vertices
## name type id
## 1 C00025 compound C00025
## 2 C00334 compound C00334
## 3 C00011 compound C00011
## 4 R00261 reaction R00261
## 6 C00026 compound C00026
## 7 C00232 compound C00232
## 9 R01648 reaction R01648
## 11 C00022 compound C00022
## 13 C00041 compound C00041
## 14 R10178 reaction R10178
## 16 C00003 compound C00003
## 17 C00001 compound C00001
## 18 C00042 compound C00042
## 19 C00004 compound C00004
## 20 C00080 compound C00080
## 21 R00713 reaction R00713
## 23 C00006 compound C00006
## 26 C00005 compound C00005
## 28 R00714 reaction R00714
##
## $edges
## from to
## 1 C00025 R00261
## 2 R00261 C00334
## 3 R00261 C00011
## 4 C00334 R01648
## 5 C00026 R01648
## 6 R01648 C00232
## 7 R01648 C00025
## 8 C00334 R10178
## 9 C00022 R10178
## 10 R10178 C00232
## 11 R10178 C00041
## 12 C00232 R00713
## 13 C00003 R00713
## 14 C00001 R00713
## 15 R00713 C00042
## 16 R00713 C00004
## 17 R00713 C00080
## 18 C00232 R00714
## 19 C00006 R00714
## 20 C00001 R00714
## 21 R00714 C00042
## 22 R00714 C00005
## 23 R00714 C00080The object returned is a list whose names are the pathway IDs submitted, and the values are lists containing two data frames (edges and vertices).
We can also get an igraph object for the a pathway (or a list of pathways):
ig <- kegg.path.conn$getPathwayIgraph(pathways[[1]])And we plot it:
vert <- igraph::as_data_frame(ig, 'vertices')
shapes <- vapply(vert[['type']], function(x) if (x == 'reaction') 'rectangle'
else 'circle', FUN.VALUE='', USE.NAMES=FALSE)
colors <- vapply(vert[['type']], function(x) if (x == 'reaction') 'yellow' else
'red', FUN.VALUE='', USE.NAMES=FALSE)
plot(ig, vertex.shape=shapes, vertex.color=colors, vertex.label.dist=1,
vertex.size=4, vertex.size2=4)
8.4 Create a decorated pathway graph picture
We will now use a KEGG pathway picture and highlight some of the enzymes and compounds on it.
For this, we first get the enzymes related to the compounds:
kegg.enz.ids <- mybiodb$entryIdsToSingleFieldValues(kegg.comp.ids,
db='kegg.compound',
field='kegg.enzyme.id')
kegg.enz.ids## [1] "4.2.1.27" "1.4.3.21" "2.3.1.74" "1.14.14.41" "2.4.1.63"
## [6] "3.2.1.21" "4.1.2.46" "4.1.2.47"define the colors we want to apply:
color2ids <- list(yellow=kegg.enz.ids, red=kegg.comp.ids)Then we call the method that builds the highlighted image and print it:
kegg.path.conn$getDecoratedGraphPicture(pathways[[1]],
color2ids=color2ids)
9 Terminate the biodb instance
When done with your biodb instance you have to terminate it, in order to ensure release of resources (file handles, database connection, etc):
mybiodb$terminate()## INFO [22:47:12.430] Closing BiodbMain instance...
## INFO [22:47:12.431] Connector "kegg.compound" deleted.
## INFO [22:47:12.439] Connector "kegg.enzyme" deleted.
## INFO [22:47:12.441] Connector "kegg.pathway" deleted.
## INFO [22:47:12.442] Connector "kegg.module" deleted.
## INFO [22:47:12.444] Connector "kegg.reaction" deleted.10 Session information
sessionInfo()## R version 4.4.0 beta (2024-04-15 r86425)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 22.04.4 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.19-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] biodbKegg_1.10.0 BiocStyle_2.32.0
##
## loaded via a namespace (and not attached):
## [1] rappdirs_0.3.3 sass_0.4.9 utf8_1.2.4
## [4] generics_0.1.3 stringi_1.8.3 RSQLite_2.3.6
## [7] hms_1.1.3 digest_0.6.35 magrittr_2.0.3
## [10] evaluate_0.23 bookdown_0.39 fastmap_1.1.1
## [13] blob_1.2.4 plyr_1.8.9 jsonlite_1.8.8
## [16] progress_1.2.3 DBI_1.2.2 tinytex_0.50
## [19] BiocManager_1.30.22 httr_1.4.7 fansi_1.0.6
## [22] XML_3.99-0.16.1 jquerylib_0.1.4 cli_3.6.2
## [25] rlang_1.1.3 chk_0.9.1 crayon_1.5.2
## [28] dbplyr_2.5.0 bit64_4.0.5 withr_3.0.0
## [31] cachem_1.0.8 yaml_2.3.8 tools_4.4.0
## [34] memoise_2.0.1 biodb_1.12.0 dplyr_1.1.4
## [37] filelock_1.0.3 curl_5.2.1 vctrs_0.6.5
## [40] R6_2.5.1 magick_2.8.3 BiocFileCache_2.12.0
## [43] lifecycle_1.0.4 stringr_1.5.1 bit_4.0.5
## [46] pkgconfig_2.0.3 pillar_1.9.0 bslib_0.7.0
## [49] glue_1.7.0 Rcpp_1.0.12 lgr_0.4.4
## [52] highr_0.10 xfun_0.43 tibble_3.2.1
## [55] tidyselect_1.2.1 knitr_1.46 igraph_2.0.3
## [58] htmltools_0.5.8.1 rmarkdown_2.26 compiler_4.4.0
## [61] prettyunits_1.2.0 askpass_1.2.0 openssl_2.1.2References
Kanehisa, Minoru, and Susumu Goto. 2000. “KEGG: Kyoto Encyclopedia of Genes and Genomes.” Nucleic Acids Research 28 (1): 27–30. https://doi.org/10.1093/nar/28.1.27.