The contributions of cartilage endplate composition and vertebral bone marrow fat to intervertebral disc degeneration in patients with chronic low back pain.

PURPOSE: The composition of the subchondral bone marrow and cartilage endplate (CEP) could affect intervertebral disc health by influencing vertebral perfusion and nutrient diffusion. However, the relative contributions of these factors to disc degeneration in patients with chronic low back pain (cLBP) have not been quantified. The goal of this study was to use compositional biomarkers derived from quantitative MRI to establish how CEP composition (surrogate for permeability) and vertebral bone marrow fat fraction (BMFF, surrogate for perfusion) relate to disc degeneration.
METHODS: MRI data from 60 patients with cLBP were included in this prospective observational study (28 female, 32 male; age = 40.0 ± 11.9 years, 19-65 [mean ± SD, min-max]). Ultra-short echo-time MRI was used to calculate CEP T2* relaxation times (reflecting biochemical composition), water-fat MRI was used to calculate vertebral BMFF, and T1ρ MRI was used to calculate T1ρ relaxation times in the nucleus pulposus (NP T1ρ, reflecting proteoglycan content and degenerative grade). Univariate linear regression was used to assess the independent effects of CEP T2* and vertebral BMFF on NP T1ρ. Mixed effects multivariable linear regression accounting for age, sex, and BMI was used to assess the combined relationship between variables.
RESULTS: CEP T2* and vertebral BMFF were independently associated with NP T1ρ (p = 0.003 and 0.0001, respectively). After adjusting for age, sex, and BMI, NP T1ρ remained significantly associated with CEP T2* (p = 0.0001) but not vertebral BMFF (p = 0.43).
CONCLUSION: Poor CEP composition plays a significant role in disc degeneration severity and can affect disc health both with and without deficits in vertebral perfusion.

Introduction

While intervertebral disc degeneration is prevalent in people without chronic low back pain (cLBP), there is a strong implication that disc degeneration causes cLBP in certain subgroups [1] and consequently, therapeutically slowing or reversing disc degeneration is a common but as yet unrealized clinical goal. Cell-based biologic therapies to regenerate the disc have demonstrated mixed efficacy in clinical trials, and a key obstacle to successful translation of these therapies is an incomplete understanding of the mechanisms of cellular homeostasis and the factors disrupting homeostasis in degenerating discs [2, 3]. One specific knowledge gap is the uncertain role of cartilage endplate (CEP) composition relative to vertebral vascularity. Both endogenous and implanted cells inside the nucleus pulposus (NP) are nourished by glucose and oxygen diffusing from the vertebral capillary bed, across the CEP, and into the NP [4]. In vitro experiments demonstrate that physiologic fluctuations in CEP composition (affecting its permeability), as well as endplate-adjacent vascularity (providing the reservoir of available nutrients), can impact NP cell viability and function [4, 5]. Despite these insights, the contributions of CEP composition and vertebral vascularity to disc degeneration in humans is unknown. Identifying the relative contributions of these two factors could clarify the etiology of disc degeneration and point to novel strategies for assessing and/or enhancing the disc’s regenerative potential.

Advances in magnetic resonance imaging (MRI) permit quantitative, noninvasive assessment of the biochemical compositions of the CEP, vertebral bone marrow, and disc. Specifically, ultra-short echo-time (UTE) MRI enables measurement of CEP T2* relaxation times, which correlate with glycosaminoglycan (GAG) content, water content, and collagen-to-GAG ratios in the CEP [6]; chemical shift encoding-based water-fat MRI (CSE-MRI) enables assessment of bone marrow fat fraction (BMFF), which is inversely proportional to the amount of hematopoietic marrow and provides a proxy for perfusion [7-10]; and T1ρ MRI enables measurement of NP T1ρ relaxation times, which correlate with proteoglycan content and provide measurements of disc degeneration severity along a continuum [11]. These advanced MRI techniques have the potential to reveal factors influencing pathologic changes to the disc-endplate complex in vivo, but have not been previously combined in a single study. Our goal was to integrate these advanced quantitative MRI techniques into a comprehensive study of CEP and vertebral bone marrow composition in relation to disc degeneration in patients with cLBP. We hypothesized that deficits in CEP and vertebral bone marrow composition would associate with more severe disc degeneration after adjusting for patient age, sex, and body mass index (BMI).

Methods

Subjects

Between January 2020 and August 2021, 84 patients with cLBP were prospectively recruited and imaged for this IRB-approved study. Written informed consent was obtained from each patient. All patients met the criteria for cLBP as defined by the National Institutes of Health Pain Consortium Research Task Force [12], assessed using the following two questions: (1) “How long has low-back pain been an ongoing problem for you?” and (2) “How often has low-back pain been an ongoing problem for you over the past six months?”. A response of greater than three months to question one, and a response of “at least half the days in the past six months” to question two, met the cLBP criteria. Major exclusion criteria were a history of lumbar spinal surgery or fracture, known disc herniation, autoimmune disorders (including ankylosing spondylitis, rheumatoid arthritis, and psoriatic arthritis), or malignancy. All recruitment and imaging was conducted at a single institution. Patient-reported measures for disability and pain were collected using the Oswestry Disability Index (ODI) and Visual Analog Scale (VAS), respectively.

MRI

MRI of the lumbar spine was performed on a GE Discovery MR750 scanner (3.0-Tesla) using an 8-channel phased-array spine coil (GE Healthcare, Waukesha, WI). MRI data used for the purposes of this study were acquired using sagittal acquisitions from a multi-echo UTE Cones sequence [13], a CSE-MRI sequence [7], a T1ρ mapping sequence [11], and standard clinical fast spin-echo sequences with T1- and T2-weighting (Supplemental Content 1).

Image analysis

To assess CEP composition, all four CEPs at the L4/5 and L5/S1 levels were manually segmented on the first UTE echo by a single trained annotator using a custom segmentation tool (IDL 8.8). The hyperintense CEP signal from the UTE images guided segmentation (Fig. 1, top row). The same annotator re-segmented a subset of 20 CEPs from five patients a minimum of four weeks after initial segmentation to assess intra-rater reliability (intraclass correlation coefficient [ICC] for resulting T2* values = 0.91). After segmentation, T2* relaxation times were computed and the 3D voxel data were rotated into a standard coordinate system (MATLAB R2020b) before a kidney-bean-shaped template was applied [14]. This algorithmic approach enabled consistent identification of the central CEP region adjacent to the NP. To account for the sensitivity of CEP T2* values to CEP orientation within the MRI bore, only CEPs oriented within ± 15° of the ‘magic angle’ (54.7°) in the sagittal plane were included, since T2* values are artificially depressed at other angles [6]. Subjects with CEP damage precluding segmentation were also excluded.

Fig. 1

Representative MRI image, tissue segmentations, and biomarker visualization for the cartilage endplate (top row), vertebral bone marrow (middle row), and intervertebral disc (bottom row). CEP = cartilage endplate, BMFF = bone marrow fat fraction, NP = nucleus pulposus

To assess vertebral bone marrow composition, the L4, L5, and S1 vertebral bodies were segmented from five mid-sagittal slices of the CSE-MRI fat-fraction images using a neural network [15], applying manual corrections as needed (Fig. 1, middle row). The segmented 3D voxel data were rotated about the shape’s principal axes into a standard coordinate system. This approach enabled isolation of the hemi-vertebral region adjacent to each endplate.

To quantify disc degeneration, the lumbar discs were automatically segmented from five mid-sagittal T1ρ images using a neural network (Fig. 1, lower row) [16]. After segmentation, T1ρ relaxation times were calculated for all voxels within each disc before the discs were rotated into a standard coordinate system in order to isolate the sub-region adjacent to each endplate. The central 40% of the disc in the anterior–posterior direction was used to represent the NP region [14]. Disc Pfirrmann grade was also scored by a radiologist [T.L.] on T2-weighted images [17].

Outcomes and statistical methods

The primary outcomes were the mean CEP T2* value in the central CEP, the mean BMFF in the hemi-vertebral body adjacent to each CEP, and the mean T1ρ value in the NP region adjacent to each CEP. Univariate linear regression was used to assess independent relationships between NP T1ρ (dependent variable), CEP T2*, and vertebral BMFF (independent variables). Logistic regression was used to assess the relationship between mean NP T1ρ and Pfirrmann grade. A mixed effects multivariable linear regression model with a random intercept accounting for multiple observations (from multiple spinal levels) per subject as well as age, sex, and BMI was used to assess the combined relationship between NP T1ρ, CEP T2*, and vertebral BMFF. Spinal level, and the interaction between level and CEP T2*, were also included in the mixed effects model. Using the statistically significant predictors of NP T1ρ from the aforementioned multivariable model, a separate multivariable mixed effects linear regression model was created to compare the effects of CEP T2*, vertebral BMFF, and their interaction on NP T1ρ. Log transformations of the data were explored but not included since the untransformed data best fit the models based on assessment of residuals.

Finally, we used published regression data [6] relating CEP T2* with tissue compositional measures to estimate the distribution of GAG, water content, and collagen-to-GAG ratio in each CEP.

Statistical analyses were conducted in JMP Pro (16.0) and two-sided p < 0.05 was considered statistically significant.

Results

Complete MRI data were successfully acquired from 80/84 prospectively enrolled patients. Following exclusion based on CEP damage (20 subjects) and orientation (110 endplates), 130 endplates from 60 patients were included in the final analysis (age = 40.0 ± 11.9 years, 28 (47%) female and 32 (53%) male, Table 1). There was a broad range of disc degeneration in the included levels: NP T1ρ ranged from 33.9–114.7 ms with mean ± SD = 62.4 ± 18.0 ms; Pfirrmann grades ranged from I–IV. NP T1ρ was significantly and negatively correlated with Pfirrmann grade (p < 0.0001), confirming that lower NP T1ρ values reflected more severe disc degeneration overall.

Table 1

Patient demographic and clinical data from the n = 60 people included in the analysis

Characteristic	n = 60 people
Age (years)	40.0 ± 11.9 (19–65)
Female, male	28 (47%), 32 (53%)
Weight (kg)	76.1 ± 13.2 (50–116)
Height (cm)	171.2 ± 12.2 (122–193)
BMI (kg/m2)	26.1 ± 4.7 (18.8–45.7)
ODI	26.1 ± 14.6 (2–74)
VAS	6.5 ± 2.4 (1–11)
Pfirrmann grade
I	7 (10%)
II	30 (43%)
III	17 (25%)
IV	15 (22%)
V	0 (0%)

Data are presented as mean ± SD (min–max) or number (percent of total)

Results from the univariate regression analysis demonstrated that CEP T2* and vertebral BMFF were each independently significantly associated with NP T1ρ (p = 0.003 and 0.0001, respectively; Fig. 2). These explanatory variables, as well as NP T1ρ, were significantly correlated with age (Table 2); after adjusting for age, sex, and BMI, NP T1ρ remained significantly associated with CEP T2* (p = 0.0001, Table 3) but not vertebral BMFF (p = 0.43). There was a statistically significant effect of level (p < 0.0001), indicating more severe disc degeneration at L5/S1 compared with L4/5 in these patients (mean NP T1ρ = 72.1 ms [95% CI: 66.7, 77.5] vs. 60.5 ms [56.5, 64.5] for L4/5 vs. L5/1, p < 0.0001). There was also a statistically significant interaction between level and CEP T2* (p = 0.0006), indicating that the relationship between CEP composition and disc degeneration was different at L4/5 compared with L5/S1 (NP T1ρ was more sensitive to variations in CEP T2* at L4/5 than at L5/S1).

Fig. 2

Scatter plots showing the independent relationships between NP T1ρ and a vertebral BMFF and b CEP T2*

Table 2

Pearson correlation coefficient quantifying the cross-correlations between NP T1ρ and the continuous explanatory variables

	NP T1ρ	CEP T2*	BMFF	Age
CEP T2*	0.259††	–
BMFF	− 0.333†††	− 0.160	–
Age	− 0.381†††	− 0.185†	0.599†††	–
BMI	0.043	0.339†††	− 0.049	− 0.100

†p < 0.05

††p < 0.01

†††p < 0.0001

Table 3

Parameter estimates (βi) and 95% confidence intervals generated from a mixed effects linear regression model accounting for multiple spinal levels per subject predicting NP T1ρ (outcome)

Term	βi estimate	Lower 95%	Upper 95%	p-value
CEP T2*	1.06	0.53	1.59	0.0001†
Level	5.80	3.76	7.84	< .0001*
CEP T2* × Level	0.83	0.37	1.30	0.0006†
Age	− 0.45	− 0.83	− 0.06	0.0245†
Vertebral BMFF	− 0.13	− 0.47	0.20	0.431
Sex	2.32	− 1.89	6.53	0.275
BMI	− 0.26	− 1.15	0.64	0.567

† indicates two-sided p < 0.05

A second multivariable regression model included the significant predictors of NP T1ρ (age and level) plus CEP T2*, vertebral BMFF, and their interaction terms. Consistent with the first multivariable model, results from the second model showed that CEP T2*, age, and spinal level were significant predictors of NP T1ρ while vertebral BMFF was not (Table 4). The interaction between CEP T2* and vertebral BMFF was not significant (p = 0.73), indicating the relationship between NP T1ρ and CEP T2* did not depend on vertebral BMFF.

Table 4

Parameter estimates (βi) and 95% confidence intervals generated from a mixed effects linear regression model accounting for multiple spinal levels per subject predicting NP T1ρ (outcome)

Term	βi estimate	Lower 95%	Upper 95%	p-value
Level	5.88	3.31	8.46	< .0001*
CEP T2*	1.07	0.51	1.63	0.0003†
CEP T2* × Level	0.83	0.32	1.34	0.0018†
Age	− 0.40	− 0.79	− 0.02	0.0395†
Vertebral BMFF	− 0.19	− 0.62	0.24	0.384
Vertebral BMFF × Level	− 0.01	− 0.30	0.27	0.922
CEP T2* × Vertebral BMFF	− 0.01	− 0.07	0.05	0.726
CEP T2* × Vertebral BMFF × Level	0.00	− 0.05	0.05	0.931

† indicates two-sided p < 0.05

Based on the range of calculated CEP T2* values along with published regression equations relating CEP T2* values to biochemical measurements [6], CEP GAG content was estimated to range from 58–95 µg/mg dry weight, water content from 49–66%, and collagen-to-GAG ratio from 4.8–10.0.

Discussion

These results provide the first in vivo evidence that lumbar disc degeneration severity in patients with cLBP relates to compositional deficits in the CEP. We found that CEPs with lower T2* values—reflecting lower CEP hydration, lower GAG content, a higher ratio of collagen-to-GAG, and lower permeability [4, 6]—were associated with more severe disc degeneration both before and after accounting for patient age, sex, BMI, and multiple measurements per subject. Conversely, higher BMFF values in the adjacent vertebral body—reflecting less hematopoietic marrow and decreased perfusion [7–10, 18]—were not associated with disc degeneration after accounting for covariates (namely, subject age). Taken together, these collective findings suggest that poor CEP composition may play a significant role in disc degeneration severity in patients with cLBP, and that CEP composition can affect disc health both with and without deficits in vertebral perfusion.

CEP composition and vertebral vascularity are important factors in disc health, though their relative importance has remained unclear because prior studies have focused on each factor alone. For example, in vitro experiments have demonstrated that variations in the amounts of aggrecan, mineral, and collagen in the CEP can impede nutrient diffusion to a degree sufficient to hinder NP cell viability and function [4, 19]. Experiments using diffusion chambers have also shown that limiting the reservoir of available nutrients can reduce NP cell viability [5]. Those prior findings demonstrate the potential mechanisms by which CEP composition and vertebral vascularity may affect NP cells; our present findings demonstrate the clinical importance of CEP composition in the context of variations in vertebral vascularity. Specifically, unlike vertebral BMFF, CEP composition was significantly associated with disc health after accounting for co-variates, which demonstrates the relative importance of poor CEP composition in disc degeneration in patients with cLBP.

These findings help clarify mechanistic factors influencing disc degeneration in a cohort of patients with non-specific cLBP, but they do not imply that CEP composition is uniquely important to patients with cLBP relative to asymptomatic subjects with disc degeneration, or that disc degeneration was the pain generator in this cohort. Since we did not study asymptomatic subjects, it remains unclear whether the relationships discovered here in cLBP patients are the same in asymptomatic subjects. Likewise, factors differentiating painful disc degeneration from asymptomatic disc degeneration as seen on MRI also remain unclear. These are important areas for future research.

Our finding that disc health strongly associates with CEP composition has implications for the clinical translation of emerging biologic therapies designed to regenerate the disc. The therapeutic potential of cell-based regenerative therapies likely depends on disc nutrient supply: mesenchymal stem cell chondrocytic differentiation, proliferation, and function depend on glucose and oxygen concentrations (the same nutrients vital to NP cell survival) and on matrix acidity (related to lactate production, an NP cellular metabolic waste product) [20]. Given the established effects of CEP composition on nutrient transport in combination with the data shown here, it follows that the efficacy of cell-based therapies may likewise depend on CEP composition. If so, CEP composition could be a possible diagnostic target to help identify patients with adequate CEP permeability to support increased cell density or metabolic activity in the NP. Similarly, poor CEP composition could also be a possible therapeutic target to enhance the disc’s regenerative potential, e.g., through enzymatic augmentation [21]. Non-invasive assessment of CEP composition with UTE MRI may thus provide a powerful tool to improve selection of patients and spinal levels likely to respond to anabolic treatment, helping to maximize the efficacy of novel biologic therapies.

The variations in biomarker values with age reported here are consistent with established age-associated biochemical changes in these tissues [18, 22, 23]. Lumbar CEP GAG and water contents in 51–67 year old subjects ranged from 43–185 µg/mg dry weight and 22–62%, respectively, in a previous study [24]. For patients with similar ages in the current study (49–65 year old), the estimated GAG and water content ranges (65–84 µg/mg dry weight and 52–63% water) are within the range previously published. A prior study using CSE-MRI [8] found that higher levels of BMFF were associated with more severe disc degeneration, consistent with the univariate regression results shown here. These findings collectively support the validity of our approach.

One strength of our study is the utilization of two recently developed neural networks [15, 16] to perform automated tissue segmentation. This enabled rapid and unbiased evaluation of the MRI biomarkers, demonstrating the relevance of neural networks to a diverse class of problems and addressing the inefficiency and subjectivity associated with manual tissue segmentation. Notably, we also demonstrated the feasibility of non-invasively assessing CEP composition with a new MRI biomarker (mean CEP T2*) in a clinical population.

This study has several limitations. First, the cross-sectional study design precludes establishment of causality. However, prior research has demonstrated the causal effects of these same CEP compositional traits on NP cell survival and function [4], which provides a mechanistic explanation for our results. A second limitation relates to the spatial resolution of the vertebral BMFF maps (voxel dimensions 1.0 × 1.0 × 4.0 mm), which does not capture subtle morphologic variations in endplate-specific vascularity which may affect nutrient supply [25]. However, the conversion from hematopoietic to fatty marrow has been shown to associate with vascular changes in both the femur and vertebral body and with age-associated changes in perfusion [9, 10, 18], suggesting that vertebral BMFF values likely provide an accurate measure of overall changes in vertebral vascularity.

In summary, we found that after adjusting for age, sex, and BMI, deficits in CEP composition—as indicated by lower T2* values—were associated with more severe disc degeneration. Conversely, lower vertebral vascularity—as indicated by higher vertebral BMFF values—was not associated with disc degeneration after adjusting for age. We conclude that poor CEP composition plays a significant role in disc degeneration severity in patients with cLBP both with and without deficits in vertebral vascularity.

Below is the link to the electronic supplementary material.

Supplementary file1 (DOCX 23 KB)