Abstract:Single-cell RNA-seq of tumor tissue has greatly strengthened our understanding of tumor development as well as the heterogeneous tissue micro-environment. Many models are dedicated to detect deferentially expressed genes among single cells of different conditions, or to discover (novel) sub-cell types via various clustering methods. However, current statistical models often neglect the zero/low expressions, and use one single modality to fit the rest of the non-zero ones. We believe that zero and low expressions are observed partly due to limited experimental resolution; on the other hand, occasionally, low expressions are falsely observed as non-zero due to the existence of incompletely degraded mRNA. In addition, the severe heterogeneity and fast dynamics of single cells makes it invalid to fit their expressions in one single modality. Here, we propose a left truncated mixture Gaussian (LTMG) distribution to uncover the multimodality in single cells’ gene expression while properly handling the zero and low expressions. LTMG is motivated by the kinetic relationships among the mRNA abundance, transcriptional regulatory signals (TRSs), and mRNA metabolism in a cell. Specifically, observations of a gene’s expression are fitted as mixture Gaussian distributions which assumes zero and low expressions as left censored data. By treating each mixing component as one state of the transcriptional regulators, we combined our LTMG model and our in-house bi-clustering algorithm to infer modules of genes co-regulated by certain (unknown) transcriptional regulatory signal, which is only shared by a subset of all the single cells. Application of the method to three high quality single-cell RNA-seq data sets of T cells extracted from liver, colon and lung cancer micro-environments, identified varied transcriptional regulations specific to subclasses of T cells. Citation Format: Changlin Wan, Wennan Chang, Xiaoyu Lu, Yifan Sun, Kaman So, Sha Cao, Xiongbin Lu, Chi Zhang. Statistical modeling of transcriptional regulatory states in single-cell RNA-Seq data of tumor and infiltrated immune cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5095.

Zero-inflated Smoothing Spline (ZISS) Models for Individual-level Single-cell Temporal Data

A comprehensive assessment of hurdle and zero-inflated models for single cell RNA-sequencing analysis

Data Smoothing Filling Method based on ScRNA-Seq Data Zero-Value Identification

RZiMM-scRNA: A regularized zero-inflated mixture model framework for single-cell RNA-seq data

A general and flexible method for signal extraction from single-cell RNA-seq data

The Poisson distribution model fits UMI-based single-cell RNA-sequencing data

TWO-SIGMA: A Novel Two-Component Single Cell Model-Based Association Method for Single-Cell RNA-seq Data

Statistics or biology: the zero-inflation controversy about scRNA-seq data

ZINBMM: a general mixture model for simultaneous clustering and gene selection using single-cell transcriptomic data

Statistical modeling of isoform splicing dynamics from RNA-seq time series data

Spline regression for zero-inflated models

Identification of cell-type-specific spatially variable genes accounting for excess zeros

Zero-preserving Imputation of Scrna-Seq Data Using Low-Rank Approximation

Variational inference of single cell time series

Zero-inflated Quantile Rank-Score Based Test (ziqrank) with Application to Scrna-Seq Differential Gene Expression Analysis

Zero-preserving imputation of single-cell RNA-seq data

scRecover: Discriminating true and false zeros in single-cell RNA-seq data for imputation

A bivariate zero-inflated negative binomial model and its applications to biomedical settings

Model-based Clustering of Zero-Inflated Single-Cell RNA Sequencing Data via the EM Algorithm

Abstract 5095: Statistical Modeling of Transcriptional Regulatory States in Single-Cell RNA-Seq Data of Tumor and Infiltrated Immune Cells