Introduction
Alzheimer’s disease (AD) is a chronic neurodegenerative disease with irreversible progression (Sperling et al., 2011; Hyman et al., 2012). Mild cognitive impairment (MCI) is the prodromal stage of AD (Petersen et al., 1999; Gauthier et al., 2006). The primary clinical manifestation of MCI is memory loss (Petersen et al., 2001). Since the progression of AD is irreversible and the treatment of AD has little effect, the detection of MCI is expected to find out the high risk of AD and further prevent its occurrence (Jacobs et al., 2013; Wang et al., 2013). However, the subtle changes of brain microstructure in MCI make it challenging to distinguish the disease from conventional radiography (Pellegrini et al., 2018). Therefore, establishing reliable biomarkers to diagnose MCI in the early stages remains challenging.
According to the biomarker modeling of AD, white matter demyelination has been proved to appear earlier during AD progression than abnormal changes of gray matter and functional connectivity (Zhuang et al., 2012; Lee et al., 2015). Recent studies have shown that genes and protein molecules ultimately cause microstructure changes in specific white matter fibers (Yu et al., 2021; Zhao et al., 2021). Thus, white matter degeneration is a valid biomarker for MCI. Diffusion tensor imaging (DTI) could detect subtle structural changes in white matter fibers, facilitating large-scale non-invasive screening (Tournier, 2019). The DTI index, mean diffusivity (MD), and fractional anisotropy (FA) describe fiber tracts’ state.
Machine learning offers a systematic approach to developing advanced, automatic, and objective classification frameworks for MCI diagnosis. The classification framework mainly includes feature extraction and classification algorithms (Rathore et al., 2017). Although there are plenty of studies on AD classification, there is insufficient research on MCI (Shatte et al., 2019). The classification accuracies for AD were all around 80%, while the accuracies for MCI were only about 60% (Wee et al., 2011; Dyrba et al., 2013, 2015a; Nir et al., 2015; Prasad et al., 2015; Dou et al., 2020). The main reason is that whole-brain changes in AD are so significant that can be classified with high accuracy by the features of white matter, gray matter, and functional connectivity. In contrast, changes in MCI are not obvious for the whole brain. However, the progression of AD is irreversible. Early diagnosis of MCI fascinates high-risk screening of AD in time.
Previous studies focused on the classification algorithm to improve classification performance. Existing researchers found that k-nearest neighbor (KNN) (Ebadi et al., 2017), random forest (RF) (Maggipinto et al., 2017; Wang et al., 2018), and support vector machine (SVM) (Cui et al., 2012; Demirhan et al., 2015; Dyrba et al., 2015a; Nir et al., 2015; Xie et al., 2015; Ahmed et al., 2017) have achieved better classification performance for MCI when taking the white matter as a feature. KNN is based on Euclid’s theorem and is classified by measuring the distance between different features. RF integrates many decision trees into a forest and combines them for predicting. SVM is based on statistical theory to solve two classification problems, mainly introducing kernel function to solve the problem of linear inseparability. However, previous studies only attached importance to algorithms while neglecting the features. In fact, selecting the appropriate modality in the data and extracting suitable features are usually more important than the underlying algorithm (Zhang et al., 2021).
The white matter feature extraction methods for MCI classification mainly include specific fiber tracking and the whole-brain white matter network (WM network). The specific fiber tracking method takes features using the fiber tracts that converge between localized brain regions and calculates the DTI index of voxels on a single fiber tract. However, this method focused on a specific fiber tract and had low classification accuracy (Nir et al., 2015; Dou et al., 2020). In the whole-brain white matter connectivity network approach, DTI images were segmented into several anatomical regions and features based on the metrics calculated from the fibers within these regions (Wee et al., 2011; Prasad et al., 2015). Recent studies had started to use multimodality for feature extraction. DTI, T1-weighted (T1w), and resting-state functional MRI (rs-fMRI) images from different modes were used to capture information from different perspectives (Hinrichs et al., 2011; Liu et al., 2014). However, most studies had concatenated the features from different modalities. The disadvantage is that all features are treated equally, and it provides no way to account for the diverse nature of features extracted from different modalities (Zhang et al., 2011, 2012; Dyrba et al., 2015b). Therefore, the urgent problem for MCI classification is to combine with different modalities before training and develop a dedicated feature fusion strategy.
To address this issue, the brain network theory was introduced. First, the brain can be seen as a network of spatially dispersed brain areas that share information continuously in functional connectivity (Sporns and Betzel, 2016). The different brain regions are connected and have synchronization within the same activation pattern. Different brain regions are temporally synchronized under the same activation pattern by the white matter pathway (Gu et al., 2015). Second, the primary clinical manifestation of MCI is memory loss (Rugg and Vilberg, 2013). The hippocampus is the crucial area for memory processing. Therefore, investigating the abnormal white matter network related to the hippocampus could promote the detection of MCI in time (Benoit and Schacter, 2015). Furthermore, previous studies have also demonstrated that white matter degeneration leads to abnormalities in the functional connectivity of the corresponding brain regions (Kleinschmidt and Vuilleumier, 2013). In the default mode network, damage to the white matter fibers such as fornix and cingulum could weaken the functional connectivity of the medial temporal lobe and precuneus to the hippocampus (Buckner and DiNicola, 2019). Therefore, we hypothesized that white matter connectivity between brain regions associated with the hippocampus could be an important biomarker for MCI recognition. Features extraction of white matter network through hippocampus-related regions would help to improve the MCI classification performance.
This study aimed to extract the effective features of white matter connectivity networks driven by regions related to the hippocampus and improve classification performance. In our study, the elderly people were divided into MCI group and normal control (NC) group. First, the DTI data were obtained, and the orientation distribution function (ODF) was calculated by constrained spherical deconvolution (CSD). Based on the ODF, the fiber tracts of the whole brain were constructed. Second, the DTI index of whole-brain white matter connectivity between each brain region under the automated anatomical labeling (AAL) was calculated to construct the MD and FA brain structure networks as features. Furthermore, the hippocampus was used as the seed point to define the regions of interest (ROIs). Brain regions with high correlation to the hippocampus in functional connectivity were obtained as ROIs compared with the AAL template. Finally, whole-brain white matter connectivity networks were screened by ROIs related to the hippocampus. Different classifiers (e.g., SVM, KNN, and RF) were used to validate the classification performance of the extracted features. The recursive feature elimination (RFE) ranked these features’ contribution to analyze the pathological mechanisms of MCI. Our study proposed a method of feature extraction to improve the MCI classification performance. The study also revealed the pathological mechanisms of MCI by the ranked contribution of features. It would provide effective early aid to MCI diagnosis.
Materials and Methods
Participants
After excluding the left-handed and other brain injury history, a total of 96 subjects met the criteria for inclusion (48 male and 48 female; age: 80.6 ± 5.4 years, mean ± std). All participants were provided written informed consent based on the Helsinki Declaration. The experimental protocol was approved by the Institutional Review Board of Tianjin University and the Ethics Committee of Chang Gung University.
Neuropsychological Behavior Testing
All subjects were tested on the clinical dementia rating (CDR) and the mini-mental status examination (MMSE) scale. The entry criteria for MCI diagnosis were as follows: (i) CDR = 0.5; (ii) 24 ≤ MMSE < 30 for well-educated subjects (education years ≥ 6) or 19 ≤ MMSE < 24 for less-educated subjects (education years < 6). The MCI group was 81.3 ± 3.6 (mean ± std) years with education duration of 5.9 ± 5.4 years (mean ± std). There was no difference between groups in age and education duration (Table 1).
TABLE 1
Enlarge this image.Table 1. Demographic and neuropsychological information of the MCI and NC groups.
Magnetic Resonance Imaging Data Acquisition
The DTI, T1w, and rs-fMRI images of each subject were acquired without personal information. The DTI with 30 diffusion encoding directions was acquired using echo planar imaging (EPI) sequence with parameters: b-value = 1,000 s/mm2, TR/TE = 11,000/104 ms, field of view = 192 × 192 mm2, matrix size = 128 × 128, slice thickness = 2 mm, voxel size = 2 × 2 × 2 mm3, number of slices = 70, and number of excitations = 1. The rs-fMRI was acquired using gradient EPI sequence with parameters: TR/TE = 2,500/27 ms, flip angle = 77°, band width = 2,400 Hz/pixel, field of view = 260 × 260 mm2, matrix size = 64 × 64, voxel size = 3.4 × 3.4 × 3.4 mm3, number of slices = 43, scan time = 360 s, and time points = 180. The T1w imaging was acquired using a 3D magnetization prepared rapid gradient EPI sequence with parameters: TR/TE = 2,530/3.5 ms; TI = 1,100 ms; field of view = 260 × 260 mm2; matrix size = 256 × 256; slice thickness = 1 mm; voxel size = 1 × 1 × 1 mm3; and number of slices = 192.
White Matter Structure Network Establishment and Feature Extraction
The feature extraction of the white matter structure network for MCI classification consists of the following steps (Figure 1):
FIGURE 1
Enlarge this image.Figure 1. Feature extraction in structural connectivity map driven by hippocampus related ROIs. Three kinds of features were extracted for MCI classification, including: (A) The whole brain WM network was acquired from diffusion MRI with all 90 × 90 AAL regions. (B) The HIP related WM network was selected by 43 × 43 ROIs from the whole brain WM network. (C) The significant HIP related WM network was acquired (*p < 0.05, **p < 0.005, and ***p < 0.0005; ranksum test with bonferroni correction).
1. MRI data preprocessing.
2. Reconstruction for white matter fiber tracts by orientation distribution function (ODF). Based on DTI, the ODF was calculated by CSD.
3. Establishment for whole-brain structure network. DTI indexes (MD and FA) of white matter connectivity between the AAL 90 brain region were calculated to establish the whole-brain structure network as the full feature for MCI classification.
4. Feature selection driven by hippocampus-related ROIs. Brain regions with significant functional connectivity to the hippocampus in rs-fMRI were selected as ROIs. The ROIs are used to screen features for the white matter structural network.
5. Machine learning for MCI recognition. Different classifiers (e.g., SVM, KNN, and RF) were used to test the features for MCI recognition.
6. Performance comparison for searching optimal features and classifiers. Classification performances were validated to prove that our feature extraction method was effective. Pathological mechanisms were analyzed with feature sort.
Magnetic Resonance Imaging Data Preprocessing
The DWI data were denoised and corrected for Gibb’s ringing using MRtrix31 and then motion-corrected, and the eddy current distortion was corrected using the eddy tool in FSL (v5.0.11).2 Next, a brain mask was constructed using the Brain Extraction Tool in FSL, and the diffusion tensor fitting was performed at each voxel within the brain mask to generate DTI index maps. The rs-fMRI data were temporal band-pass filtered (0.01–0.10 Hz) and detrended (both linear and quadratic trends). 3D geometrical displacement was used to correct for head motion. Spatial smoothing was performed with a Gaussian filter kernel (FWHM = 6 mm). The entire process was performed using the Statistical Parametric Mapping (SPM) software package,3 in which the head motion parameters, global signals, white matter signals, and cerebrospinal fluid signals were obtained and combined into a complete covariate. Covariate from functional signals was then removed using the Resting-State fMRI Data Analysis Toolkit (REST).4
Reconstruction for White Matter Fiber Tracts
To reconstruct the whole-brain white matter fiber tracts of each subject, fibers were tracked in DTI data (Supplementary Figure 1). (1) The ODF necessary for the fiber tracking algorithm can be obtained using the MRtrix software (see text footnote 1). The CSD was used with maximum spherical harmonic degree = 6 during this process. (2) The ODF was integrated into DSI Studio5 to achieve a 3D reconstruction of the whole brain’s white matter fiber connectivity. The tracking parameters were as follows: Number of tracts = 100,000, Qa_threshold = 0.04, max angle = 60, length constraint = 20.0–450.0 mm, step size = 1.00 mm, and smoothness = 0.5 mm. Seed direction was set to random, and seed position was subvoxel.
Establishment for Whole-Brain Structure Network
Since the AAL is located in the MNI152 standard space, the transfer matrixes between the standard space and native space (Supplementary Figure 2) were required to obtain the AAL brain region for each subject space. (1) The b0 image from the DTI image served as the native space. The T1w image after skull stripping in FSL served as the structure space. The MNI152 template from FSL served as the standard space. (2) With the high-resolution T1w image as the reference, FMRIB’s Linear Image Registration Tool (FLIRT) was used to obtain the transition matrix from native space to structural space (TNS) and the transition matrix from standard space to structural space (TMS). The transition matrix from native space to standard space was TNM = TNS × TMS–1. (3) The transition matrix from standard space to native space was TMN = TMS × TNS–1.
To establish the whole-brain WM network, a 90 × 90 matrix for the whole brain divided by the AAL template was calculated (Supplementary Figure 1). (1) The AAL 90 brain regions in the MNI152 standard space were transformed to each subject’s space separately by applying transfer matrix TMN. (2) For the 90 × 90 matrix, every element records the mean DTI index (MD or FA) of fiber tracts between every two of the AAL 90 brain regions. Also, the whole-brain WM networks were established, including the MD connectivity 90 × 90 matrix and the FA connectivity 90 × 90 matrix. (3) Whole-brain WM networks of 54 NC and 42 MCI were established separately as entire features for MCI recognition.
Feature Selection Driven by Hippocampus-Related Regions of Interest
The primary clinical symptom of MCI is memory loss, and the hippocampus is closely associated with memory. So, the hippocampus was chosen as a seed point to calculate the brain regions that had significant functional connectivity (Supplementary Figure 3). The concrete steps were as follows: (1) All participants’ rs-fMRI images were transferred from the native space to the standard space. (2) The left hippocampus (30, −16, −14) and right hippocampus (−30, −16, −14) with a radius of 3 mm served as the seed points to create the FC map of each participant in the MNI152 standard space. With the same time series, brain regions that activated correlation with hippocampus were calculated in software REST. The Fisher Z transformation was applied so that the results of each participant followed normal distributions. (3) After individual-level analysis, group analysis for MCI and NC was conducted separately by the second-level analysis in SPM. A one-sample T-test and a familywise error correction were employed to revise the statistics of the group analysis results. (4) Significant area (p < 0.05) with the threshold (T-value > 10, size > 27) was selected to obtain the brain regions associated with the functional connectivity of the hippocampus at the group level.
To confirm the locations of the hippocampus-related ROIs, the following steps were conducted: (1) The functional connectivity maps of the MCI and NC groups were compared with AAL templates. The percentage of overlapping voxels taken in the AAL was calculated (Supplementary Figure 3). (2) To select as many ROIs as possible in the AAL and prevent false positives, a volume percentage threshold of 10% was set. If one of the 90 regions in the AAL had a volume fraction greater than the threshold value, the brain region in the AAL was selected as an ROI (Supplementary Figure 4). The selected ROIs are listed in Supplementary Table 1, including their abbreviations. (3) Totally 43 ROIs related to the hippocampus were used to extract features of the whole-brain WM network (90 × 90 matrices), and the HIP-related WM network (43 × 43 matrices) was established. The rank sum test with Bonferroni correction of HIP-related WM network between MCI and NC was used to extract significant features. The significant HIP-related WM network was acquired (300 vectors) (Figure 1).
Machine Learning for Mild Cognitive Impairment Recognition
To train features of different WM networks for MCI classification, the process of machine learning was as follows (Figure 2): (1) The whole-brain WM network, the hippocampus-related WM network (HIP-related WM network), and the significant HIP-related WM network were separately used as features for machine learning. The recursive feature elimination (RFE) algorithm was used for both dimension reduction and features ranking. (2) The machine learning methods included SVM, KNN, and RF. For SVM, four different kernel functions were tried, namely, linear kernel (linear), polynomial kernel (polynomial), radial basis function kernel (rbf), and sigmoid kernel (sigmoid). For KNN, neighbors of 1, 3, 5, and 7 were used. For RF, 50, 100, 150, and 200 decision trees were used. (3) The hold-out method was used to evaluate classification performance; 80% of the data set was selected as the training set and the remaining 20% as the test set. This process was repeated 100 times randomly to obtain the average classification accuracy and the area under the curve (AUC). All the above experiments were based on the Scikit-Learn library.6
FIGURE 2
Enlarge this image.Figure 2. The flowchart for machine learning. Recursive feature elimination was used to search for the optimal feature subset. The classification algorithm was respectively used KNN, RF, SVM (linear, poly, rbf, sigmoid). The hold-out method was repeated 100 times randomly with 80% of the data for training and 20% for testing.
Performance Comparison for Searching Optimal Feature and Classifier
The following steps were conducted to search for optimal features and classifiers. (1) To demonstrate the effectiveness of the feature extraction method, the classification accuracies (100 times hold-out method) were statistically analyzed in the rank sum test between three features (i.e., whole-brain WM network, HIP-related WM network, and significant HIP-related WM network), and the AUC values were compared in the identical coordinate system (Figure 3). (2) To obtain a better classifier for MCI recognition, the classification accuracies (100 times hold-out method) were statistically analyzed in the rank sum test between different classifiers, and the AUC values were compared in the identical coordinate system (Figure 4). (3) To compare which biomarker is more effective for MCI recognition, the classification accuracies (100 times hold-out method) were statistically analyzed in the rank sum test between MD and FA features, and the AUC values were compared in the identical coordinate system (Figure 5). (4) To rank all features’ contribution to the classification, recursive feature elimination (RFE) was performed. Feature ranking order of the most contributing connectivity was selected based on their classification performance (Figure 6A and Table 4). A schematic illustration of degenerated white matter in MCI was made according to the ranking features (Figure 6B).
FIGURE 3
Enlarge this image.Figure 3. Classification performance comparison between different mean diffusivity (MD) feature sets for different classifier. (A) KNN. (B) RF. (C) SVM linear. (D) SVM poly. (E) SVM rbf. (F) SVM sigmoid. ‘red’: the significant HIP related WM network; ‘green’: HIP related WM network; ‘blue’: the whole brain WM network. Classification performance including mean accuracy (mean ± std; *p < 0.05, **p < 0.005, and ***p < 0.0005; ranksum test with bonferroni correction) and AUC of randomly 100 times hold out method.
FIGURE 4
Enlarge this image.Figure 4. Classification performance comparison between different algorithm through optimal feature set. Classification performance including mean accuracy (mean ± std; *p < 0.05, **p < 0.005, and ***p < 0.0005; ranksum test with Bonferroni correction) and AUC of randomly 100 times hold-out method.
FIGURE 5
Enlarge this image.Figure 5. Optimal classification performance comparison between mean diffusivity (MD) and fractional anisotropy (FA) feature sets (the significant HIP related WM network) for different classifier. (A) KNN. (B) RF. (C) SVM linear. (D) SVM poly. (E) SVM rbf. (F) SVM sigmoid. ‘red’: classification performance of mean diffusivity feature sets; ‘gray’: classification performance of fractional anisotropy feature set. Classification performance including mean accuracy (mean ± std; *p < 0.05, **p < 0.005, and ***p < 0.0005; rank sum test with Bonferroni correction) and AUC of randomly 100 times hold out method.
FIGURE 6
Enlarge this image.Figure 6. Feature ranking of MD feature (the significant HIP-related WM network) as the result of RFE in SVM rbf model. (A) Feature ranking order of the most contributing connectivity (25 features). The bar represents connectivity in different groups. “black”: limbic connectivity (15 features). “gray”: Hip-temporal connectivity (5 features). “white”: THA-frontal connectivity (5 features). (B) Schematic illustration of degenerated white matter in MCI. Nodes represented ROIs from AAL templates. Edges represented connectivity; the value represents MD variation, MCI vs. NC.
Results
Participant Characteristics
There were no significant differences in demographic information between the MCI and NC. There were only significant differences in the cognitive scales (CDR and MMSE, Table 1). This indicates that irrelevant variables’ influence was eliminated between the MCI and NC, and the results were credible.
Regions of Interest Related to Hippocampus in Resting-State Functional Magnetic Resonance Imaging
A total of 43 ROIs related to hippocampus in rs-fMRI were obtained by comparing with AAL (Supplementary Figure 4). Full names and abbreviations of the 43 ROIs were shown in Supplementary Table 1. Among them, 15 ROIs belong to the limbic lobe, 16 ROIs belong to the frontal lobe, 5 ROIs belong to the temporal lobe, 3 ROIs belong to the central region, 3 ROIs belong to the parietal lobe, and 1 ROI belongs to the occipital lobe. Most ROIs are concentrated in the limbic, frontal, and temporal lobes.
Classification Performance Promoted by Screening Features in Hippocampus-Related Regions of Interest
For the WM network of MD in the same classifier, the performances of progressively optimized features (whole-brain WM network, HIP-related WM network, and significant HIP related WM network) were significantly promoted (p < 0.05, rank sum test with Bonferroni correction) (Figure 3 and Table 2). It proved that our feature extraction method in this experiment significantly improved the performance of the MCI classification.
TABLE 2
Enlarge this image.Table 2. Classification performance based on mean diffusivity (MD) feature sets.
Comparison of Mild Cognitive Impairment Classification Performance Under Different Algorithms
The three machine learning classifiers (of which SVM contains four kernel functions) were compared separately (Figure 4). The performance in order of average classification accuracy (100 times hold-out method) was as follows: SVM rbf (ACC = 89.4%, AUC = 0.954), SVM sigmoid (ACC = 88.2%, AUC = 0.950), KNN (ACC = 86.9%, AUC = 0.920), SVM linear (ACC = 86.2%, AUC = 0.937), RF (ACC = 84.8%, AUC = 0.935), and SVM ploy (ACC = 78.5%, AUC = 0.951). The performance of the SVM rbf was significantly better than other classifiers. Therefore, training with the significant HIP-related WM network as features using SVM rbf can better recognize MCI.
Comparison of Classification Performance Under Mean Diffusivity and Fractional Anisotropy Features
For the same machine learning classifier, the performance of the feature MD was significantly better than FA (p < 0.05, rank sum test with Bonferroni correction) in the case of both using optimal features for training (Figure 5 and Table 3). Therefore, the biomarker MD is superior to FA for MCI diagnosing from a machine learning perspective.
TABLE 3
TABLE 4
Ranking the Contribution of White Matter Structural Connectivity to Mild Cognitive Impairment Classification
According to the significant HIP-related WM network’s feature ranking order through RFE (Figure 6A and Table 4), schematic illustrations of degenerated white matter in MCI were made (Figure 6B). The most common connections between MCI and NC are HIP-temporal connectivity, limbic connectivity, and THA-frontal connectivity.
Discussion
This study proposes a feature extraction method for whole-brain white matter connectivity networks driven by ROIs related to the hippocampus. The whole-brain WM network, HIP-related WM network, and significant HIP-related WM network were obtained in the process of feature extraction and optimization. Different classification algorithms, such as KNN, RF, and SVM (linear, poly, rbf, sigmoid), were used to test the classification performance. The pathological mechanisms of MCI were also revealed by RFE. Our study found that feature extraction of whole-brain white matter connectivity by hippocampus-related brain regions can significantly improve MCI classification performance. It can be summarized in three points:
1. In the feature, compared with the whole-brain WM network, the HIP-related WM network can significantly promote the MCI classification performance in machine learning.
2. In the algorithm, the classification performance of the SVM rbf is optimum when taking significant HIP-related WM networks as features.
3. In the pathology, the hippocampus- and thalamus-related white matter connectivity greatly contributed to MCI recognition. So, the method that combines MCI pathology and uses a suitable classifier can significantly improve the classification performance, while the ranking of features contribution can reveal MCI pathology in turn.
Mild Cognitive Impairment Classification Performance Significantly Promoted by Feature Extraction of Hippocampus-Related Regions of Interest
For the same classifier on features, the performance on MCI classification was significantly higher for features screening by HIP-related ROIs than the whole brain. First, in our results, the whole-brain WM network, HIP-related WM network, and significant HIP-related WM network had sequentially significantly higher classification performance (Figure 3). Second, in terms of AD pathogenesis, the degeneration of MCI first appears in specific fiber tracts and gradually spreads to the whole brain when developing to the AD stage (Daianu et al., 2015; Jones et al., 2015; Wang et al., 2016). Furthermore, in either specific fiber tracking (Nir et al., 2015; Dou et al., 2020) or whole-brain white matter connectivity network measures (Wee et al., 2011; Prasad et al., 2015) in previous studies, the accuracies of MCI classification were around 60%, while the accuracies of AD can be about 80%. It also corroborates the pathogenesis of AD from an engineering perspective. Finally, for MCI recognition, a single use of whole-brain white matter would weaken the features considering AD’s pathogenesis. Based on the pathology of MCI, our study improves the MCI classification performance effectively by feature extraction of HIP-related ROIs.
Mean Diffusivity Can Be a Valid Biomarker for Mild Cognitive Impairment Recognition
For MCI, the MD is more sensitive than FA values to reflect white matter degeneration. First, our results showed that the MD of WM networks as features outperformed FA in all the machine learning classifiers (Figure 5). Second, research showed the sensitivity of the MD index to MCI (Yu et al., 2017), while FA only begins to have statistical differences in multiple fiber tracts during AD (Mito et al., 2018). Our experiment confirmed that the MD index was more sensitive in the MCI stage from classification performance. Furthermore, FA reflects the density of white matter fibers and is influenced by the axon diameter of fiber tracts. MD represents the integrity of the myelin outer the fiber tracts. Demyelination can cause MD to increase (Tournier, 2019). It can be inferred that disintegrated myelin appeared more obvious than axon loss during MCI. Therefore, the MD index can be an effective biomarker for MCI compared with FA.
Comparison With Previous Studies
Compared with previous studies that used white matter as the feature to recognize MCI, our method using HIP-related ROIs to extract white matter features significantly improved MCI classification performance (Table 5). By comparing the classification performance between all of our classifiers, the SVM classifier with rbf kernel significantly outperformed the other classifiers (Figure 4). It can also be seen from previous studies that the SVM classifier with rbf kernel for MCI recognition performed well when taking the white matter as features (Table 5). Therefore, training with the white matter as features after HIP-related ROIs extraction can better recognize MCI when using the SVM rbf classifier.
TABLE 5
Important Hubs of the White Matter Connectivity Network for Mild Cognitive Impairment: Hippocampus and Thalamus
White matter connections related to hippocampal and thalamic contribute most to MCI classification. First, previous studies have shown that the white matter associated with the hippocampus and thalamus degenerates during MCI. The white matter microstructure between the hippocampus and medial temporal lobe changes during MCI (Teipel et al., 2016; Zhuo et al., 2016; Dumont et al., 2019). The fornix connects the hippocampus to the thalamus, and its white matter degeneration leads to decreased memory function (Christiansen et al., 2016). Degeneration of the projection tracts connecting the thalamus to the prefrontal lobes leads to a decrease in executive function (Gu and Zhang, 2019; Liu et al., 2021).
Second, in our results, features connected with the hippocampus or thalamus in the white matter connectivity network ranked high in contributions for MCI classification, and features connected with the hippocampus ranked more advanced than the thalamus (Figure 6). This result also confirmed our previous findings (Zhou et al., 2021). In our previous study, the MD of all voxels in a single fiber tract was taken as features. The fiber tract with high separability for MCI recognition passed through the hippocampus or thalamus. So, whether single fiber tract or whole-brain network had corroborated that hippocampus and thalamus were important hubs of the white matter connectivity network for MCI.
Furthermore, in our previous study, the highest average accuracy for MCI classification reached 71.0% when taking MD of all voxels in the left inferior longitudinal fasciculus. In this study, the highest average accuracy reached 89.4% by taking white matter networks related to the hippocampus. Compared with the single fiber tract, white matter networks related to the hippocampus as features can improve classification performance. Finally, the contribution of every white matter connectivity can be sorted by RFE. Thus, the hippocampus and thalamus are important hubs for MCI. Features of combinational white matter connectivity networks outperform single fiber tracts.
Limitations
This study mainly focused on feature extraction for MCI recognition. However, the valuable contributions of this study must be considered within the context of certain limitations. First, a certain amount of new data will be added as the test set alone. Data in this article will be used as the training set and validation set. Second, cross-validation will be added to obtain more accurate parameters for the classifier further to improve performance in the future. Third, deep learning will be used based on our existing findings, and generalization could be improved in the future. More subjects and multiple datasets will be acquired from different hospitals to test the generalization of the classifier. Furthermore, future studies should recruit participants with both MCI and AD. The NC, MCI, and AD should be divided to investigate the pathological mechanisms underlying AD development. Finally, the behavior scale would be added, and the support vector regression (SVR) will be used to predict MCI patients’ behavior.
Conclusion
Our study proposes a feature extraction method driven by hippocampus-related ROIs for white matter connectivity networks. In the feature extraction process, the whole-brain WM network, the HIP-related network, and the significant HIP-related network continuously optimize the performance of MCI classification. By recursive feature elimination, the pathological mechanism revealed that the hippocampus and thalamus are important hubs in white matter networks for MCI. Our results provide a valid basis for the early diagnosis of AD.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. The main codes that support the findings of this study are available on GitHub (https://github.com/ThreePoundUniverse/MCI-network-classification). Further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by Institutional Review Board of Tianjin University; Ethics Committee of Chang Gung University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
XS and YZ designed and conceptualized the research and wrote the manuscript. YZ, Y-PC, and C-PL collected the data. YZ, YC, SL, XZ, YS, and XS analyzed the data. XS, DM, and QL supervised the study. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (grant nos. 61906132 to XS, to DM, 61471263 to QL), Key Project and Team Program of Tianjin City (grant no. XC202020 to XS), Natural Science Foundation of Tianjin City (grant no. 16JCZDJC31100 to QL), Chang Gung University (grant no. BMRPC78 to Y-PC), and Tianjin University (grant no. 2020XRY-0015 to XS).
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.866230/full#supplementary-material
Footnotes
1. ^ http://www.brain.org.au/software/mrtrix/
2. ^ https://fsl.fmrib.ox.ac.uk/fsl
3. ^ http://www.fil.ion.ucl.ac.uk/spm
4. ^ http://restfmri.net/forum/REST
5. ^ http://dsi-studio.labsolver.org/
6. ^ https://scikit-learn.org/
Ahmed, O. B., Benois-Pineau, J., Allard, M., Catheline, G., and Amar, C. B. (2017). Recognition of Alzheimer’s disease and Mild Cognitive Impairment with multimodal image-derived biomarkers and Multiple Kernel Learning. Neurocomputing 220, 98–110. doi: 10.1016/j.neucom.2016.08.041
Benoit, R. G., and Schacter, D. L. (2015). Specifying the core network supporting episodic simulation and episodic memory by activation likelihood estimation. Neuropsychologia 75, 450–457. doi: 10.1016/j.neuropsychologia.2015.06.034
Buckner, R. L., and DiNicola, L. M. (2019). The brain’s default network: updated anatomy, physiology and evolving insights. Nat. Rev. Neurosci. 20, 593–608. doi: 10.1038/s41583-019-0212-7
Christiansen, K., Aggleton, J. P., Parker, G. D., O’Sullivan, M. J., Vann, S. D., and Metzler-Baddeley, C. (2016). The status of the precommissural and postcommissural fornix in normal ageing and mild cognitive impairment: an MRI tractography study. Neuroimage 130, 35–47. doi: 10.1016/j.neuroimage.2015.12.055
Cui, Y., Wen, W., Lipnicki, D. M., Beg, M. F., Jin, J. S., Luo, S., et al. (2012). Automated detection of amnestic mild cognitive impairment in community-dwelling elderly adults: a combined spatial atrophy and white matter alteration approach. Neuroimage 59, 1209–1217. doi: 10.1016/j.neuroimage.2011.08.013
Daianu, M., Jahanshad, N., Nir, T. M., Jack, C. R. Jr., Weiner, M. W., Bernstein, M. A., et al. (2015). Rich club analysis in the Alzheimer’s disease connectome reveals a relatively undisturbed structural core network. Hum. Brain Mapp. 36, 3087–3103. doi: 10.1002/hbm.22830
Demirhan, A., Nir, T. M., Zavaliangos-Petropulu, A., Jack, C. R., Weiner, M. W., Bernstein, M. A., et al. (2015). “Feature selection improves the accuracy of classifying Alzheimer disease using diffusion tensor images,” in Proceedings of the 2015 IEEE 12th International Symposium on Biomedical Imaging (ISBI), New York, NY, 126–130. doi: 10.1109/ISBI.2015.7163832
Dou, X., Yao, H., Feng, F., Wang, P., Zhou, B., Jin, D., et al. (2020). Characterizing white matter connectivity in Alzheimer’s disease and mild cognitive impairment: an automated fiber quantification analysis with two independent datasets. Cortex 129, 390–405. doi: 10.1016/j.cortex.2020.03.032
Dumont, M., Roy, M., Jodoin, P. M., Morency, F. C., Houde, J. C., Xie, Z., et al. (2019). Free Water in White Matter Differentiates MCI and AD From Control Subjects. Front. Aging Neurosci. 11:270. doi: 10.3389/fnagi.2019.00270
Dyrba, M., Barkhof, F., Fellgiebel, A., Filippi, M., Hausner, L., Hauenstein, K., et al. (2015a). Predicting Prodromal Alzheimer’s disease in subjects with mild cognitive impairment using machine learning classification of multimodal multicenter diffusion-tensor and magnetic resonance imaging data. J. Neuroimaging 25, 738–747. doi: 10.1111/jon.12214
Dyrba, M., Grothe, M., Kirste, T., and Teipel, S. J. (2015b). Multimodal analysis of functional and structural disconnection in Alzheimer’s disease using multiple kernel SVM. Hum. Brain Mapp. 36, 2118–2131. doi: 10.1002/hbm.22759
Dyrba, M., Ewers, M., Wegrzyn, M., Kilimann, I., Plant, C., Oswald, A., et al. (2013). Robust automated detection of microstructural white matter degeneration in Alzheimer’s disease using machine learning classification of multicenter DTI Data. PLoS One 8:e64925. doi: 10.1371/journal.pone.0064925
Ebadi, A., Dalboni da Rocha, J. L., Nagaraju, D. B., Tovar-Moll, F., Bramati, I., Coutinho, G., et al. (2017). Ensemble Classification of Alzheimer’s disease and mild cognitive impairment based on complex graph measures from diffusion tensor images. Front. Neurosci. 11:56. doi: 10.3389/fnins.2017.00056
Gauthier, S., Reisberg, B., Zaudig, M., Petersen, R. C., Ritchie, K., Broich, K., et al. (2006). Mild cognitive impairment. Lancet 367, 1262–1270.
Gu, L., and Zhang, Z. (2019). Exploring structural and functional brain changes in mild cognitive impairment: a whole brain ALE meta-analysis for multimodal MRI. ACS Chem. Neurosci. 10, 2823–2829. doi: 10.1021/acschemneuro.9b00045
Gu, S., Pasqualetti, F., Cieslak, M., Telesford, Q. K., Yu, A. B., Kahn, A. E., et al. (2015). Controllability of structural brain networks. Nat. Commun. 6:8414.
Hinrichs, C., Singh, V., Xu, G. F., Johnson, S. C., and Neuroimaging, A. D. (2011). Predictive markers for AD in a multi-modality framework: An analysis of MCI progression in the ADNI population. Neuroimage 55, 574–589. doi: 10.1016/j.neuroimage.2010.10.081
Hyman, B. T., Phelps, C. H., Beach, T. G., Bigio, E. H., Cairns, N. J., Carrillo, M. C., et al. (2012). National Institute on Aging-Alzheimer’s Association guidelines for the neuropathologic assessment of Alzheimer’s disease. Alzheimers Dement. 123, 1–11.
Jacobs, H. I., van Boxtel, M. P., Gronenschild, E. H., Uylings, H. B., Jolles, J., and Verhey, F. R. (2013). Decreased gray matter diffusivity: A potential early Alzheimer’s disease biomarker? Alzheimers Dement. 9, 93–97. doi: 10.1016/j.jalz.2011.11.004
Jones, D. T., Knopman, D. S., Gunter, J. L., Graff-Radford, J., Vemuri, P., Boeve, B. F., et al. (2015). Cascading network failure across the Alzheimer’s disease spectrum. Brain 139, 547–562. doi: 10.1093/brain/awv338
Kleinschmidt, A., and Vuilleumier, P. (2013). Disconnecting cognition. Curr. Opin. Neurol. 26, 333–338. doi: 10.1097/wco.0b013e328363393b
Lee, S. H., Coutu, J. P., Wilkens, P., Yendiki, A., Rosas, H. D., and Salat, D. H. (2015). TRACT-BASED ANALYSIS OF WHITE MATTER DEGENERATION IN ALZHEIMER’S DISEASE. Neuroscience 301, 79–89. doi: 10.1016/j.neuroscience.2015.05.049
Liu, F. Y., Zhou, L. P., Shen, C. H., and Yin, J. P. (2014). Multiple Kernel learning in the primal for multimodal Alzheimer’s disease classification. IEEE J. Biomed. Health 18, 984–990. doi: 10.1109/JBHI.2013.2285378
Liu, L., Jiang, H., Wang, D., and Zhao, X. F. (2021). A study of regional homogeneity of resting-state Functional Magnetic Resonance Imaging in mild cognitive impairment. Behav. Brain Res. 402, 113103. doi: 10.1016/j.bbr.2020.113103
Maggipinto, T., Bellotti, R., Amoroso, N., Diacono, D., Donvito, G., Lella, E., et al. (2017). DTI measurements for Alzheimer’s classification. Phys. Med. Biol. 62, 2361–2375. doi: 10.1088/1361-6560/aa5dbe
Mito, R., Raffelt, D., Dhollander, T., Vaughan, D. N., Tournier, J., Salvado, O., et al. (2018). Fibre-specific white matter reductions in Alzheimer’s disease and mild cognitive impairment. Brain 141, 888–902. doi: 10.1093/brain/awx355
Nir, T. M., Villalon-Reina, J. E., Prasad, G., Jahanshad, N., Joshi, S. H., Toga, A. W., et al. (2015). Diffusion weighted imaging-based maximum density path analysis and classification of Alzheimer’s disease. Neurobiol. Aging 36, S132–S140. doi: 10.1016/j.neurobiolaging.2014.05.037
Pellegrini, E., Ballerini, L., Hernandez, M., Chappell, F. M., Gonzalez-Castro, V., Anblagan, D., et al. (2018). Machine learning of neuroimaging for assisted diagnosis of cognitive impairment and dementia: A systematic review. Alzheimers Dement. 10, 519–535. doi: 10.1016/j.dadm.2018.07.004
Petersen, R. C., Doody, R., Kurz, A., Mohs, R. C., Morris, J. C., Rabins, P. V., et al. (2001). Current concepts in mild cognitive impairment. Arch. Neurol. 58, 1985–1992.
Petersen, R. C., Smith, G. E., Waring, S. C., Ivnik, R. J., Tangalos, E. G., and Kokmen, E. (1999). Mild cognitive impairment - Clinical characterization and outcome. Arch. Neurol. 56, 303–308. doi: 10.1001/archneur.56.3.303
Prasad, G., Joshi, S. H., Nir, T. M., Toga, A. W., and Thompson, P. M. (2015). Brain connectivity and novel network measures for Alzheimer’s disease classification. Neurobiol. Aging 36, S121–S131. doi: 10.1016/j.neurobiolaging.2014.04.037
Rathore, S., Habes, M., Iftikhar, M. A., Shacklett, A., and Davatzikos, C. (2017). A review on neuroimaging-based classification studies and associated feature extraction methods for Alzheimer’s disease and its prodromal stages. Neuroimage 155, 530–548. doi: 10.1016/j.neuroimage.2017.03.057
Rugg, M. D., and Vilberg, K. L. (2013). Brain networks underlying episodic memory retrieval. Curr. Opin. Neurobiol. 23, 255–260. doi: 10.1016/j.conb.2012.11.005
Shatte, A. B. R., Hutchinson, D. M., and Teague, S. J. (2019). Machine learning in mental health: a scoping review of methods and applications. Psychol. Med. 49, 1426–1448. doi: 10.1017/S0033291719000151
Sperling, R. A., Aisen, P. S., Beckett, L. A., Bennett, D. A., Craft, S., Fagan, A. M., et al. (2011). Toward defining the preclinical stages of Alzheimer’s disease: recommendations from the National Institute on Aging-Alzheimer’s Association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. 7, 280–292. doi: 10.1016/j.jalz.2011.03.003
Sporns, O., and Betzel, R. F. (2016). Modular brain networks. Annu. Rev. Psychol. 67, 613–640. doi: 10.1146/annurev-psych-122414-033634
Teipel, S., Grothe, M. J., and Alzheimer ‘s Disease Neuroimaging Initiative (2016). Does posterior cingulate hypometabolism result from disconnection or local pathology across preclinical and clinical stages of Alzheimer’s disease? Eur. J. Nucl. Med. Mol. Imaging 43, 526–536. doi: 10.1007/s00259-015-3222-3
Tournier, J. D. (2019). Diffusion MRI in the brain - Theory and concepts. Prog. Nucl. Magn. Reson. Spectrosc. 112-113, 1–16. doi: 10.1016/j.pnmrs.2019.03.001
Wang, J., Zuo, X., Dai, Z., Xia, M., Zhao, Z., Zhao, X., et al. (2013). Disrupted functional brain connectome in individuals at risk for Alzheimer’s disease. Biol. Psychiatry 73, 472–481. doi: 10.1016/j.biopsych.2012.03.026
Wang, Q., Guo, L., Thompson, P. M., Jack, C. R. Jr., Dodge, H., Zhan, L., et al. (2018). The added value of diffusion-weighted mri-derived structural connectome in evaluating mild cognitive impairment: a multi-cohort validation. J. Alzheimers Dis. 64, 149–169. doi: 10.3233/JAD-171048
Wang, T., Shi, F., Jin, Y., Yap, P.-T., Wee, C.-Y., Zhang, J., et al. (2016). Multilevel deficiency of white matter connectivity networks in Alzheimer’s disease: a diffusion MRI study with DTI and HARDI models. Neural Plast. 2016:2947136. doi: 10.1155/2016/2947136
Wee, C. Y., Yap, P. T., Li, W. B., Denny, K., Browndyke, J. N., Potter, G. G., et al. (2011). Enriched white matter connectivity networks for accurate identification of MCI patients. Neuroimage 54, 1812–1822. doi: 10.1016/j.neuroimage.2010.10.026
Xie, Y., Cui, Z., Zhang, Z., Sun, Y., Sheng, C., Li, K., et al. (2015). Identification of amnestic mild cognitive impairment using multi-modal brain features: a combined structural MRI and Diffusion Tensor Imaging Study. J. Alzheimers Dis. 47, 509–522. doi: 10.3233/JAD-150184
Yu, J., Lam, C. L. M., and Lee, T. M. C. (2017). White matter microstructural abnormalities in amnestic mild cognitive impairment: a meta-analysis of whole-brain and ROI-based studies. Neurosci. Biobehav. Rev. 83, 405–416. doi: 10.1016/j.neubiorev.2017.10.026
Yu, M., Sporns, O., and Saykin, A. J. (2021). The human connectome in Alzheimer disease - relationship to biomarkers and genetics. Nat. Rev. Neurol. 17, 545–563. doi: 10.1038/s41582-021-00529-1
Zhang, D. Q., Shen, D. G., and Alzheimer’s Disease Neuroimaging Initiative (2012). Multi-modal multi-task learning for joint prediction of multiple regression and classification variables in Alzheimer’s disease (vol 59, pg 895, 2012). Neuroimage 62, 2179–2179. doi: 10.1016/j.neuroimage.2011.09.069
Zhang, D. Q., Wang, Y. P., Zhou, L. P., Yuan, H., and Shen, D. G. (2011). Multimodal classification of Alzheimer’s disease and mild cognitive impairment. Neuroimage 55, 856–867. doi: 10.1016/j.neuroimage.2011.01.008
Zhang, Y. P., Wang, S. H., Xia, K. J., Jiang, Y. Z., and Qian, P. J. (2021). Alzheimer’s disease multiclass diagnosis via multimodal neuroimaging embedding feature selection and fusion. Inf. Fusion 66, 170–183. doi: 10.1016/j.inffus.2020.09.002
Zhao, B., Li, T., Yang, Y., Wang, X., Luo, T., Shan, Y., et al. (2021). Common genetic variation influencing human white matter microstructure. Science 372:eabf3736. doi: 10.1126/science.abf3736
Zhou, Y., Si, X., Chen, Y., Chao, Y., Lin, C. P., Li, S., et al. (2021). Hippocampus- and thalamus-related fiber-specific white matter reductions in mild cognitive impairment. Cereb. Cortex. [Epub ahead of print]. doi: 10.1093/cercor/bhab407
Zhuang, L., Sachdev, P. S., Trollor, J. N., Kochan, N. A., Reppermund, S., Brodaty, H., et al. (2012). Microstructural white matter changes in cognitively normal individuals at risk of amnestic MCI. Neurology 79, 748–754. doi: 10.1212/WNL.0b013e3182661f4d
Zhuo, J. J., Fan, L. Z., Liu, Y., Zhang, Y. C., Yu, C. S., and Jiang, T. Z. (2016). Connectivity profiles reveal a transition subarea in the parahippocampal region that integrates the anterior temporal-posterior medial systems. J. Neurosci. 36, 2782–2795. doi: 10.1523/JNEUROSCI.1975-15.2016
Yu Zhou1†, Xiaopeng Si2,3,4*†, Yi-Ping Chao5,6, Yuanyuan Chen2,3, Ching-Po Lin7, Sicheng Li2,3, Xingjian Zhang2,3, Yulin Sun2,3, Dong Ming2,3* and Qiang Li1*
* 1School of Microelectronics, Tianjin University, Tianjin, China
* 2Academy of Medical Engineering and Translational Medicine, Tianjin University, Tianjin, China
* 3Tianjin Key Laboratory of Brain Science and Neural Engineering, Tianjin University, Tianjin, China
* 4Institute of Applied Psychology, Tianjin University, Tianjin, China
* 5Graduate Institute of Biomedical Engineering, Chang Gung University, Taoyuan, Taiwan
* 6Department of Computer Science and Information Engineering, Chang Gung University, Taoyuan, Taiwan
* 7Institute of Neuroscience, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
© 2022. This work is licensed under http://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.