Moore, Raymond M.

Abstract Background Multiplexed immunofluorescence (MxIF) enables high-dimensional immune cell phenotyping and detailed characterization of the tumor immune microenvironment (TIME), but complex cyclic staining can introduce signal degradation and uncertainty in cell-level measurements. In routine practice, pathologists rely on hematoxylin and eosin (H&E) staining as the primary reference for tissue morphology and frequently cross-reference H&E to interpret MxIF findings. Methods We present a framework for aligning H&E and MxIF images at the tissue microarray (TMA) core level to support cross-modal analysis. Cell nucleus detections from each modality are used as anchor points, formulating the task as a point-set alignment problem. Coherent Point Drift (CPD) is applied for global alignment, followed by graph-matching–based refinement to improve local consistency. Evaluation on ovarian TMAs demonstrates consistent alignment performance for both restained and serial sections. We further explore the use of restained H&E as a reference for MxIF interpretation and investigate virtual H&E generation from MxIF data as a complementary approach when restained H&E is unavailable. Results The aligned images enable quantitative assessment of cross-modal feature concordance between MxIF-derived measurements and H&E morphology. The generated virtual H&E demonstrated similar cell population compared to the real H&E. Conclusions Cell-centric alignment can facilitate integrative multimodal histopathology analysis for TMA cores in spatial co-localization, cell feature validation, and virtual H&E generation.

Abstract The sentinel lymph node (SLN) is often the first site of melanoma spread and a critical immune checkpoint that determines local tumor dissemination to systemic metastasis. Although lymph nodes are central to antitumor immunity, the mechanisms that render them permissive to metastasis remain poorly understood. This study used spatial imaging to characterize immune and stromal phenotypic alteration in tumor-negative SLNs (SLN-) and tumor-positive SLNs (SLN+) from patients with cutaneous melanoma, compared to non-melanoma-associated normal lymph nodes (NLN). We performed multiplex immunofluorescence (MxIF) analysis of formalin-fixed, paraffin-embedded lymph node tissues from patients with cutaneous melanoma and control subjects, including two subjects per group (NLN, SLN-, and SLN+). Furthermore, quantitative single-cell segmentation and classification using QuPath generated median fluorescence intensities for immune, stromal, and tumor markers across LN, SLN- and SLN+ regions of interest. In addition to the initial analysis, we employed spatial neighborhood analysis to characterize markers expression profiles surrounding target immune population. This approach quantified local microenvironment based on marker normalized intensity-based co-expression pattern and stromal remodeling signatures. We grouped these markers based on the functional panels representing T-cell activation and exhaustion, myeloid polarization, tumor interaction, stromal remodeling, and apoptosis pathways providing insights into how prometastatic niche evolve within SLN. Normalized data, via centered log ratio (CLR) transformation and Gaussian Mixture Model (GMM) gating was analyzed using principal component analysis, UMAP/tSNE embedding, and hierarchical clustering in in python. Progressive immune remodeling was observed from normal NLN to SLN+, including expansion of regulatory T cells, enrichment of exhausted T cells as well as myeloid-derived suppressor cells (MDSC), and transition toward an M2-like macrophage phenotype. The spatial data from SLN- samples showed an early sign of immune suppression, with reduced cytotoxic capacity. These results show that SLN gets features of peripheral immune tolerance for melanoma early, prior to pathological dissection of metastatic spread supporting the existence of a pre-metastatic niche that facilitates lymphatic colonization of tumor cells. This work provides insights into early immune escape and may inform strategies for immunotherapy and early intervention in melanoma. Citation Format: Shankar Suman, Wendy K. Nevala, Noah A. Stueven, Raymond M. Moore, Chathu L. Atherton, Jeffrey E. Johnson, Ray Guo, James W. Jakub, Svetomir N. Markovic. Spatial remodeling of immune and stromal microenvironments in sentinel lymph nodes during melanoma progression [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 6193.

Abstract The highest risk for liver allograft rejection is within the first weeks to months after a liver transplant; however, while later-stage mechanisms are well described, the initial window at the moment of alloantigen influx (e.g., ~ 1 h post-reperfusion) remains unmapped. T-cell priming is initiated in lymph nodes, yet these sites are rarely profiled in humans and their relationship to rejection risk is unclear. To address this gap, we profiled paired pre- and ~ 1 h post-reperfusion hepatic draining lymph nodes from liver transplant recipients using multiplexed immunofluorescence and GeoMx whole-transcriptome profiling to identify features that predict rejection versus non-rejection. Early after reperfusion, CD4 + and CD8 + T cells became proximal to macrophages in rejection compared to dendritic cells in non-rejection, which also displayed an immunoregulatory-leaning profile marked by PD-1, IDO1, and CD39, suggesting that macrophages and dendritic cells differentially influence early T-cell alloimmunity. These results provide nodal targets for precise immune modulation.

Abstract Background: The tumor microenvironment (TME) plays a pivotal role in ovarian cancer (OC) progression and prognosis, while also contributing to biological differences between primary and metastatic tumors. This study examines OC TMEs to uncover associations with clinical prognosis and tumor biology by analyzing paired primary and metastatic tumors from the same patients. Methods: This study focused on 57 stage-IIIC/IV OC patients treated at the Mayo Clinic, all of whom underwent surgical debulking and adjuvant therapies. All patients provided informed consent, and the study received IRB approval. A tissue microarray (TMA) was specifically designed and constructed, containing paired primary and metastatic tissue cores from the same patients. 147 tissue cores from the constructed TMA were profiled using a 35-antibody panel, with 17 antibodies for cell phenotyping and 18 for functional state characterization. Profiling was performed using Cell DIVE, a multiplexed imaging solution at the single-cell level developed by Leica Microsystems. After cell segmentations and classifications, detailed spatial niche and regional analyses were performed at the single-cell level. Immune cell populations, as well as their TME distributions, were characterized to conduct prognostic associations and contrast primary and metastasis differences. Results: A total of 1,292,988 total cells (after excluding low-density signals and ambiguous classifications) were analyzed, with a median of 8,033 cells per core. Using spatial neighborhood analysis, 12 distinct spatial niches were identified with specific cell types and/or unique functional states. Among these, a macrophage- and fibroblast-enriched niche at the tumor-stroma interface was associated with improved overall survival (p=0.008). Another rare Treg-enriched niche, along with other immune cells, demonstrated protective effects, particularly in patients with residual macroscopic tumors after debulking surgery (p=0.05), while no significant effects were observed in those without visible residual disease (p=0.8). When comparing primary and metastatic tumors, metastatic lesions exhibited significantly increased tumor mass (p=2e-4) and a notable reduction in helper T cells (p=0.04) in tumor areas, as well as significantly increased fibroblast cells in tumor-stroma interfacing regions (p=0.002), reflecting distinct immunosuppressive and tumor-promoting shifts. Conclusion: This study reveals spatial and cellular differences in primary and metastatic OC TMEs, distribution differences linked to prognosis and surgical outcomes. These findings emphasize the potential of spatially resolved TME profiling to refine adjuvant and neoadjuvant therapeutic strategies. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 6266.

Abstract Purpose: Tumor-infiltrating lymphocytes are associated with a decreased risk of recurrence in human papillomavirus–associated oropharyngeal squamous cell carcinoma. The composition and spatial distribution of tumor-infiltrating lymphocytes and tumor-infiltrating immune cells are not well characterized. Experimental Design: Formalin-fixed, paraffin-embedded primary and lymph node (LN) tumor tissues from 10 progressors (cases) and 10 matched nonprogressors (controls) were interrogated by imaging mass cytometry. Immune, stromal, and tumor cells were quantified from selected regions of interest using machine learning. Nearest neighbors, cell–cell interactions, and niche analyses were performed. Results: In primary regions of interest, immune cell, lymphocyte, T cell, CD8+ T cell, and innate cell prevalence was significantly greater in controls. High prevalence of immune cells, lymphocytes, innate cells, and CD4+ T cells in primary tissues was significantly associated with increased time to event (TTE). Although primary and LN prevalence of T cells, CD4+ T cells, CD8+ T cells, macrophages, and tumor cells were significantly correlated, differences in LNs were neither significant nor associated with TTE. Average distances between T cells and the nearest B cells and between lymphocytes and the nearest tumor cells were decreased in control primary tissues. Interactions between B and T cells were less organized in primary tissues from cases. A niche predominantly comprising lymphocytes was associated with longer TTE. Conclusions: In human papillomavirus–associated oropharyngeal squamous cell carcinoma, immune cell subset prevalence in primary tumors differs with outcome and is associated with TTE. Interactions between B cells and both T-cell subsets are associated with longer TTE, underscoring the importance of active intratumoral immune responses in outcomes.

Abstract Motivation Multiplex immunofluorescence (MxIF) enables the quantification of multiple protein markers at a single-cell level while preserving spatial information, offering a powerful tool for studying tissue microenvironments. However, the flexibility in MxIF panel design poses challenges in standardizing cell phenotyping. Results We present SpaFlow, an efficient, customizable pipeline for unsupervised clustering and classification of MxIF data, implemented using Nextflow. SpaFlow performs quality control, clustering, and postclustering analysis on segmented and quantified MxIF data, facilitating reproducible and scalable analyses across various computing platforms. The SpaFlow pipeline integrates three clustering and classification packages—Seurat, SCIMAP, and CELESTA—each providing unique methodologies for identifying cell types based on phenotypic markers. A novel “meta-clustering” approach condenses clusters across multiple regions of interest into common meta-clusters, streamlining the cell-type identification process in large datasets. SpaFlow’s robust quality control steps, including signal summation and cell density filtering, mitigate artifacts that may impact clustering accuracy. We demonstrate the utility of SpaFlow in a case study involving 297 ovarian tumor cores, where SpaFlow successfully identified biologically meaningful cell populations, including tumor-infiltrating lymphocytes, efficiently and rapidly. Additionally, SpaFlow’s reproducibility is validated using serial tonsil sections, confirming its capability to consistently identify distinctive cell populations across matched ROIs. Availability and implementation SpaFlow is freely available with detailed documentation and examples at https://github.com/dimi-lab/SpaFlow.

Abstract Multiplexed Immunofluorescence (MxIF) enables detailed immune cell phenotyping, providing critical insights into cell behavior within the tumor immune microenvironment (TIME). However, signal integrity can be compromised due to the complex cyclic staining processes inherent to MxIF. Hematoxylin and Eosin (H&E) staining, on the other hand, offers complementary information through its depiction of cell morphology and texture patterns and is often visually cross-referenced with MxIF in clinical settings. In this study, we proposed a novel framework to align H&E and MxIF images for precise cross-modal cell feature validation. Using cell detection outputs from each modality as anchors, we formulated the multimodal image registration problem as point set alignment. Coherent Point Drift (CPD) is employed for initial alignment, followed by Graph Matching (GM) for refinement. Evaluations on ovarian cancer tissue microarrays (TMAs) demonstrate that our method achieves high alignment accuracy, enabling reliable validation of cell-level features across modalities for both restained and serial sections. Our results indicate that restained H&E enhances confidence in findings derived from MxIF. Additionally, we demonstrated the feasibility of generating high-quality virtual H&E images from MxIF data when restained H&E is unavailable, offering a viable alternative for integrated multimodal analysis.

CD36 expression in both immune and non-immune cells is known to be directly involved in cancer metastasis. Extracellular vesicles (EVs) secreted by malignant melanocytes play a vital role in developing tumor-promoting microenvironments, but it is unclear whether this is mediated through CD36. To understand the role of CD36 in melanoma, we first analyzed the SKCM dataset for clinical prognosis, evaluated the percentage of CD36 in lymphatic fluid-derived EVs (LEVs), and tested whether melanoma-derived EVs increase CD36 expression and induce M2-macrophage-like characteristics. Furthermore, we performed a multiplex immunofluorescence (MxIF) imaging analysis to evaluate the CD36 expression and its colocalization with various other cells in the lymph node (LN) of patients and control subjects. Our findings show that cutaneous melanoma patients have a worse clinical prognosis with high CD36 levels, and a higher percentage of CD36 in total LEVs were found at baseline in melanoma patients compared to control. We also found that monocytic and endothelial cells treated with melanoma EVs expressed more CD36 than untreated cells. Furthermore, melanoma-derived EVs can regulate immunosuppressive macrophage-like characteristics by upregulating CD36. The spatial imaging data show that cells in tumor-involved sentinel LNs exhibit a higher probability of CD36 expression than cells from control LNs, but this was not statistically significant. Conclusively, our findings demonstrated that CD36 plays a vital role in controlling the immunosuppressive microenvironment in the LN, which can promote the formation of a protumorigenic niche.

Abstract Background: Neuropilin 2 (NRP2) is a co-receptor that enhances signaling of growth factor and cytokine receptors and is expressed in tumor and microenvironment (TME) cells in melanoma and other tumors. Neuropilin 2 expression in tumor cells is associated with more aggressive tumors. Expression in tumor associated macrophages enhances efferocytosis and immunologically silent tumor cell death. Methods: Whole excisional lymph node biopsies were obtained from 10 treatment naïve patients with metastatic cutaneous melanoma who subsequently received anti-PD1 treatment. A single 5mm formalin-fixed paraffin embedded tissue section was used to assess a panel of 39 analytes by cyclic multiplex immunofluorescence (MxIF). Fields of view (~1mm2, n=288) were selected from pathologist-annotated regions of the tumor core, tumor-immune interface and adjacent lymphoid tissue. Pixel-based single cell segmentation and a supervised classifier approach was applied to resolve 12 distinct tumor, stromal and immune cell phenotypes and functional states (NRP2 positive or negative) in 1.5 million cells. Results: Tumor recurrence was observed in 6 patients (relapsed), whereas 4 patients were recurrence-free (non-relapsed) during follow up. Among the different cell types, NRP2 expression was more frequent in macrophages (44%), tumor cells (40%), dendritic cells (40%), followed by Tregs (31%) and T helper cells (29%), and less frequent in neutrophils (24%), cytotoxic T cells (23%) and B cells (15%). In colocalization permutation analysis, NRP2 expressing macrophages clustered away from tumor cells in non-relapsed patients while a random spatial distribution of NRP2 expressing macrophages with respect to tumor cells was observed in relapsed patients. In comparing the odds of tumor cell interactions with NRP2-expressing macrophages to those with NRP2 non-expressing macrophages, we found in the tumor:stromal interface (10% tumor cells) that cells from both non-relapsed and relapsed patients were more likely to interact with NRP2 non-expressing macrophages (OR = 0.73, 95% CI = [0.71, 0.75] and OR = 0.56, 95% CI = [0.55, 0.58] for non-relapsed and relapsed patients, respectively). However, in the tumor core (90% tumor cells), tumor cells from relapsed patients had significantly higher odds of interacting with NRP2 expressing macrophages, (OR = 1.15, 95% CI = [1.12, 1.18]), and this was not seen for non-relapsed patients (OR = 0.68, 95% CI = [0.66, 0.71]). In FOVs dominated by tumor cells (75% tumor), tumor cells from relapsed patients had 4.97 times the odds of expressing NRP2 compared to tumor cells from non-relapsed patients (95% CI = [2.23, 11.09], p = 0.0001). Conclusions: These spatial analytics suggest complex interactions among NRP2 expression in tumor and TME of melanoma determining risk for relapse with anti-PD1 treatment. Citation Format: Anastasios Dimou, Raymond M. Moore, Caitlin Ward, Alexey Leontovich, Ruifeng Guo, Chen Wang, Shankar Suman, Betty Dicke, Jill M. Schimke, Noah A. Stueven, Chathu A. Atherton, Wendy Nevala, Svetomir N. Markovic. Spatial analysis of neuropilin 2 expression in the microenvironment of melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2905.

Abstract Spread to regional lymph nodes is the most common site of metastasis from melanoma and almost all solid tumors. Lymphatic metastases begin in the first lymph nodes draining from the primary tumor, i.e., the sentinel lymph nodes (SLN), preceded by the development of the premetastatic niche (PMN) initiated by primary cutaneous melanoma. Melanoma mediators are known to play a critical role in developing the PMN through the immune subversion in SLN. The DCs are a vital component of LN immunity that directly relates to eliciting T-cell functions. Our multiplex immunofluorescence (MxIF) analysis of SLNs indicates significantly decreased levels of DCs and cytotoxic-T cells in both core and interface areas of tumors in recurred patients treated with immunotherapy compared to non-recurred melanoma patients. The goal of this study is to determine whether the lymphatic extracellular vesicles (L-EVs) of melanoma patients regulate DC tolerization in SLN in order to stimulate the development of PMN. To identify whether L-EVs stimulate DC suppression, we extracted L-EVs from the lymphatic fluid of patients with melanoma downstream of primary cutaneous melanomas, compared it with controls (non-malignant post-operative fluid collected from the lymph node dissection field, and evaluated their effect on DC maturation by challenging these L-EVs in monocytic (CD14+) and hematopoietic stem cells (CD34+) cells. We found that L-EVs derived from melanoma patients suppress differentiation of CD11c+HLA-DR+ DCs and their maturation (defined by CD83 and CD86) compared to control L-EVs in CD14+ based DC maturation assay. We also found a reduction of conventional type 1 DC (cDC1) in CD34+-based DC maturation assay those challenged with patients’ L-EVs compared to control. To understand the mechanism of DC suppression in patients' L-EVs, we also conducted a total RNA sequencing analysis of subsets of DCs from maturation experiments to evaluate differently regulated transcriptomic profiles. The gene ontology analysis indicates that patient-derived L-EVs interfere with tRNA processing, oncostatin M & IL-10 signaling, and other immunoregulatory pathways to alter the hematopoietic cell lineage pathways in the DC maturation. Furthermore, the proteomic analysis of L-EVs also indicated several differentially regulated cargo-identified proteins in patient L-EVs, which are associated with immunomodulation, including regulation of DC function compared to control L-EVs. In conclusion, our study revealed that the patient's L-EVs carry differential proteomic cargoes secreted by cutaneous melanoma, which may be involved in regulating DC differentiation and maturation to diminish their functions and support the developing PMN in SLN. Citation Format: Shankar Suman, Wendy K. Nevala, James W. Jakub, Ray Guo, Alexey A. Leontovich, Raymond M. Moore, Chathu L. Atherton, Liyi Geng, Jeffrey E. Johnson, Sarah A. McLaughlin, Svetomir N. Markovic. Melanoma-derived lymphatic extracellular vesicles suppress dendritic cell differentiation and maturation in the sentinel lymph node [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1503.