- 1Department of Electronic Engineering, Beijing National Research Center for Information Science and Technology (BNRist), Tsinghua University, Beijing, China
- 2Department of Radiology, Aerospace Center Hospital, Beijing, China
- 3Department of Neurology, Dongfang Hospital, Beijing University of Chinese Medicine, Beijing, China
In order to deeply understand the specific patterns of volume, microstructure, and functional changes in Multiple System Atrophy patients with cerebellar ataxia syndrome (MSA-c), we perform the current study by simultaneously applying structural (T1-weighted imaging), Diffusion tensor imaging (DTI), functional (BOLD fMRI) and extended Network-Based Statistics (extended-NBS) analysis. Twenty-nine MSA-c type patients and twenty-seven healthy controls (HCs) were involved in this study. First, we analyzed the whole brain changes of volume, microstructure, and functional connectivity (FC) in MSA-c patients. Then, we explored the correlations between significant multimodal MRI features and the total Unified Multiple System Atrophy Rating Scale (UMSARS) scores. Finally, we searched for sensitive imaging biomarkers for the diagnosis of MSA-c using support vector machine (SVM) classifier. Results showed significant grey matter atrophy in cerebellum and white matter microstructural abnormalities in cerebellum, left fusiform gyrus, right precentral gyrus and lingual gyrus. Extended-NBS analysis found two significant different connected components, featuring altered functional connectivity related to left and right cerebellar sub-regions, respectively. Moreover, the reduced fiber bundle counts at right Cerebellum_3 (Cbe3) and decreased fractional anisotropy (FA) values at bilateral Cbe9 were negatively associated with total UMSARS scores. Finally, the significant features at left Cbe9, Cbe1, and Cbe7b were found to be useful as sensitive biomarkers to differentiate MSA-c from HCs according to the SVM analysis. These findings advanced our understanding of the neural pathophysiological mechanisms of MSA from the perspective of multimodal neuroimaging.
Introduction
Multiple system atrophy (MSA) is a progressive neurodegenerative disorder. Its pathological feature is the deposition of alpha synuclein-positive glial cytoplasmic inclusions (GCIs) in some specific regions, such as cerebellum, striatum and olivopontine structures (Brettschneider et al., 2017; Krismer et al., 2019). MSA patients might be associated with autonomic dysfunction. It is further divided into two sub-types, according to whether the patient shows cerebellar ataxia syndrome (MSA-c) or poor levodopa-responsive parkinsonian syndrome (MSA-p) (Gilman et al., 2008). Alpha synuclein-positive GCI has been shown to be associated with the development of MSA in recent studies. However, the pathological mechanism is not well understood, and it remains to be a big challenge to accurately diagnose MSA.
Magnetic resonance imaging (MRI) has been applied to reveal the structural patterns of MSA patients in previous researches (Matsusue et al., 2009; Deguchi et al., 2015; Krismer and Wenning, 2017; Chelban et al., 2019; Dash et al., 2019). As a result, atrophy of the cerebellum, pons, middle cerebellar peduncle (MCP) and putamen was found in MSA patients, reflecting the pathological changes of the disease (Matsusue et al., 2009; Deguchi et al., 2015). Furthermore, some MRI studies reported the “putaminal slit” sign which was associated with the accumulation of iron and gliosis and the “hot-cross bun” sign which was associated with the selective depletion of myelinated transverse pontocerebellar fibers and pontine neurons in different types of MSA patients (Krismer and Wenning, 2017; Chelban et al., 2019; Dash et al., 2019). In order to improve the diagnostic accuracy, more detailed grey matter abnormality has been found in research settings, among which the cerebellar volume atrophy is widely observed (Yang et al., 2019; Lin et al., 2020). However, some MSA patients might have no abnormal structure changes on MRI (Chelban et al., 2019). It is therefore necessary to search for more sensitive biomarkers from other MRI modalities to accurately diagnose the disease.
Diffusion tensor imaging (DTI) is a noninvasive neuroimaging technique that can provide quantitative information of the diffusion characteristics of the water molecules in the brain. Using DTI, white matter changes including a significant decrease in fractional anisotropy (FA) and increase in mean diffusivity (MD) was observed in corticospinal tracts, transverse pontocerebellar fibers, pons, putamen, middle cerebellar peduncles (MCP) and cerebellum in MSA patients (Loh et al., 2011; Rulseh et al., 2016; Dash et al., 2019; Beliveau et al., 2021).
Besides the volumetric and microstructural changes of MSA, resting state functional MRI (fMRI) has revealed abnormal functional activity and connectivity patterns related to MSA. Previous studies reported the disrupted default mode network (DMN), striatal-thalamo-cortical network, cerebello-thalamo-cortical network, sensorimotor network and visual associated network in MSA patients (You et al., 2011; Wang et al., 2017; Yao et al., 2017; Ren et al., 2019; Zheng et al., 2019). Furthermore, researchers detected altered network topology and graph theory attributes by resting state graph theoretical analysis (Morisi et al., 2018; Sako et al., 2019). Previous studies provided evidence for the hypothesis that MSA may be caused by the disconnection syndrome, indicating that the accumulated alpha-synuclein may destroy several specific networks, and result in the movement disorder and other clinical symptoms (Rosskopf et al., 2018).
Most of previous studies focused on structural changes of MSA. Only a few studies reported new functional findings, which might improve diagnostic accuracy and provide novel disease progression markers. However, structural or functional analysis alone cannot provide all the necessary information to accurately diagnose MSA. Thus, a multimodal approach combining these innovative techniques seems likely to offer the best biomarker for clinical diagnosis of MSA. We speculate that a combination of structural grey matter and white matter features, DTI microstructure metrics and functional connectivity calculated from resting-state fMRI might improve diagnostic accuracy of MSA. In order to explore the patterns of brain volume, microstructural and functional changes of MSA, we simultaneously applied structural (T1-weighted imaging), DTI, functional (BOLD fMRI) and extended Network-Based Statistics (extended-NBS) analysis in this study. First of all, we analyzed the whole brain changes of volume, microstructure, and functional connectivity in MSA-c patients. Then, the correlations between multimodal MRI features and clinical measurements were explored. Finally, classification is performed using support vector machine (SVM), in order to search for sensitive imaging biomarkers for the diagnosis of MSA.
Materials and Methods
Participants
The participant collection is similar as our previous study (Zheng et al., 2020). Thirty-four right-handed MSA-c patients and twenty-nine right-handed healthy controls (HCs) matched in age and gender were recruited. All the participants were from Dongfang Hospital of Beijing University of Chinese Medicine. The recruitment criteria for MSA-c patients in this study were based on the international diagnostic criteria, which was defined by American Academy of Neurology and American Autonomic Society (Gilman et al., 2008). The study included patients diagnosed as probable MSA and without hemorrhage, infarction, tumors, trauma, or severe white matter hyperintensity. Each patient was assessed by the Unified Multiple System Atrophy Rating Scale (UMSARS), Montreal Cognitive Assessment (MoCA), Hamilton Anxiety Scale (HAMA) and Hamilton Depression Scale (HAMD). The healthy control group was included according to the following criteria: (1) no neurological or psychiatric disorders; (2) no significant cognitive decline; (3) no visual, auditory or other neurological dysfunction; (4) had not received deep brain stimulation or surgical treatment; (5) MRI scan showed no motor impairment, vascular brain injury, brain tumor, and/or significant cortical and/or subcortical atrophy. In accordance with the Declaration of Helsinki, each of the subjects gave a written informed consent. The study was approved by the Medical Research Ethical Committee of Dongfang Hospital of Beijing University of Chinese Medicine.
Among the recruited thirty-four patients, one was later identified as mis-diagnosed (patient 23) and excluded. Patient 32 did not finish clinical assessment and was discarded. During MRI scanning, two patients (patient 26 and 34) and two controls (subject 25 and 28) did not complete structural scanning and were not included in further analysis. In correlation analysis, one more subject (patient 21) was identified as an outlier and was excluded. As a result, the experiment was conducted using data collected from twenty-nine MSA-c patients and twenty-seven healthy controls. The data exclusion procedure was visualized in Supplementary Figure 1.
Data Acquisition
The image dataset was acquired using a GE 3.0T Discovery 750 scanner. 3D T1-weighted sequence parameters were as follows: repetition time (TR) = 8.2 ms, echo time (TE) = 3.17 ms, inversion time (TI) = 450 ms, flip angle (FA) = 12°, slices = 188, thickness = 1 mm, resolution = 256 × 256 matrix, voxel size = 1 mm × 1 mm× 1 mm. Diffusion tensor imaging was collected using the echo-planar imaging (EPI) sequence with 7 minutes scanning time. The imaging parameters were as follows: TR = 6,000 ms, TE = minimum, FA=12°, thickness = 3 mm, gap = 0 mm, FOV = 25.6 × 25.6 cm2, slices = 50, bandwidth = 250 kHz, resolution = 128 × 128 matrix. 64 diffusion-sensitive gradients in nonlinear directions were adopted, and the b values of the diffusion gradients were 0 s/mm2 and 1,000 s/mm2. The resting-state fMRI scanning lasted for 6 min with the following parameters: TR = 2,000 ms, TE = 30 ms, FA = 90°, FOV = 24 cm × 24 cm, resolution = 64 × 64 matrix, thickness = 3 mm, slices = 36, gap = 1 mm, voxel size = 3.75 mm × 3.75 mm × 3 mm, bandwidth = 2,232 Hz/pixel. During scanning, all participants were instructed to keep their eyes closed, move as little as possible, think of nothing in particular, and stay awake. In total, 6,480 functional images were acquired for each subject, corresponding to 180 volumes.
Data Processing
The method overview is shown in Figure 1. In order to merge multimodal features, data from each modality was first normalized to the standard Montreal Neurological Institute (MNI) 152 space with isotropic spatial resolution. We utilized the Automated Anatomical Labeling (AAL) atlas (Tzourio-Mazoyer et al., 2002) which parceled the brain and cerebellum into 116 areas for further analysis. For volumetric features extracted from T1 and DTI, values at voxels within each brain region volume were averaged as the feature associated with this brain area, and the whole brain feature was converted to a vector with 116 values, each corresponding to one brain region. For functional data, we constructed whole brain networks with 116 nodes by calculating functional connectivity between each pair of nodes. In the following subsections, we describe data processing methods for each modality.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799251/fnagi-14-799251-HTML/image_m/fnagi-14-799251-g001.jpg)
Figure 1. Method overview. Data from each modality was first preprocessed and normalized to the MNI152 standard space, producing a range of raw features. The AAL atlas was utilized to convert the raw features to feature vectors and statistical analysis was performed. We selected significant features for each brain region for further correlation as well as classification analysis. MNI, Montreal Neurological Institute; AAL, Automated Anatomical Labeling.
T1-Weighted Imaging Analysis
T1-weighted images were preprocessed using Statistical Parametric Mapping 12 (SPM12)1 software in MATLAB 2018b (Math Works, Natick, MA, United States). First, the T1-weighted images were segmented into three components, including grey matter (GM), white matter (WM) and cerebrospinal fluid (CSF). Second, we utilized the Diffeomorphic Anatomical Registration through Exponentiated Lie algebra (DARTEL) algorithm to normalize images into the standard MNI 152 space with 1mm isotropic resolution. Finally, the normalized GM and WM components were modulated for nonlinear change and smoothed using a Gaussian kernel of 8 mm full width at half-maximum (FWHM).
We utilized GM and WM components for further analysis and applied the AAL atlas to convert the raw data to feature space. The weighted grey matter volume (GMV) and weighted white matter volume (WMV) for each brain region was extracted by averaging values at each voxel within each region. Since the grey matter and white matter components produced by SPM segmentation are probability maps, the average of values at voxels within each region can be viewed as a weighted volume. As a result, the whole-brain GMV and WMV of each subject was represented by a feature vector.
Diffusion Tensor Imaging Analysis
The preprocessing of DTI data was based on FSL (Smith et al., 2004). First, brain extraction was performed by bet2 and fslmaths. Then we corrected the eddy currents and movements in diffusion data. Finally, the DTI data was normalized to a 2mm-isotropic T1 template in the MNI 152 standard space.
We used an in-house software, Multi-Modal Data Processing System (MMDPS) (Wang, 2018) and the Dipy package (Garyfallidis et al., 2014) for DTI data processing. The Constrained Spherical Deconvolution (CSD) model (Tournier et al., 2007) was firstly fitted with the preprocessed DTI data. Then the probabilistic tracking was performed by following the possible trajectory pathway of neural fibers. Based on the tracking results, we applied the AAL atlas (Tzourio-Mazoyer et al., 2002) and calculated the number of fibers going through each brain region. The tractography feature was represented by a feature vector of number of fiber bundles crossing each region (fiber bundle counts). Besides, we also calculated the fractional anisotropy (FA) and mean diffusivity (MD) with Dipy, and converted them into feature vectors using the AAL atlas by averaging values at voxels within each region.
Resting-State fMRI Analysis
The preprocessing of fMRI data was similar as in our previous research (Zheng et al., 2020). We used DPARSFA (V4.3) (Yan and Zang, 2010) and SPM12 (V6906) to preprocess data. Firstly, we discarded the first 10 time points and performed slice timing correction. Then, head motion was corrected and the data was normalized to the BOLD EPI template in the standard MNI 152 space with 2mm isotropic resolution. The image was resampled to 3-mm isotropic voxels and spatially smoothed by a 4mm full-width half-maximum (FWHM) Gaussian kernel. We further removed the linear trend and nuisance covariates, including head motions, cerebral fluid, white matter and the global signal. Finally, the signal time series were filtered and signals within 0.01-0.08 Hz were kept.
We used an in-house software, Multi-Modal Data Processing System (MMDPS) (Wang, 2018) for further data processing. The brain network nodes were defined using the AAL atlas (Tzourio-Mazoyer et al., 2002) with 116 areas. The BOLD time series was averaged within each brain region and functional connectivity was calculated as the Pearson correlation coefficient between each pair of regions. The functional connectivity is defined as the edge in the network. As a result, we built a brain network with 116 nodes for each subject.
Statistical Analysis
Two-sample t-tests were used to find significant different brain regions for grey matter volume (GMV), white matter volume (WMV), FA, MD as well as fiber bundle count results. The feature values averaged at each region volume were extracted from the two groups and group comparison was performed. All regions across the whole brain were tested and the Benjamini - Hochberg procedure was adopted to control the False Discovery Rate (FDR). The feature at a specific region was declared significant if the p-value is less than 0.001 (FDR corrected). Significant features were identified for metrics extracted from each MRI modality.
In order to evaluate clinical relativity of imaging features, we calculated Pearson correlation coefficients between significant features (p < 0.001, FDR corrected) and UMSARS-total scores of MSA-c patients. Features whose correlation p-value < 0.001 were selected for further discussion. During correlation analysis, we also inspected the distribution of significant features and clinical measurements by drawing the quantile-quantile plot (Q-Q plot) in order to examine for possible outliers in correlation. The statistical analysis, including two-sample t-tests between groups and correlation between features and clinical variables, was performed using the SciPy module, while multiple comparison correction was performed using the statsmodels module.
For functional network analysis, we utilized the extended network-based statistics (extended-NBS) (Zalesky et al., 2010; Ge et al., 2021) to evaluate functional connectivity changes in MSA patients. Firstly, the same connection was selected from the two groups and two-sample t-test was performed. This comparison was repeated for each connection in the network. As a result, each connection was associated with a t-statistic and p-value. A difference network whose edges were p-values obtained from the two-sample t-test was constructed. Then, connections whose p-value is larger than or less than a pre-defined p-value threshold were set to zero or one, respectively, producing a binary difference network. Connected components (CC) in the binary difference network were identified as a set of inter-connected regions. In order to decide whether a CC is significant, permutation test was performed. During each permutation, group label was randomly shuffled and the same procedures were performed, including comparing each connection, constructing binary difference networks, identifying CCs and storing the maximal size of CCs. The size of a connected component is defined as the number of regions within the CC. The permutation was performed 5,000 times, yielding an empirical distribution of the CC size. The empirical p-value of the original CC was calculated as the ratio of the number of permutations whose maximal CC size was larger than the original CC size, against total permutation number. Formally, let CC_sizei denote the maximal CC size at the i-th permutation, CC_size0 denote the original CC size. The empirical p-value is calculated as
Where
If the empirical p-value was less than 0.05, the original CC was declared significant. During experiment, we varied the pre-defined p-value threshold from 0.0001 to 0.001 with a step size of 0.0001. Results obtained at p-value threshold 0.0001 were taken for further discussion.
Classification
We constructed classifiers based on the significant features of each brain region in order to locate regions with biomarker potentials. The significant results produced by statistical analysis were taken as features for classification. In order to incorporate as many features as possible, we used less strict significance criteria compared to statistical analysis. Specifically, each region was associated with a list of significant features obtained from two-sample t-tests of features of each modality (p < 0.05, FDR corrected). The functional connections related to this region was also included if this region is involved in the significant CC when p-value threshold equals 0.001. For example, if region A is significant different in GMV, FA, and MD, then the feature vector of region A would contain [GMV (A), FA (A), MD (A)]. Moreover, if region A lies within the significant CC, the feature vector would further contain functional connections in the significant CC that directly related to region A. Before classification, the feature was standardized by removing the mean and scaling to unit variance.
Where x stands for the feature vector, μ stands for the mean and s stands for the standard deviation.
We utilized the linear Support Vector Machine (SVM) as classifier and performed nested leave-one-out cross validation (LOOCV) for parameter selection and testing. In the outer LOOCV, one sample was left out as test sample and other samples were fed into the inner LOOCV, where another sample was selected as validation sample. The data left was used to train a classifier with different parameters. The model parameter, C, was varied in range [0.001, 0.01, 0.1, 1, 10, 100, 500, 1000, 5000, 10000]. After selecting the best parameter using the validation sample, all samples in the inner LOOCV were utilized to train the SVM with the selected parameter and the model was tested by the left-out sample in the outer LOOCV. The leave-one-out procedure was repeated until all samples have been left out once. We calculated accuracy, sensitivity and specificity for model comparison. Their definitions are as follows.
Where TP stands for true positive, TN stands for true negative, FP stands for false positive and FN stands for false negative.
We also plotted the Receiver Operating Characteristic (ROC) curve and calculated the area under the curve (AUC) value for the classifier trained at each region. The training procedure was repeated for each brain region and the results were sorted by classification accuracy.
Results
Demographics
Twenty-nine MSA-c patients and twenty-seven healthy controls were involved in the current study. Table 1 showed the demographic and clinical information for these subjects. There was no significant difference in age and gender between MSA-c patients and healthy controls (p > 0.05).
Feature Differences and Correlations
Structural atrophy and white matter microstructural alterations in cerebellar and cerebral regions was identified in MSA-c patients. GMV at all subregions of the cerebellum were significantly decreased (p < 0.001, FDR corrected) in MSA-c patients (Figure 2A). WMV obtained from T1 showed significant decrease at bilateral Cerebellum_3 (Cbe3) and bilateral Cerebellum_10 (Cbe10) (Figure 2B). The right lingual gyrus, however, showed increased WMV (Figure 2C). On the other hand, features extracted from DTI reflected white matter alteration related to MSA-c patients. We identified reduced FA, fiber bundle count and increased MD at a range of cerebellar regions (Figures 3A–C). Increased MD in the left fusiform (FFG) and right precentral gyrus (PreCG) was also found in MSA-c patients (Figure 3D). The detailed statistics as well as significant region lists for each feature were reported in Supplementary Tables 1–5.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799251/fnagi-14-799251-HTML/image_m/fnagi-14-799251-g002.jpg)
Figure 2. Significant regions found from T1. The t-statistics of each significant region were visualized for each feature. (A) Decreased GMV in cerebellum; (B) Decreased WMV in cerebellum; (C) Increased WMV in cerebrum area. GMV, gray matter volume; WMV, white matter volume.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799251/fnagi-14-799251-HTML/image_m/fnagi-14-799251-g003.jpg)
Figure 3. Significant regions found from DTI. The t-statistics of each significant region were visualized for each feature. (A) Decreased FA in cerebellum; (B) Increased MD in cerebellum; (C) Reduced fiber bundle counts in cerebellum; (D) Increased MD in cerebrum areas. DTI, Diffusion tensor imaging; FA, fractional anisotropy; MD, mean diffusivity.
We constructed functional brain networks and compared between the two groups using extended-NBS. When varying the p-value threshold from 0.0001 to 0.001, several significant connected components whose empirical p-value was less than 0.05 were identified. The complete results were shown in Supplementary Figure 2. The first component contained increased connectivity alone, mainly consisting connections related to the right cerebellum. The second component featured connections related to the left cerebellum, and contained both increased and decreased connectivity. Another component that appeared at p-value threshold of 0.0002 contained decreased connections related to the right cerebellum, and was further merged into the second component as the p-value threshold increased to 0.0004. The first component was also incorporated into the second component when the p-value threshold increased over 0.0009. Component at p-value threshold 0.0010 was used as features for classification, which contained information provided in all three significant different components (Figure 4C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799251/fnagi-14-799251-HTML/image_m/fnagi-14-799251-g004.jpg)
Figure 4. Extended-NBS results. (A,B) Significant different connected components obtained at p-value threshold of 0.0001 with empirical p-value of 0.0020 and 0.0230, respectively; (C) Significant different connected component obtained at p-value threshold of 0.0010 with empirical p-value less than 0.0001. Red and blue lines stand for increased and decreased connectivity, respectively. Extended-NBS, Extended Network-Based Statistics.
When p-value threshold equals 0.0001, two significant different connected components were identified. The first component contained both increased and decreased connections related to sub-regions of the left cerebellum (Figure 4A). Decreased connections include FC between right middle frontal gyrus (MFG) and left Cerebellum_Crus1 (Cbe1), Cerebellum_2 (Cbe2), and FC between right inferior parietal lobe (IPL) and left Cerebellum_7b (Cbe7b). Increased connections are FC between left Cbe7b and bilateral heschl gyrus (HES), and FC between left Cbe2 and left HES. The second component contained increased connections related to sub-regions of the right cerebellum (Figure 4B). The increased connections are FC between right medial orbitofrontal cortex (mOFC) and right Cerebellum_6 (Cbe6), Cerebellum_8 (Cbe8), and FC between right anterior cingulate cortex (ACC) and right Cbe6.
Among the identified significant features, we further performed correlation analysis with UMSARS-total scores. Patient 21 was identified as an outlier by using Q-Q plot, which was excluded from all statistical and classification analysis. The outlier identification procedures were reported in the Supplementary Material. After excluding the outlier, we identified that the fiber bundle count of right Cerebellum_3 (Cbe3) and FA values at bilateral Cerebellum_9 (Cbe9) had negative significant correlation with clinical measurements (p < 0.001, Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799251/fnagi-14-799251-HTML/image_m/fnagi-14-799251-g005.jpg)
Figure 5. Correlations between the MRI significant features and UMSARS-total scores in MSA-c patients. (A) FA values of left Cbe9; (B) FA values of right Cbe9; (C) fiber bundle counts of right Cbe3. MRI, magnetic resonance imaging; UMSARS, Unified Multiple System Atrophy Rating Scale; MSA-c, cerebellum-type of multiple system atrophy; Cbe, Cerebellum; FA, fractional anisotropy.
Classification Results
For each brain region, we selected significant different features (p < 0.05, FDR corrected) from each modality, and combined functional connectivity related to this region from extended-NBS result at p-value threshold of 0.0010, forming the feature vector for classification. The linear support vector machine with embedded leave-one-out cross validation was performed to evaluate the potential value of diagnosis. The SVM analysis was repeated for each region and results were summarized in Supplementary Table 6. Here we report regions whose classification accuracy is higher than 0.90. By using significant features from all modalities, we found left Cbe9 reached the highest classification accuracy (0.928), with 1.000 sensitivity and 0.852 specificity. We also identified left Cerebellum_Crus1 (Cbe1) and left Cbe7b with 0.911 accuracy when classifying patients without white matter volume features. The classification results of these regions as well as the modalities involved were shown in Table 2. We also plotted the receiver operating characteristics (ROC) curves of these regions in Figure 6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799251/fnagi-14-799251-HTML/image_m/fnagi-14-799251-t002.jpg)
Table 2. Classification results of selected regions. The features were significant results associated with each region.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799251/fnagi-14-799251-HTML/image_m/fnagi-14-799251-g006.jpg)
Figure 6. ROC curve of selected regions when classifying with significant features of each region. (A) left Cerebellum_9; (B) left Cerebellum_Crus1; (C) left Cerebellum_7b. ROC, Receiver Operating Characteristic.
Discussion
Major Findings
By simultaneously applying structural (T1-weighted imaging), DTI, functional (BOLD fMRI) and extended-NBS analysis, we found significant grey matter atrophy in cerebellum and white matter microstructural abnormalities in cerebellum, left FFG and right PreCG and lingual gyrus. Extended-NBS analysis revealed two significant connected components, featuring altered functional connectivity related to left and right cerebellar sub-regions, respectively. Moreover, the reduced fiber bundle counts of right Cbe3 and decreased FA values of bilateral Cbe9 were negatively associated with UMSARS-total scores. Finally, we found that the significant features of left Cbe9, Cbe1 and Cbe7b could be used as sensitive biomarkers to differentiate MSA-c from HCs according to the SVM analysis. The results of this study may help us to understand the neural pathophysiological mechanism of MSA from the perspective of multimodal neuroimaging.
Structural Atrophy, White Matter Microstructural and Functional Abnormalities in MSA-c Patients in the Resting State
Significant grey matter atrophy and white matter microstructural abnormalities were found in most of the cerebellum subregions in MSA-c patients in the current study. It was shown that glial cytoplasmic inclusion body associated oligodendrocyte disease is the main pathological process of MSA (Wenning et al., 2008), and the definitive histopathological manifestation of MSA-c is mainly found in the cerebellum. In a recent study, the researchers found grey matter atrophy and white matter degeneration of cerebellum in MSA-c patients, which was consistent with the current study (Dash et al., 2019). It has been reported that the cerebellum plays an important role in the sensorimotor circuit (Timmann and Daum, 2010). Morphological changes such as structural atrophy may be caused by the primary disease process induced by pathological deposition. Additionally, increased MD in left FFG and right PreCG were found in this study. As a part of the visual association cortex, the fusiform gyrus is mainly responsible for visual association, episodic memory consolidation and visual image processing. It also participates in the regional network supporting oral declarative memory (Bremner et al., 2004; Malhi et al., 2007; Fusar-Poli et al., 2009; Tao et al., 2013; Kukolja et al., 2016). Some previous studies have shown that visual information is involved in regulating sensory and motor function (Anguera et al., 2011; Seidler and Carson, 2017; Fisher et al., 2019). Functional changes in visual-related areas have also been reported in one of our previous studies (Zheng et al., 2019). The anatomical changes in the visual-related region implied that the visual association cortex may be involved in regulating motor function following MSA, although the exact mechanism remains unclear.
Some early studies reported that the cerebellum’s main function was motor control (Andermann et al., 1975; Takagi et al., 1998), but recently many researchers reported that different subregions of cerebellum are involved in the formation of different functions, such as cognition, behavior and learning (Strick et al., 2009; Reeber et al., 2013). Consistent with previous studies, functional changes in cognitive-related cerebellar subregions (Cbe1, Cbe7b and Cbe8) were observed in current study (Desmond et al., 1997; Chen and Desmond, 2005a,b; Kirschen et al., 2005; Hayter et al., 2007). MFG and IPL as the primary hubs of the DMN (Buckner et al., 2008; Ward et al., 2014), are responsible for cognitive functions. The decreased FC between right MFG and left Cbe1, Cbe2, as well as between right IPL and left Cbe7b in the present study may suggest the cognitive dysfunction in MSA-c patients, which is partly consistent with one of our previous studies (Zheng et al., 2019). In a previous study, researchers found reduced activity in the mid-prefrontal cortex and altered functional connections to the insula, precuneus and inferior parietal in MSA patients (Yang et al., 2020). By using independent component analysis (ICA) and dual-regression analysis, researchers found the cognitive deficits may be caused by reduced cerebello-prefrontal and cerebello-amygdala functional connections in MSA (Kawabata et al., 2019). Our findings provide further evidence of cognitive impairment in MSA-c patients.
On the other hand, increased functional connectivity between left Cbe7b and bilateral HES, as well as between right mOFC and right Cbe6, Cbe8, and between right ACC and right Cbe6 were also found in the present study. Previous researches show that the ACC plays an important role in the processing of emotions (Albert et al., 2012) and is linked to mood disorders, including anxiety and depression associated with chronic pain (Barthas et al., 2015). In this study, more than 90% of the MSA patients had anxiety and depression, which may be the reason for the increased FC related to ACC. Additionally, studies have shown that individuals with non-suicidal self-injury behavior show hyperactivation in frontal areas, including the mOFC (Osuch et al., 2014). This indicates that psychological changes of MSA patients might cause functional alterations. In addition, HES belongs to the auditory cortex in humans. Studies have reported that MSA patients have a tendency to increase high-frequency auditory thresholds (Scarpa et al., 2020) and low-frequency auditory thresholds (Omichi et al., 2016; Scarpa et al., 2020), although it has also been reported that MSA patients do not experience significant auditory changes (Ikeda et al., 2013). As participants in this study did not receive audio-vestibular examination, we can only speculate that the increased functional connectivity is due to functional compensation. Further researches on auditory and vision function in MSA are needed.
Relationships Between the Multimodal Magnetic Resonance Features and Clinical Performances
In current study, negative associations were found between clinical measurements (UMSARS-total scores) and the reduced fiber bundle counts of right Cbe3 and decreased FA values of bilateral Cbe9, suggesting clinical relevance of white matter microstructural abnormalities in MSA-c. The white matter microstructural abnormalities in these regions might be used as significant imaging markers for assessing the movement disorders.
The Biomarkers to Differentiate MSA-c Patients From Healthy Controls Using Support Vector Machine Analysis
In current study, we first simultaneous applied the structural, DTI and extended-NBS analysis on MSA-c patients and HCs, and then used all significant features of each of the 116 regions to classify the two groups. Methods such as multilayer graph analysis and fusion ICA have been proposed to merge multimodal MRI data. However, multilayer graph analysis mainly works for merging brain networks constructed from different modalities (Battiston et al., 2017; Mandke et al., 2018). And fusion ICA is useful in combining contrast maps obtained from task fMRI data (Calhoun et al., 2006; Sui et al., 2013). In this work, we collected resting-state fMRI and constructed brain networks to model functional activities. Moreover, the features used in this work are not limited to networks. We consider regional features extracted from T1-weighted and DTI data, which lacks pairwise interaction information between regions. As a result, the AAL atlas was utilized to convert features from different modalities to the atlas space or the feature space. Regional features and networks were represented as feature vectors and matrices. And significant features identified from statistical analysis were used for classification training.
Results showed that the significant features of left Cbe9, Cbe1 and Cbe7b could distinguish MSA-c from HC with relative high accuracy. As previous studies reported, Cbe9 is thought to be crucial for visual guidance of movement (Glickstein et al., 1994). The function of Cbe1 is predominantly sensorimotor, whereas Cbe7b contributes to higher-level processes (Schmahmann, 1991, 1996, 2004; Schmahmann et al., 2009). Current result has an important implication since these regions could be used as valuable imaging biomarkers for the early diagnosis and prognosis assessment of MSA-c patients.
Future Considerations
This work has several limitations. First, the present study showed alterations in visual, auditory and cognitive related regions in MSA patients. However, in this study, we only focused on changes in motor function. A study including more subjects with visual and cognitive function evaluations will be conducted in the future. Next, we have only incorporated one subtype of MSA in this study. In the future, exploring MSA of different subtypes (parkinsonian and cerebellar variants) would be helpful to deepen the understanding of pathological mechanism of MSA. Moreover, a longitudinal design including the pre-symptomatic patients is of great significance to elucidate the progression of MSA. Also, it would be more clinically meaningful to consider a dynamic and longitudinal context in future studies, including the neuroplasticity mechanisms in MSA and other neurodegenerative diseases (NDDs) (Gatto et al., 2019, 2021; Gatto, 2020). We also note that novel multimodal fusion methods that combine regional features and brain networks at the same time could possibly further help to improve classification accuracy.
Conclusion
In conclusion, by simultaneously applying structural (T1-weighted imaging), DTI, functional (BOLD fMRI) and extended-NBS analysis to differentiate MSA-c from HCs, we observed significant grey matter atrophy in cerebellum and white matter microstructural abnormalities in cerebellum and several cerebrum regions, as well as altered functional connectivity between several cerebellum subregions and cerebrum regions in MSA-c patients. In addition, significant negative correlations were found between the UMSARS scores and white matter microstructural abnormalities in several cerebellum subregions of MSA-c patients. Finally, we found the significant features of left Cbe9, Cbe1 and Cbe7b could be used as useful imaging biomarkers to distinguish MSA-c from HCs according to the SVM analysis. The results of this study added new evidence for the structural atrophy, white matter microstructural and functional abnormalities of MSA, which may help us to understand the neural pathophysiological mechanism and provide potential biomarkers and new ideas for the accurate diagnosis and treatment of MSA in the future.
Data Availability Statement
The raw data analyzed in this manuscript are not publicly available due to local regulations and ethical concerns. Requests to access the processed features should be directed to ZW, wangzhiqun@126.com and WD, douwb@tsinghua.edu.cn.
Ethics Statement
The studies involving human participants were reviewed and approved by the Medical Research Ethical Committee of Dongfang Hospital of Beijing University of Chinese Medicine. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
ZW, WD, and ZC were responsible for the study concept and design. SR contributed to the acquisition of MRI data. YG devised the method framework, conducted data analysis, and assisted with data interpretation. WZ assisted with T1 data analysis and was responsible for data interpretation. YL assisted with DTI data analysis. WZ, YG, WD, and ZW drafted the manuscript. All authors provided critical revision of the manuscript for important intellectual content, critically reviewed content and approved final version for publication.
Funding
This work was supported by a grant from Tsinghua University Initiative Scientific Research Program (No. 20131089382), the National Natural Science Foundation of China (No. 61171002, 60372023, and 81571648), the Beijing Municipal Natural Science Foundation of China (No. 7182105), the Capital Special Funds for Health Improvement and Research (No. 2018-2-4202), Dongfang Hospital of Beijing University of Traditional Chinese Medicine “1166” talent training project (No. 030903010321), and Dongfang Hospital of Beijing University of Chinese Medicine Independent Selection Project (No. 2018-JYBZZ-JS124).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.799251/full#supplementary-material
Footnotes
References
Albert, J., Lopez-Martin, S., Tapia, M., Montoya, D., and Carretie, L. (2012). The role of the anterior cingulate cortex in emotional response inhibition. Hum Brain Mapp. 33, 2147–2160. doi: 10.1002/hbm.21347
Andermann, E., Andermann, F., Joubert, M., Melançon, D., Karpati, G., and Carpenter, S. (1975). Three familial midline malformtion syndromes of the central nervous system: agenesis of the corpus callosum and anterior horn-cell disease; agenesis of cerebellar vermis; and atrophy of the cerebellar vermis. Birth Def. Orig. Artic. Ser. 11, 269–293.
Anguera, J. A., Reuter-Lorenz, P. A., Willingham, D. T., and Seidler, R. D. (2011). Failure to engage spatial working memory contributes to age-related declines in visuomotor learning. J. Cogn. Neurosci. 23, 11–25. doi: 10.1162/jocn.2010.21451
Barthas, F., Sellmeijer, J., Hugel, S., Waltisperger, E., Barrot, M., and Yalcin, I. (2015). The anterior cingulate cortex is a critical hub for paininduced depression. Biol. Psychiatry 77, 236–245. doi: 10.1016/j.biopsych.2014.08.004
Battiston, F., Nicosia, V., Chavez, M., and Latora, V. (2017). Multilayer motif analysis of brain networks. Chaos: Interdiscipl. J. Nonlin. Sci. 27:047404. doi: 10.1063/1.4979282
Beliveau, V., Krismer, F., Skalla, E., Schocke, M. M., Gizewski, E. R., Wenning, G. K., et al. (2021). Characterization and diagnostic potential of diffusion tractography in multiple system atrophy. Parkinson. Relat. Disord. 85, 30–36. doi: 10.1016/j.parkreldis.2021.02.027
Bremner, J. D., Vythilingam, M., Vermetten, E., Vaccarino, V., and Charney, D. S. (2004). Deficits in hippocampal and anterior cingulate functioning during verbal declarative memory encoding in midlife major depression. Am. J. Psychiatry 161, 637–645. doi: 10.1176/appi.ajp.161.4.637
Brettschneider, J., Irwin, D. J., Boluda, S., Byrne, M. D., Fang, L., and Lee, E. B. (2017). Progression of alpha-synuclein pathology in multiple system atrophy of the cerebellar type. Neuropathol. Appl. Neurobiol. 43, 315–329. doi: 10.1111/nan.12362
Buckner, R. L., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain’s default network: anatomy, function, and relevance to disease. Ann. NY Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011
Calhoun, V. D., Adali, T., Giuliani, N., Pekar, J., Kiehl, K., and Pearlson, G. (2006). Method for multimodal analysis of independent source differences in schizophrenia: combining gray matter structural and auditory oddball functional data. Human Brain Mapp. 27, 47–62. doi: 10.1002/hbm.20166
Chelban, V., Bocchetta, M., Hassanein, S., Haridy, N. A., Houlden, H., and Rohrer, J. D. (2019). An update on advances in magnetic resonance imaging of multiple system atrophy. J. Neurol. 266, 1036–1045. doi: 10.1007/s00415-018-9121-3
Chen, S. H., and Desmond, J. E. (2005a). Temporal dynamics of cerebro-cerebellar network recruitment during a cognitive task. Neuropsychologia 43, 1227–1237. doi: 10.1016/j.neuropsychologia.2004.12.015
Chen, S. H., and Desmond, J. E. (2005b). Cerebrocerebellar networks during articulatory rehearsal and verbal working memory tasks. Neuroimage 24, 332–338. doi: 10.1016/j.neuroimage.2004.08.032
Dash, S. K., Stezin, A., Takalkar, T., George, L., Kamble, N. L., Netravathi, M., et al. (2019). Abnormalities of white and grey matter in early multiple system atrophy: comparison of parkinsonian and cerebellar variants. Eur. Radiol. 29, 716–724. doi: 10.1007/s00330-018-5594-9
Deguchi, K., Ikeda, K., Kume, K., Takata, T., Kokudo, Y., Kamada, M., et al. (2015). Significance of the hot-cross bun sign on T2*-weighted MRI for the diagnosis of multiple system atrophy. J. Neurol. 262, 1433–1439. doi: 10.1007/s00415-015-7728-1
Desmond, J. E., Gabrieli, J. D., Wagner, A. D., Ginier, B. L., and Glover, G. H. (1997). Lobular patterns of cerebellar activation in verbal working-memory and finger-tapping tasks as revealed by functional MRI. J. Neurosci. 17, 9675–9685. doi: 10.1523/JNEUROSCI.17-24-09675.1997
Fisher, Y. E., Lu, J., Alessandro, L., and Wilson, R. L. (2019). Sensorimotor experience remaps visual input to a heading-direction network. Nature 576, 121–125. doi: 10.1038/s41586-019-1772-4
Fusar-Poli, P., Placentino, A., Carletti, F., Allen, P., Landi, P., Abbamonte, M., et al. (2009). Laterality effect on emotional faces processing: ALE meta-analysis of evidence. Neurosci. Lett. 452, 262–267. doi: 10.1016/j.neulet.2009.01.065
Garyfallidis, E., Brett, M., Amirbekian, B., Rokem, A., van der Walt, S., Descoteaux, M., et al. (2014). DIPY, a library for the analysis of diffusion MRI data. Front. Neuroinform. 8:8. doi: 10.3389/fninf.2014.00008
Gatto, R. G. (2020). Molecular and microstructural biomarkers of neuroplasticity in neurodegenerative disorders through preclinical and diffusion magnetic resonance imaging studies. J. Integr. Neurosci. 19, 571–592. doi: 10.31083/j.jin.2020.03.165
Gatto, R. G., Amin, M., Finkielsztein, A., Weissmann, C., Barrett, T., Lamoutte, C., et al. (2019). Unveiling early cortical and subcortical neuronal degeneration in ALS mice by ultra-high field diffusion MRI. Amyotroph. Later. Scler. Frontotemp. Deg. 20, 549–561. doi: 10.1080/21678421.2019.1620285
Gatto, R. G., Weissmann, C., Amin, M., Angeles-López, Q. D., García-Lara, L., Castellanos, L. C. S., et al. (2021). Evaluation of early microstructural changes in the R6/1 mouse model of Huntington’s disease by ultra-high field diffusion MR imaging. Neurobiol. Aging 102, 32–49. doi: 10.1016/j.neurobiolaging.2021.02.006
Ge, Y., Yang, Z., Feng, Y., Pan, Y., and Dou, W. (2021). “Extended Network-Based Statistics for Measuring Altered Directed Connectivity Components in the Human Brain,” in In 2021 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), (Piscataway: IEEE), 3720–3727.
Gilman, S., Wenning, G. K., Low, P. A., Brooks, D. J., Mathias, C. J., Trojanowski, J. Q., et al. (2008). Second consensus statement on the diagnosis of multiple system atrophy. Neurology 71, 670–676. doi: 10.1212/01.wnl.0000324625.00404.15
Glickstein, M., Gerrits, N., Kralj-Hans, I., Mercier, B., Stein, J., and Voogd, J. (1994). Visual pontocerebellar projections in the macaque. J. Comp. Neurol. 349, 51–72. doi: 10.1002/cne.903490105
Hayter, A. L., Langdon, D. W., and Ramnani, N. (2007). Cerebellar contributions to working memory. Neuroimage 36, 943–954. doi: 10.1016/j.neuroimage.2007.03.011
Ikeda, Y., Ohta, Y., Kurata, T., Shiro, Y., Takao, Y., and Abe, K. (2013). Acoustic impairment is a distinguishable clinical feature of Asidan/SCA36. J. Neurol. Sci. 324, 109–112. doi: 10.1016/j.jns.2012.10.013
Kawabata, K., Hara, K., Watanabe, H., Bagarinao, E., Ogura, A., Masuda, M., et al. (2019). Alterations in cognition-related cerebello-cerebral networks in multiple system atrophy. Cerebellum 18, 770–780. doi: 10.1007/s12311-019-01031-7
Kirschen, M. P., Chen, S. H., Schraedley-Desmond, P., and Desmond, J. E. (2005). Load- and practice-dependent increases in cerebro-cerebellar activation in verbal working memory: an fMRI study. Neuroimage 24, 462–472. doi: 10.1016/j.neuroimage.2004.08.036
Krismer, F., and Wenning, G. K. (2017). Multiple system atrophy: insights into a rare and debilitating movement disorder. Nat. Rev. Neurol. 13, 232–243. doi: 10.1038/nrneurol.2017.26
Krismer, F., Seppi, K., Wenning, G. K., Papapetropoulos, S., Abler, V., Goebel, G., et al. (2019). Abnormalities on structural MRI associate with faster disease progression in multiple system atrophy. Parkins. Relat. Disord. 58, 23–27. doi: 10.1016/j.parkreldis.2018.08.004
Kukolja, J., Goreci, D. Y., Onur, O. A., Riedl, V., and Fink, G. R. (2016). Resting-state fMRI evidence for early episodic memory consolidation: effects of age. Neurobiol. Aging 45, 197–211. doi: 10.1016/j.neurobiolaging.2016.06.004
Lin, J., Xu, X., Hou, Y., Yang, J., and Shang, H. (2020). Voxel-based meta-analysis of gray matter abnormalities in multiple system atrophy. Front. Aging Neurosci. 12:591666. doi: 10.3389/fnagi.2020.591666
Loh, K. B., Rahmat, K., Lim, S. Y., and Ramli, N. (2011). A Hot Cross Bun sign from diffusion tensor imaging and tractography perspective. Neurol. India 59, 266–269. doi: 10.4103/0028-3886.79143
Malhi, G. S., Lagopoulos, J., Owen, A. M., Ivanovski, B., Shnier, R., and Sachdev, P. (2007). Reduced activation to implicit affect induction in euthymic bipolar patients: an fMRI study. J. Affect. Disord. 97, 109–122. doi: 10.1016/j.jad.2006.06.005
Mandke, K., Meier, J., Brookes, M. J., O’Dea, R. D., Van Mieghem, P., Stam, C. J., et al. (2018). Comparing multilayer brain networks between groups: Introducing graph metrics and recommendations. NeuroImage 166, 371–384. doi: 10.1016/j.neuroimage.2017.11.016
Matsusue, E., Fujii, S., Kanasaki, Y., Kaminou, T., Ohama, E., and Ogawa, T. (2009). Cerebellar lesions in multiple system atrophy: postmortem MR imaging-pathologic correlations. AJNR Am. J. Neuroradiol. 30, 1725–1730. doi: 10.3174/ajnr.A1662
Morisi, R., Manners, D. N., Gnecco, G., Lanconelli, N., Testa, C., Evangelisti, S., et al. (2018). Multi-class parkinsonian disorders classification with quantitative MR markers and graph-based features using support vector machines. Parkins. Relat. Disord. 47, 64–70. doi: 10.1016/j.parkreldis.2017.11.343
Omichi, R., Maeda, Y., Sugaya, A., Kataoka, Y., Kariya, S., Nagayasu, R., et al. (2016). Characteristics of audiogram configuration in multiple-system atrophy C and cortical cerebellar atrophy. Acta Otolaryngol. 136, 266–270. doi: 10.3109/00016489.2015.1104722
Osuch, E., Ford, K., Wrath, A., Bartha, R., and Neufeld, R. (2014). Functional MRI of pain application in youth who engaged in repetitive non-suicidal self-injury vs. psychiatric controls. Psychiatry Res. 223, 104–112. doi: 10.1016/j.pscychresns.2014.05.003
Reeber, S. L., Otis, T. S., and Sillitoe, R. V. (2013). New roles for the cerebellum in health and disease. Front. Syst. Neurosci. 7:83. doi: 10.3389/fnsys.2013.00083
Ren, S., Zhang, H., Zheng, W., Liu, M., Gao, F., Wang, Z., et al. (2019). Altered functional connectivity of cerebello-cortical circuit in multiple system atrophy (Cerebellar-Type). Front. Neurosci. 12:996. doi: 10.3389/fnins.2018.00996
Rosskopf, J., Gorges, M., Müller, H. P., Pinkhardt, E. H., Ludolph, A. C., and Kassubek, J. (2018). Hyperconnective and hypoconnective cortical and subcortical functional networks in multiple system atrophy. Parkins. Relat. Disord. 49, 75–80. doi: 10.1016/j.parkreldis.2018.01.012
Rulseh, A., Keller, J., Rusz, J., Syka, M., Brozova, H., Rusina, R., et al. (2016). Diffusion tensor imaging in the characterization of multiple system atrophy. Neuropsychiatr. Dis. Treat. 12, 2181–2187. doi: 10.2147/NDT.S109094
Sako, W., Abe, T., Furukawa, T., Oki, R., Haji, S., Murakami, N., et al. (2019). Differences in the intra-cerebellar connections and graph theoretical measures between Parkinson’s disease and multiple system atrophy. J. Neurol. Sci. 400, 129–134. doi: 10.1016/j.jns.2019.03.022
Scarpa, A., Cassandro, C., Vitale, C., Ralli, M., Policastro, A., Barone, P., et al. (2020). A comparison of auditory and vestibular dysfunction in Parkinson’s disease and Multiple System Atrophy. Parkins. Relat. Disord. 71, 51–57. doi: 10.1016/j.parkreldis.2020.01.018
Schmahmann, J. D. (1991). An emerging concept: the cerebellar contribution to higher function. Arch. Neurol. 48, 1178–1187. doi: 10.1001/archneur.1991.00530230086029
Schmahmann, J. D. (1996). From movement to thought: anatomic substrates of the cerebellar contribution to cognitive processing. Human Brain Mapp. 4, 174–198. doi: 10.1002/(SICI)1097-0193(1996)4:3<174::AID-HBM3>3.0.CO;2-0
Schmahmann, J. D. (2004). Disorders of the cerebellum: ataxia, dysmetria of thought, and the cerebellar cognitive affective syndrome. J. Neuropsychiatry Clin. Neurosci. 16, 367–378. doi: 10.1176/jnp.16.3.367
Schmahmann, J. D., Macmore, J., and Vangel, M. (2009). Cerebellar stroke without motor deficit: clinical evidence for motor and nonmotor domains within the human cerebellum. Neuroscience 162, 852–861. doi: 10.1016/j.neuroscience.2009.06.023
Seidler, R. D., and Carson, R. G. (2017). Sensorimotor learning: neurocognitive mechanisms and individual differences. J. Neuroeng. Rehab. 14:74.
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage 23, S208–S219. doi: 10.1016/j.neuroimage.2004.07.051
Strick, P. L., Dum, R. P., and Fiez, J. A. (2009). Cerebellum and nonmotor function. Annu. Rev. Neurosci. 32, 413–434. doi: 10.1146/annurev.neuro.31.060407.125606
Sui, J., He, H., Pearlson, G. D., Adali, T., Kiehl, K. A., Yu, Q., et al. (2013). Three-way (N-way) fusion of brain imaging data based on mCCA+jICA and its application to discriminating schizophrenia. NeuroImage 66, 119–132. doi: 10.1016/j.neuroimage.2012.10.051
Takagi, M., Zee, D. S., and Tamargo, R. J. (1998). Effects of lesions of the oculomotor vermis on eye movements in primate: saccades. J. Neurophysiol. 80, 1911–1931. doi: 10.1152/jn.1998.80.4.1911
Tao, H., Guo, S., Ge, T., Kendrick, K. M., Xue, Z., Liu, Z., et al. (2013). Depression uncouples brain hate circui. Mol. Psychiatry 18, 101–111. doi: 10.1038/mp.2011.127
Timmann, D., and Daum, I. (2010). How consistent are cognitive impairments in patients with cerebellar disorders? Behav. Neurol. 23, 81–100. doi: 10.3233/BEN-2010-0271
Tournier, J. D., Calamante, F., and Connelly, A. (2007). Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. NeuroImage 35, 1459–1472. doi: 10.1016/j.neuroimage.2007.02.016
Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI Single-Subject Brain. NeuroImage 15, 273–289. doi: 10.1006/nimg.2001.0978
Wang, N., Edmiston, E. K., Luo, X., Yang, H., Chang, M., Wang, F., et al. (2017). Comparing abnormalities of amplitude of low-frequency fluctuations in multiple system atrophy and idiopathic Parkinson’s disease measured with resting-state fMRI. Psychiatry Res. Neuroimag. 269, 73–81. doi: 10.1016/j.pscychresns.2017.09.002
Wang, Y. (2018). Data Processing System Design of Multimodal Magnetic Resonance Imaging of Human Brain. Beijing: Tsinghua University.
Ward, A. M., Schultz, A. P., Huijbers, W., Dijk, K. R., Hedden, T., and Sperling, R. A. (2014). The parahippocampal gyrus links the default-mode cortical network with the medial temporal lobe memory system. Hum. Brain Mapp. 35, 1061–1073. doi: 10.1002/hbm.22234
Wenning, G. K., Stefanova, N., Jellinger, K. A., Poewe, W., and Schlossmacher, M. G. (2008). Multiple system atrophy: a primary oligodendrogliopathy. Ann. Neurol. 64, 239–246. doi: 10.1002/ana.21465
Yan, C., and Zang, Y. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4:13.
Yang, H., Luo, X., Yu, H., Guo, M., Cao, C., Li, Y., et al. (2020). Altered resting-state voxel-level whole-brain functional connectivity in multiple system atrophy patients with cognitive impairment. Clin. Neurophysiol. 131, 54–62. doi: 10.1016/j.clinph.2019.09.026
Yang, H., Wang, N., Luo, X., Lv, H., Liu, H., Li, Y., et al. (2019). Cerebellar atrophy and its contribution to motor and cognitive performance in multiple system atrophy. Neuroimag. Clin. 23:101891. doi: 10.1016/j.nicl.2019.101891
Yao, Q., Zhu, D., Li, F., Xiao, C., Lin, X., Huang, Q., et al. (2017). Altered functional and causal connectivity of cerebello-cortical circuits between multiple system atrophy (parkinsonian type) and Parkinson’s disease. Front. Aging Neurosci. 9:266. doi: 10.3389/fnagi.2017.00266
You, H., Wang, J., Wang, H., Zang, Y. F., Zheng, F. L., Meng, C. L., et al. (2011). Altered regional homogeneity in motor cortices in patients with multiple system atrophy. Neurosci. Lett. 502, 18–23. doi: 10.1016/j.neulet.2011.07.015
Zalesky, A., Fornito, A., and Bullmore, E. T. (2010). Network-based statistic: identifying differences in brain networks. Neuroimage 53, 1197–1207. doi: 10.1016/j.neuroimage.2010.06.041
Zheng, W., Ge, Y., Ren, S., Ran, W., Zhang, X., Tian, W., et al. (2020). Abnormal static and dynamic functional connectivity of resting-state fMRI in multiple system atrophy. Aging-US 12:16341. doi: 10.18632/aging.103676
Zheng, W., Ren, S., Zhang, H., Liu, M., Zhang, Q., Chen, Z., et al. (2019). Spatial patterns of decreased cerebral blood flow and functional connectivity in multiple system atrophy (Cerebellar-Type): a combined arterial spin labeling perfusion and resting state functional magnetic resonance imaging study. Front. Neurosci. 13:777. doi: 10.3389/fnins.2019.00777
Keywords: multiple system atrophy, multimodal MRI, microstructure metrics, extended network-based statistics, support vector machine
Citation: Ge Y, Zheng W, Li Y, Dou W, Ren S, Chen Z and Wang Z (2022) Altered Brain Volume, Microstructure Metrics and Functional Connectivity Features in Multiple System Atrophy. Front. Aging Neurosci. 14:799251. doi: 10.3389/fnagi.2022.799251
Received: 21 October 2021; Accepted: 22 April 2022;
Published: 19 May 2022.
Edited by:
Jong-Min Lee, Hanyang University, South KoreaReviewed by:
Giuseppe Pontillo, University of Naples Federico II, ItalyJiajia Zhu, First Affiliated Hospital of Anhui Medical University, China
Copyright © 2022 Ge, Zheng, Li, Dou, Ren, Chen and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Weibei Dou, douwb@tsinghua.edu.cn; Zhigang Chen, chenzhigang64@126.com; Zhiqun Wang, wangzhiqun@126.com
†These authors have contributed equally to this work