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.

MCMC implementation of the optimal Bayesian classifier for non-Gaussian models: model-based RNA-Seq classification

A Bayesian model selection approach for identifying differentially expressed transcripts from RNA-Seq data

MNBC: a multithreaded Minimizer-based Naïve Bayes Classifier for improved metagenomic sequence classification

Sequential Sampling for Optimal Bayesian Classification of Sequencing Count Data

Scalable optimal Bayesian classification of single-cell trajectories under regulatory model uncertainty

Negative Binomial Additive Model for RNA-Seq Data Analysis

Nonparametric clustering of RNA-sequencing data

Non-parametric Bayesian modelling of digital gene expression data

Imaging and Spectroscopy of Arp 104: A Post-starburst Interacting Pair with Cross-Fuelling?

Bayesian Conditional Gaussian Network Classifiers with Applications to Mass Spectra Classification

A sparse negative binomial mixture model for clustering RNA-seq count data

An MCMC algorithm for base calling in sequencing-by-synthesis

A sequential Monte Carlo approach to gene expression deconvolution

Bayesian Analysis of RNA-Seq Data Using a Family of Negative Binomial Models.

Bayesian Variable Selection in Multinomial Probit Model for Classifying High-Dimensional Data

BELMM: Bayesian model selection and random walk smoothing in time-series clustering

NBLDA: negative binomial linear discriminant analysis for RNA-Seq data

Naïve Bayesian classifiers with multinomial models for rRNA taxonomic assignment

Optimal ROC-Based Classification and Performance Analysis under Bayesian Uncertainty Models

Optimal Bayesian supervised domain adaptation for RNA sequencing data

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit