Abstract
Background:
In multiple sclerosis (MS), the location of focal lesions does not always correlate with clinical symptoms, suggesting disconnection as a major pathophysiological mechanism. Resting-state (RS) functional magnetic resonance imaging (fMRI) is believed to reflect brain functional connectivity (FC) within specific neuronal networks.
Objective:
RS-fMRI was used to investigate changes in FC within two critical networks for the understanding of MS disabilities, namely, the sensory-motor network (SMN) and the default-mode network (DMN), respectively, implicated in sensory-motor and cognitive functions.
Methods:
Thirty-four relapsing–remitting (RR), 14 secondary progressive (SP) MS patients and 25 healthy controls underwent MRI at 3T, including conventional images, T1-weighted volumes, and RS-fMRI sequences. Independent component analysis (ICA) was employed to extract maps of the relevant RS networks for every participant. Group analyses were performed to assess changes in FC within the SMN and DMN in the two MS phenotypes.
Results:
Increased FC was found in both networks of MS patients. Interestingly, specific changes in either direction were observed also between RR and SP MS groups.
Conclusions:
FC changes seem to parallel patients’ clinical state and capability of compensating for the severity of clinical/cognitive disabilities.
Keywords
Introduction
Multiple sclerosis (MS) is an inflammatory/demyelinating disease of the central nervous system (CNS) that, in a proportion of cases, results in sensory-motor and cognitive disabilities. These clinical manifestations may have a strong impact on patients’ everyday life activities, affecting both their own and their relatives’ quality of life. Both macro- and microscopic white matter (WM) and grey matter (GM) abnormalities are typically observed in the brain of MS patients, 1 and the presence of such structural damage is likely to slow down brain communication within cortical and sub-cortical networks.2,3 Conventional magnetic resonance imaging (MRI) is one of the most informative tools for MS diagnosis and quantification of macroscopic damage in the CNS. However, the correlation between measures of macroscopic lesion burden and atrophy, and measures of clinical disabilities, especially those related to cognition, still remain rather poor. 4 Resting-state functional MRI (RS-fMRI) is a relatively new technique that allows us to measure, in vivo, functional connectivity (FC) as expressed by synchronization of neural activity across different brain regions, the so-called resting-state networks (RSNs).5–8 One of the most widely investigated RSNs is the so-called ‘default-mode network’ (DMN), connecting the precuneus and the posterior cingulate cortex (PCC) to the medial prefrontal cortex (medPFC) and inferior parietal regions.9–13 This network is believed to be implicated in cognitive functioning, 14 and its disruption has been consistently demonstrated in neurodegenerative conditions associated with cognitive decline.15,16 Other consistently identified RSNs include the sensory-motor network (SMN), the core (or salience) network, the visual and auditory networks, and the right- and left-lateralized parietal-frontal networks.6,10,17–19 Within the last few years, several studies have investigated changes in FC in a number of neurological and psychiatric conditions such as dementia,15,16 MS20–23 and schizophrenia. 24 Previous RS-fMRI studies in MS20–23 have either shown patterns of decreased20,22,23 or increased FC 21 within patients’ brains. The former have been interpreted as correlates of neuronal disruption, while the latter have been read as compensatory mechanisms of brain plasticity, typically detectable at early disease stages. One possible interpretation of the conflicting results reported in the literature is that the initial compensatory mechanism might be progressively outweighed by accumulating neuronal damage, and that changes in the balance between these two factors might contribute to the evolution from relapsing–remitting MS (RRMS) to secondary progressive MS (SPMS).
Here, recruiting a cohort of patients with either RRMS or SPMS, we aimed to assess whether specific changes in functional brain connectivity may account for patients clinical manifestations, with a focus on motor and cognitive disabilities. We hypothesised that clinical differences observed in MS evolution may be reflected by specific modifications of FC within networks sub-serving sensory-motor (SMN) and cognitive (DMN) functions.
Methods
Participants
For this study we recruited 48 patients with clinically definite MS (male/female (M/F) ratio=14/34; mean (standard deviation (SD)) age=38.5 (9.9) years, mean (SD) years of formal education=12.6 (3.0)) according to McDonald’s criteria. 25 Thirty-four patients had a RR phenotype (M/F ratio=8/26; mean (SD) age=35.70 (9.1) years), while the remaining 14 patients were in a SP disease course (M/F ratio=6/8; mean (SD) age=44.6 (9.1) years). Healthy volunteers had a higher education level (years=16.25, SD=2.4) than MS patients (years=12.58, SD=3.0, t(53)= −4.8, p=0.000, t-test for independent samples), while no difference in years of teaching was observed between the two MS sub-groups. Before MRI, all patients underwent an extensive clinical and neuropsychological examination. Twenty-two healthy subjects were also enrolled into the study and served as controls (M/F ratio=8/14; mean (SD) age=35.8 (11.1) years). To reduce any potential source of variability due to hemispheric dominance, all participants (patients and controls) had to be right-handed. 26
Inclusion criteria for patients were: absence of other medical conditions, no corticosteroid therapy for at least one month before entering the study, and age ranging from 20–65 years. Beyond Expanded Disability Status Scale (EDSS) score assessment, the Timed 25-Foot Walk (T25-FW; trial 1 and 2), a quantitative mobility and leg function performance test, 27 was administered to all patients. Ethical approval from the Ethics Committee of the Santa Lucia Foundation and written informed consent were obtained before study initiation.
Neuropsychological assessment
Two trained neuropsychologists explored the following cognitive domains: (a) information processing speed, through the administration of the Symbol Digit Modalities Test (SDMT) 28 ; (b) Working memory and cognitive flexibility, using the Paced Auditory Serial Addition Test at 3 s (PASAT-3) 29 ; (c) Praxis and visuo-spatial memory abilities, assessed through the Rey Complex Figure Test–immediate (RCFC) and the 20-minute delayed recall version (RCFC-Recall). 30 For the PASAT-3, we computed z-scores for the hits and errors: z-scores were considered normal for values >0 and pathological for values< –1. All scores were adjusted for sex, age and education, according to the Italian normative data. 31 Clinical and neuropsychological data are summarised in Table 1.
Principal characteristics of multiple sclerosis (MS) patients.
EDSS: Expanded Disability Status Scale; PASAT: Paced Auditory Serial Addition Test; RCFC: Rey Complex Figure Test- immediate recall; RCFC-recall: Rey Complex Figure Test-20 min delayed recall; SD: standard deviation; SDMT: Symbol Digit Modalities Test; T25-FW: Timed 25-Foot Walk.
Mann-Whitney non parametric test, for independent samples. aLevel of significance refers to Chi Square non-parametric test.
MRI scanning protocol
All imaging was obtained, in a single session, using a head-only 3.0T MRI scanner (Siemens Magnetom Allegra, Siemens Medical Solutions, Erlangen, Germany). The MRI acquisition protocol included: (a) a dual-echo turbo spin echo (TSE) (repetition time (TR)=6190 ms, echo times (TE) =12/109 ms, echo train length (ETL)=5; matrix=256×192; FOV=230×172.5 mm2; 48 contiguous 3 mm thick slices); (b) a fast fluid attenuated inversion recovery (FLAIR) (TR=8170 ms, TE=96 ms, inversion time (TI)=2100 ms; ETL=13; with same FOV, matrix and number of slices as TSE) to use as a reference for lesion identification; (c) a T2*-weighted echo-planar imaging (EPI) sensitised to BOLD contrast (TR=2080 ms, TE=30 ms, 32 axial slices parallel to AC-PC line, matrix=64×64, pixel size=3×3 mm2, slice thickness=2.5 mm, flip angle=70°) for RS-fMRI (total number of volumes=220); (d) a magnetisation prepared rapid gradient echo (MPRAGE) sequence (TR=2500 ms; TE=2.74 ms; inversion time=900 ms; flip angle=8°; matrix=256 ×208×176; slab thickness =1 mm; FoV=256×208×176 mm3). During this acquisition, subjects were instructed to keep their eyes closed, not to think of anything in particular, and not to fall asleep.
Lesion volume assessment
In each patient, T2-hyperintense lesions were identified by consensus of two independent observers on the short echo images of the TSE. Lesions were outlined on the same scan using a semi-automated local thresholding contouring technique (Jim 4.0, Xinapse System, Leicester, UK, http://www.xinapse.com/). FLAIR and long echo TSE scans were always used as a reference to increase the confidence in lesion identification. The total lesion volume (T2LL) was calculated for each patient.
RS-fMRI analysis
MRI data preprocessing was performed using statistical parametric mapping (SPM8, Wellcome Department of Imaging Neuroscience), and in-house software implemented in Matlab (MathWorks, Inc.).
MPRAGE scans were segmented to extract GM, WM and CSF probability maps for each subject. The total GM volume was computed by integrating the GM map intensity values over the whole brain. To account for subjects’ head size differences, GM-volume was expressed as percent of the total intracranial volume and yielded the GM fraction. 33
For each subject the first four volumes of the fMRI series were discarded to allow for T1 equilibration effects. Images were realigned, corrected for slice-time, normalised into Montreal Neurological Institute (MNI) space, and smoothed with a 8 mm3 Gaussian kernel. Finally, all images were filtered by a phase-insensitive band-pass filter (pass band 0.01–0.08 Hz) to reduce the effect of low frequency drift and high frequency physiological noise. A model-free analysis was employed by using independent component analysis (ICA) implemented in the Group ICA for fMRI toolbox (GIFT), in order to allow for a simultaneous separation into individual components. ICA was performed on all participants’ data grouped together. Briefly, ICA first concatenates the individual data across time, and then produces a computation of subject specific components and time courses. The toolbox performed the analysis in three steps: (a) data reduction; (b) application of the FastICA algorithm; and (c) back-reconstruction for each individual subject. The DMN and the SMN were visually identified. Full factorial analyses were performed for the SMN and DMN components. For between-group comparisons, each contrast was masked by the main effect of groups (p cluster-level uncorrected=0.005) and set at p cluster level corrected<0.05. In the patient groups only, voxel-wise correlations between FC and clinical variables were also investigated. For the SMN, the association with EDSS and T25-FW were explored, while for the DMN, we looked at the association with SDMT and PASAT-3 scores. In both cases, the interaction between the covariate of interest and the phenotype (RRMS vs SPMS) was modelled. Age and GM-fraction were always entered as covariates of no interest normalized by intra-cranial-volume 32. The threshold for statistical significance was set at p values cluster level uncorrected<0.005.
Results
Demographic, clinical and neuropsychological characteristics
Both healthy volunteers (t-Student Test, t(26)=3.9, p<0.007) and RRMS patients (t-Student Test, t(29)= −1.9, p<0.05) were younger than SPMS patients. RRMS patients had lower EDSS scores than those with SPMS (median (range) EDSS scores=1.5 (1.0–4.5) and 6.0 (3.0–7.5), respectively; Chi-square=31.7 df=13, p=0.003) and performed better at the T25-FW in both trial 1 and 2 (see Table 1, for statistical significance). Significant difference in the average T2LL was also observed between the two MS groups, showing higher T2LL in SPMS (8.5 (SD=6.8) ml), as compared to RRMS patients (6.22 (SD=8.5) ml; U(47)=152, Z= −1.95, p=0.05, Mann-Whitney Test).
The majority of our patients performed poorly on the PASAT-3, reporting scores outside the range of normality 33 with 30/34 RRMS (88%) and 14/14 SPMS (100%) performing above the 5th percentile with respect to the number of errors. Conversely, only 5.8% (for RRMS) and 7.1% (for SPMS) exhibited low scores on SDMT (below the 5th percentile), and 5.8% (for RRMS) and 14.2% (for SPMS) reported a low score at the visuo-spatial memory test (RCFC-Recall) (below the 5th percentile).
When considering all patients together, correlation analyses showed positive associations between level of motor impairment (as measured by EDSS and T25-FW trials 1 and 2), disease duration and T2LL (EDSS score vs disease duration, r=0.30; EDSS score vs T25-FW trial 1, 1 r=0.40; EDSS score vs T25-FW 2, r=0.37; T2LL vs disease duration, r=0.50; statistical threshold p<0.01, Pearson correlation).
RS-fMRI results
Among the 20 components estimated by ICA, 10 RSNs already reported by others6,10,17–19 were identified, namely, the DMN, the SMN, the core/salience, two visual components (occipital and medial), the anterior cingulate cortex/ medial frontal cortex, the temporal, the auditory, and both left and right fronto-parietal networks. The networks were identified by computing the correlation with in-house RSN templates 34 and confirmed by visual inspection.
According to our hypotheses, we selected the SMN and the DMN for FC analyses. The SMN includes the primary (Brodmann area (BA)1, 2, 3) and secondary somatosensory (BA5) areas, the primary (BA4) and the supplementary (BA6) motor cortex, and the supramarginal gyrus (SMG, BA39/40) bilaterally.6,17 The DMN includes the inferior parietal lobules (IPL, BA 39/40), the PCC (BA30, 23, 31) and precuneus (BA7), areas of the inferior, medial and superior frontal gyri (FG, BA8, 9, 10, 47) and the lateral temporal cortex (BA21) (see Figure 1 for identified components within each single group of participants).

Sensory-motor network (SMN) and default-mode network (DMN) single-group maps derived from independent component analysis (ICA) analysis in sagittal and coronal view in healthy control (HC) subjects, relapsing–remitting multiple sclerosis (RRMS) and secondary progressive multiple sclerosis (SPMS) patients.
Between-group differences of FC within the SMN
When investigating the SMN, significant increases of FC were observed in both groups of MS patients, compared against healthy controls (HCs), in the most anterior part of this network, including the pre/postcentral gyrus, extending outward to the anterior cingulate cortex (ACC, BA32). The inverse contrast (HCs>all MS patients together) did not reveal any significant difference in FC. Further, increased FC was found in the precentral gyrus (BA4) and in the ACC (BA32), in RR, as compared to SPMS patients (see Table 2 and Figure 2).
Between-group differences in functional connectivity within the sensory-motor network (SMN).
BA: Brodmann area; HCs: healthy controls; LH: left hemisphere; MS: multiple sclerosis; RH: right hemisphere; RR: relapsing–remitting; SP: secondary progressive.
T-contrasts are masked with the main effect of groups set at p=0.005. Statistical threshold: p-cluster-level uncorrected=0.05. Inverse comparisons showed no significant differences.

(a) Within the sensory-motor network (SMN), all multiple sclerosis (MS) patients, compared against healthy controls (HCs), showed increased functional brain connectivity (regions in blue) within the pre- and post-central gyrus (Brodmann area (BA) 4/6 and 3), the supramarginal gyrus (SMG, BA40), and the anterior cingulate cortex (ACC, BA 32). Also, increased functional connectivity (FC) was observed in near areas when contrasting relapsing–remitting (RR) against secondary progressive multiple sclerosis (SPMS) patients (regions in violet). The opposite contrast did not reveal any significant difference. Plot shows the BOLD signal changes in the anterior cingulate (p-corrected<0.001). FG: frontal gyri.
Between-group differences of FC within the DMN
Within the DMN, higher FC was observed in all MS patients, compared to HCs, in the right SMG (BA 40). No significant difference was found when considering the opposite contrast. When directly comparing the two MS subgroups, SPMS patients showed higher FC in the right SMG (BA 40) and in the left PCC (BA24), extending to the precuneus (BA7) (see Figure 2(b), and Table 3). An area of higher FC was also observed in the left inferior/middle temporal gyrus (BA21) of SPMS patients.
Between-group differences in functional connectivity within the default mode network (DMN).
BA: Brodmann area; HCs: healthy controls; LH: left hemisphere; MS: multiple sclerosis; RH: right hemisphere; RR: relapsing–remitting; SP: secondary progressive.
T-contrasts are masked with the main effect of groups set at p=0.005. Statistical threshold: p-cluster-level uncorrected=0.05. Inverse comparisons showed no significant differences.
Correlation analyses
Within the SMN, no significant correlations were found between FC and motor disability across the whole group. However, a significant effect of interaction was found between phenotype (RR or SPMS) and patients’ level of motor disability (as measured by T25-FW1 score) on FC in the right pre-central gyrus (i.e. inverse correlation was present in SP, but not in RRMS patients) (see Figure 3). Within the DMN, a significant positive association was found between patients’ number of mistakes at PASAT-3 (the only cognitive test in which patients scored abnormally) and FC in the ACC (p<0.005). This positive association was stronger for SP, than RRMS patients (see Figure 4).

Sensory-motor network (SMN). Interaction effect between functional connectivity (FC) in the right pre-central gyrus (Brodmann area (BA) 4/6) and the Timed 25-Foot Walk (T25-FW) trial 1 scores across the two patients’ phenotypes. The design matrix used to test the interaction is reported on the left. It includes the two groups of patients modelled separately (relapsing–remitting (RR) versus secondary progressive multiple sclerosis (SPMS)), their correspondent T25-FW trial 1 scores (as covariate of interest), and age and grey matter (GM) fraction (as covariates of no interest). In the right plot, association between individual T25-FW 1 scores and FC values derived from the brain area of interaction for the two groups of patients, the RRMS (violet circles) and the SPMS (green rhombi), is shown. This direct correlation is present in the SP, but not in the RRMS group, and it remains significant also when removing the RRMS outlier.

Default-mode network (DMN). Correlation between functional connectivity (FC) in the paracingulate and the number of errors in the Paced Auditory Serial Addition Test at 3 s (PASAT-3) across the two patients’ phenotypes. The design matrix used to test the interaction is reported on the left. It includes the two groups of patients modelled separately (relapsing–remitting (RR) versus secondary progressive multiple sclerosis (SPMS)), their number of mistakes in the PASAT-3 (as covariate of interest), and age and grey matter (GM) fraction (as covariates of no interest). In the right plot, the association between number of errors in the PASAT-3 and FC values derived from the brain area of main effect for the two groups of patients, the RR (violet circles) and the SPMS patients (green rhombi), is represented.
Discussion
In this RS-fMRI study, we investigated how MS pathology modulates brain FC across disease evolution. We focused our analyses on the SMN and the DMN, considering FC in these two networks as associated with patients’ motor and cognitive impairment. Moreover, we included patients with either RR or SPMS, considering RRMS frequently evolves to SPMS. Therefore, within the limitation of a cross-sectional design, our sample variability makes the study indirectly informative on the pathophysiological evolution of MS. As expected, SPMS compared to RRMS patients were older, had a longer disease duration, higher motor and cognitive disabilities, and higher T2-LL. When compared to HCs, MS patients showed areas of increased FC in both, the SMN and DMN. In the SMN, FC was increased in the right sensory-motor and premotor cortex and in the anterior cingulate gyrus. From the direct comparison between RR and SPMS patients, it appeared that such a positive FC modulation was driven by the former group. This suggests that, at earlier disease stages, when the motor disability is minimal, RRMS brains tend to compensate for the occurrence of tissue damage by increasing FC in the non-dominant hemisphere. As has been suggested in previous work in patients with clinically isolated syndrome (CIS), 21 compensatory mechanisms of brain plasticity may contribute in contrasting the permanent clinical effects of MS. Furthermore, when RRMS proceeds towards the SPMS phenotype the compensatory effects tend to decrease. These results are also consistent with findings showing that the re-allocation of neuronal resources might help to preserve the effects of early structural damage in other neurological conditions.15,16 Further, different patterns of FC in both networks were identified when comparing the two MS phenotypes. Within the SMN, RRMS patients showed higher FC, than those with SPMS, in frontal regions. We suggest that this represents a compensation mechanism for decline that is still not remarkable in RRMS but, instead, becomes more significant at later disease stages. 21 Indeed, a positive association between increases in FC in the pre-central gyrus and the T25-FW was observed in RRMS, but not in the SPMS patients. When considering the DMN, SPMS compared to RRMS patients showed increased FC within posterior-parietal regions, in the PCC/precuneus and in the middle and inferior temporal gyri. Among the RSNs, the DMN has by far received the most attention throughout the clinical and research community. This set of regions is typically more intensely activated during the rest and relatively de-activated during demanding tasks requiring focused attention, such as working memory tasks and visuo-spatial tasks. Again, we suggest that this hyper-activation might rely on a neural compensation mechanism which counterweights structural brain damage occurring in later MS stages. Conversely, this effect is not yet evident in RRMS, at least in our group of patients whose cognitive impairment was minimal. Furthermore, SPMS patients also revealed a more remarkable pattern of FC within the DMN throughout some temporal regions, falling out of the selected component. These areas, including the middle and inferior temporal gyrus, are involved in a number of cognitive processes, such as language and visual perception. 35
Neuropsychological assessment showed that our MS patients scored abnormally only in the PASAT-3 test, overall showing 91% of pathological scores. Correlation analysis between indices of FC and test scores showed a positive association between FC in the ACC and the PASAT-3, especially in the SPMS group. This might suggest that changes in FC (controlled here by brain atrophy) within specific networks may detect specific patterns of pathophysiological events (e.g. inflammation and neurodegeneration) which characterise relevant stages of MS evolution.
Previous studies using RS-fMRI within MS samples showed quite controversial results. Some authors observed reduced FC in the SMN 20 and in the DMN of patients with differing MS phenotypes.21,22 Another study highlighted increased FC in patients with CIS, 23 with no significant changes found in patients with RRMS. These studies show little evidence of a coherent neural organisation. However, it is possible that various pathophysiological mechanisms may differently contribute to brain modification observed in MS phenotypes and stages.
Our study has several limitations. First, we did not collect neuropsychological data in our healthy control group. Though, we suggest that normative data might somehow be representative of a useful reference point to compare our patients with. Second, the SPMS patient group was smaller than the other two groups. In addition, we focused on the DMN and SMN only. While this choice is justified by the relevance of sensory-motor and cognitive impairment in MS, we cannot exclude the presence of changes in FC in other networks and, because of the diffuse nature of the disease, this would indeed seem likely. Future work should address this question.
In conclusion, this RS-fMRI exploratory study investigated for the first time RRMS and SPMS patients together. FC changes seem to parallel patients’ capability of compensating for the severity of clinical/cognitive disabilities. Despite the limitation of a cross-sectional design, this study suggests that initial increases in SMN and DMN FC reduce when MS symptoms worsen, or a different balance between pathophysiological mechanisms induces a change of MS phenotype. Future longitudinal studies are needed however to confirm our interpretation of functional brain reorganisation across MS evolution by recruiting larger populations of patients and by following them up over time.
Footnotes
Conflicts of interest
The authors declare that there are no conflicts of interest.
Funding
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.
