Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 23 August 2022
Sec. Obesity

Causal relationships of obesity on musculoskeletal chronic pain: A two-sample Mendelian randomization study

Xiaoqing ChenXiaoqing ChenHaifeng TangHaifeng TangJinding LinJinding LinRongdong Zeng*Rongdong Zeng*
  • Department of Orthopaedics, Quanzhou First Hospital Affiliated to Fujian Medical University, Quanzhou, China

Background: The association between obesity and musculoskeletal chronic pain has attracted much attention these days; however, the causal relationship between them is uncertain. Hence, this study performed a Mendelian randomization (MR) analysis to investigate the causal effects of body mass index (BMI), waist circumference (WC), hip circumference (HC), and waist-to-hip ratio (WHR) on knee pain, hip pain, and back pain.

Materials and methods: The summary data for obesity and musculoskeletal chronic pain came from the genome-wide association study datasets. Significant and independent (p < 5 × 10−8; r2 < 0.001, kb = 10,000) single-nucleotide polymorphisms were extracted for MR analysis. The inverse variance weighted (IVW) and other methods were used for MR analysis, while sensitivity analyses were conducted to test the reliability and stability.

Results: The positive causal effects of BMI on knee pain (odds ratio (OR) = 1.049; 95% CI: 1.034 to 1.063; p = 9.88 × 10−12), hip pain (OR = 1.034; 95% CI: 1.024 to 1.044; p = 1.38 × 10−12), and back pain (OR = 1.022; 95% CI: 1.007 to 1.038; p = 0.004) were observed. WC and HC were also positively associated with knee pain (WC: OR = 1.057; 95% CI: 1.041 to 1.072; p = 1.54 × 10−13; HC: OR = 1.034; 95% CI: 1.017 to 1.052; p = 1.32 × 10−4) and hip pain (WC: OR = 1.031; 95% CI: 1.020 to 1.042; p = 2.61 × 10−8; HC: OR = 1.027; 95% CI: 1.018 to 1.035; p = 5.48 × 10−10) but not back pain. No causal relationship was found between WHR and musculoskeletal chronic pain. The results were robust according to sensitivity tests.

Conclusions: This study revealed that BMI was positively related to knee, hip, and back pain and that WC and HC were positively associated with knee and hip pain, while WHR was not related to any type of musculoskeletal chronic pain.

Introduction

Obesity is defined as a disproportionate body weight for height with an excessive accumulation of adipose tissue, which has affected over 640 million population around the world and caused a great burden on society (1). According to the classification of the World Health Organization, individuals are classified as obese when their body mass index (BMI) is over 30 kg/m2 (2). Previous studies had illustrated that individuals with obesity were at greater risk for many health problems including metabolic syndrome, type 2 diabetes, osteoarthritis, and other diseases as compared with the normal-weight population (35). Chronic pain, as defined by the International Association for the Study of Pain, refers to an unpleasant sensory and emotional experience associated with actual or potential tissue damage, containing back pain, musculoskeletal disorders, and neck pain (6, 7). Chronic pain costs US$560 to US$635 billion per year with a prevalence rate of 11% to 40% (8, 9). Therefore, it is of great significance to explore the relationships between these two heavy-burden and costly diseases.

The association between obesity and pain could be controversial, though many studies about their relationship have been performed. A prospective study including 285 patients and 191 volunteers revealed that shoulder pain was associated with obesity, and it would be helpful to treat pain by losing weight moderately (10). Amabile et al. also found that the risk of having low back pain in individuals with higher BMI was nearly twice that of those with lower BMI, and a dose–response relationship was identified between obesity and pain (11). Similarly, a cohort study in Great Britain found that BMI was associated with knee pain, with results that 19.1% of obese participants had pain symptoms (12). However, Sharon et al. performed a cross-sectional analysis including 142 subjects and revealed that no association was observed between joint pain and obesity (13). The reverse conclusions obtained by different studies might be caused by the limitations (small sample size, different races, and other existing confounders and bias) contained in the cohort and cross-sectional studies. Moreover, these studies could only find the correlation but not the causal relationship between obesity and musculoskeletal chronic pain.

Mendelian randomization (MR) is a method that takes single-nucleotide polymorphism (SNP) as instrumental variables (IVs) to explore the causal effect of exposures on outcomes (14). Due to the great development of large studies of genome-wide association studies (GWASs), MR has been a powerful alternative method for causal inference (15). Since genotypes appear before the occurrence of diseases and are largely unrelated to lifestyle or environmental factors after birth, MR could minimize the confounding factors and avoid reverse causality bias (16). This current study aims at detecting the causal effects of BMI, waist circumference (WC), hip circumference (HC), and waist-to-hip ratio (WHR) on knee, hip, and back pain with a two-sample MR method based on GWAS datasets.

Methods

Study design and genome-wide association study data source

The diagram of the study design for this MR analysis is shown in Figure 1. From the sketch map, we could learn that three assumptions are supposed to be met in the MR analysis. Firstly, the genetic variants should be closely related to the exposures. Secondly, the genetic variants are supposed to be independent of confounding factors. Lastly, the effects of the genetic variants on outcomes are only mediated by the exposures. As mentioned earlier, the SNPs were taken as IVs at a genome-wide significant level (p < 5 × 10−8) to perform the causal relationships between obesity and pain. Then linkage disequilibrium (LD) was calculated among the selected SNPs, and the SNPs with LD (while r2 > 0.001 or physical distance between them was within 10,000 kb) were excluded. Moreover, the F-statistics were also calculated, and weak IVs were removed (F < 10). The Steiger filtering test was also performed to exclude SNPs that explain more of the variance in the outcome rather than exposure.

FIGURE 1
www.frontiersin.org

Figure 1 The diagram of Mendelian randomization study. Three assumptions are required. Firstly, the genetic variants should be closely related to the exposures. Secondly, the genetic variants are supposed to be independent of confounding factors. Thirdly, the effects of the genetic variants on outcomes are only mediated by the exposures. SNPs, single-nucleotide polymorphisms.

The data used in the MR analyses originated from the GWAS summary data. Genetic data for BMI traits came from the GIANT consortium, which included 322,154 individuals and 2,554,668 SNPs of European ancestry (17). The genetic variants for waist circumference (N = 232,101), hip circumference (N = 213,038), and waist-to-hip ratio traits (N = 212,244) also originated from European population based on the GIANT consortium (18). Additionally, as for the outcomes, genetic variants for knee pain conducted by Ben et al. contained 98,704 cases and 363,153 controls from UK Biobank sources. The datasets for hip pain (52,087 patients and 409,770 controls) and back pain (118,471 cases and 343,386 controls) were also obtained from a GWAS meta-analysis of UK Biobank sources. All the genetic data were from patients with European ancestry and could be obtained from publicly available GWAS datasets (https://gwas.mrcieu.ac.uk). The MR-PRESSO test was applied to figure out and remove any potential outliers. Since the data of this study were based on existing publications and public databases, additional ethical approval or consent to participate was not required.

Statistical analyses

In this MR analysis, inverse variance weighted (IVW), weighted median, and MR-Egger were adopted to assess the causal effects of BMI, WC, HC, and WHR on knee pain, hip pain, and back pain. The IVW method analyzes each Wald ratio and provides a consistent estimate of the causal effect when all instrumental variables are valid, which is mainly used in the results with no heterogeneity or directional pleiotropy (19). The weighted median method was also used, for it could provide consistent estimates when up to 50% of the weight in the analysis originated from invalid instrumental variables (20). It should be noticed that the weighted median method would give a more accurate estimation than IVW when heterogeneity existed in the results. The MR-Egger method could identify and correct potential pleiotropy (p-value of intercept <0.05) and gives a consistent estimate (21). However, the effect size rather than the statistical significance of MR-Egger was concentrated in this current study, for the statistical power of MR-Egger was low (22). The results were reported as odds ratio (OR) and 95% confidence intervals (CIs), which provided an estimate of risk for outcome caused by each standard deviation (SD) increase in the risk factor.

Sensitivity testing

Then sensitivity analyses were conducted for the quality control of the MR results. Cochran’s Q test and I2 statistics were used to detect heterogeneity, and the intercept of the MR-Egger and MR-PRESSO tests were used to explore the directional pleiotropy. Additionally, this study performed a “leave-one-out” sensitivity test to evaluate whether the analysis was biased by a single SNP that had a particularly large horizontal pleiotropic effect. All statistical analyses were performed by the “Two-Sample MR” package (version 0.5.6) in R (version 4.1.2) software. The results were considered to be statistically significant when p < 0.05 and with heterogeneity when I2 > 50%.

Results

The extracted single-nucleotide polymorphisms for Mendelian randomization analyses

In this two-sample MR study, BMI, WC, HC, and WHR were taken as risk factors, while knee, hip, and back pain were taken as outcomes. The basic information of relevant SNPs that were selected for MR analyses is listed in Supplementary Material. A total of 69 SNPs with a mean of F = 66.78 associated with BMI and musculoskeletal chronic pain, 42 SNPs with a mean of F = 59.26 associated with WC and musculoskeletal chronic pain, 52 SNPs with a mean of F = 54.99 related to HC and musculoskeletal chronic pain, and 29 SNPs with a mean of F = 48.34 related to WHR and musculoskeletal chronic pain were included in this study.

Causal relationship of obesity on knee pain

The MR results for the causal effects of BMI, WC, HC, and WHR on knee pain are listed in Table 1, Figure 2, and Figures S1–S3. With one SD increase in BMI, the risk of knee pain was 1.049-fold as analyzed by the IVW method (95% CI: 1.034 to 1.063; p = 9.88 × 10−12). Similar results could also be acquired by weighted median (OR = 1.055; 95% CI: 1.040 to 1.070; p = 2.67 × 10−13) and MR-Egger (OR = 1.060; 95% CI: 1.018 to 1.104; p = 0.007) methods. Based on Cochran’s Q, I2, and MR-Egger intercept tests, there was heterogeneity (Q = 172.46, p = 3.04 × 10−11; I2 = 61.15%) but no pleiotropy (intercept beta = −3.00 × 10−4, p = 0.593) in the results. Then the “leave-one-out” sensitivity revealed that the causal relationship between BMI and knee pain was not changed by individual SNPs, meaning that the results were stable and reliable (Figure S1). The positive causal relationships between WC and knee pain were proved by the IVW (OR = 1.057; 95% CI: 1.041 to 1.072; p = 1.54 × 10−13), weighted median (OR = 1.064; 95% CI: 1.045 to 1.084; p = 2.51 × 10−11), and MR-Egger (OR = 1.079; 95% CI: 1.024 to 1.137; p = 0.007) methods. However, heterogeneity (Q = 72.40, p = 0.001) was observed in this causal relationship. No pleiotropy (intercept beta = −5.90 × 10−4, p = 0.416) was found in the results according to the MR-Egger intercept and MR-PRESSO tests. In the “leave-one-out” analysis, no single SNP strongly drove the overall effect of WC on knee pain (Figure S2).

TABLE 1
www.frontiersin.org

Table 1 The MR results regarding causal associations between obesity and chronic pain.

FIGURE 2
www.frontiersin.org

Figure 2 The Mendelian randomization analysis for causal effects of obesity on knee pain. (A) Forest plot about the causal effect of body mass index on knee pain. (B) Forest plot about the causal effect of waist circumference on knee pain. (C) Forest plot about the causal effect of hip circumference on knee pain. (D) Forest plot about the causal effect of waist-to-hip ratio on knee pain.

As for the HC trait, both IVW (OR = 1.034; 95% CI: 1.017 to 1.052; p = 1.32 × 10−4) and weighted median (OR = 1.029; 95% CI: 1.014 to 1.044; p = 1.34 × 10−4) methods showed positive causal effects of HC on knee pain. There was no pleiotropy (intercept beta = −6.3 × 10−4, p = 0.436), but heterogeneity (Q = 202.08, p = 7.22 × 10−20; I2 = 74.76%) was observed in the results based on Cochran’s Q, I2, and MR-Egger intercept tests. Then the “leave-one-out” analysis was conducted, revealing that the positive causal relationship was stable and reliable (Figure S3). The MR analysis for the causal relationship of WHR on knee pain was also conducted; however, no causal effect was found between them by IVW (OR = 1.003; 95% CI: 0.979 to 1.027; p = 0.822), weighted median (OR = 1.009; 95% CI: 0.984 to 1.034; p = 0.486), and MR-Egger (OR = 1.014; 95% CI: 0.904 to 1.136; p = 0. 816) methods.

Causal relationship of obesity on hip pain

Table 1, Figure 3, and Figures S4–S6 show the MR results of causal relationships between BMI, WC, HC, WHR, and hip pain. According to the results in Table 1, there was a positive causal effect of BMI on hip pain as analyzed by IVW (OR = 1.034; 95% CI: 1.024 to 1.044; p = 1.38 × 10−12), weighted median (OR = 1.034; 95% CI: 1.023 to 1.045; p = 1.60 × 10−9), and MR-Egger (OR = 1.039; 95% CI: 1.010 to 1.068; p = 0. 009) methods. Heterogeneity (Q = 132.53, p = 3.29 × 10−6) but not pleiotropy (intercept beta = −1.30 × 10−4, p = 0.739) was found in the MR results by Cochran’s Q and MR-Egger intercept tests. The “leave-one-out” analysis revealed that the result was stable and reliable, for no single SNP would influence the results (Figure S4). Additionally, with one SD higher than WC, the risk of hip pain increased by approximately 3.1% by the IVW (95% CI: 1.020 to 1.042; p = 2.61 × 10−8) and weighted median (95% CI: 1.018 to 1.045; p = 5.04 × 10−6) methods, while the risk increased by nearly 6.4% according to the MR-Egger (95% CI: 1.026 to 1.104; p = 0.002) method. However, heterogeneity (Q = 63.87, p = 0.010) existed in the results based on the Cochran’s Q test. The MR-Egger intercept and “leave-one-out” tests revealed that there was no pleiotropy (intercept beta = −9.00 × 10−4, p = 0.082) or single SNP that strongly drove the overall effect (Figure S5), which means that the result was robust.

FIGURE 3
www.frontiersin.org

Figure 3 The Mendelian randomization analysis for causal effects of obesity on hip pain. (A) Forest plot about the causal effect of body mass index on hip pain. (B) Forest plot about the causal effect of waist circumference on hip pain. (C) Forest plot about the causal effect of hip circumference on hip pain. (D) Forest plot about the causal effect of waist-to-hip ratio on hip pain.

Moreover, the IVW (OR = 1.027; 95% CI: 1.018 to 1.035; p = 5.48 × 10−10), weighted median (OR = 1.025; 95% CI: 1.013 to 1.038; p = 7.03 × 10−5), and MR-Egger (OR = 1.051; 95% CI: 1.026 to 1.077; p = 2.11 × 10−4) methods all showed a positive causal effect of HC on hip pain (Table 1). According to Cochran’s Q and MR-Egger intercept tests, there was heterogeneity (Q = 77.63, p = 0.010) but no pleiotropy (intercept beta = −7.5 × 10−4, p = 0.051) that existed in the MR result. The “leave-one-out” analysis illustrated that no single SNP strongly drove the overall effect between HC and hip pain (Figure S6). Additionally, the causal relationship between WHR and hip pain was also studied in this study. However, all the three methods yield a null causal relationship between them (IVW: OR = 1.003; 95% CI: 0.986 to 1.021; p = 0.717; weighted median: OR = 0.992; 95% CI: 0.974 to 1.011; p = 0.414; MR-Egger: OR = 1.027; 95% CI: 0.946 to 1.114; p = 0.533). No significant pleiotropy was detected in the analysis, and the “leave-one-out” test revealed a stable and reliable result.

Causal relationship of obesity on back pain

The MR results of causal relationships of BMI, WC, HC, and WHR on back pain are shown in Table 1 and Figure 4. With one SD increase in BMI trait, the incidence of back pain raised nearly 2% as analyzed by IVW (OR = 1.022; 95% CI: 1.007 to 1.038; p = 0.004) and weighted median (OR = 1.017; 95% CI: 1.002 to 1.032; p = 0.023) method. The results of Cochran’s Q and I2 tests showed the existence of heterogeneity (Q = 182.95, p = 1.06 × 10−12; I2 = 63.38). No pleiotropy (intercept beta = −1.7 × 10−4, p = 0.779) was observed in the MR analysis. In the “leave-one-out” analysis, no single SNP strongly drove the overall causal effect of BMI on back pain (Figure 4). As for the relationships between WC, HC, and back pain, causal effects were found among them by IVW (WC: OR = 1.027; 95% CI: 1.007 to 1.046; p = 0.006; HC: OR = 1.021; 95% CI: 1.006 to 1.036; p = 0.007). However, the weighted median (WC: OR = 1.011; 95% CI: 0.993 to 1.029; p = 0.234; HC: OR = 1.011; 95% CI: 0.995 to 1.027; p = 0.194) and MR-Egger (WC: OR = 1.030; 95% CI: 0.963 to 1.102; p = 0.397; HC: OR = 1.028; 95% CI: 0.982 to 1.077; p = 0.244) methods both showed a null causal effect of WC and HC on back pain. Due to the existence of heterogeneity but no pleiotropy, the weighted median method was adopted in the MR analysis; hence, the conclusions of causal effects among them were dropped. Additionally, the results revealed a null causal association between WHR and back pain by IVW (OR = 0.996; 95% CI: 0.975 to 1.018; p = 0.733). Weighted median and MR-Egger were consistent with the IVW method.

FIGURE 4
www.frontiersin.org

Figure 4 The Mendelian randomization analysis for causal effects of body mass index on back pain. (A) Scatter plot about the causal effect of body mass index on back pain. (B) Forest plot for the overall causal effects of body mass index on back pain. (C) Leave-one-out analysis for the causal effect of body mass index on back pain. (D) Funnel plot of SNPs related to body mass index and back pain. SNPs, single-nucleotide polymorphisms.

Discussion

This current research explored the causal relationships between obesity and musculoskeletal chronic pain based on the GWAS summary datasets, which might be the first study to investigate the causal effects of BMI, WC, HC, and WHR on knee, hip, and back pain by MR analysis. It revealed that BMI was positively related to knee, hip, and back pain and that WC and HC were positively associated with knee and hip pain, while WHR was not related to any type of musculoskeletal chronic pain.

Many previous studies have revealed the conclusion that being overweight or extremely obese could be associated with musculoskeletal chronic pain. In the United States, a survey concerning over 1 million individuals revealed that BMI was positively associated with daily pain and that obese people had a higher risk of being affected by daily pain (23). Another population-based longitudinal study containing 7,977 participants conducted by Rahman et al. in Finland found that central obesity defined by waist circumference was positively related to low back pain (OR = 1.40; 95% CI; 1.16 to 1.68), while general obesity measured with body mass index was positively associated with lumbar radicular pain (OR = 1.44; 95% CI: 1.12 to 1.85) (24). Similarly, Min et al. performed a research involving 40,999 individuals among community-dwelling older adults, showing that both overweight (OR = 1.166; 95% CI: 1.104 to 1.232, p < 0.01) and obesity (OR = 1.786; 95% CI: 1.530 to 2.085, p < 0.01) significantly contributed to the musculoskeletal chronic pain (25). Additionally, in another cross-sectional study with 6,524 elderly individuals in China, the researchers illustrated that the participants in overweight and obesity groups were more likely to be suffering from musculoskeletal chronic pain, compared with the normal-weight groups (26). The conclusions of these studies were consistent with our results that BMI and other obesity-related traits could causally increase musculoskeletal chronic pain.

On the contrary, there were also many scholars who found that obesity was not related to musculoskeletal chronic pain. An observational study on adolescent girls found that BMI was not associated with pain tolerance (27). Also, Paula et al. conducted a cross-sectional study containing 690 participants based in Brazil to investigate the relationship between obesity (evaluated by body mass index, bioimpedance, skinfold, arm, and abdominal circumference) and temporomandibular disorder pain in adolescents, and they discovered that all obesity-related traits were not related with temporomandibular disorder pain (28). Additionally, a matched pair study involving 1,128 female twins was performed to explore the relationships between BMI, percent body fat, WC, WHR, and low back pain, and no associations were observed between different obesity measures and low back pain after the full adjustment for genetic factors in this monozygotic within-pair case–control research (29). However, the different conclusions could be attributed to the diverse study designs and the existence of confounding factors as well as bias in observational studies. More importantly, these cross-sectional or cohort studies could only investigate the correlation but not the causal relationships between obesity and musculoskeletal chronic pain.

The association between obesity and musculoskeletal chronic pain could be explained by the mechanical loading and biochemical mechanisms (30). It was supposed that additional weight would give more pressure on the weight-bearing joints and skeletal muscle. Geoffrey et al. performed a study aiming at exploring different kinds of structural failures in intervertebral discs and found that chronic loadings would lead to a severely damaging impact on the spine, which could cause back pain (31). Similarly, another study was performed to detect the biomechanical stresses of the lower back through weight lifting tasks; it revealed that obese participants had a higher risk for back pain caused by the damage of lifting objects as compared with the normal-weight group (32). Moreover, inflammation also played a crucial part in uncovering the potential mechanism of the causal association between obesity and pain. It was widely known that obesity was considered a low-grade inflammatory disease and could produce many cytokines and adipokines through adipose tissue, which might be associated with musculoskeletal chronic pain (33). Some studies illustrated that leptin, one of the adipokines related to obesity, could increase the production of proteases and nitric oxide so as to cause low back pain (34). Another research revealed that galanin could regulate the pain threshold in obesity by central galanin receptor-1 and peripheral galanin receptor-2, though the antinociceptive effect of activating the receptors had not been clearly characterized (35). These theories could briefly explain the inner association between obesity musculoskeletal chronic pain, while more studies are need for the underlying mechanisms.

As for this Mendelian randomization analysis, there are many strengths that should be mentioned. This study investigated the causal relationships between obesity and knee, hip, and back pain with a large sample size based on GWAS datasets, which provided great persuasion for the associations. This current MR study could minimize the effects of bias and confounding factors due to the genetic variant alleles assigned randomly. Moreover, the reverse causality between obesity and musculoskeletal chronic pain that existed in the cross-sectional or cohort studies could also be avoided in this MR study, for the genetic variant exhibits earlier than diseases. Nevertheless, the limitations are also included in this study. Firstly, this MR analysis was based on patients of European ancestry, causing the obtained results to be boundedness, meaning that other races might have different results. Additionally, it was difficult for us to make further subgroup analyses for the causal effects of obesity on other sites or the degree of musculoskeletal chronic pain due to the restrictions of original data.

In conclusion, this present study found positive causal associations between BMI on knee pain, hip pain, and back pain. The waist circumference and hip circumference were also positively associated with knee pain and hip pain. However, no causal relationship was found between waist-to-hip ratio and knee, hip, and back pain. This research could provide much help for the potential mechanism and weight management in obese population with musculoskeletal chronic pain.

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

Since the data of this study were based on existing publications and public databases, additional ethical approval or consent to participate was not required.

Author contributions

XC performed the study and wrote the manuscript. HT and JL interpreted the results. RZ designed this study. All the authors contributed to the article and approved the submitted version.

Funding

The present study was supported by the Natural Science Foundation of Fujian Province (2022J011457) and Quanzhou Science and Technology Plan Project (2021N061S).

Acknowledgments

We thank the developers of the GWAS datasets.

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.971997/full#supplementary-material

Abbreviations

BMI, body mass index; MR, Mendelian randomization; SNP, single-nucleotide polymorphism; IVs, instrumental variables; GWAS, genome-wide association study; WC, waist circumference; HC, hip circumference; WHR, waist-to-hip ratio; LD, linkage disequilibrium; IVW, inverse variance weighted; OR, odds ratio; CIs, confidence intervals; SD, standard deviation.

References

1. González-Muniesa P, Mártinez-González M-A, Hu FB, Després J-P, Matsuzawa Y, Loos RJF, et al. Obesity. Nat Rev Dis Primers (2017) 3:17034. doi: 10.1038/nrdp.2017.34

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Morris A. Subtyping obesity. Nat Rev Endocrinol (2019) 15(6):316. doi: 10.1038/s41574-019-0201-y

CrossRef Full Text | Google Scholar

3. Litwin M, Kułaga Z. Obesity, metabolic syndrome, and primary hypertension. Pediatr Nephrol (2021) 36(4):825–37. doi: 10.1007/s00467-020-04579-3

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Al-Sulaiti H, Diboun I, Agha MV, Mohamed FFS, Atkin S, Dömling AS, et al. Metabolic signature of obesity-associated insulin resistance and type 2 diabetes. J Transl Med (2019) 17(1):348. doi: 10.1186/s12967-019-2096-8

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Liu Y, Ding W, Wang HL, Dai LL, Zong WH, Wang YZ, et al. Gut microbiota and obesity-associated osteoarthritis. Osteoarthritis Cartilage (2019) 27(9):1257–65. doi: 10.1016/j.joca.2019.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Cohen SP, Vase L, Hooten WM. Chronic pain: An update on burden, best practices, and new advances. Lancet (2021) 397(10289):2082–97. doi: 10.1016/S0140-6736(21)00393-7

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Murray CJL, Atkinson C, Bhalla K, Birbeck G, Burstein R, Chou D, et al. The state of us health, 1990-2010: Burden of diseases, injuries, and risk factors. JAMA (2013) 310(6):591–608. doi: 10.1001/jama.2013.13805

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Dahlhamer J, Lucas J, Zelaya C, Nahin R, Mackey S, DeBar L, et al. Prevalence of chronic pain and high-impact chronic pain among adults - united states, 2016. MMWR Morb Mortal Wkly Rep (2018) 67(36):1001–6. doi: 10.15585/mmwr.mm6736a2

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Steglitz J, Buscemi J, Ferguson MJ. The future of pain research, education, and treatment: A summary of the iom report "Relieving pain in America: A blueprint for transforming prevention, care, education, and research". Transl Behav Med (2012) 2(1):6–8. doi: 10.1007/s13142-012-0110-2

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Ozkuk K, Ates Z. The effect of obesity on pain and disability in chronic shoulder pain patients. J Back Musculoskelet Rehabil (2020) 33(1):73–9. doi: 10.3233/BMR-181384

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Dario AB, Ferreira ML, Refshauge KM, Lima TS, Ordonana JR, Ferreira PH. The relationship between obesity, low back pain, and lumbar disc degeneration when genetics and the environment are considered: A systematic review of twin studies. Spine J (2015) 15(5):1106–17. doi: 10.1016/j.spinee.2015.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Macfarlane GJ, de Silva V, Jones GT. The relationship between body mass index across the life course and knee pain in adulthood: Results from the 1958 birth cohort study. Rheumatol (Oxford) (2011) 50(12):2251–6. doi: 10.1093/rheumatology/ker276

CrossRef Full Text | Google Scholar

13. Bout-Tabaku S, Klieger SB, Wrotniak BH, Sherry DD, Zemel BS, Stettler N. Adolescent obesity, joint pain, and hypermobility. Pediatr Rheumatol Online J (2014) 12:11. doi: 10.1186/1546-0096-12-11

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol (2013) 37(7):658–65. doi: 10.1002/gepi.21758

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ (2018) 362:k601. doi: 10.1136/bmj.k601

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lawlor DA, Harbord RM, Sterne JAC, Timpson N, Davey Smith G. Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Stat Med (2008) 27(8):1133–63. doi: 10.1002/sim.3034

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature (2015) 518(7538):197–206. doi: 10.1038/nature14177

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Shungin D, Winkler TW, Croteau-Chonka DC, Ferreira T, Locke AE, Mägi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature (2015) 518(7538):187–96. doi: 10.1038/nature14132

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hartwig FP, Davies NM, Hemani G, Davey Smith G. Two-sample mendelian randomization: Avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int J Epidemiol (2016) 45(6):1717–26. doi: 10.1093/ije/dyx028

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan NA, Thompson JR. Assessing the suitability of summary data for two-sample mendelian randomization analyses using Mr-egger regression: The role of the I2 statistic. Int J Epidemiol (2016) 45(6):1961–74. doi: 10.1093/ije/dyw220

PubMed Abstract | CrossRef Full Text | Google Scholar

21. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

22. 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(2):512–25. doi: 10.1093/ije/dyv080

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Stone AA, Broderick JE. Obesity and pain are associated in the united states. Obes (Silver Spring) (2012) 20(7):1491–5. doi: 10.1038/oby.2011.397

CrossRef Full Text | Google Scholar

24. Shiri R, Falah-Hassani K, Heliovaara M, Solovieva S, Amiri S, Lallukka T, et al. Risk factors for low back pain: A population-based longitudinal study. Arthritis Care Res (Hoboken) (2019) 71(2):290–9. doi: 10.1002/acr.23710

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Qian M, Shi Y, Yu M. The association between obesity and chronic pain among community-dwelling older adults: A systematic review and meta-analysis. Geriatr Nurs (2021) 42(1):8–15. doi: 10.1016/j.gerinurse.2020.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Li J, Chen J, Qin Q, Zhao D, Dong B, Ren Q, et al. Chronic pain and its association with obesity among older adults in China. Arch Gerontol Geriatr (2018) 76:12–8. doi: 10.1016/j.archger.2018.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Goldschmidt AB, Smith KE, Bond DS, Hipwell AE, Stepp SD, Keenan K. Pain, eating behavior, and weight outcomes in adolescent girls. Int J Eat Disord (2021) 54(5):893–7. doi: 10.1002/eat.23484

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Jordani PC, Campi LB, Braido GVV, Fernandes G, Visscher CM, Gonçalves DAG. Obesity, sedentarism and tmd-pain in adolescents. J Oral Rehabil (2019) 46(5):460–7. doi: 10.1111/joor.12771

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Dario AB, Ferreira ML, Refshauge K, Sánchez-Romera JF, Luque-Suarez A, Hopper JL, et al. Are obesity and body fat distribution associated with low back pain in women? a population-based study of 1128 Spanish twins. Eur Spine J (2016) 25(4):1188–95. doi: 10.1007/s00586-015-4055-2

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Chin SH, Huang WL, Akter S, Binks M. Obesity and pain: A systematic review. Int J Obes (Lond) (2020) 44(5):969–79. doi: 10.1038/s41366-019-0505-y

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Desmoulin GT, Pradhan V, Milner TE. Mechanical aspects of intervertebral disc injury and implications on biomechanics. Spine (2020) 45(8):E457–E64. doi: 10.1097/BRS.0000000000003291

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Singh D, Park W, Hwang D, Levy MS. Severe obesity effect on low back biomechanical stress of manual load lifting. Work (2015) 51(2):337–48. doi: 10.3233/WOR-141945

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Ruiz-Fernández C, Francisco V, Pino J, Mera A, González-Gay MA, Gómez R, et al. Molecular relationships among obesity, inflammation and intervertebra disc degeneration: Are adipokines the common link? Int J Mol Sci (2019) 20(8):1–7. doi: 10.3390/ijms20082030

CrossRef Full Text | Google Scholar

34. Segar AH, Fairbank JCT, Urban J. Leptin and the intervertebral disc: A biochemical link exists between obesity, intervertebral disc degeneration and low back pain-an in vitro study in a bovine model. Eur Spine J (2019) 28(2):214–23. doi: 10.1007/s00586-018-5778-7

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Yu M, Fang P, Shi M, Zhu Y, Sun Y, Li Q, et al. Galanin receptors possibly modulate the obesity-induced change in pain threshold. Peptides (2013) 44:55–9. doi: 10.1016/j.peptides.2013.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: body mass index, waist circumference, hip circumference, waist-to-hip ratio, Mendelian randomization, chronic pain

Citation: Chen X, Tang H, Lin J and Zeng R (2022) Causal relationships of obesity on musculoskeletal chronic pain: A two-sample Mendelian randomization study. Front. Endocrinol. 13:971997. doi: 10.3389/fendo.2022.971997

Received: 17 June 2022; Accepted: 29 July 2022;
Published: 23 August 2022.

Edited by:

Luca Busetto, Università degli Studi di Padova, Italy

Reviewed by:

Xiaoqian Dang, The Second Affiliated Hospital of Xi’an Jiaotong University, China
Yi Lin, Third Affiliated Hospital of Chongqing Medical University, China

Copyright © 2022 Chen, Tang, Lin and Zeng. 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: Rongdong Zeng, zrongdong@163.com

Disclaimer: 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.