SpatialCells: automated profiling of tumor microenvironments with spatially resolved multiplexed single-cell data.

Cancer is a complex cellular ecosystem where malignant cells coexist and interact with immune, stromal and other cells within the tumor microenvironment (TME). Recent technological advancements in spatially resolved multiplexed imaging at single-cell resolution have led to the generation of large-scale and high-dimensional datasets from biological specimens. This underscores the necessity for automated methodologies that can effectively characterize molecular, cellular and spatial properties of TMEs for various malignancies. This study introduces SpatialCells, an open-source software package designed for region-based exploratory analysis and comprehensive characterization of TMEs using multiplexed single-cell data. The source code and tutorials are available at https://semenovlab.github.io/SpatialCells. SpatialCells efficiently streamlines the automated extraction of features from multiplexed single-cell data and can process samples containing millions of cells. Thus, SpatialCells facilitates subsequent association analyses and machine learning predictions, making it an essential tool in advancing our understanding of tumor growth, invasion and metastasis.

Multiplexed 3D Analysis of Immune States and Niches in Human Tissue.

Tissue homeostasis and the emergence of disease are controlled by changes in the proportions of resident and recruited cells, their organization into cellular neighbourhoods, and their interactions with acellular tissue components. Highly multiplexed tissue profiling (spatial omics) 1 makes it possible to study this microenvironment in situ , usually in 4-5 micron thick sections (the standard histopathology format) 2 . Microscopy-based tissue profiling is commonly performed at a resolution sufficient to determine cell types but not to detect subtle morphological features associated with cytoskeletal reorganisation, juxtracrine signalling, or membrane trafficking 3 . Here we describe a high-resolution 3D imaging approach able to characterize a wide variety of organelles and structures at sub-micron scale while simultaneously quantifying millimetre-scale spatial features. This approach combines cyclic immunofluorescence (CyCIF) imaging 4 of over 50 markers with confocal microscopy of archival human tissue thick enough (30-40 microns) to fully encompass two or more layers of intact cells. 3D imaging of entire cell volumes substantially improves the accuracy of cell phenotyping and allows cell proximity to be scored using plasma membrane apposition, not just nuclear position. In pre-invasive melanoma in situ 5 , precise phenotyping shows that adjacent melanocytic cells are plastic in state and participate in tightly localised niches of interferon signalling near sites of initial invasion into the underlying dermis. In this and metastatic melanoma, mature and precursor T cells engage in an unexpectedly diverse array of juxtracrine and membrane-membrane interactions as well as looser "neighbourhood" associations 6 whose morphologies reveal functional states. These data provide new insight into the transitions occurring during early tumour formation and immunoediting and demonstrate the potential for phenotyping of tissues at a level of detail previously restricted to cultured cells and organoids.

Development and validation of time-to-event models to predict metastatic recurrence of localized cutaneous melanoma.

BACKGROUND: The recent expansion of immunotherapy for stage IIB/IIC melanoma highlights a growing clinical need to identify patients at high risk of metastatic recurrence and, therefore, most likely to benefit from this therapeutic modality. OBJECTIVE: To develop time-to-event risk prediction models for melanoma metastatic recurrence. METHODS: Patients diagnosed with stage I/II primary cutaneous melanoma between 2000 and 2020 at Mass General Brigham and Dana-Farber Cancer Institute were included. Melanoma recurrence date and type were determined by chart review. Thirty clinicopathologic factors were extracted from electronic health records. Three types of time-to-event machine-learning models were evaluated internally and externally in the distant versus locoregional/nonrecurrence prediction. RESULTS: This study included 954 melanomas (155 distant, 163 locoregional, and 636 1:2 matched nonrecurrences). Distant recurrences were associated with worse survival compared to locoregional/nonrecurrences (HR: 6.21, P < .001) and to locoregional recurrences only (HR: 5.79, P < .001). The Gradient Boosting Survival model achieved the best performance (concordance index: 0.816; time-dependent AUC: 0.842; Brier score: 0.103) in the external validation. LIMITATIONS: Retrospective nature and cohort from one geography. CONCLUSIONS: These results suggest that time-to-event machine-learning models can reliably predict the metastatic recurrence from localized melanoma and help identify high-risk patients who are most likely to benefit from immunotherapy.

SpatialCells: Automated Profiling of Tumor Microenvironments with Spatially Resolved Multiplexed Single-Cell Data.

Background: Cancer is a complex cellular ecosystem where malignant cells coexist and interact with immune, stromal, and other cells within the tumor microenvironment. Recent technological advancements in spatially resolved multiplexed imaging at single-cell resolution have led to the generation of large-scale and high-dimensional datasets from biological specimens. This underscores the necessity for automated methodologies that can effectively characterize the molecular, cellular, and spatial properties of tumor microenvironments for various malignancies. Results: This study introduces SpatialCells, an open-source software package designed for region-based exploratory analysis and comprehensive characterization of tumor microenvironments using multiplexed single-cell data. Conclusions: SpatialCells efficiently streamlines the automated extraction of features from multiplexed single-cell data and can process samples containing millions of cells. Thus, SpatialCells facilitates subsequent association analyses and machine learning predictions, making it an essential tool in advancing our understanding of tumor growth, invasion, and metastasis.

Immune Profiling of Dermatologic Adverse Events from Checkpoint Blockade using Tissue Cyclic Immunofluorescence.

In this study, we demonstrate the utility of whole-slide CyCIF (tissue-based cyclic immunofluorescence) imaging for characterizing immune cell infiltrates in immune checkpoint inhibitor (ICI)-induced dermatologic adverse events (dAEs). We analyzed six cases of ICI-induced dAEs, including lichenoid, bullous pemphigoid, psoriasis, and eczematous eruptions, comparing immune profiling results obtained using both standard immunohistochemistry (IHC) and CyCIF. Our findings indicate that CyCIF provides more detailed and precise single-cell characterization of immune cell infiltrates than IHC, which relies on semi-quantitative scoring by pathologists. This pilot study highlights the potential of CyCIF to advance our understanding of the immune environment in dAEs by revealing tissue-level spatial patterns of immune cell infiltrates, allowing for more precise phenotypic distinctions and deeper exploration of disease mechanisms. By demonstrating that CyCIF can be performed on friable tissues, such as bullous pemphigoid, we provide a foundation for future studies to examine the drivers of specific dAEs using larger cohorts of phenotyped toxicity and suggest a broader role for highly multiplexed tissue imaging in phenotyping the immune mediated disease that they resemble.

Spatial intra-tumor heterogeneity is associated with survival of lung adenocarcinoma patients.

Intra-tumor heterogeneity (ITH) of human tumors is important for tumor progression, treatment response, and drug resistance. However, the spatial distribution of ITH remains incompletely understood. Here, we present spatial analysis of ITH in lung adenocarcinomas from 147 patients using multi-region mass spectrometry of >5,000 regions, single-cell copy number sequencing of ~2,000 single cells, and cyclic immunofluorescence of >10 million cells. We identified two distinct spatial patterns among tumors, termed clustered and random geographic diversification (GD). These patterns were observed in the same samples using both proteomic and genomic data. The random proteomic GD pattern, which is characterized by decreased cell adhesion and lower levels of tumor-interacting endothelial cells, was significantly associated with increased risk of recurrence or death in two independent patient cohorts. Our study presents comprehensive spatial mapping of ITH in lung adenocarcinoma and provides insights into the mechanisms and clinical consequences of GD.

UnMICST: Deep learning with real augmentation for robust segmentation of highly multiplexed images of human tissues.

Upcoming technologies enable routine collection of highly multiplexed (20-60 channel), subcellular resolution images of mammalian tissues for research and diagnosis. Extracting single cell data from such images requires accurate image segmentation, a challenging problem commonly tackled with deep learning. In this paper, we report two findings that substantially improve image segmentation of tissues using a range of machine learning architectures. First, we unexpectedly find that the inclusion of intentionally defocused and saturated images in training data substantially improves subsequent image segmentation. Such real augmentation outperforms computational augmentation (Gaussian blurring). In addition, we find that it is practical to image the nuclear envelope in multiple tissues using an antibody cocktail thereby better identifying nuclear outlines and improving segmentation. The two approaches cumulatively and substantially improve segmentation on a wide range of tissue types. We speculate that the use of real augmentations will have applications in image processing outside of microscopy.

The Spatial Landscape of Progression and Immunoediting in Primary Melanoma at Single-Cell Resolution.

Cutaneous melanoma is a highly immunogenic malignancy that is surgically curable at early stages but life-threatening when metastatic. Here we integrate high-plex imaging, 3D high-resolution microscopy, and spatially resolved microregion transcriptomics to study immune evasion and immunoediting in primary melanoma. We find that recurrent cellular neighborhoods involving tumor, immune, and stromal cells change significantly along a progression axis involving precursor states, melanoma in situ, and invasive tumor. Hallmarks of immunosuppression are already detectable in precursor regions. When tumors become locally invasive, a consolidated and spatially restricted suppressive environment forms along the tumor-stromal boundary. This environment is established by cytokine gradients that promote expression of MHC-II and IDO1, and by PD1-PDL1-mediated cell contacts involving macrophages, dendritic cells, and T cells. A few millimeters away, cytotoxic T cells synapse with melanoma cells in fields of tumor regression. Thus, invasion and immunoediting can coexist within a few millimeters of each other in a single specimen. SIGNIFICANCE: The reorganization of the tumor ecosystem in primary melanoma is an excellent setting in which to study immunoediting and immune evasion. Guided by classic histopathology, spatial profiling of proteins and mRNA reveals recurrent morphologic and molecular features of tumor evolution that involve localized paracrine cytokine signaling and direct cell-cell contact. This article is highlighted in the In This Issue feature, p. 1397.

Stepwise-edited, human melanoma models reveal mutations' effect on tumor and microenvironment.

Establishing causal relationships between genetic alterations of human cancers and specific phenotypes of malignancy remains a challenge. We sequentially introduced mutations into healthy human melanocytes in up to five genes spanning six commonly disrupted melanoma pathways, forming nine genetically distinct cellular models of melanoma. We connected mutant melanocyte genotypes to malignant cell expression programs in vitro and in vivo, replicative immortality, malignancy, rapid tumor growth, pigmentation, metastasis, and histopathology. Mutations in malignant cells also affected tumor microenvironment composition and cell states. Our melanoma models shared genotype-associated expression programs with patient melanomas, and a deep learning model showed that these models partially recapitulated genotype-associated histopathological features as well. Thus, a progressive series of genome-edited human cancer models can causally connect genotypes carrying multiple mutations to phenotype.

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging.

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

Obesity Shapes Metabolism in the Tumor Microenvironment to Suppress Anti-Tumor Immunity.

Obesity is a major cancer risk factor, but how differences in systemic metabolism change the tumor microenvironment (TME) and impact anti-tumor immunity is not understood. Here, we demonstrate that high-fat diet (HFD)-induced obesity impairs CD8+ T cell function in the murine TME, accelerating tumor growth. We generate a single-cell resolution atlas of cellular metabolism in the TME, detailing how it changes with diet-induced obesity. We find that tumor and CD8+ T cells display distinct metabolic adaptations to obesity. Tumor cells increase fat uptake with HFD, whereas tumor-infiltrating CD8+ T cells do not. These differential adaptations lead to altered fatty acid partitioning in HFD tumors, impairing CD8+ T cell infiltration and function. Blocking metabolic reprogramming by tumor cells in obese mice improves anti-tumor immunity. Analysis of human cancers reveals similar transcriptional changes in CD8+ T cell markers, suggesting interventions that exploit metabolism to improve cancer immunotherapy.

Vascular Disease and Thrombosis in SARS-CoV-2-Infected Rhesus Macaques.

The COVID-19 pandemic has led to extensive morbidity and mortality throughout the world. Clinical features that drive SARS-CoV-2 pathogenesis in humans include inflammation and thrombosis, but the mechanistic details underlying these processes remain to be determined. In this study, we demonstrate endothelial disruption and vascular thrombosis in histopathologic sections of lungs from both humans and rhesus macaques infected with SARS-CoV-2. To define key molecular pathways associated with SARS-CoV-2 pathogenesis in macaques, we performed transcriptomic analyses of bronchoalveolar lavage and peripheral blood and proteomic analyses of serum. We observed macrophage infiltrates in lung and upregulation of macrophage, complement, platelet activation, thrombosis, and proinflammatory markers, including C-reactive protein, MX1, IL-6, IL-1, IL-8, TNFα, and NF-κB. These results suggest a model in which critical interactions between inflammatory and thrombosis pathways lead to SARS-CoV-2-induced vascular disease. Our findings suggest potential therapeutic targets for COVID-19.

HDAC6 mediates an aggresome-like mechanism for NLRP3 and pyrin inflammasome activation.

Inflammasomes are supramolecular complexes that play key roles in immune surveillance. This is accomplished by the activation of inflammatory caspases, which leads to the proteolytic maturation of interleukin 1ß (IL-1ß) and pyroptosis. Here, we show that nucleotide-binding domain, leucine-rich repeat, and pyrin domain-containing protein 3 (NLRP3)- and pyrin-mediated inflammasome assembly, caspase activation, and IL-1ß conversion occur at the microtubule-organizing center (MTOC). Furthermore, the dynein adapter histone deacetylase 6 (HDAC6) is indispensable for the microtubule transport and assembly of these inflammasomes both in vitro and in mice. Because HDAC6 can transport ubiquitinated pathological aggregates to the MTOC for aggresome formation and autophagosomal degradation, its role in NLRP3 and pyrin inflammasome activation also provides an inherent mechanism for the down-regulation of these inflammasomes by autophagy. This work suggests an unexpected parallel between the formation of physiological and pathological aggregates.

Author Correction: Immunogenomic profiling determines responses to combined PARP and PD-1 inhibition in ovarian cancer.

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

SARS-CoV-2 infection protects against rechallenge in rhesus macaques.

An understanding of protective immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is critical for vaccine and public health strategies aimed at ending the global coronavirus disease 2019 (COVID-19) pandemic. A key unanswered question is whether infection with SARS-CoV-2 results in protective immunity against reexposure. We developed a rhesus macaque model of SARS-CoV-2 infection and observed that macaques had high viral loads in the upper and lower respiratory tract, humoral and cellular immune responses, and pathologic evidence of viral pneumonia. After the initial viral clearance, animals were rechallenged with SARS-CoV-2 and showed 5 log10 reductions in median viral loads in bronchoalveolar lavage and nasal mucosa compared with after the primary infection. Anamnestic immune responses after rechallenge suggested that protection was mediated by immunologic control. These data show that SARS-CoV-2 infection induced protective immunity against reexposure in nonhuman primates.

Immunogenomic profiling determines responses to combined PARP and PD-1 inhibition in ovarian cancer.

Combined PARP and immune checkpoint inhibition has yielded encouraging results in ovarian cancer, but predictive biomarkers are lacking. We performed immunogenomic profiling and highly multiplexed single-cell imaging on tumor samples from patients enrolled in a Phase I/II trial of niraparib and pembrolizumab in ovarian cancer (NCT02657889). We identify two determinants of response; mutational signature 3 reflecting defective homologous recombination DNA repair, and positive immune score as a surrogate of interferon-primed exhausted CD8 + T-cells in the tumor microenvironment. Presence of one or both features associates with an improved outcome while concurrent absence yields no responses. Single-cell spatial analysis reveals prominent interactions of exhausted CD8 + T-cells and PD-L1 + macrophages and PD-L1 + tumor cells as mechanistic determinants of response. Furthermore, spatial analysis of two extreme responders shows differential clustering of exhausted CD8 + T-cells with PD-L1 + macrophages in the first, and exhausted CD8 + T-cells with cancer cells harboring genomic PD-L1 and PD-L2 amplification in the second.

Highly multiplexed immunofluorescence images and single-cell data of immune markers in tonsil and lung cancer.

In this data descriptor, we document a dataset of multiplexed immunofluorescence images and derived single-cell measurements of immune lineage and other markers in formaldehyde-fixed and paraffin-embedded (FFPE) human tonsil and lung cancer tissue. We used tissue cyclic immunofluorescence (t-CyCIF) to generate fluorescence images which we artifact corrected using the BaSiC tool, stitched and registered using the ASHLAR algorithm, and segmented using ilastik software and MATLAB. We extracted single-cell features from these images using HistoCAT software. The resulting dataset can be visualized using image browsers and analyzed using high-dimensional, single-cell methods. This dataset is a valuable resource for biological discovery of the immune system in normal and diseased states as well as for the development of multiplexed image analysis and viewing tools.

Qualifying antibodies for image-based immune profiling and multiplexed tissue imaging.

Multiplexed tissue imaging enables precise, spatially resolved enumeration and characterization of cell types and states in human resection specimens. A growing number of methods applicable to formalin-fixed, paraffin-embedded (FFPE) tissue sections have been described, the majority of which rely on antibodies for antigen detection and mapping. This protocol provides step-by-step procedures for confirming the selectivity and specificity of antibodies used in fluorescence-based tissue imaging and for the construction and validation of antibody panels. Although the protocol is implemented using tissue-based cyclic immunofluorescence (t-CyCIF) as an imaging platform, these antibody-testing methods are broadly applicable. We demonstrate assembly of a 16-antibody panel for enumerating and localizing T cells and B cells, macrophages, and cells expressing immune checkpoint regulators. The protocol is accessible to individuals with experience in microscopy and immunofluorescence; some experience in computation is required for data analysis. A typical 30-antibody dataset for 20 FFPE slides can be generated within 2 weeks.

Phase Transitions Drive the Formation of Vesicular Stomatitis Virus Replication Compartments.

RNA viruses that replicate in the cell cytoplasm typically concentrate their replication machinery within specialized compartments. This concentration favors enzymatic reactions and shields viral RNA from detection by cytosolic pattern recognition receptors. Nonsegmented negative-strand (NNS) RNA viruses, which include some of the most significant human, animal, and plant pathogens extant, form inclusions that are sites of RNA synthesis and are not circumscribed by a membrane. These inclusions share similarities with cellular protein/RNA structures such as P granules and nucleoli, which are phase-separated liquid compartments. Here we show that replication compartments of vesicular stomatitis virus (VSV) have the properties of liquid-like compartments that form by phase separation. Expression of the individual viral components of the replication machinery in cells demonstrates that the 3 viral proteins required for replication are sufficient to drive cytoplasmic phase separation. Therefore, liquid-liquid phase separation, previously linked to organization of P granules, nucleolus homeostasis, and cell signaling, plays a key role in host-pathogen interactions. This work suggests novel therapeutic approaches to the problem of combating NNS RNA viral infections.IMPORTANCE RNA viruses compartmentalize their replication machinery to evade detection by host pattern recognition receptors and concentrate the machinery of RNA synthesis. For positive-strand RNA viruses, RNA replication occurs in a virus-induced membrane-associated replication organelle. For NNS RNA viruses, the replication compartment is a cytoplasmic inclusion that is not circumscribed by a cellular membrane. Such structures were first observed in the cell bodies of neurons from humans infected with rabies virus and were termed Negri bodies. How the replication machinery that forms this inclusion remains associated in the absence of a membrane has been an enduring mystery. In this article, we present evidence that the VSV replication compartments form through phase separation. Phase separation is increasingly recognized as responsible for cellular structures as diverse as processing bodies (P-bodies) and nucleoli and was recently demonstrated for rabies virus. This article further links the fields of host-pathogen interaction with that of phase separation.