key: cord-0030660-05dxwy1k authors: Iyer, Akshay; Hamers, Anouk A. J.; Pillai, Asha B. title: CyTOF(®) for the Masses date: 2022-04-14 journal: Front Immunol DOI: 10.3389/fimmu.2022.815828 sha: cf525941b14d67083f647b70082cd093c643cb3e doc_id: 30660 cord_uid: 05dxwy1k Mass cytometry has revolutionized immunophenotyping, particularly in exploratory settings where simultaneous breadth and depth of characterization of immune populations is needed with limited samples such as in preclinical and clinical tumor immunotherapy. Mass cytometry is also a powerful tool for single-cell immunological assays, especially for complex and simultaneous characterization of diverse intratumoral immune subsets or immunotherapeutic cell populations. Through the elimination of spectral overlap seen in optical flow cytometry by replacement of fluorescent labels with metal isotopes, mass cytometry allows, on average, robust analysis of 60 individual parameters simultaneously. This is, however, associated with significantly increased complexity in the design, execution, and interpretation of mass cytometry experiments. To address the key pitfalls associated with the fragmentation, complexity, and analysis of data in mass cytometry for immunologists who are novices to these techniques, we have developed a comprehensive resource guide. Included in this review are experiment and panel design, antibody conjugations, sample staining, sample acquisition, and data pre-processing and analysis. Where feasible multiple resources for the same process are compared, allowing researchers experienced in flow cytometry but with minimal mass cytometry expertise to develop a data-driven and streamlined project workflow. It is our hope that this manuscript will prove a useful resource for both beginning and advanced users of mass cytometry. Mass cytometry, also termed cytometry by Time-Of-Flight (CyTOF ® ), is a powerful tool for highdimensional and high-throughput single-cell assays. First introduced in 2009 by Bandura et al. (1) , mass cytometry has become an important tool in the analysis of immune cell function/activation due to its high-parameter capabilities. Since 2015, the application of mass cytometry for immunophenotyping in hematopoietic stem cell transplantation (HSCT) (2) (3) (4) (5) (6) , tumor microenvironment (TME) (7) (8) (9) (10) (11) (12) (13) and cancer immunotherapy (6, 9, (14) (15) (16) (17) has significantly expanded. Until recently, fluorescent-based (conventional) flow cytometry was the method of choice for phenotypic and functional analysis of single cells. Standard flow cytometry technologies using 4-or 5-laser data acquisition instruments allow analysis of up to 30 parameters simultaneously. The newer fluorescent-based flow cytometry machines (spectral flow cytometers) measure the total fluorescence in 1 sample and then use an unmixing technology to mathematically separate the specific fluorophore signals (18) . These data acquisition machines can process up to 50 parameters simultaneously; however, practical application typically allows a maximum of 40 parameters (19) . Due to the broader emission spectra of fluorescent probes following laser excitation, overlapping emission spectra remains a significant issue in flow cytometry. Mass cytometry replaces fluorescent labels with non-biologically available metal isotopes with concise mass spectrometry parameters, thereby overcoming the pitfalls associated with overlapping emission spectra and increasing the number of simultaneously analyzable parameters further (20, 21) . In mass cytometry, cells are incubated with a mixture of probes/ antibodies tagged with a unique non-radioactive heavy metal isotope. Single-cell suspensions are nebulized such that each droplet contains a single cell. Individual cells subsequently pass through argon (Ar) plasma, which atomizes and ionizes the sample. This converts each cell into a cloud containing ions of the elements present in or on that cell. A high-pass optic (quadrupole) removes the low-mass (mainly biologic) ions from each cloud (ions with mass below 75 Da), resulting in a cloud containing only those ions corresponding to the isotopeconjugated probes. In the Time of Flight (TOF) chamber, the ions are separated by mass-to-charge ratio. Upon encountering the detector, these ion counts are amplified and converted into electrical signals. Theoretically, 120 parameters can be studied simultaneously. However, realistically the availability of isotopes with sufficient purity as well as antibody conjugation chemistries limit applications to~60 parameters per mass cytometry panel. A single-cell technology generating even more dimensions is single-cell RNA sequencing (scRNAseq), which gives a quantitative measure of gene expression levels per cell. scRNAseq is a powerful genomic tool for dissecting cell populations. However, scRNAseq can only be run on a small number of single cells (limited mainly by increased costs), whereas mass cytometry experiments can acquire data on several times that number (in the range of 10 6 -10 7 cells), facilitating the characterization of rare cell populations. Additionally, mass cytometry can add critical functional i n fo r m a t i on t hr o u gh p ro t e i n a n a l ys e s. W h e n t h e complementary techniques of scRNAseq and mass cytometry are combined, one can rigorously phenotypically and functionally characterize diverse cell populations within a single sample. Mass cytometry can also be used to confirm data derived from scRNAseq. Considering the complexity of the TME, such a multimodal approach yields powerful data applicable to both tumor-intrinsic and tumor-extrinsic effects of immunotherapies in the TME, as well as the correlation of peripheral immune signatures with treatment response or failure or identification of new targets (22) . The purpose of this article is to detail considerations critical to designing and performing a mass cytometry experiment for immunologists and cancer biologists with limited expertise (Figure 1 ). Our target audience includes not only bench scientists and clinicians with knowledge of basic flow cytometry, but also computational scientists and immunotherapy-focused individuals working with mass cytometry datasets. For more detailed reference literature on conventional flow cytometry, we refer the reader to selected reviews, guidelines, and protocols (18, (23) (24) (25) (26) (27) . Paramount to carefully planning mass cytometry experiment design is consideration of the study goal (e.g. identifying multiple new populations in a sample, characterizing an unknown cell population, proportional comparison of multiple wellcharacterized cell populations, novel biomarker discovery, and analysis of protein expression, cell cycle or phosphorylation state, and pharmacokinetics/pharmacodynamics). The large number of parameters simultaneously analyzable by mass cytometry facilitates assessing all of these endpoints leveraging a single antibody panel. In addition, it is important to plan ahead for the sample sources to be used. Each sample source and tissue type has corresponding optimal pre-analysis sample processing and storage considerations; the researcher is referred to existing interactive resources noted elsewhere in this review to raise project-specific questions for clarification. Many experimental factors can affect mass cytometry data, including cell isolation, staining protocol, fixation, and donorspecific biological variation. All of these are described in more detail by Olsen et al. (28) Accounting for these factors, isolation and staining protocols may need several optimization rounds and unique quality controls. An important quality control (QC) approach for reproducibility and staining consistency is an internal control per sample tube, which can be achieved by sample barcoding (see Barcoding) . Another QC element when studying cytokines or transcription factor activation is stimulated versus unstimulated conditions. The sample distribution itself can contribute important QC components, since all cell types are not positive for all markers (i.e. internal positive and negative controls). Proper pairing of antibodies with metal isotopes (see Panel Design) is critical. When opting for standard pre-conjugated antibodies, the main limitation is the extent of the vendor portfolio. Alternatively, lanthanides can be conjugated de novo to purified antibodies using a Maxpar ® X8 antibody labelling kit (Fluidigm, San Francisco, CA) (29) (30) (31) . With lanthanides, a panel can contain up to 37 cellular markers/antibodies ( Figure 2A ). One limitation of this approach is that proteins such as bovine serum albumin (BSA) (often used in the buffers of purified antibodies) can bind lanthanides, resulting in failed antibody conjugation due to adsorption of the lanthanide. In addition to lanthanides, the distinct isotopes of cadmium (31), palladium (32, 33) , indium (32) , platinum (34) , and bismuth (32, 35) can be utilized, to a total of 60 distinct isotopes ( Figure 2A) . The Maxpar ® MCP9 antibody labelling kit (Fluidigm, San Francisco, CA) is specifically designed to conjugate cadmium (Cd) isotopes (31) . Most commercial Qdots used for conventional flow cytometry contain isotopes of cadmium (106) (107) (108) (109) (110) (111) (112) (113) (114) (115) (116) Cd) with either selenium ( 72-82 Se) or tellurium ( 120-130 Te), which are readily available as antibody B A FIGURE 2 | Metal isotopes utilized in mass cytometry. (A) Periodic table summarizing the elements currently available for mass cytometry experiments. Grey, possible contaminating elements; green, elements used to indicate DNA content or cell size; purple, elements available for antibody conjugations; blue, elements available as probes and tracers; orange, elements not yet explored for mass cytometry but of potential future interest. Pink font; isotopes included in QE4 calibration beads. (B) Certain isotopes work well for the use in sample 1 barcoding. Shown is an example of palladium barcoding using 6 isotopes in unique combinations of 3. This strategy generates 20 separate barcodes, allowing 20 individual samples to be combined into 1 single tube. conjugates (36, 37) . More recent options for antibody labels include streptavidin-coated gold or silver nanoparticles (38, 39) and tantalum oxide nanoparticles (40) . Sample throughput can be enhanced, costs reduced, and data quality improved by utilizing sample barcoding (33, 41) . Each experimental sample (e.g. across individuals or treatment groups) can be tagged with a unique isotope of a particular element (e.g., after which all samples are combined into 1 tube. Figure 2B demonstrates a 6-choose-3 barcoding example. Six Palladium (Pd) isotopes can be used to generate 20 different barcodes, where each barcode is created from a combination of any 3 Pd isotopes. Barcoding minimizes the possibility for intersample staining variability, reduces cell-cell doublets, and minimizes the propensity for inter-sample cross-contamination seen in serial runs across individual samples (33, (42) (43) (44) . Sample barcoding minimizes inter-sample staining variability by avoiding sample-to-sample pipetting errors and inconsistent incubation times. In a 6-choose-3 barcoding scheme, a cell-cell doublet will yield an illegal barcode (I.e. a combination of 2 existing barcodes from the scheme) with a positivity for at least 4 out of 6 isotopes, which cannot belong to a single cell event (33) . Doublets between cells within a sample cannot be detected and removed by barcoding alone. Moreover, the use of all possible Pd combinations (any combination of 6 isotopes generating 64 barcodes) can result in miscoding if one or more reagents fail, making it impossible to exclude cell-cell doublets. Sample barcoding can also be utilized to add an internal control into each tube prior to staining. This internal control generates the same results/clustering/cell proportions from one tube to the next and is therefore strongly recommended for optimal data quality and reproducibility. Internal control cells can be cryopreserved and rethawed pooled wild-type mouse splenocytes or PBMC from the same healthy donor or Vericells ® (Biolegend, San Diego, CA) (45, 46) . There are 3 main options in sample barcoding: 1. The Cell-ID ® 20-Plex palladium (Pd) barcoding kit (Fluidigm, San Francisco, CA) uses 6 distinct Pd isotopes to combine up to 20 samples per tube (41, 47, 48) . This is an intracellular method; one limitation is that samples need to be fixed prior to panel staining, so epitopes may be crosslinked in such a way that the corresponding antibody no longer recognizes the intended epitope. Therefore, protocol troubleshooting and optimization is needed prior to applying this approach. 2. Antibody-based live cell barcoding is more flexible as it does not require fixation prior to panel staining. Unique Pd, cadmium (Cd), and/or platinum (Pt) isotopes are conjugated to antibodies directed against ubiquitous epitopes such as CD45 (hematopoietic lineages), b2 microglobulin (class I MHC and CD1 isomers), or CD298 (integral membrane cationic ATPase-associated proteins), and samples are stained with different combinations of these antibodies (42, 44, 49, 50) . Pd, Cd, or Pt are ideal for live cell barcoding primarily because these isotopes are outside the CyTOF ® optimum mass range of 153 to 176 and therefore tend to be less "bright" (51) . Pd and Cd isotopes are well below and Pt is well above the 139-176 mass range of lanthanides and therefore do not influence lanthanide-based antibody detection. Pd live cell barcoding is more labor intensive than Cd, since there is no available kit for Pd conjugations. Recently, Muftuoglu et al. (52) showed that Cd-CD45 barcodes elicit higher signal intensities than Pd-CD45 barcodes, most likely attributable to superior signal resolution because MCP9 polymers used to conjugate Cd chelate a higher number of isotopes as compared to mDOTA (used to conjugate Pd). (This group also showed that it is possible to conjugate Pd isotopes to CD45 antibodies using MCP9, and that this results in an equal signal intensity as for the Cd-CD45 conjugates.) 3. Monoisotopic cisplatin-based live cell barcoding is the simplest and fastest method available. Cisplatin is used to directly label cells, without the need for antibody conjugations (53, 54) . Cisplatin, a chemotherapeutic agent, contains platinum and is available from Fluidigm as any of the following isotopes: 194 Pt, 195 Pt, 196 Pt, and 198 Pt. Thiol-reactive tellurium (TeMal) (55) or osmium and ruthenium tetroxide (56) can be added to any of the 3 barcoding strategies in order to further increase multiplexing capabilities. Mass cytometry sample staining and acquisition induces a high rate of cell loss. Therefore, starting with 800,000 -1 million cells per sample is advisable. Typically, only 50-70% of the sample can be recovered in the data; the remainder is loss due to aggregation on the walls of the spray chamber and injector (28, 57) . Of note, these numbers are based on the CyTOF Helios ® instrument (Fluidigm, San Francisco, CA). With the CyTOF2 ® instrument (Fluidigm, San Francisco, CA), cell recovery is even lower (30-40%) (20) . Of note, this cell loss inside the machine is stochastic in nature and therefore does not appear to introduce sampling bias (51) . Optimal starting cell numbers are highly dependent on the study and the planned sample staining protocol. Studies involving rare cell populations or transcription factors require a larger starting sample size for adequate rigor as compared to studies investigating prevalent subsets. During sample processing and staining, an additional cell loss of 20-30% must be considered. It is important to minimize inter-sample variability in analyzed cell number and viability in order to insure reproducible staining approaches across experiments. When processing samples, tissue digestions, freeze/thaw cycles, and incomplete fixation prior to permeabilization can introduce sampling bias by differentially affecting specific cell populations. Dead cells may compromise flow and mass cytometric data by non-specifically trapping antibodies (58) . In addition, dead cells tend to release DNA, which adheres to cells, causes cell aggregation, and increases cell doublets. Tissue digestions can also cause an overall low cell viability; achieving a high cell viability (> 80%) is important in ensuring high-quality data (59, 60) . To address low cell viability, dead cell removal kits are available [e.g. Miltenyi Biotec Inc (Auburn, CA) and STEMCELL Technologies Inc (Cambridge, MA)] for application prior to sample staining to assist in data QC. A convenient solution for low cell yields prior to sample staining is live cell barcoding. Not only can multiple low-yield samples be combined in a single tube, but a spiked-in internal control can also be added, increasing the total analysis cell number and distributing cell loss in downstream steps across both study samples and controls and thereby preserving a greater fraction of the study sample (7, 45) . For example, Winkels et al. (61) combined barcoded mouse splenocytes with mouse aorta samples, preserving a greater fraction of the murine aorta samples. Where low cell yields and/or poor viability persist, the cellular composition of these tissues can instead be studied in situ using a histologic approach. There are 2 platforms available for this purpose. In imaging mass cytometry [Hyperion ® (Fluidigm)], a laser ablates histological sections stained with metal-labelled antibodies (62, 63) . A more novel metal-based histology platform with increased speed, sensitivity, and image resolution is multiplexed ion beam imaging [MIBI (Ionpath)], which collects data through secondary ions released from the histological slide by primary ion beams (64) (65) (66) . As in conventional flow cytometry, panel design is key to mass cytometry experiment success (36, 57, 67) . The initial marker selection relies heavily on the scientist's combined biological knowledge and familiarity with statistical testing methods, varying depending on the sample type, cell type, and overall experimental objectives. Relevant biological knowledge includes that from literature and from data generated from prior RNA sequencing or conventional flow cytometry. Marker screen kits are available (68, 69) . Beyond isotope-conjugated antibodies, other probes which can be included in mass cytometry panels include tetramers (70) , carbohydrate-binding molecules (71), tellurium-based oxygen sensors (72) , inorganic nanoparticles (73), RNA probes (74, 75) , and modified nucleotides (75, 76) ( Figure 2A ). With the exception of oxygen (72), small molecules or proteins were not detectable by mass cytometry until the Nitz group developed a Tellurium-containing analog of phenylalanine, making it possible to monitor protein synthesis (77) . Poreba et al. have since developed multiple proteaseselective lanthanide-labelled probes for mass cytometry (78) . Since cells are atomized and ionized inside a mass cytometer, the resulting data lacks the Side Scatter (SSC) and Forward Scatter (FSC) parameters used for cell doublet and debris discrimination in conventional flow cytometry. Therefore, mass cytometry relies on the use of a DNA intercalator (see Sample Staining) (7) . As an alternative to the FSC parameter, Osmium Tetroxide (OsO 4 ) has been suggested as a useful tool to reconstruct cell size in mass cytometry data (79) . OsO4 is a nonpolar compound that penetrates charged membranes and can be detected directly by the mass cytometer. Good et al. have also adapted carboxy-fluorescein succinimidyl ester (CFSE)based protocols for tracking cell proliferation in mass cytometry using a metal-conjugated CFSE cross-reactive antifluorescein isothiocyanate (anti-FITC) antibody (80) . The next step in panel design is pairing an antibody or probe with a metal isotope in a manner that insures optimal signal intensity with minimal to no signal overlap. Critical considerations include: 1) isotope sensitivity range of the detection instrument, 2) intensity of surface marker expression, 3) degree of variation and patterns of expression across samples, and 4) spillover/background. 1. CyTOF ® is most sensitive in the range from atomic mass 153 to 176; therefore isotopes/mass tags within this mass range are preferable for antibodies against weakly expressed markers (51). 2. Surface marker intensity considerations follow a process akin to fluorophore-based panel design for flow cytometry; antigens/probes are first classified as either high expression (primary), medium/variable expression (secondary), or low/ unknown expression (tertiary) (81) . Antibodies with low binding affinity or directed against tertiary antigens should be paired with isotopes in the detection instrument's highsensitivity detection range. Antibodies against primary or secondary antigens need pairing with isotopes on either end of this optimal mass range. 3. The same antigen may have vastly different expression patterns depending on cell type, organ, or disease state (36) . For this reason, it is important to either have or obtain knowledge of the specific antigens in the study. Prior knowledge from conventional flow cytometry and literature will help assign antigens to the above-mentioned categories. For example, CD4 is a primary antigen and exhibits a clear bimodal expression, with clear negative and positive populations (82) . Alternatively, chemokine receptors such as CCR7 are often classified as secondary antigens and have a broad, often non-modal spectrum of expression (83). Relative to the cellular autofluorescence or channel crosstalk seen in conventional flow cytometry, sources of background are greatly reduced in mass cytometry (20, 36) . Mass cytometry background is predominantly caused by signal spillover related to instrument detection sensitivity. In a TOF analyzer ions are separated based on velocity, which in turn is determined by their mass (M) and kinetic energy. Ions of the same kind have small differences in initial position and velocity from each other resulting in slightly different detector arrival times, which is reflected in the width of the resulting mass peak. An over-abundance of the same ions causes position and velocity spreads, resulting in broader mass peaks spilling over into the adjacent mass peak (M+/-1). If an antibody against a high-expressing antigen is conjugated to a metal isotope within the high-sensitivity range of the instrument, spillover will occur due to abundance sensitivity. The second cause of spillover is oxidization of certain metal isotopes following air exposure, resulting in a background signal at 16 mass units ( 16 O) higher than the mass of the primary isotope (M+16). Oxide formation can occur in lanthanum (La), praseodymium (Pr), neodymium (Nd), and samanium (Sm)-labelled antibodies/probes. Only 7 metal isotopes form significant oxides: 139 La, [142] [143] [144] Nd, 148 Nd, and 150 Nd. Spillover matrices are available to assist in this process (84) . Before sample acquisition, CyTOF ® machine parameters are optimized to limit the 139 La oxidation to less than 3% of 155 Gd (= M+16). Finally, the largest contribution to signal overspill is isotope impurity, i.e. a contamination of a metal with one of its other isotopes. 100% purity is not pragmatically feasible for all metal isotopes. The overall principles of panel design (36) are as follows: 1. Tertiary antigens should be paired with isotopes within the CyTOF ® high-sensitivity range and primary antigens should be paired with isotopes outside this range; 2. Do not place a tertiary antigen-isotope pair in the oxide-mass (M+16) of a primary antigen-isotope pair; 3. Choose isotope tags for tertiary antigens in channels which receive no or little spillover from adjacent channels; 4. For less pure isotopes, select antigens that identify specific cell subsets (for example CD4 + and CD8 + subsets, which are mutually exclusive outside the gut and thymus); 5. Channels with high spillover can be reserved for markers to be excluded from downstream analysis (for example CD41 to gate out platelets and platelet-cell aggregates). There are many reference resources that support mass cytometry panel design, including Fluidigm's online tools, institute mass cytometry core facilities, expert collaborators, and key publications (67) . Additionally, available MaxPar panel kits (Fluidigm, San Francisco, CA) (109, 110) include a user-friendly kit-specific data analysis platform (GemStone Software, Topsham, ME). Although the list of pre-conjugated antibodies for purchase is steadily growing, customized antibody panels often require inhouse conjugations. Following panel design, unlabeled antibodies must be conjugated to selected isotopes. Fortunately, conjugation kits and published protocols are available (29, 31, (33) (34) (35) . The isotope planned for conjugation is first linked to a polymer via a chelator. Common chelators used for mass cytometry isotope conjugations include diethylene triamine pentaacetic acid (DTPA), ethylenediamine tertraacetic acid (EDTA), and 1,4,7,10-tetraazacyclododecane-1,4,7,10tetraacetic acid (DOTA). The antibody is separately modified in its hinge region by a reduction of disulfide bonds to thiols using tris-2-carboxyethyl phosphine (TCEP). Finally, the polymer and associated chelate are coupled to a thiol group of the reduced antibody. These methods can be applied to conjugate metal isotopes to IgG antibodies. Buffers containing protein or glycerol as antibody stabilizers should be avoided. For purified antibodies only available with BSA, BSA removal kits are available that can be used before proceeding to the antibody conjugations (111) . Methods for confirming successful conjugation vary depending on the protocol. This is elegantly outlined by Han et al. (29) The shelf life of antibodies conjugated "in-house" can vary significantly. If stored properly they generally maintain functionality at least 6 months from conjugation (29, 112) . Antibodies unused/stored for longer periods will require testing on the mass cytometer to confirm that the isotope remains conjugated. How these tests are best performed is explained in the section "Sample Acquisition & Data Output". Once confirmed, the antibody can be used for new cells of interest following standard titration. Results of subsequent titrations should be compared to the initial titration results to identify and troubleshoot new issues. The next crucial step is to titrate all antibodies and perform a test run involving the entire experimental protocol. This achieves the best possible signal-to-noise ratio by reducing non-specific antibody binding and spillover (113, 114) . In addition, using antibodies at non-saturated concentrations prevents ion detector saturation (33) . Similar to conventional flow cytometry, a serial dilution strategy of at least 5 dilutions is advisable. Antibodies to primary antigens (e.g. common lineage markers such as CD3, CD19, and CD11b) should be titrated individually and separately. Subsequently, antibodies directed against secondary and tertiary antigens can be titrated within the combined antibody panel, in the presence of the antibodies already optimally titrated against primary antigens. This approach has 3 advantages: 1) to titrate antibodies for staining the population of interest, 2) to enrich the signal by gating cells known to express particular secondary and/or tertiary antigens, and 3) to provide internal positive and negative controls within the titration samples. If titrating for signaling molecules such as cytokines or transcription factors, both a baseline sample and an activated sample (stimulation or treatment) are needed (115) . It is advisable to select a cell number per titration point that is comparable to the actual experiment. The antibody titers are determined by calculating the staining index, a method very similar to that used in conventional flow cytometry (116, 117) . The main difference between mass cytometry titrations and those used in flow cytometry is that the standard deviation of the negative population is essentially non-existent in mass cytometry titrations and is therefore not included in the staining index formula (118) . Finally, the performance of the panel should be tested on a few control samples prior to proceeding to valuable experimental samples. If the specificity of a signal is unclear, "metal-minusone" (MMO) controls can be used for resolution (117) . If spillover persists with MMO, Chevrier et al. developed a set of computational tools to compensate for spillover in mass cytometry data (see Compensation) (84) . In some situations, tertiary antigens are so weakly expressed that they are not cleanly discernible. Switching to a 2-step staining can augment the signal from the weakly staining primary antibody; a primary antibody conjugated to biotin, FITC, phycoerythrin (PE), or allophycocyanin (APC) is followed by a metal isotope-labelled secondary antibody or streptavidin. Intelligent mass cytometry panel design is an iterative process often requiring multiple revisions for optimization. Before sample staining, ensure that all buffers are clean by running these solutions at a 1:10,000 dilution (in ddH 2 O or Fluidigm's cell acquisition solution [CAS] ) in the machine's "solution" mode. This process is further elaborated in "Sample Acquisition & Data Output". Contaminants including barium (Ba) from dish soap used in labware cleaning; or lead (Pb), mercury (Hg) and tin (Sn) from water pipes/distilled water are a common challenge in mass cytometry (Figure 2A ) (36) . An abundance of Ba contacting the detector also damages the detector over time and decreases detector lifespan (51) . When working with patient samples, certain therapeutic reagents (cisplatin in cancer chemotherapy, gold in autoimmune therapies), or contrast reagents (Iodine, Ba) can circulate in the patient and contaminate the tissue under study, thus confounding the data ( Figure 2A ) (54, 119) . Alternative intercalators are available, including rhodium-103 ( 103 Rh) (119) . Samples often require pre-processing prior to staining. For PBMC, the anticoagulant used in the blood collection can affect specific cell types and thereby adversely impact the data (120, 121) . It is important to optimize tissue digestion protocols to minimize cell debris (reduces staining quality) and maximize viability. Multiple published protocols exist for a variety of human and mouse tissues, including tumors (43, 57, 61, (122) (123) (124) (125) (126) (127) (128) (129) (130) . In general, when attempting to enrich cells, protocols should leave the cells of interest unmanipulated (e.g., negative selection procedures for magnetic bead separation or fluorescence-activated cell sorting/FACS). Heavy metals in magnetic beads can interfere with the mass cytometer, so careful washing post-enrichment is required (131) . Immune cells from digested tissues can be enriched without antibodies or magnetic beads via density gradient isolations using agents such as Percoll ® (Cytivia, Marlborough, MA), Ficoll ® (Cytivia, Marlborough, MA), and Lympholyte ® (Cedarlane, Burlington, LC). Density gradients, can cause differential loss reduced numbers of certain cell types (for example, granulocytes after Ficoll ® gradient isolations) (132) . When working with cryopreserved samples, the effect of freezing and thawing on target epitopes needs to be tested by conventional flow cytometry prior to initiating a mass cytometry experiment. Freezing can greatly alter surface expression of certain surface antigens (CD62L or PD-1) and cytokines, due to down-regulation under cellular stress (133) (134) (135) . This can be mitigated by maintaining the samples overnight in cell culture media to allow them to equilibrate following thaw (136) . For immune cells, a common media is RPMI with 10% heatinactivated fetal bovine serum (DFBS). To measure cytokine expression potential, cells can be stimulated with phorbol myristate acetate (PMA) (an activator of NF-kb) and ionomycin. To block secretion and thus loss of intracellular signal, at the end of stimulation (2 hours minimum, maximum overnight), treatment with brefeldin A (an inhibitor of protein transport between the endoplasmic reticulum ad the Golgi apparatus) or monensin (an inhibitor of trans-Golgi transport) is required (137) (138) (139) . Cryopreservation can also differentially affect the relative frequencies of viable subpopulations upon rethaw, and this needs to be optimized to minimize selection bias (140) . When a studying phospho-proteins, because phosphorylation and dephosphorylation is a very rapid and dynamic process, capture of the phospho-protein may require exploratory conventional flow cytometry studies to codify optimal conditions and timing. Samples are often fixed prior to staining, a procedure which must be optimized to minimize selection bias (141) (142) (143) . It is critical that all antibodies in the panel be confirmed to bind target after fixation, due to epitope denaturation by fixation (144) . Occasionally, a new antibody clone must be titrated to replace the clone no longer functional post-fixation. Another approach is to forego sample fixation entirely, and block dephosphorylation by pervanadate incubation for 5 minutes prior to staining (145, 146) . Because mass cytometry target detection relies on antibody Fab' -target interaction, constant fragment (Fc) receptors need to be blocked prior to antibody staining in order to reduce false positive signal from antibody binding via their Fc region (as in conventional flow cytometry) (147) . If, however, CD16 (148) (149) (150) (151) (152) (153) and CD32 (154, 155) are critical targets, other options include 1) non-specific protein saturation with extra serum or BSA or 2) staining with anti-CD16 and anti-CD32 antibodies in a separate first step, prior to the staining with the remainder of the antibody panel. (The latter approach also leverages competitive inhibition by using the detecting antibodies for simultaneous Fab'-mediated CD16 and CD32 detection and blockade). Unlike conventional flow cytometry, in which antibody staining is performed at 4°C, mass cytometry staining can be performed at room temperature as internalization of antigens does not alter detection. There are however specific situations in which staining at 4°C is advisable, such as for myeloid populations that adhere to plastic if metabolically active above 4°C or to better preserve cell activation and viability prior to specific functional assays. After Fc blockade and live sample barcoding, every staining begins with a live/dead cell discernment step using cisplatin. Cisplatin quickly and freely diffuses into dead cells with compromised membranes and forms covalent sulfhydryl bonds with intracellular protein nucleophiles. Cisplatin is commonly applied for live-dead discrimination in mass cytometry because it a) binds covalently to cellular proteins within cells and b) stains cell membranes of viability-compromised cells to a much greater extent than live cells (54) . The first property allows cisplatin to remain bound through multiple downstream staining steps used in mass cytometry protocols, and the second property is leveraged for live-dead cell discrimination. Of note, cisplatin needs to be titrated on the sample as diffusion efficiencies vary by tissue (156) . This is followed by surface staining (36, 157) , after which many parts can be added into staining protocols, such as intracellular (36, 158, 159) , intranuclear (160, 161) , phosphostaining (141) (142) (143) , or tetramers (160) . A cryopreservation method has been developed for long-term (9 months) storage of antibody mixes. Use of such a method is highly advisable in large mass cytometry studies with multiple staining cycles to enhance staining reproducibility and data consistency across time (162) . The general order of staining steps used by our group is summarized in Figure 3 . For phospho-protein studies, a dephosphorylation block with pervanadate is inserted at the end of step III (Figure 3) . The cells are fixed with cold methanol after step IX (intracellular/ intranuclear staining) before proceeding to the phospho-staining (step XI) and the subsequent DNA intercalator (step XII) (163) . Tetramer staining (e.g., MHC or CD1 tetramers to identify cells with unique TCR specificities) can be incorporated into surface staining (160) . Multiple wash steps are incorporated after each staining step to thoroughly remove contamination from unbound reagents and minimize background. As already mentioned, in mass cytometry cells are separated from debris by addition of a DNA intercalator staining incorporated at the end of all staining steps (step XII, Figure 3 ). Natural abundance Iridium ( 191 Ir and 193 Ir) will bind to nucleic acid after cell membranes are permeabilized and the detection of both Ir isotopes allows single cells to be distinguished from debris and doublets (see Manual Gating) (7, 164) . Samples are finally washed once with cell staining buffer and counted. Cell numbers need to be kept consistent throughout the experiment and samples need to be resuspended at the optimal cell concentration for sample acquisition, improving data quality and minimizing doublets. The optimal cell concentration depends upon the specific mass cytometry instrument model and its injector type. Each injector has an optimal and maximum event rate, and event rates higher than the maximum indicated for that injector will drastically increase doublets in the resulting data. After any paraformaldehyde (PFA) fixation step, sample staining can be conveniently paused and resumed later by simply storing the sample overnight in 1.6-2% PFA at 4°C. Unlike the case with denaturation of certain conventional fluorophores, metal isotopes are not affected by prolonged fixation steps. If performing phospho-staining, a pause stop can be incorporated with 100% methanol at -80°C. When opting for completing the entire staining protocol in one workday and acquiring the samples on CyTOF ® on a subsequent day, samples can be stored in DNA intercalator (e.g. 191 Ir/ 193 Ir cocktail) overnight at 4°C. When storage needs to occur longer than overnight at either of these 2 steps, samples can be washed once with 1 mL cell staining buffer and pellets stored up to one week (165) . If stained samples need to be stored longer, Sumatoh et al. have developed a method to preserve them in 10% dimethylsulfoxide (DMSO) + 90% FBS at -80°C (165). This long-term storage method was tested on human PBMCs only and therefore may not be applicable for other types of tissues. Our group has confirmed that this protocol also works well for long-term storage of stained mouse splenocytes (data not shown). Sample acquisition on a CyTOF ® mass cytometer is performed in one of 2 modes: 1) solution mode or 2) event mode for beads and single-cell analyses. Solution mode is used for testing buffers for contaminants, verifying antibody conjugations, or re-testing older conjugated antibodies (> 6 months). A buffer or conjugated antibody is diluted 10,000 times into either ddH 2 O or CAS and acquired in solution mode. The data output of this type of acquisition includes a mass spectrum showing each element detected as a peak, accompanied by average ion counts. This method does not give the exact number of metal ions conjugated to each antibody. To generate such data, conjugated antibodies must be diluted in Tuning Solution SKU 201072 ® (Fluidigm, San Francisco, CA), which contains 6 metal isotopes at known concentrations and is used for daily calibration of the instrument. Fluidigm provides users with the necessary protocols and worksheets upon request. Immediately prior to data acquisition, at least 2 washes with ddH 2 O are required to remove residual buffer salts; these otherwise accumulate at the detector and can cause detrimental instrument drift. Tuning, the process by which the CyTOF ® instrument is calibrated before any sample acquisition, maximizes the signal intensity of the metal isotopes within the optimum range from atomic mass 153-176 while minimizing isotope oxide formation (M+16) by the inductively coupled plasma (ICP). The acquisition instrument requires re-tuning every 6 hours during prolonged data acquisition. For more detail on machine tuning and sample acquisition, we refer the reader to helpful Fluidigm machine manuals and a video by Leipold and Maecker (51, 166, 167) . The fluid tubing of a mass cytometer is of much smaller diameter than that of conventional flow cytometers (51) . This makes the CyTOF ® relatively more prone to obstruction by accumulated debris. For this reason, samples need to be passed over a double 30 mm cell strainer (Partec North America Inc., Swedesboro, NJ) immediately prior to sample acquisition. Even following filtration, there remains a periodic need to unclog the lines during sample runs. Our lab routinely uses BD FACSAria ® 50 mm sample inline filters (BD Biosciences, San Jose, CA), which perfectly fit the sample probe line of the Helios ® . When using sample barcoding, updating the CyTOF ® software with separate labels for all barcodes is strongly recommended throughout the sample run [e.g., for CD45 live barcoding, rather than a uniform label of "CD45", it is preferable to use "CD45_102" (isotope mass), "CD45_healthy" (experimental group), or "CD45_tumor" (tissue type)]. The importance of this labeling schema is clarified in "Debarcoding". Importantly, when samples sit in ddH 2 O or CAS, staining intensity and quality declines over time (168, 169) Therefore, samples should be resuspended in these media only immediately prior to data acquisition. Whether ddH 2 O or CAS is the preferred solution for acquisition depends on both the CyTOF ® model and injector type associated with that model. For the ball joint injector on CyTOF2 ® and the narrower HT injector on early Helios ® machines [the narrow bore injector (NB)], cell pellets must be resuspended in ddH 2 O. The Helios ® injector has a narrower inner diameter, resulting in smaller ion clouds (roughly one half the size of those attained with wide-bore injectors), reducing doublets and doubling the machine's event acquisition rate (cells/second) (51) . However, unintended consequences to the data such as lower median signal intensities and higher coefficients of variation (CVs) have been described using these injectors (170) . The Helios ® system introduced the wide-bore (WB) injector and CAS to address this pitfall. The WB injector has an inner diameter intermediate between the (narrower) HT injector and the (wider) CyTOF ® injector. The CAS has a higher ionic content than water and in combination with the newer WB injector resolves the data quality issues seen with narrow-bore injectors in Helios ® datasets (170) . One disadvantage of the newer WB injector, however, is that it drastically decreases the maximum event rate from 500 event/sec to 250 events/sec, requiring significantly longer sample acquisition/run times. During mass cytometry runs, signal drift occurs over time due to gradual accretion of cellular material inside the instrument and associated progressive acquisition delays (51) . Signal drift can cause variations between and even within individual data files. Since the sample acquisition speed on a CyTOF ® is quite low compared to conventional flow cytometers, experiments are often run across multiple days. The instrument is tuned each day and may be cleaned periodically within a prolonged experiment run, causing additional variance between days. Additionally, in consortia or clinical trials, data is collected at multiple sites and on multiple instruments. To correct for these variances and minimize measurement variations, samples are resuspended in EQ Four Element Calibration Beads ® (EQ 4 beads) (Fluidigm, San Francisco). EQ4 beads are polystyrene bead standards containing known relative quantities of metal isotopes from 4 metals ( 140 Ce, 151 Eu, 153 Eu, 165 Ho and 175 Lu), diluted 10-fold in either ddH 2 O or CAS and used for normalizing data within and across experiments (see Normalization & Concatenation) (171, 172) . This method for minimizing experimental variation is efficient enough that data acquired from different machines can be combined into cumulative datasets. This is important, since different machines have been shown to have discrepancies in their atomic mass sensitivity ranges (173, 174) . Recently, Liu et al. further improved the isotopic range of EQ 4 beads by adding 3 elements ( 89 Y, 115 In, and 209 Bi), resulting in a total of 7 isotopes. The authors demonstrated that the EQ4 beads did not consistently normalize isotope signals outside the mass range of 140-175, whereas the 7-element calibration bead system resolved this issue (175) . As of the time of this manuscript publication, 7-element beads are not available for purchase and interested parties must synthesize these themselves. The output data from the machine is in the form of a Flow Cytometry Standard (.fcs) file (FCS file). The FCS file structure is a standardized array with columns representing channels and rows representing events. This is used for downstream data analysis through programming languages (e.g., R) or FCS fileprocessing platforms [e.g., FCS Express ® (DeNovo Software, Pasadena, CA), Cytobank ® (Cytobank Inc, Santa Clara, CA), or FlowJo ® (FlowJo LLC, Ashland, OR)]. Generally required items are: 1) the generated FCS file; 2) a panel file listing the metals and corresponding markers; and 3) a metadata file listing the names of all FCS files used, categorizing FCS files into experimental groups (e.g., control, knockout, etc.). The panel and metadata files can be in Comma Separated Values file (CSV) or Microsoft Excel (XSLX) format. Refer to the analysis program/ package analysis instructions for guidance in selecting file format. To correct machine signal variation and thereby minimize measurement variation, "global normalization" (normalization across all FCS files) must be performed. Though this normalization is performed computationally after sample acquisition, it requires the use of bead-based standards added to the sample before sample acquisition and therefore merits advance consideration. The two algorithms for normalization are 1) MATLAB ®based bead normalization shown by Finck et al. (171) and 2) Fluidigm's bead identification and normalization (172) . The primary difference between these is that MATLAB ® normalization compares each file to the other files in the same data set, whereas Fluidigm's algorithm compares the acquisition files to a set of external values. One pitfall of the Fluidigm algorithm is that the data may be slightly reshaped to fit the external values used as compared to data normalized by the MATLAB ® method (176) . In some circumstances, a single sample must be acquired across multiple FCS files and subsequently recombined into a single data file prior to analysis. Common scenarios necessitating such an approach include: 1) a clog occurs in the sample line, sample capillary or the nebulizer requiring the acquisition to be halted and restarted 2) a particularly large sample requiring collection from multiple tubes (168, 169) . The recombining of these files into a single FCS file before analysis is termed concatenation (85) . Concatenation should be performed in order of file acquisition, to avoid introducing errors from incorrect sequencing (86) . Because events are normalized independently, the order of normalization and concatenation should not affect data quality. However, many algorithms remove artifacts and metadata during the normalization process, which can cause errors during concatenation. For this reason, unless the sample uses a large number of bead standards (> 200 hundred beads), it is advisable to first concatenate and subsequently normalize. Several different algorithms exist to aid normalization and concatenation. For the MATLAB ® algorithm, R packages such as premessa or CATALYST can be used (87, 88) . CATALYST fuses all FCS files into a single cell experiment object in R. Therefore if using CATALYST, concatenation is not required. Fluidigm's algorithmic bead normalization is performed through the instrument's software. (As of this publication, the software version is CyTOF ® 6.7). Note that each batch of EQ 4 beads has its own corresponding passport of external values, which must be updated in the CyTOF ® software whenever a new batch of EQ 4 beads is used. Though doublets have historically been considered noise or unwanted confounders, Burel et al. found that the molecular signature of T-cells in complex with other immunological cells such as monocytes or B-cells showed phenotypic differences when compared to unattached single T-cells. Moreover, they observed an increase of T-B-cell complexes in the circulation of type-1 diabetes patients (89) . The implication of these findings is that removal of doublets may in certain studies skew the data in ways that miss important biological outcomes. However, since most doublets are indeed unwanted noise, it is generally advised to study doublets separately from the singlet data analysis to determine their value to any individual experiment, to avoid compromising experimental rigor. A good QC step to visualize the effect of normalization is to plot the bead intensity over time before and after normalization (Figure 4) . In barcoding, each sample is labeled with a unique identifier and samples are mixed in a single tube before data acquisition. The output is a single FCS file consisting of the multiplexed data from all of the barcoded samples in the tube. Therefore, an extra step must be taken following data acquisition to split the output data from each barcode into its own FCS file. This process is termed "debarcoding". In debarcoding, events are stored in new FCS files corresponding to the individual barcodes, identified by the debarcoding algorithm. The number of output FCS files is equal to number of barcoded samples used. Samples are reassigned a unique identifier in the debarcoding process. During debarcoding, the single-cell deconvolution algorithm must define which channels are positive and negative for each cell. For this, a threshold value is selected, and cells above this threshold value will be assigned to their corresponding barcode sample (90) . The threshold can be selected manually or with the help of an algorithm (e.g. CATALYST provides estimation of threshold values). Rarely, high stringency in such parameters can reduce overall event number in the dataset. The vast majority of unassigned events are cell-cell doublets or debris (33) . Occasionally, cells can be sorted into incorrect barcode channels (false assignment). This false assignment rate is generally under 0.5% under uniform sample staining/ experimental conditions. A valuable internal control is to leave 1 barcode from the scheme unused, allowing an estimation of the false assignment rate for debarcoding. Both premessa and CATALYST packages in R can be used to debarcode (88, 91) . Debarcoding can also be done using the MATLAB Compiler Runtime (33) . Graphical User Interfaces (GUIs) are available for these processes. GUIs are preprogrammed interfaces that run through the code, allowing analysis without directly interacting with the code. Most of these algorithms require a debarcoding scheme alongside the FCS file, typically in the form of a table that correlates the sample identity to the element mass of the barcode isotope. This instructs the algorithm which barcode corresponds to a sample's FCS file. Formatting varies by the debarcoding resource. After normalization, concatenation and debarcoding, compensation is generally the final data pre-processing step in mass cytometry workflow. Compensation is the process whereby detection signal spillover is resolved. In conventional flow cytometry, this is a routine practice due to overlapping excitation and emission spectra of fluorophores and spillover which correlates linearly with fluorophore signal intensity (92) . Compensation subtracts the percentage of a fluorophore's spillover from the measured signal in that particular channel. Mass cytometry can have spillover, but antibody titrations and intelligent panel design can optimize the signal-to-noise ratio such that compensation is no longer necessary (67). Chevrier et al. performed a detailed analysis and found that mass cytometry spillover has a similar linear relationship with the primary signal (84) . This relationship has been applied to generate a spillover matrix allowing compensation to be performed with the CATALYST package either in R or via a web application created by the same group (87, 93) . Another helpful algorithm is CytoSpill, a statistical program that aims to minimize spillover effects. Unlike the CATALYST algorithm or conventional flow cytometry analysis programs, this algorithm does not require the use of single-cell controls, making the process significantly cheaper and easier (94, 95) . During the experiment, samples have been put through various stresses which results in the events in the FCS file being a combination of a) live cells, b) dead cells, c) cell-cell doublets, d) beads e) cell-bead doublets and f) bead-bead doublets. The goal of this step is to purge the data as much as possible of b-f, such that the output consists of live single cells. Manual gating can be performed using many established platforms from conventional flow cytometry, such as FlowJo ® or FCS Express ® . Though usually done via manual gating, some algorithmic methods such as that in the CATALYST pipeline can efficiently remove d-f. In Figure 5A , a standard mass cytometry gating strategy is shown. Here, the "event length" parameter is used in conjunction with a DNA intercalator to select single cells, followed by live/ dead gating based upon cisplatin intercalation. Though singlecell gating using this method has proven difficult, the Helios ® mass cytometer is now equipped with additional Gaussian parameters (e.g. center, offset, width, residual), which can be used to streamline data (96) . These parameters are generated in the TOF chamber. The ions enter the TOF chamber via a narrow slit ( Figure 5B ). Every 13 µs (exact interval varying by machine), a high-voltage pulse provides equal energy to all ions that have accumulated in the chamber during the interval, accelerating them across the TOF chamber and onto the detector. This highvoltage pulse is a push. As the energy provided to the ions is uniform, velocity varies by mass and ions with greater mass require longer times to reach the detector. Once these ions hit the detector, they generate an electronic pulse which when plotted and mathematically smoothened takes the shape of a bell curve (Gaussian distribution) from which the Gaussian discrimination parameters (center, offset, and width) are extracted ( Figure 5B ) (51, 97) . The Residual parameter is extracted by calculating the difference between the actual electronic pulse plot and the smoothened bell curve (97) . The improved gating strategy is shown in Figure 5C . Once the runs are pre-processed, all events are expected to represent live single cells. Before proceeding, it is worthwhile to visualize the number of cells in each sample. This can be done in R with a simple bar graph visualizing the cell number within In the Helios ® machine, the Gaussian parameters are generated in the Time-Of-Flight chamber where ions are separated on mass-to-charge ratio. Once these ions encounter the detector, an electronic pulse is generated, which is mathematically smoothened into a bell curve. Center, Offset, and Width can be extracted from this curve. (C) The added Gaussian parameters greatly improve the ability to select a single-cell population. Note that if beads are not excluded through the normalization process, they must be gated out manually before gating for single live cells. each sample and comparing across samples ( Figure 6A ). If the cell count is insufficient and the population of interest is a very rare or marginal proportion of the total cells (expected to vary by experimental conditions or cell/tissue type), the experiment may have to be rerun. Even if sufficient cells exist in the input samples, if the difference in cell numbers between runs is extremely large, such variability can affect the statistical analysis. In such cases, the cells can be randomly subsampled (randomly selected from the experimental samples based on parameter matching to controls) or the experiment repeated ensuring larger cell counts. The next step is to visualize relationships between samples. One cumulative data visualization method is Multidimensional Scaling (MDS) (98) . In the sample MDS plot showing data for murine bone-marrow derived dendritic cells (BMDCs) ( Figure 6B, upper panel) , the controls and comparisons form distinct groups, indicating that wild-type (WT) mice are similar to one another and knockout (KO) mice are similar to one another, but that there is a significant difference between WT and KO mice. In the lower MDS plot ( Figure 6B, lower panel) , while the WT are similar and form a distinct group, the KO mice are dissimilar. Furthermore, one KO clusters with the WT controls. This should prompt evaluation for experimental errors requiring re-analysis. Errors can arise from variances in staining, runs on distinct days, and human error, which can be challenging to troubleshoot. For this reason, we recommend adding a barcoded internal control in each tube as a quality control ("spiking"). Errors can be identified by variation of internal controls from expected outcomes. After critically assessing cell count and sample quality, one may proceed to analyze the specific markers for the experiment, beginning with a Non-Redundancy Score plot (NRS plot) ( Figure 6C ) (99) . This graph visually represents the variation across all samples by individual marker. Ideally, the marker variation is caused by the biologic differences between experimental groups. However, large variations in cell count or in staining intensity can also be sources of variation on the NRS plot. Markers that contribute minimal or no variation between samples can be excluded from clustering. To view the marker expression in more detail, a histogram can be used to visualize the results ( Figure 6D ). Marker expression histograms are extremely useful when plotting an internal control sample spiked into a tube (by barcoding). This provides a quality control for staining reproducibility. In Figure 6E , the internal controls are color-coded by day of staining, demonstrating that reproducibility was not optimal. After bead normalization, data must be transformed to allow proper distinction between cell populations positive and negative for each marker. The most common method is the arcsinh transformation, using a cofactor of 5 (99) . In R, CATALYST packages have functions that can perform this transformation (99, 100) . This transformation can also be accomplished directly in Cytobank ® or FlowJo ® ; however, this is lost when exporting. The next step is data clustering based on group similarities between cells. This data visualization method compartmentalizes events (cells) into groups which are then used to insure the detection, characterization, and calculation of relative abundance of the different populations in the sample. Note that not all clusters/groups represent a cell population, as some clusters might represent the same population with slightly differing marker expression. For this reason, it is important to carefully study and compare the marker expression within each cluster. Several different clustering algorithms are available for mass cytometry data, differing primarily based on the criteria they use to assign a cluster designation. The highest performing algorithms are FlowSOM, X-Shift and PhenoGraph (101) . If specifically looking for rare populations, X-Shift is preferable due to its statistical power to detect rare events (102) . While the algorithms themselves use a variety of models to cluster the data, they utilize the marker expression within each individual cell to distinguish the clusters based on phenotype. Therefore, it is important to determine in advance which cells to isolate and the potential cell types present in the samples, and design the experimental panel accordingly. The more distinguishing markers in the panel (markers that can be used to show a clear and meaningful difference in the subpopulations in the sample), the better the quality of the clustering will be to allow clear differentiation and ease of population identification (103) . As explained in detail in the previous section (Understanding Data Quality), an NRS plot ( Figure 6C ) is utilized in order to select the markers most relevant for clustering; meaning markers with a low non-redundancy score will not contribute to the cluster outcome and can therefore be excluded from clustering analysis. The clustering algorithms use the data on event phenotype (using the markers selected by the analyst) to group events into clusters or "bins". Clustering is generally performed over all FCS files to allow sample/file comparability in downstream analyses. R or Python are the recommended clustering programs, as these allow much greater control over the statistical binning process. Packages such as CATALYST or Phenograph can be used, the details of which can be found in their respective workflows (99, 104) . Note that apparently different clustering packages or software might ultimately utilize similar algorithms. For example, CATALYST uses the FlowSOM algorithm for clustering (105, 106) . Thus, it is important to study in advance the workflow being used and identify which algorithms may best suit the output data, well in advance of the experimental run. While some algorithms such as PhenoGraph will automatically select the optimum number of clusters to use, others like FlowSOM will require manually selection of the cluster number (k) (104, 105) . One method for manual selection is initial generation of a large number of clusters (~100), followed by serial merging of particular clusters until the optimum cluster number is obtained. Again, not every cluster represents a subpopulation in the sample, and it is here that applying biological understanding is important. Whether these clusters represent a novel/separate subpopulation (creating a new cluster) or are similar to each other and therefore will need to be merged must be determined based on the marker expression levels on the different clusters. All clusters that represent a will be continually merged until only clusters that are determined to be biologically different from each other remain. The same process can be performed statistically with the help of a delta plot, where the relative change under the Cumulative Distribution Function (CDF) is plotted for each cluster value. It is no longer optimal to increase the number of clusters selected (k) when doing so results in no/minimal change in the CDF ( Figure 7A) . A delta plot is generated using either the 'deltaPlotR' Package or the CATALYST package (107, 108) . When clustering initiates, a randomly generated 'seed value' is used as the starting point of the clusters. To ensure that clustering is reproducible, one should record the seed value used to initiate clustering (177) . To ensure that clustering is robust and to avoid the possibility of error due to random selection of seed value, clustering should be repeated with at least 3 random seed values and output clusters compared. These should approximate each other closely. A good QC measure to ensure clustering has been correctly performed is to include markers for a known control cell population within the sample. When opting for internal controls, clustering could be performed on the FCS files separately in order to verify similar cluster proportions between tubes. The next step is to determine whether clusters represent distinct populations and to phenotype each cluster based on statistical analysis and graphical outputs. A combination of biological and statistical expertise become highly relevant at this stage. R packages like FlowSOM, FlowCore and CATALYST offer several different visualization functions. The most common visualization approaches are heatmaps. As seen in Figure 7B , the expression for each marker is aggregated by cluster. These clusters were generated from FCS files and therefore can be applied to any group used in the analysis (in this case, WT murine controls and knockouts). Heatmaps not only help to visualize markers expressed in different clusters on the merged FCS files, but also to estimate the relative expression intensity. Marker expression levels determine to which immune lineages these clusters may belong, elucidate novel cell clusters or phenotypes, and can aid in determining whether certain clusters can be merged. The heatmap parameters can be adjusted depending on use case (e.g., comparing a particular marker's expression across groups rather than comparing prevalence of clusters, etc) (178, 179) . Mass cytometry data is of necessity multi-dimensional. Clusters, marker intensity, and differing markers complicate interpretation. To address this, dimensionality reduction algorithms can be applied to concisely visualize multiparametric data in two dimensions. t-distributed Stochastic Neighbor Embedding (tSNE), Principal Component Analysis (PCA), Isometric Feature Mapping (Isomap), and Uniform Manifold Approximation and Projection (UMAP) are all common methods for high-dimensional data visualization (99, (180) (181) (182) (183) . These facilitate easy identification of the experimental group differences by visually organizing populations using distance on the plot as a surrogate for similarity between groups. The most commonly used approaches are tSNE and UMAP. UMAP, unlike tSNE, is primarily nonparametric nonlinear dimensionality reduction algorithm originally designed to preserve global data structure with smaller data sets, though more recently parametric UMAP was demonstrated to compare favorably in performance to its onparametric counterpart while adding the benefit of learned mapping of new data. In UMAP, distances both within and between clusters are generally meaningful; whereas in tSNE, distances between clusters may not be. Therefore, UMAP is the preferred mode of visualization (181, 184) . In Figures 7C, D , a case is presented in which both tSNE and UMAP efficiently visualize the cell clusters present in the WT control and KO mouse groups. Figures have been generated in R via the CATALYST package using their standard settings. One can also use Flow Core in R, or directly generate such visualizations through Cytobank. Once clusters have been visualized and characterized, they can be annotated in R. Clustering can require hours to days to complete depending on the dimensionality reduction algorithms and the acquired total cell count per group. To prevent prolonged run times, these algorithms are ideally run on a subsample of the total population (~2000 cells). When looking for a particularly rare population, more cells can be selected for algorithm application. Because subsampling is random, a set seed function in R is used to ensure reproducibility. To ensure clusters are not lost in the sample and that subsamples are representative of the total sample, at least 3 iterations should be performed with 3 different seed values to confirm similar results (185) . Another method to visualize the relationship between clusters is the diffusion map. Clusters are arranged in order of overall similarity, creating the appearance of a gradient based on change in marker intensity per analyzed event. Diffusion maps are commonly used to study the differentiation or lineage origin of individual cells in a sample set. Another useful application of this visualization technique is for groups of cells of the same type, taken at different time points in treatment or maturation of a sample (186) . When viewing diffusion maps, it is critical to remember that association does not imply causality. Viewing 2 adjacent clusters in a diffusion map indicates that they are the most similar in the sample, but not necessarily that one has differentiated from the other. Data analytics must be tempered with biological understanding of the cell type. An algorithm commonly used for this type of analysis is Wanderlust (187) . Wanderlust assumes a linear trajectory and has significant limitations with complex datasets or in the case of multiple cell fates. For those datasets with potential bifurcating branches, Wishbone or Monocle 2 can be applied (188, 189) . There are several other resources that can be applied in place of or in addition to R coding such as FlowJo ® , FCSExpress ® , and Cytobank ® (190) (191) (192) . Several R scripts also offer GUIs making it significantly easier to analyze such data (e.g., Cytofkit package for Phenograph). However, R allows the maximum ability to finetune parameters, making it extremely attractive as an analysis platform (193) . Finally, Astrolabe (Astrolabe Diagnostics, Fort Lee, NJ) offers an interactive platform in which the researchers can interact directly with a computational biologist to analyze data (194) . Recently, automatic annotation algorithms have started to be developed [e.g., Automated Cell-type Discovery and Classification (ACDC)], which allows users to input biological parameters such as marker expression levels and apply this to identify such populations within the sample (195) . Currently, the sensitivity of metal isotope-tagged antibodies is lower than that of the most quantum-efficient fluorophores, such as phycoerythrin (PE). The main reason for this is that chelating polymers used for antibody conjugation can only accommodate a maximum of 100 metal ions, creating a ceiling on signal intensities and making it more difficult to measure extremely weakly expressed markers (20, 21) . Another major drawback of mass cytometry is the much lower acquisition flow rate as compared to flow cytometry. This is due to the dynamics in the TOF chamber, resulting in longer acquisition times. Where conventional flow cytometers can have a flow rate of up to 50,000 events/second, CyTOF2 ® only has a flow rate of 500 events/ second (250 events/second for Helios ® ) (167, 196) . Despite these low flow rates, mass cytometry machines are also quite prone to clogging due to the small diameters of the sample lines and nebulizer. The latter is easily resolved by passing samples through a cell strainer both before and during sample acquisition as explained in "Sample Acquisition & Data Output". The latest CyTOF XT ® machine from Fluidigm is capable of sensing and removing clogs automatically. In addition, sample preparation and staining require extra caution with regards to possible heavy metal contaminants such as barium from dish soap or lead from distilled water. Another limitation is the high cost of metal-tagged antibodies, conjugation kits, and other reagents. Although the costs are significantly higher as compared to flow cytometry, these are still much lower than for scRNAseq (21) . With increasing number of users and advancements in reagents and technology, these costs are expected to decrease. Finally, because cells are vaporized during the analysis, sorting populations of interest is not currently possible with mass cytometry. To address the latter, one can sort populations on a fluorescent-activated cell sorter for follow-up studies. The increased number of parameters made available by mass cytometry complicates data analysis, requiring improved bioinformatics approaches for accurate interpretation and visualization of fcs files. Data analytics requires a combination of biological, statistical, and programming knowledge, making the overall process tedious for single individuals. Though packages like ParkerICI/premessa and CATALYST in R are designed to make analysis significantly more user-friendly (88, 93), without significant R programming skills, the researcher remains largely restricted to package workflow and graphing programs. More recently new options have arisen including GUIs and web interfaces that allow analysis without necessitating R programming skills. Most notable among these is Cytobank ® (192) . Another potential limitation during manual cluster analysis is the potential for operator bias. This can be significantly mitigated by automating cluster identification through algorithms such as CITRUS (Cluster Identification, Characterization, and Regression) (197, 198) . CITRUS automatically clusters cells in each experimental group and compares them based on inputted biological parameters as well as data from publicly available datasets. Despite rapid technological advancements and vast capacities, mass cytometry has many areas ripe for future development. At the instrument level, increasing the flow rate (reducing time of analysis) without losing data quality or altering machine sensitivity are key needs. The newest mass spectrometers, have a more sensitive detector and/or improved TOF chambers than prior Helios ® mass cytometers (51, 199, 200) . Novel polymer chelators able to accommodate more than 100 metal ions at a time could increase signal intensities significantly and thus empower investigations of very weakly expressed molecules (tertiary antigens). Additional isotopes and new methods to purify existing isotopes at a requisite (>98%) purity and novel conjugation chemistries will allow mass cytometry to approach the theoretical 120-parameter capabilities of the Helios ® (Figure 2A ). With the recent development of compensation tools, stringency requirements for metal isotope purity may be relaxed (84, 95) . By merging multiple mass cytometry panels into a single clustering analysis, the computational tool CyTOFmerge is greatly expanding the number of simultaneously analyzable parameters (201) . Annotation is one of the most difficult portions of any exploratory mass cytometry analysis. It requires a depth of biological knowledge regarding the unknown cell population, which the researcher may lack at the time of the experiment. One solution is to compare the experimental dataset to available reference datasets. This approach is seen in techniques such as Seurat, a pipeline for scRNAseq. These reference maps provide information about cellular phenotype and interactions, creating a framework to compare existing data with unknown cell populations to aid in their identification and classification. Scaffold is another example of a system for immune cells, built on flow and mass cytometry data sets. Unlike for genomic data, such reference maps are scarce for mass cytometry and represent invaluable resources for future development (202) . Another much-needed application in mass cytometry is a method similar to pseudotemporal ordering used in scRNAseq. This would support immune lineage differentiation studies, allowing clusters originating from common lineages to be identified and mapped. However, distinct from pseudotemporal ordering in scRNAseq, which uses transcription data for its analysis, this would require the use of lineage markers and other characteristic surface proteins (203) . Ultimately, the greatest limitation to more widespread application of mass cytometry is the difficulty in post-acquisition data analysis and interpretation. Automatic screening and comparison of mass cytometry output variables through databases and existing publications represents a major step forward in understanding and annotating cells. Artificial Intelligence (AI) systems are likely to become standards in the near future to make increasingly accurate predictions and annotations based on the ever-increasing available data (204, 205) . Automating the entire pipeline into a one-click analysis, with selectable output visualizations and analysis factors, will make the technique more accessible for wider applications of mass cytometry in such key expanding areas as TME responses and immunotherapies. AI, AH, and AP performed the literature compilation, writing, editing, and proofreading. All authors contributed to the article and approved the submitted version. Funding support to the authors was provided by 1R01HL133462 (NIH) (AH, AP) and the Batchelor Foundation for Pediatric Research (AH, AP). Instrument and equipment support was provided by the Sylvester Comprehensive Cancer Center Flow Cytometry Shared Resource, supported by grant P30CA240139 (National Cancer Institute). The mass cytometry equipment was supported by a State of Florida Bankhead-Coley Grant Program award 8BC09. CyTOF ® is a trademark ( ™ ) and/or registered trademark ( ® ) of Fluidigm Corporation or its affiliates in the United States and/or other countries. Mass Cytometry: Technique for Real Time Single Cell Multitarget Immunoassay Based on Inductively Coupled Plasma Time-of-Flight Mass Spectrometry Mass Cytometry for the Assessment of Immune Reconstitution After Hematopoietic Stem Cell Transplantation Mass Cytometry Reveals Immune Signatures Associated With Cytomegalovirus (CMV) Control in Recipients of Allogeneic Haemopoietic Stem Cell Transplant and CMV-Specific T Cells Integrated Immune Signature Analyses Identifies Evolution of Distinct Immunoregulatory Cell Populations Which Control Alloreactivity After Allogeneic HSCT Maturation and Phenotypic Heterogeneity of Human CD4+ Regulatory T Cells From Birth to Adulthood and After Allogeneic Stem Cell Transplantation Low-Dose IL-2 Selectively Activates Subsets of CD4+ Tregs and NK Cells Validation of CyTOF Against Flow Cytometry for Immunological Studies and Monitoring of Human Cancer Clinical Trials Data-Driven Phenotypic Dissection of AML Reveals Progenitor-Like Cells That Correlate With Prognosis High-Dimensional Single-Cell Analysis Predicts Response to Anti-PD-1 Immunotherapy Phenotyping of Tumor Infiltrating Immune Cells Using Mass-Cytometry (CyTOF) Innate Immune Landscape in Early Lung Adenocarcinoma by Paired Single-Cell Analyses Commonly Occurring Cell Subsets in High-Grade Serous Ovarian Tumors Identified by Single-Cell Mass Cytometry Mass Cytometry Defines Distinct Immune Profile in Germinal Center B-Cell Lymphomas PD-1 Marks Dysfunctional Regulatory T Cells in Malignant Gliomas Immuno-Phenotyping of Tumor-Specific CD8 T Cells Using High-Dimensional Mass Cytometry Comprehensive Immune Monitoring of Clinical Trials to Advance Human Immunotherapy Single-Cell Mass Cytometry Identifies Mechanisms of R esi stance t o Im munoth erapy i n AML Spectral Flow Cytometry OMIP-069: Forty-Color Full Spectrum Flow Cytometry Panel for Deep Immunophenotyping of Major Cell Subsets in Human Peripheral Blood A Deep Profiler's Guide to Cytometry Mass Cytometry: Single Cells, Many Features Potentiality of Multiple Modalities for Single-Cell Analyses to Evaluate the Tumor Microenvironment in Clinical Specimens Flow Cytometry Controls, Instrument Setup, and the Determination of Positivity Selecting Fluorochrome Conjugates for Maximum Sensitivity Novel Full-Spectral Flow Cytometry With Multiple Spectrally-Adjacent Fluorescent Proteins and Fluorochromes and Visualization of In Vivo Cellular Movement Guidelines for the Use of Flow Cytometry and Cell Sorting in Immunological Studies Rigor and Reproducibility of Cytometry Practices for Immuno-Oncology: A Multifaceted Challenge The Anatomy of Single Cell Mass Cytometry Data Metal-Isotope-Tagged Monoclonal Antibodies for High-Dimensional Mass Cytometry Scalable Conjugation and Characterization of Immunoglobulins With Stable Mass Isotope Reporters for Single-Cell Mass Cytometry Analysis Metal-Isotope-Tagged Monoclonal Antibodies for High-Dimensional Mass Cytometry Palladium-Based Mass Tag Cell Barcoding With a Doublet-Filtering Scheme and Single-Cell Deconvolution Algorithm Platinum-Conjugated Antibodies for Application in Mass Cytometry Atomic Mass Tag of Bismuth-209 for Increasing the Immunoassay Multiplexing Capacity of Mass Cytometry Multiparameter Phenotyping of Human PBMCs Using Mass Cytometry Qdot Probes Technology Overview -Us (Accessed cited Silver Nanoparticles for the Detection of Cell Surface Antigens in Mass Cytometry Quantitative Measurement of Cell-Nanoparticle Interactions Using Mass Cytometry Tantalum Oxide Nanoparticle-Based Mass Tag for Mass Cytometry A Universal Live Cell Barcoding-Platform for Multiplexed Human Single Cell Analysis Sample Preparation for Mass Cytometry Analysis Standardization and Quality Control for High-Dimensional Mass Cytometry Studies of Human Samples Flow Cytometry Controls Multiplexed Mass Cytometry Profiling of Cellular States Perturbed by Small-Molecule Regulators Transient Partial Permeabilization With Saponin Enables Cellular Barcoding Prior to Surface Marker Staining Barcoding of Live Human Peripheral Blood Mononuclear Cells for Multiplexed Mass Cytometry Extended Live-Cell Barcoding Approach for Multiplexed Mass Cytometry Rapid Monoisotopic Cisplatin Based Barcoding for Multiplexed Mass Cytometry Tellurium-Based Mass Cytometry Barcode for Live and Fixed Cells Enhanced Multiplexing in Mass Cytometry Using Osmium and Ruthenium Tetroxide Species Method for the Analysis of the Tumor Microenvironment by Mass Cytometry: Application to Chronic Lymphocytic Leukemia Modern Flow Cytometry: A Practical Approach Flow Cytometry Quality Requirements for Monitoring of Minimal Disease in Plasma Cell Myeloma Flow Cytometry: A Blessing and a Curse Atlas of the Immune Cell Repertoire in Mouse Atherosclerosis Defined by Single-Cell RNA-Sequencing and Mass Cytometry Highly Multiplexed Imaging of Tumor Tissues With Subcellular Resolution by Mass Cytometry Multiplexed Ion Beam Imaging of Human Breast Tumors How It Works | MIBI Multiplexed Ion Beam Imaging Technology Protocols for MIBI Multiplexed Tissue Imaging Mass Cytometry Panel Optimization Through the Designed Distribution of Signal Interference MACS ® Marker Screen, Human, Version 02 | Phenotyping Assays | Kits and Support Reagents | MACS Flow Cytometry | Products | Miltenyi Biotec | USA . Available at Development of a Comprehensive Antibody Staining Database Using a Standardized Analytics Pipeline Multiplexed Peptide-MHC Tetramer Staining With Mass Cytometry Development of Mass Cytometry Methods for Bacterial Discrimination Identification of Hypoxic Cells Using an Organotellurium Tag Compatible With Mass Cytometry High-Throughput Quantitation of Inorganic Nanoparticle Biodistribution at the Single-Cell Level Using Mass Cytometry Highly Multiplexed Simultaneous Detection of RNAs and Proteins in Single Cells Simultaneous Detection of Protein and mRNA in Jurkat and KG-1a Cells by Mass Cytometry Development of a Click-Chemistry Reagent Compatible With Mass Cytometry Tellurium-Containing Phenylalanine Mimic, Allows Monitoring of Protein Synthesis In Vivo With Mass Cytometry Multiplexed Probing of Proteolytic Enzymes Using Mass Cytometry-Compatible Activity-Based Probes Cell Size Assays for Mass Cytometry Proliferation Tracing With Single-Cell Mass Cytometry Optimizes Generation of Stem Cell Memory-Like T Cells A User's Guide to Multicolor Flow Cytometry Panels for Comprehensive Immune Profiling Bimodal Down-Regulation of CD4 in Cells Expressing Human Immunodeficiency Virus Type 1 Vpu and Env Characterization of CC-Chemokine Receptor 7 Expression on Murine T Cells in Lymphoid Tissues Compensation of Signal Spillover in Suspension and Imaging Mass Cytometry Minimizing Batch Effects in Mass Cytometry Data An R-Based Reproducible and User-Friendly Preprocessing Pipeline for CyTOF Data The Challenge of Distinguishing Cell-Cell Complexes From Singlet Cells in Non-Imaging Flow Cytometry and Single-Cell Sorting An updated debarcoding tool for mass cytometry with cell type-specific and cell sample-specific stringency adjustment Fluorescence Spectral Overlap Compensation for Any Number of Flow Cytometry Parameters Ab Initio Spillover Compensation in Mass Cytometry Data CyTOForum • View Topic -The New Gaussian Parameters From Helios Automated Data Cleanup for Mass Cytometry Modern Multidimensional Scaling: Theory and Applications Differential Discovery in High-Throughput High-Dimensional Cytometry Datasets Arcsinhtransform: Create the Definition of an Arcsinh Transformation Function … in Flowcore: Flowcore: Basic Structures for Flow Cytometry Data Comparison of Clustering Methods for High-Dimensional Single-Cell Flow and Mass Cytometry Data Automated Mapping of Phenotype Space With Single-Cell Data A Beginner's Guide To Analyzing and Visualizing Mass Cytometry Data Data-Driven Phenotypic Dissection of AML Reveals Progenitor-Like Cells That Correlate With Prognosis Using Self-Organizing Maps for Visualization and Interpretation of Cytometry Data Analyzing High-Dimensional Cytometry Data Using FlowSOM Differential Discovery With CATALYST Deltaplotr : An R Package for Differential Item Functioning Analysis With Angoff's Delta Plot Topic -Stability of Maxpar Conjugated Antibodies Relevance of Antibody Validation for Flow Cytometry Titration of Mass Cytometry Reagents Titrating Complex Mass Cytometry Panels Selecting Reagents for Multicolor Flow Cytometry With BD ™ LSR II and BD FACSCanto ™ Systems Multiparameter Analysis of Stimulated Human Peripheral Blood Mononuclear Cells: A Comparison of Mass and Fluorescence Cytometry Another Step on the Path to Mass Cytometry Standardization Cell-ID Intercalator-103Rh-2000 µm Platelet Mass Cytometry: Optimization of Sample, Reagent, and Analysis Parameters Heparin Reduces Nonspecific Eosinophil Staining Artifacts in Mass Cytometry Experiments Isolation of Lamina Propria Mononuclear Cells From Murine Colon Using Collagenase E Single Cell Analysis of Human Tissues and Solid Tumors With Mass Cytometry Mass Cytometry Analysis of Immune Cells in the Brain Best Practices for Preparing a Single Cell Suspension From Solid Tissues for Flow Cytometry Isolation of Mouse Brain-Infiltrating Leukocytes for Single Cell Profiling of Epitopes and Transcriptomes Isolation and High-Dimensional Phenotyping of Gastrointestinal Immune Cells Methods for High-Dimensional Analysis of Cells Dissociated From Cryopreserved Synovial Tissue Using Mass Cytometry to Analyze the Tumor-Infiltrating Lymphocytes in Human Melanoma Phenotypic Characterization by Mass Cytometry of the Microenvironment in Ovarian Cancer and Impact of Tumor Dissociation Methods Enrichment and Detection of Antigen-Binding B Cells for Mass Cytometry On the Origin of Low-Density Neutrophils A Global Look Into Human T Cell Subsets Before and After Cryopreservation Using Multiparametric Flow Cytometry and Two-Dimensional Visualization Analysis The Effect of Cellular Isolation and Cryopreservation on the Expression of Markers Identifying Subsets of Regulatory T Cells Optimization and Limitations of Use of Cryopreserved Peripheral Blood Mononuclear Cells for Functional and Phenotypic T-Cell Characterization Standardization of Cryopreserved Peripheral Blood Mononuclear Cells Through a Resting Process for Clinical Immunomonitoring-Development of an Algorithm Monensin and Brefeldin A Differentially Affect the Phosphorylation and Sulfation of Secretory Proteins Evaluation of Monensin and Brefeldin A for Flow Cytometric Determination of Interleukin-1 Beta, Interleukin-6, and Tumor Necrosis Factor-Alpha in Monocytes But Not Monensin, Enables Flow Cytometric Detection of Interleukin-4 Within Peripheral T Cells Responding to Ex Vivo Stimulation With Chlamydia Trachomatis Cryopreservation-Related Loss of Antigen-Specific Ifng Producing CD4+ T-Cells can Skew Immunogenicity Data in Vaccine Trials: Lessons From a Malaria Vaccine Trial Substudy Cytokine-Stimulated Phosphoflow of Whole Blood Using CyTOF Mass Cytometry Intracellular Phospho-Protein Staining Techniques for Flow Cytometry: Monitoring Single Cell Signaling Events Antigen Masking During Fixation and Embedding, Dissected Flow Cytometric Analysis of Intracellular Phosphoproteins in Human Monocytes Analysis of Signaling Events by Dynamic Phosphoflow Cytometry Elimination of Erroneous Results in Flow Cytometry Caused by Antibody Binding to Fc Receptors on Human Monocytes and Macrophages Human Monocyte Heterogeneity as Revealed by High-Dimensional Mass Cytometry Human Blood Monocyte Subsets: A New Gating Strategy Defined Using Cell Surface Markers Identified by Mass Cytometry Crystallographic Data and Model Quality Harnessing CD16-Mediated NK Cell Functions to Enhance Therapeutic Efficacy of Tumor-Targeting Mabs Key Features of Gamma-Delta T-Cell Subsets in Human Diseases and Their Immunotherapeutic Implications OMIP-020: Phenotypic Characterization of Human gd T-Cells by Multicolor Flow Cytometry Preparation of Whole Bone Marrow for Mass Cytometry Analysis of Neutrophil-Lineage Cells A Regulatory Role for Fc Gamma Receptors (CD16 and CD32 Cisplatin Uptake in Macrophage Subtypes at the Single-Cell Level by LA-ICP-TOFMS Imaging Intracellular Cytokine Staining on PBMCs Using CyTOF ™ Mass Cytometry Multiplex MHC Class I Tetramer Combined With Intranuclear Staining by Mass Cytometry Stabilizing Antibody Cocktails for Mass Cytometry Phospho-Protein Staining Tips The Nightmare of a Single Cell: Being a Doublet Optimization of Mass Cytometry Sample Cryopreservation After Staining Mass Cytometry: Protocol for Daily Tuning and Running Cell Samples on a CyTOF Mass Cytometer CyTOF 2 User Manual Topic -Long CyTOF Runs in Helios Topic -Sudden Signal Intensity Drop During Helios Runs A Modified Injector and Sample Acquisition Protocol Can Improve Data Quality and Reduce Inter-Instrument Variability of the Helios Mass Cytometer Normalization of Mass Cytometry Data With Bead Standards EQ Four Element Calibration Beads-100 Ml . Available at Evaluating the Efficiency of Isotope Transmission for Improved Panel Design and a Comparison of the Detection Sensitivities of Mass Cytometer Instruments Comparison of CyTOF Assays Across Sites: Results of a Six-Center Pilot Study Metal-Encoded Polystyrene Microbeads as a Mass Cytometry Calibration/ Normalization Standard Covering Channels From Yttrium Data Scientist's Primer to Analysis of Mass Cytometry Data Random Function -RDocumentation Genetic and Environmental Determinants of Human NK Cell Diversity Revealed by Mass Cytometry NK Cell Receptor and Ligand Composition Influences the Clearance of SARS-CoV-2 Quantitative Comparison of Conventional and T-SNE-Guided Gating Analyses Evaluation of UMAP as an Alternative to T-SNE for Single-Cell Data Dimred and Coranking -Unifying Dimensionality Reduction in R Principal Components Analysis A Quantitative Framework for Evaluating Single-Cell Data Structure Preservation by Dimensionality Reduction Techniques Differential Discovery in High-Throughput High-Dimensional Cytometry Datasets Diffusion Maps for High-Dimensional Single-Cell Analysis of Differentiation Data Single-Cell Trajectory Detection Uncovers Progression and Regulatory Coordination in Human B Cell Development Wishbone Identifies Bifurcating Developmental Trajectories From Single-Cell Data Reversed Graph Embedding Resolves Complex Single-Cell Trajectories FCS Express Flow Cytometry Software. De Novo Software Cytobank: Providing an Analytics Platform for Community Cytometry Data Analysis and Collaboration Cytofkit: An Integrated Mass Cytometry Data Analysis Pipeline The Astrolabe Cytometry Platform, Solving Single-Cell Analysis Automated Cell Type Discovery and Classification Through Knowledge Transfer Supervised Machine Learning With CITRUS for Single Cell Biomarker Discovery Automated Identification of Stratifying Signatures in Cellular Subpopulations Photoionization and Photodissociation Methods in Mass Spectrometry Encyclopedia of Food Sciences and Nutrition Integrating Mass Cytometry Data Across Multiple Panels An Interactive Reference Framework for Modeling a Dynamic Immune System Trajectory-Based Differential Expression Analysis for Single-Cell Sequencing Data Predicting Cell Populations in Single Cell Mass Cytometry Data Semi-Automated Background Removal Limits Data Loss and Normalizes Imaging Mass Cytometry Data The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.Publisher's Note: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.Copyright © 2022 Iyer, Hamers and Pillai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.