|Home | About | Journals | Submit | Contact Us | Français|
We aimed to investigate functional connectivity and variability across multiple frequency bands in brain networks underlying cognitive deficits in primary-progressive multiple sclerosis (PP-MS) and to explore how they are affected by the presence of cortical lesions (CLs). We analyzed functional connectivity and variability (measured as the standard deviation of BOLD signal amplitude) in resting state networks (RSNs) associated with cognitive deficits in different frequency bands in 25 PP-MS patients (12M, mean age 50.9±10.5 years) and 20 healthy subjects (9M, mean age 51.0±9.8 years). We confirmed the presence of a widespread cognitive deterioration in PP-MS patients, with main involvement of visuo-spatial and executive domains. Cognitively impaired patients showed increased variability, reduced synchronicity between networks involved in the control of cognitive macro-domains and hyper-synchronicity limited to the connections between networks functionally more segregated. CL volume was higher in patients with cognitive impairment and was correlated with functional connectivity and variability. We demonstrate, for the first time, that a functional reorganization characterized by hypo-synchronicity of functionally-related/hyper-synchronicity of functionally-segregated large scale networks and an abnormal pattern of neural activity underlie cognitive dysfunction in PP-MS, and that CLs possibly play a role in variability and functional connectivity abnormalities.
Cognitive impairment occurs in approximately 47% of patients with primary-progressive multiple sclerosis (PP-MS)1, affecting a wide range of functions, from processing speed to executive functions and visuo-spatial memory1,2,3,4,5. Unlike the most common relapsing-remitting (RR) form of the disease, PP-MS is characterized by absence of relapses prior to clinical deterioration and by spinal cord focal and diffuse damage in presence of only few macroscopic inflammatory brain lesions6. This raises the question of what pathophysiological mechanisms underpin the development of cognitive impairment in PP-MS7. A recent study has suggested that it is related to the disruption of brain functional connectivity of the default mode network (DMN) as measured by resting-state functional MRI (rs-fMRI) and associated with microstructural abnormalities of the corpus callosum and cingulum8. However, while the relationship between cognitive-related functional connectivity abnormalities and the focal and diffuse damage of the connecting white matter (WM) tracts is well established in early and RR-MS9,10, it is not clear whether this holds true in PP-MS patients who are characterized by a relatively moderate WM lesion burden6 and a high burden of cortical lesions (CLs)11. Therefore, we postulated that, in addition to normal-appearing WM brain tissue injury, the presence of CLs might play a role in the modulation of functional connectivity in PP-MS patients.
Recent advances in rs-fMRI analysis have shown that the investigation of functional connectivity in multiple frequency bands and the assessment of blood-oxygen-level dependent (BOLD) signal amplitude variability could provide a more comprehensive characterization of brain functional correlates of cognitive impairment. Indeed, while functional connectivity of slow frequency bands (Slow-5: 0.01–0.027Hz and Slow-4: 0.027–0.073Hz)12,13 may contribute to the characterization of deep gray matter (GM) and midline cortical regions14,15,16, variability of BOLD signal amplitude across time might represent a direct index of neuronal activity. Several studies have shown that GM exhibits a higher BOLD signal amplitude than WM14,17 and that actively firing, highly connected neurons show a high spiking variability18. Thus suggesting that BOLD signal amplitude variation may reflect the variable metabolic demand of active neurons. Therefore, our aims were to: (i) identify the major domains of cognitive impairment in PP-MS through a principal component analysis (PCA) of neuropsychological tests and select resting state networks (RSNs) of interest accordingly; (ii) characterize, for the selected networks, functional connectivity and variability abnormalities associated with cognitive impairment in different frequency bands and, (iii) investigate whether the presence of CLs, which show a high prevalence in PP-MS patients11,19, has an impact on specific functional connectivity and variability changes.
Twenty-five patients with PP-MS (mean age±SD: 51.0±9.8, range: 34–63 years) were consecutively enrolled from January 1, 2012 through December 31, 2013. Inclusion criteria were: (a) age between 25 and 65 years; (b) diagnosis of clinically definite MS, according to the revised McDonald criteria20 and a PP course; (c) Expanded Disability Status Scale (EDSS) score no greater than 6.5 at screening; (d) if treated, patients had to be on stable treatment for at least 1 year. Twenty gender- and age-matched healthy volunteers (mean age±SD: 50.9±10.5, range: 32–65 years) were recruited as controls (CTRLs). For all participants, the following exclusion criteria were applied: (a) current or past history of major hematological, renal, hepatic, psychiatric disorder or neurologic disease other than MS (for the patient group); (b) contraindications to MRI. On the same day of MRI acquisition patients underwent neurological evaluation with EDSS rating.
The institutional review board of Icahn School of Medicine at Mount Sinai approved the study and all participants gave written informed consent before investigation according to the Declaration of Helsinki. All methods were performed in accordance with the relevant guidelines and regulations.
Neuropsychological evaluation was performed on the same day of the MRI session, administering the minimal assessment of cognitive function in MS (MACFIMS)21. Raw scores from MACFIMS battery were converted in z-scores22. The MACFIMS battery includes tests exploring sustained attention, concentration, processing speed and working memory (paced auditory serial addition test at 3s –PASAT; symbol digit modalities test -SDMT), verbal and visuo-spatial episodic memory (California verbal learning test II-CVLT-II; brief visuo-spatial memory test revised - BVMT-R), executive functions (Delis-Kaplan executive function scale sorting test -DKEFS; controlled oral word association test -COWAT) and spatial processing (judgment of line orientation test-JLO). Six CTRLs were non-native English speakers and did not perform MACFIMS evaluation. Raw scores were converted in z-scores. A PCA of the neuropsychogical tests (including both patients and controls) was performed using SPSS (version 20.0, IBM, Chicago, Illinois) and the first component’s eigenvalues were used to calculate a weighted average of the z-scores of the neuropsychological tests to obtain a global cognitive (GC) score. Patients were divided into cognitively impaired and cognitively preserved using a hierarchical clustering analysis of the GC score using Ward’s method with squared euclidian distance.
Participants underwent MRI on a 3.0T scanner (Achieva, Philips, The Netherlands). The MRI protocol included: (a) axial dual echo (DE) turbo spin echo (TSE) sequence: repetition time (TR)=2500ms, echo time (TE1)=10ms, TE2=80ms, field of view (FoV)=230×230mm2, matrix size=512×512, 46 contiguous 3mm-thick slices; (b) sagittal 3D T1-weighted turbo field echo (TFE) sequence: TR=7.5ms, TE=3.5ms, inversion time (TI)=900ms, flip angle=8°, voxel size=1×1×1mm3, 172 contiguous slices; (c) axial phase-sensitive inversion recovery (PSIR) sequence: TR=4500ms, TE=8ms, TI=400ms, 46 contiguous 3mm thick slices with in-plane reconstructed resolution of 0.5×0.5mm2; (d) T2* echo-planar imaging (EPI) for resting state fMRI: TR=2300ms, TE=27ms, 150 volumes; voxel size=3×3×3mm3. During resting state fMRI acquisition, subjects were asked to rest with their eyes closed without falling asleep; (f) a single-shot EPI sequence for diffusion tensor imaging: TR=8550ms, TE=89.5ms, voxel size=2×2×2mm3; b-values=0, 1000, 2000s/mm2 and 32 directions each.
Brain WM T2 and T1 lesion volumes (LV) were measured as previously described in Ghassemi and colleagues23. CLs were identified and segmented by a single observer blinded to clinical information and subject’s identity, under the supervision of a senior investigator, according to published consensus criteria24 on PSIR images (Jim 6; Xinapse Systems, England), and classified as intracortical (IC), leucocortical (LC) or juxtacortical (JC)24. JC lesions, involving only WM, were excluded from further analysis. T2 and CL probability maps were created by creating a mean image of lesion maps (sum of binarized individual maps divided by number of patients) after registration to common spatial template (see below).
Maps of GM z-scores used for subsequent fMRI analysis were calculated as described in Villain and colleagues25. Briefly, the T1-weighted images of all study subjects were segmented into GM, WM and cerebro-spinal fluid (CSF) maps with FSL VBM (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FSLVBM). The GM maps were then registered to the GM population-specific template generated from the complete image set masked with the GM population template (thresholded in order to retain only pixels having probability >30% to belong to the GM). The resulting individual GM maps were then masked with the thresholded GM population template and used to create z- score maps of GM atrophy by applying the following formula: ([subject individual GM value -control mean GM value]/control [SD] GM value). Individual z- score maps of GM atrophy and conservative GM masks, created for subsequent variability analysis by thresholding the individual GM probability maps at 70%, were finally normalized to the MNI space.
Tract-based Spatial Statistics (TBSS) was performed to compare WM microstructural damage between cognitively preserved and impaired patients. The b0 maps were used to extract the brain by using BET from FSL (www.fmrib.ox.ac.uk/fsl/). The diffusion images were co-registered and the Diffusional Kurtosis Estimator (DKE) software26 was used to obtain fractional anisotropy (FA) maps. Then, using the TBSS pipeline, individual subjects’ FA maps were registered to a common space (FMRIB58_FA standard-space) using non-linear registration with FNIRT. A mean FA skeleton was created representing the centers of all tracts common to all the subjects and then each subject’s individual aligned FA map was projected onto the skeleton.
Functional preprocessing steps were implemented in Analysis of Functional NeuroImages (AFNI)27 (http://afni.nimh.nih.gov/afni) and FSL. First, removal of the first 3 volumes of each functional dataset and slice-timing correction were performed using AFNI, followed by spatial smoothing (Gaussian kernel, FWHM=6mm) and removal of motion artifacts with FSL and ICA-AROMA28,29. In AFNI, T1-weighted and functional images for all subjects were then were normalized to the 2×2×2mm3 MNI template. The mean time-series from the WM and the CSF were regressed out. The data were then bandpass-filtered between 0.01 and 0.1Hz, the standard frequency band (SFB) thought to reflect mainly neuronal fluctuations30. Based on recent findings in healthy subjects12,14, we also focused on two separate frequency sub-bands within the SFB: Slow-5 (0.01–0.027Hz) and Slow-4 (0.027–0.073Hz).
A seed-based analysis of resting state functional connectivity was conducted using AFNI. Spheres of 6mm radius were used as seeds in regions of interest (ROI) in the dorsal attention network (DAN), right attentional network (RAN) and executive control network (ECN) according to published data31. The DAN, RAN and ECN were chosen as networks of interest corresponding to the cognitive domains found to be relevant in the PCA of the neuropsychological tests. Since the ECN, DAN and RAN show a better spatial differentiation in the frontal than in the parietal regions31,32, we decided to focus on the frontal and temporal areas in order to avoid the selection of partially overlapping seed ROIs. For the DAN, two seed-ROIs were selected in the left and right frontal eye fields (FEF) (MNI x, y, z coordinates −29,−5, 55, and 31,−5, 54); for the ECN, three seed-ROIs were used in the dorso-medial, left anterior and right anterior prefrontal cortex (dmPFC, laPFC, raPFC) (MNI x, y, z coordinates 1, 30, 44; −45, 50,−5; 46, 51,−7); for the RAN, two seed-ROIs were chosen in the right middle frontal gyrus and right middle temporal gyrus (rMFG, rMTG) (MNI x, y, z coordinates 34, 24, 44; 64, −39, −11). Mean time-series for each seed-ROI was calculated by averaging the time-series of all voxels within the sphere. Correlation analysis was performed between each seed-ROI and all other brain voxels for each frequency band (SFB, Slow-5 and Slow-4). Fisher’s r-to-z transform was used to obtain spatial z-maps with each voxel’s intensity representing its correlation strength with the seed-ROI.
Variability of the BOLD signal was measured by calculating the standard deviation of the time-series amplitude within the seed-ROIs, in each frequency band. Whole-brain variability was investigated by calculating the standard deviation within the entire GM15, using the conservative GM masks.
Results were anatomically labeled using the Talairach-Tournoux Atlas available in AFNI. Because head motion can affect BOLD time-series amplitude and is often higher in patient populations, we compared head motion (mean relative frame-to-frame displacement-FD calculated using FSL) between patients and CTRLs.
Between-group differences of gender and treatment were assessed with Fisher’s exact Test and Chi-square Test, whereas differences in other demographic, clinical and head motion parameters were tested with Mann-Whitney U Test and Kruskal-Wallis Test using the SPSS software package (version 20.0, IBM, Chicago, Illinois). Between-group differences in MRI structural parameters were tested with a multivariate analysis of variance, entering age as a covariate. Statistical significance was set at p<0.006, Bonferroni corrected for multiple comparisons.
TBSS voxel-wise non-parametric statistics were performed to compare FA spatial maps between cognitively impaired and cognitively preserved groups, using randomise with 5000 permutations, entering age as a covariate of no interest. Statistical significance was set at p<0.05, Threshold-Free Cluster Enhancement corrected for multiple comparisons.
To compare functional connectivity between patients and CTRLs, voxel-wise statistics were run in AFNI, using a two-sample t-test in each frequency band. To exclude a possible influence of regional GM atrophy on the voxel-wise comparisons of resting state activity between CTRLs and PP-MS patients, the z-scores of GM volume were used as a confounding voxel-wise covariate in the statistical analysis as described previously25. Age was also used as covariate to control for potential aging effects, as well as FD, in order to take into account the effects of head motion. In patients, voxel-wise correlations between CL/IC lesions volume and the functional connectivity maps showing significant differences in the between group comparisons were performed in AFNI, entering age, GM z-scores and FD as covariates. For both analyses, the resulting t-maps were thresholded at a corrected p-value of 0.01 (multiple-comparisons corrected using the Monte Carlo simulation as implemented in Alphasim, AFNI).
Mean variability values extracted from GM and seed regions were compared in SPSS, using an ANOVA corrected for age, GM volume and FD, Bonferroni corrected at p<0.002.
In patients, voxel-wise correlations between CL/IC lesions volume and variability in each frequency band were performed in AFNI, entering age, GM z-scores and FD as covariates and thresholding the resulting t-maps at a corrected p value of 0.01 (multiple-comparisons corrected using the Monte Carlo simulation as implemented in Alphasim, AFNI).
In patients, voxel-wise correlations between functional connectivity and variability in GM were performed in AFNI, entering age, GM z-scores and FD as covariates. All resulting t-maps were thresholded at a corrected p value of 0.01 (multiple-comparisons corrected using the Monte Carlo simulation as implemented in Alphasim, AFNI).
PP-MS patients showed a wide range of cognitive deficits (Fig. 1a). PCA of the neuropsychological tests yielded a first principal component that explained 40.26% of the variance. The main tests contributing to this component were the BVMT-R total learning (weight 0.828) and the COWAT (weight 0.727) (Fig. 1b). Therefore, our subsequent functional analysis was focused on the DAN, RAN and ECN. The DAN, involved in spatial orientation of attention33 and visuospatial working memory34,35,36, is composed of the intraparietal sulcus and the junction of the precentral and superior frontal sulcus (frontal eye field, FEF) in each hemisphere. The RAN, engaged in attention reorientation37 and visuospatial processing38, is composed of the temporal-parietal junction and ventral frontal cortex33. The ECN facilitates higher cognitive processes including flexibility of behavioral choices in response to shifting conditions and context, and includes known sites for sustained attention/working memory (dorsolateral-prefrontal cortex, lateral parietal cortex)39 and response selection (dorsomedial-frontal cortex)40.
The first component’s eigenvectors were used to calculate a weighted average of the z-scores of the neuropsychological tests to obtain a GC score. Based on the clustering analysis, 13 patients were classified as cognitively impaired and 12 patients as cognitively preserved. Significant differences in GC score were identified between cognitively impaired and preserved patients (−8.94±2.00 vs −1.90±1.88, p=0.000022) and between cognitively impaired patients and HC (−8.94±2.00 vs −1.80±4.09, p=0.000046) (Fig. 1c).
Controls, cognitively preserved, and cognitively impaired MS patients were matched for gender, and level of education. While no differences in terms of disease modifying treatment, disease duration, EDSS and WM lesion loads were detected between cognitively preserved and impaired patients, cognitively preserved patients resulted significantly older (p=0.030). Cognitively impaired patients showed a significantly higher IC lesion volume (respectively p=0.004). Comparison of cognitively preserved and impaired patients WM microstructural damage showed a widespread significant decrease in FA in the cognitively impaired group (Supplementary Fig. 2a). WM and GM lesion distribution is shown in Supplementary Fig. 2b.
No significant differences were detected between PP-MS patients and CTRLs (0.25±0.20 vs 0.19±0.10, p=0.479) and between cognitively preserved and impaired patients in terms of FD (0.29±0.25 vs 0.21±0.11, p=0.894).
The regional differences in functional connectivity between CTRLs and PP-MS patients for the DAN, RAN and the ECN seed regions, in the three frequency bands are reported in detail in Supplementary Table 1 (p=0.01, Monte Carlo corrected for multiple-comparisons).
In the SFB (Fig. 2, red panel), compared to CTRLs, PP-MS patients showed increased connectivity between the frontal DAN and regions of the DMN and between the fronto-parietal regions of the ECN and frontal areas of the RAN. Reduced connectivity was detected between the frontal ECN, cerebellum and primary somatosensory cortex.
In Slow-5 (Fig. 2, blue panel), PP-MS patients showed increased connectivity between: the frontal DAN and temporal areas of the RAN; the frontal ECN and the DMN; the frontal RAN and the frontal region of the left attentional network (LAN) and the parietal ECN component. Decreased connectivity was found between the temporal RAN component and the supplementary motor area.
In Slow-4 (Fig. 2, green panel), patients showed decreased connectivity between the temporal RAN region and the insular node of the salience network-SN. In both SFB and Slow-4 within network disconnection was present in the RAN (Fig. 2).
In the SFB (Fig. 3, red panel), patients showed decreased connectivity between the frontal ECN components, the insula and parietal RAN regions. Reduced connectivity was also present between the frontal DAN and the temporal RAN component.
In Slow-5 (Fig. 3, blue panel), patients showed increased connectivity of both the fronto-temporal RAN and the fronto-medial ECN components with the DMN. They also showed decreased connectivity between the frontal DAN and parietal LAN areas, between the frontal ECN and parietal RAN areas, and between the temporal RAN and the supplementary motor area.
In both SFB and Slow-4 (Fig. 3, green panel), patients showed increased connectivity between the frontal DAN and the cerebellum, between the frontal DAN and the cingulate gyrus, between the frontal ECN and the RAN, and within temporal components of the RAN. They also showed decreased connectivity between the frontal ECN and the primary somatosensory cortex.
Across all frequency bands, decreased connectivity was detected between the frontal ECN components and the cerebellum.
Compared to CTRLs, in the SFB (Fig. 4, red panel), patients showed increased connectivity between the temporal RAN and the frontal LAN components and within the ECN. Decreased connectivity was present between the temporal RAN region and the DMN.
In Slow-5 (Fig. 4, blue panel), patients showed decreased connectivity between the frontal components of DAN and the fronto-medial DMN, between the frontal components of ECN and the primary motor cortex and between the temporal components of RAN and the cingulate components of SN. They also showed increased connectivity between the frontal components of ECN and the fronto-medial DMN and between temporal areas of RAN, the frontal component of DMN and parietal component of LAN.
In Slow-4 (Fig. 4, green panel) patients showed decreased connectivity between the temporal component of RAN and occipital areas and within temporal components of RAN and increased connectivity between the frontal components of DAN and ECN.
In both the SFB and Slow-4, patients showed increased FC between the frontal components of the DAN and RAN.
The direct comparison of cognitively impaired and preserved patients showed, in the SFB (Fig. 5, red panel) decreased connectivity between the frontal DAN component, the insula and occipital areas.
In Slow-5 (Fig. 5, blue panel), cognitively impaired patients showed decreased FC between frontal ECN components and the temporal RAN, cingulate and occipital cortex.
In Slow-4 (Fig. 5, green panel), cognitively impaired patients showed decreased connectivity between the frontal DAN, the fronto-temporal DMN and temporal RAN regions, between frontal ECN and the cerebellum, and between the temporal and parietal regions of the RAN. They also showed increased connectivity between frontal DAN components and the cingulate cortex, and between the temporal RAN and the cerebellum.
In the SFB and Slow-5, cognitively impaired patients showed decreased connectivity within temporal and frontal components of the RAN, while in the SFB and Slow-4 increased connectivity was present within the temporal regions of the RAN.
Across all frequency bands cognitively impaired patients showed increased connectivity between the frontal components of DAN and the parietal component of RAN.
Between groups comparison of variability in GM (global variability) and seed ROIs (regional variability) showed constantly higher variability in patients than CTRLs and in cognitively impaired patients in comparison with cognitively preserved patients. Global variability was higher in PP-MS patients than CTRLs (SFB 0.43 vs 0.37; Slow-5 0.25 vs 0.20; Slow-4 0.30 vs 0.25, p<0.0001) and in cognitively impaired in comparison with cognitively preserved (SFB 0.46 vs 0.41; Slow-5 0.27 vs 0.23; Slow-4 0.31 vs 0.28, p<0.001). Significant results of regional variability analysis are reported in Table 2.
Significant correlations with connectivity were found in several regions in all frequency bands for both CL and IC lesion volume (detailed results in Supplementary Table 2). Significant correlations with variability were found for CL volumes in the cingulate gyrus and superior frontal gyrus (detailed results in Supplementary Table 3) (p 0.01, Monte Carlo corrected for multiple-comparisons).
Significant correlations between connectivity and global variability was found in several regions in all frequency bands (p 0.01, Monte Carlo corrected for multiple-comparisons). The detailed results are reported in Supplementary Table 4.
The functional connectivity analysis of the DAN, RAN and ECN confirmed the presence of a widespread cortical reorganization, involving the connection between areas pertaining to these networks and regions involved in other networks also responsible for cognitive control (i.e. LAN, DMN and SN). Such functional reorganization mainly resulted in increased between-network connectivity, with hyper-synchronization of frontal attentional networks (RAN-LAN), frontal and dorsal attentional systems (RAN-DAN), executive and attentional networks (ECN-RAN) and connections to the DMN (DAN-DMN, ECN-DMN). Disconnection resulted limited to RAN within network connectivity, connections between RAN and SN as well as connections between ECN/RAN and areas involved in motor control.
The specificity of our results in different frequency bands confirms the higher sensitivity of the lowest bands to fluctuations occurring along the midline structures (Slow-5)15,16, thus supporting the concept that resting-state BOLD oscillations exhibit a frequency-dependent, anatomically constrained spatial structure13.
When the patients were stratified on the basis of cognitive impairment and compared with CTRLs, our results disclosed the presence of a large-scale reconfiguration of brain networks interaction, with different characteristics in cognitively impaired and preserved patients. In cognitively preserved patients increased synchronicity was present between attentional networks (DAN, RAN, LAN), and within ECN. In cognitively impaired patients the connections between networks involved in the control of macro-domains (DAN, RAN and LAN for attention; ECN, SN, motor areas for stimulus selection, planning and integration) showed reduced synchronicity, with hyper-synchronicity being limited to the connections between networks functionally more segregated (ECN-RAN) and between attentional networks and cerebellum.
Functional abnormalities between large-scale neuronal networks and within DMN sub-regions have been described in RR-MS41,42 with cognitive impairment. In one study, RR-MS patients showed increased connectivity between the ECN and SN and decreased connectivity between the ECN and DMN when compared to healthy CTRLs42. In a second study, RR-MS patients without cognitive impairment showed a weaker functional connectivity in the anterior cingulate cortex and in the core of the posterior cingulate cortex (PCC), and a stronger connectivity at the periphery of the PCC41. In PP-MS patients, rather than an intra-network adaptation, we found a differential inter-network reorganization of DMN regions. Specifically, cognitively preserved patients showed increased synchronicity between both RAN and ECN and the frontal component of DMN and reduced synchronicity between temporal RAN and posterior cingulate cortex, as well as between frontal areas of DAN and DMN while cognitively impaired patients showed increased synchronicity between fronto-temporal areas of the attentional network (RAN) and fronto-temporal DMN, between frontal areas of DAN and cingulate cortex as well as between frontal ECN and frontal DMN.
Although we cannot infer from our results the beneficial vs maladaptive role of such functional reorganization43, the loss, in cognitively impaired patients, of the hyper-synchronicity between RSN involved in strictly related functions and the differential reorganization between these networks and the DMN, whose central role in cognitive systems integration is well known44, seems to suggest that these rearrangements could be responsible for the relative functional preservation observed in cognitively preserved patients.
Interestingly, in addition to the different functional connectivity reorganization, cognitively impaired patients showed higher prevalence of IC lesions and more widespread WM injury when compared with cognitively preserved patients. The disruption of physiological connections determined by diffuse WM damage and focal GM lesions in PP-MS could result in the loss of some connections and in the amplification of the connections still structurally preserved, without necessarily determining any functional advantage in terms of cognitive ability. Since in our patients the functional connectivity rearrangement was not associated with the maintenance of an adequate cognitive status, we believe it may suggest the presence of maladaptive reorganization, possibly related to the extent of WM structural damage and the presence of CLs45. Indeed, we found significant correlations between the number of CLs and functional connectivity, suggesting that the presence of lesions in the brain cortex might indeed have a role in the functional reorganization of brain networks.
Perhaps, the most novel aspect of our study is related to the variability analysis that revealed increased values in both the whole group of patients and in the two subgroups when compared with CTRLs, and in cognitively impaired when compared with cognitively preserved patients. Variability is considered an index of neural activity46,47. It increases with brain maturation47 and decreases with age showing a bidirectional pattern46. The existence of a spatial pattern, and a fluid transition of variability status across life stages seem to be a critical index of the cognitive capacity46. Although it remains speculative, the increased variability in PP-MS could reflect a dysfunctional neural activity responsible for the disruption of networks interactions. In particular, previous research suggests that there is an optimal level of variability that facilitates neural function46,48; it is therefore possible that a homogeneous increase in neuronal variability would affect the functional integration of brain networks. In alternative, the increased variability could represent a failed attempt to compensate for GM lesional damage, as suggested by the presence of a direct correlation between CLs volume and variability. In this scenario, the abnormal network integration would elicit an enhanced neuronal activity in structurally preserved GM areas.
Theoretically, variability should be independent of functional connectivity49, but recent studies have reported a close relationship between them50,51. While in healthy subjects higher variability correlates with stronger functional connectivity51, in presence of brain pathology this correlation is lost15. In our population, the two metrics showed a widespread voxelwise correlation. In patients, variability showed both direct and inverse correlations with functional connectivity in different brain areas, suggesting that a maladaptive increase in neural activity might be associated with the reorganization of functional connectivity. Our study is not without limitations. The small sample size may have limited our ability in determining significant differences in the subgroup comparisons. Most of our patients were under treatment with FDA approved disease-modifying agents; while we cannot rule out a potential treatment effect on measures of brain functional connectivity and variability or on the expression of cognitive deficits52, the lack of difference in type and frequency of treatment between the cognitively preserved and impaired patients should have minimized any confounding effects. Longitudinal studies on larger cohorts are needed to fully elucidate the role of the brain functional reorganization underlying the cognitive deficits in PP-MS and to differentiate the maladaptive from the compensatory modifications.
In conclusion, our results show for the first time that cognitive impairment in PP-MS patients is associated with increased variability and hypo-synchronicity of functionally-related/hyper-synchronicity of functionally-segregated large scale networks, which are, in turn, associated with CLs. While the variability analysis provides information on the neural activity underlying the functional connectivity abnormalities, the additional investigation of sub-frequency bands allows wider detection of functional abnormalities, with increased sensitivity along the midline structures. Hence, our novel findings contribute to improving characterization of the resting state correlates of cognitive impairment in PP-MS, and, once replicated in a larger cohort, may eventually be of use as a biomarker in clinical assessment and treatment.
How to cite this article: Petracca, M. et al. Synchronization and variability imbalance underlie cognitive impairment in primary-progressive multiple sclerosis. Sci. Rep. 7, 46411; doi: 10.1038/srep46411 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors would like to thank S. Krieger, MD, Department of Neurology, Icahn School of Medicine at Mount Sinai, C. Riley, MD, The Neurological Institute of New York, Department of Neurology, Columbia University Medical Center and J.Howard, MD, Department of Neurology, Langone Medical Center, New York University for patients’ referral; Y. Bencosme, CCRC, Department of Neurology, Icahn School of Medicine at Mount Sinai for patients’ scheduling and assistance and Nichole Bozyk for technical support. This study was supported in part by Novartis Pharmaceuticals (CFTY20DUSNC15T), National Multiple Sclerosis Society (NMSS RG 5120A3/1) and the Noto Foundation to MI and by a research fellowship from Fondazione Italiana Sclerosi Multipla (2013/B/7) to MP.
Dr. Maria Petracca has nothing to disclose. Dr. Catarina Saiote has nothing to disclose. Dr. Bender has nothing to disclose. Ms. Arias has nothing to disclose. Ms. Colleen Farrell has nothing to disclose. Dr. Magioncalda has nothing to disclose. Dr. Martino has nothing to disclose. Dr. Aaron Miller has served as a consultant and/or participant in advisory board meetings for Genzyme/sanofi-aventis, Biogen Idec, Glaxo Smith Kline, EMD Serono (Merck Serono), Novartis, ONO, Acorda, Nuron Biotech, Teva, Questcor and Accordant Health Services. He has received research support from Acorda, Novartis, Genentech, Genzyme/sanofi-aventis, Biogen Idec, Roche, and Questcor. He has served as Editor of Continuum, a continuing medical education publication of the AAN and currently serves as Editor of Continuum Audio. He is a member of the editorial board of Multiple Sclerosis and Related Disorders. He occasionally performs expert reviews of medical records or serves as an expert witness in medical malpractice cases. Dr. Northoff has nothing to disclose. Dr. Fred Lublin Sources of Funding for Research: Acorda Therapeutics, Inc.; Biogen Idec; Novartis Pharmaceuticals Corp; Teva Neuroscience, Inc.; Genzyme; Sanofi; Celgene; NIH; NMSS; Consulting Agreements/Advisory Boards/DSMB: Bayer HealthCare Pharmaceuticals; Biogen Idec; EMD Serono, Inc.; Novartis; Teva Neuroscience; Actelion; Sanofi-Aventis; Acorda; Questcor; Roche, Genentech; Celgene; Johnson & Johnson; Revalesio; Coronado Bioscience, Genzyme, MedImmune; Bristol-Myers Squibb, Xenoport, Receptos; Forward Pharma; Co-Chief Editor: Multiple Sclerosis and Related Diseases; Stock Ownership: Cognition Pharmaceuticals, Inc. Dr. Matilde Inglese has received research grants from NIH, NMSS, Novartis Pharmaceuticals Corp., Teva Neuroscience. None of the competing financial interests have influenced the presented work in any stage from conception to analysis and publication.
Author Contributions M.I. had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: M.I., G.N. and F.L. Acquisition, analysis, or interpretation of data: M.P., C.S., H.A.B., F.A., C.F., P.M., M.M., A.M., G.N., F.L., M.I. Drafting of the manuscript: M.P., C.S., M.I. Critical revision of the manuscript for important intellectual content: M.P., C.S., H.A.B., F.A., C.F., P.M., M.M., A.M., G.N., F.L., M.I. Statistical analysis: M.P., C.S., M.I. Obtained funding: M.I. Study supervision: M.I.