Single-cell sequencing reveals a senescent immune landscape in bone marrow lesions inducing articular cartilage damage in osteoarthritis

Abstract

Bone marrow lesions (BML) are early signs of osteoarthritis (OA) and are strongly correlated with the deterioration of cartilage lesions. Single-cell RNA sequencing (scRNA-seq) analyses were performed on BM from non-BML and BML areas and articular cartilage from intact and damaged areas to explore BML landscape and BML-cartilage crosstalk. We revealed the immune landscape of BM in non-BML and BML, and the transition to pro-inflammatory states of clusters in BMLs, such as classical monocytes and non-classical monocytes. Non-classical monocytes have high inflammation, OA gene signatures, and senescence scores, and are potential primary clusters promoting OA progression. Histological signs of OA related to the cellular landscape in damaged cartilage were identified, including PreFC exhaustion. The BM-cartilage crosstalk at the cell-cell interaction (CCIs) level and the TNF signal transmitted by non-classical monocytes are the critical CCIs in BML-induced cartilage damage, and PreFC is one of the primary receivers of the signal. We further validated the higher senescence level of non-classical monocyte and FC-2 in OA mice, compared with classical monocyte and PreFC, respectively. Transcription factor 7 like 2 (TCF7L2) was identified as a shared transcription factor in the senescence of monocytes and chondrocytes, facilitating the development of the senescence-associated secretory phenotype (SASP). Therefore, senescent non-classical monocytes promote BMLs and inflammation and senescence of chondrocytes by modulating BML–cartilage crosstalk in OA, with TCF7L2 serving as a regulator.

References

[1] Martel-Pelletier J, et al.. Osteoarthritis. Nat. Rev. Dis. Prim., 2016, 2 16072

[2] Hansen RT. Bone marrow lesions: plugging the holes in our knowledge using animal models. Nat. Rev. Rheumatol, 2023, 19: 429-445

[3] Driban JB, et al.. Evaluation of bone marrow lesion volume as a knee osteoarthritis biomarker-longitudinal relationships with pain and structural changes: data from the Osteoarthritis Initiative. Arthritis Res. Ther., 2013, 15: R112

[4] Raynauld J-P, et al.. Correlation between bone lesion changes and cartilage volume loss in patients with osteoarthritis of the knee as assessed by quantitative magnetic resonance imaging over a 24-month period. Ann. Rheum. Dis., 2008, 67: 683-688

[5] Davies-Tuck ML, et al.. Development of bone marrow lesions is associated with adverse effects on knee cartilage while resolution is associated with improvement-a potential target for prevention of knee osteoarthritis: a longitudinal study. Arthritis Res. Ther., 2010, 12: R10

[6] Dore D, et al.. Bone marrow lesions predict site-specific cartilage defect development and volume loss: a prospective study in older adults. Arthritis Res. Ther., 2010, 12: R222

[7] Kuttapitiya A, et al.. Microarray analysis of bone marrow lesions in osteoarthritis demonstrates upregulation of genes implicatedin osteochondral turnover, neurogenesis and inflammation. Ann. Rheum. Dis, 2017, 76: 1764-1773

[8] Bowes MA, et al.. Osteoarthritic bone marrow lesions almost exclusively colocate with denuded cartilage: a 3D study using data from the Osteoarthritis Initiative. Ann. Rheum. Dis., 2016, 75: 1852-1857

[9] Moradi K, et al.. Progression of bone marrow lesions and the development of knee osteoarthritis: osteoarthritis initiative data. Radiology, 2024, 312 e240470

[10] Findlay DM, Kuliwaba JS. Bone-cartilage crosstalk: a conversation for understanding osteoarthritis. Bone Res., 2016, 4: 16028

[11] Faust HJ, et al.. IL-17 and immunologically induced senescence regulate response to injury in osteoarthritis. J. Clin. Invest., 2020, 130: 5493-5507

[12] Smiljanovic B, et al.. Monocyte alterations in rheumatoid arthritis are dominated by preterm release from bone marrow and prominent triggering in the joint. Ann. Rheum. Dis., 2018, 77: 300-308

[13] Koushesh S, et al.. The osteoarthritis bone score (OABS): a new histological scoring system for the characterisation of bone marrow lesions in osteoarthritis. Osteoarthr. Cartil., 2022, 30: 746-755

[14] Hsueh M-F, Zhang X, Wellman SS, Bolognesi MP, Kraus VB. Synergistic roles of macrophages and neutrophils in osteoarthritis progression. Arthritis Rheumatol. Hoboken NJ, 2021, 73: 89-99

[15] Li Y, Luo W, Zhu S, Lei G. T Cells in Osteoarthritis: alterations and beyond. Front. Immunol., 2017, 8: 356

[16] Alahdal M, et al.. Potential efficacy of dendritic cell immunomodulation in the treatment of osteoarthritis. Rheumatol. Oxf. Engl., 2021, 60: 507-517

[17] Haubruck P, Pinto MM, Moradi B, Little CB, Gentek R. Monocytes, macrophages, and their potential niches in synovial joints – therapeutic targets in post-traumatic osteoarthritis?. Front. Immunol., 2021, 12: 763702

[18] Xie X, Doody GM, Shuweihdi F, Conaghan PG, Ponchel F. B-cell capacity for expansion and differentiation into plasma cells are altered in osteoarthritis. Osteoarthr. Cartil., 2023, 31: 1176-1188

[19] Salazar-Noratto GE, De Nijs N, Stevens HY, Gibson G, Guldberg RE. Regional gene expression analysis of multiple tissues in an experimental animal model of post-traumatic osteoarthritis. Osteoarthr. Cartil., 2019, 27: 294-303

[20] Fisch KM, et al.. Identification of transcription factors responsible for dysregulated networks in human osteoarthritis cartilage by global gene expression analysis. Osteoarthr. Cartil., 2018, 26: 1531-1538

[21] Kuttapitiya A, et al.. Microarray analysis of bone marrow lesions in osteoarthritis demonstrates upregulation of genes implicated in osteochondral turnover, neurogenesis and inflammation. Ann. Rheum. Dis., 2017, 76: 1764-1773

[22] Rai MF, et al.. Single cell omics for musculoskeletal research. Curr. Osteoporos. Rep., 2021, 19: 131-140

[23] Chen Y, et al.. A high-resolution route map reveals distinct stages of chondrocyte dedifferentiation for cartilage regeneration. Bone Res., 2022, 10: 38

[24] Swahn H, et al.. Senescent cell population with ZEB1 transcription factor as its main regulator promotes osteoarthritis in cartilage and meniscus. Ann. Rheum. Dis., 2023, 82: 403-415

[25] Nanus DE, et al.. Synovial tissue from sites of joint pain in knee osteoarthritis patients exhibits a differential phenotype with distinct fibroblast subsets. EBioMedicine, 2021, 72 103618

[26] Tong L, et al.. Current understanding of osteoarthritis pathogenesis and relevant new approaches. Bone Res., 2022, 10: 60

[27] Wang T, et al.. Single-cell RNA sequencing in orthopedic research. Bone Res., 2023, 11: 10

[28] Sebastian A, et al.. Single-cell RNA-Seq reveals changes in immune landscape in post-traumatic osteoarthritis. Front. Immunol., 2022, 13: 938075

[29] Lin P, et al.. Advancing skeletal health and disease research with single-cell RNA sequencing. Mil. Med. Res., 2024, 11: 33

[30] Domínguez Conde C, et al.. Cross-tissue immune cell analysis reveals tissue-specific features in humans. Science, 2022, 376: eabl5197

[31] Hu Y, et al.. Single-cell RNA-sequencing analysis reveals the molecular mechanism of subchondral bone cell heterogeneity in the development of osteoarthritis. RMD Open, 2022, 8: e002314

[32] Hao Y, et al.. Integrated analysis of multimodal single-cell data. Cell, 2021, 184: 3573-3587.e29

[33] Zhang L, et al.. Lineage tracking reveals dynamic relationships of T cells in colorectal cancer. Nature, 2018, 564: 268-272

[34] Zhang L, et al.. Single-cell analyses inform mechanisms of myeloid-targeted therapies in colon cancer. Cell, 2020, 181: 442-459.e29

[35] Bruschi M, et al.. Annexin A1 and autoimmunity: from basic science to clinical applications. Int. J. Mol. Sci., 2018, 19: 1348

[36] Kang I, Bucala R. The immunobiology of MIF: function, genetics and prospects for precision medicine. Nat. Rev. Rheumatol., 2019, 15: 427-437

[37] Toegel S, et al.. Galectin-1 couples glycobiology to inflammation in osteoarthritis through the activation of an NF-κB-regulated gene network. J. Immunol. Baltim. Md. 1950, 2016, 196: 1910-1921

[38] Raghu H, et al.. CCL2/CCR2, but not CCL5/CCR5, mediates monocyte recruitment, inflammation and cartilage destruction in osteoarthritis. Ann. Rheum. Dis., 2017, 76: 914-922

[39] Ji Q, et al.. Single-cell RNA-seq analysis reveals the progression of human osteoarthritis. Ann. Rheum. Dis., 2019, 78: 100-110

[40] Wang X, et al.. Comparison of the major cell populations among osteoarthritis, Kashin-Beck disease and healthy chondrocytes by single-cell RNA-seq analysis. Cell Death Dis., 2021, 12 551

[41] Pandey A, Bhutani N. Profiling joint tissues at single-cell resolution: advances and insights. Nat. Rev. Rheumatol., 2024, 20: 7-20

[42] Ong S-M, et al.. The pro-inflammatory phenotype of the human non-classical monocyte subset is attributed to senescence. Cell Death Dis., 2018, 9: 1-12

[43] Qu Y, et al.. A comprehensive analysis of single-cell RNA transcriptome reveals unique SPP1+ chondrocytes in human osteoarthritis. Comput. Biol. Med., 2023, 160: 106926

[44] Fan Y, et al.. Unveiling inflammatory and prehypertrophic cell populations as key contributors to knee cartilage degeneration in osteoarthritis using multi-omics data integration. Ann. Rheum. Dis., 2024, 83: 926-944

[45] Qiu X, et al.. Reversed graph embedding resolves complex single-cell trajectories. Nat. Methods, 2017, 14: 979-982

[46] Bergen V, Lange M, Peidli S, Wolf FA, Theis FJ. Generalizing RNA velocity to transient cell states through dynamical modeling. Nat. Biotechnol., 2020, 38: 1408-1414

[47] Gerwin N, et al.. Angiopoietin-like 3-derivative LNA043 for cartilage regeneration in osteoarthritis: a randomized phase 1 trial. Nat. Med., 2022, 28: 2633-2645

[48] Fu W, et al.. Nav1.7 as a chondrocyte regulator and therapeutic target for osteoarthritis. Nature, 2024, 625: 557-565

[49] Zhao Y-P, et al.. Progranulin protects against osteoarthritis through interacting with TNF-α and β-Catenin signalling. Ann. Rheum. Dis., 2015, 74: 2244-2253

[50] Yang S, et al.. NAMPT (visfatin), a direct target of hypoxia-inducible factor-2α, is an essential catabolic regulator of osteoarthritis. Ann. Rheum. Dis., 2015, 74: 595-602

[51] Xie Y, Zinkle A, Chen L, Mohammadi M. Fibroblast growth factor signalling in osteoarthritis and cartilage repair. Nat. Rev. Rheumatol., 2020, 16: 547-564

[52] Wang C, Shen J, Ying J, Xiao D, O’Keefe RJ. FoxO1 is a crucial mediator of TGF-β/TAK1 signaling and protects against osteoarthritis by maintaining articular cartilage homeostasis. Proc. Natl. Acad. Sci. USA., 2020, 117: 30488-30497

[53] Jaswal AP, et al.. BMP signaling: a significant player and therapeutic target for osteoarthritis. Osteoarthr. Cartil., 2023, 31: 1454-1468

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

[55] Rowe MA, et al.. Reduced osteoarthritis severity in aged mice with deletion of macrophage migration inhibitory factor. Arthritis Rheumatol. Hoboken NJ, 2017, 69: 352-361

[56] Arikawa T, et al.. Galectin-9 accelerates transforming growth factor beta3-induced differentiation of human mesenchymal stem cells to chondrocytes. Bone, 2009, 44: 849-857

[57] Du X, Cai L, Xie J, Zhou X. The role of TGF-beta3 in cartilage development and osteoarthritis. Bone Res., 2023, 11: 2

[58] Li D-P, et al.. LILRB2/PirB mediates macrophage recruitment in fibrogenesis of nonalcoholic steatohepatitis. Nat. Commun., 2023, 14 4436

[59] Weinmann D, et al.. Galectin-8 induces functional disease markers in human osteoarthritis and cooperates with galectins-1 and -3. Cell. Mol. Life Sci. CMLS, 2018, 75: 4187-4205

[60] Hu W, Chen Y, Dou C, Dong S. Microenvironment in subchondral bone: predominant regulator for the treatment of osteoarthritis. Ann. Rheum. Dis., 2021, 80: 413-422

[61] Saul D, et al.. A new gene set identifies senescent cells and predicts senescence-associated pathways across tissues. Nat. Commun., 2022, 13 4827

[62] Van de Sande B, et al.. A scalable SCENIC workflow for single-cell gene regulatory network analysis. Nat. Protoc., 2020, 15: 2247-2276

[63] Wang J. TFTF: An R-based integrative tool for decoding human transcription factor–target interactions. Biomolecules, 2024, 14: 749

[64] Rossi M, et al.. Pleiotropic effects of BAFF on the senescence-associated secretome and growth arrest. eLife, 2023, 12: e84238

[65] Kirsch V, Ramge J-M, Schoppa A, Ignatius A, Riegger J. In vitro characterization of doxorubicin-mediated stress-induced premature senescence in human chondrocytes. Cells, 2022, 11: 1106

[66] Zeng M, et al.. Comprehensive analysis on subchondral bone marrow lesions of human osteoarthritis by integrating bulk and single-cell transcriptomes. BMC Musculoskelet. Disord., 2023, 24 677

[67] Misharin AV, et al.. Non-classical Ly6C⁻ monocytes drive the development of inflammatory arthritis in mice. Cell Rep., 2014, 9: 591-604

[68] Wojdasiewicz P, et al.. The chemokine CX3CL1 (Fractalkine) and its receptor CX3CR1: occurrence and potential role in osteoarthritis. Arch. Immunol. Ther. Exp., 2014, 62: 395-403

[69] Canalis E, Schilling L, Denker E. TNFα has differential effects on the transcriptome profile of selected populations in murine cartilage. Osteoarthr. Cartil. Open, 2024, 6: 100528

[70] Puchner A, et al.. Non-classical monocytes as mediators of tissue destruction in arthritis. Ann. Rheum. Dis., 2018, 77: 1490-1497

[71] López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. Hallmarks of aging: an expanding universe. Cell, 2023, 186: 243-278

[72] Coryell PR, Diekman BO, Loeser RF. Mechanisms and therapeutic implications of cellular senescence in osteoarthritis. Nat. Rev. Rheumatol., 2021, 17: 47-57

[73] Fafián-Labora JA, O’Loghlen A. Classical and nonclassical intercellular communication in senescence and ageing. Trends Cell Biol., 2020, 30: 628-639

[74] Coppé J-P, Desprez P-Y, Krtolica A, Campisi J. The senescence-associated secretory phenotype: the dark side of tumor suppression. Annu. Rev. Pathol., 2010, 5: 99-118

[75] Han Z, Wang K, Ding S, Zhang M. Cross-talk of inflammation and cellular senescence: a new insight into the occurrence and progression of osteoarthritis. Bone Res., 2024, 12: 69

[76] Xie J, et al.. Cellular senescence in knee osteoarthritis: molecular mechanisms and therapeutic implication. Ageing Res. Rev., 2021, 70: 101413

[77] Wang H, et al.. Osteoclasts and osteoarthritis: novel intervention targets and therapeutic potentials during aging. Aging Cell, 2024, 23 e14092

[78] Rockel JS, et al.. Hedgehog inhibits β-catenin activity in synovial joint development and osteoarthritis. J. Clin. Invest., 2016, 126: 1649-1663

[79] Coveney CR, et al.. Role of ciliary protein intraflagellar transport protein 88 in the regulation of cartilage thickness and osteoarthritis development in mice. Arthritis Rheumatol. Hoboken NJ, 2022, 74: 49-59

[80] Malsin ES, Kim S, Lam AP, Gottardi CJ. Macrophages as a source and recipient of Wnt signals. Front. Immunol., 2019, 10: 1813

[81] Campbell TM, et al.. Mesenchymal stem cell alterations in bone marrow lesions in patients with hip osteoarthritis. Arthritis Rheumatol. Hoboken NJ, 2016, 68: 1648-1659

[82] Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat. Biotechnol., 2018, 36: 411-420

[83] Stuart T, et al.. Comprehensive integration of single-cell data. Cell, 2019, 177: 1888-1902.e21

[84] Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics J. Integr. Biol., 2012, 16: 284-287

[85] La Manno G, et al.. RNA velocity of single cells. Nature, 2018, 560: 494-498

[86] Castanza AS, et al.. Extending support for mouse data in the Molecular Signatures Database (MSigDB). Nat. Methods, 2023, 20: 1619-1620

[87] Wang P-E, et al.. Bushenhuoxue formula attenuates cartilage degeneration in an osteoarthritic mouse model through TGF-β/MMP13 signaling. J. Transl. Med., 2018, 16 72

[88] Kang D, et al.. Stress-activated miR-204 governs senescent phenotypes of chondrocytes to promote osteoarthritis development. Sci. Transl. Med., 2019, 11: eaar6659

Funding

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

Zhejiang Provincial Medical and Health Science and Technology Fund (grant no. 2024KY1223) Research Project of Zhejiang Chinese Medical University (grant no. 2023JKZKTS34)

Natural Science Foundation of Zhejiang Province (Zhejiang Provincial Natural Science Foundation)(LD22C060002)

Project of Chunyan Special Fund for Chinese Medicine Development of Zhejiang Chinese Medical University (grant no. CY202305)

RIGHTS & PERMISSIONS

The Author(s)