Real T1 relaxation time measurement and diurnal variation analysis of intervertebral discs in a healthy population of 50 volunteers

Variable flip angle T1 mapping and multi-echo T2 and T2* mapping magnetic resonance imaging sequences allow quantitative assessment of canine lumbar disc degeneration

Magnetic resonance imaging is the gold standard for diagnosing intervertebral disc (IVD) degeneration in dogs. However, published methods for quantifying severity or progression of IVD degeneration are currently limited. Mapping MRI sequences are used in humans for quantifying IVD degeneration but have rarely been applied in dogs. The objective of this prospective, method comparison study was to evaluate variable flip angle T1 mapping and multiecho T2 and T2* mapping as methods for quantifying canine lumbar IVD degeneration in twenty canine patients without clinical signs of spinal disease. Ventral and dorsal lumbar IVD widths were measured on radiographs, and lumbar IVDs were assigned a qualitative Pfirrmann grade based on standard T2-weighted sequences. T1, T2, and T2* relaxation times of the nucleus pulposus (NP) were measured on corresponding maps using manual-drawn ROIs. Strong intra- and interrater agreements were found (P < 0.01) for NP relaxation times. Radiographic IVD widths and T1, T2, and T2* mapping NP relaxation times were negatively correlated with Pfirrmann grading (P < 0.01). Significant differences in T1 NP relaxation times were found between Pfirrmann grade I and the other grades (P < 0.01). Significant differences in T2 and T2* NP relaxation times were found between grade I and the other grades and between grades II and III (P < 0.01). Findings indicated that T1, T2, and T2* MRI mapping sequences are feasible in dogs. Measured NP relaxation times were repeatable and decreased when Pfirrmann grades increased. These methods may be useful for quantifying the effects of regenerative treatment interventions in future longitudinal studies.

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.

Synthetic T2-weighted images of the lumbar spine derived from an accelerated T2 mapping sequence: Comparison to conventional T2w turbo spin echo

OBJECTIVES

To evaluate the diagnostic usefulness of synthetic T₂-weighted images of the lumbar spine derived from ten-fold undersampled k-space data using GRAPPATINI, a combination of a model-based approach for rapid T₂ and M₀ quantification (MARTINI) extended by generalized autocalibrating partial parallel acquistion (GRAPPA).

MATERIALS AND METHODS

Overall, 58 individuals (26 female, mean age 23.3 ± 8.1 years) were examined at 3 Tesla with sagittal and axial T₂w turbo spin echo (TSE) sequences compared to synthetic T_2-weighted contrasts derived at identical effective echo times and spatial resolutions. Two blinded readers graded disk degeneration and evaluated the lumbar intervertebral disks for present herniation or annular tear. One reader reassessed all studies after four weeks. Weighted kappa statistics were calculated to assess inter-rater and intra-rater agreement. Also, all studies were segmented manually by one reader to compute contrast ratios (CR) and contrast-to-noise ratios (CNR) of the nucleus pulposus and the annulus fibrosus.

RESULTS

Overall, the CR_T2w was 4.45 ± 1.80 and CR_T2synth was 4.71 ± 2.14. Both correlated (rsp = 0.768;p < 0.001) and differed (0.26 ± 1.38;p = 0.002) significantly. The CNR_T2w was 1.73 ± 0.52 and CNR_T2synth was 1.63 ± 0.50. Both correlated (rsp = 0.875;p < 0.001) and differed (-0.10 ± 0.25;p < 0.001) significantly. The inter-rater agreement was substantial to almost perfect (κ = 0.808-0.925) with the intra-rater agreement also substantial to almost perfect (κ = 0.862-0.963). The area under the curve of the receiver operating characteristics assessing disk herniation or annular tear ranged from 0.787 to 0.892.

CONCLUSIONS

This study concludes that synthetic images derived by GRAPPATINI can be used for clinical routine assessment with inter-rater and intra-rater agreements comparable to conventional T₂w TSE.

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.

Clinical implementation of accelerated T2 mapping: Quantitative magnetic resonance imaging as a biomarker for annular tear and lumbar disc herniation

OBJECTIVES

This study evaluates GRAPPATINI, an accelerated T2 mapping sequence combining undersampling and model-based reconstruction to facilitate the clinical implementation of T2 mapping of the lumbar intervertebral disc.

METHODS

Fifty-eight individuals (26 females, 32 males, age 23.3 ± 8.0 years) were prospectively examined at 3 T. This cohort study consisted of 19 patients, 20 rowers, and 19 volunteers. GRAPPATINI was conducted with the same parameters as a conventional 2D multi-echo spin-echo (MESE) sequence in 02:27 min instead of 13:18 min. Additional T2 maps were calculated after discarding the first echo (T2-WO1ST) and only using even echoes (T2-EVEN). Segmentation was done on the four most central slices. The resulting T2 values were compared for all four measurements.

RESULTS

T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of the nucleus pulposus of normal discs differed significantly from those of bulging discs or herniated discs (all p < 0.001). For the posterior annular region, only T2-GRAPPATINI showed a significant difference (p = 0.011) between normal and herniated discs. There was a significant difference between T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of discs with and without an annular tear for the nucleus pulposus (all p < 0.001). The nucleus pulposus' T2 at different degeneration states showed significant differences between all group comparisons of Pfirrmann grades for T2-GRAPPATINI (p = 0.000-0.018), T2-MESE (p = 0.000-0.015), T2-EVEN (p = 0.000-0.019), and T2-WO1ST (p = 0.000-0.015).

CONCLUSIONS

GRAPPATINI facilitates the use of T2 values as quantitative imaging biomarkers to detect disc pathologies such as degeneration, lumbar disc herniation, and annular tears while simultaneously shortening the acquisition time from 13:18 to 2:27 min.

KEY POINTS

• T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of the nucleus pulposus of normal discs differed significantly from those of discs with bulging or herniation (all p < 0.001). • The investigated T2 mapping techniques differed significantly in discs with and without annular tearing (all p < 0.001). • The nucleus pulposus' T2 showed significant differences between different stages of degeneration in all group comparisons for T2-GRAPPATINI (p = 0.000-0.018), T2-MESE (p = 0.000-0.015), T2-EVEN (p = 0.000-0.019), and T2-WO1ST (p = 0.000-0.015).

Impact of pulse sequence, analysis method, and signal to noise ratio on the accuracy of intervertebral disc T 2 measurement

Noninvasive assessments of intervertebral disc health and degeneration are critical for addressing disc degeneration and low back pain. Magnetic resonance imaging (MRI) is exceptionally sensitive to tissue with high water content, and measurement of the MR transverse relaxation time, T 2, has been applied as a quantitative, continuous, and objective measure of disc degeneration that is linked to the water and matrix composition of the disc. However, T 2 measurement is susceptible to inaccuracies due to Rician noise, T 1 contamination, and stimulated echo effects. These error generators can all be controlled for with proper data collection and fitting methods. The objective of this study was to identify sequence parameters to appropriately acquire MR data and to establish curve fitting methods to accurately calculate disc T 2 in the presence of noise by correcting for Rician noise. To do so, we compared T 2 calculated from the typical monoexponential (MONO) fits and noise corrected exponential (NCEXP) fits. We examined how the selected sequence parameters altered the calculated T 2 in silico and in vivo. Typical MONO fits were frequently poor due to Rician noise, and NCEXP fits were more likely to provide accurate T 2 calculations. NCEXP is particularly less biased and less uncertain at low SNR. This study showed that the NCEXP using sequences with data from 20 echoes out to echo times of ~300 ms is the best method for calculating T 2 of discs. By acquiring signal data out to longer echo times and accounting for Rician noise, the curve fitting is more robust in calculating T 2 despite the noise in the data. This is particularly important when considering degenerate discs or AF tissue because the SNR of these regions is lower.

Diurnal T2-changes of the intervertebral discs of the entire spine and the influence of weightlifting

The purpose was to study if (1) diurnal changes occur in the entire spine and if (2) intervertebral discs (IVDs) of weightlifters (WL) have decreased baseline T2-values in the morning as well as (3) increased diurnal changes throughout the day. This prospective cohort study investigated healthy volunteers between 2015 and 2017. WL were required to have participated in weightlifting ≥ 4×/week for ≥ 5 years, while non-weightlifters (NWL) were limited to < 2×/week for ≥ 5 years. Both groups underwent magnetic resonance imaging (MRI) of the entire spine in the morning and evening. WL were requested to perform weightlifting in-between imaging. IVD regions of interest (nucleus pulposus) were defined and T2-maps were measured. Analysis consisted of unpaired t-test, paired t-test, propensity-score matching (adjusting for age and sex), and Pearson correlation. Twenty-five individuals (15 [60.0%] males) with a mean age of 29.6 (standard deviation [SD 6.9]) years were analyzed. Both groups (WL: n = 12 versus [vs.] NWL: n = 13) did not differ demographic characteristics. Mean IVD T2-values of all participants significantly decreased throughout the day (95.7 [SD 15.7] vs. 86.4 [SD 13.9] milliseconds [ms]) in IVDs of the cervical (71.8 [SD 13.4] vs. 64.4 [SD 14.1] ms), thoracic (98.8 [SD 19.9] vs. 88.6 [SD 16.3] ms), and lumbar (117.0 [SD 23.7] vs. 107.5 [SD 21.6] ms) spine (P < 0.001 each). There were no differences between both groups in the morning (P = 0.635) and throughout the day (P = 0.681), even after adjusting for confounders. It can be concluded that diurnal changes of the IVDs occurred in the entire (including cervical and thoracic) spine. WL and NWL showed similar morning baseline T2-values and diurnal changes. Weightlifting may not negatively affect IVDs chronically or acutely.

Revisiting Radiographic L5-S1 Parallelism Using MRI T1 Mapping

Purpose:

Thirty years ago, we reported that parallel aspect of the L5-S1 disc on a lateral view of the spine might be considered to be an initial stage of disk degeneration. The current study represents an attempt to increase the validity of parallel sign on conventional radiograph using MR real T1 mapping.

Methods:

Forty-four young asymptomatic volunteers (mean age 21.6 ± 2.3) underwent lumbar spine MRI, twice the same day, morning and afternoon. Dedicated sequences using the inversion-recovery technique were used to calculate the T1 relaxation time. A region of interest (ROI) representing the nucleus pulposus was defined in each disk. The volunteers were stratified according to the presence or absence of a parallel morphology of L5-S1. Correlation between endplates angles, sacral slopes and T1 values were then evaluated.

Results:

L5-S1 space looks parallel for angles <10° (mean value 6.9° ± 1.4°). Sacral slope was lower in parallel disks (31.7 ± 4.9° vs. 40.1 ± 5.6°), showing a significant difference of 8.4° (p < 0.05). The T1 relaxation values show a significant difference between the two groups (p < 0.05) with a difference of 96 ms for the morning (1090.9 ± 33.3 ms for the parallel group and 1186.9 ± 41.2 ms for the non-parallel) and 121.9 ms for the afternoon (respectively 1004.7 ± 22.2 ms and 1126.6 ± 12.9 ms).

Conclusion:

The difference between the two groups suggests that parallel morphology of the L5-S1 disk is associated with lower water content.