The diagnostic value of morphological features of fat deposition of sacroiliac joint steatosis in axial spondyloarthritis

Abstract

Background:

Findings of fatty lesions in the context of other imaging manifestations, especially bone marrow edema and erosions can effectively assist in the diagnosis of axSpA. Chemical shift-encoded MRI is a sequence which allows for the quantification of fat signal and has been applied in the imaging evaluation of the SIJ in axSpA. The objective of this study was to investigate the diagnostic performance of morphological features of fatty lesions visualized by CSE-MRI in the imaging evaluation of SIJ in axSpA.

Methods:

Fatty lesions with morphological features (subchondral, homogeneity and distinct border) were assessed and recorded as a binary variable in each quadrant of the SIJ. Sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) were calculated for different morphological features as well as the anatomical distribution in patients with nr-axSpA and r-axSpA. T1-weighted images and CSE-MRI fat fraction maps were directly compared in the recognition of different morphological features.

Results:

Eighty-two patients [non-SpA (n = 21), nr-axSpA (n = 23), r-axSpA (n = 38)] with lower back pain (LBP) were enrolled. Presence of the three morphological features of fatty lesions had a specificity of 90.48% in axSpA. The sensitivities of being subchondral, homogeneity and distinct border were 52.17, 39.13 and 39.13% in nr-axSpA on T1-weighted images. For patients with r-axSpA, the sensitivities reached 86.84, 76.32 and 57.89%. No significant difference was found in the distribution of fatty lesions between T1-weighted images and CSE-MRI. However, CSE-MRI fat fraction maps could detect significantly more fatty lesions with homogeneity (p = 0.0412) and distinct border (p = 0.0159) than T1-weighted images in the sacroiliac joint, but not subchondral lesions (p = 0.6831).

Conclusion:

The homogeneity and distinct border are more relevant for the diagnosis of axSpA. Moreover, CSE-MRI could detect more typical morphological features of fatty lesions than T1-weighted images in showing these two features. The presence of all three features was more likely to be indicative of axSpA.

Introduction

As per the standardized definitions of MRI lesions in the sacroiliac joint of patients with axial spondyloarthritis by the ASAS MRI working group, fat metaplasia belongs in the group of structural lesions (1). Although fatty lesion is a common finding on the T1-weighted images, its relevance to the diagnosis of axSpA has been widely debated (2). According to a large population-based cohort study investigating frequency of MRI changes, fatty lesion is also very common in healthy individuals, hence its diminished specificity (3). Findings of fatty lesions are advised to be interpreted in the context of other imaging manifestations, especially bone marrow edema and erosion (4).

Three morphological features of fat infiltration have been defined so as to facilitate better characterization of such fatty lesions: (1) a distinct border around the fatty lesion; (2) homogeneity of the increased T1-weighted signal; (3) proximity of the fatty lesion to subchondral bone of SIJ (1). Previous studies investigating the diagnostic utility of these features concluded that SIJ fatty lesion per se does not facilitate the recognition of nr-axSpA, while distinct border or homogeneity showed small to moderate diagnostic utility, which was strongly associated with concomitant bone marrow edema and erosion (5, 6).

Chemical shift-encoded MRI is a sequence which allows for the quantification of fat signal (7–9) and has been applied in the imaging evaluation of the SIJ in axSpA (10). Previous studies demonstrated that the SIJ fat fraction is significantly higher in patients with axSpA, especially in r-axSpA (2, 11). Utilizing the quantification of fat fraction, CSE-MRI could generate fat fraction maps, which could be employed as a promising tool to visualize the morphological features of fatty lesions. This study aims to investigate the diagnostic utility of morphological features of fatty lesions visualized by CSE-MRI in the imaging evaluation of SIJ in axSpA.

Methods

Subjects

Participants complaining of lower back pain (LBP) who visited the outpatient clinic of the Department of Rheumatology of the Third Affiliated Hospital of Sun Yat-sen University from June 1, 2020 to January 20, 2022 were enrolled in this study. The inclusion criteria consisted of: (1) Chronic lower back pain ≥3 months; (2) Aged ≤50 years. All participants were required to meet: (a) BMI ≤30 kg/m2; (b) No history of other rheumatic immune disorders; (c) no history of malignant neoplasms; (d) no contraindications for MRI. Clinical data such as age, sex, age of onset, duration of disease, past medical history, smoking history, history of medication was collected with questionnaires. Disease activity was assessed using the Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) and the Ankylosing Spondylitis Disease Activity Score – CRP (ASDAS-CRP).

The diagnosis of each participant was determined by the expert panel, which consisted of two rheumatologists and one radiologist. Based on the clinical information and the pelvic radiographs, study subjects were classified as non-SpA (n = 21), nr-axSpA (n = 23), r-axSpA (n = 38). Patients determined as non-SpA were regarded as the control group. Grading of the radiographic sacroiliitis was conducted by the radiologist according to the modified New York criteria. Patients with previous or current treatment with biologics were excluded from this study. This study was conducted abiding by the Declaration of Helsinki and was approved by the Ethics committee of the Third Affiliated Hospital of Sun Yat-sen University. Informed consent was obtained from all the participants.

Evaluation of MR images

MRI images were acquired on 3.0 T SignaTM Architect (GE). T1 weighted, T1 fat-saturated (T1-FS) sequences and T2 fat-saturated (T2-FS) sequences in a semi-coronal orientation and T2 FS sequences in a semi-axial orientation were obtained. Scoring of fat deposition and erosions was based on T1 weighted images. The SIJ MRI scans were evaluated independently by two observers in a random order on the workstation of the institution. All the scans were evaluated by both observers twice in order to assess the intra-reader agreement. Fatty lesions with various morphological features were determined when both observers unanimously agreed on the presence of the lesions.

The following morphological features were assessed across the MRI slices: (1) Subchondral: whether the fatty lesion is located in the subchondral area under the cartilage; (2) Homogeneity: whether the signal within the fatty lesion is homogeneous; (3) Distinct border: whether the fatty lesion has a distinct contour. Assessment of the morphological features were recorded as a binary variable in each quadrant of the SIJ. Prior to the imaging evaluation, a training SIJ MRI set was established in order to reach consensus between both observers.

Statistical analysis

Software packages R version 3.6.3 (The R Project for Statistical Computing, Vienna, Austria) was used for all statistical analysis. Characteristics of the study subjects were summarized using descriptive statistics. The mean ± SD or the median and IQR were presented according to the distribution of the individual characteristics. We computed the proportions of patients with fatty lesions in at least 1, 2, 3, 4 quadrants of SIJ in different groups. Frequencies of the three morphological features were also calculated in each group.

Sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) were calculated for different morphological features as well as the anatomical distribution. Such analyses were conducted in patients with nr-axSpA and r-axSpA separately. T1-weighted images and CSE-MRI fat fraction maps were directly compared in the recognition of different morphological features. Images were defined as T1(+)CSE-MRI(+), T1(+)CSE-MRI(+), T1(−)CSE-MRI(+), T1(−)CSE-MRI(−) based on the features seen on the two sequences. The McNemar’s test was used to analyze whether there was a difference between the two sequences. Unweighted Cohen’s Kappa was calculated to compare the inter-reader and intra-reader agreement of the T1-weighted images and CSE-MRI fat fraction maps. A p value <0.05 was considered significant in all analyses.

Results

Demographic characteristics

This study included 82 patients with lower back pain. 21 patients were classified as non-SpA, while 23 were nr-axSpA and 38 were r-axSpA. Demographic and clinical characteristics were listed in Table 1. Compared with non-SpA patients, a significantly higher proportion of patients with axSpA were male. Patients with r-axSpA had a significantly longer disease duration compared with non-SpA and nr-axSpA (median 5.5 years vs. 2 years and 2 years, p < 0,001). Patients with r-axSpA also had a higher disease activity than nr-axSpA, as indicated by the ASDAS-CRP (median 2.09 vs. 1.27, p = 0.008).

Distribution of fat deposition in SIJ and its morphological features

On the T1-weighted images, the frequencies of being subchondral, homogeneity and distinct border in non-SpA patients were 23.81, 9.52 and 14.29% respectively, while the frequencies were 52.17, 39.13, 39.13% in nr-axSpA and 86.84, 76.32 and 57.89% in r-axSpA.

On the CSE-MRI fat fraction maps, the frequencies of such morphological features were 23.81, 14.29, 14.29% in non-SpA patients. For nr-axSpA patients with r-axSpA, the frequencies were 51.17, 47.83, 52.17 and 92.11%, 84.21, 73.68%, respectively.

The proportions of non-SpA patients with at least 1, 2 and 3 morphological features seen on T1-weighted images were 23.81, 14.29 and 9.52%, respectively, while in nr-axSpA and r-axSpA patients, the proportions were 65.22, 39.13, 26.09 and 89.47%, 76.32, 55.26%, respectively. In contrast, on CSE-MRI images, the proportions of non-SpA patients were 23.81, 14.29 and 14.29%. The proportions were 65.22, 52.17, 34.78 and 92.11%, 84.21 73.68% in nr-axSpA and r-axSpA patients, respectively (Shown in Table 2).

Diagnostic utilities of morphological features

The sensitivity, specificity, PPV and NPV of various morphological features pertinent to fat deposition on MRI in the diagnosis of nr-axSpA and r-axSpA were shown in Tables 3, 4. The sensitivities of being subchondral, homogeneity and distinct border were modest in nr-axSpA on T1-weighted images, reaching 52.17, 39.13 and 39.13%, respectively. For patients with r-axSpA, the sensitivities reached 86.84, 76.32 and 57.89%. Fatty lesions with all three features could only be observed in 26.09% of the nr-SpA patients, while such lesions could be found in 55.26% of r-axSpA patients. However, such typical lesions could be very specific to axSpA, with a specificity of 90.48%.

Comparison of T1-weighted images with CSE-MRI fat fraction maps

The difference between fatty lesions on T1-weighted images and CSE-MRI fat fraction maps is shown in Figure 1.

Figure 1

In order to clarify whether there was a difference between T1-weighted images and CSE-MRI fat fraction maps in the recognition of morphological features of fatty lesions, McNemar’s test was applied. Results of the comparisons could be seen in Table 5. No significant difference was found in the distribution of fatty lesions between the two MRI sequences. However, CSE-MRI fat fraction maps could detect significantly more fatty lesions with homogeneity (p = 0.0412) and distinct border (p = 0.0159) than T1-weighted images in the sacroiliac joint, but not subchondral lesions (p = 0.6831).

Results of the inter-reader agreement of the fatty lesions as well as the morphological features could be seen in Table 6, while the intra-reader agreement could be seen in Supplementary Table S1.

Discussion

Fat deposition per se has limited value for nr-axSpA patients so that it should be interpreted in conjunction with other sacroiliac joint lesions to improve the diagnostic accuracy of nr-axSpA (1). We found that there were only 39.13% of nr-axSpA patients had homogeneous or well-circumscribed fatty lesions in the T1-weighted images of nr-axSpA patients, while the feature of being located in subchondral bone could reach 52.17%. The frequency of fat deposition with one feature, two features, and three features were 65.22, 39.13, and 26.09%, respectively. Similarly, Weber et al. found that 20–30% of patients exhibited fatty lesions with three morphological features in the sacroiliac joint.

Noteworthy, fat infiltration could also be observed in the sacroiliac joints in non-SpA patients. In the T1-weighted images, the proportions of non-SpA patients with subchondral, homogeneous and clearly demarcated fatty lesions were 23.81, 9.52 and 14.29%, respectively. Meanwhile, the frequencies of fat deposition accompanied by one feature, two features, and three features were 23.81, 14.29, and 9.52% respectively, which were lower than SpA patients. Further analysis of diagnostic utility showed that the specificities of homogeneity and clear boundaries were the highest, reaching 90.48 and 85.71%, while being subchondral was only 76.19%, which reflected that the diagnostic value of being subchondral is inferior to homogeneity and distinct border.

Our results suggested that more morphological features and more quadrants with fatty lesions could improve the specificity of fat infiltration in axSpA. Fatty lesions in r-axSpA tend to have a higher fat fraction and be more extensive in the sacroiliac joint. This is due to the fact that patients with r-axSpA usually have longer disease durations and more severe sacroiliac joint destruction. Other than number of quadrants, fat deposition in r-axSpA exhibits more morphological features. 55.26% of the r-axSpA patients exhibited fatty lesions with all three morphological features. Based on these results, typical fatty lesions with all three morphological features are much more common in r-axSpA, but relatively rare in nr-axSpA patients, since the structural destruction in nr-axSpA patients tend to be mild. In clinical practice, r-axSpA patients could easily be identified with plain radiographs, since the structural changes in the sacroiliac joints could be conspicuous. Therefore, fatty lesions, even with the typical morphological features, are unlikely to offer much incremental diagnostic value. However, our results did confirm that the three morphological features, especially homogeneity and distinct border, were rather specific to axSpA, albeit less common in nr-axSpA.

The strength of CSE-MRI lies in its ability to quantify fat fractions in the region of interest (7), yet the fat fraction map generated by this sequence could also be a powerful tool to visualize the otherwise inconspicuous morphological features (12). This is the first study to utilize the fat fraction map of CSE-MRI to visualize the morphological features of fatty lesions in axSpA in comparison with traditional T1-weighted images. Although our study failed to prove that CSE-MRI could visualize more fatty lesions, observers could detect more typical morphological features in the fatty lesion, especially the feature of homogeneity and distinct border. Interestingly, such improvement in the observation of morphological features does not incur a loss in the specificity. Therefore, a CSE-MRI fat fraction map could be a helpful addition to the observation of fatty lesions in the sacroiliac joints. CSE-MRI is currently considered the state-of-the-art technique in the quantification of fat signal. Compared with other sequences such as T1-in-and-out-of-phase (IOP) MRI, CSE-MRI is a more sophisticated fat quantification approach which corrects for a number of confounding factors including T1 bias, noise bias and eddy currents. The accuracy of the fat quantification by CSE-MRI could be attributed to such intricate correction calculations, thus generating fat fraction maps enabling better observation of morphological features.

Other than CSE-MRI, a couple of other MRI sequences could also facilitate better characterization of fatty lesions, Dixon sequences rely on the difference in resonance frequency between water and fat, unlike fat suppression sequences. Hence, Dixon sequences could generate water-only, fat-only, in-phase and out-of-phase images in a single acquisition (13, 14). Kuetting et al. used multipoint Dixon technique to non-invasively distinguish chylous and nonchylous ascites and pleural effusions (15). This study showed that fat fraction was greater for chylous versus nonchylous fluids, and in vivo fat fraction was correlated with triglyceride content (15). Sequences such as volumetric interpolated breath-hold examination with fat saturation (VIBE-FS) are also promising techniques for the observation of fatty lesions (16, 17). In the future, head-to-head comparisons can be conducted to directly compare different fat-sensitive sequences in the detection and quantification of fatty lesions in SIJ.

It should be pointed out that this study had some limitations. There might be systemic bias in the study design as observers may have overestimated the diagnostic performance of CSE-MRI fat fraction maps in visualizing the morphological features and underestimated in T1-weighted images. To minimize this bias, all T1-weighted images and CSE-MRI fat fraction maps were presented in an independent, random manner, and when there was inconsistency between the two observers, a senior radiologist made the final decision. Another limitation to this study is the relatively small sample size. A cohort of 82 patients with lower back pain may not fully present the potential of CSE-MRI fat fraction maps in the detection of fatty lesions. It would be desirable to include a larger number of patients, as well as acquiring Dixon and VIBE with fat saturation images in the future.

Conclusion

Fatty lesions with morphological features of being subchondral, homogeneity and distinct border could only be found in 52.17, 39.13 and 39.13% of nr-axSpA patients. Of the three morphological features, homogeneity and distinct border are more specific to axSpA as compared to being subchondral. Typical lesions with all three morphological features could be found in 26.09 and 55.26% of nr-axSpA and r-axSpA patients, yet such lesions could provide a specificity of 90.48%. CSE-MRI fat fraction map is superior to T1-weighted images in the recognition of fatty lesions with homogeneity and distinct border, but could not visualize more fatty lesions.

Funding

This work was supported by the grants from the National Natural Science Foundation of China (81871294); the Science and Technology Planning Project of Guangdong Province, China (2019B030316004); Guangdong Clinical Research Center of Immune diseases (2020B1111170008).

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.

  Maksymowych WP Lambert RG Ostergaard M Pedersen SJ Machado PM Weber U et al . MRI lesions in the sacroiliac joints of patients with spondyloarthritis: an update of definitions and validation by the ASAS MRI working group. Ann Rheum Dis. (2019) 78:1550–8. doi: 10.1136/annrheumdis-2019-215589

  • CrossRef
  • Google Scholar

• 2.

  Chen M Yu K Hu X Van Den Berghe T Qi Y Jin B et al . Proton density fat fraction (PDFF) maps aid fat metaplasia evaluation in the sacroiliac joints in additional to T1WI: improved diagnostic accuracy in axial spondyloarthritis. Eur J Radiol. (2022) 157:110569. doi: 10.1016/j.ejrad.2022.110569

  • CrossRef
  • Google Scholar

• 3.

  Weber U Zubler V Pedersen SJ Rufibach K Lambert RG Chan SM et al . Development and validation of a magnetic resonance imaging reference criterion for defining a positive sacroiliac joint magnetic resonance imaging finding in spondyloarthritis. Arthritis Care Res. (2013) 65:977–85. doi: 10.1002/acr.21893

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  Weber U Pedersen SJ Zubler V Rufibach K Chan SM Lambert RG et al . Fat infiltration on magnetic resonance imaging of the sacroiliac joints has limited diagnostic utility in nonradiographic axial spondyloarthritis. J Rheumatol. (2014) 41:75–83. doi: 10.3899/jrheum.130568

  • CrossRef
  • Google Scholar

• 5.

  de Hooge M van den Berg R Navarro-Compán V Reijnierse M van Gaalen F Fagerli K et al . Patients with chronic back pain of short duration from the SPACE cohort: which MRI structural lesions in the sacroiliac joints and inflammatory and structural lesions in the spine are most specific for axial spondyloarthritis?Ann Rheum Dis. (2016) 75:1308–14. doi: 10.1136/annrheumdis-2015-207823

  • CrossRef
  • Google Scholar

• 6.

  Baraliakos X Ghadir A Fruth M Kiltz U Redeker I Braun J . Which magnetic resonance imaging lesions in the sacroiliac joints are most relevant for diagnosing axial spondyloarthritis? A prospective study comparing rheumatologists' evaluations with radiologists' findings. Arthritis Rheumatol. (2021) 73:800–5. doi: 10.1002/art.41595

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  Eggers H Börnert P . Chemical shift encoding-based water-fat separation methods. J Magn Reson Imaging. (2014) 40:251–68. doi: 10.1002/jmri.24568

  • CrossRef
  • Google Scholar

• 8.

  Wang X Colgan TJ Hinshaw LA Roberts NT Bancroft LCH Hamilton G et al . T(1) -corrected quantitative chemical shift-encoded MRI. Magn Reson Med. (2020) 83:2051–63. doi: 10.1002/mrm.28062

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  Kühn JP Hernando D Mensel B Krüger PC Ittermann T Mayerle J et al . Quantitative chemical shift-encoded MRI is an accurate method to quantify hepatic steatosis. J Magn Reson Imaging. (2014) 39:1494–501. doi: 10.1002/jmri.24289

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  Ren C Zhu Q Yuan H . Mono-exponential and bi-exponential model-based diffusion-weighted MR imaging and IDEAL-IQ sequence for quantitative evaluation of sacroiliitis in patients with ankylosing spondylitis. Clin Rheumatol. (2018) 37:3069–76. doi: 10.1007/s10067-018-4321-x

  • CrossRef
  • Google Scholar

• 11.

  Liu D Lin C Liu B Qi J Wen H Tu L et al . Quantification of fat metaplasia in the sacroiliac joints of patients with axial spondyloarthritis by chemical shift-encoded MRI: a diagnostic trial. Front Immunol. (2021) 12:811672. doi: 10.3389/fimmu.2021.811672

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  Navaratna R Zhao R Colgan TJ Hu HH Bydder M Yokoo T et al . Temperature-corrected proton density fat fraction estimation using chemical shift-encoded MRI in phantoms. Magn Reson Med. (2021) 86:69–81. doi: 10.1002/mrm.28669

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  Kirchgesner T Acid S Perlepe V Lecouvet F Vande Berg B . Two-point Dixon fat-water swapping artifact: lesion mimicker at musculoskeletal T2-weighted MRI. Skelet Radiol. (2020) 49:2081–6. doi: 10.1007/s00256-020-03512-x

  • CrossRef
  • Google Scholar

• 14.

  Giza SA Miller MR Parthasarathy P de Vrijer B McKenzie CA . Comparison of modified two-point Dixon and chemical shift encoded MRI water-fat separation methods for fetal fat quantification. J Magn Reson Imaging. (2018) 48:274–82. doi: 10.1002/jmri.25929

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  Kuetting D Luetkens J Fimmers R Sprinkart AM Attenberger U Pieper CC . MRI assessment of chylous and nonchylous effusions: use of multipoint Dixon fat quantification. Radiology. (2020) 296:698–705. doi: 10.1148/radiol.2020200089

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  Michaely HJ Morelli JN Budjan J Riffel P Nickel D Kroeker R et al . CAIPIRINHA-Dixon-TWIST (CDT)-volume-interpolated breath-hold examination (VIBE): a new technique for fast time-resolved dynamic 3-dimensional imaging of the abdomen with high spatial resolution. Investig Radiol. (2013) 48:590–7. doi: 10.1097/RLI.0b013e318289a70b

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  Song Q Zeng M Chen C Ma J Yun H Rao S et al . Non-contrast-enhanced magnetic resonance angiography using T2-weighted 3-dimensional fat-suppressed turbo spin echo (SPACE): diagnostic performance and comparison with contrast-enhanced magnetic resonance angiography using volume interpolated breath-hold examination in the detection of portosystemic and portohepatic collaterals. J Comput Assist Tomogr. (2012) 36:675–80. doi: 10.1097/RCT.0b013e31826999df

  • Pubmed Abstract
  • CrossRef
  • Google Scholar