ADImpute predicts unmeasured gene expression values from single cell RNA-sequencing data (dropout imputation). This R-package combines multiple dropout imputation methods. ADImpute currently supports, by default, and and two novel imputation methods: and . imputes dropouts with the average quantified expression level of the respective gene across the dataset. uses previously learnt regulatory models of gene expression to infer the expression of dropout genes from the expression of other relevant (predictive) genes in the cell. ADImpute consists of 2 fundamental functions: and .
The function allows for dropout imputation with one or more of the supported methods. In order to specify which imputation method(s) to run, pass them via the argument:
In addition to running different methods on the data, ADImpute can also determine which of these performs best for each gene and perform an “Ensemble” imputation, which combines the best performing methods for different genes. First, evaluate methods using to determine the best performing imputation method for each gene. This step sets a fraction of the quantified entries in the input data to zero, applies different imputation methods to the data and compares the imputation results to the original values. This allows ADImpute to determine which method imputes values with the lowest errors for each gene.
RPM <- NormalizeRPM(ADImpute::demo_data)
methods_pergene <- EvaluateMethods(data = RPM,
do = c("Baseline", "DrImpute", "Network"),
cores = 2, net.coef = ADImpute::demo_net)head(methods_pergene)
#> ABCE1 ABL2 ACSL3 ACTL6A ADAM9 ADD1
#> "DrImpute" "DrImpute" "Baseline" "Baseline" "Baseline" "Baseline"After determining which method performs best for each gene, the imputation can be re-done on the original data and the results of different methods combined into an ensemble:
imputed <- Impute(do = "Ensemble", method.choice = methods_pergene,
data = RPM, cores = 2, net.coef = ADImpute::demo_net)str(imputed)
#> List of 4
#> $ Baseline: num [1:696, 1:50] 8.67 9.72 10.51 11.16 9.91 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : Named chr [1:696] "ABCB5" "ABCE1" "ABL2" "ABTB1" ...
#> .. .. ..- attr(*, "names")= chr [1:696] "ABCB5__chr7" "ABCE1__chr4" "ABL2__chr1" "ABTB1__chr3" ...
#> .. ..$ : chr [1:50] "D10631_57" "D71_7" "D3en1_33" "D73_66" ...
#> $ Network : num [1:696, 1:50] 8.67 0 0 0 0 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : Named chr [1:696] "ABCB5" "ABCE1" "ABL2" "ABTB1" ...
#> .. .. ..- attr(*, "names")= chr [1:696] "ABCB5__chr7" "ABCE1__chr4" "ABL2__chr1" "ABTB1__chr3" ...
#> .. ..$ : chr [1:50] "D10631_57" "D71_7" "D3en1_33" "D73_66" ...
#> $ DrImpute: num [1:696, 1:50] 8.672 3.197 5.481 0.564 0.764 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : Named chr [1:696] "ABCB5" "ABCE1" "ABL2" "ABTB1" ...
#> .. .. ..- attr(*, "names")= chr [1:696] "ABCB5__chr7" "ABCE1__chr4" "ABL2__chr1" "ABTB1__chr3" ...
#> .. ..$ : chr [1:50] "D10631_57" "D71_7" "D3en1_33" "D73_66" ...
#> $ Ensemble: num [1:696, 1:50] 8.67 0 0 0 0 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : Named chr [1:696] "ABCB5" "ABCE1" "ABL2" "ABTB1" ...
#> .. .. ..- attr(*, "names")= chr [1:696] "ABCB5__chr7" "ABCE1__chr4" "ABL2__chr1" "ABTB1__chr3" ...
#> .. ..$ : chr [1:50] "D10631_57" "D71_7" "D3en1_33" "D73_66" ...Both the method-specific imputations and the final ensemble results are available for further examination.
Some zeros in the data correspond to genes expressed in the cell, but not captured upon sequencing - the technical dropouts - while others correspond to genes truly not expressed in the cell - the biological zeros. In order to avoid imputation of biological zeros, adapts the well-established approach of for the computation of the probability of each entry to be a technical dropout. A matrix of such probabilities, of the same size as the original data, can be provided by the user, or computed by using ’s approach, as below. To activate this option, provide a value for in the call to , as exemplified below:
str(imputed)
#> List of 3
#> $ imputations :List of 1
#> ..$ Baseline: num [1:696, 1:50] 8.67 9.72 10.51 11.16 9.91 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : Named chr [1:696] "ABCB5" "ABCE1" "ABL2" "ABTB1" ...
#> .. .. .. ..- attr(*, "names")= chr [1:696] "ABCB5__chr7" "ABCE1__chr4" "ABL2__chr1" "ABTB1__chr3" ...
#> .. .. ..$ : chr [1:50] "D10631_57" "D71_7" "D3en1_33" "D73_66" ...
#> $ zerofiltered :List of 1
#> ..$ Baseline: num [1:696, 1:50] 8.67 9.72 10.51 11.16 0 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : Named chr [1:696] "ABCB5" "ABCE1" "ABL2" "ABTB1" ...
#> .. .. .. ..- attr(*, "names")= chr [1:696] "ABCB5__chr7" "ABCE1__chr4" "ABL2__chr1" "ABTB1__chr3" ...
#> .. .. ..$ : chr [1:50] "D10631_57" "D71_7" "D3en1_33" "D73_66" ...
#> $ dropoutprobabilities: num [1:696, 1:50] NA 1 1 1 NA 0 1 0 0 1 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : Named chr [1:696] "ABCB5" "ABCE1" "ABL2" "ABTB1" ...
#> .. .. ..- attr(*, "names")= chr [1:696] "ABCB5__chr7" "ABCE1__chr4" "ABL2__chr1" "ABTB1__chr3" ...
#> .. ..$ : chr [1:50] "D10631_57" "D71_7" "D3en1_33" "D73_66" ...can also take a object as input. In this case, will result in new internal row metadata being added to the object, with the best performing methods per gene. results in new assays being added to the object. If is specified, only the results after setting biological zeros back to zero will be added to the object.
is built in a modular way, which facilitates the addition of custom functions supporting other imputation methods. Two such methods are and , with wrapper functions already contained within . To call these methods, please follow these steps: 1) install scImpute and/or SCRABBLE from their github repositories 2) clone the ADImpute repository 3) copy the lines below to the file Wrap.R in the source R/ folder of ADImpute, line #309. 4) re-load ADImpute using devtools::load_all() on ADImpute’s folder
# # call to scImpute
if('scimpute' %in% tolower(do)){
message('Make sure you have previously installed scImpute via GitHub.\n')
res <- tryCatch(ImputeScImpute(count_path, labeled = is.null(labels),
Kcluster = cell.clusters, labels = labels, drop_thre = drop_thre,
cores = cores, type = type, tr.length = tr.length),
error = function(e){ stop(paste('Error:', e$message,
'\nTry sourcing the Impute_extra.R file.'))})
imputed$scImpute <- log2( (res / scale) + pseudo.count)
}
# call to SCRABBLE
if('scrabble' %in% tolower(do)){
message('Make sure you have previously installed SCRABBLE via GitHub.\n')
res <- tryCatch(ImputeSCRABBLE(data, bulk),
error = function(e) { stop(paste('Error:', e$message,
'\nTry sourcing the Impute_extra.R file.'))})
imputed$SCRABBLE <- log2( (res / scale) + pseudo.count)
rm(res);gc()
}After these steps, and can be run with or with .
sessionInfo()
#> R version 4.4.2 (2024-10-31)
#> 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] SingleCellExperiment_1.29.1 SummarizedExperiment_1.37.0
#> [3] Biobase_2.67.0 GenomicRanges_1.59.1
#> [5] GenomeInfoDb_1.43.2 IRanges_2.41.1
#> [7] S4Vectors_0.45.2 BiocGenerics_0.53.3
#> [9] generics_0.1.3 MatrixGenerics_1.19.0
#> [11] matrixStats_1.4.1 Matrix_1.7-1
#> [13] ADImpute_1.17.0 BiocStyle_2.35.0
#>
#> loaded via a namespace (and not attached):
#> [1] sass_0.4.9 SparseArray_1.7.2 DrImpute_1.0
#> [4] lattice_0.22-6 digest_0.6.37 evaluate_1.0.1
#> [7] grid_4.4.2 fastmap_1.2.0 jsonlite_1.8.9
#> [10] backports_1.5.0 BiocManager_1.30.25 httr_1.4.7
#> [13] kernlab_0.9-33 UCSC.utils_1.3.0 codetools_0.2-20
#> [16] jquerylib_0.1.4 abind_1.4-8 cli_3.6.3
#> [19] rlang_1.1.4 crayon_1.5.3 XVector_0.47.0
#> [22] cachem_1.1.0 DelayedArray_0.33.2 yaml_2.3.10
#> [25] S4Arrays_1.7.1 tools_4.4.2 parallel_4.4.2
#> [28] BiocParallel_1.41.0 checkmate_2.3.2 GenomeInfoDbData_1.2.13
#> [31] buildtools_1.0.0 R6_2.5.1 lifecycle_1.0.4
#> [34] zlibbioc_1.52.0 bslib_0.8.0 Rcpp_1.0.13-1
#> [37] xfun_0.49 sys_3.4.3 knitr_1.49
#> [40] htmltools_0.5.8.1 snow_0.4-4 rmarkdown_2.29
#> [43] maketools_1.3.1 compiler_4.4.2ADImpute was developed in R 4.0.2, under Linux Mint 20, and tested in Linux, OS X and Windows. For further questions, please contact: [email protected]