Abstract:Advancements in RNA-sequencing have dramatically expanded our ability to study gene expression profiles of biological samples in bulk tissue and single cells. Deconvolution of bulk data with single-cell references provides the ability to study relative cell-type proportions, but most methods assume a reference is present for every cell type in bulk data. This is not true in all circumstances--cell types can be missing in single-cell profiles for many reasons. In this study, we examine the impact of missing cell types on deconvolution methods. Our experimental designs are simulation-based, using paired single-cell and single-nucleus data, since single-nucleus RNA-sequencing is able to preserve the nucleus of cell types that would otherwise be missing in a single-cell counterpart. These datasets allow us to examine the missing-cell-type phenomenon in deconvolution with realistic proportions. We apply three deconvolution methods that vary from straightforward to state-of-the-art: non-negative least squares, BayesPrism, and CIBERSORTx. We find that the performance of deconvolution methods is influenced by both the number and the similarity of missing cell types, consistent with prior results. Additionally, we find that missing cell-type profiles can be recovered from residuals using a simple non-negative matrix factorization strategy. We expect our simulation strategies and results to provide a starting point for those developing new deconvolution methods and help improve their to better account for the presence of missing cell types. Building off of our findings on simulated data, we then analyzed data from high-grade serous ovarian cancer; a tumor that has regions of highly variable levels of adipocytes dependent on the region from which it is sampled. We observe results consistent with simulation, namely that expression patterns from cell types likely to be missing appear present in residuals. Our results suggests that deconvolution methods should consider the possibility of missing cell types and provide a starting point to address this. Our source code for data simulation and analysis is freely available at .

Challenges and perspectives in computational deconvolution of genomics data

Challenges and perspectives in computational deconvolution of genomics data

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single cell RNA-sequencing datasets

Computational deconvolution of transcriptomics data from mixed cell populations

An updated State-of-the-Art Overview of transcriptomic Deconvolution Methods

Deconer: A comprehensive and systematic evaluation toolkit for reference-based cell type deconvolution algorithms using gene expression data

Deconvolution from bulk gene expression by leveraging sample-wise and gene-wise similarities and single-cell RNA-Seq data

A robust workflow to benchmark deconvolution of multi-omic data

Deconvolution algorithms for inference of the cell-type composition of the spatial transcriptome

Benchmarking second-generation methods for cell-type deconvolution of transcriptomic data

Missing cell types in single-cell references impact deconvolution of bulk data but are detectable

Fourteen years of cellular deconvolution: methodology, applications, technical evaluation and outstanding challenges

A Novel Computational Complete Deconvolution Method Using RNA-seq Data

AutoGeneS: Automatic Gene Selection Using Multi-Objective Optimization for RNA-seq Deconvolution

Benchmarking and integration of methods for deconvoluting spatial transcriptomic data

A unified model for cell-type resolution genomics from heterogeneous omics data

Advances in mixed cell deconvolution enable quantification of cell types in spatial transcriptomic data

Robust and accurate estimation of cellular fraction from tissue omics data via ensemble deconvolution

Reference-free cell type deconvolution of multi-cellular pixel-resolution spatially resolved transcriptomics data

Advances in mixed cell deconvolution enable quantification of cell types in spatially-resolved gene expression data

CATD: A reproducible pipeline for selecting cell-type deconvolution methods across tissues