Is intrinsic lumbar spine shape associated with lumbar disc degeneration? An exploratory study
BMC Musculoskeletal Disorders volume 21, Article number: 433 (2020)
Lumbar disc degeneration (LDD) is a condition associated with recurrent low back pain (LBP). Knowledge regarding effective management is limited. As a step towards the identification of risk, prognostic or potentially modifiable factors in LDD patients, the aim of this study was to explore the hypothesis that intrinsic lumbar spine shape is associated with LDD and clinical outcomes in symptomatic adults.
3 T MRI was used to acquire T2-weighted sagittal images (L1-S1) from 70 healthy controls and LDD patients (mean age 49 years, SD 11, range 31–71 years). Statistical Shape Modelling (SSM) was used to describe lumbar spine shape. SSM identified variations in lumbar shape as ‘modes’ of variation and quantified deviation from the mean. Intrinsic shape differences were determined between LDD groups using analysis of variance with post-hoc comparisons. The relationship between intrinsic shape and self-reported function, mental health and quality of life were also examined.
The first 7 modes of variation explained 91% of variance in lumbar shape. Higher LDD sum scores correlated with a larger lumbar lordosis (Mode 1 (55% variance), P = 0.02), even lumbar curve distribution (Mode 2 (12% variance), P = 0.05), larger anterior-posterior (A-P) vertebral diameter (Mode 3 (10% variance), P = 0.007) and smaller L4-S1 disc spaces (Mode 7 (2% variance), P ≤ 0.001). In the presence of recurrent LBP, LDD was associated with a larger A-P vertebral diameter (Mode 3) and a more even lumbar curvature with smaller L5/S1 disc spaces (Mode 4), which was significantly associated with patient quality of life (P = 0.002–0.04, rp = 0.43–0.61)).
This exploratory study provides new evidence that intrinsic shape phenotypes are associated with LDD and quality of life in patients. Longitudinal studies are required to establish the potential role of these risk or prognostic shape phenotypes.
Lumbar disc degeneration (LDD) is a condition associated with recurrent low back pain (LBP) , the lifetime prevalence of which may be as much as 80% . LDD is commonly evaluated using the Pfirrmann  and modified Pfirrmann grading systems , which use T2-weighted sagittal MR images of the intervertebral discs to grade the lack of distinction between the nucleus and annulus and the reduction in intervertebral disc height and signal intensity associated with LDD  .
Current treatment approaches for recurrent LBP offer small to moderate effects in terms of a sustained improvement in quality of life and disability [5, 6]. As a step towards a potential role of lumbar shape and Statistical Shape Modelling (SSM) in the identification of risk, prognostic or potentially modifiable factors in LDD patients, it is pressing to identify the phenotypes associated with LDD through the examination of observable traits. Since the external spinal curvature is not always representative of the internal spinal geometry , it seems appropriate to examine the internal architecture of the lumbar spine or intrinsic lumbar shape.
SSM is a statistical image analysis technique that, through the reduction of variables using principal components analysis, is used to describe and quantify variations in joint morphology and intrinsic shape [7, 8]. Conventionally, geometric measurements of lumbar lordosis are performed but these do not account for the distribution of lumbar curvature or morphological variability in the degenerate spine . However, SSM can describe this and has been shown to be a more reliable and precise method of characterising lumbar curvature (4% measurement error) when compared with conventional measurements, such as the Cobb angle (10% error) .
To our knowledge, SSM has not been used to examine intrinsic lumbar shape in LDD adults.
The aim of this exploratory study is to create a shape model of the lumbar spine using SSM and to test the hypothesis that intrinsic lumbar spine shape is associated with LDD and clinical outcomes in adults.
Seventy participants were recruited through advertisement from primary and secondary care between September 2015 and May 2017. Each participant provided informed consent and met strict inclusion and exclusion criteria (Table 1).
A 3 T Verio MRI scanner (Siemens Medical Systems, Erlangen, Germany) was used to acquire T2-weighted sagittal lumbar spine images (L1-L5/S1) from the asymptomatic volunteers and symptomatic participants as part of routine NHS care. The majority of scans were acquired at the same time of day in supine crook lying supported by a wedge (with shoulders, thorax and pelvis level) following a 10 min rest period .
T2-weighted sagittal images were viewed using the Picture Archiving and Communications System (PACs) (Synpase, Fujifilm Medical Systems, Tokyo, Japan). Pfirrmann and Modified Pfirrmann disc grades were determined by an experienced consultant radiologist, blinded to the demographics and clinical data of participants. LDD sum scores were calculated through the summation of Pfirrmann scores from each lumbar disc . Four groups were identified based upon the presence or absence of pain and LDD (modified Pfirrmann grade ≥ 6). Pain was determined by the presence or absence of recurrent LBP for greater than 3 months. Self-reported clinical outcome tools including the Short Form 36, Version 2 (SF-36), the Oswestry Disability Index (ODI) and the Hospital Anxiety and Depression Scale (HADS) were used to assess quality of life, function and anxiety and depression, respectively.
Statistical shape modelling (SSM)
SSM is an image processing method designed to characterise the shape of an object within a series of image . It uses principal components analysis to identify different ‘modes’ of variation (principal components) in sagittal spinal shape and computes mode scores (output variables) for each image to quantify that variation by calculating deviation from the mean of all the images [10, 11, 13].
Seventy consistent mid-sagittal slices (maximum spinal canal width) were selected from the participant’s data using Image J software and exported as Bitmap image files. A previously tested 168-point lumbar spine template [10, 14] which defined the outline of each vertebral body from the 1st lumbar to 1st sacral vertebrae (L1-S1) was recreated in custom SSM software (Shape software, University of Aberdeen) and images uploaded for analysis. In each image, key points were identified manually by the same tester and the point template fitted to the spine (Fig. 1). Intra-class correlation coefficients (ICC) demonstrated excellent reliability of key point positioning of both x (ICC = 0.99) and y (ICC = 0.98) co-ordinates. The influence of positioning and size differences between subjects were removed by scaling, translation and rotation (Procrustes transformation) of the point coordinates. The constructs of a previous model, built from sagittal MRI data (L1-S1) of 30 asymptomatic volunteers aged 20–52 (mean 29, standard deviation (SD) 9.6) years , were used to enable consistency in the shape variations described (shape ‘modes’). The original mode score distributions were normalized to have a mean of zero and unit standard deviation so that the unit of measurement was in standard deviations. Thus, the mode scores for individuals in the current cohort were calculated in relation to those of a healthy asymptomatic cohort.
Statistical analysis was undertaken using SPSS software (Version 24, IBM SPSS statistics, IBM Corp.). The normality of the data was assessed using frequency histograms, quantile-quantile plots and the Shapiro-Wilk test. To assess preliminary trends in the data, associations between LDD (LDD sum scores), age, BMI and intrinsic shape (mode scores) were established using Pearson correlations. Levene’s Test confirmed homogeneity of variance. One-way repeated-measures analysis of variance (ANOVA) was used to assess the differences between the LDD groups with Games Howell (for samples with unequal variance) or Hochberg’s T2 (for samples with equal variance but unequal sample sizes) post-hoc comparisons. In the case of unequal variance Welch’s corrected F-ratio was reported. The observed power and effect sizes were also computed (r = √ η2, where r = effect size and η2= eta squared or SSM (the between group effects)/ SST (the total amount of variance in the data)). The effects of potential confounding variables such as age, BMI and sex were explored using Pearson’s r correlation coefficient (rp) and independent sample t-tests. Results were considered significant at P < 0.05 for all tests. Missing data were excluded case-wise from the analysis and were not replaced by imputed values.
Seventy participants completed this study (31 male, 39 female, mean age 49 years (SD 11, range 31–71 years), BMI mean 26 kg m− 2 (SD 5). Ninety one percent of the variation in participant lumbar shape was explained by the first seven modes (M1 – M7) in descending order of variance (Fig. 2) . Modes contributing to greater than 1.5% of the variance observed (≤M8) were chosen according to the scree plot.
LDD sum scores significantly correlated with age (rp = 0.51, P ≤ 0.0001). There was no difference in LDD sum scores between sexes (P = 0.60). LDD sum scores did not correlate with BMI (rp = − 0.04, P = 0.76). Higher LDD sum scores were significantly associated with a more curvy or lordotic spine (corresponding to negative M1 scores, rp = − 0.28, P = 0.02), a more evenly distributed lumbar curve (negative M2 scores, rp = − 0.24, P = 0.05), larger anterior-posterior (A-P) vertebral diameter relative to vertebral height (positive M3 scores rp = 0.32, P = 0.007) and smaller L4-S1 intervertebral disc spaces (negative M7 scores, rp = − 0.56, P ≤ 0.001) (Fig. 3).
Four groups were identified based upon the presence or absence of LDD and recurrent LBP; ‘LDD pain’ (n = 24), ‘LDD no pain’ (n = 24), ‘No LDD no pain’ (n = 19) and ‘No LDD pain’ (n = 3) (Table 2). Although there was no significant difference in age between groups (F = 2.51, P = 0.07), there was a significant difference in BMI (F = 3.24, P = 0.03).
Analysis of variance revealed a significant difference in the M3 phenotype (F = 3.892, P = 0.05 ,η2 =0.13, α = 0.05, observed power = 0.72) between groups. M3 corresponded to changes in A-P vertebral diameter relative to vertebral height. Post-hoc tests showed that ‘LDD no pain’ (95% CI = -2.61 to − 0.08; P = 0.03) and ‘LDD pain’ groups (95% CI = -3.77 to − 0.27; P = 0.02) had significantly larger A-P vertebral diameters relative to vertebral height than ‘no LDD no pain’ groups (Table 3, Figs. 4 & 5). This suggested an association between M3 and LDD and/or recurrent LBP. Males had significantly higher M3 scores (P = 0.001) indicating a larger A-P vertebral diameter. Higher BMI was also associated with larger A-P vertebral diameters (positive M3 scores, rp = 0.26, P = 0.03) but not age (rp = 0.17, P = 0.20). These results supported an association between the M3 differences observed between LDD groups and BMI and sex.
Analysis of variance between groups also uncovered a significant difference in M4 (F = 3.967, P = 0.01, η2 =0.15, α = 0.05, observed power = 0.81) that was not described through direct correlation with LDD alone. M4 described variation in the evenness of the lumbar curvature occurring together with changes in L5-S1 intervertebral disc space. Differences in M4 were not explained by sex (P = 0.20), BMI (rp = 0.20, P = 0.20) or aging (rp = 0.30, P = 0.05). Post-hoc comparisons revealed that the significant differences lay between ‘LDD pain’ and ‘no LDD no pain’ groups (95% CI = -4.96 to − 0.57; P = 0.007) (Table 3, Figs. 4, 5 and 6); the ‘LDD pain’ group had a more even lumbar curvature with smaller and less wedged L5/S1 disc spaces (positive M4). This indicated that in the presence of recurrent LBP, LDD was associated with a significant difference in the M4 phenotype.
The M4 phenotype, identified as significantly different between LDD groups, was also significantly associated with quality of life in LDD patients, specifically bodily pain, vitality, social function and mental health (SF-36) (Table 4). Intrinsic shape did not correlate with measures of self-reported function (ODI), depression or anxiety in patients with recurrent LBP.
SSM has been used in recent studies to successfully detect healthy morphological changes associated with aging, functional activity, osteoporosis and osteoarthritis [15,16,17,18]. However, this is the first time that SSM has been used to explore associations between intrinsic spine shape variations and LDD. The results suggest that intrinsic lumbar shape is associated with LDD and quality of life in adults. Therefore, it seems that SSM could be used in future longitudinal research to explore factors associated with LDD modifiers, risk or prognosis.
In this current study, LDD was associated with a more lordotic or ‘curvy’ lumbar spine (M1), suggesting a direct correlation between LDD and lumbar lordosis. This association has been previously documented and has been attributed to a decrease in multifidus density . However, it is of interest in this current study that LDD, in the absence of frank neural compression, correlated with ‘curviness’ (increased lumbar lordosis), since lumbar disc herniation (disc prolapse with neural compression) is known to be characterised by a flattened lumbar lordosis . A possible explanation for this finding is that herniation with concomitant neural compression results in a protective response during which the lumbar spine flattens. The arch-model of the spine predicts that a flatter spine is stiffer due to a greater follower load and this would agree with the lower flexibility experienced by patients with advanced degenerative conditions . It is unclear, however, whether the lumbar lordosis increases as LDD progresses or whether a larger lordosis is a risk factor for LDD. Further, longitudinal studies will be necessary to explore this.
LDD also correlated with a more evenly distributed lumbar curve (M2), increased A-P relative diameter of lumbar vertebrae (M3) and smaller L4-S1 intervertebral spaces (M7). Since LDD has been previously associated with disc dehydration and reduced disc height, findings of reduced distal intervertebral spaces seem representative of the condition under investigation [22, 23]. Overall, given the sparsity of high quality, longitudinal studies, it remains challenging to determine true causality.
Analysis of variance between LDD groups established an association between LDD and larger A-P relative diameters. However, M3 was also found to be significantly associated with sex (males had significantly larger vertebral depths) and BMI. This finding is unsurprising in light of the known significant association between degeneration and genetics  and previous SSM research which has found positive correlations between M3, BMI and sex . Furthermore, since an increased risk of spondylolysis has been recently associated with smaller A-P cross sectional areas , it is possible that increasing A-P diameters may also represent a risk factor for LDD and reduced lumbar mobility.
However, the M4 phenotype, which corresponded to a greater evenness of the lumbar curvature and the smaller, less wedge-shaped, L5-S1 intervertebral spaces, was significantly associated with LDD in patients with recurrent LBP. As a preliminary step towards determining the potential for patient impact, this exploratory study also found new evidence to suggest an association between intrinsic shape phenotypes and self-reported clinical outcomes; the M4 phenotype also correlated significantly with features that affect patient quality of life, specifically bodily pain, social functioning and mental health. This finding is in agreement with previous studies in adults with degenerative changes which also found significant associations between posture and quality of life using the SF-36 [25, 26]. Therefore, the M4 phenotype warrants further investigation as a biomarker of LDD, particularly given the association with patient quality of life. In future such phenotypes could be used to stratify treatment for patients with degenerative conditions, thus ensuring that the right patient receives the right care.
In light of the established association between Modic change and low back pain , a preliminary analysis of intrinsic shape phenotypes was also undertaken using the Modic subclassification system in this cohort. Participants were assigned to each group based on the presence or absence of Modic changes and/ or pain. SSM analysis revealed that both the M3 (p = 0.03) and M4 (p = 0.03) phenotypes were similarly significant irrespective of the method of subclassification used, reinforcing the known association between LDD and Modic changes .
The assessment of supine MRI scans is a limitation in this current study. However, a standardised supine position was adopted so that routine NHS MRI scans could be directly compared with those of healthy controls. Previous SSM research has shown a correlation between intrinsic lumbar shape in supine, sitting and standing positions . This means that an element of intrinsic spine shape, or an individual’s ‘spinal signature’, is conserved despite the position adopted . It is also important to note that this signature affects the response to loading and the natural style adopted for lifting. The most lordotic spines become more lordotic under load whereas straighter spines or those with a moderate lordosis become straighter . When lifting a weight from the floor, those with curvier spines prefer to flex whereas those with straighter spines prefer to squat and find flexing difficult . Therefore, one can assume that any shape differences identified in this supine study would be further enhanced in standing and may affect the biomechanics of the lumbar spine.
It is acknowledged that the results of this exploratory, cross-sectional work are not generalisable and do not establish causality and as such should be interpreted judiciously. Firstly, the strategic recruitment of a specific number of participants to each group was not possible since the assignment to each group was dependent on the MRI outcome and strict inclusion criteria. Secondly, the results from the model used in this current study cannot be directly compared with models constructed from different images. It is possible that larger mode scores (SDs) reflected differences between the healthy adult model (mean 29 years) and the older adult cohort used in this study; LDD was found to be significantly correlated with age.
This exploratory SSM study provides new evidence that there is an association between intrinsic lumbar shape, LDD and quality of life in patients. This highlights the potential role of SSM and intrinsic shape in the identification of risk, prognostic or potentially modifiable factors in LDD patients with recurrent pain. Further longitudinal research will be required to characterise the spine over time and establish true causality.
Availability of data and materials
The datasets for this current study are available from the corresponding author upon reasonable request.
Analysis of Variance
Body Mass Index
Hospital Anxiety and Depression Scale
Intra-class Correlation Coefficients
Low Back Pain
Lumbar Disc Degeneration
Magnetic Resonance Imaging
National Health Service
Oswestry Disability Index
Statistical Shape Modelling
Teraguchi M, Yoshimura N, Hashizume H, et al. Prevalence and distribution of intervertebral disc degeneration over the entire spine in a population-based cohort: the Wakayama spine study. Osteoarthr Cartil. 2014;22(104–10):2014.
Livshits G, Popham M, Malkin I, et al. Lumbar disc degeneration and genetic factors are the main risk factors for low back pain in women: the UK twin spine study. Ann Rheum Dis. 2011;70(10):1740–5. https://doi.org/10.1136/ard.2010.137836.
Pfirrmann CW, Metzdorf A, Zanetti M, et al. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine. 2001;26(17):1873–8 published online first: 2001/09/25.
Griffith JF, Wang YX, Antonio GE, et al. Modified Pfirrmann grading system for lumbar intervertebral disc degeneration. Spine. 2007;32(24):E708–12. https://doi.org/10.1097/BRS.0b013e31815a59a0 published online first: 2007/11/17.
Gurung T, Ellard DR, Mistry D, et al. Identifying potential moderators for response to treatment in low back pain: a systematic review. Physiotherapy. 2015;101(3):243–51. https://doi.org/10.1016/j.physio.2015.01.006.
Keller A, Hayden J, Bombardier C, et al. Effect sizes of non-surgical treatments of non specific low-back pain. E Spine J. 2007;16(11):1776–88. https://doi.org/10.1007/s00586-007-0379-x.
Meakin JR, Smith FW, Gilbert FJ, et al. The effect of axial load on the sagittal plane curvature of the upright human spine in vivo. J Biomechanics. 2008;41(13):2850–4. https://doi.org/10.1016/j.jbiomech.2008.06.035 published online first: 2008/08/22.
Cootes TF, Taylor CJ.  Anatomical statistical models and their role in feature extraction. Br J Radiol;77 Spec No 2:S133-S139. doi: https://doi.org/10.1259/bjr/20343922 [published online first: 2005/01/29].
Ali AH, Cowan AB, Gregory JS, et al. The accuracy of active shape modelling and endplate measurements for characterising the shape of the lumbar spine in the sagittal plane. Comp Methods Biomechanics Biomed Engineering. 2012;15(2):167–72. https://doi.org/10.1080/10255842.2010.518962 published online first: 2012/01/25.
Meakin JR, Gregory JS, Smith FW, et al. Characterizing the shape of the lumbar spine using an active shape model: reliability and precision of the method. Spine. 2008;33(7):807–13. https://doi.org/10.1097/BRS.0b013e31816949e6.
Meakin JR, Gregory JS, Aspden RM, et al. The intrinsic shape of the human lumbar spine in the supine, standing and sitting postures: characterization using an active shape model. J Anat. 2009;215(2):206–11. https://doi.org/10.1111/j.1469-7580.2009.01102.x.
Takatalo J, Karppinen J, Taimela S, et al. Association of Abdominal Obesity with lumbar disc degeneration – a magnetic resonance imaging study. PloS one. 2013;8(2):e56244. https://doi.org/10.1371/journal.pone.0056244.
Cootes TF, Hill A, Taylor CJ, et al. Use of active shape models for locating structures in medical images. Image Vis Comput. 1994;12(6):355–65. https://doi.org/10.1016/0262-8856(94)90060-4.
Pavlova AV, Meakin JR, Cooper K, et al. The lumbar spine has an intrinsic shape specific to each individual that remains a characteristic throughout flexion and extension. Eur Spine J. 2014;23(Suppl 1):S26–32. https://doi.org/10.1007/s00586-013-3162-1 published online first: 2014/01/15.
Pavlova AV, Saunders FR, Muthuri SG, et al. Statistical shape modelling of hip and lumbar spine morphology and their relationship in the MRC National Survey of health and development. J Anat. 2017;231(2):248–59. https://doi.org/10.1111/joa.12631.
Pavlova AV, Meakin JR, Cooper K, et al. Variation in lifting kinematics related to individual intrinsic lumbar curvature: an investigation in healthy adults. BMJ Open Sport Exerc Med. 4(1):e000374. https://doi.org/10.1136/bmjsem-2018-000374 published Online First:2018/07/31.
Castro-Mateos I, Pozo JM, Cootes TF, et al. Statistical shape and appearance models in osteoporosis. Curr Osteoporos Rep. 2014;12(2):163–73. https://doi.org/10.1007/s11914-014-0206-3 published Online First: 2014/04/03.
Faber BG, Baird D, Gregson CL, et al. DXA-derived hip shape is related to osteoarthritis:findings from in the MrOS cohort. Osteoarthritis and cartilage / OARS, osteoarthritis research society. 2017;25(12):2031–8. https://doi.org/10.1016/j.joca.2017.09.006 published Online First:2017/09/25.
Kalichman L, Li L, Hunter DJ, et al. Association between computed tomography-evaluated lumbar lordosis and features of spinal degeneration, evaluated in supine position. Spine J. 2011;11(4):308–15. https://doi.org/10.1016/j.spinee.2011.02.010 published Online First: 2011/04/09.
Endo K, Suzuki H, Tanaka H, et al. Sagittal spinal alignment in patients with lumbar disc herniation. Eur Spine J. 2010;19(3):435–8. https://doi.org/10.1007/s00586-009-1240-1 published Online First: 2010/01/22.
Aspden RM. The spine as an arch. A new mathematical model. Spine. 1989;14(3):266–74 published Online First: 1989/03/01.
Adams MA, Bogduk N, Burton K, et al.  The Biomechanics of Back Pain: Churchill Livingstone Elsevier.
Barrey C, Jund J, Noseda O, et al. Sagittal balance of the pelvis-spine complex and lumbar degenerative diseases. A comparative study about 85 cases. E Spine J. 2007;16(9):1459–67. https://doi.org/10.1007/s00586-006-0294-6.
Wren TAL, Ponrartana S, Aggabao PC, et al. Increased lumbar Lordosis and smaller vertebral cross-sectional area are associated with Spondylolysis. Spine. 2018;43(12):833–8. https://doi.org/10.1097/brs.0000000000002480 published Online First: 2017/11/03.
Schwab F, Dubey A, Pagala M, Gamez L, Farcy JP. Adult scoliosis: a health assessment analysis by SF-36. Spine. 2003;28(6):602–6.
Ploumis A, Liu H, Mehbod AA, Transfeldt EE, Winter RB. A correlation of radiographic and functional measurements in adult degnerative scoliosis. SPine. 2009;34(15):1581–4. https://doi.org/10.1097/BRS.0b013e31819c94cc.
Modic MT, Ross JS. Lumbar Degenerative Disk Disease. Radiology. 2007;245(1):43–61.
Yu LP, Qian WW, Yin GY, Ren YX, Hu ZY. MRI assessment of lumbar intervertebral disc degeneration with lumbar degenerative disease using the Pfirrmann grading systems. PLoS One. 2012;7(12):e48074.
Thank you to Lesley Honeyfield (Lead Research Radiographer, Charing Cross Hospital, Imperial College Healthcare NHS Trust, London) for assisting with imaging optimisation and the acquisition of MR images and to all participants in this study. We also thank the University of Aberdeen Data Management Team for programming support for ‘Shape’ software.
JD would like to acknowledge the support of Versus Arthritis (Clinical Doctoral Fellowship award, grant number 20172) and the UKSSB and Society for Back Pain Research (Travel Fellowship award) for funding this research. Funders did not have a role in the design of the study, analysis or interpretation of the data.
Ethics approval and consent to participate
Ethical approval was obtained from the NHS Health Research Authority (NRES Committee London, Stanmore, REC reference number: 13/LO/0793). All participants in this study provided informed written consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Deane, J.A., Pavlova, A.V., Lim, A.K.P. et al. Is intrinsic lumbar spine shape associated with lumbar disc degeneration? An exploratory study. BMC Musculoskelet Disord 21, 433 (2020). https://doi.org/10.1186/s12891-020-03346-7