Skip to main content

Early remission in multiple sclerosis is linked to altered coherence of the Cerebellar Network



The development of permanent disability in multiple sclerosis (MS) is highly variable among patients, and the exact mechanisms that contribute to this disability remain unknown.


Following the idea that the brain has intrinsic network organization, we investigated changes of functional networks in MS patients to identify possible links between network reorganization and remission from clinical episodes in MS. Eighteen relapsing–remitting MS patients (RRMS) in their first clinical manifestation underwent resting-state functional MRI and again during remission. We used ten template networks, identified from independent component analysis, to compare changes in network coherence for each patient compared to those of 44 healthy controls from the Human Connectome Project test–retest dataset (two-sample t-test of pre-post differences). Combining a binomial test with Monte Carlo procedures, we tested four models of how functional coherence might change between the first clinical episode and remission: a network can change its coherence (a) with itself (“one-with-self”), (b) with another network (“one-with-other”), or (c) with a set of other networks (“one-with-many”), or (d) multiple networks can change their coherence with respect to one common network (“many-with-one”).


We found evidence supporting two of these hypotheses: coherence decreased between the Executive Control Network and several other networks (“one-with-many” hypothesis), and a set of networks altered their coherence with the Cerebellar Network (“many-with-one” hypothesis).


Given the unexpected commonality of the Cerebellar Network’s altered coherence with other networks (a finding present in more than 70% of the patients, despite their clinical heterogeneity), we conclude that remission in MS may result from learning processes mediated by the Cerebellar Network.


In multiple sclerosis (MS), the most commonly observed disease course is characterized by an alternating pattern of new or enhanced neurological symptoms (“relapses”) and phases of partial or complete recovery (“remission”). These clinical episodes often co-occur with an acute inflammatory insult, causing demyelination and neurodegeneration to the central nervous system (CNS). Such sudden injury is believed to be the source of the symptoms, as demyelination hampers signal transmission [1]. Consequently, given the network nature of the brain, which relies on complex systemwide communication to carry out tasks [2], it is reasonable to assume that any interruption of signal transmission results in suboptimal task performance, which can present as neurological symptoms [3]. Although during remission the inflammation has subsided, sclerotic, demyelinated plaques remain [4]; nevertheless, patients are often completely symptom-free, especially following early relapses [5]. Therefore, it would appear as though the brain employs mechanisms beyond lesion repair in order to regain and maintain function following lesion occurrence.

Classic work has shown that the brain is able to flexibly reorganize (generally referred to as “neuroplasticity” [6, 7]), even in elderly people [8]. Studies using functional magnetic resonance imaging (fMRI) have demonstrated that the neural correlates underlying specific tasks differ between MS patients and healthy controls (HCs), thereby suggesting neuroplastic reorganization. For example, Pantano and colleagues [9] investigated the neural correlates of hand movements in MS patients with an ipsi- as well as a contralateral hemiparesis. This study showed that, even in the absence of any clinical manifestation, the same task implicates different cortical areas in MS patients compared to HCs. Additionally, functional reorganization in MS patients has been observed at clinical recovery after pseudotumoral lesions [10]. Related observations have been reported by other similar studies [11,12,13].

One drawback of such task-based fMRI is that its findings can generally only be interpreted within the context of specific tasks that probe specific functions. However, specific tasks that uncover large-scale changes in network organization between a heterogeneous group of MS-patients and a group of HCs are not known or may not exist. Therefore, examining network properties speaks in favor of employing resting-state fMRI (rs-fMRI), which allows the analysis of network-level activity in the brain in an unconstrained manner [14]. Many investigations have looked into the reorganization of resting-state networks in MS patients, indeed demonstrating some general, fundamental differences in this organization compared to HCs (e.g. [15, 16], for a review see [17]). However, such studies do not address the mechanisms that underlie successful or failed remission.

In this study, we recruited MS patients suffering from their first clinical manifestation. Our approach sought to investigate the changes in functional coherence (a type of functional connectivity, [14]) of brain networks between the first episode, during which patients exhibited neurological symptoms, and remission (when they were symptom-free, on average four weeks after the first scan). With this experimental design, we investigated the following questions: although MS is a highly heterogeneous disease, is there nevertheless a common brain network that changes in all MS patients between a clinical episode and remission? If so, does this network change with respect to itself (1) a one-with-self change), with respect to a specific other network (2) a one-with-other change), or with respect to a variety of other networks (3) a one-with-many change)? Furthermore, is there a particular region of the brain (or set of regions) to which patients’ networks systematically change their coherence (4) a many-with-one change)? These hypotheses (Fig. 1) follow the idea of a putative general mechanism driving the recovery toward remission [18] but implicate different physiological levels with respect to the function of remission in MS.

Fig. 1
figure 1

Hypotheses. Conceptual representation of the hypotheses depicted as A brain maps and B adjacency matrices. We hypothesize that functional network coherence changes between relapse and remission can occur as changes from an original network to itself (one-with-self), to another network (one-with-other), or to a set of networks (one-with-many), or that regardless of the original network, these networks change with respect to a common network (many-with-one)




For this study, we recruited 18 relapsing–remitting MS patients (14 females, 4 males, age (mean + /– standard deviation) = 28.33 + /– 9.36 years) in their initial clinical episode (diagnosed according to the 2017 revisions of the McDonald criteria [19]: dissemination in space was demonstrated by MRI, dissemination in time by presence of oligoclonal bands). Each patient was scanned twice: once during the first episode [i.e., with ongoing presentation of symptoms (EDSS ≥ 1 in all cases)] and the second time during remission (when we expected a decrease in EDSS scores, on average four weeks after the first scan). The patients showed different clinical symptoms (Table 1). For six patients, no structured EDSS assessment was available in remission. We calculated the median EDSS scores during the first clinical episode and remission and used a Wilcoxon signed-rank test to test the null hypothesis that median ranks from both time points are equal. We calculated T2 lesion load for all but one patient (for whom no T2 FLAIR image was available from the first session due to excessive motion) in both sessions using the Lesion Prediction Algorithm of the Lesion Segmentation Toolbox [20] and assessed changes of lesion load between the two scans with a paired t-test. All patients provided written informed consent to participate in this study. All study procedures followed safety guidelines for MRI research at the University of Regensburg, complied with the 1964 Helsinki Declaration and its later amendments, and were approved by the ethics committee of the University of Regensburg.

Table 1 Demographic and clinical data of the patients

Healthy control participants

We used data from the Human Connectome Project 1200 Subjects Data Release ( as a control dataset [21, 22]. Specifically, we made use of the minimally-preprocessed rs-fMRI data, which includes complete test–retest data from 44 participants (30 females, age range = 22–35 years). Participants were scanned on two separate days, separated by an interscan interval between 3 and 11 months. All study procedures of the HCP protocol were approved by the Institutional Review Board at Washington University in St. Louis.

Neuroimaging data acquisition

Data acquisition of the 18 MS patients was carried out on a 3 Tesla Prisma scanner (Siemens, Erlangen, Germany). Resting-state functional images were acquired with a T2*-weighted EPI sequence (instruction to participants: eyes closed, repetition time (TR) = 2000 ms, echo time (TE) = 30.0 ms, flip angle (FA) = 70°, field of view (FOV) = 192 × 192 mm2, matrix size = 64 × 64, voxel resolution (VR) = 3 × 3 × 3.6 mm3, 33 axial slices (acquired ascending-interleaved), phase encoding direction = anterior–posterior, echo spacing = 0.58 ms, pixel bandwidth (BW) = 2605 Hz/Px, scan duration = 22 min 6 s, 660 volumes per run; three dummy volumes were acquired at the beginning of the scan to account for T1-saturation). For co-registration of the functional data to high-resolution native anatomical space, we acquired a T1-weighted scan using a 3D Turboflash MPRAGE sequence (TR = 1910 ms, TE = 3.67 ms, inversion time (TI) = 1040 ms, 180 sagittal slices, FA = 9°, FOV = 250 × 250 mm2, matrix size = 256 × 256, VR = 0.98 × 0.98 × 1 mm3, BW = 180 Hz/Px, scan duration = 4 min 25 s). Additionally, to quantify lesion load, we acquired a T2-weighted 2D Fluid Attenuated Inversion Recovery (FLAIR) image (TR = 9000 ms, TE = 108 ms, TI = 2500 ms, 27 axial slices, FA = 150°, FOV = 186 × 230  mm2, matrix size = 260 × 320, VR = 0.72 × 0.72 × 5.5 mm3, BW = 290 Hz/Px, scan duration = 2 min 14 s) and calculated the lesion load using an automatic lesion-segmentation method [20].

The HCP data was acquired on a 3 Tesla Connectome Skyra Scanner (instruction to participants: eyes closed, TR = 720 ms, TE = 33.1 ms, FA = 52°, FOV = 208 × 180mm2, matrix size = 104 × 90, voxel resolution = 2 × 2 × 2 mm3, 72 slices, phase encoding direction = right–left, multiband factor = 8, echo spacing = 0.58 ms, BW = 2290 Hz/Px, scan duration = 14 min 33 s, 1200 volumes per run). The full imaging protocols can be found online at

Although the scanning parameters used to acquire the patient dataset differed from those used to acquire the HCP dataset, the concern of a methodological bias is mitigated by the fact that we did not compare groups directly; rather, we compared within-group longitudinal changes between the two groups [i.e., differences (groups) of differences (time points)]. As such, given that data at both time points for a given member of each group were acquired with the same parameters, one would only expect the presence of biases if one set of scanning parameters was inherently more sensitive towards detecting connectivity-based differences over time. In this case, our method, comprising fewer acquired volumes, was less sensitive, whereas the HCP method, comprising more acquired volumes, was more sensitive [23]. Thus, the scanning parameters applied to the patient sample were more statistically conservative and therefore less likely to yield false positives.

Neuroimaging data analysis

Analysis of the acquired neuroimaging data was carried out with the FMRIB Software Library (FSL, version 6.0) [24, 25], the Connectome Workbench (version 1.4.2) [26, 27] and the CoSMoMVPA toolbox [28] for MATLAB R2018b (The Mathworks, Natick, USA). Additionally, the HCP data had already been pre-processed using tools from FreeSurfer, FSL and the Connectome Workbench at the time we used it; see [22] for details on this procedure.


Pre-processing steps included brain extraction [29], slice-time correction, and motion correction with respect to the middle volume of each run (using trilinear interpolation with 6 degrees of freedom (DOFs) [30]). We did not apply spatial smoothing, as its implementation has recently been shown to artificially increase the similarity of networks between subjects in functional connectivity analyses [31]. Additionally, we did not apply temporal filtering, which would have interfered with subsequent automatic removal of motion artifacts (see below).

Each patient’s functional scan was then linearly co-registered to the native high-resolution anatomical scan using a rigid body (6 DOFs) transformation. For statistical comparisons, we also performed a second, non-linear, co-registration of each functional scan to MNI152 2 mm standard space. Because non-linear registrations are only reasonable within the same modalities (e.g., structural, functional, etc.), we first implemented a 12 DOF affine transformation of the high-resolution anatomical scan to MNI152 2 mm space, which we then used as a starting point for the non-linear registration of the functional data to the standard space [32].

In the final pre-processing step, we corrected for head motion-induced artifacts using FSL’s Automated Removal Of Motion Artifacts (AROMA) tool, using its non-aggressive strategy [33]. This algorithm consists of three steps: first, independent component analysis (ICA) is applied to decompose the functional data into a set of spatially independent components. Second, those components that most likely relate to head motion are automatically identified based on four criteria (high-frequency content, maximum correlation with realignment parameters, edge fractions, and CSF fractions). Finally, the functional data is denoised by removing the variance explained by only the motion-related components. We performed this analysis on the functional data in native space, which we later transformed into MNI standard space using the previously generated transformation matrices.

The maps resulting from this procedure were in volumetric space, but a 2D representation of cortical data has several advantages over a 3D-representation, including enhanced alignment to the geometry of the cortex as well as increased statistical power [34]. Therefore, we transformed the preprocessed data into CIFTI files (Connectivity Informatics Technology Initiative, using the Ciftify toolbox [35], a file format which renders cortical data on the surface while maintaining volumetric information of subcortical and cerebellar data. For simplicity, we will refer to the smallest units of our data as “voxels”, comprising both surface-based and cerebellar/subcortical representations.

As a control data set, we used minimally preprocessed data (both functional and structural scans) from the HCP-retest dataset, which are also provided in CIFTI format. This data has been preprocessed following the HCP functional preprocessing pipeline [22]. The main preprocessing steps were identical between our pipeline and the HCP pipeline, encompassing removal of spatial/distortion artifacts, surface generation, registration and alignment to standard space.

Assessing changes in brain network organization between relapse and remission

The central goal of our analysis (illustrated in Figs. 2 and 3) was to investigate whether the organization of patients’ brain networks changed between relapse and remission, compared to that of a sample of HC. Using dual regression [36] we generated subject-specific component maps based on ten previously published spatial components (from ICA, Table 2) and associated time courses [14] for both patients and HCs (Fig. 2A). Here, we used a version of this template in CIFTI file format [37]. Then, each of the ten component templates was regressed (as spatial regressors in a multiple regression) into each participant’s (i.e., each patient and each of the HCs) 4D space–time dataset. This procedure resulted in a set of ten subject-specific time series, one per spatial component. The resulting time series were used as temporal predictors in a multiple regression of each participant’s 4D dataset, resulting in ten subject-specific spatial maps of parameter estimates (PE). We then smoothed the PE maps with a Gaussian kernel of 5 mm FWHM (in order to improve the spatial correspondence between the data from patients and HCs). We applied this procedure to the scan during the first clinical episode (or baseline in HC) and the scan at remission (or follow-up in HCs). As we were interested in changes in brain networks, we computed the difference between the PE maps (for each of the ten templates separately) for each patient/HC from the two scans (i.e., we calculated remission minus first clinical episode for the patients and follow-up minus baseline for the HC, Fig. 2B, upper row). We then converted the patients’ difference maps, voxelwise, to z-score maps (Fig. 2B, lower row), using the mean and standard deviation from the difference maps (follow-up—baseline) of the 44 HCs. We computed these maps for each of the ten ICA-networks separately.

Fig. 2
figure 2

Calculating statistical maps of network changes. A First, we applied dual regression for a set of ten ICA-network templates to each control’s and each patient’s data sets, both for first clinical episode (baseline) and remission (follow-up). Next, (B, upper row) we computed the difference map for each patient between the two time points. Finally, (B, lower row) we converted the patients’ difference maps to z-maps, by calculating, for each voxel, the z-score with respect to the control group

Fig. 3
figure 3

Statistical assessment. A Using a Monte-Carlo procedure, we generated null distributions to determine the expected degree of change in network organization between follow-up and baseline. B From these null distributions, we calculated empirical p-values for each voxel of the patients’ z-maps, which yielded regions reflecting a change in network organization

Table 2 Overview of the analyzed networks, according to Smith et al. [14]

In order to carry out statistical inference on the difference maps, we used a permutation testing scheme [38] based on the null hypothesis of no difference between patients and controls. The procedure was as follows (Fig. 3A): for each network, of all 62 analyzed subjects (combining 18 patients and 44 controls), we randomly selected the difference maps from 18 subjects (reflecting the number of patients in our sample) and calculated 18 z-scored difference maps, as described in the previous paragraph, from the distribution of the remaining 44 subjects (which could contain patients). From each of the resulting z-scored difference maps, we stored both the maximum and the minimum z-scores [to control for the inflated family-wise error rate (FWER)] in order to generate two separate null distributions for testing both increases and decreases of coherence. We repeated this process 10,000 times and subsequently compared the z-scored difference maps of the 18 patients (i.e. ten maps for each patient, reflecting the ten analyzed networks) to the permutation-based null distributions by counting how many values of the null distributions were more extreme than the observed z-scores in each voxel of the patients’ maps. Dividing this count by the number of permutations yielded FWER-corrected p-values (pFWER) (Fig. 3B). We considered voxels with a pFWER lower than 0.025 (i.e., 0.05, two-tailed) to indicate changes in network organization that would not be predicted by chance. Hereinafter, we refer to networks in which we observed changes between the two time points as “altered networks”. For each altered network, if we identified one or more clusters of voxels with pFWER < 0.025, we extracted those clusters using Workbench’s cifti-find-clusters algorithm. In the following sections, we refer to such clusters of voxels as regions of interest (ROIs).

Constructing adjacency matrices between altered networks and original networks

For all of the ROIs of a given altered network of a single patient, we identified to which of the ten original networks the ROIs belong, according to the ICA z-scores from the ten spatial components. This procedure yielded a 10 × 10 adjacency matrix of altered networks (rows) and original networks (columns). Thus, we counted the number of ROIs present for a given altered network and placed this count data in the corresponding column of the matrix. As we were interested in the number of patients that showed similar coherence changes, these data were binarized, so that, even if a patient had multiple ROIs reflecting coherence changes between the same two networks, this information appeared in the matrix as a one. This 10 × 10 binary adjacency matrix was populated for each patient and then summed across patients, yielding a group-level frequency matrix, which allowed us to perform group-level statistics to test our hypotheses. Such matrices were not necessarily symmetric because a pre- post change in an altered network that changed its connectivity with an original network does not imply that the original network also showed a pre- post change.

Assessing the one-with-one, one-with-many, and many-with-one hypotheses

Our hypotheses asked whether specific patterns appeared in the adjacency matrices, which would indicate systematicity between original networks and altered networks. Namely, we were interested in a one-with-one hypothesis (denoted by a single square in the adjacency matrix, Fig. 1B), a one-with-many hypothesis (denoted by a row in the adjacency matrix), and a many-with-one hypothesis (denoted by a column in the adjacency matrix). Thus, we needed to determine whether the frequencies for a single square or the pattern of frequencies for a given column or row were unlikely to occur by chance, which we carried out using a binomial test and Monte Carlo procedures.

To this end, we first calculated the expected number of patients that would contribute to a given cell of the matrix, which we obtained by summing all 100 frequencies of the group-level matrix (separately for coherence increases and decreases) and then dividing this sum by the number of elements of the matrix. This procedure yielded 1.99 patients for the increases and 2.58 patients for the decreases. Using these expected counts, we then ran a binomial test for each cell of the matrices to determine the probability of the observed number of patients, given the total number of patients that participated in the experiment. These p-values were then converted to z-scores (using a mean of 0 and standard deviation of 1) for ease of visually presenting heatmaps. The z-scores were ultimately thresholded at a value of 1.96 (i.e., p < 0.025) to be able to determine whether any column-, row, or cell-effects were present by counting how many elements in each column and row survived the threshold.

Assessing potential one-with-one effects was carried out by seeking lone cells with a p-value lower than the Bonferroni-corrected threshold (to account for the inflated family-wise error rate [FWER] from performing 100 tests in two matrices) of p < 0.025/200 (i.e., z > 3.6623). Assessing potential one-with-many or many-with-one effects (i.e., row- and column-effects, respectively) was carried out by first counting how many elements of a given row or column survived a threshold of p < 0.025 (i.e., z > 1.96) and then running a Monte Carlo procedure to determine how many elements of a row or column would be predicted by chance. This randomization procedure consisted of shuffling each patient’s binary matrix, recomputing the group-level frequency matrix, recomputing the binomial test, thresholding the resulting random p-values at p < 0.025, summing the number of elements per row and column, and storing the maximum row and column sums (also to control for the inflated FWER from ultimately testing for multiple column and row effects). We repeated this procedure 100,000 times yielding two distributions of maximum expected elements per column (once for coherence increases, once for coherence decreases) and two distributions of maximum expected elements per row (again, for coherence increases and decreases separately). We used these maximum sums to determine whether the observed number of elements comprising a column- or row-effect in the original data (following thresholding) were not likely to be predicted by chance.

Follow-up individual-level analysis

Given that our results pointed to unexpected commonalities across patients involving the Cerebellar Network and Executive Control Network, we opted to carry out a follow-up analysis to determine whether such effects were robustly present in all patients rather than merely being driven by a minority of the sample. To this end, as each ROI could belong to a different network (see above), for each patient we counted how often each “original network” appeared in the set of ROIs, ultimately yielding an individual-level probability distribution of the ten networks. These probabilities were gathered separately for coherence increases and decreases and then arcsine transformed before being submitted to a two-factor repeated-measures analysis of variance (ANOVA) with factors Coherence and Network. This procedure was repeated but instead counting how often each “altered network” appeared in order to also obtain individual-level probability distributions for the altered networks.


Our primary goal was to investigate whether there were regions in the brain that changed their coherence to functional brain networks (obtained via dual-regression ICA) between the first clinical episode (pre) and remission (post) in MS patients. We addressed this question by investigating ten networks (Table 2) for each patient, analyzing increases and decreases of functional coherence separately. To determine when pre-post changes in coherence differed from the expected degree of change present in a group of healthy controls, we employed Monte Carlo procedures, which revealed networks that differed between the first clinical episode and remission. We used these patterns of network alterations to test hypotheses regarding whether there exists systematicity across patients with respect to the networks that change and the networks to which they change.

Clinical scores, but not lesion load, change between first clinical episode and remission

The median EDSS score of our patient sample was 2.25 (interquartile range = 1.0) at first clinical episode and 1.0 (interquartile range = 0.50) at remission. The Wilcoxon signed-rank test confirmed the decrease in EDSS scores over time (z = 53.50, p = 0.0059). On the other hand, we were unable to find differences in T2 lesion load (t16 = 0.4364, p = 0.6684) between the two time points.

Many-with-one and one-with-many: coherence changes for the cerebellar and executive control networks

The Monte Carlo procedures we used to test the one-with-many (row-effects) and many-with-one (column-effects) hypotheses revealed evidence for both hypotheses. Evidence for the many-with-one hypothesis manifested with respect to the Cerebellar Network (Crb) in terms of both increased (six of ten networks present, pFWER = 9.9 × 10−6) and decreased coherence (seven of ten networks present, pFWER = 9.9 × 10−6) (Fig. 4). For the coherence increases, the networks contributing to the column-effect were the Occipital Pole (Occ, z = 2.8597, p = 0.002), Lateral Visual Areas (LVA, z = 3.906, p = 4.69 × 10−5), Sensorimotor Network (SMN, z = 3.906, p = 4.69 × 10−5), Auditory Network (Aud, z = 2.3105, p = 0.01), Executive Control Network (ECN, z = 4.4116, p = 5.13 × 10−6), and right Frontoparietal Network (rFP, z = 2.3105, p = 0.01). For the coherence decreases, the networks contributing to the column-effect were Medial Visual Areas (MVA, z = 2.356, p = 0.009), Occ (z = 3.3809, p = 3.6 × 10−4), LVA (z = 2.8749, p = 0.002), Default Mode Network (DMN, z = 2.356, p = 0.009), Crb (z = 2.356, p = 0.009) (SMN (, z = 3.8778, p = 5.3 × 10−5), and ECN (, z = 3.3809, p = 3.6 × 10−4) (Fig. 4).

Fig. 4
figure 4

Main results. A Relative frequencies depicting how many patients had a particular altered network that changed its coherence with a particular original network. B Thresholded z-scores (at z > 1.96) from inverting p-values obtained via a binomial test. Following all statistical procedures, the results provide support for the many-with-one hypothesis, as several altered networks across patients appeared to both increase and decrease their coherence with the Cerebellar Network. There was additional evidence for the one-with-many hypothesis, in that coherence increased between the ECN and a set of other networks. Column and row effects that surpassed the respective FWER-corrected threshold of p < 0.05 (two-tailed) are highlighted by the red contours

Evidence for the one-with-many hypothesis consisted of the ECN increasing its coherence with four of ten networks (pFWER =  = 2.9 × 10−4): the LVA (z = 2.31, p = 0.01), Crb (z = 4.4116, p = 5.13 × 10−6), SMN (2.31, p = 0.01), and lFP (z = 2.8597, p = 0.002) (Fig. 4). The z-scores and p-values presented in this section for the individual networks came from the binomial test.

Moreover, the binomial test did not reveal evidence for the “one-with-other” hypothesis between any networks that were not already contributing to a “one-with-many” or “many-with-one” effect.

Robust coherence changes of the cerebellar network in vast majority of patients

Our individual-level analysis revealed that the unexpected commonality of the Cerebellar Network altering its functional coherence with other networks was not driven by only a few patients. Rather, the Cerebellar Network was the most-likely (or equally-most-likely) “original network” [mean probability =  ~ 39 ± 6% (SEM)] for coherence increases (Fig. 5, left), occurring in 13 patients (i.e., ~ 72% of the sample; p = 4.31 × 10−11, from a binomial test B(18, 0.1056); Additional file 1: Figure S1). A fourteenth patient also showed coherence increases involving the Cerebellar Network, though it was not the most likely original network. With respect to the coherence decreases (Fig. 5, right), the Cerebellar Network was the most-likely original network (mean probability =  ~ 39 ± 5% [SEM]), occurring in 14 patients (i.e., ~ 78% of the sample; p = 5.82 × 10−12, from a binomial test B(18, 0.1167); Additional file 1: Figure S2). Two additional patients showed coherence decreases involving the Cerebellar Network, though it was not the most likely original network. The degree of coherence changes in the cerebellar network did not correlate with the changes in the EDSS scores (see Additional file 1). These observations were supported by the results from the ANOVA, which revealed differences between networks (F9,153 = 27.1, p = 1.49 × 10–27) and coherence directions (F1,17 = 6.15, p = 0.0239), but not an interaction between the two factors (F9,153 = 1.28, p = 0.2528). Moreover, post-hoc achieved power analyses using G*Power [39] revealed statistical power of 99.9985% for having found the effect of coherence increases involving the Cerebellar Network in 13 of 18 patients (needing only to be present in four patients, given an expected probability of 0.1056). Similarly, achieved statistical power was 99.9991% for having found the effect of coherence decreases involving the Cerebellar Network in 14 of 18 patients (needing only to be present in five patients, with an expected probability of 0.1167).

Fig. 5
figure 5

Group-averaged probabilities: original networks. The probability that each of the 18 patients’ ROIs, which showed coherence increases (red) and decreases (blue) with the altered network(s), belonged to a given one of the ten component networks. The Cerebellar Network was the most probable original network, appearing with a probability of ~ 39% for coherence increases (present in 14/18 patients) and ~ 39% for coherence decreases (present in 16/18 patients). For both panels, translucent circles reflect individual results, while opaque circles reflect the group -average. Error bars depict standard error of the mean

With respect to the individual-level probabilities of the original networks, the ANOVA revealed differences between networks (F9,153 = 3.13, p = 0.0017) and between coherence directions (F1,17 = 5.04, p = 0.038) but not an interaction between the two factors (F9,153 = 0.84, p = 0.58) (Fig. 6, Additional file 1: Figures 3, 4).

Fig. 6
figure 6

Group-averaged probabilities: altered networks. The probability that each of the 18 patients’ networks showed coherence increases (red) or decreases (blue), with respect to the sample of healthy controls. For both the coherence increases and decreases, the Sensorimotor and Executive Control Networks were most likely to be altered, appearing with a probability of ~ 20%. For both panels, translucent circles reflect individual results, while opaque circles reflect the group -average. Error bars depict standard error of the mean


“One-with-many effect”: executive control network

Our analysis revealed evidence for the “one-with-many” hypothesis, in that the functional coherence of the ECN increased with respect to the LVA, Crb, SMN, and lFP. This finding stands in partial contrast to previous work, which has shown, for example, that decreases in functional connectivity between the ECN and the DMN at one time point (as compared to controls) predict increased clinical disability (i.e., higher EDSS scores) [40]. Our findings suggest instead that a relative increase of the ECN with other networks is compatible with remission (i.e., lower EDSS scores). One important distinction, however, is that our sample consisted of patients suffering from their initial clinical episode who were scanned during that first episode and again in remission, thereby allowing us to measure changes over time within patients. In contrast, such prior work investigated relapsing–remitting (RRMS) patients that were scanned only once; as such, the authors were only able to describe differences with respect to healthy controls but not which changes may have occurred within a given patient. It is therefore possible that the connectivity behavior of the ECN differs at different stages of MS.

“Many-with-one-effect”: cerebellar network

Most critically, we found evidence for the “many-with-one” hypothesis, in that functional coherence changed (both increasing and decreasing) between a set of networks and the Cerebellar Network (Fig. 4), and that this result is rather robust at the individual level (i.e., present in more than 70% of patients), despite the clinical heterogeneity of the sample. To demonstrate the robustness of this effect, we carried out a split-half analysis in which we compared the results from the first nine patients with the results from the second nine patients (see Additional file 1: Fig. S1).

Such a finding raises the question of whether there is some capacity of the Cerebellar Network that allows RRMS patients to overcome their symptoms; this idea is supported by converging evidence that has tied functional and structural aspects of the cerebellum to the severity of MS disease progression. For example, MS patients with early cerebellar dysfunction tend to develop disability more quickly [41]. Similarly, Kutzelnigg and colleagues showed that the extent of cerebellar cortical demyelination is associated with overall MS pathology [42]. Recent cross-sectional studies investigating functional connectivity of MS patients demonstrated a relationship between symptom severity and alterations of the cerebellum’s functional connectivity [43,44,45]. Our findings extend such studies by providing additional evidence of functional changes with respect to the Cerebellar Network when measured at relapse and remission, and, more importantly, yield a new perspective that functional connectivity is altered between the Cerebellar Network and a variety of other networks, irrespective of clinical symptoms. These findings give rise to the idea that the Cerebellar Network may play a guiding role in recovery from diverse MS symptoms in the early disease course and may protect against permanent disability.

Despite these findings, the question of how the Cerebellar Network may alleviate MS symptoms remains open. One mechanism by which the Cerebellar Network might be involved in remission is that of learning. Building off of the original theory of cerebellar learning [46], the current Albus-Marr-Ito theory of cerebellar learning posits that information from sensory consequences of actions is directly conveyed to the cerebellum via the inferior olive [47, 48]. Additionally, copies of any efferent motor program are propagated to the cerebellum via collaterals from cortical neurons [49, 50]. Thus, as the cerebellum can access both efferent plans and afferent consequences, it has the potential to register when there is a deviation between the two and correct if necessary [51]. Though Marr originally thought that the cerebellum was involved only in learning motor programs, recent work has shown that the cerebellum receives input from distributed cortical areas and is involved in a broad range of cognitive functions beyond movement [52].

With respect to MS and the idea of a protective role of the cerebellum, earlier work provides evidence that such motor and/or cognitive learning may underlie the improvement of symptoms in MS patients. For example, using fMRI Saini and colleagues observed activation patterns in clinically stable MS patients that are reportedly found in healthy controls during motor learning [53]. Similarly, MS patients, whose performance on cognitive tasks successfully improved following cognitive training programs (as compared to MS patients that did not receive cognitive training), showed increased activation within the cerebellum [54] and altered functional connectivity of the cerebellum [55]. Thus, a link between learning, altered cerebellar recruitment, and MS symptom decline may explain the effectiveness of cognitive/motor training programs in MS rehabilitation [56, 57].

Limitations and future outlook

One limitation of our study concerns the interpretation of the “many-with-one” effect that we observed implicating the Cerebellar Network in processes underlying remission in MS. By definition, a “many-with-one” effect is composed of multiple “one-with-one” effects that surpassed a particular threshold and, together, form the pattern of a column. Because these columnar patterns were statistically unlikely to occur by chance, we interpreted them as evidence for the “many-with-one” hypothesis, rather than a set of isolated effects, each favoring the “one-with-one” hypothesis. Regardless, we cannot rule out the possibility that treating these effects in isolation would be the correct interpretation. For example, consider a scenario in which different patient subtypes have different altered networks that change their coherence to the Cerebellar Network. By analyzing these patients as a single group, one would observe what looks like a “many-with-one” effect, which is merely the result of the heterogeneous sample. However, our sample size is too small to differentiate subtypes. As such, this idea should be scrutinized in larger studies that can test the effectiveness of methods such as ours in identifying patient subtypes, thereby potentially improving personalized treatment strategies [58].

Despite this interpretational shortcoming of the “many-with-one” effect, there remain the facts (1) that functional coherence was commonly altered between a variety of networks and the Cerebellar Network and (2) that the Cerebellar Network was either the most likely original network or tied for most likely original network (in at least 70% of the patients for both the coherence increases and decreases). Therefore, this unexpected commonality across patients still speaks in favor of our general interpretation, providing evidence that the Cerebellar Network’s specific role in MS should be further investigated, for example, as a possible target for symptom-alleviating therapies.

Additionally, the open question regarding the meaning of increases and decreases of functional coherence of the Cerebellar Network remains unknown. Such differences in functional connectivity are believed to reflect different underlying mechanisms [59]; with respect to our findings, it is possible that the coherence increases and decreases reflect different types of plasticity at the synaptic level, specifically long-term potentiation (LTP) and long-term depression (LTD). Both of these learning mechanisms have been observed for cerebellar neurons [60,61,62,63], and LTP has been associated with symptom decline in MS patients [64, 65]. Moreover, prior work has shown that LTP increases functional connectivity in the rat cortex [66, 67]; it is thus conceivable that LTD analogously decreases functional connectivity, but this relationship has yet to be experimentally demonstrated. Consequently, alterations of cerebellar functional connectivity in MS patients may underlie neuroplastic changes that could theoretically predict symptom severity.


In summary, we investigated changes of functional connectivity in MS patients between their first clinical episode and remission using ICA-based rs-fMRI. Our analyses revealed that, the Cerebellar Network and a heterogeneous set of other networks across different MS patient tended to alter their functional coherence to one another, and that there is approximately a 39% chance of observing functional coherence changes in the Cerebellar Network for a given patient following the initial clinical episode. This finding suggests that the Cerebellar Network may be functionally compensating for structural damage to distributed brain regions, which could help explain the oft-observed mismatch between structural damage and clinical disability [68]. Furthermore, these alterations in functional coherence manifested as both increases and decreases, which may reflect neuroplastic mechanisms of learning. Within this framework, subsequent studies could investigate the cerebellum’s capacity for such learning mechanisms as a potential biomarker that predicts an individual’s clinical progression; more generally, future work could seek to understand the specific role of the cerebellum in the presence of distributed damage to the structural and functional architecture of the brain, especially at different stages of MS.

Availability of data and materials

Statistical data generated during the analyses presented in this study will be made available to readers upon reasonable request to the corresponding author, but raw/preprocessed data acquired from patients in this study cannot be shared due to the European Union’s General Data Protection Regulation.



Auditory Network


Analysis of variance


Automated Removal of Motion Artifacts




Connectivity Informatics Technology Initiative


Cerebrospinal fluid


Cerebellar Network


Default Mode Network


Degree of freedom


Executive Control Network


Expanded Disability Status Scale


Echo-planar imaging


Flip angle


Functional connectivity


Fluid-attenuated inversion recovery


Field of view


Family-wise error rate


Healthy control(s)


Human Connectome Project


Independent component analysis


Frontoparietal Network, left


Lateral Visual Areas


Long-term depression


Long-term potentiation


Medial Visual Areas


Montreal Neurological Space


Magnetization-prepared rapid acquisition gradient echo


Multiple sclerosis


Occipital Pole


Parameter estimates


Frontoparietal Network, right


Region of interest


Resting-state functional MRI


Standard error of the mean


Sensorimotor Network


Echo time


Time to repeat


Voxel resolution


  1. Lassmann H. Multiple sclerosis pathology. Cold Spring Harb Perspect Med. 2018;8:a028936.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Sporns O. Structure and function of complex brain networks. Dialogues Clin Neurosci. 2013;15:247–62.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Stam CJ. Modern network science of neurological disorders. Nat Rev Neurosci. 2014;15:683–95.

    Article  CAS  PubMed  Google Scholar 

  4. Traboulsee A. MRI relapses have significant pathologic and clinical implications in multiple sclerosis. J Neurol Sci. 2007;256:19–22.

    Article  Google Scholar 

  5. Mowry EM, Pesic M, Grimes B, Deen S, Bacchetti P, Waubant E. Demyelinating events in early multiple sclerosis have inherent severity and recovery. Neurology. 2009;72:602–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Maguire EA, Gadian DG, Johnsrude IS, Good CD, Ashburner J, Frackowiak RS, Frith CD. Navigation-related structural change in the hippocampi of taxi drivers. Proc Natl Acad Sci USA. 2000;97:4398–403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Robertson IH, Murre JMJ. Rehabilitation of brain damage: brain plasticity and principles of guided recovery. Psychol Bull. 1999;125:544–75.

    Article  CAS  PubMed  Google Scholar 

  8. Heuninckx S, Wenderoth N, Swinnen SP. Systems neuroplasticity in the aging brain: recruiting additional neural resources for successful motor performance in elderly persons. J Neurosci. 2008;28:91–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Pantano P, Iannetti GD, Caramia F, Mainero C, Di Legge S, Bozzao L, Pozzilli C, Lenzi GL. Cortical motor reorganization after a single clinical attack of multiple sclerosis. Brain. 2002;125:1607–15.

    Article  PubMed  Google Scholar 

  10. Mezzapesa DM, Rocca MA, Rodegher M, Comi G, Filippi M. Functional cortical changes of the sensorimotor network are associated with clinical recovery in multiple sclerosis. Hum Brain Mapp. 2008;29:562–73.

    Article  PubMed  Google Scholar 

  11. Braun C, Schweizer R, Elbert T, Birbaumer N, Taub E. Differential activation in somatosensory cortex for different discrimination tasks. J Neurosci. 2000;20:446–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Reddy H, Narayanan S, Matthews PM, Hoge RD, Pike GB, Duquette P, Antel J, Arnold DL. Relating axonal injury to functional recovery in MS. Neurology. 2000;54:236–236.

    Article  CAS  PubMed  Google Scholar 

  13. Reddy H, Narayanan S, Woolrich M, Mitsumori T, Lapierre Y, Arnold DL, Matthews PM. Functional brain reorganization for hand movement in patients with multiple sclerosis: defining distinct effects of injury and disability. Brain. 2002;125:2646–57.

    Article  CAS  PubMed  Google Scholar 

  14. Smith SM, Fox PT, Miller KL, Glahn DC, Fox PM, Mackay CE, Filippini N, Watkins KE, Toro R, Laird AR, Beckmann CF. Correspondence of the brain’s functional architecture during activation and rest. Proc Natl Acad Sci USA. 2009;106:13040–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Rocca MA, Valsasina P, Leavitt VM, Rodegher M, Radaelli M, Riccitelli GC, Martinelli V, Martinelli-Boneschi F, Falini A, Comi G, Filippi M. Functional network connectivity abnormalities in multiple sclerosis: Correlations with disability and cognitive impairment. Mult Scler J. 2018;24:459–71.

    Article  Google Scholar 

  16. Sbardella E, Tona F, Petsas N, Upadhyay N, Piattella M, Filippini N, Prosperini L, Pozzilli C, Pantano P. Functional connectivity changes and their relationship with clinical disability and white matter integrity in patients with relapsing-remitting multiple sclerosis. Mult Scler J. 2015;21:1681–93.

    Article  CAS  Google Scholar 

  17. Tahedl M, Levine SM, Greenlee MW, Weissert R, Schwarzbach JV. Functional connectivity in multiple sclerosis: recent findings and future directions. Front Neurol. 2018;9:1–18.

    Article  Google Scholar 

  18. Franzmeier N, Duering M, Weiner M, Dichgans M, Ewers M. Left frontal cortex connectivity underlies cognitive reserve in prodromal Alzheimer disease. Neurology. 2017;88:1054–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Thompson AJ, Banwell BL, Barkhof F, Carroll WM, Coetzee T, Comi G, Correale J, Fazekas F, Filippi M, Freedman MS, Fujihara K, Galetta SL, Hartung HP, Kappos L, Lublin FD, Marrie RA, Miller AE, Miller DH, Montalban X, Mowry EM, Sorensen PS, Tintoré M, Traboulsee AL, Trojano M, Uitdehaag BMJ, Vukusic S, Waubant E, Weinshenker BG, Reingold SC, Cohen JA. Diagnosis of multiple sclerosis: 2017 revisions of the McDonald criteria. Lancet Neurol. 2018;17:162–73.

    Article  PubMed  Google Scholar 

  20. Schmidt P, Gaser C, Arsic M, Buck D, Förschler A, Berthele A, Hoshi M, Ilg R, Schmid VJ, Zimmer C, Hemmer B, Mühlau M. An automated tool for detection of FLAIR-hyperintense white-matter lesions in multiple sclerosis. Neuroimage. 2012;59:3774–83.

    Article  PubMed  Google Scholar 

  21. Van Essen DC, Smith SM, Barch DM, Behrens TEJ, Yacoub E, Ugurbil K, WU-Minn HCP Consortium. The WU-Minn Human Connectome Project: an overview. Neuroimage. 2013;80:62–79.

    Article  PubMed  Google Scholar 

  22. Glasser MF, Sotiropoulos SN, Wilson JA, Coalson TS, Fischl B, Andersson JL, Xu J, Jbabdi S, Webster M, Polimeni JR, Van EDC, Jenkinson M, Hcp W. The minimal preprocessing pipelines for the Human Connectome Project. Neuroimage. 2013;80:105–24.

    Article  PubMed  Google Scholar 

  23. Airan RD, Vogelstein JT, Pillai JJ, Caffo B, Pekar JJ, Sair HI. Factors affecting characterization and localization of interindividual differences in functional connectivity using MRI. Hum Brain Mapp. 2016;37:1986–97.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Jenkinson M, Beckmann CF, Behrens TEJ, Woolrich MW, Smith SM. FSL. Neuroimage. 2012;62:782–90.

    Article  PubMed  Google Scholar 

  25. Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TEJ, 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.

    Article  PubMed  Google Scholar 

  26. Marcus DS, Harms MP, Snyder AZ, Jenkinson M, Wilson JA, Glasser MF, Barch DM, Archie KA, Burgess GC, Ramaratnam M, Hodge M, Horton W, Herrick R, Olsen T, McKay M, House M, Hileman M, Reid E, Harwell J, Coalson T, Schindler J, Elam JS, Curtiss SW, Van Essen DC. Human connectome project informatics: quality control, database services, and data visualization. Neuroimage. 2013;80:202–19.

    Article  PubMed  Google Scholar 

  27. Marcus DS, Harwell J, Olsen T, Hodge M, Glasser MF, Prior F, Jenkinson M, Laumann T, Curtiss SW, Van Essen DC. Informatics and data mining tools and strategies for the human connectome project. Front Neuroinform. 2011;5:1–12.

    Article  Google Scholar 

  28. Oosterhof NN, Connolly AC, Haxby JV, Rosa MJ. CoSMoMVPA : multi-modal multivariate pattern analysis of neuroimaging data in Matlab / GNU octave. Front Neuroinform. 2016;10:1–27.

    Article  Google Scholar 

  29. Smith SM. Fast robust automated brain extraction. Hum Brain Mapp. 2002;17:143–55.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Jenkinson M, Bannister P, Brady M, Smith S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage. 2002;17:825–41.

    Article  PubMed  Google Scholar 

  31. Alakörkkö T, Saarimääki H, Glerean E, Saramäki J, Korhonen O. Effects of spatial smoothing on functional brain networks. Eur J Neurosci. 2017;46:2471–80.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Andersson J, Jenkinson M, Smith S (2010): Non-linear registration, aka spatial normalisation. FMRIB Tech Rep TR07JA2.

  33. Pruim RHR, Mennes M, van Rooij D, Llera A, Buitelaar JK, Beckmann CF. ICA-AROMA: a robust ICA-based strategy for removing motion artifacts from fMRI data. Neuroimage. 2015;112:267–77.

    Article  PubMed  Google Scholar 

  34. Tucholka A, Fritsch V, Poline J-B, Thirion B. An empirical comparison of surface-based and volume-based group studies in neuroimaging. Neuroimage. 2012;63:1443–53.

    Article  PubMed  Google Scholar 

  35. Dickie EW, Anticevic A, Smith DE, Coalson TS, Manogaran M, Calarco N, Viviano JD, Glasser MF, Van Essen DC, Voineskos AN. Ciftify: a framework for surface-based analysis of legacy MR acquisitions. Neuroimage. 2019;197:818.

    Article  PubMed  Google Scholar 

  36. Beckmann CF, Mackay CE, Filippini N, Smith SM. Group comparison of resting-state FMRI data using multi-subject ICA and dual regression. NeuroImage. 2009.

    Article  PubMed  Google Scholar 

  37. Tahedl M, Schwarzbach JV. An updated and extended atlas for corresponding brain activation during task and rest. Prepr bioRxive. 2020.

    Article  Google Scholar 

  38. Nichols TE, Holmes AP. Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp. 2001;25:1–25.

    Google Scholar 

  39. Faul F, Erdfelder E, Lang A-G, Buchner A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. 2007;39:175–91.

    Article  PubMed  Google Scholar 

  40. Rocca MA, Valsasina P, Martinelli V, Misci P, Falini A, Comi G, Filippi M. Large-scale neuronal network dysfunction in relapsing-remitting multiple sclerosis. Neurology. 2012;79:1449–57.

    Article  PubMed  Google Scholar 

  41. Amato MP, Ponziani G. A prospective study on the prognosis of multiple sclerosis. Neurol Sci. 2000;21:S831-8.

    Article  CAS  PubMed  Google Scholar 

  42. Kutzelnigg A, Faber-Rod JC, Bauer J, Lucchinetti CF, Sorensen PS, Laursen H, Stadelmann C, Brück W, Rauschka H, Schmidbauer M, Lassmann H. Widespread demyelination in the cerebellar cortex in multiple sclerosis. Brain Pathol. 2007;17:38–44.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Cocozza S, Pontillo G, Russo C, Russo CV, Costabile T, Pepe A, Tedeschi E, Lanzillo R, Brescia Morra V, Brunetti A, Inglese M, Petracca M. Cerebellum and cognition in progressive MS patients: functional changes beyond atrophy? J Neurol. 2018;265:2260–6.

    Article  PubMed  Google Scholar 

  44. Dogonowski A-M, Andersen KW, Madsen KH, Sørensen PS, Paulson OB, Blinkenberg M, Siebner HR. Multiple sclerosis impairs regional functional connectivity in the cerebellum. NeuroImage Clin. 2014;4:130–8.

    Article  PubMed  Google Scholar 

  45. Sbardella E, Upadhyay N, Tona F, Prosperini L, De Giglio L, Petsas N, Pozzilli C, Pantano P. Dentate nucleus connectivity in adult patients with multiple sclerosis: functional changes at rest and correlation with clinical features. Mult Scler J. 2016;23:546–55.

    Article  Google Scholar 

  46. Marr D. A theory of cerebellar cortex. J Physiol. 1969;202:437–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Albus S. A theory of cerebellar function. Math Biosci. 1971;10:25–61.

    Article  Google Scholar 

  48. Ito M. Error detection and representation in the olivo-cerebellar system. Front Neural Circuits. 2013;7:1–8.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Brodal P. The corticopontine projection in the rhesus monkey origin and principles of organization. Brain. 1978;101:251–83.

    Article  CAS  PubMed  Google Scholar 

  50. Glickstein M, May JG, Mercier BE. Corticopontine projection in the macaque: the distribution of labelled cortical cells after large injections of horseradish peroxidase in the pontine nuclei. J Comp Neurol. 1985;235:343–59.

    Article  CAS  PubMed  Google Scholar 

  51. Ramnani N. Automatic and controlled processing in the corticocerebellar system. Prog Brain Res. 2014.

    Article  PubMed  Google Scholar 

  52. King M, Hernandez-Castillo CR, Poldrack RA, Ivry RB, Diedrichsen J. Functional boundaries in the human cerebellum revealed by a multi-domain task battery. Nat Neurosci. 2019;22:1371–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Saini S, DeStefano N, Smith S, Guidi L, Amato MP, Federico A, Matthews PM. Altered cerebellar functional connectivity mediates potential adaptive plasticity in patients with multiple sclerosis. J Neurol Neurosurg Psychiatry. 2004;75:840–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Cerasa A, Gioia MC, Valentino P, Nisticò R, Chiriaco C, Pirritano D, Tomaiuolo F, Mangone G, Trotta M, Talarico T, Bilotti G, Quattrone A. Computer-assisted cognitive rehabilitation of attention deficits for multiple sclerosis: a randomized trial with fMRI correlates. Neurorehabil Neural Repair. 2013;27:284–95.

    Article  PubMed  Google Scholar 

  55. Leavitt VM, Wylie GR, Girgis PA, DeLuca J, Chiaravalloti ND. Increased functional connectivity within memory networks following memory rehabilitation in multiple sclerosis. Brain Imaging Behav. 2014;8:394–402.

    Article  PubMed  Google Scholar 

  56. Khan F, Amatya B. Rehabilitation in multiple sclerosis: a systematic review of systematic reviews. Arch Phys Med Rehabil. 2017;98:353–67.

    Article  PubMed  Google Scholar 

  57. Prosperini L, Piattella MC, Giannì C, Pantano P. Functional and structural brain plasticity enhanced by motor and cognitive rehabilitation in multiple sclerosis. Neural Plast. 2015;2015:481574.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Drysdale AT, Grosenick L, Downar J, Dunlop K, Mansouri F, Meng Y, Fetcho RN, Zebley B, Oathes DJ, Etkin A, Schatzberg AF, Sudheimer K, Keller J, Mayberg HS, Gunning FM, Alexopoulos GS, Fox MD, Pascual-Leone A, Voss HU, Casey BJ, Dubin MJ, Liston C. Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nat Med. 2017;23:28–38.

    Article  CAS  PubMed  Google Scholar 

  59. Fox MD, Snyder AZ, Vincent JL, Corbetta M, Van Essen DC, Raichle ME. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci USA. 2005;102:9673–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Daniel H, Levenes C, Crépel F. Cellular mechanisms of cerebellar LTD. Trends Neurosci. 1998;21:401–7.

    Article  CAS  PubMed  Google Scholar 

  61. Hirano T. Depression and potentiation of the synaptic transmission between a granule cell and a Purkinje cell in rat cerebellar culture. Neurosci Lett. 1990;119:141–4.

    Article  CAS  PubMed  Google Scholar 

  62. Sakurai M. Synaptic modification of parallel fibre-Purkinje cell transmission in in vitro guinea-pig cerebellar slices. J Physiol. 1987;394:463–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Wang DJ, Da SuL, Wang YN, Yang D, Sun CL, Zhou L, Wang XX, Shen Y. Long-term potentiation at cerebellar parallel fiber-Purkinje cell synapses requires presynaptic and postsynaptic signaling cascades. J Neurosci. 2014;34:2355–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Mori F, Kusayanagi H, Nicoletti CG, Weiss S, Marciani MG, Centonze D. Cortical plasticity predicts recovery from relapse in multiple sclerosis. Mult Scler J. 2014;20:451–7.

    Article  Google Scholar 

  65. Weiss S, Mori F, Rossi S, Centonze D. Disability in multiple sclerosis: when synaptic long-term potentiation fails. Neurosci Biobehav Rev. 2014.

    Article  PubMed  Google Scholar 

  66. Álvarez-Salvado E, Pallarés V, Moreno A, Canals S. Functional MRI of long-term potentiation: imaging network plasticity. Philos Trans R Soc B Biol Sci. 2014;369:20130152.

    Article  Google Scholar 

  67. Canals S, Beyerlein M, Merkle H, Logothetis NK. Functional MRI evidence for LTP-induced neural network reorganization. Curr Biol. 2009;19:398–403.

    Article  CAS  PubMed  Google Scholar 

  68. Barkhof F. The clinico-radiological paradox in multiple sclerosis revisited. Curr Opin Neurol. 2002;15:239–45.

    Article  PubMed  Google Scholar 

  69. Liu Y, Wang H, Duan Y, Huang J, Ren Z, Ye J, Dong H, Shi F, Li K, Wang J. Functional brain network alterations in clinically isolated syndrome and multiple sclerosis: a graph-based connectome study. Radiology. 2017;282:534–41.

    Article  PubMed  Google Scholar 

  70. McIntosh AR. Towards a network theory of cognition. Neural Netw. 2000;13:861–70.

    Article  CAS  PubMed  Google Scholar 

  71. Steinman L. Immunology of relapse and remission in multiple sclerosis. Annu Rev Immunol. 2014;32:257–81.

    Article  CAS  PubMed  Google Scholar 

  72. Yeo BTT, Krienen FM, Sepulcre J, Sabuncu MR, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysio. 2011;106:1125–65.

    Article  Google Scholar 

Download references


The authors would like to thank Philipp Seidel and Regina Rockinger for assistance with data acquisition.


MT and JVS were funded by the Deutsche Multiple Sklerose Gesellschaft (DMSG), grant number 2018_DMSG_08.

Author information

Authors and Affiliations



MT, RW, SML and JVS conceived of and designed the experiment. MT, RW, ZK, DHL, and RAL recruited patients and performed clinical ratings. MT, SML, and JVS acquired data. MT analyzed the neuroimaging data. MT and SML performed statistical analyses, and created the figures, and drafted the manuscript. MT, SML, and JVS extended the conceptual framework. All authors contributed to manuscript revision and approved the final version.

Corresponding author

Correspondence to Jens V. Schwarzbach.

Ethics declarations

Ethics approval and consent to participate

All patients provided written informed consent to participate in this study. All study procedures followed safety guidelines for MRI research at the University of Regensburg, complied with the 1964 Helsinki Declaration and its later amendments, and were approved by the ethics committee of the University of Regensburg.

Consent for publication

Not applicable.

Competing interests

MT, SML, RW, ZK, and JVS declare no conflicts of interest. RAL and DHL received compensation for activities with Bayer, Biogen, Celgene, Genzyme, Merck, Novartis, Roche and TEVA as well as research support from Merck and Novartis.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Figure S1. Probability that a given network is the “original network”, whose functional coherence is changed in MS patients following their entrance into remission. Open circles are data from individual patients; closed circles depict the mean, and error bars are SEM. (A) Results from the first set of 9 patients. (B) Results from the second set of 9 patients. Note that in both cases the cerebellar network is, by far, the most likely network to change its functional coherence.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tahedl, M., Levine, S.M., Weissert, R. et al. Early remission in multiple sclerosis is linked to altered coherence of the Cerebellar Network. J Transl Med 20, 488 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Cerebellar Network
  • Functional connectivity
  • Independent component analysis
  • Multiple sclerosis
  • Remission