Treatment reduces the incidence of newly appearing multiple sclerosis lesions evolving into chronic active, slowly expanding lesions: A retrospective analysis.

BACKGROUND AND PURPOSE: Newly appearing lesions in multiple sclerosis (MS) may evolve into chronically active, slowly expanding lesions (SELs), leading to sustained disability progression. The aim of this study was to evaluate the incidence of newly appearing lesions developing into SELs, and their correlation to clinical evolution and treatment. METHODS: A retrospective analysis of a fingolimod trial in primary progressive MS (PPMS; INFORMS, NCT00731692) was undertaken. Data were available from 324 patients with magnetic resonance imaging scans up to 3 years after screening. New lesions at year 1 were identified with convolutional neural networks, and SELs obtained through a deformation-based method. Clinical disability was assessed annually by Expanded Disability Status Scale (EDSS), Nine-Hole Peg Test, Timed 25-Foot Walk, and Paced Auditory Serial Addition Test. Linear, logistic, and mixed-effect models were used to assess the relationship between the Jacobian expansion in new lesions and SELs, disability scores, and treatment status. RESULTS: One hundred seventy patients had ≥1 new lesions at year 1 and had a higher lesion count at screening compared to patients with no new lesions (median = 27 vs. 22, p = 0.007). Among the new lesions (median = 2 per patient), 37% evolved into definite or possible SELs. Higher SEL volume and count were associated with EDSS worsening and confirmed disability progression. Treated patients had lower volume and count of definite SELs (ß = -0.04, 95% confidence interval [CI] = -0.07 to -0.01, p = 0.015; ß = -0.36, 95% CI = -0.67 to -0.06, p = 0.019, respectively). CONCLUSIONS: Incident chronic active lesions are common in PPMS, and fingolimod treatment can reduce their number.

SENSE EPI reconstruction with 2D phase error correction and channel-wise noise removal.

PURPOSE: To develop a robust reconstruction pipeline for EPI data that enables 2D Nyquist phase error correction using sensitivity encoding without incurring major noise artifacts in low SNR data. METHODS: SENSE with 2D phase error correction (PEC-SENSE) was combined with channel-wise noise removal using Marcenko-Pastur principal component analysis (MPPCA) to simultaneously eliminate Nyquist ghost artifacts in EPI data and mitigate the noise amplification associated with phase correction using parallel imaging. The proposed pipeline (coined SPECTRE) was validated in phantom DW-EPI data using the accuracy and precision of diffusion metrics; ground truth values were obtained from data acquired with a spin echo readout. Results from the SPECTRE pipeline were compared against PEC-SENSE reconstructions with three alternate denoising strategies: (i) no denoising; (ii) denoising of magnitude data after image formation; (iii) denoising of complex data after image formation. SPECTRE was then tested using high b $$ b $$ -value (i.e., low SNR) diffusion data (up to b = 3000 $$ b=3000 $$ s/mm 2 $$ {}^2 $$ ) in four healthy subjects. RESULTS: Noise amplification associated with phase error correction incurred a 23% bias in phantom mean diffusivity (MD) measurements. Phantom MD estimates using the SPECTRE pipeline were within 8% of the ground truth value. In healthy volunteers, the SPECTRE pipeline visibly corrected Nyquist ghost artifacts and reduced associated noise amplification in high b $$ b $$ -value data. CONCLUSION: The proposed reconstruction pipeline is effective in correcting low SNR data, and improves the accuracy and precision of derived diffusion metrics.

Generic acquisition protocol for quantitative MRI of the spinal cord.

Quantitative spinal cord (SC) magnetic resonance imaging (MRI) presents many challenges, including a lack of standardized imaging protocols. Here we present a prospectively harmonized quantitative MRI protocol, which we refer to as the spine generic protocol, for users of 3T MRI systems from the three main manufacturers: GE, Philips and Siemens. The protocol provides guidance for assessing SC macrostructural and microstructural integrity: T1-weighted and T2-weighted imaging for SC cross-sectional area computation, multi-echo gradient echo for gray matter cross-sectional area, and magnetization transfer and diffusion weighted imaging for assessing white matter microstructure. In a companion paper from the same authors, the spine generic protocol was used to acquire data across 42 centers in 260 healthy subjects. The key details of the spine generic protocol are also available in an open-access document that can be found at https://github.com/spine-generic/protocols . The protocol will serve as a starting point for researchers and clinicians implementing new SC imaging initiatives so that, in the future, inclusion of the SC in neuroimaging protocols will be more common. The protocol could be implemented by any trained MR technician or by a researcher/clinician familiar with MRI acquisition.

Tracking White and Gray Matter Degeneration along the Spinal Cord Axis in Degenerative Cervical Myelopathy.

This study aims to determine tissue-specific neurodegeneration across the spinal cord in patients with mild-moderate degenerative cervical myelopathy (DCM). Twenty-four mild-moderate DCM and 24 healthy subjects were recruited. In patients, a T2-weighted scan was acquired at the compression site, whereas in all participants a T2*-weighted and diffusion-weighted scan was acquired at the cervical level (C2-C3) and in the lumbar enlargement (i.e., rostral and caudal to the site of compression). We quantified intramedullary signal changes, maximal canal and cord compression, white (WM) and gray matter (GM) atrophy, and microstructural indices from diffusion-weighted scans. All patients underwent clinical (modified Japanese Orthopaedic Association; mJOA) and electrophysiological assessments. Regression analysis assessed associations between magnetic resonance imaging (MRI) readouts and electrophysiological and clinical outcomes. Twenty patients were classified with mild and 4 with moderate DCM using the mJOA scale. The most frequent site of compression was at the C5-C6 level, with maximum cord compression of 38.73% ± 11.57%. Ten patients showed imaging evidence of cervical myelopathy. In the cervical cord, WM and GM atrophy and WM microstructural changes were evident, whereas in the lumbar cord only WM showed atrophy and microstructural changes. Remote cervical cord WM microstructural changes were pronounced in patients with radiological myelopathy and associated with impaired electrophysiology. Lumbar cord WM atrophy was associated with lower limb sensory impairments. In conclusion, tissue-specific neurodegeneration revealed by quantitative MRI is already apparent across the spinal cord in mild-moderate DCM before the onset of severe clinical impairments. WM microstructural changes are particularly sensitive to remote pathologically and clinically eloquent changes in DCM.

Brain microstructural and metabolic alterations detected in vivo at onset of the first demyelinating event.

In early multiple sclerosis, a clearer understanding of normal-brain tissue microstructural and metabolic abnormalities will provide valuable insights into its pathophysiology. We used multi-parametric quantitative MRI to detect alterations in brain tissues of patients with their first demyelinating episode. We acquired neurite orientation dispersion and density imaging [to investigate morphology of neurites (dendrites and axons)] and 23Na MRI (to estimate total sodium concentration, a reflection of underlying changes in metabolic function). In this cross-sectional study, we enrolled 42 patients diagnosed with clinically isolated syndrome or multiple sclerosis within 3 months of their first demyelinating event and 16 healthy controls. Physical and cognitive scales were assessed. At 3 T, we acquired brain and spinal cord structural scans, and neurite orientation dispersion and density imaging. Thirty-two patients and 13 healthy controls also underwent brain 23Na MRI. We measured neurite density and orientation dispersion indices and total sodium concentration in brain normal-appearing white matter, white matter lesions, and grey matter. We used linear regression models (adjusting for brain parenchymal fraction and lesion load) and Spearman correlation tests (significance level P ≤ 0.01). Patients showed higher orientation dispersion index in normal-appearing white matter, including the corpus callosum, where they also showed lower neurite density index and higher total sodium concentration, compared with healthy controls. In grey matter, compared with healthy controls, patients demonstrated: lower orientation dispersion index in frontal, parietal and temporal cortices; lower neurite density index in parietal, temporal and occipital cortices; and higher total sodium concentration in limbic and frontal cortices. Brain volumes did not differ between patients and controls. In patients, higher orientation dispersion index in corpus callosum was associated with worse performance on timed walk test (P = 0.009, B = 0.01, 99% confidence interval = 0.0001 to 0.02), independent of brain and lesion volumes. Higher total sodium concentration in left frontal middle gyrus was associated with higher disability on Expanded Disability Status Scale (rs = 0.5, P = 0.005). Increased axonal dispersion was found in normal-appearing white matter, particularly corpus callosum, where there was also axonal degeneration and total sodium accumulation. The association between increased axonal dispersion in the corpus callosum and worse walking performance implies that morphological and metabolic alterations in this structure could mechanistically contribute to disability in multiple sclerosis. As brain volumes were neither altered nor related to disability in patients, our findings suggest that these two advanced MRI techniques are more sensitive at detecting clinically relevant pathology in early multiple sclerosis.

White matter integrity correlates with cognition and disease severity in Fabry disease.

Cerebral white matter pathology is a common CNS manifestation of Fabry disease, visualized as white matter hyperintensities on MRI in 42-81% of patients. Diffusion tensor imaging (DTI) MRI is a sensitive technique to quantify microstructural damage within the white matter with potential value as a disease biomarker. We evaluated the pattern of DTI abnormalities in Fabry disease, and their correlations with cognitive impairment, mood, anxiety, disease severity and plasma lyso-Gb3 levels in 31 patients with genetically proven Fabry disease and 19 age-matched healthy control subjects. We obtained average values of fractional anisotropy and mean diffusivity within the white matter and performed voxelwise analysis with tract-based spatial statistics. Using a standardized neuropsychological test battery, we assessed processing speed, executive function, anxiety, depression and disease severity. The mean age (% male) was 44.1 (45%) for patients with Fabry disease and 37.4 (53%) for the healthy control group. In patients with Fabry disease, compared to healthy controls the mean average white matter fractional anisotropy was lower in [0.423 (standard deviation, SD 0.023) versus 0.446 (SD 0.016), P = 0.002] while mean average white matter mean diffusivity was higher (749 × 10-6 mm2/s (SD 32 × 10-6) versus 720 × 10-6 mm2/s (SD 21 × 10-6), P = 0.004]. Voxelwise statistics showed that the diffusion abnormalities for both fractional anisotropy and mean diffusivity were anatomically widespread. A lesion probability map showed that white matter hyperintensities also had a wide anatomical distribution with a predilection for the posterior centrum semiovale. However, diffusion abnormalities in Fabry disease were not restricted to lesional tissue; compared to healthy controls, the normal appearing white matter in patients with Fabry disease had reduced fractional anisotropy [0.422 (SD 0.022) versus 0.443 (SD 0.017) P = 0.003] and increased mean diffusivity [747 × 10-6 mm2/s (SD 26 × 10-6) versus 723 × 10-6 mm2/s (SD 22 × 10-6), P = 0.008]. Within patients, average white matter fractional anisotropy and white matter lesion volume showed statistically significant correlations with Digit Symbol Coding Test score (r = 0.558, P = 0.001; and r = -0.633, P ≤ 0.001, respectively). Average white matter fractional anisotropy correlated with the overall Mainz Severity Score Index (r = -0.661, P ≤ 0.001), while average white matter mean diffusivity showed a strong correlation with plasma lyso-Gb3 levels (r = 0.559, P = 0.001). Our findings using DTI confirm widespread areas of microstructural white matter disruption in Fabry disease, extending beyond white matter hyperintensities seen on conventional MRI. Moreover, diffusion measures show strong correlations with cognition (processing speed), clinical disease severity and a putative plasma biomarker of disease activity, making them promising quantitative biomarkers for monitoring Fabry disease severity and progression.

Disrupted principal network organisation in multiple sclerosis relates to disability.

Structural network-based approaches can assess white matter connections revealing topological alterations in multiple sclerosis (MS). However, principal network (PN) organisation and its clinical relevance in MS has not been explored yet. Here, structural networks were reconstructed from diffusion data in 58 relapsing-remitting MS (RRMS), 28 primary progressive MS (PPMS), 36 secondary progressive (SPMS) and 51 healthy controls (HCs). Network hubs' strengths were compared with HCs. Then, PN analysis was performed in each clinical subtype. Regression analysis was applied to investigate the associations between nodal strength derived from the first and second PNs (PN1 and PN2) in MS, with clinical disability. Compared with HCs, MS patients had preserved hub number, but some hubs exhibited reduced strength. PN1 comprised 10 hubs in HCs, RRMS and PPMS but did not include the right thalamus in SPMS. PN2 comprised 10 hub regions with intra-hemispheric connections in HCs. In MS, this subnetwork did not include the right putamen whilst in SPMS the right thalamus was also not included. Decreased nodal strength of the right thalamus and putamen from the PNs correlated strongly with higher clinical disability. These PN analyses suggest distinct patterns of disruptions in MS subtypes which are clinically relevant.

Characterisation of tissue-type metabolic content in secondary progressive multiple sclerosis: a magnetic resonance spectroscopic imaging study.

Proton magnetic resonance spectroscopy yields metabolic information and has proved to be a useful addition to structural imaging in neurological diseases. We applied short-echo time Spectroscopic Imaging in a cohort of 42 patients with secondary progressive multiple sclerosis (SPMS). Linear modelling with respect to brain tissue type yielded metabolite levels that were significantly different in white matter lesions compared with normal-appearing white matter, suggestive of higher myelin turnover (higher choline), higher metabolic rate (higher creatine) and increased glial activity (higher myo-inositol) within the lesions. These findings suggest that the lesions have ongoing cellular activity that is not consistent with the usual assumption of 'chronic' lesions in SPMS, and may represent a target for repair therapies.

Deep gray matter volume loss drives disability worsening in multiple sclerosis.

OBJECTIVE: Gray matter (GM) atrophy occurs in all multiple sclerosis (MS) phenotypes. We investigated whether there is a spatiotemporal pattern of GM atrophy that is associated with faster disability accumulation in MS. METHODS: We analyzed 3,604 brain high-resolution T1-weighted magnetic resonance imaging scans from 1,417 participants: 1,214 MS patients (253 clinically isolated syndrome [CIS], 708 relapsing-remitting [RRMS], 128 secondary-progressive [SPMS], and 125 primary-progressive [PPMS]), over an average follow-up of 2.41 years (standard deviation [SD] = 1.97), and 203 healthy controls (HCs; average follow-up = 1.83 year; SD = 1.77), attending seven European centers. Disability was assessed with the Expanded Disability Status Scale (EDSS). We obtained volumes of the deep GM (DGM), temporal, frontal, parietal, occipital and cerebellar GM, brainstem, and cerebral white matter. Hierarchical mixed models assessed annual percentage rate of regional tissue loss and identified regional volumes associated with time-to-EDSS progression. RESULTS: SPMS showed the lowest baseline volumes of cortical GM and DGM. Of all baseline regional volumes, only that of the DGM predicted time-to-EDSS progression (hazard ratio = 0.73; 95% confidence interval, 0.65, 0.82; p < 0.001): for every standard deviation decrease in baseline DGM volume, the risk of presenting a shorter time to EDSS worsening during follow-up increased by 27%. Of all longitudinal measures, DGM showed the fastest annual rate of atrophy, which was faster in SPMS (-1.45%), PPMS (-1.66%), and RRMS (-1.34%) than CIS (-0.88%) and HCs (-0.94%; p < 0.01). The rate of temporal GM atrophy in SPMS (-1.21%) was significantly faster than RRMS (-0.76%), CIS (-0.75%), and HCs (-0.51%). Similarly, the rate of parietal GM atrophy in SPMS (-1.24-%) was faster than CIS (-0.63%) and HCs (-0.23%; all p values <0.05). Only the atrophy rate in DGM in patients was significantly associated with disability accumulation (beta = 0.04; p < 0.001). INTERPRETATION: This large, multicenter and longitudinal study shows that DGM volume loss drives disability accumulation in MS, and that temporal cortical GM shows accelerated atrophy in SPMS than RRMS. The difference in regional GM atrophy development between phenotypes needs to be taken into account when evaluating treatment effect of therapeutic interventions. Ann Neurol 2018;83:210-222.

Spinal cord grey matter segmentation challenge.

An important image processing step in spinal cord magnetic resonance imaging is the ability to reliably and accurately segment grey and white matter for tissue specific analysis. There are several semi- or fully-automated segmentation methods for cervical cord cross-sectional area measurement with an excellent performance close or equal to the manual segmentation. However, grey matter segmentation is still challenging due to small cross-sectional size and shape, and active research is being conducted by several groups around the world in this field. Therefore a grey matter spinal cord segmentation challenge was organised to test different capabilities of various methods using the same multi-centre and multi-vendor dataset acquired with distinct 3D gradient-echo sequences. This challenge aimed to characterize the state-of-the-art in the field as well as identifying new opportunities for future improvements. Six different spinal cord grey matter segmentation methods developed independently by various research groups across the world and their performance were compared to manual segmentation outcomes, the present gold-standard. All algorithms provided good overall results for detecting the grey matter butterfly, albeit with variable performance in certain quality-of-segmentation metrics. The data have been made publicly available and the challenge web site remains open to new submissions. No modifications were introduced to any of the presented methods as a result of this challenge for the purposes of this publication.

Machine learning based compartment models with permeability for white matter microstructure imaging.

Some microstructure parameters, such as permeability, remain elusive because mathematical models that express their relationship to the MR signal accurately are intractable. Here, we propose to use computational models learned from simulations to estimate these parameters. We demonstrate the approach in an example which estimates water residence time in brain white matter. The residence time τi of water inside axons is a potentially important biomarker for white matter pathologies of the human central nervous system, as myelin damage is hypothesised to affect axonal permeability, and thus τi. We construct a computational model using Monte Carlo simulations and machine learning (specifically here a random forest regressor) in order to learn a mapping between features derived from diffusion weighted MR signals and ground truth microstructure parameters, including τi. We test our numerical model using simulated and in vivo human brain data. Simulation results show that estimated parameters have strong correlations with the ground truth parameters (R2={0.88,0.95,0.82,0.99}) for volume fraction, residence time, axon radius and diffusivity respectively), and provide a marked improvement over the most widely used Kärger model (R2={0.75,0.60,0.11,0.99}). The trained model also estimates sensible microstructure parameters from in vivo human brain data acquired from healthy controls, matching values found in literature, and provides better reproducibility than the Kärger model on both the voxel and ROI level. Finally, we acquire data from two Multiple Sclerosis (MS) patients and compare to the values in healthy subjects. We find that in the splenium of corpus callosum (CC-S) the estimate of the residence time is 0.57±0.05s for the healthy subjects, while in the MS patient with a lesion in CC-S it is 0.33±0.12s in the normal appearing white matter (NAWM) and 0.19±0.11s in the lesion. In the corticospinal tracts (CST) the estimate of the residence time is 0.52±0.09s for the healthy subjects, while in the MS patient with a lesion in CST it is 0.56±0.05s in the NAWM and 0.13±0.09s in the lesion. These results agree with our expectations that the residence time in lesions would be lower than in NAWM because the loss of myelin should increase permeability. Overall, we find parameter estimates in the two MS patients consistent with expectations from the pathology of MS lesions demonstrating the clinical potential of this new technique.

Fully automated grey and white matter spinal cord segmentation.

Axonal loss in the spinal cord is one of the main contributing factors to irreversible clinical disability in multiple sclerosis (MS). In vivo axonal loss can be assessed indirectly by estimating a reduction in the cervical cross-sectional area (CSA) of the spinal cord over time, which is indicative of spinal cord atrophy, and such a measure may be obtained by means of image segmentation using magnetic resonance imaging (MRI). In this work, we propose a new fully automated spinal cord segmentation technique that incorporates two different multi-atlas segmentation propagation and fusion techniques: The Optimized PatchMatch Label fusion (OPAL) algorithm for localising and approximately segmenting the spinal cord, and the Similarity and Truth Estimation for Propagated Segmentations (STEPS) algorithm for segmenting white and grey matter simultaneously. In a retrospective analysis of MRI data, the proposed method facilitated CSA measurements with accuracy equivalent to the inter-rater variability, with a Dice score (DSC) of 0.967 at C2/C3 level. The segmentation performance for grey matter at C2/C3 level was close to inter-rater variability, reaching an accuracy (DSC) of 0.826 for healthy subjects and 0.835 people with clinically isolated syndrome MS.

A framework for optimal whole-sample histological quantification of neurite orientation dispersion in the human spinal cord.

BACKGROUND: The complexity of fibre distributions in tissues is an important microstructural feature, now measurable in vivo by magnetic resonance imaging (MRI) through orientation dispersion (OD) indices. OD metrics have gained popularity for the characterisation of neurite morphology, but they still lack systematic validation. This paper demonstrates a framework for whole-sample histological quantification of OD in spinal cord specimens, potentially useful for validating MRI-derived OD estimates. NEW METHOD: Our methodological framework is based on (i) sagittal sectioning; (ii) Palmgren's silver staining; (iii) structure tensor (ST) analysis; (iv) directional statistics. Novel elements are the data-driven optimisation of the spatial scale of ST analysis, and a new multivariate, weighted directional statistical approach for anisotropy-informed quantification of OD. RESULTS: Palmgren's silver staining of sagittal spinal cord sections provides robust visualisation of neuronal elements, enabling OD quantification. The choice of spatial scale of ST analysis influences OD values, and weighted directional statistics provide OD maps with high contrast-to-noise. Segmentation of neurites prior to OD quantification is recommended. COMPARISON WITH EXISTING METHODS: Our framework can potentially provide OD even in demyelinating diseases, where myelin-based histology is not suitable. As compared to conventional univariate approaches, our multivariate weighted directional statistics improve the contrast-to-noise of OD maps and more accurately describe the distribution of ST metrics. CONCLUSIONS: Our framework enables practical whole-specimen characterisation of OD in the spinal cord. We recommend tuning the scale of ST analysis for optimal OD quantification, as well as neurite segmentation and weighted directional statistics, of which examples are provided herein.