Synopsis
Proton density and T2-weighted images
are frequently used in clinical MS exams. These two images can be used to obtain accurate
T2 by fitting them with prior knowledge of RF pulse shapes and refocusing
flip angles in order to compensate for indirect and stimulated echo
contributions. After demonstrating
feasibility in healthy controls, we investigate 7-year changes in T2
in subcortical grey matter in 14 relapsing remitting MS patients and related
these with disease severity and brain atrophy.Introduction
Proton
density (PD) and T2-weighted
images have been part of clinical MRI protocols for multiple sclerosis (MS) for
many years. From these images, the quantitative transverse relaxation time (T2)
may be determined to serve as an absolute parameter to examine brain changes,
rather than using weighted images. However, a simple two-point exponential fit
is not appropriate and more advanced modelling is required [1]. We examine the
feasibility of tracking disease progression and state in relapsing-remitting MS
(RRMS) deep grey matter (DGM) over 7 years using retrospective standard clinical
MRI obtained at 1.5 T, from which quantitative T2, and brain volume changes
are derived.
Methods
Imaging: Fourteen RRMS patients (age 36±7years, disease duration 9±3years and EDSS 3.0±0.9 after 7 years) were identified who underwent 1.5T MRI at baseline and 7-years follow-up scans after obtaining written consent according to institutional regulations. The clinical protocol consisted of PD and T2-weighted fast spin echo (FSE) sequences, which were used to retrospectively calculate T2. Typical scan parameters included: PD-weighted FSE sequence (TR/ESP/ETL/TE/Nav = 2000ms/6.5ms/5/13ms/2, refocusing flip angles 180° or 150°) and T2-weighted FSE sequence (TR 4000ms or 3500ms, ETL 15 or 13, TE 91ms or 88ms, refocusing flip angles 180° or 150°, Nav 2). The clinical protocol also included: T1-weighted sequence (TI/TR/ESP/TE/Nav = 2500ms/444ms/6.5ms/111ms/2, flip angle 80°). For each subject, 19 slices covering the whole brain were acquired (voxel size 0.45×0.45×5 mm3).
T2 mapping: T2 maps were computed on a pixel-by-pixel basis using a PD and a T2-weighted image, and prior knowledge of flip angles [1]. This method extracts T2 from the spin response from all echo pathways providing indirect and stimulated echo compensation (ISEC). Flip angle maps, averaged over the healthy subjects, were provided to the ISEC fitting, as previously described [1]. Additionally, T2 maps were obtained from ISEC fitting of multiecho spin echo (MESE) data, and exponential fitting of two-echo data sets.
Method validation: Seven healthy controls (age 28±6years) were imaged with the same sequences as patients. Additionally, Flip angle maps were acquired using a double angle method [2] with correction for slice profile effects [1] from two spin-echo EPI images (TR/TE=7050ms/40ms, voxel size 1.87×1.87×2mm3) with excitation angles 60° and 120°. Since PD and T2-wighted images of patients had different TR and ETL, a signal correction was made to the PD image, and the correction factor was experimentally determined from volunteers.
Analysis: Bilateral ROIs were chosen in iron-rich DGM. Global and regional brain volumes and atrophy were measured from 2D T1-weighted images using FSL-SIENAX, and FSL-SIENA [3]. Statistical analyses were performed using SPSS. MS Severity Score (MSSS) were calculated from EDSS and disease duration [4].
Results
Table 1 presents the mean T
2 values in DGM from healthy controls, obtained using the MESE-ISEC method including all echoes, two-point exponential fitting, and two-point ISEC method for FSE180° and FSE150°. The distribution of flip angles in DGM is similar between the subjects. Example PD and T
2-weighted images and the corresponding R
2 (1/T
2) map from an MS patient are shown in Fig. 1. Over 7 years, significant T
2 changes of 2-4% were observed when using two-point ISEC (Table 2), but no significant results were found when using the commonly misused standard exponential fit (not shown), which cannot account for imperfect refocusing. Multiple regression of change in T
2 of two DGM structures showed high correlation with final time point MSSS (r = 0.76, p < 0.05). The equation to predict disease severity: MSSS = -0.71*ΔT
2GP -0.48*ΔT
2CD +1.83; where ΔT
2GP and ΔT
2CD are the differences in T
2 over 7 years for globus pallidus and caudate respectively (Fig. 2). Changes in T
2 in globus pallidus correlated with change in EDSS score over time. Single time point as well as 7-year changes in T
2 in thalamus correlated with ventricular atrophy (r = -0.6, p < 0.05).
Discussion
We have retrospectively analysed 7-year T
2 changes obtained from standard clinical exams and relate these to disease severity and brain atrophy in MS patients. Though rate of percent change in T
2 over time is very low compared to that of high field [5], iron related change in DGM over time and its association with disease severity can be identified using two-point T
2 via ISEC at 1.5 T.
Conclusion
Seven years difference measurements using
two-point T
2 via ISEC from standard clinical MS exams correlated to
disease severity in MS, as well as ventricular atrophy. Therefore, the long-term
retrospective two-point T
2 and atrophy measurements at standard MRI can
be used as markers to monitor disease course in DGM in MS.
Acknowledgements
This work was supported by grants from the Multiple Sclerosis Society of Canada, and Canadian Institutes of Health Research (CIHR). Salary support for KCM was provided by scholarships from Natural Sciences and Engineering Council of Canada, and Alberta Innovates Health Solutions.References
[1]
McPhee KC, Neuroimage 118, 2015; [2]
Stollberger R, MRM 35, 1996 [3] Smith SM, Neuroimage 23, 2004, [4] Roxburgh RH, Neurology 64, 2005, [5] Uddin MN, Mult Scler 2015.