0847

Differences in distribution of MRI-based fat fraction in lower limb skeletal muscles of six different neuromuscular disorders
Harmen Reyngoudt1,2, Pierre-Yves Baudin1,2, Ericky C.A. Araujo1,2, Pierre G. Carlier3, and Benjamin Marty1,2
1NMR Laboratory, Neuromuscular Investigation Center, Institute of Myology, Paris, France, 2NMR Laboratory, CEA/DRF/IBFJ/MIRCen, Paris, France, 3CEA, DRF, Service Hospitalier Frédéric Joliot, Orsay, France

Synopsis

Quantitative MRI fat-water separation techniques such as Dixon are often used to evaluate disease progression in muscle of neuromuscular diseases. Although the mean fat fraction (FF) value per region of interest (ROI) is a valuable objective MRI biomarker, it does not reflect the variation of FF within this ROI. In this study we analyzed 6582 muscle ROIs in leg and thigh of 6 different neuromuscular diseases and analyzed, besides the mean FF, other statistical metrics such as median, standard deviation, kurtosis and skewness. The differences in FF distribution might reveal additional information about the individual patient’s disease evolution.

Introduction

Magnetic resonance imaging (MRI) is a powerful tool in the evaluation of the status and the progression of muscle destruction in neuromuscular diseases (NMDs). In the last five years, many studies in NMDs have implemented quantitative MRI-based fat-water separation techniques such as the Dixon method1. Following segmentation of the FF maps, a mean FF value per muscle, muscle group or global segment can then be reported2. Although the mean FF value per region of interest (ROI) is a valuable objective MRI biomarker, it does not reflect the variation of FF within this ROI. Indeed, with these quantitative methods, less focus is put on the actual images. Image texture analysis algorithms have been subject to research in the field of muscle imaging to elucidate whether there is additional information in the patterns that are less or not visible to the human eye3. The aim of this work was to look further than the mean FF value per ROI and to apply other simple statistical metrics, as a way to investigate the FF heterogeneity in different NMDs.

Methods - Database

For this study, we investigated 216 leg and 210 thigh MRI data sets (left and right sides), obtained in 6 different NMDs (Duchenne muscular dystrophy/DMD: n=36 for legs; dysferlinopathy/DYS: n= 43 for legs and thighs; facioscapulohumeral dystrophy/FSHD: n=21 for thighs; inclusion body myositis/IBM: n=80 for legs and thighs; immune-mediated necrotizing myositis/IMNM: n=43 for legs and thighs; late-onset Pompe disease: n=14 for legs and thighs.

Methods - Data acquisition

All data were obtained on a 3T clinical Trio/Prisma Siemens system. The quantitative MRI protocol included a 3D gradient echo sequence (Dixon) with TEs of 2.75, 3.95 and 5.15 ms, a TR of 10 ms, a flip angle of 3°, a spatial resolution of 1x1x5 mm3 and a 448x224x64 matrix size4.

Methods - Data processing/analysis

Manual segmentation was performed in 14 leg muscles and 22 thigh muscles (Fig. 1), on 5 slices, avoiding the muscle borders (fasciae, intermuscular and subcutaneous fat). The central slice was always positioned at the same anatomical level5. From the generated FF maps4, the analysis per ROI included the assessment of mean (M), median (Mdn), standard deviation (std), kurtosis (K=3=normal distribution; K<3: flatter than normal distribution; K>3: more peaked than normal distribution) and skewness (-0.5<S>0.5: normal distribution; S<-0.5: negative skew, M<Mdn; S>0.5: positive skew, M>Mdn) of the FF distribution. Histograms were generated for all ROIs. As a first analysis, we categorized ROIs into 4 FF categories: FF≤10%, 10%<FF≤30%, 30%<FF≤60% and FF>60%, which correspond to Lamminen-Mercuri (LM) scores 1, 2, 3 and 4, respectively6, which is a pseudo-quantitative scoring system based on a visual radiological assessment of T1-weighted MRI to roughly determine the degree of muscle fat replacement. Then, we also looked into smaller FF categories (e.g. 30-35%, 45-50%, 55-60%). Statistical analysis (Kruskal-Wallis) included the comparison of the FF distribution parameters between the different NMDs, per FF category (significance level: P<.005).

Results

Overall, the FF metrics were evaluated in 6582 leg and thigh muscles. For ROIs with FF≤10%, FF histograms were, logically, highly positively skewed and showed very high K values (>>3) for all 6 NMDs. Highly significant differences were observed for all metrics between (almost) all NMDs. The ROIs with FF of 10-30% were also positively skewed and showed K-values between 1 and 2.5. Here, only significant differences were found for M, Mdn and S between DYS and other NMDs. The ROIs with FF of 30-60% showed no significant differences between any of the NMDs for M and Mdn values. The significant differences found in this FF category were found for std and K, especially between DMD or DYS and the other NMDs; and histograms were rather flat (K<3) but not skewed such as the other FF categories (Fig. 2A). The ROIs with FF>60% showed predominantly differences in std between NMDs, particularly between DMD or DYS and the other NMDs (Fig. 2B). When looking more closely to the interesting FF category of 30-60%, first, at the soleus of 4 selected patients a highly similar mean FF (Fig. 3), we observed that the FF histograms corresponded to the ones in the overall analysis (Fig. 2A), except for IBM which demonstrates a highly particular FF distribution in this particular example. The next step was then to look at the FF distribution (between 30-60% FF) in all soleus muscles for the different NMDs (Fig. 4). Again we observed that the same pattern was seen for DMD and DYS, and somewhat less for IBM, as compared to the overall analysis (Fig. 2A). We see that the differences between DMD and the other NMDs are also evident when we zoom into smaller FF categories (Fig. 5).

Discussion

We observed that the muscles that are moderately affected (FF=30-60%) showed mean FF values that were the same between the investigated NMDs and that the differences in FF between these diseases could be attributed to other parameters such as std and K. Overall, the analysis of FF distribution and the heterogeneity of muscle fat replacement might reveal additional information about the individual patient’s disease evolution (i.e. natural history) and the effects of treatment7.

Acknowledgements

No acknowledgement found.

References

1Carlier PG, Marty B, Scheidegger O, et al. Skeletal muscle quantitative nuclear magnetic resonance imaging and spectroscopy as an outcome measure for clinical trials. J Neuromuscul Dis. 2016;3(1):1-28.

2Reyngoudt H, Marty B, Boisserie JM, et al. Global versus individual muscle segmentation to assess quantitative MRI-based fat fraction changes in neuromuscular diseases. Eur Radiol. 2020; doi: 10.1007/s00330-020-07487-0.

3De Certaines JD, Larcher T, Duda D, et al. Application of texture analysis to muscle MRI: 1-What kind of information should be expected from texture analysis? EPJ Nonlinear Biomed Phys. 2015;1-3:1-14.

4Azzabou N, Loureiro de Sousa P, Caldas E, Carlier PG. Validation of a generic approach to muscle water T2 determination at 3T in fat-infiltrated skeletal muscle. J Magn Reson Imaging. 2015;41(3):645-653.

5Gidaro T, Reyngoudt H, Le Louër J, et al. Quantitative nuclear magnetic resonance imaging detects subclinical changes over 1 year in skeletal muscle of GNE myopathy. J Neurol. 2020;267(1):228-238.

6Baudin PY, Marty B, Robert B, et al. Qualitative and quantitative evaluation of skeletal muscle fatty degenerative changes using whole-body Dixon nuclear magnetic resonance imaging for an important reduction of the acquisition time. Neuromuscul Disord. 2015;25(10):758-763.

7Chen G, Jiang J, Wang X, et al. Evaluation of hepatic steatosis before liver transplantation in ex vivo by volumetric quantitative PDFF-MRI. Magn Reson Med. 2020; doi: 10.1002/mrm.28592.

Figures

Fig. 1. Illustration of the individually drawn muscle ROIs across 5 slices in the thigh and the leg: AL, adductor longus, AM, adductor magnus, BF, biceps femoris, ED, extensor digitorum, GL, gastrocnemius lateralis, GM, gastrocnemius medialis, GRA, gracilis, PER, peroneus, RF, rectus femoris, SAR, sartorius, SM, semimembranosus, SOL, soleus, ST, semitendinosus, TA, tibialis anterior, TP, tibialis posterior, VI, vastus intermedius, VL, vastus lateralis, VM, vastus medialis.

Fig. 2. FF histograms for all muscles (leg and thigh) in 6 NMDs (DMD, DYS, FSHD, Pompe, IMB and IMNM). (A) ROIs with FF=30-60%. (B) ROIs with FF>60%. Values for number of pixels in ROI (#pixels), mean FF (M), median FF (Mdn), FF standard deviation (std), kurtosis (K) and skewness (S) are presented as mean ± standard deviation, of all included muscles or ROIs (indicated by n). Post-hoc statistical analysis: indices reflect significant differences (P<.005) with other diseases (1=DMD, 2=DYS, 3=FSHD, 4=Pompe, 5=IBM, 6=IMNM).

Fig. 3. FF histograms and MRI-based FF map in 4 examples of soleus muscle with very similar mean FF value, in 4 different NMDs (DMD, DYS, IMB and IMNM). Values for number of pixels in ROI (#pixels), mean FF (M), median FF (Mdn), FF standard deviation (std), kurtosis (K) and skewness (S) are presented as mean ± standard deviation, of all included muscles or ROIs (indicated by n).

Fig. 4. FF histograms for all soleus muscles in 4 NMDs (DMD, DYS, IMB and IMNM), with FF=30-60%. Values for number of pixels in ROI (#pixels), mean FF (M), median FF (Mdn), FF standard deviation (std), kurtosis (K) and skewness (S) are presented as mean ± standard deviation, of all included muscles or ROIs (indicated by n). Post-hoc statistical analysis: indices reflect significant differences (P<.005) with other diseases (1=DMD, 2=DYS, 5=IBM, 6=IMNM)

Fig. 5. FF histograms for all all muscles in 4 NMDs (DMD, DYS, IMB and IMNM), in 3 FF categories (30-35%, 45-50% and 55-60%). Values for number of pixels in ROI (#pixels), mean FF (M), median FF (Mdn), FF standard deviation (std), kurtosis (K) and skewness (S) are presented as mean ± standard deviation, of all included muscles or ROIs (indicated by n). Post-hoc statistical analysis: indices reflect significant differences (P<.005) with other diseases (1=DMD, 2=DYS, 5=IBM, 6=IMNM).

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
0847