cBioPortalData: User Start Guide
Marcel Ramos & Levi Waldron
May 20, 2022
Source:vignettes/cBioPortalData.Rmd
cBioPortalData.RmdIntroduction
The cBioPortal for Cancer Genomics website is a great resource for interactive exploration of study datasets. However, it does not easily allow the analyst to obtain and further analyze the data.
We’ve developed the cBioPortalData package to fill this
need to programmatically access the data resources available on the
cBioPortal.
The cBioPortalData package provides an R interface for
accessing the cBioPortal study data within the Bioconductor
ecosystem.
It downloads study data from the cBioPortal API (the full API specification can be found here https://cbioportal.org/api) and uses Bioconductor infrastructure to cache and represent the data.
We use the MultiAssayExperiment
(@Ramos2017-er) package to integrate,
represent, and coordinate multiple experiments for the studies availble
in the cBioPortal. This package in conjunction with
curatedTCGAData give access to a large trove of publicly
available bioinformatic data. Please see our publication here (@Ramos2020-ya).
We demonstrate common use cases of cBioPortalData and
curatedTCGAData during Bioconductor conference workshops.
Overview
Data Structures
Data are provided as a single MultiAssayExperiment per
study. The MultiAssayExperiment representation usually
contains SummarizedExperiment objects for expression data
and RaggedExperiment objects for mutation and CNV-type
data. RaggedExperiment is a data class for representing
‘ragged’ genomic location data, meaning that the measurements per sample
vary.
For more information, please see the RaggedExperiment
and SummarizedExperiment vignettes.
Identifying available studies
As we work through the data, there are some datasest that cannot be
represented as MultiAssayExperiment objects. This can be
due to a number of reasons such as the way the data is handled, presence
of mis-matched identifiers, invalid data types, etc. To see what
datasets are currently not building, we can look refer to
getStudies() with the buildReport = TRUE
argument.
cbio <- cBioPortal()## Warning in .service_validate_md5sum(api_reference_url, api_reference_md5sum, : service version differs from validated version
## service url: https://www.cbioportal.org/api/api-docs
## observed md5sum: 456be3b7a6b5871a1b41f1233ecbae85
## expected md5sum: 1615443badbeaada68463859c34f15f8
studies <- getStudies(cbio, buildReport = TRUE)
head(studies)## # A tibble: 6 × 15
## name description publicStudy groups status importDate allSampleCount
## <chr> <chr> <lgl> <chr> <int> <chr> <int>
## 1 Adrenocortica… "TCGA Adre… TRUE "PUBL… 0 2022-03-0… 92
## 2 Bladder Cance… "Whole exo… TRUE "" 0 2022-03-0… 34
## 3 Basal Cell Ca… "Whole-exo… TRUE "PUBL… 0 2022-03-0… 293
## 4 Acute Lymphob… "Comprehen… TRUE "PUBL… 0 2022-03-0… 93
## 5 Ampullary Car… "Exome seq… TRUE "PUBL… 0 2022-03-0… 160
## 6 Bladder Uroth… "Whole exo… TRUE "PUBL… 0 2022-03-0… 50
## # … with 8 more variables: readPermission <lgl>, studyId <chr>,
## # cancerTypeId <chr>, referenceGenome <chr>, pmid <chr>, citation <chr>,
## # api_build <lgl>, pack_build <lgl>The last two columns will show the availability of each
studyId for either download method (pack_build
for cBioDataPack and api_build for
cBioPortalData).
Choosing download method
There are two main user-facing functions for downloading data from the cBioPortal API.
cBioDataPackmakes use of the tarball distribution of study data. This is useful when the user wants to download and analyze the entirety of the data as available from the cBioPortal.org website.cBioPortalDataallows a more flexibile approach to obtaining study data based on the available parameters such as molecular profile identifiers. This option is useful for users who have a set of gene symbols or identifiers and would like to get a smaller subset of the data that correspond to a particular molecular profile.
Two main functions
cBioDataPack: Obtain Study Data as Zipped Tarballs
This function will access the packaged data from and return an integrative MultiAssayExperiment representation.
## Use ask=FALSE for non-interactive use
laml <- cBioDataPack("laml_tcga", ask = FALSE)
laml## A MultiAssayExperiment object of 12 listed
## experiments with user-defined names and respective classes.
## Containing an ExperimentList class object of length 12:
## [1] cna_hg19.seg: RaggedExperiment with 13571 rows and 191 columns
## [2] cna: SummarizedExperiment with 24776 rows and 191 columns
## [3] linear_cna: SummarizedExperiment with 24776 rows and 191 columns
## [4] methylation_hm27: SummarizedExperiment with 10968 rows and 194 columns
## [5] methylation_hm450: SummarizedExperiment with 10968 rows and 194 columns
## [6] mrna_seq_rpkm_zscores_ref_all_samples: SummarizedExperiment with 19720 rows and 179 columns
## [7] mrna_seq_rpkm_zscores_ref_diploid_samples: SummarizedExperiment with 19719 rows and 179 columns
## [8] mrna_seq_rpkm: SummarizedExperiment with 19720 rows and 179 columns
## [9] mrna_seq_v2_rsem_zscores_ref_all_samples: SummarizedExperiment with 20531 rows and 173 columns
## [10] mrna_seq_v2_rsem_zscores_ref_diploid_samples: SummarizedExperiment with 20440 rows and 173 columns
## [11] mrna_seq_v2_rsem: SummarizedExperiment with 20531 rows and 173 columns
## [12] mutations: RaggedExperiment with 2584 rows and 197 columns
## Functionality:
## experiments() - obtain the ExperimentList instance
## colData() - the primary/phenotype DataFrame
## sampleMap() - the sample coordination DataFrame
## `$`, `[`, `[[` - extract colData columns, subset, or experiment
## *Format() - convert into a long or wide DataFrame
## assays() - convert ExperimentList to a SimpleList of matrices
## exportClass() - save data to flat filescBioPortalData: Obtain data from the cBioPortal API
This function provides a more flexible and granular way to request a MultiAssayExperiment object from a study ID, molecular profile, gene panel, sample list.
acc <- cBioPortalData(api = cbio, by = "hugoGeneSymbol", studyId = "acc_tcga",
genePanelId = "IMPACT341",
molecularProfileIds = c("acc_tcga_rppa", "acc_tcga_linear_CNA")
)## harmonizing input:
## removing 1 colData rownames not in sampleMap 'primary'
acc## A MultiAssayExperiment object of 2 listed
## experiments with user-defined names and respective classes.
## Containing an ExperimentList class object of length 2:
## [1] acc_tcga_rppa: SummarizedExperiment with 57 rows and 46 columns
## [2] acc_tcga_linear_CNA: SummarizedExperiment with 339 rows and 90 columns
## Functionality:
## experiments() - obtain the ExperimentList instance
## colData() - the primary/phenotype DataFrame
## sampleMap() - the sample coordination DataFrame
## `$`, `[`, `[[` - extract colData columns, subset, or experiment
## *Format() - convert into a long or wide DataFrame
## assays() - convert ExperimentList to a SimpleList of matrices
## exportClass() - save data to flat filesNote. To avoid overloading the API service, the API was designed to only query a part of the study data. Therefore, the user is required to enter either a set of genes of interest or a gene panel identifier.
Clearing the cache
cBioDataPack
In cases where a download is interrupted, the user may experience a
corrupt cache. The user can clear the cache for a particular study by
using the removeCache function. Note that this function
only works for data downloaded through the cBioDataPack
function.
removeCache("laml_tcga")cBioPortalData
For users who wish to clear the entire cBioPortalData
cache, it is recommended that they use:
unlink("~/.cache/cBioPortalData/")Example Analysis: Kaplan-Meier Plot
We can use information in the colData to draw a K-M plot
with a few variables from the colData slot of the
MultiAssayExperiment. First, we load the necessary
packages:
We can check the data to lookout for any issues.
table(colData(laml)$OS_STATUS)##
## 0:LIVING 1:DECEASED
## 67 133
class(colData(laml)$OS_MONTHS)## [1] "character"Now, we clean the data a bit to ensure that our variables are of the right type for the subsequent survival model fit.
collaml <- colData(laml)
collaml[collaml$OS_MONTHS == "[Not Available]", "OS_MONTHS"] <- NA
collaml$OS_MONTHS <- as.numeric(collaml$OS_MONTHS)
colData(laml) <- collamlWe specify a simple survival model using SEX as a
covariate and we draw the K-M plot.
fit <- survfit(
Surv(OS_MONTHS, as.numeric(substr(OS_STATUS, 1, 1))) ~ SEX,
data = colData(laml)
)
ggsurvplot(fit, data = colData(laml), risk.table = TRUE)
Data update requests
If you are interested in a particular study dataset that is not currently building, please open an issue at our GitHub repository location and we will do our best to resolve the issues with either the data or the code.
We appreciate your feedback!
sessionInfo
## R version 4.2.0 (2022-04-22)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.4 LTS
##
## Matrix products: default
## BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3
## LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/liblapack.so.3
##
## 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
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] survminer_0.4.9 ggpubr_0.4.0
## [3] ggplot2_3.3.6 survival_3.3-1
## [5] cBioPortalData_2.8.1 MultiAssayExperiment_1.22.0
## [7] SummarizedExperiment_1.26.1 Biobase_2.56.0
## [9] GenomicRanges_1.48.0 GenomeInfoDb_1.32.2
## [11] IRanges_2.30.0 S4Vectors_0.34.0
## [13] BiocGenerics_0.42.0 MatrixGenerics_1.8.0
## [15] matrixStats_0.62.0 AnVIL_1.8.1
## [17] dplyr_1.0.9 BiocStyle_2.24.0
##
## loaded via a namespace (and not attached):
## [1] backports_1.4.1 BiocFileCache_2.4.0
## [3] systemfonts_1.0.4 RCircos_1.2.2
## [5] splines_4.2.0 BiocParallel_1.30.2
## [7] TCGAutils_1.16.0 digest_0.6.29
## [9] htmltools_0.5.2 fansi_1.0.3
## [11] magrittr_2.0.3 memoise_2.0.1
## [13] tzdb_0.3.0 limma_3.52.1
## [15] Biostrings_2.64.0 readr_2.1.2
## [17] vroom_1.5.7 pkgdown_2.0.3
## [19] prettyunits_1.1.1 colorspace_2.0-3
## [21] blob_1.2.3 rvest_1.0.2
## [23] rappdirs_0.3.3 textshaping_0.3.6
## [25] xfun_0.31 crayon_1.5.1
## [27] RCurl_1.98-1.6 jsonlite_1.8.0
## [29] RaggedExperiment_1.20.0 zoo_1.8-10
## [31] glue_1.6.2 GenomicDataCommons_1.20.1
## [33] gtable_0.3.0 zlibbioc_1.42.0
## [35] XVector_0.36.0 DelayedArray_0.22.0
## [37] car_3.0-13 abind_1.4-5
## [39] scales_1.2.0 futile.options_1.0.1
## [41] DBI_1.1.2 rstatix_0.7.0
## [43] Rcpp_1.0.8.3 gridtext_0.1.4
## [45] xtable_1.8-4 progress_1.2.2
## [47] bit_4.0.4 km.ci_0.5-6
## [49] httr_1.4.3 ellipsis_0.3.2
## [51] farver_2.1.0 pkgconfig_2.0.3
## [53] XML_3.99-0.9 rapiclient_0.1.3
## [55] sass_0.4.1 dbplyr_2.1.1
## [57] utf8_1.2.2 RJSONIO_1.3-1.6
## [59] labeling_0.4.2 tidyselect_1.1.2
## [61] rlang_1.0.2 AnnotationDbi_1.58.0
## [63] munsell_0.5.0 tools_4.2.0
## [65] cachem_1.0.6 cli_3.3.0
## [67] generics_0.1.2 RSQLite_2.2.14
## [69] broom_0.8.0 evaluate_0.15
## [71] stringr_1.4.0 fastmap_1.1.0
## [73] yaml_2.3.5 ragg_1.2.2
## [75] knitr_1.39 bit64_4.0.5
## [77] fs_1.5.2 survMisc_0.5.6
## [79] purrr_0.3.4 KEGGREST_1.36.0
## [81] formatR_1.12 xml2_1.3.3
## [83] biomaRt_2.52.0 compiler_4.2.0
## [85] filelock_1.0.2 curl_4.3.2
## [87] png_0.1-7 ggsignif_0.6.3
## [89] tibble_3.1.7 bslib_0.3.1
## [91] stringi_1.7.6 highr_0.9
## [93] futile.logger_1.4.3 GenomicFeatures_1.48.1
## [95] desc_1.4.1 lattice_0.20-45
## [97] Matrix_1.4-1 markdown_1.1
## [99] KMsurv_0.1-5 RTCGAToolbox_2.26.0
## [101] vctrs_0.4.1 pillar_1.7.0
## [103] lifecycle_1.0.1 BiocManager_1.30.18
## [105] jquerylib_0.1.4 data.table_1.14.2
## [107] bitops_1.0-7 rtracklayer_1.56.0
## [109] R6_2.5.1 BiocIO_1.6.0
## [111] bookdown_0.26 gridExtra_2.3
## [113] codetools_0.2-18 lambda.r_1.2.4
## [115] assertthat_0.2.1 rprojroot_2.0.3
## [117] rjson_0.2.21 withr_2.5.0
## [119] GenomicAlignments_1.32.0 Rsamtools_2.12.0
## [121] GenomeInfoDbData_1.2.8 ggtext_0.1.1
## [123] parallel_4.2.0 hms_1.1.1
## [125] grid_4.2.0 tidyr_1.2.0
## [127] rmarkdown_2.14 carData_3.0-5
## [129] restfulr_0.0.13