Multi-omics Analysis to Identify Key Immune Genes for Osteoporosis based on Machine Learning and Single-cell Analysis

Abstract

Objective: Osteoporosis is a severe bone disease with a complex pathogenesis involving various immune processes. With the in-depth understanding of bone immune mechanisms, discovering new therapeutic targets is crucial for the prevention and treatment of osteoporosis. This study aims to explore novel bone immune markers related to osteoporosis based on single-cell and transcriptome data, utilizing bioinformatics and machine learning methods, in order to provide novel strategies for the diagnosis and treatment of the disease.

Methods: Single cell and transcriptome data sets were acquired from Gene Expression Omnibus (GEO). The data was then subjected to cell communication analysis, pseudotime analysis, and high dimensional WGCNA (hdWGCNA) analysis to identify key immune cell subpopulations and module genes. Subsequently, ConsensusClusterPlus analysis was performed on the key module genes to identify different diseased subgroups in the osteoporosis (OP) training set samples. The immune characteristics between subgroups were evaluated using Cibersort, EPIC, and MCP counter algorithms. OP’s hub genes were screened using 10 machine learning algorithms and 113 algorithm combinations. The relationship between hub genes and immunity and pathways was established by evaluating the immune and pathway scores of the training set samples through the ESTIMATE, MCP-counter, and ssGSEA algorithms. Real-time fluorescence quantitative PCR (RT-qPCR) testing was conducted on serum samples collected from osteoporosis patients and healthy adults.

Results: In OP samples, the proportions of bone marrow-derived mesenchymal stem cells (BM-MSCs) and neutrophils increased significantly by 6.73% (from 24.01% to 30.74%) and 6.36% (from 26.82% to 33.18%), respectively. We found 16 intersection genes and four hub genes (DND1, HIRA, SH3GLB2, and F7). RT-qPCR results showed reduced expression levels of DND1, HIRA, and SH3GLB2 in clinical blood samples of OP patients. Moreover, the four hub genes showed positive correlations with neutrophils (0.65–0.90), immature B cells (0.76–0.92), and endothelial cells (0.79–0.87), while showing negative correlations with myeloid-derived suppressor cells (negative 0.54–0.73), T follicular helper cells (negative 0.71–0.86), and natural killer T cells (negative 0.75–0.85).

Conclusion: Neutrophils play a crucial role in the occurrence and development of osteoporosis. The four hub genes potentially inhibit metabolic activities and trigger inflammation by interacting with other immune cells, thereby significantly contributing to the onset and diagnosis of OP.

Keywords

Bioinformatics / Immunology / Machine learning / Osteoporosis / Single cells analysis

References

[1] Zhang YL, Chen Q, Zheng L, Zhang ZW, Chen YJ, Dai YC, et al. Jianpi Qingchang Bushen decoction improves inflammatory response and metabolic bone disorder in inflammatory bowel disease-induced bone loss. World J Gastroenterol. 2022; 28(13): 1315–1328.

[2] Lin S, Wu J, Chen B, Li S, Huang H. Identification of a potential MiRNA-mRNA regulatory network for osteoporosis by using bioinformatics methods: a retrospective study based on the gene expression omnibus database. Front Endocrinol (Lausanne). 2022; 13: 844218.

[3] Zhang W, Zhao W, Li W, Geng Q, Zhao R, Yang Y, et al. The imbalance of cytokines and lower levels of tregs in elderly male primary osteoporosis. Front Endocrinol (Lausanne). 2022; 13: 779264.

[4] Sapra L, Shokeen N, Porwal K, Saini C, Bhardwaj A, Mathew M, et al. Bifidobacterium longum ameliorates ovariectomy-induced bone loss via enhancing anti-osteoclastogenic and immunomodulatory potential of regulatory B cells (Bregs). Front Immunol. 2022; 13: 875788.

[5] Zhang H, Feng J, Lin Z, Wang S, Wang Y, Dai S, et al. Identification and analysis of genes underlying bone mineral density by integrating microarray data of osteoporosis. Front Cell Dev Biol. 2020; 8: 798.

[6] Marycz K, Kowalczuk A, Turlej E, Zachanowicz E, Tomaszewska A, Kulpa-Greszta M, et al. Impact of Polyrhodanine manganese ferrite binary Nanohybrids (PRHD@MnFe₂O₄) on osteoblasts and osteoclasts activities-a key factor in osteoporosis treatment. Materials (Basel). 2022; 15(11): 3990.

[7] Yu T, Xiong Y, Luu S, You X, Li B, Xia J, et al. The shared KEGG pathways between icariin-targeted genes and osteoporosis. Aging (Albany NY). 2020; 12(9): 8191–8201.

[8] Changani H, Parikh P. Molecular insights for an anti-osteoporotic properties of Litsea glutinosa on saos-2 cells: an in-vitro approach. J Ayurveda Integr Med. 2021; 13(2): 100501.

[9] Ding P, Gao C, Gao Y, Liu D, Li H, Xu J, et al. Osteocytes regulate senescence of bone and bone marrow. Elife. 2022; 11: e81480.

[10] Shukla P, Mansoori MN, Kakaji M, Shukla M, Gupta SK, Singh D. Interleukin 27 (IL-27) alleviates bone loss in estrogen-deficient conditions by induction of early growth response-2 gene. J Biol Chem. 2017; 292(11): 4686–4699.

[11] Frase D, Lee C, Nachiappan C, Gupta R, Akkouch A. The inflammatory contribution of B-lymphocytes and neutrophils in progression to osteoporosis. Cells. 2023; 12(13): 1744.

[12] Wang X, Liu X, He P, Guan K, Yang Y, Lei Y, et al. The imbalance of mitochondrial homeostasis of peripheral blood-derived macrophages mediated by MAFLD may impair the walking ability of elderly patients with osteopenia. Oxid Med Cell Longev. 2022; 2022: 5210870.

[13] Zheng L, Zhuang Z, Li Y, Shi T, Fu K, Yan W, et al. Bone targeting antioxidative nano-iron oxide for treating postmenopausal osteoporosis. Bioact Mater. 2022; 14: 250–261.

[14] Sun K, Li M, Wu Y, Wu Y, Zeng Y, Zhou S, et al. Exploring causal relationships between leukocyte telomere length, sex hormone-binding globulin levels, and osteoporosis using univariable and multivariable mendelian randomization. Orthop Surg. 2024; 16(2): 320–328.

[15] Tian L, Lu L, Meng Y. Bone marrow stromal stem cell fate decision: a potential mechanism for bone marrow adipose increase with aging-related osteoporosis. Curr Mol Med. 2023; 23(10): 1046–1057.

[16] Li CJ, Xiao Y, Sun YC, He WZ, Liu L, Huang M, et al. Senescent immune cells release grancalcin to promote skeletal aging. Cell Metab. 2021; 33(10): 1957–1973.e1956.

[17] Bi J, Zhang C, Lu C, Mo C, Zeng J, Yao M, et al. Age-related bone diseases: role of inflammaging. J Autoimmun. 2024; 143: 103169.

[18] Ho WC, Chang CC, Wu WT, Lee RP, Yao TK, Peng CH, et al. Effect of osteoporosis treatments on osteoarthritis progression in postmenopausal women: a review of the literature. Curr Rheumatol Rep. 2024; 26: 188–195.

[19] van Dalen SCM, Blom AB, Walgreen B, Slöetjes AW, Helsen MMA, Geven EJW, et al. IL-1β-mediated activation of adipose-derived mesenchymal stromal cells results in PMN reallocation and enhanced phagocytosis: a possible mechanism for the reduction of osteoarthritis pathology. Front Immunol. 2019; 10: 1075.

[20] Wang Z, Li X, Yang J, Gong Y, Zhang H, Qiu X, et al. Single-cell RNA sequencing deconvolutes the in vivo heterogeneity of human bone marrow-derived mesenchymal stem cells. Int J Biol Sci. 2021; 17(15): 4192–4206.

[21] Zhou Y, Zhu W, Zhang L, Zeng Y, Xu C, Tian Q, et al. Transcriptomic data identified key transcription factors for osteoporosis in caucasian women. Calcif Tissue Int. 2018; 103(6): 581–588.

[22] Zhou Y, Gao Y, Xu C, Shen H, Tian Q, Deng HW. A novel approach for correction of crosstalk effects in pathway analysis and its application in osteoporosis research. Sci Rep. 2018; 8(1): 668.

[23] Benisch P, Schilling T, Klein-Hitpass L, Frey SP, Seefried L, Raaijmakers N, et al. The transcriptional profile of mesenchymal stem cell populations in primary osteoporosis is distinct and shows overexpression of osteogenic inhibitors. PLoS One. 2012; 7(9): e45142.

[24] Xiao P, Chen Y, Jiang H, Liu YZ, Pan F, Yang TL, et al. In vivo genome-wide expression study on human circulating B cells suggests a novel ESR1 and MAPK3 network for postmenopausal osteoporosis. J Bone Miner Res. 2008; 23(5): 644–654.

[25] Kazezian Z, Gawri R, Haglund L, Ouellet J, Mwale F, Tarrant F, et al. Gene expression profiling identifies interferon signalling molecules and IGFBP3 in human degenerative annulus Fibrosus. Sci Rep. 2015; 5: 15662.

[26] Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016; 34(5): 525–527.

[27] Jin S, Guerrero-Juarez CF, Zhang L, Chang I, Ramos R, Kuan CH, et al. Inference and analysis of cell-cell communication using CellChat. Nat Commun. 2021; 12(1): 1088.

[28] Aibar S, Gonzalez-Blas CB, Moerman T, et al. SCENIC: single-cell regulatory network inference and clustering. Nat Methods. 2017; 14(11): 1083–1086.

[29] Morabito S, Reese F, Rahimzadeh N, Miyoshi E, Swarup V. hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data. Cell Rep Methods. 2023; 3(6): 100498.

[30] Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000; 28(1): 27–30.

[31] Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (Camb). 2021; 2(3): 100141.

[32] Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005; 102(43): 15545–15550.

[33] Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010; 26(12): 1572–1573.

[34] Kuncman Ł, Orzechowska M, Stawiski K, Masłowski M, Ciążyńska M, Gottwald L, et al. The kinetics of FMS-related tyrosine kinase 3 ligand (Flt-3L) during chemoradiotherapy suggests a potential gain from the earlier initiation of immunotherapy. Cancers (Basel). 2022; 14(16): 3844.

[35] Liu Z, Liu L, Weng S, Guo C, Dang Q, Xu H, et al. Machine learning-based integration develops an immune-derived lncRNA signature for improving outcomes in colorectal cancer. Nat Commun. 2022; 13(1): 816.

[36] Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013; 4: 2612.

[37] Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016; 17(1): 218.

[38] Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013; 14: 7.

[39] Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017; 18(1): 248–262.

[40] Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019; 47(W1): W234–W241.

[41] Zhang Y, Duan Z, Guan Y, Xu T, Fu Y, Li G. Identification of 3 key genes as novel diagnostic and therapeutic targets for OA and COVID-19. Front Immunol. 2023; 14: 1167639.

[42] Ji Q, Zheng Y, Zhang G, Hu Y, Fan X, Hou Y, et al. Single-cell RNA-seq analysis reveals the progression of human osteoarthritis. Ann Rheum Dis. 2019; 78(1): 100–110.

[43] Huber AK, Patel N, Pagani CA, Marini S, Padmanabhan KR, Matera DL, et al. Immobilization after injury alters extracellular matrix and stem cell fate. J Clin Invest. 2020; 130(10): 5444–5460.

[44] Bai J, Ge G, Wang Q, Li W, Zheng K, Xu Y, et al. Engineering stem cell recruitment and osteoinduction via bioadhesive molecular mimics to improve osteoporotic bone-implant integration. Research. 2022; 2022: 9823784.

[45] Kim S, Lee H, Hong J, et al. Bone-targeted delivery of cell-penetrating-RUNX2 fusion protein in osteoporosis model. Adv Sci. 2023; 10(28): e2301570.

[46] Yamaji M, Jishage M, Meyer C, Suryawanshi H, der E, Yamaji M, et al. DND1 maintains germline stem cells via recruitment of the CCR4-NOT complex to target mRNAs. Nature. 2017; 543(7646): 568–572.

[47] Shindo T, Doi S, Nakashima A, Sasaki K, Arihiro K, Masaki T. TGF-β1 promotes expression of fibrosis-related genes through the induction of histone variant H3.3 and histone chaperone HIRA. Sci Rep. 2018; 8(1): 14060.

[48] Serfass JM, Takahashi Y, Zhou Z, Kawasawa YI, Liu Y, Tsotakos N, et al. Endophilin B2 facilitates endosome maturation in response to growth factor stimulation, autophagy induction, and influenza a virus infection. J Biol Chem. 2017; 292(24): 10097–10111.

[49] Kondreddy V, Wang J, Keshava S, Esmon CT, Rao LVM, Pendurthi UR. Factor VIIa induces anti-inflammatory signaling via EPCR and PAR1. Blood. 2018; 131(21): 2379–2392.

[50] Capurso C, Massaro M, Scoditti E, Vendemiale G, Capurso A. Vascular effects of the Mediterranean diet part I: anti-hypertensive and anti-thrombotic effects. Vascul Pharmacol. 2014; 63(3): 118–126.

[51] Baker KS, Kopec AK, Pant A, Poole LG, Cline-Fedewa H, Ivkovich D, et al. Direct amplification of tissue factor:factor VIIa procoagulant activity by bile acids drives intrahepatic coagulation. Arterioscler Thromb Vasc Biol. 2019; 39(10): 2038–2048.

[52] Titanji K, Ofotokun I, Weitzmann MN. Immature/transitional B-cell expansion is associated with bone loss in HIV-infected individuals with severe CD4+ T-cell lymphopenia. Aids. 2020; 34(10): 1475–1483.

[53] Titanji K, Vunnava A, Sheth AN, Delille C, Lennox JL, Sanford SE, et al. Dysregulated B cell expression of RANKL and OPG correlates with loss of bone mineral density in HIV infection. PLoS Pathog. 2014; 10(10): e1004497.

[54] Tschaffon-Müller MEA, Kempter E, Steppe L, Kupfer S, Kuhn MR, Gebhard F, et al. Neutrophil-derived catecholamines mediate negative stress effects on bone. Nat Commun. 2023; 14(1): 3262.

[55] Fischer V, Haffner-Luntzer M. Interaction between bone and immune cells: implications for postmenopausal osteoporosis. Semin Cell Dev Biol. 2022; 123: 14–21.

[56] Saxena Y, Routh S, Mukhopadhaya A. Immunoporosis: role of innate immune cells in osteoporosis. Front Immunol. 2021; 12: 687037.

[57] Wampfler J, Federzoni EA, Torbett BE, Fey MF, Tschan MP. The RNA binding proteins RBM38 and DND1 are repressed in AML and have a novel function in APL differentiation. Leuk Res. 2016; 41: 96–102.

[58] Yu X, Diamond SL. Fibrin modulates shear-induced NETosis in sterile occlusive thrombi formed under Haemodynamic flow. Thromb Haemost. 2019; 119(4): 586–593.

[59] Reyes-García AML, Aroca A, Arroyo AB, García-Barbera N, Vicente V, González-Conejero R, et al. Neutrophil extracellular trap components increase the expression of coagulation factors. Biomed Rep. 2019; 10(3): 195–201.

[60] Deng Z, Ng C, Inoue K, Chen Z, Xia Y, Hu X, et al. Def6 regulates endogenous type-I interferon responses in osteoblasts and suppresses osteogenesis. Elife. 2020; 9: e59659.

[61] Wang Q, Weng H, Xu Y, Ye H, Liang Y, Wang L, et al. Anti-osteoporosis mechanism of resistance exercise in ovariectomized rats based on transcriptome analysis: a pilot study. Front Endocrinol. 2023; 14: 1162415.

[62] Fang Z, Cheng G, He M, Lin Y. CYP27A1 deficiency promoted osteoclast differentiation. PeerJ. 2023; 11: e15041.

[63] Wang X, Pei Z, Hao T, Ariben J, Li S, He W, et al. Prognostic analysis and validation of diagnostic marker genes in patients with osteoporosis. Front Immunol. 2022; 13: 987937.

[64] Sluiter C, Kettenes-van den Bosch JJ, Hop E, van der Houwen OA, Underberg WJ, Bult A. Degradation study of the investigational anticancer drug clanfenur. Int J Pharm. 1999; 185(2): 227–235.

RIGHTS & PERMISSIONS

2024 The Author(s). Orthopaedic Surgery published by Tianjin Hospital and John Wiley & Sons Australia, Ltd.