Highly multiplexed tissue imaging using repeated oligonucleotide exchange reaction

Multiplex Imaging for Cell Phenotyping of Early Human Atherosclerosis

BACKGROUND

Previous studies using animal models and cultured cells suggest that vascular smooth muscle cells (SMCs) and inflammatory cytokines are important players in atherogenesis. Validating these findings in human disease is critical to designing therapeutics that target these components. Multiplex imaging is a powerful tool for characterizing cell phenotypes and microenvironments using biobanked human tissue sections. However, this technology has not been applied to human atherosclerotic lesions and needs to first be customized and validated.

METHODS AND RESULTS

For validation, we created an 8-plex imaging panel to distinguish foam cells from SMC and leukocyte origins on tissue sections of early human atherosclerotic lesions (n=9). The spatial distribution and characteristics of these foam cells were further analyzed to test the association between SMC phenotypes and inflammation. Consistent with previous reports using human lesions, multiplex imaging showed that foam cells of SMC origin outnumbered those of leukocyte origin and were enriched in the deep intima, where the lipids accumulate in early atherogenesis. This new technology also found that apoptosis or the expression of pro-inflammatory cytokines were not more associated with foam cells than with nonfoam cells in early human lesions. More CD68+ SMCs were present among SMCs that highly expressed interleukin-1β. Highly inflamed SMCs showed a trend of increased apoptosis, whereas leukocytes expressing similar levels of cytokines were enriched in regions of extracellular matrix remodeling.

CONCLUSIONS

The multiplex imaging method can be applied to biobanked human tissue sections to enable proof-of-concept studies and validate theories based on animal models and cultured cells.

The Extra-Islet Pancreas Supports Autoimmunity in Human Type 1 Diabetes

In autoimmune Type 1 diabetes (T1D), immune cells infiltrate and destroy the islets of Langerhans - islands of endocrine tissue dispersed throughout the pancreas. However, the contribution of cellular programs outside islets to insulitis is unclear. Here, using CO-Detection by indEXing (CODEX) tissue imaging and cadaveric pancreas samples, we simultaneously examine islet and extra-islet inflammation in human T1D. We identify four sub-states of inflamed islets characterized by the activation profiles of CD8 + T cells enriched in islets relative to the surrounding tissue. We further find that the extra-islet space of lobules with extensive islet-infiltration differs from the extra-islet space of less infiltrated areas within the same tissue section. Finally, we identify lymphoid structures away from islets enriched in CD45RA + T cells - a population also enriched in one of the inflamed islet sub-states. Together, these data help define the coordination between islets and the extra-islet pancreas in the pathogenesis of human T1D.

Human BioMolecular Atlas Program (HuBMAP): 3D Human Reference Atlas Construction and Usage

The Human BioMolecular Atlas Program (HuBMAP) aims to construct a reference 3D structural, cellular, and molecular atlas of the healthy adult human body. The HuBMAP Data Portal ( https://portal.hubmapconsortium.org ) serves experimental datasets and supports data processing, search, filtering, and visualization. The Human Reference Atlas (HRA) Portal ( https://humanatlas.io ) provides open access to atlas data, code, procedures, and instructional materials. Experts from more than 20 consortia are collaborating to construct the HRA's Common Coordinate Framework (CCF), knowledge graphs, and tools that describe the multiscale structure of the human body (from organs and tissues down to cells, genes, and biomarkers) and to use the HRA to understand changes that occur at each of these levels with aging, disease, and other perturbations. The 6th release of the HRA v2.0 covers 36 organs with 4,499 unique anatomical structures, 1,195 cell types, and 2,089 biomarkers (e.g., genes, proteins, lipids) linked to ontologies. In addition, three workflows were developed to map new experimental data into the HRA's CCF. This paper describes the HRA user stories, terminology, data formats, ontology validation, unified analysis workflows, user interfaces, instructional materials, application programming interface (APIs), flexible hybrid cloud infrastructure, and demonstrates first atlas usage applications and previews.

Cross-domain information fusion for enhanced cell population delineation in single-cell spatial-omics data

Cell population delineation and identification is an essential step in single-cell and spatial-omics studies. Spatial-omics technologies can simultaneously measure information from three complementary domains related to this task: expression levels of a panel of molecular biomarkers at single-cell resolution, relative positions of cells, and images of tissue sections, but existing computational methods for performing this task on single-cell spatial-omics datasets often relinquish information from one or more domains. The additional reliance on the availability of "atlas" training or reference datasets limits cell type discovery to well-defined but limited cell population labels, thus posing major challenges for using these methods in practice. Successful integration of all three domains presents an opportunity for uncovering cell populations that are functionally stratified by their spatial contexts at cellular and tissue levels: the key motivation for employing spatial-omics technologies in the first place. In this work, we introduce Cell Spatio- and Neighborhood-informed Annotation and Patterning (CellSNAP), a self-supervised computational method that learns a representation vector for each cell in tissue samples measured by spatial-omics technologies at the single-cell or finer resolution. The learned representation vector fuses information about the corresponding cell across all three aforementioned domains. By applying CellSNAP to datasets spanning both spatial proteomic and spatial transcriptomic modalities, and across different tissue types and disease settings, we show that CellSNAP markedly enhances de novo discovery of biologically relevant cell populations at fine granularity, beyond current approaches, by fully integrating cells' molecular profiles with cellular neighborhood and tissue image information.

Clonal hematopoiesis and its impact on the aging osteo-hematopoietic niche

Clonal hematopoiesis (CH) defines a premalignant state predominantly found in older persons that increases the risk of developing hematologic malignancies and age-related inflammatory diseases. However, the risk for malignant transformation or non-malignant disorders is variable and difficult to predict, and defining the clinical relevance of specific candidate driver mutations in individual carriers has proved to be challenging. In addition to the cell-intrinsic mechanisms, mutant cells rely on and alter cell-extrinsic factors from the bone marrow (BM) niche, which complicates the prediction of a mutant cell's fate in a shifting pre-malignant microenvironment. Therefore, identifying the insidious and potentially broad impact of driver mutations on supportive niches and immune function in CH aims to understand the subtle differences that enable driver mutations to yield different clinical outcomes. Here, we review the changes in the aging BM niche and the emerging evidence supporting the concept that CH can progressively alter components of the local BM microenvironment. These alterations may have profound implications for the functionality of the osteo-hematopoietic niche and overall bone health, consequently fostering a conducive environment for the continued development and progression of CH. We also provide an overview of the latest technology developments to study the spatiotemporal dependencies in the CH BM niche, ideally in the context of longitudinal studies following CH over time. Finally, we discuss aspects of CH carrier management in clinical practice, based on work from our group and others.

Integrating multiplexed imaging and multiscale modeling identifies tumor phenotype conversion as a critical component of therapeutic T cell efficacy

Cancer progression is a complex process involving interactions that unfold across molecular, cellular, and tissue scales. These multiscale interactions have been difficult to measure and to simulate. Here, we integrated CODEX multiplexed tissue imaging with multiscale modeling software to model key action points that influence the outcome of T cell therapies with cancer. The initial phenotype of therapeutic T cells influences the ability of T cells to convert tumor cells to an inflammatory, anti-proliferative phenotype. This T cell phenotype could be preserved by structural reprogramming to facilitate continual tumor phenotype conversion and killing. One takeaway is that controlling the rate of cancer phenotype conversion is critical for control of tumor growth. The results suggest new design criteria and patient selection metrics for T cell therapies, call for a rethinking of T cell therapeutic implementation, and provide a foundation for synergistically integrating multiplexed imaging data with multiscale modeling of the cancer-immune interface. A record of this paper's transparent peer review process is included in the supplemental information.

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.

Spatial omics techniques and data analysis for cancer immunotherapy applications

In-depth profiling of cancer cells/tissues is expanding our understanding of the genomic, epigenomic, transcriptomic, and proteomic landscape of cancer. However, the complexity of the cancer microenvironment, particularly its immune regulation, has made it difficult to exploit the potential of cancer immunotherapy. High-throughput spatial omics technologies and analysis pipelines have emerged as powerful tools for tackling this challenge. As a result, a potential revolution in cancer diagnosis, prognosis, and treatment is on the horizon. In this review, we discuss the technological advances in spatial profiling of cancer around and beyond the central dogma to harness the full benefits of immunotherapy. We also discuss the promise and challenges of spatial data analysis and interpretation and provide an outlook for the future.

HCR spectral imaging: 10-plex, quantitative, high-resolution RNA and protein imaging in highly autofluorescent samples

Signal amplification based on the mechanism of hybridization chain reaction (HCR) provides a unified framework for multiplex, quantitative, high-resolution imaging of RNA and protein targets in highly autofluorescent samples. With conventional bandpass imaging, multiplexing is typically limited to four or five targets owing to the difficulty in separating signals generated by fluorophores with overlapping spectra. Spectral imaging has offered the conceptual promise of higher levels of multiplexing, but it has been challenging to realize this potential in highly autofluorescent samples, including whole-mount vertebrate embryos. Here, we demonstrate robust HCR spectral imaging with linear unmixing, enabling simultaneous imaging of ten RNA and/or protein targets in whole-mount zebrafish embryos and mouse brain sections. Further, we demonstrate that the amplified and unmixed signal in each of the ten channels is quantitative, enabling accurate and precise relative quantitation of RNA and/or protein targets with subcellular resolution, and RNA absolute quantitation with single-molecule resolution, in the anatomical context of highly autofluorescent samples.

In Vivo Stimulation of Therapeutic Antigen-Specific T Cells in an Artificial Lymph Node Matrix

T cells are critical mediators of antigen-specific immune responses and are common targets for immunotherapy. Biomaterial scaffolds have previously been used to stimulate antigen-presenting cells to elicit antigen-specific immune responses; however, structural and molecular features that directly stimulate and expand naïve, endogenous, tumor-specific T cells in vivo have not been defined. Here, an artificial lymph node (aLN) matrix is created, which consists of an extracellular matrix hydrogel conjugated with peptide-loaded-MHC complex (Signal 1), the co-stimulatory signal anti-CD28 (Signal 2), and a tethered IL-2 (Signal 3), that can bypass challenges faced by other approaches to activate T cells in situ such as vaccines. This dynamic immune-stimulating platform enables direct, in vivo antigen-specific CD8+ T cell stimulation, as well as recruitment and coordination of host immune cells, providing an immuno-stimulatory microenvironment for antigen-specific T cell activation and expansion. Co-injecting the aLN with naïve, wild-type CD8+ T cells results in robust activation and expansion of tumor-targeted T cells that kill target cells and slow tumor growth in several distal tumor models. The aLN platform induces potent in vivo antigen-specific CD8+ T cell stimulation without the need for ex vivo priming or expansion and enables in situ manipulation of antigen-specific responses for immunotherapies.

A replicating LCMV-based vaccine for the treatment of solid tumors

Harnessing the immune system to eradicate tumors requires identification and targeting of tumor antigens, including tumor-specific neoantigens and tumor-associated self-antigens. Tumor-associated antigens are subject to existing immune tolerance, which must be overcome by immunotherapies. Despite many novel immunotherapies reaching clinical trials, inducing self-antigen-specific immune responses remains challenging. Here, we systematically investigate viral-vector-based cancer vaccines encoding a tumor-associated self-antigen (TRP2) for the treatment of established melanomas in preclinical mouse models, alone or in combination with adoptive T cell therapy. We reveal that, unlike foreign antigens, tumor-associated antigens require replication of lymphocytic choriomeningitis virus (LCMV)-based vectors to break tolerance and induce effective antigen-specific CD8+ T cell responses. Immunization with a replicating LCMV vector leads to complete tumor rejection when combined with adoptive TRP2-specific T cell transfer. Importantly, immunization with replicating vectors leads to extended antigen persistence in secondary lymphoid organs, resulting in efficient T cell priming, which renders previously "cold" tumors open to immune infiltration and reprograms the tumor microenvironment to "hot." Our findings have important implications for the design of next-generation immunotherapies targeting solid cancers utilizing viral vectors and adoptive cell transfer.

Unveiling Macrophage Heterogeneity and Their Spatial Distribution Using Multiplexed Tissue Imaging

Macrophages display a high degree of phenotypic diversity and plasticity, which is influenced by their location within the tissue microenvironment. Co-Detection by Indexing (CODEX), a multiplexed imaging technique, allows the simultaneous detection of multiple membrane and cellular markers that enable the accurate identification of tissue-resident hematopoietic and non-hematopoietic cells, while conferring spatial information at a single-cell level. Here we describe the use of CODEX to visualize the phenotypic and spatial heterogeneity of murine tissue-resident macrophages in several organs, and a pipeline to characterize their cellular microenvironments and interactions.

Histological and functional characterization of 3D human skin models mimicking the inflammatory skin diseases psoriasis and atopic dermatitis

Three-dimensional (3D) human skin equivalents have emerged as valuable tools in skin research, replacing animal experimentation and precluding the need for patient biopsies. In this study, we advanced 3D skin equivalents to model the inflammatory skin diseases atopic dermatitis and psoriasis by cytokine stimulation, and were successful in integrating TH1 T cells into skin models to develop an immunocompetent 3D psoriasis model. We performed in-depth histological and functional characterization of 3D skin equivalents and validated them in terms of tissue architecture, pathological changes, expression of antimicrobial peptides and Staphylococcus aureus colonization using 3D reconstruction by multiphoton microscopy and phenotyping by highly multiplexed 'co-detection by indexing' (CODEX) microscopy. We show that our skin equivalents have a structural architecture with a well-developed dermis and epidermis, thus resembling human skin. In addition, the skin models of atopic dermatitis and psoriasis show several phenotypic features of inflammatory skin disease, including disturbed epidermal differentiation and alterations in the expression of epidermal barrier genes and antimicrobial peptides, and can be reliably used to test novel treatment strategies. Therefore, these 3D equivalents will be a valuable tool in experimental dermatological research.

T cell-mediated curation and restructuring of tumor tissue coordinates an effective immune response

Antigen-specific T cells traffic to, are influenced by, and create unique cellular microenvironments. Here we characterize these microenvironments over time with multiplexed imaging in a melanoma model of adoptive T cell therapy and human patients with melanoma treated with checkpoint inhibitor therapy. Multicellular neighborhood analysis reveals dynamic immune cell infiltration and inflamed tumor cell neighborhoods associated with CD8+ T cells. T cell-focused analysis indicates T cells are found along a continuum of neighborhoods that reflect the progressive steps coordinating the anti-tumor immune response. More effective anti-tumor immune responses are characterized by inflamed tumor-T cell neighborhoods, flanked by dense immune infiltration neighborhoods. Conversely, ineffective T cell therapies express anti-inflammatory cytokines, resulting in regulatory neighborhoods, spatially disrupting productive T cell-immune and -tumor interactions. Our study provides in situ mechanistic insights into temporal tumor microenvironment changes, cell interactions critical for response, and spatial correlates of immunotherapy outcomes, informing cellular therapy evaluation and engineering.

Hematoxylin and Eosin staining of PhenoCycler® Fusion flow cell slides

Multiplexed Imaging technologies are powerful techniques that enable ultrahigh-plex spatial phenotyping of whole tissue sections at single cell spatial resolution. Co-Detection by Indexing (CODEX) multiplexing can detect up to 100 proteins using cyclic detection of DNA conjugated antibodies applied to tissue sections. However, it is necessary to correlate multiplexed fluorescent (mIF) spatial images with Hematoxylin and Eosin (H&E) stained sections post analysis. To effectively correlate mIF spatial images with H&E morphology, an (H&E) staining protocol was developed that is directly applied to the CODEX Fusion flow-cell slide after analysis allowing for direct H&E correlation and annotation with mIF images.

NIPMAP: niche-phenotype mapping of multiplex histology data by community ecology

Advances in multiplex histology allow surveying millions of cells, dozens of cell types, and up to thousands of phenotypes within the spatial context of tissue sections. This leads to a combinatorial challenge in (a) summarizing the cellular and phenotypic architecture of tissues and (b) identifying phenotypes with interesting spatial architecture. To address this, we combine ideas from community ecology and machine learning into niche-phenotype mapping (NIPMAP). NIPMAP takes advantage of geometric constraints on local cellular composition imposed by the niche structure of tissues in order to automatically segment tissue sections into niches and their interfaces. Projecting phenotypes on niches and their interfaces identifies previously-reported and previously-unreported spatially-driven phenotypes, concisely summarizes the phenotypic architecture of tissues, and reveals fundamental properties of tissue architecture. NIPMAP is applicable to both protein and RNA multiplex histology of healthy and diseased tissue. An open-source R/Python package implements NIPMAP.

Neutrophil extracellular traps enhance S. aureus skin colonization by oxidative stress induction and downregulation of epidermal barrier genes

Staphylococcus aureus is the most common cause of bacterial skin infections in humans, including patients with atopic dermatitis (AD). Polymorphonuclear neutrophils (PMNs) are the first cells to infiltrate an infection site, where they usually provide an effective first line of defense, including neutrophil extracellular trap (NET) formation. Here, we show that infiltrating PMNs in inflamed human and mouse skin enhance S. aureus skin colonization and persistence. Mechanistically, we demonstrate that a crosstalk between keratinocytes and PMNs results in enhanced NET formation upon S. aureus infection, which in turn induces oxidative stress and expression of danger-associated molecular patterns such as high-mobility-group-protein B1 (HMGB1) in keratinocytes. In turn, HMGB1 enhances S. aureus skin colonization and persistence by promoting skin barrier dysfunctions by the downregulation of epidermal barrier genes. Using patient material, we show that patients with AD exhibit enhanced presence of PMNs, NETs, and HMGB1 in the skin, demonstrating the clinical relevance of our finding.

Integration of spatial and single-cell data across modalities with weakly linked features

Although single-cell and spatial sequencing methods enable simultaneous measurement of more than one biological modality, no technology can capture all modalities within the same cell. For current data integration methods, the feasibility of cross-modal integration relies on the existence of highly correlated, a priori 'linked' features. We describe matching X-modality via fuzzy smoothed embedding (MaxFuse), a cross-modal data integration method that, through iterative coembedding, data smoothing and cell matching, uses all information in each modality to obtain high-quality integration even when features are weakly linked. MaxFuse is modality-agnostic and demonstrates high robustness and accuracy in the weak linkage scenario, achieving 20~70% relative improvement over existing methods under key evaluation metrics on benchmarking datasets. A prototypical example of weak linkage is the integration of spatial proteomic data with single-cell sequencing data. On two example analyses of this type, MaxFuse enabled the spatial consolidation of proteomic, transcriptomic and epigenomic information at single-cell resolution on the same tissue section.

Computational Methods for Single-Cell Proteomics

Advances in single-cell proteomics technologies have resulted in high-dimensional datasets comprising millions of cells that are capable of answering key questions about biology and disease. The advent of these technologies has prompted the development of computational tools to process and visualize the complex data. In this review, we outline the steps of single-cell and spatial proteomics analysis pipelines. In addition to describing available methods, we highlight benchmarking studies that have identified advantages and pitfalls of the currently available computational toolkits. As these technologies continue to advance, robust analysis tools should be developed in tandem to take full advantage of the potential biological insights provided by these data.

A review on deep learning applications in highly multiplexed tissue imaging data analysis

Since its introduction into the field of oncology, deep learning (DL) has impacted clinical discoveries and biomarker predictions. DL-driven discoveries and predictions in oncology are based on a variety of biological data such as genomics, proteomics, and imaging data. DL-based computational frameworks can predict genetic variant effects on gene expression, as well as protein structures based on amino acid sequences. Furthermore, DL algorithms can capture valuable mechanistic biological information from several spatial "omics" technologies, such as spatial transcriptomics and spatial proteomics. Here, we review the impact that the combination of artificial intelligence (AI) with spatial omics technologies has had on oncology, focusing on DL and its applications in biomedical image analysis, encompassing cell segmentation, cell phenotype identification, cancer prognostication, and therapy prediction. We highlight the advantages of using highly multiplexed images (spatial proteomics data) compared to single-stained, conventional histopathological ("simple") images, as the former can provide deep mechanistic insights that cannot be obtained by the latter, even with the aid of explainable AI. Furthermore, we provide the reader with the advantages/disadvantages of DL-based pipelines used in preprocessing highly multiplexed images (cell segmentation, cell type annotation). Therefore, this review also guides the reader to choose the DL-based pipeline that best fits their data. In conclusion, DL continues to be established as an essential tool in discovering novel biological mechanisms when combined with technologies such as highly multiplexed tissue imaging data. In balance with conventional medical data, its role in clinical routine will become more important, supporting diagnosis and prognosis in oncology, enhancing clinical decision-making, and improving the quality of care for patients. Since its introduction into the field of oncology, deep learning (DL) has impacted clinical discoveries and biomarker predictions. DL-driven discoveries and predictions in oncology are based on a variety of biological data such as genomics, proteomics, and imaging data. DL-based computational frameworks can predict genetic variant effects on gene expression, as well as protein structures based on amino acid sequences. Furthermore, DL algorithms can capture valuable mechanistic biological information from several spatial "omics" technologies, such as spatial transcriptomics and spatial proteomics. Here, we review the impact that the combination of artificial intelligence (AI) with spatial omics technologies has had on oncology, focusing on DL and its applications in biomedical image analysis, encompassing cell segmentation, cell phenotype identification, cancer prognostication, and therapy prediction. We highlight the advantages of using highly multiplexed images (spatial proteomics data) compared to single-stained, conventional histopathological ("simple") images, as the former can provide deep mechanistic insights that cannot be obtained by the latter, even with the aid of explainable AI. Furthermore, we provide the reader with the advantages/disadvantages of the DL-based pipelines used in preprocessing the highly multiplexed images (cell segmentation, cell type annotation). Therefore, this review also guides the reader to choose the DL-based pipeline that best fits their data. In conclusion, DL continues to be established as an essential tool in discovering novel biological mechanisms when combined with technologies such as highly multiplexed tissue imaging data. In balance with conventional medical data, its role in clinical routine will become more important, supporting diagnosis and prognosis in oncology, enhancing clinical decision-making, and improving the quality of care for patients.

Organization of the human intestine at single-cell resolution

The intestine is a complex organ that promotes digestion, extracts nutrients, participates in immune surveillance, maintains critical symbiotic relationships with microbiota and affects overall health1. The intesting has a length of over nine metres, along which there are differences in structure and function2. The localization of individual cell types, cell type development trajectories and detailed cell transcriptional programs probably drive these differences in function. Here, to better understand these differences, we evaluated the organization of single cells using multiplexed imaging and single-nucleus RNA and open chromatin assays across eight different intestinal sites from nine donors. Through systematic analyses, we find cell compositions that differ substantially across regions of the intestine and demonstrate the complexity of epithelial subtypes, and find that the same cell types are organized into distinct neighbourhoods and communities, highlighting distinct immunological niches that are present in the intestine. We also map gene regulatory differences in these cells that are suggestive of a regulatory differentiation cascade, and associate intestinal disease heritability with specific cell types. These results describe the complexity of the cell composition, regulation and organization for this organ, and serve as an important reference map for understanding human biology and disease.

Concerted epithelial and stromal changes during progression of Barrett's Esophagus to invasive adenocarcinoma exposed by multi-scale, multi-omics analysis

Esophageal adenocarcinoma arises from Barrett's esophagus, a precancerous metaplastic replacement of squamous by columnar epithelium in response to chronic inflammation. Multi-omics profiling, integrating single-cell transcriptomics, extracellular matrix proteomics, tissue-mechanics and spatial proteomics of 64 samples from 12 patients' paths of progression from squamous epithelium through metaplasia, dysplasia to adenocarcinoma, revealed shared and patient-specific progression characteristics. The classic metaplastic replacement of epithelial cells was paralleled by metaplastic changes in stromal cells, ECM and tissue stiffness. Strikingly, this change in tissue state at metaplasia was already accompanied by appearance of fibroblasts with characteristics of carcinoma-associated fibroblasts and of an NK cell-associated immunosuppressive microenvironment. Thus, Barrett's esophagus progresses as a coordinated multi-component system, supporting treatment paradigms that go beyond targeting cancerous cells to incorporating stromal reprogramming.

Inflammation in multiple sclerosis: consequences for remyelination and disease progression

Despite the large number of immunomodulatory or immunosuppressive treatments available to treat relapsing-remitting multiple sclerosis (MS), treatment of the progressive phase of the disease has not yet been achieved. This lack of successful treatment approaches is caused by our poor understanding of the mechanisms driving disease progression. Emerging concepts suggest that a combination of persisting focal and diffuse inflammation within the CNS and a gradual failure of compensatory mechanisms, including remyelination, result in disease progression. Therefore, promotion of remyelination presents a promising intervention approach. However, despite our increasing knowledge regarding the cellular and molecular mechanisms regulating remyelination in animal models, therapeutic increases in remyelination remain an unmet need in MS, which suggests that mechanisms of remyelination and remyelination failure differ fundamentally between humans and demyelinating animal models. New and emerging technologies now allow us to investigate the cellular and molecular mechanisms underlying remyelination failure in human tissue samples in an unprecedented way. The aim of this Review is to summarize our current knowledge regarding mechanisms of remyelination and remyelination failure in MS and in animal models of the disease, identify open questions, challenge existing concepts, and discuss strategies to overcome the translational roadblock in the field of remyelination-promoting therapies.

Nano-DESI Mass Spectrometry Imaging of Proteoforms in Biological Tissues with High Spatial Resolution

Mass spectrometry imaging (MSI) is a powerful tool for label-free mapping of the spatial distribution of proteins in biological tissues. We have previously demonstrated imaging of individual proteoforms in biological tissues using nanospray desorption electrospray ionization (nano-DESI), an ambient liquid extraction-based MSI technique. Nano-DESI MSI generates multiply charged protein ions, which is advantageous for their identification using top-down proteomics analysis. In this study, we demonstrate proteoform mapping in biological tissues with a spatial resolution down to 7 μm using nano-DESI MSI. A substantial decrease in protein signals observed in high-spatial-resolution MSI makes these experiments challenging. We have enhanced the sensitivity of nano-DESI MSI experiments by optimizing the design of the capillary-based probe and the thickness of the tissue section. In addition, we demonstrate that oversampling may be used to further improve spatial resolution at little or no expense to sensitivity. These developments represent a new step in MSI-based spatial proteomics, which complements targeted imaging modalities widely used for studying biological systems.

Spatial profiling technologies illuminate the tumor microenvironment

The tumor microenvironment (TME) is composed of many different cellular and acellular components that together drive tumor growth, invasion, metastasis, and response to therapies. Increasing realization of the significance of the TME in cancer biology has shifted cancer research from a cancer-centric model to one that considers the TME as a whole. Recent technological advancements in spatial profiling methodologies provide a systematic view and illuminate the physical localization of the components of the TME. In this review, we provide an overview of major spatial profiling technologies. We present the types of information that can be extracted from these data and describe their applications, findings and challenges in cancer research. Finally, we provide a future perspective of how spatial profiling could be integrated into cancer research to improve patient diagnosis, prognosis, stratification to treatment and development of novel therapeutics.

High-multiplex tissue imaging in routine pathology-are we there yet?

High-multiplex tissue imaging (HMTI) approaches comprise several novel immunohistological methods that enable in-depth, spatial single-cell analysis. Over recent years, studies in tumor biology, infectious diseases, and autoimmune conditions have demonstrated the information gain accessible when mapping complex tissues with HMTI. Tumor biology has been a focus of innovative multiparametric approaches, as the tumor microenvironment (TME) contains great informative value for accurate diagnosis and targeted therapeutic approaches: unraveling the cellular composition and structural organization of the TME using sophisticated computational tools for spatial analysis has produced histopathologic biomarkers for outcomes in breast cancer, predictors of positive immunotherapy response in melanoma, and histological subgroups of colorectal carcinoma. Integration of HMTI technologies into existing clinical workflows such as molecular tumor boards will contribute to improve patient outcomes through personalized treatments tailored to the specific heterogeneous pathological fingerprint of cancer, autoimmunity, or infection. Here, we review the advantages and limitations of existing HMTI technologies and outline how spatial single-cell data can improve our understanding of pathological disease mechanisms and determinants of treatment success. We provide an overview of the analytic processing and interpretation and discuss how HMTI can improve future routine clinical diagnostic and therapeutic processes.

Integration of spatial and single-cell data across modalities with weak linkage

single-cell sequencing methods have enabled the profiling of multiple types of molecular readouts at cellular resolution, and recent developments in spatial barcoding, in situ hybridization, and in situ sequencing allow such molecular readouts to retain their spatial context. Since no technology can provide complete characterization across all layers of biological modalities within the same cell, there is pervasive need for computational cross-modal integration (also called diagonal integration) of single-cell and spatial omics data. For current methods, the feasibility of cross-modal integration relies on the existence of highly correlated, a priori "linked" features. When such linked features are few or uninformative, a scenario that we call "weak linkage", existing methods fail. We developed MaxFuse, a cross-modal data integration method that, through iterative co-embedding, data smoothing, and cell matching, leverages all information in each modality to obtain high-quality integration. MaxFuse is modality-agnostic and, through comprehensive benchmarks on single-cell and spatial ground-truth multiome datasets, demonstrates high robustness and accuracy in the weak linkage scenario. A prototypical example of weak linkage is the integration of spatial proteomic data with single-cell sequencing data. On two example analyses of this type, we demonstrate how MaxFuse enables the spatial consolidation of proteomic, transcriptomic and epigenomic information at single-cell resolution on the same tissue section.

A spatial map of human macrophage niches reveals context-dependent macrophage functions in colon and breast cancer

Tumor-associated macrophages (TAMs) display heterogeneous phenotypes. Yet the exact tissue cues that shape macrophage functional diversity are incompletely understood. Here we discriminate, spatially resolve and reveal the function of five distinct macrophage niches within malignant and benign breast and colon tissue. We found that SPP1 TAMs reside in hypoxic and necrotic tumor regions, and a novel subset of FOLR2 tissue resident macrophages (TRMs) supports the plasma cell tissue niche. We discover that IL4I1 macrophages populate niches with high cell turnover where they phagocytose dying cells. Significantly, IL4I1 TAMs abundance correlates with anti-PD1 treatment response in breast cancer. Furthermore, NLRP3 inflammasome activation in NLRP3 TAMs correlates with neutrophil infiltration in the tumors and is associated with poor outcome in breast cancer patients. This suggests the NLRP3 inflammasome as a novel cancer immunetherapy target. Our work uncovers context-dependent roles of macrophage subsets, and suggests novel predictive markers and macrophage subset-specific therapy targets.

A flexible open-source processing workflow for multiplexed fluorescence imaging based on cycles

Background  Multiplexing tissue imaging is developing as a complement for single cell analysis, bringing the spatial information of cells in tissue in addition to multiple parameters measurements. More and more commercial or home-made systems are available. These techniques allow the imaging of tens of fluorescent reporters, where the spectral overlap is solved by imaging by cycles the fluorophores using microfluidics to change the reporters between each cycle. Methods  For several systems, the acquisition system coupled to the microfluidic system is a wide field microscope, and the acquisition process is done by mosaicking to cover a large field of view, relying on image processing to obtain the data set to be analysed in intensity. The processed data set allows the identification of different populations, quite similarly to cytometry analysis, but with spatial information in addition. To obtain the final image for analysis from the raw acquisitions, several preprocessing steps are needed for inter-cycle registration, tissue autofluorescence correction or mosaicking. We propose a workflow for this preprocessing, implemented as an open source software (as a library, command line tool and standalone). Results  We exemplify the workflow on the commercial system PhenoCycler TM (formerly named CODEX®) and provide a reduced size data set for testing. Conclusions  We compare our processor with the commercially provided processor and show that we solve some problems also reported by other users.

A real-time GPU-accelerated parallelized image processor for large-scale multiplexed fluorescence microscopy data

Highly multiplexed, single-cell imaging has revolutionized our understanding of spatial cellular interactions associated with health and disease. With ever-increasing numbers of antigens, region sizes, and sample sizes, multiplexed fluorescence imaging experiments routinely produce terabytes of data. Fast and accurate processing of these large-scale, high-dimensional imaging data is essential to ensure reliable segmentation and identification of cell types and for characterization of cellular neighborhoods and inference of mechanistic insights. Here, we describe RAPID, a Real-time, GPU-Accelerated Parallelized Image processing software for large-scale multiplexed fluorescence microscopy Data. RAPID deconvolves large-scale, high-dimensional fluorescence imaging data, stitches and registers images with axial and lateral drift correction, and minimizes tissue autofluorescence such as that introduced by erythrocytes. Incorporation of an open source CUDA-driven, GPU-assisted deconvolution produced results similar to fee-based commercial software. RAPID reduces data processing time and artifacts and improves image contrast and signal-to-noise compared to our previous image processing pipeline, thus providing a useful tool for accurate and robust analysis of large-scale, multiplexed, fluorescence imaging data.

[Characterization of the tumor microenvironment by highly multiplexed microscopy]

Malignant neoplasms are highly complex ecosystems consisting of tumor cells and the tumor microenvironment (TME), which is composed of structural elements (vessels, fibroblasts, extracellular matrix) and a wide variety of infiltrating immune cell types of the innate and adaptive immune systems. The TME is the main site of tumor cell-immune cell interactions and plays a critical role in antitumoral immunity. Immunotherapies can affect the interactions between immune cell types and tumor cells in the TME, boost immune responses, and lead to tumor elimination. Novel highly multiplexed microscopy techniques, which enable the detection of more than 50 simultaneous markers in tissues, facilitate the in-depth characterization of the TME at single-cell resolution in clinically relevant samples. Detailed knowledge about the cellular and spatial composition of the TME, the specific cell types and their functional properties, and cell-cell interactions-as revealed by highly multiplexed microscopy-will improve our understanding of immunotherapies' mechanisms of action and reveal new potential therapeutic targets and predictive biomarkers.

Spatially Resolved and Highly Multiplexed Protein and RNA In Situ Detection by Combining CODEX With RNAscope In Situ Hybridization

Highly multiplexed protein and RNA in situ detection on a single tissue section concurrently is highly desirable for both basic and applied biomedical research. CO-detection by inDEXing (CODEX) is a new and powerful platform to visualize up to 60 protein biomarkers in situ, and RNAscope in situ hybridization (RNAscope) is a novel RNA detection system with high sensitivity and unprecedent specificity at a single-cell level. Nevertheless, to our knowledge, the combination of CODEX and RNAscope remained unreported until this study. Here, we report a simple and reproducible combination of CODEX and RNAscope. We also determined the cross-reactivities of CODEX anti-human antibodies to rhesus macaques, a widely used animal model of human disease.

Immunotherapy of glioblastoma explants induces interferon-γ responses and spatial immune cell rearrangements in tumor center, but not periphery

A patient-tailored, ex vivo drug response platform for glioblastoma (GBM) would facilitate therapy planning, provide insights into treatment-induced mechanisms in the immune tumor microenvironment (iTME), and enable the discovery of biomarkers of response. We cultured regionally annotated GBM explants in perfusion bioreactors to assess iTME responses to immunotherapy. Explants were treated with anti-CD47, anti-PD-1, or their combination, and analyzed by multiplexed microscopy [CO-Detection by indEXing (CODEX)], enabling the spatially resolved identification of >850,000 single cells, accompanied by explant secretome interrogation. Center and periphery explants differed in their cell type and soluble factor composition, and responses to immunotherapy. A subset of explants displayed increased interferon-γ levels, which correlated with shifts in immune cell composition within specified tissue compartments. Our study demonstrates that ex vivo immunotherapy of GBM explants enables an active antitumoral immune response within the tumor center and provides a framework for multidimensional personalized assessment of tumor response to immunotherapy.

A multiplexed epitope barcoding strategy that enables dynamic cellular phenotypic screens

Pooled genetic libraries have improved screening throughput for mapping genotypes to phenotypes. However, selectable phenotypes are limited, restricting screening to outcomes with a low spatiotemporal resolution. Here, we integrated live-cell imaging with pooled library-based screening. To enable intracellular multiplexing, we developed a method called EPICode that uses a combination of short epitopes, which can also appear in various subcellular locations. EPICode thus enables the use of live-cell microscopy to characterize a phenotype of interest over time, including after sequential stimulatory/inhibitory manipulations, and directly connects behavior to the cellular genotype. To test EPICode's capacity against an important milestone-engineering and optimizing dynamic, live-cell reporters-we developed a live-cell PKA kinase translocation reporter with improved sensitivity and specificity. The use of epitopes as fluorescent barcodes introduces a scalable strategy for high-throughput screening broadly applicable to protein engineering and drug discovery settings where image-based phenotyping is desired.

Cellular and molecular interrogation of kidney biopsy specimens

PURPOSE OF REVIEW

Traditional histopathology of the kidney biopsy specimen has been an essential and successful tool for the diagnosis and staging of kidney diseases. However, it is likely that the full potential of the kidney biopsy has not been tapped so far. Indeed, there is now a concerted worldwide effort to interrogate kidney biopsy samples at the cellular and molecular levels with unprecedented rigor and depth. This review examines these novel approaches to study kidney biopsy specimens and highlights their potential to refine our understanding of the pathophysiology of kidney disease and lead to precision-based diagnosis and therapy.

RECENT FINDINGS

Several consortia are now active at studying kidney biopsy samples from various patient cohorts with state-of-the art cellular and molecular techniques. These include advanced imaging approaches as well as deep molecular interrogation with tools such as epigenetics, transcriptomics, proteomics and metabolomics. The emphasis throughout is on rigor, reproducibility and quality control.

SUMMARY

Although these techniques to study kidney biopsies are complementary, each on its own can yield novel ways to define and classify kidney disease. Therefore, great efforts are needed in order to generate an integrated output that can propel the diagnosis and treatment of kidney disease into the realm of precision medicine.

IBEX: an iterative immunolabeling and chemical bleaching method for high-content imaging of diverse tissues

High-content imaging is needed to catalog the variety of cellular phenotypes and multicellular ecosystems present in metazoan tissues. We recently developed iterative bleaching extends multiplexity (IBEX), an iterative immunolabeling and chemical bleaching method that enables multiplexed imaging (>65 parameters) in diverse tissues, including human organs relevant for international consortia efforts. IBEX is compatible with >250 commercially available antibodies and 16 unique fluorophores, and can be easily adopted to different imaging platforms using slides and nonproprietary imaging chambers. The overall protocol consists of iterative cycles of antibody labeling, imaging and chemical bleaching that can be completed at relatively low cost in 2-5 d by biologists with basic laboratory skills. To support widespread adoption, we provide extensive details on tissue processing, curated lists of validated antibodies and tissue-specific panels for multiplex imaging. Furthermore, instructions are included on how to automate the method using competitively priced instruments and reagents. Finally, we present a software solution for image alignment that can be executed by individuals without programming experience using open-source software and freeware. In summary, IBEX is a noncommercial method that can be readily implemented by academic laboratories and scaled to achieve high-content mapping of diverse tissues in support of a Human Reference Atlas or other such applications.

FLINO: a new method for immunofluorescence bioimage normalization

MOTIVATION

Multiplexed immunofluorescence bioimaging of single-cells and their spatial organization in tissue holds great promise to the development of future precision diagnostics and therapeutics. Current multiplexing pipelines typically involve multiple rounds of immunofluorescence staining across multiple tissue slides. This introduces experimental batch effects that can hide underlying biological signal. It is important to have robust algorithms that can correct for the batch effects while not introducing biases into the data. Performance of data normalization methods can vary among different assay pipelines. To evaluate differences, it is critical to have a ground truth dataset that is representative of the assay.

RESULTS

A new immunoFLuorescence Image NOrmalization method is presented and evaluated against alternative methods and workflows. Multiround immunofluorescence staining of the same tissue with the nuclear dye DAPI was used to represent virtual slides and a ground truth. DAPI was restained on a given tissue slide producing multiple images of the same underlying structure but undergoing multiple representative tissue handling steps. This ground truth dataset was used to evaluate and compare multiple normalization methods including median, quantile, smooth quantile, median ratio normalization and trimmed mean of the M-values. These methods were applied in both an unbiased grid object and segmented cell object workflow to 24 multiplexed biomarkers. An upper quartile normalization of grid objects in log space was found to obtain almost equivalent performance to directly normalizing segmented cell objects by the middle quantile. The developed grid-based technique was then applied with on-slide controls for evaluation. Using five or fewer controls per slide can introduce biases into the data. Ten or more on-slide controls were able to robustly correct for batch effects.

AVAILABILITY AND IMPLEMENTATION

The data underlying this article along with the FLINO R-scripts used to perform the evaluation of image normalizations methods and workflows can be downloaded from https://github.com/GE-Bio/FLINO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Rhesus Macaque CODEX Multiplexed Immunohistochemistry Panel for Studying Immune Responses During Ebola Infection

Non-human primate (NHP) animal models are an integral part of the drug research and development process. For some biothreat pathogens, animal model challenge studies may offer the only possibility to evaluate medical countermeasure efficacy. A thorough understanding of host immune responses in such NHP models is therefore vital. However, applying antibody-based immune characterization techniques to NHP models requires extensive reagent development for species compatibility. In the case of studies involving high consequence pathogens, further optimization for use of inactivated samples may be required. Here, we describe the first optimized CO-Detection by indEXing (CODEX) multiplexed tissue imaging antibody panel for deep profiling of spatially resolved single-cell immune responses in rhesus macaques. This 21-marker panel is composed of a set of 18 antibodies that stratify major immune cell types along with a set three Ebola virus (EBOV)-specific antibodies. We validated these two sets of markers using immunohistochemistry and CODEX in fully inactivated Formalin-Fixed Paraffin-Embedded (FFPE) tissues from mock and EBOV challenged macaques respectively and provide an efficient framework for orthogonal validation of multiple antibody clones using CODEX multiplexed tissue imaging. We also provide the antibody clones and oligonucleotide tag sequences as a valuable resource for other researchers to recreate this reagent set for future studies of tissue immune responses to EBOV infection and other diseases.

Multiplexed imaging and automated signal quantification in formalin-fixed paraffin-embedded tissues by ChipCytometry

Deciphering the spatial composition of cells in tissues is essential for detailed understanding of biological processes in health and disease. Recent technological advances enabled the assessment of the enormous complexity of tissue-derived parameters by highly multiplexed tissue imaging (HMTI), but elaborate machinery and data analyses are required. This severely limits broad applicability of HMTI. Here we demonstrate for the first time the application of ChipCytometry technology, which has unique features for widespread use, on formalin-fixed paraffin-embedded samples, the most commonly used storage technique of clinically relevant patient specimens worldwide. The excellent staining quality permits workflows for automated quantification of signal intensities, which we further optimized to compensate signal spillover from neighboring cells. In combination with the high number of validated markers, the reported platform can be used from unbiased analyses of tissue composition to detection of phenotypically complex rare cells, and can be easily implemented in both routine research and clinical pathology.

Immune cell topography predicts response to PD-1 blockade in cutaneous T cell lymphoma

Cutaneous T cell lymphomas (CTCL) are rare but aggressive cancers without effective treatments. While a subset of patients derive benefit from PD-1 blockade, there is a critically unmet need for predictive biomarkers of response. Herein, we perform CODEX multiplexed tissue imaging and RNA sequencing on 70 tumor regions from 14 advanced CTCL patients enrolled in a pembrolizumab clinical trial (NCT02243579). We find no differences in the frequencies of immune or tumor cells between responders and non-responders. Instead, we identify topographical differences between effector PD-1⁺ CD4⁺ T cells, tumor cells, and immunosuppressive Tregs, from which we derive a spatial biomarker, termed the SpatialScore, that correlates strongly with pembrolizumab response in CTCL. The SpatialScore coincides with differences in the functional immune state of the tumor microenvironment, T cell function, and tumor cell-specific chemokine recruitment and is validated using a simplified, clinically accessible tissue imaging platform. Collectively, these results provide a paradigm for investigating the spatial balance of effector and suppressive T cell activity and broadly leveraging this biomarker approach to inform the clinical use of immunotherapies.

Tissue schematics map the specialization of immune tissue motifs and their appropriation by tumors

A schematic of a biological system, i.e., a representation of its pieces, how they are combined, and what they do, would facilitate understanding its essential organization and alteration in pathogenesis or evolution. We present a computational approach for constructing tissue schematics (TSs) from high-parameter imaging data and a biological model for interpreting them. TSs map the spatial assembly of cellular neighborhoods into tissue motifs, whose modular composition, we propose, enables the generation of complex outputs. We developed our approach in human lymphoid tissue (HLT), identifying the follicular outer zone as a potential relay between neighboring zones and a core lymphoid assembly with modifications characteristic of each HLT type. Applying the TS approach to the tumor microenvironment in human colorectal cancer identified a higher-order motif, whose mutated assembly was negatively associated with patient survival. TSs may therefore elucidate how immune architectures can be specialized and become vulnerable to reprogramming by tumors.

Strategies for Accurate Cell Type Identification in CODEX Multiplexed Imaging Data

Multiplexed imaging is a recently developed and powerful single-cell biology research tool. However, it presents new sources of technical noise that are distinct from other types of single-cell data, necessitating new practices for single-cell multiplexed imaging processing and analysis, particularly regarding cell-type identification. Here we created single-cell multiplexed imaging datasets by performing CODEX on four sections of the human colon (ascending, transverse, descending, and sigmoid) using a panel of 47 oligonucleotide-barcoded antibodies. After cell segmentation, we implemented five different normalization techniques crossed with four unsupervised clustering algorithms, resulting in 20 unique cell-type annotations for the same dataset. We generated two standard annotations: hand-gated cell types and cell types produced by over-clustering with spatial verification. We then compared these annotations at four levels of cell-type granularity. First, increasing cell-type granularity led to decreased labeling accuracy; therefore, subtle phenotype annotations should be avoided at the clustering step. Second, accuracy in cell-type identification varied more with normalization choice than with clustering algorithm. Third, unsupervised clustering better accounted for segmentation noise during cell-type annotation than hand-gating. Fourth, Z-score normalization was generally effective in mitigating the effects of noise from single-cell multiplexed imaging. Variation in cell-type identification will lead to significant differential spatial results such as cellular neighborhood analysis; consequently, we also make recommendations for accurately assigning cell-type labels to CODEX multiplexed imaging.

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.

Histological Studies of the Ventricular-Subventricular Zone as Neural Stem Cell and Glioma Stem Cell Niche

The neural stem cell niche of the ventricular-subventricular zone supports the persistence of stem and progenitor cells in the mature brain. This niche has many notable cytoarchitectural features that affect the activity of stem cells and may also support the survival and growth of invading tumor cells. Histochemical studies of the niche have revealed many proteins that, in combination, can help to reveal stem-like cells in the normal or cancer context, although many caveats persist in the quest to consistently identify these cells in the human brain. Here, we explore the complex relationship between the persistent proliferative capacity of the neural stem cell niche and the malignant proliferation of brain tumors, with a special focus on histochemical identification of stem cells and stem-like tumor cells and an eye toward the potential application of high-dimensional imaging approaches to the field.

Highly Multiplexed Phenotyping of Immunoregulatory Proteins in the Tumor Microenvironment by CODEX Tissue Imaging

Immunotherapies are revolutionizing cancer treatment by boosting the natural ability of the immune system. In addition to antibodies against traditional checkpoint molecules or their ligands (i.e., CTLA-4, PD-1, and PD-L1), therapies targeting molecules such as ICOS, IDO-1, LAG-3, OX40, TIM-3, and VISTA are currently in clinical trials. To better inform clinical care and the design of therapeutic combination strategies, the co-expression of immunoregulatory proteins on individual immune cells within the tumor microenvironment must be robustly characterized. Highly multiplexed tissue imaging platforms, such as CO-Detection by indEXing (CODEX), are primed to meet this need by enabling >50 markers to be simultaneously analyzed in single-cells on formalin-fixed paraffin-embedded (FFPE) tissue sections. Assembly and validation of antibody panels is particularly challenging, with respect to the specificity of antigen detection and robustness of signal over background. Herein, we report the design, development, optimization, and application of a 56-marker CODEX antibody panel to eight cutaneous T cell lymphoma (CTCL) patient samples. This panel is comprised of structural, tumor, and immune cell markers, including eight immunoregulatory proteins that are approved or currently undergoing clinical trials as immunotherapy targets. Here we provide a resource to enable extensive high-dimensional, spatially resolved characterization of the tissue microenvironment across tumor types and imaging modalities. This framework provides researchers with a readily applicable blueprint to study tumor immunology, tissue architecture, and enable mechanistic insights into immunotherapeutic targets.