Load the package with the library function.
library(tidyverse)
library(ggplot2)
library(dce)
set.seed(42)
We provide access to the following topological pathway databases using graphite (Sales et al. 2012) in a processed format. This format looks as follows:
dce::df_pathway_statistics %>%
arrange(desc(node_num)) %>%
head(10) %>%
knitr::kable()
database | pathway_id | pathway_name | node_num | edge_num |
---|---|---|---|---|
reactome | R-HSA-162582 | Signaling Pathways | 2488 | 62068 |
reactome | R-HSA-1430728 | Metabolism | 2047 | 85543 |
reactome | R-HSA-392499 | Metabolism of proteins | 1894 | 52807 |
reactome | R-HSA-1643685 | Disease | 1774 | 55469 |
reactome | R-HSA-168256 | Immune System | 1771 | 58277 |
panther | P00057 | Wnt signaling pathway | 1644 | 195344 |
reactome | R-HSA-74160 | Gene expression (Transcription) | 1472 | 32493 |
reactome | R-HSA-597592 | Post-translational protein modification | 1394 | 26399 |
kegg | hsa:01100 | Metabolic pathways | 1343 | 22504 |
reactome | R-HSA-73857 | RNA Polymerase II Transcription | 1339 | 25294 |
Let’s see how many pathways each database provides:
dce::df_pathway_statistics %>%
count(database, sort = TRUE, name = "pathway_number") %>%
knitr::kable()
database | pathway_number |
---|---|
pathbank | 48685 |
smpdb | 48671 |
reactome | 2406 |
wikipathways | 640 |
kegg | 323 |
panther | 94 |
pharmgkb | 90 |
Next, we can see how the pathway sizes are distributed for each database:
dce::df_pathway_statistics %>%
ggplot(aes(x = node_num)) +
geom_histogram(bins = 30) +
facet_wrap(~ database, scales = "free") +
theme_minimal()
It is easily possible to plot pathways:
pathways <- get_pathways(
pathway_list = list(
pathbank = c("Lactose Synthesis"),
kegg = c("Fatty acid biosynthesis")
)
)
lapply(pathways, function(x) {
plot_network(
as(x$graph, "matrix"),
visualize_edge_weights = FALSE,
arrow_size = 0.02,
shadowtext = TRUE
) +
ggtitle(x$pathway_name)
})
## [[1]]
##
## [[2]]
sessionInfo()
## R version 4.3.1 (2023-06-16)
## Platform: aarch64-apple-darwin20 (64-bit)
## Running under: macOS Ventura 13.6.1
##
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/lib/libRlapack.dylib; LAPACK version 3.11.0
##
## locale:
## [1] C/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
##
## time zone: America/New_York
## tzcode source: internal
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] dce_1.10.0 graph_1.80.0
## [3] cowplot_1.1.1 lubridate_1.9.2
## [5] forcats_1.0.0 stringr_1.5.0
## [7] dplyr_1.1.2 purrr_1.0.1
## [9] readr_2.1.4 tidyr_1.3.0
## [11] tibble_3.2.1 tidyverse_2.0.0
## [13] TCGAutils_1.22.0 curatedTCGAData_1.23.5
## [15] MultiAssayExperiment_1.28.0 SummarizedExperiment_1.32.0
## [17] Biobase_2.62.0 GenomicRanges_1.54.0
## [19] GenomeInfoDb_1.38.0 IRanges_2.36.0
## [21] S4Vectors_0.40.1 BiocGenerics_0.48.0
## [23] MatrixGenerics_1.14.0 matrixStats_1.0.0
## [25] ggraph_2.1.0 ggplot2_3.4.2
## [27] BiocStyle_2.30.0
##
## loaded via a namespace (and not attached):
## [1] bitops_1.0-7 httr_1.4.6
## [3] GenomicDataCommons_1.26.0 prabclus_2.3-2
## [5] Rgraphviz_2.46.0 numDeriv_2016.8-1.1
## [7] tools_4.3.1 utf8_1.2.3
## [9] R6_2.5.1 vegan_2.6-4
## [11] mgcv_1.9-0 sn_2.1.1
## [13] permute_0.9-7 withr_2.5.0
## [15] graphite_1.48.0 prettyunits_1.1.1
## [17] gridExtra_2.3 flexclust_1.4-1
## [19] cli_3.6.1 sandwich_3.0-2
## [21] labeling_0.4.2 sass_0.4.6
## [23] diptest_0.76-0 mvtnorm_1.2-2
## [25] robustbase_0.99-0 proxy_0.4-27
## [27] Rsamtools_2.18.0 FMStable_0.1-4
## [29] Linnorm_2.26.0 plotrix_3.8-2
## [31] limma_3.58.0 RSQLite_2.3.1
## [33] generics_0.1.3 BiocIO_1.12.0
## [35] gtools_3.9.4 wesanderson_0.3.6
## [37] Matrix_1.6-0 fansi_1.0.4
## [39] logger_0.2.2 abind_1.4-5
## [41] lifecycle_1.0.3 multcomp_1.4-25
## [43] yaml_2.3.7 edgeR_4.0.0
## [45] mathjaxr_1.6-0 SparseArray_1.2.0
## [47] BiocFileCache_2.10.1 Rtsne_0.16
## [49] grid_4.3.1 blob_1.2.4
## [51] promises_1.2.0.1 gdata_2.19.0
## [53] ppcor_1.1 bdsmatrix_1.3-6
## [55] ExperimentHub_2.10.0 crayon_1.5.2
## [57] lattice_0.21-8 GenomicFeatures_1.54.0
## [59] KEGGREST_1.42.0 magick_2.7.4
## [61] pillar_1.9.0 knitr_1.43
## [63] rjson_0.2.21 fpc_2.2-10
## [65] corpcor_1.6.10 codetools_0.2-19
## [67] mutoss_0.1-13 glue_1.6.2
## [69] RcppArmadillo_0.12.4.1.0 data.table_1.14.8
## [71] vctrs_0.6.3 png_0.1-8
## [73] Rdpack_2.4 mnem_1.18.0
## [75] gtable_0.3.3 kernlab_0.9-32
## [77] assertthat_0.2.1 amap_0.8-19
## [79] cachem_1.0.8 xfun_0.39
## [81] rbibutils_2.2.13 S4Arrays_1.2.0
## [83] mime_0.12 RcppEigen_0.3.3.9.3
## [85] tidygraph_1.2.3 survival_3.5-5
## [87] fastICA_1.2-3 statmod_1.5.0
## [89] interactiveDisplayBase_1.40.0 ellipsis_0.3.2
## [91] TH.data_1.1-2 tsne_0.1-3.1
## [93] nlme_3.1-162 naturalsort_0.1.3
## [95] bit64_4.0.5 progress_1.2.2
## [97] gmodels_2.18.1.1 filelock_1.0.2
## [99] bslib_0.5.0 colorspace_2.1-0
## [101] DBI_1.1.3 nnet_7.3-19
## [103] mnormt_2.1.1 tidyselect_1.2.0
## [105] bit_4.0.5 compiler_4.3.1
## [107] curl_5.0.1 rvest_1.0.3
## [109] expm_0.999-7 xml2_1.3.5
## [111] TFisher_0.2.0 ggdendro_0.1.23
## [113] DelayedArray_0.28.0 shadowtext_0.1.2
## [115] bookdown_0.34 rtracklayer_1.62.0
## [117] harmonicmeanp_3.0 sfsmisc_1.1-15
## [119] scales_1.2.1 DEoptimR_1.1-0
## [121] RBGL_1.78.0 rappdirs_0.3.3
## [123] snowfall_1.84-6.2 apcluster_1.4.10
## [125] digest_0.6.33 rmarkdown_2.23
## [127] XVector_0.42.0 htmltools_0.5.5
## [129] pkgconfig_2.0.3 highr_0.10
## [131] dbplyr_2.3.3 fastmap_1.1.1
## [133] rlang_1.1.1 shiny_1.7.4.1
## [135] farver_2.1.1 jquerylib_0.1.4
## [137] zoo_1.8-12 jsonlite_1.8.7
## [139] BiocParallel_1.36.0 mclust_6.0.0
## [141] RCurl_1.98-1.12 magrittr_2.0.3
## [143] modeltools_0.2-23 GenomeInfoDbData_1.2.10
## [145] munsell_0.5.0 Rcpp_1.0.11
## [147] viridis_0.6.3 stringi_1.7.12
## [149] zlibbioc_1.48.0 MASS_7.3-60
## [151] plyr_1.8.8 AnnotationHub_3.10.0
## [153] org.Hs.eg.db_3.17.0 flexmix_2.3-19
## [155] parallel_4.3.1 ggrepel_0.9.3
## [157] Biostrings_2.70.1 graphlayouts_1.0.0
## [159] splines_4.3.1 multtest_2.58.0
## [161] hms_1.1.3 locfit_1.5-9.8
## [163] qqconf_1.3.2 fastcluster_1.2.3
## [165] igraph_1.5.0 reshape2_1.4.4
## [167] biomaRt_2.58.0 BiocVersion_3.18.0
## [169] XML_3.99-0.14 evaluate_0.21
## [171] metap_1.8 pcalg_2.7-8
## [173] BiocManager_1.30.22 tzdb_0.4.0
## [175] tweenr_2.0.2 httpuv_1.6.11
## [177] polyclip_1.10-4 clue_0.3-64
## [179] BiocBaseUtils_1.4.0 ggforce_0.4.1
## [181] xtable_1.8-4 restfulr_0.0.15
## [183] e1071_1.7-13 later_1.3.1
## [185] viridisLite_0.4.2 class_7.3-22
## [187] snow_0.4-4 ggm_2.5
## [189] ellipse_0.4.5 memoise_2.0.1
## [191] AnnotationDbi_1.64.0 GenomicAlignments_1.38.0
## [193] cluster_2.1.4 timechange_0.2.0