Comparison of three approaches for defining nucleus pulposus and annulus fibrosus on sagittal magnetic resonance images of the lumbar spine

Morphologic Changes of the Intervertebral Disk During Growth

STUDY DESIGN

Cross-sectional.

OBJECTIVE

The aim of this study was to describe morphologic changes of the annulus fibrosus (AF) and nucleus pulposus (NP) in children during growth using magnetic resonance imaging.

SUMMARY OF BACKGROUND DATA

Little is known of intervertebral disk (IVD) maturation as opposed to degeneration, such as changes in relative AF/NP proportions and orientation during growth. Studies suggest that IVD plays a role in the etiology of pediatric spinal deformities. Therefore, understanding the morphologic development of the AF and NP during growth is key.

MATERIALS AND METHODS

An existing database of children aged 0 to 18 that had magnetic resonance imaging for indications unrelated to the spine were analyzed. The AF/NP were segmented semiautomatically from T1 to L5. The parameters: mean IVD height, cross-sectional area, slenderness (height/width ratio), volume (ratio), and relative position of the centroid of the NP within the IVD in three directions ( x , y , z ) were extracted, and compared between age, sex, and spinal level.

RESULTS

IVD height increased modestly and predominantly in the low-thoracic and lumbar spine during the first 5 to 10 years of life. Cross-sectional area and thus volume increased steadily at all levels throughout growth. IVD slenderness decreased sharply in the first years of life and remains relatively stable throughout the remainder of growth. IVDs were smaller and more slender in females, especially in the mid-thoracic spine at early adolescence. In the upper-thoracic and mid-thoracic spine the NP comprises 10% to 12% of total IVD volume during growth, this percentage increases in the low-thoracic and lumbar spine towards 20% to 25%. In the anterior-posterior direction, the position of the nucleus increasingly shifts with age, possibly in line with the developing sagittal profile of the spine.

CONCLUSION

This study describes the development of thoracic and lumbar IVDs during growth and may be used as a reference for future studies on the role of IVD in the etiology of disk-related disorders.

Lumbar intervertebral disc segmentation for computer modeling and simulation

BACKGROUND AND OBJECTIVE

The present work had as its main objective the development of a method for localizing and automatically segmenting lumbar intervertebral discs (IVD) in 3D from magnetic resonance imaging (MRI), with the goal of supporting the generation of finite element (FE) models from actual lumbar spine anatomy, by providing accurate and personalized information on the shape of the patient's IVD. The extension of the method to allow performing separate segmentations of the IVD's two main structures - annulus fibrosus (AF) and nucleus pulposus (NP) - as well as automatically detecting degenerated IVD where this distinction is no longer possible was also an objective of the work.

METHODS

The method presented here evolves from 2D segmentations in the sagittal profile using Gabor filters towards 3D segmentations. It works by detecting the spine curves and intensity regions corresponding to IVD. As so, the 2D method from Zhu et al. (2013) was partially implemented, modified and adapted to 3D use, and then tested with eight spines from two separated online datasets. The 3D adaptation was achieved by using vertebral body segmentation masks to approximate the shape of the vertebrae and to adjust the spine curves accordingly.

RESULTS

The method showed average values of 85%, 83% and 96% for the Dice coefficient, sensitivity and specificity, respectively. The method correctly identified 65 of 68 (96%) IVD as either healthy or degenerated. The method's Dice coefficient is within the range of existing 3D IVD segmentation methods in the literature (81-92%). The method took on average 6-7 s to perform a full 3D segmentation, which is well within the range of the existing methods (2 s - 19 min).

CONCLUSIONS

The developed method can be used to generate accurate 3D models of the IVD based on MRI, with AF/NP distinction and detection of marked degeneration by comparing each IVD with the remaining spine levels. Further work shall improve the method towards distinguishing between specific levels of degeneration for clinically oriented FE modeling.

Delineation of Inner Annulus Fibrosus and Nucleus Pulposus on Routine T2-weighted MR Images

STUDY DESIGN

Retrospective study of 150 IVDs.

OBJECTIVE

Assessment of costume algorithm ability to delineate the IAF and NP on routine T2 images.

SUMMARY OF BACKGROUND DATA

Central hyperintense region on T2-weighted MR images of normal lumbar IVDs represents a combination of IAF and NP. Ability to identify NP as distinct from IAF can help improve our understanding of IVD morphology in-vivo.

METHODS

Sagittal T2-weighted TSE MR images of 150 lumbar IVDs from 25 patients were analyzed. MR images were processed using a custom algorithm that markedly increased the signal intensity of structures with inherent signal intensity within 2 defined intensity thresholds. Signal intensity and contrast-to- noise ratio between outer annulus fibrosus, IAF, and NP were assessed at baseline and after processing. To assess consistency of underlying T2 differences, similar analysis was done on 108 discs from 18 patients in whom additional sagittal T2-weighted STIR images were available.

RESULTS

Following image processing, apparent IAF and NP were rendered visible in 86% and 84.3% IVDs on T2-weighted TSE and STIR images respectively. While signal intensity of these 2 regions was inherently different (P< 0.001) before processing on TSE and STIR images, their visualization was facilitated by a significant increase (P<0.001) in contrast-to-noise ratio after processing. Nonvisualization of NP was associated with disc degeneration (P<0.001).

CONCLUSION

Inherent differences exist in signal intensities of normal NP and IAF on T2-weighted MR images. Accentuating these differences using image postprocessing techniques can render these 2 structures visible.

Ultrashort time-to-echo T2* and T2* relaxometry for evaluation of lumbar disc degeneration: a comparative study

Background

To compare potential of ultrashort time-to-echo (UTE) T2* mapping and T2* values from T2*-weighted imaging for assessing lumbar intervertebral disc degeneration (IVDD),with Pfirrmann grading as a reference standard.

Methods

UTE-T2* and T2* values of 366 lumbar discs (L1/2-L5/S1) in 76 subjects were measured in 3 segmented regions: anterior annulus fibrosus, nucleus pulposus (NP), and posterior annulus fibrosus. Lumbar intervertebral discs were divided into 3 categories based on 5-level Pfirrmann grading: normal (Pfirrmann grade I),early disc degeneration (Pfirrmann grades II-III), and advanced disc degeneration (Pfirrmann grades IV-V). Regional differences between UTE-T2* and T2* relaxometry and correlation with degeneration were statistically analyzed.

Results

UTE-T2* and T2*value correlated negatively with Pfirrmann grades (P < 0.001). In NP, correlations with Pfirrmann grade were high with UTE-T2* values (r =  − 0.733; P < 0.001) and moderate with T2* values (r = -0.654; P < 0.001). Diagnostic accuracy of detecting early IVDD was better with UTE-T2* mapping than T2* mapping (P < 0.05),with receiver operating characteristic analysis area under the curve of 0.715–0.876.

Conclusions

UTE-T2* relaxometry provides another promising magnetic resonance imaging sequence for quantitatively evaluate lumbar IVDD and was more accurate than T2*mapping in the earlier stage degenerative process.

Comparison of MRI T1, T2, and T2* mapping with histology for assessment of intervertebral disc degeneration in an ovine model

An easy, reliable, and time-efficient standardized approach for assessing lumbar intervertebral disc (IVD) degeneration with relaxation times measurements in pre-clinical and clinical studies is lacking. This prospective study aims to determine the most appropriate method for lumbar IVD degeneration (IDD) assessment in sheep by comparing three quantitative MRI sequences (variable-flip-angle T1 mapping, and multi-echo T2 and T2* mapping), correlating them with Pfirrmann grading and histology. Strong intra- and interrater agreements were found for Nucleus pulposus (NP) regions-of-interest (ROI). T1, T2, and T2* mapping correlated with Pfirrmann grading and histological scoring (p < 0.05) except for the most ventral rectangular ROI on T2 maps. Correlations were excellent for all of the T1 ROIs and the T2* NP ROIs. Highly significant differences in T1 values were found between all Pfirrmann grades except between grades I/II and between grades III/IV. Significant differences were identified in the T2 and the T2* values between all grades except between grades I/III. T1, T2, and T2* relaxation times measurements of the NP are an accurate and time-efficient tool to assess lumbar IDD in sheep. Variable-flip-angle T1 mapping may be further considered as a valuable method to investigate IDD and to assess the efficacy of regenerative treatments in longitudinal studies.

Advanced Quantitative Spine Imaging

Although advanced quantitative imaging may not be currently used to any degree in the routine reporting of spinal examinations, this situation will change in the not too distant future. Advanced quantitative imaging has already allowed us to understand a great deal more regarding spinal development, marrow physiology, and disease pathogenesis. Radiologists are ideally suited to drive this research forward. To speed up this process and optimize the impact of studies reporting spine quantitative data, we should work toward universal standards on the acquisition of spine data that will allow quantitative studies to be more easily compared, contrasted, and amalgamated.

Characterization of Intervertebral Disc Changes in Asymptomatic Individuals with Distinct Physical Activity Histories Using Three Different Quantitative MRI Techniques

(1) Background: Assessments of intervertebral disc (IVD) changes, and IVD tissue adaptations due to physical activity, for example, remains challenging. Newer magnetic resonance imaging techniques can quantify detailed features of the IVD, where T2-mapping and T2-weighted (T2w) and Dixon imaging are potential candidates. Yet, their relative utility has not been examined. The performances of these techniques were investigated to characterize IVD differences in asymptomatic individuals with distinct physical activity histories. (2) Methods: In total, 101 participants (54 women) aged 25–35 years with distinct physical activity histories but without histories of spinal disease were included. T11/12 to L5/S1 IVDs were examined with sagittal T2-mapping, T2w and Dixon imaging. (3) Results: T2-mapping differentiated Pfirrmann grade-1 from all other grades (p < 0.001). Most importantly, T2-mapping was able to characterize IVD differences in individuals with different training histories (p < 0.005). Dixon displayed weak correlations with the Pfirrmann scale, but presented significantly higher water content in the IVDs of the long-distance runners (p < 0.005). (4) Conclusions: Findings suggested that T2-mapping best reflects IVD differences in asymptomatic individuals with distinct physical activity histories changes. Dixon characterized new aspects of IVD, probably associated with IVD hypertrophy. This complementary information may help us to better understand the biological function of the disc.

Differences in IVD characteristics between low back pain patients and controls associated with HIZ as revealed with quantitative MRI

Background

Magnetic resonance imaging (MRI) can provide objective continuous intervertebral disc (IVD) measures in low back pain (LBP) patients. However, there are limited studies comparing quantitative IVD measures of symptomatic and asymptomatic individuals.

Purpose

This study aimed to investigate possible differences in IVD tissue composition in patients with chronic LBP and controls using quantitative MRI and correlate IVD measures with the phenotype High-Intensity Zone (HIZ).

Methods

The lumbar spine of 25 LBP-patients (25-69y, mean 38y, 11 males) and 12 controls (25-59y, mean 38y, 7 males) was examined with T2-mapping on a 1.5T MRI scanner. The mean T2-map value and standard deviation were determined in three midsagittal IVD slices and five sub-regions dividing each IVD in the sagittal plane. The distribution of T2-map values over the IVD was also determined with histogram analysis (Δμ = distribution width).

Results

When compared to controls, patient IVDs displayed lower values for all metrics, with significant differences for the T2-map value, standard deviation (p = 0.026) and Δμ (p = 0.048). Significantly different T2-map values were found between cohorts in the region representing nucleus pulposus and the border zone between nucleus pulposus and posterior annulus fibrosus (p = 0.047–0.050). Excluding all IVDs with HIZs resulted in no significant difference between the cohorts for any of the analyzed metrics (p = 0.053–0.995). Additionally, the T2-map values were lower in patients with HIZ in comparison without HIZ.

Conclusions

Differences in IVD characteristics, measured with quantitative MRI, between LBP patients and controls were found. The T2-differences may reflect altered IVD function associated with HIZ. Future studies are recommended to explore IVD functionality in relation to HIZ and LBP.

The importance of level stratification for quantitative MR studies of lumbar intervertebral discs: a cross-sectional analysis in 101 healthy adults

PURPOSE

To investigate whether quantitative T2-times depend on lumbar intervertebral disc (IVD) level.

METHODS

The lumbar spine (Th12/L1-L5/S1) of 101 participants (53.5% female, 30.0[± 3.6]years, 173.5[± 9.6]cm and 69.9[± 13.4]kg), without history of back pain, was examined on a 3T scanner with sagittal T2-mapping. All IVDs were stratified according to Pfirrmann grade and lumbar level, with mean T2-time determined for the entire IVD volume and in five subregions of interests.

RESULTS

Significant level-dependent T2-time differences were detected, both for the entire IVD volume and its subregions. For the entire IVD volume, Pfirrmann grade 2 IVDs displayed 9-18% higher T2-times in Th12/L1 IVDs compared to L2/L3-L5/S1 IVDs (0.001 > p < 0.004) and significantly different T2-times in L1/L2-L2/L3 IVDs compared to most of the IVDs in the lower lumbar spine. In Pfirrmann grades 1, 3 and 4 IVDs, no significant level-dependent T2-time differences were observed for the entire IVD. More pronounced results were observed when comparing IVD subregions, with significant level-dependent differences also within Pfirrmann grade 1 and grade 3 IVDs. For example, in posterior IVD subregions mean T2-time was 80-82% higher in Th12/L1 compared to L3/L4-L4/L5 Pfirrmann grade 1 IVDs (p < 0.05) and 10-14% higher in L5/S1 compared to L3/L4-L4/L5 Pfirrmann grade 3 IVDs (0.02 > p < 0.001).

DISCUSSION

Significant level-dependent T2-time differences within several Pfirrmann grades, both for the entire IVD volume and for multiple IVD subregions, were shown in this large cohort study. The T2-time differences between levels existed in both non-degenerated and degenerated IVDs. These findings show the importance of stratifying for lumbar level when quantitative IVD studies are performed using T2-mapping. These slides can be retrieved under Electronic Supplementary Material.

Axial loading during MRI reveals deviant characteristics within posterior IVD regions between low back pain patients and controls

PURPOSE

To investigate differences in functional intervertebral disk (IVD) characteristics between low back pain (LBP) patients and controls using T2-mapping with axial loading during MRI (alMRI).

METHODS

In total, 120 IVDs in 24 LBP patients (mean age 39 years, range 25-69) were examined with T2-mapping without loading of the spine (uMRI) and with alMRI (DynaWell^® loading device) and compared with 60 IVDs in 12 controls (mean age 38 years, range 25-63). The IVD T2-value was acquired after 20-min loading in five regions of interests (ROI), ROI1-5 from anterior to posterior. T2-values were compared between loading states and cohorts with adjustment for Pfirrmann grade.

RESULTS

In LBP patients, mean T2-value of the entire IVD was 64 ms for uMRI and 66 ms for alMRI (p = 0.03) and, in controls, 65 ms and 65 ms (p = 0.5). Load-induced T2-differences (alMRI-uMRI) were seen in all ROIs in both patients (0.001 > p < 0.005) and controls (0.0001 > p < 0.03). In patients, alMRI induced an increase in T2-value for ROI1-3 (23%, 18% and 5%) and a decrease for ROI4 (3%) and ROI5 (24%). More pronounced load-induced decrease was detected in ROI4 in controls (9%/p = 0.03), while a higher absolute T2-value was found for ROI5 during alMRI in patients (38 ms) compared to controls (33 ms) (p = 0.04).

CONCLUSION

The alMRI-induced differences in T2-value in ROI4 and ROI5 between patients and controls most probably indicate biomechanical impairment in the posterior IVD regions. Hence, alMRI combined with T2-mapping offers an objective and clinical feasible tool for biomechanical IVD characterization that may deepen the knowledge regarding how LBP is related to altered IVD matrix composition. These slides can be retrieved under Electronic Supplementary Material.

Lumbar intervertebral discs T2 relaxometry and T1ρ relaxometry correlation with age in asymptomatic young adults

BACKGROUND

To investigate the detection of intervertebral disc (IVD) composition aging-related changes using T2 and T1ρ relaxometry in vivo in asymptomatic young adults.

METHODS

We recruited ninety asymptomatic and young adults (42 men and 48 women) between 20 and 40 years old. T2 and T1ρ lumbar spine mappings were acquired using 1.5 T magnetic resonance imaging (MRI) scanner. Two independent observers manually segmented 450 lumbar discs in all slices. They also performed sub region segmentation of annulus fibrosus (AF) and nucleus pulposus (NP) at the central MRI sagittal slices.

RESULTS

There was no difference between men and women for T2 (P=0.37) or T1ρ relaxometry (P=0.97). There was a negative correlation between age (20-40 years) and IVD T2 relaxation time of the whole disc (r=-0.30, P<0.0001), NP (r=-0.20 to -0.51, P<0.05) and posterior AF (r=-0.21 to -0.31, P<0.05) at all lumbar disc levels. There was no statistical correlation between aging and IVD T1ρ relaxation both for NP and AF.

CONCLUSIONS

T2 relaxometry detected gradual IVD dehydration in the first two decades of adulthood. We observed no significant variation of T1ρ or volumetry with aging in our study group. Our results suggest that T2 mapping may be more appropriate to detect early IVD aging changes.