Quantitative MRI to Characterize the Nucleus Pulposus Morphological and Biomechanical Variation According to Sagittal Bending Load and Radial Fissure, an ex vivo Ovine Specimen Proof-of-Concept Study

Background and context: Low back pain is a dramatic burden worldwide. Discography studies have shown that 39% of chronic low back pain patients suffer from discogenic pain due to a radial fissure of intervertebral disc. This can have major implications in clinical therapeutic choices. The use of discography is restricted because of its invasiveness and interest in it remains low as it represents a static condition of the disc morphology. Magnetic Resonance Imaging (MRI) appears to be less invasive but does not describe the biomechanical dynamic behavior of the fissure.

Purpose: We aimed to seek a quantitative MRI protocol combined with ex vivo sagittal loading to analyze the morphological and biomechanical changes of the intervertebral disc structure and stress distribution.

Study design: Proof of concept.

Methods: We designed a proof-of-concept ovine study including 3 different 3.0 T-MRI sequences (T₂-weighted, T₁ and T₂ mapping). We analyzed 3 different mechanical states (neutral, flexion and extension) on a fresh ovine spine specimen to characterize an intervertebral disc before and after puncturing the anterior part of the annulus fibrosus. We used a mark tracking method to calculate the bending angles and the axial displacements of the discal structures. In parallel, we created a finite element model to calculate the variation of the axial stress and the maximal intensity shear stress, extrapolated from our experimental boundary conditions.

Results: Thanks to an original combination of specific nuclear relaxation time quantifications (T₁, T₂) of the discal tissue, we characterized the nucleus movement/deformation into the fissure according to the synchronous mechanical load. This revealed a link between disc abnormality and spine segment range of motion capability. Our finite element model highlighted significant variations within the stress distribution between intact and damaged disc.

Conclusion: Quantitative MRI appears to provide a new opportunity to characterize intra-discal structural morphology, lesions and stress changes under the influence of mechanical load. This preliminary work could have substantial implications for non-invasive disc exploration and could help to validate novel therapies for disc treatment.

Introduction

Low back pain is the leading cause of disability worldwide (Vos et al., 2016) and represents a dramatic economic burden for western countries (Maetzel and Li, 2002; Walker et al., 2003; Dagenais et al., 2008). Back pain is frequently associated with intervertebral disc degeneration, defined as “an aberrant, cell-mediated response to progressive structural failure” (Adams and Roughley, 2006). Several pathways can lead to disc degeneration (Adams and Dolan, 2012). One of them, starts from a centrifugal (from center to periphery) and radial fissure of the annulus which alters disc stress distribution (McNally et al., 1996) and creating a stress gradient between the posterior annulus and the nucleus pulposus (Stefanakis et al., 2014). These mechanical changes can modify the cellular activities leading to dysregulation of the TIMP/MMP expression (TIMP for Tissue Inhibitor of Metallo-Proteinase and MMP for Matrix Metallo-Proteinase) (Le Maitre et al., 2004, 2007). This, in turn, leads to the acceleration of the normal ageing nucleus dehydration (Antoniou et al., 1996). Such modifications can stimulate the nociceptors naturally present in the outer third of the annulus (García-Cosamalón et al., 2010) or those which proliferate alongside the fissure (Coppes et al., 1990, 1997; Lama et al., 2018). All of these phenomena define one structural substrate of discogenic backpain.

Aiming to reverse the pathological status of the disc, multiple treatment options are available, ranging from conservative management to interventional therapies. Robust concepts supporting physiotherapy and manual therapies base their intervention on the presence of a directional preference (McKenzie, 1981; McKenzie and May, 2003; Laslett et al., 2005), which implies a dynamic disc theory. From a clinical perspective, directional preference is a direction of movement alleviating patient pain, while the other directions have no effect or worsen the pain. When observed on a patient with back pain, this type of clinical sign is specific (94%) to discogenic pain (Laslett et al., 2005) and appears to function as an effective guide for treatment (May and Aina, 2012; May et al., 2018). In addition to disc surgery, which is limited to refractory patients, interventional therapies include chemonucleolysis using chymopapain (Javid et al., 1983) – historically the first intradiscal injected drug – collagenase, chondroitinase (Ishibashi et al., 2019), ozone (Paradiso and Alexandre, 2005), radiopaque gelled ethanol (Hashemi et al., 2020), and thermocoagulation (Freeman, 2006). Other novel intradiscal therapies, such as neurotrophic growth factor (Knezevic et al., 2017), platelet-rich plasma (Li et al., 2017) and stem cell (Meisel et al., 2019) appear promising but require further and large-scale validation. Treatment indication can be outlined by correlating a clinical phenotype to a specific morphological disc profile. However, this task becomes extremely challenging when a degenerative process occurs among pain-free and healthy subjects (Brinjikji et al., 2015).

As adequate treatment requires precise diagnostics and as the standard MRI fails, one has to fall back to discography. This procedure has historically been accepted as the reference standard imaging procedure for the diagnosis of discogenic pain. By injecting a contrast agent into the central nucleus, this procedure triggers a mechanical distension of the inner annulus, reproducing the concordant patient pain (Bogduk et al., 2013). By assessing the spreading of the contrast agent with an X-Ray or a CT-scan, the clinician can identify a centrifugal nucleus radial fissure. Discography studies have demonstrated that 39% of patients with chronic low back pain present discogenic pain associated with an annular radial fissure at the level of the injured disc (Manchikanti et al., 2018). The extent of the fissure can be graded using the Dallas classification (Sachs et al., 1987). Furthermore, patient pain has been shown to be proportional to the extent of the fissure toward the periphery (Vanharanta et al., 1987; Moneta et al., 1994). However, since discography is performed on patients in prone position, it enables the physician to analyze only a static aspect of the annular fissure. It provides no information about the potential mobility of the nucleus and the dynamic evolution of the fissure. Although this procedure has been modified to avoid any degenerative process (Bogduk, 2013; McCormick et al., 2019), some studies still show that discography could provide negative side-effects (Carragee et al., 2009). Because of these definitive disadvantages, its use has been restricted to a few highly specific conditions, and substantial efforts are being placed in the development of non-invasive markers (Aprill and Bogduk, 1992).

In recent years, Thompson et al. (2009) described the connection between positive discographies on more than 2400 discs and signal intensity alterations observed via Magnetic Resonance Imaging (MRI) on T₁- and T₂-weighted scans, corresponding to inflammation of the vertebral endplate, named “MODIC sign” (Modic et al., 1984). Such signs are associated with patient pain (Thompson et al., 2009) and translate as a non-traumatic way of diagnosing discogenic backpain. However, MODIC signs describe only the inflammatory state of the endplates, and not the morphology or the dynamic behavior of the disc itself. The Pfirrmann classification is another way to describe degenerative lumbar disc (Pfirrmann et al., 2001). This classification encompasses a multi-parameter description of the disc on T₂-weighted MR images (disc height, signal intensity, etc.) allowing a grading system of the degenerative process. Several studies validate the correlation between MRI findings and histochemical composition of the intervertebral disc (Tertti et al., 1991; Benneker et al., 2005). However, none of these parameters take into account the dynamic behavior of the disc. While MRI seemed less accurate than discography to detect radial fissure (Osti and Fraser, 1992), other works suggest that successful detection is possible via MRI (Yu et al., 1988; Saifuddin et al., 1998). A specific relationship was found between radial fissure and positive discography (Aprill and Bogduk, 1992) by delineating a high intensity zone (HIZ) in the annulus fibrosus (AF), as the latter appears hyperintense compared to the nucleus pulposus (NP) on T₂-weighted MR images. However, as HIZ is frequent among pain-free subjects, controversy exists and still persists (Khan et al., 2014).

Contrarily to qualitative T₁- and T₂-weighted MRI, where the signal intensity interpretation can be affected by several experimental factors, biased perception, or scanner variability, quantitative MRI provides an absolute, objective and lower-variability characterization of the imaged anatomical structures. In intervertebral discs, typically quantified parameters are the T₁, T₂, or T_1ρ relaxation times, magnetization transfer (MT) ratio, and the apparent diffusion coefficient (ADC). These parameters can potentially be used as non-invasive biomarkers for different degeneration stages since they are related to the disc histochemical composition and condition, such as for example water, proteoglycan, or collagen content, matrix integrity, loading or aging (Tertti et al., 1991; Watanabe et al., 2007; Mwale et al., 2008; Menezes-Reis et al., 2015; Hwang et al., 2016; Galley et al., 2017; Paul et al., 2018; Enokida et al., 2020). In general, mapping the relaxation time constants provides the advantage of more accurate and reproducible differentiation between the disc tissues of NP and AF and reduced segmentation variability.

To our knowledge, there is no reference available in the literature regarding MRI capability of detecting annulus fissures and characterizing their biomechanical and morphological variations as a function of bending load. We hypothesized that we could non-invasively identify an annular radial fissure and analyze the dynamic behavior of the NP migration by using optimized T₁ and T₂ mapping protocols to assess and quantify the degree of disc fissuring, along with the degree of nucleus displacement/deformation through the fissure, depending on the spine position.

We present here the first results of an ex vivo study. Quantitative and anatomical MRI scans were performed on a loaded fresh lamb spine specimen under 3 different mechanical states (neutral, flexion, extension) in both intact and damaged conditions (anterior radial fissure for the purpose of this Proof-Of-Concept study). Our first objective was to assess whether our MRI protocol allows to detect and further characterize a radial fissure within the annulus. Our second objective was to determine the deformation/displacement fields induced in the nucleus pulposus based on the imposed bending angle of flexion/extension, in both intact and damaged conditions. Our third and last objective was to elaborate a simplified finite element model from the actual geometry and boundary conditions in order to study the contribution of the radial fissure onto the stress distribution in the nucleus.

Materials and Methods

Specimen

We used a fresh lamb spine specimen, comprising 3 vertebrae (L1-L3) and 2 intervertebral discs (2 functional units) of the upper lumbar spine (lamb fillet). The specimen was kept at 4^circC the night before the experiment. Two functional units were kept, such that a healthy control disc could be used at all times for comparison (e.g., in case of aberrant findings).

With the exception of ligamentous and articular capsule tissue, all soft tissue was otherwise removed. In the vertebral body of the 1st and 3rd vertebra, we inserted 2 MR compatible sticks of length L = 8 cm, and diameter ∅ = 6 mm. We kept the specimen hydrated at all times by wrapping it inside saline soaked gauze (Wilke et al., 1998b).

An hour and a half prior to the experiment (Adams, 1995), to avoid over-hydration (McMillan et al., 1996), we creep-loaded the specimen using an axial compressive load of 60 N, resulting in stress of 0.03 MPa. The load was applied by an elastic compressive system, which we set on the transverse vertebral process of the 1st and 3rd vertebra. To simulate an apparent compressive load resulting from muscle activity and gravity (Callaghan et al., 1998), the compressive elastic system was maintained throughout the experiment (Figure 1A).

FIGURE 1

Figure 1. Specimen preparation and boundary conditions of the specimen, in panels (A) neutral position, (B) flexion position, (C) extension position.

First, the specimen was imaged under several mechanical states with intact discs (see below). Secondly, we punctured the anterior part of the annulus using a lumbar puncture needle (18G – 1.2 mm diameter) to damage the disc. The resulting fissure represented 54% of the total height and 27% of the disc transversal diameter. We chose to damage the anterior part of the disc from its inherently easier access in comparison with the posterior part to have a perfectly reproducible lesion model. Besides, similar types of fissure have been documented using discography (Saifuddin et al., 1998) or cadaveric studies (Yu et al., 1988).

MRI Acquisitions

All images were acquired using a 3.0-Tesla imager (Siemens Healthineers, Erlangen, Germany). Room temperature was kept constant at 24°C during the entire examination, and three different MR sequences were performed:

(1) Anatomical images were acquired using a multi-slice 2D T₂-weighted Turbo Spin Echo (T₂-w TSE) sequence with fat suppression, number of averages NA = 10, echo time TE = 99 ms, repetition time TR = 3800 ms, and echo train length = 21. A total of 16 axial slices with resolution 0.31 0.31 1 mm³, and in-plane field of view (FOV) of 84 100 mm² covering the disk and part of the adjacent vertebrae were obtained in 21 min 41 s.

(2) T₂ mapping was performed using a multi-echo Spin Echo sequence with TR = 4500 ms, and NA = 1. An echo train of 32 with first echo at 30 ms and echo spacing of 30 ms was used. A single transverse slice with in-plane resolution of 0.52 0.52 mm² (FOV = 84 99 mm²) and slice thickness of 2 mm was acquired in 12 min.

(3) T₁ mapping was obtained from multiple Inversion Recovery (IR) Turbo-Spin-Echo sequences with TE/TR = 14/6000 ms, NA = 1, and inversion times of 30, 60, 130, 300, 600, 1300, 3000, 5800 ms. All 8 acquisitions were performed across a single transverse slice of resolution 0.52 0.52 mm², with corresponding FOV = 84 99 mm² and slice thickness of 2 mm, acquired in 18 min.

We used the same imaging protocol for each mechanical state of the specimen, resulting in a total of 18 scans. Axial slice orientation was chosen to provide a better visualization and a more appropriate modeling of the NP, the AF and its fissure within the disc with respect to the sagittal view.

Image Processing

Images were processed using in-house scripts and functions developed with MATLAB (MathWorks, Natick, MA, United States). T₁ and T₂ maps were obtained after normalizing the data and fitting them pixel-wise using either an exponential build-up (S₁) or decay (S₂) model:

$$S_1=a\cdot (1-2\cdot e^{t\mathrm{/}{T}_1})$$

$$S_2=b\cdot e^{-t\mathrm{/}{T}_2}$$

With S_i the normalized signal intensity in arbitrary units, t the time obtained from the echo train or inversion times, and a/b constants of the fitting models. Relaxation times are tied to the magnetic properties of the different tissue species and are expected to vary significantly between the different disc compartments.

Mechanical Loading

We controlled the boundary conditions by fixing the lower vertebra and using rubber bands in different configurations depending on the various elastic systems targeted. The first rubber band applied axial compressive loading (as described above), and the second one applied bending loading (either in extension or in flexion). We used this elastic loading system to modify the position of the specimen in the sagittal plane: flexion was created by applying a bending load on the anterior part of the sticks (Figure 1B) and extension on the posterior part (Figure 1C). Our loading protocol is described in Table 1.

TABLE 1

Table 1. Mechanical states and the corresponding estimated boundary conditions of the specimen.

Measurements

We used a mark tracking system to calculate the bending angle and the axial displacement (Figure 1A). In particular, we placed visible marks on the sticks and on the vertebra specimens, and measured their respective positions relative to each other before and during each loading step (Germaneau et al., 2016).

We calculated the bending angle using the following process:

• We determined the baseline angle made by the upper MRI-compatible sticks and the horizontal one (α₀) from the specimen in neutral loading.

• The same measurement was repeated with the bending loading, i.e., flexion or extension (α).

• We subtracted α₀ to α and obtained the corresponding bending angle. A negative value represents a flexion, and a positive value an extension.

We calculated the axial strain from the following process:

• In resting state, we determined the overall height of the vertebral specimen by measuring the distance between the two upper parts of the MRI-compatible sticks (h₀).

• The same measurement was repeated with mechanical loading (h).

• We computed the axial strain as the variation of h relative to the baseline value h₀ at resting state. A negative value represents a height reduction (i.e., a compression) and positive value a height increase.

We compared bending angles and strain values after each loading step. We evaluated the measurement uncertainty (from measurement repetition on known imposed values) at 0.1% for axial strain and 0.1° for angle variation (Germaneau et al., 2016).

Furthermore, we measured nucleus displacement and inferred its strain in the sagittal direction from the T₁ maps. Displacement on sagittal disk boundaries was calculated as the relative difference in distance between the nucleus boundaries at the reference mechanical state (step 1) and the studied state (e.g., step 3, Figure 2).

FIGURE 2

Figure 2. Method to calculate the nucleus sagittal boundary displacement (w), the degree of nucleus pulposus (NP) migration (m) within the annulus fibrosus (AF), and the nucleus sagittal strain for each mechanical step. Panel (A) is the reference mechanical loading state, i.e., step 1, panel (B) is a mechanical state with partial migration of the annulus m < 50% of the total annulus thickness, panel (C) is a mechanical state at total migration of the annulus m > 50% of the total annulus thickness. Nucleus sagittal boundary displacement is calculated for each mechanical state by w – w₀. Nucleus sagittal strain is calculated for each mechanical state by (w – w₀)/w₀. Migration m is calculated from the proportion of nucleus displacement within the annulus.

Total sagittal nucleus displacement was computed as the relative difference between the distances of the nucleus and spinal cord centroids at the reference state (step 1) and each studied step (e.g., step 4). To determine the position of each centroid, we segmented and isolated the nucleus pulposus and the spinal canal from the T₁ maps (Figure 3). The spinal cord is a fixed and identifiable structure that provides a good reference for the centroid tracking process. Nucleus displacement was obtained with an uncertainty of 0.05 mm.

FIGURE 3

Figure 3. Method to calculate global nucleus sagittal displacement. Panel (A) is the reference mechanical loading state, i.e., step 1. Panel (B) is the studied mechanical state, i.e., step 2 to 6 (Table 1). C_np is the centroid of the nucleus pulpous; C_sc is the centroid of the spinal cord; d₀ is the reference distance; d is the distance of the studied mechanical state.

We inferred the nucleus sagittal strain according to:

$$\mathrm{\epsilon }=\left(\text{w}-\text{w}0\right)/\text{w}0$$

where w0 is the distance between anterior and posterior nucleus boundaries at step 1 (Figure 2A, reference state), and w is the distance between anterior and posterior nucleus boundaries at the investigated steps (Figures 2B,C).

We further evaluated the degree of migration of the nucleus in the fissure during the different mechanical loading steps. This migration was calculated as the ratio between the displacement of the nucleus in the fissure, and the annulus thickness (Figure 2).

Finite Element Model

A finite element (FE) model was developed to perform a stress analysis from a simplified geometry extracted using MRI. Our aim was to use this specific FE model from the actual geometry and the actual boundary conditions in order to calculate the stress distribution in the nucleus and the contribution of the radial fissure. FE analysis was carried out using SolidWorks software (Dassault Systèmes Corporation, Waltham, MA, United States). The geometry of the model was defined using points extracted from the boundary of the annulus and nucleus segmented based on a threshold on the MR images acquired in the intact and fissured states (Figure 4). For this stress analysis, homogeneous behavior was considered to model materials (Shin et al., 2007), of which the properties are specified in Table 2.

FIGURE 4

Figure 4. Example of segmentation performed in the intact (top) and fissured (bottom) states, used as input for the finite element model.

TABLE 2

Table 2. Young Modulus and Poisson ratio used for the finite element model of the disc.

The boundary conditions applied to the endplates of the FE model corresponded to the loading imposed during experiments to create bending moments. For that operation, moment was applied on the upper plate of the disk and displacement was considered null on the inferior plate. From there, axial stress and intensity of shear stress induced in the nucleus could be obtained.

Results

MR Imaging Results

Examples of T₂-weighted images, T₁ and T₂ maps are presented in Figures 5, 6 for each mechanical state of the specimen, with or without radial fissure radial fissure, respectively. On Figure 5, the anterior fissure is clearly visible for each MRI sequence.

FIGURE 5

Figure 5. MR images for each mechanical loading state of the intact specimen. First row corresponds to the T₂-weighted images; second row to the computed T₁ maps, with a scale ranging from 0 ms (blue) to 1500 ms (red); and third row to the T₂ maps, with a scale ranging from 20 ms (blue) to 150 ms (red); for three different positions (A) Neutral (step 1) (B) Flexion (step 2), (C) Extension (step 3).

FIGURE 6

Figure 6. MR images for each mechanical loading state of the specimen with radial tear. First row corresponds to the T₂-weighted images; the second row to T₁ maps, with a scale ranging from 0 ms (blue) to 1500 ms (red); and third row to the T₂ maps, with a scale ranging from 20 ms (blue) to 150 ms (red); for three different positions (A) Neutral (step 6), (B) Flexion (step 5), (C) Extension (step 4).

According to the Dallas classification (Sachs et al., 1987), the analyses of the T_2–w TSE, T₁ and T₂ maps reveals a fissure shape similar to a grade II discogram in neutral and extension positions, whereas it is similar to a grade I discogram in flexion position. T₁ and T₂ maps reveal tissue infiltration within the fissure in neutral and in extension positions based on the respective relaxation times measured in the different compartments. This infiltration is no longer visible in flexion position of the specimen, and the magnetic properties of the infiltrated tissue are similar to those of the annulus fibrosus (T₁ > 800 ms, T₂ > 80 ms). These maps confirm that the nucleus deformed in accordance with the loading direction in a greater extent when the annulus is fissured.

Total Disc and Nucleus Mechanical Displacements and Strains

Axial strain, bending angle variation, nucleus boundary sagittal displacement, nucleus total sagittal displacement, and nucleus sagittal strain for each mechanical condition are summarized in Table 3. Steps 1 to 3 (Table 1) present our findings for the intact specimen at different mechanical states. Between steps 3 and 4, an experimental anterior fissure was performed at the annulus fibrosus site. Accordingly, steps 4 to 6 present the results for the damaged specimen.

TABLE 3

Table 3. Axial strain, angle variation, nucleus boundary sagittal displacement, nucleus global sagittal displacement, nucleus sagittal strain, migration of the nucleus in the annulus and cumulative time under axial compression for each mechanical loading step.

Axial strain results show an increase with time when comparing two steps with the same mechanical state, e.g., steps 1 and 6. Specimen angle results also show an increase with time. There is a considerable change in nucleus behavior before and after fissuring the annulus. Displacement of the nuclear boundaries increases with extension after the fissuring (+2.8 mm), whereas posterior displacement decreases with flexion (−1.6 mm). The same behavior is visible for nucleus global displacement (+4.6 mm for extension and −2 mm for flexion). The strongest effect of the fissure is observed on the sagittal strain (+27% in extension and +2.5% in flexion). Both flexion and extension with an intact nucleus lead to nucleus compression (negative strain values of −5.6% and −6.1% respectively). In contrast, with a fissured annulus, only flexion leads to nucleus compression (−3.1%), whereas extension causes its stretching (+21.1%).

Stress Analysis

The maps of axial stress and intensity of maximal shear stress for intact and fissured specimen are shown in Figures 7, 8, for both flexion and extension loading. For the intact specimen, we observed classical behavior with compressive and tensile components in accordance with bending loads. For the fissured specimen, we observed increased stress value in the vicinity of the fissure (around 0.2 MPa). The distribution of the intensity of maximal shear stress further exacerbates this difference by displaying increased values in front of the fissure (around 0.2 MPa), whereas maps were homogeneous in the intact disc.

FIGURE 7

Figure 7. Distribution of the axial stress map for intact (top) and fissured (bottom) disc with flexion (left) and extension (right) loading.

FIGURE 8

Figure 8. Distribution of the intensity of maximal shear stress for intact (top) and fissured (bottom) disc with flexion (left) and extension (right) loading.

Discussion

Original Findings

Our first and secondary aims were to detect radial fissure of the annulus fibrosus using quantitative MRI, and to characterize the biomechanical behavior of the fissured and intact intervertebral disc, respectively. We defined an original protocol combining qualitative and axial quantitative MRI, computed segmentation, mark tracking analysis and FE modeling. This original approach enabled characterizing the morphology and the biomechanics of the radial fissure the nucleus pulposus according to sagittal bending load and the presence / absence of an anterior radial fissure.

By comparison, classical MRI approach assess either the grade of degeneration (Pfirrmann et al., 2001), the inflammatory status of the endplate (Modic et al., 1984), the presence of a high intensity zone in the annulus (Aprill and Bogduk, 1992) or the external geometry of the disc (Fardon et al., 2014). None of these parameters are a direct observation of a radial fissure. In addition, discography allow a morphologic characterization of the fissure but it is invasive (Carragee et al., 2009) and does not explore the dynamic nature of this fissure (Bogduk, 2013).

Data from this proof-of-concept study tend to confirm the dynamic nature of the nucleus pulposus under the influence of bending loads, i.e., deformation/displacement of the nucleus pulposus away from the direction of the load. A systematic review confirmed this behavior for healthy discs, but findings were conflicting regarding pathological discs (Kolber and Hanney, 2009). We can highlight that none of the papers included in this review used a protocol similar to ours.

Our specimen was prepared according to the guidelines for spinal cadaveric studies by controlling for hydration and axial creep load (Adams, 1995; Wilke et al., 1998b). Due to experimental limitations, the axial compressive load was lower than the recommended value (60 N versus 300 N in the guidelines) and was instead applied longer in order to reach the equivalent effect (90 min versus 30 min in the guidelines).

Adjunction of quantitative imaging improves the precision of disc morphometric characterization and enables experimenters to measure both nucleus displacement with good accuracy (0.1 voxel) and strain under various loading conditions (flexion-extension). By comparing the T₁ relaxation times of the tissues at the fissure level to those of the nucleus, it appears reasonable to claim that the nucleus moves toward/deformed into the fissure under bending load influence. Our study showed that T₁ and T₂ values are relevant parameters for automated segmentation of the different disc regions. This justifies future effort in developing fast quantitative acquisitions that will offer objective metrics for the analysis of the disc biomechanics. Indeed, and as indicated by the plethora of studies in the field, there is a clear need of establishing robust quality criteria for image analyses used in clinical research and clinical trials. This task is all the more difficult as conventional MR images display only shades of gray that are very dependent on the hardware (coil profile, field strength, field homogeneity) and the operator (sequence parameters, patient positioning, signal intensity thresholds for segmentation, etc.).

When comparing two steps within the same mechanical state, e.g., steps 2 and 5, an increase of axial strain and bending angle with time was expected. This corresponds to the creep load induced by compressive axial load. However, between step 1 (intact disc in neutral state) and step 6 (fissured disc in neutral state) the angle difference is + 4°, indicating that the specimen’s position shifted from neutral toward extension even though the elastic bending system had been removed. This value is much greater than what we observed for the other two mechanical states (+ 0.7° for extension and − 0.1° for flexion). This observation, albeit not expected, could explain why MRI and mechanical results of step 6 (fissured disc in neutral state) are very close to those of step 4 (fissured disc in extension). This sizable increase of extension angle in step 6 could be a consequence of the manually performed anterior fissure and should be further explored in future studies.

Our results also revealed a link between disc abnormality and range of motion by using a non-invasive diagnostic method differing from interventional intra-discal procedures. While discography highlights the physical presence of a fissure within the annulus under fixed positional conditions, it remains impossible to observe the migration of the nucleus under postural influence (Walker et al., 2008). Using our non-invasive MRI approach, we observed and quantified the migration of the nucleus into the annulus from the quantitative maps, under conditions equivalent to the physiological loading (up to 5-6 degrees per level) usually observed for daily range of motion in humans (Alini et al., 2008). In this work, the artificial fissure induced considerable damage including the migration of the nucleus within the full thickness of the annulus (100% of migration). As a result, the nucleus strain increased with stretching in the sagittal direction inducing reduction of internal pressure. Yet, flexion partially restored not only the nucleus morphology with observed migration dropping from 100 to 21%, but also the internal pressure with a measured compression strain of 3.1% restoring internal pressure, and consequently physiological capabilities of the disc. These observations confirm the dynamic behavior changes of the fissure.

Ultimately, we used a simplified finite element model, based on geometry directly extracted from MRI data. This approach was enriched by the real boundary conditions quantified during the experiments. This dual approach enabled simultaneous measurement of stress and strain fields on the nucleus, thereby providing a complete biomechanical analysis of the normal and pathological disc. From a clinical point of view, the migration of the nucleus pulposus alongside the fissure could become a nociceptive trigger, leading patient pain. Furthermore, in this case, the singularity and its evolution observed from MRI maps can be linked to stress variation in the nucleus with an increasing of shear components.

Study Limitations

Despite being encouraging, our results suffer from several limitations. This study is a proof-of-concept study performed on only one disc. Repeating the protocol on multiple samples would allow us to assess the reproducibility of our method and to draw more robust conclusions from a statistical point of view. Another limitation could arise from the animal nature of our specimen and its difference with human discs. Indeed, quadrupedal station induces geometrical changes (O’Connell et al., 2007), variations in mechanical properties (Alini et al., 2008) and chemical composition of the disc (Zhang et al., 2014). Nonetheless, a recent review (Daly et al., 2016) concluded that ovine disc is a reasonable choice for preliminary biomechanical or injury model studies. Animal models are currently used in surgical, biomechanical (Casaroli et al., 2017) and histological studies (Schollum et al., 2008). Rather than species type, age and degeneration state of the specimen might actually have a stronger impact on the results since both of these parameters alter disc biomechanics (Adams et al., 1996). With age and degeneration, there is a decrease in the disc notochordal cell population. This triggers a reduction of proteoglycan secretion, leading in turn to a drop in nuclear hydration. Type II collagen is progressively replaced by type I collagen, leading to a more fibrous nucleus (Adams and Roughley, 2006). All of these processes impact the mechanical behavior of the nucleus according to motion and can affect the magnetic properties, such as T₁ and T₂. As our specimen comes from a young animal with no sign of degenerative disc, the results could not be generalized beyond these criteria. Another limitation comes from the type of fissure induced in our study. Our experimental setup ensured good access to the anterior part of the disc without damaging the specimen. Though rare, anterior fissures are encountered in discography (Saifuddin et al., 1998) or cadaveric studies (Yu et al., 1988). However, most of the fissures associated with discogenic backpain are either posterior or postero-lateral (Bogduk et al., 2013). Because of the different shape and thickness of the anterior annulus (Cassidy et al., 1989), results could be different with a posterior fissure. Ultimately, MRI examination time needs to be shortened to envision in vivo studies. Since ours was the first study of its kind, basic MR sequences were used in order to assess the feasibility of T₁ and T₂ maps as markers to investigate the disc biomechanics. As a consequence, repeating the current imaging protocol in the 3 different conditions leads to a 6-h examination. Such scan time did not have an impact on our specimen since it has been shown that no change in mechanical properties is observed after 20 h of testing (Wilke et al., 1998a), but faster acquisition schemes need to be implemented to allow a potential transfer to patients in vivo. We propose to implement strategies that enable simultaneous acquisitions of T₁ and T₂ (Ma et al., 2013) in future studies.

Concerning FE analysis, we used a simplified model to determine stress distribution according to imposed motion and the presence of a fracture in the annulus. Even if we used real geometry and boundary conditions read during experiments, our model had some limitations. As already performed for previous works from literature, we considered linear behaviors and isotropic properties (Shin et al., 2007; Zanjani-Pour et al., 2016). This approach was sufficient in the present work to analyze stress distribution. However, to validate the model, it would be necessary to identify the mechanical properties of tissues and to implement them in the model which could be validated from comparison of displacement and strain fields. Furthermore, we did not model the facet joints in accordance with the limited values of imposed moments (3 N.m). According to our experiments, facet joints were not activated for this magnitude of applied loads, however, they would have to be modeled for higher imposed moments.

Future directions would include replicating this study using a greater number of ovine discs, along with inducing different orientations for the radial fissure in order to better reflect clinical conditions, and applying additional movement directions, typically including the coronal and frontal plans. If results are consistent, our ambition would be to design the next phase of ex vivo studies using human cadaveric discs at different levels of degeneration and age, and thereby make it possible to transfer our MRI methods to in vivo clinical applications.

Future Clinical Implications

Once adapted to in vivo experimentations and combined with classical assessment tools (MRI, X-ray, etc.), such protocol would help clinicians to assess patients with discogenic pain and/or radicular pain due to disc herniation. Together with the directional preference concept (May and Aina, 2012; May et al., 2018), it would help to identify patients who could benefit from a physiotherapeutic mechanical treatment (McKenzie and May, 2003) from those who need surgery. We believe that this protocol would ease the clinical decision-making process, and hence optimize patient care as well as reduce health cost related to back pain.

Conclusion

In this Proof-of-Concept study, we demonstrated the possibility to characterize the morphological and biomechanical parameters of a radial fissure within an ex vivo ovine disc. To do so, we combined quantitative T₁ an T₂ mapping MRI, T₂-weighted MRI, computed segmentation, mark tracking analysis and finite element modeling.

Each MRI sequence allowed a clear and original visualization of the discal damage, i.e., radial fissure, non-invasively. The nucleus pulposus moved anteriorly with extension bending load and posteriorly with flexion bending load. After the anterior annulus was damaged, the nuclear displacement and strain increased for extension load and decreased for flexion load. The displacement and strain of the nucleus appeared to follow the direction of the fissure and the direction of the bending load.

This preliminary work, once validated on a larger scale, could have substantial applications for:

• radio-clinical non-invasive disc explorations and correlations with patient lumbar pain,

• dynamic characterization of the disc under physiological and pathological conditions.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

J-PD, TV, and PR: conceptualization. AG, MY, MS, and NS: methodology. J-PD, AG, MY, MS, and NS: software and validation. MS and NS: resources. J-PD, AG, NS, and PR: writing—original draft preparation. J-PD, TV, AG, MB, MR, NS, and PR: writing—review and editing. AG, NS, and PR: supervision and project administration. NS, MS, and PR: funding acquisition. All authors have read and agreed to the published version of the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Adams, M. A. (1995). Spine Update: mechanical testing of the spine - an appraisal of methodology. results and conclusion. Spine 20, 2151–2156.

Google Scholar

Adams, M. A., and Dolan, P. (2012). Intervertebral disc degeneration: evidence for two distinct phenotypes. J. Anat. 221, 497–506. doi: 10.1111/j.1469-7580.2012.01551.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Adams, M. A., McNally, D. S., and Dolan, P. (1996). “Stress” distributions inside intervertebral discs. the effects of age and degeneration. J. Bone Joint Surg. Br. 78, 965–972. doi: 10.1302/0301-620x78b6.1287

PubMed Abstract | CrossRef Full Text | Google Scholar

Adams, M. A., and Roughley, P. J. (2006). What is intervertebral disc degeneration, and what causes it? Spine 31, 2151–2161. doi: 10.1097/01.brs.0000231761.73859.2c

CrossRef Full Text | Google Scholar

Alini, M., Eisenstein, S. M., Ito, K., Little, C., Kettler, A. A., Masuda, K., et al. (2008). Are animal models useful for studying human disc disorders/degeneration? Eur. Spine J. 17, 2–19. doi: 10.1007/s00586-007-0414-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Antoniou, J., Steffen, T., Nelson, F., Winterbottom, N., Hollander, A. P., Poole, R. A., et al. (1996). The human lumbar intervertebral disc: evidence for changes in the biosynthesis and denaturation of the extracellular matrix with growth, maturation, ageing, and degeneration. J. Clin. Invest. 98, 996–1003. doi: 10.1172/JCI118884

PubMed Abstract | CrossRef Full Text | Google Scholar

Aprill, C., and Bogduk, N. (1992). High-intensity zone: a diagnostic sign of painful lumbar disc on magnetic resonance imaging. Br. J. Radiol. 65, 361–369. doi: 10.1259/0007-1285-65-773-361

PubMed Abstract | CrossRef Full Text | Google Scholar

Benneker, L. M., Heini, P. F., Anderson, S. E., Alini, M., and Ito, K. (2005). Correlation of radiographic and MRI parameters to morphological and biochemical assessment of intervertebral disc degeneration. Eur. Spine J. 14, 27–35. doi: 10.1007/s00586-004-0759-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bogduk, N. (2013). Practice Guidelines for Spinal Diagnostic and Treatment Procedures, 2nd Edn. San Francisco, CA: International Spine Intervention Society.

Google Scholar

Bogduk, N., Aprill, C., and Derby, R. (2013). Lumbar discogenic pain: state-of-the-art review: lumbar discogenic pain. Pain Med. 14, 813–836. doi: 10.1111/pme.12082

PubMed Abstract | CrossRef Full Text | Google Scholar

Brinjikji, W., Luetmer, P. H., Comstock, B., Bresnahan, B. W., Chen, L. E., Deyo, R. A., et al. (2015). Systematic literature review of imaging features of spinal degeneration in asymptomatic populations. Am. J. Neuroradiol. 36, 811–816. doi: 10.3174/ajnr.A4173

PubMed Abstract | CrossRef Full Text | Google Scholar

Callaghan, J. P., Gunning, J. L., and McGill, S. M. (1998). The relationship between lumbar spine load and muscle activity during extensor exercises. Phys. Ther. 78, 8–18. doi: 10.1093/ptj/78.1.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Carragee, E. J., Don, A. S., Hurwitz, E. L., Cuellar, J. M., Carrino, J., and Herzog, R. (2009). 2009 ISSLS prize winner: does discography cause accelerated progression of degeneration changes in the lumbar disc: a ten-year matched cohort study. Spine 34, 2338–2345. doi: 10.1097/BRS.0b013e3181ab5432

PubMed Abstract | CrossRef Full Text | Google Scholar

Casaroli, G., Galbusera, F., Jonas, R., Schlager, B., Wilke, H.-J., and Villa, T. (2017). A novel finite element model of the ovine lumbar intervertebral disc with anisotropic hyperelastic material properties. PloS One 12:e0177088. doi: 10.1371/journal.pone.0177088

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassidy, J. J., Hiltner, A., and Baer, E. (1989). Hierarchical structure of the intervertebral disc. Connect. Tissue Res. 23, 75–88.

Google Scholar

Coppes, M. H., Marani, E., Thomeer, R. T. W. M., and Groen, G. J. (1997). Innervation of “painful”. Lumbar Discs. Spine 22, 2342–2349. doi: 10.1097/00007632-199710150-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Coppes, M. H., Marani, E., Thomeer, R. T. W. M., Oudega, M., and Groen, G. J. (1990). Innervation of annulus fibrosis in low back pain. Lancet 336, 189–190. doi: 10.1016/0140-6736(90)91723-n

CrossRef Full Text | Google Scholar

Dagenais, S., Caro, J., and Haldeman, S. (2008). A systematic review of low back pain cost of illness studies in the United States and internationally. Spine J. 8, 8–20. doi: 10.1016/j.spinee.2007.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Daly, C., Ghosh, P., Jenkin, G., Oehme, D., and Goldschlager, T. (2016). A review of animal models of intervertebral disc degeneration: pathophysiology, regeneration, and translation to the clinic. BioMed. Res. Int. 2016:5952165. doi: 10.1155/2016/5952165

PubMed Abstract | CrossRef Full Text | Google Scholar

Enokida, S., Tanishima, S., Tanida, A., Mihara, T., Takeda, C., Yamashita, E., et al. (2020). Evaluation of age-related changes in lumbar facet joints using T2 mapping. J. Orthop. Sci. 25, 46–51. doi: 10.1016/j.jos.2019.02.017 ^∗y2019,

PubMed Abstract | Google Scholar

Fardon, D. F., Williams, A. L., Dohring, E. J., Murtagh, F. R., Gabriel Rothman, S. L., and Sze, G. K. (2014). Lumbar disc nomenclature: version 2.0. Spine J. 14, 2525–2545. doi: 10.1016/j.spinee.2014.04.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Freeman, B. J. C. (2006). IDET: a critical appraisal of the evidence. Eur. Spine J. 15, 448–457. doi: 10.1007/s00586-006-0156-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Galley, J., Maestretti, G., Koch, G., and Hoogewoud, H.-M. (2017). Real T1 relaxation time measurement and diurnal variation analysis of intervertebral discs in a healthy population of 50 volunteers. Eur. J. Radiol. 87, 13–19. doi: 10.1016/j.ejrad.2016.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Cosamalón, J., Del Valle, M. E., Calavia, M. G., García-Suárez, O., López-Muñiz, A., Otero, J., et al. (2010). Intervertebral disc, sensory nerves and neurotrophins: who is who in discogenic pain: intervertebral disc, sensory nerves and neurotrophins. J. Anat. 217, 1–15. doi: 10.1111/j.1469-7580.2010.01227.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Germaneau, A., Vendeuvre, T., Saget, M., Doumalin, P., Dupré, J. C., Brémand, F., et al. (2016). A novel approach for biomechanical spine analysis: mechanical response of vertebral bone augmentation by kyphoplasty to stabilise thoracolumbar burst fractures. J. Mech. Behav. Biomed. Mater. 59, 291–303. doi: 10.1016/j.jmbbm.2016.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Hashemi, M., Dadkhah, P., Taheri, M., Katibeh, P., and Asadi, S. (2020). Effectiveness of intradiscal injection of radiopaque gelified ethanol (DiscoGel ^®) versus percutaneous laser disc decompression in patients with chronic radicular low back pain. Korean J. Pain 33, 66–72. doi: 10.3344/kjp.2020.33.1.66

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, D., Kim, S., Abeydeera, N. A., Statum, S., Masuda, K., Chung, C. B., et al. (2016). Quantitative magnetic resonance imaging of the lumbar intervertebral discs. Quant. Imaging Med. Surg. 6, 744–755. doi: 10.21037/qims.2016.12.09

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishibashi, K., Iwai, H., and Koga, H. (2019). Chemonucleolysis with chondroitin sulfate ABC endolyase as a novel minimally invasive treatment for patients with lumbar intervertebral disc herniation. J. Spine Surg. 5, S115–S121. doi: 10.21037/jss.2019.04.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Javid, M. J., Nordby, E. J., Ford, L. T., Hejna, W. J., Whisler, W. W., Burton, C., et al. (1983). Safety and efficacy of chymopapain (Chymodiactin) in herniated nucleus pulposus with sciatica. results of a randomized, double-blind study. JAMA 249, 2489–2494.

Google Scholar

Khan, I., Hargunani, R., and Saifuddin, A. (2014). The lumbar high-intensity zone: 20 years on. Clin. Radiol. 69, 551–558. doi: 10.1016/j.crad.2013.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Knezevic, N. N., Mandalia, S., Raasch, J., Knezevic, I., and Candido, K. D. (2017). Treatment of chronic low back pain – new approaches on the horizon. J. Pain Res. 10, 1111–1123. doi: 10.2147/JPR.S132769

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolber, M. J., and Hanney, W. J. (2009). The dynamic disc model: a systematic review of the literature. Phys. Ther. Rev. 14, 181–189. doi: 10.1179/174328809X452827

CrossRef Full Text | Google Scholar

Lama, P., Le Maitre, C. L., Harding, I. J., Dolan, P., and Adams, M. A. (2018). Nerves and blood vessels in degenerated intervertebral discs are confined to physically disrupted tissue. J. Anat. 233, 86–97. doi: 10.1111/joa.12817

PubMed Abstract | CrossRef Full Text | Google Scholar

Laslett, M., Öberg, B., Aprill, C. N., and McDonald, B. (2005). Centralization as a predictor of provocation discography results in chronic low back pain, and the influence of disability and distress on diagnostic power. Spine J. 5, 370–380. doi: 10.1016/j.spinee.2004.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Maitre, C. L., Freemont, A. J., and Hoyland, J. A. (2004). Localization of degradative enzymes and their inhibitors in the degenerate human intervertebral disc. J. Pathol. 204, 47–54. doi: 10.1002/path.1608

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Maitre, C. L., Pockert, A., Buttle, D. J., Freemont, A. J., and Hoyland, J. A. (2007). Matrix synthesis and degradation in human intervertebral disc degeneration. Biochem. Soc. Trans. 35, 652–655. doi: 10.1042/BST0350652

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P., Zhang, R., and Zhou, Q. (2017). Efficacy of platelet-rich plasma in retarding intervertebral disc degeneration: a meta-analysis of animal studies. BioMed. Res. Int. 2017:7919201. doi: 10.1155/2017/7919201

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, D., Gulani, V., Seiberlich, N., Liu, K., Sunshine, J. L., Duerk, J. L., et al. (2013). Magnetic resonance fingerprinting. Nature 495, 187–192. doi: 10.1038/nature11971

PubMed Abstract | CrossRef Full Text | Google Scholar

Maetzel, A., and Li, L. (2002). The economic burden of low back pain: a review of studies published between 1996 and 2001. Best Pract. Res. Clin. Rheumatol. 16, 23–30. doi: 10.1053/berh.2001.0204

PubMed Abstract | CrossRef Full Text | Google Scholar

Manchikanti, L., Soin, A., and Beyamin, R. (2018). An update of the systematic appraisal of theaccuracy and utility of discography in chronicspinal pain. Pain Physician 1, 91–110. doi: 10.36076/ppj.2018.2.91

CrossRef Full Text | Google Scholar

May, S., and Aina, A. (2012). Centralization and directional preference: a systematic review. Man. Ther. 17, 497–506. doi: 10.1016/j.math.2012.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

May, S., Runge, N., and Aina, A. (2018). Centralization and directional preference: an updated systematic review with synthesis of previous evidence. Musculoskelet. Sci. Pract. 38, 53–62. doi: 10.1016/j.msksp.2018.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

McCormick, Z. L., Lehman, V. T., Plastaras, C. T., Walega, D. R., Huddleston, P., Moussallem, C., et al. (2019). Low-pressure lumbar provocation discography according to spine intervention society/international association for the study of pain standards does not cause acceleration of disc degeneration in patients with symptomatic low back pain: a 7-year matched cohort study. Spine 44, E1161–E1168. doi: 10.1097/BRS.0000000000003085

PubMed Abstract | CrossRef Full Text | Google Scholar

McKenzie, R. (1981). The Lumbar Spine. First. Wellington: Spinal Publication.

Google Scholar

McKenzie, R., and May, S. (2003). The Lumbar Spine: Mechanical Diagnosis and Therapy–Volume 1. New Zealand, OC: Spinal Publications.

Google Scholar

McMillan, D. W., Garbutt, G., and Adams, M. A. (1996). Effect of sustained loading on the water content of intervertebral discs: implications for disc metabolism. Ann. Rheum. Dis. 55, 880–887. doi: 10.1136/ard.55.12.880

PubMed Abstract | CrossRef Full Text | Google Scholar

McNally, D. S., Shackleford, I. M., Goodship, A. E., and Mulholland, R. C. (1996). In vivo stress measurement can predict pain on discography. Spine 21, 2580–2587.

Google Scholar

Meisel, H.-J., Agarwal, N., Hsieh, P. C., Skelly, A., Park, J.-B., Brodke, D., et al. (2019). Cell therapy for treatment of intervertebral disc degeneration: a systematic review. Global Spine J. 9, 39S–52S. doi: 10.1177/2192568219829024

PubMed Abstract | CrossRef Full Text | Google Scholar

Menezes-Reis, R., Salmon, C. E. G., Carvalho, C. S., Bonugli, G. P., Chung, C. B., and Nogueira-Barbosa, M. H. (2015). T1rho and T2 mapping of the intervertebral disk: comparison of different methods of segmentation. AJNR Am. J. Neuroradiol. 36, 606–611. doi: 10.3174/ajnr.A4125

PubMed Abstract | CrossRef Full Text | Google Scholar

Modic, M. T., Pavlicek, W., Weinstein, M. A., Boumphrey, F., Ngo, F., Hardy, R., et al. (1984). Magnetic resonance imaging of intervertebral disk disease. clinical and pulse sequence considerations. Radiology 152, 103–111. doi: 10.1148/radiology.152.1.6729099

PubMed Abstract | CrossRef Full Text | Google Scholar

Moneta, G. B., Videman, T., Kaivanto, K., Aprill, C. N., Spivey, M., Vanharanta, H., et al. (1994). Reported pain during lumbar discography as a function of anular ruptures and disc degeneration: a re-analysis of 833 discograms. Spine 19, 1968–1974.

Google Scholar

Mwale, F., Iatridis, J. C., and Antoniou, J. (2008). Quantitative MRI as a diagnostic tool of intervertebral disc matrix composition and integrity. Eur. Spine J. 17, 432–440. doi: 10.1007/s00586-008-0744-4

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Connell, G. D., Vresilovic, E. J., and Elliott, D. M. (2007). Comparison of animals used in disc research to human lumbar disc geometry. Spine 32, 328–333. doi: 10.1097/01.brs.0000253961.40910.c1

CrossRef Full Text | Google Scholar

Osti, O., and Fraser, R. (1992). MRI and discography of annular tears and intervertebral disc degeneration. a prospective clinical comparison. J. Bone Joint Surg. Br. 74, 431–435. doi: 10.1302/0301-620X.74B3.1587896

PubMed Abstract | CrossRef Full Text | Google Scholar

Paradiso, R., and Alexandre, A. (2005). The different outcomes of patients with disc herniation treated either by microdiscectomy, or by intradiscal ozone injection. Acta Neurochir. Suppl. 92, 139–142. doi: 10.1007/3-211-27458-8_30

CrossRef Full Text | Google Scholar

Paul, C. P. L., Smit, T. H., de Graaf, M., Holewijn, R. M., Bisschop, A., van de Ven, P. M., et al. (2018). Quantitative MRI in early intervertebral disc degeneration: T1rho correlates better than T2 and ADC with biomechanics, histology and matrix content. PloS One 13:e0191442. doi: 10.1371/journal.pone.0191442

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfirrmann, C. W. A., Metzdorf, A., Zanetti, M., Hodler, J., and Boos, N. (2001). Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine 26, 1873–1878. doi: 10.1097/00007632-200109010-00011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sachs, B. L., Vanharanta, H., Spivey, M. A., Guyer, R. D., Videman, T., Rashbaum, R. F., et al. (1987). Dallas discogram description - a new classification of CT/discography in low-back disorders. Spine 12, 287–294.

Google Scholar

Saifuddin, A., Emanuel, R., White, J., Renton, P., Braithwaite, I., and Taylor, B. A. (1998). An analysis of radiating pain at lumbar discography. Eur. Spine J. 7, 358–362. doi: 10.1007/s005860050090

PubMed Abstract | CrossRef Full Text | Google Scholar

Schollum, M. L., Robertson, P. A., and Broom, N. D. (2008). ISSLS prize winner: microstructure and mechanical disruption of the lumbar disc annulus: part I: a microscopic investigation of the translamellar bridging network. Spine 33, 2702–2710. doi: 10.1097/BRS.0b013e31817bb92c

PubMed Abstract | CrossRef Full Text | Google Scholar

Shin, D. S., Lee, K., and Kim, D. (2007). Biomechanical study of lumbar spine with dynamic stabilization device using finite element method. Comput.-Aided Des. 39, 559–567. doi: 10.1016/j.cad.2007.03.005

CrossRef Full Text | Google Scholar

Stefanakis, M., Luo, J., Pollintine, P., Dolan, P., and Adams, M. A. (2014). ISSLS Prize winner: mechanical influences in progressive intervertebral disc degeneration. Spine 39, 1365–1372. doi: 10.1097/BRS.0000000000000389

PubMed Abstract | CrossRef Full Text | Google Scholar

Tertti, M., Paajanen, H., Laato, M., Aho, H., Komu, M., and Kormano, M. (1991). Disc degeneration in magnetic resonance imaging. a comparative biochemical, histologic, and radiologic study in cadaver spines. Spine 16, 629–634. doi: 10.1097/00007632-199106000-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, K. J., Dagher, A. P., Eckel, T. S., Clark, M., and Reinig, J. W. (2009). Modic changes on MR images as studied with provocative diskography: clinical relevance—a retrospective study of 2457 disks. Radiology 250, 849–855. doi: 10.1148/radiol.2503080474

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanharanta, H., Sachs, B. L., Spivey, M., Guyer, R. D., Hochschuler, H., Rashbaum, R. F., et al. (1987). The Relashionship of pain provocation to lumbar disc deterioration as seen by CT/Discography. Spine 12, 295–298.

Google Scholar

Vos, T., Allen, C., Arora, M., Barber, R. M., Bhutta, Z. A., Brown, A., et al. (2016). Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet 388, 1545–1602. doi: 10.1016/S0140-6736(16)31678-6

CrossRef Full Text | Google Scholar

Walker, B. F., Muller, R., and Grant, W. D. (2003). Low back pain in australian adults: the economic burden. Asia Pac. J. Public Health 15, 79–87. doi: 10.1177/101053950301500202

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, J., El Abd, O., Isaac, Z., and Muzin, S. (2008). Discography in practice: a clinical and historical review. Curr. Rev. Musculoskelet. Med. 1, 69–83. doi: 10.1007/s12178-007-9009-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, A., Benneker, L. M., Boesch, C., Watanabe, T., Obata, T., and Anderson, S. E. (2007). Classification of intervertebral disk degeneration with axial T2 mapping. Am. J. Roentgenol. 189, 936–942. doi: 10.2214/AJR.07.2142

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilke, H.-J., Jungkunz, B., Wenger, K., and Claes, L. E. (1998a). Spinal segment range of motion as a function of in vitro test conditions: effects of exposure period, accumulated cycles, angular-deformation rate, and moisture condition. Anat. Rec. 251, 15–19.

Google Scholar

Wilke, H.-J., Wenger, K., and Claes, L. (1998b). Testing criteria for spinal implants: recommendations for the standardization of in vitro stability testing of spinal implants. Eur. Spine J. 7, 148–154. doi: 10.1007/s005860050045

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, S. W., Sether, L. A., Ho, P. S., Wagner, M., and Haughton, V. M. (1988). Tears of the anulus fibrosus: correlation between MR and pathologic findings in cadavers. AJNR Am. J. Neuroradiol. 9, 367–370.

Google Scholar

Zanjani-Pour, S., Winlove, C. P., Smith, C. W., and Meakin, J. R. (2016). Image driven subject-specific finite element models of spinal biomechanics. J. Biomech. 49, 919–925. doi: 10.1016/j.jbiomech.2016.02.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Lenart, B. A., Lee, J. K., Chen, D., Shi, P., Ren, J., et al. (2014). Histological features of endplates of the mammalian spine: from mice to men. Spine 39, E312–E317. doi: 10.1097/BRS.0000000000000174

PubMed Abstract | CrossRef Full Text | Google Scholar