A simplified method to determine tissue-water T2 from CPMG image data in fat infiltrated skeletal muscle: application in the forearm in Duchenne muscular dystrophy
Nick Zafeiropoulos1, Valeria Ricotti2, Matthew Evans1,3, Jasper Morrow3, Paul Matthews4, Robert Janiczek5, Tarek Yousry1,3, Christopher Sinclair1,3, Francesco Muntoni2, and John Thornton1,3

1Neuroradiological Academic Unit, UCL Institute of Neurology, London, United Kingdom, 2Dubowitz Neuromuscular Centre, UCL Institute of Child Health, London, United Kingdom, 3MRC Centre for Neuromuscular Diseases, London, United Kingdom, 4Imperial College London, London, United Kingdom, 5GlaxoSmithKline, London, United Kingdom

Synopsis

A simplified CPMG signal decay model was used to determine muscle-water T2 (T2m) in fat-infiltrated skeletal muscle, using a predetermined mono-exponential approximation to the fat decay component. This approach enabled the stable estimation of T2m in the forearm muscles of non-ambulant Duchenne muscular dystrophy patients and healthy controls from a multi-echo CPMG acquisition with only 12 echo-times. Values obtained were in good agreement with previous reports, and largely independent of muscle fat content.

Introduction

Muscle-water T2 (T2m) obtained by MRI relaxometry may be highly sensitive to pathophysiological changes in neuromuscular diseases: T2 increases due to oedema and inflammation as well as fat infiltration, and may decrease due to eventual muscle fibrosis [1]. In many neuromuscular diseases the final pathological pathway is fatty infiltration, and to determine T2m under these circumstances it is necessary to separate the water and fat signal components. One method proposed to achieve this is multi-exponential analysis of CPMG multi-echo imaging data (e.g.[1]–[6]). An appropriate model for the transverse-magnetization decay is required to fit to the experimental data; due to clinical scanning acquisition constraints, to obtain robust T2 estimates it is necessary in practice to restrict the number of free parameters in the fitting process. Recently a constrained tri-exponential model including bi-exponential fat components has been demonstrated to provide successful T2m estimates from 17-echo CPMG muscle MRI data [5], [6]. We developed a simpler signal decay model appropriate when fewer CPMG echo-time images are available, and demonstrated its use to estimate T2m from 12-echo CPMG measurements obtained in the forearms of a cohort of non-ambulant DMD boys and age-matched controls

Methods

The forearms of 9 healthy boys (mean age: 14.6y, range: 13-17y; mean BMI 21.5, range: 16.5-25.4) and 14 non-ambulant DMD patients (mean age: 13.3y, range: 10.8-17.3y; mean duration of non-ambulation 20.2 months, range: 4.7-41.6 and mean BMI 26.5, range: 20.8-41.7) were examined at 3T (Siemens Skyra) using a flexible matrix-coil with a non-fat suppressed multi-echo spin-echo sequence (TR=3000ms, 12 echo-times (TEs) from 10-120ms with 10ms interval, 9 x 6 mm slices, matrix 320x320, or 320x190, in-plane resolution 0.5625x0.5625 mm). The patients were examined at study baseline, 6 and 12 months. Ten forearm muscles were manually segmented on the central slice of a reference image set using ITK-snap software. The first TE image was excluded from the analysis to ensure a constant stimulated-echo-coherence background before a bi-exponential function [S(TE) = a.(1-ff).exp(-TE/T2m)+ a.ff.(0.75.exp(-TE/76)+0.25)], where a represents overall amplitude, and ff the fat-fraction (i.e. fat component amplitude as a proportion of the total signal amplitude) was fitted pixel-wise to the data using nonlinear least-squares minimization in a custom-written Matlab tool. The parameters of the second, mono-exponential component approximating the fat signal were determined separately as the mean values from 4 subcutaneous fat ROIs in 8 representative subjects. For quality control, values were excluded from the maps for pixel data which failed to meet empirically determined fit-quality criteria [R2 (goodness of fit) > 0.8, amplitude a lower confidence interval (CI) > a/2 , amplitude a < 10 x the 1st TE image amplitude, T2m-CI-width/T2m < 100%, fat fraction ff<50% and T2m<100ms]. T2m maps were generated and histogram metrics (mean, full width at quarter maximum, 75th percentile and skewness) for all the individual muscle ROIs in the cross-section combined were extracted; for comparison mean cross-sectional muscle ff was also determined independently by 3-point Dixon measurements. One-way ANOVA tested against the null hypothesis of equality of histogram metrics between patients longitudinally, and healthy volunteers.

Results

We obtained analysable T2 maps from all 23 subjects, with a reasonable number of pixel T2m values surviving fit quality-control criteria in most ROIs. Figure 1 shows histogram metrics (including the mean) for patients longitudinally, and for healthy controls. Patient cross-sectional mean T2m remained in the normative range while in general exhibiting wider distributions. ANOVA suggested significant difference in the full width at quarter maximum (p=0.0032) and skewness (p=0.0016) between patients and healthy controls. Figure 2 shows that cross-sectional mean T2m was largely independent of 3-point Dixon-determined ff for both groups. Figure 3 shows a representative control-subject T2m map, and Figure 4 a representative patient T2m map and corresponding independently obtained ff map.

Discussion

In initial exploratory analyses the tri-exponential model of Azzabou et al. [5] did not provide reliable parameter estimates with our data, presumably due to the different acquisition conditions in our case. The signal model used here was the most simple of those we investigated providing stable T2m estimates from our data. More sophisticated approaches are likely to be appropriate for analysing CPMG data with more TEs and higher SNR, particularly where determining a physically meaningful model of lipid compartmentalisation is a primary objective. Nevertheless, using our method we obtained forearm T2m values in both controls and DMD patients which are consistent with those previously reported [6], and appear largely independent of muscle fat content.

Acknowledgements

Study funding was by a grant from GSK to FM; NZ receives an EPSRC/GSK CASE PhD studentship.

References

[1] I. Arpan, S. C. Forbes, D. J. Lott, C. R. Senesac, M. J. Daniels, W. T. Triplett, J. K. Deol, H. L. Sweeney, G. A. Walter, and K. Vandenborne, “T2 mapping provides multiple approaches for the characterization of muscle involvement in neuromuscular diseases: a cross-sectional study of lower leg muscles in 5–15-year-old boys with Duchenne muscular dystrophy,” NMR Biomed., vol. 26, no. 3, pp. 320–328, 2013. [2] H. E. Kan, T. W. J. Scheenen, M. Wohlgemuth, D. W. J. Klomp, I. van Loosbroek-Wagenmans, G. W. Padberg, and A. Heerschap, “Quantitative MR imaging of individual muscle involvement in facioscapulohumeral muscular dystrophy,” Neuromuscul. Disord., vol. 19, no. 5, pp. 357–362, May 2009. [3] H. K. Kim, T. Laor, P. S. Horn, J. M. Racadio, B. Wong, and B. J. Dardzinski, “T2 Mapping in Duchenne Muscular Dystrophy: Distribution of Disease Activity and Correlation with Clinical Assessments1,” Radiology, vol. 255, no. 3, pp. 899 –908, Jun. 2010. [4] L. Yao and N. Gai, “Fat-Corrected T2 Measurement as a Marker of Active Muscle Disease in Inflammatory Myopathy,” Am. J. Roentgenol., vol. 198, no. 5, pp. W475–W481, May 2012. [5] N. Azzabou, P. Loureiro de Sousa, E. Caldas, and P. G. Carlier, “Validation of a generic approach to muscle water T2 determination at 3T in fat-infiltrated skeletal muscle,” J. Magn. Reson. Imaging, p. n/a–n/a, Mar. 2014. [6] C. Wary, N. Azzabou, C. Giraudeau, J. Le Louër, M. Montus, T. Voit, L. Servais, and P. Carlier, “Quantitative NMRI and NMRS identify augmented disease progression after loss of ambulation in forearms of boys with Duchenne muscular dystrophy,” NMR Biomed., vol. 28, no. 9, pp. 1150–1162, Sep. 2015.

Figures

Figure 1: Forearm muscle water T2 means and histogram metrics', ...'for patients at baseline, 6 months and 1 year, and healthy controls

Figure 2: Cross-sectional mean water T2 vs mean Dixon fat fraction for patients at baseline, 6 months and 1 year, and for healthy controls

Figure 3: A healthy volunteer T2m map

Figure 4a: A representative patient T2m map

Figure 4b: The corresponding Dixon ff map



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2280