Skip to main content

Vertebral bone quality score was associated with paraspinal muscles fat infiltration, but not modic classification in patients with chronic low back pain: a prospective cross-sectional study

Abstract

Background

The lumbar vertebra and paraspinal muscles play an important role in maintaining the stability of the lumbar spine. Therefore, the aim of this study was to investigate the relationship between paraspinal muscles fat infiltration and vertebral body related changes [vertebral bone quality (VBQ) score and Modic changes (MCs)] in patients with chronic low back pain (CLBP).

Methods

Patients with CLBP were prospectively collected in four hospitals and all patients underwent 3.0T magnetic resonance scanning. Basic clinical information was collected, including age, sex, course of disease (COD), and body mass index (BMI). MCs were divided into 3 types based on their signal intensity on T1 and T2-weighted imaging. VBQ was obtained by midsagittal T1-weighted imaging (T1WI) and calculated using the formula: SIL1−4/SICSF. The Proton density fat fraction (PDFF) values and cross-sectional area (CSA) of paraspinal muscles were measured on the fat fraction map from the iterative decomposition of water and fat with the echo asymmetry and least-squares estimation quantitation (IDEAL-IQ) sequences and in/out phase images at the central level of the L4/5 and L5/S1 discs.

Results

This study included 476 patients with CLBP, including 189 males and 287 females. 69% had no Modic changes and 31% had Modic changes. There was no difference in CSA and PDFF for multifidus(MF) and erector spinae (ES) at both levels between Modic type I and type II, all P values>0.05. Spearman correlation analysis showed that VBQ was weakly negatively correlated with paraspinal muscles CSA (all r values < 0.3 and all p values < 0.05), moderately positive correlation with PDFF of MF at L4/5 level (r values = 0.304, p values<0.001) and weakly positively correlated with PDFF of other muscles (all r values<0.3 and all p values<0.001). Multivariate linear regression analysis showed that age (β = 0.141, p < 0.001), gender (β = 4.285, p < 0.001) and VBQ (β = 1.310, p = 0.001) were related to the total PDFF of muscles. For MCs, binary logistic regression showed that the odds ratio values of age, BMI and COD were 1.092, 1.082 and 1.004, respectively (all p values <  0.05).

Conclusions

PDFF of paraspinal muscles was not associated with Modic classification. In addition to age and gender, PDFF of paraspinal muscles is also affected by VBQ. Age and BMI are considered risk factors for the MCs in CLBP patients.

Peer Review reports

Introduction

Chronic low back pain(CLBP) refers to low back pain for more than 3 months [1]. It affects about 13% of adults, while also causing disability, high medical and social costs [2].

Low back pain can stem from various sources in the complex lumbar spine, including the vertebra, muscles, tendons, ligaments, fascia, facet joints, and discs [3]. Meanwhile, degenerative changes in the lumbar spine are more common in CLBP patients, such as Modic changes(MCs), cartilage endplate damage [4]. In recent years, there has been increasing attention paid to MCs as a potential source of low back pain [5]. The original classification of MCs was developed by Modic in 1988 [6]. MCs are classified into three types according to their T1-weighted and T2-weighted signal intensity(SI) on magnetic resonance imaging(MRI), Type I represents bone marrow edema and inflammation, type II represents the fatty degeneration of the bone marrow, and type III represents subchondral osteosclerosis [3, 6]. Most current studies focus on the relationship between MCs and low back pain, there are also conflicting views on the relationship between the two in different studies [7,8,9,10]. But degenerative changes in the structure of the spine are co-developed, including bone, endplate, intervertebral disc, paraspinal muscles and so on [11]. Studies have shown an association between changes in paraspinal muscles mass and CLBP-related spinal pathologies [12, 13]. However, current studies have rarely explored MCs relationship with the paraspinal muscles in CLBP patients.

In addition to the MCs of vertebral end-plate, the evaluation of the vertebral body is usually reflected by the bone mineral density(BMD) obtained by dual-energy X-ray absorptiometry (DEXA) scan or computed tomography (CT) scan [14]. A scoring system for vertebral bone quality (VBQ) based on MRI has been developed in the past few years [15] and it can protect patients from the effects of radiation [16], which has shown a correlate with bone mineral density (BMD) by using quantitative computed tomography (QCT) and has certain diagnostic efficiency [17]. Moreover, BMD has been found to be associated with paraspinal muscles fat infiltration [18].

In addition to the changes in the structure of the lumbar vertebral body, the patients with CLBP also have more obvious fat infiltration of the paraspinal muscles [12, 19]. The paraspinal muscles and the vertebral body play an important role in the stability of the lumbar spine. However, to the best of our knowledge, few studies have discussed the relationship between these two structural changes in patients with CLBP by MRI. Therefore, the purpose of this study was to investigate the relationship between paraspinal muscles fat infiltration and vertebral body related changes (VBQ score and MCs) in CLBP patients.

Materials and methods

Study participants

The study was approved by the institutional review board. All procedures were in accordance with the ethical standards of the 1964 Declaration of Helsinki. Verbal informed consent was obtained from all participants of the study. Participants in the study were enrolled at four multi-center healthcare facilities between July 2021 and December 2022. The patient’s gender, age, height, weight, and course of disease(COD) were collected. Body mass index(BMI) was calculated and reported in kilograms per meter square. Inclusion criteria were as follows: (1) Patients: people with CLBP (COD ≥ 3 months); (2) The age range: Non-child patients(≥ 14 years). Exclusion criteria are as follows: (1) Patients with contraindications to MR Examination and unable to cooperate with scanning; (2) Internal source of low back pain (such as urinary calculi); (3) Lumbar trauma, tumor, infection, surgery, etc.;(4) Genetic musculoskeletal diseases and neuromuscular diseases.(5) Pregnancy; (6) Athletes and regular bodybuilders who have professional muscle training; (7) Treatment(such as medications, physiotherapy) within 3 days before scan; (8) Diabetes and other chronic diseases.

Magnetic resonance scanning

The same 3.0T MRI scanners(MR750w, GE Healthcare, Waukesha, USA) were used at all four multi-center healthcare facilities involved in the study. MRI scan sequences included routine lumbar spine sequences (Sagittal-T1WI, Sagittal-T2WI, Transverse-T2WI and Sagittal- FS T2WI), quantitative sequence(Iterative decomposition of water and fat with the echo asymmetry and least-squares estimation quantitation, IDEAL-IQ) and IDEAL sequence. The protocol was standardized across the institutions. Specific scanning parameters are listed in Table 1. An abdominal pressure band with appropriate pressure was applied to the participant’s abdomen during the scan to reduce the effect of respiratory artifacts on image quality.

Table 1 Magnetic resonance scanning parameter

Image evaluation

MCs: According to the classification criteria [20], MCs are divided into 3 types, type I: hypointense on T1WI and hyperintense on T2WI, type II: hyperintense on T1WI and isointense or mildly hyperintense on T2WI, type III: hypointense on both T1WI and T2WI. The evaluation segments included L4-S1 levels.

Vertebral bone quality (VBQ) score: According to previous research method [21], The region of interest (ROI) was placed in the L1-4 vertebra and in the L3 level cerebrospinal fluid (CSF) in the median sagittal position(Fig. 1), and the signal intensity (SI) in each ROI was recorded. VBQ was calculated using the formula: SIL1−4/SICSF.

Fig. 1
figure 1

Image of a VBQ measurement. The quasi-circular ROI represents the L1-4 vertebral body, and the elliptical ROI represents the L3 level of cerebrospinal fluid; ROI = region of interest, VBQ = vertebral bone quality

Paraspinal muscles proton density fat fraction (PDFF) and cross-sectional area (CSA) measurement:

The central level of the L4-S1 discs were the level of interest and the PDFF values were measured on the fat fraction map from the IDEAL-IQ sequence. Outline along the boundary of the left and right multifidus (MF) and erector spinae (ES) respectively, and finally calculate the average value of the left and right muscles(Fig. 2). We also calculated the average PDFF of a total of 8 muscles in two levels. In addition, We mapped the CSA of MF and ES in the in/out phase images from the IDEAL sequence by drawing the same way as before (Fig. 3).

Fig. 2
figure 2

Image of a PDFF measurement. A, B represent the L4/5 and L5/S1 levels, respectively. Green indicates ES and red indicates MF. ES = erector spinae; MF = multifidus; PDFF = proton density fat fraction

Fig. 3
figure 3

Image of a CSA measurement. A, B represent the L4/5 and L5/S1 levels, respectively. Green indicates ES and red indicates MF. ES = erector spinae; MF = multifidus; CSA = cross-sectional area

Modic changes were evaluated by two experienced radiologists who were blinded to the patients’ medical records. When they do not agree, the specific Modic classification is evaluated by the superior physician(chief physician). Images of 45 patients were randomly selected for consistency test by the same personnel, who knew nothing about the previous measurement results.

Statistical analysis

SPSS 26.0 software (IBM Corp., Armonk, NY, USA) was used for statistical analysis. Kolmogorov-Smirnov test was used to verify whether the data fit the normal distribution. The data were expressed as mean ± standard deviation or median, and the intra-class consistency test was used for data reproducibility. The t-test of two independent samples was used for the data of normal distribution, and the Mann-Whitney U test was used for the data of non-normal distribution. Spearman correlation analysis was used for correlation analysis, the r values less than 0.3 is weak, 0.3–0.6 is moderate, and greater than 0.6 is strong [22]. Multifactor linear regression analysis was used to analyze the relationship between age, gender, BMI, VBQ, COD, MCs and the total paraspinal muscles PDFF. Binary logistic regression analysis was used to explore the risk factors for MCs. As the number of Modic III types was small in this study, it was not included in the statistical analysis. P value < 0.05 was considered statistically significant.

Result

Basic information about individuals

In this study, we included a total of 476 people, including 189 males and 287 females with a median age of 48 years, a median BMI of 22.49 kg/m2, and a median disease course of 24 months (Table 2). Among the population in this study, 69% had no MCs and 31% had MCs, There were 67 individuals(14%) with MCs at the L4/5 level, while there were 110 individuals(23%) with MCs at the L5/S1 level (Table 2).

Table 2 Baseline information of cohort

Reproducibility of measurements

The intraclass coefficient of CSA was 0.911, the intraclass coefficient of PDFF was 0.847, and the intraclass coefficient of VBQ was 0.941. All P values were <0.05.

Correlation analysis

Spearman correlation analysis showed that age was low negatively correlated with CSA of MF and ES at L5/S1 level(r=-0.099,-0.094,p<0.05, respectively). Age was positively correlated with PDFF(p<0.05) in all levels of muscles, with moderate correlation with L4/5 ES PDFF(r = 0.351, p<0.05) and weak correlation with the rest(r<0.3, p<0.05) (Table 3). There was a mild to moderate positive(r = 0.284–0.328,p<0.05) association between BMI and CSA in all levels of muscles, except for ES in L5/S1(p>0.05). However, for PDFF, except for a weakly positive correlation of ES at L4/5 (r = 0.130,p<0.05), the rest were not statistically significant, and there was a weakly negative correlation(r=-0.128, p<0.05) between BMI and VBQ (Table 3). The COD was negatively correlated with CSA of all levels of muscles, and positively correlated with PDFF of all levels of muscles, but there was no statistical significance(p>0.05), however, there was a slight positive correlation (r = 0.190, p<0.001) between the COD and VBQ (Table 3). VBQ was weakly negatively correlated with paraspinal muscles CSA(all r values<0.3 and all p values<0.05), moderately positive correlation with PDFF of MF at L4/5 level(r values = 0.304, p values<0.001) and weakly positively correlated with PDFF of other muscles(all r values<0.3 and all p values<0.001) (Table 3).

Table 3 Spearman correlation analysis

Comparison of Modic type I and type II

Due to the different degrees of correlation between clinical information and lumbar spine measurement data, we first compared clinical information in order to exclude the difference between the two groups, and there was no difference in clinical information between the two groups (all p values>0.05). There was no difference in CSA and PDFF for ME and ES at both levels between Modic type I and type II (all p values>0.05) (Table 4).

Table 4 Comparison of muscle measurement parameters at L4/5 and L5/S1 levels Modic I versus II

Multifactor linear regression analysis and binary logistic regression analysis

Multivariate linear regression analysis showed that age (β = 0.141,p<0.001), gender (β = 4.285,p<0.001) and VBQ(β = 1.310,p = 0.001) were related to the total PDFF of muscles (Table 5). For MCs, binary logistic regression analysis showed that the odds ratio values of age, BMI and COD were 1.092, 1.082 and 1.004, respectively (all p values  < 0.05) (Table 6).

Table 5 Multifactor linear regression analysis of total muscles PDFF
Table 6 Binary logistic regression of risk of Modic change in CLBP petients

Discussion

In this study, we tried to explore the relationship between VBQ and MCs and paraspinal muscles in patients with CLBP. In addition, we also conducted a correlation analysis between clinical information and paraspinal muscles related parameters. We found that VBQ is associated with changes in the paraspinal muscles.

As for MCs, the current researches focused on their relationship with low back pain [7, 23]. In our study, we tried to explore the relationship between it and paraspinal muscles mass, but we found that although the PDFF of MF and ES at both levels was higher in patients with MCs type I compared with MCs type II, there was no statistical difference between them in this study. In addition, we also compared the CSA of MF and ES, and there was no statistical difference between the two types. There is currently little literature exploring the relationship between MCs and paraspinal muscles, and one study found that LBP patients with Modic type I and I/II changes had significantly higher percentage fat content compared with patients without MCs [24]. In addition, some scholars have found that at the L4/L5 and L5/S1 intervertebral disc levels, the CSA of psoas in patients with MCs is smaller than that in patients without these changes [25]. However, those studies did not explore the differences between Modic type I and II. The MCs of type I and type II are the most common patterns observed in the lumbar spine [7]. In this study, the incidence of MCs was 31%, and the amount of Modic type III was small, so it was not included in the statistical analysis. In addition, Other scholars have found that age is associated with MCs in the lower lumbar region(L4-S1) [26], and the results of this study are consistent with it, and we found that age is a risk factor(OR = 1.092). Other scholars have also pointed out that MCs are related to weight-related factors (BMI, waist circumference) [27], and this study also has similar finding. Binary logistic regression analysis showed that the OR of age, BMI and COD were all greater than 1, and the P-values were all less than 0.05, but the OR values were too small to have clinical value, which may be caused by numerous factors affecting MCs.

VBQ is based on MRI, which has been used to measure fatty infiltration within the vertebral body [15, 28]. In this study, it was found by multiple regression that there was a relationship between VBQ and total paraspinal muscles PDFF, and the PDFF increased by 1.31 units for every unit increase of VBQ. Through correlation analysis, we also found that VBQ was weakly negatively correlated with paraspinal muscles CSA, and was weakly positively correlated with PDFF of other muscles except for moderately positive correlation with PDFF of MF at L4/5 level. To our knowledge, there were few studies to explore the association of VBQ with paraspinal muscles fat infiltration and the specific causal relationship requires further longitudinal study. However, several studies have found that lower BMD may be associated with more severe paraspinal muscles fatty infiltration [18, 29, 30]and these studies may provide some reference value to our result. We also found a weak positive correlation between VBQ and age, which may indicate that the vertebral body has more fat with increased age. Interestingly, we observed a weak negative correlation between BMI and VBQ. Generally speaking, BMI and VBQ are both parameters that reflect fat to some extent, and there should be a positive linear relationship between the two. However, BMI reflects the general fat content and can’t accurately reflect the fat situation of the tissue [31], while VBQ mainly reflects the fat content of the vertebral body [28]. This may explain the relationship between the two in this study.

It has been found that the CSA of paraspinal muscles decreases with age [32, 33], this is similar to the results in our study that we found a very weak negative correlation between age and CSA of MF and ES at the L5/S1 level. In addition to the influence of age on paraspinal muscles CSA, the present studies also confirmed that age also affects paraspinal muscles fat infiltration [32, 34], which is similar to our results. We also found a low to moderate positive correlation between BMI and most paraspinal muscles’ CSA,

and only a weak positive correlation with PDFF of ES at the L4/5 level, a study of paraspinal and lower limb muscles also found that BMI was not associated with fat infiltration in paraspinal muscles [35], and it has been pointed out that BMI does not estimate the amount of fat in different parts of the body [31], which may be the reason why there was less correlation between BMI and paraspinal muscles PDFF in this study. We did not find a relationship between the COD and paraspinal PDFF, which may be due to the large span of the course of the disease in this study and the lack of detailed records of patients’ low back pain, whether it was persistent or intermittent. Finally, through multi-factor regression, we found that age and gender also had an impact on the PDFF of paraspinal muscles, which was similar to the results of Huang Y and Huang R et al [12, 32]. Of course, in order to exclude the influence of potential factors on paraspinal muscles’ PDFF, we also excluded patients with chronic diseases such as diabetes. Because diabetes is a chronic metabolic disease, which has been reported individuals with impaired glucose metabolism showed significantly higher PDFF of muscles compared with controls [36, 37].

There are some limitations to our study. First of all, because the amount of Modic III was very small in this study, statistical analysis was not included. Secondly, based on the mismatch of clinical information, we did not compare the relevant parameters of the population with and without MCs. Finally, in this study, we selected the L4-S1 levels as the level of interest, which could not represent the entire lumbar spine levels, but studies have shown that MCs and paraspinal muscles fat infiltration are more obvious at the lower lumbar spine level [26, 38].

Conclusion

There is no changes in paraspinal muscles’ CSA and PDFF between Modic type I and type II. VBQ has a positive linear relationship with total paraspinal muscles’ PDFF. In addition to age and gender, total paraspinal muscles’ PDFF is also affected by VBQ. Age has a weak negative correlation with paraspinal muscles CSA, a low to moderate positive correlation with PDFF, and a low positive correlation with VBQ. Overall, BMI had a low to moderate positive correlation with paraspinal muscles CSA, but no correlation with PDFF. Age and BMI are considered risk factors for the MCs in CLBP patients.

Data availability

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

BMD:

Bone mineral density

BMI:

Body mass index

CT:

Computed tomography

CLBP:

Chronic low back pain

COD:

Course of disease

CSA:

Cross-sectional area

CSF:

Cerebrospinal fluid

DXA:

Dual-energy X-ray absorptiometry

ES:

Erector spinae

MF:

Multifidus

IDEAL-IQ:

Iterative decomposition of water and fat with the echo asymmetry and least-squares estimation quantitation

MCs:

Modic changes

MRI:

Magnetic resonance imaging

PDFF:

Proton density fat fraction

QCT:

Quantitative computed tomography

ROI:

Region of interest

SI:

Signal intensity

VBQ:

Vertebral bone quality

References

  1. Deyo RA, Dworkin SF, Amtmann D, et al. Focus article: report of the NIH Task Force on Research Standards for Chronic Low Back Pain. Eur Spine J. 2014;23:2028–45.

    Article  PubMed  Google Scholar 

  2. Stevans JM, Delitto A, Khoja SS, et al. Risk factors Associated with Transition from Acute to Chronic Low Back Pain in US patients seeking primary care. JAMA Netw Open. 2021;4:e2037371.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Knezevic NN, Candido KD, Vlaeyen JWS, Van Zundert J, Cohen SP. Low back pain Lancet. 2021;398:78–92.

    PubMed  Google Scholar 

  4. Bailey JF, Fields AJ, Ballatori A, et al. The relationship between Endplate Pathology and patient-reported symptoms for chronic low back Pain depends on lumbar Paraspinal. Muscle Qual Spine. 2019;44:1010–7.

    Google Scholar 

  5. Hebelka H, Brisby H, Erkmar A, Lagerstrand K. Loading of the spine in low back Pain patients does not induce MRI changes in Modic lesions: a prospective clinical Study Diagnostics. (Basel). 2022;12:1815.

    Google Scholar 

  6. 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.

  7. Teraguchi M, Hashizume H, Oka H, et al. Detailed subphenotyping of Lumbar Modic Changes and Their Association with Low Back Pain in a large Population-based study: the Wakayama Spine Study. Pain Ther. 2022;11:57–71.

    Article  PubMed  Google Scholar 

  8. Annen M, Peterson C, Leemann S, Schmid C, Anklin B, Humphreys BK. Comparison of outcomes in MRI confirmed lumbar disc herniation patients with and without modic changes treated with high velocity, low amplitude spinal manipulation. J Manip Physiol Ther. 2016;39:200–9.

    Article  Google Scholar 

  9. Udby PM, Bendix T, Ohrt-Nissen S et al. Modic changes are not Associated with Long-Term Pain and Disability Spine 2019; 44: 1186–92.

  10. Saukkonen J, Määttä J, Oura P, et al. Association between Modic Changes and Low Back Pain in Middle Age. Spine. 2020;45:1360–7.

    Article  PubMed  Google Scholar 

  11. Kauppila LI. Atherosclerosis and Disc Degeneration/Low-Back Pain – a systematic review. Eur J Vasc Endovasc Surg. 2009;37:661–70.

    Article  CAS  PubMed  Google Scholar 

  12. Huang Y, Wang L, Zeng X, et al. Association of Paraspinal Muscle CSA and PDFF measurements with lumbar intervertebral disk degeneration in patients with chronic low back Pain Front Endocrinol. (Lausanne). 2022;13:792819.

    Article  Google Scholar 

  13. Cooley JR, Walker BF, Ardakani M, Kjaer E, Jensen P, Hebert TS. Relationships between paraspinal muscle morphology and neurocompressive conditions of the lumbar spine: a systematic review with meta-analysis. BMC Musculoskelet Disord. 2018;19:351.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Soliman MAR, Aguirre AO, Kuo CC, et al. Vertebral bone quality score independently predicts cage subsidence following transforaminal lumbar interbody fusion. Spine J. 2022;22:2017–23.

    Article  PubMed  Google Scholar 

  15. Chen Z, Lei F, Ye F, Yuan H, Li S, Feng D. MRI-based vertebral bone quality score for the assessment of osteoporosis in patients undergoing surgery for lumbar degenerative diseases. J Orthop Surg Res. 2023;18:257.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Ehresman J, Pennington Z, Schilling A, et al. Novel MRI-based score for assessment of bone density in operative spine patients. Spine J. 2020;20:556–62.

    Article  PubMed  Google Scholar 

  17. Salzmann SN, Okano I, Jones C, et al. Preoperative MRI-based vertebral bone quality (VBQ) score assessment in patients undergoing lumbar spinal fusion. Spine J. 2022;22:1301–8.

    Article  PubMed  Google Scholar 

  18. Zhao Y, Huang M, Serrano Sosa M, et al. Fatty infiltration of paraspinal muscles is associated with bone mineral density of the lumbar spine. Arch Osteoporos. 2019;14:99.

    Article  PubMed  Google Scholar 

  19. Mengiardi B, Schmid MR, Boos N et al. Fat content of lumbar paraspinal muscles in patients with chronic low back pain and in asymptomatic volunteers: quantification with MR spectroscopy Radiology 2006; 240: 786–92.

  20. Chen Y, Bao J, Yan Q, Wu C, Yang H, Zou J. Distribution of Modic changes in patients with low back pain and its related factors. Eur J Med Res. 2019;24:34.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Li W, Zhu H, Hua Z et al. Vertebral bone quality score as a predictor of pedicle screw loosening following surgery for degenerative lumbar Disease Spine (Phila Pa 1976) 2023; 48(23): 1635–41.

  22. Wang Q, Jia M, Jiang M, et al. Seesaw Effect between COVID-19 and Influenza from 2020 to 2023 in World Health Organization regions. Correlation Anal JMIR Public Health Surveillance. 2023;9:e44970.

    Article  Google Scholar 

  23. Bailly F, Maigne JY, Genevay S, et al. Inflammatory pain pattern and pain with lumbar extension associated with Modic 1 changes on MRI: a prospective case-control study of 120 patients. Eur Spine J. 2014;23:493–7.

    Article  CAS  PubMed  Google Scholar 

  24. Atci IB, Yilmaz H, Samanci MY, Atci AG, Karagoz Y. The prevalence of Lumbar Paraspinal Muscle Fatty Degeneration in patients with Modic type I and I/II end plate changes Asian. Spine J. 2020;14:185–91.

    Google Scholar 

  25. Arbanas J, Pavlovic I, Marijancic V, et al. MRI features of the psoas major muscle in patients with low back pain. Eur Spine J. 2013;22:1965–71.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 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.

    Article  PubMed  Google Scholar 

  27. Kuisma M, Karppinen J, Haapea M, et al. Are the determinants of vertebral endplate changes and severe disc degeneration in the lumbar spine the same? A magnetic resonance imaging study in middle-aged male workers. BMC Musculoskelet Disord. 2008;9:51.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Ehresman J, Schilling A, Yang X, et al. Vertebral bone quality score predicts fragility fractures independently of bone mineral density. Spine J. 2021;21:20–7.

    Article  PubMed  Google Scholar 

  29. Li X, Zhang Y, Xie Y, Lu R, Tao H, Chen S. Correlation between bone Mineral density (BMD) and paraspinal muscle Fat Infiltration based on QCT: a cross-sectional study. Calcif Tissue Int. 2022;110:666–73.

    Article  CAS  PubMed  Google Scholar 

  30. Yang Q, Yan D, Wang L, et al. Muscle fat infiltration but not muscle cross-sectional area is independently associated with bone mineral density at the lumbar spine. Br J Radiol. 2022;95:20210371.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Brady SRE, Urquhart DM, Hussain SM, et al. High baseline fat mass, but not lean tissue mass, is associated with high intensity low back pain and disability in community-based adults. Arthritis Res Ther. 2019;21:165.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Huang R, Pan F, Kong C, Lu S. Age- and sex-dependent differences in the morphology and composition of paraspinal muscles between subjects with and without lumbar degenerative diseases. BMC Musculoskelet Disord. 2022;23:734.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Murata Y, Nakamura E, Tsukamoto M, et al. Longitudinal study of risk factors for decreased cross-sectional area of psoas major and paraspinal muscle in 1849 individuals. Sci Rep. 2021;11:16986.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Peng X, Li X, Xu Z, et al. Age-related fatty infiltration of lumbar paraspinal muscles: a normative reference database study in 516. Chin Females Quant Imaging Med Surg. 2020;10:1590–601.

    Article  Google Scholar 

  35. Dahlqvist JR, Vissing CR, Hedermann G, Thomsen C, Vissing J. Fat replacement of Paraspinal Muscles with aging in healthy adults. Med Sci Sports Exerc. 2017;49:595–601.

    Article  PubMed  Google Scholar 

  36. Kiefer LS, Fabian J, Rospleszcz S, et al. Assessment of the degree of abdominal myosteatosis by magnetic resonance imaging in subjects with diabetes, prediabetes and healthy controls from the general population European. J Radiol. 2018;105:261–8.

    Google Scholar 

  37. Diallo TD, Rospleszcz S, Fabian J, et al. Associations of myosteatosis with disc degeneration: a 3T magnetic resonance imaging study in individuals with impaired glycaemia. J Cachexia Sarcopenia Muscle. 2023;14:1249–58.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Sollmann N, Bonnheim NB, Joseph GB, et al. Paraspinal muscle in Chronic Low Back Pain: comparison between Standard parameters and Chemical Shift Encoding-Based Water–Fat MRI. J Magn Reson Imaging. 2022;56:1600–8.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work is supported by the National Natural Science Foundation of China (No. 82260338, 82302131), and the Yunnan Fundamental Research Projects (No. 202201AC070669, 202201AU070051, 202301AS070016, 202001AY070001-200).

Author information

Authors and Affiliations

Authors

Contributions

ZZ, BH, JY and YH have made important contributions to the design, guidance and quality control of this study. ZW, DZ, ML, FP, JC and YY made important contributions to patient collection, data measurement and literature collection for this study. JC, YY, YH made important contributions to the manuscript writing and data statistics of this study. All authors approved the final version of the manuscript, contributed to the article, and approved the submitted version.

Corresponding authors

Correspondence to Juntao Yang, Zhenguang Zhang or Bo He.

Ethics declarations

Ethics approval and consent to participate

The studies approved by Research Ethics Committee of the First Affiliated Hospital of Kunming Medical University. The participants provided their verbal informed consent to participate in this study; No animal studies are presented in this manuscript; No potentially identifable human images or data is presented in this study. We solemnly promise that all the methodology of this study has been carried out in accordance with the relevant guidelines and provisions of the Declaration of Helsinki, among others.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, J., Huang, Y., Yang, Y. et al. Vertebral bone quality score was associated with paraspinal muscles fat infiltration, but not modic classification in patients with chronic low back pain: a prospective cross-sectional study. BMC Musculoskelet Disord 25, 509 (2024). https://doi.org/10.1186/s12891-024-07626-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-024-07626-4

Keywords