MRI quantification of human spine cartilage endplate geometry: Comparison with age, degeneration, level, and disc geometry

Qualitative and Quantitative MR Imaging of the Cartilaginous Endplate: A Review

The cartilaginous endplate (CEP) plays a pivotal role in facilitating the supply of nutrients and, transport of metabolic waste, as well as providing mechanical support for the intervertebral disc (IVD). Recent technological advances have led to a surge in MR imaging studies focused on the CEP. This article describes the anatomy and functions of the CEP as well as MRI techniques for both qualitative and quantitative assessment of the CEP. Effective CEP MR imaging sequences require two key features: high spatial resolution and relatively short echo time. High spatial resolution spoiled gradient echo (SPGR) and ultrashort echo time (UTE) sequences, fulfilling these requirements, are the basis for most of the sequences employed in CEP imaging. This article reviews existing sequences for qualitative CEP imaging, such as the fat-suppressed SPGR and UTE, dual-echo subtraction UTE, inversion recovery prepared and fat-suppressed UTE, and dual inversion recovery prepared UTE sequences. These sequences are employed together with other techniques for quantitative CEP imaging, including measurements of T2*, T2, T1, T1ρ, magnetization transfer, perfusion, and diffusion tensor parameters. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 2.

Combining the probabilistic finite element model and artificial neural network to study nutrient levels in the human intervertebral discs

BACKGROUND

Diffusion distance and diffusivity are known to affect nutrient transport rates, but the probabilistic analysis of these two factors remains vacant. There is a lack of effective tools to evaluate disc nutrient levels.

METHODS

Five-hundred-disc samples with different combinations of morphological and water content parameters were generated, which were used to evaluate nutrient levels in unloaded and loaded states. Spearman correlation coefficients between inputs and responses were calculated. Artificial neural networks were trained to predict nutrient concentrations based on the dataset generated by the probabilistic finite element model.

FINDINGS

In unloaded and loaded states, the minimum oxygen concentration of nucleus pulposus was negatively correlated with disc height (r = -0.83, p < 0.01 and r = -0.76, p < 0.01, respectively), and the minimum glucose concentration of annulus fibrosus was positively correlated with its water content (r = 0.68, p < 0.01 and r = 0.73, p < 0.01, respectively). The maximum lactate concentration of cartilage endplate was affected by endplate thickness (r = 0.94, p < 0.01 and r = 0.95, p < 0.01, respectively). For trained neural networks, nutrient concentrations could be well predicted, with coefficients of determination greater than 0.95 and mean absolute percentage errors less than 5 %.

INTERPRETATION

This study underscores the importance of disc height, annulus fibrosus water content, and endplate thickness in regulating nutrient levels, and precise control of these parameters should be prioritized in the design of tissue-engineered discs. Moreover, artificial neural networks might be a promising tool for evaluating nutrient levels.

The influence of geometry on intervertebral disc stiffness

Geometry plays an important role in intervertebral disc (IVD) mechanics. Previous computational studies have found a link between IVD geometry and stiffness. However, few experimental studies have investigated this link, possibly due to difficulties in non-destructively quantifying internal geometric features. Recent advances in ultra-high resolution MRI provides the opportunity to visualise IVD features in unprecedented detail. This study aimed to quantify 3D human IVD geometries using 9.4 T MRIs and to investigate correlations between geometric variations and IVD stiffness. Thirty human lumbar motion segments (fourteen non-degenerate and sixteen degenerate) were scanned using a 9.4 T MRI and geometric parameters were measured. A 1kN compressive load was applied to each motion segment and stiffness was calculated. Degeneration caused a reduction (p < 0.05) in IVD height, a decreased nucleus-annulus area ratio, and a 1.6 ± 3.0 mm inward collapse of the inner annulus. The IVD height, anteroposterior (AP) width, lateral width, cross-sectional area, nucleus-annulus boundary curvature, and nucleus-annulus area ratio had a significant (p < 0.05) influence on IVD stiffness. Linear relationships (p < 0.05, r > 0.47) were observed between these geometric features and IVD compressive stiffness and a multivariate regression model was generated to enable stiffness to be predicted from features observable on clinical imaging (stiffness, N/mm = 6062 - (61.2 × AP width, mm) - (169.2 × IVD height, mm)). This study advances our understanding of disc structure-function relationships and how these change with degeneration, which can be used to both generate and validate more realistic computational models.

Consolidating and re-evaluating the human disc nutrient microenvironment

Background

Despite exciting advances in regenerative medicine, cell-based strategies for treating degenerative disc disease remain in their infancy. To maximize the potential for successful clinical translation, a more thorough understanding of the in vivo microenvironment is needed to better determine and predict how cell therapies will respond when administered in vivo.

Aims

This work aims to reflect on the in vivo nutrient microenvironment of the degenerating IVD through consolidating what has already been measured together with investigative in silico models.

Materials and Methods

This work uses in silico modeling, underpinned by more recent experimentally determined parameters of degeneration and nutrient transport from the literature, to re-evaluate the current knowledge in terms of grade-specific stages of degeneration.

Results

Through modeling only the metabolically active cell population, this work predicts slightly higher glucose concentrations compared to previous in silico models, while the predicted results show good agreement with previous intradiscal pH and oxygen measurements. Increasing calcification with degeneration limits nutrient transport into the IVD and initiates a build-up of acidity; however, its effect is compensated somewhat by a reduction in diffusional distance due to decreasing disc height.

Discussion

This work advances in silico modeling through a strong foundation of experimentally determined grade-specific input parameters. Taken together, pre-existing measurements and predicted results suggest that metabolite concentrations may not be as critically low as commonly believed, with calcification not appearing to have a detrimental effect at stages of degeneration when cell therapies are an appropriate intervention.

Conclusion

Overall, our initiative is to provoke greater deliberation and consideration of the nutrient microenvironment when performing in vitro cell culture and cell therapy development. This work highlights urgency for robust experimental glucose measurements in healthy and degenerating IVDs, not only to validate in silico models but to significantly advance the field in fully elucidating the nutrient microenvironment and refining in vitro techniques to accelerate clinical translation.

Applications of Functionalized Hydrogels in the Regeneration of the Intervertebral Disc

Intervertebral disc degeneration (IDD) is caused by genetics, aging, and environmental factors and is one of the leading causes of low back pain. The treatment of IDD presents many challenges. Hydrogels are biomaterials that possess properties similar to those of the natural extracellular matrix and have significant potential in the field of regenerative medicine. Hydrogels with various functional qualities have recently been used to repair and regenerate diseased intervertebral discs. Here, we review the mechanisms of intervertebral disc homeostasis and degeneration and then discuss the applications of hydrogel-mediated repair and intervertebral disc regeneration. The classification of artificial hydrogels and natural hydrogels is then briefly introduced, followed by an update on the development of functional hydrogels, which include noncellular therapeutic hydrogels, cellular therapeutic hydrogel scaffolds, responsive hydrogels, and multifunctional hydrogels. The challenges faced and future developments of the hydrogels used in IDD are discussed as they further promote their clinical translation.

Diagnostic Role of Magnetic Resonance Imaging in Low Back Pain Caused by Vertebral Endplate Degeneration

Low back pain (LBP) is a common health issue worldwide with a huge economic burden on healthcare systems. In the United States alone, the cost is estimated to be $100 billion each year. Intervertebral disc degeneration is considered one of the primary causes of LBP. Moreover, the critical role of the vertebral endplates in disc degeneration and LBP is becoming apparent. Endplate abnormalities are closely correlated with disc degeneration and pain in the lumbar spine. Imaging modalities such as plain film radiography, computed tomography, and fluoroscopy are helpful but not very effective in detecting the causes behind LBP. Magnetic resonance imaging (MRI) can be used to acquire high-quality three-dimensional images of the lumbar spine without using ionizing radiation. Therefore, it is increasingly being used to diagnose spinal disorders. However, according to the American College of Radiology, current referral and justification guidelines for MRI are not sufficiently clear to guide clinical practice. This review aimed to evaluate the role of MRI in diagnosing LBP by considering the correlative contributions of vertebral endplates. The findings of the review indicate that MRI allows for fine evaluations of endplate morphology, endplate defects, diffusion and perfusion properties of the endplate, and Modic changes. Changes in these characteristics of the endplate were found to be closely correlated with disc degeneration and LBP. The collective evidence from the literature suggests that MRI may be the imaging modality of choice for patients suffering from LBP. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY STAGE: 3.

Changes in Lumbar Endplate Area and Concavity Associated With Disc Degeneration

STUDY DESIGN

Retrospective image-based analysis.

OBJECTIVE

To measure endplate three-dimensional (3D) geometry, endplate changes in vivo and to investigate correlations between disc degeneration and endplate 3D geometry dependent on symptoms of low back pain (LBP).

SUMMARY OF BACKGROUND DATA

It has been hypothesized that alteration of load transmission from the nucleus pulposus to the annulus fibrosus affects vertebral endplate geometry.

METHODS

3D surface models of inferior/superior lumbar endplates were created from computed tomography scans of n = 92 volunteers with and without LBP. Disc degeneration was evaluated using Pfirrmann scale. Concavity in both coronal and sagittal planes was assessed with the Concavity Index (unitless; larger than 1: concave; flat: 1; and less than 1: convex, respectively). Endplate area and disc height distribution were computed and the effects from demographics and spinal degeneration were sought with an analysis of variance model.

RESULTS

Both sagittal and coronal planes revealed significantly decreased concavity in those with terminal grade 5 disc degeneration (mean 0.833 ± 0.235) compared to the other grades in the cohort. Older subjects presented with larger endplate areas than the younger subjects (P = 0.0148) at L4-S1. Overall, symptomatic subjects had significantly larger endplate areas (P = 0.022), especially at the lower lumbar levels (P < 0.001). Analysis of variance showed that sex, disc level, disc degeneration grade, and disc height reached significance (P < 0.0001) as influential parameters in both Concavity Index cases.

CONCLUSION

With advancing intervertebral disc degeneration, endplates become more convex over time in both sagittal and coronal planes. Our findings implicate the endplate changes with advancing disc degeneration in the shift in load transmission from the nucleus pulposus to the annulus fibrosus, resulting in changes within the curvature of the endplates. This is also the first study to describe the direct impact of age, sex, and LBP on vertebral endplate anatomy.

LEVEL OF EVIDENCE

5.

Cartilage Endplate Thickness Variation Measured by Ultrashort Echo-Time MRI Is Associated With Adjacent Disc Degeneration

STUDY DESIGN

A magnetic resonance imaging study of human cadaver spines.

OBJECTIVE

To investigate associations between cartilage endplate (CEP) thickness and disc degeneration.

SUMMARY OF BACKGROUND DATA

Damage to the CEP is associated with spinal injury and back pain. However, CEP morphology and its association with disc degeneration have not been well characterized.

METHODS

Ten lumbar motion segments with varying degrees of disc degeneration were harvested from six cadaveric spines and scanned with magnetic resonance imaging in the sagittal plane using a T2-weighted two-dimensional (2D) sequence, a three-dimensional (3D) ultrashort echo-time (UTE) imaging sequence, and a 3D T1ρ mapping sequence. CEP thicknesses were calculated from 3D UTE image data using a custom, automated algorithm, and these values were validated against histology measurements. Pfirrmann grades and T1ρ values in the disc were assessed and correlated with CEP thickness.

RESULTS

The mean CEP thickness calculated from UTE images was 0.74 ± 0.04 mm. Statistical comparisons between histology and UTE-derived measurements of CEP thickness showed significant agreement, with the mean difference not significantly different from zero (P = 0.32). Within-disc variation of T1ρ (standard deviation) was significantly lower for Pfirrmann grade 4 than Pfirrmann grade 3 (P < 0.05). Within-disc variation of T1ρ and adjacent CEP thickness heterogeneity (coefficient of variation) had a significant negative correlation (r = -0.65, P = 0.04). The standard deviation of T1ρand the mean CEP thickness showed a moderate positive correlation (r = 0.40, P = 0.26).

CONCLUSION

This study demonstrates that quantitative measurements of CEP thickness measured from UTE magnetic resonance imaging are associated with disc degeneration. Our results suggest that variability in CEP thickness and T1ρ, rather than their mean values, may serve as valuable diagnostic markers for disc degeneration.

LEVEL OF EVIDENCE

N/A.

Axial loading during MRI induces significant T2 value changes in vertebral endplates-a feasibility study on patients with low back pain

BACKGROUND

The function of the endplate (EP) is the most important factor influencing nutritional supply to the avascular intervertebral disc (IVD). It is desired to have a non-invasive method to assess functional EP characteristics in vivo. Assessment of functional EP characteristics is important in order to understand its relation to IVD degeneration, which in turn might deepen the understanding of the pathophysiology behind low back pain (LBP). It was hypothesized that, by comparing quantitative MRI of EPs performed with conventional supine MRI (unloaded MRI) with axial loading during MRI (alMRI), dynamical properties of the EP can be displayed. The aim was therefore to investigate the feasibility of axial loading during MRI (alMRI) to instantaneously induce quantitative EP changes.

METHODS

T2 mapping of 55 vertebral EPs (L1-S1) in five LBP patients was performed during conventional supine MRI (unloaded MRI) and subsequent alMRI. With T2 mapping, the cartilaginous EP and bony EP cannot be separated; hence, the visualized EP was termed EP zone (EPZ). Each EPZ was segmented at multiple midsagittal views, generating volumetric regions of interest. EPZs demonstrating signal inhomogeneity and/or adjacent Modic changes (MC) were termed abnormal EPZs. EPZ mean T2 values were compared between unloaded MRI and alMRI, and their relationship with abnormal EPZs was determined.

RESULTS

alMRI induced significantly higher (p = 0.01) EPZ mean T2 values compared with unloaded MRI. Significantly higher mean T2 values were seen in inferior EPZs compared with superior EPZs, both with unloaded MRI (35%, p < 0.001) and with alMRI (26%, p = 0.04). Significant difference between unloaded MRI and alMRI was seen in normal (p = 0.02), but not in abnormal EPZs (p = 0.5; n = 12).

CONCLUSIONS

alMRI induces changes in human EPZ characteristics in vivo. The T2 value significantly increased in normal EPZs, with lack of such in abnormal EPZs. Combining T2 mapping with alMRI provides a clinical feasible, non-invasive method with potential to reveal biochemical behavioral patterns, thus adding another dimension of the EPZs characteristics compared with information obtained with solely unloaded MRI.

Structure-function relationships at the human spinal disc-vertebra interface

Damage at the intervertebral disc-vertebra interface associates with back pain and disc herniation. However, the structural and biomechanical properties of the disc-vertebra interface remain underexplored. We sought to measure mechanical properties and failure mechanisms, quantify architectural features, and assess structure-function relationships at this vulnerable location. Vertebra-disc-vertebra specimens from human cadaver thoracic spines were scanned with micro-computed tomography (μCT), surface speckle-coated, and loaded to failure in uniaxial tension. Digital image correlation (DIC) was used to calculate local surface strains. Failure surfaces were scanned using scanning electron microscopy (SEM), and adjacent sagittal slices were analyzed with histology and SEM. Seventy-one percent of specimens failed initially at the cartilage endplate-bone interface of the inner annulus region. Histology and SEM both indicated a lack of structural integration between the cartilage endplate (CEP) and bone. The interface failure strength was increased in samples with higher trabecular bone volume fraction in the vertebral endplates. Furthermore, failure strength decreased with degeneration, and in discs with thicker CEPs. Our findings indicate that poor structural connectivity between the CEP and vertebra may explain the structural weakness at this region, and provide insight into structural features that may contribute to risk for disc-vertebra interface injury. The disc-vertebra interface is the site of failure in the majority of herniation injuries. Here we show new structure-function relationships at this interface that may motivate the development of diagnostics, prevention strategies, and treatments to improve the prognosis for many low back pain patients with disc-vertebra interface injuries. © 2017 The Authors. Journal of Orthopaedic Research® Published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 36:192-201, 2018.

Ultrashort time to echo magnetic resonance techniques for the musculoskeletal system

Magnetic resonance (MR) imaging has been widely implemented as a non-invasive modality to investigate musculoskeletal (MSK) tissue disease, injury, and pathology. Advancements in MR sequences provide not only enhanced morphologic contrast for soft tissues, but also quantitative biochemical evaluation. Ultrashort time to echo (UTE) sequence, in particular, enables novel morphologic and quantitative evaluation of previously unseen MSK tissues. By using short minimum echo times (TE) below 1 msec, the UTE sequence can unveil short T2 properties of tissues including the deepest layers of the articular cartilage, cartilaginous endplate at the discovertebral junction, the meniscus, and the cortical bone. This article will discuss the application of UTE to evaluate these MSK tissues, starting with tissue structure, MR imaging appearance on standard versus short and ultrashort TE sequences, and provide the range of quantitative MR values found in literature.