The magnetic resonance imaging-based vertebral bone quality scoring system: A novel method to evaluate endplate changes in patients with primary single-level disk herniation and Modic changes

Abstract

Objectives: This study aimed to investigate differences in vertebral fat distribution and bone density between patients with and without Modic changes (MCs) using a magnetic resonance imaging (MRI)-based vertebral bone quality (VBQ) scoring system.

Patients and methods: In this retrospective study, 189 patients (95 males, 94 females; mean age: 54±2.2 years; range, 18 to 82 years) with primary single-level disk herniation were reviewed between June 2021 and June 2022. The patients were divided into the MC group (n=99) and the non-MC (NMC) group (n=90). The subcutaneous fat tissue thickness and bone mineral density were determined. The system consisted of two scores: the VBQ score, which reflected the fatty infiltration within the vertebral body, and the endplate bone quality (EBQ) score, which reflected the signal intensity (SI) of the upper and lower endplates. The EBQ score is a novel measurement that we introduced in this study. The VBQ and EBQ were measured and scored using MRI scans. The mean SI of the upper and lower endplates (endplate SI)/the bone marrow SI (marrow SI) was measured.

Results: There was a considerable difference in subcutaneous fat tissue thickness between the MC and NMC groups (1.40 vs. 1.16 cm, p=0.01). The EBQ scores of the L4 and L5 vertebrae and endplate SI/marrow SI of all vertebral body levels were significantly higher in the MC group.

Conclusion: The occurrence of MCs in the lumbar spine may be associated with abnormal fat distribution. The distribution of vertebral fat in patients with MCs is distributed earlier in the upper and lower endplates of the vertebral body, and this trend is not observed in patients without MC. The thickness of subcutaneous fat tissue is a key factor in the occurrence of MCs.

Introduction

The current management of low back pain (LBP) offers low to moderate improvement in pain and disability. Kjaer et al.[1] suggested that Modic changes (MCs) constitute the crucial element in the degenerative process around the disk in relation to LBP and clinical findings. Modic changes in the lumbar spine were first reported by de Roos et al.[2] in 1987. Modic et al.[3] systematically described the types of signal changes, classification criteria, and pathological changes in 1988. Modic changes are classified into three types: type 1 lesions indicate that the vascularized fibrous tissue in the bone marrow is undergoing an active degenerative process; type 2 lesions reflect fat replacement in the bone marrow and are thought to be stable with chronic fat infiltration and deposition; type 3 lesions are believed to be associated with endplate sclerosis.[3,4] Modic changes are strongly associated with disk degeneration and LBP.[5,6] However, the etiology and pathological features of MCs are not well understood and cannot explain why they occur in only a small proportion of patients with disk degeneration.

While exploring the etiology and pathogenesis of MCs, the authors found that abnormal fat metabolism may play a role in many aspects. First, abnormal fat metabolism can lead to obesity, a relatively heavy spinal load, disk degeneration, and an increased risk of endplate injury. Özcan-Ekşi et al.[7] reported that subcutaneous fat tissue thickness (SFTT) is a good predictor of severe intervertebral disk degeneration (IVDD) and MCs. Second, if the adipose tissue in the bone marrow of the vertebral body increases, it can activate some signaling pathways and inflammatory factors. Additionally, the response of myeloid tissue to inflammatory stimulation may change the structure of the bone marrow or local microenvironment, resulting in MCs in the vertebral body.[8,9]

Magnetic resonance imaging (MRI) is often used in the evaluation of patients with LBP.[10] Furthermore, MRI has a role in detecting osteoporosis, and the appearance of the bone marrow is determined by the relative amount of fat, protein, water, and bone cells on the MRI sequence. Previous studies have revealed that an additional sequence of MRI could provide additional information on the bone marrow composition.[11-13] The vertebral bone quality (VBQ) score is a popular new technique for assessing bone quality by determining fatty infiltration within the vertebral body.[14] Modic changes often occur at the vertebral endplate; therefore, the endplate bone quality (EBQ) score was introduced in this study.[15]

The signal intensity (SI) of the endplate and the vertebral body on an MRI image can reflect the degree of fat infiltration or inflammation, which may affect the bone quality and the occurrence of MCs. Therefore, the authors developed an MRI-based VBQ scoring system consisting of two scores: the VBQ score, which reflects the fatty infiltration within the vertebral body, and the EBQ score, which reflects the SI of the upper and lower endplates. The EBQ score is a novel measurement that was introduced in this study to evaluate the endplate changes associated with MCs.

This study aimed to investigate the differences in vertebral fat distribution and bone density between patients with and without MCs using the MRIbased VBQ scoring system. We hypothesized that patients with MCs would have more fat infiltration in the endplates and lower bone density than patients without MCs.

Patients and Methods

A total of 189 patients (95 males, 94 females; mean age: 54±2.2 years; range, 18 to 82 years) with primary single-segment disk herniation who were hospitalized in the Third Hospital of Hebei Medical University between June 2021 and June 2022 were retrospectively reviewed. The disk degeneration was assessed according to the Pfirrmann classification (Additional Figure 1). The patients were divided into the MC group (n=99) and the non-MC (NMC) group (n=90). The inclusion criteria for the MC group were the presence of any type (1, 2, or 3) or range (A, B, or C) of MCs on the MRI images, as described by Modic et al.[3] The MC range was classified into three categories according to the percentage of the endplate area involved in MCs: category A represented <25%, B represented 25-50%, and C represented >50%. The independent variables of this study were the presence or absence of MCs, the SFTT, and the blood lipid levels. The outcome parameters of this study were the VBQ and EBQ scores, the bone mineral density (BMD), and the mean SI of the upper and lower endplates (endplate SI) divided by the bone marrow SI (marrow SI). Subcutaneous fat tissue thickness, triglycerides, high-density lipoprotein, low-density lipoprotein, and total cholesterol levels were measured in all patients. The MRI images and BMD measurements were retrospectively analyzed. The exclusion criteriawere vertebral fracture, previous lumbar spine surgery, spinal infection, and severe spinal deformity.

ADDITIONAL FIGURE 1. Algorithm used for grading lumbar disk degeneration.

The MRI of the lumbar spine was performed with a Magnetom Avanto 1.5T (Siemens Healthineers, Erlangen, Germany) superconducting scanner. The Pfirrmann grade, SFTT, VBQ, and EBQ were obtained from the MRI data. Image evaluation was performed by an experienced radiologist and an orthopedic surgeon who were unaware of the patients’ information (T1 sagittal position: repetition time/echo time=546/9.7 msec, section thickness=4 mm; T2 sagittal position: repetition time/echo time=2,850/103 msec, section thickness=4 mm). The evaluation included an MRI scan image showing the classification and location of the MCs.[16] In addition, the degree of IVDD was assessed based on the Pfirrmann grade on the MRI T2-weighted images.[17]

The lumbar and hip dual-energy X-ray absorptiometry BMD T-scores were recorded. The BMD assessment criteria were as follows: osteopenia (T-score between −1.0 and −2.5 standard deviation [SD]) or osteoporosis (T-score <-2.5 SD) was the lowest T-score residing in the femoral or vertebral region.[17-19]

In axial T2-weighted lumbar MRI, the SFTT was measured as the vertical distance from the spinous process tip of the L1 vertebral body to the skin (Figure 1). Two authors defined this method; this technique was interpreted by the raters in a substantially consistent manner.[7]

FIGURE 1. Measurement technique of subcutaneous fat tissue thickness at the L1/L2 level on T2-weighted axial lumbar magnetic resonance imagingb

The VBQ and EBQ scores were assessed with lumbar T1-weighted MRI.[14,15] The vertebral body was excluded if the entire vertebral body was abnormal. For example, if the cerebrospinal fluid (CSF) in the L3 space was completely blocked, the region of interest (ROI) of the CSF was placed at the L2 or S1 level. The VBQ score was calculated by dividing the median SI at L1-L4 levels by the SI of the CSF.

The VBQ, EBQ, and SFTT, observed and measured by two independent researchers, were the mean values of the range and location of MCs, as well as the grading and measurement of the IVDD in statistical analysis. To assess the concordance of the various measurements between the raters, 50 patients were randomly selected, observed, and measured independently by two experienced spine surgeons.[20] Three weeks later, the same observations and measurements were made by each observer.

The VBQ score was calculated by drawing an ROI around the vertebral body on the sagittal T1-weighted MRI image, excluding the endplates and the posterior elements. Afterward, the median SI of the ROI was measured, which represented the SI of the vertebral body. Next, another ROI was drawn around the CSF on the same image, and the median SI of the ROI was measured, which represented the SI of the CSF. Finally, the SI of the vertebral body was divided by the SI of the CSF value to obtain the VBQ score. A higher VBQ score indicated a higher degree of fatty infiltration within the vertebral body (Figure 2).

FIGURE 2. Representative image of the region of interests used to calculate the vertebral bone quality score.

The EBQ score was calculated by drawing an ROI of 3 mm in width along the upper and lower endplates on the sagittal T1-weighted MRI image, excluding the cortical bone. Afterward, the median SI of the ROI was measured, which represented the SI of the endplate. Next, another ROI was drawn around the center of the vertebral body on the same image, and the median SI of the ROI was measured, which represented the SI of the bone marrow. Finally, the SI of the endplate value was divided by the SI of the bone marrow to obtain the EBQ score. A higher EBQ score indicated a higher SI of the endplate, which could reflect the presence of MCs (Figure 3).

FIGURE 3. Representative image of the region of interests used to calculate the endplate bone quality score.

Statistical analysis

The sample size for the study was determined based on a power analysis using G*Power software (Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany). A medium effect size (f=0.25), an alpha level of 0.05, and a power of 0.80 were used for the independent sample t-tests. The required sample size was estimated to be 84 patients in each group. However, a total of 189 patients were enrolled to account for possible dropouts and missing data.

Data were analyzed using IBM SPSS version 26.0 software (IBM Corp., Armonk, NY, USA). The intragroup correlation coefficient was calculated to test the reliability among and within raters. The Shapiro-Wilk test was used to assess the normality of continuous variables. Independent-sample t-tests were used for normally distributed variables and MannWhitney U tests for nonnormally distributed data. Normally distributed variables were presented as the mean ± standard deviation (SD), and nonnormally distributed variables were presented as median (range) values. The vertebral endplates and bone marrow VBQ scores were generally distributed in each group. The interobserver reliability was assessed using the intragroup correlation coefficient. The differences in VBQ and EBQ scores between the MC and NMC groups were tested by independent sample t-tests. Signal intensity differences between the upper and lower endplates of the target vertebral body and the vertebral center were analyzed by independent sample t-tests. A p-value <0.05 was considered statistically significant.

Results

There was no significant difference in sex, age, triglycerides, cholesterol, high-density lipoprotein, low-density lipoprotein, or VBQ between the two groups (p>0.05) (Table I). There was no difference in BMD between the two groups. The SFTT was higher in the MC group than in the NMC group (1.40 vs. 1.16 cm, p=0.01). Similarly, there was a significant difference in EBQ scores for the L4 and L5 vertebrae between the MC and NMC groups (L4 EBQ: 2.98±0.65 vs. 2.73±0.88, p<0.05; L5 EBQ: 3.21±1.53 vs. 2.84±0.90, p<0.05; Table II). The remaining vertebral EBQ scores were not significantly different (L1 EBQ: 2.9176±0.65 vs. 2.7970±0.91; L2 EBQ: 2.9627±0.63 vs. 2.8670±0.91; L3EBQ: 2.94±0.62 vs. 2.7676±0.92). In the MC group, the SI value within 3 mm of the endplate was higher than that at the center of the vertebral body. In contrast, in the NMC group, the SI value within 3 mm of the endplate was smaller than that at the center of the vertebral body. Therefore, the endplate SI/marrow SI of the vertebral body of each lumbar level was significantly different. The SI of the MC group was higher than that of the NMC group (Table III).

The distribution of Modic types and ranges in the MC group were compared, as shown in Tables IV and V. The majority of the MC group had type 2 MCs (87.9%), followed by type 1 (5.1%), type 3 (3%), and mixed type (4%). The most common Modic range was A (46.5%), followed by B (32.3%) and C (21.2%). The level of disk herniation and the proportion of severe IVDD between the MC and NMC groups were also compared, as shown in Tables IV and V. The most common level of disk herniation was L4/L5 in both groups (70.7% in the MC group and 76.7% in the NMC group), followed by L5/S1 (15.2% in the MC group and 15.6% in the NMC group). The proportion of severe IVDD was higher in the MC group than in the NMC group at all levels, except for L3/L4, where it was equal (57% in both groups).

In the MC group, there were five (5.1%) type 1 cases, 87 (87.9%) type 2 cases, three (3%) type 3 cases, and four (4%) mixed type cases (Figure 4). For the MC range, there were 46 cases of A, 32 cases of B, and 21 cases of C (Figure 5). At the level of disk herniation, the MCs mainly occurred at the L4/L5 level (70.7%) and L5/S1 level (15.15%).

FIGURE 4. The type of modic change.

FIGURE 5. The range of modic change.

In terms of disk herniation location, within the MC group, the majority of cases occurred at the L4/L5 level, accounting for 70.7% of cases, while the L5/S1 level was the second most common site with 15.15%. Similarly, in the NMC group, disk herniation was predominantly found at the L4/L5 level, representing 76.66% of cases, with the L5/S1 level accounting for 15.55% of cases. This indicates that both groups experienced lumbar disk herniation and disk degeneration primarily in the L4/L5 and L5/ S1 regions.

In the NMC group, the level of disk herniation was L4/5 in 79 (87%) patients and L5/S1 in 62 (68%) patients. In the MC group, there were 91 (91%) L4/L5 segments and 80 (80%) L5/S1 segments. As shown in Figure 5, the distribution of MCs was mainly concentrated in the L4/L5 and L5/S1 segments (Figure 6).

FIGURE 6. The distribution of MCs. MCs: Modic changes.

For all the above measurements, such as the VBQ score, EBQ score, and SFTT, as well as the grading of Modic change range, location, and Pfirrmann classification, the inter- and intrarater reliability coefficients were above 0.8.

Discussion

This study is the first to report the VBQ and EBQ scores of patients with MCs.[21] The results showed that the difference in EBQ between the two groups may be related to abnormal fat distribution in the vertebral body. The distribution of vertebral fat in patients with MCs is distributed earlier in the upper and lower endplates of the vertebral body, and this trend is not observed in patients without MCs. The EBQ and VBQ scores reflect the SI of the endplate and the vertebral body on MRI images, which may indicate the degree of fat infiltration or inflammation. A higher EBQ score means a higher SI of the endplate, which may suggest the presence of MCs. A lower EBQ score means a lower SI of the endplate, which may indicate a normal or sclerotic endplate. The EBQ score is different from the SI, which is a quantitative measure of the brightness of the image. The EBQ score is a ratio of the SI of the endplate to the SI of the bone marrow, which is a relative measure of the contrast between the endplate and the vertebral body. In this study, the most common level of disk herniation and MC was L4/L5, followed by L5/S1. This is in agreement with previous studies that reported a higher prevalence of MCs in the lower lumbar segments.[22] This may be explained by the higher mechanical load and stress that these segments experience, which may lead to more damage and inflammation in the endplates and disks.

Studies[23] report that endplate sclerosis occurs in all three Modic subtypes. In addition, mixed MCs were observed. Modic changes were generally more frequent in patients with LBP (43%) than in those without LBP (6%). We did not analyze the differences among Modic subtypes since we considered MCs to be a dynamic variable.

The segments in which MCs most commonly occurred in this retrospective study were concentrated at L4/L5 and L5/S1, which is consistent with findings reported in the literature.[24,25] The incidence of severe IVDD at L4/L5 and L5/S1 was 91% and 80% in the MC group and 87% and 68% in the NMC group, respectively. This is broadly consistent with reports in the literature.[26]

Albert et al.[27] found that when the lumbar disk degenerates, the stress on the functional unit of the spine redistributes, and the load on the endplate increases, leading to microfracture of the endplate. Recurrent microfractures of the endplate may be the leading cause of MC formation. Abnormal fat metabolism can lead to obesity, which aggravates stress in the functional units of the spine, and we can assess obesity and obesity-related health problems by measuring the SFTT at the upper lumbar level. Özcan-Ekşi et al.[4] reported that SFTT is a good predictor of severe IVDD and MCs. According to the literature, SFTT represents body fat content as well as being overweight.[28] This is more representative of patient weight than body mass index, with a higher risk of disk degeneration and MCs when SFTT >0.9 cm. The mean values in both groups were >0.9 cm since the inclusion criteria for this study were patients with LBP requiring surgery and severe disk degeneration. According to the results, as the thickness of subcutaneous fat increases, the occurrence of MCs also increases. In one study, 58% of endplates were ruptured at the time of 17% IVDD, suggesting that endplates are a weaker site in the functional unit of the spine.[29] In the same patients with severe IVDD and endplate microfracture, there were also many patients without MCs in the vertebral body; this phenomenon may be related to the response of vertebral bone marrow adipose tissue to inflammatory stimuli released by the intervertebral disk.

According to a previous study, a large number of cytokines, such as tumor necrosis factors, interleukin-6, and interleukin-8, accumulate in disk tissue with MCs.[8] These proinflammatory cytokines and injury-related molecules can penetrate from the intervertebral disk to the bone marrow through the endplate injury gap and react with the bone marrow of the vertebral body, resulting in changes in the local microenvironment and structure of the bone marrow of the vertebral body. Currently, fatty tissue is considered to be a highly active metabolic structure that contributes to systemic inflammation and the secretion of large amounts of proinflammatory cytokines.[30] In addition, it was reported that an increase in very-low-density lipoprotein diameter promotes the development of MCs.[31] Furthermore, an increase in blood lipids and glucose has been associated with the occurrence of MCs.[32,33]

The response of the bone marrow often depends on the composition of the bone marrow itself; if there is more fatty tissue in the bone marrow of the vertebral body, it can activate some signaling pathways and inflammatory factors. As a result, the local inflammatory reaction of the bone marrow is enhanced, resulting in alterations to the structure or regional microenvironment of the bone marrow and MCs in the vertebral body. We analyzed blood lipids from the two groups. Although none of the blood lipids were significantly different, the mean values of patients with MCs were larger than those of patients without MCs. Therefore, we tried to identify differences in vertebral fat infiltration between the two groups by measuring data using imaging-related techniques.

There was no significant difference in VBQ between the MC and NMC groups; the MCs mainly occurred near the endplate. The EBQ was used to measure the ROI of 3 mm in the upper and lower endplates; there were no significant differences in the other vertebral bodies.

Interestingly, the SI of the upper and lower endplates of each vertebral body in patients with MCs was more significant than the SI in the center of the vertebral body. In contrast, the opposite was observed in patients without MCs. Therefore, we compared the ratio of the ROI. The ratio of the ROI was the ratio of the SI in the endplate 3 mm to the central SI in the vertebral body; that is, the endplate SI/bone marrow SI in the vertebral body. The ratio of each vertebral body was significantly different, according to a t-test. Therefore, it can be concluded that the SI of the upper and lower endplates is higher in patients with MC, indicating that fat infiltration in the endplate is more evident than in NMC patients. According to relevant studies, fat distribution in the vertebral body is high in the middle and low on both sides.[34-36] Additionally, it is common in the middle, high on both sides, and can occur with aging.[37] Therefore, the abnormal distribution and metabolism of fat will occur earlier in patients with MCs. Excessive fatty tissue and inflammatory cells infiltrate the upper and lower endplates of the vertebral body. Moreover, local arteriosclerosis of small blood vessels will appear earlier and cause fat embolism, leading to earlier degeneration of the intervertebral disk.

This study has some limitations. First, this is a single-center retrospective study. Therefore, the sample size needed to be larger. The number of cases of type 1 and 3 MCs were inadequate. Second, the study focused on the phenotypic analysis of endplate changes on MRI scans rather than the study-related histopathological and mechanical analyses. Third, there were errors in SI value measurements due to subjective measures. Further studies may be needed in the future.

In conclusion, the occurrence of MCs in the lumbar spine may be associated with abnormal fat distribution, particularly in the vertebral body. The distribution of vertebral fat in patients with MCs is distributed earlier in the upper and lower endplates of the vertebral body, whereas this trend is not significant in patients without MCs. The thickness of subcutaneous fat tissue is a key factor in the occurrence of MCs.

* The two authors contributed equally to this study.

Citation: Li C, Gao W, Yao X, Tong T, Wang Y, Li W, et al. The magnetic resonance imaging-based vertebral bone quality scoring system: A novel method to evaluate endplate changes in patients with primary single-level disk herniation and Modic changes. Jt Dis Relat Surg 2024;35(2):257-266. doi: 10.52312/ jdrs.2024.1527.

The study protocol was approved by the Third Hospital of Hebei Medical University Ethics Committee (date: 25.10.2022, no: Ke 2022-115-1). The study was conducted in accordance with the principles of the Declaration of Helsinki.

A written informed consent was obtained from each patient.

Conception and design: C.J.L., W.S.G.; Administrative support: X.W.Y., T.T.; Provision of study materials or patients: Y.S.W., J.L.A.; Collection and assembly of data: L.F.W., C.J.L., W.S.G.; Data analysis and interpretation: X.W.Y., J.L.A.; Manuscript writing and final approval of manuscript: All authors.

The authors declared no conflicts of interest with respect to the authorship and/or publication of this article.

The authors received no financial support for the research and/or authorship of this article.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

References

1. Kjaer P, Korsholm L, Bendix T, Sorensen JS, Leboeuf-Yde C. Modic changes and their associations with clinical findings. Eur Spine J 2006;15:1312-9. doi: 10.1007/s00586-006-0185-x.
2. de Roos A, Kressel H, Spritzer C, Dalinka M. MR imaging of marrow changes adjacent to end plates in degenerative lumbar disk disease. AJR Am J Roentgenol 1987;149:531-4. doi: 10.2214/ajr.149.3.531.
3. Modic MT, Steinberg PM, Ross JS, Masaryk TJ, Carter JR. Degenerative disk disease: Assessment of changes in vertebral body marrow with MR imaging. Radiology 1988;166:193-9. doi: 10.1148/radiology.166.1.3336678.
4. Özcan-Ekşi EE, Yayla A, Orhun Ö, Turgut VU, Arslan HN, Ekşi MŞ. Is the distribution pattern of modic changes in vertebral end-plates associated with the severity of intervertebral disc degeneration?: A cross-sectional analysis of 527 Caucasians. World Neurosurg 2021;150:e298-e304. doi: 10.1016/j.wneu.2021.02.128.
5. Albert HB, Manniche C. Modic changes following lumbar disc herniation. Eur Spine J 2007;16:977-82. doi: 10.1007/ s00586-007-0336-8.
6. Mok FP, Samartzis D, Karppinen J, Fong DY, Luk KD, Cheung KM. Modic changes of the lumbar spine: Prevalence, risk factors, and association with disc degeneration and low back pain in a large-scale population-based cohort. Spine J 2016;16:32-41. doi: 10.1016/j.spinee.2015.09.060.
7. Özcan-Ekşi EE, Kara M, Berikol G, Orhun Ö, Turgut VU, Ekşi MŞ. A new radiological index for the assessment of higher body fat status and lumbar spine degeneration. Skeletal Radiol 2022;51:1261-71. doi: 10.1007/s00256-021-03957-8.
8. Torkki M, Majuri ML, Wolff H, Koskelainen T, Haapea M, Niinimäki J, et al. Osteoclast activators are elevated in intervertebral disks with Modic changes among patients operated for herniated nucleus pulposus. Eur Spine J 2016;25:207-16. doi: 10.1007/s00586-015-3897-y.
9. da Cruz Fernandes IM, Pinto RZ, Ferreira P, Lira FS. Low back pain, obesity, and inflammatory markers: Exercise as potential treatment. J Exerc Rehabil 2018;14:168-74. doi: 10.12965/jer.1836070.035.
10. Deyo RA, Weinstein JN. Low back pain. N Engl J Med 2001;344:363-70. doi: 10.1056/NEJM200102013440508.
11. Manenti G, Capuani S, Fusco A, Fanucci E, Tarantino U, Simonetti G. Osteoporosis detection by 3T diffusion tensor imaging and MRI spectroscopy in women older than 60 years. Aging Clin Exp Res 2013;25:S31-4. doi: 10.1007/ s40520-013-0091-0.
12. Qiu X, Fu Y, Chen J, Ye Y, Wang Z, Ming X. The correlation between osteoporosis and blood circulation function based on magnetic resonance ımaging. J Med Syst 2019;43:91. doi: 10.1007/s10916-019-1206-8.
13. Li GW, Xu Z, Chen QW, Tian YN, Wang XY, Zhou L, et al. Quantitative evaluation of vertebral marrow adipose tissue in postmenopausal female using MRI chemical shiftbased water-fat separation. Clin Radiol 2014;69:254-62. doi: 10.1016/j.crad.2013.10.005.
14. Ehresman J, Schilling A, Pennington Z, Gui C, Chen X, Lubelski D, et al. A novel MRI-based score assessing trabecular bone quality to predict vertebral compression fractures in patients with spinal metastasis. J Neurosurg Spine 2019:1-8. doi: 10.3171/2019.9.SPINE19954.
15. Jones C, Okano I, Arzani A, Dodo Y, Moser M, Reisener MJ, et al. The predictive value of a novel site-specific MRI-based bone quality assessment, endplate bone quality (EBQ), for severe cage subsidence among patients undergoing standalone lateral lumbar interbody fusion. Spine J 2022;22:1875-83. doi: 10.1016/j.spinee.2022.07.085.
16. Udby PM, Samartzis D, Carreon LY, Andersen MØ, Karppinen J, Modic M. A definition and clinical grading of Modic changes. J Orthop Res 2022;40:301-7. doi: 10.1002/ jor.25240.
17. Pfirrmann CW, Metzdorf A, Zanetti M, Hodler J, Boos N. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine (Phila Pa 1976) 2001;26:1873-8. doi: 10.1097/00007632-200109010-00011.
18. Kanis JA. Assessment of fracture risk and its application to screening for postmenopausal osteoporosis: Synopsis of a WHO report. WHO Study Group. Osteoporos Int 1994;4:368-81. doi: 10.1007/BF01622200.
19. Kanis JA, Melton LJ 3rd, Christiansen C, Johnston CC, Khaltaev N. The diagnosis of osteoporosis. J Bone Miner Res 1994;9:1137-41. doi: 10.1002/jbmr.5650090802.
20. Salzmann SN, Okano I, Jones C, Zhu J, Lu S, Onyekwere I, et al. Preoperative MRI-based vertebral bone quality (VBQ) score assessment in patients undergoing lumbar spinal fusion. Spine J 2022;22:1301-8. doi: 10.1016/j. spinee.2022.03.006.
21. Atik OŞ. Writing for Joint Diseases and Related Surgery (JDRS): There is something new and interesting in this article! Jt Dis Relat Surg 2023;34:533. doi: 10.52312/ jdrs.2023.57916.
22. Wang Y, Videman T, Battié MC. Modic changes: Prevalence, distribution patterns, and association with age in white men. Spine J 2012;12:411-6. doi: 10.1016/j. spinee.2012.03.026.
23. Jensen TS, Karppinen J, Sorensen JS, Niinimäki J, LeboeufYde C. Vertebral endplate signal changes (Modic change): A systematic literature review of prevalence and association with non-specific low back pain. Eur Spine J 2008;17:1407-22. doi: 10.1007/s00586-008-0770-2.
24. Weishaupt D, Zanetti M, Hodler J, Min K, Fuchs B, Pfirrmann CW, et al. Painful lumbar disk derangement: Relevance of endplate abnormalities at MR imaging. Radiology 2001;218:420-7. doi: 10.1148/radiology.218.2.r01fe15420.
25. Karchevsky M, Schweitzer ME, Carrino JA, Zoga A, Montgomery D, Parker L. Reactive endplate marrow changes: A systematic morphologic and epidemiologic evaluation. Skeletal Radiol 2005;34:125-9. doi: 10.1007/ s00256-004-0886-3.
26. Teichtahl AJ, Urquhart DM, Wang Y, Wluka AE, O’Sullivan R, Jones G, et al. Lumbar disc degeneration is associated with modic change and high paraspinal fat content - a 3.0T magnetic resonance imaging study. BMC Musculoskelet Disord 2016;17:439. doi: 10.1186/s12891-016-1297-z.
27. Albert HB, Kjaer P, Jensen TS, Sorensen JS, Bendix T, Manniche C. Modic changes, possible causes and relation to low back pain. Med Hypotheses 2008;70:361-8. doi: 10.1016/j. mehy.2007.05.014.
28. Berikol G, Ekşi MŞ, Aydın L, Börekci A, Özcan-Ekşi EE. Subcutaneous fat index: A reliable tool for lumbar spine studies. Eur Radiol 2022;32:6504-13. doi: 10.1007/s00330-022-08775-7.
29. Hilton RC, Ball J. Vertebral rim lesions in the dorsolumbar spine. Ann Rheum Dis 1984;43:302-7. doi: 10.1136/ ard.43.2.302.
30. Tilg H, Moschen AR. Adipocytokines: Mediators linking adipose tissue, inflammation and immunity. Nat Rev Immunol 2006;6:772-83. doi: 10.1038/nri1937.
31. Li Y, Karppinen J, Cheah KSE, Chan D, Sham PC, Samartzis D. Integrative analysis of metabolomic, genomic, and imaging-based phenotypes identify very-low-density lipoprotein as a potential risk factor for lumbar Modic changes. Eur Spine J 2022;31:735-45. doi: 10.1007/s00586-021- 06995-x.
32. Hotamisligil GS. Inflammation, metaflammation and immunometabolic disorders. Nature 2017;542:177-85. doi: 10.1038/nature21363.
33. Blaser H, Dostert C, Mak TW, Brenner D. TNF and ROS crosstalk in inflammation. Trends Cell Biol 2016;26:249-61. doi: 10.1016/j.tcb.2015.12.002.
34. Bandirali M, Di Leo G, Papini GD, Messina C, Sconfienza LM, Ulivieri FM, et al. A new diagnostic score to detect osteoporosis in patients undergoing lumbar spine MRI. Eur Radiol 2015;25:2951-9. doi: 10.1007/s00330-015- 3699-y.
35. García-Cañas R, Rodríguez-Moro C. Increased vertebral bone marrow fat content can be associated with vertebral fractures and back pain. Ann Transl Med 2022;10:850. doi: 10.21037/atm-22-3827.
36. He J, Fang H, Li X. Vertebral bone marrow fat content in normal adults with varying bone densities at 3T magnetic resonance imaging. Acta Radiol 2019;60:509-515. doi: 10.1177/0284185118786073.
37. Silva BC, Broy SB, Boutroy S, Schousboe JT, Shepherd JA, Leslie WD. Fracture risk prediction by non-BMD DXA measures: The 2015 ISCD Official Positions Part 2: Trabecular Bone Score. J Clin Densitom 2015;18:309-30. doi: 10.1016/j.jocd.2015.06.008.