Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 25 August 2025

Sec. Bone Research

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1591044

The q-Dixon sequence for MRI predicts osteoporotic vertebral compression fractures

Jing Zhang&#x;Jing ZhangQiyuan Li&#x;Qiyuan LiYao WangYao WangLi SunLi SunQingyuan ZhangQingyuan ZhangChuanping Gao*Chuanping Gao*
  • Department of Radiology, The Affiliated Hospital of Qingdao University, Qingdao, Shandong, China

Objectives: To evaluate whether q-Dixon sequence-based fat fraction (FF) values of the lumbar spine can predict osteoporotic vertebral compression fracture (OVCF) risk in older adult(s) osteoporosis patients.

Materials & methods: Thirty OVCF patients and 15 osteoporosis patients were enrolled. Areas of interest (ROIs) were manually drawn using the post-processing workstation, and FF values of the patient’s L1–L4 vertebrae (except the fractured vertebrae) were measured. The Pearson correlation test was used to analyze the correlation between the average lumbar spine FF value and bone mineral density (BMD). The receiver operating characteristic curve (ROC) and area under the curve (AUC) were used to evaluate the prediction efficiency of the FF values and to determine the best cut-off value for prediction.

Results: The average lumbar spine FF value and the FF values of the L1–L4 vertebrae in the fracture group were significantly higher than those in the non-fracture group, and there was no significant difference in BMD between the two groups. ROC analysis showed that the AUC of the average lumbar spine FF value was 0.822, with sensitivity = 73.3%, specificity = 86.7%, and cut-off value = 57.27%. Among the L1–L4 vertebrae, the FF value of L2 vertebrae had the highest AUC of 0.870, and the cutoff value was 56.62%.

Conclusion: The FF values of the lumbar spine measured by the q-Dixon sequence can help predict OVCF risk and provide complementary information to BMD measurements. The FF value of L2 vertebrae has the best prediction efficiency. However, these findings should be interpreted with caution given the relatively small sample size (n=45) and manual ROI segmentation method used in this study

Introduction

Osteoporotic vertebral compression fracture (OVCF) is the most common osteoporotic fracture in the clinic and can occur in response to slight external force under conditions of decreased bone mass and bone strength and increased bone fragility (1). According to the International Osteoporosis Foundation, approximately 50% of women and 25% of men worldwide are expected to experience at least one osteoporotic fracture during the latter half of their lives (2). As the population ages, OVCFs are increasingly affecting the older adult(s), decreasing patients’ ability to perform many routine activities of daily living and even causing patient death. Therefore, it is important to identify individuals at high risk of OVCF among osteoporosis patients, to predict fracture risk, and to initiate early intervention.

Dual energy X-ray absorptiometry (DXA) is widely used in the diagnosis and evaluation of osteoporosis (3). However, BMD is representative of only 50%–70% of bone strength, and it exposes patients to radiation (4). Therefore, alternative radiation-free imaging techniques and biomarkers need to be investigated. Despite these limitations, DXA remains the clinical gold standard due to its low cost, widespread availability, and rapid acquisition time. Studies have shown that an increase in fat cells in the bone marrow is accompanied by a decrease in osteoblasts and that fat tissue is used to “fill” the extra bone marrow space when bone mass is reduced (5). Thus, the onset of osteoporosis is associated with an increase in bone marrow fat mass.

Chemical shift-encoding water-fat imaging (Dixon) is an effective method for rapid and non-invasive assessment of the bone marrow fat content and was first proposed by Dixon in 1984 (6). Initially, the two-point Dixon sequence was proposed, wherein the hydrogen proton signals in water and fat are collected when they are in the in-phase (IP) and opposed-phase (OP), respectively, by setting different echo times and calculating the water–fat separation images. The main limitation of this method is its sensitivity to magnetic field inhomogeneity. To address this problem, Glover and Schneider proposed the three-point Dixon sequence (7). Along with the advances in MR technology, various fat content quantification sequences based on the three-point Dixon method have been developed, including Ideal IQ (GE Healthcare), mDixon (Philips Healthcare), and q-Dixon (Siemens Healthcare). In recent years, the modified Dixon sequence has been increasingly applied in the quantitative analysis of fat in vertebrae.

We hypothesize that: (1) FF values will be significantly higher in patients with OVCF compared to those without fractures, and (2) FF measurements will provide fracture risk prediction independent of BMD. In this study, the q-Dixon sequence was used to measure the fat fraction (FF) value of the lumbar spine in osteoporosis patients. We hypothesized that FF values may predict OVCF to a certain extent.

Materials and methods

Patients

This prospective study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by the institutional Review Board of the Affiliated Hospital of Qingdao University, and informed consent was obtained from all patients.

A total of 30 older adult(s) patients presenting with OVCF (fracture group) and 15 older adult(s) patients diagnosed with osteoporosis during the same period (non-fracture group) were enrolled from July 2022 to December 2023. The sample size of 45 patients (30 fracture, 15 control) was determined based on clinical availability rather than formal power calculation, which may limit the generalizability of our findings. Patients were consecutively enrolled without randomization to reflect real-world clinical practice. All participants provided written informed consent within 24 hours before undergoing MRI and BMD assessments. All patients were examined within a week of diagnosis, and MRI and BMD examinations were performed on the same day.

The inclusion criteria for the fracture group were as follows: (1) no history of fragility fractures; (2) osteoporosis confirmed by BMD measurement; and (3) fresh vertebral fracture identified on MRI. Fresh fracture was defined as acute vertebral collapse with MRI evidence of bone marrow edema (T1 hypointensity, T2/STIR hyperintensity) and symptom onset within ≤4 weeks. The exclusion criteria were as follows: (1) poor image quality; (2) number of fractured lumbar vertebrae ≥ 2; (3) compression fracture caused by high-energy trauma; (4) anti-osteoporosis therapy prior to MRI and BMD measurement; (5) presenting with other diseases that have the potential to impact bone quality; (6) and taking drugs that could affect bone metabolism.

The inclusion criteria for the non-fracture group were as follows: (1) osteoporosis was confirmed by BMD measurement; and (2) imaging studies did not show any signs of fragility fractures. The exclusion criteria were as follows: (1) poor image quality; (2) anti-osteoporosis therapy prior to MRI and BMD measurement; (3) presenting with other diseases that have the potential to impact bone quality; and (4) taking drugs that could affect bone metabolism.

Figure 1 illustrates the patient enrollment process. From an initial screening cohort of 200 eligible candidates, 45 participants met the final inclusion criteria, comprising 30 patients with acute osteoporotic vertebral compression fractures (OVCF) and 15 age-matched osteoporosis controls without vertebral fractures.

Figure 1
Flowchart detailing a study process: 200 patients were screened. After informed consent, they were divided into potential fracture and non-fracture groups, each with 100 patients. Exclusions led to the final study cohort of 30 fracture and 15 non-fracture patients. MRI and bone mineral density tests were performed, within seven days of diagnosis.

Figure 1. Flowchart of patient selection and study enrollment.

MRI examination

MRI was performed using a Prisma 3.0T MRI system (Siemens Healthcare GmbH, Erlangen, Germany) with an 18-channel body phased-array coil. The examination site was the lumbar spine. Anatomical and morphological evaluation of the vertebrae was obtained using conventional T1WI (sagittal, repetition time [TR]/echo time [TE] = 400/20 ms, field of view [FOV] = 300 × 300 mm2, matrix size = 240 × 320, slice thickness = 4.0 mm, gap = 1.0 mm, and acquisition time = 124 s) and T2WI (sagittal, TR/TE = 2700/119 ms, FOV = 300 × 300 mm2, matrix size = 288 × 384, slice thickness = 4.0 mm, gap =1.0 mm, and acquisition time = 118 s). The fat mass of the vertebrae was obtained using the q-Dixon sequence (sagittal, TR/TE = 9.40/1.29 ms, FOV = 380 × 380 mm2, matrix size = 96 × 96, slice thickness = 3.0 mm, gap = 1.0 mm, and acquisition time =108 s).

Imaging analysis

Post-processing was performed on a Syngo.via workstation. The sagittal FF map was automatically generated and transmitted to the workstation after the q-Dixon sequence. Rectangular regions of interest (ROIs) were manually drawn on the FF map by two radiologists experienced in musculoskeletal radiology (Figure 2). Manual ROI placement was performed by two experienced radiologists to accommodate anatomical variations, though this approach may introduce more subjectivity compared to automated segmentation methods. During the segmentation, cortical bone and the vertebral endplate were avoided as much as possible, and vertebrae with fractures or degenerative changes (e.g., Modic changes) were excluded. At least three vertebrae were measured for each patient. Each vertebra was measured three times and averaged. The average of the two radiologists’ measurements was taken as the final result. After measurement, the FF values of the L1–L4 vertebrae of each patient were recorded, and the average FF values of the lumbar spine were calculated.

Figure 2
Four MRI scans labeled a, b, c, and d display transverse sections of the lumbar spine. Image a shows a regular view, b is similar but appears slightly enhanced, c includes prominent spinal fluid contrast, and d highlights specific vertebral regions with outlined areas.

Figure 2. Sagittal T1-weighted (a), T2-weighted (b), Fat-sat T2-weighted (c) images and fat fraction (FF) map (d) of a 73-year-old female patient with an osteoporotic vertebral compression fracture in L1 vertebrae.

BMD measurement

Bone mineral density (BMD) was quantified using dual-energy X-ray absorptiometry (DXA; Osteosys Primus, Seoul, Korea) following standardized protocols. Scans were performed by certified technologists at two anatomical sites: (1) lumbar spine (L1-L4 vertebrae) and (2) right proximal femur. The system was calibrated daily using manufacturer-provided phantoms, with a demonstrated coefficient of variation <1% for lumbar spine measurements.For each participant:(1)Areal BMD (g/cm²) was automatically calculated for all measurable vertebrae.(2)Fractured vertebrae were systematically excluded from lumbar spine analysis.(3)T-scores were derived using the results of BMD studies in normal individuals.(4)Scan quality was verified by a musculoskeletal radiologist blinded to clinical data.

Statistical analysis

SPSS 26.0 and MedCalc 22.016 software were used for statistical calculations with a two-sided level of significance of 0.05. The Shapiro–Wilk test was used to assess the normality of the data. Student’s t-test or chi-square test was used to analyze the data with a normal distribution. The Mann–Whitney U test was used to analyze non-normally distributed data. The intraclass correlation coefficient (ICC) was used to evaluate the consistency of the FF values measured by the two radiologists. An ICC > 0.80 was considered good consistency. The Pearson correlation coefficient was used to evaluate the correlation between the FF value and BMD. The receiver operating characteristic curve (ROC) was plotted, and the area under the curve (AUC), sensitivity, and specificity were calculated to assess the predictive value of the FF value for OVCF. The cut-off value of the FF value was determined using the Youden index, which was calculated as: sensitivity + specificity − 1.

Results

Table 1 shows the general characteristics of patients in the fracture and non-fracture groups. The age range of the fracture group was 53–87 years (66.67 ± 7.30), and the age range of the non-fracture group was 52–87 years (68.33 ± 8.96). No significant differences in age, sex, height, weight, and BMI were found between the two groups (P > 0.05). In the fracture group, there were 9 fractures in the L1 vertebrae, 10 in the L2 vertebrae, 4 in the L4 vertebrae, and 7 in the thoracic vertebrae, including 1 in the T8 vertebrae, 2 in the T11 vertebrae, and 3 in the T12 vertebrae.

Table 1
www.frontiersin.org

Table 1. Patients characteristics.

The consistency in the FF values of the L1–L4 vertebrae measured by the two radiologists was excellent, with ICC values of 0.996 (95% confidence interval, 0.995–0.997). Figure 3 shows the results of the Pearson correlation analysis between the average FF value of the lumbar spine and BMD in the two groups. A significant negative correlation was observed between the average lumbar spine FF value and BMD (r = −0.6489, P < 0.001).

Figure 3
Scatter plot showing the relationship between BMD (g/cm²) on the y-axis and FF (%) on the x-axis. Data points exhibit a downward trend, with a fitted line indicating a negative correlation. Dotted lines represent confidence intervals.

Figure 3. Scatter plot of the correlation between the average FF value of the lumbar spine and BMD.

Table 2. Comparison of vertebral fat fraction (FF) values between fracture and non-fracture groups and Figure 4 show that the average lumbar spine FF value and the FF values of the L1–L4 vertebrae in the fracture group were significantly higher than those in the non-fracture group. The FF values in the two groups increased from the L1 to L4 vertebrae, and the difference was statistically significant (P < 0.05). There was no significant difference in BMD between the two groups (P > 0.05). The ROC curves of the average FF value of the lumbar spine and the FF values of the L1–L4 vertebrae are shown in Figure 5; Table 3 shows the AUC, sensitivity, specificity, cut-off value, and other parameters of the lumbar spine FF value. The FF value of the L2 vertebrae had the highest AUC of 0.870, with a sensitivity and specificity of 80%, and its cut-off value was 56.62%. The superior predictive performance at L2 (AUC=0.870) may reflect its transitional biomechanical position in the spine, experiencing both compressive and shear forces that make it particularly sensitive to marrow composition changes. The FF value of the L3 vertebrae had the highest specificity of 100%, but the AUC was only 0.752, and the cut-off value was 61.27%. The AUC of the average FF value of the lumbar spine was 0.822, with a sensitivity of 73.3%, specificity of 86.7%, and cut-off value of 57.27%.

Table 2
www.frontiersin.org

Table 2. Differences in BMD and FF values between fracture group and non-fracture group .

Figure 4
Box plot charts comparing various parameters between OVCF and HC groups. Metrics include BMD, FF, and vertebral levels L1 to L4 FF percentages. Each plot shows data distribution and interquartile ranges, with OVCF generally having wider ranges than HC.

Figure 4. Analysis of the differences between the fracture group and the non-fracture group. OVCF: fracture group; HC: non-fracture group.

Figure 5
ROC curve graph showing sensitivity versus 100% specificity for five lines: L1-FF (AUC=0.733), L2-FF (AUC=0.870), L3-FF (AUC=0.752), L4-FF (AUC=0.831), and FF (AUC=0.822). The diagonal is a reference line.

Figure 5. ROC analysis of the average FF value of the lumbar spine and the FF values of the L1–L4 vertebrae.

Table 3
www.frontiersin.org

Table 3. Performance of the parameters for predicting OVCF .

Discussion

In this study, we evaluated the ability of the q-Dixon sequence to measure the lumbar fat mass in osteoporosis patients and explored its feasibility for predicting OVCF. The q-Dixon sequence is based on chemical shifts, which exploit the differential separation of proton resonance frequencies in water and fat to achieve high-resolution and high-contrast images of tissue structure (6). This technology can collect six echoes at one time and can accurately identify the water and fat signals in human tissue using the 7-peak fat model and T2*correction, with high spatial resolution and convenient operation (8). The findings indicate that the q-Dixon sequence is effective in predicting OVCF risk, with the highest predictive efficacy in the L2 vertebral FF measurement. While we interpret elevated FF values as a predisposing factor for fracture, we acknowledge that bone marrow fat is dynamic. Although FF changes typically develop over weeks to months, post-fracture immobility or inflammatory responses could theoretically influence FF measurements even within the first week after fracture.This aligns with broader research that emphasizes the importance of accurate, early OVCF detection to prevent adverse outcomes and manage the healthcare burden.

A meta-analysis showed that BMD reduction could explain approximately 70% of the risk of osteoporotic fractures (9). For postmenopausal osteoporosis patients, a 1-standard deviation decrease in BMD was associated with a 1.5–2.0-fold increase in fracture risk (10). However, BMD has some limitations in fracture risk assessment. Most patients with osteoporotic fractures have a BMD that is not within the range of diagnosis of osteoporosis (1113). Therefore, factors related to bone mass should be considered in the assessment of patients with osteoporotic fractures. Bone marrow fat is considered an important biological marker for bone quality assessment (14, 15). As bone grows and ages, the number of adipocytes in bone marrow gradually increases, which is closely related to bone metabolism. The only difference is that bone mass continues to increase during growth and sexual maturity, whereas as the body ages, sex hormone activity decreases, resulting in bone mass loss (16). Adipocytes and osteoblasts in the bone marrow originate from the same precursor, namely, marrow mesenchymal stem cells (MSCs). MSCs have multilineage differentiation potential. Competition exists between adipocytes and osteoblasts, which are in a dynamic equilibrium state. A decrease in osteoblasts tends to be accompanied by an increase in adipocytes (17, 18). Moreover, adipocytes secrete saturated fatty acids. Saturated fatty acids have been suggested to possibly have a lipotoxic effect on osteoblasts, thereby affecting their normal function (19). Previous studies have shown that changes in the bone marrow fat precede changes in bone mass during the development of osteoporosis (20, 21).

Our findings demonstrate that FF provides complementary information to BMD, with the inverse correlation (r=-0.6489) suggesting FF captures distinct aspects of bone quality. This supports emerging paradigms of osteoporosis as both a quantitative (BMD) and qualitative (marrow composition) disorder of bone (17, 21, 22). The average lumbar spine FF value gradually increased with a decrease in BMD, reflecting the competition between adipogenesis and osteogenesis in bone marrow, in agreement with the conclusions of previous studies (22, 23). The increase in the fat content may indirectly reflect a decrease in bone mass and change in bone strength. The reasons for these results may include the following (2427): (1) the balance between adipocyte and osteoblast differentiation is disrupted, which affects bone strength; (2) saturated fatty acids secreted by adipocytes have lipotoxic effects on osteoblasts; (3) fat accumulation can replace functional hematopoietic cells and osteoblasts because of the restricted space within the bone marrow; and (4) with an increase in age, the levels of inhibitory factors such as estrogen and transforming growth factor decrease, which is not conducive to osteoblast differentiation.

The average FF value of the lumbar spine and the FF values of the L1–L4 vertebrae in patients with fracture were higher than those in patients without fractures, and there was no significant difference in BMD between the two groups. This finding is similar to that of Gassert et al. (28), who evaluated the ability of the proton density fat fraction (PDFF) of vertebral bone marrow to differentiate between osteoporotic/osteopenic patients with and without vertebral compression fractures. The results showed that patients with fractures had a significantly higher PDFF than those without fractures after adjusting for clinical factors (P < 0.001), although BMD showed no significant differences among the subgroups. This effect is not just simply caused by the strong correlation between the FF values and BMD. The FF values should be considered as an important component in the assessment of bone fragility to assess fracture risk in osteoporosis patients. Furthermore, consistent with previous studies (29, 30), the vertebral FF values increased from L1 to L4, reflecting the transformation of the bone marrow from red to yellow with age, and from peripheral to axial, but the pathophysiological relevance of this finding is unclear.

The ROC curve showed that the FF value of the L2 vertebrae had the best performance in predicting OVCF (AUC = 0.870, sensitivity = 80%, and specificity = 80%), the FF value of the L3 vertebrae had the highest specificity of 100% (AUC = 0.752). This variation in diagnostic performance across vertebral levels may reflect anatomical differences in biomechanical loading patterns, regional variations in marrow composition, or statistical variability due to sample size limitations. The AUC of the average FF value of the lumbar spine was 0.822 (sensitivity = 73.3%, specificity = 86.7%), and the cutoff value was 57.27%. The combination of multiple vertebrae was likely to result in reduced predictive sensitivity. Overall, the average FF value of the lumbar spine showed stable ability in predicting OVCF. Therefore, FF values can be considered as a complement to BMD assessment of bone strength in the risk assessment of vertebral fractures in osteoporosis patients. Yun et al. conducted a similar study and reported a lower cut-off value of 42% (31); Gassert et al. calculated a cut-off value of 44.9% (28). Potential reasons for why the cutoff value in this study is higher than those reported in similar studies may include the following. (1) Individual differences exist in the samples, such as the inclusion criteria and age range of the samples. Unlike most previous studies, which were retrospective and included a wide age range of patients, this study was a prospective study of older adult(s) patients with osteoporosis who had high levels of bone marrow fat. (2) Group differences such as sex, region, and race need to be confirmed by expanding the sample size and performing further epidemiological investigation. (3) Differences in technical metrics, such as the MR scanners used in different studies and the technical principles of the MR scanners, may lead to differences in the results.

Accurate prediction of osteoporotic fracture risk can be of significant clinical benefit in the assessment and management of osteoporosis. Beyond its predictive value, FF measurement offers unique advantages as it reflects aspects of bone quality not captured by BMD. Importantly, FF is modifiable - studies demonstrate that interventions including bisphosphonates, teriparatide, exercise, and hormonal therapy can reduce marrow fat content (17, 21, 22). This suggests FF may serve as both a predictive biomarker and a treatment response indicator. The FF value is readily available clinically and has good predictive efficacy, which can help clinicians distinguish between high- and low-risk groups of OVCF to formulate corresponding treatment plans. When a patient’s lumbar spine FF value is close to or above the cut-off value, intervention at multiple levels such as lifestyle, diet, exercise, and medication is necessary, with particular attention focused on preventing falls or injuries to avoid or delay the occurrence of OVCF.

Several limitations must be acknowledged: (1) The cross-sectional design precludes causal inferences about FF and fracture risk; (2) Generalizability may be limited by the homogeneous sample (older adult(s) Asian population) - validation in diverse ethnic groups and younger patients is needed (32); (3) Manual ROI segmentation, while clinically practical, introduces measurement variability versus automated methods; (4) The exact temporal relationship between fracture onset and FF measurement remains unclear;(5)The sample size is relatively small, and larger cohort studies are needed to validate and revise the results of this study; (6) while MRI was performed within one week of diagnosis, the exact interval between fracture onset and imaging was not recorded. This temporal relationship is important to confirm that FF changes preceded rather than resulted from the fracture. (7) Our study used DXA, which is routinely used in clinical practice to assess BMD, but a study has shown that standardized QCT was more accurate (26). Future research should focus on optimizing the measurement techniques to improve the accuracy of the data; (8)The ROIs were manually determined, which means that the shapes may not be in perfect agreement and may lead to errors. Therefore, several measurements were performed in this study, and the average was used as the final result.(9) While we controlled for age statistically, the modest sample size precluded meaningful age-stratified analyses. Future larger studies should examine age-specific FF thresholds.

Conclusion

In summary, our study shows that FF values of lumbar spine measured using the q-Dixon sequence are helpful in predicting the risk of vertebral compression fractures in osteoporosis patients. While promising, the q-Dixon sequence may become a practical tool for non-invasive assessment of fracture risk, though its clinical implementation requires validation in larger, prospective studies with standardized protocols across diverse populations.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by Medical Ethics Committee of Affiliated Hospital of Qingdao University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individuals for the publication of any potentially identifiable images or data included in this article.

Author contributions

JZ: Conceptualization, Data curation, Validation, Visualization, Writing – original draft. QL: Data curation, Formal analysis, Software, Visualization, Writing – original draft. YW: Investigation, Methodology, Resources, Validation, Writing – review & editing. LS: Conceptualization, Project administration, Supervision, Visualization, Writing – review & editing. QZ: Formal analysis, Investigation, Methodology, Visualization, Writing – review & editing. CG: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The authors declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

We thank Jenna MacArthur, PhD, from Liwen Bianji (Edanz) (www.liwenbianji.cn) for editing the English text of a draft of this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Abbreviations

AUC, Area under the curve; BMD, Bone mineral density; Dixon, Chemical shift-encoding water-fat imaging; DXA, Dual energy X-ray absorptiometry; FF, Fat fraction; ICC, Intraclass correlation coefficient; IP, In-phase; MRI, Magnetic resonance imaging; MRS, Magnetic resonance spectroscopy; MSC, Marrow mesenchymal stem cell; OP, Opposed-phase; OVCF, Osteoporotic vertebral compression fracture; PDFF, Proton density fat fraction; QCT, Quantitative computed tomography; ROC, Receiver operating characteristic curve; ROI, Region of interest.

References

1. Kutsal FY and Ergin Ergani GO. Vertebral compression fractures: Still an unpredictab le aspect of osteoporosis. Turk J Med Sci. (2021) 51:393–9. doi: 10.3906/sag-2005-315

PubMed Abstract | Crossref Full Text | Google Scholar

2. Sambrook P and Cooper C. Osteoporosis. Lancet. (2006) 367:2010–8. doi: 10.1016/S0140-6736(06)68891-0

PubMed Abstract | Crossref Full Text | Google Scholar

3. Al-Saleh Y, Sulimani R, Sabico S, Raef H, Fouda M, Alshahrani F, et al. Guidelines for osteoporosis in Saudi Arabia: recommendations from the saudi osteoporosis society. Ann Saudi Med. (2015) 35:1–12. doi: 10.5144/0256-4947.2015.1

PubMed Abstract | Crossref Full Text | Google Scholar

4. Link TM. Osteoporosis imaging: state of the art and advanced imaging. Radiology. (2012) 263:3–17. doi: 10.1148/radiol.12110462

Crossref Full Text | Google Scholar

5. Gimble JM, Robinson CE, Wu X, and Kelly KA. The function of adipocytes in the bone marrow stroma: an update. Bone. (1996) 19:421–8. doi: 10.1016/S8756-3282(96)00258-X

PubMed Abstract | Crossref Full Text | Google Scholar

6. Dixon WT. Simple proton spectroscopic imaging. Radiology. (1984) 153:189–94. doi: 10.1148/radiology.153.1.6089263

PubMed Abstract | Crossref Full Text | Google Scholar

7. Glover GH and Schneider E. Three-point Dixon technique for true water/fat decomposition with B0 inhomogeneity correction. Magn Reson Med. (1991) 18:371–83. doi: 10.1002/mrm.1910180211

PubMed Abstract | Crossref Full Text | Google Scholar

8. Hupfeld S, Pischel D, Jechorek D, Janicová A, Pech M, and Fischbach F. MRI-based fat quantification of the liver: Is it time for commercially available products? Eur J Radiol. (2021) 144:109993. doi: 10.1016/j.ejrad.2021.109993

PubMed Abstract | Crossref Full Text | Google Scholar

9. Black DM, Bauer DC, Vittinghoff E, Lui L-Y, Grauer A, Marin F, et al. Treatment-related changes in bone mineral density as a surrogate biomarker for fracture risk reduction: meta-regression analyses of individual patient data from multiple randomised controlled trials. Lancet Diabetes Endocrinol. (2020) 8:672–82. doi: 10.1016/S2213-8587(20)30159-5

PubMed Abstract | Crossref Full Text | Google Scholar

10. Kung AWC, Lee K-K, Ho AY, Tang G, and Luk KD. Ten-year risk of osteoporotic fractures in postmenopausal Chinese women according to clinical risk factors and BMD T-scores: a prospective study. J Bone Miner Res. (2007) 22:1080–7. doi: 10.1359/jbmr.070320

PubMed Abstract | Crossref Full Text | Google Scholar

11. Schuit SCE, van der Klift M, Weel AE a. M, de Laet CEDH, Burger H, Seeman E, et al. Fracture incidence and association with bone mineral density in older adult(s) men and women: the Rotterdam Study. Bone. (2004) 34:195–202. doi: 10.1016/j.bone.2003.10.001

PubMed Abstract | Crossref Full Text | Google Scholar

12. Cranney A, Jamal SA, Tsang JF, Josse RG, and Leslie WD. Low bone mineral density and fracture burden in postmenopausal women. CMAJ. (2007) 177:575–80. doi: 10.1503/cmaj.070234

PubMed Abstract | Crossref Full Text | Google Scholar

13. Stone KL, Seeley DG, Lui L-Y, Cauley JA, Ensrud K, Browner WS, et al. BMD at multiple sites and risk of fracture of multiple types: long-term results from the Study of Osteoporotic Fractures. J Bone Miner Res. (2003) 18:1947–54. doi: 10.1359/jbmr.2003.18.11.1947

PubMed Abstract | Crossref Full Text | Google Scholar

14. Fazeli PK, Horowitz MC, MacDougald OA, Scheller EL, Rodeheffer MS, Rosen CJ, et al. Marrow fat and bone–new perspectives. J Clin Endocrinol Metab. (2013) 98:935–45. doi: 10.1210/jc.2012-3634

PubMed Abstract | Crossref Full Text | Google Scholar

15. Paccou J, Hardouin P, Cotten A, Penel G, and Cortet B. The role of bone marrow fat in skeletal health: usefulness and perspectives for clinicians. J Clin Endocrinol Metab. (2015) 100:3613–21. doi: 10.1210/jc.2015-2338

PubMed Abstract | Crossref Full Text | Google Scholar

16. Aparisi Gómez MP, Ayuso Benavent C, Simoni P, Aparisi F, Guglielmi G, and Bazzocchi A. Fat and bone: the multiperspective analysis of a close relationship. Quant Imaging Med Surg. (2020) 10:1614–35. doi: 10.21037/qims.2020.01.11

PubMed Abstract | Crossref Full Text | Google Scholar

17. Zhao Y, Huang M, Ding J, Zhang X, Spuhler K, Hu S, et al. Prediction of abnormal bone density and osteoporosis from lumbar spine MR using modified dixon quant in 257 subjects with quantitative computed tomography as reference. J Magn Reson Imaging. (2019) 49:390–9. doi: 10.1002/jmri.26233

PubMed Abstract | Crossref Full Text | Google Scholar

18. Alyas F, Saifuddin A, and Connell D. MR imaging evaluation of the bone marrow and marrow infiltrative disorders of the lumbar spine. Magn Reson Imaging Clin N Am. (2007) 15:199–219, vi. doi: 10.1016/j.mric.2007.03.002

PubMed Abstract | Crossref Full Text | Google Scholar

19. Wang D, Haile A, and Jones LC. Dexamethasone-induced lipolysis increases the adverse effect of adipocytes on osteoblasts using cells derived from human mesenchymal stem cells. Bone. (2013) 53:520–30. doi: 10.1016/j.bone.2013.01.009

PubMed Abstract | Crossref Full Text | Google Scholar

20. Li G, Tang G, Liu Y, Tang R, Peng Y, and Li W. MR spectroscopy and micro-CT in evaluation of osteoporosis model in rabbits: comparison with histopathology. Eur Radiol. (2012) 22:923–9. doi: 10.1007/s00330-011-2325-x

PubMed Abstract | Crossref Full Text | Google Scholar

21. Li G-W, Chang S-X, Fan J-Z, Tian Y-N, Xu Z, and He Y-M. Marrow adiposity recovery after early zoledronic acid treatment of glucocorticoid-induced bone loss in rabbits assessed by magnetic resonance spectroscopy. Bone. (2013) 52:668–75. doi: 10.1016/j.bone.2012.11.002

PubMed Abstract | Crossref Full Text | Google Scholar

22. Liu J, Chen J-D, Li P, Liao J-W, Feng J-X, Chen Z-Y, et al. Comprehensive assessment of osteoporosis in lumbar spine using compositional MR imaging of trabecular bone. Eur Radiol. (2023) 33:3995–4006. doi: 10.1007/s00330-022-09368-0

PubMed Abstract | Crossref Full Text | Google Scholar

23. Tang R, Tang G, Hua T, Tu Y, Ji R, and Zhu J. mDIXON-Quant technique diagnostic accuracy for assessing bone mineral density in male adult population. BMC Musculoskelet Disord. (2023) 24:125. doi: 10.1186/s12891-023-06225-z

PubMed Abstract | Crossref Full Text | Google Scholar

24. Bensreti H, Alhamad DW, Gonzalez AM, Pizarro-Mondesir M, Bollag WB, Isales CM, et al. Update on the role of glucocorticoid signaling in osteoblasts and bone marrow adipocytes during aging. Curr Osteoporos Rep. (2023) 21:32–44. doi: 10.1007/s11914-022-00772-5

PubMed Abstract | Crossref Full Text | Google Scholar

25. Ambrosi TH and Schulz TJ. The emerging role of bone marrow adipose tissue in bone health and dysfunction. J Mol Med (Berl). (2017) 95:1291–301. doi: 10.1007/s00109-017-1604-7

PubMed Abstract | Crossref Full Text | Google Scholar

26. Roski F, Hammel J, Mei K, Baum T, Kirschke JS, Laugerette A, et al. Bone mineral density measurements derived from dual-layer spectral CT enable opportunistic screening for osteoporosis. Eur Radiol. (2019) 29:6355–63. doi: 10.1007/s00330-019-06263-z

PubMed Abstract | Crossref Full Text | Google Scholar

27. Karampinos DC, Melkus G, Baum T, Bauer JS, Rummeny EJ, and Krug R. Bone marrow fat quantification in the presence of trabecular bone: initial comparison between water-fat imaging and single-voxel MRS. Magn Reson Med. (2014) 71:1158–65. doi: 10.1002/mrm.24775

PubMed Abstract | Crossref Full Text | Google Scholar

28. Gasser FT, Kufner A, Gassert FG, Leondhart Y, Kronthaler S, Schwaiger BJ, et al. MR-based proton density fat fraction (PDFF) of the vertebral bone marrow differentiates between patients with and without osteoporotic vertebral fractures. Osteoporosis international: J established as result cooperation between Eur Foundation Osteoporosis Natl Osteoporosis Foundation USA. (2022) 33(2):487–96. doi: 10.1007/s00198-021-06147-3

PubMed Abstract | Crossref Full Text | Google Scholar

29. Liney GP, Bernard CP, Manton DJ, Turnbull LW, and Langton CM. Age, gender, and skeletal variation in bone marrow composition: a preliminary study at 3. 0 Tesla. J Magn Reson Imaging. (2007) 26:787–93. doi: 10.1002/jmri.21072

PubMed Abstract | Crossref Full Text | Google Scholar

30. Li X, Kuo D, Schafer AL, Porzig A, Link TM, Black D, et al. Quantification of vertebral bone marrow fat content using 3 Tesla MR spectroscopy: reproducibility, vertebral variation, and applications in osteoporosis. J Magn Reson Imaging. (2011) 33:974–9. doi: 10.1002/jmri.22489

PubMed Abstract | Crossref Full Text | Google Scholar

31. Yun JS, Lee H-D, Kwack K-S, and Park S. Use of proton density fat fraction MRI to predict the radiographic progression of osteoporotic vertebral compression fracture. Eur Radiol. (2021) 31:3582–9. doi: 10.1007/s00330-020-07529-7

PubMed Abstract | Crossref Full Text | Google Scholar

32. Şahiner Z, Güner M, and Uysal ÖS. Frailty as a predictor of vertebral fractures in osteoporosis: A clinical perspective. Jt Dis Relat Surg. (2025) 36:408–14. doi: 10.52312/jdrs

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: magnetic resonance imaging, osteoporosis, vertebral compression fracture, fat fraction, bone mineral density

Citation: Zhang J, Li Q, Wang Y, Sun L, Zhang Q and Gao C (2025) The q-Dixon sequence for MRI predicts osteoporotic vertebral compression fractures. Front. Endocrinol. 16:1591044. doi: 10.3389/fendo.2025.1591044

Received: 20 March 2025; Accepted: 01 August 2025;
Published: 25 August 2025.

Edited by:

Guanwu Li, Shanghai University of Traditional Chinese Medicine, China

Reviewed by:

Dongliang Zhang, Beijing Jishuitan Hospital, China
Tatsuya Tanaka, International University of Health and Welfare (IUHW), Japan

Copyright © 2025 Zhang, Li, Wang, Sun, Zhang and Gao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chuanping Gao, Z2FvY2h1YW5waW5nQHFkdS5lZHUuY24=

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.