Data exploration, quality control and testing in single-cell qPCR-based gene expression experiments

Abstract

Motivation: Cell populations are never truly homogeneous; individual cells exist in biochemical states that define functional differences between them. New technology based on microfluidic arrays combined with multiplexed quantitative polymerase chain reactions now enables high-throughput single-cell gene expression measurement, allowing assessment of cellular heterogeneity. However, few analytic tools have been developed specifically for the statistical and analytical challenges of single-cell quantitative polymerase chain reactions data.Results: We present a statistical framework for the exploration, quality control and analysis of single-cell gene expression data from microfluidic arrays. We assess accuracy and within-sample heterogeneity of single-cell expression and develop quality control criteria to filter unreliable cell measurements. We propose a statistical model accounting for the fact that genes at the single-cell level can be on (and a continuous expression measure is recorded) or dichotomously off (and the recorded expression is zero). Based on this model, we derive a combined likelihood ratio test for differential expression that incorporates both the discrete and continuous components. Using an experiment that examines treatment-specific changes in expression, we show that this combined test is more powerful than either the continuous or dichotomous component in isolation, or a t-test on the zero-inflated data. Although developed for measurements from a specific platform (Fluidigm), these tools are generalizable to other multi-parametric measures over large numbers of events.Availability: All results presented here were obtained using the SingleCellAssay R package available on GitHub (http://github.com/RGLab/SingleCellAssay).Contact: rgottard@fhcrc.orgSupplementary information: Supplementary data are available at Bioinformatics online.

The immune cell landscape in kidneys of lupus nephritis patients

“…Identification of genes differentially expressed between LN patients and LD controls was done using the framework proposed by McDavid et al 67 , as implemented by Seurat. P-values were corrected for multiple comparisons using the Benjamini-Hochberg method 68 .…”

Defining epithelial cell dynamics and lineage relationships in the developing lacrimal gland

“…tSNE plots were generated to aid in the 2D representation of multidimensional data independent of the clustering algorithm. The 'bimod', likelihood-ratio test for single cell gene expression was used for differential gene analysis (McDavid et al, 2013).…”

TGF-β receptor maintains CD4 T helper cell identity during chronic viral infections

“…The reciprocal analysis in which we gated first on splenic virus-specific CD4 T cells that express EOMES showed that these cells expressed equivalent amounts of PD1 ( Figure 2E), but TGFβ-RII-deficient EOMES + cells were enriched for Ly6C + PSGL1 + cells ( Figure 2F). In addition, IFN-γ-producing Th1 CD4 T cells after GP [67][68][69][70][71][72][73][74][75][76][77] peptide stimulation (from Figure 1F) expressed more granzyme B in the spleen and lung (Figure 2G). Collectively, these results demonstrate that in the absence of TGF-β signaling, virus-specific CD4 T cells showed enhanced terminal differentiation and increased expression of cytotoxic molecules.…”