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.

A mechanistic model for the negative binomial distribution of single-cell mRNA counts

Interpretable and tractable models of transcriptional noise for the rational design of single-molecule quantification experiments

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

A mathematical framework for raw counts of single-cell RNA-seq data analysis

Parameterizing Cell-to-cell Regulatory Heterogeneities Via Stochastic Transcriptional Profiles

Inferring gene regulatory networks from single-cell data: a mechanistic approach

Stochasticity or Noise in Biochemical Reactions

A deep generative model for single-cell RNA sequencing with application to detecting differentially expressed genes

Transcriptome-wide analysis of cell cycle-dependent bursty gene expression from single-cell RNA-seq data using mechanistic model-based inference

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

A deep generative model for gene expression profiles from single-cell RNA sequencing

Quantifying how post-transcriptional noise and gene copy number variation bias transcriptional parameter inference from mRNA distributions

Molecular Distributions in Gene Regulatory Dynamics

Negative binomial count splitting for single-cell RNA sequencing data

Modeling expression ranks for noise-tolerant differential expression analysis of scRNA-seq data

Uncovering the effect of RNA polymerase steric interactions on gene expression noise: analytical distributions of nascent and mature RNA numbers

Distribution shapes govern the discovery of predictive models for gene regulation

Transient power-law behaviour following induction distinguishes between competing models of stochastic gene expression

A Unified Probabilistic Modeling Framework for Eukaryotic Transcription Based on Nascent RNA Sequencing Data

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

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