Abstract:Introduction: A fundamental hallmark of cancer is that tumor cells repurpose the tissue microenvironment to promote their own survival. An increased understanding of these mechanisms may lead to improved microenvironment-directed therapies, particularly in lymphoid malignancies. In classic Hodgkin lymphoma (cHL), the rare malignant Hodgkin Reed Sternberg (HRS) cells are surrounded by a CD4+ T-cell and macrophage-rich inflammatory infiltrate. Recent multiplexed immunofluorescence studies suggest that the micron-scale niche around HRS cells is composed of distinct populations of PD-L1+ macrophages and CD4+ T cells, including regulatory CTLA4+ and LAG3+ subsets (Carey et al. Blood 2017, Patel et al. Blood 2019 and Aoki et al. Cancer Discov 2020). However, the topography of the intact tumor microenvironment of cHL requires further definition. Recent single-cell RNA sequencing studies have led to important insights into the biology of cHL; however, they do not adequately capture myeloid cells, fibroblasts, and HRS cells, likely due to the relative fragility of these cells in conventional tissue dissociation protocols. In this study, we use tandem single nucleus and spatially resolved RNA sequencing to systematically dissect the pro-tumorigenic cellular niche of cHL to define potentially targetable microenvironmental dependencies. Methods: Using single nucleus (10X Chromium V3.1 gene expression) and spatially resolved RNA sequencing (Slide-seqV2 and 10X Visium), we generated a dataset of 12 newly diagnosed cHL and 7 reactive lymphoid tissues comprising over 324,855 nuclei and 2.5 million spatially resolved gene expression profiles to construct a genome-scale single cell and spatially resolved atlas of cHL. To validate the expression of ligands/receptors on proximal cell types within the HRS cell niche at single-cell resolution, we also generated an orthogonal 1000-plex image-based spatial transcriptomic dataset (CosMx Spatial Molecular Imager). Additionally, we developed a novel Bayesian consensus tensor factorization approach (C-ZIPTF) to infer shared cancer-associated immune signatures and ligand-receptor interactions. Results: Single nucleus RNA sequencing captures a diversity of cell types and cell states in cHL, including distinct myeloid and stromal cell subsets and HRS cells. Copy number inference analysis identifies the HRS cells, and unsupervised consensus tensor factorization reveals recurrent HRS cell transcriptional states, including aberrant myeloid, stromal, neuronal, and secretory/cytokine programs. Of all clinical and histological parameters, EBV infection status showed the strongest association with HRS cell transcriptional state. Regulatory T cells, macrophages, and fibroblasts are quantitatively increased in cHL compared to reactive lymphoid tissues and additionally exhibit qualitatively altered expression programs. Using C-ZIPTF, we identify disease-defining and prognostically relevant immunosuppressive regulatory T cell and myeloid cell programs and an activated pro-fibrotic myofibroblast-like fibroblast signature. Spatial analysis provides strong evidence for a “niche” model, i.e., the spatial organization of cell types and states around HRS cells. Specifically, CD4+ T cells and myeloid cells are spatially enriched in the immediate proximity of HRS cells and, plasma cells are depleted. Moreover, the CD4+ T cells and myeloid cells exhibit a specialized pro-tumorigenic gene expression program with increased expression of multiple immunosuppressive cytokines, checkpoint molecules, and putative growth factors for HRS cells. Spatially-aware ligand-receptor interaction analysis enables systematic identification and prioritization of novel growth factors for HRS cells. Functional testing of these putative growth factors using an ex vivo 3D cell culture system confirms their growth-promoting effect on primary HRS cells. Conclusions: The HRS cells in cHL are surrounded by a distinctive pro-tumorigenic cellular niche that not only mediates immune evasion but may also directly promote tumor cell growth. These observations provide a rich set of hypotheses for preclinical model development and microenvironment-directed therapies. More broadly, our study demonstrates the power of spatial genomic approaches to deconvolute the molecular architecture of pro-tumorigenic cellular niches in tumor ecosystems.

Spatial transcriptomics exploration of the primary neuroblastoma microenvironment unveils novel paracrine interactions

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

Single-cell transcriptomics reveals shared immunosuppressive landscapes of mouse and human neuroblastoma

Medulloblastoma Spatial Transcriptomics Reveals Tumor Microenvironment Heterogeneity with High-Density Progenitor Cell Regions Correlating with High-Risk Disease

Single-cell analyses of metastatic bone marrow in human neuroblastoma reveals microenvironmental remodeling and metastatic signature

Single-cell RNA sequencing explored potential therapeutic targets by revealing the tumor microenvironment of neuroblastoma and its expression in cell death

Spatial Transcriptomic Sequencing of a DIPG-infiltrated Brainstem reveals Key Invasion Markers and Novel Ligand-Receptor Pairs contributing to Tumour to TME Crosstalk

Landscape of Bone Marrow Metastasis in Human Neuroblastoma Unraveled by Transcriptomics and Deep Multiplex Imaging

Single-cell RNA-seq reveals the transcriptional program underlying tumor progression and metastasis in neuroblastoma

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs

Neuroblastoma plasticity during metastatic progression stems from the dynamics of an early sympathetic transcriptomic trajectory

Abstract P5-06-01: Spatiotemporal gene expression analysis to identify novel tumor-neural crosstalk in TNBC brain metastases

Single-cell landscape analysis reveals distinct regression trajectories and novel prognostic biomarkers in primary neuroblastoma

Abstract 5466: In vivo characterization of neuroblastoma intratumoral heterogeneity

Genome-Scale High-Resolution Spatial Mapping of the Pro-Tumorigenic Cellular Niche in Classic Hodgkin Lymphoma

Abstract 868: A novel unsupervised machine learning model applied to neuroblastoma single-cell RNA-seq data reveals a drug-induced mesenchymal-like gene expression program

EPCO-04. SPATIAL TRANSCRIPTOMIC ANALYSIS OF MALIGNANT PERIPHERAL NERVE SHEATH TUMORS REVEALS THERAPEUTICALLY TARGETABLE MOLECULAR SIGNATURES IN REGIONS UNDERGOING HISTOPATHOLOGIC TRANSFORMATION

Angi-08. Spatial Transcriptomic Comparisons Of Glioblastoma Tissue Reveal Upregulation Of Neurodevelopmental Pathways And Synaptic Genes In Infiltrating Tumor Cells

Spatial transcriptomics reveals niche-specific enrichment and vulnerabilities of radial glial stem-like cells in malignant gliomas

Abstract 3877: Multi-modal analysis of paracrine signaling in group 3/4 medulloblastoma