Corrigendum: Assessing the causal relationship between genetically determined inflammatory biomarkers and low back pain risk: a bidirectional two-sample Mendelian randomization study
- 1Tui-Na Department, Hospital of Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
- 2Department of Acupuncture, Moxibustion, Tui-Na and Rehabilitation, Guang'an City Hospital of Traditional Chinese Medicine, Guangan, Sichuan, China
- 3Tui-Na Teaching and Research Department, College of Acupuncture and Tuina, Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
- 4Tui-Na Department, Meishan City Hospital of Traditional Chinese Medicine, Meishan, Sichuan, China
- 5Rehabilitation Department, School of Clinic Medicine & The First Affiliated Hospital of Chengdu Medical College, Chengdu, Sichuan, China
- 6Chinese Medicine Rehabilitation Department, Jiahekang Hospital, Luzhou, Sichuan, China
- 7Department of Acupuncture, Moxibustion, Tui-Na and Rehabilitation, The Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University, Luzhou, Sichuan, China
- 8Center for Traditional Chinese Medicine Prevention and Health Care, Chengdu Integrated TCM & Western Medicine Hospital, Chengdu, Sichuan, China
Background: Observational studies have suggested an association between inflammatory markers and low back pain (LBP), but the causal relationship between these factors remains uncertain.
Methods: We conducted a bidirectional two-sample Mendelian randomization analysis (MR) study to investigate whether there is a causal relationship between inflammatory markers and low back pain. We obtained genetic data for CRP, along with its upstream inflammatory markers IL-6, IL-8, and IL-10, as well as low back pain from publicly available genome-wide association studies (GWAS). We applied several MR methods, including inverse variance weighting, weighted median, MR-Egger, Wald Ratio, and MR-PRESSO, to test for causal relationships. Sensitivity analyses were also conducted to assess the robustness of the results.
Results: Our analyses utilizing the Inverse Variance Weighted (IVW) method, the MR-Egger method, and the weighted median method indicated that IL-6 may be associated with an increased risk of LBP (Effect Size: -0.009, 95% Confidence Interval: -0.013–0.006, p = 9.16e-08); however, in the reverse direction, there was no significant causal effect of LBP on inflammatory markers.
Conclusion: Our study used a Mendelian randomization approach and found that elevated IL-6 levels may reduce the risk of LBP.
1 Introduction
Low back pain (LBP) is a prevalent musculoskeletal disorder that has a profound impact on quality of life. In fact, it has been identified as the leading cause of disability and work-related absenteeism in 126 countries (1). This problem has become a significant public health issue, and it has a considerable economic burden on society (2, 3). In the United States alone, the direct and indirect costs associated with LBP surpass $100 billion annually, including medical expenses, lost wages, and decreased productivity (4). As the aging population continues to grow globally, the economic burden of LBP is expected to increase, underscoring the need to identify the risk factors for its development.
Inflammation is believed to be a key mechanism underlying LBP and spinal degeneration (3, 5). C-reactive protein (CRP) is a well-known biomarker of systemic inflammation, and previous observational studies have shown a positive correlation between its levels and the severity of LBP (6–8). Pro-inflammatory cytokines, such as interleukins, trigger the production of inflammatory markers like CRP in the liver (9). Several observational studies have found that patients with LBP have higher levels of serum IL-6 and IL-8 protein and mRNA compared to healthy controls (10–12). Additionally, IL-10 is considered an anti-inflammatory marker that inhibits the synthesis of pro-inflammatory cytokines. Observational studies have shown an inverse association between IL-10 and LBP pain severity (10, 12). However, conventional observational studies have limitations, and the association between inflammatory markers and LBP reported in previous studies may still be explained by reverse causality and residual confounding (13, 14). Therefore, the causal relationship between inflammatory markers and LBP remains uncertain.
Mendelian randomization (MR) is a widely used analytical method that employs genetic variation as an instrumental variable to produce more reliable causal estimates of the effect of risk factors on long-term exposure to disease outcomes (14, 15). The rationale for relying on MR as a more powerful method of causal inference than traditional observational studies stems from Mendel’s law and the fact that genotypes for germline genetic variation are determined at the time of conception and are often unrelated to the traditional confounders of observational studies. Previous studies have demonstrated the utility of genetic tools to elucidate the causal relationship between inflammatory markers and disease risk (16). Therefore, in this bidirectional two-sample Mendelian randomization analysis study, our aim was to investigate the causal relationship between CRP, as well as its upstream inflammatory markers IL-6, IL-8, and IL-10, and the risk of developing LBP.
2 Methods
2.1 Study design
Mendelian randomization is built upon three main assumptions. First, genetic variants are associated with the exposure. Second, there should be no association between genetic variants and confounders. Third, genetic variants influence the risk of the outcome only through the exposure, not through other pathways (17). In the present study, we employed bidirectional MR analysis to comprehensively infer the causal relationship between inflammatory markers and low back pain in both the forward and reverse directions. We have summarized our study design in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1174656/fimmu-14-1174656-HTML-r1/image_m/fimmu-14-1174656-g001.jpg)
Figure 1 Overview of the assumptions of the Mendelian randomization (MR) design and the study design. The MR design reduces residual confounding and reverse causality, thereby strengthening causal inferences about exposure-outcome associations. This is because genetic variants used as instrumental variables to study altered exposure effects are randomly assigned at the time of conception and are therefore not susceptible to confounding by environmental factors or reverse causality.
2.2 Data sources
The variants associated with CRP were obtained from a genome-wide association study (GWAS) of up to 204,402 individuals of European ancestry (18). The variants associated with IL-6 were obtained from a genome-wide association study (GWAS) of up to 67,482 individuals of European ancestry (19). Variants associated with IL-8 were derived from a GWAS of up to 21,758 individuals of European ancestry from the SCALLOP consortium (20). Genetic variants for IL-10 were collected from a GWAS of 8,293 Finnish individuals (21). The genetic association with back pain was obtained from the UK Biobank abstract data and can be found by searching for GWAS ID: ukb-b-9838 on the web at https://gwas.mrcieu.ac.uk/datasets/. The identified genetic variants were found to be associated with the self-reported experience of back pain within the past month. The phenotype was derived from responses to a touchscreen questionnaire, which asked participants: “In the last month, have you experienced any of the following symptoms that interfered with your usual activities? (You may select more than one answer).” The response options included: Headache, Facial pain, Neck or shoulder pain, Back pain, Stomach or abdominal pain, Hip pain, Knee pain, and Pain all over the body. Additional validation checks were performed to ensure the accuracy of the phenotype. Specifically, if a participant selected “Pain all over the body,” “Knee pain,” or “Neck or shoulder pain,” no additional response options were permitted. This GWAS includes 461,857 unrelated individuals of European ancestry, comprising 118,471 cases and 343,386 controls. To reduce selection bias and improve the robustness of the analysis, all datasets used in our study comprised individuals of European ancestry. All data used in this study were obtained from GWAS, with prior ethical approval and patient consent. The study protocol complied with the guidelines of the Declaration of Helsinki and was approved by the ethics committees of all participating sites. The details of these GWAS data sources are listed in Supplementary Table S1. There is no overlap between the exposure GWASs and the outcome GWASs.
2.3 Selection of genetic instruments
To identify SNPs for MR analysis, we selected SNPs associated with the exposure under genomewide significance threshold (p < 5x10^-8) as instrumental variables (IVs). These SNPs were further screened using a LD distance threshold of 10,000 kb and r2 < 0.001 to ensure independence between genetic variants (22). We calculated the F-statistic for each IV to assess the strength of its association with the exposure (23) and included IVs with an F-statistic > 10 (24). Subsequently, the palindromic SNPs were removed to ensure that their association effects were linked with the same alleles in both the exposure and outcome.
2.4 Statistical power calculation
We sought to assess the statistical power of our MR analyses through the use of an online web tool (https://sb452.shinyapps.io/power) (25). The assessment of statistical power for MR analyses was based on several parameters, including the total sample size, the significance level of 0.05/8, the proportion of variance (R2) in the exposure explained by instrumental variables, and the ratio of cases to controls.
25 Statistical analyses
To test for causal effects of exposure on outcome, we performed MR analyses using five methods, namely inverse variance weighting (IVW) (26), weighted median (27), MR-Egger (28), Wald Ratio(for IL-8, because only one SNP was available), and MR pleiotropic residuals and outliers (MR-PRESSO) (29). These five methods make different assumptions and use different strategies to deal with IVs with horizontal pleiotropy effects. The IVW method relies on the assumption that there is no pleiotropy and assumes that all SNPs are valid genetic tools (26). The weighted median method assumes that at least 50% of the IVs are valid (27). The MR-Egger method provides causal estimates even if all IVs are invalid (28). The MR-PRESSO method detects possible IV outliers by global testing and provides unbiased causal estimates by removing the identified outliers (29).
Heterogeneity in causal estimates among instrumental variables indicates a potential violation of the assumptions of MR analysis. The Cochran’s Q test was used to examine the heterogeneity in causal estimates, and we used both the causal estimates of fixed effects IVW method and MR-Egger regression to detect heterogeneity. The heterogeneities were quantified using Cochran’s Q statistics and a P-value smaller than 0.05 was considered significant heterogeneity. To assess the potential pleiotropic effects of instrumental variables, the MR-Egger regression was used. The directional horizontal pleiotropy in the causal estimates may be indicated by the intercept term in MR-Egger regression. Additionally, we performed a leave-one-out analysis where we excluded each SNP in turn and then ran MR analysis on the remaining SNPs in order to detect potentially outlying instrumental variables.
All results are expressed as effect size (ES) and corresponding 95% confidence intervals (CI). All statistical analyses were two-sided. A p-value was considered statistically significant when less than 0.006 (0.05/8 adjusted with the Bonferroni method) and was considered suggestively significant between 0.006 and 0.05. All analyses were performed using the TwoSampleMR and MR-PRESSO packages in R version 4.2.2.
3 Result
3.1 Causal effects of different inflammatory biomarkers on the risk of low back pain
To assess the impact of inflammatory markers on the risk of low back pain, we conducted an MR analysis. Initially, we identified 50, 1, 94, and 2 SNPs associated with CRP, IL-8, IL-6, and IL-10, respectively. The characteristics of these IVs are shown in Supplementary Table S2–S5. In addition, we calculated the statistical efficacy of each exposure in each cohort (CRP, 100%; IL-8, 83.3%; IL-6, 57.1%; IL-10, 4.3%) (Supplementary Table S6).
We used the IVW and MR-Egger methods to evaluate the effect of different inflammatory markers on low back pain, and the results showed IL-6 may be associated with an increased risk of LBP (ES: -0.009, 95% CI: -0.013, -0.006, p = 9.164e-08), and no significant causal relationship between other inflammatory markers on low back pain (Figure 2; Supplementary Table S6). Furthermore, Cochran’s Q test showed no heterogeneity (CRP-LBP, Q =131.904, p = 1.141; IL-6 -LBP, Q =47.016, p = 0.999; IL-10-LBP, Q = 0.215, p = 0.642). We also conducted a horizontal pleiotropy test, which indicated no directional pleiotropy (CRP-LBP, intercept = 0.0002, p = 0.480; IL-6 -LBP, intercept = 0.0001, p = 0.694). The MR-PRESSO global test also showed no horizontal pleiotropy effects (CRP-LBP, p = 0.159; IL-6-LBP, p = 0.500). To ensure that the results were not influenced by a single SNP, we conducted a leave-one-out sensitivity test for IL-6, which demonstrated that the causal effect of IL-6 on LBP did not fluctuate significantly with the absence of any single SNP (Supplementary Figure S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1174656/fimmu-14-1174656-HTML-r1/image_m/fimmu-14-1174656-g002.jpg)
Figure 2 Effect of different inflammatory markers on LBP. IL-10 did not have enough SNP to complete the weighted median, MR-Egger and MR-PRESSO. MR, Mendelian randomization; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier test; ES, effect size; CI, confidence interval.
3.2 Causal impact of low back pain on different inflammatory markers
For the MR analysis of low back pain with different inflammatory markers, we obtained SNPs associated with low back pain from the GWAS (Supplementary Tables S7-S10). We performed IVW analysis, MR-Egger method, weighted median method, and weighted mode method, but found no evidence of a causal effect of low back pain on inflammatory markers (p > 0.05) (Figure 3). Cochran’s Q test showed no heterogeneity (LBP-CRP, Q = 19.264, p = 0.115; LBP-IL6, Q = 8.536, p = 0.969; LBP-IL-10, Q = 9.637, p = 0.723; LBP-IL8, Q = 20.729, p = 0.293). All Egger regression tests were negative (LBP-CRP, intercept = 0.011, p = 0.103; LBP-IL6rα, intercept = -0.001, p = 0.919; LBP-IL10, intercept = 0.030, p = 0.198; LBP-IL8, intercept = 0.014, p = 0.273). The MR-PRESSO global test also showed no horizontal pleiotropy effects (LBP-CRP, p = 0.933; LBP-IL6rα, p = 0.360; LBP-IL10, p = 0.214; LBP-IL8, p = 0.422), indicating that our MR results were not influenced by horizontal pleiotropy. Since there are no positive results, leave-one-out sensitivity test are inapplicable here. In addition, we calculated the statistical efficacy for each outcome in each cohort. The results showed that the statistical efficacy of each exposure was 100%, thus confirming the robustness of our MR analysis. (Supplementary Table S6)
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1174656/fimmu-14-1174656-HTML-r1/image_m/fimmu-14-1174656-g003.jpg)
Figure 3 Effect of LBP on different inflammatory markers. MR, Mendelian randomization; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier test; ES, effect size; CI, confidence interval.
4 Discussion
In this study, we employed bidirectional MR analysis to investigate a potential causal relationship between different inflammatory factors and low back pain. Using the largest publicly available GWAS dataset, we found that lower IL-6 level was associated with higher risk of LBP; however, in the reverse direction, there was no significant causal effect of LBP on inflammatory markers.
Previous studies have demonstrated a strong association between them. For instance, CRP was the first inflammatory marker reported to be associated with low back pain. As early as 2005, an observational study found that CRP was linked to pain levels in patients with acute low back pain (6). IL-6 promotes the differentiation of monocytes into macrophages and activates lymphocyte maturation to help mediate the acute phase response to injury (30). Queiroz et al. found that IL-6 was associated with both the severity of LBP and disability due to LBP (31). A systematic review that included six studies found evidence of elevated CRP in patients with acute LBP, but not elevated IL-6 in patients with LBP (32). Another systematic review that included six studies found an association between elevated CRP and IL-6 and low back pain, supporting the idea that there is a positive correlation between inflammatory markers and low back pain (33). Wang et al. found higher serum IL-6 and IL-8 levels in patients with severe LBP than in mild and control subjects (12). Conversely, IL-10 is an anti-inflammatory marker, and Wang et al. found that patients with mild LBP had lower serum IL-10 levels than those with severe LBP (12). The levels of inflammatory markers were influenced by the degree of low back pain and time to relief, which may be the reason why this MR study did not find a causal effect of LBP on inflammatory biomarkers.
Previous studies have identified possible mechanisms by which inflammatory markers are involved in the development of LBP. Lumbar disc degeneration is an important cause of low back pain (34). Park et al. found that serum levels of IL-6 and IL-8 were higher in patients with lumbar disc herniation than in controls (35). Pedersen et al. found that chronic lumbar radicular pain may be associated with a sustained increase in serum pro-inflammatory substances IL-6 and IL-8 after disc herniation surgery (36). These inflammatory biomarkers are thought to promote matrix degradation, chemokine production, and cellular phenotypic changes. The resulting imbalance between catabolic and anabolic reactions leads to lumbar disc degeneration, disc herniation and nerve root pain (3).
Contrary to the findings of many previous observational studies, our results suggest that elevated levels of IL-6 may be associated with a decreased risk of LBP. This may be attributed to the potential role of an active inflammatory response in reducing the incidence of chronic low back pain. For instance, Marc et al. conducted a transcriptome-wide analysis of peripheral immune cells in 98 subjects with acute LBP and observed thousands of dynamic transcriptional changes over a 3-month period in subjects who experienced pain relief, but not in those with persistent pain (37). Their findings suggest that neutrophil-driven upregulation of the inflammatory response may prevent the transition to chronic pain in LBP patients. In a mouse pain assay, they also found that early treatment with steroids or non-steroidal anti-inflammatory drugs (NSAIDs) resulted in prolonged pain, whereas no such prolongation was observed with other analgesics. Furthermore, their analysis of pain trajectories in human subjects reporting acute back pain in the UK Biobank revealed an increased risk of pain persistence in subjects taking NSAIDs. Our study supports the findings of Marc et al., indicating that an active inflammatory response may prevent the transition to chronic low back pain in LBP patients. Additionally, the observed association between elevated IL-6 levels and a reduced risk of LBP may be explained by the fact that lumbar disc herniation is the most common cause of LBP. Weber et al. reported that serum IL-6 levels were significantly higher in patients with lumbar disc degeneration and spinal stenosis than in those with lumbar disc herniation (38).
To our knowledge, no MR studies have been reported on the causal effects of inflammatory markers on LBP or vice versa. Our study utilized multiple IVs from GWAS of inflammatory markers and LBP to increase statistical power for detecting causality, providing a more precise assessment of effect sizes.
However, our study has several limitations. Firstly, the limitations of Mendelian Randomization (MR) analysis preclude a thorough examination of the second and third assumptions and may introduce bias. For instance, several factors that contribute to low back pain may also contribute to inflammation, including mechanical trauma, obesity, and infection. Secondly, the population included in this MR analysis is of European descent and it remains to be determined whether the results can be replicated in Asian and African populations. Thirdly, we did not stratify the causal relationship between low back pain and inflammatory markers according to the duration of low back pain, despite evidence suggesting that the duration of low back pain can influence the levels of inflammatory markers in the blood. Fourth, it is important to note that the genetic association between IL-6 and IL-10 and LBP was based on a relatively small GWAS, resulting in limited statistical power (57.1% for IL6-LBP and 4.3% for IL10-LBP). Therefore, a subsequent GWAS with a larger sample size is necessary to confirm and update the findings of this study.
5 Conclusion
Our study used a Mendelian randomization approach and found that elevated IL-6 levels may reduce the risk of LBP.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials. Further inquiries can be directed to the corresponding authors.
Author contributions
Conceptualization, WL, ZL and CS; Methodology, HL, QL and CS; Formal analysis, WL, WH and JQ; Software, WL, JL and QL; Data Curation, WL and CS; Investigation, JL, YF and CL; Validation, CL, and YF; Resources, WL and YF; Visualization, WL and HL; Supervision, HL and ZL; Writing-original draft, WL, CS and QL; Writing-review and editing, CS and CL; Funding acquisition, CS; Project administration, WL and CS. All authors contributed to the article and approved the submitted version.
Funding
This project was supported by National Natural Science Foundation of China (No. 82004497) and China Postdoctoral Science Foundation (No. 2021M693788).
Acknowledgments
We thank all participants and investigators who contributed to the GWAS data used in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1174656/full#supplementary-material
References
1. Disease GBD, Injury I, Prevalence C. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet (2018) 392:1789–858. doi: 10.1016/S0140-6736(18)32279-7
2. Daly C, Ghosh P, Jenkin G, Oehme D, Goldschlager TA. Review of animal models of intervertebral disc degeneration: pathophysiology, regeneration, and translation to the clinic. BioMed Res Int (2016) 2016:5952165. doi: 10.1155/2016/5952165
3. Risbud MV, Shapiro IM. Role of cytokines in intervertebral disc degeneration: pain and disc content. Nat Rev Rheumatol (2014) 10(1):44–56. doi: 10.1038/nrrheum.2013.160
4. Katz JN. Lumbar disc disorders and low-back pain: socioeconomic factors and consequences. J Bone Joint Surg Am (2006) 88(Suppl 2):21–4. doi: 10.2106/JBJS.E.01273
5. Wuertz K, Haglund L. Inflammatory mediators in intervertebral disk degeneration and discogenic pain. Global Spine J (2013) 3(3):175–84. doi: 10.1055/s-0033-1347299
6. Stürmer T, Raum E, Buchner M, Gebhardt K, Schiltenwolf M, Richter W, et al. Pain and high sensitivity C reactive protein in patients with chronic low back pain and acute sciatic pain. Ann Rheum Dis (2005) 64(6):921–5. doi: 10.1136/ard.2004.027045
7. Gebhardt K, Brenner H, Stürmer T, Raum E, Richter W, Schiltenwolf M, et al. The course of high-sensitive C-reactive protein in correlation with pain and clinical function in patients with acute lumbosciatic pain and chronic low back pain - a 6 months prospective longitudinal study. Eur J Pain (2006) 10(8):711–9. doi: 10.1016/j.ejpain.2005.11.005
8. Rannou F, Ouanes W, Boutron I, Lovisi B, Fayad F, Macé Y, et al. High-sensitivity C-reactive protein in chronic low back pain with vertebral end-plate modic signal changes. Arthritis Rheumatol (2007) 57(7):1311–5. doi: 10.1002/art.22985
9. Khan AN, Jacobsen HE, Khan J, Filippi CG, Levine M, Lehman RA Jr, et al. Inflammatory biomarkers of low back pain and disc degeneration: a review. Ann N Y Acad Sci (2017) 1410(1):68–84. doi: 10.1111/nyas.13551
10. Uçeyler N, Rogausch JP, Toyka KV, Sommer C. Differential expression of cytokines in painful and painless neuropathies. Neurology (2007) 69(1):42–9. doi: 10.1212/01.wnl.0000265062.92340.a5
11. Licciardone JC, Kearns CM, Hodge LM, Bergamini MV. Associations of cytokine concentrations with key osteopathic lesions and clinical outcomes in patients with nonspecific chronic low back pain: results from the OSTEOPATHIC trial [published correction appears in J Am Osteopath Assoc. 2017 Jun 1;117(6):350]. J Am Osteopath Assoc (2012) 112(9):596–605. doi: 10.7556/jaoa.2012.112.9.596
12. Wang K, Bao JP, Yang S, Hong X, Liu L, Xie XH, et al. A cohort study comparing the serum levels of pro- or anti-inflammatory cytokines in patients with lumbar radicular pain and healthy subjects. Eur Spine J (2016) 25(5):1428–34. doi: 10.1007/s00586-015-4349-4
13. Phillips AN, Smith GD. How independent are "independent" effects? Relative risk estimation when correlated exposures are measured imprecisely. J Clin Epidemiol (1991) 44(11):1223–31. doi: 10.1016/0895-4356(91)90155-3
14. Smith GD, Ebrahim S. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol (2003) 32(1):1–22. doi: 10.1093/ije/dyg070
15. Smith GD, Lawlor DA, Harbord R, Timpson N, Day I, Ebrahim S. Clustered environments and randomized genes: a fundamental distinction between conventional and genetic epidemiology. PloS Med (2007) 4(12):e352. doi: 10.1371/journal.pmed.0040352
16. Hartwig FP, Borges MC, Horta BL, Bowden J, Davey Smith G. Inflammatory biomarkers and risk of schizophrenia: a 2-sample Mendelian randomization study. JAMA Psychiatry (2017) 74(12):1226–33. doi: 10.1001/jamapsychiatry.2017.3191
17. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG. EPIC- InterAct Consortium. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol (2015) 30(7):543–52. doi: 10.1007/s10654-015-0011-z
18. Ligthart S, Vaez A, Võsa U, Stathopoulou MG, de Vries PS. Genome analyses of >200,000 individuals identify 58 loci for chronic inflammation and highlight pathways that link inflammation and complex disorders. Am J Hum Genet (2018) 103(5):691–706. doi: 10.1016/j.ajhg.2018.09.009
19. Ahluwalia TS, Prins BP, Abdollahi M, Armstrong NJ, Aslibekyan S, Bain L, et al. Genome-wide association study of circulating interleukin 6 levels identifies novel loci. Hum Mol Genet (2021) 30(5):393–409. doi: 10.1093/hmg/ddab023
20. Folkersen L, Gustafsson S, Wang Q, Hansen DH, Hedman ÅK, Schork A, et al. Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals. Nat Metab (2020) 2(10):1135–48. doi: 10.1038/s42255-020-00287-2
21. Ahola-Olli AV, Würtz 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(1):40–50. doi: 10.1016/j.ajhg.2016.11.007
22. 1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, et al. A global reference for human genetic variation. Nature (2015) 526(7571):68–74. doi: 10.1038/nature15393
23. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol (2011) 40(3):740–52. doi: 10.1093/ije/dyq151
24. Brion MJ, Shakhbazov K, Visscher PM. Calculating statistical power in Mendelian randomization studies. Int J Epidemiol (2013) 42(5):1497–501. doi: 10.1093/ije/dyt179
25. Burgess S. Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol (2014) 43(3):922–9. doi: 10.1093/ije/dyu005
26. 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
27. 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(4):304–14. doi: 10.1002/gepi.21965
28. 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
29. 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 [published correction appears in Nat Genet. 2018 Aug;50(8):1196]. Nat Genet (2018) 50(5):693–8. doi: 10.1038/s41588-018-0099-7
30. 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 Rheumatol (2010) 62(7):1974–82. doi: 10.1002/art.27444
31. de Queiroz BZ, Pereira DS, Lopes RA, Felício DC, Silva JP, de Britto Rosa NM, et al. Association between the plasma levels of mediators of inflammation with pain and disability in the elderly with acute low back pain: data from the back complaints in the elders (BACE)-Brazil study. Spine (Phila Pa 1976) (2016) 41(3):197–203. doi: 10.1097/BRS.0000000000001214
32. Morris P, Ali K, Merritt M, Pelletier J, Macedo LG. A systematic review of the role of inflammatory biomarkers in acute, subacute and chronic non-specific low back pain. BMC Musculoskelet Disord (2020) 21(1):142. doi: 10.1186/s12891-020-3154-3
33. Lim YZ, Wang Y, Cicuttini FM, Hughes HJ, Chou L, Urquhart DM, et al. Association between inflammatory biomarkers and nonspecific low back pain: a systematic review. Clin J Pain (2020) 36(5):379–89. doi: 10.1097/AJP.0000000000000810
34. Moore RJ, Vernon-Roberts B, Fraser RD, Osti OL, Schembri M. The origin and fate of herniated lumbar intervertebral disc tissue. Spine (Phila Pa 1976) (1996) 21(18):2149–55. doi: 10.1097/00007632-199609150-00018
35. Park MS, Lee HM, Hahn SB, Moon SH, Kim YT, Lee CS, et al. The association of the activation-inducible tumor necrosis factor receptor and ligand with lumbar disc herniation. Yonsei Med J (2007) 48(5):839–46. doi: 10.3349/ymj.2007.48.5.839
36. Pedersen LM, Schistad E, Jacobsen LM, Røe C, Gjerstad J. Serum levels of the pro-inflammatory interleukins 6 (IL-6) and -8 (IL-8) in patients with lumbar radicular pain due to disc herniation: a 12-month prospective study. Brain Behav Immun (2015) 46:132–6. doi: 10.1016/j.bbi.2015.01.008
37. Parisien M, Lima LV, Dagostino C, El-Hachem N, Drury GL, Grant AV, et al. Acute inflammatory response via neutrophil activation protects against the development of chronic pain. Sci Transl Med (2022) 14(644):eabj9954. doi: 10.1126/scitranslmed.abj9954
Keywords: inflammatory biomarkers, low back pain, Mendelian randomization, causality, interleukin 6 (IL-6)
Citation: Li W, Lu Q, Qian J, Feng Y, Luo J, Luo C, He W, Dong B, Liu H, Liu Z and Su C (2023) Assessing the causal relationship between genetically determined inflammatory biomarkers and low back pain risk: a bidirectional two-sample Mendelian randomization study. Front. Immunol. 14:1174656. doi: 10.3389/fimmu.2023.1174656
Received: 01 March 2023; Accepted: 27 June 2023;
Published: 13 July 2023.
Edited by:
Maxim B. Freidin, Queen Mary University of London, United KingdomReviewed by:
Ian B. Stanaway, University of Washington, United StatesElizaveta Elgaeva, Institute of Cytology and Genetics (RAS), Russia
Copyright © 2023 Li, Lu, Qian, Feng, Luo, Luo, He, Dong, Liu, Liu and Su. 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: Chengguo Su, suchengguo19840804@126.com; Zhongxing Liu, 15931235@qq.com; Huahui Liu, 442121644@qq.com