Airpart : Interpretable statistical models for analyzing allelic imbalance in single-cell datasets
Preprint
Abstract
Motivation: Allelic expression analysis aids in detection of cis -regulatory mechanisms of genetic variation which produce allelic imbalance (AI) in heterozygotes. Measuring AI in bulk data lacking time or spatial resolution has the limitation that cell-type-specific (CTS), spatial-, or time-dependent AI signals may be dampened or not detected. Results: We introduce a statistical method airpart for identifying differential CTS AI from single-cell RNA-sequencing (scRNA-seq) data, or other spatially- or time-resolved datasets. airpart outputs discrete partitions of data, pointing to groups of genes and cells under common mechanisms of cis-genetic regulation. In order to account for low counts in single-cell data, our method uses a Generalized Fused Lasso with Binomial likelihood for partitioning groups of cells by AI signal, and a hierarchical Bayesian model for AI statistical inference. In simulation, airpart accurately detected partitions of cell types by their AI and had lower RMSE of allelic ratio estimates than existing methods. In real data, airpart identified differential AI patterns across cell states and could be used to define trends of AI signal over spatial or time axes. Availability: The airpart package is available as an R/Bioconductor package at https://bioconductor.org/packages/airpart .Keywords
All Related Versions
- Published version: Bioinformatics, 38 (10), 2773.
This publication has 33 references indexed in Scilit:
- The selection and function of cell type-specific enhancersNature Reviews Molecular Cell Biology, 2015
- Orchestrating high-throughput genomic analysis with BioconductorNature Methods, 2015
- De novo assembly of bacterial transcriptomes from RNA-seq dataGenome Biology, 2015
- Single-Cell RNA-Seq Reveals Dynamic, Random Monoallelic Gene Expression in Mammalian CellsScience, 2014
- Single-cell gene expression analysis reveals genetic associations masked in whole-tissue experimentsNature Biotechnology, 2013
- A powerful and flexible statistical framework for testing hypotheses of allele-specific gene expression from RNA-seq dataGenome Research, 2011
- A Statistical Framework for eQTL Mapping Using RNA‐seq DataBiometrics, 2011
- Regulation of embryonic stem cell self-renewal and pluripotency by leukaemia inhibitory factorBiochemical Journal, 2011
- Between a chicken and a grape: estimating the number of human genesGenome Biology, 2010
- Data Analysis Using Stein's Estimator and its GeneralizationsJournal of the American Statistical Association, 1975