Segmentation of brain magnetic resonance images for measurement of gray matter atrophy in multiple sclerosis patients

Natalizumab reduces loss of gray matter and thalamic volume in patients with relapsing-remitting multiple sclerosis: A post hoc analysis from the randomized, placebo-controlled AFFIRM trial

BACKGROUND

Loss of brain gray matter fractional volume predicts multiple sclerosis (MS) progression and is associated with worsening physical and cognitive symptoms. Within deep gray matter, thalamic damage is evident in early stages of MS and correlates with physical and cognitive impairment. Natalizumab is a highly effective treatment that reduces disease progression and the number of inflammatory lesions in patients with relapsing-remitting MS (RRMS).

OBJECTIVE

To evaluate the effect of natalizumab on gray matter and thalamic atrophy.

METHODS

A combination of deep learning-based image segmentation and data augmentation was applied to MRI data from the AFFIRM trial.

RESULTS

This post hoc analysis identified a reduction of 64.3% (p = 0.0044) and 64.3% (p = 0.0030) in mean percentage gray matter volume loss from baseline at treatment years 1 and 2, respectively, in patients treated with natalizumab versus placebo. The reduction in thalamic fraction volume loss from baseline with natalizumab versus placebo was 57.0% at year 2 (p < 0.0001) and 41.2% at year 1 (p = 0.0147). Similar findings resulted from analyses of absolute gray matter and thalamic fraction volume loss.

CONCLUSION

These analyses represent the first placebo-controlled evidence supporting a role for natalizumab treatment in mitigating gray matter and thalamic fraction atrophy among patients with RRMS.

CLINICALTRIALS.GOV IDENTIFIER

NCT00027300URL: https://clinicaltrials.gov/ct2/show/NCT00027300.

Grey-matter sodium concentration as an individual marker of multiple sclerosis severity

OBJECTIVE

Quantification of brain injury in patients with variable disability despite similar disease duration may be relevant to identify the mechanisms underlying disability in multiple sclerosis (MS). We aimed to compare grey-matter sodium abnormalities (GMSAs), a parameter reflecting neuronal and astrocyte dysfunction, in MS patients with benign multiple sclerosis (BMS) and non-benign multiple sclerosis (NBMS).

METHODS

We identified never-treated BMS patients in our local MS database of 1352 patients. A group with NBMS was identified with same disease duration. All participants underwent ²³Na magnetic resonance imaging (MRI). The existence of GMSA was detected by statistical analysis.

RESULTS

In total, 102 individuals were included (21 BMS, 25 NBMS and 56 controls). GMSA was detected in 10 BMS and 19 NBMS (11/16 relapsing-remitting multiple sclerosis (RRMS) and 8/9 secondary progressive multiple sclerosis (SPMS) patients) (p = 0.05). On logistic regression including the presence or absence of GMSA, thalamic volume, cortical grey-matter volume and T2-weighted lesion load, thalamic volume was independently associated with BMS status (odds ratio (OR) = 0.64 for each unit). Nonetheless, the absence of GMSA was independently associated when excluding patients with significant cognitive alteration (n = 7) from the BMS group (OR = 4.6).

CONCLUSION

Detection of GMSA in individuals and thalamic volume are promising to differentiate BMS from NBMS as compared with cortical or whole grey-matter atrophy and T2-weighted lesions.

MAGNIMS recommendations for harmonization of MRI data in MS multicenter studies

•

Sharing data from cooperative studies is essential to develop new biomarkers in MS.

•

Differences in MRI acquisition, analysis, storage represent a substantial constraint.

•

We review the state of the art and developments in the harmonization of MRI.

•

We provide recommendations to harmonize large MRI datasets in the MS field.

There is an increasing need of sharing harmonized data from large, cooperative studies as this is essential to develop new diagnostic and prognostic biomarkers. In the field of multiple sclerosis (MS), the issue has become of paramount importance due to the need to translate into the clinical setting some of the most recent MRI achievements. However, differences in MRI acquisition parameters, image analysis and data storage across sites, with their potential bias, represent a substantial constraint. This review focuses on the state of the art, recent technical advances, and desirable future developments of the harmonization of acquisition, analysis and storage of large-scale multicentre MRI data of MS cohorts. Huge efforts are currently being made to achieve all the requirements needed to provide harmonized MRI datasets in the MS field, as proper management of large imaging datasets is one of our greatest opportunities and challenges in the coming years. Recommendations based on these achievements will be provided here. Despite the advances that have been made, the complexity of these tasks requires further research by specialized academical centres, with dedicated technical and human resources. Such collective efforts involving different professional figures are of crucial importance to offer to MS patients a personalised management while minimizing consumption of resources.

Effects of Ibudilast on MRI Measures in the Phase 2 SPRINT-MS Study

OBJECTIVE

To determine whether ibudilast has an effect on brain volume and new lesions in progressive forms of multiple sclerosis (MS).

METHODS

A randomized, placebo-controlled, blinded study evaluated ibudilast at a dose of up to 100 mg over 96 weeks in primary and secondary progressive MS. In this secondary analysis of a previously reported trial, secondary and tertiary endpoints included gray matter atrophy, new or enlarging T2 lesions as measured every 24 weeks, and new T1 hypointensities at 96 weeks. Whole brain atrophy measured by structural image evaluation, using normalization, of atrophy (SIENA) was a sensitivity analysis.

RESULTS

A total of 129 participants were assigned to ibudilast and 126 to placebo. New or enlarging T2 lesions were observed in 37.2% on ibudilast and 29.0% on placebo (p = 0.82). New T1 hypointense lesions at 96 weeks were observed in 33.3% on ibudilast and 23.5% on placebo (p = 0.11). Gray matter atrophy was reduced by 35% for those on ibudilast vs placebo (p = 0.038). Progression of whole brain atrophy by SIENA was slowed by 20% in the ibudilast group compared with placebo (p = 0.08).

CONCLUSION

Ibudilast treatment was associated with a reduction in gray matter atrophy. Ibudilast treatment was not associated with a reduction in new or enlarging T2 lesions or new T1 lesions. An effect on brain volume contributes to prior data that ibudilast appears to affect markers associated with neurodegenerative processes, but not inflammatory processes.

CLASSIFICATION OF EVIDENCE

This study provides Class II evidence that for people with MS, ibudilast does not significantly reduce new or enlarging T2 lesions or new T1 lesions.

A contrast-adaptive method for simultaneous whole-brain and lesion segmentation in multiple sclerosis

Here we present a method for the simultaneous segmentation of white matter lesions and normal-appearing neuroanatomical structures from multi-contrast brain MRI scans of multiple sclerosis patients. The method integrates a novel model for white matter lesions into a previously validated generative model for whole-brain segmentation. By using separate models for the shape of anatomical structures and their appearance in MRI, the algorithm can adapt to data acquired with different scanners and imaging protocols without retraining. We validate the method using four disparate datasets, showing robust performance in white matter lesion segmentation while simultaneously segmenting dozens of other brain structures. We further demonstrate that the contrast-adaptive method can also be safely applied to MRI scans of healthy controls, and replicate previously documented atrophy patterns in deep gray matter structures in MS. The algorithm is publicly available as part of the open-source neuroimaging package FreeSurfer.

Reduced accuracy of MRI deep grey matter segmentation in multiple sclerosis: an evaluation of four automated methods against manual reference segmentations in a multi-center cohort

Background

Deep grey matter (DGM) atrophy in multiple sclerosis (MS) and its relation to cognitive and clinical decline requires accurate measurements. MS pathology may deteriorate the performance of automated segmentation methods. Accuracy of DGM segmentation methods is compared between MS and controls, and the relation of performance with lesions and atrophy is studied.

Methods

On images of 21 MS subjects and 11 controls, three raters manually outlined caudate nucleus, putamen and thalamus; outlines were combined by majority voting. FSL-FIRST, FreeSurfer, Geodesic Information Flow and volBrain were evaluated. Performance was evaluated volumetrically (intra-class correlation coefficient (ICC)) and spatially (Dice similarity coefficient (DSC)). Spearman's correlations of DSC with global and local lesion volume, structure of interest volume (ROIV), and normalized brain volume (NBV) were assessed.

Results

ICC with manual volumes was mostly good and spatial agreement was high. MS exhibited significantly lower DSC than controls for thalamus and putamen. For some combinations of structure and method, DSC correlated negatively with lesion volume or positively with NBV or ROIV. Lesion-filling did not substantially change segmentations.

Conclusions

Automated methods have impaired performance in patients. Performance generally deteriorated with higher lesion volume and lower NBV and ROIV, suggesting that these may contribute to the impaired performance.

Electronic supplementary material

The online version of this article (10.1007/s00415-020-10023-1) contains supplementary material, which is available to authorized users.

Differentiating Parkinson's disease motor subtypes using automated volume-based morphometry incorporating white matter and deep gray nuclear lesion load

Background

Periventricular leukoaraiosis may be an important pathological change in postural instability gait disorder (PIGD), a motor subtype of Parkinson's disease (PD). Clinical diagnosis of PIGD may be challenging for the general neurologist.

Purpose

To evaluate 1) the utility of a fully automated volume‐based morphometry (Vol‐BM) in characterizing imaging diagnostic markers in PD and PIGD, including, 2) novel deep gray nuclear lesion load (GMab), and 3) discriminatory performance of a Vol‐BM model construct in classifying the PIGD subtype.

Study Type

Prospective.

Subjects

In all, 23 PIGD, 21 PD, and 20 age‐matched healthy controls (HC) underwent MRI brain scans and clinical assessments.

Field Strength/Sequence

3.0T, sagittal 3D‐magnetization‐prepared rapid gradient echo (MPRAGE), and fluid‐attenuated inversion recovery imaging (FLAIR) sequences.

Assessment

Clinical assessment was conducted by a movement disorder neurologist. The MR brain images were then segmented using an automated multimodal Vol‐BM algorithm (MorphoBox) and reviewed by two authors independently.

Statistical Testing

Brain segmentation and clinical parameter differences and dependence were assessed using analysis of variance (ANOVA) and regression analysis, respectively. Logistic regression was performed to differentiate PIGD from PD, and discriminative reliability was evaluated using receiver operating characteristic (ROC) analysis.

Results

Significantly higher white matter lesion load (WMab) (P < 0.01), caudate GMab (P < 0.05), and lateral and third ventricular (P < 0.05) volumetry were found in PIGD, compared with PD and HC. WMab, caudate and putamen GMab, and caudate, lateral, and third ventricular volumetry showed significant coefficients (P < 0.005) in linear regressions with balance and gait assessments in both patient groups. A model incorporating WMab, caudate GMab, and caudate GM discriminated PIGD from PD and HC with a sensitivity = 0.83 and specificity = 0.76 (AUC = 0.84).

Data Conclusion

Fast, unbiased quantification of microstructural brain changes in PD and PIGD is feasible using automated Vol‐BM. Composite lesion load in the white matter and caudate, and caudate volumetry discriminated PIGD from PD and HC, and showed potential in classification of these disorders using supervised machine learning.

Level of Evidence: 1

Technical Efficacy: Stage 1

J. Magn. Reson. Imaging 2020;51:748–756.

Measurement of Whole-Brain and Gray Matter Atrophy in Multiple Sclerosis: Assessment with MR Imaging

Purpose To compare available methods for whole-brain and gray matter (GM) atrophy estimation in multiple sclerosis (MS) in terms of repeatability (same magnetic resonance [MR] imaging unit) and reproducibility (different system/field strength) for their potential clinical applications. Materials and Methods The softwares ANTs-v1.9, CIVET-v2.1, FSL-SIENAX/SIENA-5.0.1, Icometrix-MSmetrix-1.7, and SPM-v12 were compared. This retrospective study, performed between March 2015 and March 2017, collected data from (a) eight simulated MR images and longitudinal data (2 weeks) from 10 healthy control subjects to assess the cross-sectional and longitudinal accuracy of atrophy measures, (b) test-retest MR images in 29 patients with MS acquired within the same day at different imaging unit field strengths/manufacturers to evaluate precision, and (c) longitudinal data (1 year) in 24 patients with MS for the agreement between methods. Tissue segmentation, image registration, and white matter (WM) lesion filling were also evaluated. Multiple paired t tests were used for comparisons. Results High values of accuracy (0.87-0.97) for whole-brain and GM volumes were found, with the lowest values for MSmetrix. ANTs showed the lowest mean error (0.02%) for whole-brain atrophy in healthy control subjects, with a coefficient of variation of 0.5%. SPM showed the smallest mean error (0.07%) and coefficient of variation (0.08%) for GM atrophy. Globally, good repeatability (P > .05) but poor reproducibility (P < .05) were found for all methods. WM lesion filling technique mainly affected ANTs, MSmetrix, and SPM results (P < .05). Conclusion From this comparison, it would be possible to select a software for atrophy measurement, depending on the requirements of the application (research center, clinical trial) and its goal (accuracy and repeatability or reproducibility). An improved reproducibility is required for clinical application. ^© RSNA, 2018 Online supplemental material is available for this article.

Urgent challenges in quantification and interpretation of brain grey matter atrophy in individual MS patients using MRI

Atrophy of the brain grey matter (GM) is an accepted and important feature of multiple sclerosis (MS). However, its accurate measurement is hampered by various technical, pathological and physiological factors. As a consequence, it is challenging to investigate the role of GM atrophy in the disease process as well as the effect of treatments that aim to reduce neurodegeneration. In this paper we discuss the most important challenges currently hampering the measurement and interpretation of GM atrophy in MS. The focus is on measurements that are obtained in individual patients rather than on group analysis methods, because of their importance in clinical trials and ultimately in clinical care. We discuss the sources and possible solutions of the current challenges, and provide recommendations to achieve reliable measurement and interpretation of brain GM atrophy in MS.

SIENA-XL for improving the assessment of gray and white matter volume changes on brain MRI

In this article, SIENA-XL, a new segmentation-based longitudinal pipeline is introduced, for: (i) increasing the precision of longitudinal volume change estimation for white (WM) and gray (GM) matter separately, compared with cross-sectional segmentation methods such as SIENAX; and (ii) avoiding potential biases in registration-based methods when Jacobians are used, with a smoothing extent larger than spatial scale between tissue-interfaces, which is where atrophy usually occurs. SIENA-XL implements a new brain extraction procedure and a multi-time-point intensity equalization step before performing the final segmentation that also includes separate segmentation of deep GM structures by using FMRIB's Integrated Registration and Segmentation Tool. The detection of GM and WM volume changes with SIENA-XL was evaluated using different healthy control (HC) and multiple sclerosis (MS) MRI datasets and compared with the traditional SIENAX and two Jacobian-based approaches, SPM12 and SIENAX-JI (a version of SIENAX including Jacobian integration - JI). In scan-rescan data from HCs, SIENA-XL showed: (i) a significant decrease in error, of 50-70% when compared with SIENAX; (ii) no significant differences in error when compared with SIENAX-JI and SPM12 in a scan-rescan HC dataset that included repositioning. When tested in a HC dataset with scan-rescan both at baseline and after 1 year of follow-up, SIENA-XL showed: (i) significantly higher precision (P < 0.01) than SIENAX; (ii) no significant differences to SIENAX-JI and SPM12. Finally, in a dataset of 79 MS patients with a 2 years follow-up, SIENA-XL showed a substantial reduction of sample size, by comparison with SIENAX, SIENAX-JI, and SPM12, for detecting treatment effects of 25, 30, and 50%. Hum Brain Mapp 39:1063-1077, 2018. © 2017 Wiley Periodicals, Inc.

Evaluating the effect of multiple sclerosis lesions on automatic brain structure segmentation

In recent years, many automatic brain structure segmentation methods have been proposed. However, these methods are commonly tested with non-lesioned brains and the effect of lesions on their performance has not been evaluated. Here, we analyze the effect of multiple sclerosis (MS) lesions on three well-known automatic brain structure segmentation methods, namely, FreeSurfer, FIRST and multi-atlas fused by majority voting, which use learning-based, deformable and atlas-based strategies, respectively. To perform a quantitative analysis, 100 synthetic images of MS patients with a total of 2174 lesions are simulated on two public databases with available brain structure ground truth information (IBSR18 and MICCAI'12). The Dice similarity coefficient (DSC) differences and the volume differences between the healthy and the simulated images are calculated for the subcortical structures and the brainstem. We observe that the three strategies are affected when lesions are present. However, the effects of the lesions do not follow the same pattern; the lesions either make the segmentation method underperform or surprisingly augment the segmentation accuracy. The obtained results show that FreeSurfer is the method most affected by the presence of lesions, with DSC differences (generated - healthy) ranging from - 0.11 ± 0.54 to 9.65 ± 9.87, whereas FIRST tends to be the most robust method when lesions are present (- 2.40 ± 5.54 to 0.44 ± 0.94). Lesion location is not important for global strategies such as FreeSurfer or majority voting, where structure segmentation is affected wherever the lesions exist. On the other hand, FIRST is more affected when the lesions are overlaid or close to the structure of analysis. The most affected structure by the presence of lesions is the nucleus accumbens (from - 1.12 ± 2.53 to 1.32 ± 4.00 for the left hemisphere and from - 2.40 ± 5.54 to 9.65 ± 9.87 for the right hemisphere), whereas the structures that show less variation include the thalamus (from 0.03 ± 0.35 to 0.74 ± 0.89 and from - 0.48 ± 1.08 to - 0.04 ± 0.22) and the brainstem (from - 0.20 ± 0.38 to 1.03 ± 1.31). The three segmentation approaches are affected by the presence of MS lesions, which demonstrates that there exists a problem in the automatic segmentation methods of the deep gray matter (DGM) structures that has to be taken into account when using them as a tool to measure the disease progression.

Application of diffusional kurtosis imaging to detect occult brain damage in multiple sclerosis and neuromyelitis optica

Multiple sclerosis (MS) and neuromyelitis optica (NMO) are two common types of inflammatory demyelinating disease of the central nervous system. Early distinction of NMO from MS is crucial but quite challenging. In this study, 13 NMO spectrum disorder patients (Expanded Disability Status Scale (EDSS) of 3.0 ± 1.7, ranging from 2 to 6.5; disease duration of 5.3 ± 4.7 years), 17 relapsing-remitting MS patients (EDSS of 2.6 ± 1.4, ranging from 1 to 5.5; disease duration of 7.9 ± 7.8 years) and 18 healthy volunteers were recruited. Diffusional kurtosis imaging was employed to discriminate NMO and MS patients at the early or stable stage from each other, and from healthy volunteers. The presence of alterations in diffusion and diffusional kurtosis metrics in normal-appearing white matter (NAWM) and diffusely increased mean diffusivity (MD) in the cortical normal-appearing gray matter (NAGM) favors the diagnosis of MS rather than NMO. Meanwhile, normal diffusivities and kurtosis metrics in all NAWM as well as increases in MD in the frontal and temporal NAGM suggest NMO. Our results suggest that diffusion and diffusional kurtosis metrics may well aid in discriminating the two diseases.

Non-Local Means Inpainting of MS Lesions in Longitudinal Image Processing

In medical imaging, multiple sclerosis (MS) lesions can lead to confounding effects in automatic morphometric processing tools such as registration, segmentation and cortical extraction, and subsequently alter individual longitudinal measurements. Multiple magnetic resonance imaging (MRI) inpainting techniques have been proposed to decrease the impact of MS lesions in medical image processing, however, most of these methods make the assumption that lesions only affect white matter. Here, we propose a method to fill lesion regions using the patch-based non-local mean (NLM) strategy. The method consists of a hierarchical concentric filling strategy after identification of the lesion region. The lesion is filled iteratively, based on the surrounding tissue intensity, using an onion peel strategy. This concentric technique presents the advantage of preserving the local information and therefore the continuity of the anatomy and does not require identification of any a priori normal brain tissues. The method is first evaluated on 20 healthy subjects with simulated artificial MS lesions where we assessed our technique by measuring the peak signal-to-noise ratio (PSNR) of the images with inpainted lesion and the original healthy images. Second, in order to assess the impact of lesion filling on longitudinal image analyses, we performed a power analysis with sample size estimation to evaluate brain atrophy and ventricular growth in patients with MS. The method was compared to two different publicly available methods (FSL lesion fill and Lesion LEAP) and a more classic method, which fills the region with intensities similar to that of the surrounding healthy white matter tissue or mask the lesions. The proposed method was shown to exceed the other methods in reproducing the fidelity of healthy subject images where the lesions were inpainted. The method also improved the power to detect brain atrophy or ventricular growth by decreasing the sample size by 25% in the presence of MS lesions.

Quantifying brain tissue volume in multiple sclerosis with automated lesion segmentation and filling

Lesion filling has been successfully applied to reduce the effect of hypo-intense T1-w Multiple Sclerosis (MS) lesions on automatic brain tissue segmentation. However, a study of fully automated pipelines incorporating lesion segmentation and lesion filling on tissue volume analysis has not yet been performed. Here, we analyzed the % of error introduced by automating the lesion segmentation and filling processes in the tissue segmentation of 70 clinically isolated syndrome patient images. First of all, images were processed using the LST and SLS toolkits with different pipeline combinations that differed in either automated or manual lesion segmentation, and lesion filling or masking out lesions. Then, images processed following each of the pipelines were segmented into gray matter (GM) and white matter (WM) using SPM8, and compared with the same images where expert lesion annotations were filled before segmentation. Our results showed that fully automated lesion segmentation and filling pipelines reduced significantly the % of error in GM and WM volume on images of MS patients, and performed similarly to the images where expert lesion annotations were masked before segmentation. In all the pipelines, the amount of misclassified lesion voxels was the main cause in the observed error in GM and WM volume. However, the % of error was significantly lower when automatically estimated lesions were filled and not masked before segmentation. These results are relevant and suggest that LST and SLS toolboxes allow the performance of accurate brain tissue volume measurements without any kind of manual intervention, which can be convenient not only in terms of time and economic costs, but also to avoid the inherent intra/inter variability between manual annotations.

•

SLS and LST pipelines incorporate fully automated lesion segmentation and filling.

•

Both pipelines reduced the % of error in GM and WM volume of MS patient images.

•

Performance was similar to images with lesion annotations masked before segmentation.

•

The amount of misclassified lesion voxels was the main factor in the % of error.

•

SLS/LST can reduce time, economic costs, and the variability between manual annotations.

Effect of intramuscular interferon beta-1a on gray matter atrophy in relapsing-remitting multiple sclerosis: A retrospective analysis

BACKGROUND

Changes in gray matter (GM) volume may be a useful measure of tissue loss in multiple sclerosis (MS).

OBJECTIVES

To investigate the rate, patterns, and disability correlates of GM volume change in an MS treatment clinical trial.

METHODS

Patients (n=140) with relapsing-remitting MS were randomized to intramuscular (IM) interferon (IFN) beta-1a or placebo. Treatment effects on GM fraction (GMF) and white matter (WM) fraction (WMF) changes, differences in rates of GMF and WMF change in year one and two on treatment, and differences in atrophy rates by disease progression status were assessed retrospectively.

RESULTS

Significantly less GM atrophy (during year two), but not WM atrophy (at any point), was observed with IM IFN beta-1a compared with placebo. Pseudoatrophy effects were more apparent in WM than in GM; in year one, greater WM volume loss was observed with IM IFN beta-1a than with placebo, whereas GM volume loss was similar between groups. Risk of sustained disability progression was significantly associated with GM, but not WM, atrophy.

CONCLUSIONS

These results suggest that GMF change is more meaningful than WMF as a marker of tissue loss and may be useful to augment whole brain atrophy measurements in MS clinical trials.

Genetic associations with brain cortical thickness in multiple sclerosis

Multiple sclerosis (MS) is characterized by temporal and spatial dissemination of demyelinating lesions in the central nervous system. Associated neurodegenerative changes contributing to disability have been recognized even at early disease stages. Recent studies show the importance of gray matter damage for the accrual of clinical disability rather than white matter where demyelination is easily visualized by magnetic resonance imaging (MRI). The susceptibility to MS is influenced by genetic risk, but genetic factors associated with the disability are not known. We used MRI data to determine cortical thickness in 557 MS cases and 75 controls and in another cohort of 219 cases. We identified nine areas showing different thickness between cases and controls (regions of interest, ROI) (eight of them were negatively correlated with Kurtzke's expanded disability status scale, EDSS) and conducted genome-wide association studies (GWAS) in 464 and 211 cases available from the two data sets. No marker exceeded genome-wide significance in the discovery cohort. We next combined nominal statistical evidence of association with physical evidence of interaction from a curated human protein interaction network, and searched for subnetworks enriched with nominally associated genes and for commonalities between the two data sets. This network-based pathway analysis of GWAS detected gene sets involved in glutamate signaling, neural development and an adjustment of intracellular calcium concentration. We report here for the first time gene sets associated with cortical thinning of MS. These genes are potentially correlated with disability of MS.

Quantification of global cerebral atrophy in multiple sclerosis from 3T MRI using SPM: the role of misclassification errors

Purpose

We tested the validity of a freely available segmentation pipeline to measure compartmental brain volumes from 3T MRI in patients with multiple sclerosis (MS). Our primary focus was methodological to explore the effect of segmentation corrections on the clinical relevance of the output metrics.

Methods

Three-dimensional T1-weighted images were acquired to compare 61 MS patients to 30 age- and gender-matched normal controls (NC). We also tested the within patient MRI relationship to disability (eg, expanded disability status scale [EDSS] score) and cognition. Statistical parametric mapping v. 8 (SPM8)-derived gray matter (GMF), white matter (WMF), and total brain parenchyma fractions (BPF) were derived before and after correcting errors from T1 hypointense MS lesions and/or ineffective deep GM contouring.

Results

MS patients had lower GMF and BPF as compared to NC (P<.05). Cognitively impaired patients had lower BPF than cognitively preserved patients (P<.05). BPF was related to EDSS; BPF and GMF were related to disease duration (all P<.05). Errors caused bias in GMFs and WMFs but had no discernable influence on BPFs or any MRI-clinical associations.

Conclusions

We report the validity of a segmentation pipeline for the detection of MS-related brain atrophy with 3T MRI. Longitudinal studies are warranted to extend these results.

Subcortical brain segmentation of two dimensional T1-weighted data sets with FMRIB's Integrated Registration and Segmentation Tool (FIRST)

Brain atrophy has been identified as an important contributing factor to the development of disability in multiple sclerosis (MS). In this respect, more and more interest is focussing on the role of deep grey matter (DGM) areas. Novel data analysis pipelines are available for the automatic segmentation of DGM using three-dimensional (3D) MRI data. However, in clinical trials, often no such high-resolution data are acquired and hence no conclusions regarding the impact of new treatments on DGM atrophy were possible so far. In this work, we used FMRIB's Integrated Registration and Segmentation Tool (FIRST) to evaluate the possibility of segmenting DGM structures using standard two-dimensional (2D) T1-weighted MRI. In a cohort of 70 MS patients, both 2D and 3D T1-weighted data were acquired. The thalamus, putamen, pallidum, nucleus accumbens, and caudate nucleus were bilaterally segmented using FIRST. Volumes were calculated for each structure and for the sum of basal ganglia (BG) as well as for the total DGM. The accuracy and reliability of the 2D data segmentation were compared with the respective results of 3D segmentations using volume difference, volume overlap and intra-class correlation coefficients (ICCs). The mean differences for the individual substructures were between 1.3% (putamen) and −25.2% (nucleus accumbens). The respective values for the BG were −2.7% and for DGM 1.3%. Mean volume overlap was between 89.1% (thalamus) and 61.5% (nucleus accumbens); BG: 84.1%; DGM: 86.3%. Regarding ICC, all structures showed good agreement with the exception of the nucleus accumbens. The results of the segmentation were additionally validated through expert manual delineation of the caudate nucleus and putamen in a subset of the 3D data. In conclusion, we demonstrate that subcortical segmentation of 2D data are feasible using FIRST. The larger subcortical GM structures can be segmented with high consistency. This forms the basis for the application of FIRST in large 2D MRI data sets of clinical trials in order to determine the impact of therapeutic interventions on DGM atrophy in MS.

•

Segmentation of deep grey matter (DGM) using 2-dimensional MRI data is challenging.

•

We performed a systematic comparison between 2 and 3D data segmentation using FIRST.

•

The DGM volumes of 70 multiple sclerosis patients were statistically compared.

•

We found a good agreement of the segmentation results for the large DGM structures.

•

Our results indicate that FIRST is suitable for the segmentation of 2D data sets acquired in multicentre studies.

White matter lesion filling improves the accuracy of cortical thickness measurements in multiple sclerosis patients: a longitudinal study

Background

Previous studies have demonstrated that white matter (WM) lesions bias automated brain tissue classifications and cerebral volume measurements. However, filling WM lesions using the intensity of neighbouring normal-appearing WM has been shown to increase the accuracy of automated volume measurements in the brain. In the present study, we investigate the influence of WM lesions on cortical thickness (CTh) measures and assessed the impact of lesion filling on both cross-sectional/longitudinal and global/regional measurements of CTh in multiple sclerosis (MS) patients.

Methods

Fifty MS patients were studied at baseline as well as after three and six years of follow-up. CTh was estimated using a fully automated pipeline (CIVET) on T1-weighted magnetic resonance images data acquired at 1.5 Tesla without (original) and with WM lesion filling (filled). WM lesions were semi-automatically segmented and then filled with the mean intensity of the neighbouring voxels. For both original and filled T1 images we investigated and compared the main CIVET’s steps: tissue classification, surfaces generation and CTh measurement.

Results

On the original T1 images, the majority of WM lesion volume (72%) was wrongly classified as gray matter (GM). After lesion filling the accuracy of WM lesions classification improved significantly (p < 0.001, 94% of WM lesion volume correctly classified) as well as the WM surface generation (p < 0.0001). The mean CTh computed on the original T1 images, overall time points, was significantly thinner (p < 0.001) compared the CTh estimated on the filled T1 images. The vertex-wise longitudinal analysis performed on the filled T1 images showed an increased number of vertices in the fronto-temporal region with a significantly decrease of CTh over time compared the analysis performed on the original images.

Conclusion

These results indicate that WM lesions bias the CTh estimation both cross-sectionally as well as longitudinally. The lesion filling approach significantly improved the accuracy of the regional CTh estimation and has an impact also on the global estimation of CTh.

Microstructural changes in the striatum and their impact on motor and neuropsychological performance in patients with multiple sclerosis

Grey matter (GM) damage is a clinically relevant feature of multiple sclerosis (MS) that has been previously assessed with diffusion tensor imaging (DTI). Fractional anisotropy (FA) of the basal ganglia and thalamus might be increased in MS patients, and correlates with disability scores. Despite the established role of the striatum and thalamus in motor control, mood and cognition, the impact of DTI changes within these structures on motor and neuropsychological performance has not yet been specifically addressed in MS. We investigated DTI metrics of deep GM nuclei and their potential association with mobility and neuropsychological function. DTI metrics from 3T MRI were assessed in the caudate, putamen, and thalamus of 30 MS patients and 10 controls. Sixteen of the patients underwent neuropsychological testing. FA of the caudate and putamen was higher in MS patients compared to controls. Caudate FA correlated with Expanded Disability Status Scale score, Ambulation Index, and severity of depressive symptomatology. Putamen and thalamus FA correlated with deficits in memory tests. In contrast, cerebral white matter (WM) lesion burden showed no significant correlation with any of the disability, mobility and psychometric parameters. Our findings support evidence of FA changes in the basal ganglia in MS patients, as well as deep GM involvement in disabling features of MS, including mobility and cognitive impairment. Deep GM FA appears to be a more sensitive correlate of disability than WM lesion burden.

Accurate GM atrophy quantification in MS using lesion-filling with co-registered 2D lesion masks

Background

In multiple sclerosis (MS), brain atrophy quantification is affected by white matter lesions. LEAP and FSL-lesion_filling, replace lesion voxels with white matter intensities; however, they require precise lesion identification on 3DT1-images.

Aim

To determine whether 2DT2 lesion masks co-registered to 3DT1 images, yield grey and white matter volumes comparable to precise lesion masks.

Methods

2DT2 lesion masks were linearly co-registered to 20 3DT1-images of MS patients, with nearest-neighbor (NNI), and tri-linear interpolation. As gold-standard, lesion masks were manually outlined on 3DT1-images. LEAP and FSL-lesion_filling were applied with each lesion mask. Grey (GM) and white matter (WM) volumes were quantified with FSL-FAST, and deep gray matter (DGM) volumes using FSL-FIRST. Volumes were compared between lesion mask types using paired Wilcoxon tests.

Results

Lesion-filling with gold-standard lesion masks compared to native images reduced GM overestimation by 1.93 mL (p < .001) for LEAP, and 1.21 mL (p = .002) for FSL-lesion_filling. Similar effects were achieved with NNI lesion masks from 2DT2. Global WM underestimation was not significantly influenced. GM and WM volumes from NNI, did not differ significantly from gold-standard. GM segmentation differed between lesion masks in the lesion area, and also elsewhere. Using the gold-standard, FSL-FAST quantified as GM on average 0.4% of the lesion area with LEAP and 24.5% with FSL-lesion_filling. Lesion-filling did not influence DGM volumes from FSL-FIRST.

Discussion

These results demonstrate that for global GM volumetry, precise lesion masks on 3DT1 images can be replaced by co-registered 2DT2 lesion masks. This makes lesion-filling a feasible method for GM atrophy measurements in MS.

•

In multiple sclerosis brain atrophy measurement is affected by white matter lesions.

•

LEAP and FSL-lesion_filling replace lesion voxels with white matter intensities.

•

The two lesion-filling methods show differences.

•

2D lesion masks can be registered and used for lesion-filling on 3DT1 images.

•

This makes lesion-filling a feasible method for atrophy measurements in MS.

SELF-SUPERVISED MRI TISSUE SEGMENTATION BY DISCRIMINATIVE CLUSTERING

The study of brain lesions can benefit from a clear identification of transitions between healthy and pathological tissues, through the analysis of brain imaging data. Current signal processing methods, able to address these issues, often rely on strong prior information. In this article, a new method for tissue segmentation is proposed. It is based on a discriminative strategy, in a self-supervised machine learning approach. This method avoids the use of prior information, which makes it very versatile, and able to cope with different tissue types. It also returns tissue probabilities for each voxel, crucial for a good characterization of the evolution of brain lesions. Simulated as well as real benchmark data were used to validate the accuracy of the method and compare it against other segmentation algorithms.

Jacobian integration method increases the statistical power to measure gray matter atrophy in multiple sclerosis

Gray matter atrophy provides important insights into neurodegeneration in multiple sclerosis (MS) and can be used as a marker of neuroprotection in clinical trials. Jacobian integration is a method for measuring volume change that uses integration of the local Jacobian determinants of the nonlinear deformation field registering two images, and is a promising tool for measuring gray matter atrophy. Our main objective was to compare the statistical power of the Jacobian integration method to commonly used methods in terms of the sample size required to detect a treatment effect on gray matter atrophy. We used multi-center longitudinal data from relapsing-remitting MS patients and evaluated combinations of cross-sectional and longitudinal pre-processing with SIENAX/FSL, SPM, and FreeSurfer, as well as the Jacobian integration method. The Jacobian integration method outperformed these other commonly used methods, reducing the required sample size by a factor of 4-5. The results demonstrate the advantage of using the Jacobian integration method to assess neuroprotection in MS clinical trials.

Visual assessment of brain magnetic resonance imaging detects injury to cognitive regulatory sites in patients with heart failure

BACKGROUND

Heart failure (HF) patients exhibit depression and executive function impairments that contribute to HF mortality. Using specialized magnetic resonance imaging (MRI) analysis procedures, brain changes appear in areas regulating these functions (mammillary bodies, hippocampi, and frontal cortex). However, specialized MRI procedures are not part of standard clinical assessment for HF (which is usually a visual evaluation), and it is unclear whether visual MRI examination can detect changes in these structures.

METHODS AND RESULTS

Using brain MRI, we visually examined the mammillary bodies and frontal cortex for global and hippocampi for global and regional tissue changes in 17 HF and 50 control subjects. Significantly global changes emerged in the right mammillary body (HF 1.18 ± 1.13 vs control 0.52 ± 0.74; P = .024), right hippocampus (HF 1.53 ± 0.94 vs control 0.80 ± 0.86; P = .005), and left frontal cortex (HF 1.76 ± 1.03 vs control 1.24 ± 0.77; P = .034). Comparison of the visual method with specialized MRI techniques corroborates right hippocampal and left frontal cortical, but not mammillary body, tissue changes.

CONCLUSIONS

Visual examination of brain MRI can detect damage in HF in areas regulating depression and executive function, including the right hippocampus and left frontal cortex. Visual MRI assessment in HF may facilitate evaluation of injury to these structures and the assessment of the impact of potential treatments for this damage.

Validation of the automated method VIENA: an accurate, precise, and robust measure of ventricular enlargement

BACKGROUND

In many retrospective studies and large clinical trials, high-resolution, good-contrast 3DT1 images are unavailable, hampering detailed analysis of brain atrophy. Ventricular enlargement then provides a sensitive indirect measure of ongoing central brain atrophy. Validated automated methods are required that can reliably measure ventricular enlargement and are robust across magnetic resonance (MR) image types.

AIM

To validate the automated method VIENA for measuring the percentage ventricular volume change (PVVC) between two scans.

MATERIALS AND METHODS

Accuracy was assessed using four image types, acquired in 15 elderly patients (five with Alzheimer's disease, five with mild cognitive impairment, and five cognitively normal elderly) and 58 patients with multiple sclerosis (MS), by comparing PVVC values from VIENA to manual outlining. Precision was assessed from data with three imaging time points per MS patient, by measuring the difference between the direct (one-step) and indirect (two-step) measurement of ventricular volume change between the first and last time points. The stringent concordance correlation coefficient (CCC) was used to quantify absolute agreement.

RESULTS

CCC of VIENA with manual measurement was 0.84, indicating good absolute agreement. The median absolute difference between two-step and one-step measurement with VIENA was 1.01%, while CCC was 0.98. Neither initial ventricular volume nor ventricular volume change affected performance of the method.

DISCUSSION

VIENA has good accuracy and good precision across four image types. VIENA therefore provides a useful fully automated method for measuring ventricular volume change in large datasets.

CONCLUSION

VIENA is a robust, accurate, and precise method for measuring ventricular volume change.

A comprehensive approach to the segmentation of multichannel three-dimensional MR brain images in multiple sclerosis

Accurate classification and quantification of brain tissues is important for monitoring disease progression, measurement of atrophy, and correlating magnetic resonance (MR) measures with clinical disability. Classification of MR brain images in the presence of lesions, such as multiple sclerosis (MS), is particularly challenging. Images obtained with lower resolution often suffer from partial volume averaging leading to false classifications. While partial volume averaging can be reduced by acquiring volumetric images at high resolution, image segmentation and quantification can be technically challenging. In this study, we integrated the brain anatomical knowledge with non-parametric and parametric statistical classifiers for automatically classifying tissues and lesions on high resolution multichannel three-dimensional images acquired on 60 MS brains. The results of automatic lesion segmentation were reviewed by the expert. The agreement between results obtained by the automated analysis and the expert was excellent as assessed by the quantitative metrics, low absolute volume difference percent (36.18 ± 34.90), low average symmetric surface distance (1.64 mm ± 1.30 mm), high true positive rate (84.75 ± 12.69), and low false positive rate (34.10 ± 16.00). The segmented results were also in close agreement with the corrected results as assessed by Bland-Altman and regression analyses. Finally, our lesion segmentation was validated using the MS lesion segmentation grand challenge dataset (MICCAI 2008).

Grey matter damage in multiple sclerosis: a pathology perspective

Over the past decade, immunohistochemical studies have provided compelling evidence that gray matter (GM) pathology in multiple sclerosis (MS) is extensive. Until recently, this GM pathology was difficult to visualize using standard magnetic resonance imaging (MRI) techniques. However, with newly developed MRI sequences, it has become clear that GM damage is present from the earliest stages of the disease and accrues with disease progression. GM pathology is clinically relevant, as GM lesions and/or GM atrophy were shown to be associated with MS motor deficits and cognitive impairment. Recent autopsy studies demonstrated significant GM demyelination and microglia activation. However, extensive immune cell influx, complement activation and blood-brain barrier leakage, like in WM pathology, are far less prominent in the GM. Hence, so far, the cause of GM damage in MS remains unknown, although several plausible underlying pathogenic mechanisms have been proposed. This paper provides an overview of GM damage in MS with a focus on its topology and histopathology.

Recommendations to improve imaging and analysis of brain lesion load and atrophy in longitudinal studies of multiple sclerosis

Focal lesions and brain atrophy are the most extensively studied aspects of multiple sclerosis (MS), but the image acquisition and analysis techniques used can be further improved, especially those for studying within-patient changes of lesion load and atrophy longitudinally. Improved accuracy and sensitivity will reduce the numbers of patients required to detect a given treatment effect in a trial, and ultimately, will allow reliable characterization of individual patients for personalized treatment. Based on open issues in the field of MS research, and the current state of the art in magnetic resonance image analysis methods for assessing brain lesion load and atrophy, this paper makes recommendations to improve these measures for longitudinal studies of MS. Briefly, they are (1) images should be acquired using 3D pulse sequences, with near-isotropic spatial resolution and multiple image contrasts to allow more comprehensive analyses of lesion load and atrophy, across timepoints. Image artifacts need special attention given their effects on image analysis results. (2) Automated image segmentation methods integrating the assessment of lesion load and atrophy are desirable. (3) A standard dataset with benchmark results should be set up to facilitate development, calibration, and objective evaluation of image analysis methods for MS.

An open source multivariate framework for n-tissue segmentation with evaluation on public data

We introduce Atropos, an ITK-based multivariate n-class open source segmentation algorithm distributed with ANTs ( http://www.picsl.upenn.edu/ANTs). The Bayesian formulation of the segmentation problem is solved using the Expectation Maximization (EM) algorithm with the modeling of the class intensities based on either parametric or non-parametric finite mixtures. Atropos is capable of incorporating spatial prior probability maps (sparse), prior label maps and/or Markov Random Field (MRF) modeling. Atropos has also been efficiently implemented to handle large quantities of possible labelings (in the experimental section, we use up to 69 classes) with a minimal memory footprint. This work describes the technical and implementation aspects of Atropos and evaluates its performance on two different ground-truth datasets. First, we use the BrainWeb dataset from Montreal Neurological Institute to evaluate three-tissue segmentation performance via (1) K-means segmentation without use of template data; (2) MRF segmentation with initialization by prior probability maps derived from a group template; (3) Prior-based segmentation with use of spatial prior probability maps derived from a group template. We also evaluate Atropos performance by using spatial priors to drive a 69-class EM segmentation problem derived from the Hammers atlas from University College London. These evaluation studies, combined with illustrative examples that exercise Atropos options, demonstrate both performance and wide applicability of this new platform-independent open source segmentation tool.

Corpus callosum damage predicts disability progression and cognitive dysfunction in primary-progressive MS after five years

We aim to identify specific areas of white matter (WM) and grey matter (GM), which predict disability progression and cognitive dysfunction after five years in patients with primary-progressive multiple sclerosis (PPMS). Thirty-two patients with early PPMS were assessed at baseline and after five years on the Expanded Disability Status Scale (EDSS), and EDSS step-changes were calculated. At year five, a subgroup of 25 patients and 31 healthy controls underwent a neuropsychological assessment. Baseline imaging consisted of dual-echo (proton density and T2-weighted), T1-weighted volumetric, and diffusion tensor imaging. Fractional anisotropy (FA) maps were created, and fed into tract-based spatial statistics. To compensate for the potential bias introduced by WM lesions, the T1 volumes underwent a lesion-filling procedure before entering a voxel-based morphometry protocol. To investigate whether FA and GM volume predicted EDSS step-changes over five years and neuropsychological tests scores at five years, voxelwise linear regression analyses were performed. Lower FA in the splenium of the corpus callosum (CC) predicted a greater progression of disability over the follow-up. Lower FA along the entire CC predicted worse verbal memory, attention and speed of information processing, and executive function at five years. GM baseline volume did not predict any clinical variable. Our findings highlight the importance of damage to the interhemispheric callosal pathways in determining physical and cognitive disability in PPMS. Disruption of these pathways, which interconnect motor and cognitive networks between the two hemispheres, may result in a disconnection syndrome that contributes to long-term physical and cognitive disability.

Multivariate pattern classification of gray matter pathology in multiple sclerosis

Univariate analyses have identified gray matter (GM) alterations in different groups of MS patients. While these methods detect differences on the basis of the single voxel or cluster, multivariate methods like support vector machines (SVM) identify the complex neuroanatomical patterns of GM differences. Using multivariate linear SVM analysis and leave-one-out cross-validation, we aimed at identifying neuroanatomical GM patterns relevant for individual classification of MS patients. We used SVM to separate GM segmentations of T1-weighted three-dimensional magnetic resonance (MR) imaging scans within different age- and sex-matched groups of MS patients with either early (n=17) or late MS (n=17) (contrast I), low (n=20) or high (n=20) white matter lesion load (contrast II), and benign MS (BMS, n=13) or non-benign MS (NBMS, n=13) (contrast III) scanned on a single 1.5 T MR scanner. GM patterns most relevant for individual separation of MS patients comprised cortical areas of all the cerebral lobes as well as deep GM structures, including the thalamus and caudate. The patterns detected were sufficiently informative to separate individuals of the respective groups with high sensitivity and specificity in 85% (contrast I), 83% (contrast II) and 77% (contrast III) of cases. The study demonstrates that neuroanatomical spatial patterns of GM segmentations contain information sufficient for correct classification of MS patients at the single case level, thus making multivariate SVM analysis a promising clinical application.

The effect of hypointense white matter lesions on automated gray matter segmentation in multiple sclerosis

Previous imaging studies assessing the relationship between white matter (WM) damage and matter (GM) atrophy have raised the concern that Multiple Sclerosis (MS) WM lesions may affect measures of GM volume by inducing voxel misclassification during intensity-based tissue segmentation. Here, we quantified this misclassification error in simulated and real MS brains using a lesion-filling method. Using this method, we also corrected GM measures in patients before comparing them with controls in order to assess the impact of this lesion-induced misclassification error in clinical studies. We found that higher WM lesion volumes artificially reduced total GM volumes. In patients, this effect was about 72% of that predicted by simulation. Misclassified voxels were located at the GM/WM border and could be distant from lesions. Volume of individual deep gray matter (DGM) structures generally decreased with higher lesion volumes, consistent with results from total GM. While preserving differences in GM volumes between patients and controls, lesion-filling correction revealed more lateralised DGM shape changes in patients, which were not evident with the original images. Our results confirm that WM lesions can influence MRI measures of GM volume and shape in MS patients through their effect on intensity-based GM segmentation. The greater effect of lesions at increasing levels of damage supports the use of lesion-filling to correct for this problem and improve the interpretability of the results. Volumetric or morphometric imaging studies, where lesion amount and characteristics may vary between groups of patients or change over time, may especially benefit from this correction.

Evaluating and reducing the impact of white matter lesions on brain volume measurements

MR-based measurements of brain volumes may be affected by the presence of white matter (WM) lesions. Here, we assessed how and to what extent this may happen for WM lesions of various sizes and intensities. After inserting WM lesions of different sizes and intensities into T1-W brain images of healthy subjects, we assessed the effect on two widely used automatic methods for brain volume measurement such as SIENAX (segmentation-based) and SIENA (registration-based). To explore the relevance of partial volume (PV) estimation, we performed the experiments with two different PV models, implemented by the same segmentation algorithm (FAST) of SIENAX and SIENA. Finally, we tested potential solutions to this issue. The presence of WM lesions did not bias measurements for registration-based method such as SIENA. By contrast, the presence of WM lesions affected segmentation-based brain volume measurements such as SIENAx. The misclassification of both gray matter (GM) and WM volumes varied considerably with lesion size and intensity, especially when the lesion intensity was similar to that of the GM/WM interface. The extent to which the presence of WM lesions could affect tissue-class measures was clearly driven by the PV modeling used, with the mixel-type PV model giving a lower error in the presence of WM lesions. The tissue misclassification due to WM lesions was still present when they were masked out. By contrast, refilling the lesions with intensities matching the surrounding normal-appearing WM ensured accurate tissue-class measurements and thus represents a promising approach for accurate tissue classification and brain volume measurements.

Spatiotemporal distribution pattern of white matter lesion volumes and their association with regional grey matter volume reductions in relapsing-remitting multiple sclerosis

The association of white matter (WM) lesions and grey matter (GM) atrophy is a feature in relapsing-remitting multiple sclerosis (RRMS). The spatiotemporal distribution pattern of WM lesions, their relations to regional GM changes and the underlying dynamics are unclear. Here we combined parametric and non-parametric voxel-based morphometry (VBM) to clarify these issues. MRI data from RRMS patients with progressive (PLV, n = 45) and non-progressive WM lesion volumes (NPLV, n = 44) followed up for 12 months were analysed. Cross-sectionally, the spatial WM lesion distribution was compared using lesion probability maps (LPMs). Longitudinally, WM lesions and GM volumes were studied using FSL-VBM and SPM5-VBM, respectively. WM lesions clustered around the lateral ventricles and in the centrum semiovale with a more widespread pattern in the PLV than in the NPLV group. The maximum local probabilities were similar in both groups and higher for T2 lesions (PLV: 27%, NPLV: 25%) than for T1 lesions (PLV: 15%, NPLV 14%). Significant WM lesion changes accompanied by cortical GM volume reductions occurred in the corpus callosum and optic radiations (P = 0.01 corrected), and more liberally tested (uncorrected P < 0.01) in the inferior fronto-occipital and longitudinal fasciculi, and corona radiata in the PLV group. Not any WM or GM changes were found in the NPLV group. In the PLV group, WM lesion distribution and development in fibres, was associated with regional GM volume loss. The different spatiotemporal distribution patterns of patients with progressive compared to patients with non-progressive WM lesions suggest differences in the dynamics of pathogenesis.

Longitudinal gray matter changes in multiple sclerosis--differential scanner and overall disease-related effects

Voxel-based morphometry (VBM) has been used repeatedly in single-center studies to investigate regional gray matter (GM) atrophy in multiple sclerosis (MS). In multi-center trials, across-scanner variations might interfere with the detection of disease-specific structural abnormalities, thereby potentially limiting the use of VBM. Here we evaluated longitudinally inter-site differences and inter-site comparability of regional GM in MS using VBM. Baseline and follow up 3D T1-weighted magnetic resonance imaging (MRI) data of 248 relapsing-remitting (RR) MS patients, recruited in two clinical centers, (center1/2: n = 129/119; mean age 42.6 ± 10.7/43.3 ± 9.3; male:female 33:96/44:75; median disease duration 150 [72-222]/116 [60-156]) were acquired on two different 1.5T MR scanners. GM volume changes between baseline and year 2 while controlling for age, gender, disease duration, and global GM volume were analyzed. The main effect of time on regional GM volume was larger in data of center two as compared to center one in most of the brain regions. Differential effects of GM volume reductions occurred in a number of GM regions of both hemispheres, in particular in the fronto-temporal and limbic cortex (cluster P corrected <0.05). Overall disease-related effects were found bilaterally in the cerebellum, uncus, inferior orbital gyrus, paracentral lobule, precuneus, inferior parietal lobule, and medial frontal gyrus (cluster P corrected <0.05). The differential effects were smaller as compared to the overall effects in these regions. These results suggest that the effects of different scanners on longitudinal GM volume differences were rather small and thus allow pooling of MR data and subsequent combined image analysis.