Background The extent and clinical relevance of grey matter (GM) pathology in multiple sclerosis (MS) are increasingly recognised. GM pathology may present as focal lesions, which can be visualised using double inversion recovery (DIR) MRI, or as diffuse pathology, which can manifest as atrophy. It is, however, unclear whether the diffuse atrophy centres on focal lesions. This study aimed to determine if GM lesions and GM atrophy colocalise, and to assess their independent relationship with motor and cognitive deficits in MS.
Methods Eighty people with MS and 30 healthy controls underwent brain volumetric T1-weighted and DIR MRI at 3 T, and had a comprehensive neurological and cognitive assessment. Probability mapping of GM lesions marked on the DIR scans and voxel- based morphometry (assessing GM atrophy) were carried out. The associations of GM lesion load and GM volume with clinical scores were tested.
Results DIR-visible GM lesions were most commonly found in the right cerebellum and most apparent in patients with primary progressive MS. Deep GM structures appeared largely free from lesions, but showed considerable atrophy, particularly in the thalamus, caudate, pallidum and putamen, and this was most apparent in secondary progressive patients with MS. Very little co-localisation of GM atrophy and lesions was seen, and this was generally confined to the cerebellum and postcentral gyrus. In both regions, GM lesions and volume independently correlated with physical disability and cognitive performance.
Conclusions DIR-detectable GM lesions and GM atrophy do not significantly overlap in the brain but, when they do, they independently contribute to clinical disability.
- MULTIPLE SCLEROSIS
- IMAGE ANALYSIS
This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/
Statistics from Altmetric.com
Grey matter (GM) pathology has emerged as a significant and clinically relevant component of multiple sclerosis (MS). Advances in MRI technology have allowed the assessment of GM pathology in MS in vivo. Both GM volume loss (atrophy) and GM lesions occur in MS, and both correlate with neurological and cognitive deficits.1 However, a key question is whether they share a common pathogenic mechanism, in particular whether lesions are the cause of atrophy or are caused by an independent process that contributes to clinical outcomes.
Histopathological studies have identified extensive cortical demyelination in people with MS and in those with progressive disease, GM lesion volume may exceed that of WM lesions.2 However, it has proven difficult to detect in vivo GM lesions using conventional MRI techniques. Only a small percentage of GM lesions are identified on T2-weighted (<10%) and on fluid-attenuated inversion recovery (FLAIR) scans.3 The development of double inversion recovery (DIR) MRI has nearly doubled the detection of GM lesions.4 DIR studies demonstrate that GM lesions are spread throughout the cortex, appear early in the course of the disease and accumulate over time.5 ,6
Deep and cortical GM atrophy is now well-recognised in MS and appears to accelerate as patients with MS enter a progressive phase of the disease.1 GM atrophy is thought to mark irreversible tissue loss and it is likely to reflect a combination of neuronal morphological changes and loss, and glial abnormalities.1 Previous in vivo work provides insight into the mechanisms of GM injury in MS. Two studies have assessed the relationship between GM atrophy and WM injury, and these showed that lesions in connecting WM tracts are associated with deep GM atrophy.7 ,8 Lesions could directly cause localised atrophy, as suggested by pathological studies showing axonal transection and loss in cortical lesions.1 There has been little histopathological work examining the co-localisation of cortical atrophy and demyelination, but in the only study we are aware of, local cortical thickness did not correlate with demyelination.9 However, fixation can affect cortical thickness (and fixation time may be correlated with it9), and so it is possible that this may have obscured an association. As such, it is preferable to look for associations in fresh tissue samples or, better still, in vivo.
This study sought to clarify the spatial overlap between GM atrophy and GM lesions, as well as their independent relationship with cognitive and physical disability, in a large cohort of patients with MS and in different MS subtypes.
We recruited 80 people who fulfil the following inclusion criteria: diagnosis of clinically definite MS;10 absence of a relapse or use of corticosteroids within the preceding 4 weeks; age <65 years; and no other neurological conditions which could have influenced the pattern of GM atrophy. Thirty healthy volunteers with no known neurological disease were also studied. Patients were consecutively recruited from the National Hospital's specialist MS clinics. Status of disease modifying therapy (DMT), antidepressant for either depression or neuropathic pain, and benzodiazepine for spasms or insomnia were recorded. All participants gave written informed consent.
Clinical status was assessed using the Expanded Disability Status Scale (EDSS)11 and the MS functional composite score (MSFC) which includes walking speed on the 25-foot Timed Walk Test (25TWT), 9-hole Peg Test (9HPT) and the Paced Auditory Serial Addition Test (PASAT). Z-scores of these tests were calculated using published means and SDs. Levels of anxiety and depression were also assessed using the Hospital Anxiety and Depression Scale.12
All subjects underwent cognitive testing to assess executive function and memory. Executive function was assessed using the Hayling Sentence Completion Task13 and Stroop Task,14 from which averaged z-scores were calculated on the basis of healthy control performance, and the Symbol Digit Modalities Test (SDMT), from which the age-adjusted z-scores were calculated based on published normative values.15 Memory function was evaluated using a composite score of story recall (immediate and 30 min delay) and figure recall (immediate and 30 min delay) from the Adult Memory and Information Processing Battery,16 and word and face recognition from the Recognition Memory Test;17 z-scores based on the performance of the control sample were also calculated for these tests. When participants were not able to complete the task, their missing data would not be considered for the task's analysis; these were executive function in one patient with primary progressive MS (PPMS), 25TWT in 3 PPMS, 1 relapsing–remitting MS (RRMS) and 9 patients with secondary progressive MS (SPMS), and 9HPT in 1 RRMS and 1 patient with SPMS.
Brain MRI was performed on a Phillips 3 T Achieva TX system (Philips Healthcare, Best, the Netherlands) using a 32-channel receive-only coil. T1-weighted (TR=6.9 ms, TE=3.1 ms, TI=824.5 ms, SENSE=2, voxel size=1×1×1 mm), turbo FLAIR (TR=8000 ms, TE=125 ms, TI=2400 ms, SENSE=1.3, voxel size=1×1×3 mm), and DIR (TR=16 000 ms, TE=9.9 ms, TI=2400/325 ms, SENSE=4.16, 1×1×3 mm) sequences were acquired.
Image registration and lesion mapping
To limit the impact of WM lesions on tissues segmentations, T1-weighted (T1w) hypointense lesions were filled18 and the lesion-filled T1w images were then segmented using SPM8 (http://www.fil.ion.ucl.ac.uk/spm/). Anatomical normalisation to Montreal Neurological Institute (MNI) space was achieved via a custom diffeomorphic anatomical registration through exponentiated Lie algebra (DARTEL) template,19 generated from all subjects’ (N=110) GM tissue segmentations. To make use of MNI atlases, the DARTEL GM template was affine-registered to the MNI standard space, and each subject's T1-weighted GM was first non-linearly registered to the DARTEL template, and subsequently affine registered to MNI using the DARTEL to MNI transformation, moving it into MNI template space. This pipeline reduces the adverse effect of disease-associated brain atrophy on the registration accuracy. Spatially normalised T1-weighted images were averaged over subjects to form a study-specific T1-weighted MNI-space template. When spatially normalising the segmented GM images to MNI space (via the DARTEL and affine transformations), the images were modulated by the Jacobian determinants so that intensities represent the amount of deformation needed to normalise the images. All registrations were reviewed by SvdP to confirm their accuracy. Modulated, normalised, segmented GM images were smoothed with an 8 mm full-width at half-maximum Gaussian kernel.20
GM lesions were marked on the DIR scans using JIM (V.6.0, Xinapse Systems, Northants, UK) by two independent raters (VS and SvdP) according to consensus guidelines.21 Marking of lesions was compared between raters and a consensus was reached on all lesions. Subsequently, one in five scans were analysed by DC, after which a final consensus was reached. Total GM lesion volume was then calculated.
Each participant's DIR scan was affine registered to their T1-weighted scan; their DIR-to-T1 affine transformation, T1-to-DARTEL deformation field and DARTEL-to-MNI affine transformation were combined and used to move the binarised GM lesion mask into MNI space. All registrations were carried out in SPM8. After a nearest-neighbor interpolation, an 8 mm full-width at half-maximum Gaussian Kernel was then used to smooth the normalised lesion masks.
Two patients (one with SPMS, one with PPMS) were excluded from the subsequent imaging analysis due to inadequate registration.
Co-localisation of GM atrophy and lesions
Voxel-based morphometry (VBM) and GM lesion probability mapping (LPM) analyses were carried out in SPM8. Two types of analyses were carried out: voxel-wise comparisons and region of interest (ROI) analyses. For voxel-wise VBM comparison between groups, we used a significance level of 0.05 (family wise error (FWE) corrected). For voxel-wise LPM, we used a significance level of 0.001 (uncorrected). This threshold was used as no lesion clustering was found at 0.05 (FWE corrected). In VBM and LPM comparisons, people with MS were compared to healthy controls to assess areas significantly more affected in people with MS.
To determine co-localisation of volume loss and DIR-visible GM lesions, permutation tests were run using Randomise22 implemented in FSL (FMRIB's software library, http://www.fmrib.ox.ac.uk/fsl). For each of the MS groups, the following regression model was tested per voxel (ie, with LPM as a voxel-wise covariate; sometimes known as biological parametric mapping23):
Total intracranial volume (ICV) was estimated by summing the thresholded GM, white matter (WM) and cerebrospinal fluid volumes using the ‘get_totals’ function in SPM8.
FSL Randomise uses a non-parametric permutation-testing procedure that derives an empirical null distribution without making assumptions about the normality of the data. The tests were performed in an inclusive GM mask (the average of the GM segments over subjects, thresholded at 0.5), using a cluster-forming threshold, and an uncorrected cluster-extent-based significance level of 0.01, with 5000 permutations.
The second analysis employed an ROI approach. ROIs involved in tasks assessing clinical and cognitive functioning in MS24 ,25 were chosen a priori. These regions were the bilateral cerebellum, medial temporal lobe, postcentral gyrus, precentral gyrus, insula, prefrontal cortex and thalamus (figure 1). Masks for these regions were created by using Freesurfer to automatically segment the MNI-space cohort-specific T1-weighted template.26 ,27 For each GM ROI and in each patient, the total GM lesion volume and the GM volume were extracted from the VBM and LPM images, and co-localisation between atrophy and lesion load was assessed using linear regression analyses, correcting for age, sex, and ICV.
Associations of disability with GM lesion load and GM atrophy
To explore the associations of DIR-visible GM lesion load and GM atrophy with clinical measures (EDSS, PASAT, 9HPT, 25TWT, executive function and memory function), voxel-wise and ROI analyses were carried out. The first analysis, based on a voxel-wise factorial design, was conducted in SPM8, including age, sex and ICV, and using VBM or LPM values. In all voxel-wise comparisons, a cluster threshold of five voxels was applied.
Linear regression models used EDSS, PASAT, 9HPT, 25TWT, executive function and memory function as the dependent variables; GM lesion load and GM volume extracted from the ROIs were used as the independent variables. Age, sex and ICV were also added to the model as additional covariates. SPSS (V.21.0. Armonk, New York, USA: IBM Corp) was used to conduct this analysis. Differences between subtypes in scores, volume loss and lesion load were analysed using unpaired t tests (for VBM values) or non-parametric tests (for lesion-loads since they were non-normally distributed). When both GM lesion load and GM atrophy within a given ROI were associated with clinical performance, a linear regression model including lesion load, total GM volume, age and sex was run to determine the independent contribution of each to the outcome of interest. Shapiro-Wilk tests were used to assess normality of residuals of the linear regression analyses and unpaired t tests.
Demographics and clinical performance
Thirty of the MS group had relapsing remitting MS (RRMS), 25 had PPMS, and 25 had SPMS. Demographics are shown in table 1.
All MS subgroups differed significantly in their EDSS scores (all p<0.05), with SPMS having a higher EDSS than patients with PPMS, who in turn had a higher EDSS than patients with RRMS (see table 1). All patient groups had significantly higher levels of depression than controls (all p<0.05), but no significant difference was found between subgroups. Anxiety levels were comparable between patients and controls. People with MS performed worse than healthy controls in all clinical domains. Patients had poorer executive functioning and memory function than healthy controls (all p<0.01). Memory functioning was best in patients with RRMS, followed by patients with PPMS, who in turn performed better than the SPMS subgroup (p<0.05 and p<0.01, respectively). When correcting for depression on clinical functioning, except for the subgroup comparison of patients with SPMS with PPMS on memory function, all group comparisons remained significant on all domains.
People with MS showed significant GM atrophy compared to controls (figure 2A), predominantly in deep GM structures (thalamus, pallidum, putamen and caudate) in addition to a few small regions in the frontal (0.09 cm3), insular (0.06 cm3) and temporal lobes (0.12 cm3). Similarly, in the a priori defined ROIs (ie, cerebellum, medial temporal lobe, postcentral gyrus, precentral gyrus, insula, prefrontal cortex and thalamus, figure 1), mean GM volume was significantly smaller in patients than controls in the thalamus (p<0.001) and insula (p<0.05). Subgroup analyses, which showed that the most extensive atrophy in the deep GM regions and cortical areas was seen in SPMS, are reported in the online supplementary material.
In patients, 1476 DIR visible lesions were found throughout the GM (neocortex: 1276, cerebellum: 154, deep GM: 46), with a mean total lesion volume of 1.20 cm3 (SD=0.92 cm3). Fourteen GM lesions were identified in healthy controls (in 4 participants).
At a threshold of p<0.001, patients had significantly higher GM lesion probability only in the right cerebellar hemisphere when compared with healthy controls (0.07 cm3, tmax=3.37; figure 2B). No deep GM structures showed a significantly increased lesion probability at this threshold or at a lower threshold (p<0.01, uncorrected). However, at this lower threshold, additional cortical areas showed increased lesion probability, notably the right precentral and postcentral gyri, bilateral supplementary motor area, and bilateral temporal lobes (data not shown). People with MS had significantly more lesions in every ROI than healthy volunteers (all p<0.001). Subgroup analyses, which showed that PPMS and SPMS had regions of higher GM lesion probability than healthy controls, are detailed in the online supplementary material.
Co-localisation of GM atrophy and lesions
In the whole MS group, voxel-wise analyses showed that there was little co-localisation of GM atrophy and lesions. Overall the volume of co-localisation was much smaller (less than 1 cm3) when compared with total brain GM volumes (>500 cm3). Areas that showed a significant association between a higher probability of a voxel being lesional and smaller GM volume were scattered throughout the brain, particularly in the cerebellum (figure 2C). ROI analyses found that increased GM lesion load was significantly associated with reduced GM volume within the cerebellum (B=−17.76, p<0.001) and postcentral gyrus (B=−4.77, p=0.037). The greatest extent of co-localisation was found in the PPMS cohort, while no co-localisation was seen in any group in the deep GM regions (see online supplementary material for more details on the subgroup analyses).
Associations of GM atrophy and lesion load with clinical status
Voxel-wise models at p<0.05 (FWE corrected) in people with MS revealed an association between lower executive function scores and decreased left putamen volume (0.10 cm3), which was the only significant voxel-wise association of volume loss with any clinical metric. Furthermore, at p<0.001 (uncorrected), increased lesion probability in the cerebellum in particular was associated with lower performance in executive function, TWT speed and PASAT (table 2).
ROI analyses showed that cerebellar volume loss and increased GM lesion volume were associated with poorer functioning in all cognitive domains tested (table 3). In addition, lesions and volume loss within the postcentral gyrus were associated with cognitive performance and physical disability in almost all domains tested (ie, executive function, memory function, TWT speed, PASAT, 9HPT). Overall, lower functional scores were only weakly linked to GM lesion or atrophy clusters.
The present study assessed the distributions of DIR-visible GM lesions and GM atrophy in patients with MS, aiming to determine if they co-localise and if they both contribute to clinical outcomes. Co-localisation was at best modest, with the majority of atrophy occurring in regions that showed few GM lesions. Of the GM regions found to be consistently atrophic or to contain lesions, only the cerebellum and postcentral gyrus showed associations between both types of pathology, and physical and cognitive function.
Deep GM volume was reduced in patients compared to controls. This is in line with previous findings.1 ,28 More limited cortical atrophy was also seen in regions consistent with past literature, including the right lateral prefrontal cortex,29 the left temporal and prefrontal cortex,30 and the bilateral superior and medial frontal gyrus and orbitofrontal regions.31 In addition, small areas of cortical atrophy were seen in the frontal, insular and temporal lobes, as reported previously.5 ,32 The most extensive atrophy was seen in SPMS and a feature common to all subtypes was thalamic atrophy, which was particularly extensive in the SPMS group; this is in line with previous findings.28 ,33
GM lesions were found throughout the cerebral and cerebellar cortex, particularly in the right cerebellum, but rarely in the deep GM structures. Co-localisation between GM atrophy and GM lesions was observed only when lower statistical thresholds were used, and then too mainly in the cortical regions. It was rarely seen in the deep structures, which showed atrophy. Work by Wegner et al9 has shown pathological changes within cortical lesions and to a lesser degree, in non-lesional cortex, but did not find an association between demyelination and cortical thickness, suggesting that demyelination per se is not directly responsible for cortical atrophy. Our findings are in agreement with this and we, therefore, suggest with caveats (see the discussion of the limitations of this work below) that GM atrophy and DIR-visible GM lesion are often not directly linked. This view is supported by recent findings that cortical MTR abnormalities, which is heavily influenced by demyelination, rarely co-localises with atrophy in MS,20 further suggesting differences in pathogenic mechanisms.
Previous papers have demonstrated that the topographic distribution of GM lesions is rather similar between patient subtypes, and that PPMS and RRMS share more similarities than differences in terms of GM lesion number, volume and topographic distribution.5 Additionally, there are differences in GM atrophy between subtypes of MS which involve not only a selective involvement of brain regions but also an increased extent of atrophy in common regions, such as the thalamus, in the progressive phase of MS.34 We confirmed these findings in this paper and therefore, combined all the phenotypes to investigate the overlap between lesions and atrophy. Based on our observed differences between the subgroups, we performed a sample size calculation to estimate the number of subjects a future study might require to detect (with 80% power at 5% signiﬁcance) a difference between subgroups in the volume of the prefrontal cortex (which was the ROI with the least power to detect differences) of 8 cm3 (as observed in our study), as well as in the proportion of lesional voxels of 0.22, using the means and SDs of concentrations provided by our data and standard methods for comparisons of means.35 We found that the sample size required to detect signiﬁcant GM volume differences between subtypes of MS was 75 per group for GM volume and 39 per group for GM lesion probability. This sample size calculation can be used by future studies that will aim to investigate the differences between groups in lesion probability and GM volume.
With regard to the association between GM lesions and atrophy with clinical disability, one interesting observation from the ROI analysis was that both GM lesion volume and atrophy in the postcentral gyrus (the primary somatosensory cortex) were linked to clinical measures of physical disability and to worse performance in almost all cognitive domains. Additionally, a significant and independent association was found between cerebellar lesion load and atrophy with lower physical and cognitive performance. This fits with the known involvement of the cerebellum in sensorimotor and cognitive function.36 In addition, the postcentral gyrus and cerebellum, which are both functionally and structurally connected,36 were the only regions in which lesion burden and volume loss were linked. Further work focusing on these two regions may shed light on the mechanisms of disability in MS. Using the voxel-wise analysis, we found that deep GM atrophy was generally not linked to functional outcomes, except for an association between executive function and a small region of atrophy in the putamen. This suggests that other processes, including white matter damage and brain network dysfunction, may be involved in cognitive dysfunction.37 ,38
A number of limitations of this work should be taken into consideration. For LPM, some forms of cortical lesions may be easier to detect on DIR scans than others. DIR scans detect about a fifth of cortical lesions (18%) and less than a tenth of deep GM lesions (7%); also, while subpial lesions are the most abundant type seen in postmortem studies, these are rarely observed using DIR.4 Recent pathological work showed deep GM lesions and neurodegeneration in the deep GM that, in combination with cortical and white matter lesions, contributed to the clinical deficits;39 this insensitivity for DIR to detect deep GM lesions may partly explain the limited overlap found between lesions and atrophy in deeper GM structures. However, since the DIR-visible GM lesions represent the ‘peak’ of the real GM lesions, patients with higher DIR-visible lesions will also have higher number of real GM lesions; in this study we have investigated if these patients also have a greater amount of atrophy. Phase Sensitive Inversion Recovery MRI detects two to three times more GM lesions and so may increase the overall sensitivity of LPM analyses, but this sequence may still be less sensitive to subpial lesions than the other GM lesion subtypes.40 However, recent work has shown similar dissociations between the localisation of atrophy and demyelination, as measured by MTR,20 indicating that separations in the localisation and drivers of pathology are likely. A second limitation is that due to the low-frequency of overlapping lesions between patients we used a significance level of 0.001 (uncorrected) for the LPM analysis. While lowering the threshold increases the chance of a type I error, only one per 1000 voxels would show a false positive. In addition, the applied cluster threshold further reduces the number of false positive findings. Few studies have examined the effect of psychotropic drugs on cognition in MS. Oken et al41 found no significant effect of SSRIs and benzodiazepines on the PASAT, stroop, verbal memory or executive function test performance. These did, however, show an effect on a reaction time test. While use of psychotropic drugs may have subtly influenced our findings, the concomitant alleviation of depression, which is itself linked to cognitive impairment, is likely to have balanced out these effects. Not all patients were able to finish all clinical tasks due to disability, which led to the exclusion of their missing data from this task's analysis. It is, therefore, likely that the effects found were an underestimation of the real effects of MS as the excluded disabled patients were likely to also have the poorest cognitive function.
The present study assessed the distributions of DIR-visible GM lesions and GM atrophy in patients with different subtypes of MS and their relationship with clinical outcomes. We found that GM lesions (as seen using DIR) and GM atrophy do not usually co-localise, indicating that they are not directly spatially coupled. We also found that both GM lesions and atrophy separately contributed to disability, suggesting that the substrates of disability in MS are both pathologically and spatially heterogeneous.
The authors would like to thank all the participants in this study, as well as Daniel Altmann (London School of Hygiene and Tropical Medicine)and Carmen Tur (UCL Institute of Neurology) for help with some of the statistical analyses.
Review history and Supplementary material
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
- Data supplement 1 - Online supplement
Contributors NM, DTC, DHM, AJT and OC were involved in the conception and design of the study. NM and VS recruited and obtained consent from participants. SHPvdP, NM, VS and GRR conducted the MRI processing and statistical analyses. All authors contributed to preparation of the manuscript.
Funding This work was supported by the MS Society of Great Britain and Northern Ireland (grant numbers 892/08 and 917/09). This work was undertaken at University College London Hospitals/ University College London who received a proportion of funding from the Department of Heatlh%s National Institute for Health Research Biomedical Research Centre%s funding scheme. SHPvdP received research support from Stichting MS Research and Prinses Beatrix Fonds. NM is supported by a Wellcome Trust ISSF fellowship and a Welsh Government NISCHR Health Fellowship. GR is supported by the Medical Research Council (grant number MR/J014257/1). The Wellcome Trust Centre for Neuroimaging is supported by core funding from the Wellcome Trust (grant number 091593/Z/10/Z).
Competing interests None declared.
Patient consent Obtained.
Ethics approval UCL research ethics committee.
Provenance and peer review Not commissioned; externally peer reviewed.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.