Abstract
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 (T2-weighted, T1 and T2 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 (T1, T2) 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 T1- and T2-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 T2-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 T2-weighted MR images. However, as HIZ is frequent among pain-free subjects, controversy exists and still persists (Khan et al., 2014).
Contrarily to qualitative T1- and T2-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 T1, T2, or T1ρ 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 T1 and T2 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 4circC 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
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 T2-weighted Turbo Spin Echo (T2-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 mm3, and in-plane field of view (FOV) of 84 100 mm2 covering the disk and part of the adjacent vertebrae were obtained in 21 min 41 s.
- (2)
T2 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 mm2 (FOV = 84 99 mm2) and slice thickness of 2 mm was acquired in 12 min.
- (3)
T1 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 mm2, with corresponding FOV = 84 99 mm2 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). T1 and T2 maps were obtained after normalizing the data and fitting them pixel-wise using either an exponential build-up (S1) or decay (S2) model:
With Si 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
| Step | Mechanical loading state | Axial compression (N) | Bending force (N) | Bending moment (N.m) |
| 0 | Neutral | 0 | 0 | 0 |
| 1 | Neutral (Figure 1A) | 59.7 | 0 | 0 |
| 2 | Flexion (Figure 1B) | 59.7 | 81.1 | 2.1 |
| 3 | Extension (Figure 1C) | 59.7 | 80.9 | 2.8 |
| Fissure | ||||
| 4 | Extension (Figure 1C) | 59.7 | 80.5 | 2.8 |
| 5 | Flexion (Figure 1B) | 59.7 | 81.0 | 2.1 |
| 6 | Neutral (Figure 1A) | 59.7 | 0 | 0 |
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 (α0) from the specimen in neutral loading.
- •
The same measurement was repeated with the bending loading, i.e., flexion or extension (α).
- •
We subtracted α0 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 (h0).
- •
The same measurement was repeated with mechanical loading (h).
- •
We computed the axial strain as the variation of h relative to the baseline value h0 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 T1 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
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 T1 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
We inferred the nucleus sagittal strain according to:
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
TABLE 2
| Young modulus (MPa) | Poisson ratio | Reference | |
| Nucleus | 1 | 0.499 | Shin et al., 2007 |
| Annulus | 8.4 | 0.45 |
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 T2-weighted images, T1 and T2 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 6
According to the Dallas classification (Sachs et al., 1987), the analyses of the T2–w TSE, T1 and T2 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. T1 and T2 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 (T1 > 800 ms, T2 > 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
| Step | Axial strain1 | Specimen angle2 | Angle variation3 | Nucleus boundary sagittal displacement4 | Nucleus global sagittal displacement4 | Nucleus sagittal strain1 | Migration of nucleus in the annulus5 | Cumulative time under axial compression |
| 1 (Neutral) | −4 | 4.7 | – | – | – | – | – | lh30 |
| 2 (Flexion) | −7.3 | −8.8 | −13.5 | −2.2 | −3.6 | −5.6 | 0 | 3h |
| 3 (Extension) | −8.3 | 10.4 | 5.7 | 2.2 | −0.9 | −6.1 | 0 | 4h30 |
| Fissure | ||||||||
| 4 (Extension) | −8.8 | 11.1 | 6.4 | 5.0 | 3.7 | 21.1 | 100 | 6h30 |
| 5 (Flexion) | −7.5 | −8.9 | −13.6 | −0.6 | −1.6 | −3.1 | 21 | 8h30 |
| 6 (Neutral) | −4.6 | 8.7 | 4.0 | 4.5 | 1.5 | 20.8 | 100 | 10h |
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.
1In percentage, negative value corresponds to a decrease of the height of the specimen. 2In degrees, negative value corresponds to a flexion angle, positive value to an extension one. 3In degrees, negative value corresponds to a flexion angle, positive value to an extension one. This angle corresponds to the actual specimen angle respective to the baseline’s one (step 1 neutral). 4In millimeters, negative value corresponds to a posterior displacement, positive value to an anterior one. 5In percentage.
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 8
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 T1 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 T1 and T2 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 T1 and T2. 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 T1 and T2 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 T1 and T2 (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 T1 an T2 mapping MRI, T2-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.
Statements
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
1
AdamsM. A. (1995). Spine Update: mechanical testing of the spine - an appraisal of methodology. results and conclusion.Spine202151–2156.
2
AdamsM. A.DolanP. (2012). Intervertebral disc degeneration: evidence for two distinct phenotypes.J. Anat.221497–506. 10.1111/j.1469-7580.2012.01551.x
3
AdamsM. A.McNallyD. S.DolanP. (1996). “Stress” distributions inside intervertebral discs. the effects of age and degeneration.J. Bone Joint Surg. Br.78965–972. 10.1302/0301-620x78b6.1287
4
AdamsM. A.RoughleyP. J. (2006). What is intervertebral disc degeneration, and what causes it?Spine312151–2161. 10.1097/01.brs.0000231761.73859.2c
5
AliniM.EisensteinS. M.ItoK.LittleC.KettlerA. A.MasudaK.et al (2008). Are animal models useful for studying human disc disorders/degeneration?Eur. Spine J.172–19. 10.1007/s00586-007-0414-y
6
AntoniouJ.SteffenT.NelsonF.WinterbottomN.HollanderA. P.PooleR. 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.98996–1003. 10.1172/JCI118884
7
AprillC.BogdukN. (1992). High-intensity zone: a diagnostic sign of painful lumbar disc on magnetic resonance imaging.Br. J. Radiol.65361–369. 10.1259/0007-1285-65-773-361
8
BennekerL. M.HeiniP. F.AndersonS. E.AliniM.ItoK. (2005). Correlation of radiographic and MRI parameters to morphological and biochemical assessment of intervertebral disc degeneration.Eur. Spine J.1427–35. 10.1007/s00586-004-0759-4
9
BogdukN. (2013). Practice Guidelines for Spinal Diagnostic and Treatment Procedures, 2nd Edn. San Francisco, CA: International Spine Intervention Society.
10
BogdukN.AprillC.DerbyR. (2013). Lumbar discogenic pain: state-of-the-art review: lumbar discogenic pain.Pain Med.14813–836. 10.1111/pme.12082
11
BrinjikjiW.LuetmerP. H.ComstockB.BresnahanB. W.ChenL. E.DeyoR. A.et al (2015). Systematic literature review of imaging features of spinal degeneration in asymptomatic populations.Am. J. Neuroradiol.36811–816. 10.3174/ajnr.A4173
12
CallaghanJ. P.GunningJ. L.McGillS. M. (1998). The relationship between lumbar spine load and muscle activity during extensor exercises.Phys. Ther.788–18. 10.1093/ptj/78.1.8
13
CarrageeE. J.DonA. S.HurwitzE. L.CuellarJ. M.CarrinoJ.HerzogR. (2009). 2009 ISSLS prize winner: does discography cause accelerated progression of degeneration changes in the lumbar disc: a ten-year matched cohort study.Spine342338–2345. 10.1097/BRS.0b013e3181ab5432
14
CasaroliG.GalbuseraF.JonasR.SchlagerB.WilkeH.-J.VillaT. (2017). A novel finite element model of the ovine lumbar intervertebral disc with anisotropic hyperelastic material properties.PloS One12:e0177088. 10.1371/journal.pone.0177088
15
CassidyJ. J.HiltnerA.BaerE. (1989). Hierarchical structure of the intervertebral disc.Connect. Tissue Res.2375–88.
16
CoppesM. H.MaraniE.ThomeerR. T. W. M.GroenG. J. (1997). Innervation of “painful”. Lumbar Discs.Spine222342–2349. 10.1097/00007632-199710150-00005
17
CoppesM. H.MaraniE.ThomeerR. T. W. M.OudegaM.GroenG. J. (1990). Innervation of annulus fibrosis in low back pain.Lancet336189–190. 10.1016/0140-6736(90)91723-n
18
DagenaisS.CaroJ.HaldemanS. (2008). A systematic review of low back pain cost of illness studies in the United States and internationally.Spine J.88–20. 10.1016/j.spinee.2007.10.005
19
DalyC.GhoshP.JenkinG.OehmeD.GoldschlagerT. (2016). A review of animal models of intervertebral disc degeneration: pathophysiology, regeneration, and translation to the clinic.BioMed. Res. Int.2016:5952165. 10.1155/2016/5952165
20
EnokidaS.TanishimaS.TanidaA.MiharaT.TakedaC.YamashitaE.et al (2020). Evaluation of age-related changes in lumbar facet joints using T2 mapping.J. Orthop. Sci.2546–51. 10.1016/j.jos.2019.02.017∗y2019,
21
FardonD. F.WilliamsA. L.DohringE. J.MurtaghF. R.Gabriel RothmanS. L.SzeG. K. (2014). Lumbar disc nomenclature: version 2.0.Spine J.142525–2545. 10.1016/j.spinee.2014.04.022
22
FreemanB. J. C. (2006). IDET: a critical appraisal of the evidence.Eur. Spine J.15448–457. 10.1007/s00586-006-0156-2
23
GalleyJ.MaestrettiG.KochG.HoogewoudH.-M. (2017). Real T1 relaxation time measurement and diurnal variation analysis of intervertebral discs in a healthy population of 50 volunteers.Eur. J. Radiol.8713–19. 10.1016/j.ejrad.2016.12.001
24
García-CosamalónJ.Del ValleM. E.CalaviaM. G.García-SuárezO.López-MuñizA.OteroJ.et al (2010). Intervertebral disc, sensory nerves and neurotrophins: who is who in discogenic pain: intervertebral disc, sensory nerves and neurotrophins.J. Anat.2171–15. 10.1111/j.1469-7580.2010.01227.x
25
GermaneauA.VendeuvreT.SagetM.DoumalinP.DupréJ. C.BrémandF.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.59291–303. 10.1016/j.jmbbm.2016.02.002
26
HashemiM.DadkhahP.TaheriM.KatibehP.AsadiS. (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. Pain3366–72. 10.3344/kjp.2020.33.1.66
27
HwangD.KimS.AbeydeeraN. A.StatumS.MasudaK.ChungC. B.et al (2016). Quantitative magnetic resonance imaging of the lumbar intervertebral discs.Quant. Imaging Med. Surg.6744–755. 10.21037/qims.2016.12.09
28
IshibashiK.IwaiH.KogaH. (2019). Chemonucleolysis with chondroitin sulfate ABC endolyase as a novel minimally invasive treatment for patients with lumbar intervertebral disc herniation.J. Spine Surg.5S115–S121. 10.21037/jss.2019.04.24
29
JavidM. J.NordbyE. J.FordL. T.HejnaW. J.WhislerW. W.BurtonC.et al (1983). Safety and efficacy of chymopapain (Chymodiactin) in herniated nucleus pulposus with sciatica. results of a randomized, double-blind study.JAMA2492489–2494.
30
KhanI.HargunaniR.SaifuddinA. (2014). The lumbar high-intensity zone: 20 years on.Clin. Radiol.69551–558. 10.1016/j.crad.2013.12.012
31
KnezevicN. N.MandaliaS.RaaschJ.KnezevicI.CandidoK. D. (2017). Treatment of chronic low back pain – new approaches on the horizon.J. Pain Res.101111–1123. 10.2147/JPR.S132769
32
KolberM. J.HanneyW. J. (2009). The dynamic disc model: a systematic review of the literature.Phys. Ther. Rev.14181–189. 10.1179/174328809X452827
33
LamaP.Le MaitreC. L.HardingI. J.DolanP.AdamsM. A. (2018). Nerves and blood vessels in degenerated intervertebral discs are confined to physically disrupted tissue.J. Anat.23386–97. 10.1111/joa.12817
34
LaslettM.ÖbergB.AprillC. N.McDonaldB. (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.5370–380. 10.1016/j.spinee.2004.11.007
35
Le MaitreC. L.FreemontA. J.HoylandJ. A. (2004). Localization of degradative enzymes and their inhibitors in the degenerate human intervertebral disc.J. Pathol.20447–54. 10.1002/path.1608
36
Le MaitreC. L.PockertA.ButtleD. J.FreemontA. J.HoylandJ. A. (2007). Matrix synthesis and degradation in human intervertebral disc degeneration.Biochem. Soc. Trans.35652–655. 10.1042/BST0350652
37
LiP.ZhangR.ZhouQ. (2017). Efficacy of platelet-rich plasma in retarding intervertebral disc degeneration: a meta-analysis of animal studies.BioMed. Res. Int.2017:7919201. 10.1155/2017/7919201
38
MaD.GulaniV.SeiberlichN.LiuK.SunshineJ. L.DuerkJ. L.et al (2013). Magnetic resonance fingerprinting.Nature495187–192. 10.1038/nature11971
39
MaetzelA.LiL. (2002). The economic burden of low back pain: a review of studies published between 1996 and 2001.Best Pract. Res. Clin. Rheumatol.1623–30. 10.1053/berh.2001.0204
40
ManchikantiL.SoinA.BeyaminR. (2018). An update of the systematic appraisal of theaccuracy and utility of discography in chronicspinal pain.Pain Physician191–110. 10.36076/ppj.2018.2.91
41
MayS.AinaA. (2012). Centralization and directional preference: a systematic review.Man. Ther.17497–506. 10.1016/j.math.2012.05.003
42
MayS.RungeN.AinaA. (2018). Centralization and directional preference: an updated systematic review with synthesis of previous evidence.Musculoskelet. Sci. Pract.3853–62. 10.1016/j.msksp.2018.09.006
43
McCormickZ. L.LehmanV. T.PlastarasC. T.WalegaD. R.HuddlestonP.MoussallemC.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.Spine44E1161–E1168. 10.1097/BRS.0000000000003085
44
McKenzieR. (1981). The Lumbar Spine. First.Wellington: Spinal Publication.
45
McKenzieR.MayS. (2003). The Lumbar Spine: Mechanical Diagnosis and Therapy–Volume 1.New Zealand, OC: Spinal Publications.
46
McMillanD. W.GarbuttG.AdamsM. A. (1996). Effect of sustained loading on the water content of intervertebral discs: implications for disc metabolism.Ann. Rheum. Dis.55880–887. 10.1136/ard.55.12.880
47
McNallyD. S.ShacklefordI. M.GoodshipA. E.MulhollandR. C. (1996). In vivo stress measurement can predict pain on discography.Spine212580–2587.
48
MeiselH.-J.AgarwalN.HsiehP. C.SkellyA.ParkJ.-B.BrodkeD.et al (2019). Cell therapy for treatment of intervertebral disc degeneration: a systematic review.Global Spine J.939S–52S. 10.1177/2192568219829024
49
Menezes-ReisR.SalmonC. E. G.CarvalhoC. S.BonugliG. P.ChungC. B.Nogueira-BarbosaM. H. (2015). T1rho and T2 mapping of the intervertebral disk: comparison of different methods of segmentation.AJNR Am. J. Neuroradiol.36606–611. 10.3174/ajnr.A4125
50
ModicM. T.PavlicekW.WeinsteinM. A.BoumphreyF.NgoF.HardyR.et al (1984). Magnetic resonance imaging of intervertebral disk disease. clinical and pulse sequence considerations.Radiology152103–111. 10.1148/radiology.152.1.6729099
51
MonetaG. B.VidemanT.KaivantoK.AprillC. N.SpiveyM.VanharantaH.et al (1994). Reported pain during lumbar discography as a function of anular ruptures and disc degeneration: a re-analysis of 833 discograms.Spine191968–1974.
52
MwaleF.IatridisJ. C.AntoniouJ. (2008). Quantitative MRI as a diagnostic tool of intervertebral disc matrix composition and integrity.Eur. Spine J.17432–440. 10.1007/s00586-008-0744-4
53
O’ConnellG. D.VresilovicE. J.ElliottD. M. (2007). Comparison of animals used in disc research to human lumbar disc geometry.Spine32328–333. 10.1097/01.brs.0000253961.40910.c1
54
OstiO.FraserR. (1992). MRI and discography of annular tears and intervertebral disc degeneration. a prospective clinical comparison.J. Bone Joint Surg. Br.74431–435. 10.1302/0301-620X.74B3.1587896
55
ParadisoR.AlexandreA. (2005). The different outcomes of patients with disc herniation treated either by microdiscectomy, or by intradiscal ozone injection.Acta Neurochir. Suppl.92139–142. 10.1007/3-211-27458-8_30
56
PaulC. P. L.SmitT. H.de GraafM.HolewijnR. M.BisschopA.van de VenP. 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 One13:e0191442. 10.1371/journal.pone.0191442
57
PfirrmannC. W. A.MetzdorfA.ZanettiM.HodlerJ.BoosN. (2001). Magnetic resonance classification of lumbar intervertebral disc degeneration.Spine261873–1878. 10.1097/00007632-200109010-00011
58
SachsB. L.VanharantaH.SpiveyM. A.GuyerR. D.VidemanT.RashbaumR. F.et al (1987). Dallas discogram description - a new classification of CT/discography in low-back disorders.Spine12287–294.
59
SaifuddinA.EmanuelR.WhiteJ.RentonP.BraithwaiteI.TaylorB. A. (1998). An analysis of radiating pain at lumbar discography.Eur. Spine J.7358–362. 10.1007/s005860050090
60
SchollumM. L.RobertsonP. A.BroomN. D. (2008). ISSLS prize winner: microstructure and mechanical disruption of the lumbar disc annulus: part I: a microscopic investigation of the translamellar bridging network.Spine332702–2710. 10.1097/BRS.0b013e31817bb92c
61
ShinD. S.LeeK.KimD. (2007). Biomechanical study of lumbar spine with dynamic stabilization device using finite element method.Comput.-Aided Des.39559–567. 10.1016/j.cad.2007.03.005
62
StefanakisM.LuoJ.PollintineP.DolanP.AdamsM. A. (2014). ISSLS Prize winner: mechanical influences in progressive intervertebral disc degeneration.Spine391365–1372. 10.1097/BRS.0000000000000389
63
TerttiM.PaajanenH.LaatoM.AhoH.KomuM.KormanoM. (1991). Disc degeneration in magnetic resonance imaging. a comparative biochemical, histologic, and radiologic study in cadaver spines.Spine16629–634. 10.1097/00007632-199106000-00006
64
ThompsonK. J.DagherA. P.EckelT. S.ClarkM.ReinigJ. W. (2009). Modic changes on MR images as studied with provocative diskography: clinical relevance—a retrospective study of 2457 disks.Radiology250849–855. 10.1148/radiol.2503080474
65
VanharantaH.SachsB. L.SpiveyM.GuyerR. D.HochschulerH.RashbaumR. F.et al (1987). The Relashionship of pain provocation to lumbar disc deterioration as seen by CT/Discography.Spine12295–298.
66
VosT.AllenC.AroraM.BarberR. M.BhuttaZ. A.BrownA.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.Lancet3881545–1602. 10.1016/S0140-6736(16)31678-6
67
WalkerB. F.MullerR.GrantW. D. (2003). Low back pain in australian adults: the economic burden.Asia Pac. J. Public Health1579–87. 10.1177/101053950301500202
68
WalkerJ.El AbdO.IsaacZ.MuzinS. (2008). Discography in practice: a clinical and historical review.Curr. Rev. Musculoskelet. Med.169–83. 10.1007/s12178-007-9009-9
69
WatanabeA.BennekerL. M.BoeschC.WatanabeT.ObataT.AndersonS. E. (2007). Classification of intervertebral disk degeneration with axial T2 mapping.Am. J. Roentgenol.189936–942. 10.2214/AJR.07.2142
70
WilkeH.-J.JungkunzB.WengerK.ClaesL. 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.25115–19.
71
WilkeH.-J.WengerK.ClaesL. (1998b). Testing criteria for spinal implants: recommendations for the standardization of in vitro stability testing of spinal implants.Eur. Spine J.7148–154. 10.1007/s005860050045
72
YuS. W.SetherL. A.HoP. S.WagnerM.HaughtonV. M. (1988). Tears of the anulus fibrosus: correlation between MR and pathologic findings in cadavers.AJNR Am. J. Neuroradiol.9367–370.
73
Zanjani-PourS.WinloveC. P.SmithC. W.MeakinJ. R. (2016). Image driven subject-specific finite element models of spinal biomechanics.J. Biomech.49919–925. 10.1016/j.jbiomech.2016.02.025
74
ZhangY.LenartB. A.LeeJ. K.ChenD.ShiP.RenJ.et al (2014). Histological features of endplates of the mammalian spine: from mice to men.Spine39E312–E317. 10.1097/BRS.0000000000000174
Summary
Keywords
intervertebral disc, nucleus biomechanics, quantitative MRI, low back pain, disc degeneration, radial fissure, mechanical diagnosis and therapy
Citation
Deneuville J-P, Yushchenko M, Vendeuvre T, Germaneau A, Billot M, Roulaud M, Sarracanie M, Salameh N and Rigoard P (2021) 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. Front. Bioeng. Biotechnol. 9:676003. doi: 10.3389/fbioe.2021.676003
Received
04 March 2021
Accepted
18 May 2021
Published
09 June 2021
Volume
9 - 2021
Edited by
Marwan El-Rich, Khalifa University, United Arab Emirates
Reviewed by
Maxim Bashkuev, Charité Medical University of Berlin, Germany; Mohammad Nikkhoo, Islamic Azad University, Iran
Updates
Copyright
© 2021 Deneuville, Yushchenko, Vendeuvre, Germaneau, Billot, Roulaud, Sarracanie, Salameh and Rigoard.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Arnaud Germaneau, arnaud.germaneau@univ-poitiers.fr
This article was submitted to Biomechanics, a section of the journal Frontiers in Bioengineering and Biotechnology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.