- 1Orthopedic Research Institute, Department of Orthopedics, Sichuan University West China Hospital, Chengdu, Sichuan, China
- 2Department of Orthopedics, The Third Hospital of Mianyang, Sichuan Mental Health Center, Mianyang, Sichuan, China
Objective: To explore the genetic causal association between osteoporosis (OP) and iron status through Mendelian randomization (MR).
Methods: Publicly available genome-wide association study (GWAS) summary data were used for MR analysis with four iron status-related indicators (ferritin, iron, total iron binding capacity, and transferrin saturation) as exposures and three different types of OP (OP, OP with pathological fracture, and postmenopausal OP with pathological fracture) as outcomes. The inverse-variance weighted (IVW) method was used to analyze the genetic causal association between the four indicators of iron status and OP. The heterogeneity of MR results was determined using IVW and MR–Egger methods. The pleiotropy of MR results was determined using MR–Egger regression. A leave-one-SNP-out test was performed to determine whether the MR results were affected by a single nucleotide polymorphism (SNP). The weighted median method was conducted to further validate our results.
Results: Based on IVW, MR–Egger and weighted median models, we found no causal association between iron status (ferritin, iron, total iron binding capacity, or transferrin saturation) and OP (Pbeta > 0.05 in all models). IVW and MR–Egger analysis of OP with pathological fracture and iron status indicators showed no potential genetic causal association (Pbeta> 0.05 in the two analyses). The results of the weighted median were consistent with those of IVW (Pbeta> 0.05 in all analyses). There was no potential genetic causal association between iron status and postmenopausal OP with pathological fracture based on serum iron (Pbeta>0.05 in all models). No heterogeneity or horizontal pleiotropy was found in any of the analyses. None of the leave-one-out tests in the analyses found any SNP that could affect the results of MR.
Conclusion: Our results demonstrate that there is no genetic causal association between OP and iron status, but the effects of other factors were not excluded.
Introduction
Osteoporosis (OP) is a systemic bone disease characterized by a reduction in bone mass and a deterioration of trabecular structure, resulting in decreased bone strength and an increased risk of fragility fractures (1). It is a chronic disease that exerts significant negative impacts on the health of elderly individuals, with a prevalence of 30-40% in women and a prevalence of 10-20% in men over 50 in mainland China (2). Established risk factors for OP include advanced age, endocrine disorders, malnutrition, obesity, and the use of drugs affecting bone metabolism (3). Moreover, some genetic mutations (three single nucleotide polymorphisms [SNPs] near the TNFRSF11B gene) contribute to the increased risk of OP (4). OP is occult and tends not to be diagnosed until the occurrence of fractures (commonly in the hip, caudal vertebrae, or forearm), which seriously compromise the life of patients. In Europe, the harm caused by OP complicated with fracture is second only to lung cancer (5). However, the current treatments of OP largely rely on medications such as bisphosphonates, calcitonin, oestrogen, and oestrogen receptor agonists, which have limited therapeutic effect and considerable adverse reactions (6). Therefore, further exploration of other risk factors for OP is imperative for the treatment and prevention of OP.
As an essential element for the human body, iron is indispensable for mitochondrial function, DNA synthesis and repair, and cell survival (7). These physiological activities acquire iron mainly through three pathways: the decomposition and destruction of ageing red blood cells by macrophages, absorption of iron in food to compensate for iron loss or increased demand, and the buffering effect of the iron reserve in the liver (8). Iron metabolism in the human body is mainly realized by regulating a series of iron metabolism-related proteins such as hepcidin, ferritin, transferrin, transferrin receptor, ferroportin 1 (FPN1), and divalent metal transporter 1 (DMT1). Problems in any one of these links lead to an imbalanced iron status in the body, resulting in a series of adverse effects (8). For example, iron deficiency can result in cognitive developmental defects in children and adverse pregnancy in women (9), whereas iron overload can damage multiple organs, including the liver, heart, and pancreas (10).
At present, the impact of iron status on OP remains unclear. Accumulating studies have revealed a certain correlation between iron status and OP at the biochemical level; i.e., an imbalance in iron status may increase the risk of OP (11–13). In contrast, some animal experiments have indicated that there is no correlation between iron status and OP (14). No final conclusion has yet been reached on the correlation between iron status and OP. Hence, it is necessary to conduct a more in-depth study on the correlation between iron status and OP at the gene level.
Mendelian randomization (MR) is an analytic approach to assessing the causal association between a modifiable exposure and a clinically relevant outcome (15). MR is based on Mendel's law of inheritance (parental alleles are randomly assigned to offspring) and uses genetic variants as instrumental variables (IVs), which can improve the limitations of some observational studies such as reverse causality, confounding factors, and various biases (16). Using genetic variation as an IV for exposure, MR studies can strengthen causal inferences about exposure-outcome associations by reducing confounding factors and reverse causality (17). To ensure the reliability of the results, MR must meet three assumptions: IVs are closely related to exposures, IVs are not related to other confounding factors, and IVs affect outcomes only through exposures and not through any other pathways (18). In recent years, accumulating evidence has demonstrated the reliability of two-sample MR analysis. For example, an existing MR analysis demonstrated that high iron levels have a positive causality with gout but a negative causality with rheumatoid arthritis (RA) (17). In addition, the genetic causal association between major depressive disorder (MDD) and osteoarthritis (OA) has also been documented in a recent MR study (19). However, few studies have focused on the association between iron status and OP using MR analysis. This study selected four serum biomarkers related to iron status (ferritin, iron, total iron binding capacity, and transferrin saturation) to explore the genetic causal association between iron status (exposure) and OP (outcome) through a two-sample MR analysis.
Materials and methods
Genome-wide association study summary data for iron status-related indicators
The genome-wide association study (GWAS) summary data of iron status-related indicators was obtained from a meta-analysis of three genome-wide association studies from Iceland, the UK and Denmark of blood levels of ferritin (N = 246,139), total iron binding capacity (N = 135,430), iron (N = 163,511) and transferrin saturation (N = 131,471). Information about the samples included in the study and data processing can be found in previous studies (20).
Genome-wide association study summary data for osteoporosis
Recent large-scale GWASs and meta-analyses of OP in European populations (including OP N = 300,147, OP with pathological fracture N = 239,702, and postmenopausal OP with pathological fracture N = 173,601) were obtained from the FinnGen consortium. All three subtypes of OP were defined by the code M13 in the International Classification of Diseases, Tenth Revision (ICD-10). Detailed information on the participants, genotyping, imputation, and quality control can be found on the FinnGen website (https://www.finngen.fi/en) (21).
Instrumental variable selection
We extracted genomic single-nucleotide polymorphisms (SNPs) associated with exposure (P<5×10-8). None of the instrumental SNPs were in linkage disequilibrium (LD). We performed the clumping process (R2 < 0.001, Magna window size = 10000 kb) to eliminate LD between the SNPs. The missing SNPs in the LD control group were also deleted. SNPs with a minor allele frequency (MAF) < 0.01 were removed. By default, if the SNP for a particular request does not exist in the resulting GWAS, the SNP (agent) with the requested SNP (target) in the LD is searched. The LD agent was defined using 1000 genomes of European sample data. In addition, to test whether there was a weak tool deviation in IV, we used the F statistic (F = R2 [n-k-1]/k [1-R2]), where R2 is the variance of exposure explained by selected instrumental variables (obtained from the MR Steiger directionality test), n is the sample size, and k is the total number of variables. If the F statistic of IV is much greater than 10, it indicates that the possibility of weak instrument variable bias is very small. In addition, we eliminated some IVs that may be related to confounding factors. These confounding factors included a lack of vitamin D and calcium, menopause and a lack of exercise.
Mendelian randomization analysis
The “TwoSampleMR” package (version 0.5.6) of R language was used to analyze the genetic causal association between the four indicators of iron status and OP using three methods: inverse-variance weighted (IVW), MR–Egger, and weighted median (22).
The results were mainly based on IVW. The IVW method adopted a meta-analysis combined with Wald estimators from different SNPs. If each genetic variant can be used as an effective IV, the IVW method provides a consistent estimate of the causal effect of exposure on the outcome.
IVW and MR–Egger were performed to determine the heterogeneity of MR analysis results. Cochran’s Q statistic was adopted for IVW analysis, and Rucker’s Q statistic was used for MR–Egger analysis (23). The pleiotropy of MR results was detected using MR–Egger regression with an intercept P value > 0.05 indicating pleiotropy deficiency (24). A leave-one-SNP-out analysis was also conducted to investigate the possibility that the causal association between exposures and outcomes was driven by a single SNP (25). Moreover, we calculated the weighted median to further identify the potential causal association between exposures and outcomes.
Among the five methods of MR analysis, we mainly focus on the results of the IVW analysis. We considered that there was a genetic causal association between exposure and outcome when P<0.05 for the IVW analysis results. After correcting for multiple testing, the significance threshold of this study was P<4.17×10-3 (0.05/12=4.17×10-3). IVW and MR–Egger examined heterogeneity, and if P>0.05, there was no heterogeneity in our MR analysis. When our MR analysis results were free of heterogeneity/pleiotropy, we considered the IVW analysis results to be reliable.
Results
After removing the SNPs of incompatible alleles, the details of all independent SNPs associated with exposure are shown in Supplementary Table 1. In our study, the F statistics of the instrumental variables associated with exposure were all greater than 10, indicating that the possibility of bias in weak instrumental variables was very small.
The causal association between iron status and OP
Based on IVW and the MR–Egger model, we found that there was no causal association between iron status (ferritin, iron, total iron binding capacity, or transferrin saturation) and OP (Pbeta > 0.05 in all models) (Table 1). In addition, the results of the weighted median showed that there was no potential causal association between OP and iron status indicators (Pbeta> 0.05 in all analyses) (Table 1).
The MR–Egger intercept in the analysis showed that there was no horizontal multiplicity (MR–Egger intercept p value > 0.05) (Table 2). The MR–Egger test showed no heterogeneity in the MR analysis of iron status (ferritin, iron, total iron binding capacity, and transferrin saturation) and OP indicators (the Q-p values of the IVW and MR–Egger were both greater than 0.05) (Table 3). The leave-one-out analysis showed that the causal estimation of OP and iron status indicators was not driven by any single SNP (Supplementary Figures S1–4).
MR analyses for OP with pathological fracture
IVW and MR–Egger analysis of OP with pathological fracture and iron status indicators showed no potential genetic causal association between the two (Pbeta> 0.05 in the two analyses) (Table 1). The results of the weighted median were consistent with those of IVW (Pbeta> 0.05 in all analyses) (Table 1).
In addition, no horizontal pleiotropy was found in this analysis (all intercept P values>0.05) (Table 2). The heterogeneity analysis found no heterogeneity in this analysis, with a Q-P value > 0.05 in the two analyses (Table 3). The results of the leave-one-out method in the initial analysis and verified analysis showed that no abnormal IV affected the overall results in the two analyses (Supplementary Figures S5–8).
MR analyses for postmenopausal OP with pathological fracture
There was no potential genetic causal association between iron status and postmenopausal OP with pathological fracture based on serum iron (Pbeta>0.05 in all models) (Table 1). In addition, the results of the weighted median did not find a causal association between postmenopausal OP and pathological fracture and iron status indicators (Pbeta>0.05 in all analyses) (Table 1).
No horizontal multiplicity was found in this analysis (intercept P values in all analyses were greater than 0.05) (Table 2). The heterogeneity analysis showed no heterogeneity in the MR analysis of postmenopausal OP with pathological fracture and iron status indicators (Q-P value>0.05 in the IVW and MR–Egger analyses) (Table 3). The leave-one-out analysis showed that the causal estimation of postmenopausal OP with pathological fracture and iron status indicators was not driven by any single SNP (Supplementary Figures S9–12).
Discussion
Studies have found that iron overload-induced ferroptosis in osteoblasts can inhibit osteogenesis and promote osteoporosis (26). Iron oxide nanoparticles (IONPs) can positively regulate bone metabolism in vitro, and daily administration of IONPs can alleviate oestrogen deficiency-induced osteoporosis by removing reactive oxygen species from the body (27). Although iron is strongly associated with OP, no studies have demonstrated a genetic causal association between iron status and OP. Studying the genetic causal association between iron status and OP is of considerable significance for research on the aetiology, mechanism and treatment of OP. This is the first study to explore the causality between iron status and OP through MR analysis. No evidence of positive or negative causality between iron status and OP was found, indicating that iron status had no genetic causality with OP. However, other factors (such as the environment) might also exert potential regulatory effects on the correlation between iron status and OP.
Iron status imbalance is manifested by iron overload or iron deficiency (28). Iron overload leads to an increase in serum iron, ferritin, and transferrin saturation but a decrease in transferrin, whereas iron deficiency shows the opposite trends (29). Previous studies have shown that these iron status-related indicators are related to OP or low bone mineral density (BMD). For example, the transferrin level in the serum of patients with osteoporotic hip fracture is lower than that of normal subjects (30). A cross-sectional study based on 4,000 women aged 12-49 found that serum ferritin is negatively correlated with BMD (31). Transferrin saturation shows a negative correlation with BMD in patients with transfusion-dependent beta-thalassemia (32). A controlled clinical trial revealed that the concentration of serum iron in OP patients is significantly higher than that in healthy controls, but it is believed that iron overload is a necessary but not sufficient condition for OP because iron deficiency may also affect OP (12). Moreover, a meta-analysis demonstrated a correlation between serum iron and OP and identified a low serum iron level as a risk factor for OP (33). Some scholars have proposed that both iron overload and iron deficiency may increase the risk of OP (13). High serum iron and ferritin levels may be beneficial to elderly individuals. A study based on 262 elderly rheumatoid arthritis (RA) patients found that the serum iron level of elderly RA patients was positively correlated with BMD (34). A study on elderly individuals aged 60 and above also revealed a positive correlation between serum ferritin and BMD (35). The association between iron status and OP obtained in previous population studies may be due to confounding factors such as sex and age (31, 35).
Notably, systemic iron overload is observed in hemochromatosis protein (HFE) knockout mice, but it does not have any effect on BMD and bone microstructure (14), which is consistent with our results that found no causal association between iron status and OP, at least at the genetic level. Although a polygenic risk score study found that ferritin is genetically correlated with systemic BMD (36), the sample was from the Caucasus (our sample is European), and the research method was not adequate to explain whether there is a causal association between the two. Hence, the rationality of our research results is still valid.
Iron overload-induced reactive oxygen species (ROS) can damage DNA, proteins, and lipids, eventually leading to cell death (37) and can also promote osteoclast differentiation and proliferation via the NF-κB signalling pathway (11). Furthermore, excessive iron inhibits osteoblast viability in a concentration-dependent manner. Mild iron deficiency facilitates osteoblast viability, whereas severe iron deficiency impedes osteoblast formation (38). As the main endogenous hormone regulating iron status, hepcidin can degrade ferroportin (FPN), the sole iron exporter on the cell membrane, resulting in an increase in intracellular iron levels (39). It should be noted that hypoxia and inflammation may affect the regulation of iron status by hepcidin (40), and OP is closely related to inflammation (41). Therefore, it is reasonable to speculate that some inflammatory factors affect the correlation between iron status and OP via hepcidin. Oestrogen affects iron status via hepcidin, and menopausal women often present with iron overload (42, 43). Moreover, oestrogen can also directly affect bone metabolism (44, 45). It is suggested that oestrogen may participate in the correlation between OP and iron status. In addition, some scholars have proposed the implication of environmental factors in the relationship between iron status and OP. A tendency towards OP (decreased BMD) has been found in some malnourished young people (45). It is well established that iron deficiency is a form of malnutrition, which suggests that a poor living environment with long-term malnutrition may be a common risk factor for iron status imbalance (long-term iron deficiency) and OP. To summarize, the association between OP and iron status is complex and dependent on multiple factors, but our results at least show that OP and iron status have no genetic causal association.
Observational epidemiological studies are prone to confounding factors, reverse causation and various biases and have generated findings that have proven to be unreliable indicators of causal effects (16). However, MR studies are free from the confounding factors (as in retrospective studies) and reverse causality of traditional epidemiological approaches (46). In addition, compared with observational epidemiological studies, MR analysis does not involve high measurement costs or a large number of appropriate biospecimens (16). Therefore, MR analysis has high reliability and is widely used in many studies (17, 47).
This study excludes a causal association between iron status and OP through MR analysis based on large-scale GWAS summary data, but it cannot be denied that iron status and OP may still be related. The MR method is not without limitations. First, MR analysis is heavily dependent on the reliable associations of genetic variants with the exposure(s) of interest, which are believed to have no effect on other phenotypes that might confound the association between the exposure and disease (48). In addition, the GWAS summary datasets used in this study were not stratified by population. The genotyping errors, phenotype misclassification, and confounding factors due to population stratification may cause spurious genetic associations, which will in fact be biased instruments for MR (15).
Conclusion
This study shows that there are no positive or negative genetic causal associations between iron status and OP, but the influence of factors other than heredity cannot be ruled out.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The source of the data was a publicly available database, and no human participants were involved; hence, ethical parameters are not applicable.
Author contributions
Conception and design: JX, JM, JC, BS. Administrative support: JM, HS, BS. Provision of study materials: JX, JM, JC, CZ, YW. Collection and assembly of data: JX, LW. Data analysis and interpretation: JX, SZ, HS, YL, ML. Manuscript writing: JX, JC. Final approval of manuscript: All authors. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (grant numbers 81974347 and 81802210); the Department of Science and Technology of Sichuan Province (grant number 2021YFS0122 and 2020YFS0139). Financial support had no impact on the outcomes of this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.996244/full#supplementary-material
References
1. Armas LA, Recker RR. Pathophysiology of osteoporosis: New mechanistic insights. Endocrinol Metab Clin North Am (2012) 41(3):475–86. doi: 10.1016/j.ecl.2012.04.006
2. Yu F, Xia W. The epidemiology of osteoporosis, associated fragility fractures, and management gap in China. Arch Osteoporos (2019) 14(1):32. doi: 10.1007/s11657-018-0549-y
3. Hirschfeld PH, Kinsella, Duque. Osteosarcopenia: Where bone, muscle, and fat collide. Osteoporos Int (2017) 28(10):2781–90. doi: 10.1007/s00198-017-4151-8
4. Richards JB, Rivadeneira F, Inouye M, Pastinen TM, Spector TD. Bone mineral density, osteoporosis, and osteoporotic fractures: A genome-wide association study. Lancet (London England) (2020) (2008) 371(9623):1505–12. doi: 10.1016/S0140-6736(08)60599-1
5. Johnell O, Kanis JA. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos Int (2006) 17(12):1726. doi: 10.1007/s00198-006-0172-4
6. Chandran T, Venkatachalam I. Efficacy and safety of denosumab compared to bisphosphonates in improving bone strength in postmenopausal osteoporosis: A systematic review. Singapore Med J (2019) 60(7):364–78. doi: 10.11622/smedj.2019028
7. Bogdan AR, Miyazawa M, Hashimoto K, Tsuji K. Regulators of iron homeostasis: New players in metabolism, cell death, and disease. Trends Biochem Sci (2016) 41(3):274–86. doi: 10.1016/j.tibs.2015.11.012
8. Muckenthaler MU, Rivella S, Hentze MW, Galy B. A red carpet for iron metabolism. Cell (2017) 168(3):344–61. doi: 10.1016/j.cell.2016.12.034
9. Camaschella C. Iron-deficiency anemia. New Engl J Med (2015) 372(19):1832–43. doi: 10.1056/NEJMra1401038
11. Zhang H, Wang A, Shen G, Wang X, Liu G, Yang F, et al. Hepcidin-induced reduction in iron content and pgc-1beta expression negatively regulates osteoclast differentiation to play a protective role in postmenopausal osteoporosis. Aging (Albany NY) (2021) 13(8):11296–314. doi: 10.18632/aging.202817
12. Liu B, Liu C, Zhong W, Song M, Du S, Su J. Reduced hepcidin level features osteoporosis. Exp Ther Med (2018) 16(3):1963–7. doi: 10.3892/etm.2018.6410
13. Ledesma-Colunga MG, Weidner H, Vujic Spasic M, Hofbauer LC, Baschant U, Rauner M. Shaping the bone through iron and iron-related proteins. Semin Hematol (2021) 58(3):188–200. doi: 10.1053/j.seminhematol.2021.06.002
14. Wagner A, Alan B, Yilmaz D, Ahmad M, Liu P, Tangudu NK, et al. Despite genetic iron overload, hfe-hemochromatosis mice do not show bone loss. JBMR Plus (2019) 3(9):e10206. doi: 10.1002/jbm4.10206
15. Smith GD, Ebrahim S. 'Mendelian randomization': Can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol (2003) 32(1):1–22. doi: 10.1093/ije/dyg070
16. Davey Smith G, Hemani G. Mendelian randomization: Genetic anchors for causal inference in epidemiological studies. Hum Mol Genet (2014) 23(R1):R89–98. doi: 10.1093/hmg/ddu328
17. Yuan S, Larsson S. Causal associations of iron status with gout and rheumatoid arthritis, but not with inflammatory bowel disease. Clin Nutr (2020) 39(10):3119–24. doi: 10.1016/j.clnu.2020.01.019
18. Vanderweele TJ, Tchetgen ET, Cornelis M, Kraft P. Methodological challenges in mendelian randomization. Epidemiology (2014) 25(3):427. doi: 10.1097/EDE.0000000000000081
19. Zhang F, Rao S, Baranova A. Shared genetic liability between major depressive disorder and osteoarthritis. Bone Joint Res (2022) 11(1):12–22. doi: 10.1302/2046-3758.111.bjr-2021-0277.r1
20. Bell S, Rigas AS, Magnusson MK, Ferkingstad E, Allara E, Bjornsdottir G, et al. A genome-wide meta-analysis yields 46 new loci associating with biomarkers of iron homeostasis. Commun Biol (2021) 4(1):156. doi: 10.1038/s42003-020-01575-z
21. Canela-Xandri O, Rawlik K, Tenesa A. An atlas of genetic associations in uk biobank. Nat Genet (2018) 50(11):1593–9. doi: 10.1038/s41588-018-0248-z
22. Evans DM, Davey Smith G. Mendelian randomization: New applications in the coming age of hypothesis-free causality. Annu Rev Genomics Hum Genet (2015) 16:327–50. doi: 10.1146/annurev-genom-090314-050016
23. Shu MJ, Li JR, Zhu YC, Shen H. Migraine and ischemic stroke: A mendelian randomization study. Neurol Ther (2022) 11(1):237–46. doi: 10.1007/s40120-021-00310-y
24. Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the Mr-egger method. Eur J Epidemiol (2017) 32(5):377–89. doi: 10.1007/s10654-017-0255-x
25. Lee YH. Causal association between smoking behavior and the decreased risk of osteoarthritis: A mendelian randomization. Z Rheumatol (2019) 78(5):461–6. doi: 10.1007/s00393-018-0505-7
26. Jiang Z, Wang H, Qi G, Jiang C, Chen K, Yan Z. Iron overload-induced ferroptosis of osteoblasts inhibits osteogenesis and promotes osteoporosis: An in vitro and in vivo study. IUBMB Life (2022) 74(11):1052–69. doi: 10.1002/iub.2656
27. 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–61. doi: 10.1016/j.bioactmat.2021.11.012
28. Hentze MW, Muckenthaler MU, Andrews NC. Balancing acts. Cell (2004) 117(3):285–97. doi: 10.1016/S0092-8674(04)00343-5
29. Munoz M, Garcia-Erce JA, Remacha AF. Disorders of iron metabolism. part ii: Iron deficiency and iron overload. J Clin Pathol (2011) 64(4):287–96. doi: 10.1136/jcp.2010.086991
30. Zamir M. Nutritional assessment in elderly patients with osteoporotic hip fractures. Pakistan J Surg (2015) 31:192.
31. Lu M, Liu Y, Shao M, Tesfaye GC, Yang S. Associations of iron intake, serum iron and serum ferritin with bone mineral density in women: The national health and nutrition examination survey, 2005-2010. Calcif Tissue Int (2020) 106(3):232–8. doi: 10.1007/s00223-019-00627-9
32. Atmakusuma TD, Tenggara JB. Correlation of transferrin saturation and serum ferritin with bone mass density in adult transfusion dependent beta-thalassemia patients. J Blood Med (2021) 12:827–32. doi: 10.2147/JBM.S328547
33. Zheng J, Mao X, Ling J, He Q, Quan J. Low serum levels of zinc, copper, and iron as risk factors for osteoporosis: A meta-analysis. Biol Trace Elem Res (2014) 160(1):15–23. doi: 10.1007/s12011-014-0031-7
34. Sato H, Takai C, Kazama JJ, Wakamatsu A, Hasegawa E, Kobayashi D, et al. Serum hepcidin level, iron metabolism and osteoporosis in patients with rheumatoid arthritis. Sci Rep (2020) 10(1):9882. doi: 10.1038/s41598-020-66945-3
35. Babaei M, Bijani A, Heidari P, Hosseini SR, Heidari B. Serum ferritin levels and bone mineral density in the elderly. Caspian J Intern Med (2018) 9(3):232–8. doi: 10.22088/cjim.9.3.232
36. Liang X, Du Y, Wen Y, Liu L, Li P, Zhao Y, et al. Assessing the genetic correlations between blood plasma proteins and osteoporosis: A polygenic risk score analysis. Calcif Tissue Int (2019) 104(2):171–81. doi: 10.1007/s00223-018-0483-4
37. Li GF, Pan YZ, Sirois P, Li K, Xu YJ. Iron homeostasis in osteoporosis and its clinical implications. Osteoporos Int (2012) 23(10):2403–8. doi: 10.1007/s00198-012-1982-1
38. Toxqui L, Vaquero MP. Chronic iron deficiency as an emerging risk factor for osteoporosis: A hypothesis. Nutrients (2015) 7(4):2324–44. doi: 10.3390/nu7042324
39. Nemeth E, Tuttle MS, Powelson J, Vaughn MB, Donovan A, Ward DM, et al. Hepcidin regulates cellular iron efflux by binding to ferroportin and inducing its internalization. Science (2004) 306(5704):2090–3. doi: 10.1126/science.1104742
40. Wessling-Resnick M. Iron homeostasis and the inflammatory response. Annu Rev Nutr (2010) 30:105–22. doi: 10.1146/annurev.nutr.012809.104804
41. Lacativa PG, Farias ML. Osteoporosis and inflammation. Arq Bras Endocrinol Metabol (2010) 54(2):123–32. doi: 10.1590/s0004-27302010000200007
42. Chen B, Li GF, Shen Y, Huang XI, Xu YJ. Reducing iron accumulation: A potential approach for the prevention and treatment of postmenopausal osteoporosis. Exp Ther Med (2015) 10(1):7–11. doi: 10.3892/etm.2015.2484
43. Qian Y, Yin C, Chen Y, Zhang S, Jiang L, Wang F, et al. Estrogen contributes to regulating iron metabolism through governing ferroportin signaling Via an estrogen response element. Cell Signal (2015) 27(5):934–42. doi: 10.1016/j.cellsig.2015.01.017
44. Khosla S, Melton LJ 3rd, Riggs BL. The unitary model for estrogen deficiency and the pathogenesis of osteoporosis: Is a revision needed? J Bone Miner Res (2011) 26(3):441–51. doi: 10.1002/jbmr.262
45. Thornton D, Gordon CM. Restrictive eating disorders and skeletal health in adolescent girls and young women. Calcif Tissue Int (2017) 100(5):449–60. doi: 10.1007/s00223-016-0164-0
46. Cai Y, Zhang G, Liang J, Jing Z, Zhang R, Lv L, et al. Causal relationships between osteoarthritis and senile central nerve system dysfunction: A bidirectional two-sample mendelian randomization study. Front Aging Neurosci (2021) 13:793023. doi: 10.3389/fnagi.2021.793023
47. Zhang J. Mendelian randomization study implies causal linkage between telomere length and juvenile idiopathic arthritis in a European population. J Inflammation Res (2022) 15:977–86. doi: 10.2147/JIR.S354619
Keywords: osteoporosis, iron status, causal association, mendelian randomization, genetics
Citation: Xu J, Ma J, Chen J, Zhang S, Zheng C, Si H, Wu Y, Liu Y, Li M, Wu L and Shen B (2022) No genetic causal association between iron status and osteoporosis: A two-sample Mendelian randomization. Front. Endocrinol. 13:996244. doi: 10.3389/fendo.2022.996244
Received: 17 July 2022; Accepted: 21 November 2022;
Published: 09 December 2022.
Edited by:
Zhongjian Xie, Department of Metabolism and Endocrinology, Central South University, ChinaCopyright © 2022 Xu, Ma, Chen, Zhang, Zheng, Si, Wu, Liu, Li, Wu and Shen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bin Shen, shenbin_1971@163.com