1 The scpdata package

scpdata disseminates mass spectrometry (MS)-based single-cell proteomics (SCP) data sets formatted using the scp data structure. The data structure is described in the scp vignette.

In this vignette, we describe how to access the SCP data sets. To start, we load the scpdata package.

library("scpdata")

2 Load data from ExperimentHub

The data is stored using the ExperimentHub infrastructure. We first create a connection with ExperimentHub.

eh <- ExperimentHub()

You can list all data sets available in scpdata using the query function.

query(eh, "scpdata")
#> ExperimentHub with 12 records
#> # snapshotDate(): 2021-05-05
#> # $dataprovider: PRIDE, MassIVE, SlavovLab website
#> # $species: Homo sapiens, Mus musculus, Rattus norvegicus, Gallus gallus
#> # $rdataclass: QFeatures
#> # additional mcols(): taxonomyid, genome, description,
#> #   coordinate_1_based, maintainer, rdatadateadded, preparerclass, tags,
#> #   rdatapath, sourceurl, sourcetype 
#> # retrieve records with, e.g., 'object[["EH3899"]]' 
#> 
#>            title            
#>   EH3899 | specht2019v2     
#>   EH3900 | specht2019v3     
#>   EH3901 | dou2019_lysates  
#>   EH3902 | dou2019_mouse    
#>   EH3903 | dou2019_boosting 
#>   ...      ...              
#>   EH3906 | zhu2018NC_lysates
#>   EH3907 | zhu2018NC_islets 
#>   EH3908 | cong2020AC       
#>   EH3909 | zhu2019EL        
#>   EH6011 | liang2020_hela

Another way to get information about the available data sets is to call scpdata(). This will retrieve all the available metadata. For example, we can retrieve the data set titles along with the description to make an informed choice about which data set to choose.

info <- scpdata()
#> snapshotDate(): 2021-05-05
knitr::kable(info[, c("title", "description")])
title description
EH3899 specht2019v2 SCP expression data for monocytes (U-937) and macrophages at PSM, peptide and protein level
EH3900 specht2019v3 SCP expression data for more monocytes (U-937) and macrophages at PSM, peptide and protein level
EH3901 dou2019_lysates SCP expression data for Hela digests (0.2 or 10 ng) at PSM and protein level
EH3902 dou2019_mouse SCP expression data for C10, SVEC or Raw cells at PSM and protein level
EH3903 dou2019_boosting SCP expression data for C10, SVEC or Raw cells and 3 boosters (0, 5 or 50 ng) at PSM and protein level
EH3904 zhu2018MCP Near SCP expression data for micro-dissection rat brain samples (50, 100, or 200 ยตm width) at PSM level
EH3905 zhu2018NC_hela Near SCP expression data for HeLa samples (aproximately 12, 40, or 140 cells) at PSM level
EH3906 zhu2018NC_lysates Near SCP expression data for HeLa lysates (10, 40 and 140 cell equivalent) at PSM level
EH3907 zhu2018NC_islets Near SCP expression data for micro-dissected human pancreas samples (control patients or type 1 diabetes) at PSM level
EH3908 cong2020AC SCP expression data for Hela cells at PSM, peptide and protein level
EH3909 zhu2019EL SCP expression data for chicken utricle samples (1, 3, 5 or 20 cells) at PSM, peptide and protein level
EH6011 liang2020_hela Expression data for HeLa cells (0, 1, 10, 150, 500 cells) at PSM, peptide and protein level

To get one of the data sets (e.g. dou2019_lysates) you can either retrieve it using the ExperimentHub query function

scp <- eh[["EH3901"]]
#> see ?scpdata and browseVignettes('scpdata') for documentation
#> loading from cache
scp
#> An instance of class QFeatures containing 3 assays:
#>  [1] Hela_run_1: SingleCellExperiment with 35224 rows and 10 columns 
#>  [2] Hela_run_2: SingleCellExperiment with 35080 rows and 10 columns 
#>  [3] proteins: SingleCellExperiment with 1641 rows and 20 columns

or you can the use the built-in functions from scpdata

scp <- dou2019_lysates()
#> see ?scpdata and browseVignettes('scpdata') for documentation
#> loading from cache
scp
#> An instance of class QFeatures containing 3 assays:
#>  [1] Hela_run_1: SingleCellExperiment with 35224 rows and 10 columns 
#>  [2] Hela_run_2: SingleCellExperiment with 35080 rows and 10 columns 
#>  [3] proteins: SingleCellExperiment with 1641 rows and 20 columns

3 Data sets information

Each data set has been extensively documented in a separate man page (e.g. ?dou2019_lysates). You can find information about the data content, the acquisition protocol, the data collection procedure as well as the data sources and reference.

4 Data manipulation

For more information about manipulating the data sets, check the scp package. The scp vignette will guide you through a typical SCP data processing workflow. Once your data is loaded from scpdata you can skip section 2 Read in SCP data of the scp vignette.

Session information

R version 4.1.0 (2021-05-18)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 20.04.2 LTS

Matrix products: default
BLAS:   /home/biocbuild/bbs-3.13-bioc/R/lib/libRblas.so
LAPACK: /home/biocbuild/bbs-3.13-bioc/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        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       

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

other attached packages:
 [1] scpdata_1.0.0               ExperimentHub_2.0.0        
 [3] AnnotationHub_3.0.0         BiocFileCache_2.0.0        
 [5] dbplyr_2.1.1                QFeatures_1.2.0            
 [7] MultiAssayExperiment_1.18.0 SummarizedExperiment_1.22.0
 [9] Biobase_2.52.0              GenomicRanges_1.44.0       
[11] GenomeInfoDb_1.28.0         IRanges_2.26.0             
[13] S4Vectors_0.30.0            BiocGenerics_0.38.0        
[15] MatrixGenerics_1.4.0        matrixStats_0.58.0         
[17] BiocStyle_2.20.0           

loaded via a namespace (and not attached):
 [1] ProtGenerics_1.24.0           bitops_1.0-7                 
 [3] bit64_4.0.5                   filelock_1.0.2               
 [5] httr_1.4.2                    tools_4.1.0                  
 [7] bslib_0.2.5.1                 utf8_1.2.1                   
 [9] R6_2.5.0                      DBI_1.1.1                    
[11] lazyeval_0.2.2                withr_2.4.2                  
[13] tidyselect_1.1.1              bit_4.0.4                    
[15] curl_4.3.1                    compiler_4.1.0               
[17] DelayedArray_0.18.0           bookdown_0.22                
[19] sass_0.4.0                    rappdirs_0.3.3               
[21] stringr_1.4.0                 digest_0.6.27                
[23] rmarkdown_2.8                 XVector_0.32.0               
[25] pkgconfig_2.0.3               htmltools_0.5.1.1            
[27] highr_0.9                     fastmap_1.1.0                
[29] rlang_0.4.11                  RSQLite_2.2.7                
[31] shiny_1.6.0                   jquerylib_0.1.4              
[33] generics_0.1.0                jsonlite_1.7.2               
[35] dplyr_1.0.6                   RCurl_1.98-1.3               
[37] magrittr_2.0.1                GenomeInfoDbData_1.2.6       
[39] Matrix_1.3-3                  Rcpp_1.0.6                   
[41] fansi_0.4.2                   MsCoreUtils_1.4.0            
[43] lifecycle_1.0.0               stringi_1.6.2                
[45] yaml_2.2.1                    MASS_7.3-54                  
[47] zlibbioc_1.38.0               grid_4.1.0                   
[49] blob_1.2.1                    promises_1.2.0.1             
[51] crayon_1.4.1                  lattice_0.20-44              
[53] Biostrings_2.60.0             KEGGREST_1.32.0              
[55] knitr_1.33                    pillar_1.6.1                 
[57] glue_1.4.2                    BiocVersion_3.13.1           
[59] evaluate_0.14                 BiocManager_1.30.15          
[61] vctrs_0.3.8                   png_0.1-7                    
[63] httpuv_1.6.1                  purrr_0.3.4                  
[65] clue_0.3-59                   assertthat_0.2.1             
[67] cachem_1.0.5                  xfun_0.23                    
[69] mime_0.10                     xtable_1.8-4                 
[71] AnnotationFilter_1.16.0       later_1.2.0                  
[73] SingleCellExperiment_1.14.0   tibble_3.1.2                 
[75] AnnotationDbi_1.54.0          memoise_2.0.0                
[77] cluster_2.1.2                 ellipsis_0.3.2               
[79] interactiveDisplayBase_1.30.0