Characterization of white matter lesions in multiple sclerosis using proton density and T1-relaxation measures

Although lesion dissemination in time is a defining characteristic of multiple sclerosis (MS), there is a limited understanding of lesion heterogeneity. Currently, conventional sequences such as fluid attenuated inversion recovery (FLAIR) and T1-weighted (T1W) data are used to assess MS lesions qualitatively. Estimating water content could provide a measure of local tissue rarefaction, or reduced tissue density, resulting from chronic inflammation. Our goal was to utilize the proton spin density (PD), derived from a rapid, multi-contrast STAGE (strategically acquired gradient echo) protocol to characterize white matter (WM) lesions seen on T2W, FLAIR and T1W data. Twenty (20) subjects with relapsing-remitting MS were scanned at 3 T using T1W, T2-weighted, FLAIR and strategically acquired gradient echo (STAGE) sequences. PD and T1 maps were derived from the STAGE data. Disease severity scores, including Extended Disability Status Scale (EDSS) and Multiple Sclerosis Functional Composite (MSFC), were correlated with total, high PD and high T1 lesion volumes. A probability map of high PD regions and all lesions across all subjects was generated. Five perilesional normal appearing WM (NAWM) bands surrounding the lesions were generated to compare the median PD and T1 values in each band with the lesional values and the global WM. T1W intensity was negatively correlated with PD as expected (R = −0.87, p < 0.01, R2 = 0.756) and the FLAIR signal was suppressed for high PD volumes within the lesions, roughly for PD ≥ 0.85. The threshold for high PD and T1 regions was set to 0.905 and 1950.3 ms, respectively. High PD regions showed a high probability of occurrence near the boundary of the lateral ventricles. EDSS score and nine-hole peg test (dominant and non-dominant hand) were significantly correlated with the total lesion volume and the volumes of high PD and T1 regions (p < 0.05). There was a significant difference in PD/T1 values between the high PD/T1 regions within the lesions and the remaining lesional tissue (p < 0.001). In addition, the PD values of first NAWM perilesional band directly adjacent to the lesional boundary displayed a significant difference (p < 0.05) compared to the global WM. Lesions with high PD and T1s had the highest probability of occurrence at the boundary of the lateral ventricles and likely represent chronic lesions with significant local tissue rarefaction. Moreover, the perilesional NAWM exhibited subtly increasing PD and T1 values from the NAWM up to the lesion boundary. Unlike on the T1 maps, the perilesional band adjacent to the lesion boundary possessed significantly higher PD value than the global WM PD values. This shows that PD maps were sensitive to the subtle changes in NAWM surrounding the lesions.