In vivo gradients of thalamic damage in paediatric multiple sclerosis: a window into pathology

The thalamus represents one of the first structures affected by neurodegenerative processes in multiple sclerosis. A greater thalamic volume reduction over time, on its CSF side, has been described in paediatric multiple sclerosis patients. However, its determinants and the underlying pathological changes, likely occurring before this phenomenon becomes measurable, have never been explored. Using a multiparametric magnetic resonance approach, we quantified, in vivo, the different processes that can involve the thalamus in terms of focal lesions, microstructural damage and atrophy in paediatric multiple sclerosis patients and their distribution according to the distance from CSF/thalamus interface and thalamus/white matter interface. In 70 paediatric multiple sclerosis patients and 26 age- and sex-matched healthy controls, we tested for differences in thalamic volume and quantitative MRI metrics-including fractional anisotropy, mean diffusivity and T1/T2-weighted ratio-in the whole thalamus and in thalamic white matter, globally and within concentric bands originating from CSF/thalamus interface. In paediatric multiple sclerosis patients, the relationship of thalamic abnormalities with cortical thickness and white matter lesions was also investigated. Compared to healthy controls, patients had significantly increased fractional anisotropy in whole thalamus (f2 = 0.145; P = 0.03), reduced fractional anisotropy (f2 = 0.219; P = 0.006) and increased mean diffusivity (f2 = 0.178; P = 0.009) in thalamic white matter and a trend towards a reduced thalamic volume (f2 = 0.027; P = 0.058). By segmenting the whole thalamus and thalamic white matter into concentric bands, in paediatric multiple sclerosis we detected significant fractional anisotropy abnormalities in bands nearest to CSF (f2 = 0.208; P = 0.002) and in those closest to white matter (f2 range = 0.183-0.369; P range = 0.010-0.046), while we found significant mean diffusivity (f2 range = 0.101-0.369; P range = 0.018-0.042) and T1/T2-weighted ratio (f2 = 0.773; P = 0.001) abnormalities in thalamic bands closest to CSF. The increase in fractional anisotropy and decrease in mean diffusivity detected at the CSF/thalamus interface correlated with cortical thickness reduction (r range = -0.27-0.34; P range = 0.004-0.028), whereas the increase in fractional anisotropy detected at the thalamus/white matter interface correlated with white matter lesion volumes (r range = 0.24-0.27; P range = 0.006-0.050). Globally, our results support the hypothesis of heterogeneous pathological processes, including retrograde degeneration from white matter lesions and CSF-mediated damage, leading to thalamic microstructural abnormalities, likely preceding macroscopic tissue loss. Assessing thalamic microstructural changes using a multiparametric magnetic resonance approach may represent a target to monitor the efficacy of neuroprotective strategies early in the disease course.