Classifying cancer-associated fibroblasts-The good, the bad, and the target.

Cancer-associated fibroblasts (CAFs) are heterogeneous and ubiquitous stromal cells within the tumor microenvironment (TME). Numerous CAF types have been described, typically using single-cell technologies such as single-cell RNA sequencing. There is no general classification system for CAFs, hampering their study and therapeutic targeting. We propose a simple CAF classification system based on single-cell phenotypes and spatial locations of CAFs in multiple cancer types, assess how our scheme fits within current knowledge, and invite the CAF research community to further refine it.

Identification of an embryonic differentiation stage marked by Sox1 and FoxA2 co-expression using combined cell tracking and high dimensional protein imaging.

Pluripotent mouse embryonic stem cells (ESCs) can differentiate to all germ layers and serve as an in vitro model of embryonic development. To better understand the differentiation paths traversed by ESCs committing to different lineages, we track individual differentiating ESCs by timelapse imaging followed by multiplexed high-dimensional Imaging Mass Cytometry (IMC) protein quantification. This links continuous live single-cell molecular NANOG and cellular dynamics quantification over 5-6 generations to protein expression of 37 different molecular regulators in the same single cells at the observation endpoints. Using this unique data set including kinship history and live lineage marker detection, we show that NANOG downregulation occurs generations prior to, but is not sufficient for neuroectoderm marker Sox1 upregulation. We identify a developmental cell type co-expressing both the canonical Sox1 neuroectoderm and FoxA2 endoderm markers in vitro and confirm the presence of such a population in the post-implantation embryo. RNASeq reveals cells co-expressing SOX1 and FOXA2 to have a unique cell state characterized by expression of both endoderm as well as neuroectoderm genes suggesting lineage potential towards both germ layers.

Microenvironmental reorganization in brain tumors following radiotherapy and recurrence revealed by hyperplexed immunofluorescence imaging.

The tumor microenvironment plays a crucial role in determining response to treatment. This involves a series of interconnected changes in the cellular landscape, spatial organization, and extracellular matrix composition. However, assessing these alterations simultaneously is challenging from a spatial perspective, due to the limitations of current high-dimensional imaging techniques and the extent of intratumoral heterogeneity over large lesion areas. In this study, we introduce a spatial proteomic workflow termed Hyperplexed Immunofluorescence Imaging (HIFI) that overcomes these limitations. HIFI allows for the simultaneous analysis of > 45 markers in fragile tissue sections at high magnification, using a cost-effective high-throughput workflow. We integrate HIFI with machine learning feature detection, graph-based network analysis, and cluster-based neighborhood analysis to analyze the microenvironment response to radiation therapy in a preclinical model of glioblastoma, and compare this response to a mouse model of breast-to-brain metastasis. Here we show that glioblastomas undergo extensive spatial reorganization of immune cell populations and structural architecture in response to treatment, while brain metastases show no comparable reorganization. Our integrated spatial analyses reveal highly divergent responses to radiation therapy between brain tumor models, despite equivalent radiotherapy benefit.

Comparative transcriptomics coupled to developmental grading via transgenic zebrafish reporter strains identifies conserved features in neutrophil maturation.

Neutrophils are evolutionarily conserved innate immune cells playing pivotal roles in host defense. Zebrafish models have contributed substantially to our understanding of neutrophil functions but similarities to human neutrophil maturation have not been systematically characterized, which limits their applicability to studying human disease. Here we show, by generating and analysing transgenic zebrafish strains representing distinct neutrophil differentiation stages, a high-resolution transcriptional profile of neutrophil maturation. We link gene expression at each stage to characteristic transcription factors, including C/ebp-ß, which is important for late neutrophil maturation. Cross-species comparison of zebrafish, mouse, and human samples confirms high molecular similarity of immature stages and discriminates zebrafish-specific from pan-species gene signatures. Applying the pan-species neutrophil maturation signature to RNA-sequencing data from human neuroblastoma patients reveals association between metastatic tumor cell infiltration in the bone marrow and an overall increase in mature neutrophils. Our detailed neutrophil maturation atlas thus provides a valuable resource for studying neutrophil function at different stages across species in health and disease.

Multiplex protein imaging in tumour biology.

Tissue imaging has become much more colourful in the past decade. Advances in both experimental and analytical methods now make it possible to image protein markers in tissue samples in high multiplex. The ability to routinely image 40-50 markers simultaneously, at single-cell or subcellular resolution, has opened up new vistas in the study of tumour biology. Cellular phenotypes, interaction, communication and spatial organization have become amenable to molecular-level analysis, and application to patient cohorts has identified clinically relevant cellular and tissue features in several cancer types. Here, we review the use of multiplex protein imaging methods to study tumour biology, discuss ongoing attempts to combine these approaches with other forms of spatial omics, and highlight challenges in the field.

Cancer-associated fibroblast phenotypes are associated with patient outcome in non-small cell lung cancer.

Despite advances in treatment, lung cancer survival rates remain low. A better understanding of the cellular heterogeneity and interplay of cancer-associated fibroblasts (CAFs) within the tumor microenvironment will support the development of personalized therapies. We report a spatially resolved single-cell imaging mass cytometry (IMC) analysis of CAFs in a non-small cell lung cancer cohort of 1,070 patients. We identify four prognostic patient groups based on 11 CAF phenotypes with distinct spatial distributions and show that CAFs are independent prognostic factors for patient survival. The presence of tumor-like CAFs is strongly correlated with poor prognosis. In contrast, inflammatory CAFs and interferon-response CAFs are associated with inflamed tumor microenvironments and higher patient survival. High density of matrix CAFs is correlated with low immune infiltration and is negatively correlated with patient survival. In summary, our data identify phenotypic and spatial features of CAFs that are associated with patient outcome in NSCLC.

cytoviewer: an R/Bioconductor package for interactive visualization and exploration of highly multiplexed imaging data.

BACKGROUND: Highly multiplexed imaging enables single-cell-resolved detection of numerous biological molecules in their spatial tissue context. Interactive visualization of multiplexed imaging data is crucial at any step of data analysis to facilitate quality control and the spatial exploration of single cell features. However, tools for interactive visualization of multiplexed imaging data are not available in the statistical programming language R. RESULTS: Here, we describe cytoviewer, an R/Bioconductor package for interactive visualization and exploration of multi-channel images and segmentation masks. The cytoviewer package supports flexible generation of image composites, allows side-by-side visualization of single channels, and facilitates the spatial visualization of single-cell data in the form of segmentation masks. As such, cytoviewer improves image and segmentation quality control, the visualization of cell phenotyping results and qualitative validation of hypothesis at any step of data analysis. The package operates on standard data classes of the Bioconductor project and therefore integrates with an extensive framework for single-cell and image analysis. The graphical user interface allows intuitive navigation and little coding experience is required to use the package. We showcase the functionality and biological application of cytoviewer by analysis of an imaging mass cytometry dataset acquired from cancer samples. CONCLUSIONS: The cytoviewer package offers a rich set of features for highly multiplexed imaging data visualization in R that seamlessly integrates with the workflow for image and single-cell data analysis. It can be installed from Bioconductor via https://www.bioconductor.org/packages/release/bioc/html/cytoviewer.html . The development version and further instructions can be found on GitHub at https://github.com/BodenmillerGroup/cytoviewer .

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/ .

DNA-barcoded signal amplification for imaging mass cytometry enables sensitive and highly multiplexed tissue imaging.

Imaging mass cytometry (IMC) is a highly multiplexed, antibody-based imaging method that captures heterogeneous spatial protein expression patterns at subcellular resolution. Here we report the extension of IMC to low-abundance markers through incorporation of the DNA-based signal amplification by exchange reaction, immuno-SABER. We applied SABER-IMC to image the tumor immune microenvironment in human melanoma by simultaneous imaging of 18 markers with immuno-SABER and 20 markers without amplification. SABER-IMC enabled the identification of immune cell phenotypic markers, such as T cell co-receptors and their ligands, that are not detectable with IMC.

Multielement Z-tag imaging by X-ray fluorescence microscopy for next-generation multiplex imaging.

Rapid, highly multiplexed, nondestructive imaging that spans the molecular to the supra-cellular scale would be a powerful tool for tissue analysis. However, the physical constraints of established imaging methods limit the simultaneous improvement of these parameters. Whole-organism to atomic-level imaging is possible with tissue-penetrant, picometer-wavelength X-rays. To enable highly multiplexed X-ray imaging, we developed multielement Z-tag X-ray fluorescence (MEZ-XRF) that can operate at kHz speeds when combined with signal amplification by exchange reaction (SABER)-amplified Z-tag reagents. We demonstrated parallel imaging of 20 Z-tag or SABER Z-tag reagents at subcellular resolution in cell lines and multiple human tissues. We benchmarked MEZ-XRF against imaging mass cytometry and demonstrated the nondestructive multiscale repeat imaging capabilities of MEZ-XRF with rapid tissue overview scans, followed by slower, more sensitive imaging of low-abundance markers such as immune checkpoint proteins. The unique multiscale, nondestructive nature of MEZ-XRF, combined with SABER Z-tags for high sensitivity or enhanced speed, enables highly multiplexed bioimaging across biological scales.

Acquired resistance to anti-PD1 therapy in patients with NSCLC associates with immunosuppressive T cell phenotype.

Immune checkpoint inhibitor treatment has the potential to prolong survival in non-small cell lung cancer (NSCLC), however, some of the patients develop resistance following initial response. Here, we analyze the immune phenotype of matching tumor samples from a cohort of NSCLC patients showing good initial response to immune checkpoint inhibitors, followed by acquired resistance at later time points. By using imaging mass cytometry and whole exome and RNA sequencing, we detect two patterns of resistance¨: One group of patients is characterized by reduced numbers of tumor-infiltrating CD8+ T cells and reduced expression of PD-L1 after development of resistance, whereas the other group shows high CD8+ T cell infiltration and high expression of PD-L1 in addition to markedly elevated expression of other immune-inhibitory molecules. In two cases, we detect downregulation of type I and II IFN pathways following progression to resistance, which could lead to an impaired anti-tumor immune response. This study thus captures the development of immune checkpoint inhibitor resistance as it progresses and deepens our mechanistic understanding of immunotherapy response in NSCLC.

Cancer-associated fibroblast classification in single-cell and spatial proteomics data.

Cancer-associated fibroblasts (CAFs) are a diverse cell population within the tumour microenvironment, where they have critical effects on tumour evolution and patient prognosis. To define CAF phenotypes, we analyse a single-cell RNA sequencing (scRNA-seq) dataset of over 16,000 stromal cells from tumours of 14 breast cancer patients, based on which we define and functionally annotate nine CAF phenotypes and one class of pericytes. We validate this classification system in four additional cancer types and use highly multiplexed imaging mass cytometry on matched breast cancer samples to confirm our defined CAF phenotypes at the protein level and to analyse their spatial distribution within tumours. This general CAF classification scheme will allow comparison of CAF phenotypes across studies, facilitate analysis of their functional roles, and potentially guide development of new treatment strategies in the future.

cytoviewer: an R/Bioconductor package for interactive visualization and exploration of highly multiplexed imaging data.

Highly multiplexed imaging enables single-cell-resolved detection of numerous biological molecules in their spatial tissue context. Interactive data visualization of multiplexed imaging data is necessary for quality control and hypothesis examination. Here, we describe cytoviewer, an R/Bioconductor package for interactive visualization and exploration of multi-channel images and segmentation masks. The cytoviewer package supports flexible generation of image composites, allows side-by-side visualization of single channels, and facilitates the spatial visualization of single-cell data in the form of segmentation masks. The package operates on SingleCellExperiment, SpatialExperiment and CytoImageList objects and therefore integrates with the Bioconductor framework for single-cell and image analysis. Users of cytoviewer need little coding expertise, and the graphical user interface allows user-friendly navigation. We showcase the functionality of cytoviewer by analysis of an imaging mass cytometry dataset of cancer patients.

Multiplex imaging of breast cancer lymph node metastases identifies prognostic single-cell populations independent of clinical classifiers.

Although breast cancer mortality is largely caused by metastasis, clinical decisions are based on analysis of the primary tumor and on lymph node involvement but not on the phenotype of disseminated cells. Here, we use multiplex imaging mass cytometry to compare single-cell phenotypes of primary breast tumors and matched lymph node metastases in 205 patients. We observe extensive phenotypic variability between primary and metastatic sites and that disseminated cell phenotypes frequently deviate from the clinical disease subtype. We identify single-cell phenotypes and spatial organizations of disseminated tumor cells that are associated with patient survival and a weaker survival association for high-risk phenotypes in the primary tumor. We show that p53 and GATA3 in lymph node metastases provide prognostic information beyond clinical classifiers and can be measured with standard methods. Molecular characterization of disseminated tumor cells is an untapped source of clinically applicable prognostic information for breast cancer.

A comprehensive single-cell map of T cell exhaustion-associated immune environments in human breast cancer.

Immune checkpoint therapy in breast cancer remains restricted to triple negative patients, and long-term clinical benefit is rare. The primary aim of immune checkpoint blockade is to prevent or reverse exhausted T cell states, but T cell exhaustion in breast tumors is not well understood. Here, we use single-cell transcriptomics combined with imaging mass cytometry to systematically study immune environments of human breast tumors that either do or do not contain exhausted T cells, with a focus on luminal subtypes. We find that the presence of a PD-1high exhaustion-like T cell phenotype is associated with an inflammatory immune environment with a characteristic cytotoxic profile, increased myeloid cell activation, evidence for elevated immunomodulatory, chemotactic, and cytokine signaling, and accumulation of natural killer T cells. Tumors harboring exhausted-like T cells show increased expression of MHC-I on tumor cells and of CXCL13 on T cells, as well as altered spatial organization with more immature rather than mature tertiary lymphoid structures. Our data reveal fundamental differences between immune environments with and without exhausted T cells within luminal breast cancer, and show that expression of PD-1 and CXCL13 on T cells, and MHC-I - but not PD-L1 - on tumor cells are strong distinguishing features between these environments.

Optimizing multiplexed imaging experimental design through tissue spatial segregation estimation.

Recent advances in multiplexed imaging methods allow simultaneous detection of dozens of proteins and hundreds of RNAs, enabling deep spatial characterization of both healthy and diseased tissues. Parameters for the design of optimal multiplex imaging studies, especially those estimating how much area has to be imaged to capture all cell phenotype clusters, are lacking. Here, using a spatial transcriptomic atlas of healthy and tumor human tissues, we developed a statistical framework that determines the number and area of fields of view necessary to accurately identify all cell phenotypes that are part of a tissue. Using this strategy on imaging mass cytometry data, we identified a measurement of tissue spatial segregation that enables optimal experimental design. This strategy will enable an improved design of multiplexed imaging studies.

Imaging and Molecular Annotation of Xenographs and Tumours (IMAXT): High throughput data and analysis infrastructure.

With the aim of producing a 3D representation of tumors, imaging and molecular annotation of xenografts and tumors (IMAXT) uses a large variety of modalities in order to acquire tumor samples and produce a map of every cell in the tumor and its host environment. With the large volume and variety of data produced in the project, we developed automatic data workflows and analysis pipelines. We introduce a research methodology where scientists connect to a cloud environment to perform analysis close to where data are located, instead of bringing data to their local computers. Here, we present the data and analysis infrastructure, discuss the unique computational challenges and describe the analysis chains developed and deployed to generate molecularly annotated tumor models. Registration is achieved by use of a novel technique involving spherical fiducial marks that are visible in all imaging modalities used within IMAXT. The automatic pipelines are highly optimized and allow to obtain processed datasets several times quicker than current solutions narrowing the gap between data acquisition and scientific exploitation.

Vaccination with Designed Neopeptides Induces Intratumoral, Cross-reactive CD4+ T-cell Responses in Glioblastoma.

PURPOSE: The low mutational load of some cancers is considered one reason for the difficulty to develop effective tumor vaccines. To overcome this problem, we developed a strategy to design neopeptides through single amino acid mutations to enhance their immunogenicity. EXPERIMENTAL DESIGN: Exome and RNA sequencing as well as in silico HLA-binding predictions to autologous HLA molecules were used to identify candidate neopeptides. Subsequently, in silico HLA-anchor placements were used to deduce putative T-cell receptor (TCR) contacts of peptides. Single amino acids of TCR contacting residues were then mutated by amino acid replacements. Overall, 175 peptides were synthesized and sets of 25 each containing both peptides designed to bind to HLA class I and II molecules applied in the vaccination. Upon development of a tumor recurrence, the tumor-infiltrating lymphocytes (TIL) were characterized in detail both at the bulk and clonal level. RESULTS: The immune response of peripheral blood T cells to vaccine peptides, including natural peptides and designed neopeptides, gradually increased with repetitive vaccination, but remained low. In contrast, at the time of tumor recurrence, CD8+ TILs and CD4+ TILs responded to 45% and 100%, respectively, of the vaccine peptides. Furthermore, TIL-derived CD4+ T-cell clones showed strong responses and tumor cell lysis not only against the designed neopeptide but also against the unmutated natural peptides of the tumor. CONCLUSIONS: Turning tumor self-peptides into foreign antigens by introduction of designed mutations is a promising strategy to induce strong intratumoral CD4+ T-cell responses in a cold tumor like glioblastoma.

scQUEST: Quantifying tumor ecosystem heterogeneity from mass or flow cytometry data.

With mass and flow cytometry, millions of single-cell profiles with dozens of parameters can be measured to comprehensively characterize complex tumor ecosystems. Here, we present scQUEST, an open-source Python library for cell type identification and quantification of tumor ecosystem heterogeneity in patient cohorts. We provide a step-by-step protocol on the application of scQUEST on our previously generated human breast cancer single-cell atlas using mass cytometry and discuss how it can be adapted and extended for other datasets and analyses. For complete details on the use and execution of this protocol, please refer to Wagner et al. (2019).