Skip to main content

Genetic variants of interleukin 1B and 6 are associated with clinical outcome of surgically treated lumbar degenerative disc disease



Successfully surgically treating degenerative disc diseases can be challenging to the spine surgeons, the long-term outcome relies on both the physical and mental status of the patient before and after treatment. Numerous studies underlined the role of inflammatory cytokines – like interleukin 1B and 6 – in the development of chronic diseases such as failed back surgery syndrome (FBSS) and major depressive disorder (MDD) which alter the outcome after spinal surgery. Our aim was to evaluate the associations of IL6 and IL1B gene polymorphisms with the long-term outcome of degenerative lumbar spine surgeries.


An international genetical database (GENODISC) was combined with our institute’s clinical database to create a large pool with long term follow up data. Altogether 431 patient’s data were analysed. Patient reported outcome measures and surgical outcome was investigated in association with IL1B and IL6 SNPs with the help of ‘SNPassoc’ R genome wide association package.


Interleukin 1B variants analysis confirmed association with improvement of pain after surgery on individual SNP level and on haplotype level, moreover relationship with patient reported outcome and preoperative level of depression was found on individual SNP level. IL6 variants were associated with preoperative depression, somatization and with subsequent surgery.


Understanding the complexity of spinal surgery patients’ long-term well-being is crucial in effectively treating chronic debilitating somatic diseases and the associated mental illnesses. Further studies should investigate more comprehensively the linkage of chronic physical and mental illnesses focusing on their simultaneous treatment.

Peer Review reports


Degenerative disc disease (DDD) is a chronic and debilitating condition, which leads to loss of workdays in an active adult’s life [1]. Conservative treatments are mainly effective but occasionally surgical treatment is inevitable [2]. The surgical intervention in DDD aims to reduce pain and restore function. Surgical outcome is complex and multifactorial. It can be measured with objective (e.g.: muscle strength, vegetative functions) and subjective (patient reported outcome measures ‘PROMs’) assessment tools. Not uncommonly, despite the perfect surgical technique for an obvious pathology, the patient reports no significant improvement and continues to suffer from pain and even failed back surgery syndrome (FBSS) can develop [3].

Pain processing is a complex multifactorial pathway with different regulatory molecules. It has a biomedical component defined by tissue damage, an evaluative component which is influenced by coping mechanisms and an affective component which is altered by psychological disorders like depression and anxiety [4]. Cytokines and interleukins – such as interleukin 6 (IL6) and interleukin 1B (IL1B) - have an important role in regulating local pain response [5,6,7,8]. IL1B also act as an upregulator of other nociceptive agents and cascades such as IL6, prostaglandins, substance-P and matrix metalloprotease (MMP) 9 [9].

Single nucleotide polymorphisms (SNPs), the most common genetic variants have been identified in the genetic background of the development and outcome of several multifactorial diseases [10,11,12,13,14]. DDD is also a multifactorial entity with a strong genetic background [15] and high number of genes and their SNPs are linked to its pathomechanism [16,17,18]. Also genetic variations occurring in COL1A1, COL9a3 and VDR genes seems to be associated with the development of LDD [19]. DDD related pain is also influenced by different genes and their variations such as catechol-O-methyltransferase (COMT) [20] and β2-adrenergic receptor genes (ARDRB2) [5, 21]. According to published data the IL1 gene family (IL1A, IL1B, IL1RN) and the IL6 gene variations have connection with degenerative spinal pathologies [22,23,24], DDD related pain (low back pain, leg pain) [12, 25, 26] and the outcome of conservative treatment [27], but there is limited amount of data available about the associations of these gene variants in relation to spinal surgery outcome [28].

The current study focuses on evaluating the relationship between IL1B and IL6 gene polymorphisms and the long-term outcome of degenerative lumbar spine surgeries.

Materials and methods

Study population

Data were collected prospectively from adults (above the age of 18) who underwent routine, elective surgery for lumbar disc degeneration at one or two levels at a tertiary spine center. Prospective clinical data were linked with the subjects’ genetic data derived from the GENODISC multicenter international collaboration. Patients with minimum 2-year follow-up data were included into the final study cohort to explore the long-term outcome of the surgical procedures. Patients reoperated within 2 years due to a surgical site infection, proximal junctional kyphosis (PJK) or adjacent segment degeneration (ASD) as well as subjects undergoing either acute intervention because of neurological emergency or tumour surgery were excluded from the study. Surgeries were performed by board-certified orthopaedic surgeons or neurosurgeons specified in spinal surgery. Applied procedures included microdiscectomy, decompression and instrumented fusion (transforaminal lumbar interbody fusion or posterior fusion). All procedures were carried out using the standard median-sagittal posterior approach. All subjects signed a written consent form describing the scientific purpose of the systematic collection of their clinical and genetic data. The study was approved by the Scientific and Research Ethics Committee of the Medical Research Council Hungary (431/PI/2007).

Clinical data

Patients completed standard and validated PROMs to assess their clinical status before the surgery and during the follow-up period [29, 30]. Pain was evaluated by the 10 cm long Visual Analogue Scale. Lumbar spine related function was measured with Oswestry Disability Index (ODI). Psychologic distress was measured by evaluating the level of depression and somatisation and was assessed with the Hungarian versions of Zung Depression Scale (ZDS) [30] and the Modified Somatic Perception Questionnaire (MSPQ) (Supplementary Material), respectively. Patients were asked to rate the overall outcome of the surgery using a five-category question; “helped a lot”, “helped”, ‘helped only little”, “didn’t help”, “made things worse”. To measure global treatment outcome (GTO) a dichotomous variable was generated based on these given answers. Good outcome was defined if the patient responded by ‘helped a lot, ‘helped’ and poor in case the patient replied by ‘only little’, ‘didn’t help’, ‘made things worse’ [31, 32]. Surgical outcome was considered “good” if no re-operation was performed at the index level within 2 years and “poor” if a subsequent surgery was needed within 2 years.


DNA was extracted from venous blood or saliva samples using commercial. Five SNPs in IL1B and four SNPs in IL6 genes were selected for genotyping based on previous literature data [11, 33,34,35,36,37]. Genotyping was performed from 2007 to 2013 at the Technology Centre, Institute for Molecular Medicine Finland (FIMM), University of Helsinki using a Sequenom MassArray technology and the iPLEX Gold reagents (Sequenom Inc., San Diego, USA).


Allelic and genotype distributions, Hardy-Weinberg equilibrium, minor allele frequency (MAF) as well as associations between genetic variants and outcomes were determined and analysed using the 'SNPassoc' and 'haplo.stats' R software packages [38]. Genetic associations with preoperative and postoperative pain, disability, and psychological distress as well as global treatment and surgical outcome were investigated. Individual genotype-phenotype associations were studied in generalized linear models (GLM). Genetic subgroups with less than 4 (1%) subjects were excluded from subsequent statistical analyses. Haplotype-phenotype association was analysed applying haplo.score tests and GLM models. In haplo.score analysis, a global test of association as well as individual haplotype-specific tests are carried out using a score function. Disc herniation subgroup (patients underwent microdiscectomy) was also analysed separately to investigate the role of IL SNPs in sciatica. Significant covariates (age, gender, weight, height, preop ZDS and preop MSPQ score, type of surgery) were determined and entered into the models for each outcome. P-value less than 0.05 was considered significant.


Study population

A total of 431 subjects (all Caucasians) met the study inclusion criteria. Mean age were 52.7 (SD:13.9y) years (from 20 to 88 years) and male/female ratio was 0.6 (male:166, female:265). As the index surgery 171 patients had microdiscectomy, 22 patients had decompression, 142 patients had one level fusion and 96 patients had 2-level fusion. In the final study cohort, 44 patients required a subsequent lumbar surgery at the index level during the follow-up. Eight patients had re-discectomy or decompression, 35 required fusion and in 1 case the implants had to be removed.

Descriptive statistics of genotyping

Table 1 shows the results of the genotyping process. The genotyping success rate was more than 97% in all cases. All studied SNPs were in Hardy-Weinberg equilibrium. Two haploblock from IL1B gene were identified consisting of 2-2 SNPs (‘rs1143634-rs1143633’and ‘rs1143627-rs16944’) and no haploblock was identified on the IL6 gene as seen on Fig. 1.

Table 1 The descriptive statistics of the genotyped SNPs
Fig. 1
figure 1

Linkage disequilibrium (LD) map of IL1B (A) and IL6 (B) SNPs. Squares are colored darker if the |D’| value is high, that is, LD is strong. Empty dark squares mean |D’|=1, that is, complete LD between two single nucleotide polymorphisms

Associations of IL1B and IL6 gene variants with preoperative PROMs

In the overall population the mean±SD values of preoperative ODI score was 47.4±18.4 and the mean VAS score was 7.2±1.9, the mean ZDS was 39.6±8.1, MSPQ was 8.3±5.7. No individual SNP was associated with preoperative ODI and pain (Table 2), however both IL genes had SNPs related to the level of depression. ‘T’ allele of rs1143627 IL1B SNP was associated with higher level of depression (ZDS was 40.6±8.7, 39.2±7.3 and 38.3±8.0 in case of ‘T/T’, ‘T/C’ and ‘C/C’ genotypes, respectively, p-value=0.025 in log additive model). IL1B rs16944 ‘G’ allele carriers also showed higher level of depression (ZDS was 40.6±8.8, 39.2±7.3 and 38.0±8.0 in case of ‘G/G’, ‘A/G’ and ‘A/A’ genotypes, respectively, p-value=0.025 in log additive model). rs1143634 IL1B was associated with ZDS in an overdominant model (p=0.025, “C/T” mean ZDS±SD was 40.8±8.4 and 39.0±7.8 in case of ‘C/T’ and ‘C/C’+’T/T’ genotype groups). The ‘C’ allele of IL6 SNP rs2069835 was linked to increased level of depression (mean ZDS±SD were 39.2±7.8, 42.2±9.2, and 45.3±10.1 in case of ‘T/T’, ‘T/C’ and ‘C/C’ genotypes, respectively, p=0.003 in log-additive model).

IL6 rs2069835 was associated with the level of preoperative somatization (Mean MSPQ±SD was 8.0±5.3 and 10.2±7.0 in case of ‘T/T’ and ‘T/C’+’C/C’ genotypes respectively, p=0.010, in dominant model) (Table 2).

IL1B haplotypes were not associated with preoperative ODI, depression, somatization, and pain (data not shown).

Table 2 Associations of IL1B and IL6 gene variants with preoperative PROMs

Associations of IL1B and IL6 gene variants with postoperative outcome

Change in pain and disability

The mean overall improvement in pain intensity was 3.4±3.2 points (overall 61% improvement) in the study cohort. IL1B rs1143633 was strongly associated with the change in the reported pain at follow-up, where the ‘A’ allele carriers had the largest improvement in pain intensity (mean change±SD (%) in pain intensity was -3.7±3.3 (50%) in ‘A/G’+’A/A’ group vs -2.9 ±3.2 (40%) in ‘G/G’ genotype, p=0.00085 in dominant model) (Table 3). Another IL1B SNP (rs1143634) was associated with sciatica in the disc herniation subgroup. In this cohort, the level of preoperative pain was significantly higher in the ‘CC’ genotype (VAS=7.5±1.9, 6.6±2.2 and 6.7±2.3 for ‘C/C’, ‘C/T’ and ‘T/T’ genotypes respectively, p=0.006 in dominant model) (Fig. 2). Change in ODI score was not associated with the studied gene variants.

Table 3 Associations of IL1B and IL6 gene variants with postoperative outcome
Fig. 2
figure 2

Preoperative pain in microdiscectomy subgroup (mean ±SD)

Global treatment outcome

In the study cohort 350 patients (82%) reported good outcome while 75 patients (17%) reported poor outcome (6 patients’ data were missing). The ‘C’ allele of IL1B rs1143627 was related with better GTO (OR:1.49, p=0.049 in log-additive model) (Table 3).

Surgical outcome

In the overall population 44 patients had poor surgical outcome (10.2%). All 4 IL6 SNPs were associated with the risk of reoperation within 2 years, even after adjusting to type of index surgery. ‘G/G’ genotype of rs1800796 (OR:6.6, p=0.009, dominant model), ‘G/A’ genotype of rs2069852 (OR:5, p=0.039 in codominant model) and ‘C’ allele of rs2069835 (p=0.027, OR:1.27 in log-additive model) were associated with worse outcome. rs2069861 was associated with surgical outcome in an overdominant model (p=0.014) (Table 3).

Results of haplotype analysis

There was one haploblock in IL1B gene (rs1143634-rs1143633) which was associated with change in pain. ‘C-A’ haplotype was associated with the greater improvement in pain compared to the most common ‘C-G’ haplotype (p= 0.001) (Table 4).

Table 4 IL1B haplotype association (GLM and hapscore) with the change in pain after surgery


Number of spine surgeries because of DDD is continuously increasing. The rate of patients with poor outcome is between 5-70% in different surgical cohorts, and chronic pain condition because of failed back surgery syndrome (FBSS) is also not uncommon in this population. Understanding the pathophysiology of chronic pain conditions - such as FBSS - can lead clinicians to develop and apply new therapeutic methods in order to alleviate pain and improve the quality of life in this large patient group. The well-being of a patient is determined by multiple musculoskeletal, functional, and psychosocial factors [32]. Genetic influence on surgical outcome has been also highlighted by previous studies [28]. In the present study, polymorphisms of two interleukin (IL1B, IL6) genes in a large cohort of 431 patients who underwent elective lumbar spinal surgery for DDD were investigated in terms of the therapeutical outcome. Relationship between long-term treatment results, psychological factors, pain, and different IL gene variants were supported by individual SNP associations and haplotype analyses. Outcome of routine lumbar degenerative surgeries was analysed in different dimensions. Associations of IL gene variants with change in pain, disability as well as patient-reported global treatment outcome and need for a subsequent surgery were determined to elucidate the potential genetic influence.

IL1B variants were significantly related to the improvement in pain after the spine surgery, ‘A’ allele of rs1143633 as well as ‘C-A’ haplotype of rs1143634-rs1143633 haploblock were associated with greater improvement in pain. No other gene variant was associated with pain relieve however when we analysed the microdiscectomy subgroup we found that patients with ‘C/C’ genotype of rs1143634 had significantly higher preoperative pain compared to the other genotypes. Other IL1B variant (rs1143627) was associated with patient reported global treatment outcome, while majority of the studied IL1B variants were related to the preoperative level of depression. Interestingly IL6 variants were significantly associated with the need for a subsequent surgery during the follow-up period. The ‘C’ allele of rs2069835 IL6 SNP was associated with a higher risk for reoperation and also with increased level of preoperative depression and somatization. None of the studied gene variants were associated with preoperative spinal pain and disability level.

Number of studies supported the relationship between intervertebral disc degeneration and IL1, IL6 gene variants [23, 24, 39,40,41,42]. SNPs of these genes have been showed to be associated with the outcome of different surgical treatment [43,44,45,46], but only Moen et al. have studied the possible association of IL1 gene family and long-term outcome in patients treated because of lumbar disc herniation so far [28]. They did not find a significant relationship between rs1143627 IL1B SNP and treatment outcome, however they did not publish the genetic effect of single SNPs but their combinations on a mixed (surgically and non-surgically treated) patient groups. The same SNP (rs1143627) was found to be associated with symptomatic disc herniation [26] and with DDD associated pain [25] by others. IL1B variants have been also described in association with DDD [22, 40]. IL6 variants have not been studied related to the surgical outcome of DDD yet but they were previously associated with the process of lumbar disc degeneration [24, 39, 42].

The association between IL1B, IL6 genetic variants and the therapeutic outcome after lumbar spinal surgeries can be explained by different mechanisms:

1) Progressive degeneration process can lead to persistent spinal pain and a potential indication of a subsequent surgery. IL1B is involved in multiple pathological process of disc degeneration. It stimulates extracellular matrix degradation, accelerates cellular senescence and induces apoptosis [47]. rs1143633 in IL1B was associated with improvement of pain after surgery in our study while this SNP was found to be associated previously with the higher occurrence of disc degeneration (HIZ) [40] what can be a potential chronic pain source. IL6 variants have been also described in relation to DDD [48, 12, 42].

2) Tissue damage is often mediated through local inflammation. Inflammatory mediators such as IL6 and IL1B carry an important role in regulating and sustaining inflammation and pain. Different studies showed their potential role in disc degeneration related inflammatory process [49,50,51]. IL6 is crucial in homeostasis maintenance and host defence but its overproduction can cause the development or progression of diseases (such as pathologic pain) [52, 53]. The serum level of IL6 is increased in herniated disc which promotes upregulation of MMPs [53, 54]. Kraychete et al. also showed that patients with chronic low back pain due to disc herniation had higher level of serum IL6 [55]. The tissue level of IL6 can be related to the genetic variant of the gene. For example, rs1800796 IL6 SNP (what we found to be strongly associated with FBSS) is associated with increased promoter activity boosting the local secretion of IL6 [12, 42]. The two genes have a potential influence on each other, while IL1B is described as one of the key local inducers of IL6 production [50, 51], [56]. Not surprisingly, the variants of IL1B and IL6 genes have been associated with other chronic inflammatory conditions such as periodontitis, cancer, osteoporosis, type 2 diabetes and diabetic nephropathy [33,34,35,36, 57,58,59,60,61,62].

3) Psychological issues are also important in pain response and in the development of chronic pain. Depression and anxiety have been previously described as risk factors of DDD and poor surgical outcome after spine surgeries [4, 32, 63, 64]. Interleukin genes can significantly influence the patient’s psychological profile. Chronic inflammation and dysregulation of the immune response is a key factor in the development of major depressive disorders (MDD) [65, 66]. Patients with MDD show an abnormal profile of pro- and anti-inflammatory circulating cytokines [66,67,68,69]. In animal models of MDD, increased level of pro-inflammatory cytokines caused central serotonin depletion, hypothalamic–pituitary–adrenal (HPA) axis dysregulation, microglial activation and brain structure alteration [66]. In animal inflammatory MDD (MDD-I) models, IL1B appears to be the initial triggering complex of the inflammatory cascade both centrally and peripherally [66]. In our study, some IL1B variants were significantly associated with the preoperative level of depression. These findings are in accordance with previous report about the positive association between rs1143627 IL1B polymorphisms and MDD [70]. In accordance with our findings, Yu et al. found that the homozygotic ‘T/T’ patients of rs1143634 had a tendency of suffering from less severe depressive symptoms than ‘T/C’ homozygotes [71]. ‘T/T’ genotype of rs1143627 is reported to have a strong connection with major recurrent depression [72], in the meantime we found that patients with this particular genotype had worse scores on the depression scale. Two of the investigated SNPs, rs16944 and rs1143627 are located in the promoter region of the gene. These polymorphisms lead to altered expression of IL1B which results in local inflammation and promotes the production of MMPs [73]. A study suggested that IL1B rs16944 gene polymorphism hinder the pharmacological response in the treatment of MDD by increasing the risk of non-remission over 6 weeks of antidepressant treatment [74]. Another IL1B SNP (rs1143633) was strongly associated with postoperative pain in our study while rs1143634 was strongly associated with the preoperative pain intensity but only in the disc herniation subgroup. Previously, association of intensity of back pain and this SNP have been published in war veterans with posttraumatic stress disorder [75]. Association between rs1143633 and pain have not been published yet, however there are a few studies investigating its relationship in paediatric MDD and schizophrenia [76, 77]. rs2069861 in IL6 was associated with both depression and somatization in our cohort. Somatization is also an important factor in the development of symptomatic DDD [64]. Genetic variants of IL6 were linked to depression, somatization and anxiety in numerous studies [78,79,80,81,82].

Recently published data showed the possible role of interleukin agonist drugs in the treatment of pathological pain (e.g., chronic pain, inflammatory pain etc.) [83], therefore novel therapeutic strategies targeting IL6 or its receptors have been developed and successfully used in the treatment of selected diseases. In a paper a single intradiscal injection of tocilizumab (IL6 receptor antibody) provided short-term alleviation of discogenic pain [84]. Variants of the interleukins’ and their receptors’ genes can modify the effect of this targeted anti-inflammatory therapies, however there is no data about that so far.

In genetic association studies the sample size is highly important, as it can significantly alter the results. However, the sample size varies in the published studies, hence inconsistent genome wide association study results with non-reproducible results exist [16]. Thus, in our study we aimed to avoid sample size related study bias by using a prospective international large dataset to strengthen the findings. There are some limitations of the present study. Selection and regional population bias cannot be ruled out fully because only Caucasian patients were enrolled to the study. We did not apply any correction of the alpha-level during the genetic association testing process. We followed this method because we used a hypothesis-driven approach where effect of candidate SNPs on a phenotype was calculated. Moreover, genetic associations were tested with different statistical models (individual SNP association, haplotype analysis) to confirm the associations of the study even if the Type I error rate was not reduced. Comorbidities (e.g.: psychiatric disorders) can influence the genetic associations even if we have adjusted the statistical analyses for individual level of depression and anxiety. Therefore, study population selection bias cannot be ruled out completely. These limitations above can influence the credibility of our findings, therefore independent replications of the study are strongly recommended on different populations.

In conclusion we can state that IL1B and IL6 gene variants are associated with the psychological status and the long-term outcome of surgically treated lumbar DDD patients and these associations can be related to each other. The most plausible explanation to these findings could be linked to the major role of these cytokines in local and systemic chronic inflammation. Based on our findings and the corresponding literature advanced treatment methods could be established targeting interleukin 1B, interleukin 6 and its genes to successfully prevent/treat FBSS or even primary lumbar degenerative pathologies. On the other hand, the consideration of patient-specific genetic difference can be important to maximize the therapeutic outcome.

Availability of data and materials

The datasets generated during the current study are not publicly available due to there is still ongoing research on the study data, but raw data are available from the corresponding author upon reasonable request.


  1. Nissen M, Ikäheimo TM, Huttunen J, Leinonen V, Von Und Zu Fraunberg M. Long-term outcome of spinal cord stimulation in failed back surgery syndrome: 20 years of Experiencewith 224 consecutive patients. Neurosurgery. 2019;84(5):1011–8.

    Article  PubMed  Google Scholar 

  2. Ma XL, Zhao XW, Ma JX, Li F, Wang Y, Lu B. Effectiveness of surgery versus conservative treatment for lumbar spinal stenosis: a system review and meta-analysis of randomized controlled trials. Int J Surg. 2017;44:329–38.

    Article  PubMed  Google Scholar 

  3. Thomson S. Failed back surgery syndrome – definition, epidemiology and demographics. Br J Pain. 2013;7(1):56–9.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Strøm J, et al. Anxiety and depression in spine surgery—a systematic integrative review. Spine J. 2018;18(7):1272–85.

    Article  PubMed  Google Scholar 

  5. Crofford LJ. Chronic pain: where the body meets the brain. Trans Am Clin Climatol Assoc. 2015;126:167–83.

    PubMed  PubMed Central  Google Scholar 

  6. Bhattacharya A, Jones DNC. Emerging role of the P2X7-NLRP3-IL1β pathway in mood disorders. Psychoneuroendocrinology. 2018;98(April):95–100.

    Article  CAS  PubMed  Google Scholar 

  7. Huang SJ, Yan JQ, Luo H, Zhou LY, Luo JG. IL-33/ST2 signaling contributes to radicular pain by modulating MAPK and NF-ΚB activation and inflammatory mediator expression in the spinal cord in rat models of noncompressive lumber disk herniation. J Neuroinflammation. 2018;15(1):1–12.

    Article  CAS  Google Scholar 

  8. Milenkovic VM, Stanton EH, Nothdurfter C, Rupprecht R, Wetzel CH. The role of chemokines in the pathophysiology of major depressive disorder. Int J Mol Sci. 2019;20(9):1–17.

    Article  CAS  Google Scholar 

  9. Ren K, Torres R. Role of interleukin-1β during pain and inflammation. Brain Res Rev. 2009;60(1):57–64.

    Article  CAS  PubMed  Google Scholar 

  10. Dantzer R, O’Connor JC, Freund GG, Johnson RW, Kelley KW. From inflammation to sickness and depression: when the immune system subjugates the brain. Nat Rev Neurosci. 2008;9(1):46–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Miller AH, Maletic V, Raison CL. Inflammation and its discontents: the role of cytokines in the pathophysiology of major depression. Psiquiatr Biol. 2010;17(2):71–80.

    Article  Google Scholar 

  12. Noponen-Hietala N, et al. Genetic variations in IL6 associate with intervertebral disc disease characterized by sciatica. Pain. 2005;114(1–2):186–94.

    Article  CAS  PubMed  Google Scholar 

  13. Ge D, et al. Genetic variation in IL28B predicts hepatitis C treatment-induced viral clearance. Nature. 2009;461(7262):399–401.

    Article  CAS  PubMed  Google Scholar 

  14. Lindgren M, et al. Genetic variation in IL28B (IFNL3) and response to interferon-alpha treatment in myeloproliferative neoplasms. Eur J Haematol. 2018;100(5):419–25.

    Article  CAS  PubMed  Google Scholar 

  15. Crites Battie M, Videman T, Gibbons LE, Fisher LD, Manninen H, Gill K. 1995 volvo award in clinical sciences determinants of lumbar disc degeneration -a study relating lifetime exposures and magnetic resonance imaging findings in identical twins. Spine (Phila Pa 1976). 1995;20(24):2601–12.

    Article  Google Scholar 

  16. Biczo A, Szita J, McCall I, Varga PP, Lazary A. Association of vitamin D receptor gene polymorphisms with disc degeneration. Eur. Spine J. 2020;29(3).

  17. Ala-Kokko L. Genetic risk factors for lumbar disc disease. Annals of Medicine. 2002;34(1):42–7.

    Article  CAS  PubMed  Google Scholar 

  18. Kalichman L, Hunter DJ. The genetics of intervertebral disc degeneration. Familial predisposition and heritability estimation. Joint Bone Spine. 2008;75(4):383–7.

    Article  PubMed  Google Scholar 

  19. Toktaş ZO, et al. Association of collagen I, IX and vitamin D receptor gene polymorphisms with radiological severity of intervertebral disc degeneration in southern European ancestor. Eur Spine J. 2015;24(11):2432–41.

    Article  PubMed  Google Scholar 

  20. Baumbauer KM, et al. Contribution of COMT and BDNF genotype and expression to the risk of transition from acute to chronic low Back pain. Clin J Pain. 2020;36(6):430–9.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Skouen JS, et al. Genetic variation in the beta-2 adrenergic receptor is associated with chronic musculoskeletal complaints in adolescents. Eur J Pain (United Kingdom). 2012;16(9):1232–42.

    Article  CAS  Google Scholar 

  22. Solovieva S, et al. Intervertebral disc degeneration in relation to the COL9A3 and the IL-1β gene polymorphisms. Eur Spine J. 2006;15(5):613–9.

    Article  PubMed  Google Scholar 

  23. Virtanen IM, et al. Occupational and genetic risk factors associated with intervertebral disc disease. Spine (Phila Pa 1976). 2007;32(10):1129–34.

    Article  Google Scholar 

  24. Kelempisioti A, et al. Genetic susceptibility of intervertebral disc degeneration among young Finnish adults. BMC Med Genet. 2011;12(1):153.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Solovieva S, Leino-Arjas P, Saarela J, Luoma K, Raininko R, Riihimäki H. Possible association of interleukin 1 gene locus polymorphisms with low back pain. Pain. 2004;109(1–2):8–19.

    Article  CAS  PubMed  Google Scholar 

  26. Paz Aparicio J, et al. The IL-1β (+3953 T/C) gene polymorphism associates to symptomatic lumbar disc herniation. Eur Spine J. 2011;20(Suppl 3):383–9.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Kim DH, Lee SH, Kim KT, Yu SD. Association of interleukin-1 receptor antagonist gene polymorphism with response to conservative treatment of lumbar herniated nucleus pulposus. Spine (Phila Pa 1976). 2010;35(16):1527–31.

    Article  Google Scholar 

  28. Moen A, Schistad EI, Rygh LJ, Re C, Gjerstad J. Role of IL1A rs1800587, IL1B rs1143627 and IL1RN rs2234677 genotype regarding development of chronic lumbar radicular pain; a prospective one-year study. PLoS One. 2014;9(9):1–5.

    Article  CAS  Google Scholar 

  29. Valasek T, Varga PP, Szövérfi Z, Kümin M, Fairbank J, Lazary A. Reliability and validity study on the Hungarian versions of the oswestry disability index and the Quebec back pain disability scale. Eur Spine J. 2013;22(5):1010–8.

    Article  PubMed  PubMed Central  Google Scholar 

  30. S. A., “A Zung-féle Önértékelő Depresszió Skála (ÖDS),” Pszichodiagnosztikai Vadem. 1994. pp. 180–4.

  31. Klemencsics I, et al. Cross-cultural adaptation and validation of the Hungarian version of the Core outcome measures index for the Back (COMI Back). Eur Spine J. 2016;25(1):257–64.

    Article  PubMed  Google Scholar 

  32. Mannion AF, et al. Outcome assessment in low back pain: how low can you go? Eur Spine J. 2005;14(10):1014–26.

    Article  PubMed  Google Scholar 

  33. da Silva FRP, Vasconcelos ACCG, de Carvalho França LF, Di Lenardo D, Nascimento HMS, Vasconcelos DFP. Association between the rs1143634 polymorphism in interleukin-1B and chronic periodontitis: results from a meta-analysis composed by 54 case/control studies. Gene. 2018;668(April):97–106.

    Article  CAS  PubMed  Google Scholar 

  34. Hwang IR, et al. Effect of interleukin 1 polymorphisms on gastric mucosal interleukin 1β production in helicobacter pylori infection. Gastroenterology. 2002;123(6):1793–803.

    Article  CAS  PubMed  Google Scholar 

  35. Iacoviello L, et al. Polymorphisms of the interleukin-1β gene affect the risk of myocardial infarction and ischemic stroke at young age and the response of mononuclear cells to stimulation in vitro. Arterioscler Thromb Vasc Biol. 2005;25(1):222–7.

    Article  CAS  PubMed  Google Scholar 

  36. Kiyohara C, Horiuchi T, Takayama K, Nakanishi Y. IL1B rs1143634 polymorphism, cigarette smoking, alcohol use, and lung cancer risk in a Japanese population. J Thorac Oncol. 2010;5(3):299–304.

    Article  PubMed  Google Scholar 

  37. Omair A, Lie BA, Reikeras O, Brox JI. An association study of interleukin 18 receptor genes (IL18R1 and IL18RAP) in lumbar disc degeneration. Open Orthop J. 2012;6(1):164–71.

    Article  PubMed  PubMed Central  Google Scholar 

  38. González JR, et al. SNPassoc: an R package to perform whole genome association studies. Bioinformatics. 2007;23(5):644–5.

    Article  PubMed  Google Scholar 

  39. Eskola PJ, et al. Gender difference in genetic association between IL1A variant and early lumbar disc degeneration: a three-year follow-up. Int J Mol Epidemiol Genet. 2012;3(3):195–204.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Rajasekaran S, et al. Phenotype variations affect genetic association studies of degenerative disc disease: conclusions of analysis of genetic association of 58 single nucleotide polymorphisms with highly specific phenotypes for disc degeneration in 332 subjects. Spine J. 2013;13(10):1309–20.

    Article  CAS  PubMed  Google Scholar 

  41. Solovieva S, et al. Interleukin 1 polymorphisms and intervertebral disc degeneration. Epidemiology. 2004;15(5):626–33.

    Article  PubMed  Google Scholar 

  42. Zhu Y, et al. Association between IL4, IL6 gene polymorphism and lumbar disc degeneration in Chinese population. Oncotarget. 2017;8(51):89064–71.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Grocott HP. Genetics and outcome after cardiac surgery. J Extra Corpor Technol. 2006;38(1):49–51.

    PubMed  PubMed Central  Google Scholar 

  44. Bittar MN, et al. Interleukin 6 G-174C polymorphism influences outcome following coronary revascularization surgery. Heart Surg Forum. 2005;8(3):167–72.

    Article  Google Scholar 

  45. Yokomuro H, Yoshihara K, Watanabe Y, Shiono N, Koyama N, Takanashi Y. The variations in the immunologic features and interleukin-6 levels for the surgical treatment of cardiac myxomas. Surg Today. 2007;37(9):750–3.

    Article  CAS  PubMed  Google Scholar 

  46. Stoica AL, Stoica E, Constantinescu I, Uscatescu V, Ginghina C. Interleukin-6 and interleukin-10 gene polymorphism, endothelial dysfunction, and postoperative prognosis in patients with peripheral arterial disease. J Vasc Surg. 2010;52(1):103–9.

    Article  PubMed  Google Scholar 

  47. Wang Y, Che M, Xin J, Zheng Z, Li J, Zhang S. The role of IL-1β and TNF-α in intervertebral disc degeneration. Biomedicine and Pharmacotherapy. 2020;131. Elsevier Masson SAS, Nov. 01.

  48. Huang X, et al. Interleukin 6 (IL-6) and IL-10 promoter region polymorphisms are associated with risk of lumbar disc herniation in a northern Chinese Han population. Genet Test Mol Biomarkers. 2017;21(1):17–23.

    Article  CAS  PubMed  Google Scholar 

  49. Molinos M, Almeida CR, Caldeira J, Cunha C, Gonçalves RM, Barbosa MA. Inflammation in intervertebral disc degeneration and regeneration. Journal of the Royal Society Interface. 2015;12(104).

  50. Jin H, et al. Baicalein inhibits the IL-1β-induced inflammatory response in nucleus Pulposus cells and attenuates disc degeneration in vivo. Inflammation. 2019;42(3):1032–44.

    Article  CAS  PubMed  Google Scholar 

  51. Studer RK, Vo N, Sowa G, Ondeck C, Kang J. Human nucleus pulposus cells react to IL-6: independent actions and amplification of response to IL-1 and TNF-α. Spine (Phila Pa 1976). Apr. 2011;36(8):593–9.

    Article  Google Scholar 

  52. Wang P, Zhu F, Konstantopoulos K. Interleukin-6 synthesis in human chondrocytes is regulated via the antagonistic actions of prostaglandin (PG)E 2 and 15-deoxy-δ 12,14-PGJ 2. PLoS One. 2011;6(11):1–13.

    Article  CAS  Google Scholar 

  53. Omair A, Holden M, Lie BA, Reikeras O, Brox JI. Treatment outcome of chronic low back pain and radiographic lumbar disc degeneration are associated with inflammatory and matrix degrading gene variants: a prospective genetic association study. BMC Musculoskelet Disord. 2013;14(1):1.

    Article  CAS  Google Scholar 

  54. Deng X, Zhao F, Kang B, Zhang X. Elevated interleukin-6 expression levels are associated with intervertebral disc degeneration. Exp Ther Med. 2016;11(4):1425–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Kraychete DC, Sakata RK, Issy AM, Bacellar O, Santos-Jesus R, Carvalho EM. Serum cytokine levels in patients with chronic low back pain due to herniated disc: analytical cross-sectional study. Sao Paulo Med J. 2010;128(5):259–62.

    Article  PubMed  Google Scholar 

  56. Jimbo K, Park JS, Yokosuka K, Sato K, Nagata K. Positive feedback loop of interleukin-1beta upregulating production of inflammatory mediators in human intervertebral disc cells in vitro. J Neurosurg Spine. 2005;2(5):589–95.

    Article  PubMed  Google Scholar 

  57. Zhao B, Li X, Li R. Genetic relationship between IL-6 rs1800796 polymorphism and susceptibility to periodontitis. Immunol Investig. 2019;48(3):268–82.

    Article  CAS  Google Scholar 

  58. Kardol-Hoefnagel T, et al. A single nucleotide C3 polymorphism associates with clinical outcome after lung transplantation. Front Immunol. 2019;10, no. SEP:1–9.

  59. Zuo X, et al. Interleukin gene polymorphisms in Chinese Han population with breast cancer, a case-control study. Oncotarget. 2018;9(26):17994–8001.

    Article  PubMed  Google Scholar 

  60. Li W, et al. Gastric cancer-derived mesenchymal stromal cells trigger M2 macrophage polarization that promotes metastasis and EMT in gastric cancer. Cell Death Dis. 2019;10(12).

  61. Karstoft K, Pedersen BK. Exercise and type 2 diabetes: focus on metabolism and inflammation. Immunol Cell Biol. 2016;94(2):146–50.

    Article  CAS  PubMed  Google Scholar 

  62. Navarro-González JF, Mora-Fernández C, De Fuentes MM, García-Pérez J. Inflammatory molecules and pathways in the pathogenesis of diabetic nephropathy. Nat Rev Nephrol. 2011;7(6):327–40.

    Article  CAS  PubMed  Google Scholar 

  63. Suri P, et al. “U. S . Department of Veterans Affairs,”. 2018;17(1):4–14. Modifiable.

  64. Robertson D, Kumbhare D, Nolet P, Srbely J, Newton G. Associations between low back pain and depression and somatization in a Canadian emerging adult population. J Can Chiropr Assoc. 2017;61(2):96–105.

    PubMed  PubMed Central  Google Scholar 

  65. Dantzer R. NEUROIMMUNE interactions: from the BRAIN to the immune system and vice versa. Physiol Rev. 2018;98:477–504.

    Article  CAS  PubMed  Google Scholar 

  66. Ellul P, Boyer L, Groc L, Leboyer M, Fond G. Interleukin-1 β-targeted treatment strategies in inflammatory depression: toward personalized care. Acta Psychiatr Scand. 2016;134(6):469–84.

    Article  CAS  PubMed  Google Scholar 

  67. Miller AH, Maletic V, Raison CL. Inflammation and Its Discontents: The Role of Cytokines in the Pathophysiology of Major Depression. Biol Psychiatry. 2009;65(9):732–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Leonard B, Maes M. Mechanistic explanations how cell-mediated immune activation, inflammation and oxidative and nitrosative stress pathways and their sequels and concomitants play a role in the pathophysiology of unipolar depression. Neurosci Biobehav Rev. Feb. 2012;36(2):764–85.

    Article  CAS  PubMed  Google Scholar 

  69. Raison CL, et al. “A randomized controlled trial of the tumor necrosis factor antagonist infliximab for treatment-resistant depression: the role of baseline inflammatory biomarkers,” Arch Gen Psychiatry. 2013;70(1):31–41.

  70. Bialek K, et al. Preliminary study of the impact of single-nucleotide polymorphisms of IL-1α, IL-1β and TNF-α genes on the occurrence, severity and treatment effectiveness of the major depressive disorder. Cell Mol Neurobiol. 2020;40(6):1049–56.

    Article  PubMed  Google Scholar 

  71. Yu YWY, Chen TJ, Hong CJ, Chen HM, Tsai SJ. Association study of the interieukin-Ibeta (C-5I IT) genetic polymorphism with major depressive disorder, associated symptomatology, and antidepressant response. Neuropsychopharmacology. 2003;28(6):1182–5.

    Article  CAS  PubMed  Google Scholar 

  72. Borkowska P, et al. Interleukin-1beta promoter (-31T/C and -511C/T) polymorphisms in major recurrent depression. J Mol Neurosci. 2011;44(1):12–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Chen H, et al. Single nucleotide polymorphisms in the human interleukin-1B gene affect transcription according to haplotype context. Hum Mol Genet. 2006;15(4):519–29.

    Article  CAS  PubMed  Google Scholar 

  74. Baune BT, et al. The interleukin 1 Beta (IL1B) gene is associated with failure to achieve remission and impaired emotion processing in major depression. Biol Psychiatry. 2010;67(6):543–9.

    Article  CAS  PubMed  Google Scholar 

  75. Lončar Z, Čurić G, Havelka Meštrović A, Mićković V, Bilić M. Do IL-1B and IL-1RN modulate chronic low back pain in patients with post-traumatic stress disorder? Coll Antropol. 2013;37(4):1237–44.

    PubMed  Google Scholar 

  76. Ridout KK, et al. Interleukin 1B gene (IL1B) variation and internalizing symptoms in maltreated preschoolers. Dev Psychopathol. 2014;26:1277–87.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Kapelski P, et al. An association between functional polymorphisms of the interleukin 1 gene complex and schizophrenia using transmission disequilibrium test. Arch Immunol Ther Exp. 2016;64:161–8.

    Article  CAS  Google Scholar 

  78. Bower JE, Ganz PA, Irwin MR, Castellon S, Arevalo J, Cole SW. Cytokine genetic variations and fatigue among patients with breast cancer. J Clin Oncol. 2013;31(13):1656–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Nobis A, Zalewski D, Waszkiewicz N. Peripheral Markers of Depression. J Clin Med. 2020;9(12):3793.

  80. Won E, Kim YK. Neuroinflammation-associated alterations of the brain as potential neural biomarkers in anxiety disorders. Int J Mol Sci. 2020;21(18):1–19.

  81. Kovacs D, et al. Interleukin-6 promoter polymorphism interacts with pain and life stress influencing depression phenotypes. J Neural Transm. 2016;123(5):541–8.

    Article  CAS  PubMed  Google Scholar 

  82. Ting EY, Yang AC, Tsai SJ. “Role of interleukin-6 in depressive disorder,” Int J Mol Sci. 2020;21(6).

  83. Zhou YQ, et al. Interleukin-6: an emerging regulator of pathological pain. J Neuroinflammation. 2016;13(1):1–9.

    Article  CAS  Google Scholar 

  84. Sainoh T, et al. Single intradiscal injection of the interleukin-6 receptor antibody tocilizumab provides short-term relief of discogenic low back pain; prospective comparative cohort study. J Orthop Sci. 2016;21(1):2–6.

    Article  PubMed  Google Scholar 

Download references


Not applicable

GENODISC Consortium

4Dr Jill Urban and Prof Jeremy Fairbank, University of Oxford UK, and Chris Heywood, Team Rehab UK Ltd, Northants, Uk; Chris Heywood,

5Prof.Sarit Sivan Technion - Israel Institute of Technology, Haifa, Israel; current address ORT Braude College, Karmi’el, Israel.,

6Prof Sally Roberts Robert Jones and Agnes Hunt Orthopaedic Hospital, Oswestry Shropshire, UK,

7Dr Corneia Neidlinger-Wilke, Institute of Orthopaedc Research and Biomechanics, University of Ulm, Germany,

8Prof Jaakko Kaprio and Prof Michele Crites Battie, University of Helsinki, Finland. *current address; School of Physical Therapy,Western University, London, Ontario, mbattie@

9Prof Dimitris Kletsas, National Centre for Scientific Research ‘Demokritos’ Athens, Greece,

10Prof Keita Ito, and Prof Jaques Huyghe, Technische Universiteit Eindhoven, Eindhoven, The Netherlands

11Dr. Marco Brayda-Bruno, Istituto Ortopedico Galeazzi, Rome, Italy.

12Dr Nevenka Kregar Velikonja. Educell d.o.o Ljubljana, Slovenia,

1,3Dr Aron Lazary*, National Center for Spinal Disorders, Budapest, Hungary,

*Cosortium representative

Code availability

Not applicable


Open access funding provided by Semmelweis University. The research leading to these results received funding from the European Community’s Seventh Framework Programme (FP7, 2007–2013) under grant agreement no. HEALTH-F2-2008-201626 (GENODISC project).

Author information

Authors and Affiliations




All authors contributed to the study conception and design. Material preparation, data collection and genotyping were done by the members of GENODISC consortium. The first draft of the manuscript was written by Adam Biczo and Aron Lazary and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Aron Lazary.

Ethics declarations

Ethics approval and consent to participate

The research protocol was in accordance with the Helsinki Declaration of 1975 on human subjects testing revised in 1996. The study was approved by the Scientific and Research Ethics Committee of the Medical Research Council Hungary (431/PI/2007)

Subjects were included after signing a written informed consent. The study was approved by the Scientific and Research Ethics Committee of the Medical Research Council Hungary (431/PI/2007).

Consent for publication

Not applicable

Competing interests

The authors declare they have no financial interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Biczo, A., Bereczki, F., Koch, K. et al. Genetic variants of interleukin 1B and 6 are associated with clinical outcome of surgically treated lumbar degenerative disc disease. BMC Musculoskelet Disord 23, 774 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: