This vignette will guide you
through how accessing and manipulating the SCoPE2 data sets available
from the SingleCellMultimodal package.
SCoPE2 is a mass spectrometry (MS)-based single-cell proteomics protocol to quantify the proteome of single-cells in an untargeted fashion. It was initially developed by Specht et al. (2021).
The user can see the available data set by using the default options.
## ah_id mode file_size rdataclass rdatadateadded
## 1 EH4694 protein 33.1 Mb SingleCellExperiment 2020-09-24
## 2 EH4695 rna_assays 68.7 Mb HDF5Matrix 2020-09-24
## 3 EH4696 rna_se 0.2 Mb SingleCellExperiment 2020-09-24
## rdatadateremoved
## 1 <NA>
## 2 <NA>
## 3 <NA>Or by simply running:
## ah_id mode file_size rdataclass rdatadateadded
## 1 EH4694 protein 33.1 Mb SingleCellExperiment 2020-09-24
## 2 EH4695 rna_assays 68.7 Mb HDF5Matrix 2020-09-24
## 3 EH4696 rna_se 0.2 Mb SingleCellExperiment 2020-09-24
## rdatadateremoved
## 1 <NA>
## 2 <NA>
## 3 <NA>Currently, only the macrophage_differentiation is
available.
You can retrieve the actual data from ExperimentHub by
setting dry.run = FALSE. This example retrieves the
complete data set (transcriptome and proteome) for the
macrophage_differentiation project:
## A MultiAssayExperiment object of 2 listed
## experiments with user-defined names and respective classes.
## Containing an ExperimentList class object of length 2:
## [1] protein: SingleCellExperiment with 3042 rows and 1490 columns
## [2] rna: SingleCellExperiment with 32738 rows and 20274 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 filesThis data set has been acquired by the Slavov Lab (Specht et al. (2021)). It contains single-cell proteomics and single-cell RNA sequencing data for macrophages and monocytes. The objective of the research that led to generate the data is to understand whether homogeneous monocytes differentiate in the absence of cytokines to macrophages with homogeneous or heterogeneous profiles. The transcriptomic and proteomic acquisitions are conducted on two separate subset of similar cells (same experimental design). The cell type of the samples are known only for the proteomics data. The proteomics data was retrieved from the authors’ website and the transcriptomic data was retrieved from the GEO database (accession id: GSE142392).
For more information on the protocol, see Specht et al. (2021).
Only version 1.0.0 is currently available.
The macrophage_differentiation data set in this package
contains two assays: rna and protein.
The single-cell proteomics data contains cell type annotation
(celltype), sample preparation batch
(batch_digest and batch_sort), chromatography
batch (batch_chromatography), and the MS acquisition run
(batch_MS). The single-cell transcriptomics data was
acquired in two batches (batch_Chromium). Note that because
the cells that compose the two assays are distinct, there is no common
cell annotation available for both proteomics and transcriptomics. The
annotation were therefore filled with NAs accordingly.
## DataFrame with 21764 rows and 6 columns
## celltype batch_digest batch_sort batch_chromatography
## <character> <character> <character> <character>
## AAACCTGAGAAACCAT-1.1 NA NA NA NA
## AAACCTGAGACTAGGC-1.2 NA NA NA NA
## AAACCTGAGAGGTAGA-1.2 NA NA NA NA
## AAACCTGAGATGCGAC-1.1 NA NA NA NA
## AAACCTGAGGCTAGGT-1.1 NA NA NA NA
## ... ... ... ... ...
## i985 Macrophage Q s8 LCA10
## i986 Monocyte Q s8 LCA10
## i987 Monocyte Q s8 LCA10
## i998 Monocyte R s9 LCB3
## i999 Monocyte R s9 LCB3
## batch_MS batch_Chromium
## <character> <factor>
## AAACCTGAGAAACCAT-1.1 NA 1
## AAACCTGAGACTAGGC-1.2 NA 2
## AAACCTGAGAGGTAGA-1.2 NA 2
## AAACCTGAGATGCGAC-1.1 NA 1
## AAACCTGAGGCTAGGT-1.1 NA 1
## ... ... ...
## i985 X190321S_LCA10_X_FP9.. NA
## i986 X190321S_LCA10_X_FP9.. NA
## i987 X190321S_LCA10_X_FP9.. NA
## i998 X190914S_LCB3_X_16pl.. NA
## i999 X190914S_LCB3_X_16pl.. NAYou can extract and check the transcriptomic data through subsetting:
## class: SingleCellExperiment
## dim: 32738 20274
## metadata(0):
## assays(1): counts
## rownames(32738): MIR1302-10 FAM138A ... AC002321.2 AC002321.1
## rowData names(0):
## colnames(20274): AAACCTGAGAAACCAT-1.1 AAACCTGAGATGCGAC-1.1 ...
## TTTGTCATCGCTTAGA-1.2 TTTGTCATCGTAGATC-1.2
## colData names(1): Batch
## reducedDimNames(0):
## mainExpName: NULL
## altExpNames(0):The data is rather large and is therefore stored on-disk using the HDF5 backend. You can verify this by looking at the assay data matrix. Note that the counts are UMI counts.
## <5 x 5> sparse DelayedMatrix object of type "integer":
## AAACCTGAGAAACCAT-1.1 ... AAACCTGCAATAACGA-1.1
## MIR1302-10 0 . 0
## FAM138A 0 . 0
## OR4F5 0 . 0
## RP11-34P13.7 0 . 0
## RP11-34P13.8 0 . 0The protein assay contains MS-based proteomic data. The
data have been passed sample and feature quality control, normalized,
log transformed, imputed and batch corrected. Detailed information about
the data processing is available in another
vignette. You can extract the proteomic data similarly to the
transcriptomic data:
## class: SingleCellExperiment
## dim: 3042 1490
## metadata(0):
## assays(1): logexprs
## rownames(3042): A0A075B6H9 A0A0B4J1V0 ... Q9Y6X6 Q9Y6Z7
## rowData names(0):
## colnames(1490): i4 i5 ... i2766 i2767
## colData names(5): celltype batch_digest batch_sort batch_chromatography
## batch_MS
## reducedDimNames(0):
## mainExpName: NULL
## altExpNames(0):In this case, the protein data have reasonable size and are loaded
directly into memory. The data matrix is stored in
logexprs. We decided to not use the traditional
logcounts because MS proteomics measures intensities rather
than counts as opposed to scRNA-Seq.
## i4 i5 i7 i10 i11
## A0A075B6H9 -0.01366062 -0.1824640 0.12977307 0.08940234 0.05711272
## A0A0B4J1V0 0.13875137 0.5383824 -0.35823777 -0.10122993 -0.10688821
## A0A0B4J237 0.54897085 -0.2247036 0.50132075 -0.14652437 -0.40384721
## A0A1B0GTH6 0.05392801 -0.3811629 0.07532757 0.29708811 -0.03003732
## A0A1B0GUA6 -0.16910887 0.1542946 -0.14959632 0.14745758 0.03578250## R version 4.4.1 (2024-06-14)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 24.04.1 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/libopenblasp-r0.3.26.so; LAPACK version 3.12.0
##
## 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
##
## time zone: Etc/UTC
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] HDF5Array_1.35.0 rhdf5_2.49.0
## [3] DelayedArray_0.33.1 SparseArray_1.5.45
## [5] S4Arrays_1.5.11 abind_1.4-8
## [7] Matrix_1.7-1 RaggedExperiment_1.31.0
## [9] SingleCellExperiment_1.27.2 SingleCellMultiModal_1.19.0
## [11] MultiAssayExperiment_1.33.0 SummarizedExperiment_1.35.5
## [13] Biobase_2.67.0 GenomicRanges_1.59.0
## [15] GenomeInfoDb_1.43.0 IRanges_2.41.0
## [17] S4Vectors_0.43.2 BiocGenerics_0.53.0
## [19] MatrixGenerics_1.19.0 matrixStats_1.4.1
## [21] BiocStyle_2.35.0
##
## loaded via a namespace (and not attached):
## [1] tidyselect_1.2.1 dplyr_1.1.4 blob_1.2.4
## [4] filelock_1.0.3 Biostrings_2.75.0 fastmap_1.2.0
## [7] BiocFileCache_2.15.0 digest_0.6.37 mime_0.12
## [10] lifecycle_1.0.4 KEGGREST_1.47.0 RSQLite_2.3.7
## [13] magrittr_2.0.3 compiler_4.4.1 rlang_1.1.4
## [16] sass_0.4.9 tools_4.4.1 utf8_1.2.4
## [19] yaml_2.3.10 knitr_1.48 bit_4.5.0
## [22] curl_5.2.3 withr_3.0.2 purrr_1.0.2
## [25] sys_3.4.3 grid_4.4.1 fansi_1.0.6
## [28] ExperimentHub_2.15.0 Rhdf5lib_1.27.0 cli_3.6.3
## [31] rmarkdown_2.28 crayon_1.5.3 generics_0.1.3
## [34] httr_1.4.7 rjson_0.2.23 BiocBaseUtils_1.9.0
## [37] DBI_1.2.3 cachem_1.1.0 zlibbioc_1.51.2
## [40] AnnotationDbi_1.69.0 BiocManager_1.30.25 formatR_1.14
## [43] XVector_0.45.0 vctrs_0.6.5 jsonlite_1.8.9
## [46] bit64_4.5.2 maketools_1.3.1 magick_2.8.5
## [49] jquerylib_0.1.4 glue_1.8.0 BiocVersion_3.21.1
## [52] UCSC.utils_1.1.0 tibble_3.2.1 pillar_1.9.0
## [55] rhdf5filters_1.17.0 rappdirs_0.3.3 htmltools_0.5.8.1
## [58] GenomeInfoDbData_1.2.13 R6_2.5.1 dbplyr_2.5.0
## [61] evaluate_1.0.1 lattice_0.22-6 AnnotationHub_3.15.0
## [64] png_0.1-8 SpatialExperiment_1.15.1 memoise_2.0.1
## [67] bslib_0.8.0 Rcpp_1.0.13 xfun_0.48
## [70] buildtools_1.0.0 pkgconfig_2.0.3