Objectives A previous study identified a significant association between several single nucleotide polymorphisms (SNPs) and lumbar disc degeneration (LDD) in Indians. To validate the association between these SNPs and specific lumbar spine pathologies, we performed a case–control study in Chinese Han population.
Design An observational study.
Setting University Hospital in Nanning, China.
Participants This study included 428 patients with LDD and 400 normal controls.
Outcome measures Patients with LDD were classified into four subgroups, including disc herniation only (subgroup 1), discopathies or/and osteochondrosis associated with disc herniation (subgroup 2), spinal stenosis or/and spondylolisthesis (subgroup 3) and degenerative scoliosis (subgroup 4). This study was conducted by examining two aspects: environmental factors and SNP genotyping. The environmental factors were evaluated with a questionnaire survey including questions about body mass index, smoking habits, the physical demands of their job and exposure to vibrations. Rs1337185, rs5275, rs5277, rs7575934, rs3213718 and rs162509 were genotyped using a PCR-based invader assay.
Results The physical workload was significantly higher in patients with lumbar spine pathologies than in the normal controls (p=0.035). The genotype and allele frequencies of rs1337185 and rs162509 were significantly different between the patients with LDD and the normal controls. In rs1337185, a significant association was found between the C allele (risk allele) and the presence of disc herniation (OR=1.80; 95% CI 1.21 to 2.68; p=0.003, adjusted p=0.012) and the presence of spinal stenosis and spondylolisthesis (OR=1.92; 95% CI 1.29 to 2.89; p=0.001, adjusted p=0.004). In rs162509, the G allele represented 1.58-fold increased risk to suffer from disc herniation (OR=1.58; 95% CI 1.20 to 2.09; p=0.001, adjusted p=0.004).
Conclusion The SNPs rs1337185 in COL11A1 and rs162509 in ADAMTS5 are associated with susceptibility to LDD. The C allele of rs1337185 is risky for patients who are affected by lumbar pathologies such as disc herniation, stenosis and spondylolisthesis. The G allele of rs16250 represents a risk factor for the development of disc herniation.
- disc degeneration
- environmental factor
- single nucleotide polymorphism
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
Strengths and limitations of this study
This study was the first to validate the association between rs1337185, rs5275, rs5277, rs7575934, rs3213718 and rs162509 and specific lumbar spine pathologies.
This was also the first study to confirm an association between rs1337185 in COL11A1 and rs162509 in ADAMTS5 and lumbar disc degeneration among the Chinese Han population.
A limitation is that the sample size was relatively small.
Lumbar disc degeneration (LDD) is one of the most common musculoskeletal degenerative diseases, which may lead to lower back pain and unilateral or bilateral leg pain.1 Lower back pain affects 70%–85% of all people during their lifetime. Twenty per cent of patients with lower back pain caused by LDD require surgical treatment to relieve their prolonged or aggravated back or/and leg pain.2–4 Despite decades of research, the specific factors that contribute to LDD are still not fully understood. A variety of environmental factors, such as occupation, sporting activities, smoking and obesity, have been implicated in the aetiology and pathogenesis of LDD.5–10 Moreover, some studies have identified associations between genetic factors and LDD.11–13 This finding suggests that genetic factors may be important in determining the risk for the occurrence and development of LDD.
To date, more than 30 candidate genes have been analysed in association with LDD.14 15 The proteins encoded by these genes are classified into several categories based on their potential function in the disc: structural components, matrix turnover and organisation and inflammatory mediators.16 Recently, Rajasekaran et al 17 selected 58 single nucleotide polymorphism (SNP) from 35 candidate genes and performed a genetic association analysis between these SNPs and LDD in a cohort of young Indian adults. Their study identified significant associations between six SNPs in five genes (rs1337185 of COL11A1, rs5275 and rs5277 of COX2, rs7575934 of IL1F5, rs3213718 of CALM1 and rs162509 of ADAMTS5) and severe disc degeneration. According to the hypotheses of disease aetiology, the study of genetic associations is a powerful method for identifying the candidate genes or genomic regions that contribute to a specific disease. However, the lack of replication and false-positive findings are considered the main drawbacks of association studies.18 Thus, replicating the associations in different ethnic groups is crucial to confirm the results of association studies. Furthermore, Rajasekaran et al’s study was limited by the marginal p values, which were close to 0.05.17 To validate the association of rs1337185, rs5275, rs5277, rs7575934, rs3213718 and rs162509 with specific lumbar spine pathologies in an independent Chinese Han population, we conducted a genetic case–control study and performed a subgroup analysis of environmental factors in the development of LDD.
Materials and methods
Subjects and data collection
A total of 828 subjects were studied. Overall, 428 patients with lumbar spine pathologies (226 males and 202 females) were recruited from Spine Surgery, the First Affiliated Hospital of Guangxi Medical University between January 2012 and August 2016. The diagnosis of LDD required the following two criteria: (1) confirmed by MRI and (2) a history of lower back pain or/and leg pain for longer than 3 months. Detailed diagnoses of LDD and classification into one of four different mutually exclusive subgroups (designed subgroups 1–4) were performed by a senior spine specialist based on frontal and lateral X-ray images or CT scans and the sagittal and axial MRI images obtained using a 1.5 T MRI Achieva scanner (Philips Medical Systems, Best, The Netherlands). The classification of four subgroups was as follows: (1) subgroup 1 included 156 patients affected by disc herniation only; (2) subgroup 2 included 84 patients affected by discopathies or/and osteochondrosis associated with disc herniation; (3) subgroup 3 included 141 patients affected by spinal stenosis or/and spondylolisthesis and (4) subgroup 4 included 47 patients affected by degenerative scoliosis. The simultaneous presence of other orthopaedic diseases, such as osteoarthritis, osteoarthrosis and osteoporosis, was recorded (figure 1). Patients with an intraspinal tumour, trauma, inflammatory disease and rheumatoid arthritis were excluded from the study. Four hundred normal control subjects (224 males and 176 females) were enrolled from the Physical Examination Centre, the First Affiliated Hospital of Guangxi Medical University. In the control group, the lumbar spine MRI and CT scans were performed in 336 and 64 normal subjects, respectively. The clinical and radiographic examinations were performed by experienced orthopaedic surgeons to rule out any potential low back pain and LDD. The protocol was approved by the Ethics Review Committee of the First Affiliated Hospital of Guangxi Medical University, and signed informed consent was obtained from all subjects. All participants were Han Chinese who lived in and around the South China region.
Evaluation of environmental factors
According to the previous studies,9 the potential environmental risk factors of lumbar spine pathologies were evaluated based on the physical load assessment questionnaire and the occupational classification point standard.19 20 The collected information included body mass index (BMI), smoking habits, the physical workloads for the majority of the working years (evaluated by the following scoring system: score 0: ≤4 hours/day and ≤1 year; score 1: ≤4 hours/day and ≤2.5 years; score 2: ≤4 hours/day and 1–5 years or ≥4 hours/day and 2.5 years; score 3: ≤4 hours/day and >5 years or ≥4 hours/day and ≥2.5 years) and, over the past year, the number of hours per day spent driving or as a passenger in a motorised vehicle (exposure to vibrations).
Genomic DNA was extracted from peripheral blood leucocytes using genomic DNA isolation kits (Promega, Madison, Wisconsin, USA) according to the manufacturer’s instructions. The primers, probes and reaction conditions are available on request.
Genotyping was performed using a PCR-based invader assay with the probe sets designed. The genotyping results were obtained using an ABI 7900 sequence detection system (Applied Biosystems, Foster City, California, USA). Genotyping was done by laboratory personnel who were blinded to the status of the subjects. Further, 10% of the samples were tested twice to validate the genotyping results with 100% reproducibility. Two authors independently reviewed the genotyping results, data entry and statistical analyses.
A standard χ2-analysis was used to examine the differences in environmental risk factors, allelic frequencies and genotype distributions between the patients with LDD and controls using the SPSS software version 17.0. Considering the type I error caused by multiple testing, p values were adjusted by using Bonferroni correction (R Software 3.3.2 for Windows). The p value was adjusted to the corresponding original 5% level by dividing by the number of subgroups. Hardy-Weinberg equilibrium was tested using a goodness-of-fit χ2 test. The OR and 95% CI were calculated using the reported risk allele as a reference. Statistical significance was considered at p<0.05.
Influence of the environmental risk factors
The characteristics of the overall population of cases and controls, including age, gender, BMI, smoking habit, physical workload, exposure to vibrations and presence of other orthopaedic conditions, are shown in table 1. The patients had an average age of 49.35±15.67 years, and the control group had an average age of 45.92±13.21 years. Among the cases, there were more males (226/428, 52.8%) than females (202/428, 47.2%), whereas among controls, there were almost equal numbers of males (224/400, 56.0%) and females (176/400, 44.0%). As the study design, the simultaneous presence of other orthopaedic diseases (13.3% of cases) was presented only in the cohort of cases. Physical workload was significantly higher in patients with lumbar spine pathologies than in the normal controls (1.61±1.31 vs 1.12±1.03, p=0.035). However, neither groups showed a significant difference in BMI, smoking or exposure to vibrations (all p>0.05).
Associations between lumbar spine pathologies and putative conventional risk factors are also reported in table 1. Of the 428 patients with lumbar spine pathologies, 156 patients (36.4%) suffered from the disc herniation alone (subgroup 1), 141 patients (33.0%) suffered from spinal stenosis or/and spondylolisthesis (subgroup 3), 84 patients (19.6%) suffered from herniation associated with discopathies or/and osteochondrosis (subgroup 2) and 47 patients (11.0%) suffered from degenerative scoliosis (subgroup 4). A higher physical workload was exhibited in subgroup 3 (1.86±1.53) compared with the other three subgroups (all p<0.01). There was no statistical difference among subgroups 1, 2 and 4 (p>0.05). The four subgroups did not differ significantly regarding age, gender, BMI, smoking, exposure to vibrations or the presence of other orthopaedic conditions (all p>0.05).
Genetic association analysis
A total of 828 subjects (428 cases and 400 controls) were successfully genotyped and subjected to the statistical analysis. The distributions of alleles and genotypes for the six SNPs (rs1337185, rs5275, rs5277, rs7575934, rs3213718, rs162509) are presented in table 2. No significant deviation of genotype frequencies from the Hardy-Weinberg equilibrium was noted in the case and control groups (all p>0.05). The genotype and allele frequencies of rs1337185 and rs162509 were significantly different between the patients with lumbar spine pathologies and the normal controls.
Regarding rs1337185, a significant association was found between the C allele and the total patients. This indicated that the C allele might lead to a higher risk for lumbar spine pathologies in the Chinese population (OR=1.72; 95% CI 1.25 to 2.34; p=0.001). In accordance with the total cases, a significant association was found between the C allele and the presence of disc herniation only (subgroup 1) (OR=1.80; 95% CI 1.21 to 2.68; p=0.003). In the patients who were affected by spinal stenosis or/and spondylolisthesis (subgroup 3), the C allele was a risk factor (OR=1.92; 95% CI 1.29 to 2.89; p=0.001). After Bonferroni correction, rs1337185 showed statistical significance in abovementioned comparisons (subgroup 1 vs control, adjusted p=0.012; subgroup 3 vs control, adjusted p=0.004) (table 2).
Regarding rs162509, the risk allele (G allele) frequency of the total cases was highly significantly different from that in the controls (OR=1.38; 95% CI 1.13 to 1.69; p=0.001). Consistently, the allele distribution showed a higher frequency of the G allele in the patients with disc herniation only (subgroup 1) (OR=1.58; 95% CI 1.20 to 2.09; p=0.001, adjusted p=0.004). In the patients with degenerative scoliosis, a 1.72-fold association was found for the G allele in genetic analysis (OR=1.72; 95% CI 1.09 to 2.78; p=0.019, adjusted p=0.076). After analysis and Bonferroni correction, it revealed that rs162509 is associated with disc herniation but not with degenerative scoliosis phenotype (table 2).
Neither the genotype nor the allele frequencies of rs5275, rs5277, rs7575934 and rs3213718 were significantly different between the cases and the normal controls. Furthermore, no other significant findings in the aforementioned four SNPs were observed even when grouping subgroups 1, 2, 3 and 4.
It is generally recognised that LDD is a complex, multifactorial disease that is determined by genetic and environmental interactions.21 In this study, the patients and controls showed significant differences in the physical workload for the majority of the working years, which is consistent with previous observations.7 9 22 However, the two groups in our study did not show significant differences in other environmental factors, including age, gender, BMI, smoking and exposure to vibrations. Recently, genetic factors have been suggested to play an important role in the development and progression of disc degeneration.23 24 Two pivotal approaches have been used to map candidate genes and genome regions: linkage analysis and association study.25 Compared with linkage analysis, an association study is a more robust approach for identifying predisposition genes for common diseases and complex traits. More than 30 candidate genes related to LDD have been identified by genetic association studies. However, the most common weakness is difficulty in replicating the previous association signals.
Recently, Rajasekaran et al 17 recruited 308 Indian subjects with mild LDD and 387 Indian subjects with severe LDD and performed a genetic association analysis for 58 SNPs of 35 candidate genes. The study detected an association between severe LDD and six SNPs in five genes including rs1337185 of COL11A1, rs5275 and rs5277 of COX2, rs7575934 of IL1F5, rs3213718 of CALM1 and rs162509 of ADAMTS5. These five genes were important candidate genes related to LDD, which could be grouped into three categories based on their potential functional zones: structural genes (COL11A1, CALM1), degradative genes (ADAMTS5) and inflammatory genes (IL1F5, COX2). To the best of our knowledge, no study has validated the associations or conducted a subgroup analysis. In this study, we identified COL11A1 (rs1337185) and ADAMTS5 (rs162509) as gene polymorphisms associated with the susceptibility of patients to LDD in a Chinese Han cohort. Furthermore, we analysed four mutually exclusive subgroups of all cases having a condition in common such as a hernia, discopathy, spinal stenosis or spondylolisthesis and degenerative scoliosis. Regarding rs1337185, the C allele was associated with a higher risk for patients suffering from disc herniation only (subgroup 1) and spinal stenosis or/and spondylolisthesis (subgroup 3). Regarding rs162509, a higher frequency of the G allele was presented in the patients with disc herniation only (subgroup 1) and degenerative scoliosis (subgroup 4). After Bonferroni correction, both SNPs (rs1337185 and rs162509) showed significantly different allele and genotype frequencies between subgroup 1 and control group and rs1337185 showed statistically different allele and genotype frequencies between subgroup 3 and control group. Unfortunately, no detailed information was reported about the risk or protective alleles at rs1337185 and rs162509 in the Indian study. Thus, we did not evaluate the difference in the genetic associations between the Han Chinese population and Indian population.17
The COL11A1 gene is highly expressed in intervertebral discs and encodes the α1 chain of type XI collagen, which suggests that it is critical for intervertebral disc metabolism. Mio and colleagues26 identified a strong association between rs1676486 in the COL11A1 gene and lumbar disc herniation in Japanese patients. This study also found that COL11A1 mRNA expression in patients with lumbar disc herniation was decreased along with an increase in the severity of degeneration. This finding indicated that the functional polymorphism might produce unstable COL11A1 transcripts. Moreover, another important polymorphism in COL11A1 gene (rs1337185) was found to have significant associations with MRI-defined lumbar disc bulging in a Finnish cohort.23 In agreement with these findings, our study confirmed that a significant association exists between rs1337185 and patients with LDD. Based on the subgroup analysis, the C allele was identified as a risk allele in the patients affected by disc herniation only (subgroup 1) and spinal stenosis or/and spondylolisthesis (subgroup 3). As a member of the large family of metalloproteases, ADAMTS5 has recently received increased attention due to its role in balancing the synthesis and degradation of the extracellular disc matrix.27 28 The mRNA and protein expression of ADAMTS5 has been measured in intervertebral disc tissues.29 Furthermore, the increased expression of ADAMTS5 was correlated with degenerative changes in patients with chronic lower back pain and lumbar disc herniation.30 31 Correspondingly, Wu et al 32 reported that the genetic polymorphisms in ADAMTS5 (rs151058, rs229052 and rs162502) might be associated with susceptibility to LDD. In accordance with those previous reports, the current study identified a significant association between rs162509 in ADAMTS5 gene and LDD in Chinese Han patients. The G allele represented an approximately 1.58-fold increased risk factor for developing disc herniation. It remains plausible that COL11A1 and ADAMTS5 genes may be involved in the aetiology of LDD. However, the precise role of rs1337185 and rs162509 is still unknown. Functional analysis of the COL11A1 and ADAMTS5 genes might help elucidate the real genetic effect on the aetiopathogenesis of LDD.
The previously reported association between rs5275, rs5277, rs7575934 and rs3213718 and disc degeneration disease in young Indian adults was not replicated in our study. Several factors could have led to this lack of replication. First, the discrepancy may be related to genetic differences between the Chinese and Indian populations. Second, the previous study was limited by its marginal p value, which might have led to false positive findings. Third, a difference in phenotype selection may have been a crucial factor that could account for this conflicting result. Genetic associations may change when different phenotypes of LDD are studied.33 There are obvious discrepancies in the clinical phenotype of patients with LDD, such as annular tears, hyperintense zones, signal intensity, disc height and disc herniation. In the current study, we selected a highly specific phenotype of LDD, which was the simultaneous presence of disc degeneration on the MRI images and leg or low back pain in a population of symptomatic patients. However, the severity of LDD in the previous study was evaluated on the basis of signal changes in the lumbar discs and was graded using the total disc degeneration score. This discrepancy corresponds well with the fact that the disease course of LDD is complicated and heterogeneous.34 Therefore, we postulate that different phenotypes in the disease course of LDD may be affected by different genetic polymorphisms.
There were several limitations in our study. First, the sample size was relatively small. Moreover, stratification by pathological subgroups reduced the number of subjects compared, which may have weakened the statistical power. We would like to enlarge the group of cases in the future, to confirm our results in a larger cohort of subjects. Second, if the original data were available, we could perform a meta-analysis of other relevant cohorts. More convincing results may be provided by a large-scale meta-analysis.
The current study showed an association between rs1337185 in COL11A1 and rs162509 in ADAMTS5 and LDD predisposition in a Chinese Han population. The C allele of rs1337185 was a risk factor for patients who are affected by lumbar pathologies such as disc herniation, stenosis and spondylolisthesis. The G allele of rs16250 represented a risk factor for the development of disc herniation.
The authors thank all the participants in this study.
Review history and Supplementary material
Contributors HJ carried out the molecular genetic studies, participated in its design and coordination and drafted the manuscript. JJ participated in the molecular genetic studies and the sequence alignment. QY and XZ performed the statistical analysis and participated in the design of the study. ZX conceived of the study and participated in its design and coordination and helped to draft the manuscript.
Funding This work was supported by the Natural Science Foundation of China (81460353), Guangxi Natural Science Foundation (2015GXNSFBA139167) and Youth Science Foundation of Guangxi Medical University (GXMUYSF201329).
Competing interests None declared.
Patient consent Obtained.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement Extra data can be accessed via the Dryad data repository at http://datadryad.org/with the doi:10.5061/dryad.5cd20.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.