- 1Department of Radiology, Klinikum Rechts der Isar, School of Medicine, Technical University of Munich, Munich, Germany
- 2Department on Neuroradiology, Klinikum Rechts der Isar, School of Medicine, Technical University of Munich, Munich, Germany
- 3Department of Neuroradiology, University Hospital of Munich (LMU), Munich, Germany
Background: Quantitative magnetic resonance imaging (MRI) techniques such as chemical shift encoding-based water-fat separation techniques (CSE-MRI) are increasingly applied as noninvasive biomarkers to assess the biochemical composition of vertebrae. This study aims to investigate the longitudinal change of proton density fat fraction (PDFF) and T2* derived from CSE-MRI of the thoracolumbar vertebral bone marrow in patients that develop incidental vertebral compression fractures (VCFs), and whether PDFF and T2* enable the prediction of an incidental VCF.
Methods: In this study we included 48 patients with CT-derived bone mineral density (BMD) measurements at baseline. Patients that presented an incidental VCF at follow up (N=12, mean age 70.5 ± 7.4 years, 5 female) were compared to controls without incidental VCF at follow up (N=36, mean age 71.1 ± 8.6 years, 15 females). All patients underwent 3T MRI, containing a significant part of the thoracolumbar spine (Th11-L4), at baseline, 6-month and 12 month follow up, including a gradient echo sequence for chemical shift encoding-based water-fat separation, from which PDFF and T2* maps were obtained. Associations between changes in PDFF, T2* and BMD measurements over 12 months and the group (incidental VCF vs. no VCF) were assessed using multivariable regression models. Mixed-effect regression models were used to test if there is a difference in the rate of change in PDFF, T2* and BMD between patients with and without incidental VCF.
Results: Prior to the occurrence of an incidental VCF, PDFF in vertebrae increased in the VCF group (ΔPDFF=6.3 ± 3.1%) and was significantly higher than the change of PDFF in the group without VCF (ΔPDFF=2.1 ± 2.5%, P=0.03). There was no significant change in T2* (ΔT2*=1.7 ± 1.1ms vs. ΔT2*=1.1 ± 1.3ms, P=0.31) and BMD (ΔBMD=-1.2 ± 11.3mg/cm3 vs. ΔBMD=-11.4 ± 24.1mg/cm3, P= 0.37) between the two groups over 12 months. At baseline, no significant differences were detected in the average PDFF, T2* and BMD of all measured vertebrae (Th11-L4) between the VCF group and the group without VCF (P=0.66, P=0.35 and P= 0.21, respectively). When assessing the differences in rates of change, there was a significant change in slope for PDFF (2.32 per 6 months, 95% confidence interval (CI) 0.31-4.32; P=0.03) but not for T2* (0.02 per 6 months, CI -0.98-0.95; P=0.90) or BMD (-4.84 per 6 months, CI -23.4-13.7; P=0.60).
Conclusions: In our study population, the average change of PDFF over 12 months is significantly higher in patients that develop incidental fractures at 12-month follow up compared to patients without incidental VCF, while T2* and BMD show no significant changes prior to the occurrence of the incidental vertebral fractures. Therefore, a longitudinal increase in bone marrow PDFF may be predictive for vertebral compression fractures.
Introduction
In our rapidly aging population, vertebral compression fractures (VCFs) are increasing in frequency due to reduced bone strength as a result of osteoporosis (1). Osteoporosis is a systemic disease that leads to reduced bone mass and microarchitectural deterioration of bone tissue and as a result, low-energy (or even no evident) trauma can lead to VCFs, which affects up to 25% of all postmenopausal women (2–4). Assessment of the bone is of high importance, since reduced bone strength can be treated, reducing the risk for fractures and following complications such as pain and immobilisation (5), being a burden for the individual patient as well as for the health care and social security systems (6, 7).
In the clinical routine, quantitative assessment of the bone structure predominantly relies on dual-energy X-ray absorptiometry (DXA) and dedicated quantitative CT (QCT), which are the clinical standards for determining the bone mineral density (BMD) (6). However, quantitative MR techniques have been introduced and investigated for analysis of vertebral bone marrow, such as chemical shift encoding-based water–fat MRI (CSE-MRI). CSE-MRI is a non-invasive approach to extract precise quantitative information about water–fat compositions of the vertebral bone marrow by determining the proton density fat fraction (PDFF), showing very good concordance with histology or magnetic resonance spectroscopy (8–12). CSE-MRI therefore enables in vivo quantitative evaluation of the bone marrow without exposing the patient to ionizing radiation. PDFF measurements of the vertebral bone marrow have proved to be particularly interesting in the analysis of osteoporosis, as the composition of the vertebral bone marrow, consisting of a mixture of hematopoietic tissue and adipocytes interspersed within the trabecular bone matrix, shifts towards a higher fat content (12, 13). Previous studies demonstrate that PDFF derived by CSE-MRI represents a potential biomarker for quantification of osteoporosis, bone fragility, and assessment of osteoporotic fractures (14–16), without exposing the patient to radiation when compared to DXA and QCT. Those techniques also enable the assessment of bone marrow T2* (3, 17–19), which has been shown to correlate with the microarchitectural structure of the vertebral body (20).
However, all previous studies were cross-sectional studies that only included patients without fractures or patients already suffering from VCFs. To the best of our knowledge, there is no study that has analyzed the intraindividual dynamic of the chemical composition in the vertebral bone marrow using chemical shift encoding-based water-fat separation techniques yet, especially prior to the occurrence of a VCF. It would be of great interest, whether this technique allows for prognostic assessment in regard to risk evaluation for future VCFs. We hypothesize that there might be measurable changes in PDFF or T2* in the vertebral bone marrow prior to the occurrence of a VCF, which, as a result, could identify patients at risk and prevent imminent VCFs.
Therefore, the purpose of this study was to investigate the association between incidental VCFs and longitudinal changes in PDFF and T2* measurements prior to the occurrence of vertebral fractures in the thoracolumbar spine.
Methods
Patient selection and study design
In our institution, from June 2018 to December 2021, a total number of 478 patients received an abdominal MR imaging protocol including a multi-echo gradient-echo sequence of the thoracolumbar spine, as described below. The inclusion criteria were: i) occurrence of an incidental VCF that was ii) not present in the 2 MRI examinations prior to the VCF and iii) without evidence for a high-energy traumatic or malignant etiology (i.e. osseous metastasis or primary bone tumor). Ultimately, 12 patients (70.5 ± 7.4 years, 5 females) with incidental vertebral fractures were identified that met all the criteria (Figure 1). Those patients were frequency-matched (1:3) for age and sex to patients without incidental VCFs (n=36, 71.1 ± 8.6 years, 15 females). History of medical treatment was checked in all included patients. None of the patients had received chemotherapy prior to or during the study. All patients underwent 3T-MRI of the abdomen including axially acquired gradient echo sequences for chemical shift encoding-based water-fat separation, from which PDFF and T2* maps of the thoracolumbar spine were obtained. The first imaging session was defined as the baseline MR imaging and follow-up sessions were performed every 6 months (6.7 ± 3.4 months). The study was approved by the local institutional review board (Ethics Commission of the Medical Faculty, Technical University of Munich, Germany; Ethics proposal number 2022-433-S-SR).
Magnetic resonance imaging and T2*/PDFF assessment
Two 3T-MRI systems (one Ingenia, Philips Healthcare and one Elition, Philips Healthcare) were used for the examination of the abdomen including the lumbar spine. The patients were placed in supine position and the combination of a 16-channel torso coil array and an inbuilt table posterior 12-channel coil array was used. The imaging sessions were scheduled according to the clinical follow-up, which includes CSE-MRI for PDFF and T2* measurement. For the PDFF and T2* measurements a six-echo 3D multi-echo gradient-echo sequence was used acquiring all echoes in a single TR using bipolar gradients with the following parameters: repetition time TR/first echo time TE1/echo time step ΔTE = 7.8/1.35/1.1ms, field of view (FOV)=300 x 400 x 150 mm3, acquisition voxel size = 2 x 3 x 6 mm3, reconstruction voxel size = 1.13 x 1.13 x 6 mm3, receiver bandwidth = 1678 Hz/pixel, frequency direction = anterior/posterior (A/P), 1 average, scan time = 9.3 s. To minimize T1-bias effects a flip angle of 3° was used (21).
Complex multi-echo gradient-echo images were provided as input to the fat quantification routine provided by the vendor (mDixon Quant, Philips Healthcare). Specifically, after phase correction, a complex-based water-fat decomposition was performed using a single T2* correction and a pre-calibrated fat spectrum accounting for the presence of the multiple peaks in the fat spectrum. A seven-peak fat spectrum model was employed (22). The PDFF map was computed as the ratio of the fat signal over the sum of fat and water signals. PDFF and T2* maps were extracted.
Furthermore, the MR imaging protocol encompasses an axial and coronal T2-TSE sequence as well as axially acquired native and contrast-enhanced T1 weighted sequences with spectral fat saturation according to clinical standards.
Segmentations were performed on axial reformations (Figure 2). The vertebral bodies Th11 to L4 were manually segmented by F.T.G. and Y.L. (3 and 4 years of experience in musculoskeletal imaging). For this, cylindrical ROIs were placed in the center of each vertebra and the mean PDFF and T2* values were noted for each vertebra. The measurements from Th11 to L4 were then averaged for each examination. Fractured or otherwise altered vertebrae (e.g. after kyphoplasty, heavy degenerative alterations; n=18) were excluded.
Figure 2 Axially acquired gradient echo sequences for chemical shift encoding-based water-fat separation of a 78 year old male patient at baseline, first follow-up and at the time of fracture occurrence; the images on the right originate from the corresponding survey in sagittal orientation, which reveal an incidental VCF of the base plate of L3 (indicated by the white arrowhead). The circles centrally placed in the vertebrae represent ROIs used for the PDFF measurements.
The morphologic imaging features of the vertebral fractures (involvement of the superior and/or inferior endplate; deformity (e.g. wedge, biconcave); involvement of the posterior column; Genant classification (mild, moderate, severe)) were assessed by two board-certified radiologists (B.J.S. and A.S.G., with 11 years and 10 years of experience in musculoskeletal imaging). The radiologists ensured, in context with the clinical history, that no morphologic signs for a malignant fracture were present in any of the sequences.
In order to assess the intra- and inter-reader reproducibility error of the PDFF, T2* and BMD values, a random sample of 10 subjects was chosen and reanalyzed by the same radiologist.
QCT and BMD measurements
BMD measurements were performed in patients that received CT in a time period ± 1 month around the MR imaging sessions. CT images were acquired with one dual-layer dual-energy CT scanner (IQon Spectral CT, Philips Healthcare) and one multislice detector CT scanner (iCT 256, Philips Healthcare). CT was performed according to our routine clinical protocols: collimation, 0.6 mm; pixel spacing, 0.4/0.3 mm; pitch factor, 0.8/0.9; tube voltage (peak), 120 kV; modulated tube current, 102–132 mA. Images were reformatted in 3-mm slice thickness using a bone-specific convolution kernel.
For BMD measurements, a mid-line 15mm MPR section in sagittal reformations was created with a PACS tool (IDS7, Sectra). Cylindrical volumes of interest were manually positioned in the non-fractured lumbar vertebrae by one radiologist (G.C.F.), and mean Hounsfield Units (HU) were noted (23). Fractured or otherwise altered vertebrae (e.g. vertebrae with severe degenerative changes, vertebrae after vertebro-/kyphoplasty) were excluded. The HU values were then converted into BMD using conversion equations as previously described: (i) 0.928 x HU + 4.5 mg/cm3 for the IQon Spectral CT, (ii) 0.855 x HU + 1.172 mg/cm3 for the Philips iCT 256 (24, 25). Osteoporosis was defined as BMD < 80 mg/cm3 and osteopenia as 80 ≤ BMD ≤ 120 mg/cm3 (26).
Statistical analysis
Statistical analysis was performed with SPSS version 28 software (IBM) using a two-sided 0.05 level of significance. One-way analysis of variance (for parametric variables) and chi-square tests (for categorical variables) were used to evaluate differences in patient characteristics between patients with and without incidental fractures. Associations between changes in PDFF, T2* and BMD measurements over 12 months and the group (incidental VCF vs. no incidental VCF) were assessed using multivariable regression models adjusting for age and sex. Moreover in order to assess the validity of our hypothesis, we used mixed-effect regression models with random intercepts to test if there is a difference in the rate of change in PDFF, T2* and BMD between patients with and without incidental VCF.
Intrarater and interrater reproducibility for T2*, PDFF and trabecular BMD values were assessed by calculating the intraclass correlation coefficient and the root mean square coefficient of variation (RMSCV) of the differences between the respective measurements.
Results
In this study, 12 patients (70.5 ± 7.4 years, 5 female) with incidental VCFs were identified and were frequency-matched (1:3) for age and sex to subjects without vertebral fractures (n=36, 71.1 ± 8.6 years, 15 females) (Table 1). Of all incidental VCFs, 16.7% were located in Th12, 25% in L1, 16.7% in L2, 33.3% in L3 and 8.3% in L4.
PDFF and T2* measurements
No significant difference was detected in the average PDFF of all measured vertebrae (Th11-L4) between the group with incidental VCFs and the group without VCF at baseline, 6-month and 12-month follow-up (baseline, PDFF=45.6 ± 12.6% in the VCF group vs. PDFF=47.4 ± 10.5% in the no-VCF group; 6-month follow up, PDFF=50.0 ± 13.3% vs. PDFF=48.1 ± 10.1%; 12-month follow-up, PDFF=53.4 ± 8.7% vs. PDFF=48.6 ± 10.1%; P=0.66, P=0.65 and P=0.24 respectively) (Table 2). Moreover, there was no significant difference in T2* measurements (baseline, T2*=9.7 ± 2.1ms vs. T2*=9.1 ± 2.0ms; 6-month follow-up, T2*=10.4 ± 3.0ms vs. T2*=9.1 ± 1.8ms; 12-month follow-up T2*=10.9 ± 1.5ms vs. T2*=9.7 ± 1.8ms; P=0.35, P=0.10 and P=0.58 respectively) (Figure 3).
Table 2 PDFF, T2* and BMD values at baseline, 6-month follow-up and 12-month follow-up/time of fracture.
Figure 3 Average PDFF, T2* and BMD measurements in the groups with and without VCFs at baseline, 6-month and 12-month follow-up.
When analyzing the change of PDFF between the baseline MR examination and the 12-month follow-up, there was an average increase in vertebral PDFF of ΔPDFF=6.3 ± 3.1% in the group with incidental VCFs, which is significantly higher (P=0.03) than the change in PDFF in the group without VCFs (2.1 ± 2.5%, Figure 4; Table 2). However, there was no significant difference regarding the changes of vertebral T2* between the group with and without incidental VCF at follow up (ΔT2*=1.7 ± 1.1ms vs. ΔT2*=1.1 ± 1.3ms, P=0.31). When assessing the differences in rates of change of the PDFF and T2* values over 12 months between patients with and without incidental VCFs, there was a significant change in slope for PDFF (2.32 per 6 months, 95% confidence interval (CI) 0.31-4.32; P=0.03) but not for T2* (0.02 per 6 months, CI -0.98-0.95; P=0.90).
Figure 4 Average change of PDFF, T2* and BMD in patients with and without incidental VCF between baseline and 12-month follow-up.
BMD measurements
BMD of the spine was measured using opportunistic QCT, if CT imaging was performed in proximity to the MR examinations. There was no significant difference between the two groups regarding the BMD measurements at baseline, 6-month or 12-month follow-up (baseline, 114.1 ± 28.8mg/cm3 vs. 129.5 ± 22.8mg/cm3; 6-month follow-up, 109.8 ± 8.4mg/cm3 vs. 122.9 ± 15.8mg/cm3; 12-month follow-up, 103.4 ± 25.5mg/cm3 vs. 118.2 ± 30.1mg/cm3; P=0.21, P=0.43 and P=0.34 respectively). Moreover, when assessing the differences in BMD change over 12 months, there was no significant difference between the group with and without incidental VCF (-1.2 ± 11.3 mg/cm3 vs. -11.4 ± 24.1 mg/cm3, P=0.37). When assessing the differences in rates of change of the BMD values over 12 months between patients with and without incidental VCFs, there was no significant change in slope (-4.84 per 6 months, CI -23.4-13.7; P=0.60).
Also, according to the above-described data, there was no indication for a higher prevalence of osteoporosis in either group. At baseline, one of the patients with an incidental VCF (8.3%) showed an osteoporotic BMD (<80mg/cm3), whereas 3 patients in the group without fracture had BMD values that qualified for osteoporosis.
Intrareader reproducibility
The intrareader (ICC, 0.99 [95% CI, 0.98–0.99]) and interreader (ICC, 0.99 [95% CI, 0.97–0.99]) agreement for trabecular BMD as well as the intrareader and interreader agreement for PDFF (ICC for both, 0.98 [95% CI, 0.96–0.99]) and T2* measurements (ICC for both, 0.98 [95% CI, 0.96–0.99]) were excellent. Intrarater reproducibility, calculated with the RMSCV, was 0.4% for the BMD analysis and 0.9% for the mean PDFF analysis and 0.8% for the mean T2* analysis. The interreader agreement for PDFF, calculated with RMSCV, was 0.8%, the interreader agreement for T2* was 0.9% and the respective interreader agreement for the trabecular BMD was 0.5%.
Discussion
In this study, we investigated longitudinal changes over 12 months of PDFF and T2* measurements using chemical shift encoding-based water-fat separation for the assessment of patients with and without incidental vertebral compression fractures. Our study demonstrates that the average change of PDFF over 12 months is significantly higher in patients that develop incidental fractures at 12-month follow up compared to patients without incidental VCF, while T2* and BMD show no significant changes prior to the occurrence of the incidental vertebral fractures. Our data suggests that PDFF, when evaluated longitudinally in the same patient, might indicate changes in the vertebral structural integrity and composition very early and therefore might be a sensitive biomarker for the identification of patients at risk for incidental VCFs.
Past studies have shown that vertebral bone marrow PDFF might be a useful parameter for assessment of osteoporosis and reduced bone strength (14–16). PDFF derived from chemical shift encoding-based water-fat separation measures a map of the density of hydrogen protons attributable to fat, normalized by the total hydrogen proton density from all mobile proton species (27). PDFF enables a fairly precise estimation of the fat volume fraction as a result of the almost equal relative proton densities of fat and water (8). Therefore, it has been proposed to be a non-invasive tool for assessment of e.g. hepatic steatosis and, as previous studies demonstrate, osteoporosis or low-energy fractures (13, 18, 28, 29).
In general, osteoporosis is characterized by a reduced bone mass and microarchitectural deterioration of bone tissue. The lost bone mass in the vertebral space is not just filled with fatty bone marrow but rather includes mechanisms such as a drift in mesenchymal stem cell differentiation towards adipogenesis over osteoblastogenesis (30, 31). Multiple studies have shown that there is a negative correlation between PDFF values and BMD in vertebral bone marrow; as a result this technique enables a differentiation of patients with osteoporosis and healthy subjects (32). Another study demonstrated that PDFF allows the differentiation between osteoporotic patients with and without vertebral fractures, representing a potentially useful radiation-free tool for the assessment of bone fragility in addition to BMD measurements (16). CSE-MRI also enables T2* measurements, which have shown to be indirect measures of vertebral bone architecture, as the local magnetic field inhomogeneities induced by the trabecular bone and bone marrow interface can be measured as a shortening of the effective transverse relaxation time. T2* measurements correlated with the bone mineral density in previous studies, where T2* increases with decreasing bone density (19). However, all previous studies were cross-sectional studies that only included patients without fractures or patients already suffering from VCFs. To the best of our knowledge, there is no study that has analyzed the intraindividual dynamic of the composition of vertebral bone marrow using chemical shift encoding-based water-fat separation techniques in MR imaging yet, especially prior to the occurrence of a VCF.
In past studies opportunistic QCT was evaluated for patients at risk for recurring osteoporotic VCFs (24, 33, 34). In our patient cohort however, BMD measured by QCT and T2* were not able to predict occurrence of incidental fractures longitudinally. Also, according to our data, patients with incidental fractures do not differ from subjects without VCFs in regard to the overall average PDFF. A tendency is observable in the fracture group with a rise in average PDFF, even though it is not statistically different from the non-fracture group. Only the relative change of PDFF showed a significant difference between patients with imminent fractures compared to patients without fractures. We therefore hypothesize that PDFF may be a very sensitive biomarker that may be useful for the prediction of the occurrence of incidental fractures. Yet larger longitudinal study cohorts are needed in order to investigate these results further.
One important limitation of our study is the relatively small number of patients with VCFs, although being a frequent condition especially in older patients. However, this problem is difficult to eliminate since the occurrence of a VCF in a patient population that is not explicitly at risk for VCFs beforehand is still a stochastic event; a very large number of patients is necessary that repeatedly receives MR examinations until a sufficiently large patient number with VCFs is gathered. Another limitation is the incomplete data set for BMD measurements using QCT; as the patients already received MR examinations there was often no justifying indication for CT, which would lead to unnecessary radiation exposure. It is therefore necessary to confirm our findings in future studies with a larger data set.
In conclusion, in our study population the average change of PDFF was significantly higher in patients that developed incidental fractures compared to patients without incidental VCFs, while T2* and BMD showed no significant changes. Therefore, when evaluated longitudinally, PDFF might indicate changes in the vertebral structural integrity and composition very early and might be a sensitive biomarker for the identification of patients at risk for incidental VCFs. Yet, this needs to be further assessed in larger longitudinal study groups.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Commission of the Medical Faculty, Technical University of Munich, Germany (Ethics proposal number 2022-433-S-SR). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
Conceptualization, AG, YL; methodology, AG, YL, TB, JSK, DK; validation, TB, JSK, FTG, BS, SR, DK; formal analysis, YL, JK, LG, FTG, FGG, GF; investigation, YL, JK, LG, FTG, FGG, GF; resources, AG, RB, DK, MM; data curation, YL, JK, LG, FTG, FGG, GF; writing—original draft preparation, YL; writing—review and editing, all authors; visualization, YL, JK; supervision, AG; project administration, YL, AG; funding acquisition, AG, MM. All authors contributed to the article and approved the submitted version.
Conflict of interest
JSK is Co-Founder of BoneScreen GmbH.
The remaining 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.
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.
References
1. Melton LJ 3rd. Epidemiology of spinal osteoporosis. Spine (1997) 22(24):2S–11S. doi: 10.1097/00007632-199712151-00002
2. Mauch JT, Carr CM, Cloft H, Diehn FE. Review of the imaging features of benign osteoporotic and malignant vertebral compression fractures. AJNR Am J neuroradiology (2018) 39(9):1584–92. doi: 10.3174/ajnr.A5528
3. Kühn J-P, Hernando D, Meffert PJ, Reeder S, Hosten N, Laqua R, et al. Proton-density fat fraction and simultaneous R2* estimation as an MRI tool for assessment of osteoporosis. Eur Radiol (2013) 23(12):3432–39. doi: 10.1007/s00330-013-2950-7
4. Consensus development conference: Diagnosis, prophylaxis, and treatment of osteoporosis. Am J Med (1993) 94(6):646–50. doi: 10.1016/0002-9343(93)90218-e
5. Miller PD. Management of severe osteoporosis. Expert Opin Pharmacother (2016) 17(4):473–88. doi: 10.1517/14656566.2016.1124856
6. Hernlund E, Svedbom A, Ivergård M, Compston J, Cooper C, Stenmark J, et al. Osteoporosis in the European union: medical management, epidemiology and economic burden. A report prepared in collaboration with the international osteoporosis foundation (IOF) and the European federation of pharmaceutical industry associations (EFPIA). Arch Osteoporos (2013) 8(1-2):136. doi: 10.1007/s11657-013-0136-1
7. Nazrun AS, Tzar MN, Mokhtar SA, Mohamed IN. A systematic review of the outcomes of osteoporotic fracture patients after hospital discharge: morbidity, subsequent fractures, and mortality. Ther Clin Risk Manag (2014) 10:937–48. doi: 10.2147/TCRM.S72456
8. Reeder SB, Hu HH, Sirlin CB. Proton density fat-fraction: A standardized MR-based biomarker of tissue fat concentration. J Magn Reson Imaging (2012) 36(5):1011–4. doi: 10.1002/jmri.23741
9. Hu HH, Kan HE. Quantitative proton MR techniques for measuring fat. NMR BioMed (2013) 26(12):1609–29. doi: 10.1002/nbm.3025
10. Sollmann N, Löffler MT, Kronthaler S, Böhm C, Dieckmeyer M, Ruschke S, et al. MRI-Based quantitative osteoporosis imaging at the spine and femur. J Magnetic Resonance Imaging (2021) 54(1):12–35. doi: 10.1002/jmri.27260
11. Baum T, Yap SP, Dieckmeyer M, Ruschke S, Eggers H, Kooijman H, et al. Assessment of whole spine vertebral bone marrow fat using chemical shift-encoding based water-fat MRI. J Magn Reson Imaging (2015) 42(4):1018–23. doi: 10.1002/jmri.24854
12. Karampinos DC, Ruschke S, Dieckmeyer M, Diefenbach M, Franz D, Gersing AS, et al. Quantitative MRI and spectroscopy of bone marrow. J Magn Reson Imaging (2018) 47(2):332–53. doi: 10.1002/jmri.25769
13. Griffith JF, Yeung DKW, Antonio GE, Lee FKH, Hong AWL, Wong SYS, et al. Vertebral bone mineral density, marrow perfusion, and fat content in healthy men and men with osteoporosis: Dynamic contrast-enhanced MR imaging and MR spectroscopy. Radiology (2005) 236(3):945–51. doi: 10.1148/radiol.2363041425
14. Dieckmeyer M, Ruschke S, Cordes C, Yap SP, Kooijman H, Hauner H, et al. The need for T₂ correction on MRS-based vertebral bone marrow fat quantification: implications for bone marrow fat fraction age dependence. NMR BioMed (2015) 28(4):432–9. doi: 10.1002/nbm.3267
15. Karampinos DC, Melkus G, Baum T, Bauer JS, Rummeny EJ, Krug R. Bone marrow fat quantification in the presence of trabecular bone: initial comparison between water-fat imaging and single-voxel MRS. Magnetic resonance Med (2014) 71(3):1158–65. doi: 10.1002/mrm.24775
16. Gassert FT, Kufner A, Gassert FG, Leonhardt 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 Int 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
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(2):390–99. doi: 10.1002/jmri.26233
18. 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 shift-based water-fat separation. Clin Radiol (2014) 69(3):254–62. doi: 10.1016/j.crad.2013.10.005
19. Wu H-Z, Zhang X-F, Han S-M, Cao L, Wen J-X, Wu W-J, et al. Correlation of bone mineral density with MRI T2* values in quantitative analysis of lumbar osteoporosis. Arch Osteoporos (2020) 15(1):18. doi: 10.1007/s11657-020-0682-2
20. Leonhardt Y, Gassert FT, Feuerriegel G, Gassert FG, Kronthaler S, Boehm C, et al. Vertebral bone marrow T2* mapping using chemical shift encoding-based water-fat separation in the quantitative analysis of lumbar osteoporosis and osteoporotic fractures. Quantitative Imaging Med Surg (2021) 11(8):3715–25. doi: 10.21037/qims-20-1373
21. Karampinos DC, Yu H, Shimakawa A, Link TM, Majumdar S. T1-corrected fat quantification using chemical shift-based water/fat separation: Application to skeletal muscle. Magn Reson Med (2011) 66(5):1312–26. doi: 10.1002/mrm.22925
22. Ren J, Dimitrov I, Sherry AD, Malloy CR. Composition of adipose tissue and marrow fat in humans by 1H NMR at 7 Tesla. J Lipid Res (2008) 49(9):2055–62. doi: 10.1194/jlr.D800010-JLR200
23. Baum T, Müller D, Dobritz M, Rummeny EJ, Link TM, Bauer JS. BMD measurements of the spine derived from sagittal reformations of contrast-enhanced MDCT without dedicated software. Eur J Radiol (2011) 80(2):e140–5. doi: 10.1016/j.ejrad.2010.08.034
24. Löffler MT, Jacob A, Valentinitsch A, Rienmüller A, Zimmer C, Ryang YM, et al. Improved prediction of incident vertebral fractures using opportunistic QCT compared to DXA. Eur Radiol (2019) 29(9):4980–89. doi: 10.1007/s00330-019-06018-w
25. 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(11):6355–63. doi: 10.1007/s00330-019-06263-z
26. American College of Radiology. ACR–SPR–SSR practice parameter for the performance of musculoskeletal quantitative computed tomography (QCT) (2018). Available at: https://wwwacrorg/-/media/ACR/Files/Practice-Parameters/qctpdf?la=en (Accessed 15 October 2020).
27. Reeder SB, Sirlin CB. Quantification of liver fat with magnetic resonance imaging. Magn Reson Imaging Clin N Am (2010) 18(3):337–57. doi: 10.1016/j.mric.2010.08.013
28. Tang GY, Lv ZW, Tang RB, Liu Y, Peng YF, Li W, et al. Evaluation of MR spectroscopy and diffusion-weighted MRI in detecting bone marrow changes in postmenopausal women with osteoporosis. Clin Radiol (2010) 65(5):377–81. doi: 10.1016/j.crad.2009.12.011
29. 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 magnetic resonance Imaging JMRI (2011) 33(4):974–79. doi: 10.1002/jmri.22489
30. Dunnill MS, Anderson JA, Whitehead R. Quantitative histological studies on age changes in bone. J Pathol Bacteriol (1967) 94(2):275–91. doi: 10.1002/path.1700940205
31. Karampinos DC, Ruschke S, Gordijenko O, Grande Garcia E, Kooijman H, Burgkart R, et al. Association of MRS-based vertebral bone marrow fat fraction with bone strength in a human In vitro model. J Osteoporos (2015) 2015:152349. doi: 10.1155/2015/152349
32. Cheng X, Li K, Zhang Y, Wang L, Xu L, Liu Y, et al. The accurate relationship between spine bone density and bone marrow in humans. Bone (2020) 134:115312. doi: 10.1016/j.bone.2020.115312
33. Leonhardt Y, May P, Gordijenko O, Koeppen-Ursic VA, Brandhorst H, Zimmer C, et al. Opportunistic QCT bone mineral density measurements predicting osteoporotic fractures: A use case in a prospective clinical cohort. Front Endocrinol (2020) 11:586352–52. doi: 10.3389/fendo.2020.586352
Keywords: osteoporosis, magnetic resonance imaging, bone density, PDFF, chemical shift encoding-based water-fat separation
Citation: Leonhardt Y, Ketschau J, Ruschke S, Gassert FT, Glanz L, Feuerriegel GC, Gassert FG, Baum T, Kirschke JS, Braren RF, Schwaiger BJ, Makowski MR, Karampinos DC and Gersing AS (2022) Associations of incidental vertebral fractures and longitudinal changes of MR–based proton density fat fraction and T2* measurements of vertebral bone marrow. Front. Endocrinol. 13:1046547. doi: 10.3389/fendo.2022.1046547
Received: 16 September 2022; Accepted: 02 November 2022;
Published: 17 November 2022.
Edited by:
Xingming Shi, Augusta University, United StatesReviewed by:
Grzegorz Tatoń, Jagiellonian University, PolandRon Noah Alkalay, Beth Israel Deaconess Medical Center and Harvard Medical School, United States
Copyright © 2022 Leonhardt, Ketschau, Ruschke, Gassert, Glanz, Feuerriegel, Gassert, Baum, Kirschke, Braren, Schwaiger, Makowski, Karampinos and Gersing. 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: Yannik Leonhardt, yannik.leonhardt@tum.de