Introduction

Lumbar laminectomy is a commonly used treatment for symptomatic degenerative lumbar spinal stenosis [6]. Although the impinged nerves are decompressed and neurological symptoms, such as low back pain, sciatica, claudication, motor, sensory and reflex activity, often improve following lumbar laminectomy, it can lead to symptomatic postoperative lumbar instability or even postoperative failure of the spinal motion segment [14]. A well-known complication of lumbar laminectomy is excessive shear displacement in the intervertebral joint, leading to postoperative spondylolysis or spondylolisthesis [7]. Symptomatic clinical instability justifies reoperation to stabilize and fuse the unstable segment [8]. When residual strength and shear stiffness (SS) of the lumbar spine after laminectomy can be predicted, this may support patient selection for additional spinal stabilization. In other words, based on predicted residual shear properties, the surgeon may decide whether or not to combine laminectomy with (instrumented) fusion techniques.

Previously, we showed in an in vitro experiment that laminectomy resulted in a substantial decrease of SS and shear force to failure (SFF) of lumbar spinal segments [1]. However, the biomechanical behaviour of a spinal motion segment following laminectomy will likely also depend on disc degeneration, facet joint degeneration, Modic changes, Schmorl’s nodes, intervertebral disc and pedicle geometry, and facet joint angles.

To our best knowledge, there is a lack of information in literature, demonstrating correlations between these various anatomical and clinical parameters and the biomechanical behaviour of a spinal motion segment following lumbar laminectomy. In this study, we aim to assess the relationship between various anatomical and clinical parameters, and in vitro strength and SS of a lumbar spinal segments either untreated or following facet sparing laminectomy. A total of ten spines (Th12–L5) were used. Ten segments remained untreated (five times L2–L3 and five times L4–L5) and ten segments were treated with laminectomy (five times L2–L3 and five times L4–L5).

We hypothesized that multiple independent variables, together, determine shear biomechanics of a lumbar spinal segment either intact or treated with laminectomy. Identification of these determinants may enable prediction of shear biomechanics in the future, which may support surgical decision-making.

Methods

Specimens

Thoracolumbar spines (T12–L5) were harvested from freshly frozen (−20°C) human cadavers (mean age 72.1 years, range 53–89 years). None of the deceased subjects had any history of spinal injury, spinal surgery or spinal metastatic disease. The spines were thawed before assessment and biomechanical testing. Excessive soft tissue and muscle tissue were carefully removed, keeping the anterior and posterior longitudinal ligaments as well as the facet joints intact (Fig. 1).

Fig. 1
figure 1

Human thoracolumbar spine (T12–L5) with laminectomy at level L4

Parameters

For assessment of spines, we used clinically relevant and methodologically validated parameters of lumbar spinal degeneration as recommended by the European Spine Society [9]. Grading methods for disc degeneration with an intraclass correlation coefficient or an interobserver κ > 0.60 [5, 13, 17] were included. For facet joint degeneration, grading schemes [9] with an intraclass correlation coefficient or interobserver κ > 0.40 were used in the present study [16, 26].

Magnetic resonance imaging (MRI, Siemens© Symphony 1.5 T: Syngo MR A30, software NUMARIS/4, Berlin, Germany) of lumbar spines was performed to assess intervertebral disc degeneration according to Griffith and Pfirrmann [5, 17] and facet joint degeneration according to Weishaupt [26]. Disc degeneration, (including narrowing and osteophytes, respectively Lanes 1 and 2) [13, 27] and facet joint degeneration [16] of levels L2–L3 and L4–L5 were also assessed based on radiographs (Sedical© Digital Vet. DX-6, Arlington Heights, IL, USA). Furthermore, MRI was used to assess the presence of Modic changes [15] and Schmorl’s nodes [18] and to determine intervertebral disc and pedicle geometry and facet joint angles [2]. Disc geometry included: disc length, width, height, surface area, and volume. Disc surface area, disc volume and pedicle diameter were calculated assuming an elliptic shape (surface = 1/4π × length × width). For pedicle diameter, an average of left and right pedicles was taken for the top (L2 or L4) and bottom (L3 or L5) of each segment. Mean facet joint angle was calculated by averaging left and right angles per segmental level (L2–L3 or L4–L5) while facet joint angle differences or tropism was determined by calculating the difference between left and right facet joint angles. Segmental frontal surface area (FA), defined in cm2, bone mineral content (BMC in g) and bone mineral density (BMD in g/cm2) of lumbar spinal sections (L2–L3 and L4–L5) were measured with dual X-ray absorptiometry (DXA, Hologic© QDR 4500 Delphi DXA scanner, Waltham, MA, USA) in anteroposterior direction. All assessments were performed using Osirix software (Osirix©, version 3.8.1., Pixmeo SARL, Geneva, Switzerland).

Specimen preparation and biomechanical testing

L2–L3 and L4–L5 motion segments were isolated from each spine. Subsequently, laminectomy was performed at level L2 of five randomly chosen spines, and at level L4 of the remaining five spines. Laminectomy, analogous to standard clinical practice, was performed by removing the spinous process and part of the lamina, leaving the pars interarticularis intact. During preparation, examination, and biomechanical testing, specimens were kept hydrated using 0.9% saline-soaked gauzes. Thoracolumbar spines with bridging osteophytes, assessed on anteroposterior, lateral and oblique radiographs, were excluded from this study. After sectioning spines in L2–L3 and L4–L5 motion segments, the motion segments were potted in a casting-mould using low melting point (48°C) bismuth alloy (Cerrolow-147; 48.0% bismuth, 25.6% lead, 12% tin, 9.6% cadmium, and 4% indium). The upper and lower vertebral bodies were fixed securely into the alloy by adding screws into the vertebral body. Screw fixation was reinforced with orthopaedic bone cement (Simplex, Stryker©, Kalamazoo, MI, USA). The disc was placed parallel to the flat surface of the bismuth. Discs were placed parallel based on the visual inspection. Because muscle tissue was thoroughly and carefully removed, the intervertebral disc and corresponding endplates were clearly visible. All articulating parts were kept free. The casting mould was placed in a hydraulic materials testing machine (Instron©, model 8872, Norwood, Canada) [1, 23, 24]. The caudal vertebral body was fixed on a plateau that allowed movement in axial and transverse directions only. Transverse movements were allowed, so segments were able to find their physiological motion patterns and to correct for possible differences in embedding. Segments were loaded with a continuous axial compressive force of 1600 N [23, 24], applied using a pneumatic cylinder that had been calibrated using a load cell (Hottinger Baldwin Messtechnik©, Force Transducer Type C2, Darmstadt, Germany). Since compression was applied in a purely axial direction, bending moments were minimized. The level of compression simulated the force during bending, a condition in which high shear loading of the lumbar spinal segments typically occurs [23]. Subsequently, while maintaining the axial load, anterior shear load was applied with a constant rate of 2.0 mm/min on the casting mould containing the cranial vertebral body, until failure of the vertebral motion segment [24]. This test set-up was similar to mechanical testing by Bisschop et al. [1], van Solinge et al. [24] and van Dieën et al. [23]. An anterior shear force was used since it corresponds to the loading direction in vivo [1012, 22]. The test was stopped after hearing a crack or after a large force reduction was seen. Shear force and displacement were digitized and stored at 100 samples per second (Instron© Fast Track 2, Norwood, Canada).

For each of the 20 motion segments tested, SFF was determined. SFF was defined as the point at which maximum load was recorded in the load–displacement curves for each specimen. These data were presented previously [1]. Shear yield force (SYF) was defined as the point at which shear load caused a decrease in stiffness, i.e. a decrease in the slope of the load–displacement curve. Average SS was calculated from the load–displacement curve, between 25 and 50% of the SFF. SS was estimated by means of a least squares fit of a straight line through the data with the slope of the regression line representing stiffness. The deformation in this region was linear, with an r 2 > 0.943 (Table 1) between load and displacement for all motion segments. All analyses were performed using computer programs written in Matlab (Mathworks ©, Natick, MA, USA).

Table 1 Overview of specimens and biomechanical outcomes per tested segment

Statistical methods

Statistical analysis was performed based on two separate groups. The first group contained untreated segments (5× L2–L3 and 5× L4–L5) while the second group consisted of segments with laminectomy (5× L2–L3 and 5× L4–L5).

Independent variables were classified as: general variables, intervertebral disc geometry (MRI), pedicle geometry (MRI), facet joint orientation (MRI), bone characteristics (DXA), intervertebral disc degeneration classifications (MRI), intervertebral disc and facet joint degeneration (Radiographs), facet joint degeneration (MRI) and other (MRI). These classes of variables are specified in Table 2.

Table 2 Overview of correlations (p values, two tailed <0.05: in bold) between independent and dependent variables in untreated and treated segments

First, relations between independent and dependent variables (SS, SYF and SFF) were tested for each individual variable. For dichotomized independent variables (segment, sex, Modic changes [15] and Schmorl’s nodes [18], independent-sample t tests were used while Pearson’s coefficient of correlation was determined for continuous and ordinal values. Note that it was thus assumed that ordinal variables (Pfirrmann [17], Griffith [5], Lane 1 [13], Lane 2 [13], Wilke [27], Pathria [16] and Weishaupt [26]) represent a linear degree of severity.

When independent variables were associated with a dependent variable, here defined as independent-sample t test: p < 0.05 or as a bivariate correlation with a significance level of: p < 0.05, they were used for the combined statistical models.

Before final analysis was performed, all independent variables were checked for correlations with each other. In case a correlation >0.7 with a p < 0.05 was found, the independent variable with the strongest effect on the specific dependent variable was included in the model. Finally, backward linear regression techniques were used to create final statistical models per dependent variable per treatment group.

Results

All specimen characteristics and biomechanical properties for segments with and without laminectomy are presented in Table 1. As shown in Fig. 2, lumbar laminectomy had a substantial declining effect on SS (23.7%), SYF (41.1%) and SFF (44.3%).

Fig. 2
figure 2

Effects of lumbar laminectomy on shear biomechanics, showing a substantial decrease of shear stiffness (23.7%), shear yield force (41.1%) and shear force to failure (44.3%) following laminectomy

Table 2 gives an overview of correlations between independent and dependent variables of segments with and without laminectomy. Some of the general variables, as presented in Table 1 (sex and age for untreated segments, and segment level and sex for treated segments), were related to strength parameters (SYF and/or SFF) in both groups. In the untreated segments, SFF was found to be lower for female specimens (2.284 N for male versus 1.351 N for female). For the treated segments, SFF was also found to be lower for female specimens (1.346 N for male versus 646 N for female). SYF in treated segments proved to be level dependent (L2–L3: 428 N versus L4–L5: 851 N).

For segments treated with laminectomy, three out of five intervertebral disc geometry variables (i.e., length, surface and volume) were significantly related to biomechanical shear properties (SYF; all three, for SFF; length only). In contrast, biomechanics of untreated segments were unrelated to intervertebral disc geometry. The opposite was true for pedicle geometry and facet joint orientation. Pedicle sections and facet joint angle difference correlated significantly to respectively, SYF and SS in untreated segments but did not correlate with biomechanical outcomes in treated segments.

For both groups, bone characteristics measured with DXA, were strongly related to shear strength parameters (SYF and SFF), but not to stiffness (SS). Like intervertebral disc geometry, intervertebral disc degeneration was predictive for biomechanics (SS) of spinal segments with laminectomy. This was consistent over imaging methods and classification schemes (MRI; Pfirrmann [17] and Griffith [5] or radiographs; Wilke [27]), although not significant for radiographs (correlation: 0.558, p value: 0.094). However, in contrast to intervertebral disc geometry, these intervertebral disc degeneration scores were not related to strength (SYF and SFF). Finally, Modic changes [15] and Schmorl’s nodes [18] were not related to shear biomechanics of spinal segments with or without laminectomy.

Results of the backward linear regression, using determinants of spine biomechanics, which were identified (based on a p < 0.05) in Table 2, are presented in Table 3. All models, describing strength parameters (SYF and SFF) consisted of two independent variables.

Table 3 Overview of backward linear regression models per dependent variable in untreated and in treated segments based on significant correlation coefficients found in Table 2

SYF and SFF could accurately be predicted by the final statistical model for untreated segments (r 2 value respectively: 0.83; 0.83). Age and pedicle geometry remained in the model for SYF, while for SFF, the final model consisted of DXA parameters (frontal area and BMC) only. For segments treated with laminectomy, SYF and SFF could be predicted from independent variables with r 2 values of 0.81 (intervertebral disc volume and segment) and 0.77 (sex and intervertebral disc length), respectively. SS was less accurately predicted with only a single variable remaining in the model for both untreated segments (facet joint angle difference; r 2 = 0.50) and segments with laminectomy (degeneration score Lane-2 (osteophytes); r 2 = 0.53).

Discussion

The aim of this study was to identify parameters, that predict spinal shear properties before and after laminectomy, in order to determine which of these parameters may prognosticate spinal instability following lumbar laminectomy.

For characterization of the spinal motion segments, we used commonly applied grading systems to assess disc degeneration [5, 13, 17, 27], facet joint degeneration [9, 16, 26], Modic changes [15] and Schmorl’s nodes [18] based on MRI and radiographs. Furthermore, we measured intervertebral disc and pedicle geometry, facet joint angles [2] and bone characteristics (BMC: bone mineral content; BMD and total segmental surface area on DXA defined as frontal area: FA). These parameters all potentially affect strength and SS of the lumbar spinal motion segment before and/or after treatment with laminectomy and can be determined in clinical practice.

We showed that multiple variables are related to spinal shear properties in intact lumbar segments and lumbar segments treated with laminectomy. Statistical models with these parameters as independent variables predicted shear biomechanics, with moderate to very good accuracy with r 2 values varying from 0.50 to 0.83 (without laminectomy) and from 0.53 to 0.81 (with laminectomy). Particularly, strength parameters (SYF and SFF) in both untreated and treated segments could be predicted with good to very good accuracy. Prediction of SS was only moderately accurate.

The tests on individual variables (Table 2) showed that, for untreated segments, pedicle geometry was related to SYF and facet joint orientation to SS. In contrast, for segments with laminectomy, intervertebral disc characteristics appeared to determine shear properties. Intervertebral disc characteristics correlated to strength and disc degeneration correlated to stiffness. For both segments with and without laminectomy, DXA assessment was found to be important, although mainly for strength parameters.

SYF might be the most critical shear property, because it marks the beginning of the irreversible deformation of a spinal motion segment, signalling the appearance of the first soft tissue and or trabecular bone lesions [21]. We expect that when shear loading crosses the yield point, sub-clinical damage will occur. Such damage may, at a later stage, lead to symptomatic spondylolisthesis. Unlike SYF, SFF marks, as the description suggests, complete and irreversible failure of spinal motion segments. SFF describes an acute clinically relevant situation. Therefore, SYF and SFF represent different clinical value. In untreated segments SYF depended mainly on pedicle geometry, while SFF strongly correlated with DXA parameters (Table 2). For treated segments, both SYF and SFF were correlated with intervertebral disc geometry and DXA parameters (Table 2) and both parameters could be predicted quite accurately (Table 3).

In this study, SS was only moderately predictable (r 2 values: 0.50 and 0.53). We assumed that the intervertebral disc has a large contribution to this biomechanical parameter. This assumption was corroborated by the results (Table 2). Degenerative parameters proved to be strongly correlated with SS in treated segments. In other words, laminectomy leads to a shift in load bearing, from the pars interarticularis to the intervertebral disc. Unfortunately, we could only study the morphology and degeneration of the intervertebral disc on MRI and radiographic imaging. A more specific (histological) analysis of the state of the intervertebral disc may strengthen correlations [19], but may not be clinically applicable.

For stiffness, r 2 values of only 0.50 and 0.53 were found. As stated earlier, stiffness was mainly determined by degenerative parameters, such as disc degeneration. The fact that these parameters are based on visual assessment and have an ordinal character possibly explains their lower predictive value, compared to directly measured continuous variables such as BMD and BMC.

In our protocol, both BMD and BMC were studied. BMD is often used as a clinical parameter. However, BMC, can also be used to express the bone mineral content since it integrates information on bone density and vertebral dimensions. BMC is defined as BMD (g/cm2) multiplied by the total segmental surface area (FA) of the spinal segment (cm2) and is expressed in grams. We therefore decided to include both parameters as a factor that prognosticates instability following lumbar laminectomy.

We found a substantial difference between male and female specimen considering SFF, in both treated and untreated segments. However, considering the limited number of tested specimens, we cannot draw any conclusions from these findings.

In vivo, muscle forces are very important [20]. Muscle forces are the main generators of compression and shear forces. We simulated the effect of muscle forces on the spine using static 1600 N compressive force and an increasing shear force imposed by the material testing machine. The chosen preload of 1600 N was selected to allow for comparison with previous work [23, 24] and was a compromise between applying compression forces that are sufficiently large to simulate spinal loads that occur in vivo when large shear forces are present [1012, 22], but low enough to avoid damage due to compression forces alone [3].

One limitation of this study is that small alignment errors may have been present. Our results, however, are not likely to be very sensitive to small errors in specimen alignment. Previously, it was shown that SS and SFF were not different between specimens in neutral position and specimen in 10° of flexion [23]. Therefore, we do not expect significant changes in biomechanical outcomes when malaligning segments.

Another limitation of this study was that we did not investigate the nature of failure. Van Solinge et al. [24] investigated types of failure, occurring with shear loading. These failure mechanisms were similar to those found in clinical practice. Since our test setup was similar, we expect our segments to fail in similar fashion.

From a clinical point of view, laminectomy at a spinal segment that exhibits small intervertebral disc geometry, disc and facet joint degeneration and poor bone mineral density may need additional instrumental spinal stabilization to reduce the risk of post-operative instability. However, also pull out strength of spinal implants, proved to be dependent on bone mineral quality as measured by dual X-ray absorptiometry (DXA) [25] and this dependency needs to be taken into account when deciding on instrumentation.

Considering further research, we recommend to assess the parameters found to be predictive in a prospective or retrospective in vivo design. In addition to shear failure, further studies should also focus on other failure mechanisms of the human lumbar spine, including axial rotation [4].

Finally, while r 2 values, as we presented, may be too low to provide the sole basis for decisions upon surgical stabilization after laminectomy. Strength parameters (SYF and SFF) correlations were predicted with reasonable accuracy (r 2 values between 0.77 and 0.83). As currently surgeons decide based upon personal experience, a more informed choice might benefit this decision.

In conclusion, predictive models with moderate to good accuracy were found for SYF and SFF of human lumbar spinal segments with and without laminectomy. Significant loss of SS and strength are predicted by BMC, BMD, intervertebral disc geometry and degenerative parameters. Therefore, knowledge of a patient’s BMC, BMD, intervertebral disc geometry and the possible presence of osteophytes, might provide valuable information as predictors of the development of post-operative instability following lumbar laminectomy. Pedicle sections and facet geometry were not predictive for the possible development of post-operative instability following lumbar laminectomy.