Single-cell transcriptomic analysis of the senescent microenvironment in bone metastasis

Abstract

Bone metastasis (BM) is a mortality-related event of late-stage cancer, with non-small cell lung cancer (NSCLC) being a common origin for BM. However, the detailed molecular profiling of the metastatic bone ecosystem is not fully understood, hindering the development of effective therapies for advanced patients. In this study, we examined the cellular heterogeneity between primary tumours and BM from tissues and peripheral blood by single-cell transcriptomic analysis, which was verified using multiplex immunofluorescence staining and public datasets. Our results demonstrate a senescent microenvironment in BM tissues of NSCLC. BM has a significantly higher infiltration of malignant cells with senescent characteristics relative to primary tumours, accompanied by aggravated metastatic properties. The endothelial-mesenchymal transition involved with SOX18 activation is related to the cellular senescence of vascular endothelial cells from BM. CD4Tstr cells, with pronounced stress and senescence states, are preferentially infiltrated in BM, indicating stress-related dysfunction contributing to the immunocompromised environment during tumour metastasis to bone. Moreover, we identify the SPP1 pathway-induced cellular crosstalk among T cells, vascular ECs and malignant cells in BM, which activates SOX18 and deteriorates patient survival. Our findings highlight the roles of cellular senescence in modulating the microenvironment of BM and implicate anti-senescence therapy for advanced NSCLC patients.

References

[1] SungH, FerlayJ, SiegelRL, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209-249.

[2] GoldstrawP, BallD, JettJR, et al. Non-small-cell lung cancer. Lancet. 2011;378(9804):1727-1740.

[3] SuvaLJ, WashamC, NicholasRW, Griffin RJ. Bone metastasis: mechanisms and therapeutic opportunities. Nat Rev Endocrinol. 2011;7(4):208-218.

[4] VicentS, Perurena N, GovindanR, LecandaF. Bone metastases in lung cancer. Potential novel approaches to therapy. Am J Respir Crit Care Med. 2015;192(7):799-809.

[5] AltorkiNK, Markowitz GJ, GaoD, et al. The lung microenvironment: an important regulator of tumour growth and metastasis. Nat Rev Cancer. 2019;19(1):9-31.

[6] WoodSL, Pernemalm M, CrosbiePA, WhettonAD. The role of the tumor-microenvironment in lung cancer-metastasis and its relationship to potential therapeutic targets. Cancer Treat Rev. 2014;40(4):558-566.

[7] ZhaoE, XuH, WangL, et al. Bone marrow and the control of immunity. Cell Mol Immunol. 2012;9(1):11-19.

[8] MukaidaN, ZhangD, SasakiSI. Emergence of cancer-associated fibroblasts as an indispensable cellular player in bone metastasis process. Cancers (Basel). 2020;12(10):2896.

[9] LiX, WangL, ChenS, et al. Adverse impact of bone metastases on clinical outcomes of patients with advanced non-small cell lung cancer treated with immune checkpoint inhibitors. Thorac Cancer. 2020;11(10):2812-2819.

[10] QiangH, LeiY, ShenY, et al. Pembrolizumab monotherapy or combination therapy for bone metastases in advanced non-small cell lung cancer: a real-world retrospective study. Transl Lung Cancer Res. 2022;11(1):87-99.

[11] Del ConteA, De Carlo E, BertoliE, et al. Bone metastasis and immune checkpoint inhibitors in non-small cell lung cancer (NSCLC): microenvironment and possible clinical implications. Int J Mol Sci. 2022;23(12):6832.

[12] WangC, YuQ, SongT, et al. The heterogeneous immune landscape between lung adenocarcinoma and squamous carcinoma revealed by single-cell RNA sequencing. Signal Transduct Target Ther. 2022;7(1):289.

[13] WangZ, LiZ, ZhouK, et al. Deciphering cell lineage specification of human lung adenocarcinoma with single-cell RNA sequencing. Nat Commun. 2021;12(1):6500.

[14] KimN, KimHK, LeeK, et al. Single-cell RNA sequencing demonstrates the molecular and cellular reprogramming of metastatic lung adenocarcinoma. Nat Commun. 2020;11(1):2285.

[15] MaynardA, McCoach CE, RotowJK, et al. Therapy-induced evolution of human lung cancer revealed by single-cell RNA sequencing. Cell. 2020;182(5):1232-1251.e22.

[16] WangZ, WangY, ChangM, et al. Single-cell transcriptomic analyses provide insights into the cellular origins and drivers of brain metastasis from lung adenocarcinoma. Neuro Oncol. 2023;25(7):1262-1274.

[17] SaulD, Kosinsky RL, AtkinsonEJ, et al. A new gene set identifies senescent cells and predicts senescence-associated pathways across tissues. Nat Commun. 2022;13(1):4827.

[18] FridmanAL, Tainsky MA. Critical pathways in cellular senescence and immortalization revealed by gene expression profiling. Oncogene. 2008;27(46):5975-5987.

[19] KangTW, YevsaT, WollerN, et al. Senescence surveillance of pre-malignant hepatocytes limits liver cancer development. Nature. 2011;479(7374):547-551.

[20] WuZ, UhlB, GiresO, Reichel CA. A transcriptomic pan-cancer signature for survival prognostication and prediction of immunotherapy response based on endothelial senescence. J Biomed Sci. 2023;30(1):21.

[21] HirzT, MeiS, SarkarH, et al. Dissecting the immune suppressive human prostate tumor microenvironment via integrated single-cell and spatial transcriptomic analyses. Nat Commun. 2023;14(1):663.

[22] KenswilKJG, Pisterzi P, Sánchez-DuffhuesG, et al. Endothelium-derived stromal cells contribute to hematopoietic bone marrow niche formation. Cell Stem Cell. 2021;28(4):653-670.e11.

[23] BenadjaoudMA, Soysouvanh F, TarletG, et al. Deciphering the dynamic molecular program of radiation-induced endothelial senescence. Int J Radiat Oncol Biol Phys. 2022;112(4):975-985.

[24] WangR, SongS, QinJ, et al. Evolution of immune and stromal cell states and ecotypes during gastric adenocarcinoma progression. Cancer Cell. 2023;41(8):1407-1426.e9.

[25] ChuY, DaiE, LiY, et al. Pan-cancer T cell atlas links a cellular stress response state to immunotherapy resistance. Nat Med. 2023;29(6):1550-1562.

[26] Cicin-SainL, Messaoudi I, ParkB, et al. Dramatic increase in naive T cell turnover is linked to loss of naive T cells from old primates. Proc Natl Acad Sci USA. 2007;104(50):19960-19965.

[27] LiX, LiC, ZhangW, Wang Y, QianP, HuangH. Inflammation and aging: signaling pathways and intervention therapies. Signal Transduct Target Ther. 2023;8(1):239.

[28] LiJ, LiangC, YangKY, et al. Specific ablation of CD4(+) T-cells promotes heart regeneration in juvenile mice. Theranostics. 2020;10(18):8018-8035.

[29] LiuX, HoftDF, PengG. Senescent T cells within suppressive tumor microenvironments: emerging target for tumor immunotherapy. J Clin Invest. 2020;130(3):1073-1083.

[30] QiuX, MaoQ, TangY, et al. Reversed graph embedding resolves complex single-cell trajectories. Nat Methods. 2017;14(10):979-982.

[31] De CeccoM, ItoT, PetrashenAP, et al. L1 drives IFN in senescent cells and promotes age-associated inflammation. Nature. 2019;566(7742):73-78.

[32] VictorelliS, Salmonowicz H, ChapmanJ, et al. Apoptotic stress causes mtDNA release during senescence and drives the SASP. Nature. 2023;622(7983):627-636.

[33] LiX, ZhaoS, BianX, et al. Signatures of EMT, immunosuppression, and inflammation in primary and recurrent human cutaneous squamous cell carcinoma at single-cell resolution. Theranostics. 2022;12(17):7532-7549.

[34] JinS, Guerrero-Juarez CF, ZhangL, et al. Inference and analysis of cell-cell communication using CellChat. Nat Commun. 2021;12(1):1088.

[35] ChengJ, WuH, XieC, et al. Single-cell mapping of large and small arteries during hypertensive aging. J Gerontol A Biol Sci Med Sci. 2024;79(2):glad188.

[36] ZhangW, BadoIL, HuJ, et al. The bone microenvironment invigorates metastatic seeds for further dissemination. Cell. 2021;184(9):2471-2486.e20.

[37] FornettiJ, WelmAL, StewartSA. Understanding the bone in cancer metastasis. J Bone Miner Res. 2018;33(12):2099-2113.

[38] Welch-ReardonKM, WuN, HughesCC. A role for partial endothelial-mesenchymal transitions in angiogenesis? Arterioscler Thromb Vasc Biol. 2015;35(2):303-308.

[39] LuoH, XiaX, HuangLB, et al. Pan-cancer single-cell analysis reveals the heterogeneity and plasticity of cancer-associated fibroblasts in the tumor microenvironment. Nat Commun. 2022;13(1):6619.

[40] OlbromskiM, Podhorska-Okołów M, DzięgielP. Role of SOX protein groups F and H in lung cancer progression. Cancers (Basel). 2020;12(11):3235.

[41] DonovanP, PatelJ, DightJ, et al. Endovascular progenitors infiltrate melanomas and differentiate towards a variety of vascular beds promoting tumor metastasis. Nat Commun. 2019;10(1):18.

[42] OvermanJ, Fontaine F, MoustaqilM, et al. Pharmacological targeting of the transcription factor SOX18 delays breast cancer in mice. Elife. 2017;6:e21221.

[43] FerraraR, Naigeon M, AuclinE, et al. Circulating T-cell Immunosenescence in patients with advanced non-small cell lung cancer treated with single-agent PD-1/PD-L1 inhibitors or platinum-based chemotherapy. Clin Cancer Res. 2021;27(2):492-503.

[44] LópezP, Rodríguez-Carrio J, Martínez-ZapicoA, Caminal-MonteroL, SuarezA. Senescent profile of angiogenic T cells from systemic lupus erythematosus patients. J Leukoc Biol. 2016;99(3):405-412.

[45] SeebauerCT, Brunner S, GlockzinG, et al. Peritoneal carcinomatosis of colorectal cancer is characterized by structural and functional reorganization of the tumor microenvironment inducing senescence and proliferation arrest in cancer cells. Onco Targets Ther. 2016;5(12):e1242543.

[46] TylutkaA, Morawin B, Wawrzyniak-GramackaE, et al. Immunosenescence in aging-related vascular dysfunction. Int J Mol Sci. 2022;23(21):13269.

[47] HwangHJ, LeeYR, KangD, et al. Endothelial cells under therapy-induced senescence secrete CXCL11, which increases aggressiveness of breast cancer cells. Cancer Lett. 2020;490:100-110.

[48] HastonS, Gonzalez-Gualda E, MorsliS, et al. Clearance of senescent macrophages ameliorates tumorigenesis in KRAS-driven lung cancer. Cancer Cell. 2023;41(7):1242-1260.e6.

[49] HaoC, CuiY, HuMU, et al. OPN-a splicing variant expression in non-small cell lung cancer and its effects on the bone metastatic abilities of lung cancer cells in vitro. Anticancer Res. 2017;37(5):2245-2254.

[50] YuS, ChenM, XuL, MaoE, SunS. A senescence-based prognostic gene signature for colorectal cancer and identification of the role of SPP1-positive macrophages in tumor senescence. Front Immunol. 2023;14:1175490.

[51] LambrechtsD, Wauters E, BoeckxB, et al. Phenotype molding of stromal cells in the lung tumor microenvironment. Nat Med. 2018;24(8):1277-1289.

[52] ZhuL, CaoZ, WangS, et al. Single-cell transcriptomics reveals peripheral immune responses in anti-synthetase syndrome-associated interstitial lung disease. Front Immunol. 2022;13:804034.

[53] WolfFA, Angerer P, TheisFJ. SCANPY: large-scale single-cell gene expression data analysis. Genome Biol. 2018;19(1):15.

[54] KorsunskyI, Millard N, FanJ, et al. Fast, sensitive and accurate integration of single-cell data with harmony. Nat Methods. 2019;16(12):1289-1296.

[55] ButlerA, Hoffman P, SmibertP, PapalexiE, SatijaR. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018;36(5):411-420.

[56] McGinnisCS, MurrowLM, GartnerZJ. DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 2019;8(4):329-337.e4.

[57] PatelAP, TiroshI, TrombettaJJ, et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science. 2014;344(6190):1396-1401.

[58] Van de SandeB, FlerinC, DavieK, et al. A scalable SCENIC workflow for single-cell gene regulatory network analysis. Nat Protoc. 2020;15(7):2247-2276.

[59] AndreattaM, Carmona SJ. UCell: Robust and scalable single-cell gene signature scoring. Comput Struct Biotechnol J. 2021;19:3796-3798.

[60] YuG, WangLG, HanY, HeQY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284-287.

[61] KumarL, Futschik ME. Mfuzz: a software package for soft clustering of microarray data. Bioinformation. 2007;2(1):5-7.

[62] LeekJT, Johnson WE, ParkerHS, JaffeAE, StoreyJD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882-883.

RIGHTS & PERMISSIONS

2024 The Author(s). Cell Proliferation published by Beijing Institute for Stem Cell and Regenerative Medicine and John Wiley & Sons Ltd.