Erin F. Simonds

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.

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.

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.

Novel mutations in the inhibitory adaptor protein LNK drive JAK-STAT signaling in patients with myeloproliferative neoplasms

Dysregulated Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling due to activation of tyrosine kinases is a common feature of myeloid malignancies. Here we report the first human disease-related mutations in the adaptor protein LNK, a negative regulator of JAK-STAT signaling, in 2 patients with JAK2 V617F-negative myeloproliferative neoplasms (MPNs). One patient exhibited a 5 base-pair deletion and missense mutation leading to a premature stop codon and loss of the pleckstrin homology (PH) and Src homology 2 (SH2) domains. A second patient had a missense mutation (E208Q) in the PH domain. BaF3-MPL cells transduced with these LNK mutants displayed augmented and sustained thrombopoietin-dependent growth and signaling. Primary samples from MPN patients bearing LNK mutations exhibited aberrant JAK-STAT activation, and cytokine-responsive CD34(+) early progenitors were abnormally abundant in both patients. These findings indicate that JAK-STAT activation due to loss of LNK negative feedback regulation is a novel mechanism of MPN pathogenesis.

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.

Single-Cell Profiling Identifies Aberrant STAT5 Activation in Myeloid Malignancies with Specific Clinical and Biologic Correlates

Progress in understanding the molecular pathogenesis of human myeloproliferative disorders (MPDs) has led to guidelines incorporating genetic assays with histopathology during diagnosis. Advances in flow cytometry have made it possible to simultaneously measure cell type and signaling abnormalities arising as a consequence of genetic pathologies. Using flow cytometry, we observed a specific evoked STAT5 signaling signature in a subset of samples from patients suspected of having juvenile myelomonocytic leukemia (JMML), an aggressive MPD with a challenging clinical presentation during active disease. This signature was a specific feature involving JAK-STAT signaling, suggesting a critical role of this pathway in the biological mechanism of this disorder and indicating potential targets for future therapies.

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.

The Initial Phase of an Immune Response Functions to Activate Regulatory T Cells

An early reaction of CD4+ T lymphocytes to Ag is the production of cytokines, notably IL-2. To detect cytokine-dependent responses, naive Ag-specific T cells were stimulated in vivo and the presence of phosphorylated STAT5 molecules was used to identify the cell populations responding to IL-2. Within hours of T cell priming, IL-2-dependent STAT5 phosphorylation occurred primarily in Foxp3+ regulatory T cells. In contrast, the Ag-specific T cells received STAT5 signals only after repeated Ag exposure or memory differentiation. Regulatory T cells receiving IL-2 signals proliferated and developed enhanced suppressive activity. These results indicate that one of the earliest events in a T cell response is the activation of endogenous regulatory cells, potentially to prevent autoimmunity.

A platinum-based covalent viability reagent for single cell mass cytometry

In fluorescence‐based flow cytometry, cellular viability is determined with membrane‐impermeable fluorescent reagents that specifically enter and label plasma membrane‐compromised nonviable cells. A recent technological advance in flow cytometry uses antibodies conjugated to elemental metal isotopes, rather than to fluorophores, to allow signal detection by atomic mass spectrometry. Unhampered by the limitations of overlapping emission fluorescence, mass cytometry increases the number of parameters that can be measured in single cells. However, mass cytometry is unable to take advantage of current fluorescent viability dyes. An alternative methodology was therefore developed here in which the platinum‐containing chemotherapy drug cisplatin was used to resolve live and dead cells by mass cytometry. In a 1‐min incubation step, cisplatin preferentially labeled nonviable cells from both adherent and suspension cultures, resulting in a platinum signal quantifiable by mass cytometry. This protocol was compatible with established sample processing steps for intracellular cytometry. Furthermore, the live/dead ratios were comparable between mass‐ and fluorescence‐based cytometry. Importantly, although cisplatin is a known DNA‐damaging agent, a 1‐min “pulse” of cisplatin did not induce observable DNA damage or apoptotic responses even within 6‐h post‐exposure. Cisplatin can therefore be used as a viability reagent for a wide range of mass cytometry protocols.