OBJECTIVE Pathological studies suggest that multiple sclerosis (MS) lesions endure multiple waves of damage and repair; however, the dynamics and characteristics of these processes are poorly understood in patients living with MS. METHODS We studied 128 MS patients (75 relapsing-remitting, 53 progressive) and 72 healthy controls who underwent advanced magnetic resonance imaging and clinical examination at baseline and 2 years later. Magnetization transfer saturation and multi-shell diffusion imaging were used to quantify longitudinal changes in myelin and axon volumes within MS lesions. Lesions were grouped into 4 classes (repair, damage, mixed repair damage, and stable). The frequency of each class was correlated to clinical measures, demographic characteristics, and levels of serum neurofilament light chain (sNfL). RESULTS Stable lesions were the most frequent (n = 2,276; 44%), followed by lesions with patterns of "repair" (n = 1,352; 26.2%) and damage (n = 1,214; 23.5%). The frequency of "repair" lesion was negatively associated with disability (β = -0.04; p < 0.001) and sNfL (β = -0.16; p < 0.001) at follow-up. The frequency of the "damage" class was higher in progressive than relapsing-remitting patients (p < 0.05) and was related to disability (baseline: β = -0.078; follow-up: β = -0.076; p < 0.001) and age (baseline: β = -0.078; p < 0.001). Stable lesions were more frequent in relapsing-remitting than in progressive patients (p < 0.05), and in younger patients versus older (β = -0.07; p < 0.001) at baseline. Further, "mixed" lesions were most frequent in older patients (β = 0.004; p < 0.001) at baseline. INTERPRETATION These findings show that repair and damage processes within MS lesions occur across the entire disease spectrum and that their frequency correlates with patients disability, age, disease duration, and extent of neuroaxonal damage. ANN NEUROL 2024.

Background and Objectives In patients with multiple sclerosis (PwMS), thalamic atrophy occurs during the disease course. However, there is little understanding of the mechanisms leading to volume loss and of the relationship between microstructural thalamic pathology and disease progression. This cross-sectional and longitudinal study aimed to comprehensively characterize in vivo pathologic changes within thalamic microstructure in PwMS using advanced multiparametric quantitative MRI (qMRI). Methods Thalamic microstructural integrity was evaluated using quantitative T1, magnetization transfer saturation, multishell diffusion, and quantitative susceptibility mapping (QSM) in 183 PwMS and 105 healthy controls (HCs). The same qMRI protocol was available for 127 PwMS and 73 HCs after a 2-year follow-up period. Inclusion criteria for PwMS encompassed either an active relapsing-remitting MS (RRMS) or inactive progressive MS (PMS) disease course. Thalamic alterations were compared between PwMS and HCs and among disease phenotypes. In addition, the study investigated the relationship between thalamic damage and clinical and conventional MRI measures of disease severity. Results Compared with HCs, PwMS exhibited substantial thalamic alterations, indicative of microstructural and macrostructural damage, demyelination, and disruption in iron homeostasis. These alterations extended beyond focal thalamic lesions, affecting normal-appearing thalamic tissue diffusely. Over the follow-up period, PwMS displayed an accelerated decrease in myelin volume fraction [mean difference in annualized percentage change (MD-ApC) = −1.50; p = 0.041] and increase in quantitative T1 (MD-ApC = 0.92; p < 0.0001) values, indicating heightened demyelinating and neurodegenerative processes. The observed differences between PwMS and HCs were substantially driven by the subgroup with PMS, wherein thalamic degeneration was significantly accelerated, even in comparison with patients with RRMS. Thalamic qMRI alterations showed extensive correlations with conventional MRI, clinical, and cognitive disease burden measures. Disability progression over follow-up was associated with accelerated thalamic degeneration, as reflected by enhanced diffusion (β = −0.067; p = 0.039) and QSM (β = −0.077; p = 0.027) changes. Thalamic qMRI metrics emerged as significant predictors of neurologic and cognitive disability even when accounting for other established markers including white matter lesion load and brain and thalamic atrophy. Discussion These findings offer deeper insights into thalamic pathology in PwMS, emphasizing the clinical relevance of thalamic damage and its link to disease progression. Advanced qMRI biomarkers show promising potential in guiding interventions aimed at mitigating thalamic neurodegenerative processes.

Background and Objectives Myelin and iron play essential roles in remyelination processes of multiple sclerosis (MS) lesions. χ-separation, a novel biophysical model applied to multiecho T2*-data and T2-data, estimates the contribution of myelin and iron to the obtained susceptibility signal. We used this method to investigate myelin and iron levels in lesion and nonlesion brain areas in patients with MS and healthy individuals. Methods This prospective MS cohort study included patients with MS fulfilling the McDonald Criteria 2017 and healthy individuals, aged 18 years or older, with no other neurologic comorbidities. Participants underwent MRI at baseline and after 2 years, including multiecho GRE-(T2*) and FAST-(T2) sequences. Using χ-separation, we generated myelin-sensitive and iron-sensitive susceptibility maps. White matter lesions (WMLs), cortical lesions (CLs), surrounding normal-appearing white matter (NAWM), and normal-appearing gray matter were segmented on fluid-attenuated inversion recovery and magnetization-prepared 2 rapid gradient echo images, respectively. Cross-sectional group comparisons used Wilcoxon rank-sum tests, longitudinal analyses applied Wilcoxon signed-rank tests. Associations with clinical outcomes (disease phenotype, age, sex, disease duration, disability measured by Expanded Disability Status Scale [EDSS], neurofilament light chain levels, and T2-lesion number and volume) were assessed using linear regression models. Results Of 168 patients with MS (median [interquartile range (IQR)] age 47.0 [21.7] years; 101 women; 6,898 WMLs, 775 CLs) and 103 healthy individuals (age 33.0 [10.5] years, 57 women), 108 and 62 were followed for a median of 2 years, respectively (IQR 0.1; 5,030 WMLs, 485 CLs). At baseline, WMLs had lower myelin (median 0.025 [IQR 0.015] parts per million [ppm]) and iron (0.017 [0.015] ppm) than the corresponding NAWM (myelin 0.030 [0.012]; iron 0.019 [0.011] ppm; both p < 0.001). After 2 years, both myelin (0.027 [0.014] ppm) and iron had increased (0.018 [0.015] ppm; both p < 0.001). Younger age (p < 0.001, b = −5.111 × 10−5), lower disability (p = 0.04, b = −2.352 × 10−5), and relapsing-remitting phenotype (RRMS, 0.003 [0.01] vs primary progressive 0.002 [IQR 0.01], p < 0.001; vs secondary progressive 0.0004 [IQR 0.01], p < 0.001) at baseline were associated with remyelination. Increment of myelin correlated with clinical improvement measured by EDSS (p = 0.015, b = −6.686 × 10−4). Discussion χ-separation, a novel mathematical model applied to multiecho T2*-images and T2-images shows that young RRMS patients with low disability exhibit higher remyelination capacity, which correlated with clinical disability over a 2-year follow-up.

PURPOSE To fully characterize the orientation dependence of magnetization transfer (MT) and inhomogeneous MT (ihMT) measures in the whole white matter (WM), for both single-fiber and crossing-fiber voxels. METHODS A characterization method was developed using the fiber orientation obtained from diffusion MRI (dMRI) with diffusion tensor imaging (DTI) and constrained spherical deconvolution. This allowed for characterization of the orientation dependence of measures in all of WM, regardless of the number of fiber orientation in a voxel. Furthermore, the orientation dependence inside 31 different WM bundles was characterized to evaluate the homogeneity of the effect. Variation of the results within and between-subject was assessed from a 12-subject dataset. RESULTS Previous results for single-fiber voxels were reproduced and a novel characterization was produced in voxels of crossing fibers, which seems to follow trends consistent with single-fiber results. Heterogeneity of the orientation dependence across bundles was observed, but homogeneity within similar bundles was also highlighted. Differences in behavior between MT and ihMT measures, as well as the ratio and saturation versions of these, were noted. CONCLUSION Orientation dependence characterization was proven possible over the entirety of WM. The vast range of effects and subtleties of the orientation dependence on MT measures showed the need for, but also the challenges of, a correction method.

Introduction Multi-shell diffusion Magnetic Resonance Imaging (dMRI) data has been widely used to characterise white matter microstructure in several neurodegenerative diseases. The lack of standardised dMRI protocols often implies the acquisition of redundant measurements, resulting in prolonged acquisition times. In this study, we investigate the impact of the number of gradient directions on Diffusion Tensor Imaging (DTI) and on Neurite Orientation Dispersion and Density Imaging (NODDI) metrics. Methods Data from 124 healthy controls collected in three different longitudinal studies were included. Using an in-house algorithm, we reduced the number of gradient directions in each data shell. We estimated DTI and NODDI measures on six white matter bundles clinically relevant for neurodegenerative diseases. Results Fractional Anisotropy (FA) measures on bundles where data were sampled at the 30% rate, showed a median L1 distance of up to 3.92% and a 95% CI of (1.74, 8.97)% when compared to those obtained at reference sampling. Mean Diffusivity (MD) reached up to 4.31% and a 95% CI of (1.60, 16.98)% on the same premises. At a sampling rate of 50%, we obtained a median of 3.90% and a 95% CI of (1.99, 16.65)% in FA, and 5.49% with a 95% CI of (2.14, 21.68)% in MD. The Intra-Cellular volume fraction (ICvf) median L1 distance was up to 2.83% with a 95% CI of (1.98, 4.82)% at a 30% sampling rate and 3.95% with a 95% CI of (2.39, 7.81)% at a 50% sampling rate. The volume difference of the reconstructed white matter at reference and 50% sampling reached a maximum of (2.09 ± 0.81)%. Discussion In conclusion, DTI and NODDI measures reported at reference sampling were comparable to those obtained when the number of dMRI volumes was reduced by up to 30%. Close to reference DTI and NODDI metrics were estimated with a significant reduction in acquisition time using three shells, respectively with: 4 directions at a b value of 700 s/mm2, 14 at 1000 s/mm2, and 32 at 2000 s/mm2. The study revealed aspects that can be important for large-scale clinical studies on bundle-specific diffusion MRI.

BACKGROUND AND OBJECTIVES A subgroup of patients with multiple sclerosis (MS) presents focal paramagnetic rims at the border between cortex and white matter (juxtacortical paramagnetic rims [JPRs]). We investigated the presence of this finding in our in vivo MS cohort and explored its potential clinical relevance. Moreover, we exploited postmortem MRI of fixed whole MS brains to (1) detect those rims and (2) investigate their histologic correlation. METHODS Quantitative susceptibility mapping (QSM) and magnetization-prepared 2 rapid acquisition gradient-echo (MP2RAGE) images at 3T-MRI of 165 patients with MS from the in vivo cohort were screened for JPRs and the presence of cortical lesions. Five postmortem brains from patients with MS were imaged with 3T-MRI to obtain QSM and MP2RAGE sequences. Tissue blocks containing JPRs were excised and paraffin-embedded slices stained by immunohistochemistry for myelin basic protein (for myelin) and anti-CR3/43 (for major histocompatibility complex II-positive microglia/macrophages). DAB-Turnbull stain was performed to detect iron. RESULTS JPRs are present in approximately 10% of in vivo patients and are associated with increased cortical lesion load. One of the 5 postmortem brains showed JPRs. Histologically, JPRs correspond to an accumulation of activated iron-laden phagocytes and are associated with demyelination of the whole overlying cortical ribbon. DISCUSSION JPRs are a novel potential MRI biomarker of focal cortical demyelination, which seems related to global cortical pathology and might be useful for diagnostic and stratification purposes in a clinical setting.

Introduction Recent studies showed that the myelin of the brain changes in the life span, and demyelination contributes to the loss of brain plasticity during normal aging. Diffusion-weighted magnetic resonance imaging (dMRI) allows studying brain connectivity in vivo by mapping axons in white matter with tractography algorithms. However, dMRI does not provide insight into myelin; thus, combining tractography with myelin-sensitive maps is necessary to investigate myelin-weighted brain connectivity. Tractometry is designated for this purpose, but it suffers from some serious limitations. Our study assessed the effectiveness of the recently proposed Myelin Streamlines Decomposition (MySD) method in estimating myelin-weighted connectomes and its capacity to detect changes in myelin network architecture during the process of normal aging. This approach opens up new possibilities compared to traditional Tractometry. Methods In a group of 85 healthy controls aged between 18 and 68 years, we estimated myelin-weighted connectomes using Tractometry and MySD, and compared their modulation with age by means of three well-known global network metrics. Results Following the literature, our results show that myelin development continues until brain maturation (40 years old), after which degeneration begins. In particular, mean connectivity strength and efficiency show an increasing trend up to 40 years, after which the process reverses. Both Tractometry and MySD are sensitive to these changes, but MySD turned out to be more accurate. Conclusion After regressing the known predictors, MySD results in lower residual error, indicating that MySD provides more accurate estimates of myelin-weighted connectivity than Tractometry.

Background and Objectives Progression independent of relapse activity (PIRA) is a crucial determinant of overall disability accumulation in multiple sclerosis (MS). Accelerated brain atrophy has been shown in patients experiencing PIRA. In this study, we assessed the relation between PIRA and neurodegenerative processes reflected by (1) longitudinal spinal cord atrophy and (2) brain paramagnetic rim lesions (PRLs). Besides, the same relationship was investigated in progressive MS (PMS). Last, we explored the value of cross-sectional brain and spinal cord volumetric measurements in predicting PIRA. Methods From an ongoing multicentric cohort study, we selected patients with MS with (1) availability of a susceptibility-based MRI scan and (2) regular clinical and conventional MRI follow-up in the 4 years before the susceptibility-based MRI. Comparisons in spinal cord atrophy rates (explored with linear mixed-effect models) and PRL count (explored with negative binomial regression models) were performed between: (1) relapsing-remitting (RRMS) and PMS phenotypes and (2) patients experiencing PIRA and patients without confirmed disability accumulation (CDA) during follow-up (both considering the entire cohort and the subgroup of patients with RRMS). Associations between baseline MRI volumetric measurements and time to PIRA were explored with multivariable Cox regression analyses. Results In total, 445 patients with MS (64.9% female; mean [SD] age at baseline 45.0 [11.4] years; 11.2% with PMS) were enrolled. Compared with patients with RRMS, those with PMS had accelerated cervical cord atrophy (mean difference in annual percentage volume change [MD-APC] −1.41; p = 0.004) and higher PRL load (incidence rate ratio [IRR] 1.93; p = 0.005). Increased spinal cord atrophy (MD-APC −1.39; p = 0.0008) and PRL burden (IRR 1.95; p = 0.0008) were measured in patients with PIRA compared with patients without CDA; such differences were also confirmed when restricting the analysis to patients with RRMS. Baseline volumetric measurements of the cervical cord, whole brain, and cerebral cortex significantly predicted time to PIRA (all p ≤ 0.002). Discussion Our results show that PIRA is associated with both increased spinal cord atrophy and PRL burden, and this association is evident also in patients with RRMS. These findings further point to the need to develop targeted treatment strategies for PIRA to prevent irreversible neuroaxonal loss and optimize long-term outcomes of patients with MS.

Key Points Question Can multiple sclerosis (MS) be differentiated from a wide range of non-MS conditions showing brain white matter lesions using solely imaging biomarkers for cortical lesions (CLs) and central vein sign (CVS)? Findings In this cross-sectional study including 1051 participants, the presence of CLs had high specificity and low sensitivity, while application of the 40% CVS rule resulted in high specificity and moderate sensitivity for MS diagnosis. CVS and CLs outperformed the contribution of infratentorial, periventricular, and juxtacortical lesions in supporting the diagnosis of MS. Meaning The findings indicate that CVS and CLs may be valuable tools to increase the accuracy of MS diagnosis.

Importance Mechanisms contributing to disability accumulation in multiple sclerosis (MS) are poorly understood. Blood neurofilament light chain (NfL) level, a marker of neuroaxonal injury, correlates robustly with disease activity in people with MS (MS); however, data on the association between NfL level and disability accumulation have been conflicting. Objective To determine whether and when NfL levels are elevated in the context of confirmed disability worsening (CDW). Design, Setting, and Participants This study included 2 observational cohorts: results from the Expression, Proteomics, Imaging, Clinical (EPIC) study at the University of California San Francisco (since 2004) were confirmed in the Swiss Multiple Sclerosis Cohort (SMSC), a multicenter study in 8 centers since 2012. Data were extracted from EPIC in April 2022 (sampling July 1, 2004, to December 20, 2016) and SMSC in December 2022 (sampling June 6, 2012, to September 2, 2021). The study included 2 observational cohorts in tertiary MS centers. All participants of both cohorts with available NfL results were included in the study, and no eligible participants were excluded or declined to participate. Exposure Association between NfL z scores and CDW. Main Outcome Measures CDW was defined as Expanded Disability Status Scale (EDSS) worsening that was confirmed after 6 or more months and classified into CDW associated with clinical relapses (CDW-R) or independent of clinical relapses (CDW-NR). Visits were classified in relation to the disability worsening events into CDW(-2) for 2 visits preceding event, CDW(-1) for directly preceding event, CDW(event) for first diagnosis of EDSS increase, and the confirmation visit. Mixed linear and Cox regression models were used to evaluate NfL dynamics and to assess the association of NfL with future CDW, respectively. Results A total of 3906 EPIC visits (609 participants; median [IQR] age, 42.0 [35.0-50.0] years; 424 female [69.6%]) and 8901 SMSC visits (1290 participants; median [IQR] age, 41.2 [32.5-49.9] years; 850 female [65.9%]) were included. In CDW-R (EPIC, 36 events; SMSC, 93 events), NfL z scores were 0.71 (95% CI, 0.35-1.07; P < .001) units higher at CDW-R(-1) in EPIC and 0.32 (95% CI, 0.14-0.49; P < .001) in SMSC compared with stable MS samples. NfL elevation could be detected preceding CDW-NR (EPIC, 191 events; SMSC, 342 events) at CDW-NR(-2) (EPIC: 0.23; 95% CI, 0.01-0.45; P = .04; SMSC: 0.28; 95% CI, 0.18-0.37; P < .001) and at CDW-NR(-1) (EPIC: 0.27; 95% CI, 0.11-0.44; P < .001; SMSC: 0.09; 95% CI, 0-0.18; P = .06). Those findings were replicated in the subgroup with relapsing-remitting MS. Time-to-event analysis confirmed the association between NfL levels and future CDW-R within approximately 1 year and CDW-NR (in approximately 1-2 years). Conclusions and Relevance This cohort study documents the occurrence of NfL elevation in advance of clinical worsening and may hint to a potential window of ongoing dynamic central nervous system pathology that precedes the diagnosis of CDW.

This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion‐relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power‐law approach and a T2‐based pore size estimation technique.

Purpose To characterize the orientation dependence of magnetization transfer (MT) measures in white matter (WM) and propose a first correction method for such measures. Methods A characterization method was developed using the fiber orientation obtained from diffusion MRI (dMRI) with diffusion tensor imaging (DTI) and constrained spherical deconvolution (CSD). This allowed for characterization of the orientation dependence of measures in all of WM, regardless of the number of fiber orientation in a voxel. Furthermore, a first correction method was proposed from the results of characterization, aiming at removing said orientation dependence. Both methods were tested on a 20-subject dataset and effects on tractometry results were also evaluated. Results Previous results for single-fiber voxels were reproduced and a novel characterization was produced in voxels of crossing fibers, which seems to follow trends consistent with single-fiber results. Unwanted effects of the orientation dependence on MT measures were highlighted, for which the correction method was able to produce improved results. Conclusion Encouraging results of corrected MT measures showed the importance of such correction, opening the door for future research on the topic.