Evgenii Chekalin, Billy Zeng, Patrick Newbury, Benjamin Glicksberg, Jing Xing, Ke Liu, Dimitri Joseph, Bin Chen

Edited on October 25, 2021; Compiled on October 25, 2021

Package overview

As the field of precision medicine progresses, we start to tailor treatments for cancer patients classified not only by their clinical, but also by their molecular features. The emerging cancer subtypes defined by these features require dedicated resources to assist the discovery of drug candidates for preclinical evaluation. Voluminous cancer patient gene expression profiles have been accumulated in public databases, enabling the creation of a cancer-specific expression signature. Meanwhile, large-scale gene expression profiles of chemical compounds became available in recent years. By matching the cancer-specific expression signature to compound-induced gene expression profiles from large drug libraries, researchers can prioritize small molecules that present high potency to reverse expression of signature genes for further experimental testing of their efficacy. This approach has proven to be an efficient and cost-effective way to identify efficacious drug candidates. However, the success of this approach requires multiscale procedures, imposing significant challenges to many labs. Therefore, we present OCTAD (http://octad.org): an open workspace for virtually screening compounds targeting precise cancer patient groups using gene expression features. We have included 19,127 patient tissue samples covering 50 cancer types, and expression profiles for 12,442 distinct compounds. We will continue to include more tissue samples. We streamline all the procedures including deep-learning based reference tissue selection, disease gene expression signature creation, drug reversal potency scoring, and in silico validation. We release OCTAD as a web portal and a standalone R package to allow experimental and computational scientists to easily navigate the tool. The code and data can also be employed to answer various biological questions.

#Data included within the package.

Package includes all required data for drug repurposing OCTAD pipeline. Initialization of the pipeline starts with listing all available samples:

library(octad.db)
## Loading required package: ExperimentHub
## Loading required package: BiocGenerics
## 
## Attaching package: 'BiocGenerics'
## The following objects are masked from 'package:stats':
## 
##     IQR, mad, sd, var, xtabs
## The following objects are masked from 'package:base':
## 
##     Filter, Find, Map, Position, Reduce, anyDuplicated, aperm, append,
##     as.data.frame, basename, cbind, colnames, dirname, do.call,
##     duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
##     lapply, mapply, match, mget, order, paste, pmax, pmax.int, pmin,
##     pmin.int, rank, rbind, rownames, sapply, setdiff, table, tapply,
##     union, unique, unsplit, which.max, which.min
## Loading required package: AnnotationHub
## Loading required package: BiocFileCache
## Loading required package: dbplyr
## Welcome to octad.db package. This is a database package, pipeline available via the octad package. If you want to run the pipeline on the webserver, please refer to octad.org
#select data
suppressMessages(library(octad.db))
phenoDF=get_ExperimentHub_data("EH7274") #load data.frame with samples included in the OCTAD database. 
## see ?octad.db and browseVignettes('octad.db') for documentation
## loading from cache
head(phenoDF) #list all data included within the package
##                  sample.id sample.type                  biopsy.site cancer
## 1 GTEX-1117F-0226-SM-5GZZ7      normal       ADIPOSE - SUBCUTANEOUS normal
## 2 GTEX-1117F-0426-SM-5EGHI      normal            MUSCLE - SKELETAL normal
## 3 GTEX-1117F-0526-SM-5EGHJ      normal              ARTERY - TIBIAL normal
## 4 GTEX-1117F-0626-SM-5N9CS      normal            ARTERY - CORONARY normal
## 5 GTEX-1117F-0726-SM-5GIEN      normal     HEART - ATRIAL APPENDAGE normal
## 6 GTEX-1117F-1326-SM-5EGHH      normal ADIPOSE - VISCERAL (OMENTUM) normal
##   data.source gender read.count.file age_in_year metastatic_site tumor_grade
## 1        GTEX Female      TOIL.RData        <NA>            <NA>        <NA>
## 2        GTEX Female      TOIL.RData        <NA>            <NA>        <NA>
## 3        GTEX Female      TOIL.RData        <NA>            <NA>        <NA>
## 4        GTEX Female      TOIL.RData        <NA>            <NA>        <NA>
## 5        GTEX Female      TOIL.RData        <NA>            <NA>        <NA>
## 6        GTEX Female      TOIL.RData        <NA>            <NA>        <NA>
##   tumor_stage gain_list loss_list mutation_list subtype
## 1        <NA>                              <NA>    <NA>
## 2        <NA>                              <NA>    <NA>
## 3        <NA>                              <NA>    <NA>
## 4        <NA>                              <NA>    <NA>
## 5        <NA>                              <NA>    <NA>
## 6        <NA>                              <NA>    <NA>

Besides, the package includes examples of the output from diffExp and runsRGES functions from the octad package along with desctiption of the data:

data("res_example",package='octad.db') #load example res from octad.db
## Warning in data("res_example", package = "octad.db"): data set 'res_example'
## not found
data("sRGES_example",package='octad.db') #load example sRGES from octad.db
## Warning in data("sRGES_example", package = "octad.db"): data set
## 'sRGES_example' not found

Please note, this is a database-like package for the main package octad which can be obtained here After the package will be accepted to the bioconductor, it will be available on the bioconductor website

Session information

Here is the output of sessionInfo on the system where this document was compiled:

sessionInfo()
## R Under development (unstable) (2023-10-22 r85388)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 22.04.3 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_GB              LC_COLLATE=C              
##  [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] octad.db_1.5.0       ExperimentHub_2.11.0 AnnotationHub_3.11.0
## [4] BiocFileCache_2.11.1 dbplyr_2.4.0         BiocGenerics_0.49.0 
## 
## loaded via a namespace (and not attached):
##  [1] KEGGREST_1.43.0               xfun_0.40                    
##  [3] bslib_0.5.1                   Biobase_2.63.0               
##  [5] bitops_1.0-7                  vctrs_0.6.4                  
##  [7] tools_4.4.0                   generics_0.1.3               
##  [9] stats4_4.4.0                  curl_5.1.0                   
## [11] tibble_3.2.1                  fansi_1.0.5                  
## [13] AnnotationDbi_1.65.0          RSQLite_2.3.2                
## [15] blob_1.2.4                    pkgconfig_2.0.3              
## [17] S4Vectors_0.41.1              GenomeInfoDbData_1.2.11      
## [19] lifecycle_1.0.3               compiler_4.4.0               
## [21] Biostrings_2.71.1             httpuv_1.6.12                
## [23] GenomeInfoDb_1.39.0           htmltools_0.5.6.1            
## [25] sass_0.4.7                    RCurl_1.98-1.12              
## [27] yaml_2.3.7                    interactiveDisplayBase_1.41.0
## [29] pillar_1.9.0                  later_1.3.1                  
## [31] crayon_1.5.2                  jquerylib_0.1.4              
## [33] ellipsis_0.3.2                cachem_1.0.8                 
## [35] mime_0.12                     tidyselect_1.2.0             
## [37] digest_0.6.33                 purrr_1.0.2                  
## [39] dplyr_1.1.3                   BiocVersion_3.19.1           
## [41] fastmap_1.1.1                 cli_3.6.1                    
## [43] magrittr_2.0.3                utf8_1.2.4                   
## [45] withr_2.5.2                   filelock_1.0.2               
## [47] promises_1.2.1                rappdirs_0.3.3               
## [49] bit64_4.0.5                   rmarkdown_2.25               
## [51] XVector_0.43.0                httr_1.4.7                   
## [53] bit_4.0.5                     png_0.1-8                    
## [55] memoise_2.0.1                 shiny_1.7.5.1                
## [57] evaluate_0.22                 knitr_1.45                   
## [59] IRanges_2.37.0                rlang_1.1.1                  
## [61] Rcpp_1.0.11                   xtable_1.8-4                 
## [63] glue_1.6.2                    DBI_1.1.3                    
## [65] BiocManager_1.30.22           jsonlite_1.8.7               
## [67] R6_2.5.1                      zlibbioc_1.49.0