Linking microstructural integrity and motor cortex excitability in multiple sclerosis

Angela Radetz, Kalina Mladenova, Dumitru Ciolac, Gabriel Gonzalez-Escamilla, Vinzenz Fleischer, Erik Ellwardt, Julia Krämer, Stefan Bittner, Sven G. Meuth, Muthuraman Muthuraman, Sergiu Groppa
doi: https://doi.org/10.1101/2020.10.15.20213090

Abstract

Motor skills are frequently impaired in multiple sclerosis (MS) patients following grey (GM) and white matter (WM) damage with cortical excitability abnormalities. We performed advanced diffusion imaging for neurite orientation dispersion and density modeling and diffusion tensor imaging within the motor system of 50 MS patients and 49 age-matched healthy controls. To assess excitability, we determined resting motor thresholds using non-invasive transcranial magnetic stimulation. A hierarchical regression model revealed that lower neurite density index (NDI), suggestive for axonal loss in the GM, predicted higher motor thresholds, i.e. reduced excitability in MS patients. Furthermore, lower NDI was indicative of decreased cognitive-motor performance. Interconnected motor WM tracts of patients were characterized by overlapping clusters of lowered fractional anisotropy and NDI, with NDI exclusively capturing a higher amount of abnormally appearing voxels. Our work outlines the potential of microstructure imaging using advanced biophysical models to forecast neurodegeneration and excitability alterations in neuroinflammation.

Introduction

Impaired motor functions caused by both grey (GM) and white matter (WM) pathology are characteristic of multiple sclerosis (MS), a neurodegenerative inflammatory disorder of the central nervous system (Eshaghi et al., 2018; Filippi et al., 2012; Kurtzke, 1983). Both in animal models of MS and in patients, cortical neurons show abnormal levels of excitability, indicating a pathophysiological relevance in neuroinflammation (Caramia et al., 2004; Ellwardt et al., 2018). Excitability of neural tissue in motor brain regions can be measured in-vivo by determining the resting motor threshold with transcranial magnetic stimulation (TMS). Hereby, increased resting motor thresholds were frequently observed in MS patients (Stampanoni Bassi et al., 2020; Zipser et al., 2018). Although the driving mechanisms of altered cortical excitability are still unclear, cortical neuronal or axonal loss or altered (re-) myelination processes could play an important role (Calabrese et al., 2007; Zipser et al., 2018). Indeed, GM integrity alterations are a pathophysiological signature of MS; atrophy quantified with conventional imaging sequences is commonly used to track disease courses (Eshaghi et al., 2018; Filippi et al., 2012; Radetz et al., 2019). Whereas this measure does not allow drawing conclusions about the depicted microstructure, non-conventional imaging approaches allow an in-vivo quantification of microstructural alterations caused by neurodegenerative processes (Alexander et al., 2019). Monitoring and assessing microstructural alterations might improve the understanding of altered excitability levels. With diffusion-weighted imaging (DWI) sequences, WM microstructural integrity has been robustly quantified based on parameters derived from fitted diffusion tensors (diffusion tensor imaging [DTI]), such as fractional anisotropy (FA) (Deppe, Krämer, et al., 2016; Deppe, Tabelow, et al., 2016; Droby et al., 2015; Roosendaal et al., 2009). Whereas FA is generally sensitive to microstructural neuroinflammation-or neurodegeneration-driven tissue damage, it is relatively non-specific for the underlying pathophysiological processes (Krämer et al., 2019). FA does not distinguish restricted and hindered diffusion and is thereby biased in areas with high neurite orientation dispersion such as the GM or regions of crossing fibers (Tournier et al., 2011; Zhang et al., 2012). Advanced biophysical diffusion models that are continuously refined allow a more realistic in-vivo depiction of GM and WM microstructure including estimations of neurite density, axon diameter, neurite dispersion and cell size (Alexander et al., 2019; Lakhani et al., 2020). The neurite orientation dispersion and density imaging (NODDI) model, which is based on advanced acquisition techniques, mitigates the methodological constraints as described for FA (Zhang et al., 2012). In each voxel a neurite density (NDI) and orientation dispersion index (ODI) is computed from restricted diffusion, whereas an isotropic volume fraction (IVF) reflects cerebrospinal fluid and edema (Zhang et al., 2012). A mouse diffusion imaging study showed good agreement of the model with parallel histological quantification of microstructure. High ODI values were reported in areas of complex architecture or crossing fibers, and high NDI values in regions with low neuron content and dense fiber bundles (Wang et al., 2019). A cross-modality study of the human cortical GM suggested a strong link between NDI and the density of myelinated, and to a lesser extent of unmyelinated axons. ODI was proposed to reflect cyto- and myeloarchitecture (Fukutomi et al., 2018). Importantly, both ODI and NDI contribute to FA, thereby allowing a more specific differentiation of microstructural characteristics (Zhang et al., 2012). Here, we aimed at carrying forward the potential of the NODDI model to evaluate the microstructural integrity in the GM of primary motor cortex and address the question, whether lower excitability is related to pathological changes characterized by NDI, ODI or IVF. We further hypothesized that NODDI measures used for quantification of motor microstructural integrity are predictive for motor and cognitive function. Such assessments are of pivotal importance as the applied clinically feasible diffusion sequence could particularly forecast neurodegenerative processes that lead to progression independent of relapse activity.

As neuroinflammation with de- and remyelination, as well as neurodegeneration occurs in GM and WM, we emphasize the importance of gaining a robust depiction of GM and WM pathology specifically of the motor system for the prediction of the functional trajectories. At the whole-brain level, reduced NDI in the normal-appearing WM of MS patients compared to that of healthy controls (HC) was previously reported, indicating myelin and axonal loss (Granberg et al., 2017; Spano et al., 2018). ODI generally yielded more variable results during comparisons between MS and HC (De Santis et al., 2019; Lakhani et al., 2020; Spano et al., 2018). Alterations in ODI result from morphological changes including a combination of axonal loss, dendritic pruning and tissue reorganization. Therefore, conclusions about the pathological processes strongly depend on the original tissue architecture (Spano et al., 2018). We here compared microstructural integrity in motor WM tracts of MS patients and HC using tract-based spatial statistics (TBSS) (Smith et al., 2006). This method allows for voxel-wise group comparisons of FA and NODDI parameters in WM tracts by minimizing partial volume effects. Importantly, we investigated if NODDI values in contrast to the classical FA are more sensitive in detecting differences in microstructural integrity in MS patients versus HC. As FA increases might reflect either NDI increases, ODI decreases or a combination of both (Zhang et al., 2012), we tracked contributions of these signals to FA and evaluated if the NODDI model captures additional pathological alterations due to MS pathology. This approach allows more specific inferences about the type of microstructural pathology underlying the FA changes.

We first assessed the predictive value of motor cortex microstructure for motor excitability and cognitive-motor performance in our participants. Next, we evaluated pathological microstructure alterations in motor WM tracts of MS patients and HC using NODDI in comparison to the classical FA. We hereby examined voxel-wise intersecting group differences to disentangle potential contributions of NODDI metrics to FA.

Methods

Participants

Fifty patients with relapsing-remitting MS (RRMS) (31 female, mean age = 35.3 years, SD = 11.1), diagnosed according to the revised McDonald criteria (Thompson et al., 2018), and 49 HC (25 female, mean age = 31.5 years, SD = 9.1) participated in our study. Clinical data of patients, Nine Hole Peg Test (9HPT), Trail Making Test part A (TMT-A) and part B (TMT-B) raw scores (Reitan, 1971) and an overview of ongoing therapies at the time of MRI acquisition are depicted in Table 1. Participants were excluded from the study in case of pregnancy or relapses or systemic therapy with steroids within the month before MRI. They were also excluded in case of contraindication to TMS, i.e. previous or current susceptibility to epileptic seizures or metal implants in proximity to the cranium. Expanded Disability Status Scale (EDSS) scores (Kurtzke, 1983) of MS patients were determined by trained and certified neurologists. The study was approved by the local Ethics Committee of the State Medical Association of Rhineland-Palatinate. All participants gave written informed consent before participation and were informed about the study content in oral and written form.

View this table:
Table 1.

Overview of patients’ clinical characteristics

MRI

Scanning was performed in a 3T Siemens TrioTim MRI scanner with a 32-channel head coil (Siemens Healthcare, Erlangen, Germany). A 3D T1-weighted magnetization prepared rapid gradient echo (MP-RAGE) sequence (echo time [TE] = 2.52 ms, repetition time [TR] = 1900 ms, inversion time [TI] = 900 ms, flip angle = 9°, matrix size = 256 × 256, field of view [FOV] = 256 × 256 mm2, slice thickness = 1 mm, voxel size = 1 × 1 x 1 mm3) and a sagittal 3D turbo spin-echo fluid attenuated inversion recovery (FLAIR) sequence (TE = 388 ms, TR = 5000 ms, TI = 1800 ms, flip angle = 120°, matrix size = 256 × 258, FOV = 256 × 256 mm2, slice thickness = 1 mm, voxel size = 0.5 × 0.5 × 1 mm3) were acquired. Further, multi-shell DWI data were obtained (TE = 111 ms, TR = 10800 ms, flip angle = 90°, matrix size = 128 × 128, FOV = 2304 × 2304 mm2, slice thickness = 2 mm, voxel size = 2 × 2 x 2 mm3), with three non-zero b-values each measured in thirty unique directions (b = 900 s/mm2, b = 1800 s/mm2, b = 2700 s/mm2) in the anterior-posterior direction. Six non diffusion-weighted volumes were acquired before each change in b-value in the anterior-posterior direction. Finally, one further non diffusion-weighted volume was acquired in anterior-posterior, and one in posterior-anterior direction.

Transcranial magnetic stimulation and recordings of motor evoked potentials

Electromyography (EMG) electrodes were placed over the right first dorsal interosseous (FDI) and abductor pollicis brevis (APB). EMG signals including the motor evoked potentials (MEP) were 1000-fold amplified with a D440 amplifier from Digitimer (Fort Lauderdale, USA). Using a CED 1401 laboratory interface, signals were digitized with a sampling rate of 5 kHz (Cambridge Electronic Design, Cambridge, UK). TMS was applied using a Rapid2 Stimulator with a figure-eight coil in biphasic pulse configuration (Magstim®, Whitland, UK). The TMS coil was placed over the hand area of left primary motor cortex (M1) approximately in parallel to the central sulcus, i.e. 45-55° relative to the mid-sagittal line (Groppa et al., 2012). After each suprathreshold pulse consistently evoked an MEP in the right FDI and APB, the individual resting motor threshold was determined. To that end, pulse intensity was lowered until MEPs with peak-to-peak amplitudes of 50 μV were evoked for two out of four pulses at rest (Wilson et al., 1993).

Neuropsychological assessment

In the MS group, 31 patients underwent neuropsychological cognitive assessment of 9HPT, and 44 patients of TMT-A and TMT-B (time interval between MRI and neuropsychological testing was mean = 3.8 months, SD = 3.1 months). The 9HPT is a frequently used measure of manual dexterity that is often impaired in MS (Feys et al., 2017). Participants are asked to place and remove pegs into holes on a board as quickly as possible. We used 9HPT scores of the dominant (right) hand as we accordingly investigated the left-hemispheric M1 cortex. In the TMT, patients are instructed to connect a set of numbers and letters according to specific rules as fast as possible. Whereas we included the test due to the evaluation of motor speed, it also requires visual attention and task switching abilities. Average raw scores and standard deviations are reported in Table 1. We used z-scores of these variables by comparing the raw scores with test-specific normative data stratified for age and education.

Data preprocessing

Figure 1 provides an overview of the study pipeline.

Figure 1.
Figure 1.

Study pipeline. (1) T1-weighted images were used for reconstruction of the GM that was individually coregistered to FA, NDI, ODI and IVF maps. Motor threshold as measure of excitability was regressed on average FA, NDI, ODI and IVF within the GM of left M1 (1a). Scores of cognitive-motor performance were similarly regressed on average FA, NDI, ODI and IVF within the GM of left M1 (1b). (2) Individual FA maps of all subjects were used to compute the TBSS mean skeleton that was applied to FA, NDI, ODI and IVF for comparisons of MS and HC group. Percentage of intersecting significantly different voxels with SMATT regions were determined for each diffusion parameter. Abbreviations: FA, fractional anisotropy, NDI, neurite density index, ODI, orientation dispersion index, IVF, isotropic volume fraction, MS, multiple sclerosis, HC, healthy control, TBSS, tract-based spatial statistics, SMATT, sensorimotor area tract template, GM, grey matter, M1, primary motor cortex.

Brain segmentation

Brain segmentation based on T1-weighted images was conducted using the cross-sectional processing stream of FreeSurfer version 5.3.0 (http://surfer.nmr.mgh.harvard.edu/). Here, removal of non-brain tissue, automated Talairach transformation, cortical and subcortical segmentation, intensity normalization, tessellation of the grey and white matter boundary and automated topology correction is included (Fischl, 2012). We used the FreeSurfer-based GM segmentation for further analyses.

Preprocessing of DWI

DWI data were preprocessed using the diffusion toolbox of FSL version 5.0.9 (https://fsl.fmrib.ox.ac.uk/fsl/) by correcting susceptibility induced distortions using topup and applying eddy current and motion artifact correction using eddy. Voxel-wise diffusion tensor fitting and computation of FA was performed using FSL’s dtifit only considering the inner shell (b = 900) and the non diffusion-weighted images.

Application of the NODDI model

Zhang and colleagoues developed the NODDI model with the requirement of being sufficiently simple, while being complex enough for depiction of major characteristics of neurite morphology (Zhang et al., 2012). A further aim was a clinically feasible acquisition time below 30 minutes. We fitted the model to the DWI using the Accelerated Microstructure Imaging via Convex Optimization (AMICO) algorithm resulting in NDI, ODI and IVF maps (Daducci et al., 2015). Whereas IVF is characterized by free isotropic Gaussian diffusion, NDI and ODI are computed from restricted diffusion adapted from the orientation-dispersed cylinder model (Zhang et al., 2011, 2012). We coregistered the first non diffusion-weighted image to the brainmask obtained with FreeSurfer using FSL’s FMRIB’s Linear Image Registration Tool (FLIRT). The transformation matrix was then applied to FA, NDI, ODI and IVF maps.

Application of a motor atlas

The Human Motor Area Template (HMAT) contains six sensorimotor regions within each hemisphere (Mayka et al., 2006). As the left M1 mask of the HMAT includes both GM and WM, we first also coregistered it to the brainmask in FreeSurfer space. This allowed us to compute average values of FA, NDI, ODI and IVF exclusively within the GM of left M1 as reconstructed in FreeSurfer.

Lesion segmentation

We estimated lesion volumes with the lesion growth algorithm (Schmidt et al., 2012) implemented in the lesion segmentation toolbox (LST) version 2.0.15 (https://www.applied-statistics.de/lst.html). T2-FLAIR volumes were coregistered to the T1-weighted images and segmentation information used to compute lesion belief maps. After visual inspection, κ = 0.2 was chosen for thresholding of these maps. By growing lesions along hypointensely appearing voxels in the T2-flair images, binary lesion maps were obtained. We linearly transformed T2-flair images to the brainmask obtained with FreeSurfer using FSL’s FLIRT. The transformation matrix was applied to transform lesion maps to FreeSurfer space. The percentage of lesion load in the left M1 mask was computed and used as covariate for regression analyses.

Statistical analyses

1. Microstructural correlates of motor excitability and cognitive-motor performance

Normality was assessed visually and using the Lilliefors test in R (version 3.6.3, RStudio). For group comparisons of motor threshold, FA, NDI, ODI and IVF, t-tests were conducted and considered significant in case of p < .05. Additionally, we tested homogeneity of variance between the groups using Levene’s test. Correlation coefficients were computed in R using Pearson correlation. We used R package corrplot for depiction of correlation coefficients (Figure 2 and S1). Hierarchical backwards regression was conducted using IBM SPSS Statistics for Windows, Version 23 (IBM Corp., Armonk, NY) for the MS and HC group separately. As control variables, we entered gender and age in the model for prediction of the dependent variable motor threshold. Additionally, the percentage of lesion load in left M1 was considered as control variable for the MS model. All control variables were entered first in one block in the model. Variables of interest were FA, NDI, ODI and IVF averaged within the left M1 GM mask. Similarly, regression with backwards elimination was conducted for 9HPT, TMT-A and TMT-B scores on FA, NDI, ODI and IVF in left M1 GM in MS patients. The order of removal of variables of interest was determined based on increase in correlation coefficients with the dependent variable (Figure 2 and S1), and predictors were only retained in case of p < .05 of the regression coefficient.

Figure 2.
Figure 2.

Correlation coefficients of the correlation between left M1 GM FA, NDI, ODI, IVF and motor threshold in (A) MS and (B) HC. * p < .05. Abbreviations: M1, primary motor cortex, GM, grey matter, FA, fractional anisotropy, NDI, neurite density index, ODI, orientation dispersion index, IVF, isotropic volume fraction, MS, multiple sclerosis, HC, healthy controls.

2. White matter microstructural integrity of the motor system

First, TBSS was conducted with FA maps of both groups. This method allows voxel-wise group comparisons by non-linear registration and projection onto an alignment-invariant tract representation (Smith et al., 2006). As target for the nonlinear registration, the FMRIB58_FA_1mm standard space template as provided by FSL was used. The skeleton was thresholded to κ = 0.2 (Smith et al., 2006). The resulting skeleton was then also applied to NDI, ODI and IVF maps. Permutation test using the randomize tool was applied for statistical analyses (Winkler et al., 2014). For each group comparison of FA, NDI, ODI and IVF between HC and MS, 500 permutations were carried out with the threshold-free cluster enhancement (TFCE) option to control for multiple compairsons. A threshold of p < .05 was used to define significance (Smith & Nichols, 2009). We overlaid all tract masks of the sensorimotor area tract template (SMATT) (Archer et al., 2018) on the mean FA skeleton and computed the number of voxels in the skeleton for each SMATT region. For creation of the SMATT, HMAT regions were used as seed region with waypoints placed at the level of the posterior limb of the internal capsule and the cerebral peduncle under exclusion of transcallosal fibers (Archer et al., 2018). We then extracted the number and percentage of voxels with significant group differences intersecting with SMATT tracts for FA, NDI, ODI and IVF maps, respectively. Lastly, we computed the percentage of overlapping group differences in all diffusion parameters per region and averages over all regions.

In Supplementary Results S5, we additionally provide results of a congruent analysis regarding Johns Hopkins University (JHU) DTI-based WM atlas that also covers tracts outside the motor system (Hua et al., 2008; Mori et al., 2005; Wakana et al., 2007).

Results

HC and MS were all right-handed and did not significantly differ in age or gender distribution (both p > .05).

1a. Microstructural correlates of motor cortex excitability

We first investigated group differences of motor threshold and FA, NDI, ODI and IVF in the left M1 GM between the MS and HC groups. Excitability levels, i.e. motor thresholds were non-significantly higher in MS patients compared to HC (t(97) = 0.78, p > .05; MS group: mean [SD] = 58.06 [10.60]; HC group: mean [SD] = 56.55 [8.63]). Group comparisons in FA, NDI, ODI and IVF in the GM of left M1 all did not reveal any significant group differences (all p > .05). Variability of NDI in MS patients was significantly higher than in HC (F(1,97) = 4.25, p = .04).

Within the GM of both MS and HC, diffusion measures were significantly interrelated (Figure 2). FA and NDI were positively interrelated, but negatively correlated to ODI and IVF. In MS patients, higher NDI was significantly linked to lower motor threshold; HC exhibited a similar trend.

We first entered age and gender as potentially confounding variables into the hierarchical backward regression model for the prediction of motor threshold for MS and HC separately. Predicting motor threshold in the MS group, we additionally included lesion load in M1. As variables of interest, we considered FA, NDI, ODI and IVF in the GM of left M1, which served as the stimulation site for determination of the motor threshold. NDI in the GM of left M1 significantly predicted motor threshold in MS (F(1,48) = 7.493, p = .009) and explained 11.7% of its variance (for details of hierarchical regression analysis see Table S2). Figure 3A depicts the relationship between the two variables and the respective density estimates, showing that lower neurite density in MS is linked to higher motor threshold, i.e. lower excitability. No model was significant in the HC group, however there was similarly a trend for NDI predicting motor threshold (Figure 3B, Table S3). By including NDI of both groups and a dummy variable encoding group membership for the prediction of motor threshold, we tested significance of group-specific slopes. The coefficient of the interaction term was not significant (p > .05), indicating that slopes were not significantly different between MS patients and HC.

Figure 3.
Figure 3.

Scatterplot with regression line and density estimates of the regression of motor threshold on neurite density index in (A) patients with MS and (B) HC. Abbreviations: HC, healthy controls, MS, multiple sclerosis.

1b. Microstructural correlates of cognitive-motor performance

MS patient’s 9HPT scores were regressed on control variables age, gender and lesion load in M1 and variables of interest FA, NDI, ODI and IVF in left M1 GM. No variable of interest significantly predicted 9HPT scores (Table S4).

Equivalently, TMT z-scores were regressed on control variables age, gender and lesion load in M1 and variables of interest FA, NDI, ODI and IVF in left M1 GM. NDI in the GM of left M1 was a significant predictor of both TMT-A (F(1,42) = 8.102; p =.007; adjusted r2 = .142) and TMT-B z-scores (F(1,42) = 7.390; p =.009; adjusted r2 = .129) (for details of hierarchical regression analyses see Table S5 and S6). Lower NDI predicted worse performance in both TMT-A (Figure 4A) and TMT-B (Figure 4B).

Figure 4.
Figure 4.

Scatterplot with regression line and density estimates of the regression of (A) TMT-A (z-scores) and (B) TMT-B (z-scores) on neurite density index in patients with MS. Abbreviations: MS, multiple sclerosis, TMT-A, Trail Making Test part A, TMT-B, Trail Making Test part B.

2. White matter microstructural integrity of the motor system

First, we compared FA, NDI, ODI and IVF maps between HC and MS patients using TBSS. TFCE-based significance tests showed significantly higher FA and NDI values in HC compared to MS in a large amount of motor WM tract voxels, i.e. 6.2% and 10.9% of the mean motor skeleton, respectively (p < .05). Conversely, ODI was higher in MS patients compared to HC in 3% of the mean motor skeleton voxels (p < .05) (Figure 5). No voxels showed significantly increased NDI or decreased ODI in MS compared to HC (p > .05). No significant group differences were observed regarding IVF (p > .05).

Figure 5.
Figure 5.

Results of tract-based spatial statistics analyses. Color bar indicates TFCE-corrected p-value, thresholded at p < .05. Contrast HC > MS for FA and NDI, and MS > HC for ODI. MNI coordinates: x = 90, y = 109, z = 95. Abbreviations: FA, fractional anisotropy, NDI, neurite density index, ODI, orientation dispersion index, L, left, R, right, TFCE, threshold-free cluster enhancement.

To obtain more clarity regarding the regional effects, we computed the number and percentage of voxels showing group differences for each region of the SMATT atlas for FA, NDI and ODI maps (Table 2). Fibers exhibiting group differences in at least 10% of the region in either FA, NDI or ODI included right supplementary motor area (SMA), pre-SMA, ventral premotor cortex (PMv) and bilateral somatosensory cortex (S1) and M1.

View this table:
Table 2.

Number and percentage of voxels significantly differing between MS and HC in FA, NDI and ODI in SMATT regions. The last column contains the total number of voxels of the mean skeleton intersecting with each SMATT region. Contrast HC > MS for FA and NDI, and MS > HC for ODI, sorted decreasingly by the percentage of voxels showing significant group differences of each region, summed over FA, NDI and ODI

Next, we disentangled which combination of FA, NDI and ODI, or if one of the parameters alone contributed to observed group differences. Figure 6A depicts the percentage of all combinations and single parameters for each SMATT region, while the percentage averaged over all regions is presented in Figure 6B. Neurite density exclusively (56%) and in combination with FA (19%) accounted for the largest amount of differences, followed by ODI alone (9%).

Figure 6.
Figure 6.

Comparison of diffusion parameters contributing to observed group differences. (A) Percentage of FA, NDI and ODI separately or combined showing group differences in the tract skeleton in intersection with SMATT regions. (B) Percentage averaged over regions with any significant group difference (B). Abbreviations: FA, fractional anisotropy, NDI, neurite density index, ODI, orientation dispersion index, SMA, supplementary motor area, M1, primary motor cortex, S1, primary somatosensory cortex, PMv, ventral premotor cortex. PMd, dorsal premotor cortex.

Discussion

In order to link microstructural integrity and cortical excitability, we applied a biophysical diffusion model on the basis of a multi-shell acquisition protocol. This allowed us to obtain maps of NDI, ODI and IVF in addition to the classical DTI parameter FA within the WM and GM of the brain motor regions. We linked lower microstructural integrity in the GM to lower motor cortical excitability and lower cognitive-motor performance in MS. In the WM, we observed lower FA and NDI and higher ODI in MS patients compared to HC in all motor tracts. Voxel-wise group differences in FA strongly co-occurred with NDI differences. NDI and ODI appear sensitive for detection of neuropathological changes in MS patients undetected by FA.

1a. Microstructural correlates of motor excitability

Cortical excitability can be altered in MS patients following neuronal and axonal loss (Calabrese et al., 2007; Zipser et al., 2018). To unravel the link between microstructure and excitability, we applied NODDI, a biophysical model of microstructure. We postulate a tight excitability-microstructural integrity relationship that strengthens during the disease course. A robust quantification of this relationship could predict functional impairment and be extended to monitor the disease course or therapeutic responses.

In contrast to classical DTI-based measures such as FA that are biased in areas of pronounced neurite orientation dispersion such as the GM, this tendency is mitigated in the NODDI model that is inherently more specific in depicting the cortical microstructure (Zhang et al., 2012). We were interested in how motor cortical integrity quantified by NODDI parameters explains altered excitability as evaluated by the resting motor threshold, a unique non-invasive approach cortex excitability assessments in humans. An initial comparison between the MS and HC groups revealed that motor threshold exhibited a trending increase in MS compared to HC. Higher motor thresholds, i.e. lower excitability, in MS than HC have been reported previously, particularly for the secondary progressive disease type and in patients in the relapsing phase (Caramia et al., 2004; Vucic et al., 2012). Decreased motor cortical excitability in MS patients potentially is a consequence of neuronal loss, axonal scarcity or synaptic down-scaling of the cortico-cortical inputs (Boorman et al., 2007; Calabrese et al., 2007; Wahl et al., 2007; Zipser et al., 2018). The motor thresholds did not statistically differ between the addressed groups in our study, however showed a trending increase in MS. All diffusion measures did not significantly depict group differences, but in line with a previous study, trending lower ODI in MS compared to HC was observable (De Santis et al., 2019). A decrease in the GM ODI could indicate degeneration of single neurites, resulting in dispersion reductions (Spano et al., 2018; Zhang et al., 2012). Presumably, both a larger sample size and a longer disease duration would amplify these trending differences in motor threshold and NODDI parameters between MS and HC.

We next examined the relation of FA, NDI, ODI and IVF among each other within M1, where the stimulation for assessment of the motor threshold was conducted. FA and NDI were significantly positively interrelated, and negatively linked to ODI and IVF. A hierarchical backward regression model revealed that lower NDI predicted higher motor threshold, i.e. lower neural excitability in MS patients. This is well in line with the suggestion of axonal and myelin loss due to neuroinflammatory processes (Filippi et al., 2012; Grussu et al., 2017). Furthermore, this shows the superiority of NDI in predicting excitability over local lesion load that was included in the model as a potential confounder. The HC model was not significant, however a trend of a positive relation between NDI and motor threshold was observable here as well. We suggest that higher scattering of NDI led to a better model fit in the MS group as higher variance allows stronger correlations, while the model in general also fits to the HC group. This assumption could be tested by involving HC with more strongly scattered values, such as by including individuals of higher age that typically show higher resting motor thresholds and increasing neuronal loss (Bhandari et al., 2016; Ding et al., 2016; Merluzzi et al., 2016).

1b. Microstructural correlates of cognitive-motor performance

We further assessed if microstructural integrity of the left M1 GM can also be linked to a behavioral outcome, i.e. cognitive-motor performance in the MS group. Hierarchical regression results indicated that patients with a loss in neurite density showed a lower cognitive-motor performance as quantified by TMT-A and TMT-B. For 9HPT, the regression model did not include significant predictors of interest. However, fewer patients completed this test and a trend for a positive correlation was detected between NDI and 9HPT scores (Figure S1). Previous results demonstrated a relation between lower cognitive performance and lower GM volume in regions relevant for demands of the utilized neuropsychological test (Morgen et al., 2006), but also the whole brain (Fisniku et al., 2008; Sastre-Garriga et al., 2004). Lower magnetization transfer ratio was observed in motor regions for patients with lower 9HPT performance, likely reflecting pathological processes including inflammation, edema, demyelination and axonal loss (Khaleeli et al., 2007; Schmierer et al., 2004; van Waesberghe et al., 1999). The TMT measures multiple cognitive domains such as visuo-motor abilities and cognitive flexibility and it has been shown that precentral gyrus was activated in an fMRI-compatible TMT adaptation (Crowe, 1998; Zakzanis et al., 2005). Hence, a lowered TMT performance linked to lower NDI in left M1 as seen here appears consequential and underlines the suitability of the NODDI model for assessment of GM microstructural integrity.

2. White matter microstructural integrity of the motor system

To gain a complete depiction of motor-relevant microstructure and its pathological alterations in MS, we next investigated WM motor tracts by applying TBSS. We detected decreased FA and NDI and increased ODI in patients with MS compared to HC, where the highest percentage of voxels showing group differences was in NDI. The strongest portion of sensorimotor tracts in intersection with group differences in the TBSS FA mean skeleton included WM tracts originating in right SMA and PMv. NDI was particularly lower in bilateral M1 tracts, whereas ODI was higher mainly in right SMA and S1 tracts of MS patients. GM atrophy in MS compared to HC in precentral and postcentral regions has been reported previously (Lansley et al., 2013; Prinster et al., 2010). In a recent longitudinal study, M1 and SMA atrophy was more pronounced in patients with disability progression compared to those without after 10 years (Bergsland et al., 2018). These findings make it plausible that connecting WM tracts show abnormalities due to primary or secondary causes as well, which has been demonstrated earlier (Steenwijk et al., 2015). Also in support of this view, lower FA was observed in the corticospinal tract of MS patients, and higher mean, axial and radial diffusivity was linked to lower M1 thickness (Bergsland et al., 2015). To the best of our knowledge, NODDI has not yet been investigated systematically in sensorimotor tracts of MS patients as was done in the present work. However, a small sample-sized first NODDI study also reported decreased NDI and increased ODI in the normal-appearing WM of MS patients (Schneider et al., 2017). Higher ODI indicates fiber coherence loss and can possibly be explained by an increase in compensative axonal sprouting or branching (Timmers et al., 2015). In a recent histological study, ODI was shown to indeed reflect the neurite architecture complexity, whereas reduced NDI indicates loss of myelin and axons (Grussu et al., 2017). Regarding NDI, a decrease in MS is expected due to the demyelinating, inflammatory and neuro-axonal pathology of the disease (Filippi et al., 2012). A NODDI study that performed voxel-wise comparisons revealed reduced NDI and increased ODI in the fornix and corpus callosum of RRMS compared to HC (Spano et al., 2018). Our supplementary results applying the same analysis to JHU regions similarly yielded reduced NDI and increased ODI in MS compared to HC in these structures.

Importantly, low FA values could be caused by lower NDI, higher ODI, or a combination of both (Zhang et al., 2012), such that it is of great relevance to disentangle if group effects appear for these parameters exclusively in a voxel, or if they potentially overlap, which has not yet been systematically quantified in MS research to our knowledge. This allows a more concrete inference on the type of microstructure change in contrast to assumptions solely based on FA. On average, 56 % of voxels showing significant group differences in intersection with the mean FA skeleton were different exclusively in NDI, whereas 19 % were different simultaneously in FA and NDI, 9 % in ODI only and 7 % in FA only. This pattern was consistent among single sensorimotor WM tracts. Our findings suggest that FA is capable of detecting parts of the pathological alterations that can mainly be retraced to NDI decreases in MS. The NODDI model therefore allows to capture further microstructure changes in MS in addition to its inherently higher specificity to the underlying pathology (Zhang et al., 2012). These TBSS results support previous findings based on voxel-wise comparisons (Schneider et al., 2017; Spano et al., 2018) and further disentangle that the lower FA signal is predominantly caused by reduced neurite density in sensorimotor WM tracts of MS patients.

Limitations and prospect

We provide novel insights regarding the applicability of the NODDI model to assess microstructure in both GM and WM. Our suggestion that neurite density is also linked to motor cortical excitability in healthy participants needs to be further examined in future studies including a subject group with a higher variability in age. Furthermore, not all MS patients completed the 9HPT, which likely lowered statistical power. The link between neuropsychological test performance and NODDI parameters could in future be expanded to other cognitive domains and brain regions analyzing larger patient samples. It would be particularly interesting to observe potential effects of motor and cognitive training or medication on beneficial adaptations in microstructural integrity (Bonzano et al., 2014; Moccia et al., 2017). Due to the multiple facets of the disorder, including the specific disease type and duration, treatment, relapse rate and localization of pathologically affected regions, investigating larger sample sizes that allow stratification to control these variations is of importance. How microstructural state as evaluated by biophysical diffusion models is connected to functional adaptive and maladaptive mechanisms of brain networks would be an interesting question to revisit using functional MRI (Droby et al., 2016; Fleischer et al., 2020; Fleischer, Radetz, et al., 2019). Moreover, the integration of further structural imaging methods such as relaxation times and other diffusion models is relevant for in vivo delineation of pathological microstructure changes as are histological assessments in ex vivo samples (Alexander et al., 2019; Gracien et al., 2017; Grussu et al., 2017; Lakhani et al., 2020). A translational approach linking NODDI in different animal models of MS to human patients remains to be further pursued as was done with a cuprizone model of demyelination in mice and by applying DTI (Cerina et al., 2020; Krämer et al., 2019; Wang et al., 2019). Furthermore, differentiating the evolution of microstructure in MS compared to healthily aging individuals is of interest in longitudinal investigations (Fleischer, Koirala, et al., 2019).

Conclusion

In this work we applied a biophysical model of GM and WM microstructure in the motor system of patients with MS and age-matched HC. We hereby demonstrate that lower neurite density in left M1 is linked to decreased motor cortical excitability in patients with MS. Also on the behavioral level, decreased neurite density is indicative of impaired cognitive-motor performance. Furthermore, we showed that lower neurite density and higher orientation dispersion are characteristic in the WM of MS patients compared to HC, and that these markers are more sensitive to pathological alterations than the classical DTI measure FA. These findings suggest that advanced biophysical diffusion models are of great relevance for prediction of neurodegenerative processes, excitability alterations and disease progression.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Competing interests

S. G. Meuth has received honoraria for lecturing and travel expenses for attending meetings from Almirall, Amicus Therapeutics Germany, Bayer Health Care, Biogen, Celgene, Diamed, Genzyme, MedDay Pharmaceuticals, Merck Serono, Novartis, Novo Nordisk, ONO Pharma, Roche, Sanofi-Aventis, Chugai Pharma, QuintilesIMS, and Teva. His research is funded by the German Ministry for Education and Research (BMBF), Deutsche Forschungsgemeinschaft (DFG), Else Kröner Fresenius Foundation, German Academic Exchange Service, Hertie Foundation, Interdisciplinary Center for Clinical Studies (IZKF) Muenster, German Foundation Neurology, and by Almirall, Amicus Therapeutics Germany, Biogen, Diamed, Fresenius Medical Care, Genzyme, Merck Serono, Novartis, ONO Pharma, Roche, and Teva. J. Krämer received honoraria for lecturing from Biogen, Novartis, Merck Serono, Sanofi-Genzyme, Roche, Mylan and Teva, and financial research support from Sanofi Genzyme. S. Bittner has received honoraria and compensation for travel from Biogen Idec, Merck Serono, Novartis, Sanofi-Genzyme and Roche.

Acknowledgements

We thank all participants who took part in the study and Cheryl Ernest for proofreading the manuscript. This study was supported by the German Research Council (DFG, CRC TR-128 to S.B., S.G.M., S.G.). Parts of this research were conducted using the supercomputer Mogon and advisory services offered by Johannes Gutenberg University Mainz (hpc.uni-mainz.de), which is a member of the AHRP and the Gauss Alliance e.V. We gratefully acknowledge the computing time granted on Mogon.

Abbreviations

APB
Abductor pollicis brevis
AMICO
Accelerated Microstructure Imaging via Convex Optimization
DTI
diffusion tensor imaging
DWI
diffusion-weighted imaging
PMd
dorsal premotor cortex
TE
echo time
EMG
electromyography
EDSS
Expanded Disability Status Scale
FOV
field of view
FDI
first dorsal interosseous
FLAIR
fluid attenuated inversion recovery
FA
fractional anisotropy
GM
grey matter
HC
healthy controls
HMAT
Human Motor Area Template
TI
inversion time
IVF
isotropic volume fraction
JHU
Johns Hopkins University
L
left
LST
lesion segmentation toolbox
MEP
motor evoked potential
MS
multiple sclerosis
NDI
neurite density index
NODDI
neurite orientation dispersion and density imaging
9HPT
Nine Hole Peg Test
ODI
orientation dispersion index
M1
primary motor cortex
S1
primary somatosensory cortex
TR
repetition time
RRMS
relapsing-remitting MS
R
right
SMATT
sensorimotor area tract template
SD
standard deviation
SMA
supplementary motor area
TFCE
threshold-free cluster enhancement
TBSS
tract-based spatial statistics
TMT-A
Trail Making Test part A
TMT-B
Trail Making Test part B
TMS
transcranial magnetic stimulation
PMv
ventral premotor cortex
WM
white matter

References

  1. ↵
    Alexander, D. C., Dyrby, T. B., Nilsson, M., & Zhang, H. (2019). Imaging brain microstructure with diffusion MRI: Practicality and applications. NMR in Biomedicine, 32(4), e3841. https://doi.org/10.1002/nbm.3841
    OpenUrl
  2. ↵
    Archer, D. B., Vaillancourt, D. E., & Coombes, S. A. (2018). A template and probabilistic atlas of the human sensorimotor tracts using diffusion MRI. Cerebral Cortex, 28(5), 1685–1699. https://doi.org/10.1093/cercor/bhx066
    OpenUrl
  3. ↵
    Bergsland, N., Horakova, D., Dwyer, M. G., Uher, T., Vaneckova, M., Tyblova, M., Seidl, Z., Krasensky, J., Havrdova, E., & Zivadinov, R. (2018). Gray matter atrophy patterns in multiple sclerosis: A 10-year source-based morphometry study. NeuroImage: Clinical, 17, 444–451. https://doi.org/10.1016/j.nicl.2017.11.002
    OpenUrl
  4. ↵
    Bergsland, N., Laganà, M. M., Tavazzi, E., Caffini, M., Tortorella, P., Baglio, F., Baselli, G., & Rovaris, M. (2015). Corticospinal tract integrity is related to primary motor cortex thinning in relapsing–remitting multiple sclerosis. Multiple Sclerosis Journal, 21(14), 1771–1780. https://doi.org/10.1177/1352458515576985
    OpenUrl
  5. ↵
    Bhandari, A., Radhu, N., Farzan, F., Mulsant, B. H., Rajji, T. K., Daskalakis, Z. J., & Blumberger, D. M. (2016). A meta-analysis of the effects of aging on motor cortex neurophysiology assessed by transcranial magnetic stimulation. Clinical Neurophysiology, 127(8), 2834–2845. https://doi.org/10.1016/j.clinph.2016.05.363
    OpenUrlCrossRefPubMed
  6. ↵
    Bonzano, L., Tacchino, A., Brichetto, G., Roccatagliata, L., Dessypris, A., Feraco, P., Lopes De Carvalho, M. L., Battaglia, M. A., Mancardi, G. L., & Bove, M. (2014). Upper limb motor rehabilitation impacts white matter microstructure in multiple sclerosis. NeuroImage, 90, 107–116. https://doi.org/10.1016/j.neuroimage.2013.12.025
    OpenUrlCrossRefPubMed
  7. ↵
    Boorman, E. D., O’Shea, J., Sebastian, C., Rushworth, M. F. S., & Johansen-Berg, H. (2007). Individual Differences in White-Matter Microstructure Reflect Variation in Functional Connectivity during Choice. Current Biology, 17(16), 1426–1431. https://doi.org/10.1016/j.cub.2007.07.040
    OpenUrlCrossRefPubMedWeb of Science
  8. ↵
    Calabrese, M., Atzori, M., Bernardi, V., Morra, A., Romualdi, C., Rinaldi, L., McAuliffe, M. J. M., Barachino, L., Perini, P., Fischl, B., Battistin, L., & Gallo, P. (2007). Cortical atrophy is relevant in multiple sclerosis at clinical onset. Journal of Neurology, 254(9), 1212–1220. https://doi.org/10.1007/s00415-006-0503-6
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    Caramia, M. D., Palmieri, M. G., Desiato, M. T., Boffa, L., Galizia, P., Rossini, P. M., Centonze, D., & Bernardi, G. (2004). Brain excitability changes in the relapsing and remitting phases of multiple sclerosis: A study with transcranial magnetic stimulation. Clinical Neurophysiology, 115(4), 956–965. https://doi.org/10.1016/j.clinph.2003.11.024
    OpenUrlCrossRefPubMedWeb of Science
  10. ↵
    Cerina, M., Muthuraman, M., Gallus, M., Koirala, N., Dik, A., Wachsmuth, L., Hundehege, P., Schiffler, P., Tenberge, J.-G., Fleischer, V., Gonzalez-Escamilla, G., Narayanan, V., Krämer, J., Faber, C., Budde, T., Groppa, S., & Meuth, S. G. (2020). Myelination- and immune-mediated MR-based brain network correlates. Journal of Neuroinflammation, 17(1), 186. https://doi.org/10.1186/s12974-020-01827-z
    OpenUrl
  11. ↵
    Crowe, S. F. (1998). The differential contribution of mental tracking, cognitive flexibility, visual search, and motor speed to performance on parts A and B of the Trail Making Test. Journal of clinical psychology, 54(5), 585–591. https://doi.org/10.1002/(SICI)1097-4679(199808)54:5<585::AID-JCLP4>3.0.CO;2-K
    OpenUrlCrossRefPubMedWeb of Science
  12. ↵
    Daducci, A., Canales-Rodríguez, E. J., Zhang, H., Dyrby, T. B., Alexander, D. C., & Thiran, J.-P. (2015). Accelerated microstructure imaging via convex optimization (AMICO) from diffusion MRI data. NeuroImage, 105, 32–44. https://doi.org/10.1016/j.neuroimage.2014.10.026
    OpenUrlCrossRefPubMed
  13. ↵
    De Santis, S., Bastiani, M., Droby, A., Kolber, P., Zipp, F., Pracht, E., Stoecker, T., Groppa, S., & Roebroeck, A. (2019). Characterizing microstructural tissue properties in multiple sclerosis with diffusion MRI at 7 T and 3 T: The impact of the experimental design. Neuroscience, 403, 17–26. https://doi.org/10.1016/j.neuroscience.2018.03.048
    OpenUrlCrossRef
  14. ↵
    Deppe, M., Krämer, J., Tenberge, J.-G., Marinell, J., Schwindt, W., Deppe, K., Groppa, S., Wiendl, H., & Meuth, S. G. (2016). Early silent microstructural degeneration and atrophy of the thalamocortical network in multiple sclerosis. Human Brain Mapping, 37(5), 1866–1879. https://doi.org/10.1002/hbm.23144
    OpenUrlCrossRefPubMed
  15. ↵
    Deppe, M., Tabelow, K., Krämer, J., Tenberge, J.-G., Schiffler, P., Bittner, S., Schwindt, W., Zipp, F., Wiendl, H., & Meuth, S. G. (2016). Evidence for early, non-lesional cerebellar damage in patients with multiple sclerosis: DTI measures correlate with disability, atrophy, and disease duration. Multiple Sclerosis Journal, 22(1), 73–84. https://doi.org/10.1177/1352458515579439
    OpenUrl
  16. ↵
    Ding, X.-Q., Maudsley, A. A., Sabati, M., Sheriff, S., Schmitz, B., Schütze, M., Bronzlik, P., Kahl, K. G., & Lanfermann, H. (2016). Physiological neuronal decline in healthy aging human brain—An in vivo study with MRI and short echo-time whole-brain 1H MR spectroscopic imaging. NeuroImage, 137, 45–51. https://doi.org/10.1016/j.neuroimage.2016.05.014
    OpenUrl
  17. ↵
    Droby, A., Fleischer, V., Carnini, M., Zimmermann, H., Siffrin, V., Gawehn, J., Erb, M., Hildebrandt, A., Baier, B., & Zipp, F. (2015). The impact of isolated lesions on white-matter fiber tracts in multiple sclerosis patients. NeuroImage. Clinical, 8, 110–116. https://doi.org/10.1016/j.nicl.2015.03.003
    OpenUrl
  18. ↵
    Droby, A., Yuen, K. S. L., Muthuraman, M., Reitz, S.-C., Fleischer, V., Klein, J., Gracien, R.-M., Ziemann, U., Deichmann, R., Zipp, F., & Groppa, S. (2016). Changes in brain functional connectivity patterns are driven by an individual lesion in MS: a resting-state fMRI study. Brain Imaging and Behavior, 10(4), 1117–1126. https://doi.org/10.1007/s11682-015-9476-3
    OpenUrl
  19. ↵
    Ellwardt, E., Pramanik, G., Luchtman, D., Novkovic, T., Jubal, E. R., Vogt, J., Arnoux, I., Vogelaar, C. F., Mandal, S., Schmalz, M., Barger, Z., Ruiz de Azua, I., Kuhlmann, T., Lutz, B., Mittmann, T., Bittner, S., Zipp, F., & Stroh, A. (2018). Maladaptive cortical hyperactivity upon recovery from experimental autoimmune encephalomyelitis. Nature Neuroscience, 21(10), 1392–1403. https://doi.org/10.1038/s41593-018-0193-2
    OpenUrl
  20. ↵
    Eshaghi, A., Prados, F., Brownlee, W. J., Altmann, D. R., Tur, C., Cardoso, M. J., De Angelis, F., van de Pavert, S. H., Cawley, N., De Stefano, N., Stromillo, M. L., Battaglini, M., Ruggieri, S., Gasperini, C., Filippi, M., Rocca, M. A., Rovira, A., Sastre-Garriga, J., Vrenken, H., … MAGNIMS study group. (2018). Deep gray matter volume loss drives disability worsening in multiple sclerosis. Annals of Neurology, 83(2), 210–222. https://doi.org/10.1002/ana.25145
    OpenUrlCrossRefPubMed
  21. ↵
    Feys, P., Lamers, I., Francis, G., Benedict, R., Phillips, G., LaRocca, N., Hudson, L. D., Rudick, R., & Multiple Sclerosis Outcome Assessments Consortium. (2017). The Nine-Hole Peg Test as a manual dexterity performance measure for multiple sclerosis. Multiple Sclerosis Journal, 23(5), 711–720. https://doi.org/10.1177/1352458517690824
    OpenUrl
  22. ↵
    Filippi, M., Rocca, M. A., Barkhof, F., Brück, W., Chen, J. T., Comi, G., DeLuca, G., De Stefano, N., Erickson, B. J., Evangelou, N., Fazekas, F., Geurts, J. J. G., Lucchinetti, C., Miller, D. H., Pelletier, D., Popescu, B. F. G., Lassmann, H., & Attendees of the Correlation between Pathological MRI findings in MS workshop. (2012). Association between pathological and MRI findings in multiple sclerosis. The Lancet. Neurology, 11(4), 349–360. https://doi.org/10.1016/S1474-4422(12)70003-0
    OpenUrl
  23. ↵
    Fischl, B. (2012). FreeSurfer. NeuroImage, 62(2), 774–781. https://doi.org/10.1016/j.neuroimage.2012.01.021
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    Fisniku, L. K., Chard, D. T., Jackson, J. S., Anderson, V. M., Altmann, D. R., Miszkiel, K. A., Thompson, A. J., & Miller, D. H. (2008). Gray matter atrophy is related to long-term disability in multiple sclerosis. Annals of Neurology, 64(3), 247–254. https://doi.org/10.1002/ana.21423
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    Fleischer, V., Koirala, N., Droby, A., Gracien, R.-M., Deichmann, R., Ziemann, U., Meuth, S. G., Muthuraman, M., Zipp, F., & Groppa, S. (2019). Longitudinal cortical network reorganization in early relapsing–remitting multiple sclerosis. Therapeutic advances in neurological disorders, 12, 1–15. https://doi.org/10.1177/1756286419838673
    OpenUrl
  26. ↵
    Fleischer, V., Muthuraman, M., Anwar, A. R., Gonzalez-Escamilla, G., Radetz, A., Gracien, R.-M., Bittner, S., Luessi, F., Meuth, S. G., Zipp, F., & Groppa, S. (2020). Continuous reorganization of cortical information flow in multiple sclerosis: A longitudinal fMRI effective connectivity study. Scientific Reports, 10(1), 806. https://doi.org/10.1038/s41598-020-57895-x
    OpenUrl
  27. ↵
    Fleischer, V., Radetz, A., Ciolac, D., Muthuraman, M., Gonzalez-Escamilla, G., Zipp, F., & Groppa, S. (2019). Graph theoretical framework of brain networks in multiple sclerosis: A review of concepts. Neuroscience, 403, 35–53. https://doi.org/10.1016/j.neuroscience.2017.10.033
    OpenUrl
  28. ↵
    Fukutomi, H., Glasser, M. F., Zhang, H., Autio, J. A., Coalson, T. S., Okada, T., Togashi, K., Van Essen, D. C., & Hayashi, T. (2018). Neurite imaging reveals microstructural variations in human cerebral cortical gray matter. NeuroImage, 182, 488–499. https://doi.org/10.1016/j.neuroimage.2018.02.017
    OpenUrl
  29. ↵
    Gracien, R.-M., Reitz, S. C., Hof, S.-M., Fleischer, V., Droby, A., Wahl, M., Steinmetz, H., Groppa, S., Deichmann, R., & Klein, J. C. (2017). Longitudinal quantitative MRI assessment of cortical damage in multiple sclerosis: A pilot study. Journal of Magnetic Resonance Imaging, 46(5), 1485–1490. https://doi.org/10.1002/jmri.25685
    OpenUrl
  30. ↵
    Granberg, T., Fan, Q., Treaba, C. A., Ouellette, R., Herranz, E., Mangeat, G., Louapre, C., Cohen-Adad, J., Klawiter, E. C., Sloane, J. A., & Mainero, C. (2017). In vivo characterization of cortical and white matter neuroaxonal pathology in early multiple sclerosis. Brain: A Journal of Neurology, 140(11), 2912–2926. https://doi.org/10.1093/brain/awx247
    OpenUrlCrossRefPubMed
  31. ↵
    Groppa, S., Werner-Petroll, N., Münchau, A., Deuschl, G., Ruschworth, M. F., & Siebner, H. R. (2012). A novel dual-site transcranial magnetic stimulation paradigm to probe fast facilitatory inputs from ipsilateral dorsal premotor cortex to primary motor cortex. NeuroImage, 62(1), 500–509. https://doi.org/10.1016/j.neuroimage.2012.05.023
    OpenUrlCrossRefPubMedWeb of Science
  32. ↵
    Grussu, F., Schneider, T., Tur, C., Yates, R. L., Tachrount, M., Ianuş, A., Yiannakas, M. C., Newcombe, J., Zhang, H., Alexander, D. C., DeLuca, G. C., & Gandini Wheeler-Kingshott, C. A. M. (2017). Neurite dispersion: A new marker of multiple sclerosis spinal cord pathology? Annals of Clinical and Translational Neurology, 4(9), 663–679. https://doi.org/10.1002/acn3.445
    OpenUrl
  33. ↵
    Hua, K., Zhang, J., Wakana, S., Jiang, H., Li, X., Reich, D. S., Calabresi, P. A., Pekar, J. J., van Zijl, P. C., & Mori, S. (2008). Tract probability maps in stereotaxic spaces: Analyses of white matter anatomy and tract-specific quantification. NeuroImage, 39(1), 336–347. https://doi.org/10.1016/j.neuroimage.2007.07.053
    OpenUrlCrossRefPubMedWeb of Science
  34. ↵
    Khaleeli, Z., Cercignani, M., Audoin, B., Ciccarelli, O., Miller, D. H., & Thompson, A. J. (2007). Localized grey matter damage in early primary progressive multiple sclerosis contributes to disability. NeuroImage, 37(1), 253–261. https://doi.org/10.1016/j.neuroimage.2007.04.056
    OpenUrlCrossRefPubMedWeb of Science
  35. ↵
    Krämer, J., Brück, W., Zipp, F., Cerina, M., Groppa, S., & Meuth, S. G. (2019). Imaging in mice and men: Pathophysiological insights into multiple sclerosis from conventional and advanced MRI techniques. Progress in Neurobiology, 182, 101663. https://doi.org/10.1016/j.pneurobio.2019.101663
    OpenUrl
  36. ↵
    Kurtzke, J. F. (1983). Rating neurologic impairment in multiple sclerosis: An expanded disability status scale (EDSS). Neurology, 33(11), 1444–1444. https://doi.org/10.1212/WNL.33.11.1444
    OpenUrlCrossRefPubMed
  37. ↵
    Lakhani, D. A., Schilling, K. G., Xu, J., & Bagnato, F. (2020). Advanced Multicompartment Diffusion MRI Models and Their Application in Multiple Sclerosis. American Journal of Neuroradiology, 41(5), 751–757. https://doi.org/10.3174/ajnr.A6484
    OpenUrlAbstract/FREE Full Text
  38. ↵
    Lansley, J., Mataix-Cols, D., Grau, M., Radua, J., & Sastre-Garriga, J. (2013). Localized grey matter atrophy in multiple sclerosis: A meta-analysis of voxel-based morphometry studies and associations with functional disability. Neuroscience & Biobehavioral Reviews, 37(5), 819–830. https://doi.org/10.1016/j.neubiorev.2013.03.006
    OpenUrl
  39. ↵
    Mayka, M. A., Corcos, D. M., Leurgans, S. E., & Vaillancourt, D. E. (2006). Three-dimensional locations and boundaries of motor and premotor cortices as defined by functional brain imaging: A meta-analysis. NeuroImage, 31(4), 1453–1474. https://doi.org/10.1016/j.neuroimage.2006.02.004
    OpenUrlCrossRefPubMedWeb of Science
  40. ↵
    Merluzzi, A. P., Dean, D. C., Adluru, N., Suryawanshi, G. S., Okonkwo, O. C., Oh, J. M., Hermann, B. P., Sager, M. A., Asthana, S., Zhang, H., Johnson, S. C., Alexander, A. L., & Bendlin, B. B. (2016). Age-dependent differences in brain tissue microstructure assessed with neurite orientation dispersion and density imaging. Neurobiology of Aging, 43, 79–88. https://doi.org/10.1016/j.neurobiolaging.2016.03.026
    OpenUrlCrossRef
  41. ↵
    Moccia, M., de Stefano, N., & Barkhof, F. (2017). Imaging outcome measures for progressive multiple sclerosis trials. Multiple Sclerosis Journal, 23(12), 1614–1626. https://doi.org/10.1177/1352458517729456
    OpenUrl
  42. ↵
    Morgen, K., Sammer, G., Courtney, S. M., Wolters, T., Melchior, H., Blecker, C. R., Oschmann, P., Kaps, M., & Vaitl, D. (2006). Evidence for a direct association between cortical atrophy and cognitive impairment in relapsing–remitting MS. NeuroImage, 30(3), 891–898. https://doi.org/10.1016/j.neuroimage.2005.10.032
    OpenUrlCrossRefPubMedWeb of Science
  43. ↵
    Mori, S., Wakana, S., Van Zijl, P. C., & Nagae-Poetscher, L. M. (2005). MRI atlas of human white matter. Elsevier Science.
  44. ↵
    Prinster, A., Quarantelli, M., Lanzillo, R., Orefice, G., Vacca, G., Carotenuto, B., Alfano, B., Brunetti, A., Brescia Morra, V., & Salvatore, M. (2010). A voxel-based morphometry study of disease severity correlates in relapsing—Remitting multiple sclerosis. Multiple Sclerosis Journal, 16(1), 45–54. https://doi.org/10.1177/1352458509351896
    OpenUrl
  45. ↵
    Radetz, A., Koirala, N., Krämer, J., Johnen, A., Fleischer, V., Gonzalez-Escamilla, G., Cerina, M., Muthuraman, M., Meuth, S. G., & Groppa, S. (2019). Gray matter integrity predicts white matter network reorganization in multiple sclerosis. Human Brain Mapping. https://doi.org/10.1002/hbm.24849
  46. ↵
    Reitan, R. M. (1971). Trail making test results for normal and brain-damaged children. Perceptual and Motor Skills, 33(2), 575–581. https://doi.org/10.2466/pms.1971.33.2.575
    OpenUrlCrossRefPubMedWeb of Science
  47. ↵
    Roosendaal, S. D., Geurts, J. J., Vrenken, H., Hulst, H. E., Cover, K. S., Castelijns, J. A., Pouwels, P. J., & Barkhof, F. (2009). Regional DTI differences in multiple sclerosis patients. NeuroImage, 44(4), 1397–1403. https://doi.org/10.1016/j.neuroimage.2008.10.026
    OpenUrlCrossRefPubMedWeb of Science
  48. ↵
    Sastre-Garriga, J., Ingle, G. T., Chard, D. T., Ramió-Torrentà, L., Miller, D. H., & Thompson, A. J. (2004). Grey and white matter atrophy in early clinical stages of primary progressive multiple sclerosis. NeuroImage, 22(1), 353–359. https://doi.org/10.1016/j.neuroimage.2004.02.008
    OpenUrlCrossRefPubMedWeb of Science
  49. ↵
    Schmidt, P., Gaser, C., Arsic, M., Buck, D., Förschler, A., Berthele, A., Hoshi, M., Ilg, R., Schmid, V. J., Zimmer, C., Hemmer, B., & Mühlau, M. (2012). An automated tool for detection of FLAIR-hyperintense white-matter lesions in Multiple Sclerosis. NeuroImage, 59(4), 3774–3783. https://doi.org/10.1016/j.neuroimage.2011.11.032
    OpenUrlCrossRefPubMedWeb of Science
  50. ↵
    Schmierer, K., Scaravilli, F., Altmann, D. R., Barker, G. J., & Miller, D. H. (2004). Magnetization transfer ratio and myelin in postmortem multiple sclerosis brain. Annals of neurology, 56(3), 407–415. https://doi.org/10.1002/ana.20202
    OpenUrlCrossRefPubMedWeb of Science
  51. ↵
    Schneider, T., Brownlee, W., Zhang, H., Ciccarelli, O., Miller, D. H., & Wheeler-Kingshott, C. G. (2017). Sensitivity of multi-shell NODDI to multiple sclerosis white matter changes: A pilot study. Functional neurology, 32(2), 97. https://doi.org/10.11138/FNeur/2017.32.2.097
    OpenUrlCrossRefPubMed
  52. ↵
    Smith, S. M., Jenkinson, M., Johansen-Berg, H., Rueckert, D., Nichols, T. E., Mackay, C. E., Watkins, K. E., Ciccarelli, O., Cader, M. Z., Matthews, P. M., & Behrens, T. E. J. (2006). Tract-based spatial statistics: Voxelwise analysis of multi-subject diffusion data. NeuroImage, 31(4), 1487–1505. https://doi.org/10.1016/j.neuroimage.2006.02.024
    OpenUrlCrossRefPubMedWeb of Science
  53. ↵
    Smith, S. M., & Nichols, T. E. (2009). Threshold-free cluster enhancement: Addressing problems of smoothing, threshold dependence and localisation in cluster inference. NeuroImage, 44(1), 83–98. https://doi.org/10.1016/j.neuroimage.2008.03.061
    OpenUrlCrossRefPubMedWeb of Science
  54. ↵
    Spano, B., Giulietti, G., Pisani, V., Morreale, M., Tuzzi, E., Nocentini, U., Francia, A., Caltagirone, C., Bozzali, M., & Cercignani, M. (2018). Disruption of neurite morphology parallels MS progression. Neurology Neuroimmunology Neuroinflammation, 5(6), e502. https://doi.org/10.1212/NXI.0000000000000502
    OpenUrl
  55. ↵
    Stampanoni Bassi, M., Buttari, F., Maffei, P., De Paolis, N., Sancesario, A., Gilio, L., Pavone, L., Pasqua, G., Simonelli, I., Sica, F., Fantozzi, R., Bellantonio, P., Centonze, D., & Iezzi, E. (2020). Practice-dependent motor cortex plasticity is reduced in non-disabled multiple sclerosis patients. Clinical Neurophysiology, 131(2), 566–573. https://doi.org/10.1016/j.clinph.2019.10.023
    OpenUrl
  56. ↵
    Steenwijk, M. D., Daams, M., Pouwels, P. J.J., Balk, L., Tewarie, P. K., Geurts, J. J., Barkhof, F., & Vrenken, H. (2015). Unraveling the relationship between regional gray matter atrophy and pathology in connected white matter tracts in long-standing multiple sclerosis. Human brain mapping, 36(5), 1796–1807. https://doi.org/10.1002/hbm.22738
    OpenUrlCrossRefPubMed
  57. ↵
    Thompson, A. J., Banwell, B. L., Barkhof, F., Carroll, W. M., Coetzee, T., Comi, G., Correale, J., Fazekas, F., Filippi, M., Freedman, M. S., Fujihara, K., Galetta, S. L., Hartung, H. P., Kappos, L., Lublin, F. D., Marrie, R. A., Miller, A. E., Miller, D. H., Montalban, X., … Cohen, J. A. (2018). Diagnosis of multiple sclerosis: 2017 revisions of the McDonald criteria. The Lancet. Neurology, 17(2), 162–173. https://doi.org/10.1016/S1474-4422(17)30470-2
    OpenUrl
  58. ↵
    Timmers, I., Zhang, H., Bastiani, M., Jansma, B. M., Roebroeck, A., & Rubio-Gozalbo, M. E. (2015). White matter microstructure pathology in classic galactosemia revealed by neurite orientation dispersion and density imaging. Journal of inherited metabolic disease, 38(2), 295–304. https://doi.org/10.1007/s10545-014-9780-x
    OpenUrlCrossRef
  59. ↵
    Tournier, J.-D., Mori, S., & Leemans, A. (2011). Diffusion tensor imaging and beyond. Magnetic resonance in medicine, 65(6), 1532–1556. https://doi.org/10.1002/mrm.22924
    OpenUrlCrossRefPubMedWeb of Science
  60. ↵
    van Waesberghe, J. H., Kamphorst, W., De Groot, C. J., van Walderveen, M. A., Castelijns, J. A., Ravid, R., Lycklama à Nijeholt, G. J., van der Valk, P., Polman, C. H., Thompson, A. J., & Barkhof, F. (1999). Axonal loss in multiple sclerosis lesions: Magnetic resonance imaging insights into substrates of disability. Annals of Neurology, 46(5), 747–754. https://doi.org/10.1002/1531-8249(199911)46:5<747::aid-ana10>3.3.co;2-w
    OpenUrlCrossRefPubMedWeb of Science
  61. ↵
    Vucic, S., Burke, T., Lenton, K., Ramanathan, S., Gomes, L., Yannikas, C., & Kiernan, M. C. (2012). Cortical dysfunction underlies disability in multiple sclerosis. Multiple Sclerosis Journal, 18(4), 425–432. https://doi.org/10.1177/1352458511424308
    OpenUrl
  62. ↵
    Wahl, M., Lauterbach-Soon, B., Hattingen, E., Jung, P., Singer, O., Volz, S., Klein, J. C., Steinmetz, H., & Ziemann, U. (2007). Human Motor Corpus Callosum: Topography, Somatotopy, and Link between Microstructure and Function. Journal of Neuroscience, 27(45), 12132–12138. https://doi.org/10.1523/JNEUROSCI.2320-07.2007
    OpenUrlAbstract/FREE Full Text
  63. ↵
    Wakana, S., Caprihan, A., Panzenboeck, M. M., Fallon, J. H., Perry, M., Gollub, R. L., Hua, K., Zhang, J., Jiang, H., Dubey, P., Blitz, A., van Zijl, P., & Mori, S. (2007). Reproducibility of quantitative tractography methods applied to cerebral white matter. NeuroImage, 36(3), 630–644. https://doi.org/10.1016/j.neuroimage.2007.02.049
    OpenUrlCrossRefPubMedWeb of Science
  64. ↵
    Wang, N., Zhang, J., Cofer, G., Qi, Y., Anderson, R. J., White, L. E., & Johnson, G. A. (2019). Neurite orientation dispersion and density imaging of mouse brain microstructure. Brain Structure and Function, 224(5), 1797–1813. https://doi.org/10.1007/s00429-019-01877-x
    OpenUrl
  65. ↵
    Wilson, S. A., Thickbroom, G. W., & Mastaglia, F. L. (1993). Transcranial magnetic stimulation mapping of the motor cortex in normal subjects: The representation of two intrinsic hand muscles. Journal of the neurological sciences, 118(2), 134–144. https://doi.org/10.1016/0022-510X(93)90102-5
    OpenUrlCrossRefPubMedWeb of Science
  66. ↵
    Winkler, A. M., Ridgway, G. R., Webster, M. A., Smith, S. M., & Nichols, T. E. (2014). Permutation inference for the general linear model. NeuroImage, 92, 381–397. https://doi.org/10.1016/j.neuroimage.2014.01.060
    OpenUrlCrossRefPubMedWeb of Science
  67. ↵
    Zakzanis, K. K., Mraz, R., & Graham, S. J. (2005). An fMRI study of the trail making test. Neuropsychologia, 43(13), 1878–1886. https://doi.org/10.1016/j.neuropsychologia.2005.03.013
    OpenUrlCrossRefPubMedWeb of Science
  68. ↵
    Zhang, H., Hubbard, P. L., Parker, G. J., & Alexander, D. C. (2011). Axon diameter mapping in the presence of orientation dispersion with diffusion MRI. NeuroImage, 56(3), 1301–1315. https://doi.org/10.1016/j.neuroimage.2011.01.084
    OpenUrlCrossRefPubMedWeb of Science
  69. ↵
    Zhang, H., Schneider, T., Wheeler-Kingshott, C. A., & Alexander, D. C. (2012). NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain. NeuroImage, 61(4), 1000–1016. https://doi.org/10.1016/j.neuroimage.2012.03.072
    OpenUrlCrossRefPubMedWeb of Science
  70. ↵
    Zipser, C. M., Premoli, I., Belardinelli, P., Castellanos, N., Rivolta, D., Heidegger, T., MüllerDahlhaus, F., & Ziemann, U. (2018). Cortical excitability and interhemispheric connectivity in early relapsing–remitting multiple sclerosis studied with TMS-EEG. Frontiers in neuroscience, 12, 393. https://doi.org/10.3389/fnins.2018.00393
    OpenUrl
Back to top
Posted October 20, 2020.
Download PDF

Supplementary Material

Data/Code
Email
Linking microstructural integrity and motor cortex excitability in multiple sclerosis
Angela Radetz, Kalina Mladenova, Dumitru Ciolac, Gabriel Gonzalez-Escamilla, Vinzenz Fleischer, Erik Ellwardt, Julia Krämer, Stefan Bittner, Sven G. Meuth, Muthuraman Muthuraman, Sergiu Groppa
medRxiv 2020.10.15.20213090; doi: https://doi.org/10.1101/2020.10.15.20213090
Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Citation Tools

Subject Area