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.