Structural connectivity in multiple sclerosis and simulation of disconnection
Elisabetta Pagani1, Maria Assunta Rocca1,2, Ermelinda De Meo1, Bruno Colombo2, Mariaemma Rodegher2, Giancarlo Comi2, Andrea Falini3, and Massimo Filippi1,2

1Neuroimaging Research Unit, San Raffaele Scientific Institute, Vita-Salute San Raffaele University, Milan, Italy, 2Department of Neurology, San Raffaele Scientific Institute, Vita-Salute San Raffaele University, Milan, Italy, 3Department of Neuroradiology, San Raffaele Scientific Institute, Vita-Salute San Raffaele University, Milan, Italy

Synopsis

Aim of the study was to quantify structural connectivity integrity in multiple sclerosis (MS) patients with different clinical phenotypes, to simulate a disconnection due to T2 visible lesions and to test its effect on network based measures. Diffusion tensor MRI was obtained from 239 MS patients and 131 healthy controls; connectivity matrices were produced and then artificially disconnected based on T2 visible lesion distribution. Global and nodal network metrics were calculated for both cases. Crucial nodes of the network were found to be different in strength between MS phenotypes. Disconnection simulation highlighted the role of T2 lesions in determining structural connectivity abnormalities.

Objectives

To quantify structural connectivity integrity in multiple sclerosis (MS) patients with different clinical phenotypes. To simulate a disconnection due to T2 visible lesions and to verify its effect on network based measures.

Background

MS is a disease characterized by focal and diffuse damage to the white matter (WM). T2 lesions are disseminated in space and time. This makes the idea of modelling the brain as a network of roads and treating lesions as disconnecting events very attractive. Diffusion tensor imaging (DTI) and tractography enable to define and quantify structural damage in brain networks (1); it was already applied in MS proving that there are effects already in the early phase of the disease (2). With this study we wished to control for the variability coming from the application of tractography directly to damaged brain, to simulate the effect of lesions on the network and to explore difference between MS phenotypes.

Methods

DT and dual echo turbo spin echo (TSE) MRI scans were obtained using a 3 tesla scanner from 239 MS patients (12 with a clinically isolated syndrome [CIS], 111 relapsing-remitting [RR] MS, 45 benign MS and 71 secondary progressive [SP] MS) and 131 healthy controls (HC). DT MRI scans were processed using FSL software (3) for distortions and motion correction (topup tool), tensor estimate and non-linear registration into the MNI space (fnirt tool). From the whole group of healthy controls, 80 were selected to create a reference group, with the aim of producing an atlas of WM tracts connecting cortical areas: after averaging the transformed DT components, tractography was run from all WM voxels (Diffusion Toolkit) and Trackvis software (www.trackvis.org) (4) was used to select those fibers connecting pairs of cortical areas defined by the AAL atlas (5). This allowed us to define edges and nodal degrees of a structural connectivity matrix. For 51 HCs and 239 MS patients, after the application of the WM tract atlas, the average fractional anisotropy (FA) value of each connection was calculated and saved within the corresponding element of the connectivity matrix. A similar procedure was used to produce the matrix where lesion load within each connection considered was saved. Such a lesional matrix was used to artificially produce disconnection in the connectivity matrix, so that when the lesion load was higher than 25 mm3, the FA values was set to 0 (Figures 1-2). Global and nodal network metrics were calculated (6).

Results

Compared to HC, MS patients showed significantly decreased (p<0.001) strength, assortativity, transitivity, global efficiency and increased average path length of the whole network. The nodes with the highest strength (hubs of the network) were the same for patients and controls. Among these, the postcentral, superior parietal, precuneus and cerebellar crus 1 e 2 had a significant reduced strength in SPMS patients compared to RRMS. When the disconnected matrices were considered, more nodes, including frontal superior, orbital, middle occipital areas, the thalamus and caudate nuclei were found to be decreased in strength in SPMS patients in comparison with RRMS (Figure 3) and the calcarine, middle occipital, superior parietal regions, precuneus and paracentral lobule in SPMS compared to benign MS.

Conclusions

Global measures of structural connectivity were significantly different in MS patients compared with HC, showing an extended disruption of structural connectivity integrity in MS patients, thus confirming that MS can be considered as a disconnection syndrome. The pattern of hubs distribution appeared to be preserved in MS patients, with a preserved architecture of brain structural networks and a reduced strength of connections of the nodes identified as hubs. The simulation of disconnection due to T2 lesions allowed us to identify the role of lesions in determining abnormalities in structural connectivity not only comparing MS patients with healthy controls, but also assessing the differences across MS phenotypes.

Acknowledgements

This work has been partially supported by a grant from Fondazione Italiana Sclerosi Multipla (FISM/2011/R/19) and by a grant from Italian Ministry of Health (GR-2009-1529671).

References

(1) Hagmann P, Cammoun L, Gigandet X, Gerhard S, Grant PE, Wedeen V, Meuli R, Thiran JP, Honey CJ, Sporns O. MR connectomics: Principles and challenges. J Neurosci Methods. 2010;194(1):34-45.

(2) Li Y, Jewells V, Kim M, Chen Y, Moon A, Armao D, Troiani L, Markovic-Plese S, Lin W, Shen D. Diffusion tensor imaging based network analysis detects alterations of neuroconnectivity in patients with clinically early relapsing-remitting multiple sclerosis. Hum Brain Mapp. 2013;34(12):3376-91.

(3) Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM.Advances in functional and structural MR image analysis and implementation as FSL.Neuroimage. 2004;23 Suppl 1:S208-19.

(4)Wang R, Benner T, Sorensen AG, Wedeen VJ. Diffusion toolkit: A software package for diffusion imaging data processing and tractography. Proc Intl Soc Mag Reson Med. 2007;Vol 15.

(5) Tzourio-Mazoyer N1, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, Mazoyer B, Joliot M. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage. 2002;15(1):273-89.

(6) Rubinov M, Sporns O. Complex network measures of brain connectivity: Uses and interpretations. Neuroimage. 2010;52:1059-1069.

Figures

Figure 1. The spatial distribution of T2 visible lesions for a representative MS patient is shown in red overlaid on the FA atlas. Regions of the AAL atlas, used to define nodes of the network, are represented in different co lours.

Figure 2. Original connectivity matrix (right) of a representative MS patient and after artificial removal of connections based on T2 lesion distribution (right).

Figure 3. Average nodal strength obtained from RRMS (blue) and SPMS (red) patients, using the original connectivity matrices (upper row) and the artificially disconnected ones (bottom row). Nodes where a significant difference was found between the two groups, are identified by labels.



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