Distinct malignant cell states and myeloid glutamate signaling associated with aggressive pancreatic neuroendocrine tumors

Abstract Pancreatic neuroendocrine tumors (PNET) are rare malignancies of the endocrine pancreas with diverse clinical outcomes. While some PNETs are indolent, others are aggressive and metastasize quickly. However, clinically-relevant molecular stratification for PNET to predict outcomes and guide therapeutic decision-making is limited. Thus, there is an urgent need to understand the molecular heterogeneity of PNETs to refine prognostication and discover novel therapeutic vulnerabilities. We performed single-nucleus RNA sequencing on resected primary and metastatic PNETs ( n = 20), including two PNETs with neoadjuvant treatment. We inferred gene expression programs (GEPs) of malignant and non-malignant cells and investigated associations with clinical outcomes. Next, we inferred interactions in the tumor microenvironment (TME) and performed transwell assays for functional validation. Finally, we explored genomic and transcriptomic evolution in a unique case study of an untreated primary PNET with two asynchronous hepatic metastases. A malignant GEP enriched for neural/synaptic signaling genes was associated with worse overall survival, broad chromosomal loss of heterozygosity, and alternative lengthening of telomeres. Another malignant GEP enriched for VEGF signaling increased throughout metastatic progression in our case study. We found that macrophage-derived glutamate drives polarization towards an immunosuppressive phenotype and activates the MAPK/ERK pathway in malignant cells to increase migratory capacity. This study provides a detailed single-nucleus transcriptomic classification of malignant, stromal, and immune cell types and states in PNETs, their interactions in the TME, and associations with clinical outcomes. The refined molecular taxonomy of PNET may guide the development of more efficacious biomarkers and therapeutic strategies.