A Hitchhiker's guide to high-dimensional tissue imaging with multiplexed ion beam imaging.

Advancements in multiplexed tissue imaging technologies are vital in shaping our understanding of tissue microenvironmental influences in disease contexts. These technologies now allow us to relate the phenotype of individual cells to their higher-order roles in tissue organization and function. Multiplexed Ion Beam Imaging (MIBI) is one of such technologies, which uses metal isotope-labeled antibodies and secondary ion mass spectrometry (SIMS) to image more than 40 protein markers simultaneously within a single tissue section. Here, we describe an optimized MIBI workflow for high-plex analysis of Formalin-Fixed Paraffin-Embedded (FFPE) tissues following antigen retrieval, metal isotope-conjugated antibody staining, imaging using the MIBI instrument, and subsequent data processing and analysis. While this workflow is focused on imaging human FFPE samples using the MIBI, this workflow can be easily extended to model systems, biological questions, and multiplexed imaging modalities.

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.

Epstein-Barr Virus Orchestrates Spatial Reorganization and Immunomodulation within the Classic Hodgkin Lymphoma Tumor Microenvironment.

Classic Hodgkin Lymphoma (cHL) is a tumor composed of rare malignant Hodgkin and Reed-Sternberg (HRS) cells nested within a T-cell rich inflammatory immune infiltrate. cHL is associated with Epstein-Barr Virus (EBV) in 25% of cases. The specific contributions of EBV to the pathogenesis of cHL remain largely unknown, in part due to technical barriers in dissecting the tumor microenvironment (TME) in high detail. Herein, we applied multiplexed ion beam imaging (MIBI) spatial pro-teomics on 6 EBV-positive and 14 EBV-negative cHL samples. We identify key TME features that distinguish between EBV-positive and EBV-negative cHL, including the relative predominance of memory CD8 T cells and increased T-cell dysfunction as a function of spatial proximity to HRS cells. Building upon a larger multi-institutional cohort of 22 EBV-positive and 24 EBV-negative cHL samples, we orthogonally validated our findings through a spatial multi-omics approach, coupling whole transcriptome capture with antibody-defined cell types for tu-mor and T-cell populations within the cHL TME. We delineate contrasting transcriptomic immunological signatures between EBV-positive and EBV-negative cases that differently impact HRS cell proliferation, tumor-immune interactions, and mecha-nisms of T-cell dysregulation and dysfunction. Our multi-modal framework enabled a comprehensive dissection of EBV-linked reorganization and immune evasion within the cHL TME, and highlighted the need to elucidate the cellular and molecular fac-tors of virus-associated tumors, with potential for targeted therapeutic strategies.

Association of Baseline Tumor-Specific Neoantigens and CD8+ T-Cell Infiltration With Immune-Related Adverse Events Secondary to Immune Checkpoint Inhibitors.

PURPOSE: Recent evidence has shown that higher tumor mutational burden strongly correlates with an increased risk of immune-related adverse events (irAEs). By using an integrated multiomics approach, we further studied the association between relevant tumor immune microenvironment (TIME) features and irAEs. METHODS: Leveraging the US Food and Drug Administration Adverse Event Reporting System, we extracted cases of suspected irAEs to calculate the reporting odds ratios (RORs) of irAEs for cancers treated with immune checkpoint inhibitors (ICIs). TIME features for 32 cancer types were calculated on the basis of the cancer genomic atlas cohorts and indirectly correlated with each cancer's ROR for irAEs. A separate ICI-treated cohort of non-small-cell lung cancer (NSCLC) was used to evaluate the correlation between tissue-based immune markers (CD8+, PD-1/L1+, FOXP3+, tumor-infiltrating lymphocytes [TILs]) and irAE occurrence. RESULTS: The analysis of 32 cancers and 33 TIME features demonstrated a significant association between irAE RORs and the median number of base insertions and deletions (INDEL), neoantigens (r = 0.72), single-nucleotide variant neoantigens (r = 0.67), and CD8+ T-cell fraction (r = 0.51). A bivariate model using the median number of INDEL neoantigens and CD8 T-cell fraction had the highest accuracy in predicting RORs (adjusted r2 = 0.52, P = .002). Immunoprofile assessment of 156 patients with NSCLC revealed a strong trend for higher baseline median CD8+ T cells within patients' tumors who experienced any grade irAEs. Using machine learning, an expanded ICI-treated NSCLC cohort (n = 378) further showed a treatment duration-independent association of an increased proportion of high TIL (>median) in patients with irAEs (59.7% v 44%, P = .005). This was confirmed by using the Fine-Gray competing risk approach, demonstrating higher baseline TIL density (>median) associated with a higher cumulative incidence of irAEs (P = .028). CONCLUSION: Our findings highlight a potential role for TIME features, specifically INDEL neoantigens and baseline-immune infiltration, in enabling optimal irAE risk stratification of patients.

Graph Fourier transform for spatial omics representation and analyses of complex organs.

Spatial omics technologies are capable of deciphering detailed components of complex organs or tissue in cellular and subcellular resolution. A robust, interpretable, and unbiased representation method for spatial omics is necessary to illuminate novel investigations into biological functions, whereas a mathematical theory deficiency still exists. We present SpaGFT (Spatial Graph Fourier Transform), which provides a unique analytical feature representation of spatial omics data and elucidates molecular signatures linked to critical biological processes within tissues and cells. It outperformed existing tools in spatially variable gene prediction and gene expression imputation across human/mouse Visium data. Integrating SpaGFT representation into existing machine learning frameworks can enhance up to 40% accuracy of spatial domain identification, cell type annotation, cell-to-spot alignment, and subcellular hallmark inference. SpaGFT identified immunological regions for B cell maturation in human lymph node Visium data, characterized secondary follicle variations from in-house human tonsil CODEX data, and detected extremely rare subcellular organelles such as Cajal body and Set1/COMPASS. This new method lays the groundwork for a new theoretical model in explainable AI, advancing our understanding of tissue organization and function.

Biomarkers of pembrolizumab efficacy in advanced anal squamous cell carcinoma: analysis of a phase II clinical trial and a cohort of long-term responders.

BACKGROUND: Recent trials suggest that programmed cell death 1 (PD-1)-directed immunotherapy may be beneficial for some patients with anal squamous cell carcinoma and biomarkers predictive of response are greatly needed. METHODS: This multicenter phase II clinical trial (NCT02919969) enrolled patients with metastatic or locally advanced incurable anal squamous cell carcinoma (n=32). Patients received pembrolizumab 200 mg every 3 weeks. The primary endpoint of the trial was objective response rate (ORR). Exploratory objectives included analysis of potential predictive biomarkers including assessment of tumor-associated immune cell populations with multichannel immunofluorescence and analysis of circulating tumor tissue modified viral-human papillomavirus DNA (TTMV-HPV DNA) using serially collected blood samples. To characterize the clinical features of long-term responders, we combined data from our prospective trial with a retrospective cohort of patients with anal cancer treated with anti-PD-1 immunotherapy (n=18). RESULTS: In the phase II study, the ORR to pembrolizumab monotherapy was 9.4% and the median progression-free survival was 2.2 months. Despite the high level of HPV positivity observed with circulating TTMV-HPV DNA testing, the majority of patients had low levels of tumor-associated CD8+PD-1+ T cells on pretreatment biopsy. Patients who benefited from pembrolizumab had decreasing TTMV-HPV DNA scores and a complete responder's TTMV-HPV DNA became undetectable. Long-term pembrolizumab responses were observed in one patient from the trial (5.3 years) and three patients (2.5, 6, and 8 years) from the retrospective cohort. Long-term responders had HPV-positive tumors, lacked liver metastases, and achieved a radiological complete response. CONCLUSIONS: Pembrolizumab has durable efficacy in a rare subset of anal cancers. However, despite persistence of HPV infection, indicated by circulating HPV DNA, most advanced anal cancers have low numbers of tumor-associated CD8+PD-1+ T cells and are resistant to pembrolizumab.

MAPS: pathologist-level cell type annotation from tissue images through machine learning.

Highly multiplexed protein imaging is emerging as a potent technique for analyzing protein distribution within cells and tissues in their native context. However, existing cell annotation methods utilizing high-plex spatial proteomics data are resource intensive and necessitate iterative expert input, thereby constraining their scalability and practicality for extensive datasets. We introduce MAPS (Machine learning for Analysis of Proteomics in Spatial biology), a machine learning approach facilitating rapid and precise cell type identification with human-level accuracy from spatial proteomics data. Validated on multiple in-house and publicly available MIBI and CODEX datasets, MAPS outperforms current annotation techniques in terms of speed and accuracy, achieving pathologist-level precision even for typically challenging cell types, including tumor cells of immune origin. By democratizing rapidly deployable and scalable machine learning annotation, MAPS holds significant potential to expedite advances in tissue biology and disease comprehension.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

Genomic and Immunophenotypic Landscape of Acquired Resistance to PD-(L)1 Blockade in Non-Small-Cell Lung Cancer.

PURPOSE: Although immune checkpoint inhibitors (ICI) have extended survival in patients with non-small-cell lung cancer (NSCLC), acquired resistance (AR) to ICI frequently develops after an initial benefit. However, the mechanisms of AR to ICI in NSCLC are largely unknown. METHODS: Comprehensive tumor genomic profiling, machine learning-based assessment of tumor-infiltrating lymphocytes, multiplexed immunofluorescence, and/or HLA-I immunohistochemistry (IHC) were performed on matched pre- and post-ICI tumor biopsies from patients with NSCLC treated with ICI at the Dana-Farber Cancer Institute who developed AR to ICI. Two additional cohorts of patients with intervening chemotherapy or targeted therapies between biopsies were included as controls. RESULTS: We performed comprehensive genomic profiling and immunophenotypic characterization on samples from 82 patients with NSCLC and matched pre- and post-ICI biopsies and compared findings with a control cohort of patients with non-ICI intervening therapies between biopsies (chemotherapy, N = 32; targeted therapies, N = 89; both, N = 17). Putative resistance mutations were identified in 27.8% of immunotherapy-treated cases and included acquired loss-of-function mutations in STK11, B2M, APC, MTOR, KEAP1, and JAK1/2; these acquired alterations were not observed in the control groups. Immunophenotyping of matched pre- and post-ICI samples demonstrated significant decreases in intratumoral lymphocytes, CD3e+ and CD8a+ T cells, and PD-L1-PD1 engagement, as well as increased distance between tumor cells and CD8+PD-1+ T cells. There was a significant decrease in HLA class I expression in the immunotherapy cohort at the time of AR compared with the chemotherapy (P = .005) and the targeted therapy (P = .01) cohorts. CONCLUSION: These findings highlight the genomic and immunophenotypic heterogeneity of ICI resistance in NSCLC, which will need to be considered when developing novel therapeutic strategies aimed at overcoming resistance.

Systematic benchmarking of imaging spatial transcriptomics platforms in FFPE tissues.

Emerging imaging spatial transcriptomics (iST) platforms and coupled analytical methods can recover cell-to-cell interactions, groups of spatially covarying genes, and gene signatures associated with pathological features, and are thus particularly well-suited for applications in formalin fixed paraffin embedded (FFPE) tissues. Here, we benchmarked the performance of three commercial iST platforms on serial sections from tissue microarrays (TMAs) containing 23 tumor and normal tissue types for both relative technical and biological performance. On matched genes, we found that 10x Xenium shows higher transcript counts per gene without sacrificing specificity, but that all three platforms concord to orthogonal RNA-seq datasets and can perform spatially resolved cell typing, albeit with different false discovery rates, cell segmentation error frequencies, and with varying degrees of sub-clustering for downstream biological analyses. Taken together, our analyses provide a comprehensive benchmark to guide the choice of iST method as researchers design studies with precious samples in this rapidly evolving field.