Skip to main content

Genetic susceptibility of intervertebral disc degeneration among young Finnish adults



Disc degeneration (DD) is a common condition that progresses with aging. Although the events leading to DD are not well understood, a significant genetic influence has been found. This study was undertaken to assess the association between relevant candidate gene polymorphisms and moderate DD in a well-defined and characterized cohort of young adults. Focusing on young age can be valuable in determining genetic predisposition to DD.


We investigated the associations of existing candidate genes for DD among 538 young adults with a mean age of 19 belonging to the 1986 Northern Finland Birth Cohort. Nineteen single nucleotide polymorphisms (SNP) in 16 genes were genotyped. We evaluated lumbar DD using the modified Pfirrmann classification and a 1.5-T magnetic resonance scanner for imaging.


Of the 538 individuals studied, 46% had no degeneration, while 54% had DD and 51% of these had moderate DD. The risk of DD was significantly higher in subjects with an allele G of IL6 SNPs rs1800795 (OR 1.45, 95% CI 1.07-1.96) and rs1800797 (OR 1.37, 95% CI 1.02-1.85) in the additive inheritance model. The role of IL6 was further supported by the haplotype analysis, which resulted in an association between the GGG haplotype (SNPs rs1800797, rs1800796 and rs1800795) and DD with an OR of 1.51 (95% CI 1.11-2.04). In addition, we observed an association between DD and two other polymorphisms, SKT rs16924573 (OR 0.27 95% CI 0.07-0.96) and CILP rs2073711 in women (OR 2.04, 95% CI 1.07-3.89).


Our results indicate that IL6, SKT and CILP are involved in the etiology of DD among young adults.

Peer Review reports


Disc degeneration (DD) is a multifactorial condition that progresses with aging. With age, histological, metabolic and functional changes take place in the disc, making it stiffer and reducing its strength [1]. Although the events leading to DD are not well understood, the outcome of studies over the last decade has shown that genetic influence plays a large role in DD, together with environmental factors [2]. Known environmental risk factors for DD include body mass index (BMI) and heavy physical loading [3], although their effect is weak in comparison to heredity [4].

To date, several gene loci associated with human disc degeneration have been identified [5]. Variations in the genes involved in inflammation, extracellular matrix components, and protein metabolism have been reported as associating with DD [6, 7]. The first polymorphisms associated with DD were two variations in the Vitamin D receptor gene (TaqI and FokI) [8]. Subsequently, associations with aggrecan (ACAN) and two collagen IX genes (COL9A2, COL9A3) were reported [911]. The structural components of the intervertebral disc have been an obvious choice as candidate genes, because variations in structural genes may render discs more vulnerable to mechanical loading, causing injuries to endplates and/or annulus fibrosus, leading to loss of disc function [5]. In addition, variants in genes encoding cytokines, such as interleukin-1 (IL1) and interleukin-6 (IL6), or degradative metalloproteinases (MMPs) have been reported as being associated with DD [1215]. The genes of adipocytokines are also putative candidates for DD due to the association between DD and atherosclerosis [16]. However, of the reported associations, only a few have been verified in different ethnic populations [2].

Although genetic analyses have gained a more central role in the exploration of the underlying causes of DD, its etiology and pathogenesis still remain unclear. In general, genetic studies have suffered from small sample sizes and differences in the definition and characterization of DD. Here we investigated the association between relevant candidate gene polymorphisms and moderate DD in a well-defined and characterized cohort of young adults. Focusing on moderate DD phenotype manifesting at an early age will provide a better opportunity to determine the predisposing role of genetic polymorphisms.


Study population

The study population consisted of members of the 1986 Northern Finland Birth Cohort (NFBC 1986) with an expected date of birth between July 1, 1985 and June 30, 1986 (n = 9479) in the two northernmost provinces of Finland; Oulu and Lapland. The 2969 that lived within 100 km of the city of Oulu received a postal questionnaire in 2003. Altogether, 1987 responded, with a mean age of 18 (response rate 68%) and were invited to participate in a physical examination, which took place in 2005-2006 at a mean age of 19. The participants of the examination (n = 874) were invited to participate in lumbar spine magnetic resonance imaging (MRI) [17]. Our study population consisted of 538 subjects with both DNA sample and lumbar spine MRI. The Ethics Committee of the Northern Ostrobothnia Hospital District approved the study plan, and the study was performed according to the Declaration of Helsinki. Written informed consent was obtained from all individuals.

Magnetic resonance imaging

MRI scans were performed using a 1.5-T unit (Signa, General Electric, Milwaukee, WI) and phased array spine coil (USA Instruments, Aurora OH, USA) with the imaging protocol of a sagittal T1-weighted (440/14 [repetition time msec/echo time msec]) spin echo, a sagittal T2-weighted (3960/116) fast spin echo (FSE), and the axial T2-weighted FSE (5160/107) of the entire lumbar spine.

The degree of DD was graded from T2-weighted images using a modified Pfirrmann classification: Grade 1 (normal shape, no horizontal bands, distinction of nucleus, and annulus is clear), Grade 2 (nonhomogeneous shape with horizontal bands, some blurring between nucleus and anulus), Grade 3 (nonhomogeneous shape with blurring between nucleus and anulus, anulus shape still recognizable), Grade 4 (nonhomogeneous shape with hypointensity, anulus shape not intact and distinction between nucleus and anulus impossible, disc height usually decreased), and Grade 5 (same as Grade 4 but collapsed disc space) [17, 18]. In the modified Pfirrmann classification, hyperintensity or isointensity of the intervertebral disc to cerebrospinal fluid (CSF) ratios were not used as criteria for Grades 1 through 3, due to the fact that CSF was hyperintense compared to the discs in all the MR sequences. Moderate DD was defined as Grade 4 degeneration at one or more lumbar levels, or Grade 3 degeneration at three or more levels. Subjects with Grade 1 and Grade 2 findings and no other pathology in MRI were defined as the reference group, i.e. subjects without DD. Information on inter-rater reliability and comparison of the MRI study population and the rest of the survey respondents have been published elsewhere [17].

Candidate gene and single-nucleotide polymorphism (SNP) selection

At the time of the study conception we listed all putative DD candidate genes and polymorphisms. The current genotyping assay of choice is designed for SNP genotyping thus excluding other types of genetic variations from this study. After testing different polymerase chain reaction (PCR) multiplexes we were able to include nineteen SNPs in five separate multiplexes. These SNPs were distributed in 16 different genes (Table 1). Due to the structure of the DNA sequences at the SNP sites we were not able to include all possible polymorphisms in the analyses. Four genes with conceivable relevant biological importance and/or preliminary association in other populations (unpublished results) to the phenotype in question were included in the study (ADIPOQ: rs2241766, COL22A1:rs2292927, LEPR: rs1137100 and rs8179183, and MMP10:rs470154).

Table 1 Candidate gene information

SNP genotyping

We isolated DNA from frozen whole blood and performed the genotyping using a single base extension fluorescent method, the SNaPshot Multiplex Kit (Applied Biosystems). PCR amplification and SNaPshot extension reactions were both designed to have the same annealing temperature. The SAP/ExoI purified PCR product was processed according to the manufacturer's protocol with the exception of the halved amount (2.5 μl) of SNaPshot Ready Reaction Mix. We carried out analysis using the ABI3031xl Genetic Analyzer (Applied Biosystems), and determined the genotypes using GeneMapper software (Applied Biosystems). We validated the genotype results from a randomly chosen 10% of samples for each SNP by direct sequencing, using the abovementioned genetic analyzer with BigDye chemistry from the same manufacturer. Genotyping was blinded towards MRI status.

Statistical Analysis

We tested the potential deviation from the Hardy-Weinberg equilibrium (HWE) using the chi-square test. Prior to polymorphism association, the possible effect of gender, BMI and height on the DD status was analyzed using the SAS program (V 9.2). Crude and adjusted ORs and their 95% confidence intervals (CIs) were calculated using logistic regression with SAS. The dependent variable was DD phenotype and independent variables were genotype, allele or haplotype. Haplotypes were first estimated by SNPStats [19] and then statistically reconstructed from population genotype data, using the PHASE program and the Markov-chain method for haplotype assignments [20]. In candidate gene association studies, the questions are more specific than in whole genome association studies. Therefore, in the setting of existing a priori hypotheses, multiple testing corrections were only accounted for in genes with multiple SNPs. The Bonferroni-Holm method was used for multiple testing corrections [21].


Clinical results

Of the 538 subjects in this study, 45.7% (N = 246) did not have DD, whereas 54.3% (N = 292) had DD. Of the 292 individuals with DD, 51.4% (N = 150) had moderate DD as defined in the method section. The remaining 142 individuals had degenerative findings but were not included in the analysis. Height and BMI were associated with DD (Table 2), but after adjustment for height, BMI was no longer associated with DD. DD was more frequent among men (32.1%) than women (24.9%), but the difference was not statistically significant when adjusted for height.

Table 2 Characteristics of study population at the age of 19

Genetic associations

Overall, genotyping of all the subjects from the selected 19 SNPs was 98% successful. All the genotype frequencies analyzed were in HWE. Associations were observed in three of the 16 genes tested. Table 3 presents the genotype frequencies and respective ORs for individual SNP. The strongest associations with DD were observed for polymorphisms in IL6 (Table 3). The IL6 rs1800795 G-allele was associated with the degenerative phenotype in an additive manner; GC-genotype, OR 1.51 [0.92-2.47] and GG-genotype, OR 2.09 [1.14-3.85]. Minor allele (G) frequencies were 0.490 and 0.398 in subjects with and without DD, respectively (OR 1.45 [1.07-1.96]). The rs1800797 was similarly associated with the DD phenotype, suggesting an additive model (GA: OR 1.55 [0.95-2.54]), GG: OR 1.85 [1.01-3.39]), compared to the AA-genotype. G-allele frequencies were 0.490 and 0.409 among subjects with and without DD, respectively (OR 1.37 [1.02-1.85]). Subsequent haplotype analyses resulted in an association between the GGG haplotype (from rs1800797, rs1800796 and rs1800795, respectively) and DD (OR 1.48 [1.09-2.01]) (Table 4). From this haplotype, only one SNP (rs1800796) did not individually associate with DD. The effect of the GGG haplotype was additive, while the most frequent haplotype irrespective of the MRI findings was AGC.

Table 3 Genotype frequencies of the analyzed SNPs
Table 4 Frequencies of IL6 haplotypes and diplotypes

The GA genotype of the SKT (KIAA1217) polymorphism rs16924573 was associated with a reduced risk of DD (OR 0.27 [0.07-0.96], p = 0.024). The frequency of the GA genotype was lower among subjects with DD (2.0%), compared to subjects without DD (6.5%). The CILP SNP rs2073711 was associated with DD, but only among women. The CT/TT genotypes, compared to the CC genotype, resulted in an OR of 2.04 [1.07-3.89], and p = 0.025 (corrected p = 0.05). We did not observe statistically significant associations between the other SNPs and DD.


Disc degeneration is described as a common multifactorial and multigenic condition [3]. At least 14 genes have been shown to associate with DD, although their specific effects are not known. In this study, IL6, SKT and CILP were associated with moderate DD in a sample population of young adults who had their MRI taken at a mean age of 21 years.

IL6 is a pleiotropic cytokine which is expressed in many cells and is one of the most important mediators of inflammatory reactions in humans [22]. The role of IL6 and in particular the effect of three common promoter polymorphisms (rs1800797, rs1800796, rs1800795) has been studied previously for the presence of sciatica [23], but also in relation to DD among 12- to 14-year-old children [15]. In our population of young adults rs1800797 and rs1800795 were associated with moderate DD. The G allele of both rs1800797 and rs1800795 was the risk allele for DD (OR 1.37 [1.02-1.85] and 1.45 [1.07-1.96], respectively).

The role of IL6 was further strengthened with haplotype analyses of the promoter SNPs. IL6 haplotype GGG was associated with moderate DD (OR 1.48 [1.09-2.01]). In addition, the risk for moderate DD was higher for the carriers of the GGG/GGG diplotype (OR 2.22 [1.17-4.22]). This finding is supported by an observation that, among sciatica patients, GGGA was the risk haplotype (OR 4.80 [1.59-14.45]) [23]. Furthermore, the GGGA haplotype was a predisposing factor for future sciatica symptoms [24]. In a study of young Danish children (aged between 12 and 14) a GCG haplotype was associated with early DD (OR 6.46 [1.61-26]) [15]. This association was present only in females, and the population with DD was small, consisting of 68 individuals. The two most common haplotypes observed in that study were in accordance with the ones we observed in ours (AGC and GGG), indicating that there are no major differences in genotype frequencies between the two study populations.

Functional analyses have shown that the production of IL6 is higher in individuals carrying the GG genotypes for rs1800797 and rs1800795 variations [25]. It has also been reported that the three common promoter polymorphisms influence the transcription of IL6 with the G allele by each being responsible for the increased overall transcription of IL6 [26]. Overproduction of IL6 has been associated with a spectrum of age-related conditions, including atherosclerosis, peripheral vascular disease, coronary artery disease, and osteoporosis [27]. It is interesting to note that the same common IL6 promoter variants are also associated with the most severe forms of distal interphalangeal osteoarthritis [28]. Osteoarthritis and DD share pathological and clinical features, even though the extent to which the etiological determinants of these diseases' processes are similar is not known [29].

The association of the SKT rs16924573 variant with DD has not been previously observed. In our population, carriers of the minor allele A, which is a rare allele in the Finnish population, were at a reduced risk of DD. The importance of SKT in development of caudal intervertebral discs (IVD) in mice has been established previously [30]. The mice with homozygotely trapped SKT gene produce a deformed "kinky tail" phenotype - seemingly due to nucleus pulposus changes - during adulthood [30]. These observations led to the first study with human lumbar disc herniation patients characterized by sciatica [31]. In the study an association of SKT rs16924573 G allele with lumbar disc herniations was found (OR 1.34 [1.14-1.58]) [31]. DD, phenotype of the present study, is believed to be an important intermediate stage in symptomatic disc disease [32], which is often characterized by disc herniations either with or without sciatic pain radiating from the back to the leg [33]. Our current finding is in line with the earlier findings, however, the differences between the phenotypes in these studies need to be noted. These findings provide some support to the hypothesis that SKT could have long-term importance in the onset of symptomatic disc disease by making the discs more vulnerable for DD.

In addition, an association between the CILP variant and DD was observed. CILP encodes the cartilage intermediate protein, which inhibits the transforming growth factor β1- (TGF-β1) mediated induction of extracellular matrix proteins through direct interaction with TGF-β1 [34]. In the present study, the T allele was associated with an increased risk of DD. The TT/CT genotype compared to the CC genotype showed an OR of 2.04 [1.07-3.89]. The association between CILP and DD was, however, only found in females. This finding is in contrast with previously reported associations among Japanese subjects, which indicated the C allele as the risk allele for symptomatic lumbar disc disease [34] as well as for lumbar DD [35, 36]. However, a replication study in Finnish and Chinese populations did not find an association between CILP and symptomatic disc disease or MRI-defined DD [37]. Interestingly, the C allele was associated with reduced risk of knee OA [38], which is in accordance with our findings that the T allele, and not the C allele, is the risk allele. Although it was shown in functional studies that the C allele is responsible for increased binding and inhibition of TGF-β1 [34], the discrepancies in the association studies could simply be due to the difference in allele frequencies between European and Japanese populations

In our population, both men and women with DD were significantly taller than those without DD. The association between body height and herniated lumbar disc has previously been contradictory [39, 40]. In our earlier study, height was associated with DD on MRI among 12- to 14-year-old Danish children (N = 220) [15]. In another study among middle-aged women living in the UK (N = 796), individuals with radiographically defined DD progression over a 9-year period were taller than individuals without DD progression, but the difference (0.7 cm) was not statistically significant [41]. However, in a recent large cohort study (N = 13, 282) body height recorded at baseline was associated with higher likelihood of back surgery and low back pain during the 12-year follow-up. Although the associations reached statistical significance only among males, a similar trend was observed among females as well [42].

One possible mechanism for the role of body height as a risk factor is lumbar disc size. Larger discs might be predisposed to DD due to greater diffusion distances of nutrients into the centre of the nucleus pulposus [43]. Furthermore, higher discs might be at an increased risk of mechanical failure under external loading [44]. Finally, facet joint alterations have been found to be more frequent in taller individuals with lumbar disc herniation [45].

To date there have been many reported genetic associations in relation to DD, yet only a few have been replicated in different populations. We were not able to detect any significant associations for the majority of the polymorphisms that have been previously implicated in DD. This discrepancy may be due to differences in the definitions of the DD phenotype, sample characteristics, or true differences in the prevalence of polymorphisms among ethnic populations. The majority of these studies have been mainly conducted in adult populations, and only a few have been performed specifically in younger populations. It is well possible that a particular gene is associated with DD at only a particular age; for instance, a 5A/6A promoter polymorphism of the MMP3 gene associated with more degenerative findings in the elderly with a mean age of 73 than in younger subjects with a mean age of 21.4 [46]. It has also been demonstrated that the Trp2 allele of COL9A2 is a significant age-dependent risk factor for the development and severity of DD, increasing the risk of DD in the 40 to 49 year age group in the Trp2 carriers [47].

The strength of our study is the young age of the participants of this population-based birth cohort. The effects of environmental factors are most likely not so evident, whereas genetic factors may play a more important role in comparison to adult populations. Furthermore, in order to avoid any bias, our genotypes were defined a priori without any knowledge of the results of MRI readings. Our DD phenotype was also defined a priori before the genetic analyses and three specialists evaluated the MR images independently. Despite the initial large number of MRI scans, our sample study size is still fairly modest for studying multifactorial and multigenic conditions. A further limitation in the study was that only one SNP was included for the majority of the genes. For a comprehensive evaluation of the gene, a more extensive SNP analyses across the gene regions would be needed.

DD is a common disorder and is believed to be an important intermediate stage in symptomatic disc disease [32]. Several recent studies suggest that DD is associated with low back symptoms [4850]. The findings indicate that a higher degree of lumbar DD is related to a higher likelihood of symptoms, and the presence of moderate DD or DD at multiple levels increases the likelihood of pain [50]. Therefore, the clinical relevance of genetic studies exploring the pathomechanisms of DD is indisputable.


In summary, five variations of the 19 selected candidate SNPs were associated with susceptibility to DD in our study. Our results support the role of IL6, SKT and CILP in DD. Furthermore, we were able to show the associations in a population of young adults. Overall, genetic association studies provide information about specific gene variants on disease occurrence, but interpreting their results has proven to be a challenge [29]. The lack of association of well-known variations observed here does not necessarily contradict previous reporting. Instead our results emphasize the complexity of DD, and the important role of further studies in populations with different ethnic and age backgrounds.


  1. 1.

    Adams MA, Roughley PJ: What is intervertebral disc degeneration, and what causes it?. Spine (Phila Pa 1976). 2006, 31 (18): 2151-2161. 10.1097/01.brs.0000231761.73859.2c.

    Article  Google Scholar 

  2. 2.

    Battie MC, Videman T, Kaprio J, Gibbons LE, Gill K, Manninen H, Saarela J, Peltonen L: The Twin Spine Study: contributions to a changing view of disc degeneration. Spine J. 2009, 9 (1): 47-59. 10.1016/j.spinee.2008.11.011.

    Article  PubMed  Google Scholar 

  3. 3.

    Battie MC, Videman T: Lumbar disc degeneration: epidemiology and genetics. J Bone Joint Surg Am. 2006, 88 (Suppl 2): 3-9.

    Article  PubMed  Google Scholar 

  4. 4.

    Battie MC, Videman T, Levalahti E, Gill K, Kaprio J: Genetic and environmental effects on disc degeneration by phenotype and spinal level: a multivariate twin study. Spine (Phila Pa 1976). 2008, 33 (25): 2801-2808. 10.1097/BRS.0b013e31818043b7.

    Article  Google Scholar 

  5. 5.

    Chan D, Song Y, Sham P, Cheung KM: Genetics of disc degeneration. Eur Spine J. 2006, 15 (Suppl 3): S317-25.

    Article  PubMed  Google Scholar 

  6. 6.

    Ala-Kokko L: Genetic risk factors for lumbar disc disease. Ann Med. 2002, 34 (1): 42-47. 10.1080/078538902317338634.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Kalichman L, Hunter DJ: The genetics of intervertebral disc degeneration. Associated genes. Joint Bone Spine. 2008, 75 (4): 388-396. 10.1016/j.jbspin.2007.11.002.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Videman T, Leppavuori J, Kaprio J, Battie MC, Gibbons LE, Peltonen L, Koskenvuo M: Intragenic polymorphisms of the vitamin D receptor gene associated with intervertebral disc degeneration. Spine. 1998, 23 (23): 2477-2485. 10.1097/00007632-199812010-00002.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Kawaguchi Y, Osada R, Kanamori M, Ishihara H, Ohmori K, Matsui H, Kimura T: Association between an aggrecan gene polymorphism and lumbar disc degeneration. Spine (Phila Pa 1976). 1999, 24 (23): 2456-2460. 10.1097/00007632-199912010-00006.

    CAS  Article  Google Scholar 

  10. 10.

    Annunen S, Paassilta P, Lohiniva J, Perala M, Pihlajamaa T, Karppinen J, Tervonen O, Kroger H, Lahde S, Vanharanta H, Ryhanen L, Goring HH, Ott J, Prockop DJ, Ala-Kokko L: An allele of COL9A2 associated with intervertebral disc disease. Science. 1999, 285 (5426): 409-412. 10.1126/science.285.5426.409.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Paassilta P, Lohiniva J, Goring HH, Perala M, Raina SS, Karppinen J, Hakala M, Palm T, Kroger H, Kaitila I, Vanharanta H, Ott J, Ala-Kokko L: Identification of a novel common genetic risk factor for lumbar disk disease. JAMA. 2001, 285 (14): 1843-1849. 10.1001/jama.285.14.1843.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Solovieva S, Kouhia S, Leino-Arjas P, Ala-Kokko L, Luoma K, Raininko R, Saarela J, Riihimaki H: Interleukin 1 polymorphisms and intervertebral disc degeneration. Epidemiology. 2004, 15 (5): 626-633. 10.1097/01.ede.0000135179.04563.35.

    Article  PubMed  Google Scholar 

  13. 13.

    Karppinen J, Daavittila I, Solovieva S, Kuisma M, Taimela S, Natri A, Haapea M, Korpelainen R, Niinimaki J, Tervonen O, Ala-Kokko L, Mannikko M: Genetic factors are associated with modic changes in endplates of lumbar vertebral bodies. Spine (Phila Pa 1976). 2008, 33 (11): 1236-1241. 10.1097/BRS.0b013e318170fd0e.

    Article  Google Scholar 

  14. 14.

    Karppinen J, Solovieva S, Luoma K, Raininko R, Leino-Arjas P, Riihimaki H: Modic changes and interleukin 1 gene locus polymorphisms in occupational cohort of middle-aged men. Eur Spine J. 2009, 18 (12): 1963-1970. 10.1007/s00586-009-1139-x.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Eskola PJ, Kjaer P, Daavittila I, Solovieva S, Okuloff A, Sorensen JS, Wedderkopp Ni, Ala-Kokko L, Männikkö M, Karppinen J: Genetic risk factors of disc degeneration among 12-14-year-old Danish children: a population study. Int J Mol Epidemiol Genet. 2010, 1 (2): 158-165.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Kauppila LI: Atherosclerosis and Disc Degeneration/Low-Back Pain - A Systematic Review. Eur J Vasc Endovasc Surg. 2009, 37 (6): 661-670. 10.1016/j.ejvs.2009.02.006.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Takatalo J, Karppinen J, Niinimaki J, Taimela S, Nayha S, Jarvelin MR, Kyllonen E, Tervonen O: Prevalence of degenerative imaging findings in lumbar magnetic resonance imaging among young adults. Spine (Phila Pa 1976). 2009, 34 (16): 1716-1721. 10.1097/BRS.0b013e3181ac5fec.

    Article  Google Scholar 

  18. 18.

    Pfirrmann CW, Metzdorf A, Zanetti M, Hodler J, Boos N: Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine (Phila Pa 1976). 2001, 26 (17): 1873-1878. 10.1097/00007632-200109010-00011.

    CAS  Article  Google Scholar 

  19. 19.

    Sole X, Guino E, Valls J, Iniesta R, Moreno V: SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006, 22 (15): 1928-1929. 10.1093/bioinformatics/btl268.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68 (4): 978-989. 10.1086/319501.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Holm S: A simple sequentially rejective multiple test procedure. Scand J Statist. 1979, 65-70. 6

  22. 22.

    Samuel JM, Kelberman D, Smith AJ, Humphries SE, Woo P: Identification of a novel regulatory region in the interleukin-6 gene promoter. Cytokine. 2008, 42 (2): 256-264. 10.1016/j.cyto.2008.02.008.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Noponen-Hietala N, Virtanen I, Karttunen R, Schwenke S, Jakkula E, Li H, Merikivi R, Barral S, Ott J, Karppinen J, Ala-Kokko L: Genetic variations in IL6 associate with intervertebral disc disease characterized by sciatica. Pain. 2005, 114 (1-2): 186-194. 10.1016/j.pain.2004.12.015.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Karppinen J, Daavittila I, Noponen N, Haapea M, Taimela S, Vanharanta H, Ala-Kokko L, Mannikko M: Is the interleukin-6 haplotype a prognostic factor for sciatica?. Eur J Pain. 2008, 12 (8): 1018-1025. 10.1016/j.ejpain.2008.01.009.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Rivera-Chavez FA, Peters-Hybki DL, Barber RC, O'Keefe GE: Interleukin-6 promoter haplotypes and interleukin-6 cytokine responses. Shock. 2003, 20 (3): 218-223. 10.1097/00024382-200309000-00004.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Terry CF, Loukaci V, Green FR: Cooperative influence of genetic polymorphisms on interleukin 6 transcriptional regulation. J Biol Chem. 2000, 275 (24): 18138-44. 10.1074/jbc.M000379200.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Omoigui S: The Interleukin-6 inflammation pathway from cholesterol to aging--role of statins, bisphosphonates and plant polyphenols in aging and age-related diseases. Immun Ageing. 2007, 4: 1-10.1186/1742-4933-4-1.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Kamarainen OP, Solovieva S, Vehmas T, Luoma K, Riihimaki H, Ala-Kokko L, Mannikko M, Leino-Arjas P: Common interleukin-6 promoter variants associate with the more severe forms of distal interphalangeal osteoarthritis. Arthritis Res Ther. 2008, 10 (1): R21-10.1186/ar2374.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Ryder JJ, Garrison K, Song F, Hooper L, Skinner J, Loke Y, Loughlin J, Higgins JP, MacGregor AJ: Genetic associations in peripheral joint osteoarthritis and spinal degenerative disease: a systematic review. Ann Rheum Dis. 2008, 67 (5): 584-591.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Semba K, Araki K, Li Z, Matsumoto K, Suzuki M, Nakagata N, Takagi K, Takeya M, Yoshinobu K, Araki M, Imai K, Abe K, Yamamura K: A novel murine gene, Sickle tail, linked to the Danforth's short tail locus, is required for normal development of the intervertebral disc. Genetics. 2006, 172 (1): 445-456.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Karasugi T, Semba K, Hirose Y, Kelempisioti A, Nakajima M, Miyake A, Furuichi T, Kawaguchi Y, Mikami Y, Chiba K, Kamata M, Ozaki K, Takahashi A, Makela P, Karppinen J, Kimura T, Kubo T, Toyama Y, Yamamura K, Mannikko M, Mizuta H, Ikegawa S: Association of the tag SNPs in the human SKT gene (KIAA1217) with lumbar disc herniation. J Bone Miner Res. 2009, 24 (9): 1537-1543. 10.1359/jbmr.090314.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Videman T, Battie MC, Ripatti S, Gill K, Manninen H, Kaprio J: Determinants of the progression in lumbar degeneration: a 5-year follow-up study of adult male monozygotic twins. Spine (Phila Pa 1976). 2006, 31 (6): 671-678. 10.1097/01.brs.0000202558.86309.ea.

    Article  Google Scholar 

  33. 33.

    Frymoyer JW: Back pain and sciatica. N Engl J Med. 1988, 318 (5): 291-300. 10.1056/NEJM198802043180506.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Seki S, Kawaguchi Y, Chiba K, Mikami Y, Kizawa H, Oya T, Mio F, Mori M, Miyamoto Y, Masuda I, Tsunoda T, Kamata M, Kubo T, Toyama Y, Kimura T, Nakamura Y, Ikegawa S: A functional SNP in CILP, encoding cartilage intermediate layer protein, is associated with susceptibility to lumbar disc disease. Nat Genet. 2005, 37 (6): 607-612. 10.1038/ng1557.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Min SK, Nakazato K, Okada T, Ochi E, Hiranuma K: The cartilage intermediate layer protein gene is associated with lumbar disc degeneration in collegiate judokas. Int J Sports Med. 2009, 30 (9): 691-694. 10.1055/s-0029-1214380.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Min SK, Nakazato K, Yamamoto Y, Gushiken K, Fujimoto H, Fujishiro H, Kobayakawa Y, Hiranuma K: Cartilage intermediate layer protein gene is associated with lumbar disc degeneration in male, but not female, collegiate athletes. Am J Sports Med. 2010, 38 (12): 2552-2557. 10.1177/0363546510376714.

    Article  PubMed  Google Scholar 

  37. 37.

    Virtanen IM, Song YQ, Cheung KM, Ala-Kokko L, Karppinen J, Ho DW, Luk KD, Yip SP, Leong JC, Cheah KS, Sham P, Chan D: Phenotypic and population differences in the association between CILP and lumbar disc disease. J Med Genet. 2007, 44 (4): 285-288. 10.1136/jmg.2006.047076.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Valdes AM, Hart DJ, Jones KA, Surdulescu G, Swarbrick P, Doyle DV, Schafer AJ, Spector TD: Association study of candidate genes for the prevalence and progression of knee osteoarthritis. Arthritis Rheum. 2004, 50 (8): 2497-507. 10.1002/art.20443.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Heliovaara M: Risk factors for low back pain and sciatica. Ann Med. 1989, 21 (4): 257-264. 10.3109/07853898909149202.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Kelsey JL, Githens PB, O'Conner T, Weil U, Calogero JA, Holford TR, White AA, Walter SD, Ostfeld AM, Southwick WO: Acute prolapsed lumbar intervertebral disc. An epidemiologic study with special reference to driving automobiles and cigarette smoking. Spine (Phila Pa 1976). 1984, 9 (6): 608-613. 10.1097/00007632-198409000-00012.

    CAS  Article  Google Scholar 

  41. 41.

    Hassett G, Hart DJ, Manek NJ, Doyle DV, Spector TD: Risk factors for progression of lumbar spine disc degeneration: the Chingford Study. Arthritis Rheum. 2003, 48 (11): 3112-3117. 10.1002/art.11321.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Coeuret-Pellicer M, Descatha A, Leclerc A, Zins M: Are tall people at higher risk of low back pain surgery? A discussion on the results of a multipurpose cohort. Arthritis Care Res (Hoboken). 2010, 62 (1): 125-127. 10.1002/acr.20023.

    Article  Google Scholar 

  43. 43.

    Urban JPG, Smith S, Fairbank JCT: Nutrition of the intervertebral disc. Spine (Phila Pa 1976). 2004, 29 (23): 2700-2709. 10.1097/01.brs.0000146499.97948.52.

    Article  Google Scholar 

  44. 44.

    Natarajan RN, Andersson GB: The influence of lumbar disc height and cross-sectional area on the mechanical response of the disc to physiologic loading. Spine (Phila Pa 1976). 1999, 24 (18): 1873-1881. 10.1097/00007632-199909150-00003.

    CAS  Article  Google Scholar 

  45. 45.

    Karacan I, Aydin T, Sahin Z, Cidem M, Koyuncu H, Aktas I, Uludag M: Facet angles in lumbar disc herniation: their relation to anthropometric features. Spine (Phila Pa 1976). 2004, 29 (10): 1132-1136. 10.1097/00007632-200405150-00016.

    Article  Google Scholar 

  46. 46.

    Takahashi M, Haro H, Wakabayashi Y, Kawa-uchi T, Komori H, Shinomiya K: The association of degeneration of the intervertebral disc with 5a/6a polymorphism in the promoter of the human matrix metalloproteinase-3 gene. J Bone Joint Surg Br. 2001, 83 (4): 491-495. 10.1302/0301-620X.83B4.11617.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Jim JJ, Noponen-Hietala N, Cheung KM, Ott J, Karppinen J, Sahraravand A, Luk KD, Yip SP, Sham PC, Song YQ, Leong JC, Cheah KS, Ala-Kokko L, Chan D: The TRP2 allele of COL9A2 is an age-dependent risk factor for the development and severity of intervertebral disc degeneration. Spine. 2005, 30 (24): 2735-2742. 10.1097/01.brs.0000190828.85331.ef.

    Article  PubMed  Google Scholar 

  48. 48.

    Luoma K, Riihimaki H, Luukkonen R, Raininko R, Viikari-Juntura E, Lamminen A: Low back pain in relation to lumbar disc degeneration. Spine (Phila Pa 1976). 2000, 25 (4): 487-492. 10.1097/00007632-200002150-00016.

    CAS  Article  Google Scholar 

  49. 49.

    Cheung KM, Karppinen J, Chan D, Ho DW, Song YQ, Sham P, Cheah KS, Leong JC, Luk KD: Prevalence and pattern of lumbar magnetic resonance imaging changes in a population study of one thousand forty-three individuals. Spine (Phila Pa 1976). 2009, 34 (9): 934-940. 10.1097/BRS.0b013e3181a01b3f.

    Article  Google Scholar 

  50. 50.

    de Schepper EI, Damen J, van Meurs JB, Ginai AZ, Popham M, Hofman A, Koes BW, Bierma-Zeinstra SM: The association between lumbar disc degeneration and low back pain: the influence of age, gender, and individual radiographic features. Spine (Phila Pa 1976). 2010, 35 (5): 531-536. 10.1097/BRS.0b013e3181aa5b33.

    Article  Google Scholar 

  51. 51.

    Yang WS, Chuang LM: Human genetics of adiponectin in the metabolic syndrome. J Mol Med. 2006, 84 (2): 112-121. 10.1007/s00109-005-0011-7.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Noponen-Hietala N, Kyllonen E, Mannikko M, Ilkko E, Karppinen J, Ott J, Ala-Kokko L: Sequence variations in the collagen IX and XI genes are associated with degenerative lumbar spinal stenosis. Ann Rheum Dis. 2003, 62 (12): 1208-1214. 10.1136/ard.2003.008334.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Maeda S, Koga H, Matsunaga S, Numasawa T, Ikari K, Furushima K, Harata S, Takeda J, Sakou T, Komiya S, Inoue I: Gender-specific haplotype association of collagen alpha2 (XI) gene in ossification of the posterior longitudinal ligament of the spine. J Hum Genet. 2001, 46 (1): 1-4. 10.1007/s100380170117.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Sims AM, Timms AE, Bruges-Armas J, Burgos-Vargas R, Chou CT, Doan T, Dowling A, Fialho RN, Gergely P, Gladman DD, Inman R, Kauppi M, Kaarela K, Laiho K, Maksymowych W, Pointon JJ, Rahman P, Reveille JD, Sorrentino R, Tuomilehto J, Vargas-Alarcon G, Wordsworth BP, Xu H, Brown MA, International Genetics of Ankylosing Spondylitis: Prospective meta-analysis of interleukin 1 gene complex polymorphisms confirms associations with ankylosing spondylitis. Ann Rheum Dis. 2008, 67 (9): 1305-1309.

    Article  PubMed  Google Scholar 

  55. 55.

    Fairbrother UL, Tanko LB, Walley AJ, Christiansen C, Froguel P, Blakemore AI: Leptin receptor genotype at Gln223Arg is associated with body composition, BMD, and vertebral fracture in postmenopausal Danish women. J Bone Miner Res. 2007, 22 (4): 544-550. 10.1359/jbmr.070114.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Hirose Y, Chiba K, Karasugi T, Nakajima M, Kawaguchi Y, Mikami Y, Furuichi T, Mio F, Miyake A, Miyamoto T, Ozaki K, Takahashi A, Mizuta H, Kubo T, Kimura T, Tanaka T, Toyama Y, Ikegawa S: A functional polymorphism in THBS2 that affects alternative splicing and MMP binding is associated with lumbar-disc herniation. Am J Hum Genet. 2008, 82 (5): 1122-1129. 10.1016/j.ajhg.2008.03.013.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Cheung KM, Chan D, Karppinen J, Chen Y, Jim JJ, Yip SP, Ott J, Wong KK, Sham P, Luk KD, Cheah KS, Leong JC, Song YQ: Association of the Taq I allele in vitamin D receptor with degenerative disc disease and disc bulge in a Chinese population. Spine (Phila Pa 1976). 2006, 31 (10): 1143-1148. 10.1097/01.brs.0000216530.41838.d3.

    Article  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


We wish to thank Aira Erkkilä, Helena Satulehto and Niina Heikkinen for their excellent technical assistance.


Supported in part by grants from the Academy of Finland (121620 and 115471) and AOSPINE (AOSBRC-07-2).

Author information



Corresponding author

Correspondence to Minna Männikkö.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AK, PE and UK participated in analysis and interpretation, design and writing. JT, ID, JN, RBS, OT, SS and AO participated in analysis and interpretation and critical revising. PYPK, YQS, KMCC, DC, L-AK and M-RJ participated in conception, design and some in critical revising. JK and MM participated in conception, design of the study, writing, critical revisions and supervision. All authors have approved the final version of the manuscript to be submitted.

Anthi Kelempisioti, Pasi J Eskola contributed equally to this work.

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Kelempisioti, A., Eskola, P.J., Okuloff, A. et al. Genetic susceptibility of intervertebral disc degeneration among young Finnish adults. BMC Med Genet 12, 153 (2011).

Download citation


  • Nucleus Pulposus
  • Single Nucleotide Polymorphism
  • Disc Degeneration
  • Lumbar Disc Herniation
  • Lumbar Disc Degeneration