Single-cell and spatial transcriptomics reveals an anti-tumor neutrophil subgroup in microwave thermochemotherapy-treated lip cancer

Abstract

Microwave thermochemotherapy (MTC) has been applied to treat lip squamous cell carcinoma (LSCC), but a deeper understanding of its therapeutic mechanisms and molecular biology is needed. To address this, we used single-cell transcriptomics (scRNA-seq) and spatial transcriptomics (ST) to highlight the pivotal role of tumor-associated neutrophils (TANs) among tumor-infiltrating immune cells and their therapeutic response to MTC. MNDA⁺ TANs with anti-tumor activity (N1-phenotype) are found to be abundantly infiltrated by MTC with benefit of increased blood perfusion, and these TANs are characterized by enhanced cytotoxicity, ameliorated hypoxia, and upregulated IL1B, activating T&NK cells and fibroblasts via IL1B-IL1R. In this highly anti-tumor immunogenic and hypoxia-reversed microenvironment under MTC, fibroblasts accumulated in the tumor front (TF) can recruit N1-TANs via CXCL2-CXCR2 and clear N2-TANs (pro-tumor phenotype) via CXCL12-CXCR4, which results in the aggregation of N1-TANs and extracellular matrix (ECM) deposition. In addition, we construct an N1-TANs marker, MX2, which positively correlates with better prognosis in LSCC patients, and employ deep learning techniques to predict expression of MX2 from hematoxylin-eosin (H&E)-stained images so as to conveniently guide decision making in clinical practice. Collectively, our findings demonstrate that the N1-TANs/fibroblasts defense wall formed in response to MTC effectively combat LSCC.

References

[1] Miranda-FilhoA, BrayF. Global patterns and trends in cancers of the lip, tongue and mouth. Oral. Oncol., 2020, 102: 104551.

[2] LallasA, et al. . Clinical and dermatoscopic predictors of squamous cell carcinoma of the lips: a case-control, multicentric study. J. Eur. Acad. Dermatol. Venereol., 2022, 36: 222-227.

[3] TydingsC, SharmaKV, KimA, YarmolenkoPS. Emerging hyperthermia applications for pediatric oncology. Adv. Drug Deliv. Rev., 2020, 163-164: 157-167.

[4] IsselsRD, et al. . Effect of neoadjuvant chemotherapy plus regional hyperthermia on long-term outcomes among patients with localized high-risk soft tissue sarcoma: the EORTC 62961-ESHO 95 randomized clinical trial. JAMA Oncol., 2018, 4: 483-492.

[5] WessalowskiR, et al. . Regional deep hyperthermia for salvage treatment of children and adolescents with refractory or recurrent non-testicular malignant germ-cell tumours: an open-label, non-randomised, single-institution, phase 2 study. Lancet Oncol., 2013, 14: 843-852.

[6] DattaNR, PuricE, KlingbielD, GomezS, BodisS. Hyperthermia and radiation therapy in locoregional recurrent breast cancers: a systematic review and meta-analysis. Int J. Radiat. Oncol. Biol. Phys., 2016, 94: 1073-1087.

[7] YangZ, et al. . Thermal immuno-nanomedicine in cancer. Nat. Rev. Clin. Oncol., 2023, 20: 116-134.

[8] GuoD, et al. . Exosomes from heat-stressed tumour cells inhibit tumour growth by converting regulatory T cells to Th17 cells via IL-6. Immunology, 2018, 154: 132-143.

[9] Sorin, M. et al. Single-cell spatial landscapes of the lung tumour immune microenvironment. Nature614, 548–554 (2023).

[10] JohnsonDE, et al. . Head and neck squamous cell carcinoma. Nat. Rev. Dis. Prim., 2020, 6: 92.

[11] ChenY, McAndrewsKM, KalluriR. Clinical and therapeutic relevance of cancer-associated fibroblasts. Nat. Rev. Clin. Oncol., 2021, 18: 792-804.

[12] MoncadaR, et al. . Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat. Biotechnol., 2020, 38: 333-342.

[13] ZhuY, et al. . Rational design of biomaterials to potentiate cancer thermal therapy. Chem. Rev., 2023, 11: 7326-7378.

[14] JaillonS, et al. . Neutrophil diversity and plasticity in tumour progression and therapy. Nat. Rev. Cancer, 2020, 20: 485-503.

[15] Veglia, F. et al. Analysis of classical neutrophils and polymorphonuclear myeloid-derived suppressor cells in cancer patients and tumor-bearing mice. J. Exp. Med. 218, e20201803 (2021).

[16] GungabeesoonJ, et al. . A neutrophil response linked to tumor control in immunotherapy. Cell, 2023, 186: 1448-1464.e1420.

[17] GabrilovichDI, Ostrand-RosenbergS, BronteV. Coordinated regulation of myeloid cells by tumours. Nat. Rev. Immunol., 2012, 12: 253-268.

[18] PapayannopoulosV. Neutrophil extracellular traps in immunity and disease. Nat. Rev. Immunol., 2018, 18: 134-147.

[19] Surendran, V., Rutledge, D., Colmon, R. & Chandrasekaran, A. A novel tumor-immune microenvironment (TIME)-on-Chip mimics three dimensional neutrophil-tumor dynamics and neutrophil extracellular traps (NETs)-mediated collective tumor invasion. Biofabrication13, 10.1088/1758-5090/abe1cf (2021).

[20] Howard, A. G. et al. MobileNets: Efficient convolutional neural networks for mobile vision applications. arXiv. https://doi.org/10.48550/arXiv.1704.04861 (2017).

[21] Dosovitskiy, A., Beyer, L., Kolesnikov, A., Weissenborn, D. & Houlsby, N. An image is worth 16 x 16 words: transformers for image recognition at scale. arXiv. https://doi.org/10.48550/arXiv.2010.11929 (2020).

[22] HuJ, ShenL, AlbanieS, SunG, WuE. Squeeze-and-excitation networks. IEEE Trans. Pattern Anal. Mach. Intell., 2020, 42: 2011-2023.

[23] MantovaniA, CassatellaMA, CostantiniC, JaillonS. Neutrophils in the activation and regulation of innate and adaptive immunity. Nat. Rev. Immunol., 2011, 11: 519-531.

[24] HedrickCC, MalanchiI. Neutrophils in cancer: heterogeneous and multifaceted. Nat. Rev. Immunol., 2022, 22: 173-187.

[25] ShaulME, FridlenderZG. Tumour-associated neutrophils in patients with cancer. Nat. Rev. Clin. Oncol., 2019, 16: 601-620.

[26] GaldieroMR, et al. . Occurrence and significance of tumor-associated neutrophils in patients with colorectal cancer. Int J. Cancer, 2016, 139: 446-456.

[27] ZhangH, et al. . Tumor-infiltrating neutrophils is prognostic and predictive for postoperative adjuvant chemotherapy benefit in patients with castric cancer. Ann. Surg., 2018, 267: 311-318.

[28] PosabellaA, et al. . High density of CD66b in primary high-grade ovarian cancer independently predicts response to chemotherapy. J. Cancer Res Clin. Oncol., 2020, 146: 127-136.

[29] MahiddineK, et al. . Relief of tumor hypoxia unleashes the tumoricidal potential of neutrophils. J. Clin. Invest., 2020, 130: 389-403.

[30] GehD, et al. . Neutrophils as potential therapeutic targets in hepatocellular carcinoma. Nat. Rev. Gastroenterol. Hepatol., 2022, 19: 257-273.

[31] Fotouhi-ArdakaniN, et al. . Role for myeloid nuclear differentiation antigen in the regulation of neutrophil apoptosis during sepsis. Am. J. Respir. Crit. Care Med., 2010, 182: 341-350.

[32] MilotE, Fotouhi-ArdakaniN, FilepJG. Myeloid nuclear differentiation antigen, neutrophil apoptosis and sepsis. Front Immunol., 2012, 3: 397.

[33] KlimenkovaO, et al. . A lack of secretory leukocyte protease inhibitor (SLPI) causes defects in granulocytic differentiation. Blood, 2014, 123: 1239-1249.

[34] BroderickL, HoffmanHM. IL-1 and autoinflammatory disease: biology, pathogenesis and therapeutic targeting. Nat. Rev. Rheumatol., 2022, 18: 448-463.

[35] SimsJE, SmithDE. The IL-1 family: regulators of immunity. Nat. Rev. Immunol., 2010, 10: 89-102.

[36] PonzettaA, et al. . Neutrophils driving unconventional T cells mediate resistance against murine sarcomas and selected human tumors. Cell, 2019, 178: 346-360.e324.

[37] ChangCJ, et al. . Targeting tumor-infiltrating Ly6G(+) myeloid cells improves sorafenib efficacy in mouse orthotopic hepatocellular carcinoma. Int J. Cancer, 2018, 142: 1878-1889.

[38] TakeshimaT, et al. . Key role for neutrophils in radiation-induced antitumor immune responses: potentiation with G-CSF. Proc. Natl Acad. Sci. USA, 2016, 113: 11300-11305.

[39] SongM, et al. . Cancer-associated fibroblast-mediated cellular crosstalk supports hepatocellular carcinoma progression. Hepatology, 2021, 73: 1717-1735.

[40] Albrengues, J. et al. Neutrophil extracellular traps produced during inflammation awaken dormant cancer cells in mice. Science361, eaao4227 (2018).

[41] MeherAK, et al. . Novel role of IL (Interleukin)-1β in neutrophil extracellular trap formation and abdominal aortic aneurysms. Arterioscler Thromb. Vasc. Biol., 2018, 38: 843-853.

[42] XiaoY, et al. . Cathepsin C promotes breast cancer lung metastasis by modulating neutrophil infiltration and neutrophil extracellular trap formation. Cancer Cell, 2021, 39: 423-437.e427.

[43] WigerbladG, KaplanMJ. Neutrophil extracellular traps in systemic autoimmune and autoinflammatory diseases. Nat. Rev. Immunol., 2023, 23: 274-288.

[44] ZhangJ, et al. . Hypoxia-inducible factor-1α/interleukin-1β signaling enhances hepatoma epithelial-mesenchymal transition through macrophages in a hypoxic-inflammatory microenvironment. Hepatology, 2018, 67: 1872-1889.

[45] WangYX, et al. . Interferon-inducible MX2 is a host restriction factor of hepatitis B virus replication. J. Hepatol., 2020, 72: 865-876.

[46] MengXW, et al. . MX2: Identification and systematic mechanistic analysis of a novel immune-related biomarker for systemic lupus erythematosus. Front Immunol., 2022, 13: 978851.

[47] HaoY, et al. . Integrated analysis of multimodal single-cell data. Cell, 2021, 184: 3573-3587.e3529.

[48] HuC, et al. . CellMarker 2.0: an updated database of manually curated cell markers in human/mouse and web tools based on scRNA-seq data. Nucleic Acids Res., 2023, 51: D870-d876.

[49] BergenstråhleJ, LarssonL, LundebergJ. Seamless integration of image and molecular analysis for spatial transcriptomics workflows. BMC Genomics, 2020, 21. 482

[50] LaubacherME, AdesSE. The Rcs phosphorelay is a cell envelope stress response activated by peptidoglycan stress and contributes to intrinsic antibiotic resistance. J. Bacteriol., 2008, 190: 2065-2074.

[51] RitchieME, et al. . limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res., 2015, 43: e47.

[52] GuZ, EilsR, SchlesnerM. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics, 2016, 32: 2847-2849.

[53] Wu, T. et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation2, 100141 (2021).

[54] ColapricoA, et al. . TCGAbiolinks: an R/bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res., 2016, 44: e71.

[55] HänzelmannS, CasteloR, GuinneyJ. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinforma., 2013, 14. 7

[56] JinS, et al. . Inference and analysis of cell-cell communication using CellChat. Nat. Commun., 2021, 12. 1088

[57] BrowaeysR, SaelensW, SaeysY. NicheNet: modeling intercellular communication by linking ligands to target genes. Nat. Methods, 2020, 17: 159-162.

[58] LangfelderP, HorvathS. WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma., 2008, 9. 559

Funding

National Natural Science Foundation of China (National Science Foundation of China)(82073000)

National Science Foundation of Sichuan Province (No. 2022YFS0289), and Interdisciplinary Innovation Project of West China College of Stomatology, Sichuan University(RD-03-202004).

RIGHTS & PERMISSIONS

The Author(s)