Full spectrum flow cytometry and mass cytometry: A 32-marker panel comparison.

High-dimensional single-cell data has become an important tool in unraveling the complexity of the immune system and its involvement in homeostasis and a large array of pathologies. As technological tools are developed, researchers are adopting them to answer increasingly complex biological questions. Up until recently, mass cytometry (MC) has been the main technology employed in cytometric assays requiring more than 29 markers. Recently, however, with the introduction of full spectrum flow cytometry (FSFC), it has become possible to break the fluorescence barrier and go beyond 29 fluorescent parameters. In this study, in collaboration with the Stanford Human Immune Monitoring Center (HIMC), we compared five patient samples using an established immune panel developed by the HIMC using their MC platform. Using split samples and the same antibody panel, we were able to demonstrate highly comparable results between the two technologies using multiple data analysis approaches. We report here a direct comparison of two technology platforms (MC and FSFC) using a 32-marker flow cytometric immune monitoring panel that can identify all the previously described and anticipated immune subpopulations defined by this panel.

cyCombine allows for robust integration of single-cell cytometry datasets within and across technologies.

Combining single-cell cytometry datasets increases the analytical flexibility and the statistical power of data analyses. However, in many cases the full potential of co-analyses is not reached due to technical variance between data from different experimental batches. Here, we present cyCombine, a method to robustly integrate cytometry data from different batches, experiments, or even different experimental techniques, such as CITE-seq, flow cytometry, and mass cytometry. We demonstrate that cyCombine maintains the biological variance and the structure of the data, while minimizing the technical variance between datasets. cyCombine does not require technical replicates across datasets, and computation time scales linearly with the number of cells, allowing for integration of massive datasets. Robust, accurate, and scalable integration of cytometry data enables integration of multiple datasets for primary data analyses and the validation of results using public datasets.

Vi-Vaccinations Induce Heterogeneous Plasma Cell Responses That Associate With Protection From Typhoid Fever.

Vi-polysaccharide conjugate vaccines are efficacious against cases of typhoid fever; however, an absolute correlate of protection is not established. In this study, we investigated the leukocyte response to a Vi-tetanus toxoid conjugate vaccine (Vi-TT) in comparison with a plain polysaccharide vaccine (Vi-PS) in healthy adults subsequently challenged with Salmonella Typhi. Immunological responses and their association with challenge outcome was assessed by mass cytometry and Vi-ELISpot assay. Immunization induced significant expansion of plasma cells in both vaccines with modest T follicular helper cell responses detectable after Vi-TT only. The Vi-specific IgG and IgM B cell response was considerably greater in magnitude in Vi-TT recipients. Intriguingly, a significant increase in a subset of IgA+ plasma cells expressing mucosal migratory markers α4ß7 and CCR10 was observed in both vaccine groups, suggesting a gut-tropic, mucosal response is induced by Vi-vaccination. The total plasma cell response was significantly associated with protection against typhoid fever in Vi-TT vaccinees but not Vi-PS. IgA+ plasma cells were not significantly associated with protection for either vaccine, although a trend is seen for Vi-PS. Conversely, the IgA- fraction of the plasma cell response was only associated with protection in Vi-TT. In summary, these data indicate that a phenotypically heterogeneous response including both gut-homing and systemic antibody secreting cells may be critical for protection induced by Vi-TT vaccination.

Reversal of epigenetic aging and immunosenescent trends in humans.

Epigenetic "clocks" can now surpass chronological age in accuracy for estimating biological age. Here, we use four such age estimators to show that epigenetic aging can be reversed in humans. Using a protocol intended to regenerate the thymus, we observed protective immunological changes, improved risk indices for many age-related diseases, and a mean epigenetic age approximately 1.5 years less than baseline after 1 year of treatment (-2.5-year change compared to no treatment at the end of the study). The rate of epigenetic aging reversal relative to chronological age accelerated from -1.6 year/year from 0-9 month to -6.5 year/year from 9-12 month. The GrimAge predictor of human morbidity and mortality showed a 2-year decrease in epigenetic vs. chronological age that persisted six months after discontinuing treatment. This is to our knowledge the first report of an increase, based on an epigenetic age estimator, in predicted human lifespan by means of a currently accessible aging intervention.

A clinically meaningful metric of immune age derived from high-dimensional longitudinal monitoring.

Immune responses generally decline with age. However, the dynamics of this process at the individual level have not been characterized, hindering quantification of an individual's immune age. Here, we use multiple 'omics' technologies to capture population- and individual-level changes in the human immune system of 135 healthy adult individuals of different ages sampled longitudinally over a nine-year period. We observed high inter-individual variability in the rates of change of cellular frequencies that was dictated by their baseline values, allowing identification of steady-state levels toward which a cell subset converged and the ordered convergence of multiple cell subsets toward an older adult homeostasis. These data form a high-dimensional trajectory of immune aging (IMM-AGE) that describes a person's immune status better than chronological age. We show that the IMM-AGE score predicted all-cause mortality beyond well-established risk factors in the Framingham Heart Study, establishing its potential use in clinics for identification of patients at risk.

The anatomy of single cell mass cytometry data.

Mass cytometry enables the measurement of up to 50 features on single cell. This has catalyzed a shift toward multidimensional data analysis methods, rather than the manual gating strategies as traditionally for in flow cytometry data. This shift means that data scientists are involved in the analysis process to an increasing degree. As the data is analyzed in a more unbiased fashion, where noisy or uninformative observations are not easily excluded, a deeper knowledge of the origin, noise, and modalities of the data is therefore needed to embark on useful data analysis. In this primer, we introduce the idiosyncrasies of mass cytometry data with a focus on the technical properties of how data generated with the CyTOF® system, and the characteristics of protein expression in the cells of the hematopoietic continuum, specifically targeted toward data scientists. We also provide a comprehensive online repository of scripts, tutorials, and example data. © 2018 International Society for Advancement of Cytometry.

Comparison of CyTOF assays across sites: Results of a six-center pilot study.

For more than five years, high-dimensional mass cytometry has been employed to study immunology. However, these studies have typically been performed in one laboratory on one or few instruments. We present the results of a six-center study using healthy control human peripheral blood mononuclear cells (PBMCs) and commercially available reagents to test the intra-site and inter-site variation of mass cytometers and operators. We used prestained controls generated by the primary center as a reference to compare against samples stained at each individual center. Data were analyzed at the primary center, including investigating the effects of two normalization methods. All six sites performed similarly, with CVs for both Frequency of Parent and median signal intensity (MSI) values<30%. Increased background was seen when using the premixed antibody cocktail aliquots at each site, suggesting that cocktails are best made fresh. Both normalization methods tested performed adequately for normalizing MSI values between centers. Clustering algorithms revealed slight differences between the prestained and the sites-stained samples, due mostly to the increased background of a few antibodies. Therefore, we believe that multicenter mass cytometry assays are feasible.

Platinum-conjugated antibodies for application in mass cytometry.

Mass cytometry has overcome limitations of fluorescent single cell cytometry by allowing for the measurement of up to currently ∼40 different parameters on a single cell level. However, the cellular proteome comprises many more potential analytes, and current mass cytometry instrumentation allows for theoretically up to 121 different mass detection channels. The labeling of specific probes with appropriate metal ions is a significant hurdle for exploiting more of mass cytometry's analytical capacity. To this end, we here describe the labeling of antibody with natural abundance or isotopically purified platinum as formulated in cisplatin and circumventing the use of chelator-loaded polymers. We confirm the utility of cisplatin-antibody-conjugates for surface, intracellular, and phosphoepitope-specific immunophenotyping, as well as for application in cell surface CD45-based barcoding. Cisplatin-labeling of antibody increases the analytical capacity of the CyTOF(®) platform by two channels based on available reagents, and has the potential to add a total of six channels for detection of specific probes, thus helping to better extend the analytical mass range of mass cytometers.

Computationally efficient multidimensional analysis of complex flow cytometry data using second order polynomial histograms.

Many methods have been described for automated clustering analysis of complex flow cytometry data, but so far the goal to efficiently estimate multivariate densities and their modes for a moderate number of dimensions and potentially millions of data points has not been attained. We have devised a novel approach to describing modes using second order polynomial histogram estimators (SOPHE). The method divides the data into multivariate bins and determines the shape of the data in each bin based on second order polynomials, which is an efficient computation. These calculations yield local maxima and allow joining of adjacent bins to identify clusters. The use of second order polynomials also optimally uses wide bins, such that in most cases each parameter (dimension) need only be divided into 4-8 bins, again reducing computational load. We have validated this method using defined mixtures of up to 17 fluorescent beads in 16 dimensions, correctly identifying all populations in data files of 100,000 beads in <10 s, on a standard laptop. The method also correctly clustered granulocytes, lymphocytes, including standard T, B, and NK cell subsets, and monocytes in 9-color stained peripheral blood, within seconds. SOPHE successfully clustered up to 36 subsets of memory CD4 T cells using differentiation and trafficking markers, in 14-color flow analysis, and up to 65 subpopulations of PBMC in 33-dimensional CyTOF data, showing its usefulness in discovery research. SOPHE has the potential to greatly increase efficiency of analysing complex mixtures of cells in higher dimensions.

Multiparameter Phenotyping of Human PBMCs Using Mass Cytometry.

The standard for single-cell analysis of phenotype and function in recent decades has been fluorescence flow cytometry. Mass cytometry is a newer technology that uses heavy metal ions, rather than fluorochromes, as labels for probes such as antibodies. The binding of these ion-labeled probes to cells is quantitated by mass spectrometry. This greatly increases the number of phenotypic and functional markers that can be probed simultaneously. Here, we review topics that must be considered when adapting existing flow cytometry panels to mass cytometry analysis. We present a protocol and representative panels for surface phenotyping and intracellular cytokine staining (ICS) assays.

Barcoding of live human peripheral blood mononuclear cells for multiplexed mass cytometry.

Mass cytometry is developing as a means of multiparametric single-cell analysis. In this study, we present an approach to barcoding separate live human PBMC samples for combined preparation and acquisition on a cytometry by time of flight instrument. Using six different anti-CD45 Ab conjugates labeled with Pd104, Pd106, Pd108, Pd110, In113, and In115, respectively, we barcoded up to 20 samples with unique combinations of exactly three different CD45 Ab tags. Cell events carrying more than or less than three different tags were excluded from analyses during Boolean data deconvolution, allowing for precise sample assignment and the electronic removal of cell aggregates. Data from barcoded samples matched data from corresponding individually stained and acquired samples, at cell event recoveries similar to individual sample analyses. The approach greatly reduced technical noise and minimizes unwanted cell doublet events in mass cytometry data, and it reduces wet work and Ab consumption. It also eliminates sample-to-sample carryover and the requirement of instrument cleaning between samples, thereby effectively reducing overall instrument runtime. Hence, CD45 barcoding facilitates accuracy of mass cytometric immunophenotyping studies, thus supporting biomarker discovery efforts, and it should be applicable to fluorescence flow cytometry as well.

Phenotyping of Live Human PBMC using CyTOF™ Mass Cytometry.

Single-cell analysis has become an method of importance in immunology. Fluorescence flow cytometry has been a major player. However, due to issues such as autofluorescence and emission spillover between different fluorophores, alternative techniques are being developed. In recent years, mass cytometry has emerged, wherein antibodies labeled with metal ions are detected by ICP-MS. In order for a cell to be seen, a metal in the mass window must be present; there is no analogous parameter to forward or side scatter. The current mass window selected is approximately AW 103-196, which includes the lanthanides used for most antibody labeling, as well as iridium and rhodium for DNA intercalators. In this protocol, we use a cocktail of antibodies labeled with MAXPAR metal-chelating polymers to surface-stain live PBMC that have been previously cryopreserved. Many of these markers were taken from a standard fluorescence phenotyping panel (Maecker et al., 2012). No intracellular antibodies are used. We use a CyTOF™ (Cytometry by Time-Of-Flight) mass cytometer to acquire the ICP-MS data. Subsequent analysis of the dual count signal data using FlowJo software allows for cell types to be analyzed based on the dual count signal in each mass channel. The percentage of each cell type is determined and reported as a percent of the parent cell type.

The Split Virus Influenza Vaccine rapidly activates immune cells through Fcγ receptors.

Seasonal influenza vaccination is one of the most common medical procedures and yet the extent to which it activates the immune system beyond inducing antibody production is not well understood. In the United States, the most prevalent formulations of the vaccine consist of degraded or "split" viral particles distributed without any adjuvants. Based on previous reports we sought to determine whether the split influenza vaccine activates innate immune receptors-specifically Toll-like receptors. High-dimensional proteomic profiling of human whole-blood using Cytometry by Time-of-Flight (CyTOF) was used to compare signaling pathway activation and cytokine production between the split influenza vaccine and a prototypical TLR response ex vivo. This analysis revealed that the split vaccine rapidly and potently activates multiple immune cell types but yields a proteomic signature quite distinct from TLR activation. Importantly, vaccine induced activity was dependent upon the presence of human sera indicating that a serum factor was necessary for vaccine-dependent immune activation. We found this serum factor to be human antibodies specific for influenza proteins and therefore immediate immune activation by the split vaccine is immune-complex dependent. These studies demonstrate that influenza virus "splitting" inactivates any potential adjuvants endogenous to influenza, such as RNA, but in previously exposed individuals can elicit a potent immune response by facilitating the rapid formation of immune complexes.

Genetic and environmental determinants of human NK cell diversity revealed by mass cytometry.

Natural killer (NK) cells play critical roles in immune defense and reproduction, yet remain the most poorly understood major lymphocyte population. Because their activation is controlled by a variety of combinatorially expressed activating and inhibitory receptors, NK cell diversity and function are closely linked. To provide an unprecedented understanding of NK cell repertoire diversity, we used mass cytometry to simultaneously analyze 37 parameters, including 28 NK cell receptors, on peripheral blood NK cells from 5 sets of monozygotic twins and 12 unrelated donors of defined human leukocyte antigen (HLA) and killer cell immunoglobulin-like receptor (KIR) genotype. This analysis revealed a remarkable degree of NK cell diversity, with an estimated 6000 to 30,000 phenotypic populations within an individual and >100,000 phenotypes in the donor panel. Genetics largely determined inhibitory receptor expression, whereas activation receptor expression was heavily environmentally influenced. Therefore, NK cells may maintain self-tolerance through strictly regulated expression of inhibitory receptors while using adaptable expression patterns of activating and costimulatory receptors to respond to pathogens and tumors. These findings further suggest the possibility that discrete NK cell subpopulations could be harnessed for immunotherapeutic strategies in the settings of infection, reproduction, and transplantation.

Mass cytometry: protocol for daily tuning and running cell samples on a CyTOF mass cytometer.

In recent years, the rapid analysis of single cells has commonly been performed using flow cytometry and fluorescently-labeled antibodies. However, the issue of spectral overlap of fluorophore emissions has limited the number of simultaneous probes. In contrast, the new CyTOF mass cytometer by DVS Sciences couples a liquid single-cell introduction system to an ICP-MS. Rather than fluorophores, chelating polymers containing highly-enriched metal isotopes are coupled to antibodies or other specific probes. Because of the metal purity and mass resolution of the mass cytometer, there is no "spectral overlap" from neighboring isotopes, and therefore no need for compensation matrices. Additionally, due to the use of lanthanide metals, there is no biological background and therefore no equivalent of autofluorescence. With a mass window spanning atomic mass 103-203, theoretically up to 100 labels could be distinguished simultaneously. Currently, more than 35 channels are available using the chelating reagents available from DVS Sciences, allowing unprecedented dissection of the immunological profile of samples. Disadvantages to mass cytometry include the strict requirement for a separate metal isotope per probe (no equivalent of forward or side scatter), and the fact that it is a destructive technique (no possibility of sorting recovery). The current configuration of the mass cytometer also has a cell transmission rate of only ~25%, thus requiring a higher input number of cells. Optimal daily performance of the mass cytometer requires several steps. The basic goal of the optimization is to maximize the measured signal intensity of the desired metal isotopes (M) while minimizing the formation of oxides (M+16) that will decrease the M signal intensity and interfere with any desired signal at M+16. The first step is to warm up the machine so a hot, stable ICP plasma has been established. Second, the settings for current and make-up gas flow rate must be optimized on a daily basis. During sample collection, the maximum cell event rate is limited by detector efficiency and processing speed to 1000 cells/sec. However, depending on the sample quality, a slower cell event rate (300-500 cells/sec) is usually desirable to allow better resolution between cells events and thus maximize intact singlets over doublets and debris. Finally, adequate cleaning of the machine at the end of the day helps minimize background signal due to free metal.

Development of mass cytometry methods for bacterial discrimination.

Fluorescent flow cytometry has become the method of choice for interrogation of bacterial populations at the single-cell level. However, limitations of this technique include issues of dynamic range, spectral overlap, photobleaching, and overall low signal intensity due to the small size of bacteria. The recent development of mass cytometry allows single-cell analysis with the resolution of inductively coupled plasma mass spectrometry, facilitating multiparametric analysis. Using a combination of a metal-based membrane stain and lectins conjugated to lanthanide-chelating polymers, we demonstrate that individual Escherichia coli cells can be differentiated based on their cell surface polysaccharides using mass cytometry. The model E. coli system involves evaluation of three different surface polysaccharides using element-tagged concanavalin A and wheat germ agglutinin lectins. Finally, this technique enabled experiments designed to follow the export of O-antigen substituted lipopolysaccharide in a conditional mutant. These studies revealed that the culture responds as a uniform population and that lipopolysaccharide export is approximately 10 times faster than the logarithmic bacterial doubling time.