Altered anterior default mode network dynamics in progressive multiple sclerosis ================================================================================ * Giulia Bommarito * Anjali Tarun * Younes Farouj * Maria Giulia Preti * Maria Petracca * Amgad Droby * Mohamed Mounir El Mendili * Matilde Inglese * Dimitri Van De Ville ## Abstract **Background** Modifications in brain function remain relatively unexplored in progressive multiple sclerosis (PMS), despite their potential to provide new insights into the pathophysiology of this disease stage. **Objectives** To characterize the dynamics of functional networks at rest in patients with PMS, and the relation with clinical disability. **Methods** Thirty-two patients with PMS underwent clinical and cognitive assessment. The dynamic properties of functional networks, retrieved from transient brain activity, were obtained from patients and 25 healthy controls (HC). Sixteen HC and 19 patients underwent a one-year follow-up clinical and imaging assessment. Differences in the dynamic metrics between groups, their longitudinal changes, and the correlation with clinical disability were explored. **Results** PMS patients, compared to HC, showed a reduced dynamic functional activation of the anterior default mode network (aDMN) and its opposite-signed coactivation with the executive-control network, at baseline and follow-up. Processing speed and visuo-spatial memory negatively correlated to aDMN dynamic activity. The anti-coupling between aDMN and auditory/sensory-motor network, temporal-pole/amygdala or salience networks were differently associated to separate cognitive domains. **Conclusion** Patients with PMS presented an altered aDMN functional recruitment and anti-correlation with ECN. The aDMN dynamic functional activity and interaction with other networks explained cognitive disability. Keywords * progressive multiple sclerosis * resting-state functional MRI * innovation-driven co-activation patterns * default mode network * clinical disability ## Introduction Progressive multiple sclerosis (PMS) is characterized by a gradual worsening in clinical disability, occurring from the onset of the disease in the case of primary progressive (PP) form, or after a relapsing-remitting (RR) phase in secondary progressive (SP) multiple sclerosis (MS). PP and SP MS phenotypes are considered different aspects of a unique spectrum1 and this disease stage, compared to RR MS, presents with a slower, chronic, compartmentalized inflammation and a prominent neurodegeneration, together with a depleted ability to respond to damage. The contribution of imaging to the definition of PMS ranges from the identification of biomarkers, as atrophy, to the detection, in vivo, of hallmarks such as cortical lesions and diffuse white matter damage1,2. Nonetheless, while a growing number of reports is revealing the clinical significance of the networks function disruption occurring in MS, studies targeting brain functional activity in PMS are relatively lacking. Previous studies of brain function at rest in patients with PMS reported altered network activity or connectivity involving the default mode network (DMN)3, attentional, executive control and sensory-motor networks4–6. However, the discordant results and the limited number of studies does not allow a unique interpretation of the changes in functional networks occurring in this phase of the disease and the hypothesis of an exhaustion of plasticity characterizing PMS is mainly inferred by studies on RR patients at later stages. Moreover, to our knowledge, only a longitudinal study focusing on the sensory-motor functional network has been performed, so far7. As emerging for other neurodegenerative conditions, defining the functional disruption occurring in PMS could contribute to the understanding of pathophysiology, to extricate the changes due to the disease or to brain aging and to pinpoint potential treatments8. Beside conventional analysis, new approaches have recently emerged to explore the time-varying nature of brain functional activity at rest9. In healthy subjects, the investigation of the dynamic properties of rs-FC has revealed its ability in better grasping behavioral aspects, compared to the classic time-averaged connectivity10. The dynamics of rs-FC have also been investigated in patients at the earlier stages of MS, emphasizing the complex functional changes occurring in this disorder, and their relationship to the cognitive status11–16. Moreover, the ability of dynamic rs-FC in better capturing cognitive decline has been showed also in neurodegenerative disorders17. Multiple methodologies have been developed to study the variations over time in functional brain activity. Among them, the characterization of patterns of functional co-activation among regions has been demonstrated able to tackle the pathophysiological processes occurring in neurological or neuropsychiatric disorders, and to parallel the clinical improvement after therapeutic procedures18–20 In this study, we hypothesized that the dynamics of functional networks, investigated by means of co-activation patterns, would be altered in patients with PMS and would explain clinical disability. Specifically, we investigated: i) the dynamics of functional networks in patients with PMS, compared to healthy controls (HC); ii) their changes after one-year of follow-up (FU), and iii) their role in explaining cognitive impairment, assessed with the Brief International Cognitive Assessment for Multiple Sclerosis, and clinical disability. ## Materials and Methods ### Participants Forty-eight patients with either PP or SP MS and 26 healthy controls (HC) were prospectively recruited. Inclusion and exclusion criteria are detailed in the Supplementary Information (SI). After image quality check and motion correction (detailed in SI), the final sample included 57 subjects: 32 patients with PMS (mean age: 50.8±9.9 years, 20 women, 18 patients with PP MS and 14 patients with SP MS) and 25 HC (mean age: 47.3±8.7 years, 10 women). All subjects underwent MRI and, on the same day, patients underwent clinical evaluation, including expanded disability status scale (EDSS), timed-25-foot-walking (T25FW) test, 9-hole-peg-test (9HPT), and neuropsychological assessment with the symbol digit modalities test (SDMT), California verbal learning test second edition (CVLT-II) and the brief visuospatial memory test revised (BVMT-R). For the CVLT-II and the BVMT-R the first five or three recall trials were considered, respectively. Raw scores were transformed into Z-scores according to21. Impairment for a single test was assessed on the level of 1.5 SD, and patients were defined as cognitively impaired (CI) or preserved (CP) based on scoring outside the normal range in one or more of the tests22. Thirty-five (16 HC and 19 patients, 11 with PP MS and 8 with SP MS) out of the 57 subjects included at baseline underwent MRI and clinical evaluation after 1-year ([mean ±SD]= 11±2 months). Subjects who were included at FU or dropped out from the study did not differ in terms of demographic, clinical and imaging variable (see Table SI 1). The study was approved by the Institutional Review Board of the Icahn School of Medicine at Mount Sinai, and all the subjects gave written informed consent to participate. ### MRI acquisition All the participants underwent MRI on the same 3.0□T scanner (Magnetom Skyra, Erlangen, Siemens, Germany) using a 32-channel head coil, with a constant protocol at baseline and FU, including T2-weighted, T1-weighted, and a single shot gradient echo planar imaging sequence for the resting state. Sequence parameters are reported in the SI. ### Structural MRI analysis T2 and T1 lesions were segmented, and lesion volumes (T2 LV and T1 LV, respectively) obtained using a semiautomatic technique (Jim 7, Xinapse Systems, Northants, UK). The T1-weighted images of PMS patients underwent a lesion in-painting procedure prior to the brain extraction and anatomical structures’ segmentation23. The baseline normalized brain volume (BV), grey matter (GMV) and white matter volumes (WMV), and the FU percentage brain volume changes were estimated using SIENAX and SIENA, respectively24. ### Functional resting state MRI preprocessing and iCAPs analysis Resting state functional magnetic resonance imaging (rs-fMRI) data were preprocessed using SPM12 (FIL, UCL, UK), complementary functions of the Data Processing Assistant for Resting-State fMRI (DPARSF)25 and the Individual Brain Atlases using Statistical Parametric Mapping toolboxes26. Preprocessing steps included: realignment and discarding of the first 10 volumes, resulting in 390 time-points, co-registration of the T1-weighted and functional images, segmentation of the structural images into white matter (WM), grey matter and cerebrospinal fluid (CSF). The regression of the nuisance variables, including the average WM and CSF signals and six motion parameters, was performed using DPARSF toolbox. A spatial smoothing (Gaussian kernel of a 6 mm full width at half maximum) was finally applied. The innovation-driven-coactivation-patterns (iCAPs) analysis is a novel technique that allows to derive a set of whole-brain spatial patterns of regions whose activity simultaneously increases or decreases, thus characterized by similar functional dynamic behavior. These are obtained by detecting the *transients* or changes in rs-fMRI activity time-courses and each iCAP represents a brain functional *state* or network. The procedure was performed using a publicly available implementation ([https://c4science.ch/source/iCAPs/](https://c4science.ch/source/iCAPs/)). For a comprehensive explanation of the methodology, please refer to27, to SI, and Figure 1. The iCAPs retrieved from the baseline sample were fitted into the FU sample, as detailed in the SI. From now on, the terms “iCAPs”, “networks” or “states” will be used interchangeably to describe the spatial patterns of functional activations derived from the iCAPs analysis. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/28/2020.11.26.20238923/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/04/28/2020.11.26.20238923/F1) Figure 1. Top: fMRI time-courses undergo a deconvolution step to obtain activity-inducing signals. Innovation signals, obtained by computing the temporal derivative of the activity-inducing signals, enter a two-step thresholding procedure (temporal and spatial) and the resulting frames undergo the temporal clustering step to extract iCAPs. Bottom: the total duration for each iCAP was computed considering both the positive and negative durations, i.e. the time-points when that iCAP was active (after thresholding using a cut-off of |1|). The overlap between 2 iCAPs was measured with the jaccard score, obtained dividing the number of time-points when the two iCAPs were simultaneously active for the number of time-points when at least one of the two iCAPs was active. The couplings and anti-couplings represented the time-points when the 2 iCAPs co- (de)activated (A and C) or showed an opposite co(de)-activation (B), respectively. ### iCAPs temporal properties For each subject at baseline and FU, we obtained iCAPs’ temporal characteristics: i) the *duration* of overall functional activation, a metric describing the amount of sustained activity (or functional “engagement”) for that specific iCAP; for each pair of iCAPs ii) the *coupling*: same-signed coactivation, and iii) the *anti-coupling*: opposite-signed coactivation (Figure 1). Couplings and anti-couplings reflect the amount of overlap between the activity, in the same or the opposite signs respectively, of each iCAPs pair, and are closely related to the networks’ static functional connectivity. ### Statistical analysis Demographic and clinical parameters were compared using unpaired and paired t-test between HC and patients with PMS or between values at baseline and FU, respectively. Analysis of covariance (ANCOVA) was used to compare iCAPs duration between HC and PMS patients and between CI and CP patients, using age, gender and education as covariates. A false discovery rate (FDR) correction according to28 was applied. To compare FU and baseline iCAPs temporal properties, a paired t-test within each group and a repeated measure ANCOVA analysis, considering the group x time interaction, were used. Further analyses on coupling/anti-coupling features were computed only for iCAPs whose duration was significantly different in PMS patients compared to HC. ICAPs temporal metrics were correlated with T2 and T1 LV using Spearman’s rank correlation coefficient. Partial least squares correlation (PLSC) analysis was used to investigate patterns of correlation between clinical parameters and the temporal properties of iCAPs. Age, gender and education were stepwise regressed from each variable, except for the normalized cognitive tests Z scores. The toolbox used to run PLSC analysis is publicly available ([https://miplab.epfl.ch/index.php/software/PLS](https://miplab.epfl.ch/index.php/software/PLS)). Since aDMN temporal properties resulted altered in patients with PMS compared to HC (see Results), we assessed their relationship with clinical disability. Grouped PLSC analyses were performed to assess the different contribution of aDMN metrics to cognitive test scores in CI or CP patients. Moreover, we performed an exploratory analysis to assess any relationship between all iCAPs durations and clinical test scores. ## Data availability Data supporting the findings of this study will be shared upon request. ## Results ### Demographic and clinical data Demographic data for the two groups and clinical data for the patients with PMS are reported in Table 1. Table SI 1 reports data for the subjects who underwent FU and drop out subjects. No significant differences in demographic variables were observed between HC and PMS, while BV, WMV and GMV were significantly lower in PMS patients. PMS patients did not differ in terms of any clinical parameter of disability between baseline and FU (Table SI 2). Twenty-four patients were classified as CI at baseline, and the PMS group showed low normalized scores at the BVMT and SDMT tests (see also SI and Figure SI 3). View this table: [Table 1.](http://medrxiv.org/content/early/2021/04/28/2020.11.26.20238923/T1) Table 1. Demographic and clinical data at baseline and follow-up. ### Spatial properties of iCAPs Eleven iCAPs with different spatial distribution were obtained and considered for subsequent analysis. The networks represented in each iCAP are shown in Figure 2 and named based on the observed spatial pattern and its overlap with well-known functional networks (see Table SI 3). Specifically, they represented the deep grey matter (DGM), auditory/sensory-motor (Aud/SM), primary visual (PrimVis), executive control (ECN), anterior DMN (aDMN), salience network (SAL), temporo-parietal/language (TempPar/Lan), secondary visual (SecVis), precuneus/posterior DMN (pCun/pDMN), amygdala/temporal pole (Amy/TP) and orbito-frontal cortex (OFC). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/28/2020.11.26.20238923/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/04/28/2020.11.26.20238923/F2) Figure 2. Spatial patterns of the 11 iCAPs retrieved from the analysis on HC and patients with PMS. Under each iCAP are reported the average consensus, and the number of innovation frames assigned to it. Coordinates refer to the Montreal Neurological Institute space. DGM: deep grey matter; Aud/SM: the auditory/sensory-motor; PrimVis: primary visual; ECN: executive control network; aDMN: anterior DMN; SAL: salience; TempPar/Lan: temporo-parietal/language; SecVis: secondary visual; pCun/pDMN: precuneus/posterior DMN; Amy/TP: amygdala/temporal pole; OFC: orbito-frontal cortex. ### Temporal properties of iCAPs: group differences #### Baseline The overall duration of the aDMN was significantly reduced in patients with PMS compared to HC (*F*(1,52)=10.99, p= .002, Figure 3A, Table SI 5). A significantly reduced anti-coupling between aDMN and ECN was also observed in patients with PMS compared to HC (*F*(1,52)=7.2, p=.010, Figure 3B, Table SI 6 and 7). No significant differences were detected when iCAPs temporal metrics were compared between CI and CP patients (Table SI 8, 9 and 10). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/28/2020.11.26.20238923/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/04/28/2020.11.26.20238923/F3) Figure 3. *A*: bar plots of the duration of the 11 iCAPs in HC and patients with PMS. *B*: anti-coupling between ECN and aDMN, resulted significantly different between the HC and PMS patients. DGM: deep grey matter; Aud/SM: the auditory/sensory-motor; PrimVis: primary visual; ECN: executive control network; aDMN: anterior DMN; SAL: salience; TempPar/Lan: temporo-parietal/language; SecVis: secondary visual; pCun/pDMN: precuneus/posterior DMN; Amy/TP: amygdala/temporal pole; OFC: orbito-frontal cortex. *p<0.05, ** surviving FDR correction. #### FU The longitudinal analysis did not reveal differences in iCAPs durations in PMS patients or HC at FU vs baseline while the aDMN-ECN coupling significantly increased in PMS patients at FU compared to baseline (*t*(18)=-3.86, *p*= .001, Figure SI 4). No significant group x time interaction was found for the 11 networks durations or for the couplings/anti-couplings involving the aDMN. A group effect was confirmed at FU for aDMN duration (*F*(1,30)=14.5, *p*= .001) and the anti-coupling between aDMN and ECN (*F*(1,30)=12.0, *p*= .002). Meanwhile, a group effect emerged for Aud/SM (*F*(1,30)=13.3, *p*= .001) and the aDMN-Aud/SM anti-coupling (*F*(1,30)=4.2, *p*= .049). ### Temporal properties of iCAPs and structural MRI metrics No significant associations were found between aDMN duration at baseline and T2 or T1 LV. The aDMN-pCun/pDMN anti-coupling and the coupling between aDMN and primVis were positively correlated to T1 LV and T2 LV (*r*(27)= .40, *p*= .04) (*r*(27)= .41, *p*= .03, respectively). No correlations were found between aDMN temporal metrics and volumetric parameters. ### iCAPs temporal properties and clinical parameters at baseline #### aDMN duration No significant results were found for EDSS or motor scores. Cognitive scores were associated to aDMN durations (p=0.030, Figure 4). Specifically, a negative correlation was observed for the SDMT and BVMT-R. No significant results emerged from the grouped PLSC analysis (CI and CP patients considered separately). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/28/2020.11.26.20238923/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2021/04/28/2020.11.26.20238923/F4) Figure 4. Partial least squares correlation analyses performed to assess the relationship between aDMN duration (top) or anti-couplings (bottom) and cognitive test Z scores. In yellow are shown the robust weights after bootstrapping, and the error bars represent bootstrap intervals. DGM: deep grey matter; Aud/SM: the auditory/sensory-motor; PrimVis: primary visual; ECN: executive control network; aDMN: anterior DMN; SAL: salience; TempPar/Lan: temporo-parietal/language; SecVis: secondary visual; pCun/pDMN: precuneus/posterior DMN; Amy/TP: amygdala/temporal pole; OFC: orbito-frontal cortex; BVMT-R: brief visual memory test revised; CVLT-II: California verbal learning test (second edition); SDMT: symbol digit modalities test. #### aDMN anti-couplings A significant association between aDMN anti-couplings and cognitive performance was found (p=0.005, Figure 4) in the whole patients group. In particular, SDMT and CVLT-II were related to aDMN-Aud/SM, aDMN-Amy/TempP, and aDMN-SAL, with the result mainly driven by positive correlations between aDMN-Aud/SM or aDMN-SAL and CVLT-II score and the negative correlations between aDMN-Aud/SM and SDMT and between aDMN-Amy/TempP and CVLT-II. When the grouped PLSC was performed, separating CI and CP patients, we found that aDMN anti-couplings were related to cognitive performance in the CI, but not in the CP group (p= .019, Figure SI 5). #### All iCAPs durations This exploratory analysis additionally revealed that the durations of DGM, Aud/SM, TempPar/LAN and pCun/pDMN were significantly correlated to EDSS and T25FW (*p*= .027, Figure SI 6). ### Reproducibility analyses Analyses to verify the reproducibility of the results on study sub-samples are reported in the SI. ## Discussion In this work, we characterized the changes in the dynamics of brain functional patterns occurring in patients with PMS patients, showing that they present changes in the functional activation of the aDMN, and in its interplay with ECN. Such changes were detected also over one year FU, in a sample of patients who did not show a significant worsening in clinical disability. Moreover, the dynamic temporal characteristics of the aDMN, namely its reduced sustained activity and the interaction with Aud/SM, SAL and Amy/TempPole, explained cognitive disability in patients with PMS. Regarding the methodology, we chose to apply this specific data-driven approach since allows to overcome the limitation of the temporal and spatial segregation of functional networks, which characterizes other techniques utilized to investigate functional dynamics, thus better grasping the dynamic interplay between networks over the time of a rs-fMRI sequence. Moreover, these spatial patterns of functional activity have been showed to be reliable across subjects29. As first contribution, this study revealed that the aDMN overall activity duration and its anti-correlation with ECN activity were reduced in patients with PMS. Previous studies showed a reduced aDMN conventional functional activity, related to cognitive disability, in PMS3, and highlighted its functional changes in other neurodegenerative disorders30. The ECN and DMN are usually oppositely engaged and involved in the switching between the “idling” vs “task-oriented” condition. Besides, together with the SAL, they have been associated to clinical phenotypes in neurodegenerative and psychiatric conditions. Here, we confirm the central role of aDMN dysfunction in the progressive phase of MS and advance the hypothesis an unbalanced and less dynamic interaction with ECN could be linked to the underlying neurodegenerative process. Intriguingly, the anti-coupling between aDMN and pCun/pDMN as well as its coupling with the primVis network were correlated to lesion load, endorsing the connection between structural damage and functional disruption. A second relevant finding was the absence of significant changes in functional dynamics over time, accompanied by the lack of clinical worsening. Although we cannot exclude this result could be impacted by attrition bias or learning effect for cognitive tests scores, previous studies in PMS propend for slow changes in clinical and structural imaging metrics31,32 and the longitudinal variation of brain functional activity might be as gradual. Nonetheless, further longitudinal studies with larger sample sizes are needed to address whether more subtle brain functional modifications occur in the late progressive phase. As last result, we found a strong relationship between altered functional dynamics of the aDMN and cognitive disability. The DMN plays a pivotal role in the cognitive disability in patients with MS3,6,33 and previous studies on the dynamic rs-FC in patients with MS highlighted its role in explaining cognitive impairment, including processing speed11,13. Our results indicated that an increased engagement of the aDMN was associated to a worse cognitive performance. Besides, the cognitive domains expressed by the three tests were differently explained by the interaction between aDMN and the other networks. Namely, a better performance in verbal memory, the cognitive domain less affected in our patients’ population, was associated to an increased anti-coupling between aDMN and both Aud/SM and the SAL networks, suggesting the synchronous opposite-signed activation contributes to ensure a better cognitive functioning, while the opposite was observed for the interaction between aDMN and Amy/TempP. Besides, the aDMN-Aud/SM reduced segregation correlated with a better processing speed, suggesting a compensatory role for the aDMN, likely involving its synchronicity with sensory-motor regions34. Data from literature suggests an influence of functional dynamics during rest on the subsequent performance during tasks35, hence a possible interpretation of the reduced aDMN dynamics as compensatory, in an attempt to adapt to cognitive tasks demand. Another interpretation entails the aDMN interplay with other networks and the reduced overall duration as an attempt to respond to structural damage reducing the dynamic fluidity to privilege efficiency. When a grouped analysis was performed, the aDMN anti-couplings explained the cognitive status in the CI patients, supporting the idea that changes in behavioral performance are better seized by the dynamic interaction among functional networks. In particular, it emerged an additional relationship between a reduced aDMN-pCun/pDMN anti-coupling and better visuo-spatial memory and processing speed. The latter also correlated to a reduced opposite-signed activity of aDMN with ECN and increased anti-correlation with the higher visual network. An additional exploratory analysis showed a relationship between T25FW and EDSS and the Aud/SM, the TempPar, DGM and pCun/pDMN activity duration. Although none of these networks temporal properties significantly differed, compared to HC, their relative changes in functional engagement might underlie the behavioral modifications characterizing disability in progressive MS. This study comes with several limitations. First, the small sample size and we cannot exclude type II errors. Second, the high rate of drops-out, which could have biased the results and determined the lack of significant changes we observe over one-year FU. Third, learning effects could have affected the longitudinal results, thus the absence of a significant worsening in cognitive tests at FU compared to baseline. Fourth, motor and cognitive test scores were not available for the HC group, which urges caution in the interpretation of the correlation analyses. Fifth, the use of a multi-band acceleration factor of 7 that could have introduced noise in the rs-fMRI data. Sixth, the lack of an external validation, although the analyses on the sub-samples of our study population proved stable and reproducible results. Seventh, the absence of a population of patients at the earlier stages of the disease, which prevents us to conclude that our findings are specific of the progressive stages of MS. In conclusion, this study characterized functional brain dynamics in patients with PMS, revealing a reduced engagement of the aDMN *state*, though not detrimental for the cognitive status which is additionally related to the dynamic interplay between aDMN and other functional networks. Future studies are needed to confirm these findings and to address the coupling between functional modifications and local or diffuse structural damage through the different disease stages. ## Supporting information Supplementary Information [[supplements/238923_file06.pdf]](pending:yes) ## Data Availability Data supporting the findings of this study will be shared upon request. ## Disclosures Giulia Bommarito, Anjali Tarun, Younes Farouj, Maria Giulia Preti, Maria Petracca, Amgad Droby, Mohamed Mounir El Mendili, and Dimitri Van De Ville have nothing to disclose. Matilde Inglese is currently receiving a grant from Teva Neuroscience. ## Funding G. Bommarito was supported by a research fellowship FISM-Fondazione Italiana Sclerosi Multipla, Cod.: 2017/B/2 and financed or co-financed with the ‘5 per mille’ public funding. A. Tarun and Y. Farouj were supported by the Swiss National Science Foundation under the Project Grant 205321-163376. ## Abbreviations PMS : progressive multiple sclerosis PP : primary progressive SP : secondary progressive fMRI : functional MRI rs-fMRI : resting state functional MRI DMN : default mode network rs-FC : resting state functional connectivity iCAP : innovation driven co-activation pattern BOLD : blood-oxygenated-level dependent FU : follow-up HC : healthy control EDSS : expanded disability status scale T25FW : timed 25-foot walking 9HPT : 9-hole peg test SDMT : symbol digit modalities test CVLT-II : California verbal learning test second edition BVMT-R : brief visuospatial memory test revised WM : white matter CSF : cerebrospinal fluid FD : framewise displacement TA : total activation MNI : Montreal Neurological Institute ANCOVA : analysis of covariance FDR : false discovery rate PLSC : partial least squares correlation DGM : deep grey matter Aud/SM : the auditory/sensory-motor PrimVis : primary visual ECN : executive control network aDMN : anterior DMN SAL : salience TempPar/Lan : temporo-parietal/language SecVis : secondary visual pCun/pDMN : precuneus/posterior DMN Amy/TP : amygdala/temporal pole OFC : orbito-frontal cortex. * Received November 26, 2020. * Revision received April 27, 2021. * Accepted April 28, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Lassmann H. Pathogenic mechanisms associated with different clinical courses of multiple sclerosis. Front Immunol. 2019;10:1–14. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2019.00858&link_type=DOI) 2. 2.Filippi M, Preziosa P, Langdon D, et al. Identifying Progression in Multiple Sclerosis: New Perspectives. Ann Neurol. 2020;88:438–452. 3. 3.Rocca MA, Valsasina P, Absinta M, et al. Default-mode network dysfunction and cognitive impairment in progressive MS. Neurology. United States; 2010;74:1252–1259. 4. 4.Petracca M, Saiote C, Bender HA, et al. Synchronization and variability imbalance underlie cognitive impairment in primary-progressive multiple sclerosis. Sci Rep. England; 2017;7:46411. 5. 5.Rocca MA, Valsasina P, Leavitt VM, et al. Functional network connectivity abnormalities in multiple sclerosis: Correlations with disability and cognitive impairment. Mult Scler J. 2018;24:459–471. 6. 6.Basile B, Castelli M, Monteleone F, et al. Functional connectivity changes within specific networks parallel the clinical evolution of multiple sclerosis. Mult Scler J. 2014;20:1050– 1057. 7. 7.Droby A, Fleysher L, Petracca M, et al. Lower cortical gamma-aminobutyric acid level contributes to connectivity in sensory-motor inter-connected regions in progressive MS. Mult Scler Relat Disord [online serial]. Elsevier B.V.; 2020;43:102183. Accessed at: [https://doi.org/10.1016/j.msard.2020.102183](https://doi.org/10.1016/j.msard.2020.102183). 8. 8.Chard DT, Alahmadi AAS, Audoin B, et al. Mind the gap: from neurons to networks to outcomes in multiple sclerosis. Nat Rev Neurol. Epub 2021. 9. 9.Preti MG, Bolton TA, Van De Ville D. The dynamic functional connectome: State-of-the-art and perspectives. Neuroimage. 2017;160:41–54. 10. 10.Liégeois R, Li J, Kong R, et al. Resting brain dynamics at different timescales capture distinct aspects of human behavior. Nat Commun [online serial]. Springer US; 2019;10. Accessed at: [http://dx.doi.org/10.1038/s41467-019-10317-7](http://dx.doi.org/10.1038/s41467-019-10317-7). 11. 11.van Geest Q, Douw L, van ‘t Klooster S, et al. Information processing speed in multiple sclerosis: Relevance of default mode network dynamics. NeuroImage Clin [online serial]. Elsevier; 2018;19:507–515. Accessed at: [https://doi.org/10.1016/j.nicl.2018.05.015](https://doi.org/10.1016/j.nicl.2018.05.015). 12. 12.d’Ambrosio A, Valsasina P, Gallo A, et al. Reduced dynamics of functional connectivity and cognitive impairment in multiple sclerosis. Mult Scler. England; Epub 2019 Mar.:1352458519837707. 13. 13.Eijlers AJC, Wink AM, Meijer KA, Douw L, Geurts JJG, Schoonheim MM. Reduced Network Dynamics on Functional MRI Signals Cognitive Impairment in Multiple Sclerosis. Radiology. United States; 2019;292:449–457. 14. 14.van Geest Q, Hulst HE, Meijer KA, Hoyng L, Geurts JJG, Douw L. The importance of hippocampal dynamic connectivity in explaining memory function in multiple sclerosis. Brain Behav. 2018;8:1–12. 15. 15.Rocca MA, Hidalgo de La Cruz M, Valsasina P, et al. Two-year dynamic functional network connectivity in clinically isolated syndrome. Mult Scler J. Epub 2019.:1–14. 16. 16.Tijhuis FB, Broeders TAA, Santos FAN, et al. Dynamic functional connectivity as a neural correlate of fatigue in multiple sclerosis. NeuroImage Clin [online serial]. Elsevier Inc.; 2021;29:102556. Accessed at: [https://doi.org/10.1016/j.nicl.2020.102556](https://doi.org/10.1016/j.nicl.2020.102556). 17. 17.Fu Z, Caprihan A, Chen J, et al. Altered static and dynamic functional network connectivity in Alzheimer’s disease and subcortical ischemic vascular disease: shared and specific brain connectivity abnormalities. Hum Brain Mapp. 2019;40:3203–3221. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.24591&link_type=DOI) 18. 18.Tuleasca C, Bolton TAW, Régis J, et al. Normalization of aberrant pretherapeutic dynamic functional connectivity of extrastriate visual system in patients who underwent thalamotomy with stereotactic radiosurgery for essential tremor: A resting-state functional MRI study. J Neurosurg. 2020;132:1792–1801. 19. 19.Zöller D, Sandini C, Karahanoğlu FI, et al. Large-Scale Brain Network Dynamics Provide a Measure of Psychosis and Anxiety in 22q11.2 Deletion Syndrome. Biol Psychiatry Cogn Neurosci Neuroimaging. 2019;4:881–892. 20. 20.Zhuang X, Walsh RR, Sreenivasan K, Yang Z, Mishra V, Cordes D. Incorporating spatial constraint in co-activation pattern analysis to explore the dynamics of resting-state networks: An application to Parkinson’s disease. Neuroimage [online serial]. Elsevier Ltd; 2018;172:64–84. Accessed at: [https://doi.org/10.1016/j.neuroimage.2018.01.019](https://doi.org/10.1016/j.neuroimage.2018.01.019). 21. 21.Parmenter BA, Testa SM, Schretlen DJ, Weinstock-Guttman B, Benedict RHB. The utility of regression-based norms in interpreting the minimal assessment of cognitive function in multiple sclerosis (MACFIMS). J Int Neuropsychol Soc. Ecole Polytechnique Fédérale de Lausanne; 2010;16:6–16. 22. 22.Dusankova JB, Kalincik T, Havrdova E, Benedict RHB. Cross cultural validation of the minimal assessment of cognitive function in multiple sclerosis (MACFIMS) and the brief international cognitive assessment for multiple sclerosis (BICAMS). Clin Neuropsychol. 2012;26:1186–1200. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/13854046.2012.725101&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23034066&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F28%2F2020.11.26.20238923.atom) 23. 23.Battaglini M, Jenkinson M, De Stefano N. Evaluating and reducing the impact of white matter lesions on brain volume measurements. Hum Brain Mapp. 2012;33:2062–2071. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.21344&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21882300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F28%2F2020.11.26.20238923.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307593400004&link_type=ISI) 24. 24.Smith SM, Zhang Y, Jenkinson M, et al. Accurate, robust, and automated longitudinal and cross-sectional brain change analysis. Neuroimage. United States; 2002;17:479–489. 25. 25.Chao-Gan Y, Yu-Feng Z. DPARSF: A MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front Syst Neurosci. 2010;4:1–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnsys.2010.00145&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20204156&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F28%2F2020.11.26.20238923.atom) 26. 26.Aleman-Gomez Y, Melie-García L, Valdés-Hernandez P. IBASPM: toolbox for automatic parcellation of brain structures. 12th Annu Meet Organ Hum Brain Mapp. 2006. p. 11–15. 27. 27.Karahanoglu FI, Van De Ville D. Transient brain activity disentangles fMRI resting-state dynamics in terms of spatially and temporally overlapping networks. Nat Commun. England; 2015;6:7751. 28. 28.Benjamini Y, Hochberg Y, Benjamini Yoav HY. Benjamini and Y FDR.pdf [online]. J. R. Stat. Soc. Ser. B 1995. p. 289–300. Accessed at: [http://www.stat.purdue.edu/~doerge/BIOINFORM.D/FALL06/BenjaminiandYFDR.pdf%5Cn](http://www.stat.purdue.edu/~doerge/BIOINFORM.D/FALL06/BenjaminiandYFDR.pdf) [http://engr.case.edu/ray\_soumya/mlrg/controlling\_fdr\_benjamini95.pdf](http://engr.case.edu/ray_soumya/mlrg/controlling_fdr_benjamini95.pdf). 29. 29.Lurie DJ, Kessler D, Bassett DS, et al. Questions and controversies in the study of time-varying functional connectivity in resting fMRI. Netw Neurosci. 2020;4:30–69. 30. 30.Mormino EC, Smiljic A, Hayenga AO, et al. Relationships between beta-amyloid and functional connectivity in different components of the default mode network in aging. Cereb Cortex. 2011;21:2399–2407. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhr025&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21383234&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F28%2F2020.11.26.20238923.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000294808800020&link_type=ISI) 31. 31.Camp SJ, Stevenson VL, Thompson AJ, et al. A longitudinal study of cognition in primary progressive multiple sclerosis. Brain. 2005;128:2891–2898. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awh602&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16049040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F28%2F2020.11.26.20238923.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233667500015&link_type=ISI) 32. 32.Cole PhD JH, Raffel MD J, Friede PhD T, et al. Longitudinal Assessment of Multiple Sclerosis with the Brain-Age Paradigm. Ann Neurol. 2020;88:93–105. 33. 33.Hawellek DJ, Hipp JF, Lewis CM, Corbetta M, Engel AK. Increased functional connectivity indicates the severity of cognitive impairment in multiple sclerosis. Proc Natl Acad Sci U S A. 2011;108:19066–19071. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA4LzQ3LzE5MDY2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDQvMjgvMjAyMC4xMS4yNi4yMDIzODkyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 34. 34.Silva PHR, Spedo CT, Barreira AA, Leoni RF. Symbol Digit Modalities Test adaptation for Magnetic Resonance Imaging environment: A systematic review and meta-analysis. Mult Scler Relat Disord [online serial]. Elsevier B.V.; 2018;20:136–143. Accessed at: [https://doi.org/10.1016/j.msard.2018.01.014](https://doi.org/10.1016/j.msard.2018.01.014). 35. 35.Gonzalez-Castillo J, Bandettini PA. Task-based dynamic functional connectivity: Recent findings and open questions. Neuroimage [online serial]. Elsevier Ltd; 2018;180:526–533. Accessed at: [https://doi.org/10.1016/j.neuroimage.2017.08.006](https://doi.org/10.1016/j.neuroimage.2017.08.006).