Identification of cell types in multiplexed in situ images by combining protein expression and spatial information using CELESTA.

Advances in multiplexed in situ imaging are revealing important insights in spatial biology. However, cell type identification remains a major challenge in imaging analysis, with most existing methods involving substantial manual assessment and subjective decisions for thousands of cells. We developed an unsupervised machine learning algorithm, CELESTA, which identifies the cell type of each cell, individually, using the cell's marker expression profile and, when needed, its spatial information. We demonstrate the performance of CELESTA on multiplexed immunofluorescence images of colorectal cancer and head and neck squamous cell carcinoma (HNSCC). Using the cell types identified by CELESTA, we identify tissue architecture associated with lymph node metastasis in HNSCC, and validate our findings in an independent cohort. By coupling our spatial analysis with single-cell RNA-sequencing data on proximal sections of the same specimens, we identify cell-cell crosstalk associated with lymph node metastasis, demonstrating the power of CELESTA to facilitate identification of clinically relevant interactions.

CODEX multiplexed tissue imaging with DNA-conjugated antibodies.

Advances in multiplexed imaging technologies have drastically improved our ability to characterize healthy and diseased tissues at the single-cell level. Co-detection by indexing (CODEX) relies on DNA-conjugated antibodies and the cyclic addition and removal of complementary fluorescently labeled DNA probes and has been used so far to simultaneously visualize up to 60 markers in situ. CODEX enables a deep view into the single-cell spatial relationships in tissues and is intended to spur discovery in developmental biology, disease and therapeutic design. Herein, we provide optimized protocols for conjugating purified antibodies to DNA oligonucleotides, validating the conjugation by CODEX staining and executing the CODEX multicycle imaging procedure for both formalin-fixed, paraffin-embedded (FFPE) and fresh-frozen tissues. In addition, we describe basic image processing and data analysis procedures. We apply this approach to an FFPE human tonsil multicycle experiment. The hands-on experimental time for antibody conjugation is ~4.5 h, validation of DNA-conjugated antibodies with CODEX staining takes ~6.5 h and preparation for a CODEX multicycle experiment takes ~8 h. The multicycle imaging and data analysis time depends on the tissue size, number of markers in the panel and computational complexity.

Highly multiplexed tissue imaging using repeated oligonucleotide exchange reaction.

Multiparameter tissue imaging enables analysis of cell-cell interactions in situ, the cellular basis for tissue structure, and novel cell types that are spatially restricted, giving clues to biological mechanisms behind tissue homeostasis and disease. Here, we streamlined and simplified the multiplexed imaging method CO-Detection by indEXing (CODEX) by validating 58 unique oligonucleotide barcodes that can be conjugated to antibodies. We showed that barcoded antibodies retained their specificity for staining cognate targets in human tissue. Antibodies were visualized one at a time by adding a fluorescently labeled oligonucleotide complementary to oligonucleotide barcode, imaging, stripping, and repeating this cycle. With this we developed a panel of 46 antibodies that was used to stain five human lymphoid tissues: three tonsils, a spleen, and a LN. To analyze the data produced, an image processing and analysis pipeline was developed that enabled single-cell analysis on the data, including unsupervised clustering, that revealed 31 cell types across all tissues. We compared cell-type compositions within and directly surrounding follicles from the different lymphoid organs and evaluated cell-cell density correlations. This sequential oligonucleotide exchange technique enables a facile imaging of tissues that leverages pre-existing imaging infrastructure to decrease the barriers to broad use of multiplexed imaging.

Coordinated Cellular Neighborhoods Orchestrate Antitumoral Immunity at the Colorectal Cancer Invasive Front.

Antitumoral immunity requires organized, spatially nuanced interactions between components of the immune tumor microenvironment (iTME). Understanding this coordinated behavior in effective versus ineffective tumor control will advance immunotherapies. We re-engineered co-detection by indexing (CODEX) for paraffin-embedded tissue microarrays, enabling simultaneous profiling of 140 tissue regions from 35 advanced-stage colorectal cancer (CRC) patients with 56 protein markers. We identified nine conserved, distinct cellular neighborhoods (CNs)-a collection of components characteristic of the CRC iTME. Enrichment of PD-1+CD4+ T cells only within a granulocyte CN positively correlated with survival in a high-risk patient subset. Coupling of tumor and immune CNs, fragmentation of T cell and macrophage CNs, and disruption of inter-CN communication was associated with inferior outcomes. This study provides a framework for interrogating how complex biological processes, such as antitumoral immunity, occur through concerted actions of cells and spatial domains.

Author Correction: Ultra-high throughput single-cell analysis of proteins and RNAs by split-pool synthesis.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Ultra-high throughput single-cell analysis of proteins and RNAs by split-pool synthesis.

Single-cell omics provide insight into cellular heterogeneity and function. Recent technological advances have accelerated single-cell analyses, but workflows remain expensive and complex. We present a method enabling simultaneous, ultra-high throughput single-cell barcoding of millions of cells for targeted analysis of proteins and RNAs. Quantum barcoding (QBC) avoids isolation of single cells by building cell-specific oligo barcodes dynamically within each cell. With minimal instrumentation (four 96-well plates and a multichannel pipette), cell-specific codes are added to each tagged molecule within cells through sequential rounds of classical split-pool synthesis. Here we show the utility of this technology in mouse and human model systems for as many as 50 antibodies to targeted proteins and, separately, >70 targeted RNA regions. We demonstrate that this method can be applied to multi-modal protein and RNA analyses. It can be scaled by expansion of the split-pool process and effectively renders sequencing instruments as versatile multi-parameter flow cytometers.

Profiling myelodysplastic syndromes by mass cytometry demonstrates abnormal progenitor cell phenotype and differentiation.

BACKGROUND: We sought to enhance the cytometric analysis of myelodysplastic syndromes (MDS) by performing a pilot study of a single cell mass cytometry (MCM) assay to more comprehensively analyze patterns of surface marker expression in patients with MDS. METHODS: Twenty-three MDS and five healthy donor bone marrow samples were studied using a 34-parameter mass cytometry panel utilizing barcoding and internal reference standards. The resulting data were analyzed by both traditional gating and high-dimensional clustering. RESULTS: This high-dimensional assay provided three major benefits relative to traditional cytometry approaches: First, MCM enabled detection of aberrant surface maker at high resolution, detecting aberrancies in 27/31 surface markers, encompassing almost every previously reported MDS surface marker aberrancy. Additionally, three previously unrecognized aberrancies in MDS were detected in multiple samples at least one developmental stage: increased CD321 and CD99; and decreased CD47. Second, analysis of the stem and progenitor cell compartment (HSPCs), demonstrated aberrant expression in 21 of the 23 MDS samples, which were not detected in three samples from patients with idiopathic cytopenia of undetermined significance. These immunophenotypically abnormal HSPCs were also the single most significant distinguishing feature between clinical risk groups. Third, unsupervised clustering of high-parameter MCM data allowed identification of abnormal differentiation patterns associated with immunophenotypically aberrant myeloid cells similar to myeloid derived suppressor cells. CONCLUSIONS: These results demonstrate that high-parameter cytometry methods that enable simultaneous analysis of all bone marrow cell types could enhance the diagnostic utility of immunophenotypic analysis in MDS.

Proliferation tracing with single-cell mass cytometry optimizes generation of stem cell memory-like T cells.

Selective differentiation of naive T cells into multipotent T cells is of great interest clinically for the generation of cell-based cancer immunotherapies. Cellular differentiation depends crucially on division state and time. Here we adapt a dye dilution assay for tracking cell proliferative history through mass cytometry and uncouple division, time and regulatory protein expression in single naive human T cells during their activation and expansion in a complex ex vivo milieu. Using 23 markers, we defined groups of proteins controlled predominantly by division state or time and found that undivided cells account for the majority of phenotypic diversity. We next built a map of cell state changes during naive T-cell expansion. By examining cell signaling on this map, we rationally selected ibrutinib, a BTK and ITK inhibitor, and administered it before T cell activation to direct differentiation toward a T stem cell memory (TSCM)-like phenotype. This method for tracing cell fate across division states and time can be broadly applied for directing cellular differentiation.

Deep Profiling of Mouse Splenic Architecture with CODEX Multiplexed Imaging.

A highly multiplexed cytometric imaging approach, termed co-detection by indexing (CODEX), is used here to create multiplexed datasets of normal and lupus (MRL/lpr) murine spleens. CODEX iteratively visualizes antibody binding events using DNA barcodes, fluorescent dNTP analogs, and an in situ polymerization-based indexing procedure. An algorithmic pipeline for single-cell antigen quantification in tightly packed tissues was developed and used to overlay well-known morphological features with de novo characterization of lymphoid tissue architecture at a single-cell and cellular neighborhood levels. We observed an unexpected, profound impact of the cellular neighborhood on the expression of protein receptors on immune cells. By comparing normal murine spleen to spleens from animals with systemic autoimmune disease (MRL/lpr), extensive and previously uncharacterized splenic cell-interaction dynamics in the healthy versus diseased state was observed. The fidelity of multiplexed spatial cytometry demonstrated here allows for quantitative systemic characterization of tissue architecture in normal and clinically aberrant samples.

Three-dimensional intact-tissue sequencing of single-cell transcriptional states.

Retrieving high-content gene-expression information while retaining three-dimensional (3D) positional anatomy at cellular resolution has been difficult, limiting integrative understanding of structure and function in complex biological tissues. We developed and applied a technology for 3D intact-tissue RNA sequencing, termed STARmap (spatially-resolved transcript amplicon readout mapping), which integrates hydrogel-tissue chemistry, targeted signal amplification, and in situ sequencing. The capabilities of STARmap were tested by mapping 160 to 1020 genes simultaneously in sections of mouse brain at single-cell resolution with high efficiency, accuracy, and reproducibility. Moving to thick tissue blocks, we observed a molecularly defined gradient distribution of excitatory-neuron subtypes across cubic millimeter-scale volumes (>30,000 cells) and a short-range 3D self-clustering in many inhibitory-neuron subtypes that could be identified and described with 3D STARmap.

GateFinder: projection-based gating strategy optimization for flow and mass cytometry.

Motivation: High-parameter single-cell technologies can reveal novel cell populations of interest, but studying or validating these populations using lower-parameter methods remains challenging. Results: Here, we present GateFinder, an algorithm that enriches high-dimensional cell types with simple, stepwise polygon gates requiring only two markers at a time. A series of case studies of complex cell types illustrates how simplified enrichment strategies can enable more efficient assays, reveal novel biomarkers and clarify underlying biology. Availability and implementation: The GateFinder algorithm is implemented as a free and open-source package for BioConductor: https://nalab.stanford.edu/gatefinder. Supplementary information: Supplementary data are available at Bioinformatics online.

Publisher Correction: High-resolution myogenic lineage mapping by single-cell mass cytometry.

In the version of this Article originally published, the name of author Andrew Tri Van Ho was coded wrongly, resulting in it being incorrect when exported to citation databases. This has been corrected, though no visible changes will be apparent.

Single-cell mass cytometry reveals distinct populations of brain myeloid cells in mouse neuroinflammation and neurodegeneration models.

Neuroinflammation and neurodegeneration may represent two poles of brain pathology. Brain myeloid cells, particularly microglia, play key roles in these conditions. We employed single-cell mass cytometry (CyTOF) to compare myeloid cell populations in the experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis, the R6/2 model of Huntington's disease (HD) and the mutant superoxide dismutase 1 (mSOD1) model of amyotrophic lateral sclerosis (ALS). We identified three myeloid cell populations exclusive to the CNS and present in each disease model. Blood-derived monocytes comprised five populations and migrated to the brain in EAE, but not in HD and ALS models. Single-cell analysis resolved differences in signaling and cytokine production within similar myeloid populations in EAE compared to HD and ALS models. Moreover, these analyses highlighted α5 integrin on myeloid cells as a potential therapeutic target for neuroinflammation. Together, these findings illustrate how neuropathology may differ between inflammatory and degenerative brain disease.

Single-cell developmental classification of B cell precursor acute lymphoblastic leukemia at diagnosis reveals predictors of relapse.

Insight into the cancer cell populations that are responsible for relapsed disease is needed to improve outcomes. Here we report a single-cell-based study of B cell precursor acute lymphoblastic leukemia at diagnosis that reveals hidden developmentally dependent cell signaling states that are uniquely associated with relapse. By using mass cytometry we simultaneously quantified 35 proteins involved in B cell development in 60 primary diagnostic samples. Each leukemia cell was then matched to its nearest healthy B cell population by a developmental classifier that operated at the single-cell level. Machine learning identified six features of expanded leukemic populations that were sufficient to predict patient relapse at diagnosis. These features implicated the pro-BII subpopulation of B cells with activated mTOR signaling, and the pre-BI subpopulation of B cells with activated and unresponsive pre-B cell receptor signaling, to be associated with relapse. This model, termed 'developmentally dependent predictor of relapse' (DDPR), significantly improves currently established risk stratification methods. DDPR features exist at diagnosis and persist at relapse. By leveraging a data-driven approach, we demonstrate the predictive value of single-cell 'omics' for patient stratification in a translational setting and provide a framework for its application to human cancer.

Commonly Occurring Cell Subsets in High-Grade Serous Ovarian Tumors Identified by Single-Cell Mass Cytometry.

We have performed an in-depth single-cell phenotypic characterization of high-grade serous ovarian cancer (HGSOC) by multiparametric mass cytometry (CyTOF). Using a CyTOF antibody panel to interrogate features of HGSOC biology, combined with unsupervised computational analysis, we identified noteworthy cell types co-occurring across the tumors. In addition to a dominant cell subset, each tumor harbored rarer cell phenotypes. One such group co-expressed E-cadherin and vimentin (EV), suggesting their potential role in epithelial mesenchymal transition, which was substantiated by pairwise correlation analyses. Furthermore, tumors from patients with poorer outcome had an increased frequency of another rare cell type that co-expressed vimentin, HE4, and cMyc. These poorer-outcome tumors also populated more cell phenotypes, as quantified by Simpson's diversity index. Thus, despite the recognized genomic complexity of the disease, the specific cell phenotypes uncovered here offer a focus for therapeutic intervention and disease monitoring.

Scalable multi-sample single-cell data analysis by Partition-Assisted Clustering and Multiple Alignments of Networks.

Mass cytometry (CyTOF) has greatly expanded the capability of cytometry. It is now easy to generate multiple CyTOF samples in a single study, with each sample containing single-cell measurement on 50 markers for more than hundreds of thousands of cells. Current methods do not adequately address the issues concerning combining multiple samples for subpopulation discovery, and these issues can be quickly and dramatically amplified with increasing number of samples. To overcome this limitation, we developed Partition-Assisted Clustering and Multiple Alignments of Networks (PAC-MAN) for the fast automatic identification of cell populations in CyTOF data closely matching that of expert manual-discovery, and for alignments between subpopulations across samples to define dataset-level cellular states. PAC-MAN is computationally efficient, allowing the management of very large CyTOF datasets, which are increasingly common in clinical studies and cancer studies that monitor various tissue samples for each subject.

High-resolution myogenic lineage mapping by single-cell mass cytometry.

Muscle regeneration is a dynamic process during which cell state and identity change over time. A major roadblock has been a lack of tools to resolve a myogenic progression in vivo. Here we capitalize on a transformative technology, single-cell mass cytometry (CyTOF), to identify in vivo skeletal muscle stem cell and previously unrecognized progenitor populations that precede differentiation. We discovered two cell surface markers, CD9 and CD104, whose combined expression enabled in vivo identification and prospective isolation of stem and progenitor cells. Data analysis using the X-shift algorithm paired with single-cell force-directed layout visualization defined a molecular signature of the activated stem cell state (CD44+/CD98+/MyoD+) and delineated a myogenic trajectory during recovery from acute muscle injury. Our studies uncover the dynamics of skeletal muscle regeneration in vivo and pave the way for the elucidation of the regulatory networks that underlie cell-state transitions in muscle diseases and ageing.

SESSION INTRODUCTION.

Recent technological developments allow gathering single-cell measurements across different domains (genomic, transcriptomics, proteomics, imaging etc). Sophisticated computational algorithms are required in order to harness the power of single-cell data. This session is dedicated to computational methods for single-cell analysis in various biological domains, modelling of population heterogeneity, as well as translational applications of single cell data.

Diabetes-linked transcription factor HNF4α regulates metabolism of endogenous methylarginines and ß-aminoisobutyric acid by controlling expression of alanine-glyoxylate aminotransferase 2.

Elevated levels of circulating asymmetric and symmetric dimethylarginines (ADMA and SDMA) predict and potentially contribute to end organ damage in cardiovascular diseases. Alanine-glyoxylate aminotransferase 2 (AGXT2) regulates systemic levels of ADMA and SDMA, and also of beta-aminoisobutyric acid (BAIB)-a modulator of lipid metabolism. We identified a putative binding site for hepatic nuclear factor 4 α (HNF4α) in AGXT2 promoter sequence. In a luciferase reporter assay we found a 75% decrease in activity of Agxt2 core promoter after disruption of the HNF4α binding site. Direct binding of HNF4α to Agxt2 promoter was confirmed by chromatin immunoprecipitation assay. siRNA-mediated knockdown of Hnf4a led to an almost 50% reduction in Agxt2 mRNA levels in Hepa 1-6 cells. Liver-specific Hnf4a knockout mice exhibited a 90% decrease in liver Agxt2 expression and activity, and elevated plasma levels of ADMA, SDMA and BAIB, compared to wild-type littermates. Thus we identified HNF4α as a major regulator of Agxt2 expression. Considering a strong association between human HNF4A polymorphisms and increased risk of type 2 diabetes our current findings suggest that downregulation of AGXT2 and subsequent impairment in metabolism of dimethylarginines and BAIB caused by HNF4α deficiency might contribute to development of cardiovascular complications in diabetic patients.