Networks of microstructural damage predict disability in multiple sclerosis

Abstract

BackgroundNetwork-based measures are emerging MRI markers in multiple sclerosis (MS). We aimed to identify networks of white (WM) and grey matter (GM) damage that predict disability progression and cognitive worsening using data-driven methods.MethodsWe analysed data from 1836 participants with different MS phenotypes (843 in a discovery cohort and 842 in a replication cohort). We calculated standardised T1-weighted/T2-weighted (sT1w/T2w) ratio maps in brain GM and WM, and applied spatial independent component analysis to identify networks of covarying microstructural damage. Clinical outcomes were Expanded Disability Status Scale worsening confirmed at 24 weeks (24-week confirmed disability progression (CDP)) and time to cognitive worsening assessed by the Symbol Digit Modalities Test (SDMT). We used Cox proportional hazard models to calculate predictive value of network measures.ResultsWe identified 8 WM and 7 GM sT1w/T2w networks (of regional covariation in sT1w/T2w measures) in both cohorts. Network loading represents the degree of covariation in regional T1/T2 ratio within a given network. The loading factor in the anterior corona radiata and temporo-parieto-frontal components were associated with higher risks of developing CDP both in the discovery (HR=0.85, p<0.05 and HR=0.83, p<0.05, respectively) and replication cohorts (HR=0.84, p<0.05 and HR=0.80, p<0.005, respectively). The decreasing or increasing loading factor in the arcuate fasciculus, corpus callosum, deep GM, cortico-cerebellar patterns and lesion load were associated with a higher risk of developing SDMT worsening both in the discovery (HR=0.82, p<0.01; HR=0.87, p<0.05; HR=0.75, p<0.001; HR=0.86, p<0.05 and HR=1.27, p<0.0001) and replication cohorts (HR=0.82, p<0.005; HR=0.73, p<0.0001; HR=0.80, p<0.005; HR=0.85, p<0.01 and HR=1.26, p<0.0001).ConclusionsGM and WM networks of microstructural changes predict disability and cognitive worsening in MS. Our approach may be used to identify patients at greater risk of disability worsening and stratify cohorts in treatment trials.