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.

The shape of gene expression distributions matter: how incorporating distribution shape improves the interpretation of cancer transcriptomic data

Deciphering gene expression patterns using large-scale transcriptomic data and its applications

Identification of Gene Expression Pattern Related to Breast Cancer Survival Using Integrated TCGA Datasets and Genomic Tools.

A Novel Approach to Detect Differentially Expressed Genes from Count-Based Digital Databases by Normalizing with Housekeeping Genes.

On the gene expression landscape of cancer

Spatial transcriptomics inferred from pathology whole-slide images links tumor heterogeneity to survival in breast and lung cancer

Distribution shapes govern the discovery of predictive models for gene regulation

A Normalization-Free and Nonparametric Method Sharpens Large-Scale Transcriptome Analysis and Reveals Common Gene Alteration Patterns in Cancers

Common Human Cancer Genes Discovered by Integrated Gene-Expression Analysis

Expression patterns of small numbers of transcripts from functionally-related pathways predict survival in multiple cancers

Estimating the number of available states for normal and tumor tissues in gene expression space

The Prognostic Role of Genes with Skewed Expression Distribution in Lung Adenocarcinoma.

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

Gene Co-expression Is Distance-Dependent in Breast Cancer

Gene signatures for cancer research: A 25-year retrospective and future avenues

Perfect genetic biomarkers for cancer from a fresh view of gene dysregulation

SurvODE: Extrapolating Gene Expression Distribution for Early Cancer Identification

Incorporating prior information in gene expression network-based cancer heterogeneity analysis

Zipf's law in importance of genes for cancer classification using microarray data.

Characterization of Expression-Based Gene Clusters Gives Insights into Variation in Patient Response to Cancer Therapies

Normal tissue transcriptional signatures for tumor-type-agnostic phenotype prediction