- Research article
- Open Access
- Open Peer Review
Mild (not severe) disc degeneration is implicated in the progression of bilateral L5 spondylolysis to spondylolisthesis
BMC Musculoskeletal Disordersvolume 19, Article number: 98 (2018)
Spondylolytic (or lytic) spondylolisthesis is often associated with disc degeneration at the index-level; however, it is not clear if disc degeneration is the cause or the consequence of lytic spondylolisthesis. The main objective of this computed tomography based finite element modelling study was to examine the role of different grades of disc degeneration in the progression of a bilateral L5-lytic defect to spondylolisthesis.
High-resolution computed tomography data of the lumbosacral spine from an anonymised healthy male subject (26 years old) were segmented to build a 3D-computational model of an INTACT L1-S1 spine. The INTACT model was manipulated to generate four more models representing a bilateral L5-lytic defect and the following states of the L5-S1 disc: nil degeneration (NOR LYTIC), mild degeneration (M-DEG LYTIC), mild degeneration with 50% disc height collapse (M-DEG-COL LYTIC), and severe degeneration with 50% disc height collapse(S-COL LYTIC). The models were imported into a finite element modelling software for pre-processing, running nonlinear-static solves, and post-processing of the results.
Compared with the baseline INTACT model, M-DEG LYTIC model experienced the greatest increase in kinematics (Fx range of motion: 73% ↑, Fx intervertebral translation: 53%↑), shear stresses in the annulus (Fx anteroposterior: 163%↑, Fx posteroanterior: 31%↑), and strain in the iliolumbar ligament (Fx: 90%↑). The S-COL LYTIC model experienced a decrease in mobility (Fx range of motion: 48%↓, Fx intervertebral translation: 69%↓) and an increase in normal stresses in the annulus (Fx Tensile: 170%↑; Fx Compressive: 397%↑). No significant difference in results was noted between M-DEG-COL LYTIC and S-COL LYTIC models.
In the presence of a bilateral L5 spondylolytic defect, a mildly degenerate index-level disc experienced greater intervertebral motions and shear stresses compared with a severely degenerate index-level disc in flexion and extension bending motions. Disc height collapse, with or without degenerative changes in the stiffness properties of the disc, is one of the plausible re-stabilisation mechanisms available to the L5-S1 motion segment to mitigate increased intervertebral motions and shear stresses due to a bilateral L5 lytic defect.
Isthmic spondylolysis (or lytic defect) is characterised by a bony defect of the pars interarticularis and occurs most commonly as a bilateral defect in the L5 vertebra [1, 2]. The mechanism for the progression of the lytic defect to spondylolisthesis remains unclear despite viable qualitative theories supported with clinical radiographic evidence [3,4,5,6,7,8,9].
Abnormal spinopelvic morphology and orientation measured through parameters such as pelvic incidence, sacral slope, pelvic tilt may create a biomechanical environment leading to the development of a pars defect and its progression to spondylolisthesis [3,4,5,6, 9]. In children and adolescents, the progression of the defect to spondylolisthesis is attributed to growth plate injury which may perpetuate to epiphyseal ring separation [8, 10]. Slippage in skeletally immature spines is not attributed to any disc degeneration, but rather growth plate lesions .
Lytic spondylolisthesis is often associated with disc degeneration at the olisthetic segment; however, it is not clear whether disc degeneration is the cause or consequence of lytic spondylolisthesis . Clinical evidence suggests that disc degeneration occurs more rapidly and prematurely in the presence of lytic spondylolisthesis [11, 12]. Szypryt et al. (1988) compared disc degeneration at the olisthetic segment in 40 consecutive patients with an age-and-sex matched control group of 40 asymptomatic volunteers. The authors observed that over the age of 25 years, the prevalence of disc degeneration was significantly higher in the spondylolytic group compared with the control group . Seitsalo et al. (1997) in their long-term follow-up study of operatively and conservatively treated isthmic spondylolisthesis patients reported a strong correlation between the severity of disc degeneration and the degree of slippage in the conservatively treated cohort. Progression to spondylolisthesis in adults is almost always attributed to disc degeneration [11, 13]. Case studies, however, have cited the rare progression of pars defect to vertebral slippage in adults without any associated disc degeneration, adding to the ambiguity around the role of disc degeneration in the pathomechanism for slippage progression [14, 15].
The intervertebral disc plays an important role in resisting shear and torsional forces acting on the spinal column during physiological bending motions. Biomechanical studies on cadaveric lumbar spines with healthy intervertebral discs have shown increased kinematics at the index-level following a bilateral lytic defect [16,17,18]. The loss of posterior tension band in the lumbar spine following the defect may redirect excessive shear and torsional forces on to the index-level disc, which may predispose the disc to accelerated or premature degeneration . With progressive disc degeneration, the shear and torsional stiffness of the disc may get compromised, which in addition to the loss of posterior tension band, may perpetuate to spondylolisthesis. However, from a biomechanical standpoint, it remains unclear how different grades of disc degeneration may be implicated in the progression of a lytic defect to spondylolisthesis. Yong-Hing and Kirkaldy-Willis (1983) suggested that a degenerative progression occurs in the intervertebral disc from a state of hypermobility in the early stages of degeneration to hypomobility in the later stages . Disc degeneration grades III and IV (Thomson’s grading) characterised by nuclear clefts, radial and concentric tears in the annulus fibrosus have been reported to cause kinematic instability in lumbar motion segments, whereas grade V degeneration characterised by disc height collapse and osteophyte formation resulted in biomechanical stability in the segments .
Here we use computed tomography (CT) data from a healthy human subject and nonlinear finite element modelling (FE) to assemble three-dimensional models of the lumbosacral spine to investigate the role of different grades of disc degeneration in the progression of a bilateral L5 lytic defect to spondylolisthesis. We hypothesised that a mildly degenerated disc is not able to arrest hypermobility following a bilateral lytic defect in the L5 vertebra and therefore the defect is disposed to progress to spondylolisthesis; however, in a severely degenerated disc, hypomobility is sufficient to stabilise the motion segment and limit the progression of the defect to spondylolisthesis.
Image segmentation and 3D model generation from CT data
Prior approval from UNSW Human Research Ethics Advisory Panel was obtained for the use of retrospectively acquired CT data from a human subject (UNSW-NRR-HC16754). High-resolution lumbosacral spine CT data (437 axial cuts, 512 × 512 pixel resolution, voxel dimensions: 0.30 × 0.30 × 0.50 mm) from an anonymised healthy male subject (26 years old) were obtained in DICOM (Digital Imaging and Communications in Medicine) file format from Carl Bryan Radiology (St. George Private Hospital, Sydney, Australia). Using previously published protocols, the CT data were imported into image processing software Avizo Standard (vers. 8.1, FEI Visualization Sciences Group, Hillsboro, USA) for the segmentation of anatomical regions of interest, and subsequent generation and refinement of surface and volumetric meshes . The segmented CT data were manipulated to generate a bilateral L5 spondylolytic defect model by deleting pixels from the L5 pars region (NOR LYTIC), creating an approximately 2-mm (mm) wide fracture gap between the bony fragments. Three more L5 spondylolysis models were built representing: mild disc degeneration (M-DEG LYTIC), mild disc degeneration with a disc height collapse (M-DEG-COL LYTIC), and severe disc degeneration (S-COL LYTIC) as illustrated in Fig. 1.
The intervertebral disc is a composite structure comprising the annulus fibrosus and the nucleus pulposus, which degenerates accordingly with morphological and material property changes to both. In accordance with published material property data for normal and degenerated discs, the shear modulus for the annulus fibrosus (AF) was increased (mild: + 25%, severe: + 50%) and the Poisson’s ratio for the nucleus pulposus (NP) was decreased (mild: − 13%, severe: − 34%) compared with the INTACT state in the FE models representing L5-S1 disc degeneration [23,24,25]. Disc height collapse was modelled as a 50% reduction in disc height measured in the mid-sagittal plane. The segmented CT data were manipulated to represent disc height collapse by moving the sacrum superiorly to reduce sagittal mid-disc height at the L5-S1 level from 14.35 mm to 6.64 mm while preserving the L5-S1 facet articulation.
Modelling annulus Fibres and ligaments
Volumetric meshes in Nastran file format (.nas) representing different variants of the lumbosacral spine model were imported into the FE modelling software Strand7 (vers. 2.4.6, Strand7 Pty. Ltd., Sydney, Australia). The annulus fibrosus was modelled per previously published protocols as a composite structure comprising concentric layers (n = 4) of criss-cross collagen fibres embedded within a homogenous ground substance, with the ends of the fibres rigidly anchored in the superior and inferior endplates . Ligaments were modelled using cylindrical beam elements, with the attachment and insertion sites based on previously recorded anatomical observations and published protocols [22, 26]. The type and number of elements used in assembling the five FE models are presented in Table 1.
Facet joint articulations and pars Interarticularis fracture gap
The compressive load bearing characteristics of the bony articulating pillars at each facet joint was modelled using nonlinear Point Contact-Tension elements in Strand7 (n = 5 per joint), which were normally oriented and uniformly distributed over the bony articulating surfaces. The bilateral pars interarticularis fracture gap in the NOR LYTIC, M-DEG LYTIC, M-DEG-COL LYTIC, and S-COL LYTIC models was connected using nonlinear Point Contact-Zero Gap elements (n = 25 each side, compressive stiffness only) to allow for load transfer between the fractured fragments in the event of gap closure during simulated bending motions.
Loads and boundary constraints
In all the models, a centre node on the anterior surface of the sacral mass was fixed in all translational and rotational degrees of freedom (DOFs). Bending motions were simulated using a cross-beam construct accommodated on the L1 superior endplate by means of a surface cap, both of which were assigned material properties of stainless steel (E = 200 GPa, ν =0.25). A force couple was applied to the anterior and posterior ends of the cross-beam to simulate flexion (Fx) and extension (Ex) bending. The models were loaded in pure unconstrained moments (without any compressive preload) with step-wise increments in load from 1.0 N-metre (Nm) to 10 Nm. The pre-processed FE models were solved for geometry, material, and boundary nonlinearities using the nonlinear static solver in Strand7.
Material property values
Previously published material property values calibrated against in vitro biomechanical testing data were assigned to brick elements representing the intervertebral disc, beam elements representing the primary ligaments, and non-linear point contact elements representing facet articulation between the bony pillars (see Additional file 1). 
All the results were analysed at peak Fx and Ex bending loads (10 Nm) in the five FE models.
L5-S1 range of motion (RoM)
The L5-S1 Fx and Ex RoM results are presented in Fig. 2. From the baseline INTACT state, the greatest increase in Fx and Ex RoMs was observed in the M-DEG LYTIC model (Fx: 7.2° (INTACT) to 12.4° (M-DEG LYTIC); Ex: 7.0° (INTACT) to 9.5° (M-DEG LYTIC)).
L5-S1 Interpedicular kinematics
The L5-S1 interpedicular kinematics in Fx and Ex were evaluated per published protocols [16, 27]. The in-plane component of the interpedicular travel (IPT) vector was projected onto the L5-S1 mid-discal plane to evaluate the relative posteroanterior translatory motions between the pedicles on the adjacent vertebrae.
The L5-S1 mid-discal plane projections of the IPT vector in Fx and Ex in the five models are presented in Fig. 3. From the baseline INTACT state, the greatest increase in IPT projections on to the mid-disc plane in Fx and Ex were observed in the M-DEG LYTIC model (Fx: 3.5 mm (INTACT) to 5.4 mm (M-DEG LYTIC); Ex: 1.9 mm (INTACT) to 2.7 mm (M-DEG LYTIC)).
Normal stresses in the L5-S1 AF
Colour coded compressive (−) and tensile (+) stress distribution in the L5-S1 AF in the mid-axial plane of the L5-S1 disc during Fx and Ex is presented in Fig. 4. The average values for the compressive and tensile stresses in the L5-S1 AF during Fx and Ex is presented in Fig. 5.
The greatest increase in the average normal stress from the baseline INTACT state was observed in M-DEG-COL LYTIC and S-COL LYTIC models. During Fx, average compressive stress increased from 0.11 (± 0.20) MPa in the INTACT state to 0.54 (± 1.22) MPa in the M-DEG-COL LYTIC model and 0.56 (± 1.23) MPa in the S-COL LYTIC model. During Ex, average compressive stress increased from 0.15 (± 0.21) MPa in the INTACT state to 0.68 (± 1.60) MPa in the M-DEG-COL LYTIC model and 0.69 (± 1.57) MPa in the S-COL LYTIC model.
Shear stresses in the L5-S1 AF
Colour coded shear stress distribution in the L5-S1 AF in the mid-axial plane of the L5-S1 disc during Fx and Ex is presented in Fig. 6. The average values for the posteroanterior and anteroposterior shear stress in the L5-S1 AF during Fx and Ex is presented in Fig. 7.
The greatest increase in shear stresses from the baseline INTACT state was observed in M-DEG LYTIC models. During Fx, average anteroposterior shear stress increased from 0.015 (± 0.016) MPa in the INTACT state to 0.04 (± 0.046) MPa in the M-DEG LYTIC model. During Ex, average posteroanterior stress increased from 0.015 (± 0.014) MPa in the INTACT state to 0.029 (± 0.029) MPa in the M-DEG LYTIC model.
Axial strain in iliolumbar ligament (ILL) fibres
Table 2 Shows mean axial strain in the ILL fibres during peak Fx load. The ILL fibres did not experience any positive strain during peak ex load.
The main objective of this biomechanical study was to examine the role of disc degeneration in the progression of a bilateral L5 spondylolytic defect to spondylolisthesis. The natural course of isthmic spondylolisthesis has been associated with disc degeneration and spontaneous stabilisation of the olisthetic segment, but the question whether disc degeneration is the cause or the consequence of vertebral slippage remains open .
The results from the present study confirm our hypothesis that in the presence of a bilateral L5 lytic defect, motion segments with a mildly degenerate index-level disc are at a greater risk for vertebral slippage compared to motion segments with a severely degenerate index-level disc. Compared with the baseline INTACT model, M-DEG LYTIC model experienced the greatest increase in kinematics (Fx ROM: 73% ↑; Fx IPT mid-discal projection: 53%↑), shear stresses in the annulus (Fx anteroposterior: 163%↑; Fx posteroanterior: 31%↑), and average strain in the ILL fibres (Fx: 90%↑). The S-COL LYTIC model experienced a decrease in mobility (Fx ROM: 48% ↓; Fx IPT mid-discal projection: 69%↓) and an increase in normal stresses in the annulus (Fx Tensile: 170%↑; Fx Compressive: 397%↑) compared with the INTACT model. The results are corroborated by previous clinical findings for conservatively treated L5-S1 lytic defect patients which revealed immobility of the motion segment with severe degeneration of the index-level disc . In the light of Kirkaldy-Willis model for degenerative changes in the lumbar spine and results from the present study, we posit that severe disc degeneration of the index-level disc is one of the restabilisation mechanisms available to the L5-S1 motion segment in order to mitigate increased intervertebral motions and shear stresses due to a bilateral L5 lytic defect . The results also indicate that the disc height collapse and not the degenerative changes in the stiffness properties of the disc per se, play a predominant role in the restabilisation of the L5-S1 segment. With a 50% loss in disc height, the L5-S1 segment became less mobile, offloaded the ILL fibres, and experienced lower shear stresses compared with the intact disc height models. No significant difference in results was noted between M-DEG-COL LYTIC and S-DEG-COL LYTIC models. The collapsed disc height models, however, experienced high normal stresses compared with the intact disc height models suggesting that loss in disc height redistributed loads within the motion segment, decreasing the shear component of the intradiscal forces and increasing the normal component.
Over the age of 25 years, the prevalence of index-level disc degeneration in lytic defect patients was found to be significantly greater compared with an age-and-sex matched control group . Although progression of a lytic defect to spondylolisthesis in adults depends on numerous factors such as spinopelvic balance, lower-lumbar muscle strength, repetitive flexion-extension activities, strength of the iliolumbar ligament; the most important structure resisting the vertebral slippage is the index-level disc [4, 6, 29, 30]. The health (or lack thereof) of the index-level disc is therefore critical to resisting the vertebral slippage. Two distinct pathways may be available to a L5-S1 segment with a bilateral lytic defect and a mildly degenerate index-level disc:
Within the limits of biomechanical stability provided by the active and passive spinal elements, the index-level disc may progressively lose height resulting in a decrease in intervertebral motions and redirection of the intradiscal shear stresses to normal stresses, until biomechanical stabilisation is achieved.
Beyond the limits of biomechanical stability provided by the active and passive spinal elements, the defect may progress to spondylolisthesis followed by a loss in disc height, in small and closed-loop increments, until biomechanical stabilisation is achieved.
From the onset of a bilateral L5 lytic defect, the exact duration of unstable and restabilisation phases may vary from patient to patient. Why some L5-S1 motion segments will adopt one pathway over the other may depend on genetic, anatomical, and lifestyle factors, and remains to be clinically explored. If the slippage progression is relatively rapid and intractable radicular pain associated with nerve root entrapment is evident, surgical intervention may be necessary to achieve biomechanical stability .
Progressive degeneration of the intervertebral disc is characterised by a cascade of cellular, biochemical, morphological, and functional changes. The sequence of these changes and the interplay between them are largely unknown. In the present study, a rather simplified approach was used to model progressive degeneration of the L5-S1 disc, altering only the structure (height) and material properties. The material properties of the AF and NP tissues were altered for different grades of disc degeneration in accordance with published data . Whilst the shear modulus for the AF tissues was increased with progressive degeneration, Poisson’s ratio for the NP tissues was decreased to simulate the loss in incompressibility with progressive degeneration.  A 50% loss in mid-sagittal plane disc height was modelled to simulate a collapsed disc.
Some limitations to this study were noted. Poroelasticity and degeneration related structural defects (Schmorl’s node, annular tears, nuclear clefts, and rim lesions) in the disc were not modelled, which may have a bearing on the results. Instead of assembling multiple FE models representing a gradual change in the disc height, only a 50% collapse in disc height was modelled. The AF was modelled as a composite structure with only four concentric layers of collagen fibres (without translamellar bridges) embedded within a homogenous ground substance.
In the presence of a bilateral L5 spondylolytic defect, a mildly degenerate index-level disc experienced greater intervertebral motions and shear stresses compared with a severely degenerate index-level disc during flexion and extension bending. Disc height collapse, with or without severe degenerative changes in the stiffness properties of the disc, is one of the plausible restabilisation mechanism available to the L5-S1 motion segment to mitigate increased intervertebral motions and shear stresses due to a bilateral L5 lytic defect.
Anterior Longitudinal Ligament
Digital Imaging and Communications in Medicine
Degrees of Freedom
An intact finite element model of L1-S1 spine without any defect or degeneration
First Lumbar Vertebra
Fifth Lumbar Vertebra
A Spondylolytic Defect
- M-DEG LYTIC:
A finite element model of L1-S1 spine with a bilateral L5 lytic defect and mild degenerative changes in the stiffness properties of the L5-S1 disc
- M-DEG-COL LYTIC:
A finite element model of L1-S1 spine with a bilateral L5 lytic defect, mild degenerative changes in the stiffness properties of the L5-S1 disc and 50% disc height collapse
- NOR LYTIC:
A finite element model of L1-S1 spine with a bilateral L5 lytic defect and a normal L5-S1 disc
Posterior Longitudinal Ligament
Range of Motion
- S-COL LYTIC:
A finite element model of L1-S1 spine with a bilateral L5 lytic defect, severe degenerative changes in the stiffness properties of the L5-S1 disc and 50% disc height collapse
University of New South Wales
Fredrickson BE, Baker D, McHolick WJ, Yuan HA, Lubicky JP. The natural history of spondylolysis and spondylolisthesis. J Bone Joint Surg Am. 1984;36:699–707.
Leone A, Cianfoni A, Cerase A, Magarelli N, Bonomo L. Lumbar spondylolysis: a review. Skelet Radiol. 2011;40:683–700.
Hresko MT, Labelle H, Roussouly P, Berthonnaud E. Classification of high-grade spondylolistheses based on pelvic version and spine balance: possible rationale for reduction. Spine (Phila Pa 1976). 2007;32:2208–13.
Labelle H, Mac-Thiong JM, Roussouly P. Spino-pelvic sagittal balance of spondylolisthesis: a review and classification. Eur Spine J. 2011;20(Suppl 5):641–6.
Mac-Thiong JM, Labelle H. A proposal for a surgical classification of pediatric lumbosacral spondylolisthesis based on current literature. Eur Spine J. 2006;15:1425–35.
Roussouly P, Gollogly S, Berthonnaud E, Labelle H, Weidenbaum M. Sagittal alignment of the spine and pelvis in the presence of L5-s1 isthmic lysis and low-grade spondylolisthesis. Spine (Phila Pa 1976). 2006;31:2484–90.
Sairyo K, et al. Vertebral forward slippage in immature lumbar spine occurs following epiphyseal separation and its occurrence is unrelated to disc degeneration: is the pediatric spondylolisthesis a physis stress fracture of vertebral body? Spine (Phila Pa 1976). 2004;29:524–7.
Sairyo K, Katoh S, Sakamaki T, Inoue M, Komatsubara S, Ogawa T, et al. A review of the pathomechanism of forward slippage in pediatric spondylolysis: the Tokushima theory of growth plate slippage. J Med Investig. 2015;62:11–8.
Vialle R, Ilharreborde B, Dauzac C, Lenoir T, Rillardon L, Guigui P. Is there a sagittal imbalance of the spine in isthmic spondylolisthesis? A correlation study. Eur Spine J. 2007;16:1641–9.
Ikata T, Miyake R, Katoh S, Morita T, Murase M. Pathogenesis of sports-related spondylolisthesis in adolescents. Radiographic and magnetic resonance imaging study. Am J Sports Med. 1996;24:94–8.
Szypryt EP, Twining P, Mulholland C, Worthington BS. The prevalence of disc degeneration associated with neural arch defects of the lumbar spine assessed by magnetic resonance imaging. Spine (Phila Pa 1976). 1989;14:977–81.
Seitsalo S, Schlenzka Z, Poussa M, Österman K. Disc degeneration in young patients with isthmic spondylolisthesis treated operatively or conservatively: a long-term follow-up. Eur Spine J. 1996:393–7.
Floman Y. Progression of lumbosacral isthmic spondylolisthesis in adults. Spine (Phila Pa 1976). 2000;25:342–7.
Lyras DN, Tilkeridis K, Stavrakis T. Progression of spondylolysis to isthmic spondylolisthesis in an adult without accompanying disc degeneration: a case report. Acta Orthop Belg. 2008;74:141–4.
Stone AT, Tribus CB. Acute progression of spondylolysis to isthmic spondylolisthesis in an adult. Spine (Phila Pa 1976). 2002;27:E370–2.
Chamoli U, Chen AS, Diwan AD. Interpedicular kinematics in an in vitro biomechanical assessment of a bilateral lumbar spondylolytic defect. Clin Biomech (Bristol, Avon). 2014;29:1108–15.
Mihara H, Onari K, Cheng BC, David SM, Zbelick TA. The biomechanical effects of spondylolysis and its treatment. Spine (Phila Pa 1976). 2003;28:235–8.
Fan J, Yu GR, Liu F, Zhao J, Zhao WD. A biomechanical study on the direct repair of spondylolysis by different techniques of fixation. Orthop Surg. 2010;2:46–51.
Wiltse LL. The effect of the common anomalies of the lumbar spine upon disc degeneration and low back pain. Orthop Clin North Am. 1971;2:569–82.
Yong-Hing K, Kirkaldy-Willis WH. The pathophysiology of degenerative disease of the lumbar spine. Orthop Clin North Am. 1983;14:491–504.
Tanaka N, An HS, Lim T, Fujiwara A, Jeon C, Haughton VM. The relationship between disc degeneration and flexibility of the lumbar spine. Spine J. 2001;1:47–56.
Ramakrishna VAS, Chamoli U, Viglione LL, Tsafnat N, Diwan ADD. The role of sacral slope in the progression of a bilateral Spondylolytic defect at L5 to spondylolisthesis: a biomechanical investigation using finite element analysis. Global Spine Journal. 2017; https://doi.org/10.1177/2192568217735802.
Malandrino A, Planell JA, Lacroix D. Statistical factorial analysis on the poroelastic material properties sensitivity of the lumbar intervertebral disc under compression, flexion and axial rotation. J Biomech. 2009;42:2780–8.
Nguyen AM, Johannessen W, Yoder JH, Wheaton AJ, Vresilovic EJ, Borthakur A, Elliott DM. Noninvasive quantification of human nucleus pulposus pressure with use of T1rho-weighted magnetic resonance imaging. J Bone Joint Surg Am. 2008;90:796–802.
Wang Y, Chen H, Zhang L, Liu J, Wang Z. Influence of degenerative changes of intervertebral disc on its material properties and pathology. Chin J Traumatol. 2012;15:67–76.
Behrsin JF, Briggs CA. Ligaments of the lumbar spine: a review. Surg Radiol Anat. 1988;10:211–9.
Chamoli U, Chen AS, Diwan AD. Letter to the editor regarding the article " Interpedicular travel in the evaluation of spinal implants: an application in posterior dynamic stabilization" by D.J Cook, MS Yeager, and BC Cheng: Spine 2012; 37(11): 923–931. Spine (Phila Pa 1976). 2014;39(11):921.
Kirkaldy-Willis WH, Farfan HF. Instability of the lumbar spine. Clin Orthop Relat Res. 1982:110–23.
Motley G, Nyland J, Jacobs J, Caborn DNM. The pars Interarticularis stress reaction, spondylolysis, and spondylolisthesis progression. J Athl Train. 1998;33:351–8.
Aihara T, Takahashi K, Yamagata M, Moriya H, Shimada Y. Does the iliolumbar ligament prevent anterior displacement of the fifth lumbar vertebra with defects of the pars? J Bone Joint Surg Br. 2000;82:846–50.
The authors thank Emma Wheatley, Karen Ruut, and Carl Bryant from CJ Bryant Radiology, St. George Private Hospital, Sydney (Australia) for providing access to de-identified patient CT data.
This work was supported by a research grant from AO Spine Foundation Asia Pacific (AOSAUNZ(R) 2016–01) to ADD, a Taste of Summer Research Scholarship from UNSW Australia to VR, an International Postgraduate Research Scholarship from UNSW Australia and DIISRTE to UC, and internal research funds from Spine Service.
Availability of data and materials
The datasets used during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Prior approval from the University of New South Wales Human Research Ethics Advisory Panel was obtained for the use of retrospectively acquired computed tomography data from a human subject (UNSW-NRR-HC16754). A written consent was obtained from the patient for the use of their de-identified medical records for research purposes.
Consent for publication
A written consent was obtained from the patient for the use of their de-identified medical records for research purposes.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Mild (not severe) disc degeneration is implicated in the progression of L5 spondylolysis to spondylolisthesis. This supplementary information file contains data on material property values assigned to various elements used in building the five FE models of lumbosacral spine in this study. (DOCX 36 kb)