Abstract:Abstract Background Single-cell RNA-Sequencing (scRNA-Seq) has provided single-cell level insights into complex biological processes. However, the high frequency of gene expression detection failures in scRNA-Seq data make it challenging to achieve reliable identification of cell-types and Differentially Expressed Genes (DEG). Moreover, with the explosive growth of single-cell data using 10x genomics protocol, existing methods will soon reach the computation limit due to scalability issues. The single-cell transcriptomics field desperately need new tools and framework to facilitate large-scale single-cell analysis. Results In order to improve the accuracy, robustness, and speed of scRNA-Seq data processing, we propose a generalized zero-inflated negative binomial mixture model, “JOINT,” that can perform probability-based cell-type discovery and DEG analysis simultaneously without the need for imputation. JOINT performs soft-clustering for cell-type identification by computing the probability of individual cells, i.e. each cell can belong to multiple cell types with different probabilities. This is drastically different from existing hard-clustering methods where each cell can only belong to one cell type. The soft-clustering component of the algorithm significantly facilitates the accuracy and robustness of single-cell analysis, especially when the scRNA-Seq datasets are noisy and contain a large number of dropout events. Moreover, JOINT is able to determine the optimal number of cell-types automatically rather than specifying it empirically. The proposed model is an unsupervised learning problem which is solved by using the Expectation and Maximization (EM) algorithm. The EM algorithm is implemented using the TensorFlow deep learning framework, dramatically accelerating the speed for data analysis through parallel GPU computing. Conclusions Taken together, the JOINT algorithm is accurate and efficient for large-scale scRNA-Seq data analysis via parallel computing. The Python package that we have developed can be readily applied to aid future advances in parallel computing-based single-cell algorithms and research in various biological and biomedical fields.

Joint CC and Bimax: A Biclustering Method for Single-Cell RNA-Seq Data Analysis

Joint learning dimension reduction and clustering of single-cell RNA-sequencing data

Identifying Cell Subpopulations and Their Genetic Drivers from Single-Cell RNA-Seq Data Using a Biclustering Approach

jSRC: a flexible and accurate joint learning algorithm for clustering of single-cell RNA-sequencing data

UniBic: Sequential Row-Based Biclustering Algorithm for Analysis of Gene Expression Data.

CAbiNet: joint clustering and visualization of cells and genes for single-cell transcriptomics

CDSKNN XMBD : a Novel Clustering Framework for Large-Scale Single-Cell Data Based on a Stable Graph Structure

BREM-SC: A Bayesian Random Effects Mixture Model for Joint Clustering Single Cell Multi-omics Data

JOINT for large-scale single-cell RNA-sequencing analysis via soft-clustering and parallel computing

QUBIC2: A novel biclustering algorithm for large-scale bulk RNA-sequencing and single-cell RNA-sequencing data analysis

Exploring High-throughput Biomolecular Data with Multiobjective Robust Continuous Clustering

A new and effective two-step clustering approach for single cell RNA sequencing data

Microarray Data Biclustering with Multi-objective Immune Optimization Algorithm.

An Effective Biclustering-Based Framework for Identifying Cell Subpopulations From scRNA-seq Data

A probabilistic model-based bi-clustering method for single-cell transcriptomic data analysis

CellBIC: bimodality-based top-down clustering of single-cell RNA sequencing data reveals hierarchical structure of the cell type

nsDCC: dual-level contrastive clustering with nonuniform sampling for scRNA-seq data analysis

SCSMD: Single Cell Consistent Clustering based on Spectral Matrix Decomposition

scMCs: a framework for single-cell multi-omics data integration and multiple clusterings

A New Binary Biclustering Algorithm Based on Weight Adjacency Difference Matrix for Analyzing Gene Expression Data.

A Deep Matrix Factorization Based Approach for Single-Cell RNA-seq Data Clustering.