Background

Quantitative susceptibility mapping (QSM) identifies iron accumulation and myelin loss in smoldering white matter lesions (WMLs). Yet, QSM may be also used to provide a broader understanding of focal and diffuse MS pathology.

Objectives

To study QSM features across WMLs, to assess myelin and axonal loss in WMLs with different QSM features and to quantify QSM pathology in normal-appearing white and cortical grey matter (NAWM, NAGM).

Methods

Ninety-one MS patients (62 RRMS, 29 PMS) and 72 healthy controls (HC) underwent QSM, myelin water imaging (MWI) and multishell diffusion at 3T MRI. In WMLs, cortical lesions (CLs), NAWM and NAGM, we extracted mean QSM, myelin water fraction (MWF) and neurite density index (NDI). WMLs were classified into 5 groups according to their appearance on 3D-EPI QSM: (i) isointense; (ii) with hyperintense rim, Rim+ (iii); with hypointense rim relative to the lesion core, hypo Rim; (iv) hyperintense; (v) hypointense. Mann-Whitney and Kruskal-Wallis test with Dunn’s correction for multiple comparison were used to compare (a) lesion types and (b) specific lesions vs all other WMLs. Voxel-wise comparisons of NAWM QSM were performed using Threshold-Free Cluster Enhancement (TFCE) clustering. Cortical analysis of QSM NAGM and GM-HC was performed using FreeSurfer and compared using a General Linear model (GLM).

Results

Of 1136 WMLs in QSM maps, we detected: (i) 314 (27.6%), (ii) 183 (16.1%), (iii) 16 (1.41%), (iv) 577 (50.8%) and (v) 46 (4.05%) WML. All WML exhibited lower NDI than NAWM and WM-HC (P<0.0001). Isointense lesions exhibited higher NDI (P=0.0115) and MWF (P<0.0001) than other WMLs. Rim + and hyperintense lesions exhibited lower MWF than NAWM and WM-HC (P<0.0001). Rim + lesions showed lower MWF and NDI than other WML types (P<0.001). Hypo Rim+ lesions and hypointense lesions exhibited higher MWF than other WMLs (P=0.0006, P<0.05). Hyperintense lesions exhibited lower MWF than other WMLs types (P<0.01) except Rim+ lesions. TFCE and vertex-wise cortical surface analysis showed areas throughout the NA tissue, where QSM is either lower or higher compared to healthy tissue in HC and in PMS compared to RMS (P<0.01).

Conclusions

QSM is sensitive to diffuse and focal pathology with various myelin and axonal characteristics. We hypothesize that isointense WMLs show high repair activity, hypointense WMLs are remyelinated lesions and hyperintense WMLs are chronic inactive lesions. MRI-histopathology work is ongoing to confirm these findings.

Background

Microstructural biophysical models reconstructed from advanced diffusion MRI (dMRI) data provide quantitative measures (qMs), which inform about the brain tissue microenvironment, based on different assumptions.

Objectives

To compare the sensitivity of available qMs to focal pathology in multiple sclerosis (MS), and to explore which qMs– or combinations of qMs – are best correlated with patients’ disability.

Methods

dMRI (1.8 mm isotropic resolution, 149 directions, b-values were 0, 700, 1000, 2000, 3000 s/mm²) was acquired from 67 relapsing-remitting and 33 progressive MS patients (median EDSS: 2.5). The qMs for the isotropic and intra-axonal compartments were derived from the following available models: Ball and Stick, NODDI, SMT-NODDI, MCMDI, NODDIDA, DIAMOND, Microstructure Bayesian approach (MB) and microstructure fingerprinting. In total, 13 qMs were included and subject-wise normalized within brain tissue (nqMs).

To identify the nqMs sensitive to focal pathology, an attention-based convolutional neural network (aCNN) was built to (a) classify randomly sampled WM lesion and perilesional WM patches and (b) generate attention weights (AWs) representing the relative importance of the qMs in the classification. Twenty patients were randomly selected in the test dataset (709 lesion patches and 746 perilesional WM patches), and the rest were in the cross-validation (CV) dataset (2925 lesion patches and 3176 perilesional WM patches). The performance metric was the area under the receiver operating characteristic curve (AUC). Because of the correlation between the nqMs, which may influence the relative AWs, we performed 10-fold CV and selected the nqMS that most contributed to the classification.

To assess which nqMS – or combination of nqMS was best correlated with EDSS, we used Spearman’s correlation coefficient (ρ) with two-sided 20000 permutation tests and followed by Bonferroni correction.

Results

The test AUC was 0.911 indicating the aCNN learned the right AWs to differentiate lesions and perilesional WM. The most discriminating nqMs included isotropic and intra-axonal compartments from MB, the neural density index (NDI) from the NODDI and the intra-axonal compartment from MCMDI.

The sum of isotropic and intra-axonal compartments of the MB (sMB) showed the strongest correlation with EDSS (ρ=-0.40,corr. p<0.0001) followed by the sum of sMB and NDI (ρ=-0.30,corr. p<0.05), and the sum of sMB and intra-axonal compartment from MCMDI (ρ=-0.32,corr. p<0.05). None of the selected nqMs as a single measure and their other combinations correlated with EDSS.

Conclusions

By performing aCNN-aided selection of the openly available WM quantitative measures, we have identified the measures most sensitive to MS focal pathology; furthermore, we have derived a new contrast that – by combining the measures of isotropic and intracellular diffusion – strongly correlated with patients’ disability.

Background

Myelin water imaging (MWI) and neurite orientation dispersion and density imaging (NODDI) provide sensitive surrogate markers of myelin and axonal content in lesions and normal-appearing tissue. However, to date, there is scarce information about the relationship of these measures with (i) disability; and (ii) the axonal damage specific biomarker serum neurofilament light chain (sNfL).

Objectives

To explore the correlation of MWI and NODDI measures in MS lesions and in normal-appearing (NA) brain tissue with disability and sNfL.

Methods

Ninety-one MS patients (62 relapsing-remitting MS-RRMS and 29 progressive MS-PMS) underwent MWI and NODDI. Mean myelin water fraction (MWF) and neurite density index (NDI) were extracted in white matter lesions (WMLs), cortical lesions (CLs), NA white matter (NAWM) and cortical NA gray matter (CNAGM). For sNfL, a logarithmic transformation was applied to comply with normality assumption. Correlation studies between MRI measures, sNfL and EDSS were performed using linear models, with age and gender as covariates. The models were performed for the whole sample and for patients with clinical deficits only (EDSS >1).

Results

MWF and NDI did not correlate with EDSS when the entire cohort was considered (P>0.05). However, for those patients with clinical deficits (EDSS> 1), NDI in WMLs was associated with EDSS (NDI: P<0.01, beta=-10.00; N=74). We also found that MWF and NDI in WMLs were related to sNfL (MWF: P<0.01, beta=0.13; NDI: P<0.01, beta=-3.60). Again, this correlation was stronger in patients with EDSS>1 (MWF: P<0.01, beta=0.13; NDI: P <0.01, beta=-3.60).

Conclusions

Imaging surrogate markers of myelin and axon pathology in WML – and not in CLs and NA tissues - are correlated with disability and sNfL. Interestingly, associations between those imaging markers and disability/sNFL were more evident in patients with clinical deficits as compared to those without neurological deficits.

Background

The pathogenesis of neurodegeneration in multiple sclerosis (MS) is multifactorial and the determinants of brain atrophy rates are not completely understood.

Objectives

To investigate the association between annualized atrophy rate (AAR) of multiple brain measures (regional cortical thickness (CTh), volumes of basal ganglia, thalamus, white matter, gray matter, brain and brain parenchymal fraction (BPF)) and: (1) annualized rate of new and enlarging white matter lesions (WMLs); (2) annualized rate of resolved WMLs; (3) occurrence of progression independent of relapse activity (PIRA) during follow-up.

Methods

We included 1573 1.5T or 3T brain MRI scans from 378 patients of the Swiss MS Cohort Study (331 relapsing-remitting MS (RRMS), 27 clinically isolated syndrome (CIS), 11 secondary-progressive MS (SPMS), 9 primary-progressive MS (PPMS); 70% female; median age: 41.9 yrs; disease duration: 8.3 yrs; EDSS: 2.0; follow-up time: 4.0 yrs). Longitudinal changes in WMLs were obtained using an automated prototype (LeMan-PV). Brain volumes and CTh AARs were obtained using FreeSurfer longitudinal pipeline (v6.0) after WMLs filling. In patients fulfilling PIRA an EDSS progression had to be confirmed ≥6 months after the index event. Multivariable generalized linear models were used to model the association between AAR (dependent variable) and independent variables (1-3), correcting for age, sex, disease duration and baseline EDSS. p-values were adjusted for Bonferroni multiple comparison correction; for vertex-wise CTh analysis, Monte Carlo Z simulation was performed (cluster threshold p<0.05).

Results

We found positive associations between annualized rate of new and enlarging WMLs and (i) CTh AAR of 8 extensive clusters (bilateral frontal, temporal and occipital regions and right insula, all p<0.01) and (ii) AAR of: caudate bilaterally (p=0.02), white matter volume, brain volume and BPF (p<0.001 for all).

We also found a negative association between annualized rate of resolved WMLs and CTh AAR in 3 cortical clusters (right insula, precentral area and anterior cingulate region, all p<0.05); no associations with AAR of volumes emerged.

57 patients fulfilled PIRA whereas 295 experienced no EDSS progression events: no significant differences in AAR measures were found between these two groups.

Conclusions

In a large cohort of MS patients, with a median follow-up of 4 years, local radiological inflammatory and reparative activity were associated with AAR in multiple brain regions. PIRA did not seem to be related to increased AAR in any of the regions studied.

Background

Microstructural biophysical models reconstructed from advanced diffusion MRI (dMRI) data provide quantitative measures (qMs), which inform about the brain tissue microenvironment, based on different assumptions.

Objectives

To compare the sensitivity of available qMs to focal pathology in multiple sclerosis (MS), and to explore which qMs– or combinations of qMs – are best correlated with patients’ disability.

Methods

dMRI (1.8 mm isotropic resolution, 149 directions, b-values were 0, 700, 1000, 2000, 3000 s/mm²) was acquired from 67 relapsing-remitting and 33 progressive MS patients (median EDSS: 2.5). The qMs for the isotropic and intra-axonal compartments were derived from the following available models: Ball and Stick, NODDI, SMT-NODDI, MCMDI, NODDIDA, DIAMOND, Microstructure Bayesian approach (MB) and microstructure fingerprinting. In total, 13 qMs were included and subject-wise normalized within brain tissue (nqMs).

To identify the nqMs sensitive to focal pathology, an attention-based convolutional neural network (aCNN) was built to (a) classify randomly sampled WM lesion and perilesional WM patches and (b) generate attention weights (AWs) representing the relative importance of the qMs in the classification. Twenty patients were randomly selected in the test dataset (709 lesion patches and 746 perilesional WM patches), and the rest were in the cross-validation (CV) dataset (2925 lesion patches and 3176 perilesional WM patches). The performance metric was the area under the receiver operating characteristic curve (AUC). Because of the correlation between the nqMs, which may influence the relative AWs, we performed 10-fold CV and selected the nqMS that most contributed to the classification.

To assess which nqMS – or combination of nqMS was best correlated with EDSS, we used Spearman’s correlation coefficient (ρ) with two-sided 20000 permutation tests and followed by Bonferroni correction.

Results

The test AUC was 0.911 indicating the aCNN learned the right AWs to differentiate lesions and perilesional WM. The most discriminating nqMs included isotropic and intra-axonal compartments from MB, the neural density index (NDI) from the NODDI and the intra-axonal compartment from MCMDI.

The sum of isotropic and intra-axonal compartments of the MB (sMB) showed the strongest correlation with EDSS (ρ=-0.40,corr. p<0.0001) followed by the sum of sMB and NDI (ρ=-0.30,corr. p<0.05), and the sum of sMB and intra-axonal compartment from MCMDI (ρ=-0.32,corr. p<0.05). None of the selected nqMs as a single measure and their other combinations correlated with EDSS.

Conclusions

By performing aCNN-aided selection of the openly available WM quantitative measures, we have identified the measures most sensitive to MS focal pathology; furthermore, we have derived a new contrast that – by combining the measures of isotropic and intracellular diffusion – strongly correlated with patients’ disability.