Sean C. Bendall

Single-Cell Mass Cytometry of Differential Immune and Drug Responses Across a Human Hematopoietic Continuum

Simultaneous measurement of more than 30 properties in individual human cells is used to characterize signaling in the immune system. Flow cytometry is an essential tool for dissecting the functional complexity of hematopoiesis. We used single-cell “mass cytometry” to examine healthy human bone marrow, measuring 34 parameters simultaneously in single cells (binding of 31 antibodies, viability, DNA content, and relative cell size). The signaling behavior of cell subsets spanning a defined hematopoietic hierarchy was monitored with 18 simultaneous markers of functional signaling states perturbed by a set of ex vivo stimuli and inhibitors. The data set allowed for an algorithmically driven assembly of related cell types defined by surface antigen expression, providing a superimposable map of cell signaling responses in combination with drug inhibition. Visualized in this manner, the analysis revealed previously unappreciated instances of both precise signaling responses that were bounded within conventionally defined cell subsets and more continuous phosphorylation responses that crossed cell population boundaries in unexpected manners yet tracked closely with cellular phenotype. Collectively, such single-cell analyses provide system-wide views of immune signaling in healthy human hematopoiesis, against which drug action and disease can be compared for mechanistic studies and pharmacologic intervention.

viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia

New high-dimensional, single-cell technologies offer unprecedented resolution in the analysis of heterogeneous tissues. However, because these technologies can measure dozens of parameters simultaneously in individual cells, data interpretation can be challenging. Here we present viSNE, a tool that allows one to map high-dimensional cytometry data onto two dimensions, yet conserve the high-dimensional structure of the data. viSNE plots individual cells in a visual similar to a scatter plot, while using all pairwise distances in high dimension to determine each cells location in the plot. We integrated mass cytometry with viSNE to map healthy and cancerous bone marrow samples. Healthy bone marrow automatically maps into a consistent shape, whereas leukemia samples map into malformed shapes that are distinct from healthy bone marrow and from each other. We also use viSNE and mass cytometry to compare leukemia diagnosis and relapse samples, and to identify a rare leukemia population reminiscent of minimal residual disease. viSNE can be applied to any multi-dimensional single-cell technology.

IGF and FGF cooperatively establish the regulatory stem cell niche of pluripotent human cells in vitro.

Distinctive properties of stem cells are not autonomously achieved, and recent evidence points to a level of external control from the microenvironment. Here, we demonstrate that self-renewal and pluripotent properties of human embryonic stem (ES) cells depend on a dynamic interplay between human ES cells and autologously derived human ES cell fibroblast-like cells (hdFs). Human ES cells and hdFs are uniquely defined by insulin-like growth factor (IGF)- and fibroblast growth factor (FGF)-dependence. IGF 1 receptor (IGF1R) expression was exclusive to the human ES cells, whereas FGF receptor 1 (FGFR1) expression was restricted to surrounding hdFs. Blocking the IGF-II/IGF1R pathway reduced survival and clonogenicity of human ES cells, whereas inhibition of the FGF pathway indirectly caused differentiation. IGF-II is expressed by hdFs in response to FGF, and alone was sufficient in maintaining human ES cell cultures. Our study demonstrates a direct role of the IGF-II/IGF1R axis on human ES cell physiology and establishes that hdFs produced by human ES cells themselves define the stem cell niche of pluripotent human stem cells.

Extracting a Cellular Hierarchy from High-dimensional Cytometry Data with SPADE

The ability to analyze multiple single-cell parameters is critical for understanding cellular heterogeneity. Despite recent advances in measurement technology, methods for analyzing high-dimensional single-cell data are often subjective, labor intensive and require prior knowledge of the biological system. To objectively uncover cellular heterogeneity from single-cell measurements, we present a versatile computational approach, spanning-tree progression analysis of density-normalized events (SPADE). We applied SPADE to flow cytometry data of mouse bone marrow and to mass cytometry data of human bone marrow. In both cases, SPADE organized cells in a hierarchy of related phenotypes that partially recapitulated well-described patterns of hematopoiesis. We demonstrate that SPADE is robust to measurement noise and to the choice of cellular markers. SPADE facilitates the analysis of cellular heterogeneity, the identification of cell types and comparison of functional markers in response to perturbations.

Cytometry by Time-of-Flight Shows Combinatorial Cytokine Expression and Virus-Specific Cell Niches within a Continuum of CD8+ T Cell Phenotypes

Cytotoxic CD8(+) T lymphocytes directly kill infected or aberrant cells and secrete proinflammatory cytokines. By using metal-labeled probes and mass spectrometric analysis (cytometry by time-of-flight, or CyTOF) of human CD8(+) T cells, we analyzed the expression of many more proteins than previously possible with fluorescent labels, including surface markers, cytokines, and antigen specificity with modified peptide-MHC tetramers. With 3-dimensional principal component analysis (3D-PCA) to display phenotypic diversity, we observed a relatively uniform pattern of variation in all subjects tested, highlighting the interrelatedness of previously described subsets and the continuous nature of CD8(+) T cell differentiation. These data also showed much greater complexity in the CD8(+) T cell compartment than previously appreciated, including a nearly combinatorial pattern of cytokine expression, with distinct niches occupied by virus-specific cells. This large degree of functional diversity even between cells with the same specificity gives CD8(+) T cells a remarkable degree of flexibility in responding to pathogens.

A deep profiler's guide to cytometry

In recent years, advances in technology have provided us with tools to quantify the expression of multiple genes in individual cells. The ability to measure simultaneously multiple genes in the same cell is necessary to resolve the great diversity of cell subsets, as well as to define their function in the host. Fluorescence-based flow cytometry is the benchmark for this; with it, we can quantify 18 proteins per cell, at >10 000 cells/s. Mass cytometry is a new technology that promises to extend these capabilities significantly. Immunophenotyping by mass spectrometry provides the ability to measure >36 proteins at a rate of 1000 cells/s. We review these cytometric technologies, capable of high-content, high-throughput single-cell assays.

Multiplexed mass cytometry profiling of cellular states perturbed by small-molecule regulators

Mass cytometry facilitates high-dimensional, quantitative analysis of the effects of bioactive molecules on human samples at single-cell resolution, but instruments process only one sample at a time. Here we describe mass-tag cellular barcoding (MCB), which increases mass cytometry throughput by using n metal ion tags to multiplex up to 2n samples. We used seven tags to multiplex an entire 96-well plate, and applied MCB to characterize human peripheral blood mononuclear cell (PBMC) signaling dynamics and cell-to-cell communication, signaling variability between PBMCs from eight human donors, and the effects of 27 inhibitors on this system. For each inhibitor, we measured 14 phosphorylation sites in 14 PBMC types at 96 conditions, resulting in 18,816 quantified phosphorylation levels from each multiplexed sample. This high-dimensional, systems-level inquiry allowed analysis across cell-type and signaling space, reclassified inhibitors and revealed off-target effects. High-content, high-throughput screening with MCB should be useful for drug discovery, preclinical testing and mechanistic investigation of human disease.

Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis

Acute myeloid leukemia (AML) manifests as phenotypically and functionally diverse cells, often within the same patient. Intratumor phenotypic and functional heterogeneity have been linked primarily by physical sorting experiments, which assume that functionally distinct subpopulations can be prospectively isolated by surface phenotypes. This assumption has proven problematic, and we therefore developed a data-driven approach. Using mass cytometry, we profiled surface and intracellular signaling proteins simultaneously in millions of healthy and leukemic cells. We developed PhenoGraph, which algorithmically defines phenotypes in high-dimensional single-cell data. PhenoGraph revealed that the surface phenotypes of leukemic blasts do not necessarily reflect their intracellular state. Using hematopoietic progenitors, we defined a signaling-based measure of cellular phenotype, which led to isolation of a gene expression signature that was predictive of survival in independent cohorts. This study presents new methods for large-scale analysis of single-cell heterogeneity and demonstrates their utility, yielding insights into AML pathophysiology.

Multiplexed ion beam imaging of human breast tumors

Immunohistochemistry (IHC) is a tool for visualizing protein expression that is employed as part of the diagnostic workup for the majority of solid tissue malignancies. Existing IHC methods use antibodies tagged with fluorophores or enzyme reporters that generate colored pigments. Because these reporters exhibit spectral and spatial overlap when used simultaneously, multiplexed IHC is not routinely used in clinical settings. We have developed a method that uses secondary ion mass spectrometry to image antibodies tagged with isotopically pure elemental metal reporters. Multiplexed ion beam imaging (MIBI) is capable of analyzing up to 100 targets simultaneously over a five-log dynamic range. Here, we used MIBI to analyze formalin-fixed, paraffin-embedded human breast tumor tissue sections stained with ten labels simultaneously. The resulting data suggest that MIBI can provide new insights into disease pathogenesis that will be valuable for basic research, drug discovery and clinical diagnostics.

Normalization of mass cytometry data with bead standards.

Mass cytometry uses atomic mass spectrometry combined with isotopically pure reporter elements to currently measure as many as 40 parameters per single cell. As with any quantitative technology, there is a fundamental need for quality assurance and normalization protocols. In the case of mass cytometry, the signal variation over time due to changes in instrument performance combined with intervals between scheduled maintenance must be accounted for and then normalized. Here, samples were mixed with polystyrene beads embedded with metal lanthanides, allowing monitoring of mass cytometry instrument performance over multiple days of data acquisition. The protocol described here includes simultaneous measurements of beads and cells on the mass cytometer, subsequent extraction of the bead‐based signature, and the application of an algorithm enabling correction of both short‐ and long‐term signal fluctuations. The variation in the intensity of the beads that remains after normalization may also be used to determine data quality. Application of the algorithm to a one‐month longitudinal analysis of a human peripheral blood sample reduced the range of median signal fluctuation from 4.9‐fold to 1.3‐fold.