
==== Front
Medicine (Baltimore)
Medicine (Baltimore)
MD
Medicine
0025-7974
1536-5964
Lippincott Williams & Wilkins Hagerstown, MD

MD-D-24-06249
00004
10.1097/MD.0000000000039521
3
4400
Research Article
Observational Study
Causal relationship between systemic circulatory inflammatory regulators and intervertebral disc degeneration: A bidirectional 2-sample Mendelian randomization study
Chen Zi-Xuan MD 120863313@QQ.com
ab
Xu Bo PhD 619835857@qq.com
a
Huang Ze-Ling MD 1127062302@qq.com
ab
Liu Yu-Jiang MD 601895480@qq.com
a
Li Yu-Wei PhD lyw97538@126.com
a
Lu Bin-Jie MD lbj900212@sina.com
a
Lin Jun MD 490330850@QQ.com
ab
Zhang Xian-Da MD 2946146575@qq.com
ac
https://orcid.org/0009-0007-2891-8955
Shen Xiao-Feng PhD a*
a Department of Orthopedics, Suzhou TCM Hospital Affiliated to Nanjing University of Chinese Medicine, Suzhou, Jiangsu, China
b Nanjing University of Chinese Medicine, Nanjing, Jiangsu, China
c Anhui University of Chinese Medicine, Hefei, Anhui, China.
* Correspondence: Xiao-Feng Shen, Department of Orthopedics, Suzhou TCM Hospital Affiliated to Nanjing University of Chinese Medicine, Suzhou, Jiangsu 215003, China (e-mail: 29240818@qq.com).
06 9 2024
06 9 2024
103 36 e3952105 6 2024
12 7 2024
09 8 2024
Copyright © 2024 the Author(s). Published by Wolters Kluwer Health, Inc.
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the Creative Commons Attribution License 4.0 (CCBY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In the context of the development of intervertebral disc degeneration (IDD), inflammatory mediators play a pivotal role. Nevertheless, due to the influence of the inflammatory microenvironment, the causal relationship between specific inflammatory mediators and the development of IDD remains uncertain. The understanding of the causal relationship between inflammatory mediators and IDD is of great importance in preventing and delaying disc degeneration in the future. We utilized genetic data concerning systemic circulating inflammatory regulators obtained from a Genome-Wide Association Study (GWAS) analyzing 41 serum cytokines in a cohort of 8293 individuals from Finland. The genetic data for IDD were derived from the most recent GWAS summary statistics conducted within the FinnGen consortium, encompassing 37,636 IDD cases and 270,964 controls. Our analysis employed bidirectional 2-sample Mendelian randomization (MR) techniques, which included several MR methods such as MR Egger, weighted median, inverse variance weighted, weighted mode, and simple mode. Additionally, the MR-PRESSO method was employed to identify horizontal pleiotropy, heterogeneity was quantified using the Cochran Q statistic, and MR-Egger intercept analysis was performed to assess pleiotropy. We established causal relationships between 3 specific inflammatory factors and IDD. Elevated levels of MIP-1β (OR = 0.956, 95% CI: −0.08 to −0.006; P = .02) and IFN-G (OR = 0.915, 95% CI: −0.16 to −0.02; P = .01) expression were associated with a reduced risk of IDD. Conversely, genetic susceptibility to IDD was linked to a decrease in IL-13 levels (OR = 0.967, 95% CI: −0.063 to −0.004; P = .03). In this study, we have identified inflammatory factors that exhibit a causal relationship with the onset and progression of IDD, as supported by genetic predictions.

genetic relationships
infection susceptibility
inflammatory regulators
intervertebral disc degeneration
Mendelian randomization
National Natural Science Foundation of China 10.13039/501100001809 82174399 Xiao-Feng ShenNatural Science Foundation of Jiangsu Province 10.13039/501100004608 BK20211084 Xiao-Feng ShenCollaborative Innovation Center for Modern Science and Technology and Industrial Development of Jiangxi Traditional Medicine 10.13039/501100012248 ZD202230 Xiao-Feng ShenConstruction Project of Key Laboratory of Traditional Chinese Medicine Orthopedics and TraumatologyJSDW202253, SZS2022019 Xiao-Feng ShenOPEN-ACCESSTRUE
SDCT
==== Body
pmc1. Introduction

Lower back pain (LBP) stands as the most prevalent musculoskeletal disorder, ranking first in terms of disability.[1] Intervertebral disc degeneration (IDD) constitutes a significant cause of LBP, with approximately 40% of chronic LBP cases being associated with IDD.[2] In recent years, degenerative spinal conditions have triggered severe LBP among middle-aged and even younger individuals, thus evolving into a public health concern. This phenomenon greatly impacts the quality of life and imposes a substantial economic burden on society. The global prevalence of LBP stands at 9.4%, and the associated healthcare costs have surged by 10% over the past 43 years.[3] Nevertheless, the pathogenic mechanisms of IDD remain unclear, and a deeper understanding of its etiology is paramount for the prevention and treatment of LBP.

Inflammatory environments are believed to play a pivotal role in the progression of IDD,[4] yet their mechanisms are complex and not entirely elucidated. In specific circumstances, the inflammatory microenvironment plays a crucial role in maintaining a stable environment within the intervertebral disc (IVD). For instance, there exists a delicate balance between IL-1β and its inhibitor, IL-1Ra, within the normal IVD to sustain the homeostasis of the extracellular matrix (ECM).[5] Simultaneously, IL-1β and TNF-α, at particular concentrations, cooperatively uphold the equilibrium between synthesis and degradation of the extracellular matrix.[6] However, various factors can disrupt the biomechanical structure and cellular homeostasis within the normal IVD, stimulating the upregulation of multiple inflammatory factors, thus giving rise to an inflammatory microenvironment. Furthermore, the IVD’s own inflammatory microenvironment recruits additional exogenous inflammatory cells, further exacerbating the severity of IDD.[7] These pieces of evidence suggest that inflammation serves as both a cause of IDD and a consequence of a vicious cycle following IDD’s onset.

Nonetheless, current research faces various limitations: Replicating the complex inflammatory microenvironment within the IVD presents a significant challenge. It is difficult to eliminate the interplay between individual factors. Observing individual inflammatory factor levels at specific time points in clinical settings may not provide a comprehensive reflection of the overall inflammatory state during the disease process. Hence, there is an immediate need for an accurate assessment of the relationship between circulating inflammatory factors and the risk of IDD.

Mendelian randomization (MR) leverages genetic variations as instrumental variables (IVs) to estimate causal relationships between exposures, intermediates, and outcomes, thereby mitigating biases stemming from residual confounding and reverse causality, which can be encountered by conventional epidemiological approaches.[8] Numerous prior investigations have successfully established causal links between circulating inflammatory factors and various elements.[9–11] Consequently, in this 2-sample MR analysis, conducted bidirectionally, our primary aim is to explore the causal connection between systemic circulatory inflammatory regulators and the development of IDD. Revealing the relationship between the 2 at the genetic level has important implications for future targeted conservative treatment of degenerative discs or accelerated rehabilitation of patients after surgery.

2. Materials and methods

2.1. Study design

The overall workflow of this study is depicted in the figure. MR analysis requires meeting 3 assumptions: a strong association between IVs and exposure, IVs independence from confounders, and IVs relevance exclusively to the outcome through the exposure.[12] This study utilizes bidirectional 2-sample MR analysis to thoroughly examine the causal relationship between systemic inflammatory regulators and IDD. Our research design is summarized in Figure 1.

Figure 1. Research design.

2.2. Data sources

Both datasets employed in this MR analysis are extracted from summary data accessible to the public, originating from Genome-Wide Association Studies (GWAS). We obtained the most recent GWAS summary statistics for IDD from the FinnGen consortium R9 version (r9.finngen.fi), comprising 37,636 cases and 270,964 controls.[13] Data for 41 circulating inflammatory regulators were extracted from GWAS involving 8293 individuals across 3 distinct cohorts.[14]

2.3. Single nucleotide polymorphism filtering

The core assumptions of MR require that all single nucleotide polymorphisms (SNPs) act as strong and independent predictors of the exposure, demonstrating genome-wide statistical significance. In our study, which examined 41 inflammatory factors and IDD, we chose a relatively lenient yet highly significant P-value threshold of 5 × 10−6. This decision was made because a more stringent threshold of 5 × 10−8 would exclude a substantial number of SNPs, making the study unfeasible.[15]

To address linkage disequilibrium, we applied an R2 < 0.001 and a physical distance of 10,000 kb. Palindrome SNPs were removed from the analysis. Additionally, we considered genetic variants with an F-value (calculated as F = β²_exposure/SE²_exposure) below 10 as weak IVs and excluded them from the study. These methodological steps strengthened the association between IVs and exposure factors, thereby enhancing the scientific validity and credibility of our MR analyses.

2.4. Bidirectional MR analysis

We employed 5 MR methods to calculate the relationship between the impact of SNPs on the outcome and their influence on exposure, consolidating the results obtained for each SNP to thoroughly evaluate the causal connection between exposure and outcome. The 5 MR methods include Inverse Variance Weighted (IVW), MR-Egger regression, weighted median method, simple mode, weighted mode. Among these, the results obtained through the IVW method serve as the primary reference in this study,[16] with P < .05 considered statistically significant evidence. We detected horizontal pleiotropy through a global test using the MR-PRESSO method and corrected potential pleiotropic outliers by removing SNP outliers when necessary.[17] Finally, we quantified heterogeneity through the Cochran Q statistic and conducted MR-Egger intercept analysis to test for pleiotropy, with P > .05 serving as a reference value indicating no significant heterogeneity and pleiotropy.[18,19]

MR analysis, including 2-sample and multivariable MR, was performed using the TwoSampleMR package in R software (version 4.2.2).

3. Results

3.1. Causal effects of 41 circulating inflammatory regulators on IDD

Out of the initial 41 systemic circulatory inflammatory regulators, a total of 428 SNPs were selected. The primary prediction method, IVW MR, demonstrated that increased expression levels of 2 inflammatory regulators, MIP-1β (OR = 0.956, CI: −0.08 to −0.006; P = .02) and IFN-G (OR = 0.915, CI: −0.16 to −0.02; P = .01), were associated with a decreased risk of IDD. Although not exhibiting statistical significance, the remaining 4 MR analysis methods provided some level of support for the aforementioned results (see Table S1, Supplemental Digital Content, http://links.lww.com/MD/N529; Figure 2 to 6). The Cochran Q method did not detect heterogeneity (P > .05), and both MR-Egger and MR-PRESSO tests did not reveal horizontal pleiotropy (P > .05). The “leave-one-out” method involves systematically removing each SNP 1 at a time, calculating the combined effect of the remaining SNPs, and observing whether the results vary after each SNP is excluded. If significant changes in the results occur when a particular SNP is removed, it suggests that this specific SNP has a substantial impact on the overall findings. To ensure that the results are not unduly influenced by any single SNP, we conducted a “leave-one-out” sensitivity test. This analysis revealed that the overall results remained largely consistent after excluding each SNP, demonstrating the robustness and reliability of our findings.

Figure 2. Effect of different exposures on outcomes.

Figure 3. Forest plot of the association between systemic circulating inflammatory modulators and IDD. (A) IFN-G and IDD; (B) MIP-1β and IDD; (C) IDD and IL-13.

Figure 4. Funnel plots of the association between systemic circulating inflammatory regulators and IDD. (A) IFN-G versus IDD; (B) MIP-1β versus IDD; (C) IDD versus IL-13.

Figure 5. Leave-one-out sensitivity analysis of the relationship between systemic circulating inflammatory modulators and IDD. (A) IFN-G and IDD; (B) MIP-1β and IDD; (C) IDD and IL-13.

Figure 6. Scatter plots of the association between systemic circulating inflammatory modulators and IDD. (A) IFN-G versus IDD; (B) MIP-1β versus IDD; (C) IDD versus IL-13.

3.2. Causal effects of IDD on 41 circulating inflammatory regulators

To assess the reverse causal effects, we extracted 135 SNPs associated with IDD, with significance levels at P < 5 × 10−6. As some SNPs were not applicable to certain cytokines, varying numbers of SNPs were used for different cytokines, which can be found in Table S2, Supplemental Digital Content, http://links.lww.com/MD/N529. Using the IVW method, We observed that an increased genetic predisposition to IDD correlated with a reduction in the levels of 2 factors: GROA (OR = 0.961, 95% CI: −0.071 to −0.008; P = .02) and IL-13 (OR = 0.967, 95% CI: −0.063 to −0.004; P = .03) (Figs. 2-6). As the analysis using the simple mode in GROA yielded results contrary to the other 4 methods, we decided to exclude GROA. Sensitivity analysis was also performed, which did not detect heterogeneity or horizontal pleiotropy (P > .05), and the leave-one-out sensitivity test did not reveal variation in the overall outcome.

4. Discussion

This study employed a 2-sample MR analysis, utilizing data from the FinnGen project and GWAS data, to assess the causal association between 41 circulating inflammatory regulators and IDD. The study results indicated that higher levels of MIP-1β and IFN-γ were associated with a lower risk of IDD, while the reverse causal analysis revealed an association between genetic susceptibility to IDD and decreased levels of IL-13.

Inflammation consistently plays a significant role in the initiation and progression of IDD.[20] Various factors, both intrinsic and extrinsic to the intervertebral disc, can create an inflammatory microenvironment. These factors include microcrystals, degradation products of the ECM, mechanical loading, and immune cell infiltration. They stimulate disc cells to produce various inflammatory mediators, such as cytokines, chemokines, NLRP3 inflammasomes, prostaglandin E2 (PGE2), IFN-γ.[21] These inflammatory mediators contribute to the pathology of IDD in multiple ways. They promote the expression and activity of ECM-degrading enzymes, leading to structural and functional damage to the ECM. [22] Additionally, they induce senescence, apoptosis, and abnormal proliferation of disc cells[7]; Inflammatory mediators also increase neurovascular generation and the production of pain-related factors, which cause nerve fiber growth and innervation within the intervertebral disc.[23] Therefore, regulating the inflammatory microenvironment and restoring the balance of cytokines are effective strategies for controlling IDD.

Previous observational studies have indicated that the expression of inflammatory factors such as IL-1β and TNF-α significantly increases in IDD, and their expression levels are positively correlated with patient age and the severity of IDD.[24–26] Furthermore, a higher expression of the NLRP3 inflammasome in IDD can mediate the production of various inflammatory cytokines, further promoting IDD development.[25]

However, these observational studies might be subject to confounding factors or reverse causality that could distort the true causal relationship. For example, Weber et al[27] found that serum levels of IL-6 in IDD patients showed no significant association with the severity of IDD as indicated by MRI. These findings emphasize the complex relationship between inflammation and IDD.

MIP-1β, also known as CCL4, is a CC chemokine that plays a role in inflammation and host defense mechanisms by interacting with its specific receptor CCR5.[28] Previous studies have suggested that MIP-1β expression levels are significantly elevated in severely degenerated NP tissue compared to mildly degenerated NP tissue.[29,30] IFN-γ, also known as IFN-gamma, is 1 of the main inducers of M0 macrophage polarization.[31] Gorth et al[32] found that transgenic mice with overexpressed TNF-α exhibited early spontaneous IVD protrusion and elevated levels of IL-6, TNF-α, and IFN-γ upon NP exposure. Furthermore, research has shown that IFN-γ can increase the expression of TNF-α, IL-1β, COX-2, MMP-3, MMP-13, ADAMT-4, and ADAMT-5 in IVD cells.[33] Interleukin-13 (IL-13) is an anti-inflammatory factor produced by Th2 cells, macrophages, and mast cells, and it participates in tissue remodeling and fibrosis through activation of the JAK/STAT signaling pathway.[34] IL-13 receptors are expressed on macrophages and can inhibit M1 polarization while promoting M2 polarization, serving anti-inflammatory and tissue repair functions.[35] Currently, there is limited research on the relationship between IL-13 and IDD.

In contrast to many previous observational studies, our study suggests that higher levels of IFN-γ and MIP-1β may lead to a decreased genetic risk for IDD, which could be attributed to the immunoregulation and tissue repair abilities of inflammatory factors on degenerated intervertebral discs. Yu et al[36] proposed that the loss of immune privilege caused by intervertebral disc lesions could be a significant factor in the development of IDD and sciatica. While the infiltration of inflammatory cells, growth factors, and cytokines may promote IDD, the functions of specific immune cells may aid in the regeneration of damaged disc tissues and the restoration of immune privilege. It has been reported that IFN-γ-induced inflammation of mesenchymal stem cells may improve their ability to suppress T-cell proliferation, serving an anti-inflammatory and tissue repair function.[37] In the tumor microenvironment, IFN-γ and MIP-1β have shown dual roles in both antitumor and pro-tumor effects.[38,39] Whether these roles are relevant to the process of IDD remains an avenue for further research. In the assessment of reverse causality, we found that the progression of IDD is associated with decreased levels of IL-13. We speculate that during the IDD process, the destabilization of the intracellular environment leads to increased expression of pro-inflammatory factors such as TNF-α and IL-1β, while concurrently resulting in a decrease in the expression of anti-inflammatory factors such as IL-13, ultimately creating an inflammatory microenvironment within the IVD. Therefore, activating the expression of anti-inflammatory factors like IL-13 within the IVD may be a strategy to slow the progression of IDD.

Establishing a causal relationship between inflammatory factors and IDD holds significant potential for advancing targeted therapies and preventive strategies. Understanding these relationships enables the development of precision medicine approaches, including the use of biologic therapies such as monoclonal antibodies and gene therapy to modulate specific inflammatory mediators like IL-1β, TNF-α, and IFN-γ. These targeted therapies promise improved efficacy and reduced side effects compared to nonspecific anti-inflammatory treatments. Additionally, existing drugs that target these inflammatory pathways can be repurposed, accelerating the availability of new treatments. Early identification of genetic markers for heightened inflammatory responses facilitates risk stratification and early intervention, crucial for preventing significant degeneration. This allows for personalized preventive measures, including lifestyle modifications and early pharmacologic interventions tailored to high-risk individuals. Public health strategies can also be informed by these insights, guiding recommendations for physical activity, occupational health, and nutrition to reduce the incidence of IDD. Research and clinical applications benefit from an enhanced understanding of IDD pathophysiology. This understanding leads to better-designed clinical trials that utilize identified biomarkers as endpoints, improving their accuracy and relevance. Moreover, a holistic treatment approach that integrates anti-inflammatory strategies with physical therapy, surgical interventions, and other treatments can provide more comprehensive care plans for patients. Overall, recognizing the causal role of inflammation in IDD opens new avenues for innovative therapies and preventive measures. This understanding not only advances the development of targeted and effective medical interventions but also supports the creation of personalized and population-based strategies to manage and prevent IDD, ultimately improving patient outcomes and quality of life.

This study represents the first application of MR to examine the link between inflammation and IDD. This method helps overcome the problems of reverse causality bias that are common in traditional observational studies. Furthermore, we imposed strict study restrictions to reduce heterogeneity and pleiotropy in the outcome measures, providing reliable causal relationship estimates. Another notable aspect is the wide availability of serum samples in clinical studies, making our findings amenable to further validation and application. While the typical genome-wide significance P-value threshold is set at 5 × 10−8, our study used a relatively lenient threshold of 5 × 10−6 to explore more possibilities.

MR is a robust method for inferring causality, yet it is imperative to recognize its inherent limitations. MR relies on key assumptions: genetic variants must exhibit a strong association with the exposure (relevance), SNPs should be independent of confounders influencing both the exposure and the outcome (independence), and genetic variants should impact the outcome solely through the exposure, not via alternative pathways (exclusion restriction). Violations of these assumptions can result in biased estimates due to weak instruments, residual confounding, and pleiotropy, wherein a single genetic variant affects multiple traits. Despite genetic randomization, unmeasured confounding persists as a concern, as pleiotropy and population stratification – variations in allele frequencies across subpopulations – can introduce bias if not properly controlled. Furthermore, measurement error in the exposure can cause attenuation bias, underestimating the true effect of the exposure on the outcome, underscoring the necessity of accurate exposure measurement. The generalizability of MR study results is also limited, given the potential influence of gene-environment interactions or population-specific allele frequencies, thus necessitating caution when extrapolating findings to diverse demographic groups. Lastly, large sample sizes are often essential to detect small effects, particularly when genetic variants account for a minor proportion of the variance in the exposure; inadequate sample sizes may result in insufficient power and inconclusive outcomes.

5. Conclusion

In conclusion, this 2-sample MR study suggests a causal relationship between inflammation and IDD. Higher levels of MIP-1β and IFN-γ appear to reduce the risk of IDD, and the progression of IDD is associated with a decrease in IL-13 levels. These findings provide new insights into strategies for preventing and delaying IDD progression.

Author contributions

Funding acquisition: Xiao-Feng Shen.

Supervision: Xiao-Feng Shen.

Writing – original draft: Zixuan Chen.

Writing – review & editing: Bo Xu, Ze-Ling Huang, Yu-Jiang Liu, Yu-Wei Li, Bin-Jie Lu, Jun Lin, Xian-Da Zhang, Xiao-Feng Shen.

Supplementary Material

Abbreviations:

ADAMT-4 a disintegrin and metalloproteinase with thrombospondin motifs 4

ADAMT-5 a disintegrin and metalloproteinase with thrombospondin motifs 5

CCR5 C-C chemokine receptor type 5

CI confidence interval

COX-2 cyclooxygenase-2

ECM extracellular matrix

GWAS Genome-Wide Association Study

IDD intervertebral disc degeneration

IFN-G interferon gamma

IL-13 interleukin-13

IL-1Ra interleukin 1 receptor antagonist

IL-1β interleukin 1 beta

IVD intervertebral disc

IVs instrumental variables

IVW inverse variance weighted

LBP lower back pain

MIP-1β macrophage inflammatory protein 1 beta

MMP-13 matrix metalloproteinase-13

MMP-3 matrix metalloproteinase-3

MR Mendelian randomization

MR-PRESSO Mendelian Randomization Pleiotropy RESidual Sum and Outlier

MSCs mesenchymal stem cells

NLRP3 NOD-, LRR-, and pyrin domain-containing protein 3

OR odds ratio

PGE2 prostaglandin E2

SNP single nucleotide polymorphism

STAT signal transducer and activator of transcription

TNF-α tumor necrosis factor alpha

This study was supported by the National Natural Science Foundation of China (82174399), Natural Science Foundation of Jiangsu Province (BK20211084), Science and Technology Development Program Project of Jiangsu Provincial Administration of Traditional Chinese Medicine (ZD202230), Construction Project of Key Laboratory of Traditional Chinese Medicine Orthopedics and Traumatology (JSDW202253, SZS2022019).

Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.

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.

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Supplemental Digital Content is available for this article.

How to cite this article: Chen Z-X, Xu B, Huang Z-L, Liu Y-J, Li Y-W, Lu B-J, Lin J, Zhang X-D, Shen X-F. Causal relationship between systemic circulatory inflammatory regulators and intervertebral disc degeneration: A bidirectional 2-sample Mendelian randomization study. Medicine 2024;103:36(e39521).

This manuscript was previously posted to Research Square: doi: https://doi.org/10.21203/rs.3.rs-3619712/v2.
==== Refs
References

[1] GBD Disease I Injury C Prevalence. 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.30496104
[2] Li Y Zhang H Zhu D . Notochordal cells: a potential therapeutic option for intervertebral disc degeneration. Cell Prolif. 2023;57 :e13541.37697480
[3] Yi J Zhou Q Huang J Niu S Ji G Zheng T . Lipid metabolism disorder promotes the development of intervertebral disc degeneration. Biomed Pharmacother. 2023;166 :115401.37651799
[4] Wuertz K Haglund L . Inflammatory mediators in intervertebral disk degeneration and discogenic pain. Global Spine J. 2013;3 :175–84.24436868
[5] Le Maitre CL Freemont AJ Hoyland JA . The role of interleukin-1 in the pathogenesis of human intervertebral disc degeneration. Arthritis Res Ther. 2005;7 :R732–45.15987475
[6] 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.17688691
[7] Molinos M Almeida CR Caldeira J Cunha C Gonçalves RM Barbosa MA . Inflammation in intervertebral disc degeneration and regeneration. J R Soc Interface. 2015;12 :20150429.26040602
[8] Smith GD Ebrahim S . “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32 :1–22.12689998
[9] Qin J Zhang L Ke B Liu T Kong C Jin C . Causal relationships between circulating inflammatory factors and IgA vasculitis: a bidirectional Mendelian randomization study. Front Immunol. 2023;14 :1248325.37753071
[10] Yeung CHC Schooling CM . Systemic inflammatory regulators and risk of Alzheimer’s disease: a bidirectional Mendelian-randomization study. Int J Epidemiol. 2021;50 :829–40.33313759
[11] Li H Li M Liu C . Causal effects of systemic inflammatory regulators on chronic kidney diseases and renal function: a bidirectional Mendelian randomization study. Front Immunol. 2023;14 :1229636.37711613
[12] Burgess S Thompson SG . Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32 :377–89.28527048
[13] Kurki MI Karjalainen J Palta P . FinnGen. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. 2023;613 :508–18.36653562
[14] Ahola-Olli AV Würtz P Havulinna AS . Genome-wide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors. Am J Hum Genet. 2017;100 :40–50.27989323
[15] Wang Q Shi Q Lu J Wang Z Hou J . Causal relationships between inflammatory factors and multiple myeloma: a bidirectional Mendelian randomization study. Int J Cancer. 2022;151 :1750–9.35841389
[16] Lawlor DA Harbord RM Sterne JA Timpson N Davey Smith G . Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27 :1133–63.17886233
[17] 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 :693–8.29686387
[18] Bowden J Del Greco MF Minelli C . A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. 2017;36 :1783–802.28114746
[19] 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.26050253
[20] Khan AN Jacobsen HE Khan J . Inflammatory biomarkers of low back pain and disc degeneration: a review. Ann N Y Acad Sci. 2017;1410 :68–84.29265416
[21] Chen X Wang Z Deng R Yan H Liu X Kang R . Intervertebral disc degeneration and inflammatory microenvironment: expression, pathology, and therapeutic strategies. Inflamm Res. 2023;72 :1811–28.37665342
[22] Purmessur D Cornejo MC Cho SK Hecht AC Iatridis JC . Notochordal cell-derived therapeutic strategies for discogenic back pain. Global Spine J. 2013;3 :201–18.24436871
[23] Johnson ZI Schoepflin ZR Choi H Shapiro IM Risbud MV . Disc in flames: roles of TNF-α and IL-1β in intervertebral disc degeneration. Eur Cell Mater. 2015;30 :104–16; discussion 116.26388614
[24] Altun I . Cytokine profile in degenerated painful intervertebral disc: variability with respect to duration of symptoms and type of disease. Spine J. 2016;16 :857–61.26975459
[25] Chen ZH Jin SH Wang MY . Enhanced NLRP3, caspase-1, and IL- 1β levels in degenerate human intervertebral disc and their association with the grades of disc degeneration. Anat Rec (Hoboken). 2015;298 :720–6.25284686
[26] Weiler C Nerlich AG Bachmeier BE Boos N . Expression and distribution of tumor necrosis factor alpha in human lumbar intervertebral discs: a study in surgical specimen and autopsy controls. Spine (Phila Pa 1976). 2005;30 :44–53; discussion 54.15626980
[27] Weber KT Alipui DO Sison CP . Serum levels of the proinflammatory cytokine interleukin-6 vary based on diagnoses in individuals with lumbar intervertebral disc diseases. Arthritis Res Ther. 2016;18 :3.26743937
[28] Bondue A Jao SC Blanpain C Parmentier M LiWang PJ . Characterization of the role of the N-loop of MIP-1 beta in CCR5 binding. Biochemistry. 2002;41 :13548–55.12427015
[29] Li Z Wang X Pan H . Resistin promotes CCL4 expression through toll-like receptor-4 and activation of the p38-MAPK and NF-κB signaling pathways: implications for intervertebral disc degeneration. Osteoarthritis Cartilage. 2017;25 :341–50.27737814
[30] Zhang J Li Z Chen F . TGF-β1 suppresses CCL3/4 expression through the ERK signaling pathway and inhibits intervertebral disc degeneration and inflammation-related pain in a rat model. Exp Mol Med. 2017;49 :e379.28935976
[31] Zhu L Zhao Q Yang T Ding W Zhao Y . Cellular metabolism and macrophage functional polarization. Int Rev Immunol. 2015;34 :82–100.25340307
[32] Gorth DJ Shapiro IM Risbud MV . Transgenic mice overexpressing human TNF-α experience early onset spontaneous intervertebral disc herniation in the absence of overt degeneration. Cell Death Dis. 2018;10 :7.30584238
[33] Yang H Liu B Liu Y . Secreted factors from intervertebral disc cells and infiltrating macrophages promote degenerated intervertebral disc catabolism. Spine (Phila Pa 1976). 2019;44 :E520–9.30540714
[34] Junttila IS . Tuning the cytokine responses: an update on interleukin (IL)-4 and IL-13 receptor complexes. Front Immunol. 2018;9 :888.29930549
[35] Hershey GK . IL-13 receptors and signaling pathways: an evolving web. J Allergy Clin Immunol. 2003;111 :677–90; quiz 691.12704343
[36] Yu P Mao F Chen J . Characteristics and mechanisms of resorption in lumbar disc herniation. Arthritis Res Ther. 2022;24 :205.35999644
[37] Cassano JM Schnabel LV Goodale MB Fortier LA . Inflammatory licensed equine MSCs are chondroprotective and exhibit enhanced immunomodulation in an inflammatory environment. Stem Cell Res Ther. 2018;9 :82.29615127
[38] Gocher AM Workman CJ Vignali DAA . Interferon-γ: teammate or opponent in the tumour microenvironment? Nat Rev Immunol. 2022;22 :158–72.34155388
[39] Mukaida N Sasaki SI Baba T . CCL4 signaling in the tumor microenvironment. Adv Exp Med Biol. 2020;1231 :23–32.32060843
