Small cell lung cancer (SCLC) is an aggressive malignancy that includes subtypes defined by differential expression of ASCL1, NEUROD1, and POU2F3 (SCLC-A, -N, and -P, respectively). To define the heterogeneity of tumors and their associated microenvironments across subtypes, we sequenced 155,098 transcriptomes from 21 human biospecimens, including 54,523 SCLC transcriptomes. We observe greater tumor diversity in SCLC than lung adenocarcinoma, driven by canonical, intermediate, and admixed subtypes. We discover a PLCG2-high SCLC phenotype with stem-like, pro-metastatic features that recurs across subtypes and predicts worse overall survival. SCLC exhibits greater immune sequestration and less immune infiltration than lung adenocarcinoma, and SCLC-N shows less immune infiltrate and greater T-cell dysfunction than SCLC-A. We identify a profibrotic, immunosuppressive monocyte/macrophage population in SCLC tumors that is particularly associated with the recurrent, PLCG2-high subpopulation.

The response to tumor-initiating inflammatory and genetic insults can vary among morphologically indistinguishable cells, suggesting as yet uncharacterized roles for epigenetic plasticity during early neoplasia. To investigate the origins and impact of such plasticity, we performed single-cell analyses on normal, inflamed, premalignant, and malignant tissues in autochthonous models of pancreatic cancer. We reproducibly identified heterogeneous cell states that are primed for diverse, late-emerging neoplastic fates and linked these to chromatin remodeling at cell-cell communication loci. Using an inference approach, we revealed signaling gene modules and tissue-level cross-talk, including a neoplasia-driving feedback loop between discrete epithelial and immune cell populations that was functionally validated in mice. Our results uncover a neoplasia-specific tissue-remodeling program that may be exploited for pancreatic cancer interception.

Radiotherapy is a mainstay cancer therapy whose antitumor effects partially depend on T cell responses. However, the role of Natural Killer (NK) cells in radiotherapy remains unclear. Here, using a reverse translational approach, we show a central role of NK cells in the radiation-induced immune response involving a CXCL8/IL-8–dependent mechanism. In a randomized controlled pancreatic cancer trial, CXCL8 increased under radiotherapy, and NK cell positively correlated with prolonged overall survival. Accordingly, NK cells preferentially infiltrated irradiated pancreatic tumors and exhibited CD56 ^dim -like cytotoxic transcriptomic states. In experimental models, NF-κB and mTOR orchestrated radiation-induced CXCL8 secretion from tumor cells with senescence features causing directional migration of CD56 ^dim NK cells, thus linking senescence-associated CXCL8 release to innate immune surveillance of human tumors. Moreover, combined high-dose radiotherapy and adoptive NK cell transfer improved tumor control over monotherapies in xenografted mice, suggesting NK cells combined with radiotherapy as a rational cancer treatment strategy.

Factor analysis decomposes single-cell gene expression data into a minimal set of gene programs that correspond to processes executed by cells in a sample. However, matrix factorization methods are prone to technical artifacts and poor factor interpretability. We address these concerns with Spectra, an algorithm that combines user-provided gene programs with the detection of novel programs that together best explain expression covariation. Spectra incorporates existing gene sets and cell-type labels as prior biological information, explicitly models cell type and represents input gene sets as a gene–gene knowledge graph using a penalty function to guide factorization toward the input graph. We show that Spectra outperforms existing approaches in challenging tumor immune contexts, as it finds factors that change under immune checkpoint therapy, disentangles the highly correlated features of CD8⁺ T cell tumor reactivity and exhaustion, finds a program that explains continuous macrophage state changes under therapy and identifies cell-type-specific immune metabolic programs.

Locoregional failure (LRF) in patients with breast cancer post-surgery and post-irradiation is linked to a dismal prognosis. In a refined new model, we identified ectonucleotide pyrophosphatase/phosphodiesterase 1/CD203a (ENPP1) to be closely associated with LRF. ENPP1hi circulating tumor cells (CTC) contribute to relapse by a self-seeding mechanism. This process requires the infiltration of polymorphonuclear myeloid-derived suppressor cells and neutrophil extracellular trap (NET) formation. Genetic and pharmacologic ENPP1 inhibition or NET blockade extends relapse-free survival. Furthermore, in combination with fractionated irradiation, ENPP1 abrogation obliterates LRF. Mechanistically, ENPP1-generated adenosinergic metabolites enhance haptoglobin (HP) expression. This inflammatory mediator elicits myeloid invasiveness and promotes NET formation. Accordingly, a significant increase in ENPP1 and NET formation is detected in relapsed human breast cancer tumors. Moreover, high ENPP1 or HP levels are associated with poor prognosis. These findings unveil the ENPP1/HP axis as an unanticipated mechanism exploited by tumor cells linking inflammation to immune remodeling favoring local relapse. Significance: CTC exploit the ENPP1/HP axis to promote local recurrence post-surgery and post-irradiation by subduing myeloid suppressor cells in breast tumors. Blocking this axis impairs tumor engraftment, impedes immunosuppression, and obliterates NET formation, unveiling new opportunities for therapeutic intervention to eradicate local relapse and ameliorate patient survival. This article is highlighted in the In This Issue feature, p. 1171