Gaohua Liu1†
Wenkang Chen- 1Institute of Clinical Medicine, The First Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, Hunan, China
- 2Department of Hepatobiliary Surgery, The First Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, Hunan, China
- 3Department of Hematology, Fujian Medical University Union Hospital, Fuzhou, Fujian, China
- 4Speciality of Sports Medicine in Department of Orthopaedics, The First Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, Hunan, China
Introduction: The relationship between intervertebral disc degeneration (IVDD) and osteoporosis (OP), diagnosed primarily using bone mineral density (BMD), remains unclear so far. The present study, therefore, aimed to investigate the potential relationship between osteoporosis and intervertebral disc degeneration using Mendelian randomization and genome-wide association analyses. Specifically, the impact of bone mineral density on the development of intervertebral disc degeneration was evaluated.
Materials and methods: The genome-wide association studies (GWAS) summary data of OP/BMDs and IVDD were collected from the FinnGen consortium, the GEFOS consortium, and MRC-IEU. The relationship between IVDD and OP was then explored using TSMR. The inverse-variance weighted (IVW) method was adopted as the primary effect estimate, and the reliability and stability of the results were validated using various methods, including MR-Egger, weighted median, simple mode, weighted mode, and MR-PRESSO.
Results: No significant causal relationship was observed between OP and IVDD (IVW, P > 0.05) or between femoral neck BMD (FA-BMD) and IVDD when OP and FA-BMD were used as exposures. However, increased levels of total body BMD (TB-BMD) and lumbar spine BMD (LS-BMD) were revealed as significant risk factors for IVDD (TB-BMD: IVW, OR = 1.201, 95% CI: 1.123–1.284, P = 8.72 × 10−8; LS-BMD: IVW, OR = 1.179, 95% CI: 1.083–1.284, P = 1.43 × 10−4). Interestingly, both heel BMD (eBMD) and femur neck BMD (FN-BMD) exhibited potential causal relationships (eBMD: IVW, OR = 1.068, 95% CI: 1.008–1.131, P = 0.0248; FN-BMD, IVW, OR = 1.161, 95% CI: 1.041–1.295, P = 0.0074) with the risk of IVDD. The reverse MR analysis revealed no statistically causal impact of IVDD on OP and the level of BMD (P > 0.05).
Conclusion: OP and the level of FA-BMD were revealed to have no causal relationship with IVDD. The increased levels of TB-BMD and LS-BMD could promote the occurrence of IVDD. Both eBMD and FN-BMD have potential causal relationships with the risk of IVDD. No significant relationship exists between IVDD and the risk of OP. Further research is warranted to comprehensively comprehend the molecular mechanisms underlying the impact of OP and BMD on IVDD and vice versa.
1 Introduction
Intervertebral disc degeneration (IVDD) serves as the pathological basis for various spinal degenerative diseases that contribute to disability and reduce the quality of life of the affected patients (1). The intervertebral disc comprises the nucleus pulposus, the cartilage endplate, and the annulus fibrosus. Each of these three main components of the intervertebral disc is mostly composed of collagen and proteoglycan, both of which are crucial for imparting critical qualities to the disc (2). Intervertebral disc degeneration (IVDD) is a pathological condition characterized by a progressive decline in the levels of proteoglycans and the water content inside the nucleus pulposus. This degenerative condition is well-recognized as the primary contributor to the development of lower back pain (3). In severe cases, the destroyed discs compress the spinal nerve roots located at L4-S2, which often results in sciatica (4).
Osteoporosis, often referred to as OP, is a multifaceted skeletal disorder characterized by a loss in bone mass and impairment of bone microarchitecture, which results in heightened vulnerability to fractures and increased fragility of the bones (5). A key feature of osteoporosis is the reduction in bone mineral density (BMD). Dual-energy X-ray absorptiometry (DXA) is currently the primary method of establishing the clinical diagnosis of osteoporosis (6). OP is prevalent globally, affecting nearly 200 million individuals, particularly postmenopausal women (7). The important risk factors reported for OP include genetics, living habits, and medical history (5). The timely identification and management of osteoporosis (OP) is important as it would significantly contribute to preventing bone fractures and enhance the overall quality of life of the affected patients.
The correlation between osteoporosis and disc degeneration remains debatable in the scientific community to date. Liang et al. (8) demonstrated through a retrospective study that reduced vertebral bone mineral density (BMD), osteoporosis, and the exacerbation of disc degeneration were correlated. On the other hand, Kaiser et al. (9) demonstrated, through a prospective study, that higher trabecular BMD is negatively associated with disc height loss. Interestingly, certain retrospective and prospective studies have also reported no correlation between BMD and the progression of IVDD (10–13). These inconsistent conclusions could be originating from the confounders or covariates, which impact the association between BMD and IVDD.
IVDD and osteoporosis share a few genetic mechanisms and signaling pathways. Inflammatory cytokines, including interleukin-1 (IL-1), IL-6, and tumor necrosis factor-α (TNF-α), influence the development, progression, and severity of intervertebral disc degeneration (IVDD) (14–16) and OP (17–19). The matrix metallopeptidase 9 (MMP-9), which is an important proteolytic enzyme, is also closely associated with IVDD and OP. In addition, several same signaling pathways, including PI3K-AKT (20, 21), WNT/β-catenin (22, 23), and NF-κB (24, 25), are associated with intervertebral disc degeneration (IVDD) and also reported to significantly affect the pathogenesis of osteoporosis.
While numerous studies have reported the association between IVDD and OP, the direct causal relationship remains obscure. Mendelian randomization (MR) is a novel methodology that is being increasingly utilized to explore the causal associations between modifiable exposures and various illnesses or features (26, 27). MR allows for effectively overcoming certain limitations, especially unfeasible causality in the randomized controlled trials (RCT) and inevitable confounding bias or reverse causality in observational studies (26). MR also aids in distinguishing causal pathways from risk factors that are challenging to randomize or prone to measurement errors (28).
Our study indicates no discernible causal association between OP and the risk of intervertebral disc degeneration (IVDD) or between FA-BMD and the risk of IVDD. The elevated levels of TB-BMD and LS-BMD could, however, contribute to the development of intervertebral disc degeneration (IVDD). In addition, both eBMD and FN-BMD were revealed to have a putative causal relationship with the risk of intervertebral disc degeneration (IVDD). Notably, no substantial correlation was revealed between intervertebral disc degeneration (IVDD) and osteoporosis (OP) or among the bone mineral density across various anatomical locations and different age cohorts.
2 Methods
2.1 Data resources and sample information
The study design is depicted in Figure 1. The study was conducted by following the STROBE-MR guidelines (29). The data used in the present study were the summary genome-wide association studies (GWAS) data of intervertebral disc degeneration (IVDD) and osteoporosis, which were obtained from the FinnGen collaboration. The IVDD dataset comprised 20,001 cases and 164,692 controls, whereas the osteoporosis dataset comprised 3,203 cases and 209,575 controls (30). The genome-wide association study (GWAS) data on lumbar spine bone mineral density (LS-BMD), femur neck bone mineral density (FN-BMD), and forearm bone mineral density (FA-BMD) collected by the Genetic Factors for Osteoporosis (GEFOS) consortium were also used. The dataset comprised 25,509 cases for FN-BMD, 12,906 cases for LS-BMD, and 563 cases for FA-BMD (31, 32). In addition, heel bone mineral density (eBMD) data were obtained from the MRC-IEU, and the dataset comprised 265,627 individuals. The data for TB-BMD for different age groups were obtained from the reports of a large GWAS meta-analysis study (32). The TB-BMD, LS-BMD, FN-BMD, and FA-BMD assessments were conducted using dual-energy X-ray absorptiometry (DXA) (6). The eBMD assessments were conducted using ultrasonography (33). A reverse analysis was conducted to obtain the genetic instruments for OP and bone mineral density (BMD) at different anatomical sites, and their impact on intervertebral disc degeneration (IVDD) was determined. IVDD was diagnosed using the International Classification of Diseases 10th Revision (ICD-10) M51, ICD-9 722, and ICD-8 275, with removed ICD-9 7220|7224|7227|7228Av and ICD-8 7250. Detailed information of these GWAS summary data is provided in Supplementary Table 1.
Figure 1 The flowchart for the two-sample bidirectional Mendelian randomization analysis. The blue line represents the Mendelian randomization analysis of the causal relationship of intervertebral disc degeneration (IVDD) with bone mineral density (BMD) or osteoporosis. The red line represents the Mendelian randomization analysis of the causal relationship of bone mineral density or osteoporosis with IVDD. The following three assumptions had to be considered: the selected instrument variants (SNPs) should be significantly associated with the exposure; the variants should be independent of any exposure–outcome relationship confounders; the variants affect the outcome only via the exposure.
2.2 MR analysis
2.2.1 Selection of genetic variants
In order to address the issue of weak instrument bias, a genome-wide significance threshold of P < 5 × 10−8 was used as the default criterion for identifying the single-nucleotide polymorphisms (SNPs), whereas the threshold of minor allele frequency (MAF) > 0.01 was used as the genetic instruments for osteoporosis (OP), bone mineral density (BMD), and intervertebral disc degeneration (IVDD). Next, the clumping threshold (the genetic distance = 10,000 kb and the SNP linkage disequilibrium value r2 < 0.001) was established using the data from the 1000 genomes project for European ancestry (34) to resolve the linkage disequilibrium (LD) issue based on the identified SNPs.
Diabetes, BMI, and smoking are significant prognostic factors for both OP and IVDD. Several scholars have proposed that the risk factors for diabetes mellitus (35, 36) and body mass index (BMI) are associated with IVDD (37). Moreover, a retrospective study indicated that cigarette smoking accelerates the development of cervical disc degeneration (38). Interestingly, a higher level of education was reported as a risk factor for cervical disc degeneration, regardless of age differences among the respondents (39). OP is also reported to be closely associated with diabetes mellitus, BMI (40), smoking (41, 42), and the level of education (43). BMI and smoking are considered to have significant effects on OP, particularly in postmenopausal women. Accordingly, to ensure obtaining statistically meaningful outcomes, the present study used PhenoScannerV2 to identify and exclude certain confounding single-nucleotide polymorphisms (SNPs) associated with diabetes, body mass index (BMI), smoking, and social status in our study.
The strength of the IVs was assessed by computing the F-statistics using the formula F = R2 × (N – 2)/(1 – R2). In the formula, R2 denotes the proportion of variation in the exposure variable explained by the IV, and N denotes the sample size of the original GWAS used as the outcome variable (44). The R2 for each IV was computed using the formula R2 = [2 × EAF × (1 – EAF) × beta2)/[(2 × EAF × (1 – EAF) × beta2] + [2 × EAF × (1 – EAF) × N × (SE × beta2)]. In this formula, EAF denotes the effect of allele frequency, beta represents the estimated genetic effect on the outcome, N denotes the sample size of the genome-wide association study (GWAS), and SE denotes the standard error of the genetic effect (45). The IVs with the F statistics over 10 were considered credible instrumental variables and selected for the subsequent Mendelian randomization study.
2.2.2 Sensitivity analysis and statistical analysis
The potential causal relationship between intervertebral disc abnormalities and osteoporosis was investigated in the present study primarily using the inverse variance weighted (IVW) method, as described in a previous report (46). The Cochrane’s Q test was conducted to evaluate the heterogeneity among the single-nucleotide polymorphisms (SNPs). Furthermore, the strength and reliability of the primary results obtained in the previous steps were assessed using various sensitivity analyses, including the weighted-median approach and the IVW method. Moreover, the possible presence of directional pleiotropy was evaluated using MR-Egger regression (47). The slope of MR-Egger regression indicates the causal estimates that have been adjusted for pleiotropy, whereas the intercept’s value provides an estimation of the extent of pleiotropy. In addition, the Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) test was conducted as a supplementary approach to identify and address the horizontal pleiotropic outliers (48). A global test of heterogeneity was conducted through a regression analysis of the relationships between the single-nucleotide polymorphisms (SNPs) and the outcomes while considering the connections between the SNPs and the exposures. The observed distance of each SNP from the regression was then compared with the anticipated distance according to the null hypothesis of no pleiotropy. The robustness of the obtained results was assessed using “leave-one-out” analyses, which involved systematically eliminating one single SNP at a time. This was followed by a reanalysis of the Mendelian randomization (MR) results using the IVW approach with the remaining instrumental variables (IVs).
2.3 Calculation of statistical power
The statistical power was evaluated using the mRnd website (https://shiny.cnsgenomics.com/mRnd/) (49). The key determinants of statistical power were the size of the sample for the result and the extent to which the genetic instrument accounted for the variance in the exposure variable.
All analyses were performed using R (version 4.1.2) and the related R packages (TwoSampleMR and forestplot). A P-value of <0.005 (0.05/10) indicated strong evidence of a causal association after the Bonferroni correction threshold was applied. A P-value ranging between 0.05 and 0.005 was considered suggestive evidence for a potential causal association.
3 Results
3.1 Causal effects of osteoporosis and bone mineral density on IVDD
All the IVs that were ultimately selected exhibited F-statistic values of over 10. Details of the method adopted to screen out the single-nucleotide polymorphisms (SNPs) are provided in Supplementary Table 2. In the analysis of the effect of eBMD on the risk of IVDD, rs7816131 was excluded because of its robust association with IVDD (P = 8.3 × 10−9), which is not in accordance with the third assumption of MR analysis. Finally, 2 SNPs of OP, 80 SNPs of TB-BMD, 22 SNPs of LS-BMD, 20 SNPs of FN-BMD, 3 SNPs of FA-BMD, and 331 SNPs of eBMD were used as IVs to analyze the relationship between OP/BMDs and the risk of intervertebral disc degeneration (IVDD). Comprehensive data on the IVs for BMD are provided in Supplementary Table 3. The variation explained by these IVs was 0.02% for OP, 9.05% for TB-BMD, 2.1% for FN-BMD, 1.65% for FA-BMD, 2.27% for LS-BMD, and 14.27% for eBMD. The statistical power values determined in the Mendelian randomization analyses of OP and BMD in terms of their effects on IVDD are presented in Supplementary Table 5.
The MR analysis did not yield statistically significant evidence of the causal effects of OP (IVW, P > 0.05) and FA-BMD (IVW, P > 0.05) on IVDD (Figure 2). However, positive correlations of the BMDs (TB-BMD, LS-BMD, and eBMD) with the risk of IVDD were revealed. The results indicated a significant association between the genetically enhanced BMD levels and a higher susceptibility to IVDD (TB-BMD: IVW, OR = 1.201, 95% CI: 1.123–1.284, P = 8.72 × 10−8; LS-BMD: IVW, OR = 1.179, 95% CI: 1.083–1.284, P = 1.43 × 10−4; Figure 2). The different results may be attributed to too few amounts of SNPs associated with OP and FA-BMD which served as exposures. Thus, the results of the MR-Egger and (or) MR-PRESSO test are not available. Moreover, heterogeneity (P for Cochrane’s Q in IVW = 0.049 or 0.005 when OP and FA-BMD were served as exposures, Table 1) may exert an influence on this difference (Table 1).
Figure 2 Causal effect of OP and BMD at different anatomical sites on IVDD. IVDD, intervertebral disc degeneration; OP, osteoporosis; TB-BMD, total body bone mineral density; FN-BMD, femoral neck bone mineral density; FA-BMD, forearm bone mineral density; LS-BMD, lumbar spine bone mineral density; eBMD, heel bone mineral density; IVW, inverse variance weighted; nsnp, number of single nucleotide polymorphisms; CI, confidence interval.
Table 1 MR sensitivity analyses of IVDD and BMD at different sites and in different age groups.
MR sensitivity analysis of TB-BMD and IVDD exerted heterogeneous (P for Cochrane’s Q in IVW = 0.007, Table 1) and no horizontal pleiotropy (P for MR-Egger intercept = 0.171, Table 1) when TB-BMD was served as the exposure. The MR-Egger analysis exhibited less statistical power, as seen by the lack of a significant P value and broader confidence intervals when compared with the IVW technique (50). Therefore, the MR-PRESSO Outlier Test, another horizontal pleiotropy test, was performed. The SNP rs4846580 was identified as an outlier using the MR-PRESSO Outlier Test, with an observed residual sum of squares (RSSobs) value of 1.92 × 10−3 (P < 0.08). Therefore, this SNP was removed during the analysis of the relationship between TB-BMD and the risk of IVDD. The results of the MR analysis indicated no statistically significant association between TB-BMD and IVDD (IVW, OR = 1.073, 95% CI: 1.014–1.135, P = 1.47 × 10−2) when the analysis was repeated with the remaining five single-nucleotide polymorphisms (SNPs) (Supplementary Figure 5). Notably, neither heterogeneity nor horizontal pleiotropy was indicated (P for Cochrane’s Q in IVW = 0.349; P for MR-Egger intercept = 0.016, P for MR-PRESSO Global Test = 0.37, Table 1) in the MR sensitivity analysis of LS-BMD with IVDD.
Interestingly, both eBMD and FN-BMD were revealed to have potential causal relationships (eBMD: IVW, OR = 1.068, 95% CI: 1.008–1.131, P = 0.0248; FN-BMD, IVW, OR = 1.161, 95% CI: 1.041–1.295, P = 7.40 × 10−3; Figure 2) with the risk of IVDD according to statistical significance standard (P value threshold) defined in the Method section. Neither heterogeneity nor horizontal pleiotropy was indicated in the directional pleiotropy in the MR sensitivity analysis of FN-BMD with IVDD (P for Cochran’s Q in IVW = 0.118, MR-Egger intercept = −0.029, P = 0.097; P for MR-PRESSO Global Test = 0.37, Table 1). However, horizontal pleiotropy and heterogeneity (P for MR-PRESSO Global Test <0.001, P for Cochrane’s Q in IVW <0.0001, Table 1) were indicated in the MR sensitivity analysis of eBMD with IVDD, although P for the MR-Egger intercept is more than 0.05 (Table 1). This may be attributed to too many amounts of SNPs associated with eBMD that served as the exposures. The scatter plots for the effect sizes of the SNPs for OP and TB-BMD at different anatomical sites on IVDD are presented in Supplementary Figure 1. The forest plots, leave-one-out analysis plots, and funnel plots for the causal effect of BMD on IVDD are depicted in Supplementary Figures 2-4. When the leave-one-out analysis was conducted using the IVW method, most of the determined correlations remained unchanged even when considering a single SNP associated with bone mineral density (BMD).
3.2 Causal effect of IVDD on BMD at different anatomical sites or in different age groups
Six SNPs (rs3010043, rs4473430, rs3135840, rs6470763, rs4284332, and rs17487277) were selected as IVs for determining the causal effect of IVDD on the risk of OP, TB-BMD, FN-BMD, FA-BMD, and LS-BMD. Notably, rs6470763 was directly associated with heel BMD (beta = 0.021, P = 7.60 × 10−11), which is not in accordance with the third assumption of MR analysis. Therefore, this SNP was removed from the analysis of the relationship between intervertebral disc degeneration (IVDD) and the risk of eBMD. All IVs that were finally selected for the analysis had F-statistic values over 10. Detailed information on the IVs strongly associated with IVDD is provided in Supplementary Table 4. These IVs accounted for a variance of 0.11% or 1.20% in the analysis of the effect of IVDD on BMD.
The MR analysis revealed that IVDD was not associated with BMD at different anatomical sites (Figure 3) and in different age groups (Figure 4). Neither heterogeneity nor horizontal pleiotropy (all P for Cochran’s Q in IVW, MR-Egger intercept, and MR-PRESSO Global Test are more than 0.05) was indicated in the MR sensitivity analysis of IVDD with BMD at different sites, and in different age groups, excluded heterogeneity was found in the analysis of IVDD with TB-BMD and TB-BMD (age 0–15) (P for Cochrane’s Q in IVW is less than 0.05). Otherwise, the scatter plots for the effect sizes of SNPs for IVDD’s relationship to OP and TB-BMD at different anatomical sites are depicted in Supplementary Figures 6, 7. The forest plots, leave-one-out analysis plots, and funnel plots for the causal effect of IVDD on BMDs are depicted in Supplementary Figures 8, 9.
Figure 3 Causal effect of IVDD on OP and BMD at different anatomical sites. IVDD, intervertebral disc degeneration; OP, osteoporosis; TB-BMD, total body bone mineral density; FN-BMD, femoral neck bone mineral density; FA-BMD, forearm bone mineral density; LS-BMD, lumbar spine bone mineral density; eBMD, heel bone mineral density; IVW, inverse variance weighted; nsnp, number of single-nucleotide polymorphisms; CI, confidence interval.
Figure 4 Causal effect of IVDD on BMD in different age groups. TB-BMD, total body bone mineral density; IVW, inverse variance weighted; nsnp, number of single-nucleotide polymorphisms; OR, odds ratio; CI, confidence interval.
4 Discussion
The present study revealed no significant evidence to support the causal effect of OP and FA-BMD on the risk of IVDD. Moreover, the positive relationships of TB-BMD and LS-BMD with the risk of IVDD were revealed. In addition, both eBMD and FN-BMD exhibited potential causal relationships with the risk of IVDD. In reverse MR analysis, no causal effects of IVDD on the risk of OP and the change in bone mineral density were revealed.
The scientific community continues to debate whether bone mineral density (BMD) affects IVDD and, if so, how. While certain studies have reported osteoporosis as a causal factor of IVDD, others have considered it a protective factor against IVDD. Otherwise, several studies demonstrate there is no association between OP and IVDD (12, 51).
It has long been suggested that diminished bone quality is associated with the gradual deterioration of the endplates and the development of spondylosis, which ultimately results in heightened disc degeneration (52). Nevertheless, the incidence of intervertebral disc degeneration (IVDD) is lower among individuals with a low bone mineral density (BMD), even though these people are more susceptible to vertebral body fractures (53–55). According to this concept, it is postulated that osteoporosis, a degenerative and incapacitating disorder associated with aging and affecting a significant population globally (56), potentially delays the onset of intervertebral disc degeneration (IVDD). Functional investigations aimed at elucidating the correlation between the genetic factors influencing bone mineral density (BMD) and intervertebral disc degeneration (IVDD) are expected to provide valuable insights into this association and facilitate the discovery of novel treatment approaches.
Mechanistically, some scholars believe that osteoporosis leads to lower vertebral BMD, aggravates disc load, reduces the supply of nutrients to the intervertebral disc (57), and increases inflammatory factors that cause disc degeneration (15, 58). Others consider osteoporosis to delay disc degeneration that the vertebral body with lower BMD has a loose bony microstructure, which allows the vascular buds that nourish the cartilage endplate to grow better, and the increase in the number of vascular buds to enrich the blood supply of the cartilage endplate and can better provide nutrients to the intervertebral discs, thereby delaying IVDD (59, 60). In terms of genetic predictions, our results support the latter opinion that increased levels of TB-BMD and LS-BMD are the risk factors for IVDD, whereas eBMD and FN-BMD were revealed as the potential risk factors for IVDD. On the other hand, the results of the present study revealed no causal effect of OP on IVDD when just two SNPs were selected as IVs. These different results could be attributed to the different number of SNPs associated with the exposure and types of exposure variables. Furthermore, no statistically significant correlation was observed between intervertebral disc degeneration (IVDD) and osteoporosis (OP) or among the bone mineral density at different anatomical locations and in different age groups.
The present MR study offers several advantages. The use of Mendelian randomization effectively mitigated the potential influences of confounding bias and reverse causation. In addition, several sensitivity analyses were conducted to assess the robustness of the three Mendelian randomization assumptions, which reduced the likelihood of spurious findings. The instrumental variable weighted (IVW) method, which was adopted as the principal approach in the present study, resulted in a superior statistical power compared with other Mendelian randomization (MR) approaches, particularly the MR-Egger approach (50). Consequently, the MR-Egger analysis led to lower statistical power, as evident in the lack of a significant P value and broader confidence intervals compared with the IVW technique. The findings of the present study provided additional support for the necessity of maintaining a constant beta direction across all magnetic resonance (MR) procedures. Furthermore, harmonized data were consistently identified and rectified using the MR-PRESSO method (48). This ensured the absence of horizontal pleiotropy throughout the Mendelian randomization (MR) analysis and enhanced the reliability of the obtained findings.
Nonetheless, similar to other studies, the present study also had certain limitations. The first limitation is that the study participants were exclusively of European origin, which may cause the findings of this study not to be generalizable to individuals of other nations, such as those of African or East Asian descent. Furthermore, the comprehensive elimination of pleiotropy was challenging due to the limited understanding of the overall biological functionality of these instrumental variations. In addition, while the findings of the present study indicated the existence of potential causal relationships between bone mineral density at various anatomical sites and the risk of intervertebral disc degeneration (IVDD), a further comprehensive investigation of the intricate underlying processes is nonetheless warranted. Lastly, it is difficult to validate the findings of the present study in wet lab assays.
In summary, the present study revealed no substantial causal relationship, either directly from the causative impact of OP or indirectly via FA-BMD, on IVDD. The elevation in the levels of TB-BMD (total body bone mineral density) and LS-BMD (lumbar spine bone mineral density) could, however, contribute to the development of intervertebral disc degeneration (IVDD). Both eBMD and FN-BMD exhibited a putative causal relationship with the risk of intervertebral disc degeneration (IVDD). Furthermore, the association of intervertebral disc degeneration (IVDD) with osteoporosis (OP) and bone mineral density (BMD) was not statistically significant. Therefore, further investigation is required to elucidate the molecular mechanisms underlying the impact of bone mineral density (BMD) on intervertebral disc degeneration (IVDD), which will facilitate the precise treatment of these two bone degenerative diseases by orthopedic surgeons.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
GL: Investigation, Validation, Writing – review & editing, Conceptualization, Data curation, Formal analysis, Writing – original draft, Visualization, Software, Project administration, Methodology. HZ: Writing – review & editing, Supervision, Resources, Funding acquisition, Conceptualization. MC: Writing – original draft, Visualization, Validation, Methodology, Formal analysis, Data curation, Conceptualization. WC: Funding acquisition, Supervision, Writing – review & editing, Investigation, Validation.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Natural Science Foundation in Hunan Province (2022JJ40391) and the Health Department of Hunan Province (B202304078506).
Acknowledgments
We thank participants and staff of the FinnGen consortium, the UK Biobank, the MRC Integrative Epidemiology Unit, and the European Bioinformatics Institute.
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.2024.1298531/full#supplementary-material
Abbreviations
GEFOS, the Genetic Factors for Osteoporosis; MRC-IEU, The Medical Research Council Integrative Epidemiology Unit; IVDD, intervertebral disc degeneration; OP, osteoporosis; TB-BMD, total body bone mineral density; FN-BMD, femoral neck bone mineral density; FA-BMD, forearm bone mineral density; LS-BMD, lumbar spine bone mineral density; eBMD, heel bone mineral density; NA, not available; IVW, inverse variance weighted; nsnp, number of single-nucleotide polymorphism; OR, odd ratios; CI, confidence interval.
References
1. Li Z, Chen X, Xu D, Li S, Chan MTV, Wu WKK. Circular RNAs in nucleus pulposus cell function and intervertebral disc degeneration. Cell proliferat. (2019) 52:e12704. doi: 10.1111/cpr.12704
2. Zehra U, Tryfonidou M, Iatridis JC, Illien-Jünger S, Mwale F, Samartzis D. Mechanisms and clinical implications of intervertebral disc calcification. Nat Rev Rheumatol. (2022) 18:352–62. doi: 10.1038/s41584-022-00783-7
3. Knezevic NN, Candido KD, Vlaeyen JWS, Van Zundert J, Cohen SP. Low back pain. Lancet (London England). (2021) 398:78–92. doi: 10.1016/S0140-6736(21)00733-9
5. Lane NE. Epidemiology, etiology, and diagnosis of osteoporosis. Am J obstetrics gynecol. (2006) 194:S3–11. doi: 10.1016/j.ajog.2005.08.047
6. Anam AK, Insogna K. Update on osteoporosis screening and management. Med Clinics North America. (2021) 105:1117–34. doi: 10.1016/j.mcna.2021.05.016
7. Reginster JY, Burlet N. Osteoporosis: a still increasing prevalence. Bone. (2006) 38:S4–9. doi: 10.1016/j.bone.2005.11.024
8. Liang X, Liu Q, Xu J, Ding W, Wang H. Hounsfield unit for assessing bone mineral density distribution within cervical vertebrae and its correlation with the intervertebral disc degeneration. Front Endocrinol. (2022) 13:920167. doi: 10.3389/fendo.2022.920167
9. Kaiser J, Allaire B, Fein PM, Lu D, Jarraya M, Guermazi A, et al. Correspondence between bone mineral density and intervertebral disc degeneration across age and sex. Arch osteoporosis. (2018) 13:123. doi: 10.1007/s11657-018-0538-1
10. Yang Z, Griffith JF, Leung PC, Lee R. Effect of osteoporosis on morphology and mobility of the lumbar spine. Spine. (2009) 34:E115–121. doi: 10.1097/BRS.0b013e3181895aca
11. Salo S, Leinonen V, Rikkonen T, Vainio P, Marttila J, Honkanen R, et al. Association between bone mineral density and lumbar disc degeneration. Maturitas. (2014) 79:449–55. doi: 10.1016/j.maturitas.2014.09.003
12. Pan J, Lu X, Yang G, Han Y, Tong X, Wang Y. Lumbar disc degeneration was not related to spine and hip bone mineral densities in Chinese: facet joint osteoarthritis may confound the association. Arch osteoporosis. (2017) 12:20. doi: 10.1007/s11657-017-0315-6
13. Geng J, Wang L, Li Q, Huang P, Liu Y, Blake GM, et al. The association of lumbar disc herniation with lumbar volumetric bone mineral density in a cross-sectional chinese study. Diagn (Basel Switzerland). (2021) 11:938. doi: 10.3390/diagnostics11060938
14. Séguin CA, Pilliar RM, Roughley PJ, Kandel RA. Tumor necrosis factor-alpha modulates matrix production and catabolism in nucleus pulposus tissue. Spine. (2005) 30:1940–8. doi: 10.1097/01.brs.0000176188.40263.f9
15. Le Maitre CL, Hoyland JA, Freemont AJ. Catabolic cytokine expression in degenerate and herniated human intervertebral discs: IL-1beta and TNFalpha expression profile. Arthritis Res Ther. (2007) 9:R77. doi: 10.1186/ar2275
16. Shamji MF, Setton LA, Jarvis W, So S, Chen J, Jing L, et al. Proinflammatory cytokine expression profile in degenerated and herniated human intervertebral disc tissues. Arthritis rheumatism. (2010) 62:1974–82. doi: 10.1002/art.27444
17. Zheng SX, Vrindts Y, Lopez M, De Groote D, Zangerle PF, Collette J, et al. Increase in cytokine production (IL-1 beta, IL-6, TNF-alpha but not IFN-gamma, GM-CSF or LIF) by stimulated whole blood cells in postmenopausal osteoporosis. Maturitas. (1997) 26:63–71. doi: 10.1016/S0378-5122(96)01080-8
18. Ni Y, Li H, Zhang Y, Zhang H, Pan Y, Ma J, et al. Association of IL-6 G-174C polymorphism with bone mineral density. J Bone mineral Metab. (2014) 32:167–73. doi: 10.1007/s00774-013-0477-2
19. Liao L, Su X, Yang X, Hu C, Li B, Lv Y, et al. TNF-α Inhibits FoxO1 by Upregulating miR-705 to Aggravate Oxidative Damage in Bone Marrow-Derived Mesenchymal Stem Cells during Osteoporosis. Stem Cells (Dayton Ohio). (2016) 34:1054–67. doi: 10.1002/stem.2274
20. Yang S, Zhang F, Ma J, Ding W. Intervertebral disc ageing and degeneration: The antiapoptotic effect of oestrogen. Ageing Res Rev. (2020) 57:100978. doi: 10.1016/j.arr.2019.100978
21. Ma Y, Ran D, Zhao H, Song R, Zou H, Gu J, et al. Cadmium exposure triggers osteoporosis in duck via P2X7/PI3K/AKT-mediated osteoblast and osteoclast differentiation. Sci Tot Environ. (2021) 750:141638. doi: 10.1016/j.scitotenv.2020.141638
22. Cheng CH, Chen LR, Chen KH. Osteoporosis due to hormone imbalance: an overview of the effects of estrogen deficiency and glucocorticoid overuse on bone turnover. Int J Mol Sci. (2022) 23:1376. doi: 10.3390/ijms23031376
23. Hao Y, Ren Z, Yu L, Zhu G, Zhang P, Zhu J, et al. p300 arrests intervertebral disc degeneration by regulating the FOXO3/Sirt1/Wnt/β-catenin axis. Aging Cell. (2022) 21:e13677. doi: 10.1111/acel.13677
24. Rachner TD, Khosla S, Hofbauer LC. Osteoporosis: now and the future. Lancet (London England). (2011) 377:1276–87. doi: 10.1016/S0140-6736(10)62349-5
25. Tao H, Li W, Zhang W, Yang C, Zhang C, Liang X, et al. Urolithin A suppresses RANKL-induced osteoclastogenesis and postmenopausal osteoporosis by, suppresses inflammation and downstream NF-κB activated pyroptosis pathways. Pharmacol Res. (2021) 174:105967. doi: 10.1016/j.phrs.2021.105967
26. 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: JASN. (2016) 27:3253–65. doi: 10.1681/ASN.2016010098
27. Swanson SA, Tiemeier H, Ikram MA, Hernán MA. Nature as a trialist?: deconstructing the analogy between mendelian randomization and randomized trials. Epidemiol (Cambridge Mass.). (2017) 28:653–9. doi: 10.1097/EDE.0000000000000699
28. Bowden J, Holmes MV. Meta-analysis and Mendelian randomization: A review. Res synthesis Methods. (2019) 10:486–96. doi: 10.1002/jrsm.1346
29. Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomization: the STROBE-MR statement. Jama. (2021) 326:1614–21. doi: 10.1001/jama.2021.18236
30. Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. (2023) 613:508–18. doi: 1038/s41586-022-05473-8
31. 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. doi: 10.1038/nature14878
32. Medina-Gomez C, Kemp JP, Trajanoska K, Luan J, Chesi A, Ahluwalia TS, et al. Life-course genome-wide association study meta-analysis of total body BMD and assessment of age-specific effects. Am J Hum Genet. (2018) 102:88–102. doi: 10.1016/j.ajhg.2017.12.0
33. Howard GM, Nguyen TV, Harris M, Kelly PJ, Eisman JA. Genetic and environmental contributions to the association between quantitative ultrasound and bone mineral density measurements: a twin study. J Bone mineral rRes: Off J Am Soc Bone Mineral Res. (1998) 13:1318–27. doi: 10.1359/jbmr.1998.13.8.1318
34. Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, et al. A map of human genome variation from population-scale sequencing. Nature. (2010) 467:1061–73. doi: 10.1038/nature09534
35. Fabiane SM, Ward KJ, Iatridis JC, Williams FM. Does type 2 diabetes mellitus promote intervertebral disc degeneration? Eur Spine J. (2016) 25:2716–20. doi: 10.1007/s00586-016-4612-3
36. Teraguchi M, Yoshimura N, Hashizume H, Yamada H, Oka H, Minamide A, et al. Progression, incidence, and risk factors for intervertebral disc degeneration in a longitudinal population-based cohort: the Wakayama Spine Study. Osteoarthritis Cartilage. (2017) 25:1122–31. doi: 10.1016/j.joca.2017.01.001
37. Hangai M, Kaneoka K, Kuno S, Hinotsu S, Sakane M, Mamizuka N, et al. Factors associated with lumbar intervertebral disc degeneration in the elderly. Spine J. (2008) 8:732–40. doi: 10.1016/j.spinee.2007.07.392
38. Chen Z, Li X, Pan F, Wu D, Li H. A retrospective study: Does cigarette smoking induce cervical disc degeneration? Int J Surg. (2018) 53:269–73. doi: 10.1016/j.ijsu.2018.04.004
39. Markotić V, Zubac D, Miljko M, Šimić G, Zalihić A, Bogdan G, et al. Level of education as a risk factor for extensive prevalence of cervical intervertebral disc degenerative changes and chronic neck pain. Cent Eur J Public Health. (2017) 25:245–50. doi: 10.21101/cejph.a4897
40. Tang G, Feng L, Pei Y, Gu Z, Chen T, Feng Z. Low BMI, blood calcium and vitamin D, kyphosis time, and outdoor activity time are independent risk factors for osteoporosis in postmenopausal women. Front Endocrinol (Lausanne). (2023) 14:1154927. doi: 10.3389/fendo.2023.1154927
41. Brook JS, Balka EB, Zhang C. The smoking patterns of women in their forties: their relationship to later osteoporosis. Psychol Rep. (2012) 110:351–62. doi: 10.2466/13.18.PR0.110.2.351-362
42. Trevisan C, Alessi A, Girotti G, Zanforlini BM, Bertocco A, Mazzochin M, et al. The impact of smoking on bone metabolism, bone mineral density and vertebral fractures in postmenopausal women. J Clin Densitom. (2020) 23:381–9. doi: 10.1016/j.jocd.2019.07.007
43. Tabor E, Grodzki A, Pluskiewicz W. Higher education and better knowledge of osteoporosis improve bone health in Polish postmenopausal women. Endokrynol Pol. (2022) 73:831–6. doi: 10.5603/EP.a2022.0055
44. Burgess S, Thompson SG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. (2011) 40:755–64. doi: 10.1093/ije/dyr036
45. Papadimitriou N, Dimou N, Tsilidis KK, Banbury B, Martin RM, Lewis SJ, et al. Physical activity and risks of breast and colorectal cancer: a Mendelian randomisation analysis. Nat Commun. (2020) 11:597. doi: 10.1038/s41467-020-14389-8
46. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758
47. 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
48. Verbanck M, Chen C-Y, 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
49. Brion M-JA, Shakhbazov K, Visscher PM. Calculating statistical power in Mendelian randomization studies. Int J Epidemiol. (2012) 42:1497–501. doi: 10.1093/ije/dyt179
50. Lin Z, Deng Y, Pan W. Combining the strengths of inverse-variance weighting and Egger regression in Mendelian randomization using a mixture of regressions model. PloS Genet. (2021) 17:e1009922. doi: 10.1371/journal.pgen.1009922
51. Wang YX, Kwok AW, Griffith JF, Leung JC, Ma HT, Ahuja AT, et al. Relationship between hip bone mineral density and lumbar disc degeneration: a study in elderly subjects using an eight-level MRI-based disc degeneration grading system. J magnetic resonance Imag: JMRI. (2011) 33:916–20. doi: 10.1002/jmri.22518
52. Yong-Hing K, Kirkaldy-Willis WH. The pathophysiology of degenerative disease of the lumbar spine. Orthopedic Clinics North America. (1983) 14:491–504. doi: 10.1016/S0030-5898(20)31329-8
53. Miyakoshi N, Itoi E, Murai H, Wakabayashi I, Ito H, Minato T. Inverse relation between osteoporosis and spondylosis in postmenopausal women as evaluated by bone mineral density and semiquantitative scoring of spinal degeneration. Spine. (2003) 28:492–5. doi: 10.1097/01.BRS.0000048650.39042.58
54. Nanjo Y, Morio Y, Nagashima H, Hagino H, Teshima R. Correlation between bone mineral density and intervertebral disk degeneration in pre- and postmenopausal women. J Bone mineral Metab. (2003) 21:22–7. doi: 10.1007/s007740300004
55. Mattei TA. Osteoporosis delays intervertebral disc degeneration by increasing intradiscal diffusive transport of nutrients through both mechanical and vascular pathophysiological pathways. Med Hypotheses. (2013) 80:582–6. doi: 10.1016/j.mehy.2013.01.030
56. Cooper C, Campion G, Melton LJ 3rd. Hip fractures in the elderly: a world-wide projection. Osteoporosis international: J established as result cooperation between Eur Foundation Osteoporosis Natl Osteoporosis Foundation USA. (1992) 2:285–9. doi: 10.1007/BF01623184
57. Ou-Yang L, Lu GM. Dysfunctional microcirculation of the lumbar vertebral marrow prior to the bone loss and intervertebral discal degeneration. Spine. (2015) 40:E593–600. doi: 10.1097/BRS.0000000000000834
58. Risbud MV, Shapiro IM. Role of cytokines in intervertebral disc degeneration: pain and disc content. Nat Rev Rheumatol. (2014) 10:44–56. doi: 10.1038/nrrheum.2013.160
59. Turgut M, Uslu S, Uysal A, Yurtseven ME, Ustün H. Changes in vascularity of cartilage endplate of degenerated intervertebral discs in response to melatonin administration in rats. Neurosurgical Rev. (2003) 26:133–8. doi: 10.1007/s10143-003-0259-8
Keywords: osteoporosis, bone mineral density, intervertebral disc degeneration, Mendelian randomization, genome-wide association studies
Citation: Liu G, Zhang H, Chen M and Chen W (2024) Causal relationship between intervertebral disc degeneration and osteoporosis: a bidirectional two-sample Mendelian randomization study. Front. Endocrinol. 15:1298531. doi: 10.3389/fendo.2024.1298531
Received: 21 September 2023; Accepted: 02 April 2024;
Published: 30 April 2024.
Edited by:
Xiaonan Liu, Johns Hopkins Medicine, United StatesReviewed by:
Dipti Deshpande, Kirkland and Ellis, United StatesJi Tu, University of New South Wales, Australia
Vladimir Palicka, University Hospital Hradec Kralove, Czechia
Copyright © 2024 Liu, Zhang, Chen and Chen. 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: Wenkang Chen, 2018011264@usc.edu.cn
†These authors have contributed equally to this work