Current issue
Archive
Manuscripts accepted
About the Journal
Editorial office
Editorial board
Section Editors
Abstracting and indexing
Subscription
Contact
Ethical standards and procedures
Most read articles
Instructions for authors
Article Processing Charge (APC)
Regulations of paying article processing charge (APC)
OSTEOPOROSIS / CLINICAL RESEARCH
Exploring causal correlations between inflammatory response-related genes and osteoporosis: a Multi-Omics Mendelian Randomization Study
1
Affiliated Hospital of Shandong University of Traditional Chinese Medicine, Jinan, China
Submission date: 2024-10-30
Final revision date: 2025-05-13
Acceptance date: 2025-05-17
Online publication date: 2025-06-25
Corresponding author
KEYWORDS
TOPICS
ABSTRACT
Introduction:
The relationship between the inflammatory response (IR) and osteoporosis (OP) has been the subject of extensive research; however, their genetic link remains unclear. This study used IR-related genes as instrumental variables (IVs) to represent IR, while summary data of OP served as the outcome to explore their genetic relationship.
Material and methods:
IR-related genes were retrieved from the GeneCards database. OP transcriptome datasets were collected from the Gene Expression Omnibus (GEO) database and meta-analyzed to identify differentially expressed genes (DEGs) related to IR in OP. Genetic proxy instruments for IR-related genes were derived from studies of corresponding gene expression (n = 31,684) and DNA methylation (n = 1,980) quantitative trait loci (eQTLs and mQTLs), respectively. Aggregated data for OP (1,351 OP cases and 209,313 controls) were extracted from the largest genome-wide association study (GWAS) of OP. We integrated QTL data with OP GWAS data to estimate their genetic associations using summary data-based Mendelian randomization analysis (SMR). Additionally, Bayesian colocalization analysis was employed to reveal the potential relationships between IR gene expression and inflammatory factors, as well as various hormones. Finally, to further validate whether the statistical evidence provided in GWAS comprised true-positive findings, a replication study (1,955 cases and 278,169 controls) was conducted here through genotype-phenotype associations.
Results:
A meta-analysis of four datasets identified 115 IR-related DEGs in OP out of 612 IR-related genes. Through SMR analysis, we found that the expression levels of two IR-related genes were associated with OP risk. Specifically, elevated gene expression levels of FAS (odds ratio (OR) = 1.094; 95% confidence interval (CI) = 0.892–1.341; false discovery rate (FDR) = 0.034) increased the risk of OP. Conversely, increased expression levels of CHUK decreased the risk of OP (OR = 0.518; 95% CI = 0.424–0.637; FDR = 0.039). Colocalization analysis identified potential interactions between the FAS gene and estradiol (PP.H4 = 0.95) as well as interleukin-1α (IL-1α) (PP.H4 = 0.65). Potential interactions were also observed between the CHUK gene and growth hormone (PP.H4 = 0.59) as well as macrophage inflammatory protein-1α (MIP-1α) (PP.H4 = 0.62). In addition, consistent results were observed in the replication study, further demonstrating the reliability of our findings.
Conclusions:
This multi-omics integration study revealed a genetic link between IR and OP, as represented by IR-related genes, and provided new insights into the potential pathogenic mechanisms of OP. Additionally, these identified candidate genes offer avenues for future targeted functional studies aimed at developing appropriate therapeutic interventions and preventing OP.
The relationship between the inflammatory response (IR) and osteoporosis (OP) has been the subject of extensive research; however, their genetic link remains unclear. This study used IR-related genes as instrumental variables (IVs) to represent IR, while summary data of OP served as the outcome to explore their genetic relationship.
Material and methods:
IR-related genes were retrieved from the GeneCards database. OP transcriptome datasets were collected from the Gene Expression Omnibus (GEO) database and meta-analyzed to identify differentially expressed genes (DEGs) related to IR in OP. Genetic proxy instruments for IR-related genes were derived from studies of corresponding gene expression (n = 31,684) and DNA methylation (n = 1,980) quantitative trait loci (eQTLs and mQTLs), respectively. Aggregated data for OP (1,351 OP cases and 209,313 controls) were extracted from the largest genome-wide association study (GWAS) of OP. We integrated QTL data with OP GWAS data to estimate their genetic associations using summary data-based Mendelian randomization analysis (SMR). Additionally, Bayesian colocalization analysis was employed to reveal the potential relationships between IR gene expression and inflammatory factors, as well as various hormones. Finally, to further validate whether the statistical evidence provided in GWAS comprised true-positive findings, a replication study (1,955 cases and 278,169 controls) was conducted here through genotype-phenotype associations.
Results:
A meta-analysis of four datasets identified 115 IR-related DEGs in OP out of 612 IR-related genes. Through SMR analysis, we found that the expression levels of two IR-related genes were associated with OP risk. Specifically, elevated gene expression levels of FAS (odds ratio (OR) = 1.094; 95% confidence interval (CI) = 0.892–1.341; false discovery rate (FDR) = 0.034) increased the risk of OP. Conversely, increased expression levels of CHUK decreased the risk of OP (OR = 0.518; 95% CI = 0.424–0.637; FDR = 0.039). Colocalization analysis identified potential interactions between the FAS gene and estradiol (PP.H4 = 0.95) as well as interleukin-1α (IL-1α) (PP.H4 = 0.65). Potential interactions were also observed between the CHUK gene and growth hormone (PP.H4 = 0.59) as well as macrophage inflammatory protein-1α (MIP-1α) (PP.H4 = 0.62). In addition, consistent results were observed in the replication study, further demonstrating the reliability of our findings.
Conclusions:
This multi-omics integration study revealed a genetic link between IR and OP, as represented by IR-related genes, and provided new insights into the potential pathogenic mechanisms of OP. Additionally, these identified candidate genes offer avenues for future targeted functional studies aimed at developing appropriate therapeutic interventions and preventing OP.
REFERENCES (68)
1.
Lane JM, Russell L, Khan SN. Osteoporosis. Clin Orthop Relat Res 2000; 372: 139-50.
2.
Odén A, McCloskey EV, Kanis JA, Harvey NC, Johansson H. Burden of high fracture probability worldwide: secular increases 2010-2040. Osteoporos Int 2015; 26: 2243-8.
3.
Sigurdsson G, Halldorsson BV, Styrkarsdottir U, Kristjansson K, Stefansson K. Impact of genetics on low bone mass in adults. J Bone Miner Res 2008; 23: 1584-90.
5.
Xue Y, Zhou Y, Li C, Zhang J, Liu F, Shi R. Causal relationship between Interleukin-27 expression levels and osteoporosis: a bidirectional mendelian randomization study. BMC Musculoskelet Disord 2024; 25: 680.
6.
Rifas L, Weitzmann MN. A novel T cell cytokine, secreted osteoclastogenic factor of activated T cells, induces osteoclast formation in a RANKL-independent manner. Arthritis Rheum 2009; 60: 3324-35.
7.
Amarasekara DS, Kim S, Rho J. Regulation of osteoblast differentiation by cytokine networks. Int J Mol Sci 2021; 22: 2851.
8.
Zhang R, Peng S, Zhu G. The role of secreted osteoclastogenic factor of activated T cells in bone remodeling. Jpn Dent Sci Rev 2022; 58: 227-32.
9.
Rong K, Lang Y, Zhou Y, et al. Risk genetic variants (IL-10) for osteoporosis in han population from Northwest China. J Inflamm Res 2023; 16: 1091-102.
10.
Ke K, Sul OJ, Chung SW, Suh JH, Choi HS. Lack of NOD2 attenuates ovariectomy-induced bone loss via inhibition of osteoclasts. J Endocrinol 2017; 235: 85-96.
11.
Xie J, Ma R, Xu X, et al. Identification of genetic association between mitochondrial dysfunction and knee osteoarthritis through integrating multi-omics: a summary data-based Mendelian randomization study. Clin Rheumatol 2024; 43: 3487-96.
12.
Xie N, Xie J, Wang Z, et al. The role of calcium, 25-hydroxyvitamin D, and parathyroid hormone in irritable bowel syndrome: a bidirectional two-sample mendelian randomization study. Nutrients 2022; 14: 5109.
13.
Wu Y, Zeng J, Zhang F, et al. Integrative analysis of omics summary data reveals putative mechanisms underlying complex traits. Nat Commun 2018; 9: 918.
14.
Zhu Z, Zhang F, Hu H, et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat Genet 2016; 48: 481-7.
15.
Lei SF, Wu S, Li LM, et al. An in vivo genome wide gene expression study of circulating monocytes suggested GBP1, STAT1 and CXCL10 as novel risk genes for the differentiation of peak bone mass. Bone 2009; 44: 1010-4.
16.
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: 668.
17.
Benisch P, Schilling T, Klein-Hitpass L, 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: e45142.
18.
Võsa U, Claringbould A, Westra HJ, et al. Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression. Nat Genet 2021; 53: 1300-10.
19.
Ollier W, Sprosen T, Peakman T. UK Biobank: from concept to reality. Pharmacogenomics 2005; 6: 639-46.
20.
Kurki MI, Karjalainen J, Palta P, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature 2023; 613: 508-18.
21.
Xu S, Li X, Zhang S, et al. Oxidative stress gene expression, DNA methylation, and gut microbiota interaction trigger Crohn’s disease: a multi-omics Mendelian randomization study. BMC Med 2023; 21: 179.
22.
Walker MD, Shane E. Postmenopausal osteoporosis. N Engl J Med 2023; 389: 1979-91.
23.
Yamada A, Arakaki R, Saito M, Kudo Y, Ishimaru N. Dual role of Fas/FasL-mediated signal in peripheral immune tolerance. Front Immunol 2017; 8: 403.
24.
Galli G, Vacher P, Ryffel B, Blanco P, Legembre P. Fas/CD95 signaling pathway in damage-associated molecular pattern (DAMP)-sensing receptors. Cells 2022; 11: 1438.
25.
Yan H, Hong Y, Cai Y. Association between FAS gene -670 A/G and -1377 G/A polymorphisms and the risk of autoimmune diseases: a meta-analysis. Biosci Rep 2020; 40: BSR20191197.
26.
Huang D, Xiao J, Deng X, et al. Association between Fas/FasL gene polymorphism and musculoskeletal degenerative diseases: a meta-analysis. BMC Musculoskelet Disord 2018; 19: 137.
27.
Glavan BJ, Holden TD, Goss CH, et al. Genetic variation in the FAS gene and associations with acute lung injury. Am J Respir Crit Care Med 2011; 183: 356-63.
28.
Zhang L, Cheng J, Su H, Wang Z, Dai W. Diagnostic value of circulating bone turnover markers osteocalcin, cathepsin K, and osteoprotegerin for osteoporosis in middle-aged and elderly postmenopausal women. Arch Med Sci 2024; 20: 1727-30.
29.
Wu X, Pan G, McKenna MA, Zayzafoon M, Xiong WC, McDonald JM. RANKL regulates Fas expression and Fas-mediated apoptosis in osteoclasts. J Bone Miner Res 2005; 20: 107-16.
30.
Udagawa N, Koide M, Nakamura M, et al. Osteoclast differentiation by RANKL and OPG signaling pathways.
32.
Zhao G, Lyu J, Huang X, et al. The role and underlying mechanism of dental pulp stem cell-derived exosomal miR-31 in the treatment of osteoarthritis by targeting mTOR to enhance chondrocyte autophagy levels. Arch Med Sci 2024; 20: 1680-94.
33.
Morikawa N, Kato Y, Takeshita N, Shimizu Y. Pharmacological characterization of AS2690168, a novel small molecule RANKL signal transduction inhibitor. Eur J Pharmacol 2022; 924: 174941.
34.
Xiong J, Piemontese M, Thostenson JD, Weinstein RS, Manolagas SC, O’Brien CA. Osteocyte-derived RANKL is a critical mediator of the increased bone resorption caused by dietary calcium deficiency. Bone 2014; 66: 146-54.
35.
Chen T, Wang H, Jiang C, Lu Y. PKD1 alleviates oxidative stress-inhibited osteogenesis of rat bone marrow-derived mesenchymal stem cells through TAZ activation. J Cell Biochem 2021; 122: 1715-25.
36.
Astudillo P, Ríos S, Pastenes L, Pino AM, Rodríguez JP. Increased adipogenesis of osteoporotic human-mesenchymal stem cells (MSCs) characterizes by impaired leptin action. J Cell Biochem 2008; 103: 1054-65.
37.
Chen K, Qiu P, Yuan Y, et al. Pseurotin A inhibits osteoclastogenesis and prevents ovariectomized-induced bone loss by suppressing reactive oxygen species. Theranostics 2019; 9: 1634-50.
38.
Agidigbi TS, Kim C. Reactive oxygen species in osteoclast differentiation and possible pharmaceutical targets of ROS-mediated osteoclast diseases. Int J Mol Sci 2019; 20: 3576.
39.
Liu HD, Ren MX, Li Y, et al. Melatonin alleviates hydrogen peroxide induced oxidative damage in MC3T3-E1 cells and promotes osteogenesis by activating SIRT1. Free Radic Res 2022; 56: 63-76.
40.
Gürler EB, Çilingir-Kaya ÖT, Peker Eyüboglu I, et al. Melatonin supports alendronate in preserving bone matrix and prevents gastric inflammation in ovariectomized rats. Cell Biochem Funct 2019; 37: 102-12.
41.
Rippo MR, Babini L, Prattichizzo F, et al. Low FasL levels promote proliferation of human bone marrow-derived mesenchymal stem cells, higher levels inhibit their differentiation into adipocytes. Cell Death Dis 2013; 4: e594.
42.
Wang L, Liu S, Zhao Y, et al. Osteoblast-induced osteoclast apoptosis by fas ligand/FAS pathway is required for maintenance of bone mass. Cell Death Differ 2015; 22: 1654-64.
43.
Harding O, Holzer E, Riley JF, Martens S, Holzbaur ELF. Damaged mitochondria recruit the effector NEMO to activate NF-B signaling. Mol Cell 2023; 83: 3188-204.e7.
44.
Li X, Hu Y. Attribution of NF-B activity to CHUK/IKK-involved carcinogenesis. Cancers 2021; 13: 1411.
45.
Takeda K, Takeuchi O, Tsujimura T, et al. Limb and skin abnormalities in mice lacking IKKalpha. Science 1999; 284: 313-6.
46.
Antonia RJ, Hagan RS, Baldwin AS. Expanding the view of IKK: new substrates and new biology. Trends Cell Biol 2021; 31: 166-78.
47.
Walkon LL, Strubbe-Rivera JO, Bazil JN. Calcium overload and mitochondrial metabolism. Biomolecules 2022; 12: 1891.
48.
Suh J, Kim NK, Shim W, et al. Mitochondrial fragmentation and donut formation enhance mitochondrial secretion to promote osteogenesis. Cell Metab 2023; 35: 345-60.e7.
49.
Zhai Q, Chen X, Fei D, et al. Nanorepairers rescue inflammation-induced mitochondrial dysfunction in mesenchymal stem cells. Adv Sci 2022; 9: e2103839.
50.
Eijken M, Meijer IMJ, Westbroek I, et al. Wnt signaling acts and is regulated in a human osteoblast differentiation dependent manner. J Cell Biochem 2008; 104: 568-79.
51.
Liu F, Yuan Y, Bai L, et al. LRRc17 controls BMSC senescence via mitophagy and inhibits the therapeutic effect of BMSCs on ovariectomy-induced bone loss. Redox Biol 2021; 43: 101963.
52.
Ma W, Li C. Enhancing postmenopausal osteoporosis: a study of KLF2 transcription factor secretion and PI3K-Akt signaling pathway activation by PIK3CA in bone marrow mesenchymal stem cells. Arch Med Sci 2024; 20: 918-37.
53.
Liang Y, Li W, Li X, Nan J. Involvement of PI3K/Akt/-catenin signaling in schisandrin B-mitigated bone deterioration in an experimental rat model of estrogen deficiency. Arch Med Sci 2023; 19: 1520-9.
54.
Yang X, Jiang T, Wang Y, Guo L. The role and mechanism of SIRT1 in resveratrol-regulated osteoblast autophagy in osteoporosis rats. Sci Rep 2019; 9: 18424.
55.
Zheng Z, Zhang X, Huang B, et al. Site-1 protease controls osteoclastogenesis by mediating LC3 transcription. Cell Death Differ 2021; 28: 2001-18.
56.
Shao B, Fu X, Yu Y, Yang D. Regulatory effects of miRNA 181a on FasL expression in bone marrow mesenchymal stem cells and its effect on CD4+T lymphocyte apoptosis in estrogen deficiency induced osteoporosis. Mol Med Rep 2018; 18: 920-30.
57.
Tella SH, Gallagher JC. Prevention and treatment of postmenopausal osteoporosis. J Steroid Biochem Mol Biol 2014; 142: 155-70.
58.
Yang S, Zhu L, Luo F, Jiang R, Yu L. Association between female reproductive characteristics and bone mineral density: evidence from Mendelian randomization. Arch Med Sci 2024; 20: 2063-6.
59.
Boyce BF, Aufdemorte TB, Garrett IR, Yates AJ, Mundy GR. Effects of interleukin-1 on bone turnover in normal mice. Endocrinology 1989; 125: 1142-50.
60.
Landin-Wilhelmsen K, Nilsson A, Bosaeus I, Bengtsson BA. Growth hormone increases bone mineral content in postmenopausal osteoporosis: a randomized placebo-controlled trial. J Bone Miner Res 2003; 18: 393-405.
61.
Kukita T, Nomiyama H, Ohmoto Y, et al. Macrophage inflammatory protein-1 alpha (LD78) expressed in human bone marrow: its role in regulation of hematopoiesis and osteoclast recruitment. Lab Investig J Tech Methods Pathol 1997; 76: 399-406.
62.
Kawao N, Kaji H. Interactions between muscle tissues and bone metabolism. J Cell Biochem 2015; 116: 687-95.
63.
Chen J, Xu F, Ruan X, et al. Therapeutic targets for inflammatory bowel disease: proteome-wide Mendelian randomization and colocalization analyses. EBioMedicine 2023; 89: 104494.
64.
Gudjonsson A, Gudmundsdottir V, Axelsson GT, et al. A genome-wide association study of serum proteins reveals shared loci with common diseases. Nat Commun 2022; 13: 480.
65.
Shi J, Wang Z, Yi M, et al. Evidence based on Mendelian randomization and colocalization analysis strengthens causal relationships between structural changes in specific brain regions and risk of amyotrophic lateral sclerosis. Front Neurosci 2024; 18: 1333782.
66.
Rasooly D, Peloso GM, Pereira AC, et al. Genome-wide association analysis and Mendelian randomization proteomics identify drug targets for heart failure. Nat Commun 2023; 14: 3826.
67.
Arvanitis M, Tayeb K, Strober BJ, Battle A. Redefining tissue specificity of genetic regulation of gene expression in the presence of allelic heterogeneity. Am J Hum Genet 2022; 109: 223-39.
68.
Emilsson V, Ilkov M, Lamb JR, et al. Co-regulatory networks of human serum proteins link genetics to disease. Science 2018; 361: 769-73.
Share
RELATED ARTICLE
| eISSN: | 1896-9151 |
| ISSN: | 1734-1922 |
We process personal data collected when visiting the website. The function of obtaining information about users and their behavior is carried out by voluntarily entered information in forms and saving cookies in end devices. Data, including cookies, are used to provide services, improve the user experience and to analyze the traffic in accordance with the Privacy policy. Data are also collected and processed by Google Analytics tool (more).
You can change cookies settings in your browser. Restricted use of cookies in the browser configuration may affect some functionalities of the website.
You can change cookies settings in your browser. Restricted use of cookies in the browser configuration may affect some functionalities of the website.
