Quantitative T2 and atrophy in multiple sclerosis: A retrospective 7-year study using standard clinical brain images
Md. Nasir Uddin1, Kelly C. McPhee2, Gregg Blevins3, and Alan H. Wilman1

1Department of Biomedical Engineering, University of Alberta, Edmonton, AB, Canada, 2Department of Physics, University of Alberta, Edmonton, AB, Canada, 3Division of Neurology, Department of Medicine, University of Alberta, Edmonton, AB, Canada

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 T2 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 T2-weighted images and the corresponding R2 (1/T2) map from an MS patient are shown in Fig. 1. Over 7 years, significant T2 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 T2 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*ΔT2GP -0.48*ΔT2CD +1.83; where ΔT2GP and ΔT2CD are the differences in T2 over 7 years for globus pallidus and caudate respectively (Fig. 2). Changes in T2 in globus pallidus correlated with change in EDSS score over time. Single time point as well as 7-year changes in T2 in thalamus correlated with ventricular atrophy (r = -0.6, p < 0.05).

Discussion

We have retrospectively analysed 7-year T2 changes obtained from standard clinical exams and relate these to disease severity and brain atrophy in MS patients. Though rate of percent change in T2 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 T2 via ISEC at 1.5 T.

Conclusion

Seven years difference measurements using two-point T2 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 T2 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.

Figures

Figure 1: Proton density (a), T2-weighted images (b) and the corresponding R2 (1/T2) map (scaled 0 to 17 s-1) (c) containing deep grey matter obtained from a 31 year old MS patient at 1.5 T.

Figure 2: Predictor of disease severity in MS using multiple regression of the 7-year T2 difference measurements with MS severity scores (MSSS). The regression line equations: MSSS (after 7-years) = -0.71*ΔT2GP-0.48*ΔT2CD +1.83; where ΔT2GP and ΔT2CD are the differences in T2 over 7 years for globus pallidus and caudate respectively.

Table 1: Transverse relaxation times (T2) obtained from seven healthy controls using MESE via ISEC, two- point ISEC and two-point exponential fitting

Table 2: Transverse relaxation times (T2) in ms obtained from fourteen patients with MS. Percent changes over time are reported with effect size.



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