Probability state modeling of memory CD8+ T-cell differentiation

Multi-site reproducibility of a human immunophenotyping assay in whole blood and peripheral blood mononuclear cells preparations using CyTOF technology coupled with Maxpar Pathsetter, an automated data analysis system

High-dimensional mass cytometry data potentially enable a comprehensive characterization of immune cells. In order to positively affect clinical trials and translational clinical research, this advanced technology needs to demonstrate a high reproducibility of results across multiple sites for both peripheral blood mononuclear cells (PBMC) and whole blood preparations. A dry 30-marker broad immunophenotyping panel and customized automated analysis software were recently engineered and are commercially available as the Fluidigm® Maxpar® Direct™ Immune Profiling Assay™. In this study, seven sites received whole blood and six sites received PBMC samples from single donors over a 2-week interval. Each site labeled replicate samples and acquired data on Helios™ instruments using an assay-specific acquisition template. All acquired sample files were then automatically analyzed by Maxpar Pathsetter™ software. A cleanup step eliminated debris, dead cells, aggregates, and normalization beads. The second step automatically enumerated 37 immune cell populations and performed label intensity assessments on all 30 markers. The inter-site reproducibility of the 37 quantified cell populations had consistent population frequencies, with an average %CV of 14.4% for whole blood and 17.7% for PBMC. The dry reagent coupled with automated data analysis is not only convenient but also provides a high degree of reproducibility within and among multiple test sites resulting in a comprehensive yet practical solution for deep immune phenotyping.

Quantifying the Dynamics of Hematopoiesis by In Vivo IdU Pulse-Chase, Mass Cytometry, and Mathematical Modeling

We present a new method to directly quantify the dynamics of differentiation of multiple cellular subsets in unperturbed mice. We combine a pulse-chase protocol of 5-iodo-2'-deoxyuridine (IdU) injections with subsequent analysis by mass cytometry (CyTOF) and mathematical modeling of the IdU dynamics. Measurements by CyTOF allow for a wide range of cells to be analyzed at once, due to the availability of a large staining panel without the complication of fluorescence spillover. These are also compatible with direct detection of integrated iodine signal, with minimal impact on immunophenotyping based on the surface markers. Mathematical modeling beyond a binary classification of surface marker abundance allows for a continuum of cellular states as the cells transition from one state to another. Thus, we present a complete and robust method for directly quantifying differentiation at the systemic level, allowing for system-wide comparisons between different mouse strains and/or experimental conditions. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Abnormal T Cell Frequencies, Including Cytomegalovirus-Associated Expansions, Distinguish Seroconverted Subjects at Risk for Type 1 Diabetes

We analyzed T cell subsets from cryopreserved PBMC obtained from the TrialNet Pathway to Prevention archives. We compared subjects who had previously seroconverted for one or more autoantibodies with non-seroconverted, autoantibody negative individuals. We observed a reduced frequency of MAIT cells among seroconverted subjects. Seroconverted subjects also possessed decreased frequencies of CCR4-expressing CD4 T cells, including a regulatory-like subset. Interestingly, we found an elevation of CD57+, CD28–, CD127–, CD27– CD8 T cells (SLEC) among seroconverted subjects that was most pronounced among those that progressed to disease. The frequency of these SLEC was strongly correlated with CMV IgG abundance among seroconverted subjects, associated with IA-2 levels, and most elevated among CMV+ seroconverted subjects who progressed to disease. Combined, our data indicate discrete, yet profound T cell alterations are associated with islet autoimmunity among at-risk subjects.

High-Dimensional Modeling for Cytometry: Building Rock Solid Models Using GemStone™ and Verity Cen-se'™ High-Definition t-SNE Mapping

This chapter outlines how to approach the complex tasks associated with designing models for high-dimensional cytometry data. Unlike gating approaches, modeling lends itself to automation and accounts for measurement overlap among cellular populations. Designing these models is now easier because of a new technique called high-definition t-SNE mapping. Nontrivial examples are provided that serve as a guide to create models that are consistent with data.

Impaired Subset Progression and Polyfunctionality of T Cells in Mice Exposed to Methamphetamine during Chronic LCMV Infection

Methamphetamine (METH) is a widely used psychostimulant that severely impacts the host’s innate and adaptive immune systems and has profound immunological implications. T cells play a critical role in orchestrating immune responses. We have shown recently how chronic exposure to METH affects T cell activation using a murine model of lymphocytic choriomeningitis virus (LCMV) infection. Using the TriCOM (trinary state combinations) feature of GemStone^™ to study the polyfunctionality of T cells, we have analyzed how METH affected the cytokine production pattern over the course of chronic LCMV infection. Furthermore, we have studied in detail the effects of METH on splenic T cell functions, such as cytokine production and degranulation, and how they regulate each other. We used the Probability State Modeling (PSM) program to visualize the differentiation of effector/memory T cell subsets during LCMV infection and analyze the effects of METH on T cell subset progression. We recently demonstrated that METH increased PD-1 expression on T cells during viral infection. In this study, we further analyzed the impact of PD-1 expression on T cell functional markers as well as its expression in the effector/memory subsets. Overall, our study indicates that analyzing polyfunctionality of T cells can provide additional insight into T cell effector functions. Analysis of T cell heterogeneity is important to highlight changes in the evolution of memory/effector functions during chronic viral infections. Our study also highlights the impact of METH on PD-1 expression and its consequences on T cell responses.

Computational flow cytometry: helping to make sense of high-dimensional immunology data

Recent advances in flow cytometry allow scientists to measure an increasing number of parameters per cell, generating huge and high-dimensional datasets. To analyse, visualize and interpret these data, newly available computational techniques should be adopted, evaluated and improved upon by the immunological community. Computational flow cytometry is emerging as an important new field at the intersection of immunology and computational biology; it allows new biological knowledge to be extracted from high-throughput single-cell data. This Review provides non-experts with a broad and practical overview of the many recent developments in computational flow cytometry.

Stochastic Measurement Models for Quantifying Lymphocyte Responses Using Flow Cytometry

Adaptive immune responses are complex dynamic processes whereby B and T cells undergo division and differentiation triggered by pathogenic stimuli. Deregulation of the response can lead to severe consequences for the host organism ranging from immune deficiencies to autoimmunity. Tracking cell division and differentiation by flow cytometry using fluorescent probes is a major method for measuring progression of lymphocyte responses, both in vitro and in vivo. In turn, mathematical modeling of cell numbers derived from such measurements has led to significant biological discoveries, and plays an increasingly important role in lymphocyte research. Fitting an appropriate parameterized model to such data is the goal of these studies but significant challenges are presented by the variability in measurements. This variation results from the sum of experimental noise and intrinsic probabilistic differences in cells and is difficult to characterize analytically. Current model fitting methods adopt different simplifying assumptions to describe the distribution of such measurements and these assumptions have not been tested directly. To help inform the choice and application of appropriate methods of model fitting to such data we studied the errors associated with flow cytometry measurements from a wide variety of experiments. We found that the mean and variance of the noise were related by a power law with an exponent between 1.3 and 1.8 for different datasets. This violated the assumptions inherent to commonly used least squares, linear variance scaling and log-transformation based methods. As a result of these findings we propose a new measurement model that we justify both theoretically, from the maximum entropy standpoint, and empirically using collected data. Our evaluation suggests that the new model can be reliably used for model fitting across a variety of conditions. Our work provides a foundation for modeling measurements in flow cytometry experiments thus facilitating progress in quantitative studies of lymphocyte responses.

Algorithmic Tools for Mining High-Dimensional Cytometry Data

The advent of mass cytometry has led to an unprecedented increase in the number of analytes measured in individual cells, thereby increasing the complexity and information content of cytometric data. Although this technology is ideally suited to the detailed examination of the immune system, the applicability of the different methods for analyzing such complex data is less clear. Conventional data analysis by manual gating of cells in biaxial dot plots is often subjective, time consuming, and neglectful of much of the information contained in a highly dimensional cytometric dataset. Algorithmic data mining has the promise to eliminate these concerns, and several such tools have been applied recently to mass cytometry data. We review computational data mining tools that have been used to analyze mass cytometry data, outline their differences, and comment on their strengths and limitations. This review will help immunologists to identify suitable algorithmic tools for their particular projects.

Human B-cell and progenitor stages as determined by probability state modeling of multidimensional cytometry data

BACKGROUND

Human progenitor and B-cell development is a highly regulated process characterized by the ordered differential expression of numerous cell-surface and intracytoplasmic antigens. This study investigates the underlying coordination of these modulations by examining a series of normal bone marrow samples with the method of probability state modeling or PSM.

RESULTS

The study is divided into two sections. The first section examines B-cell stages subsequent to CD19 up-regulation. The second section assesses an earlier differentiation stage before and including CD19 up-regulation. POST-CD19 ANTIGENIC UP-REGULATION: Statistical analyses of cytometry data derived from sixteen normal bone marrow specimens revealed that B cells have at least three distinct coordinated changes, forming four stages labeled as B1, B2, B3, and B4. At the end of B1; CD34 antigen expression down-regulates with TdT while CD45, CD81, and CD20 slightly up-regulate. At the end of B2, CD45 and CD20 up-regulate. At the end of B3 and beginning of B4; CD10, CD38, and CD81 down-regulate while CD22 and CD44 up-regulate. PRE-CD19 ANTIGENIC UP-REGULATION: Statistical analysis of ten normal bone marrows revealed that there are at least two measurable coordinated changes with progenitors, forming three stages labeled as P1, P2, and P3. At the end of P1, CD38 up-regulates. At the end of P2; CD19, CD10, CD81, CD22, and CD9 up-regulate while CD44 down-regulates slightly.

CONCLUSIONS

These objective results yield a clearer immunophenotypic picture of the underlying cellular mechanisms that are operating in these important developmental processes. Also, unambiguously determined stages define what is meant by "normal" B-cell development and may serve as a preliminary step for the development of highly sensitive minimum residual disease detection systems. A companion article is simultaneously being published in Cytometry Part A that will explain in further detail the theory behind PSM. Three short relevant videos are available in the online supporting information for both of these papers.

Methods for discovery and characterization of cell subsets in high dimensional mass cytometry data

The flood of high-dimensional data resulting from mass cytometry experiments that measure more than 40 features of individual cells has stimulated creation of new single cell computational biology tools. These tools draw on advances in the field of machine learning to capture multi-parametric relationships and reveal cells that are easily overlooked in traditional analysis. Here, we introduce a workflow for high dimensional mass cytometry data that emphasizes unsupervised approaches and visualizes data in both single cell and population level views. This workflow includes three central components that are common across mass cytometry analysis approaches: (1) distinguishing initial populations, (2) revealing cell subsets, and (3) characterizing subset features. In the implementation described here, viSNE, SPADE, and heatmaps were used sequentially to comprehensively characterize and compare healthy and malignant human tissue samples. The use of multiple methods helps provide a comprehensive view of results, and the largely unsupervised workflow facilitates automation and helps researchers avoid missing cell populations with unusual or unexpected phenotypes. Together, these methods develop a framework for future machine learning of cell identity.

Automated analysis of flow cytometric data for measuring neutrophil CD64 expression using a multi-instrument compatible probability state model

BACKGROUND

Leuko64™ (Trillium Diagnostics) is a flow cytometric assay that measures neutrophil CD64 expression and serves as an in vitro indicator of infection/sepsis or the presence of a systemic acute inflammatory response. Leuko64 assay currently utilizes QuantiCALC, a semiautomated software that employs cluster algorithms to define cell populations. The software reduces subjective gating decisions, resulting in interanalyst variability of <5%. We evaluated a completely automated approach to measuring neutrophil CD64 expression using GemStone™ (Verity Software House) and probability state modeling (PSM).

METHODS

Four hundred and fifty-seven human blood samples were processed using the Leuko64 assay. Samples were analyzed on four different flow cytometer models: BD FACSCanto II, BD FACScan, BC Gallios/Navios, and BC FC500. A probability state model was designed to identify calibration beads and three leukocyte subpopulations based on differences in intensity levels of several parameters. PSM automatically calculates CD64 index values for each cell population using equations programmed into the model. GemStone software uses PSM that requires no operator intervention, thus totally automating data analysis and internal quality control flagging. Expert analysis with the predicate method (QuantiCALC) was performed. Interanalyst precision was evaluated for both methods of data analysis.

RESULTS

PSM with GemStone correlates well with the expert manual analysis, r(2) = 0.99675 for the neutrophil CD64 index values with no intermethod bias detected. The average interanalyst imprecision for the QuantiCALC method was 1.06% (range 0.00-7.94%), which was reduced to 0.00% with the GemStone PSM. The operator-to-operator agreement in GemStone was a perfect correlation, r(2) = 1.000.

CONCLUSION

Automated quantification of CD64 index values produced results that strongly correlate with expert analysis using a standard gate-based data analysis method. PSM successfully evaluated flow cytometric data generated by multiple instruments across multiple lots of the Leuko64 kit in all 457 cases. The probability-based method provides greater objectivity, higher data analysis speed, and allows for greater precision for in vitro diagnostic flow cytometric assays.