Sha-Sha Tao1,2,3†
Peng Wang4†
Xin-Yi Wang5,6†
Kang-Jia Yin1,2
Xiao-Ke Yang7Zhi-Xin Wang7De-Guang Wang8*
Hai-Feng Pan1,2*- 1Department of Epidemiology and Biostatistics, School of Public Health, Anhui Medical University, Hefei, Anhui, China
- 2Inflammation and Immune Mediated Diseases Laboratory of Anhui Province, Hefei, Anhui, China
- 3Experimental Teaching Center for Preventive Medicine, School of Public Health, Anhui Medical University, Hefei, Anhui, China
- 4Teaching Center for Preventive Medicine, School of Public Health, Anhui Medical University, Hefei, Anhui, China
- 5Department of Radiation Oncology, The First Affiliated Hospital of Anhui Medical University, Hefei, Anhui, China
- 6The First Clinical Medical College, Anhui Medical University, Hefei, Anhui, China
- 7Department of Rheumatology and Immunology, The First Affiliated Hospital of Anhui Medical University, Hefei, China
- 8Department of Nephrology, The Second Affiliated Hospital of Anhui Medical University, Hefei, Anhui, China
Background: Polyunsaturated fatty acids (PUFAs) are closely related to osteoporosis. To test their causal relationship, we conducted a Mendelian randomization (MR) analysis.
Methods: We analyzed the causal relationship between four PUFAs measures, n-3 PUFAs (n-3), n-6 PUFAs (n-6), the ratio of n-3 PUFAs to total fatty acids (n-3 pct), and the ratio of n-6 PUFAs to n-3 PUFAs (n-6 to n-3), and five measures of osteoporosis, including estimated bone mineral density (eBMD), forearm (FA) BMD, femoral neck (FN) BMD, lumbar spine (LS) BMD, and fracture, using two-sample MR analysis. In order to verify the direct effect between PUFAs and BMD, we chose interleukin-6 (IL-6), tumor necrosis factor-β (TNF-β), and bone morphogenetic proteins 7 (BMP-7), three markers or cytokines strongly related to BMD, as possible confounding factors, and analyzed the possible causal relationships between them and PUFAs or BMD by MR. Inverse variance weighting (IVW), MR-Egger, weighted and weighted median were conducted. MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) and MR-Egger regression methods were used to evaluate the potential pleiotropy of instrumental variables (IVs) and outliers were identified by MR-PRESSO. Cochran’s Q statistic was used to detect the heterogeneity among IVs. Leave-one-out sensitivity analysis was used to find SNPs that have a significant impact on the results. All results were corrected by the Bonferroni correction.
Results: The IVW results showed that n-3 PUFAs (OR = 1.030, 95% CI: 1.013, 1.047, P = 0.001) and n-6 PUFAs (OR = 1.053, 95% CI: 1.034, 1.072, P < 0.001) were positively correlated with eBMD, while n-6 to n-3 (OR = 0.947, 95% CI: 0.924, 0.970, P < 0.001) were negatively correlated with eBMD. These casual relationships still existed after Bonferroni correction. There were positive effects of n-3 PUFAs on FA BMD (OR = 1.090, 95% CI: 1.011, 1.176, P = 0.025) and LS BMD (OR = 1.056, 95% CI: 1.011, 1.104, P = 0.014), n-3 pct on eBMD (OR = 1.028, 95% CI: 1.002, 1.055, P = 0.035) and FA BMD (OR = 1.090, 95% CI: 1.011, 1.174, P = 0.025), n-6 to n-3 on LS BMD (OR = 1.071, 95% CI: 1.021, 1.124, P = 0.005); negative effects of n-3 pct on fracture (OR = 0.953, 95% CI: 0.918, 0.988, P = 0.009) and n-6 to n-3 on FA BMD (OR = 0.910, 95% CI: 0.837, 0.988, P = 0.025). However, these causal effects all disappeared after Bonferroni correction (all P > 0.0025). None of IL-6, TNF-β, and BMP-7 had a causal effect on PUFA and BMD simultaneously (all P > 0.05).
Conclusion: Evidence from this MR study supports the genetically predicted causal effects of n-3, n-6, n-3 pct, and n-6 to n-3 on eBMD. In addition, n-3 not only associate with FA BMD and LS BMD through its own level and n-6 to n-3, but also link to fracture through n-3 pct.
Introduction
Osteoporosis is one of the most common systemic bone diseases in clinic, which is mainly manifested with the degradation of bone tissue microenvironment, the reduction of bone mineral density (BMD), the increase of bone fragility, and the high risk of fracture (1, 2). Osteoporosis mostly occurs in the postmenopausal women. Worldwide, fractures caused by osteoporosis can be as high as 9 million times a year (3). With the increased life expectancy of the global population and the improvement of medical and health conditions, the prevention and treatment of osteoporosis and fractures is still a common public health challenge and health care problem in the world today (4–7). It has been suggested that 20–40% of the risk of osteoporosis is caused by environmental factors such as nutrition (8).
Polyunsaturated fatty acids (PUFAs) are important immune nutrients, which play a role in nutritional treatment of a variety of diseases, including cancers, inflammatory diseases and osteoporosis (9–20). PUFAs are straight chain fatty acids with two or more double bonds, and the length of carbon chain is 18–22 carbon atoms. In PUFAs, the first unsaturated double bond is located between the 3rd and 4th carbon atoms starting from the methyl end, which is called omega-3 (n-3) PUFAs, and between the 6th and 7th carbon atoms, which is called omega-6 (n-6) PUFAs (15). Mounting evidence has suggested that PUFAs are involved in osteolysis, bone formation, bone development, bone metabolism, and metabolic bone diseases including osteoporosis and may be beneficial for skeletal health (21–29).
Recently, the relationship between PUFAs and osteoporosis has attracted a lot of attention. PUFAs intake and the ratio of n-6 to n-3 PUFAs are reported to be associated with BMD in humans (27, 30–32). Some studies found that total PUFAs, n-3 PUFAs, n-6 PUFAs intake can increase BMD, decrease the risk of fractures and beneficial for osteoporosis (30–34). A low n-6 to n-3 PUFAs ratio was proposed to be beneficial for the bone quality of rats (35). However, conclusions about the relationship between n-3 PUFAs, n-6 PUFAs, and osteoporosis are inconsistent. It was also observed that the intake of n-3 PUFAs or the ratio of n-6 to n-3 PUFAs are not associated with osteoporotic fractures, while the intake of n-6 PUFAs is positively associated with an elevated risk of fracture (32). Other study also indicated that increased intake of PUFAs is associated with greater perimenopausal femoral neck (FN) BMD loss (36). In postmenopausal women, the lower intake of marine n-3 PUFAs and the higher intake of n-6 PUFAs were observed that can decrease the risk of total fracture (33).
Clarifying the relationship between various PUFAs levels and BMD can not only further study the prevention, health care and treatment strategies of osteoporosis patients, but also can estimate the fracture risk of osteoporosis patients. It is imperative to evaluate the causal relationship between PUFAs level and BMD. However, from the above evidence, conclusions from observational studies and randomized controlled trials (RCT) about the relationship between PUFAs and osteoporosis are inconsistent, as well as the causal relationship between PUFAs and osteoporosis remains obscure.
Genome wide association study (GWAS) shows that BMD is a trait controlled by multiple genes, and is easy to be affected by environmental factors and has the tendency of family aggregation. Some genetic determinants of fractures are regulated by lower BMD (37–40). Mendelian randomization (MR) is an advanced research using the law of free combination, which uses genetic variation as instrumental variable (IVs) (41). The free combination of alleles can reduce the confounding effect of environmental factors; the one-way relationship that genes affect traits but diseases do not change genotypes can effectively avoid reverse causal bias (42). Compared with RCT, MR can extract useful information from the existing GWAS database on a large scale, with stronger statistical ability and wider coverage (43).
Materials and methods
Study design
In this study, we used two-sample MR studies to assess the causal relationships between PUFAs and osteoporosis. We selected four measures of PUFAs, including the circulating level of n-3 PUFAs (n-3), the circulating level of n-6 PUFAs (n-6), the ratio of n-3 fatty acids to total fatty acids (n-3 pct), and the ratio of n-6 PUFAs to n-3 PUFAs (n-6 to n-3). For osteoporosis, estimated BMD (eBMD), forearm (FA) BMD, femoral neck (FN) BMD, lumbar (LS) BMD, and fracture were selected as measures.
In order to verify the direct effect between PUFAs and BMD, we chose interleukin-6 (IL-6) (44), tumor necrosis factor-β (TNF-β) (45) and bone morphogenetic proteins 7 (BMP-7) (46), three markers or cytokines strongly related to BMD, as possible confounding factors, and analyzed the possible causal relationships between them and PUFAs or BMD by MR.
Data sources and single-nucleotide polymorphism selection
Genome wide association study of polyunsaturated fatty acids
The GWAS data on PUFAs came from the latest and largest public GWAS analysis by Borges MC et al. (47). In the study, the GWAS data of n-3, n-6, n-3 pct, and n-6 to n-3 were generated from 114,999 UK Biobank participants of European ancestry using BOLT-LMM (v2.3) (Supplementary Table 1) (48).
Genome wide association studyof osteoporosis
Genetic associations of eBMD and fracture were obtained from the largest public GWAS of Morris et al. (1), in which eBMD GWAS data were retrieved from 426,824 individuals of European ancestry in the UK Biobank and fracture GWAS data were retrieved from 416,795 UK Biobank European ancestry participants (ncases = 53,184 and ncontrols = 373,611) (Supplementary Table 1).
The research by Zheng et al. (49) provided GWAS data about FA BMD, FN BMD, and LS BMD. In the research, the GWAS data for the three BMDs were obtained from 10,805, 49,988, and 44,731 individuals of European ancestry, respectively (Supplementary Table 1).
Genome wide association study of possible confounding factors
Genetic associations of IL-6 and TNF-β were obtained from the public GWAS of Ahola-Olli et al. (50), and the GWAS data for both were retrieved from 8,293 individuals of European ancestry. The GWAS data of BMP-7 was obtained from 3,301 individuals of European ancestry (51) (Supplementary Table 1).
Single-nucleotide polymorphism selection
We screened single-nucleotide polymorphisms (SNPs) strongly related to exposure factors (P < 5 × 10–8) from the exposure GWAS, used clustering process (R2 < 0.001 and clumping distance = 10,000 kb) to eliminate linkage disequilibrium (LD) between SNPs, and excluded SNPs with minor allele frequency (MAF < 0.01), which ensured that the end result was undisturbed and feasible. The selected SNPs were matched with the outcome GWAS, and if SNP can not be found in outcome GWAS, its proxy SNP with high LD (r2 > 0.8) was used instead. Finally, other SNPs were selected as IVs after the palindrome SNPs were removed.
Statistical analyses
We used four complementary and mutually corroborative methods to analyze the causal relationship between PUFAs and osteoporosis, including inverse variance weighting (IVW), MR-Egger, weighted and weighted median, among which IVW is the main analysis method. The weighted median estimator serves as an unbiased causal effect estimate when up to 50% of the instruments are invalid, by estimating the causal effect as the median of the weighted ratio estimates (52). At the same time, we used MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) and MR-Egger regression methods to evaluate the potential level pleiotropy of IVs (53, 54). Meanwhile, MR-PRESSO can also find abnormal values in IVs. After removing the abnormal values, MR-PRESSO and MR-Egger tests were performed again until there was no horizontal pleiotropic SNP in all IVs. Then, we applied Cochran’s Q statistic to detect and quantify the heterogeneity among IVs (55). Leave-one-out sensitivity analysis was used to find and eliminate SNPs that have a significant impact on the results, so as to ensure the robustness of causal relationship estimation. There were four exposures (n-3, n-6, n-3pct, and n-6 to n-3) and five outcomes (eBMD, FA, FN, LS, and fracture) in this study, therefore the Bonferroni method was conducted to correct for multiple comparisons and the P-value was less than 0.0025 (0.05 was divided by 4 × 5) (56, 57). All statistical analyses were performed using the packages “TwoSampleMR” and “MRPRESSO” in R version 4.1.1.
Results
Instrumental variables selection
Instrumental variables of polyunsaturated fatty acids selection
After clumping process, there were 52 SNPs, 63 SNPs, 41 SNPs, and 41 SNPs strongly associated (P < 5 × 10–8) with the circulating level of n-3 PUFAs, the circulating level of n-6 PUFAs, n-3 pct, n-6 to n-3, respectively and no LD was screened out. MAF of all above SNPs were not less than 0.01. The main information of SNPs is listed in Supplementary Tables 2–5.
Instrumental variables of possible confounding factors selection
After clumping process, there were 2 SNPs and 2 SNPs strongly associated (P < 5 × 10–8) with TNF-β and BMP-7, respectively. There was no SNPs available at the genome wide significance threshold (P < 5 × 10–8) of IL-6, so we relaxed the significance threshold to P < 5 × 10–7, and 2 SNPs were strongly associated (P < 5 × 10–7) with IL-6. Among the above selected SNPs, no LD was found, and MAF of these SNPs were not less than 0.01. The main information of SNPs is provided in Supplementary Tables 6, 7.
Causal relationship between polyunsaturated fatty acids and osteoporosis
n-3 polyunsaturated fatty acids on osteoporosis
From the IVW results, n-3 PUFAs had positive effects on eBMD (OR = 1.030, 95% CI: 1.013, 1.047, P = 0.001), FA BMD (OR = 1.090, 95% CI: 1.011, 1.176, P = 0.025), and LS BMD (OR = 1.056, 95% CI: 1.011, 1.104, P = 0.014), however, after further Bonferroni correction, only the effect on eBMD remained. No causal effects of n-3 PUFAs on FN BMD (OR = 0.988, 95% CI: 0.950, 1.028, P = 0.557) and fracture (OR = 0.982, 95% CI: 0.938, 1.028, P = 0.433) were observed (Figure 1 and Supplementary Table 8).
Figure 1. Mendelian randomization (MR) estimate results of polyunsaturated fatty acids (PUFAs) on outcomes. n-3 pct, the ratio of n-3 fatty acids to total fatty acids; n-6 to n-3, the ratio of n-6 PUFAs to n-3 PUFAs; BMD, bone mineral density; eBMD, estimated BMD; FA, forearm; FN, femoral neck; LS, lumbar.
n-6 polyunsaturated fatty acids on osteoporosis
The results of IVW showed a positive effect of n-6 PUFAs on eBMD (OR = 1.053, 95% CI: 1.034, 1.072, P < 0.001), which still persisted after Bonferroni correction. However, no causal associations were found between n-6 PUFAs and FA BMD (OR = 1.022, 95% CI: 0.925, 1.128, P = 0.673), FN BMD (OR = 1.015, 95% CI: 0.962, 1.071, P = 0.584), LS BMD (OR = 1.013, 95% CI: 0.951, 1.078, P = 0.688), and fracture (OR = 1.014, 95% CI: 0.970, 1.059, P = 0.545) (Figure 1 and Supplementary Table 9).
n-3 pct on osteoporosis
The results of IVW revealed that n-3 pct had positive causal relationships with eBMD (OR = 1.028, 95% CI: 1.002, 1.055, P = 0.035) and FA BMD (OR = 1.090, 95% CI: 1.011, 1.174, P = 0.025) as well as a negative causal relationship with fracture (OR = 0.953, 95% CI: 0.918, 0.988, P = 0.009). However, these causal effects all disappeared after Bonferroni correction. There were no causal effects of n-3 pct on FN BMD (OR = 0.992, 95% CI: 0.951, 1.036, P = 0.732) and LS BMD (OR = 1.025, 95% CI: 0.969, 1.085, P = 0.384) (Figure 1 and Supplementary Table 10).
n-6 to n-3 on osteoporosis
The results of IVW indicated that n-6 to n-3 had a negative effect on eBMD (OR = 0.947, 95% CI: 0.924, 0.970, P < 0.001), which remain persisted after Bonferroni correction and had a negative relationship with FA BMD (OR = 0.910, 95% CI: 0.837, 0.988, P = 0.025) while a positive relationship with LS BMD (OR = 1.071, 95% CI: 1.021, 1.124, P = 0.005), which both disappeared after Bonferroni correction. No causal effects of n-6 to n-3 on FN BMD (OR = 1.002, 95% CI: 0.959, 1.046, P = 0.945) and fracture (OR = 1.035, 95% CI: 0.979, 1.095, P = 0.225) were observed (Figure 1 and Supplementary Table 11).
Causal relationship between possible confounding factors and polyunsaturated fatty acids
From the IVW results, BMP-7 had a negative effect on n-3 pct (OR = 0.967, 95% CI: 0.937, 0.998, P = 0.038) and a positive effect on n-6 to n-3 (OR = 1.035, 95% CI: 1.002, 1.069, P = 0.036). However, no causal effect was observed of BMP-7 on n-3 PUFAs and n-6 PUFAs (all P > 0.05). Moreover, no causal relationship was found between IL-6, TNF-β and all four outcomes (all P > 0.05) (Figure 2 and Supplementary Table 12).
Figure 2. Mendelian randomization (MR) estimate results of possible confounding factors on polyunsaturated fatty acids (PUFAs). n-3 pct, the ratio of n-3 fatty acids to total fatty acids; n-6 to n-3, the ratio of n-6 PUFAs to n-3 PUFAs; IL, interleukin; TNF, tumor necrosis factor; BMP, bone morphogenetic protein.
Causal relationship between possible confounding factors and osteoporosis
From the IVW results, no causal relationship was found between IL-6, TNF-β, BMP-7 and eBMD, FA, FN, LS as well as fracture (all P > 0.05) (Figure 3 and Supplementary Table 13).
Figure 3. Mendelian randomization (MR) estimate results of possible confounding factors on BMD. BMD, bone mineral density; eBMD, estimated BMD; FA, forearm; FN, femoral neck; LS, lumbar; IL, interleukin; TNF, tumor necrosis factor; BMP, bone morphogenetic protein.
Pleiotropy and sensitivity analysis
The heterogeneity test did not find any heterogeneity between selected IVs of n-3 PUFAs and n-6 PUFAs. Heterogeneity was not observed between the IVs of n-3 pct on eBMD and on FA BMD, while was found among IVs of n-3 pct on FN BMD, on fractrue and on LS BMD. No heterogeneities were found between IVs of n-6 to n-3, except for IVs of n-6 to n-3 on FN BMD (Q = 50.008, P = 0.022). No heterogeneity was found between the IVs of IL-6, TNF-β and BMP-7.
MR Pleiotropy RESidual Sum and Outlier global test defined SNPs with horizontal pleiotropy as outliers which were listed in Supplementary Tables 2–7. After removing the outliers, MR-Egger regression and MR-PRESSO global test were used to verify that there was no horizontal pleiotropy between IVs and results. Leave-one-out analysis suggested that the outcomes were not caused by any SNPs. Supplementary Tables 8–13 and Supplementary Figures 1–16 show the results of pleiotropy and sensitivity analysis.
Discussion
As the IVW results shown, after Bonferroni correction n-3 PUFAs and n-6 PUFAs were still positively correlated with eBMD, while n-6 to n-3 were negatively correlated with eBMD, which provides new evidence to support the relationship between PUFAs and osteoporosis.
In this study, no causal effect was observed of IL-6 and TNF-β on four measures of PUFAs and five measures of osteoporosis, suggesting that the causal effect of four measures of PUFAs on BMD was not affected by IL-6 and TNF-β, moreover IL-6 and TNF-β could not affect BMD by affecting PUFAs. As we all know, only in situations where the total effect, direct effect and indirect effect all act in the same direction, the exposure having indirect effect can be identified as a mediator (58). BMP-7 has a causal effect on-3 pct and n-6 to n-3, but has no causal effect on BMD, suggesting that BMP-7 can not affect BMD through PUFAs and is also not a mediator in the causal effect of PUFAs on BMD.
Recently, the relationship between PUFAs and osteoporosis has attracted a lot of attention. Some studies have found that intake of total PUFA, n-3 PUFAs, n-6 PUFAs can increase BMD and be beneficial to osteoporosis (30–34), which is mutually confirmed with our results.
These may attribute to the profitable role of PUFAs in bone formation, absorption, development and metabolism, and n-3 PUFAs can also regulate bone health by increasing osteoblast activity and decreasing osteoclast activity, promoting intestinal calcium absorption and mineral deposition during bone development (25). The intake of n-3 PUFAs has also been observed to be associated with increased bone regeneration, improved bone microstructure and strength (59–61). Osteoporosis increases the apoptosis, adipogenic differentiation, and levels of RANKL and sclerostin of bone marrow mesenchymal stem cells and osteoblasts (62). Bone mineral loss is the result of an imbalance between osteoblastic bone formation and osteoclast bone resorption. As an important n-3 PUFAs, DHA is a lipid component specific to the osteoblast membrane, which induces extensive lipid remodeling in mesenchymal stem cells, resulting in more stable membrane microdomains and thus enhanced osteogenic differentiation (63). It has been reported that dietary n-3 can reduce osteoclast formation and bone loss in ovariectomized mice (64). In rats, taking fish oil can also inhibit alveolar bone absorption and osteoclast activity (65).
In this study, n-6 to n-3 was observed to be negatively related to eBMD, which was consistent with previous studies. In previous study, a low n-6 to n-3 PUFAs ratio has also been proposed to be beneficial for the bone quality of rats (35). It has been widely documented that reducing the n-6 to n-3 PUFAs ratio can prevent bone mineral loss and prostaglandin(PG) E2 production in animal and in vitro cell culture experiments (66). Kelly et al. (67) also proposed that the high proportion of n-6 to n-3 PUFAs may be one of the important reasons for the increased risk of obesity and osteoporosis. According to the IVW results, the relationship between n-3 PUFAs, n-6 PUFAs, or the ratio of n-6 to n-3 PUFAs and fractures was not observed, while n-3 pct was negatively related to fractures, which suggests that it may not be the level of PUFAs who affected the fracture, but n-3 pct. These results, together with the n-6 to n-3 results, suggest that n-3 may influence BMD not only by itself or by its ratio to n-6, but also by its ratio to total fatty acids. The single item has little significance, but its proportion has important diagnostic value, which is not uncommon in clinical practice. For example, neutrophils to lymphocytes ratio (NLR) is an inflammatory index, which plays a role in the prognosis evaluation of sepsis and various diseases. A higher NLR may indicate more serious infection and worse prognosis (68, 69). The ratio of soluble fms-like tyrosine kinase 1 (sFlt-1) to placental growth factor (PlGF) can be used as a monitoring indicator of preeclampsia (PE), which is a kind of hypertensive disorder complicating pregnancy and seriously endangers the health of mothers and infants, with unpredictable outcomes (70).
Moreover, the IVW results showed that FA BMD and LS BMD, which were less affected by weight bearing, were more susceptible to n-3 and n-6 to n-3 than FN BMD which was more affected by weight bearing. Fat is digested, decomposed and metabolized into glycerin and PUFAs (71). Like vitamins and minerals, PUFAs are closely related to bone health through various ways. The distribution and accumulation of adipose tissue is extremely important for bone health. At the same time, adipose tissue secretions such as leptin, adiponectin, estrogen and osteocalcin can also act on bones. Multiple studies have shown that BMI is positively correlated with BMD, with lower BMI being thought to increase the risk of osteoporosis, while higher body weight (even obesity) protects bones (72–75). However, in recent years, this “obesity paradox” has been challenged like never before. Fat-rich bone marrow may be the cause of osteoporosis, especially in postmenopausal women (76). In 2011, the UK Fracture Liaison Service first reported that the incidence of obesity in postmenopausal women with fragility fractures can be as high as 27% (77). Two other studies were also pointed out that visceral fat was significantly associated with bone loss (78, 79). Excessive body fat, especially abdominal fat, produces inflammatory cytokines that stimulate increased bone marrow lipogenesis, increased bone resorption, decreased bone strength, and decreased bone mass (80).
Advantages and disadvantages
This study has several advantages. First of all, as a MR study, this study investigated the causal associations between different types of PUFAs and BMD in different parts, with detailed classification and comprehensive research. Secondly, we not only studied the level of PUFAs, but also the proportion of PUFAs, and all the results verified each other, and the conclusions were unified. Third, we used a strict Bonferroni correction, thus the conclusions are robust.
At the same time, this study also has some limitations. First, all GWAS data in this study were from European populations, and the representativeness of the results to the entire population remains to be determined. Second, the relationship between different doses of PUFAs and BMD has not been studied, and more detailed quantitative experiments are needed. In addition, the mediating effects of obesity, BMD and other factors still need further research.
Conclusion
This MR study establishes that n-3, n-6, n-3 pct, and n-6 to n-3 are causally associated with eBMD. In addition, n-3 not only associated with FA BMD and LS BMD through its own level and n-6 to n-3, but also associated with fracture through n-3 pct. Our findings provide new clue to further reveal the pathogenic role and therapeutic potential of PUFAs in osteoporosis.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
H-FP and D-GW conceived the present idea and were responsible for the design of the study. S-ST performed the statistical analysis and manuscript writing. PW and X-YW participated in acquisition of data and data analysis. K-JY, X-KY, and Z-XW participated in acquisition of data. All authors have read and approved the submitted version for publication.
Funding
This study was funded by grants from the National Natural Science Foundation of China (81872687), Anhui Provincial University Natural Science Research Project (KJ2021A0230), and Research Fund of Anhui Institute of Translational Medicine (2021zhyx-B04).
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/fnut.2022.1014847/full#supplementary-material
References
1. Morris JA, Kemp JP, Youlten SE, Laurent L, Logan JG, Chai RC, et al. An atlas of genetic influences on osteoporosis in humans and mice. Nat Genet. (2019) 51:258–66. doi: 10.1038/s41588-018-0302-x
2. Liu J, Curtis EM, Cooper C, Harvey NC. State of the art in osteoporosis risk assessment and treatment. J Endocrinol Invest. (2019) 42:1149–64. doi: 10.1007/s40618-019-01041-6
3. Johnell O, Kanis JA. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos Int. (2006) 17:1726–33.
4. Reid IR. A broader strategy for osteoporosis interventions. Nat Rev Endocrinol. (2020) 16:333–9. doi: 10.1038/s41574-020-0339-7
5. Khosla S, Hofbauer LC. Osteoporosis treatment: recent developments and ongoing challenges. Lancet Diabetes Endocrinol. (2017) 5:898–907. doi: 10.1016/S2213-8587(17)30188-2
6. Compston JE, McClung MR, Leslie WD. Osteoporosis. Lancet. (2019) 393:364–76. doi: 10.1016/S0140-6736(18)32112-3
7. Hurley DL, Khosla S. Update on primary osteoporosis. Mayo Clin Proc. (1997) 72:943–9. doi: 10.4065/72.10.943
8. Weaver CM, Gordon CM, Janz KF, Kalkwarf HJ, Lappe JM, Lewis R, et al. The National Osteoporosis Foundation’s position statement on peak bone mass development and lifestyle factors: a systematic review and implementation recommendations. Osteoporos Int. (2016) 27:1281–386. doi: 10.1007/s00198-015-3440-3
9. Vega OM, Abkenari S, Tong Z, Tedman A, Huerta-Yepez S. Omega-3 polyunsaturated fatty acids and lung cancer: nutrition or pharmacology? Nutr Cancer. (2021) 73:541–61. doi: 10.1080/01635581.2020.1761408
10. Nasir A, Bullo MMH, Ahmed Z, Imtiaz A, Yaqoob E, Jadoon M, et al. Nutrigenomics: epigenetics and cancer prevention: a comprehensive review. Crit Rev Food Sci Nutr. (2020) 60:1375–87. doi: 10.1080/10408398.2019.1571480
11. Ghamarzad Shishavan, N, Masoudi S, Mohamadkhani A, Sepanlou SG, Poustchi H, Hekmatdoost A, et al. The association of dietary intake and plasma fatty acid panel in pancreatic cancer patients: results from golestan cohort study. Nutr Health. (2022) 11:2601060221114712. doi: 10.1177/02601060221114712
12. Wawrzyniak P, Noureddine N, Wawrzyniak M, Lucchinetti E, Krämer SD, Rogler G, et al. Nutritional lipids and mucosal inflammation. Mol Nutr Food Res. (2021) 65:e1901269. doi: 10.1002/mnfr.201901269
13. Engelen M, Jonker R, Sulaiman H, Fisk HL, Calder PC, Deutz NEP. ω-3 polyunsaturated fatty acid supplementation improves postabsorptive and prandial protein metabolism in patients with chronic obstructive pulmonary disease: a randomized clinical trial. Am J Clin Nutr. (2022) 116:686–98. doi: 10.1093/ajcn/nqac138
14. Jiang H, Wang L, Wang D, Yan N, Li C, Wu M, et al. Omega-3 polyunsaturated fatty acid biomarkers and risk of type 2 diabetes, cardiovascular disease, cancer, and mortality. Clin Nutr. (2022) 41:1798–807. doi: 10.1016/j.clnu.2022.06.034
15. Wall R, Ross RP, Fitzgerald GF, Stanton C. Fatty acids from fish: the anti-inflammatory potential of long-chain omega-3 fatty acids. Nutr Rev. (2010) 68:280–9. doi: 10.1111/j.1753-4887.2010.00287.x
16. Fritsche K. Fatty acids as modulators of the immune response. Annu Rev Nutr. (2006) 26:45–73. doi: 10.1146/annurev.nutr.25.050304.092610
17. Kang JX, Liu A. The role of the tissue omega-6/omega-3 fatty acid ratio in regulating tumor angiogenesis. Cancer Metastasis Rev. (2013) 32:201–10. doi: 10.1007/s10555-012-9401-9
18. Boston PF, Bennett A, Horrobin DF, Bennett CN. Ethyl-EPA in Alzheimer’s disease–a pilot study. Prostaglandins Leukot Essent Fatty Acids. (2004) 71:341–6. doi: 10.1016/j.plefa.2004.07.001
19. Zhang W, Wang H, Zhang H, Leak RK, Shi Y, Hu X, et al. Dietary supplementation with omega-3 polyunsaturated fatty acids robustly promotes neurovascular restorative dynamics and improves neurological functions after stroke. Exp Neurol. (2015) 272:170–80. doi: 10.1016/j.expneurol.2015.03.005
20. Orchard TS, Pan X, Cheek F, Ing SW, Jackson RD. A systematic review of omega-3 fatty acids and osteoporosis. Br J Nutr. (2012) 107(Suppl. 2):S253–60. doi: 10.1017/S0007114512001638
21. Watkins BA, Shen CL, Allen KG, Seifert MF. Dietary (n-3) and (n-6) polyunsaturates and acetylsalicylic acid alter ex vivo PGE2 biosynthesis, tissue IGF-I levels, and bone morphometry in chicks. J Bone Miner Res. (1996) 11:1321–32. doi: 10.1002/jbmr.5650110917
22. Watkins BA, Lippman HE, Le Bouteiller L, Li Y, Seifert MF. Bioactive fatty acids: role in bone biology and bone cell function. Prog Lipid Res. (2001) 40:125–48. doi: 10.1016/S0163-7827(00)00016-3
23. Sakaguchi K, Morita I, Murota S. Eicosapentaenoic acid inhibits bone loss due to ovariectomy in rats. Prostaglandins Leukot Essent Fatty Acids. (1994) 50:81–4. doi: 10.1016/0952-3278(94)90151-1
24. Li Y, Seifert MF, Ney DM, Grahn M, Grant AL, Allen KG, et al. Dietary conjugated linoleic acids alter serum IGF-I and IGF binding protein concentrations and reduce bone formation in rats fed (n-6) or (n-3) fatty acids. J Bone Miner Res. (1999) 14:1153–62. doi: 10.1359/jbmr.1999.14.7.1153
25. Lau BY, Cohen DJ, Ward WE, Ma DW. Investigating the role of polyunsaturated fatty acids in bone development using animal models. Molecules. (2013) 18:14203–27. doi: 10.3390/molecules181114203
26. Hutchins-Wiese HL, Picho K, Watkins BA, Li Y, Tannenbaum S, Claffey K, et al. High-dose eicosapentaenoic acid and docosahexaenoic acid supplementation reduces bone resorption in postmenopausal breast cancer survivors on aromatase inhibitors: a pilot study. Nutr Cancer. (2014) 66:68–76. doi: 10.1080/01635581.2014.847964
27. Weiss LA, Barrett-Connor E, von Mühlen D. Ratio of n-6 to n-3 fatty acids and bone mineral density in older adults: the rancho bernardo study. Am J Clin Nutr. (2005) 81:934–8. doi: 10.1093/ajcn/81.4.934
28. Griel AE, Kris-Etherton PM, Hilpert KF, Zhao G, West SG, Corwin RL. An increase in dietary n-3 fatty acids decreases a marker of bone resorption in humans. Nutr J. (2007) 6:2. doi: 10.1186/1475-2891-6-2
29. Claassen N, Coetzer H, Steinmann CM, Kruger MC. The effect of different n-6/n-3 essential fatty acid ratios on calcium balance and bone in rats. Prostaglandins Leukot Essent Fatty Acids. (1995) 53:13–9. doi: 10.1016/0952-3278(95)90078-0
30. Coetzee M, Haag M, Kruger MC. Effects of arachidonic acid, docosahexaenoic acid, prostaglandin E(2) and parathyroid hormone on osteoprotegerin and RANKL secretion by MC3T3-E1 osteoblast-like cells. J Nutr Biochem. (2007) 18:54–63. doi: 10.1016/j.jnutbio.2006.03.002
31. Högström M, Nordström P, Nordström A. N-3 Fatty acids are positively associated with peak bone mineral density and bone accrual in healthy men: the NO2 Study. Am J Clin Nutr. (2007) 85:803–7. doi: 10.1093/ajcn/85.3.803
32. Martínez-Ramírez MJ, Palma S, Martínez-González MA, Delgado-Martínez AD, de la Fuente C, Delgado-Rodríguez M. Dietary fat intake and the risk of osteoporotic fractures in the elderly. Eur J Clin Nutr. (2007) 61:1114–20. doi: 10.1038/sj.ejcn.1602624
33. Orchard TS, Cauley JA, Frank GC, Neuhouser ML, Robinson JG, Snetselaar L, et al. Fatty acid consumption and risk of fracture in the Women’s Health Initiative. Am J Clin Nutr. (2010) 92:1452–60. doi: 10.3945/ajcn.2010.29955
34. Dou Y, Wang Y, Chen Z, Yu X, Ma D. Effect of n-3 polyunsaturated fatty acid on bone health: a systematic review and meta-analysis of randomized controlled trials. Food Sci Nutr. (2022) 10:145–54. doi: 10.1002/fsn3.2655
35. Yang L, Yang C, Chu C, Wan M, Xu D, Pan D, et al. Beneficial effects of monounsaturated fatty acid-rich blended oils with an appropriate polyunsaturated/saturated fatty acid ratio and a low n-6/n-3 fatty acid ratio on the health of rats. J Sci Food Agric. (2022) 102:7172–85. doi: 10.1002/jsfa.12083
36. Macdonald HM, New SA, Golden MH, Campbell MK, Reid DM. Nutritional associations with bone loss during the menopausal transition: evidence of a beneficial effect of calcium, alcohol, and fruit and vegetable nutrients and of a detrimental effect of fatty acids. Am J Clin Nutr. (2004) 79:155–65. doi: 10.1093/ajcn/79.1.155
37. Richards JB, Zheng HF, Spector TD. Genetics of osteoporosis from genome-wide association studies: advances and challenges. Nat Rev Genet. (2012) 13:576–88. doi: 10.1038/nrg3228
38. Trajanoska K, Rivadeneira F. The genetic architecture of osteoporosis and fracture risk. Bone. (2019) 126:2–10. doi: 10.1016/j.bone.2019.04.005
39. Yang TL, Shen H, Liu A, Dong SS, Zhang L, Deng FY, et al. A road map for understanding molecular and genetic determinants of osteoporosis. Nat Rev Endocrinol. (2020) 16:91–103. doi: 10.1038/s41574-019-0282-7
40. Song YQ, Hu SD, Lin X, Meng XH, Wang X, Zhang YH, et al. Identification of PDXDC1 as a novel pleiotropic susceptibility locus shared between lumbar spine bone mineral density and birth weight. J Mol Med. (2022) 100:723–34. doi: 10.1007/s00109-021-02165-0
41. Smith GD, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. (2003) 32:1–22. doi: 10.1093/ije/dyg070
42. Zheng J, Baird D, Borges MC, Bowden J, Hemani G, Haycock P, et al. Recent developments in mendelian randomization studies. Curr Epidemiol Rep. (2017) 4:330–45. doi: 10.1007/s40471-017-0128-6
43. Sekula P, Del Greco MF, Pattaro C, Köttgen A. Mendelian randomization as an approach to assess causality using observational data. J Am Soc Nephrol. (2016) 27:3253–65. doi: 10.1681/ASN.2016010098
44. Krela-Kaźmierczak I, Szymczak-Tomczak A, Łykowska-Szuber L, Wysocka E, Michalak M, Stawczyk-Eder K, et al. Interleukin 6, osteoprotegerin, sRANKL and bone metabolism in inflammatory bowel diseases. Adv Clin Exp Med. (2018) 27:449–53. doi: 10.17219/acem/75675
45. Kim H, Chun S, Ku SY, Suh CS, Choi YM, Kim JG. Association between polymorphisms in tumor necrosis factor (TNF) and TNF receptor genes and circulating TNF, soluble TNF receptor levels, and bone mineral density in postmenopausal Korean women. Menopause. (2009) 16:1014–20. doi: 10.1097/gme.0b013e3181a039c8
46. Chen Y, Ma B, Wang X, Zha X, Sheng C, Yang P, et al. Potential functions of the BMP family in bone, obesity, and glucose metabolism. J Diabetes Res. (2021) 2021:6707464. doi: 10.1155/2021/6707464
47. Borges MC, Haycock PC, Zheng J, Hemani G, Holmes MV, Davey Smith G, et al. Role of circulating polyunsaturated fatty acids on cardiovascular diseases risk: analysis using Mendelian randomization and fatty acid genetic association data from over 114,000 UK Biobank participants. BMC Med. (2022) 20:210. doi: 10.1186/s12916-022-02399-w
48. Loh PR, Kichaev G, Gazal S, Schoech AP, Price AL. Mixed-model association for biobank-scale datasets. Nat Genet. (2018) 50:906–8. doi: 10.1038/s41588-018-0144-6
49. Zheng HF, Forgetta V, Hsu YH, Estrada K, Rosello-Diez A, Leo PJ, et al. Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture. Nature. (2015) 526:112–7.
50. Ahola-Olli AV, Wurtz P, Havulinna AS, Aalto K, Pitkänen N, Lehtimäki T, et al. Genome-wide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors. Am J Hum Genet. (2017) 100:40–50. doi: 10.1016/j.ajhg.2016.11.007
51. Sun BB, Maranville JC, Peters JE, Stacey D, Staley JR, Blackshaw J, et al. Genomic atlas of the human plasma proteome. Nature. (2018) 558:73–9. doi: 10.1038/s41586-018-0175-2
52. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
53. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
54. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
55. Greco MF, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. (2015) 34:2926–40. doi: 10.1002/sim.6522
56. Sedgwick P. Multiple hypothesis testing and Bonferroni’s correction. BMJ. (2014) 349:g6284. doi: 10.1136/bmj.g6284
57. Curtin F, Schulz P. Multiple correlations and Bonferroni’s correction. Biol Psychiatry. (1998) 44:775–7. doi: 10.1016/S0006-3223(98)00043-2
58. Carter AR, Sanderson E, Hammerton G, Richmond RC, Davey Smith G, Heron J, et al. Mendelian randomisation for mediation analysis: current methods and challenges for implementation. Eur J Epidemiol. (2021) 36:465–78. doi: 10.1007/s10654-021-00757-1
59. Martyniak K, Wei F, Ballesteros A, Meckmongkol T, Calder A, Gilbertson T, et al. Do polyunsaturated fatty acids protect against bone loss in our aging and osteoporotic population? Bone. (2021) 143:115736. doi: 10.1016/j.bone.2020.115736
60. Sharma T, Mandal CC. Omega-3 fatty acids in pathological calcification and bone health. J Food Biochem. (2020) 44:e13333. doi: 10.1111/jfbc.13333
61. Jørgensen HS, Eide IA, Jenssen T, Åsberg A, Bollerslev J, Godang K, et al. Marine n-3 polyunsaturated fatty acids and bone mineral density in kidney transplant recipients: a randomized, placebo-controlled trial. Nutrients. (2021) 13:2361. doi: 10.3390/nu13072361
62. 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:e45142. doi: 10.1371/journal.pone.0045142
63. Levental KR, Surma MA, Skinkle AD, Lorent JH, Zhou Y, Klose C, et al. ω-3 polyunsaturated fatty acids direct differentiation of the membrane phenotype in mesenchymal stem cells to potentiate osteogenesis. Sci Adv. (2017) 3:eaao1193. doi: 10.1126/sciadv.aao1193
64. Sun D, Krishnan A, Zaman K, Lawrence R, Bhattacharya A, Fernandes G. Dietary n-3 fatty acids decrease osteoclastogenesis and loss of bone mass in ovariectomized mice. J Bone Miner Res. (2003) 18:1206–16. doi: 10.1359/jbmr.2003.18.7.1206
65. Iwami-Morimoto Y, Yamaguchi K, Tanne K. Influence of dietary n-3 polyunsaturated fatty acid on experimental tooth movement in rats. Angle Orthod. (1999) 69:365–71.
66. Albertazzi P, Coupland K. Polyunsaturated fatty acids. Is there a role in postmenopausal osteoporosis prevention? Maturitas. (2002) 42:13–22. doi: 10.1016/s0378-5122(02)00022-1
67. Kelly OJ, Gilman JC, Kim Y, Ilich JZ. Long-chain polyunsaturated fatty acids may mutually benefit both obesity and osteoporosis. Nutr Res. (2013) 33:521–33. doi: 10.1016/j.nutres.2013.04.012
68. Wang K, Wang X, Du J, Liu C, Jiang Y, Zhang H, et al. Relationship between changes in the course of COVID-19 and ratio of neutrophils-to-lymphocytes and related parameters in patients with severe vs. common disease. Epidemiol Infect. (2021) 149:e81. doi: 10.1017/S0950268821000674
69. Templeton AJ, McNamara MG, Šeruga B, Aneja P, Ocaña A, Leibowitz-Amit R, et al. Prognostic role of neutrophil-to-lymphocyte ratio in solid tumors: a systematic review and meta-analysis. J Natl Cancer Inst. (2014) 106:dju124. doi: 10.1093/jnci/dju124
70. Zeisler H, Llurba E, Chantraine F, Vatish M, Staff AC, Sennström M, et al. Predictive value of the sFlt-1:PlGF ratio in women with suspected preeclampsia. N Engl J Med. (2016) 374:13–22. doi: 10.1056/NEJMoa1414838
71. de Carvalho C, Caramujo MJ. The various roles of fatty acids. Molecules. (2018) 23:2583. doi: 10.3390/molecules23102583
72. Poet JL, Galinier Pujol A, Tonolli Serabian I, Conte Devolx B, Roux H. Lumbar bone mineral density in anorexia nervosa. Clin Rheumatol. (1993) 12:236–9. doi: 10.1007/BF02231534
73. Asomaning K, Bertone-Johnson ER, Nasca PC, Hooven F, Pekow PS. The association between body mass index and osteoporosis in patients referred for a bone mineral density examination. J Womens Health. (2006) 15:1028–34. doi: 10.1089/jwh.2006.15.1028
74. Sukumar D, Schlussel Y, Riedt CS, Gordon C, Stahl T, Shapses SA. Obesity alters cortical and trabecular bone density and geometry in women. Osteoporos Int. (2011) 22:635–45. doi: 10.1007/s00198-010-1305-3
75. Turcotte AF, O’Connor S, Morin SN, Gibbs JC, Willie BM, Jean S, et al. Association between obesity and risk of fracture, bone mineral density and bone quality in adults: a systematic review and meta-analysis. PLoS One. (2021) 16:e0252487. doi: 10.1371/journal.pone.0252487
76. Menagh PJ, Turner RT, Jump DB, Wong CP, Lowry MB, Yakar S, et al. Growth hormone regulates the balance between bone formation and bone marrow adiposity. J Bone Miner Res. (2010) 25:757–68. doi: 10.1359/jbmr.091015
77. Compston J. Obesity and fractures in postmenopausal women. Curr Opin Rheumatol. (2015) 27:414–9. doi: 10.1097/BOR.0000000000000182
78. Jankowska EA, Rogucka E, Medraś M. Are general obesity and visceral adiposity in men linked to reduced bone mineral content resulting from normal ageing? A population-based study. Andrologia. (2001) 33:384–9. doi: 10.1046/j.1439-0272.2001.00469.x
79. Tseng YH, Huang KC, Liu ML, Shu WT, Sheu WH. Association between metabolic syndrome (MS) and bone mineral loss: a cross-sectional study in Puli Township in Taiwan. Arch Gerontol Geriatr. (2009) 49(Suppl. 2):S37–40. doi: 10.1016/S0167-4943(09)70011-1
Keywords: polyunsaturated fatty acids, osteoporosis, bone mineral density, genetic association, causal relationship, Mendelian randomization
Citation: Tao S-S, Wang P, Wang X-Y, Yin K-J, Yang X-K, Wang Z-X, Wang D-G and Pan H-F (2022) Causal effect of polyunsaturated fatty acids on bone mineral density and fracture. Front. Nutr. 9:1014847. doi: 10.3389/fnut.2022.1014847
Received: 09 August 2022; Accepted: 21 November 2022;
Published: 08 December 2022.
Edited by:
Sorayya Kheirouri, Tabriz University of Medical Sciences, IranReviewed by:
Hans Demmelmair, Ludwig Maximilian University of Munich, GermanyUndurti Narasimha Das, UND Life Sciences LLC, United States
Copyright © 2022 Tao, Wang, Wang, Yin, Yang, Wang, Wang and Pan. 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: De-Guang Wang, wangdeguang@ahmu.edu.cn; Hai-Feng Pan, panhaifeng1982@sina.com, panhaifeng@ahmu.edu.cn
†These authors have contributed equally to this work and share first authorship