Single-cell allele-specific gene expression analysis

Abstract

Allele-specific expression is traditionally studied by bulk RNA sequencing, which measures average gene expression across cells. Single-cell RNA sequencing (scRNA-seq) allows the comparison of expression distribution between the two alleles of a diploid organism, and characterization of allele-specific bursting. Here we describe SCALE, a bioinformatic and statistical framework for allele-specific gene expression analysis by scRNA-seq. SCALE estimates genome-wide bursting kinetics at the allelic level while accounting for technical bias and other complicating factors such as cell size. SCALE detects genes with significantly different bursting kinetics between the two alleles, as well as genes where the two alleles exhibit non-independent bursting processes. Here, we illustrate SCALE on a mouse blastocyst single-cell dataset with step-by-step demonstration from the upstream bioinformatic processing to the downstream biological interpretation of SCALE’s output.