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

Changlin Wan, Wennan Chang, Yu Zhang, Fenil Shah, Xiaoyu Lu, Yong Zang, Anru Zhang, Sha Cao, Melissa L. Fishel, Qin Ma, Chi Zhang

Research output: Contribution to journalArticle

2 Scopus citations

Abstract

A key challenge in modeling single-cell RNA-seq data is to capture the diversity of gene expression states regulated by different transcriptional regulatory inputs across individual cells, which is further complicated by largely observed zero and low expressions. We developed a left truncated mixture Gaussian (LTMG) model, from the kinetic relationships of the transcriptional regulatory inputs, mRNA metabolism and abundance in single cells. LTMG infers the expression multi-modalities across single cells, meanwhile, the dropouts and low expressions are treated as left truncated. We demonstrated that LTMG has significantly better goodness of fitting on an extensive number of scRNA-seq data, comparing to three other state-of-the-art models. Our biological assumption of the low non-zero expressions, rationality of the multimodality setting, and the capability of LTMG in extracting expression states specific to cell types or functions, are validated on independent experimental data sets. A differential gene expression test and a co-regulation module identification method are further developed. We experimentally validated that our differential expression test has higher sensitivity and specificity, compared with other five popular methods. The co-regulation analysis is capable of retrieving gene co-regulation modules corresponding to perturbed transcriptional regulations. A user-friendly R package with all the analysis power is available at https://github.com/zy26/LTMGSCA.

Original languageEnglish (US)
Pages (from-to)e111
JournalNucleic acids research
Volume47
Issue number18
DOIs
StatePublished - Oct 10 2019

ASJC Scopus subject areas

  • Genetics

Fingerprint Dive into the research topics of 'LTMG: a novel statistical modeling of transcriptional expression states in single-cell RNA-Seq data'. Together they form a unique fingerprint.

  • Cite this