Advanced high-grade serous ovarian cancer (HGSC) is an aggressive disease that accounts for 70% of all ovarian cancer deaths. Nevertheless, 15% of patients diagnosed with advanced HGSC survive more than 10 years. The elucidation of predictive markers of these long-term survivors (LTS) could help identify therapeutic targets for the disease, and thus improve patient survival rates. To investigate the stromal heterogeneity of the tumor microenvironment (TME) in ovarian cancer, we used spatial transcriptomics to generate spatially resolved transcript profiles in treatment-naïve advanced HGSC from LTS and short-term survivors (STS) and determined the association between cancer-associated fibroblasts (CAF) heterogeneity and survival in patients with advanced HGSC. Spatial transcriptomics and single-cell RNA-sequencing data were integrated to distinguish tumor and stroma regions, and a computational method was developed to investigate spatially resolved ligand–receptor interactions between various tumor and CAF subtypes in the TME. A specific subtype of CAFs and its spatial location relative to a particular ovarian cancer cell subtype in the TME correlated with long-term survival in patients with advanced HGSC. Also, increased APOE-LRP5 cross-talk occurred at the stroma-tumor interface in tumor tissues from STS compared with LTS. These findings were validated using multiplex IHC. Overall, this spatial transcriptomics analysis revealed spatially resolved CAF-tumor cross-talk signaling networks in the ovarian TME that are associated with long-term survival of patients with HGSC. Further studies to confirm whether such cross-talk plays a role in modulating the malignant phenotype of HGSC and could serve as a predictive biomarker of patient survival are warranted. Significance: Generation of spatially resolved gene expression patterns in tumors from patients with ovarian cancer surviving more than 10 years allows the identification of novel predictive biomarkers and therapeutic targets for better patient management. See related commentary by Kelliher and Lengyel, p. 1383

Perineural invasion (PNI) is a well-established factor of poor prognosis in multiple cancer types¹, yet its mechanism remains unclear. Here we provide clinical and mechanistic insights into the role of PNI and cancer-induced nerve injury (CINI) in resistance to anti-PD-1 therapy. Our study demonstrates that PNI and CINI of tumour-associated nerves are associated with poor response to anti-PD-1 therapy among patients with cutaneous squamous cell carcinoma, melanoma and gastric cancer. Electron microscopy and electrical conduction analyses reveal that cancer cells degrade the nerve fibre myelin sheets. The injured neurons respond by autonomously initiating IL-6- and type I interferon-mediated inflammation to promote nerve healing and regeneration. As the tumour grows, the CINI burden increases, and its associated inflammation becomes chronic and skews the general immune tone within the tumour microenvironment into a suppressive and exhaustive state. The CINI-driven anti-PD-1 resistance can be reversed by targeting multiple steps in the CINI signalling process: denervating the tumour, conditional knockout of the transcription factor mediating the injury signal within neurons (Atf3), knockout of interferon-α receptor signalling (Ifnar1^−/−) or by combining anti-PD-1 and anti-IL-6-receptor blockade. Our findings demonstrate the direct immunoregulatory roles of CINI and its therapeutic potential.

Intraductal papillary mucinous neoplasms (IPMN) of the pancreas are bona fide precursor lesions of pancreatic ductal adenocarcinoma (PDAC). The most common subtype of IPMNs harbors a gastric foveolar-type epithelium, and these low-grade mucinous neoplasms are harbingers of IPMNs with high-grade dysplasia and cancer. The molecular underpinning of gastric differentiation in IPMNs is unknown, although identifying drivers of this indolent phenotype might enable opportunities for intercepting progression to high-grade IPMN and cancer. We conducted spatial transcriptomics on a cohort of IPMNs, followed by orthogonal and cross-species validation studies, which established the transcription factor NKX6-2 as a key determinant of gastric cell identity in low-grade IPMNs. Loss of NKX6-2 expression is a consistent feature of IPMN progression, while reexpression of Nkx6-2 in murine IPMN lines recapitulates the aforementioned gastric transcriptional program and glandular morphology. Our study identifies NKX6-2 as a previously unknown transcription factor driving indolent gastric differentiation in IPMN pathogenesis. Significance: Identification of the molecular features driving IPMN development and differentiation is critical to prevent cancer progression and enhance risk stratification. We used spatial profiling to characterize the epithelium and microenvironment of IPMN, which revealed a previously unknown link between NKX6-2 and gastric differentiation, the latter associated with indolent biological potential. See related commentary by Ben-Shmuel and Scherz-Shouval, p. 1768. This article is highlighted in the In This Issue feature, p. 1749

Inactivating TP53 mutations leads to a loss of function of p53, but can also often result in oncogenic gain-of-function (GOF) of mutant p53 (mutp53) proteins which promotes tumor development and progression. The GOF activities of TP53 mutations are well documented, but the mechanisms involved remain poorly understood. Here, we study the mutp53 interactome and find that by targeting minichromosome maintenance complex components (MCMs), GOF mutp53 predisposes cells to replication stress and chromosomal instability (CIN), leading to a tumor cell-autonomous and cyclic GMP–AMP synthase (cGAS)-stimulator of interferon genes (STING)-dependent cytosolic DNA response that activates downstream non-canonical nuclear factor kappa light chain enhancer of activated B cell (NC-NF-κB) signaling. Consequently, GOF mutp53-MCMs-CIN-cytosolic DNA-cGAS-STING-NC-NF-κB signaling promotes tumor cell metastasis and an immunosuppressive tumor microenvironment through antagonizing interferon signaling and regulating genes associated with pro-tumorigenic inflammation. Our findings have important implications for understanding not only the GOF activities of TP53 mutations but also the genome-guardian role of p53 and its inactivation during tumor development and progression.

Most platforms used for the molecular reconstruction of the tumor–immune microenvironment (TIME) of a solid tumor fail to explore the spatial context of the three-dimensional (3D) space of the tumor at a single-cell resolution, and thus lack information about cell–cell or cell–extracellular matrix (ECM) interactions. To address this issue, a pipeline which integrated multiplex spatially resolved multi-omics platforms was developed to identify crosstalk signaling networks among various cell types and the ECM in the 3D TIME of two FFPE (formalin-fixed paraffin embedded) gynecologic tumor samples. These platforms include non-targeted mass spectrometry imaging (glycans, metabolites, and peptides) and Stereo-seq (spatial transcriptomics) and targeted seqIF (IHC proteomics). The spatially resolved imaging data in a two- and three-dimensional space demonstrated various cellular neighborhoods in both samples. The collection of spatially resolved analytes in a voxel (3D pixel) across serial sections of the tissue was also demonstrated. Data collected from this analytical pipeline were used to construct spatial 3D maps with single-cell resolution, which revealed cell identity, activation, and energized status. These maps will provide not only insights into the molecular basis of spatial cell heterogeneity in the TIME, but also novel predictive biomarkers and therapeutic targets, which can improve patient survival rates.

Comprehensive investigation of CD8+ T cells in acute myeloid leukemia (AML) is essential for developing immunotherapeutic strategies beyond immune checkpoint blockade. Herein, we performed single-cell RNA profiling of CD8+ T cells from 3 healthy bone marrow donors and 23 newly diagnosed (NewlyDx) and 8 relapsed/refractory (RelRef) patients with AML. Cells coexpressing canonical exhaustion markers formed a cluster constituting <1% of all CD8+ T cells. We identified two effector CD8+ T-cell subsets characterized by distinct cytokine and metabolic profiles that were differentially enriched in NewlyDx and RelRef patients. We refined a 25-gene CD8-derived signature correlating with therapy resistance, including genes associated with activation, chemoresistance, and terminal differentiation. Pseudotemporal trajectory analysis supported enrichment of a terminally differentiated state in CD8+ T cells with high CD8-derived signature expression at relapse or refractory disease. Higher expression of the 25-gene CD8 AML signature correlated with poorer outcomes in previously untreated patients with AML, suggesting that the bona fide state of CD8+ T cells and their degree of differentiation are clinically relevant. Immune clonotype tracking revealed more phenotypic transitions in CD8 clonotypes in NewlyDx than in RelRef patients. Furthermore, CD8+ T cells from RelRef patients had a higher degree of clonal hyperexpansion associated with terminal differentiation and higher CD8-derived signature expression. Clonotype-derived antigen prediction revealed that most previously unreported clonotypes were patient-specific, suggesting significant heterogeneity in AML immunogenicity. Thus, immunologic reconstitution in AML is likely to be most successful at earlier disease stages when CD8+ T cells are less differentiated and have greater capacity for clonotype transitions.