Graphical Models for Zero-Inflated Single Cell Gene Expression

Abstract

Bulk gene expression experiments relied on aggregations of thousands of cells to measure the average expression in an organism. Advances in microfluidic and droplet sequencing now permit expression profiling in single cells. This study of cell-to-cell variation reveals that individual cells lack detectable expression of transcripts that appear abundant on a population level, giving rise to zero-inflated expression patterns. To infer gene co-regulatory networks from such data, we propose a multivariate Hurdle model. It is comprised of a mixture of singular Gaussian distributions. We employ neighborhood selection with the pseudo-likelihood and a group lasso penalty to select and fit undirected graphical models that capture conditional independences between genes. The proposed method is more sensitive than existing approaches in simulations, even under departures from our Hurdle model. The method is applied to data for T follicular helper cells, and a high-dimensional profile of mouse dendritic cells. It infers network structure not revealed by other methods; or in bulk data sets. An R implementation is available at https://github.com/amcdavid/HurdleNormal .

Fixed error in software URL

References in corpus (6)

• High-dimensional graphs and variable selection with the Lasso
• High-dimensional Ising model selection using ${\ell }_1$-regularized logistic regression
• Variable selection with error control: Another look at Stability Selection
• Graphical Models for Zero-Inflated Single Cell Gene Expression
• Exact block-wise optimization in group lasso and sparse group lasso for linear regression
• Vector-Space Markov Random Fields via Exponential Families

Cited by in corpus (2)

• Graphical Models for Zero-Inflated Single Cell Gene Expression
• Approaches for benchmarking single-cell gene regulatory network inference methods