Abstract:In immunology studies, flow cytometry is a commonly used multivariate single-cell assay. One key goal in flow cytometry analysis is to detect the immune cells responsive to certain stimuli. Statistically, this problem can be translated into comparing two protein expression probability density functions (pdfs) before and after the stimulus; the goal is to pinpoint the regions where these two pdfs differ. Further screening of these differential regions can be performed to identify enriched sets of responsive cells. In this paper, we model identifying differential density regions as a multiple testing problem. First, we partition the sample space into small bins. In each bin, we form a hypothesis to test the existence of differential pdfs. Second, we develop a novel multiple testing method, called TEAM (Testing on the Aggregation tree Method), to identify those bins that harbor differential pdfs while controlling the false discovery rate (FDR) under the desired level. TEAM embeds the testing procedure into an aggregation tree to test from fine- to coarse-resolution. The procedure achieves the statistical goal of pinpointing density differences to the smallest possible regions. TEAM is computationally efficient, capable of analyzing large flow cytometry data sets in much shorter time compared with competing methods. We applied TEAM and competing methods on a flow cytometry data set to identify T cells responsive to the cytomegalovirus (CMV)-pp65 antigen stimulation. With additional downstream screening, TEAM successfully identified enriched sets containing monofunctional, bifunctional, and polyfunctional T cells. Competing methods either did not finish in a reasonable time frame or provided less interpretable results. Numerical simulations and theoretical justifications demonstrate that TEAM has asymptotically valid, powerful, and robust performance. Overall, TEAM is a computationally efficient and statistically powerful algorithm that can yield meaningful biological insights in flow cytometry studies.

High-Speed Automatic Characterization of Rare Events in Flow Cytometric Data

Ultra High-Throughput Multiparametric Imaging Flow Cytometry: Towards Diffraction-Limited Sub-Cellular Detection

Ultrahigh-Throughput, Real-Time Flow Cytometry for Rare Cell Quantification from Whole Blood

flowAI: automatic and interactive anomaly discerning tools for flow cytometry data

Expanding the use of clustering and dimensionality reduction in high parameter flow cytometry data through machine learning for novel samples.

High-throughput fluorescence lifetime imaging flow cytometry

Machine Learning for Flow Cytometry Data Analysis

From big flow cytometry datasets to smart diagnostic strategies: The EuroFlow approach

Automatic Phenotyping Using Exhaustive Projection Pursuit

Computational analysis of high-throughput flow cytometry data

Flow cytometry data analysis: Recent tools and algorithms

Statistical file matching of flow cytometry data

TEAM: A MULTIPLE TESTING ALGORITHM ON THE AGGREGATION TREE FOR FLOW CYTOMETRY ANALYSIS

An Algorithmic Pipeline for Analyzing Multi-parametric Flow Cytometry Data

Data analysis in flow cytometry: The future just started

FlowAtlas: an interactive tool for high-dimensional immunophenotyping analysis bridging FlowJo with computational tools in Julia

Flow Cytometry: The Next Revolution

Coarsened mixtures of hierarchical skew normal kernels for flow cytometry analyses

High-speed Cell Recognition Algorithm for Ultrafast Flow Cytometer Imaging System

A Scalable Pipeline for High-Throughput Flow Cytometry

Development, application and computational analysis of high-dimensional fluorescent antibody panels for single-cell flow cytometry