Abstract:High-throughput gene expression data are often obtained from pure or complex (heterogeneous) biological samples. In the latter case, data obtained are a mixture of different cell types and the heterogeneity imposes some difficulties in the analysis of such data. In order to make conclusions on gene expresssion data obtained from heterogeneous samples, methods such as microdissection and flow cytometry have been employed to physically separate the constituting cell types. However, these manual approaches are time consuming when measuring the responses of multiple cell types simultaneously. In addition, exposed samples, on many occasions, end up being contaminated with external perturbations and this may result in an altered yield of molecular content. In this paper, we model the heterogeneous gene expression data using a Bayesian framework, treating the cell type proportions and the cell-type specific expressions as the parameters of the model. Specifically, we present a novel sequential Monte Carlo (SMC) sampler for estimating the model parameters by approximating their posterior distributions with a set of weighted samples. The SMC framework is a robust and efficient approach where we construct a sequence of artificial target (posterior) distributions on spaces of increasing dimensions which admit the distributions of interest as marginals. The proposed algorithm is evaluated on simulated datasets and publicly available real datasets, including Affymetrix oligonucleotide arrays and national center for biotechnology information (NCBI) gene expression omnibus (GEO), with varying number of cell types. The results obtained on all datasets show a superior performance with an improved accuracy in the estimation of cell type proportions and the cell-type specific expressions, and in addition, more accurate identification of differentially expressed genes when compared to other widely known methods for blind decomposition of heterogeneous gene expression data such as Dsection and the nonnegative matrix factorization (NMF) algorithms. MATLAB implementation of the proposed SMC algorithm is available to download at https://github.com/moyanre/smcgenedeconv.git.

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

Computational deconvolution of transcriptomics data from mixed cell populations

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

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

A Novel Computational Complete Deconvolution Method Using RNA-seq Data

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

CDSeq: A Novel Complete Deconvolution Method for Dissecting Heterogeneous Samples Using Gene Expression Data.

A robust workflow to benchmark deconvolution of multi-omic data

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

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

Comprehensive evaluation of deconvolution methods for human brain gene expression

Challenges and perspectives in computational deconvolution of genomics data

DeCovarT, a multidimensional probalistic model for the deconvolution of heterogeneous transcriptomic samples

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

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

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

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

Unraveling the complexity: understanding the deconvolutions of RNA-seq data

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

A sequential Monte Carlo approach to gene expression deconvolution