An end-to-end workflow for multiplexed image processing and analysis.

Multiplexed imaging enables the simultaneous spatial profiling of dozens of biological molecules in tissues at single-cell resolution. Extracting biologically relevant information, such as the spatial distribution of cell phenotypes from multiplexed tissue imaging data, involves a number of computational tasks, including image segmentation, feature extraction and spatially resolved single-cell analysis. Here, we present an end-to-end workflow for multiplexed tissue image processing and analysis that integrates previously developed computational tools to enable these tasks in a user-friendly and customizable fashion. For data quality assessment, we highlight the utility of napari-imc for interactively inspecting raw imaging data and the cytomapper R/Bioconductor package for image visualization in R. Raw data preprocessing, image segmentation and feature extraction are performed using the steinbock toolkit. We showcase two alternative approaches for segmenting cells on the basis of supervised pixel classification and pretrained deep learning models. The extracted single-cell data are then read, processed and analyzed in R. The protocol describes the use of community-established data containers, facilitating the application of R/Bioconductor packages for dimensionality reduction, single-cell visualization and phenotyping. We provide instructions for performing spatially resolved single-cell analysis, including community analysis, cellular neighborhood detection and cell-cell interaction testing using the imcRtools R/Bioconductor package. The workflow has been previously applied to imaging mass cytometry data, but can be easily adapted to other highly multiplexed imaging technologies. This protocol can be implemented by researchers with basic bioinformatics training, and the analysis of the provided dataset can be completed within 5-6 h. An extended version is available at https://bodenmillergroup.github.io/IMCDataAnalysis/ .

Breast tumor microenvironment structures are associated with genomic features and clinical outcome.

The functions of the tumor microenvironment (TME) are orchestrated by precise spatial organization of specialized cells, yet little is known about the multicellular structures that form within the TME. Here we systematically mapped TME structures in situ using imaging mass cytometry and multitiered spatial analysis of 693 breast tumors linked to genomic and clinical data. We identified ten recurrent TME structures that varied by vascular content, stromal quiescence versus activation, and leukocyte composition. These TME structures had distinct enrichment patterns among breast cancer subtypes, and some were associated with genomic profiles indicative of immune escape. Regulatory and dysfunctional T cells co-occurred in large 'suppressed expansion' structures. These structures were characterized by high cellular diversity, proliferating cells and enrichment for BRCA1 and CASP8 mutations and predicted poor outcome in estrogen-receptor-positive disease. The multicellular structures revealed here link conserved spatial organization to local TME function and could improve patient stratification.

A quantitative analysis of the interplay of environment, neighborhood, and cell state in 3D spheroids.

Cells react to their microenvironment by integrating external stimuli into phenotypic decisions via an intracellular signaling network. To analyze the interplay of environment, local neighborhood, and internal cell state effects on phenotypic variability, we developed an experimental approach that enables multiplexed mass cytometric imaging analysis of up to 240 pooled spheroid microtissues. We quantified the contributions of environment, neighborhood, and intracellular state to marker variability in single cells of the spheroids. A linear model explained on average more than half of the variability of 34 markers across four cell lines and six growth conditions. The contributions of cell-intrinsic and environmental factors to marker variability are hierarchically interdependent, a finding that we propose has general implications for systems-level studies of single-cell phenotypic variability. By the overexpression of 51 signaling protein constructs in subsets of cells, we also identified proteins that have cell-intrinsic and cell-extrinsic effects. Our study deconvolves factors influencing cellular phenotype in a 3D tissue and provides a scalable experimental system, analytical principles, and rich multiplexed imaging datasets for future studies.

The single-cell pathology landscape of breast cancer.

Single-cell analyses have revealed extensive heterogeneity between and within human tumours1-4, but complex single-cell phenotypes and their spatial context are not at present reflected in the histological stratification that is the foundation of many clinical decisions. Here we use imaging mass cytometry5 to simultaneously quantify 35 biomarkers, resulting in 720 high-dimensional pathology images of tumour tissue from 352 patients with breast cancer, with long-term survival data available for 281 patients. Spatially resolved, single-cell analysis identified the phenotypes of tumour and stromal single cells, their organization and their heterogeneity, and enabled the cellular architecture of breast cancer tissue to be characterized on the basis of cellular composition and tissue organization. Our analysis reveals multicellular features of the tumour microenvironment and novel subgroups of breast cancer that are associated with distinct clinical outcomes. Thus, spatially resolved, single-cell analysis can characterize intratumour phenotypic heterogeneity in a disease-relevant manner, with the potential to inform patient-specific diagnosis.

Imaging mass cytometry and multiplatform genomics define the phenogenomic landscape of breast cancer.

Genomic alterations shape cell phenotypes and the structure of tumor ecosystems in poorly defined ways. To investigate these relationships, we used imaging mass cytometry to quantify the expression of 37 proteins with subcellular spatial resolution in 483 tumors from the METABRIC cohort. Single-cell analysis revealed cell phenotypes spanning epithelial, stromal and immune types. Distinct combinations of cell phenotypes and cell-cell interactions were associated with genomic subtypes of breast cancer. Epithelial luminal cell phenotypes separated into those predominantly impacted by mutations and those affected by copy number aberrations. Several features of tumor ecosystems, including cellular neighborhoods, were linked to prognosis, illustrating their clinical relevance. In summary, systematic analysis of single-cell phenotypic and spatial correlates of genomic alterations in cancer revealed how genomes shape both the composition and architecture of breast tumor ecosystems and will enable greater understanding of the phenotypic impact of genomic alterations.

Analysis of the Human Kinome and Phosphatome by Mass Cytometry Reveals Overexpression-Induced Effects on Cancer-Related Signaling.

Kinase and phosphatase overexpression drives tumorigenesis and drug resistance. We previously developed a mass-cytometry-based single-cell proteomics approach that enables quantitative assessment of overexpression effects on cell signaling. Here, we applied this approach in a human kinome- and phosphatome-wide study to assess how 649 individually overexpressed proteins modulated cancer-related signaling in HEK293T cells in an abundance-dependent manner. Based on these data, we expanded the functional classification of human kinases and phosphatases and showed that the overexpression effects include non-catalytic roles. We detected 208 previously unreported signaling relationships. The signaling dynamics analysis indicated that the overexpression of ERK-specific phosphatases sustains proliferative signaling. This suggests a phosphatase-driven mechanism of cancer progression. Moreover, our analysis revealed a drug-resistant mechanism through which overexpression of tyrosine kinases, including SRC, FES, YES1, and BLK, induced MEK-independent ERK activation in melanoma A375 cells. These proteins could predict drug sensitivity to BRAF-MEK concurrent inhibition in cells carrying BRAF mutations.

In-Depth Characterization of Monocyte-Derived Macrophages using a Mass Cytometry-Based Phagocytosis Assay.

Phagocytosis is a process in which target cells or particles are engulfed and taken up by other cells, typically professional phagocytes; this process is crucial in many physiological processes and disease states. The detection of targets for phagocytosis is directed by a complex repertoire of cell surface receptors. Pattern recognition receptors directly detect targets for binding and uptake, while opsonic and complement receptors detect objects coated by soluble factors. However, the importance of single and combinatorial surface marker expression across different phenotypes of professional phagocytes is not known. Here we developed a novel mass cytometry-based phagocytosis assay that enables the simultaneous detection of phagocytic events in combination with up to 40 other protein markers. We applied this assay to distinct monocyte derived macrophage (MDM) populations and found that prototypic M2-like MDMs phagocytose more E. coli than M1-like MDMs. Surface markers such as CD14, CD206, and CD163 rendered macrophages phagocytosis competent, but only CD209 directly correlated with the amount of particle uptake. Similarly, M2-like MDMs also phagocytosed more cancer cells than M1-like MDMs but, unlike M1-like MDMs, were insensitive to anti-CD47 opsonization. Our approach facilitates the simultaneous study of single-cell phenotypes, phagocytic activity, signaling and transcriptional events in complex cell mixtures.

A Map of Human Type 1 Diabetes Progression by Imaging Mass Cytometry.

Type 1 diabetes (T1D) results from the autoimmune destruction of insulin-producing ß cells. A comprehensive picture of the changes during T1D development is lacking due to limited sample availability, inability to sample longitudinally, and the paucity of technologies enabling comprehensive tissue profiling. Here, we analyzed 1,581 islets from 12 human donors, including eight with T1D, using imaging mass cytometry (IMC). IMC enabled simultaneous measurement of 35 biomarkers with single-cell and spatial resolution. We performed pseudotime analysis of islets through T1D progression from snapshot data to reconstruct the evolution of ß cell loss and insulitis. Our analyses revealed that ß cell destruction is preceded by a ß cell marker loss and by recruitment of cytotoxic and helper T cells. The approaches described herein demonstrate the value of IMC for improving our understanding of T1D pathogenesis, and our data lay the foundation for hypothesis generation and follow-on experiments.

Compensation of Signal Spillover in Suspension and Imaging Mass Cytometry.

The advent of mass cytometry increased the number of parameters measured at the single-cell level while decreasing the extent of crosstalk between channels relative to dye-based flow cytometry. Although reduced, spillover still exists in mass cytometry data, and minimizing its effect requires considerable expert knowledge and substantial experimental effort. Here, we describe a novel bead-based compensation workflow and R-based software that estimates and corrects for interference between channels. We performed an in-depth characterization of the spillover properties in mass cytometry, including limitations defined by the linear range of the mass cytometer and the reproducibility of the spillover over time and across machines. We demonstrated the utility of our method in suspension and imaging mass cytometry. To conclude, our approach greatly simplifies the development of new antibody panels, increases flexibility for antibody-metal pairing, opens the way to using less pure isotopes, and improves overall data quality, thereby reducing the risk of reporting cell phenotype artifacts.

Simultaneous Multiplexed Imaging of mRNA and Proteins with Subcellular Resolution in Breast Cancer Tissue Samples by Mass Cytometry.

To build comprehensive models of cellular states and interactions in normal and diseased tissue, genetic and proteomic information must be extracted with single-cell and spatial resolution. Here, we extended imaging mass cytometry to enable multiplexed detection of mRNA and proteins in tissues. Three mRNA target species were detected by RNAscope-based metal in situ hybridization with simultaneous antibody detection of 16 proteins. Analysis of 70 breast cancer samples showed that HER2 and CK19 mRNA and protein levels are moderately correlated on the single-cell level, but that only HER2, and not CK19, has strong mRNA-to-protein correlation on the cell population level. The chemoattractant CXCL10 was expressed in stromal cell clusters, and the frequency of CXCL10-expressing cells correlated with T cell presence. Our flexible and expandable method will allow an increase in the information content retrieved from patient samples for biomedical purposes, enable detailed studies of tumor biology, and serve as a tool to bridge comprehensive genomic and proteomic tissue analysis.

histoCAT: analysis of cell phenotypes and interactions in multiplex image cytometry data.

Single-cell, spatially resolved omics analysis of tissues is poised to transform biomedical research and clinical practice. We have developed an open-source, computational histology topography cytometry analysis toolbox (histoCAT) to enable interactive, quantitative, and comprehensive exploration of individual cell phenotypes, cell-cell interactions, microenvironments, and morphological structures within intact tissues. We highlight the unique abilities of histoCAT through analysis of highly multiplexed mass cytometry images of human breast cancer tissues.

An Immune Atlas of Clear Cell Renal Cell Carcinoma.

Immune cells in the tumor microenvironment modulate cancer progression and are attractive therapeutic targets. Macrophages and T cells are key components of the microenvironment, yet their phenotypes and relationships in this ecosystem and to clinical outcomes are ill defined. We used mass cytometry with extensive antibody panels to perform in-depth immune profiling of samples from 73 clear cell renal cell carcinoma (ccRCC) patients and five healthy controls. In 3.5 million measured cells, we identified 17 tumor-associated macrophage phenotypes, 22 T cell phenotypes, and a distinct immune composition correlated with progression-free survival, thereby presenting an in-depth human atlas of the immune tumor microenvironment in this disease. This study revealed potential biomarkers and targets for immunotherapy development and validated tools that can be used for immune profiling of other tumor types.

Influence of node abundance on signaling network state and dynamics analyzed by mass cytometry.

Signaling networks are key regulators of cellular function. Although the concentrations of signaling proteins are perturbed in disease states, such as cancer, and are modulated by drug therapies, our understanding of how such changes shape the properties of signaling networks is limited. Here we couple mass-cytometry-based single-cell analysis with overexpression of tagged signaling proteins to study the dependence of signaling relationships and dynamics on protein node abundance. Focusing on the epidermal growth factor receptor (EGFR) signaling network in HEK293T cells, we analyze 20 signaling proteins during a 1-h EGF stimulation time course using a panel of 35 antibodies. Data analysis with BP-R2, a measure that quantifies complex signaling relationships, reveals abundance-dependent network states and identifies novel signaling relationships. Further, we show that upstream signaling proteins have abundance-dependent effects on downstream signaling dynamics. Our approach elucidates the influence of node abundance on signal transduction networks and will further our understanding of signaling in health and disease.

Systems-level analysis of mechanisms regulating yeast metabolic flux.

Cellular metabolic fluxes are determined by enzyme activities and metabolite abundances. Biochemical approaches reveal the impact of specific substrates or regulators on enzyme kinetics but do not capture the extent to which metabolite and enzyme concentrations vary across physiological states and, therefore, how cellular reactions are regulated. We measured enzyme and metabolite concentrations and metabolic fluxes across 25 steady-state yeast cultures. We then assessed the extent to which flux can be explained by a Michaelis-Menten relationship between enzyme, substrate, product, and potential regulator concentrations. This revealed three previously unrecognized instances of cross-pathway regulation, which we biochemically verified. One of these involved inhibition of pyruvate kinase by citrate, which accumulated and thereby curtailed glycolytic outflow in nitrogen-limited yeast. Overall, substrate concentrations were the strongest driver of the net rates of cellular metabolic reactions, with metabolite concentrations collectively having more than double the physiological impact of enzymes.

Extending the limits of quantitative proteome profiling with data-independent acquisition and application to acetaminophen-treated three-dimensional liver microtissues.

The data-independent acquisition (DIA) approach has recently been introduced as a novel mass spectrometric method that promises to combine the high content aspect of shotgun proteomics with the reproducibility and precision of selected reaction monitoring. Here, we evaluate, whether SWATH-MS type DIA effectively translates into a better protein profiling as compared with the established shotgun proteomics. We implemented a novel DIA method on the widely used Orbitrap platform and used retention-time-normalized (iRT) spectral libraries for targeted data extraction using Spectronaut. We call this combination hyper reaction monitoring (HRM). Using a controlled sample set, we show that HRM outperformed shotgun proteomics both in the number of consistently identified peptides across multiple measurements and quantification of differentially abundant proteins. The reproducibility of HRM in peptide detection was above 98%, resulting in quasi complete data sets compared with 49% of shotgun proteomics. Utilizing HRM, we profiled acetaminophen (APAP)(1)-treated three-dimensional human liver microtissues. An early onset of relevant proteome changes was revealed at subtoxic doses of APAP. Further, we detected and quantified for the first time human NAPQI-protein adducts that might be relevant for the toxicity of APAP. The adducts were identified on four mitochondrial oxidative stress related proteins (GATM, PARK7, PRDX6, and VDAC2) and two other proteins (ANXA2 and FTCD). Our findings imply that DIA should be the preferred method for quantitative protein profiling.

Long-term exposure to bis(2-ethylhexyl)phthalate (DEHP) inhibits growth of guppy fish (Poecilia reticulata).

Bis(2-ethylhexyl)phthalate (DEHP) is a widely used plasticizer that is a commonly found contaminant of aquatic environments. However, little is known about the long-term effects of DEHP on fish development, as previous studies yielded conflicting results and mostly investigated the effects of concentrations higher than those found in natural habitats. We thus aimed to investigate the effects of DHEP (i) at concentrations present in the environment, and (ii) under conditions that might accentuate any deleterious consequences (larvae rather than adult fish, use of higher temperature). Different concentrations of DEHP (0.1-10 microg l(-1)rpar; applied continuously for 91 days were tested on guppy fish that were less than one week old at the beginning of the treatment. As early as 14 days after the start of exposure, guppies treated with 10 microg l(-1) DEHP showed significantly reduced body length as compared with control fish. The inhibitory effect of DEHP was concentration-dependent and increased with time, leading to a maximal reduction in body length of 15 and 40% at 1 and 10 microg l(-1) DEHP, respectively. The effect was even more pronounced for body weight, which was diminished by up to 40 and 70% at 1 and 10 microg l(-1) DEHP, respectively. The reduction in growth was still significant at 91 days of DEHP treatment, whereas the Fulton's condition factor was unaffected. While DEHP significantly blocked growth in both male and female guppies, no shift in the sexual development was observed. These data show that DEHP, at concentrations present in aquatic environments, can profoundly affect development in fish.