Multimodal magnetic resonance longitudinal study on the deep gray matter in multiple sclerosis patients with teriflunomide.

Disease-modifying therapies (DMTs) are used in an increasing number of patients with multiple sclerosis (MS). However, whether DMTs have intrinsic effects on deep gray matter (DGM) microstructure and atrophy is still poorly understood. In this study, we described the quantitative susceptibility values (QSV) and diffusion kurtosis imaging (DKI) metrics of DGM in relapsing-remitting MS (RRMS) patients and their association with cognitive deficits. We recruited 62 patients with RRMS receiving DMTs and 30 patients with RRMS not receiving DMTs underwent MRI on a 3T scanner. Fractional anisotropy (FA), kurtosis fractional anisotropy (KFA), mean diffusivity (MD), mean kurtosis (MK), QSV and volumes of bilateral caudate nucleus (CAU), amygdala (AMYG), putamen (PUT), hippocampus (Hipp), globus pallidus (GP) and thalamus (THA) were measured. Correlation analysis was performed between those image indexes with longitudinal significant changes and clinical neurological scores, including Expanded Disability Status Scale (EDSS), Digit Span Testand (DST), Symbol Digit Modalities Test (SDMT), Mini-Mental State Examination (MMSE) and Montreal Cognitive Assessment (MoCA). Significant longitudinal increases in FA, KFA and MK values were found in both groups in bilateral CAU, AMYG, PUT, Hipp, GP and THA (all p < 0.005). MD values of the right of CAU in the two groups were significant longitudinal increase (p = 0.009, p = 0.047); MD values of the right of GP (p = 0.042), the left of THA (p = 0.003), the right of THA (p = 0.001) in treated MS were significant longitudinal decrease; There were no significant longitudinal changes between treated and untreated groups in normalized deep gray matter volume. For QSV, longitudinal increase in the right of PUT (p = 0.022) in the treated MS group and in the left of Hipp (p = 0.045) in the untreated MS group. The QSV and DKI measures were highly correlated with cognitive and disability tests. The treated RRMS patients showed different longitudinal changes of MD value and QSV with untreated in several DGM regions, and these differences were correlated with cognitive and microstructural integrity.

Alterations of brain structural and functional connectivity networks topology and decoupling in pediatric myelin oligodendrocyte glycoprotein antibody-associated disease.

OBJECTIVE: To investigate the alteration in structural and functional connectivity networks (SCN and FCN) as well as their coupling in pediatric myelin oligodendrocyte glycoprotein antibody-associated disease (MOGAD), and determine if these properties could serve as potential biomarkers for the disease. MATERIALS AND METHODS: In total of 32 children with MOGAD and 30 age- and sex-matched healthy controls (HC) were employed to construct the SCN and FCN, respectively. The graph-theoretical analyses of the global properties, node properties of the 90 brain nodes, and the structural-functional connectivity (SC-FC) coupling of the two networks were performed. The graph-theoretical properties that exhibited significant differences were analyzed using partial correlation analysis in conjunction with the clinical scales, including the expanded disability status scale (EDSS), modified Rankin scale (mRS), and pediatric cerebral performance category (PCPC) of the MOGAD group. Subsequently, a machine learning model was developed to discriminate between MOGAD and the HC group, aiming to explore the potential of these properties as biomarkers. RESULTS: The SCN of the MOGAD group exhibited aberrant global properties, including an increased characteristic path length (Lp) and a decreased global efficiency (Eg), along with reduced nodal properties such as degree centrality (Dc), nodal efficiency (Ne), and local efficiency in multiple nodes. The FCN of the MOGAD group only exhibited decreased Dc, Ne, and betweenness centrality in two nodes of nodal properties. Besides, MOGAD showed a significant decrease in SC-FC coupling compared to the HC group. The analysis of partial correlation revealed significant correlations between several properties and the scales of EDSS and mRS in the MOGAD group. The machine learning method was used to extract six features and establish the model, achieving a classification accuracy of 82.3% for MOGAD. CONCLUSIONS: Pediatric MOGAD showed a more pronounced impairment in the SCN along with decoupling of SC-FC. Both partial correlation analysis and discriminant modeling suggest that alterations in brain network properties have the potential as biomarkers for assessing brain damage in MOGAD.

Cross-sectional and longitudinal evaluation of white matter microstructure damage and cognitive correlations by automated fibre quantification in relapsing-remitting multiple sclerosis patients.

The purpose of this study was to characterize whole-brain white matter (WM) fibre tracts by automated fibre quantification (AFQ), capture subtle changes cross-sectionally and longitudinally in relapsing-remitting multiple sclerosis (RRMS) patients and explore correlations between these changes and cognitive performance A total of 114 RRMS patients and 71 healthy controls (HCs) were enrolled and follow-up investigations were conducted on 46 RRMS patients. Fractional anisotropy (FA), mean diffusion (MD), axial diffusivity (AD), and radial diffusivity (RD) at each node along the 20 WM fibre tracts identified by AFQ were investigated cross-sectionally and longitudinally in entire and pointwise manners. Partial correlation analyses were performed between the abnormal metrics and cognitive performance. At baseline, compared with HCs, patients with RRMS showed a widespread decrease in FA and increases in MD, AD, and RD among tracts. In the pointwise comparisons, more detailed abnormalities were localized to specific positions. At follow-up, although there was no significant difference in the entire WM fibre tract, there was a reduction in FA in the posterior portion of the right superior longitudinal fasciculus (R_SLF) and elevations in MD and AD in the anterior and posterior portions of the right arcuate fasciculus (R_AF) in the pointwise analysis. Furthermore, the altered metrics were widely correlated with cognitive performance in RRMS patients. RRMS patients exhibited widespread WM microstructure alterations at baseline and alterations in certain regions at follow-up, and the altered metrics were widely correlated with cognitive performance in RRMS patients, which will enhance our understanding of WM microstructure damage and its cognitive correlation in RRMS patients.

Interpretable and Intuitive Machine Learning Approaches for Predicting Disability Progression in Relapsing-Remitting Multiple Sclerosis Based on Clinical and Gray Matter Atrophy Indicators.

RATIONALE AND OBJECTIVES: To investigate whether clinical and gray matter (GM) atrophy indicators can predict disability in relapsing-remitting multiple sclerosis (RRMS) and to enhance the interpretability and intuitiveness of a predictive machine learning model. MATERIALS AND METHODS: 145 and 50 RRMS patients with structural MRI and at least 1-year follow-up Expanded Disability Status Scale (EDSS) results were retrospectively enrolled and placed in the discovery and external test cohorts, respectively. Six clinical and radiomics feature-based machine learning classifiers were trained and tested to predict disability progression in the discovery cohort and validated in the external test set. Partial dependence plot (PDP) analysis and a Shiny web application were conducted to enhance the interpretability and intuitiveness. RESULTS: In the discovery cohort, 98 patients had disability stability, and 47 patients were classified as having disability progression. In the external test set, 35 patients were disability stable, and 15 patients had disability progression. Models trained with both clinical and radiomics features (area under the curve (AUC), 0.725-0.950) outperformed those trained with clinical (AUC, 0.600-0.740) or radiomics features only (AUC, 0.615-0.945). Among clinical+ radiomics feature models, the logistic regression (LR) classifier-based model performed best, with an AUC of 0.950. Only the radiomics feature-only models were applied in the external test set due to the data collection problem and showed fair performance, with AUCs ranging from 0.617 to 0.753. PDP analysis showed that female patients and those with lower volume, surface area, and symbol digit modalities test (SDMT) scores; greater mean curvature and age; and no disease modifying therapy (DMT) had increased probabilities of disease progression. Finally, a Shiny web application (https://lauralin1104.shinyapps.io/LRshiny/) was developed to calculate the risk of disability progression. CONCLUSION: Interpretable and intuitive machine learning approaches based on clinical and GM atrophy indicators can help physicians predict disability progression in RRMS patients for clinical decision-making and patient management.

Aberrant white matter microstructure detected by automatic fiber quantification in pediatric myelin oligodendrocyte glycoprotein antibody-associated disease.

BACKGROUND AND OBJECTIVES: Myelin oligodendrocyte glycoprotein antibody-associated diseases (MOGAD) is an idiopathic inflammatory demyelinating disorder in children, for which the precise damage patterns of the white matter (WM) fibers remain unclear. Herein, we utilized diffusion tensor imaging (DTI)-based automated fiber quantification (AFQ) to identify patterns of fiber damage and to investigate the clinical significance of MOGAD-affected fiber tracts. METHODS: A total of 28 children with MOGAD and 31 healthy controls were included in this study. The AFQ approach was employed to track WM fiber with 100 equidistant nodes defined along each tract for statistical analysis of DTI metrics in both the entire and nodal manner. The feature selection method was used to further screen significantly aberrant DTI metrics of the affected fiber tracts or segments for eight common machine learning (ML) to evaluate their potential in identifying MOGAD. These metrics were then correlated with clinical scales to assess their potential as imaging biomarkers. RESULTS: In the entire manner, significantly reduced fractional anisotropy (FA) was shown in the left anterior thalamic radiation, arcuate fasciculus, and the posterior and anterior forceps of corpus callosum in MOGAD (all p < 0.05). In the nodal manner, significant DTI metrics alterations were widely observed across 37 segments in 10 fiber tracts (all p < 0.05), mainly characterized by decreased FA and increased radial diffusivity (RD). Among them, 14 DTI metrics in seven fiber tracts were selected as important features to establish ML models, and satisfactory discrimination of MOGAD was obtained in all models (all AUC > 0.85), with the best performance in the logistic regression model (AUC = 0.952). For those features, the FA of left cingulum cingulate and the RD of right inferior frontal-occipital fasciculus were negatively and positively correlated with the expanded disability status scale (r = -0.54, p = 0.014; r = 0.43, p = 0.03), respectively. CONCLUSION: Pediatric MOGAD exhibits extensive WM fiber tract aberration detected by AFQ. Certain fiber tracts exhibit specific patterns of DTI metrics that hold promising potential as biomarkers.

Radiomics derived from T2-FLAIR: the value of 2- and 3-classification tasks for different lesions in multiple sclerosis.

Background: White matter (WM) lesions can be classified into contrast enhancement lesions (CELs), iron rim lesions (IRLs), and non-iron rim lesions (NIRLs) based on different pathological mechanism in relapsing-remitting multiple sclerosis (RRMS). The application of radiomics established by T2-FLAIR to classify WM lesions in RRMS is limited, especially for 3-class classification among CELs, IRLs, and NIRLs. Methods: A total of 875 WM lesions (92 CELs, 367 IRLs, 416 NIRLs) were included in this study. The 2-class classification was only performed between IRLs and NIRLs. For the 2- and 3-class classification tasks, all the lesions were randomly divided into training and testing sets with a ratio of 8:2. We used least absolute shrinkage and selection operator (LASSO), reliefF algorithm, and mutual information (MI) for feature selection, then eXtreme gradient boosting (XGBoost), random forest (RF), and support vector machine (SVM) were used to establish discrimination models. Finally, the area under the curve (AUC), accuracy, sensitivity, specificity, and precision were used to evaluate the performance of the models. Results: For the 2-class classification model, LASSO classifier with RF model showed the best discrimination performance with the AUC of 0.893 (95% CI: 0.838-0.942), accuracy of 0.813, sensitivity of 0.833, specificity of 0.781, and precision of 0.851. However, the 3-class classification model of LASSO with XGBoost displayed the highest performance with the AUC of 0.920 (95% CI: 0.887-0.950), accuracy of 0.796, sensitivity of 0.839, specificity of 0.881, and precision of 0.846. Conclusions: Radiomics models based on T2-FLAIR images have the potential for discriminating among CELs, IRLs, and NIRLs in RRMS.

Radiomics models based on cortical damages for identification of multiple sclerosis with cognitive impairment.

BACKGROUND: Cognitive impairment (CI) is a common symptom in multiple sclerosis (MS) patients. Cortical damages can be closely associated with cognitive network dysfunction and clinically significant CI in MS. So, in this study, We aimed to develop a radiomics model to efficiently identify the MS patients with CI based on clinical data and cortical damages. METHODS: One hundred and eighteen patients with MS were divided into CI and normal cognitive (NC) cohorts (62/56) as defined by the Montreal Cognitive Assessment (MoCA). All participants were randomly divided into train and test sets with a ratio of 7:3. The radiomic features were selected by using the least absolute shrinkage and selection operator (LASSO) method. The discrimination models were built with the support vector machines (SVM) by the clinical data, radiomic features, and merge data, respectively. And the patients were further divided according to each cognitive domain including memory, visuospatial, language, attention and executive, and each domain model was applied by the most suitable classifier. RESULTS: A total of 2298 features were extracted, of which 36 were finally selected. The merge model showed the greatest performance with the area under the curve (AUC) of 0.86 (95 % confidence interval: 0.81-0.91), accuracy (ACC) of 0.78, sensitivity of 0.79 and specificity of 0.77 in test cohort. However, although the visuospatial domain model showed the highest AUC of 0.71 (95 % confidence interval: 0.61-0.81) among five domain models, other domain models did not meet satisfactory results with a relatively low AUC, ACC, sensitivity and specificity. CONCLUSIONS: The radiomics model based on clinical data and cortical damages had a great potential to identify the MS patients with CI for clinical cognitive assessment.

Multimodal Investigation of Deep Gray Matter Nucleus in Patients with Multiple Sclerosis and Their Clinical Correlations: A Multivariate Pattern Analysis Study.

Deep gray matter (DGM) nucleus are involved in patients with multiple sclerosis (MS) and are strongly associated with clinical symptoms. We used machine learning approach to further explore microstructural alterations in DGM of MS patients. One hundred and fifteen MS patients and seventy-one healthy controls (HC) underwent brain MRI. The fractional anisotropy (FA), mean diffusivity (MD), quantitative susceptibility value (QSV) and volumes of the caudate nucleus (CN), putamen (PT), globus pallidus (GP), and thalamus (TH) were measured. Multivariate pattern analysis, based on a machine-learning algorithm, was applied to investigate the most damaged regions. Partial correlation analysis was used to investigate the correlation between MRI quantitative metrics and clinical neurological scores. The area under the curve of FA-based classification model was 0.83, while they were 0.93 for MD and 0.81 for QSV. The Montreal cognitive assessment scores were correlated with the volume of the DGM and the expanded disability status scale scores were correlated with the MD of the GP and PT. The study results indicated that MS patients had involvement of DGM with the CN being the most affected. The atrophy of DGM in MS patients mainly affected cognitive function and the microstructural damage of DGM was mainly correlated with clinical disability.

Increased cortical lesion load contributed to pathological changes beyond focal lesion in cortical gray matter of multiple sclerosis: a diffusion kurtosis imaging analysis.

Biomarkers specific to cortical gray matter (cGM) pathological changes of multiple sclerosis (MS) are desperately needed to better understand the disease progression. The cGM damage occurs in cortical lesion (CL) and normal-appearing cGM (NAcGM) areas. While the association between CL load and cGM damage has been reported, little is known about how different CL types, i.e. intracortical lesion (ICL) and leukocortical lesion (LCL) would be associated with cGM damage. In our study, relapsing-remitting MS patients and healthy controls were divided into 4 groups according to CL load level. NAcGM diffusion kurtosis imaging (DKI)/diffusion tensor imaging (DTI) values and cGM volume (cGMV) were used to characterize the pathological changes in cGM. Univariate general linear model was used for group comparisons and stepwise regression analysis was used to assess the effects of ICL volume and LCL volume on NAcGM damage. We found peak values in DKI/DTI values, cGMV and neuropsychological scores in high CL load group. Kurtosis fractional anisotropy (KFA) was the most sensitive in characterizing NAcGM damage, and LCL volume related more to NAcGM damage. Our findings suggested KFA could become a surrogate biomarker to cGM damage, and LCL might be the main factor in whole brain NAcGM damage.

Enlarged choroid plexus related to iron rim lesions and deep gray matter atrophy in relapsing-remitting multiple sclerosis.

BACKGROUND: Choroid plexus (CP) is considered to be linked to inflammation of multiple sclerosis (MS), but its connection with markers of inflammation in vivo in MS is unclear, the markers such as lesions load and brain atrophy, particularly the white matter lesions (WMLs) edge surrounded by an iron rim, termed as iron rim lesions (IRLs). PURPOSE: To investigate the association between CP volume and brain lesions load, especially IRLs load and atrophy in MS, and its relationship with clinical characteristics. METHODS: 3.0 T brain MRI images were acquired from 99 relapsing-remitting MS (RRMS) and 60 healthy controls (HCs) to obtain the volumes of CP, whole brain and lesions. Volumes were expressed as a ratio of intracranial volume. Expanded Disability Status Scale (EDSS), Montreal Cognitive Assessment (MoCA) and Symbol Digit Modalities Test (SDMT) were used to assess the severity of disability and cognitive function. Student's t-test and Multivariable regression analyses were performed to evaluate the difference of CP volumes between RRMS and HC and the association between CP volume and lesions load, brain volumes and clinical scale scores in RRMS. RESULTS: CP volume was 30% larger in patients with RRMS than HCs (p < 0.001) and was 20% larger in patients with IRLs than those without IRLs (p = 0.007). Moreover, the larger CP volume was related to greater WMLs volume in the whole RRMS (r = 0.46, p < 0.001). Further analysis in patients with IRLs showed a positive correlation between CP volume and WMLs volume (r = 0.45, p = 0.003), and IRLs volume (r = 0.51, p < 0.001). Meanwhile, enlarged CP was related to lower volumes in the whole brain (r = -0.30, p = 0.006), deep gray matter (r = -0.51, p < 0.001) and most regional deep gray matter nuclei (except amygdala), but no correlation with cortical lesions or cortex volume (both p > 0.05). In addition, CP volume was significantly higher in patients with cognitive impairment than those with cognitive preservation by MoCA scores (p = 0.011); the larger CP volume was associated with higher EDSS scores (r = 0.25, p = 0.014) and lower SDMT Z scores in RRMS (r = -0.26, p = 0.014). CONCLUSION: The enlargement of CP in RRMS had close correlations with inflammatory lesions, especially IRLs and deep gray matter atrophy, but not the cortex. Meanwhile, the larger CP volume was associated with higher disability and lower cognitive scores. CP volume may be a surrogate imaging marker for MS disease activity.

Microstructural alterations in different types of lesions and their perilesional white matter in relapsing-remitting multiple sclerosis based on diffusion kurtosis imaging.

BACKGROUND AND OBJECTIVES: In multiple sclerosis (MS), contrast enhancement lesions and chronic active lesions have been demonstrated to have different degrees of inflammation. Accordingly, they exist different degrees of tissue damage, one is short and acute, and another is slow and longstanding. This study aimed to explore whether diffusion parameters can differentiate different types of lesions, and investigate the microstructural damage between different types of MS lesions by using diffusion magnetic resonance imaging (dMRI) and its correlation with clinical biomarkers of disability and cognitive states. METHODS: We retrospectively identified 77 contrast enhancement lesions (CELs), 384 iron rim lesions (IRLs), 393 non-iron rim lesions (NIRLs), their corresponding perilesional white matter (PLWM), and 68 normal-appearing white matter (NAWM) from 68 relapsing-remitting MS (RRMS). Additionally, 44 white matter in healthy controls (WM in HCs) were also enrolled in this study. The DTI and DKI parameters were measured in the above white matter, including kurtosis fractional anisotropy (KFA), fractional anisotropy (FA), mean kurtosis (MK), and mean diffusivity (MD). All the patients were assessed with the Digital Span Test (DST), the Symbol Digit Modalities Test (SDMT), the Mini-Mental State Examination (MMSE), the Montreal Cognitive Assessment (MoCA), and the Expanded Disability Status Scale (EDSS). RESULTS: The lowest KFA, FA, MK values and the highest MD values were found in CELs, followed by IRLs, NIRLs, NAWM, and WM in HCs. In KFA and FA values, there were significant differences between each type of lesion, as well as each type of PLWM (P < 0.05). The MK values of CELs and IRLs were significantly lower than NIRLs, but inversely for MD (P < 0.05). There were no differences between CELs and IRLs for MK (P = 1) and MD (P = 0.261). The results of MK and MD values in CELs-PLWM and IRLs-PLWM were similar to the CELs and IRLs. There were no significant differences between NAWM and WM in HCs in all the enrolled diffusion parameters (P >0.05) and the FA values between NIRLs-PLWM and NAWM or between NIRLs-PLWM and WM in HCs were no significant differences (P >0.05). The KFA and MD values in IRLs-PLWM (r =0.443, P =0.021; r =-0.518, P =0.006) were correlated with the DST scores and the KFA of CELs-PLWM (r =0.396, P =0.041) was correlated with SDMT scores. CONCLUSION: Our findings demonstrate that the KFA values have the potential to distinguish different types of MS white matter tissues. Furthermore, the diffusion parameters can reflect the microstructure abnormalities in different MS lesions and might help us better understand the pathological mechanism and lesion evolution.

Alterations in White Matter Fiber Tracts Characterized by Automated Fiber-Tract Quantification and Their Correlations With Cognitive Impairment in Neuromyelitis Optica Spectrum Disorder Patients.

Objectives: To investigate whether patients with neuromyelitis optica spectrum disorder (NMOSD) have tract-specific alterations in the white matter (WM) and the correlations between the alterations and cognitive impairment. Materials and Methods: In total, 40 patients with NMOSD and 20 healthy controls (HCs) who underwent diffusion tensor imaging (DTI) scan and neuropsychological scale assessments were enrolled. Automated fiber-tract quantification (AFQ) was applied to identify and quantify 100 equally spaced nodes of 18 specific WM fiber tracts for each participant. Then the group comparisons in DTI metrics and correlations between different DTI metrics and neuropsychological scales were performed. Results: Regardless of the entire or pointwise level in WM fiber tracts, patients with NMOSD exhibited a decreased fractional anisotropy (FA) in the left inferior fronto-occipital fasciculus (L_IFOF) and widespread increased mean diffusion (MD), axial diffusivity (AD), and radial diffusivity (RD), especially for the thalamic radiation (TR), corticospinal tract (CST), IFOF, inferior longitudinal fasciculus (ILF), superior longitudinal fasciculus (SLF) [p < 0.05, false discovery rate (FDR) correction], and the pointwise analyses performed more sensitive. Furthermore, the negative correlations among MD, AD, RD, and symbol digit modalities test (SDMT) scores in the left TR (L_TR) were found in NMOSD. Conclusion: Patients with NMOSD exhibited the specific nodes of WM fiber tract damage, which can enhance our understanding of WM microstructural abnormalities in NMOSD. In addition, the altered DTI metrics were correlated with cognitive impairment, which can be used as imaging markers for the early identification of NMOSD cognitive impairment.