Abstract:Motivation: The gene expression profile of a cell dictates its function in molecular processes, and can be used to probe its health status. This represents a step forward in the deep characterization of diseases such as cancer and may lead to breakthroughs in their treatment. The technology used to measure the gene expression of isolated cells, single-cell RNA-seq (scRNA-seq), has emerged in the last decade as a key enabler of this progress. However, the use of existing methods for dimensionality reduction, clustering and differential expression is limited by the specificities of the data obtained from scRNA-seq experiments, where technical factors may confound analyses of the true biological signal and contribute to spurious results. To overcome this issue, a possible approach is designing probabilistic generative models of the data with hidden variables encoding different underlying processes. Results: We propose two novel probabilistic models for scRNA-seq data: modified probabilistic count matrix factorization (m-pCMF) and Bayesian zero-inflated negative binomial factorization (ZINBayes). These build upon previous models in the literature while leveraging scalable Bayesian inference via variational methods. We show that the proposed methods are competitive with the state-of-the-art models for robust dimensionality reduction in modern data sets, and improve upon the current best Bayesian model for small numbers of cells. The results show that building probabilistic models of latent variables which encode domain knowledge and using variational inference constitute a promising approach to analyse scRNA-seq data in a scalable way. Availability: m-pCMF and ZINBayes are publicly available as Python packages at https://github.com/pedrofale/, along with the code to reproduce all the results. Contact: susanavinga@tecnico.ulisboa.pt

MGPfact : A Model-Based Factorization Method for scRNA Data Unveils Bifurcating Transcriptional Modules Underlying Cell Fate Determination

MGPfact : A Model-Based Factorization Method for scRNA Data Unveils Bifurcating Transcriptional Modules Underlying Cell Fate Determination

Functional Inference of Gene Regulation Using Single-Cell Multi-Omics

Principled feature attribution for unsupervised gene expression analysis

FISHFactor: A Probabilistic Factor Model for Spatial Transcriptomics Data with Subcellular Resolution

scTrace+: enhance the cell fate inference by integrating the lineage-tracing and multi-faceted transcriptomic similarity information

LTMG: a novel statistical modeling of transcriptional expression states in single-cell RNA-Seq data

Manifold learning analysis suggests strategies to align single-cell multimodal data of neuronal electrophysiology and transcriptomics

Manifold learning-based methods for analyzing single-cell RNA-sequencing data

Extraction of biological signals by factorization enables the reliable analysis of single-cell transcriptomics

Single-cell analysis via manifold fitting: A framework for RNA clustering and beyond

GMFGRN: a matrix factorization and graph neural network approach for gene regulatory network inference

Topological identification and interpretation for single-cell gene regulation elucidation across multiple platforms using scMGCA

Scalable probabilistic matrix factorization for single-cell RNA-seq analysis

BANMF-S: a blockwise accelerated non-negative matrix factorization framework with structural network constraints for single cell imputation

scMerge leverages factor analysis, stable expression, and pseudoreplication to merge multiple single-cell RNA-seq datasets

Dissecting and improving gene regulatory network inference using single-cell transcriptome data

Probabilistic modeling of bifurcations in single-cell gene expression data using a Bayesian mixture of factor analyzers

Comprehensive single-cell RNA-seq analysis using deep interpretable generative modeling guided by biological hierarchy knowledge

A fast and efficient count-based matrix factorization method for detecting cell types from single-cell RNAseq data

A deep generative model for multi-view profiling of single-cell RNA-seq and ATAC-seq data