Abstract
There is mounting evidence regarding the role of impairment in neuromodulatory networks for neurodegenerative diseases, such as Parkinson’s and Alzheimer’s disease. However, the role of neuromodulatory networks in multiple sclerosis (MS) has not been assessed. We applied resting-state functional connectivity and graph theory to investigate the changes in the functional connectivity within neuromodulatory networks including the serotonergic, noradrenergic, cholinergic, and dopaminergic systems in MS. Twenty-nine MS patients and twenty-four age- and gender-matched healthy controls performed clinical and cognitive assessments including the expanded disability status score, symbol digit modalities test, and Hamilton Depression rating scale. We demonstrated a diffuse reorganization of network topography (P < 0.01) in serotonergic, cholinergic, noradrenergic, and dopaminergic networks in patients with MS. Serotonergic, noradrenergic, and cholinergic network functional connectivity derangement was associated with disease duration, EDSS, and depressive symptoms (P < 0.01). Derangements in serotonergic, noradrenergic, cholinergic, and dopaminergic network impairment were associated with cognitive abilities (P < 0.01). Our results indicate that functional connectivity changes within neuromodulatory networks might be a useful tool in predicting disability burden over time, and could serve as a surrogate endpoint to assess efficacy for symptomatic treatments.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
Focal demyelination in multiple sclerosis (MS) leads to the disruption of the axon–myelin unit and sub-sequential neuronal damage [1]. Loss of pre-synaptic neurons causes reduced neurotransmitter release at the synaptic cleft resulting in disrupted functional connectivity among distant brain regions [2]. Post-mortem studies in MS reported a reduced synaptic density and impaired axonal transport, synaptic plasticity, glutamate neurotransmission, and glutamate homoeostasis [3]. Demyelination, axonal loss, and synaptic pathology cause a dysfunctional communication among brain regions. Neurons throughout the brain communicate through the release of neurotransmitters, subsequently modulating neuronal activity. Several studies have outlined neurotransmitter dysfunction in neurodegenerative disorders, such as Parkinson’s, Alzheimer’s, and Huntington’s disease, involving mainly the dopaminergic, serotonergic, and cholinergic networks. Network impairment in such neurodegenerative disorders has been selectively associated with neurological symptoms such as cognitive, motor, and non-motor symptoms [4,5,6,7,8]. Tools to measure impaired neurotransmitter networks integrity require either the application of PET or complex MRI methodology with contrast injection. For this reason, since MS is not primarily a neurotransmitter disorder, the precise involvement of neuromodulatory networks in MS has not been fully elucidated and existing evidence is scarce. Previous studies reported that a reduced glutamate and gamma-aminobutyric acid concentration throughout several brain areas in MS patients measured through magnetic resonance spectroscopy in the brain correlated with clinical disability [9,10,11].
Resting-state functional magnetic resonance imaging (fMRI) is able to explore connectivity in predefined networks as for example the default mode network, the visual network, and the motor network. In MS, patients could experience either a decrease or an increase in functional connectivity [12,13,14,15,16]. In addition, resting-state fMRI might also depict network topography changes, through graph theoretical analyses. Graph theory is a mathematical representation of the brain connections, consisting of a set of nodes (corresponding to cerebral regions) and edges (corresponding to functional or structural connection) interposed between them. The topography of the network can be described using different measures such as clustering coefficient, path length, and efficiency measures to characterize the ‘small-world’ properties of brain networks or degree and betweenness centrality to identify the crucial areas within the network [17].
Recent evidence suggests that motor and cognitive decline in neurodegenerative disorders, such as Parkinson’s and Alzheimer’s disease, selectively depends upon neurotransmission integrity of dopaminergic and cholinergic networks, respectively [18, 19]. We hypothesized that demyelination and axonal loss randomly occurring throughout the brain, in MS patients, might cause dysfunctional activity and, hence, impaired functional connectivity (FC) in neuromodulatory networks. Furthermore, clinical outcomes in MS, such as depressive symptoms or cognitive impairment, could depend on the integrity of selected neuromodulatory networks as has been observed in other neurological disorders, such as the selective loss of dopaminergic neurons in Parkinson’s disease leading to motor, or the loss of cholinergic neurons in Alzheimer’s disease leading to cognitive dysfunction [18, 19]. Therefore, we aimed at investigating FC and network topography changes within the serotonergic, noradrenergic, cholinergic, and dopaminergic networks in MS, and to explore their relationship with both physical and cognitive disability.
Materials and methods
Subjects
Twenty-nine patients with a definite diagnosis of relapsing–remitting MS according to the 2010 revised McDonalds criteria [20] and twenty-four healthy controls (HCs) were recruited from the specialist clinics at King’s College Hospital NHS Foundation Trust and through advertisement. All participants had no history of other neurological or psychiatric disorders, and were successfully screened to undertake MRI scanning according to the standard MRI scanning safety criteria (https://www.mrisafety.com). For each subject, clinical disability was assessed using the Expanded Disability Status Scale (EDSS). Multiple Sclerosis Severity Score (MSSS) was calculated as measure of disease severity [21]. Cognitive function was assessed with the symbol digit modalities test (SDMT). This test was chosen among other cognitive tests, because it showed a high sensibility and specificity in detecting cognitive impairment in MS [22]. The Hamilton depression rating scale (HRSD) was used to assess neuropsychiatric symptoms. The study was approved by the London-South East Research Ethics Committee. Written informed consent was obtained from all study’s subjects prior to inclusion in the study in accordance to the Declaration of Helsinki. Data will be made available via a request to the authors after a formal data sharing agreement.
MRI acquisition
MRI images were acquired using a 3-T Magnetom Biograph mMR PET/MR hybrid scanner (Siemens, Erlangen, Germany). MRI sequences included three-dimensional magnetization-prepared rapid gradient‐echo T1 (three-D-T1w MPRAGE, TR [repetition time]: 1700 ms, TE [echo time]: 2.63 ms, TI [inversion time]: 900 ms, flip angle: 9°, voxel size: 1.1 × 1.1 × 1.1 mm), turbo spin-echo T2-weighed images (T2w TSE, TR: 3200 ms, TE: 409 ms, TI: 900 ms, voxel size: 1 × 1 × 1mm); and fast fluid-attenuated inversion recovery (FLAIR, TR: 5000 ms, TE: 499 ms, TI: 1800 ms, voxel size: 1 × 1 × 1mm). In addition, resting-state fMRI scans were acquired using an echo-planar imaging sequence (TR: 3000 ms, TE: 30 ms, voxel size: 3 × 3 × 3mm, echo-planar imaging factor: 64, field of view: 192 mm, number of slices: 36; number of volumes: 240, flip angle: 90°).
MR image processing
For each subject, T1 and T2 lesions were manually defined through Analyze medical imaging software (version 12, Mayo Foundation AnalyzeDirect).
T1 lesions mask was used to perform a ‘lesion filling’ procedure [23] on the 3D-T1w MPRAGE. We applied SIENAX [24] to estimate cortical gray matter volume, normalized for subject head size. Multi-atlas propagation with enhanced registration approach (MAPER) [25] method was applied to create the subject-specific brain segmentation. This robust, automated technique segments subjects’ structural T1w MRI into 95 anatomical regions with high accuracy and is applicable even to subjects with significant cortical atrophy and/or ventriculomegaly.
Structural and fMRI data were further processed analyzed through SPM12 software (Statistical Parametric Mapping version 12; https://www.fil.ion.ucl.ac.uk/spm/software/spm12/) and the CONN toolbox version 17a [26]. Briefly, T1 lesion-filled scans and gray matter regions of interest map were normalized into the Montreal Neurological Institute (MNI) 152 space atlas. FMRI data were realigned and co-registered to the structural images, normalized to the MNI-152 space atlas, and then smoothed, using 6 mm full width half max isotropic Gaussian kernels band-pass filtering was performed with a frequency window of 0.008–0.09 Hz. Signal contributions from micro-head movements were corrected using the image ‘scrubbing’ ART method implemented in CONN toolbox. Outlier time points were included as covariates in denoising, and first-level general linear model along with motion parameters.
Networks’ definition
We performed a region-of-interest-based seed-to-seed connectivity analysis within four selected neuromodulatory networks: the serotonergic, the noradrenergic, the cholinergic, and the dopaminergic network. Network maps were created including anatomic regions of the central nervous system with afferent and efferent neurons using selected neuromodulatory molecules. To maximize the precision of the anatomical localization, we selected brain regions from subject-specific brain segmentation with the MAPER approach [25]. This approach allows brain segmentation in native space. When unavailable on the subject-specific brain segmentation, brain regions were selected from standard atlases to complete the network. For the detailed description of networks’ anatomy, please see Fig. 1 and Supplementary Appendix A-1. Briefly, the serotonergic network was created following Strac and colleagues, and Jacobs and Azmitia anatomical description [27, 28]. The cholinergic network was created following Selden and colleagues anatomical description, including four divisions: the Ch1-2–3 division, the Ch4 medial division, the Ch4 lateral perisylvian division, and the Ch4 lateral capsula division [29]. For the noradrenergic network, taking into account Samuels and colleagues anatomical description and highly overlapping regions within the serotonergic network (noradrenergic network includes the locus coeruleus, pedunculopontine nucleus, and ventral tegmental area in addition to the serotonergic network), we decided to merge the two networks creating a serotonergic–noradrenergic network to avoid ROI overlapping [30]. The dopaminergic network was created following Tziortzi and colleagues anatomical description, including three divisions: the dopaminergic executive division, the dopaminergic limbic division, and the dopaminergic motor division excluding the ventral tegmental area already included in the serotonergic–noradrenergic network to avoid ROI overlapping [31]. For each region of interest, we verified that the volume for each region was greater than the full width half max of the scanner resolution to avoid blurring effect. For each subject, we created a region of interest-to-region of interest serotonergic, cholinergic, noradrenergic, and dopaminergic functional connectivity matrix in MNI space.
Graph theory measures
For each network matrices, we performed a network topographic analysis through CONN functional connectivity toolbox [26]. We used a set of functional connectivity indices, namely the global efficiency, the local efficiency, the betweenness centrality, the average path length, the clustering coefficient, and the degree [17]. Global efficiency is a measure of network integration describing the information flow over the whole network. Local efficiency is the distance between the nodes and is a measure of network integration. Betweenness centrality is a measure of network centrality; it corresponds to the fraction of all the shortest paths in a network over a given node. Nodes with higher betweenness centrality participate in a large number of short paths. Average path length corresponds to the shortest distance between pairs of nodes, and is a measure of network efficiency of information transfer within the network. An increased average path length reflects the inability of the network to transfer information in parallel. Clustering coefficient is a measure of network organization reflecting the number of connections between the neighbor’s nodes. Degree is a measure of centrality and reflects the number of connections for each node to all other nodes in the network. A higher degree means a greater influence of the corresponding region in the network. To extract the best indices of network organization, network edges were adjusted with a threshold for a cost higher than 0.15 on a two-sided test for the adjacency matrix [32].
Statistics
Statistical analyses were performed through CONN functional connectivity toolbox using a seed-to-seed approach in SPM. Between-group analysis was performed using two-tailed t tests to compare the resting-state functional connectivity and graph theory measures change between MS patients and HCs for each network. Linear regression analysis was also performed to investigate the association between clinical disability and both the resting-state functional connectivity and graph theory measures for each neurotransmitter network for MS patients. Linear regression was performed only for those networks that showed a difference at the between-group analysis. Age and gender were included in the analysis as covariates for the between-group comparison, while age, gender, cortical gray matter volume and T2 lesions volume were used as covariate for the linear regression analyses with the EDSS, SDMT, and HRSD as dependent variables. Only age and gender were included as covariates for the linear regression analyses with disease duration and disease severity as dependent variable, considering the close interaction between these two latter variables and both lesion volume and cortical atrophy. Given the exploratory nature of the study, a threshold of P < 0.01 uncorrected was applied for between-groups and linear regression data analyses. Notwithstanding the exploratory nature of the study, we also applied the Benjamini–Hochberg procedure to control for the false discovery rate when evaluating correlation matrices in the FC analysis, setting the false discovery rate to 10% [33].
Results
Clinical assessments
MS patients performed worse than HCs in the SDMT (P < 0.001) and scored higher in HRSD (P < 0.001). Clinical, demographic, and imaging characteristics of the HCs and MS patients are summarized in Table 1.
Functional connectivity in serotonergic–noradrenergic network
MS patients showed decreased resting-state functional connectivity between different cerebellum regions (P = 0.006) and between the cerebellum and both the right amygdala and right thalamus (P < 0.01) compared with HCs. Patients also showed increased resting-state functional connectivity between the dorsal raphe and both the left and right superior frontal gyrus (P < 0.001), between the dorsal raphe and the left middle frontal gyrus, (P = 0.002), and between the left globus pallidus and the right hippocampus (P = 0.001) compared with HCs (Table 2 and Fig. 2).
Graph theory analysis revealed that MS patients showed increased betweenness centrality in the cerebellum (t-value = 2.84; P = 0.007) and a decreased betweenness centrality in the brainstem (t value = − 3.03; P = 0.004) for the serotonergic–noradrenergic network (Table 3 and Fig. 3).
Functional connectivity in cholinergic network
Ch 1-2-3 division
Regarding Ch1-2-3 division of the cholinergic network, MS patients showed no differences for both functional connectivity and graph theory evaluation compared to HCs.
Ch 4 medial division
Functional connectivity evaluation for the Ch4 medial division of the cholinergic network showed no difference between HCs and MS patients. Graph theory analysis revealed that MS patients have a higher local efficiency (t value = 0.36; P = 0.006) and a higher clustering coefficient (t-value = 0.34; P = 0.007) compared with HCs for the Ch4 medial division of the cholinergic network (Table 3 and Fig. 3).
Ch 4 lateral perisylvian division
MS patients have a similar functional connectivity for the Ch4 lateral perisylvian division of the cholinergic network compared with HCs, whereas they showed a decreased betweenness centrality in the right insula (t-value = -3.02; P = 0.004) and a decreased average path length in the left insula (t value = -2.94; P = 0.005) compared with HCs when evaluating the network topography through graph theory (Table 3 and Fig. 3).
Ch 4 lateral capsula division
Resting-state fMRI evaluation revealed that MS patients, compared with HCs, have reduced functional connectivity between the Ch4 nucleus and both the middle and superior frontal gyrus (all P < 0.006), and between the left posterior cingulate cortex and the left cuneus (t-value = − 2.91; P = 0.005). MS patients also showed an increased connectivity between the left middle temporal gyrus and the right middle temporal gyrus (t-value = − 3.3; P = 0.002), between the left middle temporal gyrus and both the left and right pre- and post-central gyrus (all P < 0.006), and between the right fusiform gyrus and the left inferior frontal gyrus (t-value = 2.73; P = 0.009) (Table 2 and Fig. 2).
Graph theory assessment showed that MS patients have a reduced global efficiency in the right superior frontal gyrus (t-value = − 2.85; P = 0.006), a reduced average path length in the left middle frontal gyrus (t-value = − 3.4; P = 0.001), in the right and left postcentral gyrus (t-value = − 2.79; P = 0.008), and in the left inferior frontal gyrus (t-value = − 2.76; P = 0.008), and a reduced degree in the right superior frontal gyrus (t-value = − 2.85, P = 0.006) compared with HCs for the Ch4 lateral capsula division of the cholinergic network. MS patients also showed an increased global efficiency in the left middle frontal gyrus (t-value = 2.73; P = 0.009) (Table 3 and Fig. 3).
Functional connectivity in dopaminergic network
Executive division
MS patients showed no differences for both functional connectivity and graph theory evaluation for the executive division of the dopaminergic network.
Limbic division
MS patients showed no differences for both functional connectivity and graph theory evaluation for the executive division of the dopaminergic network.
Motor division
MS patients showed no differences in resting-state functional connectivity when compared with HCs for the motor division of the dopaminergic network. Conversely, for the same network, graph theory analysis revealed that MS patients showed a decreased degree in the right middle frontal gyrus compared with HCs (t-value = − 2.86; P = 0.006) (Table 3 and Fig. 3).
Clinical correlation analysis
Disease duration correlated with changes in the FC within the serotonergic–noradrenergic network (between cerebellum and right superior frontal gyrus, t-value = − 3.54, amygdala right and dorsal raphe, t-value = − 3.49, cerebellum and left inferior frontal gyrus, t-value = 3.54, right putamen and left globus pallidus, t-value = -2.92, cerebellum and right globus pallidus, t-value = − 3.29, ventral tegmental area and cerebellum, t-value = 3.54; all P < 0.01; Supplementary Table B-1), and the Ch4 lateral capsula division of the cholinergic network (between left inferior frontal gyrus and right amygdala, t-value = 3.67, right parahippocampal gyrus and left occipital lobe, t-value = − 2.87; all P < 0.01; Supplementary Table B-1). Similarly, disease duration also correlated with changes in the graph theory measures in the serotonergic–noradrenergic network, and the Ch4 lateral capsula and Ch4 lateral perisylvian divisions of the cholinergic network (all P < 0.01; Supplementary Table B-2). EDSS was associated with changes in FC within the serotonergic–noradrenergic network (between cerebellum and left middle frontal gyrus, t-value = − 3.75, right hippocampus and left putamen, t-value = − 3.04, right putamen and right hippocampus, t-value = - 3.52; all P < 0.01; Supplementary Table B-1), the Ch4 lateral capsula division of the cholinergic network (between left superior temporal gyrus and left cuneus, t-value = 3.18, and Ch4 nucleus and left angular gyrus, t-value = 3.25, right amygdala and right lingual gyrus, t-value = − 3.80; all P < 0.01; Supplementary Table B-1). EDSS was also associated with changes in the graph theory measures for the serotonergic–noradrenergic network, for the Ch4 lateral capsula division and Ch4 lateral perisylvian division of the cholinergic network, and for the motor division of the dopaminergic network (all P < 0.01; Supplementary Table B-2).
Disease severity, assessed with MSSS, correlated with changes in the FC in the serotonergic–noradrenergic network (between dorsal raphe and right amygdala, t-value = 3.60, left thalamus and left putamen, t-value = − 3.12, within the cerebellum, t-value = − 2.93; all P < 0.01; Supplementary Table B-1). Disease severity also correlated with changes in the graph theory measures for the serotonergic–noradrenergic network, for the Ch4 lateral perisylvian division of the cholinergic network, and for the motor division of the dopaminergic network (all P < 0.01) (Supplementary Table B-2).
Cognitive function, assessed with the SDMT, correlated with changes in the resting-state FC in the serotonergic–noradrenergic network (between right caudate nucleus and right inferior frontal gyrus, t-value = 4.25, cerebellum and left putamen, t-value = − 6.44, cerebellum and left hippocampus, t-value = − 4.04, right hippocampus and right caudate nucleus, t-value = 5.37, left hypothalamus and right putamen, t-value = 4.12, right inferior frontal gyrus and left superior frontal gyrus, t-value = 4.20, left putamen and left caudate nucleus, t-value = 4.04; right thalamus and cerebellum, t-value = − 4.19; locus coeruleus and right hippocampus, t-value = − 4.16; ventral tegmental area and left superior frontal gyrus, t-value = 6.71; all P < 0.01, Supplementary Table B-1), the Ch4 lateral capsular division of the cholinergic network (between right angular gyrus and left fusiform gyrus, t-value = − 4.52, right fusiform gyrus and right precentral gyrus, t-value = − 4.24, right posterior cingulate cortex and left occipital lobe, t-value = 6.60; all P < 0.01; Supplementary Table B-1), and the limbic division of the dopaminergic network (between left anterior cingulate cortex and left medial orbital gyrus, t-value = 5.77, left medial orbital gyrus and left anterior cingulate cortex, t-value = 3.90; all P < 0.01; Supplementary Table B-1). SDMT also correlated with changes in the graph theory measures for the serotonergic–noradrenergic network and for the Ch4 lateral perisylvian division of the cholinergic network (all P < 0.01; Supplementary Table B-2).
Depressive symptoms, assessed through the HRSD, were associated with changes in the resting-state FC for the serotonergic–noradrenergic network (between cerebellum and left nucleus accumbens, t-value = − 4.75, cerebellum and right hypothalamus, t-value = − 4.17, brainstem and left hypothalamus, t-value = − 5.53, cerebellum and right amygdala, t-value = 10.05, cerebellum and left hippocampus, t-value = − 4.14, cerebellum and right thalamus, t-value = 5.63; cerebellum and locus coeruleus, t-value = − 4.17; cerebellum and ventral tegmental area, t-value = 4.47; all P < 0.01; Supplementary Table B-1), and the Ch4 lateral capsula division of the cholinergic network (between Ch4 nucleus and right angular gyrus, t-value = 4.05, right amygdala and Postcentral gyrus left, t-value = 4.56, left amygdala and left supramarginal gyrus, t-value = 4.07; all P < 0.01; Supplementary Table B-1). HRSD correlated with changes in the graph theory measures for the serotonergic–noradrenergic network (all P < 0.01; Supplementary Table B-2).
Discussion
Using a resting-state fMRI, we demonstrated FC changes and network topography derangements in serotonergic, noradrenergic, cholinergic, and dopaminergic systems. We highlighted that these changes has direct clinical implications resulting in both physical and cognitive disability in MS patients.
In MS, changes in network connectivity and topography within serotonergic, noradrenergic, cholinergic, and dopaminergic networks were not completely unexpected. First, the random occurrence of demyelinating lesions throughout the brain may cause white matter tract disconnection and brain reorganization, resulting in impaired FC. Considering the many long-range and diffuse projections to both cortical and subcortical regions from the relatively small cell populations that produce modulatory neurotransmitters (e.g., dorsal Raphe, Ch4 nucleus, and locus coeruleus), it is reasonable to expect that these neurochemical networks would be highly susceptible to white matter pathologies. Furthermore, given the wide range of physiological processes that these modulatory transmitter populations are implicated in it, quantification of functional changes in their associated networks could well provide a guide to the development and progression of clinical symptoms in MS. However, pathological changes lead to structural and functional reorganization of the brain to compensate for the loss of functions. Brain reorganization in MS is a dynamic process, likely characterized by increased FC as a compensatory mechanism during early stages and reduced ability to compensate for progressive brain damage at later stages, subsequently resulting in decreased FC [14]. Therefore, we cannot completely rule out the role of spatial reorganization of network anatomy in MS patients following brain damages. However, to partly overcome this issue, we used subject-specific regions of interest map, which allow us to take into account the subject anatomical variability and brain changes that might occur during the disease course, for example brain atrophy. Although ROI selection was performed to avoid brain regions overlapping, between the neuromodulatory networks due to the fact that in one single cortical area, we might find neuronal subpopulation belonging to different neuromodulatory networks, we cannot completely disentangle the role of each sub-population to the connectivity within the explored network. Notwithstanding this limitation, our exploratory study might drive new studies aimed at evaluating neuromodulatory network changes in MS taking advantage of newly developed PET tracers for neurotransmitters which are better equipped to dissociate neuromodulatory components in, for example, the dorsal raphe or cholinergic nuclei. Similarly, the reported alterations in FC might not necessarily imply a selective derangement in neuromodulatory networks, since they could reflect changes in non-neuromodulatory networks partly overlapping in terms of brain regions. Nonetheless these limitations, our results are consistent with the previous results. Hesse and colleagues reported reduced brain serotonin transporter levels in the cingulate cortex, the thalamus, and the insula, and a higher transporter availability in the orbitofrontal cortex in MS patients compared to HCs using PET molecular imaging with the radiotracer [11C]DASB [34]. In our study, we found reduced functional connectivity within the cerebellum, and between the cerebellum and both the amygdala and the thalamus, whereas there was an increased connectivity between the dorsal raphe and the frontal regions. In addition, we also outlined a decreased betweenness centrality for the brainstem in MS compared to HCs, suggesting that ascending projections from this structure are damaged and have a reduced activity in controlling the cortical regions of the network.
Differently for the serotonergic network, very little is known about noradrenalin deficits in MS. Noradrenalin concentration in the brain of MS patients was found to be either increased [35] or decreased [36]. The primary source of noradrenalin in the brain is the locus coeruleus [37]. Both inflammation and neurodegeneration have been described within the locus coeruleus of MS patients and the extent of neuronal damage shown to correlate with cognitive disability [38, 39]. Therefore, the reduced centrality of the brainstem in the noradrenalin network, reported in here, might be related to an overall reduced activity of the locus coeruleus with reduced function of the long-range projections heading to frontal regions and adaptive short-path intra-cerebellum adaptive connections, which increase the dysfunctional centrality of the cerebellum.
In our study, we also reported a functional disconnection between cortical and subcortical regions modulated by cholinergic neurons, with direct implications for cognition. Kimura and colleagues already reported tissue damage in the white matter bundles of the cholinergic network using a semi-quantitative measure and a correlation with cognitive impairment [40]. Finally, we also reported altered FC in the dopaminergic network in MS patients compared to HCs, showing a decreased connectivity between the ventral tegmental area and the orbital cortex and a decreased degree in the middle frontal gyrus. Up to now, no studies have specifically assessed dopaminergic network in MS. However, indirect report outlined the occurrence of white matter lesions over the bundles connecting the substantia nigra, the ventral tegmental areas, the striatum, the limbic areas, and the prefrontal cortex with changes in the local metabolic activity in the prefrontal cortex and in the striatum [41,42,43]. These damages may underpin FC changes and network reorganization, and are associated with cognitive impairment and fatigue [44,45,46].
Regarding clinical outcome measures, disease duration, disease severity, and physical disability were also associated at different extent with FC and network topography changes within serotonergic, cholinergic, dopaminergic, and noradrenergic network. We also reported a link between networks derangement and SDMT scores, suggesting a role for cognitive impairment. However, the SDMT is a single test assessing attention and information processing speed, thus exploring only a single cognitive domain. A more extensive cognitive battery, assessing different cognitive domains, such as the Brief Repeatable Battery of Neuropsychological Tests [47], could better elucidate the relationship between cognitive impairment and networks derangements. The association between network impairment and disease burden was evident especially for those networks with a wider distribution throughout the brain. This might be possibly due to a lower probability of demyelinating lesions over the connecting fiber bundles and subsequent disruption of neuron communication. Moreover, it is worth mentioning that both disease modifying and concomitant treatments, especially those aiming at providing relief from depressive symptoms, might impact on FC in neuromodulatory networks. The small sample in the present study prevents us from being able to draw conclusions on this topic. Further studies with larger samples, and targeted study design, are warranted to better evaluate the impact of disease modifying and concomitant therapies on FC changes in neuromodulatory networks. Therefore, the overall neurotransmitter network impairment might be used as a marker of brain damage and could be useful to stratify patients according to the severity of the disease.
In addition, we reported an association between depressive symptoms and both serotonergic and noradrenergic network topography reorganization. Our findings are in line with the reduction in the monoamine network activity, namely noradrenalin and serotonin, in major depressive disorder [48] and in MS [34]. In addition, we found that depression in MS was associated with an increased connectivity throughout the cholinergic network. This is in line with the previous findings describing an association between depressive symptoms and over-activity of the cholinergic system throughout the brain [49, 50]. Taken together, these findings suggest that depressive symptoms in MS are highly dependent upon neurotransmission integrity and an integrated therapeutic approach using drugs targeting different neurotransmitters might be useful. Moreover, depressive symptoms in MS include a broad spectrum of symptoms such as decreased concentration, fatigue, and sexual disinterest. Understanding better the neuronal basis of each of these symptoms might help clinicians towards more targeted symptomatic treatment choices to obtain a higher treatment efficacy.
Conclusions
Through utilization of resting-state fMRI, and adopting an innovative method for reconstructing neuromodulatory networks maps in a subject-specific fashion, we revealed a diffuse FC impairment and topography derangements within the serotonergic, noradrenergic, cholinergic, and dopaminergic networks in MS with direct impact on clinical disability. This is the first study investigating neurotransmitter pathology in MS. Considering that resting-state fMRI might not disentangle FC changes specifically due to neuro-modulatory damages from those related to non-neuromodulatory connections, further studies assessing the impairment for each individual network using selective PET tracers are highly warranted to better understand the role of each network in disease pathophysiology and drive towards clinical applications of these findings. Neurotransmitter impairment might be a predictor of disability burden over time and might be used ideally to assess treatment efficacy. Moreover, confirmatory studies might also lead to the development of new symptomatic treatment based on neurotransmitter modulation.
Data availability
Data will be made available via a request to the authors after a formal data sharing agreement.
References
Ciccarelli O, Barkhof F, Bodini B, De Stefano N, Golay X, Nicolay K, Pelletier D, Pouwels PJ, Smith SA, Wheeler-Kingshott CA, Stankoff B, Yousry T, Miller DH (2014) Pathogenesis of multiple sclerosis: insights from molecular and metabolic imaging. Lancet Neurol 13(8):807–822. https://doi.org/10.1016/S1474-4422(14)70101-2
Palop JJ, Chin J, Mucke L (2006) A network dysfunction perspective on neurodegenerative diseases. Nature 443(7113):768–773. https://doi.org/10.1038/nature05289
Dutta R, Chang A, Doud MK, Kidd GJ, Ribaudo MV, Young EA, Fox RJ, Staugaitis SM, Trapp BD (2011) Demyelination causes synaptic alterations in hippocampi from multiple sclerosis patients. Ann Neurol 69(3):445–454. https://doi.org/10.1002/ana.22337
Caminiti SP, Tettamanti M, Sala A, Presotto L, Iannaccone S, Cappa SF, Magnani G, Perani D, Alzheimer's Disease Neuroimaging I (2017) Metabolic connectomics targeting brain pathology in dementia with Lewy bodies. J Cereb Blood Flow Metab 37(4):1311–1325. https://doi.org/10.1177/0271678X16654497
Klein JC, Eggers C, Kalbe E, Weisenbach S, Hohmann C, Vollmar S, Baudrexel S, Diederich NJ, Heiss WD, Hilker R (2010) Neurotransmitter changes in dementia with Lewy bodies and Parkinson disease dementia in vivo. Neurology 74(11):885–892. https://doi.org/10.1212/WNL.0b013e3181d55f61
Pagano G, Niccolini F, Politis M (2017) The serotonergic system in Parkinson's patients with dyskinesia: evidence from imaging studies. J Neural Transm (Vienna). https://doi.org/10.1007/s00702-017-1823-7
Niccolini F, Pagano G, Fusar-Poli P, Wood A, Mrzljak L, Sampaio C, Politis M (2018) Striatal molecular alterations in HD gene carriers: a systematic review and meta-analysis of PET studies. J Neurol Neurosurg Psychiatry 89(2):185–196. https://doi.org/10.1136/jnnp-2017-316633
Roy R, Niccolini F, Pagano G, Politis M (2016) Cholinergic imaging in dementia spectrum disorders. Eur J Nucl Med Mol Imaging 43(7):1376–1386. https://doi.org/10.1007/s00259-016-3349-x
Bhattacharyya PK, Phillips MD, Stone LA, Bermel RA, Lowe MJ (2013) Sensorimotor cortex gamma-aminobutyric acid concentration correlates with impaired performance in patients with MS. AJNR Am J Neuroradiol 34(9):1733–1739. https://doi.org/10.3174/ajnr.A3483
Cawley N, Solanky BS, Muhlert N, Tur C, Edden RA, Wheeler-Kingshott CA, Miller DH, Thompson AJ, Ciccarelli O (2015) Reduced gamma-aminobutyric acid concentration is associated with physical disability in progressive multiple sclerosis. Brain 138(Pt 9):2584–2595. https://doi.org/10.1093/brain/awv209
Muhlert N, Atzori M, De Vita E, Thomas DL, Samson RS, Wheeler-Kingshott CA, Geurts JJ, Miller DH, Thompson AJ, Ciccarelli O (2014) Memory in multiple sclerosis is linked to glutamate concentration in grey matter regions. J Neurol Neurosurg Psychiatry 85(8):833–839. https://doi.org/10.1136/jnnp-2013-306662
Faivre A, Rico A, Zaaraoui W, Crespy L, Reuter F, Wybrecht D, Soulier E, Malikova I, Confort-Gouny S, Cozzone PJ, Pelletier J, Ranjeva JP, Audoin B (2012) Assessing brain connectivity at rest is clinically relevant in early multiple sclerosis. Mult Scler 18(9):1251–1258. https://doi.org/10.1177/1352458511435930
Roosendaal SD, Schoonheim MM, Hulst HE, Sanz-Arigita EJ, Smith SM, Geurts JJ, Barkhof F (2010) Resting state networks change in clinically isolated syndrome. Brain 133(Pt 6):1612–1621. https://doi.org/10.1093/brain/awq058
Rocca MA, Valsasina P, Martinelli V, Misci P, Falini A, Comi G, Filippi M (2012) Large-scale neuronal network dysfunction in relapsing-remitting multiple sclerosis. Neurology 79(14):1449–1457. https://doi.org/10.1212/WNL.0b013e31826d5f10
Rocca MA, Valsasina P, Leavitt VM, Rodegher M, Radaelli M, Riccitelli GC, Martinelli V, Martinelli-Boneschi F, Falini A, Comi G, Filippi M (2018) Functional network connectivity abnormalities in multiple sclerosis: correlations with disability and cognitive impairment. Mult Scler 24(4):459–471. https://doi.org/10.1177/1352458517699875
d'Ambrosio A, Hidalgo de la Cruz M, Valsasina P, Pagani E, Colombo B, Rodegher M, Comi G, Filippi M, Rocca MA (2017) Structural connectivity-defined thalamic subregions have different functional connectivity abnormalities in multiple sclerosis patients: implications for clinical correlations. Hum Brain Mapp 38(12):6005–6018. https://doi.org/10.1002/hbm.23805
Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52(3):1059–1069. https://doi.org/10.1016/j.neuroimage.2009.10.003
O'Keeffe GW, Sullivan AM (2018) Evidence for dopaminergic axonal degeneration as an early pathological process in Parkinson's disease. Parkinsonism Relat Disord 56:9–15. https://doi.org/10.1016/j.parkreldis.2018.06.025
Hampel H, Mesulam MM, Cuello AC, Farlow MR, Giacobini E, Grossberg GT, Khachaturian AS, Vergallo A, Cavedo E, Snyder PJ, Khachaturian ZS (2018) The cholinergic system in the pathophysiology and treatment of Alzheimer's disease. Brain 141(7):1917–1933. https://doi.org/10.1093/brain/awy132
Polman CH, Reingold SC, Banwell B, Clanet M, Cohen JA, Filippi M, Fujihara K, Havrdova E, Hutchinson M, Kappos L, Lublin FD, Montalban X, O'Connor P, Sandberg-Wollheim M, Thompson AJ, Waubant E, Weinshenker B, Wolinsky JS (2011) Diagnostic criteria for multiple sclerosis: 2010 revisions to the McDonald criteria. Ann Neurol 69(2):292–302. https://doi.org/10.1002/ana.22366
Roxburgh RH, Seaman SR, Masterman T, Hensiek AE, Sawcer SJ, Vukusic S, Achiti I, Confavreux C, Coustans M, le Page E, Edan G, McDonnell GV, Hawkins S, Trojano M, Liguori M, Cocco E, Marrosu MG, Tesser F, Leone MA, Weber A, Zipp F, Miterski B, Epplen JT, Oturai A, Sorensen PS, Celius EG, Lara NT, Montalban X, Villoslada P, Silva AM, Marta M, Leite I, Dubois B, Rubio J, Butzkueven H, Kilpatrick T, Mycko MP, Selmaj KW, Rio ME, Sa M, Salemi G, Savettieri G, Hillert J, Compston DA (2005) Multiple sclerosis severity score: using disability and disease duration to rate disease severity. Neurology 64(7):1144–1151. https://doi.org/10.1212/01.WNL.0000156155.19270.F8
Van Schependom J, D'Hooghe MB, Cleynhens K, D'Hooge M, Haelewyck MC, De Keyser J, Nagels G (2014) The Symbol Digit Modalities Test as sentinel test for cognitive impairment in multiple sclerosis. Eur J Neurol 21(9):1219–1225. https://doi.org/10.1111/ene.12463(e1271–1212)
Chard DT, Jackson JS, Miller DH, Wheeler-Kingshott CA (2010) Reducing the impact of white matter lesions on automated measures of brain gray and white matter volumes. J Magn Reson Imaging 32(1):223–228. https://doi.org/10.1002/jmri.22214
Smith SM, Zhang Y, Jenkinson M, Chen J, Matthews PM, Federico A, De Stefano N (2002) Accurate, robust, and automated longitudinal and cross-sectional brain change analysis. Neuroimage 17(1):479–489
Heckemann RA, Keihaninejad S, Aljabar P, Rueckert D, Hajnal JV, Hammers A, Alzheimer's Disease Neuroimaging I (2010) Improving intersubject image registration using tissue-class information benefits robustness and accuracy of multi-atlas based anatomical segmentation. Neuroimage 51(1):221–227. https://doi.org/10.1016/j.neuroimage.2010.01.072
Whitfield-Gabrieli S, Nieto-Castanon A (2012) Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect 2(3):125–141. https://doi.org/10.1089/brain.2012.0073
Svob Strac D, Pivac N, Muck-Seler D (2016) The serotonergic system and cognitive function. Transl Neurosci 7(1):35–49. https://doi.org/10.1515/tnsci-2016-0007
Jacobs BL, Azmitia EC (1992) Structure and function of the brain serotonin system. Physiol Rev 72(1):165–229. https://doi.org/10.1152/physrev.1992.72.1.165
Selden NR, Gitelman DR, Salamon-Murayama N, Parrish TB, Mesulam MM (1998) Trajectories of cholinergic pathways within the cerebral hemispheres of the human brain. Brain 121(Pt 12):2249–2257
Samuels ER, Szabadi E (2008) Functional neuroanatomy of the noradrenergic locus coeruleus: its roles in the regulation of arousal and autonomic function part I: principles of functional organisation. Curr Neuropharmacol 6(3):235–253. https://doi.org/10.2174/157015908785777229
Tziortzi AC, Haber SN, Searle GE, Tsoumpas C, Long CJ, Shotbolt P, Douaud G, Jbabdi S, Behrens TE, Rabiner EA, Jenkinson M, Gunn RN (2014) Connectivity-based functional analysis of dopamine release in the striatum using diffusion-weighted MRI and positron emission tomography. Cereb Cortex 24(5):1165–1177. https://doi.org/10.1093/cercor/bhs397
Alaerts K, Geerlings F, Herremans L, Swinnen SP, Verhoeven J, Sunaert S, Wenderoth N (2015) Functional organization of the action observation network in autism: a graph theory approach. PLoS ONE 10(8):e0137020. https://doi.org/10.1371/journal.pone.0137020
Haynes W (2013) Benjamini-Hochberg method. In: Dubitzky W, Wolkenhauer O, Cho K-H, Yokota H (eds) Encyclopedia of systems biology. Springer, New York, New York, p 78. https://doi.org/10.1007/978-1-4419-9863-7_1215
Hesse S, Moeller F, Petroff D, Lobsien D, Luthardt J, Regenthal R, Becker GA, Patt M, Thomae E, Seese A, Meyer PM, Bergh FT, Sabri O (2014) Altered serotonin transporter availability in patients with multiple sclerosis. Eur J Nucl Med Mol Imaging 41(5):827–835. https://doi.org/10.1007/s00259-013-2636-z
Barkhatova VP, Zavalishin IA, Askarova L, Shavratskii V, Demina EG (1998) Changes in neurotransmitters in multiple sclerosis. Neurosci Behav Physiol 28(4):341–344
Markianos M, Koutsis G, Evangelopoulos ME, Mandellos D, Karahalios G, Sfagos C (2009) Relationship of CSF neurotransmitter metabolite levels to disease severity and disability in multiple sclerosis. J Neurochem 108(1):158–164. https://doi.org/10.1111/j.1471-4159.2008.05750.x
Benarroch EE (2009) The locus ceruleus norepinephrine system: functional organization and potential clinical significance. Neurology 73(20):1699–1704. https://doi.org/10.1212/WNL.0b013e3181c2937c
Polak PE, Kalinin S, Feinstein DL (2011) Locus coeruleus damage and noradrenaline reductions in multiple sclerosis and experimental autoimmune encephalomyelitis. Brain 134(Pt 3):665–677. https://doi.org/10.1093/brain/awq362
Gadea M, Martinez-Bisbal MC, Marti-Bonmati L, Espert R, Casanova B, Coret F, Celda B (2004) Spectroscopic axonal damage of the right locus coeruleus relates to selective attention impairment in early stage relapsing-remitting multiple sclerosis. Brain 127(Pt 1):89–98. https://doi.org/10.1093/brain/awh002
Kimura Y, Sato N, Ota M, Maikusa N, Maekawa T, Sone D, Enokizono M, Sugiyama A, Imabayashi E, Matsuda H, Okamoto T, Yamamura T, Sugimoto H (2017) A structural MRI study of cholinergic pathways and cognition in multiple sclerosis. eNeurological Sci 8:11–16. https://doi.org/10.1016/j.ensci.2017.06.008
Roelcke U, Kappos L, Lechner-Scott J, Brunnschweiler H, Huber S, Ammann W, Plohmann A, Dellas S, Maguire RP, Missimer J, Radu EW, Steck A, Leenders KL (1997) Reduced glucose metabolism in the frontal cortex and basal ganglia of multiple sclerosis patients with fatigue: a 18F-fluorodeoxyglucose positron emission tomography study. Neurology 48(6):1566–1571
Engstrom M, Flensner G, Landtblom AM, Ek AC, Karlsson T (2013) Thalamo-striato-cortical determinants to fatigue in multiple sclerosis. Brain Behav 3(6):715–728. https://doi.org/10.1002/brb3.181
Finke C, Schlichting J, Papazoglou S, Scheel M, Freing A, Soemmer C, Pech LM, Pajkert A, Pfuller C, Wuerfel JT, Ploner CJ, Paul F, Brandt AU (2015) Altered basal ganglia functional connectivity in multiple sclerosis patients with fatigue. Mult Scler 21(7):925–934. https://doi.org/10.1177/1352458514555784
Pardini M, Bonzano L, Mancardi GL, Roccatagliata L (2010) Frontal networks play a role in fatigue perception in multiple sclerosis. Behav Neurosci 124(3):329–336. https://doi.org/10.1037/a0019585
Genova HM, Rajagopalan V, Deluca J, Das A, Binder A, Arjunan A, Chiaravalloti N, Wylie G (2013) Examination of cognitive fatigue in multiple sclerosis using functional magnetic resonance imaging and diffusion tensor imaging. PLoS ONE 8(11):e78811. https://doi.org/10.1371/journal.pone.0078811
Vytlacil J, Kayser A, Miyakawa A, D'Esposito M (2014) An approach for identifying brainstem dopaminergic pathways using resting state functional MRI. PLoS ONE 9(1):e87109. https://doi.org/10.1371/journal.pone.0087109
Boringa JB, Lazeron RH, Reuling IE, Ader HJ, Pfennings L, Lindeboom J, de Sonneville LM, Kalkers NF, Polman CH (2001) The brief repeatable battery of neuropsychological tests: normative values allow application in multiple sclerosis clinical practice. Mult Scler 7(4):263–267. https://doi.org/10.1177/135245850100700409
Schildkraut JJ (1965) The catecholamine hypothesis of affective disorders: a review of supporting evidence. Am J Psychiatry 122(5):509–522. https://doi.org/10.1176/ajp.122.5.509
Pytka K, Dziubina A, Mlyniec K, Dziedziczak A, Zmudzka E, Furgala A, Olczyk A, Sapa J, Filipek B (2016) The role of glutamatergic, GABA-ergic, and cholinergic receptors in depression and antidepressant-like effect. Pharmacol Rep 68(2):443–450. https://doi.org/10.1016/j.pharep.2015.10.006
Pytka K, Podkowa K, Rapacz A, Podkowa A, Zmudzka E, Olczyk A, Sapa J, Filipek B (2016) The role of serotonergic, adrenergic and dopaminergic receptors in antidepressant-like effect. Pharmacol Rep 68(2):263–274. https://doi.org/10.1016/j.pharep.2015.08.007
Funding
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Author information
Authors and Affiliations
Contributions
MP and AC: conceptualization. BG, PB, and ES: data acquisition. AC, ZC, SPC, and HW: data analysis. AC, HW, PB, SPC, ES, SCRW, AH, and MP: manuscript drafting.
Corresponding author
Ethics declarations
Conflicts of interest
Antonio Carotenuto, Heather Wilson, Beniamino Giordano, Silvia P. Caminiti, Zachary Chappell, Steven C.R. Williams, and Alexander Hammers declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this abstract. Dr Eli Silber and Peter Brex received funding to attend lectures, advisory board, and for lecturing from Biogen, Merck, Teva, Genzyme, and Roche. Professor Marios Politis research is supported by Parkinson’s UK, Lily and Edmond J. Safra Foundation, Michael J Fox Foundation (MJFF) for Parkinson’s research, and CHDI Foundation.
Ethical approval
The study was approved by the London-South East Research Ethics Committee.
Informed consent
Written informed consent was obtained from all study’s subjects prior to inclusion in the study in accordance to the Declaration of Helsinki.
Electronic supplementary material
Below is the link to the electronic supplementary material.
415_2020_9806_MOESM1_ESM.docx
Supplementary file1 includes an in-depth explanation on the methodology we used to set the network maps and additional tables containing ROI-based connectivity analysis (appendix A-1 and Table B-1 and B-2) and correlations with clinical outcome measures (DOCX 95 kb)
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 http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Carotenuto, A., Wilson, H., Giordano, B. et al. Impaired connectivity within neuromodulatory networks in multiple sclerosis and clinical implications. J Neurol 267, 2042–2053 (2020). https://doi.org/10.1007/s00415-020-09806-3
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00415-020-09806-3